You're viewing our updated article page. If you need more time to adjust, you can return to the old layout.

ORIGINAL RESEARCH article

Front. Med., 09 February 2026

Sec. Hepatobiliary Diseases

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1646323

Serum-urine metabolic integration via UPLC-QTOF/MS uncovers shared pathway biomarkers for cirrhosis diagnosis

  • 1. First Affiliated Hospital of Baotou Medical College, Baotou, China

  • 2. Department of Pharmacy, Baotou Medical College, Inner Mongolia University of Science and Technology, Baotou, China

Article metrics

View details

284

Views

28

Downloads

Abstract

Introduction:

Liver cirrhosis is the terminal stage of chronic liver disease, which is marked by high morbidity and mortality in its advanced phases. Although liver biopsy still serves as the gold - standard diagnostic method, the detection of serum and urine metabolites holds great promise for the identification of cirrhosis.

Methods:

Untargeted metabolomics analysis was carried out using ultra - performance liquid chromatography coupled with quadrupole time - of - flight mass spectrometry (UPLC - QTOF/MS). We compared the serum and urine metabolic profiles between 30 healthy individuals and 28 liver cirrhosis patients to screen for biomarkers associated with liver cirrhosis.

Results:

A total of 55 endogenous metabolites showed dysregulation in serum, and 51 did so in urine. Four shared differential metabolites—glycoursodeoxycholic acid, urobilin, glycocholic acid, and urobilinogen—were identified in both biofluids. Pathway enrichment analysis revealed three co - regulated metabolic pathways: tryptophan metabolism, glycerophospholipid metabolism, and porphyrin metabolism (p < 0.05).

Discussion:

This study delineates the distinct metabolic signatures of cirrhosis and proposes a diagnostic strategy based on dual - biofluid analysis. The intersectional biomarkers and pathways elucidate the mechanisms linking bile acid homeostasis and hemoprotein catabolism to cirrhotic progression, offering a noninvasive approach for clinical detection.

1 Introduction

Chronic liver diseases (CLDs) affect approximately 800 million individuals worldwide, leading to roughly 2 million fatalities per year (1, 2). The progression of CLDs can be ascribed to a variety of factors, such as alcohol misuse, obesity and other metabolic disorders, autoimmune hepatitis, and viral hepatitis (3). Cirrhosis represents the end-stage consequence of chronic hepatitis and can advance to diffuse cirrhosis. In patients with cirrhosis, the normal hepatic architecture is supplanted by regenerative nodules, and in severe instances, it may culminate in liver failure. The early stage of cirrhosis is frequently asymptomatic yet potentially reversible (4). Initial assessments for cirrhosis encompass serological assays for viral hepatitis, quantification of ferritin and transferrin saturation, abdominal ultrasonography, a complete blood count, liver function evaluations, and determination of the prothrombin time/international normalized ratio (5). The therapeutic objectives for cirrhosis are to prevent its onset, decompensation, and mortality. Non-selective beta-blockers are frequently employed (6, 7). The management of ascites entails the use of diuretics, salt restriction, and antibiotic therapy. Liver biopsy still serves as the gold-standard method for the detection of cirrhosis. Non-invasive assessments are particularly valuable in identifying either the early or advanced stages of cirrhosis. Serum-based biomarkers for cirrhosis have been established (8).

Metabolomics, a crucial domain in “omics” research, is centered on the high-throughput identification and quantitative assessment of small-molecule metabolites (with a molecular weight < 1,500 Da) in the metabolome (9, 10). The progress of analytical technologies and bioinformatics has elevated metabolomics to a key position as a systems biology tool, propelling its extensive application as an integrated diagnostic strategy in clinical and biomedical research (11, 12). Yang (13) uncovered the distinctive metabolic features of cancer cachexia via serum and urine metabolomics and developed a diagnostic model. Blood, being a vital and easily accessible biofluid, is of great significance in clarifying the pathogenesis and progression of human diseases (14). Urine represents an optimal biofluid for disease research owing to its non-invasive collection and relatively simpler composition compared with other bodily fluids. The convenience of serial sampling facilitates the longitudinal monitoring of disease progression and the assessment of therapeutic response (15). In this study, untargeted metabolomics was utilized to characterize the serum and urine metabolites in patients with cirrhosis. Our findings provide mechanistic underpinnings for non-invasive diagnosis of cirrhosis and identification of therapeutic targets.

2 Materials and methods

2.1 Study participants

We recruited patients with cirrhosis from the First Affiliated Hospital of Baotou Medical College between February and July 2023, and obtained clinical data from 28 patients with cirrhosis. The diagnosis of cirrhosis included liver biopsy, imaging studies, coagulation tests, complete blood count, and complications of decompensated cirrhosis (ascites, gastrointestinal bleeding, sepsis, hepatic encephalopathy, and hepatorenal syndrome). The serum and urine samples of 30 healthy people without disease were obtained at the same time. Inclusion criteria of cirrhosis: (1) the cirrhosis group had clinical symptoms such as loss of appetite, anorexia, fatigue, and discomfort or pain in the liver region on admission; (2) the healthy group did not have any liver disease; (3) complete medical history, aged 18 years or older, and no cognitive impairment. Exclusion criteria include: (1) severe organic diseases; (2) diagnosis of liver cancer or other malignant tumors; (3) recent use of drugs that affect liver function indexes. The Medical Ethics Committee of the First Affiliated Hospital of Baotou Medical College approved this study (Ethical approval number: 2022026). Written informed consent was obtained from all participants. General data (age, gender, height, weight, body mass index, etc.) and laboratory examinations were collected for all subjects. All methods were performed in accordance with the relevant guidelines and regulations.

2.2 Untargeted metabolomics

