Job strain in relation to body mass index: pooled analysis of 160 000 adults from 13 cohort studies

Nyberg ST, Heikkilä K, Fransson EI, Alfredsson L, De Bacquer D, Bjorner JB, Bonenfant S, Borritz M, Burr H, Casini A, Clays E, Dragano N, Erbel R, Geuskens GA, Goldberg M, Hooftman WE, Houtman IL, Jöckel K-H, Kittel F, Knutsson A, Koskenvuo M, Leineweber C, Lunau T, Madsen IEH, Magnusson Hanson LL, Marmot MG, Nielsen ML, Nordin M, Oksanen T, Pentti J, Rugulies R, Siegrist J, Suominen S, Vahtera J, Virtanen M, Westerholm P, Westerlund H, Zins M, Ferrie JE, Theorell T, Steptoe A, Hamer M, Singh-Manoux A, Batty GD, Kivimäki M, for the IPD-Work Consortium (Finnish Institute of Occupational Health, Helsinki, Finland; Institute of Environmental Medicine, Karolinska Institutet, Stockholm, Sweden; Jönköping University, Jönköping, Sweden; Ghent University, Ghent, Belgium; National Research Centre for the Working Environment, Copenhagen, Denmark; Versailles-Saint Quentin University, Versailles, France; Centre for Research in Epidemiology and Population Health, Villejuif, France; Bispebjerg University Hospital, Copenhagen, Denmark; Centre for Maritime Health and Safety, Esbjerg, Denmark; Université Libre de Bruxelles, Brussels, Belgium; University Duisburg-Essen, Essen, Germany; West-German Heart Center Essen, University Duisburg-Essen, Essen, Germany; TNO, Hoofddorp, the Netherlands; Mid Sweden University, Sundsvall, Sweden; University of Helsinki, Helsinki, Finland; Stockholm University, Stockholm, Sweden; University College London, London, UK; Bispebjerg University Hospital, Copenhagen, Denmark; Umeå University, Umeå, Sweden; Finnish Institute of Occupational Health, Turku, Finland; Harvard School of Public Health, Boston, MA, USA; University of Copenhagen, Copenhagen, Denmark; University of Düsseldorf, Düsseldorf, Germany; University of Turku, Turku; Folkhälsan Research Center, Helsinki; Turku University Hospital, Turku, Finland; Uppsala University, Uppsala; Karolinska Institutet, Stockholm, Sweden; University of Bristol, Bristol; University of Edinburgh, Edinburgh, UK; and University of Helsinki, Helsinki, Finland). Job strain in relation to body mass index: pooled analysis of 160 000 adults from 13 cohort studies. J Intern Med 2012; 272: 65–73. Background Evidence of an association between job strain and obesity is inconsistent, mostly limited to small-scale studies, and does not distinguish between categories of underweight or obesity subclasses. Objectives To examine the association between job strain and body mass index (BMI) in a large adult population. Methods We performed a pooled cross-sectional analysis based on individual-level data from 13 European studies resulting in a total of 161 746 participants (49% men, mean age, 43.7 years). Longitudinal analysis with a median follow-up of 4 years was possible for four cohort studies (n = 42 222). Results A total of 86 429 participants were of normal weight (BMI 18.5–24.9 kg m−2), 2149 were underweight (BMI < 18.5 kg m−2), 56 572 overweight (BMI 25.0–29.9 kg m−2) and 13 523 class I (BMI 30–34.9 kg m−2) and 3073 classes II/III (BMI ≥ 35 kg m−2) obese. In addition, 27 010 (17%) participants reported job strain. In cross-sectional analyses, we found increased odds of job strain amongst underweight [odds ratio 1.12, 95% confidence interval (CI) 1.00–1.25], obese class I (odds ratio 1.07, 95% CI 1.02–1.12) and obese classes II/III participants (odds ratio 1.14, 95% CI 1.01–1.28) as compared with participants of normal weight. In longitudinal analysis, both weight gain and weight loss were related to the onset of job strain during follow-up. Conclusions In an analysis of European data, we found both weight gain and weight loss to be associated with the onset of job strain, consistent with a ‘U’-shaped cross-sectional association between job strain and BMI. These associations were relatively modest; therefore, it is unlikely that intervention to reduce job strain would be effective in combating obesity at a population level.

Objectives. To examine the association between job strain and body mass index (BMI) in a large adult population.
Methods. We performed a pooled cross-sectional analysis based on individual-level data from 13 European studies resulting in a total of 161 746 participants (49% men, mean age, 43.7 years). Longitudinal analysis with a median follow-up of 4 years was possible for four cohort studies (n = 42 222).
Results. A total of 86 429 participants were of normal weight (BMI 18.5-24.9 kg m )2 ), 2149 were underweight (BMI < 18.5 kg m )2 ), 56 572 overweight (BMI 25.0-29.9 kg m )2 ) and 13 523 class I (BMI 30-34.9 kg m )2 ) and 3073 classes II ⁄ III (BMI ‡ 35 kg m )2 ) obese. In addition, 27 010 (17%) participants reported job strain. In cross-sectional analyses, we found increased odds of job strain Introduction Obesity and job strain (i.e. stress at work) are major public health issues in modern societies, potentially contributing to a range of health-related outcomes, such as reduced quality of life, disability, cardiovascular and cerebrovascular diseases and depression [1][2][3]. According to recent European Union estimates, stress is cited as a factor in half of all lost working days and thus represents a substantial cost in terms of human distress and impaired economic performance [4]. There may be a link between job strain and body mass index (BMI) [5][6][7][8][9][10][11][12] -the most commonly utilized measure of adiposity -as stress might contribute to an unhealthy lifestyle [5], such as physical inactivity [6] and a poor diet [7], which in turn could induce weight gain. Other mechanisms are also plausible. Conversely, psychosocial stress may reduce appetite leading to weight loss [8][9][10]. In addition to stress being a risk factor for weight change, there is a suggestion that this relationship might be bi-directional. Obesity, for instance, may reduce work capacity [11], increasing the risk of feelings of stress (the reverse causation hypothesis). Finally, given its association with both increased weight and exposure to stressful work conditions [12], it is also likely that socio-economic disadvantage may have an important role in these relationships (the common cause hypothesis).
To date, empirical evidence for an association between job strain (or other forms of work stress) and BMI has been inconsistent, with findings of a posi-tive (more stress, higher BMI) [8,[13][14][15][16][17], an inverse (more stress, lower BMI) [18,19] or no relation [20,21]. Small sample sizes in most of these studies may have contributed to the mixed results. This low study power has also led to an inability to distinguish between categories of underweight or different classes of obesity. To enable more precise characterization of the association between job strain and BMI than in previous studies, we pooled data from 13 independent cohort studies, resulting in an individual-level meta-analysis of 161 746 men and women.

Study population and design
This study is part of the individual-participantdata meta-analysis in working populations (IPD-Work) consortium of European cohort studies. A collaboration of five studies was established at a workshop in London, UK on 8 November 2008, since when a further eight cohort studies have been included.
In this study, we pooled data from 13 prospective cohort studies (see Table 1 for full names): from Belgium (Belstress), Denmark (DWECS, IPAW, PUMA), Finland (FPS, HeSSup), France (Gazel), Germany (HNR), the Netherlands (POLS), Sweden (SLOSH, WOLF-N, WOLF-S) and the UK (Whitehall II). Details of the design, recruitment, measurements and ethical approval of the participating studies are presented in the Appendix S1. Participants with complete data on BMI, job strain, sex and age were included in these analyses, yielding a sample of 78 487 employed men and 83 259 employed women (mean age 43.7 years at study entry). Characteristics of these studies and the participants are shown in Table 1.

Assessment of BMI
Body mass index was calculated using the usual formula: weight in kilograms divided by height in metres squared. Participants with missing values for weight or height were excluded (n = 2220; 1.4% of all participants). To avoid a few potentially unreliable measurements unduly affecting the results, participants with BMI values <15 or >50 kg m )2 were excluded from the analysis (n = 100; 0.1%). We classified BMI into five categories according to World Health Organization (WHO) recommendations [22]. Participants with a BMI < 18.5 kg m )2 were categorized as underweight, those with a BMI between 18.5 and <25 kg m )2 were classified as normal weight, those with a BMI between 25 and <30 kg m )2 as overweight and, following the WHO international classification of adult obesity [22], we included two subcategories of obesity: class I (BMI 30 to <35 kg m )2 ) and class II and III combined (BMI ‡ 35 kg m )2 ).

Definition of job strain
According to the job strain model -the most widely tested model of work stress -job strain arises when an employee simultaneously has high psychological job demands and a low level of control over work [23]. In the included studies, job strain was assessed using participant-completed questionnaires. All questions within the job demand and job control scales required responses in Likert-type formats. Mean response scores for job demand items and for job control items were calculated for each participant. An unfavourable (high) level of job demand was denoted by a score above the study-specific median, whereas an unfavourable (low) level of job control was defined as a score below the studyspecific median score. We defined job strain as this combination of the two variables (i.e. a high level of job demand and a low level of job control). All other combinations of job demand and control, including levels equal to the median values, were defined as nonjob strain. Participants with missing data on more than half of the items of job demand or job control were excluded (n = 1714; 1.1%).

