Figures
Abstract
Physical exercise is a tool to prevent and treat some of the chronic diseases affecting the world’s population. A mechanism through which exercise could exert beneficial effects in the body is by provoking alterations to the gut microbiota, an environmental factor that in recent years has been associated with numerous chronic diseases. Here we show that physical exercise performed by women to at least the degree recommended by the World Health Organization can modify the composition of gut microbiota. Using high-throughput sequencing of the 16s rRNA gene, eleven genera were found to be significantly different between active and sedentary women. Quantitative PCR analysis revealed higher abundance of health-promoting bacterial species in active women, including Faecalibacterium prausnitzii, Roseburia hominis and Akkermansia muciniphila. Moreover, body fat percentage, muscular mass and physical activity significantly correlated with several bacterial populations. In summary, we provide the first demonstration of interdependence between some bacterial genera and sedentary behavior parameters, and show that not only does the dose and type of exercise influence the composition of gut microbiota, but also the breaking of sedentary behavior.
Citation: Bressa C, Bailén-Andrino M, Pérez-Santiago J, González-Soltero R, Pérez M, Montalvo-Lominchar MG, et al. (2017) Differences in gut microbiota profile between women with active lifestyle and sedentary women. PLoS ONE 12(2): e0171352. https://doi.org/10.1371/journal.pone.0171352
Editor: Suryasarathi Dasgupta, Harvard Medical School, UNITED STATES
Received: August 22, 2016; Accepted: January 19, 2017; Published: February 10, 2017
Copyright: © 2017 Bressa et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Raw data in fastq format are available from the NCBI Biosample database (accession number: PRJNA350839).
Funding: Ministerio de Economía y Competitividad. Ramón y Cajal Grant 2012_11910. www.mineco.gob.es. Universidad Europea de Madrid 2015UEM06, 2015UEM24 and Cátedra ASISA Grant 2015UEM46. www.universidadeuropea.es. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Sedentary lifestyle is associated with a high incidence of chronic diseases such as cardiovascular disease, cancer and diabetes. Physical exercise is a powerful preventative and treatment intervention that is known to be effective in generating metabolic and immune health benefits [1]. The gut microbiota is essential for processing dietary components and has a major role in shaping the immune system [2]. Accordingly, the role of gut microbiota in determining host health and development of disease has been gaining interest in the last decade; however, whether gut microbiota is a target for the therapeutic benefits of exercise remains unknown [3].
Dysbiosis or imbalance in gut microbiota has been associated with many diseases, among which are ulcerative colitis, Crohn's disease, colon cancer [4,5], metabolic syndrome [6], type I and type II diabetes [7,8], cardiovascular disease [9], allergy, asthma, eczema and autism [10]. While it remains unclear in many cases whether the imbalance in gut microbiota is a cause or a consequence of the disease, interventions using probiotics and/or prebiotics are able to ameliorate symptoms of some diseases [11], indicating that manipulation of gut microbiota could be a viable therapeutic strategy [12]. Accordingly, there is a concerted effort to determine the composition of the core healthy gut microbiota and the factors that modify it [13,14]. Along this line, a very recent study that analyzed the gut microbiota of 1135 individuals by deep sequencing found the influence of gut microbiota on 31 intrinsic factors, 12 diseases, 19 drug groups, 4 smoking categories and 60 dietary factors [14]. Although a total of 207 factors were evaluated in this work, physical exercise as a factor that may influence gut microbiota composition was not evaluated. Several studies in experimental models have addressed the relationship between gut microbiota composition and physical exercise. In 2007, Bäckhed et al. suggested the existence of a microbiota-muscle axis that protects mice from obesity [15]. Consistent with this notion, the bacterial abundance of Clostridiaceae and Bacteroidaeae families and Ruminoccocus genus were found to be negatively associated with blood lactate levels in exercised animals, whereas a positive association was found for Oscillospira genus [16]. Evans et al. observed that changes in the Firmicutes/Bacteroidetes ratio exerted by exercise were inversely proportional to the distance traveled by animals [17]. Furthermore, voluntary wheel running was found to induce changes in the abundance of 2150 taxa and was able to offset the decrease in gut microbial richness in mice induced by exposure to polychlorinated biphenyls [18], and also increased Lactobacillus genus and Blautia coccoides-Eubacterium rectale group in rats [19]. Moreover, it seems that the changes exerted by physical exercise depend on the physiological state of the individual. Accordingly, forced exercise increased microbiota richness in obese, hypertensive and normal rats, which was dependent on their physiological state [16]. Additionally, alterations to the microbiota exerted by exercise in rats on high-fat diet were different to those produced in rats on normal diet [20], and those produced in diabetic mice were different to those produced in their control counterparts [21]. Collectively, these findings indicate that modulation of the microbiota by exercise depends not only on the physiological state of the individual, but also on the diet. Finally, it has been observed that exercise induces more effective changes in the microbiota in juvenile rats than in adult rats [22]. The only human study carried out to date that addresses alterations in microbiota by physical exercise compared the gut microbiota of healthy male controls with that of professional rugby players [21]. A greater diversity of microbiota species was found in rugby players, which would indicate a greater state of health in athletes since loss of diversity has been associated with several diseases [23,24], while increased diversity is associated with a better state of health in older people [25]. Nevertheless, in the human study [21], the large difference between diets of athletes and healthy controls (differences in energy, carbohydrates [both complex and simple], protein, fat and fiber) makes it difficult to isolate the effects of exercise against other divergent factors.
The main objective of the present work was to perform an observational study comparing the composition of the gut microbiota between individuals who did not practice any physical exercise and those who at the very least practiced exercise at the minimun dose recommended by the World Health Organization (WHO) [26]. Because hormonal status and sex are factors that influence gut microbiota, the study was conducted in premenopausal women.
Materials and methods
Subject characteristics
Forty premenopausal women were recruited according to inclusion criteria: age 18–40 years and body mass index (BMI) 20–25 kg/m2. Criteria for physical exercise were as follows: sedentary women who did not perform the minimal exercise established as healthy by WHO, ie 3 days of exercise per week for 30 minutes at a moderate intensity [26]; and for active women, those who performed at least 3 hours of physical exercise per week. Exclusion criteria were any kind of pathology, previous gastrointestinal surgery, antibiotics intake during three months prior to the study, smoking, prebiotics, probiotics, vegetarian or vegan, nutritional or ergogenic complements, pregnancy or lactation. All subjects were Caucasian. The Ethics Committee for Clinical Research of the Ramón y Cajal Hospital (Madrid, Spain) approved the study. Written and informed consent was obtained from all participants. The study is registered in Clinicaltrials.gov with the accession number NCT02901912.
Physical activity and sedentary time determination
Levels of physical activity performed by both groups of women were measured with an ActiSleep V.3.4.2 accelerometer. (Actigraph, Manufacturing Technology Inc., Shalimar, FL). This allowed us to verify in an objective manner the presorting of sedentary and active individuals as it determines body movements in all three axes. Acceleration, energy expenditure, intensity of physical activity and body position was recorded. Results were analyzed with Actilife software (Actilife6, Actigraph). Recordings were performed for 7 days (5 weekdays and 2 weekend days). Participants were instructed to wear the accelerometer on the right wrist all day except when they had a shower or performed pool activities. Data were considered valid when they contained a record of at least five valid days including at least one weekend day, and at least 10 hours of activity were recorded. The time-sampling interval (epoch) was every minute during the 7 days. According to previous reports, a cut-off of <100 counts/minute was chosen to categorize sedentary time, which typically includes activities such as sitting or working quietly (eg, reading or typing on a computer) [27]. A sedentary bout was defined as 1 or more consecutive minutes with less than 100 counts/minute. A sedentary break was considered when there was an interruption in sedentary time of at least 1 minute with values above 100 counts/minute. Physical activity was considered as light when the activity count was 100–1951 counts/minute, and moderate-to-vigorous when counts were 1952–5724 counts/minute according to the methodology of Freedson et al. [28].
Food frequency questionnaire
Dietary pattern characterization of the participants was carried out using a food frequency questionnaire (FFQ) with 97 food items. The FFQ was given to participants to complete the day they donated the fecal sample. Data from the FFQ were analyzed using Dietsource software 3.0 (Novartis, Barcelona, Spain) to obtain the total energy ingested (in kcals) of proteins, fat, carbohydrates, fiber and ethanol intake.
Body composition measurement
Body composition was evaluated by dual-energy X ray absorptiometry (DEXA) (GE Healthcare, Madison, WI). Participants were provided with a gown and were placed in the middle of the table in a supine position. The measures of body composition assessed in the study were as follows: total body fat mass, estimated visceral adipose tissue (VAT), total muscle mass as well as fat and muscle mass distribution in the trunk and extremities. The following indices were calculated using the obtained values: adiposity index (AI) = total fat/height2; muscular mass index (MMI) = total muscle mass/height2 and appendicular muscular mass index (AppMMI) = muscle mass in arms+legs/height2.
Stool collection, DNA extraction and sequencing
Participants were provided with a plastic device to collect stool samples, which were stored at -80°C until extraction. DNA was extracted using the commercial Power Soil DNA Kit (Mobio, CA) with a bead-beating homogenizer (Bullet Blender® Storm, Next Advance, NY). The concentration and purity of DNA was measured using the Quant-iT PicoGreen dsDNA Assay Kit (ThermoFisher Scientific, Waltham, MA).
Sequencing and bioinformatics
Microbiota composition analysis of fecal samples was performed by amplifying the hypervariable regions V3 and V4 of the 16s rRNA gene using the primer pair 5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3' and 5'-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3’ [29] to generate amplicons of 459 bp that were sequenced on a MiSeq Illumina platform (Illumina, San Diego, CA). Sequence outputs were analyzed using the Quantitative Insights into Microbial Ecology (QIIME) program, version 1.9.1 [30], using QIIME default parameters except for split library demultiplexing (phred quality threshold of 20 and better). The 16s pair-end reads were assembled using the script multiple_join_paired_ends.py, which joins forward and reverse demultiplexed reads. The output file was processed for quality filtering by split_libraries_fastq.py. High quality sequences were grouped into Operational Taxonomic Units (OTUs) with a sequence identity threshold of 97%, and taxonomy was assigned by interrogating the high quality sequences with the Greengenes database (13_5). Beta-diversity was evaluated by calculating weighted and unweighted Unifrac distances [31]. To study alpha-diversity, Shannon and Simpson diversity indices and Chao1 rarefaction estimators were calculated. Raw data in fastq format are available with the accession number PRJNA350839in the NCBI Biosample database (https://www.ncbi.nlm.nih.gov/bioproject/350839).
Enterotype classification
Participants were classified into “enterotypes” according to the relative abundance of genera and Partitioning Around Medoids (PAM) using the Jensen-Shannon distance based on described methodology [32]. The results were assessed for the optimal number of clusters using the R package NbClust [33] and the Calinski-Harabasz index [34].
Quantitative PCR analysis
Quantitative PCR (qPCR) analysis was carried out on a CFX Connect™ Real-Time PCR Detection System (BioRad, Barcelona, Spain) using SYBR Green I chemistry (BioRad) in 20 µL reactions containing one microliter (1–10 ng) of DNA template and 200 nM of primers. Cycling parameters were 95 °C for 10 minutes, followed by 40 cycles of 95°C for 15 seconds, 1 minute at the established annealing temperature and 72°C for 45 seconds. Subsequently, melting curve analysis was performed, in which fluorescence was measured as the temperature increased from 50°C to 95°C. For bacterial quantification, standard curves were made using serial dilutions of a known DNA concentration corresponding to a known number of bacteria for each group or bacterial species. The calculation of the percentage for each bacterial species or group was performed considering the quantity of bacteria obtained per gram of stool obtained with a universal primer pair as 100% [35]. For Bifidobacterium longum and Roseburia hominis detection, primers were manually designed considering the conserved regions in the 16s rRNA gene and the hypervariable regions to avoid cross-reaction with phylogenetically close species. To do this, the nucleotide sequences in FASTA format (.fas) were obtained from the NCBI nucleotide database. Primer sets were tested in silico with Testprime 1.0 against the small subunit ribosomal RNA SILVA database [29] and with MFE-primer-2.0 report [36]. Neither primer pair set matched to other bacterial species, while Bifidobacterium longum primers matched to the longum and infantis subspecies. The efficiency of PCR amplification for both primer pairs was calculated from the slope and intercept of a standard curve prepared using serial dilutions of DNA extracted from Bifidobacterium longum (DMS 20219) and Roseburia hominis (DMS 16839) acquired from the Leibniz Institute DSMZ-German Collection of Microorganism and Cell Culture. Primer sequences, product size, annealing temperatures and efficiencies are listed in S1 Table. Melting curves obtained from DNA and fecal samples are shown in S1 Fig.
Fecal water content
To obtain the fecal water content, 500 mg of each fecal sample was weighed and then lyophilized in a freeze-drier (Christ Alpha 1–2 LD, Osterode am Harz, Germany). Water content was expressed as the percentage of weight loss of stool samples.
Enzymatic fecal activity
Fecal enzymatic activity was determined using the API-ZYM system (Biomerieux, Marcy l'Etoile, France) according to the method described by Delgado et al. for fecal samples [37] and the manufacturer´s instructions. The API-ZYM system is a semi-quantitative method that determines 19 enzyme activities: alkaline phosphatase, acid phosphatase, phosphohydrolase, esterase, lipase, lipase-esterase, leucine arylamidase, valine arylamidase, cystine arylamidase, trypsin, α-chymotrypsin, α-galactosidase, β-galactosidase, β-glucuronidase, α-glucosidase, β-glucosidase, N-acetyl-β-glucosaminidase, α-mannosidase and α-fucosidase.
Statistical analysis
Statistical analysis was carried out using QIIME version 1.9.1, SPSS software 20.0 (SPSS, Chicago, IL) and the R statistical package 3.3.1. Variable normal distribution was assessed using the Shapiro-Wilk test. Comparisons of variables between sedentary and active subjects were analyzed by Student's t-test (two-tailed). Those variables that did not meet normality were analyzed using the Mann-Whitney test. Bivariate correlation analyses were assessed using the Spearman correlation coefficient (Spearman Rho) using R, employing diet, body composition physical activity parameters and those genera that were present at least in the 80% of the samples. Linear regression analysis was performed for those microbial taxa that were differentially significant between sedentary and active subjects. The stepwise backward elimination method was employed using the following covariates: diet variables (total energy, protein, carbohydrates, lipids, fiber, ethanol, vegetables, cereals, dairy products, fruits, processed meat, beer and coffee), fecal water, body composition parameters (BMI, percentage of body fat, AI, VAT, MMI and AppMMI) and physical activity parameters (moderate physical activity, vigorous physical activity, average hourly kilocalories expended, total time in sedentary bouts, average time per sedentary bout, maximum time per sedentary bout, total time in sedentary breaks, average time per sedentary break, number of sedentary breaks, maximum time per sedentary break and total energy expenditure). Multiple regression models were built when multiple variables were predictive. Principal Component Analyses (PCoA) of community structure (β-diversity) using the unweighted and weighted distance metric were generated by QIIME [30], visualized by EMPeror [30] and analyzed by permutational multivariate analysis of variance (PERMANOVA) using the script compare_categories.py. Significance was set initially at p<0.05. In multivariate analyses (diet, microbial taxa [family and genus], correlation analysis and fecal enzymatic activities), p values were corrected with the Benjamini-Hochberg false discovery rate.
Results
Participant’s characteristics
Forty women were enrolled in the present study: 19 active (ACT) and 21 sedentary (SED). Anthropometric characteristics were similar between groups and no significant differences in age (ACT = 30.7±5.9 years vs SED = 32.2±8.7 years), weight (ACT = 63.45±11.2 kg vs SED = 61.4±9.0 kg), height (ACT = 165.6±8.65 cm vs SED = 163.7±6.3 cm) and body mass index (BMI) (ACT = 24.4±4.5 vs SED = 22.9±3.0) were found between groups.
Body composition analysis
Analysis of body composition parameters indicated significant differences in almost all adiposity and muscle parameters between sedentary and active women (Table 1).
Dietary habits
Data on annual food consumption were collected with a FFQ and analyzed using Dietsource 3.0. Diets of both sedentary and active women conformed to the parameters of a balanced diet (Table 2), where the appropriate percentage of macronutrients is considered as follows: 10–15% proteins, 50–60% carbohydrates and 30–35% lipids [38]. No significant differences in macronutrient consumption and total energy were found between groups; however, the quantity of fiber ingested by the active group was significantly higher than that of the sedentary group (p = 0.005) (Table 2), while no differences were found for ethanol intake. Regarding the results by food groups, there was a significantly greater consumption of fruits (p = 0.027) and vegetables (p = 0.037) by the active group, while the sedentary group ingested significantly greater amounts of processed meats (p = 0.007) (Table 2).
Physical activity
When physical activity was compared between groups, there was a significant difference in the percentage of moderate physical activity in the active group (29.96±5.61%) versus the sedentary group (22.82±3.64%) (p = 0.039), while there were no significant differences in light physical activity measurements. The average expenditure was higher in the active group (145±31 kcal/h) than in the sedentary group (105±22 kcal/h), and consequently the total energy expenditure was also higher for active women (3432±911 kcals) than for sedentary women (2523±982 kcals) (p = 0.048). No significant differences were found for the remaining determined variables. Differences in sedentary behavior were also reflected in the variables “maximum time per sedentary break on working days” (SED = 691±96 min vs ACT = 515±91 min; p = 0.005), “average time per sedentary break on working days” (SED = 81.15±14.71 min vs ACT = 58.82±15.89 min; p = 0.038), and in the variable “total time in sedentary bouts on working days” (SED = 311±21 min vs ACT = 367±35 min; p = 0.027).
Microbiota metagenomic analysis
The average number of reads per sample was 26,890. Rarefaction curves based on observed species, Chao1, and phylogenetic distance whole tree measures were virtually saturated, suggesting sufficient sequencing depth (S2 Fig). No significant differences between the active and sedentary groups were detected in alpha diversity parameters (Observed, Chao1) or Shannon diversity index, indicating that species richness was similar between both groups (S2 Fig). Regarding β-diversity, no significant differences were detected in the Principal Coordinates Analysis (PCoA) based on unweighted Unifrac distance metrics (Fig 1A) (p = 0.141, pseudo-F = 1.15) or on weighted Unifrac distance metrics (Fig 1B) (p = 0.052, pseudo-F = 2.61). A total of 15 phyla were detected, in order of presence: Bacteroidetes (54%), Firmicutes (44%), Proteobacteria (0.96%), Tenericutes (0.39%), Verrucomicrobia (0.11%), Euryarchaeota (0.08%), Actinobacteria (0.07%), Lentisphaerae (0.06%), Cyanobacteria (0.050%), Spirochaetes (0.04%), Fusobacteria (0.014%), Elusimicrobia (0.009%), Synergistetes (0.007%), kTM7 (0.003%), and Acidobacteria (0.0001%). Acidobacteria (2 subjects), Elusimicrobia (2 subjects) and Spirochaetes (2 subjects) phyla were detected only in sedentary subjects. At the phylum level, no significant differences were observed between the two groups, although there was a trend towards a higher presence of Firmicutes (p = 0.085) and a lower presence of Bacteroidetes (p = 0.076) in active women (Table 3). Likewise, the relative abundance of Firmicutes:Bacteroidetes ratio did not significantly differ between groups (p = 0.115). Significant differences were found in two bacterial families, Barnesiellaceae (p = 0.001) and Odoribacteraceae (p = 0.009), with a higher presence in sedentary women (Table 3). At the genus level, there were significant differences in eleven genera: Bifidobacterium, Barnesiellaceae, Odoribacter, Paraprevotella, Turicibacter, Clostridiales, Coprococcus, Ruminococcus, and two unknown genera of Ruminococcaceae family (Fig 2). Given the importance of some bacterial species in health, the presence of Bifidobacterium longum (Fig 3A), Faecalibacterium prausnitzii (Fig 3B), Roseburia hominis (Fig 3C), Akkermansia muciniphila (Fig 3D) was measured by qPCR. Analyses revealed a more significant abundance of F. prautznnii (p = 0.029) (Fig 3B), R. hominis (p = 0.005) (Fig 3C) and A. muciniphila (p = 0.002) (Fig 3D) in active than in sedentary women.
Principal Coordinates Analysis (PcoA) plots of unweighted (A) and weighted (B) Unifrac distance metrics obtained from sequencing the 16s rRNA gene in fecal samples. Axes represent percentage of data explained by each coordinate dimension.
Only statistical significant genera with an adjusted value p<0.001 are represented.
(A) Bifidobacterium longum, (B) Faecalibacterium prausnitzii, (C) Roseburia hominis, (D) Akkermansia muciniphila. Data were log transformed and analyzed by t-test * p<0.05 **p<0.01.
Enterotype classification
PAM analysis favored partitioning into two clusters (Fig 4A), the Bacteroides and Prevotella enterotypes, whereas the Ruminococcus enterotype was not present in our samples (Fig 4B). Seventy-seven percent of the samples were within the Bacteroides enterotype and the remaining 23% belonged to the Prevotella enterotype, in which 66% of the population was from the active group.
Enterotype classification A) k refers to clustering using the Calinski-Harabasz index B) PcoA plots using the Jensen-Shannon distance.
Correlation of gut microbial taxa, body composition and physical activity parameters
Bivariate correlation analysis was performed to study the associations between different variables and their influence on the microbiota. Body composition parameters did not correlate significantly with the percentage of Firmicutes or Bacteroidetes phyla and the Firmicutes/Bacteroidetes (F/B) ratio. At the family level no correlations were found. At the genus level, Odoribacter was positively associated with BFP, AI and VAT (ρ = 0.43, p = 0.001; ρ = 0.44, p = 0.008; ρ = 0.40; p = 0.018, respectively), whereas Turicibacter was negatively associated with BMI (ρ = -0.41, p = 0.016) and Haemophilus was negatively correlated with BFP, AI and VAT (ρ = -0.40, p = 0.02; ρ = -0.43, p = 0.01; ρ = -0.51, p = 0.002, respectively) (Fig 5). Regarding genera associated with muscular parameters, Faecalibacterium was positively associated with MMI (ρ = 0.421, p = 0.013) and AppMI (ρ = 0.416, p = 0.008) and Coprococcus and Lachnospiraceae unclassified1 were positively correlated with AppMMI (ρ = 0.40, p = 0.002; ρ = 0.43, p = 0.001, respectively) (Fig 5). Interestingly, the minimum time in sedentary bouts was significantly associated with the number of species and Shannon and Simpson indices (Fig 5), whereas light physical activity was negatively associated with both indices (Fig 5). Among other genera, Odoribacter was positively associated with the minimum time per sedentary bout, whereas Paraprevotella correlated with the average time per sedentary bout (ρ = 0.52, p = 0.002), Desulfovibrionaceae unclassified correlated with the maximal time per sedentary break (ρ = 0.53, p = 0.004), and Akkermansia with the maximal time per sedentary bout (ρ = 0.42, p = 0.03) (Fig 5).
Multiple regression analysis
Multiple regression analysis was performed to establish the influence of the different variables to the differences observed at genera level between the active and sedentary groups. None of the variables considered for physical activity and body composition had an impact on Turicibacter, Bifidobacterium, Odoribacter, Ruminococcus and Ruminococcus unclassified2 genera, whereas the unclassified genus of the Barnesiellaceae family and Haemophilus were predicted by the percentage of body fat, and Coprococcus and the Ruminococcae family unclassified1 genus were partially explained by muscular mass parameters. Interestingly, the unclassified genus of the Desulfovibrionaceae family and the Paraprevotella genus were predicted by sedentary behavior (Table 4).
Fecal enzymatic activity
Bacterial enzymes in the intestine are involved in the metabolism of foods and drugs, and changes to the composition of intestinal bacterial populations may result in modifications to the bacterial enzyme profile that could impact host health status. Accordingly, enzymatic activities of fecal samples were compared between the active and sedentary groups. The activity of cysteine aminopeptidase in feces was significantly higher in the active group than in the sedentary group (6.6-fold higher, p = 0.033), whereas no differences were detected for other enzyme activities. In correlation analysis, cysteine aminopeptidase activity was negatively correlated with the presence of Bacteroides (ρ = -0.39, p = 0.013), whereas α-fucosidase activity was positively correlated with the presence of Bifidobacterium and Odoribacter genera (ρ = 0.52, p = 0.0005; ρ = 0.45, p = 0.003, respectively) and alkaline phosphatase was negatively correlated with Desulfovibrio genus (ρ = -0.42, p = 0.007).
Discussion
Physical exercise is a low-cost strategy to prevent a wide variety of diseases. One area where physical exercise may have beneficial effects is the gut microbiota, an environmental factor whose imbalance or dysbiosis is associated with various types of diseases including inflammatory bowel disease (IBD), colon cancer, cardiovascular diseases and diabetes, among others. We sought to test whether exercise practiced at a dose recommended by WHO was sufficient to alter the composition of the gut microbiota in women. Our findings indicate that whereas physical exercise did not produce significant changes to the microbiota diversity/richness, sedentary parameters (i.e., sedentary time and breaks) correlated with microbiota richness. Our results agree with those of Welly et al., who also failed to find an increase in microbiota richness in an experimental model of exercise performed spontaneously by animals [39]. By contrast, an increase in microbiota diversity has been reported in professional rugby players [40], and also in animal models in which exercise was performed voluntarily, at high intensity, or forced [16,18,41]. This discrepancy in findings may reflect different factors, such as the modality of exercise (voluntary or forced), intensity, and time in life when exercise is performed [22,42]. In a study comparing voluntary versus forced exercise on gut microbiota, forced exercise was found to produce more marked effects on microbiota diversity than voluntary exercise [39]. A greater diversity of bacterial species has been related to better health conditions [25], and low diversity has been linked to obesity and several diseases such as diabetes [43], IBD [44], colon cancer [45], and autism [45]. Accordingly, it would be of great interest to establish the type and dose of exercise that would be sufficient to increase microbial diversity in the gut.
A decrease in the Bacteroidetes/Firmicutes ratio has been related to obesity, suggesting that the normalization of this ratio might prevent obesity-associated diseases [46]. The effect of exercise on Bacteroidetes and Firmicutes in animal models is still controversial. Evans et al. detected an increase in Bacteroidetes and a decrease in Firmicutes that was reflected in an increased B/F ratio [15], whereas Kang et al. [18] reported a decrease in Bacteroidetes and an increase in Firmicutes, reflecting a decrease in the B/F ratio. In the present study, no changes were detected in Bacteroidetes and Firmicutes; however, we detected a trend for a lower presence of the Bacteroidetes population in the active group, which is consistent with findings in professional athletes [40]. Based on the study by Arumugam et al. [32], we performed an enterotype classification of our samples and we found that only 23% of the population belonged to the Prevotella enterotype, with twice as many active than sedentary women within this enterotype. The Prevotella enterotype has been linked to long-term consumption of carbohydrate-rich diets, whereas the Bacteroides enterotype has been related to a diet rich in protein and animal fat [47]; however, we did not find any relationship between enterotypes and the carbohydrate or protein content of the participant’s diets. Recently, this enterotype classification scheme has been questioned as statistical analysis of a large number of samples of microbiota communities has shown that aside from these discrete clusters, continuous gradients (enterogradients) of the dominant genera can be found that could be better considered as biomarkers of the microbiota, and might more accurately associate with disease or health markers than a discrete enterotype cluster [48,49].
Among all the genera studied, the abundance of eleven of them was significantly different between the active and sedentary group, with Paraprevotella and an unclassified genus of the Desulfovibrionaceae family specifically associated with sedentarism parameters, while the remaining genera where largely associated with diet parameters. Diet is one of the known factors that most influences gut microbiota and to ascertain its influence in human studies is challenging [14]. Individuals who perform high levels of exercise often have special diets [40]. In our study, the physical exercise load was the minimum recommended by WHO for adults aged 18–64 years, which clearly does not require an extreme diet. Nonetheless, as exercise and diet often go hand in hand, an active lifestyle is frequently associated with a high consumption of fruits and vegetables, whereas sedentarism is associated with the consumption of high-calorie and fatty foods. Indeed, exercise interventions in human populations have resulted in an improvement in diet habits [50–52]. Although the diets were similar in our study regarding total carbohydrates, protein and fat content, significant differences were observed for fiber (higher in the active group) and processed meat (higher in the sedentary group). Moreover, we found an inverse correlation between fat intake and muscle parameters, and between fiber intake and body fat composition, which prevented multiple regression analysis of dietary factors and exercise-related factors together because of collinearity problems. According to studies conducted in rodents, females who performed physical exercise had a lower proportion of the Turicibacteraceae family and Turicibacter genus [17,42]. Cereal consumption has been shown to modify Turicibacter abundance in microbiota in animals [53], and therefore alterations in the Turicibacter genus and its relationship with exercise in humans deserves further research. Somewhat surprisingly, a marked increase in the Haemophilus genus, which includes some species related to pathogenicity, was detected in active women. Whether a sizeable presence of Haemophilus in gut microbiota can have a positive or negative effect on health remains unknown. An increase in the abundance of Haemophilus has been described in patients with multiple sclerosis [54], and has also been associated with colorectal carcinoma in patients with adenoma [55]. Conversely, a decrease in Haemophilus spp. has been linked to the development of type II diabetes [56], and a complete depletion of Haemophilus spp. was detected in patients with rheumatoid arthritis and negatively correlated with the level of serum autoantibodies [57].
The alterations in the levels of an unclassified genus of the Desulfovibrionaceae family and the Paraprevotella genus were partially explained by the disruption of sedentary lifestyle. To our knowledge, this is the first demonstration of a correlation between the presence of certain bacterial taxa and the disruption of sedentary lifestyle. The Paraprevotella genus is one of four comprising the Prevotellaceae family, and is characterized by the production of succinic and acetic acid as major fermentation products. Only two species have been described from this genus and its potential effects on human health are unknown [58]. In line with previous results, the Coprococcus genus was more abundant in active women [42]. A low presence of Coprococcus has been linked to IBD and a lower treatment response in pediatric patients [59], and also to decreased resistance to Campylobacter infection [60], and to an increased severity of atopic disease [61]. Coprococcus is a butyrate-producing genus, which could explain some exercise health-promoting effects [62,63]. An unclassified genus of Barnesiellaceae was also modified by exercise. The implication of these bacterial genera in the intestinal microbiota and their possible role in the human body await further investigation.
Species analysis by qPCR showed that the abundance of R. hominis, A. muciniphila and F. prausnitzii was significantly higher in active women than in their sedentary peers. These species have been related to health-promoting effects. A high abundance of A. muciniphila has been previously described in the microbiota of athletes [40], while low levels have been linked to metabolic disorders (obesity, metabolic syndrome and type II diabetes) in patients with IBD [64]. In animal studies, the presence of A. muciniphila has been associated with the existence of other beneficial bacteria such as Bifidobacterium [65]. Although a similar finding was found in our study, with the Bifidobacterium genus more abundant in the active group, an interventional study would be required to test if there is a co-increment of Bifidobacterium and A. muciniphila as the exercise regimen increases. Both R. hominis and F. prausnitzii have a beneficial effect on health as they produce butyrate, which has a positive impact on intestinal health and lipid metabolism [66,67]. F. prausnitzii also produces metabolites with anti-inflammatory action [68]. As described, some bacterial genera were associated with the break of sedentary behavior. This may be important for health because time in sedentary behavior has been associated with the increased risk of metabolic syndrome, type 2 diabetes and cardiovascular disease [69], whereas breaks in sedentary time have been linked to benefits in triglyceride and glucose levels, waist circumference and body mass index [70].
The determination of fecal enzymatic activities revealed a dramatic (6.6-fold) increase in cysteine aminopeptidase activity in active women versus sedentary women. Bacterial cysteine aminopeptidases releases N-terminal cysteine residues from polypeptides and proteins [71], which can serve as energy sources for bacteria that metabolize the protein or for other bacterial species [72]. While bacterial gut protein metabolism is proportionally less significant than carbohydrate metabolism, it is nontheless important because the reaction products produced, including short-chain fatty acids, ammonia, phenols and dicarboxylic acids, can influence host health [72]. Although little is known about specific enzymatic activities of microbiota bacterial groups, cysteine aminopeptidase activity is widely found in the Lactobacillales order [73]; however, we did not find significant differences at any level within this group of bacteria. Further research is needed to understand how changes in microbiota profiles provoke changes in the metabolic profile of the commensal microbiota.
In conclusion, the results presented here together with previous findings show that physical activity modulates the microbiota profile. The dose of exercise that has been considered in our study is the minimum that is recommended by WHO and can thus be easily implemented by individuals in their daily life. Our results indicate that physical activity performed at low doses but continuously can increase the abundance of health-promoting bacteria (Bifidobacterium spp, R. hominis, A. muciniphila and F. prausnitzii) in the microbiota. Although we failed to find an increase in microbiota richness, we found an inverse association between sedentary parameters and microbiota richness (number of species, and Shannon and Simpson indices), which might indicate that perhaps not only is the dose and mode of exercise important, but also the pattern of exercise, such as breaks in sedentary time, avoiding long periods of inactivity in daily life, to induce changes in gut microbiota.
Supporting information
S1 Fig.
Melting curve dissociation analysis of Bifidobacterium longum (A) and Roseburia hominis (B).
https://doi.org/10.1371/journal.pone.0171352.s001
(TIF)
S2 Fig. Observed species, Chao1, and phylogenetic distance whole tree rarefaction curves.
Plotted red line is the mean±standard deviation of active women samples and blue line is the mean±standard deviation of sedentary women samples.
https://doi.org/10.1371/journal.pone.0171352.s002
(TIF)
S1 Table. Primers and amplification conditions for real-time PCR.
https://doi.org/10.1371/journal.pone.0171352.s003
(DOCX)
Author Contributions
- Conceptualization: ML MP.
- Data curation: CB MBA RGS.
- Formal analysis: ML DM.
- Funding acquisition: ML MBA RGS.
- Investigation: JPS MGML JLMM RD DM.
- Methodology: ML MP.
- Project administration: ML.
- Resources: ML JLMM RD DM.
- Software: CB MBA RGS.
- Supervision: ML.
- Validation: MGML JLMM.
- Visualization: ML CB MBA.
- Writing – original draft: ML.
- Writing – review & editing: ML MP RGS.
References
- 1. Huang CJ, Zourdos MC, Jo E, Ormsbee MJ (2013) Influence of physical activity and nutrition on obesity-related immune function. ScientificWorldJournal 2013: 752071. pmid:24324381
- 2. Wu HJ, Wu E (2012) The role of gut microbiota in immune homeostasis and autoimmunity. Gut Microbes 3: 4–14. pmid:22356853
- 3. Flint HJ, Duncan SH, Scott KP, Louis P (2015) Links between diet, gut microbiota composition and gut metabolism. Proc Nutr Soc 74: 13–22. pmid:25268552
- 4. Tamboli CP, Neut C, Desreumaux P, Colombel JF (2004) Dysbiosis in inflammatory bowel disease. Gut 53: 1–4. pmid:14684564
- 5. Sobhani I, Tap J, Roudot-Thoraval F, Roperch JP, Letulle S, Langella P, et al. (2011) Microbial dysbiosis in colorectal cancer (CRC) patients. PloS one 6: e16393. pmid:21297998
- 6. Delzenne NM, Neyrinck AM, Bäckhed F, Cani PD (2011) Targeting gut microbiota in obesity: effects of prebiotics and probiotics. Nature Reviews Endocrinology 7: 639–646. pmid:21826100
- 7. Giongo A, Gano KA, Crabb DB, Mukherjee N, Novelo LL, Casella G, et al. (2011) Toward defining the autoimmune microbiome for type 1 diabetes. ISME J 5: 82–91. pmid:20613793
- 8. Larsen N, Vogensen FK, van den Berg FWJ, Nielsen DS, Andreasen AS, Pedersen BK, et al. (2010) Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PloS one 5: e9085. pmid:20140211
- 9. Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, Dugar B, et al. (2011) Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature 472: 57–63. pmid:21475195
- 10. Sekirov I, Russell SL, Antunes LC, Finlay BB (2010) Gut microbiota in health and disease. Physiological Reviews 90: 859–904. pmid:20664075
- 11. Marchesi JR, Adams DH, Fava F, Hermes GD, Hirschfield GM, Hold G, et al. (2016) The gut microbiota and host health: a new clinical frontier. Gut 65: 330–339. pmid:26338727
- 12. Brahe LK, Astrup A, Larsen LH (2016) Can We Prevent Obesity-Related Metabolic Diseases by Dietary Modulation of the Gut Microbiota? Advances in nutrition (Bethesda, Md) 7: 90–101.
- 13. Falony G, Joossens M, Vieira-Silva S, Wang J, Darzi Y, Faust K, et al.(2016) Population-level analysis of gut microbiome variation. Science (New York, NY) 352: 560–564.
- 14. Zhernakova A, Kurilshikov A, Bonder MJ, Tigchelaar EF, Schirmer M, Vatanen T, et al. (2016) Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science (New York, NY) 352: 565–569.
- 15. Backhed F, Manchester JK, Semenkovich CF, Gordon JI (2007) Mechanisms underlying the resistance to diet-induced obesity in germ-free mice. Proc Natl Acad Sci U S A 104: 979–984. pmid:17210919
- 16. Petriz BA, Castro AP, Almeida JA, Gomes CP, Fernandes GR, Kruger RH, et al. (2014) Exercise induction of gut microbiota modifications in obese, non-obese and hypertensive rats. BMC Genomics 15: 511. pmid:24952588
- 17. Evans CC, LePard KJ, Kwak JW, Stancukas MC, Laskowski S, Dougherty J, et al. (2014) Exercise prevents weight gain and alters the gut microbiota in a mouse model of high fat diet-induced obesity. PLoS One 9: e92193. pmid:24670791
- 18. Choi JJ, Eum SY, Rampersaud E, Daunert S, Abreu MT, Toborek M. (2013) Exercise attenuates PCB-induced changes in the mouse gut microbiome. Environmental health perspectives 121: 725–730. pmid:23632211
- 19. Queipo-Ortuño MI, Seoane LM, Murri M, Pardo M, Gomez-Zumaquero JM, Cardona F, et al. (2013) Gut microbiota composition in male Rat models under different nutritional status and physical activity and its association with serum leptin and ghrelin levels. PloS one 8: e65465. pmid:23724144
- 20. Kang SS, Jeraldo PR, Kurti A, Miller ME, Cook MD, Whitlock K, et al. (2014) Diet and exercise orthogonally alter the gut microbiome and reveal independent associations with anxiety and cognition. Mol Neurodegener 9: 36. pmid:25217888
- 21. Lambert JE, Myslicki JP, Bomhof MR, Belke DD, Shearer J, Reimer RA (2015) Exercise training modifies gut microbiota in normal and diabetic mice. Appl Physiol Nutr Metab: 1–4.
- 22. Mika A, Van Treuren W, González A, Herrera JJ, Knight R, Fleshner M (2015) Exercise is more effective at altering gut microbial composition and producing stable changes in lean mass in juvenile versus adult male F344 rats. PloS one 10: e0125889. pmid:26016739
- 23. Ott SJ, Schreiber S (2006) Reduced microbial diversity in inflammatory bowel diseases. Gut 55: 1207.
- 24. Chang JY, Antonopoulos DA, Kalra A, Tonelli A, Khalife WT, Schmidt TM, et al. (2008) Decreased diversity of the fecal Microbiome in recurrent Clostridium difficile-associated diarrhea. The Journal of infectious diseases 197: 435–438. pmid:18199029
- 25. Claesson MJ, Jeffery IB, Conde S, Power SE, O’Connor EM, Cusack S, et al. (2012) Gut microbiota composition correlates with diet and health in the elderly. Nature 488: 178–184. pmid:22797518
- 26.
World Health O (2010) Global recommendations on physical activity for health. Geneva: WHO Library Cataloguing-in-Publication Data.
- 27. Healy GN, Wijndaele K, Dunstan DW, Shaw JE, Salmon J, Zimmet PZ, et al. (2008) Objectively measured sedentary time, physical activity, and metabolic risk: the Australian Diabetes, Obesity and Lifestyle Study (AusDiab). Diabetes Care 31: 369–371. pmid:18000181
- 28. Freedson PS, Melanson E, Sirard J (1998) Calibration of the Computer Science and Applications, Inc. accelerometer. Medicine and science in sports and exercise 30: 777–781. pmid:9588623
- 29. Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, et al. (2013) Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res 41: e1. pmid:22933715
- 30. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. (2010) QIIME allows analysis of high-throughput community sequencing data. Nature methods 7: 335–336. pmid:20383131
- 31. Lozupone C, Knight R (2005) UniFrac: a new phylogenetic method for comparing microbial communities. Applied and Environmental Microbiology 71: 8228–8235. pmid:16332807
- 32. Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, et al. (2011) Enterotypes of the human gut microbiome. Nature 473: 174–180. pmid:21508958
- 33. Charrad M, Ghazzali N, Boiteau V, Niknafs A, Charrad MM (2014) Package ‘NbClust’. J Stat Soft 61: 1–36.
- 34. Caliński T, Harabasz J (1974) A dendrite method for cluster analysis. Communications in Statistics-theory and Methods 3: 1–27.
- 35. Brankatschk R, Bodenhausen N, Zeyer J, Burgmann H (2012) Simple absolute quantification method correcting for quantitative PCR efficiency variations for microbial community samples. Applied and Environmental Microbiology 78: 4481–4489. pmid:22492459
- 36. Qu W, Zhou Y, Zhang Y, Lu Y, Wang X, Zhao D, et al. (2012) MFEprimer-2.0: a fast thermodynamics-based program for checking PCR primer specificity. Nucleic Acids Res 40: W205–208. pmid:22689644
- 37. Delgado S, Ruas-Madiedo P, Suárez A, Mayo B (2006) Interindividual differences in microbial counts and biochemical-associated variables in the feces of healthy Spanish adults. Digestive diseases and sciences 51: 737–743. pmid:16614997
- 38. Martínez Álvarez JR, Villarino Marín A, Iglesias Rosado C, de Arpe Muñoz C, Gómez Candela C (2010) Recomendaciones de alimentación para la población española. Nutrición Clínica y Dietética Hospitalaria 30: 4–14.
- 39. Welly RJ, Liu TW, Zidon TM, Rowles JL 3rd, Park YM, Smith TN, et al. (2016) Comparison of Diet vs. Exercise on Metabolic Function & Gut Microbiota in Obese Rats. Medicine and science in sports and exercise.
- 40. Clarke SF, Murphy EF, O'Sullivan O, Lucey AJ, Humphreys M, Hogan A, et al. (2014) Exercise and associated dietary extremes impact on gut microbial diversity. Gut: gutjnl-2013-306541.
- 41. Denou E, Marcinko K, Surette MG, Steinberg GR, Schertzer JD (2016) High-intensity exercise training increases the diversity and metabolic capacity of the mouse distal gut microbiota during diet-induced obesity. American journal of physiologyEndocrinology and metabolism 310: E982–993.
- 42. Allen JM, Berg Miller ME, Pence BD, Whitlock K, Nehra V, Gaskins HR, et al. (2015) Voluntary and forced exercise differentially alters the gut microbiome in C57BL/6J mice. Journal of applied physiology (Bethesda, Md: 1985) 118: 1059–1066.
- 43. Knip M, Siljander H (2016) The role of the intestinal microbiota in type 1 diabetes mellitus. Nature Reviews Endocrinology.
- 44. Huttenhower C, Kostic AD, Xavier RJ (2014) Inflammatory bowel disease as a model for translating the microbiome. Immunity 40: 843–854. pmid:24950204
- 45. Ahn J, Sinha R, Pei Z, Dominianni C, Wu J, Shi J, et al. (2013) Human gut microbiome and risk for colorectal cancer. Journal of the National Cancer Institute 105: 1907–1911. pmid:24316595
- 46. Ley RE, Backhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI (2005) Obesity alters gut microbial ecology. Proc Natl Acad Sci U S A 102: 11070–11075. pmid:16033867
- 47. Wu GD, Chen J, Hoffmann C, Bittinger K, Chen YY, Keilbaugh SA, et al. (2011) Linking long-term dietary patterns with gut microbial enterotypes. Science 334: 105–108. pmid:21885731
- 48. Gorvitovskaia A, Holmes SP, Huse SM (2016) Interpreting Prevotella and Bacteroides as biomarkers of diet and lifestyle. Microbiome 4: 15. pmid:27068581
- 49. Knights D, Ward TL, McKinlay CE, Miller H, Gonzalez A, McDonald D, et al. (2014) Rethinking "enterotypes". Cell Host Microbe 16: 433–437. pmid:25299329
- 50. Endoh K, Kuriki K, Kasezawa N, Tohyama K, Goda T (2016) Impact of Interactions Between Self-Reported Psychological Stress and Habitual Exercise on the Dietary Intake of Japanese Men and Women: a Large-Scale Cross-Sectional Study. Asian Pac J Cancer Prev 17: 2007–2017. pmid:27221888
- 51. Baumann S, Toft U, Aadahl M, Jorgensen T, Pisinger C (2015) The long-term effect of screening and lifestyle counseling on changes in physical activity and diet: the Inter99 Study—a randomized controlled trial. Int J Behav Nutr Phys Act 12: 33. pmid:25886540
- 52. Viitasalo A, Eloranta AM, Lintu N, Vaisto J, Venalainen T, Kiiskinen S, et al. (2016) The effects of a 2-year individualized and family-based lifestyle intervention on physical activity, sedentary behavior and diet in children. Prev Med 87: 81–88. pmid:26915641
- 53. Zhong Y, Nyman M, Fåk F (2015) Modulation of gut microbiota in rats fed high‐fat diets by processing whole‐grain barley to barley malt. Molecular nutrition & food research 59: 2066–2076.
- 54. Chen J, Chia N, Kalari KR, Yao JZ, Novotna M, Soldan MM, et al. (2016) Multiple sclerosis patients have a distinct gut microbiota compared to healthy controls. Scientific reports 6: 28484. pmid:27346372
- 55. Kasai C, Sugimoto K, Moritani I, Tanaka J, Oya Y, Inoue H, et al. (2016) Comparison of human gut microbiota in control subjects and patients with colorectal carcinoma in adenoma: terminal restriction fragment length polymorphism and next-generation sequencing analyses. Oncology reports 35: 325–333. pmid:26549775
- 56. Zhang X, Shen D, Fang Z, Jie Z, Qiu X, Zhang C, et al. (2013) Human gut microbiota changes reveal the progression of glucose intolerance. PloS one 8: e71108. pmid:24013136
- 57. Zhang X, Zhang D, Jia H, Feng Q, Wang D, Liang D, et al. (2015) The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nature medicine.
- 58. Morotomi M, Nagai F, Sakon H, Tanaka R (2009) Paraprevotella clara gen. nov., sp. nov. and Paraprevotella xylaniphila sp. nov., members of the family ‘Prevotellaceae’isolated from human faeces. International Journal of Systematic and Evolutionary Microbiology 59: 1895–1900. pmid:19567577
- 59. Shaw KA, Bertha M, Hofmekler T, Chopra P, Vatanen T, Srivatsa A, et al. (2016) Dysbiosis, inflammation, and response to treatment: a longitudinal study of pediatric subjects with newly diagnosed inflammatory bowel disease. Genome Medicine 8: 1.
- 60. Kampmann C, Dicksved J, Engstrand L, Rautelin H (2016) Composition of human faecal microbiota in resistance to Campylobacter infection. Clinical Microbiology and Infection 22: 61. e61–61. e68.
- 61. Nylund L, Nermes M, Isolauri E, Salminen S, Vos WM, Satokari R (2015) Severity of atopic disease inversely correlates with intestinal microbiota diversity and butyrate‐producing bacteria. Allergy 70: 241–244. pmid:25413686
- 62. Pryde SE, Duncan SH, Hold GL, Stewart CS, Flint HJ (2002) The microbiology of butyrate formation in the human colon. FEMS microbiology letters 217: 133–139. pmid:12480096
- 63. Leonel AJ, Alvarez-Leite JI (2012) Butyrate: implications for intestinal function. Current opinion in clinical nutrition and metabolic care 15: 474–479. pmid:22797568
- 64. Derrien M, Belzer C, de Vos WM (2016) Akkermansia muciniphila and its role in regulating host functions. Microbial pathogenesis.
- 65. Org E, Parks BW, Joo JW, Emert B, Schwartzman W, Kang EY, et al. (2015) Genetic and environmental control of host-gut microbiota interactions. Genome research 25: 1558–1569. pmid:26260972
- 66. Canani RB, Costanzo MD, Leone L, Pedata M, Meli R, Calignano A, (2011) Potential beneficial effects of butyrate in intestinal and extraintestinal diseases. World J Gastroenterol 17: 1519–1528. pmid:21472114
- 67. Machiels K, Joossens M, Sabino J, De Preter V, Arijs I, Eeckhaut V, et al. (2014) A decrease of the butyrate-producing species Roseburia hominis and Faecalibacterium prausnitzii defines dysbiosis in patients with ulcerative colitis. Gut 63: 1275–1283. pmid:24021287
- 68. Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermudez-Humaran LG, Gratadoux JJ, et al. (2008) Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proceedings of the National Academy of Sciences of the United States of America 105: 16731–16736. pmid:18936492
- 69. Hamilton MT, Hamilton DG, Zderic TW (2007) Role of low energy expenditure and sitting in obesity, metabolic syndrome, type 2 diabetes, and cardiovascular disease. Diabetes 56: 2655–2667. pmid:17827399
- 70. Healy GN, Dunstan DW, Salmon J, Cerin E, Shaw JE, Zimmet PZ, et al. (2008) Breaks in sedentary time: beneficial associations with metabolic risk. Diabetes care 31: 661–666. pmid:18252901
- 71. Gonzales T, Robert-Baudouy J (1996) Bacterial aminopeptidases: properties and functions. FEMS Microbiol Rev 18: 319–344. pmid:8703509
- 72. Bernalier-Donadille A (2010) [Fermentative metabolism by the human gut microbiota]. Gastroenterol Clin Biol 34 Suppl 1: S16–22.
- 73. Christensen JE, Dudley EG, Pederson JA, Steele JL (1999) Peptidases and amino acid catabolism in lactic acid bacteria. Antonie Van Leeuwenhoek 76: 217–246. pmid:10532381