Air Pollution and Glucose Metabolism: An Analysis in Non-Diabetic Participants of the Heinz Nixdorf Recall Study.
Exposure to particulate matter (PM) in ambient air is a major environmental health risk, accounting for an estimated 3.1 million deaths and 3.1% of global disability-adjusted life years lost worldwide in 2010 (Lim et al. 2010). Short- and long-term exposure to these inhalable particles can aggravate respiratory and cardiovascular conditions, increase hospital admissions, and result in mortality from cardiovascular and respiratory diseases as well as from lung cancer (WHO 2006). Emerging evidence suggests that exposure to PM may also increase the risk of type 2 diabetes mellitus (T2DM) in the general population (Balti et al. 2014). Although specific pathophysiological mechanisms that might contribute to associations between PM and T2DM are unclear at present, one plausible hypothesis is that inhalation of particulate matter causes oxidative stress in the lungs that can lead to systemic inflammation, inflammation of adipose tissue, and insulin resistance. Together, these intermediate outcomes may contribute to a diabetogenic metabolism and eventually lead to the onset of T2DM (Franklin et al. 2015; Rajagopalan and Brook 2012).
In an effort to clarify these pathways, several epidemiological studies have explored whether higher air pollution exposure is associated with elevated blood glucose levels, a potential sign of increased insulin resistance. Short-term (days to weeks; Peng et al. 2016), medium-term (weeks to months; Peng et al. 2016; Sade et al. 2015, 2016), and long-term ([greater than or equal to] 1 year; Cai et al. 2017; Chuang et al. 2011; Liu et al. 2016; Ward-Caviness et al. 2015; Wolf et al. 2016) exposure studies have shown positive associations between a variety of air pollution measures and blood glucose levels. Nevertheless, study findings have been inconsistent for specific pollutants [e.g., particulate matter (PM), nitrogen dioxide (N[O.sub.2]), sulfur dioxide (S[O.sub.2])], and the latency period required for a cause-effect relationship has received little attention.
Because blood glucose measures are subject to high intrapersonal variability, glycated hemoglobin A1c (HbA1c), a biomarker that reflects average blood glucose levels over the previous 6-8 wk, is a useful instrument for assessing glucose levels and potential insulin resistance. At present, the epidemiological studies that have evaluated associations between long-term outdoor air pollution exposures (PM, N[O.sub.2], S[O.sub.2]) and HbA1c level show mixed results (Chuang et al. 2011; Honda et al. 2017; Liu et al. 2016; Tamayo et al. 2014, 2016; Wolf et al. 2016). The few studies investigating medium-term exposure periods, which may be the more relevant period for HbA1c levels, have also yielded mixed results (Sade et al. 2015,2016).
Additionally, few studies have examined these associations in persons without diabetes (Brook et al. 2013; Chen et al. 2016; Honda et al. 2017; Kim and Hong 2012; Peng et al. 2016; Wolf et al. 2016), an important study group for bettering our understanding of how air pollution exposure may play a role in the early development of diabetes. We examined whether exposure to medium-term (28-, 91-d mean) air pollution [particulate matter with aerodynamic diameter [less than or equal to] 2.5 [micro]m (P[M.sub.2.5]), particulate matter with aerodynamic diameter [less than or equal to] 10 [micro]m (P[M.sub.10]), N[O.sub.2], and accumulation mode particle number (P[N.sub.AM])] is associated with blood glucose and HbA1c levels in nondiabetic participants of the German population-based prospective Heinz Nixdorf Recall (HNR) study using data from two examination times. Because blood glucose levels vary more with daily changes than HbA1c levels do, we hypothesized that 28-d mean exposure windows would be more strongly associated with blood glucose measures, whereas 91-d mean exposure windows would be more strongly associated with HbA1c levels.
This study was conducted using data from the baseline (2000-2003) and first follow-up (2006-2008) examinations of the HNR study, an ongoing prospective population-based cohort study located in three adjacent cities (Bochum, Essen, and Mulheim) within the highly urbanized German Ruhr area. The study design has been described in detail elsewhere (Schmermund et al. 2002; Stang et al. 2005). Briefly, potential participants between the ages of 45 and 75 were identified through random sampling of local residency lists and were recruited via letter or telephone. Subjects were not eligible for study inclusion if they were institutionalized, had died or moved away at time of recruitment, had an incorrect or nonexistent address, did not speak a sufficient level of German, could not be interviewed because of severe illness, or were pregnant. In all, 4,814 participants were enrolled into the HNR study between December 2000 and August 2003 (recruitment efficacy proportion: 55.8%; Stang et al. 2005), and 4,157 participants returned for a follow-up examination between 2006 and 2008. Assessment at both examinations included a self-administered questionnaire, face-to-face interviews, clinical examinations, and comprehensive laboratory tests following standard protocols. The study was approved by the institutional ethics committees of the University of Duisburg-Essen and the University Hospital of Essen and adhered to strict internal and external quality assurance protocols. All participants gave informed consent.
Air pollution. In this study, we used the validated, time-dependent, three-dimensional European Air Pollution Dispersion (EURAD) chemistry transport model (Buns et al. 2012; Hass et al. 1993; Memmesheimer et al. 2004) to estimate daily mass concentrations of P[M.sub.2.5], P[M.sub.10], N[O.sub.2], and P[N.sub.AM] (particle number with aerodynamic diameter between 0.1-1.0 [micro]m; Nonnemacher et al. 2014). The EURAD model is a multilayer, multigrid model system that simulates the transport, chemical transformation, and deposition of tropospheric constituents (Buns et al. 2012). It employs four sequential nesting grid sizes (125 km, 25 km, 5 km, 1 km) for assigning exposure to Europe, central Europe, North Rhine-Westphalia in Germany, and the Ruhr area (Duisburg-Mulheim-Essen-Bochum), respectively (Buns et al. 2012; Memmesheimer et al. 2004). Specific information concerning the input data for the EURAD model have been published previously (Memmesheimer et al. 2004). The EURAD model produces hourly estimates for each grid square for a set of chemical compounds and a set of volatile organic compounds. Additionally, an assimilation process for P[M.sub.10] and N[O.sub.2] concentrations integrating measurement information from routine state-operated monitoring was conducted (measurement data on P[M.sub.2.5] and P[N.sub.AM] for this study region and period were only available from individual measurement campaigns and therefore were not used for the assimilation process). P[N.sub.AM] estimates from the EURAD model have been validated against measurements obtained between January 2011 and December 2014 by the Institute of Energy and Environmental Technology (IUTA) at its measuring station in Mulheim-Styrum using a TSI 3926 scanning mobility particle sizer spectrometer (size range: 0.014-0.750 [micro]m; TSI Inc.; for more details, see Birmili et al. 2016). For the accumulation mode (0.1-1.0 [micro]m), Pearson correlation coefficients between the model estimates and measured averages were calculated by season (spring: 1 March-31 May; summer: 1 June-31 August; fall: 1 September-30 November; winter: 1 December-29 February) and ranged between 0.51 and 0.61, with the highest correlation occurring during winter and fall seasons (H. Jakobs, unpublished data, 2016).
The HNR study area covers approximately 600 [km.sup.2] in the Ruhr area of Germany. Each participant in the HNR was assigned daily mean P[M.sub.2.5], P[M.sub.10], N[O.sub.2], and P[N.sub.AM] concentrations from the 1 [km.sup.2] grid cell corresponding to his/her given residential address (ArcView, version 9.2, ESRI; Hennig et al. 2016; Nonnemacher et al. 2014). Residential mean exposure values were calculated for each participant for the 28- and 91-d before each examination. Additional means for 1-, 2-, 3-, 7-, 14-, 45-, 60-, 75-, 105-, 120-, and 182-d exposure windows were calculated for use in exploring the temporal shape of the association of interest. Interquartile range (IQR) values were calculated for each exposure using all observations across both examinations.
Meteorological variables. Meteorological data for humidity and temperature were calculated as averages for the same 28- and 91-d time periods before examination using the Mesoscale Meteorological Model component of the EURAD model (Memmesheimer et al. 2004). In the modeling process, meteorological exposures were estimated for each 1 [km.sup.2] grid using centralized data from the German weather service (Deutscher Wetterdienst) center in Essen, Germany. Means for 1-, 2-, 3-, 7-, 14-, 45-, 60-, 75-, 105-, 120-, and 182-d meteorological windows were also calculated.
Noise measures. Long-term road noise was modeled for the year 2007 at facade points according to European Union Directive 2002/49/EC (European Commission 2002) and was calculated as the weighted 24-h mean ([L.sub.den]). Residential exposures were then assigned to participants using the maximum noise value at the most exposed facade point for the building of residence or, if building-specific information was lacking, the maximum value in a 10-m buffer around each participant's address. Noise was included as a categorical variable with the following groups: <45, 45-50, 50-55,55-60,60-65,65-70,70-75, and [greater than or equal to] 75 dB (A).
Traffic indicators. Proximity (in meters) of a participant's residence to a major road, as defined by a traffic count of >26,000 vehicles/d, was calculated using official digitized maps with a precision of at least 0.5 m. The reference line was the median strip between oncoming traffic lanes.
Blood Glucose and HbA1c Measures
Blood glucose levels (mg/dL) were measured at each examination with participants being advised to fast before the visit. Information on time since last meal (hours) was also collected to facilitate classification of each test as random or fasting. Serum samples were centrifuged immediately after being obtained and were analyzed enzymatically by the hexokinase method (Burrin and Price, 1985). HbA1c levels (percentage points, p.p.) were also measured at each examination using an automated nephelometer (BN-II; Dade Behring, Inc.). All analyses were performed in the central laboratory of the University Hospital of Essen following a standard procedure.
Definition of Covariates
Individual socioeconomic status (SES) was defined as years of education in four categories ([less than or equal to] 10, 11-13, 14-17, and [greater than or equal to] 18 years) according to the International Standard Classification of Education (UNESCO 1997). Neighborhood unemployment rates (as percentages) were obtained from local census authorities for each administrative neighborhood (median size: 11,263 inhabitants) for the years of the baseline examination (2000-2003; Dragano et al. 2009). Smoking status was defined as current, former (>1 year since quitting), or never smoker. Cumulative smoking exposure was assessed for former and current smokers using pack-years and accounting for periods of nonsmoking. Exposure to environmental tobacco smoke (ETS; yes/no) was defined as regular passive exposure to smoke at home, work, or other location. Physical activity (yes/no) was assessed as participation in regular sporting activities at least once a week for a minimum of 30 min. Alcohol consumption was obtained through a dietary questionnaire (0, 1-3, 4-6, >6 drinks per week). Anthropometric measurements (height, weight) were measured at examinations according to standard protocols, and body mass index (BMI; kg/[m.sup.2]) was calculated. Nutrition status was assessed using a dietary pattern index created by incorporating consumption frequency of 13 food items and diet quality classifications used in previous studies (Winkler and Doring 1995, 1998). Possible scores ranged from 0 to 26, with 26 representing an ideal diet, and were categorized into quartiles (<10,11-12,13-15, >15). Diabetes mellitus status (yes/no) was classified as positive if the participant reported a physician diagnosis or was taking an antihyperglycemic drug. No differentiation was made by diabetes type (1 or 2). For medication use, participants were asked to bring all medication taken in the last seven days to each examination, where confirmed use was assigned using the WHO Anatomical Therapeutic Chemical classification system (ATC; Collaborating Centre for Drug Statistics Methology, 2014). Statin use (yes/no) was thus assigned. Seasonality at the time of blood draw was modeled using Fourier series terms cos (2p X doy =365.25) and sin(2p X doy=365.25), where doy represents day of year (Peng et al. 2016). Updated information on all baseline characteristics was obtained at the follow-up study visit, with the exception of neighborhood unemployment rate and cumulative smoking exposure.
A total of 8,971 observations were collected over the two examinations of the HNR study. Observations were excluded from the analysis if the participant had a diagnosis of diabetes or was taking antihyperglycemic medication at the time of the examination ([n.sub.obs] = 973). Participant observations from the baseline examination were not excluded if participants first reported a diabetes diagnosis or antihyperglycemic medication at follow-up. Additionally, 890 observations were excluded owing to missing variable information at baseline ([n.sub.obs] = 560) or at the follow-up examination ([n.sub.obs] = 330). Of these, air pollution exposure ([n.sub.obs] = 152), HbA1c measure ([n.sub.obs] = 141), information on statin use ([n.sub.obs] = 314), alcohol consumption ([n.sub.obs] = 195), and nutrition index ([n.sub.obs] = 139) were most commonly missing. After exclusions, the main analyses for blood glucose and HbA1c included 7,108 total observations from 4,176 participants.
We evaluated the association between air pollution exposures (P[M.sub.2.5], P[M.sub.10], N[O.sub.2], P[N.sub.AM]) and blood sugar measures (blood glucose, HbA1c) using a linear mixed-effects regression model with a random participant intercept to account for the correlation between individuals at baseline and at the follow-up examination. Separate 28- and 91-d mean exposure models were run for each air pollutant. Model estimates were calculated per IQR increase of exposure and are presented as mean exposure point estimates with 95% confidence intervals (CIs). Because significant nonlinearity was present for age, temperature, and humidity, these variables were modeled using restricted cubic splines with four knots. The normality of the model residuals was checked in a pre-analysis.
Three models (Models 1-3) of increasing covariate adjustment were conducted for each air pollutant exposure. Updated covariate information was used for all follow-up observations, with the exception of neighborhood unemployment and cumulative smoking (information only available at baseline). In Model 1, we adjusted for temperature, humidity, and examination. In Model 2, we additionally adjusted for age, season, nutrition index, smoking status, and BMI. For blood glucose analyses, time since last meal was also included. Covariates included for adjustment in Model 2 were identified using a directed acyclic graph (DAG; see Figure S1; Textor et al. 2011). To assess whether residual confounding was present, we further adjusted in Model 3 (the main model) for sex, education level, and additional lifestyle factors, including physical activity, alcohol consumption, statin use, ETS exposure, cumulative smoking exposure, and neighborhood unemployment. Covariates included in Model 3 were selected based on inclusion in prior studies (Peng et al. 2016; Wolf et al. 2016).
To explore whether associations between air pollution exposure and glucose metabolism varied according to exposure time window, we ran the main model analysis using a range of exposure windows before examination (1-, 2-, 3-, 7-, 14-, 28-, 45-, 60-, 75-, 91-, 105-, 120-, and 182-d means) for each pollutant. We chose the 28- and 91-d exposures a priori for the main analyses and treated the additional time points as an exploratory analysis.
For P[N.sub.AM], we repeated analyses using multipollutant models adjusted for P[M.sub.2.5], P[M.sub.10], or N[O.sub.2] (as continuous variables) to determine whether associations with this exposure metric were independent of associations with more traditional air pollution markers. Spearman correlation coefficients were calculated between all air pollutant exposures (28- and 91-d means).
Analyses were also conducted to investigate whether different associations were apparent in observations designated as prediabetic ([n.sub.obs] = 4,114) based on a fasting blood glucose level between 100 mg/dL and 125 mg/dL or an HbA1c level between 5.7% and 6.5% (American Diabetes Association 2016) compared with observations having no indication of prediabetes. We also looked at whether associations differed by working status (<15 hours/week vs. [greater than or equal to] 15 hours/week), because persons who worked less had a greater likelihood of spending their time at home and thus theoretically more accurate exposure assignments.
To evaluate the robustness of our results to the inclusion of traffic noise exposure, we added noise ([L.sub.den]) and proximity to major roads to the main model individually as well as jointly. We also repeated analyses after excluding participants with a high blood glucose test (fasting blood glucose [greater than or equal to] 126 mg/dL or non-fasting blood glucose [greater than or equal to] 200 mg/dL) as well as participants who had a diagnosis of diabetes or who were using antihyperglycemic medication at the time of examination. The blood glucose analyses were also run after excluding participants who were not fasting at the time of blood draw ([n.sub.obs] = 2,139; time since last meal <8 hours) to limit any effects due to incomplete fasting or to differences in fasting between examinations of the same individual.
We evaluated potential effect modification by sex (male, female), smoking status (current/former/never), examination (baseline/follow-up), age (<65/ [greater than or equal to] 65 y), season (spring/summer/fall/winter), and regular physical activity (yes/no) through addition of a multiplicative interaction term between air pollution exposure (28- and 91-d means; continuous) and the covariate of interest into the main analysis model. Interactions were considered significant if the corresponding likelihood ratio test yielded a p-value <0.05.
All statistical analyses were conducted in R version 3.3.0 (R Core Team).
A total of 7,108 blood glucose and HbA1c measurements were collected from 4,176 participants in the HNR study without diabetes and with complete covariate data (Table 1). Among fasting participants, mean blood glucose levels [[+ or -] standard deviation (SD)] were 107.97 [+ or -] 17.4 mg/dL and 107.6 [+ or -] 15.8 mg/dL at baseline and at the follow-up examination, respectively. Glucose levels were similar for observations from nonfasting participants (104.6 [+ or -] 17.0 mg/dL at baseline and 103.9 [+ or -] 14.5 mg/dL at follow-up), possibly explained by the fact that although these participants did not fast the complete 8 h, the mean fasting time was still approximately 4.5 h. Mean HbA1c levels ([+ or -] SD) at baseline and at follow-up were 5.40 [+ or -] 0.64 p.p. and 5.52 [+ or -] 0.52 p.p., respectively. Compared with information reported by participants at the baseline examination ([n.sub.obs] = 3,857), participants at the follow-up examination ([n.sub.obs] = 3,251) reported higher consumption of alcoholic drinks per week and lower exposure to ETS and were more likely to be former smokers and to report statin use. For all air pollutants, mean exposure levels decreased between visits (e.g., from 17.4 to 15.9 [micro]g/[m.sup.3] for 28-d P[M.sub.2.5]; Table 2). P[N.sub.AM] exposure levels were moderately correlated with other air pollutants (p = 0.53, 0.51, and 0.33 for 28-d P[M.sub.2.5], P[M.sub.10], and N[O.sub.2], respectively; see Table S1).
Participants excluded because of missing data (covariate, exposure, and/or outcome) at the baseline examination were more likely to have attended their examination in the winter, to report exposure to ETS, to have [less than or equal to] 10 years of education, and to report being physically inactive than participants with complete data from the baseline examination (see Table S2). Participants excluded because of missing data at the follow-up examination were more likely to have attended their examination in the spring and to be former smokers than participants with complete data at follow-up.
Main Analysis: Air Pollution and Glucose Metabolism
In Model 1, P[M.sub.2.5] (28- and 91-d), P[M.sub.10] (28-d), and P[N.sub.AM] (28- and 91-d) exposures were positively associated with blood glucose (Table 3). Similarly, P[M.sub.2.5], P[M.sub.10], and P[N.sub.AM] exposures (28 d and 91 d) were positively associated with HbA1c. Further adjustment for additional covariates identified using a DAG (Model 2) slightly reduced the point estimates for blood glucose but yielded similar estimates for HbA1c.
In Model 3, IQR increases in 28-d mean exposure to P[M.sub.2.5], P[M.sub.10], and P[N.sub.AM] were positively associated with blood glucose levels (Table 3). Significant positive associations were also seen for 91-d mean exposure to P[M.sub.2.5]. For N[O.sub.2], neither 28- nor 91-d exposure levels were clearly associated with blood glucose levels. After evaluation of the association between a range of short-and medium-term exposures and blood glucose, we observed the strongest positive associations with 28-, 45-, and 60-d exposure windows for PM and P[N.sub.AM] (Figure 1). Associations between N[O.sub.2] and blood glucose levels were close to the null for most exposure windows, although nonsignificant negative associations were estimated for the longest time windows (75- to 182-d).
For HbA1c, significant positive associations were observed for 28-d and 91-d mean exposure to P[M.sub.2.5], P[M.sub.10], and P[N.sub.AM] using Model 3 (Table 3), with stronger point estimates being observed for 91-d than 28-d mean exposures. For N[O.sub.2], neither 28- nor 91-d exposure levels were associated with HbA1c. Using the range of exposure windows, we observed the strongest positive associations with slightly longer exposure windows for P[M.sub.2.5], P[M.sub.10], and P[N.sub.AM] (Figure 2) than for glucose. Associations between N[O.sub.2] and HbA1c levels were null for all exposure windows.
In the multipollutant models, the association between P[N.sub.AM] and blood glucose was attenuated when P[M.sub.2.5] was included in the model but remained after adjustment for P[M.sub.10] and N[O.sub.2] (see Figure S2). All associations between P[N.sub.AM] and HbA1c were robust to adjustment for P[M.sub.2.5], P[M.sub.10], and N[O.sub.2] exposure (see Figure S2).
When considering nondiabetic and prediabetic participants, associations for both blood glucose and HbA1c were attenuated among prediabetic participants ([n.sub.obs] = 3,881), but they remained in the same direction as those in the main analysis (see Figures S3 and S4). Estimates did not differ greatly by working status ([n.sub.obs] = 4,624 working part-time or less), when additional observations were excluded based on a broader definition of diabetes ([n.sub.obs] = 423 additionally excluded), or when analyses were limited to only fasting participants ([n.sub.obs] = 4,969). Model estimates were also very robust to single and simultaneous addition of proximity of the participant's residence to a major road and chronic noise exposure ([n.sub.obs] = 7,048).
For blood glucose, we observed evidence of effect modification by examination for P[N.sub.AM] (28-, 91-d) with associations that were positive for the baseline examination but close to the null at the follow-up (p-values<0.01; see Figures S5 and S6). We also observed significant effect modification by season for P[M.SUP.2.5] (91-d; p-value = 0.01), N[O.sub.2] (91-d; p-value = 0.02), and P[N.sub.AM] (28-d; p-value <0.01), with no clear seasonal pattern across pollutants. In analyses of HbA1c, we observed evidence of effect modification by season for all air pollutants, with stronger associations between air pollutants and HbA1c during the spring season (see Figures S7 and S8). Effect modification by examination was also present for 91-d P[M.sup.2.5], P[M.sub.10], and P[N.sub.AM], with stronger associations apparent at the baseline examination [e.g., 0.09 (95% CI: 0.06, 0.13) vs. 0.06 (95% CI: 0.03,0.09)] for P[M.sub.2.5]). Additionally, stronger associations were seen among former smokers than among current smokers and nonsmokers for P[M.sub.10] (28- and 91-d p-values: 0.04 and 0.02, respectively) and P[N.sub.AM] (91-d; p-value = 0.02). Although not statistically significant, the same pattern was also observed for P[M.sup.2.5] and N[O.sub.2].
In the present study, we showed that higher levels of residential medium-term exposure to P[M.sub.2.5], P[M.sub.10], and P[N.sub.AM] were positively associated with blood glucose levels and HbA1c among nondiabetic adults. These associations were robust to extensive adjustment for lifestyle factors, personal characteristics, traffic noise, and meteorological covariates. Importantly, we were able to observe a consistent temporal pattern of associations between air pollution exposure and glucose metabolism measures across a range of short- and medium-term exposure windows. The slightly different patterns observed for blood glucose and HbA1c are biologically plausible because more recent exposures were most strongly associated with blood glucose, whereas longer exposure windows were most strongly associated with HbA1c levels. To our knowledge, this is the first study to also examine the effects of particle number concentration (accumulation mode), a submicron particle metric of growing interest in the health community, on blood glucose levels and HbA1c using a spatiotemporal exposure model.
As described in previous reviews (Liu et al. 2013; Rajagopalan and Brook 2012; Rao et al. 2015), several pathophysiologic pathways may explain the associations between air pollution exposures and glucose metabolism. Inhalation of air pollution has been shown to induce low-grade oxidative stress and inflammation in the lungs (Rajagopalan and Brook 2012), and these effects do not remain localized in the lung but can lead to systemic oxidative stress and inflammation, including in the adipose tissue (Rao et al. 2015; Sun et al. 2005, 2009). Systemic and adipose tissue inflammation have been observed to lead to impaired insulin signaling in the pathways that moderate glucose metabolism (Haberzettl et al. 2016; Rao et al. 2015).
Evidence from several epidemiologic studies supports this connection between air pollution and glucose metabolism given that researchers have observed associations between short- and medium-term air pollution exposure and increased inflammation (Rajagopalan and Brook 2012) as well as worsening insulin resistance (Brook et al. 2013; Chen et al. 2016; Kelishadi et al. 2009; Kim and Hong 2012). Similarly, ongoing research suggests that medium-term exposures to air pollution during pregnancy may be related to abnormal glucose regulation (Fleisch et al. 2014; Lu et al. 2017; Robledo et al. 2015) and incidence of gestational diabetes mellitus among pregnant women (Hu et al. 2015). The results from the present study provide additional evidence that short- and medium-term air pollution exposures may represent biologically relevant exposure windows for influencing glucose metabolism measures. These associations with glucose metabolism may help explain emerging epidemiologic data that there is an association between long-term air pollution exposure and incidence of diabetes mellitus (Balti et al. 2014).
Comparison with Prior Studies
The observed association between medium-term exposure to PM and P[N.sub.AM] and blood glucose concentration has been previously shown for PM measures (Kim and Hong 2012; Peng et al. 2016). In a recent study among nondiabetic participants in the Normative Aging Study using a hybrid spatiotemporal prediction model, Peng et al. (2016) observed a significant positive association between 28-d P[M.sup.2.5] exposure and fasting blood glucose (0:89 mg/dL per 3.12 [micro]g/[m.sup.3] increase). Chen et al. (2016) also observed associations between centrally measured short-term P[M.sup.2.5] exposures and fasting blood glucose among a cohort of Mexican Americans with high risk of diabetes (e.g., 0.70 mg/dL per 5 [micro]g/[m.sup.3] 7-d mean P[M.sup.2.5]). In contrast to several previous studies (Chen et al. 2016; Chuang et al. 2011; Honda et al. 2017; Kim and Hong 2012; Sade et al. 2015; Wolf et al. 2016), we did not observe an association between N[O.sub.2] exposure and blood glucose level, for which we have no clear explanation.
We specifically investigated nondiabetic participants, an important group for evaluating how air pollution exposure may influence early disease mechanisms that lead to a diabetogenic metabolism and eventually to T2DM in the general population. Three other studies including persons of mixed diabetes status have reported positive associations between air pollution exposures and blood glucose levels (Sade et al. 2015, 2016; Ward-Caviness et al. 2015). However, effects of diet, physical activity, and medication on blood glucose levels in participants with diabetes may overshadow relatively small effects of air pollution, suggesting that these effects might be easier to detect among metabolically healthy persons with intact glucose metabolism.
An important and novel finding of our study is that we found a positive association between medium-term PM and P[N.sub.AM] exposure and HbA1c levels in persons with no known metabolic disorder. This finding is of pathophysiological importance because it suggests that air pollution exposure may lead to a diabetogenic metabolism, manifesting itself as increased glucose levels over periods of several weeks to months. Of the two prior studies among nondiabetic participants, one reported no associations with P[M.sub.2.5], N[O.sub.2], or a variety of other exposure metrics (Wolf et al. 2016), whereas the second reported an association with N[O.sub.2] but not with P[M.sup.2.5] (Honda et al. 2017). Although all three studies were conducted in Western populations, the divergence in the results may be due to the difference in exposure windows used. With our spatiotemporal exposure model, we were able to choose biologically motivated medium-term (28- and 91-d mean) exposures before the blood draw for our analysis, whereas Wolf et al. (2016) utilized a temporally static long-term mean, and Honda et al. (2017) used 1- to 5-y moving average exposure windows. From a biological standpoint, medium-term exposures (averaged over weeks to months) would be more likely to be associated with HbA1c than exposures averaged over longer periods of time given that HbA1c levels in individual red blood cells increase with age, and the lifespan of a typical red blood cell is approximately 115 d (Franco 2012).
In recent years, experimental and epidemiologic evidence has emerged suggesting that submicron and ultrafine (UFP; aerodynamic diameter <0:1 [micro]m) particles may have significant adverse health effects owing to their large surface area-to-mass ratio, their ability to generate reactive oxygen species, and their small size, which enables them to penetrate deeply into the respiratory system (Li et al. 2016). Nevertheless, epidemiologic studies linking exposure to submicron particles or UFPs with health effects are limited because they are not routinely measured by monitoring networks, and their heterogeneous distribution in the atmosphere makes them challenging to model and measure. To our knowledge, only three studies have investigated whether UFP exposure is associated with diabetes-related outcomes. One small study, conducted among 58 Danish households of unknown diabetes status, observed that 48-h exposure to indoor UFPs was associated with HbA1c levels, but it did not observe any connection between 48-h outdoor UFP exposure (assigned from a central monitoring station) and HbA1c (Karottki et al. 2014). In a study of near-highway households, no association was observed between annual UFP exposure and prevalent diabetes (Lietal. 2017). In contrast, a recent Canadian cohort study by Bai et al. (2018) observed a positive association between annual ambient UFP exposure and risk of diabetes. In the present study, we considered the accumulation mode rather than ultrafine particles specifically and found positive associations between both short-term (28-d) and medium-term (91-d) outdoor P[N.sub.AM] exposure and two measures of glucose metabolism that merit further investigation.
Although the increases in blood glucose and HbA1c levels observed in the present study (e.g., increase in HbA1c of 0.15 p. p. per 10 [micro]g/[m.sup.3] P[M.sub.2.5]) are small, and the exposure range in the Ruhr area of Germany (e.g., 91-d mean P[M.sub.2.5] IQR: 4.0 [micro]g/[m.sup.3]) is low, air pollution exposure is ubiquitous and unavoidable for all members of the population. Even small increases in blood glucose and HbA1c levels may increase the risk of cardiometabolic disease (Zhang et al. 2010) as well as of cardiovascular events and mortality (Cohen et al. 2009; Schottker et al. 2016). Thus, even these small increases per unit of air pollution exposure may have relevant and wide-reaching impacts on health worldwide.
Strengths and Limitations
This study has several strengths. The HNR study is a well-characterized, population-based cohort with extensive covariate data, which allowed us to adjust for potential confounders as well as for the co-occurring exposure of long-term traffic noise. Additionally, the use of a spatiotemporal model allowed us to estimate a range of short- and medium-term exposures. The consistency and the temporal pattern across short- to medium-term exposure windows also support the hypothesis that a true association exists between air pollution exposure and glucose metabolism measures among nondiabetic adults.
Nevertheless, there are several limitations of this study. Because the EURAD model estimates urban background exposures on a 1 [km.sup.2] grid and not directly at participants' residences, some exposure misclassification, particularly in traffic-related air pollution, likely exists. We expect this misclassification to be nondifferential and thus, in principle, to bias our results towards the null. Additionally, it is possible that a background air pollution estimate for the 1 [km.sup.2] around a residence may not represent a good estimate of overall exposure because people typically do not spend the entirety of their time at home. Assimilated estimates of P[M.sub.2.5] and P[N.sub.AM] are un-available in this study area owing to a small number of measurements; thus, error in these exposure estimates is likely greater than for P[M.sub.10] and N[O.sub.2].
Medium-term exposure to PM and P[N.sub.AM] was positively associated with blood glucose and HbA1c levels, consistent with an adverse effect on glucose regulation in nondiabetic adults. Further studies investigating PM components and source-specific particulate matter as well as possible mediation pathways are needed to better understand the mechanisms by which air pollution exposure influences blood glucose and HbA1c levels in metabolically healthy persons.
We are indebted to the investigative group and the study personnel of the Heinz Nixdorf Recall study. We also thank the North Rhine-Westphalia State Agency for Nature, Environment and Consumer Protection for providing emission and land use data for North Rhine-Westphalia. We thank A. Buschka for her data management work and U. Quass from the Institute of Energy and Environmental Technology for his help in describing the validation study. The Heinz Nixdorf Recall study was funded by the Heinz Nixdorf Foundation, Germany [Chairman: M. Nixdorf; Past Chairman: Dr. jur. G. Schmidt (deceased)]. The HNR was also supported by grants from the German Research Council [Deutsche Forschungsgemeinschaft (DFG); ER 155/6-1, ER 155/6-2, SI 236/8-1, and SI 236/9-1] and the Kulturstiftung Essen, Germany. The AirFlamm III study was supported by the German Research Council (DFG; HO 3314/4-3). This work was completed in partial fulfillment of the requirements of the Ph.D. thesis for S.L.
American Diabetes Association. 2016. Erratum. Classification and diagnosis of diabetes. Sec. 2. In Standards of Medical Care in Diabetes-2016. Diabetes Care 2016;39(Suppl. 1):S13-S22. Diabetes Care 39(9):1653 S13-S22, PMID: 26696675, https://doi.org/10.2337/dc16-S005.
Bai L, Chen H, Hatzopoulou M, Jerrett M, Kwong JC, Burnett RT, et al. 2018. Exposure to Ambient Ultrafine Particles and Nitrogen Dioxide and Incident Hypertension and Diabetes. Epidemiology [Epub ahead of print], PMID: 29319630, https://doi.org/10.1097/EDE.0000000000000798.
Balti EV, Echouffo-Tcheugui JB, Yako YY, Kengne AP. 2014. Air pollution and risk of type 2 diabetes mellitus: A systematic review and meta-analysis. Diabetes Res Clin Pract 106(2):161-172, PMID: 25262110, https://doi.org/10.1016/j.diabres. 2014.08.010.
Birmili W, Weinhold K, Rasch F, Sonntag A, Sun J, Merkel M, et al. 2016. Long-term observations of tropospheric particle number size distributions and equivalent black carbon mass concentrations in the German Ultrafine Aerosol Network (GUAN). Earth Syst Sci Data 8(2):355-382, https://doi.org/10.5194/essd-8-355-2016.
Brook RD, Xu X, Bard RL, Dvonch JT, Morishita M, Kaciroti N, et al. 2013. Reduced metabolic insulin sensitivity following sub-acute exposures to low levels of ambient fine particulate matter air pollution. Sci Total Environ 448:66-71, PMID: 22901427, https://doi.org/10.1016/j.scitotenv.2012.07.034.
Buns C, Klemm O, Wurzler S, Hebbinghaus H, Steckelbach I, Friesel J, et al. 2012. Comparison of four years of air pollution data with a mesoscale model. Atmos. Res 118:404-417, https://doi.org/10.1016/j.atmosres.2012.07.009.
Burrin JM, Price CP. 1985. Measurement of blood glucose. Ann Clin Biochem 22(Part 4):327-342, PMID: 3898972, https://doi.org/10.1177/000456328502200401.
Cai Y, Hansell AL, Blangiardo M, Burton PR, de Hoogh K, Doiron D, et al. 2017. Long-term exposure to road traffic noise, ambient air pollution, and cardiovascular risk factors in the HUNT and Lifelines cohorts. Eur. Heart J 44:1-8, PMID: 28575405, https://doi.org/10.1093/eurheartj/ehx263.
Chen Z, Salam MT, Toledo-Corral C, Watanabe RM, Xiang AH, Buchanan TA, et al. 2016. Ambient air pollutants have adverse effects on insulin and glucose homeostasis in Mexican Americans. Diabetes Care 39(4):547-554, PMID: 26868440, https://doi.org/10.2337/dc15-1795.
Chuang KJ, Yan YH, Chiu SY, Cheng TJ. 2011. Long-term air pollution exposure and risk factors for cardiovascular diseases among the elderly in Taiwan. Occup Environ Med 68(1):64-68, PMID: 20833756, https://doi.org/10.1136/oem.2009.052704.
Cohen BE, Barrett-Connor E, Wassel CL, Kanaya AM. 2009. Association of glucose measures with total and coronary heart disease mortality: does the effect change with time? The Rancho Bernardo Study. Diabetes Res Clin Pract 86(1):67-73, PMID: 19671481, https://doi.org/10.1016/j.diabres.2009.07.006.
Dragano N, Hoffmann B, Stang A, Moebus S, Verde PE, Weyers S, et al. 2009. Subclinical coronary atherosclerosis and neighbourhood deprivation in an urban region. Eur J Epidemiol 24(1):25-35, PMID: 18931923, https://doi.org/10. 1007/s10654-008-9292-9.
European Commission. 2002. Directive 2002/49/EC of the European Parliament and of the Council of 25 June 2002 relating to the assessment and management of environmental noise. Off J Eur Commun L18918.72002:12-25. http://eur-lex.europa. eu/legal-content/EN/TXT/?uri=CELEX:32002L0049 [accessed 22 May 2017].
Fleisch AF, Gold DR, Rifas-Shiman SL, Koutrakis P, Schwartz JD, Kloog I, et al. 2014. Air pollution exposure and abnormal glucose tolerance during pregnancy: the Project Viva cohort. Environ Health Perspect 122(4):378-383, PMID: 24508979, https://doi.org/10.1289/ehp.1307065.
Franco RS. 2012. Measurement of red cell lifespan and aging. Transfus Med Hemother 39(5):302-307, PMID: 23801920, https://doi.org/10.1159/000342232.
Franklin BA, Brook R, Arden Pope C. 2015. Air pollution and cardiovascular disease. Curr Probl Cardiol 40(5):207-238, PMID: 25882781, https://doi.org/10.1016/ j.cpcardiol.2015.01.003.
Haberzettl P, O'Toole TE, Bhatnagar A, Conklin DJ. 2016. Exposure to fine particulate air pollution causes vascular insulin resistance by inducing pulmonary oxidative stress. Environ Health Perspect 124(12):1830-1839, PMID: 27128347, https://doi.org/10.1289/EHP212.
Hass H, Ebel A, Feldmann H, Jakobs HJ, Memmesheimer M. 1993. Evaluation studies with a regional chemical transport model (EURAD) using air quality data from the EMEP monitoring network. Atmos Environ A Gen Top 27(6):867-887, https://doi.org/10.1016/0960-1686(93)90007-L.
Hennig F, Sugiri D, Tzivian L, Fuks K, Moebus S, Jockel KH, et al. 2016. Comparison of land-use regression modeling with dispersion and chemistry transport modeling to assign air pollution concentrations within the Ruhr area. Atmosphere (Basel) 7(3):48, https://doi.org/10.3390/atmos7030048.
Honda T, Pun VC, Manjourides J, Suh H. 2017. Associations between long-term exposure to air pollution, glycosylated hemoglobin and diabetes. Int J Hyg Environ Health 220(7):1124-1132. PMID: 28712959, https://doi.org/10.1016/j. ijheh.2017.06.004.
Hu H, Ha S, Henderson BH, Warner TD, Roth J, Kan H, et al. 2015. Association of atmospheric particulate matter and ozone with gestational diabetes mellitus. Environ Health Perspect 123(9):853-859, PMID: 25794412, https://doi.org/10. 1289/ehp.1408456.
Karottki DG, Beko G, Clausen G, Madsen AM, Andersen ZJ, Massling A, et al. 2014. Cardiovascular and lung function in relation to outdoor and indoor exposure to fine and ultrafine particulate matter in middle-aged subjects. Environ Int 73:372-381, PMID: 25233101, https://doi.org/10.1016/j.envint.2014.08.019.
Kelishadi R, Mirghaffari N, Poursafa P, Gidding SS. 2009. Lifestyle and environmental factors associated with inflammation, oxidative stress and insulin resistance in children. Atherosclerosis 203(1):311-319, PMID: 18692848, https://doi.org/10.1016/j. atherosclerosis.2008.06.022.
Kim JH, Hong YC. 2012. GSTM1, GSTT1, and GSTP1 polymorphisms and associations between air pollutants and markers of insulin resistance in elderly Koreans. Environ Health Perspect 120(10):1378-1384, PMID: 22732554, https://doi.org/10. 1289/ehp.1104406.
Li N, Georas S, Alexis N, Fritz P, Xia T, Williams MA, et al. 2016. A work group report on ultrafine particles (American Academy of Allergy, Asthma & Immunology): Why ambient ultrafine and engineered nanoparticles should receive special attention for possible adverse health outcomes in human subjects. J Allergy Clin Immunol 138(2):386-396, PMID: 27130856, https://doi.org/10.1016/j.jaci.2016.02.023.
Li Y, Lane KJ, Corlin L, Patton AP, Durant JL, Thanikachalam M, et al. 2017. Association of long-term near-highway exposure to ultrafine particles with cardiovascular diseases, diabetes and hypertension. Int J Environ Res Public Health 14(5):E461-16, PMID: 28445425, https://doi.org/10.3390/ijerph14050461.
Lim SS, Vos T, Flaxman AD, Danaei G, Shibuya K, Adair-Rohani H, et al. 2012. A comparative risk assessment of burden of disease and injury attributable to 67 riskfactors and riskfactor clusters in 21 regions, 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet 380(9859):2224-2260, PMID: 23245609, https://doi.org/10.1016/S0140-6736(12)61766-8.
Liu C, Yang C, Zhao Y, Ma Z, Bi J, Liu Y, et al. 2016. Associations between long-term exposure to ambient particulate air pollution and type 2 diabetes prevalence, blood glucose and glycosylated hemoglobin levels in China. Environ Int 92-93:416-421, PMID: 27148900, https://doi.org/10.1016/j.envint.2016.03.028.
Liu C, Ying Z, Harkema J, Sun Q, Rajagopalan S. 2013. Epidemiological and experimental links between air pollution and type 2 diabetes. Toxicol Pathol 41(2):361-373, PMID: 23104765, https://doi.org/10.1177/0192623312464531.
Lu MC, Wang P, Cheng TJ, Yang CP, Yan YH. 2017. Association of temporal distribution of fine particulate matter with glucose homeostasis during pregnancy in women of Chiayi City, Taiwan. Environ Res 152:81-87, PMID: 27743970, https://doi.org/10. 1016/j.envres.2016.09.023.
Memmesheimer M, Friese E, Ebel A, Jakobs HJ, Feldmann H, Kessler C, et al. 2004. Long-term simulations of particulate matter in Europe on different scales using sequential nesting of a regional model. IJEP 22(1/2):108-132, https://doi.org/10. 1504/IJEP.2004.005530.
Nonnemacher M, Jakobs H, Viehmann A, Vanberg I, Kessler C, Moebus S, et al. 2014. Spatio-temporal modelling of residential exposure to particulate matter and gaseous pollutants for the Heinz Nixdorf Recall Cohort. Atmos Environ 91:15-23, https://doi.org/10.1016/j.atmosenv.2014.03.052.
Peng C, Bind M-AC, Colicino E, Kloog I, Byun H-M, Cantone L, et al. 2016. Particulate air pollution and fasting blood glucose in non-diabetic individuals: associations and epigenetic mediation in the Normative Aging Study, 2000-2011. Environ. Health Perspect 124(11):1715-1721, PMID: 27219535, https://doi.org/10.1289/EHP183.
Rajagopalan S, Brook RD. 2012. Air pollution and type 2 diabetes: Mechanistic insights. Diabetes 61(12):3037-3045, PMID: 23172950, https://doi.org/10.2337/db12-0190.
Rao X, Patel P, Puett R, Rajagopalan S. 2015. Air pollution as a riskfactor for type 2 diabetes. Toxicol Sci 143(2):231-241, PMID: 25628401, https://doi.org/10.1093/ toxsci/kfu250.
Robledo CA, Mendola P, Yeung E, Mannisto T, Sundaram R, Liu D, et al. 2015. Preconception and early pregnancy air pollution exposures and risk of gestational diabetes mellitus. Environ Res 137:316-322, PMID: 25601734, https://doi.org/10. 1016/j.envres.2014.12.020.
Sade MY, Kloog I, Liberty IF, Katra I, Novack L, Novack V. 2015. Air pollution and serum glucose levels: a population-based study. Medicine (Baltimore) 94(27): e1093, PMID: 26166095, https://doi.org/10.1097/MD.0000000000001093.
Sade MY, Kloog I, Liberty IF, Schwartz J, Novack V. 2016. The association between air pollution exposure and glucose and lipids levels. J Clin Endocrinol Metab 101(6):2460-2467, PMID: 27218271, https://doi.org/10.1210/jc.2016-1378.
Schmermund A, Mohlenkamp S, Stang A, Gronemeyer D, Seibel R, Hirche H, et al. 2002. Assessment of clinically silent atherosclerotic disease and established and novel risk factors for predicting myocardial infarction and cardiac death in healthy middle-aged subjects: Rationale and design of the Heinz Nixdorf RECALL Study. Am Heart J 144(2):212-218, PMID: 12177636, https://doi.org/10.1067/mhj.2002.123579.
Schottker B, Rathmann W, Herder C, Thorand B, Wilsgaard T, Njolstad I, et al. 2016. HbA1c levels in non-diabetic older adults - No J-shaped associations with primary cardiovascular events, cardiovascular and all-cause mortality after adjustment for confounders in a meta-analysis of individual participant data from six cohort studies. BMC Med 14:26, PMID: 26867584, https://doi.org/10.1186/s12916016-0570-1.
Stang A, Moebus S, Dragano N, Beck EM, Mohlenkamp S, Schmermund A, et al. 2005. Baseline recruitment and analyses of nonresponse of the Heinz Nixdorf Recall study: identifiability of phone numbers as the major determinant of response. Eur J Epidemiol 20(6):489-496, PMID: 16121757, https://doi.org/10. 1007/s10654-005-5529-z.
Sun Q, Wang A, Jin X, Natanzon A, Duquaine D, Brook RD, et al. 2005. Acceleration of atherosclerosis and vascular inflammation in an animal model. JAMA 294(23):30033010, PMID: 16414948, https://doi.org/10.1001/jama.294.23.3003.
Sun Q, Yue P, Deiuliis JA, Lumeng CN, Kampfrath T, Mikolaj MB, et al. 2009. Ambient air pollution exaggerates adipose inflammation and insulin resistance in a mouse model of diet-induced obesity. Circulation 119(4):538-546, PMID: 19153269, https://doi.org/10.1161/CIRCULATI0NAHA.108.799015.
Tamayo T, Rathmann W, Kramer U, Sugiri D, Grabert M, Holl RW. 2014. Is particle pollution in outdoor air associated with metabolic control in type 2 diabetes? PLoS One 9(3):e91639, PMID: 24619127, https://doi.org/10.1371/journal.pone.0091639.
Tamayo T, Rathmann W, Stahl-Pehe A, Landwehr S, Sugiri D, Kramer U, et al. 2016. No adverse effect of outdoor air pollution on HbA1c in children and young adults with type 1 diabetes. Int J Hyg Environ Health 219(4-5):349-355, PMID: 26935923, https://doi.org/10.1016/j.ijheh.2016.02.002.
Textor J, Hardt J, Knuppel S. 2011. DAGitty: A graphical tool for analyzing causal diagrams. Epidemiology 22(5):745, PMID: 21811114, https://doi.org/10.1097/EDE. 0b013e318225c2be.
UNESCO. 1997. International standard classification of education ISCED 1997. http:// www.unesco.org/education/information/nfsunesco/doc/isced_1997.htm [accessed 30 March 2017].
Ward-Caviness CK, Kraus WE, Blach C, Haynes CS, Dowdy E, Miranda ML, et al. 2015. Association of roadway proximity with fasting plasma glucose and metabolic risk factors for cardiovascular disease in a cross-sectional study of cardiac catheterization patients. Environ Health Perspect 123(10):1007-1014, PMID: 25807578, https://doi.org/10.1289/ehp.1306980.
WHO (World Health Organization). 2006. WHO Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide-Global update 2005. http:// www.who.int/phe/health_topics/outdoorair/outdoorair_aqg/en/index.html [accessed 29 September 2016].
WHO Collaborating Centre for Drug Statistics Methodology. 2014. ATC/DDD Index 2014. http://www.whocc.no/atc_ddd_index/[accessed 16 March 2017].
Winkler G, Doring A. 1995. Kurzmethoden zur Charakterisierung des Ernahrungsmusters: Einsatz und Auswertung eines Food-Frequency-Fragebogens. Ernahrungs-Umschau 42(8):289-291.
Winkler G, Doring A. 1998. Validation of a short qualitative food frequency list used in several German large scale surveys. Z Ernahrungswiss 37(3):234-241, PMID: 9800314, https://doi.org/10.1007/PL00007377.
Wolf K, Popp A, Schneider A, Breitner S, Hampel R, Rathmann W, et al. 2016. Association between long-term exposure to air pollution and biomarkers related to insulin resistance, subclinical inflammation, and adipokines. Diabetes 65(11):3314-3326, PMID: 27605624, https://doi.org/10.2337/db15-1567.
Zhang X, Gregg EW, Williamson DF, Barker LE, Thomas W, Bullard KM, et al. 2010. A1C level and future risk of diabetes: A systematic review. Diabetes Care 33(7):1665-1673, PMID: 20587727, https://doi.org/10.2337/dc09-1939.
Sarah A. Lucht, (1) Frauke Hennig, (1) Clara Matthiessen, (1) Simone Ohlwein, (1) Andrea Ichs, (2,3) Susanne Moebus, (4) Karl-Heinz Jockel, (4) Hermann Jakobs, (5) and Barbara Hoffmann (1)
(1) Environmental Epidemiology Group, Institute of Occupational, Social and Environmental Medicine, Medical Faculty, Heinrich-Heine University Dusseldorf, Dusseldorf, Germany
(2) Institute for Health Services Research and Health Economics, Centre for Health and Society, Medical Faculty, Heinrich-Heine University Dusseldorf, Dusseldorf, Germany
(3) Institute for Health Services Research and Health Economics, German Diabetes Center (DDZ), Dusseldorf, Germany
(4) Institute of Medical Informatics, Biometry and Epidemiology (IMIBE), University Hospital Essen, University of Duisburg-Essen, Essen, Germany
(5) Rhenish Institute for Environmental Research (RIU), University of Cologne, Cologne, Germany
Address correspondence to S. Lucht, Environmental Epidemiology Group, Institute for Occupational, Social and Environmental Medicine, Medical Faculty, Heinrich-Heine University Dusseldorf, Postfach 101007; 40001 Dusseldorf. Telephone: 49 (0) 211586729114. Email: sarah.lucht@ uni-duesseldorf.de
Supplemental Material is available online (https://doi.org/10.1289/EHP2561).
The authors declare that they have no actual or potential competing financial interests.
Received 21 July 2017; Revised 26 February 2018; Accepted 26 February 2018; Published 3 April 2018.
Note to readers with disabilities: EHP strives to ensure that all journal content is accessible to all readers. However, some figures and Supplemental Material published in EHP articles may not conform to 508 standards due to the complexity of the information being presented. If you need assistance accessing journal content, please contact firstname.lastname@example.org. Our staff will work with you to assess and meet your accessibility needs within 3 working days.
Caption: Figure 1. Associations between an interquartile range (IQR) increase in mean air pollution exposure and blood glucose level (mg/dL) using a range of short-and medium-term exposure windows before blood draw. Linear mixed-effects regression models with random participant intercepts were run using the main model (covariates provided in Table 3) and are shown stratified by air pollutant (N[O.sub.2], nitrogen dioxide; P[M.sub.2.5], particulate matter with aerodynamic diameter [less than or equal to] 2.5 [micro]m; P[M.sub.10], particulate matter with aerodynamic diameter [less than or equal to] 10 [micro]m; P[N.sub.AM], accumulation mode particle number concentration). Error bars represent the 95% confidence interval for each point estimate. IQR values are provided in Table 2.
Caption: Figure 2. Associations between an interquartile range (IQR) increase in mean air pollution exposure and glycated hemoglobin A1c (HbA1c; percentage points) using a range of short- and medium-term exposure windows before blood draw. Linear mixed-effects models with random participant intercepts were run using the main model (covariates provided in Table 3) and are shown stratified by air pollutant (N[O.sub.2], nitrogen dioxide; P[M.sub.2.5], particulate matter with aerodynamic diameter [less than or equal to] 2.5 [micro]m; P[M.sub.10], particulate matter with aerodynamic diameter [less than or equal to] 10 [micro]m; P[N.sub.AM], accumulation mode particle number concentration). Error bars represent the 95% confidence interval for each point estimate. IQR values are provided in Table 2.
Table 1. Demographic characteristics of the Heinz Nixdorf Recall study participant observations included in this study, stratified by examination time (baseline: 2000-2003; follow-up: 2006-2008). Baseline (n = 3,857) Mean [+ or -] SD Variable or n (%) Age (y) 59.4 [+ or -] 7.8 HbA1c (percentage points) 5.4 [+ or -] 0.6 Blood glucose (mg/dL) 106.7 [+ or -] 17.3 Time since last meal (h) 10.1 [+ or -] 4.7 Fasting Status (yes) 2,442 (63.3) BMI (kg/[m.sup.2]) 27.7 [+ or -] 4.5 Neighborhood unemployment 12 4 [+ or -] 3 4 (%) Cumulative smoking 27.9 [+ or -] 27.7 (pack-years) (a) Sex (male) 1,879 (48.7) Regular physical activity (yes) 2,149 (55.7) Statin use (yes) 378 (9.8) ETS exposure (yes) 1,378 (35.7) Season at blood draw Spring 1,032 (26.8) Summer 1,098 (28.5) Fall 917 (23.8) Winter 810 (21.0) Smoking status Never smoker 1,634 (42.4) Former smoker 1,321 (34.2) Current smoker 902 (23.4) Education level [less than or equal to] 10 years 397 (10.3) 11-13 years 2,154 (55.8) 14-17 years 883 (22.9) [greater than or equal to] 18 years 423 (11.0) Nutrition index <10 650 (16.9) 11-12 1,281 (33.2) 13-15 880 (22.8) >15 1,046 (27.1) Alcoholic drinks per week 0 1,877 (48.7) 1-3 693 (18.0) 4-6 378 (9.8) >6 909 (23.6) Noise density (dB) <45 633 (16.5) 45-50 873 (22.8) 50-55 783 (20.4) 55-60 490 (12.8) 60-65 431 (11.3) 65-70 424 (11.1) 70-75 164 (4.3) >75 31 (0.8) Missing 28 (0.7) Follow-up (n = 3,251) Mean [+ or -] SD Variable or n (%) Age (y) 64.0 [+ or -] 7.7 HbA1c (percentage points) 5.5 [+ or -] 0.5 Blood glucose (mg/dL) 106.8 [+ or -] 15.6 Time since last meal (h) 12.0 [+ or -] 4.5 Fasting Status (yes) 2,527 (77.8) BMI (kg/[m.sup.2]) 27.9 [+ or -] 4.7 Neighborhood unemployment -- (%) Cumulative smoking -- (pack-years) (a) Sex (male) 1,551 (47.7) Regular physical activity (yes) 1,911 (58.8) Statin use (yes) 551 (16.9) ETS exposure (yes) 807 (24.8) Season at blood draw Spring 897 (27.6) Summer 731 (22.5) Fall 797 (24.5) Winter 826 (25.4) Smoking status Never smoker 1,409 (43.3) Former smoker 1,266 (38.9) Current smoker 576 (17.7) Education level [less than or equal to] 10 years 304 (9.4) 11-13 years 1,830 (56.3) 14-17 years 730 (22.5) [greater than or equal to] 18 years 387 (11.9) Nutrition index <10 462 (14.2) 11-12 1,073 (33.0) 13-15 737 (22.7) >15 979 (30.9) Alcoholic drinks per week 0 1,170 (36.0) 1-3 670 (20.6) 4-6 248 (7.6) >6 1,163 (35.8) Noise density (dB) <45 547 (17.0) 45-50 767 (23.8) 50-55 641 (19.9) 55-60 419 (13.0) 60-65 345 (10.7) 65-70 349 (10.8) 70-75 133 (4.1) >75 18 (0.6) Missing 32 (1.0) Note: --, no information was collected at that particular examination point; BMI, body mass index; ETS, environmental tobacco smoke; HbA1c, glycated hemoglobin A1c; SD, standard deviation. (a) Among current or former smokers only. Table 2. Summary statistics for residential 28- and 91-day mean exposure levels at the baseline (2000-2003) and follow-up (2006-2008) examinations from the EURAD model. Baseline (n = 3,857) Exposure Mean [+ or -] SD 28-Day Exposures [PM.sub.2.5] ([micro]g/[m.sup.3]) 17.4 [+ or -] 3.8 [PM.sub.10] ([micro]g/[m.sup.3]) 20.9 [+ or -] 4.9 N[O.sub.2] ([micro]g/[m.sup.3]) 39.9 [+ or -] 12.0 [PN.sub.AM] (n/mL) 3,710.2 [+ or -] 1,908.3 Temperature ([degrees]C) 10 5 [+ or -] 6 6 Humidity (%) 6 6 [+ or -] 2 2 91-Day Exposures [PM.sub.2.5] ([micro]g/[m.sup.3]) 17.6 [+ or -] 2.7 [PM.sub.10] ([micro]g/[m.sup.3]) 21.1 [+ or -] 4.0 N[O.sub.2] ([micro]g/[m.sup.3]) 40.6 [+ or -] 10.7 [PN.sub.AM] (n/mL) 3,744.9 [+ or -] 1,343.6 Temperature ([degrees]C) 10 0 [+ or -] 5 9 Humidity (%) 6.4 [+ or -] 1.9 Baseline (n = 3,857) Exposure Range 28-Day Exposures [PM.sub.2.5] ([micro]g/[m.sup.3]) 10.3-31.2 [PM.sub.10] ([micro]g/[m.sup.3]) 11.1-44.3 N[O.sub.2] ([micro]g/[m.sup.3]) 18.9-104.2 [PN.sub.AM] (n/mL) 888.4-15,540.0 Temperature ([degrees]C) -1.9-24.3 Humidity (%) 2.9-11.1 91-Day Exposures [PM.sub.2.5] ([micro]g/[m.sup.3]) 12.3-27.3 [PM.sub.10] ([micro]g/[m.sup.3]) 13.2-36.3 N[O.sub.2] ([micro]g/[m.sup.3]) 20.7-76.7 [PN.sub.AM] (n/mL) 1,662.0-9,775.0 Temperature ([degrees]C) 0.1-20.1 Humidity (%) 3.7-10.0 Follow-up (n = 3,251) Exposure Mean [+ or -] SD 28-Day Exposures [PM.sub.2.5] ([micro]g/[m.sup.3]) 15.9 [+ or -] 5.0 [PM.sub.10] ([micro]g/[m.sup.3]) 19.3 [+ or -] 6.0 N[O.sub.2] ([micro]g/[m.sup.3]) 37.8 [+ or -] 10.9 [PN.sub.AM] (n/mL) 3,011.2 [+ or -] 1,530.1 Temperature ([degrees]C) 10.7 [+ or -] 5.8 Humidity (%) 6.5 [+ or -] 1.9 91-Day Exposures [PM.sub.2.5] ([micro]g/[m.sup.3]) 15 9 [+ or -] 3 0 [PM.sub.10] ([micro]g/[m.sup.3]) 19 5 [+ or -] 4 3 N[O.sub.2] ([micro]g/[m.sup.3]) 38.0 [+ or -] 9.0 [PN.sub.AM] (n/mL) 3,018.9 [+ or -] 842.8 Temperature ([degrees]C) 10 7 [+ or -] 5 2 Humidity (%) 6.5 [+ or -] 1.7 Follow-up (n = 3,251) Exposure Range IQR 28-Day Exposures [PM.sub.2.5] ([micro]g/[m.sup.3]) 7.5-33.7 5.7 [PM.sub.10] ([micro]g/[m.sup.3]) 8.7-40.6 7.4 N[O.sub.2] ([micro]g/[m.sup.3]) 15.6-80.5 15.7 [PN.sub.AM] (n/mL) 691.5-9,961.0 2,142.3 Temperature ([degrees]C) 0.8-24.7 11.5 Humidity (%) 3.6-10.9 3.8 91-Day Exposures [PM.sub.2.5] ([micro]g/[m.sup.3]) 9.6-25.2 4.0 [PM.sub.10] ([micro]g/[m.sup.3]) 10.6-34.9 5.5 N[O.sub.2] ([micro]g/[m.sup.3]) 18.1-65.9 15.0 [PN.sub.AM] (n/mL) 1,457.0-6,781.0 1,352.7 Temperature ([degrees]C) 3.5-20.5 10.6 Humidity (%) 4.3-10.0 3.5 Note: EURAD, European Air Pollution Dispersion; IQR, interquartile range; N[O.sub.2], nitrogen dioxide; PM2.5, particulate matter with aerodynamic diameter [less than or equal to] 2.5 [micro]m; [PM.sub.10], particulate matter with aerodynamic diameter [less than or equal to] 10 [micro]m; [PN.sub.AM], accumulation mode particle number; SD, standard deviation. Table 3. Estimated association (95% confidence interval) between interquartile range increase in 28- and 91-day mean air pollution and HbA1c and blood glucose (nobs = 7,108). Blood glucose (mg/dL) (a) 28-Day estimate 91-Day estimate Exposure (95% CI) (95% CI) [PM.sub.2.5] Model 1 (b) 1.04 (0.50, 1.58) 0.84 (0.15, 1.52) Model 2 (c) 0.97 (0.44, 1.50) 0.78 (0.01, 1.55) Model 3 (Main) (d) 0.91 (0.38, 1.44) 0.81 (0.05, 1.58) [PM.sub.10] Model 1 (b) 0.77 (0.21, 1.34) 0.37 (-0.24, 0.98) Model 2 (c) 0.70 (0.15, 1.26) 0.21 (-0:44, 0.86) Model 3 (Main) (d) 0.59 (0.04, 1.14) 0.10 (-0.54, 0.75) N[O.sub.2] Model 1 (b) 0.45 (-0.28, 1.18) 0.19 (-0.80, 1.17) Model 2 (c) 0.48 (-0.31, 1.13) -0.23 (-1.29,0.83) Model 3 (Main) (d) 0.27 (-0:45, 0.99) -0.41 (-1.46, 0.65) P[N.sub.AM] Model 1 (b) 0.74 (0.17, 1.31) 0.86 (0.29, 1.44) Model 2 (c) 0.75 (0.17, 1.32) 0.75 (0.19, 1.32) Model 3 (Main) (d) 0.64 (0.07, 1.21) 0.67 (0.10, 1.24) HbA1c (percentage points) 28-Day estimate 91-Day estimate Exposure (95% CI) (95% CI) [PM.sub.2.5] Model 1 (b) 0.03 (0.02, 0.05) 0.04 (0.01,0.06) Model 2 (c) 0.03 (0.01, 0.05) 0.07 (0.04, 0.10) Model 3 (Main) (d) 0.03 (0.01, 0.05) 0.07 (0.04, 0.10) [PM.sub.10] Model 1 (b) 0.04 (0.02, 0.06) 0.03 (0.00, 0.05) Model 2 (c) 0.04 (0.02, 0.06) 0.04 (0.02, 0.06) Model 3 (Main) (d) 0.04 (0.02, 0.06) 0.04 (0.02, 0.06) N[O.sub.2] Model 1 (b) 0.02 (-0.01, 0.04) -0:02 (-0:05, 0.02) Model 2 (c) 0.01 (-0.01, 0.04) 0.00 (-0:03, 0.04) Model 3 (Main) (d) 0.01 (-0:02, 0.04) 0.00 (-0.04, 0.04) P[N.sub.AM] Model 1 (b) 0.04 (0.02, 0.06) 0.09 (0.07,0.11) Model 2 (c) 0.03 (0.01, 0.06) 0.09 (0.07,0.11) Model 3 (Main) (d) 0.03 (0.01, 0.05) 0.09 (0.07,0.11) Note: All models were linear mixed-effects regression models with random participant intercepts. Interquartile range values are provided in Table 2. CI, confidence interval; HbA1c, glycated hemoglobin A1c; N[O.sub.2], nitrogen dioxide; PM2.5, particulate matter with aerodynamic diameter <2.5 [micro]m; [PM.sub.10], particulate matter with aerodynamic diameter <10 [micro]m; [PN.sub.AM], accumulation mode particle number concentration. (a) All blood glucose models were additionally adjusted for time since last meal (hours). (b) Adjusted for humidity, temperature, and examination. (c) Model 1 + age, smoking status, nutrition index, season, and BMI. (d) Model 2 + physical activity, alcohol consumption, statin use, and exposure to environmental tobacco smoke at each visit; and pack-years of smoking and neighborhood unemployment at baseline.
|Printer friendly Cite/link Email Feedback|
|Author:||Lucht, Sarah A.; Hennig, Frauke; Matthiessen, Clara; Ohlwein, Simone; Ichs, Andrea; Moebus, Susanne;|
|Publication:||Environmental Health Perspectives|
|Date:||Apr 1, 2018|
|Previous Article:||Bias Amplification in Epidemiologic Analysis of Exposure to Mixtures.|
|Next Article:||Prenatal Organophosphate Pesticide Exposure and Traits Related to Autism Spectrum Disorders in a Population Living in Proximity to Agriculture.|