The 100 μL sample was transferred to an eppendorf tube, and then 400 μL of extract (methanol: acetonitrile = 1:1) was added. The mixture was vortexed and mixed for 30s, and sonicated for 10 min under an ice water bath, and then left at −40 °C for 1 h. The serum samples were centrifuged at 12000 rpm for 15 min at 4 °C, and the supernatant was taken for detection. Equal volumes of supernatant from all samples were mixed to make quality control (QC) samples. The QC samples were scanned at intervals of 10 samples during the collection process. Systematic errors were corrected with the quality gap between QC samples. The target compounds were separated by chromatography using a Vanquish ultra-high performance liquid chromatographer (Thermo Fisher Scientific). The parameters for setting the liquid phase gradient are as follows: 0.0–0.8 min, 2%B; 0.8–2.8 min, 2–70%B; 2.8–5.3 min, 70–90%B; 5.3–5.9 min, 90–100% B; 5.9–7.5 min, 100% B; 7.5–7.6 min, 100–2% B; 7.6–10.0 min, 2% B. The flow rate was 0.3 mL/min. Phase A of liquid chromatography consisted of 25 mmol/L ammonium acetate and 25 mmol/L ammonia, and phase B consisted of acetonitrile. Sample plate temperature: 4 °C, injection volume: 2 μL. Primary and secondary mass spectrometry data were acquired using a mass spectrometer (Orbitrap Exploris 120) controlled by Xcalibur (V4.4, Thermo). After converting the raw data into mzXML format, peak identification, extraction, alignment, and integration were performed. When integrating standardized datasets, probabilistic quotient normalization (PQN) was applied to ensure inter-sample comparability and mitigate scaling variability (16). By comparing the spectral intensities of a sample with those of other samples within its neighborhood, the probability quotient was calculated and the relative intensity was adjusted. The kyoto encyclopedia of genes and genomes (KEGG) database and human metabolome database (HMDB) were used to annotate metabolites, and the molecular weight data (m/z) were matched with the data in the database. Metabolites with a mass difference of less than 10 ppm between the observed and database values were annotated, and the molecular formula of metabolites were further identified and validated through isotope distribution measurements. Serum metabolic profile analysis was performed for all metabolites. Metabolites of variable importance in projection (VIP) > 1 and p < 0.05 were screened to draw the volcano map. The differential metabolites were compared with the information on the HMDB website1 to identify endogenous metabolites that met the requirements, and a heat map of metabolite expression clusters was generated. The endogenous differential metabolites were plotted using the MetaboAnalyst.2 We screened metabolic pathways according to impact > 0.02 and raw p < 0.05, and observed the key metabolites that affect cirrhosis within these metabolic pathways.

2.3 Statistical analysis

IBM SPSS Statistics 26.0 software was used to process the data of patients in this study. The enumeration data were expressed as n/%, and we used the chi-square test for calculations. The measurement data were expressed as x ± s, and we used the t and F test for calculations. To mitigate false positives arising from multiple testing, all statistical tests were corrected for the false discovery rate using the Benjamini-Hochberg procedure.

3 Results

3.1 Changes in serum and urine metabolic profiles

The separation trends between groups were observed using score plots. PCA (Principal Component Analysis) and PLS- DA (Partial Least Squares-Discriminant Analysis) were performed on the metabolites of the healthy control group and the disease group. Meanwhile, 200 times of 7-fold cross-validation was carried out on the PLS-DA results to determine whether the model was overfitted (Table 1).

Table 1

