Baseline Characteristics of Patients
In our study, a total of 80 patients were enrolled: 20 were healthy subjects and 60 were CKD patients. The average age of 20 healthy control (12 females, 8 males) was 55.15±13.77 years. The average age of 60 CKD patients (30 females, 30 males) was 55.18±14.80 years. In the total 60 CKD patients, the average age of stage 3 CKD patients (n=22) was 58.27±14.37 years, the average age of stage 4 CKD patients (n=17) was 54.35±15.09 years and the average age of stage 5 CKD patients (n=21) was 52.62±15.13 years. The parameters of healthy subjects included age, BMI, blood pressure, Scr, eGFR, BUN, UA, ALB, CRP, serum cholesterol, triglyceride, ALT, AST, K, Ca, P, Na ,CL and haemoglobin. Detailed informations about healthy people were listed in Table 1. The 60 CKD patients were divided into three stages: stage 3, stage 4, and stage 5. The clinical parameters of CKD patients included age, BMI, blood pressure, Scr, eGFR, BUN, UA, ALB, 24-hUTP CRP, serum cholesterol, triglyceride, ALT, AST, K, Ca, P, Na ,CL and haemoglobin. The levels of Scr were significantly increased from CKD 3 to CKD 5 (138.89±33.26 ,259.15±66.25, and 647.43±279.31, P = 0.000). eGFR was significantly decreased from CKD 3 to CKD 5 (43.57±8.66, 21.65±3.84, and 7.87±3.08, P = 0.000). The levels of BUN were significantly increased from CKD 3 to CKD 5 (9.38±3.79, 15.88±3.36, and 26.17±9.29, P = 0.000). The levels of Ca significantly decreased from CKD 3 to CKD 5 (2.27±0.23, 2.16±0.17, and 2.03±0.27, P = 0.004). The levels of P were significantly increased from CKD 3 to CKD 5 (1.28±0.32, 1.42±0.29, and 2.04±0.75, P = 0.000). The levels of haemoglobin decreased from CKD 3 to CKD 5 (12.33±2.12, 10.75±1.47, and 9.15±1.85, P = 0.000).
Gut Microbiota Structure and Alpha Diversity Analysis
A total of 80 patient samples were enrolled to test gut microbiota. The Venn graph results indicated that the gut microbiota in the two groups (healthy subjects and CKD patients) had 2351 OTUs in total (Figure 1). Among the two groups, the healthy subjects group had 1076 OTUs and the CKD group had 2259 OTUs. The rarefaction curve tended to be flat, indicating that the sequencing data volume was reasonable (Figure 2). The rank abundance curve showed that the sample size was sufficient (Figure 3). Within a certain range, with the increase of the sample size, the curve rises sharply, indicating that the sample size is insufficient and needs to be increased; On the contrary, the flatit carve indicates that the sample size is sufficient for data analysis. More than 10 samples are meaningful.
Next, We evaluated the alpha diversity index of the healthy subjects group, CKD 3 group, CKD 4 group, and CKD 5 group. As shown in Figure 4, Shannon analyses performed that the diversity of gut microbiota in the CKD 3 group was lower than that in the other groups in our study. Therefore, decreased diversity of gut microbiota might be the characters of CKD 3. It might be of enormus significance for indicating the progression and development of CKD.
Taxonomic Analysis
We investigated the distribution of phylum, class, order, family, and genus indifferent groups with taxonomic analysis. The phyla in different groups included bacteroidota, firmicutes, proteobacteria, actinobacteriota, fusobacteriota, verrucomicrobiota, desulfobacterota and cyanobacteria(Figure 5(A)). The classes in different groups included bacteroidia, clostridia, gammaproteobacteria, negativicutes, actinobacteria, fusobacteriia, bacilli, verrucomicrobiae, coriobacteriia, alphaproteobacteria and desulfovibrionia(Figure 5(B)). The orders in different groups included bacteroidales, lachnospirales, lachnospirales, enterobacterales, bifidobacteriales, selenomonadales, fusobacteriales, lactobacillales, verrucomicrobiales, christensenellales, vurkholderiales, pasteurellales, acidaminococcales, erysipelotrichales, coriobacteriales, rhodospirillales, tissierellales, clostridiales and monoglobales(Figure 5(C)). The families in different groups included bacteroidaceae, lachnospiraceae, prevotellaceae, ruminococcaceae, enterobacteriaceae, bifidobacteriaceae, fusobacteriaceae, veillonellaceae, tannerellaceae, selenomonadaceae, lactobacillaceae, akkermansiaceae, christensenellaceae, sutterellaceae, oscillospiraceae, pasteurellaceae, rikenellaceae, acidaminococcaceae and streptococcaceae(Figure 5(D)). The genera in different groups included bacteroides, faecalibacterium, prevotella, blautia, roseburia, megamonas, klebsiella, shigella, agathobacter, bifidobacterium, ruminococcus, parabacteroides, subdoligranulum, dorea, alistipes, phascolarctobacterium, lachnospira, lachnoclostridium, fusicatenibacter, fusobacterium and megasphaera(Figure 5(E)). Furthermore, taxonomic analysis in different samples were conducted. The phyla included firmicutes, bacteroidota proteobacteria and actinobacteriota(Figure 6(A)). The classes included clostridia, bacteroidia, negativicutes, gammaproteobacteria, bacilli, actinobacteria and coriobacteriia(Figure 6(B)). The orders included lachnospirales, bacteroidales, selenomonadales, tissierellales, lactobacillales and enterobacterales(Figure 6(C)). The families in different groups included lachnospiraceae, bacteroidaceae, veillonellaceae, peptostreptococcaceae and enterobacteriaceae(Figure 6(D)). The genera in different groups included bacteroides, blautia, roseburia, veillonella, lachnoclostridium, klebsiella, clostridioides, lactobacillus and parasutterella(Figure 6(E)).
Figure 6 Taxonomic analyses in all samples. (A) The phyla in all samples. (B) The classes in all samples. (C) The orders in all samples. (D) The families in all samples. (E) The genera in all samples.The figure showed the species information with relative abundance above 1%.
Diversity of gut microbiota among different CKD patients
We divided 60 CKD patients into the male group (n=30) and the female group (n=30). There were 1741 OTUs in the male group and 1778 OTUs in the female group, and 1440 OTUs were shared(Figure7). Then, by evaluating the α diversity index of gut microbiota between the male and female groups, it was found that the species and diversity of gut microflora in the male group were lower than that in the female group(Figure8).
Secondly, we divided CKD patients into three groups according to 24hUTP, and evaluated α diversity index(Figure 9). The results showed that the diversity index of gut microflora was highest when the 24hUTP was less than 1g/d, and was higher than the 24hUTP 1g/d-3.5g/d group and the 24hUTP>3.5g/d group.
We evaluated the association between clinical parameters and gut microbiota in 60 patients with CKD. The top 20 genera in 60 CKD patients were calculated by taxonomic analysis. The top 20 genera included bacteroides, faecalibacterium, prevotella, blautia, escherichia-Shigella, bifidobacterium, roseburia, fusobacterium, parabacteroides, agathobacter, akkermansia, megamonas, klebsiella, uncultered, ruminococcus_torques_group, ruminococcus, subdoligranulum, megasphaera, lachnoclostridium, alistipes(Figure10). The results performed that Scr was positively correlated with bacteroides(r=0.271, P=0.036) and blautia(r=0.435, P=0.0005), while was negatively correlated with bifidobacterium(r=-0.325, P=0.011). BUN was positively correlated with blautia(r=0.461, P=0.0002), while was negatively correlated with prevotella(r=-0.271, P=0.036), bifidobacterium(r=-0.349, P=0.006) and megamonas(r=-0.368, P=0.004). UA was negatively correlated with subdoligranulum(r=-0.286, P=0.027). ALB was negatively correlated with bacteroides(r=-0.308, P=0.017). 24hUTP was positively correlated with blautia(r=0.337, P=0.009). Serum Cholesterol was positively correlated with bacteroides(r=0.307, P=0.021), blautia(r=0.284, P=0.034), parabacteroides(r=0.296, P=0.027) and megasphaera(r=0.294, P=0.028). Haemoglobin was positively correlated with prevotella(r=0.326, P=0.011) and megamonas(r=0.369, P=0.0037), while was negatively correlated with bacteroides(r=-0.297, P=0.021), blautia(r=-0.331, P=0.0098) and klebsiella(r=-0.270, P=0.037). eGFR was negatively correlated with bacteroides(r=-0.258, P=0.046) and blautia(r=-0.45, P=0.0003), while was positively correlated with bifidobacterium(r=0.351, P=0.006), megamonas(r=0.279, P=0.0308). Some kinds of gut microbiota including bacteroides, blautia, parabacteroides, megasphaera and klebsiella might be detrimental factors in CKD, while other kinds of gut microbiota including bifidobacterium, prevotella, megamonas and subdoligranulum might be beneficial factors in CKD.