Roux-en-Y gastric bypass surgery changes fungal and bacterial microbiota in morbidly obese patients—A pilot study

Abstract

The Roux-en-Y gastric bypass (RYGB) remains the most effective treatment for morbidly obese patients to lower body weight and improve glycemic control. There is recent evidence that the mycobiome (fungal microbiome) can aggravate disease severity in a number of diseases including inflammatory bowel disease (IBD), irritable bowel syndrome (IBS) and hepatitis; moreover, a dysbiotic fungal microbiota has been reported in the obese. We characterized fungal and bacterial microbial composition in fecal samples of 16 morbidly obese patients before and three months after RYGB surgery and compared with nine healthy controls. We found that RYGB surgery induced a clear alteration in structure and composition of the gut fungal and bacterial microbiota. Beta diversity analysis revealed significant differences in bacterial microbiota between obese patients before surgery and healthy controls (P < 0.005) and a significant, unidirectional shift in RYGB patients after surgery (P < 0.001 vs. before surgery). In contrast, there was no significant difference in fungal microbiota between groups but individually specific changes after RYGB surgery. Interestingly, RYGB surgery induced a significant reduction in fungal alpha diversity namely Chao1, Sobs, and Shannon diversity index (P<0.05, respectively) which contrasts the trend for uniform changes in bacteria towards increased richness and diversity post-surgery. We did not observe any inter-kingdom relations in RYGB patients but in the healthy control cohort and there were several correlations between fungi and bacteria and clinical parameters (P<0.05, respectively) that warrant further research. Our study identifies changes in intestinal fungal communities in RYGB patients that are distinct to changes in the bacterial microbiota.

About the Authors:

Robert E. Steinert

Roles Conceptualization, Investigation, Project administration, Supervision, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Department of Surgery, University Hospital Zürich, Zürich, Switzerland

ORCID logo http://orcid.org/0000-0001-7884-5036

Ateequr Rehman

Roles Data curation, Formal analysis, Methodology

Affiliation: Institute of Clinical Molecular Biology, Christian-Albrechts-University of Kiel, Kiel, Germany

Everton Job Souto Lima

Roles Data curation, Formal analysis, Methodology, Software, Visualization

Affiliation: Microbiology and Systems Biology, The Netherlands Organization for Applied Scientific Research (TNO), Zeist, The Netherlands

Valeria Agamennone

Roles Data curation, Methodology, Project administration

Affiliation: Microbiology and Systems Biology, The Netherlands Organization for Applied Scientific Research (TNO), Zeist, The Netherlands

Frank H. J. Schuren

Roles Data curation, Methodology, Project administration, Validation

Affiliation: Microbiology and Systems Biology, The Netherlands Organization for Applied Scientific Research (TNO), Zeist, The Netherlands

Daniel Gero

Roles Investigation, Software, Visualization

Affiliation: Department of Surgery, University Hospital Zürich, Zürich, Switzerland

ORCID logo http://orcid.org/0000-0003-2941-9801

Phillip Schreiner

Roles Investigation

Affiliation: Department of Surgery, University Hospital Zürich, Zürich, Switzerland

René Vonlanthen

Roles Writing – review & editing

Affiliation: Department of Surgery, University Hospital Zürich, Zürich, Switzerland

Aiman Ismaeil

Roles Writing – review & editing

Affiliations Department of Surgery, University Hospital Zürich, Zürich, Switzerland, Department of Surgery, Aswan University, Tingar, Egypt

Stefanos Tzafos

Roles Writing – review & editing

Affiliation: Department of Surgery, University Hospital Zürich, Zürich, Switzerland

Hanna Hosa

Roles Writing – review & editing

Affiliation: Department of Surgery, University Hospital Zürich, Zürich, Switzerland

Diana Vetter

Roles Writing – review & editing

Affiliation: Department of Surgery, University Hospital Zürich, Zürich, Switzerland

Benjamin Misselwitz

Roles Writing – review & editing

Affiliation: Division of Visceral Surgery and Medicine, Inselspital Bern and Bern University, Bern, Switzerland

Marco Bueter

Roles Conceptualization, Data curation, Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing

Affiliation: Department of Surgery, University Hospital Zürich, Zürich, Switzerland

Introduction

