Supplementary Fig. 1
Bacterial Sequencing and processing. (A) Around 33.91 Gbases sequences were generated by Illumina NextSeq using 2×150 paired-end reads. (B) Each sample yielded a median of 1.9 Gbases. Sequenced were trimmed and stitched reads mapped to the human genome were removed (42.12%, 0.05%). After quality control, the median number of quality-filtered reads per samples was 10,164,901.
dmj-2019-0202-s003.pdf
Supplementary Fig. 2
Bacterial summary taxonomic composition. High quality reads classified using Greengenes v. 13_8 as the reference database. For each taxonomic rank, we aggregated operational taxonomic units then plotted the relative abundance of the most prevalent genera for healthy and type 2 diabetes mellitus (T2DM) with relation to diet, high fiber diet (HFD) and low fiber diet (LFD).
dmj-2019-0202-s004.pdf
Supplementary Fig. 3
Bacterial summary taxonomic composition. High quality reads classified using Greengenes v. 13_8 as the reference database. For each taxonomic rank, we aggregated operational taxonomic units then plotted the relative abundance of the most prevalent genera for type 2 diabetes mellitus categorized based on glycosylated hemoglobin (HbA1c) levels and diet fiber contents, high fiber diet (HFD) and low fiber diet (LFD).
dmj-2019-0202-s005.pdf
Fig. 1Evaluation of community composition among samples. Panel shows non-metric multidimensional scaling (NMDS) ordination. Ordination plots were categorized by color according to the patient health status (Group 1) and by shape according to their diet (Group 2). T2DM, type 2 diabetes mellitus; LFD, low fiber diet; HFD, high fiber diet.
Fig. 2Difference of functional profiles among samples. Panel shows difference in pathway richness (number of unique pathways) for (A) healthy and (B) type 2 diabetes mellitus (T2DM). Diabetic patients had higher richness of functional genes than that of healthy subjects (P<0.001) irrespective to differences in diet. LFD, low fiber diet; HFD, high fiber diet.
Fig. 3Differential abundance testing of functional groups. (A) Panel shows pathways with significant relative abundance (P<0.001). (B) Log fold difference for significantly enriched pathways for type 2 diabetes mellitus (T2DM) and T2DM with high glycosylated hemoglobin (HbA1c) are PWY-5838 30 fold, PWY-2941 0.9 fold, and CRNFORCAT-PWY −0.9 fold in contrast to healthy and T2DM with low HbA1c.
Fig. 4Covariance network of bacterial genera and metabolic pathways. Panel a shows a force-driven network based on the predicted covariances (analysis of variance with P<0.05) between the genera and metabolic pathways identified as statistically altered in type 2 diabetes mellitus as compared to healthy.