Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Traffic-related air pollution and spectacles use in schoolchildren

  • Payam Dadvand ,

    payam.dadvand@isglobal.org

    Affiliations ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain, Pompeu Fabra University, Barcelona, Spain, Ciber on Epidemiology and Public Health (CIBERESP), Madrid, Spain

  • Mark J. Nieuwenhuijsen,

    Affiliations ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain, Pompeu Fabra University, Barcelona, Spain, Ciber on Epidemiology and Public Health (CIBERESP), Madrid, Spain

  • Xavier Basagaña,

    Affiliations ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain, Pompeu Fabra University, Barcelona, Spain, Ciber on Epidemiology and Public Health (CIBERESP), Madrid, Spain

  • Mar Alvarez-Pedrerol,

    Affiliations ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain, Pompeu Fabra University, Barcelona, Spain, Ciber on Epidemiology and Public Health (CIBERESP), Madrid, Spain

  • Albert Dalmau-Bueno,

    Affiliations ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain, Pompeu Fabra University, Barcelona, Spain, Ciber on Epidemiology and Public Health (CIBERESP), Madrid, Spain

  • Marta Cirach,

    Affiliations ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain, Pompeu Fabra University, Barcelona, Spain, Ciber on Epidemiology and Public Health (CIBERESP), Madrid, Spain

  • Ioar Rivas,

    Affiliations ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain, Pompeu Fabra University, Barcelona, Spain, Ciber on Epidemiology and Public Health (CIBERESP), Madrid, Spain, Institute of Environmental Assessment and Water Research (IDAEA), Spanish National Research Council (CSIC), Barcelona, Spain

  • Bert Brunekreef,

    Affiliations Institute for Risk Assessment Sciences, Division Environmental Epidemiology, Utrecht University, Utrecht, The Netherlands, Julius Center for Health Sciences and Primary Care, University Medical Center Utrecht, Utrecht, The Netherlands

  • Xavier Querol,

    Affiliation Institute of Environmental Assessment and Water Research (IDAEA), Spanish National Research Council (CSIC), Barcelona, Spain

  • Ian G. Morgan,

    Affiliations ARC Centre of Excellence in Vision Science, Research School of Biology, College of Medicine, Biology and Environment, Australian National University, Canberra, Australia, Department of Preventive Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China

  • Jordi Sunyer

    Affiliations ISGlobal, Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain, Pompeu Fabra University, Barcelona, Spain, Ciber on Epidemiology and Public Health (CIBERESP), Madrid, Spain, IMIM-Parc Salut Mar, Spain

Abstract

Purpose

To investigate the association between exposure to traffic-related air pollution and use of spectacles (as a surrogate measure for myopia) in schoolchildren.

Methods

We analyzed the impact of exposure to NO2 and PM2.5 light absorbance at home (predicted by land-use regression models) and exposure to NO2 and black carbon (BC) at school (measured by monitoring campaigns) on the use of spectacles in a cohort of 2727 schoolchildren (7–10 years old) in Barcelona (2012–2015). We conducted cross-sectional analyses based on lifelong exposure to air pollution and prevalent cases of spectacles at baseline data collection campaign as well as longitudinal analyses based on incident cases of spectacles use and exposure to air pollution during the three-year period between the baseline and last data collection campaigns. Logistic regression models were developed to quantify the association between spectacles use and each of air pollutants adjusted for relevant covariates.

Results

An interquartile range increase in exposure to NO2 and PM2.5 absorbance at home was respectively associated with odds ratios (95% confidence intervals (CIs)) for spectacles use of 1.16 (1.03, 1.29) and 1.13 (0.99, 1.28) in cross-sectional analyses and 1.15 (1.00, 1.33) and 1.23 (1.03, 1.46) in longitudinal analyses. Similarly, odds ratio (95% CIs) of spectacles use associated with an interquartile range increase in exposures to NO2 and black carbon at school was respectively 1.32 (1.09, 1.59) and 1.13 (0.97, 1.32) in cross-sectional analyses and 1.12 (0.84, 1.50) and 1.27 (1.03, 1.56) in longitudinal analyses. These findings were robust to a range of sensitivity analyses that we conducted.

Conclusion

We observed increased risk of spectacles use associated with exposure to traffic-related air pollution. These findings require further confirmation by future studies applying more refined outcome measures such as quantified visual acuity and separating different types of refractive errors.

Introduction

Myopia is the most common refractive error of vision, currently affecting about one-fifth of the world’s population (~1.5 billion people) [14]. Once considered a purely genetic condition, it is now increasingly recognized as having a multifactorial etiology, with both genetic and environmental factors involved [1, 57]. During the past few decades, there has been a notable increase in the global prevalence of myopia, representing an alarming epidemic worldwide [14]. Although the reason(s) for this increasing trend are yet to be established, such a rapid increase can be suggestive for a more important contribution of non-genetic and environmental factors in the causation of refractive errors [1, 5, 8].

