Inter- and intra-individual variations in seasonal and daily stabilities of the human gut microbiota in Japanese

Relationships between human gut microbiota, dietary habits, and health/diseases are the subject of epidemiological and clinical studies. However, the temporal stability and variability of the bacterial community in fecal samples remain unclear. In this study, middle-aged Japanese male and female volunteers (n = 5 each) without disease were recruited from the Sakura Diet Study. Fecal samples and lifestyle information were collected in every quarter and at each defecation for 7 continuous days. Next-generation sequencing of 16S rDNA and hierarchical clustering showed no time trend and intra-individual differences in both fecal sample sets. Significant inter-individual variations in seasonal and daily fecal sample sets were detected for 24 and 23 out of 39 selected dominant genera (>0.1 % of the total human gut microbiota; occupation rate >85 %), respectively. Intra- to inter-individual variance ratios in 26 and 35 genera were significantly <1.0 for seasonal and daily stabilities. Seasonal variation in fermented milk consumption might be associated with Bifidobacterium composition, but not with Lactobacillus. For most of the dominant genera in the human gut microbiota, inter-individual variations were significantly larger than intra-individual variations. Further studies are warranted to determine the contributions of human gut microbiota to nutritional metabolism, health promotion, and prevention/development of diseases.


Introduction
The human intestinal bacterial community is diverse (more than 1000 species) and varies between individuals (Turnbaugh et al. 2009;Lahti et al. 2014). Human gut microbiota is thought to play important roles in association with dietary habits and diseases such as obesity (Ley et al. 2006), inflammatory bowel disease (Walters et al. 2014), and cancer (Schwabe and Jobin 2013). Human intestinal microbiota has been assessed by using terminal restriction fragment length polymorphism analysis (Nagashima et al. 2003); however, identifying each bacterial species is difficult, as it requires the sequencing of all corresponding fragments in a sample. Recently, the use of next-generation sequencing (NGS) led to new findings on specific genera and species of human intestinal bacteria (Shokralla et al. 2012;Quail et al. 2012) with a relatively higher sensitivity and accuracy with respect to bacterial identification.
In recent case-control studies, significant differences in some intestinal bacteria were found between healthy controls and patients with colorectal cancer (Wang et al. 2012;Wu et al. 2013;Ahn et al. 2013). However, consistent findings between studies have not been reported regarding the risk of colorectal cancer. The subject-specific microbiota is thought to be stable for a long period (Rajilić-Stojanović et al. 2012;Martínez et al. 2013;Faith et al. 2013), whereas the composition of each human intestinal bacterium within the microbiota has been demonstrated to change rapidly with a diet rich in fermented milk (Veiga et al. 2014) and animal-based products (David et al. 2014). Thus, the temporal stability and variance remain unclear, rendering bacterial analysis in epidemiological and clinical studies.
Using two sets of seasonal and daily fecal samples from middle-aged Japanese male and female volunteers (n = 5 each) without disease recruited from the Sakura Diet Study, seasonal and daily stabilities were analyzed by NGS of 16S rDNA and hierarchical clustering. We then examined the relative contribution of the inter-and intra-individual variance to each human intestinal bacterium.

Study subjects
Fecal samples were collected from 10 middle-aged Japanese subjects (five men and five women) from January 2013 to March 2014. The subjects were systematically recruited (based on age, gender, residence, and non-severe diseases such as heart and cerebrovascular diseases, and cancer) as volunteers from participants in the Sakura Diet Study in Shizuoka, located in central Japan. They had lived in the area for at least 1 year. Briefly, we orally explained to them the purpose of the study and obtained their signed informed consent for participation in this study. Lifestyle information, including dietary habits, was collected using a scientifically validated questionnaire (Tokudome et al. 2005). During 1-year interval between the first and the second administration of the food frequency questionnaires, the following samples were systematically collected each season: a 3-day dietary record, blood, urine, saliva, feces, and green tea samples as well as blood pressure, height, and body weight. The study was approved by the Ethics Committee, University of Shizuoka .

Questionnaires on defecation and collection of fecal samples
We asked the subjects to collect fecal samples and fill out questionnaires on defecation times in every quarter (four times, representing seasonal samples), and at all times of defecation on 7 continuous days after the fourth season (more than three times in 1 week, considering ones with constipation). For both sets of samples, soon after defecation, each fecal sample was individually suspended by the study subjects in 100 mM Tris-HCl (pH 9), 40 mM EDTA, 4 M guanidine thiocyanate (protein denaturant to inhibit bacterial growth), and 0.001 % bromothymol, as previously described (Fukuda and Fujita 2014;Shiozaki et al. 2014). To clarify inter-and intra-individual variations of human gut microbiota, we collected seasonal and daily fecal sample sets from the same subjects.

DNA extraction from fecal samples
Fecal solids in the suspension were broken down by using FastPrep 24 Instrument (MP Biomedicals, Santa Ana, CA, USA) with zirconia beads at 5 m/s for 2 min. Bacterial DNA extraction from 200 μL of the suspension was performed by using a Magtration System 12GC (Precision System Science, Japan), with MagDEA DNA 200 (Precision System Science) as a reagent for the automatic nucleic acid extraction.

NGS analysis of bacterial community structure in feces
In this study, a series of representative bacteria in the human gut microbiota was analyzed by previously described NGS of 16S rDNA methods (Takahashi et al. 2014) using the following primers (for the V3-V4 region of 16S rDNA of prokaryotes): forward primer 5′-AATGATACGGCGAC CACCGAGATCTACACXXXXXXXXACACTCTTTCCC TACACGACGCTCTTCCGATCTCCTACGGGNBGC ASCAG-3′, where Xs represent the sample-specific 8-bp barcode sequences (CTCTCTAT, TATCCTCT, GTAAGG AG, ACTGCATA, AAGGAGTA, CTAAGCCT, CGTCT AAT, TCTCTCCG, TCGACTAG, and TTCTAGCT) and reverse primer 5′-CAAGCAGAAGACGGCATACGAGAT Z Z Z Z Z Z Z Z G T G AC T G G AG T T C AG AC G T G T G CTCTTCCGATCTGACTACNVGGGTATCTAATCC-3′, where Zs represent the sample-specific 8-bp barcode sequences (TCGCCTTA, CTAGTACG, TTCTGCCT, GCT  CAGGA, AGGAGTCC, CATGCCTA, GTAGAGAG, CAG  CCTCG, TGCCTCTT, TCCTCTAC, TCATGAGC, and  CCTGAGAT); the underlined sequences represent the PCR primer region (Pro341F and Pro805R). Sequencing was Table 1 Possible confounding variables regarding dietary consumption and defecation conditions for each subject (S1 to 10) a,b NS not significant a Seasonal and daily stabilities were defined as stabilities in the human gut microbiota for 1 year (i.e., 4 times representing each season) and 1 week (i.e., 7 continuous days), respectively b Information on food intake and defecation and fecal samples were collected on defecation times in every quarter during 1 year, and at all defecation times during 1 week. We asked volunteers to provide information about such conditions 1 week and 1 day before their defecation, respectively. Thus, regarding daily stability, the data were summarized as the values for 1 week. The results on defecation represent the median (upper) and inter-quartile range (bottom) based on the data for each season during 1 year and 1 week (4 and at least 3 times) c The stool volume was expressed as the number (e.g., 0.5, 1, and 2) of Japanese standard chicken egg size "S" (i.e., 40-52 g). The data represent the median (upper) and inter-quartile range (bottom) at different times (from 3 to 7 times) d The numbers indicate the following: 1: watery, 2: muddy, 3: soft, 4: banana shape, 5: hard, or 6: very hard e The letters represent the following: A: yellow, B: green ocher, C: brown, D: blackish brown, E: slightly blackish, or F: blackish f Shape/consistency: the number of days with defecation in 1 week g The number of defecation times in 1 week h The data represent the median (upper) and inter-quartile range (bottom)

Seasonal stability
Fermented milk (the beverages or yoghurt) (day/week) conducted using a paired-end and modified to 2 × 300-bp cycle run on an Illumina MiSeq sequencing system (Illumina, San Diego, CA, USA) and MiSeq Reagent Kit version 3 (600 Cycle) chemistry. Paired-end sequencing with read lengths of 301 bp was performed. After demultiplexing, a clear overlap in the paired-end reads was observed.
The method of quality filtering of sequences was as follows: only reads that had quality value (QV) scores of ≥20 for more than 99 % of the sequence were extracted for further analysis.

16S rDNA-based taxonomic analysis
Bacterial identification from sequences was performed using the Metagenome@KIN analysis software (World Fusion, Japan) and the TechnoSuruga Lab Microbial Identification database DB-BA 9.0 (TechnoSuruga Laboratory, Japan). Regarding the composition (%) of human intestinal bacteria, we focused on genera representing >0.1 % of the total human gut microbiota, considering the measurement precision, which was approximately >85 % in this study. Based on the analysis of a series of 20 pooled standard samples (i.e., a mixture of our study samples) measured within one run (for 2 days), intra-assay coefficients of variation were <5.0, <10.0, and <15.0 % for 6, 7, and 9 genera, respectively, out of the 22 selected human intestinal bacteria (>0.2 % of the total human gut microbiota), except for Odoribacter. Inter-assay coefficients of variation based on replicate analyses of a total of 60 pooled standard samples over three runs were <7.0 % (<5.0 % for 13 genera out of these), except for Bifidobacterium and Akkermansia. For the remaining genera (including Collinsella, 0.1-0.2 % of the total human gut microbiota), however, the intraand inter-assay coefficients were >15.0, and 2.3-14.4 %, respectively.

Statistical analyses
Considering the measurement precision, in the stability analyses, hierarchical clustering for the selected 17 human intestinal bacterium (>1.0 % of the total human gut microbiota) was performed using GeneMaths software (Applied Maths, Belgium) to obtain an overview of the similarity between 16S rDNA genomic profiles in the subjects. For clustering, the similarity between the profiles was calculated using Pearson's correlation coefficient (r p ), while linkage was calculated using the Ward method. Each variation in the profiles was visualized by different lengths of branches (i.e., r p ) in the clustering tree for seasonal and daily fecal sample sets. For mean values of each human intestinal bacterium over four seasons (v_ four seasons %), r p (only for >1.0 % of the total human gut microbiota, considering r p based on hierarchical clustering analysis) and Spearman's correlation coefficients (r s ) are shown with each genus for 1-3 continuous days (v_ 1d , v_ 2d , and v_ 3d %, in the order). However, regarding r p and r s for daily fecal sample set, the two mean values of 1-3 and 7 days could not be used as independent variables. Cochran's Q-test was used to assess changes in fermented milk intake. For the selected 39 human intestinal bacteria, the inter-and intra-individual differences from each composition of the human intestinal microbiota were analyzed by one-way repeated measures ANOVA. Using a single imputation method, each mean value was used to replace missing values for subjects who could not provide their fecal samples because of constipation. In case of subjects with two or more defecations in a day, mean value of all defecations in the day was used. Friedman test, as a nonparametric method, was also performed, considering the small sample size, i.e., 40 seasonal and 72 daily fecal samples obtained from the 10 subjects. The total variance was partitioned by ANOVA into the following two sources: (1) inter-individual variance in the subject's stationary composition and (2) residual variance, which measured intra-individual variance in either the seasonal or daily fecal sample set (Ogawa et al. 1999). We also calculated the number of days needed to estimate the true value with 95 % confidence intervals within 10 and 20 % of their true mean (Beaton et al. 1979). These analyses were performed with SPSS version 18 (IBM Corporation, Chicago, IL, USA).

Subject characteristics
Means [standard deviation (SD)] of age and BMI were 37.2 (2.5) and 38.2 (10.0) years and 26.4 (3.5) and 23.1 (1.6) kg/ m 2 in men and women, respectively. Fermented milk intake and defecation conditions of each subject are presented in Table 1. Significant seasonal variation was found for fermented milk intake (times/week) (p = 0.022). Probiotics, prebiotics, or antibiotics were not used by the recruited Fig. 1 Double-hierarchical clustering of the 17 selected dominant genera and individual fecal samples collected for each season. As reference, the scale for the genus (representing >1.0 % of the total human gut microbiota) was attached on the above panel of the figure because it was not being outputted in the original figure. The genera and distances between them are depicted. On the right and left, samples are expressed by combined codes for the 10 subjects (S1-S10) and four seasons (1-4), and distances between the samples are depicted. Using the Ward method, each distance matrix between the genera or the samples is shown as Pearson's correlation coefficient (r p ) volunteers. Regarding daily fecal samples, three subjects had constipation and two subjects had frequent defecation. The fecal volume, shape/consistency, and color did not change over the week.

Inter-and intra-individual variations for each intestinal human bacterium
We grouped each 16S rDNA genomic profile by using hierarchical clustering for the selected 17 genera representing >1.0 % of the total human gut microbiota (Figs. 1, 2). Regarding both seasonal and daily stabilities, no time trend was observed for each subject. Even if they were adjacent, consistent continuous orders (i.e., "season 1-4" or "day 1-7") were not observed. Obviously, r p were greater for daily than for seasonal stability. For both sets of samples, Prevotella and Bacteroides were classified in the two largest clusters, at the right and left, respectively (Figs. 1 and 2). Moreover, except Bifidobacterium, most r s between v_ four seasons and v_ 2d or v_ 3d % were similar or greater than those between v_ four seasons and v_ 1d % (e.g., r s = 0.38_0.89, 0.46_0.62, and 0.32_0.53, and 0.64_0.91, 0.78_0.95, and 0.79_0.95 for v_ 1d to v_ 3d % in Bifidobacterium and Bacteroides, respectively). Most of r p were similar to the corresponding r s , but r p in Bifidobacterium, Prevotella, Megamonas, and Roseburia were apparently greater (data not shown).

Relative contributions of intra-(A) and inter-individual variance (B) in each human intestinal bacterium
Relative contributions of intra-(A) and inter-individual variance (B) in each genus and the coefficients of within-person variance (CV w ) and between-person variance (CV b ) are presented in Tables 4 and 5. We also calculated the number of days of fecal sample collection required to estimate the true composition within 10 and 20 % of their true mean. In seasonal and daily fecal sample sets, A/B ratios for 26 and 35 out of the 39 selected genera were <1.0, except for Bacillus, Megasphaera, Mitsuokella, Sporobacter, Brevundimonas, and Mesorhizobium. In the two fecal sample sets, the mean values of each genus were similar, except for Prevotella, Mitsuokella, and Brevundimonas (1.5 vs. 3.5, 0.1 vs. 0.4, and 0.1 vs. 1.1 %, respectively). Compared with CV b , the values of CV w were apparently greater for the seasonal fecal sample set than for the daily one, and most CV b values were >100 % for both fecal sample sets. The former values in compositions of Bifidobacterium and Dorea were >2 times greater in the seasonal fecal sample set than in the daily one, and the latter values in those of 28 genera were >100 % in both or either fecal sample set. According to relatively greater latter values, it might be difficult to accurately estimate the true means, but significant or greater "differences between subjects" (i.e., inter-individual variance) were supported by the results presented in Tables 2, 3, 4, and 5.

