Kidney Res Clin Pract > Epub ahead of print
Kim, Oh, Ji, Choi, Oh, Suh, Choi, Bae, Ma, Paik, and Kim: Systematic metabolomics study in the serum and urine of a mouse model of Fabry disease

Abstract

Background

Fabry disease (FD) is an X-linked lysosomal disorder caused by α-galactosidase A enzyme activity deficiency. Although glycosphingolipid analogs have been identified in the plasma or urine of patients with FD, there is a limited understanding of altered metabolomics profiles beyond the globotriaosylceramide accumulation in FD.

Methods

Metabolomics study was performed for monitoring of biomarker and altered metabolism related with disease progression in serum and urine from male α-galactosidase A knockout mice and age-matched wild-type mice at 20 and 40 weeks. Profiling analysis for metabolites, including organic acids, amino acids, fatty acids, kynurenine pathway metabolites, and nucleosides in the serum and urine was performed using gas chromatography-tandem mass spectrometry and liquid chromatography-tandem mass spectrometry combined with star symbol patterns and partial least squares discriminant analysis (PLS-DA).

Results

A total of 27 and 23 metabolites from the serum and urine of FD mice were distinguished from those of wild-type mice, respectively, based on p-value (<0.05) and variable importance in projection scores (>1.0) of PLS-DA. In the serum, metabolites of the glutathione, glutathione disulfide, citrulline, and kynurenine pathways that are related to oxidative stress, nitric oxide biosynthesis, and inflammation were increased, whereas those involved in pyruvate and tyrosine metabolism and the tricarboxylic acid cycle were altered in the 20- and 40-week-old urine of FD model mice.

Conclusion

Altered metabolic signatures associated with disease progression by oxidative stress, inflammation, nitric oxide biosynthesis, and immune regulation in the early and late stages of FD.

Introduction

Fabry disease (FD, OMIM number 301500) is an X-linked inherited lysosomal storage disease caused by mutations of the GLA gene. This biological phenomenon leads to the accumulation of globotriaosylceramide (Gb3) in podocytes, cardiac myocytes, endothelial cells, and neural cells because of defective or absent activity of the α-galactosidase A (α-Gal A) enzyme [1]. This may lead to life-threatening complications of hypertrophic cardiomyopathy, cerebrovascular events, and kidney injury [2]. Moreover, another sub-metabolite of globotriaosylsphingosine (lyso-Gb3), which is a deacylated form of Gb3, has been used as a biomarker for diagnosis and correlation with phenotype. Recently, many Gb3 isoforms and lyso-Gb3 analogs have been found in the plasma and urine of patients with FD [36].
Metabolomics is a global and multifaceted analysis of screening small molecules altered in biological samples [7]. Recent studies have approached the alteration of metabolic pathways in FD using a systematic and comprehensive metabolic analysis [8,9]. Previous metabolomics studies to identify new Gb3 or lyso-Gb3-related biomarkers in the plasma and urine of patients with FD have been conducted using a mass spectrometry metabolic approach, but these were focused on screening accumulation of glycosphingolipid analogs [46,10]. Although Ducatez et al. [11] identified 13 specific metabolites in the plasma of patients with FD from the French Fabry cohort through a network-based analysis, most of the targeted metabolites were glycerophospholipids, amino acids (AAs), and amines.
FD is associated with inflammation and oxidative stress. A recent study showed that α-Gal A–deficient podocyte dysregulated proteins involved in thermogenesis, lysosomal trafficking, metabolic activity, and cell-cell interaction and cycle [12]. In addition, the accumulation of Gb3 is related to an increase in reactive oxygen species (ROS) in endothelial cells, which may affect metabolic imbalance and mitochondrial dysfunction linked to the tricarboxylic acid cycle (TCA) cycle [1315]. Thus, in this study, we focused on the analysis on systematic metabolomics including organic acids (OAs) along with TCA cycle metabolites, AAs, nucleosides (NSs), kynurenine pathway metabolites (KYNs), and fatty acids (FAs) as final products of lipid metabolism in the serum and urine of 20-week-old (20w) and 40-week-old (40w) FD mice compared to age-matched wild-type (WT) mice to reflect disease progression [16]. Our findings are expected to give a clue to the potential biomarker beyond the glycosphingolipids for predicting disease progression and provide therapeutic strategies in patients with FD.

Methods

Animals

