Ambient air quality and the risk for Chronic Obstructive Pulmonary Disease among Metro Manila Development Authority traffic enforcers in Metro Manila: An exploratory study.

BACKGROUND: Air pollution and poor ambient air quality are significantly related to multiple health risks. One associated disease is chronic obstructive pulmonary disease (COPD), a preventable disease with several contributing factors and one of the leading causes of morbidity/mortality locally and globally. A potentially high-risk population are traffic enforcers who are constantly exposed to air pollution. In the Philippines, the MMDA has the widest coverage in traffic management. The study determined the risk of COPD among Metro Manila Development Authority (MMDA) traffic enforcers in relation to ambient air quality level, as well as identified other factors that increase the risk of developing COPD.
METHODS: Fifty-two MMDA traffic enforcers deployed in PM2.5 air quality sensor areas in Metro Manila from 2016 to 2018 were recruited through stratified sampling. The International Primary Airways Guidelines (IPAG) questionnaire was utilized to measure risk of COPD. Respiratory health and working history were obtained through questionnaires. Department of environment and natural resources provided PM2.5 ambient air quality data which aided in the construction of the Exposure-Month Index. Ordinal logistic regression was used to examine the association of PM2.5 together with the relevant factors and the risk of COPD.
RESULTS: We found statistically significant associations between PM2.5 and COPD among high risk category [odds risk (OR): 1.24, 95% confidence interval (CI): 1.07-1.44]. Age (Moderate, OR: 1.16, 95% CI: 0.98-1.38 and High, OR: 10.06, 95% CI: 4.02-25.17) and chest pain (Moderate, OR: 68.65, 95% CI: 1.71-2.75 × 103) were potential risk factors, whereas body mass index (BMI) (OR: 0.05, 95% CI: 0.01-0.53) exhibited protective effect.
CONCLUSIONS: Exposure to PM2.5 was associated with an increased risk of COPD among high-risk category MMDA traffic enforcers. Age and chest pain were potential risk factors to risk of COPD, whereas BMI exhibited a potential protective effect. Results of this study can be used for clinical management of high-risk populations, such that of MMDA traffic enforcers.

Introduction

Globally, air pollution has negatively affected air quality and has become a threat to public health. In fact, 4.2 million deaths annually were estimated to be due to ambient air pollution. Morbidity and mortality from cardiovascular and respiratory causes, in particular, have been observed to be associated with particulate air pollution., Among the various criteria air pollutants, particulate matter (PM) is found to be responsible for the majority of associated health risks., An increasing emphasis was thus placed on PM with a size less than 2.5 microns in diameter (PM2.5) because of its ability to reach distal regions of the lung.

One of the diseases associated with poor air quality is chronic obstructive pulmonary disease (COPD), which is a preventable yet commonly undiagnosed and life-threatening respiratory condition. It is characterized by persistent limitation of airflow and enhanced chronic inflammatory response in the airways to noxious particles or gases. Mainstay symptoms of COPD are chronic, progressive dyspnea, cough, sputum, or mucus production—these usually become apparent after 40–50 years of age.

The primary cause of COPD remains to be exposure to tobacco smoke. However, recent scientific developments revealed the substantial role of PM2.5 on the exacerbation of COPD.8, 9, 10, 11, 12, 13 Other modifiable risk factors include obesity and occupational exposures to organic and inorganic dusts, chemicals, fumes, and irritants. Meanwhile, non-modifiable factors include increasing age, male gender, lung development abnormalities, bronchial hyperreactivity, long-term asthma, childhood respiratory infections, and family history of respiratory diseases.,,

In the Philippines, the Department of Health recognizes COPD as the seventh leading cause of mortality. In 2007, an estimated 14% of Metro Manila's population alone reportedly had COPD. Several risk factors are present in the country—the most significant of which are air pollution (96%) in urban settings and biomass exposure (87%) in rural settings. A significant contributor to air pollution are motor vehicles—this is especially true in Metro Manila where there is a 2.5 million annual average daily traffic. In fact, the region's average PM2.5 level (30.44 μg/m3) from 2016 to 2018 lies beyond the recommended annual PM2.5 level (20 μg/m3) by the World Health Organization. With considerably high ambient PM2.5 levels, outdoor occupations are mostly hazardous in nature given the prolonged exposure.

