Gut Mycobiome in Patients With Chronic Kidney Disease Was Altered and Associated With Immunological Profiles

Objectives: Mounting evidence suggests that bacterial dysbiosis and immunity disorder are associated with patients with chronic kidney disease (CKD), but the mycobiome is beginning to gain recognition as a fundamental part of our microbiome. We aim to characterize the profile of the mycobiome in the gut of CKD patients and its correlation to serum immunological profiles.

Methods and materials: Ninety-two CKD patients and sex–age–body mass index (BMI)–matched healthy controls (HCs) were recruited. Fresh samples were collected using sterile containers. ITS transcribed spacer ribosomal RNA gene sequencing was performed on the samples. An immunoturbidimetric test was used to assess the serum levels of immunological features.

Results: The CKD cohort displayed a different microbial community from that in the HC cohort according to principal coordinate analysis (PCoA). (P=0.001). The comparison of the two cohorts showed that the CKD cohort had significantly higher gut microbial richness and diversity (P<0.05). The CKD cohort had lower abundances of Candida, Bjerkandera, Rhodotorula, and Ganoderma compared to the HC cohort, while it had higher Saccharomyces (P<0.05). However, the microbial community alteration was inconsistent with the severity of kidney damage in patients, as only patients in CKD stage 1~3 had differed microbial community concerning for HCs based on PCoA (P<0.05). The serum concentration of the kappa light chain in CKD patients was positively associated with Saccharomyces, whereas the it was negatively associated with Ganoderma (P<0.05).

Conclusions: Not only was gut mycobiome dysbiosis observed in CKD patients, but the dysbiosis was also associated with the immunological disorder. These findings suggest that therapeutic strategies targeting gut mycobiome might be effective.

Introduction

Chronic kidney disease (CKD) is a kidney impairment that is irreversible and ultimately leads to end-stage renal disease and renal replacement therapy (e.g., transplantation or dialysis). Globally, 1.2 million people died from CKD in 2017, and the all-age mortality rate from CKD increased by 41.5% between 1990 and 2017 (1). Although the fact of the climbing incidence of CKD has been recognized as a public health problem, its etiologies are not currently clear. As the management of CKD includes treating the underlying etiologies, it is urgent to explore the potential pathology associated with the disease.

It is well established that CKD is associated with bacterial dysbiosis in the gut (2–8). The existence of gut microbiome alterations such as a reduction in bacterial diversity has been related to CKD (3, 4, 6, 7). Moreover, several potential pathogenic bacteria, such as Enterobacteriaceae, Escherichia-Shigella, Klebsiella, and Pseudomonas increased in CKD patients (3, 6, 7), while bacteria with probiotic features, including Akkermansia, Bifidobacterium, Blautia, Faecalibacterium, Prausnitzii, Lactobacillaceae, and Roseburia, depleted (3–5, 7).

Kidneys contribute to immune homeostasis, and the components of the immune system mediate renal disease and play a central role in the progression of CKD (9). For instance, serum complement 3 (C3), as partially originated from the kidney, was not only increased in CKD patients but also associated with a decreased estimated glomerular filtration rate (eGFR) (10). Additionally, an elevated serum C-reactive protein (CRP) level has been linked to CKD (11–13).

The kidneys are frequently targeted by pathogenic immune responses against renal autoantigens or by local manifestations of systemic autoimmunity (14). The interaction between the microbiome and immunity in health and disease might drive the pathogenesis of renal disorders (15). For example, our previous study demonstrated that the depletion of Actinobacteria and Bifidobacterium in CKD patients predicted the serum high level of free light chain lambda (FLC λ) (3).

Invasive fungal infections are a major challenge in the management of patients with renal dysfunction (16). It is reported that fungal infections complicate the course of 4%–7% of CKD patients, with a mortality of over 65% (17). Traditionally, fungi have been considered pathogens that commonly affect any organs of the body, including the colon. However, beyond the bacterial community, the human gut also harbors a large number of fungi (18). Gut fungi are causally implicated in microbiome assembly and immune development. Accumulating findings highlighted that the human mycobiome can strongly influence the host immune system (19, 20).

Given that CKD patients are usually complicated with fungal infection (16), and persistent immunity disorder is a common feature of CKD and inflammation is the body’s immune system response to an irritant (including bacteria, viruses, or fungi) (21, 22), we hypothesize that the gut mycobiome of individuals with and without CKD differs, and the differences correlate with their renal damage severity, along with patients’ immunological characteristics.

Methods and Materials

Participant Recruitment

