Relationship Between Exposure to Dusty Days, Temperature, and Wind Speed and the Prevalence of Multiple Sclerosis in Khuzestan Province, Iran

authors:

avatar Maryam Dastoorpour ORCID 1 , avatar Hossein Mirshekarpour 2 , avatar Dadkhoda Sofi 3 , avatar Khosro Jamebozorgi 4 , *

Department of Biostatistics and Epidemiology, School of Public Health, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
Department of Radiology, Afzalipour Hospital, School of Medicine, Kerman University of Medical Sciences, Kerman, Iran
Department of Internal Medicine, Amir-Al-Momenin Hospital, School of Medicine, Zabol University of Medical Sciences, Zabol, Iran
Department of Neurology, Amir-Al-Momenin Hospital, School of Medicine, Zabol University of Medical Sciences, Zabol, Iran

How To Cite? Dastoorpour M, Mirshekarpour H, Sofi D, Jamebozorgi K. Relationship Between Exposure to Dusty Days, Temperature, and Wind Speed and the Prevalence of Multiple Sclerosis in Khuzestan Province, Iran. Jundishapur J Chronic Dis Care. 2025;14(1):e155650. https://doi.org/10.5812/jjcdc-155650.

Abstract

Background:

Various factors, including environmental influences, may contribute to the onset of multiple sclerosis (MS). Khuzestan province is notable for being one of the hottest and most dust-storm-affected regions in the Middle East.

Objectives:

This ecological study aimed to examine the relationship between exposure to dusty days, temperature, and wind speed and the prevalence of MS in cities across Khuzestan province, southern Iran.

Methods:

This ecological research was conducted in 2019 in Khuzestan province and utilized available data on MS prevalence, dusty days, temperature, and wind speed. The data were obtained from the Khuzestan MS Center database, the Environmental Protection Agency, and the Meteorological Organization of Khuzestan province. A multi-pollutant generalized additive model (GAM) was applied to evaluate the relationship between the independent variables and MS prevalence. Data analysis was performed using R software version 3.4.3.

Results:

The average age of the patients was 31.40 ± 8.94 years, with a female-to-male sex ratio of 3.05. In 2019, the average annual prevalence of MS across the 27 counties in Khuzestan province was 42 per 100,000 people. The GAM analysis identified a significant positive association between dusty days and MS prevalence [relative risk (RR) = 1.113, 95% CI: 1.090 - 1.137, P < 0.001] and a significant negative association between average wind speed and MS prevalence (RR = 0.999, 95% CI: 0.9990 - 0.9997, P < 0.001).

Conclusions:

The findings suggest that an increase in dusty days is linked to higher MS prevalence, whereas higher average wind speeds are associated with lower MS prevalence. This study contributes to a deeper understanding of the role of environmental factors in the prevalence of MS and its management.

1. Background

While multiple sclerosis (MS) is the most common cause of non-traumatic neurological disability in young adults (1), environmental factors are increasingly recognized for their role in disease onset, particularly in regions like Khuzestan. The prevalence of this disease is twice as high in women as in men, and the average age of onset is 30 years (2). It is assumed that, in addition to genetics, environmental factors may play a significant role in the etiology of the disease, as geographical location has a substantial impact on the prevalence of MS. The MS is considered uncommon in tropical and subtropical regions (1). According to statistics, Africa and Asia have the lowest prevalence rates, but recent studies suggest an increase in prevalence in these regions (3).

Although Iran has historically had a low prevalence of MS, it has recently become a country with a high prevalence, and the numbers continue to rise (4, 5). In 2013, while the global prevalence of MS was reported as 73 per 100,000, Iran reported a prevalence of 60 per 100,000 (6). In 2015, Dehghani et al. observed an increasing trend in all provinces of Iran over a five-year study period (7). Several review studies have reported MS prevalence rates in Iran ranging from 5.3 to 89 per 100,000 and incidence rates ranging from 7 to 148.1 per 100,000 (4, 6).

The exact cause of the increase in MS prevalence remains uncertain. The role of environmental factors in MS susceptibility has been explored in limited studies. Among these factors, the effect of air temperature on the disease is noteworthy. A study by Ebrahimi and Sedighi in various cities of Kerman province found a linear relationship between air temperature and MS prevalence, with higher prevalence reported in colder climates (8). Similarly, Spelman, studying 30 countries worldwide while considering geographical latitude, demonstrated a pronounced seasonal recurrence pattern with an annual sinusoidal trend. This pattern peaked in early spring and declined in autumn in the northern hemisphere, with extensions observed in the southern hemisphere (9).