Mice were maintained in a 12-hour light/dark cycle and had free access to standard chow and tap water. B6; 129-Glatm1Kul/J mice (α-Gal A knockout mice) were purchased from the Jackson Laboratory (JAX stock #003535). The male α-Gal A knockout (hemizygote) mice used in this study were selected by PCR-based genotyping using ear punch DNA samples. The primers used to amplify the GLA gene were as follows: oIMR5947 5’-AGG TCC ACA GCA AAG GAT TG-3’; oIMR5948 5’-GCA AGT TGC CCT CTG ACT TC-3’; and oIMR7415 5’-GCC AGA GGC CAC TTG TGT AG-3’. We used 20w and 40w α-Gal A knockout mice (n = 6 and 7 mice, respectively). Age-matched WT mice without a GLA mutation were used as healthy controls (n = 7 for each age group). Serum was collected from the retro-orbital sinus and centrifuged at 845 ×g for 10 minutes. Urine samples were collected in metabolic cages to examine the metabolites 2 days before sacrifice. The function data of the kidney and liver, and immunohistochemical staining in 20w and 40w FD mice and WT mice were shown in Supplementary Fig. 1A–E (available online).
All animal experiments were approved by the Institutional Animal Care and Use Committee of Chonnam National University Medical School and were conducted in accordance with the institution’s guidelines for experimental animal care and use (No. CNUHI-ACUC-23014).

Chemicals and reagents

All standards used and internal standards (ISs) for simultaneous profiling analyses of various metabolites, trimethylamine (TEA), and ammonium formate were purchased from Sigma-Aldrich and Tokyo Chemical Industry. Distilled water (DW), methanol, and acetonitrile (ACN) were purchased from J.T. Baker Inc. Formic acid and acetic acid were purchased from Wako Pure Chemical. Diethyl ether (DEE), ethyl acetate (EA), toluene, dichloromethane, and sodium chloride were purchased from Kanto Chemical. Tetrahydrofuran was purchased from Fisher Scientific Korea. Sulfuric acid and sodium hydroxide (NaOH) were purchased from Daejung Reagents Chemicals. O-Methoxyamine hydrochloride and N-methyl-N-(tert-butyldimethylsilyl) trifluoroacetamide (MTBSTFA) + 1% tert-butyldimethylchlorosilane were obtained from Thermo Scientific. All other chemicals were analytical grade.

Gas chromatography-tandem mass spectrometry (GC-MS/MS)

Metabolic profiling analyses of OAs and FAs were conducted by GCMS-TQ8040 (Shimadzu Corp.) equipped with an Ultra-2 (5% phenyl-95% methylpolysiloxane bonded phase; 25 × 0.20-mm inner diameter, 0.11 μm film thickness) cross-linked capillary column (Agilent Technologies) and interfaced with a triple quadrupole mass spectrometer. The profiling analyses were performed in multiple reaction monitoring (MRM) mode in electron impact ionization mode.

Liquid chromatography-tandem mass spectrometry (LC-MS/MS)

Metabolic profiling analyses of AAs, KYNs, and NSs were performed by a Shimadzu Nexera UPLC system (Shimadzu Corp.) coupled with an LCMS-8050 triple quadrupole mass spectrometer (Shimadzu Corp.). For chromatographic separation, an Intrada AA column (150 × 2 mm, 3 μm; Intakt) for AA profiling, a Zorbax Eclipse XDB-C18 (100 × 4.6 mm, 3.5 μm; Agilent Technologies) for KYN profiling, and a Hydro-RP C18 (150 × 4.6 mm, 4 μm; Phenomenex) for NS profiling were used. The MS system was operated in electrospray ionization mode and the profiling analyses were performed in MRM mode.

Sample preparation for organic acid and fatty acid profiling analyses in serum and urine by GC-MS/MS

Serum and urine samples from FD model mice were subjected to methoxime (MO)/tert-butyldimethylsilyl (TBDMS) derivatives for OA and FA profiling analysis by GC-MS/MS, as described in our previous studies [1719]. Briefly, for deproteinization, ACN was added to serum (40 μL) and urine (20 μL) including ISs (13C2-succinic acid, 3,4-methoxybenzoic acid, lauric-d2-acid, and pentadecanoic acid), and they were mixed and centrifuged at 12,300 ×g for 3 minutes. Then, it was transferred to DW, the aqueous phase was adjusted to pH ≥12 with 5-M NaOH, and the MO reaction was performed with O-methoxyamine hydrochloride (1 mg) at 60 ℃ for 60 minutes. The aqueous phase was then acidified to pH ≤2 with 10% sulfuric acid and saturated with sodium chloride. Liquid-liquid extraction was conducted twice with 3 mL of DEE and 2 mL of EA. After 5 μL of TEA was added, the extracts were evaporated to dryness under a smooth stream of nitrogen at 40 ℃. To form the TBDMS derivative, the extracts were reacted with toluene (10 μL) and MTBSTFA (20 μL) for 60 minutes at 60 ℃. Then, 1 μL was injected for GC-MS/MS analysis.

Sample preparation for amino acid, kynurenine pathway metabolite, and nucleoside profiling analyses in serum and urine by LC-MS/MS

For deproteinization, ACN was added to serum (30 μL) and urine (10 μL) including ISs (13C1-phenylalanine and 13C1-leucine for AA profiling, 3,4-dimethoxybenzoic acid for KYN profiling, and 3-deazauridine for NS profiling), and they were mixed and centrifuged at 12,300 ×g for 3 minutes. After the supernatant was transferred to a Spin-X centrifuge filter tube and centrifuged at 12,300 ×g for 3 minutes, 1 μL was injected for LC-MS/MS analysis according to our previous method [20].

Star symbol pattern recognition analysis

The levels of OAs, FAs, AAs, KYNs, and NSs were calculated based on their standard calibration curves. The mean concentrations of each metabolite of the 20w and 40w FD groups were normalized to the corresponding mean values of the WT groups (20w and 40w, respectively). Each normalized value was plotted as a line emanating from a common central point. Star symbol patterns, which are used to readily observe metabolic alterations through a visualized format, were drawn with the normalized values using Microsoft Excel [21,22].

Statistical analysis

Data from animal experiments are presented as the mean ± standard deviation. After converting concentration data to log10-transformed data, a statistical comparison for the univariate analysis was performed using the Wilcoxon rank-sum test to distinguish discriminatory features between the two groups. The significant comparison results are indicated by a p-value of <0.05. Multivariate statistical analyses of metabolomics data, including partial least squares discriminant analysis (PLS-DA), hierarchical clustering heatmap analysis, univariate receiver operating characteristic (ROC) curve analysis, and pathway analysis, were performed using MetaboAnalyst 5.0 (https://www.metaboanalyst.ca) [2325]. The cross validation of PLS-DA was confirmed with parameters, a correlation coefficient (R2), a cross-validation correlation coefficient (Q2), and an accuracy.

Results

Altered metabolic profiles and univariate analysis

Serum

Among 207 metabolites in our database according to the standard solutions, 91 metabolites including 17 OAs, 20 FAs, 36 AAs, seven KYNs, and 11 NSs were determined in the 20w and 40w serum of FD model (Supplementary Table 1, available online). The 20w and 40w FD groups were normalized to the corresponding 20w and 40w WT groups, respectively. Star symbol patterns were drawn to visually monitor altered metabolism and discrimination of WT and FD states from the metabolic profiles (Fig. 1). As a result of the Wilcoxon rank-sum test, eight metabolites, including glutathione disulfide (GSSG), isoleucine, and N4-acetylcytidine, showed significant changes (p < 0.05) between the 20w WT group and the 20w FD group, and 24 metabolites, including 5’-deoxy-5’-methylthioadenosine (MTA), 5,6-dihydrouridine, and N4-acetylcytidine, showed significant changes between the 40w WT group and the 40w FD group.

Urine

A total of 99 metabolites including 52 OAs, 24 AAs, seven KYNs, and 16 NSs were determined in the 20w and 40w urine of FD model mice (Supplementary Table 2, available online). The 20w and 40w FD groups were normalized to the corresponding 20w and 40w WT groups, respectively. Star symbol patterns were drawn to visually monitor altered metabolism and discrimination of WT and FD states from the metabolic profiles (Fig. 2). As a result of the Wilcoxon rank-sum test, nine metabolites, including pyruvic acid, 2-hydroxybutyric acid, and glycolic acid, showed significant changes (p < 0.05) between the 20w WT group and the 20w FD group, and 19 metabolites, including 3-methyladipic acid, ethylmalonic acid, and malonic acid, showed significant changes between the 40w WT group and the 40w FD group.

Multivariate analyses

Serum

To select metabolites that could discriminate between the WT group and the FD group, 27 metabolites (glutamic acid, 1-methylguanosine, MTA, N4-acetylcytidine, arachidonic acid, docosahexaenoic acid, docosatetraenoic acid, tryptophan, pipecolic acid, α-aminoadipic acid, creatine, citrulline, lysine, picolinic acid, 5-hydroxyindoleacetic acid, serotonin, 5,6-dihydrouridine, pseudouridine, N2-methylguanosine, 5-methylcytidine, N2,N2-dimethylguanosine, isoleucine, aspartic acid, GSSG, and cytidine) were considered key metabolites from serum of the 20w and 40w model mice; these were chosen based on p-values (<0.05) and variable importance in projection (VIP) scores (>1.0) from PLS-DA analysis (Table 1). Statistical analyses were performed with the 27 metabolites selected. The WT and FD groups of 20w and 40w model mice were separated in PLS-DA (Fig. 3A, D), but were not classified in heatmap analysis (Supplementary Fig. 2A, B; available online). In the serum of 20w model mice, the R2, Q2, and accuracy of PLS-DA were 0.88, 0.63, and 0.93, respectively, while in the serum of 40w model mice, they were 0.70, 0.58, and 0.87, respectively.

Urine

To select metabolites that could discriminate between the WT and FD groups, the 23 metabolites (2-hydroxybutyric acid, 3-hydroxypropionic acid, 1-methyladenosine, N6-methyladenosine, kynurenic acid, phenylacetic acid, malonic acid, ethylmalonic acid, isovaleryglycine, methylsuccinic acid, hexanoylglycine, adipic acid, 3-methyladipic acid, vanillic acid, 3-indolecarboxylic acid, vanillylmandelic acid, picolinic acid, 5-hydroxyindoleacetic acid, xanthosine, pyruvic acid, glycolic acid, α-ketoglutaric acid, and 4-hydroxyphenylpyruvic acid) were considered as key metabolites from urine of 20w and 40w model mice; these were chosen based on p-values (<0.05) and VIP scores (>1.0) from the PLS-DA analysis (Table 2). Statistical analyses were performed with the 23 metabolites selected. The WT group and the FD group of 20w and 40w model were separated in PLS-DA (Fig. 4A, D), and were classified in heatmap analysis, except for one case (Supplementary Fig. 2C, D; available online). In the urine of 20w model mice, the R2, Q2, and accuracy of PLS-DA were 0.99, 0.54, and 0.70, respectively, while in the urine of 40w model mice, they were 0.99, 0.81, and 0.90, respectively.

Identification of differential metabolites

Serum

Univariate ROC analysis, a practical tool for evaluating the performance of diagnostic tests, was performed on the 27 metabolites. 11 Metabolites in the 20w serum and 23 metabolites in the 40w serum showed an area under the curve (AUC) values (>0.8) (Table 1). Particularly, four metabolites including isoleucine (AUC = 0.95, log2 fold change [FC] = 0.50), aspartic acid (AUC = 0.93, log2FC = –0.60), GSSG (AUC = 1.00, log2FC = –1.49), and cytidine (AUC = 0.86, log2FC = 0.45) in the 20w serum were significantly altered (p < 0.05, VIP score > 1.0, AUC > 0.8) and thus may represent altered metabolites of an early diagnosis of FD. In addition, 19 metabolites including cis-aconitic acid (AUC = 0.88, log2FC = –0.35), arachidonic acid (AUC = 0.90, log2FC = 0.76), docosahexaenoic acid (AUC = 0.90, log2FC = 0.42), docosatetraenoic acid (AUC = 0.84, log2FC = 0.39), tryptophan (AUC = 0.84, log2FC = –0.57), pipecolic acid (AUC = 0.90, log2FC = –0.93), α-aminoadipic acid (AUC = 0.88, log2FC = –1.33), threonine (AUC = 0.84, log2FC = –0.39), creatine (AUC = 0.90, log2FC = –0.47), citrulline (AUC = 0.94, log2FC = –0.67), lysine (AUC = 0.86, log2FC = –0.54), picolinic acid (AUC = 0.86, log2FC = –0.31), 5-hydroxyindoleacetic acid (AUC = 0.86, log2FC = –0.71), serotonin (AUC = 0.92, log2FC = 0.44), 5,6-dihydrouridine (AUC = 0.98, log2FC = 1.25), pseudouridine (AUC = 0.90, log2FC = 0.64), N2-methylguanosine (AUC = 0.90, log2FC = 0.96), 5-methylcytidine (AUC = 0.90, log2FC = 0.56), and N2,N2-dimethylguanosine (AUC = 0.90, log2FC = 0.64) in the 40w serum were significantly altered (p < 0.05, VIP score > 1.0, AUC > 0.8) and thus may represent altered metabolites of progressive FD. Importantly, four metabolites (glutamic acid, 1-methylguanosine, MTA, and N4-acetylcytidine) in both 20w and 40w serum showed changes in equal directions (p < 0.05, VIP score > 1.0, AUC > 0.8). Glutamic acid increased more from 20w to 40w, and 1-methylguanosine, MTA, and N4-acetylcytidine decreased more from 20w to 40w. These metabolites could be interpreted as serum altered metabolites of progressive FD. Furthermore, GSSG was significantly increased in the 20w serum, but this change was not present in the 40w serum.

Urine

Univariate ROC analysis was performed on the 23 metabolites. Nine metabolites in the 20w urine and 19 metabolites in the 40w urine showed AUC values (>0.8) (Table 2). In particular, four metabolites including pyruvic acid (AUC = 1.00, log2FC =–0.89), glycolic acid (AUC = 0.97, log2FC = –0.36), α-ketoglutaric acid (AUC = 0.86, log2FC = –0.95), and 4-hydroxyphenylpyruvic acid (AUC = 0.86, log2FC = –1.02) in the 20w urine were significantly altered (p < 0.05, VIP score > 1.0, AUC > 0.8) and thus may represent altered metabolites of an early diagnosis of FD. In addition, 14 metabolites including phenylacetic acid (AUC = 0.86, log2FC = –0.62), malonic acid (AUC = 0.94, log2FC = –0.77), ethylmalonic acid (AUC = 0.96, log2FC = –1.03), isovaleryglycine (AUC = 0.88, log2FC = –0.58), methylsuccinic acid (AUC = 0.94, log2FC = –0.38), hexanoylglycine (AUC = 0.84, log2FC = –1.17), adipic acid (AUC = 0.88, log2FC = –0.92), 3-methyladipic acid (AUC = 1.00, log2FC = –1.07), vanillic acid (AUC = 0.86, log2FC = 0.87), 3-indolecarboxylic acid (AUC = 0.94, log2FC = –0.31), vanillylmandelic acid (AUC = 0.88, log2FC = –0.43), picolinic acid (AUC = 0.86, log2FC = –2.34), 5-hydroxyindoleacetic acid (AUC = 0.90, log2FC = –0.34), and xanthosine (AUC = 0.94, log2FC = –0.81) in the 40w urine were significantly altered (p < 0.05, VIP score > 1.0, AUC > 0.8) and thus may represent altered metabolites of progressive FD. Importantly, four metabolites (2-hydroxybutyric acid, 3-hydroxypropionic acid, 1-methyladenosine, and N6-methyladenosine) in both 20w and 40w urine showed changes in equal directions. More specifically, 2-hydroxybutyric acid and 3-hydroxypropionic acid increased in both the 20w and 40w urine, while 1-methyladenosine and N6-methyladenosine decreased in both the 20w and 40w urine. Thus, these metabolites could be interpreted as altered metabolites of progressive FD in urine. Interestingly, kynurenic acid was significantly decreased in the 20w urine, while it was significantly increased in the 40w urine, showing a change in the opposite direction.

Metabolic pathway analysis

Serum

To confirm which metabolic pathways the 27 metabolites were associated with, pathway analysis was performed. Alanine, aspartate, and glutamate metabolism, arginine and proline metabolism, arginine biosynthesis, and D-glutamine and D-glutamate metabolism were evaluated as key pathways in the 20w serum (p < 0.05, impact score > 0.1) (Fig. 3C). In addition, arginine and proline metabolism, arginine biosynthesis, D-glutamine and D-glutamate metabolism, tryptophan metabolism, arachidonic acid metabolism, and lysine degradation were evaluated as key pathways in the 40w serum (p < 0.05, impact score > 0.1) (Fig. 3F).

Urine

To confirm which metabolic pathways the 23 metabolites were associated with, pathway analysis was performed. Pyruvate metabolism, glycolysis/gluconeogenesis, and the TCA cycle were evaluated as a key pathway in the 20w urine (p < 0.05, impact score > 0.1) (Fig. 4C). Tyrosine metabolism (p = 0.221, impact score = 0.08) and the TCA cycle (p = 0.61, impact score = 0.10) were found to be related pathways in the 40w urine (Fig. 4F).

Discussion

In the present study, we conducted metabolic profiling analysis to identify 91 metabolites in the serum and 99 metabolites in the urine of mice with FD. Multivariate analysis was used to identify differential metabolites. We selected 27 key metabolites associated with FD that were different from those found in the serum of WT mice (p < 0.05, VIP score > 1.0). Specifically, eight metabolites were altered at 20w, and 23 metabolites were altered at 40w (p < 0.05, VIP score > 1.0, AUC > 0.8). In addition, 23 metabolites were found to be different in the urine of mice with FD compared with WT mice (p < 0.05, VIP score > 1.0). Especially, nine metabolites were altered at 20w and 19 metabolites were altered at 40w (p < 0.05, VIP score, > 1.0, AUC > 0.8). Pathway analysis revealed that tryptophan, arachidonic acid metabolism, arginine biosynthesis, and glutamine and glutamate metabolism were affected in the serum of mice with FD. Our findings could be differentiated from previous human studies in that we performed systemic metabolomics analyses including OAs, kynurenine, or NSs in serum and urine at early and late age points to find the potential biomarker beyond Gb3 and lyso-Gb3 analogs.
Regarding the statistical analysis conducted on the serum of mice with FD, a total of eight metabolites including four AAs and four NSs were significantly altered in mice with FD at 20w (p < 0.05, VIP score > 1.0, AUC > 0.8). The kidneys are involved in the synthesis and exchange of several AAs between different organs [26]. Glutathione metabolism is considered the most important antioxidant activity and various physiologic functions, including detoxification, modulation of redox-regulated signal transduction, storage and transport of cysteine, regulation of cell proliferation and immune responses, and metabolism of leukotrienes and prostaglandins [27]. Under oxidative stress, glutathione is synthesized from L-glutamate, L-cysteine, and glycine through glutathione synthetic enzymes and is oxidized to GSSG. A recent study showed that at the time of diagnosis, patients with FD exhibit altered glutathione metabolism and elevated glutathione levels [28]. Similar to the results of this previous study, we found that the levels of glutathione and GSSG were higher in mice with FD than in WT mice at 20w, indicating that an increase in glutathione metabolism could be an adaptive response to oxidative stress. However, the GSSG level in mice with FD was found to be similar to that of WT mice at 40w. These results may be indicative of decreased activity of glutathione peroxidase and GSSG reductase as well as an increase in ROS production, which could lead to a decrease in the GSSG level with the progression of FD [29]. Moreover, increased concentrations of ROS reduce the bioavailability of nitric oxide (NO) due to its rapid oxidative degradation and dysfunction of endothelial NO synthase, the enzyme responsible for producing NO in endothelial cells [30]. In particular, the accumulation of Gb3 in the vascular endothelium of FD is linked to an increase in the production of ROS [15,31]. Citrulline is a non-essential AA that is produced as an intermediate product of the urea cycle. It is also a byproduct of NO synthesis by endothelial NO synthase in the presence of tetrahydrobiopterin as a cofactor [30]. In this light, the increase in citrulline observed in mice with FD at 40w may be due to an increase in the byproducts of NO production, which could serve as a defense mechanism against oxidative stress.
Furthermore, the kynurenine pathway of tryptophan catabolism is known to be dysregulated in conditions associated with inflammation, such as cardiovascular disease, atherosclerosis, and chronic kidney disease [32,33]. A recent deep plasma targeted proteomic profiling and network analysis study also showed that several inflammation-related proteins are dysregulated in patients with FD [34]. Notably, our findings demonstrated that the levels of metabolites such as tryptophan, picolinic acid, and serotonin, which are involved in the kynurenine pathway, are altered in the serum of mice with FD at 40w. In addition, our final pathway analyses of the serum suggest that various pathways including those involved in tryptophan metabolism, arachidonic acid metabolism, arginine biosynthesis, proline metabolism, and glutamine and glutamate metabolism can be affected in mice with FD with high oxidative stress, inflammation, and immune intolerance.
Urinary NSs and deoxynucleosides are metabolites that result from oxidative DNA damage or RNA degradation [35]. Degraded normal NSs can be salvaged as nucleotide components, while the modified NSs are excreted into the urine. Therefore, identifying metabolic changes of metabolites in NSs in the urine in FD is important. Metabolites associated with NSs metabolism showed higher concentrations in the urine of mice with FD than those with WT mice at 40w. This finding suggests that urinary excretion of modified NSs were increased in the later stages of FD in mice. Interestingly, picolinic acid, which is another tryptophan metabolite of the downstream kynurenine pathway, was significantly increased in the urine at 40w in mice with FD. Picolinic acid has been implicated in neuroprotection and immune modulation [36,37]. Furthermore, the levels of α-ketoglutaric acid, which is a dicarboxylate OA and a key intermediate metabolite in the TCA cycle involved in maintaining mitochondrial metabolism homeostasis [38,39], were found to be increased in the urine of mice with FD collected at 20w compared with the levels in the urine of WT mice. This increased level of α-ketoglutaric acid in the urine may indicate mitochondrial disorders. Although substantial differences in AAs were observed at 40w between mice with FD and WT mice, there were no statistically significant differences in the urine collected at 20w and 40w.
In conclusion, this systematic metabolomics study revealed that 27 metabolites in the serum and 23 metabolites in the urine were associated with FD in mice at 20w and 40w. These altered metabolites were associated with oxidative stress, inflammation, NO biosynthesis, and immune regulation. Through a better understanding of the differential metabolic pathways involved in FD mice, we may be able to identify potential biomarkers to predict FD progression and provide insight into disease prognosis by assessing metabolic responses to treatment in early and late stages. Furthermore, these findings may provide novel therapeutic strategies to protect and reverse end-organ damage. Future studies are needed to validate the alteration of metabolic profiles related to oxidative stress, inflammation, and NO biosynthesis in the FD patients cohort.

Notes

Conflicts of interest

All authors have no conflicts of interest to declare.

Funding

This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2018R1D1A1B07042999), the NRF funded by the Korea government (MSIT) (NRF-2019R1A2C2086276), and the NRF grant funded by the Korea government (MSIT) (RS-2023-00217317).

Data sharing statement

The data presented in this study are available from the corresponding author upon reasonable request.

Authors’ contributions

Conceptualization: CSK, MJP, SWK

Data curation, Formal analysis: CSK, SO, MJ, BC, TRO, SHS, HSC, EHB, SKM

Funding acquisition: CSK, SWK

Methodology: CSK, SO, MJ, BC, MJP

Supervision: MPJ, SWK

Writing–original draft: CSK, SO

Writing–review & editing: MJ, BC, MJP, SWK

All authors read and approved the final manuscript.

Figure 1.
Star symbol plots in the serum of WT and FD groups at 20 weeks (A–E) and 40 weeks (F–J). (A, F) Organic acids, (B, G) fatty acids, (C, H) amino acids, (D, I) kynurenine pathway metabolites, and (E, J) nucleosides. The numbers on rays correspond to those in Supplementary Table 1 (available online).
FD, Fabry disease; WT, wild-type.
*p < 0.05.
j-krcp-23-218f1.jpg
Figure 2.

Star symbol plots in the urine of WT and FD groups at 20 weeks (A–E) and 40 weeks (F–J). (A, F) Organic acids, (B, G) fatty acids, (C, H) amino acids, (D, I) kynurenine pathway metabolites, and (E, J) nucleosides.

The numbers on rays correspond to those in Supplementary Table 2 (available online).
FD, Fabry disease; WT, wild-type.
*p < 0.05.
j-krcp-23-218f2.jpg
Figure 3.

Multivariate analyses of serum from 20w and 40w FD model mice performed with the 27 selected metabolites.

PLS-DA scores plot (A: 20w, D: 40w), VIP scores of PLS-DA (B: 20w, E: 40w), and pathway analysis (C: 20w, F: 40w).
FD, Fabry disease; GSSG, glutathione disulfide; MTA, 5’-deoxy-5’-methylthioadenosine; PLS-DA, partial least squares discriminant analysis; VIP, variable importance in projection; WT, wild-type; 20w, 20-week-old; 40w, 40-week-old.
j-krcp-23-218f3.jpg
Figure 4.

Multivariate analyses of urine from 20w and 40w FD model mice performed with the 23 selected metabolites.

PLS-DA scores plot (A: 20w, D: 40w), VIP scores of PLS-DA (B: 20w, E: 40w), and pathway analysis (C: 20w, F: 40w).
FD, Fabry disease; PLS-DA, partial least squares discriminant analysis; TCA, tricarboxylic acid cycle; VIP, variable importance in projection; WT, wild-type; 20w, 20-week-old; 40w, 40-week-old.
j-krcp-23-218f4.jpg
Table 1.
Levels of the 27 selected metabolites, p-value, VIP scores, and AUC values in serum of 20w and 40w FD mice models
No. Metabolite 20w serum concentration (ng/μL, mean ± SD)
40w serum concentration (ng/μL, mean ± SD)
WT FD Normalized valuea p-valueb VIPc AUCd WT FD Normalized valuea p-valueb VIPc AUCd
1 Glutamic acid 6.5 ± 0.8 10.8 ± 3.4 1.65 0.008 1.97 0.93 12.1 ± 3.5 24.5 ± 15.9 2.03 0.03 1.34 0.86
2 1-Methylguanosine 0.014 ± 0.002 0.009 ± 0.003 0.67 0.04 1.97 0.86 0.010 ± 0.004 0.004 ± 0.001 0.45 0.007 1.69 0.92
3 MTA 0.013 ± 0.004 0.008 ± 0.003 0.61 0.04 1.74 0.86 0.015 ± 0.006 0.004 ± 0.001 0.26 0.001 2.04 1.00
4 N4-Acetylcytidine 0.19 ± 0.04 0.11 ± 0.03 0.59 0.005 2.17 0.95 2.9 ± 1.0 1.6 ± 0.2 0.54 0.002 1.77 0.96
5 Isoleucine 12.5 ± 1.9 8.8 ± 0.9 0.70 0.005 2.18 0.95 8.6 ± 0.7 9.5 ± 2.2 1.11 0.62 0.55 0.59
6 Aspartic acid 1.8 ± 0.2 2.7 ± 0.6 1.52 0.008 2.04 0.93 4.0 ± 0.8 5.1 ± 2.7 1.29 0.71 0.60 0.57
7 GSSG 1.8 ± 0.2 5.1 ± 2.8 2.81 0.001 2.15 1.00 0.77 ± 0.24 0.77 ± 0.38 1.01 0.62 0.16 0.59
8 Cytidine 0.43 ± 0.08 0.32 ± 0.12 0.73 0.04 1.57 0.86 0.34 ± 0.15 0.32 ± 0.19 0.96 0.21 0.70 0.71
9 cis-Aconitic acid 0.47 ± 0.02 0.46 ± 0.01 0.98 0.37 0.81 0.67 0.10 ± 0.02 0.13 ± 0.02 1.28 0.02 1.41 0.88
10 Arachidonic acid 120.4 ± 54.4 141.2 ± 29.3 1.17 0.45 0.83 0.64 149.1 ± 41.9 88.2 ± 17.2 0.59 0.01 1.67 0.90
11 Docosahexaenoic acid 52.8 ± 14.7 61.9 ± 13.4 1.17 0.30 0.86 0.69 81.8 ± 13.2 61.1 ± 11.1 0.75 0.01 1.49 0.90
12 Docosatetraenoic acid 17.3 ± 6.9 20.9 ± 5.0 1.21 0.10 1.03 0.79 28.1 ± 4.8 21.5 ± 2.7 0.76 0.04 1.55 0.84
13 Tryptophan 8.5 ± 1.8 10.9 ± 3.2 1.27 0.30 1.11 0.69 9.7 ± 3.2 14.4 ± 2.2 1.48 0.04 1.50 0.84
14 Pipecolic acid 0.60 ± 0.38 1.1 ± 1.0 1.75 0.45 0.82 0.64 0.61 ± 0.22 1.2 ± 0.5 1.91 0.01 1.43 0.90
15 α-Aminoadipic acid 1.0 ± 1.0 0.93 ± 0.54 0.91 0.84 0.15 0.55 1.0 ± 0.6 2.5 ± 1.2 2.52 0.02 1.43 0.88
16 Threonine 13.2 ± 3.0 14.2 ± 4.7 1.07 0.95 0.21 0.52 20.8 ± 3.3 27.2 ± 5.6 1.31 0.04 1.31 0.84
17 Creatine 17.3 ± 3.5 22.2 ± 4.2 1.28 0.05 1.42 0.83 34.4 ± 6.6 47.7 ± 9.9 1.39 0.01 1.47 0.90
18 Citrulline 7.5 ± 1.6 9.1 ± 3.0 1.22 0.37 0.77 0.67 9.9 ± 1.9 15.8 ± 3.5 1.60 0.004 1.72 0.94
19 Lysine 20.8 ± 4.5 25.6 ± 12.2 1.23 0.84 0.47 0.55 68.3 ± 19.2 99.2 ± 22.5 1.45 0.03 1.46 0.86
20 Picolinic acid 0.017 ± 0.005 0.017 ± 0.004 0.96 0.84 0.13 0.55 0.030 ± 0.005 0.037 ± 0.005 1.24 0.03 1.44 0.86
21 5-Hydroxyindoleacetic acid 0.027 ± 0.013 0.032 ± 0.007 1.18 0.73 0.76 0.57 0.071 ± 0.018 0.12 ± 0.05 1.64 0.03 1.34 0.86
22 Serotonin 0.55 ± 0.05 0.66 ± 0.12 1.21 0.23 1.50 0.71 2.2 ± 0.4 1.6 ± 0.17 0.74 0.007 1.63 0.92
23 5,6-Dihydrouridine 0.22 ± 0.05 0.22 ± 0.06 0.98 >0.99 0.14 0.50 0.19 ± 0.07 0.08 ± 0.02 0.42 0.001 1.85 0.98
24 Pseudouridine 0.74 ± 0.19 0.59 ± 0.18 0.80 0.10 1.11 0.79 0.51 ± 0.14 0.33 ± 0.05 0.64 0.01 1.59 0.90
25 N2-Methylguanosine 0.014 ± 0.001 0.010 ± 0.004 0.72 0.07 1.76 0.81 0.011 ± 0.005 0.006 ± 0.001 0.52 0.01 1.60 0.90
26 5-Methylcytidine 0.016 ± 0.003 0.012 ± 0.003 0.72 0.07 1.70 0.81 0.018 ± 0.005 0.012 ± 0.001 0.68 0.01 1.55 0.90
27 N2,N2-Dimethylguanosine 0.005 ± 0.001 0.004 ± 0.002 0.80 0.14 1.25 0.76 0.007 ± 0.002 0.005 ± 0.001 0.64 0.01 1.58 0.90

AUC, area under the curve; FD, Fabry disease; GSSG, glutathione disulfide; MTA, 5’-deoxy-5’-methylthioadenosine; SD, standard deviation; VIP, variable importance in projection; WT, wild-type; 20w, 20-week-old; 40w, 40-week-old.

a Values normalized to the corresponding mean value of the WT concentration.

b p-value calculated by Wilcoxon rank-sum test.

c VIP score of partial least squares discriminant analysis.

d Values of univariate area under the receiver operating characteristic curve.

Table 2.
Levels of the 23 selected metabolites, p-value, VIP scores, and AUC values in urine of 20w and 40w FD models
No. Metabolite 20w urine concentration (pg/ng creatinine, mean ± SD)
40w urine concentration (pg/ng creatinine, mean ± SD)
WT FD Normalized valuea p-valueb VIPc AUCd WT FD Normalized valuea p-valueb VIPc AUCd
1 2-Hydroxybutyric acid 7.9 ± 0.9 10.4 ± 1.3 1.33 0.004 2.51 0.97 4.2 ± 0.8 6.8 ± 1.8 1.62 0.007 1.63 0.92
2 3-Hydroxypropionic acid 86.1 ± 12.1 135.8 ± 46.6 1.58 0.04 2.02 0.86 54.7 ± 11.2 71.7 ± 12.8 1.31 0.04 1.29 0.84
3 1-Methyladenosine 1.2 ± 1.1 0.22 ± 0.14 0.18 0.02 2.27 0.92 0.33 ± 0.21 0.14 ± 0.02 0.43 0.007 1.52 0.92
4 N6-Methyladenosine 1.2 ± 1.1 0.23 ± 0.15 0.18 0.02 2.27 0.92 0.31 ± 0.20 0.14 ± 0.03 0.45 0.02 1.43 0.88
5 Kynurenic acid 30.3 ± 4.0 25.9 ± 3.0 0.85 0.04 1.77 0.86 7.4 ± 2.0 9.8 ± 1.9 1.33 0.04 1.35 0.84
6 Pyruvic acid 92.6 ± 12.5 171.4 ± 55.0 1.85 0.002 2.58 1.00 69.3 ± 25.2 66.5 ± 23.9 0.96 >0.99 0.09 0.51
7 Glycolic acid 400.2 ± 46.3 512.7 ± 59.2 1.28 0.004 2.40 0.97 258.4 ± 66.3 294.1 ± 39.8 1.14 0.26 0.82 0.69
8 α-Ketoglutaric acid 3,483.1 ± 1,327.8 6,723.5 ± 3,410.5 1.93 0.04 1.86 0.86 3,331.2 ± 861.8 2,876.8 ± 1,536.7 0.86 0.62 0.66 0.59
9 4-Hydroxyphenylpyruvic acid 387.6 ± 148.6 786.4 ± 317.7 2.03 0.04 1.72 0.86 606.9 ± 323.3 648.5 ± 322.7 1.07 0.62 0.18 0.59
10 Phenylacetic acid 27.6 ± 15.6 27.7 ± 10.2 1.00 0.82 0.27 0.56 8.3 ± 1.6 12.8 ± 3.3 1.54 0.03 1.54 0.86
11 Malonic acid 80.4 ± 9.3 79.6 ± 14.7 0.99 0.94 0.21 0.53 41.6 ± 12.6 70.8 ± 8.3 1.70 0.004 1.87 0.94
12 Ethylmalonic acid 25.0 ± 4.6 23.7 ± 2.9 0.95 0.59 0.43 0.61 10.4 ± 2.5 21.1 ± 7.3 2.04 0.002 1.80 0.96
13 Isovaleryglycine 1,206.3 ± 176.1 1,198.5 ± 211.9 0.99 0.94 0.12 0.53 612.8 ± 178.5 913.1 ± 406.4 1.49 0.02 1.21 0.88
14 Methylsuccinic acid 41.7 ± 5.5 40.6 ± 3.0 0.97 0.82 0.30 0.56 22.5 ± 4.1 29.2 ± 2.3 1.30 0.004 1.60 0.94
15 Hexanoylglycine 840.1 ± 162.8 928.0 ± 261.9 1.10 0.59 0.42 0.61 329.4 ± 190.9 741.3 ± 488.5 2.25 0.04 1.46 0.84
16 Adipic acid 259.9 ± 36.1 229.0 ± 44.5 0.88 0.31 1.21 0.69 101.5 ± 45.5 191.7 ± 72.4 1.89 0.02 1.54 0.88
17 3-Methyladipic acid 7.4 ± 0.7 7.0 ± 1.5 0.95 0.94 0.70 0.53 3.5 ± 0.5 7.4 ± 2.5 2.09 0.001 1.89 1.00
18 Vanillic acid 11.2 ± 3.9 7.8 ± 3.4 0.70 0.24 1.26 0.72 9.6 ± 3.0 5.3 ± 3.1 0.55 0.03 1.39 0.86
19 3-Indolecarboxylic acid 18.3 ± 2.6 16.5 ± 2.3 0.90 0.49 1.09 0.64 9.0 ± 1.2 11.1 ± 1.2 1.24 0.004 1.55 0.94
20 Vanillylmandelic acid 1.3 ± 0.2 1.2 ± 0.2 0.94 0.59 0.76 0.61 0.61 ± 0.13 0.82 ± 0.18 1.35 0.02 1.34 0.88
21 Picolinic acid 5.0 ± 4.5 5.0 ± 7.7 1.01 0.59 0.56 0.61 0.52 ± 0.21 2.6 ± 2.8 5.06 0.03 1.43 0.86
22 5-Hydroxyindoleacetic acid 34.9 ± 7.8 32.3 ± 2.5 0.93 0.59 0.57 0.61 12.9 ± 2.3 16.3 ± 1.6 1.27 0.01 1.49 0.90
23 Xanthosine 6.3 ± 3.1 8.4 ± 3.4 1.33 0.31 0.95 0.69 6.1 ± 1.5 10.7 ± 2.5 1.75 0.004 1.77 0.94

AUC, area under the curve; FD, Fabry disease; SD, standard deviation; VIP, variable importance in projection; WT, wild-type; 20w, 20-week-old; 40w, 40-week-old.

a Values normalized to the corresponding mean value of the WT concentration.

b p-value calculated by Wilcoxon rank-sum test.

c VIP score of partial least squares discriminant analysis.

d Values of univariate area under the receiver operating characteristic curve.

References

1. Zarate YA, Hopkin RJ. Fabry’s disease. Lancet 2008;372:1427–1435.
crossref pmid
2. Eng CM, Fletcher J, Wilcox WR, et al. Fabry disease: baseline medical characteristics of a cohort of 1765 males and females in the Fabry Registry. J Inherit Metab Dis 2007;30:184–192.
crossref pmid pdf
3. Boutin M, Lavoie P, Abaoui M, Auray-Blais C. Tandem mass spectrometry quantitation of lyso-gb3 and six related analogs in plasma for Fabry disease patients. Curr Protoc Hum Genet 2016;90:17.23.1–17.23.9.
crossref pmid pdf
4. Lavoie P, Boutin M, Auray-Blais C. Multiplex analysis of novel urinary lyso-Gb3-related biomarkers for Fabry disease by tandem mass spectrometry. Anal Chem 2013;85:1743–1752.
crossref pmid
5. Auray-Blais C, Boutin M. Novel gb(3) isoforms detected in urine of Fabry disease patients: a metabolomic study. Curr Med Chem 2012;19:3241–3252.
crossref pmid
6. Manwaring V, Boutin M, Auray-Blais C. A metabolomic study to identify new globotriaosylceramide-related biomarkers in the plasma of Fabry disease patients. Anal Chem 2013;85:9039–9048.
crossref pmid
7. Hu L, Liu J, Zhang W, et al. Functional metabolomics decipher biochemical functions and associated mechanisms underlie small-molecule metabolism. Mass Spectrom Rev 2020;39:417–433.
crossref pmid pdf
8. Boutin M, Auray-Blais C. Metabolomic discovery of novel urinary galabiosylceramide analogs as Fabry disease biomarkers. J Am Soc Mass Spectrom 2015;26:499–510.
crossref pmid pdf
9. Beneduci A, Cuccurullo M, Pontoni G, Chidichimo G, Capasso G. Perspectives of 1H-NMR-based urinary metabonomics in Fabry disease. J Nephrol 2010;23 Suppl 16:S213–S220.
pmid
10. Auray-Blais C, Boutin M, Gagnon R, Dupont FO, Lavoie P, Clarke JT. Urinary globotriaosylsphingosine-related biomarkers for Fabry disease targeted by metabolomics. Anal Chem 2012;84:2745–2753.
crossref pmid
11. Ducatez F, Mauhin W, Boullier A, et al. Parsing Fabry disease metabolic plasticity using metabolomics. J Pers Med 2021;11:898.
crossref pmid pmc
12. Jehn U, Bayraktar S, Pollmann S, et al. α-Galactosidase a deficiency in Fabry disease leads to extensive dysregulated cellular signaling pathways in human podocytes. Int J Mol Sci 2021;22:11339.
crossref pmid pmc
13. Müller KB, Galdieri LC, Pereira VG, Martins AM, D’Almeida V. Evaluation of oxidative stress markers and cardiovascular risk factors in Fabry disease patients. Genet Mol Biol 2012;35:418–423.
crossref pmid pmc
14. Ravarotto V, Simioni F, Carraro G, Bertoldi G, Pagnin E, Calò LA. Oxidative stress and cardiovascular-renal damage in Fabry disease: is there room for a pathophysiological involvement? J Clin Med 2018;7:409.
crossref pmid pmc
15. Shen JS, Meng XL, Moore DF, et al. Globotriaosylceramide induces oxidative stress and up-regulates cell adhesion molecule expression in Fabry disease endothelial cells. Mol Genet Metab 2008;95:163–168.
crossref pmid pmc
16. Ohshima T, Murray GJ, Swaim WD, et al. alpha-Galactosidase A deficient mice: a model of Fabry disease. Proc Natl Acad Sci U S A 1997;94:2540–2544.
pmid pmc
17. Hur H, Paik MJ, Xuan Y, et al. Quantitative measurement of organic acids in tissues from gastric cancer patients indicates increased glucose metabolism in gastric cancer. PLoS One 2014;9:e98581.
crossref pmid pmc
18. Paik MJ, Cho EY, Kim H, et al. Simultaneous clinical monitoring of lactic acid, pyruvic acid and ketone bodies in plasma as methoxime/tert-butyldimethylsilyl derivatives by gas chromatography-mass spectrometry in selected ion monitoring mode. Biomed Chromatogr 2008;22:450–453.
crossref pmid
19. Paik MJ, Kim KR. Sequential ethoxycarbonylation, methoximation and tert-butyldimethylsilylation for simultaneous determination of amino acids and carboxylic acids by dual-column gas chromatography. J Chromatogr A 2004;1034:13–23.
crossref pmid
20. Seo C, Kim SH, Lee HS, et al. Metabolomic study on bleomycin and polyhexamethylene guanidine phosphate-induced pulmonary fibrosis mice models. Metabolomics 2019;15:111.
crossref pmid pdf
21. Kim Y, Kim SH, Oh SJ, et al. Metabolomic analysis of organic acids, amino acids, and fatty acids in plasma of Hanwoo beef on a high-protein diet. Metabolomics 2020;16:114.
crossref pmid pdf
22. Seo C, Na M, Jang J, et al. Monitoring of altered amino acid metabolic pattern in rat urine following intraperitoneal injection with γ-hydroxybutyric acid. Metabolomics 2018;14:111.
crossref pmid pdf
23. Choi RY, Ji M, Lee MK, Paik MJ. Metabolomics study of serum from a chronic alcohol-fed rat model following administration of defatted Tenebrio molitor larva fermentation extract. Metabolites 2020;10:436.
crossref pmid pmc
24. Ji M, Jo Y, Choi SJ, et al. Plasma metabolomics and machine learning-driven novel diagnostic signature for non-alcoholic steatohepatitis. Biomedicines 2022;10:1669.
crossref pmid pmc
25. Kim HY, Lee HS, Kim IH, et al. Comprehensive targeted metabolomic study in the lung, plasma, and urine of PPE/LPS-induced COPD mice model. Int J Mol Sci 2022;23:2748.
crossref pmid pmc
26. van de Poll MC, Soeters PB, Deutz NE, Fearon KC, Dejong CH. Renal metabolism of amino acids: its role in interorgan amino acid exchange. Am J Clin Nutr 2004;79:185–197.
crossref pmid
27. Santangelo F, Witko-Sarsat V, Drüeke T, Descamps-Latscha B. Restoring glutathione as a therapeutic strategy in chronic kidney disease. Nephrol Dial Transplant 2004;19:1951–1955.
crossref pmid
28. Biancini GB, Jacques CE, Hammerschmidt T, et al. Biomolecules damage and redox status abnormalities in Fabry patients before and during enzyme replacement therapy. Clin Chim Acta 2016;461:41–46.
crossref pmid
29. Kim JW, Kim HW, Nam SA, et al. Human kidney organoids reveal the role of glutathione in Fabry disease. Exp Mol Med 2021;53:1580–1591.
crossref pmid pmc pdf
30. Förstermann U. Nitric oxide and oxidative stress in vascular disease. Pflugers Arch 2010;459:923–939.
crossref pmid pdf
31. Biancini GB, Vanzin CS, Rodrigues DB, et al. Globotriaosylceramide is correlated with oxidative stress and inflammation in Fabry patients treated with enzyme replacement therapy. Biochim Biophys Acta 2012;1822:226–232.
crossref pmid
32. Wang Q, Liu D, Song P, Zou MH. Tryptophan-kynurenine pathway is dysregulated in inflammation, and immune activation. Front Biosci (Landmark Ed) 2015;20:1116–1143.
crossref pmid pmc
33. Cheng Y, Li Y, Benkowitz P, Lamina C, Köttgen A, Sekula P. The relationship between blood metabolites of the tryptophan pathway and kidney function: a bidirectional Mendelian randomization analysis. Sci Rep 2020;10:12675.
crossref pmid pmc pdf
34. Tebani A, Barbey F, Dormond O, et al. Deep next-generation proteomics and network analysis reveal systemic and tissue-specific patterns in Fabry disease. Transl Res 2023;258:47–59.
crossref pmid
35. Patejko M, Struck-Lewicka W, Siluk D, Waszczuk-Jankowska M, Markuszewski MJ. Urinary nucleosides and deoxynucleosides. Adv Clin Chem 2018;83:1–51.
crossref pmid
36. Ding K, McGee-Lawrence ME, Kaiser H, et al. Picolinic acid, a tryptophan oxidation product, does not impact bone mineral density but increases marrow adiposity. Exp Gerontol 2020;133:110885.
crossref pmid pmc
37. Rapisarda A, Pastorino S, Massazza S, Varesio L, Bosco MC. Antagonistic effect of picolinic acid and interferon-gamma on macrophage inflammatory protein-1alpha/beta production. Cell Immunol 2002;220:70–80.
pmid
38. Yuan Y, Zhu C, Wang Y, et al. α-Ketoglutaric acid ameliorates hyperglycemia in diabetes by inhibiting hepatic gluconeogenesis via serpina1e signaling. Sci Adv 2022;8:eabn2879.
crossref pmid pmc
39. Yang Y, He P, Hou Y, Liu Z, Zhang X, Li N. Osmundacetone modulates mitochondrial metabolism in non-small cell lung cancer cells by hijacking the glutamine/glutamate/α-KG metabolic axis. Phytomedicine 2022;100:154075.
crossref pmid


ABOUT
BROWSE ARTICLES
EDITORIAL POLICY
FOR CONTRIBUTORS
Editorial Office
#301, (Miseung Bldg.) 23, Apgujenog-ro 30-gil, Gangnam-gu, Seoul 06022, Korea
Tel: +82-2-3486-8736    Fax: +82-2-3486-8737    E-mail: registry@ksn.or.kr                

Copyright © 2024 by The Korean Society of Nephrology.

Developed in M2PI

Close layer