Next Article in Journal
Comparative Lipidomic Study of Human Milk from Different Lactation Stages and Milk Formulas
Next Article in Special Issue
Mediterranean Diet Effects on Type 2 Diabetes Prevention, Disease Progression, and Related Mechanisms. A Review
Previous Article in Journal
Impact of Treatment with RUTF on Plasma Lipid Profiles of Severely Malnourished Pakistani Children
Previous Article in Special Issue
Adherence to the 2018 World Cancer Research Fund/American Institute for Cancer Research Recommendations and Breast Cancer in the SUN Project
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mediterranean Diet, Physical Activity and Gut Microbiome Composition: A Cross-Sectional Study among Healthy Young Italian Adults

by
Francesca Gallè
1,†,
Federica Valeriani
2,†,
Maria Sofia Cattaruzza
3,
Gianluca Gianfranceschi
2,
Renato Liguori
4,
Martina Antinozzi
3,
Beatriz Mederer
5,
Giorgio Liguori
1 and
Vincenzo Romano Spica
2,*
1
Department of Movement Sciences and Wellbeing, University of Naples “Parthenope”, 80133 Naples, Italy
2
Department of Movement, Human and Health Sciences, University of Rome “Foro Italico”, 00135 Rome, Italy
3
Department of Public Health and Infectious Diseases, Sapienza University, 00185 Rome, Italy
4
Department of Sciences and Technologies, University of Naples “Parthenope”, 80143 Naples, Italy
5
Department of Didactics of Language and Literatura, University of Granada, 18150 Granada, Spain
*
Author to whom correspondence should be addressed.
These Authors contributed equally to this work.
Submission received: 29 June 2020 / Revised: 16 July 2020 / Accepted: 19 July 2020 / Published: 21 July 2020
(This article belongs to the Special Issue Diet Quality and Human Health)

Abstract

:
Background. This cross-sectional study aimed to explore the microbial composition of the gut and its possible association with the Mediterranean diet (MD) after adjusting for demographic and anthropometric characteristics in a sample of healthy young Italian adults. Methods. Gut microbiota, demographic information, and data on adherence to MD and physical activity (PA) habits were collected in a sample of 140 university students (48.6% males, mean age 22.5 ± 2.9) with a mean body mass index (BMI) of 22.4 ± 2.8 kg/m2 (15.2–33.8) and a mean PA level of 3006.2 ± 2973.6 metabolic equivalent (MET)-minutes/week (148–21,090). Results. A high prevalence of Firmicutes and Bacteroidetes was found in all the fecal samples. Significant dissimilarities in the microbiota composition were found on the basis of MD adherence and PA levels (p = 0.001). At the genus level, Streptococcus and Dorea were highly abundant in overweight/obese individuals, Ruminococcus and Oscillospira in participants with lower adherence to MD, and Lachnobacterium in subjects with low levels of PA (p = 0.001). A significantly higher abundance of Paraprevotella was shown by individuals with lower BMI, lower MD adherence, and lower PA levels (p = 0.001). Conclusions. This study contributes to the characterization of the gut microbiome of healthy humans. The findings suggest the role of diet and PA in determining gut microbiota variability.

1. Background

Recent evidence has shown that gut microbiota, depending on its composition, may influence human health and diseases [1,2] and that it can be influenced by host and environmental factors [3,4,5]. In healthy adults, two bacterial phyla have been recognized as predominant: the Gram-positive Firmicutes and the Gram-negative Bacteroidetes [6]. However, the Firmicutes/Bacteroidetes (F/B) ratio has been observed to change with the nutritional status, age, and gender [3,5,7,8,9]. In particular, it was shown to be associated with the body mass index (BMI), since it tends to be higher in obese people and to decrease with weight loss [7,10,11].
Diet and physical activity (PA) are likely to play a role—in particular, diet is a fundamental factor in determining which nutrients gut microbes can use for their biological processes, and the metabolic products of those processes may have important impacts on health and disease [12]. In particular, the Mediterranean diet (MD), due to its inclusion of non-refined cereals, fruits and vegetables, olive oil, and red wine, may favor the presence of some bacteria in the gut, such as Bifidobacterium and Lactobacobacillus, which are typically health-related [13]. MD seems to be the best pattern to beneficially modulate gut microbiota biodiversity and, consequently, prevent several disorders, including gastrointestinal or even neurological disease, through the gut-brain axis [14,15,16]. The relationship between PA and the gut microbiota is still debated, but evidence shows that exercise may foster specific gut microbial profiles, especially among lean individuals [17,18].
A deeper understanding of the characteristics associated with gut microbiota composition in healthy people is necessary to improve knowledge about its role in physiological and pathological processes. Thus, the present study aimed to explore the bacterial composition of fecal samples and its possible association with dietary and PA habits in a sample of healthy young Italian adults. In particular, the BMI and adherence to MD patterns were used as indicators of diet, while weekly activities were used to assess the habitual levels of PA.

2. Materials and Methods

2.1. Setting and Participants

The study protocol was approved by the institutional review board of the University of Rome “Foro Italico” (number CAR 10/2018) and tested in a previous pilot study [19]. Students attending the University of Naples “Parthenope” and University of Rome “La Sapienza” were invited to participate to the study at the end of their lessons. During the invitation, a researcher explained the aim and the procedures of the study and guaranteed anonymity in the collection and treatment of personal information. Those who agreed to adhere were asked to sign a written informed consent in accordance with the standards of the Helsinki Declaration.

2.2. Questionnaires

