Front. Nutr. Frontiers in Nutrition Front. Nutr. 2296-861X Frontiers Media S.A. 10.3389/fnut.2022.843076 Nutrition Original Research Disparity of Gut Microbiota Composition Among Elite Athletes and Young Adults With Different Physical Activity Independent of Dietary Status: A Matching Study Xu Yongjin 1 † Zhong Fei 1 † Zheng Xiaoqian 1 Lai Hsin-Yi 2 3 4 Wu Chunchun 1 Huang Cong 1 5 * 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

Edited by: Jianhua Wu, University of Leeds, United Kingdom

Reviewed by: Huan Li, Lanzhou University, China; Francesco Suriano, University of Gothenburg, Sweden

 *Correspondence: Cong Huang cohuang@zju.edu.cn

This article was submitted to Nutritional Immunology, a section of the journal Frontiers in Nutrition

†These authors have contributed equally to this work

 18 03 2022 2022 9 843076 24 12 2021 27 01 2022 Copyright © 2022 Xu, Zhong, Zheng, Lai, Wu and Huang. 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.

 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.

 gut microbiota dietary status physical activity elite athlete inflammation matching study 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.

Schematic presentation of the flow of participants through screening.

 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).

Characteristics of the participants.

Characteristics Athlete group HPA group LPA group
(n = 22) (n = 22) (n = 22)
Age, year (SD) 21.55 ± 2.42 21.27 ± 2.47 21.27 ± 2.51
Male (%) 12 (54.5) 12 (54.5) 12 (54.5)
BMI, kg/m2 (SD) 20.67 ± 5.26 22.28 ± 2.71 21.89 ± 2.81
PA, METs· hour/week 138.50 ± 108.26 59.13 ± 29.96 20.25 ± 13.17
Muscle mass, kg (SD) 27.45 ± 6.51 23.31 ± 6.44 20.52 ± 6.02
Fat mass, kg (SD) 9.70 ± 4.41 12.99 ± 4.76 15.06 ± 6.72
Fat percentage, % (SD) 14.30 ± 5.26 20.47 ± 6.67 24.08 ± 8.46
Grip strength, kg (SD)
  Left 37.34 ± 9.81 33.41 ± 10.42 29.22 ± 9.33
  Right 39.15 ± 10.58 34.84 ± 11.70 29.71 ± 7.77
  Standing long jump, m (SD) 2.46 ± 0.29 1.89 ± 0.32 1.77 ± 0.34
  Sit-and-reach, cm (SD) 18.26 ± 8.68 14.54 ± 9.07 8.33 ± 7.82
One-leg standing with eyes closed, sec. (SD)
  Left 30.42 ± 30.07 24.67 ± 23.21 24.39 ± 26.41
  Right 35.13 ± 31.24 24.72 ± 23.27 15.57 ± 10.36
Harvard step test (SD)
  Completed, n (%) 21 (95.45) 15 (68.18) 7 (31.82)
  Harvard step index 96.56 ± 18.43 76.75 ± 18.81 56.39 ± 24.75
Daily diet, g/d (SD)
  Cereals 164.07 ± 142.82 122.70 ± 111.51 44.09 ± 29.50
  Vegetables and fruits 320.38 ± 225.33 201.82 ± 128.26 120.08 ± 107.19
  Meat 642.65 ± 438.48 319.17 ± 207.61 390.34 ± 339.50
  Beans and nuts 31.06 ± 36.27 18.41 ± 22.43 34.55 ± 46.32
  Fatty food 62.35 ± 74.98 39.47 ± 49.02 37.95 ± 45.93
  Alcohol 41.36 ± 99.61 10.98 ± 29.67 24.31 ± 101.23

PA, physical activity; HPA, high physical activity; LPA, low physical activity; BMI, body mass index; SD, Standard deviation.

 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.

Difference in alpha diversity of gut microbiota composition among athletes and young adults.

Athlete group HPA group LPA group
(n = 22) (n = 22) (n = 22)
Alpha diversity
mean ±SEM mean ±SEM mean ±SEM p-value
sobs† 175.14 ± 12.50 165.27 ± 8.30 169.23 ± 10.29 0.934
chao† 202.94 ± 13.55 197.07 ± 9.56 200.53 ± 13.24 0.925
Ace† 207.91 ± 13.63 196.60 ± 9.16 194.87 ± 12.59 0.712
Shannon† 2.52 ± 0.16 2.60 ± 0.15 2.64 ± 0.14 0.924
Simpson† 0.20 ± 0.03 0.19 ± 0.03 0.19 ± 0.03 0.936
Coverage‡ 1.00 ± 0.00 1.00 ± 0.00 1.00 ± 0.00 0.744

HPA, high physical activity; LPA, low physical activity. SEM, Standard error of the mean.

Kruskal–Wallis test;

One-way analysis of variance.

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).

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.

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).

Phylotypes with significant difference among athletes and young adults before and after adjusting dietary status.

Non-adjusted model† Dietary status adjusted model‡
Taxa specific Athlete group HPA group LPA group Athlete group HPA group LPA group
(n = 22) (n = 22) (n = 22) (n = 22) (n = 22) (n =22)
Mean ±SEM Mean ±SEM Mean ±SEM p-value Mean ±SEM Mean ±SEM Mean ±SEM p-value
Phylum Lentisphaerae 0.004 ± 0.004 0.000 ± 0.000 0.000 ± 0.000 0.045 0.005 ± 0.002 0.000 ± 0.002 −0.001 ± 0.002 0.242
Class Erysipelotrichia 0.143 ± 0.036 0.110 ± 0.028 0.221 ± 0.051# 0.040 0.129 ± 0.045 0.115 ± 0.042 0.231 ± 0.044 0.148
Lentisphaeria 0.004 ± 0.004 0.000 ± 0.000 0.000 ± 0.000 0.045 0.005 ± 0.002 0.000 ± 0.002 −0.001 ± 0.002 0.242
Order Erysipelotrichales 0.143 ± 0.036 0.110 ± 0.028 0.221 ± 0.051# 0.040 0.129 ± 0.045 0.115 ± 0.042 0.231 ± 0.044 0.148
Victivallales 0.004 ± 0.004 0.000 ± 0.000 0.000 ± 0.000 0.045 0.005 ± 0.002 0.000 ± 0.002 −0.001 ± 0.002 0.242
Family Clostridiaceae 0.144 ± 0.070 0.023 ± 0.008* 0.026 ± 0.014 0.029 0.180 ± 0.049 0.014 ± 0.046* −0.002 ± 0.048* 0.029
Erysipelotrichaceae 0.143 ± 0.036 0.110 ± 0.028 0.221 ± 0.051# 0.040 0.129 ± 0.045 0.115 ± 0.042 0.231 ± 0.044 0.148
Victivallaceae 0.004 ± 0.004 0.000 ± 0.000 0.000 ± 0.000 0.045 0.005 ± 0.002 0.000 ± 0.002 −0.001 ± 0.002 0.242
Genus Bilophila 0.100 ± 0.032 0.101 ± 0.029 0.249 ± 0.053 0.034 0.113 ± 0.044 0.076 ± 0.041 0.261 ± 0.043*## 0.011
Faecalicoccus 0.037 ± 0.024 0.005 ± 0.002 0.092 ± 0.038## 0.008 0.022 ± 0.028 0.010 ± 0.026 0.103 ± 0.028* 0.050
Megamonas 11.670 ± 3.898 5.154 ± 4.276* 2.239 ± 0.776 0.020 11.254 ± 3.871 6.970 ± 3.622 0.839 ± 3.794 0.206
Parasutterella 0.257 ± 0.144 0.969 ± 0.276* 0.703 ± 0.344 0.026 0.360 ± 0.316 1.029 ± 0.295 0.539 ± 0.309 0.290
Romboutsia 0.134 ± 0.034 0.062 ± 0.017* 0.111 ± 0.056 0.038 0.127 ± 0.046 0.071 ± 0.043 0.109 ± 0.045 0.661
Victivallis 0.004 ± 0.004 0.000 ± 0.000 0.000 ± 0.000 0.045 0.005 ± 0.002 0.000 ± 0.002 −0.001 ± 0.002 0.242
Species Arcobacter_halophilus 0.000 ± 0.000 0.001 ± 0.001 0.000 ± 0.000 0.045 0.000 ± 0.000 0.001 ± 0.000 0.000 ± 0.000 0.350
Bilophila_wadsworthia 0.100 ± 0.032 0.101 ± 0.029 0.249 ± 0.053 0.034 0.113 ± 0.044 0.076 ± 0.041 0.261 ± 0.043*## 0.011
Campylobacter_jejuni 0.046 ± 0.045 0.000 ± 0.000 0.000 ± 0.000 0.045 0.072 ± 0.030 −0.007 ± 0.028 −0.020 ± 0.030* 0.103
Clostridium_spiroforme 0.000 ± 0.000 0.002 ± 0.001* 0.000 ± 0.000 0.017 0.000 ± 0.001 0.002 ± 0.001 0.000 ± 0.001 0.139
Megamonas_funiformis 11.065 ± 3.721 5.092 ± 4.240* 2.109 ± 0.733 0.019 10.661 ± 3.773 6.879 ± 3.530 0.726 ± 3.698 0.216
Megamonas_rupellensis 0.125 ± 0.046 0.004 ± 0.003* 0.041 ± 0.032 0.036 0.125 ± 0.038 0.005 ± 0.035* 0.039 ± 0.037 0.087
Paraprevotella_xylaniphila 0.100 ± 0.062 0.000 ± 0.000 0.000 ± 0.000 0.045 0.038 ± 0.035 −0.015 ± 0.033 0.077 ± 0.034 0.167
Parasutterella_excrementihominis 0.253 ± 0.144 0.969 ± 0.276* 0.699 ± 0.344 0.023 0.357 ± 0.316 1.030 ± 0.295 0.533 ± 0.309 0.285
Romboutsia_sedimentorum 0.134 ± 0.034 0.062 ± 0.017* 0.111 ± 0.056 0.038 0.127 ± 0.046 0.071 ± 0.043 0.109 ± 0.045 0.661
Rothia_aeria 0.000 ± 0.000 0.004 ± 0.000 0.000 ± 0.000 0.045 0.000 ± 0.000 0.000 ± 0.000 0.000 ± 0.000 0.166
Victivallis_vadensis 0.004 ± 0.003 0.000 ± 0.000 0.000 ± 0.000 0.045 0.005 ± 0.002 0.000 ± 0.002 −0.001 ± 0.002 0.242

HPA, high physical activity; LPA, low physical activity; SEM, Standard error of the mean.

p < 0.05 when HPA or LPA group compared with Athlete group.

p < 0.05 when LPA group compared with HPA group;

p < 0.01 when LPA group compared with HPA group.

Kruskal–Wallis test;

Analysis of covariance.

 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).

Blood biomarkers with significant difference among athletes and young adults before and after adjusting dietary status.

Non-adjusted model Dietary status adjusted modelΨ
Blood parameters Athlete group HPA group LPA group Athlete group HPA group LPA group
(n = 22) (n = 22) (n = 22) (n = 22) (n = 22) (n =22)
Mean ±SEM Mean ±SEM Mean ±SEM p-value Mean ±SEM Mean ±SEM Mean ±SEM p-value
White body cell, (109/L)† 4.96 ± 0.20 5.53 ± 0.23 6.45 ± 0.41* 0.014 5.18 ± 0.34 5.57 ± 0.32 6.19 ± 0.33 0.148
Basophils, (109/L)† 0.027 ± 0.002 0.039 ± 0.003* 0.034 ± 0.003 0.022 0.028 ± 0.003 0.038 ± 0.003* 0.034 ± 0.003 0.111
Eosinophils, (109/L)† 0.136 ± 0.020 0.136 ± 0.015 0.142 ± 0.019 0.886 0.130 ± 0.022 0.141 ± 0.020 0.144 ± 0.021 0.902
Neutrophils, (109/L)‡ 2.30 ± 0.15 2.93 ± 0.18* 3.32 ± 0.27* 0.004 2.46 ± 0.23 2.94 ± 0.22 3.15 ± 0.23 0.147
Lymphocytes, (109/L)† 2.19 ± 0.10 2.10 ± 0.10 2.62 ± 0.17 0.062 2.25 ± 0.15 2.13 ± 0.14 2.54 ± 0.15 0.151
Monocytes, (109/L)‡ 0.303 ± 0.016 0.326 ± 0.022 0.343 ± 0.023 0.392 0.311 ± 0.024 0.330 ± 0.022 0.331 ± 0.023 0.818
Immunoglobulin G, (g/L)‡ 11.29 ± 0.26 12.07 ± 0.37 11.97 ± 0.47 0.285 11.18 ± 0.44 12.14 ± 0.41 12.00 ± 0.43 0.279
Immunoglobulin M, (g/L)‡ 1.46 ± 0.09 1.34 ± 0.09 1.52 ± 0.12 0.412 1.46 ± 0.12 1.34 ± 0.11 1.53 ± 0.12 0.490
C-reactive protein, (g/L)† 0.81 ± 0.22 1.14 ± 0.53 0.74 ± 0.10 0.122 0.97 ± 0.39 1.27 ± 0.37 0.45 ± 0.39 0.324

HPA, high physical activity; LPA, low physical activity; SEM, Standard error of the mean.

p < 0.05 when HPA or LPA group compared with Athlete group.

Kruskal–Wallis test in non-adjusted model;

One-way analysis of variance in non-adjusted model;

Analysis of covariance in dietary status adjusted model.

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).

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, zf358235743@zju.edu.cn.

 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.

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

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.

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.

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 Mougeot JLC Stevens CB Morton DS Brennan MT Mougeot FB. Oral microbiome and cancer therapy-induced oral mucositis. JNCI Monographs. (2019) 2019:lgz002.10.1093/jncimonographs/lgz00231425594 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. 10.1016/j.cgh.2013.03.03823643636 Scott KP Gratz SW Sheridan PO Flint HJ Duncan SH. The influence of diet on the gut microbiota. Pharmacol Res. (2013) 69:52–60. 10.1016/j.phrs.2012.10.02023147033 Fenneman AC Rampanelli E Yin YS Ames J Blaser MJ Fliers E . Gut microbiota and metabolites in the pathogenesis of endocrine disease. Biochem Soc Trans. (2020) 48:915–31. 10.1042/BST2019068632412045 Cani PD Jordan BF. Gut microbiota-mediated inflammation in obesity: a link with gastrointestinal cancer. Nat Rev Gastroenterol Hepatol. (2018) 15:671–82. 10.1038/s41575-018-0025-629844585 Brandsma E Kloosterhuis NJ Koster M Dekker DC Gijbels MJ Van Der Velden S . A proinflammatory gut microbiota increases systemic inflammation and accelerates atherosclerosis. Circ Res. (2019) 124:94–100. 10.1161/CIRCRESAHA.118.31323430582442 Maslowski KM Mackay CR. Diet, gut microbiota and immune responses. Nat Immunol. (2011) 12:5–9. 10.1038/ni0111-533435203 Wu HJ Wu E. The role of gut microbiota in immune homeostasis and autoimmunity. Gut Microbes. (2012) 3:4–14. 10.4161/gmic.1932022356853 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. 10.3390/ijerph1720761833086688 Wang J Wang Y Zhang X Liu J Zhang Q Zhao Y . Gut microbial dysbiosis is associated with altered hepatic functions and serum metabolites in chronic hepatitis B patients. Front Microbiol. (2017) 8:2222. 10.3389/fmicb.2017.0222229180991 Zhou C Zhao H. Xiao XY Guo RJ Wang Q Chen H . Metagenomic profiling of the pro-inflammatory gut microbiota in ankylosing spondylitis. J Autoimm. (2020) 107:102360. 10.1016/j.jaut.2019.10236031806420 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. 10.1111/nyas.1366029528123 Ischander M Zaldivar F Eliakim A Nussbaum E Dunton G Leu SY . Physical activity, growth, and inflammatory mediators in bmi-matched female adolescents. Med Sci Sports Exerc. (2007) 39:1131–8. 10.1249/mss.0b013e318053e7a217596781 Slavish DC Jones DR Smyth JM Engeland CG Song S McCormick NM . Positive and negative affect and salivary markers of inflammation among young adults. Int J Behav Med. (2020) 27:282–93. 10.1007/s12529-019-09795-231222631 Clarke SF Murphy EF O'Sullivan O Lucey AJ Humphreys M Hogan A . Exercise and associated dietary extremes impact on gut microbial diversity. Gut. (2014) 63:1913–20. 10.1136/gutjnl-2013-30654125021423 Zhong F Wen X Yang M Lai HY Momma H Cheng L . Effect of an 8-week exercise training on gut microbiota in physically inactive older women. Int J Sports Med. (2021) 42:610–23. 10.1055/a-1301-701133321523 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. 10.1038/s41575-018-0061-230262901 Bolte LA Vila AV Imhann F Collij V Gacesa R Peters V . Long-term dietary patterns are associated with pro-inflammatory and anti-inflammatory features of the gut microbiome. Gut. (2021) 70:1287–98. 10.1136/gutjnl-2020-32267033811041 Bibiloni MDM Julibert A Argelich E Aparicio-Ugarriza R Palacios G Pons A . Western and mediterranean dietary patterns and physical activity and fitness among spanish older adults. Nutrients. (2017) 9:704. 10.3390/nu907070428684703 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. 10.3390/nu1205136032397433 Zhu Q Jiang S Du G. Effects of exercise frequency on the gut microbiota in elderly individuals. Microbiologyopen. (2020) 9:e1053. 10.1002/mbo3.105332356611 Piercy KL Troiano RP Ballard RM Carlson SA Fulton JE Galuska DA . The physical activity guidelines for Americans. Jama. (2018) 320:2020–8. 10.1001/jama.2018.1485430418471 Craig CL Marshall AL Sjöström M Bauman AE Booth ML Ainsworth BE . International physical activity questionnaire: 12-country reliability and validity. Med Sci Sports Exerc. (2003) 35:1381–95. 10.1249/01.MSS.0000078924.61453.FB12900695 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. 10.13325/j.cnki.acta.nutr.sin.2011.05.012 Magoč T Salzberg SL. Flash: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. (2011) 27:2957–63. 10.1093/bioinformatics/btr50721903629 Edgar RC. Uparse: highly accurate otu sequences from microbial amplicon reads. Nat Methods. (2013) 10:996–8. 10.1038/nmeth.260423955772 Edgar RC Haas BJ Clemente JC Quince C Knight R. Uchime improves sensitivity and speed of chimera detection. Bioinformatics. (2011) 27:2194–200. 10.1093/bioinformatics/btr38121700674 Caporaso JG Kuczynski J Stombaugh J Bittinger K Bushman FD Costello EK . Qiime allows analysis of high-throughput community sequencing data. Nat Methods. (2010) 7:335–6. 10.1038/nmeth.f.30320383131 Edgar RC. Search and clustering orders of magnitude faster than blast. Bioinformatics. (2010) 26:2460–1. 10.1093/bioinformatics/btq46120709691 Schloss PD Westcott SL Ryabin T Hall JR Hartmann M Hollister EB . Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. (2009) 75:7537–41. 10.1128/AEM.01541-0919801464 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. 10.1016/j.chemosphere.2006.04.05116737727 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. 10.4103/aca.ACA_94_1931621677 Segata N Izard J Waldron L Gevers D Miropolsky L Garrett WS . Metagenomic biomarker discovery and explanation. Genome Biol. (2011) 12:R60. 10.1186/gb-2011-12-6-r6021702898 Kulecka M Fraczek B Mikula M Zeber-Lubecka N Karczmarski J Paziewska A . The composition and richness of the gut microbiota differentiate the top polish endurance athletes from sedentary controls. Gut Microbes. (2020) 11:1374–84. 10.1080/19490976.2020.175800932401138 Han M Yang K Yang P Zhong C Chen C Wang S . Stratification of athletes' gut microbiota: the multifaceted hubs associated with dietary factors, physical characteristics, and performance. Gut Microbes. (2020) 12:1–18. 10.1080/19490976.2020.184299133289609 Fart F Rajan SK Wall R Rangel I Ganda-Mall JP Tingo L . Differences in gut microbiome composition between senior orienteering athletes and community-dwelling older adults. Nutrients. (2020) 12:2610. 10.3390/nu1209261032867153 Zhang N Ju Z Zuo T. Time for food: the impact of diet on gut microbiota and human health. Nutrition. (2018) 51:80–5. 10.1016/j.nut.2017.12.00529621737 Tan J McKenzie C Vuillermin PJ Goverse G Vinuesa CG Mebius RE . Dietary fiber and bacterial scfa enhance oral tolerance and protect against food allergy through diverse cellular pathways. Cell Rep. (2016) 15:2809–24. 10.1016/j.celrep.2016.05.04727332875 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. 10.1016/j.chom.2018.05.01229902436 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. 10.1002/mbo3.119434180599 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. 10.1186/s12970-019-0290-y31053143 Mohr AE Jäger R Carpenter KC Kerksick CM Purpura M Townsend JR . The athletic gut microbiota. J Int Soc Sports Nutr. (2020) 17:1–33. 10.1186/s12970-020-00353-w32398103 David LA Maurice CF Carmody RN Gootenberg DB Button JE Wolfe BE . Diet rapidly and reproducibly alters the human gut microbiome. Nature. (2014) 505:559–63. 10.1038/nature1282024336217 Feng Z Long W Hao B Ding D Ma X Zhao L . A human stool-derived bilophila wadsworthia strain caused systemic inflammation in specific-pathogen-free mice. Gut Pathog. (2017) 9:59. 10.1186/s13099-017-0208-729090023 Ye Z Zhang N Wu C Zhang X Wang Q Huang X . A metagenomic study of the gut microbiome in Behcet's disease. Microbiome. (2018) 6:135. 10.1186/s40168-018-0520-630077182 Devkota S Wang Y Musch MW Leone V Fehlner-Peach H Nadimpalli A . Dietary-fat-induced taurocholic acid promotes pathobiont expansion and colitis in Il10-/- mice. Nature. (2012) 487:104–8. 10.1038/nature1122522722865 Forbes JD Chen CY Knox NC Marrie RA El-Gabalawy H de Kievit T . A comparative study of the gut microbiota in immune-mediated inflammatory diseases-does a common dysbiosis exist? Microbiome. (2018) 6:221. 10.1186/s40168-018-0603-430545401 Ribaldone DG Pellicano R Actis GC. Inflammation: a highly conserved, janus-like phenomenon—a gastroenterologist'perspective. J Mol Med. (2018) 96:861–71. 10.1007/s00109-018-1668-z29987405 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. 10.1093/aje/153.3.24211157411 Brolinson PG Elliott D. Exercise and the immune system. Clin Sports Med. (2007) 26:311–9. 10.1016/j.csm.2007.04.01117826186 Petersen AMW Pedersen BK. The anti-inflammatory effect of exercise. J Appl Physiol. (2005) 98:1154–62. 10.1152/japplphysiol.00164.200415772055 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. 10.1113/jphysiol.2009.17951519752112 Booth FW Roberts CK Laye MJ. Lack of Exercise is a major cause of chronic diseases. Compr Physiol. (2012) 2:1143. 10.1002/cphy.c11002523798298 Handschin C Spiegelman BM. The role of exercise and Pgc1α in inflammation and chronic disease. Nature. (2008) 454:463–9. 10.1038/nature0720618650917 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. 10.1038/nri215617717539 Qin C Zhou L Hu Z Zhang S Yang S Tao Y . Dysregulation of immune response in patients with coronavirus 2019 (Covid-19) in Wuhan, China. Clin Infect Dis. (2020) 71:762–8. 10.1093/cid/ciaa24832161940 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. 10.1016/j.aanat.2015.09.00226434754 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. 10.1016/S0735-1097(98)00093-X9581711 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. 10.5603/FM.a2017.010129064542 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. 10.1249/01.MSS.0000058440.28108.CC12673139 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. 10.1080/02701367.2014.93040725141092 Belda J Ricart S Casan P Giner J Bellido-Casado J Torrejon M . Airway inflammation in the elite athlete and type of sport. Br J Sports Med. (2008) 42:244–8. 10.1136/bjsm.2007.03633517711871 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. 10.1089/jir.2011.003821916608 Muniz Pedrogo DA Chen J Hillmann B Jeraldo P Al-Ghalith G Taneja V . 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. 10.1093/ibd/izy31830321331 Luti S Modesti A Modesti PA. Inflammation, peripheral signals and redox homeostasis in athletes who practice different sports. Antioxidants. (2020) 9:65. 10.3390/antiox911106533143147 Lin CL Hsu YJ Ho HH Chang YC Kuo YW Yeh YT . 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. 10.3390/nu1207197232630786 Zhao X Zhang Z Hu B Huang W Yuan C Zou L. Response of Gut Microbiota to Metabolite Changes Induced by Endurance Exercise. Front Microbiol. (2018) 9:765. 10.3389/fmicb.2018.0076529731746