 Research
 Open Access
 Open Peer Review
 Published:
Structural equation modeling of parasympathetic and sympathetic response to traffic air pollution in a repeated measures study
Environmental Health volume 12, Article number: 81 (2013)
Abstract
Background
Trafficrelated air pollution has been associated to a range of adverse health impacts, including decreased heart rate variability (HRV). The association between trafficrelated pollution and HRV, however, has varied by trafficrelated or HRV marker as well as by study, suggesting the need for a more comprehensive and integrative approach to examining air pollutionmediated biological impacts on these outcomes. In a Bayesian framework, we examined the effect of traffic pollution on HRV using structural equation models (SEMs) and looked at effect modification by participant characteristics.
Methods
We studied measurements of 5 HRV markers [high frequency (HF), low frequency (LF), 5min standard deviation of normaltonormal intervals (SDNN), square root of the mean squared differences of successive normaltonormal intervals (rMSSD), and LF/HF ratio (LF/HF)] for 700 elderly men from the Normative Aging Study. Using SEMs, we fit a latent variable for traffic pollution that is reflected by levels of carbon monoxide, nitrogen monoxide, nitrogen dioxide, and black carbon (BC) to estimate its effect on latent variable for parasympathetic tone that included HF, SDNN and rMSSD, and the sympathetic tone marker, LF/HF. Exposure periods were assessed using 4, 24, 48, 72hour moving average previsit. We compared our main effect findings using SEMs with those obtained using linear mixed models.
Results
Traffic pollution was not associated with mean parasympathetic tone and LF/HF for all examined moving averages. In Bayesian linear mixed models, however, BC was related to increased LF/HF, an inter quartile range (IQR) increase in BC was associated with a 6.5% (95% posterior interval (PI): 0.7%, 14.2%) increase in mean LF/HF 24hours later. The strongest association observed was for the 4hour moving average (10.1%; 95% PI: 3.0%, 17.6%). The effect of traffic on parasympathetic tone was stronger among diabetic as compared to nondiabetic participants. Specifically, an IQR increase in traffic pollution in the 48hr prior to the clinic visit was associated with a 44.3% (95% PI: 67.7%, 4.2%) lower mean parasympathetic tone among diabetics, and a 7.7% (95% PI: 18.0%, 41.4%) higher mean parasympathetic tone among nondiabetics.
Conclusions
BC was associated with adverse changes LF/HF in the elderly. Traffic pollution may decrease parasympathetic tone among diabetic elderly.
Background
Studies in air pollution epidemiology have correlated increased levels of trafficrelated pollution to increased cardiovascular morbidity and mortality and have also recognized trafficrelated pollution as an important risk factor [1]. While the physiological mechanisms of this association have not been fully explained, existing scientific evidence suggests that altered cardiac autonomic control plays a significant role [2, 3]. For instance, Schwartz et al. in 2005 found associations between trafficrelated pollution exposures and disturbances of autonomic control of the heart as measured through heart rate variability (HRV), a measure of naturally occurring beattobeat interval in heart rate [4]. This association between trafficrelated pollution and HRV, however, has varied by trafficrelated pollutants or HRV markers as well as by study, suggesting the need for a more comprehensive and integrative approach to examining air pollutionmediated biological impacts on these outcomes and in addition to consider multiple trafficrelated pollutants and multiple health markers simultaneously.
We used structural equation models (SEMs) to address this issue. SEMs, which represent a family of statistical techniques, allow one to estimate association among multiple latent variables, for example, the relation between traffic pollution and parasympathetic tone. For our analyses, we used latent variables to conceptualize cardiac parasympathetic tone and traffic pollution. These latent variables, traffic and parasympathetic tone, are correlated with measured trafficrelated pollutants and markers of HRV, respectively, however they are not observed directly. The latent variables correspond to what is “common” among the parameters measured and do not preclude in part different biological processes. Several studies have used SEMs that assessed sourcespecific health effects of air pollution [5, 6] and in health effects of methyl mercury [7] and lead on neurodevelopment [8]. Notably, in models with 3 or more surrogates of the latent exposure, SEMs reduce the attenuation due to measurement error which is an ongoing concern in air pollution studies [9]. This is a typical phenomenon in measurement error corrections, it tends to deattenuate point estimates coupled with wider confidence intervals due to reflecting additional uncertainty as a result of measurement error correction. The hope is that the total mean squared error (bias squared + variance) is lower than uncorrected estimate.
In this paper, we develop in a Bayesian framework SEMs to examine separately the impact of shortterm changes in traffic pollution on parasympathetic tone, and a sympathetic tone marker, the low frequency HRV (LF) to high frequency HRV (HF) ratio (LF/HF), among participants in the Normative Aging Study (NAS). In addition, we look at effect modification by participant characteristics and also compare and contrast the main results of the SEMs to more standard linear mixed models in the Bayesian linear mixed models (BLMMs) framework.
Methods
Study population
The NAS is an ongoing longitudinal study of aging launched in 1963 by the Veterans Administration [10]. 2,280 communitydwelling healthy men living in the Greater Boston area were enrolled between 1963 and 1968. Every 3–5 years, after an overnight fast and abstention from smoking, participants visited the clinic for an extensive physical examination, laboratory tests, blood collection, and a selfadministered questionnaire on medical history, food intake, medication usage, alcohol consumption, smoking history, and other factors that could affect health.
Outcome data
The cardiac rhythm was measured for approximately 7 minutes. From November 2000 to December 2009, square root of the mean squared differences of successive normaltonormal intervals (rMSSD), 5min standard deviation (SD) of normaltonormal intervals (SDNN), HF, LF, and LF/HF measurements were obtained during each participant’s regularly scheduled visit [11, 12]. To date, HRV measurements have been made for 700 participants, in some cases on multiple occasions (214 participants with one, 216 with two, 259 with three, and 11 with four measurements). A total of 1467 valid measurements were included in our analyses. All participants gave written consent and this study had Institutional Review Board approval.
Exposure and weather data
We measured concentrations of ambient BC continuously at the Harvard University Countway Library stationary ambient monitoring site that is located < 1 km from the clinical laboratory where subjects were examined. Hourly ambient carbon monoxide (CO), nitrogen monoxide (NO), and nitrogen dioxide (NO_{2}) concentrations were obtained from the Massachusetts Department of Environmental Protection (MDEP) for monitoring stations within 20 km of the clinical laboratory, including those in Boston (Bremen, Kenmore Square, NorthEnd, and Roxbury), Waltham and Lynn MA. The pollutant exposure was averaged across all the MDEP monitoring stations; data from the MDEP monitoring stations were used to calculate the mean hourly concentrations of each pollutant. We obtained data on dew point temperature and ambient temperature from the first order National Weather Service station at Boston Logan airport. We calculated apparent temperature (ATemp), a human discomfort index [13], as: ATemp = −2.653 + (0.994 × AT) + (0.0153 × DPT^{2}), where DPT is the dew point temperature in Celsius and AT is the ambient temperature in Celsius.
Statistical analysis
We examined the association between latent traffic pollution and latent parasympathetic tone and between latent traffic pollution and sympathetic tone marker, LF/HF, using SEMs in a Bayesian framework that account for repeated measures. SEMs consist of two key components; a structural model part which shows the association dependencies between two latent constructs, between a latent construct and a measured variable or between two measured variables, and a measurement model part which shows the relation between latent constructs and their indicators/markers [9, 14]. HF, SDNN, rMSSD, and LF/HF were log transformed to satisfy the standard SEM assumption that all continuous variables are approximately normally distributed. Descriptive statistics and correlations between HRV markers, trafficrelated pollutants and covariates were calculated.
Structural equation model with repeated measures in Bayesian framework
We examined associations between trafficrelated pollution and markers of parasympathetic tone and sympathetic tone using SEMs, parametric Directed Acyclic Graphs (DAGs) that contain specified paths connecting observed and latent variables [9, 14]. SEMs in Bayesian framework consist of (1) distributional assumptions on the latent variables, (2) prior distributions for the parameters, and (3) probability models for the observed variables given latent variables and model parameters. SEMs also specify linear models relating the unknown means of the latent traffic pollution, sympathetic tone marker LF/HF and latent parasympathetic tone variables to fixed covariates (e.g. confounding factors), other latent variables, and a subject level random intercept to accommodate the repeated measures design of the study. For the effect of traffic pollution on both LF/HF and the latent parasympathetic tone, exposure periods were assessed using 4hr, 24hr, 48hr and 72hr moving average previsit, given results from previous studies that show effects at similar exposure windows [4, 11].
For the measurement part of the SEM, we specified measurement models describing the relationships between observed measures and latent variables. The relationship between the trafficrelated observed measures (BC, CO, NO, and NO_{2}) and the latent exposure variable Traffic was first modeled (Figure 1).
For moving average concentration previsit measurement of participants i = 1, …I, of trafficrelated pollutants j = 1, …J, on participant’s scheduled time visits t = 1, …T,
where X _{ ijt } is the normally distributed trafficrelated exposure variable for moving average concentration previsit measurement of participant i of trafficrelated pollutant j on participant’s scheduled time visit t, with unknown mean μ _{ ijt } ^{X} and precision (1/variance) τ _{ j } ^{X}. The mean exposure of trafficrelated pollutant j, μ _{ ijt } ^{X}, was related to the latent exposure variable Traffic _{ it } through a linear model,
where Traffic _{ it } is the normally distributed latent exposure variable Traffic for participant i for scheduled time visit t, with unknown mean μ _{ it } ^{Traffic} and precision (1/variance) τ ^{Traffic},
The mean Traffic exposure, μ _{ it } ^{Traffic}, was related to confounding factors, through a linear model,
where θ (t) is the calendar day of the year the participant visited the clinic, the sine and cosine terms capture seasonal trends in traffic pollution, ATemp is apparent temperature averaged across the same exposure window, and residual (ATemp) is the residual of a model regressing ATemp against sine and cosine terms of the calendar day of the year in the study. This approach was taken because temperature data was highly seasonal, and this avoids substantial collinearity between the seasonal terms and temperature, which would otherwise be present in the model. Furthermore, Eq. (4) looks like it doesn’t vary with participant i but it does vary because the scheduled time visits t are unique to that participant.
The relationship between the latent outcome variable, and its observed markers was also modeled. The HRV markers HF, SDNN and rMSSD were used as the markers for the latent Parasympathetic tone variable. For HRV measurement of participants i = 1, …I, of HRV parasympathetic tone markers k = 1, …K, on participant’s scheduled time visits t = 1, …T,
where Y _{ ikt } represents the normally distributed markers of parasympathetic tone for HRV measurement of participant i of parasympathetic tone marker k on participant’s scheduled time visit t, with unknown mean μ _{ ikt } ^{Y} and precision (1/variance) τ _{ k } ^{Y}. For parasympathetic tone marker k, the mean outcome, μ _{ ikt } ^{Y}, was related to the latent outcome variable Parasympathetic tone through a linear model,
where Parasympathetic tone _{ it } is the normally distributed latent outcome variable Parasympathetic tone for scheduled time visit t for participant i, with unknown mean μ _{ it } ^{Parasympathetic tone} and precision (1/variance) τ ^{Parasympathetic tone},
Furthermore, for LF/HF measurement of participants i = 1, …I, on participant’s scheduled time visits t = 1, …T,
where LF/HF_{ it } is the normally distributed sympathetic tone marker outcome for participant i for scheduled time visit t, with unknown mean μ _{ it } ^{LF/HF} and precision (1/variance) τ ^{LF/HF}.
The mean Parasympathetic tone outcome, μ _{ it } ^{Parasympathetic tone} was related to the latent exposure variable Traffic _{ it }, confounding factors, predictors of the outcome variable, and a subject level random intercept, a _{ i }, through a linear model,
In addition, the mean LF/HF outcome, μ _{ it } ^{LF/HF}, was also related to the latent exposure variable Traffic _{ it }, predictors of the outcome variable, confounding factors, and a subject level random intercept, b _{ i }, through a linear model,
For Eqs. (9) and (10), the fixed covariates (X _{ 1 } , …, X _{ m }) included age, fasting blood glucose, body mass index (BMI), smoking status (ever/never), room temperature, sitting mean arterial blood pressure, ≥ 2 servings of alcohol per day (yes/no), and cardiac medication usage (beta blocker, angiotensin converting enzyme inhibitor, and calcium channel blocker). These covariates were chosen a priori as potentially important predictors of both Parasympathetic tone and LF/HF. To account for seasonal variation in Parasympathetic tone and LF/HF, a function of calendar date {sine [2πθ(t)/365.24] + cosine [2πθ(t)/365.24], where θ (t) is the calendar day of the year the participant visited the clinic} was used. To account for temperature variation in both Parasympathetic tone and LF/HF, we regressed ATemp on sine [2πθ(t)/365.24] + cosine [2πθ(t)/365.24] and used the linear [Residual (ATemp)_{t}] and quadratic [Residual^{2} (ATemp)_{t}] residuals of the equation as part of the fixed effects of the model. We also included a subject level random intercept, a _{ i } for Eq. (9) and b _{ i } for Eq. (10), to represent subjectspecific permanent effects for each participant to accommodate the repeated measures design of the study. Moreover, Eqs. (9) and (10) implicitly depend on participant i because the scheduled time visits t depend on participant i. To evaluate obesity as an effect modifier, subjects were classified into two groups according to BMI (obese: BMI ≥ 30 vs. nonobese: BMI < 30). We also assessed effect modification by a history of diabetes (diabetic vs. nondiabetic) and by the four possible combinations of obesity and history of diabetes (diabetic and obese, diabetic and nonobese, nondiabetic and obese, and nondiabetic and nonobese). We included latent interaction term constructs between the dichotomized modifier variable and traffic pollution in Equations (9) and (10).
For model identifiability, the location of all the latent variables was centered and set to 0 (i.e. δ _{ 0 } , φ _{ 0 } and γ _{ 0 } = 0). A scale must be specified for each latent variable because latent variables were constructed from multiple measured parameters, which have different scales of variation. We chose HF as the reference scale for Parasympathetic tone and BC as the reference scale for Traffic pollution. In order to set the reference scale for Parasympathetic tone and Traffic and for identifiability and interpretation of these latent variables, α_{1} ^{HF} and λ_{1} ^{BC}, were also constrained to be equal to 1 [6]. For the DAG diagram of the SEMs, see Figure 1.
Linear mixed models in Bayesian framework
In order to compare the results from the SEM in Bayesian framework with repeated measures to those obtained from a standard univariate analysis, associations between trafficrelated air pollutant BC and mean markers of parasympathetic tone (HF) and sympathetic tone (LF/HF) were estimated using linear mixed models with random subjectspecific intercepts in Bayesian linear mixed models (BLMMs) framework. Additionally, we compared the results from the above BLMMs to Frequentist framework estimates (Frequentist linear mixed models, LMMs) obtained from LME in R (http://www.rproject.org). The same fixed and random covariates used in SEMs were included in all the LMMs and BLMMs. To allow more direct comparison for the SEMs, BLMMs, and LMMs, effect size estimates were reported for an interquartile range (IQR) increase in BC (reference pollutant of latent variable Traffic pollution).
We specified a model for the measured health effect outcomes, HF and LF/HF. The means of both HF and LF/HF for participant i on participant’s scheduled time visit t, E[HF _{ it } ] and E[LF/HF _{ it } ], respectively, was related to BC, confounding factors and predictors of the health effect outcomes, and a subject level random intercept, c _{ i } or d _{ i }, through a linear model,
In terms of Bayesian inference, results were quantified and characterized in terms of the strength of the following hypothesized relationships: (1) inverse relationships between BC and HRV marker HF, and between Traffic pollution and Parasympathetic tone, and (2) positive relationships between BC and the HRV marker LF/HF, and between Traffic pollution and sympathetic tone marker, LF/HF. In addition, results were presented as probability statements, providing evidence that coefficients describing Traffic pollution (γ _{ 1 } and δ _{ 1 }) or BC (β _{ 1 } and ρ _{ 1 }) are either negative or positive depending on the hypothesized relationship. See Additional file 1 for details for implementation particulars and description of prior distributions of parameters, and see Additional file 2 for the WinBUGS code.
Results
700 NAS eligible study participants with valid HRV measurements available were included for analysis. Subjects were male with a mean age of 75.0 years (SD = 6.7 years), were mostly ever cigarette smokers (69.5%), and were mostly overweight with a mean BMI of 28.1 kg/m^{2} (SD = 4.2 kg/m^{2}). Table 1 shows the descriptive statistics for ambient air pollutants and HRV measures. Traffic pollutant concentrations, averaged over 24 hours before blood collection, were strongly correlated with CO and NO_{2} having the highest correlation (Spearman correlation coefficient, ρ = 0.63). The correlation between HRV markers was highest for HF and rMSSD (ρ = 0.93). The LF/HF marker is negatively correlated with the rest of the HRV markers.
Traffic pollution and parasympathetic tone measurement models
Table 2 tabulates the factor loadings (λ _{ 1j } coefficients of Equation 2) of the Traffic Pollution measurement model and the corresponding variance of the Traffic latent (1/τ ^{Traffic} of Equation 3) and its marker variables (1/τ _{ j } ^{X} of Equation 1). Moreover, for the measurement model for the relation of the Parasympathetic tone latent variable to its various marker variables, the factor loadings (α _{ 1k } coefficients of Equation 6) and the variance of the Parasympathetic tone (1/τ ^{Parasympathetic tone} of Equation 7) and its marker variables (1/τ _{ k } ^{Y} of Equation 5) were also summarized in Table 2. Across all the different moving averages, the results of the factor loading estimates (α _{ 1k }) for the relation of Parasympathetic Tone to its markers were very similar. In using BC as the reference marker, for the 4hr moving average exposure, the highest reliability measure of Traffic was CO [coefficient as factor loading of 0.65; 95% Posterior Interval (PI): 0.59 to 0.72], while the lowest reliability measure of Traffic was NO_{2} (0.01; 95% PI: 0.01 to 0.01). Furthermore, for the 24hr moving average exposure, CO predominantly represented Traffic (0.95; 95% PI: 0.86 to 1.05) while NO_{2} was still the lowest reliability measure of Traffic (0.02; 95% PI: 0.02 to 0.02) (Figure 2A). In addition, in using HF as the reference marker for the 24hr moving average exposure, Parasympathetic tone was strongly represented by rMSSD with a coefficient as factor loading of 0.56 (95% PI: 0.55 to 0.57) (Figure 2B). In addition, for the 4hr moving average exposure, the variance of latent Parasympathetic Tone variable was higher (1.720; 95% PI: 1.535 to 1.921) compared to the variance of the latent Traffic variable (0.137; 95% PI: 0.112 to 0.165). For the HRV markers, for a 4hr moving average exposure, the variance of HF was the highest (0.309; 95% PI: 0.273 to 0.346) and the variance of rMSSD was the lowest of the 3 markers (0.052; 95% PI: 0.043 to 0.062). Of the 4 trafficrelated pollutants, the variance of the 4hr moving average BC was the highest (0.394; 95% PI: 0.364 to 0.425), while NO_{2} was the lowest in variance (4.3E5; 95% PI: 4.0E5 to 4.7E5) for all the pollutants.
Effect of traffic pollution on parasympathetic tone and sympathetic tone marker LF/HF
Table 3 shows the estimated percent change in parasympathetic tone per IQR increase in BC, the reference pollutant of the latent traffic pollution, for various hourlymoving averages. For the 48hr previsit moving average, a change of −3.0% (95% PI: 24.7% to 25.4%) in parasympathetic tone (reference marker: HF) was associated with an IQR increase in BC (reference pollutant of traffic pollution exposure). Results showed that the 4hr posterior distribution gave the smallest variance compared to other hourly moving averages. Table 3 also presents the estimated percent change in LF/HF, a marker of sympathetic tone, per IQR increase in BC (reference pollutant of latent traffic pollution) for various hourlymoving averages. We found a change of 6.4% (95%PI: 9.5% to 24.3%) in LF/HF associated with an IQR increase in BC (reference pollutant of traffic pollution) 48hr previsit. The distributions for SEM also showed that the 4hr posterior distribution gave the smallest variance compared to other hourly moving averages and as the averaging period increases from 4hr to 72hr, the variability or spread of the distribution also increased, as evidenced by the wider posterior probability distributions. Furthermore, in both the SEM and BLMM approaches, the posterior probabilities characterized and quantified strong positive relationships in all the exposure windows (4hr, 24hr, 48hr, and 72hr) between Traffic pollution and sympathetic tone marker, LF/HF, and between BC and HRV marker LF/HF (posterior probabilities ≥ 0.78).
Effect modification
We assessed whether being diabetic, obese, or being both obese and diabetic, modified the effect of traffic pollution on parasympathetic tone and sympathetic tone marker, LF/HF. An IQR increase in BC (reference pollutant of latent traffic pollution) 48hr previsit was associated with a −44.3% (95% PI: 67.7% to −4.2%) change in mean parasympathetic tone (reference marker: HF) among participants with diabetes, while for the nondiabetic participants the association was a 7.7% change in mean parasympathetic tone (95% PI: 18.0% to 41.4%). The effect of traffic pollution on parasympathetic tone was also stronger for diabetic versus nondiabetic and obese versus nonobese participants (for details, see Table 4). Furthermore, among obese participants, the association of parasympathetic tone with traffic pollution was stronger for diabetic participants versus nondiabetic participants at all exposure periods (Figure 3A). However, among diabetic participants, the association was similar for obese participants versus nonobese participants across all examined exposure windows, except for the 72hr moving average where the effect was two times stronger for obese versus nonobese participants (Figure 3B). In addition, Figure 3C compares the effect of traffic on parasympathetic tone among participants that are both diabetic and obese versus participants that are both nondiabetics and nonobese and the effects were stronger for participants that are both diabetic and obese (for the effect of traffic on sympathetic tone marker (LF/HF) effect modification results, see Additional file 3 for details).
Model comparison
Ordinary LMMs were fitted to doublecheck the estimates and posterior intervals from the BLMMs. The results from BLMMs and LMMs produced almost identical estimates and 95% intervals (results not shown here). Figure 4 compares the posterior mean and 95% PI from SEMs to the mean and 95% PI from BLMMs. For Parasympatheric tone, results generally showed larger effect estimates (more negative posterior means) for BLMMs as compared to SEMs. For example in a 4hr previsit exposure, an IQR increase in BC was associated with changes of −4.8% (95% PI: 14.8 to 6.3%) in HF for BLMM as compare to a 0.1% change (95% PI: 19.0 to 25.0%) in Parasympathetic tone for an IQR increase in BC (reference pollutant for traffic pollution) for SEM. For sympathetic tone, comparison of the SEMs and BLMMs results showed that the effect estimates (posterior means) were all positive for both models. In contrast to that observed for the SEMs, BLMM showed a stronger association between 4hr BC and LF/HF, with an increase of 10.1% (95% PI: 3.0 to 17.6%) in LF/HF per IQR increase in BC. Further, estimated variances of the SEMs were larger compared to the variances of the BLMMs, as evidenced by the wider estimated posterior intervals for SEMs as compared to BLMMs (see Table 3 for details).
Discussion
In using SEMs, we found traffic pollution, as measured by BC, CO, NO, and NO_{2}, to have null association with parasympathetic tone, as measured by HF, SDNN, and rMSSD, and weak positive association with sympathetic tone marker, LF/HF, at all the exposure windows. Our null association of traffic and parasympathetic tone, and the weak positive association of traffic and LF/HF are in agreement with the previous study of Park and coworkers that found no significant association between HRV markers and trafficrelated pollutants, BC, NO_{2}, and CO for any of the exposure averaging periods [11]. In addition, the previous study of Weichenthal and et al. did not also observe significant associations between BC and HRV [15] and from a controlled crossover study of diesel exhaust exposure [16]. However, our results are contrary to those from previous studies that showed associations between markers of trafficrelated pollutants and HRV markers [4, 17–22]. For example, Schwartz and coworkers found that exposure to trafficrelated pollutant (e.g. BC) was associated with decreased SDNN and rMMSD response and associated with increased LF/HF response [4]. In addition, impacts assessed using the Bayesian linear mixed models (BLMMs) were associated with decreased parasympathetic tone marker, HF, and associated with increased sympathetic tone marker, LF/HF, and were consistent across all examined exposure windows, as evidenced by their respective estimated means (Table 3). Moreover, in the linear mixed modeling approach, we found a strong positive relationship between increased LF/HF and exposure to BC at the 4hr (Posterior probability that ρ _{ 1 } > 0 = 0.99) and 24hr (Posterior probability ρ _{ 1 } > 0 = 0.96) exposure windows, which was in agreement with the Schwartz et al. study that also showed a strong and positive association between BC and LF/HF [4].
The differences in the result we observed between the SEM approach and the Bayesian linear mixed model approach implies that the effect of traffic on parasympathetic is not the same as the effect of BC on HF (Table 3). The dissimilarity in results may be explained by the difference between traffic in general, which includes, CO, NOx from gas combustion vehicles, with modest BC emitted, and diesel particulate matter that is most closely associated with the BC variable. In 2010, SpiraCohen and coworkers published a study that looked at the analyses of urban elemental carbon samples, although the study was not in Boston, the study found that over 90% was due to diesel emissions in particular [23]. The results from the Bayesian linear mixed model approach and the SEM approach would also imply the higher toxicity of diesel traffic emissions versus gasoline powered vehicle emissions. Our results showed that the examination of parasympathetic effects of traffic pollution is not needed for all the participants. As for the sympathetic tone outcome, in the SEM approach, traffic pollution was associated with increased LF/HF at all exposure windows. Moreover, in the linear mixed model approach, BC was also associated with increased LF/HF with the 4hr and 24hr exposure windows being strongly associated (Posterior probability that ρ _{ 1 } > 0 is greater than 0.95).
Several works in the past have used factor analysis, a key component and the measurement model part of structural equation models, to relate multiple exposure variables with health variables. For example, the work of the EPA Workshop on Source Apportionment of Particulate Matter Health Effects summarized the results of time series analyses of health effects with source factors [24]. However, in our approach, we used both the measurement model part and the structural model part of SEMs. We created latent variables for traffic (exposure) and parasympathetic tone (outcome) and modeled the association dependencies between these two latent constructs.
Correspondingly, our observation of effect modification by diabetic status in the association between traffic pollution and parasympathetic tone may be due to a higher likelihood of diabetics to have decreased parasympathetic tone than nondiabetics [25] and/or to their greater baseline oxidative stress, which has been shown to mediate the autonomic effects of particles [26–28]. SEMs at all exposure windows showed a consistent pattern in the acute estimated effect of traffic on parasympathetic tone in diabetics. These results are supported by those from a populationbased study that showed a high prevalence of decreased parasympathetic tone in type 2 diabetic participants [25]. Additionally, we observed evidence of effect modification by obesity on the association between traffic pollution exposure and decreased parasympathetic tone. Cardiac autonomic nervous system and obesity are related, Poirier et al. in 2006 reported that an increase in body weight is associated with a decline in parasympathetic tone and accompanied by a rise in mean heart rate [29]. This may explain why we found a stronger association between traffic pollution and parasympathetic tone among obese participants as compared to nonobese participants. In addition, when we stratified the participants to obese and nonobese, we still observed evidence of effect modification by diabetic status on the effect of traffic pollution on parasympathetic tone. However, no effect modification by obesity status was observed when we stratified the participants to diabetic and nondiabetic. Since only 8% of the participants were both obese and diabetic, it is unlikely in this case that obesity is acting as a surrogate for diabetes. Furthermore, we also observed evidence of effect modification by being both obese and diabetic versus nonobese and nondiabetic on the association between traffic pollution exposure and decreases parasympathetic tone (Table 4).
Our study to our knowledge is the first to demonstrate the simultaneous impacts of multiple trafficrelated pollutants on parasympathetic tone and sympathetic tone marker LF/HF, cardiovascular health outcomes that have the capability to predict cardiac morbidity and sudden cardiac death. Our study showed that multiple markers of traffic pollution and health impacts could be analyzed concurrently using SEMs. Moreover, in using Bayesian methods whether SEMs or linear mixed models, the estimates are easily interpreted and are useful quantities, an advantage of our proposed approach. Our study in a Bayesian framework that utilized SEMs was able to integrate data and to account for variations in one probabilistic framework. Additionally, SEMs have the ability to analyze unbalanced longitudinal data or repeated measures, with multiple outcomes and/or exposure markers, and benefit from a reduction in measurement error bias. In addition, our study used a random subject intercept which means that the contrasts are a mixture of within and between subjects; for this reason, confounding factor bias could be small relative to that in a purely crosssectional study design [30]. The comprehensive structural equation modeling approach, however, is not required for study designs that only have one exposure variable and one measured outcome variable, for which linear regression model and linear mixed model will still be the favored method of analysis for crosssectional studies, and repeated and longitudinal studies, respectively. Nonetheless, more studies that look at the parallel use of the two approaches (linear mixed models and SEMs) are needed to verify and validate our results.
There are a number of limitations in our study. In our study, exposures were estimated using concentrations measured at several ambient monitoring sites for CO, NO and NO_{2} and at a single ambient monitoring site for BC, which will not capture spatial variation in air pollutant concentrations. Because participants lived up to 22 km on median straight line distance from our ambient monitoring site, this spatial variation may result in exposure measurement error that would probably be nondifferential and could bias the results in either direction [31]. Additionally, bias due to unmeasured or residual confounding and predictor factors cannot be ruled out, and the potential for larger variance and wider posterior interval should also be considered if SEMs and not linear mixed models are to be used in air pollution epidemiologic studies [9, 32]. Lastly, since the study population consists of a high percentage of older ever smoker males who are mostly veterans and predominately white that have particular occupational exposure; the results may not be generalizable to women, younger individuals, or other ethnic and racial groups. The effect of traffic pollution on parasympathetic tone and LF/HF on these other populations should be tackled in future studies.
Conclusions
In conclusion, this study documented the null association between exposure to traffic pollution and parasympathetic tone and the weak association between traffic pollution exposure and sympathetic tone marker, LF/HF, health effect outcomes that have the capability to predict cardiovascular disease, among older men. In addition, BC is adversely associated with LF/HF in the elderly and that traffic pollution may decrease parasympathetic tone among diabetic elderly. The study also showed that SEMs in a Bayesian framework would be a plausible alternative method for unbalanced longitudinal and repeated measures studies in air pollution epidemiology that look at the cumulative effect of traffic pollution exposure on multiple health effect outcomes.
Abbreviations
 AT:

Ambient temperature
 ATemp:

Apparent temperature
 BC:

Black carbon
 BLMM:

Bayesian linear mixed model
 BMI:

Body mass index
 CO:

Carbon monoxide
 DAG:

Directed acyclic graph
 DM:

Diabetes mellitus
 DPT:

Dew point temperature
 HF:

High frequency
 HRV:

Heart rate variability
 IQR:

Inter quartile range
 LMM:

Frequentist linear mixed model
 LF:

Low frequency
 LF/HF:

Low frequency/high frequency ratio
 NAS:

Normative Aging Study
 NO2:

Nitrogen dioxide
 NO:

Nitrogen monoxide
 PI:

Posterior interval
 rMSSD:

Square root of the mean squared differences of successive normaltonormal intervals
 SDNN:

5minutes standard deviation of normaltonormal intervals
 SEM:

Structural equation model.
References
 1.
Brook RD, Rajagopalan S, Pope CA, Brook JR, Bhatnagar A, DiezRoux AV, Holguin F, Hong Y, Luepker RV, Mittleman MA, Peters A, Siscovick D, Smith SC, Whitsel L, Kaufman JD: Particulate matter air pollution and cardiovascular disease: An update to the scientific statement from the American Heart Association. Circulation. 2010, 121: 23312378. 10.1161/CIR.0b013e3181dbece1.
 2.
LuttmannGibson H, Suh HH, Coull BA, Dockery DW, Sarnat SE, Schwartz J, Stone PH, Gold DR: Systemic inflammation, heart rate variability and air pollution in a cohort of senior adults. Occup Environ Med. 2010, 67: 625630. 10.1136/oem.2009.050625.
 3.
Pope CA, Hansen ML, Long RW, Nielsen KR, Eatough NL, Wilson WE, Eatough DJ: Ambient particulate air pollution, heart rate variability, and blood markers of inflammation in a panel of elderly subjects. Environ Health Perspect. 2004, 112: 339345. 10.1289/ehp.112a339.
 4.
Schwartz J, Litonjua A, Suh H, Verrier M, Zanobetti A, Syring M, Nearing B, Verrier R, Stone P, MacCallum G: Traffic related pollution and heart rate variability in a panel of elderly subjects. Thorax. 2005, 60: 455461. 10.1136/thx.2004.024836.
 5.
Nikolov MC, Coull BA, Catalano PJ, Godleski JJ: An informative Bayesian structural equation model to assess sourcespecific health effects of air pollution. Biostatistics. 2007, 8: 609624. 10.1093/biostatistics/kxl032.
 6.
McBride SJ, Norris GA, Williams RW, Neas LM: Bayesian hierarchical modeling of cardiac response to particulate matter exposure. J Expo Sci Environ Epidemiol. 2011, 21: 7491. 10.1038/jes.2009.58.
 7.
Grandjean P, BudtzJorgensen E, Jorgensen PJ, Weihe P: Umbilical cord mercury concentration as biomarker of prenatal exposure to methylmercury. Environ Health Perspect. 2005, 113: 905908. 10.1289/ehp.7842.
 8.
SÃ¡nchez BN, BudtzJÃ¸rgensen E, Ryan LM, Hu H: Structural Equation Models. J Am Stat Assoc. 2005, 100: 14431455. 10.1198/016214505000001005.
 9.
Bollen K: Structural Equations with Latent Variables. 1989, New York City: Wiley, 2:10–39, 5:159–168, 6:179:183, 8:319–322
 10.
Bell B, Rose C, Damon A: The Normative Aging Study: an interdisciplinary and longitudinal study of health and aging. Aging Human Develop. 1972, 3: 417.
 11.
Park SK, O’Neill MS, Vokonas PS, Sparrow D, Schwartz J: Effects of air pollution on heart rate variability: the VA normative aging study. Environ Health Perspect. 2005, 113: 304309. 10.1289/ehp.113a304.
 12.
Pope CA, Eatough DJ, Gold DR, Pang Y, Nielsen KR, Nath P, Verrier RL, Kanner RE: Acute exposure to environmental tobacco smoke and heart rate variability. Environ Health Perspect. 2001, 109: 711716.
 13.
Kalkstein LS, Valimont KM: An Evaluation of Summer Discomfort in the United State Using a Relative Climatological Index. Bull Am Meteorol Soc. 1986, 67: 842848. 10.1175/15200477(1986)067<0842:AEOSDI>2.0.CO;2.
 14.
Pearl J: Causal Diagrams for Empirical Research. Biometrika. 1995, 82: 669688.
 15.
Weichenthal S, Kulka R, Dubeau A, Martin C, Wang D, Dales R: Trafficrelated air pollution and acute changes in heart rate variability and respiratory function in urban cyclists. Environ Health Perspect. 2011, 119: 13738. 10.1289/ehp.1003321.
 16.
Peretz A, Kaufman JD, Trenga CA, Allen J, Carlsten C, Aulet MR, Adar SD, Sullivan JH: Effects of diesel exhaust inhalation on heart rate variability in human volunteers. Environ Res. 2008, 107: 17884. 10.1016/j.envres.2008.01.012.
 17.
Adar SD, Gold DR, Coull BA, Schwartz J, Stone PH, Suh H: Focused Exposures to Airborne Traffic Particles and Heart Rate Variability in the Elderly. Epidemiology. 2007, 18: 95103. 10.1097/01.ede.0000249409.81050.46.
 18.
Liao D, Duan Y, Whitsel EA, Zheng ZJ, Heiss G, Chinchilli VM, Lin HM: Association of higher levels of ambient criteria pollutants with impaired cardiac autonomic control: a populationbased study. Am J Epidemiol. 2004, 159: 76877. 10.1093/aje/kwh109.
 19.
Chan CC, Chuang KJ, Shiao GM, Lin LY: Personal exposure to submicrometer particles and heart rate variability in human subjects. Environ Health Perspect. 2004, 112: 10637. 10.1289/ehp.6897.
 20.
Gold DR, Litonjua A, Schwartz J, Lovett E, Larson A, Nearing B, Allen G, Verrier M, Cherry R, Verrier R: Ambient pollution and heart rate variability. Circulation. 2000, 101: 126773. 10.1161/01.CIR.101.11.1267.
 21.
Liao D, Creason J, Shy C, Williams R, Watts R, Zweidinger R: Daily variation of particulate air pollution and poor cardiac autonomic control in the elderly. Environ Health Perspect. 1999, 107: 5215. 10.1289/ehp.99107521.
 22.
Zanobetti A, Gold DR, Stone PH, Suh HH, Schwartz J, Coull BA, Speizer FE: Reduction in heart rate variability with traffic and air pollution in patients with coronary artery disease. Environ Health Perspect. 2010, 118: 324330.
 23.
SpiraCohen A, Chen LC, Kendall M, Sheesley R, Thurston GD: Personal exposures to trafficrelated particle pollution among children with asthma in the South Bronx, NY. J Expo Sci Environ Epidemiol. 2010, 20: 446456. 10.1038/jes.2009.34.
 24.
Thurston GD, Ito K, Mar T, Christensen WF, Eatough DJ, Henry RC, Kim E, Laden F, Lall R, Larson TV: Workgroup Report: Workshop on Source Apportionment of Particulate Matter Health EffectsIntercomparison of Results and Implications. Environ Health Perspect. 2005, 113: 17681774. 10.1289/ehp.7989.
 25.
Liao D, Cai J, Brancati FL, Folsom A, Barnes RW, Tyroler HA, Heiss G: Association of vagal tone with serum insulin, glucose, and diabetes mellitus–The ARIC Study. Diabetes Res Clin Pract. 1995, 30: 211221. 10.1016/01688227(95)011900.
 26.
Chahine T, Baccarelli A, Litonjua A, Wright RO, Suh H, Gold DR, Sparrow D, Vokonas P, Schwartz J: Particulate air pollution, oxidative stress genes, and heart rate variability in an elderly cohort. Environ Health Perspect. 2007, 115: 16171622. 10.1289/ehp.10318.
 27.
Schwartz J, Park SK, O’Neill MS, Vokonas PS, Sparrow D, Weiss S, Kelsey K: GlutathioneStransferase M1, obesity, statins, and autonomic effects of particles: genebydrugbyenvironment interaction. Am J Respir Crit Care Med. 2005, 172: 15291533. 10.1164/rccm.2004121698OC.
 28.
Rhoden CR, Wellenius GA, Ghelfi E, Lawrence J, GonzalezFlecha B: PMinduced cardiac oxidative stress and dysfunction are mediated by autonomic stimulation. Biochim Biophys Acta. 2005, 1725: 305313. 10.1016/j.bbagen.2005.05.025.
 29.
Poirier P, Giles TD, Bray GA, Hong Y, Stern JS, PiSunyer FX, Eckel RH: Obesity and cardiovascular disease: pathophysiology, evaluation, and effect of weight loss: an update of the, American Heart Association Scientific Statement on Obesity and Heart Disease from the Obesity Committee of the Council on Nutrition, Physical Activity, and Metabolism. Circulation. 1997, 2006 (113): 898918.
 30.
Fitzmaurice GM, Laird NM, Ware JH: Applied Longitudinal Analysis. 2004, Hoboken, NJ: WileyInterscience, 402408. Volume 14
 31.
Sarnat JA, Brown KW, Schwartz J, Coull BA, Koutrakis P: Ambient gas concentrations and personal particulate matter exposures: implications for studying the health effects of particles. Epidemiology. 2005, 16: 385395. 10.1097/01.ede.0000155505.04775.33.
 32.
Hoyle RH: The Structural Equation Modeling Approach: Basic concepts and fundamental issues. Structural equation modeling: Concepts, issues, and applications. Edited by: Hoyle RH. 1995, Thousand Oaks, CA: Sage, 115. Volume 1
Acknowledgments
This work was supported by grants from the National Institute of Environmental Health Sciences (ES01466301A2, ES015172, ES015774, and PO1 ES09825) and from the U.S. Environmental Protection Agency (EPA R827353 and R83241601).
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
ESB designed the study with support from HHS, JDS, BAC, GAW and PSV, conducted the analysis under the guidance of BAC, and wrote the manuscript with assistance from HHS. JDS, BAC, GAW, and PSV reviewed the drafts. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Baja, E.S., Schwartz, J.D., Coull, B.A. et al. Structural equation modeling of parasympathetic and sympathetic response to traffic air pollution in a repeated measures study. Environ Health 12, 81 (2013) doi:10.1186/1476069X1281
Received
Accepted
Published
DOI
Keywords
 Bayesian
 Diabetes
 HRV
 Obesity
 Parasympathetic response
 Structural equation models
 Sympathetic response and Traffic air pollution
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.