- 1Department of Sports and Exercise Science, Zhejiang University, Hangzhou, China
- 2Department of Neurology and Research Center of Neurology in Second Affiliated Hospital, Key Laboratory of Medical Neurobiology of Zhejiang Province, Interdisciplinary Institute of Neuroscience and Technology, Zhejiang University School of Medicine, Hangzhou, China
- 3Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering and Instrument Science, Zhejiang University School of Medicine, Hangzhou, China
- 4Department of Neurology of Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou, China
- 5Department of Medicine and Science in Sports and Exercise, Tohoku University Graduate School of Medicine, Sendai, Japan
Objective: This study aimed to investigate the disparity of gut microbiota among elite athletes and young adults with different physical activity independent of dietary status.
Methods: In Hangzhou, China, an age and sex matching study was conducted between April and May 2021. A total of 66 Chinese young adults were recruited in this study and divided into an elite athlete group, physically active group, and physically inactive group. Fecal samples were collected to assess gut microbiota composition. Dietary status was measured using a food-frequency questionnaire. Comparisons in gut microbiota and blood biomarkers among three groups were analyzed by using the analysis of covariance.
Results: The findings depicted a tendency to form clusters for beta diversity among three groups, while no significant difference was observed in both alpha and beta diversity. In the multiple analysis model, by adjusting dietary status, a significantly higher abundance of Clostridiaceae (p = 0.029) and Megamonas_rupellensis (p = 0.087) was observed in elite athletes compared to that in general young adults. Furthermore, inflammation-related bacteria such as Bilophila (p = 0.011) and Faecalicoccus (p = 0.050) were enriched in physically inactive young adults compared to two other groups. Pearson's correlation analysis showed a positive association between Bilophila and circulating white body cell count (r = 0.332, p = 0.006) and its subtypes including neutrophils (r = 0.273, p = 0.027), and lymphocytes (r = 0.327, p = 0.007). Megamonas_rupellensis has been shown associated positively with serum lymphocytes levels (r = 0.268, p = 0.03). Although no significant differences were observed, the elite athletes tended to have lower levels of blood biomarkers of immunity within a normal range, which may reflect a better immune function.
Conclusion: This matching study indicated that physically inactive young adults are more likely to have a lower immune function and a higher abundance of pro-inflammatory gut bacteria than elite athletes and physically active young adults. Dietary status should be considered as an important factor that may affect the association of physical activity with immune function and gut microbiota.
Introduction
Humans are superorganisms composed of human bodies and commensal microbiota (1). The human gut microbiota is approximately 100 trillion organisms, outnumbering the human cells by an estimated 10-fold (2–4). Increasing evidence suggests that there is a close relationship between gut microbiota composition, inflammation (5, 6), and immunity (7, 8). A systematic review showed that gut microbiota is thought to contribute to subacute systemic inflammatory and it may also be influenced by outcome, thereby reinforcing the disease symptoms (9). Some gut microbiota may be associated with pro-inflammatory effects, while others have anti-inflammatory properties. As potential anti-inflammatory microbiota, Lachnospiraceae and Ruminococcus are thought to prevent or ameliorate chronic hepatitis B (10). Conversely, pro-inflammatory gut microbiota, such as Bacteroides fragilis, were highly enriched in the gut of inflammatory arthritis (11). In addition, the Bacteroides fragilis toxin was also reported to be pro-inflammatory and carcinogenic (12). Thus, gut microbiota could be considered an important indicator of infectious and metabolic diseases.
It is well-known that people with high physical activity levels tend to have better immune function and lower levels of inflammation, which may help them improve health status and prevent diseases. Evidence showed that physically active adolescent girls had lower circulating IL-6 concentration compared to their sedentary counterparts (13). Moreover, a recent study reported that higher physical activity levels may be associated with lower salivary C-reactive protein in young adults (14). Studies regarding the relationship of exercise habits with inflammation-related gut microbiota are emerging. The new literature enables a better understanding of the underlying mechanisms of the association between physical activity and inflammation. Higher proportions of Akkermansiaceae family and Akkermansia genus in elite rugby athletes were detected than in overweight or obese young adults (15). Meanwhile, our previous randomized controlled trial showed that 8 weeks of combined aerobic and resistance exercise not only resulted in improved physical function but also increased abundance of anti-inflammatory bacteria and decreased pro-inflammatory bacteria in physically inactive older women (16). Even with these findings, the evidence exploring the distribution of inflammation-related gut microbiota in people with different physical activity levels is limited, and further studies are warranted.
Conversely, dietary status and nutrient intake are very important and non-negligible indicators when discussing physical activity levels and gut microbiota. Diet is considered to be a pivotal factor affecting the genes and composition of gut microbiota, of which the number, type, and balance of macronutrients (fats, proteins, and carbohydrates) all have a significant impact on gut microbiota (17). In addition, data indicated that the consumption of nuts, oily fish, fruit, vegetables, and cereals is linked to a higher abundance of anti-inflammatory bacteria (18). Meanwhile, dietary status was also shown to be associated with physical activity levels. Previous studies indicated that habitual Mediterranean dietary pattern is associated with more time spent on leisure-time physical activity among older adults (19). Momentary physical activity was concomitant with momentary consumption of both healthy and unhealthy dietary intake in college students (20). Thus, it is hypothesized that dietary status may play a contributory role in the association of physical activity with gut microbiota. However, most previous observational studies and exercise intervention trials have failed to confirm the relationship between physical activity and gut microbiota composition concerning independent of dietary status (15, 16, 21).
Therefore, this matching study aims to investigate the differences in the diversity and composition of gut microbiota among elite athletes and the general young adults with different physical activity levels, as well as explore whether the association between physical activity and gut microbiota are independent of dietary status.
Materials and Methods
Participants
An age and sex matching study was conducted between April 2021 and May 2021 in Hangzhou, China, and aimed to examine the difference of gut microbiota, anthropometry, life-behavior habits, and physical performance between athletes and physically active or inactive young adults. All elite athletes are engaged in track and field events (including 400 m, 1,500 m, high jump, long jump, discus, javelin, and pole vault), and had at least 8 years of professional training experience as well as participated in official sports events held by the Chinese Athletics Association or International Association of Athletics Federation.
Thirty young elite athletes were included in this study according to the following inclusion criteria: (1) 18–25 years old; (2) no history of cardiovascular diseases; (3) no history of congenital diseases; (4) no history of physical dysfunction. Participants from the control group included general young adults who met the physical activity recommendation levels according to the guidelines of the American College of Sports Medicine (22). A total of 64 general college students were invited to participate in this study. Among them, 34 students engaged in more than 150 min of moderate-intensity physical activity every week, and 30 students did not engage in any moderate or vigorous-intensity physical activity per day. Physical activity was assessed using the International Physical Activity Questionnaire (IPAQ), which has good reliability and validity for assessing physical activity in healthy adults (23). Thus, 94 elite athletes and college students were allowed to enter the matching process. These participants were then divided into three groups (Athlete group; High physical activity group, HPA; and Low physical activity group, LPA) according to equal proportion principle with age (±24 months) and sex. Finally, a total of 66 young adults were enrolled in the present study (Figure 1). The study protocol and all amendments were approved by the Ethics Committee of the Department of Psychology and Behavioral Sciences, Zhejiang University. All individuals provided written informed consent to participate in the study.
Dietary Evaluation
In the present study dietary information was collected by food-frequency questionnaire (FFQ25) which was demonstrated to have good reproducibility and validity (24). The FFQ25 primarily investigates the food intake of participants over the past week and contains 25 categories of foods. The frequency of intake of each food is divided into nine levels, from never eating to at least three times a day. The amount of food consumed is divided into five levels, from no more than 50 to at least 250 g a meal. Further, the frequency and amount categories were calculated as times per day by using a midpoint (e.g., “3–4 times per week” was used in the calculations as 3.5 times per week). Finally, all kinds of foods were divided into six groups: cereals, vegetables and fruits, meat, beans and nuts, fatty food, and alcohol.
Feces Samples Collection and DNA Extraction
The fecal samples were collected by the participants themselves, and the sample collection supplies and methods are uniformly provided and guided by an experimenter before the start of the sample collection. All samples were stored in dry ice at −80°C and transported to Huada Gene Detection Center for DNA extraction. The DNA of the intestinal microbial community was extracted using MagPure Stool DNA KF kit B (Magen, China). The Qubit® dsDNA BR Assay kit (Invitrogen, USA) was used to quantify the DNA with a Qubit Fluorometer. The regions V3–V4 of gut microbiota 16S rRNA genes were amplified using the degenerate PCR primers, 341F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). Illumina adapter, pad, and linker sequences were added to both forward and reverse primers. A 50 μL reaction containing a 30-ng template, fusion PCR primer, and PCR master mix was used for PCR enrichment. The following were the PCR cycling conditions; 94°C for 3 min; subsequently, 30 cycles of 94°C for 30 s; 56°C for 45 s; 72°C for 45 s, and final extension at 72°C for 10 min. AmpureXP beads were used to purify the PCR products, which were then eluted in the Elution buffer. The Agilent 2,100 bioanalyzer (Agilent, USA) was used to qualify the libraries. Validated libraries were sequenced on the IlluminaMiSeq platform (BGI, Shenzhen, China) according to Illumina's standard procedures, and 2 × 300 bp paired-end reads were generated. The raw reads were filtered to detach adaptors, low-quality, and ambiguous bases to get the tags (25). In addition, The tags were clustered into OTUs using UPARSE software (v7.0.1090) (26). UCHIME (v4.2.40) was used to compare the chimera sequence with the Gold database (27) and classify the representative sequences of OUT by QIIME v1.8.0 (28). All tags were compared with OTU using the USEARCH global (29), and a statistical table of OTU abundance for each sample was obtained. We used MOTHUR (v1.31.2) (30) and QIIME (v1.8.0) (28) to evaluate alpha and beta diversity, respectively. Bar and Heat maps of different taxonomic levels were drawn using R package v3.4.1 and the R package “gplots,” respectively. Linear discriminant analysis effect size (LefSe) was used to conducted LEfSe cluster or LDA analysis. In addition, the Wilcox-test or Kruskal-Test was used as statistical analysis for significant species.
Blood Parameters Measurements
Blood samples were collected for the measurement of biomarkers of immune function and inflammation. The laser flow cytometry was used to measure white blood cell count and the respective subtypes including basophils, eosinophils, neutrophils, lymphocytes, and monocytes. Immunoturbidimetry was used to determine C-reactive protein and immunoglobulin (Ig) G and IgM. The hematology analysis was completed by Adicon Clinical Laboratories, Hangzhou, China. In our study, 0.1 g/L was used as if the value of C-reactive protein had been measured <0.2 g/L occurs (31).
Statistical Analysis
All data analyses were performed using SPSS Statistics version 23.0 (SPSS Inc., Chicago, IL, USA). The Kolmogorov-Smirnov test was used to assess the normal distribution of variables before statistical analysis. One-way analysis of variance with normally distributed data or the Kruskal-Wallis test with non-normally distributed data was used to compare the composition of gut microbiota and blood biomarkers of immune function among elite athletes and general college students. In addition, analysis of covariance (ANCOVA) was used to evaluate whether the difference of gut microbiota and blood immune indicators among participants is independent of dietary status which included weekly intake of cereals, vegetables and fruits, meat, beans and nuts, fatty food, and alcohol. The ANCOVA is usually applied to test whether the independent variable still influences the dependent variable after the influence of the covariate has been removed (32). Pearson's correlation was used to examine the relationship between blood immune or inflammatory biomarkers and gut microbiota composition in all participants. The continuous variables were expressed as mean ± SD or SEM as appropriate, and percentages for categorical variables. A p-value < 0.05 was considered as an indication of a significant difference.
Results
Characteristics of Participants
A total of 66 young adults were enrolled in the present study. Basic information of all participants is shown in Table 1. There was no significant difference in the body mass index among the three groups. However, the Athlete group had higher muscle mass (p = 0.002) and lower fat mass (p = 0.002) than other groups. In terms of physical activity, there was a significant difference in total physical activity per week among the three groups (p < 0.001). As shown in Supplementary Table S1, only the HPA group had job-related physical activity, and the HPA group also participated for a longer total time in transportation and housework physical activity per week. The total time of participating in moderate and high-intensity physical activities in the Athlete group was longer than that in other groups per week (estimated to be higher than in the HPA group by 4-fold and in the LPA group by 150-fold) concerning leisure-time physical activity according to IPAQ. However, the LPA group participated for a longer total time in walking and sitting per week. With regard to physical function, the performance of both hands grip strength (p = 0.008 for the left hand, p = 0.002 for right hand), standing long jump (p < 0.001), sit-and-reach (p < 0.001), one-leg standing (right) with eyes closed (p = 0.007), and Harvard step test (p < 0.001) in the Athlete group were significantly better than the LPA group. Meanwhile, the Athlete group performed significantly better than the HPA group in the standing long jump (p < 0.001) and Harvard step test (p = 0.003). In addition, there were significant differences in the sit-and-reach (p = 0.019) and Harvard step test (p = 0.002) in the HPA group compared with the LPA group. In terms of diet, athletes tended to have higher intakes of almost all food items than the other two groups, except beans and nuts. Specifically, the intake of cereals in the Athlete (p < 0.001) and HPA groups (p = 0.017) was significantly higher than that in the LPA group. Furthermore, the intake of vegetables and fruits, as well as meat in the Athlete group was significantly higher than in the HPA group (p = 0.018, p < 0.001, respectively) and the LPA group (p < 0.001, p < 0.001, respectively).
Gut Microbial Diversity
There were a total of 9,512,790 (9.5 million) 16S rRNA reads from all sixty-three isolated fecal DNA. After quality filtering, 136,836 effective sequences were collected from each fecal sample and the effective rate is 95.80%. After removal of singletons, a total of 11,212 (Athlete = 3,853, HPA = 3,636, LPA = 3,723) operational taxonomic units (OTUs) were identified using 97% sequence similarity. The average number of observed OTUs was 170 per sample. No significant differences in the OTUs were observed among the three groups (p = 0.934).
Alpha diversity is the analysis of species diversity in a single microbiota sample, including observed species index, richness (Chao and Ace index), diversity (Shannon and Simpson index), and good coverage index (30). Table 2 shows that there was no statistically significant difference in alpha diversity among the Athlete, HPA, and LPA groups.
Table 2. Difference in alpha diversity of gut microbiota composition among athletes and young adults.
Regarding beta diversity, no significant differences were observed in the Principal Coordinates Analysis (PCoA) based on unweighted-Unifrac distance metrics (p = 0.845) or weighted-Unifrac distance metrics (p = 0.372) (Supplementary Figure S1). However, the beta weighted-Unifrac heatmap showed a trend of the cluster among the three groups (Supplementary Figure S2). The taxonomy of the three groups of gut microbiomes is different. Linear discriminant analysis effect size (LEfSe) showed that the Athlete group enriched nine taxa, the HPA group enriched one taxon, and the LPA group enriched two taxa (p < 0.05). Specifically, the most abundant microbiota of the Athlete and LPA group was the Firmicutes, while the Proteobacteria in the HPA group (Supplementary Figure S3). These bacterial taxa indicated significant enrichment due to the linear discriminant analysis (LDA) score >2 or <-2 (33).
Gut Microbiota Composition at Phylum and Class Levels
The gut microbiota composition at the phylum level among three groups was represented in Figure 2A. In total, 16 phyla were observed in the intestines of different groups. Of these, 12 phyla were observed in Athlete and HPA groups, while 14 phyla were observed in the LPA group. The abundant phyla were the Bacteroidetes (52.53 in the Athlete group, 55.35 in the HPA group, and 62.81% in the LPA group) and the Firmicutes (43.99 in the Athlete group, 39.67 in the HPA group, and 32.14% in the LPA group) among three groups. However, there were no significant differences in the abundance of the Bacteroidetes and the Firmicutes among the three groups (p = 0.278, p = 0.164, respectively, data not shown). Similarly, the relative abundance of Firmicutes: Bacteroidetes ratio was evaluated, and no difference was observed among the three groups (p = 0.417, data not shown).
Figure 2. Differences of gut microbiota relative abundance at the phylum, class, order, and family levels among three groups. Relative abundance represents the proportion of each species in each group. (A) phylum level; (B) class level; (C) order level; and (D) family level. Athlete, Athlete group; HPA, High physical activity group; LPA, Low physical activity group.
Figure 2B shows the gut microbiota composition at the class level among the three groups. Bacteroidia (52.63 in the Athlete group, 55.34 in the HPA group, and 62.71% in the LPA group), Clostridia (25.68 in the Athlete group, 28.46 in the HPA group, and 24.03% in the LPA group), and Negativicutes (17.84 in the Athlete group, 10.39 in the HPA group, and 7.62% in the LPA group) were the dominant genus among the three groups.
Gut Microbiota Composition at Order and Family Levels
At the order level, Bacteroidales (52.53 in the Athlete group, 55.34 in the HPA group, and 62.71% in the LPA group), Clostridiales (25.76 in the Athlete group, 28.46 in the HPA group, and 24.03% in the LPA group), Selenomonadales (17.84 in the Athlete group, 10.39 in the HPA group, and 7.62% in the LPA group) were the dominant genus among the three groups (Figure 2C). For family level, Prevotellaceae (29.60 in the Athlete group, 25.32 in the HPA group, and 28.61% in the LPA group), Bacteroidaceae (20.88 in the Athlete group, 26.81 in the HPA group, and 31.92% in the LPA group), Ruminococcaceae (12.28 in the Athlete group, 14.67 in the HPA group, and 11.56% in the LPA group) were the dominant genus among the three groups (Figure 2D).
Gut Microbiota Composition at Genus and Species Levels
Figure 3A indicated the gut microbiota composition at the genus level among the three groups. Prevotella (20.88 in Athlete group, 26.81 in the HPA group, and 31.92% in the LPA group), Bacteroides (24.96 in Athlete group, 25.01 in the HPA group, and 27.41% in the LPA group), Faecalibacterium (6.86 in Athlete group, 10.57 in the HPA group and 6.88% in LPA group) and Megamonas (11.67 in the Athlete group, 5.15 in the HPA group and 2.24% in the LPA group) were the dominant genus among three groups.
Figure 3. Differences of gut microbiota relative abundance at the genus and species level among three groups. Relative abundance represents the proportion of each species in each group; (A) genus level; (B) species level; Athlete, Athlete group; HPA,High physical activity group; LPA, Low physical activity group.
At the species level, Prevotella_copri (21.80 in the Athlete group, 24.64 in the HPA group, and 26.11% in the LPA group), Faecalibacterium_prausnitzii (6.86 in the Athlete group, 10.57 in the HPA group, and 6.88% in the LPA group), Bacteroides_vulgatus (7.74 in Athlete group, 6.78 in HPA group, and 6.73% in LPA group), Bacteroides_plebeius (3.32% in the Athlete group, 6.87 in the HPA group, and 8.61% in the LPA group), and Megamonas_funiformis (11.07 in the Athlete group, 5.09 in the HPA group, and 2.11% in the LPA group) were the dominant genus among three groups (Figure 3B).
Differences in Gut Microbiota Composition Among Athletes and Young Adults With Different Physical Activity Levels
A total of 25 taxa of gut microbiota were observed significant differences among the groups in the non-adjusted model. In the Athlete group, the abundance of Lentisphaerae (Victivallaceae belongs to Lentisphaerae; p = 0.045), Clostridiaceae (p = 0.029), Megamonas (p = 0.020), Romboutsia (p = 0.038), Campylobacter_jejuni (p = 0.045), and Paraprevotella_xylaniphila (p = 0.045), were higher. In the LPA group, the abundance of Erysipelotrichia (p = 0.040), Bilophila (p = 0.034), and Faecalicoccus (p = 0.008) were higher. In addition, the abundance of Parasutterella (p = 0.026), Arcobacter_halophilus (p = 0.045), Parasutterella_excrementihominis (p = 0.023) and Clostridium_spiroforme (p = 0.017) were higher in the HPA group.
Moreover, we re-analyzed the bacteria with differences between groups after dietary adjustment, and the results show that there were no group differences in most bacteria after adjusted diet intake. Only the abundance of Clostridiaceae (p = 0.029), Bilophila (p = 0.011), Faecalicoccus (p = 0.050), and Bilophila_wadsworthia (p = 0.011) were still different between groups. In addition, although no differences were observed between the three groups, two bacteria (Campylobacter_jejuni and Megamonas_rupellensis) showed differences among participants. Campylobacter_jejuni was higher in the Athlete group than in the LPA group (p = 0.050), and Megamonas_rupellensis sp. was higher in the Athlete group than in the HPA group (p = 0.029) (Table 3).
Table 3. Phylotypes with significant difference among athletes and young adults before and after adjusting dietary status.
Association of Serum Immune Function Biomarkers With Gut Microbiota
Concerning immune function biomarkers, a significantly higher absolute neutrophil count (p = 0.004) and white blood cell count (p = 0.014) were observed in the LPA group than in other groups in the non-adjusted model. In addition, the HPA group had a higher absolute basophil count (p = 0.022). When adjusting for dietary status, almost all these significant differences disappeared and only a significantly higher absolute count of basophils was observed in the HPA group than in the athlete group (p = 0.037) (Table 4).
Table 4. Blood biomarkers with significant difference among athletes and young adults before and after adjusting dietary status.
Pearson's correlation was conducted to further explore the relationship between gut microbiota, inflammation, and immunity (Figure 4). It was shown that Bilophila genus (Bilophila_wadsworthia) had positive linear correlation with the absolute neutrophil (r = 0.273, p=0.027), lymphocyte count (r = 0.327, p = 0.007), and white blood cell count (r = 0.332, p = 0.006). At the species level, Megamonas_rupellensis was also positively correlated with the absolute lymphocyte count (r = 0.268, p = 0.030).
Figure 4. Association between blood biomarkers of immune function and gut microbiota with Pearson's correlation analysis. Athlete, Athlete group; HPA, High physical activity group; LPA, Low physical activity group. *p < 0.05; **p < 0.01.
Discussion
This matching study investigated the characteristics of gut microbiota at all levels (including phylum, class, order, family, genus, and species) among elite athletes and general young adults with different physical activity levels. Beta diversity tended to cluster between groups with different physical activity levels, while no significant difference was observed in both alpha and beta diversity. Meanwhile, a higher abundance of several bacteria was observed in elite athletes or physically inactive young adults. Moreover, it is inferred that dietary status may be a key factor to be considered when analyzing the association between physical activity and gut microbiota.
In this study, the elite athletes had a significantly higher level of total physical activity and moderate-to-vigorous intensity physical activity, as well as lower sitting time compared to general young adults. It is worth noting that these gut microbiological differences between people with different levels of physical activity have also been observed in previous studies. Specifically, the Bacteroides and Prevotella were more abundant in top Polish endurance athletes compared to sedentary participants (34). The elite athletes and youth non-elite athletes have different taxonomical functional and phenotypic compositions of the gut microbial community. Clostridiales and Faecalibacterium were particularly enriched in the elite athletes (35). Senior orienteers were found to have a more homogeneous microbiota and more abundant Faecalibacterium prausnitzii than the community-dwelling older adults (36). Indeed, our results indicate partial disparity of gut microbiota composition among elite athletes and young adults with different physical activity, particularly after adjusting dietary status.
Diet is an important factor affecting gut microbiota (17). Some previous studies confirmed that the different dietary ingredients have a substantial impact on the gut microbiota (37). Specifically, high-fiber feeding can reshape gut microbiota and promote the release of short-chain fatty acids (SCFAs), which play a role in maintaining the normal functions of the innate and adaptive immune system (38, 39). The intervention of low-calorie weight loss on a high-fat diet induced the growth of bile-resistant bacteria and reduction of bacteria associated with inflammation in humans (40). Moreover, in athletes, high-protein diets may be negatively associated with the diversity of gut microbiota, and a decreased relative abundance of short-chain fatty acid-producing commensal bacteria was observed in a high protein low carbohydrates diet for athletes in resistance sports (41). It is worth noting that athletes often consume a diet that differs from the general population, including in dietary diversity and total energy consumption. The eating habits of athletes are even different between different events (e.g., increased protein intake in resistance-trained athletes or carbohydrate intake in endurance athletes) (42). Previous studies have shown that there is a rapid change in the composition of the gut microbial and increased abundance of Alistipes, Bilophila, and Bacteroides after consuming a high-fat/protein diet for 5 days (43). Thus, there is a need to consider the influence of an athletic diet when comparing the differences in gut microbiota among athletes of different sports, as well as between athletes and the general population. In our study, 19 taxa of gut microbiota showed no statistical differences, the remaining 6 taxa (Clostridiaceae, Bilophila, Faecalicoccus, Bilophila_wadsworthia, Campylobacter_jejuni, Megamonas_rupellensis) observed significant differences among the groups after adjusting the covariates of dietary intake. These changes show that dietary intake is also a particularly important factor in regulating microbiota composition of the gut between elite athletes and general young adults.
Previous studies indicated that some abundant gut microbiota in physically inactive young adults of this study may be closely associated with inflammation and immunity. For example, a high abundance of Bilophila_wadsworthia (belong to Bilophila genus) caused a systemic inflammatory response in mice that included elevated IL-6 (44), which was enriched in patients with Behcet's disease (45) and can promote a Th1-mediated immune response in dietary-fat-induced colitis (46). In addition, a recent study demonstrated that Faecalicoccus was considered to be an important factor in the classification of subjects in Crohn's disease and enriched in patients with Crohn's disease (47). Our study showed that a higher abundance of Bilophila and Faecalicoccus was observed in young adults with insufficient physical activity. It is worth noting that Bilophila genus was positively associated with the counts of circulating white blood cells and its subtypes such as neutrophils and lymphocytes. The bacterium Bilophila_wadsworthia can erode the mucus layer on colon surfaces and allow more bacterial flora to approach lining cells, then resulting in inflammation (48). These findings suggest that physical inactivity may be correlated to pro-inflammatory gut bacteria.
It is well-known that physical inactivity is associated with high-grade systemic inflammation (49) and low immune function (50). However, the mechanism by which physical inactivity induces a reduction of immune function has not been fully elucidated. One possible underlying mechanism is that physical inactivity or low physical activity may lead to impaired ability to store fat and inflammation of subcutaneous adipose tissue (51, 52). In addition, the subjects with low physical activity have a higher susceptibility to infection compared to regular moderate exercise subjects (53). The reason may be that the growth hormone, cortisol, adrenaline, prolactin, and other factors showed a decreased release, leading to attenuated immunomodulatory effects in the subjects with physical inactivity (54).
It is well known that, the white blood cell, comprising neutrophils, eosinophils, basophils, lymphocytes, and monocytes, play an important role in the implementation of immune function and anti-infection in humans (55). The leukocyte count can increase in response to infections. Evidence indicated that patients with severe COVID-19 infection tend to have a higher leukocyte count in the Chinese population (56). Moreover, the elevated leukocyte count, particularly neutrophil and monocyte count were reliable inflammatory markers within the normal range, even in physically healthy individuals (57–59). Thus, it may reduce the body's inflammation and the risk of inflammation-induced diseases if the number of leukocytes can be kept at an appropriate level. Studies have reported regular exercise, particularly in more moderate-to-vigorous physical activity, can have an anti-inflammatory effect and reduce the odds of elevated white blood cell count (51, 60, 61). Although previous studies focused on elite athletes who practice water sports (62), cyclists who participate in intense endurance exercise (63) have higher levels of inflammation and low immune function. The reason for these findings could be that exercise-induced immunosuppression makes elite athletes more susceptible to infection symptoms after short-term acute exercise. Our results showed that athletes participating in regular exercise training have better immune function reflected by lower leukocyte count. Conversely, the Clostridiaceae, a bacterium associated with inflammation (64) was more abundant in athletes. As another bacterium relatively abundant in athletes, Megamonas_rupellensis also showed a positive correlation with circulating lymphocytes. These results may be interpreted as injury-induced local inflammation when the athletes are susceptible to sports injury. Evidence showed that muscle injury in basketball players can lead to increased levels of inflammation (65). Due to the lack of data on acute or chronic sports injury among athletes, these speculations could not be confirmed in this study and need to be further clarified in subsequent work.
We have to admit that there are some limitations to this study. First, our elite athletes only include track and field athletes, and the lack of participants from other sports may affect the results of the gut microbiota in different groups, hence similar findings may not necessarily be observed. Second, there was no data on the nutrient supplementation of participants, particularly in elite athletes, hence we cannot rule out the effect of these factors on main observation results (66). Third, we did not analyze the metabolites of gut microbes, which might normally be directly involved in the human metabolism process (67), hence the in-depth interpretation of results was limited.
Conclusion
These findings indicated that several gut bacteria are inversely associated with immune function and tended to be abundant in physically inactive people, while elite athletes are more likely to have a better immune function and specific gut bacteria composition than general young adults. It is important to note that dietary intakes are potential confounders that may affect the observed association of physical activity with gut microbiota and immune function. Further studies are warranted to determine the correlation between physical activity and gut microbiota. Moreover, considering diet, the effects and mechanisms of long-term habits with moderate and vigorous-intensity physical activity on immune indices and gut microbiota need to be further clarified.
Data Availability Statement
The datasets presented in this article are not readily available without the consent and permission of the study participants. Requests to access the datasets should be directed to Fei Zhong, emYzNTgyMzU3NDMmI3gwMDA0MDt6anUuZWR1LmNu.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Department of Psychology and Behavioral Sciences, Zhejiang University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
YX, XZ, CW, and CH: designing the experiment. YX, FZ, XZ, and H-YL: acquisition of subjects and data. YX, FZ, XZ, H-YL, CW, and CH: statistical analyses and interpretation of data and revision of the manuscript. YX, FZ, XZ, and CH: drafting the manuscript. CH: study supervision. All authors have revised and approved the final version.
Funding
This study was supported by the fundamental research funds (204201*172220192) for the Central Universities, China, and the Hundred Talents Program funding (188020*194221802/004/001) from Zhejiang University, China.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are extremely grateful to all the participants in this study for allowing their data to be published. We would also like to express sincere gratitude to Li Shu and Ping Wang for their assistance with participants' recruitment and data collection.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.843076/full#supplementary-material
Supplementary Figure S1. Principal Coordinate Analysis (PCoA) of gut microbiota. (A) unweighted-Unifrac distance metrics; (B) weighted-Unifrac distance metrics. Each point represents a sample, athlete group in blue, high physical activity group in orange, and low physical activity group in green. The two coordinates are plotted and further indicated the percentage of variability on the axis.
Supplementary Figure S2. Beta weighted-Unifrac heatmap. The Gradient color in the heatmap represents the degree of difference of each sample among athletes, high physical activity group, and low physical activity group. The distance between samples indicated increases gradually with the color from blue to red.
Supplementary Figure S3. Differential taxon features were analyzed by linear discriminant analysis (LDA) effect size (LEfSe) analysis among three groups. (A) showing differentially abundant bacterial taxa. (B) indicating LDA scores. Bacterial taxa indicated significantly high when LDA score >2 or LDA score <-2. Athlete, Athlete group; HPA, High physical activity group; LPA, Low physical activity group.
References
1. Mougeot JLC, Stevens CB, Morton DS, Brennan MT, Mougeot FB. Oral microbiome and cancer therapy-induced oral mucositis. JNCI Monographs. (2019) 2019:lgz002.doi: 10.1093/jncimonographs/lgz002
2. Wu GD, Lewis JD. Analysis of the human gut microbiome and association with disease. clinical gastroenterology and hepatology: the official clinical practice. J Am Gastroenterol Assoc. (2013) 11:38. doi: 10.1016/j.cgh.2013.03.038
3. Scott KP, Gratz SW, Sheridan PO, Flint HJ, Duncan SH. The influence of diet on the gut microbiota. Pharmacol Res. (2013) 69:52–60. doi: 10.1016/j.phrs.2012.10.020
4. Fenneman AC, Rampanelli E, Yin YS, Ames J, Blaser MJ, Fliers E, et al. Gut microbiota and metabolites in the pathogenesis of endocrine disease. Biochem Soc Trans. (2020) 48:915–31. doi: 10.1042/BST20190686
5. Cani PD, Jordan BF. Gut microbiota-mediated inflammation in obesity: a link with gastrointestinal cancer. Nat Rev Gastroenterol Hepatol. (2018) 15:671–82. doi: 10.1038/s41575-018-0025-6
6. Brandsma E, Kloosterhuis NJ, Koster M, Dekker DC, Gijbels MJ, Van Der Velden S, et al. A proinflammatory gut microbiota increases systemic inflammation and accelerates atherosclerosis. Circ Res. (2019) 124:94–100. doi: 10.1161/CIRCRESAHA.118.313234
7. Maslowski KM, Mackay CR. Diet, gut microbiota and immune responses. Nat Immunol. (2011) 12:5–9. doi: 10.1038/ni0111-5
8. Wu HJ, Wu E. The role of gut microbiota in immune homeostasis and autoimmunity. Gut Microbes. (2012) 3:4–14. doi: 10.4161/gmic.19320
9. Al Bander Z, Nitert MD, Mousa A, Naderpoor N. The gut microbiota and inflammation: an overview. Int J Environ Res Public Health. (2020) 17:7618. doi: 10.3390/ijerph17207618
10. Wang J, Wang Y, Zhang X, Liu J, Zhang Q, Zhao Y, et al. Gut microbial dysbiosis is associated with altered hepatic functions and serum metabolites in chronic hepatitis B patients. Front Microbiol. (2017) 8:2222. doi: 10.3389/fmicb.2017.02222
11. Zhou C, Zhao H., Xiao XY, Guo RJ, Wang Q, Chen H, et al. Metagenomic profiling of the pro-inflammatory gut microbiota in ankylosing spondylitis. J Autoimm. (2020) 107:102360. doi: 10.1016/j.jaut.2019.102360
12. Deniz E, Guo R. Needle in a haystack: bacteroides fragilis polysaccharide a as the archetypical symbiosis factor. Ann N Y Acad Sci. (2018) 1417:116–29. doi: 10.1111/nyas.13660
13. Ischander M, Zaldivar F, Eliakim A, Nussbaum E, Dunton G, Leu SY, et al. Physical activity, growth, and inflammatory mediators in bmi-matched female adolescents. Med Sci Sports Exerc. (2007) 39:1131–8. doi: 10.1249/mss.0b013e318053e7a2
14. Slavish DC, Jones DR, Smyth JM, Engeland CG, Song S, McCormick NM, et al. Positive and negative affect and salivary markers of inflammation among young adults. Int J Behav Med. (2020) 27:282–93. doi: 10.1007/s12529-019-09795-2
15. Clarke SF, Murphy EF, O'Sullivan O, Lucey AJ, Humphreys M, Hogan A, et al. Exercise and associated dietary extremes impact on gut microbial diversity. Gut. (2014) 63:1913–20. doi: 10.1136/gutjnl-2013-306541
16. Zhong F, Wen X, Yang M, Lai HY, Momma H, Cheng L, et al. Effect of an 8-week exercise training on gut microbiota in physically inactive older women. Int J Sports Med. (2021) 42:610–23. doi: 10.1055/a-1301-7011
17. Zmora N, Suez J, Elinav E. You are what you eat: diet, health and the gut microbiota. Nat Rev Gastroenterol Hepatol. (2019) 16:35–56. doi: 10.1038/s41575-018-0061-2
18. Bolte LA, Vila AV, Imhann F, Collij V, Gacesa R, Peters V, et al. Long-term dietary patterns are associated with pro-inflammatory and anti-inflammatory features of the gut microbiome. Gut. (2021) 70:1287–98. doi: 10.1136/gutjnl-2020-322670
19. Bibiloni MDM, Julibert A, Argelich E, Aparicio-Ugarriza R, Palacios G, Pons A, et al. Western and mediterranean dietary patterns and physical activity and fitness among spanish older adults. Nutrients. (2017) 9:704. doi: 10.3390/nu9070704
20. Maher JP, Harduk M, Hevel DJ, Adams WM, McGuirt JT. Momentary physical activity co-occurs with healthy and unhealthy dietary intake in African American college freshmen. Nutrients. (2020) 12:1360. doi: 10.3390/nu12051360
21. Zhu Q, Jiang S, Du G. Effects of exercise frequency on the gut microbiota in elderly individuals. Microbiologyopen. (2020) 9:e1053. doi: 10.1002/mbo3.1053
22. Piercy KL, Troiano RP, Ballard RM, Carlson SA, Fulton JE, Galuska DA, et al. The physical activity guidelines for Americans. Jama. (2018) 320:2020–8. doi: 10.1001/jama.2018.14854
23. Craig CL, Marshall AL, Sjöström M, Bauman AE, Booth ML, Ainsworth BE, et al. International physical activity questionnaire: 12-country reliability and validity. Med Sci Sports Exerc. (2003) 35:1381–95. doi: 10.1249/01.MSS.0000078924.61453.FB
24. Gao J. Fei JQ, Jiang LJ, Yao WQ, Lin B, Guo HW. Assessment of the reproducibility and validity of a simple food-frequency questionnaire used in dietary pattern studies. Acta Nutr Sin. (2011) 33:452–6. doi: 10.13325/j.cnki.acta.nutr.sin.2011.05.012
25. Magoč T, Salzberg SL. Flash: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. (2011) 27:2957–63. doi: 10.1093/bioinformatics/btr507
26. Edgar RC. Uparse: highly accurate otu sequences from microbial amplicon reads. Nat Methods. (2013) 10:996–8. doi: 10.1038/nmeth.2604
27. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. Uchime improves sensitivity and speed of chimera detection. Bioinformatics. (2011) 27:2194–200. doi: 10.1093/bioinformatics/btr381
28. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. Qiime allows analysis of high-throughput community sequencing data. Nat Methods. (2010) 7:335–6. doi: 10.1038/nmeth.f.303
29. Edgar RC. Search and clustering orders of magnitude faster than blast. Bioinformatics. (2010) 26:2460–1. doi: 10.1093/bioinformatics/btq461
30. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. (2009) 75:7537–41. doi: 10.1128/AEM.01541-09
31. Helsel DR. Fabricating data: how substituting values for non-detects can ruin results, and what can be done about it. Chemosphere. (2006) 65:2434–9. doi: 10.1016/j.chemosphere.2006.04.051
32. Mishra P, Singh U, Pandey CM, Mishra P, Pandey G. application of student's t-test, analysis of variance, and covariance. Ann Card Anaesth. (2019) 22:407–11. doi: 10.4103/aca.ACA_94_19
33. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. (2011) 12:R60. doi: 10.1186/gb-2011-12-6-r60
34. Kulecka M, Fraczek B, Mikula M, Zeber-Lubecka N, Karczmarski J, Paziewska A, et al. The composition and richness of the gut microbiota differentiate the top polish endurance athletes from sedentary controls. Gut Microbes. (2020) 11:1374–84. doi: 10.1080/19490976.2020.1758009
35. Han M, Yang K, Yang P, Zhong C, Chen C, Wang S, et al. Stratification of athletes' gut microbiota: the multifaceted hubs associated with dietary factors, physical characteristics, and performance. Gut Microbes. (2020) 12:1–18. doi: 10.1080/19490976.2020.1842991
36. Fart F, Rajan SK, Wall R, Rangel I, Ganda-Mall JP, Tingo L, et al. Differences in gut microbiome composition between senior orienteering athletes and community-dwelling older adults. Nutrients. (2020) 12:2610. doi: 10.3390/nu12092610
37. Zhang N, Ju Z, Zuo T. Time for food: the impact of diet on gut microbiota and human health. Nutrition. (2018) 51:80–5. doi: 10.1016/j.nut.2017.12.005
38. Tan J, McKenzie C, Vuillermin PJ, Goverse G, Vinuesa CG, Mebius RE, et al. Dietary fiber and bacterial scfa enhance oral tolerance and protect against food allergy through diverse cellular pathways. Cell Rep. (2016) 15:2809–24. doi: 10.1016/j.celrep.2016.05.047
39. Makki K, Deehan EC, Walter J, Backhed F. The impact of dietary fiber on gut microbiota in host health and disease. Cell Host Microbe. (2018) 23:705–15. doi: 10.1016/j.chom.2018.05.012
40. Jaagura M, Viiard E, Karu-Lavits K, Adamberg K. Low-carbohydrate high-fat weight reduction diet induces changes in human gut microbiota. Microbiologyopen. (2021) 10:e1194. doi: 10.1002/mbo3.1194
41. Jang LG, Choi G, Kim SW, Kim BY, Lee S, Park H. The combination of sport and sport-specific diet is associated with characteristics of gut microbiota: an observational study. J Int Soc Sports Nutr. (2019) 16:1–10. doi: 10.1186/s12970-019-0290-y
42. Mohr AE, Jäger R, Carpenter KC, Kerksick CM, Purpura M, Townsend JR, et al. The athletic gut microbiota. J Int Soc Sports Nutr. (2020) 17:1–33. doi: 10.1186/s12970-020-00353-w
43. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. (2014) 505:559–63. doi: 10.1038/nature12820
44. Feng Z, Long W, Hao B, Ding D, Ma X, Zhao L, et al. A human stool-derived bilophila wadsworthia strain caused systemic inflammation in specific-pathogen-free mice. Gut Pathog. (2017) 9:59. doi: 10.1186/s13099-017-0208-7
45. Ye Z, Zhang N, Wu C, Zhang X, Wang Q, Huang X, et al. A metagenomic study of the gut microbiome in Behcet's disease. Microbiome. (2018) 6:135. doi: 10.1186/s40168-018-0520-6
46. Devkota S, Wang Y, Musch MW, Leone V, Fehlner-Peach H, Nadimpalli A, et al. Dietary-fat-induced taurocholic acid promotes pathobiont expansion and colitis in Il10-/- mice. Nature. (2012) 487:104–8. doi: 10.1038/nature11225
47. Forbes JD, Chen CY, Knox NC, Marrie RA, El-Gabalawy H, de Kievit T, et al. A comparative study of the gut microbiota in immune-mediated inflammatory diseases-does a common dysbiosis exist? Microbiome. (2018) 6:221. doi: 10.1186/s40168-018-0603-4
48. Ribaldone DG, Pellicano R, Actis GC. Inflammation: a highly conserved, janus-like phenomenon—a gastroenterologist'perspective. J Mol Med. (2018) 96:861–71. doi: 10.1007/s00109-018-1668-z
49. Geffken DF, Cushman M, Burke GL, Polak JF, Sakkinen PA, Tracy RP. Association between physical activity and markers of inflammation in a healthy elderly population. Am J Epidemiol. (2001) 153:242–50. doi: 10.1093/aje/153.3.242
50. Brolinson PG, Elliott D. Exercise and the immune system. Clin Sports Med. (2007) 26:311–9. doi: 10.1016/j.csm.2007.04.011
51. Petersen AMW, Pedersen BK. The anti-inflammatory effect of exercise. J Appl Physiol. (2005) 98:1154–62. doi: 10.1152/japplphysiol.00164.2004
52. Pedersen BK. The diseasome of physical inactivity–and the role of myokines in muscle–fat cross talk. J Physiol. (2009) 587(Pt 23):5559–68. doi: 10.1113/jphysiol.2009.179515
53. Booth FW, Roberts CK, Laye MJ. Lack of Exercise is a major cause of chronic diseases. Compr Physiol. (2012) 2:1143. doi: 10.1002/cphy.c110025
54. Handschin C, Spiegelman BM. The role of exercise and Pgc1α in inflammation and chronic disease. Nature. (2008) 454:463–9. doi: 10.1038/nature07206
55. Ley K, Laudanna C, Cybulsky MI, Nourshargh S. Getting to the site of inflammation: the leukocyte adhesion cascade updated. Nat Rev Immunol. (2007) 7:678–89. doi: 10.1038/nri2156
56. Qin C, Zhou L, Hu Z, Zhang S, Yang S, Tao Y, et al. Dysregulation of immune response in patients with coronavirus 2019 (Covid-19) in Wuhan, China. Clin Infect Dis. (2020) 71:762–8. doi: 10.1093/cid/ciaa248
57. Chmielewski PP, Boryslawski K, Chmielowiec K, Chmielowiec J, Strzelec B. The association between total leukocyte count and longevity: evidence from longitudinal and cross-sectional data. Ann Anat. (2016) 204:1–10. doi: 10.1016/j.aanat.2015.09.002
58. Mehta JL, Saldeen TG, Rand K. Interactive role of infection, inflammation and traditional risk factors in atherosclerosis and coronary artery disease. J Am Coll Cardiol. (1998) 31:1217–25. doi: 10.1016/S0735-1097(98)00093-X
59. Chmielewski PP, Strzelec B. Elevated leukocyte count as a harbinger of systemic inflammation, disease progression, and poor prognosis: a review. Folia Morphol. (2018) 77:171–8. doi: 10.5603/FM.a2017.0101
60. King DE, Carek P, Mainous AG III, Pearson WS. Inflammatory markers and exercise: differences related to exercise type. Med Sci Sports Exerc. (2003) 35:575–81. doi: 10.1249/01.MSS.0000058440.28108.CC
61. Loprinzi PD, Richart SM. White blood cell counts mediate the effects of physical activity on prostate-specific antigen levels. Res Q Exerc Sport. (2014) 85:409–13. doi: 10.1080/02701367.2014.930407
62. Belda J, Ricart S, Casan P, Giner J, Bellido-Casado J, Torrejon M, et al. Airway inflammation in the elite athlete and type of sport. Br J Sports Med. (2008) 42:244–8. doi: 10.1136/bjsm.2007.036335
63. Nieman DC, Konrad M, Henson DA, Kennerly K, Shanely RA, Wallner-Liebmann SJ. Variance in the acute inflammatory response to prolonged cycling is linked to exercise intensity. J Interferon Cytokine Res. (2012) 32:12–7. doi: 10.1089/jir.2011.0038
64. Muniz Pedrogo DA, Chen J, Hillmann B, Jeraldo P, Al-Ghalith G, Taneja V, et al. An increased abundance of clostridiaceae characterizes arthritis in inflammatory bowel disease and rheumatoid arthritis: a cross-sectional study. Inflamm Bowel Dis. (2019) 25:902–13. doi: 10.1093/ibd/izy318
65. Luti S, Modesti A, Modesti PA. Inflammation, peripheral signals and redox homeostasis in athletes who practice different sports. Antioxidants. (2020) 9:65. doi: 10.3390/antiox9111065
66. Lin CL, Hsu YJ, Ho HH, Chang YC, Kuo YW, Yeh YT, et al. Bifidobacterium longum subsp. longum olp-01 supplementation during endurance running training improves exercise performance in middle- and long-distance runners: a double-blind controlled trial. Nutrients. (2020) 12:1972. doi: 10.3390/nu12071972
Keywords: gut microbiota, dietary status, physical activity, elite athlete, inflammation, matching study
Citation: Xu Y, Zhong F, Zheng X, Lai H-Y, Wu C and Huang C (2022) Disparity of Gut Microbiota Composition Among Elite Athletes and Young Adults With Different Physical Activity Independent of Dietary Status: A Matching Study. Front. Nutr. 9:843076. doi: 10.3389/fnut.2022.843076
Received: 24 December 2021; Accepted: 27 January 2022;
Published: 18 March 2022.
Edited by:
Jianhua Wu, University of Leeds, United KingdomCopyright © 2022 Xu, Zhong, Zheng, Lai, Wu and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Cong Huang, Y29odWFuZyYjeDAwMDQwO3pqdS5lZHUuY24=
†These authors have contributed equally to this work