Next Article in Journal
Healthcare Workers and Manual Patient Handling: A Pilot Study for Interdisciplinary Training
Next Article in Special Issue
County-Wide Mortality Assessments Attributable to PM2.5 Emissions from Coal Consumption in Taiwan
Previous Article in Journal
Prediction Model for Dry Eye Syndrome Incidence Rate Using Air Pollutants and Meteorological Factors in South Korea: Analysis of Sub-Region Deviations
Previous Article in Special Issue
Attributable Risk to Assess the Health Impact of Air Pollution: Advances, Controversies, State of the Art and Future Needs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Ultrafine Particles in Ambient Air on Primary Health Care Consultations for Diabetes in Children and Elderly Population in Ljubljana, Slovenia: A 5-Year Time-Trend Study

by
Vesna Viher Hrženjak
1,
Andreja Kukec
2,3,
Ivan Eržen
2,3 and
Dalibor Stanimirović
3,*
1
National Laboratory of Health, Environment and Food, Prvomajska 1, 2000 Maribor, Slovenia
2
Faculty of Medicine, University of Ljubljana, Vrazov trg 2, 1000 Ljubljana, Slovenia
3
National Institute of Public Health, Trubarjeva 2, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(14), 4970; https://doi.org/10.3390/ijerph17144970
Submission received: 25 May 2020 / Revised: 26 June 2020 / Accepted: 8 July 2020 / Published: 10 July 2020
(This article belongs to the Special Issue Burden of Disease Attributable to Air Pollution)

Abstract

:
Epidemiological studies indicate that exposure to ultrafine particles (UFP) in ambient air represents an important environmental public health issue. The aim of this study was to determine the association between UFP in ambient air and the daily number of consultations in the primary health care unit due to diabetes mellitus in children and elderly population of the Municipality of Ljubljana. A 5-year time-trend ecological study was carried out for the period between 1 January 2013 and 31 December 2017. The daily number of primary health care consultations due to diabetes mellitus among children and elderly population was observed as the health outcome. Daily mean UFP concentrations (different size from 10 to 100 nm) were measured and calculated. Poisson regression analysis was used to investigate the association between the observed outcome and the daily UFP, particulate matter fine fraction (PM2.5), and particulate matter coarse fraction (PM10) concentrations, adjusted to other covariates. The results show that the daily number of consultations due to diabetes mellitus were highly significantly associated with the daily concentrations of UFP (10 to 20 nm; p ≤ 0.001 and 20 to 30 nm; p ≤ 0.001) in all age groups and in the elderly population. In observed the population of children, we did not confirm the association. Findings indicate that specified environmental challenges should be addressed by comprehensive public health strategies leading to the coordinated cross-sectoral measures for the reduction of UFP in ambient air and the mitigation of adverse health effects.

1. Introduction

Air pollution and diabetes mellitus are critical public health issues worldwide [1]. Air pollution is a problem that is faced by all societies at all levels of development; especially at risk is the urban population [2]. Air pollution leads to 8.8 million premature deaths each year globally [3] and is the second leading course of deaths from noncommunicable diseases, including especially ischaemic heart disease, stroke, chronic obstructive pulmonary disease, and lung cancer [4]. Evidence is emerging on possible links between air pollution and other health outcomes, such as neurodevelopmental disorders, cognitive impairment and chronic disease conditions, such as diabetes [5]. Diabetes is a metabolic disorder caused by genetic and environmental factors, which lead to insufficient insulin secretion. The prevalence of diabetes has rapidly increased over the last few decades, and it is considered as one of the major contributors to the global burden of disease and premature deaths [6,7,8,9].
Currently, the most recognized adverse effects of air pollution on human health are due to the action of particulate matter (PM), nitrogen dioxide (NO2), and ground-level ozone (O3) [10]. PM, mainly its fine fraction (PM2.5), is a prime public health concern [4]. PM health effects have been investigated in many studies. Results from original studies have also been combined in reviews and meta-analyzes [11,12,13,14,15,16,17]. These studies have found that PM in ambient air has been linked to all-cause mortality, cardiovascular diseases, and respiratory diseases. Some population groups are more affected by air pollution than others, such as older people, children, and those with pre-existing health conditions, who are more vulnerable. Toxicological evidence suggests that PM may lead to adverse effects on human health through inflammation in the respiratory or circulatory system, systemic oxidative stress, or by activating the autonomic nervous system [18]. Over the last few years, a growing number of toxicological studies have also suggested an association between PM2.5 and diabetes through an inflammatory pathway leading to endothelial dysfunction, immune response alterations in visceral adipose tissues, and endoplasmic reticulum stress leading to alterations in insulin transduction, insulin sensitivity, and glucose metabolism [19,20,21,22,23]. A growing number of epidemiological studies have been trying to explore the associations between particulate air pollution and diabetes, and also some reviews and meta-analyzes concerning this topic have been published [1,24,25,26,27,28]. Several of them concluded that air pollution is likely a risk factor for diabetes. Current cumulative evidence appears to suggest that type 2 diabetes mellitus-related biomarkers increase with increasing exposure duration and concentrations of air pollutants [1]. However, the results from diverse epidemiological studies are not fully consistent, and no systematic reviews have yet been published to determine whether these claims are substantiated and verified.
In the toxicity of particles, size and morphology play a key role. Division of particle size most commonly draws a distinction between coarse and fine particles—fine particles (PM2.5) are smaller than 2.5 µm (AED ≤ 2.5 μm), and coarse particles (PM10) are larger (AED ≤ 10 μm). Ultrafine particles (UFP) are defined as particles smaller than 0.1 μm (AED ≤ 0.1 μm) [29]. In the particle measurements, they are the prime contributor to the number concentrations with more than 80% of the total airborne particle number in urban and industrial environments, while they contribute little to the mass [29,30,31]. They are mainly generated by road traffic, the burning of fossil fuels and biomass, and by industrial emissions, but can also originate from natural sources [29,32,33,34].
Due to their small size, UFP are believed to exert higher toxicity than larger particles. UFP have a higher surface area per unit mass. Large relative surface area absorbs relatively higher amounts of toxic materials increasing the exchange of toxic material with the physiological medium. The lower aerodynamic diameter of UFP allows them to penetrate into the deepest parts of the respiratory tract (alveoli) where there are no flagella, and where the exchange of gases between the inhaled air and the blood takes place. They are eliminated less efficiently from the respiratory tract, which increases their residence time and facilitates their accumulation. After that, they cross the alveolar barrier, enter the bloodstream, and they can travel throughout the body, passing even through the placenta [29,35,36].
A growing number of studies evaluate the role of UFP exposure on adverse health outcomes. Most of what we know about possible health effects from UFP exposure comes from research done on animals [29]. Strong and consistent evidence from animal studies indicates long-term exposure to UFP is related to negative effects on the brain, nervous, and respiratory system [29,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55]. A few animal studies have reported effects on reproduction and development [29,56,57,58,59,60]. The body of epidemiological literature on the association between health and UFP exposure is limited but growing [29,61,62,63]. Increasing evidence suggests UFP can cause adverse health outcomes, especially on respiratory and cardiovascular morbidity and mortality. The evidence suggests adverse short-term associations with inflammatory and cardiovascular changes, which may be at least partly independent of other pollutants [62], but it is not clear whether short- and long-term exposure in humans can induce adverse health effects to the same extent as those observed in animals [32,36]. For other studied health outcomes, such as the diabetes mellitus, the evidence on independent health effects of UFP remains inconclusive or insufficient [62], which is related to the fact that actual data on UFP exposure are rather undependable. The aim of our study was to determine the association between UFP in ambient air and the daily number of consultations in the primary health care unit due to diabetes mellitus in children and elderly population of the municipality of Ljubljana (MOL), adjusted for other co-pollutants, seasonal, and meteorological factors.

2. Materials and Methods

2.1. Study Design, Period, and Area of Observation

Ecological time-trend design was used to assess the association between UFP in ambient air and the daily number of consultations in the primary health care unit due to diabetes mellitus in children and the elderly population of MOL. The unit of observation was a single day of the observed period between 1 January 2013 and 31 December 2017.
The study was carried out in an administrative unit of the MOL, which is the largest city and capital of Slovenia, with a population between approximately 282,741 and 288,250 people over the period 2013–2017 (Figure 1). It has an area of 275 km2 and its population density varies from 1028.1 to 1048.2 people per km2 over the period 2013–2017 [64]. It is located at 46°03′20″ N 14°30′30″ E and has an average elevation of 295 m. It has a continental climate (bordering on a subtropical humid climate), which is characterized by moderately cold winters and warm summers. The average yearly temperature in 2018 was 12.5 °C (January 2018 average: 4.8 °C and July 2018 average: 22.3 °C). Yearly precipitation in 2018 was approximately 1400 mm and the average annual percentage of humidity was 75.0%. The economic activity in the MOL is very dynamic, with around 50,000 companies operating mainly in the pharmaceutical, petrochemical, and food industries. Other economic activities include banking, finance, transport, construction, crafts, services, and tourism. Due to air quality perspective, MOL is a highly polluted municipality in Slovenia with PM10. The reasons are geographical characteristics (poorly ventilated basin), meteorological conditions (temperature inversions are frequently present in winter and autumn season), and largely road transport factors, especially in its narrower urban area with heavy traffic [65].

2.2. Observed Population

