Associations of self-reported physical activity and depression in 10,000 Irish adults across harmonised datasets: a DEDIPAC-study

Background Depression is a prevalent, debilitating, and often recurrent mood disorder for which successful first-line treatments remains limited. The purpose of this study was to investigate the cross-sectional associations between self-reported physical activity (PA) and depressive symptoms and status among Irish adults, using two existing datasets, The Irish Longitudinal Study on Ageing (TILDA) and The Mitchelstown Cohort Study. Methods The two selected databases were pooled (n = 10,122), and relevant variables were harmonized. PA was measured using the short form International Physical Activity Questionnaire. Depressive symptoms were measured by the Center for Epidemiologic Studies Depression (CES-D) questionnaire. Participants were classified as meeting World Health Organization moderate-to-vigorous PA (MVPA) guidelines or not, and divided into tertiles based on weekly minutes of MVPA. A CES-D score of ≥16 indicated elevated depressive symptoms. Data collection were conducted in 2010–2011. Results Significantly higher depressive symptoms were reported by females (7.11 ± 7.87) than males (5.74 ± 6.86; p < 0.001). Following adjustment for age, sex, BMI, and dataset, meeting the PA guidelines was associated with 44.7% (95%CI: 35.0 to 52.9; p < 0.001) lower odds of elevated depressive symptoms. Compared to the low PA tertile, the middle and high PA tertiles were associated with 25.2% (95%CI: 8.7 to 38.6; p < 0.01) and 50.8% (95%CI: 40.7 to 59.2; p < 0.001) lower odds of elevated depressive symptoms, respectively. Conclusion Meeting the PA guidelines is associated with lower odds of elevated depressive symptoms, and increased volumes of MVPA are associated with lower odds of elevated depressive symptoms. Electronic supplementary material The online version of this article (10.1186/s12889-018-5702-4) contains supplementary material, which is available to authorized users.


