Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Spatial environmental factors predict cardiovascular and all-cause mortality: Results of the SPACE study

  • Michael B. Hadley ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing

    michael.hadley@mountsinai.org

    Affiliation Icahn School of Medicine at Mount Sinai, New York, New York, United States of America

  • Mahdi Nalini,

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliations Digestive Disease Research Center, Digestive Disease Research Institute, Tehran University of Medical Sciences, Tehran, Iran, Cardiovascular Research Center, Kermanshah University of Medical Sciences, Kermanshah, Iran

  • Samrachana Adhikari,

    Roles Conceptualization, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

    Affiliation New York University Grossman School of Medicine, New York, New York, United States of America

  • Jackie Szymonifka,

    Roles Formal analysis, Writing – review & editing

    Affiliation New York University Grossman School of Medicine, New York, New York, United States of America

  • Arash Etemadi,

    Roles Conceptualization, Writing – review & editing

    Affiliation Metabolic Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, Maryland, United States of America

  • Farin Kamangar,

    Roles Conceptualization, Formal analysis, Methodology, Writing – review & editing

    Affiliation Department of Biology, School of Computer, Mathematical, and Natural Sciences, Morgan State University, Baltimore, Maryland, United States of America

  • Masoud Khoshnia,

    Roles Conceptualization, Data curation, Writing – review & editing

    Affiliation Golestan University of Medical Sciences, Gorgan, Golestan, Iran

  • Tyler McChane,

    Roles Formal analysis

    Affiliation Icahn School of Medicine at Mount Sinai, New York, New York, United States of America

  • Akram Pourshams,

    Roles Conceptualization, Data curation, Writing – review & editing

    Affiliation Digestive Disease Research Center, Digestive Disease Research Institute, Tehran University of Medical Sciences, Tehran, Iran

  • Hossein Poustchi,

    Roles Conceptualization, Data curation

    Affiliation Digestive Disease Research Center, Digestive Disease Research Institute, Tehran University of Medical Sciences, Tehran, Iran

  • Sadaf G. Sepanlou,

    Roles Conceptualization, Data curation, Formal analysis, Writing – review & editing

    Affiliation Digestive Disease Research Center, Digestive Disease Research Institute, Tehran University of Medical Sciences, Tehran, Iran

  • Christian Abnet,

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Metabolic Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, Maryland, United States of America

  • Neal D. Freedman,

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Metabolic Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, Maryland, United States of America

  • Paolo Boffetta,

    Roles Conceptualization, Formal analysis, Methodology, Writing – review & editing

    Affiliations Stony Brook Cancer Center, Stony Brook University, Stony Brook, New York, United States of America, Department of Medical and Surgical Sciences, University of Bologna, Bologna, Italy

  • Reza Malekzadeh,

    Roles Conceptualization, Data curation, Methodology, Writing – review & editing

    Affiliation Digestive Disease Research Center, Digestive Disease Research Institute, Tehran University of Medical Sciences, Tehran, Iran

  •  [ ... ],
  • Rajesh Vedanthan

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing

    Affiliation New York University Grossman School of Medicine, New York, New York, United States of America

  • [ view all ]
  • [ view less ]

Abstract

Background

Environmental exposures account for a growing proportion of global mortality. Large cohort studies are needed to characterize the independent impact of environmental exposures on mortality in low-income settings.

Methods

We collected data on individual and environmental risk factors for a multiethnic cohort of 50,045 individuals in a low-income region in Iran. Environmental risk factors included: ambient fine particular matter air pollution; household fuel use and ventilation; proximity to traffic; distance to percutaneous coronary intervention (PCI) center; socioeconomic environment; population density; local land use; and nighttime light exposure. We developed a spatial survival model to estimate the independent associations between these environmental exposures and all-cause and cardiovascular mortality.

Findings

Several environmental factors demonstrated associations with mortality after adjusting for individual risk factors. Ambient fine particulate matter air pollution predicted all-cause mortality (per μg/m3, HR 1.20, 95% CI 1.07, 1.36) and cardiovascular mortality (HR 1.17, 95% CI 0.98, 1.39). Biomass fuel use without chimney predicted all-cause mortality (reference = gas, HR 1.23, 95% CI 0.99, 1.53) and cardiovascular mortality (HR 1.36, 95% CI 0.99, 1.87). Kerosene fuel use without chimney predicted all-cause mortality (reference = gas, HR 1.09, 95% CI 0.97, 1.23) and cardiovascular mortality (HR 1.19, 95% CI 1.01, 1.41). Distance to PCI center predicted all-cause mortality (per 10km, HR 1.01, 95% CI 1.004, 1.022) and cardiovascular mortality (HR 1.02, 95% CI 1.004, 1.031). Additionally, proximity to traffic predicted all-cause mortality (HR 1.13, 95% CI 1.01, 1.27). In a separate validation cohort, the multivariable model effectively predicted both all-cause mortality (AUC 0.76) and cardiovascular mortality (AUC 0.81). Population attributable fractions demonstrated a high mortality burden attributable to environmental exposures.

