Skip to main content

Greater exposure to PM2.5 and PM10 was associated with lower corneal nerve measures: the Maastricht study - a cross-sectional study

Abstract

Background

Epidemiological and toxicological studies indicate that increased exposure to air pollutants can lead to neurodegenerative diseases. To further confirm this relationship, we evaluated the association between exposure to ambient air pollutants and corneal nerve measures as a surrogate for neurodegeneration, using corneal confocal microscopy.

Methods

We used population-based observational cross-sectional data from The Maastricht Study including N = 3635 participants (mean age 59.3 years, 51.6% were women, and 19.9% had type 2 diabetes) living in the Maastricht area. Using the Geoscience and hEalth Cohort COnsortium (GECCO) data we linked the yearly average exposure levels of ambient air pollutants at home address-level [particulate matter with diameters of ≤ 2.5 µm (PM2.5), and ≤ 10.0 µm (PM10), nitrogen dioxide (NO2), and elemental carbon (EC)]. We used linear regression analysis to study the associations between Z-score for ambient air pollutants concentrations (PM2.5, PM10, NO2, and EC) and Z-score for individual corneal nerve measures (corneal nerve bifurcation density, corneal nerve density, corneal nerve length, and fractal dimension).

Results

After adjustment for potential confounders (age, sex, level of education, glucose metabolism status, corneal confocal microscopy lag time, inclusion year of participants, smoking status, and alcohol consumption), higher Z-scores for PM2.5 and PM10 were associated with lower Z-scores for corneal nerve bifurcation density, nerve density, nerve length, and nerve fractal dimension [stβ (95% CI): PM2.5 -0.10 (-0.14; -0.05), -0.04 (-0.09; 0.01), -0.11 (-0.16; -0.06), -0.20 (-0.24; -0.15); and PM10 -0.08 (-0.13; -0.03), -0.04 (-0.09; 0.01), -0.08 (-0.13; -0.04), -0.17 (-0.21; -0.12)], respectively. No associations were found between NO2 and EC and corneal nerve measures.

Conclusions

Our population-based study demonstrated that exposure to higher levels of PM2.5 and PM10 are associated with higher levels of corneal neurodegeneration, estimated from lower corneal nerve measures. Our results suggest that air pollution may be a determinant for neurodegeneration assessed in the cornea and may impact the ocular surface health as well.

Peer Review reports

Introduction

Air pollution significantly impacts health, with an estimated seven million annual deaths in 2021, according to the World Health Organization (WHO) [1]. Efforts have been made to improve air quality by imposing stricter guidelines set by the WHO [1].

Air pollution is a complex mixture of gasses, solid particles, and liquid droplets. The pollutants that have the strongest supporting evidence regarding public health issues encompass particulate matter (PM), ozone, nitrogen dioxide (NO2), sulfur dioxide, and carbon monoxide [1]. The broad spectrum complications attributed to air pollution extends beyond respiratory issues [2] to ocular health [3] and neurodegeneration [4]. Previous studies investigated the impact of air pollution ambient level exposure on alterations of the ocular surface [5], dry eye syndrome [6], allergic conjunctivitis [7, 8], keratoconus [9] and corneal morphology [3, 10,11,12]. Moreover, studies have shown that air pollution can cause neuroinflammation and neurodegeneration [13].

Despite the cornea's direct vulnerability to environmental factors, there is currently a notable gap in our understanding of the association between ambient air pollutant exposure and corneal nerve morphology in well-characterized study populations. The integrity of the corneal nerves can serve as a valuable window for observing corneal health [14], and central and peripheral neurodegenerative disorders [15,16,17,18,19,20]. Studying this association is important because alterations in corneal sensitivity due to corneal neurodegeneration can lead to dry eye, delayed wound healing after injury, infections, and persistent epithelial defects [21] which may result in severe vision loss or even blindness [22]. Moreover, as postulated in the ticking clock hypothesis, the onset and/or progression of early neurodegeneration can be delayed by reducing exposure to potentially modifiable risk factors [23, 24]. To the best of our knowledge, the association between exposure to air pollutants and neurodegeneration assessed in the cornea has not been studied. In the present paper, we focus on the exposure to air pollutants as a determinant of morphological changes of corneal nerve fibers assessed using the in vivo, non-invasive and sensitive ophthalmic technique of corneal confocal microscopy [25].

Using data from The Maastricht Study we examined the associations of exposure to ambient air pollution estimated at home address with corneal nerve fiber measures, as assessed using corneal confocal microscopy.

Methods

Study population and design

We used data from the Maastricht Study, an observational prospective population-based cohort study enriched with type 2 diabetes individuals, for efficiency reasons. The rationale and methodology have been described previously [26]. In brief, the study focuses on the etiology, pathophysiology, complications and comorbidities of type 2 diabetes mellitus, and is characterized by an extensive phenotyping approach. Eligible for participation were all individuals aged between 40 and 75 years living in the southern part of the Netherlands. Participants were recruited from the Maastricht area, as shown in Fig. 1, through mass media campaigns and from the municipal registries and the regional Diabetes Patient Registry via mailings. Recruitment was stratified according to known type 2 diabetes status [26]. We had cross-sectional data for N = 9187 participants who were included for the baseline survey between November 2010 and October 2020. The examinations for each participant were performed within a time window of 3 months. Estimates of address-level yearly average levels of ambient air pollutants were provided by the Geoscience and Health Cohort Consortium (GECCO) [27]. The study was approved by the institutional medical ethical committee (NL31329.068.10) and the Minister of Health, Welfare and Sports of the Netherlands (permit 131088–105234-PG). All participants gave written informed consent [26].

