Next Article in Journal
Breastfeeding in the Community—How Can Partners/Fathers Help? A Systematic Review
Previous Article in Journal
Arsenic Concentrations and Dietary Exposure in Rice-Based Infant Food in Australia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Ambient Air Pollution on Rheumatoid Arthritis Disease Activity Score Index

by
Ahmad Alsaber
1,*,
Jiazhu Pan
1,*,
Adeeba Al-Herz
2,
Dhary S. Alkandary
3,
Adeeba Al-Hurban
3,
Parul Setiya
4 and
on behalf of the KRRD Group
1
Department of Mathematics and Statistics, University of Strathclyde, Glasgow G1 1XH, UK
2
Department of Rheumatology, Al-Amiri Hospital, P.O. Box 4077, Safat 13041, Kuwait
3
Department of Earth and Environmental Sciences, Kuwait University, P.O. Box 25944, Safat 1320, Kuwait
4
Department of Agrometerology, Govind Ballabh Pant University of Agriculture & Technology, Pantnagar, Uttarakhand 263145, India
*
Authors to whom correspondence should be addressed.
Kuwait Registry for Rheumatic Diseases (KRRD) Group Members: Khulood Saleh, Adel Al-Awadhi, Adeeba Al-Herz, Waleed Al-Kandari, Eman Hasan, Aqeel Ghanem, Mohammed Hussain, Yaser Ali, Ebrahim Nahar, Ahmad Alenizi, Sawsan Hayat, Fatemah Abutiban, Ali Aldei, Hebah Alhajeri, Ali Al-Sughayer, Naser Alhadhood, Khaled Mokaddam, Ahmad Khadrawy, Youssef Bartella.
Int. J. Environ. Res. Public Health 2020, 17(2), 416; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17020416
Submission received: 19 October 2019 / Revised: 30 December 2019 / Accepted: 3 January 2020 / Published: 8 January 2020
(This article belongs to the Section Environmental Health)

Abstract

:
Rheumatoid arthritis (RA) is a chronic autoimmune of an unknown etiology. Air pollution has been proposed as one of the possible risk factors associated with disease activity, although has not been extensively studied. In this study, we measured the relationship between exposure to air pollutants and RA activity. Data on RA patients were extracted from the Kuwait Registry for Rheumatic Diseases (KRRD). Disease activity was measured using disease activity score with 28 examined joints (DAS-28) and the Clinical Disease Activity Index (CDAI) during their hospital visits from 2013 to 2017. Air pollution was assessed using air pollution components (PM 10 , NO 2 , SO 2 , O 3 , and CO). Air pollution data were obtained from Kuwait Environmental Public Authority (K-EPA) from six different air quality-monitoring stations during the same period. Multiple imputations by the chained equations (MICE) algorithm were applied to estimate missing air pollution data. Patients data were linked with air pollution data according to date and patient governorate address. Descriptive statistics, correlation analysis, and linear regression techniques were employed using STATA software. In total, 1651 RA patients with 9875 follow-up visits were studied. We detected an increased risk of RA using DAS-28 in participants exposed to SO 2 and NO 2 with β = 0 . 003 (95% CI: 0.0004–0.005, p < 0 . 01 ) and β = 0 . 003 (95% CI: 0.002–0.005, p < 0 . 01 ), respectively, but not to PM 10 , O 3 , and CO concentrations. Conclusively, we observed a strong association between air pollution with RA disease activity. This study suggests air pollution as a risk factor for RA and recommends further measures to be taken by the authorities to control this health problem.

1. Introduction

Rheumatoid arthritis (RA) is a chronic autoimmune disease that mainly affects the joints causing inflammation, pain, and difficulty to use the joints. Although the exact cause is unknown, many genetic and environmental factors have been linked to the disease. Exposure to chemicals has previously been proposed as a possible, if not the main cause of the disease [1].
Several studies suggested that exposure of air pollution may increase the risk of RA [2,3]. Furthermore, epidemiological evidence indicates a significant association between the risk of RA and exposure to environmental factors, such as cigarette smoke, dioxin, noise, and traffic-related air pollution [4,5,6].
Pollutants with the strongest evidence for public health concern include particulate matter (PM), ozone (O 3 ), nitrogen dioxide (NO 2 ), and sulfur dioxide (SO 2 ) [7]. However, few studies have followed adequate methodologies correlating meteorological variables with RA. In this context, further investigation concerning the impact of air pollution on the risk of developing RA is still necessary.
With respect to RA, there is a need to investigate the impact of air pollution on RA through elaborate research. Air pollutants are part of the environmental components resulting from dust storms and fossil fuel combustion, which determine RA symptoms and worsen the overall disease. Additionally, gases such as SO 2 , NO 2 , carbon monoxide (CO), and O 3 are the other main pollutants that cause RA [8]. More in general, extensive investigations have been performed by several researchers about the impact of ambient air pollution on human health [9,10]. For instance, SO 2 resulting from the combustion of fossil fuels with high sulfur content is considered one of the most common pollutants with the worst impact on air quality. The combustion of fossil fuels containing high sulfur content causes the release of sulfur dioxide into the atmosphere. Several researchers have adequately documented the harmful effects of long-term exposure to high levels of SO 2 on overall health [11,12]; therefore, increased emission of pollutants into the atmosphere may potentially result in several adverse health effects, including RA.
Many composite indices for RA progress measurement are actually available: for example, the Disease Activity Score (DAS) with 28 examined joints (DAS-28) is one common RA index that has been extensively employed to identify the disease progress level for RA patients [13,14,15]. The Disease Activity Score was developed to measure and assess RA disease activity in daily clinical practice, clinical trials, and long-term observational studies [16]. The second RA index is the Clinical Disease Activity index (CDAI) [17,18], which is used to assess disease activity. The CDAI was developed to provide physicians and patients with simple and more understandable instruments.
In the present study, we aimed to investigate whether exposure to ambient air pollution (i.e., PM with aerodynamic diameter < 10 µm (PM 10 ), NO 2 , SO 2 , O 3 , and CO) is associated with an increased risk of RA using the DAS-28 and CDAI indices.

2. Materials and Methods

2.1. Patients Data Source (Kuwait Registry for Rheumatic Diseases—KRRD)

The State of Kuwait is a small country with a total area equal to 17,818 km 2 located on the far west side of the Asian continent. The total population of Kuwait is around 4.6 million, distributed into six main governorates [19]. All RA patients represented in this study were officially registered from the Kuwait Registry for Rheumatic Diseases (KRRD) from 2013 until the end of 2017. The KRRD is a national registry listing adult patients with rheumatic diseases. Patients who fulfilled the American College of Rheumatology (ACR) criteria for RA [20] registered from January 2013 to December 2017 were included in the study. The RA information data were collected from the rheumatology departments of four major government hospitals in Kuwait based on patient visits. The selected hospitals are mainly distributed in different governorates covering the ethnic diversity of the Kuwaiti population. The KRRD, from which this study originated, was approved by the Ethics Committees of the Faculty of Medicine at Kuwait University, and the Ministry of Health. Additionally, official consent was obtained from all represented patients enrolled in the registry [21].

2.2. Calculating RA Indices