Participants were asked to report through a questionnaire their age, gender, weight, height, and particular diet regimens (i.e., vegetarian or vegan), which did not imply exclusion. Intercurrent chronic diseases, pregnancy, food intolerances, concomitant infections (and the associated use of antibiotics and probiotics), or gastrointestinal surgical procedures in the 3 previous months were also investigated and represented exclusion criteria. The 9-items questionnaire developed by Martínez-González et al. and the short version of the International Physical Activity Questionnaire (IPAQ) were used to assess the dietary habits and the habitual PA levels of the students, respectively [20,21]. The first is a short screener employed to assess the adherence to cardioprotective MD, derived by the 136-item questionnaire validated in a sample of Spanish adults [22]. It allows to attribute a value of 1 to each reported dietary habit corresponding to the MD pattern, resulting in a maximum total score of 9, which accounts for the highest adherence level. The IPAQ total score is expressed in metabolic equivalent (MET)-minutes/week.

2.3. Analysis of Fecal Samples

A fecal swab (Copan Italia S.P.A., Brescia, Italy) was given to each participant, together with the instructions for the collection of a stool sample; participants were asked to deliver the swab on a planned day within two hours since collection. Samples were then stored at 4–8 °C in a refrigerated container and were taken within 24 h to the laboratory of the University of Rome “Foro Italico”, where they were processed with a previously validated protocol for DNA extraction from fecal traces [23]. Feces were weighed prior to the extraction, and DNA was purified and normalized; the libraries for next-generation sequencing (NGS) were prepared according to the 16S Metagenomic Sequencing Library Preparation Guide (part# 15044223 rev A; Illumina, San Diego, CA, USA). The PCR amplicons were obtained using primers containing overhang adapters: Ba27F 5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGAGAGTTTGATCCTGGCTCAG-3′ and Ba338R 5′-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGTGCTGCCTCCCGTAGGAGT-3′ [24,25]. Tagged PCR products were generated using primer pairs with unique barcodes through two-step PCR. In this strategy, target primers containing overhang adapters were used in the first PCR reaction to amplify the target gene, and that product was then used in the second PCR using primers-containing barcodes. Each amplification reaction had a total volume of 25 μL, containing 12.5 μL of KAPA HiFi Hot Start Ready Mix (Roche, Pleasanton, CA, USA), 5 μL of each primer (1 μM), and 2 μL of template DNA. Reactions were carried out on a Techne® TC-PLUS thermocycler (VWR International, LLC, Radnor, PA, USA). Following amplification, 5 μL of PCR product from each reaction was used for agarose gel (1%) electrophoresis to confirm amplification. The final concentration of cleaned DNA amplicon was determined using the Qubit PicoGreen dsDNA BR assay kit (Invitrogen, Grand Island, NY, USA) and validated on a Bioanalyzer DNA 1000 chip (Agilent, Santa Clara, CA, USA). Libraries were prepared using the MiSeq Reagent Kit Preparation Guide (Illumina, San Diego, CA, USA). Raw sequence data was processed using an in-house pipeline that was built on the Galaxy platform and incorporated various software tools to evaluate the quality of the raw sequence data (FASTA/Q Information tools, Mothur). All datasets were rigorously screened to remove low-quality reads (short reads >200 nt, zero-ambiguous sequences). Demultiplexing was performed to remove PhiX sequences and sort sequences; moreover, to minimize sequencing errors and ensure sequence quality, the reads were trimmed based on the sequence quality score using Btrim (an average quality score of 30 from the ends, and remove reads that are less any 200 bp after end-trimming) [26]. OTUs (operational taxonomic units) were clustered at a 97% similarity level, and final OTUs were generated based on the clustering results, and taxonomic annotations of individual OTUs were based on representative sequences using RDP’s 16S Classifier 2.5. Observed OTUs were defined as observed species. A level of 97% sequence identity is often chosen as representative of a species and 95% for a genus. The sequence reads were analyzed, also, in the cloud environment BaseSpace through the 16S Metagenomics app (version 1.0.1; Illumina®): the taxonomic database used was the Illumina-curated version (May 2013 release of the Greengenes Consortium Database) [27]. The raw sequencing data were submitted to NCBI Sequence Read Archive with the project accession number of PRJNA630035.

2.4. Statistical Analysis

A descriptive analysis was carried out on the collected information. Weight and height values were used to calculate the BMI. Age, BMI, adherence to MD, and levels of habitual PA were expressed as mean ± SD and ranges; numbers and percentages of participants for each gender and BMI category as defined by the World Health Organization standards [28] were also reported. Participants were also grouped by MD scores (lower or equal/higher than the median value obtained in the whole sample) and by PA levels (low, moderate, or high) [29].
Relative abundances of community members were determined with rarefied data and summarized at each taxonomic level. The proportion of the gut microbiome at each taxonomic rank, such as phylum, order, class, family, and genus, was determined using the RDP classifier and the Greengenes Database. The standard pipeline for 16S rDNA amplicon analysis revealed that the bacteria fell into 29 phylotypes (similarity level = 97%). Briefly, the relative abundance of the bacterial community was determined on 29 phylotypes (% abundance = (number of sequences for phylum or genera/total sequences for sample) × 100). The subsequent mathematical analyses excluded the 24 rarest phyla—that is, those with a relative abundance below 0.1 % in all the samples. A total of 23,388,649 sequences were generated from 140 samples. The number of sequences for each sample ranged from 42,072 to 437,525, leading to the identification of 630 OTUs defined at 97% identity. Alpha and beta diversity were calculated using EstimateS software at a level of 97% sequence similarity. Regarding alpha diversity, the Shannon index and equitability index at the species level were computed [30,31]. Principal coordinates analysis (PCoA) was performed using the METAGENassist platform [32] in order to investigate the dissimilarity between the groups. To assess the sequencing depth, alpha rarefaction plots were done using the software Mothur (version 1.31.1, www.mothur.org) and R (version 3.1.3, www.R-project.org)) with packages “ggplot2” and “vegan” (R Core team 2013). The R package was employed for comparative analyses of taxonomic and functional microbiota composition. Heatmaps and clustering analyses were based on the Pearson’s correlation as a measure of the distance and Ward’s method using the METAGENassist platform [32].
Student’s t-test was used to compare the bacterial diversity and relative abundance of the main bacterial phyla and genus (accounting for the 90% of fecal bacterial component) of samples from participants with different BMI (included in the underweight/normal weight and overweight/obese groups) and levels of MD adherence. ANOVA with a Bonferroni post hoc test was employed to evaluate differences among participants grouped by habitual PA levels and to compare lifestyle subgroups obtained by merging MD adherence and PA. The statistical significance of factors potentially contributing to compositional differences among microbiota samples was tested also with the PRIMER software (version 7) by using the nonparametric permutation Analysis Of Similarity (ANOSIM) function and the default setting [33]. In order to identify possible determinants of the microbiome composition, two linear regression analyses were performed by considering age, gender, BMI, MD adherence, and PA level as independent variables and the Shannon index or F/B ratio as the outcomes. A stepwise procedure was used in order to adjust for possible confounders. A value of p < 0.05 was considered statistically significant. Data were analyzed with IBM SPSS version 25 for Windows (SPSS, Chicago, IL, USA).

