1. Background
2. Objectives
3. Methods
3.1. Study Population
3.2. Anthropometric and Biochemical Parameters
3.3. NMR Acquisition and Data Processing
3.4. Statistical Analysis
3.4.1. The Anthropometric and Clinical Experiments
3.4.2. Metabolomics Experiments
3.5. Metabolite Identification
3.6. Metabolic Pathway Analysis
4. Results
4.1. Clinical Patient Characteristics
| Control (N = 36) | NAFLD (N = 37) | P Value | |
|---|---|---|---|
| Gender, M/F | 13/23 | 17/20 | 0.478 |
| Age, y | 35.61 ± 1.72 | 38.55 ± 1.46 | 0.199 |
| Height, cm | 164.7 ± 1.46 | 164.45 ± 1.9 | 0.914 |
| Weight, kg | 66.35 ± 1.76 | 79.26 ± 2.07 | 0.000 |
| BMI, kg/m2 | 24.44 ± 0.57 | 29.46 ± 0.81 | 0.000 |
| FBS, mg/dL | 91.05 ± 1.11 | 92.19 ± 2.05 | 0.628 |
| TC, mg/dL | 170.82 ± 5.56 | 194.02 ± 7.29 | 0.014 |
| TG, mg/dL | 109.02 ± 8.49 | 147.83 ± 12.22 | 0.011 |
| LDL-C, mg/dL | 100.91 ± 4.24 | 103.97 ± 4.77 | 0.633 |
| HDL-C, mg/dL | 45.8 ± 1.21 | 45.01 ± 1.68 | 0.708 |
| ALT, IU/L | 21.51 ± 1.79 | 41.72 ± 4.69 | 0.000 |
| AST, IU/L | 18.11 ± 0.98 | 28.77 ± 2.75 | 0.001 |
| ALP, U/L | 160.42 ± 7.14 | 191.33 ± 9.24 | 0.01 |
| T Bili, mg/dL | 0.79 ± 0.04 | 0.83 ± 0.06 | 0.613 |
| D Bili, mg/dL | 0.21 ± 0.01 | 0.24 ± 0.01 | 0.202 |
Abbreviations: ALP, alkaline phosphatase; ALT, alanine aminotransferase; AST, aspartate aminotransferase; D Bili, direct bilirubin; FBS, fasting blood sugar; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; T Bili, total bilirubin; TC, total cholesterol; TG, triglyceride.
aValues are presented as mean ± SE.
bComparison between two groups (control vs. NAFLD patients) using the independent t-test and Fisher's exact test (for categorical variables such as gender).
4.2. Discrimination Between NAFLD and Control Groups Using Multivariate Analysis
Multivariate statistical analysis from NMR-based metabolic profiling. A, PCA score plot after excluding outliers; B, OPLS-DA score plot of the two study groups (yellow circle, control; purple triangle, NAFLD patients); C, OPLS-DA loading plot demonstrating discriminant variables between the two study groups (the variables on the right are elevated chemical shifts and the variables on the left are reduced chemical shifts in NAFLD patients compared with control group, related metabolites are tabulated in Table 2). D, ROC curve analysis was performed to assess the diagnostic performance of the model, AUC value was 0.99.
4.3. Metabolite Identification and Pathway Analysis of Altered Profiles
| Metabolite | δ 1H (p.p.m.)a | Fold Change | Direction of Variationb | Metabolic Pathway |
|---|---|---|---|---|
| 1- Butyric acid | 0.89, 2.13 | 1.9 | ↓ | Gut microbiome-derived metabolism |
| 2- Propionic acid | 1.05, 2.15 | 2.29 | ↓ | Gut microbiome-derived metabolism |
| 3- Isovaleric acid | 2.05 | 2.05 | ↓ | Gut microbiome-derived metabolism |
| 4- Glutamine | 2.45 | 2.36 | ↓ | Amino acid metabolism |
| 5- Glycine | 3.55 | 1.67 | ↓ | Amino acid metabolism |
| 6- Serine | 3.83 | 1.5 | ↓ | Amino acid metabolism |
| 7- Glycochenodeoxycholic acid | 0.83, 2.01, 2.19 | 1.83 | ↑ | Bile acid biosynthesis |
| 8- Taurocholic acid | 0.81, 0.85, 1.95, 3.63 | 2.17 | ↑ | Bile acid biosynthesis |
| 9- Glycocholic acid | 1.25, 1.83 | 4.15 | ↑ | Bile acid biosynthesis |
| 10- Deoxycholic acid | 1.13, 2.21 | 4.46 | ↑ | Bile acid biosynthesis |
| 11- Valine | 1.01 | 1.68 | ↑ | BCAA metabolism |
| 12- Isoleucine | 1.27, 1.45, 1.97 | 2.29 | ↑ | BCAA metabolism |
| 13- Succinic acid | 2.41 | 1.82 | ↑ | TCA cycle |
| 14- Isocitric acid | 2.97 | 1.68 | ↑ | TCA cycle |
| 15- 2-Ketoglutaric acid | 2.99 | 1.67 | ↑ | TCA cycle |
| 16- Trimethylamine | 2.89 | 2.01 | ↑ | Gut microbiome-derived metabolism |
| 17- Proline | 1.99 | 1.8 | ↑ | Amino acid metabolism |
| 18- Hydroxyproline | 3.37 | 1.82 | ↑ | Amino acid metabolism |
| 19- Tyrosine | 3.17 | 2.77 | ↑ | Amino acid metabolism |
Abbreviations: BCAA, branched-chain amino acid; TCA: tricarboxylic acid cycle.
aChemical shift scale of NMR signal used for quantification of metabolites
bIncreased or decreased metabolites in NAFLD group compared with control group.
| Pathway Name | Match Statusa | P Value | FDRb |
|---|---|---|---|
| 1- Aminoacyl-tRNA biosynthesis | 7/75 | 8.0379E-7 | 6.4303E-5 |
| 2- Primary bile acid biosynthesis | 4/47 | 3.9933E-4 | 0.011655 |
| 3- Citrate cycle (TCA cycle) | 3/20 | 4.3707E-4 | 0.011655 |
| 4- Alanine, aspartate and glutamate metabolism | 3/24 | 7.6061E-4 | 0.015212 |
| 5- Methane metabolism | 3/34 | 0.002139 | 0.030465 |
| 6- Propanoate metabolism | 3/35 | 0.0023279 | 0.030465 |
| 7- D-glutamine and D-glutamate metabolism | 2/11 | 0.003113 | 0.030465 |
| 8- Nitrogen metabolism | 3/39 | 0.0031861 | 0.030465 |
| 9- Butanoate metabolism | 3/40 | 0.0034273 | 0.030465 |
| 10- Glyoxylate and dicarboxylate metabolism | 3/50 | 0.0064674 | 0.048247 |
| 11- Cyanoamino acid metabolism | 2/16 | 0.0066339 | 0.048247 |
aThe total number of metabolites in each pathway/the number of identified metabolites in each pathway.
bThe false discovery rate.
Metabolic pathway analysis. MetaboAnalyst 3 displaying changed metabolic pathways in serum from NAFLD group compared with control group. The homosapiens library and the default algorithms were used for enrichment and topology analyses. Statistics for pathways with major change based on the p-value and FDR (P value < 0.05 and FDR < 0.05). These pathways include: a, aminoacyl-tRNA biosynthesis; b, primary bile acid biosynthesis; c, citrate cycle TCA cycle); d, alanine; aspartate and glutamate metabolism; e, methane metabolism; f, propanoate metabolism; g, D-Glutamine and D-glutamate metabolism; h, nitrogen metabolism; i, butanoate metabolism; j, glyoxylate and dicarboxylate metabolism; k, cyanoamino acid metabolism.

