- Open Access
- Open Peer Review
Longitudinal associations of long-term exposure to ultrafine particles with blood pressure and systemic inflammation in Puerto Rican adults
Environmental Healthvolume 17, Article number: 33 (2018)
Few longitudinal studies have examined the association between ultrafine particulate matter (UFP, particles < 0.1 μm aerodynamic diameter) exposure and cardiovascular disease (CVD) risk factors. We used data from 791 adults participating in the longitudinal Boston Puerto Rican Health Study (Massachusetts, USA) between 2004 and 2015 to assess whether UFP exposure was associated with blood pressure and high sensitivity C-reactive protein (hsCRP, a biomarker of systemic inflammation).
Residential annual average UFP exposure (measured as particle number concentration, PNC) was assigned using a model accounting for spatial and temporal trends. We also adjusted PNC values for participants’ inhalation rate to obtain the particle inhalation rate (PIR) as a secondary exposure measure. Multilevel linear models with a random intercept for each participant were used to examine the association of UFP with blood pressure and hsCRP.
Overall, in adjusted models, an inter-quartile range increase in PNC was associated with increased hsCRP (β = 6.8; 95% CI = − 0.3, 14.0%) but not with increased systolic blood pressure (β = 0.96; 95% CI = − 0.33, 2.25 mmHg), pulse pressure (β = 0.70; 95% CI = − 0.27, 1.67 mmHg), or diastolic blood pressure (β = 0.55; 95% CI = − 0.20, 1.30 mmHg). There were generally stronger positive associations among women and never smokers. Among men, there were inverse associations of PNC with systolic blood pressure and pulse pressure. In contrast to the primary findings, an inter-quartile range increase in the PIR was positively associated with systolic blood pressure (β = 1.03; 95% CI = 0.00, 2.06 mmHg) and diastolic blood pressure (β = 1.01; 95% CI = 0.36, 1.66 mmHg), but not with pulse pressure or hsCRP.
We observed that exposure to PNC was associated with increases in measures of CVD risk markers, especially among certain sub-populations. The exploratory PIR exposure metric should be further developed.
Long-term exposure to fine particulate matter (PM2.5, < 2.5 μm aerodynamic diameter) has been associated with cardiovascular disease (CVD) risk factors, such as increased blood pressure (BP) and concentrations of biomarkers of systemic inflammation, as well as increased risk of hypertension [1,2,3,4,5,6]. Less is known, however, about the impact of the smallest size fraction of PM, ultrafine PM (UFP, < 0.1 μm aerodynamic diameter), on these indicators. Compared to larger size fractions, UFP has a larger total deposition fraction, can penetrate more deeply into the lungs, has greater total surface area with which to interact with epithelial cells, is more likely to cross biological barriers, and can induce oxidative stress more readily [7,8,9,10,11,12,13,14,15].
A major source of UFP in urban areas is motor vehicle exhaust and there is high spatial and temporal variability of UFP compared to other sizes of particulate matter, especially near roadways . Due in part to the challenges of modeling UFP [17,18,19,20,21,22], there is little epidemiological literature on the health consequences of long-term UFP exposure. In a cross-sectional study of the health effects of annual average particle number concentration (PNC) among adults in the Community Assessment of Freeway Exposure and Health Study residing in several near-highway and urban background areas of the greater Boston area, we found positive associations with high sensitivity C-reactive protein (hsCRP), a biomarker of systemic inflammation . The few previously published longitudinal analyses of the long-term health effects of UFP considered UFP modeled with larger spatial resolutions (between 200 m and 4 km). One study found that UFP mass was associated with increased ischemic heart disease mortality . Other studies of long-term exposure to PNC have found that PNC was positively associated with sub-clinical markers of atherosclerosis and was inconsistently associated with biomarkers of inflammation and with respiratory outcomes [25,26,27].
Given the small number of studies on the health effects of long-term UFP exposure, we investigated the relationship between UFP and cardiovascular risk factors in the prospective Boston Puerto Rican Health Study (BPRHS) . In this population, proximity to traffic was previously associated with changes in hsCRP levels over two years . For the present study, our primary objective was to assess whether ambient residential UFP number concentrations at a fine spatial scale (≤20 m resolution) were associated with BP and hsCRP levels over six years. We also explored a novel exposure metric, the particle inhalation rate (PIR, particles inhaled/time), that may more closely approximate the biologically-relevant dose of UFP.
The BPRHS is a prospective cohort study of 1499 individuals designed to investigate the risk factors for chronic disease among Puerto Ricans living in eastern Massachusetts . Briefly, participants were recruited through door-to-door enumeration and through community approaches from census tracts in the greater Boston area with at least 10 Hispanics aged 45–75 years. Individuals were eligible for inclusion in the BPRHS if they were 45–75 years old at baseline, were able to answer questions in English or Spanish, and self-identified as being Puerto Rican. Participants were excluded if they had plans to move within two years or if they had low cognitive function as measured by the Mini Mental State Examination (scores ≤10). We restricted our analyses to those participants who lived within a 1000 m buffer of our air pollution monitoring area at any study visit (n = 791, Fig. 1).
This study was approved by the Institutional Review Boards at Tufts Medical Center, Northeastern University, and the University of Massachusetts Lowell. All participants provided written informed consent.
Participants were visited up to three times over approximately six years (visit one between 2004 and 2009, visit two between 2006 and 2011, and visit three between 2011 and 2015). The mean time between visit one and visit two was 2.2 years while the mean time between visit two and visit three was 4.1 years. Trained Spanish-English bilingual interviewers administered questionnaires in participants’ homes. Educational attainment was analyzed as a dichotomous variable (>8th grade/≤8th grade), based on the distribution of educational attainment in this population. Smoking status was assessed as current, former, or never (< 100 cigarettes smoked). Poverty status was determined by comparing participants’ total self-reported annual household income to the U.S. Census Bureau annual thresholds . Medication use (prescription and over-the-counter) was assessed by visual inspection of medications. Physical activity was assessed using a modified Paffenbarger questionnaire of the Harvard Alumni Activity Survey which has been previously validated in an elderly Puerto Rican population [31, 32]. Validated scales were used to assess psychological acculturation and perceived stress [33, 34]. Participants were considered depressed if they reported taking medication for depression or if they scored ≥16 on the Center for Epidemiology Studies Depression Scale [35,36,37].
Height and weight were measured in duplicate. Body mass index (BMI) was calculated as kg/m2. A trained interviewer measured participants’ seated systolic blood pressure (SBP) and diastolic blood pressure (DBP) with an electronic sphygmomanometer (Dinamap™ Model 8260, Critikon, Tampa, FL), in duplicate, three times. The second and third sets of readings were averaged . Participants were considered hypertensive if they had a SBP ≥140 mmHg, a DBP ≥90 mmHg, self-reported a diagnosis of hypertension, or if they were taking medication for hypertension. Pulse pressure (PP) was calculated as the difference between SBP and DBP.
A phlebotomist obtained blood samples after each study visit (> 81% had fasting samples at each visit; visits generally occurred in the morning). Cholesterol, triglycerides, serum glucose, and hsCRP concentrations were measured as described previously . We excluded the top 1% of hsCRP values (≥47 mg/L) since extremely high hsCRP values are associated with acute infections . Triglyceride and hsCRP concentrations were natural log transformed due to their skewed distributions. Participants were considered diabetic if their glucose concentration was ≥126 mg/dL, they were taking medications for diabetes, or they self-reported diabetes.
Geolocation of participants’ residences
Participants’ residential addresses at each study visit were geocoded using ArcMap . Addresses that could not be parcel-matched using the Boston parcel data were geocoded using Google Earth and publicly accessible site maps of housing developments. We geocoded 97% of participants’ residential locations at their second study visit (41% matched automatically) and 96% at their third study visit (39% matched automatically). A randomly selected subset of 12% of the geocoded locations for participants who moved was independently checked. All addresses were geocoded to the same parcel and the mean difference in position was less than 10 m.
Exposure to UFP was estimated using a land-use regression model based on a previously published model . Details on model development and validation are provided in the Additional file 1 (Part 1; see Table S1 for the full model). Briefly, the model was built using measured PNC from mobile and stationary platforms (Fig. 1) , meteorological data, and distances from specific roadways and bus routes (model-adjusted r2 = 0.37; see Additional file 1, Part 1) [42, 43]. The r2 is similar to other hourly PNC models specific to metropolitan Boston . The model predicted hourly ambient PNC at each participant’s address with ≤20 m resolution. For all PNC estimates before 2012, the hourly model was back-extrapolated using meteorological data collected during the time period of interest. This was possible because PNC regression models are largely stable over time  and models of traffic-related air pollutants, including PNC, perform reasonably well when back-extrapolated [46, 47]. Each participant was then assigned an ambient annual average PNC corresponding to the 365 days immediately preceding each of their study visits.
Calculation of particle inhalation rate
The particle inhalation rate (PIR, number of particles inhaled/h) was estimated for each participant as the product of the annual average PNC estimate (particles/L) and the hourly respiratory volume (tidal volume * breaths/h = L of air inhaled/h). We used published estimates for age- and sex-specific minute respiratory volume (L of air inhaled/min-kg) adjusted for weight and physical activity  together with data on how many hours per typical weekday and weekend day participants engaged in various levels of physical activity (lying down, sitting, light activity, moderate activity, and vigorous activity; for algorithm see Additional file 1, Part 2).
For each outcome (SBP, DBP, PP, ln(hsCRP)), we developed two different multilevel linear models to consider the longitudinal associations of 1) PNC and 2) PIR with the levels of the outcomes across study visits. All models controlled for age, included a random intercept for each participant, and used data from every completed visit. All modeling was performed in Stata v14 . To facilitate direct comparisons between the PNC and PIR models, we scaled results to the inter-quartile range (IQR; 4.6 thousand particles/mL and 6.2 billion particles inhaled/h, respectively).
We used a multi-stage process to select covariates. From an initial set of potential covariates identified by a literature review, we first included variables in the models if they were 1) associated with the outcome (p < 0.15) in bivariate analysis, and 2) either associated with the air pollution measure (p < 0.15) or changed the effect estimate for the air pollution measure by ≥10% in models that included only the air pollution measure and the potential covariate. Variables were retained if they were associated with the outcome in the multivariate model (p < 0.15) and if they did not introduce problems with collinearity based on variance inflation factors and correlation coefficients. We then assessed the effect of adding other variables that were considered potentially important based on the literature but had not met our initial inclusion criteria. If these variables were not associated with the outcome (p < 0.15) and did not materially change the effect estimates for PNC or PIR, they were excluded. Model residuals and model fit statistics were examined after the addition of each new covariate. Time varying predictors considered during the covariate selection process included age, BMI, high-density lipoprotein (HDL) cholesterol, low-density lipoprotein (LDL) cholesterol, ln(triglycerides), diabetes, hypertension medication, smoking, anxiety medication, perceived stress, psychological acculturation, marital status, and physical activity. Time invariant predictors considered included sex, educational attainment, and year of baseline visit. Year of baseline visit was included if it was a significant predictor to account for annual trends not captured within the PNC model. Physical activity and sex were not considered as confounders in the PIR models, as they were used in the calculation of the exposure. Variables assessed only at the third study visit (e.g., secondhand smoke exposure, family history of hypertension) were included only in sensitivity analyses and were assumed to be time-invariant. For each model, we checked collinearity and intra-class correlations. We also checked the normality and homoscedasticity of the residual errors.
Based on evidence from previous studies , we examined effect modification by sex, medication use, family history of hypertension (for BP), family history of CVD (for hsCRP), diabetes, smoking, employment status at baseline, physical activity, age, and BMI. To account for the high prevalence of baseline hypertension and cardiovascular disease, we also considered the effect of using the baseline measure of the outcome as a covariate in models that only used outcome data from the second and third study visits. Additionally, we conducted sensitivity analyses of the main models excluding participants who did not complete all three study visits, including participants with hsCRP values >99th percentile and excluding participants who, at baseline, self-reported at least one previous heart attack or stroke, had hypertension, had high baseline hsCRP (> 3 mg/L), or who died before their third study visit (n = 50).
Sixty-nine percent of participants were female and about half had attained more than an eighth grade education (Table 1). At each study visit, more than 70% of participants reported a household income below 120% of the federal poverty line  and only 22.1% of participants were employed at baseline. At baseline, 44.4% of participants reported never smoking, 23.7% were current smokers, and 31.9% were former smokers. Although the mean age at baseline was 57.1 years (standard deviation = 7.4; Table 1), 11.2% of participants had suffered at least one heart attack or stroke, 72.3% had hypertension, 46.0% had diabetes, 37.3% took antilipidemic (statin) medications, and 65.7% showed depressive symptomology. The participants included in the present analysis were similar to the larger BPRHS population .
Exposure distributions are summarized in Table 2 and spatial distributions are shown in Fig. 1. The distributions of both PNC and PIR were fairly stable across study visits (Table 2). PNC and PIR values were significantly correlated (r = 0.356, p < 0.001).
Overall, long-term exposure to PNC was not associated with SBP (β = 0.96; 95% CI = − 0.33, 2.25 mmHg per 4600 particles/mL; Fig. 2). Sex modified the effect of PNC on SBP. Among females, PNC was positively associated with SBP (β = 1.84; 95% CI = 0.21, 3.48 mmHg) while for males, PNC was inversely but not significantly associated with SBP (β = − 1.09; 95% CI = − 3.16, 0.98 mmHg). Additionally, PNC was positively associated with SBP among never smokers (β = 2.20; 95% CI = 0.04, 4.37 mmHg), but not among current or former smokers (Additional file 1: Table S2). When controlling for the baseline level of SBP, PNC was associated with changes in SBP (β = 1.66; 95% CI = 0.17, 3.14 mmHg). When using the PIR in place of PNC, exposure was associated with SBP (β = 1.03; 95% CI = 0.00, 2.06 mmHg per 6.2 billion particles inhaled/h). These associations were stronger among generally healthier participants (Additional file 1: Table S3). The associations between PIR and BP were attenuated when people with hypertension at baseline and, separately, people who dropped out before their third study visit were excluded (β = 1.12, 95% CI = − 0.49, 2.73 mmHg; β = 0.36, 95% CI = − 0.93, 1.64 mmHg, respectively).
There was no association overall between PNC and DBP (β = 0.55; 95% CI = − 0.20, 1.30 mmHg per 4600 particles/mL), except among never smokers (β = 1.32; 95% CI = 0.19, 2.46 mmHg; Fig. 2). Additionally, excluding people taking respiratory medications made the association between PNC and DBP levels significant (β = 1.03, 95% CI = 0.16, 1.89 mmHg). Nevertheless, controlling for the baseline level of DBP did not substantially affect the effect estimates (β = 0.15; 95% CI = − 0.72, 1.02 mmHg). In contrast, the PIR was associated with DBP (β = 1.01; 95% CI = 0.36, 1.66 mmHg per 6.2 billion particles inhaled/h). The positive associations between PIR and DBP were strongest among males, participants employed at baseline, and older participants (Additional file 1: Table S3). The associations between PIR and DBP were attenuated when excluding people with hypertension at baseline and, separately, excluding people who dropped out before their third study visit (β = 1.01, 95% CI = − 0.05, 2.07 mmHg; β = 0.64, 95% CI = − 0.18, 1.46 mmHg, respectively).
Similarly to the trends for SBP, PNC was not associated with PP overall (β = 0.70; 95% CI = − 0.27, 1.67 mmHg per 4600 particles/mL) but there was effect modification by sex (Fig. 2). Among females, PNC was positively associated with PP (β = 1.73; 95% CI = 0.52, 2.93 mmHg) while for males, PNC was inversely associated with PP (β = − 1.62; 95% CI = − 3.18, − 0.05 mmHg). The positive associations between PNC and PP were generally stronger among never smokers, participants who were not employed at baseline, the least physically active participants, older participants, non-diabetics, and participants with a family history of hypertension (Additional file 1: Table S2). When controlling for the baseline level of PP, PNC was associated with changes in PP (β = 1.69; 95% CI = 0.66, 2.72 mmHg). When using the PIR in place of PNC, however, exposure was not associated with PP (β = 0.14; 95% CI = − 0.62, 0.90 mmHg per 6.2 billion particles inhaled/h).
In contrast to the blood pressure measures, exposure to PNC was associated with 6.8% higher hsCRP overall (95% CI = − 0.3, 14.0% per 4600 particles/mL; Fig. 2). The associations were strongest among never smokers, people without a family history of CVD, and people taking statins (Additional file 1: Table S2). Including individuals with a hsCRP value ≥47, excluding individuals with a baseline hsCRP value of ≥3 mg/L, excluding individuals with a non-fasting baseline blood draw, or excluding individuals who dropped out before their third study visit all attenuated the association between PNC and levels of hsCRP (β = 6.3%, 95% CI = − 1.3, 14.0%; β = 2.8%, 95% CI = − 5.8, 11.4%; β = 5.6%, 95% CI = − 2.1, 13.4%; β = 5.0%, 95% CI = − 3.7, 13.6%, respectively). Additionally, in models controlling for the baseline level of hsCRP, PNC was not associated with changes in hsCRP (β = − 1.7; 95% CI = − 9.9, 6.5%). The PIR was also not associated with hsCRP overall (β = − 4.0; 95% CI = − 9.4, 1.3% per 6.2 billion particles inhaled/h), although there was an inverse association between PIR and hsCRP among diabetics, participants without a family history of CVD, and participants who were not taking statins (Additional file 1: Table S3).
Sensitivity analyses excluding people who died or excluding those who had a previous heart attack or stroke did not materially change the results (data not shown). Excluding hypertension medication as a covariate from the models did not materially change any of the results (data not shown).
Our study is one of only a few longitudinal studies to consider the health effects of long-term exposure to UFP, and it is the first to do so with UFP measured at high spatial and temporal resolution. Our primary finding is that annual average PNC exposure was positively associated with hsCRP but was not associated with the blood pressure measures overall. Sex and smoking status modified the associations of PNC and the outcomes with generally stronger positive associations among women and never smokers. In sub-populations where PNC affected blood pressure, the effect estimates were modest in terms of clinical significance. An IQR increase in PNC exposure was associated with increases in blood pressure approximately equivalent to the increase of blood pressure seen with an additional year of age .
The association between PNC and hsCRP was somewhat stronger than the associations with blood pressure despite the fact that 56% of participants had hsCRP values > 3 mg/L at baseline (indicative of high cardiovascular risk) . Each IQR increase in PNC was associated with a 7% higher mean hsCRP concentration (95% CI = − 0.3, 14.0%). This positive association is consistent with the only other longitudinal study examining the association between PNC and hsCRP . Similarly, a cross-sectional study among a different population in the Boston metropolitan area found positive associations between annual average PNC exposure and hsCRP . Furthermore, our finding that the associations of PNC with hsCRP and with blood pressure were stronger among never smokers compared to current or former smokers is consistent with the PM literature . This might be explained by the constant low-grade inflammation that smokers experience .
Although toxicological evidence suggests that UFP may exert cardiovascular effects, likely through mechanisms mediated by oxidative stress, inflammation, and endothelial dysfunction [55,56,57], the epidemiologic evidence is less clear regarding associations between UFP and cardiovascular risk factors. The vast majority of the previous epidemiological literature has only considered acute health effects or short-term exposure to UFP and the results have been inconsistent [58,59,60,61,62,63,64,65,66]. Our reported results do not account for potential acute changes in the outcome measures due to changes in short-term exposures. While it is possible that there are modest acute effects, a recent longitudinal study suggested that short-term effects did not substantially change the effect estimates for long-term exposures to PNC . Moreover, within our data, including short-term effects could affect the performance of the models due to the high correlation between the short- and long-term exposures. Furthermore, if short-term effects aggregate to the long-term effects, including the short-term effects may not be conceptually valid as the short-term effects would be on the causal pathway .
One of the innovations of our study was the use of an exploratory exposure metric, the PIR. We developed an algorithm to calculate PIR in a longitudinal study based on an established framework  and using validated age-, sex-, weight- and physical activity level-specific estimates for hourly respiratory volume . In a previous study, PIR was calculated by multiplying PNC by the amount of air inhaled per minute . That study found different exposure trends compared to studies that did not account for inhalation [69, 70]. Neglecting inhalation rate may thus introduce exposure misclassification.
Nevertheless, there are fundamental differences in what PNC and PIR attempt to quantify. For 3% of observations, participants had a high PIR (>75th percentile) but low PNC exposure (<25th percentile) and a high average inhalation rate (>75th percentile). For these people, the particle deposition fraction and clearance rate may be the most relevant factors in relation to health effects [7, 71]. This is particularly important for certain sub-populations who have higher inhalation rates, such as males and those with greater physical activity. In our study, we found some evidence that associations between the PIR and CVD risk factors were stronger among males than females. Additionally, as may be expected if there were a true association between UFP and DBP levels, the association between the PIR and DBP in males was stronger than the association between PNC and DBP in males.
We might expect people who inhale more particles per hour to have increased susceptibility to high ambient concentrations of air pollutants, especially if individuals are exercising in highly polluted areas . It is also possible that these people may be less vulnerable to health effects from air pollution if the higher PIR is due to greater levels of physical activity since exercise is associated with better cardiovascular health and the beneficial effect of exercise tends to outweigh the negative effect of exposure to air pollution [72, 73]. This was partially supported as the effect estimates for PNC and PIR were slightly stronger among the most sedentary participants in certain models. This trend was not entirely consistent, however, and in some cases, healthier individuals in our study appeared to be more susceptible. It is also possible that it was easier for us to observe associations among people who are not taking medications, such as statins, that could counteract the negative effects of PM exposure [74, 75].
Further work is needed to refine our PIR algorithm. In particular, it would be useful to refine the metric based on indoor measurements of UFP concentrations or based on an accounting of how much time participants spend in different indoor and outdoor micro-environments. Additionally, we may be over-accounting for weight or physical activity since including BMI as a covariate in the models for hsCRP changed the direction of the associations. This trend was not apparent in the BP models. Finally, it might be interesting to see if using age-, sex-, weight- and physical activity level estimates for hourly respiratory volume specific to this population (rather than from the EPA Exposure Factors Handbook) changes the results. While this may be expected if the health status of our participants differs substantially from that of the reference population, limiting our analyses to only those participants who were not taking respiratory medications and who did not have extremely high physical activity levels did not substantially change the effect estimates for the PIR. We believe that further developing the PIR metric is worthwhile since it addresses a critical step on the exposure pathway and may reduce exposure misclassification.
Our study had several limitations. One was the temporal mismatch between our exposure monitoring and participant visits. Participant visits occurred between 2004 and 2015 while we monitored UFP concentrations from December 2011 through November 2013. During this time, UFP emissions could have changed. Nevertheless, much of the temporal variability in PNC exposure is explained by meteorological conditions and we have historical data for these parameters [76, 77]. Additionally, modeled estimates compared reasonably well to PNC measurements at a stationary site for the years in which the model was back-extrapolated (see Additional file 1, Part 1). Annual trends were also accounted for in the health association models if the term for the year participants started the study was a significant predictor.
Other limitations of our exposure assessment include our assumptions that spatial variability was constant with time, that the model is valid up to 1000 m from the monitoring area, and that participants’ residential annual average PNC was representative of their overall personal exposure. Although we do not have time-activity data available in the BPRHS, we have shown in previous work that accounting for time-activity could reduce exposure misclassification . We also assumed that observed associations were due to long-term UFP exposure without accounting for potential interactions with, or independent effects of, other traffic-related pollutants or traffic-related noise. It was not possible to account for other traffic-related pollutants because we do not currently have exposure models for any other pollutants. Additionally, while our assumption that a one year averaging period represents a biologically relevant time-frame is in accordance with much of the related literature on long-term exposure to traffic-related air pollutants [2, 4], it is possible that we did not capture the critical averaging window for UFP. It is also possible that short-term effects would change the effect estimates for the long-term exposures in this study even though they did not in a previous longitudinal study .
Limitations also included the substantial attrition and resulting potential for selection bias. While only 49% of the 791 participants contributed data at all three time points, the baseline characteristics of participants who stayed in the study were similar to those who dropped out and most of the results were not substantially affected by excluding individuals who completed fewer than three study visits. Furthermore, our results may not be generalizable to healthier populations. All of our study participants were Puerto Rican, most had low socioeconomic status, and most had at least one chronic health condition at baseline. Other potential limitations of our analysis include the large number of comparisons and resulting possibility that some significant findings are due to chance, relatively low exposure contrast across the study population possibly limiting our ability to find true associations, the potential for misclassification of covariates such as diabetes, and the potential for residual confounding.
We found that both PNC and PIR were associated with biomarkers of CVD risk over six years, although the trends were not entirely consistent. The PIR is a novel exposure metric that accounts for differential inhalation rate. As our study is among the first to address these questions, future work is needed to validate these findings.
Body mass index
Boston Puerto Rican Health Study
Diastolic blood pressure
High sensitivity C-reactive protein
Particle inhalation rate
Particle number concentration
Systolic blood pressure
Ultrafine particulate matter
Brook RD, Rajagopalan S, Pope CA, Brook JR, Bhatnagar A, Diez-Roux AV, et al. Particulate matter air pollution and cardiovascular disease an update to the scientific statement from the American Heart Association. Circulation. 2010;121:2331–78.
Chuang K-J, Yan Y-H, Chiu S-Y, Cheng T-J. Long-term air pollution exposure and risk factors for cardiovascular diseases among the elderly in Taiwan. Occup Environ Med. 2010;oem.2009.052704.
Coogan PF, White LF, Jerrett M, Brook RD, Su JG, Seto E, et al. Air pollution and incidence of hypertension and diabetes mellitus in black women living in Los Angeles. Circulation. 2012;125:767–72.
Fuks K, Moebus S, Hertel S, Viehmann A, Nonnemacher M, Dragano N, et al. Long-term urban particulate air pollution, traffic noise, and arterial blood pressure. Environ Health Perspect. 2011;119:1706–11.
Hennig F, Fuks K, Moebus S, Weinmayr G, Memmesheimer M, Jakobs H, et al. Association between source-specific particulate matter air pollution and hs-CRP: local traffic and industrial emissions. Environ Health Perspect. 2014;122:703–10.
Zhang Z, Laden F, Forman JP, Hart JE. Long-term exposure to particulate matter and self-reported hypertension: a prospective analysis in the nurses’ health study. Environ Health Perspect. 2016;124:1414–20.
Lippmann M, Yeates DB, Albert RE. Deposition, retention, and clearance of inhaled particles. Br J Ind Med. 1980;37:337–62.
Jaques PA, Kim CS. Measurement of total lung deposition of inhaled ultrafine particles in healthy men and women. Inhal Toxicol. 2000;12:715–31.
Shukla A, Timblin C, BeruBe K, Gordon T, McKinney W, Driscoll K, et al. Inhaled particulate matter causes expression of nuclear factor (NF)- κ B–related genes and oxidant-dependent NF- κ B activation In Vitro. Am J Respir Cell Mol Biol. 2000;23:182–7.
Donaldson K, Brown D, Clouter A, Duffin R, MacNee W, Renwick L, et al. The pulmonary toxicology of ultrafine particles. J Aerosol Med Off J Int Soc Aerosols Med. 2002;15:213–20.
Li N, Sioutas C, Cho A, Schmitz D, Misra C, Sempf J, et al. Ultrafine particulate pollutants induce oxidative stress and mitochondrial damage. Environ Health Perspect. 2003;111:455–60.
Oberdörster G, Sharp Z, Atudorei V, Elder A, Gelein R, Kreyling W, et al. Translocation of inhaled ultrafine particles to the brain. Inhal Toxicol. 2004;16:437–45.
Kreyling WG, Semmler-Behnke M, Möller W. Ultrafine particle–lung interactions: does size matter? J Aerosol Med. 2006;19:74–83.
Stoeger T, Reinhard C, Takenaka S, Schroeppel A, Karg E, Ritter B, et al. Instillation of six different ultrafine carbon particles indicates a surface area threshold dose for acute lung inflammation in mice. Environ Health Perspect. 2006;114:328–33.
Rissler J, Swietlicki E, Bengtsson A, Boman C, Pagels J, Sandström T, et al. Experimental determination of deposition of diesel exhaust particles in the human respiratory tract. J Aerosol Sci. 2012;48:18–33.
Karner AA, Eisinger DS, Niemeier DA. Near-roadway air quality: synthesizing the findings from real-world data. Environ Sci Technol. 2010;44:5334–44.
Durant JL, Ash CA, Wood EC, Herndon SC, Jayne JT, Knighton WB, et al. Short-term variation in near-highway air pollutant gradients on a winter morning. Atmospheric Chem Phys Print. 2010;10:5599–626.
Grundström M, Hak C, Chen D, Hallquist M, Pleijel H. Variation and co-variation of PM10, particle number concentration, NOx and NO2 in the urban air – relationships with wind speed, vertical temperature gradient and weather type. Atmos Environ. 2015;120:317–27.
Kumar P, Fennell P, Britter R. Effect of wind direction and speed on the dispersion of nucleation and accumulation mode particles in an urban street canyon. Sci Total Environ. 2008;402:82–94.
Padró-Martínez LT, Patton AP, Trull JB, Zamore W, Brugge D, Durant JL. Mobile monitoring of particle number concentration and other traffic-related air pollutants in a near-highway neighborhood over the course of a year. Atmos Env. 2012;61:253–64.
Pattinson W, Longley I, Kingham S. Using mobile monitoring to visualise diurnal variation of traffic pollutants across two near-highway neighbourhoods. Atmos Environ. 2014;94:782–92.
Zhu Y, Kuhn T, Mayo P, Hinds WC. Comparison of daytime and nighttime concentration profiles and size distributions of ultrafine particles near a major highway. Environ Sci Technol. 2006;40:2531–6.
Lane KJ, Levy JI, Scammell MK, Peters JL, Patton AP, Reisner E, et al. Association of modeled long-term personal exposure to ultrafine particles with inflammatory and coagulation biomarkers. Environ Int. 2016;92–93:173–82.
Ostro B, Hu J, Goldberg D, Reynolds P, Hertz A, Bernstein L, et al. Associations of Mortality with Long-Term Exposures to Fine and Ultrafine Particles, Species and Sources: Results from the California Teachers Study Cohort. Environ Health Perspect [Internet]. 2015;123:549-556. Available from: http://ehp.niehs.nih.gov/1408565. [cited 2015 Mar 18]. doi:https://doi.org/10.1289/ehp.1408565.
Aguilera I, Dratva J, Caviezel S, Burdet L, de Groot E, Ducret-Stich RE, et al. Particulate matter and subclinical atherosclerosis: associations between different particle sizes and sources with carotid intima-media thickness in the SAPALDIA study. Environ Health Perspect. 2016;
Viehmann A, Hertel S, Fuks K, Eisele L, Moebus S, Möhlenkamp S, et al. Long-term residential exposure to urban air pollution, and repeated measures of systemic blood markers of inflammation and coagulation. Occup Environ Med. 2015;72:656–63.
Weichenthal S, Bai L, Hatzopoulou M, Van Ryswyk K, Kwong JC, Jerrett M, et al. Long-term exposure to ambient ultrafine particles and respiratory disease incidence in in Toronto, Canada: a cohort study. Environ Health. 2017;16:64.
Tucker KL, Mattei J, Noel SE, Collado BM, Mendez J, Nelson J, et al. The Boston Puerto Rican health study, a longitudinal cohort study on health disparities in Puerto Rican adults: challenges and opportunities. BMC Public Health. 2010;10:107.
Rioux CL, Tucker KL, Mwamburi M, Gute DM, Cohen SA, Brugge D. Residential traffic exposure, pulse pressure, and C-reactive protein: consistency and contrast among exposure characterization methods. Environ Health Perspect. 2010;118:803–11.
US Census Bureau. Poverty Thresholds [Internet]. 2017 [cited 2017 Jul 10]. Available from: https://www.census.gov/data/tables/time-series/demo/income-poverty/historical-poverty-thresholds.html.
Paffenbarger RS, Wing AL, Hyde RT. Physical activity as an index of heart attack risk in college alumni. Am J Epidemiol. 1978;108:161–75.
Tucker KL, Bermudez OI, Castaneda C. Type 2 diabetes is prevalent and poorly controlled among Hispanic elders of Caribbean origin. Am J Public Health. 2000;90:1288–93.
Tropp LR, Erkut S, Coll CG, Alarcón O, Vázquez García HA. Psychological acculturation: development of a new measure for Puerto Ricans on the U.S. mainland. Educ Psychol Meas. 1999;59:351–67.
Cohen S, Kamarck T, Mermelstein R. A global measure of perceived stress. J Health Soc Behav. 1983;24:385.
Miller TQ, Markides KS, Black SA. The factor structure of the CES-D in two surveys of elderly Mexican Americans. J Gerontol B Psychol Sci Soc Sci. 1997;52:S259–69.
Mościcki EK, Locke BZ, Rae DS, Boyd JH. Depressive symptoms among Mexican Americans: the Hispanic health and nutrition examination survey. Am J Epidemiol. 1989;130:348–60.
Tucker KL, Falcon LM, Bianchi LA, Cacho E, Bermudez OI. Self-reported prevalence and health correlates of functional limitation among Massachusetts elderly Puerto Ricans, Dominicans, and non-Hispanic white neighborhood comparison group. J Gerontol A Biol Sci Med Sci. 2000;55:M90–7.
Clyne B, Olshaker JS. The C-reactive protein. J Emerg Med. 1999;17:1019–25.
Rioux CL, Gute DM, Brugge D, Peterson S, Parmenter B. Characterizing urban traffic exposures using transportation planning tools: an illustrated methodology for health researchers. J urban health. Bull N Y Acad Med. 2010;87:167–88.
Patton AP, Collins C, Naumova EN, Zamore W, Brugge D, Durant JL. An hourly regression model for ultrafine particles in a near-highway urban area. Environ Sci Technol. 2014;48:3272–80.
Simon MC, Hudda N, Naumova EN, Levy JI, Brugge D, Durant JL. Comparisons of traffic-related ultrafine particle number concentrations measured in two urban areas by central, residential, and mobile monitoring. Atmos Environ. 2017;169:113–27.
National Climatic Data Center. Land-Based Station Data [Internet]. 2015 [cited 2015 Mar 16]. Available from: http://www.ncdc.noaa.gov/data-access/land-based-station-data.
MassGIS. MassGIS Datalayers [Internet]. Adm. Finance. 2012 [cited 2015 Aug 13]. Available from: http://www.mass.gov/anf/research-and-tech/it-serv-and-support/application-serv/office-of-geographic-information-massgis/datalayers/layerlist.html
Patton AP, Zamore W, Naumova EN, Levy JI, Brugge D, Durant JL. Transferability and generalizability of regression models of ultrafine particles in urban neighborhoods in the Boston area. Environ Sci Technol. 2015;49:6051–60.
Montagne DR, Hoek G, Klompmaker JO, Wang M, Meliefste K, Brunekreef B. Land use regression models for ultrafine particles and black carbon based on short-term monitoring predict past spatial variation. Environ Sci Technol. 2015;49:8712–20.
Abernethy RC, Allen RW, McKendry IG, Brauer M. A land use regression model for ultrafine particles in Vancouver, Canada. Environ Sci Technol. 2013;47:5217–25.
Levy I, Levin N, Yuval, Schwartz JD, Kark JD. Back-extrapolating a land use regression model for estimating past exposures to traffic-related air pollution. Environ Sci Technol. 2015;49:3603–10.
US EPA. Metabolically Derived Human Ventilation Rates: A Revised Approach Based Upon Oxygen Consumption Rates (Final Report) 2009 [Internet]. 2009 [cited 2015 Jul 27]. Available from: http://cfpub.epa.gov/ncea/cfm/recordisplay.cfm?deid=202543.
StataCorp. Stata Statistical Software: Release 14 [Internet]. College Station, TX: StataCorp LP; 2015 [cited 2016 Oct 31]. Available from: http://www.stata.com/support/faqs/resources/citing-software-documentation-faqs/
Poverty Research [Internet]. ASPE. 2015 [cited 2018 Jan 30]. Available from: https://aspe.hhs.gov/poverty-research
Franklin SS, Gustin W, Wong ND, Larson MG, Weber MA, Kannel WB, et al. Hemodynamic patterns of age-related changes in blood pressure. Circulation. 1997;96:308–15.
Pearson TA, Mensah GA, Alexander RW, Anderson JL, Cannon RO, Criqui M, et al. Markers of inflammation and cardiovascular disease application to clinical and public health practice: a statement for healthcare professionals from the Centers for Disease Control and Prevention and the American Heart Association. Circulation. 2003;107:499–511.
Pope CA 3rd, Burnett RT, Thun MJ, Calle EE, Krewski D, Ito K, et al. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution. J Am Med Assoc. 2002;287:1132–41.
Yasue H, Hirai N, Mizuno Y, Harada E, Itoh T, Yoshimura M, et al. Low-grade inflammation, Thrombogenicity, and Atherogenic lipid profile in cigarette smokers. Circ J. 2006;70:8–13.
Upadhyay S, Stoeger T, George L, Schladweiler MC, Kodavanti U, Ganguly K, et al. Ultrafine carbon particle mediated cardiovascular impairment of aged spontaneously hypertensive rats. Part Fibre Toxicol. 2014;11:36.
Araujo JA. Particulate air pollution, systemic oxidative stress, inflammation, and atherosclerosis. Air Qual Atmosphere Health. 2010;4:79–93.
Schmid O, Möller W, Semmler-Behnke M, Ferron GA, Karg E, Lipka J, et al. Dosimetry and toxicology of inhaled ultrafine particles. Biomarkers. 2009;14:67–73.
Rückerl R, Ibald-Mulli A, Koenig W, Schneider A, Woelke G, Cyrys J, et al. Air pollution and markers of inflammation and coagulation in patients with coronary heart disease. Am J Respir Crit Care Med. 2006;173:432–41.
Zeka A, Sullivan JR, Vokonas PS, Sparrow D, Schwartz J. Inflammatory markers and particulate air pollution: characterizing the pathway to disease. Int J Epidemiol. 2006;35:1347–54.
Rückerl R, Greven S, Ljungman P, Aalto P, Antoniades C, Bellander T, et al. Air pollution and inflammation (interleukin-6, C-reactive protein, fibrinogen) in myocardial infarction survivors. Environ Health Perspect. 2007;115:1072–80.
Yue W, Schneider A, Stölzel M, Rückerl R, Cyrys J, Pan X, et al. Ambient source-specific particles are associated with prolonged repolarization and increased levels of inflammation in male coronary artery disease patients. Mutat Res. 2007;621:50–60.
Delfino RJ, Staimer N, Tjoa T, Polidori A, Arhami M, Gillen DL, et al. Circulating biomarkers of inflammation, antioxidant activity, and platelet activation are associated with primary combustion aerosols in subjects with coronary artery disease. Environ Health Perspect. 2008;116:898–906.
Delfino RJ, Tjoa T, Gillen DL, Staimer N, Polidori A, Arhami M, et al. Traffic-related Air Pollution and Blood Pressure in Elderly Subjects With Coronary Artery Disease. Epidemiol Camb Mass [Internet]. 2010;21(3). Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3872093/. [cited 2016 Dec 31];21. doi:https://doi.org/10.1097/EDE.0b013e3181d5e19b.
Hertel S, Viehmann A, Moebus S, Mann K, Bröcker-Preuss M, Möhlenkamp S, et al. Influence of short-term exposure to ultrafine and fine particles on systemic inflammation. Eur J Epidemiol. 2010;25:581–92.
Rich DQ, Zareba W, Beckett W, Hopke PK, Oakes D, Frampton MW, et al. Are ambient ultrafine, accumulation mode, and fine particles associated with adverse cardiac responses in patients undergoing cardiac rehabilitation? Environ Health Perspect. 2012;120:1162–9.
Chung M, Wang DD, Rizzo AM, Gachette D, Delnord M, Parambi R, et al. Association of PNC, BC, and PM2.5 measured at a central monitoring site with blood pressure in a predominantly near highway population. Int J Environ Res Public Health. 2015;12:2765–80.
Knol AB, de Hartog JJ, Boogaard H, Slottje P, van der Sluijs JP, Lebret E, et al. Expert elicitation on ultrafine particles: likelihood of health effects and causal pathways. Part Fibre Toxicol. 2009;6:19.
Int Panis L, de Geus B, Vandenbulcke G, Willems H, Degraeuwe B, Bleux N, et al. Exposure to particulate matter in traffic: a comparison of cyclists and car passengers. Atmos Environ. 2010;44:2263–70.
Kaur S, Nieuwenhuijsen M, Colvile R. Personal exposure of street canyon intersection users to PM2.5, ultrafine particle counts and carbon monoxide in Central London, UK. Atmos Environ. 2005;39:3629–41.
Boogaard H, Borgman F, Kamminga J, Hoek G. Exposure to ultrafine and fine particles and noise during cycling and driving in 11 Dutch cities. Atmos Environ. 2009;43:4234–42.
Leikauf LG. Toxic responses of the respiratory system. In: Klaassen C, Watkins JB, editors. Casarett Doulls Essent Toxicol. New York: McGraw-Hill; 2010.
Sharman JE, Cockcroft JR, Coombes JS. Cardiovascular implications of exposure to traffic air pollution during exercise. QJM Int J Med. 2004;97:637–43.
Thompson PD, Buchner D, Piña IL, Balady GJ, Williams MA, Marcus BH, et al. Exercise and physical activity in the prevention and treatment of atherosclerotic cardiovascular disease. Circulation. 2003;107:3109–16.
Schwartz J, Park SK, O’Neill MS, Vokonas PS, Sparrow D, Weiss S, et al. Glutathione-S-transferase M1, obesity, statins, and autonomic effects of particles. Am J Respir Crit Care Med. 2005;172:1529–33.
Miyata R, Hiraiwa K, Cheng JC, Bai N, Vincent R, Francis GA, et al. Statins attenuate the development of atherosclerosis and endothelial dysfunction induced by exposure to urban particulate matter (PM10). Toxicol Appl Pharmacol. 2013;272:1–11.
Noble CA, Mukerjee S, Gonzales M, Rodes CE, Lawless PA, Natarajan S, et al. Continuous measurement of fine and ultrafine particulate matter, criteria pollutants and meteorological conditions in urban El Paso, Texas. Atmos Environ. 2003;37:827–40.
Kozawa KH, Winer AM, Fruin SA. Ultrafine particle size distributions near freeways: effects of differing wind directions on exposure. Atmos Environ. 2012;63:250–60.
Lane KJ, Kangsen Scammell M, Levy JI, Fuller CH, Parambi R, Zamore W, et al. Positional error and time-activity patterns in near-highway proximity studies: an exposure misclassification analysis. Environ Health. 2013;12:75.
We would like to acknowledge the Boston Puerto Rican Health Study team including Lori Lyn, Shilpa Saklani, and Esther Carver. We are also grateful for the assistance of Kevin Lane, Wig Zamore, Cara Goodman, and Allison Patton. Particle number concentration data used for model validation were collected by the Harvard T.H. Chan School of Public Health with support from USEPA (RD 83479801) and NIEHS (PO1ES009825).
Funding for this work was provided by Tufts University Department of Civil and Environmental Engineering, Tufts Institute of the Environment, the National Heart, Lung, and Blood Institute (P01 AG023394 and P50 HL105185), the National Institute of Environmental Health Sciences (ES015462), the National Science Foundation (0966093), and a Philanthropic Educational Organization Esther Garrett Edgerton Endowed Scholar Award. No funder had any role in the design of the study, collection of the data, analysis of the data, interpretation of the results, or writing the manuscript.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from Doug Brugge (firstname.lastname@example.org) on reasonable request.
Ethics approval and consent to participate
This study was approved by the Institutional Review Boards at Tufts Medical Center, Northeastern University (protocol #10–04-12), and the University of Massachusetts Lowell. All participants provided written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Part 1. Exposure assessment and Table S1. PNC model used for exposure assignment; Part 2. Inhalation rate adjustment; Part 3. Table S2. Longitudinal associations with an IQR increase in PNC (4600 particles/mL) and Table S3. Longitudinal associations with an IQR increase in PIR (6.2 billion particles inhaled/h). (DOCX 62 kb)