No. Rt(min) m/z Formula Metabolite HMDB ID VIP p value Trend
DIS vs. CON
1 0.79 241.03 C6H12N2O4S2 L-cystine HMDB0000192 1.00 0.01 up
2 4.30 498.29 C26H45NO6S taurodeoxycholic acid HMDB0000896 2.67 0.00 up
3 4.18 448.31 C26H43NO5 glycoursodeoxycholic acid HMDB0000708 1.89 0.00 up
4 4.18 899.63 C26H43NO5 glycohyodeoxycholic acid HMDB0304944 1.54 0.00 up
5 8.46 585.27 C33H36N4O6 bilirubin HMDB0000054 1.23 0.00 up
6 3.17 182.05 C8H9NO4 4-pyridoxic acid HMDB0000017 1.08 0.02 up
7 3.10 188.18 C9H21N3O N1-acetylspermidine HMDB0001276 1.50 0.00 up
8 4.60 591.32 C33H42N4O6 urobilin HMDB0004160 1.69 0.00 up
9 3.9 464.30 C26H43NO6 glycocholic acid HMDB0000138 2.10 0.00 up
10 4.10 498.29 C26H45NO6S chenodeoxycholyltaurine HMDB0242411 2.68 0.00 up
11 4.21 407.28 C24H40O5 allocholic acid HMDB0000505 1.45 0.00 up
12 3.36 318.19 C15H27NO6 2,4-dimethylhexanedioylcarnitine HMDB0241047 1.71 0.00 up
13 4.11 383.33 C27H42O N-[(3a,5b,7b)-7-hydroxy-24-oxo-3-(sulfooxy)cholan-24-yl]-glycine HMDB0002409 1.79 0.00 up
14 3.96 627.38 C33H57O9P PA(10:0/20:4(5Z,8Z,11Z,13E)-OH(15S)) HMDB0262659 1.42 0.00 up
15 5.16 389.26 C24H38O4 12 alpha-hydroxy-3-oxo-5beta-cholan-24-oic acid HMDB0062742 1.19 0.00 up
16 7.25 674.44 C36H68NO8P PC(14:1(9Z)/14:1(9Z)) HMDB0007900 2.37 0.00 up
17 4.16 627.38 C33H57O9P PA(10:0/20:4(5Z,8Z,10E,14Z)-OH(12S)) HMDB0262661 1.63 0.00 up
18 4.29 253.14 C11H18N4O3 histidylvaline HMDB0028898 1.01 0.00 up
19 4.21 625.36 C33H55O9P PA(10:0/20:4(5Z,8Z,11Z,13E) + =O(15)) HMDB0262691 1.38 0.00 up
20 4.18 531.30 C27H48O8S 5b-cyprinol sulfate HMDB0006888 1.98 0.00 up
21 1.25 203.15 C8H18N4O2 asymmetric dimethylarginine HMDB0001539 1.07 0.00 up
22 3.64 448.31 C26H43NO6 glutamic acid HMDB0000148 2.01 0.00 up
23 3.97 589.30 C33H42N4O6 urobilinogen HMDB0004158 1.99 0.00 up
24 3.95 462.27 C26H45NO7S taurocholic acid HMDB0000036 2.83 0.00 up
25 2.26 181.05 C9H10O4 hydroxyphenyllactic acid HMDB0000755 1.05 0.00 up
26 6.14 431.31 C27H42O4 7 alpha-hydroxy-3-oxo-4-cholestenoate HMDB0012458 1.10 0.00 up
27 4.11 471.24 C24H40O7S chenodeoxycholic acid 3-sulfate HMDB0002586 1.64 0.00 up
28 5.07 433.24 C21H39O7P LysoPA(18:2(9Z,12Z)/0:0) HMDB0007856 1.73 0.00 down
29 6.13 436.28 C21H44NO6P LysoPE(P-16:0/0:0) HMDB0011152 1.73 0.00 down
30 5.78 500.28 C25H44NO7P LysoPE(20:4(8Z,11Z,14Z,17Z)/0:0) HMDB0011518 1.01 0.00 down
31 3.53 277.16 C15H22N2O3 phenylalanylisoleucine HMDB0028998 1.91 0.00 down
32 5.72 524.28 C27H44NO7P LysoPE(22:6(4Z,7Z,10Z,13Z,16Z,19Z)/0:0) HMDB0011526 1.02 0.00 down
33 0.73 146.05 C5H9NO4 L-4-hydroxyglutamate semialdehyde HMDB0006556 1.62 0.00 down
34 3.29 160.08 C10H12N2O serotonin HMDB0000259 1.22 0.00 down
35 5.25 518.32 C26H48NO7P LysoPC(18:3(9Z,12Z,15Z)/0:0) HMDB0010388 1.50 0.00 down
36 5.82 570.35 C30H52NO7P LysoPC(22:5(7Z,10Z,13Z,16Z,19Z)/0:0) HMDB0010403 1.52 0.00 down
37 5.59 546.34 C28H52NO7P LysoPC(20:3(8Z,11Z,14Z)/0:0) HMDB0010394 1.68 0.00 down
38 5.05 313.24 C18H34O4 octadecanedioic acid HMDB0000782 1.74 0.00 down
39 0.92 124.01 C2H7NO3S taurine HMDB0000251 1.06 0.00 down
40 3.59 313.15 C18H20N2O3 phenylalanylphenylalanine HMDB0013302 1.82 0.00 down
41 0.78 280.09 C8H21NO6P+ glycerophosphocholine HMDB0000086 2.23 0.00 down
42 5.18 468.31 C22H46NO7P LysoPC(14:0/0:0) HMDB0010379 1.73 0.00 down
43 0.86 132.08 C4H6N2O2 dihydrouracil HMDB0000076 1.02 0.00 down
44 0.72 132.03 C4H7NO4 L-aspartic acid HMDB0000191 1.78 0.00 down
45 6.04 496.34 C24H50NO7P LysoPC(16:0/0:0) HMDB0010382 1.68 0.00 down
46 6.42 480.34 C24H50NO6P LysoPC(P-16:0/0:0) HMDB0010407 1.68 0.00 down
47 3.82 120.08 C8H11NO 2-hydroxyphenethylamine HMDB0001065 1.34 0.00 down
48 6.20 303.23 C20H32O3 8-HETE HMDB0004679 2.09 0.00 down
49 7.75 304.24 C20H32O2 arachidonic acid HMDB0001043 1.23 0.00 down
50 9.39 546.35 C26H54NO7P LysoPC(0:0/18:0) HMDB0011128 1.80 0.00 down
51 3.76 199.01 C8H8O4S 4-vinylphenol sulfate HMDB0062775 1.77 0.00 down
52 3.59 120.08 C10H12N2O3 kynurenine HMDB0000684 1.44 0.00 down
53 3.93 308.19 C17H25NO4 4-phenylbutanoylcarnitine HMDB0241867 1.82 0.00 down
54 4.82 267.12 C14H20O5 3-carboxy-4-methyl-5-pentyl-2-furanpropanoic acid HMDB0061643 1.77 0.00 down
55 3.50 454.29 C21H44NO7P LysoPE(16:0/0:0) HMDB0011503 1.93 0.00 down

Differential metabolites in the serum of patients with cirrhosis.

The PCA results showed that there was a significant deviation (8.4%) between the SCON group and the SDIS group along the PCA2 axis (Figure 1A). However, the difference between the UCON group and the UDIS group was not obvious (Figure 1B). The PLS-DA results indicated that there was a significant deviation (21.79%) between the SCON group and the SDIS group along the PC1 axis (Figure 2A). For the urine samples, there was also a significant deviation (12.49%) between the UCON group and the UDIS group along the PC1 axis (Figure 2C). Based on the permutation test results in panels B and D of Figure 2, the R2 and Q2 values of the original model are significantly higher than the distributions of the R2 and Q2 values of the permuted models. This indicates that there is a significant correlation between the grouping of samples (DIS vs. CON) and metabolites. It further supports our hypothesis that there are unique metabolic characteristics in cirrhosis patients (Figures 2B,D).

Figure 1

Two PCA scatter plots labeled A and B show different group distributions. Plot A compares groups SCON (red) and SDIS (blue) with PCA1 at eighty-five point fifty-three percent and PCA2 at eight point four percent. Plot B compares UCON (red) and UDIS (blue) with PCA1 at seventy-two point forty-five percent and PCA2 at seven point thirty-one percent. Each plot includes ellipses encapsulating data points, with a legend indicating group colors.

PCA score plots of serum and urine metabolomic profiles in patients with cirrhosis: (A) serum; (B) urine.