Interpretation

Several environmental factors predicted cardiovascular and all-cause mortality, independent of each other and of individual risk factors. Mortality attributable to environmental factors represents a critical opportunity for targeted policies and programs.

Introduction

Environmental factors contribute significantly to global mortality [14]. In 2019, environmental hazards were responsible for an estimated 11.3 million deaths, of which 5.1 million were from cardiovascular disease (CVD) [5, 6].

A growing list of environmental factors present particular risks to cardiovascular health [14]. Ambient fine particulate matter air pollution (PM2.5) from traffic, industry, fires, and dust is a risk factor for all-cause mortality, cardiovascular mortality, ischemic heart disease (IHD), and stroke [7, 8]. In 2019, ambient PM2.5 ranked seventh among all health risk factors for mortality, responsible for 4.14 million deaths, of which 2.47 million were from CVD [5, 6]. Similarly, household air pollution from inefficient stoves and solid fuels was responsible for 2.31 million deaths (1.07 million from CVD) [5, 6]. Proximity to traffic pollution and noise is associated with increased rates of adverse CVD events, particularly IHD and stroke [4, 7, 9]. Distance to health care services affects access to preventive care, tertiary care, and emergent percutaneous coronary interventions (PCI) [10, 11]. Socioeconomic environment is an independent risk factor for CVD, even after controlling for individual socioeconomic status [12]. Population density has been shown to be both positively and negatively correlated with CVD in different environments [13, 14]. Exposures to artificial light at night cause circadian dysregulation and have been associated with ischemic heart disease outcomes [4, 15, 16]. Finally, land use (e.g., greenspaces, mixed land use) may predict CVD through association with physical activity, social engagement, and access to health services [17, 18]. Together, these variables act directly and indirectly to precipitate cardiovascular disease and mortality (Fig 1).

thumbnail
Fig 1. Conceptual framework.

Impact of eight SEFs on individual risk factors and CV and all-cause mortality.

https://doi.org/10.1371/journal.pone.0269650.g001

Understanding the relationships between environmental risk factors and health is a critical step towards designing targeted policies and programs to reduce the immense burden of attributable disease. To date, most investigations of environmental factors have studied single risk factors in high-income settings. We therefore developed a spatial environmental model to provide a better understanding of the independent associations between multiple spatial environmental factors and mortality, within a low-income population in the middle-income country of Iran.

Methods

Study setting

We analyzed data of participants enrolled in the Golestan Cohort Study (GCS) in Golestan Province, Iran, a middle-income country with diverse ethnicities and lifestyles. In both Iran and Golestan, CVD is the leading cause of death and disability [19].

Study participants

The GCS enrolled 50,045 individuals (28,811 females and 21,234 males) across northeastern Golestan from 2004 to 2008 [20]. Participants ranged in age from 40 to 75 years in order to capture individuals with higher rates of non-communicable disease, particularly esophageal cancer. Approximately 80% were enrolled from 326 rural villages ranging in size from 20 to 150 residents. The remaining 20% were selected randomly from Gonbad City, the second-largest city in Golestan with a population of approximately 130,000. Exclusion criteria were: unwillingness to participate for any reason; being a temporary resident; or having a previous diagnosis of upper gastrointestinal cancer. Among those selected to enroll, participation rates were approximately 80% for women and 65% for men. Participants were followed-up actively every 12 months with a follow-up success rate was 99%. Study methods were approved by ethics review committees of the Tehran University of Medical Sciences, the International Agency for Research on Cancer, and the National Cancer Institute. All participants signed a written informed consent at enrollment. The full study protocol is publicly available [20].

Individual characteristics

Individual baseline characteristics were collected at time of enrollment in the GCS [20]. Participants were interviewed by a physician in their native language and completed a detailed lifestyle questionnaire and physical exam. The following covariates were included in our analysis: age, sex, ethnicity, marital status, education, socioeconomic status, waist and hip circumference, physical activity, medical history (IHD, stroke, diabetes, hypertension), and substance use (tobacco, alcohol, and opium). Details on the socioeconomic status score can be found in S1 File.

