Skip to main content

Dysbiosis of the gut microbiome is associated with CKD5 and correlated with clinical indices of the disease: a case–controlled study



Chronic kidney disease (CKD) is a universal chronic disease in China. The balance of the gut microbiome is highly crucial for a healthy human body, especially for the immune system. However, the relationship between the gut microbiome and CKD has not yet been clarified.


A total of 122 patients were recruited for this study. Among them, 24 patients were diagnosed with CKD5 but did not receive hemodialysis therapy, 29 patients were diagnosed with CKD5 and received hemodialysis therapy and 69 were matched healthy controls. The gut microbiome composition was analyzed by a 16S rRNA (16S ribosomal RNA) gene-based sequencing protocol. High-performance liquid chromatography–electrospray ionization-tandem mass spectrometry (HPLC/ESI-MS/MS) technology was used to evaluate the levels of microbiome-related protein-binding uremic toxins level, indoxyl sulfate (IS) and p-cresyl sulfate (PCS), in the patients.


We compared the gut microbiome results of 122 subjects and established a correlation between the gut microbiome and IS and PCS levels. The results indicated that alpha and beta diversity were different in patients with CKD5 than in the healthy controls (p < 0.01). In comparison to healthy controls, CKD5 patients exhibited a significantly higher relative abundance of Neisseria (p < 0.001), Lachnoclostridium (p < 0.001) and Bifidobacterium (p < 0.001). Faecalibacterium (p < 0.001) displayed a notably lower relative abundance for CKD5 patients both with and without hemodialysis than for controls. It was also found that the concentrations of IS and PCS were correlated with the gut microbiome.


Our results indicate that CKD5 patients both with and without hemodialysis had dysbiosis of the gut microbiome and that this dysbiosis was associated with an accumulation of IS and PCS. These results may support further clinical diagnosis to a great extent and help in developing potential probiotics to facilitate the treatment of CKD5.


With technological advancement, the relationship between the composition of the human microbiome, including the gut microbiome and corresponding diseases has been explored gradually and studied extensively [1]. The intestinal tract, which is the biggest digestive organ in a human body, harbors trillions of bacteria whose genomes contain 100 times more gens than that of humans [2]. Human intestinal tract fosters a number of phyla including Firmicutes, Bacteroidetes, Actinobacteria and Proteobacteria [3]. Recently, it was reported that variation in the gut microbiome was connected with chronic kidney disease (CKD) [4]. Chronic kidney disease has become a public health issue of global concern due to its increasing prevalence, high treatment cost and numerous complications. According to statistical analyses, the prevalence of CKD is 10.8% in China and 8–16% worldwide, and this number is growing [5, 6]. A study covering Europe, America and Asia has shown that the prevalence of CKD in elderly people over 64 years old is 23.4% to 35.8% [7]. As the disease progresses, kidney function is gradually lost, and eventually, the patient can receive only dialysis or kidney transplantation. Regardless of either of the treatment plans, the cost is quite large, and the quality of life of patients is greatly affected. Meanwhile, cardiovascular complications are one of the most important causes of death in patients with chronic kidney disease. In patients with stage 5 chronic kidney disease, there are many harmful substances known as uremic toxins, which are one of the main causes of cardiovascular complications. It is well known that protein-binding uremic toxin production increases in patients with CKD5, of which the most representative examples are IS and PCS. These substances have been proven to promote renal fibrosis, the risk of cardiovascular disease, insulin resistance and chronic kidney disease mineral and bone disorder (CKD-MBD) [8,9,10], and they could be influenced by the composition of the gut microbiome [11]. At the same time, the gut microbiome could affect the health of the kidney and thus even the whole body by the kidney-gut axis [12]. Some studies of changes in the gut microbiome in CKD patients have been reported [13], but the exact relationship between the gut microbiome and CKD5, especially in Chinese people, has not been validated clearly.

Therefore, in this study, we analyzed the gut microbiome of CKD5 (with and without hemodialysis) patients and healthy controls by sequencing the 16S ribosomal RNA (rRNA) gene. The results revealed effective biomarkers of the disease while revealing significant information about variation in the gut microbiome based on structure, composition and function. Some specific microbial biomarkers of CKD5 were evaluated and correlated with certain clinical indices.


Study subjects

Following approval by the Medical Ethics Committee of Shandong Provincial Qianfoshan Hospital (2017-S-015), we registered on the Chinese Clinical Trial Registry (ChiCTR1800015221) and recruited 122 individuals at Shandong Provincial Qianfoshan Hospital based on the diagnostic criteria for CKD in adults outlined in the KDIGO [14]. Among these individuals, 24 patients suffered from CKD5 but did not receive hemodialysis (CKD5-NHD group, C), 29 patients were diagnosed with CKD5 and had received regular dialysis therapy for more than 3 years (three times per week, CKD5-HD group, H), and 69 were healthy controls (normal group, N). All patients received regular phosphorus binder treatment, and all individuals resided in local areas and had similar eating habits. All individuals were fully informed of the experimental content and signed an informed consent form. Basic information from the clinical registration included the following: (1) general information: name, age, gender, contact information and place of birth; (2) history of medication: use of antibiotics, prebiotics/probiotics, medicinal activated carbon and other similar drugs within 1 month; and (3) clinical information: estimated glomerular filtration rate (eGFR) (according to the CKD-EPI formulation), serum creatinine (SCr), blood urea nitrogen (BUN), the presence of constipation, diarrhea or other types of intestinal diseases, family genetic disease or other diseases and infectious diseases, such as viral hepatitis.