3. Results

About 450 undergraduates attending classes were invited to take part in the study. Of a total of 244 (54.2%) students who initially agreed to participate, 140 (57.3%) completed the questionnaire, consigned the fecal swab on the planned date, and fulfilled the inclusion criteria (Figure 1).
Table 1 shows the main characteristics of the total sample. No particular diet regimens were reported by the participants.
Out of 28 bacterial phyla detected in the fecal samples, the highest abundances were registered for Firmicutes (61.6 ± 14.6) and Bacteroidetes (30.7 ± 13.3) (Figure 2).
Table 2 shows the differences in the Shannon index, Firmicutes, and Bacteroidetes relative abundances and F/B ratios between the BMI and MD score groups and among the PA groups, with corresponding p-values. No significant differences in variability nor in the phyla relative abundance were found.
In line with these results, also, the ANOSIM test yielded no significant dissimilarity for the BMI groups (R = −0.011, p = 0.5). In the regression analyses, only the male gender was significantly associated with the Shannon index (p = 0.02, odds ratio (OR) −0.08, confidence interval (CI) 95% −0.16–−0.01) and with the F/B ratio (p = 0.03, OR 0.50, CI 95% 0.08–0.93).
The global structure of the gut microbiota composition, as well as grouping patterns based on the MD score or PA level, are shown in Figure 3a and Figure 3b, respectively. As for the MD score, it is possible to observe slight differences also supported by the ANOSIM analysis (R value = 0.120, p = 0.001). The dissimilarity among the PA groups emerged from the principal component analysis, as shown in Figure 3b. The results of the ANOSIM indicated a slight dissimilarity among the three groups (R value = 0.350, p = 0.001).
Merging declared lifestyles, we obtained subgroups defined by the MD score and PA level. In the comparison of the abovementioned parameters among these lifestyle subgroups, no significant differences were detected (Table 3). The Bonferroni post hoc test also revealed no significant differences among the subgroups.
Plotting the correlation between the genera and adherence to the MD, an influence of the diet on the structure of the microbiota could be detected, although not strong in terms of absolute values (Figure 4). Bacteroidetes and Prevotella genera are the most represented in all fecal samples (19.8 ± 11.1 and 7 ± 13, respectively). However, no significant difference was detected between the MD score groups.
Table 4 shows the significant differences detected in the comparison of the genera relative abundances between the BMI and MD score categories and among the PA groups. The whole set of comparisons are reported in Supplementary Tables S1–S3. The Streptococcus and Dorea genera were highly represented in overweight/obese individuals, while Ruminococcus and Oscillospira were abundant among participants with lower adherence to MD. The Lachnobacterium genus was more abundant in subjects with low levels of PA. Furthermore, it is possible to observe some overlapping in these profiles. In particular, significantly higher amounts of Megasphaera, Dialister, and Lachnobacterium were found in overweight/obese and low-active subjects; a significantly higher abundance of Paraprevotella was shown also by individuals with lower BMI, lower MD adherence, and lower PA levels. The post-hoc Bonferroni test showed significant differences in genera abundance between the low PA group and the others.

4. Discussion