Fig. 1
figure 1

Geographical Area of Participant Recruitment: Maastricht Region

Assessment of corneal confocal microscopy measurements

We utilized the Heidelberg retina tomograph III (HRT3) along with the Rostock cornea module (Heidelberg Engineering, Germany) for imaging the corneal nerves [25] of the left eye exclusively, due to logistical considerations. Prior to the assessment, we applied oxybuprocaine hydrochloride 0.4 mg/ml eye drops (Minims; Bausch & Lomb, France) to both eyes to prevent blinking, and applied carbomer 2 mg/g eye gel (Vidsic Bausch & Lomb) to both eyes to ensure optimal contact between the cornea and the applanation cap [20].

We performed large-scale corneal confocal microscopy imaging (1600 × 1600 μm). We employed the U-net-based convolutional neural network for fully automated tracing of the corneal nerves [28] (example can be found in Supplemental Fig. S1). Traced corneal nerve images were reviewed manually to ensure quality criteria were met (example can be found in Supplemental Fig. S2), and the location of the captured images was assessed based on the orientation of the corneal nerves (Supplemental Fig. S3). Additional details regarding this process can be found in the supplementary methods.

Our analysis included the following indices: corneal nerve bifurcation density (number of bifurcation points or branching points per mm2), corneal nerve density (total number of corneal nerve fibers, including main fibers and branches, per mm2; 'main nerve fibers' denoting the largest and most prominent nerve fibers), corneal nerve length (total length of corneal nerve fibers in mm, including both main fibers and branches, per mm2), and corneal nerve fractal dimension (a measure of nerve structure complexity, unit-less). Additional information can be found in the supplementary methods.

Assessment of air pollution

The assessment of air pollution exposure cannot be captured by a single pollutant due to the diversity of sources and dispersion patterns. Consequently, we selected particulate matter with diameters of ≤ 2.5 µm (PM2.5), and ≤ 10.0 µm (PM10), nitrogen dioxide (NO2), and elemental carbon (EC) as key pollutants for evaluation, given their regulatory significance and health implications [29].

The Institute for Public Health and the Environment (RIVM) models and maps the average concentrations of these pollutants every year, using a grid resolution of 25 m. These maps were derived from 1 km resolution nationwide background concentration maps, which were further refined with local traffic data. In essence, the nationwide background maps incorporated dispersion models that included data on emissions from industrial, vehicular, and household sources both within the Netherlands and from neighboring regions, as well as meteorological and chemical information [30]. The local traffic data, sourced from the Dutch Nationaal Samenwerkingsprogramma Luchtkwaliteit, were integrated with the national background maps using two specific models, one for urban roads and another for highways in more open areas [31, 32].

Air pollution concentrations were then computed at 9 million data points across the Netherlands and interpolated to a 25 m resolution raster map using ordinary Kriging, the data are publicly available [33]. The GECCO linked these air pollution concentrations to residential addresses and exported the data into a tabular format. Within GECCO, estimated concentrations of air pollutants were available for the years 2013 to 2020 for PM2.5, PM10, and EC, or 2014 to 2020 for NO2. While the model predictions and absolute concentrations for NO2 and PM2.5 generally aligned well with measured values [34], the interpretation of quality metrics is complex due to the use of measurements in model calibration. The models for elemental carbon are more recent and less thoroughly calibrated, potentially leading to misclassification of absolute concentrations. The EC maps are based on scaled emissions of PM2.5 and are calibrated using black smoke measurements that are converted to estimated EC values [35]. However, since a potential misclassification is likely consistent across The Maastricht Study population, comparisons of relative differences between groups remain valid.

We linked the estimated annual concentrations of air pollutants to The Maastricht Study participants’ home address at baseline, specifically for the year of the corneal observation. This approach ensures that the exposure assessment is contemporaneous with the corneal nerve assessment. For example: If corneal nerve assessment was conducted in 2015, we linked it to the air pollution levels of 2015.

Covariates

We determined glucose metabolism status (normal glucose metabolism, prediabetes, type 2 diabetes and other types of diabetes) according to the WHO 2006 criteria, based on a 75-g oral glucose tolerance test and use of glucose-lowering medication [36]. We assessed educational level (low, intermediate, high), occupational level (low, intermediate, high, self-employee, or other), household income level [37], smoking status (never, former, current), alcohol use (none, low, high), history of ocular disorders (corneal diseases or uveitis), and use of contact lenses by questionnaire. The Dutch Healthy Diet index [38] was based on a validated food frequency questionnaire [39]. Urbanity was determined by the average "environmental address density" within a center (summarized to neighborhoods) and classified as (very strong urban (> 2500 addresses per km2), strong urban (1500 – 2500 addresses per km2), moderately urban (1000 – 1500 addresses per km2), limited urban (500 – 1000 addresses per km2), non-urban (< 500 addresses per km2).

Statistical analyses

All statistical analyses were conducted using Statistical Package for Social Sciences version 25.0 (IBM SPSS, USA). A significance level of P < 0.05 was adopted for all analyses. Population characteristics and corneal nerve measurements were outlined for the entire study cohort and categorized by tertiles of PM2.5.

It is important to note that higher levels of air pollution are considered detrimental to health, and similarly, a decrease in corneal nerve measures is indicative of adverse outcomes.

Main analyses

We employed multivariable linear regression analyses to investigate the associations of Z-scores for PM2.5, PM10, NO2, and EC with Z-scores for corneal nerve bifurcation density, corneal nerve density, corneal nerve length, and fractal dimension.

Corneal confocal microscopy assessments were conducted from April 2013 to October 2020. Previous participants of the Maastricht Study who had been involved before the commencement of corneal confocal microscopy (i.e., prior to April 2013) were re-invited for corneal confocal microscopy examinations (n = 950; 26.1%). Among these participants, there was a median interval of 5.2 years, referred to as the 'visit interval' or 'lag time', between corneal confocal microscopy assessments and the initial baseline covariate measurements (such as glucose metabolism status). We incorporated adjustments for this lag time in our analysis.

We implemented various levels of adjustments in our analysis. Model 1 presents unadjusted results. Model 2 was adjusted for age, sex, level of education, glucose metabolism status, corneal confocal microscopy lag time (the interval between baseline measurements and corneal confocal microscopy scans), and inclusion year (presented as a categorical variable) of participants to control for potential period effects. Additionally, Model 3 included adjustments for smoking status and alcohol consumption. These variables were selected as they are significant potential confounders. All associations are presented as standardized regression coefficients along with their corresponding 95% confidence intervals (CI).

Additional analyses

To ensure the robustness of our findings, we performed several additional analyses. Firstly, we examined interaction effects with sex and glucose metabolism status by incorporating interaction terms into Model 3. Secondly, we repeated the analyses with additional adjustments for lifestyle factors (dietary intake excluding alcohol and physical activity), ocular variables (corneal diseases, uveitis), use of glasses or contact lenses, and urbanicity. These adjustments were not included in the primary analyses due to missing data for a substantial number of participants. Thirdly, we conducted additional analyses where adjustments were made for the location of captured corneal nerve images (encoded as dummy variables for center, semi-center, and inferior whorl), and the orientation of the corneal nerve fibers (treated as a continuous variable). This was done because the characteristics of corneal nerves vary based on the location of the captured image (e.g., corneal nerve density differs between the inferior whorl, center, and semi-center). Finally, we substituted level of education with either occupational status or income equivalent.

Results

Figure 2 shows an overview of selection of the study population. Table 1 and Supplementary Table S1 show the general participant characteristics according to tertiles of composite Z-score for PM2.5. The study population consisted of 3635 participants with an average age of 59.3 years (± 8.7), with 51.6% being women. The population was divided into tertiles based on PM2.5 concentration: Tertile 1 (high), Tertile 2 (middle), and Tertile 3 (low), each comprising approximately 1200 individuals.

Fig. 2
figure 2

Shows the selection of participants for inclusion in the analyses

Table 1 General study population characteristics in the study population with complete data on PM2.5 and according to tertiles of PM2.5

The general characteristics of participants included in the study were found to be akin to those with missing data (see Supplementary Table S2).

The associations between air pollutants and corneal nerve measures are shown in Table 2 and Fig. 3. According to model 3, higher Z-scores for PM2.5 and PM10 were associated with lower Z-scores for corneal nerve bifurcation density, nerve density, nerve length, and nerve fractal dimension [stβ (95% CI): PM2.5 -0.10 (-0.14; -0.05), -0.04 (-0.09; 0.01), -0.11 (-0.16; -0.06), -0.20 (-0.24; -0.15); and PM10 -0.08 (-0.13; -0.03), -0.04 (-0.09; 0.01), -0.08 (-0.13; -0.04), -0.17 (-0.21; -0.12)], respectively.

Table 2 Associations of Z-scores for PM2.5, PM10, NO2, and EC with Z-scores for corneal nerve bifurcation density, corneal nerve density, corneal nerve length, and corneal nerve fractal dimension
Fig. 3
figure 3

Associations of Z-score for PM2.5, PM10, NO2, and EC with Z-score for bifurcation density, corneal nerve density, corneal nerve length, and corneal nerve fractal dimension. Standardized regression coefficients (stβ) represent the differences in corneal nerve measures in SD per SD greater measure of air pollutants. For PM2.5, PM10, and EC, 1 SD corresponds to 39.9 number of branches/mm2 for corneal nerve bifurcation density, 24.3 number of main and branches/mm2 for corneal nerve density, 4.4 mm/mm2 for corneal nerve length, 0.1 (unit-less) for corneal nerve fractal dimension, 1.4 µg/m3 for PM2.5, 1.7 µg/m3 for PM10, 3.1 µg/m3 for NO2, and 0.1 µg/m3 for EC. Variables entered into model: age, sex, level of education (high, medium, low), glucose metabolism status (prediabetes and type 2 diabetes versus normal glucose metabolism), corneal confocal microscopy lag time, inclusion year of participants, smoking status (never, former, current), and alcohol consumption status (none, low, high). Asterisks indicate values that are statistically significant (P < 0.05). Abbreviations: CI, confidence interval; SD, standard deviation

No statistically significant associations were found between Z-scores for NO2 and EC with Z-scores for corneal nerve bifurcation density, nerve density, nerve length, and nerve fractal dimension [stβ (95% CI): PM2.5 -0.01 (-0.04; 0.03), 0.01 (-0.03; 0.04), -0.02 (-0.05; 0.01), -0.02 (-0.05; 0.02); and PM10 -0.01 (-0.04; 0.03), 0.01 (-0.03; 0.04), -0.02 (-0.05; 0.02), -0.001 (-0.03; 0.03)], respectively; Fig. 3 and Table 2.

Sex and glucose metabolism status did not modify the associations of Z-scores for PM2.5, PM10, NO2, and EC with composite Z score for corneal nerve measures [p value for the interaction term (potential determinant × sex; potential determinant × prediabetes versus normal glucose metabolism; potential determinant × type 2 diabetes versus normal glucose metabolism) were all < 0.05); Supplementary Table S3 and S4.

Quantitatively similar results were observed in a range of sensitivity analyses. First, associations were generally comparable to the main results after additional adjustment for dietary intake and physical activity (model 4), for ocular variables (corneal diseases, uveitis, model 5), use of glasses or contact lenses (model 6), and urbanicity (model 7) (Supplemental Table S5). Second, associations were not altered after additional adjustment for location of captured corneal nerve images (model 8), or for the orientation of the corneal nerve fibers (model 9) (Supplemental Table S6). Last, associations remained similar after replacement of level of education with occupational status (model 10) or income level (model 11) (Supplemental Table S7).

Discussion

The primary findings of this cross-sectional observational study conducted at the population level are as follows. First, after adjustment for a range of potential confounders, higher levels of PM2.5 and PM10 were associated with lower corneal nerve measures. Second, NO2 and EC were not associated with corneal nerve measures. Third, the associations were not significantly different between men and women or between different glucose metabolism strata. Further, our results remained similar after additional adjustment for lifestyle factors, corneal diseases, uveitis, use of glasses or contact lenses, and urbanicity.

To the best of our knowledge, our study represents the inaugural large-scale population-based investigation unveiling novel associations between exposure to air pollutants and corneal nerve morphology. This unique association highlights a potential relationship between specific air pollutant concentrations and neurodegeneration. No earlier study has included corneal nerve measures as an outcome to investigate ocular surface disease or as a marker for neurodegeneration.

Gayraud et al. previously revealed an association between air pollutants and retinal nerve fiber layer thickness in the Alienor study population. In this study of 683 participants, they showed that higher exposure to PM2.5 and black carbon for 10 years was significantly associated with a faster retinal nerve fiber layer thinning during the 11-year follow-up. Moreover, no statistically significant association was found with NO2 [40], in line with our findings. Our findings also align with other research regarding the impact of air pollution on the ocular surface [3]. Studies have demonstrated the harmful effects of exposure to air pollutants, which may cause irritation and inflammation, resulting in conjunctivitis [3]. Air pollution may also cause tear film instability inducing the expression of inflammatory mediators in the tear film leading to dry eye disease [3]. Jurkiewicz et al. also proposed that fine particulate matter may be a risk factor for keratoconus [9].

Our results are consistent with studies demonstrating associations between air pollution exposure and neurodegeneration [41]. A systematic review and meta-analysis demonstrated a significant association between PM2.5 exposure and lower general cognition, verbal fluency, and executive function in individuals aged 40 and above [42]. Furthermore, Sakhvidi et al. showed that exposure to PM2.5 was associated with lower cognitive function in the French CONSTANCES cohort, n = 220,000 people (aged 18–69 years) [43]. Moreover in a recent longitudinal study, late life exposure to PM2.5 and NO2 was associated with Alzheimer's disease-related neurodegeneration, specifically, medial temporal lobe atrophy [44]. Indeed, recent studies have shown an association between corneal nerve loss and cognitive function in patients with mild cognitive impairment and dementia [15, 45, 46] and corneal nerve loss in subjects with mild cognitive impairment predicted the development of dementia [47]. Thus the corneal nerve loss observed could act as a surrogate marker for neurodegeneration in the brain [48].

The significant associations observed between PM2.5 and PM10 and corneal neurodegeneration, in contrast to the absence of associations with NO2 and EC, may be attributed to distinct mechanisms by which particulate matter affects corneal nerves. There is a strong body of data showing that different components of air pollutants, including particulate matter, carbon and nitrogen oxides have a varying impact on the nervous system and neurological disorders [49]. Longitudinal studies evaluating alterations over time are needed to provide insights into the temporal aspects of these differing associations.

As previously detailed by our research team [20] the corneal nerve metrics observed in our study population diverge from those in previous studies [50] due to several factors: (1) the assessment of corneal nerve variables in larger images, spanning an area of up to 1600 × 1600 µm, which partly encompassed the inferior whorl in certain composite images, characterized by distinct architecture, distribution, and density of corneal nerve fibers [51]; (2) the utilization of a different deep learning model (U-Net-based convolutional neural network) for fully automated tracing and analysis of corneal nerves [52]; nd (3) the definition of corneal nerve density and length differing from other studies (total number and length of corneal nerves, encompassing both main and branches, per mm2).

Our study possesses several strengths. First, to the best of our knowledge, it represents the first extensive population-based observational investigation seeking to establish the association between different air pollutants and corneal nerve morphology. Second, we adjusted for a broad range of confounders and our additional analyses generally yielded consistent findings. Third, we computed z-scores for air pollution and corneal nerve measures (continuous variables) to enable comparisons of associations between determinants and outcomes across various characteristics and units. Lastly, we revealed consistent associations with a range of different corneal nerve measures (bifurcation density, nerve density, nerve length, and nerve fractal dimension).

Our study also had certain limitations. First, due to its observational cross-sectional design, caution should be exercised in drawing any causal conclusions. However, the likelihood of reverse causation is deemed low, and we meticulously adjusted for numerous potential confounders. Secondly, despite accounting for a wide range of confounders, we cannot entirely rule out bias stemming from unmeasured confounding factors (such as workplace environment and indoor air pollution). Third, we were not able to account for change in address over time, hence we may have spuriously estimated ambient air pollution levels for those individuals who moved, potentially leading to an underestimation of associations under study. Fourth, there is a possibility that we underestimated associations due to the restricted range of distribution, as our analysis only encompassed participants residing in the Maastricht area. since corneal confocal microscopy measurements commenced in April 2013, individuals who engaged in the Maastricht Study before the initiation of corneal confocal microscopy assessment (i.e., before April 2013; n = 950; 26.1%) had a median visit interval of 5.2 years [4.9–6.1]. Although we adjusted for visit interval in our analyses to yield a more accurate estimate of the true association. Finally, this study encompassed individuals of European descent aged between 40 and 75 years residing in the Maastricht area. Further investigations are warranted to ascertain the generalizability of our findings to diverse populations and regions.

Our findings may have implications for public health policy. The significant association between exposure to PM2.5 and PM10 and corneal neurodegeneration suggests a potential need to consider enhancing air quality standards in accordance with the current European Union (EU) guidelines. Notably, our associations persist even at pollution levels below the EU's recommended thresholds, which could support discussions around revising these standards, as suggested by the WHO in 2021. In Maastricht, where our study was conducted, average concentrations of PM2.5 and PM10 (12.2 and 19.0 μg/m3, respectively) surpass WHO-recommended levels (5 and 15 μg/m3, respectively) (1) but fall below EU thresholds (25 and 40 μg/m3, respectively) (European Commission, 2008).

In conclusion, the present population-based study demonstrated that exposure to higher levels of PM2.5 and PM10 are associated with higher levels of corneal neurodegeneration, estimated from lower corneal nerve measures. Our results suggest that air pollution may be a determinant to corneal neurodegeneration and impact the ocular surface health.

Availability of data and materials

Data from the Maastricht Study are available to any researcher who meets the criteria for access to confidential data https://demaastrichtstudie.app/data-dictionary/. The corresponding author may be contacted to request access. Environmental exposure data can be obtained via www.gecco.nl.

Data Availability

No datasets were generated or analysed during the current study.

Abbreviations

PM2.5 :

Particulate Matter 2.5

PM10 :

Particulate Matter 10

NO2 :

Nitrogen Dioxide

CI:

Confidence interval

SD:

Standard deviation

References

  1. Taylor L. WHO cuts air pollution limits to save millions of lives. BMJ. 2021;374:n2349.

  2. Manisalidis I, Stavropoulou E, Stavropoulos A, Bezirtzoglou E. Environmental and health impacts of air pollution: a review. Front Public Health. 2020;8:14.

    Article  Google Scholar 

  3. Lin C-C, Chiu C-C, Lee P-Y, Chen K-J, He C-X, Hsu S-K, Cheng K-C. The adverse effects of air pollution on the eye: a review. Int J Environ Res Public Health. 2022;19(3):1186.

    Article  CAS  Google Scholar 

  4. You R, Ho Y-S, Chang RC-C. The pathogenic effects of particulate matter on neurodegeneration: a review. J Biomed Sci. 2022;29(1):15.

    Article  Google Scholar 

  5. Saxena R, Srivastava S, Trivedi D, Anand E, Joshi S, Gupta SK. Impact of environmental pollution on the eye. Acta Ophthalmol Scand. 2003;81(5):491–4.

    Article  CAS  Google Scholar 

  6. Mo Z, Fu Q, Lyu D, Zhang L, Qin Z, Tang Q, et al. Impacts of air pollution on dry eye disease among residents in Hangzhou, China: a case-crossover study. Environ Pollut. 2019;246:183–9.

    Article  CAS  Google Scholar 

  7. Mimura T, Ichinose T, Yamagami S, Fujishima H, Kamei Y, Goto M, et al. Airborne particulate matter (PM2. 5) and the prevalence of allergic conjunctivitis in Japan. Sci Total Environ. 2014;487:493–9.

    Article  CAS  Google Scholar 

  8. Hong J, Zhong T, Li H, Xu J, Ye X, Mu Z, et al. Ambient air pollution, weather changes and outpatient visits for allergic conjunctivitis: A retrospective registry study. Sci Rep. 2016;6(1):23858.

    Article  CAS  Google Scholar 

  9. Jurkiewicz T, Marty A-S. Correlation between Keratoconus and Pollution. Ophthalmic Epidemiol. 2021;28(6):495–501.

    Article  CAS  Google Scholar 

  10. Yu D, Deng Q, Wang J, Chang X, Wang S, Yang R, et al. Air Pollutants are associated with dry eye disease in urban ophthalmic outpatients: a prevalence study in China. J Transl Med. 2019;17(1):46.

    Article  Google Scholar 

  11. Hao R, Zhang M, Zhao L, Liu Y, Sun M, Dong J, et al. Impact of Air Pollution on the Ocular Surface and Tear Cytokine Levels: A Multicenter Prospective Cohort Study. Front Med (Lausanne). 2022;9:909330.

  12. Kashiwagi K, Iizuka Y. Effect and underlying mechanisms of airborne particulate matter 2.5 (PM2.5) on cultured human corneal epithelial cells. Sci Rep. 2020;10(1):19516.

    Article  CAS  Google Scholar 

  13. Alhussaini AR, Aljabri MR, Al-Harbi ZT, Abdulrahman Almohammadi G, Al-Harbi TM, Bashir S. Air pollution and its adverse effects on the central nervous system. Cureus. 2023;15(5):e38927.

    Google Scholar 

  14. Shaheen BS, Bakir M, Jain S. Corneal nerves in health and disease. Surv Ophthalmol. 2014;59(3):263–85.

    Article  Google Scholar 

  15. Ponirakis G, Al Hamad H, Sankaranarayanan A, Khan A, Chandran M, Ramadan M, et al. Association of corneal nerve fiber measures with cognitive function in dementia. Ann Clin Transl Neurol. 2019;6(4):689–97.

    Article  Google Scholar 

  16. Misra SL, Kersten HM, Roxburgh RH, Danesh-Meyer HV, McGhee CN. Corneal nerve microstructure in Parkinson’s disease. J Clin Neurosci. 2017;39:53–8.

    Article  Google Scholar 

  17. Ferrari G, Grisan E, Scarpa F, Fazio R, Comola M, Quattrini A, et al. Corneal confocal microscopy reveals trigeminal small sensory fiber neuropathy in amyotrophic lateral sclerosis. Front Aging Neurosci. 2014;6:278.

    Article  Google Scholar 

  18. Bitirgen G, Akpinar Z, Malik RA, Ozkagnici A. Use of corneal confocal microscopy to detect corneal nerve loss and increased dendritic cells in patients with multiple sclerosis. JAMA Ophthalmol. 2017;135(7):777–82.

    Article  Google Scholar 

  19. So WZ, Qi Wong NS, Tan HC, Yu Lin MT, Yu Lee IX, Mehta JS, Liu YC. Diabetic corneal neuropathy as a surrogate marker for diabetic peripheral neuropathy. Neural Regen Res. 2022;17(10):2172–8.

    Article  CAS  Google Scholar 

  20. Mokhtar SBA, van der Heide FCT, Oyaert KAM, van der Kallen CJH, Berendschot T, Scarpa F, et al. (Pre)diabetes and a higher level of glycaemic measures are continuously associated with corneal neurodegeneration assessed by corneal confocal microscopy: the Maastricht Study. Diabetologia. 2023;66(11):2030–41.

    Article  CAS  Google Scholar 

  21. Vereertbrugghen A, Galletti JG. Corneal nerves and their role in dry eye pathophysiology. Exp Eye Res. 2022;222:109191.

    Article  CAS  Google Scholar 

  22. Tidke SC, Tidake P. A review of corneal blindness: causes and management. Cureus. 2022;14(10):e30097.

    Google Scholar 

  23. Haffner SM, Stern MP, Hazuda HP, Mitchell BD, Patterson JK. Cardiovascular risk factors in confirmed prediabetic individuals: does the clock for coronary heart disease start ticking before the onset of clinical diabetes? JAMA. 1990;263(21):2893–8.

    Article  CAS  Google Scholar 

  24. Wong M-S, Gu K, Heng D, Chew S-K, Chew L-S, Tai ES. The Singapore impaired glucose tolerance follow-up study: does the ticking clock go backward as well as forward? Diabetes Care. 2003;26(11):3024–30.

    Article  Google Scholar 

  25. Allgeier S, Zhivov A, Eberle F, Koehler B, Maier S, Bretthauer G, et al. Image reconstruction of the subbasal nerve plexus with in vivo confocal microscopy. Invest Ophthalmol Vis Sci. 2011;52(9):5022–8.

    Article  Google Scholar 

  26. Schram MT, Sep SJ, van der Kallen CJ, Dagnelie PC, Koster A, Schaper N, et al. The Maastricht study: an extensive phenotyping study on determinants of type 2 diabetes, its complications and its comorbidities. Eur J Epidemiol. 2014;29(6):439–51.

    Article  Google Scholar 

  27. Timmermans EJ, Lakerveld J, Beulens JWJ, Boomsma DI, Kramer SE, Oosterman M, et al. Cohort profile: the Geoscience and Health Cohort Consortium (GECCO) in the Netherlands. BMJ Open. 2018;8(6):e021597.

    Article  Google Scholar 

  28. Scarpa F, Colonna A, Ruggeri A. Multiple-image deep learning analysis for neuropathy detection in corneal nerve images. Cornea. 2020;39(3):342–7.

    Article  Google Scholar 

  29. Organization WH. WHO global air quality guidelines: particulate matter (PM2.5 and PM10), ozone, nitrogen dioxide, sulfur dioxide and carbon monoxide. Geneva: World Health Organization; 2021. Available from: https://www.thelancet.com/action/showPdf?pii=S2542-5196%2823%2900258-9.

  30. Hoogerbrugg R HS, Huitema M, Siteur K, Smeets W, Soenario I, Visser S, de Vries WJ, Wichink Kruit RJ. Large-scale concentration maps for the Netherlands. 2023 report. the Netherlands: RIVM; 2023. Contract No.: 2023–0113.

  31. Wesseling J BR, Bezemer A. An efficient modeling system for nation-wide compliance testing. 14th Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes; Kos. 2011.

  32. Lagali NS, Allgeier S, Guimarães P, Badian RA, Ruggeri A, Köhler B, et al. Reduced corneal nerve fiber density in type 2 diabetes by wide-area mosaic analysis. Invest Ophthalmol Vis Sci. 2017;58(14):6318–27.

    Article  CAS  Google Scholar 

  33. Atlas EH. Explore and discover your living environment: Environmental Health Atlas. Available from: https://www.atlasleefomgeving.nl/en.

  34. Wesseling JNL, Hoogerbrugge R. Gemeten en berekende concentraties stikstof(di)oxiden en fijnstof in de periode 2010 t/m 2015 (Update) : Een test van de standaardrekenmethoden 1 en 2. the Netherlands: RIVM; 2018 2018–12–19. Contract No.: 2016–0106.

  35. Wesseling JPvdS, E.; Hoogerbrugge R.; Velders, G.J.M.; de Vries, W.J.; van den Elshout, S.; Ameling, M. Elemental Carbon (EC) concentrations in the Netherlands. the Netherlands: RIVM; 2015. Contract No.: 2014–0128.

  36. Alberti KGMM, Zimmet PZ. Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus. Provisional report of a WHO consultation. Diabetic Med. 1998;15(7):539–53.

    Article  CAS  Google Scholar 

  37. Qi Y, Koster A, van Boxtel M, Köhler S, Schram M, Schaper N, Stehouwer C, Bosma H. Adulthood Socioeconomic Position and Type 2 Diabetes Mellitus—A Comparison of Education, Occupation, Income, and Material Deprivation: The Maastricht Study. Int J Environ Res Public Health. 2019;16(8):1435. https://doi.org/10.3390/ijerph16081435.

  38. Looman M, Feskens EJ, de Rijk M, Meijboom S, Biesbroek S, Temme EH, et al. Development and evaluation of the Dutch healthy diet index 2015. Public Health Nutr. 2017;20(13):2289–99.

    Article  Google Scholar 

  39. van Dongen MC, Wijckmans-Duysens NEG, den Biggelaar LJ, Ocké MC, Meijboom S, Brants HA, et al. The Maastricht FFQ: Development and validation of a comprehensive food frequency questionnaire for the Maastricht study. Nutrition. 2019;62:39–46.

    Article  Google Scholar 

  40. Gayraud L, Mortamais M, Schweitzer C, de Hoogh K, Cougnard-Grégoire A, Korobelnik JF, et al. Association of long-term exposure to ambient air pollution with retinal neurodegeneration: the prospective Alienor study. Environ Res. 2023;232:116364.

    Article  CAS  Google Scholar 

  41. Calderón-Garcidueñas L, Leray E, Heydarpour P, Torres-Jardón R, Reis J. Air pollution, a rising environmental risk factor for cognition, neuroinflammation and neurodegeneration: The clinical impact on children and beyond. Revue Neurologique. 2016;172(1):69–80.

    Article  Google Scholar 

  42. Thompson R, Smith RB, Karim YB, Shen C, Drummond K, Teng C, Toledano MB. Air pollution and human cognition: a systematic review and meta-analysis. Sci Total Environ. 2023;859:160234.

    Article  CAS  Google Scholar 

  43. Zare Sakhvidi MJ, Yang J, Lequy E, Chen J, de Hoogh K, Letellier N, et al. Outdoor air pollution exposure and cognitive performance: findings from the enrolment phase of the CONSTANCES cohort. Lancet Planet Health. 2022;6(3):e219–29.

    Article  Google Scholar 

  44. Petkus AJ, Salminen L, Wang X, Driscoll I, Millstein J, Beavers D, et al. Alzheimer's related neurodegeneration mediates air pollution effects on medial temporal lobe atrophy. Medrxiv. 2023:2023.11. 29.23299144.

  45. Al-Janahi E, Ponirakis G, Al Hamad H, Vattoth S, Elsotouhy A, Petropoulos IN, et al. Corneal nerve and brain imaging in mild cognitive impairment and dementia. J Alzheimers Dis. 2020;77(4):1533–43.

    Article  CAS  Google Scholar 

  46. Ponirakis G, Hamad HA, Khan A, Petropoulos IN, Gad H, Chandran M, et al. Loss of corneal nerves and brain volume in mild cognitive impairment and dementia. Alzheimer’s Dement: Transl Res Clin Intervent. 2022;8(1):e12269.

    Article  Google Scholar 

  47. Ponirakis G, Al Hamad H, Omar DA, Petropoulos IN, Khan A, Gad H, et al. Corneal nerve loss predicts dementia in patients with mild cognitive impairment. Ann Clin Transl Neurol. 2023;10(4):599–609.

    Article  CAS  Google Scholar 

  48. Chaitanuwong P, Singhanetr P, Chainakul M, Arjkongharn N, Ruamviboonsuk P, Grzybowski A. Potential ocular biomarkers for early detection of Alzheimer’s disease and their roles in artificial intelligence studies. Neurol Ther. 2023;12(5):1517–32.

    Article  Google Scholar 

  49. Fu C, Kuang D, Zhang H, Ren J, Chen J. Different components of air pollutants and neurological disorders. Front Public Health. 2022;10:959921.

    Article  Google Scholar 

  50. Gad H, Petropoulos IN, Khan A, Ponirakis G, MacDonald R, Alam U, Malik RA. Corneal confocal microscopy for the diagnosis of diabetic peripheral neuropathy: A systematic review and meta-analysis. J Diabetes Investig. 2022;13(1):134–47.

    Article  Google Scholar 

  51. Patel DV, McGhee CN. Mapping of the normal human corneal sub-Basal nerve plexus by in vivo laser scanning confocal microscopy. Invest Ophthalmol Vis Sci. 2005;46(12):4485–8.

    Article  Google Scholar 

  52. Colonna A, Scarpa F, Ruggeri A. Segmentation of corneal nerves using a U-net-based convolutional neural network: first international workshop, COMPAY 2018, and 5th International Workshop, OMIA 2018, Held in Conjunction with MICCAI 2018, Granada, Spain, September 16 - 20, 2018, Proceedings. 2018. p. 185-92.

Download references

Acknowledgements

The authors gratefully acknowledge the ZIO Foundation (Vereniging Regionale HuisartsenZorg Heuvelland) or their support of the Maastricht Study. We are indebted to the participants for their invaluable contribution and willingness to take part in this study.

Web references

European Commission. "EU air quality standards." 2008. Accessed March 21, 2024. https://environment.ec.europa.eu/topics/air/air-quality/eu-air-quality-standards_en.

Funding

This study was supported by the European Regional Development Fund via OP-Zuid, the Province of Limburg, the Dutch Ministry of Economic Affairs (grant 31O.041), Stichting de Weijerhorst (Maastricht, the Netherlands), the Pearl String Initiative Diabetes (Amsterdam, the Netherlands), the Cardiovascular Center (Maastricht, the Netherlands), the School for Cardiovascular Diseases (Maastricht, the Netherlands), the School of Mental Health and Neuroscience (Maastricht, the Netherlands), the School for Public Health and Primary Care (Maastricht, the Netherlands), the School for Nutrition and Translational Research in Metabolism (Maastricht, the Netherlands), Stichting Annadal (Maastricht, the Netherlands), Health Foundation Limburg (Maastricht, the Netherlands), Perimed (Järfälla, Sweden) and unrestricted grants from Janssen-Cilag (Tilburg, the Netherlands), Novo Nordisk (Alphen aan den Rijn, the Netherlands) and Sanofi-Aventis(Gouda, the Netherlands). GECCO is financially supported by The Netherlands Organisation for Scientific Research (NWO), The Netherlands Organisation for Health Research and Development (ZonMw), and Amsterdam UMC. More information on GECCO can be found on www.GECCO.nl.

Author information

Authors and Affiliations

Authors

Contributions

SBAM contributed to the conception and design of this study, participated in acquisition of data, analysed and interpreted data, drafted the manuscript (with AK, CDAS, CABW, CJHvdK and TTJMB), revised the manuscript critically for important intellectual content, and provided final approval of the version to be published. SBAM is also the guarantor of this work, and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analyses. JV, JA, JS, FS, AC, FCTvdH, MMJvG, HB, SK, NCS, AAK, RMMAN, MG, JL, and RM contributed to the conception and design of the study, revised the manuscript critically for important intellectual content, and provided final approval of the version to be published.

Corresponding author

Correspondence to Sara B. A. Mokhtar.

Ethics declarations

Ethics approval and consent to participate

The Maastricht study was approved by the institutional medical ethical committee (NL31329.068.10) and the Minister of Health, Welfare and Sports of the Netherlands (permit 131088–105234-PG). All participants gave written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mokhtar, S.B.A., Viljoen, J., van der Kallen, C.J.H. et al. Greater exposure to PM2.5 and PM10 was associated with lower corneal nerve measures: the Maastricht study - a cross-sectional study. Environ Health 23, 70 (2024). https://doi.org/10.1186/s12940-024-01110-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12940-024-01110-1

Keywords