The increase in the global prevalence of myopia has coincided with the rapid and ongoing increase in the population residing in urban areas where the prevalence of myopia is consistently reported to be higher than rural areas [1, 4, 911]. The higher prevalence of myopia in urban areas could suggest that urban lifestyle such as more near-work (i.e. tasks such as reading book and working with computer that need sustained gaze on a close object) or less time spent outdoor and/or urban-related environmental factors contribute to the pathogenesis of these conditions. Air pollution is the main environmental stressor in urban areas, and is responsible for most of the global burden of disease due to environmental causes [12]. Exposure to traffic-related air pollution is associated with a wide range of adverse health outcomes, with the lungs being one of the most commonly affected organs, mainly because of their constant direct exposure to air pollutants. Similarly, the eyes are directly exposed to air pollution, making them a prime target organ for the adverse effects of such an exposure. In addition to the short-term effects of air pollution on the eye, such as irritation of the ocular surface and its accompanying symptoms and complaints, chronic exposure to air pollution has been associated with long-lasting ocular conditions such as dry eye disease [13] and cataract [14]. Although air pollution could induce myopia through systemic inflammation and oxidative stress (as discussed later in the text), to date no studies have reported on the potential effect of air pollution on the development of myopia.

The aim of this analysis was to investigate the association between exposure to traffic-related air pollution and use of spectacles in schoolchildren. We considered use of spectacles as a surrogate for refractive errors of vision and specifically myopia because when low visual acuity increases during childhood, this is particularly likely to be associated with the onset of myopia [1517]. Our choice of primary schoolchildren to test our hypothesized association was in line with most studies (e.g. [1823]) of the environmental determinants of refractive errors that have focused on early years of primary school as a suitable window of exposure because it is a period when environmental factors have greatest opportunity to affect rapidly changing eyes.

Materials and methods

Study setting and participants

We undertook this study in Barcelona, Spain, a port city situated on the Northeastern part of the Iberian Peninsula. Air pollution concentrations in Barcelona are among the highest in Europe, partly attributed to high traffic density and large proportion (~50%) of diesel-powered vehicles, relatively low precipitation, high population density (about 16,000/km2), and an urban landscape characterized by 5–6 story buildings and narrow streets, which reduces the dispersion of pollutants [24, 25].

This study was based on data collected by the BRain dEvelopment and Air polluTion ultrafine particles in scHool childrEn (BREATHE) project, which aimed to evaluate the impact of air pollution exposure on neurobehavioral development in primary schoolchildren. As described in detail previously [26, 27], of the 416 schools in Barcelona, 40 schools were initially selected to obtain maximum contrast in traffic-related air pollution levels (i.e. NO2) of which 39 accepted to participate and were included in the study. Participating schools were similar to the remaining schools in Barcelona in terms of the neighborhood socioeconomic vulnerability index (0.46 versus 0.50, Kruskal—Wallis test p = 0.57) and NO2 levels (51.5 versus 50.9 μg/m3, Kruskal—Wallis test p = 0.72).

We invited all schoolchildren (n = 5,019) without special needs in the 2nd to 4th grades (7–10 years) of participating schools to participate through letters and/or presentations in schools for parents, of which 2,897 (58%) agreed to participate in BREATHE. All participants had been in the same school for more than six months (and 98% more than one year) before the beginning of the study. All parents or guardians signed the informed consent and the BREATHE project was approved (No. 2010/41221/I) by the Clinical Research Ethical Committee of the Parc de Salut MAR, Barcelona, Spain.

Outcome and covariate data

We considered the use of spectacles as a surrogate for myopia and applied it as a binary (yes/no) outcome variable. Data on the use of spectacles reported by parents were collected twice: once in the first (i.e. baseline) data collection campaign during 2012 and once during the last data collection campaign during 2015. Sociodemographic data including child’s sex and age and parental ethnicity and indicators of socioeconomic status such as educational achievement and employment status together with data on pregnancy period and childhood were obtained from parents through questionnaires.

Air pollution exposure

We assessed exposure to NO2, and particulate matter with aerodynamic diameter ≤ 2.5 μm (PM2.5) light absorption (hereafter referred to as PM2.5 absorbance, a proxy for black carbon (BC)) at residential addresses and to NO2 and BC at schools. These pollutants have been used extensively as markers of air pollution generated by traffic.

Residential air pollution levels.

We utilized an established spatiotemporal exposure assessment framework based on temporally-adjusted spatial estimates of air pollutant levels by land use regression (LUR) models developed as part of the European Study of Cohorts for Air Pollution Effects (ESCAPE) [28, 29]. These models could predict 72% and 83% of variation in annual (2009) levels of NO2 and PM2.5 absorbance, respectively, across the Barcelona [30]. By temporally adjusting (ratio method) the LUR spatial estimates, we were able to predict the ambient pollutant levels at the geocoded home address of each study participant for the periods between 1) her/his birth and baseline data collection campaign (hereafter referred to as lifelong exposure) and 2) between the baseline and last data collection campaigns (hereafter referred to as prospective exposure). Further details of our applied spatiotemporal exposure assessment framework and LUR models have been reported before [2830]. As part of BREATHE questionnaire, the participants were asked to report their current address together with their previous addresses and the period they spent in each address.

School air pollution levels.

As described in detail previously [26, 31], air pollution levels at each BREATHE school were measured twice during one-week campaigns separated by six months, once during the warm and once during the cold seasons of the year 2012. Air samples were collected in a classroom (i.e. indoor) at a height between 0.7 and 1.5 m above floor level, which is at the eye level of pupils aged 7–9 years and is also the height at which they would usually inhale. Weekly averaged NO2 concentrations were measured by Gradko Environmental passive dosimeters. BC concentrations were measured using the MicroAeth AE51 (AethLabs). Considering that the air pollution sampling in different schools were conducted during different weeks in each campaign period, we deseasonalized the monitored air pollution levels using the levels (during the corresponding sampling week for each school) measured by a background air pollution monitoring station in Barcelona to remove temporal fluctuation in background levels from our analyses using a method that has been reported previously [31].

Data analysis

We used cross-sectional and longitudinal frameworks to analyze the association between air pollution exposure and spectacles use. For the cross-sectional analyses, we applied spectacles use at baseline (i.e. prevalent spectacles use) as the outcome variable and lifelong exposure to residential air pollution as the main exposure variable. The longitudinal framework was based on the incidence of spectacles use and exposure to air pollution during the three-year period (2012–2015) between the baseline and last data collection campaigns. Accordingly, we developed a binary variable indicating whether the participant started to use spectacles during this period and we used this variable as the outcome variable together with prospective exposure to residential air pollution as the main exposure variable for the longitudinal analyses. For the exposure to air pollution at schools, we used the annual estimate (2009) for both cross-sectional and longitudinal analyses.

Because of the multilevel nature of the data (children at schools), we used logistic mixed effects models with prevalent/incident spectacles use (one at a time) as the outcome variable, estimates of lifelong/prospective exposure to each pollutant at home and school (one at a time) as the fixed effect predictor, and school as random effect. For the longitudinal analyses, we excluded those participants using spectacles at baseline since they could not be considered to be at the risk of using spectacles during the course of the follow-up. The analyses were adjusted for a number of covariates identified a priori: age (at the time of the baseline data collection for the cross-sectional analyses and the time of the last data collection campaign for the longitudinal analyses), sex, paternal and maternal ethnicities (European or non-European), prematurity (yes/no) [3234], child’s exposure to environmental tobacco smoke (yes/no) [34, 35], child’s average screen time per week, child annual total time (hours) of playing in green spaces, and indicators of socioeconomic status (SES) at both individual and area levels. In particular, we used weekly screen time as a surrogate for ‘near-work’ and child annual total time of playing in green spaces as a surrogate for outdoor activity which have been associated with the risk of myopia [1]. Screen time and green space playing time were reported by parents as the average time (separately for working days and weekends) the child would spend on watching TV or playing game on videogame console or computer and the average time the child would spend playing in green spaces separately for the working days and weekends during school period and summer holidays. Paternal and maternal educational achievements (primary school, secondary school, or university) was used as the indicator of individual-level SES and Urban Vulnerability Index [36], a measure of neighborhood SES at the census tract (median area of 0.08 km2 for the study region) was applied as the indicator of area-level SES. The odds ratios (ORs) were reported for an interquartile range (IQR) increase in each pollutant based on all study participants separately for cross-sectional and longitudinal analyses.

Results

Of 2727 (94.1%) BREATHE participants with available data on spectacles use who were included in the cross-sectional analyses, 359 (13.2%) used glasses at baseline. As presented in Table 1, those participants using spectacles were more likely to be girls (p = 0.07), older (p<0.01), and spending less time playing outdoors in green spaces (p = 0.06).

thumbnail
Table 1. Descriptiona of characteristics of the study participants.

https://doi.org/10.1371/journal.pone.0167046.t001

Of BREATHE participants included in the cross-sectional analyses, 1812 (66.5%) were followed till the last data collection campaign and included in the longitudinal analyses. The last data collection campaign identified 155 incident cases of spectacles use (i.e. 155 BREATHE participants started to use spectacles during our three-year follow-up period between the first and last data collection campaigns). Compared to participants included in the cross-sectional analyses, those included in the longitudinal analyses were younger (p<0.01) at baseline (S1 Table) which was expected because older children were supposed to finish primary school and move to high school before our last data collection campaign (This was the main reason for lost to follow-up in our study). The participants included in the longitudinal analyses also reported less screen time (p<0.01) and were more likely to be of European descent (p<0.01).

The Spearman’s correlation coefficients between different exposures are presented in S2 Table. While levels of different pollutants at school and at home were strongly correlated, the correlation between school and home levels of each pollutant was weak to moderate.

Cross-sectional analyses

An IQR increase in NO2 level at home was associated with 16% (95% confidence intervals (CIs): 3%, 29%) increase in spectacles use (Table 2). Similarly, we observed an increase in the risk of using spectacles associated with one-IQR increase in exposure to PM2.5 absorbance at home but the association was marginally statistically significant.

thumbnail
Table 2. Median (InterQuartile Range, IQR) of air polltants and adjusteda odds ratio (95% confidence intervals) of the use of spectacles associated with one Inter-Quartile Range (IQR) increase in exposure to each pollutant.

https://doi.org/10.1371/journal.pone.0167046.t002

An IQR increase in NO2 level at school was associated with 32% (95% CIs: 9%, 59%) increase in the risk of wearing spectacles. BC exposure at school was also associated with spectacles use but the association did not attain statistical significance.

Longitudinal analyses

For residential NO2 exposure we observed identical association in terms of direction and strength with that of the cross-sectional analysis (Table 2); however, as expected (because of smaller sample size), the 95% confidence intervals (CI) were wider and include one (p-value = 0.06). The longitudinal association for residential exposure to PM2.5 absorbance became stronger compared to that of cross-sectional analyses and attained statistical significance (Table 2). Similarly, the association for school BC exposure was stronger and statistically significant in the longitudinal analyses (Table 2). On the other hand, the association for NO2 exposure at school became weaker and lost its statistical significant in longitudinal analyses.

The estimates for 10 unit increase in NO2 (μg/m3) and one unit increase in PM2.5 Absorbance (10−5/m3) and Black Carbon (μg/m3) for cross-sectional and longitudinal analyses are presented in S3 Table.

Sensitivity analyses

Further adjustment of analyses for maternal smoking during pregnancy [33, 37, 38], breastfeeding [34], neighborhood socioeconomic status (Urban Vulnerability Index) of the school, parental employment status (unemployed, employee, or self-employed), parental marital status, and child’s height [34] did not change our findings notably (Data not shown). However, after adjustment of cross-sectional analyses for parental employment status, the association for residential PM2.5 absorbance became stronger and attained statistical significance (OR: 1.13, 95% CIs: 1.00, 1.29). Limiting the cross-sectional analyses of residential air pollution to those who had not moved since birth (n = 1689, 61.9%) did not result in a considerable change in direction and strength of the associations, but, as expected, the confidence intervals for NO2 became wider and the association became nearly statistically significant (OR: 1.14, 95% CI: 0.99, 1.32, p-value = 0.07). Additionally, we did not observe any statistically significant effect modification by child’s sex, maternal education or neighborhood SES for our associations.

Discussion

Interpretation of results

To our knowledge, this is the first study to evaluate the association between air pollution exposure and the use of spectacles (a marker of refractive errors of vision and specifically myopia). We observed an increased likelihood of spectacles use associated with higher exposure to traffic-related air pollutants that were generally consistent for exposures at home and at school and in cross-sectional and longitudinal analytical frameworks. These findings were also robust against a range of sensitivity analyses that we conducted.

Further adjustment of our analyses for SES indicators other than educational attainments of parents and residential neighborhood SES such as parental employment status, marital status and school neighborhood SES did not result in a notable change in our findings. Moreover, the indicators of SES (parental educational attainment and neighborhood SES) were not associated with the risk of spectacles use (Table 1). These observations could suggest that our analyses were not likely to have been influenced by residual SES confounding.

While the results of residential exposure to air pollution were consistent in cross-sectional and longitudinal analyses, for air pollution exposure at school we observed a difference between findings of these analytical frameworks. The findings for the cross-sectional analyses of school exposures need to be interpreted with caution because most of the study participants were recruited in their first years of primary schools, thus for them the length of exposure to school air pollution before reporting spectacles use at baseline could have been too short to be able to induce refractive errors. On the other hand, the consistency of our findings (in terms of direction and strength of associations) for the exposure to air pollution at school in the longitudinal analyses with those of residential exposure to air pollution in both cross-sectional and longitudinal analyses could offer us more confidence about these findings.

Available evidence and potential underlying mechanisms

We are not aware of any previous epidemiological studies on our investigated association between air pollution and refractive errors; therefore, it is not possible to compare our findings with those of others. There are also no available animal model or in vitro or in vivo studies evaluating the direct link between exposure to air pollution and refractive errors of vision. However, our findings are consistent with a number of previous observations. A recent study has reported a higher risk of “near visual difficulty” (defined as difficulty in seeing and recognizing an object at arm’s length) associated with having a cooking stove in the same room as sleeping area which could be an indicator of indoor air pollution [39].

Local and systemic inflammation and oxidative stress are the most established mechanisms for the adverse health effects of air pollution. Systemic inflammatory diseases often affect different parts of the eyes, including the sclera, cornea, vitreous, and retina, resulting in intraocular inflammatory conditions such as uveitis and retinal vasculitis [40]. Similarly, systemic inflammation induced by long-term exposure to air pollution can prompt changes in retinal microvasculature, including narrowing of arteriolar diameters and widening of venular diameters, which in turn have been associated with arteriolar damage, endothelial dysfunction, and intraocular inflammation [4145]. The latter can induce myopia in the eyes by affecting the optical power of the lens [46] and/or impairing retinal neuroactivity [47], which regulates the axial growth of eyes early in life. Inadequate or excessive axial growth of eyes relative to the ocular refractive power is one of the known causes of myopia. In addition to inducing systemic inflammation, exposure to air pollution generates local inflammation on the ocular surface [13, 48]. Animal studies have shown that eye surface inflammation can infiltrate into the eye [49], resulting in retinal inflammation [50], which in turn can affect its neuroactivity and regulation of axial length growth.

In addition to inflammation, air pollution induces oxidative stress in the eyes, which has been reported to be involved in a number of conditions such as cataract, uveitis, age-related macular degeneration, glaucoma, and various types of retinopathy [51, 52]. Oxidative stress has been shown to impair the release of dopamine from retinal cells [53], which plays a critical role in regulating the axial growth of the eye [1, 8, 54]. Therefore, oxidative stress-related impairment of dopamine release from retinal cells could be one explanation for the impact of air pollution on myopia.

Furthermore, air pollution has been implicated as a risk factor for dry eye disease by inducing instability of tear film, ocular surface inflammation, epithelial differentiation and hyperplasia of goblets cells [13, 48]. Dry eye disease has been shown to result in reduction in corneal thickness [55] and irregularities in corneal surface [56] which can ultimately lead to impaired visual acuity [57, 58]. Dry eye disease has also been reported to increase oxidative stress in conjunctival epithelium [59] which could contribute to the aforementioned oxidative stress pathway.

Limitations of study

The generalizability of our findings might have been affected by selection bias in that those participants participated in BREATHE might be different from those not participated with respect to SES. The Urban Vulnerability Index of the schools was not associated with school participation rate (Spearman’s correlation coefficient = -0.11, p-value = 0.51); this might suggest that the SES was less likely to be a major predictor of participating in the study. A part from non-participation, we had additional loss to follow-up from baseline to the end of the study. Those participants dropped out of the study during the follow-up period were different from those followed until the end of the study in terms of age, screen time and ethnicity, and maternal education which could have introduced selection bias in our findings. As the differing characteristics are not expected to be associated to residential air pollution levels, we do not expect this loss to follow-up to bias the results. We used eyeglasses as a surrogate for myopia. However, the use of eyeglasses by our participants could also have been due to other refractive errors such as hyperopia and astigmatism that have different pathogenesis. This could be more relevant for our cross-sectional analyses based on the prevalent eyeglasses use. On the other hand, considering the age of our study participants, we expect that the incident reduction in visual acuity during the course of our longitudinal study is more strongly associated with myopia than other refractive errors [1517]. Therefore, the findings of our longitudinal analyses were more likely to be relevant to myopia compared to those of cross-sectional analyses. Moreover, children with less severe refractive problems who did not use spectacles were not considered as having refractive errors in our study which could have biased our findings towards null. Likewise, we did not obtain information on use of contact lenses and the resulting outcome misclassification could have biased our estimates towards null. Furthermore, we did not have data on refractive errors in parents enabling us to address the genetic contribution in our analyses. However, we do not have any reason to assume a differential exposure to air pollution for children of parents with and without refractive error. Moreover, by temporally adjusting the LUR spatial estimates of pollutant levels, we effectively assumed that the city spatial surface and the spatial distribution of pollutants remained unchanged over the study period. Previous studies in Europe have shown the stability of these spatial contrasts over a long period [60, 61]. Moreover, we are not aware of any major change in traffic flow, land use, or emissions profiles occurred between the year of LUR model construction and the study period. Accordingly, we observed a strong correlation (Spearman’s correlation coefficient of 0.81) between modeled NO2 levels at school using LUR models and measured NO2 levels at schools during BREATHE campaigns, assuring us about the long-term validity the assigned home exposure levels.

Conclusions

We found an increase in the likelihood of myopia (as surrogated by spectacles use) in association with exposure to traffic-related air pollution at home and at school. Because of the aforementioned limitations, this study might not be able to establish a causal link; however, considering the consistent pattern of our observed associations for school and residential exposures while they were weakly correlated, the consistency of our findings based on cross-sectional and longitudinal analyses, and the robustness of these associations to several sensitivity analyses, we are convinced that our findings merits further investigation.

Currently, about half of the world population reside in cities, and it is predicted that by 2030 around 70% of global population will live in urban areas [62] where myopia is more prevalent. According to a recently published report by the WHO, most of the world’s cities (mainly in developing countries) are currently in breach of its guidelines on air pollution levels [63]. Uncorrected refractive errors are a major contributor to the global burden of disease accounting for more than 11 million disability-adjusted life-years (DALYs) [64]. Taking into account such a considerable burden of refractive errors with the ongoing rise in their prevalence worldwide, an adverse impact of air pollution on these conditions, if established by future studies, offers policymakers an evidence base for developing policies and implementing targeted interventions in order to slow down and ideally reverse the current ongoing rise in prevalence of eye refractive errors. Such an impact could also open a whole new area in our understanding of the causes of eye refractive errors in general and myopia in particular as well as adverse health effects of air pollution which are of great importance for research community in various disciplines. Further animal, in vivo, and in vitro studies are required to elucidate potential pathways underling such an impact, if any. We advise future epidemiological studies to apply more refined outcome measures such as quantified visual acuity and to separate different types of refractive errors.

Supporting information

S1 Table. Description of characteristics of the study participants included in cross-sectional and longitudinal analyses.

https://doi.org/10.1371/journal.pone.0167046.s001

(DOCX)

S2 Table. Spearman’s correlation coefficient between estimates of air pollution for the cross-sectional analyses.

https://doi.org/10.1371/journal.pone.0167046.s002

(DOCX)

S3 Table. Adjusted odds ratio (95% confidence intervals) of the use of spectacles associated with 10 unit increase in NO2 (μg/m3) and one unit increase in PM2.5 Absorbance (10−5/m3) and Black Carbon (μg/m3).

https://doi.org/10.1371/journal.pone.0167046.s003

(DOCX)

Acknowledgments

We are grateful to all the families and schools participating in the study for their altruism and their collaboration. The research leading to these results has received funding from the European Research Council under the ERC Grant Agreement number 268479 –the BREATHE project. Payam Dadvand is funded by a Ramón y Cajal fellowship (RYC-2012-10995) awarded by the Spanish Ministry of Economy and Competitiveness.

Author Contributions

  1. Conceptualization: PD MJN BB XQ IGM JS.
  2. Data curation: MC IR MAP.
  3. Formal analysis: ADB XB.
  4. Funding acquisition: JS BB.
  5. Investigation: MAP MC JS.
  6. Methodology: PD XB.
  7. Project administration: MAP JS.
  8. Resources: MAP.
  9. Supervision: IGM.
  10. Writing – original draft: PD.
  11. Writing – review & editing: PD MJN XB MAP ADB MC IR BB XQ IGM JS.

References

  1. 1. Morgan IG, Ohno-Matsui K, Saw S-M. Myopia. Lancet. 2012;379(9827):1739–48. pmid:22559900
  2. 2. Dolgin E. The myopia boom. Nature. 2015;519(7543):276–8. pmid:25788077
  3. 3. Holden B, Sankaridurg P, Smith E, Aller T, Jong M, He M. Myopia, an underrated global challenge to vision: where the current data takes us on myopia control. Eye. 2014;28(2):142–6. pmid:24357836
  4. 4. Pan C-W, Ramamurthy D, Saw S-M. Worldwide prevalence and risk factors for myopia. Ophthalmic Physiol Opt. 2012;32(1):3–16. pmid:22150586
  5. 5. Goldschmidt E, Jacobsen N. Genetic and environmental effects on myopia development and progression. Eye. 2014;28(2):126–33. pmid:24357837
  6. 6. Wojciechowski R. Nature and nurture: the complex genetics of myopia and refractive error. Clin Genet. 2011;79(4):301–20. pmid:21155761
  7. 7. Morgan I, Rose K. How genetic is school myopia? Prog Retin Eye Res. 2005;24(1):1–38. pmid:15555525
  8. 8. Lougheed T. Myopia: the evidence for environmental factors. Environmental health perspectives. 2014;122(1):A12. pmid:24380886
  9. 9. Ip JM, Rose KA, Morgan IG, Burlutsky G, Mitchell P. Myopia and the Urban Environment: Findings in a Sample of 12-Year-Old Australian School Children. Invest Ophthalmol Vis Sci. 2008;49(9):3858–63. pmid:18469186
  10. 10. Wolfram C, Höhn R, Kottler U, Wild P, Blettner M, Bühren J, et al. Prevalence of refractive errors in the European adult population: the Gutenberg Health Study (GHS). Br J Ophthalmol 2014;98(7):857–61. pmid:24515986
  11. 11. Katz J, Tielsch JM, Sommer A. Prevalence and risk factors for refractive errors in an adult inner city population. Invest Ophthalmol Vis Sci. 1997;38(2):334–40. pmid:9040465
  12. 12. Lim SS, Vos T, Flaxman AD, Danaei G, Shibuya K, Adair-Rohani H, et al. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2013;380(9859):2224–60.
  13. 13. Torricelli AAM, Novaes P, Matsuda M, Alves MR, Monteiro MLR. Ocular surface adverse effects of ambient levels of air pollution. Arq Bras Oftalmol. 2011;74(5):377–81. pmid:22184003
  14. 14. Bruce N, Perez-Padilla R, Albalak R. Indoor air pollution in developing countries: a major environmental and public health challenge. Bull World Health Organ. 2000;78(9):1078–92. pmid:11019457
  15. 15. Leone J, Mitchell P, Morgan I, Kifley A, Rose K. Use of visual acuity to screen for significant refractive errors in adolescents: Is it reliable? Arch Ophthalmol. 2010;128(7):894–9. pmid:20625052
  16. 16. Cumberland PM, Peckham CS, Rahi JS. Inferring myopia over the lifecourse from uncorrected distance visual acuity in childhood. Br J Ophthalmol. 2007;91(2):151–3. pmid:17020900
  17. 17. Xiang F, He M, Zeng Y, Mai J, Rose KA, Morgan IG. Increases in the prevalence of reduced visual acuity and myopia in Chinese children in Guangzhou over the past 20 years. Eye. 2013;27(12):1353–8. pmid:24008929
  18. 18. Tong L, Saw S-M, Carkeet A, Chan W-Y, Wu H-M, Tan aD. Prevalence Rates and Epidemiological Risk Factors for Astigmatism in Singapore School Children. Optom Vis Sci. 2002;79(9):606–13. pmid:12322931
  19. 19. Guo Y, Liu LJ, Xu L, Tang P, Lv YY, Feng Y, et al. Myopic Shift and Outdoor Activity among Primary School Children: One-Year Follow-Up Study in Beijing. PLoS ONE. 2013;8(9):e75260. pmid:24086484
  20. 20. Saw S-M, Shankar A, Tan S-B, Taylor H, Tan DTH, Stone RA, et al. A Cohort Study of Incident Myopia in Singaporean Children. Invest Ophthalmol Vis Sci. 2006;47(5):1839–44. pmid:16638989
  21. 21. Williams C, Northstone K, Howard M, Harvey I, Harrad RA, Sparrow JM. Prevalence and risk factors for common vision problems in children: data from the ALSPAC study. Br J Ophthalmol. 2008;92(7):959–64. pmid:18480306
  22. 22. Guggenheim JA, Northstone K, McMahon G, Ness AR, Deere K, Mattocks C, et al. Time Outdoors and Physical Activity as Predictors of Incident Myopia in Childhood: A Prospective Cohort Study. Invest Ophthalmol Vis Sci. 2012;53(6):2856–65. pmid:22491403
  23. 23. He M, Xiang F, Zeng Y, Mai J, et al. Effect of time spent outdoors at school on the development of myopia among children in china: A randomized clinical trial. JAMA. 2015;314(11):1142–8. pmid:26372583
  24. 24. Ajuntament de Barcelona. Dades bàsiques de mobilitat 2011 [in Catalan] Barcelona2012 [cited 2013 22 February]. http://w110.bcn.cat/Mobilitat/Continguts/Documents/Fitxers/dadesbasiques2011complert.pdf.
  25. 25. Amato F, Pandolfi M, Escrig A, Querol X, Alastuey A, Pey J, et al. Quantifying road dust resuspension in urban environment by multilinear engine: a comparison with PMF2. Atmospheric Environment. 2009;43(17):2770–80.
  26. 26. Sunyer J, Esnaola M, Alvarez-Pedrerol M, Forns J, Rivas I, López-Vicente M, et al. Association between Traffic-Related Air Pollution in Schools and Cognitive Development in Primary School Children: A Prospective Cohort Study. PLoS Med. 2015;12(3):e1001792. pmid:25734425
  27. 27. Dadvand P, Nieuwenhuijsen MJ, Esnaola M, Forns J, Basagaña X, Alvarez-Pedrerol M, et al. Green spaces and cognitive development in primary schoolchildren. Proc Natl Acad Sci U S A. 2015;112(26):7937–42. pmid:26080420
  28. 28. Beelen R, Hoek G, Vienneau D, Eeftens M, Dimakopoulou K, Pedeli X, et al. Development of NO2 and NOx land use regression models for estimating air pollution exposure in 36 study areas in Europe—the ESCAPE project. Atmos Environ. 2013;72:10–23.
  29. 29. Eeftens M, Beelen R, de Hoogh K, Bellander T, Cesaroni G, Cirach M, et al. Development of land use regression models for PM2.5, PM2.5 absorbance, PM10 and PMcoarse in 20 European study areas; results of the ESCAPE project. Environ Sci Technol. 2012;46(20):11195–205. pmid:22963366
  30. 30. Dadvand P, Figueras F, Basagaña X, Beelen R, Martinez D, Cirach M, et al. Ambient Air Pollution and Preeclampsia: a Spatiotemporal Analysis. Environ Health Perspect. 2013;121(11–12):1365–71. pmid:24021707
  31. 31. Rivas I, Viana M, Moreno T, Pandolfi M, Amato F, Reche C, et al. Child exposure to indoor and outdoor air pollutants in schools in Barcelona, Spain. Environ Int. 2014;69:200–12. pmid:24875803
  32. 32. Chen T-C, Tsai T-H, Shih Y-F, Yeh P-T, Yang C-H, Hu F-C, et al. Long-term Evaluation of Refractive Status and Optical Components in Eyes of Children Born Prematurely. Invest Ophthalmol Vis Sci. 2010;51(12):6140–8. pmid:20688740
  33. 33. Rahi JS, Cumberland PM, Peckham CS. Myopia Over the Lifecourse: Prevalence and Early Life Influences in the 1958 British Birth Cohort. Ophthalmology. 2011;118(5):797–804. pmid:21185080
  34. 34. Beuerman RW, Saw SM, Tan DTH, Wong T-Y. Myopia: animal models to clinical trials. Singapore: World Scientific Publishing Co. Pte. Ltd.; 2010.
  35. 35. El-Shazly AA. Passive smoking exposure might be associated with hypermetropia. Ophthalmic Physiol Opt. 2012;32(4):304–7. pmid:22670861
  36. 36. Spanish Ministry of Public Works. Atlas of Urban Vulnerability in Spain. Methodology and contents. Madrid: Spanish Ministry of Public Works; 2012 [cited 2016 27 May]. http://www.fomento.gob.es/NR/rdonlyres/40668D5E-26B6-4720-867F-286BD55E1C6B/135960/20160201METODOLOGIAATLASVULNERABILIDAD2001Y2011.pdf.
  37. 37. Borchert MS, Varma R, Cotter SA, Tarczy-Hornoch K, McKean-Cowdin R, Lin JH, et al. Risk Factors for Hyperopia and Myopia in Preschool Children: The Multi-Ethnic Pediatric Eye Disease and Baltimore Pediatric Eye Disease Studies. Ophthalmology. 2011;118(10):1966–73. pmid:21856013
  38. 38. Fernandes M, Yang X, Li JY, Cheikh Ismail L. Smoking during pregnancy and vision difficulties in children: a systematic review. Acta Ophthalmologica. 2014; n/a–n/a.
  39. 39. Freeman EE, Zunzunegui M-V, Kouanda S, Aubin M-Je, Popescu ML, Miszkurka M, et al. Prevalence and risk factors for near and far visual difficulty in Burkina Faso. Ophthalmic Epidemiol. 2010;17(5):301–6. pmid:20868256
  40. 40. McCluskey P, Powell RJ. The eye in systemic inflammatory diseases. The Lancet. 2004;364(9451):2125–33.
  41. 41. Adar SD, Klein R, Klein BEK, Szpiro AA, Cotch MF, Wong TY, et al. Air Pollution and the Microvasculature: A Cross-Sectional Assessment of In Vivo Retinal Images in the Population-Based Multi-Ethnic Study of Atherosclerosis (MESA). PLoS Med. 2010;7(11):e1000372. pmid:21152417
  42. 42. Wong TY, Amirul Islam FM, Klein R, Klein BEK, Cotch MF, Castro C, et al. Retinal Vascular Caliber, Cardiovascular Risk Factors, and Inflammation: The Multi-Ethnic Study of Atherosclerosis (MESA). Invest Ophthalmol Vis Sci. 2006;47(6):2341–50. pmid:16723443
  43. 43. Nguyen TT, Wang JJ, Sharrett AR, Islam FMA, Klein R, Klein BEK, et al. Relationship of Retinal Vascular Caliber With Diabetes and Retinopathy: The Multi-Ethnic Study of Atherosclerosis (MESA). Diabetes Care. 2008;31(3):544–9. pmid:18070990
  44. 44. Kolodjaschna J, Berisha F, Lung S, Schaller G, Polska E, Jilma B, et al. LPS-induced microvascular leukocytosis can be assessed by blue-field entoptic phenomenon. Am J Physiol Heart Circ Physiol. 2004;287(2):H691–H4. pmid:15016626
  45. 45. Louwies T, Panis LI, Kicinski M, De Boever P, Nawrot TS. Retinal microvascular responses to short-term changes in particulate air pollution in healthy adults. Environ Health Perspect. 2013;121(9):1011–6. pmid:23777785
  46. 46. Herbort CP, Papadia M, Neri P. Myopia and Inflammation. J Ophthalmic Vis Res. 2011;6(4):270–83. pmid:22454750
  47. 47. Bronzetti E, Artico M, Kovacs I, Felici LM. Expression of neurotransmitters and neurotrophins in neurogenic inflammation of the rat retina. Eur J Histochem. 2009;51(4):251–60.
  48. 48. Tau J, Novaes P, Matsuda M, Tasat DR, Saldiva PH, Berra A. Diesel Exhaust Particles Selectively Induce Both Proinflammatory Cytokines and Mucin Production in Cornea and Conjunctiva Human Cell Lines. Invest Ophthalmol Vis Sci. 2013;54(7):4759–66. pmid:23722391
  49. 49. Williams RN, Paterson CA, Eakins KE, Bhattacherjee P. Quantification of ocular inflammation: Evaluation of polymorphonuclear leucocyte infiltration by measuring myeloperoxidase activity. Curr Eye Res. 1982;2(7):465–70. pmid:6303695
  50. 50. Chinnery HR, McLenachan S, Binz N, Sun Y, Forrester JV, Degli-Esposti MA, et al. TLR9 Ligand CpG-ODN Applied to the Injured Mouse Cornea Elicits Retinal Inflammation. Am J Pathol. 2012;180(1):209–20. pmid:22085974
  51. 51. Saccà SC, Roszkowska AM, Izzotti A. Environmental light and endogenous antioxidants as the main determinants of non-cancer ocular diseases. Mutat Res. 2013;752(2):153–71. pmid:23337404
  52. 52. Wakamatsu TH, Dogru M, Tsubota K. Tearful relations: oxidative stress, inflammation and eye diseases. Arq Bras Oftalmol. 2008;71:72–9. pmid:19274416
  53. 53. Njie-Mbye YF, Kulkarni-Chitnis M, Opere CA, Barrett A, Ohia SE. Lipid peroxidation: pathophysiological and pharmacological implications in the eye. Front Physiol. 2013;4:366. pmid:24379787
  54. 54. Rymer J, Wildsoet CF. The role of the retinal pigment epithelium in eye growth regulation and myopia: A review. Vis Neurosci. 2005;22(03):251–61.
  55. 55. Liu Z, Pflugfelder SC. Corneal Thickness Is Reduced in Dry Eye. Cornea. 1999;18(4):403–7. pmid:10422850
  56. 56. Liu Z, Pflugfelder SC. Corneal surface regularity and the effect of artificial tears in aqueous tear deficiency. Ophthalmology. 1999;106(5):939–43. pmid:10328393
  57. 57. Goto E, Yagi Y, Matsumoto Y, Tsubota K. Impaired functional visual acuity of dry eye patients. Am J Ophthalmol. 2002;133(2):181–6. pmid:11812420
  58. 58. International Dry Eye WorkShop. The epidemiology of dry eye disease: report of the epidemiological subcommittee of the international dry eye workshop. Ocul Surf. 2007;5(2):93–107. pmid:17508117
  59. 59. Macri A, Scanarotti C, Bassi A, Giuffrida S, Sangalli G, Traverso C, et al. Evaluation of oxidative stress levels in the conjunctival epithelium of patients with or without dry eye, and dry eye patients treated with preservative-free hyaluronic acid 0.15% and vitamin B12 eye drops. Graefes Arch Clin Exp Ophthalmol. 2014;
  60. 60. Eeftens M, Beelen R, Fischer P, Brunekreef B, Meliefste K, Hoek G. Stability of measured and modelled spatial contrasts in NO2 over time. Occup Environ Med. 2011;68(10):765–70. pmid:21285243
  61. 61. Cesaroni G, Porta D, Badaloni C, Stafoggia M, Eeftens M, Meliefste K, et al. Nitrogen dioxide levels estimated from land use regression models several years apart and association with mortality in a large cohort study. Environ Health. 2012;11:48. pmid:22808928
  62. 62. Martine G, Marshall A. State of world population 2007: unleashing the potential of urban growth. New York: UNFPA; 2007.
  63. 63. World Health Organization. Ambient (outdoor) air quality and health. Fact sheet N°313. 2014;http://www.who.int/mediacentre/factsheets/fs313/en/.
  64. 64. Vos T, Barber RM, Bell B, Bertozzi-Villa A, Biryukov S, Bolliger I, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 301 acute and chronic diseases and injuries in 188 countries, 1990–2013;2013: a systematic analysis for the Global Burden of Disease Study 2013. Lancet. 2015;386(9995):743–800.