Skip to main content
JAMA Network logoLink to JAMA Network
. 2022 Apr 13;79(6):550–559. doi: 10.1001/jamapsychiatry.2022.0609

Association Between Physical Activity and Risk of Depression

A Systematic Review and Meta-analysis

Matthew Pearce 1, Leandro Garcia 1,2, Ali Abbas 1, Tessa Strain 1, Felipe Barreto Schuch 3, Rajna Golubic 1,4, Paul Kelly 5, Saad Khan 6, Mrudula Utukuri 6, Yvonne Laird 7,8, Alexander Mok 1,9, Andrea Smith 1, Marko Tainio 1,10,11, Søren Brage 1, James Woodcock 1,
PMCID: PMC9008579  PMID: 35416941

This systematic review and meta-analysis examines data from published prospective studies of adults to gauge the dose-response association between physical activity and incident depression.

Key Points

Question

What is the dose-response association between physical activity and incident depression in adults?

Findings

This systematic review and meta-analysis of 15 prospective studies including more than 2 million person-years showed an inverse curvilinear association between physical activity and incident depression, with greater differences in risk at lower exposure levels. Adults meeting physical activity recommendations (equivalent to 2.5 h/wk of brisk walking) had lower risk of depression, compared with adults reporting no physical activity.

Meaning

In this study, relatively small doses of physical activity were associated with substantially lower risks of depression.

Abstract

Importance

Depression is the leading cause of mental health–related disease burden and may be reduced by physical activity, but the dose-response relationship between activity and depression is uncertain.

Objective

To systematically review and meta-analyze the dose-response association between physical activity and incident depression from published prospective studies of adults.

Data Sources

PubMed, SCOPUS, Web of Science, PsycINFO, and the reference lists of systematic reviews retrieved by a systematic search up to December 11, 2020, with no language limits. The date of the search was November 12, 2020.

Study Selection

We included prospective cohort studies reporting physical activity at 3 or more exposure levels and risk estimates for depression with 3000 or more adults and 3 years or longer of follow-up.

Data Extraction and Synthesis

Data extraction was completed independently by 2 extractors and cross-checked for errors. A 2-stage random-effects dose-response meta-analysis was used to synthesize data. Study-specific associations were estimated using generalized least-squares regression and the pooled association was estimated by combining the study-specific coefficients using restricted maximum likelihood.

Main Outcomes and Measures

The outcome of interest was depression, including (1) presence of major depressive disorder indicated by self-report of physician diagnosis, registry data, or diagnostic interviews and (2) elevated depressive symptoms established using validated cutoffs for a depressive screening instrument.

Results

Fifteen studies comprising 191 130 participants and 2 110 588 person-years were included. An inverse curvilinear dose-response association between physical activity and depression was observed, with steeper association gradients at lower activity volumes; heterogeneity was large and significant (I2 = 74%; P < .001). Relative to adults not reporting any activity, those accumulating half the recommended volume of physical activity (4.4 marginal metabolic equivalent task hours per week [mMET-h/wk]) had 18% (95% CI, 13%-23%) lower risk of depression. Adults accumulating the recommended volume of 8.8 mMET hours per week had 25% (95% CI, 18%-32%) lower risk with diminishing potential benefits and higher uncertainty observed beyond that exposure level. There were diminishing additional potential benefits and greater uncertainty at higher volumes of physical activity. Based on an estimate of exposure prevalences among included cohorts, if less active adults had achieved the current physical activity recommendations, 11.5% (95% CI, 7.7%-15.4%) of depression cases could have been prevented.

Conclusions and Relevance

This systematic review and meta-analysis of associations between physical activity and depression suggests significant mental health benefits from being physically active, even at levels below the public health recommendations. Health practitioners should therefore encourage any increase in physical activity to improve mental health.

Introduction

Depression is the leading cause of mental health–related disease burden and a major cause of disability worldwide, affecting approximately 280 million people and accounting for more than 47 million disability-adjusted life-years in 2019.1 Depression is also associated with premature mortality from other illnesses2 and suicide.3

Prevention of depression requires effective interventions, including modification of established risk factors.4 Narrative reviews have concluded that physical activity can prevent future depression.5,6 One meta-analysis of prospective studies reported that compared with people with low levels of physical activity, those with higher levels had 17% (95% CI, 12%-21%) lower odds of developing depression,7 while another meta-analysis reported 21% (95% CI, 18%-25%) lower odds when synthesizing 106 associations from 65 studies using diverse exposure definitions.8 To our knowledge, no study has yet synthesized the evidence to describe the strength or shape of the association by conducting a dose-response meta-analysis using harmonized exposure estimates.

Estimating the dose-response relationship between physical activity and any health outcome using meta-analysis is challenging because of the diversity of assessment and inconsistent reporting. Most often, the harmonization of different estimates of physical activity exposure is achieved using binary categorization of low vs high activity, but this approach results in loss of information and does not tell us about the variation in risk across a range of physical activity doses. In contrast, previous work has shown that by using more detailed exposure harmonization, it is possible to investigate the shape of the association between physical activity and type 2 diabetes using dose-response meta-analytical techniques.9 Examining the association between the dose of physical activity and depression in this way allows approximation of the reduction in risk associated with different levels of activity. By combining estimates of reduction in risk with population prevalence estimates of activity, it is possible to quantify the population burden of depression related to insufficient physical activity and the potential public health impact of activity-related interventions.

The aim of this systematic review and meta-analysis was to investigate the dose-response association between physical activity and depression. We also assessed the potential population changes in depression that may be preventable by higher physical activity levels.

Methods

The study protocol is available at PROSPERO (CRD42018095507) and followed Meta-analysis of Observational Studies in Epidemiology (MOOSE) reporting guidelines.

Eligibility Criteria

We included prospective cohort studies of adults (≥18 years of age) that reported any dimension of physical activity at 3 or more exposure levels and reported risk estimates for depression. Studies were excluded if the sample size was fewer than 3000 participants (to limit small-study effects indicative of publication bias) or if the follow-up period was less than 3 years (to minimize reverse causality bias from undiagnosed depression at baseline). Full details are provided in eTable 1 in the Supplement. Studies were eligible if the exposure included leisure-time physical activity, either alone or combined with other activity domains. The outcome of interest was depression, including (1) presence of major depressive disorder indicated by self-report of physician diagnosis, registry data, or diagnostic interviews using DSM criteria or International Statistical Classification of Diseases and Related Health Problems, Tenth Revision, codes F32 through F33; (2) elevated depressive symptoms established using validated cutoffs for a depressive screening instrument, eg, Center for Epidemiologic Studies Depression scale.

Search and Selection

We searched PubMed, SCOPUS, Web of Science, PsycINFO, and the reference lists of systematic reviews and articles retrieved from the search or known to the authors. Specific search strings are available in eMethods 1 in the Supplement. We considered peer-reviewed articles published in academic journals until November 12, 2020. No language limits were set. Titles and abstracts and subsequently full-text articles were screened by 2 independent reviewers for eligibility, with disagreements resolved by a third reviewer (eFigure 1 in the Supplement).

Data Extraction

Data extraction was carried out using a standardized extraction form and conducted independently by 2 researchers, with disagreements resolved by a third researcher. For each exposure category, we extracted information to estimate physical activity volume if not reported directly, number of depression cases, number of participants and person-years of follow-up, and the effect estimate with 95% CIs. Effect estimates from the most adjusted model were used. Results reported separately for men and women, or for multiple cohorts within a study, were extracted as separate associations.

When the data required for exposure harmonization or meta-analysis were not reported in the original articles, nor obtainable through other publications from the same cohort, we contacted authors before imputing critical information if necessary (eMethods 2 in the Supplement).

Physical Activity Exposure Harmonization

We harmonized reported exposure levels to a common metric of physical activity volume in marginal metabolic equivalent task hours per week (mMET-h/wk), reflecting the energy expended above the resting metabolic rate (1 MET). Combinations of harmonization techniques were used depending on what was reported or obtained from authors and availability of validation work (details of the original exposure data, harmonization methods, and harmonization flowchart are available in an Open Science Framework10 repository).

Frequency, Duration, and Intensity Assumptions

Studies that described physical activity exposure as frequency and duration were converted to a measure of weekly duration. If session duration was not provided, a session duration of 0.75 hours was assumed. Weekly durations were then converted to activity volume in mMET hours per week by multiplying by intensity values of 1.5 mMET for light, 3.5 mMET for moderate and moderate to vigorous, and 7.0 mMET for vigorous physical activity.11 A score of 3.5 mMET corresponds to the midpoint of the range for moderate-intensity activity (≥2 to <5 mMETs). Vigorous was scored at 7 mMET because activity of this intensity contributes twice as much toward meeting physical activity guidelines compared with moderate activity (the World Health Organization advises 150-300 minutes of moderate intensity activity or 75-150 minutes of vigorous intensity activity per week12). Light activities are those in the range of 0.5 to 2 mMETs; however, these were scored slightly higher than the midpoint at 1.5 mMET to reflect the intensity of behaviors commonly referred to as light activity in questionnaires (eg, light housework, light gardening).

Converting Energy Expenditure to MET Units

One study reported energy expenditure without adjustment for body weight (kcal/wk). Weight was calculated using body mass index and height from national survey data.13 The exposure in kcal per week was then divided by derived weight and further converted to MET hours using 1 kcal/kg = 1 MET hour.

Removing the Resting Energy Expenditure Component of Physical Activity

For studies reporting physical activity volume and corresponding duration data, 1 MET hour per week was subtracted for each hour of activity undertaken; this is equivalent to the original studies having multiplied activity duration with net (mMET) rather than gross (MET) intensity values. If duration for each volume level was not reported or obtainable from authors, we used a prediction equation derived from studies where both volume and duration were available (eMethods 3 in the Supplement).

Removing Occupational Component of Aggregate Physical Activity Estimates

Nonoccupational physical activity was used as the exposure in all but 1 study14 that used total physical activity, including leisure-time, domestic, occupational, and transport domains. Given no domain-specific activity results were available for this study, we assumed occupational activity to be 40 hours per week at 1.25 METs (or 0.25 mMETs, ie, equivalent to sedentary office work) and estimated nonoccupational physical activity by subtracting this assumed occupational volume.

Analyses

We conducted dose-response meta-analyses for outcomes with at least 4 independent studies available. For studies reporting only sex-stratified results, strata-specific risks were combined using fixed-effect meta-analysis. Sex-stratified meta-analyses were not conducted because only 2 sets of sex-stratified results were available.

A 2-stage random-effects meta-analysis was used to model the dose-response association between physical activity and depression. First, study-specific associations were estimated using generalized least-squares regression.15,16 Second, we estimated the pooled association by combining the study-specific coefficients using restricted maximum likelihood.17 We assumed the dose-response associations were nonlinear and modeled them by fitting restricted cubic splines with 3 knots at the 0, 37.5th, and 75th percentiles of person-years. If the statistical model was unable to converge for these knots, we increased the percentile for the upper knot until it did.

Estimated Population Risk

To provide a population perspective of the relative importance of the estimated dose-response associations, potential impact fractions (PIFs) were calculated based on the exposure prevalence in the populations of included cohorts.18 PIFs were calculated for 3 exposure levels based on the World Health Organization recommendations for adults12: the minimum recommended level of 8.8 mMET hours per week (volume equivalent to approximately 2.5 h/wk of physical activity at moderate intensity of 3.5 marginal METs), the level recommended for additional health benefits of 17.5 mMET hours per week, as well as 4.4 mMET hours per week (half the minimum recommended level).

Sensitivity and Subgroup Analysis

We conducted 2 sensitivity analyses. First, we assigned 2.5 and 6 mMET to moderate and vigorous physical activity, respectively, when intensity had to be assumed for physical activity exposure harmonization; for studies requiring assumptions about duration of sessions, we assigned a 0.5-hour duration (rather than 0.75-hour duration). Second, we tested the use of knots at the 0, 42.5th, and 85th percentiles in the cubic spline functions.

We conducted subgroup analyses by depression diagnosis method using 2 outcome subtypes: (1) major depression from registry data, self-report of physician diagnosis, or diagnostic interview; (2) use of a screening tool for elevated depressive symptoms with validated cutoffs.

Analysis of Heterogeneity

Cochran Q test and I2 statistic were used to assess heterogeneity of effect measures between studies. We also conducted random-effect meta-regressions (at the physical activity level of 8.8 mMET-h/wk) using restricted maximum likelihood to identify how much of the variance in effect measures was due to proportion of men and women (women only, men only, mixed), outcome ascertainment (registry data, self-report of physician diagnosis, diagnostic interview, Center for Epidemiologic Studies Depression scale, other depressive symptoms scale), duration of follow-up (above or below median of 8.5 years), exclusion of prevalent cases at baseline (yes, no), handling other morbidities at baseline (exclusion of participants, statistical adjustment), and our exposure harmonization approaches (measurement unit conversion, measurement unit conversion plus occupational physical activity assumption, duration/intensity/frequency assumptions). See eTable 2 in the Supplement for more details. We conducted a leave-one-out sensitivity analysis to investigate the influence of each study on the overall effect-size (at the physical activity level of 8.8 mMET-h/wk).

Software, Data, and Code Availability

Analyses were performed using R, version 4.0.5, and the dosresmeta package,19 version 2.0.1. An interactive interface to visualize dose-response curves was developed using the Shiny package, version 1.0.5. Codes for all analyses and the interactive interface are available in a repository on GitHub20 (see README file).

Results

Literature Search

The systematic literature search yielded 19 175 results following removal of duplicates. After excluding 18 827 records based on title and abstract screening, 348 full-text articles were reviewed. Full-text screening identified 15 eligible publications reporting 15 associations including 191 130 participants contributing 28 806 incident depression events and 2 110 588 person-years. Approximately 64% of participants in the studies were women. All but 1 of the included studies originated in high-income countries: 6 from the United States,21,22,23,24,25,26 6 from Europe,27,28,29,30,31,32 1 from Australia,33 1 from Japan,34 and 1 study that included data from India, Ghana, Mexico, and Russia.14 Study characteristics are described in Table 1.14,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43

Table 1. Study Characteristics and Exposure Harmonization.

Source (country) Study name Participants (cases), No.a Sex Age range at baseline, y Follow-up time, y Outcome subtype Outcome ascertainment Exposure harmonizationb
Mikkelsen et al,27 2010 (Denmark) Copenhagen City Heart Study 8804 (436) Women 20-93 ≤26 Major depression Registry Used duration category midpoints. Multiplied by intensity (mMET) for light (1.5) and moderate (3.5) activity.
Paffenbarger et al,24 1994 (United States) Harvard Alumni Health Study 10 201 (387) Men 35-74 23-27 Major depression Self-report of physician diagnosis Interpolated to get median BMI for sample. Mean height for men aged 35-74 y, 173 cm. Used these data to get weight by category. Divided midpoint of absolute activity volume by weight to give kcal/kg/wk. which is MET-h/wk. Equation used to remove resting metabolic rate component to give mMET-h/wk.
Ten Have et al,32 2011 (the Netherlands) Netherlands Mental Health Survey and Incidence Study 3998 (244) Both (49% women) 18-64 3 Major depression Electronic CIDI35 Used duration category midpoints. Multiplied by intensity (mMET) for moderate (3.5) activity.
Chang et al,25 2016 (United States) Nurses’ Health Study 21 728 (3945) Women ≥65 ≤10 Major depression Self-report of depressive symptoms, use of antidepressants, physician diagnosis Used duration category midpoints. Multiplied by intensity (mMET) for moderate (3.5) activity.
Cabello et al,14 2017 (Ghana, India, Mexico, Russia) World Health Organization’s Study on Global Aging and Adult Health 5970 (594) Both (67% women) ≥18 3-8 Major depression World Mental Health Survey Initiative version of the CIDI,36 self-report of physician diagnosis, medication, or treatment Assigned MET-h/wk based on scoring described and used equation to remove resting metabolic rate component. Instrument does not capture sedentary behaviors, so no volume removed for assumed occupational physical activity.
Fernandez-Montero et al,31 2020 (Spain) Seguimiento Universidad de Navarra 15 488 (870) Both (60% women) 18-64 Mean (SD): 8.5 (4.4) Major depression Self-report of physician diagnosis Used volume category midpoints. Equation used to remove resting metabolic rate component to give mMET-h/wk.
Hallgren et al,29 2019 (Sweden) Swedish National March Cohort Study 20 594 (289) Both (65% women) ≥18 ≤13 Major depression Registry Used duration category midpoints. Multiplied by intensity (mMET) for moderate (3.5) activity.
España-Romero et al,21 2013 (United States) Aerobics Center Longitudinal Study 5110 (641) Both (20% women) 20-83 Mean (range): 6.1 (1-12) Elevated depressive symptoms 10-Item CES-D scale37 Used volume category midpoints. Equation used to remove resting metabolic rate component to give mMET-h/wk.
Camacho et al,23 1991 (United States) Alameda County Study 3664 (733) Men ≥20 9 Elevated depressive symptoms Human Population Laboratory depression index38 Article reported activity index 0-14 based on frequency of 4 behaviors; 14 represents often engaging in sports, swimming or walking, doing exercises, and gardening. The mMET-h/wk for a score of 14 was calculated by assuming frequency and duration and by assigning mMET based on description of activities. Midpoints of the index ranges of the 3 categories (0-4 [2], 5-8 [6.5], 9-14 [11.5]) were taken and assigned mMET-h/wk relative to the known score of an index of 14.
Pavey et al,33 2013 (Australia) Australian Longitudinal Study on Women’s Health 12 094 (2419) Women 45-53 3 Elevated depressive symptoms 10-Item CES-D scale37 Used volume category midpoints. Equation used to remove resting metabolic rate component to give mMET-h/wk.
Wise et al,22 2006 (United States) Black Women’s Health Study 35 224 (9465) Women 21-69 3-5 Elevated depressive symptoms 20-Item CES-D scale39 Used duration category midpoints. Multiplied by intensity (mMET) for vigorous (7) activity.
Hamer et al,28 2009 (England) English Longitudinal Study of Aging 4323 (348) Both (53% women) ≥50 4 Elevated depressive symptoms 8-Item CES-D scale40 Used frequency category midpoints. Assumed session duration of 0.75 h. Multiplied by intensity (mMET) for light (1.5), moderate (3.5), and vigorous (7) activity.
Kuwahara et al,34 2015 (Japan) Japan Epidemiology Collaboration on Occupational Health Study 29 082 (6177) Both (15% women) 20-64 Mean: 4.7 Elevated depressive symptoms 13-Item cohort specific scale41 Used volume category midpoints. Equation used to remove resting metabolic rate component to give mMET-h/wk.
Harvey et al,30 2018 (Norway) Health Study of Nord-Trøndelag County 22 564 (1578) Both (50% women) ≥20 9-13 Elevated depressive symptoms 14-Item Hospital Anxiety and Depression scale42 Used duration category midpoints. Multiplied by intensity (mMET) for moderate (3.5) activity.
Hughes et al,26 2019 (United States) Health Professionals Follow-up Study 6311 (788) Men 40-75 25 Elevated depressive symptoms 5-Item Mental Health Inventory43 Used volume category means. Equation used to remove resting metabolic rate component to give mMET-h/wk.

Abbreviations: BMI, body mass index; CES-D, Center for Epidemiologic Studies Depression scale; CIDI, Composite International Diagnostic Interview; mMET, marginal metabolic equivalent of task.

a

Cases of incident depression.

b

Volume category refers to the description of physical activity volume in the article, eg, <10, 10-20, and >20 MET-h/wk. Duration category refers to the description of physical activity duration in the article, eg, 0, 1-3, and >3 h/wk. In both cases, these ranges had to be converted to point estimates for the meta-analysis (further details available in Open Science Framework10 repository).

Dose-Response Analyses

The majority (78%) of participants reported exposure levels below 17.5 mMET hours per week, with almost all (95%) data below 35 mMET hours per week; participants in studies of elevated depressive symptoms tended to be less active than those in studies of major depression (eFigure 2 in the Supplement). Figure 1 shows an inverse, curvilinear dose-response association between physical activity and depression, with greater differences in risk in the lower-dose region. Relative to adults not reporting any activity, those accumulating half the recommended volume of physical activity (4.4 mMET-h/wk) had 18% (95% CI, 13%-23%) lower risk of depression. Adults accumulating the recommended volume of 8.8 mMET hours per week had 25% (95% CI, 18%-32%) lower risk with diminishing potential benefits and higher uncertainty observed beyond that exposure level (Table 2). Interactive dose-response curves and exposure distributions are available online.44 The dose-response curves for major depression and elevated depressive symptoms showed similar curvilinear relationships (Table 2 and Figure 2).

Figure 1. Association Between Physical Activity and Incidence of Depression.

Figure 1.

Dark line represents the meta-analytical dose-response curve (constrained to be linear beyond upper knot at 75% of person-years). Shaded area displays 95% CI. Vertical dotted lines indicate knots at the 37.5th and 75th percentiles of person-years. I2 = 73.7%; P < .001. Interactive dose-response curves and exposure distributions are available online.44

Table 2. Relative Risk and Potential Impact Fractions of Incident Depression, Major Depression, and Elevated Depressive Symptoms at 3 Physical Activity Levelsa.

Risk association by activity volume [mMET-h/wk], RR (95% CI) Population impact by activity volume [mMET-h/wk], PIF (95% CI), %
4.4 8.8 17.5 4.4 8.8 17.5
Depression 0.82 (0.77-0.87) 0.75 (0.68-0.82) 0.72 (0.64-0.81) 6.38 (4.25-8.63) 11.53 (7.69-15.43) 13.89 (8.44-19.25)
Major depression 0.83 (0.75-0.92) 0.75 (0.64-0.87) 0.74 (0.61-0.88) 2.97 (1.27-4.91) 7.28 (3.36-11.44) 8.04 (2.38-13.82)
Elevated depressive symptoms 0.80 (0.73-0.88) 0.73 (0.64-0.84) 0.70 (0.59-0.84) 9.45 (5.19-13.86) 14.44 (7.88-20.92) 17.01 (8.39-25.24)

Abbreviations: mMET, marginal metabolic equivalent of task; PIF, potential impact fraction; RR, relative risk.

a

Risks are expressed relative to accumulating 0 mMET-h/wk.

Figure 2. Associations Between Physical Activity and Incidence of Major Depression and Elevated Depressive Symptoms.

Figure 2.

Dark lines represent the meta-analytical dose-response curve (constrained to be linear beyond upper knot at 75% of person-years). Shaded area displays 95% CIs. Vertical dotted lines indicate knots at the 37.5th and 75th percentiles of person-years. A, Major depression I2 = 54.2%; P = .01. B, Elevated depressive symptoms I2 = 81.3%; P < .001.

Estimated Population Risk

Based on PIF analyses of prevalences of physical activity and depression outcomes in the included cohorts, 11.5% (95% CI, 7.7%-15.4%) of incident depression could have been prevented if all adults had achieved at least 8.8 mMET hours per week of physical activity (Table 2). PIFs were approximately twice as high for elevated depressive symptoms than for major depression at both 8.8 and 17.5 mMET hours per week and approximately 3 times higher at 4.4 mMET hours per week (Table 2).

Sensitivity Analyses

Using alternative assumptions of 0.5 hours for session duration and 1 MET lower-intensity values for moderate and vigorous physical activity (only when those assumptions were necessary) or alternative placement of knots for the splines did not materially alter dose-response associations (eTable 3 in the Supplement) or PIFs (eTable 4 in the Supplement).

Analysis of Heterogeneity

Differences between proportion of women, exposure harmonization method, outcome ascertainment, duration of follow-up, exclusion of prevalent cases at baseline, or managing other morbidities at baseline did not significantly explain variance in the effect measures of associations between physical activity and depression (eFigure 3 in the Supplement). The leave-one-out sensitivity analysis did not identify any outliers (eFigure 4 in the Supplement).

Discussion

This study reports the first dose-response meta-analysis of associations between physical activity and incident depression to our knowledge. Our results show an inverse curvilinear association with the greatest differences in risk observed between low doses of physical activity, suggesting most benefits are realized when moving from no activity to at least some. Accumulating an activity volume equivalent to 2.5 hours of brisk walking per week was associated with 25% lower risk of depression, and at half that dose, risk was 18% lower compared with no activity. Only minor additional benefits were observed at higher activity levels.

Previous meta-analyses have shown lower risks of depression among adults reporting high vs low physical activity7 and suggested a dose-response association using meta-regression without quantifying the association.8 Our study directly models how risk of depression varies across the physical activity exposure range using a continuous scale rather than cruder categories. We also found that even small volumes of activity were beneficial but go further by quantifying differences in risk for these doses. Our findings therefore have important new implications for health practitioners making lifestyle recommendations, especially to inactive individuals who may perceive the current recommended target as unrealistic.

The associations we observed are likely explained by more than 1 mechanism. Proposed pathways include acute neuroendocrine and inflammatory responses to activity such as activation of the endocannabinoid system (“runner’s high”)45 and longer-term adaptations, including changes in the brain’s neural architecture.46 Psychosocial and behavioral explanations have also been suggested, including improved physical self-perceptions and body image, more social interactions, and the personal development of coping strategies.46 The social aspect of activity participation may operate even at relatively low doses, consistent with the dose-response curve we observed. The role of the environment as a potential moderator of the association between physical activity and depression should also be considered. For example, the use of green space is associated with lower risk of depression,47 with mediation analysis suggesting only part of the association is explained by physical activity.48 Conversely, noise pollution49 and neighborhood deprivation50 might diminish the mental health benefits of activity. Such contextual factors may have contributed to the high level of heterogeneity we observed between studies. The above mechanisms may not operate to the same degree across different types, frequencies, intensities, and contexts of activity,51 and measurement of these activity subdimensions may also differ between studies, potentially resulting in higher heterogeneity. Future work should therefore explore the shape of the dose-response relationship for these aspects of activity in addition to total volume.

In subgroup analyses using results from different cohorts for each outcome, we observed that the dose response with physical activity was similar for major depression and elevated depressive symptoms, as has been reported using binary exposure expressions.7 However, PIFs for elevated depressive symptoms were higher; this is because participants in those studies were generally less active than participants in the studies examining major depression. Examining depressive symptoms on a continuous scale instead of using dichotomized outcome variables may provide further insight regarding the benefits of physical activity to alleviate depressive symptoms, but no study that we reviewed reported results in this way.

To conduct nonlinear dose-response analyses, we only included studies with at least 3 levels of physical activity exposure, a requirement that limited the number of eligible results compared with previous meta-analyses.7,8 Physical activity exposures in primary studies were classified using a variety of methods and expressed in different formats, and this could have resulted in excluding even more studies. A key strength of the present meta-analytical work is our extensive exposure harmonization, maximizing the inclusion of the existing evidence. We did observe high heterogeneity between study estimates, but differences in method of exposure harmonization did not explain variation in effect sizes. Duration of follow-up, proportions of men and women, original analytical choices relating to the handling of other morbidities or prevalent depression at baseline, and method of depression diagnosis were not statistically significant predictors of effect size.

It is still possible that the associations observed in the present analyses could overestimate the role of physical activity because of reverse causality. That is, depression at baseline results in lower physical activity in participants who then go on to record depression during follow-up. For depression diagnosed at baseline, this bias is often mitigated by excluding these prevalent cases from analyses. Two of the studies included in our meta-analysis did not exclude participants with depression at baseline or adjust for these statistically,24,33 and our subgroup analysis indicated (although not statistically significantly so) that effect sizes of this subgroup were stronger than those from studies that excluded prevalent cases. Reverse causality bias may also result from depression that is undiagnosed at baseline, particularly given the often recurring nature of the condition.52 To limit risk of bias from undiagnosed depression, we only included studies with at least 3 years of follow-up. This does not rule out the possibility that physical activity may be causally associated with depression on a shorter time scale. Two studies included in our full-text screening were excluded solely based on follow-up time. A study of 108 000 South Korean adults reported an inverse association for physical activity in the exposure range approximating 7 to 40 mMET hours per week, but no association at very high activity levels.53 A study of 4600 Irish adults reported lower (albeit not statistically significant) odds of depression in those meeting physical activity guidelines compared with those who did not.54

Longer follow-up time limits the risk of reverse causality affecting the results but can also introduce regression dilution bias owing to exposure measurement error caused by true variation in physical activity behavior during follow-up.55 For example, in a recent study, mortality associations estimated using 28 years’ follow-up with a single measure of physical activity at baseline were 2- to 3-fold weaker than those estimated using repeated measures.56 We did not include results using more than 1 measure of physical activity because none were identified that met our inclusion criteria. However, we did observe weaker associations in studies with longer follow-up times, although again, this meta-regression result was not statistically significant. Given the recurrent nature of depression and engagement in physical activity as a treatment for this condition,57 use of repeated measures may be particularly important and would strengthen the evidence base for evaluating the dose-response relationship. Further work is required to understand the effects of different strategies for managing reverse causality and exposure measurement error in studies of depression.

Limitations

A limitation of these analyses is that included studies only used self-reported measures of physical activity, which are subject to recall and social-desirability biases. Our analyses also included relatively limited data at higher physical activity doses. Studies using device-based measures of activity capturing a wider range of exposure and with longer follow-up of incident depression are therefore warranted. There were insufficient studies with stratified results to examine sex, age, or geographical subgroups, and notably, lower- and middle-income countries were underrepresented in included studies.

Conclusions

This meta-analysis found an association between physical activity and incident depression. This suggests substantial mental health benefits can be achieved at physical activity levels even below the public health recommendations, with additional benefit for meeting the minimum recommended target but limited extra benefit beyond that. Assuming causality, 1 in 9 cases of depression might have been prevented if everybody in the population was active at the level of current health recommendations.

Supplement.

eMethods 1. Search strings

eMethods 2. Data imputation procedures

eMethods 3. Estimating the resting component of energy expenditure

eTable 1. Study inclusion and exclusion criteria

eTable 2. Study characteristics for analysis of heterogeneity

eTable 3. Relative risks using alternate assumptions

eTable 4. Potential impact fractions using alternate assumptions

eFigure 1. Study screening and selection flowchart

eFigure 2. Distribution of marginal MET hours per week

eFigure 3. Subgroup analysis

eFigure 4. Leave-one-out sensitivity analysis

eReferences

References

  • 1.Abbafati C, Abbas KM, Abbasi-Kangevari M, et al. ; GBD 2019 Diseases and Injuries Collaborators . Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2020;396(10258):1204-1222. doi: 10.1016/S0140-6736(20)30925-9 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 2.Walker ER, McGee RE, Druss BG. Mortality in mental disorders and global disease burden implications: a systematic review and meta-analysis. JAMA Psychiatry. 2015;72(4):334-341. doi: 10.1001/jamapsychiatry.2014.2502 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 3.Patel V, Chisholm D, Parikh R, et al. Addressing the burden of mental, neurological, and substance use disorders: key messages from Disease Control Priorities, 3rd edition. Lancet. 2016;387(10028):1672-1685. doi: 10.1016/S0140-6736(15)00390-6 [DOI] [PubMed] [Google Scholar]
  • 4.Salazar de Pablo G, Solmi M, Vaquerizo-Serrano J, et al. Primary prevention of depression: an umbrella review of controlled interventions. J Affect Disord. 2021;294:957-970. doi: 10.1016/j.jad.2021.07.101 [DOI] [PubMed] [Google Scholar]
  • 5.Teychenne M, Ball K, Salmon J. Physical activity and likelihood of depression in adults: a review. Prev Med. 2008;46(5):397-411. doi: 10.1016/j.ypmed.2008.01.009 [DOI] [PubMed] [Google Scholar]
  • 6.Mammen G, Faulkner G. Physical activity and the prevention of depression: a systematic review of prospective studies. Am J Prev Med. 2013;45(5):649-657. doi: 10.1016/j.amepre.2013.08.001 [DOI] [PubMed] [Google Scholar]
  • 7.Schuch FB, Vancampfort D, Firth J, et al. Physical activity and incident depression: a meta-analysis of prospective cohort studies. Am J Psychiatry. 2018;175(7):631-648. doi: 10.1176/appi.ajp.2018.17111194 [DOI] [PubMed] [Google Scholar]
  • 8.Dishman RK, McDowell CP, Herring MP. Customary physical activity and odds of depression: a systematic review and meta-analysis of 111 prospective cohort studies. Br J Sports Med. 2021;55(16):926-934. doi: 10.1136/bjsports-2020-103140 [DOI] [PubMed] [Google Scholar]
  • 9.Smith AD, Crippa A, Woodcock J, Brage S. Physical activity and incident type 2 diabetes mellitus: a systematic review and dose-response meta-analysis of prospective cohort studies. Diabetologia. 2016;59(12):2527-2545. doi: 10.1007/s00125-016-4079-0 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 10.Open Science Framework . Physical activity and risk of depression: dose-response meta-analysis of prospective studies. Published 2022. Accessed January 21, 2022. https://osf.io/rc9nq/?view_only=7b516d9693d94c1d971512033a54176e
  • 11.Ainsworth BE, Haskell WL, Herrmann SD, et al. 2011 Compendium of physical activities: a second update of codes and MET values. Med Sci Sports Exerc. 2011;43(8):1575-1581. doi: 10.1249/MSS.0b013e31821ece12 [DOI] [PubMed] [Google Scholar]
  • 12.Bull FC, Al-Ansari SS, Biddle S, et al. World Health Organization 2020 guidelines on physical activity and sedentary behaviour. Br J Sports Med. 2020;54(24):1451-1462. doi: 10.1136/bjsports-2020-102955 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 13.Najjar MF, Rowland M. Anthropometric reference data and prevalence of overweight, United States, 1976-80. Vital Health Stat 11. 1987;(238):1-73. [PubMed] [Google Scholar]
  • 14.Cabello M, Miret M, Caballero FF, et al. The role of unhealthy lifestyles in the incidence and persistence of depression: a longitudinal general population study in four emerging countries. Global Health. 2017;13(1):18. doi: 10.1186/s12992-017-0237-5 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 15.Greenland S, Longnecker MP. Methods for trend estimation from summarized dose-response data, with applications to meta-analysis. Am J Epidemiol. 1992;135(11):1301-1309. doi: 10.1093/oxfordjournals.aje.a116237 [DOI] [PubMed] [Google Scholar]
  • 16.Orsini N, Li R, Wolk A, Khudyakov P, Spiegelman D. Meta-analysis for linear and nonlinear dose-response relations: examples, an evaluation of approximations, and software. Am J Epidemiol. 2012;175(1):66-73. doi: 10.1093/aje/kwr265 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 17.Gasparrini A, Armstrong B, Kenward MG. Multivariate meta-analysis for non-linear and other multi-parameter associations. Stat Med. 2012;31(29):3821-3839. doi: 10.1002/sim.5471 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 18.Zapata-Diomedi B, Barendregt JJ, Veerman JL. Population attributable fraction: names, types and issues with incorrect interpretation of relative risks. Br J Sports Med. 2018;52(4):212-213. doi: 10.1136/bjsports-2015-095531 [DOI] [PubMed] [Google Scholar]
  • 19.Crippa A, Orsini N. Multivariate dose-response meta-analysis: the dosresmeta R package. J Stat Softw. 2016;72(code snippet 1):1-15. doi: 10.18637/jss.v072.c01 [DOI] [Google Scholar]
  • 20.GitHub . Meta-analysis repository. Published 2022. Accessed February 9, 2022. https://github.com/meta-analyses/meta-analysis
  • 21.España-Romero V, Artero EG, Lee DC, et al. A prospective study of ideal cardiovascular health and depressive symptoms. Psychosomatics. 2013;54(6):525-535. doi: 10.1016/j.psym.2013.06.016 [DOI] [PubMed] [Google Scholar]
  • 22.Wise LA, Adams-Campbell LL, Palmer JR, Rosenberg L. Leisure time physical activity in relation to depressive symptoms in the Black Women’s Health Study. Ann Behav Med. 2006;32(1):68-76. doi: 10.1207/s15324796abm3201_8 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 23.Camacho TC, Roberts RE, Lazarus NB, Kaplan GA, Cohen RD. Physical activity and depression: evidence from the Alameda County Study. Am J Epidemiol. 1991;134(2):220-231. doi: 10.1093/oxfordjournals.aje.a116074 [DOI] [PubMed] [Google Scholar]
  • 24.Paffenbarger RS Jr, Lee I-M, Leung R. Physical activity and personal characteristics associated with depression and suicide in American college men. Acta Psychiatr Scand Suppl. 1994;377:16-22. doi: 10.1111/j.1600-0447.1994.tb05796.x [DOI] [PubMed] [Google Scholar]
  • 25.Chang SC, Pan A, Kawachi I, Okereke OI. Risk factors for late-life depression: a prospective cohort study among older women. Prev Med. 2016;91:144-151. doi: 10.1016/j.ypmed.2016.08.014 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 26.Hughes KC, Gao X, Molsberry S, Valeri L, Schwarzschild MA, Ascherio A. Physical activity and prodromal features of Parkinson disease. Neurology. 2019;93(23):e2157-e2169. doi: 10.1212/WNL.0000000000008567 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 27.Mikkelsen SS, Tolstrup JS, Flachs EM, Mortensen EL, Schnohr P, Flensborg-Madsen T. A cohort study of leisure time physical activity and depression. Prev Med. 2010;51(6):471-475. doi: 10.1016/j.ypmed.2010.09.008 [DOI] [PubMed] [Google Scholar]
  • 28.Hamer M, Molloy GJ, de Oliveira C, Demakakos P. Leisure time physical activity, risk of depressive symptoms, and inflammatory mediators: the English Longitudinal Study of Ageing. Psychoneuroendocrinology. 2009;34(7):1050-1055. doi: 10.1016/j.psyneuen.2009.02.004 [DOI] [PubMed] [Google Scholar]
  • 29.Hallgren M, Nguyen TTD, Lundin A, et al. Prospective associations between physical activity and clinician diagnosed major depressive disorder in adults: a 13-year cohort study. Prev Med. 2019;118(118):38-43. doi: 10.1016/j.ypmed.2018.10.009 [DOI] [PubMed] [Google Scholar]
  • 30.Harvey SB, Øverland S, Hatch SL, Wessely S, Mykletun A, Hotopf M. Exercise and the prevention of depression: results of the HUNT cohort study. Am J Psychiatry. 2018;175(1):28-36. doi: 10.1176/appi.ajp.2017.16111223 [DOI] [PubMed] [Google Scholar]
  • 31.Fernandez-Montero A, Moreno-Galarraga L, Sánchez-Villegas A, et al. Correction to: Dimensions of leisure-time physical activity and risk of depression in the “Seguimiento Universidad de Navarra” (SUN) prospective cohort. BMC Psychiatry. 2020;20(1):132. doi: 10.1186/s12888-020-02549-5 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 32.Ten Have M, de Graaf R, Monshouwer K. Physical exercise in adults and mental health status findings from the Netherlands mental health survey and incidence study (NEMESIS). J Psychosom Res. 2011;71(5):342-348. doi: 10.1016/j.jpsychores.2011.04.001 [DOI] [PubMed] [Google Scholar]
  • 33.Pavey TG, Peeters G, Bauman AE, Brown WJ. Does vigorous physical activity provide additional benefits beyond those of moderate? Med Sci Sports Exerc. 2013;45(10):1948-1955. doi: 10.1249/MSS.0b013e3182940b91 [DOI] [PubMed] [Google Scholar]
  • 34.Kuwahara K, Honda T, Nakagawa T, et al. Associations of leisure-time, occupational, and commuting physical activity with risk of depressive symptoms among Japanese workers: a cohort study. Int J Behav Nutr Phys Act. 2015;12(1):119. doi: 10.1186/s12966-015-0283-4 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 35.Smeets RM, Dingemans PMA. Composite International Diagnostic Interview (CIDI), version 1.1. World Health Organization; 1993. [Google Scholar]
  • 36.Kessler RC, Ustün TB. The World Mental Health (WMH) Survey Initiative Version of the World Health Organization (WHO) Composite International Diagnostic Interview (CIDI). Int J Methods Psychiatr Res. 2004;13(2):93-121. doi: 10.1002/mpr.168 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 37.Andresen EM, Malmgren JA, Carter WB, Patrick DL. Screening for depression in well older adults: evaluation of a short form of the CES-D (Center for Epidemiologic Studies Depression Scale). Am J Prev Med. 1994;10(2):77-84. doi: 10.1016/S0749-3797(18)30622-6 [DOI] [PubMed] [Google Scholar]
  • 38.Kaplan GA, Roberts RE, Camacho TC, Coyne JC. Psychosocial predictors of depression: prospective evidence from the human population laboratory studies. Am J Epidemiol. 1987;125(2):206-220. doi: 10.1093/oxfordjournals.aje.a114521 [DOI] [PubMed] [Google Scholar]
  • 39.Radloff LS. The CES-D scale: a self-report depression scale for research in the general population. Appl Psychol Meas. 1977;1(3):385-401. doi: 10.1177/014662167700100306 [DOI] [Google Scholar]
  • 40.Karim J, Weisz R, Bibi Z, ur Rehman S. Validation of the eight-item center for epidemiologic studies depression scale (CES-D) among older adults. Curr Psychol. 2015;34(4):681-692. doi: 10.1007/s12144-014-9281-y [DOI] [Google Scholar]
  • 41.Yakura N. Verification of the validity of depressive symptom scale based on the existing health questionnaire. Asian Pac J Dis Manag. 2009;3(1):21-26. doi: 10.7223/apjdm.3.21 [DOI] [Google Scholar]
  • 42.Zigmond AS, Snaith RP. The hospital anxiety and depression scale. Acta Psychiatr Scand. 1983;67(6):361-370. doi: 10.1111/j.1600-0447.1983.tb09716.x [DOI] [PubMed] [Google Scholar]
  • 43.Berwick DM, Murphy JM, Goldman PA, Ware JE Jr, Barsky AJ, Weinstein MC. Performance of a five-item mental health screening test. Med Care. 1991;29(2):169-176. doi: 10.1097/00005650-199102000-00008 [DOI] [PubMed] [Google Scholar]
  • 44.MRC Epidemiology Unit . Physical activity meta-analyses interface. Published 2022. Accessed February 9, 2022. https://shiny.mrc-epid.cam.ac.uk/meta-analyses-physical-activity/
  • 45.Brellenthin AG, Crombie KM, Hillard CJ, Koltyn KF. Endocannabinoid and mood responses to exercise in adults with varying activity levels. Med Sci Sports Exerc. 2017;49(8):1688-1696. doi: 10.1249/MSS.0000000000001276 [DOI] [PubMed] [Google Scholar]
  • 46.Kandola A, Ashdown-Franks G, Hendrikse J, Sabiston CM, Stubbs B. Physical activity and depression: towards understanding the antidepressant mechanisms of physical activity. Neurosci Biobehav Rev. 2019;107:525-539. doi: 10.1016/j.neubiorev.2019.09.040 [DOI] [PubMed] [Google Scholar]
  • 47.Dzhambov AM, Hartig T, Tilov B, Atanasova V, Makakova DR, Dimitrova DD. Residential greenspace is associated with mental health via intertwined capacity-building and capacity-restoring pathways. Environ Res. 2019;178:108708. doi: 10.1016/j.envres.2019.108708 [DOI] [PubMed] [Google Scholar]
  • 48.van den Berg MM, van Poppel M, van Kamp I, et al. Do physical activity, social cohesion, and loneliness mediate the association between time spent visiting green space and mental health? Environ Behav. 2019;51(2):144-166. doi: 10.1177/0013916517738563 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 49.Dzhambov AM, Lercher P. Road traffic noise exposure and depression/anxiety: an updated systematic review and meta-analysis. Int J Environ Res Public Health. 2019;16(21):E4134. doi: 10.3390/ijerph16214134 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 50.Wickham S, Taylor P, Shevlin M, Bentall RP. The impact of social deprivation on paranoia, hallucinations, mania and depression: the role of discrimination social support, stress and trust. PLoS One. 2014;9(8):e105140. doi: 10.1371/journal.pone.0105140 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 51.Lubans D, Richards J, Hillman C, et al. Physical activity for cognitive and mental health in youth: a systematic review of mechanisms. Pediatrics. 2016;138(3):e20161642. doi: 10.1542/peds.2016-1642 [DOI] [PubMed] [Google Scholar]
  • 52.Burcusa SL, Iacono WG. Risk for recurrence in depression. Clin Psychol Rev. 2007;27(8):959-985. doi: 10.1016/j.cpr.2007.02.005 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 53.Kim SY, Park JH, Lee MY, Oh KS, Shin DW, Shin YC. Physical activity and the prevention of depression: a cohort study. Gen Hosp Psychiatry. 2019;60:90-97. doi: 10.1016/j.genhosppsych.2019.07.010 [DOI] [PubMed] [Google Scholar]
  • 54.McDowell CP, Dishman RK, Hallgren M, MacDonncha C, Herring MP. Associations of physical activity and depression: results from the Irish Longitudinal Study on Ageing. Exp Gerontol. 2018;112:68-75. doi: 10.1016/j.exger.2018.09.004 [DOI] [PubMed] [Google Scholar]
  • 55.Strain T, Wijndaele K, Sharp SJ, Dempsey PC, Wareham N, Brage S. Impact of follow-up time and analytical approaches to account for reverse causality on the association between physical activity and health outcomes in UK Biobank. Int J Epidemiol. 2020;49(1):162-172. doi: 10.1093/ije/dyz212 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 56.Lee DH, Rezende LFM, Ferrari G, et al. Physical activity and all-cause and cause-specific mortality: assessing the impact of reverse causation and measurement error in two large prospective cohorts. Eur J Epidemiol. 2021;36(3):275-285. doi: 10.1007/s10654-020-00707-3 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 57.Stubbs B, Vancampfort D, Smith L, Rosenbaum S, Schuch F, Firth J. Physical activity and mental health. Lancet Psychiatry. 2018;5(11):873. doi: 10.1016/S2215-0366(18)30343-2 [DOI] [PubMed] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Supplement.

eMethods 1. Search strings

eMethods 2. Data imputation procedures

eMethods 3. Estimating the resting component of energy expenditure

eTable 1. Study inclusion and exclusion criteria

eTable 2. Study characteristics for analysis of heterogeneity

eTable 3. Relative risks using alternate assumptions

eTable 4. Potential impact fractions using alternate assumptions

eFigure 1. Study screening and selection flowchart

eFigure 2. Distribution of marginal MET hours per week

eFigure 3. Subgroup analysis

eFigure 4. Leave-one-out sensitivity analysis

eReferences


Articles from JAMA Psychiatry are provided here courtesy of American Medical Association

RESOURCES