Research December 2016 | Volume 124 | Issue 12
Individual and Neighborhood Socioeconomic Status and the Association between Air Pollution and Cardiovascular Disease
Gloria C. Chi,1 Anjum Hajat,2 Chloe E. Bird,3 Mark R. Cullen,4 Beth Ann Griffin,5 Kristin A. Miller,1 Regina A. Shih,5 Marcia L. Stefanick,6 Sverre Vedal,2 Eric A. Whitsel,7,8 and Joel D. Kaufman1,2
PDF Version (370 KB)
Background: Long-term fine particulate matter (PM2.5) exposure is linked with cardiovascular disease, and disadvantaged status may increase susceptibility to air pollution-related health effects. In addition, there are concerns that this association may be partially explained by confounding by socioeconomic status (SES).
Objectives: We examined the roles that individual- and neighborhood-level SES (NSES) play in the association between PM2.5 exposure and cardiovascular disease.
Methods: The study population comprised 51,754 postmenopausal women from the Women’s Health Initiative Observational Study. PM2.5 concentrations were predicted at participant residences using fine-scale regionalized universal kriging models. We assessed individual-level SES and NSES (Census-tract level) across several SES domains including education, occupation, and income/wealth, as well as through an NSES score, which captures several important dimensions of SES. Cox proportional-hazards regression adjusted for SES factors and other covariates to determine the risk of a first cardiovascular event.
Results: A 5 μg/m3 higher exposure to PM2.5 was associated with a 13% increased risk of cardiovascular event [hazard ratio (HR) 1.13; 95% confidence interval (CI): 1.02, 1.26]. Adjustment for SES factors did not meaningfully affect the risk estimate. Higher risk estimates were observed among participants living in low-SES neighborhoods. The most and least disadvantaged quartiles of the NSES score had HRs of 1.39 (95% CI: 1.21, 1.61) and 0.90 (95% CI: 0.72, 1.07), respectively.
Conclusions: Women with lower NSES may be more susceptible to air pollution-related health effects. The association between air pollution and cardiovascular disease was not explained by confounding from individual-level SES or NSES.
Citation: Chi GC, Hajat A, Bird CE, Cullen MR, Griffin BA, Miller KA, Shih RA, Stefanick ML, Vedal S, Whitsel EA, Kaufman JD. 2016. Individual and neighborhood socioeconomic status and the association between air pollution and cardiovascular disease. Environ Health Perspect 124:1840–1847; http://dx.doi.org/10.1289/EHP199
Address correspondence to G.C. Chi, 1959 NE Pacific St., Box 357236, Seattle, WA 98195 USA. Telephone: (626) 872-3007. E-mail: email@example.com
This work was supported by grants from the U.S. Environmental Protection Agency (EPA, grants RD831697 and RD-83479601-0 to the University of Washington), the National Institutes of Health (NIH), U.S. Department of Health and Human Services (DHHS), National Institute of Environmental Health Sciences (grants F31ES025475, P50ES015915, K24ES013195, and 1K99ES023498-01 to the University of Washington and grant R01-ES012238 to the University of North Carolina at Chapel Hill), and the NIH/DHHS/National Heart, Lung, and Blood Institute (grant R01-HL084425 to the RAND Corporation). The WHI (Women’s Health Initiative) program is funded by the National Heart, Lung, and Blood Institute (grants HSN268201100046C, HHSN268201100001C, HHSN268201100002C, HSN268201100003C, HHSN268201100004C, and HHSN271201100004C).
The contents of this publication are solely the responsibility of the grantee and do not necessarily represent the official views of the U.S. EPA or the NIH. Further, the U.S. EPA does not endorse the purchase of any commercial products or services mentioned in the publication.
The authors declare they have no actual or potential competing financial interests.
Received: 14 June 2015
Revised: 3 November 2015
Accepted: 19 April 2016
Published: 3 May 2016
Note to readers with disabilities: EHP strives to ensure that all journal content is accessible to all readers. However, some figures and Supplemental Material published in EHP articles may not conform to 508 standards due to the complexity of the information being presented. If you need assistance accessing journal content, please contact firstname.lastname@example.org. Our staff will work with you to assess and meet your accessibility needs within 3 working days.
- Supplemental Material PDF (816 KB)
Note to readers with disabilities: EHP has provided a 508-conformant table of contents summarizing the Supplemental Material for this article (see below) so readers with disabilities may determine whether they wish to access the full, nonconformant Supplemental Material. If you need assistance accessing journal content, please contact email@example.com. Our staff will work with you to assess and meet your accessibility needs within 3 working days.
Related EHP Articles
Comment on “Individual and Neighborhood Socioeconomic Status and the Association between Air Pollution and Cardiovascular Disease”
Response to “Comment on ‘Individual and Neighborhood Socioeconomic Status and the Association between Air Pollution and Cardiovascular Disease’”
Large population studies have routinely demonstrated that exposure to air pollution is associated with increased risk of cardiovascular morbidity and mortality (Brook et al. 2010). Low-socioeconomic status (SES) has also consistently been identified as a risk factor for cardiovascular disease (CVD) (Elo 2009). In addition, SES putatively co-varies with the spatial distribution of air pollution (Hajat et al. 2015). In this study, we set out to address two different questions regarding the role of SES in the air pollution–CVD relationship.
Our first aim was to answer a substantive question of whether individuals with low SES are more susceptible to the effects of air pollution on CVD. This question is crucial in informing air quality standards sufficient to protect the health of sensitive groups. We addressed this question by testing whether individual or neighborhood SES are effect modifiers of the air pollution–CVD relationship. Extant literature provides mixed evidence of effect modification by SES on the association between air pollution and health outcomes, including CVD. Low-SES individuals may be more susceptible to adverse effects of air pollution because they have poorer health resulting from reduced material resources, have higher psychosocial stress, and exhibit more individuals risk factors such as unhealthy behaviors and lifestyles (Elo 2009; O’Neill et al. 2003).
Our second question addresses the important methodological problem of whether SES confounds the association between air pollution and CVD. Confounding by SES is particularly concerning because low SES is a strong risk factor for CVD (Elo 2009) and also co-varies spatially with air pollution. Some North American studies have reported that communities with low SES are more likely to be exposed to higher concentrations of air pollution (Hajat et al. 2015), whereas, European research has been mixed (Hajat et al. 2015). The direction of confounding may depend on how SES co-varies with air pollution in the study population.
Epidemiological studies of air pollution health effects commonly include some measures of SES, such as individual-level education or income, but few incorporate multiple levels of SES. Both individual-level SES and neighborhood-level SES (NSES) are independently related to PM2.5 (Chaix et al. 2006; Hajat et al. 2013). NSES may have greater impact on estimated associations of air pollution on mortality than does individual-level SES (Næss et al. 2007). Not controlling for both levels of SES may lead to potential residual confounding. In addition, individual-level and contextual NSES may increase susceptibility to air pollution-related health outcomes via different processes. For instance, individual poor health status (e.g., diabetes and obesity) may act in synergy with air pollution exposure to promote CVD (O’Neill et al. 2003). On the neighborhood level, contextual factors such as lower-housing stock may increase co-exposure of other pollutants to induce greater susceptibility to harmful effects of air pollution (O’Neill et al. 2003).
SES has been characterized as a multidimensional concept often operationalized by measuring three specific domains: education, occupation, and income/wealth—each having different effects at various times in the life course (Elo 2009). For example, some have proposed that education affects health by promoting accumulation of knowledge regarding health-promoting behaviors and technologies and by improving decision-making and problem-solving skills (Elo 2009). As for occupational class, those working in higher occupational class positions tend to have less exposure to potentially harmful chemicals and pollutants present in the workplace (O’Neill et al. 2003). Lastly, income and wealth are financial resources that enable access to health-generating resources (such as good quality housing in safe neighborhoods) and access to higher quality health care (Elo 2009). The three domains are also interrelated because educational attainment influences subsequent occupation and income. Thus, epidemiological studies involving SES should include measures that reflect various domains of SES.
Unlike many data sets, the Women’s Health Initiative (WHI) Observational Study has available a unique set of SES indicators on both the individual and neighborhood levels that span the three domains of SES: education, occupation, and income/wealth. This is a substantial methodological improvement in measuring SES. In this article, we examine the roles that individual-level SES and NSES play in the association between fine particulate matter (PM2.5) and incident cardiovascular events in the longitudinal WHI Observational Study.
The WHI Observational Study enrolled 93,676 postmenopausal women 50–79 years old from 40 centers throughout the United States between 1993 and 1998. Eligible women included those who provided written informed consent, who planned to stay in the area, and who were free of conditions that might interfere with follow-up. The study design and participant characteristics were described previously (Women’s Health Initiative Study Group 1998; Hays et al. 2003; Langer et al. 2003). In the observational study, 93,676 women were screened at baseline to obtain information on demographics, lifestyle, medical history, cardiovascular risk factors, and anthropometric and blood pressure measurements (Langer et al. 2003). Annual mailed questionnaires collected updates on health outcomes. In this analysis, participants were followed from baseline until the end of follow-up of the main cohort in September 2005 (mean follow-up 7.6 years). Participants of this study were postmenopausal women > 50 years old at baseline who were better educated than women of the same age in the U.S. general population, limiting our ability to generalize our results to the overall population.
The current analysis was restricted to participants free of CVD (myocardial infarction, congestive heart failure, coronary revascularization, and stroke) at baseline with at least one PM2.5 prediction resolved to the street (i.e., higher resolution than ZIP code centroid) over the study period. Of the 93,676 participants, 18,576 had CVD and 2,006 had missing CVD status at baseline and were excluded. A further 17,115 participants were excluded due to missing covariates. Of the remaining 55,979 women, 1,493 had completely missing PM2.5 predictions and 2,732 had geocodes not resolved to the street. These categories were not mutually exclusive, and our final analytic sample included 51,754 women. Those excluded had lower SES and were more likely to be nonwhite, smokers, and diabetics. However, the analytic sample is generally representative of the baseline sample of women who were free of CVD, except for having more white participants than the baseline sample (86.3% vs. 83.6%). All standardized mean differences comparing important measures (including exposure and CVD incidence) of complete cases to the original sample of eligible participants were < 0.1 (see Table S1). Thus, all important measures were well-balanced between the complete cases and original sample. The proportion of missing data ranged from 0.05% to 10.6%.
All known participant home addresses over the follow-up period were geocoded. For each address, the point-specific annual average PM2.5 concentration was predicted using U.S. Environmental Protection Agency’s (EPA) Air Quality System (AQS) and Interagency Monitoring of Protected Visual Environments (IMPROVE) monitoring data for the year 2000 and used to represent ambient PM2.5 concentrations at that address over the entire follow-up. The year 2000 was selected because it represented an early year of complete national PM2.5 monitoring and a representative year of the follow-up period. Relative concentrations of particulate pollution were largely consistent for study locations throughout the study period (Miller et al. 2007). In addition, analysis from the American Cancer Society’s study demonstrated that PM2.5 was strongly correlated between sites during a 20-year period and the hazard function was not time-dependent—suggesting that fine particulate matter measured at any point over the study period is a reasonable surrogate for long-term particulate matter exposure (Abrahamowicz et al. 2003). Likelihood-based ambient point-specific PM2.5 predictions at participant residences were obtained using a regionalized national universal kriging model that included over 200 geographic covariates reduced via partial least squares techniques (Sampson et al. 2013). This approach resulted in a high level of cross-validated accuracy of prediction with an overall R2 of 0.88.
To calculate time-varying PM2.5 exposure, data were split on each time that a first cardiovascular event occurred, generating multiple records for each address for each participant. For each record, PM2.5 exposure was calculated as an average of the current and all previous PM2.5 predictions weighted by time spent at each residence. Splitting the data allowed us to calculate exposures that incorporated information before a cardiovascular event or censorship. The exposure was only time-varying in the sense that it incorporated residential history, but not time-varying in calendar time as all predictions were estimated for the year 2000.
The outcome of interest was time from enrollment until first cardiovascular event, which included myocardial infarction, stroke, death from coronary heart disease, and death from cerebrovascular disease. The WHI identified CVD outcomes through annual follow-up questionnaires to participants. Outcomes were ascertained via local and central review and adjudication of medical records by trained physicians (Curb et al. 2003). Deaths were ascertained via proxy reports and data linkage with the National Death Index of the National Center for Health Statistics. Physician adjudicators reviewed all available records for deaths including hospitalization records, autopsy records, and death certificate diagnoses (Curb et al. 2003). See Supplemental Material, “Women’s Health Initiative Classification Criteria for Cardiovascular Disease Events,” for further details about WHI criteria for classification of CVD events. Institutional review boards at the University of Washington and the Fred Hutchinson Cancer Research Center approved the study.
We assessed three distinct SES domains: education, occupation, and income/wealth, which have varied effects on health (Elo 2009). Individual-level SES characteristics were obtained from the baseline questionnaire and included education, family income, and occupation. We included four categories for education (less than high school, high school/GED/trade school, some college/associate degree, and bachelor’s degree or higher), five categories for family income (< $20,000; $20,000–$34,999; $35,000–$49,999; $50,000–$74,999; and ≥ $75,000), and four categories for occupation (managerial/professional, technical/sales/administrative, service/labor, and homemaker).
On the neighborhood level, we had more measures available and included the corresponding measures of percent of adults 25 years and older with a high school degree (education), percent of civilian population 16 years and older with professional/managerial/executive occupations (occupation), median family income (income), and percent of families above the poverty line (income). Data from the 2000 Census was used to assess baseline NSES at the tract level, a unit of geography small enough to be considered a reasonable proxy for neighborhood (Soobader et al. 2006). We also included median home value of owner-occupied housing units as a surrogate of wealth. Income and wealth are not surrogates for one another and both may influence health (Braveman et al. 2005). For instance, wealth may buffer consequences of temporary income loss (e.g., due to unemployment). More importantly, however, wealth may be a better indicator of economic SES among older adults because a) income and occupation become less important for retired individuals and b) accumulated financial assets such as home ownerships become more significant (Pollack et al. 2007). An individual-level measure of wealth was not available in this data set. Finally, we had available an NSES score that was previously related to incident coronary heart disease in this cohort (Bird et al. 2009). The NSES score is a composite measure of six Census tract-level variables that was created from a confirmatory factor analysis examining 12 theoretically relevant measures and was only available in metropolitan statistical areas (Dubowitz et al. 2008). This index was composed of a) percent of adults 25 years and older with less than a high school education, b) percent male unemployment, c) percent of households with income below the poverty line, d) percent of households receiving public assistance, e) percent of households with children headed only by a female, and f) median household income. Values of the NSES score during intercensal years were interpolated, and participants were assigned baseline values based on their year of enrollment. Higher values on the score indicate less deprivation.
Individual-level SES indicators had weak to moderate correlations with each other and with NSES indicators (correlation coefficients range from 0.14 to 0.36) (see Table S2). NSES indicators exhibited stronger correlation with each other (correlation coefficients from 0.45 to 0.85); even so, NSES indicators represent distinct domains of SES that putatively affect health via distinct mechanisms and pathways (Elo 2009).
The relationship between long-term annual average PM2.5 exposure and time from enrollment until incident cardiovascular events was assessed using Cox proportional hazards models. The following baseline characteristics were controlled for as potential confounders: age, race/ethnicity, diabetes, hypertension, hypercholesterolemia, smoking (smoking status, cigarettes per day, years smoked), and body mass index. Analyses were stratified by 5-year age categories, body mass index (five categories), and diabetes status for a more thorough adjustment. Race/ethnicity was condensed into a binary variable for white not of Hispanic origin and a group including American Indian/Alaskan Native, Asian/Pacific Islander, black, Hispanic, and unknown race/ethnicity due to the small numbers of participants in the latter category.
Although the use of explanatory variables at both the individual and neighborhood levels suggests a multilevel approach, multilevel Cox regression models are often computationally intensive and cumbersome (Goldstein 1995). Therefore, this study utilized the more tractable marginal method which uses traditional estimation. To obtain estimates of standard errors and p-values unbiased by geographic clustering of individuals, we adjusted the variances of these coefficients using a sandwich estimator (Lee et al. 1992; Lin 1994).
Effect modification by each individual-level SES and NSES indicator was investigated by fitting multiplicative interaction terms for different levels of the SES variable with PM2.5. Separate interaction models were fit for each individual-level SES and NSES variable. Models adjusted for all individual-level SES and NSES variables in addition to other baseline covariates. The model for the composite NSES score included all individual-level SES variables and adjustment covariates but no other NSES variables. Joint tests were conducted to simultaneously test all interaction terms for the SES indicator in question. The Benjamini–Hochberg method was used to control the false discovery rate at 5% (Benjamini and Hochberg 1995). Statistical analyses were performed using Stata (release 13; Stata Statistical Software).
To observe potential confounding by SES variables, we fitted separate models for each SES variable (individual or neighborhood level) and adjusted for all non-SES covariates. Individual-level SES variables were included as categorical variables. All NSES variables were measured continuously, but split into quartiles and included as factor variables in the analyses to allow for non-linear relationships between these measures and time to CVD event. All SES models were adjusted separately (not sequentially). We then proceeded to fit three combinations of SES variables: all individual-level SES only, all NSES only, and both.
Hypertension may lie along the causal pathway between air pollution and CVD. In sensitivity analyses, hypertension was removed from models. We also evaluated whether after adjustment for SES (individual level first and then contextual) there was residual confounding from individual behavioral factors. This is pertinent to large cohort studies using administrative data that lack individual variables. To address the possibility of selection bias due to complete case analysis, missing values in SES variables and adjustment covariates were multiply imputed (see Supplemental Material, “Multiple Imputation”). However, these analyses were run using baseline PM2.5 instead of a time-weighted average PM2.5 due to issues of computational feasibility. In addition, cross-level interactions were explored, looking at the following categories: low SES in both levels, low SES in one level and high SES in the other, and high SES in both levels (see Supplemental Material, “Cross-level Interaction”).
Our analytic sample included 51,754 women with 387,840 women-years of follow-up. Mean age at enrollment was 63 years. Most participants were non-Hispanic whites (86.3%) and were never or past smokers (52.4% and 41.4%, respectively) (Table 1). In general, subject characteristics were similar across different categories of first PM2.5 prediction, although those in the highest exposure quartile tended to have fewer non-Hispanic whites and lower NSES (Tables 1 and 2). Those who experienced CVD events tended to have less education and lower income, were less likely to work in managerial or professional positions, and more likely to live in lower-NSES neighborhoods (see Table S3).
We observed 1,737 cardiovascular events. The number of events in each quartile of first available PM2.5 prediction is shown in Table 3. The highest number of events was observed in the highest quartile of PM2.5. The overall mean concentration of all PM2.5 observations was 12.7 μg/m3 (SD, 2.9; interquartile range, 4.1); the minimum was 2.2 μg/m3 and the maximum was 25.1 μg/m3. Figure S1 shows a scatterplot of first available PM2.5 predictions by NSES score with a locally weighted scatterplot smoothing curve. In general, areas with lower NSES tended to experience slightly higher levels of PM2.5.
Exposure to PM2.5 was significantly associated with risk of cardiovascular events. After adjustment for age, race/ethnicity, smoking, body mass index, diabetes, hypertension, and hypercholesterolemia, a 5 μg/m3 higher exposure to PM2.5 was associated with a 12% higher risk of cardiovascular event [hazard ratio (HR): 1.12; 95% confidence interval (CI): 1.00, 1.25; Table 4]. Further adjustment for individual-level SES or NSES (singly or combined) did not change the HR materially. In Table 4, each adjustment listed is separate and not sequentially related to the adjustment above it. For example, the individual income model only adjusted for individual income and all non-SES covariates. The fully adjusted model, which included all potential confounders and all individual-level SES and NSES variables except for the NSES score, had an HR of 1.13 (95% CI: 1.02, 1.25).
The associations of PM2.5 with CVD events by categories of individual SES and NSES variables are shown in Figures 1 and 2, respectively. None of the individual-level SES variables significantly modified the association between PM2.5 and CVD events. Although those with the lowest individual income (< $20,000) had an HR of 1.30 (95% CI: 1.12, 1.52), it was not significantly different from other income categories. There is evidence of statistically significant effect modification by the NSES score (2-sided p = 0.008), median home value (2-sided p < 0.001), and percentage above poverty (2-sided p = 0.013) after accounting for multiple comparisons. Those in the most disadvantaged quartile of the NSES score had an HR of 1.39 (95% CI: 1.21, 1.61); whereas, those in the least disadvantaged quartile had an HR of 0.90 (95% CI: 0.72, 1.07). Similarly, those in the lowest quartile of median home value had an HR of 1.40 (95% CI: 1.24, 1.58) compared to those in the highest quartile with an HR of 0.87 (95% CI: 0.77, 0.99). Furthermore, significant positive associations between PM2.5 exposure and CVD risk were observed in the most disadvantaged quartiles of all NSES variables examined. HRs tended to decrease as NSES increased, and this trend is consistent across categories for multiple neighborhood-level indicators, unlike for individual-level SES indicators.
Figure 1 – Estimated hazard ratios and 95% confidence intervals for time to first cardiovascular event associated with 5 μg/m3 higher exposure to PM2.5 according to levels of individual socioeconomic status (SES) and p-values for interactions. Models adjusted for age, race/ethnicity, smoking, body mass index, diabetes, hypertension, hypercholesterolemia, and all individual-level SES and neighborhood-level SES indicators except for the neighborhood SES score.
Figure 2 – Estimated hazard ratios and 95% confidence intervals for time to first cardiovascular event associated with 5 μg/m3 higher exposure to PM2.5 according to levels of neighborhood socioeconomic status (NSES) and p-values for interactions. Models adjusted for age, race/ethnicity, smoking, body mass index, diabetes, hypertension, hypercholesterolemia, and all individual-level socioeconomic status and NSES indicators except for the NSES score. All NSES variables were grouped into quartiles, ranging from lowest NSES (most deprived) to highest NSES (least deprived). The NSES score model adjusted for individual-level SES indicators but no other NSES indicators.
In sensitivity analyses, removal of hypertension did not change our results materially (see Table S4). In the model only adjusted for age, the estimated HR was 1.13 (95% CI: 1.01, 1.28). Additional control for SES indicators in age-only models did not change HRs meaningfully. There was no evidence of confounding by SES indicators even without controlling for individual risk factors and no indication of residual confounding by individual factors. Multiple imputation of missing covariates, including SES variables, did not change the results materially (see Tables S5 and S6). Looking at cross-level interactions, having both low individual education and low NSES (any characteristic) did not confer greater vulnerability than having low SES on only one level (see Figure S2). However, there is evidence that having low individual income and low NSES (any characteristic) conferred greater risk of CVD than having high SES on at least one level.
Results corroborate previous studies that exposure to long-term PM2.5 is a risk factor for CVD, and this association cannot be explained by confounding by individual-level SES or NSES. Furthermore, the association was stronger for women residing in lower-SES neighborhoods.
Our results of effect modification by NSES are consistent with the hypothesis that those with low SES may be disproportionately affected by the adverse health effects of air pollution. Researchers documented that individuals with low SES and racial minorities experience higher exposure to air pollution (Hajat et al. 2015) and also suffer from worse health outcomes resulting from poverty and psychosocial stress in poor communities (Diez Roux et al. 2004). The combination of greater exposure to air pollution, poorer health, and fewer resources to cope with the effects of air pollution may result in increased susceptibility to air pollution-related health outcomes (O’Neill et al. 2003). We see stronger effect modification for neighborhood property values compared to median household income, which might stem from property values being a better reflection of SES for older individuals. In addition, it is possible that property values may be higher in communities of more owners (vs. renters), and these owners may be more invested in the long term, which could contribute to contextual factors such as neighborhood stability or investment, and other social processes not captured by income or administrative data.
The lack of evidence for effect modification by individual-level SES suggests that neighborhood-level processes may increase susceptibility to air pollution-related CVD. First, NSES is on the same spatial scale as air pollution, and empirical evidence shows that the association between individual-level SES and PM2.5 is often weaker than that observed between NSES and air pollution (Hajat et al. 2013). Secondly, macro-level contextual factors, such as racial-residential segregation are hypothesized to differentially distribute exposures to environmental hazards and to concentrate poverty (Gee and Payne-Sturges 2004; Morello-Frosch and Lopez 2006). Disadvantaged neighborhood environments may be working through the stress pathway to impact health (Diez Roux and Mair 2010) making residents more susceptible to the health effects of PM2.5.
Exposure measurement error, particulate matter infiltration, dose reduction, and subject time activity patterns may differ according to individual-level SES or NSES and could explain part of our findings of effect modification by NSES. Higher rates of PM2.5 infiltration have been reported for lower-SES individuals (Hystad et al. 2009), which may be explained by decreased use of air conditioning and older and poorer housing quality among low-SES individuals. Thus, using ambient exposures would systematically underestimate true exposures for lower-SES persons compared to those with higher SES, which would be consistent with the direction of effect modification observed in this study. Furthermore, the health-motivated individuals among those with more resources may use their resources not only to seek cleaner residential areas but also to reduce background exposures (e.g., by better air conditioning). Thus, some of the effect modification may actually represent true dose reduction in those with high NSES.
Our findings of no confounding by individual-level SES in this cohort are consistent with studies reporting small changes in relative risk estimates after adjustment for SES including education and income (Brochu et al. 2011; Dockery et al. 1993; Pope et al. 2002). However, research in Canada suggested that NSES positively confounded the relationship between particulate air pollution and mortality, where adjustment for several NSES variables changed risk estimates more than 10% (Jerrett et al. 2005). The WHI cohort has higher SES relative to the United States as a whole, and relatively small SES variability in our data could explain the lack of confounding by SES in our results. In addition, while individual-level SES is associated with CVD outcomes, it is not strongly associated with exposure, and the converse is true for NSES (NSES is associated with exposure but not strongly associated with CVD outcomes)—hence, neither served as strong confounders. Either individual-level SES or NSES may still be an important confounder in populations where associations between SES and air pollution and SES and CVD are large.
NSES may also lie along the causal pathway between air pollution and CVD. For example, poor air quality due to a polluting facility and increasing traffic may change the attractiveness of a neighborhood, causing higher SES individuals to move away, lower-SES individuals to move in, and/or home values to decline. The resulting lower NSES of the neighborhood may cause negative changes in the neighborhood’s social, physical, and built environments that could result in adverse health effects for residents. If so, NSES would be in the causal pathway and should be dealt with in an analytically appropriate manner. Understanding the directionality of NSES and air quality is difficult especially in a multisite study of this nature where different processes are likely occurring in different places (Saha and Mohai 2005).
This study is consistent with the findings in the analysis by Miller et al. (HR = 1.24; 95% CI: 1.09, 1.4; per 10 μg/m3 higher PM2.5) (Miller et al. 2007) and analysis in a Health Effects Institute report (HR = 1.25; 95% CI: 1.09, 1.44; per 10 μg/m3 higher PM2.5) (Vedal et al. 2013) which used the same cohort of women. Our fully adjusted HR is 1.29 (95% CI: 1.04, 1.59) when scaled to a 10 μg/m3 higher exposure to PM2.5. However, the current analysis is primarily focused (unlike the prior analyses) on the very important methodological problem of disentangling air pollution exposures from putatively co-varying socioeconomic factors. The first study (Miller et al. 2007) assessed NSES in a sensitivity analysis and found no evidence of confounding by U.S. Census-derived measures aggregated to ZIP code level, but did not assess NSES as a potential effect modifier. The current study also uses improved exposure assessment that is resolved to participant addresses and incorporates residential history. The first study assigned exposures based on nearest monitor to participant homes. The second utilized baseline PM2.5 predictions at geocoded addresses, but included sensitivity analyses that incorporated exposures based on residential location in the 1 or 2 years before event or death (or corresponding years in subjects with no events). The current study uses a time-weighted average PM2.5 that incorporates all residential history. Finally, the current study had 2 additional years of follow-up compared to the first study.
There are several limitations to this study. First, selection bias may arise from conducting a complete case analysis. However, important measures including exposure and CVD incidence were well-balanced between complete cases and eligible participants in the original sample (see Table S1), indicating that our analytic sample is generally representative of the baseline sample of women who were CVD free. There was also no individual measure of wealth available, and we were not able to assess either effect modification or confounding by individual wealth. In addition, sensitivity analyses that multiply imputed all SES variables and adjustment covariates were not materially different from the main analysis. Another issue is unobserved confounding by self-selection into neighborhoods. However, our estimates are robust to adjustment for many demographic, socioeconomic, and health characteristics that may correlate with self-selection. We also do not have measures of indoor air pollution, which may better reflect true exposures in an older population which spends more time indoors. Future studies are needed to assess the effect of this measurement error. Our SES measures were not adjusted for variation in cost of living and housing in different regions, which may lead to measurement error in our SES assessment in a national cohort. The effect of this error is likely to be location- and population-specific. However, while developing the NSES score, the authors conducted sensitivity analysis to adjust for differences in cost of living, and it did not have an effect on the NSES score. Further research is needed to assess the importance of adjusting for cost of living in the measurement of individual-level SES. Future work is also needed to test for sex differences and the impact of NSES on the association between PM2.5 and CVD events in samples that are more representative of the United States.
This study has several strengths. First, analyses were conducted using a large sample size and a long follow-up time. Second, outcomes were adjudicated based on protocol-based review of medical records, thereby reducing outcome misclassification. Third, we were able to resolve PM2.5 exposures to the level of the individual’s residence based on geocodes and a state-of-the-art fine-scale modeling framework, reducing exposure misclassification. In addition, this study is among the few to investigate the roles of both individual-level SES and NSES in different domains on the association between air pollution and cardiovascular disease, which is an important methodological improvement on prior attempts to measure SES. Finally, the study examines a range of women of moderate income residing across a wide range of NSES, which gives considerable ability to assess the impact of NSES. Many past studies did not have data on both.
We investigated the role that SES plays in the association between PM2.5 and CVD. We found that individual-level SES and NSES did not confound the positive association between PM2.5 and CVD in this cohort. Furthermore, risk estimates were higher for women living in more disadvantaged neighborhoods. Our findings contribute to the understanding of confounding by SES in air pollution health effects research and support an evolving understanding of the synergistic adverse effects of air pollution and socioeconomic factors.
Abrahamowicz M, Schopflocher T, Leffondré K, du Berger R, Krewski D. 2003. Flexible modeling of exposure-response relationship between long-term average levels of particulate air pollution and mortality in the American Cancer Society study. J Toxicol Environ Health A 66:1625–1654, doi: 10.1080/15287390306426.
Bird CE, Shih RA, Eibner C, Griffin BA, Slaughter ME, Whitsel EA, et al. 2009. Neighborhood socioeconomic status and incident coronary heart disease among women [Abstact]. J Gen Intern Med 24(suppl 1):S127.
Braveman PA, Cubbin C, Egerter S, Chideya S, Marchi KS, Metzler M, et al. 2005. Socioeconomic status in health research: one size does not fit all. JAMA 294:2879–2888, doi: 10.1001/jama.294.22.2879.
Brochu PJ, Yanosky JD, Paciorek CJ, Schwartz J, Chen JT, Herrick RF, et al. 2011. Particulate air pollution and socioeconomic position in rural and urban areas of the Northeastern United States. Am J Public Health 101(suppl 1):S224–S230, doi: 10.2105/AJPH.2011.300232.
Brook RD, Rajagopalan S, Pope CA III, Brook JR, Bhatnagar A, Diez-Roux AV, et al. 2010. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation 121:2331–2378, doi: 10.1161/CIR.0b013e3181dbece1.
Chaix B, Gustafsson S, Jerrett M, Kristersson H, Lithman T, Boalt Å, et al. 2006. Children’s exposure to nitrogen dioxide in Sweden: investigating environmental injustice in an egalitarian country. J Epidemiol Community Health 60:234–241.
Curb JD, McTiernan A, Heckbert SR, Kooperberg C, Stanford J, Nevitt M, et al. 2003. Outcomes ascertainment and adjudication methods in the Women’s Health Initiative. Ann Epidemiol 13(9 Suppl):S122–S128.
Diez Roux AV, Borrell LN, Haan M, Jackson SA, Schultz R. 2004. Neighbourhood environments and mortality in an elderly cohort: results from the cardiovascular health study. J Epidemiol Community Health 58:917–923.
Diez Roux AV, Mair C. 2010. Neighborhoods and health. Ann NY Acad Sci 1186:125–145, doi: 10.1111/j.1749-6632.2009.05333.x.
Dockery DW, Pope CA III, Xu X, Spengler JD, Ware JH, Fay ME, et al. 1993. An association between air pollution and mortality in six U.S. cities. N Engl J Med 329:1753–1759, doi: 10.1056/NEJM199312093292401.
Dubowitz T, Heron M, Bird CE, Lurie N, Finch BK, Basurto-Dávila R, et al. 2008. Neighborhood socioeconomic status and fruit and vegetable intake among whites, blacks, and Mexican Americans in the United States. Am J Clin Nutr 87:1883–1891.
Gee GC, Payne-Sturges DC. 2004. Environmental health disparities: a framework integrating psychosocial and environmental concepts. Environ Health Perspect 112:1645–1653, doi: 10.1289/ehp.7074.
Hajat A, Diez-Roux AV, Adar SD, Auchincloss AH, Lovasi GS, O’Neill MS, et al. 2013. Air pollution and individual and neighborhood socioeconomic status: evidence from the Multi-Ethnic Study of Atherosclerosis (MESA). Environ Health Perspect 121:1325–1333, doi: 10.1289/ehp.1206337.
Hajat A, Hsia C, O’Neill MS. 2015. Socioeconomic disparities and air pollution exposure: a global review. Curr Environ Health Rep 2:440–450, doi: 10.1007/s40572-015-0069-5.
Langer RD, White E, Lewis CE, Kotchen JM, Hendrix SL, Trevisan M. 2003. The Women’s Health Initiative Observational Study: baseline characteristics of participants and reliability of baseline measures. Ann Epidemiol 13(9 suppl):S107–S121.
Lee EW, Wei LJ, Amato DA, Leurgans S. 1992. Cox-type regression analysis for large numbers of small groups of correlated failure time observations. In: Survival Analysis: State of the Art (Klein JP, Goel PK, eds). Dordrecht, Netherlands:Springer Netherlands, 237–247.
Miller KA, Siscovick DS, Sheppard L, Shepherd K, Sullivan JH, Anderson GL, et al. 2007. Long-term exposure to air pollution and incidence of cardiovascular events in women. N Engl J Med 356:447–458, doi: 10.1056/NEJMoa054409.
Morello-Frosch R, Lopez R. 2006. The riskscape and the color line: examining the role of segregation in environmental health disparities. Environ Res 102:181–196, doi: 10.1016/j.envres.2006.05.007.
O’Neill MS, Jerrett M, Kawachi I, Levy JI, Cohen AJ, Gouveia N, et al. 2003. Health, wealth, and air pollution: advancing theory and methods. Environ Health Perspect 111:1861–1870, doi: 10.1289/ehp.6334.
Sampson PD, Richards M, Szpiro AA, Bergen S, Sheppard L, Larson TV, et al. 2013. A regionalized national universal kriging model using Partial Least Squares regression for estimating annual PM2.5 concentrations in epidemiology. Atmos Environ 75:383–392.
Vedal S, Campen MJ, McDonald JD, Larson TV, Sampson PD, Sheppard L, et al. 2013. National Particle Component Toxicity (NPACT) initiative report on cardiovascular effects. Res Rep Health Eff Inst 178:5–8.