The Roux-en-Y gastric bypass (RYGB) remains the most effective treatment for morbidly obese patients to lower body weight and improve glycemic control [1]. Several interrelated mechanism have been suggested for this including alterations in the postprandial secretion of gastrointestinal hormones [2, 3], increased delivery of bile acid salts to distal L-cells [2, 4], faster gastric transit and changes in vagal nerve signaling [5] and, more recently, changes in gut microbiota composition [6]. The gut microbiota is a dense and highly diverse microbial community composed of bacteria, fungi, phages, viruses, protozoa, and archaea that inhabits the gastrointestinal tract and is intimately involved in host physiology [7, 8]. An impairment of the gut microbiota, referred to as ‘dysbiosis’ and characterized by changes in its composition and function, is commonly observed in gastrointestinal and systemic diseases [9, 10]. In the obese, bacterial dysbiosis is characterized by compositional changes in specific bacterial groups and low bacterial gene richness and related functional pathways [11–13]. RYGB surgery restores bacterial richness and diversity and substantially alters the abundance of several bacterial groups [14–16]. It has been speculated that these alterations may support weight loss and metabolic improvements seen after surgery.

The fungal components of the gut microbiota have hardly been explored in health and disease. Fungi comprise less than 2% of total gut microorganisms, with only about 60 genera and 180 species [17]. A fungal cell is, however, about >100-fold larger than a bacterial cell and thus, represents substantially greater biomass than suggested by the number of available genomes [17]. Mycobiome dysbiosis has recently been implicated in the pathogenesis of inflammatory bowel diseases (IBD) [18], irritable bowel syndrome (IBS) [19] and hepatitis [20] and Rodríguez et al. [21] observed an increased abundance of Ascomycota, Saccharomycetes, Dipodascaceae and Saccharomycetaceae in the obese [21]. To date, it has not been explored whether RYGB surgery has an impact on fungal microbiota and whether potential changes are related to the metabolic improvements postoperatively.

Here we used an amplicon-based barcoding approach to characterize fungal and bacterial microbiota in fecal samples of morbidly obese patients before and three months after RYGB surgery and healthy, normal-weight controls. We applied a variety of analyses to assess differences in microbial community composition, possible inter-kingdom correlations, and relationships with clinical patient parameters.

Materials and methods

Patients and sample collection

Sixteen RYGB patients were recruited at the Department of Surgery, University Hospital Zurich, Switzerland. Each participant gave written, informed consent before participating in the study. All patients underwent complete evaluation before the operation and during follow-up, including medications, nutritional behavior, anthropometric and clinical parameters, and blood sampling for laboratory tests. RYGB patients received surgery between July and September 2017. Fecal samples were collected the day before surgery and three months after, using the deoxyribose nucleic acid (DNA) Genotek Omnigene Gut Collection kits (DNA Genotek, Ontario, Canada) following the manufacturer’s instructions. The kits were returned to a −80 °C freezer within 24 h of collection and stored until analysis.

A separate cohort of Dutch healthy normal weight subjects was used as a baseline. This approach was implied to see the direction of microbial changes after RYGB surgery towards, or against a healthy microbiome profile. Normal-weight controls were nine non-smokers, non-drug using individuals from the Netherlands who did not regularly experience intestinal problems and were not under treatment for conditions or diseases of the digestive tract. Healthy subjects had not taken antibiotics for at least seven months before sampling, nor had they used other medications in the previous two weeks. Each participant gave written, informed consent before participating in the study.

Subject characteristics are presented in Table 1. All research was carried out in compliance with the Helsinki Declaration, and for the Swiss RYGB cohort, it was approved by the ethics committee Zurich, Switzerland (BASEC Nr. 2012–0260) while for the healthy Dutch controls it was approved by the ethics committee Brabant, The Netherlands.

[Figure omitted. See PDF.]

Table 1. Subject characteristics.

https://doi.org/10.1371/journal.pone.0236936.t001

DNA isolation, amplicon PCR and illumina sequencing

DNA was extracted from fecal samples, as described previously [22]. DNA samples were then used for the 16S rRNA gene and Internal Transcribed Spacer (ITS) region sequencing to characterize bacterial and fungal communities (for details, see S1 File).

Microbiome (16S and ITS) data analysis

The 16S rRNA gene and ITS sequencing data were pre-processed, analyzed, and classified using Mothur (v. 1.36.1) [23]. Low-quality regions of the 16S rRNA gene sequences were trimmed using Btrim [24] with a sliding window size of 5 nucleotide and an average quality score of 25. For the ITS sequences, cutadapt was used for trimming [25] with a quality cutoff of 20. Paired-end reads were merged, and unique sequences were filtered based on length, with no ambiguous bases allowed. Chimeric sequences were identified and removed using UCHIME (Edgar 2011). Unique 16S rRNA gene sequences were aligned to the bacterial SILVA SEED reference database (Release 128, available at https://mothur.org/wiki/Silva_reference_files). Taxonomic assignment was performed using the Ribosomal Database Project (RDP) naïve Bayesian classifier with a confidence threshold of 60% and 1000 iterations [26] and the mothur-formatted version of the RDP training set v.14 (trainset14_032015). Unique ITS sequences were classified using the UNITE reference database v.6 (mothur release, available at https://unite.ut.ee/repository.php), with a confidence threshold of 60% and 1000 iterations. 16S rRNA gene and ITS sequence data have been submitted to the European Nucleotide Archive (accession no. PRJEB38917).

Statistical analysis

Statistical analysis was conducted using the open source software R, version 3.5.1 [27], and figures were generated using ggplot2 [28]. Descriptive statistics were used for demographic variables of the RYGB and heathy control cohort. To test for significant differences between groups, comparisons used Student’s independent and paired t-test. Differences in bacterial and fungal community composition between the sample groups were assessed at the phylum, class, family and genus level using the DESeq2 package (Love 2014). The Benjamini-Hochberg (BH) false discovery rate (FDR) correction was applied to the resulting p-values [29]. Microbial alpha diversity estimations are known to be influenced by sequencing depth. Thus, to rarify this uneven sequencing depths, subsampling (bacterial 11514 seqeunces per samples and mycobiome 98 sequences per sample) approach was implemented to calculate four diversity indices.: (i) Abundance-based Coverage Estimator (ACE) metric, (ii) Chao1, (iii) observed species, and (iv) Shannon’s diversity. Analysis of variance (ANOVA) was performed to test global differences in diversity among groups, in addition, paired t-tests was applied to compare RYGB patinets before and after surgery. For beta diversity analysis, due to uneven sampling depth (library size), the microbiome data was normalized using the variance stabilizing transformation (VST) method in the DESeq2 package. The transformed data was analyzed with PerMANOVA using the adonis function of the package vegan in R (http://cc.oulu.fi/~jarioksa/softhelp/vegan/html/adonis.html) to determine community composition variation across groups at different taxonomic levels. P-value adjustment was done via p.adjust function using the Benjamini & Hochberg method. Spearman correlation analysis was used to explore the relationship between changes in bacterial and fungal communities, and between these changes and clinical patient parameters. The Benjamin-Hochberg method was used for adjustment of p-values [29]. Ordination analysis was performed using Non-metric Multidimensional Scaling (NMDS) to investigate patterns in community composition.

Results

RYGB surgery changes mycobiome diversity

Bacterial and fungal alpha-diversity was comparable between in RYGB patients before surgery and the healthy control cohort (Fig 1A and 1B). After surgery, a trend for an increase in bacterial alpha-diversity was observed when compared to patients before surgery in measures of estimated alpha-diversity (ACE and Chao1). There was, however, no noticeable change in the number of observed bacterial species (Sobs) and Shannon diversity (Fig 1A). In contrast, for fungi we observed a significant reduction in estimated (Chao1, P = 0.03) and observed (Sobs, P = 0.01) species richness and Shannon diversity (P = 0.01, Fig 1B).

[Figure omitted. See PDF.]

Fig 1.

Alpha diversity of observed species of bacterial (A) and fungal (B) microbiota composition in patients before and after RYGB surgery and the healthy control cohort. Analysis of variance (ANOVA) was performed to test global differences in diversity among groups, in addition, paired t-tests was applied to compare RYGB patients before and after surgery. Non-metric Multi Dimensional Scaling (NMDS) plots were generated using Bray-Curtis distances (Beta diversity) for bacterial (C) and fungal (D) microbiota genera. Data were analyzed with PerMANOVA to determine community composition variation across groups at different taxonomic levels. P-value adjustment was done via p.adjust function using the Benjamini & Hochberg method.

https://doi.org/10.1371/journal.pone.0236936.g001

For bacterial composition, PerMANOVA suggested significant differences between RYGB patients before surgery and the healthy control cohort (P = 0.005). The effect of surgery on microbial composition and structure was very distinct and unidirectionally resulting in significant differences between patients before and after surgery (P = 0.001). Patients after surgery were, however, still distinct from the healthy control cohort (P = 0.001, Fig 1C). For fungal microbiota, NMDS analysis revealed no differences between RYGB patients before surgery and the control cohort. However, in patients after surgery, we observed a clear alteration in composition and structure when compared to before surgery that in contrast to the homogeneous shift in bacterial beta-diversity was specific to individuals and, therefore, not significantly different between groups (Fig 1D).

Bacterial microbiota composition changes after surgery

The microbiota of RYGB patients and normal-weight controls was dominated by Firmicutes, specifically by members of the Clostridia class. Several significant differences in normalized abundance at all taxonomical levels were observed between groups (Fig 2 and S1 Fig). At the phylum level, pre-surgery patients had significantly higher levels of Firmicutes and Actinobacteria and lower levels of Verrucomicrobia when compared to healthy controls (P ≤ 0.05, respectively). Post-surgery, abundances of these bacterial groups were not significantly different anymore, indicating a rebiosis of microbiota (Fig 2A and 2B). Proteobacteria levels were significantly different between all groups (P ≤ 0.05, respectively) with healthy controls showing the lowest levels and patients after surgery the highest (Fig 2A and 2B). At the genus level, patients before surgery had significantly higher levels of Blautia and Roseburia when compared with healthy controls (P ≤ 0.05, respectively) which decreased post-surgery in the direction of a healthy microbiome (Fig 2C and 2D). Abundances of Faecalibacterium and Bifidobacterium also decreased after surgery. However, this decrease was beyond abundances observed in healthy controls (Fig 2C and 2D). In line with observations at family and class level (S1 Fig), aerotolerant Streptococcus was significantly increased after surgery (P ≤ 0.05, respectively).

[Figure omitted. See PDF.]

Fig 2. Normalized mean abundance ± SEM (left) and relative abundance counts (right) for major bacteria at phylum (A and B) and genus (C and D) level in RYGB patients before and 3 month after surgery and in healthy controls.

Minor bacteria are grouped under ‘Other’. Differential abundance test was done by DESeq2 and considered to be significant at p <0.05, where (1) indicates differentially between control and after surgery, (2) between control and before surgery, and (3) between before surgery and after surgery.

https://doi.org/10.1371/journal.pone.0236936.g002

RYGB induces fungal microbiota changes at individual level

As mentioned, changes in fungal microbiota composition in response RYGB surgery were observed on an individual level. There was, however, no unidirectional change towards, or against, a healthy microbiome (Fig 1). This effect resulted in a lack of significant changes in normalized abundance at all taxonomical levels (Fig 3 and S2 Fig). Ascomycota and Basidiomycota were the most dominant fungal phyla in all subjects, and RYGB surgery induced a change of these fungal groups towards a healthy microbial profile (Fig 3A and 3B). Candida, Saccharomyces, and Pichia were among the most abundant fungi at the genus level with decreases in Candida and Saccharomyces and increases in Pichia after surgery (Fig 3C and 3D). When compared to healthy controls, changes in Saccharomyces were towards, and changes in Candida and Pichia, against, a healthy mycobiome. Cladosporium was less abundant in patients before surgery when compared with healthy controls. However, surgery increased its abundances towards a healthy mycobiome. In contrast, Cryptococcus was observed to be less abundant before surgery when compared with healthy controls, and this trend increased after surgery (Fig 3C and 3D).

[Figure omitted. See PDF.]

Fig 3. Normalized mean abundance ± SEM (left) and relative abundance counts (right) for major fungi at phylum (A and B) and genus (C and D) level in RYGB patients before and 3 month after surgery and in healthy controls.

Minor fungi are grouped under ‘Other’. There were no significant differences between subject groups.

https://doi.org/10.1371/journal.pone.0236936.g003

Associations between fungal and bacterial microbiota and clinical parameters

To further understand the interactions between bacterial and fungal groups in healthy subjects and RYGB patients, correlation analysis was performed in all three study groups separately. In healthy subjects, we found six positive correlations (Fig 4A) and three negative correlations while we failed to identify such inter-kingdom correlation in RYGB subjects before and after surgery.

[Figure omitted. See PDF.]

Fig 4. Heatmap of Spearman’s pairwise correlation coefficients between fungi and bacterial microbiota at genus level (A) and with metabolic parameters before and after surgery (B).

Only the top 20 abundant bacterial and top 8 fungi are shown. Blue squares indicate positive correlations, and red squares indicate inverse correlations. The shading of squares indicates the magnitude of the association; darker shades are more strongly associated than lighter shades. P < 0.05 significance denoted by *. BMI, body mass index; HbA1c, haemoglobin A1c; HDL, high-density-lipoprotein; LDL, low-density lipoprotein, CRP, C-reactive protein.

https://doi.org/10.1371/journal.pone.0236936.g004

We also identified several significant correlations between fungal and bacterial genera and clinical patient parameters before and after surgery (Fig 4B). For fungal genera before surgery, Saccharomyces and Pichia correlated negatively with fasting glucose and Cladosporium negatively with BMI. After surgery, we lost these correlations but identified a negative correlation between Pichia and triglycerides and a positive correlation between Debaryomyces and fasting glucose (Fig 4B). For bacterial genera before surgery, Blautia correlated positively with HDL cholesterol and negatively with BMI, while Streptococcus correlated negatively with triglycerides. After surgery, we were again unable to find any of those correlations. However, there were positive correlations between Roseburia and cholesterol (total and LDL), as well as between one unclassified Lachnospiraceae bacterium and Escherichia/Shigella and LDL cholesterol. Furthermore, Streptococcus correlated negatively with HDL cholesterol (Fig 4B).

Discussion

We characterized fungal and bacterial microbial composition in fecal samples of morbidly obese patients before and 3 months after RYGB surgery and compared with a healthy control cohort. We found clear but highly individualized changes in the fungal kingdom with a significant reduction in mycobiome richness and diversity in patients after surgery when compared to before surgery. This effect strikingly contrasts the well-described and uniform changes in bacteria towards increased richness and diversity post-surgery.

Consistent with many other studies, we found the gut mycobiome was dominated by fungi from the Ascomycota and Basidiomycota phyla [18, 19, 21] while Candida, Saccharomyces, and Cladosporium were the main genera [30–32]. Ascomycota appeared more prevalent in obese patients before surgery than after surgery and normal-weight controls while at genus level similar trends were observed for Saccharomyces and Debaryomyces. These observations are in line with a previous study comparing the fungal microbiome of 52 obese and non-obese patients [21], possibly indicating an effect of the obese phenotype on fungal microbiota composition. We also found a tendency for higher levels of Candida in obese patients before surgery than after surgery supporting a link of Candida with metabolic changes and intestinal inflammation as previously observed for diabetes [33] and IBD [34].

There was no significant difference in fungal alpha-diversity between obese patients before surgery and the healthy control cohort used in this study. Other studies in patients with IBD and chronic hepatitis showed increased alpha-diversity when compared with healthy subjects [20, 35, 36] while one recent study also observed reduced alpha-diversity in IBD patients [18]. Rodríguez et al. [21] observed no significant difference in fungal alpha-diversity between obese and non-obese subjects which is consistent with our findings. Interestingly, we found that RYGB surgery resulted in a significant reduction in observed and estimated fungal richness and diversity. This clearly contrasts the many reports of increased bacterial alpha-diversity after RYGB surgery [14–16, 37–39] although in the current study we were not able to replicate these earlier observations but only observed a trend for increased alpha-diversity measures (ACE and Chao1). Perhaps this is due to the limited number of subjects in our study (n = 16) when compared to earlier studies (n = 31 and 61) [14, 37], and/or different sample times after surgery. In fact, it is likely that the gut microbiota still changes between 3 month and 1 year after surgery to dynamically adapt to alterations in host metabolism (weight loss, fat mass loss, etc.). Thus, further investigation is warranted to better understand changes in bacterial and fungal microbiome over time and whether there are opposite changes of fungal and bacterial diversity in relation to host disease and pathogenesis.

With regard to the bacterial microbiota, we found the typical increase in the relative abundance of Proteobacteria after surgery [14–16, 37–39]. While due to their proinflammatory properties Proteobacteria are generally not considered to be beneficial [40, 41], there was no increase in CRP levels after surgery in our cohort (Table 1). This Proteobacteria increase might reflect the increase in oxygen availability in the large intestine after surgery, favoring facultative anaerobes such as Escherichia, Klebsiella, and Pseudomonas [42, 43]. Along the same lines, we observed a decrease in obligatory anaerobic gram-positive bacterial groups such as Blautia, Roseburia and Faecalibacterium from the phylum Firmicutes as well as Bifidobacterium from the phylum Actinobacteria [15, 16, 39]. It is possible that these compositional changes are also linked to (i) reduced gastric acid secretion after surgery [44], (ii) reduction in total energy intake with modifications in nutrient composition and (iii) increased availability of undigested nutrients in the large intestine due to the anatomic alteration after surgery.

Interestingly, RYGB surgery resulted in individualized changes in the fungal microbiota and not in unidirectional changes, as observed for bacterial microbiota. One explanation is that fungi are aerobic and thus less sensitive to the increase in oxygen availability after surgery. This may limit unidirectional changes as seen with bacteria that are highly sensitive to such gastrointestinal conditions. Moreover, while it is known that eating habits can influence the bacterial gut microbiota, there may be an even stronger influence of daily diet on the gut mycobiome. This is supported by a study by Hallen-Adams et al. [45], which analyzed the composition of the gut mycobiome at two time points each from 24 participants and found detection of the same fungus occurred in less than 20% of the time. Also David et al. [46] conducted a controlled diet study and found while bacterial composition showed a clear response to nutrient availability (carbohydrates/fiber vs. proteins and fats), fungal composition appeared to be driven by food colonization with the same species of fungi detected in participant fecal samples and in cheese fed to those participants. In fact, sequencing methods do not distinguish between DNA contributed from live or dead cells. Thus, it is not possible to differentiate fungi colonizing the gut from transients derived from diet and/or environment. Among the most dominant genera in our cohorts were Saccharomyces and Debaryomyces which are both commonly found in foods such as mushrooms, cheese, bread and beer.

We also performed correlation analyses to identify inter-kingdom relations and associations of bacterial and fungal microbiota with anthropometric and metabolic parameters. Several studies suggested a direct interaction of fungi and bacteria. In IBD, Sokol et al. [18] reported disease-specific patterns for inter-kingdom networks with increased strenghts of correlations between fungi and bacteria in ulcerative colitis but weaker correlations in Crohn´s disease. We observed several significant correlations between bacteria and fungi in healthy, normal-weight controls but not in RYGB patients, which may indicate persistence of a dysbiotic microbiota state in the obese before and after surgery. Moreover, we found that Pichia was negatively correlated with fasting glucose before surgery and with triglycerides after surgery, while Streptococcus was negatively correlated with triglycerides before surgery and with HDL after surgery. These metabolic and microbial changes might be independent results of the surgical intervention. However, it warrants further research to investigate if the correlated fungal and bacterial groups might also directly modulate metabolic parameters independent of surgery. Such direct effects could be tested in a RYGB rat model with a targeted fungicide/antibiotic intervention.

There are some limitations of our study. 1) In this pilot study, we only included 16 RYGB patients and 9 healthy control subject, thus, whether our data are representative of larger cohorts remains to be proven. 2) The healthy subject cohort which was used as a baseline to interpret the direction of changes in our main patient cohort was sampled at a distinct place using a different sampling method and was not adjusted for age and gender. Since geography including different eating habits, sampling methodology, age and gender can all impact microbiome measurements, caution is required when directly comparing these two cohorts. 3) Fungal and bacterial microbial composition was characterized only before and three months after surgery. While this allows to better understand the immediate effects of surgery, it does not provide information on the dynamic changes that may occur later in the postoperative phase as an adaption to changes in host metabolism such as weight loss, fat mass loss, etc. 4)) All current sequencing technologies only asses DNA signatures and thus cannot distinguish between DNA contributed from colonizing cells or fungi or transients derived from the diet. 5) Associations between specific fungi and 3-month clinical data such as cholesterol values need to be confirmed with a longer follow-up.

