The role of physical activity in the development of first cardiovascular disease event: a tree-structured survival analysis of the Danish ADDITION-PRO cohort
Cardiovascular Diabetology volume 17, Article number: 126 (2018)
Ambiguity exists in relation to the role of physical activity (PA) for cardiovascular disease (CVD) risk reduction. We examined the interplay between PA dimensions and more conventional CVD risk factors to assess which PA dimensions were associated with the first CVD event and whether subgroup differences exist.
A total of 1449 individuals [median age 65.8 (IQR: 61.2, 70.7) years] with low to high risk of type 2 diabetes and free from CVD from the Danish ADDITION-PRO study were included for survival analysis. PA was measured by individually calibrated heart rate and movement sensing for 7 consecutive days. The associations of different PA dimensions (PA energy expenditure, time spent in light-, moderate- and vigorous intensity PA), sedentary time and other conventional CVD risk factors with the first CVD event were examined by tree-structured survival analysis. Baseline information was linked to data on the first CVD event (ischemic heart disease, ischemic stroke, heart failure, atrial flutter/fibrillation and atherosclerotic disease) and mortality obtained from Danish registers.
During a median follow-up time of 5.5 (IQR: 5.1–6.1) years, a total of 201 individuals (13.9%) developed CVD. Overall CVD incidence rate was 2.6/100 person-years. PA energy expenditure above 43 kJ/kg/day was associated with lower rates of CVD events among participants ≤ 70 years and with HbA1c ≤ 5.7% (39 mmol/mol), systolic blood pressure ≤ 156 mmHg and albumin creatinine ratio ≤ 70 (incidence rates 0.0–0.8/100 person-years).
Any type of PA resulting in increased PA energy expenditure may over time be the best prevention strategy to uphold reduced risk of CVD.
As in the general population, cardiovascular disease (CVD) is a major contributor to morbidity and mortality in individuals with pre-diabetes and type 2 diabetes. However, in these high-risk individuals the risk of developing CVD is greatly increased [1, 2].
Increasing physical activity (PA) levels have shown a substantial favourable impact on type 2 diabetes progression in high-risk individuals [3, 4]. However, the evidence regarding the role of PA in the prevention of CVD in different subgroups is less clear. Observational studies have demonstrated that greater amounts of PA are associated with lower rates of incident and recurrent CVD [5, 6]. However, this is not unambiguously supported by results from large randomized controlled trials of lifestyle interventions in subgroups at high CVD risk . Despite significant improvements of CVD risk factors after increases in PA, these do not seem to convey any appreciable reduction in CVD event rates . The reasons for this inconsistency are not yet fully understood. In contrast, one randomized controlled study found that PA reduced the incidence of CVD among individuals with pre-diabetes after 23 years .
Furthermore, the majority of observational studies of the role of PA and the development of CVD have looked only at single measures of PA. Knowledge on how different dimension of PA such as intensity, frequency, duration, and total energy expenditure are associated to CVD risk is therefore limited. In addition, these studies have mainly been conducted in healthy individuals leaving questions unanswered as to whether reported associations are generalizable to high-risk populations. The potential impact of PA on CVD risk may, to an extent, be mediated or modified by other CVD risk factors such as glucose intolerance, obesity, blood pressure, smoking, sex and age but this has rarely been explored. Earlier studies have not examined the potential interplay between different dimensions of PA and CVD risk factors, probably because the conventional statistical methods such as linear regression analyses are primarily focused on mean effects in the entire population and most often do not take into account interactions between multiple explanatory variables . One of the reasons is that assessment of interactions, using the traditional regression models, requires pre-specification of the interaction terms. As the number of explanatory variables in the model increases, the number of possible interactions that can be investigated is large and leads to a complicated model that can be difficult to fit and interpret. Newer tree-based regression techniques provide one approach which account for potential heterogeneity within the population being studied. One of the main advantages of tree-based methods are their data driven exploration of complex interactions among multiple risk factors without any a priori specification. Subgroups are thereby identified by the presence or absence of specific risk factors and who differ in disease risk .
Because PA patterns and CVD risk factors are represented in various combinations in individuals, looking into subgroup differences could unveil which role different dimensions of PA in daily living play in cardiovascular health. We recently identified subgroup characteristics and different activity dimensions associated with type 2 diabetes risk . Insight into which combinations of CVD risk/protective factors that are associated to higher or lower CVD incidence could likewise help defining subgroups of the population to which specific interventions, including PA, could be tailored and targeted. We hypothesized that different dimensions of PA may be associated with CVD risk in different ways. Moreover, we speculated whether the association between PA dimensions and CVD are the same in every individual or whether subgroup differences exist. Thus, we applied a tree-based approach with an extension to survival data, tree-structured survival analysis (TSSA)  with the purpose to examine the interplay between objectively measured PA dimensions and more conventional CVD risk factors in a population with low to high risk for type 2 diabetes , to assess which PA dimensions are associated with risk of first CVD event and whether subgroup differences exist.
Study design and participants
Participants for this prospective observational study were recruited between 2009 and 2011 from the Danish arm of the population-based ADDITION-Europe study (ADDITION-DK), a stepwise screening program for type 2 diabetes in primary care . ADDITION-PRO was the follow-up health examination of individuals without type 2 diabetes at screening but with low to high risk of type 2 diabetes defined by the screening procedure in ADDITION-DK. In brief, the stratification of type 2 diabetes risk in ADDITION-DK was based on a Danish diabetes risk score questionnaire . Participants were asked to indicate known risk factors for type 2 diabetes including age, sex, BMI, known hypertension, family history of type 2 diabetes, gestational diabetes and leisure time PA. Participants with a risk score ≥ 5 points (maximum 15 points) were invited to continue in the stepwise screening programme, which included random blood glucose and glycated haemoglobin A1c (HbA1c) testing, a fasting blood glucose test, and an oral glucose tolerance test. World Health Organisation criteria were used to diagnose type 2 diabetes . The sampling frame for the ADDITION-PRO study consisted of participants classified in groups of increasing type 2 diabetes risk according to the type 2 diabetes risk score and glycaemic status: low type 2 diabetes risk (< 5 points on type 2 diabetes risk score); high type 2 diabetes risk (≥ 5 points on type 2 diabetes risk score) with (i) normoglycemia; (ii) isolated impaired fasting glucose; (ii) isolated impaired glucose tolerance; and (iv) combined impaired fasting glucose and impaired glucose tolerance. In total, 2082 individuals agreed to participate in the ADDITION-PRO health examination which serves as baseline for the current study. The ADDITION-PRO study focuses on the aetiology, pathophysiology, complications and comorbidities of type 2 diabetes and is characterized by an extensive phenotyping. The study has been described previously .
Measurements and definitions
Cardiometabolic risk factors
In the ADDITION-PRO cohort, cardiometabolic risk factors were assessed at baseline in 2009–2011 and included biochemical, anthropometric, and behavioural assessments as part of the larger health examination procedure. The detailed methods of measurements or calculations of HbA1c, total cholesterol, HDL-cholesterol, LDL-cholesterol, triglycerides, plasma creatinine, urinary albumin creatinine ratio and estimated glomerular filtration rate as well as body mass index (BMI) and diastolic and systolic blood pressure have been previously published .
Information on age and sex were obtained from the unique Danish civil registration number. Information on smoking, occupation and medication use was collected from self-administered questionnaires at baseline. Smoking was categorized as current smokers and non-smokers (including prior smokers). Occupation was categorized into those who were employed and those who were unemployed (including housewives, retired individuals and individuals receiving government provisions). Medication was categorized into CVD protective medications: use of any kind of anti-diabetic drug, anti-hypertensive drug, lipid-lowering drug, and/or aspirin.
Physical activity was measured using a combined accelerometer and heart rate monitor (ActiHeart, CamNTech, Cambridge, UK) . The monitor measures uniaxial acceleration and heart rate independently. The monitor was fitted horizontally on the participant’s chest with two standard electrocardiogram electrodes (Maxensor, Alton, UK), one at the lower part of the sternum and the other one on the same horizontal level, on the left side, as laterally as possible. Each participant wore a monitor continuously for seven consecutive days. An 8-min submaximal step test was performed to individually calibrate heart rate to PA intensity . For participants who did not complete an individual calibration due to physical impairment or certain cardiovascular conditions, e.g. angina pectoris, valid calibration tests in the rest of the sample (n = 1046) were used to derive a group calibration equation adjusted for age, sex, beta-blocker use, and sleeping heart rate, for the translation of heart rate into activity intensity. Heart rate data collected during the free-living period was processed using noise classification followed by Gaussian robust regression , and average activity intensity was estimated using a branched equation framework .
The energy cost of any daily PA beyond sleep expressed as PA energy expenditure (PAEE) (kJ/kg/day), in addition to sedentary time excluding self-reported sleep, and time spent in light-, moderate- and vigorous PA were derived, whilst minimizing diurnal information bias . Time spent at various metabolic equivalent of task (MET) levels represent the PA intensity distribution across a 24 h spectrum . Sedentary time was defined as an intensity of < 1.5 METs, light PA as 1.5–3.0 METs, moderate-to-vigorous PA (MVPA) as ≥ 3.0 METs and vigorous PA as an intensity ≥ 6.0 METs, using a fixed value of 20.35 J/ml O2 × 3.5 ml O2/min/kg to define 1 MET. To indicate the degree of activity accumulation occurring in bouts, time spent in continuous MVPA lasting ≥ 10 min was derived as well as two categorical variables of meeting Danish MVPA-guidelines : 30 min of MVPA per day regardless of bout duration and 30 min of MVPA per day in bouts of at least 10 min. To ensure the best PA estimates, only measures from participants with a minimum of 48 h of ActiHeart wear time were considered valid for the present analysis.
Ascertainment of cardiovascular disease events
All Danish residents have a unique civil registration number recorded in the Danish Civil Registration System . The registers in Denmark are nationwide and cover all residents. Clinical and PA data were linked to data on diagnoses and procedures from the Danish National Patient Register  using this civil registration number. Similarly, data on death and emigration was obtained from the Central Person Register . We defined CVD as a composite outcome of events of ischemic heart disease, ischemic stroke, heart failure, atrial fibrillation/flutter and atherosclerotic disease. See Additional file 1: Table S1 for the specific International Classification of Diseases codes and Danish procedure codes used.
All the participants were followed from baseline of the study (date of the ADDITION-PRO clinical health examination) until first event of CVD, emigration, death or the end of the follow-up period (date of register extraction; May 4, 2016), whichever came first. Participants with a previous CVD diagnosis in the registers at baseline were excluded (n = 400). Also, participants with less than 48 h of monitor wear or in whom PA measurement failed were excluded (n = 233), leaving a total of 1449 individuals with 7720 person-years to be included in the present analysis.
Descriptive characteristics of the study sample at baseline by CVD status at end of follow-up were summarized as medians with inter quartile ranges (IQR) for non-normal data, means with standard deviations (SD) for normally distributed data or as numbers and percentages. To examine differences in characteristics between the groups, we conducted χ2 tests for categorical variables and t-tests for normally distributed data. In non-normal data, differences were tested using the Man-Whitney U test.
The interplay between PA and conventional CVD risk factors as well as the association between PA dimensions and CVD incidence was assessed by tree-structured survival analysis (TSSA) (“party” package in R) [12, 27]. TSSA extends tree-based methods to survival data taking into account time to event. Briefly, a survival tree is constructed using recursive partitioning . The recursive partitioning algorithm identifies the risk factor and the split in this factor which gives the maximal difference in CVD event rates between the two resulting subgroups. This procedure is applied recursively until the tree has been grown to an optimal number of “terminal leaves”. The terminal leaves comprise subgroups in the study population, characterized by a different sequence of classifications by the risk factors included in the model.
In the TSSA, the following dimensions of PA were considered as risk factors for CVD: (a) overall PAEE (b) light PA (c) MVPA (d) MVPA accumulated in bouts ≥ 10 min and (e) vigorous PA (f) meeting Danish MVPA-guidelines of 30 min of MVPA per day (bouted and non-bouted). Also, sedentary time was considered a risk factor for CVD. Additional conventional CVD risk factors included for analysis were: sex, age categorized in 5-year age bands (45–49, 50–54, …, 75–79, years) and smoking, in addition to BMI, systolic and diastolic blood pressure, HbA1c, total-, HDL- and LDL-cholesterol and triglycerides as continuous variables. Beyond their role as markers for kidney disease risk, albumin creatinine ratio and estimated glomerular filtration rate are well-established powerful independent risk factors for cardiovascular morbidity and mortality across populations with or at risk of type 2 diabetes. We included these renal biomarkers as continuous variables in the model.
Because a large proportion of the study population was using CVD protective medication at baseline, we did a sensitivity analysis, deriving a TSSA including medication use as a potential risk factor for incident CVD.
Incidence rates (IR) per 100 person-years with 95% confidence intervals (CI) for the composite CVD outcome were calculated for the entire population and for the final subgroups that emerged from the TSSA by Poisson regression analyses using log-person time as the offset variable.
Data management was performed in SAS version 9.4. Statistical analyses were performed in R version 3.1.3 (The R Foundation for Statistical Computing, http://www.R-project.org).
Characteristics of the study population
Population age ranged from 45 to 79 years with men constituting 51% of the population. Median follow-up was 5.5 years (IQR: 5.1–6.1), during which 201 (13.9%) individuals developed CVD, 57 died and 3 emigrated. The overall incidence rate for CVD was 2.6 per 100 person-years. Ischemic heart disease totalled 43% of the CVD events, ischemic stroke was 20%, atrial fibrillation/flutter 30%, heart failure was 1%, and atherosclerotic disease was 6% of the CVD events.
Participants who developed CVD during follow-up were older, had a larger proportion of men and a higher proportion of participants in anti-hypertensive and/or aspirin treatment compared to participants who did not develop CVD (Table 1). There were no baseline differences in the proportions of smokers or participants with normal glucose tolerance, pre-diabetes or type 2 diabetes between participants who did and did not develop CVD. The metabolic profile among participants who developed CVD was less favourable with higher measures of BMI, systolic blood pressure and HbA1c. Among participants who did not develop CVD, total PAEE was higher and more time was spent in both light PA, MVPA and MVPA bouts of more than 10 min compared to participants who developed CVD. None of the participants accumulated any time in vigorous PA. Participants who developed CVD spent more time sedentary (≈ 30 min) than participants who did not develop CVD.
Tree-structured survival analysis
The TSSA retained eight factors associated with incident CVD: age, HbA1c, systolic blood pressure, sex, albumin creatinine ratio, HDL-cholesterol, PAEE and estimated glomerular filtration rate (Fig. 1). Age, HbA1c and systolic blood pressure were the strongest risk factors associated with incident CVD as these were ranked highest in the tree. Interactions between the retained risk factors associated with CVD incidence and splits in these defined nine subgroups. To interpret the tree, risk factor splits and branches are followed through the tree from the top “root node” which contains all the participants until a “terminal leaf” constituting a subgroup is reached. For example, risk factors that are associated with CVD and define subgroup 1 are age (≤ 70 years), HbA1c [≤ 5.7% (39 mmol/mol)], systolic blood pressure (≤ 156 mmHg), albumin creatinine ratio (≤ 70) and PAEE (≤ 43 kJ/kg/day) while for subgroup 7 and 8 age (≤ 70 years), HbA1c [> 5.7% (39 mmol/mol)], and sex was associated with CVD, with men having the highest incidence rates for CVD.
Among the analysed PA dimensions, only PAEE was associated with incident CVD and solely among participants aged ≤ 70 years, and having HbA1c ≤ 5.7% (39 mmol/mol), systolic blood pressure ≤ 156 mmHg and albumin creatinine ratio ≤ 70. The optimal split in PAEE was at 43 kJ/kg/day resulting in two branches. One branch contained participants with PAEE levels ≤ 43 kJ/kg/day [subgroup 1; IR 3.0 (CI 1.7–4.8)] and the second included participants with PAEE > 43 kJ/kg/day. In the high PAEE subgroup, estimated glomerular filtration rate was also associated with CVD development and thereby formed two further subgroups [subgroups 2 and 3; IR 0.8 (CI 0.1–6.5) and 0.0 (CI: NA)].
After adding CVD protective medications use as risk factors for CVD in the subsidiary TSSA (sensitivity analysis), no dimension of PA was associated with CVD. Risk factors most strongly associated with incident CVD were age, use of aspirin, HbA1c, use of any anti-hypertensive drugs and HDL-cholesterol (Additional file 2: Figure S1).
We followed 1449 Danish men and women free of CVD at baseline for a median time of 5.5 years regarding first CVD events. We investigated the role of various dimensions of PA and the interplay with other CVD risk factors associated with first CVD event and identified the dimension of PA with potential benefit to cardiovascular health. The TSSA identified nine subgroups through varying combinations of risk factors and they corresponded to three levels of CVD risk. The lowest risk groups (subgroup 2, 3 and 4; IR 0.0–0.8) were characterized by younger age, higher PA level and a relatively good CVD risk factor profile (Tables 2 and 3). Five moderate CVD risk groups (subgroup 1, 6, 7, 8; IR 1.8–3.8) had varying PA characteristic and varying levels of different CVD risk factors. Finally, one high risk group (subgroup 5; IR 9.9) was characterized primarily by levels of systolic and diastolic blood pressure, LDL-cholesterol and triacylglycerol which were the highest of all of the nine groups identified by the TSSA.
Interestingly, we found PAEE to be the sole activity dimension with an association to incident CVD, independent of other CVD risk factors, and only among individuals with a relatively healthy CVD profile; this subgroup was comprised of individuals aged ≤ 70 years with HbA1c levels ≤ 5.7% (39 mmol/mol), systolic blood pressure ≤ 156 mmHg and albumin creatinine ratio ≤ 70 (Fig. 1).
Physical activity, CVD risk factors and modulation of CVD risk
Cardiovascular disease is a multifactorial disease involving many risk factors. PA positively affects the same risk factors that contribute to cardiovascular risk such as BMI, blood pressure, lipids, insulin sensitivity, blood coagulability and cardiac function (increased myocardial oxygen supply, improved myocardial contraction, and electrical stability) . Detailed molecular mechanisms of the PA-mediated prevention of CVD are not fully uncovered. However, recently it has been shown that modulations of microRNAs in an immediate response to PA are able to induce significant cardio-protection . Adherence to a physically active lifestyle over the long term may therefore modify risk factors for CVD and subsequently prime individuals for a more favourable CVD profile which in turn could lead to a lower cardiovascular risk.
Differential impact of physical activity on CVD risk in subgroups
We considered possible explanations for the lack of significant associations between any of the PA dimensions and incident CVD in the other subgroups with relatively unfavourable CVD profiles. Our diabetes risk-stratified study population makes a direct comparison with other study populations challenging. Nevertheless, our results corroborate those studies examining the impact of PA on CVD risk. Previous evidence does support more favourable effects of PA on CVD risk among low risk compared to high risk subgroups . In the Look AHEAD Study there was no significant effect of a lifestyle intervention including PA on CVD events . However, using a recursive partitioning approach as we did, Baum and colleagues recently examined whether the intervention effect differed between different subgroups of study participants. Findings showed that individuals with low HbA1c and a good general health experienced the most significant reduction in CVD events from the intervention . These findings support the notion that the impact of PA in subgroups with unfavourable CVD profiles is attenuated and thus may emphasize the importance of being physically active when CVD risk is still low. Another reason why we found no PA dimension associated with CVD in the subgroups with unfavorable CVD profiles, could be that individuals in these subgroups show homogeneous PA patterns. Lack of variation in PA patterns (and therefore lack of an independent explanatory effect) would not allow for any of the PA dimensions to further discriminate CVD event risk between mutually exclusive subgroups. Therefore, a potential CVD benefit from any of the PA dimensions cannot be ruled out in these subgroups.
There is little evidence about the potential interaction of medication use in relation to PA and CVD. While some studies find no differential effect of PA between user and non-user of CVD protective medication , others report of non-significant associations between PA and CVD after adjustment for CVD drug use . However, PA is known to be highly and inversely correlated to the use of CVD protective medication . In our study, the subgroups with the highest incidence of CVD were comprised of individuals who on average were more sedentary and with higher proportions of individuals who received various CVD protective medications compared to the other subgroups. As a result CVD protective medication use may be a proxy for physical inactivity, which could explain why inclusion of medication use in the statistical model resulted in no association of PA with CVD incidence.
Implications of increasing physical activity energy expenditure and for future physical activity guidelines
Overall PA levels in the study population as a whole were low to very low. In accordance with studies in populations similar to ours, we also found that more than 70% of the time awake was spent sedentary and only 3% of the time was spent on MVPA . Although the study population accumulated a median of 28 min/day of MVPA, the large fraction of time spent sedentary resulted in an overall low total PAEE (median 30 kJ/kg/day). As such high levels of sedentary behaviour may coexist with meeting recommended levels of MVPA. However, accumulating evidence suggests that sedentary behaviour, independent of PA levels, may be associated with an increased risk of cardiovascular disease and a variety of other health problems . Based on the evidence for the benefit of regular PA in primary and secondary prevention of CVD, the Danish and international PA guidelines for older adults (age ≥ 65) (and for adults aged 50–64 with clinically significant chronic conditions) emphasize participation in minimum 30 min/day of MVPA in bouts of more than 10 min [23, 39,40,41]. However, individuals with or at risk of type 2 diabetes generally have low levels of MVPA and spend a large portion of their day being sedentary . Due to various barriers known to limit and interfere with the ability or motivation to adhere to PA recommendations, high risk individuals do not achieve increases in MVPA of sufficient magnitude to confer health benefit [43, 44]. For example, individuals with uncomplicated type 2 diabetes exhibit an attenuated increase in stroke volume during PA attributed to impaired left ventricular filling at higher heart rates . However, the view that PA has to be of moderate-to-vigorous intensity in order to yield cardiovascular risk reduction has been debated [46,47,48]. Our findings suggest that increases in overall PAEE may have health benefits regardless of how that increase is achieved. Although higher intensity PA may offer superior cardiovascular benefits in terms of improved cardiorespiratory fitness , it may not be a prerequisite for CVD risk reduction. In fact, regular PA of even low-to-moderate intensity, typical of everyday life has shown favourable effects on cardio metabolic risk factors including HbA1C as well as CVD events [6, 49, 50]. We found that PAEE above 43 kJ/kg/day was associated with lower rates of CVD events. In a comparable study population using the same PA monitor to validate a four-category PA index, PAEE above 43 kJ/kg/day categorized individuals as being moderately active . This activity level typically represents a lifestyle for an individual with sedentary work-conditions but who are engaged in leisure time PA of > 3.5 and ≤ 7.0 h/week . As such, our finding of an association between PAEE and CVD incidence has important implications for public health. In most populations, PAEE is mainly composed of “daily activities” across the energy expenditure spectrum, rather than regular MVPA sessions. Targeting increases in PAEE can be achieved by seeking out opportunities to be more physically active in all types of activities in daily living (e.g. walking, household tasks, stair climbing). Because these activities need not necessarily be of high intensity or accumulated in a structured manner, it may be a more palatable message and feasible strategy to managing CVD risk especially for those who perhaps struggle to increase their PA level through MVPA. This argument of increasing PAEE and not just MVPA is reflected in the inclusion of sedentary behaviour recommendations and the often included “more is better” message in several recently updated PA guidelines .
Strengths and limitations
Using the ActiHeart monitor for objective PA measurement in the ADDITION-PRO study overcame the limitations, bias and poor validity associated with self-report [17, 53,54,55]. In addition, the individually calibrated and combined measures from heart rate and accelerometry has shown more accurate measures of PAEE in comparison to measures obtained by solely heart rate or solely accelerometry [18, 20, 56]. The validity and reliability of the ActiHeart monitor has been successfully tested within various populations in several studies [17, 57, 58]. The ActiHeart monitor is limited by its relatively lower feasibility and costs when used in large-scale epidemiological studies. The individual calibrations of the monitor, data cleaning and processing are time consuming and require trained staff to obtain good quality measures. Furthermore, our measure of free-living PA takes into account all types and intensities of activities performed in daily life and not only structured exercise and related structured activities. This is important in a study population that spends a considerable amount of time in low intensity, unplanned, and unstructured activities that are difficult to quantify by self-report. However, the interpretation of our results is somewhat limited by the single baseline measurement of PA. Repeated PA measures would allow for examination of whether changes in PA per se are associated with CVD. In addition, PA measures may be overestimated among participants who did not perform a step-test but for whom a group calibration value was assigned, since these participants were less healthy participants compared to the rest of the sample.
The complex interplay between PA and CVD risk factors in the development of CVD argue for studying the association of PA and several CVD risk factors with incident CVD together rather than adjusting for these factors to study the other. We applied a novel approach using TSSA to assess the possible heterogeneity of the role of PA dimensions affecting CVD development. TSSA has highly comparable metrics compared to conventional regression and conveys results in a visual and an easily interpretable manner . One constraint for the present study is collinearity between the modelled PA dimensions. However, collinearity is a bigger concern with the use of conventional regression models where correlation makes the estimated regression coefficients more difficult to interpret and increases their variability. In the TSSA, at each step, the most predictive or informative explanatory variable is selected to split a node of the tree. After splitting on the variable, its explanatory power is depleted, and other variables may be selected in subsequent steps. If two or more highly correlated explanatory variables are almost equally informative, however, the most statistically significant one is selected. More important, given the data-driven nature of the TSSA, external validation of the results is essential. Furthermore, while we have shown that the TSSA provided insight into the interplay between PA dimensions and other CVD risk factors and their association to CVD in specific subgroups, we did not perform complementary regression analyses. This is because the goal of our research was not to compare results yielded by the two statistical models.
Finally, we used register data to obtain information on CVD events. The unique Danish Civil Registration Number assigned to all Danish citizens allowed for accurate ascertainment of CVD events with no loss to follow-up. However, we also acknowledge the relative short follow-up period as one of the limitations of this study. Furthermore, during follow-up, 57 study participants died from causes other than CVD. However, because this number constitutes < 5% of the population, the effect of competing risk from non-CVD death was considered insignificant in this study.
The effect of PA in prevention of CVD events remains to be fully determined. As the number of ageing individuals and people with pre-diabetes and type 2 diabetes is predicted to increase, the prevention of CVD becomes even more urgent. Taken together, our prospective findings suggest that among various PA dimension, only PAEE is associated with the development of CVD. The fact that we found an association only in the subgroups at lowest risk of CVD may possibly indicate that in order to effectively reduce cardiovascular risk and events over time, any type of PA that increases PAEE should be implemented before a high risk state for CVD develops. Understanding the heterogeneity of the role of PA is important to guide personalized PA recommendations and translation of PA prevention studies to clinical and public health settings. Future prospective and clinical trials are needed to unravel and firmly establish the presumed CVD protective role of different PA dimensions, especially at any course and level of type 2 diabetes risk.
body mass index
glycated hemoglobin A1c
metabolic equivalent of task
moderate to vigorous physical activity
physical activity energy expenditure
tree-structured survival analysis
Ford ES, Zhao G, Li C. Pre-diabetes and the risk for cardiovascular disease: a systematic review of the evidence. J Am Coll Cardiol. 2010;55:1310–7.
Barr ELM, Zimmet PZ, Welborn TA, Jolley D, Magliano DJ, Dunstan DW, Cameron AJ, Dwyer T, Taylor HR, Tonkin AM, et al. Risk of cardiovascular and all-cause mortality in individuals with diabetes mellitus, impaired fasting glucose, and impaired glucose tolerance. Aust Diabetes Obes Lifestyle Study (AusDiab). 2007;116:151–7.
Laaksonen DE, Lindstrom J, Lakka TA, Eriksson JG, Niskanen L, Wikstrom K, Aunola S, Keinanen-Kiukaanniemi S, Laakso M, Valle TT, et al. Physical activity in the prevention of type 2 diabetes: the Finnish diabetes prevention study. Diabetes. 2005;54:158–65.
Pan XR, Li GW, Hu YH, Wang JX, Yang WY, An ZX, Hu ZX, Lin J, Xiao JZ, Cao HB, et al. Effects of diet and exercise in preventing NIDDM in people with impaired glucose tolerance. The Da Qing IGT and Diabetes Study. Diabetes Care. 1997;20:537–44.
Sofi F, Capalbo A, Cesari F, Abbate R, Gensini GF. Physical activity during leisure time and primary prevention of coronary heart disease: an updated meta-analysis of cohort studies. Eur J Cardiovasc Prev Rehabil. 2008;15:247–57.
Yates T, Haffner SM, Schulte PJ, Thomas L, Huffman KM, Bales CW, Califf RM, Holman RR, McMurray JJ, Bethel MA, et al. Association between change in daily ambulatory activity and cardiovascular events in people with impaired glucose tolerance (NAVIGATOR trial): a cohort analysis. Lancet. 2014;383:1059–66.
Koivula RW, Tornberg ÅB, Franks PW. Exercise and diabetes-related cardiovascular disease: systematic review of published evidence from observational studies and clinical trials. Curr DiabRep. 2013;13:372–80.
Li G, Zhang P, Wang J, An Y, Gong Q, Gregg EW, Yang W, Zhang B, Shuai Y, Hong J, et al. Cardiovascular mortality, all-cause mortality, and diabetes incidence after lifestyle intervention for people with impaired glucose tolerance in the Da Qing Diabetes Prevention Study: a 23-year follow-up study. Lancet Diabetes Endocrinol. 2014;2:474–80.
Moons KG, Altman DG, Reitsma JB, Ioannidis JP, Macaskill P, Steyerberg EW, Vickers AJ, Ransohoff DF, Collins GS. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): explanation and elaboration. Ann Intern Med. 2015;162:W1–73.
Lemon SC, Roy J, Clark MA, Friedmann PD, Rakowski W. Classification and regression tree analysis in public health: methodological review and comparison with logistic regression. Ann Behav Med. 2003;26:172–81.
Amadid H, Johansen NB, Bjerregaard AL, Vistisen D, Faerch K, Brage S, Lauritzen T, Witte DR, Sandbaek A, Jorgensen ME. Physical activity dimensions associated with impaired glucose metabolism. Med Sci Sports Exerc. 2017;49:2176–84.
Segal MR. Features of tree-structured survival analysis. Epidemiology. 1997;8:344–6.
Johansen NB, Hansen AL, Jensen TM, Philipsen A, Rasmussen SS, Jorgensen ME, Simmons RK, Lauritzen T, Sandbaek A, Witte DR. Protocol for ADDITION-PRO: a longitudinal cohort study of the cardiovascular experience of individuals at high risk for diabetes recruited from Danish primary care. BMC Public Health. 1078;2012:12.
Lauritzen T, Griffin S, Borch-Johnsen K, Wareham NJ, Wolffenbuttel BH, Rutten G. The ADDITION study: proposed trial of the cost-effectiveness of an intensive multifactorial intervention on morbidity and mortality among people with Type 2 diabetes detected by screening. Int J Obes Relat Metab Disord. 2000;24(Suppl 3):S6–11.
Glumer C, Carstensen B, Sandbaek A, Lauritzen T, Jorgensen T, Borch-Johnsen K. A Danish diabetes risk score for targeted screening: the Inter99 study. Diabetes Care. 2004;27:727–33.
Organization WH. Definition, diagnosis and classification of diabetes mellitus and its complications: report of a WHO consultation. Part 1, Diagnosis and classification of diabetes mellitus. (Surveillance DoND ed). DoND ed. Geneva: World Health Organization; 1999.
Brage S, Brage N, Franks PW, Ekelund U, Wareham NJ. Reliability and validity of the combined heart rate and movement sensor Actiheart. Eur J Clin Nutr. 2005;59:561–70.
Brage S, Ekelund U, Brage N, Hennings MA, Froberg K, Franks PW, Wareham NJ. Hierarchy of individual calibration levels for heart rate and accelerometry to measure physical activity. J Appl Physiol. 1985;2007(103):682–92.
Stegle O, Fallert SV, MacKay DJ, Brage S. Gaussian process robust regression for noisy heart rate data. IEEE Trans Biomed Eng. 2008;55:2143–51.
Brage S, Brage N, Franks PW, Ekelund U, Wong MY, Andersen LB, Froberg K, Wareham NJ. Branched equation modeling of simultaneous accelerometry and heart rate monitoring improves estimate of directly measured physical activity energy expenditure. J Appl Physiol. 1985;2004(96):343–51.
Brage S, Westgate K, Franks PW, Stegle O, Wright A, Ekelund U, Wareham NJ. Estimation of free-living energy expenditure by heart rate and movement sensing: a doubly-labelled water study. PLoS ONE. 2015;10:e0137206.
McArdle WD. Exercise physiology: energy, nutrition, and human performance. 4th ed. Philadelphia: Lea & Febiger; 1996.
Physical activity recommendations for adults. https://www.sst.dk/en/health-and-lifestyle/physical-activity/recommendations/recommendations-for-adults. Accessed Jan 2018.
Pedersen CB, Gotzsche H, Moller JO, Mortensen PB. The Danish civil registration system. A cohort of eight million persons. Dan Med Bull. 2006;53:441–9.
Lynge E, Sandegaard JL, Rebolj M. The Danish national patient register. Scand J Public Health. 2011;39:30–3.
Ebbesen AV. The creation of the central person registry in Denmark. In: Gram C, Rasmussen P, Østergaard SD, editors. History of Nordic computing 4: 4th IFIP WG 97 Conference, HiNC 4, Copenhagen, Denmark, August 13–15; 2014, Revised Selected Papers. Cham: Springer International Publishing; 2015. p. 49–57.
Package ‘party’ A laboratory for recursive partytioning. http://party.r-forge.r-project.org/. Accessed Jan 2018.
Zhang H, Burton S. Recursive partitioning and applications. 2nd ed. New York: Springer; 2010.
Paffenbarger RSLI. Exercise and fitness. In: Manson JERP, Gaziano JM, Hennekens CH, editors. Prevention of myocardial infarction. New York: Oxford University Press; 1996. p. 172–202.
Lew JK, Pearson JT, Schwenke DO, Katare R. Exercise mediated protection of diabetic heart through modulation of microRNA mediated molecular pathways. Cardiovasc Diabetol. 2017;16:10.
Fritz T, Caidahl K, Krook A, Lundstrom P, Mashili F, Osler M, Szekeres FL, Ostenson CG, Wandell P, Zierath JR. Effects of Nordic walking on cardiovascular risk factors in overweight individuals with type 2 diabetes, impaired or normal glucose tolerance. Diabetes Metab Res Rev. 2013;29:25–32.
Wing RR, Bolin P, Brancati FL, Bray GA, Clark JM, Coday M, Crow RS, Curtis JM, Egan CM, Espeland MA, et al. Cardiovascular effects of intensive lifestyle intervention in type 2 diabetes. N Engl J Med. 2013;369:145–54.
Baum A, Scarpa J, Bruzelius E, Tamler R, Basu S, Faghmous J. Targeting weight loss interventions to reduce cardiovascular complications of type 2 diabetes: a machine learning-based post hoc analysis of heterogeneous treatment effects in the Look AHEAD trial. Lancet Diabetes Endocrinol. 2017;5:808–15.
Stamatakis E, Hamer M, Primatesta P. Cardiovascular medication, physical activity and mortality: cross-sectional population study with ongoing mortality follow-up. Heart. 2009;95:448–53.
Bertoni AG, Whitt-Glover MC, Chung H, Le KY, Barr RG, Mahesh M, Jenny NS, Burke GL, Jacobs DR. The association between physical activity and subclinical atherosclerosis: the multi-ethnic study of atherosclerosis. Am J Epidemiol. 2009;169:444–54.
Williams PT, Franklin B. Vigorous exercise and diabetic, hypertensive, and hypercholesterolemia medication use. Med Sci Sports Exerc. 1933;2007:39.
Henson J, Yates T, Biddle SJH, Edwardson CL, Khunti K, Wilmot EG, Gray LJ, Gorely T, Nimmo MA, Davies MJ. Associations of objectively measured sedentary behaviour and physical activity with markers of cardiometabolic health. Diabetologia. 2013;56:1012–20.
Healy GN, Dunstan DW, Salmon J, Shaw JE, Zimmet PZ, Owen N. Television time and continuous metabolic risk in physically active adults. Med Sci Sports Exerc. 2008;40:639–45.
Nelson ME, Rejeski WJ, Blair SN, Duncan PW, Judge JO, King AC, Macera CA, Castaneda-Sceppa C. Physical activity and public health in older adults: recommendation from the American College of Sports Medicine and the American Heart Association. Med Sci Sports Exerc. 2007;39:1435–45.
Winzer EB, Woitek F, Linke A. Physical activity in the prevention and treatment of coronary artery disease. J Am Heart Assoc Cardiovasc Cerebrovasc Dis. 2018;7:e007725.
Zachariah G, Alex A. Exercise for prevention of cardiovascular disease: evidence-based recommendations. J Clin Prev Cardiol. 2017;6:109–14.
Morrato EH, Hill JO, Wyatt HR, Ghushchyan V, Sullivan PW. Physical activity in U.S. adults with diabetes and at risk for developing diabetes, 2003. Diabetes Care. 2007;30:203–9.
Avery L, Flynn D, van Wersch A, Sniehotta FF, Trenell MI. Changing physical activity behavior in type 2 diabetes: a systematic review and meta-analysis of behavioral interventions. Diabetes Care. 2012;35:2681–9.
Lidegaard LP, Schwennesen N, Willaing I, Faerch K. Barriers to and motivators for physical activity among people with Type 2 diabetes: patients’ perspectives. Diabet Med. 2016;33:1677–85.
Wilson GA, Wilkins GT, Cotter JD, Lamberts RR, Lal S, Baldi JC. Impaired ventricular filling limits cardiac reserve during submaximal exercise in people with type 2 diabetes. Cardiovasc Diabetol. 2017;16:160.
Sparling PB, Howard BJ, Dunstan DW, Owen N. Recommendations for physical activity in older adults. BMJ. 2015;350:h100.
de Souto Barreto P. Global health agenda on non-communicable diseases: has WHO set a smart goal for physical activity? BMJ. 2015;350:h23.
Smith L, Ekelund U, Hamer M. The potential yield of non-exercise physical activity energy expenditure in public health. Sports Med. 2015;45:449–52.
Grace A, Chan E, Giallauria F, Graham PL, Smart NA. Clinical outcomes and glycaemic responses to different aerobic exercise training intensities in type II diabetes: a systematic review and meta-analysis. Cardiovasc Diabetol. 2017;16:37.
Wen CP, Wai JP, Tsai MK, Yang YC, Cheng TY, Lee MC, Chan HT, Tsao CK, Tsai SP, Wu X. Minimum amount of physical activity for reduced mortality and extended life expectancy: a prospective cohort study. Lancet. 2011;378:1244–53.
Ekelund U, Palla L, Brage S, Franks PW, Peters T, Balkau B, Diaz MJ, Huerta JM, Agnoli C, Arriola L, et al. Physical activity reduces the risk of incident type 2 diabetes in general and in abdominally lean and obese men and women: the EPIC-InterAct Study. Diabetologia. 2012;55:1944–52.
Health Council of the Netherlands, Physical activity guidelines 2017. Hague; 2017. https://www.gezondheidsraad.nl/sites/default/files/grpublication/physical_activity_guidelines_2017_0.pdf. Accessed Aug 2018.
Wahid A, Manek N, Nichols M, Kelly P, Foster C, Webster P, Kaur A, Friedemann Smith C, Wilkins E, Rayner M, et al. Quantifying the association between physical activity and cardiovascular disease and diabetes: a systematic review and meta-analysis. J Am Heart Assoc Cardiovasc Cerebrovasc Dis. 2016;5:e002495.
Besson H, Brage S, Jakes RW, Ekelund U, Wareham NJ. Estimating physical activity energy expenditure, sedentary time, and physical activity intensity by self-report in adults. Am J Clin Nutr. 2010;91:106–14.
Shephard RJ. Limits to the measurement of habitual physical activity by questionnaires. Br J Sports Med. 2003;37:197–206 (discussion 206).
Crouter SE, Churilla JR, Bassett DR Jr. Accuracy of the Actiheart for the assessment of energy expenditure in adults. Eur J Clin Nutr. 2008;62:704–11.
Assah FK, Ekelund U, Brage S, Wright A, Mbanya JC, Wareham NJ. Accuracy and validity of a combined heart rate and motion sensor for the measurement of free-living physical activity energy expenditure in adults in Cameroon. Int J Epidemiol. 2011;40:112–20.
Villars C, Bergouignan A, Dugas J, Antoun E, Schoeller DA, Roth H, Maingon AC, Lefai E, Blanc S, Simon C. Validity of combining heart rate and uniaxial acceleration to measure free-living physical activity energy expenditure in young men. J Appl Physiol. 1985;2012(113):1763–71.
Colombet I, Ruelland A, Chatellier G, Gueyffier F, Degoulet P, Jaulent MC. Models to predict cardiovascular risk: comparison of CART, multilayer perceptron and logistic regression. In: Proceedings of the AMIA Symposium. 2000:156–160.
HA and DV contributed to study design. HA performed the analysis and wrote the manuscript. DV helped with analyzing data and contributed to revising the article. MJ, AS, NJ, AB, SB, KF and DW contributed with interpretation of data and revising the content of the manuscript. HA is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. The authors declare that the results of the study are presented clearly, honestly, and without fabrication, falsification, or inappropriate data manipulation. All authors read and approved the final manuscript.
The ADDITION-Denmark study was supported by the National Health Services in the counties of Copenhagen, Aarhus, Ringkøbing, Ribe, and Southern Jutland in Denmark; the Danish Council for Strategic Research; the Danish Research Foundation for General Practice; Novo Nordisk Foundation; the Danish Centre for Evaluation and Health Technology Assessment; the Diabetes Fund of the National Board of Health; the Danish Medical Research Council; and the Aarhus University Research Foundation. Additionally, the ADDITION-PRO study was funded by an unrestricted grant from the European Foundation for the Study of Diabetes/Pfizer for Research into Cardiovascular Disease Risk Reduction in Patients with Diabetes (74550801), the Danish Council for Strategic Research, and internal research and equipment funds from Steno Diabetes Center.
We thank Research Centre for Prevention and Health, Capital Region, Denmark for the opportunity given to Hanan Amadid to carry out a research stay at the research center. The supervision and guidance in register data management received during her stay was essential for the present study.
A.B., K.F, M.J., T.L., D.V. and D.W. own shares in Novo Nordisk A/S. No other potential competing interest relevant to this article were reported for the remaining authors.
Availability of data and materials
The data that support the findings of this study are available from the Danish Data Protection Agency and Danish Health Data Agency (j-No: 2012-58-0004 and FSEID-00002001) but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Danish Data Protection Agency and Danish Health Data Agency (j-No: 2012-58-0004 and FSEID-00002001).
Consent for publication
Ethics approval and consent to participate
Ethical approval for the ADDITION-PRO study was obtained from the ethics committee of the Central Denmark Region (Reference No. 20080229) and the study was conducted in accordance with the Helsinki Declaration. All participants provided oral and written informed consent before participating in the ADDITION-PRO study. According to Danish law, ethical approval and patient consent is not required for the register-part of this study. Access and use of the described register data are approved by the Danish Data Protection Agency and Danish Health Data Agency (j-No: 2012-58-0004 and FSEID-00002001).
H.A. has received scholarship funding from the Danish Ministry of Higher Education and Science, the Danish Heart Foundation, the Novo Nordisk Foundation and Aarhus University. D.W. is supported by the Danish Diabetes Academy, which is funded by the Novo Nordisk Foundation. K.F. is supported by the Novo Nordisk Foundation. The work of S.B is supported by UK Medical Research Council [MC_UU_12015/3].
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Amadid, H., Johansen, N.B., Bjerregaard, AL. et al. The role of physical activity in the development of first cardiovascular disease event: a tree-structured survival analysis of the Danish ADDITION-PRO cohort. Cardiovasc Diabetol 17, 126 (2018). https://doi.org/10.1186/s12933-018-0769-x