Trajectories of Physical Activity Predict the Onset of Depressive Symptoms but Not Their Progression: A Prospective Cohort Study
This prospective, community-based study examined trajectories of physical activity from childhood to adulthood and whether these trajectories contributed to depressive symptoms in adulthood to a greater degree than adulthood physical activity. Participants () were from the ongoing Cardiovascular Risk in Young Finns Study which started in 1980. Depressive symptoms were measured with Beck Depression Inventory (BDI-II) in 2012, and physical activity was assessed from 1980 to 2011 with self-reports. Analyses were adjusted for age, sex, childhood negative emotionality, socioeconomic factors, previous depressive symptoms, social support, body mass index, and smoking status (1980–2007). Highly, moderately, and lightly physically active trajectory groups were identified. Highly active participants reported lower levels of depressive symptoms compared to lightly active ones () and compared to moderately active ones (). Moderately active participants had less symptoms than lightly active ones (). High levels of adulthood physical activity associated with lower levels of depressive symptoms (). The findings did not withstand adjustment for previous depressive symptoms (). Lifelong physical activity trajectories or adulthood physical activity was not associated with the progression of depressive symptoms in adulthood. Thus, physical activity history does not contribute to the progression of the depressive symptoms to a greater degree than adulthood physical activity.
Depressive disorders affect 340 million people (approximately 9%) globally regardless of gender, culture, or ethnicity . From 1990 till present time, depression has been rated one of the leading causes of years lived with disability (YLDs) . Depression associates with a wide range of subsequent health outcomes. For instance, depression predicts future coronary events and cardiac deaths in healthy individuals and in those with established cardiovascular disease (CHD) [3, 4]. Depression has also been shown to be associated with several other physical diseases, including cancer and osteoporosis , as well as with severe psychiatric disorders such as schizophrenia . Furthermore, depression relates to low educational attainment and reduced work productivity . Identifying factors that contribute to depression and developing tools that can be applied in early prevention of depression have been prominently highlighted .
Even though some conflicting findings concerning the association between physical activity and depression exist [8, 9], prior evidence provides support for the value of physical activity in reducing depressive symptoms in both healthy and clinical populations . Literature has also indicated that even low doses of physical activity may protect against depression . The relation between physical activity and depression may also be bidirectional, as physical activity may alleviate depressive symptoms, but these symptoms may also decrease the likelihood of initiating physical activity [8, 11]. Physical activity has also been related to neurobiological functioning, which may preempt depressive mood . In addition, physical activity has been shown to be associated with the development of individual qualities such as self-confidence, emotional self-regulation , and capacity of social bonding , which, in turn, may be protective against depression [13, 14].
Several studies have indicated that regular physical activity associates with reduced risk of depression in all age groups from early childhood to late adulthood [8, 10, 15]. There is also evidence that childhood physical activity may lead to maintenance of lifetime physical activity patterns conducive to well-being (e.g., ). These patterns, however, have a tendency to decline towards later adulthood . Many previous findings regarding the association between physical activity and depression derive from cross-sectional designs . Previous longitudinal or intervention designs cover relatively short follow-up phases of physical activity [10, 18]. Prospective studies examining the association between lifelong physical activity and depressive symptoms are very rare .
It has been stated that more sophisticated methods are needed when studying the development and etiology of depression . It has also been shown that people’s health behaviors may have various developmental tendencies instead of a one growth trajectory [20, 21], and some of these behavioral profiles may be more detrimental to health than others. For instance, it is theoretically possible that physical activity’s decline towards later life  may start earlier and/or be steeper among specific groups than in others, and thus the different physical activity profiles may associate with distinct health outcomes. The need for examinations assessing the development and decline of physical activity has been acknowledged [17, 22]. It is also essential to study whether long-term physical activity contributes to health over and above the concurrent physical activity. Furthermore, identifying people, whose behavioral tendencies associate with decreased health, might be useful for professionals who aim to adjust lifestyle interventions to targeted groups. It is also important to gain information of the behavioral determinants and contributors of good health .
Growth mixture modeling provides an appropriate framework for studying developmental processes [20, 24]. The key advantage of such modeling is that it allows for the estimation of interindividual variability in intraindividual patterns of change over time [20, 24]. Current approaches to growth modeling have also been regarded as flexible for instance in terms of including partially missing data, unequally spaced measurement points, nonnormally distributed or discretely scaled repeated measures, and linearly or nonlinearly shaped trajectories . These issues, or some of these, typically emerge in developmental research. It has been denoted that these issues present challenges for more traditional methods .
Extensive adjustment for potential confounders has been lacking in many previous studies regarding physical activity and health outcomes [16, 18]. Development of depression often depends on the interplay of multiple psychological and lifestyle related factors  that ideally should be considered when examining associations between health behaviors and depression [25–30]. Child’s early emotional experiences (e.g., negative emotionality) may have effects on the development of depression later in life . Childhood family’s socioeconomic status has been shown to be associated with the development of health behaviors and well-being . In adulthood, occurrence of depressive episodes has been shown to predict subsequent depression . It has been proven that depression is related to age , and women tend to report higher levels of depression than men . Adulthood socioeconomic position, experiences of social support, body mass index, and smoking status are also associated with depression in adulthood [26, 29, 30].
We examined the potential heterogeneity in physical activity trajectories in relation to depressive symptoms over a 30-year period from childhood to adulthood in a population-based sample with six cohorts and including eight study waves. To our knowledge, this is the first study assessing physical activity and depressive symptoms in such design. The analyses were adjusted for several risk factors for depression [25–30]. Physical activity measurements were performed during participants’ childhood and adulthood (from the age of 9 to 49), and depressive symptoms were assessed in participants’ adulthood (participants aged from 35 to 50). To gain a comprehensive picture of the association between physical activity and depressive symptoms, the association was studied cross-sectionally, with respect to change and with respect to long-term trajectories. The objectives of this study were (1) to explore the potentially distinct trajectories of physical activity from childhood to adulthood, (2) to examine whether physical activity was associated with cross-sectional, short-term, and long-term changes in depressive symptoms in adulthood, and (3) to examine whether the physical activity trajectories contributed to depressive symptoms in adulthood to a greater degree than adulthood physical activity.
2.1. Study Design and Participants
The study participants were from the ongoing prospective Cardiovascular Risk in Young Finns Study that began in 1980 . The original sample consisted of 3596 children and adolescents (83.2% of those invited, 1832 females and 1764 males) from six birth cohorts (aged 3, 6, 9, 12, 15, and 18). To obtain a representative sample, Finland was divided into five areas based on the locations of universities with medical schools (Helsinki, Kuopio, Oulu, Tampere, and Turku), and the participants were randomly selected based on their social security numbers from nearby urban and rural areas. Informed consent was requested from each participant (or from the parents of small children), and the study was approved by the local ethics committees. The study was conducted according to Declaration of Helsinki and American Psychological Association’s ethical principles.
The sample was followed in 8 waves, 1983, 1986, 1989, 1992, 1997, 2001, 2007-2008, and 2012, in which medical, psychological, and physical activity studies were conducted. Physical activity from childhood to middle adulthood was assessed in 1980, 1983, 1986, 1989, 1992, 2001, 2007, and 2011, response rate ranging from 53.1% to 72.8% () of the original study participants (Table 1). Depressive symptoms were measured in 2012 and 47.9% () of the original study subjects participated in the examination (Table 1).
Previous studies of sample attrition have shown that there has not been systematic selection bias regarding study participants’ medical profiles or physical activity [31, 32], but some selective attrition with respect to personality and depressive symptoms exists [33, 34]. Participants who were less self-directed and less agreeable and had higher levels of neuroticism as well as depressive symptoms had discontinued the study more often than others [33, 34].
2.2.1. Physical Activity
Physical activity was assessed with self-administered questionnaires. Information concerning 3- and 6-year-old children’s (born in 1974 and 1977) physical activity levels is missing from this study, because they were not able to self-report their physical activity levels in 1980. Participants, from whom physical activity was assessed, were aged from 9 to 18 in 1980 and from 34 to 49 in 2011.
From 1980 to 1989, physical activity questionnaires consisted of five questions focusing on the intensity and frequency of participants’ leisure time physical activity, participation in sports-club training, participation in sports competitions, and participants’ usual way of spending leisure time . From 1992 till present time, questionnaires consisted of five questions as well, assessing the intensity and frequency of leisure time physical activity, hours spent on physical activity per week, average duration of a physical activity session, and participation in organized physical activity . From 1980 to 1992, the answers for the questions were coded into 3 categories (ranging from 1 to 3), excluding the items that assessed participation in sports competitions (1980–1989) and participation in organized sports (1992) which had a response range from 1 to 2. From 2001 to 2011, all responses to the questions were coded into 3 categories (response scale ranging from 1 to 3). A sum score (physical activity index) of question responses was created for each participant each year (Table 1), higher scores reflecting higher physical activity level. The index has been found reliable and valid .
2.2.2. Depressive Symptoms
Participants’ depressive symptoms were assessed with Beck Depression Inventory II (BDI-II) . These symptoms were measured in 2012 when the participants were aged from 35 to 50. BDI-II consists of 21 symptoms with a severity range from 0 (no symptoms) to 3 (severe level of depressive symptoms). A sum score of all items was computed for each participant (Table 1), and no missing items were allowed. The reliability estimate (Cronbach’s α) for the depressive symptom scores was >0.90. BDI-II has demonstrated to be a valid instrument [36–38], and it has been regarded as an acknowledged standard in the measurement of depressive mood [36–39]. It is applicable in clinical and nonclinical contexts, including in general populations [36–39]. BDI-II correlates highly with its earlier versions, including modified BDI [36, 38], which has also been considered as a valid measure for assessing depressive symptoms in general populations . Furthermore, BDI-II correlates well with other widely used scales for depression . The instrument has been designed and also demonstrated to be a useful screening tool for potential depressed cases [36, 38].
Childhood, adulthood, and general (age, sex, and body mass index) covariates were controlled for in this study (Table 1) [25–30]. Participants’ negative emotionality  was reported by the primary caretaker via six questions reflecting participants’ behavior in childhood (e.g., “The child hits/kicks other children “accidentally””), on a scale from 1 (true) to 2 (not true), and average of the items was calculated for each participant. As some of the participants were adolescents in 1980, their caretakers responded to this question retrospectively. Participants’ parents’ socioeconomic status was assessed via two indices, educational and income levels . Parents’ educational level was determined via educational information collected from participants’ mothers’ and fathers’ (1 = below 9 years/comprehensive school, 2 = 9 to 12 years/secondary school, and 3 = over 12 years/academic education). If parents’ educational information differed, we based parental educational status on the information collected from the parent with the higher educational level. If educational information was available for only one parent, family’s educational status was determined using his/her educational information. Family’s income level was rated in an 8-point scale [1 =< 15 000 marks (~2523 euros) and 8 => 100 000 marks (~16819 euros)].
Symptoms of participants’ adulthood depression were determined in 1992, 1997, 2001, and 2007 via a modified version of Beck Depression Inventory, referred to as a modified BDI [34, 40]. Items of the measure were rated in a 5-point scale, and average of the items was computed each year for each participant. Participants’ socioeconomic status (2007) was determined via two indices; education was assessed via a 3-category scale (1 = comprehensive school, 2 = secondary school, and 3 = academic level) and income level via an 8-point scale (1 =< 10 000 euros and 8 => 70 000 euros). Additionally, participants’ experiences of social support, body mass index, and smoking status measured in 2007 were controlled for [29, 30]. Social support was determined via a 12-question inventory  using a 5-point scale, and a mean score of the items was calculated for each participant. Participants’ smoking status was examined via a 5-category scale (1 = smokes a cigarette per day or more, 2 = smokes once in a week, 3 = smokes less than once in a week, 4 = has quitted smoking, and 5 = has never smoked).
2.3. Statistical Analyses
Physical activity questionnaires, which were designed for children and adolescents (1980–1989) and adults (1992–2011) differed slightly in their content. To assure that the findings of the present study were based on changes in physical activity and not due to a measurement artifact, a confirmatory factor model was used to examine whether the physical activity indices consisting of five indicator variables had measurement and structural invariance over time [42, 43]. Weighted least squares means and variance adjusted (WLSMV) estimation was used for all analyses . The goodness of fit for scalar invariance was assessed with comparative fit index (CFI), Tucker-Lewis index (TLI), and root-mean square error of approximation index (RMSEA). Factor scores derived from this examination were used in subsequent analyses.
Within the growth mixture modeling framework, Latent Class Growth Analysis (LCGA) was used to explore the trajectories of physical activity from childhood to adulthood. LCGA captures information about developmental processes at inter- and intraindividual levels, detecting subpopulations with distinct growth trajectories . The determination of the number of subgroups for physical activity was based on Akaike’s Information Criterion (AIC) . In addition, the determination of the groups was based on the classification quality estimations and practical considerations [20, 45]. Within the LCGA model, the average temporal trajectories in the physical activity groups were modeled by regression equations, in which both the linear and quadratic terms were tested for the independent variable (time).
The associations between physical activity factor scores (assessed from age 9 to 49) and depressive symptoms (participants aged from 35 to 50) were first examined cross-sectionally and longitudinally with linear regression analyses. Due to the potential multiple testing problem, Bonferroni-corrected values () were used in determining the significant associations. Thereafter, the associations between the physical activity trajectory groups and depressive symptoms measured in 2012 were examined with analyses of variance, and post hoc tests were also performed (Bonferroni’s method). Furthermore, we examined the longitudinal associations between physical activity levels assessed in participants’ adulthood (2007, including participants aged from 30 to 45) and depressive symptoms (2012) using a linear regression. Due to the number of missing values, the variance analyses and the regression analyses in which the adulthood physical activity (2007) was used as a predictor were performed in another dataset which was imputed using the expectation-maximization (EM) algorithm . Analyses were performed in statistical software programs Mplus (version 7.1 and version 7.2), IBM SPSS (version 21), and Stata (version 13). values of <0.05 were considered significant.
Descriptives of the original sample () are shown in Table 1. Supplementary Tables and (see Supplementary Material available online at http://dx.doi.org/10.1155/2016/8947375) provide descriptives of the sample in the complete () and imputed () data, respectively. Although the scalar invariance model for physical activity did not demonstrate strong factorial invariance over time, the fit for partial scalar invariance model was adequate (CFI = 0.90, TLI = 0.90, and RMSEA = 0.047), given partial invariance of the threshold parameters. For RMSEA, values <0.05 indicate a very close model fit, and CFI and TLI values close to 0.90 denote an adequate fit . Since the partial scalar invariance model was considered acceptable, factor scores were predicted for each subject to be used in subsequent analyses (Supplementary Table ).
The prerequisites of growth mixture modeling (e.g., ) were met satisfactorily. LCGA suggested that a three-factor (group) solution was the best fitting model for the data based on AIC indices (AIC = 48915.44, 48915.39, 48901.43, 48902.85, and 48907.08 for 1, 2, 3, 4, and 5 groups, resp.; a lower AIC implies a better model). The classification quality of the three-factor model was also adequate based on the average probability estimates (group 1 = 0.74, group 2 = 0.86, and group 3 = 0.71). Values >0.70 delineate that the group consists of individuals with similar patterns of change . Linear parameter estimates were found for the physical activity groups (Supplementary Table ). Figure 1 shows the central tendencies of the three groups, lightly (), moderately (), and highly () physically active groups by age. Based on the estimated marginal means, participants’ physical activity levels remained relatively unchanged from childhood to adulthood in each group, although these levels appeared to diminish minimally towards middle adulthood in all participants (Figure 1).
The regression analyses performed in the original sample () indicated that low levels of physical activity factor scores were, in most examined ages, associated with higher levels of depressive symptoms in participants’ adulthood (), although some of the associations attenuated when Bonferroni-corrected values () were applied (Table 2).
The analyses of variance in the original sample () indicated that the physical activity trajectory groups predicted symptoms of depression [ (2,1719 = 8.12, , adjusted = 0.01)]. Post hoc tests showed that highly physically active group had lower levels of depressive symptoms than lightly active group (mean difference = −3.26; , 95%CI: −5.25 to −1.26). Highly physically active participants had lower levels of depression in comparison to moderately active ones (mean difference = −1.92, , 95%CI: −3.66 to −0.18). Also moderately active participants had lower levels of depressive symptoms than lightly active participants (mean difference = −1.33, , 95%CI: −2.48 to −0.19). Based on the analyses of variance after adjusting for the covariates [25–30], the associations attenuated to nonsignificance [ (2, 631) = 2.13, , adjusted = 0.25]. In particular the previous symptoms of depression measured in 1992, 1997, 2001, and 2007 attenuated the associations (for details, see Supplementary Table ).
We also examined the unadjusted associations in the data with full information on all study variables (), in which case the association between the physical activity trajectory groups and depressive symptoms became only marginally significant [ (2,645) = 2.76, , adjusted = 0.01]. As the sample attrition might have affected this association (i.e., by reducing the statistical power), the analyses were also performed in a dataset which was imputed using EM-algorithm (). This method was applied as the Little’s MCAR test  confirmed that the data were not missing completely at random ( = 3903.02, = 2923, ). In the imputed dataset, physical activity groups predicted the symptoms of depression [ (2,3561) = 16.74, , adjusted = 0.01]. Post hoc tests indicated that highly physically active group had lower levels of depressive symptoms than lightly active group (mean difference = −2.74; , 95%CI: −3.91 to −1.56) (Figure 2). Highly physically active participants had lower levels of depressive symptoms in comparison to moderately active ones (mean difference = −1.59, , 95%CI: −2.60 to −0.57) (Figure 2). Also moderately active participants had lower levels of depressive symptoms than lightly active participants (mean difference = −1.15, , 95%CI: −1.81 to −0.49) (Figure 2). When the covariates [25–30] were adjusted for, the results became nonsignificant [ (2,3547) = 0.53, , adjusted = 0.47]. In particular the previous symptoms of depression measured in 1997 and 2001 attenuated the associations (for details, see Supplementary Table ).
Thereafter, we examined the longitudinal associations between adulthood physical activity (assessed in 2007, participants’ aged from 30 to 45) and symptoms of depression (assessed in 2012). In the original sample, adulthood physical activity was associated with depressive symptoms ( = −1.10, , 95%CI: −1.58 to −0.61, adjusted = 0.01). When the covariates [25–30] were adjusted for, the association attenuated to nonsignificance ( = −0.17, , 95%CI: −0.85 to 0.51, adjusted = 0.24). In particular the previous symptoms of depression assessed in 2001 and 2007 attenuated the association (for details, see Supplementary Table ). Thereafter, the analyses were performed in a sample with full information on all study variables (), in which case the adulthood physical activity was not associated with depressive symptoms ( = −0.55, , 95%CI: −1.31 to 0.21, adjusted = 0.002). Due to the sample attrition, the association was studied also in the imputed data (), and the results showed that the adulthood physical activity was associated with decreased levels of depressive symptoms ( = −1.07, , 95%CI: −1.36 to −0.78, adjusted = 0.01). When the covariates [25–30] were adjusted for, the association attenuated to nonsignificance ( = 0.06, , 95%CI: −0.16 to 0.28, adjusted = 0.47). In particular the symptoms of depression assessed in 2001 and 2007 attenuated the association (for details, see Supplementary Table ).
This study examined whether distinct trajectories of lifelong physical activity existed in the data and whether physical activity was related to depressive symptoms in adulthood. We also studied whether the lifelong physical activity trajectories contributed to the outcome to a greater degree than adulthood physical activity. This inspection was important, because the information whether lifelong trajectories contribute to depressive symptoms over and above the concurrent (adulthood) physical activity is lacking. LCGA revealed three distinct groups, the lightly, moderately, and highly physically active groups. Physical activity levels remained relatively similar from childhood to adulthood in each group, although the groups’ physical activity levels decreased slightly towards middle adulthood. These results are in line with previous studies demonstrating the decline of physical activity with age .
High physical activity associated with lower level of depressive symptoms in adulthood, which is in accord with previous research [10, 18]. The mechanisms behind the observed association can be both physiological and psychological. Physical activity may alleviate depressive symptoms through neurobiological alterations . Physical activity may also relate to enhanced self-confidence, emotion regulation skills, and social capacities that are associated with positive mood [13, 14].
The association between high physical activity and lower levels of depressive symptoms in adulthood, however, disappeared when covariates (age, sex, childhood negative emotionality, parents’ socioeconomic status, participants’ previous symptoms of depression, participants’ socioeconomic status, social support, body mass index, and smoking status) were taken into account. The explanatory power of the fully adjusted models was substantially higher in comparison to the unadjusted ones (46% higher in the imputed data). These same results were found when lifelong physical activity trajectories and adulthood physical activity (assessed in 2007) were used as predictors for depressive symptoms in adulthood. Specifically, physical activity did not predict lower levels of depressive symptoms in adulthood when the experiences of previous depression were taken into account. The bidirectional nature of the association between physical activity and depression has been documented [8, 11]. Hence, it is theoretically possible that the childhood symptoms of depression preceding our first measurement have decreased the likelihood of initiating physical activity [8, 11] or that of maintaining adequate level of physical activity .
Summarizing, the study indicated that lifelong physical activity trajectories or adulthood physical activity levels were not associated with the progression of depressive symptoms in adulthood. The study suggests that lifelong physical activity history does not contribute to the progression of the depressive symptoms to a greater degree than adulthood physical activity.
4.1. Limitations and Strengths
Our results need to be interpreted in light of the following limitations. Participants did not provide information regarding each variable across the measurement years, which diminished the complete-sample size considerably. However, we attempted to control the potential attrition bias via an imputation method. The measurements concerning the participants’ physical activity focused only on the self-reported leisure time physical activity, not total physical activity or energy expenditure. Also the depressive symptoms were assessed with self-administered questionnaires, and thus the possibility of subjective bias cannot be excluded. However, this is common in epidemiological studies by necessity. To our knowledge, comparably long follow-up studies using diagnostic interviews do not exist. Furthermore, early childhood depression was not assessed in the study, but this deficiency was retrieved to an extent by controlling for the childhood negative emotionality .
The strengths of the study were the prospective, population-based study design, relatively large sample size, use of LCGA, and utilization of an extensive set of covariates. In LCGA, the main advantage is that it allows a researcher to model developmental processes at inter- and intraindividual levels . Such modeling has been shown to be especially useful regarding the development of health behaviors [20, 21], and the need for further studies has been recognized . However, in the case of this study, the traditional cross-sectional and longitudinal studies yielded to very similar results. This was not evident a priori, and therefore our study adds evidence regarding the etiology of depression  by showing that physical activity trajectories appear not to play a special role over and above the adulthood physical activity.
Furthermore, depressive symptoms were studied with a well-validated instrument, BDI-II [36–39].
We were able to study the whole variance of participants’ depressive symptoms instead of categorical diagnoses. This is important when attempting to find preventive tools for depression, because functional impairment is much more strongly associated with symptom severity than with diagnostic symptom count . Also the people who have experienced symptoms of depression are at risk of getting a diagnosis .
This study identified three distinct physical activity groups, the lightly, moderately, and highly physically active ones. Each group’s physical activity levels remained relatively unchanged from childhood to adulthood, although the levels tended to diminish slightly towards later adulthood. Physical activity was associated with depressive symptoms in adulthood. Highly physically active participants from childhood to adulthood had lower levels of depressive symptoms in adulthood compared to lightly physically active ones. Furthermore, participants’ adulthood physical activity assessed in 2007 was associated with decreased levels of depressive symptoms in adulthood. The associations between physical activity and depressive symptoms disappeared when the preexisting symptoms of depression were controlled for, indicating that physical activities did not associate with the progression of depressive symptoms. Thus, the study suggests that lifelong physical activity history does not contribute to the progression of the depressive symptoms to a greater degree than adulthood physical activity. Obtaining information of mental health history might benefit clinicians and other professionals in evaluating the role of physical activity in well-being.
The sponsors had no role in study design, collection, analysis, interpretation of the data, writing the report, and decision to submit the report for publication.
The authors declare that there is no conflict of interests regarding the publication of this paper.
The Young Finns Study has received financial support from the Academy of Finland, Social Insurance Institution of Finland, Kuopio, Tampere and Turku University Hospital Medical Funds, Juho Vainio Foundation, Paavo Nurmi Foundation, Finnish Foundation of Cardiovascular Research and Finnish Cultural Foundation, Sigrid Juselius Foundation, Tampere Tuberculosis Foundation, and Emil Aaltonen Foundation. This study was also supported by Academy of Finland [Grant nos. 258711, 265869, and 258578 (Liisa Keltikangas-Järvinen, Laura Pulkki-Råback, Kaisa Kaseva, Mirka Hintsanen, and Tom Rosenström)], Signe and Ane Gyllenberg Foundation (Liisa Keltikangas-Järvinen, Mirka Hintsanen, and Laura Pulkki-Råback), Emil Aaltonen Foundation (Mirka Hintsanen), Alli Paasikivi Foundation (Mirka Hintsanen), Juho Vainio Foundation (Laura Pulkki-Råback), Finnish Ministry of Education and Culture (Tuija Tammelin, Xiaolin Yang, and Kaisa Kaseva), Finnish Cultural Foundation (Christian Hakulinen), and Urheiluopistosäätiö (Kaisa Kaseva). Irina Lisinen and Mauri Niiniaho are gratefully acknowledged for data management.
The Supplementary Material consists of six tables. Supplementary Tables 1-3 present the descriptives of the data. Supplementary Table 4 presents the parameter estimates of the physical activity trajectory groups. Supplementary Tables 5-6 present the regression analyses between physical activity and symptoms of depression.
P. C. Dinas, Y. Koutedakis, and A. D. Flouris, “Effects of exercise and physical activity on depression,” Irish Journal of Medical Science, vol. 180, no. 2, pp. 319–325, 2011.View at: Publisher Site | Google Scholar
T. Vos, A. D. Flaxman, M. Naghavi et al., “Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010,” The Lancet, vol. 380, no. 9859, pp. 2163–2196, 2010.View at: Publisher Site | Google Scholar
R. M. Carney and K. E. Freedland, “Depression, mortality, and medical morbidity in patients with coronary heart disease,” Biological Psychiatry, vol. 54, no. 3, pp. 241–247, 2003.View at: Publisher Site | Google Scholar
D. M. Clarke and K. C. Currie, “Depression, anxiety and their relationship with chronic diseases: a review of the epidemiology, risk and treatment evidence,” Medical Journal of Australia, vol. 190, no. 7, pp. 54–60, 2009.View at: Google Scholar
K. W. Sax, S. M. Strakowski, P. E. Keck Jr., V. H. Upadhyaya, S. A. West, and S. L. Mcelroy, “Relationships among negative, positive, and depressive symptoms in schizophrenia and psychotic depression,” British Journal of Psychiatry, vol. 168, no. 1, pp. 68–71, 1996.View at: Publisher Site | Google Scholar
J. M. Donohue and H. A. Pincus, “Reducing the societal burden of depression: a review of economic costs, quality of care and effects of treatment,” PharmacoEconomics, vol. 25, no. 1, pp. 7–24, 2007.View at: Publisher Site | Google Scholar
R. F. Muñoz, P. Cuijpers, F. Smit, A. Z. Barrera, and Y. Leykin, “Prevention of major depression,” Annual Review of Clinical Psychology, vol. 6, pp. 181–212, 2010.View at: Publisher Site | Google Scholar
S. M. Pereira, M. Geoffroy, and C. Power, “Investigation of bidirectional associations between depressive symptoms and physical activity from age 11 to 50 years in the 1958 British Birth Cohort,” The Lancet, vol. 382, supplement 80, 2013.View at: Publisher Site | Google Scholar
D. Kritz-Silverstein, E. Barrett-Connor, and C. Corbeau, “Cross-sectional and prospective study of exercise and depressed mood in the elderly: The Rancho Bernardo study,” American Journal of Epidemiology, vol. 153, no. 6, pp. 596–603, 2001.View at: Publisher Site | Google Scholar
M. Teychenne, K. Ball, and J. Salmon, “Physical activity and likelihood of depression in adults: a review,” Preventive Medicine, vol. 46, no. 5, pp. 397–411, 2008.View at: Publisher Site | Google Scholar
M. Azevedo Da Silva, A. Singh-Manoux, E. J. Brunner et al., “Bidirectional association between physical activity and symptoms of anxiety and depression: the Whitehall II study,” European Journal of Epidemiology, vol. 27, no. 7, pp. 537–546, 2012.View at: Publisher Site | Google Scholar
I. Helmich, A. Latini, A. Sigwalt et al., “Neurobiological alterations induced by exercise and their impact on depressive disorders,” Clinical Practice and Epidemiology in Mental Health, vol. 6, pp. 115–125, 2010.View at: Publisher Site | Google Scholar
L. Crust, “A review and conceptual re-examination of mental toughness: implications for future researchers,” Personality and Individual Differences, vol. 45, no. 7, pp. 576–583, 2008.View at: Publisher Site | Google Scholar
E. McAuley, B. Blissmer, D. X. Marquez, G. J. Jerome, A. F. Kramer, and J. Katula, “Social relations, physical activity, and well-being in older adults,” Preventive Medicine, vol. 31, no. 5, pp. 608–617, 2000.View at: Publisher Site | Google Scholar
J. M. Ruuskanen and I. Ruoppila, “Physical activity and psychological well-being among people aged 65 to 84 years,” Age and Ageing, vol. 24, no. 4, pp. 292–296, 1995.View at: Publisher Site | Google Scholar
R. Malina, “Tracking of physical activity across the lifespan,” President's Council on Physical Fitness and Sports Research Digest, vol. 3, no. 14, pp. 1–8, 2001.View at: Google Scholar
R. Telama and X. Yang, “Decline of physical activity from youth to young adulthood in Finland,” Medicine and Science in Sports and Exercise, vol. 32, no. 9, pp. 1617–1622, 2000.View at: Google Scholar
X. Yang, M. Hirvensalo, M. Hintsanen et al., “Longitudinal associations between changes in physical activity and depressive symptoms in adulthood: the Young Finns study,” International Journal of Behavioral Medicine, vol. 21, no. 6, pp. 908–917, 2014.View at: Publisher Site | Google Scholar
P. De Jonge and A. M. Roest, “Depression and cardiovascular disease: the end of simple models,” British Journal of Psychiatry, vol. 201, no. 5, pp. 337–338, 2012.View at: Publisher Site | Google Scholar
B. O. Muthén and L. K. Muthén, “Integrating person-centered and variable-centered analyses: growth mixture modeling with latent trajectory classes,” Alcoholism: Clinical and Experimental Research, vol. 24, no. 6, pp. 882–891, 2000.View at: Publisher Site | Google Scholar
N. R. Riggs, C.-P. Chou, C. Li, and M. A. Pentz, “Adolescent to emerging adulthood smoking trajectories: when do smoking trajectories diverge, and do they predict early adulthood nicotine dependence?” Nicotine & Tobacco Research, vol. 9, no. 11, pp. 1147–1154, 2007.View at: Publisher Site | Google Scholar
S. C. Duncan, T. E. Duncan, L. A. Strycker, and N. R. Chaumeton, “A cohort-sequential latent growth model of physical activity from ages 12 to 17 years,” Annals of Behavioral Medicine, vol. 33, no. 1, pp. 80–89, 2007.View at: Publisher Site | Google Scholar
N. M. Peel, R. J. McClure, and H. P. Bartlett, “Behavioral determinants of healthy aging,” American Journal of Preventive Medicine, vol. 28, no. 3, pp. 298–304, 2005.View at: Publisher Site | Google Scholar
P. J. Curran, K. Obeidat, and D. Losardo, “Twelve frequently asked questions about growth curve modeling,” Journal of Cognition and Development, vol. 11, no. 2, pp. 121–136, 2010.View at: Publisher Site | Google Scholar
K. A. Dodge, “Social-cognitive mechanisms in the development of conduct disorder and depression,” Annual Review of Psychology, vol. 44, no. 1, pp. 559–584, 1993.View at: Publisher Site | Google Scholar
B. Galobardes, M. Shaw, D. A. Lawlor, J. W. Lynch, and G. D. Smith, “Indicators of socioeconomic position (part 1),” Journal of Epidemiology and Community Health, vol. 60, no. 1, pp. 7–12, 2006.View at: Publisher Site | Google Scholar
J. Karsten, C. A. Hartman, J. H. Smit et al., “Psychiatric history and subthreshold symptoms as predictors of the occurrence of depressive or anxiety disorder within 2 years,” British Journal of Psychiatry, vol. 198, no. 3, pp. 206–212, 2011.View at: Publisher Site | Google Scholar
J. Mirowsky, “Age and the gender gap in depression,” Journal of Health and Social Behavior, vol. 37, no. 4, pp. 362–380, 1996.View at: Publisher Site | Google Scholar
A. Allgöwer, J. Wardle, and A. Steptoe, “Depressive symptoms, social support, and personal health behaviors in young men and women,” Health Psychology, vol. 20, no. 3, pp. 223–227, 2001.View at: Publisher Site | Google Scholar
F. S. Luppino, L. M. de Wit, P. F. Bouvy et al., “Overweight, obesity, and depression: a systematic review and meta-analysis of longitudinal studies,” Archives of General Psychiatry, vol. 67, no. 3, pp. 220–229, 2010.View at: Publisher Site | Google Scholar
O. T. Raitakari, M. Juonala, T. Rönnemaa et al., “Cohort profile: the cardiovascular risk in young Finns study,” International Journal of Epidemiology, vol. 37, no. 6, pp. 1220–1226, 2008.View at: Publisher Site | Google Scholar
X. Yang, R. Telama, M. Hirvensalo et al., “Leadership component of type A behavior predicts physical activity in early midlife,” International Journal of Behavioral Medicine, vol. 19, no. 1, pp. 48–55, 2012.View at: Publisher Site | Google Scholar
T. Rosenström, M. Jokela, C. R. Cloninger et al., “Associations between dimensional personality measures and preclinical atherosclerosis: the cardiovascular risk in Young Finns study,” Journal of Psychosomatic Research, vol. 72, no. 5, pp. 336–343, 2012.View at: Publisher Site | Google Scholar
T. Rosenström, M. Jokela, S. Puttonen et al., “Pairwise measures of causal direction in the epidemiology of sleep problems and depression,” PLoS ONE, vol. 7, no. 11, Article ID e50841, 2012.View at: Publisher Site | Google Scholar
R. Telama, J. Viikari, and I. Välimäki, “Atherosclerosis precursors in Finnish children and adolescents. X. Leisure-time physical activity,” Acta Paediatrica Scandinavica. Supplement, vol. 74, no. 318, pp. 169–180, 1985.View at: Publisher Site | Google Scholar
A. Beck, R. Steer, and G. Brown, Manual for the Beck Depression Inventory-II, Psychological Corporation, San Antonio, Tex, USA, 1996.
E. A. Storch, J. W. Roberti, and D. A. Roth, “Factor structure, concurrent validity, and internal consistency of the beck depression inventory—second edition in a sample of college students,” Depression and Anxiety, vol. 19, no. 3, pp. 187–189, 2004.View at: Publisher Site | Google Scholar
Y.-P. Wang and C. Gorenstein, “Psychometric properties of the Beck Depression Inventory-II: a comprehensive review,” Revista Brasileira de Psiquiatria, vol. 35, no. 4, pp. 416–431, 2013.View at: Publisher Site | Google Scholar
A. Beck, R. Steer, and G. Brown, BDI-II—Beck Depression Inventory®-II: Finnish Standardization, Psychological Corporation and Psykologien Kustannus Oy, Helsinki, Finland, 2004.
A. T. Beck, C. H. Ward, M. Mendelson, J. Mock, and J. Erbaugh, “An inventory for measuring depression,” Archives of General Psychiatry, vol. 4, pp. 561–571, 1961.View at: Publisher Site | Google Scholar
G. Zimet, N. Dahlem, S. Zimet et al., “The multidimensional scale of perceived social support,” Journal of Personality Assessment, vol. 52, no. 1, pp. 30–41, 1988.View at: Publisher Site | Google Scholar
R. J. Vandenberg and C. E. Lance, “A review and synthesis of the measurement invariance literature: suggestions, practices, and recommendations for organizational research,” Organizational Research Methods, vol. 3, no. 1, pp. 4–70, 2000.View at: Publisher Site | Google Scholar
L. Muthén and B. Muthén, Mplus User's Guide, Muthén & Muthén, Los Angeles, Calif, USA, 6th edition, 1998–2010.
S. I. Vrieze, “Model selection and psychological theory: a discussion of the differences between the Akaike information criterion (AIC) and the Bayesian information criterion (BIC),” Psychological Methods, vol. 17, no. 2, pp. 228–243, 2012.View at: Publisher Site | Google Scholar
H. Andruff, N. Carraro, A. Thompson et al., “Latent class growth modelling: a tutorial,” Tutorials in Quantitative Methods for Psychology, vol. 5, no. 1, pp. 11–24, 2009.View at: Google Scholar
A. P. Dempster, N. M. Laird, and D. B. Rubin, “Maximum likelihood from incomplete data via the EM algorithm,” Journal of the Royal Statistical Society, vol. 39, no. 1, pp. 1–38, 1977.View at: Google Scholar | MathSciNet
G. W. Cheung and R. B. Rensvold, “Evaluating goodness-of-fit indexes for testing measurement invariance,” Structural Equation Modeling. A Multidisciplinary Journal, vol. 9, no. 2, pp. 233–255, 2002.View at: Publisher Site | Google Scholar | MathSciNet
R. J. Little, “A test of missing completely at random for multivariate data with missing values,” Journal of the American Statistical Association, vol. 83, no. 404, pp. 1198–1202, 1988.View at: Publisher Site | Google Scholar | MathSciNet
J. Karsten, C. A. Hartman, J. Ormel, W. A. Nolen, and B. W. J. H. Penninx, “Subthreshold depression based on functional impairment better defined by symptom severity than by number of DSM-IV symptoms,” Journal of Affective Disorders, vol. 123, no. 1–3, pp. 230–237, 2010.View at: Publisher Site | Google Scholar