ACCELEROMETER-BASED SEDENTARY BEHAVIOUR AND PHYSICAL ACTIVITY ARE ASSOCIATED WITH THE GUT MICROBIOTA IN 8507 INDIVIDUALS FROM THE POPULATION-BASED SCAPIS
=========================================================================================================================================================

* Gabriel Baldanzi
* Sergi Sayols-Baixeras
* Elin Ekblom-Bak
* Örjan Ekblom
* Koen F. Dekkers
* Ulf Hammar
* Diem Nguyen
* Shafqat Ahmad
* Ulrika Ericson
* Daniel Arvidsson
* Mats Börjesson
* Peter J. Johanson
* J Gustav Smith
* Göran Bergström
* Lars Lind
* Gunnar Engström
* Johan Ärnlöv
* Beatrice Kennedy
* Marju Orho-Melander
* Tove Fall

## ABSTRACT

**Objective** Population-based studies investigating the relationship between physical activity and the gut microbiota composition have mainly relied on self-reported activity, potentially influenced by reporting bias. Here, we investigated associations of accelerometer-based sedentary behaviour and physical activity with the gut microbiota composition and functional profile in the large Swedish CArdioPulmonary bioImage Study.

**Methods** In 8507 participants aged 50-65, the proportion of time in sedentary (SED), moderate-intensity (MPA), and vigorous-intensity (VPA) physical activity were estimated with hip-worn accelerometer. The gut microbiota was profiled using shotgun metagenomics of fecal samples. We fitted multivariable regression models, and adjusted for sociodemographic, lifestyle, and technical covariates while also accounting for multiple testing.

**Results** Overall, SED and MPA were associated with microbiota species in opposite directions. For example, the strongest positive regression coefficient for MPA and the strongest negative for SED were with *Prevotella copri*, a plant-polysaccharide-degrading bacteria. Species associated with VPA aligned with the MPA associations, although with clear discrepancies. For instance, *Phocaeicola vulgatus* was negatively associated with MPA, while the association with VPA was non-significant and in the positive direction. Additional adjustment for dietary variables or adiposity attenuated some of the associations. For the functional profile, MPA and VPA were generally associated with lower capacity for amino acid degradation.

**Conclusion** Our findings suggest that sedentary behaviour and physical activity are associated with a similar set of gut microbiota species and functions, but in opposite directions. Furthermore, the intensity of physical activity may have specific effects on certain species of the gut microbiota.

## INTRODUCTION

Physical activity has well-established health benefits, including reduced risk for cardiovascular disease, type 2 diabetes,1,2 and psychiatric conditions like depression.3 Conversely, sedentary behaviour, defined as sitting or non-sleep lying activities with low energy expenditure, is associated with increased risk of type 2 diabetes and cardiovascular mortality.4–6 Some studies have though indicated that the risks attributed to sedentary behaviour could be substantially attenuated by physical activity at high intensities.7–9

The gut microbiota is a community of microorganisms within the gastrointestinal tract that interacts with the host.10 Evidence suggests that the gut microbiota plays a role in the development of type 2 diabetes and cardiovascular diseases.11–13 Moreover, the gut microbiota may influence brain homeostasis through the microbiota-gut-brain axis, which includes microbe-produced neurotransmitters and other molecules that can affect the central nervous system via neuronal pathways or the immune system.14

Regular physical activity may affect the gut microbiota through various mechanisms, including modulation of the gut immune system, reduction in the intestinal transit time15 and splanchnic blood flow, transient increase in the intestinal permeability,16 and modulation of the enterohepatic circulation of bile acids.17 Smaller intervention studies in specific populations have reported changes to the gut microbiota composition after structured exercise, with a decrease in *Clostridium* and *Blautia* and an increase in *Bifidobacterium* and *Dorea*.18–20 In population-based studies, self-reported moderate and vigorous-intensity physical activity have been associated with higher gut microbiota diversity21, and sedentary behaviour with increased abundance of *Roseburia hominis* and *Erysipelatoclostridium* species.22 However, these studies had limited microbiota taxonomic resolution and used self-reported physical activity information, which may be affected by reporting bias.23 Therefore, there is a need for larger population-based studies that combine sensor-based physical activity assessment with gut microbiota profiled in higher taxonomic resolution. Here, we aimed to identify associations of accelerometer-based sedentary behaviour and physical activity with the gut microbiota analysed with deep shotgun metagenomics, in 8507 participants, using cross-sectional data from the Swedish CArdioPulmonary BioImage Study (SCAPIS).

## METHODS

### Study Population

The SCAPIS cohort includes 30 154 women and men aged 50 to 65 enrolled between 2013 and 2018 from six different regions in Sweden.24 Baseline investigation included wear of a hip-accelerometer for 7 days. Participants from Malmö and Uppsala were also invited to provide faecal samples for metagenomic analysis (n = 9831). The Swedish Ethical Review Authority approved the SCAPIS study (DNR 2010-228- 31M) and the present study (DNR 2018-315 with amendment 2020-06597). All participants provided written informed consent.

### Accelerometer data processing

Accelerometer data was processed using a previously described protocol25 with the added exclusion of registrations during estimated bedtime. For details, see supplemental methods. In brief, participants were instructed to wear the ActiGraph (Pensacola, USA) tri-axial accelerometer over the hip for 7 consecutive days, except during sleep and water-based activities. Raw data was transformed into counts per minute (cpm) over 60s epochs. Sedentary time (SED), low-intensity physical activity (LIPA), moderate-intensity physical activity (MPA), and vigorous-intensity physical activity (VPA) were defined as <200 cpm, 200 – 2689 cpm, 2690 – 6166 cpm, and ≥6167 cpm, respectively. Non-wear time was defined as periods of ≥60 min with zero counts, with intervals of maximum two minutes of 0 – 199 cpm. A valid day was defined as a day with >10 hours of wear time. We excluded 467 participants who had <4 valid days. The percentage of time in SED, MPA, or VPA were calculated by dividing the time spent in each activity by the total wear time.25 Because SED and LIPA were highly negatively correlated (Spearman’s correlation = −0.95, supplemental figure 1), we chose to exclude LIPA from the subsequent analyses.

### Faecal metagenomics

The gut microbiota was assessed through metagenomic analyses of faecal samples using a previously described protocol.26 In summary, faecal samples were collected at home, kept in the home freezer until the second study visit, and stored at −80°C until shipped to Clinical Microbiomics A/S (Copenhagen, Denmark) for DNA extraction, library preparation, sequencing with Ilumina Novaseq 6000 (Illumina, CA, USA), bioinformatics processing, and taxonomic annotation. Metagenomic species were defined by binning of co-abundance genes, as previously described.27 Alpha diversity, a metric of the diversity of species and the homogeneity of their abundance within a sample28 was estimated using the Shannon diversity index. Beta diversity, a metric of the composition dissimilarity between samples, was estimated using the Bray-Curtis dissimilarity.29 These estimations were performed with R package *vegan* on sequence data that was rarefied to 210 430 read-pairs for all samples to account for the differences in sequencing depth. The functional potential of the gut microbiota was defined by the abundance of genes of the manually curated gut metabolic modules30 (GMM) and microbiota-gut-brain modules14 (MGB). GMMs capture the microbiota metabolic potential and anaerobic fermentation capacity30, while MGBs comprise the capacity to degrade or produce potentially neuroactive compounds.14 To calculate the abundance of respective modules, we used the R package Omixer-RPM v.0.3.231 considering a minimum module coverage of 66.6%. Analyses were focused on species identity, GMMs, and MGBs with a relative abundance >0.01% in ≥1% of individuals, resulting in 1336 species, 76 GMMs, and 28 MGBs for further analysis. Before statistical analysis, species, GMM, and MGB were log(x+1) transformed, where x denotes the relative abundance. For details, see supplemental methods.

### Covariates

Covariate information was obtained from the SCAPIS questionnaire, and anthropometric measurements and fasting plasma samples collected during study site visit. Mean daily intake of alcohol, fibre, added sugar, protein, carbohydrate, fat, and total energy were estimated from the food frequency questionnaire.32,33 Fibre, protein, carbohydrate, fat, and added sugar intake were transformed to percentages of non-alcohol energy intake. We categorized smoking status as current, former, or non-smoker, and highest achieved education level as incomplete compulsory, complete compulsory, secondary, or university education. Country of birth was grouped as Scandinavia (Sweden, Denmark, Norway, or Finland), non-Scandinavian Europe, Asia, and other countries. Information on medications prescribed for hypertension, type 2 diabetes, dyslipidaemia, depression, and anxiety within six months before the first site visit was retrieved from the Swedish Prescribed Drug Register (supplemental methods). Proton-pump inhibitor usage was defined as a measurable level of omeprazole or pantoprazole metabolites in plasma. Antibiotic use was defined as a dispensed prescription (Anatomical Therapeutic Chemical code J01) up to three months before the first visit.

### Statistical method

#### Model specification

Based on current literature21,25,34, we created a directed acyclic graph prior to the analysis phase using the DAGitty tool ([www.daggity.net](http://www.daggity.net), supplemental figure 2), and selected potential confounders for the main model based on the d-separation criteria.35 The main model comprised age, sex, alcohol intake, smoking, education, country of birth, study site, and month of accelerometer wear. We also adjusted for technical variation covariates including total wear time, percentage of wear time on weekend, and faecal DNA extraction plate. To identify associations independent of dietary intake or adiposity, we explored two additional models: (1) Diet-adjusted model: additional adjustment for total energy intake and intake from added sugars, proteins, carbohydrates, and fibre, or (2) BMI/WHR-adjusted model: additional adjustment for body mass index (BMI) and waist-hip ratio (WHR). Statistical analyses were conducted in R version 4.1.1 ([http://www.r-project.org](http://www.r-project.org)).

#### Alpha and beta diversity

For alpha diversity, we applied linear regression models with Shannon index as the dependent variable and SED, MPA, and VPA separately as independent variables, while adjusting for covariates. Regression coefficients presented represent standard-deviation changes in SED, MPA, or VPA. For beta diversity, we used distance-based redundancy analysis to estimate the proportion of the interindividual gut microbiota dissimilarity explained by SED, MPA, and VPA. The analyses accounted for the covariates by using the parameter “condition” in the function “capscale” (R package *vegan*) and p-values were calculated based on 9999 permutations.

#### Species, gut metabolic modules, and microbiota-gut-brain modules

To identify microbiota features associated with physical activity overall, we applied a series of linear regression models with each species, GMM, and MGB introduced as the dependent variable, and SED, MPA, and VPA jointly introduced as independent variables together with the main model covariates. We used the likelihood ratio test to compare these models to an alternative model that only contained the covariates. To control for multiple testing, we applied the Benjamini-Hochberg method with a 5% false discovery rate and reported significance as q-values.36 The species, GMMs, and MGBs identified (q-value <0.05) where subsequently evaluated in models where SED, MPA, and VPA were introduced separately as an independent variable while adjusting for main model covariates, as well as for the diet-adjusted model and the BMI/WHR-adjusted model. To identify associations potentially caused by single influential observations, we calculated dfbetas for the main model. If removing the observation with the highest dfbeta resulted in a change in the direction of the association or a p-value ≥0.05, we discarded the association. Lastly, remaining species were examined in two sensitivity analyses: (1) additional adjustment for use of proton-pump inhibitors and medication for hypertension, type 2 diabetes, dyslipidaemia, anxiety, and/or depression; (2) exclusion of participants that were prescribed antibiotic treatment within the last three months.

## RESULTS

### Study population

The study population consisted of 4197 participants from Malmö and 4310 participants from Uppsala with valid accelerometer data, faecal metagenomics data, and complete information on the main model covariates (Table 1).

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/06/04/2023.06.01.23290817/T1)

Table 1. 
Characteristics of the study population stratified by study site.

### Sedentary behaviour and physical activity associations with alpha diversity and beta diversity

We found that SED was associated with lower alpha diversity (β = −0.020, p-value = 1.2×10-5), whereas MPA and VPA were associated with higher alpha diversity (MPA: β =0.024, p-value = 5.5×10-8; VPA: β = 0.032, p-value = 1.7×10-12). Similar results were found in the diet-adjusted model. In the BMI/WHR-adjusted model, only MPA and VPA were associated with alpha diversity (MPA: β = 0.017, p-value = 1.3×10-4; VPA: β = 0.022, p-value = 1.4×10-6).

To evaluate the contribution of sedentary behaviour and physical activity on the interindividual variation of the gut microbiota composition, we performed distance-based redundancy analysis on the Bray-Curtis dissimilarly matrix. After conditioning on the main model covariates, the adjusted R2 for SED was 0.083%, for MPA 0.084%, and for VPA 0.081%. The adjusted R2 for SED, MPA and VPA jointly was 0.188% in the main model, 0.141% in the diet-adjusted model, and 0.154% in the BMI/WHR-adjusted model (all p-values = 1×10-4). For comparison, the R2 for BMI without accounting for other variables was 0.76%, for fibre intake 0.39%, and for the dietary variables combined 0.69%.

### Sedentary behaviour and physical activity were associated with a large number of the species present in the human gut

We identified 345 species associated with physical activity (supplemental table 1). Among these, SED was associated with 224 species (q-value <0.05), MPA with 245 species, and VPA with 228 species. However, we found that six associations for SED, 16 associations for MPA, and 31 associations for VPA were caused by single influential observations (supplemental table 2). Overall, associations for SED and MPA mirrored each other (figure 1). The Pearson correlation between the main model regression coefficients for SED and MPA was −0.92, while between SED and VPA was −0.67. Amongst the largest coefficients, SED was positively, and MPA and VPA were negatively associated with *Blautia obeum* (unique identifier HG3A.0001) and *Ruminococcus torques*. Additionally, SED was negatively and MPA was positively associated with *Prevotella copri* and an unclassified *Eubacteriales* species (identifier HG3A.0100). The largest positive coefficients for VPA were with four unclassified *Eubacteriales* species (HG3A.0125, HG3A.0100, HG3A.0062, HG3A.0162).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/04/2023.06.01.23290817/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/06/04/2023.06.01.23290817/F1)

Figure 1. 
Top 20 associations based on the absolute regression coefficients for SED, MPA, and VPA with gut microbiota species. Effect estimates show changes in the log(relative abundance+1) of species by standard-deviation changes in SED, MPA, or VPA. Main model: adjustment for age, sex, alcohol intake, smoking, education, country of birth, study site, month of accelerometer wear, total accelerometer wear time, percentage of wear time on weekend, and fecal DNA extraction plate. Diet model: main model with additional adjustment for total energy intake, and percentage of energy intake from carbohydrates, protein, fibers, and added sugars. BMI/WHR model: main model with additional adjustment for BMI and waist-hip ratio. Circles are the regression coefficients and bars represent the 95% confidence intervals. Filled circles are associations with q-values <0.05. SED: percentage of time in sedentary behaviour; MPA: percentage of time in moderate-intensity physical activity; VPA: percentage of time in vigorous-intensity physical activity.

Regression coefficients for MPA and VPA were in largely consistent in direction and magnitude (Pearson correlation = 0.67). However, for specific species, there were marked differences in the association with MPA or VPA. For instance, *P. copri* and *Faecalibacterium prausnitzii* (HG3A.0025) were among the largest positive coefficients for MPA, but were not associated with VPA. Other differences include *Phocaeicola vulgatus* and *Bacteroides uniformis*, which were negatively associated with MPA, but not with VPA (figure 1).

The main model and the diet-adjusted model regression coefficients were highly correlated (Pearson correlation >0.99 for SED, MPA, and VPA), although generally attenuated in diet-adjusted model (supplemental table 2 and supplemental figure 3). The BMI/WHR-adjusted model also produced similar coefficients for MPA and VPA compared to the main model coefficients (Pearson correlation >0.97). However, regression coefficients for SED from the BMI/WHR model and from the main model were not as strongly correlated (Pearson correlation = 0.89). Certain main model associations for SED were substantially attenuated in the BMI/WHI model, especially for *B. obeum* (HG3A.0001) which was no longer associated with SED (figure 1 and supplemental table 2). On the other hand, the negative association between SED and *F. prausnitzii* (HG3A.0025) became stronger in the BMI/WHR model.

To investigate the potential effect of medication usage on the associations identified in the main model (q-value <0.05), we performed further adjustment for use of proton-pump inhibitors and medications for hypertension, type 2 diabetes, dyslipidaemia, anxiety, and/or depression. After this adjustment, SED was associated (p-value < 0.05) with 187 of the 224 species associated in the main model, MPA with 205 of 245, and VPA with 169 of 228 species. In the sensitivity analysis removing 498 participants who had used any antibiotic in the last 3 months, SED was associated (p-value < 0.05) with 215 species, MPA with 221, and VPA with 195 (supplemental table 3). The direction of the associations did not change.

### Sedentary behaviour and physical activity are associated with the gut microbiota functional potential

The functional potential was determined by the abundance of previously curated functional modules, which include modules of carbohydrate, amino acid, or lipid degradation,30 and modules of the microbiota-gut-brain axis.14 We identified 74 functional modules associated with physical activity (supplemental table 4). Out of 17 modules of carbohydrate degradation, SED was associated with lower or higher abundance of 12 modules in the main model (figure 2 and supplemental table 5). The strongest negative association was between SED and degradation of arabinoxylan, which is a dietary fibre. Even in the diet-adjusted model, SED was associated with degradation of arabinoxylan. Out of the 26 modules of amino acid degradation, SED was associated with higher abundance of 17 in the main model. In the diet-adjusted model, SED was associated with higher abundance of 15 modules of amino acid degradation, whereas MPA was associated with lower abundance of 14, and VPA with lower abundance of 7.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/04/2023.06.01.23290817/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/06/04/2023.06.01.23290817/F2)

Figure 2. 
Associations of SED, MPA and VPA with functional modules of carbohydrate, amino acid, and lipid degradation. Effect estimates show changes in the log(relative abundance+1) of modules by standard-deviation change in SED, MPA, or VPA. Main model: adjustment for age, sex, alcohol intake, smoking, education, country of birth, study site, month of accelerometer wear, total accelerometer wear time, percentage of wear time on weekend, and fecal DNA extraction plate. Diet model: main model with additional adjustment for total energy intake, and percentage of energy intake from carbohydrates, protein, fibers, and added sugars. BMI/WHR model: main model with additional adjustment for BMI and waist-hip ratio. Circles are the regression coefficients and bars represent the 95% confidence intervals. Filled circles are associations with q-values <0.05. SED: percentage of time in sedentary behaviour; MPA: percentage of time in moderate-intensity physical activity; VPA: percentage of time in vigorous-intensity physical activity.

Among the modules of the microbiota-gut-brain axis, the largest main model coefficient was found for heat-shock protein ClpB, which was negatively associated with SED, and positively associated with MPA and VPA (figure 3 and supplemental table 5). Comparing main model and BMI/WHR-adjusted model associations, we observed marked differences for the butyrate synthesis I module. In the main model, VPA was positively associated with butyrate synthesis I, while no association was detected for SED or MPA. In the BMI/WHR-adjusted model, we could not detect the association for VPA; instead, SED was positively and MPA was negatively associated with this module. No association was driven by a single influential observation.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/04/2023.06.01.23290817/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/06/04/2023.06.01.23290817/F3)

Figure 3. 
Associations of SED, MPA and VPA with microbiota-gut-brain modules and other functional modules. Effect estimates show changes in the log(relative abundance+1) of modules by standard-deviation change in SED, MPA, or VPA. Main model: adjustment for age, sex, alcohol intake, smoking, education, country of birth, study site, month of accelerometer wear, total accelerometer wear time, percentage of wear time on weekend, and fecal DNA extraction plate. Diet model: main model with additional adjustment for total energy intake, and percentage of energy intake from carbohydrates, protein, fibers, and added sugars. BMI/WHR model: main model with additional adjustment for BMI and waist-hip ratio. Circles are the regression coefficients and bars represent the 95% confidence intervals. Filled circles are associations with q-values <0.05. SED: percentage of time in sedentary behaviour; MPA: percentage of time in moderate-intensity physical activity; VPA: percentage of time in vigorous-intensity physical activity.

## DISCUSSION

In this largest-to-date population-based study of physical activity with gut microbiome encompassing 8507 individuals, we found that, from the 1336 species investigated, 345 (25.8%) were associated with physical activity. Overall, SED and MPA were associated with the same set of species but with regression coefficients in opposite directions. Furthermore, MPA and VPA had concordant associations with gut microbiota species, although with notable exceptions. Similar results were observed after adjustment for dietary variables, but some associations for SED were markedly affected by the adjustment for BMI/WHR. Additionally, SED, MPA, and VPA were associated with specific microbial functions, some potentially involved in the microbiota-gut-brain axis.

### The associations of sedentary behaviour and physical activity with gut microbiota species

Among the strongest findings with species-level annotation, MPA was associated with higher abundance of *P. copri* and *F. prausnitzii* (HG3A.0241 and HG3A.0025). Time spent in exercise has previously been associated with higher *P. copri* abundance in cyclists.37 Nevertheless, the associations between *P. copri* and health outcomes have been inconsistent, with some studies suggesting a role in fibre fermentation and improved glucose metabolism,38 while others reported associations with insulin resistance39 and hypertension.40 MPA and VPA had similar associations, but there were some clear exceptions. MPA was associated with higher abundance of *P. copri* and *F. prausnitzii* (HG3A.0025), and lower abundance of *B. uniformis* and *P. vulgatus*, while we could not detect the same associations for VPA. We found that VPA was positively associated with *Roseburia homininis* and *F. prausnitzii* (HG3A.0241), but negatively associated with another subspecies of *F. prausnitzii* (HG3A.0010). In general, *F. prausnitzii* is considered to have anti-inflammatory properties41 and individuals with type 2 diabetes have been reported to have a lower abundance of this bacteria.42 A higher abundance of *F. prausnitzii* and *R. hominis* has also been reported in active premenopausal women, based on accelerometer assessment.43 Both *F. prausnitzii* and *R. hominis* are important producers of butyrate, which is a main energy source for colonocytes and critical for gut homeostasis.44 A previous study in a small sample of young adults has also reported an association between cardiorespiratory fitness and fecal butyric acid.45

In the present study, SED was associated with higher abundance of *B. obeum* (HG3A.0001 and HG3A.0009), *R. torques*, *P. vulgatus,* and *Escherichia coli*. An increased abundance of the *Escherichia/Shigella* taxon was also reported in sedentary individuals in the Healthy Life in an Urban Setting cohort (HELIUS, N = 1334), the largest population-based study on physical activity and gut microbiota until the present study.22 A higher abundance of *E. coli* has been observed in individuals with atherosclerotic cardiovascular disease46 and users of the anti-diabetes medication metformin,47 suggesting a link between physical activity, gut microbiota, and cardiometabolic diseases. Due to different microbiota profiling methods, our findings cannot be directly compared with the HELIUS study. We confirmed though the positive association between *P. copri* and physical activity but found a positive association between VPA and *R. hominis*, which was more abundant in sedentary individuals in the HELIUS study.

The median VPA in our study sample was only 0.06%. A previous study has suggested that as little as 15 min/week of VPA reduces the risk for all-cause mortality by 18%;48 therefore, effects of VPA could be expected even with little time spent in these activities. Vigorous-intensity activities include running and high-intensity sports, while moderate-intensity activities cover a wider range of activities, including brisk walking and bicycling at low speed, but also household chores.49 Future studies could further investigate the associations of subcategories of MPA with the gut microbiota.

### Adjustment for BMI and WHR had a pronounced effect on the associations observed between SED and species of the gut microbiota

In the models additionally adjusted for dietary variables or for BMI/WHR, we observed attenuated coefficients. The largest effect was observed for associations with SED after adjustment for BMI/WHR. This attenuation could be caused by BMI acting as a confounder or a mediator. However, for 14 species, including *F. prausnitzii* (HG3A.0025), the association with SED became stronger after BMI/WHR adjustment. This could be due to a negative confounding effect of adiposity, which masked the true associations in the main model, or that the adjustment for BMI/WHR led to collider bias. The greater impact of adjusting for BMI/WHR on SED coefficients than on MPA or VPA coefficients aligns with previous Mendelian randomization analyses, which have demonstrated stronger bidirectional effects between SED and BMI than between physical activity time and BMI.50,51

### Sedentary behaviour associated with higher gut microbiota capacity to degrade amino acids and to produce certain short-chain fatty acids

With regards to the functional potential of the gut microbiota, we found that SED was associated with lower capacity for fibre degradation and higher capacity for amino acid degradation. Likewise, a previous study has reported an increased abundance of carbohydrate degradation pathways in athletes.52 These findings could be due to lower fibre content in the diet of sedentary individuals, which we aimed to address by adjusting for fibre intake. A lower availability of fermentable carbohydrates in the distal gut can result in a reduction in saccharolytic bacteria and an increase in proteolytic bacteria.30 Intervention studies with standardized diet would be needed to disentangle physical activity associations from associations due to differences dietary intakes.

In our main model results, MPA was negatively associated with modules for short-chain fatty acids synthesis, more specifically acetate and propionate synthesis, while SED was positively associated with these modules. Conversely, six weeks of exercise training has been reported to increase the abundance of bacterial genes involved in propionate synthesis.53 The same study also reported a BMI-dependent exercise-induced increase in bacterial genes of butyrate synthesis. In our study, VPA was positively associated with butyrate synthesis in the main model, while in the BMI/WHR-adjusted model, it was SED that was positively associated with this same module. For other modules in the microbiota-gut-brain axis, SED was negatively and MPA positively associated with the module for the heat-shock protein ClpB that has been suggested to influence appetite regulation.54 SED was also associated with higher abundance of the GABA synthesis module and higher abundance of *E. coli*, one of the main GABA-producing bacteria in the gut.55 Although lower plasma levels of GABA have been described in individuals with depression,56 recent studies found elevated levels in a mixed sample of medicated and non-medicated individuals with major depressive disorder.57,58

### Strengths and limitations

The strengths of this study are the accelerometer-based assessment of sedentary behaviour and physical activity phenotypes, the large sample of participants from the general population, and the high taxonomic resolution microbiome data. Moreover, we had access to comprehensive information on potential confounders. Some limitations apply. One concern is whether the associations described reflect the lifestyle of health-conscious individuals. Despite covariate adjustments, it is implausible to capture all dimensions of dietary intake and residual confounding may remain. Accelerometers measure absolute physical activity intensity, but relative intensity could be more clinically relevant. Standardized accelerometer cut-offs can misclassify low and high fitness individuals. Estimating relative intensity would though require data on individual maximal capacity, such as maximal oxygen uptake.59 Additionally, social desirability bias and adherence to the study instructions could affect accelerometer-based assessment. However, these misclassifications would be non-differential. The accelerometer needed to be removed during water-based activity and may also underestimate physical activity intensity during cycling, upper body activities, and weight-lifting.25 Our study was conducted in a Swedish population aged 50-65, thus generalizability to other populations is limited. Additionally, we cannot assess the direction of the associations using cross-sectional data. It is suggested that gut microbiota may favour the practice of exercise by enhancing the enjoyment of physical activity60 or improving the host performance.61 The gut microbiota can also affect adiposity,62 which is negatively associated with physical activity. Lastly, the use of relative abundances incur on the issues of compositional data, which can produce false-positive associations.63

In summary, sedentary behaviour and physical activity were associated with a large number of gut microbiota species and functional modules. Our findings can be used to guide research on the interplay between physical activity, the gut microbiota composition, and health outcomes.

## Supporting information

Supplemental tables 1-5 [[supplements/290817_file02.xlsx]](pending:yes)

Supplemental figures 1-3 [[supplements/290817_file03.pdf]](pending:yes)

Supplemental methods [[supplements/290817_file04.pdf]](pending:yes)

## Data Availability

The de-identified metagenomic sequences can be found in the European Nucleotide Archive under the accession code PRJEB51353. The statistical analyses R codes can be found at [https://github.com/MolEpicUU/physicalactivity\_gut](https://github.com/MolEpicUU/physicalactivity_gut). The individual-level data underlying this article were provided by the SCAPIS study and are not shared publicly due to confidentiality. Data will be shared upon reasonable request to the corresponding author only after permission by the Swedish Ethical Review Authority ([https://etikprovningsmyndigheten.se](https://etikprovningsmyndigheten.se)) and by the SCAPIS Data access board ([https://www.scapis.org/data-access/](https://www.scapis.org/data-access/)).

[https://github.com/MolEpicUU/physicalactivity\_gut](https://github.com/MolEpicUU/physicalactivity_gut)

## Contributorship

JGS, GBergström, LL, GE, JÄ, MO-M and TF obtained the funding for the study. GBaldanzi, SS-B, EEB, ÖK, UH, BK, and TF planned and designed the study. GBaldanzi carried out the statistical analyses with contribution from SS-B, KFD, and UH. GBaldanzi wrote the first version of the manuscript with support from SS-B, BK, and TF. All authors contributed with the critical interpretation of the results and revision of the manuscript.

## Funding/support

Financial support was obtained in the form of grants from the European Research Council [ERC-STG- 2018-801965 (TF); ERC-CoG-2014-649021 (MO-M); ERC-STG-2015-679242 (JGS)], the Swedish Heart-Lung Foundation [Hjärt-Lungfonden, 2019-0505 (TF); 2018-0343 (JÄ); 2020-0711 (MO-M); 2020-0173 (GE); 2019-0526 (JGS)], the Swedish Research Council [VR, 2019-01471 (TF), 2018-02784 (MO-M), 2019-01015 (JÄ), 2020-00243 (JÄ), 2019-01236 (GE), 2021-02273 (JGS), 2018-02837 (EXODIAB 2009-1039)], the Swedish Research Council for Sustainable Development [FORMAS, 2020- 00989 (SA)], EASD/Novo Nordisk (SA), Göran Gustafsson foundation [2016 (TF)], Axel and Signe Lagerman’s foundation (TF), the A.L.F. governmental grant 2018-0148 (MO-M), The Novo Nordic Foundation NNF20OC0063886 (MO-M), The Swedish Diabetes Foundation DIA 2018-375 (MO-M), the Skandia Risk&Hälsa, 2023 (ÖE, EEB), and governmental funding of clinical research within the Swedish National Health Service (JGS). The study was also supported by the Swedish Foundation for Strategic Research (LUDC-IRC 15-0067). Funding for the SCAPIS study was provided by the Swedish Heart-Lung Foundation, the Knut and Alice Wallenberg Foundation, the Swedish Research Council and VINNOVA (Sweden’s Innovation agency), the University of Gothenburg and Sahlgrenska University Hospital, Karolinska Institutet and Stockholm County council, Linköping University and University Hospital, Lund University and Skåne University Hospital, Umeå University and University Hospital, Uppsala University and University Hospital.

## Competing interests

J.Ä. has served on the advisory boards for Astella, AstraZeneca and Boehringer Ingelheim and has received lecturing fees from AstraZeneca and Novartis, all unrelated to the present work. Remaining authors declare no competing interests.

## Data Sharing

The de-identified metagenomic sequences can be found in the European Nucleotide Archive under the accession code “PRJEB51353”. The statistical analyses R codes can be found at [https://github.com/MolEpicUU/physicalactivity\_gut](https://github.com/MolEpicUU/physicalactivity_gut). The individual-level data underlying this article were provided by the SCAPIS study and are not shared publicly due to confidentiality. Data will be shared upon reasonable request to the corresponding author only after permission by the Swedish Ethical Review Authority ([https://etikprovningsmyndigheten.se](https://etikprovningsmyndigheten.se)) and by the SCAPIS Data access board ([https://www.scapis.org/data-access/](https://www.scapis.org/data-access/)).

## Acknowledgments

We acknowledge the SCAPIS board for enabling the current study, along with the Swedish Heart-Lung Foundation, the main funding body of SCAPIS. The computations and data handling were enabled by resources in project sens2019512 provided by the National Academic Infrastructure for Supercomputing in Sweden (NAISS) and the Swedish National Infrastructure for Computing (SNIC) at Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX), funded by the Swedish Research Council through grant agreement no. 2018-0597 and no. 2022-06725.

*   Received June 1, 2023.
*   Revision received June 1, 2023.
*   Accepted June 4, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## REFERENCES

1.  1.Aune, D., Norat, T., Leitzmann, M., Tonstad, S. & Vatten, L. J. Physical activity and the risk of type 2 diabetes: a systematic review and dose-response meta-analysis. Eur. J. Epidemiol. 30, 529–542 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10654-015-0056-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26092138&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

2.  2.Lear, S. A. et al. The effect of physical activity on mortality and cardiovascular disease in 130 000 people from 17 high-income, middle-income, and low-income countries: the PURE study. Lancet Lond. Engl. 390, 2643–2654 (2017).
    
    

3.  3.Brown, W. J., Ford, J. H., Burton, N. W., Marshall, A. L. & Dobson, A. J. Prospective Study of Physical Activity and Depressive Symptoms in Middle-Aged Women. Am. J. Prev. Med. 29, 265–272 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amepre.2005.06.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16242588&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000232809500003&link_type=ISI) 

4.  4.Patterson, R. et al. Sedentary behaviour and risk of all-cause, cardiovascular and cancer mortality, and incident type 2 diabetes: a systematic review and dose response meta-analysis. Eur. J. Epidemiol. 33, 811–829 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10654-018-0380-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29589226&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

5.  5.Wilmot, E. G. et al. Sedentary time in adults and the association with diabetes, cardiovascular disease and death: systematic review and meta-analysis. Diabetologia 55, 2895–2905 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-012-2677-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22890825&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309561600004&link_type=ISI) 

6.  6.Teychenne, M., Costigan, S. A. & Parker, K. The association between sedentary behaviour and risk of anxiety: a systematic review. Bmc Public Health 15, 513 (2015).
    
    

7.  7.Ekelund, U. et al. Do the associations of sedentary behaviour with cardiovascular disease mortality and cancer mortality differ by physical activity level? A systematic review and harmonised meta-analysis of data from 850 060 participants. Br. J. Sports Med. 53, 886–894 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiYmpzcG9ydHMiO3M6NToicmVzaWQiO3M6OToiNTMvMTQvODg2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDQvMjAyMy4wNi4wMS4yMzI5MDgxNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

8.  8.Ekelund, U. et al. Does physical activity attenuate, or even eliminate, the detrimental association of sitting time with mortality? A harmonised meta-analysis of data from more than 1 million men and women. The Lancet 388, 1302–1310 (2016).
    
    

9.  9.Ekelund, U. et al. Joint associations of accelero-meter measured physical activity and sedentary time with all-cause mortality: a harmonised meta-analysis in more than 44 000 middle-aged and older individuals. Br. J. Sports Med. 54, 1499–1506 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiYmpzcG9ydHMiO3M6NToicmVzaWQiO3M6MTA6IjU0LzI0LzE0OTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNC8yMDIzLjA2LjAxLjIzMjkwODE3LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

10. 10.Qin, J. et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464, 59–65 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature08821&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20203603&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275117500033&link_type=ISI) 

11. 11.Molinaro, A. et al. Imidazole propionate is increased in diabetes and associated with dietary patterns and altered microbial ecology. Nat. Commun. 11, 5881 (2020).
    
    

12. 12.Witkowski, M., Weeks, T. L. & Hazen, S. L. Gut Microbiota and Cardiovascular Disease. Circ. Res. 127, 553–570 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/CIRCRESAHA.120.316242&link_type=DOI) 

13. 13.Pinto, Y. et al. Gestational diabetes is driven by microbiota-induced inflammation months before diagnosis. Gut gutjnl-2022-328406 (2023) doi:10.1136/gutjnl-2022-328406.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjg6IjcyLzUvOTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDQvMjAyMy4wNi4wMS4yMzI5MDgxNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

14. 14.Valles-Colomer, M. et al. The neuroactive potential of the human gut microbiota in quality of life and depression. Nat. Microbiol. 4, 623–632 (2019).
    
    

15. 15.Stanich, P. P., Peck, J., Murphy, C., Porter, K. M. & Meyer, M. M. Physical activity during video capsule endoscopy correlates with shorter bowel transit time. Endosc. Int. Open 5, E856–E860 (2017).
    
    

16. 16.van Wijck, K. et al. Exercise-induced splanchnic hypoperfusion results in gut dysfunction in healthy men. PloS One 6, e22366 (2011).
    
    

17. 17.Wertheim, B. C. et al. Physical activity as a determinant of fecal bile acid levels. Cancer Epidemiol. Biomark. Prev. Publ. Am. Assoc. Cancer Res. Cosponsored Am. Soc. Prev. Oncol. 18, 1591–1598 (2009).
    
    

18. 18. Donati Zeppa, S., et al. Nine weeks of high-intensity indoor cycling training induced changes in the microbiota composition in non-athlete healthy male college students. J. Int. Soc. Sports Nutr. 18, 74 (2021).
    
    

19. 19.Kern, T. et al. Structured exercise alters the gut microbiota in humans with overweight and obesity— A randomized controlled trial. Int. J. Obes. 44, 125–135 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41366-019-0440-y&link_type=DOI) 

20. 20.Motiani, K. K. et al. Exercise Training Modulates Gut Microbiota Profile and Improves Endotoxemia. Med. Sci. Sports Exerc. 52, 94–104 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/MSS.0000000000002112&link_type=DOI) 

21. 21.Manor, O. et al. Health and disease markers correlate with gut microbiome composition across thousands of people. Nat. Commun. 11, 5206 (2020).
    
    

22. 22.Houttu, V. et al. Physical Activity and Dietary Composition Relate to Differences in Gut Microbial Patterns in a Multi-Ethnic Cohort-The HELIUS Study. Metabolites 11, 858 (2021).
    
    

23. 23.Adams, S. A. et al. The effect of social desirability and social approval on self-reports of physical activity. Am. J. Epidemiol. 161, 389–398 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwi054&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15692083&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000227145500011&link_type=ISI) 

24. 24.Bergström, G. et al. Prevalence of Subclinical Coronary Artery Atherosclerosis in the General Population. Circulation 144, 916–929 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/circulationaha.121.055340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

25. 25.Ekblom-Bak, E. et al. Accelerometer derived physical activity patterns in 27.890 middle-aged adults: The SCAPIS cohort study. Scand. J. Med. Sci. Sports 32, 866–880 (2022).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

26. 26.Dekkers, K. F. et al. An online atlas of human plasma metabolite signatures of gut microbiome composition. Nat. Commun. 13, 5370 (2022).
    
    

27. 27.Nielsen, H. B. et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat. Biotechnol. 32, 822–828 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.2939&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24997787&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

28. 28.Spellerberg, I. F. & Fedor, P. J. A tribute to Claude Shannon (1916–2001) and a plea for more rigorous use of species richness, species diversity and the ‘Shannon–Wiener’ Index. Glob. Ecol. Biogeogr. 12, 177–179 (2003).
    
    

29. 29.Bray, J. R. & Curtis, J. T. An Ordination of the Upland Forest Communities of Southern Wisconsin. Ecol. Monogr. 27, 325–349 (1957).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/1942268&link_type=DOI) 

30. 30.Vieira-Silva, S. et al. Species-function relationships shape ecological properties of the human gut microbiome. Nat. Microbiol. 1, 16088 (2016).
    
    

31. 31.Darzi, Y., Falony, G., Vieira-Silva, S. & Raes, J. Towards biome-specific analysis of meta-omics data. ISME J. 10, 1025–1028 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ismej.2015.188&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26623543&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

32. 32.Christensen, S. E. et al. Relative Validity of Micronutrient and Fiber Intake Assessed With Two New Interactive Meal-and Web-Based Food Frequency Questionnaires. J. Med. Internet Res. 16, e59 (2014).
    
    

33. 33.Christensen, S. E. et al. Two New Meal-and Web-Based Interactive Food Frequency Questionnaires: Validation of Energy and Macronutrient Intake. J. Med. Internet Res. 15, e109 (2013).
    
    

34. 34.Aasmets, O., Krigul, K. L., Lüll, K., Metspalu, A. & Org, E. Gut metagenome associations with extensive digital health data in a volunteer-based Estonian microbiome cohort. Nat. Commun. 13, 869 (2022).
    
    

35. 35.Pearl, J. Causal Diagrams for Empirical Research. Biometrika 82, 669–688 (1995).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomet/82.4.669&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995TT22800001&link_type=ISI) 

36. 36.Benjamini, Y. & Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 57, 289–300 (1995).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.2517-6161.1995.tb02031.x)&link_type=DOI) 

37. 37.Petersen, L. M. et al. Community characteristics of the gut microbiomes of competitive cyclists. Microbiome 5, 98 (2017).
    
    

38. 38.Kovatcheva-Datchary, P. et al. Dietary Fiber-Induced Improvement in Glucose Metabolism Is Associated with Increased Abundance of Prevotella. Cell Metab. 22, 971–982 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2015.10.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26552345&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

39. 39.Pedersen, H. K. et al. Human gut microbes impact host serum metabolome and insulin sensitivity. Nature 535, 376–381 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature18646&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27409811&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

40. 40.Li, J. et al. Gut microbiota dysbiosis contributes to the development of hypertension. Microbiome 5, 14 (2017).
    
    

41. 41.Breyner, N. M. et al. Microbial Anti-Inflammatory Molecule (MAM) from Faecalibacterium prausnitzii Shows a Protective Effect on DNBS and DSS-Induced Colitis Model in Mice through Inhibition of NF-κB Pathway. Front. Microbiol. 8, 114 (2017).
    
    

42. 42.Qin, J. et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 490, 55–60 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11450&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23023125&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309446800031&link_type=ISI) 

43. 43.Bressa, C. et al. Differences in gut microbiota profile between women with active lifestyle and sedentary women. PLOS ONE 12, e0171352 (2017).
    
    

44. 44.Donohoe, D. R. et al. The Microbiome and Butyrate Regulate Energy Metabolism and Autophagy in the Mammalian Colon. Cell Metab. 13, 517–526 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2011.02.018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21531334&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000290291300006&link_type=ISI) 

45. 45.Estaki, M. et al. Cardiorespiratory fitness as a predictor of intestinal microbial diversity and distinct metagenomic functions. Microbiome 4, 42 (2016).
    
    

46. 46.Jie, Z. et al. The gut microbiome in atherosclerotic cardiovascular disease. Nat. Commun. 8, 845 (2017).
    
    

47. 47.Wu, H. et al. Metformin alters the gut microbiome of individuals with treatment-naive type 2 diabetes, contributing to the therapeutic effects of the drug. Nat. Med. 23, 850–858 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm.4345&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28530702&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

48. 48.Ahmadi, M. N. et al. Vigorous physical activity, incident heart disease, and cancer: how little is enough? Eur. Heart J. 43, 4801–4814 (2022).
    
    

49. 49.Ainsworth, B. E. et al. 2011 Compendium of Physical Activities: A Second Update of Codes and MET Values. Med. Sci. Sports Exerc. 43, 1575–1581 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/MSS.0b013e31821ece12&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21681120&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000292773000025&link_type=ISI) 

50. 50.Carrasquilla, G. D., García-Ureña, M., Fall, T., Sørensen, T. I. A. & Kilpeläinen, T. O. Mendelian randomization suggests a bidirectional, causal relationship between physical inactivity and adiposity. eLife 11, e70386 (2022).
    
    

51. 51.Wang, Z. et al. Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention. Nat. Genet. 54, 1332–1344 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-022-01165-1&link_type=DOI) 

52. 52.Barton, W. et al. The microbiome of professional athletes differs from that of more sedentary subjects in composition and particularly at the functional metabolic level. Gut 67, 625–633 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjg6IjY3LzQvNjI1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDQvMjAyMy4wNi4wMS4yMzI5MDgxNy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

53. 53.Allen, J. M. et al. Exercise Alters Gut Microbiota Composition and Function in Lean and Obese Humans. Med. Sci. Sports Exerc. 50, 747–757 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/MSS.0000000000001495&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

54. 54.Tennoune, N. et al. Bacterial ClpB heat-shock protein, an antigen-mimetic of the anorexigenic peptide α-MSH, at the origin of eating disorders. Transl. Psychiatry 4, e458 (2014).
    
    

55. 55.Strandwitz, P. et al. GABA-modulating bacteria of the human gut microbiota. Nat. Microbiol. 4, 396– 403 (2019).
    
    

56. 56.Petty, F., Kramer, G. L., Gullion, C. M. & Rush, A. J. Low plasma gamma-aminobutyric acid levels in male patients with depression. Biol. Psychiatry 32, 354–363 (1992).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0006-3223(92)90039-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1420649&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JW49300006&link_type=ISI) 

57. 57.Bhandage, A. K. et al. Depression, GABA, and Age Correlate with Plasma Levels of Inflammatory Markers. Int. J. Mol. Sci. 20, 6172 (2019).
    
    

58. 58.Pan, J.-X. et al. Diagnosis of major depressive disorder based on changes in multiple plasma neurotransmitters: a targeted metabolomics study. Transl. Psychiatry 8, 1–10 (2018).
    
    

59. 59.Warner, A., Vanicek, N., Benson, A., Myers, T. & Abt, G. Agreement and relationship between measures of absolute and relative intensity during walking: A systematic review with meta-regression. PloS One 17, e0277031 (2022).
    
    

60. 60.Dohnalová, L. et al. A microbiome-dependent gut-brain pathway regulates motivation for exercise. Nature 612, 739–747 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-022-05525-z&link_type=DOI) 

61. 61.Scheiman, J. et al. Meta-omics analysis of elite athletes identifies a performance-enhancing microbe that functions via lactate metabolism. Nat. Med. 25, 1104–1109 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-019-0485-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 

62. 62.Turnbaugh, P. J. et al. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444, 1027–1031 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature05414&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17183312&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F04%2F2023.06.01.23290817.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000242971100045&link_type=ISI) 

63. 63.Gloor, G. B., Macklaim, J. M., Pawlowsky-Glahn, V. & Egozcue, J. J. Microbiome Datasets Are Compositional: And This Is Not Optional. Front. Microbiol. 8, (2017).