Objectives We investigate the independent and interacting long-term associations of occupational physical activity (OPA) and sport physical activity (SpPA) with the incidence of coronary heart disease (CHD) and cardiovascular diseases (CVD; CHD plus ischaemic stroke) in North Italian male workers.
Methods 3574 employed men aged 25–64 years, free of CVD at baseline, recruited in three population-based and one factory-based cohorts, were included in the analysis. The Baecke Questionnaire was used to assess OPA and SpPA in ‘minutes per week’ of moderate or vigorous PA. We estimated the associations between different domains of PA and the endpoints, adjusting for major CVD risk factors, using Cox models.
Results During a median follow-up of 14 years, 135 and 174 first CHD and CVD events, fatal and non-fatal, occurred. Compared with the intermediate OPA tertile, the HRs for CHD among low and high OPA workers were 1.66 (95% CI 1.06 to 2.59) and 1.18 (0.72 to 1.94), respectively (P value=0.07). Decreasing trends in CHD and CVD rates across increasing levels of SpPA were also found, with an HR for CVD of 0.68 (0.46 to 0.98) for intermediate/recommended SpPA compared with poor SpPA. We also found a statistically significant SpPA-OPA interaction, and the protective effect of SpPA was only found among sedentary workers, for both endpoints. Conversely, high OPA workers with intermediate/recommended SpPA levels had increased CHD and CVD rates compared with the poor SpPA category.
Conclusions Our results provide further evidence on the health paradox of OPA, with higher CVD rates among workers with intense PA at work. Moreover, the protective effect on CVDs of SpPA is prominent in sedentary workers, but it attenuates and even reverses in moderate and strenuous OPA workers.
- physical activity
- cardiovascular diseases
- preventive medicine
- longitudinal studies
Statistics from Altmetric.com
The 2015 Global Burden of Disease Study estimated that 15.2 million people died worldwide from coronary heart disease (CHD) and stroke.1 In Europe alone, CHD accounts for 1.8 million deaths annually.2 Physical activity (PA) is known as an important modifiable protective factor for the development of CHD and as such it has been included as an integral component of more updated cardiovascular prevention guidelines.3 4 Since 2010, the American Heart Association (AHA) recommends adults to perform at least 150 min per week of moderate-intensity activity or at least 75 min per week of vigorous-intensity activity, or at least 150 min per week of moderate and vigorous-intensity activities combined. These recommendations were integrated into the WHO global recommendations on PA,5 and allow classifying adults into three classes: poor, including those who practice none PA; intermediate, including those who practice PA without reaching the recommended minutes per week; and ideal, including those who practice PA for the recommended minutes per week.6
While most studies consistently indicated that high levels of leisure time physical activity (LTPA) decrease the risk of cardiovascular disease (CVD),7–22 the relationship between occupational PA (OPA) and CVD has shown contradictory results. A few studies have reported both intermediate and high levels of OPA as protective factors for CHD and stroke9 15 16 or intermediate levels of OPA only as a protective factor for CVD mortality and incidence.11 13 Conversely, other studies indicated that moderate and high levels of OPA increased the risk of major cardiovascular events.14 17 21 Finally, other studies did not find any association.12 19 Some studies have also estimated the interplay between LTPA and OPA on the CVD risk15 17 19–22 without univocal results, in particular concerning the effect of SpPA within groups exposed to intense OPA.
The aim of this study is to investigate the independent and interacting associations between OPA and SpPA categories, the latter based on AHA recommendations, and the incidence of coronary events in North Italian working men, over a follow-up of a median of 14 years.
Material and methods
This study is a pooled analysis of three population-based (the Multinational Monitoring of Trends and Determinants in Cardiovascular Disease (MONICA)—Brianza surveys II and III and the Pressioni Arteriose Monitorate e Loro Associazioni (PAMELA)—Study) and one factory-based (the Study of Employed in the Municipality of Milan (SEMM)—study) prospective cohorts of men aged 25 to 64 years, working as salaried employees (executive, white and blue-collar workers).23 24 No written informed consent signed by participants was required at time of recruitment. The study cohorts summed up to 3995 men who were free of CVD and currently employed at baseline. Women were not included in these analyses due to the low number of events during follow-up: 37 and 27 first coronary and ischaemic stroke events, respectively. A total of 421 men (10.5%) were excluded from the analysis due to missing information on CVD risk factors (n=175) or on PA (n=246). The final sample size consisted of 3574 male workers.
Baseline risk factors assessment and follow-up procedures
Information on adopted procedures to measure CVD risk factors at baseline and to detect endpoints has been published23–25; we summarise here major methodological aspects. The MONICA, PAMELA and SEMM baseline risk factors assessments were carried out by a unique team between 1989 and 1996. The overall participation rate was not <65%, with small variations among the four cohorts. Cardiovascular risk factors were collected adhering to the standardised procedures and quality standards of the WHO-MONICA Project (http://www.thl.fi/publications/monica/manual/index.htm). In the PAMELA and in the SEMM studies, risk factors were measured based on MONICA-like procedures (see online supplementary material).
Supplementary file 1
PA was measured using the Baecke Questionnaire,26 a known reliable and valid tool,27 to assess habitual PA. The questionnaire consists of 16 items from which an OPA score (eight items) and an SpPA index (four items) were derived. The OPA items investigated participant’s main occupation; his/her self-rating of the work’s vigour as compared with workers of the same age; the frequency of sitting, standing, walking, lifting and sweating at work; and the frequency of fatigue after work. The OPA score was categorised in tertiles, with cut-offs at 2.5 and 3.125 points (on a scale ranging from 1 to 5). Based on the method suggested by Atherosclerosis Risk in Communities (ARIC) investigators,10 SpPA items were converted in ‘minutes per week’ of moderate or vigorous PA based on metabolic equivalent of task,28 over the number of months annually performing the activities. The ‘minutes per week’ variable was then categorised into SpPA categories, according to AHA CVD health guidelines6: poor SpPA (0 min/week of activity); intermediate SpPA (1–149 min/week moderate or 1–74 min/week vigorous or 1–149 min/week moderate plus vigorous activity); and recommended SpPA (≥150 min/week moderate or ≥75 min/week vigorous or ≥150 moderate plus vigorous activity).
The same follow-up procedures were adopted for all cohorts, and people were followed until the end of 2008. The study endpoints were (1) a first acute coronary event as myocardial infarction, acute coronary syndrome or coronary revascularisation whichever occurred first; and (2) a first major cardiovascular event, including in addition a first ischaemic stroke or carotid endarterectomy. Loss-to-follow-up rate was <4%, with no differences across study cohorts.
We estimated the age-adjusted mean (prevalence) of major CVD risk factors by OPA and SpPA categories from linear (logistic) regression models, and tested differences among groups using the Wald χ2 test. A generalised logits model was used whenever the CVD risk factor had more than two levels. To estimate the association of OPA and SpPA categories with the incidence of cardiovascular events, we used Cox proportional hazards models, including both SpPA and OPA indexes and (1) age, cohort type (population vs factory-based), educational level (low, intermediate and high); and (2) age, cohort type, educational level, total and high-density lipoprotein (HDL) cholesterol, systolic blood pressure, current smoking, diabetes, alcohol intake and body mass index (BMI). Since the association between SpPA categories and CVD risk was linear in the ARIC study10 as well as in our preliminary analysis on event rates, we tested for statistical trend by additionally running each model with SpPA treated as a continuous variable. To investigate the effect of the length of follow-up time on CVD risk, we plotted Kaplan-Meier survival curves across OPA and SpPA categories. Finally, we explored the protective role of SpPA among workers engaged in different levels of OPA by adding relevant SpPA-OPA interaction terms in age and CVD risk factor-adjusted Cox regression models, and formally testing the presence of additive interaction on the log(HR) scale with a Wald χ2 test. The analyses were performed using the Statistical Analysis System (V.9.4 release; SAS Institute, Cary, North Carolina, USA), while figures were drawn using R (V.3.1.3, 2015, The R Foundation for Statistical Computing, Wien, Austria).
Among the included 3574 men, in a median follow-up time of 14 years (IQR 12.9–15.9 years) 135 first fatal and non-fatal CHD events, and 174 CVD events, occurred. Age-adjusted demographic characteristics and CVD risk factor levels at baseline for OPA and SpPA categories are shown in tables 1 and 2, respectively. Workers in high OPA were younger, less educated, with higher prevalence of current smokers and more than 50 gr per day alcohol consumers. People in the recommended SpPA category showed a better CVD risk factor profile, with lower mean levels of total cholesterol and BMI, higher median levels of HDL cholesterol and less current smokers.
The age, cohort and education-adjusted HRs and 95% CIs for the incidence of major CHD and CVD events across OPA and SpPA categories are shown in table 3. As these rates and HRs are mutually adjusted for the other type of PA, their estimates of the associations between two exposures (OPA and SpPA) with the endpoints are reciprocally independent. Keeping the intermediate OPA as the reference category, low OPA showed a significant higher age and educational class-adjusted CHD risk with HRs of 1.66 (95% CI 1.06 to 2.59) and high OPA an increased risk of 1.18 (0.72 to 1.94), but not reaching statistical significance. These results were not substantially modified when adjusting for other CHD risk factors, and were confirmed, with minor attenuations of HRs, for CVD as the endpoint.
Progressively decreased CHD risks were detected in the intermediate and the recommended SpPA categories, in comparison to the poor SpPA category (χ2 test for trend borderline statically significant at 0.08), with an overall reduction of 28% in CHD risk when considering the two categories combined. Here, the adjustments for other CVD risk factors reduced the associations, as most of them can be considered as direct or indirect mediators. The results were consistent for the CVD endpoint, with the notable observation that the reduced risk across SpPA categories turned statistically significant (χ2 test for trend, P=0.04) and the HR when considering the two intermediate and recommended categories combined was 0.68 (0.46 to 0.98).
The Kaplan-Meyer survival curves show that the reduced CHD (figure 1A) and CVD (figure 1B) risks of intermediate and recommended SpPA in comparison to the poor SpPA category started from the very beginning of the follow-up period, and increased progressively. Similarly, the excessive risk of the low OPA category in comparison to the two other categories was also clearly detectable for both endpoints (figure 1C,D), but with the notable exception that the risks for poor OPA started to reveal later in the follow-up period (3–5 years after baseline). This latency time may explain some of the discrepancies of the results in the literature, when the follow-up time is short.
Table 4 reports the risk factor-adjusted HRs (with 95% CI) for SpPA in different OPA categories, and the formal test for SpPA-OPA interaction on CHD and CVD risks. The analyses were carried out combining the intermediate and recommend SpPA categories to be able to reach sufficient statistical power to test for statistical interaction. In low OPA workers, characterised mainly by sedentary work activities, prominent protective effects of recommended or intermediate SpPA were seen, with consistent reduced HRs of 0.45 (95% CI 0.24 to 0.87) and 0.45 (95% CI 0.25 to 0.82) for CHD and CVD, respectively, for the two active SpPA categories combined. These decreased HRs for higher levels of SpPA disappeared in the intermediate OPA category, and even reversed in the high OPA categories with values of 1.84 (0.88 to 3.87) for CHD and 1.66 (0.87 to 3.14) for CVD. The heterogeneity tests and the interaction tests were statically significant at an α level of 0.05, for both endpoints. The detrimental effect of SpPA among high OPA workers was even more pronounced when recommended SpPA levels were met (online supplementary table S1) and when taking into account commuting PA (walking or cycling to work for more than 30 min/day; online supplementary table S2).
We found long-lasting (over a median 14 years of follow-up) and prominent protective effects of performing SpPA on CHD and CVD risks in adult North Italian working men. These reductions, when intermediate and recommended levels of sport PA were combined, mounted to 28% for CHD and 32% for CVD, the latter being statistically significant. Similar results in men were reported by most of the previous primary and secondary studies for LTPA.10 14 17 18 21 22 For women, the results are less consistent, with a reduced CHD risk in a Finnish study,15 and with mostly null effects reported by another study.18 Furthermore, the protective effects in our study were only partially reduced when adjusting for behavioural and metabolic risk factors. In our data, lower prevalence of BMI and current cigarette smokers and means of total cholesterol with higher HDL cholesterol levels were found among those performing recommended levels of SpPA, indicating potential mediation effects of these factors. Other studies adjusted for different sets of demographic characteristics and CVD risk factors,14 and this might explain in part the heterogeneity of their results.
We focused on three categories of SpPA, based on AHA recommended levels,3 4 which can be derived from the Baecke Questionnaire, and are well characterised in terms of frequency, minutes per week and intensity. Most of other papers used LTPA, which includes a combination of sport and recreational (walking and bicycling during leisure time and for going to work) domains of PA.15 18 21 A few more recent studies better characterised PA in terms of frequency (times per week) and duration (minutes per session).17 22 The variety of LTPA definitions may explain in part the differences in estimating the size of the effects, as PA levels may have different meanings in different geographical, gender and age groups,10 and a non-quantitative scoring may not grasp such differences. The Baecke Questionnaire includes items exploring the recreational domain of LTPA. In our data, these items showed a poor internal correlation (Cronbach α coefficient <0.4), modest predictive ability (dissimilar single item–endpoints associations) and reduced construct validity (only the item ‘During leisure time I walk’ loads the factor with an acceptable coefficient of 0.72), indicating that walking was the single, main determinant of the LPTA score (data not shown). Due to this poor performance of these items, we excluded recreational LTPA from our analysis.
Our study found a statistically significant 66% increase in CHD incidence rates among workers with low levels of OPA, when compared with the intermediate OPA category and adjusting for SpPA as well as for other major CVD risk factors. We also estimated a 18% increase in both CHD and CVD risks among workers with higher levels of OPA, neither statistically significant. A recent meta-analysis12 reported similar results for the higher but not for the lower levels of OPA. The discrepancies may be explained by the inclusion of different proportion of sedentary workers in the examined samples, as well as by the conjunct consideration of sport/LTPA (as discussed later in this paper). The higher CVD risk in sedentary people (mainly due to the work component) was also found in other studies.9 11 On the other hand, a Danish study21 did not find a higher risk in ischaemic heart disease mortality in sedentary work activities, but the different endpoint and the effect of Myocardial Infarction (MI) survival should be acknowledged. A Swedish study19 did not find an increase in risk of myocardial infarction in high OPA. Again, reasons for these discrepancies may be due to the different methods adopted for the assessment of OPA,10 14 and also to the differences in the length of follow-up, as can be seen in figure 1, where OPA differences in risk start to show up after 3–5 years of follow-up.
In our study, the protective effect of AHA recommended and intermediate SpPA levels was strongly dependent of OPA level, with a prominent risk reduction among sedentary workers (55%), and an opposite relationship for workers with higher level of OPA increasing the risk up to 66% for CVD and 84% for CHD. These findings suggest that the strongest effect of SpPA can be detected among people involved in mainly sedentary jobs, in particular when recommended SpPA levels are met and when commuting PA is limited (online supplementary tables S1 and S2). Since these jobs mainly characterise white collars, the residual confounding of social class should be acknowledged, although we adjusted for the educational level. On the other hand, performing heavy physical working tasks and engaging in high levels of sport activities during leisure time may lead to premature exhaustion and fatigue, and generate an overloading of the cardiovascular system, which over years with insufficient recovery can lead to CVD.29 Our online supplementary material findings on workers meeting recommended SpPA levels and with intense daily commuting PA strengthen our confidence to the main results, and provide further support to the overloading hypothesis. Moreover, this opposite effect of high SpPA among people with high levels of OPA may be a further explanation of the paradoxical adverse effect of OPA on CVD risk. In table 5, we report the effect of high versus low LTPA/SpPA in different OPA categories from published studies which addressed the interplay between OPA and LTPA/SpPA in working men. One study19 was not included as published results do not allow deriving comparable calculations. A Belgian study,22 which adopted similar OPA assessment tools and endpoint (CHD incidence), did find the same results. Our study just extends these findings on a longer follow-up period. A Finnish15 and an Israeli17 study found a prominent protective effect of LTPA in low OPA categories only. Conversely, another Danish study, which used fatal ischaemic heart diseases only as the endpoint, found protective effects of LTPA over all OPA categories.21 Finally, Wang et al20 measured OPA adopting the relative aerobic strain or workload method, which takes into account both the absolute energy expenditure of the tasks and the workers’ individual aerobic capacity. In this study, the effect of LTPA was modest both in low and high OPA categories. Moreover, a recent report has shown that the beneficial effect of LTPA increasing heart rate variability reduces when combined with higher OPA levels.30 Anyhow, more collaborative studies are needed to investigate the effects of LTPA/SpPA among workers with high OPA, adopting common and standardised methodologies.
Our study presents some limitations. First, PA was calculated using a self-administered questionnaire, once on the entire follow-up period. We can thus assume some bias due to misclassifications of the exposures: in general, HRs of higher levels of both sport PA and OPA may be underestimated as people getting old moderate the former and are reallocated at work and reduce strenuous tasks. Although we adjusted for major risk factors, we do not have data on nutrition as a potential confounder. Another limitation is the focus on men only, as the low number of events among women determined very unstable effect estimates, in particular for the OPA-SpPA interaction analyses. As mentioned, results in men and women can be heterogeneous because of both different levels of PA and different associations between PA domains and the CVD endpoint, so combining the two gender groups may be misleading and can produce confused findings. Finally, our data do not allow to further investigate the separate effects of walking and lifting at work on the CVD risk, as the large majority of high OPA workers (67%) reported both frequent walking and lifting.
In conclusion, our study conducted on North Italian male workers suggests a prominent protective effect of AHA recommended and intermediate sport PA levels on CVDs among sedentary workers, and an opposite effect in people involved in physically demanding working tasks. The second finding particularly needs further investigations in other cohorts to assess consistency of the results. Due to the increasing number of sedentary workers in postindustrialised countries, interventions to promote SpPA in working populations could be an important factor to reduce the burden of CVDs. In addition, future CVD prevention guidelines4 should consider tailoring the recommendations on the amount of PA during leisure time to the levels of PA at work.
What is already known on this subject?
While consistent evidence indicates that elevated and intermediate levels of leisure time physical activity decrease the risk of cardiovascular diseases (CVD), the relationship between occupational physical activity (OPA) levels and CVD has showed contrasting results.
What might this study add?
Our study, conducted on North Italian male working cohorts, confirms the protective effect of both American Heart Association recommended and intermediate levels of sport activities, with a reduced CVD risk of 32%, as well as a paradoxical effect of higher CVD risks in low and high levels of OPA when compared with the intermediate level.
The protective effect of recommended and intermediate sport physical activity (SpPA) levels strongly depends on OPA levels, with a prominent risk reduction among sedentary workers (55%), and an opposite relationship for workers with higher level of OPA increasing the risk up to 66% for CVD and 84% for coronary heart disease.
How might this impact on clinical practice?
The increasing number of sedentary workers in postindustrialised societies asks for specific interventions to promote SpPA even at the workplace. If our results were confirmed in larger studies, the CVD prevention recommendation of SpPA ought to consider the level of OPA.
Contributors MMF originally conceived the research question. MMF, MR and GV drafted the manuscript and MMF and GV prepared the revisions. MMF is the PI of the SEMM study. MMF and GC are co-PIs of the MONICA Brianza cohorts and GG is PI for the PAMELA study. GV was responsible for the statistical analyses. AH, MC, EC, RB, GG and GC critically revised the text, with particular reference to the introduction and to the discussion sections.
Funding This work and the latest activities of the Brianza, PAMELA and SEMM cohorts have been funded by the Health Administration of the Lombardia Region (grants nos. 17155/2004 and 10800/2009). The follow-up was partially supported with grants from the Italian Ministry of Health (grant 2012/597) and it was carried out in collaboration with the Department for Cardiovascular Dysmetabolic and Ageing –Associated Diseases of the Istituto Superiore di Sanità in Rome.
Competing interests None declared.
Ethics approval The MONICA-Brianza and the PAMELA studies were approved by the Ethical Committee of the University Hospital of Monza (Italy). The SEMM study was approved by the Ethical Committee of the University Hospital of Varese (Italy).
Provenance and peer review Not commissioned; externally peer reviewed.
Collaborators The Cohorts Collaborative Study in Northern Italy (CCSNI) Research Group: MONICA-Brianza cohorts: PIs: Professor Marco Mario Ferrario, Professor Giancarlo Cesana; Baseline data collection: Professor R Sega, Dr M T Gussoni, Dr F Duzioni, Sig L Bolognesi; Data quality assessment: Dr G Veronesi, Dr C Fornari, Dr L Bertù. PAMELA study: PI: Professor Guido Grassi. Professor R Sega and Professor G Mancia (former PIs). Baseline data collection and follow-up activities: Professor R Sega, Dr M Bombelli. Data quality assessment: Dr R Facchetti. SEMM Study: PI: Professor Marco M Ferrario. Baseline data collection: Dr A Vivaldi. Data quality assessment: Dr G Veronesi, Dr C Fornari, Dr L Bertù. Laboratory and biobanking: Professor P Brambilla, Dr S Signorini.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.