Synthesis and characterization of phthalyl pullulan nanoparticles (PPNs)
The chemical reaction scheme of phthalyl pullulan is shown in Fig. 1A. The content of phthalate groups in the phthalyl pullulan was confirmed by measurement of 1H-NMR spectroscopy and was estimated by determining the ratio of protons of phthalic acid to sugar protons as the same method previously reported [15]. Also, the morphologies of the PPNs observed by SEM were spherical shapes as shown in Fig. 1B. The internalization of PPNs was confirmed by a confocal laser scattering microscope (CLSM). It was found that the PPNs were internalized into LP as shown in Fig. 1C.
To examine any internal changes in LP by PPNs, LP was treated with PPNs or pullulan, and its commonly secreted SCFAs such as acetate (C2) (Fig. S1A), propionate (C3) (Fig. S1B), and butyrate (C4) (Fig. S1C) were analyzed in vitro. It was found that the total SCFA amount in the culture medium of LP was not changed when treated PPNs whereas a significant increase was observed when treated pullulan (Fig. 1D).
The antimicrobial activities of the PPNs-internalized LP against pathogens EC and Listeria monocytogenes (LM) were checked by co-culture assay and agar diffusion test. The results showed that the antimicrobial activity of the PPNs-internalized LP was much higher than that of untreated LP or pullulan itself by co-culture assay and agar diffusion test (Fig. 1E, F).
Physiological changes in host
The mouse feeding experiment was performed to evaluate the effect of synbiotics against dysbiosis as mentioned in Fig. 2A. The effect of antibiotics treatment to mimic dysbiosis was determined by the extent of decrease of colonies on LB agar compared to that of a saline-treated group (Fig. S2).
Groups fed diets supplemented with PPNs increased body weight at the end of the feeding experiment compared with other groups (Fig. 2B). In particular, the highest body weights were obtained for the T5 groups which were supplemented with LP/P/PPNs combination, while the body weights of the T1 group were decreased on the endpoint. In the case of feed intake, the values were higher in the PPNs-supplemented group (T4 and T5) than in the others (Fig. 2C).
There were also significant changes in colon length and cecum weight by pro-/synbiotics treatment (Fig. 2D, E). The T4 group had the longest colon length and all groups supplemented with synbiotics (T3, T4, and T5) had a longer colon length compared with the T1 group (Fig. 2F), showing that the pattern of the result was similar to that of the body weight. Likewise, a noticeable increase in values of the cecum weights was only observed in the T1 group, to whom only EC was administered, whereas others fed with probiotics or synbiotics showed lower cecum weight similar to that of the C group.
To assess gut barrier restoration by pro-/synbiotics, the levels of serum endotoxin were measured. The T1 group fed only with EC, showed the highest values whereas the levels were significantly lower in the T4 and T5 groups fed with LP/PPNs and LP/P/PPNs respectively (Fig. 2G). Likewise, when injecting fluorescein isocyanate (FITC)-dextran into the guts of mice to determine intestinal permeability, a similar tendency was observed that the T1 group showed the highest level of FITC-dextran in serum while significantly lower in synbiotics-fed groups including T4 and T5 (Fig. S3).
Effects of synbiotics on the gut microbiota
To determine the effects of PPNs on gut microbiota, both culture-in-/dependent analysis were conducted using intestinal contents and their genomic DNA (gDNA).
Firstly, viable cells of coliform bacteria including EC and lactic acid bacteria (LAB) were enumerated by plating intestinal contents onto MacConkey agar and MRS agar respectively. The mean value of viable coliform bacteria of the T1 group was approximately 6 log10 (CFU/mg of intestinal content) whereas other groups treated pro-/synbiotics were under 4 log10 (CFU/mg of intestinal content) (Fig. 3A). Interestingly, the groups treated PPNs (T4 and T5) showed much lower coliform bacteria than T2 and T3.
On the other hand, the LAB amount of groups tended to be contrary to the result of coliform (Fig. 3B). The T1 group (2-3 log10 CFU/mg) was lowest among groups whereas the values of pro-/synbiotics groups (T2, T3, T4 and T5) were all over 4 log10 (CFU/mg) (Fig. 3B). Also, the values of those groups were significantly increased than that of the C group.
To cross-check the results of culture-dependent analysis and investigate overall changes in the gut microbial community by pro-/synbiotics, gDNA-based analysis such as quantitative PCR (qPCR) and 16S rRNA sequencing was performed. By qPCR, similar results were observed in the levels of enteropathogenic Escherichia coli (intimin) (Fig. S4A) and LABs (Lactobacillus spp., Fig. S4B; Bifidobacterium spp., Fig. S4C).
Meanwhile, microbial community dynamics were explored based on 16S rRNA sequencing. Observed OTUs (operational taxonomic units), a microbial richness index, were higher in the order of T1, C, T2, T3, T4 and T5 (Fig. 3C). Other alpha diversity indices such as Shannon (diversity index; Fig. S5A) and Simpson (evenness index; Fig. S5B) had similar patterns that the lowest in the T1 and highest in the T5 were commonly observed. To examine the effect of PPNs as synbiotic partners on the microbial richness, the groups were reorganized by PPNs treatment (C, T1, T2, T3 vs. T4, T5). Interestingly, observed OTUs were higher when treating PPNs (Fig. S5C).
Principal coordinate analysis (PCoA) based on unweighted (Fig. 3D; R2=0.28, P<0.001) and weighted (Fig. 3E; R2=0.85, P<0.0001) UniFrac distances revealed that the gut microbiota was altered by synbiotics treatment, indicating their significant influence by Adonis test. In the PCoA plot, samples were clustered into three distinct groups (C vs. T1 vs T2, T3, T4 and T5). The samples of T2, T3, T4 and T5 were placed between the C and T1 samples, each of which is also distinguished from each other.
Next, the relative abundance of the microbial taxa in each group was compared and it was found that several phyla and genera appeared to be at quite different levels. At the phylum level, all groups shared the following 13 phyla: Actinobacteria, Bacteroidetes, Cyanobacteria, Deferribacteres, Euryarchaeota, Firmicutes, Fusobacteria, Lentisphaerae, Proteobacteria, Spirochaetes, TM7, Tenericutes, and Verrucomicrobia (Table S2). Three dominant phyla, containing more than 95% of total 16S rRNA gene sequences, were Bacteroidetes, Firmicutes, and Proteobacteria at C and T1 group, while Bacteroidetes, Firmicutes, and Tenericutes at T2, T3, T4, and T5 group (Fig. 3F). In particular, Proteobacteria was more abundant in the T1 group than C group, while it was significantly reduced in T2, T3, T4, and T5 groups (Fig. S6).
At the genus level, the gut microbiota of the six groups shared 102 genera (Table S2). Three dominant genera containing more than 55% of total 16S rRNA gene sequences were as follows: 1) C group: an unclassified genus of family S24-7, Helicobacter, and Odoribacter; 2) T1 group: unclassified genera of family Enterobacteriaceae, Erysipelotrichaceae, and, Lachnospiraceae; 3) T2 group: Oscillospira, unclassified genera of family Lachnospiraceae, Rikenellaceae; 4) T3 and T4 group: Oscillospira, the unclassified genus of family Lachnospiraceae, and unclassified genus of family Ruminococcaceae; 4) T5 group: unclassified genera of family Lachnospiraceae, Rikenellaceae, and Ruminococcaceae (Fig. S7A). In particular, Lactobacillus was more abundant in T4 and T5 than other groups (Fig. S7B). Likewise, Bifidobacterium was significantly more abundant in T4 fed with LP/PPN than other groups, which is similar to the results of qPCR (Fig. S7C). Also, Faecalibacterium and unclassified genus of Veillonellaceae family showed significantly higher abundance in the T5 group than others (Fig. S7D, E).
Taken together, supplementing synbiotics, especially LP/PPNs or LP/P/PPNs, modulated gut microbiota in increasing microbial richness, diversity. Concurrently, the relative abundances of several bacteria were different among groups.
Predicted effects of synbiotics on the gut metagenome
To predict the functions of the gut metagenome of each group, the abundances of Kyoto Encyclopedia of Genes and Genome (KEGG) pathways were predicted by PICRUSt software and assessed its prediction accuracy by the Nearest Sequenced Taxon Index (NSTI) scores. The average NSTI score of C, T1, T2, T3, T4, and T5 were 0.18 (±0.02), 0.07 (±0.03), 0.12 (±0.01), 0.17 (±0.01), 0.15 (±0.01), and 0.15 (±0.01), respectively, which were in similar with other mammal microbiota studies [27, 28]. Subsequently, linear discriminant analysis (LDA) effect size (LEfSe) analysis was performed to determine KEGG pathways whose abundances were different among groups.
The effects of synbiotics on the KEGG were predicted by comparing T1 and T5 groups. It was found that several significantly different KEGG pathways were identified between the two groups (Fig. 4A). For example, ‘Metabolism’, ‘Amino acid metabolism’, ‘Replication and repair’, ‘Cellular processes’ were predicted at significantly higher levels in the T5 group whereas ‘Infectious diseases’, ‘Lipopolysaccharide biosynthesis’, ‘Bacterial secretion system’, ‘Membrane transport’ were predicted at significantly higher levels in T1 group. A similar result was observed when compared with T1 and T4 (Fig. S8A).
The effects of PPNs on the intestinal microbiome were also predicted by comparing T2 and T4 groups (Fig. 4B). T4 groups showed higher levels within ‘Membrane transport’, ‘Transporters’, ‘ABC transporters’, ‘Carbohydrate metabolism’, ‘Transcription factors’, ‘Transcription’ whereas T2 groups showed higher levels within ‘Lipopolysaccharide biosynthesis’, ‘Pore’s ion channels’, ‘Folding sorting and degradation’, ‘Glycan biosynthesis and metabolism’ in T2 group, respectively. The results suggest that the administration of PPNs with probiotics may upregulate gene expression of the gut microbiome to induce the internalization of outer PPNs. Also, the PPNs are suggested to downregulate the biosynthesis of glycan such as lipopolysaccharide (LPS).