This study provides data on the composition of the gut microbiota in a Mediterranean population of young adults with different habits regarding MD and PA.
The results confirm the higher prevalence of Firmicutes and Bacteroidetes as prevalent phyla in the gut microbiota of healthy humans and further support the a role for diet and PA in determining the variability of the gut microbiota. Indeed, significant dissimilarities in fecal microbiota compositions were found between the MD groups and among the PA groups.
No significant differences in Firmicutes and Bacteroidetes abundances were observed on the basis of the BMI classification, as reported by Peters et al. [34]. As for the MD adherence, some studies performed on adults from Italy as well as other countries of the Mediterranean basin have shown the link between the MD and gut microbiota [13,35,36]. In line with the study by Mitsou et al. [13], we did not find significant differences related to diet in the relative amounts of the main phyla. However, our results were not in agreement with those from De Filippis et al. and Garcia-Mantrana et al., who reported an association between higher F/B ratios and a lower adherence to the MD or a high intake of animal-associated foods [35,36].
Nevertheless, at the genus level, several important differences were highlighted in our sample. In particular, as previously reported, the Streptococcus genus was highly abundant among overweight/obese participants, so as the Dorea and Dialister, which were also mainly represented among low-active individuals [37,38]. Interestingly, Dorea and Dialister were found to be related with insulin secretion and fasting blood glucose, suggesting their role in type 2 diabetes development in overweight-obese individuals [38]. Unlike what others have observed, Oscillospira was higher in subjects with a low adherence to MD patterns [37]. Lactic acid bacteria, such as the Lactobacillus and Lactococcus genera, were more represented in individuals who reported a higher adherence to the MD. This may be related to the dietary polyphenols, such as those present in olive oil, that represent substrates for Lactobacteria, mainly Lactobacillus strains, that are known to contributee in the healthy maintenance of the gut microbiota balance by modulating the oxidative status of the intestinal barrier, the inflammatory processes, and the immune response [39]. Lachnospira was more represented among individuals with higher adherence to the MD, while Paraprevotella was more represented in those with lower MD scores. The literature reports that vegetable-based diets tend to correlate with the amount of short-chain fatty acids (SCFA), which are produced by gut bacteria through dietary fiber fermentation and have fundamental effects on the immunity, metabolism, and inflammation processes of the host [38]. Lachnospira among Firmicutes, so as Prevotella among Bacteroidetes, seem to be involved in the fermentation of fibers, leading to a higher SCFA production. Even if the abundance of Firmicutes is generally associated with a higher intake of animal-associated foods, several members of this phylum were associated with SCFA production and were reported to be more abundant in individuals with higher adherence to the MD patterns and higher consumption of vegetables [17,39,40,41].
Our results did not show any correlation between Prevotella abundance and adherence to the MD pattern. Otherwise, the abundance of Paraprevotella genus was found to be significantly higher in the group with lower MD adherence, so as in less-active people, in accordance with previous studies [41,42]. Differences in the microbiota composition detected among the PA groups seem to confirm the role of exercise in determining the gut colonization by bacteria capable of SCFA production [17]. Furthermore, in our study, the relative abundance of Ruminoccoccus was found to be higher in subjects with a lower level of adherence to MD, and this is in line with the work of Martinez-Medina et al., showing an association between this genus and the consumption of a diet rich in saturated fat and sugars [43]. The Megasphaera genus was found to be more abundant in low-active subjects, in accordance to the findings from Liang et al., who reported a higher abundance of this genus in lower-level athletes in comparison with higher-level ones [40]. The higher abundance of Lachnobacterium in low-active participants is a novel observation that needs to be explored in depth. Interestingly, indeed, we observed several significant differences in genera abundance in the inactive group respect to the others, underlining a possible influence of physical activity in the presence of Lachnobacterium within the microbiota composition.
The analysis of biodiversity indexes and F/B ratio did not show any statistical significance among the different lifestyle subgroups. This could be due to the small sample size. Larger studies should be performed to verify these findings. In addition to the size of the sample, this study presents other limitations. First of all, the BMI, adherence to the MD, and the PA levels were calculated on the basis of self-reported information, and this could have led to inaccuracy or bias. Furthermore, we used a short and easy questionnaire to assess the MD adherence to facilitate the compliance of participants, but this did not allow us to perform a detailed analysis related to the quality and the quantity of nutrient intakes, e.g., fibers or animal proteins. However, our study offers an additional description of the microbiota composition in young adults from a Mediterranean population in association with their lifestyles and further studies are needed to characterize in depth the gut microbiota biodiversity in relation to several individual, social, and behavioral factors.

5. Conclusions

This study is a contribution to the characterization of the gut microbiome of healthy young Italian adults. It provides a picture of the relative abundance of the main bacterial phyla in the human gut and suggests the influence of diet and PA on the gut microbiota composition. From the perspectives of preventive medicine and public health, the results of this study underline the importance of promoting healthy behaviors in the population.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6643/12/7/2164/s1, Table S1: Mean abundance values ± SD of bacterial genera from fecal samples categorized per BMI category of participants. * p < 0.05, Table S2: Mean abundance values ± SD of bacterial genera from fecal samples categorized per MD score of participants. * p < 0.05, Table S3: Mean abundance values ± SD of bacterial genera from fecal samples categorized per PA level of participants. * p < 0.05.

Author Contributions

Project administration and funding acquisition: F.G. and G.L. Conceptualization, data curation, formal analysis, methodology, investigation, and writing—original draft preparation: F.G., F.V., and M.S.C. Investigation: G.G., R.L., B.M. and M.A. Supervision and writing—review and editing: F.G., V.R.S., G.L., and M.S.C. Writing—review and editing: B.M. All authors have read and agreed to the published version of the manuscript.

Funding

University of Naples Parthenope funding for individual research, years 2016–17. Recipients: F.G. and G.L.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Grenham, S.; Clarke, G.; Cryan, J.F.; Dinan, T.G. Braingut-microbe communication in health and disease. Front. Physiol. 2011, 2, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Belkaid, Y.; Hand, T.W. Role of the microbiota in immunity and inflammation. Cell 2014, 157, 121–141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Mariat, D.; Firmesse, O.; Levenez, F.; Guimarăes, V.; Sokol, H.; Doré, J.; Corthier, G.; Furet, J.P. The Firmicutes/Bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. 2009, 9, 123. [Google Scholar] [CrossRef] [PubMed]
  4. Walker, A.W.; Ince, J.; Duncan, S.H.; Webster, L.M.; Holtrop, G.; Ze, X.; Brown, D.; Stares, M.D.; Scott, P.; Bergerat, A.; et al. Dominant and diet-responsive groups of bacteria within the human colonic microbiota. ISME J. 2011, 5, 220–230. [Google Scholar] [CrossRef]
  5. Mahnic, A.; Rupnik, M. Different host factors are associated with patterns in bacterial and fungal gut microbiota in Slovenian healthy cohort. PLoS ONE 2018, 13, e0209209. [Google Scholar] [CrossRef] [Green Version]
  6. Arumugam, M.; Raes, J.; Pelletier, E.; Le Paslier, D.; Yamada, T.; Mende, D.R.; Fernandes, G.R.; Tap, J.; Bruls, T.; Batto, J.; et al. Enterotypes of the human gut microbiome. Nature 2011, 473, 174–180. [Google Scholar] [CrossRef]
  7. Koliada, A.; Syzenko, G.; Moseiko, V.; Budovska, L.; Puchkov, K.; Perederiy, V.; Gavalko, Y.; Dorofeyev, A.; Romanenko, M.; Tkach, S.; et al. Association between body mass index and Firmicutes/Bacteroidetes ratio in an adult Ukrainian population. BMC Microbiol. 2017, 17, 120. [Google Scholar] [CrossRef] [Green Version]
  8. Zhernakova, A.; Kurilshikov, A.; Bonder, M.J.; Tigchelaar, E.F.; Schirmer, M.; Vatanen, T.; Mujagic, Z.; Vila, A.V.; Falony, G.; Vieira-Silva, S.; et al. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science 2016, 352, 565–569. [Google Scholar] [CrossRef] [Green Version]
  9. Odamaki, T.; Kato, K.; Sugahara, H.; Hashikura, N.; Takahashi, S.; Xiao, J.; Abe, F.; Osawa, R. Age-related changes in gut microbiota composition from newborn to centenarian: A cross-sectional study. BMC Microbiol. 2016, 16, 90. [Google Scholar] [CrossRef] [Green Version]
  10. Sweeney, T.E.; Morton, J.M. The human gut microbiome: A review of the effect of obesity and surgically induced weight loss. JAMA Surg 2013, 148, 563–569. [Google Scholar] [CrossRef] [Green Version]
  11. Mathur, R.; Barlow, G.M. Obesity and the microbiome. Expert Rev. Gastroenterol. Hepatol. 2015, 9, 1087–1099. [Google Scholar] [CrossRef] [PubMed]
  12. Gentile, C.L.; Weir, T.L. The gut microbiota at the intersection of diet and human health. Science 2018, 362, 776–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mitsou, E.K.; Kakali, A.; Antonopoulou, S.; Mountzouris, K.C.; Yannakoulia, M.; Panagiotakos, D.B.; Kyriacou, A. Adherence to the Mediterranean diet is associated with the gut microbiota pattern and gastrointestinal characteristics in an adult population. Br. J. Nutr. 2017, 117, 1645–1655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rinninella, E.; Cintoni, M.; Raoul, P.; Lopetuso, L.R.; Scaldaferri, F.; Pulcini, G.; Miggiano, G.A.D.; Gasbarrini, A.; Mele, M.C. Food Components and Dietary Habits: Keys for a Healthy Gut Microbiota Composition. Nutrients 2019, 11, 2393. [Google Scholar] [CrossRef] [Green Version]
  15. Di Liddo, R.; Piccione, M.; Schrenk, S.; Dal Magro, C.; Cosma, C.; Padoan, A.; Contran, N.; Scapellato, M.L.; Pagetta, A.; Spica, V.R.; et al. S100B as a new fecal biomarker of inflammatory bowel diseases. Eur. Rev. Med. Pharm. Sci 2020, 24, 323–332. [Google Scholar]
  16. Mayer, E.A.; Tillisch, K.; Gupta, A. Gut/brain axis and the microbiota. J. Clin. Investig. 2015, 125, 926–938. [Google Scholar] [CrossRef] [PubMed]
  17. Allen, J.M.; Mailing, L.J.; Niemiro, G.M.; Moore, R.; Cook, M.D.; White, B.A.; Holscher, H.D.; Woods, J.A. Exercise alters gut microbiota composition and function in lean and obese humans. Med. Sci. Sports Exerc. 2018, 50, 747–757. [Google Scholar] [CrossRef]
  18. Gallè, F.; Valeriani, F.; Cattaruzza, M.S.; Ubaldi, F.; Romano Spica, V.; Liguori, G.; WDPP, Working Group on Doping Prevention Project; GSMS-SItI. Exploring the association between physical activity and gut microbiota composition: A review of current evidence. Ann. Ig. 2019, 31, 582–589. [Google Scholar]
  19. Valeriani, F.; Gallè, F.; Cattaruzza, M.S.; Antinozzi, M.; Gianfranceschi, G.; Postiglione, N.; Romano Spica, V.; Liguori, G. Are nutrition and physical activity associated with gut microbiota? A pilot study on a sample of healthy young adults. Ann. Ig. 2020, 32. [epub ahead of print]. [Google Scholar]
  20. Martínez-González, M.A.; Fernández-Jarne, E.; Serrano-Martínez, M.; Wright, M.; Gomez-Gracia, E. Development of a short dietary intake questionnaire for the quantitative estimation of adherence to a cardioprotective Mediterranean diet. Eur. J. Clin. Nutr. 2004, 58, 1550–1552. [Google Scholar] [CrossRef] [Green Version]
  21. Craig, C.L.; Marshall, A.L.; Sjöström, M.; Bauman, A.E.; Booth, M.L.; Ainsworth, B.E.; Pratt, M.; Ekelund, U.; Yngve, A.; Sallis, J.F.; et al. International physical activity questionnaire: 12-country reliability and validity. Med. Sci. Sports Exerc. 2003, 35, 1381–1395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Martin-Moreno, J.M.; Boyle, P.; Gorgojo, L.; Maisonneuve, P.; Fernandez-Rodriguez, J.C.; Salvini, S.; Willett, W.C. Development and validation of a food frequency questionnaire in Spain. Int. J. Epidemiol. 1993, 22, 512–519. [Google Scholar] [CrossRef] [PubMed]
  23. Valeriani, F.; Agodi, A.; Casini, B.; Cristina, M.L.; D’Errico, M.M.; Gianfranceschi, G.; Liguori, G.; Liguori, R.; Mucci, N.; Mura, I.; et al. Potential testing of reprocessing procedures by real-time polymerase chain reaction: A multicenter study of colonoscopy devices. Am. J. Infect. Control. 2018, 46, 159–164. [Google Scholar] [CrossRef] [Green Version]
  24. Valeriani, F.; Crognale, S.; Protano, C.; Gianfranceschi, G.; Orsini, M.; Vitali, M.; Romano Spica, V. Metagenomic analysis of bacterial community in a travertine depositing hot spring. New Microbiol. 2018, 41, 126–135. [Google Scholar] [PubMed]
  25. Valeriani, F.; Protano, C.; Gianfranceschi, G.; Leoni, E.; Galasso, V.; Mucci, N.; Vitali, M.; Romano Spica, V. Microflora Thermarum Atlas project: Biodiversity in thermal spring waters and natural SPA pools. Water Sci. Technol. Water Supply 2018, 18, 1472–1483. [Google Scholar] [CrossRef]
  26. Kong, Y. Btrim: A fast, lightweight adapter and quality trimming program for next-generation sequencing technologies. Genomics 2011, 98, 152–153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Wang, S.; Hou, W.; Dong, H.; Jiang, H.; Huang, L.; Wu, G.; Zhang, C.; Song, Z.; Zhang, Y.; Ren, H.; et al. Control of Temperature on Microbial Community Structure in Hot Springs of the Tibetan Plateau. PLoS ONE 2013, 8, 62901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. World Health Organization. Body Mass Index—BMI. Available online: http://www.euro.who.int/en/health-topics/disease-prevention/nutrition/a-healthy-lifestyle/body-mass-index-bmi (accessed on 3 June 2020).
  29. IPAQ Group. Guidelines for Data Processing and Analysis of the International Physical Activity Questionnaire (IPAQ)—Short and Long Forms. Available online: www.ipaq.ki.se (accessed on 3 June 2020).
  30. Colwell, R.K.; Chao, A.; Gotelli, N.J.; Lin, S.Y.; Mao, C.X.; Chazdon, R.L.; Longino, J.T. Models and estimators linking individual-based and sample-based rarefaction, extrapolation, and comparison of assemblages. J. Plant. Ecol. 2012, 5, 3–21. [Google Scholar] [CrossRef] [Green Version]
  31. Magurran, A.E. Measuring Biological Diversity; Wiley-Blackwell: Oxford, UK, 2016. [Google Scholar]
  32. Arndt, D.; Xia, J.; Liu, Y.; Zhou, Y.; Guo, A.C.; Cruz, J.A.; Sinelnikov, I.; Budwill, K.; Nesbø, C.L.; Wishart, D.S. METAGENassist: A comprehensive web server for comparative metagenomics. Nucleic Acids Res. 2012, 40, 88–95. [Google Scholar] [CrossRef] [Green Version]
  33. Clark, K.R. Gorley RN PRIMER v6: User Manual/Tutorial; Primer-e: Plymouth, UK, 2006. [Google Scholar]
  34. Peters, B.A.; Shapiro, J.A.; Church, T.R.; Miller, G.; Trinh-Shevrin, C.; Yuen, E.; Friedlander, C.; Hayes, R.B.; Ahn, J. A taxonomic signature of obesity in a large study of American adults. Sci. Rep. 2018, 8, 9749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. De Filippis, F.; Pellegrini, N.; Vannini, L.; Jeffery, I.B.; La Storia, A.; Laghi, L.; Serrazanetti, D.I.; Di Cagno, R.; Ferrocino, I.; Lazzi, C.; et al. High-level adherence to a Mediterranean diet beneficially impacts the gut microbiota and associated metabolome. Gut 2016, 65, 1812–1821. [Google Scholar] [CrossRef] [PubMed]
  36. Garcia-Mantrana, I.; Selma-Royo, M.; Alcantara, C.; Collado, M.C. Shifts on Gut Microbiota Associated to Mediterranean Diet Adherence and Specific Dietary Intakes on General Adult Population. Front. Microbiol. 2018, 9, 890. [Google Scholar]
  37. Naderpoor, N.; Mousa, A.; Gomez-Arango, L.F.; Barrett, H.L.; Dekker Nitert, M.; de Courten, B. Faecal Microbiota Are Related to Insulin Sensitivity and Secretion in Overweight or Obese Adults. J. Clin. Med. 2019, 8, 452. [Google Scholar] [CrossRef] [Green Version]
  38. Luisi, M.; Lucarini, L.; Biffi, B.; Rafanelli, E.; Pietramellara, G.; Durante, M.; Vidali, S.; Provensi, G.; Madiai, S.; Gheri, C.F.; et al. Effect of Mediterranean Diet Enriched in High Quality Extra Virgin Olive Oil on Oxidative Stress, Inflammation and Gut Microbiota in Obese and Normal Weight Adult Subjects. Front. Pharm. 2019, 10, 1366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Tannock, G.W.; Liu, Y. Guided dietary fibre intake as a means of directing short-chain fatty acid production by the gut microbiota. J. R. Soc. N. Z. 2020, 50, 434–455. [Google Scholar] [CrossRef]
  40. Liang, R.; Zhang, S.; Peng, X.; Yang, W.; Xu, Y.; Wu, P.; Chen, J.; Cai, Y.; Zhou, J. Characteristics of the gut microbiota in professional martial arts athletes: A comparison between different competition levels. PLoS ONE 2019, 14, e0226240. [Google Scholar] [CrossRef]
  41. Bressa, C.; Bailén-Andrino, M.; Pérez-Santiago, J.; González-Soltero, R.; Pérez, M.; Montalvo-Lominchar, M.G.; Maté-Muñoz, J.L.; Domínguez, R.; Moreno, D.; Larrosa, M. Differences in gut microbiota profile between women with active lifestyle and sedentary women. PLoS ONE 2017, 12, e0171352. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Shikany, J.M.; Demmer, R.T.; Johnson, A.J.; Fino, N.F.; Meyer, K.; Ensrud, K.E.; Lane, N.E.; Orwoll, E.S.; Kado, D.M.; Zmuda, J.M.; et al. Association of dietary patterns with the gut microbiota in older, community-dwelling men. Am. J. Clin. Nutr. 2019, 110, 1003–1014. [Google Scholar] [CrossRef] [PubMed]
  43. Martinez-Medina, M.; Denizot, J.; Dreux, N.; Robin, F.; Billard, E.; Bonnet, R.; Darfeuille-Michaud, A.; Barnich, N. Western diet induces dysbiosis with increased E coli in CEABAC10 mice, alters host barrier function favouring AIEC colonisation. Gut 2014, 63, 116–124. [Google Scholar] [CrossRef] [PubMed]
