Mean 1-year black carbon concentration (A) and mean intraocular pressure (B) across allelic risk scores for endothelial function, metal processing, and oxidative stress pathways. Error bars show SEs.
Difference in mean intraocular pressure for a 1–interquartile range (IQR) increase (1 IQR = 0.22 μg/m3) in 1 year of black carbon exposure was compared between participants with low pathway-specific allelic risk scores and those with high scores in 3 model frameworks. Tier 3 indicates the main fully adjusted mixed-effects model; tier 4, sensitivity analysis model using the tier 3 model adjusted for glaucoma disease status; and medication, sensitivity analysis using the tier 3 model adjusted for glaucoma eyedrops and hypertension disease status replaced with antihypertension medication. P values were corrected for false discovery rate.
Customize your JAMA Network experience by selecting one or more topics from the list below.
Nwanaji-Enwerem JC, Wang W, Nwanaji-Enwerem O, et al. Association of Long-term Ambient Black Carbon Exposure and Oxidative Stress Allelic Variants With Intraocular Pressure in Older Men. JAMA Ophthalmol. 2019;137(2):129–137. doi:10.1001/jamaophthalmol.2018.5313
What is the association between long-term ambient black carbon exposure and intraocular pressure?
In this analysis of 911 follow-up visits by 419 older men enrolled in a longitudinal cohort study, the association of long-term black carbon exposure with intraocular pressure was greater in individuals with a high oxidative stress allelic score when compared with individuals with a low score.
Ambient black carbon exposure may be a risk factor for increased intraocular pressure, particularly in individuals susceptible to other biological oxidative stressors.
Elevated intraocular pressure is a major risk factor for glaucoma, a leading cause of irreversible blindness worldwide. Environmental air pollution has been suggested as a potential contributor to elevated intraocular pressure; however, no studies have demonstrated such an association to date.
To investigate the association of long-term ambient black carbon exposure with intraocular pressure in community-dwelling older adults.
Design, Setting, and Participants
This population-based analysis, conducted from October 18, 2017, through March 22, 2018, used data from the all-male, New England–based Normative Aging Study of the US Department of Veterans Affairs. The analysis included 419 older men with a total of 911 follow-up study visits between January 1, 2000, and December 30, 2011. Intraocular pressure was measured by Goldmann applanation tonometry during the study visits. Validated spatiotemporal models were used to generate 1-year black carbon exposure levels at the addresses of the participants.
Main Outcomes and Measures
An independently developed genetic score approach was used to calculate allelic risk scores for 3 pathways associated with black carbon toxicity: endothelial function, oxidative stress, and metal processing. The associations among black carbon exposure, allelic risk scores, and intraocular pressure were explored using linear mixed-effects models.
All 419 participants were men with a mean (SD) age of 75.3 (6.9) years. The mean (SD) 1-year black carbon exposure was 0.51 (0.18) μg/m3, and the mean (SD) intraocular pressure for the left eye was 14.1 (2.8) mm Hg and for the right eye was 14.1 (3.0) mm Hg. Of the 911 visits, 520 (57.1%) had a high endothelial function allelic risk score, 644 (70.7%) had a high metal-processing allelic risk score, and 623 (68.4%) had a high oxidative stress allelic risk score. In fully adjusted linear mixed-effects models, the association of black carbon with intraocular pressure was greater in individuals with a high oxidative stress allelic score (β = 0.36; 95% CI, 0.003-0.73) compared with individuals with a low score (β = −0.35; 95% CI, −0.86 to 0.15).
Conclusions and Relevance
Ambient black carbon exposure may be a risk factor for increased intraocular pressure in individuals susceptible to other biological oxidative stressors. If additional studies confirm these results, monitoring ambient black carbon exposure and physiological oxidative stress may prevent the development and progression of intraocular pressure–related disease.
Black carbon (BC), a byproduct of combustion processes, is among the most pervasive ambient particles in the world.1 Exposure to BC has a well-documented association with numerous adverse health outcomes, but the association of BC with increased systemic blood pressure (BP) is among the most extensively studied.2-4 Existing studies both replicate the BC-BP association and provide evidence suggesting that biological processes, such as endothelial function, metal processing, and oxidative stress, may explain elements of the association.5-8
Intraocular pressure (IOP) is a major modifiable risk factor for ocular diseases, including glaucoma, a leading cause of blindness in the world.9,10 Evidence suggests that elevated systemic BP and elevated IOP share some physiologic mechanisms.11-15 To our knowledge, studies of the association of ambient BC with IOP in a population-based sample have not yet been performed. Understanding this association might provide another reason for reducing the contribution of BC to environmental pollution, especially if subsequent studies add adequate strength to such an association.16,17
We hypothesized that BC exposure is associated with increased IOP and that biological processes involved in the BC-BP association are also implicated in the BC-IOP association. The latter hypothesis was evaluated by testing if genetic variants in 3 biological pathways implicated in ambient particle–related disease (endothelial function, metal processing, and oxidative stress) could modify the association of BC with IOP.5 The present study was conducted from October 18, 2017, to March 22, 2018. We tested these hypotheses by using longitudinal data from the all-male cohort of the US Department of Veterans Affairs (VA) Normative Aging Study (NAS). Since its establishment in 1963, the NAS has collected a multitude of health-related variables throughout participants’ lifetimes with the aim of conducting analyses to better understand the aging process.18 The present study is in line with that NAS aim.
The NAS is an ongoing longitudinal study of aging in male volunteers. In 1963, it began to recruit male participants from the Greater Boston area who were free of any chronic disease.18 Although now a closed cohort, the NAS participants return for follow-up study visits every 3 to 5 years since recruitment. These on-site, follow-up visits involve comprehensive physical examinations and detailed collection of data on lifestyle factors that may affect health, including physical activity, diet, smoking habits, alcohol intake, and medications. In this cohort, dropout has been less than 1% per year and predominantly occurs when participants move out of the study area. In the present study, the sample (419 men with 911 study visits) consisted of NAS participants who had all available exposure, outcome, sociodemographic, genetic, and clinical data. Of these 419 participants, 124 (29.6%) had 1 study visit, 135 (32.2%) had 2 study visits, 123 (29.4%) had 3 study visits, and 37 (8.8%) had 4 study visits between January 1, 2000, and December 30, 2011. At recruitment, all NAS participants provided written informed consent to the VA Institutional Review Board and were at least aged 18 years; this initial consent covered subsequent studies such as ours. For this current study, human participant approval was granted by the VA and Harvard T. H. Chan School of Public Health Institutional Review Boards.
In the morning of each study visit, predilated IOP of right and left eyes was measured by an ophthalmologist using Goldmann applanation tonometry (the current criterion standard).19 Measurements could range from 0 to 30 mm Hg, but our study sample had a range of 5 to 28 mm Hg. For each participant, we calculated the IOP between both eyes to obtain a mean IOP.
We focused on long-term (1-year) BC exposure because it has previously been associated with BP and has consistently been associated with other aging-related diseases.20-23 We began by using a validated spatiotemporal land-use regression model to generate daily BC exposure estimates (in micrograms per cubic meter) at the 1 × 1-km area resolution on the basis of participants’ residences.24 The spatiotemporal model was based on information from daily BC concentrations at a central monitor, daily BC concentrations at 83 monitoring sites throughout the Greater Boston area, land use (eg, traffic density), meteorological conditions (eg, wind speed), and other descriptors (eg, day of the week). In the training data set, the prediction model had a high multivariate coefficient of determination (R2) of 0.83, and the mean correlation between estimated values and observed BC levels in 4 out-of-sample validation samples was moderate (R2 = 0.59). Mean daily exposures for the 365 days prior to the day of each participant’s NAS visit were calculated to obtain the 1-year exposure estimate.
We used genotyping data from the NAS data set and a novel allelic risk score method that was independently developed by Bind et al5 to investigate interactions between environmental exposures and the biological pathways of endothelial function, metal processing, and oxidative stress. Bind et al5 developed the scores using the known biological functionality of genetic polymorphisms25 and independent outcomes representative of each pathway (augmentation index for endothelial function, patella lead concentration for metal processing, and 8-hydroxy-2′-deoxyguanosine for oxidative stress). Each score represents a pathway-specific allelic risk profile, and scores are dichotomized as high or low allelic risk according to score distributions in the study sample of interest (ie, high: ≥median; low: <median). Further details on the allelic scores can be found in the original publication.5
We used linear mixed-effects models with data from all patient visits to determine the association of 1-year BC levels with mean IOP. All covariates were determined a priori on the basis of existing literature regarding BC and IOP.3,26-28 Next, we used a tiered framework to build the final model. Tier 1 adjusted for chronological age (continuous) and year of visit (continuous). Tier 2 made additional adjustments for mean 1-year temperature (continuous address-specific satellite measurements29), cumulative cigarette pack-years (continuous), smoking status (current, former, or never), season of visit (spring [March-May], summer [June-August], fall [September-November], and winter [December-February]), body mass index (calculated as weight in kilograms divided by height in meters squared; lean [<25], overweight [25-30], or obese [>30]), alcohol intake (yes or no; ≥2 drinks daily), systolic BP (continuous), diastolic BP (continuous), fasting blood glucose level (continuous), C-reactive protein level (continuous), and total serum cholesterol level (continuous). Tier 3 paralleled tier 2 but replaced fasting blood glucose level, systolic and diastolic BP, and cholesterol level with the dichotomized disease variables of diabetes, hypertension, and ischemic heart disease. To examine the robustness of tier 3, we created a fourth tier of covariates. Tier 4 used the tier 3 covariates but made an additional adjustment for glaucoma defined as follows:
either eye’s cup-disc ratio of 0.7 or higher;
asymmetric cup-disc ratio (the difference of 2 eyes’ cup-disc ratio ≥0.2);
any eye’s cup-disc ratio of 0.6 or higher, with either disc hemorrhage or visual field defect; or
vision loss resulting from nerve fiber layer loss.
Tier 3 resulted in the best model fit (lowest Akaike information criterion) and was used as the final fully adjusted model.
For descriptive purposes, we ran a preliminary model of BC and the 3 allelic risk scores as joint predictors of mean IOP using the fully adjusted framework. We tested whether any of the allelic risk scores substantially modified the association of BC with mean IOP. False discovery rate correction was conducted with the Benjamini-Hochberg procedure.30 To further examine the robustness of our effect modification findings, we ran sensitivity analyses using the tier 4 covariates without alterations and then with the glaucoma and hypertension disease status replaced with glaucoma eyedrops and any antihypertensive medication.
All statistical analyses were performed using R, version 3.4.1 (R Core Team). Linear mixed-effects models were run using the lme function of the nlme R package31 and included a random participant-specific intercept to account for the correlation between repeated outcome measures (ie, multiple visits for a participant). The lme function uses a t distribution/statistic to calculate a 2-sided P value. P < .05 was considered to be statistically significant.
The demographic and clinical characteristics of participants across all study visits are presented in Table 1. All 419 participants were men, with a mean (SD) age of 75.3 (6.9) years. The mean (SD; interquartile range [IQR]) 1-year BC exposure was 0.51 (0.18; 0.22) μg/m3, and the mean (SD) IOP for the left eye was 14.1 (2.8) mm Hg and for the right eye was 14.1 (3.0) mm Hg. Most participants had completed at least 12 years of formal education (as indicated at 674 visits [74.0%]), were overweight or obese (711 visits [78.0%]), and were former smokers (595 visits [65.3%]). In this study sample, the prevalence of diabetes was 183 participant visits (20.1%); glaucoma, 86 (9.4%); hypertension; 700 (76.8%), and ischemic heart disease, 297 (32.6%).
Among the 911 visits, 520 (57.1%) showed evidence that the participant had a high endothelial function allelic risk score; 644 (70.7%), a high metal-processing allelic risk score; and 623 (68.4%), a high oxidative stress allelic risk score. The pathway-specific genetic variants used to generate the allelic risk scores are presented in Table 2. Participants with a low endothelial function or a low oxidative stress allelic risk score did not differ in mean 1-year BC exposures from their counterparts with high allelic risk scores (Figure 1A). Participants with a low metal-processing allelic risk score had a higher mean 1-year BC exposure (mean [SE], 0.53 [0.01] μg/m3; P = .02) when compared with individuals with a high metal-processing allelic risk score (mean [SE], 0.50 [0.007] μg/m3). No differences in mean IOP were found between individuals with low and those with high allelic risk scores across the 3 pathways (Figure 1B).
Table 3 presents the results of the tiered framework we used to examine the direct associations of 1-year BC exposure with mean IOP. We did not observe a direct association of IQR increases in 1-year BC exposure with mean IOP in the tier 1 model, which was adjusted for age and year of visit (β = 0.07; 95% CI, −0.23 to 0.37; P = .64). In addition, we did not observe an association after adjusting for important covariates and confounders in the tier 2 model (β = 0.13; 95% CI, −0.18 to 0.45; P = .40), the final fully adjusted tier 3 model (β = 0.14; 95% CI, −0.18 to 0.45; P = .40), or the tier 4 model (β = 0.13; 95% CI, −0.19 to 0.44; P = .42). Again, tier 3 was selected as the final fully adjusted framework because it reflected the best model fit (ie, lowest Akaike information criterion, including important lifestyle and disease covariates). No association with BC was observed when we modeled BC exposure and the 3 allelic risk scores as joint predictors of mean IOP. Moreover, none of the 3 allelic risk scores had direct associations with mean IOP.
Figure 2 depicts the modifying association of the 1-year BC exposure with mean IOP. Of all of the allelic risk scores for each of the 3 pathways (endothelial function, metal processing, and oxidative stress), only the oxidative stress allelic score changed the association of BC with IOP (false discovery rate–corrected P = .03 for interaction). Specifically, the association was greater in individuals who had high oxidative stress allelic risk scores (β = 0.36; 95% CI, 0.003-0.73) when compared with individuals who had low oxidative stress allelic risk scores (β = −0.35; 95% CI, −0.86 to 0.15). Our results remained unchanged in the sensitivity analysis using the tier 4 model (ie, glaucoma added as a variable in the fully adjusted tier 3 model) (Figure 2). In the sensitivity analysis, in which the variables for glaucoma eyedrops and antihypertensive medications were used instead of glaucoma and hypertension disease status, again high oxidative stress allelic risk scores (β = 0.36; 95% CI, –0.008 to 0.72) were compared with low oxidative stress allelic risk scores (β = −0.34; 95% CI, –0.84 to 0.17).
In the present study, we describe the novel association among long-term BC exposure, oxidative stress polymorphisms, and IOP in a cohort of community-dwelling men. These associations persisted even after adjusting for important disease states, such as inflammation, glaucoma, obesity, ischemic heart disease, and diabetes. These results suggest that ambient BC exposure may be a risk factor for elevated IOP in individuals at risk for elevated oxidative stress that is independent of discernable disease diagnoses. To our knowledge, this is the first study to report the association of an ambient air pollutant with intraocular pressure on a population level. As such, this study may point to the potential need to broaden the factors considered when evaluating and managing elevated IOP.
It has long been postulated that environmental air pollution may affect eye health; however, evidence supporting this theory is limited. The present study addresses this research gap. Although we considered genetic polymorphisms from 3 different pathways (endothelial function, metal processing, and oxidative stress), we observed findings only in oxidative stress genetic risk. In addition, the polymorphisms used to generate these genetic scores were selected using the independent outcomes of pathways related to the polymorphisms. In the case of the oxidative stress allelic risk score, this score was developed from polymorphisms associated with serum levels of 8-hydroxy-2′-deoxyguanosine,5 which is one of the most common and quantifiable DNA adducts that reflects cellular oxidative stress.32 Thus, the oxidative stress score, although still a measure of genetic risk, is associated with a tangible, physiologic measure of oxidative stress.
When individuals with high or low oxidative stress allelic risk scores were compared, we detected a moderate difference in mean IOP (0.73 mm Hg) for an IQR increase in 1 year of BC exposure. These results suggest that BC exposure susceptibility may be biologically important. Independent of ambient air pollution, increased oxidative stress or reduced antioxidant capacity has been linked to elevated IOP and glaucoma pathologic findings.33,34 For example, a 2017 study compared blood and aqueous humor samples in 96 patients with glaucoma with samples from 64 age-matched healthy controls. The researchers observed higher levels of protein carbonyls and other oxidative stress–related molecular measures in the blood and more prominently in the aqueous humor from individuals with glaucoma.35 Another study of 531 individuals reported an association between systemic oxidative stress and IOP. The researchers first quantified biological antioxidant potential in the participants’ serum and then observed that individuals with the lowest biological antioxidant potential were more likely to have higher IOP.36 The conclusions reached by these observational studies are further supported by experiments in animal models, which provide a greater mechanistic understanding of the association of oxidative stress with IOP physiologic features.37,38 Additional support for these conclusions can be derived from human studies that suggest that oxidative stress–reducing activities (eg, diet, exercise) may curtail some IOP-related conditions.39,40 Exposure to BC has been identified as a source of localized and systemic oxidative stress8; thus, BC may be acting through some of the same oxidative stress physiologic functions already associated with glaucoma and IOP. Still, further research is necessary to confirm our findings and to better establish this association.
We found a novel gene-environment interaction associated with IOP, which is consistent with findings from other studies investigating gene-environment interactions and smoking. A 2011 study reported interactions between the NOS3 gene rs7830 single-nucleotide polymorphism and smoking, which can be viewed as a personal form of air pollution. Specifically, CC homozygous individuals who were past or current smokers were at a greater risk of developing primary open-angle glaucoma when compared with CC homozygous individuals who were never smokers. However, no association with smoking was found in carriers of the single-nucleotide polymorphism A variant.41 Like the present study, the study by Kang et al41 revealed an association that may have been missed if gene-environment interactions had not been considered. Moreover, it also exposed a more nuanced but important point: Not all environmental exposures may be directly associated with IOP-related eye diseases. Future research into why some exposures are directly, indirectly, or not associated with eye diseases will be useful for better understanding disease pathogenesis and potentially safeguarding vulnerable populations from disease onset.42
The strengths of the present study include the use of novel genetic pathway tools and a longitudinal cohort as well as the repeated measures of ambient pollutant exposures, IOP, covariates, and potential confounders. However, the study has some limitations. First, we used an allelic risk score approach that does not provide genome-wide resolution of the 3 biological pathways and that allots polymorphisms equal weights when they may have different functional magnitudes.5 Still, the variants we used were developed in an independent study and have been used in other studies since they were first described.43,44 Moreover, the variants are representative of their respective biological pathways. Second, we used a validated spatiotemporal model to estimate 1-year BC levels at the participant home. Given that most participants are retired and spend most of their time at home, these estimates were believed to be good proxies of the participants’ personal exposure. Moreover, any nondifferential misclassification bias from this approach of estimating personal exposures is likely to attenuate statistical associations rather than steer them away from the null.45,46 Last, this study cohort consisted of older white men residing in a lightly polluted environment, the lower bound of the 95% CI for the observed association in the high oxidative stress genetic risk group was close to zero, and the design did not reflect causality. These final points reflect the need for additional studies involving other environments, other demographics, and larger groups to confirm our findings more broadly.
This study highlights the potential contributions of gene-environment interactions to the complex physiologic functions of IOP-related disease. This finding could be particularly important because most of the global burden of eye disease is in the developing world.47 The differences in the global prevalence of eye disease are often attributed to genetics alone or the difficulties with health care access that affect early disease detection, management, and treatment.48 These factors are highly relevant but also note that the developing world often faces a larger burden of environmental pollution.49 Annual BC exposure levels greater than 5 times the levels observed in our study have been reported in India and China.50,51 Of importance, the present associations, although novel, do not prove causality and should not affect policy in isolation. Whether these findings persist in more diverse populations experiencing greater pollution and in study designs that can demonstrate causality will be interesting to see. If these future studies substantiate this association, integrated initiatives (ie, combining environmental improvement, socioeconomic outreach, and targeted pharmaceutical interventions) may prove useful for future policy or public health initiatives aimed at addressing the global burden of eye disease.
Accepted for Publication: September 5, 2018.
Corresponding Author: Jamaji C. Nwanaji-Enwerem, PhD, Department of Environmental Health, Harvard T. H. Chan School of Public Health and MD-PhD Program, Harvard Medical School, 665 Huntington Ave, Building 1, Room G7, Boston, MA 02115 (firstname.lastname@example.org).
Published Online: November 8, 2018. doi:10.1001/jamaophthalmol.2018.5313
Author Contributions: Dr J. C. Nwanaji-Enwerem had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Concept and design: O. Nwanaji-Enwerem, Baccarelli, Park.
Acquisition, analysis, or interpretation of data: J. C. Nwanaji-Enwerem, Wang, Vokonas, Baccarelli, Weisskopf, Herndon, Wiggs, Park, Schwartz.
Drafting of the manuscript: J. C. Nwanaji-Enwerem, Vokonas.
Critical revision of the manuscript for important intellectual content: J. C. Nwanaji-Enwerem, Wang, O. Nwanaji-Enwerem, Baccarelli, Weisskopf, Herndon, Wiggs, Park, Schwartz.
Statistical analysis: J. C. Nwanaji-Enwerem, Wang, Schwartz.
Obtained funding: J. C. Nwanaji-Enwerem, Vokonas, Baccarelli, Schwartz.
Administrative, technical, or material support: Wang, O. Nwanaji-Enwerem, Vokonas, Baccarelli, Park.
Supervision: Vokonas, Baccarelli, Weisskopf.
Conflict of Interest Disclosures: Dr Weisskopf reported receiving grants from the National Institutes of Health (NIH) during the conduct of the study; consulting fees from Spaulding Rehabilitation Hospital; and compensation from Springer for editor duties. Dr Wiggs reported receiving grants from the National Eye Institute during the conduct of the study. Dr Schwartz reported receiving grants from the NIH during the conduct of the study. No other disclosures were reported.
Funding/Support: This work was funded in part by grants R01ES021733, R01ES025225, R01ES027747, ES015172, ES014663, ES020010, P30ES009089, P30ES000002, K01ES016587, and P30ES017885 from the National Institute of Environmental Health Sciences; grants RD832416 and RD83587201 from the Environmental Protection Agency; grant 2T32HL007118-41 from the National Heart, Lung, and Blood Institute; grant T42-OH008455 from the National Institute for Occupational Safety and Health; Ruth L. Kirschstein National Research Service Award 1 F31AG056124-01A1 from the NIH/National Institute on Aging (Dr J. C. Nwanaji-Enwerem); and contract 53-K06-510 from the US Department of Agriculture, Agricultural Research Service. The US Department of Veterans Affairs Normative Aging Study, supported by the Cooperative Studies Program/Epidemiologic Research and Information Center and the US Department of Veterans Affairs, is a research component of the Massachusetts Veterans Epidemiology Research and Information Center.
Role of the Funder/Sponsor: The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Disclaimer: The views expressed herein are those of the authors and do not reflect the official policy or position of the US Department of Veterans Affairs or any of the funding agencies.
Additional Contributions: Emilia J. Ling, SM, and Brett P. Salazar, BSc, both from Stanford University School of Medicine, assisted with the promotional image. These individuals were not compensated for their contributions.
Create a personal account or sign in to: