1. Introduction
More than one trillion microbes are found in the human body [1], and they encode 3.3 million genes, which is 150-fold higher than the number of genes encoded by the human genome [2]. By creating a stable environment for symbiotic interactions with the host, these microbes play critical metabolic roles in maintaining host health and homeostasis [2,3,4,5]. In chronic kidney disease (CKD), however, the influx of urea, uric acid, and oxalate from the gastrointestinal tract increases as kidney function deteriorates, which can perturb the balance between commensals and pathogens and lead to gut dysbiosis [2,6]. In parallel, the proliferation of urease-producing bacteria may result in the loss of intestinal barrier integrity and a sequential increase of bacterial translocation to the systemic circulation, leading to the deterioration of kidney function [2,7,8].
Gut dysbiosis due to kidney impairment is specifically associated with an alteration in uremic toxin metabolism [9,10]. Various uremic toxins, including p-cresol and indole derivatives, produced by dietary protein metabolism in humans have been identified [11]. The levels of p-cresol and indole derivatives in serum are negatively correlated with kidney function and these protein-bound solutes are inefficiently eliminated by dialysis [12,13,14]. They are also associated with decreased kidney function, cardiovascular disease, and mortality in patients with pre-dialytic CKD [11,12,15,16]. In addition to such amino acid-derived metabolites, trimethylamine-N-oxide (TMAO), which is produced by the microbial degradation of dietary quaternary amines, including choline, betaine, and carnitine, is an important uremic toxin related to the progression of CKD, atherosclerosis, heart attack, and metabolic syndrome [17,18,19].
Although various studies have reported the relationship between gut dysbiosis and kidney function impairment, as well as the involvement of gut microbiota in uremic toxin metabolism [8,19,20], few studies have examined the consequential linkages among gut microbiota, uremic toxins, and kidney dysfunction. This study was conducted to analyze the interrelationships among the gut microbiota, uremic metabolites known to be produced in the gut, and kidney function impairment of patients in various stages of CKD.
2. Materials and Methods 2.1. Study Participants This study was approved by the Medical Ethics Committee of Seoul National University Hospital (IRB number: 1808-153-967) and complied with the Declaration of Helsinki. The study included only participants who provided informed consent and agreed to submit fecal specimens to the human stool repository (IRB number: 1802-062-921). The human stool repository includes samples collected from patients who underwent percutaneous kidney biopsy (IRB number: 1508-046-694) and from kidney transplant donors and recipients prior to transplantation (IRB number: 1703-062-839). We also obtained serum samples at the same time as feces collection from registered participants. The serum samples for this study were provided by the Seoul National University Hospital Human Biobank, a member of the National Biobank of Korea. All samples derived from the National Biobank of Korea were obtained with informed consent under institutional review-board-approved protocols. 2.2. Clinical Information on Study Participants We collected demographic information from the study participants, including details on age, sex, height, and weight, and established whether any of these individuals had comorbidities, including hypertension and diabetes mellitus, as determined from their clinical and medication histories compiled in an electronic medical record system. Participants who had used antibiotics within one month were excluded. Furthermore, given that gut microbiota may be affected by the systemic inflammation status, we also collected data on the serum levels of highly sensitive C-reactive protein, which is used as a marker of inflammation. Kidney function was assessed by measuring serum blood urea nitrogen and creatinine levels at the time of fecal sample collection. The presence of hematuria and amount of proteinuria were quantified based on urine microscopic examination and random urine protein to creatinine ratio determinations, respectively. Plasma hemoglobin and serum albumin concentrations were examined to evaluate the anemia and nutritional status, respectively. The etiologies of CKD were established by pathological confirmation, imaging studies, or clinical diagnosis. The estimated glomerular filtration rate was calculated using the Chronic Kidney Disease Epidemiology Collaboration calculation formula according to CKD staging. We divided the study participants into the following groups: kidney donors without evidence of kidney disease as healthy controls; patients with stage 1 and 2 CKD as mild CKD; patients with stage 3 and 4 CKD, and stage 5 CKD without dialysis, as moderate to severe CKD; and patients with stage 5 CKD requiring dialysis as end-stage renal disease (ESRD). 2.3. Stool DNA Extraction and MiSeq Sequencing
The stool samples collected from all study participants were immediately stored in a deep freezer at −80 °C. Stool DNA extraction was performed using a QIAamp® Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions [21]. The extracted DNA was used to amplify the V4–5 variable regions of the 16S rRNA gene. Amplification was performed in accordance with the MiSeq system protocol for preparing a 16S metagenomics sequencing library (Illumina, Inc., San Diego, CA, USA). The amplicons of each sample were purified using Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA), and the purified amplicons were quantified using a PicoGreen dsDNA Assay kit (Invitrogen, Carlsbad, CA, USA). Equimolar concentrations of each library were pooled and sequenced using the Illumina MiSeq system (250-base pair paired ends) according to the manufacturer’s instructions.
2.4. Sequence Data Analysis
For microbiota analysis, the obtained sequence reads were analyzed using the Microbial Genomics Module of CLC genomic workbench v. 11.0.1 (Qiagen, Aarhus, Denmark). Briefly, raw sequences were merged, and sequences with short read lengths (<400 base pairs of merged reads) or low-quality score and chimeric reads were removed using the USEARCH pipeline v. 11.0.667 (http://www.drive5.com/usearch). Primer sequences were removed from the merged sequences, and filtered sequences were subsequently clustered into operational taxonomic units (OTUs) based on 97% sequence identity. The taxonomic positions of representative sequences in each OTU were assigned based on comparisons with the EzTaxon-e reference database [22]. To compare diversity indices among samples, sequence read numbers were normalized by random subsampling and the indices were calculated using Mothur [23]. Principal coordinate analysis (PCoA) based on Bray–Curtis distances was performed using Calypso to compare microbiota compositions among samples [24]. The functional roles of microbiota were predicted using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) [25], and the co-occurrence networks for microbiota in each group were inferred based on Spearman correlation matrices and selected according to q < 0.05 (Benjamini–Hochberg-corrected). Networks were constructed for significantly different direct interactions among genera in the different patient groups, with the visualization of networks and calculations performed using CoNet with Cytoscape (v. 3.4.0).
2.5. Serum Metabolite Analysis
Serum specimens collected from all study participants were stored in a −180 °C nitrogen tank. To measure metabolites that are representative of uremic toxins, we selected as target metabolites p-cresyl sulfate, p-cresyl glucuronide, indoxyl sulfate, and TMAO based on the findings of previous studies [11,26,27]. The serum concentrations of these four metabolites were determined by liquid chromatography–tandem mass spectrometry, as described previously [28,29]. Twenty-microliter samples were placed in microtubes, followed by the addition of 20 μL of acetonitrile containing an internal standard and 500 μL of acetonitrile containing 0.1% formic acid. These sample mixtures were vortexed for 30 s, followed by centrifugation at 13,000 rpm for 5 min. The resulting supernatants were transferred into injection vials and subjected to liquid chromatography–tandem mass spectrometry. We used an Agilent 1260 Infinity high-performance liquid chromatography system in conjunction with an API 4000 QTRAP mass spectrometry system (Agilent Technologies, Santa Clara, CA, USA). From the chromatograms, we calculated the area ratios of p-cresyl sulfate, p-cresyl glucuronide, indoxyl sulfate, TMAO, and the internal standard and determined the concentrations of the metabolites using previously prepared calibration curves. For p-cresyl glucuronide, the minimum detection limit was 5 ng/mL, and values were designated as zero when measurements were below this value.
2.6. Statistical Analysis For the baseline characteristics, continuous variables are expressed as the means and standard deviations and categorical variables are expressed as percentages. Differences between samples were evaluated using the Mann–Whitney U and Kruskal–Wallis tests. The p-value for trends was calculated using the Stata module “nptrend,” which is an extension of the Mann–Whitney U test that can be used to perform nonparametric tests for trends across ordered groups; p < 0.05 was considered to indicate statistical significance. Correlations between microbiota and uremic toxins were determined by linear regression. The Benjamini–Hochberg false discovery rate (FDR) was applied to correct for multiple testing and FDR-adjusted p-values of less than 0.05 were considered as significant. Statistical analyses were performed using R v. 3.5.0 (R Core Team), Stata v. 15.1 (StataCorp, College Station, TX, USA), and GraphPad Prism v. 8.1.1 (GraphPad, Inc., San Diego, CA, USA). 3. Results 3.1. Comparisons of Baseline Characteristics and Serum Uremic Metabolites according to CKD Group
We analyzed samples collected from a total of 149 participants, among whom there were 46 controls and 103 patients with CKD. The patients with CKD comprised 36 subjects with mild CKD, 32 with moderate to severe CKD, and 35 with dialysis-dependent ESRD. The baseline characteristics, laboratory data, and uremic metabolite concentrations were compared among groups (Table 1). We found no significant differences among groups with respect to the distributions of age and sex. Patients with ESRD were, however, significantly leaner than patients in the control and other CKD groups (p = 0.003). The participants with lower kidney function had more comorbidities, including hypertension and diabetes mellitus (p < 0.001 and 0.020, respectively). In line with expectations, we observed that plasma hemoglobin and serum albumin levels were negatively correlated with CKD severity (p < 0.001, in both cases). The proportion of patients with anemia increased significantly with worsening CKD severity (p < 0.001). However, serum levels of highly sensitive C-reactive protein did not differ significantly among the four groups.
Measurements of the four uremic metabolites p-cresyl sulfate, p-cresyl glucuronide, indoxyl sulfate, and TMAO in the sera revealed clear positive correlations with CKD severity (Kruskal–Wallis, p < 0.001), although there were no significant differences in the levels of these metabolites between participants in the control and mild CKD groups (Supplementary Figure S1).
3.2. Differences in Microbiota Composition according to CKD Group
A total of 1,796 OTUs were detected in all patients. To compare the diversity indices among samples, read numbers were normalized to 6,400 by random subsampling, and we accordingly detected no significant differences in bacterial diversity among the different patient groups (p > 0.05; Figure 1a). However, we detected a larger number of observed OTUs in the ESRD group compared to the control and mild CKD groups, respectively (p < 0.05; Figure 1b). PCoA based on Bray–Curtis distances did not clearly distinguish differences in the microbiota in different CKD groups (Figure 1c). There were also no differences among the groups with respect to the proportions of the different phyla, although Bacteroidetes, Firmicutes, and Proteobacteria were dominant phyla in the gut microbiota of all groups (Figure 1d). A comparison of the microbiota among groups at the genus level (Figure 1e) revealed 18 main genera, defined as those constituting more than 1% of the total microbiota in samples. In most samples, Bacteroides and Prevotella were the dominant genera (with averages of 28.5% and 15.0% in samples, respectively). When we evaluated the differences in genera among groups based on multiple group comparisons, the proportions of Alistipes, Oscillibacter, Lachnospira, Veillonella, and Dialister were shown to be significantly different among the four groups (p < 0.05; Figure 2). The proportions of these genera in the microbiota of the control group were significantly different from those in the moderate to severe CKD group, but not from those in the mild CKD group. Furthermore, the proportions of Alistipes, Oscillibacter, Lachnospira, and Veillonella in the mild CKD group differed from those in the moderate to severe CKD group. In contrast, we detected no significant differences in the proportions of Alistipes, Oscillibacter, Lachnospira, Veillonella, and Dialister in the moderate to severe CKD and ESRD groups. The relative abundance of Alistipes and Oscillibacter was increased with the progression in CKD severity (tests for trends, p = 0.001 and 0.016, respectively), whereas the abundance of Lachnospira, Veillonella, and Dialister decreased with increasing CKD severity (tests for trends, p = 0.019, 0.012, and p < 0.001, respectively).
To determine whether mutual interactions among the aforementioned five genera and other gut microbes differed according to CKD severity, we performed co-occurrence network analysis. We accordingly found that as kidney function decreased, the network appeared to become more active and complex in terms of both positive and negative interactions of gut microbiota in these five genera with other microbes (Figure 3). Oscillibacter and Veillonella were found to show the highest values of betweenness centrality in networks and appeared to act as hubs in the microbial networks constructed for patients with moderate or higher stage CKD. Oscillibacter showed positive correlations with other genera, whereas Veillonella species were negatively correlated. These results indicate that the relative abundance of Oscillibacter associated with kidney function impairment gives rise to the co-occurrence of other gut microbiota.
3.3. Microbiota-Related Uremic Toxins
We subsequently analyzed differences in the relative abundance of gut microbiota according to the levels of uremic toxins in the sera of patients and controls. The genera found to be associated with the levels of uremic toxins are shown on a logarithmic scale in Table 2. According to the FDR-adjusted p-values (q-values) in multivariable linear regression, p-cresyl sulfate showed a significant association with six major genera, with positive correlations observed with Alistipes, Oscillibacter, and Subdoligranulum (q < 0.001, q < 0.001, and q = 0.023, respectively) and negative correlations with Lachnospira, Veillonella, and Megamonas (q = 0.039, q = 0.014, and q = 0.034, respectively). Serum p-cresyl glucuronide and indoxyl sulfate levels were positively correlated with Alistipes (q = 0.010 and q = 0.035, respectively) and Oscillibacter (q = 0.001 and q = 0.037, respectively). Serum TMAO levels were associated only with Oscillibacter (q = 0.006). Unexpectedly, we found Oscillibacter to be associated with all four measured uremic toxins, although each of these metabolites is derived from different parent compounds.
3.4. Predicted Functional Analysis of Gut Microbiota among the CKD Groups
The functions of the gut microbiota were predicted based on PICRUSt analysis and a comparison of pathways based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) categories among the different CKD groups. At KEGG Ortholog (KO) level 3, four pathways were predicted to be significantly different among groups (q < 0.05; Supplementary Table S1). Among these significant pathways, KEGG terms relating to “Pyruvate metabolism” and “Methane metabolism” were predicted to increase with decreasing kidney function, whereas “Riboflavin metabolism” was predicted to decrease with kidney function deterioration (Figure 4a).
Given that pyruvate metabolism represents a key intersection in the network of various metabolic pathways, we focused on this pathway and examined the contribution of major genera to pyruvate metabolism (Figure 4b). We accordingly found that only the proportional contributions of Oscillibacter and Veillonella showed significant differences among the patient groups (p = 0.0011 and p = 0.0014, respectively; Figure 4c), and as CKD severity increased, the contribution of Oscillibacter to pyruvate metabolism increased, whereas that of Veillonella decreased.
As Oscillibacter showed some evidence of an association with uremic metabolites and CKD severity as well as differences among CKD groups regarding the proportional contribution to pyruvate metabolism, we performed a detailed analysis of the contribution of Oscillibacter to each orthologous gene involved in pyruvate metabolism. Among all orthologs, we found that functional orthologs of the E1 component of pyruvate dehydrogenase (K00161, K00162) showed the highest proportional contribution of Oscillibacter (33.0% for K00161 and 33.1% for K00162). Among all genera, Oscillibacter also showed a high proportional contribution (21.3%) to the pyruvate dehydrogenase E2 component (K00627). These findings indicate that Oscillibacter is a primary source of pyruvate dehydrogenase, a key enzyme involved in glucose oxidation that converts pyruvate to acetyl-CoA. The proportional contributions of Oscillibacter to each orthologous gene involved in pyruvate metabolism compared with those of other genera are shown in Supplementary Table S2.
Although comparisons among the different patient groups revealed no significant differences in the contributions of Oscillibacter to each of pyruvate metabolism orthologs according to CKD severity, we found that the contribution to pyruvate dehydrogenase E1 component decreased by approximately 10% in CKD groups compared with the control group (for K00161: 40.1% in the control, 27.0% in mild CKD, 32.2% in moderate to severe CKD, and 31.4% in ESRD; for K00162: 40.2% in the control, 27.0% in mild CKD, 32.2% in moderate to severe CKD, and 31.4% in ESRD). In contrast, we found that the contribution of Oscillibacter to lactate dehydrogenase, which converts pyruvate to lactate, was significantly increased (Kruskal–Wallis p = 0.003; q-value between control and ESRD = 0.009; Figure 5). In addition, we observed significant increases in the contributions to other collateral metabolic pathways, among which was an increase in oxaloacetate and formate production at the expense of a direct generation of acetyl-CoA via pyruvate dehydrogenase (K01958: pyruvate carboxylase, p < 0.001; K00656: formate C-acetyltransferase, p = 0.009; Supplementary Table S2 and Supplementary Figure S2).
4. Discussion
The gut microbiota as a source of uremic toxin accumulation in patients with CKD has gained attention. CKD affects bacterial fermentation processes, including colonic transit time and colon microbiota composition, resulting in changes in the microenvironment of the colon [2,6]. However, the integrating effect of CKD on the gut microbiota and their associated metabolites have not been characterized in detail. In this study, we investigated the associations between specific gut microbes, including Oscillibacter, and levels of uremic metabolites, along with changes in predicted metabolic pathways related to kidney dysfunction in patients with CKD.
Although several studies have analyzed the gut microbiota in patients with CKD, most studies focused on the differences between patients with ESRD and healthy controls. For example, Vaziri et al. identified a significant elevation in the relative abundance of 190 OTUs in patients with ESRD compared with healthy controls [30]. Similarly, using a bacterial culture method based on diluted stool samples, Hida et al. observed elevated contents of several gut bacteria in patients on hemodialysis compared to healthy controls [31]. In contrast to these previous studies, we focused on patients with pre-dialytic CKD, who were subclassified based on their estimated glomerular filtration rate to examine serial changes in gut microbiota according to kidney functional impairment. We accordingly identified specific gut microbes (Alistipes, Oscillibacter, Lachnospira, Veillonella, and Dialister), the relative abundances of which showed continuous change concomitant with changes in kidney function. Given that previous studies only investigated differences in OTUs at the family level or differences in OTUs between control and subjected with ESRD [30,31], it is difficult to properly compare the changes in OTUs discovered in previous studies and those in the current study. However, in the study of Varizi et al. [30], the class Clostridia which includes Oscillibacter, was reported to show a significant elevation in patients with ESRD compared with control subjects, thereby indicating that our findings are at least partially consistent with those reported previously.
Of the bacterial taxa showing significant differences among the CKD groups, Oscillibacter and Veillonella appear to have important roles as network hubs for the other microbes in the analysis of network co-occurrence. Particularly, Oscillibacter showed significant positive associations with diverse genera in patients with advanced CKD, indicating that the bacterial species in this genus interact with those in other bacteria in the gut microbiota of patients with advanced CKD. Moreover, we found that Oscillibacter showed common associations with the four uremic metabolites p-cresyl sulfate, p-cresyl glucuronate, indoxyl sulfate, and TMAO, which are derived from a diverse range of parent compounds, including phenols, indoles, and quaternary amines. These broad effects on uremic metabolites and other microbes suggest that Oscillibacter play a pivotal role in creating a favorable inflammatory environment that facilitates the proliferation and activation of a number of uremic toxin-producing pathogens. Although they are poorly represented in culture collections, the genus Oscillibacter was detected in human gut microbiota related to some pathologic state [32,33,34]. Elevated Oscillibacter abundance had been found in patients with stroke and closely related to gut permeability and host inflammation [34,35]. However, the physiological role of Oscillibacter with respect to kidney disease has not been reported. Veillonella, another core microbiota associated with CKD stages, has not been studied in gut microbiota from patients with kidney disease, although they were previously reported to be decreased in saliva and oral swab samples of patients with CKD [36].
In the present study, we conducted predicted functional analyses to detect clues indicating the physiological role of the microbiota in the production of toxins associated with kidney dysfunction. We found that as kidney function deteriorated, there was an increase in the microbial contributions to pyruvate metabolism, particularly with regards to the proportional contribution of Oscillibacter-related genes. Moreover, we observed that the elevated abundance of Oscillibacter and their contribution to pyruvate metabolism were weighted toward anaerobic glycolysis (represented by lactate production) rather than to aerobic glycolysis as kidney function deteriorated. A recent metabolomics study revealed that citric acid metabolism in the tricarboxylic acid cycle, a subsequent pathway of pyruvate, is the most altered metabolic pathway in patients with nondiabetic CKD stages 3–4 compared to in healthy controls [37]. Furthermore, the genus Oscillibacter is significantly associated with serum lactate levels according to an animal study and the genus Veillonella is related to the metabolism of lactate to propionate in athletes [38,39]. The changes in lactate metabolism mediated by these microbes promote changes in intestinal pH, a decrease in which has previously been shown to disrupt the growth of certain members of the gut microbiota, thereby modifying microbial and metabolic interactions [40]. Based on these processes, changes in the gut microbiota may influence the production and absorption of uremic metabolites by gut intraluminal environmental disturbance, and vice versa. Further comprehensive studies of the association between changes in local and systemic metabolic pathways and specific microbes in patients with renal disease are needed.
Although the present study has notable strengths, such as including an entire spectrum of CKD stages, there were also limitations. First, we did not establish any clear causal relationships among microbiota, microbial metabolites, and CKD. Although we demonstrated associations between microbes and specific metabolites within the microbiota, the causal mechanisms should be determined in further studies. Second, the number of participants registered in each group was relatively small. However, compared with previous studies of CKD and microbiomes, we evaluated a similar or larger number of participants. Furthermore, although we obtained comparative values at certain times during CKD progression, we did not observe serial changes in renal function over time. Finally, we did not assess interactions between the human host and observed microbiota as immunological factors or perform gene expression analyses. However, this study provides a possible association between gut microbiota and renal functions. 5. Conclusions In this study, we detected significant correlations among the gut microbiota, uremic metabolites, and renal functions in patients with CKD, providing insights into the role of the gut microbiome in the progression of kidney disease. Our findings indicate that the deterioration in renal function observed in patients with CKD is related to increases in the relative abundance of Oscillibacter within the gut microbiota, the bacteria of which interact with other gut microbes. Based on our observations, Oscillibacter may create a favorable environment for the production of various uremic metabolites, and such action may be related to the altered contribution of Oscillibacter to the pyruvate metabolism pathway. Further studies are needed to confirm our results and perform more in-depth analyses of the various interrelationships among the gut microbiome, microbial metabolites, and CKD progression.
Enlarge this image.Figure 1. Comparison of diversity and taxonomy of gut microbiota according to chronic kidney disease (CKD) severity. (a) Comparison of Shannon diversity indices of gut microbiota among CKD groups. (b) Comparison of number of observed operational taxonomic units (OTUs) of gut microbiota among CKD groups (*, p < 0.05; **, p < 0.01). (c) Principal coordinates analysis (PCoA) based on Bray-Curtis distances of gut microbiota among healthy controls (red), patients with mild CKD (blue), patients with moderate CKD (gray), and patients with end-stage renal disease (ESRD, yellow). The first two axes of the PCoA plot are represented by principal coordinate axis 1 (PCoA1) and principal coordinate axis 2 (PCoA2). (d) Comparison of microbiota composition among CKD groups at the phylum level. red color, Bacteroidetes; blue color, Firmicutes; yellow color, Proteobacteria; gray color, and other phyla with mean relative abundances <1% of total abundance in samples. (e) Comparison of microbiota composition among CKD groups at the genus level. The heatmap plot shows mainly detected genera with mean relative abundances >1% of total abundance in samples; values <1%, unclassified or unidentified, are classified as 'others.' The abundances of all genera were plotted after conversion to a binary logarithmic scale.
Enlarge this image.Figure 2. Significantly different genera according to renal function in chronic kidney disease (CKD) patients. We determined bacterial genera showing significant differences among patient groups. (a) Genera showing significant increasing trends according to CKD group. (b) Genera showing significant decreasing trends according to CKD group. q-values were determined using the Benjamini-Hochberg method based on the p-values obtained using Mann-Whitney U tests. Single and double asterisks for group comparisons indicate q < 0.05 and q < 0.005, respectively. ns, not significant.
Enlarge this image.Figure 3. Co-occurrence network for significantly different genera associated with different chronic kidney disease (CKD) groups. The network shows only significant correlations with q < 0.05. Green lines denote a positive correlation between microbiota, whereas red lines denote a negative correlation. The size of circles indicates the relative abundance of each genus.
Enlarge this image.Figure 4. Predicted functional analysis. (a) Comparison of significantly different metabolic pathways among chronic kidney disease (CKD) groups based on predicted functional analysis. Only pathways showing a significant difference among groups were plotted (see Supplementary Table S1 for a complete list of functional pathways). Single and double asterisks for group comparisons indicate q < 0.05 and q < 0.005, respectively. (b) Comparison of the total contribution of genera to the pyruvate metabolism pathway among groups. Only major genera showing a mean proportional contribution >1% among all samples are listed. (c) Comparison of the proportional contributions of Oscillibacter (upper) and Veillonella to the pyruvate metabolism pathway among patient groups. q-values were determined using the Benjamini-Hochberg method based on the p-values obtained using Mann-Whitney U tests. *, q < 0.05; **, q < 0.005.
Enlarge this image.Figure 5. Schematic diagram showing the contribution of Oscillibacter to pyruvate metabolism according to renal function. Flow charts were derived based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Small bar charts represent a comparison of the contribution of Oscillibacter to different orthologous genes according to the patient group. q-values were determined using the Benjamini-Hochberg method based on the p-values obtained using Mann-Whitney U tests. Asterisks in the bar charts indicate significance (q-value) compared with the control group. ns, not significant; *, q < 0.05.
| Variables | Control | Mild CKD | Moderate to Severe CKD | ESRD | |
|---|---|---|---|---|---|
| Total N = 149 | N = 46 | N = 36 | N = 32 | N = 35 | p |
| Clinical parameters | |||||
| Age (y) | 47.0 ± 10.8 | 49.8 ± 15.1 | 52.4 ± 11.9 | 48.9 ± 12.2 | 0.251 |
| Male sex (%) | 16 (34.8) | 21 (58.3) | 17 (53.1) | 15 (42.9) | 0.15 |
| Body mass index (%) | 23.3 ± 3.0 | 24.6 ± 3.5 | 23.6 ± 3.2 | 21.8 ± 4.3 | 0.003 |
| Diabetes mellitus (%) | 0 (0) | 3 (8.3) | 6 (18.8) | 6 (17.1) | 0.02 |
| Hypertension (%) | 3 (6.5) | 21 (58.3) | 22 (68.8) | 26 (74.3) | <0.001 |
| Blood urea nitrogen (mg/dL) | 12.1 ± 2.9 | 14.4 ± 3.9 | 43.7 ± 25.7 | 45.8 ± 15.8 | <0.001 |
| Serum creatinine (mg/dL) | 0.7 ± 0.2 | 0.8 ± 0.2 | 3.8 ± 2.5 | 7.8 ± 2.6 | <0.001 |
| CKD-EPI eGFR (mL/min/1.73 m2) | 101.6 ± 19.0 | 98.3 ± 26.1 | 25.2 ± 17.2 | 7.2 ± 2.5 | <0.001 |
| Urine RBC (number /HPF) | <0.001 | ||||
| 0 | 24 (52.2) | 4 (11.1) | 7 (21.9) | NA | |
| 1–4 | 17 (37.0) | 7 (19.4) | 9 (28.1) | NA | |
| 5≤ | 5 (10.9) | 25 (69.4) | 16 (50.0) | NA | |
| Urine protein/creatinine ratio | 0.05 ± 0.03 | 3.6 ± 3.4 | 3.3 ± 3.4 | NA | <0.001 |
| Plasma hemoglobin (g/dL) | 13.8 ± 1.3 | 12.9 ± 1.6 | 11.1 ± 2.0 | 10.3 ± 1.6 | <0.001 |
| Anemia (%) | 3 (6.5) | 13 (36.1) | 24 (75.0) | 30 (85.7) | <0.001 |
| Serum albumin (mg/dL) | 4.4 ± 0.3 | 3.5 ± 0.7 | 3.8 ± 0.5 | 3.8 ± 0.4 | <0.001 |
| Serum C-reactive protein (mg/dL) | 0.1 ± 0.4 | 0.2 ± 0.3 | 0.7 ± 1.3 | 0.3 ± 0.8 | 0.321 |
| Etiology of CKD (biopsy proven/clinical diagnosis) | <0.001 | ||||
| Diabetes mellitus | NA | 0 | 4 (0/4) | 6 (0/6) | |
| Hypertension | NA | 0 | 1 (1/0) | 1 (0/1) | |
| Glomerulonephritis | NA | 35 (35/0) | 21 (18/3) | 14 (4/10) | |
| Polycystic kidney | NA | 0 | 4 (0/4) | 3 (0/3) | |
| Others | NA | 1 | 2 | 11 | |
| Serum uremic metabolites | |||||
| P-cresyl sulfate (ug/mL) | 9.5 ± 10.8 | 7.00 ± 8.7 | 63.2 ± 56.0 | 111.6 ± 87.0 | <0.001 |
| P-cresyl glucuronide * (ng/mL) | 18.2 ± 18.0 | 19.8 ± 19.3 | 114.5 ± 110.1 | 746.7 ± 880.5 | <0.001 |
| Indoxyl sulfate (ug/mL) | 0.7 ± 0.4 | 0.7 ± 0.6 | 7.3 ± 7.6 | 26.0 ± 17.8 | <0.001 |
| TMAO (ug/mL) | 0.6 ± 1.1 | 0.8 ± 1.2 | 4.9 ± 5.9 | 13.9 ± 17.4 | <0.001 |
* Calculated excluding samples measured below the minimum measurement limit of p-cresyl glucuronide (5 ng/mL). Abbreviations: CKD, chronic kidney disease; eGFR, estimated glomerular filtration rate; RBC, red blood cell; HPF, high-power field; TMAO, trimethylamine N-oxide; NA, not available.
| Predictors | Regression Coefficient | Standard Error | Adjusted R2 | p | FDR |
|---|---|---|---|---|---|
| p-cresyl sulfate (log) | |||||
| Alistipes | 0.207 | 0.050 | 0.100 | <0.001 | <0.001 |
| Oscillibacter | 0.238 | 0.045 | 0.155 | <0.001 | <0.001 |
| Lachnospira | –0.105 | 0.043 | 0.033 | 0.016 | 0.039 |
| Veillonella | –0.079 | 0.027 | 0.050 | 0.004 | 0.014 |
| Subdoligranulum | 0.234 | 0.085 | 0.042 | 0.007 | 0.023 |
| Megamonas | –0.092 | 0.036 | 0.036 | 0.012 | 0.034 |
| p-cresyl glucuronate (log) * | |||||
| Prevotella | –0.023 | 0.009 | 0.034 | 0.014 | 0.062 |
| Alistipes | 0.189 | 0.057 | 0.062 | 0.001 | 0.010 |
| Oscillibacter | 0.213 | 0.053 | 0.094 | < 0.001 | 0.001 |
| Lachnospira | –0.104 | 0.049 | 0.024 | 0.034 | 0.103 |
| Subdoligranulum | 0.224 | 0.097 | 0.028 | 0.023 | 0.081 |
| Indoxyl sulfate (log) | |||||
| Alistipes | 0.126 | 0.043 | 0.048 | 0.004 | 0.035 |
| Oscillibacter | 0.112 | 0.041 | 0.043 | 0.007 | 0.037 |
| Lachnospira | –0.091 | 0.036 | 0.034 | 0.014 | 0.058 |
| Subdoligranulum | 0.177 | 0.073 | 0.032 | 0.016 | 0.058 |
| TMAO (log) | |||||
| Prevotella | –0.016 | 0.007 | 0.027 | 0.026 | 0.080 |
| Alistipes | 0.104 | 0.046 | 0.027 | 0.026 | 0.080 |
| Oscillibacter | 0.145 | 0.042 | 0.067 | 0.001 | 0.006 |
| Lachnospira | –0.087 | 0.039 | 0.026 | 0.027 | 0.080 |
| Dialister | –0.107 | 0.043 | 0.034 | 0.014 | 0.080 |
Only variables with p-values under 0.05 are listed in this table while FDR was calculated with all major genera. Listed variables are in descending order according to relative abundance. * Samples measured below the minimum measurement limit (5 ng/mL) are calculated as imputation to zero. Abbreviations: FDR, false discovery rate; TMAO, trimethylamine N-oxide.
Supplementary Materials
The following are available online at https://www.mdpi.com/2076-2607/8/6/907/s1, Figure S1. Serum levels of four uremic metabolites according to CKD groups; Figure S2. Hypothetical pathway diagram using the KEGG module associated with changes in Oscillibacter contribution according to renal function; Table S1. Differentially represented KEGG categories predicted by PICRUSt (level 3) in 4 groups; Table S2. Contribution of Oscillibacter to functional orthologs associated with the pyruvate metabolism pathway.
Author Contributions
B.-S.K. and H.L. contributed to the study conception. J.E.K., J.I.P., H.C., and H.L. contributed to the study design. J.E.K., J.I.P., H.J.C., S.H.Y., J.P.L., D.K.K., K.W.J., and Y.S.K. contributed to data acquisition. J.E.K., H.-E.K., M.-J.K., B.-Y.K. and B.-S.K. contributed to data analysis. J.E.K., H.-E.K., J.I.P., H.J.C., B.-S.K., and H.L. contributed to drafting of the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding
This research was supported by the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (grant number: 2019R1A2C2011465, 2015M3C9A2054342).
Acknowledgments
The biospecimens for this study were provided by the Seoul National University Hospital Human Biobank, a member of the National Biobank of Korea, which is supported by the Ministry of Health and Welfare.
Conflicts of Interest
The authors declare no conflicts of interest.
1. Glendinning, L.; Free, A. Supra-organismal interactions in the human intestine. Front. Cell. Infect. Microbiol. 2014, 4, 47.
2. Mahmoodpoor, F.; Rahbar Saadat, Y.; Barzegari, A.; Ardalan, M.; Zununi Vahed, S. The impact of gut microbiota on kidney function and pathogenesis. Biomed. Pharmacother. 2017, 93, 412-419.
3. Vijay-Kumar, M.; Aitken, J.D.; Carvalho, F.A.; Cullender, T.C.; Mwangi, S.; Srinivasan, S.; Sitaraman, S.V.; Knight, R.; Ley, R.E.; Gewirtz, A.T. Metabolic syndrome and altered gut microbiota in mice lacking Toll-like receptor 5. Science 2010, 328, 228-231.
4. Turnbaugh, P.J.; Ley, R.E.; Mahowald, M.A.; Magrini, V.; Mardis, E.R.; Gordon, J.I. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 2006, 444, 1027-1031.
5. Round, J.L.; Mazmanian, S.K. The gut microbiota shapes intestinal immune responses during health and disease. Nat. Rev. Immunol. 2009, 9, 313-323.
6. Anders, H.J.; Andersen, K.; Stecher, B. The intestinal microbiota, a leaky gut, and abnormal immunity in kidney disease. Kidney Int. 2013, 83, 1010-1016.
7. de Almeida Duarte, J.B.; de Aguilar-Nascimento, J.E.; Nascimento, M.; Nochi, R.J., Jr. Bacterial translocation in experimental uremia. Urol. Res. 2004, 32, 266-270.
8. Ramezani, A.; Raj, D.S. The gut microbiome, kidney disease, and targeted interventions. J. Am. Soc. Nephrol. 2014, 25, 657-670.
9. Smith, E.A.; Macfarlane, G.T. Enumeration of human colonic bacteria producing phenolic and indolic compounds: Effects of pH, carbohydrate availability and retention time on dissimilatory aromatic amino acid metabolism. J. Appl. Bacteriol. 1996, 81, 288-302.
10. Wong, J.; Piceno, Y.M.; DeSantis, T.Z.; Pahl, M.; Andersen, G.L.; Vaziri, N.D. Expansion of urease- and uricase-containing, indole- and p-cresol-forming and contraction of short-chain fatty acid-producing intestinal microbiota in ESRD. Am. J. Nephrol. 2014, 39, 230-237.
11. Mair, R.D.; Sirich, T.L.; Meyer, T.W. Uremic Toxin Clearance and Cardiovascular Toxicities. Toxins (Basel) 2018, 10, 226.
12. Wu, I.W.; Hsu, K.H.; Lee, C.C.; Sun, C.Y.; Hsu, H.J.; Tsai, C.J.; Tzen, C.Y.; Wang, Y.C.; Lin, C.Y.; Wu, M.S. p-Cresyl sulphate and indoxyl sulphate predict progression of chronic kidney disease. Nephrol. Dial. Transplant. 2011, 26, 938-947.
13. Vanholder, R.; Meert, N.; Schepers, E.; Glorieux, G.; Argiles, A.; Brunet, P.; Cohen, G.; Drueke, T.; Mischak, H.; Spasovski, G.; et al. Review on uraemic solutes II-Variability in reported concentrations: Causes and consequences. Nephrol. Dial. Transplant. 2007, 22, 3115-3121.
14. Martinez, A.W.; Recht, N.S.; Hostetter, T.H.; Meyer, T.W. Removal of P-cresol sulfate by hemodialysis. J. Am. Soc. Nephrol. 2005, 16, 3430-3436.
15. Melamed, M.L.; Plantinga, L.; Shafi, T.; Parekh, R.; Meyer, T.W.; Hostetter, T.H.; Coresh, J.; Powe, N.R. Retained organic solutes, patient characteristics and all-cause and cardiovascular mortality in hemodialysis: Results from the retained organic solutes and clinical outcomes (ROSCO) investigators. BMC Nephrol. 2013, 14, 134.
16. Niwa, T. Role of indoxyl sulfate in the progression of chronic kidney disease and cardiovascular disease: Experimental and clinical effects of oral sorbent AST-120. Ther. Apher. Dial. 2011, 15, 120-124.
17. Velasquez, M.T.; Ramezani, A.; Manal, A.; Raj, D.S. Trimethylamine N-Oxide: The Good, the Bad and the Unknown. Toxins (Basel) 2016, 8, 326.
18. Wright, A.T. Gut commensals make choline too. Nat. Microbiol. 2019, 4, 4-5.
19. Barrows, I.R.; Ramezani, A.; Raj, D.S. Gut Feeling in AKI: The Long Arm of Short-Chain Fatty Acids. J. Am. Soc. Nephrol. 2015, 26, 1755-1757.
20. Meijers, B.K.; Claes, K.; Bammens, B.; de Loor, H.; Viaene, L.; Verbeke, K.; Kuypers, D.; Vanrenterghem, Y.; Evenepoel, P. p-Cresol and cardiovascular risk in mild-to-moderate kidney disease. Clin. J. Am. Soc. Nephrol. 2010, 5, 1182-1189.
21. Li, M.; Gong, J.; Cottrill, M.; Yu, H.; de Lange, C.; Burton, J.; Topp, E. Evaluation of QIAamp DNA Stool Mini Kit for ecological studies of gut microbiota. J. Microbiol. Methods 2003, 54, 13-20.
22. Yoon, S.H.; Ha, S.M.; Kwon, S.; Lim, J.; Kim, Y.; Seo, H.; Chun, J. Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 2017, 67, 1613-1617.
23. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 2009, 75, 7537-7541.
24. Zakrzewski, M.; Proietti, C.; Ellis, J.J.; Hasan, S.; Brion, M.J.; Berger, B.; Krause, L. Calypso: A user-friendly web-server for mining and visualizing microbiome-environment interactions. Bioinformatics 2017, 33, 782-783.
25. Langille, M.G.; Zaneveld, J.; Caporaso, J.G.; McDonald, D.; Knights, D.; Reyes, J.A.; Clemente, J.C.; Burkepile, D.E.; Vega Thurber, R.L.; Knight, R.; et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 2013, 31, 814-821.
26. Lin, C.J.; Chen, H.H.; Pan, C.F.; Chuang, C.K.; Wang, T.J.; Sun, F.J.; Wu, C.J. p-Cresylsulfate and indoxyl sulfate level at different stages of chronic kidney disease. J. Clin. Lab. Anal. 2011, 25, 191-197.
27. Liabeuf, S.; Glorieux, G.; Lenglet, A.; Diouf, M.; Schepers, E.; Desjardins, L.; Choukroun, G.; Vanholder, R.; Massy, Z.A.; European Uremic Toxin Work, G. Does p-cresylglucuronide have the same impact on mortality as other protein-bound uremic toxins? PLoS ONE 2013, 8, e67168.
28. Kanemitsu, Y.; Asaji, K.; Matsumoto, Y.; Tsukamoto, H.; Saigusa, D.; Mukawa, C.; Tachikawa, T.; Abe, T.; Tomioka, Y. Simultaneous quantitative analysis of uremic toxins by LC-MS/MS with a reversed-phase/cation-exchange/anion-exchange tri-modal mixed-mode column. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2017, 1068-1069, 1-8.
29. Shu, C.; Chen, X.; Xia, T.; Zhang, F.; Gao, S.; Chen, W. LC-MS/MS method for simultaneous determination of serum p-cresyl sulfate and indoxyl sulfate in patients undergoing peritoneal dialysis. Biomed. Chromatogr. 2016, 30, 1782-1788.
30. Vaziri, N.D.; Wong, J.; Pahl, M.; Piceno, Y.M.; Yuan, J.; DeSantis, T.Z.; Ni, Z.; Nguyen, T.H.; Andersen, G.L. Chronic kidney disease alters intestinal microbial flora. Kidney Int. 2013, 83, 308-315.
31. Hida, M.; Aiba, Y.; Sawamura, S.; Suzuki, N.; Satoh, T.; Koga, Y. Inhibition of the accumulation of uremic toxins in the blood and their precursors in the feces after oral administration of Lebenin, a lactic acid bacteria preparation, to uremic patients undergoing hemodialysis. Nephron 1996, 74, 349-355.
32. Mondot, S.; Kang, S.; Furet, J.P.; Aguirre de Carcer, D.; McSweeney, C.; Morrison, M.; Marteau, P.; Dore, J.; Leclerc, M. Highlighting new phylogenetic specificities of Crohn's disease microbiota. Inflamm. Bowel Dis. 2011, 17, 185-192.
33. Claesson, M.J.; O'Sullivan, O.; Wang, Q.; Nikkila, J.; Marchesi, J.R.; Smidt, H.; de Vos, W.M.; Ross, R.P.; O'Toole, P.W. Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PLoS ONE 2009, 4, e6669.
34. Lam, Y.Y.; Ha, C.W.; Campbell, C.R.; Mitchell, A.J.; Dinudom, A.; Oscarsson, J.; Cook, D.I.; Hunt, N.H.; Caterson, I.D.; Holmes, A.J.; et al. Increased gut permeability and microbiota change associate with mesenteric fat inflammation and metabolic dysfunction in diet-induced obese mice. PLoS ONE 2012, 7, e34233.
35. Claesson, M.J.; Jeffery, I.B.; Conde, S.; Power, S.E.; O'Connor, E.M.; Cusack, S.; Harris, H.M.; Coakley, M.; Lakshminarayanan, B.; O'Sullivan, O.; et al. Gut microbiota composition correlates with diet and health in the elderly. Nature 2012, 488, 178-184.
36. Hu, J.; Iragavarapu, S.; Nadkarni, G.N.; Huang, R.; Erazo, M.; Bao, X.; Verghese, D.; Coca, S.; Ahmed, M.K.; Peter, I. Location-Specific Oral Microbiome Possesses Features Associated With CKD. Kidney Int. Rep. 2018, 3, 193-204.
37. Hallan, S.; Afkarian, M.; Zelnick, L.R.; Kestenbaum, B.; Sharma, S.; Saito, R.; Darshi, M.; Barding, G.; Raftery, D.; Ju, W.; et al. Metabolomics and Gene Expression Analysis Reveal Down-regulation of the Citric Acid (TCA) Cycle in Non-diabetic CKD Patients. EBioMedicine 2017, 26, 68-77.
38. Petriz, B.A.; Castro, A.P.; Almeida, J.A.; Gomes, C.P.; Fernandes, G.R.; Kruger, R.H.; Pereira, R.W.; Franco, O.L. Exercise induction of gut microbiota modifications in obese, non-obese and hypertensive rats. BMC Genom. 2014, 15, 511.
39. Rothschild, D.; Weissbrod, O.; Barkan, E.; Kurilshikov, A.; Korem, T.; Zeevi, D.; Costea, P.I.; Godneva, A.; Kalka, I.N.; Bar, N.; et al. Environment dominates over host genetics in shaping human gut microbiota. Nature 2018, 555, 210-215.
40. Ilhan, Z.E.; Marcus, A.K.; Kang, D.W.; Rittmann, B.E.; Krajmalnik-Brown, R. pH-Mediated Microbial and Metabolic Interactions in Fecal Enrichment Cultures. mSphere 2017, 2, e00047-17.
Ji Eun Kim1,2, Hyo-Eun Kim3, Ji In Park4, Hyunjeong Cho5, Min-Jung Kwak6, Byung-Yong Kim6, Seung Hee Yang7, Jung Pyo Lee7,8, Dong Ki Kim1,7, Kwon Wook Joo1,7, Yon Su Kim1,7, Bong-Soo Kim9,* and Hajeong Lee1,*
1Department of Internal Medicine, Seoul National University Hospital, Seoul 03080, Korea
2Department of Internal Medicine, Korea University Guro Hospital, Seoul 08308, Korea
3Seoul National University Hospital Biomedical Research Institute, Seoul 03082, Korea
4Department of Internal Medicine, Kangwon National University Hospital, Chuncheon 24295, Korea
5Department of Internal Medicine, Chungbuk National University Hospital, Cheongju 28644, Korea
6ChunLab, Inc., Seoul 06725, Korea
7Kidney Research Institute, Seoul National University, Seoul 08826, Korea
8Department of Internal Medicine, Seoul National University Boramae Medical Center, Seoul 07061, Korea
9Department of Life Science, Multidisciplinary Genome Institute, Hallym University, Chuncheon 24252, Korea
*Authors to whom correspondence should be addressed.
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Cho, Hyunjeong
; Min-Jung, Kwak; Byung-Yong, Kim; Yang, Seung Hee; Lee, Jung Pyo; Kim, Dong Ki; Kwon, Wook Joo; Yon Su Kim
; Bong-Soo, Kim
; Lee, Hajeong