Corrigendum: Human gut microbiome across different lifestyles: from hunter-gatherers to urban populations
- 1Centro de Investigación en Políticas, Población y Salud, Universidad Nacional Autónoma de México, Mexico City, Mexico
- 2Posgrado en Ciencias Biológicas, Facultad de Medicina, Universidad Nacional Autónoma de México, Mexico City, Mexico
- 3Departamento de Investigación en Tabaquismo y EPOC, Instituto Nacional de Enfermedades Respiratorias Ismael Cosìo Villegas, México City, Mexico
- 4Laboratorio Nacional de Ciencias de la Sostenibilidad (LANCIS), Instituto de Ecología, Universidad Nacional Autónoma de México, Mexico City, Mexico
- 5Departamento de Microbiología y Parasitología, Facultad de Medicina, Universidad Nacional Autónoma de México, Mexico City, Mexico
Human lifestyle and its relationship with the human microbiome has been a line of research widely studied. This is because, throughout human history, civilizations have experienced different environments and lifestyles that could have promoted changes in the human microbiome. The comparison between industrialized and non-industrialized human populations in several studies has allowed to observe variation in the microbiome structure due to the population lifestyle. Nevertheless, the lifestyle of human populations is a gradient where several subcategories can be described. Yet, it is not known how these different lifestyles of human populations affect the microbiome structure on a large scale. Therefore, the main goal of this work was the collection and comparison of 16S data from the gut microbiome of populations that have different lifestyles around the world. With the data obtained from 14 studies, it was possible to compare the gut microbiome of 568 individuals that represent populations of hunter-gatherers, agricultural, agropastoral, pastoral, and urban populations. Results showed that industrialized populations present less diversity than those from non-industrialized populations, as has been described before. However, by separating traditional populations into different categories, we were able to observe patterns that cannot be appreciated by encompassing the different traditional lifestyles in a single category. In this sense, we could confirm that different lifestyles exhibit distinct alpha and beta diversity. In particular, the gut microbiome of pastoral and agropastoral populations seems to be more similar to those of urban populations according to beta diversity analysis. Beyond that, beta diversity analyses revealed that bacterial composition reflects the different lifestyles, representing a transition from hunters-gatherers to industrialized populations. Also, we found that certain groups such as Bacteoidaceae, Lanchospiraceae, and Rickenellaceae have been favored in the transition to modern societies, being differentially abundant in urban populations. Thus, we could hypothesize that due to adaptive/ecological processes; multifunctional bacterial groups (e.g., Bacteroidaceae) could be replacing some functions lost in the transition to modern lifestyle.
Introduction
The human microbiome has been proposed as a fundamental key for maintaining health, including the regulation of the immune system, food digestion, synthesis of vitamins, and protection against pathogens (Selber-Hnativ et al., 2017). Specifically for the gut microbiome, it has been documented that alterations or dysbiosis in the microbial communities have consequences on human health, and have been associated with inflammatory bowel disease, diabetes, and arthritis (Ding et al., 2019). Thus, given the directed or undirected consequences that changes in the microbiome have for human health, it is crucial to understand what factors affect the diversity and composition of the human microbiome in order to get insights into the relationship between human health and microbiomes.
In particular, human lifestyle and its relationship with the gut microbiome have been widely studied (De Filippo et al., 2010; Yatsunenko et al., 2012; Schnorr et al., 2014; Clemente et al., 2015; Obregon-Tito et al., 2015; Jha et al., 2018; Hansen et al., 2019). Throughout human history, civilizations have experienced different environments and lifestyles, which could have led to changes in the microbial communities in the human gut (Jha et al., 2018). For example, the transition from hunter-gatherers to industrialized populations (urbanization) represents important modifications in diet, hygiene practices, health care (including antibiotics), access to resources, and the way that humans interact with the environment (Warinner et al., 2015; Jha et al., 2018). The urbanization process has also promoted globalized supply chains with an increase in the consumption of dairy products, processed foods, and high fat (Jennings et al., 2015). Moreover, it is known that urban populations are less exposed to natural habitats than rural ones, modifying the interaction with microorganisms in the environment (Callewaert et al., 2020). Altogether, the factors involved in the different lifestyles of human populations drive important shifts in the composition, diversity, and structure of the gut microbiome, which in turn have consequences in human health.
Previous studies have found that the gut microbiomes associated with industrialized populations are less diverse and have different microbial compositions than those from traditional lifestyles (Clemente et al., 2015). For instance, it has been shown several times that the microbiomes from human populations living in industrialized areas are dominated by the genus Bacteroides, while hunter-gatherer populations exhibit a high abundance of the genera Prevotella and Treponema (Schnorr et al., 2014; Obregon-Tito et al., 2015). It is important to highlight that such genera are commonly associated with particular functions in the gut microbiome. Specifically, Bacteroides has been associated with high-fat consumption (Singh et al., 2020), whereas Prevotella and Treponema have been linked to high-fiber consumption (De Filippo et al., 2010; Ayeni et al., 2018). Metagenomic analysis also has revealed that the microbial communities associated with industrialized populations present more antibiotic resistance genes than traditional populations (Obregon-Tito et al., 2015; Wibowo et al., 2021). More recent investigations also have found that traditional lifestyles have a higher abundance of mobile genetic elements, probably as adaptations to the seasonal variations (Wibowo et al., 2021).
Although the changes in abundance of microbial groups associated with particular functions seem to be consistent in the studies that compare industrialized and non-industrialized populations (Mancabelli et al., 2017) much remains to be understood between the relationship of the human gut microbiome and human lifestyles. That is, most large-scale studies (meta-analysis) have tried to show the differences between industrialized and non-industrialized populations without taking into account that, within traditional lifestyles such as hunter-gatherer, pastoralist, and agricultural populations, different human practices could be reflected in the composition of the gut microbiome. That is to say, we still don’t know if there are differences in the gut microbiome of populations that carry out activities of domestication of plants, animals, or in the populations that depend on the resources that are available in the environment without an established production system.
Despite the importance of testing how the urbanization process and different lifestyles can impact the composition of the human gut microbiome, no systematic analysis of the reported data exists to date that take into account the different traditional lifestyles as distinct categories (hunter-gatherers, pastoral, agropastoral, agricultural, and urban populations). Thus, we propose a meta-analysis approach that allows us to integrate a large amount of data and identify patterns among different studies that are constructed under the same theoretical framework. The main goal of this meta-analysis was the comparison of the gut microbiome of populations that have different lifestyles across different geographical regions. For this meta-analysis, we analyzed publicly available sequences of 16S rRNA gene from different studies comparing the gut microbiome of human populations. The samples include meta-amplicon sequences of gut microbiomes of hunter-gatherers, pastoral, agropastoral, agricultural, and urban populations.
The collected data allowed us to test if the urbanization process and lifestyle can impact the composition of the gut microbiome. We calculated alpha and beta diversity indices and compared them according to the population’s lifestyle. Also, we identified deferentially abundant microbes between the microbial gut microbiome of human populations with different lifestyles that could be biomarkers. Based on the Phylogenetic Diversity Index we found that, in general, industrialized populations present less diversity than those from non-industrialized populations. Moreover, the PCoA with unweighted UniFrac distance and Beta Dispersion analyses confirm that different lifestyles exhibit distinct beta diversity, representing the transition from hunter-gatherers to urban populations.
Materials and Methods
Data Collection
To compare different lifestyles around the world, raw sequences of 16S rRNA gene of human gut microbiome were downloaded and processed. The analysis included sequences of gut microbiome from hunter-gatherers (Jha et al., 2018; Smits et al., 2018; Hansen et al., 2019; Rubel et al., 2020), pastoral populations (Angelakis et al., 2016; Li H. et al., 2018; Hansen et al., 2019; Rubel et al., 2020), agropastoral populations (Quagliariello et al., 2019; Rubel et al., 2020), traditional agricultural populations (Yatsunenko et al., 2012; Jha et al., 2018; Ruggles et al., 2018) and people who live in urban areas (Yatsunenko et al., 2012; Aly et al., 2016; Angelakis et al., 2016; Bressa et al., 2017; Fujio-Vejar et al., 2017; Girard et al., 2017; McDonald et al., 2018; Quagliariello et al., 2019). Information about geographical region, sex, and age of each individual are available at Supplementary Table 1.
Sample Selection Criteria and Establishment of Lifestyle
To control variables that have been documented as factors that influence the human gut microbiome, exclusion criteria were established for the selection of samples and construction of the database. We took a minimum number of six samples per population with individuals over 18 years old. Individuals reported with metabolic diseases were removed from the analysis. In the case of studies where healthy controls were analyzed against individuals with some type of disease, only healthy individuals were selected. Studies that employed sequencing techniques other than the Illumina platform were also excluded. For this analysis, all 16S gene regions were included (Supplementary Table 1). To categorize human populations lifestyles we used the classification of Food and Agricultural Organization (FAO), and information reported in the individual analyzed studies. We categorized the lifestyles as follows: (a) Hunters-gatherers: people whose main activity is foraging and hunting resources that are around their environments (Galvin and Beeton, 2013), (b) Pastoral: populations whose main activity is raising livestock (Little, 2015; FAO, 2016), (c) Agropastoral: populations that mix small agriculture and livestock (Little, 2015; FAO, 2016), and (d) Agricultural: populations that practice traditional agriculture mainly crop-based and in some cases with small livestock (Edwards et al., 1988).
Data Processing
16S rRNA raw sequences (.fastq files) were downloaded from public repositories (e.g., NCBI, MG-Rast, Github) and analyzed with the Quantitative Insights Into Microbial Ecology (QIIME2) platform (Bolyen et al., 2019) version 2020.2. Processing was done with single-end sequences and individually for each article dataset. Sequences denoising and quality filtering, including primers and adapters trimming, were performed with the DADA2 plug-in (Callahan et al., 2016). We established a cutoff between 100 and 150 nt for all sequences depending of their quality and primers position. The taxonomic assignment of the Amplicon Sequence Variants (ASVs) was carried out independently for each article dataset with the “sklearn-base” classifier using the latest version of Greengenes 13_8 database (McDonald et al., 2012) and with a 99% identity threshold. Samples with less than 7,000 sequences after quality filtering were excluded from the analyses.
Since we analyzed different variable regions of the 16S rRNA gene, the results of each data set were joined in a single feature table and then processed with the Fragment-Insertion QIIME2 plug-in (Janssen et al., 2018) using the reference database sepp-refs-gg-13-8 with default parameters. This plug-in is based on SATé-Enabled Phylogenetic Placement (SEPP) and has been reported as a solution when analyzing different variable regions of the 16S gene (Janssen et al., 2018). The process resulted in a phylogenetic tree that was used to filter the feature table for the subsequent analyses. All analyses related to sequence data processing were performed in QIIME2 (v.2020-2). Finally, all ASVs that match with mitochondria and chloroplast were removed and the number of sequences was standardized by applying a rarefaction of 7,000 (700 permutations) sequences per sample with R’s phyloseq package (McMurdie and Holmes, 2013).
Statistical Analyses
Alpha Diversity Analysis
To compare the gut microbiome diversity of populations living different lifestyles, we calculated the Phylogenetic diversity (PD) and Shannon-Wiener indexes. The indexes were calculated with the “vegan” and “PhyloMeasures” packages in R (Oksanen et al., 2008; Tsirogiannis and Sandel, 2016). Statistical significant differences in alpha diversity between microbial communities were calculated with the Wilcoxon rank-sum test using “vegan” package in R (Oksanen et al., 2019).
Beta Diversity Analysis
For beta diversity, a Principal Coordinate Analysis (PCoA) with unweighted UniFrac was performed at the family level with “phyloseq” R package (McMurdie and Holmes, 2013). Permutation analysis of variance test (PERMANOVA) with unweighted UniFrac matrix (999 permutations) was performed to detect global differences between lifestyles with “vegan” R package. In addition, pairwise differences in beta diversity were analyzed with a multivariate dispersion analysis, using the betadisper function with the unweighted UniFrac distance matrix in “vegan” R package (Oksanen et al., 2019). Then, a pairwise Wilcoxon rank-sum test was employed to detect statistical significant differences in beta dispersion among lifestyles. Moreover, we employed different approaches to detect differentially abundant groups. First, pairwise comparison of differentially abundant groups between lifestyles was carried out with “DESeq2” package (Love et al., 2014), using “poscounts” size factor estimation and default settings. Bacterial groups with a p-adjusted < 0.01 were taken for the analysis (Supplementary Figures 3, 4). Secondly, a Linear Discriminant Analysis (LDA) with effect size (LEfSe) was performed at the family level in the MicrobiomeAnalyst platform (Dhariwal et al., 2017). In LEfSe, only features with an LDA score higher than 2 and a p-adjusted < 0.01 were used for the analysis. The bacterial groups that were differentially abundant in both analyses (DESeq2 and LEfSe) were represented through an alluvial graphic plotted in R. Stacked bar plots of the 20 most abundant families were performed with “phyloseq” R packages (McMurdie and Holmes, 2013).
Finally, in order to corroborate the patterns found at the family level, all the diversity analyses described above were also carried out with the same parameters at the genus level (Supplementary Figures 2, 5–7).
Results
Study Populations
To analyze the gut microbiome of populations living different lifestyles, we search for publicly available sequences of the 16S gene data. According to the exclusion criteria, we collected data from 14 studies with a total of 568 individuals. These data represent three populations of hunters-gatherers (n = 187), four pastoral populations (n = 67), two agropastoral populations (n = 53), six populations of agriculturists (n = 99), and 10 populations that live in urban areas (n = 162). We included in the study the v1-v2, v3-v4, v4, and v4-v5 16S variable regions (Supplementary Table 1).
Alpha Diversity Across Different Lifestyles
Based on the Phylogenetic Diversity index (PD), the gut microbiome of urban populations presented the lowest values of alpha diversity, being significantly less diverse than all other lifestyles (Figure 1). Also, it is important to highlight that within traditional lifestyles there were no significant differences in terms of alpha diversity (Figure 1). In contrast, the Shannon-Winner index showed that hunter-gatherers and agropastoral populations were significantly more diverse than pastoral, agricultural, and urban populations (Supplementary Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/843170/fmicb-13-843170-HTML-r1/image_m/fmicb-13-843170-g001.jpg)
Figure 1. Alpha diversity of human gut microbiome across lifestyles. Boxplot of Phylogenetic diversity analysis (PD) between hunter-gatherers, agricultural, agropastoral, pastoral, and urban populations. Asterisks denote statistical significant differences given by Wilcoxon rank-sum test (*p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001).
Beta Diversity Across Different Lifestyles
Regarding beta diversity, the PCoA with unweighted UniFrac distance at the family level showed a the microbiome structure is slightly different according to lifestyles of hosts and explained almost 40% of the total variation (axis1 = 16.1% and axis2 = 23.6%) (Figures 2A,B). These differences were corroborated with the PERMANOVA test (p = 0.001, R2 = 0.21508, F = 38.5) and the beta dispersion analysis, which showed that all lifestyles were significantly different (Figure 2C). Also in the PCoA figure, we appreciated a shift of diversity from hunter-gatherers to urban populations. For instance, while hunter-gatherers and urban population are in opposite sites of the of the plot, pastoral, agropastoral, and agricultural populations hosts represent intermediate states of microbiome structure. Finally, it is worth to mention that we found the same patterns of beta diversity at genus level (Supplementary Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/843170/fmicb-13-843170-HTML-r1/image_m/fmicb-13-843170-g002.jpg)
Figure 2. Beta diversity analysis of human gut microbiome between lifestyles. (A) Principal Coordinates Analysis (PCoA) with unweighted UniFrac distance at the family level. (B) Boxplots of the distribution of each lifestyle across axis 1 of the PCoA. (C) Boxplot of the Multivariate Beta Dispersion analysis comparison of different lifestyles. Asterisks denote statistical significant differences given by Wilcoxon rank-sum test (*p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001).
On the other hand, both the DEseq2 and LEfse analyses at the family level (Figure 3C and Supplementary Figures 3, 4) revealed some patterns in urban and traditional lifestyles that are summarized in the alluvial plot (Figure 3A) and that were also recovered at the genus level (Supplementary Figures 5–7). For instance, we found that the bacterial families Bacteroidaceae, Bifidobacteriaceae, and Rikenellaceae were differentially abundant in urban and pastoral populations compared to the other traditional lifestyles (Figures 3A–C and Supplementary Figures 3, 4). Contrary, Prevotellaceae, Paraprevotellaceae, Succinivibrionaceae, and Spirochaetaceae, were more abundant in hunter-gatherers, agricultural, and agropastoral populations (Figure 3A). Besides, in LEfse analysis, we found that hunter-gatherers showed a particularly high abundance of Spirochaetaceae while agricultural populations present a high abundance of Succinivibrionaceae and Prevotellacea. Contrary, for urban populations were Bacteroidaceae, Lachnospiraceae, Bifidobacteriaceae, Rikenellaceae, and Enterobacteriaceae (Figure 3B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/843170/fmicb-13-843170-HTML-r1/image_m/fmicb-13-843170-g003.jpg)
Figure 3. Differentially abundant bacterial families between lifestyles. (A) Alluvial plot representing the differentially abundant families among lifestyles according to LEfSe and DEseq2 analyses. (B) Linear Discriminant Analysis with Effect Size (LEfSe) at the family level between lifestyles. Only features with LDA score higher than 2.0 and p < 0.05 were included. (C) Stacked barplot of relative abundance of the top 20 most abundant bacterial families in each lifestyle.
Discussion
Different studies have suggested that human lifestyle is involved in the diversity and composition of the human gut microbiome (De Filippo et al., 2010; Yatsunenko et al., 2012; Schnorr et al., 2014; Clemente et al., 2015; Martínez et al., 2015). However, few studies have tested this hypothesis through systematic collection and analysis of information from populations living in different lifestyles. To determine if, traditional lifestyles and urban populations show previously reported diversity patterns (out of the categorization of industrialized and non-industrialized populations), we performed a meta-analysis focused on collecting information and analyzing 16S data of the human gut microbiome.
From Traditional Lifestyles to Urban Populations
As previously reported, we found that populations with traditional lifestyles exhibit significantly higher alpha diversity than populations living in urban areas or with some degree of industrialization (Figure 1). This finding supports the hypothesis that the urbanization process involves a series of changes related to the modification of the alpha and beta diversity of the human gut microbiome.
Further, by recognizing that non-industrialized lifestyles populations can be separated into sub categories new patterns arise. For instance, we found that some categories such as pastoralist showed less alpha diversity than all other traditional lifestyles (Figure 1). This result has been previously reported in pastoral populations compared to hunter-gatherers and agropastoral populations (Hansen et al., 2019). Historically, pastoralism was practiced in regions where the lands were poor to carry out other activities like agriculture, therefore animal products were considered as the main source of food (e.g., milk, meat, and other animal-derived products) (Little, 2015). It has been shown that a high-fat diet is associated with a decrease in alpha diversity, while a diet rich in fiber is generally associated with increased alpha diversity (Wan et al., 2019). Thus, the limited variety of food resources coupled with a high-fat diet could be influencing the low alpha diversity in pastoral populations, highlighting that a traditional lifestyle does not necessarily lead to a greater diversity of the gut microbiome. Overall, considering the findings described above, it makes sense that populations that are known to have a high-fiber diet like agricultural populations showed the highest values for alpha diversity, whereas populations with higher intakes of animal products like urban and pastoral populations presented the lowest alpha diversity (Figure 1).
On the other hand, by separating the populations into different lifestyles we observed changes in the microbiome composition that could be associated with the transition of lifestyles. Most changes are related with bacterial groups [at family (Figure 3) and genus level (Supplementary Figure 2)] that are known to degrade plant polysaccharides and are associated with a high-fiber diet, such as Prevotellaceae (g_Prevotella), Paraprevotellaceae, Spirochetaceae (g_Treponema), Clostridiaceae, and Succinivibrionaceae (g_Succinivibrio) (De Filippo et al., 2010; Ormerod et al., 2016; Guo et al., 2020; Tett et al., 2021). These polysaccharides-degrading groups were significantly abundant in all traditional lifestyles but less present in pastoral populations and urban populations (Figure 3). In terms of human health, it is important to mention that members of the family Prevotellaceae (g_Prevotella) and Succinivibrionaceae (g_Succinivibrio) are able to ferment complex polysaccharides to produce short-chain fatty acids (SCFAs) that play an important role in colon health (Alemao et al., 2020). Contrary, it seems that the microbiomes associated with urban populations have suffered a depletion of these bacterial families but an enrichment of families correlated with high-fat and high-animal protein diets, such as Bacteroidaceae, Lachnospiraceae, and Rikenellaceae (Okeke et al., 2014; Lecomte et al., 2015; Vacca et al., 2020). It is important to highlight that these groups also have metabolic functions associated with the degradation of complex polysaccharides (Flint and Duncan, 2014; Vacca et al., 2020). However, it has been reported that in the absence of nutrients found in high-fiber food (galacto- and fructo-oligosaccharides, resistant starch), some bacterial groups adapt to obtain energy through digestion of the mucus barrier (Alemao et al., 2020). This is the case of Bacteroidaceae, Lachnospiraceae, and Rikenellaceae that have been reported as main mucosal-sugar generalists (Pereira et al., 2020). These groups probably have to adapt when urban populations transition from a high-fiber diet to a high-fat diet. Nevertheless, alterations of mucosa membranes and the lack of foods rich in fiber have been associated with immune disorders such as inflammatory bowel disease, obesity, and diabetes (Alemao et al., 2020). Moreover, Lachnospiraceae and Bacteroidaceae families have been also related to obesity, colon cancer, and diabetes (Alemao et al., 2020; Vacca et al., 2020), which are considered diseases of western lifestyles (Carrera-Bastos et al., 2011) and nowadays more than a billion people in the world has been affected by metabolic syndrome (Saklayen, 2018).
Opposite to Prevotellaceae, the Bacteroidaceae group seems to have been favored in the transition to modern lifestyles. While in hunter-gatherers Bacteroidaceae is almost absent, in the transition to plant domestication (agricultural populations) the abundance of this group increases (Figure 3A). However, in the transition to animal domestication (pastoral populations), we appreciate a significant increase of such group, culminating in modern societies where Bacteroidaceae was the differentially abundant group (Figure 3B and Supplementary Figure 2). An important characteristic of some members of the Bacteroidaceae (e.g., Bacteroides) and Rikenellaceae family is the capability to tolerate bile, which is linked to higher consumption of animal products (Graf, 2014; Tomova et al., 2019). Additionally, metatranscriptomic analyses have shown that Bacteroidaceae is involved in almost all bacterial functional genes categories in the active gut microbiome (Yap et al., 2014), including the degradation of complex carbohydrates (Flint et al., 2012). In this sense, the Bacteroidaceae group could have increased their abundance due to the loss of bacterial groups that are not tolerant to high bile content, in addition to its capability to perform different metabolic functions in the human intestine. Moreover, members of the Bacteroidaceae (g_Bacteroides) family as well as other genera such as Akkermansia [found highly abundant in urban populations in this study (Supplementary Figure 2)] have been adapted to degrade the mucus barrier in the gut. This probably resulted in a drastic increase of Bacteroidaceae abundance in the transition to modern lifestyle.
Similarly, animal domestication could lead to the increase of other groups such as the Bifidobacterium. In particular, the analysis at the genus level revealed that pastoral and urban populations presented significantly more abundance of this group than all other lifestyles (Supplementary Figure 2). Bifidobacterium is dominant in the gut of healthy breast-fed infants but in adulthood their abundance decreases (Arboleya et al., 2016). Nevertheless, the consumption of dairy products in adulthood has promoted the enrichment of such groups, since various species of Bifidobacterium utilize lactose as their main source of energy (Li X. et al., 2018). Additionally, has been reported that Bifidobacterium can inhibit the growth of Treponema species in the gut (Belkacemi et al., 2020) as we found in the differentially abundant analyses at the genus level (Supplementary Figure 2). Treponema genus is a member of the family Spirochaetaceae that is associated with high-fiber diet (Angelakis et al., 2019). This could explain why the abundance of the family Spirochaetaceae is considerably low in urban and pastoral populations, whereas in hunter-gatherers and agricultural populations the abundance of Spirochaetaceae is higher than Bifidobacteriaceae (Figure 3A). Altogether, these findings could support the antagonistic interaction between Treponema-Bifidobacteirum.
The PCoA analysis performed by lifestyle at both family and genus level (Figure 2A and Supplementary Figure 2) showed that the different hosts’ lifestyles are recovered in the microbiome structure. Particularly, the pattern of change in the microbiome structure seems to recover the transition of the human populations, from the passage of hunter-gatherers to the domestication of plants, followed by the domestication of animals and finally the modern urban populations. However, the observed differences in microbiome composition according to lifestyles (Figure 3C) could be mainly driven by changes in diet composition (e.g., high-fat or high-fiber diets). For example, populations that have a higher consumption of animal-derived products, such as urban, pastoral and agropastoral populations were closer in the PCoA. Contrary, the populations that are known to have a high-fiber diet such as agricultural and hunter-gatherers were closer to each other but more distant to urban populations. Furthermore, we can appreciate that between the categories associated with a high-fat diet (pastoral and agropastoral) there are some differences in the intakes of animal products that seem to be reflected in the PCoA. For instance, pastoral populations that have a higher fat diet compared to other traditional lifestyles were closer or even overlapped with urban populations, followed by agropastoral populations. In the case of agropastoral populations, they have mixed practices between agriculture and pastoralism, which results in higher intakes of foods rich in fiber compared to pastoral or urban populations. That could be the reason why agropastoral populations were shown as an intermediate state between pastoral and agricultural populations in the PCoA. In this context, all the results obtained in this meta-analysis seem to corroborate the hypothesis that, over other variables, diet could have a primary role in shaping the gut microbiome (De Filippo et al., 2010).
Conclusion
Our results support the hypothesis that populations that have undergone a process of urbanization present less diversity in the gut microbiome. According to our meta-analysis based on 16S gene sequences, we could confirm that beyond the industrialized and non-industrialized categories, the different lifestyles established in this study reflect specific structural patterns in the microbial communities of the human gut. Moreover, the major differences between lifestyles were found in the proportion of groups associated with high-fiber and high-fat diets, indicating that diet plays a major role in the composition of human gut microbiome. In particular, the transition from traditional to modern lifestyle seems to impact the gut microbiome by decreasing diversity and promoting the growth of particular microbial groups. That is the case of Bacteroidaceae family which represents one of the most abundant groups in urbanized populations and is known to be involved in almost all bacterial functional genes categories. Thus, we could hypothesize that due to adaptive/ecological processes; multifunctional bacterial groups (e.g., Bacteroidaceae) could be replacing some functions lost in the transition to modern lifestyle. Finally, it is important to mention that the integration of more detailed and precise information of populations and individual samples in future gut microbiome studies is crucial to perform meta-analyses that can integrate and analyze other variables related to the individuals.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
SR-P, AE, RM-E, and RC conceived and designed the project. SR-P and AH-T performed the bioinformatics and statistical analyses. SR-P, AE, RM-E, RC, AH-T, and MN-D wrote the manuscript. RC supervised the project and led the team. All authors discussed the results and commented on the manuscript.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This manuscript is part of the doctoral project of SR-P with the doctoral scholarship supported by Consejo Nacional de Humanidades, Ciencia y Tecnología (CONAHCYT) (CVU: 776199).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
This manuscript is part of the requirements for obtaining a Doctoral degree at the Posgrado en Ciencias Biologicas, UNAM of SR-P. SR-P acknowledges the scholarship supported by Consejo Nacional de Humanidades, Ciencia y Tecnología (CONAHCYT) (CVU: 776199). The authors also acknowledge Rodrigo García Herrera for his support in HPC computing and LANCIS-Instituto de Ecología, UNAM for computing resources.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.843170/full#supplementary-material
Supplementary Table 1 | Metadata associated to the analyzed samples (sample accession, lifestyle, 16S variable region, geographic information, sex, and age).
References
Alemao, C. A., Budden, K. F., Gomez, H. M., Rehman, S. F., Marshall, J. E., Shukla, S. D., et al. (2020). Impact of diet and the bacterial microbiome on the mucous barrier and immune disorders. Allergy Eur. J. Allergy Clin. Immunol. 76, 714–734. doi: 10.1111/all.14548
Aly, A. M., Adel, A., El-Gendy, A. O., Essam, T. M., and Aziz, R. K. (2016). Gut microbiome alterations in patients with stage 4 hepatitis C. Gut Pathog. 8:42. doi: 10.1186/s13099-016-0124-2
Angelakis, E., Bachar, D., Yasir, M., Musso, D., Djossou, F., Gaborit, B., et al. (2019). Treponema species enrich the gut microbiota of traditional rural populations but are absent from urban individuals. New Microbes New Infect. 27, 14–21. doi: 10.1016/j.nmni.2018.10.009
Angelakis, E., Yasir, M., Bachar, D., Azhar, E. I., Lagier, J. C., Bibi, F., et al. (2016). Gut microbiome and dietary patterns in different Saudi populations and monkeys. Sci. Rep. 6, 1–9. doi: 10.1038/srep32191
Arboleya, S., Watkins, C., Stanton, C., and Ross, R. P. (2016). Gut bifidobacteria populations in human health and aging. Front. Microbiol. 7:1204. doi: 10.3389/fmicb.2016.01204
Ayeni, F. A., Biagi, E., Rampelli, S., Fiori, J., Soverini, M., Audu, H. J., et al. (2018). Infant and Adult Gut Microbiome and Metabolome in Rural Bassa and Urban Settlers from Nigeria. Cell Rep. 23, 3056–3067. doi: 10.1016/j.celrep.2018.05.018
Belkacemi, S., Tidjani Alou, M., Million, M., Levasseur, A., Khelaifia, S., and Raoult, D. (2020). Prevalence of Treponema species in the Gut Microbiome is Linked to Bifidobacterium sp. and Bacteroides sp. Res. Sq. [Prepint]. 1–28.
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Bressa, C., Bailén-Andrino, M., Pérez-Santiago, J., González-Soltero, R., Pérez, M., Montalvo-Lominchar, M. G., et al. (2017). Differences in gut microbiota profile between women with active lifestyle and sedentary women. PLoS One 12:e0171352. doi: 10.1371/journal.pone.0171352
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Callewaert, C., Ravard Helffer, K., and Lebaron, P. (2020). Skin Microbiome and its Interplay with the Environment. Am. J. Clin. Dermatol. 21, 4–11. doi: 10.1007/s40257-020-00551-x
Carrera-Bastos, P., Fontes-Villalba, M., O’Keefe, J. H., Lindeberg, S., and Cordain, L. (2011). The western diet and lifestyle and diseases of civilization. Res. Rep. Clin. Cardiol. 2, 15–35. doi: 10.2147/rrcc.s16919
Clemente, J. C., Pehrsson, E. C., Blaser, M. J., Sandhu, K., Gao, Z., Wang, B., et al. (2015). The microbiome of uncontacted Amerindians. Sci. Adv. 1:e1500183. doi: 10.1126/sciadv.1500183
De Filippo, C., Cavalieri, D., Di Paola, M., Ramazzotti, M., Poullet, J. B., Massart, S., et al. (2010). Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc. Natl. Acad. Sci. U.S.A. 107, 14691–14696. doi: 10.1073/pnas.1005963107
Dhariwal, A., Chong, J., Habib, S., King, I. L., Agellon, L. B., and Xia, J. (2017). MicrobiomeAnalyst: a web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. Nucleic Acids Res. 45, 180–187. doi: 10.1093/nar/gkx295
Ding, R. X., Goh, W. R., Wu, R. N., Yue, X. Q., Luo, X., Khine, W. W. T., et al. (2019). Revisit gut microbiota and its impact on human health and disease. J. Food Drug Anal. 27, 623–631. doi: 10.1016/j.jfda.2018.12.012
Edwards, P., Pullin, R. S., and Gartner, J. (1988). Research and Education for the Development of Integrated Crop-Livestock-Fish Farming Systems in the Tropics. (New York, NY: International Center for Living Aquatic Resources Management), 5–9. doi: 10.1016/0308-521x(90)90008-e
FAO (2016). Guidelines for the Enumeration of Nomadic and Semi-Nomadic (Transhumant) Livestock. Food and Agriculture Organization of the United Nations. Available online at: https://www.fao.org/3/ca6397en/ca6397en.pdf (accessed May 22, 2021).
Flint, H. J., and Duncan, S. H. (2014). “Bacteroides and Prevotella,” in Encyclopedia of Food Microbiology, 2nd Edn, eds C. A. Batt and M. L. Tortorello (Oxford: Academic Press), 203–208. doi: 10.1016/B978-0-12-384730-0.00031-8
Flint, H. J., Scott, K. P., Duncan, S. H., Louis, P., and Forano, E. (2012). Microbial degradation of complex carbohydrates in the gut. Gut Microbes 3, 289–306. doi: 10.4161/gmic.19897
Fujio-Vejar, S., Vasquez, Y., Morales, P., Magne, F., Vera-Wolf, P., Ugalde, J. A., et al. (2017). The gut microbiota of healthy Chilean subjects reveals a high abundance of the phylum Verrucomicrobia. Front. Microbiol. 8:1221. doi: 10.3389/fmicb.2017.01221
Galvin, K. A., and Beeton, T. (2013). “Hunter-gatherer societies, ecological impact of,” in Encyclopedia of Biodiversity, 2nd Edn, ed. S. A. Levin (Fort Collins, CO: Elsevier Ltd), 162–165. doi: 10.1016/B978-0-12-384719-5.00072-1
Girard, C., Tromas, N., Amyot, M., and Shapiro, B. J. (2017). Gut Microbiome of the Canadian Arctic Inuit. mSphere 2:e00297-16. doi: 10.1128/msphere.00297-16
Graf, J. (2014). “The Family Rikenellaceae,” in The Prokaryotes: Other Major Lineages of Bacteria and The Archaea, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 857–859. doi: 10.1007/978-3-642-38954-2
Guo, P., Zhang, K., Ma, X., and He, P. (2020). Clostridium species as probiotics: potentials and challenges. J. Anim. Sci. Biotechnol. 11:24. doi: 10.1186/s40104-019-0402-1
Hansen, M. E. B., Rubel, M. A., Bailey, A. G., Ranciaro, A., Thompson, S. R., Campbell, M. C., et al. (2019). Population structure of human gut bacteria in a diverse cohort from rural Tanzania and Botswana. Genome Biol. 20:16. doi: 10.1186/s13059-018-1616-9
Janssen, S., Mcdonald, D., Gonzalez, A., Navas-molina, J. A., Jiang, L., Xu, Z. Z., et al. (2018). Phylogenetic Placement of Exact Amplicon Sequences. mSystems 3:e00021-18.
Jennings, S., Cottee, J., Curtis, T., and Miller, S. (2015). Food in an Urbanised World: The Role of City Region Food Systems. Available online at: https://edepot.wur.nl/408621 (accessed May 22, 2021).
Jha, A. R., Davenport, E. R., Gautam, Y., Bhandari, D., Tandukar, S., Ng, K. M., et al. (2018). Gut microbiome transition across a lifestyle gradient in Himalaya. PLoS Biol. 16:e2005396. doi: 10.1371/journal.pbio.2005396
Lecomte, V., Kaakoush, N. O., Maloney, C. A., Raipuria, M., Huinao, K. D., Mitchell, H. M., et al. (2015). Changes in gut microbiota in rats fed a high fat diet correlate with obesity-associated metabolic parameters. PLoS One 10:e0126931. doi: 10.1371/journal.pone.0126931
Li, H., Li, T., Li, X., Wang, G., Lin, Q., and Qu, J. (2018). Gut Microbiota in Tibetan herdsmen reflects the degree of urbanization. Front. Microbiol. 9:1745. doi: 10.3389/fmicb.2018.01745
Li, X., Yin, J., Zhu, Y., Wang, X., Hu, X., Bao, W., et al. (2018). Effects of whole milk supplementation on gut microbiota and cardiometabolic biomarkers in subjects with and without lactose malabsorption. Nutrients 10:1403. doi: 10.3390/nu10101403
Little, M. A. (2015). “Pastoralism,” in Basics in Human Evolution, ed. M. P. Muehlenbein (San Antonio, TX: Elsevier Inc), 337–347. doi: 10.1016/B978-0-12-802652-6.00024-4
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Mancabelli, L., Milani, C., Lugli, G. A., Turroni, F., Ferrario, C., van Sinderen, D., et al. (2017). Meta-analysis of the human gut microbiome from urbanized and pre-agricultural populations. Environ. Microbiol. 19, 1379–1390. doi: 10.1111/1462-2920.13692
Martínez, I., Stegen, J. C., Maldonado-Gómez, M. X., Eren, M. A., Siba, P. M., Greenhill, A. R., et al. (2015). The Gut Microbiota of Rural Papua New guineans: composition, diversity patterns, and ecological processes. Cell Rep. 11, 527–538. doi: 10.1016/j.celrep.2015.03.049
McDonald, D., Hyde, E., Debelius, J. W., Morton, J. T., Gonzalez, A., Ackermann, G., et al. (2018). American gut: an open platform for citizen-science microbiome research. mSystems 3:e00031-18. doi: 10.1101/277970
McDonald, D., Price, M. N., Goodrich, J., Nawrocki, E. P., Desantis, T. Z., Probst, A., et al. (2012). An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6, 610–618. doi: 10.1038/ismej.2011.139
McMurdie, P. J., and Holmes, S. (2013). Phyloseq: an R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
Obregon-Tito, A. J., Tito, R. Y., Metcalf, J., Sankaranarayanan, K., Clemente, J. C., Ursell, L. K., et al. (2015). Subsistence strategies in traditional societies distinguish gut microbiomes. Nat. Commun. 6:6505. doi: 10.1038/ncomms7505
Okeke, F., Roland, B. C., and Mullin, G. E. (2014). The Role of the gut Microbiome in the Pathogenesis and Treatment of Obesity. Glob. Adv. Heal. Med. 3, 44–57. doi: 10.7453/gahmj.2014.018
Oksanen, J., Blanchet, F. G., Michael, F., Kindt, R., Legendre, P., McGlinn, D., et al. (2019). Vegan: Community Ecology Package, 1–285. Available online at: https://cran.r-project.org/package=vegan
Oksanen, J., Kindt, R., Legendre, P., O’Hara, B., Simpson, G. L., Solymos, P. M., et al. (2008). The vegan package. Community Ecol. Packag 1–182. Available online at: https://cran.r-project.org/package=vegan (accessed January 21, 2021).
Ormerod, K. L., Wood, D. L. A., Lachner, N., Gellatly, S. L., Daly, J. N., Parsons, J. D., et al. (2016). Genomic characterization of the uncultured Bacteroidales family S24-7 inhabiting the guts of homeothermic animals. Microbiome 4:36. doi: 10.1186/s40168-016-0181-2
Pereira, F. C., Wasmund, K., Cobankovic, I., Jehmlich, N., Herbold, C. W., Lee, K. S., et al. (2020). Rational design of a microbial consortium of mucosal sugar utilizers reduces Clostridiodes difficile colonization. Nat. Commun. 11:5104. doi: 10.1038/s41467-020-18928-1
Quagliariello, A., Di Paola, M., De Fanti, S., Gnecchi-Ruscone, A. G., Martinez-Priego, L., Pérez-Villaroya, D., et al. (2019). Gut microbiota composition in Himalayan and Andean populations and its relationship with diet, lifestyle and adaptation to the high-altitude environment. J. Anthropol. Sci. 96, 189–208. doi: 10.4436/JASS.97007
Rubel, M. A., Abbas, A., Taylor, L. J., Connell, A., Tanes, C., Bittinger, K., et al. (2020). Lifestyle and the presence of helminths is associated with gut microbiome composition in Cameroonians. Genome Biol. 21:122. doi: 10.1186/s13059-020-02020-4
Ruggles, K. V., Wang, J., Volkova, A., Contreras, M., Noya-Alarcon, O., Lander, O., et al. (2018). Changes in the Gut Microbiota of Urban Subjects during an Immersion in the Traditional Diet and Lifestyle of a Rainforest Village. mSphere 3:e00193-18. doi: 10.1128/msphere.00193-18
Saklayen, M. G. (2018). The global epidemic of the metabolic syndrome. Curr. Hypertens. Rep. 20:12. doi: 10.1007/s11906-018-0812-z
Schnorr, S. L., Candela, M., Rampelli, S., Centanni, M., Consolandi, C., Basaglia, G., et al. (2014). Gut microbiome of the Hadza hunter-gatherers. Nat. Commun. 5:3654. doi: 10.1038/ncomms4654
Selber-Hnativ, S., Rukundo, B., Ahmadi, M., Akoubi, H., Al-Bizri, H., Aliu, A. F., et al. (2017). Human gut microbiota: toward an ecology of disease. Front. Microbiol. 8:1265. doi: 10.3389/fmicb.2017.01265
Singh, R. P., Halaka, D. A., Hayouka, Z., and Tirosh, O. (2020). High-Fat Diet Induced Alteration of Mice Microbiota and the Functional Ability to Utilize Fructooligosaccharide for Ethanol Production. Front. Cell. Infect. Microbiol. 10:376. doi: 10.3389/fcimb.2020.00376
Smits, S. A., Leach, J., Sonnenburg, E. D., Gonzalez, C. G., Lichtman, J. S., Reid, G., et al. (2018). Seasonal cycling in the gut microbiome of the Hadza hunter-gatherers of Tanzania. Science 357, 802–806. doi: 10.1126/science.aan4834
Tett, A., Pasolli, E., Masetti, G., Ercolini, D., and Segata, N. (2021). Prevotella diversity, niches and interactions with the human host. Nat. Rev. Microbiol. 19, 585–599. doi: 10.1038/s41579-021-00559-y
Tomova, A., Bukovsky, I., Rembert, E., Yonas, W., Alwarith, J., Barnard, N. D., et al. (2019). The effects of vegetarian and vegan diets on gut microbiota. Front. Nutr. 6:47. doi: 10.3389/fnut.2019.00047
Tsirogiannis, C., and Sandel, B. (2016). PhyloMeasures: a package for computing phylogenetic biodiversity measures and their statistical moments. Ecography 39, 709–714. doi: 10.1111/ecog.01814
Vacca, M., Celano, G., Calabrese, F. M., Portincasa, P., Gobbetti, M., and De Angelis, M. (2020). The controversial role of human gut lachnospiraceae. Microorganisms 8:573. doi: 10.3390/microorganisms8040573
Wan, Y., Wang, F., Yuan, J., Li, J., Jiang, D., Zhang, J., et al. (2019). Effects of dietary fat on gut microbiota and faecal metabolites, and their relationship with cardiometabolic risk factors: a 6-month randomised controlled-feeding trial. Gut 68, 1417–1429. doi: 10.1136/gutjnl-2018-317609
Warinner, C., Speller, C., Collins, M. J., and Lewis, C. M. (2015). Ancient human microbiomes. J. Hum. Evol. 79, 125–136. doi: 10.1016/j.jhevol.2014.10.016
Wibowo, M. C., Yang, Z., Borry, M., Hübner, A., Huang, K. D., Tierney, B. T., et al. (2021). Reconstruction of ancient microbial genomes from the human gut. Nature 594, 234–239. doi: 10.1038/s41586-021-03532-0
Yap, G. C., Hong, P. Y., and Lee, B. W. (2014). Microflora of the Intestine: the natural microflora of humans. Encycl. Food Microbiol. Second Ed. 2, 634–638. doi: 10.1016/B978-0-12-384730-0.00207-X
Keywords: meta-analysis, lifestyles, 16S r RNA, gut microbiome, bacterial diversity
Citation: Rosas-Plaza S, Hernández-Terán A, Navarro-Díaz M, Escalante AE, Morales-Espinosa R and Cerritos R (2022) Human Gut Microbiome Across Different Lifestyles: From Hunter-Gatherers to Urban Populations. Front. Microbiol. 13:843170. doi: 10.3389/fmicb.2022.843170
Received: 25 December 2021; Accepted: 25 February 2022;
Published: 26 April 2022.
Edited by:
Andreas Ernst Schwiertz, MVZ Institute of Microecology, GermanyReviewed by:
Juan Antonio Ugalde, Universidad del Desarrollo, ChileAndrea Quagliariello, University of Padua, Italy
Copyright © 2022 Rosas-Plaza, Hernández-Terán, Navarro-Díaz, Escalante, Morales-Espinosa and Cerritos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: René Cerritos, renecerritos@gmail.com