In conclusion, we found a clear alteration in structure and composition of fungal and bacterial microbiota in response to RYGB surgery. Beta diversity analysis revealed significant differences in bacterial microbiota between obese patients before surgery and healthy controls (P < 0.005) and a significant, unidirectional shift in RYGB patients after surgery (P < 0.001 vs. before surgery). In contrast, there was no significant difference in fungal microbiota between groups but individually specific changes after RYGB surgery. Interestingly, RYGB surgery induced a significant reduction in fungal alpha diversity namely Chao1, Sobs, and Shannon diversity index (P<0.05, respectively) which contrasts the trend for uniform changes in bacteria towards increased richness and diversity post-surgery. We did not observe any inter-kingdom relations in RYGB patients but in the healthy control cohort and there were several correlations between fungi and bacteria and clinical parameters (P<0.05, respectively) that warrant further research.

Supporting information

[Figure omitted. See PDF.]

S1 File. Supplementary material and methods.

https://doi.org/10.1371/journal.pone.0236936.s001

(DOCX)

S1 Fig. Normalized mean abundance +/- SEM (left) and relative abundance counts (right) for major bacteria at class (A and B) and family (C and D) level in RYGB patients before and 3 month after surgery and in healthy controls.

Minor bacteria are grouped under ‘Other’. Differential abundance test was done by DESeq2 and considered to be significant at p <0.05, where (1) indicates differentially between control and after surgery, (2) between control and before surgery, and (3) between before surgery and after surgery.

