Skip to main content

Association of physical activity and air pollution exposure with the risk of type 2 diabetes: a large population-based prospective cohort study

Abstract

Background

The interplay between physical activity (PA) and air pollution in relation to type 2 diabetes (T2D) remains largely unknown. Based on a large population-based cohort study, this study aimed to examine whether the benefits of PA with respect to the risk of T2D are moderated by exposure to air pollution.

Methods

UK Biobank participants (n = 359,153) without diabetes at baseline were included. Information on PA was obtained using the International Physical Activity Questionnaire short form. Exposure to air pollution, including PM2.5, PMcoarse (PM2.5−10), PM10, and NO2, was estimated from land use regression models. Cox regression models were used to estimate the hazard ratios (HRs) and 95% confidence intervals (95% CIs).

Results

During a median of 8.9 years of follow-up, 13,706 T2D events were recorded. Compared with a low PA level, the HRs for the risk of T2D among individuals with moderate and high PA were 0.82 (95% CI, 0.79–0.86) and 0.73 (95% CI, 0.70–0.77), respectively. Compared with low levels of air pollution, the HRs for risk of T2D for high levels of air pollution (PM2.5, PMcoarse, PM10, and NO2) were 1.19 (1.14–1.24), 1.06 (1.02–1.11), 1.13 (1.08–1.18), and 1.19 (1.14–1.24), respectively. There was no effect modification of the associations between PA and T2D by air pollution (all P-interactions > 0.05). The inverse associations between PA and T2D in each air pollution stratum were generally consistent (all P for trend < 0.05).

Conclusion

A higher PA and lower air pollution level were independently associated with a lower risk of T2D. The beneficial effects of PA on T2D generally remained stable among participants exposed to different levels of air pollution. Further studies are needed to replicate our findings in moderately and severely polluted areas.

Peer Review reports

Background

The global burden of type 2 diabetes (T2D) has increased over the past few decades, and its prevention is a public health priority [1]. Physical activity (PA) has been suggested to play an important role in the prevention of T2D and related morbidities [2,3,4]. In contrast, exposure to air pollution has been associated with an elevated risk of T2D [5,6,7]. Because PA increases the respiration rate, the intake of polluted air during outdoor PA may increase considerably and potentially intensify the detrimental health effects of air pollution. Therefore, it is of great public health interest to examine the tradeoffs between the health benefits of PA and the intensified harmful effects of air pollution during PA on T2D.

To date, only two studies have investigated the joint effect of PA and air pollution on diabetes [8, 9]. In a study of older Korean adults, Kim et al. [8] found that the benefits of PA on diabetes outweighed the risks related to ambient particulate matter (PM) exposure. This study examined only the weekly frequency of moderate to vigorous PA, which might not fully represent the amount of PA. Another study showed that habitual physical activity can reduce the risk of diabetes regardless of the levels of PM ≤ 2.5 μm (PM2.5) exposure [9]. To our knowledge, no study has examined the combined effects of PA and exposure to traffic-related air pollution, such as nitrogen dioxide (NO2), on the risk of T2D to date.

To address this knowledge gap, we used data from the UK Biobank, a large-scale prospective, population-based cohort study, to analyze associations of the frequency and duration of PA, long-term exposure to air pollution, and their joint effect with the risk of T2D. We assessed exposure to PM2.5, PM ≤ 10 μm (PM10), PM 2.5–10 μm (PMcoarse), and NO2 to comprehensively examine the effects of air pollution on T2D. Furthermore, due to the suggestions of previous studies that genetic variations may modify associations between environmental factors and the risk of T2D [10, 11], we also examined interactions between PA or air pollution and genetic risk for T2D in this study.

Methods

Study design and participants

Between April 2007 and December 2010, the UK Biobank recruited 502,536 participants aged 40–69 years who attended one of 22 assessment centers across England, Wales, and Scotland [12]. Participants completed a touch screen questionnaire, had physical measurements taken, and provided biological samples, as described in detail elsewhere [13]. Ethical approval for the UK Biobank study was obtained from the North West Multi-centre Research Ethics Committee (06/MRE08/65), and all participants provided written informed consent.

We excluded 1299 participants who subsequently withdrew from the study, 132,920 participants who had incomplete information (99,907 with missing PA data and 33,013 with missing information for any of the air pollution exposures), and 9164 participants who had T2D at baseline. Our primary analyses included 359,153 participants. Participants (n = 7968) with missing genetic data were excluded from the genetic analyses. A flowchart of the study sample selection process is presented in Supplementary Fig. 1. There were no significant differences in baseline characteristics between the included individuals and the total population of the UK Biobank (Supplementary Table 1S).

Assessment of physical activity

PA assessment was based on self-reports at baseline and obtained using the International Physical Activity Questionnaire (IPAQ) short form, which includes 6 questions about the frequencies and durations of three types of activities (walking, moderate-intensity activities, and vigorous activities) performed in a typical week [14]. Data were analyzed in accordance with the IPAQ scoring protocol. PA was computed in metabolic equivalent of task minutes per week (MET-min/wk), which could effectively measure the overall PA level of the participants and was then categorized as low (< 600 MET-min/wk), moderate (600–3000 MET-min/wk), or high (> 3000 MET-min/wk) PA [15]. Both categorical and continuous PA variables were used for the data analyses.

Assessment of air pollution