Figure 1. CONSORT diagram for the sample generation.
Figure 1. CONSORT diagram for the sample generation.
Nutrients 12 02164 g001
Figure 2. Box plot showing the relative abundance for the dominant phyla determined by the Global Alignment Sequence Taxonomy assignments. The boundary of the box closest to zero indicates the 25th percentile, the line within the box represents the median, and the boundary of the box farthest from zero indicates the 75th percentile. Whiskers above and below the box indicate the 10th and 90th percentiles. The arithmetic mean is indicated by x. Phyla less than 1% of the sequences were not reported.
Figure 2. Box plot showing the relative abundance for the dominant phyla determined by the Global Alignment Sequence Taxonomy assignments. The boundary of the box closest to zero indicates the 25th percentile, the line within the box represents the median, and the boundary of the box farthest from zero indicates the 75th percentile. Whiskers above and below the box indicate the 10th and 90th percentiles. The arithmetic mean is indicated by x. Phyla less than 1% of the sequences were not reported.
Nutrients 12 02164 g002
Figure 3. Principal coordinates analysis (PCoA) of the weighted Bray–Curtis for 16S rDNA data related to (a) the Mediterranean diet (MD) score: 1: low adherence (score ranging from 0 to 5) and 2: higher adherence to MD (score > 5) and (b) habitual physical activity (PA) level: 1: low, 2: moderate, and 3: high. Percentages on the axes represent the proportion of the variation explained by the 2 first components of the PCoA.
Figure 3. Principal coordinates analysis (PCoA) of the weighted Bray–Curtis for 16S rDNA data related to (a) the Mediterranean diet (MD) score: 1: low adherence (score ranging from 0 to 5) and 2: higher adherence to MD (score > 5) and (b) habitual physical activity (PA) level: 1: low, 2: moderate, and 3: high. Percentages on the axes represent the proportion of the variation explained by the 2 first components of the PCoA.
Nutrients 12 02164 g003
Figure 4. Heatmap and clustering of individual gut microbiota samples for taxonomic composition (genus level) in relation to the MD adherence. Clustering was based on Pearson distances and Ward’s linkage method. Only the genera above 2% abundance in at least one sample are depicted.
Figure 4. Heatmap and clustering of individual gut microbiota samples for taxonomic composition (genus level) in relation to the MD adherence. Clustering was based on Pearson distances and Ward’s linkage method. Only the genera above 2% abundance in at least one sample are depicted.
Nutrients 12 02164 g004
Table 1. Main characteristics of the total samples.
Table 1. Main characteristics of the total samples.
Participants
n = 140
Age
mean value ± SD (range)22.5 ± 2.9 (18–36)
Gender
n (%)
males68 (48.6)
females72 (51.4)
BMI kg/m2
mean value ± SD (range)22.4 ± 2.8 (15.2–33.8)
BMI category
n (%)
under weight7 (5.0)
normal weight106 (75.7)
over weight24 (17.1)
obese3 (2.2)
MD adherence score
mean value ± SD (median; range)5.3 ± 1.6 (5; 2–9)
MD adherence level
n (%)
≤median73 (52.1)
>median67 (47.9)
Habitual PA
MET-minutes/week
mean value ± SD (range)3006.2 ± 2973.6 (148–21,090)
PA level
n (%)
low17 (12.1)
moderate57 (40.7)
high66 (47.1)
BMI, Body Mass Index, PA, physical activity, MD, Mediterranean diet, and MET, metabolic equivalent.
Table 2. Comparisons of the Shannon index, Firmicutes, Bacteroidetes, and the Firmicutes/Bacteroidetes (F/B) ratio of the BMI, MD score, and PA level groups, with corresponding p-values.
Table 2. Comparisons of the Shannon index, Firmicutes, Bacteroidetes, and the Firmicutes/Bacteroidetes (F/B) ratio of the BMI, MD score, and PA level groups, with corresponding p-values.
BMIpa
Underweight/
Normal Weight
Overweight/
Obese
Shannon Index2.5 ± 0.22.5 ± 0.20.77
Firmicutes58.9 ± 13.161.1 ± 8.70.47
Bacteroidetes33.4 ± 10.431.8 ± 8.90.54
F/B ratio2.1 ± 1.12.2 ± 1.30.56
MD scorepa
Variable0–56–9
Shannon Index2.6 ± 0.22.6 ± 0.20.96
Firmicutes58.9 ± 13.059.8 ± 11.90.68
Bacteroidetes33.6 ± 11.132.6 ± 9.10.58
F/B ratio2.1 ± 1.22.1 ± 1.00.94
PA levelpb
LowmoderateHigh
Shannon Index2.6 ± 0.12.6 ± 0.22.5 ± 0.20.46
Firmicutes58.3 ± 16.058.1 ± 12.060.6 ± 11.80.55
Bacteroidetes32.7 ± 14.834.6 ± 9.431.9 ± 9.40.38
F/B ratio2.3 ± 1.51.9 ± 0.92.2 ± 1.10.27
a Student’s t-test and b ANOVA.
Table 3. Comparisons of the Shannon Index, Firmicutes, Bacteroidetes, and F/B ratio of the subgroups categorized on the basis of MD adherence and PA levels, with corresponding p-values.
Table 3. Comparisons of the Shannon Index, Firmicutes, Bacteroidetes, and F/B ratio of the subgroups categorized on the basis of MD adherence and PA levels, with corresponding p-values.
Habitual PA LevelMD Adherence Scorep
0–5
n = 73
6–9
n = 67
Shannon Indexlow active2.6 ± 0.12.5 ± 0.2
moderately active2.6 ± 0.22.6 ± 0.20.47
highly active2.5 ± 0.22.6 ± 0.2
Firmicuteslow active57.1 ± 17.163.8 ± 9.1
moderately active59.0 ± 11.556.6 ± 12.90.51
highly active59.7 ± 12.661.1 ± 11.5
Bacteroideteslow active32.8 ± 16.131.8 ± 8.4
moderately active33.9 ± 9.135.6 ± 9.90.61
highly active33.4 ± 10.631.1 ± 8.6
F/B ratiolow active2.3 ± 1.62.1 ± 0.7
moderately active1.9 ± 1.01.8 ± 0.90.81
highly active2.1 ± 1.32.1 ± 1.0
Table 4. Relative abundance of significantly different genera in BMI, MD score, and PA groups.
Table 4. Relative abundance of significantly different genera in BMI, MD score, and PA groups.
GenusBMIpa
Underweight/
Normal Weight
Overweight/
Obese
Selenomonas0.23 ± 0.380.38 ± 1.300.02
Megasphaera0.12 ± 0.290.32 ± 0.870.001
Streptococcus1.66 ± 2.084.02 ± 6.690.001
Dorea1.16 ± 2.292.44 ± 4.250.001
Lachnobacterium0.27 ± 0.740.68 ± 2.430.007
Jannaschia0.19 ± 0.310.27 ± 0.500.02
Dialister0.35 ± 0.761.07 ± 2.340.001
Eubacterium0.71 ± 2.191.61 ± 3.420.01
Paraprevotella0.19 ± 0.300.08 ± 0.130.01
MD scorepa
Genus0–56–9
Lachnospira0.37 ± 1.200.75 ± 2.280.02
Oscillospira2.18 ± 1.901.73 ± 1.250.03
Lactobacillus0.29 ± 0.521.26 ± 4.520.002
Ruminococcus9.51 ± 6.697.63 ± 4.900.002
Lactococcus0.22 ± 0.530.52 ± 1.300.01
Veillonella0.09 ± 0.050.14 ± 0.170.001
Paraprevotella0.23 ± 0.380.13 ± 0.220.008
PA levelpb
Genuslowmoderatehigh
Megasphaera0.55 ± 1.08 10.12 ± 0.240.11 ± 0.370.002
Lachnobacterium1.30 ± 3.20 10.25 ± 0.410.17 ± 0.250.001
Dialister1.53 ± 2.69 10.40 ± 0.800.24 ± 0.580.001
Paraprevotella0.32 ± 0.49 20.21 ± 0.290.12 ± 0.250.037
a Student’s t-test and b ANOVA with Bonferroni post hoc test: 1 significantly different from moderate and high PA groups and 2 significantly different from the high PA group.

Share and Cite

MDPI and ACS Style

Gallè, F.; Valeriani, F.; Cattaruzza, M.S.; Gianfranceschi, G.; Liguori, R.; Antinozzi, M.; Mederer, B.; Liguori, G.; Romano Spica, V. Mediterranean Diet, Physical Activity and Gut Microbiome Composition: A Cross-Sectional Study among Healthy Young Italian Adults. Nutrients 2020, 12, 2164. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12072164

AMA Style

Gallè F, Valeriani F, Cattaruzza MS, Gianfranceschi G, Liguori R, Antinozzi M, Mederer B, Liguori G, Romano Spica V. Mediterranean Diet, Physical Activity and Gut Microbiome Composition: A Cross-Sectional Study among Healthy Young Italian Adults. Nutrients. 2020; 12(7):2164. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12072164

Chicago/Turabian Style

Gallè, Francesca, Federica Valeriani, Maria Sofia Cattaruzza, Gianluca Gianfranceschi, Renato Liguori, Martina Antinozzi, Beatriz Mederer, Giorgio Liguori, and Vincenzo Romano Spica. 2020. "Mediterranean Diet, Physical Activity and Gut Microbiome Composition: A Cross-Sectional Study among Healthy Young Italian Adults" Nutrients 12, no. 7: 2164. https://0-doi-org.brum.beds.ac.uk/10.3390/nu12072164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop