ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 27 September 2022

Sec. Microbiome in Health and Disease

Volume 12 - 2022 | https://doi.org/10.3389/fcimb.2022.639624

Altered gut microbiota and gut-derived p-cresyl sulfate serum levels in peritoneal dialysis patients

  • 1. Department of Nephrology, Zhongshan Hospital, Fudan University, Shanghai, China

  • 2. Shanghai Institute of Kidney Disease and Dialysis, Shanghai, China

  • 3. Shanghai Key Laboratory of Kidney disease and Blood Purification, Shanghai, China

  • 4. Shanghai Medical Center of Kidney, Shanghai, China

Abstract

Peritoneal dialysis (PD) is a renal replacement therapy for end-stage renal disease. Gut microbiota-derived uremic solutes, indoxyl sulfate (IS), p-cresyl sulfate (PCS), and trimethylamine-N-oxide (TMAO) accumulate in PD patients. The objective was to explore the gut microbiota and their influence on uremic toxins in PD patients and healthy controls (HC). Fecal samples were collected from PD patients (n = 105) and HC (n = 102). 16S rRNA gene regions were sequenced for gut microbiota analysis. IS, PCS, and TMAO levels were measured using HPLC-MS. PD patients exhibited lower alpha diversity and altered gut microbiota composition compared to HC. At the genus level, PD patients showed increased abundance of opportunistic pathogenic bacteria, and decreased abundance of beneficial bacteria. Three Operational Taxonomic Units discriminated PD patients from HC. Phenylalanine metabolism increased in PD, whereas tryptophan metabolism was unaltered. Low serum PCS did not necessarily mean healthier due to the loss of alpha diversity, increased Proteobacteria and opportunistic pathogenic bacteria. High serum PCS was mainly caused by elevated p-cresol-producing bacteria, enriched amino acid related enzymes, and enhanced sulfur metabolism, rather than declined residual renal function. In patients with different urine volumes, the gut microbiota alpha diversity and composition were unaltered, but serum IS and TMAO were significantly elevated in anuric patients. In conclusion, the gut microbiota abundance, composition, and function were altered in PD patients, which increased the PCS levels. We provided a better understanding of the microbiota-metabolite-kidney axis in PD patients. Targeting certain bacteria could decrease the PCS levels, whereas preserving the residual renal function could reduce the IS and TMAO levels.

Introduction

The gut microbiota is the largest microecosystem in our body (). Gut microbiota dysbiosis is related to many diseases including acute kidney injury and chronic kidney disease (CKD) (; ). Peritoneal dialysis (PD) is an effective renal replacement therapy administered to approximately 11% of end-stage renal disease patients undergoing dialysis. The long-term frequent inflow and outflow of glucose dialysate in the peritoneal cavity and uremia may change the abundance and composition of the gut microbiota. A few studies found that PD patients have altered gut microbial composition (Wang et al., 2012; Stadlbauer et al., 2017). However, the differences in the gut microbiota between PD patients and healthy controls (HC) are not well known.

Indoxyl sulfate (IS), P-cresyl sulfate (PCS), and trimethylamine-N-oxide (TMAO) are gut microbiota-derived uremic toxins, which can contribute to the progression and cardiovascular complications of CKD and PD (Zeisel and Warrier, 2017; ). Since IS and PCS are protein-bound uremic toxins, they are hardly removed through PD (Vanholder et al., 1999). Uremic toxins accumulate in PD patients due to reduced renal function and limited removal of uremic toxins through dialysis (Viaene et al., 2014). We hypothesized that the gut microbiota in PD could accelerate the generation of uremic toxins.

In this study, we aimed to compare the gut microbiota of PD patients with that of HC by conducting16S rRNA gene sequencing of fecal samples. We also investigated the correlations between gut microbiota, serum uremic toxins, and urine volumes in PD patients.

Methods

Participants

In this case-control study, we enrolled 105 continuous ambulatory PD patients using glucose dialysate (Dianeal, Baxter) with dialysis vintage > 6 months, treated in Zhongshan Hospital from February 2018 to July 2019, and 102 healthy controls (HC) who got their annual physical examination done in our hospital. All participants were Han nationality and older than 18 years old. Exclusion criteria included: (1) gastrointestinal diseases; (2) any specific drug use in the previous three months, including antibiotics, probiotics, prebiotics, synbiotics, proton pump inhibitors, and immunosuppressive agents; (3) severe liver diseases, hepatitis; (4) tumor; (5) other immunological or autoimmune disorders. The study was approved by the Ethical Committee of Zhongshan Hospital, Fudan University (Approval No.: B2017-108R). All participants provided written informed consent.

Fecal and serum sample collection

Stool samples were collected from the participants by spontaneous evacuation, transported to our laboratory on dry ice, divided into four equal parts of 200 mg, and immediately stored at −80°C for further analysis. Blood was collected in the morning pre-prandial. After centrifugation (3000×g, 10 min), the serum was immediately frozen at −80°C until use.

Clinical assessment

Serum creatinine, blood urea nitrogen, uric acid, albumin, hemoglobin, triglyceride, total cholesterol, calcium, phosphorus, and C-reactive protein (CRP) levels were measured using standard methods followed in Zhongshan Hospital. Body mass index (BMI) was calculated using the formula: BMI = kg/m2. The CKD-EPI equation was used to calculate the estimated glomerular filtration rate (eGFR) (). IS, PCS, and TMAO levels were measured using high performance liquid chromatography mass spectrometry (HPLC-MS) as previously described (). The normalized protein catabolic rate (nPCR) was calculated using the PD adequest 2.0 software (Baxter Healthcare, Norfolk, UK) as a measure of daily protein intake. Kt/V (kidney) was used to evaluate the renal excretory function, and Kt/V (peritoneum) was used to monitor PD treatment efficiency based on urea clearance (Wang and Wang, 2015). A questionnaire was administered to PD patients (n = 40) and HC (n = 20) to collect information about the number of bowel movements per week and the Bristol Stool Scale assessing constipation (harder stool).

16S rRNA microbial profiling analysis and functional annotation

Bacterial DNA was extracted from stool samples using the E.Z.N.A.® Stool DNA Kit following the manufacturer’s instructions (Omega Bio-tek, Inc., GA). The V3-V4 variable regions of the microbial 16S rRNA gene were amplified using the following primers: 341F (CCTACGGGNGGCWGCAG) and 805R (GACTACHVGGGTATCTAATCC). The products from different samples were indexed, mixed at equal ratios, and sequenced using an Illumina Miseq platform (Illumina Inc., USA) by Shanghai Mobio Biomedical Technology Co. Ltd, China ().

After raw data extraction using USEARCH 8.0, an average of 35727 ± 8763 reads per sample was obtained. The sequences were classified into Operational Taxonomic Units (OTUs) clusters based on 97% similarity (3% divergence). Each 16S rRNA gene sequence was annotated using the RDP Classifier (http://rdp.cme.msu.edu/), referring to the Silva (SSU123) 16S rRNA database.

Alpha diversities are presented as the Shannon, Chao, and Ace diversity indices using the “vegan” package in R. Principal Coordinate Analysis (PCoA) plots and non-metric multidimensional scaling (NMDS) analysis were generated to visualize the unweighted UniFrac distances using the QIIME pipeline in R. The differential taxa between the PD and HC groups were detected using the linear discriminant analysis (LDA) effect size (LEfSe) (LDA score = 2 as the cut-off value).

Kyoto Encyclopedia of Genes and Genomes (KEGG) () analysis was used to predict functional pathway information as categorized using the phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) algorithm ().

Statistical analysis

All continuous data are presented as mean ± SD. Categorical variables are expressed as percentages. The student’s t-test or one-way ANOVA with the Bonferroni multiple comparison test was adopted to compare continuous variables. The Chi-squared test was used to compare categorical variables. Univariate correlation was assessed using Spearman’s correlation analyses. A P-value of <0.05 (two tailed) was considered statistically significant. SPSS (version 23.0) and GraphPad Prism (version 7.0) were used for statistical analysis.

Results

Differences in the serum metabolomes of PD patients and HC

The participants in the PD and HC groups were matched for age, gender, and BMI. Compared with the HC, the renal function in the PD group was significantly decreased as shown by the serum creatinine and eGFR. IS, PCS, and TMAO levels were remarkably increased in PD patients. The inflammatory biomarker CRP was also enhanced in PD (Table 1).

Table 1

CharacteristicPD (n=105)HC (n=102)P-value
Age, years56.79 ± 13.7454.20 ± 11.250.139
Gender, male/female, n52/5351/511.000
Body mass index, kg/m223.02 ± 4.1622.90 ± 2.130.801
Serum creatinine, μmol/L857.46 ± 257.2774.29 ± 14.39<0.001
Blood urea nitrogen, mmol/L18.38 ± 4.924.51 ± 0.92<0.001
Uric acid, umol/L371.01 ± 96.69344.07 ± 70.550.023
eGFR, ml/min/1.73m25.13 ± 1.6690.70 ± 13.67<0.001
Albumin, g/L35.00 ± 4.8943.53 ± 2.71<0.001
RBC, ×1012/L3.52 ± 0.764.61 ± 0.49<0.001
Hemoglobin, g/L105.67 ± 19.80139.86 ± 16.86<0.001
WBC, ×109/L7.22 ± 2.195.55 ± 1.11<0.001
Platelet, ×109/L205.53 ± 80.08201.66 ± 46.10.672
Triglyceride, mmol/L1.94 ± 1.161.66 ± 0.950.061
Total cholesterol, mmol/L4.38 ± 1.134.78 ± 0.900.005
Calcium, mmol/L2.48 ± 1.982.28 ± 0.100.301
Phosphorus, mmol/L1.57 ± 0.451.17 ± 0.15<0.001
PCS, μg/mL20.68 ± 15.540.63 ± 0.81<0.001
IS, μg/mL27.00 ± 14.810.35 ± 0.26<0.001
TMAO, μg/mL5.14 ± 3.450.19 ± 0.28<0.001
CRP, mg/L7.53 ± 12.881.00 ± 1.56<0.001
Bristol Stool Scale score4.1 ± 1.324.2 ± 0.770.755
Number of weekly bowel movements11.50 ± 5.227.05 ± 2.01<0.001
Hypertension, n (%)100 (95.24)NO
Diabetes, n (%)23 (21.90)NO

Clinical characteristics of peritoneal dialysis patients and healthy controls.

eGFR, estimated glomerular filtration rate; RBC, red blood cells; WBC, white blood cells; PCS, p−cresyl sulfate; IS, Indoxyl sulfate; TMAO, trimethylamine-N-oxide; CRP, C-reactive protein. Values were presented as mean ± SD or n (%).

Decreased gut microbiota alpha-diversity in PD patients

The rarefaction curve showed that with the increasing number of sequences or samples, the number of OTUs reached a plateau, meaning the data were qualified and the sample size was sufficient (Figures 1A, B). The Venn diagram showed that 110 OTUs were unique for PD, 406 OTUs were unique for HC, and 1364 OTUs were shared by both groups (Figure 1C). The OTUs were significantly reduced in PD patients compared to those of HC (P<0.0001; Figure 1D). The Shannon (P <0.05), Chao (P <0.0001), and Ace indices (P <0.0001) were significantly lower in PD patients compared to that in HC (Figures 1E–G).

Figure 1

Altered gut microbiota composition in PD patients

To compare the beta diversity, the PCoA and NMDS analysis based on unweighted UniFrac distance analysis, displayed the microbiome space between the samples. A significant difference was observed between the PD and HC groups with the distribution of the OTUs from PC1 (30.98%), PC2 (20.20%), and PC3 (6.73%) (P<0.001, Figure 2A). The NMDS analysis also revealed a separation trend of PD and HC from NMDS1 and NMDS2 (P<0.001, Figure 2B).

Figure 2

Differences in gut microbiota abundance between PD and HC

The top 50 differentially abundant OTUs between PD patients and HC were presented in a heatmap (Supplementary Figure S1). Cladograms generated by LEfSe indicated the differences in the predominant bacteria between the PD and HC groups (Supplementary Figure S2). Furthermore, most gut microbiota fall into three phyla: Firmicutes, Bacteroidetes and Proteobacteria in PD (~50%, 28%, 14%, respectively) and HC (~62%, 24%, 7%, respectively) (Figure 3A). Proteobacteria (P <0.01), Fusobacteria (P <0.001), and Bacteria unclassified (P <0.001) were significantly increased, while Firmicutes (P <0.001), Tenericutes (P <0.001), and Cyanobacteria (P <0.001) were significantly decreased in PD patients compared to their levels in HC (Figure 3B). At the genus level, 16 and 27 genera were significantly increased and decreased, respectively, in PD compared to their levels in HC (LDA >2). We found a significant enrichment of opportunistic pathogens including Bacteroides, Escherichia Shigella, and Flavonifractor. Furthermore, a depletion in the abundance of beneficial bacteria, such as short-chain fatty acids (SCFA)-producing bacteria, including acetate-producing Bifidobacterium, butyrate-producing Faecalibacterium, Subdoligranulum, Roseburia, Erysipelotrichaceae, Lachnospiraceae, and Clostridiales () was observed in PD patients compared to that in HC (Figure 3C). Moreover, p-cresol-producing bacteria were enriched in PD, including Fusobacterium, Clostridium, and Eubacterium; Bacteroides and Fusobacterium are indole-producing bacteria (). They could lead to an increased production of PCS and IS. At the species levels, the similar results were observed (Supplementary Figure S3). Meanwhile, Blautia (r = 0.262, P = 0.008) and Fusobacterium (r = 0.233, P = 0.019) were associated with CRP. Thus, PD status influenced the relative quantity of certain gut microbiota and strengthened the importance of bacteria in systemic inflammation or accumulation of uremic toxins in PD patients.

Figure 3

Predictive function analysis

KEGG pathway analysis was conducted to predict bacterial functions in the PD and HC groups. At level three, we found that 22 predicted gut microbial functions, including oxidative phosphorylation, fructose and mannose metabolism, and pentose phosphate pathway, were significantly increased; while 34 functions, including DNA repair and recombination proteins, glycolysis and gluconeogenesis, were significantly decreased in PD patients compared with those of HC (P <0.05, LDA >2. Figure 4; Supplementary Data S1). The long-term effects of high glucose dialysate in the peritoneal cavity of PD patients may explain the disordered glucose metabolism. Notably, bacterial motility proteins, bacterial chemotaxis, flagellar assembly, and peptidoglycan biosynthesis were reduced in PD patients (Supplementary Data S1). Although IS was a tryptophan metabolite, no significant difference in tryptophan metabolism was observed between the two groups (P >0.05). Furthermore, the enriched phenylalanine metabolism (P = 0.043, LDA < 2) in PD confirmed that gut microbiota play a more important role in the metabolism of PCS than that of IS.

Figure 4

Validity of gut microbiota in discriminating PD

To detect unique OTUs markers of PD patients, we conducted fivefold cross-validation on a random forest model between 105 PD patients and 102 HC samples. Three OTUs (Pseudobutyrivibrio, Lachnospiraceae, and Erysipelotrichaceae) were selected as the optimal markers of PD with minor CV error (Figure 5A). In the training cohort, the probability of disease (POD) value calculated using the three optimal OTUs was significantly higher in PD patients (approached 1) compared to that of HC (approached 0, Figure 5B). The receiver operating characteristics (ROC) curve showed that the area under the curve (AUC) of the POD value was 0.9231 (Figure 5C). In the testing cohort, the POD value was also significantly increased in PD (Figure 5D), and the AUC was 0.8452 (Figure 5E). Although PD does not require specific gut microbiota as diagnostic markers, the high value of AUC suggested that PD could be distinguished from HC by the three OTUs with high specificity and sensitivity. Measuring the variation of the 3 OTUs described in Figure 5A may help to monitor the effectiveness of PD treatment (such as probiotics or synbiotics).

Figure 5

Correlation between the gut microbiome and PCS in PD patients

In PD patients, alpha-diversity, as measured by Shannon, Chao, Ace indices, and observed OTUs, was positively correlated with PCS in PD patients (Table 2; Figure 6A). There was a weak correlation between alpha-diversity and IS or phosphorus levels (r <0.3). The alpha-diversity did not correlate with TMAO levels and other variables. PCS was associated with 74 OTUs, IS with 20 OTUs, and TMAO with six OTUs in PD patients (Figure 6B). Among them, 14 OTUs were both associated with PCS and IS. Hence, the degree of microbiota disorder in PD patients was more closely related to PCS than that to IS and TMAO.

Table 2

Shannon indexChao indexAce indexObserved OTUs
rP-valuerP-valuerP-valuerP-value
Age0.1070.2780.1350.1700.0820.4070.1520.121
Dialysis vintage0.0050.963-0.0450.6510.0010.989-0.1050.288
BMI-0.1530.1200.0470.6370.0220.8260.0460.639
PCS0.3060.0020.512<0.0010.435<0.0010.509<0.001
IS0.1060.2800.2590.0080.2060.0350.2120.030
TMAO-0.0440.6580.1100.2630.0280.7760.0270.783
eGFR0.0390.694-0.0290.770-0.0500.6140.0170.862
24h urine volumns-0.0350.7230.1010.3030.0480.6290.1350.170
Creatinine-0.0310.7560.0210.8280.0200.839-0.0280.776
Calcium0.0410.6750.0930.3440.0900.3630.0950.337
Phosphorus-0.2270.020-0.0620.528-0.0800.420-0.0760.439
Serum albumin0.0090.9270.1000.3120.0980.3200.0950.336
CRP-0.0200.8430.0500.6200.0640.5270.0290.771

Correlation between alpha diversity and variables.

OTUs, operational taxonomic units; BMI, body mass index; PCS, p-cresyl sulfate; IS, Indoxyl sulfate; TMAO, trimethylamine-N-oxide; eGFR, estimated glomerular filtration rate; CRP, C-reactive protein.

Figure 6

The metabolic pathway of PCS is shown in Figure 6D. To clarify the relationship between PCS and gut microbiota, we divided the PD patients into PCS-low (L) (PCS ≤ 18.7 μg/ml, n = 53) and PCS-high (H) (PCS > 18.7 μg/ml, n = 52) subgroups according to the PCS median concentration and compared the clinical data in these subgroups (Table 3). Here, two subgroups had similar protein intake as evaluated by nPCR and albumin. No difference was observed in the eGFR, urine volumes, Kt/V (kidney), and Kt/V (peritoneum), implying similar renal excretory function and PD treatment efficiency. The comorbidities and medicinal regimens also showed no significant difference between two subgroups (Supplementary Table S1).

Table 3

CharacteristicPCS-L (n=53)PCS-H (n=52)P-value
Age, years55.51 ± 12.7458.1 ± 14.690.337
Gender, male/female, n28/2524/280.560
Body mass index, kg/m223.64 ± 4.2322.38 ± 4.030.117
Dialysis vintage, month27.62 ± 24.3229.63 ± 33.360.724
Kt/V (kidney)0.55 ± 0.560.63 ± 0.660.522
Kt/V (pritoneum)1.50 ± 0.471.55 ± 0.370.554
24h urine volumn, L0.60 ± 0.530.73 ± 0.660.279
eGFR, ml/min/1.73m24.92 ± 1.635.34 ± 1.680.195
Serum creatinine, μmol/L901.13 ± 273.60812.94 ± 233.760.079
Blood urea nitrogen, mmol/L18.37 ± 5.2418.39 ± 4.610.981
Uric acid, μmol/L374.34 ± 75.91367.62 ± 114.740.723
Hemoglobin, g/L104.83 ± 21.44106.52 ± 18.140.664
Serum albumin, g/L34.51 ± 4.6635.50 ± 5.100.301
nPCR0.89 ± 0.200.95 ± 0.280.208
Triglyceride, mmol/L2.07 ± 1.271.80 ± 1.020.245
Total cholesterol, mmol/L4.36 ± 1.214.39 ± 1.050.923
Calcium, mmol/L2.25 ± 0.182.72 ± 2.800.226
Phosphorus, mmol/L1.63 ± 0.411.51 ± 0.490.168
CRP, mg/L6.21 ± 12.248.73 ± 13.420.327
Number of observed OTUs239.19 ± 49.53292.81 ± 73.65<0.001
Shannon index3.15 ± 0.463.37 ± 0.590.040
Chao index310.65 ± 63.26376.02 ± 87.74<0.001
Ace index324.62 ± 68.99382.68 ± 84.66<0.001
PCS, μg/mL8.39 ± 5.3633.2 ± 12.09<0.001
IS, μg/mL22.97 ± 13.2331.11 ± 15.320.004
TMAO, μg/mL5.02 ± 3.255.26 ± 3.660.723

Clinical data of PCS-L and PCS-H subgroups.

eGFR, estimated glomerular filtration rate; PCS, p-cresyl sulfate; nPCR, normalized protein catabolic rate; CRP, C-reactive protein; PCS, p-cresyl sulfate; IS, Indoxyl sulfate; TMAO, trimethylamine-N-oxide. Values were presented as mean ± SD or n.

Notably, the alpha diversity of gut microbiota was significantly reduced in the PCS-L subgroup compared to PCS-H subgroup (Table 3). A Venn diagram showed that 111 OTUs were unique in the PCS-L subgroup, and 325 OTUs were unique in the PCS-H subgroup (Figure 6C), suggesting that patients with low serum PCS were not necessarily healthier due to a diminished bacterial abundance. Gut microbiota composition was also significantly altered in the two subgroups as demonstrated using PCoA and NMDS analysis (P<0.001, Supplementary Figure S4). The top three abundant phyla were Firmicutes, Bacteroidetes, and Proteobacteria in the PCS-L (~49%, 29%, 16%, respectively) and PCS-H subgroups (~52%, 27%, 11%, respectively, Figure 6E). Interestingly, the ratio of the three phyla in the PCS-H subgroup was much similar to that of HC. The significantly increased genera in PCS-H might contribute to increased PCS (Figure 6F). The p-cresol-producing bacteria Ruminococcus, Roseburia, and Anaerococcus enriched in PCS-H could accelerate PCS biosynthesis. Remarkably, the opportunistic pathogenic bacterium Escherichia Shigella was elevated in PCS-L. KEGG pathway analysis predicted bacterial functions in two subgroups at levels 1–3 (Supplementary Figure S5). As expected, amino acid-related enzymes (P = 0.003, LDA>2) and sulfur metabolism (P = 0.037, LDA>2) were increased in the PCS-H subgroup, which can explain the enhanced serum PCS and IS levels.

Furthermore, comparisons of the gut microbiota in terms of the IS and TMAO levels, showed that gut microbiome abundance and composition did not differ between the IS-low and IS-high subgroups divided by the median IS concentration (25.4 μg/ml). No difference was observed between the TMAO-low and TMAO-high subgroups divided by the median TMAO concentration (4.1 μg/ml), either. Thus, the disordered gut microbiota in PD was associated with serum PCS, rather than with IS or TMAO.

Gut microbiome and urine volume correlations in PD patients

The mean of endogenous creatinine clearance rate and urea clearance rate, or urine volume was used to estimate residual renal function (RRF) (). RRF is usually assessed by urine volumes in uremic patients. PCS levels showed no correlation with urine volumes (r = 0.054, P = 0.585). We assigned 105 PD patients into three subgroups according to urine volumes and investigated the levels of uremic toxins. Anuria was defined as 24-h urine volumes <100 mL. PD patients with 24-h urine volumes of <100 mL (n = 34), 100–1000 mL (n = 31), and >1000 mL (n = 40) were assigned to the U1, U2, and U3 groups, respectively. IS was found to be remarkably increased in the U1 group compared to U2 and U3 groups (both P <0.0001, Figure 7A). TMAO levels were higher in the U1 group compared to U3 group (P <0.05, Figure 7C). However, the PCS levels did not significantly differ among three groups (P >0.05, Figure 7B).

Figure 7

Significant differences in the alpha and beta diversities were not detected among the U1, U2, and U3 groups (Supplementary Data S2; Supplementary Figure S6). A Venn diagram showed that 108, 69 and 149 OTUs were unique in the U1, U2 and U3 groups, respectively (Figure 7D). At the genus level, the levels of Flavonifractor (r = −0.321, P <0.05) and Bilophila (r = −0.231, P <0.05) gradually decreased with the increasing urine volumes in the U1, U2, and U3 groups (Figure 7E). Whether Flavonifractor and Bilophila can produce IS or TMAO has never been reported. Hence, there was no correlation between urine volumes and microbiota diversity.

Discussion

Many studies have concentrated on the disordered gut microbiota in CKD stages 1–5 or in hemodialysis patients, however, the relationship between gut microbiota and uremic toxins in PD patients has been rarely reported.

In this study, first, we compared the gut microbiota between PD patients and HC (summarized in Supplementary Figure S7). We found that the alpha diversity was significantly decreased in PD patients compared to HC, consistent with previous studies (; Stadlbauer et al., 2017). Many factors such as dietary restrictions, drugs, dialysis treatment, and malnutrition may inhibit the growth of certain bacteria. The genera enriched in PD patients were mainly opportunistic pathogens, and p-cresol and indole-producing bacteria, proving a correlation between uremic toxins and gut microbiota in PD. SCFAs including acetate, propionate, and butyrate have protective effects in inflammatory responses, stabilizing the gut barrier, and attenuating bacterial translocation (; Routy et al., 2018). Supplying SCFAs could reduce local and systemic inflammation and improve renal function (; ). The decreased acetate producer Bifidobacterium and butyrate producers Faecalibacterium and Subdoligranulum, suggest a diminished ability to generate SCFAs and a potential gut barrier dysfunction and systemic inflammation in PD patients. Furthermore, three OTUs could accurately differentiate PD patients from HC.

Second, we compared the gut microbiota in different PD subgroups. Increased PCS levels were directly associated with cardiovascular events and mortality, and the inflammatory marker interleukin-6 and oxidative stress marker glutathione peroxidase (; Rossi et al., 2014). The gut microbiota directly increased serum PCS levels with increased alpha diversity, altered microbiota composition, higher abundance of p-cresol-producing bacteria, enriched amino acid related enzymes, and elevated sulfur metabolism. A previous study showed that PCS and TMAO levels were positively correlated with the Shannon index in HC (Wilmanski et al., 2019). We proved that the higher alpha diversity was only associated with the serum microbial metabolite PCS in PD patients. Targeting p-cresol-producing bacteria may diminish PCS in PD patients resulting in therapeutic benefits. Combined supplementation with prebiotics and probiotics could reduce serum PCS in CKD patients (Rossi et al., 2016). However, a randomized crossover study found that inulin-type fructan intervention could not decline PCS in PD patients (). Thus, it is imperative to reduce PCS levels by targeting certain gut microbiota in PD patients.

We found that low serum PCS was partly due to reduced alpha diversity and abundance of p-cresol-producing bacteria. Microbial diversity could be influenced by inflammatory and nutritional status to some extent (), and primarily by intestinal physiology (). Although reduced alpha diversity leads to lower PCS levels, it could be related to diseases such as inflammatory bowel disease (; Wilmanski et al., 2019). The phylum Proteobacteria was significantly enriched in patients with low PCS levels. Increased Proteobacteria is a potential diagnostic signature of microbial dysbiosis and risk of inflammation and disease (Shin et al., 2015; ). The elevated opportunistic pathogenic bacteria also associated with the occurrence of peritonitis (; ). Those suggested that patients with lower PCS levels were not always healthier. Their gut microbiota should be analyzed. If their low PCS were due to a lack of α-diversity or disordered bacteria, they may be worse than patients with high PCS.

Finally, we investigated whether decreased renal function could increase serum PCS levels in anuric PD patients. Some researchers found that serum PCS increased with the decreased eGFR or residual renal function in CKD and PD patients (; Viaene et al., 2014). In contrast, others found that increasing serum PCS did not parallel with declining residual renal function in PD patients (; ), while IS and TMAO were reported to have a strong correlation with eGFR (; ). Consistently, we found that PCS was not correlated with urine volumes in PD patients, but IS and TMAO were still related to residual renal function, suggesting the validity of our enrolled patients. The kidney clearance rate of IS was three times that of PCS (). Serum PCS levels did not change with urine volumes possibly due to its lower kidney clearance rate compared to IS, as well as the different diets and life styles among three groups. The gut microbiota did not significantly change among the U1, U2, and U3 groups. We propose that kidney excretion has a greater influence on serum IS and TMAO levels than their generation by gut microbiota. A recent study showed that impaired kidney function was the main contributor to the increased serum IS levels at different CKD stages (). Our findings proved that increased IS and TMAO levels were mainly due to the loss of residual renal function in PD patients. The need to focus more on the production mechanisms of PCS rather than on its renal elimination has been emphasized through this study.

Our study had several limitations. Compared with 16S rRNA analysis, shotgun metagenome sequencing can provide a more accurate description of the gut microbiota composition and function. Furthermore, we did not include fiber intake and lifestyle assessments or measured the PCS level in urine, stool, and dialysate, which may have limited the interpretations of our results.

Conclusion

The gut microbiota abundance, composition, and function were altered in PD patients. Our study provided novel insights into the relationship between gut microbiota and uremic toxins in PD patients. The gut microbiota might provide therapeutic possibilities for preventing the generation of PCS. Moreover, preserving residual renal function could reduce the IS and TMAO levels.

Funding

This work was funded by the National Natural Science Foundation of China (81970667) and the Development Plan of Top Young Talents in Shanghai (2018).

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The raw Illumina read data for all samples were deposited in the European Bioinformatic Institute database under the accession code PRJNA682853 and PRJNA687563.

Ethics statement

The studies involving human participants were reviewed and approved by Zhongshan Hospital, Fudan University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

MB designed the experiments and wrote the manuscript. MB and PZ collected samples and analyzed data. SG was involved in analyzing the data. JZ and JJ helped interpret data. XD and XY designed and supervised all experiments, and edited the manuscript. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.639624/full#supplementary-material

References

  • 1

    Andrade-OliveiraV.AmanoM. T.Correa-CostaM.CastoldiA.FelizardoR. J.de AlmeidaD. C.et al. (2015). Gut bacteria products prevent AKI induced by ischemia-reperfusion. J. Am. Soc. Nephrol.26 (8), 18771888. doi: 10.1681/asn.2014030288

  • 2

    BammensB.EvenepoelP.VerbekeK.VanrenterghemY. (2005). Time profiles of peritoneal and renal clearances of different uremic solutes in incident peritoneal dialysis patients. Am. J. Kidney Dis.46 (3), 512519. doi: 10.1053/j.ajkd.2005.05.016

  • 3

    CaoX. S.ChenJ.ZouJ. Z.ZhongY. H.TengJ.JiJ.et al. (2015). Association of indoxyl sulfate with heart failure among patients on hemodialysis. Clin. J. Am. Soc. Nephrol.10 (1), 111119. doi: 10.2215/cjn.04730514

  • 4

    CaporasoJ. G.LauberC. L.WaltersW. A.Berg-LyonsD.HuntleyJ.FiererN.et al. (2012). Ultra-high-throughput microbial community analysis on the illumina HiSeq and MiSeq platforms. Isme. J.6 (8), 16211624. doi: 10.1038/ismej.2012.8

  • 5

    CosolaC.RocchettiM. T.CupistiA.GesualdoL. (2018). Microbiota metabolites: Pivotal players of cardiovascular damage in chronic kidney disease. Pharmacol. Res.130, 132142. doi: 10.1016/j.phrs.2018.03.003

  • 6

    Crespo-SalgadoJ.VehaskariV. M.StewartT.FerrisM.ZhangQ.WangG.et al. (2016). Intestinal microbiota in pediatric patients with end stage renal disease: a Midwest pediatric nephrology consortium study. Microbiome4 (1), 50. doi: 10.1186/s40168-016-0195-9

  • 7

    De AndrésJ.ManzanoS.GarcíaC.RodríguezJ. M.Espinosa-MartosI.JiménezE. (2018). Modulatory effect of three probiotic strains on infants' gut microbial composition and immunological parameters on a placebo-controlled, double-blind, randomised study. Benef. Microbes9 (4), 573584. doi: 10.3920/bm2017.0132

  • 8

    De SordiL.KhannaV.DebarbieuxL. (2017). The gut microbiota facilitates drifts in the genetic diversity and infectivity of bacterial viruses. Cell Host Microbe22 (6), 801808.e803. doi: 10.1016/j.chom.2017.10.010

  • 9

    DoğukanA.OymakO.TaskapanH.CinarS.TokgözB.UtasC. (2000). Shigella sonnei peritonitis in a patient on CAPD. Perit. Dial Int.20 (6), 806.

  • 10

    DuY.TangG.YuanW. (2020). Suppression of HDAC2 by sodium butyrate alleviates apoptosis of kidney cells in db/db mice and HG−induced NRK−52E cells. Int. J. Mol. Med.45 (1), 210222. doi: 10.3892/ijmm.2019.4397

  • 11

    EmalD.RampanelliE.StrooI.ButterL. M.TeskeG. J.ClaessenN.et al. (2017). Depletion of gut microbiota protects against renal ischemia-reperfusion injury. J. Am. Soc. Nephrol.28 (5), 14501461. doi: 10.1681/asn.2016030255

  • 12

    FengX.YangX.YiC.GuoQ.MaoH.JiangZ.et al. (2014). Escherichia coli peritonitis in peritoneal dialysis: the prevalence, antibiotic resistance and clinical outcomes in a south China dialysis center. Perit. Dial Int.34 (3), 308316. doi: 10.3747/pdi.2013.00012

  • 13

    GrypT.De PaepeK.VanholderR.KerckhofF. M.Van BiesenW.Van de WieleT.et al. (2020). Gut microbiota generation of protein-bound uremic toxins and related metabolites is not altered at different stages of chronic kidney disease. Kidney Int.97 (6), 12301242. doi: 10.1016/j.kint.2020.01.028

  • 14

    KanehisaM.FurumichiM.TanabeM.SatoY.MorishimaK. (2017). KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res.45 (D1), D353d361. doi: 10.1093/nar/gkw1092

  • 15

    LangilleM. G.ZaneveldJ.CaporasoJ. G.McDonaldD.KnightsD.ReyesJ. A.et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol.31 (9), 814821. doi: 10.1038/nbt.2676

  • 16

    LeveyA. S.StevensL. A.SchmidC. H.ZhangY. L.CastroA. F.FeldmanH. I.et al. (2009). A new equation to estimate glomerular filtration rate. Ann. Intern. Med.150 (9), 604612. doi: 10.7326/0003-4819-150-9-200905050-00006

  • 17

    LinT. Y.HungS. C.Association of subjective global assessment of nutritional status with gut microbiota in hemodialysis patients: A case-control study. Nephrol. Dial Transplant. (2021) 36 (6), 11041111. doi: 10.1093/ndt/gfaa019

  • 18

    LiL.XiongQ.ZhaoJ.LinX.HeS.WuN.et al. (2020). Inulin-type fructan intervention restricts the increase in gut microbiome-generated indole in patients with peritoneal dialysis: a randomized crossover study. Am. J. Clin. Nutr.111 (5), 10871099. doi: 10.1093/ajcn/nqz337

  • 19

    LuW.RenC.HanX.YangX.CaoY.HuangB. (2018). The protective effect of different dialysis types on residual renal function in patients with maintenance hemodialysis: A systematic review and meta-analysis. Med. (Baltimore)97 (37), e12325. doi: 10.1097/md.0000000000012325

  • 20

    MeijersB. K.ClaesK.BammensB.de LoorH.ViaeneL.VerbekeK.et al. (2010). P-cresol and cardiovascular risk in mild-to-moderate kidney disease. Clin. J. Am. Soc. Nephrol.5 (7), 11821189. doi: 10.2215/cjn.07971109

  • 21

    PascalV.PozueloM.BorruelN.CasellasF.CamposD.SantiagoA.et al. (2017). A microbial signature for crohn's disease. Gut66 (5), 813822. doi: 10.1136/gutjnl-2016-313235

  • 22

    PelletierC. C.CroyalM.EneL.AguesseA.Billon-CrossouardS.KrempfM.et al. (2019). Elevation of trimethylamine-N-Oxide in chronic kidney disease: Contribution of decreased glomerular filtration rate. Toxins (Basel)11 (11), 635. doi: 10.3390/toxins11110635

  • 23

    PhamN. M.RechtN. S.HostetterT. H.MeyerT. W. (2008). Removal of the protein-bound solutes indican and p-cresol sulfate by peritoneal dialysis. Clin. J. Am. Soc. Nephrol.3 (1), 8590. doi: 10.2215/cjn.02570607

  • 24

    PoesenR.ViaeneL.VerbekeK.ClaesK.BammensB.SprangersB.et al. (2013). Renal clearance and intestinal generation of p-cresyl sulfate and indoxyl sulfate in CKD. Clin. J. Am. Soc. Nephrol.8 (9), 15081514. doi: 10.2215/cjn.00300113

  • 25

    ReeseA. T.DunnR. R. (2018). Drivers of microbiome biodiversity: A review of general rules, feces, and ignorance. mBio9 (4), e01294–18. doi: 10.1128/mBio.01294-18

  • 26

    RenZ.FanY.LiA.ShenQ.WuJ.RenL.et al. (2020). Alterations of the human gut microbiome in chronic kidney disease. Adv. Sci. (Weinh)7 (20), 2001936. doi: 10.1002/advs.202001936

  • 27

    RizzattiG.LopetusoL. R.GibiinoG.BindaC.GasbarriniA. (2017). Proteobacteria: A common factor in human diseases. BioMed. Res. Int.2017, 9351507. doi: 10.1155/2017/9351507

  • 28

    RossiM.CampbellK. L.JohnsonD. W.StantonT.VeseyD. A.CoombesJ. S.et al. (2014). Protein-bound uremic toxins, inflammation and oxidative stress: A cross-sectional study in stage 3-4 chronic kidney disease. Arch. Med. Res.45 (4), 309317. doi: 10.1016/j.arcmed.2014.04.002

  • 29

    RossiM.JohnsonD. W.MorrisonM.PascoeE. M.CoombesJ. S.ForbesJ. M.et al. (2016). Synbiotics easing renal failure by improving gut microbiology (SYNERGY): A randomized trial. Clin. J. Am. Soc. Nephrol.11 (2), 223231. doi: 10.2215/cjn.05240515

  • 30

    RoutyB.GopalakrishnanV.DaillèreR.ZitvogelL.WargoJ. A.KroemerG. (2018). The gut microbiota influences anticancer immunosurveillance and general health. Nat. Rev. Clin. Oncol.15 (6), 382396. doi: 10.1038/s41571-018-0006-2

  • 31

    ShinN. R.WhonT. W.BaeJ. W. (2015). Proteobacteria: Microbial signature of dysbiosis in gut microbiota. Trends Biotechnol.33 (9), 496503. doi: 10.1016/j.tibtech.2015.06.011

  • 32

    StadlbauerV.HorvathA.RibitschW.SchmerböckB.SchilcherG.LemeschS.et al. (2017). Structural and functional differences in gut microbiome composition in patients undergoing haemodialysis or peritoneal dialysis. Sci. Rep.7 (1), 15601. doi: 10.1038/s41598-017-15650-9

  • 33

    VanholderR.De SmetR.LesafferG. (1999). P-cresol: A toxin revealing many neglected but relevant aspects of uraemic toxicity. Nephrol. Dial Transplant.14 (12), 28132815. doi: 10.1093/ndt/14.12.2813

  • 34

    ViaeneL.MeijersB. K.BammensB.VanrenterghemY.EvenepoelP. (2014). Serum concentrations of p-cresyl sulfate and indoxyl sulfate, but not inflammatory markers, increase in incident peritoneal dialysis patients in parallel with loss of residual renal function. Perit. Dial Int.34 (1), 7178. doi: 10.3747/pdi.2012.00276

  • 35

    WangI. K.LaiH. C.YuC. J.LiangC. C.ChangC. T.KuoH. L.et al. (2012). Real-time PCR analysis of the intestinal microbiotas in peritoneal dialysis patients. Appl. Environ. Microbiol.78 (4), 11071112. doi: 10.1128/aem.05605-11

  • 36

    WangL.WangT. (2015). Adequacy of peritoneal dialysis: Kt/V revisited. Eur. Rev. Med. Pharmacol. Sci.19 (7), 12721275.

  • 37

    WilmanskiT.RappaportN.EarlsJ. C.MagisA. T.ManorO.LovejoyJ.et al. (2019). Blood metabolome predicts gut microbiome α-diversity in humans. Nat. Biotechnol.37 (10), 12171228. doi: 10.1038/s41587-019-0233-9

  • 38

    ZeiselS. H.WarrierM. (2017). Trimethylamine n-oxide, the microbiome, and heart and kidney disease. Annu. Rev. Nutr.37, 157181. doi: 10.1146/annurev-nutr-071816-064732

Summary

Keywords

Uremic toxins, Peritoneal Dialysis, residual renal function (RRF), Gut micobiota, p-cresyl sulfate (PCS)

Citation

Bao M, Zhang P, Guo S, Zou J, Ji J, Ding X and Yu X (2022) Altered gut microbiota and gut-derived p-cresyl sulfate serum levels in peritoneal dialysis patients. Front. Cell. Infect. Microbiol. 12:639624. doi: 10.3389/fcimb.2022.639624

Received

10 December 2020

Accepted

07 September 2022

Published

27 September 2022

Volume

12 - 2022

Edited by

Alain Stintzi, University of Ottawa, Canada

Reviewed by

Ferdinand Xiankeng Choong, Karolinska Institutet (KI), Sweden; Maria Teresa Rocchetti, University of Foggia, Italy

Updates

Copyright

*Correspondence: Xiaofang Yu, ; Xiaoqiang Ding,

†These authors have contributed equally to this work

This article was submitted to Microbiome in Health and Disease, a section of the journal Frontiers in Cellular and Infection Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics