Abstract Simple Summary Human gut microbiome plays an important role for health. This study was thus aimed to analyze the predominant species and metabolic routes involved in short-chain fatty acids (SCFAs) production in the human gut microbiome after treatment with copra meal hydrolysate (CMH). Using integrative metagenomics, key predominant bacterial species and metabolic routes involved in cooperative microbiome networks in relation to SCFAs biosynthesis were identified. This suggests that CMH becomes a potential prebiotic diet for modulating and maintaining the gut microbiome implicated in human health. Abstract Gut microbiome plays an essential role in host health, and there is interest in utilizing diet to modulate the composition and function of microbial communities. Copra meal hydrolysate (CMH) is commonly used as a natural additive to enhance health. However, the gut microbiome is largely unknown at species level and is associated with metabolic routes involving short-chain fatty acids (SCFAs). In this study, we aimed to analyze, using integrative metagenomics, the predominant species and metabolic routes involved in SCFAs production in the human gut microbiome after treatment with CMH. The effect of CMH treatment on the Thai gut microbiome was demonstrated using 16S rRNA genes with whole-metagenome shotgun (WMGS) sequencing technology. Accordingly, these results revealed that CMH has potentially beneficial effects on the gut microbiome. Twelve predominant bacterial species, as well as their potential metabolic routes, were involved in cooperative microbiome networks under sugar utilization (e.g., glucose, mannose, or xylose) and energy supply (e.g., NADH and ATP) in relation to SCFAs biosynthesis. These findings suggest that CMH may be used as a potential prebiotic diet for modulating and maintaining the gut microbiome. To our knowledge, this is the first study to reveal the predominant bacterial species and metabolic routes in the Thai gut microbiome after treatment with potential prebiotics. Keywords: copra meal hydrolysate, human gut microbiome, metagenomics, metabolism, short-chain fatty acids 1. Introduction The human gut microbiome contains highly diverse microbial communities that play essential roles in human health. The gut microbiome is influenced by several factors, such as host genetics, diet, lifestyle, medication, and environment [[38]1,[39]2]. In addition, consumption of specific nutritional substances, such as prebiotics, can modulate gut microbial composition [[40]3,[41]4]. Prebiotics include diverse types of oligosaccharides that are subjected to microbial fermentation in the gastrointestinal (GI) tract, and they affect microbial communities and their functions, offering energy and maintaining gut homeostasis [[42]5,[43]6]. Manno-oligosaccharides (MOS) consist of a linear chain of mannose and have gained great interest as a prebiotic. MOS can be derived from mannan-rich plants, such as copra meal, a by-product of coconut milk and coconut oil processing. In Thailand, 25 million metric tons of copra meal is produced annually [[44]7]. By enzymatic hydrolysis using β-mannanase, a copra meal hydrolysate (CMH) was obtained as a source for further MOS production. CMH is stable under human gastrointestinal tract conditions. CMH is relatively stable in the small intestine and is easily fermented by lactobacilli and bifidobacteria [[45]8]. Very recently, the impact of CMH on gastrointestinal symptoms and gut microbiome was demonstrated to have a positive relationship with health using 16S rRNA gene sequencing [[46]9]. In the context of taxonomy and metabolic function-associated pathways, the human gut microbiome is largely unknown at the species level, as are the associated metabolic routes involved in short-chain fatty acids (SCFAs). Therefore, this study aimed to analyze the predominant species and metabolic routes involved in SCFAs production in the human gut microbiome after treatment with CMH using integrative metagenomics. Different treatments, such as baseline, placebo, or CMH, on the Thai gut microbiome were demonstrated using 16S rRNA genes and whole-metagenome shotgun (WMGS) sequencing technology. DNA extraction from fecal samples was initially performed, and sequencing data were then analyzed using bioinformatics and systems biology tools and databases for taxonomic profiles and metabolic functions of the gut microbiome. Integrative analysis with 16S rRNA genes and WMGS data revealed a number of predominant bacterial species, as well as a list of potential metabolic functions and associated routes involved in SCFAs production after treatment with CMH. This is the first study to reveal that CMH is a potential prebiotic capable of modulating and maintaining the gut microbiome implicated in human health. 2. Materials and Methods 2.1. Participants and Fecal Sample Collection For participants, thirty-seven Thai adults—from Bangkok and its closer city, Thailand—aged between 18–45 years under a BMI of 18.5–24.0 kg/m^2 were registered with double-blinded, placebo-controlled trials in this cohort. Concerning recruitment, they were informed at King Chulalongkorn Memorial Hospital, Bangkok, Thailand, under stringent inclusion and exclusion criteria i.e., dietary intake, age, and status of health. Notably, the participants had no intestinal diseases or diarrheal symptoms in the months prior to sampling and none of the patients had a family history of colorectal cancer. In addition, the participants had to not receive antibiotics within at least three months as well as probiotics, prebiotics, and synbiotics at least one month before sampling. Participants with allergies to coconuts or food intolerance were excluded. The clinical and demographic characteristics of the study participants in the cohort were described according to Sathitkowitchai et al. (2021) [[47]9] (see [48]Supplementary File S1). For fecal sample collection, participants were assigned to four different groups in a total of 74 samples: baseline CMH (bCMH for 17 samples), baseline placebo (bPB for 20 samples), treatment with CMH (tCMH for 17 samples), and treatment with placebo (tPB for 20 samples). Remarkably, the baseline means that they were not treated with either CMH or PB. In addition, treatment meant that they were subjected to CMH or PB for 21 days. Fecal samples (20 g) were immediately collected at the time of defecation and placed into a collection tube in a cooler bag. The fecal samples were stored at –80 °C for further analysis. This study was approved by the Thai Clinical Trials Registry (trial identification number TCTR20190426003) and the Ethics Committee of King Chulalongkorn Memorial Hospital, Bangkok, Thailand (IRB No. 388/61). All methods were performed in accordance with relevant guidelines and regulations. Written informed consent was obtained from all participants. 2.2. DNA Extraction and Metagenome Sequencing All 74 samples were subjected to DNA extraction using a modified method described by Nakphaichit et al. (2014) [[49]10]. Subsequently, 16S rRNA gene sequencing was performed according to the method described by Sathitkowitchai et al. (2021) [[50]9]. Of the 74 samples, 20 were prepared for WMGS sequencing, annotation, and analysis. Initially, fecal samples were centrifuged for 2 min at 13,000 ×g. The pellet was washed twice, centrifuged in phosphate-buffered saline solution (PBS, 1 mL) at 13,000 ×g for 5 min, and then suspended in 900 µL PBS. The supernatant was then discarded. Total DNA purification with quality control assessments and metagenome sequencing was performed according to the method described by Raethong et al. (2021) [[51]11]. Clean reads, which were not mapped to the human genome, were obtained as the metagenomic dataset of the gut microbiome. 2.3. Microbial Taxonomic Analysis and Functional Annotation of 16S rRNA Gene Sequencing Data and WMGS Datasets Initially, 16S rRNA gene sequence datasets (i.e., 17 bCMH samples, 20 bPB samples, 17 tCMH samples, and 20 tPB samples) were analyzed for microbial taxonomy and function using QIIME 2, PICRUSt2 pipeline package (v2.1.0_b) [[52]12], and MetGEMs toolbox (v1.0) [[53]13]. Exploring the MetGEMs toolbox, Core Function was used to investigate the metabolic functions, the KO IDs in each sample were rank-transformed, and the geometric means of KO IDs of each group were then computed. To analyze WMGS datasets, 20 samples out of 74 samples were randomly selected (5 bCMH samples, 5 bPB samples, 5 tCMH samples, and 5 tPB samples). Alpha- and beta-diversity analyses were initially performed using a vegan package [[54]14] in the R program (version 2.5-6). For alpha-diversity, the observed species and Shannon diversity indices were used to calculate species richness and abundance. Statistical differences in diversity indices between possible pairwise comparison groups were identified using the Wilcoxon rank-sum test. Beta-diversity was computed using Bray–Curtis distances with the metaMDS function in the vegan R package [[55]14]. Beta-diversity analysis was performed to investigate differences between microbial communities across possible pairwise comparison groups. Differences in beta-diversity were visualized through non-metric multidimensional scaling (NMDS) ordination using the ggplot2 R package [[56]15]. Beta-diversity was tested for inference by permutational multivariate analysis of variance (PERMANOVA), as implemented in the ADONIS function from the vegan R package [[57]14], using permutations equal to 999. To test the difference in microbial composition between two or more groups, analysis of similarities (ANOSIM) was also evaluated based on the Bray–Curtis dissimilarity using the vegan R package [[58]14]. For microbial taxonomic analysis, MetaPhlAn (v3.0) [[59]16] was used to detect and quantify individual species with a library of clade-specific markers (mpa_v30_CHOCOPhlAn_201901) as database and to further generate overall metagenome profiles, for example, phyla, families, genera, species, and strain levels. The relative abundance of the microbial taxonomic level was plotted using the ggplot2 package [[60]15] in the R program (v.3.5.3). The Wilcoxon rank-sum test was used to identify statistically significant differences (p < 0.1) between possible pairwise comparison groups. The relative abundance of gene families and pathways was determined using HUMAnN 3. The gene profile was summed into KO IDs and normalized using the built-in script in HUMAnN 3 [[61]17]. To further determine the metabolic function of the gut microbiome, the obtained KO IDs from the 16S rRNA gene sequence and WMGS dataset were searched against the KEGG database using KEGG mapper in KEGG mapping tools [[62]18,[63]19,[64]20]. Once KEGG Orthology (KO) IDs were identified, they were classified into six functional categories: metabolism, genetic information processing, cellular processes, environmental information processing, organism systems, and human diseases. Further, the functional and pathway enrichment analysis was explored under a distinct up-directional p-value < 0.05 for the potential group, for example, pairwise comparison between tCMH and tPB. 2.4. Identification of Predominant Bacterial Species, Potential Metabolic Functions and Associated Routes Involved in Treatment with CMH Using Integrative Analysis To identify the predominant bacterial species, potential metabolic functions, and associated routes when treated with CMH, taxonomic profiles, and metabolic functions obtained from 16S rRNA gene sequences and WMGS datasets were integrated. After analyzing taxonomy and metabolic functional as well as pathway enrichments, the predominant bacterial species and enriched KO IDs under positive median mode and p-value < 0.05 were then considered. Based on the KEGG database, the targets of enriched KO IDs across predominant species were mapped to metabolic pathways, for example, SCFA biosynthesis, using mapper in KEGG mapping tools. Furthermore, literature mining and manual curation were performed to uncover a number of predominant bacterial species, as well as a list of potential metabolic functions and associated routes. 3. Results and Discussion 3.1. Assessment of 16S rRNA Gene Datasets on Microbial Composition and Metabolic Function After treatment of the Thai gut microbiome with CMH (tCMH), the microbial composition in the human gut was assessed using 16S rRNA gene datasets. Four phyla (Firmicutes, Bacteroidetes, Actinobacteria, and Proteobacteria) were most commonly identified. Firmicutes showed the highest relative abundance in the bacterial community, accounting for 80.6%. At the family level, tCMH was found in the top five bacterial families (Lachnospiraceae, Bacteroidaceae, Erysipelotrichaceae, Prevotellaceae, and Rikenellaceae) under the positive median mode ([65]Figure 1A). Ruminococcaceae was reduced in the tCMH group. These results are consistent with those reported by Sathitkowitchai et al., 2021 [[66]9]. Across the other pairwise groups, for example, at the phyla level, the changes in the gut microbiome were similar (see [67]Supplementary File S1). Figure 1. [68]Figure 1 [69]Open in a new tab Microbial composition and metabolic function after treatment with CMH, using 16S rRNA gene sequences datasets. (A) Microbial composition changes at family level. (B) Top 15 associated pathways in carbohydrate metabolism. Note: tPB and tCMH represent treatment with placebo and CMH, respectively. Metabolic function prediction of the gut microbiome focused on metabolism ([70]Figure 1B). The top three associated pathways involved in carbohydrate metabolism were the major enrichments. Considering the positive median mode under statistical significance (p ≤ 0.05), glycolysis/gluconeogenesis, propanoate metabolism, C5-Branched dibasic acid metabolism were targeted in relation to the central carbon metabolism and energy supply. There were target of enzymes e.g., phosphoenolpyruvate carboxykinase (ATP) (EC:4.1.1.49), 3-isopropylmalate dehydrogenase (EC:1.1.1.85), pyruvate, orthophosphate dikinase (EC:2.7.9.1), aldose 1-epimerase (EC:5.1.3.3), fructose-1,6-bisphosphatase II (EC:3.1.3.11), fructose-1,6-bisphosphatase I (EC:3.1.3.11), fructose-1,6-bisphosphatase III (EC:3.1.3.11), fructose-1,6-bisphosphatase I/sedoheptulose-1,7-bisphosphatase (EC:3.1.3.11/3.1.3.37), fructose-1,6-bisphosphatase II/sedoheptulose-1,7-bisphosphatase (EC:3.1.3.11 3.1.3.37), fructose 1,6-bisphosphate aldolase/phosphatase (EC:4.1.2.13/3.1.3.11), pyruvate kinase (EC:2.7.1.40) pyruvate kinase isozymes R/L (EC:2.7.1.40), methylglyoxal synthase (EC:4.2.3.3) (see [71]Supplementary File S1). 3.2. Assessment of Gemicrobial Diversity and Composition from WMGS Datasets Raw WMGS datasets were assessed for 1141.74 Megabases (Mb) of all samples ([72]Table 1). After discarding low-quality reads, adapter and human genome contaminants, a total of 1132.41 Mb of clean reads were retrieved with an average effective rate of 99.18% (see [73]Supplementary File S2) and used for further analysis. Table 1. Summary of WMGS datasets for each group. Groups Raw Reads (Mb) Clean Reads (Mb) Effective Rates * (%) Baseline bPB 273.40 ^# 271.29 99.23 bCMH 276.56 ^# 273.71 98.97 Treatment tPB 301.76 299.58 99.28 tCMH 290.02 287.83 99.24 Total 1141.74 1132.41 99.18 [74]Open in a new tab Note: * Means of effective rates are presented. ^# Data was taken from Raethong et al. 2021, [[75]11]. In addition, bPB, bCMH, tPB, tCMH represent baseline placebo, baseline CMH, treatment with placebo, treatment with CMH, respectively. Alpha- and beta-diversity analyses were performed to determine the microbial diversity across all possible groups. The observed species and Shannon diversity indices showed that each group was comparable, as shown in [76]Figure 2A,B. Based on the Wilcoxon rank-sum test, there were no statistically significant differences in alpha-diversity among all groups. Regarding beta diversity, an NMDS plot of the Bray–Curtis dissimilarity index is shown in [77]Figure 2C. The results also showed that beta-diversity was not significantly different across all groups (ANOSIM analysis: R = –0.236, p > 0.05; ADONIS analysis: R^2 = 0.127; p-value > 0.05). This indicates that the diversity between the baseline and treatment groups was comparable in the context of taxonomic richness and abundance. Figure 2. [78]Figure 2 [79]Open in a new tab Comparison of microbial diversity across all possible groups. (A,B) Boxplots of the alpha-diversity in contexts of observed species and Shannon index of each pairwise group, respectively; (C) Beta-diversity plot of non-metric multidimensional scaling (NMDS) ordination in context of Bray–Curtis dissimilarity among each group. Taxonomic diversity was assessed at species level. In addition, bPB, bCMH, tPB, and tCMH represent baseline placebo, baseline CMH, treatment with placebo, and treatment with CMH, respectively. To initially assess the taxonomic profiles of the gut microbiome, the WMGS datasets were analyzed for all taxa. As presented in [80]Figure 3A,B, we found that four phyla (Firmicutes, Actinobacteria, Proteobacteria, and Bacteroidetes) were comparable across the baselines (bPB and bCMH) and treatments (tPB and tCMH). Among these four phyla, Firmicutes showed the highest relative abundance in the bacterial community for all groups, accounting for an average of 60.99%. The family Lachnospiraceae, belonging to the genus Roseburia in the phylum Firmicutes, was dominant at a high abundance (23.73%). Considering all groups under the top 15 bacterial families (see [81]Supplementary File S2), five dominant families ([82]Figure 3C,D) were found: Bifidobacteriaceae, Enterobacteriaceae, Eubacteriaceae, Lachnospiraceae, and Ruminococcaceae. These results are consistent with those of La-ongkham et al. (2020) [[83]21] in the context of the core taxonomy of the Thai gut microbiome. Figure 3. [84]Figure 3 [85]Open in a new tab Taxonomic profiles of gut microbiome using WMGS datasets. (A,B) Boxplot of relative abundance of microbial phyla across baselines and treatments, respectively. (C,D) Stacked column graph showing the relative abundance of taxonomic families across baselines and treatments, respectively. In addition, bPB, bCMH, tPB, and tCMH represent baseline placebo, baseline CMH, treatment with placebo, and treatment with CMH, respectively. Focusing on tCMH, five alternative dominant families—Ruminococcaceae, Coriobacteriaceae, Erysipelotrichaceae, Veillonellaceae, and Acidaminococcaceae—were observed. Of these families, eight potential species were identified ([86]Supplementary File S2): Acidaminococcus intestini, Agathobaculum butyriciproducens, Collinsella aerofaciens, Collinsella stercoris, Faecalibacterium prausnitzii, Holdemanella biformis, Ruminococcus bicirculans, and Ruminococcus bromii ([87]Figure 4). Notably, two species (p-value < 0.1), A. intestini and A. butyriciproducens, were significantly involved in metabolic functions. Figure 4. [88]Figure 4 [89]Open in a new tab Comparison of relative abundance between tPB and tCMH on eight potential species, using WMGS datasets. For those marked with an asterisk (*), significant difference between tCMH and tPB groups is shown by different letters (p-value < 0.1). 3.3. Assignment of Metabolic Function Underlying KO IDs from WMGS Datasets Concerning on all bacterial families achieved from taxonomic profiles (see [90]Supplementary File S2), metabolic functional analysis was performed using KEGG. A total of 5421 KO identifiers (KO IDs) were identified. Across six functional categories, metabolism was shown to have the highest number of KO IDs (1,712 KO IDs), followed by genetic information processing (365 KO IDs), environmental information processing (550 KO IDs), cellular processes (177 KO IDs), organism systems (17 KO IDs), and human diseases (56 KO IDs), as shown in [91]Figure 5A and [92]Supplementary File S2. As observed, 530 out of 1,712 KO IDs were mainly found in carbohydrate metabolism. Genes involved in glycolysis/gluconeogenesis (72 KO IDs), amino sugar and nucleotide sugar metabolism (49 KO IDs), pentose and glucuronate interconversions (47 KO IDs), fructose and mannose metabolism (46 KO IDs), starch and sucrose metabolism (43 KO IDs), glyoxylate and dicarboxylate metabolism (42 KO IDs), pyruvate metabolism (39 KO IDs), galactose metabolism (38 KO IDs), butanoate metabolism (35 KO IDs), propanoate metabolism (31 KO IDs), TCA cycle (30 KO IDs), pentose-phosphate pathway (30 KO IDs), ascorbate and aldarate metabolism (13 KO IDs), inositol phosphate metabolism (11 KO IDs), and C5-branched dibasic acid metabolism (4 KO IDs) ([93]Figure 5B and [94]Supplementary File S2). Figure 5. [95]Figure 5 [96]Open in a new tab Functional categories underlying KO IDs from WMGS datasets of gut microbiome (A) Main functional categories using KEGG (B) Sub-metabolic categories in carbohydrate metabolism. Metabolic function and pathway enrichment analyses were performed for carbohydrate metabolism. Interestingly, we found that the TCA cycle and pentose and glucuronate interconversions were significantly enriched (p < 0.05) ([97]Figure 5B and [98]Figure 6A; see [99]Supplementary File S2). Figure 6. [100]Figure 6 [101]Open in a new tab Heatmaps represent the enriched metabolic pathways underlying KO IDs (A) KO IDs involved in TCA cycle and pentose and glucuronate interconversions (B) Bacterial species underlying significant KO IDs. Note: Color maps display fourth root of RPKM value. Upon tCMH, four KO IDs were identified: altronate hydrolase ([102]K01685, EC:4.2.1.7), 2-oxoglutarate ferredoxin oxidoreductase subunit delta ([103]K00176, EC:1.2.7.3), l-iditol 2-dehydrogenase ([104]K00008, EC:1.1.1.14), and d-lyxose ketol-isomerase (K09988, EC:5.3.1.15) ([105]Figure 6A) (see [106]Supplementary File S2). These enzymes are essential for 2-oxoglutarate and succinate formation, leading to the biosynthesis of SCFAs. Succinate is a key precursor and plays an important role in either propionate or acetyl-CoA formation, which can then be converted to acetate or butyrate formation [[107]22,[108]23,[109]24]. After SCFA-producing species mapping with the target of four KO IDs, 11 species were identified: Acidaminococcus intestini, Anaerostipes hadrus, Bacteroides dorei, Bacteroides massiliensis, Bacteroides vulgatus, Clostridium saccharolyticum, Desulfovibrio piger, Escherichia coli, Eubacterium siraeum, Roseburia hominis, and Roseburia intestinalis ([110]Figure 6B and [111]Supplementary File S2). These results are consistent with those of Nogal et al. (2021), who showed that SCFAs are produced by enteric and acetogenic commensal microbes, such as Bacteroides spp. and Clostridium spp. [[112]22]. Exhibiting functions involved in succinate formation towards SCFA biosynthesis, for example, propionate, butyrate, and acetate, four out of 11 species (i.e., A. hadrus, R. hominis, A. intestine, and E. siraeum) were interested ([113]Figure 6B) [[114]25,[115]26]. Altogether, our results agree well with earlier reports on propionate/butyrate-producing microbes, such as A. intestini, E. coli, Bacteroides spp., Anaerostipes, Eubacterium and Roseburia [[116]27]. 3.4. Identifying Predominant Bacterial Species and Potential Metabolic Routes Using Integrative Metagenomics To identify the predominant bacterial species, the results obtained from the 16S rRNA gene sequence datasets ([117]Figure 1B) and the WMGS datasets ([118]Figure 6B) were integrated. Post-tCMH, we found 12 predominant species: A. intestini and A. butyriciproducens, A. hadrus, B. dorei, B. massiliensis, B. vulgatus, C. saccharolyticum, D. piger, E. coli, E. siraeum, R. hominis, and R. intestinalis. According to literature supports for all identified species [[119]28], very interestingly, they were SCFA-producing species. The results of mapping the enriched KO IDs identified from these SCFA-producing species onto SCFA biosynthesis are shown in [120]Figure 7. Indeed, it revealed five metabolic pathways—glycolysis/gluconeogenesis, TCA cycle, pentose and glucuronate interconversions, C5-branched dibasic acid metabolism, and propanoate metabolism—in relation to SCFA biosynthesis. Figure 7. [121]Figure 7 [122]Open in a new tab Potential species and related metabolic routes involved in SCFAs biosynthetic pathways. Note: KEGG was used as a pathway database. Among SCFAs, the metabolic routes involved in acetate, propionate, and butyrate production were identified ([123]Figure 7), as they are commonly found in the human gut [[124]23]. In glycolysis/gluconeogenesis, we found 11 KO IDs across five key enzymes involved in nutrient utilization and ATP supply: pyruvate kinase (EC:2.7.1.40), phosphoenolpyruvate carboxykinase (ATP) (EC:4.1.1.49), pyruvate, and phosphate dikinase (EC:2.7.9.1). These enzymes are important for generating key precursors such as oxaloacetate and pyruvate for acetyl-COA towards acetate and butyrate formation [[125]4,[126]28]. Among SCFAs, acetate is the most abundant SCFA, and it is essential for the growth of other microbes in the human gut. Butyrate is also the main energy source for human colonocytes and enterocytes, and can activate intestinal gluconeogenesis for beneficial effects on glucose and energy homeostasis. In the TCA cycle, we also found NADH supply via only one KO ID and one key enzyme, 2-oxoglutarate synthase (EC:1.2.7.3), across A. intestini, B. dorei, B. massiliensis, B. vulgatus, D. piger, and E. siraeum. This enzyme is interesting because it produces succinyl-COA, which is a precursor for further propionate formation [[127]22,[128]29]. Consistent with other populations from earlier reports, propionate is commonly produced by Bacteroidetes, including Bacteroides, for example, B. vulgatus through the succinate pathway [[129]30] when treatment with CMH [[130]31]. For pentose and glucuronate interconversions, we observed three KO IDs and three key enzymes involved in NADH supply: l-iditol 2-dehydrogenase (EC: EC:1.1.1.14). This enzyme is widely distributed and has been described in archaea and bacteria across five bacterial species (i.e., A. hadrus, C. saccharolyticum, E. coli, R. hominis and R. intestinalis). This enzyme is widely distributed in the bacteria. It acts on a number of sugar alcohols, such as D-xylitol, which increases the concentration of propionate, as seen in E. coli [[131]32]. For the remaining C5-branched dibasic acid metabolism and propanoate metabolism, we observed two KO IDs and two key enzymes involved in NADH and precursor supply (i.e., 1,2 propanediol) for propionate, 3-isopropylmalate dehydrogenase (EC:1.1.1.85), and methylglyoxal synthase (EC:4.2.3.3), respectively. As previously reported, propanediol-associated pathway carriers almost exclusively belong to the family Lachnospiraceae, mainly from the genera Ruminococcus and Blautia [[132]33,[133]34]. In addition, propionate can be converted to glucose via intestinal gluconeogenesis. Propionate has antibacterial and anti-inflammatory functions that protect the human intestine against pathogens [[134]28]. Beyond integrative metagenomics throughout this study, the predominant species and potential metabolic routes were revealed to be cooperative microbiome networks under different sugar utilizations (e.g., glucose, mannose, or xylose) in relation to SCFA biosynthesis. Taken with tCMH, these results suggest that CMH is a potential prebiotic for human gut microbiome modulation and maintenance. 4. Conclusions Our integrative metagenomics approach under different treatment regimens was enabled using 16S rRNA gene and WMGS sequencing technology. The data were analyzed using bioinformatics and systems biology. Upon treatment with CMH, we positively identified twelve predominant bacterial species, potential metabolic functions and associated routes. Cooperative microbiome networks under sugar utilization (e.g., glucose, mannose, or xylose) and energy supply (i.e., ATP and NADH) engage in central carbon metabolism for biosynthesizing SCFAs, such as propionate, butyrate, and acetate. These findings suggest that CMH may be a potential prebiotic for modulating the gut microbiome. Therefore, this study sheds light on the gut microbiome–metabolism axis implicated in human health. Acknowledgments