Visual Overview
Abstract
Elevated levels of fine particulate matter <2.5 µm in aerodynamic diameter (PM2.5) are associated with increased risk of cardiovascular outcomes and death, but their association with risk of CKD and ESRD is unknown. We linked the Environmental Protection Agency and the Department of Veterans Affairs databases to build an observational cohort of 2,482,737 United States veterans, and used survival models to evaluate the association of PM2.5 concentrations and risk of incident eGFR <60 ml/min per 1.73 m2, incident CKD, eGFR decline ≥30%, and ESRD over a median follow-up of 8.52 years. County-level exposure was defined at baseline as the annual average PM2.5 concentrations in 2004, and separately as time-varying where it was updated annually and as cohort participants moved. In analyses of baseline exposure (median, 11.8 [interquartile range, 10.1–13.7] µg/m3), a 10-µg/m3 increase in PM2.5 concentration was associated with increased risk of eGFR<60 ml/min per 1.73 m2 (hazard ratio [HR], 1.21; 95% confidence interval [95% CI], 1.14 to 1.29), CKD (HR, 1.27; 95% CI, 1.17 to 1.38), eGFR decline ≥30% (HR, 1.28; 95% CI, 1.18 to 1.39), and ESRD (HR, 1.26; 95% CI, 1.17 to 1.35). In time-varying analyses, a 10-µg/m3 increase in PM2.5 concentration was associated with similarly increased risk of eGFR<60 ml/min per 1.73 m2, CKD, eGFR decline ≥30%, and ESRD. Spline analyses showed a linear relationship between PM2.5 concentrations and risk of kidney outcomes. Exposure estimates derived from National Aeronautics and Space Administration satellite data yielded consistent results. Our findings demonstrate a significant association between exposure to PM2.5 and risk of incident CKD, eGFR decline, and ESRD.
Elevated levels of fine particulate matter of <2.5 µm in aerodynamic diameter (PM2.5) are associated with increased risk of death, reduced life expectancy, increased risk of cardiovascular disease, and stroke, as well as a host of other adverse health outcomes.1–13 Studies on the relationship of PM2.5 and the kidney are limited.
Experimental laboratory evidence suggests that exposure to deep exhaust particles leads to disturbances in renal hemodynamics, promotes oxidative stress, inflammation, and DNA damage in renal tissue, exacerbates AKI, and further promulgates chronic renal injury in murine models14,15
Data on the relationship of air pollution and kidney disease in humans is very scarce. Hendryx described higher kidney disease mortality in coal mining Appalachian areas in the United States and suggested that this increase in kidney disease mortality burden may reflect environmental exposure to particulate matter.16 In a cross-sectional analysis of 1103 Boston-area patients hospitalized for ischemic stroke, Lue et al.17 showed that residential proximity to major roads is associated with reduced eGFR. In seminal findings, Mehta et al.18 examined the relationship between long-term exposure to PM2.5 and longitudinal changes in eGFR in a regional cohort of 669 older men who were participants in the Boston-based Veterans Administration Normative Aging Study, and found that higher 1-year PM2.5 exposure was associated with lower eGFR and an additional annual decrease in eGFR. In a large Chinese study, Xu et al.19 examined an 11-year (for indication) renal biopsy series of 71,151 native kidney biopsies and reported that long-term exposure to high levels of PM2.5 was associated with increased odds of membranous nephropathy, but not other glomerular diseases.
The sum of the experimental and clinical findings provides biologic plausibility and supports the hypothesis that environmental exposure to elevated levels of PM2.5 is associated with increased risk of kidney disease. However, whether exposure to elevated levels of PM2.5 is associated with increased risk of development of CKD and CKD progression has not been examined in large national longitudinal epidemiologic studies. Identification of air pollution as a potential contributor to kidney disease will inform national and global burden of disease estimates, stimulate further policy discussions on the importance of curbing air pollution on health and disease, and inform the public about the hazards of air pollution. In this work, we built a longitudinal national cohort of 2,482,737 United States veterans and characterized the relationship of PM2.5 and risk of incident CKD, and progression to ESRD.
Results
There were 2,482,737 cohort participants followed for a median of 8.52 years (interquartile range, 8.04–8.80). PM2.5 concentrations ranged from 5.0 to 9.1 µg/m3, 9.2 to 11.0 µg/m3, 11.1 to 12.6, and 12.7 to 22.1 µg/m3 in quartiles 1, 2, 3, and 4, respectively. Overall cohort participants were mostly white men (Table 1). Cohort participants living in counties in the highest quartile of PM2.5 concentrations were more likely to be black, more likely to have hypertension, diabetes mellitus, cardiovascular disease, and higher eGFR at time of cohort entry (Table 1). They were also more likely to be current or former smokers (Table 1). Counties in the highest quartile of PM2.5 concentrations had significantly higher population density. There was a gradual increase in incident rate of eGFR <60 ml/min per 1.73 m2, CKD, eGFR decline ≥30%, and ESRD across quartiles of county level PM2.5 concentrations (Table 2). Adjusted survival curves by PM2.5 quartiles are presented in Figure 1, A–D.
Demographic and health characteristics of overall study cohort and according to quartiles of annual average county PM2.5 concentrations
Incident rate of renal outcomes in overall study cohort and according to quartiles of annual average county PM2.5 concentrations
Adjusted survival curves by PM2.5 quartiles. (A) Incident eGFR <60 ml/min per 1.73 m2, (B) incident CKD, (C) eGFR decline ≥30%, and (D) ESRD. Survival curves are adjusted for age, race, sex, and T0 eGFR.
Exposure to Ambient Fine Particulate Matter and the Risk of Incident eGFR <60 ml/min per 1.73 m2 and Incident CKD
We examined the association of PM2.5 concentrations and the risk of incident eGFR <60 ml/min per 1.73 m2 among those who had no history of eGFR<60 ml/min per 1.73 m2 before time of cohort entry. Where exposure was defined at baseline as the annual average PM2.5 concentrations in 2004, a 10-µg/m3 increase in PM2.5 concentration was associated with increased risk of incident eGFR <60 ml/min per 1.73 m2 (hazard ratio [HR], 1.21; 95% confidence interval [95% CI], 1.14 to 1.29) (Table 3). In analyses where exposure was time-varying throughout duration in the cohort (exposure for each cohort participant was matched with his/her county of residence at a given time and the county’s average annual particulate matter concentration for that year), the risk was also increased (HR, 1.25; 95% CI, 1.17 to 1.34) (Table 3). In analyses where exposure was considered in quartiles, compared with cohort participants in the lowest quartile (quartile 1) of PM2.5 concentrations, those in quartile 2, 3, and 4 had graded increased risk of incident eGFR <60 ml/min per 1.73 m2 in analyses considering baseline 2004 exposure and time-varying exposure (Supplemental Table 1). Spline analyses suggested a linear relationship between PM2.5 concentrations and risk incident eGFR <60 ml/min per 1.73 m2 (P for nonlinearity =0.90); the linear representation of the association of PM2.5 concentrations and risk is depicted in Figure 2A.
Risk of renal outcomes for every 10 μg/m3 increase in PM2.5 concentrations:
Analyses of risk of renal outcomes by PM2.5 concentrations (PM2.5 of 5.7 μg/m3 served as a reference) with PM2.5 probability distribution in the background. (A) Risk of incident eGFR <60 ml/min/1.73 m2 (P for nonlinearity =0.90). (B) Risk of incident CKD (P for nonlinearity =0.90). (C) Risk of eGFR decline ≥30% (P for nonlinearity =0.84). (D) Risk ESRD (P for nonlinearity =0.47). Models adjusted for age, race, sex, cancer, cardiovascular disease, chronic lung disease, diabetes mellitus, hyperlipidemia, hypertension, T0 eGFR, BMI, smoking status, angiotensin-converting enzyme inhibitor/angiotensin receptor blocker use, county population density, number of outpatient eGFR measurements, number of hospitalizations, and county percent in poverty.
We evaluated the risk of incident CKD (defined as two eGFR measurements <60 ml/min per 1.73 m2 at least 90 days apart) in a subcohort of people with at least two eGFR measurements separated by at least 90 days during follow-up and who had no history of eGFR <60 ml/min per 1.73 m2 before time of cohort entry. The results were consistent (Figure 2B, Supplemental Table 1, Table 3).
Exposure to Ambient Fine Particulate Matter and the Risk of eGFR Decline and Progression to ESRD
In the overall cohort, an increase in PM2.5 of 10 µg/m3 was associated with increased risk of eGFR decline ≥30% in analyses considering baseline exposure (HR, 1.28; 95% CI, 1.18 to 1.39), and time-varying exposure (HR, 1.36; 95% CI, 1.26 to 1.46) (Table 3). Compared with the lowest quartile (quartile 1) of PM2.5 concentrations, quartiles 2, 3, and 4 exhibited a graded increase in risk of eGFR decline ≥30% (Supplemental Table 1). There was a linear relationship between PM2.5 concentrations and risk of eGFR decline ≥30% (Figure 2C) (P for nonlinearity =0.84). The results were consistent in analyses considering the outcome of ESRD (Figure 2D, Supplemental Table 1, Table 3).
PM2.5 and Risk of Adverse Kidney Outcomes Using National Aeronautics and Space Administration Data
We additionally considered PM2.5 estimates derived from the National Aeronautics and Space Administration (NASA) space-borne satellite sensors as an alternative data source to define ambient PM2.5 exposure levels. Analyses considering the NASA satellites remote sensing data yielded results consistent with those shown using exposure data obtained from ground-based air monitoring stations (Table 4).
Risk of renal outcomes for every 10 μg/m3 increase in PM2.5 concentration using NASA Socioeconomic Data and Applications Center data for exposure levels
Population Attributable Fraction
Population attributable fraction (PAF) represents the proportional reduction in population disease that would occur if exposure to PM2.5 was reduced to the Environmental Protection Agency’s (EPA) recommended levels of 12 µg/m3. PAF for incident eGFR <60 ml/min per 1.73 m2 and incident CKD were 1.51% (95% CI, 1.44 to 1.58), and 2.06% (95% CI, 1.96 to 2.17), respectively. PAF for eGFR decline ≥30%, and ESRD were 2.21% (95% CI, 2.14 to 2.30), and 2.31% (95% CI, 1.86 to 2.75), respectively. Our estimate of the national burden of CKD attributable to elevated levels of PM2.5 exceeding the EPA standard (where the theoretical minimum risk exposure level [TMREL] was set at the EPA standard of 12 µg/m3) in the contiguous United States was 44,793 incident cases per year (95% uncertainty interval [95% UI], 42,716 to 46,869). The national burden of ESRD attributable to PM2.5 levels in excess of EPA standards was 2438 incident cases per year (95% UI, 1963 to 2902). We conducted analyses where we defined TMREL on the basis of the methodologies of the Global Burden of Disease (GBD) studies; this TMREL was assigned on the basis of a uniform distribution of PM2.5 from 2.4 to 5.9 μg/m3 representing exposure values between the minimum and fifth percentiles of exposure distributions from outdoor air pollution cohort studies.20–22 Using the GBD definition of TMREL, the estimate of the national burden of incident CKD attributable to air pollution was 337,032 (95% UI, 207,976 to 466,087) and the national burden of incident ESRD attributable to air pollution was 13,537 (95% UI, 8377 to 18,696).
Maps depicting the geographic distribution of the burden of incident CKD where the theoretical minimum was set at the EPA standard, and according to the GBD methodologies are presented in Figures 3, A and B, respectively.
Geographic distribution of the national burden of incident CKD attributable to air pollution in the United States. (A) Risk attributable to exposure levels of PM2.5 above the EPA recommended concentration of 12 μg/m3. (B) Risk attributable to exposure levels of PM2.5 above a uniform distribution between 2.4 and 5.9 µg/m3. Burden is per 100,000 population.
Sensitivity Analyses
In order to test different distance thresholds (and spatial resolutions) for exposure definition, we assigned PM2.5 exposure levels to each cohort participant on the basis of the nearest air monitoring station in those who were within 30, 10, and 5 miles of an air monitoring station (Supplemental Table 2); the results were consistent and yielded slightly higher risk estimates in time-varying analyses where the maximum distance from an air monitoring station was <5 miles (Supplemental Table 2). In order to evaluate the robustness of study findings in the most populated counties in the United States, the analyses were restricted to the 100 counties with highest numbers of population, and results were consistent in terms of direction and magnitude of risk estimates (Supplemental Table 3). Because variation in regional characteristics may confound the association of PM2.5 and kidney disease, we developed strategies to evaluate the association in analyses, considering city-adjusted and within-city models (as described in Supplemental Material); the results suggest that within metropolitan areas (city-adjusted or within-city models), higher PM2.5 concentrations were associated with higher risk of renal outcomes (Supplemental Table 4). Between-city risk estimates were slightly higher than within-city estimates (Supplemental Table 4). Because characteristics of geographies might confound the association between PM2.5 and risk of kidney outcomes,23 we curated the County Health Rankings datasets24,25 and built analyses additionally controlling for 55 United States county-level variables in six domains, including demographics, physical environment, social and economic factors, health behaviors, clinical care, and health outcomes; the results remained consistent (Supplemental Table 5). We repeated the analyses using expanded and more sensitive definitions of hypertension and diabetes: the definition of hypertension included relevant diagnostic codes and average systolic BP (treated as a continuous variable) in the year before time zero (T0); the definition of diabetes included relevant diagnostic codes as well as use of diabetic medications (including oral hypoglycemic agents and insulin) and hemoglobin A1c levels >6.4%. The results remained stable to this challenge of using expanded definitions of key drivers of kidney disease (Supplemental Table 6). The results were reproduced in analyses considering alternative renal outcomes, including odds of rapid eGFR decline (eGFR slope <−5 ml/min per 1.73 m2 per year), and risk of ESRD or eGFR decline ≥50% (Supplemental Table 7). As a measure of calibration, we examined the association of PM2.5 and risk of death where a priori observations suggest that an association is expected (positive control).6,10 Our results show a significant association between PM2.5 concentrations and risk of death (Supplemental Table 8). We also considered the outcome of myocardial infarction as an additional positive control3,26; the results were consistent with published literature6,10 in that an increase in PM2.5 was associated with increased risk of myocardial infarction (Supplemental Table 9). Results of sensitivity analyses for the competing risk of death were consistent with those shown in primary analyses27 (Supplemental Table 10).
Negative control is a valuable complement to other epidemiologic methods and serves to identify and resolve both suspected and unsuspected sources of spurious causal inference, including confounding, mismeasurements, and other biases, design, or analytic flaws.28 Ambient air sodium concentration is one of the parameters measured by air monitoring stations. There is no biologic basis to support an association between levels of sodium concentrations in the air and risk of adverse renal outcomes; this renders ambient air sodium a suitable negative control.28 We therefore, as a negative exposure control, tested the association between ambient air sodium levels and the risk of renal outcomes, and the results show a vanishingly weak or nonsignificant association in models, considering baseline exposure in the year 2004 and time-varying exposure (Supplemental Table 11). An analysis considering the association between air sodium levels and risk of death also shows no association in models considering baseline exposure (HR, 1.00; 95% CI, 1.00 to 1.01) and those considering time-varying exposure (HR, 1.00; 95% CI, 0.99 to 1.00) (Supplemental Table 11). Additional details on sensitivity analyses are presented in Supplemental Material.
In formal interaction analyses, race did not modify the association of PM2.5 and risk of kidney outcomes; the association was more pronounced among women and in those below the median age and in those below the median body mass index (BMI) of the overall cohort (Supplemental Figure 1).
Discussion
In a large national cohort of United States veterans, we observed a linear relationship between PM2.5 concentrations and risk of incident CKD and progression to ESRD. The results were consistent where baseline exposure was defined as the annual average PM2.5 concentrations in the year 2004, and where exposure was time-varying to reflect movement of cohort participants and changes in PM2.5 concentrations over the years. Furthermore, we examined a range of kidney outcomes including development of kidney disease, kidney function decline (eGFR decline ≥30%), and the terminal outcome of ESRD. The results consistently showed a linear relationship between PM2.5 levels and risk of kidney outcomes. The results were robust in sensitivity analyses including the examination of different distance thresholds from an air monitoring station, and analyses evaluating the association within metropolitan areas. The results were also consistent in analyses using ambient PM2.5 estimates derived from NASA’s satellite data. Ambient air sodium concentrations (used as a negative control) were not associated with increased risk of adverse renal outcomes. The constellation of findings suggests that chronic exposure to fine particulate matter air pollution is a significant risk factor for the development and progression of kidney disease.
Although air quality has significantly improved in the United States over the past several decades, air quality remains suboptimal in many parts of the country and in multiple geographies around the world. In our analyses, the risk of CKD and its progression was most pronounced at the highest levels of fine particulate matter concentrations. However, analyses where PM2.5 concentrations were categorized in quartiles suggest a graded relationship; spline analyses (and tests of nonlinearity) endorse a linear relationship where risk starts to increase at relatively low concentrations of particulate matter that are well below the recommended levels by the World Health Organization (10 µg/m3) and the EPA (12 µg/m3). We used a conservative approach to estimate the attributable burden of CKD and ESRD in the contiguous United States, and our results suggest a small but significant number of incident cases of CKD (44,793 per year) and ESRD (2438 per year) are attributable to levels of particulate matter air pollution exceeding the EPA recommended level of 12 µg/m3. A lower and more stringent TMREL of PM2.5 between 2.4 and 5.9 μg/m3 yielded a higher estimate of incident CKD (337,032 per year) and incident ESRD (13,537 per year). The findings provide a quantitative assessment of the potential reduction in burden of CKD that is achievable with improvement in air quality in the United States, and suggest the need for a broader assessment of the global burden of kidney disease attributable to air pollution.
Multiple studies reported substantial geographic variation in the burden of CKD and ESRD in the United States and globally; the variation was persistent after accounting for diabetes, hypertension, and obesity conditions generally considered to be major drivers of kidney disease, suggesting that variation in burden of disease is likely due to factors other than these traditional drivers.23,29 It has been hypothesized that some of the geographic variation in kidney disease burden may be due to environmental factors.23 Our study results suggest that particulate matter air pollution is an important, but yet unrecognized risk factor for kidney disease and its progression, which may explain some of the geographic variation in kidney disease burden in the United States. Examination of the effect of particulate matter on risk of incident kidney disease and its progression outside the United States, and particularly in areas with much higher levels of particulate matter air pollution is warranted.
Three distinct hypotheses have been proposed to explain the epidemiologic observations of a relationship between PM2.5 and cardiovascular outcomes; these may also be pertinent in the evaluation of renal outcomes.30 The first suggests that inhaled particles provoke pulmonary inflammation which may then lead to systemic inflammation.30 The second suggests that the mechanism involves pollutant-induced disturbances in the lung autonomic nervous system.30 The third (and most plausible hypothesis) is predicated on the premise that airborne particulates enter the bloodstream where they may then interact with tissue components to promote the observed pathologic effects30,31; the latter is supported by emerging evidence suggesting that inhaled inert gold nanoparticles not only enter the bloodstream of healthy adult volunteers, but are detected in the urine within minutes after exposure, providing a proof of concept that inhaled nanoparticles get filtered and excreted by the kidney.31 These three hypotheses provide contextual background to evaluate the experimental and clinical findings describing the extrapulmonary effect of particulate matter air pollution, where it has been reported that exposure to elevated levels of PM2.5 is associated with increased inflammatory mediators (including TNF-α, IL-6, and plasminogen activator inhibitor-1), oxidative stress,32–34 increased atherosclerotic plaque area, and exaggerated vasoconstrictor responses to phenylephrine and serotonin.35 Evidence suggests that increased PM2.5 concentrations are associated with significant decrease in flow-mediated dilatation,36,37 increases in systolic BP and pulse pressure,38–40 and disturbances in the hypothalamic-pituitary-adrenal axis.41 Emerging evidence also suggests that exposure to ambient air pollutants can lead to metabolic disturbances, including glucose intolerance, decreased insulin sensitivity, higher blood lipid concentrations, weight gain, and increased risk of diabetes mellitus.42–44 It is plausible that one or more of these mechanistic pathways may explain the association described here.
This study has several limitations. Cohort participants were United States veterans and mostly older white men; therefore, the findings may not be generalizable to other populations. Although we accounted for known confounders, the possibility of residual confounding due to either unknown or unmeasured confounders cannot be completely excluded. We specifically note that the counties with the highest measures of PM2.5 had higher population density and greater percentage of citizens living in poverty. Although we took care to develop analyses for the 100 most populated counties, built city-adjusted and within-city analyses, and also accounted for United States county-level characteristics23 in domains including demographics, physical environment, social and economic conditions, health behaviors, clinical care, and health outcomes, it is possible that individual differences in genomic makeup, dietary habits, physical activity, other environmental or occupational attributes, variation in exposure to heavy metals, or other factors not captured in our analyses might explain the described results. To define covariates we relied on VA administrative databases, and although we used comorbidity definitions validated for use in administrative datasets, misclassification or inaccurate measurement of predictor variable is not impossible. Our analyses did not consider the composition of PM2.5, which may vary by region and over time; however, estimates using nonspecific PM2.5 mass alone may underestimate the total effect of PM2.5 on health outcomes and bias the results toward the null hypothesis.10 Our datasets did not contain information on time spent in traffic or outdoors, and did not include information on potential exposure to indoor air pollutants, which may result in misclassification of exposure. The study has a number of strengths, including the large national cohort of veterans who are recipients of care in a single integrated network of health care systems designed to reduce variation in care practices, and the analyses benefited from the merging of large databases including those of the VA, EPA, Census Bureau, and other data sources. We followed our cohort participants for a median of 8.52 years (interquartile range, 8.04–8.80), designed analytic strategies using time-varying exposure (to capture updated exposure levels as participants moved over the years), and evaluated a range of well defined chronic kidney outcomes, including development of kidney disease, CKD progression, and the terminal outcome of ESRD. Our analytic approach to examine within-city and city-adjusted effect reduces concern about confounding due to variation in regional characteristics, and potential geographic variation in composition and toxic content of PM2.5. The analyses also considered NASA’s satellite data as alternative data source to define PM2.5 exposure. The analytic strategies also included the development and testing of positive and negative controls to detect possible hidden bias. In summary, our results demonstrate a significant association between PM2.5 concentrations and risk of development of kidney disease, and its progression to ESRD. Effort to improve air quality might ease the burden of kidney disease in the United States and globally.
Concise Methods
Cohort Participants
We selected users of the Veterans Affairs (VA) Healthcare System, using data from the US Department of Veterans Affairs, who had at least one outpatient eGFR measurement between October 1, 2003 and September 30, 2004 and no prior history of ESRD, and designated the date of last eGFR measurement in this time period as T0 (n=2,751,717). Patients were further chosen on having at least one eGFR measurement after T0 (n=2,680,431), and were followed until September 30, 2012 or death. Participants were then limited to those who had data on PM2.5 (n=2,628,465) and data on all covariates, yielding an analytic cohort of 2,482,737 (Figure 4). The study was approved by the Institutional Review Board of the VA Saint Louis Health Care System (Saint Louis, MO).
Flow diagram of cohort participant inclusion.
Data Sources
Department of Veterans Affairs datasets were used to obtain patient demographics, inpatient and outpatient data, laboratory information, vital signs, and medications.45,46 Details on the VA datasets used are provided in Supplemental Material.47–51 Data from the US Renal Database System (USRDS) was used to supplement ESRD status information. The Center for Disease Control’s (CDC) National Environmental Public Health Tracking Network furnished annual particulate matter estimates for the contiguous United States that originate from Community Multiscale Air Quality modeled output.52,53 The Community Multiscale Air Quality System uses, for counties with air monitoring stations, the measures from the EPA’s Air Quality System. EPA data also provided information on sodium levels, as well as the latitude and longitude of said monitoring stations.54 The NASA Socioeconomic Data and Applications Center Global Annual PM2.5 grids from Moderate Resolution Imaging Spectroradiometer, Multi-angle Imaging Spectroradiometer, and Sea-Viewing Wide Field-of View Sensor aerosol optical depth remote space-borne satellite sensing data provided an additional source of ambient PM2.5 estimates at the 10×10 km resolution.55,56 National United States–based estimates of incident rates of CKD and treated ESRD were obtained from the CDC CKD Surveillance Project and the 2016 USRDS Annual Data Report, respectively.57,58 Data on county-level poverty, population, population density, metropolitan statistical areas, and zip code centroid were obtained from the US Census Bureau. More detailed description of data sources is provided in Supplemental Material.
Exposure Assessment
The primary predictor variable for analyses was PM2.5 concentrations. Cohort participants were assigned geographic location, which may have varied over time, on the basis of their county information contained in outpatient or inpatient data closest to but before a given time point. Using annual monitor and modeled data, exposure (in micrograms per cubic meter) was defined as (1) the annual average in year 2004, where a patient’s geographic location was designated as location at T0 (used in baseline models); and (2) time varying where geographic location was updated as cohort participants moved (and average annual exposure was matched to their updated geographic location at any specific time). In all primary analyses, unless otherwise indicated, measures correspond to a 10-μg/m3 increase in PM2.5.9 PM2.5 exposure was additionally categorized into quartiles defined by the county-level distribution. In the time-varying model, quartiles were independently defined by exposure distributions among the counties in each given year. PM2.5 was alternatively defined by NASA’s SEDAC Global Annual PM2.5 grids from MODIS, MISR and SeaWiFS aerosol optical depth data through linkage to participant’s zip code of residence.55,56 Further details are provided in Supplemental Material.
Ascertainment of Outcomes
Outcomes evaluated included the risk of incident eGFR <60 ml/min per 1.73 m2, the risk of incident CKD where CKD was defined as two eGFR measurements <60 ml/min per 1.73 m2 at least 90 days apart,47 time until ≥30% decline in eGFR from eGFR at T0, and time until ESRD.59 Patients were censored after onset of ESRD, for all outcomes other than ESRD, and at time of death or end of study follow-up. The date of first ESRD services was ascertained through linkage of VA and USRDS databases. Outpatient eGFR was used in the assessment of all outcomes except for ESRD. eGFR was calculated using the four-variable abbreviated CKD Epidemiology Collaboration equation on the basis of age, race, sex, and serum creatinine.60
Covariates
Covariate selection was on the basis of factors that could potentially confound the association of PM2.5 and kidney disease outcomes, and was informed by prior studies.3,23,47,61–64 Baseline covariates were ascertained from October 1, 1999 until cohort entry (T0). Covariates included age, race, sex, cancer, cardiovascular disease, chronic lung disease, diabetes mellitus, hyperlipidemia, hypertension, T0 eGFR, BMI, smoking status, angiotensin-converting enzyme inhibitor/angiotensin receptor blocker use, county population density, number of outpatient eGFR measurements, number of hospitalizations, and county percent in poverty. Details of covariate definitions are provided in Supplemental Material. Covariates were treated as continuous variables where appropriate, unless otherwise indicated.
Statistical Analyses
Demographic and clinical characteristics of the overall cohort and by county PM2.5 quartile are presented as frequency (percent) for categorical variables, and as mean (standard deviation) or median (interquartile range) for continuous variables if normally or non-normally distributed, respectively. PM2.5 quartiles for Table 1 and baseline models were derived from the distribution of 2004 monitor and modeled PM2.5 among counties. Age, race, sex, and eGFR adjusted incident rates are presented, and were standardized to the distribution of these variables in the overall cohort. Survival curves for PM2.5 quartiles, adjusted for age, race, sex, and eGFR, are presented. Cox proportional hazard models were adjusted for covariates as described. In order to account for intracounty correlation, a robust sandwich estimator was used. Analyses were performed in baseline models and time varying models, where PM2.5 exposure was treated as time varying. Patients were excluded from time-varying analyses if their county of residence had missing data at any time during follow-up. The inclusion of time-varying models in the analytic strategies was in consideration of the fact that cohort participants may have moved over time, but also importantly in recognition of the fact that the Clean Air Act and other environmental policies at state and local levels have resulted in reduction in PM2.5 over the years.65 The primary analyses were repeated utilizing PM2.5 exposures derived from NASA satellite data. Effect modification of the association between PM2.5 exposure and renal outcomes by age, race, sex, and BMI was examined through the addition of interaction terms in separate baseline models using above and below median age, black and nonblack, men and women, and above and below median BMI interactions.
Cubic spline analyses were performed.66 A Wald chi-squared test for nonlinearity of spline terms indicated that, for all outcomes, there was no statistical evidence of deviation from linearity. Further information on the spline analysis is provided in Supplemental Material. A representation of the linear relationship of PM2.5 and the HRs from proportional hazard regression models are provided, where the PM2.5 concentration of 5.7 µg/m3 was used as the reference. Distribution histograms of PM2.5 are included in the background of these graphs.
PAF is presented as a measure of the proportion of the outcome in the population attributable to PM2.5 exposure above the EPA standard of 12 μg/m3. Further details on calculation of PAF and attributable burden of disease are provided in Supplemental Material.
Missing data were not imputed. In analyses, a 95% CI of an HR that does not include unity was considered statistically significant. In all analyses, a P value ≤0.05 was considered statistically significant. All statistical analyses were done using SAS Enterprise Guide version 7.1 and SAS 9.4 (SAS Institute, Cary, NC).
Sensitivity Analyses
To test robustness of study findings, we undertook a number of sensitivity analyses as described in Supplemental Material.
Disclosures
None.
Acknowledgments
Support for Veterans Affairs (VA)/Centers for Medicare and Medicaid Services data are provided by the Department of Veterans Affairs, Veterans Health Administration, Office of Research and Development, Health Services Research and Development, and VA Information Resource Center (Project Number/Data Use Agreement ID Al-Aly-01).
This work was funded by a grant from the US Department of Veterans Affairs (to Z.A.-A.).
The funders of this study had no role in study design, collection, analysis, and interpretation of data, writing the report, and the decision to submit the report for publication. The contents do not represent the views of the US Department of Veterans Affairs or the US Government.
Footnotes
Published online ahead of print. Publication date available at www.jasn.org.
This article contains supplemental material online at http://jasn.asnjournals.org/lookup/suppl/doi:10.1681/ASN.2017030253/-/DCSupplemental.
- Copyright © 2018 by the American Society of Nephrology