Risk of Incident Diabetes in Relation to Long-term Exposure to Fine Particulate Matter in Ontario, Canada

Background: Laboratory studies suggest that fine particulate matter (≤ 2.5 µm in diameter; PM2.5) can activate pathophysiological responses that may induce insulin resistance and type 2 diabetes. However, epidemiological evidence relating PM2.5 and diabetes is sparse, particularly for incident diabetes. Objectives: We conducted a population-based cohort study to determine whether long-term exposure to ambient PM2.5 is associated with incident diabetes. Methods: We assembled a cohort of 62,012 nondiabetic adults who lived in Ontario, Canada, and completed one of five population-based health surveys between 1996 and 2005. Follow-up extended until 31 December 2010. Incident diabetes diagnosed between 1996 and 2010 was ascertained using the Ontario Diabetes Database, a validated registry of persons diagnosed with diabetes (sensitivity = 86%, specificity = 97%). Six-year average concentrations of PM2.5 at the postal codes of baseline residences were derived from satellite observations. We used Cox proportional hazards models to estimate the associations, adjusting for various individual-level risk factors and contextual covariates such as smoking, body mass index, physical activity, and neighborhood-level household income. We also conducted multiple sensitivity analyses. In addition, we examined effect modification for selected comorbidities and sociodemographic characteristics. Results: There were 6,310 incident cases of diabetes over 484,644 total person-years of follow-up. The adjusted hazard ratio for a 10-µg/m3 increase in PM2.5 was 1.11 (95% CI: 1.02, 1.21). Estimated associations were comparable among all sensitivity analyses. We did not find strong evidence of effect modification by comorbidities or sociodemographic covariates. Conclusions: This study suggests that long-term exposure to PM2.5 may contribute to the development of diabetes.

