Changes in the Human Gut Microbiome during Dietary Supplementation with Modified Rice Bran Arabinoxylan Compound

1. Introduction

Communities of microorganisms reside within humans, with the gut being the most densely colonised site [1]. Microbial communities that cohabit the intestinal tract, known as gut microbiota, can be beneficial, neutral, or detrimental to the host. Gut microbiota can impact human health by mediating physiological homeostasis through immune function, digestion, vitamin synthesis, and pathogen colonisation [2]. The effects are exerted through interfaces between gut microbiota and intestinal epithelial cells, the immune system, and dietary intake [2,3,4]. Gut microbiota diversity and abundance is associated with health through processes such as reduced inflammation and disorders linked to inflammatory events [3].

Microbial diversity is required for healthy gut function, with a loss of diversity associated with an increased risk of disease [5]. In healthy adults, the gut microbiota comprises eight dominant phyla including Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria, Fusobacteria, and Verrucomicrobia. The Firmicutes and Bacteroidetes together represent 90% of gut microbiota [6,7]. Bacteria in these phyla perform diverse roles in regulating host health [7]. Bifidobacterium species comprise only 2% of total gut microbiota but play a significant role in the breakdown of complex carbohydrates [8,9], protecting the host against pathogens through competitive exclusion, modulating the immune system, and providing vitamins and other nutrients for the host [10]. The genus Lactobacillus includes species that produce short-chain fatty acids (SCFAs) and regulate intestinal transit [11,12]. Production of SCFAs by these bacteria can improve insulin sensitivity and protect against diet-induced obesity [9]. Other bacteria of importance to human health include Roseburia, Faecalibacterium, Ruminococcus, and Bacteroides [13,14,15].

The presence of microbes and their interaction with the host and other gut microbes in the gastrointestinal tract can lead to adverse health outcomes [1]. Metabolic products from gut microbiota have been linked to an increased risk of major adverse cardiovascular events [16] and type 2 diabetes [17]. Supporting an optimal composition of gut microbiota is vital for human health [18]. The gut microbiota in healthy individuals includes shared common microbiota, with variations across age, ethnicity, and socioeconomic status [19,20,21,22,23]. Although relatively stable throughout adulthood, variations in the common microbial genera are observed, with major shifts identified in older adults [23,24,25,26]. Composition of the gut microbiome is also affected by lifestyle [27,28]. Cigarette smoking and alcohol consumption impact the gut microbiome with changes identified with the cessation of smoking, and heavy alcohol consumption [27,28].

Diet has also been shown to be a significant factor affecting gut microbial composition [18,29]. Functional foods aim to improve health and wellbeing and have been shown to interact with gut microbiota [30,31,32]. Prebiotics are non-digestible dietary ingredients that stimulate the growth of microbes in the gut after fermentation [33]. Recent research identifies prebiotics that enhance microbial diversity to promote health and defend against dysbiosis of gut microbiota [4,5]. Arabinoxylans (AXs) are the non-digestible fibre of cereal grains, including wheat, rice, rye, maize, and sorghum [9]. The structures of AXs have been shown to influence substrate fermentation and degradation by gut microbiota [34]. Francios et al. [35] and Kjolbaek et al. [36] showed that supplementation with wheat AX/AX oligosaccharides in humans led to changes in the microbiota profile including increased beneficial bacteria such as Bifidobacterium. Lachnospiraceae have also been shown to be increased in AX-fed mice [37]. The extent and mechanism of the effects of AXs on gut microbiota differs due to variations in the chemical structure and molecular weight [38]. Recent in vivo research has focused on linear AXs due to their bioavailability and simple structures [9,39].

Rice bran arabinoxylan compound (RBAC), known as BioBran, MGN-3, and Ribraxx, is a modified AX compound [40]. RBAC consists of a water-soluble hemicellulose-$\beta$ fraction (degree of polymerisation approximately 200), partially decomposed by enzymes extracted from a cultured medium of Lentinula edodes (shiitake mushrooms) [40]. However, knowledge of its effect on gut microbiota, particularly in healthy adults, is lacking. Therefore, the purpose of this clinical study was to determine the changes in gut microbiota in response to dietary RBAC supplementation in healthy adults. Additional observational analysis was also performed to assess the abundance of the beneficial bacteria Bifidobacterium, Lactobacillus, Roseburia, Ruminococcaceae, and Faecalibacterium.

2. Results

2.1. Participant characteristics and dietary intake