Figure 2

Four-panel image displaying principal component analysis (PCA) and predictive plots. Panel A shows PCA of classes SCON and SDIS, with red and blue data points. Panel B illustrates correlation versus value for the same classes, featuring a linear trend. Panel C presents PCA of classes UCON and UDIS with overlapping clusters. Panel D displays their correlation versus value, with a similar linear pattern. Each plot has axes labeled with respective PCA components, correlation, and value, along with R² and Q² statistics.

PLS-DA analysis of metabolomics in patients with cirrhosis: (A) PLS-DA score plot of serum samples; (B) Permutation test for the PLS-DA model of serum samples, R2 = 0.726, Q2 = −0.349; (C) PLS-DA score plot of urine samples; (D) Permutation test for the PLS-DA model of urine samples, R2 = 0.8041, Q2 = −0.331.

3.2 Screening of differential metabolites in serum and urine

V-plots score graph was constructed with log2(FC) on the x-axis and -log10(P) on the y-axis for all metabolites to visualize the overall distribution of differential metabolites. Potential differential metabolites associated with the disease were screened with VIP > 1 and p < 0.05. Compared with the healthy control group, in the serum of the liver cirrhosis group, the contents of 94 metabolites were up-regulated, and those of 123 metabolites were down-regulated (Figure 3A). In urine, the contents of 97 metabolites were up-regulated, and those of 98 metabolites were down-regulated (Figure 3B). The metabolites obtained from the V-plots score graph were further screened on the HMDB website (see Footnote 1) to identify potential metabolic biomarkers associated with the disease. Heat maps of metabolite contents showed that 55 endogenous metabolites in the serum of liver cirrhosis patients were affected (Figure 4), and 51 endogenous metabolites in urine were affected (Figure 5). Among them, 4 metabolites were differentially affected in both serum and urine of liver cirrhosis patients: glycoursodeoxycholic acid, urobilin, glycocholic acid, and urobilinogen. The contents of these four metabolites were significantly up-regulated in the body fluid samples of liver cirrhosis patients (Table 2).

Figure 3

Two scatter plots, labeled A and B, show data comparisons. Plot A, "SDISVSSCON," displays points where genes are significantly upregulated (red), downregulated (blue), or show no difference (gray). Plot B, "UDISVSUCON," uses the same color scheme. Both plots depict log two fold change on the x-axis and negative log ten p-value on the y-axis, indicating statistical significance.

Metabolomics volcano map of patients with cirrhosis: (A) serum samples and (B) urine samples.

Figure 4

Heatmap illustrating metabolite concentrations across 60 samples (SCON and SDS groups), with values ranging from -2 to 2, indicated by colors from blue to red. Metabolites are listed on the right, samples at the bottom. Dendrograms on the left show hierarchical clustering.

Heat map of metabolite expression in serum of patients with cirrhosis.

Figure 5

Heatmap illustrating the expression levels of various metabolites across different samples, labeled horizontally as UCON and UDS. Each cell is colored based on expression value—from blue (low) to red (high). Dendrograms on the sides indicate hierarchical clustering of samples and metabolites, emphasizing patterns in the data distribution.

Heat map of metabolite expression in urine of patients with cirrhosis.

Table 2