https://doi.org/10.1371/journal.pone.0236936.s002

(TIFF)

S2 Fig. Normalized mean abundance +/- SEM (left) and relative abundance counts (right) for major fungi at class (A and B) and family (C and D) level in RYGB patients before and 3 month after surgery and in healthy controls.

Minor fungi are grouped under ‘Other’. There were no significant differences between subject groups.

https://doi.org/10.1371/journal.pone.0236936.s003

(TIFF)

Citation: Steinert RE, Rehman A, Souto Lima EJ, Agamennone V, Schuren FHJ, Gero D, et al. (2020) Roux-en-Y gastric bypass surgery changes fungal and bacterial microbiota in morbidly obese patients—A pilot study. PLoS ONE 15(7): e0236936. https://doi.org/10.1371/journal.pone.0236936

References

1. Adams TD, Davidson LE, Litwin SE, Kim J, Kolotkin RL, Nanjee MN, et al. Weight and Metabolic Outcomes 12 Years after Gastric Bypass. N Engl J Med. 2017;377(12):1143–55. pmid:28930514

2. Steinert RE, Feinle-Bisset C, Asarian L, Horowitz M, Beglinger C, Geary N. Ghrelin, CCK, GLP-1, and PYY(3–36): Secretory Controls and Physiological Roles in Eating and Glycemia in Health, Obesity, and After RYGB. Physiol Rev. 2017;97(1):411–63. pmid:28003328