Covariates
Covariates included in the analyses were age, sex, socio-economic status (SES; high, intermediate and low) and smoking status (current smoker versus nonsmoker). Participants with missing values for either age or sex were excluded from all analyses (n = 367; 0.2%). A more detailed description of the assessment of covariates is presented in the Appendix S1.

Statistical analyses
We examined individual-level data from nine studies. For a further four, we provided syntax and instructions for statistical analysis, as the study investigators chose to carry out their own analyses. One-stage and two-stage meta-analyses of individual participant data [24][25][26] were performed. In the cross-sectional analysis, we used two-stage meta-analysis to include all cohort studies irrespective of whether individual-level or aggregate data were available from the study.
For each study, effect estimates and their standard errors were estimated using logistic regression (the first stage); these study-specific results were then pooled using random-effects meta-analysis (the second stage) [27]. We calculated summary odds ratios and their 95% confidence intervals (CIs) for job strain in individuals who were categorized as underweight, overweight or obese (classes I and II ⁄ III), compared with individuals of normal weight. We adjusted the odds ratios for sex, age, SES and smoking. To test the associations between BMI and both job demand and job control, we computed summary mean difference in demand and control scores between BMI categories using linear regression. Heterogeneity amongst study-specific estimates was assessed using the I 2 statistic [28]. In a sensitivity analysis, we ran the analyses separately for studies in which individuallevel data were available for pooled analysis. Additionally, to examine measurement method as a source of heterogeneity, we ran these analyses separately for studies with measured height and weight and for those with self-reported values.
In order to examine subgroup differences and longitudinal associations, we used a one-stage meta-analysis pooling all available individual-level data into one data set. We tested for possible interactions of BMI category, sex and age group (>50 vs. £50 years) by including an interaction term (BMI*covariate) in the model using a mixed-effects logistic regression model with study as the random effect. In four studies (Belstress, FPS, HeSSup and Whitehall II), BMI and job strain components had been re-measured approximately 4 years apart thus allowing us to examine the longitudinal associations between job strain and BMI categories in this subgroup of cohorts.
To define job strain at follow-up, we used the same study-specific cut-off points that were used at baseline. These studies allowed us to examine a series of subsidiary questions. (i) Does exposure to job strain predict obesity amongst nonobese participants, and is an association with obesity stronger for those with repeated exposure to job strain (test of a dose-response association)? (ii) Are both weight gain (change from nonobese to obese between baseline and followup) and weight loss (change from obese to nonobese during the same period) related to the onset of job strain at follow-up? (iii) Does obesity at baseline predict the onset of job strain at follow-up (test of reverse causation)? (iv) Does SES at baseline predict obesity and job strain at follow-up and are the associations between job strain and obesity attenuated in a stratified analysis within the three strata of SES (test of the common cause hypothesis)?
Models were fitted with PROC GENMOD, PROC GLIMMIX and PROC MIXED in sas 9 or spss 17. The meta-analysis was conducted using r (version 2.11; library Meta, http://www.r-project.org). More details about statistical analysis can be found in the Appendix S1.  Table 1.

Amongst
Job demand, job control and obesity Figure S1 shows a forest plot of the mean differences in job demand score in each BMI category relative to the normal weight group. In an age-and sex-adjusted model (model 1), no association was observed between BMI category and job demand score. After further control for SES (model 2), there was some suggestion of a dose-related link, with higher job demand being associated with a higher risk of obesity, although all point estimates included zero. This positive relation was also seen in the longitudinal analysis of incident obesity (age-, sex-and SES-adjusted odds ratio for top versus bottom quintile 1.14, 95% CI 0.99-1.32, Table S1). Figure S2 shows a corresponding forest plot for job control and BMI categories. In age-and sex-adjusted analyses (model 1), job control was slightly lower amongst underweight, overweight and obese participants compared with their normal-weight counterparts. However, after adding SES to the multivariable model, with the exception of the underweight group, all these differences were statistically nonsignificant (model 2).
Job strain and obesity Figure 1 shows a forest plot of the random-effect summary odds ratios for job strain in each BMI category (study-specific results are provided in Fig. S3-S5). In an age-and sex-adjusted model (model 1), there was a suggestion of a 'U'-shaped relation: the greatest risk of job strain was evident in the underweight and obese groups, whilst the risk was lowest in the normalweight group. Thus, the odds ratio for job strain was 1. , but values remained statistically significant for the two obesity categories. Further adjustment for smoking had essentially no effect on these estimates.

Longitudinal associations
Amongst the participants who were nonobese at baseline, low versus high SES at baseline was related to the risk of subsequent obesity, with an age-and sex-adjusted odds ratio of 1.54 (95% CI 1.35-1.76). Table 2 shows the longitudinal associations between job strain and obesity at follow-up in this population. These analyses are based on four cohort studies with a median (interquartile range) follow-up of 4 (1) years. Job strain at baseline alone or at both baseline and follow-up was not associated with obesity at follow-up. Similarly, in further analyses, change in BMI during follow-up did not differ between initially nonobese participants with and without job strain at baseline (age-, sex-and SES-adjusted mean difference in BMI change )0.02, 95% CI )0.06 to 0.02 kg m )2 , P = 0.46), or between those with and without job strain at baseline and follow-up (mean difference )0.04, 95% CI )0.10 to 0.02 kg m )2 , P = 0.22) (data not shown). However, new exposure to job strain at follow-up was associated with becoming obese at follow-up (odds ratio compared with no job strain at baseline and follow-up 1.18, 95% CI 1.02-1.36) ( Table 2). When we examined this relationship within the three strata of SES, the results were essentially unchanged. Table 3 shows the longitudinal analysis relating BMI with job strain at follow-up amongst participants without job strain at baseline. Low SES at baseline was a strong predictor of job strain at follow-up (odds ratio 2.93, 95% CI 2.64-3.24), but baseline BMI categories were not associated with subsequent job strain (no support for the reverse causation hypothesis). Becoming obese was associated with a raised risk of job strain at follow-up (odds ratio 1.18, 95% CI 1.02-1.36). This was also evident within all strata of SES although all CI values included unity. Change from obese to nonobese was also associated with an increased odds of job strain at follow-up (odds ratio 1.31, 95% CI 1.03-1.68 compared with nonobese at baseline and follow-up), a finding replicated at low and intermediate levels of SES, although these analyses were restricted by low numbers (only five cases of incident job strain in the high-SES group).

Sensitivity analyses
We found no statistical evidence to suggest that the cross-sectional association between job strain and obesity varied between particip ants younger and older than 50 years of age (P for interaction = 0.36) or between men and women (P for interaction = 0.35). Furthermore, the results described earlier remained largely unchanged after exclusion of the four studies that did not share individual-level data or when the analyses were performed separately for clinically measured versus self-reported BMI. Adjustment for the length of follow-up had essentially no effect on the longitudinal association estimates.

Discussion
The aim of this analysis of pooled data from approximately 160 000 adults in 13 European studies was to describe the association between job strain and BMI  in greater detail than has previously been possible.
The results show a 'U'-shaped association between the two factors, with the proportion of employees with job strain highest in the underweight and obese groups, and the lowest risk of job strain in individuals of normal weight. From two narrative reviews, based on studies with smaller numbers, the authors found no consistent cross-sectional association between work stress and BMI [29,30]. However, in these analyses, the stressed and nonstressed participants were compared in terms of mean BMI, making it difficult to detect higher levels of stress amongst both underweight and obese individuals.
Our longitudinal analysis shows that changes in job strain and BMI category tend to co-occur. First, we found that change from no job strain at study baseline to job strain at follow-up is correlated with change from nonobese at baseline to obese at followup, a finding also apparent when we stratified analysis for each socio-economic group. Second, we found that change from no job strain at baseline to job strain at follow-up was also associated with reduction in weight (from obese to nonobese), again largely independently of SES. Thus, both weight gain and weight loss were associated with the onset of job strain, a finding which is consistent with the 'U'-shaped crosssectional association between job strain and BMI category.
We found little direct evidence to suggest that job strain is a causal risk factor for weight gain. First, the association was substantially reduced after adjustment for SES; second, baseline job strain did not predict change in BMI or the risk for obesity in longitudinal analysis. Finally, repeated measurements of job strain provided no evidence of dose-response associations between job strain and BMI or obesity. These findings are in agreement with those of previous studies. In a study-based meta-analysis of 8514 participants, Wardle and colleagues found no clear evidence for a longitudinal association between job strain and BMI (correlation coefficient 0.014, 95% CI )0.002 to 0.031, P = 0.09) [31]. This is in agreement with data from Japanese [32], Swiss [33], Swedish [34] and Finnish [35] studies which reported no association between job strain ⁄ work stress and change in adiposity. It has been suggested that the effect of job strain on change in BMI might differ between subgroups of individuals [10,34], or the relationship may be limited to waist circumference [32]. However, BMI and waist circumference are strongly correlated [36] implying that a predictive association should also be seen for BMI if job strain was a strong predictor of waist circumference. Some studies have examined associations between other indicators of work stress (e.g. job insecurity or iso-strain) and weight change but with inconsistent findings [8,15,16,33,37].
Considering the reverse causation hypothesis, there was no evidence to suggest that obesity confers an increased risk of job strain. The fact that neither a direct causal effect nor the reverse causation hypothesis were supported by the results of our longitudinal analyses raises the possibility that common causes might underlie the apparent association between the onset of job strain and weight change. In cross-sectional age-and sex-adjusted analyses, the excess odds of job strain were approximately 20% in obese class I individuals and 30% for those in obese classes II ⁄ III; however, adjustment for SES attenuated these estimates to 7% and 14%, respectively. This attenuation suggests that socio-economic adversity is likely to at least partially explain the association between job strain and obesity. In the longitudinal analyses, similar associations between the onset of job strain and weight change were observed within socio-economic groups which means that these associations are unlikely to be solely explained by SES, but other, yet unknown, factors may also be involved. Further research is needed to confirm this. It may be that adverse life events and the onset of psychiatric disorders, particularly depressive symptoms, contribute to the association between the onset of job strain and weight gain, as these factors are known to affect weight control and reporting of job strain [38]. Previous research suggests a robust association between nonintentional weight loss, being underweight and increased mortality [39,40], which is largely attributable to a pre-existing physical illness. This explanation might also apply in the present study with preexisting physical morbidity potentially underlying the associations between weight loss, being underweight and job strain.
Our study has several important strengths. First, to our knowledge, this is the first study in which the association between BMI and job strain was studied across the entire BMI distribution; that is, including underweight individuals as well as two subcategories of obese individuals. Secondly, the analysis covers multiple study populations from several countries, increasing the generalizability of the findings. Given that the sample size was larger than in any prior study, the likelihood of random error influencing our results is also lower than in previous studies. Thirdly, we defined work stress based on the job strain model, which is the most widely used though not the only conceptualization in this area of research. However, some limitations should also be noted. Apart from SES and smoking, we did not examine the role of potential mediating or confounding factors. Despite data harmonization, variation in the assessment of job strain and SES between studies may have contributed to inaccuracy of the estimates. Data harmonization also meant that the measures of job strain and SES used in this study might not be optimally adjusted for the specific contexts of each participating study, potentially contributing to underestimation of the associations. On the other hand, using studyspecific measurements, as in previous analyses, may introduce information bias and overestimation of the associations.
In summary, data from 13 European cohort studies show a cross-sectional 'U'-shaped association between job strain and being either obese or underweight, and corresponding longitudinal associations between the onset of job strain and both weight gain and weight loss. As these associations were relatively modest in terms of absolute effect size and not necessarily causal, our data do not suggest that interventions to reduce job strain would be effective in combating obesity at a population level. However, early screening for job strain and obesity in the workplace may inform appropriate treatment strategies or lifestyle changes to prevent adverse health outcomes associated with these conditions, such as work disability and depressive disorders.

Conflict of interest statement
No conflict of interest was declared.

Supporting Information
Additional Supplementary Information may be found in the online version of this article: Figure S1. Summary estimates for the association between BMI categories and job demand score. Figure S2. Summary estimates for the association between BMI categories and job control score. Figure S3. Study-specific meta-analysis of job strain amongst underweight (BMI < 18.5 kg m )2 ) versus normal weight (BMI 18.5-24.9 kg m )2 ) adults. Figure S4. Study-specific meta-analysis of job strain amongst overweight (BMI 25.0-29.9 kg m )2 ) versus normal weight (BMI 18.5-24.9 kg m )2 ) adults. Figure S5. Study-specific meta-analysis of job strain amongst obese (BMI ‡ 30 kg m )2 ) versus normal weight (BMI 18.5-24.9 kg m )2 ) adults. Table S1. Longitudinal associations between job demand and job control quintiles and incident obesity in four studies with repeat data.
Appendix S1. Study design, recruitment of participants and measurements in the 13 participating studies.
Please note: Wiley-Blackwell are not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author of the article.