Type 2 diabetes mellitus (T2DM) is a growing public health problem in Bulgaria. While individual and lifestyle determinants have been researched; till date there has been no study on environmental risks such as road traffic, noise, and air pollution. As a first step toward designing a large-scale population-based survey, we aimed at exploring the overall associations of prevalent T2DM with exposures to road traffic, noise, and air pollution. A total of 513 residents of Plovdiv city, Bulgaria were recruited. Individual data on self-reported doctor-diagnosed T2DM and confounding factors were linked to objective and self-rated exposure indicators. Logistic and log-link Poisson regressions were conducted. In the fully adjusted logistic models, T2DM was positively associated with exposures to L den 71-80 dB (odds ratio (OR) = 4.49, 95% confidence interval (CI): 1.38, 14.68), fine particulate matter (PM) 2.5 25.0-66.8 μg/m 3 (OR = 1.32, 95% CI: 0.28, 6.24), benzo alpha pyrene 6.0-14.02 ng/m 3 (OR = 1.76, 95% CI: 0.52, 5.98) and high road traffic (OR = 1.40, 95% CI: 0.48, 4.07). L den remained a significant risk factor in the: Poisson regression model. Other covariates with consistently high multivariate effects were age, gender, body mass index, family history of T2DM, subjective sleep disturbance, and especially bedroom location. We concluded that residential noise exposure might be associated with elevated risk of prevalent T2DM. The inferences made by this research and the lessons learned from its limitations could guide the designing of a longitudinal epidemiological survey in Bulgaria.
Keywords: Air pollution, noise, particulate matter (PM), polycyclic aromatic hydrocarbons, road traffic, type 2 diabetes mellitus (T2DM)
|How to cite this article:|
Dzhambov AM, Dimitrova DD. Exposures to road traffic, noise, and air pollution as risk factors for type 2 diabetes: A feasibility study in Bulgaria. Noise Health 2016;18:133-42
|How to cite this URL:|
Dzhambov AM, Dimitrova DD. Exposures to road traffic, noise, and air pollution as risk factors for type 2 diabetes: A feasibility study in Bulgaria. Noise Health [serial online] 2016 [cited 2019 May 21];18:133-42. Available from: http://www.noiseandhealth.org/text.asp?2016/18/82/133/181996
Diabetes mellitus is a growing public health problem in Bulgaria. According to the International Diabetes Federation, Bulgaria was one the European countries with relatively high prevalence of diabetes - 7.63% of those aged 20-79 years in 2013.  In 2012, 9.6% of the population had diabetes, with 26.1% of those still being undiagnosed.  That proportion had increased since 2006 when the prevalence was 7.9%.  Nearly 90% of all cases in 2012 could be attributed to type 2 diabetes mellitus (T2DM) (International Classification of Diseases (ICD) code - 10: E11). Additionally, 3.7% had prediabetes.  Given the low socioeconomic standard in the country and the dysfunctional healthcare system with limited resources hampering the access of patients to quality health services,  various preventive measures to reduce the burden of diabetes are gaining more and more prominence. However, lifestyle and genetic factors cannot fully explain the pandemic of cardiometabolic diseases worldwide.  While the high prevalence of insufficient physical activity, obesity, and metabolic syndrome in Bulgaria makes healthy lifestyle a logical target for preventive interventions,  health-enhancing behaviors will not tackle the risks associated with environmental determinants such as noise and air pollution.
These exposures are believed to have additive but distinct pathways of increasing the risk of T2DM. With respect to air pollution, sufficient data has been accumulated and has been synthetized in recently published meta-analysis.  Several biological mechanisms underlying the effect on T2DM have been discussed: Endothelial dysfunction, decreased peripheral glucose uptake, hepatic insulin resistance, systemic inflammation, and oxidative stress.
Conversely, only a few studies have explored the risk associated with noise pollution, their results are conflicting and their methodologies are largely incompatible.  Nevertheless, significantly higher risk of T2DM has been reported among those exposed to high residential noise.  The detrimental effects of noise are probably mediated through stress response, disruption of normal sleep patterns, and chronic sleep deprivation.  To this moment, there has been no previous study in Bulgaria on the environmental risk factors for T2DM. However, given the limited resources for environmental research in the country, and the fact that medical institutions give priority to clinical investigations, prior to designing a large scale longitudinal survey we believe such evidence is warranted, exploring its feasibility, supplementing basic effect size estimates that could be used to determine the required sample size and give us a feel of the anticipated practical importance of our findings; furthermore, possible methodological pitfalls could be highlighted and addressed at this preliminary stage in order to avoid additional costs and sources of bias in a population-based survey. Therefore, we aimed at exploring the overall associations of prevalent T2DM with objective and self-rated indicators of exposure to road traffic, noise, and air pollution in Plovdiv, Bulgaria.
| Methods|| |
Design and study population
This was a cross-sectional small scale observational survey. It was carried out in the period July-November, 2014 in the city of Plovdiv, Bulgaria. It is the second largest city in the country with a population of about 341 thousand people and a territory of 101.98 km 2 ,  chosen for the present study as it is characterized by noise and air pollution exceeding the acceptable standards. , Likewise, because more than 80% of diabetes-related mortality occurs in low- and middle-income countries,  this makes Plovdiv, Bulgaria an appropriate setting for our research.
A questionnaire comprising 59 closed and open-ended questions were administered in order to explore the associations between various neighborhood characteristics and several health endpoints. Only some of these questions were relevant to the objectives of the present study. In order to ensure sufficient sample size, two data collection procedures were employed - a nonprobability snowball sampling and a field sampling. Current residence in Plovdiv, Bulgaria and being ≥18 years of age were mandatory criteria for inclusion in the survey. People reporting hearing loss or uncorrected hearing impairment were excluded.
Snowball sampling started with a small group (n = 57) of acquaintances, colleagues, relatives, and students who were contacted either in person or via e-mail; they filled the questionnaire and were then asked to recruit members of their social network who became first wave participants. Each participant had to recruit up to 10 others to ensure continuity of the procedure. Because snowball sampling is a nonprobability approach and generates biased estimates due to the initial selection of seeds, no more than 50% of the predetermined number of participants were to be collected via this procedure. Simultaneously with it, the lead researcher carried out random field interviews among residents of Plovdiv, Bulgaria. As a first step, 41 neighborhood blocks were selected on an interactive city map. The researcher then visited those neighborhoods starting from the one closest to the center of the city. Participants were approached in local green spaces, at the street, in front of their homes or in local shopping centers, and they were invited to participate in the survey. They were informed that their residential satisfaction was surveyed. The interviews took about 10 min. If ,for some reason, the interview was discontinued before completion of the questionnaire, the latter was still coded in the dataset if sufficient data on exposure or outcome variables were provided, or if other variables were available that could be used to impute the missing values (e.g., self-reported instead of objective exposure indicator). No financial incentives were offered.
The necessary sample size was calculated for a multivariate logistic regression model based on the most reliable exposure variable with lowest rate of missing values anticipated, i.e., self-reported traffic intensity (STI), which is considered a proxy for both road traffic air pollution and noise.  The sample size necessary to detect an OR of 2.11 for prevalent T2DM among people exposed to "extreme traffic" in comparison to those exposed to "none/very rare traffic" was 558. 
A total of 249 questionnaires were distributed via snowball sampling, and 213 were returned completed (response rate = 85.5%); 1,906 people were approached during the field sampling and of those only 368 agreed to participate (response rate = 19.3%). Main reasons for nonparticipation were lack of time and/or fear of revealing personal information.
Because of its noninvasive and observational nature and adult population the study did not undergo approval by an Ethics Committee.  All participants were assured of anonymity and participation was voluntary. Answering the questionnaire was taken to constitute informed consent.
Each participant answered a question about their address in Plovdiv, Bulgaria. For those reporting an exact address, it was geocoded using BatchGeocode (BatchGeo LLC, 2014)  and Google Earth Pro TM v. 7.1 (Google Inc., Mountain View, CA).  When, due to privacy and confidentiality, concerns participants had provided only the name of the street they lived at, the addresses were matched to the center point of the street, according to Agay-Shay et al.  Addresses not matched automatically or mismatched were coded manually. When possible, the building of residence was identified based on the street name and description of the building's characteristics (number of apartment floor, type of housing, presence of a backyard garden, etc.). Relevant exposure data were assigned after overlaying the air pollution and traffic noise maps onto the map of geocoded addresses.
Road traffic is the dominant transportation source of noise in Plovdiv, Bulgaria with railway, aircraft, and industrial noise levels being less influential. L den - defined as "average" noise levels during daytime, evening, and nighttime, applying a 5-dB penalty to noise in the evening and a 10-dB penalty to noise during the night - was chosen for a noise indicator. Combined L den data from all traffic sources were elicited from the official strategic noise maps (10 m × 10 m grid, 4 m height) of Plovdiv, Bulgaria (http://www.webcitation.org/6biU8dDj0) created in 2009 in compliance with the Environmental Noise Directive 2002/49/EC.  Noise levels were assessed at the coordinates of the residential address unless the participant had indicated the orientation of the apartment. L den was modelled with LimA v. 5 (Brüel and Kjær, Nærum, Denmark). Further description of the noise mapping process is available elsewhere. 
Noise annoyance was derived from an 11-point scale: "To what extent are you disturbed, annoyed, or irritated by noise when you are at home?" ("0, not at all" to "10, very much"). Lifetime exposure to loud noise was measured using two questions: "For how long during your lifetime have you lived in a place where noise was loud enough to disturb normal conversations?" and "For how long during your lifetime have you worked at a place where noise was loud enough to disturb normal conversations?" The need to speak in raised voices might be considered a proxy for exposure to at least 66 dB. 
Air pollution exposure
The objective indicators for air pollution in the neighborhood that we had information on, were average annual all-source fine particulate matter (PM) 2.5 levels, a pollutant derived from fossil-fuel combustion, and average annual all-source benzo alpha pyrene (BaP) levels, a polycyclic aromatic hydrocarbon formed in the process of incomplete combustion of organic material. Due to insufficient resources, traffic count data could not be collected for each residential street, nor could we take field samples at residential addresses. Therefore, information on calculated levels was extracted from an official municipality source reporting PM 2.5 and BaP pollution maps.  (http://www.webcitation.org/6biSbgjwd) Pollutant dispersion modelling was done for 2011 with SELMA GIS 9.28 (Lohmeyer GmbH and Co. KG). Available data on PM 2.5 levels had been categorized: "0.0-17.5 μg/m 3", "17.5-20.3 μg/m 3", "20.3-25.0 μg/m 3", "25.0-40.0 μg/m 3", and "40.0-66.8 μg/m 3"; and for BaP levels as follows: "0.0-1.0" ng/m 3 , "1.0-3.75" ng/m 3 , "3.75-6.0" ng/m 3 , "6.0-9.75" ng/m 3 , "9.75-14.02" ng/m 3 .  The spatial resolution of the visualized pollutant levels (≈400 m × 400 m) was lower than the noise map resolution, therefore, PM 2.5 and BaP levels were still assigned to the participants even if they had reported only the name of their residential neighborhood, provided that exposure levels did not vary within the whole neighborhood.
The effect of subjective all-source air pollution in the neighborhood was assessed in order to capture a wider range of air pollutants and psychophysical aspects of air pollution. People's self-rated air pollution has its place in epidemiology when adequate assessment of long-term exposure is not possible.  Moreover, some of the deleterious effects of environmental pollution on health could be mediated through psychological mechanisms. , Perceived air pollution was derived from three 11-point-scale-based questions: "How would you rate the air quality in your neighborhood?" ("0, very low" to "10, very high", reverse coded question); "According to you, how severe is the air pollution in your neighborhood?" ("0, not at all" to "10, extremely"); and "To what extent are you annoyed by the air pollution in your neighborhood?" ("0, not at all" to "10, very much").
We anticipated high percentage of missing values on the objective noise exposure and air pollution indicators, because Bulgarians often refuse to report their residential address.  Therefore, we introduced a variable representing STI as a proxy for traffic counts, air pollution, and noise.  STI question asked: "How would you describe the road that your home is located at and its traffic? Please, rate the traffic intensity based on comparisons with other streets of Plovdiv, Bulgaria." ("very rare/no traffic", "moderately busy street", "considerably busy street", "heavy traffic", "extremely busy street/extreme traffic").
Assessment of type 2 diabetes mellitus
T2DM status was determined by self-reported doctor diagnosis: "Have you been diagnosed with any of the following diseases/conditions?" with one of the response option being "Type 2 diabetes".
Assessment of confounding factors
Additional information was collected on:
Age, gender, ethnicity, highest educational level attained, marital status, occupation, self-rated socioeconomic status.
Body Mass Index (kg/meters 2 ), sport activity (hours per week), number of wholegrain meals, alcohol beverages, meat/meat product meals, soda drinks, and sweet treats per week during the past year. When the participant had indicated the number of workouts per week rather than their overall duration, each workout was assumed to last 30 min on average.
Pack-years of smoking [(number of cigarettes a day*number of years of smoking)/20]
Family history of T2DM: "Have any of your first-degree relatives been diagnosed with any of the following diseases/conditions?", with one of the response option being "Type 2 diabetes".
Noise sensitivity was assessed with the Noise Sensitivity Questionnaire Short Form (NoiSeQSF). It measures multidimensional noise sensitivity on five subscales with three items in each. 
Duration of residence at the present address, number of nights per week sleeping with open bedroom windows, bedroom location.
Sleep disturbance was derived from an 11-point scale question: "How good is your sleep at home?" ("0, very good" to "10, cannot sleep at all").
The type of sampling was coded as well.
The dataset was screened for missing rates. The variable with the highest percentage of missing values (40.9%) was L den due to underreporting of residential addresses. Next, data were screened for univariate normality (graphical analysis, D'Agostino-Pearson K 2 test, and Jarque-Bera test), and for univariate outliers (modified outlier labelling rule).  Outliers were to be retained if they were considered honest answers or if the cause for the outliers could not be determined. Otherwise they were winsorised.
Initial descriptive statistics were computed. Likert scale-like variables with ≥5 categories were included in parametric tests.  Welch's t-test and Pearson's chi-squared test/Fisher's Exact Test/Fisher-Freeman-Halton Test were conducted to assess the association between T2DM and other variables. The main analyses were based on multivariate logistic and log-linear Poisson regressions (distribution "Poisson", link "Log", method "Hybrid", scale parameter "Fixed value", covariance matrix "Robust estimator"). , The latter is a preferable regression method when the aim is to approximate the relative risk (RR) in cross-sectional studies. However, because logistic regression is still widely used and formed the basis of our power analysis, all models are primarily reported in a logistic regression framework.
Three adjusted models were tested:
Model 1: Objectively measured noise and air pollution exposure (L den + PM 2.5 + BaP);
Model 2: Self-rated exposure to road traffic, noise, and air pollution (STI);
Model 3: Lifetime noise exposure (residential + occupational exposure to loud noise).
After inspection of the distribution of T2DM, cases across the 5-dB L den categories in the dataset (with missing data), L den was dichotomized ("71-80 dB" versus "51-70 dB") in the regression models in order to enhance the contrast between the exposure categories, to minimize exposure misclassification bias rising from the limitations of the noise map, the fact that some residential addresses were matched to the center point of the street, and the unequal distribution of T2DM cases across the 5-dB categories. PM 2.5 ("25.0-66.8 μg/m 3" versus "0.0-25.0 μg/m 3") and BaP ("6.0-14.02 ng/m 3" vs. "0.0-6.0 ng/m 3") were dichotomized as well in order to account for the small number of T2DM cases across their initial five categories. For the same reason, the categories of STI were collapsed ("very rare/no traffic" and "moderately busy street" → "low traffic exposure"; "considerably busy street" → "moderate traffic exposure"; "heavy traffic" and "extremely busy street/extreme traffic" → "high traffic exposure").
The regression models were run on 50 imputed datasets in order to increase the precision and power of the analyses. Multiple imputation yields accurate estimates in <50% missing rates. , This method has proven to be valid alternative for estimating missing rates in epidemiological surveys of environmental risk factors.  The exact adjustment set for the direct effects models was based on directed acyclic graphs (DAG) using DAGitty v. 2.3 (http://dagitty.net/). A recent simulation study showed that a priori theoretical identification of important confounders in logistic regression should always be preferred to empirical strategies.  Possible moderation effect of the type of sampling procedure was tested by computing interaction terms and evaluating them at the relaxed P < 0.2 level. Multicollinearity was tested. In sensitivity analyses, the results were stratified by participants' characteristics when it was appropriate. Post hoc statistical power was computed for the logistic regression models using G*Power v. 188.8.131.52. 
Results were considered statistically significant at P < 0.05 (two-tailed) or when 95% confidence intervals did not overlap zero. Data were processed with Statistical Package for the Social Sciences (SPSS) [International Business Machines (IBM) Corporation released 2010. IBM SPSS Statistics for Windows, Version 19.0. Armonk, New York: IBM Corporation].
| Results|| |
A total of 581 questionnaires were available for processing, but after excluding 68 questionnaires due to an unacceptably high percentage of missing data (>50%) on all key variables, or due to meeting exclusion criteria, we were left with 513 individual cases for further analyses.
Participants' mean age was 36.45 years (SD = 15.39, range: 18-83 years) and 185 participants (36.1%) were males. Most participants were ethnic Bulgarians (n = 436, 85.0%), had upper secondary (n = 257, 50.1%) or bachelor/master (n = 236, 46.0%) educational level, were married/had a spouse (n = 318, 62.0%), were employed (n = 286, 55.8%) and rated their socioeconomic status as "middle" (n = 352, 68.6%). There were 35 cases (6.8%) of T2DM in the sample.
According to official statistics for Plovdiv, Bulgaria  our sample was only partially representative of the studied population. That is, it was skewed toward women, younger people, those with middle socioeconomic status and higher education. The prevalence of diagnosed diabetes in the population is about 7.1% with 83% of those being >50 years of age and 56.7% being male.  In our sample, the prevalence of T2DM was similar, however, with a larger proportion of male cases and about 57.1% >50 years of age.
Participants differed between the two types of sampling [See supplementary Table S1 [Additional file 1]]. People in the snowball sample were younger, more likely to be female, with a lower educational level, single, and studying. The two groups did not differ significantly on T2DM prevalence, but did on road traffic and BaP exposures. Most of these differences were due to the fact that we had recruited more young people at our university campus than in the field sample, and those students had younger peers in their social network. Therefore, the type of sampling was tested for moderation effects in the multivariate models.
According to [Table 1], people with T2DM were more likely to be male and significantly older than those without. In addition, both groups differed on socioeconomic classes and occupations. Lifestyle factors such as dietary habits, alcohol consumption, and smoking were significantly associated with T2DM. Some environmental exposures were higher among those with T2DM. Correlations between the exposure indicators are given in [Supplementary Table S2 [Additional file 2]].
|Table 1: Participants' characteristics by their type 2 diabetes mellitus status|
Click here to view
L den , PM 2.5 , BaP, and type 2 diabetes mellitus
In the univariate analyses, there were higher odds for T2DM associated with L den , PM 2.5 , and BaP-7.19 (95% CI: 2.79, 18.53), 1.25 (95% CI: 0.51, 3.09), and 1.15 (95% CI: 0.50, 2.67), respectively. After these three exposure indicators were adjusted for each other and additionally controlled for important a priori confounders, the elevated odds of T2DM persisted in the multivariate model, but only the effect of L den reached statistical significance [See [Table 2]]. The model explained about 45% of the variance in T2DM, but the observed power of the test was high only for L den (100%), whereas for PM 2.5 (11%) and BaP (33%), it was low enough to account for the nonsignificant results. Interaction terms between the type of sampling and L den , PM 2.5 , or BaP were nonsignificant. The ORs associated with L den in men and women were 5.55 (95% CI: 0.93, 33.19) and 6.42 (95% CI: 0.64, 64.67), respectively. When the analysis was restricted to participants living at the address for ≥10 years the OR associated with L den remained similar (OR = 4.54, 95% CI: 1.05, 19.65) but increased for PM 2.5 (OR = 2.08, 95% CI: 0.27, 16.26) and BaP (OR = 2.27, 95% CI: 0.54, 9.46). In the Poisson regression model, the risks were relative risk (RR) = 3.08 (95% CI: 1.24, 7.62), RR = 1.29 (95% CI: 0.49, 3.41) and 1.39 (95% CI: 0.60, 3.18) for L den , PM 2.5 , and BaP, respectively.
|Table 2: Multivariate associations between Lden, PM2.5 and BaP and type 2 diabetes mellitus (logistic regression)|
Click here to view
Road traffic exposure and type 2 diabetes mellitus
People exposed to high road traffic had a crude OR = 3.60 (95 % CI: 1.57, 8.29). In the adjusted logistic model, however, high road traffic exposure was associated with nonsignificant effect (OR = 1.40, 95% CI: 0.48, 1.07) [See [Table 3]]. The model explained about 39% of the variance in T2DM and had low statistical power (15%). The odds associated with high traffic exposure among men and women were OR = 1.09 (95% CI: 0.25, 4.75) and OR = 2.18 (95% CI: 0.26, 18.55), respectively. Poisson regression model yielded RR = 1.19 (95% CI: 0.54, 2.62) associated with high traffic exposure.
|Table 3: Multivariate associations between self-reported traffic intensity and type 2 diabetes mellitus (logistic regression)|
Click here to view
Lifetime exposure to loud noise and type 2 diabetes mellitus
Only the univariate effect of lifetime occupational exposure was significant (OR = 2.13, 95% CI: 1.45, 3.12). Neither of the two lifetime exposures was associated with elevated multivariate odds of T2DM [See [Table 4]]. Stratification by gender yielded ORs of 0.63 (95% CI: 0.20, 1.94) and 1.03 (95% CI: 0.51, 2.07) in men and 0.38 (95% CI: 0.04, 3.51) and 0.29 (0.02, 5.11) in women. In the Poisson regression, the RR was 0.69 (95% CI: 0.34, 1.41) for residential and 0.91 (95% CI: 0.66, 1.27) for occupational exposure.
|Table 4: Multivariate associations between duration of lifetime residential and occupational exposure to loud noise (per 1|
interquartile range increase) and type 2 diabetes mellitus (logistic regression)
Click here to view
| Discussion|| |
To our knowledge, this was the first exploration of the feasibility of a large scale epidemiological study on traffic, noise, and air pollution as environmental risk factors for T2DM in Bulgaria. Overall, although our results are somewhat weak due to methodological limitations, they suggested that future research in the country is feasible. People exposed to L den 71-80 dB might have significantly higher odds of prevalent T2DM. Exposures to high road traffic, PM 2.5 25.0-66.8 μg/m 3 , and BaP 6.0-14.02 ng/m 3 were suspected risk factors as well, although they failed statistical significance. Lifetime exposure to loud noise was not associated with T2DM. These effects were higher among women, although nonsignificant. Given the systematic overestimation of the effects when ORs are applied as if they were RRs and the ongoing debate on the appropriateness of the two estimates; , we additionally looked at Poisson regressions, which estimated statistically significantly higher risk of T2DM associated with L den and nonsignificant risks associated with PM 2.5 and BaP. A nonsignificantly elevated risk was seen as well for participants exposed to high traffic.
The fact that our self-rated exposure indicators were not significant risk factors can partially be attributed to the low power of some tests or to a temporal instability of participants' perception of exposure. When we tested noise sensitivity and noise annoyance in a multivariate context, they were not significant or influential covariates. If noise sensitivity that has been linked to cardiovascular health  does not independently rise or modify the risk of T2DM, this would mean that highly sensitive groups will benefit just as much from physical reduction of noise exposure. On the contrary, if personal sensitivity increases the risk on its own and is not contingent on the actual exposure, then even reducing noise levels might not be enough to sufficiently decrease morbidity rates, while additional psychosocial approaches might help. These speculations prompt further scientific inquiry into the subjective exposure indicators and individual environmental sensitivity.
The RR of T2DM associated with exposure to L den 71-80 dB is lower than that estimated in a recent meta-analysis,  which might be due to the weak conclusions of that meta-analysis (due to methodological discrepancies in the primary studies), the limitations of our survey and/or the different social and environmental context in a middle-income country such as Bulgaria. The effect of PM 2.5 was positive, even given the low power of our test and it is likely to become significant in a further study as is the case in the literature.  Moreover, significant effects might have been missed due to the low spatial resolution of the air pollution maps we used. As for BaP, the mechanisms underlying its alleged effects are still unclear, although increase in adipose mass, production of reactive oxygen species, and inflammation are suspected pathways.  Other authors have found links between urinary concentrations of polycyclic aromatic hydrocarbons and diabetes, glucose homeostasis, and metabolic syndrome. ,, We cannot compare our results to previous findings because, as far as we are aware, the direct relationship between atmospheric BaP levels and T2DM has not been established. The significant and strong effect of sleeping in a bedroom with noisy façade is in line with previous work. Recently, Pirrera et al.  concluded that "averaged sound pressure level as noise indicator in a field context is insufficient to draw conclusions with respect to noise-induced sleep disturbances" and that "bedroom location is not only a major determinant for noise assessment, but also for sleep outcomes". Nighttime noise, however, was not tested in our models due to its almost perfect linear correlation with L den .
Methodological considerations and future research
Our study has several noteworthy limitations. First, its cross-sectional design precludes us from inferring causal links between exposure and outcome variables. However, as it has been argued with respect to cardiovascular effects of noise,  it is unlikely that people with T2DM will purposefully move to noisy or air polluted locations. Moreover, the effect of L den did not change in long-term residents and, due to their poor socio-economic standard, most Bulgarians are unlikely to have moved to a different address simply because of the traffic. Some of the covariates were assessed only for the last year (e.g., dietary habits), whereas others reflected lifetime exposure (e.g., pack-years of smoking). Thus, the diet, for example, might have been modified in compliance with the prescribed nonpharmacological treatment of T2DM. This would explain why people with T2DM consumed less sweet treats and soda beverages. As an attractive alternative design requiring fewer resources, we propose the cross-sectional cohort study, where a representative population is sampled cross-sectionally and then thorough exposure and outcome history is collected over a specified time period.  Regardless of the study design, however, random sampling will be essential.
Another set of inherent methodological pitfalls are sociological. While we fail to provide sufficient evidence for definitively confirming or rejecting some risk factors due to the smallish sample size, even in representative epidemiological studies, one should refrain from interpreting the 95% CIs solely in terms of statistical significance.  Furthermore, the power was sufficient for the effect of L den . Regarding missing data, imputing 50 datasets should have addressed those adequately. , The low response rate alone, on the other hand, cannot be held as an evidence of nonparticipation bias:  Bulgarians have poor survey culture and might be reluctant to participate in surveys on environmental pollutants.  As for the nonrepresentativeness of our sample, it is a common misconception that random sampling is mandatory in risk assessment; if the sole purpose is to estimate the relative risk, it is unnecessary unless the factor making the sample nonrandom acts as an effect modifier.  In the future, offering some form of financial incentives will probably increase participation rates, while public popularization of environmental health concerns might get Bulgarians more involved.
Regarding our exposure and outcome variables, there is a possibility for misclassification bias, but it is likely to be nondifferential. Noise and air pollution levels were modelled for 2009 and 2011, respectively. An update of the strategic noise map of Plovdiv, Bulgaria by recalculating and redefining the acoustic model of the street network might be necessary as well because its accuracy is currently partially outside the acceptable range defined by the Environmental Noise Directive 2002/49/EC.  The fact that only major arterial streets were sampled to calibrate the model might be skewing the exposure. On the other hand, the pathogenesis of T2DM is long, therefore, exposure levels referring to a past period are arguably more informative than current exposure. Matching addresses to the center point of the street is problematic as well as is the fact that lifetime residential history was not collected. By dichotomizing L den , PM 2.5 , and BaP, however, we should have reduced these imprecisions. Stratification of the analyses by duration of residence was not always feasible because we had very few cases in the subgroups (i.e., T2DM was reported mostly by long-term residents).
With respect to STI, it measured traffic intensity on the residential street, but for those living on small side streets close to a major road artery and vice versa adjusting for the distance to the closest major road will improve the accuracy of the models. Finally, our PM 2.5 and BaP indicators were somewhat crude proxies for air pollution, because they were derived from air pollution maps with large grid size, rather than sampled and verified at each residential address. Using readily available noise and air pollution maps was associated with some degree of exposure misclassification and is not substantiated for future studies aiming at estimation of precise exposure-response relationships, because taking all exposure levels <60 dB as one reference category loses information. They were not intended for the purposes of epidemiological research, making them unsuitable for generating high quality evidence at lower exposure levels. Enrolment of acoustic engineers and officers from the Regional Health Inspections will facilitate field sampling of noise and air pollution, calculating, and calibrating land-use regression models. Extensive occupational and residential address history should be acquired as well, and the type of domestic heating should be assessed. A look at the effects of wider range of air pollutants and further study of polycyclic aromatic hydrocarbons, by assessing both their concentrations in air and urinary metabolites, are warranted.
The fact that lifetime exposure to loud noise did not raise the risk of T2DM but L den 71-80 dB did, is probably due to an inherent limitation of our design and the inaccuracy in the self-reported exposure measure. The observed range of values on this variable and the few people ever exposed might be the causes of the flattened association with T2DM. Likewise, the question on speech disturbance made no distinction between transportation noise and other sources such as domestic equipment, toddlers, noisy neighbors, etc. As for occupational exposure to loud noise, the fact that previous studies failed to establish a link with T2DM  as well suggests that there might be a qualitative difference between noise exposure at work and at home. Provided that people typically spend more active time at work, nighttime noise mediated through sleep disturbance might, in fact, be driving the effect. Disruption of sleep patterns has been recognized as a risk factor for metabolic disorders and diabetes. , Our data additionally supports the importance of bedroom location, sleeping with open bedroom windows, and sleep disturbance as risk factors for T2DM. Therefore, occupational noise might simply be less influential because people typically do not sleep at work. Moreover, they might have developed some psychological and behavioral coping strategies in occupational settings. Finally, because we did not specifically target industrial workers, the variability in this variable might have additionally been too low.
T2DM status was assessed via self-reported doctor diagnosis but about 26.1% of the diabetes cases remain undiagnosed.  Further, we failed to assess the date of diagnosis, therefore, the temporal links with noise and air pollution cannot be disentangled. An emphasis should be put on obtaining valid biomedical information about metabolic control indicators and accompanying diseases through clinical examinations or medical documentation. The use of more sophisticated and available questionnaires in Bulgaria for the assessment of some covariates such as sleep quality,  physical activity, and diet  is recommended as well.
| Conclusions|| |
Residential noise exposure was associated with significantly higher odds of prevalent T2DM. Although there were suspected effects of other self-rated and objective exposures, they failed statistical significance. Bedroom location, sleeping in a room with open windows, and self-rated sleep disturbance were associated with elevated risks as well.
The inferences made by this small scale research and the lessons learned from its limitations could guide the designing of a population-based epidemiological survey in Bulgaria. Future environmental health research calls for acquisition of valid exposure and biomedical data, involving the community and further investigating both objective and self-rated exposure indicators. Our results might be of interest not only to Bulgarian environmentalists but to the wider international research community as well. Despite the sources of bias in our study and its crude design, it adds to the limited body of evidence of the metabolic effects of environmental noise. By providing various exposure indicators and reporting both ORs and RRs, we facilitate comparisons with previous findings. Moreover, results from South-Eastern Europe are still scant.
The authors would like to thank all the contributors who recruited other participants during the snowball sampling and all participants who made this study possible. This study received no external funding.
Financial support and sponsorship
Conflicts of interest
The authors declare that, with respect to this study, they do not have any potential conflict of interest.
| References|| |
International Diabetes Federation. IDF Diabetes Atlas Update Poster. 6 th
ed. Brussels, Belgium: International Diabetes Federation, 2013. p. 114.
Borisova A-M, Shinkov A, Vlahov J, Dakovska L, Blajeva E, Todorov T. Prevalence of diabetes mellitus and prediabetes in Bulgaria today. Endocrinologya 2012;17:172-82.
Borissova A-M, Shinkov A, Kovatcheva R, Vlahov J, Dakovska L, Todorov T. Changes in the prevalence of diabetes mellitus in Bulgaria (2006-2012). Clin Med Insights Endocrinol Diabetes 2015;8:41-5.
Atanasova E, Pavlova M, Groot W. Out-of-pocket patient payments for public health care services in Bulgaria. Front Public Health 2015;3:175.
Baillie-Hamilton PF. Chemical toxins: A hypothesis to explain the global obesity epidemic. J Altern Complement Med 2002;8:185-92.
Stefanov TS, Temelkova-Kurktschiev TS. The metabolic syndrome in Bulgaria. Folia Med (Plovdiv) 2011;53:5-14.
Eze IC, Hemkens LG, Bucher HC, Hoffmann B, Schindler C, Künzli N, et al
. Association between ambient air pollution and diabetes mellitus in Europe and North America: Systematic review and meta-analysis. Environ Health Perspect 2015;123:381-9.
Dzhambov AM. Long-term noise exposure and the risk for type 2 diabetes: A meta-analysis. Erratum. Noise Health 2015;17:123.
Noise Observation and Information Service for Europe (NOISE) [website]. Copenhagen, European Environment Agency. Available from: http://noise.eionet.europa.eu/
. [Last accessed on 2014 Jul 22].
Takuchev N, Vasileva I, Petrova S. Dispersion modeling of the air pollution, emitted by the traffic in the transport tunnel under the old town of Plovdiv, Bulgaria. Ecol Balk 2014;6:73-86.
Mathers CD, Loncar D. Projections of global mortality and burden of disease from 2002 to 2030. PLoS Med 2006;3:e442.
Dzhambov AM. Validity of self-reported traffic intensity as a proxy for road traffic counts and noise. Noise Control Engr J 2015;63:11-9.
Heidemann C, Niemann H, Paprott R, Du Y, Rathmann W, Scheidt-Nave C. Residential traffic and incidence of Type 2 diabetes: The German health interview and examination surveys. Diabet Med 2014;31:1269-76.
Dzhambov AM, Dimitrova DD. Validating a short Bulgarian version of a psychometric instrument for multidimensional noise sensitivity assessment. Folia Med (Plovdiv) 2014;56:116-125.
Cinnamon J, Schuurman N. Injury surveillance in low-resource settings using geospatial and social web technologies. Int J Health Geogr 2010;9:25.
Kamadjeu R. Tracking the polio virus down the Congo river: A case study on the use of Google Earth in public health planning and mapping. Int J Health Geogr 2009;8:4.
Agay-Shay K, Peled A, Crespo AV, Peretz C, Amitai Y, Linn S, et al
. Green spaces and adverse pregnancy outcomes. Occup Environ Med 2014;71:562-9.
European Union. 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 2002;L189:12-25.
Dzhambov AM, Dimitrova DD, Turnovska TH. Improving traffic noise simulations using space syntax: Preliminary results from two roadway systems. Arh Hig Rada Toksikol 2014;65:259-72.
Lazarus H. Prediction of verbal communication is noise - A review: Part 1. App Acoust 1986;19:439-64.
Municipality of Plovdiv. Program for reaching the normative levels of fine particulate matter under 2.5 microns (PM2,5) and polycyclic aromatic hydrocarbons (PAH) in the atmospheric air on the territory of Municipality of Plovdiv. Plovdiv: Municipality of Plovdiv; 2011.
Weiland SK, Mundt KA, Rückmann A, Keil U. Self-reported wheezing and allergic rhinitis in children and traffic density on street of residence. Ann Epidemiol 1994;4:243-7.
McCarron P, Harvey I, Brogan R, Peters TJ. Self reported health of people in an area contaminated by chromium waste: Interview study. BMJ 2000;320:11-5.
Hoaglin DC, Iglewicz B. Fine-tuning some resistant rules for outlier labeling. J Am Stat Assoc 1987;82:1147-9.
Norman G. Likert scales, levels of measurement and the "laws" of statistics. Adv Health Sci Educ Theory Pract 2010;15:625-32.
Duijnhoven RG, Knol MJ. Logistic Regression and Odds Ratios as Means to Adjust for Baseline Incomparability›s in Randomised Controlled Trials: Description of the Disadvantages, Alternatives and Frequency of Use. [Master Thesis]. Utrecht: University of Utrecht; 2012.
McNutt LA, Wu C, Xue X, Hafner JP. Estimating the relative risk in cohort studies and clinical trials of common outcomes. Am J Epidemiol 2003;157:940-3.
Choi YJ, Nam CM, Kwak MJ. Multiple imputation technique applied to appropriateness ratings in cataract surgery. Yonsei Med J 2004;45:829-37.
Schlomer GL, Bauman S, Card NA. Best practices for missing data management in counseling psychology. J Couns Psychol 2010;57:1-10.
Beyea J, Stellman SD, Teitelbaum S, Mordukhovich I, Gammon MD. Imputation method for lifetime exposure assessment in air pollution epidemiologic studies. Environ Health 2013;12:62.
Lee PH. Should we adjust for a confounder if empirical and theoretical criteria yield contradictory results? A simulation study. Sci Rep 2014;4:6085.
Faul F, Erdfelder E, Lang AG, Buchner A. GFNx01Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods 2007;39:175-91.
Dzhambov AM, Dimitrova DD. Psychometric properties of the Bulgarian translation of noise sensitivity scale short form (NSS-SF): Implementation in the field of noise control. Noise Health 2014;16:361-7.
Cook TD. Advanced statistics: Up with odds ratios! A case for odds ratios when outcomes are common. Acad Emerg Med 2002;9:1430-4.
A′Court C, Stevens R, Heneghan C. Against all odds? Improving the understanding of risk reporting. Br J Gen Pract 2012;62:e220-3.
Heinonen-Guzejev M, Vuorinen HS, Mussalo-Rauhamaa H, Heikkilä K, Koskenvuo M, Kaprio J. The association of noise sensitivity with coronary heart and cardiovascular mortality among Finnish adults. Sci Total Environ 2007;372:406-12.
Hu H, Kan H, Kearney GD, Xu X. Associations between exposure to polycyclic aromatic hydrocarbons and glucose homeostasis as well as metabolic syndrome in nondiabetic adults. Sci Total Environ 2014;505:56-64.
Khalil A, Villard PH, Dao MA, Burcelin R, Champion S, Fouchier F, et al
. Polycyclic aromatic hydrocarbons potentiate high-fat diet effects on intestinal inflammation. Toxicol Lett 2010;196:161-7.
Alshaarawy O, Zhu M, Ducatman AM, Conway B, Andrew ME. Urinary polycyclic aromatic hydrocarbon biomarkers and diabetes mellitus. Occup Environ Med 2014;71:437-41.
Pirrera S, De Valck E, Cluydts R. Field study on the impact of nocturnal road traffic noise on sleep: The importance of in-and outdoor noise assessment, the bedroom location and nighttime noise disturbances. Sci Total Environ 2014;500-501:84-90.
Babisch W, Pershagen G, Selander J, Houthuijs D, Breugelmans O, Cadum E, et al
. Noise annoyance--a modifier of the association between noise level and cardiovascular health? Sci Total Environ 2013; 452-453:50-7.
Hudson JI, Pope HG Jr, Glynn RJ. The cross-sectional cohort study: An underutilized design. Epidemiology 2005;16:355-9.
Poole C. Low P
-values or narrow confidence intervals: Which are more durable? Epidemiology 2001;12:291-4.
Kypri K, Stephenson S, Langley J. Assessment of nonresponse bias in an internet survey of alcohol use. Alcohol Clin Exp Res 2004;28:630-4.
Woodward M. Cohort studies: Analytical considerations. In: Woodward M, editor. Epidemiology: Study Design and Data Analysis. 3 rd
ed. Boca Raton, Florida: Chapman and Hall/CRC; 2013. p. 169.
Municipality of Plovdiv. Acoustic monitoring and specific traffic counting at locations used in the development of "Strategic noise map of Plovdiv agglomeration". Plovdiv: Municipality of Plovdiv; 2013.
Knutson KL, Van Cauter E. Associations between sleep loss and increased risk of obesity and diabetes. Ann N Y Acad Sci 2008;1129:287-304.
Spiegel K, Tasali E, Leproult R, Van Cauter E. Effects of poor and short sleep on glucose metabolism and obesity risk. Nat Rev Endocrinol 2009;5:253-61.
Buysse DJ, Reynolds CF 3 rd
, Monk TH, Berman SR, Kupfer DJ. The Pittsburgh sleep quality index (PSQI): A new instrument for psychiatric practice and research. Psychiatry Res 1989;28:193-213.
Atanasova V, Gatseva P, Bivolarska A, Fronas G. Body mass index and food frequency intake of foreign medical students. Trakia J Sci 2014;12(Suppl 1):367-70.
Angel M Dzhambov
PhD student, Department of Hygiene and Ecomedicine, Faculty of Public Health, Medical University of Plovdiv, No. 15-A Vassil Aprilov Boulevard, Plovdiv - 4002
Source of Support: None, Conflict of Interest: None
[Table 1], [Table 2], [Table 3], [Table 4]