The annual average concentrations of PM2.5, PMcoarse, PM10, and NO2 were calculated using a land use regression (LUR) model developed by the European Study of Cohorts and Air Pollution Effects (ESCAPE) [16] and linked to the geocoded residential addresses of UK Biobank participants. The LUR model calculated the spatial variations in the annual average air pollutant concentrations at the participants’ home addresses, which were provided at the baseline visit, using predictor variables obtained from a geographic information system such as traffic, land use, and topography. Because of the use of a high-resolution European map [17], annual concentration data for PM10 and NO2 were available for several years (2007 and 2010 for PM10 and 2005, 2006, 2007, and 2010 for NO2) in the UK Biobank, therefore, we averaged the obtained values to obtain the air pollutant concentrations of PM10 and NO2. All other particulate matter (PM2.5 and PMcoarse) exposure data were available for a single year only (2010) in the UK Biobank. Participants were divided into 3 categories on the basis of tertile cutoffs for each air pollution concentration: low PM2.5 (< 9.5 µg/m3), moderate PM2.5 (9.5–10.3 µg/m3), and high PM2.5 (> 10.3 µg/m3); low PMcoarse (< 5.9 µg/m3), moderate PMcoarse (5.9–6.4 µg/m3), and high PMcoarse (> 6.4 µg/m3); low PM10 (< 18.4 µg/m3), moderate PM10 (18.4–19.9 µg/m3), and high PM10 (> 19.9 µg/m3); and low NO2 (< 24.6 µg/m3), moderate NO2 (24.6–31.5 µg/m3), and high NO2 (> 31.5 µg/m3). Both categorical and continuous air pollution exposure data were used in the data analyses.

Definition of genetic risk score

We created a genetic risk score (GRS) for T2D using 102 single nucleotide polymorphisms (SNPs), which passed quality control, based on a previous study (Supplementary Table 1S) [18]. A weighted method was used to calculate the T2D GRS. We calculated the sum of the number of associated alleles (0, 1, or 2); each SNP was weighted by the strength of its association with T2D in a previous genome-wide association study [19]. The T2D GRSs showed a normal distribution (Supplementary Fig. 2S), and a higher T2D GRS indicated a higher genetic predisposition to T2D. The T2D GRS was then divided into tertiles to stratify individuals into high, intermediate, and low genetic risk categories. Detailed information about genotyping, imputation, and quality control in the UK Biobank study has been described previously [20].

Ascertainment of type 2 diabetes

The prevalence and incidence of diabetes were assessed based on the UK Biobank algorithms for the diagnosis of diabetes [21]. Incident T2D was ascertained using hospital inpatient records containing data on admissions and diagnoses obtained from the Hospital Episode Statistics for England, Scottish Morbidity Record data for Scotland, and the Patient Episode Database for Wales [22]. Participants with T2D were defined by the International Classification of Diseases, 10th revision (ICD-10) code E11. Follow-up time was defined as the time from the date of attendance until the date of first diagnosis, February 28, 2017, for Scotland, or February 25, 2018, for Wales and England, whichever occurred first. Detailed information on the ascertainment of T2D is available online at http://content.digital.nhs.uk/services.

Assessment of covariates

Potential confounders were selected according to the previously published literature [8, 23, 24]. We used the baseline touch screen questionnaire to assess several potential confounders: age, sex, race, education, household income, smoking status, alcohol consumption, body mass index (BMI), vegetable intake, fruit intake, family history of diabetes, comorbidities (hypertension, cardiovascular disease [CVD], depression, and cancer), genotyping chip and first 10 principal components of ancestry. BMI was calculated by dividing a participant’s weight by the square of his or her height in meters (kg/m2). Hypertension was defined as a self-reported history of hypertension, a systolic blood pressure ≥ 140 mmHg, a diastolic blood pressure ≥ 90 mmHg, or antihypertensive medication use. Supplementary Table 3S includes the coding of the variables under the assessment of the covariates. Details of these measurements are provided on the website of the UK Biobank (www.ukbiobank.ac.uk).

Statistical analyses

Baseline characteristics of the participants are summarized across T2D status as means (standard deviations [SDs]) for continuous variables and numbers (percentages) for categorical variables. To impute the missing covariate values (all covariates had < 3% of the missing values), we used multiple imputation by chained equations using the R package “mice” to impute the missing covariate values [25]. Cox proportional hazard models were constructed to calculate hazard ratios (HRs) and 95% confidence intervals (CIs) for associations of PA, air pollution (PM2.5, PMcoarse, PM10, or NO2), genetic risk, and their combination with incident T2D. The proportional hazards assumption, assessed using the Schoenfeld residuals technique [26], was satisfied.

We ran three models that included an increasing number of covariates: Model 1 included age and sex (men or women); Model 2 (multivariable-adjusted model) was adjusted as in Model 1 but also included race (white, Asian, black, Chinese, mixed, or other race), education (degree or no degree), household income (<£18,000, £18,000-£30,999, £31,000-£51,999, £52,000-£100,000, or >£100,000), smoking status (current, former, or never), alcohol consumption (current, former, or never), BMI, fruit and vegetable intake (< 2 or ≥ 2 servings per day), family history of diabetes (yes or no), hypertension (yes or no), CVD (yes or no), depression (yes or no), and cancer (yes or no); Model 3 was further mutually adjusted for PA and air pollution. PA and each air pollution variable were treated as continuous variables, and HRs were calculated per 600 MET-min/wk difference in PA, per 5 µg/m3 difference in PM2.5 and PMcoarse, and per 10 µg/m3 difference in PM10 and NO2. Stratified analyses were conducted to examine the associations of PA in each air pollution stratum. To investigate the joint associations of PA and air pollution with T2D, the participants were then classified into 9 groups according to the categories of PA and each individual air pollutant with reference to the participants with low PA and high air pollutant exposure levels. For analyses of the genetic data, Model 3 was additionally adjusted for the genotyping chip and the first 10 principal components of ancestry. Interactions between PA, individual air pollutants and T2D GRS were assessed with the likelihood ratio test.