The observed population consisted of all individuals, residing in MOL temporarily and permanently, who visited the primary health care units of the Community Health Center (CHC) Ljubljana due to diabetes mellitus in the period of observation.

2.3. Data Collection

2.3.1. Observed Health Outcomes

Routinely collected health data were obtained from the health information system of the CHC Ljubljana. The daily number of primary health care consultations due to diabetes mellitus (WHO International Classification of Diseases, version 10: E10-E14) [66] among all observed populations, children (from birth to 10 years old) and the elderly (60 year +) population, were observed as the health outcomes.

2.3.2. Observed UFP and PM Data and Meteorological Data

Within the project “Ultrafine particles—an evidence based contribution to the development of regional and European environmental and health policy” (UFIREG) UFP number concentrations (PNC) measurements in different size ranges from 10 until 800 nm using custom-made mobility particle size spectrometers in five cities, including Ljubljana, were conducted. In Ljubljana UFP number concentrations were collected on a measuring site located in the urban environment of MOL (Hajdrihova street 19, Ljubljana). Kindergarten, primary school and multistory buildings were situated in the vicinity of the measuring site. The measuring site was situated 50 m from the nearest low-frequented road, but it was not in the direct vicinity of highly-frequented roads. The railway station is situated near the measuring site (approximately 1 km). The measuring site was considered as representative for the substantial part of the population of MOL by the UFIREG study group. The project started in July 2011 and ended in December 2014 [67]. After the end of the project UFP number concentrations measurements in Ljubljana were continued for another 5 years. UFP number concentrations were measured in 5 min time intervals and then calculated to hourly and daily average concentrations for five size ranges: UFP size from 0.01 to 0.02 μm, from 0.02 to 0.03 μm, from 0.03 to 0.05 μm, from 0.05 to 0.07 μm, and from 0.07 to 0.10 μm as particle number/cm3. For this study, data on daily average concentrations for the outlined five size ranges between 1 January 2013 and 31 December 2017 were used.
Daily average PM10 and PM2.5 concentrations (in µg/m3), the daily mean value of air temperature (in °C), and relative humidity (as %) were obtained at the fixed measuring station in Ljubljana, which is a part of the National automated network for monitoring air quality. The data were obtained in agreement with the Slovenian Environment Agency [68].

2.3.3. Data on Seasonal Factors and Pollen Concentration

The following seasonal factors were considered in the analysis: season of the year (spring, summer, autumn, winter), working day (yes-weekday/no-weekend), and influenza season (yes/no). Data on influenza season were obtained from annual reports (Epidemiological surveillance of communicable diseases in Slovenia) of the National Institute of Public Health of the Republic of Slovenia [69].
Data on pollen concentration (grains of pollen/m3 over a 24-h period) of different allergenic plants (cypress family, grass family, ragweed, birch) were obtained from measuring the station in Ljubljana for the period of observation. The data were acquired from a 2-h collection of pollen concentrations from the database located at the National Laboratory of Health, Environment and Food [70].

2.4. Statistical Analysis

2.4.1. Data Description

The distributions of observed health outcomes, selected air pollution data, and meteorological parameters were statistically described by non-parametric typical statistical values (mean, standard deviation, minimum, maximum, 1st, median, and 3rd quartile). The temporal variability of observed health outcomes, selected air pollution data and pollen concentration were presented in sequence charts [71].

2.4.2. Correlation Analysis

The Pearson correlation coefficient was used to assess the ecological correlation a) between exposure (daily average (24-h) UFP (five size ranges), average 24-h PM10 and PM2.5 concentrations), and outcome variables (daily number of consultations due to diabetes mellitus, separate for three age groups—all age groups, children, and the elderly), and b) between the daily average 24-h number concentration of UFP and the average 24-h concentration of PM10 and PM2.5 and observed meteorological conditions.

2.4.3. Temporal Relationship Analysis

In the temporal relationship analysis, the daily number of consultations due to diabetes mellitus (separate for all age groups, children, and the elderly) were considered as the observed health outcomes. As the explanatory factors, the daily average 24-h number concentration of UFP and the average 24-h concentration of PM10 and PM2.5 were considered. Data on seasonal factors, pollen concentration, and meteorological parameters were considered as covariates. Respectively, lags from zero up to five days from exposure to the consultation were examined to determine the amount of time between exposure and effect for observed pollutants in ambient air and pollen data.
Poisson regression analysis was used to investigate the association between the observed outcome and the daily UFP, PM2.5, and PM10 concentrations, adjusted to other covariates. Poisson regression is a regression technique used when the outcome is a count variable or a rate [72].
The modelling procedure was performed in two stages. In the first stage, non-adjusted models were built by only relating the observed outcomes and explanatory factors. In the second stage, adjusted models were built by adding the explanatory factor to a core covariate model that included seasonal parameters, meteorological data, and pollen data. The interpretable end-result was the incidence rate ratio (IRR). It was presented together with its lower and upper 95% confidence interval (CI). A p-value of 0.05 or less was considered as statistically significant in all statistical tests. All statistical analyses were carried out by using the SPSS software platform, version 25.0 (SPSS Inc., Chicago, IL, USA).
The protocol of the study was approved by the National Medical Ethics Committee (0120-109/2017-5 KME 71/03/17, 2017). In addition, consent from the Community Health Center Ljubljana for using data was obtained (820-1-4/2017, 2017). All data were depersonalized and provided for use in this study in an aggregated form.

3. Results

3.1. Data Description

The daily number of consultations due to diabetes mellitus: all age groups, children, and the elderly were available for all 1.826 days of the study period. There were 373/1.826 (20.4%) days with no consultations for observed health outcomes. More detailed information is available in Table 1. The temporal variability of daily consultations for observed health outcomes is presented with sequence diagrams in Figure 2a–c.
The average 24-h number concentration of UFP data for individual size ranges (UFP0.01–0.02, UFP0.02–0.03, UFP0.03–0.05, UFP0.05–0.07, and UFP0.07–0.10) in the observed period was measured for 1.655/1.826 (90.6%) days. The average 24-h concentration of PM10 and PM2.5 in the observed period was available for 1.765/1.826 (96.6%) days and for 1.742/1.826 (95.4%) days. More detailed information is available in Table 1. Temporal variability of daily UFP and PM10 and 2.5 data is presented with sequence diagrams in Figure 3a,b. In all size ranges of UFP and PM10 and 2.5, the daily number of concentrations were highest in the winter season.
Complete data for the daily meteorological factors (air temperature and relative humidity) were available for all 1.826 days of the study period (Table 1). Daily average pollen concentration was the highest in the spring season (Figure 3c).

3.2. Correlation Analysis

The Pearson correlation coefficient between average 24-h UFP0.01–0.02 concentrations and daily number of consultations due to diabetes mellitus in the total observed population (separate for all age groups, children, and the elderly) ranged from 0.033 to 0.212. With the exception of children (p = 0.183), a statistically significant correlation was confirmed for all age groups and the elderly (p ≤ 0.001).
The Pearson correlation coefficient between average 24-h UFP0.02–0.03 concentrations and the daily number of consultations due to diabetes mellitus in the total observed population (separate for all age groups, children, and the elderly) ranged from 0.047 to 0.254. With the exception of children (p = 0.058), a statistically significant correlation was confirmed for all age groups and the elderly (p ≤ 0.001).
The Pearson correlation coefficient between average 24-h UFP0.03–0.05 concentrations and daily number of consultations due to diabetes mellitus in the total observed population (separate for all age groups, children, and the elderly) ranged from 0.054 to 0.265. A statistically significant correlation was confirmed for all age groups, the elderly (p ≤ 0.001), and children (p = 0.029).
The Pearson correlation coefficient between average 24-h UFP0.05–0.07 concentrations and the daily number of consultations due to diabetes mellitus in the total observed population (separate for all age groups, children, and the elderly) ranged from 0.034 to 0.187. With the exception of children (p = 0.171), a statistically significant correlation was confirmed for all age groups and the elderly (p ≤ 0.001).
The Pearson correlation coefficient between average 24-h UFP0.07–0.10 concentrations and the daily number of consultations due to diabetes mellitus in the total observed population (separate for all age groups, children, and the elderly) ranged from 0.010 to 0.114. With the exception of children (p = 0.688), a statistically significant correlation was confirmed for all age groups and the elderly (p ≤ 0.001).
The Pearson correlation coefficient between average 24-h PM2.5 concentrations and the daily number of consultations due to diabetes mellitus in the total observed population (separate for all age groups, children, and the elderly) ranged from 0.006 to 0.030. No statistically significant correlation for all age groups (p = 0.209), the elderly (p = 0.222), and children (p = 0.817) was found.
The Pearson correlation coefficient between average 24-h PM10 concentrations and the daily number of consultations due to diabetes mellitus in the total observed population (separate for all age groups, children, and the elderly) ranged from 0.000 to 0.036. No statistically significant correlation for all age groups (p = 0.126), the elderly (p = 0.143), and children (p = 0.996) was found.
The correlation analysis showed a statistically significant correlation between the daily average 24-h temperature and the daily average 24-h number concentration of UFP0.02–0.03 (r = −0.229, p ≤ 0.001), UFP0.03–0.05 (r = −0.295, p ≤ 0.001), UFP0.05–0.07 (r = −0.372, p ≤ 0.001), and UFP0.07–0.10 (r = −0.438, p ≤ 0.001), and the average 24-h concentration of PM10 (r = −0.482, p ≤ 0.001) and PM2.5 (r = −0.557, p ≤ 0.001). The correlation between the daily average 24-h temperature and the daily average 24-h number concentration of UFP0.01–0.02 was not statistically significant (r = −0.014, p = 0.577).
The daily average 24-h number concentration of UFP0.01–0.02 (r = −0.131, p ≤ 0.001), UFP0.03–0.05 (r = 0.091, p ≤ 0.001), UFP0.05–0.07 (r = 0.183, p ≤ 0.001), and UFP0.07–0.10 (r = 0.249, p ≤ 0.001), and the average 24-h concentration of PM10 (r = 0.230, p ≤ 0.001) and PM2.5 (r = 0.296, p ≤ 0.001) have a statistically significant correlation with the daily average 24-h relative humidity. We did not confirm correlation between the daily average 24-h number concentration of UFP0.02–0.03 and the daily average 24-h relative humidity (r = 0.020, p = 0.415).