Traffic enforcers had the highest exposure, who were found to have an increased risk of developing respiratory and cardiac diseases.22, 23, 24 In Metro Manila, the Metro Manila Development Authority (MMDA) remains to be the agency with the widest coverage in traffic management. On the job, MMDA traffic enforcers work on 8-hour shifts on weekdays. Given the nature and scope of the work, the MMDA may then adequately represent traffic enforcers with a high risk of exposure to high levels air pollution in Metro Manila.

In response to a growing need for evidence linking air quality to health outcomes in the country, specifically among high risk groups, this study explored the association between ambient air quality and risk of COPD among MMDA traffic enforcers in Metro Manila.

Methods

We utilized a cross-sectional study design in elucidating the risk of developing COPD due to ambient air quality. The list of variables with their corresponding units, definitions, type, and sources of information are summarized in Table S1.

Ethical approval

Ethical clearance was obtained from the ASMPH Ethical Research Board with Ethical Approval number: SMPH 2019 Group 02. Permission to conduct the study in the selected MMDA satellites was also coordinated with the main headquarters in Makati City and each MMDA satellite head prior to each visit. Prior to assessment, all participants were briefed and given informed written consent forms (see Supplementary Document 3), detailing the purpose and methods of the study and their rights as study participants. Gathered data were kept confidential throughout the entire study, in compliance with the Data Privacy Act of 2012.

Independent variable

For the independent variable, an Exposure-Month Index (EMI) which is calculated as the average PM2.5 exposure of the participant in the month/s where s/he was deployed. Since traffic enforcers are assigned to different areas in varying time periods, EMI captures the cumulative exposure of the traffic enforcers, which is dependent on the length of time the subject stayed in the location. The EMI was assessed with the aid of a working history form (see Supplementary Document 1) and was computed as follows:whereby, k represents the work location; PM2.5 represents the PM2.5 concentration in the work location; duration is the length of time the traffic officer was assigned in that location. The summation of the duration weighted PM2.5 is the exposure index () for traffic officer y.

PM 2.5 levels were obtained from air quality ground monitoring data from the department of environment and natural resources-air quality management section (DENR-AQMS) measured via BAM (Thermos, Met 1) and TEOM (Thermos) sensors from 2016 to 2018 (see Table S2 for PM 2.5 Yearly Average).

Dependent variable

For the dependent variable, risk of COPD was measured using the International Primary Airways Guidelines (IPAG) risk score, which was later transformed into a categorical variable with three levels: namely, low, moderate, and high risk. The IPAG Questionnaire, highlighted in Appendix 1 of Sichletidis, Spyratos, Papaioannou et al, is the most widely validated screening questionnaire aimed at determining risk of COPD (sensitivity = 64.5%; specificity = 65.2%). This tool includes questions regarding respiratory symptoms, smoking history, allergies, age, and anthropometrics. There are a total of 38 points, with a cut-off score of greater than 17 points as high risk of COPD. The questionnaire was translated to Tagalog for comprehensibility and incorporated in the Respiratory Health Questionnaire (see Supplementary Document 2).

Study sites

The study was conducted in nine selected MMDA satellite stations with an active PM2.5 air quality sensor installed and monitored by the DENR-AQMS. Fig. 1 illustrates the coverage area of PM2.5 sensors within a 5-kilometer radius and the location of the different MMDA satellites. Three-year PM2.5 levels in these sites were obtained from DENR-AQMS and processed with the assistance of environmental health experts at the ASMPH Center for Research and Innovation (ACRI), Environmental Health unit.

Fig. 1

Map of Metro Manila showing air quality sensors' area coverage (in circles) and selected MMDA satellites.

Study population

