Long-term exposure to PM2.5 and cardiovascular disease incidence and mortality in an Eastern Mediterranean country: findings based on a 15-year cohort study
Environmental Health volume 20, Article number: 112 (2021)
Evidence concerning the impact of long-term exposure to fine Particulate Matter ≤2.5 μm (PM2.5) on Cardio-Vascular Diseases (CVDs) for those people subject to ambient air pollution in developing countries remains relatively scant. This study assessed the relationship of 15-year PM2.5 exposure with cardiovascular incidence and mortality rate in Isfahan province, Iran.
The cohort comprised 3081 participants over 35 years old who were free of CVDs. They were selected through multi-stage cluster sampling in Isfahan, Iran. PM2.5 exposure was determined separately for each individual via satellite-based spatiotemporal estimates according to their residential addresses. In this context, CVD is defined as either fatal and non-fatal Acute Myocardial Infarctions (AMI) or stroke and sudden cardiac death. The incidence risk for CVD and the ensuing mortality was calculated based on the average PM2.5 exposure within a study period of 15 years using the Cox proportional hazards frailty model upon adjusting individual risk factors. The mean annual rate of PM2.5 and the follow-up data of each residential area were combined.
Mean three-year PM2·5 exposure for the cohort was measured at 45.28 μg/m3, ranging from 20.01 to 69.80 μg/m3. The median time period for conducting necessary follow-ups was 12.3 years for the whole population. Notably, 105 cardiovascular and 241 all-cause deaths occurred among 393,786 person-months (27 and 61 per 100,000 person-months, respectively). In well-adjusted models, 10 μg/m3 increase in PM2.5 corresponded to a 3% increase in the incidence rate of CVDs [0.95 CI = 1.016, 1.036] (in case of p = 0.000001 per 10 μg/m3 increase in PM2.5, the Hazard Ratio (HR) for AMI and Ischemic Heart Disease (IHD) was 1.031 [0.95 CI = 1.005, 1.057] and 1.028 [0.95 CI = 1.017, 1.039]), respectively. No consistent association was observed between PM2.5 concentration and fatal CVD (fatal AMI, fatal stroke, SCD (Sudden Cardiac Death)) and all-cause mortality.
Results from analyses suggest that the effect of PM2.5 on cardiovascular disease occurrence was stronger in the case of older people, smokers, and those with high blood pressure and diabetes. The final results revealed that long-term exposure to ambient PM2.5 with high concentrations positively correlated with IHD incidence and its major subtypes, except for mortality. The outcome accentuates the need for better air quality in many countries.
The current concentration of ambient air pollution may cause detrimental health effects including a high rate of mortality and morbidity, partly induced by CVD. Based on epidemiological findings specific to western countries, elevated particles with a mass median aerodynamic diameter ≤ 2.5 μm (PM2.5) consistently lead to CVD incidence and death . The fine particulate matter calculated as 2.5 μg/m3 or less is in close association with myocardial infarctions, thrombosis, and stroke . The mean of air pollutant concentrations in the USA  and European  countries is much lower than that in Eastern Asia . It should be noted that Eastern Asia is usually bombarded with frequent Asian dust events from industrial facilities and the Gobi Desert in Central Asia, reportedly producing over 20% of the total global dust emissions through the westerlies [6, 7].
The dominance of vast deserts in the surrounding areas of Isfahan, Iran together with the significant number of large industries in its suburbs constitute the reasons why the city is one of the most polluted cities in Iran . High PM2.5 levels in Isfahan could be arguably attributed to the growing emergence of deserts and mines, of lead and zinc, as well as energy conversion systems such as power plants and oil refineries in immediate surroundings. In addition, in Isfahan province, wind direction changes per season, meaning that with the exception of summer, the winds generally blow from the west to the east. A significant impact on the air quality was observed in Isfahan which inevitably resulted from the adaptation of wind currents to the degraded areas caused by the activities of gypsum, clay, sand, and mines as well as a wide range of related industries in more than 12,900 ha (6800 ha of gypsum and 6100 ha of clay and sand). The concentration of suspended particles in the air was the key to the process. Furthermore, studies that have investigated the amount of dust and suspended particles in Isfahan’s main stations suggest that the operations of local centers have yielded more than 30% increase in dust volume in the summer. As mentioned above, such particles are produced by summer winds in the east .
Over 90% of the world’s population inhabit areas whose air quality exceeds the mean annual PM2.5 standard defined at 10 μg/m3 by the World Health Organization (WHO) .
Furthermore, such concentration-responsive detrimental impacts of long-term PM2.5 exposure on cardiovascular mortality rates in countries with higher air pollution levels have not been thoroughly investigated. Therefore, it is crucial to carefully scrutinize the effects of such exposure and develop fitting preventive strategies to reduce potential risks for the community. Given the accessibility of data about individual exposure to air pollution and the significance of other individual CVD risk factors, the stated problem has been addressed in the present study.
Assessments carried out in line with the Global Burden of Disease initiative demonstrate that high systolic blood pressure, smoking, high-sodium diet, and ambient particulate matter pollution are among the leading four risk factors contributing to deaths and disability-adjusted life years across the world .
Nearly 50% of these estimated attributable deaths result from ischemic heart disease and stroke, to which countries with low and medium incomes and high outdoor PM2·5 levels are primarily subject .
Most air pollution studies enjoying a large cohort did not make appropriate adjustments for individual-level risk factors (e.g., Social Economic Status (SES), health-related behaviors, or healthcare access) to the extent that they failed to properly appreciate the associations between PM2·5 and CVD.
The aging population and fast-paced urbanization are increasingly adding fuel to the prevalence of cardiovascular diseases in Iran.
Statistics retrieved from the Iranian Journal of Cardiovascular Nursing point to the cardiovascular diseases (incident cases) as the most frequent cause of disability in most countries including Iran in 2013 . To be specific, the total number of cardiovascular diseases in the present study was 441, with 351 cases (79.6%) in urban and 90 (20.4%) in rural areas and it followed a steadily rising trend.
A precise assessment of PM2.5 exposure could be quite beneficial to medical studies with the help of satellite data. Satellite observations measure the Earth’s optical radiations affected by aerosols. There are physical and acceptable statistical models for calculating PM2.5 density via satellite observations .
Backed by the datasets obtained from high-quality satellite-based PM2.5 estimates with a 1 × 1 km spatial resolution as well as the prediction records of Cardiovascular Disease Risk in Iran (Isfahan Cohort Study (ICS) - an established prospective project), the present study attempts to discover the potential relationship between long-term PM2.5 exposure and CVD-induced mortality and morbidity with respect to total CVD and its subtypes. This is reportedly the first study that has investigated the longitudinal effect of PM2.5 on cardiovascular diseases in Iran. In addition, it investigates the possibility that PM2.5 greater than a specific value (PM2.5 > 56.61) be closely associated with increased risks of CVD-induced morbidity and mortality. Subgroup analysis was carried out to study the PM2.5-induced morbidity and mortality under cardiovascular subtypes.
It is hoped that the findings of the present study would enhance conventional understandings about the global impacts of PM2·5-induced air pollution on CVDs given that they capture different population communities subject to a wide range of PM2·5 concentrations on a global scale and include standardized objective measures for CVD risk factors.
Materials and methods
Isfahan Cohort Study (ICS) was employed to investigate CVD incidence and the corresponding risk factors involved in the case of Iranian population. ICS is a population-based, longitudinal study that targets 6504 adults aging equal to or greater than 35 years at baseline, living in urban and rural areas from three counties in central Iran (Isfahan, Arak, and Najafabad) who had participated in the baseline survey of a community trial for CVD prevention and control, entitled Isfahan Healthy Heart Program (IHHP). They were recruited from January 2 to September 28, 2001. The baseline survey of the IHHP was conducted in a representative population of adults aging ≥19 years who were living in urban and rural areas of Isfahan, Arak, and Najafabad. ICS comprised 3081 adults living in 37 clusters. The participant selection method was multi-stage random sampling. The study population was first classified in terms of participants’ area of residence (urban or rural), and the regional population distribution was obtained from the national population census conducted in 1999. Census blocks were randomly selected from each county with a probability of selection proportional to the expected number of households, which were divided into approximately 1000 households. In each cluster, approximately 5-10% of the households were randomly selected for enumeration. One eligible 19-year-old subject was randomly selected from each house, provided they were mentally stable, unpregnant, and Iranian in nationality. The response rate for home interviews was 98%. However, only 95% attended the examination clinic. Because individual cardiovascular health was not included in the IHHP sampling criteria, participants initially included 108 cases (2.8) with a history of MI, stroke, or heart failure who were excluded from the ICS survey. IHHP samples were derived among different age and sex groups to reflect the fair age-sex distribution of the community. The prevalence of cardiovascular risk factors in the control area was estimated to be 0.2, and the sample size for IHHP had 90% strength to detect a relative risk of 0.75 at a significance level of 0.05. Arranged in the form of clustering, the total sample size in each region comprised 4828 people. To perform a series of follow-ups concerning the losses, we managed to employ 6300 people from each region (12,600 people). A total of 12,514 people were included in the initial survey. Individual selection for ICS ethical certification was obtained from the ethics committee of the Isfahan Cardiovascular Research Center, i.e., a cooperative center at World Health Organization. Of the 12,514 individuals reviewed in the baseline survey, the total number of 6640 adults aging ≥35 years were enrolled into the ICS . A full and detailed description of the study design was already given in .
Given that our study solely focuses on Isfahan province, the samples from Isfahan and Najafabad were included in the analysis, only. A total of 3081 adults from both rural and urban areas with complete information were initially enrolled for baseline examinations in 2001.
Multi-stage cluster random sampling was employed to select participants in terms of gender, age, and residential status (i.e., urban or rural). By using a formula for calculating the sample size, 3081 of them were selected in 2001.
Following an initial review of samples over the years 2001 to 2016, the follow-up procedure (telephone interview) for the candidates was performed every 2 years. Initial structured interviews were conducted based on a questionnaire with three main questions. “Is he/she alive?”, “has he/she been hospitalized for any reason? (With special focus on cardiovascular and cerebrovascular events), and “did the participant experience any of the following five neurological symptoms (hemiparesis, dysarthria, facial asymmetry, imbalance, and transient intermittent blindness)?” In the event of death, hospitalization, or neurological symptoms, the date of the event and the physician’s diagnosis must be taken into account. In case of any event, the questionnaire was reviewed along with relevant health records. In the case of out-of-hospital deaths, death certificates from the provincial mortality database were obtained and a verbal autopsy was performed by a specialist nurse over the course of a secondary interview in the presence of the surviving family members. A pre-defined questionnaire composed of such matters as medical history, signs, and symptoms of the patient prior to his/her death was made use of in the verbal autopsy .
The rate of ICS loss in the follow-up process was 6.4 and 3.9% in the second and third phases, respectively. Moreover, Arak province is situated in the central region of northwest of Isfahan and due to such problems as distance and difficulty in taking samples in rural areas, it was removed from the second phase of the ICS. Then, only Isfahan province was considered in this study to achieve complete follow-ups. Details of the loss to follow-up information and other issues are accessible in the following references [16, 17].
Measurements of the confounding factors
Having signed the informed written consent forms, the participants under study were interviewed by trained health professionals for 30 min. A valid questionnaire consisting of questions about demographic characteristics, socioeconomic statuses, behaviors, attitudes, skills, and knowledge about chronic, non-communicable, and other related diseases was employed. The participants’ lifestyle behaviors (smoking, physical activities, and eating habits), measurement of blood pressure, and anthropometric parameters were taken into consideration in compliance with standard protocols [18, 19] and calibrated instruments were used in this regard . In terms of smoking status, subjects were categorized as smokers (at least one cigarette per day or at least one cigarette per day in the past) and non-smokers . Participants’ fasting blood samples (10 ml) were obtained so as to gauge their total cholesterol (mg/dl), HDL-cholesterol (HDL-C), LDL-cholesterol (LDL-C), and triglycerides (mg/dl). The ratio of LDL to HDL was divided into three categories: normal (< 5.4), borderline (5.4-7), and abnormal (> 7) . Diabetes mellitus is defined as a condition with fasting blood glucose ≥126 mg/dl or when the patient is under anti-diabetic treatment. Impaired glucose tolerance is identified if the 2-h glucose levels are equal to and above 140 and below 200 mg/dl .
Obesity, overweight, and normal weight constitute the Body Mass Index (BMI): BMI ≥30 kg/m2, 25 ≤ BMI < 29.9 kg/m2, and BMI < 24.9 kg/m2 .
Hypertension was defined as a condition with blood pressure ≥ 140/90 mmHg or when the patients were under anti-hypertensive treatment .
A standard Food Frequency Questionnaire (FFQ) was employed to determine the participants’ dietary behaviors with appropriate validity and reliability .
Food Frequency Questionnaire includes dishes that are commonly consumed by Iranians. For each food item, participants were asked to report the amount and frequency of consumption in the past year. The total diet score is calculated as the sum of the scores given for each of the 12 food groups. Those subjects with a diet consisting of cereals, hydrogenated vegetable oils, red meat, and processed meat are considered healthy. The unhealthy group consumed sweets and pizza as well. Therefore, the total diet score for each person could range between 0 and 12. Participants with a total diet score of 8-10 are classified as healthy, while people with a diet score below 8 are classified as poor .
Physical activity was measured according to the International Physical Activity Questionnaire (IPAQ), the validity and reliability of which had already been reviewed and reported [26,27,28], and it was divided into three categories: low (less than 600 Met-min/Week), medium (at least 600 Met-min/Week), and high (at least 3000 Met-min/Week) .
This research was endorsed by the Isfahan University of Medical Sciences and approved by the Health Sciences Research Ethics Board, the code number for which is: 399099. The written consent forms were signed by each participant prior to data collection.
Outdoor PM2·5 air pollution exposure assessment
A satellite-based spatiotemporal model with a spatial resolution of 1 × 1 km was utilized to measure ambient PM2.5 levels. PM2.5 modeling refers to Land Use Regression (LUR) models that have recently drawn much attention because of their critical role in understanding and assessing air pollution [14, 30].
MODIS (Moderate Resolution Imaging Spectroradiometer) aerosol products, namely MOD04, and MAIAC (Multi-Angle Implementation of Atmospheric Correction) are multispectral sensors onboard two Terra and Aqua satellites and they were utilized in the present study. The satellites overpass the study area at approximately 10:30 and 13:30 local times every day . The statistical model proposed by Kong et al.  was employed to convert the obtained statistical observations to PM2.5. Thus, the seasonal localized regression coefficients were recalculated for our area based on the model presented by Kong et al.  and suggestions from Li et al. .
The data associated with ground stations in the city of Isfahan were taken from the Department of Environment for the reference years 2013 and 2018. Cloudy days, as well as days with wind speeds of more than 10 m/s, were put aside in compliance with the reports of the Iran Meteorological Organization for Isfahan. On such days, instantaneous satellite measurements significantly differ from ground-based measurements on a daily average, an issue that causes uncertainty.
The mean, standard deviation, and maximum PM2.5 concentration were calculated, and outlier data were deleted at a confidence level of 99%. In both of our references (Kong et al.  and Li et al. ), the seasonal division of data yielded better results than annual calculation. The seasons were divided according to those recommendations and based on the local climate (see Additional file 1).
In the method proposed by Kong et al. , calculating the average AOD in an area of 50 × 50 km was recommended so as to eliminate the effects of Land uses, although this averaging would make the PM2.5 estimations coarse in terms of spatial resolution. However, Li et al.  implemented the LUR model to control Land-use effects based on the variables such as population density, distance to industrial areas, and the area slope. Thus, they keep the spatial resolution, especially when considering the area within a city, as in our case.
Then, we use seasonal linear regression functions between AOD and PM2.5 which can be represented as PM2.5 = a * AOD + b, where a and b are regression coefficients. Regression coefficients are rounded to significant digits according to their estimation error (see Additional file 2).
A three-year rolling average of PM2·5 estimates was available in 2001, 2007, and 2013 with an approximate resolution of 1 × 1 km. Patients were enrolled in 2001 and their information was recorded and followed up until the end of 2015. PM2·5 levels in the participants’ residences were measured in 2001, 2007, and 2013, with an approximate resolution of 1 × 1 km. The values for the days of each year and subsequently, those for the 3 years were averaged. The mean of PM2.5 values for the years before the event was calculated for each participant, i.e., if individuals survived before 2007, only the amount of PM2.5 in 2001 was calculated; the same holds for the following: if they survived before 2013, the average PM2.5 for the years 2001 and 2007 was assessed; and if they survived before 2015, the average PM2.5 for 2001, 2007, and 2013 was measured. Indeed, the estimated exposure for the study participants was the mean concentration of PM2.5 over their surviving years during the cohort period. The obtained result was subsequently inputted in the model as a three-year rolling average and descriptive information about the three-year rolling average of PM2·5 variable was eventually provided.
The reported cases of Acute Myocardial Infarction (AMI) and stroke were monitored by the Isfahan Cardiovascular Research Center. In the absence of admission data in the registration database, well-trained nurses reviewed hospital-related medical records. The diagnosis of acute MI was based on the presence of at least two of the following criteria: (a) typical chest pain lasting more than 30 min with ST elevation> 0.1 mV in at least two adjacent electro cardiograph leads and (b) increased serum levels of cardiac biomarkers [34, 35]. Heart death was defined as death within 1 h of disease onset, cardiac arrest, or sudden collapse exhibiting no symptoms 41 h in advance. Ischemic Heart Disease (IHD) includes definite or probable MI, Unstable Angina (UA), and sudden cardiac death. In addition, the definition of stroke employed in this study to help diagnose stroke is in line with the definition presented by the World Health Organization. CVD was perceived as a combination of IHD and stroke. Although the doctors’ diagnoses were taken into account in the hospital, final decisions were made independently by the board . Total CVD, fatal CVD, non-fatal CVD, all-cause mortality, stroke, AMI, and IHD were separately included in different models.
Description of the baseline characteristics of the participants was issued based on means, standard deviation for continuous measurements, and percentages of categorical variables. The present study used the frailty model as an extended form of multivariate-adjusted Cox proportional hazards to estimate Hazard Ratios (HRs) and Confidence Interval of 95% (CIs) in the case of incidence of total CVD and its subtypes per 10 μg/m3 increase in long-term PM2.5 concentrations. The residential areas were considered as a random effect in models (according to the census blocks). The census blocks were randomly selected from each county with a probability of selection proportional to the expected number of households, which were divided into approximately 1000 households. Nearly 5–10% of the households were randomly selected for enumeration in each cluster.
Three models were applied to explore the relationship between long-term exposure to PM2.5 and CVD incidence or all-cause mortality including progressively adjusted covariates that were measured at baseline.
Model 1 incorporates age, gender, geographic area (urban or rural), and residential areas as random effects.
Model 2 is inclusive of individual risk factors namely smoking status, physical activity, dietary behavior, obesity, Social Economic Status (SES), and hypertension.
Model 3 is further adjusted in terms of diabetes, total cholesterol, triglycerides, LDL-to-HDL ratio, fatality, and history of heart disease in the family. This model can also be run by residence clustering as random effect.
To investigate individual covariates in association with total CVD, fatal CVD and stroke or AMI as two important subcategories of CVD, the two following hierarchical cox models were implemented.
As a basic model, Model 1 includes each covariate. Model 2 is adjusted by all the remaining covariates and residence clusters as random effect. The p-values presented in these models are the interaction results between PM2.5 exposure and individual/clinical variables.
The R 4.0.3 software was used to conduct the required analysis.
Over a 15-year course of following up with 3081 patients, 241 deaths occurred among which 105 death cases resulted from total CVD (fatal AMI, fatal stroke, and SCD). Also, there were 336 non-fatal CVD events (non-fatal incidence of AMIs and strokes). The mean (SD) three-year PM2.5 concentration in the study was 45.28 μg/m3 (11.63), ranging from a minimum of 20.01 μg/m3 in Najafabad to a maximum of 69.80 μg/m3 in Isfahan. There was a substantial variation in some individual/clinical variables by PM2·5 levels (Table 1). For example, 2429 (78.8%) of the study participants were urban dwellers, among whom 983 (95.8%) were exposed to the lowest PM2·5 tertiles, 825 (80.4%) to the middle, and 621(60.3%) to the highest PM2·5. In addition, the difference between ever smokers and cholesterol variables in different tertiles was significant, while the difference between the other variables in different tertiles was insignificant.
Model results for all-cause mortality and incidence related to CVDs are presented in Table 2. In Model 1, 10 μg/m3 increase in PM2.5 was significantly associated with 1.024 HR (0.95 CI = 1.015, 1.033) for CVD, which increased to 1.025 (0.95 CI = 1.016, 1.034) in Model 2. When the other variables were added (Model 3), the HR value increased to 1.026 (0.95 CI = 1.016, 1.036) for CVD events. In the fully adjusted models (Model 3), the most significant association was observed for non-fatal CVD with an HR of 1.029 (0.95 CI = 1.015, 1.034), for AMIs with an HR of 1.031 (0.95 CI = 1.005, 1.057), and for IHDs with an HR of 1.028 (0.95 CI = 1.017, 1.039). No risk association was observed between PM2.5 and fatal CVD.
The relationships of PM2.5 with fatal CVD, CVD, AMI, and stroke events according to the individual/clinical variables are summarized in Tables 3 and 4. The results illustrate that the effect of heterogeneity (or interaction with PM2.5) concerning different ages, gender, smoking status, high blood pressure, and diabetes status was statistically significant for different events.
For example, per 10 μg/m3 increase in PM2·5 for those subjects in the age of > 60 years, the HR was 1.046 for fatal CVD, 1.017 for CVD events, and 1.031 for stroke events. Moreover, the male gender determined HR of 1.006 for CVD events. Ever-smokers were subject to HR of 1.006 for CVD events and those with diabetes were subject to HR of 1.033 for fatal CVD, 1.016 for CVD events, 1.025 for AMI events, and 1.024 for stroke events (P for heterogeneity < 0·05). Overall, the results in Table 3 suggest that the effect of PM2.5 on cardiovascular disease was stronger among older people, smokers, and those with high blood pressure and diabetes.
According to the tertiles models, increasing the risk of CVD and stroke events was presented only for PM2.5 > 51.61 μg/m3.
Figure 1 shows the cumulative risk ratio and its associated confidence interval for different outcomes, based on which an upward trend could be perceived.
In Isfahan province, long-term outdoor PM2·5 exposure results from major CVD events. The findings of the present study offer new information on the relationship between ambient PM2·5 exposure and both CVD incidence and mortality in a wide range of PM2.5 concentrations (20.01 to 69.80 μg/m3) in different urban/rural areas while adjusting a comprehensive set of individual confounders.
Particulate Matter is of particular importance due to its specific characteristics such as composition and size distribution. This pollutant is characterized by a very large surface area and it could be adsorbed by many diverse organic materials such as polycyclic aromatic hydrocarbons, nitro-polycyclic aromatic hydrocarbons, heavy metals, pathogens, and radioactive materials. It contains fine particles penetrating the lower respiratory system and the blood and migrating to other organs, even to the brain .
According to the WHO recommendation, long-term exposure to extreme levels of PM2.5 over 10-25 μg/m3 could potentially impair the coagulation process, reduce inflammation, damage blood vessels, and eventually cause cardiovascular disease. Several pathways that could help justify the strong link between PM2.5 and cardiovascular diseases were identified [37, 38].
In Europe (Barcelona), PM2.5 concentration in the air of urban areas has been decreasing, and reducing air pollution in the metropolitan area of Barcelona would bring about substantial health and economic benefits .
In the present study, data sets obtained from environmental monitoring stations specific to a limited time period were provided and the correlation between satellite information and the data sets was found reasonable (R2 > 0.60).
A research study was conducted in Tehran to investigate the impact of temporal changes on the concentration of PM2.5 and a wide range of diseases. Reportedly, long-term exposure to PM2.5 caused deaths, estimated to be ranging between 24.5 and 36.2% out of total deaths induced by cerebrovascular diseases (stroke) and between 19.8 and 24.1% by IHD. They also found that IHD-induced deaths constituted the largest number of deaths resulting from long-term exposure to PM2.5. Overall, long-term exposure to PM2.5 significantly contributes to mortality in Tehran . Our results are consistent with those reported in this Iranian study, although we have found the significant effects only in the case of CVD incidence rather than CVD mortality .
Chen et al. (2008) conducted a systematic review of the relationship between long-term exposure to ambient pollution and chronic diseases. They found that such PM2·5 exposure would increase the risk of cardiovascular mortality by nearly 12–14% per 10-μg/m3 increase in PM2·5, independent of age, gender, and geographic region . Our results for CVD-related deaths turned out to be insignificant statistically due to the smaller sample size as compared to the results of different reviewed studies .
Inconsistency was found between our study and other recently conducted meta-analyses in which 53 studies on long-term PM2.5 and mortality in more than 150,000 adults with different SES levels were reviewed . It was demonstrated that in the case of all-cause mortality, a 1 μg/m3 increase in PM2.5 was significantly associated with a 1.29% increase in all-age and all-cause mortality (95%CI 1.09-1.50) at a mean exposure of 10 μg/m3. The percentage of increase in the rate of cardiovascular mortality was significant at 1.46% (95% CI, 1.25–1.67) with a mean exposure of 10 μg/m3.
Moreover, our study observed that the mean (SD) of 3-year PM2·5 concentrations was 45.28 (11.63) μg/m3 with the HR value of 1.026 (0.95 CI = 1.016, 1.036) per 10 μg/m3 increase in PM2·5 for all CVD events (see Table 2). Of note, Model 3 is a complete model that considers both confounders adjustment and residence clustering as a random effect. These estimates were considerably smaller than those cited in the meta-regression analyses due to the smaller sample size of our study in comparison to the result of different reviewed studies. However, some researchers have employed different modeling approaches (e.g., space-time exposure models) instead of land use regression method or fixed monitoring stations and considered the zip-codes of their study subjects controlling for SES levels; therefore, their approach estimated a larger number of deaths, unlike our approach which estimated a lower rate of mortality due, mainly, to the smaller sample size .
Outcome of the adjusted models achieved by Hystad et al.  demonstrated that during a median follow-up period of 9·3 years (between Jan 1, 2003 and July 14, 2018) targeting 157,436 adults belonging to 747 communities in 21 countries, a 10 μg/m3 increase in PM2·5 was associated with an increased risk for cardiovascular disease events (hazard ratio 1·05 [95% CI 1·03–1·07]), myocardial infarction (1·03 [1·00–1·05]), stroke (1·07 [1·04–1·10]), and cardiovascular disease mortality (1·03 [1·00–1·05]). They concluded that long-term outdoor PM2·5 concentrations were associated with increased risks of cardiovascular disease in adults aging 35–70 years. In terms of the risk of contracting cardiovascular disease, our research shares similarities to the above-mentioned research study.
In another study carried out by Abdolahnejad et al., Relative Risk (RR) pointed to the increased risk of exposure to pollutants, as obtained through time-series experiments. In Isfahan, upon every 10 μg/m3 increase in the concentration of the pollutant, the value of RR per increase in the total mortality induced by PM10 was 0.8% . Our results are inconsistent with those of Abdolahnejad et al. and this divergence results from different modeling methods and follow-up periods.
Hvidtfeldt UA et al. put a greater emphasis on exposure data from central air monitoring stations and followed a Danish cohort of 49,564 individuals from the years 1993–1997 to 2015. Among the 895,897 annual follow-up cases, 10,193 deaths induced by all causes occurred, of which 2319 were CVD related. They observed that the relationship between long-term exposure to PM2.5 and CVD mortality was significant (hazard ratio of 1.03; 95% confidence interval: 1.01, 1.05) . Our results were not significant because of the longer follow-up period and the larger number of participants .
Zhang et al. investigated a retrospective cohort consisting of 39,054 subjects from four cities in northern China for the mortality induced by all-cause and specific cardiovascular diseases over the course of the years 1998 to 2009 in order to assess the cardiovascular effects of long-term exposure to high-level concentrations of inhalable particulate. Their analyses suggested that the effects of PM10 on cardiovascular mortality were more pronounced in males, smokers, and people with a higher socioeconomic status. Long-term exposure to PM10 increases mortality from cardiovascular disease, especially from ischemic heart disease ; this finding is consistent with our presentation.
Chen et al. performed a random-effects meta-analysis. This study included cohort and case-control studies on outdoor air pollution in human populations using individual-level data. In addition to natural-cause mortality, they evaluated mortality from circulatory diseases (IHD and cerebrovascular disease (stroke)). They concluded that there was enough evidence to suggest that both PM2.5 and PM10 were associated with increased mortality from all causes and cardiovascular disease . This result is not in agreement with that of our study because the aforementioned research study carried out different types of outdoor exposure and checked PM2.5 and PM10 simultaneously .
In the present paper, following each 10 μg/m3 increase in PM2.5, HR for fatal CVD was measured at 0.998 (0.95 CI = 0.979, 1.017). Therefore, no evidence of risk associations was observed between PM2.5 increase and fatal CVD. Findings acquired from Tseng et al. confirmed the results of the present study. A cohort study was carried out on the potential associations between PM2.5 and cardiovascular mortality with 43,227 patients in Taiwan. Participants were followed up from the year 1989 to 2008 and their vital status was determined based on death records. In this study, no evidence of an increased risk of all-cause or cardiovascular mortality with a prolonged exposure to PM2.5 was detected , similar to our study.
Regional characteristics of Isfahan are the reasons for a higher percentage of PM2.5 concentration in nearby rural areas. The occurrence of this phenomenon could be attributed to the proximity of these areas to deserts and swamps. Deserts are severely affected by wind erosion and due to the prevailing wind direction in these areas blowing from the east to the west, many suspended particles are transported to the western side, inevitably causing air pollution. The Gavkhouni Wetland, i.e., the endpoint of Zayanderud (a famous river in Isfahan) which is a place to trap suspended particles, has run dry due to reduced rainfall and improper farming and irrigation methods and has turned into a source of many dust particles in the studied period.
Furthermore, sand mines and brick kilns situated in the east and northeast of Isfahan are essential sources of the production of dust particles and their spread to the air of the city. Similarly, anomalies could be observed as the average concentration of suspended particles in different city areas. Such anomalies could be caused by various factors such as industrial activities in the suburbs, vehicle and urban traffic, proximity of the area to the landfill, reduced vegetation, and green spaces [52, 53].
The main objective of the study of Khusfi et al. was to investigate the effects of atmospheric conditions and physical characteristics of the Earth on the seasonal variations in dust concentration in the semi-arid regions of Central Iran Zone (CIZ) in a typical period (2001–2008 and 2009–2016). It was found that the activity of Sand and Dust Storms (SDS) in 2009-2016 was greater than that in the previous period (2001-2008), i.e., a finding consistent with our study results .
Indeed, according to the results presented in Table 1, 78.8% of the cohort population inhabited the urban areas. However, in the third tertile, this number was 60.3% (urban areas in the third tertile of PM2.5 decreased in number and more rural areas were included), and the difference between the two mentioned values was significant [52, 53].
Moreover, in rural areas, villagers are exposed to smoke from burning agricultural waste. These results were similarly obtained using a cross-sectional study performed through random selection of three villages in Isfahan in order to evaluate the respiratory effects of this smoke-induced air pollution. It was shown that in symptomatic cases before and after rice burning, the total particulate matters doubled during the burning time. In this study, all clinical and spirometry changes were statistically significant between study subjects .
Based on our results, only CVD events and stroke, as its essential subtype, had a significantly strong dose-response relationship; HRs related to the 3rd tertile of PM2·5 were1.652 (0.95 CI = 1.269, 2.150) for total CVD events and 1.791 (0.95 CI = 1.036, 3.097) for stroke events. Taken together, our study results support the evidence of a significant increase in total CVD risk, especially for stroke, at high PM2·5 concentrations.
In addition to the diverse geographical population and extensive PM concentration range, the ICS is unique in its depth of individual variables to adjust the potential confounding factors.
Besides the long-term follow-up (15 years), ICS considered a range of sociodemographic, behavioral, metabolic, and clinical variables that affect CVD, an issue that most large cohort studies examining air pollution did not measure .
From the results, one could conclude that long-term outdoor PM2·5 is a significant contributing factor in CVD in Iran. While it is unlikely that outdoor PM2·5 concentrations would vary substantially over such small areas, some exposure misclassifications exist that could bias estimates towards the null.
We compared PM2.5 concentrations and CVD events across different areas. While residual confounding cannot be ruled out, we adjusted more individual CVD risk factors than previous studies; the analyses controlling for unmeasured factors among the centers using random-effects demonstrated larger estimates of the effect. We could not examine specific causes of non-fatal CVD due to the smaller number of events; however, these analyses can be done in the future with additional follow-ups and more events.
The present research enjoys some merits including being a longitudinal study with a 15-year follow-up, sufficient number of samples, inclusive outcomes, sufficient number of urban and rural population exposed to high PM2.5 concentrations for the whole study period, objective measurement of a comprehensive suite of individual CVD risk factors, standardized collection of data on household and community characteristics, and prospective recording of fatal and non-fatal events that were evaluated through standard definitions [15, 34, 35].
The present study was subject to the following limitations. First, there was an insufficient scope data on indoor air pollution sources such as household use of solid fuels. Second, exposure to ambient PM2.5 was estimated based on the participants’ residential addresses regardless of their daily activities and location changes during days. Third, acute (i.e., daily) variations in PM2.5 exposures and their impacts could not be controlled. Fourth, we did not consider local environmental monitoring data because our study was a 15-year cohort. Existing local environmental monitoring data for these 15 years were not entirely available; local environmental monitoring centers were limited; and the distance from the patients’ residential neighborhoods was far away from the stations. However, satellite data with a 1 × 1 km spatial resolution included the patients’ residential addresses which were more accurate while local environmental monitoring data were not satisfactory in terms of quality.
The present study demonstrated that long-term ambient PM2.5 exposure was associated with the increased risk of contracting CVDs including AMI, IHD, and strokes. Our data support the assumption that PM2.5 is an important risk factor in CVD.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available upon request to the corresponding author.
- PM2.5 :
Particles with a mass median aerodynamic diameter ≤ 2.5 μm
Acute Myocardial Infarction
Ischemic heart disease
Coronary heart disease
Sudden cardiac death
Unstable angina pectoris
Social Economic Status
Di Q, Wang Y, Zanobetti A, Wang Y, Koutrakis P, Choirat C, et al. Air pollution and mortality in the Medicare population. N Engl J Med. 2017;376(26):2513–22.
Rhinehart ZJ, Kinnee E, Essien UR, Saul M, Guhl E, Clougherty JE, et al. Association of fine particulate matter and risk of stroke in patients with atrial fibrillation. JAMA Netw Open. 2020;3(9):e2011760.
USEPA. National ambient air quality standards for particulate matter. Final rule. Fed Reg. 2013;78:3086–287 40 CFR Parts 50, 51, 52, 53, and 58.
European Environment Agency. Air quality in Europe - 2017 report. 2017. http://www.eea.europa.eu//publications/air-quality-in-europe-2017.
Yim SH, Gu Y, Shapiro MA, Stephens B. Air quality and acid deposition impacts of local emissions and transboundary air pollution in Japan and South Korea. Atmospheric Chemistry and Physics. 2019;19(20):13309–23.
Tanaka TY, Chiba M. A numerical study of the contributions of dust source regions to the global dust budget. Glob Planet Chang. 2006;52(1-4):88–104.
Chang SC, Lee CT. Assessment of PM10 enhancement by yellow sand on the air quality of Taipei, Taiwan in 2001. Environ Monit Assess. 2007;132(1-3):297–309.
Modarres R, Dehkordi AK. Daily air pollution time series analysis of Isfahan City. Int J Environ Sci Technol. 2005;2(3):259–67.
Talebi SM, Tavakoli T, Ghinani A. Levels of PM10 and its chemical composition in the atmosphere of the city of Isfahan. Iran J Chem Eng. 2008;5(3):62–7.
World Health Organization. Ambient air pollution: World Health Organization website; 2019. https://www.who.int/airpollution/ambient/en/. Accessed 10 Apr 2019.
Zhou M, Wang H, Zeng X, Yin P, Zhu J, Chen W, et al. Mortality, morbidity, and risk factors in China and its provinces, 1990–2017: a systematic analysis for the global burden of disease study 2017. Lancet. 2019;394(10204):1145–58.
Hystad P, Larkin A, Rangarajan S, AlHabib KF, Avezum Á, Calik KB, et al. Associations of outdoor fine particulate air pollution and cardiovascular disease in 157 436 individuals from 21 high-income, middle-income, and low-income countries (PURE): a prospective cohort study. Lancet Planet Health. 2020;4(6):e235–45.
Shahsavari S, Nazari F, Karimyar Jahromi M, Sadeghi M. Epidemiologic study of hospitalized cardiovascular patients in Jahrom hospitals in 2012-2013. Iran J Cardiovasc Nurs. 2013;2(2):14–21.
Gulliver J, Morley D, Dunster C, McCrea A, van Nunen E, Tsai M-Y, et al. Land use regression models for the oxidative potential of fine particles (PM2.5) in five European areas. Environ Res. 2018;160:247–55.
Sarrafzadegan N, Talaei M, Sadeghi M, Kelishadi R, Oveisgharan S, Mohammadifard N, et al. The Isfahan cohort study: rationale, methods and main findings. J Hum Hypertens. 2011;25(9):545–53.
Sarrafzadegan N, Hassannejad R, Roohafza H, Sadeghi M, Talaei M, Oveisgharan S, et al. A 10-year Isfahan cohort on cardiovascular disease as a master plan for a multi-generation non-communicable disease longitudinal study: methodology and challenges. J Hum Hypertens. 2019;33(11):807–16.
Talaei M, Sarrafzadegan N, Sadeghi M, Oveisgharan S, Marshall T, Thomas GN, et al. Incidence of cardiovascular diseases in an Iranian population: the Isfahan Cohort Study. Arch Iran Med. 2013;16(3):138–44.
North American Association for the Study of Obesity, National Heart, Lung, Blood Institute, NHLBI Obesity Education Initiative. The practical guide: identification, evaluation, and treatment of overweight and obesity in adults. National Institutes of Health, National Heart, Lung, and Blood Institute, NHLBI Obesity Education Initiative, North American Association for the Study of Obesity; 2000.
Chobanian AV. National heart, lung, and blood institute joint national committee on prevention, detection, evaluation, and treatment of high blood pressure; national high blood pressure education program coordinating committee: the seventh report of the joint national committee on prevention, detection, evaluation, and treatment of high blood pressure: the JNC 7 report. JAMA. 2003;289:2560–72.
National Cholesterol Education Program (US). Expert Panel on Detection, Treatment of High Blood Cholesterol in Adults. Third report of the National Cholesterol Education Program (NCEP) Expert Panel on detection, evaluation, and treatment of high blood cholesterol in adults (Adult Treatment Panel III). The Program; 2002.
World Health Organization. WHO draft protocol and manual of operations population survey for cardiovascular disease risk factors in the Eastern Mediterranean Region. Alexandria: World Health Organization; 1995. p. 35.
World Health Organization. Obesity: preventing and managing the global epidemic: report of a WHO consultation on obesity, Geneva, 3-5 June 1997. World Health Organization; 1998.
Working Group on Risk and High Blood Pressure. An epidemiological approach to describing risk associated with blood pressure levels. Hypertension. 1985;7:641–51.
Mohammadifard N, Sajjadi F, Maghroun M, Alikhasi H, Nilforoushzadeh F, Sarrafzadegan N. Validation of a simplified food frequency questionnaire for the assessment of dietary habits in Iranian adults: Isfahan Healthy Heart Program, Iran. ARYA Atheroscler. 2015;11(2):139.
Mohammadifard N, Talaei M, Sadeghi M, Oveisegharan S, Golshahi J, Esmaillzadeh A, et al. Dietary behaviorss and mortality from cardiovascular disease: Isfahan Cohort Study. Eur J Clin Nutr. 2017;71(2):252–8.
Vasheghani-Farahani A, Tahmasbi M, Asheri H, Ashraf H, Nedjat S, Kordi R. The Persian, last 7-day, long form of the international physical activity questionnaire: translation and validation study. Asian J Sports Med. 2011;2(2):106.
Moghaddam MB, Aghdam FB, Jafarabadi MA, Allahverdipour H, Nikookheslat SD, Safarpour S. The Iranian Version of International Physical Activity Questionnaire (IPAQ) in Iran: content and construct validity, factor structure, internal consistency and stability. World Appl Sci J. 2012;18(8):1073–80.
Klishadi R, Khosravi A, Famouri F, Sadeghi M, Shirani S. Assessment of physical activity of adolescents in Isfahan. J Shahrekord Univ Med Sci. 2001;3(2):55–66.
Talaei M, Rabiei K, Talaei Z, Amiri N, Zolfaghari B, Kabiri P, et al. Physical activity, sex, and socioeconomic status: a population based study. ARYA Atheroscler. 2013;9(1):51.
Wu C-D, Chen Y-C, Pan W-C, Zeng Y-T, Chen M-J, Guo YL, et al. Land-use regression with long-term satellite-based greenness index and culture-specific sources to model PM2.5 spatial-temporal variability. Environ Pollut. 2017;224:148–57.
Nordio F, Kloog I, Coull BA, Chudnovsky A, Grillo P, Bertazzi PA, et al. Estimating spatio-temporal resolved PM10 aerosol mass concentrations using MODIS satellite data and land use regression over Lombardy, Italy. Atmos Environ. 2013;74:227–36.
Kong L, Xin J, Zhang W, Wang Y. The empirical correlations between PM2. 5, PM10 and AOD in the Beijing metropolitan region and the PM2. 5, PM10 distributions retrieved by MODIS. Environ Pollut. 2016;216:350–60.
Li R, Ma T, Xu Q, Song X. Using MAIAC AOD to verify the PM2. 5 spatial patterns of a land use regression model. Environ Pollut. 2018;243:501–9.
Luepker RV, Apple FS, Christenson RH, Crow RS, Fortmann SP, Goff D, et al. Case definitions for acute coronary heart disease in epidemiology and clinical research studies: a statement from the AHA Council on Epidemiology and Prevention; AHA Statistics Committee; World Heart Federation Council on Epidemiology and Prevention; the European Society of Cardiology Working Group on Epidemiology and Prevention; Centers for Disease Control and Prevention; and the National Heart, Lung, and Blood Institute. Circulation. 2003;108(20):2543–9.
Committee Members, Braunwald E, Antman EM, Beasley JW, Califf RM, Cheitlin MD, et al. ACC/AHA guideline update for the management of patients with unstable angina and non–ST-segment elevation myocardial infarction—2002: summary article: a report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines (Committee on the Management of Patients With Unstable Angina). Circulation. 2002;106(14):1893–900.
Mokhtari M, Miri M, Mohammadi A, Khorsandi H, Hajizadeh Y, Abdolahnejad A. Assessment of air quality index and health impact of PM10, PM2. 5 and SO2 in Yazd, Iran. J Maz Univ Med Sci. 2015;25(131):14–23.
WHO. World-Health-Organization. WHO Air Quality Guidelines for particular matter, ozone, dioxide and sulfur dioxide. Global update 2005. Summary of risk assessment. WHO Publications, Germany. 2005;74:75.
Doyle K. Pollution particles damage blood vessels, may lead to heart disease. New York: Reuters; 2016.
Pérez L, Sunyer J, Künzli N. Estimating the health and economic benefits associated with reducing air pollution in the Barcelona metropolitan area (Spain). Gac Sanit. 2009;23(4):287–94.
McGuinn LA, Ward-Caviness CK, Neas LM, Schneider A, Diaz-Sanchez D, Cascio WE, et al. Association between satellite-based estimates of long-term PM2.5 exposure and coronary artery disease. Environ Res. 2016;145:9–17.
Liang F, Liu F, Huang K, Yang X, Li J, Xiao Q, et al. Long-term exposure to fine particulate matter and cardiovascular disease in China. J Am Coll Cardiol. 2020;75(7):707–17.
Madrigano J, Kloog I, Goldberg R, Coull BA, Mittleman MA, Schwartz J. Long-term exposure to PM2. 5 and incidence of acute myocardial infarction. Environ Health Perspect. 2013;121(2):192–6.
Sahu SK, Sharma S, Zhang H, Chejarla V, Guo H, Hu J, et al. Estimating ground level PM2. 5 concentrations and associated health risk in India using satellite based AOD and WRF predicted meteorological parameters. Chemosphere. 2020;255:126969.
Faridi S, Shamsipour M, Krzyzanowski M, Künzli N, Amini H, Azimi F, et al. Long-term trends and health impact of PM2.5 and O3 in Tehran, Iran, 2006–2015. Environ Int. 2018;114:37–49.
Chen H, Goldberg MS, Villeneuve PJ. A systematic review of the relation between long-term exposure to ambient air pollution and chronic diseases. Rev Environ Health. 2008;23(4):243–97.
Vodonos A, Awad YA, Schwartz J. The concentration-response between long-term PM2. 5 exposure and mortality; a meta-regression approach. Environ Res. 2018;166:677–89.
Abdolahnejad A, Jafari N, Mohammadi A, Miri M, Hajizadeh Y, Nikoonahad A. Cardiovascular, respiratory, and total mortality ascribed to PM10 and PM2. 5 exposure in Isfahan, Iran. J Educ Health Promot. 2017;6:109.
Hvidtfeldt UA, Geels C, Sørensen M, Ketzel M, Khan J, Tjønneland A, et al. Long-term residential exposure to PM2. 5 constituents and mortality in a Danish cohort. Environ Int. 2019;133:105268.
Zhang LW, Chen X, Xue XD, Sun M, Han B, Li CP, et al. Long-term exposure to high particulate matter pollution and cardiovascular mortality: a 12-year cohort study in four cities in northern China. Environ Int. 2014;62:41–7.
Chen J, Hoek G. Long-term exposure to PM and all-cause and cause-specific mortality: a systematic review and meta-analysis. Environ Int. 2020;143:105974.
Tseng E, Ho WC, Lin MH, Cheng TJ, Chen PC, Lin HH. Chronic exposure to particulate matter and risk of cardiovascular mortality: cohort study from Taiwan. BMC Public Health. 2015;15(1):1–9.
Mazloumi A, Alaminia Z, Fadaei M. Study of suspended solids concentration based on air quality index (AQI) in Isfahan. In: The third national conference on environmental and agricultural research in Iran, Hamadan; 2015.
Gauvin S, Reungoat P, Cassadou S, Dechenaux J, Momas I, Just J, et al. Contribution of indoor and outdoor environments to PM2. 5 personal exposure of children—VESTA study. Sci Total Environ. 2002;297(1-3):175–81.
Khusfi ZE, Khosroshahi M, Roustaei F, Mirakbari M. Spatial and seasonal variations of sand-dust events and their relation to atmospheric conditions and vegetation cover in semi-arid regions of central Iran. Geoderma. 2020;365:114225.
Golshan M, Faghihi M, Roushan-Zamir T, Marandi MM, Esteki B, Dadvand P, et al. Early effects of burning rice farm residues on respiratory symptoms of villagers in suburbs of Isfahan, Iran. Int J Environ Health Res. 2002;12(2):125–31.
This study was funded by Isfahan Cardiovascular Research Institute.
Ethics approval and consent to participate
Study protocols were reviewed and approved by the Research and Ethics Council of Isfahan University of Medical Sciences (code number 399099) and conformed to the Declaration of Helsinki.
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.
About this article
Cite this article
Jalali, S., Karbakhsh, M., Momeni, M. et al. Long-term exposure to PM2.5 and cardiovascular disease incidence and mortality in an Eastern Mediterranean country: findings based on a 15-year cohort study. Environ Health 20, 112 (2021). https://doi.org/10.1186/s12940-021-00797-w
- Outdoor air pollution
- Cardiovascular diseases
- Survival models
- Cox proportional hazards frailty models