Ecological correlation between diabetes hospitalizations and fine particulate matter in Italian provinces

Fine particulate matter air pollution (PM2.5) is associated with an increased risk for cardiovascular events and this evidence is well supported in the literature [1]. Several studies demonstrate that during periods of high ambient air pollution, diabetic patients have increased susceptibility for vascular reactivity and doubled rates of hospital admission for heart disease [2, 3]. It has been suggested that some biological mechanisms related to air pollution exposure and cardiovascular diseases may also be involved in the onset of type 2 diabetes [4]. A recent experimental study on mice exposed to PM2.5 for 10 months (equivalent to a human exposure period of ~ 40 years) shows that this exposure was sufficient to cause insulin resistance and increased the type 2 diabetes prevalence and susceptibility to cardiovascular diseases [5].

There are now at least few epidemiologic studies showing some association between air pollution and diabetes, including a recent meta-analysis [6]. An ecological study has found that diabetes prevalence among adults in USA was greater in areas with higher PM2.5 concentrations [7]. This association was strong and the increased risk of diabetes was present even among areas that are below the EPA legal limits for PM2.5 [7]. In a German perspective cohort study, women exposed to higher levels of traffic-related air pollution (NO₂ and PM) developed type 2 diabetes at a higher rate than controls with a risk increase over 16 years of follow-up by 15 % per 1 interquartile range of traffic-related PM exposure [8]. Another cohort study conducted in the Denmark suggested that the risk for diabetes was weakly and positively associated with increasing mean levels of traffic-related air pollution at the residence [9]. In contrast, associations between diabetes and exposure to particulate matter in the year before diagnosis, were not found in the Nurses’ Health Study [10]. However, an association between incident diabetes and residential proximity to a road was statistically significant [10].

In Italy there are well established evidences of the effect of increased air pollution levels on mortality for cardiovascular and respiratory diseases [11, 12]. Exposures to PM10 were positively associated with cardiovascular hospital admissions in 9 Italian cities [13]. A large cohort study on the most populated city of Italy (Rome) showed that long term exposure to NO2 and PM2.5 was associated with ischemic heart diseases, cardiovascular diseases and lung cancer [14]. Despite diabetes prevalence and mortality have been linked to its common risk factors on the Italian territory [15], there are not yet evidences if diabetes is also correlated with levels of air pollution.

Here we present a cross-sectional study assessing the ecological association between diabetes hospitalizations in Italy and ambient PM2.5 concentration at the province level. The relationship is adjusted for census level covariates for common risk factors, socio economic factors and for appropriateness of hospital admission with diabetes. All analyses were carried out from datasets retrieved from official and freely available databases.

Type 2 diabetes mellitus has been shown to be linked with some risk factors like gender [22], ageing [23], overweight or obesity [24], physical inactivity [25], ethnicity [26], genetic factors [27] and others [28]. More recently, several studies have also suggested a potential role of particulate matter in air pollution, in the development of diabetes. Biological plausible mechanisms include, among others, increased systemic oxydative stress and inflammatory responses to air pollutants (reviewed elsewhere [29, 30]). In animal models, exposure to PM2.5 resulted in adipose inflammatory changes (pro-inflammatory to anti-inflammatory macrophage ratio and insulin signalling abnormalities) which could lead to insulin resistance [5, 31]. Recent research identified cellular mechanisms that are implied in insulin resistance development caused by joint exposure to high fat diet and high levels of PM2.5 [32]. Human studies have given mixed evidences of the association between particulate matter and systemic inflammation and/or insulin resistance possibly due to differences in study populations, exposure levels and specific markers examined [2, 33–35].

In our analysis we used PCA to reduce the (correlated) variables describing the various risk and socioeconomic factors into a synthetic variable (PC1) for subsequent modelling with linear regression. Obesity, overweight, physical inactivity, cigarette smoking were all positively correlated with PC1 while educational level and GDP were inversely related. As expected, diabetes hospital discharge was positively related with PC1 and inversely related to hospitalization appropriateness in both male and female models. Regarding particulate matter pollution, our ecological analysis showed that diabetes hospitalizations (of 45 years or older age) for men increase by 3.5 % and women increase by approximately 4 % each unit increase of PM2.5.

Epidemiological studies examining the association between air pollution from traffic emissions and incidence or prevalence of diabetes are scarce and produced non consistent results. In the Toronto area, [4] used hospital discharges to estimate prevalence and found only for women an increased prevalence by 17 % across the interquartile range of NO2. Contrary to those findings, a cross-sectional ecological study conducted in the Netherlands [36] found no associations between diabetes prevalence (assessed through self reported previous diagnosis of the diabetic condition by a doctor) and exposure to traffic related pollution (though there were some indications for a relation with traffic in a 250 m buffer from the residential address of participants). A more recent study in the US found a 1 % diabetes prevalence (estimated from self reporting of diabetes diagnosis by participants) with an increase by 10 μg/m³ of PM2.5, after adjusting for many risk factors [7]. By comparison, our estimated increase per 10 μg/m³ PM2.5 is much higher (35 % for men and 40 % for women). However, the PM2.5 concentration range examined in [7] is lower, having a range between 2.5–17.7 μg/m³ (with median = 11 μg/m³) while in our study we had a larger range of PM2.5 (11–32 μg/m³) and a higher median value (18.68 μg/m³) that may explain the difference in diabetes increase. Furthermore, diabetes prevalence in [7] is based on self reported previous diagnosis of the diabetic condition by a doctor collected through a population survey, while we based our study on hospital discharge data. Prevalence estimates based on self reporting from individuals such as in surveys may miss those individuals that are yet not aware of the disease. Also hospital discharge data with diabetes as principal diagnosis possibly underestimate the true population prevalence due to the comorbidity of diabetes with cardiovascular conditions that are often the primary diagnosis of hospitalization. It is also true that many general practitioners nowadays screen for insulin resistance those (especially if of older age) individuals carrying some risk factors, therefore increasing their likelihood of hospitalization.

To date we are aware of 5 cohort studies that examined the association between air pollution and incidence of type 2 diabetes. [8] found a significant increase of the hazard for incident diabetes of 15–42 % per interquartile range of traffic-related exposure (PM10, NO2) in a women cohort from the highly industrialized Ruhr district in Western Germany. More recently, in Ontario (Canada), [37] found an 11 % increased diabetes incidence rate by a 10 μg/m3 increased PM2.5 level in a cohort of more than 60,000 participants. Contrary to those findings, [38] did not found significant association between incident diabetes and PM2.5, while they found an increase by 25 % of incident diabetes with interquartile range increase of NO2 (12.4 ppb) in a cohort of black women living in Los Angeles. No associations of incident diabetes with PM2.5 or PM10 were found in the analysis of 2 prospective cohorts in the US [10], although an association with distance to road (used as proxy for exposure to traffic related pollution) was reported. Finally another recent cohort study conducted in Denmark found a borderline statistical association between confirmed cases of diabetes and NO2 levels [9]. Notably, the risk was higher between non smokers and physically active people.

In Italy, hospital use by diabetic patients is substantially greater than that by the non-diabetic population [39]. The persistence of sub-optimal care can be responsible for an increased risk of complications, which in turn determines an increased demand for hospital care [40]. Therefore, excess hospitalization and longer hospital stays for patients with diabetic complications can be seen as an indicator of the quality of primary diabetes care [41]. In our ecological analysis the indicator variable of hospital length of stay for diabetes with complications (routinely used as appropriateness index by the Ministry of Health to monitor primary health care) was highly significant. Therefore, we advocate that an index of appropriate hospitalization should always be considered in ecological analyses dealing with hospital discharges. Previous analysis of hospital admissions for diabetes in Italy in different years has shown a decreasing trend from 2005 (13.7 per1000 inhabitants) to 2009 (12.0 per 1000 inhabitants). While this effect may be connected with increased appropriate hospitalization of diabetic patients and better management of primary cares, still Italian regions exhibit large differences in hospital discharge [19].

Ecological analyses are affected by well known limitations and our study does not allow us to conclude about causal relationship between higher fine particulate matter levels and increased hospital discharge diabetes in the Italian territory. It is important to remember that linear relationship at ecological aggregated level could not hold at the individual level and/or might be non linear. Our analysis was based on official and validated data and covers roughly 2 thirds of the Italian population and all major urban areas. Additionally, we used a cross sectional approach by adopting a 3 year averaging of all variables to increase the statistical strength, as suggested in a similar study [7]. We used covariates available from official surveys designed for providing meaningful data at the regional level of aggregation therefore assuming that the average characteristics are representative of the population in that area. Furthermore, we assumed that average exposure of individuals can be estimated from PM2.5 levels from all ground monitoring stations on a given territory. In sensitivity analysis we entered in models different metrics for annual PM2.5 measured levels (mean, median, 95 percentile and maximum) and levels estimated from a dynamic model with 4 × 4 km of resolution and still resulted significant terms. Like others [7], we also note that it is possible that the hospitals attracting more diabetic patients are placed in most polluted areas. However, when we entered in the models a term correcting for geographical macro-region (North East, North West, South) aiming in capturing residual unexplained geographical variation in hospital discharges, still the PM2.5 term remained significant. Finally, we could have omitted variables that could explain part of the observed relationship between PM2.5 and diabetes hospital discharges such as NO2 concentrations that are reported to have the same ecological relationship with diabetes prevalence.

In conclusion, we have found a significant ecological relationship between sex and age standardised hospital discharge with diabetes as principle diagnosis and PM2.5 concentrations in Italian provinces, once that covariates have been accounted for. The relationship was robust to different means of estimating PM2.5 exposure. A large portion of the variance of diabetes hospitalizations was linked to differences of primary care appropriateness between Italian regions and this variable should routinely be included in ecological analyses of hospitalizations.