Metabolomic Biomarkers in the Diagnosis of Non-Alcoholic Fatty Liver Disease

authors:

avatar Saeed Chashmniam 1 , avatar Mohsen Ghafourpour 2 , avatar Azam Rezaei Farimani 3 , 4 , avatar Ali Gholami 3 , 5 , avatar Bibi Fatemeh Nobakht Motlagh Ghoochani 3 , 4 , *

Department of Chemistry, Sharif University of Technology, Tehran, Iran
22 Bahman Hospital, Neyshabur University of Medical Sciences, Neyshabur, Iran
Noncommunicable Diseases Research Center, Neyshabur University of Medical Sciences, Neyshabur, Iran
Department of Basic Medical Sciences, Neyshabur University of Medical Sciences, Neyshabur, Iran
Department of Epidemiology and Biostatistics, School of Public Health, Neyshabur University of Medical Sciences, Neyshabur, Iran
Hepatitis Monthly: Vol.19, issue 9; e92244
published online: September 15, 2019
article type: Research Article
received: April 14, 2019
revised: July 16, 2019
accepted: August 31, 2019
DOI: https://doi.org/10.5812/hepatmon.92244

how to cite: Chashmniam S , Ghafourpour M , Rezaei Farimani A , Gholami A, Nobakht Motlagh Ghoochani B F. Metabolomic Biomarkers in the Diagnosis of Non-Alcoholic Fatty Liver Disease. Hepat Mon. 2019;19(9):e92244. https://doi.org/10.5812/hepatmon.92244.

Abstract

Background:

Nonalcoholic fatty liver disease (NAFLD) is the most abundant chronic liver disorder, because racial and ethnic differences may influence prevalence and severity of NAFLD.

Objectives:

This metabolomic study was conducted to identify the metabolic biomarkers and determine the mechanism of progress of NAFLD in Iranian patients.

Methods:

Serum samples were collected from 75 participants (37 healthy controls and 38 patients with NAFLD) after an overnight fast. The metabolome of all samples were determined by nuclear magnetic resonance (NMR) and were compared by multivariate statistical analysis.

Results:

Totally, 19 metabolomic biomarkers were identified by NMR. Compared to healthy controls, NAFLD patients had increased serum concentrations of glycochenodeoxycholic acid, taurocholic acid, glycocholic acid, deoxycholic acid, valine, isoleucine, succinic acid, isocitric acid, 2-ketoglutaric acid, trimethylamine, proline, hydroxyproline and tyrosine, while the concentrations of butyric acid, propionic acid, isovaleric acid, glutamine, glycine, and serine decreased.

Conclusions:

A robust set of biomarkers for diagnosis of NAFLD was established. Serum metabolomics biomarkers revealed changes in some amino acids and their derivatives, bile acids, short chain fatty acids, and tricarboxylic acid cycle intermediates in subjects with NAFLD compared to healthy controls. These markers could be used as indicators regarding the efficacy of therapeutic interventions.

1. Background

Nonalcoholic fatty liver disease (NAFLD) is the most abundant chronic liver disorder (1), affecting 6.3% - 33% of the population worldwide (2). The prevalence of NAFLD is nearly 12–24% among Asian countries (3). Variety of NAFLD incidence may be attributed to the differences in age, gender, ethnicity, location, genetic effects, and lifestyle modification (4, 5). Studies showed that NAFLD may progress to both cirrhosis and hepatocellular carcinoma (HCC) (6, 7). Therefore, early diagnosis of NAFLD may be important to prevent these disorders. Diagnostic approaches for detection of NAFLD include serum liver enzyme measurement, imaging analysis such as sonography and magnetic resonance imagining (MRI) and liver biopsy (8). The last one is the gold standard for diagnosis of NAFLD, but this method is invasive and expensive and lots of subjects with suspected NAFLD are reluctant to use it as a diagnostic approach. Other traditional non-invasive methods mentioned above are not precise and sometimes cannot detect the patients with NAFLD (9). Thus, these challenges motivated researchers to find a non-invasive method with high accuracy for early detection of NAFLD.

In the recent decades, “omics” technologies such as metabolomics, the newest member in the “omics” approaches have provided an opportunity to find diagnostic biomarkers, clarify the underling mechanisms, and evaluate the treatment of various diseases using the nuclear magnetic resonance (NMR) and gas/liquid chromatography mass spectrometry (GC/LC-MS) (10-12). The “omics” methods offer an advantage over traditional approaches that is detection of several molecules in a single assay versus investigation of one molecule in each assay.

2. Objectives

The present study was conducted to evaluate whether NMR-based metabolomics may detect patients with NAFLD compared to healthy controls in Iranian population.

3. Methods

3.1. Study Population

Metabolomic analyses were performed on 37 NAFLD patients and on 36 healthy controls. Subjects with NAFLD were enrolled in the study (between January and August 2017) from 22 bahman hospital (Neyshabur, Iran). The NAFLD patients had been newly diagnosed based on hepatic ultrasound. Among 37 NAFLD patients, 8 patients had hyperlipidemia and 3 patients had high blood pressure; therefore, these subjects received medications associated with their disorders. The control group included the volunteers who underwent medical checkups and were verified to have normal liver function. Healthy controls were matched for age and gender with the patient group. To verify the absence of steatosis, all healthy controls also performed hepatic ultrasound.

Blood sample was taken under fasting conditions form each participant in the morning. Serum samples were separated and were stored in aliquots at -80°C until use. The inclusion criteria were: (i) age between 22 - 63 years old; (ii) non-smoker and non-alcoholic people (intake of alcohol for women and men was less than 20 and 30 mg/d, respectively) and exclusion criteria were: (i) having acute or chronic liver disorders including viral and autoimmune hepatitis, alcoholic liver disease and other metabolic liver disorders such as hemochromatosis and Wilson’s disease; (ii) having other medical disorders such as diabetes mellitus, hyper/hypothyroidism, cancer, cardiovascular diseases (iii) being pregnant and lactating women.

The present study protocol was approved by the Ethics Committee of Neyshabur University of Medical Sciences (IR.NUMS.REC.1395.33). A consent form was obtained from all the participants before enrollment in the study.

3.2. Anthropometric and Biochemical Parameters

Body height and weight were measured with an accuracy 0.1 cm and 0.1 kg by a digital stadiometer (Model BSM 370, Seoul, Korea) and InBody (Model 770 device, Seoul, Korea), respectively. Body weight index was calculated by weight/(height)2. All blood specimens were taken in the morning after an overnight fast, then serum samples were separated and biochemical parameters were immediately measured, the rest of serum samples were stored at -80°C until used for metabolomics analyses.

3.3. NMR Acquisition and Data Processing

Preparation of serum for NMR and measurement were performed as explained previously (13, 14). Serum samples were thawed on ice and mixed with 10 % D2O (deuterium oxide, 99.9% D, Aldrich Chemicals Company) for the NMR lock signal, then transferred to 5 mm NMR tubes. NMR analysis was carried out at department of Chemistry at Sharif University of Technology, Tehran, Iran using a Bruker DRX500 MHz spectrometer operating at 500.13 MHz, equipped with 5mm high-quality NMR tubes (Sigma Aldrich., RSA). The 1H NMR analysis were performed at a constant temperature of 298 K and in blinded to the two groups. For all serum samples the Carr-Purcell-Meiboom-Gill (CPMG) spin-echo pulse sequence, π/2 - tD - π - tD experiments were used in order to detect small molecular weight metabolites (15). The 1D spectra acquisition parameters were a 10.5 µs 90° pulse, 128 transients collected into 32 k data points with a spectral width of 8389.26 Hz, an acquisition time of 1.95 seconds and a relaxation delay of 2 seconds. Prior to Fourier transformation, an exponential function corresponding to a line-broadening factor of 0.3 Hz was applied to the free induction decay (FID).

All spectra was referred to solvent and manually phased and baseline corrected by using the XWINNMR software (version 3.5, Bruker Spectrospin Ltd). All spectra were imported to MATLAB (version 6.5.1, The MathWorks, Cambridge, UK) and binned into 408 bins (an equal 0.02 ppm bin width) over a chemical shift in the range between 0.2 and 10 ppm using ProMetab software (version prometab-v3-3) (16) in MATLAB. Serum spectra regions containing water (4.2 - 5.5 ppm) were discarded. H NMR spectra were normalized and aligned by ProMetab in MATLAB. The final data were imported to SIMCA (SIMCA 14.0, Umetrics, Umeå, Sweden) for multivariate statistical analysis.

3.4. Statistical Analysis

3.4.1. The Anthropometric and Clinical Experiments

Comparisons of clinical and anthropometric data were preformed using SPSS 16.0 (SPSS, Inc., Chicago, IL) between control and patient groups. Normality distribution for each variable was evaluated by the Kolmogorov-Smirnov test. Categorical variables between the two groups were compared using Fisher’s exact test. The independent t-test was used for continues variables with normally distributed data.

3.4.2. Metabolomics Experiments

Multivariate statistical analyses were carried out by SIMCA (SIMCA 14.0, Umetrics, Umeå, Sweden). All variables (chemical shifts along with confounding factors) were imported into the SIMCA. First, unsupervised principal component analysis (PCA) was used to identify clusters and outliers in the dataset. Then, supervised orthogonal projections to latent structures-discriminant analysis (OPLS-DA) was applied to recognize differentiation chemical shifts between healthy controls and NAFLD patients. The model quality was described by three parameters (R2X, R2Y and Q2). Moreover, the model prediction performance was evaluated by receiver-operating characteristic (ROC) curve that generated from seven-fold cross validation obtained by SIMCA.

3.5. Metabolite Identification

The metabolites were recognized based on literatures and on-line databases such as human metabolome database (HMDB) (http://hmdb.ca/) and biological magnetic resonance databank (BMRB) (http://bmrb.wisc.edu) according to their chemical shifts obtained by OPLS-DA model.

3.6. Metabolic Pathway Analysis

Metabolite pathway analysis were performed using MetaboAnalyst 3.0. The discriminant metabolites between healthy controls and NAFLD patients were imported into the MetaboAnalyst software (17) in order to seek significant metabolic pathways involved in this disease. The P value and false discovery rate (FDR) < 0.05 were considered statistically significant.

4. Results

4.1. Clinical Patient Characteristics

In this study, we analyzed the metabolic profiles of 75 serum samples collected from 37 healthy controls and 38 patients with NAFLD for NMR analysis (two samples were identified as outliers and deleted from analysis) that their anthropometric and biochemical parameters were summarized in Table 1. Of the 37 participants, 19 were with grade 1 and 18 with grade 2. In the NAFLD and control groups, the gender ratio showed no significant difference. Age and height were similar in the NAFLD and control groups, whereas body weight and BMI were significantly higher in the NAFLD than in the control group. Also, NAFLD patients had significantly higher levels of serum total cholesterol, triglycerides and hepatic enzymes including alkaline phosphatase, aspartate aminotransferase, and alanine aminotransferase compared to the healthy controls (Table 1).

Table 1.

Anthropometric and Clinical Characteristics of Participantsa,b

Control (N = 36)NAFLD (N = 37)P Value
Gender, M/F13/2317/200.478
Age, y35.61 ± 1.7238.55 ± 1.460.199
Height, cm164.7 ± 1.46164.45 ± 1.90.914
Weight, kg66.35 ± 1.7679.26 ± 2.070.000
BMI, kg/m224.44 ± 0.5729.46 ± 0.810.000
FBS, mg/dL91.05 ± 1.1192.19 ± 2.050.628
TC, mg/dL170.82 ± 5.56194.02 ± 7.290.014
TG, mg/dL109.02 ± 8.49147.83 ± 12.220.011
LDL-C, mg/dL100.91 ± 4.24103.97 ± 4.770.633
HDL-C, mg/dL45.8 ± 1.2145.01 ± 1.680.708
ALT, IU/L21.51 ± 1.7941.72 ± 4.690.000
AST, IU/L18.11 ± 0.9828.77 ± 2.750.001
ALP, U/L160.42 ± 7.14191.33 ± 9.240.01
T Bili, mg/dL0.79 ± 0.040.83 ± 0.060.613
D Bili, mg/dL0.21 ± 0.010.24 ± 0.010.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

The obtained matrix containing the numbers of samples and 409 variables (408 chemical shifts along with BMI as confounding factor) was imported into SIMCA and analyzed using unsupervised PCA for detecting outliers. The PCA scoter plot showed two serum samples were found to be outside the 99% Hotelling’s T2 confidence limit and were removed for further analysis (Appendix A and B). After removing two outliers, the final PCA was performed (Figure 1A: R2X = 0.857, Q2 = 0.729). Then, the OPLS-DA was performed to identify potential metabolic biomarkers between the two study groups (Figure 1B: R2X = 0.67, R2Y = 0.982, Q2 = 0.924). The loading plot of OPLS-DA model has been illustrated in Figure 1C. This loading plot demonstrated the chemical shifts of the NMR spectra that are responsible for clustering in OPLS-DA score plot. To validate the OPLS-DA model ROC curve was used, the area under the ROC curve (AUC) value was 0.99 for the OPLS-DA model (Figure 1D).

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.
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.
Figure 1.

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

Our findings showed nineteen different metabolites that were recognized based on the literatures and on-line databases. Chemical shifts and signal multiplicity were used for identification of metabolites. Compared to healthy controls, NAFLD patients had increased serum concentrations of glycochenodeoxycholic acid, taurocholic acid, glycocholic acid, deoxycholic acid, valine, isoleucine, succinic acid, isocitric acid, 2-ketoglutaric acid, trimethylamine, proline, hydroxyproline and tyrosine, while the concentrations of butyric acid, propionic acid, isovaleric acid, glutamine, glycine, and serine decreased (Table 2).

Table 2.

List of Small Molecular Weight Metabolites That Play an Important Role in Discrimination of NAFLD Patients and Healthy Controls Using NMR

Metaboliteδ 1H (p.p.m.)aFold ChangeDirection of VariationbMetabolic Pathway
1- Butyric acid0.89, 2.131.9↓Gut microbiome-derived metabolism
2- Propionic acid1.05, 2.152.29↓Gut microbiome-derived metabolism
3- Isovaleric acid2.052.05↓Gut microbiome-derived metabolism
4- Glutamine2.452.36↓Amino acid metabolism
5- Glycine3.551.67↓Amino acid metabolism
6- Serine3.831.5↓Amino acid metabolism
7- Glycochenodeoxycholic acid0.83, 2.01, 2.191.83↑Bile acid biosynthesis
8- Taurocholic acid0.81, 0.85, 1.95, 3.632.17↑Bile acid biosynthesis
9- Glycocholic acid1.25, 1.834.15↑Bile acid biosynthesis
10- Deoxycholic acid1.13, 2.214.46↑Bile acid biosynthesis
11- Valine1.011.68↑BCAA metabolism
12- Isoleucine1.27, 1.45, 1.972.29↑BCAA metabolism
13- Succinic acid2.411.82↑TCA cycle
14- Isocitric acid2.971.68↑TCA cycle
15- 2-Ketoglutaric acid2.991.67↑TCA cycle
16- Trimethylamine2.892.01↑Gut microbiome-derived metabolism
17- Proline1.991.8↑Amino acid metabolism
18- Hydroxyproline3.371.82↑Amino acid metabolism
19- Tyrosine3.172.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 analysis was performed on metabolites detected via NMR analysis using the MetaboAnalyst 3.0 servers. This software showed eleven metabolic pathways changed in serum of patients with NAFLD. These pathways included aminoacyl-tRNA biosynthesis, primary bile acid biosynthesis, citrate cycle (TCA cycle), alanine, aspartate and glutamate metabolism, methane metabolism, propanoate metabolism, D-glutamine and D-glutamate metabolism, nitrogen metabolism, butanoate metabolism, glyoxylate and dicarboxylate metabolism, cyanoamino acid metabolism (Figure 2 and Table 3).

Table 3.

Statistics for Pathways with Major Change Based on the P Value and False Discovery Rate

Pathway NameMatch StatusaP ValueFDRb
1- Aminoacyl-tRNA biosynthesis7/758.0379E-76.4303E-5
2- Primary bile acid biosynthesis4/473.9933E-40.011655
3- Citrate cycle (TCA cycle)3/204.3707E-40.011655
4- Alanine, aspartate and glutamate metabolism3/247.6061E-40.015212
5- Methane metabolism3/340.0021390.030465
6- Propanoate metabolism3/350.00232790.030465
7- D-glutamine and D-glutamate metabolism2/110.0031130.030465
8- Nitrogen metabolism3/390.00318610.030465
9- Butanoate metabolism3/400.00342730.030465
10- Glyoxylate and dicarboxylate metabolism3/500.00646740.048247
11- Cyanoamino acid metabolism2/160.00663390.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.
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.
Figure 2.

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.

5. Discussion

NAFLD is associated with the obesity and alterations in life style affecting many body systems (17). Nevertheless, the pathogenesis is multifactorial and not clearly realized. A few studies have been conducted to find the metabolic biomarkers of this disease. However, because racial and ethnic differences may influence NAFLD prevalence and severity (18, 19), this metabolomic study was performed to detect the metabolic biomarkers and determine the mechanism of progress of NAFLD in Iranian patients.

The present study demonstrated that, 19 metabolites were altered in the patients with NAFLD compared to the control group as discussed below.

5.1. Amino Acid Metabolism

In the present study, a change was observed in the levels of amino acids including, glutamine, aspartate, proline, tyrosine, serine, glycine, isoleucine, leucine, and valine in the patients with NAFLD compared to healthy controls. These findings are not astonishing because most proteins and amino acids are metabolized in the liver. Glutamine, serine and glycine were observed to decrease in the patients with NAFLD compared to healthy controls. These amino acids are as precursors of the most important antioxidants, glutathione (GSH), directly or indirectly (20-22). Our results are in agreement with prior studies that reported reduced levels of these amino acids in plasma of patients with NAFLD (23, 24). However, one study showed that plasma serine did not change significantly in NAFLD patients in comparison with control group (25). Glutamine has been shown to decreases the oxidative stress (OS) through increasing the GSH biosynthesis and ameliorate the hepatic steatosis in the liver of NAFLD rats (26). Besides, Mardinoglu et al. showed that, serine supplementation decreases the hepatic steatosis in humans (27), and another group indicated that, supplementation with glycine might protect against NASH in rats (28).

Our findings indicated high levels of proline and hydroxyproline in the serum of the patients with NAFLD compared to healthy controls. These metabolites are the main components of the collagen protein and are potential biomarkers of elevated hepatic collagen (29). In line with our results, increased serum level of proline has been shown in subjects with NAFLD (30). Also previous studies showed that patients with alcoholic and non-alcoholic liver diseases had high mRNA expression of procollagen compared to the controls (31). Another study showed that, the synthesis of hepatic collagen was elevated in these two groups of patients (32).

Increased levels of tyrosine have been often reported in the serum of patients with NAFLD (23, 33). This is consistent with previous studies, where the aromatic amino acids were associated with the prevalence of liver diseases (34). This finding is not surprising because the liver metabolizes aromatic amino acids, and lipid accumulation in it causes impairment in hepatic metabolism of these amino acids. Down-regulation of aromatic amino acid transporter (SLC16A10) is another possible explanation, as it uptakes these amino acids from blood to liver, in NASH patients (35). Inconsistent with our findings, Mukherjee et al. (25) reported that the levels of plasma tyrosine significantly decreased in patients with NAFLD compared the to control group, but this group did not indicate any explanation for their results.

Our results are in line with previous investigations, which often found that the serum levels of branched chain amino acids (BCAA: isoleucine, leucine, and valine) increased in the subjects with obesity and/or NAFLD (36-38). The previous studies demonstrated that the activity of mitochondrial BCAA-catabolizing enzymes, branched-chain aminotransferase and branched-chain α-ketoacid dehydrogenase (BDKDH), decreased in obese population (39). Because BMI was significantly higher in the patients with NAFLD than the control group; therefore, this cofounding factor was adjusted. Similarly, after adjusting for BMI, an increase in the level of BCAA was shown.

5.2. TCA Cycle

The elevation of circulating BCAA and/or their catabolic intermediate products are involved in the up-regulation of TCA cycle in the liver of these patients (38, 40, 41). As observed, TCA cycle intermediate levels (succinate, α-ketoglutarate and isocitric acid) elevated in the patients with NAFLD compared to the healthy controls in the present study. Moreover, reduced activities of mitochondrial respiratory chain (MRC) complexes demonstrated in this malady can cause OS and mitochondrial dysfunction in the patients with NAFLD (42, 43).

5.3. Gut Microbiome-Derived Metabolites

Dysbiosis, as disruption of the normal gut microflora, has been observed in the NAFLD and metabolic disorders (44). Interestingly, in this study, some metabolites including, short chain fatty acids (SCFAs) (decreased levels of propionic acid, butyric acid, and isovaleric acid), methylamines (increased levels of trimethylamine) and primary and secondary bile acids produced by the fermentation of gut microbiota (44) were observed to alter in the serum of patients with NAFLD compared to the control group.

SCFAs are organic acids with less than 7 carbon atoms and play a main role in attenuating the obesity and dyslipidemia (45, 46). Ruminococcaceae is one of the gut bacteria physiological function of which is producing the SCFAs (47). Previous reports demonstrated a lower proportion of Ruminococcaceae in the patients with NASH compared to healthy subjects (48). Christensenellaceae, another family associated to low BMI in humans were positively correlated with total SCFAs (49). Decreased levels of SCFAs might be because of significantly higher BMI in the patients with NAFLD than the control group. However, after adjusting for BMI, a decrease in the level of SCFAs was also shown.

In the present study, a significant increase was found in the level of primary (glycocholic acid, taurocholic acid, glycochenodeoxycholate) and secondary bile acids (deoxycholic acid (DCA)) in the serum of the patients with NAFLD compared to the control group. Our results are in line with the previous investigations, which showed a higher level of bile acids (BAs) in subjects with hyperlipidemia and NAFLD (24, 50, 51), resulting from (a) increases in synthesis of BAs; (b) defective BA transports. Previous studies indicated that the mRNA expression of farnesoid X receptor (FXR) was down-regulated in patients with NAFLD compared to the controls (52). This protein suppresses cholesterol 7 alpha-hydroxylase (CYP7A1) expression. An increase in the hepatic expression of CYP7A1, the rate-limiting enzyme synthesizing BAs from cholesterol, results in the elevated synthesis of primary bile acids.

As bile acids from hepatic-portal circulation are removed by the liver; therefore, triglycerides accumulation in the liver may cause liver dysfunction and increased levels of bile acids in the circulation in NAFLD patients. Given that NAFLD patients in this study are obese, thus, it can be said that, the obesity causes defective hepatic BA transport and BMI is negatively related to the expression of liver BA transporters (53). Dysbiosis may be another reason for decreased level of secondary bile acids in these patients, as an imbalance in gut flora. Primary BAs are converted into secondary ones through deconjugation and dehydroxylation by intestinal flora. Alterations in gut microbiota have been observed in both NAFLD and obesity (54, 55). For instance, Clostridium cluster XI and XIVa producing DCA, a hydrophobic secondary BA which generates ROS and causes DNA damage, from primary bile acids have been reported to elevate in the obese mice (55).

The present study had some limitations. The diet of two study groups was not evaluated, and also the metabolites were identified by one method (NMR), thus it is recommended to verify our results other approaches.

5.4. Conclusions

The present study was the first investigation assessed the global metabolomics in the Iranian NAFLD patients using the NMR. A total of 19 significantly altered metabolites were detected in the patients compared to the controls. It was found that, some amino acids and their derivatives, bile acids, SCFAs, and TCA cycle intermediates were significantly changed. In brief, our findings demonstrated that, a combination of NMR-based metabolomics along with chemometric techniques offers an effective method for detection of markers of NAFLD with a high accuracy. Our results showed that, precursors of glutathione, glutamine, serine and glycine reduced in NAFLD patients. Therefore, supplementation with these amino acids might be useful for treatment of NAFLD patients. Furthermore, decreased levels of SCFAs could be because of altered levels of gut microbiota, as a result; supplementation with probiotics is another recommendation for treatment of these patients.

Acknowledgements

References

  • 1.

    Younossi Z, Henry L. Contribution of alcoholic and nonalcoholic fatty liver disease to the burden of liver-related morbidity and mortality. Gastroenterology. 2016;150(8):1778-85. [PubMed ID: 26980624]. https://doi.org/10.1053/j.gastro.2016.03.005.

  • 2.

    Chalasani N, Younossi Z, Lavine JE, Diehl AM, Brunt EM, Cusi K, et al. The diagnosis and management of non-alcoholic fatty liver disease: Practice Guideline by the American Association for the Study of Liver Diseases, American College of Gastroenterology, and the American Gastroenterological Association. Hepatology. 2012;55(6):2005-23. [PubMed ID: 22488764]. https://doi.org/10.1002/hep.25762.

  • 3.

    Fan JG, Saibara T, Chitturi S, Kim BI, Sung JJ, Chutaputti A, et al. What are the risk factors and settings for non-alcoholic fatty liver disease in Asia-Pacific? J Gastroenterol Hepatol. 2007;22(6):794-800. [PubMed ID: 17498218]. https://doi.org/10.1111/j.1440-1746.2007.04952.x.

  • 4.

    Schwimmer JB, McGreal N, Deutsch R, Finegold MJ, Lavine JE. Influence of gender, race, and ethnicity on suspected fatty liver in obese adolescents. Pediatrics. 2005;115(5):e561-5. [PubMed ID: 15867021]. https://doi.org/10.1542/peds.2004-1832.

  • 5.

    Mansour-Ghanaei F, Joukar F, Mobaraki SN, Mavaddati S, Hassanipour S, Sepehrimanesh M. Prevalence of non-alcoholic fatty liver disease in patients with diabetes mellitus, hyperlipidemia, obesity and polycystic ovary syndrome: A cross-sectional study in north of Iran. Diabetes Metab Syndr. 2019;13(2):1591-6. [PubMed ID: 31336526]. https://doi.org/10.1016/j.dsx.2019.03.009.

  • 6.

    Matteoni CA, Younossi ZM, Gramlich T, Boparai N, Liu YC, McCullough AJ. Nonalcoholic fatty liver disease: A spectrum of clinical and pathological severity. Gastroenterology. 1999;116(6):1413-9. [PubMed ID: 10348825]. https://doi.org/10.1016/s0016-5085(99)70506-8.

  • 7.

    Marrero JA, Fontana RJ, Su GL, Conjeevaram HS, Emick DM, Lok AS. NAFLD may be a common underlying liver disease in patients with hepatocellular carcinoma in the United States. Hepatology. 2002;36(6):1349-54. [PubMed ID: 12447858]. https://doi.org/10.1053/jhep.2002.36939.

  • 8.

    Mansour-Ghanaei R, Mansour-Ghanaei F, Naghipour M, Joukar F. Biochemical markers and lipid profile in nonalcoholic fatty liver disease patients in the PERSIAN Guilan cohort study (PGCS), Iran. J Family Med Prim Care. 2019;8(3):923-8. [PubMed ID: 31041226]. [PubMed Central ID: PMC6482810]. https://doi.org/10.4103/jfmpc.jfmpc_243_18.

  • 9.

    Hadizadeh F, Faghihimani E, Adibi P. Nonalcoholic fatty liver disease: Diagnostic biomarkers. World J Gastrointest Pathophysiol. 2017;8(2):11-26. [PubMed ID: 28573064]. [PubMed Central ID: PMC5437499]. https://doi.org/10.4291/wjgp.v8.i2.11.

  • 10.

    Zhang A, Sun H, Wang X. Saliva metabolomics opens door to biomarker discovery, disease diagnosis, and treatment. Appl Biochem Biotechnol. 2012;168(6):1718-27. [PubMed ID: 22971835]. https://doi.org/10.1007/s12010-012-9891-5.

  • 11.

    Nobakht BF, Aliannejad R, Rezaei-Tavirani M, Arefi Oskouie A, Naseri MT, Parastar H, et al. NMR- and GC/MS-based metabolomics of sulfur mustard exposed individuals: A pilot study. Biomarkers. 2016;21(6):479-89. [PubMed ID: 26984270]. https://doi.org/10.3109/1354750X.2016.1153725.

  • 12.

    Nobakht M Gh BF, Aliannejad R, Rezaei-Tavirani M, Taheri S, Oskouie AA. The metabolomics of airway diseases, including COPD, asthma and cystic fibrosis. Biomarkers. 2015;20(1):5-16. [PubMed ID: 25403491]. https://doi.org/10.3109/1354750X.2014.983167.

  • 13.

    Nobakht BF M Gh, Arefi Oskouie A, Rezaei-Tavirani M, Aliannejad R, Taheri S, Fathi F, et al. NMR spectroscopy-based metabolomic study of serum in sulfur mustard exposed patients with lung disease. Biomarkers. 2017;22(5):413-9. [PubMed ID: 27319271]. https://doi.org/10.1080/1354750X.2016.1203995.

  • 14.

    Ghoochani B, Aliannejad R, Oskouie AA, Rezaei-Tavirani M, Kalantari S, Naseri MT, et al. Metabolomics diagnostic approach to mustard airway diseases: a preliminary study. Iran J Basic Med Sci. 2018;21(1):59-69. [PubMed ID: 29372038]. [PubMed Central ID: PMC5776438]. https://doi.org/10.22038/IJBMS.2017.23792.5982.

  • 15.

    Brown FF, Campbell ID, Kuchel PW, Rabenstein DC. Human erythrocyte metabolism studies by 1H spin echo NMR. FEBS Lett. 1977;82(1):12-6. [PubMed ID: 21099]. https://doi.org/10.1016/0014-5793(77)80875-2.

  • 16.

    Viant MR. Improved methods for the acquisition and interpretation of NMR metabolomic data. Biochem Biophys Res Commun. 2003;310(3):943-8. [PubMed ID: 14550295]. https://doi.org/10.1016/j.bbrc.2003.09.092.

  • 17.

    Xia J, Broadhurst DI, Wilson M, Wishart DS. Translational biomarker discovery in clinical metabolomics: An introductory tutorial. Metabolomics. 2013;9(2):280-99. [PubMed ID: 23543913]. [PubMed Central ID: PMC3608878]. https://doi.org/10.1007/s11306-012-0482-9.

  • 18.

    Bambha K, Belt P, Abraham M, Wilson LA, Pabst M, Ferrell L, et al. Ethnicity and nonalcoholic fatty liver disease. Hepatology. 2012;55(3):769-80. [PubMed ID: 21987488]. [PubMed Central ID: PMC3278533]. https://doi.org/10.1002/hep.24726.

  • 19.

    Rich NE, Oji S, Mufti AR, Browning JD, Parikh ND, Odewole M, et al. Racial and ethnic disparities in nonalcoholic fatty liver disease prevalence, severity, and outcomes in the united states: A systematic review and meta-analysis. Clin Gastroenterol Hepatol. 2018;16(2):198-210 e2. [PubMed ID: 28970148]. [PubMed Central ID: PMC5794571]. https://doi.org/10.1016/j.cgh.2017.09.041.

  • 20.

    Zhou X, He L, Wu C, Zhang Y, Wu X, Yin Y. Serine alleviates oxidative stress via supporting glutathione synthesis and methionine cycle in mice. Mol Nutr Food Res. 2017;61(11). [PubMed ID: 28759161]. https://doi.org/10.1002/mnfr.201700262.

  • 21.

    Yu JC, Jiang ZM, Li DM. Glutamine: A precursor of glutathione and its effect on liver. World J Gastroenterol. 1999;5(2):143-6. [PubMed ID: 11819414]. [PubMed Central ID: PMC4688527]. https://doi.org/10.3748/wjg.v5.i2.143.

  • 22.

    McCarty MF, O'Keefe JH, DiNicolantonio JJ. Dietary glycine is rate-limiting for glutathione synthesis and may have broad potential for health protection. Ochsner J. 2018;18(1):81-7. [PubMed ID: 29559876]. [PubMed Central ID: PMC5855430].

  • 23.

    Gaggini M, Carli F, Rosso C, Buzzigoli E, Marietti M, Della Latta V, et al. Altered amino acid concentrations in NAFLD: Impact of obesity and insulin resistance. Hepatology. 2018;67(1):145-58. [PubMed ID: 28802074]. https://doi.org/10.1002/hep.29465.

  • 24.

    Kalhan SC, Guo L, Edmison J, Dasarathy S, McCullough AJ, Hanson RW, et al. Plasma metabolomic profile in nonalcoholic fatty liver disease. Metabolism. 2011;60(3):404-13. [PubMed ID: 20423748]. [PubMed Central ID: PMC2950914]. https://doi.org/10.1016/j.metabol.2010.03.006.

  • 25.

    Mukherjee S, Vaidyanathan K, Vasudevan DM, Das SK. Role of plasma amino acids and gaba in alcoholic and non-alcoholic fatty liver disease-a pilot study. Indian J Clin Biochem. 2010;25(1):37-42. [PubMed ID: 23105881]. [PubMed Central ID: PMC3453006]. https://doi.org/10.1007/s12291-010-0007-0.

  • 26.

    Lin Z, Cai F, Lin N, Ye J, Zheng Q, Ding G. Effects of glutamine on oxidative stress and nuclear factor-kappaB expression in the livers of rats with nonalcoholic fatty liver disease. Exp Ther Med. 2014;7(2):365-70. [PubMed ID: 24396406]. [PubMed Central ID: PMC3881322]. https://doi.org/10.3892/etm.2013.1434.

  • 27.

    Mardinoglu A, Bjornson E, Zhang C, Klevstig M, Soderlund S, Stahlman M, et al. Personal model-assisted identification of NAD(+) and glutathione metabolism as intervention target in NAFLD. Mol Syst Biol. 2017;13(3):916. [PubMed ID: 28254760]. [PubMed Central ID: PMC5371732]. https://doi.org/10.15252/msb.20167422.

  • 28.

    Zhou X, Han D, Xu R, Wu H, Qu C, Wang F, et al. Glycine protects against high sucrose and high fat-induced non-alcoholic steatohepatitis in rats. Oncotarget. 2016;7(49):80223-37. [PubMed ID: 27784003]. [PubMed Central ID: PMC5348315]. https://doi.org/10.18632/oncotarget.12831.

  • 29.

    Lai YS, Chen WC, Kuo TC, Ho CT, Kuo CH, Tseng YJ, et al. Mass-spectrometry-based serum metabolomics of a C57BL/6J mouse model of high-fat-diet-induced non-alcoholic fatty liver disease development. J Agric Food Chem. 2015;63(35):7873-84. [PubMed ID: 26262841]. https://doi.org/10.1021/acs.jafc.5b02830.

  • 30.

    Chen Y, Li C, Liu L, Guo F, Li S, Huang L, et al. Serum metabonomics of NAFLD plus T2DM based on liquid chromatography-mass spectrometry. Clin Biochem. 2016;49(13-14):962-6. [PubMed ID: 27211699]. https://doi.org/10.1016/j.clinbiochem.2016.05.016.

  • 31.

    Annoni G, Weiner FR, Colombo M, Czaja MJ, Zern MA. Albumin and collagen gene regulation in alcohol- and virus-induced human liver disease. Gastroenterology. 1990;98(1):197-202. [PubMed ID: 2293578]. https://doi.org/10.1016/0016-5085(90)91310-3.

  • 32.

    Murawaki Y, Kato S, Hirayama C. Hepatic collagen synthesis in patients with alcoholic and nonalcoholic liver disease. Gastroenterol Jpn. 1991;26(4):465-71. [PubMed ID: 1916155]. https://doi.org/10.1007/bf02782815.

  • 33.

    Jin R, Banton S, Tran VT, Konomi JV, Li S, Jones DP, et al. Amino acid metabolism is altered in adolescents with nonalcoholic fatty liver disease-an untargeted, high resolution metabolomics study. J Pediatr. 2016;172:14-19 e5. [PubMed ID: 26858195]. [PubMed Central ID: PMC5321134]. https://doi.org/10.1016/j.jpeds.2016.01.026.

  • 34.

    Morgan MY, Marshall AW, Milsom JP, Sherlock S. Plasma amino-acid patterns in liver disease. Gut. 1982;23(5):362-70. [PubMed ID: 7076013]. [PubMed Central ID: PMC1419690]. https://doi.org/10.1136/gut.23.5.362.

  • 35.

    Mariotta L, Ramadan T, Singer D, Guetg A, Herzog B, Stoeger C, et al. T-type amino acid transporter TAT1 (Slc16a10) is essential for extracellular aromatic amino acid homeostasis control. J Physiol. 2012;590(24):6413-24. [PubMed ID: 23045339]. [PubMed Central ID: PMC3533202]. https://doi.org/10.1113/jphysiol.2012.239574.

  • 36.

    Takashina C, Tsujino I, Watanabe T, Sakaue S, Ikeda D, Yamada A, et al. Associations among the plasma amino acid profile, obesity, and glucose metabolism in Japanese adults with normal glucose tolerance. Nutr Metab (Lond). 2016;13:5. [PubMed ID: 26788116]. [PubMed Central ID: PMC4717594]. https://doi.org/10.1186/s12986-015-0059-5.

  • 37.

    Newgard CB, An J, Bain JR, Muehlbauer MJ, Stevens RD, Lien LF, et al. A branched-chain amino acid-related metabolic signature that differentiates obese and lean humans and contributes to insulin resistance. Cell Metab. 2009;9(4):311-26. [PubMed ID: 19356713]. [PubMed Central ID: PMC3640280]. https://doi.org/10.1016/j.cmet.2009.02.002.

  • 38.

    Sunny NE, Kalavalapalli S, Bril F, Garrett TJ, Nautiyal M, Mathew JT, et al. Cross-talk between branched-chain amino acids and hepatic mitochondria is compromised in nonalcoholic fatty liver disease. Am J Physiol Endocrinol Metab. 2015;309(4):E311-9. [PubMed ID: 26058864]. [PubMed Central ID: PMC4537921]. https://doi.org/10.1152/ajpendo.00161.2015.

  • 39.

    She P, Van Horn C, Reid T, Hutson SM, Cooney RN, Lynch CJ. Obesity-related elevations in plasma leucine are associated with alterations in enzymes involved in branched-chain amino acid metabolism. Am J Physiol Endocrinol Metab. 2007;293(6):E1552-63. [PubMed ID: 17925455]. [PubMed Central ID: PMC2767201]. https://doi.org/10.1152/ajpendo.00134.2007.

  • 40.

    Sunny NE, Parks EJ, Browning JD, Burgess SC. Excessive hepatic mitochondrial TCA cycle and gluconeogenesis in humans with nonalcoholic fatty liver disease. Cell Metab. 2011;14(6):804-10. [PubMed ID: 22152305]. [PubMed Central ID: PMC3658280]. https://doi.org/10.1016/j.cmet.2011.11.004.

  • 41.

    Satapati S, Kucejova B, Duarte JA, Fletcher JA, Reynolds L, Sunny NE, et al. Mitochondrial metabolism mediates oxidative stress and inflammation in fatty liver. J Clin Invest. 2015;125(12):4447-62. [PubMed ID: 26571396]. [PubMed Central ID: PMC4665800]. https://doi.org/10.1172/JCI82204.

  • 42.

    Perez-Carreras M, Del Hoyo P, Martin MA, Rubio JC, Martin A, Castellano G, et al. Defective hepatic mitochondrial respiratory chain in patients with nonalcoholic steatohepatitis. Hepatology. 2003;38(4):999-1007. [PubMed ID: 14512887]. https://doi.org/10.1053/jhep.2003.50398.

  • 43.

    Begriche K, Massart J, Robin MA, Bonnet F, Fromenty B. Mitochondrial adaptations and dysfunctions in nonalcoholic fatty liver disease. Hepatology. 2013;58(4):1497-507. [PubMed ID: 23299992]. https://doi.org/10.1002/hep.26226.

  • 44.

    Leung C, Rivera L, Furness JB, Angus PW. The role of the gut microbiota in NAFLD. Nat Rev Gastroenterol Hepatol. 2016;13(7):412-25. [PubMed ID: 27273168]. https://doi.org/10.1038/nrgastro.2016.85.

  • 45.

    Zhou D, Pan Q, Xin FZ, Zhang RN, He CX, Chen GY, et al. Sodium butyrate attenuates high-fat diet-induced steatohepatitis in mice by improving gut microbiota and gastrointestinal barrier. World J Gastroenterol. 2017;23(1):60-75. [PubMed ID: 28104981]. [PubMed Central ID: PMC5221287]. https://doi.org/10.3748/wjg.v23.i1.60.

  • 46.

    Jin CJ, Sellmann C, Engstler AJ, Ziegenhardt D, Bergheim I. Supplementation of sodium butyrate protects mice from the development of non-alcoholic steatohepatitis (NASH). Br J Nutr. 2015;114(11):1745-55. [PubMed ID: 26450277]. https://doi.org/10.1017/S0007114515003621.

  • 47.

    Berni Canani R, Sangwan N, Stefka AT, Nocerino R, Paparo L, Aitoro R, et al. Lactobacillus rhamnosus GG-supplemented formula expands butyrate-producing bacterial strains in food allergic infants. ISME J. 2016;10(3):742-50. [PubMed ID: 26394008]. [PubMed Central ID: PMC4817673]. https://doi.org/10.1038/ismej.2015.151.

  • 48.

    Raman M, Ahmed I, Gillevet PM, Probert CS, Ratcliffe NM, Smith S, et al. Fecal microbiome and volatile organic compound metabolome in obese humans with nonalcoholic fatty liver disease. Clin Gastroenterol Hepatol. 2013;11(7):868-75 e1-3. [PubMed ID: 23454028]. https://doi.org/10.1016/j.cgh.2013.02.015.

  • 49.

    Dong Y, Fei P, Han Y, Guo L. Characterization of fecal microbiota, short-chain fatty acids and lactic acid concentrations in 5 - 8-year-old children with cow milk protein allergy. Iran J Pediatr. 2018;28(4). e64638. https://doi.org/10.5812/ijp.64638.

  • 50.

    Pennington CR, Ross PE, Bateson MC, Bouchier IA. Serum bile acids in patients with hyperlipidaemia. J Clin Pathol. 1978;31(1):58-62. [PubMed ID: 627619]. [PubMed Central ID: PMC476721]. https://doi.org/10.1136/jcp.31.1.58.

  • 51.

    Ferslew BC, Xie G, Johnston CK, Su M, Stewart PW, Jia W, et al. Altered bile acid metabolome in patients with nonalcoholic steatohepatitis. Dig Dis Sci. 2015;60(11):3318-28. [PubMed ID: 26138654]. [PubMed Central ID: PMC4864493]. https://doi.org/10.1007/s10620-015-3776-8.

  • 52.

    Yang ZX, Shen W, Sun H. Effects of nuclear receptor FXR on the regulation of liver lipid metabolism in patients with non-alcoholic fatty liver disease. Hepatol Int. 2010;4(4):741-8. [PubMed ID: 21286345]. [PubMed Central ID: PMC2994619]. https://doi.org/10.1007/s12072-010-9202-6.

  • 53.

    Haeusler RA, Camastra S, Nannipieri M, Astiarraga B, Castro-Perez J, Xie D, et al. Increased bile acid synthesis and impaired bile acid transport in human obesity. J Clin Endocrinol Metab. 2016;101(5):1935-44. [PubMed ID: 26684275]. [PubMed Central ID: PMC4870845]. https://doi.org/10.1210/jc.2015-2583.

  • 54.

    Henao-Mejia J, Elinav E, Jin C, Hao L, Mehal WZ, Strowig T, et al. Inflammasome-mediated dysbiosis regulates progression of NAFLD and obesity. Nature. 2012;482(7384):179-85. [PubMed ID: 22297845]. [PubMed Central ID: PMC3276682]. https://doi.org/10.1038/nature10809.

  • 55.

    Yoshimoto S, Loo TM, Atarashi K, Kanda H, Sato S, Oyadomari S, et al. Obesity-induced gut microbial metabolite promotes liver cancer through senescence secretome. Nature. 2013;499(7456):97-101. [PubMed ID: 23803760]. https://doi.org/10.1038/nature12347.