No. Rt(min) m/z Formula Metabolite HMDB ID VIP p value Trend
DIS vs. CON
1 4.64 316.25 C17H33NO4 decanoylcarnitine HMDB0000651 1.02 0.03 down
2 7.22 313.27 C19H38O4 MG(0:0/16:0/0:0) HMDB0011533 1.63 0.00 up
3 4.41 302.23 C16H31NO4 6-methyloctanoylcarnitine HMDB0241048 1.36 0.00 down
4 2.87 222.08 C11H13NO4 N-acetyl-L-tyrosine HMDB0000866 1.93 0.00 up
5 4.58 330.26 C18H35NO4 8-methyldecanoylcarnitine HMDB0240799 1.16 0.02 down
6 4.12 302.23 C16H31NO4 5-methyloctanoylcarnitine HMDB0241049 1.43 0.00 down
7 4.22 591.32 C33H42N4O6 urobilin HMDB0004160 2.02 0.00 up
8 4.52 655.28 C36H38N4O8 coproporphyrin III HMDB0000570 1.66 0.00 up
9 4.29 655.28 C36H38N4O8 coproporphyrin I HMDB0000643 1.67 0.00 up
10 2.59 209.09 C10H12N2O3 L-kynurenine HMDB0000684 1.17 0.00 up
11 3.80 206.08 C11H11NO3 indole-3-lactic acid HMDB0000671 1.43 0.00 up
12 1.43 204.13 C9H18NO4 o-acetylcarnitine HMDB0000201 1.05 0.03 up
13 3.30 244.15 C12H21NO4 2-tiglylcarnitine HMDB0241658 1.63 0.00 down
14 4.05 464.30 C26H43NO6 glycocholic acid HMDB0000138 1.52 0.00 up
15 1.24 130.09 C6H11NO2 pipecolic acid HMDB0000070 1.10 0.03 up
16 3.93 130.07 C4H9N3O2 creatine HMDB0000064 1.03 0.01 up
17 1.67 262.16 C12H23NO5 5-hydroxypentanoylcarnitine HMDB0241655 1.09 0.00 down
18 4.49 491.30 C28H44O7 24(28)-dehydromakisterone HMDB0302988 2.02 0.00 up
19 3.55 318.19 C15H27NO6 3-methylheptanedioylcarnitine HMDB0241046 1.28 0.00 up
20 3.66 340.18 C17H25NO6 deca-3,5,7-trienedioylcarnitine HMDB0241129 2.36 0.00 down
21 2.55 302.16 C14H23NO6 hept-4-enedioylcarnitine HMDB0241688 1.53 0.00 down
22 3.81 204.07 C11H11NO3 indolelactic acid HMDB0000671 2.02 0.00 up
23 3.79 329.21 C21H30O4 corticosterone HMDB0001547 1.75 0.00 down
24 2.53 194.05 C9H9NO4 alpha-hydroxyhippuric acid HMDB0002404 1.06 0.02 down
25 0.85 264.12 C9H17N3O6 creatine riboside HMDB0240254 1.44 0.00 up
26 1.51 314.12 C14H19NO7 tyramine glucuronide HMDB0010328 1.12 0.03 up
27 7.89 297.28 C19H36O2 cis-10-nonadecenoic acid HMDB0013622 3.32 0.00 up
28 3.35 318.19 C15H27NO6 octanedioylcarnitine HMDB0241733 1.11 0.00 up
29 5.94 437.26 C21H43O7P lysoPA(18:0/0:0) HMDB0007854 1.63 0.00 up
30 3.74 539.25 C27H41O9P PA(22:6(4Z,7Z,10Z,12E,16Z,19Z)-OH(14)/2:0) HMDB0266564 2.18 0.00 down
31 3.96 627.37 C33H57O9P PA(20:4(5Z,8Z,11Z,13E)-OH(15S)/10:0) HMDB0262660 1.38 0.00 up
32 2.80 340.10 C15H17NO8 5-hydroxy-6-methoxyindole glucuronide HMDB0010363 1.69 0.00 down
33 4.72 386.29 C21H39NO5 3-hydroxy-cis-5-tetradecenoylcarnitine HMDB0013330 1.22 0.00 up
34 3.92 398.25 C21H35NO6 tetradeca-9,11-dienedioylcarnitine HMDB0241391 1.20 0.00 up
35 4.01 329.21 C21H30O4 reichsteins Substance S HMDB0000015 1.77 0.00 down
36 5.65 435.24 C21H41O7P LysoPA(18:1(9Z)/0:0) HMDB0007855 1.27 0.00 up
37 3.71 184.10 C9H15NO3 3-hepteneoylglycine HMDB0094729 1.52 0.00 down
38 4.35 199.13 C11H20O3 2-hydroxyundec-8-enoic acid HMDB0340607 1.38 0.00 up
39 1.62 260.15 C12H21NO5 4-oxopentanoylcarnitine HMDB0241664 1.06 0.01 down
40 3.91 298.20 C16H27NO4 nona-3,6-dienoylcarnitine HMDB0241762 1.82 0.00 down
41 4.13 157.12 C9H18O2 isopropyl hexanoate HMDB0040430 1.01 0.00 down
42 3.31 305.03 C11H14O8S 4-hydroxy-5-(dihydroxyphenyl)-valeric acid-O-sulfate HMDB0059978 1.36 0.02 down
43 2.03 262.16 C12H23NO5 3-hydroxyisovalerylcarnitine HMDB0061189 1.06 0.01 down
44 3.97 593.33 C33H42N4O6 urobilinogen HMDB0004158 1.53 0.00 up
45 3.69 316.21 C16H29NO5 6-hydroxynon-4-enoylcarnitine HMDB0241753 1.06 0.01 down
46 1.56 240.11 C9H13N5O3 dihydrobiopterin HMDB0000038 1.11 0.00 down
47 3.48 300.18 C15H25NO5 3-hydroxyocta-2,5-dienoylcarnitine HMDB0241723 1.29 0.00 down
48 3.75 331.23 C21H30O3 17 alpha-Hydroxyprogesterone HMDB0000374 1.21 0.00 down
49 4.27 414.30 C26H43NO5 glycoursodeoxycholic acid HMDB0000708 1.12 0.02 up
50 4.10 204.09 C8H15NO6 N-acetylmannosamine HMDB0001129 1.16 0.02 up
51 3.48 182.08 C9H13NO3 normetanephrine HMDB0000819 1.43 0.00 down

Differential metabolites in the urine of patients with cirrhosis.

3.3 Analysis of differential metabolic pathways in serum and urine

Pathway analysis of differential metabolites was performed on the MetaboAnalyst6.0 website (see Footnote 2). The key metabolic pathways influencing liver fibrosis were screened out based on raw p < 0.05 and impact > 0.02. It was found that liver cirrhosis affected 11 metabolic pathways (Figure 6A) in serum samples, including glycerophospholipid metabolism, alanine, aspartate and glutamate metabolism, arginine biosynthesis, tryptophan metabolism, beta-alanine metabolism, pantothenate and CoA biosynthesis, arachidonic acid metabolism, taurine and hypotaurine metabolism, porphyrin metabolism, primary bile acid biosynthesis, pyrimidine metabolism. Cirrhosis affects seven metabolic pathways (Figure 6B) in urine samples, includingarginine and proline metabolism, porphyrin metabolism, amino sugar and nucleotide sugar metabolism, tryptophan metabolism, pentose and glucuronate interconversions, steroid hormone biosynthesis, glycerophospholipid metabolism. Tryptophan metabolism, glycerophospholipid metabolism, and porphyrin metabolism were co-regulated by both serum and urine.

Figure 6

Graphs A and B depict metabolic pathways with pathway impact on the x-axis and negative logarithm of the p-value on the y-axis. In Graph A, glycerophospholipid metabolism and pantothenate and CoA biosynthesis are highlighted. Graph B emphasizes tryptophan metabolism and glycerophospholipid metabolism. Color and size of circles indicate significance and impact, with red and larger circles showing more significance and impact.

Metabolic pathways affected by patients with cirrhosis: (A) serum and (B) urine.

4 Discussion

The development and progression of liver diseases are intricately associated with individuals’ dietary patterns and lifestyles. Liver cirrhosis has multifaceted etiologies, predominantly stemming from hepatitis viruses, alcohol intake, and bile stasis. Portal hypertension and impaired liver function serve as critical indicators for the diagnosis of liver cirrhosis. The liver assumes a pivotal role in metabolism, and analyzing the serum and urine of patients with liver cirrhosis aids in the identification of metabolic biomarkers related to the condition. To enhance the accuracy and efficiency of our investigation, we implemented stringent participant selection criteria. Inclusion was restricted to patients with clinically diagnosed cirrhosis, while individuals with severe comorbidities or recent use of medications known to affect metabolic profiles were excluded. This rigorous screening protocol ensured cohort homogeneity, thereby allowing clear observation of cirrhosis-associated metabolic alterations despite the limited sample size (n = 28).

Our research findings indicate that in the serum of patients with liver cirrhosis, the levels of bilirubin, glycocholate, L-glutamate, urobilinogen, taurocholate, and 7α-hydroxy-3-oxo-4-cholestenoate are significantly upregulated, whereas those of serotonin, taurine, glycerophosphocholine, L-aspartate, arachidonate, kynurenine, and dihydrouracil are significantly downregulated. Moreover, glycoursodeoxycholic acid, urobilin, glycocholic acid, and urobilinogen are the metabolites that show significant increases in both the serum and urine of patients with liver cirrhosis. The liver plays an important role in metabolic processes. Serum samples from patients with cirrhosis were analyzed to find serum biomarkers associated with cirrhosis. Bilirubin in the blood is primarily produced by reticuloendothelial cells in the spleen (17). Bilirubin binds to human peroxisome proliferator-activated receptor α (PPAR α), which contributes to the reduction of hepatic fat accumulation and the alleviation of obesity and metabolic dysfunctions (18, 19). Nevertheless, excessively elevated bilirubin levels (>150 μM) may trigger the responses of pruritus receptors, and pruritus serves as the initial manifestation of cholestasis (20). The binding of bilirubin to albumin is frequently employed to predict the long-term prognosis of patients with hepatocellular carcinoma (21) and shows a significant correlation with the histological stage of patients with primary biliary cirrhosis (22). Bilirubin exhibits reactive oxygen species scavenging activity and immunosuppressive effects (23). However, an excessively high level of bilirubin may act as an indicator of cirrhosis.

The study found that alterations in the function of the glutamate-nitric oxide-cGMP pathway in cirrhosis can cause changes in the nervous system, giving rise to hepatic encephalopathy. In this process, nitric oxide activates soluble guanylate cyclase, leading to increased expression levels of nitric oxide in the cerebral cortex and thereby affecting the patients’ neurological function (24). The dysregulation of bile acid metabolism and its subsequent accumulation in the liver result in progressive liver injury and fibrosis. Cirrhosis can cause bile duct rupture, leading to bile acid leakage. Therefore, the accumulation of bile acids in the blood is associated with cirrhosis (25). Primary bile acids are synthesized within hepatocytes via cholesterol oxidation. Subsequently, they are conjugated with glycine or taurine and then excreted into the gallbladder by the bile salt export pump (26). Studies have found that taurocholate promotes the activation of hepatic stellate cells through the S1PR2/p38 MAPK/YAP signaling pathway (27). The levels of amino acids also changed significantly in patients with cirrhosis. Taurine and L-aspartate are non-essential amino acids, and their levels decreased significantly in patients with cirrhosis. Taurine is mainly synthesized in the liver and kidneys. It can reduce lipid peroxidation products, alleviate inflammation, and prevent calcium accumulation. The deficiency of taurine in hepatocytes leads to severe liver injury and triggers compensatory hepatocyte proliferation, which is closely related to the development of cirrhosis (28). L-ornithine L-aspartate has been utilized for the prevention and treatment of hepatic encephalopathy in cirrhotic patients (29).

Our study offers additional evidence indicating that the deficiency of L-aspartate is linked to the development of cirrhosis. Research has revealed that the level of tryptophan significantly increased in cirrhotic rats (30). Kynurenine, a product of tryptophan catabolism, is associated with signaling within the host microbiome, immune cell responses, and neuronal excitability. A decrease in the total activity of tryptophan 2,3-dioxygenase in liver tissue impedes the conversion of tryptophan to kynurenine. Consequently, a reduced level of kynurenine can serve as an indicator of cirrhosis. Serotonin induces the contraction and proliferation of smooth muscle cells and stimulates endothelial cells to release vasodilator substances. The level of serotonin is implicated in diseases such as hypertension, primary pulmonary hypertension, and cirrhosis (31). Our study further validates the close association between serotonin and the development of cirrhosis. Tryptophan metabolism, glycerophospholipid metabolism, and porphyrin metabolism were identified as cirrhosis-associated pathways detected in both serum and urine. Tryptophan metabolism is closely linked to gut microbiota. As an essential amino acid acquired exclusively through dietary intake, tryptophan plays a central role in metabolism. Within the gut, tryptophan is metabolized into 5-hydroxytryptamine (5-HT, serotonin), kynurenine, and various indole derivatives, demonstrating significant associations with the pathogenesis and progression of obesity, diabetes, non-alcoholic fatty liver disease, and atherosclerosis (32).