Characteristics of the study participants and their dietary intake from the screening survey and Australian Eating Survey® (AES) are summarised in Table 1. The participant group consisted of balanced numbers of males and females (sex ratio 1:1). Ages ranged from 21 to 56 years, with most participants (n=8) being 30 years or younger (average age of 26.5 years), and two participants above 30 years (average age of 46.5 years). Nine participants did not smoke cigarettes, and six participants consumed alcohol. The dietary analysis showed that most participants did not follow a particular dietary type and were considered omnivores, however, two participants were either vegan or pescatarian. Australian Recommended Food Score (ARFS) was calculated as an average from AES results obtained before intervention commencement and at the completion of the study. For the sampling, 80 faecal samples were collected, with 3 faecal samples lost in transit. A total of 77 samples were analysed for this study.

2.2. DNA quality

There were 14,883,521 high-quality 16S rRNA sequences obtained with a mean of 173,579 forward and reverse reads per sample from 77 samples in addition to 8 repeats for quality control purposes.

2.3. Phylogenetic taxonomy of the gut microbiota

The gut microbiota taxonomy for all samples from each participant were visualised at the phylum level using taxa bar plots (Figure 1). This graph shows that Bacteroidetes (green) and Firmicutes (purple) were the two dominant phyla across all samples. Furthermore, the plot shows that the taxa of gut microbiota in participants vary from each other. However, the ratio of Firmicutes/ Bacteroidetes fluctuated over the time points and was not statistically significant.

2.4. Associations of alpha diversity of gut microbiota with explanatory factors

The rarefaction plot generated using Shannon’s index [41] showed a maximum depth of 47,647 based on the median frequency value from the frequency per sample results. This is sufficient for analysis as a levelling out on the Y-axis was observed at this depth, indicating additional sequences beyond this depth would unlikely result in additional observed features.

The association of alpha diversity metrics with explanatory factors using Kruskal-Wallis analysis was performed on sequences pooled according to explanatory factors (Table 2). These factors were participant, sex, age group, cigarette smoking, alcohol consumption, and ARFS group. Association of alpha diversity was also determined for time point (1-8), experimental phase (baseline, low dose, washout, high dose, and post-intervention), RBAC dosage (0g/day, 1g/day, and 3g/day), and interventional period (baseline to experimental period including washout and post-intervention). Participant was the only factor that showed a statistically significant effect on alpha diversity in Shannon’s evenness (p=9.66e^-8) and Faith’s PD (p=0.002). Statistically significant changes in Shannon’s evenness were shown for alcohol consumption (p=7.92e^-6), cigarette smoking (p=0.008), sex (p=0.003), and ARFS group (p=0.032). There were no significant changes in alpha diversity across age group, time point, experimental phase, RBAC dosage, and interventional period.

2.5. Associations of beta diversity of gut microbiota with explanatory factors

Associations of beta diversity with explanatory factors were determined using PERMANOVA (Table 3). Participant sequence data were grouped according to the explanatory factor being tested using the same technique as alpha diversity analysis. Statistically significant differences in beta diversities were found across participant (p=0.001), sex (p=0.001), alcohol consumption (p=0.001), ARFS group (p=0.001), and age group (p≤0.013). Cigarette smoking was significantly associated with all matrices (p=0.001) except weighted unifrac (p=0.1). The analysis for time point, experimental phase, RBAC dosage, and interventional period did not show statistically significant associations with any beta diversity measurements.

An EMPeror plot of Bray Curtis beta diversity showing clustering according to individual participants is demonstrated in Figure 2. These association analyses showed that participant is the major explanatory factor associated with both alpha and beta diversity patterns. This finding suggests that individuality is a critical determinant of microbial diversity.

These findings suggest the requirement to analyse beta diversity for each participant to determine whether there were statistically significant longitudinal changes from baseline across the experimental period. Individual beta diversity for each participant across samples (1-8) was calculated using weighted unifrac and Bray Curtis methods (Table 4).

The explanatory factors used for this analysis were RBAC dosage, experimental phase, and interventional period. In P02, the interventional period was significantly associated with weighted unifrac (p=0.01) and Bray Curtis (p=0.012) beta diversity methods. In P03, the experimental phase and interventional period were significantly associated with both methods (p≤0.052). P04 observed associated changes in RBAC dosage for both methods (p=0.035 and p=0.018, respectively), as well as the experimental phase with Bray Curtis (p=0.044). In P07, the experimental phase (p=0.025) and interventional period (p=0.025) were significantly associated with Bray Curtis. In P08, RBAC dosage was significantly associated with both methods (p<0.05), and experimental phase and interventional period for Bray Curtis only. P09 observed that dosage and experimental phase were significantly associated (p<0.05) with both methods. It is further noted that P05 and P10 showed marginally significant values in Bray Curtis method for the experimental phase.

The longitudinal findings are illustrated using the weighted unifrac beta diversity distance matrices across the 8 time points in EMPeror plots for selected participants (Figure 3). For P03, all timepoints after baseline (red points) show a significant change. For P03, all timepoints after baseline (red points) show a significant change. The graph for P04 illustrates an increase from baseline in the low dose interventional period (orange points), but no change in the high dose (blue points) or washout (green points) periods. P01 shows no change across any timepoints from the baseline. Technical replicates are presented on the graph and are shown on their associated time point alongside the sample. Please refer to Figure S1 of the supplementary materials for the plots for all participants.

2.6. Detection of the association of microbial diversity with interactions of explanatory factors

Further investigation was carried out to determine whether multiple variables could explain the variation in participants’ gut microbiota and reduce false positives in the beta diversity analysis. The calculated unifrac distance matrix was used to perform an Adonis PERMANOVA test on participant and other explanatory factors, including RBAC dosage, experimental phase, time point, and interventional period. The results showed that the variable participant (R²=0.796, p=0.001) explains approximately 80% of variance in the beta diversity. of the variance in the beta diversity. Regression analysis with linear models showed a statistical significance for the experimental phase (R²=0.019, p=0.026) and intervention (R²=0.012, p=0.011). Hence, after participants’ variation was accounted for, the experimental and interventional periods explain approximately 2% of the variance observed in beta diversity with the remaining 18% of variance unaccounted for.

2.7. Significant taxonomic changes from baseline to interventional period

Taxonomic changes from baseline to the interventional period were statistically assessed using ANCOM. Data for all statistically significant changes with a W value of 5 or above are listed in Table 5. Associated changes were compared from baseline to the total interventional period unless otherwise specified. These results show statistically significant taxonomic changes (from baseline) detected down to a species level for 7 participants - P01, P02, P03, P07, P08, P09 and P10.

In P01, a marked increase in the Hungatella genus with a W value of 40 was revealed and less marked changes across other genera and families were observed. In P02, decreases in Anaerococcus, Corynebacterium, and Finegoldia genera were detected with W values ranging between 48 to 64. Large declines in Dialister (W=94) and Gastranaerophilales (W=86) were observed in P03. A large decrease in Erysipelatoclostridiaceae (W=102) was displayed in P07. In P08, there was a decrease in the Eubacterium siraeum group (W=35) and lower significant changes in other bacterial genera. P09 revealed large increases in multiple bacterial genera, including Eubacterium siraeum, (W=38). This participant also showed smaller increases in Eubacterium hallii group and Stomatobaculum (W=5 for each) and further smaller increases in Anaerococcus, Ruminococcus, Megasphaera, Solobacterium (W=7-8) in either low or high RBAC doses. P10 shows a large decrease in Prevotella (W=20) and smaller changes across other groups. Some changes were only detected during the washout stage, as observed in P01.

3. Discussion

The purpose of this clinical study was to determine changes in the gut microbial composition in response to dietary RBAC supplementation in healthy adults. Our results revealed that alpha and beta diversity were not associated with experimental phase, interventional period, RBAC dosage, or time point. This suggests that the dietary supplementation did not significantly moderate the gut microbiota composition. However, the absence of supplementation effect on the gut microbiota composition should be further explicated. Firstly, the unknown species are not accounted for in alpha diversity measurements. This indicates that the true richness and evenness of the microbial environment was not fully examined in the present study. Consequently, the results on alpha diversity may be misrepresented [42]. Secondly, the finding that participant is a significant determinant for alpha and beta diversity suggests that the individual microbiome requires separate assessment alongside the explanatory factors during RBAC supplementation on gut microbiota.

The present research findings revealed differences in alpha and beta diversity of gut microbiota within individual participants. Gut microbial composition presents a profile that is distinct to the individual, with moderating factors such as genetic, environment, and lifestyle [43]. Previous research on twins has indicated heritable common gut microbiota, with cohabitating twins sharing the highest number of strains [44]. However, differences in the microbiota between twins were also shown in previous studies, indicating non-genetic factors are also at play [43,44]. Since genetics and environment are known factors influencing gut microbiota composition, it was not unexpected that the microbial composition between participants significantly differed in the present study. Compounding these factors is the individual lifestyle, which further impacts microbiota diversity among participants [45].