To examine the robustness of the primary findings, we conducted a series of sensitivity analyses. First, to examine the possibility of reverse causation bias, we excluded participants who developed T2D within the first two years of follow-up. Second, we excluded participants who had T2D-related diseases (CVD, cancer and hypertension) to eliminate potential comorbidity effects. Third, we restricted the analyses to participants with no missing covariate data. Fourth, employment status (working, retired, unemployed, other) of participants was further adjusted in the models. Fifth, participants were divided into 2 categories based on the WHO air quality guidelines for each air pollutant to examine the associations air pollution with incident T2D. Finally, we restricted the analyses to participants of European ancestry to test the association between T2D -GRS and incident T2D.

All analyses were performed using R software, version 4.0.2 (R Development Core Team, Vienna, Austria). A p value less than 0.05 was considered statistically significant in all analyses.

Results

Baseline characteristics

The baseline characteristics of the included participants (n = 359,153) are provided in Table 1. Overall, participants had a mean (SD) age of 56.3 (8.1) years, and 52.8% were women. Most participants (50.4%) performed a moderate volume of physical activity (600–3000 MET-min/wk). Participants who engaged in more PA were more likely to be former smokers and have a lower prevalence of hypertension, CVD, and depression. The mean (SD) annual average concentrations of PM2.5, PMcoarse, PM10, and NO2 were 9.97 (1.06) µg/m3, 6.42 (0.90) µg/m3, 19.29 (1.96) µg/m3, and 29.14 (9.30) µg/m3, respectively. The pearson correlation coefficients of the air pollutants were shown in Supplementary Table 4S.

Table 1 Baseline characteristics of participants

Associations of PA and air pollution with incident T2D

During a median of 8.9 (interquartile range: 8.2–9.5) years of follow-up, 13,706 T2D events were recorded. Table 2 shows the associations of PA and individual air pollution variables with incident T2D. We found that a higher level of PA was associated with a decreased risk of T2D after adjusting for a series of covariates, including air pollution. Compared with the low PA group, the moderate PA and high PA groups had adjusted HRs of 0.82 (95% CI, 0.79–0.86) and 0.73 (95% CI, 0.70–0.77), respectively. In contrast, higher air pollution levels were associated with an increased risk of T2D after adjusting for a series of covariates, including PA. Compared with the low air pollution group, the adjusted HRs for T2D in the moderate and high air pollution groups were 1.07 (95% CI,1.02–1.11) and 1.19 (95% CI,1.14–1.24) for PM2.5, 1.02 (95% CI, 0.98–1.07) and 1.06 (95% CI, 1.02–1.11) for PMcoarse, 1.06 (95% CI, 1.02–1.11) and 1.13 (95% CI, 1.08–1.18) for PM10, and 1.08 (95% CI, 1.04–1.13) and 1.19 (95% CI, 1.14–1.24) for NO2, respectively. Moreover, we found significant trends for the associations of incident T2D across the categories of PA and all air pollution variables (Table 2).

Table 2 Association of physical activity or air pollution with risk of incident type 2 diabetes

Joint effect of PA and air pollution on T2D

Table 3 indicates the associations between PA and incident T2D stratified by individual air pollution variables. Stratified analyses indicated that PA was inversely associated with the risk of T2D in each air pollution group. High levels of PA were associated with a 25-28% lower risk for T2D (HR between 0.72 and 0.75) than low PA levels at different levels of each air pollutant. No significant interactions between PA and air pollution were observed (Table 3 and Supplementary Table 5S; all P-interactions > 0.05).

Table 3 Risk of incident type 2 diabetes according to physical activity category within each of air pollution category

Figure 1 presents the joint associations of PA and air pollution with the risk of T2D. The analyses indicated that participants in the high PA and low air pollution groups had the lowest risk of T2D. Compared to participants in the low PA and high air pollution groups, the HR for T2D among participants with high PA in the low PM2.5 group was 0.61 (95% CI: 0.57–0.67), in the low PMcoarse group was 0.69 (95% CI: 0.64–0.75), in the low PM10 group was 0.66 (95% CI: 0.61–0.72), and in the low NO2 group was 0.63 (95% CI: 0.58–0.69). The inverse associations between PA and T2D in each air pollution stratum were generally consistent (all P for trend < 0.05).

Fig. 1
figure 1

Joint associations of physical activity and PM2.5 (A), PM10 (B), PMcoarse (C), and NO2 (D) with the incidence of type 2 diabetes. Abbreviations: CI, confidence interval; HR, hazard ratio; PA, physical activity. The results were obtained after adjusting for age, sex, race, education, household income, smoking status, alcohol consumption, body mass index, fruit and vegetable intake, family history of diabetes, hypertension, cardiovascular disease, depression, and cancer

The sensitivity analyses showed no substantial changes when we excluded participants who developed T2D during the first two years of follow-up (Supplementary Fig. 3S), those for whom covariate data were missing (Supplementary Fig. 4S), those who had T2D-related diseases (CVD, cancer and hypertension) (Supplementary Fig. 5S), further adjusted the employment status in the models (Supplementary Fig. 6S), or divided participants according to the WHO air quality guidelines (Supplementary Table 6S).