3. Peterli R, Steinert RE, Woelnerhanssen B, Peters T, Christoffel-Courtin C, Gass M, et al. Metabolic and hormonal changes after laparoscopic Roux-en-Y gastric bypass and sleeve gastrectomy: a randomized, prospective trial. Obesity surgery. 2012;22(5):740–8. pmid:22354457

4. Pournaras DJ, Glicksman C, Vincent RP, Kuganolipava S, Alaghband-Zadeh J, Mahon D, et al. The Role of Bile After Roux-en-Y Gastric Bypass in Promoting Weight Loss and Improving Glycaemic Control. Endocrinology. 2012;153(8):3613–9. pmid:22673227

5. Sinclair P, Brennan DJ, le Roux CW. Gut adaptation after metabolic surgery and its influences on the brain, liver and cancer. Nature reviews Gastroenterology & hepatology. 2018;15(10):606–24.

6. Furet JP, Kong LC, Tap J, Poitou C, Basdevant A, Bouillot JL, et al. Differential adaptation of human gut microbiota to bariatric surgery-induced weight loss: links with metabolic and low-grade inflammation markers. Diabetes. 2010;59(12):3049–57. pmid:20876719

7. Crommen S, Simon MC. Microbial Regulation of Glucose Metabolism and Insulin Resistance. Genes (Basel). 2017;9(1).

8. Clemente JC, Ursell LK, Parfrey LW, Knight R. The impact of the gut microbiota on human health: an integrative view. Cell. 2012;148(6):1258–70. pmid:22424233

9. Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490(7418):55–60. pmid:23023125

10. Joossens M, Huys G, Cnockaert M, De Preter V, Verbeke K, Rutgeerts P, et al. Dysbiosis of the faecal microbiota in patients with Crohn’s disease and their unaffected relatives. Gut. 2011;60(5):631–7. pmid:21209126

11. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480–4. pmid:19043404

12. 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. pmid:23985870

13. Aron-Wisnewsky J, Prifti E, Belda E, Ichou F, Kayser BD, Dao MC, et al. Major microbiota dysbiosis in severe obesity: fate after bariatric surgery. Gut. 2018.

14. Kong LC, Tap J, Aron-Wisnewsky J, Pelloux V, Basdevant A, Bouillot JL, et al. Gut microbiota after gastric bypass in human obesity: increased richness and associations of bacterial genera with adipose tissue genes. Am J Clin Nutr. 2013;98(1):16–24. pmid:23719559

15. Guo Y, Huang ZP, Liu CQ, Qi L, Sheng Y, Zou DJ. Modulation of the gut microbiome: a systematic review of the effect of bariatric surgery. European journal of endocrinology / European Federation of Endocrine Societies. 2018;178(1):43–56.

16. Magouliotis DE, Tasiopoulou VS, Sioka E, Chatedaki C, Zacharoulis D. Impact of Bariatric Surgery on Metabolic and Gut Microbiota Profile: a Systematic Review and Meta-analysis. Obesity surgery. 2017;27(5):1345–57. pmid:28265960

17. Mukherjee PK, Sendid B, Hoarau G, Colombel JF, Poulain D, Ghannoum MA. Mycobiota in gastrointestinal diseases. Nature reviews Gastroenterology & hepatology. 2015;12(2):77–87.

18. Sokol H, Leducq V, Aschard H, Pham HP, Jegou S, Landman C, et al. Fungal microbiota dysbiosis in IBD. Gut. 2017;66(6):1039–48. pmid:26843508

19. Botschuijver S, Roeselers G, Levin E, Jonkers DM, Welting O, Heinsbroek SEM, et al. Intestinal Fungal Dysbiosis Is Associated With Visceral Hypersensitivity in Patients With Irritable Bowel Syndrome and Rats. Gastroenterology. 2017;153(4):1026–39. pmid:28624575

20. Chen Y, Chen Z, Guo R, Chen N, Lu H, Huang S, et al. Correlation between gastrointestinal fungi and varying degrees of chronic hepatitis B virus infection. Diagn Microbiol Infect Dis. 2011;70(4):492–8. pmid:20846815

21. Mar Rodriguez M, Perez D, Javier Chaves F, Esteve E, Marin-Garcia P, Xifra G, et al. Obesity changes the human gut mycobiome. Scientific reports. 2015;5:14600. pmid:26455903

22. Ladirat SE, Schols HA, Nauta A, Schoterman MH, Keijser BJ, Montijn RC, et al. High-throughput analysis of the impact of antibiotics on the human intestinal microbiota composition. J Microbiol Methods. 2013;92(3):387–97. pmid:23266580

23. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41. pmid:19801464

24. Kong Y. Btrim: a fast, lightweight adapter and quality trimming program for next-generation sequencing technologies. Genomics. 2011;98(2):152–3. pmid:21651976

25. Marcel M. Cutadapt Removes Adapter Sequences From High-Throughput Sequencing Reads. EMBnetjournal. 2011;17(1):10–2.

26. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7. pmid:17586664

27. R Core Team (2018). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/.

28. Wickham H. Elegant Graphics for Data Analysis. Springer-Verlag New York. 2016.

29. Yoav BHY. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological). 1995;57(1).

30. Huseyin CE, Rubio RC, O’Sullivan O, Cotter PD, Scanlan PD. The Fungal Frontier: A Comparative Analysis of Methods Used in the Study of the Human Gut Mycobiome. Front Microbiol. 2017;8:1432. pmid:28824566

31. Gouba N, Raoult D, Drancourt M. Plant and fungal diversity in gut microbiota as revealed by molecular and culture investigations. PLoS One. 2013;8(3):e59474. pmid:23555039

32. Strati F, Di Paola M, Stefanini I, Albanese D, Rizzetto L, Lionetti P, et al. Age and Gender Affect the Composition of Fungal Population of the Human Gastrointestinal Tract. Front Microbiol. 2016;7:1227. pmid:27536299

33. Gosiewski T, Salamon D, Szopa M, Sroka A, Malecki MT, Bulanda M. Quantitative evaluation of fungi of the genus Candida in the feces of adult patients with type 1 and 2 diabetes—a pilot study. Gut Pathog. 2014;6(1):43. pmid:25328543

34. Standaert-Vitse A, Sendid B, Joossens M, Francois N, Vandewalle-El Khoury P, Branche J, et al. Candida albicans colonization and ASCA in familial Crohn’s disease. The American journal of gastroenterology. 2009;104(7):1745–53. pmid:19471251

35. Ott SJ, Kuhbacher T, Musfeldt M, Rosenstiel P, Hellmig S, Rehman A, et al. Fungi and inflammatory bowel diseases: Alterations of composition and diversity. Scandinavian journal of gastroenterology. 2008;43(7):831–41. pmid:18584522

36. Li Q, Wang C, Tang C, He Q, Li N, Li J. Dysbiosis of gut fungal microbiota is associated with mucosal inflammation in Crohn’s disease. J Clin Gastroenterol. 2014;48(6):513–23. pmid:24275714

37. Aron-Wisnewsky J, Prifti E, Belda E, Ichou F, Kayser BD, Dao MC, et al. Major microbiota dysbiosis in severe obesity: fate after bariatric surgery. Gut. 2019;68(1):70–82. pmid:29899081

38. Tremaroli V, Karlsson F, Werling M, Stahlman M, Kovatcheva-Datchary P, Olbers T, et al. Roux-en-Y Gastric Bypass and Vertical Banded Gastroplasty Induce Long-Term Changes on the Human Gut Microbiome Contributing to Fat Mass Regulation. Cell Metab. 2015;22(2):228–38. pmid:26244932

39. Ejtahed HS, Angoorani P, Hasani-Ranjbar S, Siadat SD, Ghasemi N, Larijani B, et al. Adaptation of human gut microbiota to bariatric surgeries in morbidly obese patients: A systematic review. Microb Pathog. 2018;116:13–21. pmid:29306011

40. Cani PD, Amar J, Iglesias MA, Poggi M, Knauf C, Bastelica D, et al. Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes. 2007;56(7):1761–72. pmid:17456850

41. Mukhopadhya I, Hansen R, El-Omar EM, Hold GL. IBD-what role do Proteobacteria play? Nature reviews Gastroenterology & hepatology. 2012;9(4):219–30.

42. Celiker H. A new proposed mechanism of action for gastric bypass surgery: Air hypothesis. Med Hypotheses. 2017;107:81–9. pmid:28915970

43. Hartman AL, Lough DM, Barupal DK, Fiehn O, Fishbein T, Zasloff M, et al. Human gut microbiome adopts an alternative state following small bowel transplantation. Proc Natl Acad Sci U S A. 2009;106(40):17187–92. pmid:19805153

44. Mason EE, Munns JR, Kealey GP, Wangler R, Clarke WR, Cheng HF, et al. Effect of gastric bypass on gastric secretion. Am J Surg. 1976;131(2):162–8. pmid:1251957

45. Hallen-Adams HE, Kachman SD, Kim J, Legge RM, Martínez I. Fungi inhabiting the healthy human gastrointestinal tract: a diverse and dynamic community. Fungal Ecology. 2015;15:9–17.

46. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505(7484):559–63. pmid:24336217