At baseline, individuals enrolled in the GCS were, on average, 52.1 years of age (SD 8.9 years) and 58% female. The majority was married (88%) and illiterate (70%). A minority had a history of IHD or stroke (6%), diabetes (7%), hypertension (20%), tobacco use (22%), opium use (17%), or alcohol use (3%). Three-quarters were of Turkmen ethnicity and the remainder were of other ethinc groups (e.g., Persians, Baluchis, Qizilbash), consistent with the prevalence of ethnic groups in the sampled region. The derivation and validation cohorts were similar in the distribution of risk factors (Table 1).

Mortality data

In the GCS, nearly 100% of deaths are captured via documents collected from participants, hospital files, and verbal autopsy questionnaires, which are reviewed by at least two independent internists to define diagnoses according to the International Classification of Diseases, 10th Revision (ICD-10) codes [21, 22]. All-cause mortality included all reported deaths. Cardiovascular mortality included deaths attributable to IHD (ICD-10 codes I20-25), cerebrovascular disease (I60-69), cardiac arrest (I46), congestive heart failure (I50), hypertensive diseases (I10-15), chronic rheumatic heart diseases (I05-09), pulmonary heart disease (I26-28), and other cardiovascular system diseases not otherwise specified. Complete case analysis was used for missing data.

Individual geocodes

Each individual was assigned a geocode (latitude and longitude) based on residence location. Details on geocode assignment can be found in S1 File.

Spatial environmental factors

We developed exposure variables for eight prespecified spatial environmental factors (SEFs) using a combination of GCS data and publicly available datasets [20, 2328]. These variables were chosen based on a review of the literature on environmental risks for cardiovascular disease [14], as well as data availability. The SEFs were: ambient air pollution [24, 25], household fuel use and ventilation [20], socioeconomic environment [23], proximity to traffic (within 100m of a minor highway or within 500m of a major highway) [20], distance to percutaneous coronary intervention centers [20], population density [26], nighttime light exposure [27], and land use [28]. Environmental exposures were assigned according to year of enrollment. Some potential environmental hazards were not included because no data was available for the study region (e.g., noise pollution; toxins in water and food). All SEFs were included in the final multivariable model. Data sources and related methodologies are summarized in S1 Table in S1 File. Please see Data Sharing Statement for information on how to access GCS-specific data.

Derivation and validation cohorts

The GCS dataset was randomized into derivation and validation cohorts, stratified to ensure mortality was balanced in both groups, following established methodology [29]. The derivation cohort (90%) was used to construct the multivariable spatial environmental model; the validation cohort (10%) was used subsequently to test the model’s predictive value. The derivation cohort contained 45,042 individuals and 5996 deaths, of which 2733 were cardiovascular deaths. The validation cohort contained 5003 individuals and 655 deaths, of which 286 were cardiovascular deaths. A summary of the study design is shown in S1 Fig in S1 File.

Statistical analyses

We developed spatial environmental survival models to study the association of each SEF on the hazard of mortality, for both all-cause and cardiovascular mortality. All eight SEFs were tested simultaneously in the multivariable survival model, in order to adjust for each other. Additionally, we adjusted for common individual risk factors, including sex, age, individual socioeconomic status, anthropometric measures, and history of cardiovascular disease, hypertension, diabetes, and smoking. We used a spatial random effects survival model [3032] (e.g., shared frailty model) to model time to mortality. For an individual j in geocode i with censoring time t given covariates xij,, the survival function in the frailty model is:

Here, β is a vector of regression coefficients and νi is a frailty parameter (or the random effect) for geocode i and is assumed to have a gaussian distribution. Environmental exposures demonstrated low correlation at the geocode level reducing the risk of collinearity and model overfitting. Computations were performed in R (packages ‘survival’ and ‘spBayesSurv’) [33]. We report the exponentiated coefficients as the estimate of hazard ratio along with 95% confidence interval of the estimates.

To account for spatial dependence, as a sensitivity analysis we also fitted the Bayesian survival models that adjust for spatial autocorrelation based on distance between geocodes [34, 35]. Gaussian random field priors were specified on the frailty parameters to allow for autocorrelation between neighboring geocodes.

We validated the model using recent methodologies for external validation [36]. First, we used a chi-squared test to determine whether the addition of SEFs improved the model’s goodness-of-fit beyond traditional risk factors in the derivation cohort. Next, time-dependent area-under-the-curve analysis was used to test how well the spatial frailty model predicted all-cause and cardiovascular mortality in the novel validation cohort.