RA disease activity scores are measured using two different indices: DAS-28 and CDAI. The DAS-28 is the sum of four outcome parameters: TJC28, the number of tender joints (0–28); SJC28, the number of swollen joints (0–28); ESR, the erythrocyte sedimentation rate (in mm/h) (C-reactive protein (CRP) may be used as an alternative to ESR in the calculation); and GH, the patient global health assessment (from 0 = best to 100 = worst) (see Equation (1)).
The second index is The Clinical Disease Activity Index (CDAI). CDAI takes into account the following items: TJC28, the number of tender joints (0–28); SJC28, the number of swollen joints (0–28); PaGH, the patient global health assessment (from 0 = best to 10 = worst); and PrGH, the care provider global health assessment (from 0 = best to 10 = worst) (see Equation (2)).
DAS 28 = 0 . 56 × T J C 28 + 0 . 28 × S J C 28 + 0 . 70 × ln ( E S R O r C R P ) + 0 . 014 × G H
CDAI = T J C 28 + S J C 28 + P a G H + Pr G H

2.3. Ambient Air Pollutants Data (Environmental Public Authority of Kuwait—K-EPA)

Pollutants data (PM 10 , NO 2 , SO 2 , O 3 , and CO) were obtained from six fixed monitoring stations run by the Environmental Public Authority of Kuwait (K-EPA). The air pollutant measurement sampling was from 1 January 2013 to 31 December 2017 based on hourly observations.
The pollutant data were distributed throughout the residential areas where stations were measuring different parameters including PM 10 , NO 2 , SO 2 , O 3 , and CO. The average concentration time of PM 10 , NO 2 , and SO 2 was twenty-four hours, and eight-hour time-average concentrations of O 3 and CO were collected for each measurement place, relying on 75% of the contributing values being present as the calculated averages.
According to Johnson et al. [22], the Air Quality Index (AQI) is defined as a measure of the condition of air relative to the requirements of one or more biotic species or to any human need [22,23]. The AQI is divided into categories, in which they are numbered, and each slot is marked with a color code. This provides a scale from a healthy level of zero to a very hazardous level of above 300 as a health risk indicator associated with air quality.
In this study, the air quality was assessed using the AQI developed by Al-Shayji et al. [24] for the State of Kuwait, based on the guidelines proposed by the United States Environmental Protection Agency (USEPA) [25]. The AQI is an index for reporting the day-to-day air quality. It gives details about the cleanliness of ambient air. The following equation was used to convert from pollutant concentration to AQI:
I p = I h i g h I l o w C h i g h C l o w C p C l o w + I l o w ,
where I p is the Air Quality Index for the given pollutant, C p is the pollutant concentration, C l o w is the concentration breakpoint that is C p , C h i g h is the concentration breakpoint that is C p , I l o w is the index breakpoint corresponding to C l o w , and I h i g h is the index breakpoint corresponding to C h i g h [26] (see Table 1).

2.4. Air Pollution Data Processing and Treatment

AQI data were examined by checking the normality assumption and detecting for any possible outliers before any statistical analysis or testing between the variables were done. About 5.8%, 1.6%, 48.4%, 5.3%, and 6.5% of data for PM 10 , NO 2 , SO 2 , O 3 , and CO were missing, respectively. Multiple imputations were performed to improve the accuracy of AQI prediction, where a final estimate was composed of the outputs of several multivariate fill-in methods [27,28].
The imputation process was determined by considering log transformation for better prediction to estimate and filling the missing data for each pollutant. Finally, AQI with RA information were matched with patients’ hospital location and date of patient visit. All AQI data for every pollutant were aggregated from hourly to daily observation.

2.5. Matching Procedure between Patients and AQI

All data management and combination was conducted using R Studio Version 1.1.463 running R 3.5.1 GUI 1.70 [29] software. Various R packages were used to clean, match, and combine the two datasets, including p l y r [30], d p l y r [31], t i d y r [32], and s t r i n g r [33].
The matching procedure was done using a developed R code to match between RA patient information and AQI monitoring station using date and governorate variables for both KRRD and K-EPA. As mentioned above, air pollution information was taken from six different stations distributed into all six governorates in the state of Kuwait. The matching procedure was conditioned on the date of patient visit with the date of the daily average AQI using a developed R code grouped by governorate physical address (e.g., if a patient lived in the Ahmadi governorate, the AQI information that came from the Ahmadi monitoring station was added to their visit information after matching the same date; see Figure 1).

2.6. Statistical Analysis

In the current study, means, standard deviations (SDs), and percentages were used to summarize and compare RA characteristics between the governorate levels. To estimate the association between the pollutants and RA indices, hierarchal linear model (HLM) analysis was performed using a regression approach. Pearson correlation test was performed to highlight the significant associations between AQI for NO 2 , CO, PM 10 , SO 2 , and O 3 and with the RA indices (DAS-28 and CDAI). Regression analyses were performed separately for DAS-28 and CDAI as response variables. Four regression models were estimated to highlight the most significant variables associated with the response variables (DAS-28 or CDAI). Because the variables belong to two different databases, we started with the first regression (M1) that measured the association between patients demographics (e.g., gender, RA disease duration, nationality, governorate, and comorbidity) with the response variables (DAS-28 or CDAI), and then we estimated the second regression model (M2) that measured the association between rheumatoid factors after adjusting for gender, RA disease duration, nationality, governorate, and comorbidity with the response variables (DAS-28 or CDAI). Models 1 and 2 were estimated from KRRD database. Then, we estimated Models 3 and 4 to highlight if there is any association between air pollution with disease activity indices (DAS-28 or CDAI) after merging the EPA with KRRD databases. Model 3 (M3) measured the direct effect from air pollutants to disease activity indices (DAS-28 or CDAI), whereas Model 4 (M4) measured the association between air pollutants to disease activity indices after adjusting for the rheumatoid factors mentioned in Model 2 (M2) (e.g., comorbidity, treatment class, swollen, tender, etc.). For better data fit, model comparison techniques using deviance scores were implemented to confirm the best choice of model [34,35]. All statistical procedures were performed using Stata 15.1 SE version software (StataCorp, College Station, TX, USA).
Model 1 (M1) was made to explain the influence of demographic variables (Disease Duration, Gender, Governorate, Nationality, Comorbidity, and Treatment Class) on the response variables (DAS-28 and CDAI). Model 2 (M2) was made to determine the effect of RA factors (Swollen, Tender, RF (rheumatoid factor), ACPA, Patient Global Assessment, and Physician Global Assessment) plus demographics on the response variables. Model 3 (M3) was made to estimate the relationship between the increase of NO 2 , CO, PM 10 , SO 2 , and O 3 concentrations and the response variables (DAS-28 and CDAI). Model 4 (M4) was made to explain the effect of AQI in terms of NO 2 , CO, PM 10 , SO 2 , and O 3 with RA factors (Swollen, RF, ACPA (anti-cyclic citrullinated peptide), ESR, and CRP) to indicate RA disease activity.

3. Results

The data of RA patients visits were obtained from KRRD, and the air pollution data were obtained from K-EPA. The analysis was performed during the period from 2013 to the end of 2017. There were 1651 RA patients with 9875 follow-up visits and 13,152 daily air pollution records. Because of the matching process to combine the data from the air pollution dataset with the RA patient visits from the KRRD dataset, the final dataset had to meet the matching conditions (matching based on date and governorate) with a total of 9875 records.
Table 2 shows some information about RA patient characteristics group by governorate location in the state of Kuwait. From the results, most patients were from Fawaniya governorate (n = 4378 visits; 44.3%). Most of the patients belonged to the local country with Kuwaiti nationality (n = 5783 visits; 58.6%). Females accounted for the majority of total visits (n = 6008; 60.8%). The average RA disease duration for all patients was 9.82 years with a SD of 6.48 years.
Most of the patients were positive for rheumatoid factor (RF; n = 6881; 74.6%) and positive for anti-cyclic citrullinated peptide (ACPA; n = 4934; 60.5%). The majority of RA patients presented co-morbidities (e.g., hypertension, hyperlipidemia, diabetes mellitus, chronic kidney disease, coronary artery disease, cancer, or any other illness; n = 5393; 54.6%). From the results, most of the consumed drugs were biologics (n = 5214; 52.8%).
In Table 2, we can see that most RA patients visited Al-Amiri hospital, with total visits n = 5051 (51.1%). Al-Amiri hospital is the only hospital in the Kuwait City governorate. The second-most patients were from Farwaniya hospital, with total visits numbering n = 3981 (40.3%). Farwaniya Hospital is the only public general hospital in Al Farwaniyah governorate. In Kuwait, there are six governorates, and each governorate has only one hospital.
With regard to disease activity, results in Table 2 present the clinical features of all RA patient visits. The average DAS-28 score was 2.67 with SD 1.26, and the average score of CDAI was 6.24 with a SD of 9.96, both indicating a low disease activity. The average and SD scores for ESR and CRP were x ¯ = 27 . 19 mm/h and SD = 21 . 79 and x ¯ = 6 . 32 mg/L and SD = 4 . 85 , respectively, which were both within the normal ranges according to our laboratory. Moreover, the average and SD for swollen and tender joints for all RA patient visits were x ¯ = 0 . 69 and SD = 2 . 26 and x ¯ = 2 . 87 and SD = 5 . 60 , respectively.
Table 3 shows the average air pollutant concentrations using AQI scores. The mean and SD AQIs for PM 10 , CO, NO 2 , O 3 , and SO 2 were 167 . 62 ± 214 . 27 , 1 . 28 ± 0 . 61 , 47 . 98 ± 26 . 64 , 17 . 76 ± 8 . 94 , and 15 . 87 ± 17 . 66 , respectively. The mean exposure levels for PM 10 , CO, NO 2 , O 3 , and SO 2 in Group 2 were 158 . 51 ± 68 . 02 , 1 . 31 ± 0 . 44 , 42 . 74 ± 18 . 83 , 18 . 17 ± 7 . 56 , and 13 . 94 ± 12 . 04 ) , respectively. Figure 2 shows the monthly AQI average time series for PM 10 , CO, NO 2 , O 3 , and SO 2 .
Moreover, all pollutants distribution are positively skewed (i.e., the means are higher than the medians for all pollutants). However, log transformation was employed to all pollutants in the regression model for better and quality parameter estimation. Logarithmic transformation pulls extreme values in the pollutants into the normal distribution. In the regression modeling, the logarithmic transformation is not considered, as linear explanation of the regression model is made easier to transform all other variables in the model with the same logarithmic transformation. However, this is not necessary as long as all other variables achieved the normality assumption.
Correlation analysis using Pearson correlation test was conducted to highlight the significant relationship between the pollutants (PM 10 , NO 2 , SO 2 , O 3 , and CO). Table 4 shows significant positive correlations between NO 2 and CO ( r p = 0.22), NO 2 and SO 2 ( r p = 0.51), and O 3 and PM 10 ( r p = 0.08). Significant negative correlations were discovered between O 3 and NO 2 ( r p = −0.12), PM 10 and NO 2 ( r p = −0.12), PM 10 and SO 2 ( r p = −0.03), and PM 10 and CO ( r p = −0.05).
Table 4 presents the Pearson correlation coefficients between different air pollutants and RA variables. For the score of RA disease activity using the DAS-28 index, the correlation results showed a positive significant correlation with exposure of SO 2 using AQI ( r p = 0 . 07 ), and the same results were returned with the with the exposure to NO 2 using AQI ( r p = 0 . 07 ). As for particular pollutants, only Hart et al. [36] provided evidence of elevated risks for NO 2 and SO 2 , especially in terms of seronegative RA. Other pollutants (PM 10 , CO, and O 3 ) did not show any significant correlation. For the CDAI, the correlation results showed a positive significant correlation with exposure to SO 2 using AQI ( r p = 0 . 10 ), and the same results were returned with the exposure to NO 2 ( r p = 0 . 11 ). Other pollutants (PM 10 , CO, and O 3 ) did not show any significant correlation with CDAI. For PM 10 , Hart et al. [36] reported an elevated odds ratio (OR), which failed to reach statistical significance. The effects of both P M 2 . 5 and PM 10 were non-significant across the analyses, but were consistently more pronounced for seronegative RA.
Table 5 and Table 6 present the hierarchal linear model (HLM) analysis using multiple linear regression. The first model (M1) demonstrates the effect of patient demographics on RA disease activity, where DAS-28 and CDAI were the response variables (see Table 5 and Table 6). For both DAS-28 and CDAI, the demographics explained the disease activity scores ( R D A S 28 2 = 0 . 034 and R C D A I 2 = 0 . 030 ). RA Disease Duration was not significant for DAS-28, but it had a significant association with CDAI (see Table 6 (M1)). Model 2 (M2) shows that RA factors (Swollen, Tender, RF, ACPA, Patient Global Assessment, Physician Global Assessment, ESR, and CRP) plus demographics affected DAS-28 ( R 2 = 0 . 865 ) and CDAI ( R 2 = 0 . 924 ). In the CDAI model (M2), Gender and RA Disease Duration were not significant in explaining the CDAI. Model 3 (M3) demonstrated and highlighted the effects of gaseous pollutants (PM 10 , NO 2 , SO 2 , O 3 , and CO) using AQI on RA disease activity; only SO 2 and NO 2 were significant risk factors for RA patients using the information of DAS-28 ( R 2 = 0 . 007 ) and CDAI ( R 2 = 0 . 015 ). The final model demonstrated the effect of gaseous air pollutants with RA factors (Swollen, RF, ACPA, ESR, and CRP) on RA disease activity. The AQI of NO 2 and SO 2 still showed positive associations with disease activity performance of RA. The positive effects of NO 2 in Model 4 (M4) were β = 0 . 003 (95% CI: 0.002–0.005) and β = 0 . 048 (95% CI: 0.030–0.066) for DAS-28 and CDAI, respectively (e.g., for a 1 μ g / m 3 increase in daily concentration of NO 2 , DAS-28 index will increase by 0.002 (95%CI:0.002–0.005) and CDAI index will increase by 0.048 (95% CI: 0.030–0.066)), whereas, for SO 2 , the results showed a positive significant effect with β = 0 . 003 (95% CI: 0.0004–0.005) and β = 0 . 044 (95% CI: 0.018–0.070) for DAS-28 and CDAI, respectively (e.g., for 1 μ g / m 3 increase in daily concentration of SO 2 , DAS-28 index will increase by 0.003 (95%CI: 0.0004–0.005) and CDAI index will increase by 0.044 (95% CI: 0.018–0.070)).

4. Discussion

Air pollution is a major concern for human health, since it is known to trigger and/or induce several pathologies and subsequently to increase morbidity and mortality rates, particularly in the Middle Eastern countries such as Kuwait. Therefore, air pollution control is a crucial element that should be prioritized by governments. According to the World Health Organization (WHO), six major air pollutants including NO 2 , CO, PM 10 , SO 2 , O 3 , and lead ( P b ) were identified to be primary and secondary pollutants with health risks. Many studies have found connections between particulates in the air and rates of hospitalization, chronic obstructive pulmonary disease, and restricted activity due to illness [7].
RA is considered as the most common chronic systemic auto-immune disease affecting joints, musculoskeletal apparatus, and fibrous tissues [37], whose incidence is expected to follow a positive trend during the following years [38].
Our findings show a significant relationship between the concentration of certain air pollutants and RA disease activity: particularly, our results show that increased SO 2 and NO 2 concentrations may increase the risk of developing RA disease activity using the index DAS-28 or CDAI. These results agree with those reported in other studies [36,39].
Concerning PM 10 , O 3 , and CO, our results do not show any effect being significant enough for considering them as risk factors of RA disease activity in terms of DAS-28 index; however, O 3 was found to be a risk factor together with SO 2 and NO 2 in terms of predicting RA disease activity using CDAI. A study from Korea reported an increased risk for developing RA in adults exposed to increased O 3 and CO concentrations [40].
For better understanding, the major sources of air pollutants in Kuwait are oil refineries, traffic, and power plants (mostly using fossil fuels), which are thought to be the main sources of SO 2 and NO 2 in Kuwait City [41]. Indeed, urban pollution results from various sources, including industrial activities and traffic within cities: more specifically, the most prevalent pollutants generated from road transport are NO 2 , CO, SO 2 , volatile organic compounds (VOCs), and P M [42].
Particularly, Kuwait—with a population of over four million people and a fleet of more than two million vehicles, both growing rapidly—is experiencing increases in traffic volumes, trip frequency, and trip length [43]; therefore, the poor air quality outdoors is becoming a major concern, especially for people living in Kuwait City. Kuwait is a relatively small country, with a high quantity of fixed and mobile sources of different pollutants affecting air quality, especially with a high level of road traffic. As a major oil supplier in the region, fossil fuel combustion products emitted by power stations and water distillation plants are liable for the elevated atmospheric SO 2 concentrations, whereas VOCs and NO 2 concentrations are mainly increased by automobiles [41].
In the present study, the relation between air pollutants and RA disease activity was measured using several regression models, to ensure that this association was still present even after the addition of RA factors that were highly significant for DAS-28 and CDAI. From Model 1 for both DAS-28 and CDAI, demographics controlled the disease activity level; these results comply with patient-reported outcomes used by the National Databank for Rheumatic Diseases (NDB) [44,45].
The results of Model 1 did not show any significant evidence concerning the effect of disease duration on the disease activity level as measured by the DAS-28; this also agrees with another study [46]. Nonetheless, the results of Models 3 and 4 confirmed the existence of a significant association between exposure to SO 2 and NO 2 and increase of RA disease activity: more specifically, in Model 3, where the effect of air pollutants was presented without adding RA variables, SO 2 and NO 2 showed significant relationships with RA’s disease activity; in Model 4, where other RA factors were included, SO 2 and NO 2 remained risk factors for RA disease activity level.
The Swedish Epidemiological Investigation examined the impact of prolonged exposure to air pollution on the probability of having RA: from this research, there was no evidence of a higher risk of RA resulting from exposure to PM; additionally, the overall risk of RA mildly increased due to atmospheric pollutants. This research confirmed the presence of a higher risk of RA following increases in NO 2 and SO 2 concentrations in ambient air [36].
Among these various fields of research, a study using distance-to-road as a proxy component for exposure to traffic-associated pollution reported a substantial increase of RA risk in those residing within shorter distances from the road [4]. Moreover, a Canadian nested case-control study reported an augmented risk of developing RA following exposure to O 3 [5]. In another case-control study, a significant increase in the risk of RA was reported in line with increasing concentrations of NO 2 and SO 2 [36]. Furthermore, other studies suggested an increased risk of RA in subjects who were exposed to NO 2 , particularly in women [1], as well as positive associations between O 3 , CO, and NO 2 and RA incidence [47].
This study has some strengths such as combining patients records with air pollution concentration to initiate a complete dataset that could be used for future academic studies. In addition, dealing with missing values using MICE algorithm increased the accuracy of the estimated regression coefficients. This study used KRRD database that includes data on RA patients in Kuwait with all previous records. One of the Limitations in the study is related to the records of the follow-up visits. As the data were extracted from a registry, the number of hospital visits is not equal for all patients. It ranges from 1 to 49 visits. Because of this limitation, time series analysis during patients visits could not be performed to estimate lag effect between air pollution and RA disease activity. Single lag day effect or moving average of several previous days lag effect could not be investigated in this study because of the data layout and the study duration is very short (from 2013 to 2017). However, it could be developed by improving time series features in the future.

5. Conclusions

Air pollution was significantly associated with disease activity scores in RA patients. NO 2 and SO 2 were found to be significant risk factors for RA activity, with 7% positive correlation with disease activity indices DAS-28 and CDAI. Future research could also be based on time-series analysis by employing univariate or multivariate time series analysis. It is also recommended that researchers classify the data on air pollution and disease activity score using a cluster technique and perform an adequate cluster analysis on the data.

6. Ethical Approval

The authors obtained ethical approval from the Ethics Committees of The Faculty of Medicine, University of Kuwait; The Ministry of Health in Kuwait; and the Kuwait Institute for Medical Specialization (KIMS).

Author Contributions

A.A. developed the research methodology, analyzed the data using STATA, and finished writing the manuscript; J.P. contributed for review and supervision; A.A.-H. (Adeeba Al-Herz) provided significant details on RA characteristics and provided the database as she is the founder of KRRD system including her revision to the manuscript; D.S.A. defined air pollutants and performed the AQI calculation; A.A.-H. (Adeeba Al-Hurban) carried out the review of the environmental literature works with provision of the K-EPA data; and P.S. analyzed the data and made all figures, graphs, and tables. All authors have read and agreed to the published version of the manuscript.

Funding

This is an investigator study from KRRD registry. The KRRD registry is supported by unrestricted grants from Pfizer.

Acknowledgments

