Analysis of Serum Bile Acid Profile Characteristics and Identification of New Biomarkers in Lean Metabolic Dysfunction-Associated Fatty Liver Disease Based on LC-MS/MS

Abstract

Objectives:

Plasma bile acid (BA) has been widely studied as pathophysiological factors in chronic liver disease. But the changes of plasma BA level in lean metabolic dysfunction-associated fatty liver disease (MAFLD) remains unclear. Here, we clarified the BA metabolic characteristics of lean MAFLD and explored its significance and mechanism as a marker.

Methods:

We employed ultra-performance liquid chromatography tandem mass spectrometry based on BA metabonomics to characterize circulating bile acid in lean MAFLD patients. Explore its significance as serum biomarkers by further cluster analysis, functional enrichment analysis, and serum concentration change analysis of differential BAs. Evaluation of diagnostic value of differential BAs by ROC analysis.

Results:

A total of 65 BAs were detected and 17 BAs were identified which showed different expression in the lean-MAFLD group compared with the normal group. Functional annotation and enrichment analysis of KEGG and HMDB showed that differential BAs were mainly related to bile acid biosynthesis, bile secretion, cholesterol metabolism, and familial hypercholangitis, involving diseases including but not limited to cirrhosis, hepatocellular carcinoma, chronic active hepatitis, colorectal cancer, acute liver failure, and portal vein obstruction. ROC analysis displayed that the 6 BA metabolites (GCDCA-3S, GUDCA-3S, CDCA-3S, NCA, TCDCA, and HDCA) exhibited well differential diagnostic ability in discriminating between lean MAFLD patients and normal individuals with an area under the curve (AUC) ⩾0.85.

Conclusions:

We delineated the characteristics of BA level in patients with lean MAFLD, and identified 6 potential plasma BA biomarkers of lean MAFLD.

Plain Language Summary

Analysis of serum bile acid profile characteristics and identification of new biomarkers in fatty liver disease accompanied by metabolic abnormalities in people with normal weight based on the technology of high-resolution mass spectrometry

Objectives: The physique of lean MAFLD patient is normal or even leaner. They often does not pay enough attention to the onset of fatty liver disease. Plasma bile acids (BAs) have been extensively studied as pathophysiological actors in chronic liver disease. But the changes of plasma BA level in fatty liver disease accompanied by metabolic abnormalities in people with normal weight remains unclear. Here, we clarified the BA metabolic characteristics of lean MAFLD and explored its significance and mechanism as a marker. Methods: we employed an advanced mass spectrometry technology to characterize circulating bile acid in lean lean MAFLD patients. To explore its significance as a marker by bioinformatics methods, such as cluster analysis, functional enrichment analysis, and relative content change analysis of differential BAs. Evaluation diagnostic accuracy and determine threshold points of BAs by Receiver Operating Characteristic analysis. Results: A total of 65 BAs were detected and 17 BAs were identified which showed different expression in the lean MAFLD group compared with the normal group. Bioinformatics analysis showed that differential BAs were mainly related to bile acid biosynthesis, bile secretion, cholesterol metabolism, and familial hypercholangitis, involving diseases including but not limited to cirrhosis, hepatocellular carcinoma, chronic active hepatitis, colorectal cancer, acute liver failure, and portal vein obstruction. ROC analysis displayed that the six BA metabolites (GCDCA-3S, GUDCA-3S, CDCA-3S, NCA, TCDCA and HDCA) exhibited well differential diagnostic ability in discriminating between lean MAFLD patients and normal individuals with an area under the curve (AUC) ≥ 0.85. Conclusions: We delineated the characteristics of BA level in patients with lean MAFLD, and identified six potential plasma BA biomarkers of lean MAFLD. This strategy provided broad clinical application prospects for disease assessment.

Keywords

Metabolic dysfunction-associated fatty liver disease bile acid profile LC-MS/MS metabonomics serum biomarkers

Introduction

Metabolic dysfunction-associated fatty liver disease (MAFLD) once called non-alcoholic fatty liver disease (NAFLD). In addition to fatty denaturation and lipid accumulation in hepatocytes, NAFLD is often accompanied by cardiovascular and metabolic diseases such as hypertension, atherosclerosis, obesity, and diabetes.

Given the importance of metabolic abnormalities in this disease, NAFLD has now been renamed as MAFLD in order to reflect the disease characteristics and risks more appropriately.¹ With the prevalence of metabolic syndrome (insulin resistance, obesity, diabetes, etc.), MAFLD has become the main cause of chronic liver disease worldwide in the past decades. The progression of MAFLD can be described by the level of inflammatory activity and fibrosis stage, instead of simple fatty liver and steatohepatitis. Without therapeutic intervention, some MAFLD patients may subsequently progress to cirrhosis and eventually develop into hepatocellular carcinoma (HCC). The pathogenesis of MAFLD progression remains to be fully elucidated. MAFLD development depends on multiple hepatic insults via several different pathways.² Among them, intestine-liver axis interaction and abnormal BA metabolism cannot be ignored.

In HFC diet induced obese mice model, A muciniphila efficiently increased mitochondrial oxidation and BA metabolism in the intestine-liver axis by regulating the metabolism of L-aspartic acid, and improved steatosis and inflammation in MAFLD.³ Vitamin C and vitamin D3 may be an effective method for treating MAFLD by regulating gut microbiota and BA metabolism via the intestine-liver axis, which may be a potential drug target for future MAFLD interventions.⁴ The latest research has found that 3-succinylated cholic acid derived from Bacteroides uniformis strains has been shown to be negatively correlated with liver injury in MAFLD patients.⁵ BA are the main component of bile. It can not only promote the digestion and absorption of lipids, but also have important physiological signals and metabolic regulation effects. The quantitative detection of the BA profiles plays an important role in assessing liver disease.^6,7 There were significant differences in BA profiles among patients with chronic liver disease of different causes which suggested that BA profiles has clinical potential in distinguishing liver injury types.⁸ Compared with healthy individuals, there was no difference in total bile acid (TBA) levels among NAFLD patients, but BA composition changed significantly.⁹ Animal experiments showed significant changes in BA levels in the enterohepatic circulation of non-alcoholic steatohepatitis, and the effects could be corrected by diet.¹⁰ A certain study has shown that primary BA is associated with follow-up liver-related events among NAFLD patients, suggesting that BA metabolism may predict the prognosis of NAFLD.¹¹ Together, increasing evidence suggested a close relationship between BA plasma concentrations and MAFLD patients.

It was noteworthy that “lean” or “non-obese” MAFLD has been extensively reported among in Asian populations in recent years and there were few studies on BA metabolism in lean MAFLD since its redefinition in 2020. MAFLD patients in the advanced stage were usually accompanied by lobular inflammation and liver fibrosis. Here, we assessed the changes in serum BA spectrum inMAFLD progression patients with normal BMI with the aim of explore appropriate serum BA markers to diagnose lean MAFLD by LC-MS/MS technology and metabonomics.

Materials and Methods

Sample collection

This study was a retrospective study. From January to June 2023, serum samples of 15 diagnosed lean MAFLD patients (lean-MAFLD group) and 15 healthy controls (Normal group) were collected from Nanjing Jinling Hospital. The diagnosis of MAFLD was made according to the diagnostic criterion proposed in 2020, emphasized the role of systemic metabolic disorders in triggering liver diseases.¹ In this study, enrolled MAFLD patients with normal BMI (18.5-23.9), serum TG ⩾ 1.70 mmol/L, and hepatic steatosis through abdominal ultrasound scan, were accompanied by diabetes or arterial hypertension for at least 5 years and failure of receive standardized medical treatment within the past 6 to 12 months. The exclusion criteria were as follows: Other chronic liver disease (including, but not limited to viral hepatitis, autoimmune liver diseases and malignancy); medical history of drugs affecting liver function; hazardous alcohol intake. Patients were screened for the presence of lobular inflammation and fibrosis by liver biopsy upon additional informed consent. Histopathological diagnosis of all liver tissue samples was conducted by experienced pathologists. The Scheuer scoring system was adopted as the histological standard of Liver inflammation (G0~G4) and fibrosis (S0~S4). The grades of liver inflammation were classified into the following 5 stages: G0, no inflammation; G1, portal inflammation; G2, mild piecemeal necrosis or acidophil bodies; G3, moderate piecemeal necrosis and severe focal cell damage; and G4, widely bridging necrosis and lobular structural abnormalities. Liver fibrosis was scored as follows: S0, no fibrosis; S1, portal tracts expansion and portal fibrosis without septa; S2, normal lobular structural and portal fibrosis with rare septa; S3, disordered lobular structure and numerous septa without cirrhosis; and S4, early or confirmed cirrhosis. Hepatic steatosis was scored on a 4-point scale (F0~F3). The percentage of fat in liver parenchyma: F0, <5%; F1, 5%~33%; F2, 34%~66%; F3, >66%. Fasting serum from 5 mL peripheral venous blood was collected from each participant and stored in a −80° C refrigerator.

The research protocol was approved by the ethics committees of Nanjing University School of Medicine Affiliated Jinling Hospital. This study was conducted according to the ethical guidelines of the Declaration of Helsinki. Participants willingly agreed to participate in the study and written consents were taken.

BAs detection and quantification by LC-MS/MS

Bile acids contents were detected by MetWare (http://www.metware.cn/) based on the AB Sciex QTRAP 6500 LC-MS/MS platform.

Sample preparation and extraction Samples (50 μL) were extracted with 200 μL methanol/acetonitrile(v/v = 2:8). 10 μL 13 kinds internal standard mixed solution (1 μg/mL) was added into the extract as internal standards (IS) for the quantification. Put the samples at −20°C for 10 minutes to precipitated protein. Then centrifugation for 10 minutes (12 000 r/min, and 4°C), the supernatant was transferred to clean plastic microtubes. The extracts were evaporated to dryness, reconstituted in 100 μL 50% methanol (V/V) for further LC-MS/MS analysis.

The HPLC and ESI-MS/MS conditions based on references.^12,13

HPLC Conditions The sample extracts were analyzed using an LC-ESI-MS/MS system (UHPLC, ExionLC™ AD, https://sciex.com.cn/; MS, Applied Biosystems 6500 Triple Quadrupole, https://sciex.com.cn/). The analytical conditions were as follows, HPLC: column, Waters ACQUITY UPLC HSS T3 C18 (100 mm × 2.1 mm i.d., 1.8 µm); solvent system, water with 0.01% acetic acid and 5 mmol/L ammonium acetate (A), acetonitrile with 0.01% acetic acid (B); The gradient was optimized at 5% to 40%B in 0.5 minute, then increased to 50% B in 4 minutes, then increased to 75% B in 3 minutes, and then 75% to 95% in 2.5 minutes, washed with 95% B for 2 minutes, finaly ramped back to 5% B (12-14 minutes); flow rate, 0.35 mL/min; temperature, 40°C; injection volume: 3 μL.

ESI-MS/MS Conditions Linear ion trap (LIT) and triple quadrupole (QQQ) scans were acquired on a triple quadrupole-linear ion trap mass spectrometer (QTRAP), QTRAP^® 6500+ LC-MS/MS System, equipped with an ESI Turbo Ion-Spray interface, operating in negative ion mode and controlled by Analyst 1.6.3 software (Sciex). The ESI source operation parameters were as follows: ion source, ESI-; source temperature 550°C; ion spray voltage (IS) −4500 V; curtain gas (CUR) was set at 35 psi, respectively. Bile acids were analyzed using scheduled multiple reaction monitoring (MRM). Data acquisitions were performed using Analyst 1.6.3 software (Sciex). Multiquant 3.0.3 software (Sciex) was used to quantify all metabolites. Mass spectrometer parameters including the declustering potentials (DP) and collision energies (CE) for individual MRM transitions were done with further DP and CE optimization. A specific set of MRM transitions were monitored for each period according to the metabolites eluted within this period.

Data analysis

Unsupervised PCA (principal component analysis) was performed by statistics function prcomp within R (www.r-project.org). The data was unit variance scaled before unsupervised PCA.

Significantly regulated BAs between groups were determined by Variable Importance in Projection (VIP) and absolute LogFC (fold change). VIP values were extracted from OPLS-DA result, which also contain score plots and permutation plots, was generated using R package MetaboAnalystR. The data was log transform (log) and mean centering before OPLS-DA. A permutation test (200 permutations) was performed in order to avoid overfitting. Screening criteria are as follows:

For comparison between 2 groups:

(1) BA metabolites with VIP > 1 were selected. The VIP value represents the influence intensity of the inter group differences of corresponding BA metabolites on the classification and discrimination of each group of samples in the model. It is generally believed that BA metabolites with VIP > 1 have significant differences.

(2) BA metabolites with fold change ⩾2 and fold change ⩽ 0.5 were selected . If the difference in BA metabolites between the lean-MAFLD group and the Normal group is more than twice or less than 0.5, the difference is considered significant.

The diagnostic power of the different BAs were evaluated by the receiver operating characteristic (ROC) analysis by Sangerbox3.0 (http://www.sangerbox.com/tool ).

Public databases

KEGG database (http://www.kegg.jp/kegg/compound/) and HMDB database (http://www.hmdb.ca/) were used for functional annotation and the enrichment analysis.

Results

Statistics of clinical characteristics of the study cohort

Participants in this study included 15 lean MAFLD patients and 15 healthy people. The demographic, clinical, and laboratory characteristics of the 2 groups of participants were shown in Table 1.

Table 1.

Statistics of clinical characteristics of the study cohort.

Vairables	Normal (n = 15)	lean-MAFLD (n = 15)	P value
Sex(F/M)	8/7	9/6	.713
Age (median ± SD, range)	40.80 ± 7.69 (25-53)	44.80 ± 9.61 (24-64)	.219
BMI (kg/m², median ± SD)	21.35 ± 1.39	21.97 ± 1.27	.208
TG (mmol/L, median ± SD )	1.01 ± 0.36	2.31 ± 0.53	.000
Hypertension (Y/N)	0/15	6/9	.006
Diabetes (Y/N)	0/15	12/3	.000
Steatosis grade <5%/5%~33%/34%~66%/>66%	15/0/0/0	0/2/9/4	.000
G0/1/2/3/4	15/0/0/0/0	0/5/9/1/0	.000
S0/1/2/3/4	15/0/0/0/0	0/4/8/3/0	.000

Screening of BAs metabolites with significant differences

Using UPLC-MS/MS quantitative BA detection method, a total of 65 differential BA metabolites were detected. The standard curve linear equations and correlation coefficients of BA were shown in the Supplemental Table 1. The quantitative results of BA concentration were shown in Supplemental Table 2. OPLS-DA model showed that the BA metabolites in the lean-MAFLD group and the Normal group were clearly separated (Figure 1). The resulting R²X, R²Y, and Q² of permutation test model validation were 0.481, 0.8, and 0.236. By the cutoff of fold change ⩾2 or ⩽0.5, and VIP ⩾ 1.0, 17 differential BAs were identified in patients with lean-MAFLD compared to Normal, of which 14 were upregulated and 3 were downregulated (Table 2 and Figure 2A). The upregulated BAs were 6-ketoLCA, CDCA, GCDCA-3S, GUDCA-3S, CDCA-3S, NCA, CDCA-3Gln, GCA, GCDCA, TUDCA, Tω-MCA, TCA, THDCA, TCDCA, while the downregulated BAs were HDCA, 23-DCA, MDCA. In order to observe the variation of BAs content, we used normalization (Unit Variance Scaling, UV Scaling) to normalize the original content of differential metabolites identified by applying the screening criteria, and drew the heat map with the R software package (Figure 2B). The relative contents of differential BAs between lean-MAFLD group and Normal group were shown in Figure 3.

Figure 1.

Establishment and verification of OPLS-DA model. (A) OPLS-DA Score Chart of lean-MAFLD versus Normal. (B) OPLS-DA S-plot of lean-MAFLD versus Normal. The abscissa represents the covariance of the principal component and metabolite, and the ordinate represents the correlation coefficient of the principal component and metabolite. The closer the metabolites are to the upper right and lower left corners, the more significant the difference is. The red dots indicate that the metabolites have VIP values greater than or equal to 1, and the green dots indicate that the metabolites have VIP values less than 1. (C) The model had good prediction ability. Orange represents the random grouping model R2Y, and purple represents the random grouping Model Q2, the values represented by the black arrows are the R2X, R2Y, and Q2 values of the original model.

Table 2.

BAs with significant differences in lean-MAFLD versus normal.

Index	Compounds	Class	VIP	Log2FC	Type
6-ketoLCA	5-β-Cholanic Acid-3α-ol-6-one	BAs	1.52	1.32	up
CDCA	Chenodeoxycholic acid	BAs	1.09	1.33	up
GCDCA-3S	Glycochenodeoxycholic acid 3 sulfate disodium salt	BAs	1.92	1.85	up
GUDCA-3S	Glycoursodeoxycholic acid 3 sulfate sodium	BAs	1.52	3.00	up
CDCA-3S	chenodeoxycholic acid3-sulfate disodium salt	BAs	1.74	1.61	up
NCA	norcholic acid	BAs	2.12	1.66	up
CDCA-3Gln	Chenodeoxycholic acid-3-β-D-glucuronide	Glucuronide	1.25	1.19	up
GCA	Glycocholic acid	Glyco	1.74	1.74	up
GCDCA	Glycochenodeoxycholic acid	Glyco	1.27	1.09	up
TUDCA	Tauroursodeoxycholic acid	Tauro	1.06	1.60	up
Tω-MCA	Tauro-ω-muricholic acid sodium salt	Tauro	1.18	1.31	up
TCA	Taurocholic acid	Tauro	1.81	1.23	up
THDCA	Taurohyodeoxycholic acid (sodium salt)	Tauro	1.13	1.14	up
TCDCA	Taurochenodeoxycholic acid	Tauro	1.84	1.91	up
HDCA	Hyodeoxycholic acid	BAs	1.41	−1.08	down
23-DCA	Nor-Deoxycholic acid	BAs	1.49	−1.54	down
MDCA	murideoxycholic acid	BAs	1.10	−2.06	down

Figure 2.

Seventeen significant different BAs in lean-MAFLD versus Normal. (A) Bar chart of fold change. Red represents up-regulation, and green represents down-regulation. (B) Heat map analysis. Red blocks indicate upregulation, whereas green indicates down-regulation.

Figure 3.

Relative contents of differential BAs between lean-MAFLD and Normal. The box in the middle represents the quartile range, the thin black line extending from it represents the 95% confidence interval, the black horizontal line in the middle is the median, and the external shape represents the distribution density of the data.

Metabolic pathway analysis

KEGG pathway enrichment analysis of differential BAs found that the plasma between the lean-MAFLD group and Normal group were different in the metabolism pathways of primary bile acid biosynthesis, metabolic pathways, bile secretion, cholesterol metabolism, taurine, and hypotaurine metabolism (Figure 4A and B and Table 3). The detailed annotation information of the BAs with significant differences screened in the KEGG database was shown in Supplemental Table 3.

Figure 4.

Classification and enrichment of KEGG/HMDB pathway of differential BAs. (A) KEGG classification diagrams. The ordinate is the name of the KEGG metabolic pathway, and the abscissa is the number of BAs from the annotation to the pathway and the proportion of the number of BAs to the total number of annotated metabolites. (B) KEGG enrichment analysis graphs. (C) Enrichment of HMDB pathway of differential BAs. The color of bubbles represents the value of adjusted P value, and the size of bubbles represents the number of counts.

Table 3.

KEGG annotation results for differential Bas.

Index	Compounds	cpd_ID	kegg_map
6-ketoLCA	5-β-Cholanic Acid-3α-ol-6-one	-	-
CDCA	Chenodeoxycholic acid	C02528	ko00120, ko01100, ko04976
GCDCA-3S	Glycochenodeoxycholic acid 3 sulfate disodium salt	C15559	-
GUDCA-3S	Glycoursodeoxycholic acid 3 sulfate sodium	-	-
CDCA-3S	Chenodeoxycholic acid 3-sulfate disodium salt	-	-
NCA	Norcholic acid	-	-
CDCA-3Gln	Chenodeoxycholic acid-3-β-D-glucuronide	-	-
GCA	Glycocholic acid	C01921	ko00120, ko01100, ko04976, ko04979
GCDCA	Glycochenodeoxycholic acid	C05466	ko00120, ko04976, ko04979
TUDCA	Tauroursodeoxycholic acid	C16868	-
Tω-MCA	Tauro-ω-muricholic acid sodium salt	-	-
TCA	Taurocholic acid	C05122	ko00120, ko00430, ko01100, ko04976, ko04979
THDCA	Taurohyodeoxycholic acid (sodium salt)	-	-
TCDCA	Taurochenodeoxycholic acid	C05465	ko00120, ko04976, ko04979
HDCA	Hyodeoxycholic acid	C15517	-
23-DCA	Nor-Deoxycholic acid	-	-
MDCA	Murideoxycholic acid	C15515	-

ko00120 Primary bile acid biosynthesis.

ko01100 Metabolic pathways.

ko04976 Bile secretion.

ko04979 Cholesterol metabolism.

ko00430 Taurine and hypotaurine metabolism.

HMDB is a widely used database of more than 40 000 endogenous metabolites and more than 5000 protein (DNA) information associated with it. The HMDB sub-database SMPDB also provides human metabolism, metabolic disease pathways, metabolite signals, and drug activity pathways. Since the SMPDB only opens the image of Primary pathways download, the pathway enrichment analysis was only performed for the Primary pathway. HMDB Primary pathway enrichment analysis of differential BAs found that the plasma between the lean-MAFLD group and Normal group were different in Bile Acid Biosynthesis, Congenital Bile Acid Synthesis Defect Type II, Congenital Bile Acid Synthesis Defect Type III, Familial Hypercholanemia (FHCA), Zellweger Syndrome, Cerebrotendinous Xanthomatosis (CTX), 27-Hydroxylase Deficiency (Table 4 and Figure 4C). The detailed annotation information of the BAs with significant differences screened in SMPDB was shown in Supplemental Table 4.

Table 4.

HMDB annotation results for differential Bas.

Index	Compounds	hmdb_ID	primary_SMPDB_ID
6-ketoLCA	5-β-Cholanic acid-3α-ol-6-one	-	-
CDCA	Chenodeoxycholic acid	HMDB00518	SMP0000035; SMP0000314; SMP0000318; SMP0000317; SMP0000316; SMP0000315; SMP0000720
GCDCA-3S	Glycochenodeoxycholic acid 3 sulfate disodium salt	HMDB0002496	-
GUDCA-3S	Glycoursodeoxycholic acid 3 sulfate sodium	HMDB0002409	-
CDCA-3S	Chenodeoxycholic acid3-sulfate disodium salt	-	-
NCA	Norcholic acid	-	-
CDCA-3Gln	Chenodeoxycholic acid-3-β-D-glucuronide	HMDB02430	-
GCA	Glycocholic acid	HMDB00138	SMP0000035; SMP0000314; SMP0000318; SMP0000317; SMP0000316; SMP0000315; SMP0000720
GCDCA	Glycochenodeoxycholic acid	HMDB00637	SMP0000035; SMP0000314; SMP0000318; SMP0000317; SMP0000316; SMP0000315; SMP0000720
TUDCA	Tauroursodeoxycholic acid	HMDB00874	-
Tω-MCA	Tauro-ω-muricholic acid sodium salt	-	-
TCA	Taurocholic acid	HMDB00036	SMP0000035; SMP0000314; SMP0000318; SMP0000317; SMP0000316; SMP0000315; SMP0000720
THDCA	Taurohyodeoxycholic acid (sodium salt)	-	-
TCDCA	Taurochenodeoxycholic acid	HMDB00951	SMP0000035; SMP0000314; SMP0000318; SMP0000317; SMP0000316; SMP0000315; SMP0000720
HDCA	Hyodeoxycholic acid	HMDB0000733	-
23-DCA	Nor-Deoxycholic acid	-	-
MDCA	Murideoxycholic acid	HMDB00811	-

SMP0000035 Bile Acid Biosynthesis.

SMP0000314 Congenital Bile Acid Synthesis Defect Type II.

SMP0000318 Congenital Bile Acid Synthesis Defect Type III.

SMP0000317 Familial Hypercholanemia (FHCA).

SMP0000316 Zellweger Syndrome.

SMP0000315 Cerebrotendinous Xanthomatosis (CTX).

SMP0000720 27-Hydroxylase Deficiency.

Differential BA metabolites associated with disease

According to HMDB database, the disease information associated with differential BAs was found. The results were shown in Table 5. Related diseases included, but were not limited to, liver cirrhosis, hepatocellular carcinoma, chronic active hepatitis, colorectal cancer, acute liver failure, portal vein obstruction, etc.

Table 5.

Differential BAs associated diseases.

Index	Compound name	Hmdb diseases
HDCA	Hyodeoxycholic acid	Primary biliary cirrhosis
6-ketoLCA	5-β-Cholanic acid-3α-ol-6-one	-
23-DCA	Nor-deoxycholic acid	-
MDCA	Murideoxycholic acid	-
CDCA	Chenodeoxycholic acid	Cystic fibrosis \| Biliary atresia \| Cirrhosis \| Hepatocellular carcinoma \| Primary biliary cirrhosis
GCDCA-3S	Glycochenodeoxycholic acid 3 sulfate disodium salt	-
GUDCA-3S	Glycoursodeoxycholic acid 3 sulfate sodium	-
CDCA-3S	Chenodeoxycholic acid3-sulfate disodium salt	-
NCA	Norcholic acid	-
CDCA-3Gln	Chenodeoxycholic acid-3-β-D-glucuronide	-
GCA	Glycocholic acid	Hepatobiliary diseases \| Celiac disease \| Biliary atresia \| Choledochal cysts \| Intrahepatic biliary hypoplasia \| Neonatal hepatitis \| Alpha-1-antitrypsin deficiency \| Cystic fibrosis \| Galactosemia type 1 \| Chronic active hepatitis \| Glycogen storage disease \| Acute liver failure \| Portal vein obstruction \| Wilson’s disease \| Colorectal cancer \| Crohn’s disease \| Ulcerative colitis
GCDCA	Glycochenodeoxycholic acid	Colorectal cancer \| Biliary atresia
TUDCA	Tauroursodeoxycholic acid	Colorectal cancer
Tω-MCA	Tauro-ω-muricholic acid sodium salt	-
TCA	Taurocholic acid	Hepatocellular carcinoma \| Cirrhosis \| Colorectal cancer \| Crohn’s disease \| Ulcerative colitis \| Metastatic melanoma \| Biliary atresia
THDCA	Taurohyodeoxycholic acid (sodium salt)	-
TCDCA	Taurochenodeoxycholic acid	Biliary atresia \| Colorectal cancer \| Crohn’s disease \| Ulcerative colitis

Predictive value of bile acid spectrum subcomponents on lean MAFLD

All the identified differential BAs were used to perform an ROC analysis. Six of the BAs with an AUC > 0.85 (GCDCA-3S, GUDCA-3S, CDCA-3S, NCA, TCDCA, HDCA) displayed a high discriminating capacity (Figure 5). At the cut-off value less than 109.6, 65.8, 2.31, 2.03, 44.61 ng/mL for GCDCA-3S, GUDCA-3S, CDCA-3S, NCA, TCDCA, the sensitivity, and the specificity were 86.7%, 80%, 100%, 93.3%, 86.7% and 86.7%, 93.3%, 66.7%, 80%, 86.7%, respectively. At the cut-off value higher than 4.54 ng/mL for HDCA, the sensitivity and the specificity were 93.3% and 80%, respectively.

Figure 5.

Predictive value of bile acid spectrum subcomponents on lean-MAFLD. GCDCA-3S yielded an AUC of 0.89 (95% CI: 0.76-1.00); GUDCA-3S yielded an AUC of 0.92 (95% CI: 0.83-1.00); CDCA-3S yielded an AUC of 0.89 (95% CI: 0.78-1.00); NCA yielded an AUC of 0.94 (95% CI: 0.86-1.00); TCDCA yielded an AUC of 0.92 (95% CI: 0.81-1.00); HDCA yielded an AUC of 0.86 (95% CI: 0.70-1.00).

Discussion

The increase of MAFLD incidence rate is closely related to the global obesity crisis¹⁴; however, up to 10%-20% of MAFLD patients have a normal body mass index (BMI).^15,16 The “lean” or “non-obese” MAFLD patients usually do not have systemic obesity but are related to visceral fat accumulation. It seemed that the pathogenesis of lean MAFLD is at least partially related to the changes in gut microbiota composition.¹⁷ The intestinal microbiota could alter the metabolism of BA, cause the onset of MAFLD by regulating the activation of the farnesome X receptor (FXR), thereby affect the homeostasis of fat and glucose. Lean people with MAFLD was associated with comorbidities (such as: diabetes, hypertension, metabolic syndrome, and cardiovascular disease) and cardiovascular mortality,¹⁸ as assessed by histology¹⁹ and retrospective study.²⁰ Other studies challenged this view.^21,22 Although the lean MAFLD phenotype typically manifests as a milder form of disease, it may also exhibit a wide range of histopathological features related to metabolic dysfunction-associated steatohepatitis (MASH), including steatosis, hepatocyte ballooning degeneration, lobular inflammation, and liver fibrosis.^22,23 The physique of lean MAFLD patient is normal or even leaner. They often do not pay enough attention to the onset of fatty liver disease. So far, specific biomarkers to diagnose lean MAFLD are missing.

In this study, 17 differentially expressed BAs were identified via BA metabolomics between the lean-MAFLD group and the Normal group. KEGG analysis showed that the 17 differentially expressed BAs were mainly enriched in metabolic pathways, bile secretion, cholesterol metabolism, taurine, and hypotaurine metabolism that have been reported to be related to glycolipid metabolism disorder and oxidative stress. It provides clues for further comprehending the role of BAs in the occurrence and progression of lean MAFLD.

Our results showed that serum bile acids in lean MAFLD patients are mainly upregulated compared to normal individuals. The results of this study were similar to those of previous studies, which indicated that the concentrations of glycochenodeoxycholate, glycocholate, and taurocholate, such as GCDCA, CDCA, GCA, TCA, TCDCA, etc, were significantly higher in subjects with non-alcoholic steatohepatitis (NASH).^24
-26 In this study, ROC curve analysis shows that GCDCA-3S, GUDCA-3S, CDCA-3S, NCA, TCDCA, and HDCA have a high predictive value for hepatic injury in lean MAFLD patients. Many studies demonstrated numerous interesting connections between changes in plasma BA profile and the severity of NAFLD. In more severe forms of NAFLD, plasma levels of GUDCA, GCDCA, and TCDCA were associated with fibrosis and lobular inflammation in NAFLD and NASH.²⁷ NCA can promote the progression and immune invasion of HCC by inhibiting the FXR signaling pathway,²⁸ which may play a certain role in the progression of MAFLD into HCC. CDCA, the major hydrophobic primary bile acids, was significantly modulate the expression of the human miRNA and mRNA in primary human hepatocytes and has a related impact on the function of gene networks involved in BA lipid, and drug metabolism.²⁹ In addition, excessive CDCA may be an endogenous danger factor, which activates NLRP3 inflammasome and triggers liver inflammation during cholestasis.³⁰ We speculate that excessive CDCA may induce inflammation and fibrosis in MAFLD.

In this study, HDCA level decreased in lean MAFLD compared with healthy people. Similar findings indicate that a group of gut microbiota-modified BAs, hyodeoxycholic acid (HDCA) species, are negatively correlated with the presence and severity of NAFLD.^31,32 HDCA also influences cholesterol and glucose homeostasis in LDL receptor-null mice. The method of supplementing HDCA could inhibit intestinal cholesterol absorption in mice, reduce plasma VLDL/IDL/LDL cholesterol levels, improve high-density lipoprotein function, as well as lower obesity and blood glucose levels.³³ It means that HDCA may be a protective indicator.

It should be mentioned that this study has a number of limitations. The sample size was relatively small and needs to be increased in further research. Furthermore, the power analysis for sample size calculation was not done. Also, we diagnosed MAFLD patients based on the novel consensus of 2020 diagnostic criteria, metabolic at-risk criteria only include diabetes, hypertension and hyperlipidemia. Finally, this study did not differentiate the severity of lean MAFLD patients.

Conclusion

In summary, we have demonstrated significant changes in composition of bile acids in lean MAFLD. GCDCA-3S, GUDCA-3S, CDCA-3S, NCA, TCDCA, and HDCA were screened and identified as potential serum biomarkers for diagnosing lean MAFLD by LC-MS/MS technique. This strategy provided broad clinical application prospects for disease assessment.

Supplemental Material

sj-docx-1-end-10.1177_11795514241282253 – Supplemental material for Analysis of Serum Bile Acid Profile Characteristics and Identification of New Biomarkers in Lean Metabolic Dysfunction-Associated Fatty Liver Disease Based on LC-MS/MS

Supplemental material, sj-docx-1-end-10.1177_11795514241282253 for Analysis of Serum Bile Acid Profile Characteristics and Identification of New Biomarkers in Lean Metabolic Dysfunction-Associated Fatty Liver Disease Based on LC-MS/MS by Bing Wang, Fei Zhang, Hong Qiu, Yujie He, Haotian Shi and Yuerong Zhu in Clinical Medicine Insights: Endocrinology and Diabetes

Footnotes

Acknowledgements

The authors would like to thank to pathology department, for pathological evaluation of liver tissue.

Declarations

ORCID iD

Bing Wang

Supplemental Material

Supplemental material for this article is available online.

References

Eslam

Newsome

Sarin

, et al. A new definition for metabolic dysfunction-associated fatty liver disease: an international expert consensus statement. J Hepatol. 2020;73:202-209.

Juanola

Martínez-López

Francés

Gómez-Hurtado

Non-alcoholic fatty liver disease: metabolic, genetic, epigenetic and environmental risk factors. Int J Environ Res Public Health. 2021;18:5227.

Rao

Kuang

, et al. Gut Akkermansia muciniphila ameliorates metabolic dysfunction-associated fatty liver disease by regulating the metabolism of L-aspartate via gut-liver axis. Gut Microbes. 2021;13:1-19.

Chen

Zhao

Mei

, et al. Vitamin C and vitamin D3 alleviate metabolic-associated fatty liver disease by regulating the gut microbiota and bile acid metabolism via the gut-liver axis. Front Pharmacol. 2023;14.

Nie

Luo

Wang

, et al. Gut symbionts alleviate MASH through a secondary bile acid biosynthetic pathway. Cell. 2024;187:2717-2734.e33.

Bing

YL.

The role of bile acid metabolism in the occurrence and development of NAFLD. Front Mol Biosci. 2022;9:1089359.

Puri

Daita

Joyce

, et al. The presence and severity of nonalcoholic steatohepatitis is associated with specific changes in circulating bile acids. Hepatology. 2018;67:534-548.

Sang

Wang

Zhou

, et al. Bile acid profiles are distinct among patients with different etiologies of chronic liver disease. J Proteome Res. 2021;20:2340-2351.

Caussy

Hsu

Singh

, et al. Serum bile acid patterns are associated with the presence of NAFLD in twins, and dose-dependent changes with increase in fibrosis stage in patients with biopsy-proven NAFLD. Aliment Pharmacol Ther. 2019;49:183-193.

10.

Gillard

Clerbaux

Nachit

, et al. Bile acids contribute to the development of non-alcoholic steatohepatitis in mice. JHEP Rep. 2022;4:100387.

11.

Wegermann

Howe

Henao

, et al. Serum bile acid, vitamin E, and serotonin metabolites are associated with future liver-related events in nonalcoholic fatty liver disease. Hepatol Commun. 2021;5:608-617.

12.

Song

Sun

Xiao

Dihydromyricetin prevents obesity via regulating bile acid metabolism associated with the farnesoid X receptor in ob/ob mice. Food Funct. 2022;13:2491-2503.

13.

Zhou

, et al. Plasma oxidative lipidomics reveals signatures for sepsis-associated acute kidney injury. Clin Chim Acta. 2023;551:117616.

14.

Chen

Wang

Morgan

Robson

Annett

Obesity, non-alcoholic fatty liver disease and hepatocellular carcinoma: current status and therapeutic targets. Front Endocrinol. 2023;14:1148934.

15.

Zou

Yeo

, et al. Global prevalence, incidence, and outcomes of non-obese or lean non-alcoholic fatty liver disease: a systematic review and meta-analysis. Lancet Gastroenterol Hepatol. 2020;5:739-752.

16.

Zheng

Rios

, et al. Global epidemiology of lean non-alcoholic fatty liver disease: a systematic review and meta-analysis. J Gastroenterol Hepatol. 2020;35:2041-2050.

17.

Duarte

SMB

Stefano

Miele

, et al. Gut microbiome composition in lean patients with NASH is associated with liver damage independent of caloric intake: A prospective pilot study. Nutr Metab Cardiovasc Dis. 2018;28:369-384.

18.

Golabi

Paik

Fukui

, et al. Patients with lean nonalcoholic fatty liver disease are metabolically abnormal and have a higher risk for mortality. Clin Diabetes. 2019;37:65-72.

19.

Denkmayr

Feldman

Stechemesser

, et al. Lean patients with non-alcoholic fatty liver disease have a severe histological phenotype similar to obese patients. J Clin Med. 2018;7:562.

20.

Zou

Yeo

Nguyen

, et al. Prevalence. Characteristics and mortality outcomes of obese, nonobese and lean NAFLD in the United States, 1999-2016. J Intern Med. 2020;288:139-151.

21.

Hagström

Nasr

Ekstedt

, et al. Risk for development of severe liver disease in lean patients with nonalcoholic fatty liver disease: A long-term follow-up study. Hepatol Commun. 2018;2:48-57.

22.

Leung

Loong

Wei

, et al. Histological severity and clinical outcomes of nonalcoholic fatty liver disease in nonobese patients. Hepatology. 2017;65:54-64.

23.

Lim

Chun

Kim

, et al. Fibrotic burden in the liver differs across metabolic dysfunction-associated fatty liver disease subtypes. Gut Liver. 2023;17:610-619.

24.

Kalhan

Guo

Edmison

, et al. Plasma metabolomic profile in non-alcoholic fatty liver disease. Metabolism. 2011;60:404-413.

25.

Rivera-Andrade

Petrick

Alvarez

, et al. Circulating bile acid concentrations and non-alcoholic fatty liver disease in Guatemala. Aliment Pharmacol Ther. 2022;56:321-329.

26.

Chávez-Talavera

Haas

Grzych

Tailleux

Staels

Bile acid alterations in nonalcoholic fatty liver disease, obesity, insulin resistance and type 2 diabetes: what do the human studies tell?

Curr Opin Lipidol. 2019;30:244-254.

27.

Nimer

Choucair

Wang

, et al. Bile acids profile, histopathological indices and genetic variants for non-alcoholic fatty liver disease progression. Metabolism. 2021;116:154457.

28.

Gong

Qin

, et al. Norcholic acid promotes tumor progression and immune escape by regulating farnesoid X receptor in hepatocellular carcinoma. Front Oncol. 2021;11:711448.

29.

Krattinger

Boström

Lee

SML

, et al. Chenodeoxycholic acid significantly impacts the expression of miRNAs and genes involved in lipid, bile acid and drug metabolism in human hepatocytes. Life Sci. 2016;156:47-56.

30.

Gong

Zhou

Zhao

, et al. Chenodeoxycholic acid activates NLRP3 inflammasome and contributes to cholestatic liver fibrosis. Oncotarget. 2016;7:83951-83963.

31.

Kuang

Wang

, et al. Hyodeoxycholic acid alleviates non-alcoholic fatty liver disease through modulating the gut-liver axis. Cell Metab. 2023;35:1752-1766.e8.

32.

Zhong

Gao

, et al. Hyodeoxycholic acid ameliorates nonalcoholic fatty liver disease by inhibiting RAN-mediated PPARα nucleus-cytoplasm shuttling. Nat Commun. 2023;14:5451.

33.

Shih

Shaposhnik

Meng

, et al. Hyodeoxycholic acid improves HDL function and inhibits atherosclerotic lesion formation in LDLR-knockout mice. FASEB J. 2013;27:3805-3817.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.04 MB