Finally, we calculated the population attributation fraction (PAF) for all categorical predictors in the multivariable model. PAFs incorporate both the hazard and prevalence of a risk factor to assess the fraction of total disease risk in the population that would be eliminated (or added) if the risk factor were eliminated from the population (i.e., if all individuals with that risk factor were moved to the reference category) [37].

Patient and public involvement

Local residents, physicians, elders, religious leaders, and university physicians were deeply involved in the design and implementation of the Golestan Cohort Study. Details can be found in S1 File.

Results

Distribution of spatial environmental factors

Across all geocodes, average annual fine particulate matter air pollution exposures were 33.5 μg/m3 for the 5 years prior to enrollment (167.7 μg/m3 cumulatively, SD 17.5 μg/m3). Most households burned kerosene fuel (71%), of which 42% had a chimney for ventilation. A total of 7% of households used biomass fuels (typically wood or dung burned indoors for cooking or heating), of which 81% had a chimney. The remainder of households providing responses used either gas (12%) or mixed fuels (9%). One third of participants (34%) lived close to major highways. Distances to the nearest percutaneous coronary intervention center averaged 92.2 km (SD 37.4 km). The socioeconomic status score varied with location, with lower scores concentrated in the northeast. Local population density averaged 1732 persons per square kilometer (SD 3069). The intensity of nighttime light averaged 22.7 (SD 22.2) on a NOAA intensity metric. Both population density and nighttime light exposure were highest in Golestan’s central agricultural valley. Most households were located amidst cropland (57%) or urban settings (25%), with the remainder located among shrubland (9%), grassland (9%), or barren earth (<1%). The derivation and validation cohorts were similar in the distribution of all SEFs (Table 2). The spatial distribution of the SEFs is illustrated in Fig 2.

thumbnail
Fig 2. Spatial models for six spatial environmental factors (SEFs) across Golestan Province, Iran.

These models were used to assign environmental exposures to individuals based on their location of residence. Methodologies are described in S1 Table in S1 File. A. Land use based on satellite imagery [28]. B. Proportion of households burning biomass fuels without a chimney. Estimates derived from Gaussian process regression (kriging) of average fuel use patterns for each village geocode [20]. C. Proximity to traffic, defined as within 100m of a minor highway or within 500m of a major highway [20]. D. Proximity to percutaneous coronary intervention centers [20]. E. Average socioeconomic status index score. Estimates derived from kriging of median socioeconomic score for each geocode [23]. F. Average annual intensity of light-at-night using satellite imagery [27]. Models developed in ArcGIS [38]. All displayed data is either in the public domain, collected by the investigators, or for illustrative purposes only.

https://doi.org/10.1371/journal.pone.0269650.g002

thumbnail
Table 2. Individual exposures to spatial environmental factors.

https://doi.org/10.1371/journal.pone.0269650.t002

Mortality data

In the derivation cohort (n = 45,042), there were 2733 cardiovascular deaths and 5996 all-cause deaths. In the validation cohort (n = 5003), there were 286 cardiovascular deaths and 655 all-cause deaths. The mean follow-up time in the derivation and validation cohorts was 10.2 years and 10.9 years, respectively. Mean time to all-cause death in the derivation and validation cohorts was 6.2 years and 6.1 years, respectively. Mean time to cardiovascular death in both cohorts was 5.8 years. In both derivation and validation cohorts, all-cause mortality was 13% and cardiovascular mortality was 6% over the follow-up period.

Individual characteristics predict mortality

Individual characteristics predicted all-cause and cardiovascular mortality in the derivation cohort (Table 3). The following variables demonstrated increased hazard for both all-cause and cardiovascular mortality: older age, male gender, being unmarried, lower socioeconomic status, illiteracy, higher waist circumference (adjusted for hip circumference), lower hip circumference (adjusted for waist circumference), lower physical inactivity, tobacco use, opium use, and history of hypertension, diabetes, IHD, and stroke. Turkmen ethnicity was associated with increased hazard for all-cause mortality, but not for cardiovascular mortality.

Spatial environmental factors predict mortality

We identified several SEFs that predicted mortality in the multivariable derivation model (Table 4). (Univariate estimates are provided in S2 Table in S1 File).

thumbnail
Table 4. Hazard ratios for eight spatial environmental factors.

https://doi.org/10.1371/journal.pone.0269650.t004

In the multivariable model, exposure to outdoor air pollution predicted all-cause mortality (per μg/m3, HR 1.20, 95% CI 1.07 to 1.36) and cardiovascular mortality (per μg/m3, HR 1.17, 95% CI 0.98 to 1.39).

Exposure to household air pollution also demonstrated predictive power. Specifically, the use of household kerosene fuel use without a chimney compared to gas predicted cardiovascular mortality (HR 1.19, 95% CI 1.01 to 1.41) and all-cause mortality (HR 1.09, 95% CI 0.97 to 1.23). Additionally, the use of biomass fuel without a chimney compared to gas predicted both cardiovascular mortality (HR 1.36, 95% CI 0.99 to 1.87) and all-cause mortality (HR 1.23, 95% CI 0.99 to 1.53).

Greater distance to PCI centers also was associated with increased hazard of both cardiovascular mortality (per 10 km, HR 1.02, 95% CI 1.004 to 1.03) and all-cause mortality (per 10 km, HR 1.01, 95% CI 1.004 to 1.02). Proximity to traffic also increased the hazard of cardiovascular mortality (HR 1.13, 95% CI 1.01 to 1.27).

The remaining SEFs—neighborhood socioeconomic status, local population density, nighttime light, and land use—did not demonstrate relationships with either all-cause or cardiovascular mortality.

Spatial environmental factors add predictive power

The addition of the eight SEFs improved the predictive power of the model beyond traditional risk factors. Goodness of fit statistics (chi-squared statistics and corresponding p-values) were used to compare the derivation model incorporating environmental variables to the reduced model. For both cardiovascular and all-cause mortality, the addition of SEFs improved goodness of fit (chi-squared(df): 57.5(16), p<0.001 for cardiovascular mortality; 54.3(16), p<0.001 for all-cause mortality).

Stratification by sex

A sex-stratified analysis found no effect modification of the relationship between SEFs and either all-cause or cardiovascular mortality.

Adjusting for spatial autocorrelation

After adjustments for spatial autocorrelation, there was no meaningful change in the magnitude or direction of the point estimates for the eight SEFs. The spatial autocorrelation model was consistent with the original frailty estimates.

Model validation

The model incorporating SEFs was tested on a novel validation cohort. In this validation cohort, the spatial frailty model effectively predicted mortality with time-dependent areas-under-the-curve of 0.76 for all-cause mortality and 0.81 for cardiovascular mortality.

Population attributable fractions

PAFs for common individual risk factors were characteristically large, including hypertension (CV mortality 0.42; all-cause mortality 0.25), diabetes (CV 0.41; all-cause 0.35), and current tobacco use (CV 0.16; all-cause 0.12). PAFs were substantial for SEFs, including ambient PM2.5 (highest quartile, CV mortality 0.14; all-cause mortality 0.15), biomass fuel use without chimney (CV 0.24; all-cause 0.15), kerosene fuel use without chimney (CV 0.15; all-cause 0.07), and proximity to traffic (CV 0.11; all-cause 0.03). These PAFs are illustrated in S2 Fig in S1 File, along with PAFs for several common individual risk factors. PAFs for all categorical predictors can be found in S3 Table in S1 File.

Discussion

Statement of principal findings

In this analysis of the Golestan Cohort Study, several SEFs predicted cardiovascular or all-cause mortality in the spatial survival model after adjusting for individual risk factors. When applied to a novel validation cohort, the model effectively predicted both cardiovascular and all-cause mortality.

Ambient air pollution.

In our model, ambient air pollution levels were associated with all-cause and cardiovascular mortality after adjusting for other environmental risk factors in the multivariable model. Additionally, PAFs identified a large burden of mortality attributable to ambient air pollution. This is consistent with the existing literature demonstrating associations between ambient PM2.5 and both cardiovascular and all-cause mortality [37].

Household air pollution.

We found that household kerosene fuel use without a chimney was associated with both cardiovascular mortality and all-cause mortality. This is consistent with a previous study of the GCS that found associations between cumulative kerosene exposure and both cardiovascular and all-cause mortality [39]. Additionally, we observed modest evidence of an association between the use of biomass fuel without a chimney and both cardiovascular and all-cause mortality. This is consistent with the existing literature on the cardiovascular effects of indoor burning of solid biomass [3, 5, 6]. Additionally, PAF calculations illustrate a high burden of cardiovascular and all-cause mortality attributable to biomass and kerosene fuel burning.

Proximity to traffic.