Joint effect of PA or air pollution and T2D GRS on T2D

Genetic data were available for 351,185 participants in this study. In the multivariable-adjusted model, compared with participants with a low T2D GRS, those with an intermediate (HR: 1.45, 95% CI: 1.38–1.52) or high (HR: 2.20, 95% CI: 2.10–2.30) T2D GRS had an increased risk of T2D (Supplementary Table 7S). No substantial changes of the associations when we restricted the analyses to participants of European ancestry (Supplementary Table 7S). In the joint effect analyses, PA and each air pollution variable were significantly associated with the risk of T2D independent of T2D GRS (Fig. 2). There was no significant interaction between PA or each type of air pollutant and T2D GRS (all P-interactions > 0.05), indicating that the associations with PA and each type of air pollutant did not vary substantially on the basis of genetic risk.

Fig. 2
figure 2

Joint associations of physical activity or air pollution and genetic risk with the incidence of type 2 diabetes. Abbreviations: CI, confidence interval; HR, hazard ratio. The results were obtained after adjusting for age, sex, race, education, household income, smoking status, alcohol consumption, body mass index, fruit and vegetable intake, family history of diabetes, hypertension, cardiovascular disease, depression, cancer, genotyping, the first 10 genetic principal components, and air pollution or physical activity. Individuals in the low PA or low air pollution group were used as the reference group (*). P-interaction describes the interactions between PA or air pollution and the genetic risk of type 2 diabetes. The interactions between PA or air pollution and T2D GRS were assessed with the likelihood ratio test by including an interaction term between PA (low, moderate, and high PA were coded as 0, 1 and 2, respectively) or air pollution (low, moderate, and high air pollution were coded as 0, 1 and 2, respectively) and T2D GRS (low, intermediate, and high PA were coded as 0, 1 and 2, respectively) in the multivariable-adjusted model

Discussion

In this large population-based cohort study involving 359,153 individuals, we identified the following important findings: (1) higher PA and lower air pollution exposure (PM2.5, PMcoarse, PM10, and NO2) levels were associated with a lower risk of T2D after adjusting for a series of covariates, including a mutual adjustment for PA or air pollution; (2) the inverse association between PA and incident T2D generally remained stable among participants exposed to different levels of each air pollutant (all P-interactions > 0.05); and (3) there was no significant interaction between PA or air pollution and genetic risk (all P-interactions > 0.05), and higher PA and lower air pollution exposure levels were associated with a lower risk of T2D regardless of genetic risk.

Our significant finding of a low risk of developing T2D related to PA was consistent with existing evidence [27, 28]. Notably, a recent meta-analysis summarized 55 cohort studies and reported a 28% lower risk of T2D among participants with high PA levels than among participants with low PA levels [27]; similarly, our study indicated a 27% (HR: 0.73, 95% CI: 0.70–0.77) lower risk of T2D among participants with high PA levels than among those with low PA levels. Several mechanisms might explain the effect of PA on T2D. For example, PA might reduce the risk of T2D by increasing cardiorespiratory fitness, improving lipid levels and endothelial function [29], reducing glycosylated hemoglobin (HbA1c) levels and improving insulin sensitivity [30]. The activation of anti-inflammatory signaling pathways may be another potential mechanism underlying the effect of PA on T2D [31, 32]. In addition, we found that long-term exposure of individuals to air pollutants, including PM2.5, PMcoarse, PM10, and NO2, was associated with a higher risk of T2D, providing further evidence of a positive association between long-term exposure to air pollution and T2D [5,6,7, 33, 34]. The adverse effects of air pollution on HbA1c and fasting glucose concentrations have been well documented in previous studies [33]. In addition, air pollution was associated with increased levels of systemic inflammation [35, 36] and oxidative stress [37], which may increase the risk of T2D [38, 39]. Overall, our study further supported the need to establish measures to increase PA and tackle air pollution, which might contribute to reducing the burden of diabetes.

Whether the health benefits of PA are moderated by long-term exposure to air pollution remains in dispute. Notably, a few studies found significant interactions between air pollution and PA, including evidence that air pollution counteracted the benefits of PA on asthma [40] and stroke [41]. In contrast, several studies revealed no significant interactions between PA and long-term exposure to air pollution on lung function/respiratory diseases [42, 43], hypertension [24], myocardial infarction [44], and mortality [45]. Other studies have also shown that although PA did not interact with air pollution, it counteracted the hazardous effect of air pollution on blood pressure [46] and atherosclerotic cardiovascular disease [47]. This discrepancy may be partially attributed to the differences in sample sizes, study durations, and the levels of air pollution in the study regions[24]. Moreover, recent studies found that the benefits of PA for diabetes did not interact with fine particulate matter and outweighed the risks related to ambient particulate matter exposure [8, 9]. Consistent with those findings, we found in the current study that the long-term benefits of PA with respect to T2D were not significantly moderated by exposure to air pollution. Furthermore, our finding that the long-term benefits of PA for T2D were not moderated by exposure to high levels of NO2 is novel. The current study, therefore, may indicate that the effects of long-term exposure to air pollution and PA on T2D are independent of each other, with the benefits of PA not being reduced by exposure to air pollution, including PM2.5, PMcoarse, PM10, and NO2.