The exclusion criteria included the following: (1) long-term constipation, chronic diarrhea or other types of intestinal diseases; (2) antibiotic, probiotic, prebiotic or synbiotic use 4 weeks before sample collection; (3) collection of specimens, including but not limited to polyethylene glycol electrolyte dispersants, enemas and other laxative drugs; (4) family hereditary kidney disease; and (5) viral hepatitis or other infectious diseases. The basic information for the patients and healthy controls is shown in Table 1.

Table 1 Statistics of the basic data collected from the patients and healthy controls

Biospecimen collection, DNA extraction and sequencing

Fecal specimens were collected in sterilized 2-ml tubes containing pure ethanol on ice, immediately frozen (within 30 min) and stored at − 80 °C until analysis. Genomic DNA was extracted using the cetyl trimethylammonium bromide (CTAB) method [15]. A NanoDrop 2000 (Thermo Electron Corporation, USA) spectrophotometer was used to determine the concentration of the extracted DNA. The V1–V2 regions of the 16S rRNA gene were amplified and sequenced on an Illumina HiSeq 2500 system. According to the reported bacterial 16S rRNA gene sequences in GenBank, we used DNAMAN V6 (Chinese version) software (Lynnon Biosoft Company, USA) to compare the homology of the bacteria and chose the highly conserved V1–V2 region to design a pair of primers (F: 5′-AGAGTTTGATCMTGGCTCAG-3′ and R: 5′-GCTGCCTCCCGTAGGAGT-3′) according to the basic principles of primer design [16]. PCR amplification was performed on the target fragments. Follow-up experiments were carried out according to the sequencing manual.

16S rRNA gene sequence analysis

The 16S sequencing paired-end data set was joined and quality filtered using the Laser FLASH method [17]. All sequences were analyzed using the Quantitative Insights into Microbial Ecology (QIIME, version 1.9.1) software suite [18]. Sequences were clustered against the Greengenes (13.8 release) ribosomal database’s 97% reference data set ( Sequences that remained unmatched with any of the entries in this reference were subsequently clustered into de novo OTUs at 97% similarity with the UCLUST algorithm. Taxonomy was assigned to all OTUs using the RDP classifier [19] within QIIME and the Greengenes reference data set. Rarefaction and rank abundance curves were calculated from OTU tables using alpha diversity and rank abundance scripts within the QIIME pipeline. The hierarchical clustering based on population profiles of the most common and abundant taxa was performed using UPGMA (unweighted pair group method with arithmetic mean, also known as average linkage) clustering on the distance matrix of OTU abundance. This method resulted in a Newick-formatted tree, which was obtained utilizing the QIIME package. Furthermore, Calypso online software (version 8.20, and R software (version 3.5.1) were used to analyze alpha diversity (Shannon, ACE, Chao1), beta diversity [unweighted/weighted UniFrac, principal coordinate analysis (PCoA)] and correlations between biomarkers and clinical indices. Linear discriminant analysis (LDA) and effect size (LEfSe) were performed with the online software Galaxy (

IS and PCS determination and analysis

In the morning, approximately 4 ml of fasting (for at least 8 h before the blood collection) blood per individual was collected into a vacuum blood collection tube. The collected blood samples were allowed to stand for 30 min and then centrifuged at 1300×g for 10 min at room temperature using a conventional centrifuge. The supernatant (serum) of the sample after centrifugation was placed in a 1.5 ml cryotube. Then, 50.0 mg of an IS standard (purity: 99.8%, APExBIO Company, USA) was accurately weighed and diluted with methanol in a 50 ml volumetric flask (1 mg/ml). Next, 1.0 ml of the above solution was accurately pipetted, placed in a 10 ml volumetric flask, and diluted to the mark (100 μg/ml) with methanol. The solution was then diluted with physiological saline to an IS standard solution of 0.5, 1, 2, 10, 20, and 50 μg/ml. After that, 10.0 mg of hydrochlorothiazide standard (purity: 99.5%, APExBIO Company) was accurately weighed, diluted with methanol in a volumetric flask of 100 ml, and diluted to 0.05 μg/ml with methanol. The same method was used to obtain a PCS standard solution, and the internal standard for PCS was warfarin (purity: 99.2%, APExBIO Company) rather than hydrochlorothiazide. All solutions were stored at 4 °C in the dark. The participants’ sera were accurately absorbed in 100 µl and placed in a 1.5 ml centrifuge tube, the 500 µl internal standard (0.05 µg/ml hydrochlorothiazide or warfarin methanol solution) was precisely added and then vortexed for 1 min and centrifuged at 4500×g for 10 min at 4 °C. Next, 100 µl of the supernatant was mixed with an equal volume of water, placed in an auto injection bottle, and injected 5 µl, and the IS content in human serum was calculate by the internal standard method. HPLC/ESI-MS/MS conditions were as follows. (1) Chromatographic conditions: Shimadzu Shim-pack XR-ODS (75 mm * 2.0 mm, 2.1 μm); mobile phase A (0.1% aqueous formic acid, containing 2 mM ammonium acetate): B (gradient elution with acetonitrile); flow rate of 0.3 ml/min; column temperature of 40 °C; and injection volume of 5 μl. (2) Mass spectrometry conditions: the ion source was an electrospray ionization (ESI) source; the detection method was negative ion detection; the scanning mode was multireaction monitoring (MRM); the injection voltage was − 4500 V; the ion source temperature was 550 °C; the air curtain gas flow rate was 30 Psi; the atomizing gas flow rate was 55 Psi; the collision gas flow rate was 8 Psi; and the scanning time was 0.2 s. The ion conditions for quantification and characterization are shown in Tables 2 and 3.

Table 2 IS and internal standard mass spectrometry ion conditions
Table 3 PCS and internal standard mass spectrometry ion conditions

The IS and PCS standard solutions of 0.5, 1, 2, 10, 20, and 50 μg/ml were used separately and processed according to the abovementioned operation steps for HPLC analysis. A standard curve (weighting coefficient was 1/C) was fitted by a linear regression method with the ratio of IS (PCS) to the internal standard peak area (Y) versus IS (PCS) concentration (C, μg/ml).

Statistical analysis

The clinical characteristics of the subjects are represented as the median (quartile interval) [M(IQR)] and the relative numbers, depending on the data and distribution types. Depending on the characteristics of the measurement data, one-way ANOVA was applied to multiple groups. The data from multiple groups of classified data were compared with the R*C Chi square test. Diversity categorization of alpha and beta diversity was defined in the OTU table to a sequencing depth of 20,000 per sample. Moreover, alpha diversity was determined using the Mann–Whitney U test, and beta diversity was acquired by ANOSIM (analysis of similarities). LEfSe combines the Kruskal–Wallis test or pairwise Wilcoxon rank sum test with linear discriminant analysis (LDA), whose threshold value on the logarithmic LDA score equals 3.0. Spearman’s rank correlation method was used to analyze the relationship between the microbiome and clinical indices. Analyses were performed using the SPSS statistical package (version 17.0), R software (version 3.5.1), Galaxy online software, Calypso online software (version 8.20) and GraphPad Prism 7 software. p values less than 0.05 were considered statistically significant.


To establish the gut microbiome characteristics of patients with CKD5, we conducted 16S rDNA gene sequencing to analyze 122 fecal samples from 122 individuals (24 CKD5 patients without dialysis, 29 CKD5 patients with dialysis and 69 healthy controls). The clinical characteristics of the subjects are summarized in Table 1. No significant differences were noticed in either age or gender among the three groups. Through preprocessing of the sequencing data, 6,577,239 high-quality sequences (Phred ≥ Q30) with an average of 548,103 sequences per sample were obtained.

Gut microbiome composition of the CKD5-NHD, CKD5-HD and healthy control subjects

To characterize the community attributes of the gut microbiome of individuals with CKD5, relative taxon abundances of the phyla Firmicutes and Bacteroidetes of the CKD5-NHD, CKD5-HD and healthy control groups were compared (Fig. 1a, b). Venn diagrams distinguished between the gut microbiomes of the three groups (Fig. 1c). Within the healthy control group, 227 unique species were found. Similarly, 5 unique species were found in the CKD5-NHD group, while 8 unique species were found in the CKD5-HD group (details shown in Additional file 1).

Fig. 1

Comparison of relative taxon abundances of the phyla Firmicutes and Bacteroidetes among the CKD5-NHD, CKD5-HD, and healthy control subjects. a Comparison of relative abundances of the phylum Firmicutes among the CKD5-NHD, CKD5-HD, and healthy control subjects. b Comparison of relative abundances of the phylum Bacteroidetes among the CKD5-NHD, CKD5-HD, and healthy control subjects. c Venn diagram. C CKD5-NHD group, H CKD5-HD group, N normal group

Variations in the gut microbiome among the three groups

We analyzed the alpha and beta diversity to explore the variations in the gut microbiome among the three groups. Significant differences were found among the healthy control, CKD5-NHD and CKD5-HD groups with respect to microbial diversity (Shannon index, p < 0.001) and microbial abundance (ACE index, p < 0.001; Chao1 index, p < 0.001). Microbial diversity and abundance were higher in the healthy control group than in both CKD5 groups (Fig. 2). In addition, we calculated the unweighted and weighted UniFrac distances to investigate the correlation of the gut microbiome between the CKD5 groups and the healthy control group. These calculations revealed that the healthy control group was significantly different from the CKD5-NHD group and CKD5-HD group (Fig. 3). Together, the results indicated a significant variation in the gut microbiome between the CKD5 disease and healthy control groups.

Fig. 2

Phylogenetic diversity of the gut microbiome among the CKD5-NHD, CKD5-HD and healthy control subjects. Box plots depict microbiome diversity and abundance differences according to the Shannon (a), ACE (b) and Chao1 (c) indices among the CKD5-NHD, CKD5-HD and healthy control subjects. C CKD5-NHD group, H CKD5-HD group, N normal group

Fig. 3

PCoA of the microbiome among the CKD5-NHD, CKD5-HD and healthy control subjects. a Weighted UniFrac PCoA. b Unweighted UniFrac PCoA. C CKD5-NHD group, H CKD5-HD group, N normal group

Bacterial taxonomic differences among the three groups

To select biomarkers from the disease groups, we used LEfSe to analyze and distinguish the composition of the gut microbiome between the CKD5 and healthy control groups. The gut microbiome of the CKD5-NHD group was characterized by a dominance of Lachnoclostridium, Neisseria, etc., whereas the microbiome in the healthy control group was dominated by the genus Faecalibacterium, etc. (p < 0.05, Fig. 4). The CKD5-HD group’s microbiome was also characterized by a dominance of Lachnoclostridium, Neisseria, etc., whereas the microbiome in the healthy control group displayed a dominance of Faecalibacterium, etc. (p < 0.05, Fig. 5). The images of Figs. 4a and 5a are two cut-off images from the integral LEfSe graph to show the major different taxa at the genus level. The integral graphs are shown in Additional files 2 and 3.

Fig. 4

Characteristics of the microbial community composition in the CKD5-NHD and healthy control groups. a The major differentially abundant taxa between CKD5-NHD and healthy control subjects (LDA score above 3.0), which was generated from LEfSe analysis. b The enriched taxa in the CKD5-NHD and healthy control gut microbiomes are represented in the cladogram. The central point represents the root of the tree (bacteria) and each ring represents the next lower taxonomic level (order to genus: D_3 order, D_4 family, D_5 genus). The diameter of each circle represents the relative abundance of the taxon. ce Relative abundance of taxa between the CKD5-NHD and healthy control subjects was compared, and all the p values were less than 0.0001 (Mann–Whitney U test). C CKD5-NHD group, N normal group

Fig. 5

Characteristics of the microbial community composition in the CKD5-HD and healthy control groups. a The major differentially abundant taxa between CKD5-HD and healthy control subjects (LDA score above 3.0), which was generated from LEfSe analysis. b The enriched taxa in the CKD5-HD and healthy control gut microbiomes are represented in the cladogram. The central point represents the root of the tree (bacteria), and each ring represents the next lower taxonomic level (order to genus: D_3 order, D_4 family, D_5 genus). The diameter of each circle represents the relative abundance of the taxon. ce Relative abundance of taxa between the CKD5-HD and healthy control subjects was compared (Mann–Whitney U test). H CKD5-HD group, N normal group

The gut microbiome of CKD5 patients was associated with clinical indices

With HPLC/ESI-MS/MS technology, the IS and PCS concentrations and chromatogram (shown in Additional file 4) of serum samples in each group were obtained. The details of IS and PCS concentrations are shown in Table 4. The Spearman’s rank correlation coefficient method was used to evaluate the correlation between each subject’s gut microbiome and clinical index parameters, including eGFR, SCr, BUN, IS and PCS. The details are shown in Fig. 6. Strong correlations (correction r > 0.6 or r < − 0.6, p < 0.05) were found among 39 taxa and the 5 clinical indices in CKD5-NHD subjects (Fig. 6a). Similarly, a relationship between the biomarkers (41 taxa) and clinical indices (5 indices) in CKD5-HD subjects was established (Fig. 6b). IS, PCS, BUN and CR exhibited a significant positive correlation with the genera Enterococcus and Bifidobacterium, whereas eGFR showed a significant negative correlation with the two bacterial genera. However, no significant correlations (|r| < 0.4, p < 0.05) were found between the gut microbiome and clinical indices in the normal group. The details are described in the figure below.

Table 4 The IS and PCS concentrations in each group
Fig. 6

Heat map of Spearman’s correlation analysis between the gut microbiome of the three groups and the clinical indices. a Heat map of Spearman’s correlation analysis between the gut microbiome of CKD5-NHD subjects and the clinical indices. b Heat map of Spearman’s correlation analysis between the gut microbiome of CKD5-HD subjects and the clinical indices. c Heat map of Spearman’s correlation analysis between the gut microbiome of the normal group subjects and the clinical indices. *p < 0.05, **p < 0.01, ***p < 0.001


In this study, an analysis of gut microbiome composition was performed for stool samples from subjects affected by CKD5 and matched healthy controls. It was found that the gut microbiome of patients with both CKD5-NHD and CKD5-HD demonstrated lower microbial richness as well as distinct composition compared to those of the healthy control group subjects. It was further revealed that the gut microbiome might be associated with certain enterogenous protein-binding uremic toxins. Our pilot study reported that the composition of the gut microbiome was correlated with CKD5 with and without hemodialysis.

Several studies have demonstrated that dysbiosis of human microbiome could be primarily caused by inflammatory disorders throughout the body [20,21,22,23,24]. In our study, we observed that the microbial richness was obviously lower in both the CKD5-NHD and CKD5-HD groups than in the healthy control group. A number of studies have suggested that greater overall diversity implies better health [25]. The assembled results indicate that an opportunistic pathogen may colonize patients suffering from CKD5 disease.

With respect to gut the microbiome composition, the relative abundance of Lachnoclostridium (p < 0.001) and Neisseria (p < 0.001) was significantly higher in both the CKD5-NHD and CKD5-HD groups than in the health control group. There were some differences between our results and the discovery of Vaziri et al. [13]. With regard to these differences, race, dietary habit, living environment and some other aspects may play critical roles. In our study, Neisseria levels were found to be elevated in the setting of esophageal inflammation, which was correlated with inflammatory disorders in pancreatic disease [26,27,28,29]. As we all know, patients with CKD are under the condition of a microinflammatory state. Neisseria may aggravate this pathological condition. Several studies also demonstrated that reduced bone volume of fat diet-induced bone loss in male mice was associated with the presence of Lachnoclostridium [30]. There were situations of dyslipidemia and mineral bone disorder (CKD-MBD) in patients with CKD5; Lachnoclostridium may be one of the culprits. Interestingly, Bifidobacterium, usually known as a kind of probiotic, was found to be significantly higher in the disease groups than in the healthy group. As we know, butyrate (BT) is one of the short chain fatty acids (SCFAs) that can maintain the health of the intestinal tract [31]. Adequate butyrate has been proven to be critical for sustaining the health of the gut microbiome. Butyrate not only maintains the integrity of the intestinal epithelium but also provides energy for respiration of intestinal epithelial cells and regulates immune responses of the intestinal tract [32, 33]. For instance, BT develops an antiinflammatory effect by inhibiting the recruitment and proinflammatory activity of immune cells, such as neutrophils. It could activate the differentiation of colonic Treg cells to suppress inflammation [34]. However, when the abundance of some bacteria, such as Bifidobacterium, increased, the abundance of bacteria that generate butyrate was significantly reduced which may damage epithelial cells and deteriorate the microinflammatory state. Vaziri et al. [35,36,37] found that ESRD patients exhibited significant expansion of bacterial families possessing urease, which resulted in the accumulation of some uremic toxins, such as IS and PCS, and the destruction of tight junctions. Alterations in the gut microbiome in both the CKD5-NHD and CKD5-HD groups may be due to intestinal microecological environmental changes. These changes may occur for several of the following reasons. The generally poor nutritional status of patients with injured renal function makes it difficult for a microbiome that is beneficial to the body to obtain sufficient nutritional substrates, which inhibits the growth of advantageous bacterial groups [38]. Because of restraint of the advantageous bacterial groups, intestinal peristalsis slows, intestinal microvilli are injured, and the intestinal clearance ability is reduced, which provides opportunities for pathogenic bacteria to contact and adhere to mucous membranes; this possibility is also an important reason for the changes in the microbiome structure [39]. Renal dysfunction leads to the accumulation of waste material, and thus the intestinal mucosa becomes ischemic, the mucosal barrier is damaged, and the permeability of the mucous membrane increases; moreover, the growth environment for the bacteria changes due to underlying inflammation [40]. Our results showed an increase in these three genera of bacteria, which indicates that these bacteria could be a factor that affects renal function.

The gut microbiome plays a critical role in substance metabolism and influences the essential diagnosis and treatment of various pivotal diseases, such as cancer, diabetes and CKD. In this study, our results indicated that the gut microbiome was associated with the disease indices. IS and PCS were found to have a significant positive association with Enterococcus (p < 0.01) and Bifidobacterium (p < 0.001). It is known that IS and PCS are secreted and excreted primarily through the organic anion transporter (OAT) in renal tubules. The renal excretion function of CKD patients is impaired, resulting in a reduction in the excretion of the substance and its accumulation in the body [41]. Simultaneously, due to the structure and function of the microbiome disturbance, overgrowth of Enterococcus, Bifidobacterium and other phenolic- and indole-producing bacteria occurred in CKD patients, resulting in increased levels of IS, PCS and other harmful metabolites and further leading to the accumulation of these two substances. Studies have shown that IS and PCS promote renal fibrosis, myocardial apoptosis, and insulin resistance and affect osteoblast function [42,43,44,45,46,47,48]. BUN and CR are the two main indices that reflect kidney function and are filtrated by the glomerulus and partly excreted by the renal tubules. With renal (the decline of eGFR) and gut bacteria metabolism dysfunction, the concentrations of BUN and CR were increased. Due to the disordered internal environment, CKD patients usually have an intestinal microecological imbalance (i.e., proliferation of a large number of harmful bacteria and a reduction in beneficial bacteria), resulting in a state of long-term microinflammation within the body. Disordered microbiomes contribute to the state of CKD by causing the accumulation of IS, PCS, BUN, CR and other uremic toxins, thereby aggravating the illness and increasing the risk of a variety of complications, such as cardiovascular disease, and eventually increasing all-cause mortality. Therefore, modifying the gut microbiome by prebiotics may be a novel way to improve the disease state [49, 50].


This study indicates that variation in the gut microbiome is associated with CKD5 disease. It also suggests the potential pathogens active in CKD. These results may be essential for clinical diagnosis and could be further used to develop potential probiotics that could facilitate the treatment of CKD. However, it remains unknown whether renal disorders cause dysbiosis or vice versa. Further studies need to scrutinize and explore the mechanisms between the gut microbiome and CKD using modeled organisms.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.



chronic kidney disease


ribosomal RNA


high-performance liquid chromatography–electrospray ionization-tandem mass spectrometry


indoxyl sulfate


p-Cresyl sulfate


chronic kidney disease mineral and bone disorder


estimated glomerular filtration rate


serum creatinine


blood urea nitrogen


cetyl trimethylammonium bromide


principal coordinate analysis


linear discriminant analysis


linear discriminant analysis effect size


short chain fatty acids


organic anion transporter


  1. 1.

    Gilbert JA, Blaser MJ, Caporaso JG, Jansson JK, Lynch SV, Knight R. Current understanding of the human microbiome. Nat Med. 2018;24(4):392–400.

    CAS  Article  Google Scholar 

  2. 2.

    De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, Collini S, Pieraccini G, Lionetti P. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci USA. 2010;107(33):14691–6.

    Article  Google Scholar 

  3. 3.

    Dethlefsen L, McFall-Ngai M, Relman DA. An ecological and evolutionary perspective on human-microbe mutualism and disease. Nature. 2007;449(7164):811–8.

    CAS  Article  Google Scholar 

  4. 4.

    Edwards JK. Microbiology: altered gut microbiota in CKD. Nat Rev Nephrol. 2016;12:126.

    Google Scholar 

  5. 5.

    Zhang L, Wang F, Wang L, et al. Prevalence of chronic kidney disease in China: a cross-sectional survey. Lancet. 2012;379(9818):815.

    Article  Google Scholar 

  6. 6.

    Jha V, Garcia-Garcia G, Iseki K, et al. C hronic kidney disease: global dim ension and perspectives. Lancet. 2013;382(9888):260–72.

    Article  Google Scholar 

  7. 7.

    Zhang QL, Rothenbacher D. Prevalence of chronic kidney disease in population-based studies: systematic review. BMC Public Health. 2008;8(1):117–29.

    Article  Google Scholar 

  8. 8.

    Lee SB, Kalluri R. Mechanistic connection between inflammation and fibrosis. Kidney Int. 2010;78(119):S22–6.

    Article  Google Scholar 

  9. 9.

    Huang Y, Wongamorntham S, Kasting J, McQuillan D, Owens RT, Yu L, Noble NA, Border W. Renin increases mesangial cell transforming growth factor-beta1 and matrix proteins through receptor-mediated, angiotensin II-independent mechanisms. Kidney Int. 2006;69(1):105–13.

    CAS  Article  Google Scholar 

  10. 10.

    Schepers E, Meert N, Glorieux G, et al. P-cresylsulphate, the main in vivo metabolite of p-cresol, activates leucocyte free radical production. Nephrol Dial Transplant. 2007;22(2):592–6.

    CAS  Article  Google Scholar 

  11. 11.

    Stockler-Pinto MB, Soulage CO, Borges Natália A, Cardozo LFMF, Dolenga CJ, Nakao LS, et al. From bench to the hemodialysis clinic: protein-bound uremic toxins modulate nf-κb/nrf2 expression. Int Urol Nephrol. 2018;50(2):347–54.

    CAS  Article  Google Scholar 

  12. 12.

    The Kidney-Gut Axis. Implications for Nutrition Care. J Renal Nutr. 2015;25(5):S1051227615000199.

    Google Scholar 

  13. 13.

    Vaziri ND, Wong J, Pahl M, Piceno YM, Yuan J, DeSantis TZ, Ni Z, Nguyen TH, Andersen GL. Chronic kidney disease alters intestinal microbial flora. Kidney Int. 2013;83(2):308–15.

    Article  Google Scholar 

  14. 14.

    Levin A, Stevens PE, Bilous RW, Coresh J, De Francisco ALM, De Jong PE, et al. Kidney disease: improving global outcomes (KDIGO). KDIGO 2012 clinical practice guideline for evaluation and management of chronic kidney disease. Kidney Int Suppl. 2013;3(1):1–150.

    Article  Google Scholar 

  15. 15.

    Wang X, Zhang L, Wang Y, Liu X, Zhang H, Liu Y, Shen N, Yang J, Gai Z. Gut microbiota dysbiosis is associated with Henoch-Schonlein Purpura in children. Int Immunopharmacol. 2018;58:1–8.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Moreno C, Romero J, Espejo RT. Polymorphism in repeated 16S rRNA genes is a common property of type strains and environmental isolates of the genus Vibrio. Microbiology. 2002;148(4):1233–9.

    CAS  Article  Google Scholar 

  17. 17.

    Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.

    CAS  Article  Google Scholar 

  18. 18.

    Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6.

    CAS  Article  Google Scholar 

  19. 19.

    Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM, Tiedje JM. The ribosomal database project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009;37(Database issue):D141–5.

    CAS  Article  Google Scholar 

  20. 20.

    Huck O, Al-Hashemi J, Poidevin L, Poch O, Davideau JL, Tenenbaum H, Amar S. Identification and characterization of microrna differentially expressed in macrophages exposed to porphyromonas gingivalis infection. Infect. Immunity. 2017;85(3):e00771-16.

    Article  Google Scholar 

  21. 21.

    Pere-Vedrenne C, Prochazkova-Carlotti M, Rousseau B, He W, Chambonnier L, Sifre E, Buissonniere A, Dubus P, Megraud F, Varon C, Menard A. The cytolethal distending toxin subunit cdtb of helicobacter hepaticus promotes senescence and endoreplication in xenograft mouse models of hepatic and intestinal cell lines. Front Cell Infect Microbiol. 2017;7:268.

    Article  Google Scholar 

  22. 22.

    Zitvogel L, Ayyoub M, Routy B, Kroemer G. Microbiome and anticancer immunosurveillance. Cell. 2016;165(2):276–87.

    CAS  Article  Google Scholar 

  23. 23.

    Rajagopala SV, Vashee S, Oldfield LM, Suzuki Y, Venter JC, Telenti A, Nelson KE. The human microbiome and cancer. Cancer Prev. Res. 2017;10(4):226–34.

    Article  Google Scholar 

  24. 24.

    Meng S, Chen B, Yang J, Wang J, Zhu D, Meng Q, Zhang L. Study of microbiomes in aseptically collected samples of human breast tissue using needle biopsy and the potential role of in situ tissue microbiomes for promoting malignancy. Front Oncol. 2018;8:318.

    Article  Google Scholar 

  25. 25.

    Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013;500(7464):541–6.

    Article  Google Scholar 

  26. 26.

    Benitez AJ, Hoffmann C, Muir AB, Dods KK, Spergel JM, Bushman FD, Wang ML. Inflammation-associated microbiota in pediatric eosinophilic esophagitis. Microbiome. 2015;3:23.

    Article  Google Scholar 

  27. 27.

    Macfarlane S, Furrie E, Macfarlane GT, Dillon JF. Microbial colonization of the upper gastrointestinal tract in patients with Barrett’s esophagus. Clin Infect Dis. 2007;45(1):29–38.

    Article  Google Scholar 

  28. 28.

    Si J, Lee C, Ko G. Oral microbiota: microbial biomarkers of metabolic syndrome independent of host genetic factors. Front Cell Infect. Microbiol. 2017;7:516.

    Article  Google Scholar 

  29. 29.

    Memba R, Duggan SN, Ni Chonchubhair HM, Griffin OM, Bashir Y, O’Connor DB, Murphy A, McMahon J, Volcov Y, Ryan BM, Conlon KC. The potential role of gut microbiota in pancreatic disease: a systematic review. Pancreatology. 2017;17(6):867–74.

    CAS  Article  Google Scholar 

  30. 30.

    Mccabe LR, Irwin R, Tekalur A, et al. Exercise prevents high fat diet-induced bone loss, marrow adiposity and dysbiosis in male mice. Bone. 2018;118:20–31.

    Article  Google Scholar 

  31. 31.

    Zhao J, Nian L, Kwok LY, Sun T, Zhao J. Reduction in fecal microbiota diversity and short-chain fatty acid producers in Methicillin-resistant Staphylococcus aureus infected individuals Endocrine as revealed by PacBio single molecule, real-time sequencing technology. Eur J Clin Microbiol Infect Dis. 2017;36(8):1463–72.

    CAS  Article  Google Scholar 

  32. 32.

    Pryde SE, Duncan SH, Hold GL, Stewart CS, Flint HJ. The microbiology of butyrate formation in the human colon. FEMS Microbiol Lett. 2002;217(2):133.

    CAS  Article  Google Scholar 

  33. 33.

    Rose DJ, DeMeo MT, Keshavarzian A, Hamaker BR. Influence of dietary fiber on inflammatory bowel disease and colon cancer: importance of fermentation pattern. Nutr Rev. 2007;65(2):51–62.

    Article  Google Scholar 

  34. 34.

    Furusawa Y, Obata Y, Fukuda S, Endo TA, Nakato G, Takahashi D, Nakanishi Y, Uetake C, Kato K, Kato T, Takahashi M, Fukuda NN, Murakami S, Miyauchi E, Hino S, Atarashi K, Onawa S, Fujimura Y, Lockett T, Clarke JM, Topping DL, Tomita M, Hori S, Ohara O, Morita T, Koseki H, Kikuchi J, Honda K, Hase K, Ohno H. Commensal microbederived butyrate induces the differentiation of colonic regulatory T cells. Nature. 2013;504(7480):446–50.

    CAS  Article  Google Scholar 

  35. 35.

    Wong J, Piceno YM, DeSantis TZ, Pahl M, Andersen GL, Vaziri ND. 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(3):230–7.

    CAS  Article  Google Scholar 

  36. 36.

    Vaziri ND, Yuan J, Norris K. Role of urea in intestinal barrier dysfunction and disruption of epithelial tight junction in chronic kidney disease. Am J Nephrol. 2013;37(1):1–6.

    CAS  Article  Google Scholar 

  37. 37.

    Vaziri ND, Yuan J, Rahimi A, Ni Z, Said H, Subramanian VS. Disintegration of colonic epithelial tight junction in uremia: a likely cause of CKD-associated inflammation. Nephrol Dial Transplant. 2012;27(7):2686–93.

    CAS  Article  Google Scholar 

  38. 38.

    Guarner C, González-Navajas JM, Sánchez E, Soriando G, Francés R, Chiva M, Zapater P, Benlloch S, Muñoz C, Pascual S, Balanzó J, Pérez-Mateo M, Such J. The detection of bacterial DNA in blood of rats with CCl4-induced cirrhosis with ascites represents episodes of bacterial translocation. Hepatology. 2006;44(3):633.

    CAS  Article  Google Scholar 

  39. 39.

    Pérez-Paramo M, Muñoz J, Albillos A, Freile I, Portero F, Santos M, Ortiz-Berrocal J. Effect of propranolol on the factors promoting bacterial translocation in cirrhotic rats with ascites. Hepatology. 2000;31(1):43–8.

    Article  Google Scholar 

  40. 40.

    Wiest R, Garciatsao G. Bacterial translocation (BT) in cirrhosis. Hepatology. 2005;41(3):422–33.

    CAS  Article  Google Scholar 

  41. 41.

    Nigam SK, Wu W, Bush KT, et al. Handling of drugs, metabolites, and uremic toxins by kidney proximal tubule drug transporters. Clin J Am Soc Nephrol. 2015;10(11):2039–49.

    CAS  Article  Google Scholar 

  42. 42.

    Sun CY, Chang SC, Wu MS. Suppression of Klotho expression by protein-bound uremic toxins is associated with increased DNA methyltransferase expression and DNA hypermethylation. Kidney Int. 2012;81(7):640–50.

    CAS  Article  Google Scholar 

  43. 43.

    Shimizu H, Bolati D, Higashiyama Y, Nishijima F, Shimizu K, Niwa T. Indoxyl sulfate upregulates renal expression of MCP-1 via production of ROS and activation of NF-κB, p53, ERK, and JNK in proximal tubular cells. Life Sci. 2012;90(13–14):525–30.

    CAS  Article  Google Scholar 

  44. 44.

    Shimizu H, Yisireyili M, Higashiyama Y, Nishijima F, Niwa T. Indoxyl sulfate upregulates renal expression of ICAM-1 via production of ROS and activation of NF-kappaB and p53 in proximal tubular cells. Life Sci. 2013;92(2):143–8.

    CAS  Article  Google Scholar 

  45. 45.

    Liabeuf S, Barreto DV, Barreto FC, Meert N, Glorieux G, Schepers E, Temmar M, Choukroun G, Vanholder R, Massy ZA, European Uraemic Toxin Work Group (EUTox). Free p-cresylsulphate is a predictor of mortality in patients at different stages of chronic kidney disease. Nephrol Dial Transplant. 2010;25(4):1183–91.

    CAS  Article  Google Scholar 

  46. 46.

    Yang K, Xu X, Nie L, Xiao T, Guan X, He T, Yu Y, Liu L, Huang Y, Zhang J, Zhao J. Indoxyl sulfate induces oxidative stress and hypertrophy in cardiomyocytes by inhibiting the AMPK/UCP2 signaling pathway. Toxicol Lett. 2015;234(2):110–9.

    CAS  Article  Google Scholar 

  47. 47.

    Koppe Laetitia, Pillon Nicolas J, Vella Roxane E, Croze Marine L, Pelletier Caroline C, Chambert Stéphane, Massy Ziad, Glorieux Griet, Vanholder Raymond, Dugenet Yann, Soula Hédi A, Fouque Denis, Soulage Christophe O. p-Cresyl Sulfate promotes insulin resistance associated with CKD. J Am Soc Nephrol JASN. 2013;24(1):88–99.

    CAS  Article  Google Scholar 

  48. 48.

    Tanaka H, Iwasaki Y, Yamato H, Mori Y, Komaba H, Watanabe H, Maruyama T, Fukagawa M. p-Cresyl sulfate induces osteoblast dysfunction through activating JNK and p38 MAPK pathways. Bone. 2013;56(2):347–54.

    CAS  Article  Google Scholar 

  49. 49.

    Kieffer DA, Piccolo BD, Vaziri ND, Liu S, Lau WL, Khazaeli M, Nazertehrani S, Moore ME, Marco ML, Martin RJ, Adams SH. Resistant starch alters gut microbiome and metabolomic profiles concurrent with amelioration of chronic kidney disease in rats. Am J Physiol Renal Physiol. 2016;310(9):F857–71.

    CAS  Article  Google Scholar 

  50. 50.

    Lau WL, Vaziri ND, Nunes ACF, Comeau AM, Langille MGI, England W, Khazaeli M, Suematsu Y, Phan J, Whiteson K. The phosphate binder ferric citrate alters the gut microbiome in rats with chronic kidney disease. J Pharmacol Exp Ther. 2018;367(3):452–60.

    CAS  Article  Google Scholar 

Download references


We thank Dr. Qunye Zhang and Dr. Chao Lin for helpful discussions. Miss Menghan Jiang, Mr. Bin Chen, and Mr. Meng Wang all helped collect the samples and perform 16S rRNA gene sequencing.


Not applicable.

Author information




YL recruited all participants, collected raw data from the mice, performed all DNA extractions, performed the statistical analyses and was a major contributor in writing the manuscript. XS and LZ processed the 16S rRNA gene sequencing data and data analysis and were also major contributors in writing the manuscript. YWL recruited the patients and helped explain some medical questions. MS performed the HPLC/ESI-MS/MS and analyzed the data. CL and YG helped perform 16S rRNA gene sequencing and mouse feeding. ZW and DX designed and supervised the project. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Dongmei Xu or Zunsong Wang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Medical Ethics Committee of Shandong Provincial Qianfoshan Hospital (2017-S-015) and registered on the Chinese Clinical Trial Registry (ChiCTR1800015221).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1.

The unique species in three groups.

Additional file 2.

LEfSe graph of group C and N.

Additional file 3.

LEfSe graph of group H and N.

Additional file 4: Fig. S1.

IS chromatogram of the CKD5-HD group. Fig. S2. PCS chromatogram of the CKD5-HD group. Fig. S3. IS chromatogram of the CKD5-NHD group. Fig. S4. PCS chromatogram of the CKD5-NHD group.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Su, X., Zhang, L. et al. Dysbiosis of the gut microbiome is associated with CKD5 and correlated with clinical indices of the disease: a case–controlled study. J Transl Med 17, 228 (2019).

Download citation


  • Gut microbiome
  • Chronic kidney disease
  • 16S rRNA gene sequencing
  • Indoxyl sulfate
  • p-Cresyl sulfate