Furthermore, a study analyzing monthly average meteorological data revealed significant seasonal variations in ambient temperature, relative humidity, and wind speed. The monthly onset distribution of MS in the study group was highest in May and June (12.2% each) and lowest in September (3.7%) (10). However, contrary findings were reported in Israel, where Saaroni et al. observed no significant effect of meteorological parameters, season, month, or specific days on MS relapse occurrence (11).

Air pollution is another critical environmental factor linked to MS. In Tehran, Heydarpour et al. demonstrated that long-term exposure to air pollutants, including PM10, SO2, NO, NO2, and NOx, may act as risk factors for MS (12). Supporting this, a similar study conducted in Italy in 2016 suggested that air pollution might contribute to MS onset and relapse. Increased exposure to PM10 was associated with higher hospitalization rates for the disease (13).

Although the exact cause of MS remains unknown, a combination of factors, including environmental influences, infections, genetic predispositions, immune system dysfunctions, viruses, toxins, metabolic issues, and substances like nitric oxide, may play a role in its development.

2. Objectives

The 2013 annual report by the World Health Organization identified Ahvaz as the most polluted city in the world, with an annual mean PM10 concentration of 372 μg/m³ (14-16). Due to industrialization, global warming, and Khuzestan province's exposure to dusty weather conditions, this study aimed to quantitatively assess the relationship between exposure to dusty days, temperature, and wind speed and the prevalence of MS, focusing on identifying critical thresholds. Epidemiological information can assist in the rapid identification, diagnosis, and control of the consequences of this disease (17).

3. Methods

Khuzestan province, located in southwest Iran along the Persian Gulf coast, is a key center for the country’s oil and gas extraction activities. Covering an area of 64,055 square kilometers and housing a population of over 4,711 million (based on the 2016 census), it is the fifth most populous province in Iran. Geographically, Khuzestan is situated at coordinates 31.33˚N and 48.69˚E (18).

This ecological study was conducted in Khuzestan province in 2019 using available data. Ecological studies are particularly suitable for evaluating population-level associations between environmental factors and disease prevalence.

The dataset included the prevalence of MS, the number of dusty days, the average daily temperature, and wind speed for the 27 counties of Khuzestan province, Iran, all collected in 2019.

Participants were identified using data from the Khuzestan MS Center, affiliated with Ahvaz Jundishapur University of Medical Sciences, to ensure a representative sample of the MS population in the region. This center serves as a referral hub for MS patients in Khuzestan province, with initial diagnoses made by neurologists. Patient evaluations at the center are overseen by neurologists, with final diagnoses based on clinical and radiological evidence. Data accuracy was ensured through cross-referencing multiple databases and validation by trained personnel.

Patient information, including age, age at diagnosis, disease duration, initial symptoms, the Expanded Disability Status Scale (EDSS), gender, education, employment status, ethnicity, marital status, and MS subtype, was extracted from the Khuzestan MS Center database. This database covered records from the association's establishment in 2005 to the end of 2019. Prevalence estimates of MS by county in 2019 were calculated using population data adjusted for natural population growth rates between census years.

The data regarding the number of dusty days, average daily temperature (in degrees Celsius), and daily wind speed (in meters per second) were obtained from the Environmental Protection Agency and the Meteorological Organization of Khuzestan province. These datasets were collected from the beginning to the end of 2019 for various counties within the province. Subsequently, the daily values of these variables were averaged annually for each county. Finally, the prevalence rates of MS for 2019 were correlated with the annual averages of the number of dusty days, temperature, and wind speed across the 27 counties.

3.1. Statistical Analysis

Descriptive statistics, including frequency, mean, standard deviation, median, and interquartile range (IQR), were reported. A multi-pollutant generalized additive model (GAM) was employed to estimate the relative risk (RR) of the effects of dusty days, temperature, and wind speed on the prevalence of MS. Generalized additive models represent a powerful class of models for capturing nonlinear effects of continuous quantitative variables in regression models where the response follows an exponential family distribution (19). These models are an extension of generalized linear models (GLMs) and offer a high degree of flexibility. Due to their ability to model nonlinear trends, GAMs have been widely used in many studies (20-22). In fact, GAMs allow the integration of variables in a non-parametric manner using smooth functions, such as splines (19).

The assumptions for using GAMs are as follows:

1. Independence of observations.

2. The response variable should follow a distribution from the exponential family (e.g., normal, binomial, Poisson).

3. The effects of the predictors are assumed to be additive and are modeled through separate smooth functions for each predictor.

4. The smoothing functions (e.g., splines) used in the model should be smooth.

5. The link function that connects the predictors to the response variable must be appropriate (23).

The standard formulation of a GAM (24) can be stated as follows:

gEY=β0+f1X1+f2X2++fnXn

- Y is the response variable.

- E[Y] is the expected value of the response variable.

- g is the link function that relates the expected value of Y to the predictors (for example, normal, binomial, or Poisson distributions).

- β0 is the intercept.

- fi (Xi) are smooth functions applied to the predictor variables Xi (e.g., splines).

- n represents the number of predictor variables.

In this ecological study, the response variable consisted of count data (the number of MS cases per 100,000 people), thus following a Poisson distribution. However, due to the unequal mean and variance in the response variable data, the quasi-Poisson regression model was utilized.

The GAM model provides the capability to detect nonlinear relationships between variables. In this study, a semiparametric regression approach was employed, which models the effects of the predictors using a combination of parametric and nonparametric regression techniques while retaining essential data characteristics. This method allows for more precise calculation of coefficients and the effects of independent variables, resulting in a better fit to the data through the use of nonlinear and flexible smoothing functions, such as splines. The GAM utilized penalized cubic regression splines to model the nonlinear relationships between predictors and the response variable (23).

The model used in this study is outlined as follows:

ηX= Bj1 ×Temperature+ Bj2×Dusty days+Bj3 ×Wind speed+ j=1qSj(Dusty days)
ηX=0.0001516 ×Temperature+0.1071761 × Dusty days - 0.0006243×Wind speed +j=1qSj(Dusty days) 

Where:

- Bj is the estimated coefficient of covariate variable j.

- Sj is a smooth function of the covariate "dusty days" (j = 1, …., q).

- q is the number of spline functions in the GAM model.

In this analysis, the unit of measurement used was the increase of one unit in the level of all independent variables, which was applied to calculate the RR for the outcome (the prevalence of MS per 100,000 people). The unconnected lines above and below the curves represent the 95% confidence interval. The 'mgcv' package (23) in R software version 3.4.3 was employed for the GAM analysis. A significance level of less than 0.05 was considered.

4. Results

In this study, a total of 2,676 patients with MS were registered at the Khuzestan MS Center from its establishment in 2005 until the end of December 2019, all of whom were included in the analysis. The mean age of the MS patients was 31.40 ± 8.94 years. The sex ratio (female to male) was 3.05, indicating that approximately two-thirds of the patients were female. The mean age at MS diagnosis was 29.19 ± 9.16 years. Notably, 75% of the patients were diagnosed before the age of 35.75. Of the patients, 30% were female, with the majority holding associate's and bachelor's degrees (35.10%), being married (62.10%), and ethnically belonging to the Lor group (34.50%). Regarding the type of MS, relapsing-remitting MS (RRMS) was the most prevalent, while only 3.5% had primary progressive MS (PPMS) and 4.9% had secondary progressive MS (SPMS) (Table 1).

Table 1.

Distribution of Demographic and Clinical Variables in Patients with MS (N = 2676) a

VariablesValues
Age (y)
Mean ± SD31.40 ± 8.94
Median ± IQR30.00 ± 11
EDSS
Mean ± SD1.857 ± 1.7846
Median ± IQR1.500 ± 1.5
Duration of MS disease (y)
Mean ± SD1.95 ± 3.658
Median ± IQR0.0 ± 2
Age-incidence of MS
Mean ± SD29.19 ± 9.16
Median ± IQR28.00 ± 12
Time between onset of symptoms and definitive diagnosis (months)
Mean ± SD14.84 ± 31.601
Median ± IQR2.00 ± 12
Gender
Female2015 (75.3)
Male661 (24.7)
The first sensory sign
No1684 (62.9)
Yes992 (37.1)
The first motor sign
No2214 (82.7)
Yes462 (17.3)
The first equilibrium sign
No2332 (87.1)
Yes344 (12.9)
The first visual sign
No1663 (62.1)
Yes1013 (37.9)
The first urinary disorder sign
No2603 (97.3)
Yes73 (2.7)
The first fatigue sign
No2550 (95.3)
Yes126 (4.7)
Education
Illiterate144 (5.4)
Under diploma679 (25.4)
Diploma782 (29.2)
Post-diploma and bachelor's degree940 (35.1)
Master’s degree and higher131 (4.9)
Job
Employed697 (26.0)
Unemployed749 (28.0)
Disabled91 (3.4)
Housewife1139 (42.6)
Race
Persian891 (33.3)
Turkish61 (2.3)
Kurd63 (2.4)
Lor922 (34.5)
Arab704 (26.3)
Others35 (1.3)
Marital status
Single931 (34.8)
Married1661 (62.1)
Divorced62 (2.3)
Widower/widow22 (0.8)
MS type
RR-MS2041 (76.3)
PP-MS93 (3.5)
SP-MS132 (4.9)
RP-MS25 (0.9)
DEVIC67 (2.5)
CIS318 (11.9)

Overall, in 2019, the average annual prevalence of MS in Khuzestan province was 42 per 100,000 people. Table 2 presents the distribution of MS prevalence according to geographical classification in Khuzestan province. As shown, the lowest prevalence of MS was observed in Bavi county, with a rate of 4 cases per 100,000 individuals. Conversely, the highest prevalence of MS occurred in Behbahan county, with a rate of 112 cases per 100,000 individuals, followed by Ahvaz county with a rate of 109 cases per 100,000 individuals in 2019 (Table 2, Figure 1). Additionally, the mean temperature, number of dusty days, and wind speed (expressed in kilometers per hour) across Khuzestan province in 2019 were recorded as 25.25 ± 2.16 (ranging from 20 to 30), 47.26 ± 15.7 (ranging from 18 to 78), and 39.37 ± 13.13 (ranging from 14.30 to 60.80), respectively.

Table 2.

Distribution of Annual Average Temperature and Wind Speed, Frequency of Dusty Days, and Prevalence of Multiple Sclerosis per 100,000 Individuals by County in Khuzestan Province, 2019

Names MS Prevalence per 100,000 IndividualsTemperature (°C)Dusty DaysWind Speed (km/h)
Bavi4247824.2
Karun1225.26630.1
Haftgel14274648.5
Andika15274648.5
Hamidiyeh17232556.74
Shadegan17255046.5
Gotvand20202443.8
Shush20224632.1
Ramshir20256030.1
Dasht-e Azadegan22283624.9
Dezful29241860.8
Izeh33247824.2
Hendijan34242556.74
Andimeshk3423.93746.7
Baghmalek39305224.9
Masjedsoleyman41263445.3
Lali4227.55224.9
Bandar-e Mahshahr5427.43043.5
Abadan5527.55955.9
Khorramshahr5524.56057.3
Omidiyeh55254946.5
Shushtar66254946.5
Ramhormoz6624.54943.8
Aghajari6828.43614.3
Hoveyzeh7722.75421.7
Ahvaz10925.26630.1
Behbahan112265134.5
Distribution of multiple sclerosis (MS) prevalence (per 100,000 individuals) in Khuzestan province by geographical classification
Distribution of multiple sclerosis (MS) prevalence (per 100,000 individuals) in Khuzestan province by geographical classification

The results of the GAM model indicated that there was no statistically significant association between average temperature and the prevalence of MS (P-value = 0.994). However, a significant positive statistical association was observed between the number of dusty days and the prevalence of MS in Khuzestan province (P-value < 0.001). For each unit increase in the number of dusty days, the prevalence of MS increased by 11 cases per 100,000 individuals (RR = 1.113, P-value < 0.001). Conversely, a significant negative statistical association was observed between average wind speed and the prevalence of MS. For each unit increase in average wind speed, the prevalence of MS decreased by 1 case per 100,000 individuals (RR = 0.999, P-value < 0.001) (Table 3).

Table 3.

The Association Between Temperature, Dusty Days, and Wind Speed with the Prevalence of Multiple Sclerosis in Khuzestan Province, 2019

VariablesCoefficientStandard ErrorRR95 % CIP-Value
Temperature (°C)0.0000.0201.0000.961 - 1.0410.994
Dusty days0.1070.0111.1131.090 - 1.137< 0.001 a
Wind speed (km/h)-0.0010.0000.9990.9990 - 0.9997< 0.001 a

5. Discussion

Recent studies in Iran have investigated the role of air pollution as a significant factor in the occurrence of autoimmune diseases, including MS (7, 25, 26). A study conducted between 2008 and 2016 in Isfahan found that the level of the Air Quality Index (AQI) was associated with the EDSS of MS patients, suggesting that air pollution may play a role in both the severity of MS (27).

Outside of Iran, Jeanjean et al. (28), in Strasbourg, France, conducted a study between 2000 and 2009 based on meteorological parameters, school holidays, and public holidays. They reported that the impact of PM10 on the onset and recurrence of MS was more significant in the winter season than in the summer season, which is consistent with our study's findings that an increase in the number of dusty days is associated with an increase in the prevalence of MS in Khuzestan province. Similar positive associations between particulate matter air pollution and MS have been reported in other studies conducted in Iran, Italy, and Finland (12, 13, 29).

In a cross-sectional study by Tateo et al. in Padua, Italy, the results indicated a strong association between the prevalence of MS and annual PM2.5 concentration levels (r = 0.81, P < 0.001). Regression analysis further confirmed this relationship, revealing that MS cases were significantly associated with average PM2.5 levels (β = 0.11, P < 0.001) (30). In the United States, a case-control study by Lavery et al. showed that particulate matter (PM2.5) (OR = 3.96, 95% CI: 1.42 - 11.1, P < 0.01), along with three other criteria pollutants (SO₂, CO, and lead), was significantly associated with an increased likelihood of pediatric MS (31).

A case-crossover study by Jeanjean et al. in Paris, France, indicated that elevated NO2 and PM10 exposure were linked to an increased risk of MS relapse during the cold season (OR = 1.08, 95% CI: 1.03 - 1.14; OR = 1.06, 95% CI: 1.01 - 1.11, respectively). Moreover, a significant association was found between ozone (O3) exposure and MS relapse risk during the hot season (OR = 1.16, 95% CI: 1.07 - 1.25) (28). Similarly, a cross-sectional study by Ashtari et al. in Isfahan, Iran, revealed that the AQI level correlated with the degree of complete remission following the first attack of MS, with an odds ratio of 1.005 (95% CI: 1.004 - 1.006) (27). In an ecological study by Heydarpour et al. conducted from 2003 to 2013 in Tehran, Iran, the findings indicated that MS cases exhibited a clustered distribution throughout the city. Additionally, a significant difference in exposure to pollutants such as PM10, SO2, NO2, and NOx (P < 0.001) was found when comparing MS cases to control subjects (12).

Collectively, these studies underscore the significant association between particulate matter and MS, suggesting that exposure to air pollution, particularly particulate matter, may contribute to the prevalence and severity of this neurological disease.

Among the environmental factors considered, some studies have focused on the role of temperature, sunlight, and other meteorological parameters. In this regard, Vojinovic et al., in the city of Nis, Serbia, conducted a 5-year study and reported a significant positive correlation between total cloudiness and daily cloudiness with the number of MS relapses, although no relationship was found between disease activity and direct sunlight exposure in terms of hours, days, or sunny months (32). In another study, Tremlett in southern Tasmania, Australia, observed a lower frequency of MS relapses in summer compared to winter, based on data collected from January 2002 to April 2005 (33). Similarly, Ma and Zhang in China, from January 2002 to December 2012, associated MS onset/relapse with geographical latitude, temperature, and sunlight intensity. They observed the highest onset/relapse of MS in winter (134 cases) and the lowest in summer (97 cases). According to them, environmental factors, especially sun exposure, had the most significant impact on disease onset and relapse, while gender, age of onset, and disease duration were not significant (34). On the other hand, in a study by Abella-Corral et al. conducted from 1997 to 2002, a significant difference in the prevalence of MS was reported between summer (high prevalence, peaking in June) and winter (low prevalence, lowest in December) (35).

Some studies have examined the seasonal and monthly correlation between meteorological parameters and the prevalence and relapse of MS. For instance, Fundora-Hernandez et al. in Cuba, between April 2004 and November 2007, reported the highest frequency of MS prevalence in the second quarter of the year, from April to June, attributing it to the transition from winter to summer and the variable thermal oscillation during that period (36). Additionally, Ogawa et al. in Japan, in 2004, investigated the correlation between temperature and MS exacerbation. They divided the 12 months of the year into 6 groups based on the monthly average temperature and found that the total number of disease attacks was higher in the warmest months (July and August) and the coldest months (January and February) compared to other months (37).

Recently, in a population-based study involving three countries with high MS prevalence (Canada, Italy, and Norway), behavior and lifestyle related to sunlight exposure in summer and winter were examined. The results indicated that the risk of MS was higher among individuals who spent more time indoors during their childhood and protected themselves from the sun during the limited periods they spent outdoors. However, this protection did not affect the risk of MS among individuals who spent a considerable amount of time outdoors in the summer and a moderate amount of time outdoors in the winter (38). These findings contrast with the results of our study regarding the role of temperature in the prevalence of MS. Our study did not observe a significant statistical correlation between average temperature and MS prevalence. Comparing our findings with those of the aforementioned studies is somewhat challenging, as they used different meteorological criteria related to the prevalence or exacerbation of MS, likely due to the substantial differences in climate between Khuzestan province and the other regions studied.

A significant limitation of this study is its reliance on aggregated data, which restricts the applicability of the findings to individual cases. Additionally, the exposure levels captured by air pollution and meteorological stations may not accurately reflect personal exposure. Since the study was ecological and population-based, it was not possible to account for individual confounders such as smoking, education, and other relevant variables. One of the strengths of this study is that, at the time it was conducted, no prior research had explored the connection between exposure to dusty days, temperature, wind speed, and the prevalence of MS in Khuzestan province, Iran.

Our study results showed that an increase in dusty days is associated with a higher prevalence of MS, highlighting the need for further investigation into the environmental factors affecting MS. These findings also suggest that public health interventions in Khuzestan province should prioritize monitoring and managing environmental factors, particularly during dusty seasons. Finally, future research should explore the mechanisms by which environmental factors, especially particulate matter, influence MS risk and consider longitudinal studies to assess changes over time.

Acknowledgements

References

  • 1.

    Ascherio A. Environmental factors in multiple sclerosis. Expert Rev Neurother. 2013;13(12 Suppl):3-9. [PubMed ID: 24289836]. https://doi.org/10.1586/14737175.2013.865866.

  • 2.

    Hirtz D, Thurman DJ, Gwinn-Hardy K, Mohamed M, Chaudhuri AR, Zalutsky R. How common are the "common" neurologic disorders? Neurology. 2007;68(5):326-37. [PubMed ID: 17261678]. https://doi.org/10.1212/01.wnl.0000252807.38124.a3.

  • 3.

    Papathanasopoulos P, Gourzoulidou E, Messinis L, Georgiou V, Leotsinidis M. Prevalence and incidence of multiple sclerosis in western Greece: a 23-year survey. Neuroepidemiology. 2008;30(3):167-73. [PubMed ID: 18382116]. https://doi.org/10.1159/000122334.

  • 4.

    Azami M, YektaKooshali MH, Shohani M, Khorshidi A, Mahmudi L. Epidemiology of multiple sclerosis in Iran: A systematic review and meta-analysis. PLoS One. 2019;14(4). e0214738. [PubMed ID: 30964886]. [PubMed Central ID: PMC6456231]. https://doi.org/10.1371/journal.pone.0214738.

  • 5.

    Sahraian MA, Sahebkar M, Dehghani R, Derakhshan-Jazari M, Kazami-Moghaddam V, Kouchaki E. Multiple sclerosis-A disease on a dramatically rising trend in Iran: Review of possible reasons. Iran J Neurol. 2017;16(1):34-40. [PubMed ID: 28717432]. [PubMed Central ID: PMC5506754].

  • 6.

    Dehghani R, Kazemi Moghaddam V. [Investigation of the possible causes of the increased prevalence of multiple sclerosis in Iran: a review]. Pars J Med Sci. 2015;13(2):17-25. FA. https://doi.org/10.29252/jmj.13.2.3.

  • 7.

    Dehghani R, Yunesian M, Sahraian MA, Gilasi HR, Kazemi Moghaddam V. The Evaluation of Multiple Sclerosis Dispersal in Iran and Its Association with Urbanization, Life Style and Industry. Iran J Public Health. 2015;44(6):830-8. [PubMed ID: 26258096]. [PubMed Central ID: PMC4524308].

  • 8.

    Ebrahimi HA, Sedighi B. Prevalence of multiple sclerosis and environmental factors in Kerman province, Iran. Neurol Asia. 2013;18(4):385-9.

  • 9.

    Spelman T, Gray O, Trojano M, Petersen T, Izquierdo G, Lugaresi A, et al. Seasonal variation of relapse rate in multiple sclerosis is latitude dependent. Ann Neurol. 2014;76(6):880-90. [PubMed ID: 25283272]. https://doi.org/10.1002/ana.24287.

  • 10.

    Salvi F, Bartolomei I, Smolensky MH, Lorusso A, Barbarossa E, Malagoni AM, et al. A seasonal periodicity in relapses of multiple sclerosis? A single-center, population-based, preliminary study conducted in Bologna, Italy. BMC Neurol. 2010;10:105. [PubMed ID: 21040535]. [PubMed Central ID: PMC2988761]. https://doi.org/10.1186/1471-2377-10-105.

  • 11.

    Saaroni H, Sigal A, Lejbkowicz I, Miller A. Mediterranean weather conditions and exacerbations of multiple sclerosis. Neuroepidemiology. 2010;35(2):142-51. [PubMed ID: 20571282]. https://doi.org/10.1159/000313446.

  • 12.

    Heydarpour P, Amini H, Khoshkish S, Seidkhani H, Sahraian MA, Yunesian M. Potential impact of air pollution on multiple sclerosis in Tehran, Iran. Neuroepidemiology. 2014;43(3-4):233-8. [PubMed ID: 25501708]. https://doi.org/10.1159/000368553.

  • 13.

    Angelici L, Piola M, Cavalleri T, Randi G, Cortini F, Bergamaschi R, et al. Effects of particulate matter exposure on multiple sclerosis hospital admission in Lombardy region, Italy. Environ Res. 2016;145:68-73. [PubMed ID: 26624240]. [PubMed Central ID: PMC4797990]. https://doi.org/10.1016/j.envres.2015.11.017.

  • 14.

    Dastoorpoor M, Sekhavatpour Z, Masoumi K, Mohammadi MJ, Aghababaeian H, Khanjani N, et al. Air pollution and hospital admissions for cardiovascular diseases in Ahvaz, Iran. Sci Total Environ. 2019;652:1318-30. [PubMed ID: 30586817]. https://doi.org/10.1016/j.scitotenv.2018.10.285.

  • 15.

    Masoumi K, Haddadzadeh Shoushtari M, Forouzan A, Asgari Darian A, Dastoorpoor M, Ebrahimzadeh P, et al. Rainfall-Associated Bronchospasm Epidemics: The Epidemiological Effects of Air Pollutants and Weather Variables. Can Respir J. 2017;2017:9252069. [PubMed ID: 29089817]. [PubMed Central ID: PMC5635285]. https://doi.org/10.1155/2017/9252069.

  • 16.

    Dastoorpoor M, Idani E, Goudarzi G, Khanjani N. Acute effects of air pollution on spontaneous abortion, premature delivery, and stillbirth in Ahvaz, Iran: a time-series study. Environ Sci Pollut Res Int. 2018;25(6):5447-58. [PubMed ID: 29214476]. https://doi.org/10.1007/s11356-017-0692-9.

  • 17.

    Olek M. Multiple Sclerosis: Etiology, Diagnosis, and New. Treatment Strategies. Springer: Berlin, Germany; 2007.

  • 18.

    Dastoorpoor M, Khodadadi N, Majdinasab N, Soltani M, Ghanbari S. The Time-trend of Multiple Sclerosis Incidence and Prevalence in Khuzestan Province, Iran. Caspian J Neurol Sci. 2022;8(2):104-16. https://doi.org/10.32598/CJNS.8.29.4.

  • 19.

    Terzi Y, Cengiz MA. Using of generalized additive model for model selection in multiple poisson regression for air pollution data. Sci Res Essays. 2009;4(9):867-71.

  • 20.

    Dominici F, McDermott A, Zeger SL, Samet JM. On the use of generalized additive models in time-series studies of air pollution and health. Am J Epidemiol. 2002;156(3):193-203. [PubMed ID: 12142253]. https://doi.org/10.1093/aje/kwf062.

  • 21.

    Moradgholi A, Ahmadi Angali K, Akhoond MR, Goudarzi GH, Dastoorpour M. [The Severity of the Relationship Between Daily air Pollution and Cardiovascular Deaths in Ahvaz, Iran- Using Generalized Additive Models(GAMs) for Seven Years During March 2008 - March 2015]. Armaghan-e-Danesh. 2019;24(4):638-54. FA. https://doi.org/10.52547/armaghanj.24.4.638.

  • 22.

    Dastoorpoor M, Khanjani N, Moradgholi A, Sarizadeh R, Cheraghi M, Estebsari F. Prenatal exposure to ambient air pollution and adverse pregnancy outcomes in Ahvaz, Iran: a generalized additive model. Int Arch Occup Environ Health. 2021;94(2):309-24. [PubMed ID: 32936369]. https://doi.org/10.1007/s00420-020-01577-8.

  • 23.

    Wood SN. Generalized Additive Models. New York, US: Taylor & Francis; 2017. https://doi.org/10.1201/9781315370279.

  • 24.

    Hastie TJ. Generalized Additive Models. In: Chambers JM, Hastie TJ, editors. Statistical Models in S. England, UK: Routledge; 2017. p. 249-307. https://doi.org/10.1201/9780203738535-7.

  • 25.

    Vojdani A. A Potential Link between Environmental Triggers and Autoimmunity. Autoimmune Dis. 2014;2014:437231. [PubMed ID: 24688790]. [PubMed Central ID: PMC3945069]. https://doi.org/10.1155/2014/437231.

  • 26.

    Esmaeil Mousavi S, Heydarpour P, Reis J, Amiri M, Sahraian MA. Multiple sclerosis and air pollution exposure: Mechanisms toward brain autoimmunity. Med Hypotheses. 2017;100:23-30. [PubMed ID: 28236843]. https://doi.org/10.1016/j.mehy.2017.01.003.

  • 27.

    Ashtari F, Esmaeil N, Mansourian M, Poursafa P, Mirmosayyeb O, Barzegar M, et al. An 8-year study of people with multiple sclerosis in Isfahan, Iran: Association between environmental air pollutants and severity of disease. J Neuroimmunol. 2018;319:106-11. [PubMed ID: 29526408]. https://doi.org/10.1016/j.jneuroim.2018.02.019.

  • 28.

    Jeanjean M, Bind MA, Roux J, Ongagna JC, de Seze J, Bard D, et al. Ozone, NO(2) and PM(10) are associated with the occurrence of multiple sclerosis relapses. Evidence from seasonal multi-pollutant analyses. Environ Res. 2018;163:43-52. [PubMed ID: 29426027]. [PubMed Central ID: PMC5886008]. https://doi.org/10.1016/j.envres.2018.01.040.

  • 29.

    Oikonen M, Laaksonen M, Laippala P, Oksaranta O, Lilius EM, Lindgren S, et al. Ambient air quality and occurrence of multiple sclerosis relapse. Neuroepidemiology. 2003;22(1):95-9. [PubMed ID: 12566960]. https://doi.org/10.1159/000067108.

  • 30.

    Tateo F, Grassivaro F, Ermani M, Puthenparampil M, Gallo P. PM2.5 levels strongly associate with multiple sclerosis prevalence in the Province of Padua, Veneto Region, North-East Italy. Mult Scler. 2019;25(13):1719-27. [PubMed ID: 30270719]. https://doi.org/10.1177/1352458518803273.

  • 31.

    Lavery AM, Waubant E, Casper TC, Roalstad S, Candee M, Rose J, et al. Urban air quality and associations with pediatric multiple sclerosis. Ann Clin Transl Neurol. 2018;5(10):1146-53. [PubMed ID: 30349849]. [PubMed Central ID: PMC6186930]. https://doi.org/10.1002/acn3.616.

  • 32.

    Vojinovic S, Savic D, Lukic S, Savic L, Vojinovic J. Disease relapses in multiple sclerosis can be influenced by air pollution and climate seasonal conditions. Vojnosanit Pregl. 2015;72(1):44-9. [PubMed ID: 26043590]. https://doi.org/10.2298/vsp140121030v.

  • 33.

    Tremlett H, van der Mei IA, Pittas F, Blizzard L, Paley G, Mesaros D, et al. Monthly ambient sunlight, infections and relapse rates in multiple sclerosis. Neuroepidemiology. 2008;31(4):271-9. [PubMed ID: 18971584]. https://doi.org/10.1159/000166602.

  • 34.

    Ma J, Zhang X. [The relationship between season/latitude and multiple sclerosis]. Zhonghua Nei Ke Za Zhi. 2015;54(11):945-8. ZH. [PubMed ID: 26759213].

  • 35.

    Abella-Corral J, Prieto JM, Dapena-Bolano D, Iglesias-Gomez S, Noya-Garcia M, Lema M. [Seasonal variations in the outbreaks in patients with multiple sclerosis]. Rev Neurol. 2005;40(7):394-6. ES. [PubMed ID: 15849671].

  • 36.

    Fundora-Hernandez H, Dorta-Contreras AJ, Socias-Lopez M, Fraga-Santana S, Villatoro-Rodriguez SN, Padilla-Docal B, et al. Seasonal distribution and evolving forms of multiple sclerosis patients diagnosed from April 2004 to November 2007. Arq Neuropsiquiatr. 2009;67(3A):661-3. [PubMed ID: 19722045]. https://doi.org/10.1590/s0004-282x2009000400016.

  • 37.

    Ogawa G, Mochizuki H, Kanzaki M, Kaida K, Motoyoshi K, Kamakura K. Seasonal variation of multiple sclerosis exacerbations in Japan. Neurol Sci. 2004;24(6):417-9. [PubMed ID: 14767690]. https://doi.org/10.1007/s10072-003-0200-5.

  • 38.

    Magalhaes S, Pugliatti M, Riise T, Myhr KM, Ciampi A, Bjornevik K, et al. Shedding light on the link between early life sun exposure and risk of multiple sclerosis: results from the EnvIMS Study. Int J Epidemiol. 2019;48(4):1073-82. [PubMed ID: 30561654]. [PubMed Central ID: PMC6693814]. https://doi.org/10.1093/ije/dyy269.