Ability of bifidobacteria to metabolize chitin-glucan and its impact on the gut microbiota

Scientific Reports, Apr 2019

Chitin-glucan (CG) represents a natural carbohydrate source for certain microbial inhabitants of the human gut and may act as a prebiotic for a number of bacterial taxa. However, the bifidogenic activity of this substrate is still unknown. In the current study, we evaluated the ability of chitin-glucan to influence growth of 100 bifidobacterial strains belonging to those species commonly identified within the bifidobacterial communities residing in the infant and adult human gut. Such analyses were coupled with transcriptome experiments directed to explore the transcriptional effects of CG on Bifidobacterium breve 2L, which was shown to elicit the highest growth performance on this natural polysaccharide. In addition, an in vivo trial involving a rat model revealed how the colonization efficiency of this bifidobacterial strain was enhanced when the animals were fed with a diet containing CG. Altogether our analyses indicate that CG is a valuable novel prebiotic compound that may be added to the human diet in order to re-establish/reinforce bifidobacteria colonization in the mammalian gut.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://www.nature.com/articles/s41598-019-42257-z.pdf

Ability of bifidobacteria to metabolize chitin-glucan and its impact on the gut microbiota

Ability of bifidobacteria to metabolize chitin-glucan and its impact on the gut microbiota Giulia Alessandri Christian Milani Sabrina Duranti Leonardo Mancabelli Thibaut Ranjanoro Salvatore Modica Luca Carnevali Rosario statello Francesca Bottacini Francesca t urroni Maria Cristina ossiprandi Andrea sgoi fo Douwe van si nderen Marco Ventura Published: xx xx xxxx Chitin-glucan (CG) represents a natural carbohydrate source for certain microbial inhabitants of the human gut and may act as a prebiotic for a number of bacterial taxa. However, the bifidogenic activity of this substrate is still unknown. In the current study, we evaluated the ability of chitin-glucan to influence growth of 100 bifidobacterial strains belonging to those species commonly identified within the bifidobacterial communities residing in the infant and adult human gut. Such analyses were coupled with transcriptome experiments directed to explore the transcriptional effects of CG on Bifidobacterium breve 2L, which was shown to elicit the highest growth performance on this natural polysaccharide. In addition, an in vivo trial involving a rat model revealed how the colonization efficiency of this bifidobacterial strain was enhanced when the animals were fed with a diet containing CG. Altogether our analyses indicate that CG is a valuable novel prebiotic compound that may be added to the human diet in order to re-establish/reinforce bifidobacteria colonization in the mammalian gut. - Chitin-glucan (CG) is a high-purity bio-polymer composed of two distinct polysaccharides represented by chitin (?-1,4-poly-N-acetyl-D-glucosamine) and ?-1,3-D-glucan in a ratio ranging from 20:80 to 40:60 (w/w)1. CG is typically extracted from Aspergillus niger mycelial cell walls. Furthermore, another natural sources of CG are represented by different fungal species and yeasts due to the presence of this bio-polymer in their inner cell walls2,3. Recently, CG has been indicated as food supplements fixing the maximum consumption rate to five gr per day for an average person4. However, mammalian enzymes are unable to degrade CG and therefore following ingestion this glycopolymer will arrive in its intact form in the large intestine where it may influence growth and/or metabolic activity of different members of the gut microbiota. In this context, there is growing scientific evidence of possible prebiotic effects elicited by CG towards various microorganisms of the mammalian gut5,6. Recently, an in vitro study performed in a dynamic gut simulator (SHIME) illustrated the effect of CG on microbiota composition and activity, leading to a decrease of the Firmicutes/Bacteroidetes ratio and an increase in the relative abundance of Roseburia spp1, which could suggest a potential role exploited by CG in shaping the gut microbiota through a prebiotic effect as previously displayed for other complex carbohydrates7. The term ?prebiotic? includes compounds, such as non-digestible (i.e. cannot be metabolized or degraded by the host) carbohydrates, that are selectively metabolized by beneficial gut bacteria. Prebiotic treatment is a dietary strategy by which the gastrointestinal microbiota can be modified, both in composition and/or activity, for the purpose of conferring health benefits to the host8. In fact, it has been demonstrated that prebiotics can reduce symptoms associated with inflammatory bowel disease, the duration of infectious and antibiotic-associated diarrhea, and the risk of cardiovascular diseases. Furthermore, prebiotics may also promote satiety and weight loss, preventing obesity and enhancing protective effects against the onset of colon cancer9. 1Department of Veterinary Medical Science, University of Parma, Parma, italy. 2Laboratory of Probiogenomics, Department of chemistry, Life Sciences, and environmental Sustainability, University of Parma, Parma, i taly. 3KitoZyme, Herstal, Belgium. 4Stress Physiology Laboratory, Department of chemistry, Life Sciences and environmental Sustainability, University of Parma, Parma, italy. 5APc Microbiome institute and School of Microbiology, Bioscience institute, national University of ireland, cork, i reland. Giulia Alessandri, christian Milani and Sabrina Duranti contributed equally. correspondence and requests for materials should be addressed to M.V. (email: ) Bifidobacteria are very prevalent and abundant human gut microbiota members, especially during the first months following birth10, though their numbers decrease following weaning and in elderly. Bifidobacterial abundance in the human microbiota is markedly reduced following gastrointestinal diseases11, suggesting that this taxon plays a positive role in the promotion of host health12. Prebiotics, such as dietary fibers, have been used to counteract the reduction of bifidobacterial abundance in the human gut13. So far, a very limited scientific evidence is available regarding the possible prebiotic effect of CG towards members of the Bifidobacterium genus6. In the current report, we evaluated possible bifidogenic features of CG towards 100 different bifidobacterial strains using an in vitro approach. Such analyses revealed the presence of a highly CG-responsive bifidobacterial strain, Bifidobacterium breve 2L, whose transcriptome when cultivated on CG was investigated in detail by RNAseq experiments. Furthermore, in vivo trials in a rat model fed with CG and B. breve 2L clearly demonstrated that CG increases the abundance of B. breve 2L and several other members of the mammalian gut microbiota. Results and Discussion Evaluation of prebiotic capability of chitin-glucan toward bifidobacteria. One hundred bifidobacterial strains previously isolated from the human environment (reported in Table?1), were evaluated for their ability to grow on CG, as unique carbon source. These strains were isolated from infant feces (Bifidobacterium bifidum, Bifidobacterium breve and Bifidobacterium pseudocatenulatum) as well as from gastrointestinal tract of adult (Bifidobacterium adolescentis, Bifidobacterium angulatum, Bifidobacterium catenulatum, Bifidobacterium dentium, Bifidobacterium longum, Bifidobacterium pseudolongum subsp. globosum). In addition, the bifidobacterial collection used in this project included strains belonging to Bifidobacterium animalis subsp. lactis, isolated from various commercial products sold as probiotics. Despite the fact that CG consists of two different polysaccharides, i.e. chitin (?-1,4-poly-N-acetyl-D-glucosamine) and ?-1,3-D-glucan, we focused our study on the biopolymer CG due to its commercial relevance as an extract of Aspergillus niger. Furthermore, we have used the Tyndallization procedure, which operates a series moderate heat treatment steps separated by rapid cooling stages, in order to achieve a sterilization of CG-based medium without causing extensive hydrolysis of CG as might occur when using high-level temperature treatments such as autoclaving. Growth assays were conducted on a modified MRS growth media, i.e. MRS w/o glu + CG (0.5% (w/v) together with a positive control (complete MRS) as well as a negative control (MRS w/o glu). Experiments were performed in triplicates for each growth media. We observed a minimal growth also on MRS w/o glu (?103 cells/ml), which is probably linked to the presence of components such as meat and yeast extract providing extra organic carbon sources. However, the final cell count of positive control and MRS w/o glu + CG were calculated by subtracting the cells count obtained on MRS w/o glu. Interestingly, all strains were shown to exhibit growth using MRS w/o glu + CG as the unique carbon source. Remarkably, almost all assayed strains displayed growth levels comparable with those obtained in complete MRS (MRS + glucose), i.e. comprised between 107?109 cells/mL (Table?1). Notably, this may be caused by partial degradation of CG in simpler sugars, e.g. glucose monomers from glucan, due to the heat treatment included in the modified Tyndallization approach employed to sterilize the MRS w/o glu + CG medium. Only in case of strains belonging to the B. adolescentis and B. longum species, the observed growth levels were significantly different (p-value < 0.05) (Fig.?1). In detail, B. adolescentis strains exhibited a final cell count of 7.59 ? 107 cells/ mL, and 4.19 ? 108 cells/mL, in MRS w/o glu + CG and in complete MRS, respectively. Regarding B. longum, these strains were able to reach an average final cell count of 4.40 ? 108 cells/mL in MRS w/o glu + CG, which is a lower cell yield compared to the final cell count in complete MRS (1.05 ? 109 cells/mL). This finding indicates that these bifidobacterial species have a reduced ability to utilize CG as compared to the other assessed bifidobacterial strains. In contrast, B. breve and B. bifidum strains were shown to reach the highest cell numbers on this substrate (cell numbers ranging from 108 to 109 cells/mL), and in some cases surpassing growth yields obtained in complete MRS (Fig.?1). These two bacterial species are typically isolated from infants, and consequently expected to be metabolically adapted to degrade host-specific glycans such as mucin and host glycan constituents like N-acetylglucosamine14,15. Notably, CG is partially composed of polymerized ?-1,4-poly-N-acetyl-D-glucosamine, which may explain the reason for the more robust metabolic efficiency of these species toward CG (compared to the other examined bifidobacterial species). For those bifidobacterial strains showing an identical or higher growth level in MRS w/o glu + CG vs. MRS, growth performances were further evaluated using a different approach, i.e. viable cell count on MRS Agar (Table?1). Moreover, three bifidobacterial strains, i.e., B. bifidum LMG 11041, B. breve 31L and B. dentium LMG 11405, were selected as controls, since these strains showed reduced growth levels in MRS w/o glu + CG when compared to the cell numbers reached when grown on MRS. Interestingly, all counts obtained through the plating method were lower than those obtained employing the Thoma cell counting chamber. This may not be surprising, since the Thoma chamber utilizes an indirect counting method and is unable to distinguish between living or devitalized cells, while the plating method exclusively evaluates viable cells. Notably, all strains grown using CG as the sole carbon source and then plated on MRS Agar were shown to be able to utilize the substrate tested. Bacterial counts ranged from 2.09 ? 104 CFU/mL for B. breve 31L cultures to 1.33 ? 108 CFU/mL for B. breve 2L cultures. Thus, B. breve 2L was the strain showing the highest growth performance when CG was provided as the sole carbon source. For this reason, B. breve 2L was chosen as a model bifidobacterial strain to dissect the genetic repertoire responsible for efficient CG metabolism. Identification of genes induced by CG in the genome of B. breve 2L. In order to identify the genes of B. breve 2L responsible for CG metabolism, we evaluated the transcriptomes of this strain when cultivated on CG by RNAseq analyses. The average transcriptome profile corresponding to the three replicates of B. breve 2L cultivated on MRS w/o glu + CG (0.5%(w/v) as the unique carbon source was compared with the average transcriptome profile of the three replicates of the positive control (B. breve 2L grown on complete MRS). Sequencing Strainsa B. adolescentis 1901B B. adolescentis 1902B B. adolescentis 1903B B. adolescentis 1904B B. adolescentis 22L B. adolescentis 236B B. adolescentis 382B B. adolescentis 388B B. adolescentis 42B B. adolescentis 487B B. adolescentis 532B B. adolescentis 548B B. adolescentis 59B B. adolescentis 61B B. adolescentis 65B B. adolescentis 679B B. adolescentis 703B B. adolescentis 70B B. adolescentis 711B B. adolescentis 723B B. adolescentis 731B B. adolescentis 734B B. adolescentis ATCC 15703 B. adolescentis LMG 10733 B. adolescentis LMG 10734 B. adolescentis LMG 18897 B. angulatum LMG 11039 B. animalis subsp. lactis Bb-12 B. animalis subsp. lactis DSM 10140 B. bifidum 156B B. bifidum 324B B. bifidum 361B B. bifidum 85B B. bifidum LMG 11041 B. bifidum LMG 11582 B. bifidum LMG 11583 B. bifidum LMG 13195 B. bifidum LMG 13200 B. bifidum PRL2010 B. breve 12L B. breve 2L B. breve 31L B. breve 687B B. breve 689B B. breve 691B B. breve LMG 13208 B. catenulatum 1231B B. catenulatum 1232B B. catenulatum 1233B B. catenulatum 1234B B. catenulatum LMG 11043 B. dentium 125B B. dentium 181B B. dentium 183B B. dentium 369B Continued Average Agar plate (CFU/mL) Agar plate SD Strainsa Origin B. dentium LMG 11405 Oral cavity B. gallicum LMG 11596 Colonoscopic sample B. longum 123B Colonoscopic sample B. longum 134B Colonoscopic sample B. longum 159B Colonoscopic sample B. longum 207B Colonoscopic sample B. longum 220B Colonoscopic sample B. longum 224B Colonoscopic sample B. longum 229B Colonoscopic sample B. longum 296B Adult stool sample B. longum 314B Colonoscopic sample B. longum 319B Colonoscopic sample B. longum 340B Adult stool sample B. longum 346B Adult stool sample B. longum 350B Adult stool sample B. longum 351B Adult stool sample B. longum 397B Adult stool sample B. longum 419B Adult stool sample B. longum 428B Adult stool sample B. longum 432B Adult stool sample B. longum 433B Adult stool sample B. longum 434B Adult stool sample B. longum 442B Adult stool sample B. longum 447B Adult stool sample B. longum 451B Adult stool sample B. longum 499B Colonoscopic sample B. longum 553B Colonoscopic sample B. longum 606B Colonoscopic sample B. longum 633B Adult stool sample B. longum 707B Adult stool sample B. longum 71B Adult stool sample B. longum 743B Adult stool sample B. longum 861B Adult stool sample B. longum 908B Colonoscopic sample B. longum subsp. longum LMG 13197 Colonoscopic sample B. pseudocatenulatum 202B Colonoscopic sample B. pseudocatenulatum 263B Adult stool sample B. pseudocatenulatum 289B Adult stool sample B. pseudocatenulatum 318B Colonoscopic sample B. pseudocatenulatum LMG 10505 Infant stool sample B. pseudolongum subsp. globosum 555B Colonoscopic sample B. pseudolongum subsp. globosum 685B Adult stool sample B. pseudolongum subsp. globosum 686B Colonoscopic sample LBM.pGseu1d1o5l9o6ngum subsp. globosum Bovine rumen B. stercoris JCM 15918 Adult stool sample reads of B. breve 2L grown on CG as well as on glucose were mapped on the genome sequence of B. breve 2L. Subsequently, evaluation of RPKM (reads per kilobase per million mapped reads) values for each gene revealed high expression (ranging from 1464 to 24685 RPKM) of genes predicted to encode carbohydrate transporters (Table?S1, annotation highlighted in red). Notably, B. breve 2L grown in MRS w/o glu + CG (0.5% (w/v) showed increased transcription, ranging from 8.1 to 106-fold, of various genes, some of which were predicted to be involved in carbohydrate internalization (Table?S2, annotation highlighted in red). Furthermore, significant transcriptional up-regulation, ranging from 8.2-fold to 14.6-fold, was observed for three genes encoding enzymes putatively involved in the hydrolysis of glycosidic linkages between hexose sugars (Table?S2, annotation highlighted in blue). A detailed scrutiny of the upregulated B. breve 2L genes upon its cultivation on MRS w/o glu + GC, revealed that the majority of these genes were organized in seven loci, which are predicted to encode carbohydrate transport systems (Table?S3, annotation highlighted in red), or are coding for putative glycosyl hydrolases (Table?S3, annotation highlighted in blue). Altogether, these results indicate that the inclusion of CG in the growth medium as the unique carbohydrate source, modulates the expression of genes encoding enzymes toward the degradation and metabolism of CG. Notably, CG was also observed to increase transcription of tad genes (Table?S3, annotation highlighted in green) responsible for the synthesis and assembly of the Type IVb pilus locus, which has been shown in another B. breve strain to mediate the colonization and persistence of bifidobacterial cells in the mammalian gut16?18. Evaluation of the colonization of B. breve 2L in rats following GC treatment. In order to evaluate if CG is able to modulate colonization of B. breve 2L in the mammalian gut, we performed an in vivo trial using Groningen rats (Rattus norvegicus). The trial consisted of three groups of animals, one receiving a daily inoculum of approximately 109 CFU of B. breve 2L, i.e., Breve2L, a second group fed with standard diet supplemented with 10% CG, i.e., CG, and another one treated with the same amount of B. breve 2L cells plus 10% CG, i.e., CG + Breve2L (Fig.?2). The putative physiological effect of CG on the Body Weight (BW) of the animals was recorded and compared to the Food Intake (FI) (Fig.?2b,c). Two-way ANOVA for repeated measures provided a significant effect of CG over time for BW changes (p-value < 0.05) and FI/BW (p-value < 0.01). BW increment was significantly higher in the Breve2L group as compared to the CG + Breve2L group at the end of the intervention period (T2) (p-value < 0.05). FI/BW ratio was significantly higher in the Breve2L group compared to either the CG (p-value < 0.05) or CG + Breve2L groups (p-value = 0.01), after the first week of the intervention period (T1). These data suggest that CG stimulates appetite while reducing/limiting body weight increment. In this context, characterization of the modulatory effects of CG toward the mammalian gut microbiota composition is pivotal to understand the physiological and metabolic aspects responsible for the reduction of body weight. In order to assess the number of B. breve 2L in the fecal samples of the animals enrolled in this study, we applied a qPCR approach (Fig.?3). These analyses highlighted that rats of the CG groups did not show a statistically significant increment of B. breve 2L at the end of the experiment (T3) compared to the other time points. In contrast, a statistically significant increased load of B. breve 2L (p-value < 0.05) was observed at time points T1 (1.14E + 03 CFU/gr) and T2 (1.20E + 03 CFU/gr) compared to T0 (3.30E + 02 CFU/gr) in the Breve2L group, but there were no significant differences with respect to T3 (9.58E+ 02 CFU/gr). These data suggest that the daily administration of B. breve 2L allows a transient colonization of the bacterial species as its concentration decreases after the washout week (T3). In this context, when rats were fed with both CG and B. breve 2L, the abundance of the B. breve strain was higher in T2 (4.37E + 03) and T3 (2.85E + 03) as compared to that observed for rats fed with only B. breve 2L, i.e. 1.20E + 03 CFU/gr at T2 and 9.58E + 02 CFU/gr at T3. Thus, supplementation of CG in the standard diet appeared to cause an increase in the abundance of B. breve 2L, thus enhancing gut colonization/ persistence of this strain. Characterization of CG effects on the rat gut microbiota composition. Evaluation of the gut microbiota composition of the animals enrolled in this study was performed by 16 S gene rRNA microbial profiling analyses on fecal samples collected during the trial. Analysis of sequencing data produced a total of 4,581,503 quality-filtered reads with an average of 47,724 reads per sample (Table?S4). Evaluation of alpha-diversity, i.e. the biodiversity, of the collected rat fecal samples was performed through analysis of rarefaction curves constructed with 10 sub-sampling of the whole sequenced datasets. Notably, alpha-diversity data revealed, as expected, that at T0 all animal groups possessed similar microbiota diversity (t-test p-values > 0.05) (Fig.?3). T0 samples were used as ?Control? non-treated group for analysis of T1, T2 and T3 time points. Intriguingly, analysis of data collected for T1 and T2 revealed that CG supplementation appeared to reduce the microbiota diversity, as observed for CG and CG + Breve2L compared to T0 and Breve2L sample groups (Fig.?4). These observations were confirmed by statistical analysis by means of t-test at 20,000 reads of CG and CG + Breve2L with respect to the same animal at T0, resulting in p-values < 0.05 (Fig.?4). Furthermore, data collected at T3 showed that cessation of CG supplementation does not reverse the biodiversity of the CG and CG + Breve2L groups to pre-treatment levels. These observations indicate that CG modulates the microbiota through selection of specific bacterial taxa, with subsequent reduction of the overall gut microbiota diversity, which persists (at least for the period tested) following termination of CG supplementation. Similar data were previously observed for other prebiotic compounds7. Beta-diversity analysis did not reveal any statistically significant differences between sample groups at T0 (PERMANOVA p-value > 0.05) and confirmed the modulatory effect exerted by CG at T1 and T2 (Fig.?5). In fact, PERMANOVA statistical analysis between CG and CG + Breve2L groups at T1 and T2 with respect to T0 datasets resulted in all cases in p-values < 0.05 (Fig.?5). Notably, cessation of CG supplementation causes the CG and CG + Breve2L sample groups to cluster together with Breve2L and Control datasets at T3 (Fig.?5). Thus, beta-diversity data confirms the modulatory effect exerted by CG on the animals? gut microbiota composition. Taxonomic composition of all samples was reconstructed at phylum and genus levels in order to detail the impact of CG on the gut microbiota (Supplementary File?1). A t-test between the relative abundance of each profiled bacterial genus observed in the CG, Breve2L and CG + Breve2L groups was performed and compared between the various time points. Furthermore, statistically significant taxa observed for each comparison were mapped in order to precisely identify those bacterial genera whose relative abundance was increased or decreased upon CG and/or B. breve 2 L supplementation (Fig.?6). Intriguingly, 10 and 12 taxa showed increased (when compared to T0) relative abundance during CG supplementation at T1 and T2, respectively, followed by a decrease in abundance at T3, i.e. after cessation of CG supplementation (Fig.?6). In this context, it is worth mentioning that Prevotella 1, U. m. of Prevotellaceae family and Eubacterium ventriosum group (Lachnospiraceae family) increased their abundance by 100%, 145.4% and 2072.5% in T2 when compared to T0 (Fig.?6) (Supplementary File?1). In contrast, 42 and 38 taxa showed decreased relative abundance at T1 and T2, respectively, followed by an increase in abundance at T3 (Fig.?6) (Supplementary File?1). The observed decreased relative abundance of Bacteroides as well as Ruminococcus (?15.5% and ?67.2% in T2 when compared to T0) and the simultaneous increase in relative abundance of Prevotella suggest that CG promotes a shift of the gut microbiota towards a Type 2 enterotype, i.e. Prevotella-driven enterotype19. Intriguingly, CG was observed to induce increased relative abundance of Akkermansia , a taxon with health-promoting activities20 and reduced relative abundance of Peptoclostridium, a genus encompassing opportunistic pathogens21. Analysis of the impact of B. breve 2L and CG administration on the animal microbiota could not be accurately performed through 16S rRNA gene microbial profiling due to the low relative abundance of bifidobacteria in the assessed fecal samples. For this reason, we performed a detailed cataloguing of bifidobacterial communities through the use of a previously published bifidobacterial ITS profiling approach22,23, allowing a detailed cataloguing of the bifidobacterial population down to the (sub)species level. Cataloguing of the rat fecal bifidobacterial communities residing. In order to detail the impact of CG supplementation on the bifidobacterial population colonizing the rat gut, we performed bifidobacterial ITS profiling of DNA extracted from all fecal samples collected in this study. Sequencing produced a total of 316,080 reads, with an average of 3,293 reads per sample (Table?S5). In order to detail the modulatory effects of CG, we performed inspection of (sub)species-level profiles and statistical analysis through t-test of CG, Breve2L as well as CG + Breve2L groups of samples at T1, T2 and T3 compared to T0 (Supplementary File?2). As expected, an increase in the relative abundance of B. breve was observed in the Breve2L and CG + Breve2L groups at T1 and T2, which was followed by a relative decrease in B. breve 2L abundance at T3 (Fig.?7). Moreover, data of samples constituting the CG group revealed that CG supplementation induces a statistically significant increase in relative abundance (p-value < 0.05) of eight bifidobacterial species at T2 (Supplementary File?2). Furthermore, when CG was supplemented in combination with B. breve 2L, this bifidobacterial species (and most likely the B. breve 2L strain) was shown to reach a higher relative abundance as compared to those animals only receiving B. breve 2L (Fig.?7), as also indicated by qPCR results (Fig.?2). These results indicate that CG exerts a species-specific modulation of the bifidobacterial population harbored by the rat gut. Conclusion Evaluation of the potential prebiotic features of chitin-glucan toward bifidobacteria under in vitro conditions highlighted the ability of 100 bifidobacterial strains to use this substrate as its sole carbon source. The bifidobacterial species that was most effective in utilizing CG for growth was shown to be B. breve. Notably, these results reflect the ability of this typical infant gut colonizer to utilize host-produced glycans, whose molecular structure partially resembles that of CG. RnaSeq transcriptomics data was performed for the B. breve strain showing the best CG utilization capabilities, i.e. B. breve 2L, pointing out a modulation of the B. breve 2L genes involved in the transport and metabolism of hexose sugars. The in vivo trials in a rat model clearly supported the notion that CG exploits a clear bifidogenic effect on bifidobacteria and in particular of small number of bifidobacterial species such as B. breve. These findings reinforce the potential of GC in modulating and shaping the bifidobacterial communities especially in ecological conditions where bifidobacteria are depleted (e.g., very often associated with metabolic disorders or gut diseases). In this context, the capability of bifidobacteria to use the CG as carbon source and the subsequent degradation of this bio-polymer in simpler derivatives, i.e. chitooligosaccharides (COS), may support the notion that CG act as a prebiotic. Indeed, it has been largely demonstrated that COS are able to inhibit the growth of pathogenic bacteria24,25. Moreover, CG consumption was shown to reduce body weight increment in rats, thus pointing at CG as an interesting novel prebiotic for prevention and treatment of obesity. Such in vivo data should be further confirmed by clinical trials performed in human beings consuming CG or CG-based products (e.g., symbiotic products) in those categories of individuals where bifidobacterial abundance is naturally low (e.g., in the elderly) or is depleted as a consequence of metabolic disorder (e.g., constipation) or diseases (auto-immune diseases) as well as antibiotic therapy. It is also arguable that cross-feeding interactions might be established by the different members of the bifidobacterial communities as well as with the different members of the human gut microbiota for the complete metabolism of CG. This was previously shown in co-participated trophic interactions15,26 where a partner partially metabolizes GC in favor to another microorganism that is genetically incapable to utilize this substrate. Thus, cross-feeding might represent another valuable way exploited by CG to induce a more general prebiotic effect in the human gut. Methods Strains and culture conditions. Bifidobacterium strains used in this study are listed in Table?1. Strains were routinely grown anaerobically in De Man, Rogosa, Sharpe (MRS) medium (Scharlau) containing 2% glucose (w/v), which was supplemented with 0.05% L-cysteine-HCl and incubated at 37 ?C for 24 h. Anaerobic conditions were achieved by the use of an anaerobic cabinet (Ruskin), in which the atmosphere consisted of 17% CO2, 80% N2, and 2.99% H2. Ethical statement. All experimental procedures and protocols involving animals were approved by the Veterinarian Animal Care and Use Committee of Parma University (approved protocol 370/2018-PR) and conducted in compliance with the European Community Council Directives dated 22 September 2010 (2010/63/UE). Chitin-glucan growth assay. CG (KitoZyme, Belgium) was added to MRS without glucose (MRS w/o glu) at a final concentration of 0.5% (w/v), as previously assessed through growth experiments involving various bifidobacteria on different carbon sources14. This growth medium was termed MRS w/o glu+ CG and subjected to sterilization employing a modified Tyndallization approach consisting of two thermal cycles at 80?C for 30 minutes each alternated with cooling in ice. For all growth tests, cells were recovered from an overnight culture and turbidity was measured at 600 nm, using a biophotometer (Eppendorf). A growth tube containing 6 mL of MRS w/o glu + CG was inoculated with active viable bacterial cells diluted to an OD600nm of ~1.0, obtaining a final inoculum with an OD600nm of ~0.1. Cultures were grown in biologically independent triplicates and the resulting growth data sets were expressed as the means from these replicates. Moreover, positive (MRS) and negative (MRS w/o glu) growth controls were performed. Cultures were incubated under anaerobic conditions at 37 ?C for 24 h. Cell growth was monitored using a Thoma cell counting chamber (Herka). For some bifidobacterial strains, cell growth was also monitored by viable cell count in MRS Agar. For this purpose, strains were inoculated as mentioned above in MRS w/o glu + CG. Then, 1mL of each grown strain was serially diluted in PBS (Phosphate buffered saline) and plated on MRS agar. Plates were incubated under anaerobic conditions at 37 ?C for 48 h. Bacterial growth was assessed by colony counting. RNA-Seq transcriptomic analysis and identification of genes induced by CG. Bacterial cells were recovered from an overnight culture and turbidity was measured at 600 nm, using a biophotometer (Eppendorf). A growth tube containing 40 mL of complete MRS (reference condition) as well as 40 ml of MRS w/o glu CG was inoculated with viable bacterial cells diluted to an OD600nm of ~1.0, obtaining a final inoculum with an OD600nm of ~0.1. All growth conditions were performed by incubation in anaerobic cabinet at 37 ?C. Following inoculation, 1 0 growth was monitored and when an OD600nm value between 0.6 and 0.8 (exponential phase) cells were centrifuged at 6000 rpm for 5 min. Finally, prior to RNA extraction cells were frozen at ?80 ?C. Growth assays were carried out in triplicate. Total RNA was isolated from B. breve 2L cultures grown in MRS medium (Scharlau, Italy) as well as MRS w/o glu + GC. The obtained cell pellet was resuspended in 1ml of QIAZOL (Qiagen, United Kingdom) and placed in a tube containing 0.8 g of glass beads (diameter, 106 ?m; Sigma). Cells were lysed by shaking the mix on a Precellys 24 homogenizer (Bertin instruments, France). The mixture was then centrifuged at 13,000rpm for 15 min, and the RNA-containing upper phase was recovered. RNA was further purified using RNeasy mini kit (Qiagen, UK) as reported in the manufacturer?s instructions. RNA quality was checked by a Tape station 2200 (Agilent Technologies, USA) analysis. RNA concentration and purity were evaluated by Picodrop microliter spectrophotometer (Picodrop, UK). For RNA sequencing, 2 ?g of total RNA was treated to remove ribosomal RNA by the Ribo-Zero Magnetic Kit (Illumina), followed by purification of the rRNA-depleted sample by ethanol precipitation. RNA was processed according to the manufacturer?s protocol. The efficacy of rRNA depletion was checked by a Tape station 2200 (Agilent Technologies). Then, 500ng of rRNA-depleted RNAs was fragmented using Bioruptor NGS ultrasonicator (Diagenode, USA) followed by size evaluation using Tape station 2200 (Agilent Technologies). A whole transcriptome library was constructed using the TruSeq Stranded RNA LT Kit (Illumina). Samples were loaded into a Flow cell V2 75 cycles (Illumina) as reported by the technical support guide. Sequencing reads were mapped with Burrows-Wheeler Aligner (BWA)27 to the genomic sequence of B. breve 2L (NCBI accession number AWUG00000000.1). Reads Per Kilobase per Million mapped reads (RPKM) of each gene were assessed using Artemis software28. Animal housing. Experiments involved 5-month-old male wild-type Groningen rats (Rattus norvegicus). After weaning, rats were housed in same sex sibling groups in rooms under humidity (50 ? 10%) and temperature-controlled conditions (22 ? 2 ?C), a 12-h light-dark cycle (lights on at 7 a.m.), and with food and water available ad libitum. Experimental design of the in vivo trials. From the initiation of the experiments, rats were housed individually in polymethyl methacrylate (Plexiglas?) cages (39 cm ? 23 cm ? 15 cm). The first week represented an acclimatization period, during which rats continued to consume a standard chow diet supplemented with an oral administration of 500 ? l of sucrose solution (2%) in order to adapt the rats to drink from the syringe. Maintaining their habits, rats could represent the negative control of themselves, acting as the baseline for subsequent microbiota analyses29. For the following two weeks (14 days), rats (n = 24) were randomized to 3 groups: a first group fed with standard diet supplemented with CG [90% standard diet (w/w) + 10% KiOnutrime-CG from KitoZyme, Belgium; CG group], a second group fed with a standard diet and an oral treatment with B. breve 2L (Breve2L group) and a third group fed with standard diet supplemented with CG (same composition of the CG group) and oral treatment with B. breve 2L (CG + Breve2L group) (Table?2). The CG group rats during the intervention period were orally inoculated with 500 ? l of sucrose solution (2%) in order to maintain the same condition for the three groups except for the experimental variables. Finally, during the last week of the wash-out period, all animals returned to the standard chow diet. Standard diet consisted of 54.61% nitrogen-free extract (mainly represented by starch and hemicellulose), 5.54% fibers, 19.42% protein, 11.09% water, 2.58% lipids, and 6.76% ash (non-organic mineral matter) (3.9 kcal/g; 4RF21, Mucedola, Italy). The percentage of supplementation of the substrate CG to standard chow diet was fixed at 10% (w/w) as previously illustrated6. Food intake (FI) and body weight (BW) were measured daily and weekly, respectively, and the BW changes and FI were calculated as previously described7. Evaluation of Bifidobacterium breve cell numbers by qPCR. Quantitative PCR (qPCR) was assessed using the species-specific primers BiBre1 (5?-CCGGATGCTCCATCACAC-3?) and BiBre2 (5?-ACAAAGTGCCTTGCTCCCT-3?). qPCR was performed using GoTaq qPCR Master Mix (Promega, USA) on a CFX96 system (BioRad, CA, USA) following previously described protocols30. PCR products were detected with SYBR green fluorescent dye and amplified according to the following protocol: one cycle of 95?C for 2 min, followed by 40 cycles of 95 ?C for 3 s and 56 ?C for 30 s. The melting curve was 65?C to 95 ?C with increments of 0.5 ?C/s. In each run, negative controls (no DNA) were included. A standard curve was built using the CFX96 software (BioRad). Fecal bacterial DNA extraction and 16S rRNA/ITS microbial profiling. Fecal samples were sub jected to DNA extraction using the QIAmp DNA Stool Mini Kit following the manufacturer?s instructions (Qiagen). Partial 16S rRNA gene sequences were amplified from extracted DNA using primer pair Probio_Uni/ Probio_Rev, which targets the V3 region of the 16S rRNA gene sequence31. Partial ITS sequences were amplified from extracted DNA using the primer pair Probio-bif_Uni/Probio-bif_Rev, which targets the spacer region between the 16S rRNA and the 23S rRNA genes within the ribosomal RNA (rRNA) locus22. Illumina adapter overhang nucleotide sequences were added to the partial 16S rRNA gene-specific amplicons and to the generated ITS amplicons of 200 bp, which were further processed using the 16S Metagenomic Sequencing Library Preparation Protocol (Part No. 15044223 Rev. B?Illumina). Amplifications were carried out using a Verity Thermocycler (Applied Biosystems). The integrity of the PCR amplicons was analyzed by electrophoresis on a 2200 TapeStation Instrument (Agilent Technologies, USA). DNA products obtained following PCR-mediated amplification of the 16S rRNA gene sequences were purified by a magnetic purification step involving the Agencourt AMPure XP DNA purification beads (Beckman Coulter Genomics GmbH, Bernried, Germany) in order to remove primer 1 1 Samples WT1A ? WT8A WT1B ? WT8B WT1C ? WT8C WT1D ? WT8D WT9A ? WT16A WT9B ? WT16B WT9C ? WT16C WT9D ? WT16D WT17A ? WT24A WT17B ? WT24B WT17C ? WT24C WT17D ? WT24D Type of intervention Time points CG group Breve2L group CG + Breve2L group dimers. DNA concentration of the amplified sequence library was determined by a fluorimetric Qubit quantification system (Life Technologies, USA). Amplicons were diluted to a concentration of 4 nM, and 5 ? L quantities of each diluted DNA amplicon sample were mixed to prepare the pooled final library. 16S rRNA gene and ITS sequencing were performed using an Illumina MiSeq sequencer with MiSeq Reagent Kit v3 chemicals. Metagenomics analyses. Following sequencing, the.fastq files were processed using a custom script based on the QIIME software suite32. Paired-end read pairs were assembled to reconstruct the complete Probio_Uni/ Probio_Rev and Probio-bif_Uni/Probio-bif_Rev amplicons. Quality control retained sequences with a length between 140 and 400 bp and mean sequence quality score >20 while sequences with homopolymers >7 bp and mismatched primers were omitted. In order to calculate downstream diversity measures (alpha and beta diversity indices, Unifrac analysis), 16S rRNA Operational Taxonomic Units (OTUs) were defined at ?99% sequence homology using uclust33 and OTUs with less than 10 sequences were filtered. ITS Operational Taxonomic Units (OTUs) were defined at 100% sequence homology using uclust33. All reads were classified to the lowest possible taxonomic rank using QIIME32 and a reference dataset from the SILVA database34 for 16S rRNA data or an updated version of the bifidobacterial ITS database22 for ITS data. Biodiversity of the samples (alpha-diversity) were calculated with Chao1 and Shannon indexes. Similarities between samples (beta-diversity) were calculated by unweighted uniFrac35. The range of similarities is calculated between the values 0 and 1. PCoA representations of beta-diversity were performed using QIIME32. Statistical analysis. For in vitro trials, statistical significance between means was analyzed using the Student?s t-test. Values are expressed as the means ? standard errors from three experiments. Statistical calculations were performed using the software program GraphPad Prism 5 (La Jolla, CA, USA). For in vivo experiments, two-way ANOVA for repeated measures with ?group? as between-subject factor (three levels: CG group, Breve2L group and CG + Bbreve2L group) was performed for: (i) BW changes, with ?time? as within-subject factor (three levels: T1, T2, T3); (ii) FI-to-BW ratio, with ?time? as within-subject factor (four levels: week1, week2, week3 and week4). Follow up analysis was performed using Student?s t-test, with a Bonferroni correction for multiple comparisons. Furthermore, the difference between the rarefaction curves and relative abundance of taxa, as well as the differential abundance of bacterial genera were statistically analysed by t-test. All statistical analyses were performed through SPSS software (www.ibm.com/software/it/analytics/spss/). Data Deposition. Raw sequences of 16S rRNA gene profiling and bifidobacterial ITS profiling as well as RNAseq data are accessible through SRA study accession numbers SRP164382 and SRP164394. 1 2 Acknowledgements This study was supported by KitoZyme, Belgium. Financial support has been provided by a grant from the Walloon Region (ADIPOSTOP Project, convention 7366). This research benefited from the HPC (High Performance Computing) facility of the University of Parma, Italy. We furthermore thank GenProbio srl for financial support of the Laboratory of Probiogenomics. Author Contributions G.A., C.M. and S.D. designed and performed the experiments and wrote the manuscript. G.A. and S.D. performed the experiments. C.M., L.M. and F.B. performed the bioinformatics analyses. L.C. and R.S. performed the in vivo experiments. T.R. and S.M. provided the chitin-glucan charbohydrate. M.C.O. and S.M. contributed to the manuscript preparation. F.T. and A.S. participated in the design of the study. D.v.S. participated and supervised the study. M.V. conceived the study, participated in its design and coordination, and contributed to the manuscript preparation. All authors reviewed the manuscript. All authors read and approved the final manuscript. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-42257-z. Competing Interests: The authors declare no competing interests. 1 3 Publisher?s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article?s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article?s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. 1. Marzorati , M. et al. High-fiber and high-protein diets shape different gut microbial communities, which ecologically behave similarly under stress conditions, as shown in a gastrointestinal simulator . Molecular nutrition & food research 61 , https://doi. org/10.1002/mnfr.201600150 ( 2017 ). 2. Araujo , D. , Freitas , F. , Sevrin , C. , Grandfils , C. & Reis , M. A. Co-production of chitin-glucan complex and xylitol by Komagataella pastoris using glucose and xylose mixtures as carbon source . Carbohydr Polym 166 , 24 - 30 , https://doi.org/10.1016/j.carbpol. 2017 . 02 .088 ( 2017 ). 3. Meichik , N. R. & Vorob'ev , D. V. Chitin-glucan complex in cell walls of the Peltigera aphthosa lichen . Prikl Biokhim Mikrobiol 48 , 340 - 345 ( 2012 ). 4. Commission , E. Commission Decision of 2 February 2011 authorising the placing on the market of a chitin-glucan from Aspergillus niger as a novel food ingredient under Regulation (EC) No 258/97 of the European Parliament and of the Council . Official Journal of the European Union ( 2011 ). 5. Berecochea-Lopez , A. et al. Fungal chitin-glucan from Aspergillus niger efficiently reduces aortic fatty streak accumulation in the high-fat fed hamster, an animal model of nutritionally induced atherosclerosis . Journal of agricultural and food chemistry 57 , 1093 - 1098 , https://doi.org/10.1021/jf803063v ( 2009 ). 6. Neyrinck , A. M. et al. Dietary modulation of clostridial cluster XIVa gut bacteria (Roseburia spp .) by chitin-glucan fiber improves host metabolic alterations induced by high-fat diet in mice . J Nutr Biochem 23 , 51 - 59 , https://doi.org/10.1016/j.jnutbio. 2010 . 10 .008 ( 2012 ). 7. Ferrario , C. et al. How to Feed the Mammalian Gut Microbiota: Bacterial and Metabolic Modulation by Dietary Fibers . Frontiers in microbiology 8 , 1749, https://doi.org/10.3389/fmicb. 2017 . 01749 ( 2017 ). 8. Gibson , G. R. & Roberfroid , M. B . Dietary modulation of the human colonic microbiota: introducing the concept of prebiotics . J Nutr 125 , 1401 - 1412 , https://doi.org/10.1093/jn/125.6. 1401 ( 1995 ). 9. Slavin , J. Fiber and prebiotics: mechanisms and health benefits . Nutrients 5 , 1417 - 1435 , https://doi.org/10.3390/nu5041417 ( 2013 ). 10. Turroni , F. et al. Diversity of bifidobacteria within the infant gut microbiota . PloS one 7 , e36957, https://doi.org/10.1371/journal. pone. 0036957 ( 2012 ). 11. Duranti , S. et al. Elucidating the gut microbiome of ulcerative colitis: bifidobacteria as novel microbial biomarkers . FEMS microbiology ecology 92 , https://doi.org/10.1093/femsec/fiw191 ( 2016 ). 12. Ventura , M. , Turroni , F. , Lugli , G. A. & van Sinderen , D. Bifidobacteria and humans: our special friends, from ecological to genomics perspectives . Journal of the science of food and agriculture 94 , 163 - 168 , https://doi.org/10.1002/jsfa.6356 ( 2014 ). 13. Collins, S. & Reid , G. Distant Site Effects of Ingested Prebiotics. Nutrients 8 , https://doi.org/10.3390/nu8090523 ( 2016 ). 14. Turroni , F. et al. Genome analysis of Bifidobacterium bifidum PRL2010 reveals metabolic pathways for host-derived glycan foraging . Proceedings of the National Academy of Sciences of the United States of America 107 , 19514 - 19519 , https://doi.org/10.1073/ pnas.1011100107 ( 2010 ). 15. Egan , M. , O 'Connell Motherway , M. , Ventura , M. & van Sinderen , D. Metabolism of sialic acid by Bifidobacterium breve UCC2003 . Applied and environmental microbiology 80 , 4414 - 4426 , https://doi.org/10.1128/AEM.01114- 14 ( 2014 ). 16. O 'Connell Motherway , M. et al. Functional genome analysis of Bifidobacterium breve UCC2003 reveals type IVb tight adherence (Tad) pili as an essential and conserved host-colonization facto . r. Proceedings of the National Academy of Sciences of the United States of America 108 , 11217 - 11222 , https://doi.org/10.1073/pnas.1105380108 ( 2011 ). 17. Ventura , M. , Turroni , F. , Motherway , M. O. , MacSharry , J. & van Sinderen , D. Host-microbe interactions that facilitate gut colonization by commensal bifidobacteria . Trends in microbiology 20 , 467 - 476 , https://doi.org/10.1016/j.tim. 2012 . 07 .002 ( 2012 ). 18. Turroni , F. et al. Glycan Utilization and Cross-Feeding Activities by Bifidobacteria . Trends Microbiol , https://doi.org/10.1016/j. tim. 2017 . 10 .001 ( 2017 ). 19. Arumugam , M. et al. Enterotypes of the human gut microbiome . Nature 473 , 174 - 180 , https://doi.org/10.1038/nature09944 ( 2011 ). 20. Collado , M. C. , Derrien , M. , Isolauri , E., de Vos , W. M. & Salminen , S. Intestinal integrity and Akkermansia muciniphila, a mucindegrading member of the intestinal microbiota present in infants, adults, and the elderly . Applied and environmental microbiology 73 , 7767 - 7770 , https://doi.org/10.1128/AEM.01477- 07 ( 2007 ). 21. Milani , C. et al. Gut microbiota composition and Clostridium difficile infection in hospitalized elderly individuals: a metagenomic study . Sci Rep 6 , 25945 , https://doi.org/10.1038/srep25945 ( 2016 ). 22. Milani , C. et al. Evaluation of bifidobacterial community composition in the human gut by means of a targeted amplicon sequencing (ITS) protocol . FEMS microbiology ecology 90 , 493 - 503 , https://doi.org/10.1111/ 1574 - 6941 .12410 ( 2014 ). 23. Milani , C. et al. Unveiling bifidobacterial biogeography across the mammalian branch of the tree of life . ISME J, https://doi. org/10.1038/ismej. 2017 . 138 ( 2017 ). 24. Okawa , Y. , Kobayashi , M. , Suzuki , S. & Suzuki , M. Comparative study of protective effects of chitin, chitosan, and N-acetyl chitohexaose against Pseudomonas aeruginosa and Listeria monocytogenes infections in mice . Biol Pharm Bull 26 , 902 - 904 ( 2003 ). 25. Simunek , J. , Koppova , I. , Filip , L. , Tishchenko , G. & Belzecki , G. The antimicrobial action of low-molar-mass chitosan, chitosan derivatives and chitooligosaccharides on bifidobacteria . Folia Microbiol (Praha) 55 , 379 - 382 , https://doi.org/10.1007/s12223-010- 0063- 0 ( 2010 ). 26. Egan , M. et al. Cross-feeding by Bifidobacterium breve UCC2003 during co-cultivation with Bifidobacterium bifidum PRL2010 in a mucin-based medium . Bmc Microbiol 14 , 282, https://doi.org/10.1186/s12866-014-0282- 7 ( 2014 ). 27. Li , H. & Durbin , R. Fast and accurate short read alignment with Burrows-Wheeler transform . Bioinformatics 25 , 1754 - 1760 , https:// doi.org/10.1093/bioinformatics/btp324 ( 2009 ). 28. Rutherford , K. et al. Artemis: sequence visualization and annotation . Bioinformatics 16 , 944 - 945 ( 2000 ). 29. Turroni , F. et al. Deciphering bifidobacterial-mediated metabolic interactions and their impact on gut microbiota by a multi-omics approach . ISME J 10 , 1656 - 1668 , https://doi.org/10.1038/ismej. 2015 . 236 ( 2016 ). 30. Milani , C. et al. Bifidobacteria exhibit social behavior through carbohydrate resource sharing in the gut . Sci Rep 5 , 15782 , https://doi. org/10.1038/srep15782 ( 2015 ). 31. Milani , C. et al. Assessing the fecal microbiota: an optimized ion torrent 16S rRNA gene-based analysis protocol . PloS one 8 , e68739, https://doi.org/10.1371/journal.pone. 0068739 ( 2013 ). 32. Caporaso , J. G. et al. QIIME allows analysis of high-throughput community sequencing data . Nature Methods 7 , 335 - 336 , https:// doi.org/10.1038/ Nmeth .F. 303 ( 2010 ). 33. Edgar , R. C. Search and clustering orders of magnitude faster than BLAST . Bioinformatics 26 , 2460 - 2461 , https://doi.org/10.1093/ bioinformatics/btq461 ( 2010 ). 34. Quast , C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools . Nucleic Acids Res 41 , D590 -596, https://doi.org/10.1093/nar/gks1219 ( 2013 ). 35. Lozupone , C. & Knight , R. UniFrac: a new phylogenetic method for comparing microbial communities . Appl Environ Microb 71 , 8228 - 8235 , https://doi.org/10.1128/Aem.71.12. 8228 - 8235 . 2005 ( 2005 ). ? The Author(s) 2019


This is a preview of a remote PDF: https://www.nature.com/articles/s41598-019-42257-z.pdf

Giulia Alessandri, Christian Milani, Sabrina Duranti, Leonardo Mancabelli, Thibaut Ranjanoro, Salvatore Modica, Luca Carnevali, Rosario Statello, Francesca Bottacini, Francesca Turroni, Maria Cristina Ossiprandi, Andrea Sgoifo, Douwe van Sinderen, Marco Ventura. Ability of bifidobacteria to metabolize chitin-glucan and its impact on the gut microbiota, Scientific Reports, 2019, DOI: 10.1038/s41598-019-42257-z