Gut microbiome and mycobiome in inflammatory bowel disease patients with Clostridioides difficile infection
Si Yu1† 
Xiaomeng Ge2† Hui Xu1† 
Bei Tan1 
Bowen Tian1 Yujie Shi1 Yimin Dai1 Yue Li1* 
Songnian Hu3,4* 
Jiaming Qian1- 1Department of Gastroenterology, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, China
- 2Microbial Resource and Big Data Center, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
- 3State Key Laboratory of Microbial Resources, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
- 4University of Chinese Academy of Sciences, Beijing, China
Background: Clostridium difficile infection (CDI) is common in patients with inflammatory bowel disease (IBD) and has been reported as a risk factor for poor outcome. However, gut microbiome and mycobiome of IBD patients with CDI have been barely investigated. This study aimed to assess the gut microbiome and mycobiome in IBD patients with CDI.
Methods: We collected fecal samples from patients with active IBD and concomitant CDI (IBD-CDI group, n=25), patients with active IBD and no CDI (IBD-only group, n=51), and healthy subjects (HC, n=40). Patients’ characteristics including demographic data, disease severity, and medication history were collected. Metagenomic sequencing, taxonomic and functional analysis were carried out in the samples.
Results: We found that the bacterial alpha diversity of the IBD-CDI group was decreased. The bacterial and fungal beta diversity variations between IBD patients and HC were significant, regardless of CDI status. But the IBD-CDI group did not significantly cluster separately from the IBD-only group. Several bacterial taxa, including Enterococcus faecium, Ruminococcus gnavus, and Clostridium innocuum were overrepresented in the IBD-CDI group. Furthermore, IBD patients with CDI were distinguished by several fungal taxa, including overrepresentation of Saccharomyces cerevisiae. We also identified functional differences in IBD patients with CDI include enrichment of peptidoglycan biosynthesis. The network analysis indicated specific interactions between microbial markers in IBD-CDI patients.
Conclusion: IBD patients with CDI had pronounced microbial dysbiosis. Gut micro-ecological changes in IBD patients with CDI might provide insight into the pathological process and potential strategies for diagnosis and treatment in this subset of patients.
1 Introduction
Inflammatory bowel disease (IBD) is a chronic inflammatory disorder of the gastrointestinal tract, including two main subtypes, ulcerative colitis (UC) and Crohn’s disease (CD). Dysbiosis of gut microflora is well recognized as an important factor in the pathogenesis and development of IBD (Pickard et al., 2017; Franzosa et al., 2019). Clostridioides difficile (C. difficile) is a Gram-positive, sporulating, anaerobic bacillus, widely distributed in human gut as an opportunistic pathogen (Czepiel et al., 2019). C. difficile infection (CDI) was reported to be the major cause of healthcare-associated infections around the world, and was also strongly associated with distinct alternations in the gut microbiota (Lamendella et al., 2018; Rodríguez et al., 2020; Berkell et al., 2021). Severe and recurrent CDI patients can benefit from gut microbiota reconstruction by fecal microbiota transplantation (FMT) (Hvas et al., 2019). IBD with concomitant CDI is a common condition with a prevalence up to 32% (Rodríguez et al., 2020). IBD patients with CDI were demonstrated to have more severe disease course, poor outcome, and an increase in morbidity and mortality (D’Aoust et al., 2017; Li et al., 2018).
Recent studies have tried to investigate the gut microbiome and in IBD patients with CDI. Dysbiotic microbiome and alternations of certain gut microorganisms were revealed (Sokol et al., 2018; Bushman et al., 2020; Hellmann et al., 2020; Mahnic et al., 2022). However, there is still a limited number of studies exploring the composition and changes of gut microbiome in IBD patients with CDI. Besides, most of these studies were performed with 16S rRNA gene sequencing, while metagenomic data could provide both the taxonomic and functional information in high resolution (Maccaferri et al., 2011). In addition, all the current studies have been focused on bacterial dysbiosis, while there is more and more evidence that nonbacterial microbes, such as fungi, might play unique and important roles in IBD pathogenesis and disease activity. Knowledge of the gut mycobiome might provide new insights for IBD therapeutic strategies and prognosis prediction (Chehoud et al., 2015; Sokol et al., 2017; Imai et al., 2019; Lemoinne et al., 2020; Frau et al., 2021; Ventin-Holmberg et al., 2021; Guzzo et al., 2022; Underhill and Braun, 2022; Li et al., 2022b). Nevertheless, gut mycobiome in IBD patients with CDI have never been described before.
Therefore, we conducted the study with the following aims: (1) to investigate the taxonomic and functional alternations of gut microbiome in IBD patients with CDI with metagenome sequencing; (2) to describe the gut mycobiome in IBD patients with CDI.
2 Materials and methods
2.1 Study design and sample collection
All participants were consecutively recruited from IBD patients who were hospitalized for active disease in Peking Union Medical College Hospital (Beijing, China) from November 2017 to May 2020. Diagnosis of IBD was established based on the clinical, radiological, endoscopic and histological criteria. Diagnosis of CDI was based on the presence of diarrhea and positive results of toxin A/B enzyme immunoassays and/or nucleic acid amplification test (Kucharzik et al., 2021). Participants who had taken metronidazole and/or vancomycin within a month before index hospitalization were excluded. Patients who had received fecal microbiota transplantation and/or colectomy were also excluded. Healthy volunteers who were physically fit and did not have any gastrointestinal disease, infectious disease, or other chronic disease were recruited as the healthy controls. Fecal samples of patients were collected once a patient was hospitalized and had received neither any antibiotics for CDI nor escalated medication for IBD. Fecal samples were collected in sterile collection tubes and were immediately transferred to a −80°C freezer for further analysis. Research approval was obtained from the Ethics Committees of Peking Union Medical College Hospital (approval no. JS-1494). All patients provided informed consent. The study conformed with the principles in the Declaration of Helsinki.
2.2 Shotgun sequencing and metagenomic analyses
Metagenomic DNA libraries were sequenced on Illumina platform and 150 bp paired-end reads were generated. Low-quality reads were removed with Trimmomatic (Bolger et al., 2014). Low-quality samples with a high proportion of low-quality reads and human-contaminating reads were removed in further analysis. Taxonomic classifications were assigned to metagenomic reads using MetaPhlAn3 (Beghini et al., 2021). Functional profile was obtained using HUMAnN2 (Abubucker et al., 2012). Discriminative microbial markers in each study group was identified with linear discriminant analysis effect size (Lefse) (Segata et al., 2011). Comparison and identification of differential pathways were performed using STAMP (Parks et al., 2014). Ecological network of gut microbiota in the three groups was analyzed using SparCC (Friedman and Alm, 2012). Cytoscape (version 3.9.1) was used for visualization of the networks.
2.3 Statistical analysis
R (version 4.2.1) was used for statistical analysis. Continuous variables with a non-Gaussian distribution are presented as the median and interquartile range (IQR) and were compared using the Mann-Whitney U test. Categorical data are presented as counts and percentages and were compared using chi-square test or Fisher’s exact test. Multiple group statistics was performed using Kruskal- Wallis test, and the post-hoc test was performed with Dunn’s multiple comparison test. The p values were corrected using the Benjamini and Hochberg procedure to control for the false discovery rate. P < 0.05 was considered to be significant.
3 Results
3.1 Clinical characteristics
A total of 169 participants including 129 IBD patients and 40 healthy controls were recruited. and 116 of them were included in the study population after quality control. The study population consist of three groups: patients with IBD and no CDI (IBD-only, n = 51), patients with IBD and CDI (IBD-CDI, n = 25), and healthy subjects (HC, n = 40). Clinical characteristics of the patients were shown in Supplementary Table 1. 50 patients were diagnosed with UC (65.8%), and 26 patients were diagnosed with CD (34.2%). Distribution of the diagnostic types was comparable between the IBD-CDI and IBD-only groups (Supplementary Table 1, p = 0.588). A majority of patients had moderate-to-severe disease (76.6% of the UC patients according to the Truelove-Witts criteria; 72.0% of the CD patients according to the Crohn disease activity index). Most patients had prior 5-aminosalicylic acid exposure (77.0%), while only a fraction of patients had been exposed to corticosteroid (36.8%), immunosuppressant (17.1%), biologics (14.5%), and antibiotics (18.4%) within a month before the index hospitalization. The two groups of patients were comparable, notably in terms of age, lifestyles, body mass index, diagnostic subtype, disease severity, laboratory data, and medication history. However, patients in the IBD-only group were more frequently male, and had more often perianal lesion than patients in the IBD-CDI group.
3.2 Gut microbiome profiles of the study groups
Compared with HC samples, the alpha diversity (assessed by the Shannon index) was significantly decreased in IBD-only and IBD-CDI groups. No significant difference of alpha diversity was observed between IBD-only and IBD-CDI groups (Figure 1A). The beta diversity was significantly different among all IBD patients compared with controls (assessed by the Bray-Curtis dissimilarity index, PERMNOVA P = 0.001). However, the beta diversity did not show a clustering of samples according to C. difficile infection status (Figure 1B). At the species level, the species counts of gut microbiota were decreased in both the IBD-only and IBD-CDI groups, compared with HC (Figure 1C). There were 258 common species among the three groups (Figure 1D; Supplementary Table 2).
Figure 1 Biodiversity of gut microbiome in the HC, IBD-only, and IBD-CDI groups. (A) The alpha diversity was measured in the Shannon index. (B) PCoA plots of the three study groups. (C) Box plots of the number of species among the three groups. (D) The venn diagram of the common and unique species among the three groups. ****P < 0.0001.
In the three study groups, the gut microbiome was dominated by Firmicutes, Bacteroidetes, Actinobacteria, and Proteobacteria at the phylum level (Figure 2A; Supplementary Table 3), and Lachnospiraceae, Bacteroidaceae, Bifidobacteriaceae, Ruminococcaceae, Enterobacteriaceae, Enterococcaceae, and Prevotellaceae at the family level (Figure 2B; Supplementary Table 3). The most dominant genera identified were Bacteroides, Bifidobacterium, Blautia, Faecalibacterium, Enterococcus, Prevotella, Escherichia, Roseburia, Lactobacillus, and Eubacterium (Figure 2C; Supplementary Table 3). Specifically, relative abundance of Faecalibacterium and Bacteroides were decreased, and the Enterococcus and Blautia abundance were increased in the IBD-CDI group compared with the other two groups. The HC group had higher relative abundance of Roseburia compared with the other two groups of patients (all P < 0.05, Figures 2D–F; Supplementary Table 3).
Figure 2 Gut microbiome composition in the HC, IBD-only, and IBD-CDI groups. (A-C) Global composition at the phylum, family, and genus level. (D-F) Box plots of the top 10 genera sorted by the mean values of relative abundance of the three groups.
Microbial markers among healthy control and IBD patients with or without CDI were identified with lefse analysis (Figure 3). According to the cladogram plot, the differential taxa in the IBD-CDI group included the phylum Proteobacteria, the order Enterobacterales, the class Gammaproteobacteria, and the family Enterococcaceae and Enterobacteriaceae (Figure 3A; Supplementary Table 4). At the species level, the IBD-CDI group was characterized by the enrichment of Enterococcus faecium (E. faecium), Ruminococcus gnavus (R. gnavus), Clostridium innocuum (C.innocuum), C. difficile, and Bacteroides vulgatus (Figure 3B; Supplementary Table 4). The random forest classifier was used to identify the species that best discriminated the three groups of samples from each other. As shown in Figure 3C, the two species enriched in the IBD-CDI group, C. difficile and R. gnavus, were identified as two important microbial markers.
Figure 3 Microbial marker analysis. (A) The cladogram plot of the three study groups. (B) LDA score distribution bar plot for marker species. Species with LDA score >4.0 and P < 0.05 were regarded as significantly discriminative markers. (C) Analysis with the machine learning algorithm random forest. Mean decrease accuracy indicated the importance of the markers. Only taxa present in >20% of the samples were taken into account.
3.3 Gut mycobiome alternations in IBD patients with CDI
We then assessed the gut fungal microbiota to identify the mycobiome alternations in IBD patients with CDI. The alpha diversity of fungi in the IBD-only and IBD-CDI groups was significantly higher than that in the HC group. The analysis of alpha diversity showed no difference between IBD-only and IBD-CDI groups (Figure 4A; Supplementary Table 5). To further explore the equilibrium between bacteria and fungi diversity in the gut, we calculated the fungi-to-bacteria diversity ratio. This ratio was significantly increased in patients with IBD-only or IBD-CDI compared with HC. Similarly, no distinct difference was observed between the IBD-only and IBD-CDI groups (Figure 4B). A significant difference of fungal beta diversity was observed between samples from HC and patients with IBD-only or IBD-CDI (Figure 4C). The results of lefse analysis showed several mycobiome alterations associated with IBD-CDI, including the enrichments of the order Saccharomycetales, the class Saccharomycetes, the families Saccharomycetaceae and Debaryomycetaceae, the genera Saccharomyces and Candida, and the species Saccharomyces cerevisiae (S. cerevisiae) and Candida albicans (C. albicans) (Figure 4D; Supplementary Table 5). As showed in Figures 4E, F, the relative abundance of Saccharomyces cerevisiae and Candida albicans increased in the sequence of HC, IBD-only, and IBD-CDI.
Figure 4 Gut mycobiome alternations among the HC, IBD-only, and IBD-CDI groups. (A) The fungal alpha diversity. (B) Fungi-to-bacteria alpha diversity ratio. (C) Fungal beta diversity analysis with the PCoA plots of the three study groups. (D) Lefse analysis of the mycobiome alterations. Taxa with LDA score >4.0 and P < 0.05 were regarded to be significantly discriminative. Only taxa present in >20% of the samples were taken into account. (E) Box plot of relative abundance of Saccharomyces cerevisiae. (F) Box plot of relative abundance of Candida albicans. ***P < 0.001. ****P < 0.0001.
3.4 Functional signatures of IBD patients with CDI
Functional analysis of the three study groups was conducted with HUMANn2. Principal component analysis (PCA) was used to demonstrate the functional composition. The HC samples clustered together, while the IBD-only and IBD-CDI samples were widely scattered (Figure 5A). A total of 275 differential pathways among the HC, IBD-only, and IBD-CDI groups were identified (Supplementary Table 5). In post-hoc analysis, 22 pathways with significantly different relative abundance between IBD-only and IBD-CDI were identified (Figure 5B; Supplementary Table 5). Four of these pathways were related to peptidoglycan metabolism: PWY-5265 (peptidoglycan biosynthesis II [staphylococci]), PWY-6471 (peptidoglycan biosynthesis IV [Enterococcus faecium]), PWY0-1586 (peptidoglycan maturation [meso-diaminopimelate containing]), and PWY-6470 (peptidoglycan biosynthesis V [beta-lactam resistance]). The mean relative abundance of PWY-5265 and PWY-6471 in patients with IBD-CDI were higher than those with IBD-only or HC. Patients with IBD-CDI had decreased abundance of PWY0-1586 and increased abundance of PWY-6470 compared with patients with IBD-only (Figure 5C). We further conducted stratified analysis of these four pathways to identify the metagenomically contributed species (Figure 5D; Supplementary Table 5). Enterococcus faecium was the dominant contributor in PWY-5265 in the IBD-CDI group. Additionally, Enterococcus faecium and Clostridioides difficile were the main contributors in PWY-6470 and PWY-6471 in these patients. PWY0-1586 were not dominated by any single species.
Figure 5 Functional analysis(A) Principal component analysis of the functional composition. (B) Fold changes of the significantly different pathways between IBD-CDI and IBD-only groups. Only pathways present in >20% of the samples were taken into account. Data was normalized with z score in the heatmap. (C) Bar chart of the four pathways involved with peptidoglycan metabolism. (D) Stratified analysis of the four pathways by species.
3.5 Correlation networks of the study groups
SparCC analysis was used to construct the correlation network of the three study groups (Friedman and Alm, 2012). Compared with HC and IBD-only, the network of IBD-CDI had more nodes and edges. But the HC network had higher average number of neighbors than the IBD-only and IBD-CDI networks (Figure 6; Supplementary Table 6). Several microbial markers were present in the network, such as Ruminococcus gnavus and Clostridium innocuum. Ruminococcus gnavus was present in both the IBD-only and IBD-CDI network. In the IBD-only network, it was positively correlated with Bacteroides ovatus, and was negatively correlated with Bacteroides thetaiotaomicron. In the IBD-CDI network, it was positively related with Clostridium innocuum and Intestinibacter bartlettii, and was negatively correlated with Lactobacillus salivarius (L. salivarius). Clostridium innocuum was presented in the center of the IBD-CDI network, and was also negatively correlated with Lactobacillus salivarius.
Figure 6 Correlation network analysis. (A) The HC network. (B) The IBD-only network. (C) The IBD-CDI network. Only species present in >20% of the samples were taken into account, and only the statistically significantly correlations were shown.
4 Discussion
CDI is a common complication of IBD. C. difficile coinfection is associated with longer hospital stay, escalated medication, and higher colectomy rates of IBD patients. Investigation of gut microecology in IBD patients with CDI may facilitate the understanding of the disease, and provide new perspectives for treatments. In this study, we showed that IBD-CDI patients had a pronounced microbial dysbiosis and specific alternations in the gut mycobiome and functional profile. Bacterial dysbiosis in gut of IBD patients with CDI, including decreased alpha diversity and altered composition, were consistent with previous reports (Sokol et al., 2018; Mahnic et al., 2022). In keeping with previous results, clustering of the samples in IBD-only and IBD-CDI indicated heterogeneity regardless of C. difficile coinfection status (Bushman et al., 2020). In addition to confirmation, new data about the bacterial microbiota alternations were also provided. We found a decreased proportion of Faecalibacterium and an increased proportion of Enterococcus in IBD patients with CDI compared with IBD-only. Faecalibacterium has been reported as a probiotic, which had anti-inflammatory potential via butyrate production and regulation of butyrate-mediated inflammatory pathways (Roychowdhury et al., 2018; Zhou et al., 2018; Xu et al., 2020; Li et al., 2022a). Depletion of Faecalibacterium has been observed in pediatric CD patients with CDI (Hellmann et al., 2020). In our study, relative abundance of Faecalibacterium was significantly decreased in IBD-only patients compared with HC, and this alternation was even enhanced when comparing IBD-CDI and HC, suggesting a more drastic micro-ecological changes with both medical complications.
Several studies have shown that gut mycobiome shifts existed between IBD patients and healthy population. For the first time, we report the mycobiome alternations in IBD patients with CDI. Gut mycobiome in IBD patients with CDI was characterized by an increased alpha diversity compared with healthy controls, and by certain taxonomic alternations, including a markedly increased abundance of S. cerevisiae. Chehoud et al. have reported decreased fungal alpha diversity in IBD patients (Chehoud et al., 2015). However, more recent studies reported no significant difference between IBD patients and healthy people in fungal diversity (Imai et al., 2019; Lemoinne et al., 2020; Li et al., 2022b). Furthermore, primary sclerosing cholangitis, another common comorbidity of IBD, was reported with elevated biodiversity compared with healthy people (Lemoinne et al., 2020). In accordance with previous studies, we found that the fungi-to-bacteria diversity ratio was decreased in IBD patients compared with HC (Sokol et al., 2017; Lemoinne et al., 2020), suggesting disturbed balance between fungal and bacterial flora. S. cerevisiae has long been correlated with IBD, since the anti-S. cerevisiae antibodies (ASCA) were specifically elevated in CD patients and might help in diagnosis and outcome prediction of CD patients (Main et al., 1988; Quinton et al., 1998; Forcione et al., 2004; Reese et al., 2006). S. cerevisiae was detected in the fecal samples of pediatric CD patients (Lewis et al., 2015), and its abundance was decreased in UC patients and active IBD patients compared with healthy controls (Sokol et al., 2017; Li et al., 2022b). Overall, current knowledge of gut mycobiome is limited, and the results remain inconsistent. Whether gut mycobiome contributes to pathogenesis and development of IBD-CDI remains to be explore.
In the present study, we demonstrated that E. faecium was one of the differential microbial markers in IBD patients with CDI, and was involved in the peptidoglycan metabolism pathways, which were enriched in this group of patients. Both E. faecium and C. difficile were nosocomial pathogens and shared similar risk factors (Gerding, 1997). As reported, vancomycin-resistant E. faecium was detected in around 17%~19% patients with CDI (Garbutt et al., 1999; Ozsoy and Ilki, 2017). Vancomycin-resistant E. faecium colonization was reported as a risk factor of C. difficile-associated disease in immunocompromised population (Trifilio et al., 2013). Glycopeptide-resistant E. faecium can produce modified peptidoglycan precursors, which block the target of the antimicrobial agents (Hendrickx et al., 2013). However, there is also studies indicating that some commensal strains of E. faecium can facilitate probiotic bacteria activity against C. difficile pathogenesis by secreting peptidoglycan hydrolase and enhancing host immunity (Zheng et al., 2016; Kim et al., 2019). Moreover, this peptidoglycan remodeling activity might also benefit patients receiving immunotherapy (Griffin et al., 2021). In addition to peptidoglycan metabolism, previous study has suggested that both E. faecium and C. difficile could affect the antibiotic resistant activities of the antimicrobial resistant (Kunishima et al., 2019).
The other two microbial markers in IBD patients with CDI, C.innocuum and R. gnavus, were also implied in the network analysis. C.innocuum is a gram-positive, spore-forming, anaerobic bacterium (Cherny et al., 2021). C.innocuum frequently coexist with C. difficile, a previous study indicated that C.innocuum might cause antibiotic-associated diarrhoea with clinical manifestation similar to CDI (Chia et al., 2018; Cherny et al., 2022). Besides, C.innocuum can lead to creeping fat in Crohn’s disease (Ha et al., 2020), and is associated with poor clinical outcomes in UC patients (Le et al., 2022). Another biomarker, R. gnavus, was positively correlated with C.innocuum in IBD-CDI network. R. gnavus can utilize the mucin and fucosylated glycans in the gastrointestinal tract, and it can also produce an inflammatory polysaccharide, which might induce secretion of tumor necrosis factor alpha by dendritic cells (Crost et al., 2013; Henke et al., 2019; Henke et al., 2021). In addition, high levels of R. gnavus were observed in the donors of UC patients who failed the FMT, confirming a proinflammatory role of R. gnavus in the gut of IBD patients (Fuentes et al., 2017). Notably, both C.innocuum and R. gnavus were negatively correlated with L. salivarius, a long recognized probiotics. Previous studies showed that probiotic strategies containing L. salivarius had beneficial effects on IBD patients (Palumbo et al., 2016; Ghavami et al., 2020). A recent study suggested that DL-endopeptidase-producing L. salivarius could promote the restoration of the pattern-recognition receptor NOD2 in mice with CD, and further exerted potent anti-colitis effects (Gao et al., 2022). Therefore, probiotic supplement with L. salivarius and regulation of the balance between C.innocuum, R. gnavus, and L. salivarius might be potential therapeutic options for IBD patients with CDI.
This study has several limitations. First, it is a single-centre study, with relatively small sample size, which may introduce bias about the composition and alternations of gut microflora. But our sample size is larger than former studies in IBD patients with CDI (Sokol et al., 2018; Bushman et al., 2020). Second, the study is cross-sectional, samples were only collected before the initiation of treatment. Therefore, this study cannot provide information about dynamics of gut microbiota. Longitudinal analysis is warranted in the future. Lastly, the study population were heterogenous, and the prior medication exposure varied among patients. Future analysis to address these limitations would be useful.
In summary, this study showed that IBD patients with CDI had pronounced microbial dysbiosis, with decreased biodiversity and certain bacterial and fungal alternations. Several functional changes, such as peptidoglycan metabolism, were enriched in this subset of patients. We also found signatures of interactions between microbial markers in these patients, which might provide new insights for microbiota-based treatment in IBD-CDI patients.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics statement
The studies involving human participants were reviewed and approved by The Ethics Committees of Peking Union Medical College Hospital. The patients/participants provided their written informed consent to participate in this study.
Author contributions
SY, XG, and HX contributed equally to this study. YL and HX conceived and designed the study. SY, XG, HX, BWT, YS, BT, and YD performed the experiments. SY, XG, and HX analysed the data and wrote the manuscript. YL, SH, and JQ provided general guidance and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by CAMS Innovation Fund for Medical Sciences (CIFMS) from Chinese Academy of Medical Sciences 2021-I2M-1-062; Beijing Municipal Natural Science Foundation (7212078); and National key clinical specialty construction project ZK108000.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2023.1129043/full#supplementary-material
Supplementary Table 1 | Clinical profile of the patients
Supplementary Table 2 | Species list of the study groups
Supplementary Table 3 | Taxonomy composition of the study groups
Supplementary Table 4 | Microbial markers with significant differences among the study groups
Supplementary Table 5 | Gut mycobiome analysis of the study groups
Supplementary Table 6 | Functional analysis of the study groups
Supplementary Table 7 | Network analysis of the study groups
References
Abubucker, S., Segata, N., Goll, J., Schubert, A. M., Izard, J., Cantarel, B. L., et al. (2012). Metabolic reconstruction for metagenomic data and its application to the human microbiome. PloS Comput. Biol. 8 (6), e1002358. doi: 10.1371/journal.pcbi.1002358
Beghini, F., McIver, L. J., Blanco-Miguez, A., Dubois, L., Asnicar, F., Maharjan, S., et al. (2021). Integrating taxonomic, functional, and strain-level profiling of diverse microbial communities with bioBakery 3. Elife 10, e65088. doi: 10.7554/eLife.65088
Berkell, M., Mysara, M., Xavier, B. B., van Werkhoven, C. H., Monsieurs, P., Lammens, C., et al. (2021). Microbiota-based markers predictive of development of clostridioides difficile infection. Nat. Commun. 12 (1), 2241. doi: 10.1038/s41467-021-22302-0
Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics 30 (15), 2114–2120. doi: 10.1093/bioinformatics/btu170
Bushman, F. D., Conrad, M., Ren, Y., Zhao, C., Gu, C., Petucci, C., et al. (2020). Multi-omic analysis of the interaction between clostridioides difficile infection and pediatric inflammatory bowel disease. Cell Host Microbe 28 (3), 422–433 e427. doi: 10.1016/j.chom.2020.07.020
Chehoud, C., Albenberg, L. G., Judge, C., Hoffmann, C., Grunberg, S., Bittinger, K., et al. (2015). Fungal signature in the gut microbiota of pediatric patients with inflammatory bowel disease. Inflammation Bowel Dis. 21 (8), 1948–1956. doi: 10.1097/MIB.0000000000000454
Cherny, K. E., Muscat, E. B., Balaji, A., Mukherjee, J., Ozer, E. A., Angarone, M. P., et al. (2022). Association between clostridium innocuum and antibiotic-associated diarrhea in adults and children: A cross-sectional study and comparative genomics analysis. Clin. Infect. Dis. 20, ciac483. doi: 10.1093/cid/ciac483
Cherny, K. E., Muscat, E. B., Reyna, M. E., Kociolek, L. K. (2021). Clostridium innocuum: Microbiological and clinical characteristics of a potential emerging pathogen. Anaerobe 71, 102418. doi: 10.1016/j.anaerobe.2021.102418
Chia, J. H., Wu, T. S., Wu, T. L., Chen, C. L., Chuang, C. H., Su, L. H., et al. (2018). Clostridium innocuum is a vancomycin-resistant pathogen that may cause antibiotic-associated diarrhoea. Clin. Microbiol. Infect. 24 (11), 1195–1199. doi: 10.1016/j.cmi.2018.02.015
Crost, E. H., Tailford, L. E., Le Gall, G., Fons, M., Henrissat, B., Juge, N. (2013). Utilisation of mucin glycans by the human gut symbiont ruminococcus gnavus is strain-dependent. PloS One 8 (10), e76341. doi: 10.1371/journal.pone.0076341
Czepiel, J., Drozdz, M., Pituch, H., Kuijper, E. J., Perucki, W., Mielimonka, A., et al. (2019). Clostridium difficile infection: Review. Eur. J. Clin. Microbiol. Infect. Dis. 38 (7), 1211–1221. doi: 10.1007/s10096-019-03539-6
D'Aoust, J., Battat, R., Bessissow, T. (2017). Management of inflammatory bowel disease with clostridium difficile infection. World J. Gastroenterol. 23 (27), 4986–5003. doi: 10.3748/wjg.v23.i27.4986
Forcione, D. G., Rosen, M. J., Kisiel, J. B., Sands, B. E. (2004). Anti-saccharomyces cerevisiae antibody (ASCA) positivity is associated with increased risk for early surgery in crohn's disease. Gut 53 (8), 1117–1122. doi: 10.1136/gut.2003.030734
Franzosa, E. A., Sirota-Madi, A., Avila-Pacheco, J., Fornelos, N., Haiser, H. J., Reinker, S., et al. (2019). Gut microbiome structure and metabolic activity in inflammatory bowel disease. Nat. Microbiol. 4 (2), 293–305. doi: 10.1038/s41564-018-0306-4
Frau, A., Ijaz, U. Z., Slater, R., Jonkers, D., Penders, J., Campbell, B. J., et al. (2021). Inter-kingdom relationships in crohn's disease explored using a multi-omics approach. Gut Microbes 13 (1), 1930871. doi: 10.1080/19490976.2021.1930871
Friedman, J., Alm, E. J. (2012). Inferring correlation networks from genomic survey data. PloS Comput. Biol. 8 (9), e1002687. doi: 10.1371/journal.pcbi.1002687
Fuentes, S., Rossen, N. G., van der Spek, M. J., Hartman, J. H., Huuskonen, L., Korpela, K., et al. (2017). Microbial shifts and signatures of long-term remission in ulcerative colitis after faecal microbiota transplantation. ISME J. 11 (8), 1877–1889. doi: 10.1038/ismej.2017.44
Gao, J., Zhao, X., Hu, S., Huang, Z., Hu, M., Jin, S., et al. (2022). Gut microbial DL-endopeptidase alleviates crohn's disease via the NOD2 pathway. Cell Host Microbe 30 (10), 1435–1449.e1439. doi: 10.1016/j.chom.2022.08.002
Garbutt, J. M., Littenberg, B., Evanoff, B. A., Sahm, D., Mundy, L. M. (1999). Enteric carriage of vancomycin-resistant enterococcus faecium in patients tested for clostridium difficile. Infect. Control Hosp. Epidemiol. 20 (10), 664–670. doi: 10.1086/501562
Gerding, D. N. (1997). Is there a relationship between vancomycin-resistant enterococcal infection and clostridium difficile infection? Clin. Infect. Dis. 25 (Suppl 2), S206–S210. doi: 10.1086/516247
Ghavami, S. B., Yadegar, A., Aghdaei, H. A., Sorrentino, D., Farmani, M., Mir, A. S., et al. (2020). Immunomodulation and generation of tolerogenic dendritic cells by probiotic bacteria in patients with inflammatory bowel disease. Int. J. Mol. Sci. 21 (17), 6266. doi: 10.3390/ijms21176266
Griffin, M. E., Espinosa, J., Becker, J. L., Luo, J. D., Carroll, T. S., Jha, J. K., et al. (2021). Enterococcus peptidoglycan remodeling promotes checkpoint inhibitor cancer immunotherapy. Science 373 (6558), 1040–1046. doi: 10.1126/science.abc9113
Guzzo, G. L., Andrews, J. M., Weyrich, L. S. (2022). The neglected gut microbiome: Fungi, Protozoa, and bacteriophages in inflammatory bowel disease. Inflammation Bowel Dis. 28 (7), 1112–1122. doi: 10.1093/ibd/izab343
Ha, C. W. Y., Martin, A., Sepich-Poore, G. D., Shi, B., Wang, Y., Gouin, K., et al. (2020). Translocation of viable gut microbiota to mesenteric adipose drives formation of creeping fat in humans. Cell 183 (3), 666–683.e617. doi: 10.1016/j.cell.2020.09.009
Hellmann, J., Andersen, H., Fei, L., Linn, A., Bezold, R., Lake, K., et al. (2020). Microbial shifts and shorter time to bowel resection surgery associated with c. difficile in pediatric crohn's disease. Inflammation Bowel Dis. 26 (8), 1212–1221. doi: 10.1093/ibd/izz263
Hendrickx, A. P., van Schaik, W., Willems, R. J. (2013). The cell wall architecture of enterococcus faecium: from resistance to pathogenesis. Future Microbiol. 8 (8), 993–1010. doi: 10.2217/fmb.13.66
Henke, M. T., Brown, E. M., Cassilly, C. D., Vlamakis, H., Xavier, R. J., Clardy, J. (2021). Capsular polysaccharide correlates with immune response to the human gut microbe ruminococcus gnavus. Proc. Natl. Acad. Sci. U S A. 118, (20). doi: 10.1073/pnas.2007595118
Henke, M. T., Kenny, D. J., Cassilly, C. D., Vlamakis, H., Xavier, R. J., Clardy, J. (2019). Ruminococcus gnavus, a member of the human gut microbiome associated with crohn's disease, produces an inflammatory polysaccharide. Proc. Natl. Acad. Sci. U S A. 116 (26), 12672–12677. doi: 10.1073/pnas.1904099116
Hvas, C. L., Dahl Jørgensen, S. M., Jørgensen, S. P., Storgaard, M., Lemming, L., Hansen, M. M., et al. (2019). Fecal microbiota transplantation is superior to fidaxomicin for treatment of recurrent clostridium difficile infection. Gastroenterology 156 (5), 1324–1332.e1323. doi: 10.1053/j.gastro.2018.12.019
Imai, T., Inoue, R., Kawada, Y., Morita, Y., Inatomi, O., Nishida, A., et al. (2019). Characterization of fungal dysbiosis in Japanese patients with inflammatory bowel disease. J. Gastroenterol. 54 (2), 149–159. doi: 10.1007/s00535-018-1530-7
Kim, B., Wang, Y. C., Hespen, C. W., Espinosa, J., Salje, J., Rangan, K. J., et al. (2019). Enterococcus faecium secreted antigen a generates muropeptides to enhance host immunity and limit bacterial pathogenesis. Elife 8, e45343. doi: 10.7554/eLife.45343
Kucharzik, T., Ellul, P., Greuter, T., Rahier, J. F., Verstockt, B., Abreu, C., et al. (2021). ECCO guidelines on the prevention, diagnosis, and management of infections in inflammatory bowel disease. J. Crohns Colitis. 15 (6), 879–913. doi: 10.1093/ecco-jcc/jjab052
Kunishima, H., Ishibashi, N., Wada, K., Oka, K., Takahashi, M., Yamasaki, Y., et al. (2019). The effect of gut microbiota and probiotic organisms on the properties of extended spectrum beta-lactamase producing and carbapenem resistant enterobacteriaceae including growth, beta-lactamase activity and gene transmissibility. J. Infect. Chemother. 25 (11), 894–900. doi: 10.1016/j.jiac.2019.04.021
Lamendella, R., Wright, J. R., Hackman, J., McLimans, C., Toole, D. R., Bernard Rubio, W., et al. (2018). Antibiotic treatments for clostridium difficile infection are associated with distinct bacterial and fungal community structures. mSphere 3 (1), e00572–17. doi: 10.1128/mSphere.00572-17
Le, P. H., Chiu, C. T., Yeh, P. J., Pan, Y. B., Chiu, C. H. (2022). Clostridium innocuum infection in hospitalised patients with inflammatory bowel disease. J. Infect. 84 (3), 337–342. doi: 10.1016/j.jinf.2021.12.031
Lemoinne, S., Kemgang, A., Ben Belkacem, K., Straube, M., Jegou, S., Corpechot, C., et al. (2020). Fungi participate in the dysbiosis of gut microbiota in patients with primary sclerosing cholangitis. Gut 69 (1), 92–102. doi: 10.1136/gutjnl-2018-317791
Lewis, J. D., Chen, E. Z., Baldassano, R. N., Otley, A. R., Griffiths, A. M., Lee, D., et al. (2015). Inflammation, antibiotics, and diet as environmental stressors of the gut microbiome in pediatric crohn's disease. Cell Host Microbe 18 (4), 489–500. doi: 10.1016/j.chom.2015.09.008
Li, X. V., Leonardi, I., Putzel, G. G., Semon, A., Fiers, W. D., Kusakabe, T., et al. (2022b). Immune regulation by fungal strain diversity in inflammatory bowel disease. Nature 603 (7902), 672–678. doi: 10.1038/s41586-022-04502-w
Li, H. B., Xu, M. L., Xu, X. D., Tang, Y. Y., Jiang, H. L., Li, L., et al. (2022a). Faecalibacterium prausnitzii attenuates CKD via butyrate-renal GPR43 axis. Circ. Res. 131 (9), e120–e134. doi: 10.1161/circresaha.122.320184
Li, Y., Xu, H., Xu, T., Xiao, M., Tang, H., Wu, D., et al. (2018). Case-control study of inflammatory bowel disease patients with and without clostridium difficile infection and poor outcomes in patients coinfected with c. difficile and cytomegalovirus. Dig. Dis. Sci. 63 (11), 3074–3083. doi: 10.1007/s10620-018-5230-1
Maccaferri, S., Biagi, E., Brigidi, P. (2011). Metagenomics: key to human gut microbiota. Dig. Dis. 29 (6), 525–530. doi: 10.1159/000332966
Mahnic, A., Pintar, S., Skok, P., Rupnik, M. (2022). Gut community alterations associated with clostridioides difficile colonization in hospitalized gastroenterological patients with or without inflammatory bowel disease. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.988426
Main, J., McKenzie, H., Yeaman, G. R., Kerr, M. A., Robson, D., Pennington, C. R., et al. (1988). Antibody to saccharomyces cerevisiae (bakers' yeast) in crohn's disease. Bmj 297 (6656), 1105–1106. doi: 10.1136/bmj.297.6656.1105
Ozsoy, S., Ilki, A. (2017). Detection of vancomycin-resistant enterococci (VRE) in stool specimens submitted for clostridium difficile toxin testing. Braz. J. Microbiol. 48 (3), 489–492. doi: 10.1016/j.bjm.2016.12.012
Palumbo, V. D., Romeo, M., Marino Gammazza, A., Carini, F., Damiani, P., Damiano, G., et al. (2016). The long-term effects of probiotics in the therapy of ulcerative colitis: A clinical study. BioMed. Pap. Med. Fac. Univ. Palacky Olomouc. Czech Repub. 160 (3), 372–377. doi: 10.5507/bp.2016.044
Parks, D. H., Tyson, G. W., Hugenholtz, P., Beiko, R. G. (2014). STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30 (21), 3123–3124. doi: 10.1093/bioinformatics/btu494
Pickard, J. M., Zeng, M. Y., Caruso, R., Núñez, G. (2017). Gut microbiota: Role in pathogen colonization, immune responses, and inflammatory disease. Immunol. Rev. 279 (1), 70–89. doi: 10.1111/imr.12567
Quinton, J. F., Sendid, B., Reumaux, D., Duthilleul, P., Cortot, A., Grandbastien, B., et al. (1998). Anti-saccharomyces cerevisiae mannan antibodies combined with antineutrophil cytoplasmic autoantibodies in inflammatory bowel disease: Prevalence and diagnostic role. Gut 42 (6), 788–791. doi: 10.1136/gut.42.6.788
Reese, G. E., Constantinides, V. A., Simillis, C., Darzi, A. W., Orchard, T. R., Fazio, V. W., et al. (2006). Diagnostic precision of anti-saccharomyces cerevisiae antibodies and perinuclear antineutrophil cytoplasmic antibodies in inflammatory bowel disease. Am. J. Gastroenterol. 101 (10), 2410–2422. doi: 10.1111/j.1572-0241.2006.00840.x
Rodríguez, C., Romero, E., Garrido-Sanchez, L., Alcaín-Martínez, G., Andrade, R. J., Taminiau, B., et al. (2020). Microbiota insights in clostridium difficile infection and inflammatory bowel disease. Gut Microbes 12 (1), 1725220. doi: 10.1080/19490976.2020.1725220
Roychowdhury, S., Cadnum, J., Glueck, B., Obrenovich, M., Donskey, C., Cresci, G. A. M. (2018). Faecalibacterium prausnitzii and a prebiotic protect intestinal health in a mouse model of antibiotic and clostridium difficile exposure. JPEN J. Parenter Enteral Nutr. 42 (7), 1156–1167. doi: 10.1002/jpen.1053
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12 (6), R60. doi: 10.1186/gb-2011-12-6-r60
Sokol, H., Jegou, S., McQuitty, C., Straub, M., Leducq, V., Landman, C., et al. (2018). Specificities of the intestinal microbiota in patients with inflammatory bowel disease and clostridium difficile infection. Gut Microbes 9 (1), 55–60. doi: 10.1080/19490976.2017.1361092
Sokol, H., Leducq, V., Aschard, H., Pham, H. P., Jegou, S., Landman, C., et al. (2017). Fungal microbiota dysbiosis in IBD. Gut 66 (6), 1039–1048. doi: 10.1136/gutjnl-2015-310746
Trifilio, S. M., Pi, J., Mehta, J. (2013). Changing epidemiology of clostridium difficile-associated disease during stem cell transplantation. Biol. Blood Marrow Transplant. 19 (3), 405–409. doi: 10.1016/j.bbmt.2012.10.030
Underhill, D. M., Braun, J. (2022). Fungal microbiome in inflammatory bowel disease: a critical assessment. J. Clin. Invest. 132, (5). doi: 10.1172/JCI155786
Ventin-Holmberg, R., Eberl, A., Saqib, S., Korpela, K., Virtanen, S., Sipponen, T., et al. (2021). Bacterial and fungal profiles as markers of infliximab drug response in inflammatory bowel disease. J. Crohns Colitis. 15 (6), 1019–1031. doi: 10.1093/ecco-jcc/jjaa252
Xu, J., Liang, R., Zhang, W., Tian, K., Li, J., Chen, X., et al. (2020). Faecalibacterium prausnitzii-derived microbial anti-inflammatory molecule regulates intestinal integrity in diabetes mellitus mice via modulating tight junction protein expression. J. Diab. 12 (3), 224–236. doi: 10.1111/1753-0407.12986
Zheng, A., Luo, J., Meng, K., Li, J., Bryden, W. L., Chang, W., et al. (2016). Probiotic (Enterococcus faecium) induced responses of the hepatic proteome improves metabolic efficiency of broiler chickens (Gallus gallus). BMC Genomics 17, 89. doi: 10.1186/s12864-016-2371-5
Keywords: inflammatory bowel disease, Clostridioides difficile infection, metagenomics, gut microbiome, mycobiome
Citation: Yu S, Ge X, Xu H, Tan B, Tian B, Shi Y, Dai Y, Li Y, Hu S and Qian J (2023) Gut microbiome and mycobiome in inflammatory bowel disease patients with Clostridioides difficile infection. Front. Cell. Infect. Microbiol. 13:1129043. doi: 10.3389/fcimb.2023.1129043
Received: 21 December 2022; Accepted: 24 January 2023;
Published: 06 February 2023.
Edited by:
Xiaoyan Wang, Department of Gastroenterology, Central South University, ChinaReviewed by:
Min Zhi, The Sixth Affiliated Hospital of Sun Yat-sen University, ChinaHui Yan, The University of Texas Health Science Center at San Antonio, United States
Copyright © 2023 Yu, Ge, Xu, Tan, Tian, Shi, Dai, Li, Hu and Qian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yue Li, liyue@pumch.ac.cn; Songnian Hu, husn@im.ac.cn
†These authors have contributed equally to this work