The ethics committee of the Affiliated Wuxi No.2 Hospital of Nanjing Medical University approved this study (Ref. 2018051). Informed consent was provided by all subjects before their inclusion in the study. As sex, age, and BMI are important variables affecting the gut microbiome (23–25), we recruited ninety-two CKD patients and sex–age–BMI–matched healthy controls (HCs) were recruited from December 2018 to February 2020. The inclusion and exclusion criteria were based on our recent study on the bacterial microbiome in CKD patients (3). The diagnostic criteria for CKD include a decreased eGFR [<60 ml/min/1.73 m²) or an evidence of kidney damage such as albuminuria (albumin excretion rate ≥30 mg/24 h; urinary albumin– creatinine ratio (UACR) ≥30 mg/g], urine sediment abnormalities, electrolytes, and other abnormalities due to tubular disorders, abnormalities detected by histology, and structural abnormalities detected by the imaging or history of kidney transplantation (26). Patients with CKD who have never undergone hemodialysis were recruited. CKD patients were separated into five subgroups based on the stages of renal function at recruitment (27): stage 1, normal or high eGFR CKD (eGFR ≥ 90 ml/min/1.73m²); stage 2, mild CKD (eGFR ≥ 60~89 ml/min/1.73m²); stage 3, moderate CKD (30~59 ml/min/1.73m²); stage 4, severe CKD (15~29 ml/min/1.73m²); and stage 5, end-stage CKD (<15 ml/min/1.73m²). Subjects with kidney damage and positive urinary protein, or eGFR less than 90 ml/min/1.73 m², current illness (including diabetes and hypertension), elevated body temperature, or white blood cell count or serum C-reactive protein were excluded from the HC cohort. In addition, participants with acute intercurrent disease and infections, diarrhea, kidney transplantation, pregnancy, and breastfeeding and those who used antibiotics, probiotics, or immunosuppressive drugs within 30 days before enrollment were excluded from the present study.

Sample Collection and Procession

Fecal samples were collected in a sterile container by the patient. Approximately 30 mg of feces were transferred to a sterile bottle containing 500 μl lysis buffer containing Tris 0.1mol/L (pH 8.0), 2 mM Pleas, and 2% sodium dodecyl sulfate (SDS) (Guhe Health. com., Hangzhou, China) by a technician using sterile techniques and stored at −80°C until further processing.

The blood samples for detecting renal function and immunological features were collected on the day of fecal sample collection. An immunoturbidimetric test was used to assess the serum levels of immunological features on the day of sample collection (AU5400; Beckman Coulter, Pasadena, CA, USA). Information on clinical manifestations; concurrent diseases, such as diabetes, hypertension, and hyperlipidemia; and the current medication usage were assessed by reviewing clinical records, medical interviews, and face-to-face interviews. A Chinese version of the Food Frequency Questionnaire was used to assess food intake to measure the effects of confounding on the microbial community in CKD patients (28).

Fungal DNA Extraction

The DNeasy PowerSoil Pro Kit was used to isolate fungal genomic DNA from fecal samples according to the manufacturer’s instructions in a biological safety cabinet (QIAGEN, Hilden, Germany), with additional glass-bead beating steps performed using a Mini-Beadbeater (FastPrep; Thermo Electron, Boston, Waltham MA, USA). The amount of DNA was determined using a NanoDrop ND-1000 spectrophotometer (Thermo Electron). The integrity and size of DNA were verified by electrophoresis on a 1.0% agarose gel containing 0.5 mg/ml ethidium bromide. All DNA samples were stored at −20°C prior to further analysis.

The ITS regions were amplified using ITS1F (5’-CTTGGTCATTTAGAGGAAGTAA-3’) and ITS2 (2043R; 5’-GCTGCGTTCTTCATCGATGC-3’) primers (29). All PCR reactions were performed using Phusion High-Fidelity PCR Master Mix (Thermo Scientific Inc., Waltham, MA, USA) according to the manufacturer’s protocol and approximately 50 ng of extracted DNA per reaction. Thermocycling conditions were set at 98°C for 15 s for 1 cycle, then at 98°C for 15 s, 58°C for 15 s, then 72°C for 15 s for 30 cycles, followed by a final extension at 72°C for 1 min. Negative DNA extraction samples (lysis buffer and kit reagents only) were amplified and sequenced as contamination controls. Amplified products were purified using Agencourt AMPure XP beads (1 volume; Beckman Coulter, Pasadena, CA, USA) and samples were run on a 1% agarose gel in order to size-select gel slices of approximately 430 bp. The amount of DNA was determined using a Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA). Sequencing was performed with 2 × 150 bp on a Novaseq 6000 platform (Illumina Inc., San Diego, CA, USA).

Bioinformatic Analysis

The ITS sequence dataset was merged and demultiplexed into per-sample data using QIIME (V1.9.1) with default parameters (30). Raw sequencing reads with exact matches to the barcodes were assigned to their respective samples and identified as valid sequences. The low-quality sequences were filtered through the following criteria: the sequences of a specified length of <150 bp, sequences of average Phred scores of <20, sequences of ambiguous bases, and sequences of mononucleotide repeats of >8 bp. Paired-end reads were assembled using Vsearch (V2.4.4; -fastq_mergepairs -fastq_minovlen 0). Operational taxonomic unit (OTU) picking included dereplication (-derep_full length), cluster (-cluster_fast, -id 0.97), and the detection of chimeras (-uchime_ref). A representative sequence was selected from each OTU using default parameters. OTU taxonomic classification was conducted by Vsearch searching the representative sequences set against the UNITE 12_11 (https://unite.ut.ee) database (31).

An OTU table was further generated to record the abundance of each OTU in each sample and the taxonomy of the OTUs. A minimum library size was chosen to rarefy the OTUs in our present study, as it is critical to normalize the OTU table to eliminate any bias due to differences in the sampling sequencing depth. In addition, total sum scaling was applied to transform the OTU table into relative abundance by dividing the number of total reads of each sample. OTUs containing less than 0.01% of total sequences across all samples were discarded.

Sequence data analysis was performed using the QIIME and R packages (V3.2.0). OTU-level alpha richness and diversity indices, including Chao1, Shannon, and Simpson, were calculated using the OTU table. Beta diversity analysis was performed to investigate the structural variation of fungal communities across samples using Bray–Curtis metrics and visualized via principal coordinate analysis (PCoA) based on the permutational multivariate analysis of variance (PERMANOVA) calculated by the ‘adonis’ function.

Statistical Analysis

Pearson’s chi-square or Fisher’s exact tests were used with categorical variables; Student’s t-test and ANOVA were used on normalized continuous variables and the Wilcoxon rank-sum test on non-normal continuous variables. The P-value was adjusted for multiple comparisons using the Benjamini–Hochberg (BH) false discovery rate (FDR). Spearman correlation analysis was performed on the abundant bacterial genera (>1% relative abundances) and immunological features that differed between the CKD and HCs.

Results

Demographics

We assessed the mycobiome and serum immunological features of a total of 92 CKD patients and 92 sex-age-BMI matched HCs (Table 1). Both cohorts were 50% (n=46) women and 50% (n=46) men. The CKD cohort had higher levels of urinary protein, serum creatinine, urea nitrogen, uric acid, and eGFR (P<0.05 for all comparisons). Of the 92 CKD patients, the prevalence of renal function from stage 1~5 ranged from 23.91% to 17.39%. The CKD cohort had significantly elevated CRP, serum κ and λ light chain, and rheumatoid factor (Table 2; P<0.05). The medication usages of antihypertensive/glucocorticoid/hypoglycemic/hypolipidemic agents were not confounding factors of the difference of immunological features between the CKD patients and controls since the medication users did not differ in their immunological states from non-users (Tables S1–S4). As the nutrient intake was compared between the two cohorts, no significant difference was found (Table S5; P<0.05 for all nutrients); thus, they were not listed as confounding factors in the downstream analysis.

TABLE 1

Table 1 Characteristics of participants.

TABLE 2

Table 2 Comparison of immunological status between CKD and HC cohort.

Gut Mycobiome Was Altered in CKD Patients

To test whether the bladder microbiome differs between CKD patients and HCs, we first assessed the microbial community structure using all OTUs presenting in each sample. The PCoA of Bray–Curtis dissimilarities revealed a differential clustering between the CKD patients and HCs (Figure 1A; R² =0.086, P=0.001), reflecting a dysbiosis fungal gut microbiome in CKD patients. As only the participants in the CKD cohort were being administered medications at recruitment, we compared the CKD patients who were being administered antihypertensive/glucocorticoid/hypoglycemic/hypolipidemic agents and non-users using PCoA to clarify the confounding effects of medication usage. In addition, no significant differences were found between the users and non-users (Figures S1A, B; P>0.05), indicating that the abovementioned medication administrations were not confounding factors when we compare the microbiome between CKD patients and controls. As the mounting of studies demonstrated that gender plays a role in the human bacterial microbiome in the gut (32–34), we compared the bacterial community between the male and female subjects from either the CKD group and controls. However, we noticed that there were no differences between the men and women in both groups (Figure S2; P>0.05).

FIGURE 1

Figure 1 Microbial community, diversity, composition, and differed genera in the cohorts of CKD and HC. (A) PCoA based on Bray–Curtis distances at the OTU level showed different microbial compositions between groups of CKD patients and HCs (P < 0.05). Permutational multivariate analysis of variance (PERMANOVA) was performed for statistical comparisons of samples in the two cohorts. P-value was adjusted by the Benjamini–Hochberg FDR. (B) Bacterial richness and diversity measured by Chao1, Shannon, and Simpson were calculated at the microbial OTU level. The CKD patients had significantly higher levels of bacterial richness and diversity. The Wilcoxon rank-sum test was performed and adjusted by the Benjamini–Hochberg FDR. ** indicates P < 0.01. (C) Microbial profile at the phylum and genus levels. Sankey plot representing the overall gut mycobiome composition and corresponding abundance area for CKD patients and HCs. The taxonomic classification levels of phylum and genus are displayed. The top ten most abundant genera and their affiliated phyla are shown in the Sankey plot. (D) Microbial genera that were differentially abundant between CKD patients and HCs. Only the genera with above 1% are displayed. P-value was calculated using the Wilcoxon rank-sum test and adjusted by the Benjamini and Hochberg FDR. *, **, and *** indicate P < 0.05, P < 0.01, and P < 0.001, respectively.

When the bacterial richness and diversity were assessed, the CKD patients demonstrated significantly higher levels of Chao1, Shannon, and Simpson indices compared to those in the HC cohort (Figure 1B; P<0.05).

As shown in Figure 1C, either the CKD cohort or HC cohort was dominated by Ascomycota (28.39% vs. 35.64%) and Basidiomycota (11.61% vs. 13.04%) at the phylum level, but they did not show significant differences (P>0.05). At the genus level, four identified fungi were predominant (above 1% abundance) in the CKD cohort, including Candida (9.95%), Saccharomyces (5.60%), Aspergillus (2.94%), and Debaryomyces (2.63%). In addition, we noticed that there were more fungal genera with an abundance above 1% in the HC cohort, such as Candida (24.21%), Aspergillus (2.65%), Ganoderma (2.52%), Rhodotorula (2.08%), Auricularia (1.33%), Malassezia (1.20%), and Saccharomyces (1.12%).

When the fungal genus with an abundance above 1% was compared, the CKD cohort had significantly lower abundance of Candida, Bjerkandera, Rhodotorula, and Ganoderma compared to the HC cohort, while it had a higher level of Saccharomyces (Figure 1D; P<0.05). At the fungal species level, we did not observe any species displaying a significant difference between the CKD patients and healthy subjects (P>0.05).

Gut Mycobiome Fluctuated With CKD Stages

The CKD patients were divided into five subgroups based on CKD stages (29). Ten patients in each CKD stage and 10 healthy subjects who were individually sex–age–BMI matched were selected (Table S6). We noticed that the microbial community in patients with normal or high eGFR CKD/moderate CKD/end-stage CKD was significantly different from HCs; patients with normal or high eGFR CKD differed from mild/moderate CKD (Figure 2A; P<0.05). When the bacterial richness and diversity were assessed, only Chao1 showed significantly higher levels in patients with normal or high eGFR compared to that in patients with moderate CKD and HCs (Figure 2B; P<0.05). When microbial genera were compared among groups, Saccharomyces was significantly elevated in patients with normal or high eGFR CKD and patients with moderate CKD compared to HCs. In addition, patients with moderate CKD had significantly increased Saccharomyces than that in the patients with end-stage CKD (Figure 2C; P<0.05).

FIGURE 2

Figure 2 Microbial community, diversity, composition, and Saccharomyces in the subgroups according to CKD stages and HC. (A) PCoA based on Bray–Curtis distances at the OTU level showed different microbial compositions between the subgroups of patients’ renal function damage and HCs. Permutational multivariate analysis of variance (PERMANOVA) was performed for statistical comparisons of samples in the two cohorts. Patients with normal- or high-eGFR CKD/moderate CKD/end-stage CKD showed different microbial communities compared to HCs/mild CKD/moderate CKD (P < 0.05). P-value was adjusted by the Benjamini and Hochberg FDR. (B) Bacterial richness and diversity measured by Chao1, Shannon, and Simpson were calculated at the microbial OTU level. Chao1 showed significantly higher in normal or high eGFR CKD in relation to HCs/moderate CKD (P < 0.05). Wilcoxon rank-sum test was performed and adjusted by the Benjamini and Hochberg FDR. * and ** indicate P < 0.05 and P < 0.01, respectively. (C) Comparison of the abundances of Saccharomyes in CKD patients from normal or high eGFR CKD to end-stage CKD and HC. *, **, and *** indicate P < 0.05, P < 0.01, and P < 0.001, respectively.

Immunological Features Were Altered in CKD Patients and Associated With Gut Mycobiome

When the serum immunological features were compared, the CKD patients had significantly higher levels of serum CRP, FLC κ and FLC λ, and rheumatoid factor (Table 2; P<0.05). Next, we assessed the associations using Spearman correlation analysis between the immunological features (except for the rheumatoid factor) and fungal genera that were altered in the CKD cohort. We noticed that the CKD-enriched genus, Saccharomyces, was positively correlated to the enriched level of FLC κ in CKD patients (Figure 3, P<0.05). However, the CKD-depleted genera, Bjerkandera and Ganoderma, were negatively correlated to FLC κ. In addition, the reduction of the microbial genus, Bjerkandera, was negatively correlated to the enhanced expression of FLC λ in CKD patients. The decreased abundance of Ganoderma was also negatively correlated to the increased level of serum CRP (Figure 3, P<0.05).

FIGURE 3

Figure 3 Gut mycobiome correlation to immunological profiles. Spearman correlation analysis was performed on the abundant bacterial genera (>1% relative abundances) that displayed a significant difference between CKD patients and HCs and the disease profiles that showed significant difference between CKD patients and HCs. The correlations of two variables with values of P < 0.05 are displayed. *, and ** indicate P < 0.05 and P < 0.01 respectively.

As the rheumatoid factor was detectable only in 5/92 CKD patients and was not detected in any healthy subjects, we did not perform a correlation analysis between the microbial genus and the rheumatoid factor.

Discussion

As the association of the gut bacterial microbiome with the health state has gained attention and its correlation to immunity has been extensively explored in recent years, we firstly explored the fungal profile in CKD patients’ gut and its associations with serum immunological features. Like the bacterial microbiome, the mycobiome is also strongly affected by environmental factors (35), such as dietary exposure (36). Therefore, we collected and compared the nutrient intake of participants in our present study. However, at the current sample size, we did not observe a difference in the food intake between CKD patients and HCs.

Similar to the bacterial profile in CKD patients (2–8), the fungal community was altered in patients in our present study. However, the alteration was fluctuated in various clinical CKD stages. The CKD patients with normal or high eGFR, moderate CKD, and end-stage CKD displayed a different fungal community with respect to healthy subjects, whereas the patients with mild CKD and severe CKD did not exhibit a difference from that in the HC cohort. These findings also appeared in our previous study on the bacterial microbiome in CKD patients (3). The restoring process that appeared in patients, including the fungal and bacterial microbiome, should be examined using animal models.

The CKD patients possessed higher levels of microbial richness and a diversity of mycobiome, which was dissimilar to several previous bacterial microbiome studies in CKD patients (3, 4, 6, 7). When we divided the CKD patients into five subgroups based on their clinical stages of CKD, we found that the higher bacterial richness was only observed in the stage 1 patients with respect to HCs. In addition, the stage 1 patients also demonstrated higher bacterial richness than that in stage 3 patients. Ren Z. et al. compared the bacterial diversity in the gut in CKD patients in which the patients were grouped into three subgroups, such as group A (stage 1~2 patients), group B (stage 3~4 patients), and group C (stage 5 patients) (7). Additionally, their study demonstrated that there were no significant differences among the three subgroups (7). We cannot conclude that our findings were dissimilar to the Ren Z. et al. study (7) as we used a different grouping method that might influence the findings.

Reports using next-generation sequencing have found diverse fungal communities in all sections of the human gut, consisting mainly of the phyla Ascomycota and Basidiomycota (37–39). Similarly, our present study also found that Ascomycota and Basidiomycota dominated in the participant’s gut, either CKD patients or HCs.

It is worthy to note that Candida was the most prevalent fungus in both CKD and HC cohorts, which was significantly depleted in the patients. A similar finding was observed in Sciavilla, P and his colleague’s study (40). When they compared the fungal microbiome using culturomics in the gut in patients with irritable bowel syndrome, Candida spp. exhibited a higher frequency in the patient group compared to control (40). Di Paola M also reported that patients with Crohn’s disease had higher levels of isolated Candida spp. compared to the healthy subjects (41). Therefore, the role of Candida spp. in the fungal community in the gut should be extensively explored using animal models in the future. As gut wall inflammation is the key pathology in patients with irritable bowel syndrome and Crohn’s disease (42), and it is also common in CKD patients, it is necessary to investigate the role of Candida spp. in the gut of CKD patients.

Previous mycobiome studies frequently reported that Saccharomyces was one of the most prevalent fungi in the human gut (18, 37, 39, 43, 44). Our present study also demonstrated that it was one of most abundant fungi in the participants’ gut, especially in the CKD patients. Moreover, the high level of Saccharomyces was responsible for the high concentration of serum FLC κ in CKD patients, as the abundance of Saccharomyces was positively associated with the level of FLC κ in patients. Saccharomyces spp. is a probiotic whose clinical efficacy, anti-inflammatory, and immunomodulatory effects are supported by previous extensive studies (45, 46). As FLCs, such as FLC κ, accumulate in the circulation, their concentrations progressively increase with renal damage in CKD patients (47).

We noticed that Bjerkandera was declined in CKD patients and it was negatively correlated to the levels of FLC κ and λ in patients. It is necessary to explore the immunological function of Bjerkandera in the human gut, as it has been rarely reported by human gut microbiome studies.

In the present study, Ganoderma decreased in CKD patients’ gut and it was negatively linked to the levels of CRP and FLC κ. Ganoderma spp., a component of edible mushroom, has potential beneficial effects. Thus, the increase of Saccharomyces in the patients’ gut and its positive link to FLC κ might be a protective response in CKD progression (48, 49). For example, Chang and colleagues have shown that Ganoderma spp. reduced obesity by modulating the gut microbiome (48). Wu et al. also demonstrated that the administration of Ganoderma spp. increased the abundance of beneficial bacteria in the gut of mice (49). Therefore, the decrease of Ganoderma in CKD patients and its negative associations with CRP and FLC κ suggest that Ganoderma might play a potential role of immune modulation in the onset and progression of CKD.

Renal diseases are associated with an imbalanced bacterial microbiome and the loss of immune homeostasis (2–9). Coincidentally, our present study demonstrated that dysbiosis of the gut mycobiome is evident in CKD patients and it is associated with the immunological disorders. It is especially noteworthy that Saccharomyces and Ganoderma, the components of probiotics (45), were linked with patients’ immunological profiles. For example, a previous study reported that the administration of Saccharomyces boulardii to patients with HIV lowered the concentrations of some gut bacterial species, such as the Clostridiaceae family, which were correlated with systemic levels of bacterial translocation and inflammation markers (45). Another study demonstrated that Reishi mushroom increased the gut Bacteroides/Firmicutes ratio and promoted the growth of anti-inflammatory and short-chain fatty acid (SCFA)–producing bacteria (50). Thus, the findings in our present study shed light on mycobiome-based therapy in CKD patients.

There are limitations to our present study. On one hand, we cannot remove the bias led by the co-occurrence of diabetes and hypertension in CKD patients as renal damage is one of the most common complications of diabetes and hypertension (51,52). On the other hand, although we compared the composition of microbiome using subgroups based on the patients’ clinical stages, the small sample size might skew the findings. A future study with a large and multicenter population is needed to confirm the association between the severity of renal damage and the mycobiome.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/sra/?term=PRJNA647266.

Ethics Statement

The studies involving human participants were reviewed and approved by the ethics committee of the Affiliated Wuxi Second Hospital of Nanjing Medical University (Ref. 2018051). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

Conceptualization: LL, FL, YZ, FY, and NF. Methodology: JH, SW, YW, YF, JS, LH, CG, PJ, YT, WG, and FY. Software: FL, JH, and WG. Validation: NF. Writing: JH, FL, YZ, and YG. Supervision: NF. Funding acquisition: NF. Project administration: FL, JH, LH, PJ, JS, CG, YF, FY, and YW. All authors contributed to the article and approved the submitted version.

Funding

Wuxi “Taihu Talents Program” Medical and Health High-level Talents Project (THRCJH20200901); Wuxi “key medical discipline construction” Municipal Clinical Medical Center (municipal public health center) Project (LCYXZX202103).

Conflict of Interest

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

Publisher’s Note

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

Acknowledgments

We gratefully acknowledge the volunteers who participated in our study.

Supplementary Material

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

Supplementary Figure 1 | The effects of medication usageon CKD patients’ gut mycobiome. PCoA based on Bray-Curtis distances at OUT-level did not show different microbial compositions between CKD patients who took antihypertensive agent and those who had not take antihypertensive agent (A); who took glucocorticoid agent and those did not take glucocorticoid agent (B); who took hypoglycemic agent and those did not take hypoglycemic agent (C); who took hypolipidemic agent and those did not take hypolipidemic agent (D); Permutational multivariate analysis of variance (PERMANOVA) was performed for statistical comparisons of samples using different levels of hydroxychloroquine. P-value was adjusted by Benjamini and Hochberg false discovery rate.

References

1. Bikbov B, Purcell C A, Levey A S, Smith M, Abdoli A, Abebe M, et al. Global, Regional, and National Burden of Chronic Kidney Disease, 1990-2017: A Systematic Analysis for the Global Burden of Disease Study 2017. Lancet (2020) 395:709–33. doi: 10.1016/S0140-6736(20)30045-3

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Jiang S, Xie S, Lv D, Wang P, He H, Zhang T, et al. Alteration of the Gut Microbiota in Chinese Population With Chronic Kidney Disease. Sci Rep (2017) 7:2870. doi: 10.1038/s41598-017-02989-2

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Liu F, Xu X, Chao L, Chen K, Shao A, Sun D, et al. Alteration of the Gut Microbiome in Chronic Kidney Disease Patients and Its Association With Serum Free Immunoglobulin Light Chains. Front Immunol (2021) 12:609700. doi: 10.3389/fimmu.2021.609700

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Li F, Wang M, Wang J, Li R, Zhang Y. Alterations to the Gut Microbiota and Their Correlation With Inflammatory Factors in Chronic Kidney Disease. Front Cell Infect Microbiol (2019) 9:206. doi: 10.3389/fcimb.2019.00206

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Wang X, Yang S, Li S, Zhao L, Hao Y, Qin J, et al. Aberrant Gut Microbiota Alters Host Metabolome and Impacts Renal Failure in Humans and Rodents. GUT (2020) 69:2131–42. doi: 10.1136/gutjnl-2019-319766

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Hu X, Ouyang S, Xie Y, Gong Z, Du J. Characterizing the Gut Microbiota in Patients With Chronic Kidney Disease. POSTGRAD Med (2020) 132:495–505. doi: 10.1080/00325481.2020.1744335

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Ren Z, Fan Y, Li A, Shen Q, Wu J, Ren L, et al. Alterations of the Human Gut Microbiome in Chronic Kidney Disease. Adv Sci (Weinh) (2020) 7:2001936. doi: 10.1002/advs.202001936

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Lecamwasam A, Nelson TM, Rivera L, Ekinci EI, Saffery R, Dwyer KM. Gut Microbiome Composition Remains Stable in Individuals With Diabetes-Related Early to Late Stage Chronic Kidney Disease. Biomedicines (2020) 9:19. doi: 10.3390/biomedicines9010019

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Tecklenborg J, Clayton D, Siebert S, Coley SM. The Role of the Immune System in Kidney Disease. Clin Exp Immunol (2018) 2:142–50. doi: 10.1111/cei.13119

CrossRef Full Text | Google Scholar

10. Bao X, Borne Y, Muhammad IF, Schulz CA, Persson M, Orho-Melander M, et al. Complement C3 and Incident Hospitalization Due to Chronic Kidney Disease: A Population-Based Cohort Study. BMC Nephrol (2019) 20:61. doi: 10.1186/s12882-019-1248-7

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Adejumo O, Okaka E, Okwuonu C, Iyawe I, Odujoko O. Serum C-Reactive Protein Levels in Pre-Dialysis Chronic Kidney Disease Patientsin Southern Nigeria. Ghana Med J (2016) 1:31–8. doi: 10.4314/gmj.v50i1.5

CrossRef Full Text | Google Scholar

12. Abraham G, Sundaram V, Sundaram V, Mathew M, Leslie N, Sathiah V. C-Reactive Protein, a Valuable Predictive Marker in Chronic Kidney Disease. Saudi J Kidney Dis Transpl (2009) 20:811–5.

PubMed Abstract | Google Scholar

13. Yeo E, Hwang J, Park JE, Choi YJ, Huh KB, Kim WY. Tumor Necrosis Factor (TNF-α) and C-Reactive Protein (CRP) are Positively Associated With the Risk of Chronic Kidney Disease in Patients With Type 2 Diabetes. ymj (2010) 51:519–25. doi: 10.3349/ymj.2010.51.4.519

CrossRef Full Text | Google Scholar

14. Kurts C, Panzer U, Anders HJ, Rees AJ. The Immune System and Kidney Disease: Basic Concepts and Clinical Implications. Nat Rev Immunol (2013) 13:738–53. doi: 10.1038/nri3523

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Felix K, Richard B J, Richard AF. Immunity, Microbiota and Kidney Disease. Nat Rev Nephrol (2019) 15:263–74. doi: 10.1038/s41581-019-0118-7

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Gandhi BV, Bahadur MM, Dodeja H, Aggrwal V, Thamba A, Mali M. Systemic Fungal Infections in Renal Diseases. J POSTGRAD Med (2005) 51 Suppl 1:S30–6.

PubMed Abstract | Google Scholar

17. Jha V, Chugh S, Chugh KS. Infections in Dialysis and Transplant Patients in Tropical Countries. Kid Int (2000) 57:S85–93. doi: 10.1046/j.1523-1755.2000.07415.x

CrossRef Full Text | Google Scholar

18. Nash AK, Auchtung TA, Wong MC, Smith DP, Gesell JR, Ross MC, et al. The Gut Mycobiome of the Human Microbiome Project Healthy Cohort. MICROBIOME (2017) 5:153. doi: 10.1186/s40168-017-0373-4

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zhang I, Pletcher SD, Goldberg AN, Barker BM, Cope EK. Fungal Microbiota in Chronic Airway Inflammatory Disease and Emerging Relationships With the Host Immune Response. Front Microbiol (2017) 8:2477. doi: 10.3389/fmicb.2017.02477

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Fan D, Coughlin LA, Neubauer MM, Kim J, Kim MS, Zhan X, et al. Activation of HIF-1alpha and LL-37 by Commensal Bacteria Inhibits Candida Albicans Colonization. Nat Med (2015) 21:808–14. doi: 10.1038/nm.3871

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Yilmaz MI, Carrero JJ, Axelsson J, Lindholm B, Stenvinkel P. Low-Grade Inflammation in Chronic Kidney Disease Patients Before the Start of Renal Replacement Therapy: Sources and Consequences. Clin Nephrol (2007) 68:1–9. doi: 10.5414/CNP68001

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Chen L, Deng H, Cui H, Fang J, Zuo Z, Deng J, et al. Inflammatory Responses and Inflammation-Associated Diseases in Organs. Oncotarget (2018) 9:7204–18. doi: 10.18632/oncotarget.23208

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Holmes RP, Assimos DG. The Impact of Dietary Oxalate on Kidney Stone Formation. Urol Res (2004) 32:311. doi: 10.1007/s00240-004-0437-3

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Aleman F, Valenzano DR. Microbiome Evolution During Host Aging. PLoS Pathog (2019) 15:e1007727. doi: 10.1371/journal.ppat.1007727

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Gao X, Zhang M, Xue J, Huang J, Zhuang R, Zhou X, et al. Body Mass Index Differences in the Gut Microbiota Are Gender Specific. Front Microbiol (2018) 9:1250. doi: 10.3389/fmicb.2018.01250

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Perez-Gomez MV, Bartsch LA, Castillo-Rodriguez E, Fernandez-Prado R, Fernandez-Fernandez B, Martin-Cleary C, et al. Clarifying the Concept of Chronic Kidney Disease for Non-Nephrologists. Clin Kidney J (2019) 12:258–61. doi: 10.1093/ckj/sfz007

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Levey AS, Stevens LA, Coresh J. Conceptual Model of CKD: Applications and Implications. Am J Kidney Dis (2009) 53:S4–S16. doi: 10.1053/j.ajkd.2008.07.048

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Zhao WH, Huang ZP, Zhang X, He L, Willett W, Wang JL, et al. Reproducibility and Validity of a Chinese Food Frequency Questionnaire. BioMed Environ Sci (2010) 23:1–38. doi: 10.1016/S0895-3988(11)60014-7

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Sun B, Xia Y, Garber PA, Amato KR, Gomez A, Xu X, et al. Captivity Is Associated With Gut Mycobiome Composition in Tibetan Macaques (Macaca Thibetana). Front Microbiol (2021) 12:665853. doi: 10.3389/fmicb.2021.665853

PubMed Abstract | CrossRef Full Text | Google Scholar

30. White JR, Maddox C, White O, Angiuoli SV, Fricke WF. CloVR-ITS: Automated Internal Transcribed Spacer Amplicon Sequence Analysis Pipeline for the Characterization of Fungal Microbiota. Microbiome (2013) 1:6. doi: 10.1186/2049-2618-1-6

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Langarica-Fuentes A, Zafar U, Heyworth A, Brown T, Fox G, Robson GD. Fungal Succession in an in-Vessel Composting System Characterized Using 454 Pyrosequencing. FEMS Microbiol Ecol (2014) 88:296–308. doi: 10.1111/1574-6941.12293

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gomez A, Luckey D, Taneja V. The Gut Microbiome in Autoimmunity: Sex Matters. Clin Immunol (2015) 159:154–62. doi: 10.1016/j.clim.2015.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Razavi AC, Potts KS, Kelly TN, Bazzano LA. Sex, Gut Microbiome, and Cardiovascular Disease Risk. Biol Sex Differ (2019) 10:29. doi: 10.1186/s13293-019-0240-z

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Min Y, Ma X, Sankaran K, Ru Y, Chen L, Baiocchi M, et al. Sex-Specific Association Between Gut Microbiome and Fat Distribution. Nat Commun (2019) 10:2408. doi: 10.1038/s41467-019-10440-5

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Fiers WD, Leonardi I, Iliev ID. From Birth and Throughout Life: Fungal Microbiota in Nutrition and Metabolic Health. Annu Rev Nutr (2020) 40:323–43. doi: 10.1146/annurev-nutr-013120-043659

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Shankar J. Food Habit Associated Mycobiota Composition and Their Impact on Human Health. Front Nutr (2021) 8:773577. doi: 10.3389/fnut.2021.773577

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Hoffmann C, Dollive S, Grunberg S, Chen J, Li H, Wu GD, et al. Archaea and Fungi of the Human Gut Microbiome: Correlations With Diet and Bacterial Residents. PLoS One (2013) 8:e66019. doi: 10.1371/journal.pone.0066019

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Suhr MJ, Hallen-Adams HE. The Human Gut Mycobiome: Pitfalls and Potentials–a Mycologist's Perspective. MYCOLOGIA (2015) 107:1057–73. doi: 10.3852/15-147

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Mar RM, Perez D, Javier CF, Esteve E, Marin-Garcia P, Xifra G, et al. Obesity Changes the Human Gut Mycobiome. Sci Rep (2015) 5:14600. doi: 10.1038/srep14600

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Sciavilla P, Strati F, Di Paola M, Modesto M, Vitali F, Cavalieri D, et al. Gut Microbiota Profiles and Characterization of Cultivable Fungal Isolates in IBS Patients. Appl Microbiol Biotechnol (2021) 105:3277–88. doi: 10.1007/s00253-021-11264-4

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Di Paola M, Rizzetto L, Stefanini I, Vitali F, Massi-Benedetti C, Tocci N, et al. Comparative Immunophenotyping of Saccharomyces Cerevisiae and Candida Spp. Strains From Crohn's Disease Patients and Their Interactions With the Gut Microbiome. J Transl Autoimmun (2020) 3:100036. doi: 10.1016/j.jtauto.2020.100036

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Lau WL, Kalantar-Zadeh K, Vaziri ND. The Gut as a Source of Inflammation in Chronic Kidney Disease. NEPHRON (2015) 130:92–8. doi: 10.1159/000381990

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Shah S, Locca A, Dorsett Y, Cantoni C, Ghezzi L, Lin Q, et al. Alterations of the Gut Mycobiome in Patients With MS. EBIOMEDICINE (2021) 71:103557. doi: 10.1016/j.ebiom.2021.103557

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Hoarau G, Mukherjee PK, Gower-Rousseau C, Hager C, Chandra J, Retuerto MA, et al. Bacteriome and Mycobiome Interactions Underscore Microbial Dysbiosis in Familial Crohn's Disease. MBIO (2016) 7:e01250–16. doi: 10.1128/mBio.01250-16

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Villar-Garcia J, Guerri-Fernandez R, Moya A, Gonzalez A, Hernandez JJ, Lerma E, et al. Impact of Probiotic Saccharomyces Boulardii on the Gut Microbiome Composition in HIV-Treated Patients: A Double-Blind, Randomised, Placebo-Controlled Trial. PLoS One (2017) 12:e173802. doi: 10.1371/journal.pone.0173802

CrossRef Full Text | Google Scholar

46. Canonici A, Pellegrino E, Siret C, Terciolo C, Czerucka D, Bastonero S, et al. Saccharomyces Boulardii Improves Intestinal Epithelial Cell Restitution by Inhibiting Alphavbeta5 Integrin Activation State. PLoS One (2012) 7:e45047. doi: 10.1371/journal.pone.0045047

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Kopple JD, Massry SG, Kalantar ZK. Uremic Toxicity. In: Nutritional Management of Renal Disease. London: Academic Press (2013). 49–77. doi: 10.1016/B978-0-12-391934-2.00004-7

CrossRef Full Text | Google Scholar

48. Chang CJ, Lin CS, Lu CC, Martel J, Ko YF, Ojcius DM, et al. Ganoderma Lucidum Reduces Obesity in Mice by Modulating the Composition of the Gut Microbiota. Nat Commun (2015) 6:7489. doi: 10.1038/ncomms8489

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Wu X, Cao J, Li M, Yao P, Li H, Xu W, et al. An Integrated Microbiome and Metabolomic Analysis Identifies Immunoenhancing Features of Ganoderma Lucidum Spores Oil in Mice. Pharmacol Res (2020) 158:104937. doi: 10.1016/j.phrs.2020.104937

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Li MY, Yu LL, Zhao JX, Zhang H, Chen W, Zhai QX, et al. Role of Dietary Edible Mushrooms in the Modulation of Gut Microbiota. J Funct FOODS (2021) 83:104538. doi: 10.1016/j.jff.2021.104538

CrossRef Full Text | Google Scholar

51. Alicic RZ, Rooney MT, Tuttle KR. Diabetic Kidney Disease: Challenges, Progress, and Possibilities. Clin J Am Soc Nephrol (2017) 12:2032–45. doi: 10.2215/CJN.11491116

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Griffin KA. Hypertensive Kidney Injury and the Progression of Chronic Kidney Disease. Hypertension (2017) 70:687–94. doi: 10.1161/HYPERTENSIONAHA.117.08314

PubMed Abstract | CrossRef Full Text | Google Scholar