3.3. Relationship Analysis

The results of the non-adjusted Poisson regression models showed that the daily number of consultations due to diabetes mellitus by all age groups were statistically significantly associated with UFP0.01–0.02_lag0 and UFP0.02–0.03_lag0 (p ≤ 0.001) and with PM2.5_lag2 (p ≤ 0.001) and with PM10_lag0 (p ≤ 0.001). The results of the adjusted Poisson regression models showed that the daily number of consultations due to diabetes mellitus by all age groups were statistically significantly associated with UFP0.01–0.02_lag0 and UFP0.02–0.03_lag0 (p ≤ 0.001). More detailed information is available in Table 2.
In the observed children population, we confirmed a statistically significantly association between the daily number of consultations due to diabetes mellitus and the daily concentration of UFP0.02–0.03_lag0 (p = 0.045) and UFP0.03–0.05_lag0 (p = 0.023). The association was confirmed in case of non-adjusted Poisson regression models (Table 3). With adjusted Poisson regression models, we did not find a statistically significantly association between the daily number of consultations due to diabetes mellitus and the daily concentration of UFP, PM2.5 or PM10 (Table 3).
In the observed elderly population, we found a statistically significantly association between the daily number of consultations due to diabetes mellitus and the 24-h number concentration of UFP0.01–0.02_lag0 and UFP0.02–0.03_lag0 (p ≤ 0.001) and with PM2.5_lag2 (p ≤ 0.001) and with PM10_lag0 (p ≤ 0.001). The association was confirmed in case of non-adjusted Poisson regression models. The results of adjusted Poisson regression models showed that the daily number of consultations due to diabetes mellitus by the elderly population were statistically significantly associated with UFP0.01–0.02_lag0 and UFP0.02–0.03_lag0 (p ≤ 0.001). More detailed information is available in Table 4.

4. Discussion

In this study, we investigated the association between daily UFP concentrations in the ambient air and the daily number of consultations due to diabetes mellitus in primary health care units of the CHC in the MOL. The results of the study showed that the daily number of consultations due to diabetes mellitus were highly significantly associated with the daily concentrations of UFP (10 to 20 nm; p ≤ 0.001 and 20 to 30 nm; p ≤ 0.001) on the day of exposure in all age groups and in the elderly population. In the observed population of children, the association was not statistically significant. The reason for this is most probably due to the fact that children suffer mainly from Type 1 diabetes, which is genetically determined. In adults, who suffer mainly from Type 2 diabetes, lifestyle and other factors (overweight, low physical activity) play an important role in the development of diabetes. However, as already put forward, the study of the impacts of environmental factors is becoming increasingly important in the field. The associations were confirmed in adjusted models that included UFP number concentrations for five different size ranges, PM2.5 and PM10, seasonal parameters, meteorological data, and pollen data. Our findings indicated that the effects of UFP in the ambient air on the daily number of consultations due to diabetes mellitus are independent from those caused by other air pollutants, seasonal, or meteorological conditions.
Despite evidence regarding associations between PM2.5 and diabetes [19,20,21,22,23,24,25,26,27,28], a limited number of studies have investigated effects of UFP on diabetes mellitus [73,74,75,76]. The UFIREG study from Central Europe found associations between short-term exposure to UFP, PNC, and PM and hospital admissions caused by diabetes [73]. Due to their small size, UFP are believed to exert higher toxicity than larger particles and are associated with effects not only restricted to the lung. Three biological pathways that are assumed not to act independently are discussed. Inhalation of UFP may lead to the release of proinflammatory mediators or vasoactive molecules from lung cells causing systemic oxidative stress and inflammation, UFP deposited in the pulmonary tree may be associated with imbalance of the autonomic nervous system, and UFP and their constituents can be translocated into the blood causing endothelial dysfunction. Systemic oxidative stress and inflammation, imbalance of the autonomic nervous system, and endothelial dysfunction can lead to insulin resistance and, therefore, can promote the progression of type 2 diabetes [73]. A population-based cohort study in Toronto, Canada found each interquartile change in exposure to UFP was associated with an increased risk of diabetes [74]. One epidemiology study reported a short-term association between the prediabetic marker HbA1c and indoor, but not outdoor, UFP exposure [75]. A study of the association of long-term near-highway exposure to UFP with cardiovascular diseases, diabetes, and hypertension did not observe the association of TAA-PNC (time activity adjusted annual average UFP exposure) with diabetes [76]. These inconsistent conclusions might be explained by the subgroup of populations being targeted and the different methodologies used. The characteristic of these research studies is that they mainly covered a short observed period and the adjusted covariates were few. They examined the association between mortality and hospital admissions due to selected disease groups and UFP ambient air pollution. To our knowledge, our research is the first that examined the association between the number of first visits to the doctor at the primary health care level due to selected groups of diseases and UFP ambient air pollution. Consequently, this study addressed another important issue, namely the occurrence of less pronounced changes in population health associated with UFP ambient air pollution that do not affect mortality or require hospital admission. Further high-quality studies are needed to expand and validate the conclusions.

4.1. Limitations of the Study

The conducted study has some potential limitations. The possibility to establish associations between different air pollutants (including UFP) and health outcomes depends on how well the exposure can be quantified. The most obvious limitation of this study is that measured concentrations are obtained from one measuring station, which potentially creates exposure misinterpretation concerning the level of exposure for some individuals among the population. In this study, we used one single measurement site, assuming temporal variability at this site reflected the temporal variability of the exposure in the population, which was sufficient in order to estimate short-time variations in pollutant concentrations in time-series studies [77,78,79]. In addition, UFP have been shown to have a higher spatial variability than other air pollutants, therefore additional exposure misclassification might be assumed. We are aware that many factors are constantly present that affect the emissions, transformation, and dispersion of pollutants in ambient air, and thus contribute to the spatial variability that exists within urban areas. However, a study in Dresden, Germany, showed low spatial variability in UFP among urban background stations and suggested that short-time UFP exposure for the average population might be adequately represented by one measurement site if the location was chosen carefully [79]. Therefore, we assumed that the data obtained from a central measuring point, due to the location at the bottom of the MOL basin and in the city center, are a satisfactory approximation of the concentration of pollutants in ambient air for the entire urban area of MOL.
A frequent problem in using routinely collected environmental data for time-trend studies is a lack of a certain share of the measurements. In our study only a small amount of data were missing for all measured air pollutants (from 9.4% for UFP data to 3.4% for PM10 data). Missing data were mainly due to the calibration of instruments or filter blockage of the measuring device.
In the toxicity of PM, the size plays a key role; nonetheless, previous studies suggested that the effects of PM on diabetes mellitus might differ with chemical composition as well [80,81,82,83]. This could be an additional limitation of the study since only the size of PM was observed, and their chemical composition was not considered. These issues certainly raise questions for future researchers on the effects of UFP on public health and suggest that forthcoming experiments should address the chemical composition of UFP.
When using medical data, we have to take into account that the daily number of consultations due to selected diagnosis might be misestimated as well. On one hand, we assume that most subjects have a chosen personal physician near their place of residence. We are aware that this is not always the case. Therefore, individuals who have permanent residence in the MOL, but have a selected personal physician outside MOL, were excluded from the research. On the other hand, an additional limitation of such collected health data is the fact that some permanent residents of the MOL who have also a selected personal doctor there, due to daily or even weekly migrations can spend most of the time outside their place of residence, and therefore their exposure to ambient air pollutants can be different. Another limitation can be related to the relatively small population and consequently less daily number of consultations due to diabetes mellitus in primary health care units of the CHC.
As in other similar studies where data from the medical records were used, we also encountered a certain degree of data inaccuracy, which was related to coding. Some similar problems in using medical data in Slovenia were observed in a project finished in 2012 [84]. On one side, there were no problems with the completeness of data collection. One the other side, there may be differences in coding the primary diagnosis in different primary health care units of the CHC, and since only the primary diagnosis was considered, the number of cases might be misestimated.
One could argue that not all steps in a time series regression study proposed by Bhaskaran and co-authors [85] were done. However, to a certain extent we have followed these steps, while the others (a) periodic functions of time (sine/cosine functions); (b) flexible spline functions of time; (c) consider possible non-linear associations as in other regression contexts; (d) diagnostic plots based on deviance residuals; and (e) multiple sensitivity analyses changing key modelling decisions, were out of the scope of this study and were planned to be addressed in further analyses.
Despite certain limitations, this study is surely relevant, because the association between the daily UFP concentration and the daily number of consultations due to diabetes mellitus in CHC was not yet investigated in Slovenia. Further epidemiological studies are needed to establish the exact exposure assessment and to explore the impact of UFP on specific population groups, especially the more vulnerable ones and those that are more exposed.

4.2. Strenghts of the Study

In contrast, our study has several strengths. First, we had an extensive dataset with continuous measurements of many air pollutants and an extensive dataset about the daily number of consultations in primary health care units of the CHC Ljubljana. This allowed us to study the associations of UFP with different health outcomes in adjusted models that included UFP number concentrations for five different size ranges, PM2.5 and PM10, seasonal parameters, meteorological data, and pollen data. It was complex, extensive, and demanding research, due to the specific environment and a large number of factors that can critically affect the results of the research, and the large amount of data that needed to be processed and analyzed. Particularly it was very challenging to obtain a large amount of health data from the information system of the CHC. This action included the preparation of the command script for the extraction of data on the number of consultations, including selected diagnoses and the age of patients for a five-year period, and later the editing, cleaning, and consolidation of data by forming age classes. Large amounts of data made it possible to conclude that the effects of UFP in the ambient air on the daily number of consultations due to diabetes mellitus were independent from those caused by other air pollutants, seasonal, or meteorological conditions.
Second, the study is one of few time-trend studies in Slovenia which linked routinely collected health and environmental data with the aim of examining whether there was an association between the number of first visits to a doctor at the primary health care level and ambient air pollution [86,87,88,89].
Furthermore, this is the first study where the effects of UFP on diabetes were studied over a period of 5 years. Several studies already concluded that air pollution is likely a risk factor for diabetes. In Slovenia, as in many countries, we are seeing an increase in the prevalence of diabetes, which has become an important public health problem with accompanying morbidity and consequences for patients, families, and the community. Compared to European and OECD countries, Slovenia is in the middle zone in terms of the prevalence of diabetes. The European Health Interview Survey (EHIS) conducted in 2014 estimated the prevalence of diabetes among the adult population at 6.9%. The same result was also shown in 2016 by the Health-Related Behavioral Style Survey [90]. According to the International Diabetes Federation (IDF) Atlas, the national prevalence in adults in 2019 in Slovenia was 7.8%, and it was similar in neighboring countries (Italy 8.3%, Hungary 9.3%, Austria 9.7%, Croatia 6.8%) and Europe (8.9%) [91]. Assuming that (similarly to European countries) in Slovenia about 15% of people with diabetes are treated without medication, i.e., only with adequate physical activity, diet, and weight loss, it is estimated that in Slovenia in 2017 about 131,000 people were diagnosed with diabetes. The number of all people with diabetes (identified and undiagnosed) could range from 154,000 to 175,000 [90].
The results of the study could provide a valuable platform for formulating recommendations to environmental and health policymakers in Slovenia. Findings imply the adoption of an appropriate legislative framework for routine UFP measurement in the national network for ambient air quality and conveying this information to competent environmental and health institutions, and the general population. Extensive dissemination of the results and their suitable interpretation might be used for raising the awareness of the professional and lay public about the health risks associated with the presence of UFP in ambient air and the possible implementation of preventive and mitigation measures, both at the level of the population and the individual.

5. Conclusions

The aim of this study was to determine the association between UFP in ambient air and the daily number of consultations in the primary health care unit due to diabetes mellitus in all age groups, in children, and the elderly population of the MOL. The results of the study showed that the daily number of consultations due to diabetes mellitus were highly significantly associated with the daily concentrations of UFP (10 to 20 nm; p ≤ 0.001 and 20 to 30 nm; p ≤ 0.001) on the day of exposure in all age groups and in the elderly population. In the observed population of children, the association was not statistically significant. Our findings indicated that effects of UFP in the ambient air on the daily number of consultations due to diabetes mellitus were independent from those caused by other air pollutants, seasonal, or meteorological conditions.
The presented study addresses some important challenges in the field of public health, which have obviously been rather disregarded in view of the research carried out so far and underrated in terms of the enacted environmental health policies by Slovenia and other countries. We have established an association between the daily UFP (10 to 30 nm) concentration and the daily number of consultations for diabetes among all age groups and the elderly population. In view of that, the results of the study indicate that there is a compelling necessity for coordinated public health activities in sense of the reduction of harmful environmental factors in the MOL and the broader region. Unfortunately, the consideration of this issue in previous campaigns of the national, as well as regional or municipal, authorities has been partial and without wider socio-economic measures, such as the actions for the prevention of air pollution and the mobilization of the public for this cause.
It is clear that isolated measures of close social communities (countries, regions, cities, etc.), even if successful, will not be able to eliminate the negative effects of UFP on health. Air pollution is just one of the detrimental global consequences of uncontrolled industrialization, the lack of a sustainable model for economic development, excessive consumption, and the general ignorance of a large part of the population. All indications are that the current socio-economic arrangements concerning environmental issues urgently need radical redefinition. Some events in recent years suggest that this shift has already begun in individual segments of global society and may indicate an alternative path to less air pollution and better public health.

Author Contributions

Conceptualization, V.V.H., I.E. and D.S.; Data curation, A.K. and D.S.; Formal analysis, A.K. and D.S.; Investigation, V.V.H. and I.E.; Methodology, V.V.H., A.K., I.E. and D.S.; Software, D.S.; Supervision, I.E. and D.S.; Validation, A.K., I.E. and D.S.; Visualization, A.K.; Writing—original draft, V.V.H., A.K., I.E. and D.S.; Writing—review and editing, V.V.H., A.K., I.E. and D.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Authors of this study would like to thank the members of UFIREG project, colleagues from the Slovenian Environment Agency (ARSO), the National Laboratory of Health, Environment and Food (NLZOH), and the Community Health Center Ljubljana for the provision of environmental, meteorological, and health data, and for their help and support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, Y.; Xu, L.; Shan, Z.; Teng, W.; Han, C. Association between air pollution and type 2 diabetes: An updated review of the literature. Ther. Adv. Endocrinol. Metab. 2019, 10, 1–15. [Google Scholar] [CrossRef] [PubMed]
  2. WHO. Air Pollution. Available online: https://www.who.int/health-topics/air-pollution#tab=tab_1 (accessed on 13 May 2020).
  3. Lelieveld, J.; Klingmüller, K.; Pozzer, A.; Pöschl, U.; Fnais, M.; Daiber, A.; Münzel, T. Cardiovascular disease burden from ambient air pollution in Europe reassessed using novel hazard ratio functions. Eur. Heart. J. 2019, 40, 1590–1596. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. WHO Regional Office for Europe. Noncommunicable diseases and air pollution. In Proceedings of the WHO European high-level Conference on Noncommunicable diseases, Ashgabat, Turkmenistan, 9–10 April 2019.
  5. WHO. Review of Evidence on Health Aspects of Air Pollution—REVIHAAP Project; Technical report; WHO Regional Office for Europe: Copenhagen, Denmark, 2013; Available online: http://www.euro.who.int/__data/assets/pdf_file/0004/193108/REVIHAAP-Final-technical-report-final-version.pdf?ua=1 (accessed on 13 May 2020).
  6. Polonsky, K.S. The past 200 years in diabetes. N. Engl. J. Med. 2012, 367, 1332–1340. [Google Scholar] [CrossRef]
  7. Nolan, C.J.; Damm, P.; Prentki, M. Type 2 diabetes across generations: From pathophysiology to prevention and management. Lancet 2011, 378, 169–181. [Google Scholar] [CrossRef]
  8. Balakumar, P.; Maung, U.K.; Jagadeesh, G. Prevalence and prevention of cardiovascular disease and diabetes mellitus. Pharm. Res. 2016, 11, 600–609. [Google Scholar] [CrossRef] [PubMed]
  9. Forouzanfar, M.H.; Afshin, A.; Alexander, L.T.; Andersn, H.R.; Bhutta, Z.A.; Biryukov, S.; Brauer, M.; Burnett, R.; Cercy, K.; Charlson, F.J.; et al. Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990–2015: A systematic analysis for the Global Burden of Disease Study 2015. Lancet 2016, 388, 1659–1724. [Google Scholar] [CrossRef] [Green Version]
  10. Air Quality in Europe—2019 Report; EEA Report No 10/2019; European Environment Agency: Copenhagen, Denmark, 2019; Available online: https://www.eea.europa.eu/publications/air-quality-in-europe-2019 (accessed on 13 May 2020).
  11. Hoek, G.; Krishnan, R.M.; Beelen, R.; Peters, A.; Ostro, B.; Brunekreef, B.; Kaufman, J.D. Long-term air pollution exposure and cardio-respiratory mortality: A review. Environ. Health 2013, 12. [Google Scholar] [CrossRef] [Green Version]
  12. Lu, F.; Xu, D.; Cheng, Y.; Dong, S.; Guo, C.; Jiang, X.; Zheng, X. Systematic review and meta-analysis of the adverse health effects of ambient PM2.5 and PM10 pollution in the Chinese population. Environ. Res. 2015, 136, 196–204. [Google Scholar] [CrossRef]
  13. Mustafic, H.; Jabre, P.; Caussin, C.; Murad, M.H.; Escolano, S.; Tafflet, M.; Perier, M.-C.; Marijon, E.; Vernerey, D.; Empana, J.-P.; et al. Main air pollutants and myocardial infarction: A systematic review and meta-analysis. J. Am. Med. Assoc. 2012, 307, 713–721. [Google Scholar] [CrossRef]
  14. Shah, A.S.; Langrish, J.P.; Nair, H.; McAllister, D.A.; Hunter, A.L.; Donaldson, K.; Newby, D.E.; Mills, N.L. Global association of air pollution and heart failure: A systematic review and meta-analysis. Lancet 2013, 382, 1039–1048. [Google Scholar] [CrossRef] [Green Version]
  15. Raaschou-Nielsen, O.; Andersen, Z.J.; Beelen, R.; Samoli, E.; Staffogia, M.; Weinmayr, G.; Hoffmann, B.; Fischer, P.; Nieuwenhuijsen, M.J.; Brunekreef, B.; et al. Air pollution and lung cancer incidence in 17 European cohorts: Prospective analyses from the European Study of Cohorts for Air Pollution Effects (ESCAPE). Lancet Oncol. 2013, 14, 813–822. [Google Scholar] [CrossRef]
  16. Cesaroni, G.; Forastiere, F.; Stafoggia, M.; Andersen, Z.J.; Badaloni, C.; Beelen, R.; Caracciolo, B.; de Faire, U.; Erbel, R.; Eriksen, K.T.; et al. Long term exposure to ambient air pollution and incidence of acute coronary events: Prospective cohort study and meta-analysis in 11 European cohorts from the ESCAPE Project. BMJ 2014, 348, f7412. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Beelen, R.; Raaschou-Nielsen, O.; Stafoggia, M.; Andersen, Z.J.; Weinmayr, G.; Hoffmann, B.; Wolf, K.; Samoli, E.; Fischer, P.; Nieuwenhuijsen, M.; et al. Effects of long-term exposure to air pollution on natural-cause mortality: An analysis of 22 European cohorts within the multicentre ESCAPE project. Lancet 2014, 383, 785–795. [Google Scholar] [CrossRef]
  18. Brook, R.D.; Rajagopalan, S.; Pope, C.A.; Brook, J.R.; Bhatnagar, A.; Diez-Roux, A.V.; Holguin, F.; Hong, Y.; Luepker, R.V.; Mittleman, M.A.; et al. Particulate matter air pollution and cardiovascular disease an update to the scientific statement from the American Heart Association. Circulation 2010, 121, 2331–2378. [Google Scholar] [CrossRef] [Green Version]
  19. Rajagopalan, S.; Brook, R.D. Air pollution and type 2 diabetes: Mechanistic insights. Diabetes 2012, 61, 3037–3045. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, C.; Ying, Z.; Harkema, J.; Sun, Q.; Rajagopalan, S. Epidemiological and experimental links between air pollution and type 2 diabetes. Toxicol. Pathol. 2013, 41, 361–373. [Google Scholar] [CrossRef] [Green Version]
  21. Liu, C.; Xu, X.; Bai, Y.; Tse-Yao, W.; Rao, X.; Wang, A.; Sun, L.; Ying, Z.; Gushchina, L.; Maiseyeu, A.; et al. Air pollution-mediated susceptibility to inflammation and insulin resistance: Influence of CCR2 pathways in mice. Environ. Health Perspect. 2014, 122, 17–26. [Google Scholar] [CrossRef] [Green Version]
  22. Zheng, Z.; Xu, X.; Zhang, X.; Wang, A.; Zhang, C.; Hüttemann, M.; Grossman, L.I.; Chen, L.C.; Rajagopalan, S.; Sun, Q.; et al. Exposure to ambient particulate matter induces a NASH-like phenotype and impairs hepatic glucose metabolism in an animal model. J. Hepatol. 2013, 58, 148–154. [Google Scholar] [CrossRef] [Green Version]
  23. Sun, Q.; Yue, P.; Deiuliis, J.A.; Lumeng, C.N.; Kampfrath, T.; Mikolaj, M.B.; Cai, Y.; Ostrwski, M.C.; Lu, B.; Parthasarathy, S.; et al. Ambient air pollution exaggerates adipose inflammation and insulin resistance in a mouse model of diet-induced obesity. Circulation 2009, 119, 538–546. [Google Scholar] [CrossRef] [Green Version]
  24. Park, S.K.; Wang, W. Ambient Air Pollution and Type 2 Diabetes: A Systematic Review of Epidemiologic Research. Curr. Environ. Health Rep. 2014, 1, 275–286. [Google Scholar] [CrossRef] [Green Version]
  25. Wang, B.; Xu, D.; Jing, Z.; Liu, D.; Yan, S.; Wang, Y. Effect of long-term exposure to air pollution on type 2 diabetes mellitus risk: A systemic review and meta-analysis of cohort studies. Eur. J. Endocrinol. 2014, 171, R173–R182. [Google Scholar] [CrossRef] [Green Version]
  26. Eze, I.C.; Hemkens, L.G.; Bucher, H.C.; Hoffmann, B.; Schindler, C.; Künzli, N.; Schilowski, T.; Probst-Hensch, N.M. Association between ambient air pollution and diabetes mellitus in Europe and North America: Systematic review and meta-analysis. Environ. Health Perspect. 2015, 123, 381–389. [Google Scholar] [CrossRef]
  27. He, D.; Wu, S.; Zhao, H.; Qiu, H.; Fu, Y.; Li, X.; He, Y. Association between particulate matter 2.5 and diabetes mellitus: A meta-analysis of cohort studies. J. Diabetes Investig. 2017, 8, 687–696. [Google Scholar] [CrossRef]
  28. Yang, B.Y.; Qian, Z.M.; Li, S.; Chen, G.; Bloom, M.S.; Elliott, M.; Syberg, K.W.; Heinrich, J.; Markevych, I.; Wang, S.Q.; et al. Ambient air pollution in relation to diabetes and glucose-homoeostasis markers in China: A cross-sectional study with findings from the 33 Communities Chinese Health Study. Lancet Planet Health 2018, 2, e64–e73. [Google Scholar] [CrossRef] [Green Version]
  29. Diaz, E.; Marien, K.; Manahan, L.; Fox, J. Summary of Health Research on Ultrafine Particles; Report; Washington State Department of Health: Seattle, WA, USA, December 2019. Available online: https://www.doh.wa.gov/Portals/1/Documents/4000/334-454.pdf (accessed on 13 May 2020).
  30. Hussein, T.; Puustinen, A.; Aalto, P.P.; Mäkelä, J.M.; Hämeri, K.; Kulmala, M. Urban aerosol number size distributions. Atmos. Chem. Phys. 2004, 4, 391–411. [Google Scholar] [CrossRef] [Green Version]
  31. Shi, J.P.; Evans, D.E.; Khan, A.A.; Harrison, R.M. Sources and concentration of nanoparticles (<10 nm diameter) in the urban atmosphere. Atmos. Environ. 2001, 35, 1193–1202. [Google Scholar]
  32. Baldauf, R.W.; Devlin, R.B.; Gehr, P.; Giannelli, R.; Hassett-Sipple, B.; Jung, H.; Martini, G.; McDonald, J.; Sacks, J.; Walker, K. Ultrafine particle metrics and research considerations: Review of the 2015 UFP workshop. Int. J. Environ. Res. Public Health 2016, 13, 1054. [Google Scholar] [CrossRef]
  33. Manigrasso, M.; Natale, C.; Vitali, M.; Protano, C.; Avino, P. Pedestrians in Traffic Environments: Ultrafine Particle Respiratory Doses. Int. J. Environ. Res. Public Health 2017, 14, 288. [Google Scholar] [CrossRef] [Green Version]
  34. Marcias, G.; Fostinelli, J.; Catalani, S.; Uras, M.; Sanna, A.M.; Avataneo, G.; De Palma, G.; Fabbri, D.; Paganelli, M.; Lecca, L.I.; et al. Composition of Metallic Elements and Size Distribution of Fine and Ultrafine Particles in a Steelmaking Factory. Int. J. Environ. Res. Public Health 2018, 15, 1192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kreyling, W.G.; Semmler-Behnke, M.; Möller, W. Ultrafine particle-lung interactions: Does size matter? J. Aerosol. Med. 2006, 19, 74–83. [Google Scholar] [CrossRef]
  36. Li, N.; Georas, S.; Alexis, N.; Fritz, P.; Xia, T.; Williams, M.A.; Horner, E.; Nel, A. A work group report on ultrafine particles (AAAAI) why ambient ultrafine and engineered nanoparticles should receive special attention for possible adverse health outcomes in humans. J. Allergy Clin. Immunol. 2016, 138, 386–396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Li, N.; Sioutas, C.; Cho, A.; Schmitz, D.; Misra, C.; Sempf, J.; Wang, M.; Oberley, T.; Froines, J.; Nel, A. Ultrafine particulate pollutants induce oxidative stress and mitochondrial damage. Environ. Health Perspect. 2003, 111, 455–460. [Google Scholar] [CrossRef] [PubMed]
  38. Bhargava, A.; Tamrakar, S.; Aglawe, A.; Lad, H.; Srivastava, R.K.; Mishra, D.K.; Tiwri, R.; Chaudhury, K.; Goryacheva, I.Y.; Mishra, P.K. Ultrafine particulate matter impairs mitochondrial redox homeostasis and activates phosphatidylinositol 3-kinase mediated DNA damage responses in lymphocytes. Environ. Pollut. 2018, 234, 406–419. [Google Scholar] [CrossRef] [PubMed]
  39. Cheng, H.; Saffari, A.; Sioutas, C.; Forman, H.J.; Morgan, T.E.; Finch, C.E. Nano-scale particulate matter from urban traffic rapidly induces oxidative stress and inflammation in olfactory epithelium with concomitant effects on brain. Environ. Health Perspect. 2016, 124, 1537–1546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Tyler, C.R.; Zychoowski, K.E.; Snchez, B.N.; Rivero, V.; Lucas, S.; Herbert, G.; Liu, J.; Irshad, H.; McDonald, J.D.; Bleske, B.E. Surface area-dependence of gas-particle interactions influences pulmonary and neuroinflammatory outcomes. Part. Fibre Toxicol. 2016, 13, 64. [Google Scholar] [CrossRef] [Green Version]
  41. Zhang, H.; Liu, H.; Davies, K.J.A.; Siutas, C.; Finch, C.E.; Morgan, T.E.; Forman, H.J. Nrf2-regulated phase II enzymes are induced by chronic ambient nanoparticle exposure in young mice with age-related impairments. Free Radic. Biol. Med. 2012, 52, 2038–2046. [Google Scholar] [CrossRef] [Green Version]
  42. Seagrave, J.; Campen, M.J.; McDonald, J.D.; Mauderly, J.L.; Rohr, A.C. Oxidative stress, inflammation, and pulmonary function assessment in rats exposed to laboratory-generated pollutant mixtures. J. Toxicol. Environ. Health 2008, 71, 1352–1362. [Google Scholar] [CrossRef]
  43. Weissenberg, A.; Sydlik, U.; Peuschel, H.; Schroeder, P.; Schneider, M.; Schins, R.P.F.; Abel, J.; Unfried, K. Reactive oxygen species as mediators of membrane-dependent signaling induced by ultrafine particles. Free Radic. Biol Med. 2010, 49, 597–605. [Google Scholar] [CrossRef]
  44. Morgan, T.E.; Davis, D.A.; Iwata, N.; Tanner, J.A.; Snyder, D.; Ning, Z.; Kam, W.; Hsu, Y.-T.; Winkler, J.W.; Chen, J.-C.; et al. Glutamatergic neurons in rodent models respond to nanoscale particulate urban air pollutants in vivo and in vitro. Environ. Health Perspect. 2011, 119, 1003–1009. [Google Scholar] [CrossRef]
  45. Kleinman, M.T.; Araujo, J.; Nel, A.; Siutas, C.; Campbell, A.; Cong, P.Q.; Li, H.; Bondy, S.C. Inhaled ultrafine particulate matter affects CNS inflammatory processes and may act via MAP kinase signaling pathways. Toxicol. Lett. 2008, 178, 127–130. [Google Scholar] [CrossRef] [Green Version]
  46. Allen, J.L.; Liu, X.; Weston, D.; Prince, L.; Oberdörster, G.; Finkelstein, J.N.; Johnston, C.J.; Cory-Slechta, D.A. Developmental exposure to concentrated ambient ultrafine particulate matter air pollution in mice results in persistent and sex-dependent behavioral neurotoxicity and glial activation. Toxicol. Sci. 2014, 140, 160–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Woodward, N.C.; Pakbin, P.; Saffari, A.; Shirmohammadi, F.; Haghani, A.; Sioutas, C.; CCacciottolo, M.; Morgan, T.E.; Caleb, E.F. Traffic-related air pollution impact on mouse brain accelerates myelin and neuritic aging changes with specificity for CA1 neurons. Neurobiol. Aging 2017, 53, 48–58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Klocke, C.; Allen, J.L.; Sobolewski, M.; Mayer-Pröschel, M.; Blum, J.L.; Lauterstein, D.; Zelikoff, J.T.; Cory-Slechta, D.A. Neuropathological Consequences of Gestational Exposure to Concentrated Ambient Fine and Ultrafine Particles in the Mouse. Toxicol. Sci. 2017, 156, 492–508. [Google Scholar] [CrossRef] [PubMed]
  49. Allen, J.L.; Conrad, K.; Oberdörster, G.; Johnston, C.J.; Sleezer, B.; Cory-Slechta, D.A. Developmental exposure to concentrated ambient particles and preference for immediate reward in mice. Environ. Health Perspect. 2013, 121, 32–38. [Google Scholar] [CrossRef] [PubMed]
  50. Davis, D.A.; Bortolato, M.; Godar, S.C.; Sander, T.K.; Iwata, N.; Pakbin, P.; Shih, J.C.; Berhane, K.; McConnell, R.; Siutas, C.; et al. Prenatal exposure to urban air nanoparticles in mice causes altered neuronal differentiation and depression-like responses. PLoS ONE 2013, 8, e64128. [Google Scholar] [CrossRef] [PubMed]
  51. Fonken, L.K.; Xu, X.; Weil, Z.M.; Chen, G.; Sun, Q.; Rajagopalan, S.; Nelson, R.J. Air pollution impairs cognition, provokes depressive-like behaviors and alters hippocampal cytokine expression and morphology. Mol. Psychiatry 2011, 16, 973, 987–995. [Google Scholar] [CrossRef] [Green Version]
  52. Suzuki, T.; Oshio, S.; Iwata, M.; Saburi, H.; Odagiri, T.; Udagawa, T.; Sugawara, I.; Umezawa, M.; Takeda, K. In utero exposure to a low concentration of diesel exhaust affects spontaneous locomotor activity and monoaminergic system in male mice. Part. Fibre Toxicol. 2010, 7. [Google Scholar] [CrossRef] [Green Version]
  53. Li, N.; Harkema, J.R.; Lewandowski, R.P.; Wang, M.; Bramble, L.A.; Gookin, G.R.; Ning, Z.; Kleinman, M.T.; Sioutas, C.; Nel, A.E. Ambient ultrafine particles provide a strong adjuvant effect in the secondary immune response: Implication for traffic-related asthma flares. Am. J. Physiol. Lung Cell Mol. Physiol. 2010, 299, 374–383. [Google Scholar] [CrossRef]
  54. Kleinman, M.; Hamade, A.; Meacher, D.; Oldham, M.; Sioutas, C.; Chakrabarti, B.; Stram, D.; Froines, J.R.; Cho, A.K. Inhalation of concentrated ambient particulate matter near a heavily trafficked road stimulates antigen-induced airway responses in mice. J. Air Waste Manag. Assoc. 2005, 55, 1277–1288. [Google Scholar] [CrossRef] [Green Version]
  55. Saleh, Y.; Antherieu, S.; Dusautoir, R.; Alleman, L.Y.; Sotty, J.; De Sousa, C.; Platel, A.; Perdrix, E.; Riffault, V.; Fronval, I.; et al. Exposure to Atmospheric Ultrafine Particles Induces Severe Lung Inflammatory Response and Tissue Remodeling in Mice. Int. J. Environ. Res. Public Health 2019, 16, 1210. [Google Scholar] [CrossRef] [Green Version]
  56. Li, C.; Li, X.; Jigami, J.; Hasegawa, C.; Suzuki, A.K.; Zhang, Y.; Fujitani, Y.; Nagaoka, K.; Watanabe, G.; Taya, K. Effect of nanoparticle-rich diesel exhaust on testosterone biosynthesis in adult male mice. Inhal. Toxicol. 2012, 24, 599–608. [Google Scholar] [CrossRef] [PubMed]
  57. Laurent, O.; Hu, J.; Li, L.; Cockburn, M.; Escobedo, L.; Kleeman, M.J.; Wu, J. Sources and contents of air pollution affecting term low birth weight in Los Angeles County, California, 2001–2008. Environ. Res. 2014, 134, 488–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Laurent, O.; Hu, J.; Li, L.; Kleeman, M.J.; Bartell, S.M.; Cockburn, M.; Ecsobedo, L.; Wu, J. A statewide nested case-control study of preterm birth and air pollution by source and composition: California, 2001–2008. Environ. Health Perspect. 2016, 124, 1479–1486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Chang, Y.C.; Cole, T.B.; Costa, L.G. Prenatal and early-life diesel exhaust exposure causes autism-like behavioral changes in mice. Part. Fibre Toxicol. 2018, 15, 18. [Google Scholar] [CrossRef]
  60. Bourdon, M.; Torres-Rovira, L.; Monniaux, D.; Faure, C.; Levy, R.; Tarrade, A.; Rousseau-Ralliard, D.; Chavatte-Palmer, P.; Jolivet, G. Impact of a gestational exposure to diesel exhaust on offspring gonadal development: Experimental study in the rabbit. J. Dev. Orig. Health Dis. 2018, 9, 519–529. [Google Scholar] [CrossRef] [PubMed]
  61. Integrated Science Assessment (ISA) For Particulate Matter; Final Report; U.S. Environmental Protection Agency: Washington, DC, USA, December 2009. Available online: https://cfpub.epa.gov/ncea/risk/recordisplay.cfm?deid=216546 (accessed on 14 May 2020).
  62. Ohlwein, S.; Kappeler, R.; Kutlar Joss, M.; Künzli, N.; Hoffmann, B. Health effects of ultrafine particles: A systematic literature review update of epidemiological evidence. Int. J. Public Health 2019, 64, 547–559. [Google Scholar] [CrossRef]
  63. HEI Review Panel on Ultrafine Particles. Understanding the Health Effects of Ambient Ultrafine Particles. HEI Perspectives 3; Health Effects Institute: Boston, MA, USA, 2013; Available online: https://www.healtheffects.org/publication/understanding-health-effects-ambient-ultrafine-particles (accessed on 14 May 2020).
  64. Population by Age and Sex, Municipalities, Slovenia, Half-Yearly. The Statistical Office of the Republic of Slovenia (SURS). Available online: https://pxweb.stat.si/SiStatDb/pxweb/en/10_Dem_soc/10_Dem_soc__05_prebivalstvo__10_stevilo_preb__20_05C40_prebivalstvo_obcine/05C4002S.px/. (accessed on 15 May 2020).
  65. SOER. Air Pollution. Slovenian Environment Agency. Available online: https://www.arso.gov.si/en/soer/air_pollution.html (accessed on 15 May 2020).
  66. WHO. International Statistical Classification of Diseases and Related Health Problems 10th Revision (ICD-10). Available online: https://icd.who.int/browse10/2010/en (accessed on 15 May 2020).
  67. UFIREG. Data Collection and Methods; Report; UFIREG: Dresden, Germany, 2014; Available online: http://www.umwelt.sachsen.de/umwelt/download/UFIREG_Handbook_2014_final.pdf (accessed on 15 May 2020).
  68. Slovenian Environment Agency. Available online: https://www.arso.gov.si (accessed on 14 November 2019).
  69. Epidemiološko Spremljanje Nalezljivih Bolezni (Epidemiological Surveillance of Communicable Diseases in Slovenia). Nacionalni Inštitut za Javno Zdravje (National Institute of Public Health of the Republic of Slovenia). Available online: https://www.nijz.si/sl/epidemiolosko-spremljanje-nalezljivih-bolezni-letna-in-cetrtletna-porocila (accessed on 15 May 2020).
  70. National Laboratory of Health, Environment and Food (Nacionalni Laboratorij za Zdravje, Okolje in Hrano). Zbirka Dve-Urnih Koncentracij Cvetnega Prahu (Two-Hour Collection of Pollen Concentrations). Available online: https://www.nlzoh.si (accessed on 1 August 2019).
  71. Morgenstern, H.; Thomas, D. Principles of study design in environmental epidemiology. Environ. Health Perspect. 1993, 101, 23–38. [Google Scholar]
  72. De Souza Tadano, Y.; Ugaya, C.M.L.; Franco, A.T. Methodology to assess air pollution impact on human health using the generalized linear model with Poisson Regression. In Air Pollution-Monitoring, Modelling and Health; Khare, M., Ed.; In Tech: Rijeka, Croatia, 2012. [Google Scholar]
  73. Lanzinger, S.; Schneider, A.; Breitner, S.; Stafoggia, M.; Erzen, I.; Dostal, M.; Pastorkova, A.; Bastian, S.; Cyrys, J.; Zscheppang, A.; et al. Ultrafine and Fine Particles and hospital admission in Central Europe. Results from the UFIREG study. Am. J. Respir. Crit Med. 2016, 194, 1233–1241. [Google Scholar] [CrossRef] [Green Version]
  74. Bai, L.; Chen, H.; Hatzopoulou, M.; Jerrett, M.; Kwong, J.C.; Burnett, R.T.; van Donkelaar, A.; Copes, R.; Martin, R.V.; Van Ryswyk, K.; et al. Exposure to Ambient Ultrafine Particles and Nitrogen Dioxide and Incident Hypertension and Diabetes. Epidemiology 2018, 29, 323–332. [Google Scholar] [CrossRef]
  75. Karottki, D.G.; Bekö, G.; Clausen, G.; Madsen, A.M.; Andersen, Z.J.; Massling, A.; Ketzel, M.; Ellermann, T.; Lund, R.; Sigsgaard, T.; et al. Cardiovascular and lung function in relation to outdoor and indoor exposure to fine and ultrafine particulate matter in middle-aged subjects. Environ. Int. 2014, 73, 372–381. [Google Scholar] [CrossRef] [Green Version]
  76. Li, Y.; Lane, K.J.; Corlin, L.; Patton, A.P.; Durant, J.L.; Thanikachalam, M.; Woodin, M.; Wang, M.; Brugge, D. Association of Long-Term Near-Highway Exposure to Ultrafine Particles with Cardiovascular Diseases, Diabetes and Hypertension. Int. J. Environ. Res. Public Health 2017, 14, 461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Johansson, C.; Norman, M.; Gidhagen, L. Spatial & temporal variations of PM10 and particle number concentrations in urban air. Environ. Monit. Assess. 2007, 127, 477–487. [Google Scholar]
  78. Cyrys, J.; Pitz, M.; Heinrich, J.; Wichmann, H.-E.; Peters, A. Spatial and temporal variation of particle number concentration in Augsburg, Germany. Sci. Total Environ. 2008, 401, 168–175. [Google Scholar] [CrossRef] [PubMed]
  79. Birmili, W.; Tomsche, L.; Sonntag, A.; Opelt, C.; Weinhold, K.; Nordmann, S.; Schmidt, W. Variability of aerosol particles in the urban atmosphere of Dresden (Germany): Effects of spatial scale and particle size. Meteorol. Z. 2013, 22, 195–211. [Google Scholar] [CrossRef]
  80. Weinmayr, G.; Hennig, F.; Fuks, K.; Nonnemacher, M.; Jacobs, H.; Möhlenkamp, S.; Erbel, R.; Jöckel, K.-H.; Hoffman, B.; Moebus, S. Long-term exposure to fine particulate matter and incidence of type 2 diabetes mellitus in a cohort study: Effects of total and traffic-specific air pollution. Environ. Health 2015, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Strak, M.; Janssen, N.; Beelen, R.; Schmitz, O.; Vaartjes, I.; Karssenberg, D.; van den Brink, C.; Bots, M.L.; Dijst, M.; Brunekreef, B.; et al. Long-term exposure to particulate matter, NO2 and the oxidative potential of particulates and diabetes prevalence in a large national health survey. Environ. Int. 2017, 108, 228–236. [Google Scholar] [CrossRef]
  82. Sun, S.; Qiu, H.; Ho, K.-F.; Tian, L. Chemical components of respirable particulate matter associated with emergency hospital admissions for type 2 diabetes mellitus in Hong Kong. Environ. Int. 2016, 97, 93–99. [Google Scholar] [CrossRef] [Green Version]
  83. Zanobetti, A.; Franklin, M.; Koutrakis, P.; Schwartz, J. Fine particulate air pollution and its components in association with cause-specific emergency admissions. Environ. Health 2009, 8, 58. [Google Scholar] [CrossRef] [Green Version]
  84. Kukec, A.; Zaletel-Kragelj, L.; Bizjak, M.; Fink, R.; Jereb, G.; Košnik, M.; Poljšak, B.; Zadnik, V.; Farkaš-Lainščak, J. A Study of Linkage of Environmental and Health Data in Zasavje as a Model Study to Support the Development and Implementation of Cross-Sectoral Policies on the Environment and Health; Final Report; Chair of Public Health: Ljubljana, Slovenia, 2012. (In Slovene) [Google Scholar]
  85. Bhaskaran, K.; Gasparrini, A.; Shakoor, H.; Smeeth, L.; Armstrong, B. Time series regression studies in environmental epidemiology. Int. J. Epidemiol. 2013, 42, 1187–1195. [Google Scholar] [CrossRef]
  86. Šimac, N. Air Pollution by Ozone in Nova Gorica Region—Assessment of the Impact on Human Health; Specialistic Study Thesis; Regional Institute of Public Health: Nova Gorica, Slovenia, 2008. (In Slovene) [Google Scholar]
  87. Šimac, N.; Hladnik, M.; Zaletel-Kragelj, L. The impact of temperature on tropospheric ozone in the Nova Gorica region. Zdr. Var. 2011, 50, 121–130. [Google Scholar] [CrossRef] [Green Version]
  88. Rems-Novak, M.M. Effects of air pollution with ozone on primary health care consultations for respiratory tract diseases in children in Koper Municipality. Master’s Thesis, University of Nova Gorica, Graduate School of Environmental Sciences, Nova Gorica, Slovenia, 2013. [Google Scholar]
  89. Kukec, A.; Eržen, I.; Farkaš, J.; Zaletel Kragelj, L. Impact of air pollution with PM10 on primary health care consultations for respiratory diseases in children in Zasavje, Slovenia: A time-trend study. Zdr. Var. 2014, 53, 55–68. [Google Scholar] [CrossRef]
  90. Paulin, S.; Nadrag, P.; Kelšin, N.; Truden-Dobrin, P.; Korošec, A.; Simončič, M.; Zaletel, M.; Zupanc, I.; Eržen, I.; Zaletal, J.; et al. Ožji Nabor Kazalnikov za Spremljanje Obvladovanja Sladkorne Bolezni v Sloveniji; Nacionalni Inštitut za Javno Zdravje: Ljubljana, Slovenia, 2019; Available online: https://www.nijz.si/ (accessed on 24 June 2020).
  91. IDF Diabetes Atlas. Ninth Edition. 2019. Available online: https://www.diabetesatlas.org (accessed on 25 June 2020).
Figure 1. A map of Slovenia surrounded by bordering countries as well as the location of the capital city Ljubljana with the borders of the municipality, hydrography, and roads.
Figure 1. A map of Slovenia surrounded by bordering countries as well as the location of the capital city Ljubljana with the borders of the municipality, hydrography, and roads.
Ijerph 17 04970 g001
Figure 2. Temporal variability of the daily number of consultations due to diabetes mellitus: (a) all age groups; (b) children; (c) the elderly in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
Figure 2. Temporal variability of the daily number of consultations due to diabetes mellitus: (a) all age groups; (b) children; (c) the elderly in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
Ijerph 17 04970 g002
Figure 3. Temporal variability of: (a) average 24-h number concentration of UFP (number/cm3); (b) average 24-h concentration of PM10 and PM2.5 (µg/m3); (c) daily average pollen concentration (number of particles/m3) in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
Figure 3. Temporal variability of: (a) average 24-h number concentration of UFP (number/cm3); (b) average 24-h concentration of PM10 and PM2.5 (µg/m3); (c) daily average pollen concentration (number of particles/m3) in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
Ijerph 17 04970 g003
Table 1. Descriptive statistics for the daily number of consultations due to diabetes mellitus and ultrafine, fine, and coarse particulate matter (PM) and meteorological parameters and pollen concentration in the municipality of Ljubljana (MOL), Slovenia between 1 January 2013 and 31 December 2017.
Table 1. Descriptive statistics for the daily number of consultations due to diabetes mellitus and ultrafine, fine, and coarse particulate matter (PM) and meteorological parameters and pollen concentration in the municipality of Ljubljana (MOL), Slovenia between 1 January 2013 and 31 December 2017.
nMeanSDMin.Q1MedianQ3Max.
Daily number of consultations due to diabetes mellitus: all age groups, children, elderly
All age groups182611.219.350.001.0011.0018.0042.0
Children18260.040.210.000.000.000.002.00
Elderly18268.177.190.001.008.0013.0033.00
Average 24-h number concentration of UFP (number/cm3)
UFP0.01–0.0216551230.53478.2816.50915.301172.301482.504370.40
UFP0.02–0.0316551086.93434.3088.10809.501031.301308.204769.40
UFP0.03–0.0516551517.34605.96247.501113.001420.301819.204446.20
UFP0.05–0.0716551039.73467.17198.10732.200932.401224.203446.70
UFP0.07–0.1016551070.34578.18181.20696.400924.301251.404207.90
Average 24-h concentration of PM10 and PM2.5 (µg/m3)
PM10176525.4616.622.0015.0020.0030.00126.00
PM2.5174220.6715.253.0011.0016.0024.00114.00
Average 24-h of temperature (°C) and relative humidity (%)
Temperature182612.068.107−7.905.8012.5018.3230.00
Relative humidity182675.9913.6631.0065.0076.0088.00100.00
Daily average pollen concentration (number of particles/m3)
Pollen concentration18268.9927.360.000.001.007.00323.00
Legend: n—total number of data; SD—standard deviation; Q1—the first quartile; Q3—the third quartile; UFP—ultrafine particles; PM10—particulate matter 10 micrometers or less in diameter; PM2.5—particulate matter 2.5 micrometers or less in diameter.
Table 2. Association between the daily number of consultations due to diabetes mellitus by all age groups and the average 24-h number concentration of UFP (number/cm3); average 24-h concentration of PM10 and PM2.5 (µg/m3) in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
Table 2. Association between the daily number of consultations due to diabetes mellitus by all age groups and the average 24-h number concentration of UFP (number/cm3); average 24-h concentration of PM10 and PM2.5 (µg/m3) in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
PollutantsModel 1 aModel 2 b
IRR95% CIpIRR95% CIp
UFP0.01–0.02_lag01.0011.0001.003<0.0011.0011.0001.003<0.001
UFP0.02–0.03_lag01.0011.0001.003<0.0011.0011.0001.003<0.001
UFP0.03–0.05_lag01.0001.0001.000<0.0011.0001.0001.0000.133
UFP0.05–0.07_lag01.0001.0001.000<0.0011.0001.0001.0000.590
UFP0.07–0.10_lag01.0001.0001.000<0.0011.0001.0001.0000.689
PM2.5_lag21.0021.0021.003<0.0011.0011.0001.0020.130
PM10_lag01.0021.0011.003<0.0011.0011.0001.0020.606
Legend: a non-adjusted model; b adjusted model for season of the year (spring, summer, autumn, winter), working day, influenza season, pollen concentration, air temperature, and relative humidity. IRR—incident rate ratio; CI—confidence interval; pp-value of 0.05 or less was considered as statistically significant; UFP—ultrafine particles; PM10—particulate matter 10 micrometers or less in diameter; PM2.5—particulate matter 2.5 micrometers or less in diameter.
Table 3. Association between the daily number of consultations due to diabetes mellitus by the population of children and the average 24-h number concentration of UFP (number/cm3); average 24-h concentration of PM10 and PM2.5 (µg/m3) in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
Table 3. Association between the daily number of consultations due to diabetes mellitus by the population of children and the average 24-h number concentration of UFP (number/cm3); average 24-h concentration of PM10 and PM2.5 (µg/m3) in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
PollutantsModel 1 aModel 2 b
IRR95% CIpIRR95% CIp
UFP0.01–0.02_lag01.0001.0001.0010.1641.0001.0001.0010.663
UFP0.02–0.03_lag01.0001.0001.0010.0451.0001.0001.0010.377
UFP0.03–0.05_lag01.0001.0001.0010.0231.0001.0001.0010.237
UFP0.05–0.07_lag01.0001.0001.0010.1541.0001.0001.0010.534
UFP0.07–0.10_lag01.0001.0001.0000.6761.0000.9991.0010.873
PM2.5_lag21.0040.9881.0180.5711.0040.9851.0210.668
PM10_lag01.0000.9851.0130.9960.9970.9781.0130.699
Legend: a non-adjusted model; b adjusted model for season of the year (spring, summer, autumn, winter), working day, influenza season, pollen concentration, air temperature, and relative humidity. IRR—incident rate ratio; CI—confidence interval; pp-value of 0.05 or less was considered as statistically significant; UFP—ultrafine particles; PM10—particulate matter 10 micrometers or less in diameter; PM2.5—particulate matter 2.5 micrometers or less in diameter.
Table 4. Association between the daily number of consultations due to diabetes mellitus by the elderly population and the average 24-h number concentration of UFP (number/cm3); average 24-h concentration of PM10 and PM2.5 (µg/m3) in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
Table 4. Association between the daily number of consultations due to diabetes mellitus by the elderly population and the average 24-h number concentration of UFP (number/cm3); average 24-h concentration of PM10 and PM2.5 (µg/m3) in the MOL, Slovenia between 1 January 2013 and 31 December 2017.
PollutantsModel 1 aModel 2 b
IRR95% CIpIRR95% CIp
UFP0.01–0.02_lag01.0011.0001.003<0.0011.0011.0001.003<0.001
UFP0.02–0.03_lag01.0011.0001.003<0.0011.0011.0001.003<0.001
UFP0.03–0.05_lag01.0001.0001.000<0.0011.0001.0001.0000.065
UFP0.05–0.07_lag01.0001.0001.000<0.0011.0001.0001.0000.419
UFP0.07–0.10_lag01.0001.0001.000<0.0011.0001.0001.0000.560
PM2.5_lag21.0031.0021.004<0.0011.0011.0001.0030.057
PM10_lag01.0021.0011.003<0.0011.0000.9991.0010.795
Legend: a non-adjusted model; b adjusted model for season of the year (spring, summer, autumn, winter), working day, influenza season, pollen concentration, air temperature, and relative humidity. IRR—incident rate ratio; CI—confidence interval; pp-value of 0.05 or less was considered as statistically significant; UFP—ultrafine particles; PM10—particulate matter 10 micrometers or less in diameter; PM2.5—particulate matter 2.5 micrometers or less in diameter.