Serum and urine samples offer advantages such as simple operation, short processing time, good repeatability, and low cost, facilitating rapid disease diagnosis. However, in clinical applications, various factors that may affect the results need to be carefully taken into account. Before sample collection, health education for patients should be enhanced, and standardizing sample collection and storage procedures is essential. Moreover, it is necessary to improve the operational proficiency and comprehensive capabilities of laboratory technicians to eliminate the influence of subjective factors on test results. It should be noted that a limitation of this study is the absence of a validation dataset. Constrained by factors such as patient availability, geographical distribution, and ethical approval, we were unable to assemble an independent validation cohort. Although alternative measures like strict sample screening, quality control, and multiple-testing correction were implemented, the lack of external validation may affect the generalizability of our findings. Future research should aim to address this limitation by including larger, multi-center validation cohorts to enhance the robustness and clinical applicability of the results. In this study, metabolomics data were not employed to construct models for predicting the likelihood or severity of fibrosis. Despite conducting a detailed characterization of the urinary and serum metabolomes, we did not conduct further in-depth analysis of these datasets to establish predictive frameworks. The development of models using serum metabolomics alone, urine metabolomics alone, or integrated serum-urine data presents a significant opportunity for clinical translation. Accurately predicting fibrosis progression is of critical importance for early diagnosis and therapeutic intervention in clinical practice. Future studies should concentrate on leveraging these metabolomic profiles to develop prediction models with greater clinical utility. For instance, machine learning algorithms could be utilized to integrate multi-metabolite features from both serum and urine, and then establish robust fibrosis classifiers that can be validated in larger clinical cohorts. Simultaneously, this study did not take into account the metabolite ratios between serum and urine as potential indicators of disease status. Metabolite ratios may offer better pathophysiological insights than individual metabolite concentrations, as changes in relative abundance across biological matrices often reflect the underlying disease mechanisms. In specific pathological conditions, key metabolite ratios may show progressive changes that correlate with disease progression. Future research should prioritize exploring the relationships between serum-urine metabolite ratios and clinical disease states, which may potentially uncover novel diagnostic or prognostic biomarkers to guide therapeutic strategies.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by Medical Ethics Committee of the First Affiliated Hospital of Baotou Medical College (Ethical approval number: 2022026). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

XL: Conceptualization, Funding acquisition, Formal analysis, Writing – review & editing, Writing – original draft, Data curation. RW: Data curation, Formal analysis, Conceptualization, Methodology, Writing – review & editing, Writing – original draft, Software. HZ: Data curation, Methodology, Conceptualization, Writing – original draft. RL: Methodology, Software, Writing – original draft. HC: Supervision, Conceptualization, Writing – review & editing, Writing – original draft, Funding acquisition, Visualization, Resources. SS: Data curation, Funding acquisition, Conceptualization, Validation, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the National Natural Science Foundation of China (81641137, 81760782; and 82060784); the 11th “Grassland Talents “Talent Project of Inner Mongolia Autonomous Region, China [(2021)08]; the “Grassland Talents “Youth Innovation and Entrepreneurship Talent Project of Inner Mongolia Autonomous Region, China (Q2017046); the Natural Science Foundation of Inner Mongolia Autonomous Region, China (2023MS08022 and 2019MS08189); the Project of Baotou Medical College Innovation Team Development, China (bycxtd-02); Inner Mongolia Natural Science Foundation Project (2023MS08029); the Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region (NJYT24079); Science and Technology Program of Baotou Health Commission (wsjkwkj005); Clinical medicine +X multidisciplinary joint Scientific research Fund project of Baotou Medical College (BYJJ-DXK 2022016).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1646323/full#supplementary-material

Abbreviations

CLD, Chronic liver diseases; HMDB, human metabolome database; KEGG, kyoto encyclopedia of genes and genomes; PCA, Principal Component Analysis; PLS-DA, Partial Least Squares-Discriminant Analysis; QC, quality control; UPLC-QTOF/MS, ultra-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry; VIP, variable importance in projection.

References

  • 1.

    Byass P . The global burden of liver disease: a challenge for methods and for public health. BMC Med. (2014) 12:159. doi: 10.1186/s12916-014-0159-5,

  • 2.

    Marcellin P Kutala BK . Liver diseases: a major, neglected global public health problem requiring urgent actions and large-scale screening. Liver Int. (2018) 38:26. doi: 10.1111/liv.13682,

  • 3.

    Cheemerla S Balakrishnan M . Global epidemiology of chronic liver disease. Clin Liver Dis. (2021) 17:36570. doi: 10.1002/cld.1061,

  • 4.

    Ginés P Quintero E Arroyo V Terés J Bruguera M Rimola A et al . Compensated cirrhosis: natural history and prognostic factors. Hepatology. (1987) 7:1228. doi: 10.1002/hep.1840070124,

  • 5.

    Smith A Baumgartner K Bositis C . Cirrhosis: diagnosis and management. Am Fam Physician. (2019) 100:75970.

  • 6.

    Cromer M Wilcox CM Shoreibah M . Beta-blockers and cirrhosis: striking the right balance. Am J Med Sci. (2024) 367:22834. doi: 10.1016/j.amjms.2024.01.009,

  • 7.

    Turco L Reiberger T Vitale G La Mura V . Carvedilol as the new non-selective beta-blocker of choice in patients with cirrhosis and portal hypertension. Liver Int: Official J Int Association for the Study of the Liver. (2023) 43:118394. doi: 10.1111/liv.15559,

  • 8.

    Chen Q Mei L Zhong R Han P Wen J Han X et al . Serum liver fibrosis markers predict hepatic decompensation in compensated cirrhosis. BMC Gastroenterol. (2023) 23:317. doi: 10.1186/s12876-023-02877-2,

  • 9.

    Raikou VD . Renoprotective strategies. World J Nephrol. (2024) 13:89637. doi: 10.5527/wjn.v13.i1.89637,

  • 10.

    Wishart DS . Current Progress in computational metabolomics. Brief Bioinform. (2007) 8:27993. doi: 10.1093/bib/bbm030,

  • 11.

    Bujak R Struck-Lewicka W Markuszewski MJ Kaliszan R . Metabolomics for laboratory diagnostics. J Pharm Biomed Anal. (2015) 113:10820. doi: 10.1016/j.jpba.2014.12.017,

  • 12.

    Chen DQ Guo Y Li P . New insights into a novel metabolic biomarker and therapeutic target for chronic kidney disease. IMNA. (2024) 11:00019:e24. doi: 10.1097/IMNA-D-24-00019

  • 13.

    Yang QJ Zhao JR Hao J Li B Huo Y Han YL et al . Serum and urine metabolomics study reveals a distinct diagnostic model for cancer cachexia. J Cachexia Sarcopenia Muscle. (2018) 9:7185. doi: 10.1002/jcsm.12246,

  • 14.

    James EL Parkinson EK . Serum metabolomics in animal models and human disease. Curr Opin Clin Nutr Metab Care. (2015) 18:47883. doi: 10.1097/MCO.0000000000000200,

  • 15.

    Zhang A Sun H Wu X Wang X . Urine metabolomics. Clin Chim Acta. (2012) 414:659. doi: 10.1016/j.cca.2012.08.016,

  • 16.

    Brix F Demetrowitsch T Jensen-Kroll J Zacharias HU Szymczak S Laudes M et al . Evaluating the effect of data merging and Postacquisition normalization on statistical analysis of untargeted high-resolution mass spectrometry based urinary metabolomics data. Anal Chem. (2024) 96:3340. doi: 10.1021/acs.analchem.3c01380,

  • 17.

    Creeden JF Gordon DM Stec DE Hinds TD Jr . Bilirubin as a metabolic hormone: the physiological relevance of low levels. Am J Physiol Endocrinol Metab. (2021) 320:E191207. doi: 10.1152/ajpendo.00405.2020,

  • 18.

    Hinds TD Jr Stec DE . Bilirubin, a cardiometabolic signaling molecule. Hypertension. (2018) 72:78895. doi: 10.1161/HYPERTENSIONAHA.118.11130,

  • 19.

    Hinds TD Stec DE . Bilirubin safeguards cardiorenal and metabolic diseases: a protective role in health. Curr Hypertens Rep. (2019) 21:87. doi: 10.1007/s11906-019-0994-z,

  • 20.

    Meixiong J Vasavda C Green D Zheng Q Qi L Kwatra SG et al . Identification of a bilirubin receptor that may mediate a component of cholestatic itch. eLife. (2019) 8:e44116. doi: 10.7554/eLife.44116,

  • 21.

    Wang J Zhang Z Yan X Li M Xia J Liu Y et al . Albumin-bilirubin (ALBI) as an accurate and simple prognostic score for chronic hepatitis B-related liver cirrhosis. Dig Liver Dis. (2019) 51:11728. doi: 10.1016/j.dld.2019.01.011,

  • 22.

    Chan AW Chan RC Wong GL Wong VW Choi PC Chan HL et al . New simple prognostic score for primary biliary cirrhosis: albumin-bilirubin score. J Gastroenterol Hepatol. (2015) 30:13916. doi: 10.1111/jgh.12938,

  • 23.

    Vitek L Hinds TD Jr Stec DE Tiribelli C . The physiology of bilirubin: health and disease equilibrium. Trends Mol Med. (2023) 29:31528. doi: 10.1016/j.molmed.2023.01.007,

  • 24.

    Rodrigo R Felipo V . Brain regional alterations in the modulation of the glutamate-nitric oxide-cGMP pathway in liver cirrhosis. Role of hyperammonemia and cell types involved. Neurochem Int. (2006) 48:4727. doi: 10.1016/j.neuint.2005.10.014,

  • 25.

    Liberal R Grant CR . Cirrhosis and autoimmune liver disease: current understanding. World J Hepatol. (2016) 8:115768. doi: 10.4254/wjh.v8.i28.1157,

  • 26.

    de Aguiar Vallim TQ Tarling EJ Edwards PA . Pleiotropic roles of bile acids in metabolism. Cell Metab. (2013) 17:65769. doi: 10.1016/j.cmet.2013.03.013,

  • 27.

    Yang J Tang X Liang Z Chen M Sun L . Taurocholic acid promotes hepatic stellate cell activation via S1PR2/p38 MAPK/YAP signaling under cholestatic conditions. Clin Mol Hepatol. (2023) 29:46581. doi: 10.3350/cmh.2022.0327,

  • 28.

    Warskulat U Borsch E Reinehr R Heller-Stilb B Mönnighoff I Buchczyk D et al . Chronic liver disease is triggered by taurine transporter knockout in the mouse. FASEB J. (2006) 20:5746. doi: 10.1096/fj.05-5016fje,

  • 29.

    Goh ET Stokes CS Sidhu SS Vilstrup H Gluud LL Morgan MY . L-ornithine L-aspartate for prevention and treatment of hepatic encephalopathy in people with cirrhosis. Cochrane Database Syst Rev. (2018) 5:CD012410. doi: 10.1002/14651858.CD012410.pub2,

  • 30.

    Ohta Y Uemura M Saito K Sasaki E Ishiguro I . Relationship between the level of serum L-tryptophan and its hepatic uptake and metabolism in rats with carbon tetrachloride-induced liver cirrhosis. Amino Acids. (1996) 10:36978. doi: 10.1007/BF00805864,

  • 31.

    Culafic DM Mirkovic DS Vukcevic MD Rudic JS . Plasma and platelet serotonin levels in patients with liver cirrhosis. World J Gastroenterol. (2007) 13:57503. doi: 10.3748/wjg.v13.i43.5750,

  • 32.

    Luo Z Liu Y Wang X Fan F Yang Z Luo D . Exploring tryptophan metabolism: the transition from disturbed balance to diagnostic and therapeutic potential in metabolic diseases. Biochem Pharmacol. (2024) 230:116554. doi: 10.1016/j.bcp.2024.116554,

Summary

Keywords

cirrhosis, metabolic markers, serum metabolomics, UPLC-QTOF/MS, urine metabolomics

Citation

Li X, Wang R, Zhou H, Li R, Chang H and Shi S (2026) Serum-urine metabolic integration via UPLC-QTOF/MS uncovers shared pathway biomarkers for cirrhosis diagnosis. Front. Med. 12:1646323. doi: 10.3389/fmed.2025.1646323

Received

13 June 2025

Revised

07 December 2025

Accepted

11 December 2025

Published

09 February 2026

Volume

12 - 2025

Edited by

Antonio Riva, Roger Williams Institute of Liver Studies (King's College London & Foundation for Liver Research), United Kingdom

Reviewed by

Jiangxin Wang, Shenzhen University, China

Chandrima Gain, University of California, Los Angeles, United States

Updates

Copyright

*Correspondence: Songli Shi, ; Hong Chang,

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics