Article Text
Abstract
Objective Although 30 min/day of moderate-intensity physical activity is suggested for preventing type 2 diabetes (T2D), the current recommendations exclusively rely on self-reports and rarely consider the genetic risk. We examined the prospective dose-response relationships between total/intensity-specific physical activity and incident T2D accounting for and stratified by different levels of genetic risk.
Methods This prospective cohort study was based on 59 325 participants in the UK Biobank (mean age=61.1 years in 2013–2015). Total/intensity-specific physical activity was collected using accelerometers and linked to national registries until 30 September 2021. We examined the shape of the dose-response association between physical activity and T2D incidence using restricted cubic splines adjusted for and stratified by a polygenic risk score (based on 424 selected single nucleotide polymorphisms) using Cox proportional hazards models.
Results During a median follow-up of 6.8 years, there was a strong linear dose-response association between moderate-to-vigorous-intensity physical activity (MVPA) and incident T2D, even after adjusting for genetic risk. Compared with the least active participants, the HRs (95% CI) for higher levels of MVPA were: 0.63 (0.53 to 0.75) for 5.3–25.9 min/day, 0.41 (0.34 to 0.51) for 26.0–68.4 min/day and 0.26 (0.18 to 0.38) for >68.4 min/day. While no significant multiplicative interaction between physical activity measures and genetic risk was found, we found a significant additive interaction between MVPA and genetic risk score, suggesting larger absolute risk differences by MVPA levels among those with higher genetic risk.
Conclusion Participation in physical activity, particularly MVPA, should be promoted especially in those with high genetic risk of T2D. There may be no minimal or maximal threshold for the benefits. This finding can inform future guidelines development and interventions to prevent T2D.
- physical activity
- endocrine system
- genetics
- epidemiology
- preventive medicine
Data availability statement
The data analyzed in this study are available from the UK Biobank website with approved access. Data could be obtained upon direct application to the UK Biobank Study.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
WHAT IS ALREADY KNOWN ON THIS TOPIC
The WHO recommends at least 30 min of moderate-intensity physical activity per day for the prevention of type 2 diabetes (T2D). However, such recommendations almost exclusively rely on studies using self-reported measures of physical activity, which is subject to bias. The dose-response association between device-measured physical activity, particularly light-intensity physical activity, and T2D is rarely explored. Furthermore, while genetics play a major role in the onset of T2D, the interaction between genetic risk and physical activity on T2D is less understood.
WHAT THIS STUDY ADDS
First, there was a strong linear dose-response inverse association between accelerometer-measured physical activity and T2D even when adjusting for the genetic risk. The magnitude of the association appeared larger than in studies using self-reported physical activity. Second, this association was stronger for moderate-to vigorous-intensity physical activity (MVPA), while the association with light-intensity physical activity was weak and inconsistent. Third, we found no multiplicative interaction between physical activity and genetic risk but a significant additive interaction with MVPA. The association between total physical activity/MVPA and T2D was similar across genetic risk strata but the absolute risk reduction from MVPA was the largest for those with high genetic risk.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
In contrast to existing evidence primarily based on self-reported physical activity, our findings imply that there may be no threshold for the ‘optimal amount’ of physical activity for T2D prevention and the association may be stronger than previously estimated. These findings can inform future guidelines development and interventions to prevent T2D. Physical activity, specifically MVPA, is beneficial, especially in those with high genetic risk, and should be promoted as a priority strategy for T2D prevention.
Introduction
Diabetes remains, to date, a global public health concern.1 In 2021, there were 537 million adults living with diabetes worldwide, with a prevalence of 10.5%. The direct health expenditures associated with diabetes during the same year ascended to US$966 billion.2 Type 2 diabetes (T2D) accounts for approximately 90% of all types of diabetes mellitus.
Physical activity is a first-line strategy for the prevention and management of T2D.3 4 Health authorities recommended at least 30 min of moderate-intensity physical activity (MPA) on most days of the week or 150 min accumulated per week to prevent the onset of T2D.1 4 However, such recommendations almost exclusively rely on studies using self-reported measures of physical activity, which is subject to bias. In addition, the role of low-intensity physical activity (LPA) on T2D remains unknown. As a result, the 2018 Physical Activity Guidelines Advisory Committee Scientific Report called for more research on device-based physical activity and the potential health benefits of LPA.5
Although the important role genetics plays in the aetiology of T2D is well established, evidence on whether the genetic risk of T2D could be attenuated by physical activity, and whether the protective effects of physical activity on T2D incidence differ by genetic risk remains scarce. Recent advancements in genome-wide association analyses have led to the identification of more than 400 single nucleotide polymorphisms (SNPs) for T2D,6 7 allowing for in-depth investigation of the interactions between physical activity and genetic risk on incident T2D.
Capitalising on a large sample of adults with available accelerometry, genetic and prospective health data, this study aimed to examine (1) the prospective dose-response relationships between total and intensity-specific physical activity and incident T2D while accounting for genetic risk; and (2) the association between total and intensity-specific physical activity and incident T2D across different levels of genetic risk, and the associations between combinations of genetic risk and physical activity levels and incident T2D.
Methods
Study design and participants
We used data from the UK Biobank (reference number: 63454), a prospective cohort study of over 500 000 participants aged 40–69 years from the UK, with detailed information reported elsewhere.8 Briefly, between the years 2006 and 2010, baseline assessments were carried out in 22 assessment centres comprising a self-completed touch-screen questionnaire in addition to physical and functional measures and collection of a fresh blood sample. Additionally, between February 2013 and December 2015, a subsample of 103 712 participants wore an accelerometer on their dominant wrist for seven consecutive 24-hour days.9
For the current study, we excluded participants with a diagnosis of prevalent diabetes (n=4634; online supplemental table 1), cardiovascular diseases (CVDs) or cancer (n=18 545) before accelerometer measurement. We also excluded those with no genetic data (n=1537) and those with missing values on covariates (n=8119) (online supplemental figure 1).
Supplemental material
Measures
Accelerometer-measured physical activity
Raw device-based physical activity data were collected using an accelerometer (Axivity AX3; designed by Open Lab, Newcastle University) and processed following the protocol described elsewhere.10 11 The Euclidean Norm Minus One (ENMO) in milligravity (mg) was used to quantify the acceleration related to the movement registered.12 A validated tailored machine learning model using balanced random forests with Hidden Markov models was then used to identify different movement behaviours in 30 s time windows, including moderate-to-vigorous-intensity physical activity (MVPA), LPA, sedentary behaviour and sleep.9 The total volume of physical activity was summarised as the average ENMO while awake (defined as 06:00 to 22:00), while daily MVPA and LPA were defined as the sum of all minutes spent on waking behaviours at ≥3 and 1.5–2.9 metabolic equivalent of tasks (METs) in a typical day.
Genetic risk
Genetic risk was estimated by applying a polygenic risk score (PRS) including 424 selected SNPs identified from genome-wide association analyses.6 7 The detailed genotyping process, imputation and quality control in the UK Biobank have been described elsewhere.13 A PRS was calculated using the following formula with each SNP weighted: PRS=β1×SNP1+β2×SNP2+ …+βn×SNPn, where βn was the relative effect size and SNPn is the risk allele number of each SNP (online supplemental table 2). A higher PRS indicates a higher genetic predisposition to developing T2D. PRSs were categorised by tertile into low, intermediate or high genetic risk.
Incident T2D
Prevalent diabetes was identified based on the modified Biobank algorithms by Eastwood et al,14 via hospital inpatient and outpatient records, primary care data, death registration, self-reported medical history and medication, as well as biochemical examination for blood glucose (if random GLU ≥11.1 mmol/L) and glycated haemoglobin (if HbA1c ≥6.5%). Incident T2D was defined as developing T2D (both fatal and nonfatal) after completing an accelerometer assessment up until the censoring date (30 September 2021 for England, 31 July 2021 for Scotland and 28 February 2018 for Wales), ascertained by using the code E11 from the International Classification of Diseases 10th Revision (ICD-10) (online supplemental table 1).7 15 We calculated the follow-up time from the accelerometer assessment to the time of T2D diagnosis, death, loss to follow-up or censorship, whichever occurred first.
Covariates
Potential confounders were selected based on an a priori developed acyclic graph (online supplemental figure 2). Sociodemographic characteristics included age (continuous in years), self-reported gender (men; women), ethnicity (white European; non-white European), educational attainment (no qualifications; other qualifications than college/university degree; college or university degree), household income (<18 000; 18 000–30 999; 31 000–51 999; >52 000 £/year), Townsend deprivation index (quartiles), employment status (unemployed; employed; retired) and assessment centres (22 categories). Additional covariates included non-physical activity lifestyle factors: smoking status (never; former; current), alcohol consumption (never or special occasions only; one to three times a month; once or twice a week; three or four times a week; daily or on most days), healthy diet score (quartiles) based on a previously validated dietary index (online supplemental table 3) and body mass index (BMI; kg/m2).16 Pre-existing chronic conditions were defined as either a self-reported history of hypertension, dyslipidaemia and depression at the baseline or diagnosed diseases through health records between the time of baseline assessment (the year 2006–2010) and accelerometer measurement (the year 2013–2015). Finally, we adjusted for accelerometer-related variables, including total wear days and seasonality.
Statistical analysis
Complete case analysis was used in this study. We presented descriptive statistics in counts and percentages for categorical variables, and means and SD for continuous variables. The continuous dose-response analyses assessed the shape of associations between total/intensity-specific physical activity and incident T2D, with data trimmed at the 5th and 95th percentile of the exposure distribution.17 Restricted cubic splines were used to allow for potential nonlinearity with three knots placed at the 10th, 50th (reference) and 90th percentiles. We assumed linearity for values of the bottom and top 10%. Departure from linearity was examined by a Wald test.
Total physical activity, MVPA and LPA were then categorised into 4-level variables based on the 10th, 50th and 90th percentiles of data distribution, separately. We conducted multivariable Cox proportional hazards models to estimate the HRs and 95% CIs for incident T2D according to the categories of physical activity exposures, using age as the underlying timescale. The proportional hazards assumption was checked using the Schoenfeld residuals, and no violation was found.
Five sequential models were built: Model 1 accounted for age as the timescale and adjusted for gender and accelerometer-related variables as covariates. Model 2 additionally adjusted for other sociodemographic characteristics and Model 3 further adjusted for non-physical activity lifestyle factors and pre-existing chronic conditions. We further adjusted for PRS, genotyping array and the first 10 principal components of ancestry in Model 4. For the intensity-specific physical activity model, we mutually adjusted for MVPA and LPA in Model 5. We used Model 4 as the main model for total physical activity and Model 5 as the main model for intensity-specific physical activity.
We tested PRS as a potential effect modifier by adding a multiplicative interaction term (total physical activity×PRS in Model 4; MVPA×PRS and LPA×PRS in Model 5 one at a time), followed by stratified analysis by PRS tertiles (low, intermediate and high risk). Considering the different baseline hazards in models for each stratum, we also examined the joint association of physical activity and genetic risk by creating (4×3) mutually exclusive categories combining physical activity levels and PRS with the highest risk combination (ie, highest PRS and lowest physical activity level) as the reference.
We conducted several sensitivity analyses. First, to reduce the risk of reverse causation, we conducted a sensitivity analysis by left truncating the first 2 years of the follow-up period. Second, considering that the current accelerometer processing algorism (eg, setting 06:00 to 22:00 as waking time) may not accurately classify physical activity for shift workers, we conducted a sensitivity analysis excluding shift workers (n=5111). Third, we reran the models defining awake time as 07:00 to 21:00 and 08:00 to 20:00. Fourth, considering that BMI could be a potential confounder or a mediator, and the current main models treated BMI as a mediator and therefore did not adjust for it, we conducted an additional sensitivity analysis adjusting for BMI as a confounder.18 19 Lastly, as a post-hoc sensitivity analysis, we performed Aalen’s additive hazards models to test the additive interaction effects.20 The statistical analysis plan was pre-registered on Open Science Framework (Identifier: DOI 10.17605/OSF.IO/9ST7J). We conducted all analyses in SAS (V.9.4) and R (V.3.6.0). Results were reported according to STrengthening the Reporting of OBservational studies in Epidemiology Checklists (online supplemental table 11).
Equity, diversity and inclusion statement
Our research team included two women and three men, four authors are people of colour. All authors are at the early-stage or mid-stage of their careers. The UK Biobank sample was population-based (including a broad range of sociodemographic characteristics), but not population-representative, which we acknowledged as a limitation in the Discussion section.
Results
The final analytical sample included 59 325 participants without diabetes, CVDs or cancer at the time of accelerometer assessment, with valid physical activity, genetic data and all covariates (online supplemental figure 1).
Table 1 presents the baseline characteristics of the analytical sample by levels of total physical activity. Overall, 44.0% of the participants were men; 97.1% were of white European descent; 46.6% had a college or university degree, over 60% were employed and a similar proportion had a household income of over 31 000 £/year. The prevalence of current smoking, poor diet, overweight and obesity, hypertension, dyslipidaemia and depression was the highest among the group with the lowest level of total physical activity. Additional information on the baseline characteristics of the sample by different levels of MVPA and LPA are displayed in online supplemental tables 4 and 5, respectively.
Physical activity and incident T2D
We observed 884 incident T2D cases during a median follow-up period of 6.8 years (IQR=6.3–7.3). Table 2 shows the HRs and 95% CIs for the association between total physical activity/MVPA/LPA and incident T2D.
For overall physical activity, compared with the least active participants, those in the most active group had an 80% lower risk of incident T2D in Model 1. The association was attenuated after the covariates were introduced. In the final model with PRS adjusted (Model 4), the inverse relationship between total physical activity and incident T2D remained significant with the most active group having a 68% lower risk of incident T2D. When examining intensity-specific physical activity, MVPA had a strong association with T2D (HR=0.26; 95% CI=0.18 to 0.38 in the most active group), even after adjusting for LPA and all other covariates. However, for LPA, we only found a statistically significant association with incident T2D at the highest level (top 10%). The dose-response analysis revealed linear relationships between all physical activity exposures and incident T2D (figure 1).
Stratified and joint analysis
There is a strong positive association between genetic risk score and incident T2D (HR=1.41; 95% CI=1.16 to 1.71 for the intermediate risk group and HR=2.43; 95% CI=2.04 to 2.90 for the high genetic risk group, as compared with low genetic risk group) (online supplemental table 6). Stratified by genetic risk, within each category, total physical activity and MVPA were consistently associated with incident T2D, but LPA was not (table 3). For example, compared with those with the lowest level of total physical activity, participants with the highest level had a 66%, 76% and 61% lower risk of T2D in the high, intermediate and low genetic risk groups. There was no significant multiplicative interaction based on p value and stratified analyses. However, we found a significant additive interaction with MVPA (p=0.011), but not with total physical activity (p=0.113) or LPA (p=0.524). Specifically, we observed absolute risk differences of 3.8, 2.1 and 1.0 cases per 100 000 person-years between the most and least active groups based on MVPA in those with high, intermediate and low genetic risk, respectively (online supplemental figure 3). Within each genetic risk subgroup, we observed stronger associations between MVPA and T2D compared with those between LPA and T2D.
Figure 2 shows the association between joint physical activity/PRS categories and T2D. Overall, results based on the total physical activity and MVPA showed similar patterns. In both cases, high genetic risk and the highest total physical activity/MVPA combination were associated with a lower risk of incident T2D than low genetic risk and the lowest total physical activity/MVPA combination. In contrast, associations between LPA/PRS combinations and T2D were less clear—those with high genetic risk and the lowest three LPA categories had similar HRs and all the other combinations had similarly smaller HRs.
Sensitivity analysis
Left truncating the first 2 years of follow-up did not appreciably change the associations between total physical activity/MVPA and incident T2D (online supplemental table 6), but the highest level of LPA was no longer associated with incident T2D. Excluding shift workers (n=5111) and redefining awake time as 07:00 to 21:00 or 08:00 to 20:00 did not substantially alter any estimates (online supplemental tables 7–9). Furthermore, the relationships between total/intensity-specific physical activity and incident T2D were attenuated after adjusting for BMI, and no inverse relationship was found between LPA and T2D incidence after accounting for BMI (online supplemental table 10).
Discussion
Based on the largest accelerometry dataset paired with the most comprehensive genetic risk indicators to our knowledge, our study showed that higher levels of total physical activity and particularly MVPA were strongly associated with a lower risk of developing T2D, independent of genetic risk. Considering that physical inactivity is a modifiable risk factor in most cases, our findings provide insights for informing public health recommendations to prevent T2D. Particularly, we found a linear dose-response association with a considerably larger magnitude than studies using self-reported physical activity. The association was the strongest when physical activity was operationalised as MVPA, followed by total physical activity, and the association between LPA and T2D was the least consistent and weakest. Finally, although we did not find a significant effect modification by genetic risk on a multiplicative scale, we found a significant additive interaction indicating a larger absolute risk reduction among those with the highest genetic risk.
Our results confirm previously observed protective effects of physical activity on the onset of T2D, regardless of genetic risk. In our study, 68.4 min/day of MVPA (ie, 90th percentile) was associated with a 74% lower risk of developing T2D when compared with participants who spent less than 5.2 min in MVPA (ie, 10th percentile). This association is of a larger magnitude than previous studies based on self-reported physical activity,3 which was also observed for mortality outcomes.21 This is likely explained by the measurement errors in self-reported physical activity that tend to bias the association towards the null.22 23
Remarkably, a linear dose-response pattern with no minimal or maximal threshold was identified. This finding contrasts previous evidence supporting a curvilinear dose-response association whereby the association became weaker at higher levels of physical activity.5 For example, Kyu and colleagues’ dose-response meta-analysis including 55 prospective cohort studies with self-reported physical activity found a clear threshold effect,24 where the largest risk reduction was observed when physical activity was increased from 600 to 3600 MET-min/week (equivalent to approximately from 150 to 900 min/week of MPA). Beyond this threshold, increasing levels of physical activity were associated with little improvements in diabetes outcomes. However, reporting biases, particularly at the higher levels of physical activity, may have potentially distorted the observed associations. Our results indicate that individuals should be encouraged to be as physically active as possible to maximise the benefits.
When examining intensity-specific physical activity, the association between MVPA and incident T2D was strong and significant at every level when compared with the bottom 10% MVPA level (≤5.2 min/week). In contrast, only participants in the top 10% LPA level (over 7.4 hours/day) had a statistically significantly lower risk of developing T2D and the association became non-significant in several sensitivity analyses. Due to the difficulty of measuring LPA using self-reported instruments, the evidence base for LPA and health outcomes is less developed than that for MVPA. The large amount of LPA required to reduce the risk of T2D may, at least in part, be explained by the interaction of the insulin system with different intensities of physical activity. A position statement by the American Diabetes Association pointed out that 20 min of MVPA can help elevate glucose uptake and enhance insulin action, while far more LPA is needed to achieve a similar effect.4 Evidence from our study implies that, when possible, it may be more effective for T2D prevention programmes to focus on promoting MVPA.
We did not observe a multiplicative interaction between physical activity and genetic risk using 424 SNPs among the UK Biobank participants. However, we found a significant additive interaction with MVPA, suggesting the largest absolute risk reduction in those with a high genetic predisposition to T2D. Interestingly, our results showed that participants with a high genetic risk and the highest total physical activity and MVPA (ie, upper 10th percentile, at least 68.4 min/day of MVPA) had a lower risk of incident T2D than those with a low genetic risk and the lowest total physical activity/MVPA (ie, bottom 10th percentile, less than 5.2 min/day of MVPA). These results underscore the importance of being physically active for T2D prevention, particularly in those with high genetic risk. Our finding is in line with two studies using 19 SNPs in Poland and 49 SNPs in nine European countries.25 26 In contrast, Klimentidis and his colleagues found that the protective effect of physical activity in those with higher genetic risk was weaker compared with those with lower genetic risk using a PRS comprised of 65 SNPs in the USA.27 The difference between studies is likely due to different measurements of physical activity (ie, our study used accelerometers while previous studies used self-reported physical activity) and genetic risk (our study captured more T2D-related SNPs than previous studies).28
To our knowledge, this investigation is the first population-based prospective cohort study to explore the relationship between accelerometer-measured physical activity and incident T2D while accounting for genetic risk. Our study has several strengths. First, while most previous studies used self-reported physical activity measures, which are prone to recall and social desirability biases, we used accelerometers that monitor physical activity of different intensities continuously through 24-hour days. Second, we used machine learning approache to quantify the volume of intensity-specific physical activity instead of using pre-defined cut points, which are subject to misclassifications, and have shown substantial agreement with wearable camera ground-truth data.29 30 Furthermore, the machine learning model used in the current study was trained and validated in a realistic free-living environment, which covered a wider variety of activities than in previous lab-based studies.31 Third, we included 424 selected SNPs for the genetic risk of T2D, which is the most comprehensive genetic study on physical activity and T2D available to date.
This study is subject to limitations. Physical activity was only measured at a single time point which could not capture long-term patterns of physical activity. Additionally, participants of the UK Biobank study are not nationally representative of the general UK population; however, a previous study found that the association between risk factors and health outcomes derived from the study may be generalisable.32 Lastly, since the accelerometry data were collected between 2013 and 2015, the follow-up time is relatively short.
In conclusion, our results reinforced the importance of physical activity for T2D prevention, particularly for those with a high genetic risk. Future guideline development should consider physical activity intensity and the dose-response nature of the association. The promotion of physical activity, particularly MVPA, should be a priority strategy for T2D prevention.
Data availability statement
The data analyzed in this study are available from the UK Biobank website with approved access. Data could be obtained upon direct application to the UK Biobank Study.
Ethics statements
Patient consent for publication
Ethics approval
This study involves human participants and was approved by North West Multi-centre Research Ethics Committee. Participants gave informed consent to participate in the study before taking part.
Acknowledgments
We gratefully acknowledge the commitment and dedication of the participants of the UK Biobank.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Twitter @SusanLuo1, @DrMelodyDing
LC and DD contributed equally.
Contributors ML, LC and DD conceptualised and designed the study. ML analysed the data with the help of CY, BDPC, LC and DD, and drafted the study with DD. All the authors reviewed the results, provided input during critical revisions and approved the submitted manuscript. All authors have final responsibility for the decision to submit for publication.
Funding DD is supported by a National Health and Medical Research Council Investigator Grant (# 2009254). LC is supported by the Young Elite Scientists Sponsorship Program by China Association for Science and Technology (grant number: YESS20210143). BDPC is supported by the Government of Andalusia, Research Talent Recruitment Program (EMERGIA #2020/00158).
Competing interests None declared.
Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.