Share and Cite

MDPI and ACS Style

Viher Hrženjak, V.; Kukec, A.; Eržen, I.; Stanimirović, D. Effects of Ultrafine Particles in Ambient Air on Primary Health Care Consultations for Diabetes in Children and Elderly Population in Ljubljana, Slovenia: A 5-Year Time-Trend Study. Int. J. Environ. Res. Public Health 2020, 17, 4970. https://doi.org/10.3390/ijerph17144970

AMA Style

Viher Hrženjak V, Kukec A, Eržen I, Stanimirović D. Effects of Ultrafine Particles in Ambient Air on Primary Health Care Consultations for Diabetes in Children and Elderly Population in Ljubljana, Slovenia: A 5-Year Time-Trend Study. International Journal of Environmental Research and Public Health. 2020; 17(14):4970. https://doi.org/10.3390/ijerph17144970

Chicago/Turabian Style

Viher Hrženjak, Vesna, Andreja Kukec, Ivan Eržen, and Dalibor Stanimirović. 2020. "Effects of Ultrafine Particles in Ambient Air on Primary Health Care Consultations for Diabetes in Children and Elderly Population in Ljubljana, Slovenia: A 5-Year Time-Trend Study" International Journal of Environmental Research and Public Health 17, no. 14: 4970. https://doi.org/10.3390/ijerph17144970

APA Style

Viher Hrženjak, V., Kukec, A., Eržen, I., & Stanimirović, D. (2020). Effects of Ultrafine Particles in Ambient Air on Primary Health Care Consultations for Diabetes in Children and Elderly Population in Ljubljana, Slovenia: A 5-Year Time-Trend Study. International Journal of Environmental Research and Public Health, 17(14), 4970. https://doi.org/10.3390/ijerph17144970

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop