Gut Microbiota Community and Its Assembly Associated with Age and Diet in Chinese Centenarians
Gut Microbiota Community and Its Assembly Associated with Age and Diet in Chinese Centenarians
Journal of Microbiology and Biotechnology. 2015. Sep, 25(8): 1195-1204
Copyright © 2015, The Korean Society For Microbiology And Biotechnology
  • Received : October 07, 2014
  • Accepted : April 01, 2015
  • Published : September 28, 2015
Export by style
Cited by
About the Authors
Fang, Wang
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Ting, Yu
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Guohong, Huang
Food and Biological Technology Department, Guangxi Polytechnic, Nanning 530226, P.R. China
Da, Cai
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Xiaolin, Liang
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Haiyan, Su
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Zhenjun, Zhu
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Danlei, Li
College of Life Sciences, Tianjin Normal University, Tianjin 300387, P.R. China
Yang, Yang
College of Life Science and Technology, Guangxi University, Nanning 530004, P.R. China
Peihong, Shen
College of Life Science and Technology, Guangxi University, Nanning 530004, P.R. China
Ruifeng, Mao
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Lian, Yu
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Mouming, Zhao
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China
Quanyang, Li
College of Light Industry and Food Engineering, Guangxi University, Nanning 530004, P.R. China

Increasing evidence suggests that gut microbiota underpin the development of health and longevity. However, our understanding of what influences the composition of this community of the longevous has not been adequately described. Therefore, illumina sequencing analysis was performed on the gut microbiota of centenarians (aged 100-108 years; RC) and younger elderlies (aged 85-99 years; RE) living in Bama County, Guangxi, China and the elderlies (aged 80-92 years; CE) living in Nanning City, Guangxi, China. In addition, their diet was monitored using a semiquantitative dietary questionary (FFQ 23). The results revealed the abundance of Roseburia and Escherichia was significantly greater, whereas that of Lactobacillus, Faecalibacterium, Parabacteroides, Butyricimonas, Coprococcus, Megamonas, Mitsuokella, Sutterella, and Akkermansia was significantly less in centenarians at the genus level. Both clustering analysis and UniFraq distance analysis showed structural segregation with age and diet among the three populations. Using partial least square discriminate analysis and redundancy analysis, we identified 33 and 34 operational taxonomic units (OTUs) as key OTUs that were significantly associated with age and diet, respectively. Age-related OTUs were characterized as Ruminococcaceae, Clostridiaceae, and Lachnospiraceae, and the former two were increased in the centenarians; diet-related OTUs were classified as Bacteroidales, Lachnospiraceae , and Ruminococcaceae . The former two were deceased, whereas the later one was increased, in the high-fiber diet. The age and high-fiber diet were concomitant with changes in the gut microbiota of centenarians, suggesting that age and high-fiber diet can establish a new structurally balanced architecture of gut microbiota that may benefit the health of centenarians.
Bama Yao Autonomous County, Guangxi, China is the fifth longevous village in the world, as identified at the Thirteenth International Natural Medicine Meeting [1] . According to the Population Census of China in 2010 [34] , centenarians accounted for 36.5 out of 1 × 10 5 people and the longevous phenomenon in the Bama County is more obvious compared with previous reports [40 , 41] . Although longevity is associated with many factors, including heredity, diet, society, and the environment [10] , the gut microbiota is considered a possible determining factor [4 , 14 , 40 , 41] .
Almost inconceivable numbers of microorganisms (10 times greater than the total human eukaryotic cells) inhabit the human intestine [2] , where they anaerobically digest food ingredients undigested in the upper gastrointestinal tract, resulting in improved nutrient absorption [36] , accelerated physiological chemical transformation, enhanced host self-repair, and regulated immunity [2] .
Colonization of the gut microbiota is considered to be deeply rooted in the evolution of human beings [20] . Recent studies on the gut microbiota in human beings have found that their population is associated with numerous factors, including age [40 , 41] , region [40 , 41] , disease [21] , and antibiotic use [12] . Researchers have made efforts to search for gut microbiota structure possibly related with age and region by comparing their abundance among centenarians, younger elderly, and adults, based on culture- or PCR-based techniques, even though they only provided incomplete information on microbial diversity [14 , 40 , 41] . Furthermore, age-related differences were also observed in centenarians living in Italy, using the Human Intestinal Tract Chip (HITChip) and quantitative PCR assays [4] , revealing high amounts of opportunistic pathogens and a relocation of Firmicutes population. These findings suggest that the two main factors affecting the microbiota structure of healthy subjects are age and region, where the latter includes lifestyle and dietary factors. However, only a few studies have explored the association of the gut microbiota in centenarians with age, and even fewer have investigated the correlation of the composition of gut microbiota with diet.
In the current study, illumina sequencing was used to investigate the differences in gut microbiota among centenarians, Bama younger elderly, and Nanning elderly and to pinpoint the association between different ages and diets and a specific combination of populations in the gut microbiota. Our findings may help to identify specific microbiota signatures for Chinese centenarians and their relationships between diet and longevity.
Materials and Methods
- Sampling
A total of 24 subjects were classified according to their age and region as follows: Group RC ( n = 8), aged 100-108 years, from the Jiazhuan village of Bama Yao Autonomous County, Guangxi, China, a Bama suburb; Group RE ( n = 8), aged 85-99 years, living in the Jiazhuan and Xishan village of Bama; and Group CE ( n = 8), the urban elderly, aged 80-92 years, living in Xixiangtang and Langdong district of Nanning City, Guangxi, China (Table S1). The healthy elderly populations were selected from the same socially homogeneous growth environment and eating diet. A detailed clinical examination and medical history were obtained from these elderly’s family members, which showed all participants were healthy. None of participants had had gastrointestinal tract disorders, hypertension, diabetes, or other systemic diseases. All participants provided written informed consent prior to study entering and had not taken antibiotics or probiotics one month before the sampling dates. The study was approved by the Ethics Committee of Bama Yao Autonomous County government of Guangxi Province, China.
Fecal samples were collected in sterile containers with ice box by the subjects themselves or their family members, and immediately stored at -20℃. The samples were then transported to the laboratory by a team member within 1-2 h of sampling and preserved at -80℃ until extraction of genomic DNA.
The habitual long-term dietary information was obtained from all of the participants or their family members using a semiquantitative dietary questionary (FFQ 23) (Table S2).
- DNA Extraction
The genomic DNA was extracted according to the manufacturer’s instructions of the Stool DNA Kit (TIANamp Beijing, China) with slight modification as follows: (i) bacteria cells were split at 95℃ for 10 min; (ii) additional 4 µl of RNase A was added in the stool sample; and (iii) samples were incubated at 70℃ for 30 min. DNA integrity and size were verified by 0.5% agarose gel electrophoresis and DNA concentrations were determined using NanoDrop spectrophotometry (NanoDrop Technologies, USA). The DNA was stored at -20℃.
- PCR Amplification of V4 Region of Bacterial 16S rRNA Gene and Illumina Sequencing
The bacteria genomic DNA was amplified with the 515F and 806R primers specific for the V4 hypervariable regions of the 16S rDNA gene. In order to distinguish each sample and yield accurate phylogenetic and taxonomic information, the reverse primer contained a 6 bp error-correcting barcode unique to each sample. All DNAs were amplified by following a protocol described previously [6] . Twenty-four samples were sequenced on an Illumina MiSeq platform according to the manufacturer’s recommendations, provided by Beijing Novogene Genomics Technology Co. Ltd (China).
- Bioinformatic Analysis
Pairs of reads from the original DNA fragments were merged by using FLASH [25] , which was designed to merge pairs of reads when original DNA fragments were shorter than twice the read length. Sequencing reads were assigned to each sample according to the unique barcode and were analyzed with the QIIME [5] software package (Quantitative Insights Into Microbial Ecology) and UPARSE pipeline [17] . In brief, the reads were filtered by QIIME quality filters using default setting for Illumina processing, and operational taxonomic units (OTUs) were picked using the UPARSE pipeline. A total of 802,799 ± 100,349, 539,472 ± 67,434, and 393,005 ± 49,125 sequences from the stool samples of centenarians, Bama younger elderly, and Nanning elderly, respectively, were selected, and the sequences trimmed to 253 bp and assigned to OTUs at 97% similarity ( Table 2 ). A representative sequence was picked for each OTU and assigned to taxonomic data using the RDP classifier [37] . The diversity and composition of the bacterial communities were determined using the protocol described by Caporaso et al . [6] . In addition, both weighted and unweighted UniFrac were calculated with QIIME to cluster samples using principal coordinate analysis (PCoA) and unweighted pair group method with arithmetic mean (UPGMA).
The average percentage of protein, fat, carbohydrate, and fiber in centenarians, Bama younger elderly, and Nanning elderly.
PPT Slide
Lager Image
aIndicates statistical significance for both group RC and group CE bIndicates statistical significance for both group RE and group CE. cIndicates statistical significance for both group RC and group RE.
Illumina sequencing data summary.
PPT Slide
Lager Image
The OTU number, richness estimates (Chao 1), and diversity estimates (Shannon index) were calculated at 3% distance. aIndicates statistical significance for both group RC and group CE. bIndicates statistical significance for both group RE and group CE.
- Multivariate Statistical Analysis
OTUs were identified using partial least square discriminate analysis (PLS-DA) using Simca-p+11.5. Variable importance in projection (VIP) reflects the variables with the most significant contribution in discriminating. Those OTUs with the value of VIP higher than 1.5 were considered as key OTUs. Correlations between key OTUs and age or diet were examined based on t-value analysis obtained from redundancy analysis (RDA) and preformed using Canoco 4.5 (Biometrics, The Netherlands). Using these key OTUs, a heatmap was constructed with HemI 1.0.
Statistical analyses were performed with SPSS 19.0 for Windows (SPSS Inc., USA) to determine the statistical differences among the three groups using parametric methods (ANOVA) and nonparametric statistical methods (Mann–Whitney U-test and Kruskal–Wallis test).
Results and Discussion
- Characterization of Dietary Habits
To characterize the diet habits of the 24 elders, a survey with food frequency questionnaire (FFQ23) (Tables S2 and S3) was performed. The average percentage of fiber in the Bama diet was 2.3% in RC and 2.3% in RE, significantly higher than that of 1.8% ( p < 0.05) in CE. The average percentage of fat in the diet was 5.5% in RC, significantly higher than that of 3.9% ( p < 0.05) in CE, and the average percentage of carbohydrate in RE was significantly higher than that of 13.5% ( p < 0.05) and 12.6% ( p < 0.05) in RC and CE, respectively ( Table 1 ). Although the average percentage of fat and carbohydrate were different between centenarians and Bama young elderly, compared with that in the urbanized world, the fiber average content was similar among the elderly living in Bama village. In addition, their diets were rich in legumes, cereals, polenta, red meat, and vegetables (dark and light vegetable), all of which were completely self-produced, cultivated and harvested, and considered as high-fiber diet. Separately, the diets of Nanning elderly mainly included dairy products, poultry, egg, fish, cereals, dark vegetables, and fruits, which were purchased from a large local farmers’ market and was considered as urbanized diet.
- Statistical Characteristics of Illumina Sequencing Data
We observed that the gut microbiota of centenarians was more diverse than that of younger elderly. The observed species matrix, which was simplified as the count of the unique OTU numbers, showed an obvious increase in RC and the Chao 1 index was significantly higher in RC and RE ( p < 0.01), although the species diversity in a community, as estimated by the Shannon index, was similar among RC, RE, and CE ( Table 2 ). For most samples, no curves plateaued with the sequencing (Figs. S1 and S2), while the Shannon diversity started to level off, indicating that although new phylotypes could be found in additional sequencing, most diversity had been well captured with the sequencing effort (Fig. S3). Gut microbial community is a mini-ecosystem whose diversity is regarded as a key health indicator of healthy individuals and affected by the health status of the hosts [35] . It has been reported that reduced biodiversity of the gut microbiota is associated with disease status, including obesity and antibiotic treatment [12 , 35] . In contrast, previous studies have suggested that the microbial diversity of gut microbiota was significantly reduced [3 , 4 , 26] in centenarians and younger elderly although they are healthy.
- Distinctive Fecal Microbial Communities Associated with Centenarians
Taxonomically, 14 bacterial phyla were observed in all samples by using Illumina sequencing of the 16S rRNA V4 region ( Fig. 1 ), and more than 94.89% of the sequences in all of the centenarians and younger elderly were found belonging to the three most populated bacterial phyla, namely Firmicutes, Bacteroidetes, and Proteobacteria, in agreement with previous studies showing that such phyla accounted for the majority of gut microbiota in the elderly [4 , 29] . Bacteroidetes was the most predominant phylum in all samples, followed by Firmicutes and Proteobacteria . The relative abundance of the Bacteroidetes was greater in CE than in RE ( p < 0.05), whereas no significant difference was found in other dominant phyla among the three groups.
PPT Slide
Lager Image
Relative abundance of the most abundant bacterial phyla in RC1-RC8 (centenarians), RE1-RE8 (Bama younger elderly), and CE1-CE8 (Nanning elderly).
The relative abundance of Enterobacteriaceae, S24-7, Methanobacteriaceae, and Comamonadaceae was found to be significantly greater in group RC than group CE, whereas the other seven families Verrucomicrobiaceae, Veillonellaceae, Rikenellaceae, Porphyromonadaceae, Barnesiellaceae, Odoribacteraceae, and Alcaligenaceae were significantly lower in group RC than group CE ( p < 0.05; Table 3 ). In addition, whereas the relative abundance of Verrucomicrobiaceae, Rikenellaceae, Lactobacillaceae, and Barnesiellaceae were significantly greater in group RE than in group RC, no significant difference was found in other dominant families between the two groups ( p < 0.05).
List of families that were significantly different among RC, RE, and CE.
PPT Slide
Lager Image
Abbreviations: RC, Centenarians; RE, Bama younger elderly; CE, Nanning elderly. aIndicates statistical significance for both group RC and group CE. bIndicates statistical significance for both group RE and group CE. cIndicates statistical significance for both group RC and group RE.
The relative abundance of two genera, Roseburia and Escherichia , was found to be significantly higher in RC than CE ( p < 0.05; Table 4 ), whereas the other seven genera Parabacteroides, Butyricimonas, Coprococcus, Megamonas, Mitsuokella, Sutterella, and Akkermansia were significantly lower in RC than CE ( p < 0.05), and Lactobacillus, Coprococcus, Faecalibacterium, Mitsuokella, and Akkermansia were significantly lower in RC than RE ( p < 0.05). Moreover, Methanobrevibacter, Lactobacillus, Escherichia, and Akkermansia were significantly higher in RE than in CE ( p < 0.05), whereas Butyricimonas, Lachnospira, Megamonas, Mitsuokella, and Sutterella were significantly lower in RE than in CE ( p < 0.05)
List of genera that were significantly different among RC, RE, and CE.
PPT Slide
Lager Image
Abbreviations: RC, Centenarians; RE, Bama younger elderly; CE, Nanning elderly. aIndicates statistical significance for both group RC and group CE. bIndicates statistical significance for both group RE and group CE. cIndicates statistical significance for both group RC and group RE.
We observed that another important characteristic of gut microbiota in centenarians is the structural change in butyrate-producing bacteria in the phylum Firmicutes . Among of them, Clostridium cluster IV (represented by Faecalibacterium ) and Clostridium cluster XIVa (represented by Roseburia ) are considered beneficial for energy metabolism and normal development of colonic epithelial cells [28] , the reduced amounts of which have been noticed in frail elderly and patients with colorectal cancer and liver cirrhosis [7 , 38 , 42] . The abundance of genus Faecalibacterium of Clostridium cluster IV was less in centenarians than in young elderly, which is in accordance with a previous research [4] . In contrast, compared with younger elderly, we observed that the abundance of Roseburia of Clostridium cluster XIVa was greater in centenarians ( p < 0.05). Many researchers have found the importance of Roseburia in the production of butyrate and conjugated linoleic acids [13 , 15] . A previous study also found that the relative high abundance of Clostridium coccoides and Eubacterium rectale of Clostridium cluster XIVa is clearly greater in centenarians than Bama younger elderly based on the analysis using culture technology [40] . However, a significant reduction of bacteria of Clostridium cluster XIVa in Italian centenarians has been noticed [4] . The discrepancy of Clostridium cluster IV and Clostridium cluster XIVa in gut microbiota of centenarians from different geographic provenience may be related to different dietary habits.
We also observed that Bacteroidetes was more commonly detected in centenarians than in young elderly from the same area, although lower levels of Rikenellaceae and Porphyromonadaceae were seen in centenarians ( p < 0.05). Bacteroidetes species in the gut microbiota play an important role in carbohydrate fermentation, catabolism of polysaccharides, bile acids, and cholesterol, and utilization of amino acids and proteins, and could affect human health [32] . Some evidence showed that Bacteroidetes in elderly individuals tends to be dominated, differing from that of younger individuals [8 , 26] , and consistent with our comparison results. However, an age-related reduction in Bacteroidetes was found in Italian elderly and centenarians [4 , 26] . In addition, consistent with the previous results [4 , 14] , Proteobacteria was more commonly detected in centenarians, such as Enterobacteriaceae was more abundant in centenarians who lived in the ancient communities that usually had poor hygiene environments and thus could favor the spread of these species.
The family Methanobacteriaceae was not present in the Nanning elderly, but its level was more abundant in centenarians. The influence of diet on the gut microbiota has been found in previous studies. For example, fermentation of complex polysaccharides results in the production of primarily acetate, propionate, and butyrate, normally in a ratio of 3:1:1 [33] . The high-fiber diet may support Methanobacteria strains, which utilize acetic acids and produce methane in centenarians and in the environment of ancient communities where intestinal archaea can be dominant. Conversely, the relative abundance of Megamonas and Mitsuokell within the family Veillonellaceae was obviously lower in centenarians ( p < 0.05). The enrichment of family Veillonellaceae in gut micobiota, typical small intestine bacteria [24] , may be bowel movements at a slower pace [7] . Hence, the decrease of Veillonellaceae in centenarians suggested their bowel movements at a regular pace. Our results also showed that the level of Akkermansia , a mucindegrading bacterium, was relatively low in centenarians compared with that in younger elderly ( p < 0.05), which was in line with a previous study [9] , but contrary to another study [4] , showing that the level of A. muciniphila was more abundant in aged people. Hence, whether the decrease of Akkermansia in centenarians is correlated to host aging and longevity would require further investigation.
- Structural Comparison of Gut Microbiota among RC, RE, and CE with Multivariate Statistics
Differences in the total bacterial community at the single sample level were assessed by a net separation of the three groups produced by the RC, RE, and CE samples clustering and a complete linkage hierarchical clustering according to their bacterial genera as found by the RDP classifier ( Fig. 2 A). It was noteworthy that the samples could be grouped into two distinct clusters. Each cluster contained subjects of the same dietary habit (ignoring CE1 and RE1), indicating the effect of dietary habit on the gut microbiota community. It was also noteworthy that one of the clusters was separated, because the RE elderly and RC elderly belonged to different age groups. The results indicated a dominant role of age in shaping the composition of gut microbiota. Similar results were also obtained in a PCoAbased UniFraq distance analysis ( Fig. 2 B).
PPT Slide
Lager Image
16S rRNA gene surveys reveal a clear separation of the three elderly populations investigated. (A) Dendrogram obtained with complete linkage hierarchical clustering of the samples from centenarians (RC1-RC8), Bama younger elderly (RE1-RE8), and Nanning elderly (CE1-CE8) based on the total OTUs. (B) Principal coordinate analysis (PCoA) plot based on the weighted UniFrac metric.
Supporting these observations, associations between the structural shift of gut microbiota with different age and diet in human beings, including infants [18 , 19 , 27] , elderly [8] , and residents with different diets [11 , 30] , have been highlighted by several earlier studies. In addition, alterations of gut microbiota also associated with age have been suggested in centenarians in comparison with other adults and younger elderly using culture-dependent [14 , 41] and -independent molecular approaches [4 , 40] .
- Correlations of Relative Abundance of Bacterial Groups with Age and Diet
PLS-DA, a multivariate analysis method, was used to identify the key OTUs responsible for differences in the gut microbiota. A total of 54 key OTUs with a VIP > 1.5 were identified and used to identify the specific bacteria with intrinsic difference in gut microbiota compositions among groups with different age and diet. Redundancy analysis using the key OTUs as the species variable and age and region as the environmental variables showed that both the constrained ordination models by age and diet were significant ( p = 0.002 and p = 0.002).
The RDA indicated that 33 OTUs of gut microbiota in the three groups were significantly correlated with age ( Fig. 3 A; p < 0.05). Nineteen OTUs were more abundant and 14 OTUs were less in RC. Most of these OTUs were classified to the families Ruminococcaceae (12 OTUs), Lachnospiraceae (6 OTUs), and Clostridiaceae (4 OTUs). OTUs from family Ruminococcaceae were identified as the key specific gut microbiota, among which seven OTUs were highly abundant and five OTUs less abundant in centenarians, and OTUs from Clostridiaceae disappeared in younger elders. In addition, there were 25 OTUs that were significantly different between centenarians and younger elderly (Table S4).
PPT Slide
Lager Image
Heatmap constructed using the key OTUs of elders that correlated with age and diet. The order, family, and genus names of the key OTUs are shown on the left at the corresponding position. Individuals are shown at the top of the heatmap. (A) Relative abundances of the 33 OTUs identified as key variables for the differentiation of microbiota of centenarians (RC1-RC8) and younger elderly (RE1-RE8 and CE1-CE8). (B) Relative abundance of the 34 OTUs identified as key variables for the differentiation of microbiota of high-fiber diet (RC1-RC8 and RE1-RE8) and urbanized diet (CE1-CE8).
As shown in Fig. 3 B, the abundance of 34 OTUs were significantly correlated with the diet of the elderly ( p < 0.05). The abundance of 14 of these OTUs was greater, and of 20 was less, in the fiber diet group than in urbanized diet group. Most of these OTUs were distributed across order Bacteroidales (8 OTUs) and families Lachnospiraceae (7 OTUs) and Ruminococcaceae (6 OTUs). The abundance of these selected phylotypes in Bacteroidales and Lachnospiraceae was lower in the high-fiber diet group, but the OTUs of Ruminococcaceae were higher in the high-fiber diet group. Twenty-six of these OTUs were significantly different between the high-fiber diet and urbanized diet groups (Table S5).
The family Ruminococcaceae of Clostridium cluster IV has been known to degrade dietary fiber and produce SCFA [22] , and the higher abundance of Ruminococcaceae in subjects on diets supplemented with resistant starch had also been described by Salonen et al . [31] , which was in line with our current study. Additionally, OTUs classified to Lachnospiraceae of Clostridium cluster XIVa were found to contribute to the structural differences between gut microbiota of different age and diet groups based on the results of Heatmap. The health-promoting functions of Lachnospiraceae include participating in carbohydrate fermentation into short-chain fatty acids, CO 2 , and H 2 , resulting in increasing nutrients for the host and modulating colonic pH [16 , 39] . Although having no significant difference among the three groups, OTUs related to Lachnospiraceae were less detected in high-fiber diet. Hence, the peculiar compositional layout of butyrate-producing bacteria associated with age and diet in centenarians may contribute to modulating the structure of the gut microbiota to a more balanced state.
OTUs closely related to Clostridiaceae were more abundant in the gut microbiota of centenarians. Although the prevalence of Clostridiaceae is positively associated with infancy food allergy, there is still no sufficient evidence about the role of Clostridiaceae in pathogenicity [21] . In addition, in the case of high-fiber diet, Bacteroidetes was found to be less, although it has always been regarded as the most abundant in polysaccharides-rich diet [11] . Bacteria within Bacteroidetes were adept in degrading resistant starch in the human large intestine [23] and dominant species vary in the ability to utilize different types of starch substrates [33] . Therefore, we speculated that the decreased OTUs of Bacteroidetes in high-fiber diet group might originate from different types of dietary carbohydrates.
Our results point to a new balanced gut microbiota architecture of Chinese centenarian induced by different ages and diets. Given the potential key role of a high-fiber diet in mediating longevity, it will be of interest to the general public and become a possible guide of dietary interventions for health and longevity in the future.
The present work was supported by grants from the National Natural Science Foundation of China (No. 31371762), Innovation Project of Guangxi Graduate Education (No. YCBZ2014019), and the “Bagui Scholars Distinguished Professor” Special Project.We thank Dr. Bo Wu (Guangxi University) for providing the experiment platform and guidance for the work.We have no conflicts of interest to declare.
1996 Use of geographic information systems in epidemiology (GIS-Epi). Epidemiol. Bull. 17 1 - 6
Backhed F , Ley RE , Sonnenburg JL , Peterson DA , Gordon JI 2005 Host-bacterial mutualism in the human intestine. Science 307 1915 - 1920    DOI : 10.1126/science.1104816
Biagi E , Candela M , Fairweather-Tait S , Franceschi C , Brigidi P 2012 Ageing of the human metaorganism: the microbial counterpart. Age 34 247 - 267    DOI : 10.1007/s11357-011-9217-5
Biagi E , Nylund L , Candela M , Ostan R , Bucci L , Pini E 2010 Through ageing, and beyond: gut microbiota and inflammatory status in seniors and centenarians. PLoS One 5 1 - 14    DOI : 10.1371/annotation/df45912f-d15c-44ab-8312-e7ec0607604d
Caporaso J G , Kuczynski J , Stombaugh J , Bittinger K , Bushman FD , Costello EK 2010 QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7 335 - 336    DOI : 10.1038/nmeth.f.303
Caporaso JG , Lauber CL , Walters WA , Berg-Lyons D , Lozupone CA , Turnbaugh PJ 2011 Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. USA 108 4516 - 4522    DOI : 10.1073/pnas.1000080107
Chen Y , Yang F , Lu H , Wang B , Chen Y , Lei D 2011 Characterization of fecal microbial communities in patients with liver cirrhosis. Hepatology 54 562 - 572    DOI : 10.1002/hep.24423
Claesson MJ , Cusack S , O’Sullivan O , Greene-Diniz R , de Weerd H , Flannery E 2011 Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc. Natl. Acad. Sci. USA 108 4586 - 4591    DOI : 10.1073/pnas.1000097107
Collado MC , Derrien M , Isolauri E 2007 Intestinal integrity andAkkermansia muciniphila, a mucin-degrading member of the intestinal microbiota present in infants, adults, and the elderly. Appl. Environ. Microbiol. 73 7767 - 7770    DOI : 10.1128/AEM.01477-07
Davinelli S , Willcox DC , Scapagnini G 2012 Extending healthy ageing: nutrient sensitive pathway and centenarian population. Immun. Ageing 9 1 - 7    DOI : 10.1186/1742-4933-9-9
De Filippo C , Cavalieri D , Di Paola M , Ramazzotti M , Poullet JB , Massart S 2010 Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc. Natl. Acad. Sci. USA 107 14691 - 14696    DOI : 10.1073/pnas.1005963107
Dethlefsen L , Huse S , Sogin ML , Relman DA 2008 The pervasive effects of an antibiotic on the human gut microbiota, as revealed by deep 16S rRNA sequencing. PLoS Biol. 6 2383 - 2400    DOI : 10.1371/journal.pbio.0060280
Devillard E , McIntosh FM , Duncan SH , Wallace RJ 2007 Metabolism of linoleic acid by human gut bacteria: different routes for biosynthesis of conjugated linoleic acid. J. Bacteriol. 189 2566 - 2570    DOI : 10.1128/JB.01359-06
Drago L , Toscano M , Rodighiero V , De Vecchi E , Mogna G 2012 Cultivable and pyrosequenced fecal microflora in centenarians and young subjects. J. Clin. Gastroenterol. 46 S81 - S84    DOI : 10.1097/MCG.0b013e3182693982
Duncan SH , Belenguer A , Holtrop G , Johnstone AM , Flint HJ , Lobley GE 2007 Reduced dietary intake of carbohydrates by obese subjects results in decreased concentrations of butyrate and butyrate-producing bacteria in feces. Appl. Environ. Microbiol. 73 1073 - 1078    DOI : 10.1128/AEM.02340-06
Duncan SH , Louis P , Flint HJ 2007 Cultivable bacterial diversity from the human colon. Lett. Appl. Microbiol. 44 343 - 350    DOI : 10.1111/j.1472-765X.2007.02129.x
Edgar RC 2013 UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10 996 - 998    DOI : 10.1038/nmeth.2604
Fan W , Huo G , Li X , Yang L , Duan C 2014 Impact of diet in shaping gut microbiota revealed by a comparative study in infants during the first six months of life. J. Microbiol. Biotechnol. 24 133 - 143    DOI : 10.4014/jmb.1309.09029
Ku H-J , Lee J-H 2014 Development of a novel long-range 16S rRNA universal primer set for metagenomic analysis of gastrointestinal microbiota in newborn infants. J. Microbiol. Biotechnol. 24 812 - 822
Ley RE , Hamady M , Lozupone C , Turnbaugh PJ , Ramey RR , Bircher JS 2008 Evolution of mammals and their gut microbes. Science 320 1647 - 1651    DOI : 10.1126/science.1155725
Ling ZX , Li ZL , L iu X , Cheng YW , Luo YQ , Tong XJ 2014 Altered fecal microbiota composition associated with food allergy in infants. Appl. Environ. Microbiol. 80 2546 - 2554    DOI : 10.1128/AEM.00003-14
Louis P , Duncan SH , McCrae SI , Millar J , Jackson MS , Flint HJ 2004 Restricted distribution of the butyrate kinase pathway among butyrate-producing bacteria from the human colon. J. Bacteriol. 186 2099 - 2106    DOI : 10.1128/JB.186.7.2099-2106.2004
Macfarlane GT , Englyst HN 1986 Starch utilization by the human large intestinal microflora. J. Appl. Bacteriol. 60 195 - 201    DOI : 10.1111/j.1365-2672.1986.tb01073.x
Macfarlane S , Furrie E , Macfarlane GT , Dillon JF 2007 Microbial colonization of the upper gastrointestinal tract in patients with Barrett's esophagus. Clin. Infect. Dis. 45 29 - 38    DOI : 10.1086/518578
Magoc T , Salzberg SL 2011 FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27 2957 - 2963    DOI : 10.1093/bioinformatics/btr507
Mueller S , Saunier K , Hanisch C , Norin E , Alm L , Midtvedt T 2006 Differences in fecal microbiota in different European study populations in relation to age, gender, and country: a cross-sectional study. Appl. Environ. Microbiol. 72 1027 - 1033    DOI : 10.1128/AEM.72.2.1027-1033.2006
Palmer C , Bik EM , DiGiulio DB , Relman DA , Brown PO 2007 Development of the human infant intestinal microbiota. PloS Biol. 5 1556 - 1573    DOI : 10.1371/journal.pbio.0050177
Pryde SE , Duncan SH , Hold GL , Stewart CS , Flint HJ 2002 The microbiology of butyrate formation in the human colon. FEMS Microbiol. Lett. 217 133 - 139    DOI : 10.1111/j.1574-6968.2002.tb11467.x
Rajilic-Stojanovic M , Heilig HGHJ , Molenaar D , Kajander K , Surakka A , Smidt H , de Vos WM 2009 Development and application of the human intestinal tract chip, a phylogenetic microarray: analysis of universally conserved phylotypes in the abundant microbiota of young and elderly adults. Environ. Microbiol. 11 1736 - 1751    DOI : 10.1111/j.1462-2920.2009.01900.x
Ruengsomwong S , Korenori Y , Sakamoto N , Wannissorn B , Nakayama J , Nitisinprasert S 2014 Senior Thai fecal microbiota comparison between vegetarians and non-vegetarians using PCR-DGGE and real-time PCR. J. Microbiol. Biotechnol. 24 1026 - 1033    DOI : 10.4014/jmb.1310.10043
Salonen A , Lahti L , Salojarvi J , Holtrop G , Korpela K , Duncan SH 2014 Impact of diet and individual variation on intestinal microbiota composition and fermentation products in obese men. ISME J. 8 2218 - 2230    DOI : 10.1038/ismej.2014.63
Salyers AA 1984 Bacteroides of the human lower intestinal tract. Annu. Rev. Microbiol. 38 293 - 313    DOI : 10.1146/annurev.mi.38.100184.001453
Scott KP , Gratz SW , Sheridan PO , Flint HJ , Duncan SH 2013 The influence of diet on the gut microbiota. Pharmacol. Res. 69 52 - 60    DOI : 10.1016/j.phrs.2012.10.020
Shi S , Liu S , Huang Q , Huang F , Wei Y 2012 Guangxi Zhuang Autonomous County’s 2010 Population Census China Statistics Press Beijing 222 - 240
Turnbaugh PJ , Hamady M , Yatsunenko T , Cantarel BL , Duncan A , Ley RE 2009 A core gut microbiome in obese and lean twins. Nature 457 480 - 487    DOI : 10.1038/nature07540
van Hylckama Vlieg JE , Veiga P , Zhang C , Derrien M , Zhao L 2011 Impact of microbial transformation of food on health-from fermented foods to fermentation in the gastrointestinal tract. Curr. Opin. Biotechnol. 22 211 - 219    DOI : 10.1016/j.copbio.2010.12.004
Wang Q , Garrity GM , Tiedje JM , Cole JR 2007 Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73 5261 - 5267    DOI : 10.1128/AEM.00062-07
Wang T , Cai G , Qiu Y , Fei N , Zhang M , Pang X 2012 Structural segregation of gut microbiota between colorectal cancer patients and healthy volunteer. ISME J. 6 320 - 329    DOI : 10.1038/ismej.2011.109
Wong JMW , de Souza R , Kendall CWC , Emam A , Jenkins DJA 2006 Colonic health: fermentation and short chain fatty acids. J. Clin. Gastroenterol. 40 235 - 243    DOI : 10.1097/00004836-200603000-00015
Zhao L , Qiao X , Zhu J , Zhang X , Jiang J , Hao Y , Ren F 2011 Correlations of fecal bacterial communities with age and living region for the elderly living in Bama, Guangxi, China. J. Microbiol. 49 186 - 192    DOI : 10.1007/s12275-011-0405-x
Zhao L , Xu W , Ibrahim SA , Jin J , Feng J , Jiang J 2011 Effects of age and region on fecal microflora in elderly subjects living in Bama, Guangxi, China. Curr. Microbiol. 62 64 - 70    DOI : 10.1007/s00284-010-9676-4
Zwielehner J , Liszt K , Handschur M , Lassl C , Lapin A , Haslberger AG 2009 Combined PCR-DGGE fingerprinting and quantitative-PCR indicates shifts in fecal population sizes and diversity ofBacteroides, bifidobacteria andClostridiumcluster IV in institutionalized elderly. Exp. Gerontol. 44 440 - 446    DOI : 10.1016/j.exger.2009.04.002