MMDA traffic enforcers satisfying any of the following criteria were excluded from the study population: (1) were not currently assigned to the field, (2) were assigned to the field for more than 6 months in an area that does not have a PM 2.5 air quality sensor, (3) were newly assigned to the field from January 2019 onwards, and (4) previously diagnosed with COPD before January 2016. Study sample consisted of 52 MMDA traffic enforcers (prevalence = 0.883%, power = 80, design effect = 1, confidence interval = 95%) assigned to field jobs for at least 6 months.

Data collection procedure

Pre-testing of questionnaires was done to three participants in the EDSA-MMDA satellite for the purpose of training the researchers and testing the data collection instruments. During actual data collection, a minimum of three researchers went per study site and requested for available participants who were currently not on duty. All data gathered were transcribed and encrypted to ensure confidentiality and, afterwards, will be permanently deleted once the study has concluded.

Statistical analyses

Data obtained was encoded in Excel version 2016 (Microsoft, USA). We utilized an ordinal logistic regression to assess association of the risk of COPD and PM2.5, after adjusting for several covariates. A bivariable analysis was utilized to examine the association of IPAG variable and the covariates. Covariates which exhibited statistical significance were then included to the full model in the multivariable regression. The full model was then subjected to a stepwise multivariable regression via backward elimination, whereby variables which did not show statistical significance in the full model were dropped one-by-one. The final model favored the inclusion of body mass index (BMI), age, symptoms of chest pain and PM2.5. The subsequent interpretations of the odds risk (OR) for moderate and high categories are relative to the low category, which is set as the base level of comparison in the ordinal logistic regression. All statistical analyses were done using R Statistical Programming (R Development Core Team, USA).

Results

In Table 1, we observed strong statistically significant association of the IPAG categories with age (High, OR: 1.48, 95% CI: 1.05–2.07), pack years (Moderate, OR:1.14, 95% CI: 0.03–50.94 and High, OR:1.20, 95% CI: 1.05–1.38), chest pain (Moderate, OR: 19.50, 95% CI: 6.97–54.57). Whereas, BMI (High, OR: 0.76, 95% CI: 0.57–1.02), years of service (High, OR: 1.11, 95% CI: 0.98–1.25), past respiratory disease (High, OR: 6.17, 95% CI: 0.78–48.64) and PM2.5 (High, OR: 1.06, 95% CI: 0.99–1.13) exhibited a relatively moderate statistical significance. All statistically significant variables exhibited a positive association with IPAG categories, except for BMI, which indicated a possible negative association.

Table 1

Bivariable, ordinal logistic regression results of International Primary Airways Guidelines categories with age.

Variables	Moderate		High
	OR (95% CI)	P values	OR (95% CI)	P values
BMIa	0.90 (0.00–915)	0.36	0.76 (0.57–1.02)	0.07b
Agea	1.09 (0.00–7.5 × 107)	0.14	1.48 (1.05–2.07)	0.02b
Current Smoker	0.31 (0.09–1.03)	0.31	1.04 (0.16–6.94)	0.97
Past Smoker	0.51 (0.07–4.01)	0.51	1.54 (0.12–19.5)	0.74
Years of Servicea	1.09 (0.13–9.09)	0.12	1.11 (0.98–1.25)	0.10b
Pack Yearsa	1.14 (0.03–50.94)	0.03b	1.20 (1.05–1.38)	0.01b
Difficulty of Breathing	1.55 (0.47–5.07)	0.64	2.07 (14.18–0.46)	0.46
Sneezing	0.87 (0.26–2.86)	0.88	1.16 (0.17–7.72)	0.88
Colds	1.55 (0.55–4.39)	0.64	0.78 (0.08–7.76)	0.83
Atopy	–	–	2.31 (0.21–26.05)	0.50
Cough	0.83 (0.29–2.33)	0.87	1.03 (0.10–10.53)	0.98
Chest Paina	19.50 (6.97–54.57)	0.01b	4.88 (0.36–66.41)	0.23
Shortness of Breath	0.39 (0.12–1.27)	0.40	1.29 (0.19–8.61)	0.80
Current Respiratory Disease	3.90 (1.39–10.91)	0.30	4.87 (0.36–66.41)	0.23
Past Respiratory Diseasea	4.62 (1.43–15.00)	0.13	6.17 (0.78–48.64)	0.08b
Symptoms in family	2.56 (0.35–18.92)	0.31	5.11 (0.52–49.80)	0.16
Family history of respiratory disease	1.35 × 104 (4.68 × 103–3.92 × 104)	0.82	0.32 (0.03–3.11)	0.33
PM2.5a	1.03 (0.00–2.93 × 103)	0.32	1.06 (0.99–1.13)	0.10b

BMI: Body Mass Index; OR: Odds Ratio; 95% CI: 95% Confidence Interval; PM2.5: particulate matter with a size of 2.5 microns or less in diameter.

Reference level of COPD risk is low.

Variables which were included in the full model of the subsequent multivariable, multinomial analysis.

P value ≤ 0.10 indicate statistical significance.

Further multivariable, ordinal logistic regression analysis favored a simpler model with the inclusion of BMI (High, OR: 0.05, 95% CI: 0.01–0.53), age (Moderate, OR: 1.16, 95% CI: 0.98–1.38 and High, OR: 10.06, 95% CI: 4.02–25.17), chest pain (Moderate, OR: 68.65, 95% CI: 1.71–2.75 × 103) and PM2.5 (High, OR: 1.24, 95% CI: 1.07–1.44) as shown in Table 2. While age and chest pain exhibited statistical significance in the high and moderate categories, respectively, results should be interpreted with caution due to the substantially wide confidence intervals. We observe a consistent statistically significant and negative association of BMI with IPAG category.

Table 2

Multivariate, ordinal logistic regression results of final model.

Variable	Moderate		High
	OR (95% CI)	P values	OR (95% CI)	P values
BMI	0.87 (0.65–1.18)	0.37	0.05 (0.01–0.53)	0.01a
Age	1.16 (0.98–1.38)	0.09a	10.06 (4.02–25.17)	<0.001a
Chest Pain	68.65 (1.71–2.75 × 103)	0.03a	0.12 (6.33 × 10−7 – 2.28 × 104)	0.73
PM2.5	1.01 (0.94–1.08)	0.87	1.24 (1.07–1.44)	0.01a

BMI: Body Mass Index; OR: Odds Ratio; 95% CI: 95% Confidence Interval; PM2.5: particulate matter with a size of 2.5 microns or less in diameter.

P value ≤ 0.10 indicate statistical significance.

Discussion

In this study, we examined the effects of ambient air pollution on the risk of COPD and observed statistically significant association between PM2.5 exposure and the risk of COPD among MMDA traffic enforcers, after adjusting for relevant covariates. Several risk and protective factors, aside from PM2.5, were also observed.

PM2.5 and COPD

We observed statistically significant association between PM2.5 exposure and the risk of COPD. A unit increase in PM2.5 exposure was associated with a statistically significant OR in high-risk category (OR: 1.24, 95% CI: 1.07–1.44) relative to the low-risk category, in Table 2. Among high risk populations, such as traffic enforcers, studies have indicated a decline in the workers’ lung function due to traffic-generated PM2.5. In India, taxi motorcyclists were observed to have substantial risk in developing restrictive syndrome, a form of impairment in the pulmonary function. In Zagreb and Croatia, Zuskin et al noted the higher risks in developing respiratory symptoms for bus drivers and mechanics. Epidemiological studies have demonstrated an association between short-term ambient PM2.5 exposure and lung function decline. Short-term exposure to ambient air pollutants are said to result in oxidative stress and pulmonary and systemic inflammatory responses such as reduction in the ciliary activity in the airways and acute decrease of pulmonary function. Schikowski et al argued that if repeated acute exacerbations of COPD were considered a factor in the development of COPD, then the causal role of air pollution in aggravation of COPD is substantial. That being said, the effect of long-term exposure to PM2.5 is slowly but considerably being investigated and has been additionally associated with increased morbidity and mortality of COPD. A recent longitudinal cohort study by Guo et al found that PM2.5 was consistently associated with reduction in lung function with respect to four parameters—namely FEV1, FVC, MMEF, and FEV1/FVC ratio—as well as increased risk of COPD development. Prolonged and chronic exposure may potentially induce pulmonary inflammation, decrease lung function, and cause emphysematous changes. Likewise, this study supports the reported specificity and sensitivity of IPAG Questionnaire as the most widely validated screening questionnaire for COPD.

COPD risk factors: age and chest pain

Age was positively associated with the risk of COPD in the high-risk category (OR: 10.06, 95% CI: 4.02–25.17) compared to the low-risk category; Table 2. The prevalence and manifestations of COPD increases with age for both men and women throughout most of the lifespan, especially among those above 40–50 years old. Since COPD is a progressive lung disease, symptoms worsen over time and often do not appear until there is significant lung damage. Aside from age, we also observed that symptom of chest pain was associated with statistically significant, yet with highly uncertain effects estimates. HajGhanbari et al noted that the hyper-expanded and relatively rigid chest wall, which affects the thoracic articulations in a hyper extended position, may contribute to thoracic pain experienced by people with COPD. Whereas, Bordoni et al posits the role of costal muscles, which face an alteration of their physiologic length (shortening or lengthening) caused by the change of costal biomechanics, influenced in turn by lung movements. However, due to the relatively uncertain effects estimate of chest pain, further studies are warranted.

COPD protective factor: BMI

Aside from risk factors, we observed statistically significant protective factor of BMI on the risk of COPD. A unit increase in BMI was associated with a decrease in the risk for COPD among high-risk category (OR: 0.05, 95%CI: 0.01–0.53). Galesanu et al observed similar findings noting that overweight/obese patients had better lung function and a larger midthigh muscle cross-sectional area obtained using computed tomography (MTCSACT) than those with normal BMI as well as significantly higher peak work rate than patients with normal BMI. A likely mechanism proposed by Divo et al posits the role of adipose tissue plays a modulating role in disease expression in the two (normal and obese) groups as the disease begins and progresses. Sun et al, on the other hand, postulates a possible reverse-causation, where increased COPD disease activity or severity leads to weight loss and cachexia. We invite the readers to read the systematic review by Sun et al for a detailed discussion regarding the inverse association of BMI and COPD.

The study recognizes the following limitations. First, this study only included ambient air quality data from designated PM2.5 sensors of the DENR-AQMS. Indoor air pollution was not quantified. Other gaseous pollutants (e.g. nitric oxide) aside from PM2.5 were not taken into account. Second, the developed Exposure-Month Index incorporates only the three most recent years of the participants' working history. Additionally, with a cross-sectional study design, exposure and outcomes were assessed simultaneously, thereby only capturing a snapshot of the participants’ health. Third, although the IPAG Questionnaire is a widely used screening tool, it does not account other risk factors of COPD. Additional information regarding respiratory health as well as candidate risk factors is encouraged. Fourth, because of the number of associated risk factors, a larger sample size may be more appropriate for future studies. Fifth, this study does not aim to clinically diagnose COPD. The clinical diagnosis of COPD requires a more extensive history taking and physical examination by a certified and experienced medical practitioner. Other subsequent tests may also aid in the accurate diagnosis of COPD, such as a post-bronchodilator spirometry assessment, arterial blood gas, and chest radiograph, all of which our study did not include.

In conclusion, exposure to PM2.5 was associated with an increased risk of COPD among high-risk category MMDA traffic enforcers. While age and symptom of chest pain may potentially pose a risk in developing COPD, caution should be exercised due to substantial uncertainty in the effect estimates. On the other hand, BMI exhibited a potential protective effect on the risk of COPD. Results of this study can be used for clinical management of high-risk populations, such that of MMDA traffic enforcers.

Conflicts of interest

None.