The authors would like to express thanks for the support given by the KRRD and Kuwait EPA for the data collection. Without their contribution and support, this research would not be possible.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chang, K.H.; Hsu, C.C.; Muo, C.H.; Hsu, C.Y.; Liu, H.C.; Kao, C.H.; Chen, C.Y.; Chang, M.Y.; Hsu, Y.C. Air pollution exposure increases the risk of rheumatoid arthritis: A longitudinal and nationwide study. Environ. Int. 2016, 94, 495–499. [Google Scholar] [CrossRef] [PubMed]
  2. Shen, L.; Zhang, H.; Zhou, X.; Liu, R. Association between polymorphisms of interleukin 12 and rheumatoid arthritis associated biomarkers in a Chinese population. Cytokine 2015, 28, 363–367. [Google Scholar] [CrossRef] [PubMed]
  3. Solus, J.F.; Chung, C.P.; Oeser, A.; Li, C.; Rho, Y.H.; Bradley, K.M.; Kawai, V.K.; Smith, J.R.; Stein, C.M. Genetics of serum concentration of IL-6 and TNFα in systemic lupus erythematosus and rheumatoid arthritis: A candidate gene analysis. Clin. Rheumatol. 2015, 34, 1375–1382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hart, J.E.; Laden, F.; Puett, R.C.; Costenbader, K.H.; Karlson, E.W. Exposure to traffic pollution and increased risk of rheumatoid arthritis. Environ. Health Perspect. 2009, 117, 1065–1069. [Google Scholar] [CrossRef] [Green Version]
  5. De Roos, A.J.; Koehoorn, M.; Tamburic, L.; Davies, H.W.; Brauer, M. Proximity to traffic, ambient air pollution, and community noise in relation to incident rheumatoid arthritis. Environ. Health Perspect. 2014, 122, 1075–1080. [Google Scholar] [CrossRef] [Green Version]
  6. Kobayashi, S.; Okamoto, H.; Iwamoto, T.; Toyama, Y.; Tomatsu, T.; Yamanaka, H.; Momohara, S. A role for the aryl hydrocarbon receptor and the dioxin TCDD in rheumatoid arthritis. Rheumatology 2008, 47, 1317–1322. [Google Scholar] [CrossRef] [Green Version]
  7. Organization, W.H. Air Quality Management: Air Quality Guidelines; World Health Organization: Geneva, Switzerland, 1999; Available online: www.who.int/peh/air/airqualitygd.htm (accessed on 16 April 2019).
  8. Sun, G.; Hazlewood, G.; Bernatsky, S.; Kaplan, G.G.; Eksteen, B.; Barnabe, C. Association between air pollution and the development of rheumatic disease: A systematic review. Int. J. Rheumatol. 2016, 2016, 5356307. [Google Scholar] [CrossRef] [Green Version]
  9. Bernatsky, S.; Smargiassi, A.; Barnabe, C.; Svenson, L.W.; Brand, A.; Martin, R.V.; Hudson, M.; Clarke, A.E.; Fortin, P.R.; van Donkelaar, A.; et al. Fine particulate air pollution and systemic autoimmune rheumatic disease in two Canadian provinces. Environ. Res. 2016, 146, 85–91. [Google Scholar] [CrossRef]
  10. Tobón, G.J.; Youinou, P.; Saraux, A. The environment, geo-epidemiology, and autoimmune disease: Rheumatoid arthritis. Autoimmun. Rev. 2010, 9, A288–A292. [Google Scholar] [CrossRef]
  11. Seinfeld, J.H. Air Pollution: Physical and Chemical Fundamentals; McGraw-Hill: New York, NY, USA, 1975. [Google Scholar]
  12. Scott, H.M.; Soskolne, C.L.; Martin, S.W.; Ellehoj, E.A.; Coppock, R.W.; Guidotti, T.L.; Lissemore, K.D. Comparison of two atmospheric-dispersion models to assess farm-site exposure to sour-gas processing-plant emissions. Prev. Vet. Med. 2003, 57, 15–34. [Google Scholar] [CrossRef]
  13. Prevoo, M.; Van’T Hof, M.; Kuper, H.; Van Leeuwen, M.; Van De Putte, L.; Van Riel, P. Modified disease activity scores that include twenty-eight-joint counts development and validation in a prospective longitudinal study of patients with rheumatoid arthritis. Arthritis Rheum. 1995, 38, 44–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. van Gestel, A.M.; Haagsma, C.J.; van Riel, P.L. Validation of rheumatoid arthritis improvement criteria that include simplified joint counts. Arthritis Rheum. 1998, 41, 1845–1850. [Google Scholar] [CrossRef]
  15. Fransen, J.; Creemers, M.; Van Riel, P. Remission in rheumatoid arthritis: Agreement of the disease activity score (DAS28) with the ARA preliminary remission criteria. Rheumatology 2004, 43, 1252–1255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Van Riel, P. The development of the disease activity score (DAS) and the disease activity score using 28 joint counts (DAS28). Clin. Exp. Rheumatol. 2014, 32, S65–S74. [Google Scholar]
  17. Smolen, J.; Breedveld, F.; Schiff, M.; Kalden, J.; Emery, P.; Eberl, G.; Van Riel, P.; Tugwell, P. A simplified disease activity index for rheumatoid arthritis for use in clinical practice. Rheumatology 2003, 42, 244–257. [Google Scholar] [CrossRef] [Green Version]
  18. Martins, F.M.; da Silva, J.A.P.; Santos, M.J.; Vieira-Sousa, E.; Duarte, C.; Santos, H.; Costa, J.A.; Pimentel-Santos, F.M.; Cunha, I.; Cunha Miranda, L.; et al. DAS28, CDAI and SDAI cut-offs do not translate the same information: Results from the Rheumatic Diseases Portuguese Register Reuma. pt. Rheumatology 2014, 54, 286–291. [Google Scholar] [CrossRef] [Green Version]
  19. For Civil Information (PACI). Directory of Population and Labor Force. Available online: https://www.paci.gov.kw/Home.aspx (accessed on 23 May 2019).
  20. Aletaha, D.; Neogi, T.; Silman, A.J.; Funovits, J.; Felson, D.T.; Bingham, C.O., III; Birnbaum, N.S.; Burmester, G.R.; Bykerk, V.P.; Cohen, M.D.; et al. 2010 rheumatoid arthritis classification criteria: An American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum. 2010, 62, 2569–2581. [Google Scholar] [CrossRef]
  21. Al-Herz, A.; Al-Awadhi, A.; Saleh, K.; Al-Kandari, W.; Hasan, E.; Ghanem, A.; Abutiban, F.; Alenizi, A.; Hussain, M.; Ali, Y.; et al. A comparison of rheumatoid arthritis patients in Kuwait with other populations: Results from the KRRD registry. J. Adv. Med. Med. Res. 2016, 14, 1–11. [Google Scholar] [CrossRef]
  22. Johnson, M.; Isakov, V.; Touma, J.; Mukerjee, S.; Özkaynak, H. Evaluation of land-use regression models used to predict air quality concentrations in an urban area. Atmos. Environ. 2010, 44, 3660–3668. [Google Scholar] [CrossRef]
  23. Johnson, D.L.; Ambrose, S.H.; Bassett, T.J.; Bowen, M.L.; Crummey, D.E.; Isaacson, J.S.; Johnson, D.N.; Lamb, P.; Saul, M.; Winter-Nelson, A.E. Meanings of environmental terms. J. Environ. Qual. 1997, 26, 581–589. [Google Scholar] [CrossRef]
  24. Al-Shayji, K.; Lababidi, H.; Al-Rushoud, D.; Al-Adwani, H. Development of a fuzzy air quality performance indicator. Kuwait J. Sci. Eng. 2008, 35, 101–126. [Google Scholar]
  25. Fitz-Simons, T. Guideline for Reporting of Daily Air Quality: Air Quality Index (AQI); Technical Report; Environmental Protection Agency, Office of Air Quality Planning: Washington, DC, USA, 1999. [Google Scholar]
  26. Reporting, A.Q.I.; Rule, F. Protection Agency. 1999. Available online: https://www3.epa.gov/airnow/40cfrpt58_aqi-reporting.pdf (accessed on 7 July 2019).
  27. Junninen, H.; Niska, H.; Tuppurainen, K.; Ruuskanen, J.; Kolehmainen, M. Methods for imputation of missing values in air quality data sets. Atmos. Environ. 2004, 38, 2895–2907. [Google Scholar] [CrossRef]
  28. Schafer, J.L. Analysis of Incomplete Multivariate Data; Chapman and Hall/CRC: London, UK, 1997. [Google Scholar]
  29. Core, R.; Team, R. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013; Available online: http://www.R-project.org (accessed on 16 May 2019).
  30. Wickham, H. The split-apply-combine strategy for data analysis. J. Stat. Softw. 2011, 40, 1–29. [Google Scholar] [CrossRef] [Green Version]
  31. Wickham, H.; Francois, R. dplyr: A Grammar of Data Manipulation (Version 0.8.3) [Software]. 2019. Available online: http://CRAN.R-project.org/package=dplyr (accessed on 7 June 2019).
  32. Wickham, H. tidyr: Easily Tidy Data with Spread () and Gather () Functions. R Package. Version 1.0.0. 2019. Available online: http://CRAN.R-project.org/package=tidyr (accessed on 7 June 2019).
  33. Wickham, H. stringr: Make it easier to work with strings. R Package Vers. 0.6 2012, 2, 96–97. [Google Scholar]
  34. Nelder, J.A.; Wedderburn, R.W. Generalized linear models. J. R. Stat. Soc. 1972, 135, 370–384. [Google Scholar] [CrossRef]
  35. Dobson, A.J.; Barnett, A.G. An Introduction to Generalized Linear Models; Chapman and Hall/CRC: London, UK, 2008. [Google Scholar]
  36. Hart, J.E.; Källberg, H.; Laden, F.; Costenbader, K.H.; Yanosky, J.D.; Klareskog, L.; Alfredsson, L.; Karlson, E.W. Ambient air pollution exposures and risk of rheumatoid arthritis. Arthritis Care Res. 2013, 65, 1190–1196. [Google Scholar] [CrossRef] [Green Version]
  37. Gabriel, S.E.; Crowson, C.S.; O’Fallon, M. The epidemiology of rheumatoid arthritis in Rochester, Minnesota, 1955–1985. Arthritis Rheum. 1999, 42, 415–420. [Google Scholar] [CrossRef]
  38. Chaudhari, K.; Rizvi, S.; Syed, B.A. Rheumatoid arthritis: Current and future trends. Nat. Rev. Drug Discov. 2016, 15, 305–306. [Google Scholar] [CrossRef]
  39. Hart, J.E.; Källberg, H.; Laden, F.; Bellander, T.; Costenbader, K.H.; Holmqvist, M.; Klareskog, L.; Alfredsson, L.; Karlson, E.W. Ambient air pollution exposures and risk of rheumatoid arthritis: Results from the Swedish EIRA case–control study. Ann. Rheum. Dis. 2013, 72, 888–894. [Google Scholar] [CrossRef]
  40. Shin, J.; Lee, J.; Lee, J.; Ha, E.H. Association between Exposure to Ambient Air Pollution and Rheumatoid Arthritis in Adults. Int. J. Environ. Res. Public Health 2019, 16, 1227. [Google Scholar] [CrossRef] [Green Version]
  41. Al-Awadhi, J. Measurement of air pollution in Kuwait City using passive samplers. Atmos. Clim. Sci. 2014, 4, 253–271. [Google Scholar] [CrossRef] [Green Version]
  42. Hankey, S.; Marshall, J.D.; Brauer, M. Health impacts of the built environment: Within-urban variability in physical inactivity, air pollution, and ischemic heart disease mortality. Environ. Health Perspect. 2011, 5, 247–253. [Google Scholar] [CrossRef]
  43. Al-Mutairi, N.; Koushki, P. Potential contribution of traffic to air pollution in the State of Kuwait. Am. J. Environ. Sci. 2009, 5, 218–222. [Google Scholar] [CrossRef]
  44. Covic, T.; Tyson, G.; Spencer, D.; Howe, G. Depression in rheumatoid arthritis patients: Demographic, clinical, and psychological predictors. J. Psychosom. Res. 2006, 60, 469–476. [Google Scholar] [CrossRef] [PubMed]
  45. Godha, D.; Shi, L.; Mavronicolas, H. Association between tendency towards depression and severity of rheumatoid arthritis from a national representative sample: The Medical Expenditure Panel Survey. Curr. Med. Res. Opin. 2010, 26, 1685–1690. [Google Scholar] [CrossRef] [PubMed]
  46. Gonzalez-Alvaro, I.; Ortiz, A.; Garcia-Vicuna, R.; Balsa, A.; Pascual-Salcedo, D.; Laffon, A. Increased serum levels of interleukin-15 in rheumatoid arthritis with long-term disease. Clin. Exp. Rheumatol. 2003, 21, 639–642. [Google Scholar]
  47. Jung, C.R.; Hsieh, H.Y.; Hwang, B.F. Air pollution as a potential determinant of rheumatoid arthritis: A population-based cohort study in Taiwan. Epidemiology 2017, 28, S54–S59. [Google Scholar] [CrossRef]
Figure 1. Matching Procedures to combine Air Quality Index (AQI) information with Kuwait Registry for Rheumatic Diseases (KRRD) patient profile records using date and governorate address information). K-EAPL, Environmental Public Authority of Kuwait; RA, rheumatoid arthritis.
Figure 1. Matching Procedures to combine Air Quality Index (AQI) information with Kuwait Registry for Rheumatic Diseases (KRRD) patient profile records using date and governorate address information). K-EAPL, Environmental Public Authority of Kuwait; RA, rheumatoid arthritis.
Ijerph 17 00416 g001
Figure 2. Air quality index (AQI) ambients for six governorate monitoring stations in Kuwait from 2013 to 2017.
Figure 2. Air quality index (AQI) ambients for six governorate monitoring stations in Kuwait from 2013 to 2017.
Ijerph 17 00416 g002
Table 1. Kuwait Air Quality Index. PM 10 : particulate matter with aerodynamic diameter ≤10 μ m.
Table 1. Kuwait Air Quality Index. PM 10 : particulate matter with aerodynamic diameter ≤10 μ m.
CategoriesAQI,
Sub-Index
O 3 (ppm),
8-h
PM 10 (µg/m 3 ),
24-h
CO (ppm),
24-h
SO 2 (ppm),
24-h
NO 2 (ppm),
24-h
I low I high I low I high I low I high I low I high I low I high I low I high
Good0–500.0–0.030.0–900.0–4.00.0–0.030.0–0.03
Moderate51–1000.031–0.0690.1–350.04.1–8.00.031–0.060.04–0.05
Unhealthy (1)101–1500.061–0.092350.1–431.18.1–11.70.061–0.1820.06–0.30
Unhealthy (2)151–2000.093–0.124431.4–512.511.8–15.40.183–0.3040.31–0.55
Very Unhealthy201–3000.125–0.374512.6–675.015.5–30.40.305–0.6040.56–1.04
Hazardous301–5000.375–0.504675.1–100030.5–50.40.605–1.0041.05–2.04
Table 2. RA patient visit demographic and clinical features groups by governorate—2012 to 2017.
Table 2. RA patient visit demographic and clinical features groups by governorate—2012 to 2017.
[ALL]AhmadiFarwaniyaHawallyJahraKuwait CityMubarak
n = 9875n = 356n = 4378n = 1272n = 226n = 3007Alkabeer n = 636
Gender:
  male3867 (39.2%)73 (20.5%)2656 (60.7%)285 (22.4%)82 (36.3%)653 (21.7%)118 (18.6%)
  female6008 (60.8%)283 (79.5%)1722 (39.3%)987 (77.6%)144 (63.7%)2354 (78.3%)518 (81.4%)
Nationality:
  Kuwaitis5783 (58.6%)328 (92.1%)1499 (34.2%)773 (60.8%)145 (64.2%)2462 (81.9%)576 (90.6%)
  non-Kuwaitis4092 (41.4%)28 (7.87%)2879 (65.8%)499 (39.2%)81 (35.8%)545 (18.1%)60 (9.43%)
Visited Hospital:
  Amiri5051 (51.1%)294 (82.6%)364 (8.31%)774 (60.8%)120 (53.1%)2955 (98.3%)544 (85.5%)
  Farwaniya3981 (40.3%)0 (0.00%)3976 (90.8%)5 (0.39%)0 (0.00%)0 (0.00%)0 (0.00%)
  Jahra111 (1.12%)0 (0.00%)2 (0.05%)3 (0.24%)106 (46.9%)0 (0.00%)0 (0.00%)
  Mubarak732 (7.41%)62 (17.4%)36 (0.82%)490 (38.5%)0 (0.00%)52 (1.73%)92 (14.5%)
Disease Duration 3 9.82 (6.48)13.3 (9.11)9.39 (5.79)10.4 (7.03)10.2 (6.45)9.42 (5.78)11.4 (9.60)
Comorbidity 2 :
  Yes5393 (54.6%)226 (63.5%)1658 (37.9%)802 (63.1%)123 (54.4%)2157 (71.7%)427 (67.1%)
  No4482 (45.4%)130 (36.5%)2720 (62.1%)470 (36.9%)103 (45.6%)850 (28.3%)209 (32.9%)
Treatment Class:
  Biologics5214 (52.8%)327 (91.9%)1342 (30.7%)762 (59.9%)131 (58.0%)2124 (70.6%)528 (83.0%)
  cDMARDs 4 4661 (47.2%)29 (8.15%)3036 (69.3%)510 (40.1%)95 (42.0%)883 (29.4%)108 (17.0%)
RF 1 :
  Positive6881 (74.6%)235 (72.8%)3148 (76.3%)817 (71.0%)177 (93.2%)2042 (72.0%)462 (76.5%)
  Negative2348 (25.4%)88 (27.2%)976 (23.7%)334 (29.0%)13 (6.84%)795 (28.0%)142 (23.5%)
ACPA:
  Positive4934 (60.5%)102 (33.6%)2665 (70.3%)593 (63.1%)73 (61.3%)1205 (48.1%)296 (60.0%)
  Negative3216 (39.5%)202 (66.4%)1125 (29.7%)347 (36.9%)46 (38.7%)1299 (51.9%)197 (40.0%)
Patient GA1.64 (2.36)1.54 (2.34)1.02 (1.85)2.60 (2.69)2.09 (2.52)1.95 (2.56)2.39 (2.69)
Physician GA1.05 (1.77)1.06 (1.82)0.72 (1.50)1.63 (2.02)1.58 (2.21)1.13 (1.78)1.64 (2.18)
DAS-282.67 (1.26)1.85 (1.35)2.70 (1.21)2.77 (1.29)3.04 (1.39)2.61 (1.22)2.79 (1.44)
CDAI6.24 (9.96)4.64 (8.89)4.83 (8.56)8.31(10.72)9.45 (14.25)6.78 (10.49)9.00 (11.53)
ESR27.19 (21.79)15.12 (16.82)30.24 (23.10)26.77 (20.30)30.06 (19.33)23.97 (19.54)27.80 (24.04)
CRP6.32 (4.85)4.32 (3.91)7.29 (4.89)4.53 (4.53)6.36 (4.42)6.08 (4.68)5.38 (4.87)
Swollen Joints0.69 (2.26)0.34 (1.57)1.08 (2.60)0.53 (1.97)0.95 (3.63)0.26 (1.67)0.57 (1.99)
Tender Joints2.87 (5.60)1.72 (4.32)2.02 (4.18)3.55 (6.21)4.82 (8.39)3.46 (6.36)4.54 (7.13)
1 RF, rheumatoid factor; ACPA, anti-cyclic citrullinated peptide antibody. 2 Comorbidity (e.g., hypertension, hyperlipidemia, diabetes mellitus, etc.); 3 Disease Duration, RA disease duration by years. 4 DMARDs, conventional disease modifying anti-rheumatic drugs. GA, global assessment.
Table 3. Distribution of Kuwait ambient air pollution exposure using AQI during 2012–2017.
Table 3. Distribution of Kuwait ambient air pollution exposure using AQI during 2012–2017.
Air PollutantAhmadi (n = 356)Farwaniya (n = 4378)Hawally (n = 1272)Jahra (n = 226)Kuwait City (n = 3007)Mubarak Alkabeer (n = 636)ALL (n = 9875)
PM 10
  min17.10820.34621.94818.83712.1385.4215.421
  25th 1 75.83976.31280.69267.41771.11474.88674.965
  median112.623120.779113.58692.788116.081109.917113.586
  75th 2 180.403186.145180.581151.762193.121184.282186.568
  max549.419511.826588.494545.789577.706585.077588.494
  mean (SD 3 )142.36 ± 99.58146.47 ± 94.20145.69 ± 99.23123.64 ± 95.00146.81 ± 102.42142.40 ± 102.79144.87 ± 100.64
CO
  min0.2400.2070.0870.0420.0870.2920.042
  25th0.9380.9450.9840.8541.0060.9780.975
  median1.3671.3381.3371.1511.3691.3801.346
  75th1.6791.6031.6721.4551.6791.7281.672
  max4.8944.7014.4715.1228.1435.2878.143
  mean (SD)1.37 ± 0.621.33 ± 0.571.40 ± 0.661.14 ± 0.641.40 ± 0.661.46 ± 0.701.39 ± 0.66
NO 2
  min9.0809.0805.3468.2295.3465.3465.346
  25th25.76828.20027.20829.12727.31932.03727.391
  median35.76236.53735.81052.42036.79544.74437.355
  75th54.21853.15051.40967.86652.05468.21054.552
  max137.960135.878134.123107.279208.670207.557208.670
  mean (SD)42.85 ± 24.1342.01 ± 20.7041.29 ± 20.6050.92 ± 25.5443.00 ± 22.7652.09 ± 27.9143.74 ± 23.13
O 3
  min4.0514.0514.0515.8873.4764.8773.476
  25th10.96511.03010.85112.45710.58111.32810.851
  median15.21515.37215.10418.19214.70914.16414.985
  75th20.87422.24020.45125.17220.45119.50720.759
  max54.26269.68280.65637.53988.62357.33088.623
  mean (SD)17.04 ± 8.7518.53 ± 11.2717.58 ± 10.5218.71 ± 7.2416.85 ± 10.0816.26 ± 7.5117.19 ± 9.90
SO 2
  min0.0031.0001.0000.6650.0031.0000.003
  25th4.2085.2934.8755.4354.4907.3334.875
  median8.3338.0008.59413.7927.99314.0838.727
  75th14.14617.29217.16924.58316.74622.94617.504
  max121.833121.833121.83376.875111.917127.875127.875
  mean (SD)13.15 ± 15.4613.26 ± 14.2214.18 ± 16.1017.90 ± 16.5213.39 ± 14.6018.62 ± 17.0714.26 ± 15.42