Background
Depression is a prevalent, debilitating, and often recurrent mood disorder [1,2] that affects over 300 million people globally [1] and is the second leading cause of global disability [3]. The financial burden associated with depression continues to grow, due largely to increasing treatment costs and wider societal costs of lost productivity [4]. Depression is the mental disorder with the highest cost implications in Europe, accounting for over €90 billion in direct and indirect cost across Europe [4].
Successful treatment of depression with traditional first-line treatments, such as psychotherapy and pharmacotherapy, remains limited [5]. Exercise has been show to be a viable treatment for depression and depressive symptoms [6,7]; however, evidence indicates that approximately 63% of people with mental health problems, including depression, do not seek specialist treatment and only half of those who do appear to improve [8]. Thus, there is a continued need to identify and investigate cost-effective prevention and treatment strategies that are effective and more accessible.
The salutary benefits of regular physical activity (PA) for mental health are well-established [9] and without the physiologic side effects and costs associated with antidepressants and psychotherapy [10]. In the context of depression, evidence supports: (i) significantly lower PA levels among adults with depression compared to adults without depression [11,12]; (ii) inverse associations between PA and onset of depression [13,14]; and (iii) that habitual PA may moderate treatment response [15]. However, many studies have not used validated instruments to capture important features of PA such as type, frequency, intensity, and duration, and the use of inconsistent self-report measures of PA between studies has limited the ability to examine dose-response relationships between PA and depression [13]. Despite this, large-scale studies exist that have measured PA and depression, though this research question often was not the primary focus of the original research [16]. Furthermore, in studies that did address PA as potential correlate of depression, different methodological approaches have been adopted, limiting the possibility to generalize the findings [13].
The cost of gathering such data is significant and research funders often emphasize maximising the potential of existing data resources through data sharing and secondary analysis [17][18][19]. The research and knowledge potential of existing datasets can be further exploited through dataset pooling and subsequent harmonisation of variables [17,[20][21][22][23]. Reported benefits of pooling include increased statistical power, potential for subgroup analysis and comparative research, increased exposure heterogeneity and generalisability, and increased cost-efficiency of research programmes [17,[24][25][26].
Thus, the aim of the current investigation was to harmonise and integrate two existing datasets in order to increase statistical power and generalisability in investigating and quantifying the association between meeting PA guidelines and elevated depressive symptoms, and to explore whether larger associations are observed with increased PA, in a large cross-sectional population sample of Irish adults. It is hypothesised that meeting PA guidelines will be associated with lower odds of elevated depressive symptoms, with lower odds of elevated depressive symptoms observed with greater volumes of PA.

Methods
This study used Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) recommendations to guide reporting [27].

Selection of participating studies
The platform for the current research was the DEterminants of DIet and Physical ACtivity Knowledge Hub (DEDIPAC-KH), a pan-European joint programme initiative active from 2013 to 2016. The focus of DEDIPAC was to better understand how individual and contextual factors influence diet, PA, and sedentary behaviour [28,29]. Within DEDIPAC, a compendium of relevant and accessible datasets regarding the listed health risk behaviours and their determinants was developed [30]. Datasets for inclusion in the compendium were identified by: i) a search of the CORDIS project platform, which is the European Commission's primary public repository and portal to disseminate information on all EU-funded research projects and their results, ii) an examination of existing recent literature reviews and noting the nature of datasets used, and iii) recommendation from the DEDIPAC network and expert consultation. The compendium provides information on the study, the study population, the determinants and health behaviours measured and level of accessibility. The compendium is accessible through https://www.dedipac.eu.
On examination of the compendium, 23 datasets were noted that measured both PA behaviour and a psychological determinant of some nature (e.g. motivation, perceived competence, anxiety, depression etc.). An evaluation of the available data dictionaries which were either publicly available or secured directly from the targeted dataset owners revealed five relevant datasets (three Irish, one Dutch, and one multi-nation) which had potential for harmonisation (i.e., PA and depression measured using a same or similar methodology). Target variables (including relevant covariates) were agreed and listed by three authors (CMcD, CMacD, MPH). The dataset owners were subsequently approached for access to the specific target variables. Of these, all agreed to participate; however, the multi-national study was inaccessible within the timeframe of the study. Therefore, the Dutch study was also excluded due to its smaller sample and nationality that differed from the remaining three studies. Another study (the Cork Children's Lifestyle Survey) was subsequently excluded as it assessed lifetime diagnosis of depression rather than current depressive symptoms. Appropriate data access agreements were signed as required. Only cross-sectional data for datasets were available. Thus, two datasets, The Irish Longitudinal Study on Ageing (TILDA) and The Mitchelstown Cohort Study together consisting of 10,122 Irish adults, were harmonised for the present study.

Participating datasets
The Irish Longitudinal Study on Ageing (TILDA) is a large (n = 8329) prospective cohort study that assesses the social, economic, and health circumstances of community dwelling adults aged ≥50 years and their partners of any age in Ireland [31]. For the current study, the first wave of data collected in 2011 was used. PA was measured using the short International Physical Activity Questionnaire (IPAQ-SF) [32], in which respondents are asked to report the number of days and the duration of the vigorous, moderate, and walking activities they undertook during the last seven days. Depression was measured with the Center for Epidemiologic Studies Depression (CES-D) questionnaire [33].
The Mitchelstown Cohort Study provided data (n = 1804) from the second phase of the Cork and Kerry Diabetes and Heart Disease Study [34], designed to assess diabetes and heart disease among 55-74 year-old adults. Data collection began in 2010. PA was measured with the IPAQ-SF while depression was measured using the CES-D.

Harmonisation process
A target list of relevant variables that may influence the association between PA and depression, based on a logical, theoretical, or prior empirical relation with PA and/or depression, and that had potential for data harmonisation (i.e., same or similar measure) was generated from the two included datasets. Appropriate variable harmonisation methods for PA, depression, and participant age, sex, and body mass index (BMI) were agreed by a six-person expert consensus group from within DEDIPAC (Table 1). These variables were then harmonised to form a large PA and mental health integrated dataset.

Data analysis
PA data were processed according to the IPAQ guidelines for data processing and analysis [35], and respondents were categorised as meeting or not meeting World Health Organization PA guidelines (i.e., ≥150 min.week − 1 of MVPA, ≥75 min.week − 1 of vigorous PA, or ≥ 500 MET.min.week − 1 ) [36]. As the PA data were positively skewed, within individual datasets and across the integrated dataset participants were divided into tertiles based on weekly minutes of MVPA (i.e., 0 min.week − 1 , 10-295 min.week − 1 , and ≥ 300 min.week − 1 ). CES-D scores range from 0 to 60, and individuals at risk for clinical depression were indicated by scores ≥16 [33]; this cut-off score has demonstrated 100% sensitivity and 87.6% specificity among adults [37]. Respondents above this cut-off score are referred to as having elevated depressive symptoms throughout the current manuscript as they are at increased risk of depression; however, meeting the cut-off score does not necessarily indicate the presence of clinically diagnosed depression. BMI was categorized as underweight, normal weight, overweight, and obese using established thresholds of < 18.5 kg.m − 2 , 18.5-24.99 kg.m − 2 , 25-29.99 kg.m − 2 , and ≥ 30 kg.m − 2 , respectively [38]. Statistical analyses were conducted using SPSS Version 22.0 (Armonk, NY: IBM Corp.).

Within dataset analyses
For descriptive purposes, Chi-square tests examined differences in proportions of elevated depressive symptom status, ten-year age categories, sex, BMI categories, and meeting PA guidelines between datasets. For significant Chi-square tests, Z tests were calculated for column proportions for each row in the Chi-square contingency table and adjusted using a Bonferroni correction to provide further detail regarding PA and depression across the different levels of covariates [39]. Within each dataset, binomial logistic regression quantified crude and adjusted associations between PA and elevated depressive symptoms. Covariates in adjusted models were age, sex, and BMI.

Integrated dataset analysis
Chi-square tests examined differences in elevated depressive symptom status, age categories, sex, and BMI categories, between individuals meeting and not meeting PA guidelines. Binomial logistic regression quantified crude and adjusted associations between PA and elevated depressive symptom status. Covariates in adjusted models were age, sex, BMI, and dataset. Likelihood ratio tests examined covariate significance. One-way ANO-VAs followed by Bonferroni-corrected post-hoc tests quantified differences in depressive symptoms between those meeting and not meeting PA guidelines, PA tertiles, and sexes. Hedges' g effect sizes and associated 95% confidence intervals (95%CI) were calculated to quantify the magnitude of differences in depressive symptom scores between meeting PA guidelines, PA tertiles, and sexes [40]. Effect sizes of 0.2, 0.5, and 0.8 can been used as cut-offs for small, medium, and large effects, respectively [41].

Participant characteristics within datasets
Harmonised data characteristics are presented in Table 2.  Table 2.
Briefly, a significantly greater proportion of Mitchelstown participants reported elevated depressive symptoms, did not meet PA guidelines, and were aged 50-59 years and male (all p < 0.05).

Participant characteristics across integrated dataset
Participant characteristics by PA guideline adherence in the integrated dataset are presented in Table 4. Elevated depressive symptom status (χ 2 (1, N = 10,133) = 94.60, p < 0.001), sex (χ 2 (1, N = 10,133) = 243.30, p < 0.001), age (χ 2 (4, N = 10,122) = 168.63, p < 0.001), and BMI (χ 2 (3, N = 7818) = 22.03, p < 0.001) significantly differed between meeting and not meeting PA guidelines. Briefly, people not meeting PA guidelines were significantly more likely to report elevated depressive symptoms and be female, aged 70-79 or ≥ 80 years, and obese (all p < 0.05). Table 3 presents crude and adjusted odds ratios for the associations between meeting PA guidelines and PA tertiles, and elevated depressive symptoms across the integrated dataset. Odds ratios for covariates in adjusted analyses are reported in Additional file 1: Table S1. Significant crude associations remained robust following adjustment for age, sex, BMI, and dataset. Compared to not meeting PA guidelines, meeting the guidelines was associated with 44.3% (95% CI: 34.5 to 52.6; p < 0.001) lower odds of elevated depressive symptoms. Age, sex, and BMI were significant covariates (all p < 0.005). Significant crude associations remained robust following adjustment for age, sex, BMI, and dataset. The middle and high PA tertiles were associated with 26.4% (95% CI: 10.7 to 39.4; p < 0.01) and 51.1% (95% CI: 40.9 to 59.6; p < 0.001) lower odds of elevated depressive symptoms, respectively. Age, sex, and BMI were significant covariates (all p < 0.05).

Discussion
Using secondary analyses of pooled data from large-scale Irish studies, this work contributes to deepen the  Different subscript letters indicate a subset of each category whose column proportions differ statistically significantly at the .05 level BMI Body mass index, PA Participant characteristics knowledge on the association between meeting PA guidelines and elevated depressive symptoms, and a potential dose-response in the PA-depression relationship, in a European country. A previous umbrella systematic review of psychological determinants of PA across the life course highlighted that the negative association between PA and depression/depressive symptoms' presents limited and no conclusive evidence [42]. The present study supports associations between meeting PA guidelines and lower odds of depression, and a dose response in the PA-depression relationship. Approximately 10.9% of the current sample reported elevated depressive symptoms, supporting previous findings in worldwide [2] and Irish [43] samples. Females were significantly more likely to report higher depressive symptoms (g = 0.18), which is consistent with previous evidence [44]. Adjusted analyses, both within and across datasets, support the association between not meeting MVPA guidelines and elevated depressive symptoms, consistent with previous findings regarding the beneficial effects of meeting MVPA guidelines on depressive symptoms [13]. The 44% lower odds of reporting elevated depressive symptoms among people meeting MVPA guidelines is consistent with previous findings [9,45].
A potential dose-response relationship was observed in the present study. Compared to the low PA tertile, the middle and high PA tertiles were associated with 26 and 51% lower odds of elevated depressive symptoms, respectively, such that participation in the highest volume of MVPA appeared to confer more benefits than a more moderate volume. Previous findings are limited and contradictory, and a clear dose-response relationship between PA and depression has not been readily apparent [13,15]. Research by Dunn et al. (2005) [46] observed similar findings to ours in that significantly lower post-intervention depression scores among adults with mild to moderate major depressive disorder whose total energy expenditure was 17.5 kcal/kg/week for 12-weeks compared to those whose total energy expenditure was 7.0 kcal/kg/week [46]. Further, those who routinely engage in higher levels of PA have been shown to respond better to depression treatment than those who are less active [15]. Collectively, these findings support the need to engage in at least a moderate level of PA, with the potential for increased benefits with increased PA.
Although not examined in the present study, there is also evidence to suggest that compared to engaging in no PA, engaging in volumes of MVPA less than recommended levels [13], and increasing levels of light intensity PA [47], may also protect against depression. This may have important public health implications as, although the current study strongly supports the benefits of meeting PA guidelines, the wider research appears to indicate that any amount of PA is better than none.
The iterative steps for effective retrospective data harmonisation published by Fortier and colleagues (2017) [24] underpinned the harmonisation process within this study (e.g. definition of question, detailed characteristics regarding potential datasets, expert input, valid data input, etc.). In concordance with an effective data harmonisation process, the dataset owners were involved from an early stage and access to detailed information on the data collection procedures adopted by the two studies included was available. The implementation of such methodological rigour optimised the quality of the harmonisation process and, thereby the quality and scientific utility of the final database generated.

Strengths and limitations
The present study benefits from a large sample size, a PA measure that measured important components of PA, such as frequency, duration, and intensity of PA, direct examination of meeting PA guidelines, and following rigorous harmonisation procedures. However, several limitations exist. There is the possibility of differences in sampling methodology across studies, while some differences between studies can be impossible to reconcile. For example, including other lifestyle and health behaviours, socio-economic status, and social support as covariates may influence the association between PA level and depression, yet differences in how they were measured made it impossible to include it as a covariate in the present analyses. Because of this, the field will benefit from a more standardised set of measurements across studies. Over-reporting of PA is frequently encountered in epidemiological research partly due to influences such as social desirability and recall bias [48] and may have led to respondents being misclassified as meeting PA guidelines in the present study; however, there is no evidence to suggest that reporting error should differ according to mental health status. Finally, due to the cross-sectional design of this study, causation cannot be deduced. Nonetheless, the current findings are to the authors' knowledge the first to provide support for the benefits of PA for depression, with indications of a dose-response relation, among a large, Irish sample.

Conclusions
Inverse associations between meeting PA guidelines and odds of reporting elevated depressive symptoms were observed, with stronger associations observed with increased volumes of MVPA. Future prospective cohort studies should examine this potential dose-response relationship and focus on potential moderators and mediators of the association between PA and depression in large, European-based samples. Further, more large-scale prospective studies in large European-based samples using objective measure of PA are required to establish cause and effect and to further elucidate the role of potential moderators and mediators of the PA-depression relationship.

Additional file
Additional file 1: Table S1. Odds ratios (OR) and 95% confidence intervals (CI) derived from binominal logistic regression analyses as indicators of association between physical activity (PA) and covariates and depressive symptoms. The sponsors played no role in study design, methods, subject recruitment, data collection, analysis or preparation of the paper.

Availability of data and materials
The TILDA data that support the findings of this study are publicly available in the Irish Social Science Data Archive [http://www.ucd.ie/issda/data/tilda/]. The Mitchelstown data that support the findings of this study are available from Dr. Janas Harrington, University College Cork, Cork, Ireland but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Additional characteristics and degree of accessibility to the original datasets pooled and harmonized in the current study are defined within the DEDIPAC Compendium of datasets (https://www.dedipac.eu/compendium).

Authors' contributions
CMcD, CMacD, and MH had full access to all data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Study concept and design: CMcD, CMacD, and MH. Acquisition, analysis or interpretation of data: CMcD, CMacD, MH, JL, JB, CD, and JH. Drafting of the manuscript: CMcD, CMacD, MH, AC, LC, CD, AL, and FCML. Critical revision of the manuscript for important intellectual content: All authors. Intellectual content: All authors. Statistical analysis: CMcD and MH. All authors read and approved the final manuscript.

Ethics approval and consent to participate
The data collection procedures used in these original studies were in accordance with the ethical standards of the respective institutional research committees, the 1964 Helsinki Declaration, and its later amendments. Informed consent was obtained from all participants included in the original studies and all data was strictly confidential and anonymous. In the context of this ethical approval, access and permission to utilise the data from the original studies was received from the dataset owners at Trinity College Dublin and University College Cork, Ireland.

Competing interests
The authors declare that they have no competing interests.