Article Text
Abstract
Objective To examine whether occupational physical activity changes predict future body mass index (BMI) changes.
Methods This longitudinal cohort study included adult participants attending ≥3 consecutive Tromsø Study surveys (examinations 1, 2 and 3) from 1974 to 2016 (N=11 308). If a participant attended >3 surveys, the three most recent surveys were included. Occupational physical activity change (assessed by the Saltin-Grimby Physical Activity Level Scale) was computed from the first to the second examination, categorised into persistently inactive (n=3692), persistently active (n=5560), active to inactive (n=741) and inactive to active (n=1315). BMI change was calculated from the second to the third examination (height being fixed at the second examination) and regressed on preceding occupational physical activity changes using analysis of covariance adjusted for sex, birth year, smoking, education and BMI at examination 2.
Results Overall, BMI increased by 0.84 kg/m2 (95% CI 0.82 to 0.89). Following adjustments as described previously, we observed no differences in BMI increase between the occupational physical activity change groups (Persistently Inactive: 0.81 kg/m2, 95% CI 0.75 to 0.87; Persistently Active: 0.87 kg/m2, 95% CI 0.82 to 0.92; Active to Inactive: 0.81 kg/m2, 95% CI 0.67 to 0.94; Inactive to Active: 0.91 kg/m2, 95% CI 0.81 to 1.01; p=0.25).
Conclusion We observed no prospective association between occupational physical activity changes and subsequent BMI changes. Our findings do not support the hypothesis that occupational physical activity declines contributed to population BMI gains over the past decades. Public health initiatives aimed at weight gain prevention may have greater success if focusing on other aspects than occupational physical activity.
- epidemiology
- public health
- longitudinal studies
- preventive medicine
- physical work
Statistics from Altmetric.com
Key messages
What is already known about this subject?
The inconclusive results from observational studies on occupational physical activity change and body mass index (BMI) gain may be due to methodological issues.
What are the new findings?
Occupational physical activity declines were not prospectively associated with BMI gains in this large population-based sample.
How might this impact on policy or clinical practice in the foreseeable future?
Public health initiatives aimed at weight gain prevention may have greater success if focusing on other aspects than occupational physical activity.
Introduction
Excessive adiposity and weight gain arise from an imbalance between energy intake and expenditure.1 Increased energy intake is likely the main driver for population weight gains,2 but declines in physical activity levels may also contribute.1 3 At the population level, it may be easier to prevent weight gain by increasing physical activity levels than changing food habits.1 Although the evidence for a prospective association between physical activity and weight gain is limited by methodological challenges,4 higher levels of physical activity are reported to prevent weight gain at the population level.5
Energy expenditure contribution from occupational physical activity is considered higher than that from leisure-time physical activity.3 6 Since leisure-time physical activity appears stable over the past decades and occupational physical activity has declined in western countries,3 7–10 lower levels of occupational physical activity, rather than leisure-time physical activity, may contribute to population gains in weight.3 11 12
Studies assessing the association between occupational physical activity and body mass index (BMI) or weight show conflicting results.11–16 Some studies reported no association between baseline occupational physical activity and future BMI change11 13–16; however, baseline physical activity does not take the reciprocal relationship of changing weight and physical activity into account (ie, physical activity level at baseline may change over time to follow-up, which may be related or unrelated to weight change).4 Other studies computed change scores for both occupational physical activity and BMI and reported conflicting results12 17; however, without adjusting for previous physical activity or BMI/weight at baseline, this represents a cross-sectional analysis of change scores (ie, it is as likely that physical activity change leads to weight change as vice versa), and thus, the direction of the association is unexamined.4
To overcome these methodological challenges, the aim of this study was to assess whether changes in occupational physical activity predicted future changes in BMI over a 40-year period in a large cohort of Norwegian adults examined at three time points with ~6 years of follow-up between each time point.
Methods
Design
The Tromsø Study is an ongoing population-based cohort study in the municipality of Tromsø, Norway, which includes seven repeated surveys with high attendance (%): 1974 (Tromsø 1) (83%), 1979–1980 (Tromsø 2) (85%), 1986–1987 (Tromsø 3) (81%), 1994–1995 (Tromsø 4) (77%), 2001 (Tromsø 5) (79%), 2007–2008 (Tromsø 6) (66%) and 2015–2016 (Tromsø 7) (65%). The cohort includes invited participants from total birth cohorts and random samples of inhabitants in the Tromsø municipality.10 18 Tromsø 1 included only men, while Tromsø 2–7 included both sexes (details described elsewhere (Tromsø 1–618 and Tromsø 710). In this study, we included participants attending at least three consecutive surveys (hereafter examinations 1–3). We computed change in physical activity from examination 1 to 2 followed by change in BMI and weight from examination 2 to 3. Consequently, the follow-up period for physical activity change from examination 1 to 2 and BMI change from examination 2 to 3 were 6–7 years (mean: 6.5 years) for all included participants. Inclusion criteria were information on (1) physical activity at examination 1 and 2, and height and weight at examination 2 and 3; (2) educational level and smoking habits at examination 2; and (3) not pregnant at examination 2 and/or 3. If participants attended more than three consecutive surveys, data from the three most recent surveys were included in the main analyses (overall cohort), while one participant could be included in multiple period-specific samples (Tromsø 1–3: 1974–1987, Tromsø 2–4: 1979–1995, Tromsø 3–5: 1986–2001, Tromsø 4–6: 1994–2008 and Tromsø 5–7: 2001–2016). The layout for the analyses is illustrated in figure 1.
Participants
A flowchart illustrates the selection of participants for our samples (online supplemental figure 1). In short, the overall cohort comprised 11 308 participants with their three most recent attendances. The period-specific sample sizes were as follows: Tromsø 1–3 (1974–1987): n=3570, Tromsø 2–4 (1979–1995): n=9679, Tromsø 3–5 (1986–2001): n=3827, Tromsø 4–6 (1994–2008): n=2212 and Tromsø 5–7 (2001–2016): n=1146). Each individual was eligible for inclusion in multiple period-specific samples. Some participants were excluded due to missing confounders: Tromsø 1–3 (1974–1987): n=512, Tromsø 2–4 (1979–1995): n=595, Tromsø 3–5 (1986–2001): n=15, Tromsø 4–6 (1994–2008): n=39 and Tromsø 5–7 (2001–2016): n=20 (online supplemental figure 1).
Supplemental material
The descriptive characteristics at examination 2 for the overall cohort and period-specific samples are presented in table 1. Tromsø 1 (1974) included only men; thus, the Tromsø 1–3 (1974–1987) sample included only men. All other cohorts are well balanced on sex distribution. Across period-specific samples, age distribution increases, current smokers decrease and educational-level increase (table 1).
Patient and public involvement
There was no public involvement in the design or implementation of this study. The Tromsø 7 advisory board included patient (University Hospital of North Norway) and public (Norwegian Health Association, Tromsø municipality) representatives, and some participants were invited as ambassadors during data collection where they actively contributed to participant recruitment.
Physical activity
Physical activity was measured using the Saltin-Grimby Physical Activity Level Scale (SGPALS) questionnaire19 20 for occupational and leisure-time physical activity (leisure time during the last 12 months) (four hierarchical levels), slightly modified compared with the original SGPALS from 196819 (differences described in online supplemental file 1; the SGPALS layout is presented in online supplemental table 1). For the occupational SGPALS, those reporting rank 1, predominantly sedentary work, were considered inactive; those reporting rank 2, sitting or standing work with some walking; rank 3, walking and some handling of material; or rank 4, heavy manual work, were considered active (online supplemental table 1). Similar inactive/active categorisation was used for the leisure time SGPALS (online supplemental table 1). The occupational SGPALS has shown acceptable reliability21 and an ability to rank participants compared with accelerometry.22
Change in occupational and leisure time SGPALS was computed as (1) persistently inactive (reporting rank 1 at examinations 1 and 2), (2) persistently active (rank ≥2 at examinations 1 and 2), (3) active to inactive (rank ≥2 at examination 1 and rank 1 at examination 2) and (4) inactive to active (rank 1 at examination 1 and rank ≥2 at examination 2).
The occupational time SGPALS was used in all surveys of the Tromsø Study, while the leisure time SGPALS was used in all except Tromsø 4 (1994–1995). In Tromsø 5 (2001), the leisure time SGPALS was answered by those under 70 years.
BMI and weight
Weight and height were measured in light clothing and expressed as kilogram (kg) and metre (m). BMI at examination 2 was calculated as weight divided by the square height (kg/m2). To eliminate the effect of possible height loss between examination 2 and 3, change in BMI at examination 3 was calculated as weight at examination 3 divided by the square height at examination 2. BMI change is our primary outcome, while weight change results are secondary outcomes (online supplemental tables 2 and 3 and 5–9).
Confounders and effect modifiers
Our selected confounders were sex, birth year, smoking, education and baseline BMI/weight (at examination 2). Effect modifiers included the abovementioned confounders, in addition to leisure-time physical activity change. Smoking (from questionnaire) was categorised into (1) current smoker, (2) previous smoker and (3) never smoker. Years of education (from questionnaire) were reported in Tromsø 2 (1979–1980), Tromsø 3 (1986–1987) and Tromsø 5 (2001), which we categorised into (1) primary school (<10 years), (2) high school (10–12 years), (3) university<4 years (13–15 years) and (4) university≥4 years (≥16 years). A five-group alternative for education was reported in Tromsø 4 (1994–1995) and Tromsø 6 (2007–2008), including the four aforementioned groups and a fifth named ‘technical school 2 years senior high’ (eg, craftsman, plumber, electrician and carpenter), which we categorised as high school. All confounders included in the models were retrieved from examination 2.
Statistical analyses
We used paired t-tests to assess whether participants changed BMI and weight from examination 2 to 3. We used analysis of covariance to assess whether physical activity changes from examination 1 to 2 predicted BMI or weight changes from examination 2 to 3 as overall and in strata of sex, birth year, smoking, education and leisure-time physical activity change, with adjustment for sex, birth year, smoking, education and BMI or weight at examination 2. Q–Q plots confirmed change in BMI and weight from examination 2 to 3 to not deviate from normal distribution. The Levene test of equality variance confirmed homogeneity of variance across occupational physical activity change groups (all p>0.07). We assessed interaction effects between occupational physical activity change and potential effect modifiers (sex, birth year, smoking, education and leisure-time physical activity change from examination 1 to 2) in the overall cohort. For sensitivity analyses, we computed occupational physical activity change into six groups: (1) Persistently Inactive, (2) Persistently Active, (3) active but decreasing (rank 4 or 3→3 or 2), (4) active and increasing (rank 2 or 3→3 or 4), (5) Active to Inactive and (6) Inactive to Active. Data are shown as mean and 95% CIs unless otherwise stated. We used the Statistical Package for Social Sciences V.26 for all statistical analyses.
Results
The participants in the overall cohort and period-specific samples increased their BMI from examination 2 to 3 (all p<0.01) (table 2). Weight change results are shown in online supplemental table 2).
Change in BMI by change in occupational physical activity
Changes in BMI by occupational physical activity change, overall and by strata of sex, birth year, smoking, education and leisure-time physical activity changes are presented in table 3. We observed no differences in BMI change from examination 2 to 3 by occupational physical activity changes from examination 1 to 2 (Persistently Inactive: 0.81 kg/m2, 95% CI 0.75 to 0.87; Persistently Active: 0.87 kg/m2, 95% CI 0.82 to 0.92; Active to Inactive: 0.81 kg/m2, 95% CI 0.67 to 0.94; Inactive to Active: 0.91 kg/m2, 95% CI 0.81 to 1.01; p=0.25), which was consistent in stratified analyses (all p≥0.054) (table 3).
We found no interaction effects of potential effect modifiers for the association between occupational physical activity changes and BMI changes (sex: p=0.87, smoking status: p=0.64, education: p=0.25 and leisure-time physical activity changes: p=0.24), except by birth year (p=0.01).
Overall and stratified weight change results for the overall cohort are found in online supplemental table 3; we found no differences in weight change from examination 2 to 3 by occupational physical activity change from examination 1 to 2 (all p≥0.049).
In the sensitivity analyses, where we computed occupational physical activity change into six groups (1) Persistently Inactive, (2) Persistently Active, (3) Active but Decreasing (rank 4 or 3 to 3 or 2), (4) Active and Increasing (rank 2 or 3 to 3 or 4), (5) Active to Inactive and (6) Inactive to Active, the results generally remained unchanged (overall analysis: p=0.15); however, some differences were observed in some strata analyses (birth year, born ≤1929: p=0.03; education, high school: p=0.04, university≥4 years: p=0.049; and leisure-time physical activity changes; PA: p=0.003) (online supplemental table 4). We found no interaction in the association between occupational physical activity change and BMI change (sex: p=0.21, smoking: p=0.59, education: p=0.88, leisure-time physical activity change: p=0.12, except by birth year: p=0.04).
We observed no differences in BMI change by occupational physical activity change in any period-specific sample (table 4): (1) there were no differences in BMI change from Tromsø 2 (1979–1980) to Tromsø 3 (1986–1987) between the physical activity change groups from Tromsø 1 (1974) to Tromsø 2 (1979–1980) (p=0.68), (2) BMI change from Tromsø 3 (1986–1987) to Tromsø 4 (1994–1995) between the physical activity change groups from Tromsø 2 (1979–1980) to Tromsø 3 (1986–1987) (p=0.50), (3) BMI change from Tromsø 4 (1994–1995) to Tromsø 5 (2001) between the physical activity change groups from Tromsø 3 (1986–1987) to Tromsø 4 (1994–1995) (p=0.90), (4) BMI change from Tromsø 5 (2001) to Tromsø 6 (2007–2008) between the physical activity change groups from Tromsø 4 (1994–1995) to Tromsø 5 (2001) (p=0.98), (5) BMI change from Tromsø 6 (2007–2008) to Tromsø 7 (2015–2016) between the physical activity change groups from Tromsø 5 (2001) to Tromsø 6 (2007–2008) (p=20). Stratified analyses for the period-specific samples are presented in online supplemental tables 5–9). We observed no differences in BMI or weight change by occupational physical activity change in any strata analysis (all p≥0.13; except Tromsø 2–4 (1979–1995) sample, ≥4 years university education: p≤0.04; online supplemental table 8).
Discussion
In this large Norwegian population-based prospective study over four decades, we found no association between occupational physical activity changes and future BMI and weight changes.
Most previous longitudinal studies examined the association between baseline occupational physical activity and future BMI change,13–16 which do not account for the reciprocal temporal changes in physical activity and BMI.4 Two studies assessed changes in both occupational physical activity and BMI where one found lower occupational physical activity to be associated with weight gain,12 while one found no association.17 Without adjustment for previous physical activity levels, the direction of association and thus indication of causality, remains uncertain.4 Our study corroborates the findings of a recent study by Dobson et al,23 which regressed trajectories of self-reported BMI (ie, weight and height) on physical work exertion trajectories over nine time points in Canadian adults and showed no association between physical work exertion change and BMI trajectories, except for higher odds of being in a very obese trajectory (from 36 to 40 kg/m2 at follow-up) compared to a reference normal weight trajectory (22 to 24 kg/m2) with no higher odds of being in other BMI trajectories among those who decreased their physical work exertion compared with those who sustained low physical work exertion.23 Our study expands the work by Dobson et al 23 by using measured weight and height on both examinations and non-dichotomized BMI change as the outcome. Consequently, with higher accuracy in the outcome,24 the observed magnitudes in the association between occupational physical activity change and BMI change can be interpreted with higher confidence.4
As we did not adjust for energy intake due to unavailable data, our results may be influenced by residual confounding. Nevertheless, a previous study estimated that increasing physical activity energy expenditures of about 100 kcal/day would be sufficient for weight gain prevention at the population level,25 indicating that equivalent decreases would result in weight gain. This is similar to the estimated lower energy expenditure deriving from declines in occupational physical activity.3 As leisure time physical activity influences energy expenditure, one could hypothesise that occupational physical activity decline is only hazardous for those being physically inactive in leisure time. However, we observed no effect modification by leisure time physical activity changes.
It has been suggested that achieving energy balance and weight stability is easier at higher energy turnover.1 For example, energy intake increased by 500 kcal/day from the 1970s to 2000s in the USA, and 110–150 min of walking per day is needed to compensate for this increase.26 Consequently, as 150 min of walking per day is up to seven times higher than the current minimum recommendations for physical activity (150 min/week),27 and considering that one out of three adults in Western high-income countries fail to meet the recommendations,28 it is unlikely that the physical activity volume performed by the general population is sufficiently high to prevent weight gain.29
As occupational physical activity energy expenditure is dependent on activity duration, the effect of occupational physical activity on weight gain prevention may be influenced by whether individuals work full or part time. Thus, as we did not adjust for full-time and part-time work due to unavailable data, this may also have introduced residual confounding. However, these energy expenditure differences may in reality be small. For example, heavy manual labour workers are estimated to work at ~30% to 35% of maximal oxygen uptake over an 8 hours work day,30 which can be a sufficient volume to compensate the 500 kcal/day energy intake increase.26 However, few individuals in the Tromsø Study report heavy manual labour (~8% in 1979–1980,~2% in 2015–201610). In contrast, most occupational physical activities in the Tromsø Study changed from standing and walking to sitting,10 which is consistent with some cohorts.3 11 12 The energy expenditure difference while sitting compared with standing is estimated to be 54 kcal over 6 hours (ie, 72 kcal over 8 hours),31 which is unlikely to have any apparent effect on weight gain.
Some study cohorts in Southern Europe include a substantially larger proportion of heavy manual labour workers (Portugal, 37%32; Spain, Barcelona, 68%17); however, this is not consistent (Madrid, Spain: 2%,33 Italy: 8%34). Consequently, the generalisability of our findings may be limited to Northern/Central European8–10 and North American3 11 high-income countries. Studies examining weight gain prevention in heavy manual labour workers may be a future research target.
In our study, 741 (7%) participants are categorised as ‘Active to Inactive’, while 1315 (12%) participants were categorised as ‘Inactive to Active’ (table 3), indicating that more individuals increased their occupational physical activity level in our cohort. However, this is due to our crude categorisation of physical activity change; in our sensitivity analysis, 1315 (12%) are categorised as active but decreasing (rank 4 or 3→rank 3 or 2) (online supplemental table 4), where these are categorised as ‘Persistently Active’ in our main analysis (rank ≥2→rank ≥2) (table 3). Thus, the consistent pattern of declining occupational physical activity levels as in previous studies3 7–10 is confirmed in our study.
Our results indicate that occupational physical activity declines play a minor, if any, role in the observed population gain in BMI and weight. Consequently, public health initiatives aimed at weight gain prevention may have greater success by focusing on other aspects than occupational physical activity, for example, intake of energy dense food.2 26
The association between physical activity and BMI gain may also be reversed and/or bidirectional.4 High body weight appears causally associated with lower levels of physical activity when examining these associations using a Mendelian randomisation approach.35 However, intuitively, leisure-time physical activity is self-regulated while occupational physical activity is less controllable by the individual. Whether individuals regulate their occupational physical activity level depending on their BMI gain is questionable.
Strengths
First, as population gains in BMI have gradually increased over decades,36 the long follow-up time (~6 years) between each examination allowed us to examine whether occupational physical activity has contributed to BMI gain in this cohort.4 Second, by computing change in physical activity followed by change in BMI (accounting for previous physical activity level), we are able to interpret the direction of the association with more certainty.4 Third, by merging our period-specific samples to an overall cohort, we had higher power to examine multiple potential effect modifiers (table 3). For example, one warranted effect modification to be elucidated in associations between occupational physical activity and health outcomes is sex.37 Although we found differences in BMI gain by sex, we observed no effect modification of the associations by sex. Fourth, we used measured weight and height to calculate BMI as our outcome, which are more valid than self-reported weight and height,24 likely influenced by social desirability bias. Finally, the efforts to recruit representative samples and the high attendance in the Tromsø Study surveys indicate high representability of the population.18
Limitations
We categorised self-reported physical activity into crude groups, which have introduced misclassification, as described previously. Thus, we may have missed potential energy expenditure changes deriving from physical activity that could influence energy balance. However, crude groups of self-reported physical activity are valuable for categorisation of population levels of physical activity,38 and the SGPALS categorisations have previously shown associations with multiple health outcomes, suggesting predictive validity of the instrument.20 Moreover, our findings were unaltered when occupational physical activity change was categorised into six groups.
The recall and social desirability bias associated with self-reported physical activity likely results in over-reporting of physical activity levels,39 which is also demonstrated in office workers.40 Over-reporting of physical activity underestimates or overestimates the effect magnitude between physical activity and health outcomes.4 However, self-reported physical activity is currently the only instrument available in long-term ongoing cohort studies.4 Finally, as we did not adjust our models for energy intake and full-time/part-time work due to unavailable data, our results may be influenced by residual confounding.
Conclusion
We observed no association between changes in occupational physical activity and subsequent changes in BMI. Our findings do not support the hypothesis that occupational physical activity declines contributed to population gains in BMI over the past decades. Public health initiatives aimed at weight gain prevention may have greater success if focusing on other aspects than occupational physical activity.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Twitter @edvardhsagelv, @Ulf_Ekelund, @lailahopstock, @mariusfimland, @MorsethBente
Contributors EHS, BM, UE and LAH designed the study; EHS carried out statistical analyses; OL and TW provided statistical expertise; EHS drafted the manuscript; all authors interpreted the study results and contributed with manuscript revisions and approved the final version of the manuscript.
Funding The work of EHS is funded by Population Studies in the High North (Befolkningsundersøkelser i Nord: BiN), an internally funded project by UiT The Arctic University of Norway (no grant/award number). The remaining authors are funded by their respective positions/tenures.
Competing interests None declared.
Patient consent for publication Not required.
Ethics approval All participants in Tromsø 4–7 provided written informed consent, and the present study was approved by the Regional Ethics Committee for Medical Research (ref. 2016/758410).
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data may be obtained from a third party and are not publicly available. The data that support the findings of this study are available from the Tromsø Study, but restrictions apply to the availability of these data, which were used under licence for the current study, and so are not publicly available. The data can be made available from the Tromsø Study upon application to the Data and Publication Committee for the Tromsø Study (see www.tromsostudy.com).
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.