In our model, residence close to highways increased the hazard of cardiovascular mortality. Although this hazard ratio was only 1.13, many individuals live near major roadways, resulting in a PAF of 0.11 for cardiovascular mortality. Importantly, this association was observed after adjusting for ambient air pollution and socioeconomic environment in the multivariable model. This suggests that the calculated hazard may be due chiefly to traffic-related noise or to another unmeasured variable associated with proximity to traffic. These findings are consistent with prior studies that have demonstrated associations between proximity to traffic and cardiovascular events, particulary IHD and stroke [7, 9].

Distance to health care services.

We found that greater distance to PCI centers increased the hazard of both cardiovascular and all-cause mortality. These results align with a previous cohort study demonstrating that distance to hospital is an independent predictor of mortality among patients with incident MI in the community [11].

Several other SEFs did not demonstrate clear relationships in the multivariable model: neighborhood socioeconomic status, local population density, nighttime light, and land use. Low socioeconomic status neighborhood environment previously was shown to predict cardiovascular events [12, 40], and did so in our univariate model for cardiovascular mortality, but not in the multivariable environmental model. This suggests that socioeconomic environment may be a proxy for other environmental risk factors (e.g., air pollution) that were included in our model.

Strengths and limitations

A key strength of our model was the simultaneous testing of multiple environmental risk factors. To date, most studies of environmental risk look at a single environmental risk factor against a background of individual characteristics. This can lead to confounding by other unmeasured environmental factors. In this study, we incorporate a diversity of spatially resolved SEFs in a prospective model predicting cardiovascular mortality. Additionally, the study was performed in a rural, low-income setting, helping to bridge a gap in the medical literature.

This study also benefitted from the use of the GCS dataset, which included a large sample size and excellent follow-up rate, as well as systematic methods to identify mortality and attributable causes of death. Additionally, our spatial random effects survival model controlled for possible spatial dependency in the data.

Our study had several limitations. First, exposures were assigned according to village- or neighborhood-level geocodes, rather than specific home addresses (due to human subjects-related privacy considerations). This could result in exposure misclassification and bias our findings towards the null. However, given the small sizes of villages, we estimate the average distance between assigned geocodes and true home address to be less than 500 meters. With the exception of proximity to traffic, none of our modeled environmental exposures vary dramatically over this distance.

Second, the spatial environmental factors in our model were assessed at the year of enrollment. This may misclassify exposure for participants exposed to SEFs prior to enrollment, for participants migrating to new locations with different exposures, or if the exposure varied over time. Similarly, given that environmental exposures were assigned at enrollment, our model does not account for acute exposures that may result in acute cardiovascular events (e.g. acute air pollution exposure triggering coronary plaque rupture). Again, we would expect these factors to bias our findings towards the null.

Third, household fuel use and ventilation were used as proxies for air pollution exposure. Although not optimal, this method has been used previously to study air pollution exposures in Golestan [39], as well as in studies on the association between household air pollution and childhood mortality [41].

Fourth, the socioeconomic environment variable used wealth indices from study participants and could be inaccurate if the mean wealth index at each geocode is not representative of the community.

Finally, although we have adjusted for eight SEFs in the environmental model, there remains the possibility of residual confounding from unmeasured SEFs (e.g. climate; temperature variation; noise pollution; toxins in water and food) and individual risk factors (e.g. dyslipidemia).

Conclusions

We tested the prospective associations between eight environmental factors and cardiovascular and all-cause mortality in a large cohort in a low-income setting. Our findings demonstrate that the burden of disease attributable to the environment may be as large as traditional cardiovascular risk factors, and thus represents a critical opportunity for targeted policies and programs. Furthermore, these findings illustrate the utility and feasibility of incorporating environmental data in survival models, even in low-income settings.

A growing literature illustrates how health care providers, governments, and charities can identify and intervene on environmental exposures at the individual and population levels [8, 42, 43]. We anticipate that these findings and analytic approach will stimulate further studies to promote better health for populations and the environment worldwide.

Acknowledgments

The authors appreciate statistical support provided by Alan Chen, Department of Population Health, NYU Grossman School of Medicine.

References

  1. 1. Bhatnagar A. Environmental determinants of cardiovascular disease. Circulation research. 2017;121(2):162–80. pmid:28684622
  2. 2. Cosselman KE, Navas-Acien A, Kaufman JD. Environmental factors in cardiovascular disease. Nat Rev Cardiol. 2015;12(11):627–42. pmid:26461967.
  3. 3. Landrigan PJ, Fuller R, Acosta NJR, Adeyi O, Arnold R, Basu NN, et al. The Lancet Commission on pollution and health. Lancet. 2017. pmid:29056410.
  4. 4. Münzel T, Hahad O, Sørensen M, Lelieveld J, Duerr GD, Nieuwenhuijsen M, et al. Environmental risk factors and cardiovascular diseases: a comprehensive review. Cardiovascular Research. 2021. pmid:34609502
  5. 5. Murray CJ, Aravkin AY, Zheng P, Abbafati C, Abbas KM, Abbasi-Kangevari M, et al. Global burden of 87 risk factors in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019. The Lancet. 2020;396(10258):1223–49. pmid:33069327
  6. 6. Institute for Health Metrics and Evaluation (2021). Global Health Data Exchange. ghdx.healthdata.org. Accessed July 20, 2021.
  7. 7. Brook RD, Rajagopalan S, Pope CA 3rd, 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(21):2331–78. pmid:20458016.
  8. 8. Kaufman JD, Elkind MS, Bhatnagar A, Koehler K, Balmes JR, Sidney S, et al. Guidance to reduce the cardiovascular burden of ambient air pollutants: a policy statement from the American Heart Association. Circulation. 2020;142(23):e432–e47. pmid:33147996
  9. 9. Rosenbloom JI, Wilker EH, Mukamal KJ, Schwartz J, Mittleman MA. Residential proximity to major roadway and 10-year all-cause mortality after myocardial infarction. Circulation. 2012;125(18):2197–203. pmid:22566348
  10. 10. Nallamothu BK, Bates ER, Wang Y, Bradley EH, Krumholz HM. Driving times and distances to hospitals with percutaneous coronary intervention in the United States: implications for prehospital triage of patients with ST-elevation myocardial infarction. Circulation. 2006;113(9):1189–95. pmid:16520425.
  11. 11. Wei L, Lang CC, Sullivan FM, Boyle P, Wang J, Pringle SD, et al. Impact on mortality following first acute myocardial infarction of distance between home and hospital: cohort study. Heart. 2008;94(9):1141–6. pmid:17984217
  12. 12. Diez Roux AV, Merkin SS, Arnett D, Chambless L, Massing M, Nieto FJ, et al. Neighborhood of residence and incidence of coronary heart disease. N Engl J Med. 2001;345(2):99–106. pmid:11450679.
  13. 13. Lindroth M, Lundqvist R, Lilja M, Eliasson M. Cardiovascular risk factors differ between rural and urban Sweden: the 2009 Northern Sweden MONICA cohort. BMC Public Health. 2014;14:825. pmid:25106746
  14. 14. Kulshreshtha A, Goyal A, Dabhadkar K, Veledar E, Vaccarino V. Urban-rural differences in coronary heart disease mortality in the United States: 1999–2009. Public Health Rep. 2014;129(1):19–29. pmid:24381356
  15. 15. Morris CJ, Purvis TE, Hu K, Scheer FA. Circadian misalignment increases cardiovascular disease risk factors in humans. Proceedings of the National Academy of Sciences. 2016;113(10):E1402–E11. pmid:26858430
  16. 16. Sun S, Cao W, Ge Y, Ran J, Sun F, Zeng Q, et al. Outdoor light at night and risk of coronary heart disease among older adults: a prospective cohort study. European Heart Journal. 2021;42(8):822–30. pmid:33205210
  17. 17. Rosso AL, Grubesic TH, Auchincloss AH, Tabb LP, Michael YL. Neighborhood amenities and mobility in older adults. Am J Epidemiol. 2013;178(5):761–9. pmid:23666814
  18. 18. Nieuwenhuijsen MJ. Influence of urban and transport planning and the city environment on cardiovascular disease. Nature reviews cardiology. 2018;15(7):432–8. pmid:29654256
  19. 19. Ahmadi A, Soori H, Mehrabi Y, Etemad K. Spatial analysis of myocardial infarction in Iran: National report from the Iranian myocardial infarction registry. J Res Med Sci. 2015;20(5):434–9. pmid:26487871
  20. 20. Pourshams A, Khademi H, Malekshah AF, Islami F, Nouraei M, Sadjadi AR, et al. Cohort Profile: The Golestan Cohort Study—a prospective study of oesophageal cancer in northern Iran. Int J Epidemiol. 2010;39(1):52–9. pmid:19332502
  21. 21. Khademi H, Etemadi A, Kamangar F, Nouraie M, Shakeri R, Abaie B, et al. Verbal autopsy: reliability and validity estimates for causes of death in the Golestan Cohort Study in Iran. PLos one. 2010;5(6):e11183. pmid:20567597
  22. 22. Organization WH. International statistical classification of diseases and related health problems: World Health Organization; 2004.
  23. 23. Islami F, Kamangar F, Nasrollahzadeh D, Aghcheli K, Sotoudeh M, Abedi-Ardekani B, et al. (2009). Socio-economic status and oesophageal cancer: results from a population-based case-control study in a high-risk area. Int J Epidemiol 38(4):978–88. pmid:19416955
  24. 24. van Donkelaar A, Martin RV, Brauer M, Boys BL. Use of satellite observations for long-term exposure assessment of global concentrations of fine particulate matter. Environ Health Perspect. 2015;123(2):135–43. pmid:25343779
  25. 25. Atmospheric Composition Analysis Group (2019). Global Estimates of All Composition Surface PM2.5. Accessed December 10, 2019. https://sites.wustl.edu/acag/datasets/surface-pm2-5/#V5.GL.02
  26. 26. Laboratory. ORN. LandScan. http://web.ornl.gov/sci/landscan/. 2017. http://web.ornl.gov/sci/landscan/
  27. 27. NOAA National Centers for Environmental Information (2016). Nighttime lights composites for the years 1992–2013 from DMSP satellite data. Accessed December 10, 2019. http://ngdc.noaa.gov/eog/dmsp/downloadV4composites.html
  28. 28. Friedl, M. & D Sulla-Menashe (2019). MCD12Q1 MODIS/Terra+Aqua Land Cover Type Yearly L3 Global 500m SIN Grid V006 [Data set]. NASA EOSDIS Land Processes DAAC. Accessed December 12, 2019. https://lpdaac.usgs.gov/products/mcd12q1v006/
  29. 29. Zhang Z, Cortese G, Combescure C, Marshall R, Lee M, Lim HJ, et al. Overview of model validation for survival regression model with competing risks using melanoma study data. Annals of translational medicine. 2018;6(16). pmid:30364028
  30. 30. Banerjee S, Carlin BP, Gelfand AE. Hierarchical modeling and analysis for spatial data: Crc Press; 2014.
  31. 31. Mitra R. Nonparametric bayesian inference in biostatistics. New York, NY: Springer Science+Business Media; 2015. pages cm p.
  32. 32. Auchincloss AH, Gebreab SY, Mair C, Diez Roux AV. A review of spatial methods in epidemiology, 2000–2010. Annu Rev Public Health. 2012;33:107–22. pmid:22429160
  33. 33. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria. https://www.R-project.org/
  34. 34. Zhou H, Hanson T, Zhang J. spbayessurv: Fitting Bayesian spatial survival models using r. arXiv preprint arXiv:170504584. 2017.
  35. 35. Zhou H, Hanson T. A unified framework for fitting Bayesian semiparametric models to arbitrarily censored survival data, including spatially referenced data. Journal of the American Statistical Association. 2018;113(522):571–81.
  36. 36. Royston P, Altman DG. External validation of a Cox prognostic model: principles and methods. BMC medical research methodology. 2013;13(1):33. pmid:23496923
  37. 37. Rockhill B, Newman B, Weinberg C. Use and misuse of population attributable fractions. American journal of public health. 1998;88(1):15–9. pmid:9584027
  38. 38. ESRI 2019. ArcGIS Desktop: Release 10.7. Redlands, CA: Environmental Systems Research Institute.
  39. 39. Mitter SS, Vedanthan R, Islami F, Pourshams A, Khademi H, Kamangar F, et al. Household Fuel Use and Cardiovascular Disease Mortality: Golestan Cohort Study. Circulation. 2016;133(24):2360–9. pmid:27297340
  40. 40. Sundquist K, Winkleby M, Ahlen H, Johansson SE. Neighborhood socioeconomic environment and incidence of coronary heart disease: a follow-up study of 25,319 women and men in Sweden. Am J Epidemiol. 2004;159(7):655–62. pmid:15033643.
  41. 41. Ezeh OK, Agho KE, Dibley MJ, Hall JJ, Page AN. The effect of solid fuel use on childhood mortality in Nigeria: evidence from the 2013 cross-sectional household survey. Environmental Health. 2014;13(1):113. pmid:25514998
  42. 42. Hadley MB, Baumgartner J, Vedanthan R. Developing a Clinical Approach to Air Pollution and Cardiovascular Health. Circulation. 2018;137(7):725–42. pmid:29440198.
  43. 43. Rajagopalan S, Brauer M, Bhatnagar A, Bhatt DL, Brook JR, Huang W, et al. Personal-Level Protective Actions Against Particulate Matter Air Pollution Exposure: A Scientific Statement From the American Heart Association. Circulation. 2020;142(23):e411–e31. pmid:33150789