The present results also revealed significant association between alpha (Shannon’s evenness) and beta diversity with alcohol consumption, cigarette smoking, and diet (Australian Recommended Food Score, ARFS). Alcohol consumption has been shown to affect gut microbiota and has been linked to substantial losses in diversity, including reductions in Bacteroidetes and Firmicutes [46]. Alcohol consumption also correlates with decreased connectivity of the microbial network and subsequent alteration of gut microbiome composition [28]. Cigarette smoking causes taxonomic changes in gut microbiota, specifically decreased Bacteroidetes and increased Firmicutes, due to nicotine exposure [47]. Moreover, after smoking cessation, increased alpha diversity of gut microbiota was observed [27]. Due to the toxic effects on gut mucosa, cigarette smoking has also been linked to gut microbial dysbiosis [47]. The present study identified lifestyle including diet and subsequent ARFS, alcohol consumption, and cigarette smoking are important factors that impact gut microbiota in healthy adults.

Five of the six participants in the present study that were categorised in the highest ARFS groups revealed significant changes in beta diversity. Although there is limited research on how the ARFS is associated with the gut microbiome, Aslam et al. [48] reported that higher ARFS correlated with the consumption of a more diverse diet, with differences shown in beta diversity compared to limited diets with lower ARFS. Hence, the association of ARFS with changes in beta diversity may be explained by the differences between the consumption of core foods (e.g. grains and water) and discretionary foods (e.g. fried products and packaged sweets) in the diet [49]. The increased consumption of core foods and low discretionary foods also shows a higher proportion of complex carbohydrates that benefits gut microbiota. This may explain the significant changes observed in beta diversity for participants with a higher ARFS in the present study [49].

Age group as an explanatory factor was revealed to be associated with significant changes in beta diversity. Considerable shifts in the microbiome have previously been shown during infancy, puberty, to the later stage of life (>75 years of age) [21]. The stability of the microbiome is reached during adulthood (18-25 years of age) and remains relatively constant until approximately 75 years of age when a loss of diversity occurs [21]. However, Odamaki et al. [23] observed significant differences between adult clusters at the ages of 33 (cluster 1) and 42 (cluster 2), with higher levels of Bacteroidetes observed in the younger age group (cluster 1). This is consistent with the significant differences in beta diversity observed in the present study for participants ≤30 years and >30 years of age.

Associations between sex, and alpha and beta diversity were observed. Sex is a known modulator of the gut microbiome through the actions of sex hormones [50]. This is supported by recent studies showing that sex is a prime contributor to microbial diversity with increased alpha diversity observed in females compared to males [51,52]. Furthermore, sex has also been linked to responses of gut microbiota to diet, anti-microbial effects, and obesity [50,53]. The present findings affirm that age and sex significantly influence gut microbiota composition within healthy adults.

Dietary RBAC supplementation was shown to change the abundance of beneficial bacteria including Lactobacillus, Roseburia, and Ruminococcaceae, although these changes were not observed in all participants. Some participants exhibited trends in taxonomic changes during the interventional period. For example, participant 9 (P09) showed increases across multiple genera. These beneficial bacteria play a vital role in gut health through the production of SCFAs as essential energy sources for colonic enterocytes [12,14], the provision of vitamins [10], and possess anti-inflammatory properties [13]. This finding suggests that dietary RBAC supplementation may be used as a prebiotic for regulating these beneficial bacteria.

The increases in beneficial bacteria observed in some participants may also be associated with the ARFS. Diets with less diversity have been associated with lower bacterial variety [54]. Of the four participants who exhibited taxonomic changes in beneficial bacteria, three had the lowest ARFS. Hence, this finding suggests a potential association between RBAC supplementation responses and ARFS. This indicates that RBAC supplementation may influence the growth of beneficial gut microbiota of people with diets that are lower in core foods and higher in discretionary foods. In addition to increased beneficial bacteria, there were also reductions in bacteria known to be opportunistic pathogens, including Anaerococcus and Corynebacterium [55,56]. However, these taxonomic changes are difficult to interpret due to the difference in gut microbiota between individuals and how rapidly the gut microbiome can change in response to daily diet and lifestyle factors [45].