Interestingly, our study found that the associations of PA with T2D remained consistent among participants exposed to different individual air pollution levels, indicating that PA may decrease the risk of T2D, be it among people with relatively high or low levels of air pollution exposure. However, the exact reasons for the stable protective effect of PA on T2D regardless of air pollution have not been clarified. One potential hypothesis to explain this observation is that the additional air pollutants inhaled during PA account for only a small fraction of the total inhaled dose of air pollution [48] and, therefore, are not sufficient to increase the risk of T2D. It is also possible that health benefits due to increased PA levels generally outweigh the risks related to increases in inhaled air pollution doses during physical activity or exercise [49]. Thus, our study highlights an increase in PA as a potentially effective measure for the prevention of the incidence and progression of T2D, regardless of air pollution exposure.

Consistent with previous studies [50, 51], the T2D GRS was significantly positively associated with the risk of T2D in the current study. However, we did not observe a significant interaction between PA or air pollution and T2D GRS on the risk of T2D, suggesting that higher PA levels or reducing air pollution exposure may protect against T2D, regardless of the genetic risk profile.

Strengths and limitations

The strengths of this study included a large sample size and outcome events, the inclusion of information on a wealth of potential confounders, highly accurate diabetes diagnoses by UK Biobank algorithms, confirmation of T2D events by medical record review, and uniform data collection protocols, which reduced measurement error and thus reduced biased estimates. Another major novelty of the current study is the examination of the interactions of PA and air pollution with genetic risk.

However, there were several limitations in this study that should be noted. First, we did not distinguish between indoor and outdoor PA, and indoor air pollution and relocation of participants during the study period were not considered in this study, which may contribute to exposure misclassification. Second, the assessment of exposure to air pollution relied on residence locations and did not completely reflect personal exposure. Given the inherent imprecision of the spatially derived exposure levels, the air pollution assessment was likely subject to nondifferential misclassification that may have attenuated our results. Third, air pollution exposure in our study was available at only certain periods in time. Although we made a reasonable assumption that the spatial contrast in air pollution exposure was relatively stable in the UK over these years, the possibility of exposure misclassification cannot be excluded. Fourth, our study was conducted in a low-pollution area. Further studies in areas with moderate and severe pollution are needed to examine the applicability of our findings. Fifth, although this study tried to identify incident T2D based on hospital inpatient records, there are no biomarkers (e.g., HbA1c and fasting blood glucose) to assess the status of T2D. Furthermore, in observational studies, the possibility of residual confounding factors due to imprecise measurements or unknown factors cannot be excluded. Therefore, although we carefully adjusted for various confounders in our analyses, the associations might have been affected by unknown factors. Finally, reverse causality might exist in our study, although the results remained unchanged when we excluded participants with outcome events that occurred during the first two years of follow-up.

Conclusion

In summary, the results of this large-scale prospective cohort study showed that higher PA and lower air pollution levels were independently associated with a lower risk of T2D, and the benefits of PA for T2D generally remained stable among participants exposed to different levels of each air pollution variable, including PM2.5, PMcoarse, PM10, and NO2. Our findings suggested that PA should be promoted to prevent T2D among people with both relatively high and low levels of air pollution exposure. Further studies are needed to validate our findings in regions with moderate and severe pollution levels.

Availability of data and materials

The UK Biobank data are available from the UK Biobank upon request (www.ukbiobank.ac.uk/).

Abbreviations

BMI:

Body mass index

CIs:

Confidence intervals

CVD:

Cardiovascular disease

ESCAPE:

European Study of Cohorts and Air Pollution Effects

GRS:

Genetic risk score

HbA1c:

Glycosylated hemoglobin

HRs:

Hazard ratios

ICD-10:

International Classification of Diseases, 10th revision

IPAQ:

International Physical Activity Questionnaire

LUR:

Land use regression

MET:

Metabolic equivalent of task

NO2:

Nitrogen dioxide

PA:

Physical activity

PM:

Particulate matter

SD:

Standard deviation

SNPs:

Single nucleotide polymorphisms

T2D:

Type 2 diabetes

References

  1. Lascar N, Brown J, Pattison H, Barnett AH, Bailey CJ, Bellary S. Type 2 diabetes in adolescents and young adults. Lancet Diabetes Endocrinol. 2018;6(1):69–80.

    Article  Google Scholar 

  2. Zheng Y, Ley SH, Hu FB. Global aetiology and epidemiology of type 2 diabetes mellitus and its complications. Nat Rev Endocrinol. 2018;14(2):88–98.

    Article  Google Scholar 

  3. Balducci S, D’Errico V, Haxhi J, Sacchetti M, Orlando G, Cardelli P, Vitale M, Bollanti L, Conti F, Zanuso S, et al. Effect of a Behavioral Intervention Strategy on Sustained Change in Physical Activity and Sedentary Behavior in Patients With Type 2 Diabetes: The IDES_2 Randomized Clinical Trial. JAMA. 2019;321(9):880–90.

    Article  Google Scholar 

  4. Hu G, Jousilahti P, Barengo NC, Qiao Q, Lakka TA, Tuomilehto J. Physical activity, cardiovascular risk factors, and mortality among Finnish adults with diabetes. Diabetes Care. 2005;28(4):799–805.

    Article  Google Scholar 

  5. Andersen ZJ, Raaschou-Nielsen O, Ketzel M, Jensen SS, Hvidberg M, Loft S, Tjønneland A, Overvad K, Sørensen M. Diabetes incidence and long-term exposure to air pollution: a cohort study. Diabetes Care. 2012;35(1):92–8.

    Article  CAS  Google Scholar 

  6. Coogan PF, White LF, Jerrett M, Brook RD, Su JG, Seto E, Burnett R, Palmer JR, Rosenberg L. Air pollution and incidence of hypertension and diabetes mellitus in black women living in Los Angeles. Circulation. 2012;125(6):767–72.

    Article  CAS  Google Scholar 

  7. GBD 2019 Risk Factors Collaborators. Global burden of 87 risk factors in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet (London, England). 2020, 396(10258):1223–1249.

  8. Kim SR, Choi D, Choi S, Kim K, Lee G, Son JS, Kim KH, Park SM. Association of combined effects of physical activity and air pollution with diabetes in older adults. Environ Int. 2020;145:106161.

    Article  CAS  Google Scholar 

  9. Guo C, Yang HT, Chang LY, Bo Y, Lin C, Zeng Y, Tam T, Lau AKH, Hoek G, Lao XQ. Habitual exercise is associated with reduced risk of diabetes regardless of air pollution: a longitudinal cohort study. Diabetologia. 2021;64(6):1298–308.

    Article  Google Scholar 

  10. Eze IC, Imboden M, Kumar A, von Eckardstein A, Stolz D, Gerbase MW, Künzli N, Pons M, Kronenberg F, Schindler C, et al. Air pollution and diabetes association: modification by type 2 diabetes genetic risk score. Environ Int. 2016;94:263–71.

    Article  CAS  Google Scholar 

  11. Zanuso S, Sacchetti M, Sundberg CJ, Orlando G, Benvenuti P, Balducci S. Exercise in type 2 diabetes: genetic, metabolic and neuromuscular adaptations. A review of the evidence. Br J Sports Med. 2017;51(21):1533–8.

    Article  Google Scholar 

  12. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, Downey P, Elliott P, Green J, Landray M. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779.

    Article  Google Scholar 

  13. Palmer LJ. UK Biobank: bank on it. Lancet (London England). 2007;369(9578):1980–2.

    Article  Google Scholar 

  14. Celis-Morales CA, Lyall DM, Anderson J, Iliodromiti S, Fan Y, Ntuk UE, Mackay DF, Pell JP, Sattar N, Gill JM. The association between physical activity and risk of mortality is modulated by grip strength and cardiorespiratory fitness: evidence from 498 135 UK-Biobank participants. Eur Heart J. 2017;38(2):116–22.

    Google Scholar 

  15. Lear SA, Hu W, Rangarajan S, Gasevic D, Leong D, Iqbal R, Casanova A, Swaminathan S, Anjana RM, Kumar R, et al. The effect of physical activity on mortality and cardiovascular disease in 130,000 people from 17 high-income, middle-income, and low-income countries: the PURE study. Lancet (London, England). 2017;390(10113):2643–54.

    Article  Google Scholar 

  16. Beelen R, Raaschou-Nielsen O, Stafoggia M, Andersen ZJ, Weinmayr G, Hoffmann B, Wolf K, Samoli E, Fischer P, Nieuwenhuijsen M, et al. Effects of long-term exposure to air pollution on natural-cause mortality: an analysis of 22 European cohorts within the multicentre ESCAPE project. Lancet (London England). 2014;383(9919):785–95.

    Article  CAS  Google Scholar 

  17. Vienneau D, de Hoogh K, Bechle MJ, Beelen R, van Donkelaar A, Martin RV, Millet DB, Hoek G, Marshall JD. Western European land use regression incorporating satellite- and ground-based measurements of NO2 and PM10. Environ Sci Technol. 2013;47(23):13555–64.

    Article  CAS  Google Scholar 

  18. Scott RA, Scott LJ, Mägi R, Marullo L, Gaulton KJ, Kaakinen M, Pervjakova N, Pers TH, Johnson AD, Eicher JD, et al. An Expanded Genome-Wide Association Study of Type 2 Diabetes in Europeans. Diabetes. 2017;66(11):2888–902.

    Article  CAS  Google Scholar 

  19. Vetter C, Dashti HS, Lane JM, Anderson SG, Schernhammer ES, Rutter MK, Saxena R, Scheer F. Night Shift Work, Genetic Risk, and Type 2 Diabetes in the UK Biobank. Diabetes Care. 2018;41(4):762–9.

    Article  Google Scholar 

  20. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, Motyer A, Vukcevic D, Delaneau O, O’Connell J, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562(7726):203–9.

    Article  CAS  Google Scholar 

  21. Eastwood SV, Mathur R, Atkinson M, Brophy S, Sudlow C, Flaig R, de Lusignan S, Allen N, Chaturvedi N. Algorithms for the Capture and Adjudication of Prevalent and Incident Diabetes in UK Biobank. PLoS One. 2016;11(9):e0162388.

    Article  Google Scholar 

  22. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, Downey P, Elliott P, Green J, Landray M, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779.

    Article  Google Scholar 

  23. Howell NA, Tu JV, Moineddin R, Chen H, Chu A, Hystad P, Booth GL. Interaction between neighborhood walkability and traffic-related air pollution on hypertension and diabetes: The CANHEART cohort. Environ Int. 2019;132:104799.

    Article  CAS  Google Scholar 

  24. Guo C, Zeng Y, Chang LY, Yu Z, Bo Y, Lin C, Lau AK, Tam T, Lao XQ. Independent and Opposing Associations of Habitual Exercise and Chronic PM(2.5) Exposures on Hypertension Incidence. Circulation. 2020;142(7):645–56.

    Article  CAS  Google Scholar 

  25. Schafer JL. Multiple imputation: a primer. Stat Methods Med Res. 1999;8(1):3–15.

  26. Grambsch PM, Therneau TM. Proportional hazards tests and diagnostics based on weighted residuals. Biometrika. 1994;81(3):515–26.

    Article  Google Scholar 

  27. Kyu HH, Bachman VF, Alexander LT, Mumford JE, Afshin A, Estep K, Veerman JL, Delwiche K, Iannarone ML, Moyer ML, et al. Physical activity and risk of breast cancer, colon cancer, diabetes, ischemic heart disease, and ischemic stroke events: systematic review and dose-response meta-analysis for the Global Burden of Disease Study 2013. BMJ (Clinical research ed). 2016;354:i3857.

    Google Scholar 

  28. Aune D, Norat T, Leitzmann M, Tonstad S, Vatten LJ. Physical activity and the risk of type 2 diabetes: a systematic review and dose-response meta-analysis. Eur J Epidemiol. 2015;30(7):529–42.

    Article  Google Scholar 

  29. Colberg SR, Sigal RJ, Yardley JE, Riddell MC, Dunstan DW, Dempsey PC, Horton ES, Castorino K, Tate DF. Physical Activity/Exercise and Diabetes: A Position Statement of the American Diabetes Association. Diabetes Care. 2016;39(11):2065–79.

    Article  Google Scholar 

  30. Arakawa K. Antihypertensive mechanism of exercise. J Hypertens. 1993;11(7):H45.

    Article  Google Scholar 

  31. Donath MY, Shoelson SE. Type 2 diabetes as an inflammatory disease. Nat Rev Immunol. 2011;11(2):98–107.

    Article  CAS  Google Scholar 

  32. Donath MY. Targeting inflammation in the treatment of type 2 diabetes: time to start. Nat Rev Drug Discovery. 2014;13(6):465–76.

    Article  CAS  Google Scholar 

  33. Thiering E, Heinrich J. Epidemiology of air pollution and diabetes. Trends Endocrinol Metab. 2015;26(7):384–94.

    Article  CAS  Google Scholar 

  34. Pearson JF, Bachireddy C, Shyamprasad S, Goldfine AB, Brownstein JS. Association between fine particulate matter and diabetes prevalence in the U. S. Diabetes Care. 2010;33(10):2196–201.

    Article  Google Scholar 

  35. Pope CA 3rd, Bhatnagar A, McCracken JP, Abplanalp W, Conklin DJ, O’Toole T. Exposure to Fine Particulate Air Pollution Is Associated With Endothelial Injury and Systemic Inflammation. Circul Res. 2016;119(11):1204–14.

    Article  CAS  Google Scholar 

  36. Li H, Cai J, Chen R, Zhao Z, Ying Z, Wang L, Chen J, Hao K, Kinney PL, Chen H, et al. Particulate Matter Exposure and Stress Hormone Levels: A Randomized, Double-Blind, Crossover Trial of Air Purification. Circulation. 2017;136(7):618–27.

    Article  CAS  Google Scholar 

  37. Chuang KJ, Chan CC, Su TC, Lee CT, Tang CS. The effect of urban air pollution on inflammation, oxidative stress, coagulation, and autonomic dysfunction in young adults. Am J Respir Crit Care Med. 2007;176(4):370–6.

    Article  CAS  Google Scholar 

  38. Odegaard AO, Jacobs DR Jr, Sanchez OA, Goff DC Jr, Reiner AP, Gross MD. Oxidative stress, inflammation, endothelial dysfunction and incidence of type 2 diabetes. Cardiovasc Diabetol. 2016;15:51.

    Article  Google Scholar 

  39. Brooks-Worrell BM, Palmer JP. Setting the Stage for Islet Autoimmunity in Type 2 Diabetes: Obesity-Associated Chronic Systemic Inflammation and Endoplasmic Reticulum (ER) Stress. Diabetes Care. 2019;42(12):2338–46.

    Article  CAS  Google Scholar 

  40. Paulin LM, Hansel NN. Physical Activity and Air Pollution Exposures in the Urban Environment. Am J Respir Crit Care Med. 2016;194(7):786–7.

    Article  Google Scholar 

  41. Lin H, Guo Y, Di Q, Zheng Y, Kowal P, Xiao J, Liu T, Li X, Zeng W, Howard SW, et al. Ambient PM(2.5) and Stroke: Effect Modifiers and Population Attributable Risk in Six Low- and Middle-Income Countries. Stroke. 2017;48(5):1191–7.

    Article  CAS  Google Scholar 

  42. Fisher JE, Loft S, Ulrik CS, Raaschou-Nielsen O, Hertel O, Tjønneland A, Overvad K, Nieuwenhuijsen MJ, Andersen ZJ. Physical Activity, Air Pollution, and the Risk of Asthma and Chronic Obstructive Pulmonary Disease. Am J Respir Crit Care Med. 2016;194(7):855–65.

    Article  CAS  Google Scholar 

  43. Fuertes E, Markevych I, Jarvis D, Vienneau D, de Hoogh K, Antó JM, Bowatte G, Bono R, Corsico AG, Emtner M, et al. Residential air pollution does not modify the positive association between physical activity and lung function in current smokers in the ECRHS study. Environ Int. 2018;120:364–72.

    Article  CAS  Google Scholar 

  44. Kubesch NJ, Therming Jørgensen J, Hoffmann B, Loft S, Nieuwenhuijsen MJ, Raaschou-Nielsen O, Pedersen M, Hertel O, Overvad K, Tjønneland A, et al. Effects of leisure-time and transport-related physical activities on the risk of incident and recurrent myocardial infarction and interaction with traffic-related air pollution: a cohort study. J Am Heart Assoc. 2018;7(15).

    Article  CAS  Google Scholar 

  45. Sun S, Cao W, Qiu H, Ran J, Lin H, Shen C, Siu-Yin Lee R, Tian L. Benefits of physical activity not affected by air pollution: a prospective cohort study. Int J Epidemiol. 2020;49(1):142–52.

    Article  Google Scholar 

  46. Avila-Palencia I, Laeremans M, Hoffmann B, Anaya-Boig E, Carrasco-Turigas G, Cole-Hunter T, de Nazelle A, Dons E, Götschi T, Int Panis L, et al. Effects of physical activity and air pollution on blood pressure. Environ Res. 2019;173:387–96.

    Article  CAS  Google Scholar 

  47. Tu R, Hou J, Liu X, Li R, Dong X, Pan M, Mao Z, Huo W, Chen G, Guo Y, et al. Physical activity attenuated association of air pollution with estimated 10-year atherosclerotic cardiovascular disease risk in a large rural Chinese adult population: A cross-sectional study. Environ Int. 2020;140:105819.

    Article  CAS  Google Scholar 

  48. Rojas-Rueda D, de Nazelle A, Tainio M, Nieuwenhuijsen MJ. The health risks and benefits of cycling in urban environments compared with car use: health impact assessment study. BMJ (Clinical research ed). 2011;343:d4521.

    Article  Google Scholar 

  49. Andersen ZJ, de Nazelle A, Mendez MA, Garcia-Aymerich J, Hertel O, Tjønneland A, Overvad K, Raaschou-Nielsen O, Nieuwenhuijsen MJ. A study of the combined effects of physical activity and air pollution on mortality in elderly urban residents: the Danish Diet, Cancer, and Health Cohort. Environ Health Perspect. 2015;123(6):557–63.

    Article  Google Scholar 

  50. Langenberg C, Sharp SJ, Franks PW, Scott RA, Deloukas P, Forouhi NG, Froguel P, Groop LC, Hansen T, Palla L, et al. Gene-lifestyle interaction and type 2 diabetes: the EPIC interact case-cohort study. PLoS Med. 2014;11(5):e1001647.

    Article  Google Scholar 

  51. Ma H, Li X, Zhou T, Sun D, Liang Z, Li Y, Heianza Y, Qi L. Glucosamine Use, Inflammation, and Genetic Susceptibility, and Incidence of Type 2 Diabetes: A Prospective Study in UK Biobank. Diabetes Care. 2020;43(4):719–25.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the UK Biobank participants. This research was conducted using the UK Biobank resource under application number 43795.

Funding

This work was supported by the National Natural Science Foundation of China (82204115, 82103931, 82003443, and 81973109), the Project Supported by the Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme (2019), the Construction of High-level University of Guangdong (G820332010, G618339167, and G618339164), the Guangdong Basic and Applied Basic Research Foundation (2022A1515012085, 2021A1515011629, and 2021A1515110230), the Guangzhou Science and Technology Project (202002030255), and the Young Elite Scientists Sponsorship Program by CAST (2021QNRC001 and 2019QNRCO01). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

Z-HL and W-FZ are joint first authors, contributed to the statistical analyses, and had primary responsibility for writing the manuscript. Z-HL and W-FZ contributed equally to this article. CM, VBK, and XG directed the study. XRZ, VCC, W-QS, QC, X-MW, Q-MH, XC, H-LY, and M-CC contributed to the data cleaning. CM, DS, P-DZ, DL, Y-JZ, and P-LC contributed to the analysis or interpretation of the data. All authors critically reviewed the manuscript for important intellectual content. CM is the study guarantor. The corresponding author (CM) attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Chen Mao.

Ethics declarations

Ethics approval and consent to participate

The UK Biobank received ethical approval from the research ethics committee (REC reference for UK Biobank 11/NW/0382), and participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: 

Table 1S. Distribution of population characteristics in included individuals and the total population of UK Biobank. Table 2S. Single-nucleotide polymorphisms used to build the genetic risk score for type 2 diabetes. Table 3S. The coding of the variables under assessment of the covariates in the UK Biobank study. Table 4S. Pearson correlation coefficients among the five air pollutants. Table 5S. Risk of incident type 2 diabetes according to air pollution within each of physical activity category. Table 6S. Association of air pollution with risk of incident type 2 diabetes according to WHO air quality guidelines. Table 7S. Hazard ratios for type 2 diabetes across tertiles categories of type 2 diabetes genetic risk score. Figure 1S. Flow chart of participants enrolment. Figure 2S. Distribution of the polygenic risk score for type 2 diabetes. Figure 3S. Joint associations of physical activity and air pollutants with incidence of type 2 diabetes after excluding participants with type 2 diabetes within 2 years of baseline. Figure 4S. Joint associations of physical activity and air pollutants with incidence of type 2 diabetes after excluding participants with missing data for covariates.Figure 5S. Joint associations of physical activity and air pollutants with incidence of type 2 diabetes after excluding participants who had type 2 diabetes related diseases. Figure 6S. Joint associations of physical activity and air pollutants with incidence of type 2 diabetes after adjusting the employment status. Figure 7S. Joint associations of physical activity or air pollution and genetic risk with the incidence of type 2 diabetes after excluding participants of non-European ancestry.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, ZH., Zhong, WF., Zhang, XR. et al. Association of physical activity and air pollution exposure with the risk of type 2 diabetes: a large population-based prospective cohort study. Environ Health 21, 106 (2022). https://doi.org/10.1186/s12940-022-00922-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12940-022-00922-3

Keywords