1 25th, lower quartile (25th percentile). 2 75th, upper quartile (75th percentile). 3 SD: standard deviation.
Table 4. Correlation analysis between rheumatoid arthritis disease factors and AQI for SO 2 , NO 2 , CO, O 3 , and PM 10 .
Table 4. Correlation analysis between rheumatoid arthritis disease factors and AQI for SO 2 , NO 2 , CO, O 3 , and PM 10 .
DAS-28CDAINO 2 O 3 SO 2 COPM 10 SwollenTenderESR
DAS-28
CDAI0.77 ****
NO 2 0.07 ****0.11 ****
O 3 0.000.00−0.12 ****
SO 2 0.07 ****0.10 ****0.51 ****−0.09 ****
CO−0.010.020.22 ****0.020.07 ****
PM 10 0.00−0.02−0.12 ****0.08 ****−0.03 *−0.05 **
Swollen0.50 ****0.60 ****0.010.010.010.00−0.02
Tender0.72 ****0.93 ****0.13 ****0.010.11 ****0.03−0.010.42 ****
ESR0.65 ****0.20 ****0.00−0.020.04 *−0.04 *0.020.16 ****0.17 ****
CRP0.28 ****0.02 *0.010.020.01−0.01−0.010.11 ****0.02 *0.37 ****
Note: * p < 0.1; ** p < 0.05; *** p < 0.01; **** p < 0.001.
Table 5. Coefficients estimated by hierarchal linear model (HLM) for DAS-28 (standard error in parentheses).
Table 5. Coefficients estimated by hierarchal linear model (HLM) for DAS-28 (standard error in parentheses).
Dependent Variable
DAS-28
(M1)(M2)(M3)(M4)
Gender (male)−0.213 ***−0.040 ***
(−0.268, −0.157)(−0.064, −0.017)
RA Disease Duration−0.002−0.004 ***
(−0.006, 0.002)(−0.006, −0.003)
Nationality (non-Kuwaitis)0.272 ***0.022
(0.214, 0.331)(−0.007, 0.051)
Governorate (Farwaniya)0.807 ***0.299 ***
(0.666, 0.947)(0.239, 0.358)
Governorate (Hawally)0.852 ***0.277 ***
(0.704, 1.000)(0.214, 0.341)
Governorate (Jahra)1.143 ***0.254 ***
(0.933, 1.352)(0.150, 0.359)
Governorate (Kuwait City)0.744 ***0.290 ***
(0.606, 0.882)(0.232, 0.348)
Governorate (Mubarak Alkabeer)0.955 ***0.175 ***
(0.793, 1.117)(0.106, 0.244)
Comorbidity (Yes)0.060 **−0.051 ***
(0.007, 0.114)(−0.074, −0.029)
Treatment Class (cDMARDs) 0.064 ***
(0.036, 0.092)
Swollen 0.090 *** 0.226 ***
(0.085, 0.095) (0.208, 0.244)
Tender 0.099 ***
(0.096, 0.101)
RF (Positive) 0.035 *** 0.004
(0.010, 0.060) (−0.078, 0.085)
ACPA (Positive) 0.008 0.007
(−0.015, 0.031) (−0.063, 0.078)
Patient Global Assessment 0.097 ***
(0.088, 0.105)
Physician Global Assessment 0.014 **
(0.002, 0.025)
ESR 0.028 *** 0.035 ***
(0.027, 0.028) (0.034, 0.037)
CRP 0.017 *** 0.001
(0.015, 0.020) (−0.006, 0.009)
NO 2 0.003 **0.003 ***
(0.001, 0.005)(0.002, 0.005)
O 3 0.0020.003
(−0.002, 0.006)(−0.001, 0.006)
SO 2 0.004 ***0.003 **
(0.001, 0.007)(0.0004, 0.005)
CO −0.051−0.001
(−0.114, 0.012)(−0.053, 0.052)
PM 10 0.00020.00003
(−0.0002, 0.001)(−0.0003, 0.0004)
Constant1.845 ***1.029 ***2.586 ***1.506 ***
(1.701, 1.988)(0.966, 1.093)(2.435, 2.738)(1.358, 1.654)
R 2 0.0340.8650.0070.488
Adjusted R 2 0.0330.8650.0060.486
Note: ** p < 0.05; *** p < 0.01.
Table 6. Coefficients estimated through HLM for CDAI (standard error in parentheses).
Table 6. Coefficients estimated through HLM for CDAI (standard error in parentheses).
Dependent Variable
CDAI
(M1)(M2)(M3)(M4)
Gender (male)−0.748 ***−0.078
(−1.187, −0.309)(−0.218, 0.061)
RA Disease Duration−0.052 ***0.005
(−0.082, −0.021)(−0.005, 0.015)
Nationality (non-Kuwaitis)1.209 ***0.255 ***
(0.747, 1.671)(0.081, 0.429)
Governorate (Farwaniya)0.066−1.199 ***
(−1.044, 1.176)(−1.550, −0.848)
Governorate (Hawally)3.408 ***0.526 ***
(2.241, 4.576)(0.152, 0.900)
Governorate (Jahra)4.662 ***−0.559 *
(3.008, 6.315)(−1.179, 0.062)
Governorate (Kuwait City)1.947 ***−0.135
(0.858, 3.036)(−0.476, 0.207)
Governorate (Mubarak Alkabeer)4.430 ***−0.064
(3.150, 5.709)(−0.474, 0.345)
Comorbidity (Yes)0.984 ***0.147 **
(0.564, 1.405)(0.014, 0.281)
Treatment Class (cDMARDs) −0.059
(−0.225, 0.107)
Swollen 1.145 *** 3.035 ***
(1.114, 1.175) (2.848, 3.222)
Tender 1.423 ***
(1.410, 1.435)
RF (Positive) 0.025 −0.066
(−0.123, 0.173) (−0.901, 0.770)
ACPA (Positive) 0.226 *** −0.594
(0.090, 0.363) (−1.320, 0.132)
Patient Global Assessment 1.067 ***
(1.059, 1.075)
Physician Global Assessment 0.871 ***
(0.861, 0.881)
ESR 0.019 *** 0.085 ***
(0.016, 0.022) (0.068, 0.103)
CRP −0.050 *** −0.171 ***
(−0.064, −0.037) (−0.246, −0.095)
NO 2 0.040 ***0.048 ***
(0.022, 0.058)(0.030, 0.066)
O 3 0.0270.039 **
(−0.008, 0.062)(0.003, 0.074)
SO 2 0.044 ***0.044 ***
(0.018, 0.070)(0.018, 0.070)
CO −0.0620.185
(−0.601, 0.476)(−0.358, 0.729)
PM 10 0.0010.0004
(−0.003, 0.004)(−0.003, 0.004)
Constant4.537 ***1.322 ***5.306 ***2.540 ***
(3.404, 5.671)(0.948, 1.697)(4.006, 6.606)(1.017, 4.064)
R 2 0.0300.9250.0150.299
Adjusted R 2 0.0290.9250.0140.297
Note: * p < 0.1; ** p < 0.05; *** p < 0.01.

Share and Cite

MDPI and ACS Style

Alsaber, A.; Pan, J.; Al-Herz, A.; Alkandary, D.S.; Al-Hurban, A.; Setiya, P.; on behalf of the KRRD Group. Influence of Ambient Air Pollution on Rheumatoid Arthritis Disease Activity Score Index. Int. J. Environ. Res. Public Health 2020, 17, 416. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17020416

AMA Style

Alsaber A, Pan J, Al-Herz A, Alkandary DS, Al-Hurban A, Setiya P, on behalf of the KRRD Group. Influence of Ambient Air Pollution on Rheumatoid Arthritis Disease Activity Score Index. International Journal of Environmental Research and Public Health. 2020; 17(2):416. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17020416

Chicago/Turabian Style

Alsaber, Ahmad, Jiazhu Pan, Adeeba Al-Herz, Dhary S. Alkandary, Adeeba Al-Hurban, Parul Setiya, and on behalf of the KRRD Group. 2020. "Influence of Ambient Air Pollution on Rheumatoid Arthritis Disease Activity Score Index" International Journal of Environmental Research and Public Health 17, no. 2: 416. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17020416

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop