4.1. Comparison of clinical parameters between two groups
The age and gender were matched between two groups. Serum creatinine, urea nitrogen and uric acid were increased in IgAN group(P < 0.05). The level of eGFR showed no significance between normal control and IgAN groups. Total Cholesterol, triglyceride and LDL levels as well as leukocytes and neutrophils were significantly higher in the IgAN group compared with normal controls (P < 0.05). The level of serum albumin was significantly decreased in IgAN group(P < 0.05). Lymphocyte, hemoglobin, platelet and, glucose showed no significance between control and IgAN groups (Table 1).
Table 1
Comparison of clinical parameters between control and IgAN groups
Clinical parameters
|
Normal control
|
IgAN
|
P value
|
WBC(109/L)
|
5.69 ± 1.35
|
6.83 ± 1.78
|
0.012
|
Neutrophils(109/L)
|
3.34 ± 1.24
|
4.26 ± 1.36
|
0.014
|
Lymphocytes(109/L)
|
1.84 ± 0.44
|
1.93 ± 0.71
|
0.549
|
Hemoglobin(g/L)
|
140.34 ± 16.64
|
133.93 ± 17.83
|
0.137
|
Platelets(109/L)
|
230.00 (192.00, 283.50)
|
230.00 (201.50, 259.00)
|
0.938
|
Creatinine (µmol/L)
|
61.43 ± 13.10
|
80.29 ± 28.55
|
0.003
|
Urea nitrogen (mol/L)
|
4.63 ± 1.04
|
5.81 ± 1.82
|
0.003
|
Uric acid (µmol/L)
|
297.57 ± 64.69
|
370.72 ± 86.03
|
0.001
|
Albumin(g/L)
|
45.35 ± 2.52
|
40.91 ± 4.59
|
< 0.001
|
Total cholesterol(mmol/L)
|
4.35 ± 0.58
|
4.93 ± 1.05
|
0.009
|
Triglycerides (mmol/L)
|
1.01 (0.77, 1.33)
|
1.31 (1.01, 2.18)
|
0.007
|
Low-density lipoprotein (mmol/L)
|
2.51 ± 0.43
|
2.90 ± 0.77
|
0.025
|
Glucose(mmol/L)
|
4.71 ± 0.43
|
4.77 ± 0.62
|
0.580
|
eGFR(mL/min·1.73 m2)
|
103.89 (99.40,114.34)
|
101.95 (72.70, 123.44)
|
0.380
|
24-h urine protein (mg/24h)
|
--
|
774.75 (366.18, 2030.60)
|
--
|
4.2. Fecal short-chain fatty acid levels between two groups
The SCFAs such as acetic acid, propionic acid, butyric acid, isobutyric acid, valeric acid, isovaleric acid, and caproic acid were quantified by GC/MS (Table 2). The levels of acetic acid, propionic acid, butyric acid, isobutyric acid and caproic acid in the IgAN patients were significantly reduced compared with control group(P < 0.05).
Table 2
Comparison of SCFAs between control and IgAN groups
SCFAs
|
Control group
|
IgAN group
|
P value
|
Acetic acid (mmol/g)
|
79.56(66.31,135.83)
|
36.38(30.10, 65.92)
|
< 0.001
|
Propionic acid(mmol/g)
|
6.98(3.37,8.13)
|
4.16(3.05,5.47)
|
0.013
|
Butyric acid(mmol/g)
|
1.45(1.12,2.14)
|
0.66(0.49,0.96)
|
< 0.001
|
Isobutyric acid(mmol/g)
|
7.17 ± 5.26
|
0.68 ± 0.16
|
< 0.001
|
Valeric acid(mmol/g)
|
0.28(0.16,0.36)
|
0.19(0.09,0.37)
|
0.129
|
Isovaleric acid(mmol/g)
|
0.53(0.39,0.76)
|
0.39(0.28,0.67)
|
0.095
|
Caproic acid(mmol/g)
|
0.95 ± 0.45
|
0.54 ± 0.12
|
< 0.001
|
4.3. Correlation analysis of fecal SCFAs and clinical parameters
As shown in Table 3, butyric acid(r=-0.336, P = 0.010) and isobutyric acid(r=-0.298, P = 0.022) were negatively correlated with urea acid; butyric acid(r=-0.316, P = 0.016) was negatively correlated with urea nitrogen; caproic acid(r=-0.415, P = 0.025) showed negative correlation with 24-h urine protein level. Meanwhile, a strong inverse relationship between SCFA and serum lipids. Propionic acid, butyric acid and isobutyric acid showed negative correlation with total cholesterol (P < 0.05). Butyric acid, isobutyric acid and valeric acid showed negative correlation with triglyceridec. Butyric acid and isobutyric acid was negatively correlated with low-density lipoprotein(P < 0.05). Interestingly, serum albumin level showed a positive correlation with acetic acid, isobutyric acid and caproic acid(P < 0.05).
Table 3
Correlations between fecal SCFAs and clinical parameters
Variable 1
|
Variable 2
|
r
|
P value
|
Acetic acid
|
Albumin
|
0.337
|
0.010
|
Propionic acid
|
Total cholesterol
|
-0.264
|
0.045
|
Butyric acid
|
Urea nitrogen
|
-0.316
|
0.016
|
Low-density lipoprotein
|
-0.390
|
0.002
|
Total cholesterol
|
-0.386
|
0.003
|
Uric acid
|
-0.336
|
0.010
|
Triglyceride
|
-0.390
|
0.002
|
Isobutyric acid
|
Albumin
|
0.410
|
0.001
|
Urea nitrogen
|
-0.299
|
0.022
|
Total cholesterol
|
-0.306
|
0.020
|
Triglyceride
|
-0.320
|
0.014
|
Low-density lipoprotein
|
-0.279
|
0.034
|
Uric acid
|
-0.298
|
0.022
|
Valeric acid
|
Triglyceride
|
-0.277
|
0.036
|
Caproic acid
|
Albumin
|
0.430
|
0.001
|
24-h urine protein
|
-0.415
|
0.025
|
4.4. Characteristics of gut microbiota in IgAN patients
The α-diversity, including Chao1 (P = 0.0945, Fig. 1A), observed species (P = 0.051, Fig. 1B) and Simpson (P = 0.067, Fig. 1C), were approaching to the significant difference between two groups. Shannon index was different between two groups (P = 0.033, Figs. 1D). The β-diversity was significantly different between the IgAN and normal control groups by Non-Metric Multi-Dimensional Scaling(NMDS) (Stress = 0.132)(Fig. 1E) .
Overall, 18 bacterial phyla were recovered across the samples. The main phyla were Firmicutes (53.40%), Bacteroidetes (36.12%), Proteobacteria (6.89%), Actinobacteria (3.28%) in IgAN group, the same as those of Firmicutes (51.73%), Bacteroidetes (43.55%), Proteobacteria (3.08%), Actinobacteria (0.67%) in control group. The relative abundance of Actinobacteria in IgAN was higher than that in controls with significantly statistical difference(P = 0.013).
The LEfSe analysis (Supplementary Figure S1) was used to identify the metagenomic biomarker by way of class comparison. The results of LEfSe showed that 179 bacteria taxa had biologically consistent and statistically significant differences between two groups. The 91 bacteria taxa were more abundant in normal control group, and 88 bacteria taxa were more abundant in IgA nephropathy group.
4.5. Correlations between SCFAs and microbial indexes
The correlations of SCFAs and microbial indexes detected in the IgAN patients were analyzed (Fig. 2A). Acetic acid was positively associated with c_Clostridia(r = 0.357, P = 0.008), o_Clostridiales(r = 0.357, P = 0.008) and g_Eubacterium_coprostanoligenes_group(r = 0.283, P = 0.036). Butyric acid was positively associated with g_Alistipes (r = 0.278, P = 0.040). Isobutyric acid was positively associated with g_Lachnospiraceae_NK4A136_group and f_Ruminococcaceae(P < 0.05). Valeric acid was positively associated with g_Lactobacillus(r = 0.300, P = 0.026). c_Alphaproteobacteria, o__Rhizobiales, f_Rhizobiaceae, f__Enterococcaceae, g__Intestinibacter, g__Enterococcus, g_Megamonas and g__Ruminococcaceae_UCG-002 were negatively correlated with SCFAs (P < 0.05)(Fig. 2A).
The LDA (Figure F) was used to identify the metagenomic biomarker by way of class comparison. Compared with control group, the relative abundance of c_Clostridia, o_Clostridiales, g_Lachnospiraceae_NK4A136_group, g_Ruminococcaceae_UCG_010, g_Alistipes and g_Eubacterium_coprostanoligenes_group were significantly decreased in the IgAN group. Meanwhile, the relative abundance of c_Alphaproteobacteria, o__Rhizobiales, f_Rhizobiaceae, f__Enterococcaceae, g__Intestinibacter, g__Enterococcus, g_Megamonas and g__Ruminococcaceae_UCG-002 were increased significantly in the IgAN group(Fig. 2B) .
Discussion
This study is the first to examine fecal SCFAs in patients with IgAN, and we found significant differences in SCFAs between the IgAN group and normal controls, with significantly lower levels of acetic acid, propionic acid, butyric acid, isobutyric acid and caproic acid in IgAN. There were few studies consistent with our study, where propionic acid and butyric acid levels were significantly reduced in mice with diabetic nephropathy[18]. Wang et al. found significantly lower fecal SCFAs in CKD 5 than in patients with CKD stages 1–4, and negatively correlated with the renal function[19]. Wu et al. found significantly lower serum SCFAs levels in patients with CKD stages 4–5 compared to the normal group [20].
In our study there was the dysbiosis of gut microbiota in IgAN patients compared to controls, exemplified by the results of α-diversity, β-diversity and LEfSe analysis. De Angelis et al[21] reported that the lower microbial diversity was found in IgAN patients(P group and NP group) estimated by rarefaction, Chao 1and Shannon diversity index. In agreement with our study, Hu X et al[22] also demonstrated that the community richness of fecal microbiota in the IgAN patients was significantly lower than that in the healthy controls.
The reason for the significant decrease of fecal SCFA in IgAN patients, was related to the decreased relative abundances of SCFA-producing bacteria compared to the normal group. The bacteria taxa of c_Clostridia, o_Clostridiales, and g_Eubacterium_coprostanoligenes_ group were positively correlated with acetic acid, and g_Alistipes was positively correlated with butyric acid. Isobutyric acid was positively associated with g_ Lachnospiraceae_ NK4A136_group and f_Ruminococcaceae. Furthermore, The relative abundance of those were significantly decreased in IgAN group compared to control group. Especially, Clostridia, Alistipes and Ruminococcaceae were comfirmed as important SCFA-producing bacteria [23, 24]. Clostridium was reported to produce acetate by the Wood-Ljungdahl pathway[25]. Lachnospiraceae was demonstrated to use lactate and acetate to produce butyrate[26]. Baxter NT et al[27]found that gut microbiota with an increase in Ruminococcus or Clostridium were more likely to yield higher butyrate concentrations.
The relative abundance of c_Alphaproteobacteria, o__Rhizobiales, f_Rhizobiaceae, f__Enterococcaceae, g__Intestinibacter, g__Enterococcus, g_Megamonas and g__Ruminococcaceae_UCG-002 were negatively correlated with SCFAs. SCFAs often exhibited broad-spectrum antimicrobial activity, mainly due to its ability to penetrate and destroy microorganisms. It has been demonstrated that SCFA has both bactericidal and antibacterial functions against oral microorganisms. In vitro experiments have demonstrated that butyric acid can inhibit Salmonella infection [28].SCFAs are relatively inert to the bacteria that produce them, but can effectively inhibit the growth of other bacteria. [29]In addition to its direct effect on the intestinal barrier, SCFA can lower the pH in the intestinal lumen, which can directly promote the growth of commensal bacteria and inhibit the proliferation of pathogenic bacteria [30].
In this study, the negative correlations between SCFAs and kidney injury-related indicators were found. Butyric acid was negatively correlated with urea acid and urea nitrogen; caproic acid showed negative correlation with 24-h urine protein level. Wang et al. found that the serum acetate and butyrate levels in CKD 5 were significantly lower than those in CKD 1–4 patients [19]. Wu et al. found that compared with the normal group, the serum propionic acid level of CKD4-5 patients was significantly lower, and it can be used to distinguish patients with severe renal impairment from the normal group [20].
The role of short-chain fatty acids in IgAN remained unclear. Previous studies have confirmed that SCFAs modulated inflammation both in intestinal and extra-intestinal environments and possessed multifarious effects against inflammatory bowel disease and allergic airway disease by decreasing inflammatory response. [31, 32]Meanwhile, a strong relationship between IgAN and intestinal inflammation was reported, such as inflammatory bowel diseases [33] and coeliac disease [34].Qin DAI et al, found that P-aIgA1 (aggregated IgA1 from IgAN patients) promoted the proliferation of human renal mesangial cells (HMCs), and markedly increased the protein levels of HDACⅠin the cells. P-aIgA1 markedly increased the protein levels of Col1a1 and PAI-1, as well as pSmad2/3 and pStat3 in the cells. Pretreatment with SCFAs inhibitor of HDAC inpartially reversed P-aIgA1-induced cell proliferation and extracellular matrix synthesis in HMCs[35].
Meanwhile, a strong inverse relationship between SCFAs and blood lipids. Propionic acid, butyric acid and isobutyric acid showed negative correlations with total cholesterol (P < 0.05). Butyric acid, isobutyric acid and valeric acid showed negative correlation with triglyceridec. Butyric acid and isobutyric acid was negatively correlated with low-density lipoprotein(P < 0.05). Previous studies have found that SCFAs also have a role in regulating blood lipids. Propionic acid and butyric acid inhibited isoprenaline- and adenosine deaminase-stimulated lipolysis in the presence of phosphodiesterase (PDE3) inhibitors as well as propionic acid and butyric acid inhibited basal and insulin-stimulated de novo lipogenesis, which was associated with increased phosphorylation, thereby inhibiting the activity of HMG-CoA, the rate-limiting enzyme for fatty acid synthesis. Moreover, SCFAs have an effect on fat storage and mobilization in rat primary adipocytes, thus possibly contributing to healthier adipocyte contracts and improving energy metabolism and reducing circulating free fatty acids [36].
There are several limitations in this study. First, we found that gut microbiota and SCFAs were involved in IgAN, but the possible mechanisms and pathways have not been elucidated. Second, the sample size was relative small, and the multiple confounders were not controlled for all correlation analyses. Our data should be interpreted with caution, and further studies are needed to verify the role of SCFA produced by the intestinal flora in IgAN.