Diabetes mellitus and its associated macrovascular (Kannel and McGee 1979;Manson et al. 1991a) and microvascular (Watkins and Edmonds 1983) complications constitute a serious threat to global human health and welfare (Wild et al. 2004). The burden of diabetes relates particularly to type 2 diabetes, which accounts for 90-95% of cases globally (American Diabetes Association 2012). Although obesity (Colditz et al. 1995), diet (Hu et al. 2001), and physical inactivity (Manson et al. 1991b) have been identified as important risk factors for type 2 diabetes, there is increasing evidence that insulin resistance, the underlying hallmark and pathophysiological mechanism leading to type 2 diabetes, can be aggravated by factors that promote inflammatory responses (Hotamisligil 2006).
It is only recently that ambient air pollution has been implicated in the etiology of type 2 diabetes (Brook et al. 2008). Exposure to air pollution has been associated with cardiovascular-related mortality and morbidity (Brook et al. 2010;Chen et al. 2008). There is also evidence that persons with diabetes are particularly susceptible to the acute effects of air pollution (Goldberg et al. 2006;O'Neill et al. 2005). In a mouse model, exposure to fine particulate matter (particles with an aerodynamic diameter ≤ 2.5 µm; PM 2.5 ) increased blood glucose and induced adipose inflammation and insulin resistance (Sun et al. 2009). This animal study provides a potential biological basis for the link between air pollution and diabetes.
Because of the ubiquitous nature of exposure to air pollution, even a modest effect of air pollution on increasing the risk of diabetes may pose a large public health burden. To date, only four epidemiological studies have investigated the relationship between air pollution and incident diabetes (Andersen et al. 2012;Coogan et al. 2012;Krämer et al. 2010;Puett et al. 2011). Two reported positive associations between incident diabetes and PM 2.5 and traffic-related pollutants such as nitrogen oxides (NO x ) and nitrogen dioxide (NO 2 ) (Coogan et al. 2012;Krämer et al. 2010). However, another study reported no association between incident diabetes and PM 2.5 (Puett et al. 2011), and the fourth reported only a small increase in diabetes associated with NO 2 (Andersen et al. 2012). Of the four cohort studies, one was based on a general population (Andersen et al. 2012), two included women only (Coogan et al. 2012;Krämer et al. 2010), and one comprised female nurses and male health professionals (Puett et al. 2011).
Given that few epidemiological studies have examined associations between air pollution and diabetes, and because little is known about the association in the general population, we conducted a population-based cohort study of long-term exposure to PM 2.5 and incident diabetes in Ontario, Canada.

Materials and Methods
Study design and population. The study was designed as a follow-up of respondents from Ontario, Canada, to the 1996/1997 cycle of National Population Health Survey of all Canadians (Statistics Canada 2010a) and the 2000/2001, 2002, 2003, and 2005 cycles of the Canadian Community Health Surveys among Canadians ≥ 12 years of age (Statistics Canada 2010b). These population-based surveys collected information related to health status, health care utilization, and determinants of health for the Canadian population in all provinces and territories, excluding fulltime members of the Canadian military, institutional residents, and individuals living on Indian Reserves, Crown Lands, and certain remote regions. The same questions were used across surveys. The response rates of the surveys in Ontario varied from 73.4 to 92.8%, depending on the year (Statistics Canada 2010a, 2010b). The surveys have been used in health research, such as estimating the burden of cardiovascular diseases in Canada (Manuel et al. 2003).
The study population for the present analysis comprised all respondents who, at the time of survey, resided in Ontario, were ≥ 35 years of age, were registered with Ontario's provincial health insurance plan, provided informed consent to share and link their responses to provincial health administrative data, and were free of diabetes. In addition, we restricted our analyses to Canadian-born individuals, because immigrants to Canada have unknown prior exposures and tend to have better health and health behaviors (McDonald and Kennedy 2004), but are more likely to live in areas with higher ground-level concentrations of PM 2.5 than those born in Canada (Villeneuve et al. 2011). A total of 62,012 participants were included in the study and they were followed-up from the time of survey until 31 December 2010.
The study was approved by the Research Ethics Board of Sunnybrook Health Sciences Centre in Toronto, Ontario, Canada.
Ascertainment of diabetes and other comorbidities. We used the Ontario Diabetes Database (Toronto, ON, CA), a validated registry of diabetics in Ontario, to identify cohort members with and without diabetes (Hux et al. 2002;Lipscombe and Hux 2007). This database was developed using hospital discharge abstracts from the Canadian Institute for Health Information and physician service claims from the Ontario Health Insurance Plan database (Toronto, ON, CA). The hospital discharge database captures all hospital admissions among Ontario residents, whereas the claims database includes claims from approximately 98% of Ontario physicians because of the universal nature of the provincial medicare system (Chan 2000;Lipscombe and Hux 2007 (Hux et al. 2002). Gestational diabetes was excluded. This algorithm has been validated using chart review and shown to have high sensitivity (86%) and specificity (97%) for identifying persons with diabetes (Hux et al. 2002). More recently, a study evaluating medical records of 22 family practitioners and their patients in southwestern Ontario confirmed the high sensitivity and specificity of this algorithm (Harris et al. 2010). Once included in the database, individuals remain in it until death or termination of Ontario health coverage. The eligibility of cohort members for health insurance and their vital status through the followup period were assessed using data from the Registered Persons Database (Toronto, ON, CA), a registry of all Ontario residents who have a health insurance number.
We linked the 62,012 cohort members to the diabetes database using encrypted health insurance numbers. We defined incident diabetes as cases diagnosed between the time of entry into the cohort and the end of followup (31 December 2010). Prevalent cases of diabetes diagnosed before the baseline survey were excluded from the analyses.
Assessment of long-term exposure to PM 2.5. To assess long-term exposure to ambient air pollution, we used satellite-based estimates of surface concentrations of PM 2.5 (van Donkelaar et al. 2010). The satellite-based concentrations were derived from aerosol optical depth data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Multiangle Imaging Spectroradiometer (MISR) instruments onboard the National Aeronautics and Space Administration (NASA)'s Terra satellite (van Donkelaar et al. 2010). Using satellite data collected from 1 January 2001 through 31 December 2006, we derived long-term average concentrations of PM 2.5 at a resolution of approximately 10 × 10 km and covered all North America below 70°N, which includes all of Ontario. These satellite-based long-term average concentrations of PM 2.5 have been shown to correlate well with ground measurements at fixed-site stations across North America (Pearson correlation coefficient r = 0.77, n = 1057) (van Donkelaar et al. 2010). These satellite-based estimates have been applied previously to estimate associations between longterm exposure to air pollution and mortality (Crouse et al. 2012;Villeneuve et al. 2011) and to estimate the global burden of illnesses due to air pollution (Brauer et al. 2012). The remote sensing methodology used to estimate ground-level ambient concentrations of PM 2.5 have been described in detail elsewhere (van Donkelaar et al. 2010). We assigned exposure by linking the exposure surface concentrations of PM 2.5 to subjects' residences at cohort entry using six-character postal codes. Six-character postal codes in urban areas represent the centroid of the blocks in which the cohort members lived.
Potential confounding variables. From the self-reported health surveys, we extracted information on marital status, race/ethnicity (white, black, Asian, Arab, Latin American, other), education attainment (less than high school, high school, some postsecondary, and university), and household income adequacy (lowest income, lower middle income, middle income, upper middle income, and upper income). Household income adequacy is an index used by Statistics Canada that accounts for total household income and household size (Statistics Canada 2010a, 2010b). Because 98% of the cohort was classified as white, we dichotomized race/ethnicity as white or nonwhite.
In addition, we extracted baseline survey data for height and weight to calculate body mass index (BMI; kilograms per meter squared). We also obtained information on smoking status (never smoker, current smoker, former smoker), alcohol consumption (more than once a month, less than once a month, former drinker, never drank), daily consumption of fruits and vegetables (< 5 servings/ day, ≥ 5 servings/day), physical activity (≥ 3.0, 1.5-2.9, < 1.5 kcal/kg/day of energy expenditure for leisure activities), and urban/rural residence. Urban areas are those continuously built-up areas having a population ≥ 1,000 and a population density ≥ 400/km 2 based on current census population counts (Statistics Canada 2010a, 2010b).
Using the Canadian Census data for 1996(Statistics Canada 2012a, 2012b, 2012c, we created contextual variables at the census tract level for the proportion of population ≥ 15 years of age with less than high school education, unemployment rate, and mean household income. Because census tracts are not defined for rural areas, we derived these three variables for rural residents according to their census subdivision, which is the next higher geographic unit, usually volume 121 | number 7 | July 2013 • Environmental Health Perspectives representing a municipality or equivalent. We assigned the values of contextual variables derived from the 1996/1997 census for individuals who entered the cohort in 1996; from the 2001 census for individuals who entered in 2000/2001, 2002, or 2003; and from the 2006 census for individuals who entered in 2005. Census tracts are small and relatively homogeneous geographic units that usually comprise a population of 2,500-8,000.
To control for regional-scale spatial patterns in the incidence of diabetes that could be explained by factors other than pollution, we created a dichotomous indicator variable classifying Ontario into southern and northern regions based on the 14 Ontario Local Health Integrated Networks. The Local Health Integrated Networks are responsible for planning, integrating, and funding various local health care services in Ontario.
Statistical analysis. We used a stratified Cox proportional hazards model with strata defined as single-year age groups, cycle of survey, and region (south/north). We excluded 1,936 individuals with missing data for BMI, and we created for all other covariates a separate category for missing values, leaving 60,076 cohort members in all analyses.
The outcome was the incident diagnosis date of diabetes as indicated in the diabetes database. Follow-up time was measured in days, calculated from the date of interview until the date of incident diabetes, or death, or were no longer eligible for provincial health insurance, or until the end of follow-up (31 December 2010).
We modeled associations between PM 2.5 and incident diabetes, adjusting for sex, marital status, education, household income adequacy, race/ethnicity, BMI (a linear term and a quadratic term), physical activity, smoking, drinking, diet, urban residency, hypertension at baseline, area-level unemployment, education, and mean household income at baseline. In separate analyses we also controlled for other comorbidities including congestive heart failure, acute myocardial infarction, COPD, and asthma.
We routinely tested for deviations from the proportional hazards assumption by assessing whether the cross-product of each variable with the natural logarithm of the time variable was statistically significant (α = 0.05). We also verified the assumption of linearity for PM 2.5 and all other continuous variables by using natural cubic spline functions with two and three degrees of freedom. We examined plots of the concentration-response curves and used the Akaike Information Criterion (AIC) to assess the relative goodness of fit for these models. A difference of > 4 AIC points is considered to prefer one model over another (Leffondre et al. 2002). Because there was no evidence of departure from linearity for the relation of PM 2.5 and diabetes, we report adjusted hazard ratios (HR) and 95% CIs for every 10-µg/m 3 increase of PM 2.5 (referred to as HR 10 ). We chose an increment of 10 µg/m 3 to facilitate comparisons of our findings with other studies (Coogan et al. 2012;Puett et al. 2011).
In addition, we investigated potential effect modification by age, sex, BMI, education, race/ ethnicity, household income adequacy, physical activity, smoking, and comorbidities by assessing whether the interaction term that was the cross-product of each variable with PM 2.5 value was statistically significant. Sensitivity analyses. Our main analysis focused on exposure assigned according to residential postal codes at cohort entry. To assess the impact of residential mobility on the effect estimates, we performed three sensitivity analyses by restricting the followup period to first 2 or 5 years, respectively; restricting to 55,708 participants who had lived at their baseline address for at least 5 years before baseline; and modeling timeweighted exposure since cohort entry until the event, with weights for each participant defined by the time spent at each address. For the latter analysis, we obtained residence for each subject and each year from the Registered Persons Database (Toronto, ON, CA) for the period 1996-2010, and derived annual estimates of PM 2.5 exposure by assigning the 6-year mean concentration of PM 2.5 at each annual postal code.
We also performed additional sensitivity analyses by restricting to participants who had used health services 1-2 years before the baseline because of a concern that the frequency of health care utilization might influence the likelihood of detecting diabetes; excluding participants with missing information on diet because two of the five surveys (1996 and 2002; 27% of the study population) did not include dietrelated questions; excluding participants from the 2002 cycle of health survey because this survey had the lowest response rate (73.4%); and restricting the analysis to southern Ontario, where 83% of the cohort lived.
We next investigated whether the association might be influenced by spatial dependence among study subjects. We fitted the Cox model with a frailty (random effect) term for Ontario Local Health Integration Networks to allow for the possibility that the effect estimate for diabetes may vary from network to network in the estimation of the main effect and its variance. A gamma distribution for the frailties was assumed, with an exchangeable correlation structure within network. We compared models with and without a frailty term using the AIC. We repeated this analysis by using a frailty term for grids from PM 2.5 exposure surface (10 × 10 km) as a random effect.
Furthermore, we examined whether associations changed over time by additionally adjusting the Cox model for annual mean concentration of PM 2.5 in Ontario (as a linear term) for the period 1996-2010. We derived the long-term trend of PM 2.5 using data from all air quality monitoring stations across Ontario that operated for at least half of the study period [see Supplemental Material, Figure S1 (http:// dx.doi.org/10.1289/ehp.1205958)]. We also tested for inter actions between time periods (1996-2000 and 2007-2010) and PM 2.5 .

Results
The cohort comprised 484,644 person-years of observations. Mean (± SD) follow-up was 8 ± 3.2 years. The mean age of the cohort at time of entry was 54.9 years (Table 1). Fortyfive percent of the cohort members were men, 63% were married, 24% were current smokers, and 54% were either overweight or obese (BMI ≥ 25 kg/m 2 ). In addition, 29% of the cohort had hypertension, 10% had COPD, and 3% had congestive heart failure at baseline. Average unemployment among the census tracts was 7%, and mean household income was approximately Can$62,300.
Associations between diabetes and PM 2.5. For every 10-µg/m 3 increase of PM 2.5 , the hazard ratio for incident diabetes (HR 10 ) was 1.08 (95% CI: 0.99, 1.17) adjusting for sex and stratifying on age, survey year, and region (Table 2). Adjustment for education, household income, and BMI strengthened the association between PM 2.5 and diabetes (HR 10 _ ∼ 1.12). Additionally controlling for physical activity, smoking, and other individual-level covariates did not appreciably change the HR. In the fully adjusted model that included all individual-level covariates, contextual covariates, and comorbid conditions, the HR 10 was 1.11 (95% CI: 1.02, 1.21). Modeling PM 2.5 using natural splines did not improve model fit according to AIC relative to the model that assumed linearity for PM 2.5 (AIC = 57441 for the model with a linear term of PM 2.5 ; AIC = 57442 for the model with natural spline), suggesting a potential log-linear relationship between PM 2.5 and incident diabetes [see Supplemental Material, Figure S2 (http://dx.doi.org/10.1289/ehp.1205958)].

Sensitivity analyses.
Restricting the analy sis to the first 2 and 5 years of followup slightly increased the effect estimates for diabetes (Table 4). Using time-weighted exposure as an alterative exposure metric did not result in appreciable difference in the estimate of association. In addition, the risk estimate was insensitive to exclusion of participants who had recently moved to their baseline addresses before the study, those who had not used health services within the previous 2 years, those who enrolled in the 2002 survey, or those who lived in northern Ontario. Additionally, we excluded participants with missing information on diet, yielding an HR 10 of 1.10 (95% CI: 0.96, 1.23). Furthermore, adding a frailty term to allow for random effects according to Ontario Local Health Integration Networks or grids from PM 2.5 exposure surface had little impact on the HR or its variance. Last, there was no strong evidence of variation in the estimated effects of PM 2.5 over time, and the test for interactions between periods and PM 2.5 was not statistically significant (p-interactions = 0.95 and 0.35).

Discussion
In this population-based cohort study of 62,012 adults in Ontario, we found that long-term exposure to PM 2.5 was associated with an increased risk of incident diabetes after controlling for various individual and neighborhood covariates. The estimate of association was insensitive to various sensitivity analyses. Additionally, we did not find strong evidence for effect modification by selected comorbidities and sociodemographic covariates.
This study provides evidence for the association between incident diabetes and PM 2.5 . Few studies have investigated the relationship between incident diabetes and air pollution. In a study of 1,775 women in the Ruhr district of Germany, Krämer et al. (2010) reported an adjusted HR of 1.27 (95% CI: 1.09, 1.48) for incident diabetes for every interquartile-range (IQR) increase of PM 2.5 (IQR = 0.4/10 5 m using absorbance-based measurement) and 1.42 (95% CI: 1.16, 1.73) per IQR of NO 2 (IQR = 8 ppb). A second study of 4,204 African-American women in Los Angeles, California, reported an HR of 1.25 (95% CI: 1.07, 1.46) per IQR of NO x (12.4 ppb) and 1.63 (95% CI: 0.78, 3.44) for a 10-µg/m 3 increase in PM 2.5 (Coogan et al. 2012).
In contrast, there was little evidence of an association between incident diabetes and particulate pollutants [PM 2.5 , PM 2.5-10 , PM 10 (PM with diameter 2.5-10 µm and ≤ 10 µm)] in a cohort of female nurses (n = 74,412) or a cohort of male health professionals (n = 15,048) living in metropolitan areas of the northeastern and midwestern United States, though residential proximity to roadways was associated with diabetes among the nurses (0-49 m vs. ≥ 200 m: HR = 1.14; 95% CI: 1.03, 1.27) (Puett et al. 2011). In Table 4. Sensitivity analyses a for the associations of incident diabetes with every 10-µg/m 3 increase of PM 2.5 .  a study of 51,818 participants in Denmark (Andersen et al. 2012), an IQR increase in NO 2 (2.6 ppb) was not associated with all cases of diabetes (HR = 1.00; 95% CI: 0.97, 1.04), but was associated with cases who were identified using a more strict case definition (HR = 1.04; 95% CI: 1.00, 1.08). Finally, there was a 4% increase in the adjusted odds (95% CI: 0%, 8%) of diabetes prevalence with each parts per billion increase of NO 2 among women in a cross-sectional study of respiratory clinic patients in Hamilton and Toronto, Ontario, but there was no association among men (Brook et al. 2008).

Sensitivity analysis
Although associations with diabetes in previous studies were relatively consistent for NO 2 , associations with PM 2.5 were not. The inconsistency may be attributable to chance, differences in population characteristics, misclassification of exposure (resulting from different exposure methods across studies), or inherent differences in the toxicological properties of the pollutants. In these prior studies, positive associations were reported mostly among women (Coogan et al. 2012;Krämer et al. 2010;Puett et al. 2011), and we also found a stronger association among women than men. This stronger association could be attributable to chance, but it may also reflect smaller errors in exposure assessment for women, because women tend to spend more time in and around home than men (Brook et al. 2008). Physiological and lifestyle differences may also contribute to the difference in effect estimates between men and women (Brook et al. 2008). Associations also varied by age, with stronger effects of PM 2.5 among those above and below 50-65 years of age at baseline. However, our data are inadequate to investigate this pattern further.
Although increasing BMI is a strong risk factor for diabetes (Colditz et al. 1995), our results do not support the idea that being overweight or obese would enhance susceptibility to the effects of air pollution on diabetes. We did not find strong evidence suggesting that other comorbidities materially altered the association between PM 2.5 and diabetes, though power to detect differences was limited. Because inflammation is a key feature of common chronic diseases such as COPD (Rana et al. 2004) and an important pathophysiological response to PM 2.5 exposure, further investigation is needed on whether the chronic state of inflammation in patients with these conditions may heighten their susceptibility to diabetes as a consequence of PM 2.5 exposure.
To our knowledge, this is the largest study of incident diabetes in a population-based cohort to date. According to the Canadian Census in 2001 (Statistics Canada 2012b), this cohort is representative of the Canadian-born population ≥ 35 years of age in Ontario (mean age = 53.3, men = 48%, married = 65%, and white = 99%). We obtained extensive individual information on known risk factors, and the diagnosis of diabetes was based on a validated registry with very high sensitivity and specificity (Hux et al. 2002). Finally, the use of satellite-based long-term average estimates of PM 2.5 ensures virtually complete spatial coverage of PM 2.5 among all cohort members. The satellite-based estimates have been shown to correlate well with ground-based measurements (van Donkelaar et al. 2010). The ambient level of PM 2.5 in Ontario (annual mean in 2000, 11.2 µg/m 3 ) was much lower than average exposures in previous studies conducted in Los Angeles, California (annual mean PM 2.5 in 2000, 20.7 µg/m 3 ) (Coogan et al. 2012) and in the Ruhr district of Germany (annual mean PM 2.5 in 2002, 22.4 µg/m 3 ) (Krämer et al. 2010).
This study has several limitations that should be considered. First, we could not differentiate between type 1 and type 2 diabetes. However, given that type 2 diabetes accounts for > 90% of all diabetes cases globally and that all cohort members were ≥ 35 years of age at entry (mean, ~ 55 years), the vast majority of incident diabetes in this cohort is expected to be type 2 diabetes (American Diabetes Association 2012).
Second, we could not identify undiagnosed cases of diabetes in the cohort. Although incomplete diagnosis of cohort members is a potential limitation of this study, effect estimates were virtually unchanged when we restricted the analysis to subjects who had used health care services during the 1-2 years before baseline as a proxy indicator of health care utilization, which may be related to the diagnosis of diabetes. Because of universal health care in Ontario, it is expected that we may have underestimated the true effects because this measurement error was likely independent of the exposure.
Third, the spatial pattern in exposure used was derived for the period 2001 to 2006 only. However, the spatial gradients of ambient PM 2.5 in Ontario remained stable during the follow-up period (1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010) and that variability in the concentrations of PM 2.5 is primarily spatial in nature and not temporal [see Supplemental Material, p. 4 (http://dx.doi. org/10.1289/ehp.1205958)]. Studies conducted in diverse locations across the United States have also demonstrated long-term stability in the spatial patterns of PM 2.5 (Jerrett et al. 2005;Miller et al. 2007;Pope et al. 2002). We therefore expect that the spatial contrasts in PM 2.5 during 2001-2006 is a reasonable representation of longer-term spatial exposure to PM 2.5 in Ontario. The spatial resolution of PM 2.5 exposure surface (i.e., 10 × 10 km), however, meant that we were unable to estimate associations at finer spatial scale. We also did not have information on daily activity. To assess the impact of residential mobility, we performed various sensitivity analyses, which did not result in appreciable differences in risk estimates. Also, our analyses did not consider the mixture of air pollutants to which subjects may have been exposed.
Fourth, we did not have family history of diabetes or occupational exposure to dust/ fumes in the surveys, and information on potential confounding variables was obtained at baseline only. Twenty-seven percent of the cohort had missing information on diet which was modeled as a separate category. However, there was little change in the association between PM 2.5 and diabetes when we excluded members with missing information on diet.
This study suggests that long-term exposure to PM 2.5 may contribute to the development of diabetes. A plausible biological mechanism linking exposure to PM 2.5 with diabetes may be indirect effects mediated through systemic proinflammatory and oxidative responses (Brook et al. 2008). Sun et al. (2009) found that exposing mice to ambient PM 2.5 at the concentration of 72 µg/m 3 for 6 hr each day over 10 weeks exacerbated insulin resistance by enhancing systemic inflammatory response and inflammation in adipose tissue. Other plausible mechanisms include autonomic nervous system imbalance and endothelial dysfunction that may be triggered by PM 2.5 , which would in turn induce vasoconstriction and result in reduced insulin sensitivity (Brook et al. 2008;Coogan et al. 2012).

Conclusions
In summary, we investigated the association between long-term exposure to PM 2.5 and the risk of incident diabetes in a large cohort in Ontario, Canada. Results from this study support a possible relationship between PM 2.5 and diabetes.