In six participants, significant changes in beta diversity were associated with RBAC dosage, experimental phase, and/or interventional period. Previous research has indicated that the composition and diversity of the gut microbiome can determine whether participants respond to interventions including dietary modification for health-related purposes [57]. This is further extended to other interventions such as treatments for cervical cancer and responsiveness to exercise for prediabetes treatment [57,58,59]. The present study used multiple analyses to classify participants into two groups, responders and non-responders to RBAC supplementation [60]. Four participants (P02, P03, P07, and P09) observed significant changes in all three classification criteria of individual beta diversity, weighted unifrac EMPeror plots, and taxonomic analysis, thus making them responders to RBAC supplementation. Two participants (P04 and P08) revealed changes in two criteria, hence they were probable responders to RBAC supplementation, however further supplementation and analysis is required to determine whether these participants can be classified as responders. Two participants (P01 and P10) only showed microbial taxonomic changes, therefore were unlikely to be responders to the supplementation. Further two participants (P05 and P06) showing no statistical changes across the analyses were considered non-responders to the RBAC supplementation.

Given the sample size in the present study, further research with a larger sample size is required. As demonstrated in the present study, research aiming to investigate the effects of dietary RBAC supplementation on gut microbiota may require controlling for individual characteristics as confounding factors, such as age. Also individualised doses according to body weight may be more appropriate to assess the impact of supplementation and dosage on gut microbiota. A further limitation of this study was the short timeframe of the intervention. An extended period of supplementation may reveal further changes in gut microbiota composition. Maintaining a food journal by participants to track the dietary intake and physical activity could be correlated with gut microbiota composition [45].

4. Materials and Methods

4.1. Participant recruitment and diet tracking

Ten healthy adult volunteers were recruited through University presentations and social media posts. The inclusion criteria comprised healthy adults aged 18-60 years old with informed consent. Exclusion criteria encompassed chronic health conditions and reported gut dysbiosis, ongoing medication use, history of antibiotic use in the three months prior to taking the first dietary supplement, and pregnancy or breastfeeding. Participants were required to maintain their usual diet with no major changes (e.g., changes from omnivore to vegan or vegetarian).

Baseline demographic data were collected in a screening survey before study commencement. Dietary intake was tracked by participants completing an online food frequency questionnaire, the AES, before study commencement and at study completion. Data collated from the AES detailed the participants’ macronutrient intake as well as the percentage of discretionary and core food intakes. An ARFS was obtained from survey results for each participant [61]. These scores were categorised into very low (<23), low (24-29), medium (30-33, high (34-40), and very high (>40) levels of core food intakes in the diet.

4.2. Design and intervention

The study intervention required participants to supplement their regular diet with two doses of RBAC, a low dose (1g per day for 6 consecutive weeks) and a high dose (3g per day for 6 consecutive weeks). Each RBAC sachet contained 1g of RBAC as the active ingredient plus microcrystalline cellulose (0.50g), modified starch (0.26g), dextrin (0.20g), and tricalcium phosphate (0.04g) as excipients. The net weight in each sachet amounted to 2g. Participants were instructed to take 1 sachet per day for 6 weeks (low dose), then 3 sachets per day for 6 weeks (high dose), with a washout period of 3 weeks between these doses. Two baseline faecal samples were collected to establish the individual baseline, and one post-intervention faecal sample was collected to assess for ongoing effects. Implementation of this design resulted in a total of 8 time points for faecal sample collection over 21 weeks. This interventional design was implemented as per the schedule shown in Table 6.

4.3. Faecal sample collection, preservation, and DNA extraction

Faecal samples were collected by participants using a faecal swab collection and preservation system (cat. 45670-B, Norgen Biotek Corp Thorold, Ontario, Canada). Nucleic acid within the collected swab remains stable for up to 2 years with DNA stored in the buffer at room temperature [62]. Participants were instructed to swab their faeces at the start of the sampling week as per the sample schedule and place the swab in the provided vial containing buffer. Swabs were either collected or posted from participants and stored in a -80°C freezer. Faecal samples were subjected to DNA extraction using the microbiome DNA isolation kit (cat. 64100, Norgen Biotek Corp, Thorold, Ontario, Canada). For the lysate preparation, swabs were briefly vortexed before 500 $\mu$L of the sample faecal/buffer mixture was transferred to a DNAase-free microcentrifuge tube. Lysis buffers E and A were added as per product procedure and samples were incubated in a 65°C water bath for 10 minutes. Column binding and column washing were performed as per product procedure. DNA elution was completed twice to improve total yield.

4.4. Amplification and sequencing of extracted microbial DNA

The concentration and purity of the extracted DNA were determined for each sample using the Nanodrop 2000 spectrophotometer (ThermoFisher Scientific, Waltham, Massachusetts, USA). DNA concentration and purity (260/280 ratio) were noted for sequencing submission. As per MyTaq Red Mix procedure, PCRs were performed on randomly selected extractions to ensure the extracted DNA was primarily of microbial origin and suitable for downstream analysis. Primers (0341-CCTACGGGNGGCWGCAG and 0785-GACTACHVGGGTATCTAATCC 16S rRNA target primers) were diluted to a 1:5 ratio [63]. A master mix was prepared and PCR was conducted as per MyTaq Red Mix insert procedure and electrophoresis using E-Gel precast Agarose Electrophoresis System 1%. Ten (10) $\mu$L of sample and 20 $\mu$L of a 1-kb DNA ladder (Sigma-Aldrich-D0428) were used, two wells (10 and 11) were loaded as a positive and negative control, respectively.

Thirty (30) $\mu$L of the extracted DNA samples were prepared in a round bottom 96 well microplate, sealed and transported to Ramaciotti Centre for Genomics (UNSW, Sydney, NSW, Australia). DNA amplification and sequencing were performed on the Illumina MiSeq platform, targeting the gene region 16S rRNA V3-V4 using primers 341-forward and 805-reverse.

4.5. Microbiota profiling and statistical analysis

Microbiota profiling and statistical analysis were performed using the QIIME2 platform [64]. Data were imported to QIIME2 following the Casava 1.8 paired end demultiplexed fastq instructions. The DADA2 plugin was used to find amplicon sequences, remove chimeric sequences, and minimise noise created by spurious operational taxonomic units [65]. The trim and truncating parameters (from 5’ – 3’) were run with truncating-length-forward 285 and truncating-length-reverse 220, and used for data analysis due to sequencing quality cut-offs and sequencing depths.

Phylogenetic taxonomy was assigned to the sequencing data using VSEARCH with QIIME-compatible SILVA release 138 SSURef NR99 full-length region sequences and 138 SSURef NR99 full-length region taxonomy [66]. Taxonomic changes were determined using the analysis of the composition of microbiota with bias correction (ANCOM) method, where the statistical significance was determined using a W value [67]. A one-way ANOVA was performed on the data to determine significant changes over time.

Alpha diversity (richness and evenness) and beta diversity (change in compositions over time) analyses were completed using QIIME2 q2-diversity plugin based on a phylogenetic diversity tree. This was generated with a sampling depth of 26430 determined by 56.32% of retained features in 100% of samples. Alpha diversity applying Shannon’s evenness and Faith’s phylogenetic diversity (Faith’s PD) were calculated using Kruskal-Wallis H test, with the data being pooled and ranked from smallest to largest (one-way ANOVA on ranks).

Beta diversity was calculated using four different methods: unweighted unique fraction metric (unifrac), weighted unifrac, Bray Curtis, and Jaccard. Statistical analysis of beta diversity was based on the permutational multivariate analysis of variance (PERMANOVA). Analyses were performed to identify associations of beta diversity with explanatory factors and the RBAC supplementation. The Adonis PERMANOVA test for beta diversity significance in QIIME2 was used to determine significant factors that affect beta diversity through fitting polynomial regression models [68].

Participants were classified as responders to RBAC intervention if they exhibited significant changes across three analysis methods (with p<0.05). This included individual beta diversity analysis, changes to beta diversity visualised with Earth Microbiome Project (EMPeror) plots across experimental phases, and alterations in taxonomic compositions. Multiple analyses were performed to reduce the likelihood of false positives due to the limited sample size.

5. Conclusions

Overall, the present study revealed that individual factors such as age and lifestyle significantly influence alpha and beta diversity of human gut microbiota. Dietary supplementation of RBAC was shown to influence the beta diversity of gut bacteria of most participants. Based on individual changes in the gut microbiota profile, participants were classified as responders or non-responders. RBAC supplementation was shown to influence the abundance of the beneficial bacteria, Bifidobacterium, Lactobacillus, Roseburia, Ruminococcaceae, and Faecalibacterium in five participants. These findings suggest that fermentation and degradation of RBAC by gut microbiota may be highly individualised. This evidence advocates further research controlling for participant characteristics such as age and lifestyle during dietary supplementation on the gut microbiota in healthy adults.