Discussion
In the seasonal and daily fecal sample sets, 39 genera representing >0.1 % of total bacteria were identified in the human gut microbiota and occupied approximately >85 % of the total human gut microbiota. In the two dendrograms for 17 genera that represent >1.0 % of the total human gut microbiota, no time trend of fecal sample collection was found and intra-individual variations were observed in Fig. 2 Double-hierarchical clustering of the 17 selected dominant genera and individual fecal samples collected for 7 continuous days. As reference, the scale for the genus (representing >1.0 % of the total human gut microbiota) was attached on the above panel of the figure because it was not being outputted in the original figure. The genera and distances between them are depicted. On the right and left, samples are expressed by combined codes for the 10 subjects (S1-S10) and 7 days (day 1-7), and distances between the samples are also depicted. Using the Ward method, each distance matrix between the genera or the samples is shown as Pearson's correlation coefficient (r p ) Table 2 Compositions (%) of 39 selected dominant genera in the human gut microbiota of the 10 subjects (S1-S10), based on seasonal stability Phylum Subject (S1-S10) Genus S1 S2 S3 S4 S5 S6 Actinobacteria Bifidobacterium 4.0 ± 4.2 a 8.6 ± 6.9 8.8 ± 3.6 1.6 ± 1.3 2.0 ± 2.8 6.2 ± 5.8 Collinsella Bacteroides 11.4 ± 1.3 24.0 ± 7.2 25.6 ± 4.6 10.4 ± 4.6 17.5 ± 3.3 9.7 ± 3.5 Blautia 9.7 ± 3.5 14.5 ± 1.3 15.2 ± 3.3 14.7 ± 5.1 22.2 ± 3.9 9.8 ± 1.9 Catenibacterium Eubacterium 5.7 ± 0.4 2.9 ± 0.5 4.8 ± 1.3 7.4 ± 1.5 5.4 ± 1.1 4.8 ± 0.8 Faecalibacterium 7.4 ± 2.9 4.7 ± 1.5 2.0 ± 1.2 7.0 ± 1.5 6.7 ± 2.9 5.2 ± 2.2 Ruminococcus 7.1 ± 1.1 6.5 ± 1.9 3.2 ± 1.7 7.7 ± 0.4 9.2 ± 3.0 10.1 ± 0.7 Others 12.9 ± 2.9 9.7 ± 0.5 8.2 ± 1.2 25.7 ± 3.8 9.8 ± 2.4 30.2 ± 4.0    Others 4.0 ± 0.4 7.9 ± 3.5 38.9 ± 7.9 6.9 ± 1.7 0.003 0.215 a The individual mean ± SD was calculated using all values for 1 year b According to one-way repeated measures ANOVA, the residual was defined as "Season." The two variables, "Subject" and "Season," correspond to "inter-" and "intra-" individual variations, respectively c Not detected Table 3 Compositions (%) of 39 selected dominant genera in the human gut microbiota for the 10 subjects (S1-S10), based on daily stability Phylum Subject (S1-S10) Genus S1 S2 S3 S4 S5 S6 Actinobacteria Bifidobacterium 4.0 ± 4.2 a 8.6 ± 6.9 8.8 ± 3.6 1.6 ± 1.3 2.0 ± 2.8 6.2 ± 5.8 Collinsella 5.1 ± 1.6 4.7 ± 1.1 0.1 ± 0.1 0.4 ± 0.6 3.8 ± 0.7 3.7 ± 0.2 Bacteroides 11.4 ± 1.3 24.0 ± 7.2 25.6 ± 4.6 10.4 ± 4.6 17.5 ± 3.3 9.7 ± 3.5 Blautia 9.7 ± 3.5 14.5 ± 1.3 15.2 ± 3.3 14.7 ± 5.1 22.2 ± 3.9 9.8 ± 1.9 Catenibacterium Eubacterium 5.7 ± 0.4 2.9 ± 0.5 4.8 ± 1.3 7.4 ± 1.5 5.4 ± 1.1 4.8 ± 0.8 Faecalibacterium 7.4 ± 2.9 4.7 ± 1.5 2.0 ± 1.2 7.0 ± 1.5 6.7 ± 2.9 5.2 ± 2.2 Ruminococcus 7.1 ± 1.1 6.5 ± 1.9 3.2 ± 1.7 7.7 ± 0.4 9.2 ± 3.0 10.1 ± 0.7 Others 12.9 ± 2.9 9.7 ± 0.5 8.2 ± 1.2 25.7 ± 3.8 9.8 ± 2.4 30.2 ± 4.0  In one-way repeated measures ANOVA, the residual was defined as "Day." The two variables "Subject" and "Day" correspond to "inter-" and "intra-" individual variations, respectively c Not detected Table 4 Relative contributions of intra-and inter-individual variance in 39 selected dominant genera, coefficient of within-person variance (CV w ) and between-person variance (CV b ), and the number of days ("Days") of fecal sample collection required to estimate the true values within 10 and 20 % of their true mean, based on seasonal stability a The inter-individual (A) variation represents variation between individual subjects, and intra-individual (B) variation represents variation in season and residual b The mean of the composition (%) among the 10 subjects (all 40 samples) c The number of days of fecal sample collection required to estimate the values within 10 and 20 % of their true mean with 95 % confidence    seasonal and daily stabilities. However, compared with the intra-individual variation, a greater inter-individual variation was found to be statistically significant for almost all 39 selected dominant genera. This suggest that "differences between subjects" were detected based on their human gut bacterial community structure. Our data indicated that changes in the gut human microbiota composition in term of Bifidobacterium, but not Lactobacillus, corresponded to seasonal variations in fermented milk consumption. Dietary assessment studies indicated that inter-individual variance was smaller than intra-individual variance in term of dietary intake of foods and nutrients (Ogawa et al. 1999;Tokudome et al. 2002). Similarly, among American-European individuals, A/B ratios have been reported to be smaller for dietary intake of macronutrients such as protein and carbohydrates, but not for that of minerals and vitamins. In our study, A/B ratios were also smaller for the 39 selected dominant genera, excluding few of them. Thus, similar to dietary assessments, bacterial community structure in the human gut allows the assessment of the intestinal microbiota composition for such dominant selected genera (and species; see supplemental Tables S1 to S4). However, the number of days of dietary information or fecal sample collection for both assessment of dietary foods/nutrient intake and bacterial community structure were not small to estimate the true values within 10 and 20 % of their true means. Even if r s were adequately higher for some of the human intestinal bacteria, it might be difficult to precisely estimate point values for an individual subject, but not mean values in a group of subjects, because of CV w greater values. This study was performed by using bacterial genus data to assess seasonal and daily stabilities of the human intestinal bacteria, but further investigation is needed to determine the relationships between bacterial species data, dietary habits, and health/diseases.
Large-scale projects such as the Human Microbiome Project and Meta-HIT have made substantial progress toward understanding the symbiotic relationships between the human gut microbiota and their hosts (The Human Microbiome Project Consortium 2012; Qin et al. 2012). According to cross-national cluster analyses, three enterotypes have been identified; i.e., Bacteroides, Prevotella, and Ruminococcus (Arumugam et al. 2011). Bacteroides and Prevotella compositions have been reported to be associated with host meals rich in animal proteins, several amino acids, and saturated fats as well as host diet featuring carbohydrates and monosaccharide (Wu et al. 2011;David et al. 2014). Moreover, higher consumption of fermented milk or administration of probiotics and prebiotics has been demonstrated to increase Bifidobacterium composition (Saito et al. 2002;Matsumoto et al. 2010;Petry et al. 2012). This study also found that changes in Bifidobacterium composition might correspond with seasonal changes in terms of fermented milk consumption, and this is related to greater values of intra-individual variation and greater A/B ratios for Bifidobacterium. However, further studies are needed to clarify the association between fermented milk intake and the proportions of Bifidobacterium and Lactobacillus in the microbiota for each season. In epidemiological studies, therefore, the composition of human intestinal bacteria has been shown to associate with change based on food ingestion or administration of supplements.
Individually, unique compositions of human intestinal bacterium were shown in 24 fecal samples derived from three collections of eight subjects without disease (Nam et al. 2011). Additionally, differences in some of human intestinal bacterium were found between healthy controls and patients with Crohn's disease or ulcerative colitis (Gillevet et al. 2010). Through host inflammation, differences between human intestinal bacteria have been commonly shown to relate to the pathogenesis of cancer of the colon, stomach, breast, gallbladder, esophagus, and pancreas (Sheflin et al. 2014;Ohtani et al. 2014). In an animal study, the mechanism involving inflammation suggests that pathogenic species such as Bacteroides fragilis promote colorectal cancer and that Fusobacterium species are related to tumor progression (Tjalsma et al. 2012;McCoy et al. 2013).
This study has some limitations. The number of subjects in the study was small (n = 10), but 40 and 72 fecal samples were systematically collected to examine seasonal and daily stabilities of human intestinal bacteria, respectively. Bacterial community structure was assessed using the two sets of fecal samples from the same subjects. Inter-and intraindividual variations were examined only for a specific age group with equal gender ratio. However, taking into account the age and gender distribution, we appropriately recruited the subjects. In general, bacterial community structure was evaluated bacterial composition (%) because it was difficult to accurately evaluate "bacterial numbers per fecal dry or wet weight." The data on specific species with >0.1 % of the total human gut microbiota are shown, whereas that on genera (and species) with <0.1 % were not shown, considering the measurement limitation for some of them. Regarding species classification, in this study, 83 species with >0.1 % of the total human gut microbiota were identified (i.e., 75 and 77, including 69 common, in seasonal and daily fecal sample sets, respectively), and a statistically significant greater inter-individual variation was observed for approximately 65 % of these (supplemental Tables S1 to S4). Eighteen identified dominant species within >1.0 % represented >45 % of the total human gut microbiota.
In conclusion, using the two sets of seasonal and daily fecal sample sets, we determined the following indices on human gut bacterial community structure: the intra-and inter-individual variance, the ratios, CV w (%), CV b (%), and the number of days of fecal sample collection required to estimate the compositions within 10 and 20 % of their true means. Compared with the intra-individual variance, a statistically significant greater inter-individual variation ("difference between persons") was found for the 39 selected dominant genera. Our findings can be helpful to interpret human bacterial contribution to the role of human gut microbiota in nutritional metabolism, health promotion, and prevention/development of diseases in epidemiological and clinical studies.