2.1. Study Design and Study Population:
In this retrospective (historical) cohort study, the study population was all people exposed to mustard gas in sardasht in 1987 based on the approval of the VMAF (the Veterans and Martyr Affair Foundation). VMAF is the legal authority in Iran that provides pension and health services to war survivor. This foundation registered all survivors of war and then regularly provides long-term care services to them and their families by maintaining excellence in health care services through professional collaboration [6, 8]. Therefore, the inclusion criterion was the approval of the Veterans and Martyr Affair Foundation of Iran to contact of the persons with the chemical agent of mustard gas and the exclusion criterion was the occurrence of premature death. According to the contract in this study, premature death refers to death at the time of exposure (during the first 14 days after exposure).
2.2. Data Collection:
Data related to health status (death / life) and other information related to the subjects were obtained by referring to the Veterans and Martyr Affair Foundation of Iran and reviewing the VMAF database up to March 20, 2019. Normally, all information related to the deaths of war survivors in Iran is transferred to the Veterans and Martyr Affair Foundation of Iran, where a death commission consisting of various specialists and experts investigates the causes of death. The members of this committee are aware of the possible chronic effects of mustard gas and determine the causes of death based on the medical record and death certificate.
In addition to data on the health status (death / life) of the subjects, other information included gender (male / female), age of the person at the time of exposure, severity of late complication (no symptom, mild symptom, moderate symptom, severe symptom), organ of the late complication (lung, Eyes, skin) beside evacuation and hospital admission status (EA, NEA, None) were also taken from the mentioned database and recorded. In order to classified subject according to evacuation and hospital admission, based on the evidence of exposure to MG, the data were categorized into 3 groups: (1) evacuated and admitted (EA), confirmed history of exposure in the affected geographic region with evacuation and hospital admission; (2) not evacuated or admitted (NEA), confirmed history of exposure without evacuation or hospital admission; and (3) undocumented cases, not classified in the other 2 groups but showing late clinical manifestations suggestive of exposure. The classification criterion for the severity of the late complication which periodically examined and record by a specialized medical team was based on the protocol compiled by a panel of experts at the Veterans and Martyr Affair Foundation of Iran. Based on this protocol, the complications in these people are divided into four categories: no symptom, mild, moderate and severe symptom. The panel included pulmonologists, psychiatrists, neurologists, ophthalmologists, dermatologists, cardiologists, forensic specialists and nephrologists, reviewed the relevant documents and evidence and categorized them [6]. Definition of late complications in lung, skin and eye were as follow:
Lungs:
Categorization of complication in lung following MG exposure was defined using spirometry results of forced expiratory volume in the first second (FEV1) as follow: no symptoms, FEV1 higher than 80%; mild, FEV1 higher than 70% but lower than or equal to 80%; moderate, FEV1 higher than 50% but lower than or equal to 70%; and severe, FEV1 lower than or equal to 50% with normal arterial blood gas or severe tracheabronchomalacia [6, 9].
Skin:
Category of skin complications subsequent exposure to MG was determined as follows: Mild category was defined as the presence of an MG-specific scar (pigmentation, vascular, or trophic changes) across less than 5%of the body surface; the presence of such a scar plus mild dermatitis and mild xerosis or mild-to-moderate autoimmune diseases; and concurrent medical documents plus severe autoimmune diseases or any mild dermatitis. The moderate category was defined as the presence of an MG-specific scar affecting between 5%and 20%of the body surface2; presence of an MG-specific scar with significant trophic changes in the genital or anus areas; the presence of an MG-specific scar of any size with or without concurrent medical documents plus severe dermatitis covering less than 50% of the body surface; concurrent medical documents with severe dermatitis or xerosis covering less than 50% of the body surface; and the presence of an MG-specific scar of any size with or without concurrent medical documents plus any severe autoimmune disease affecting the skin covering more than 50% of the body surface. The severe category was defined as the presence of an MG-specific scar covering more than 20%of the body surface; the presence of an MG-specific scar of any size with or without concurrent medical documents with severe dermatitis or xerosis covering more than 50% of the body surface; and concurrent medical documents with severe dermatitis or xerosis covering less than 50% of the body surface [6, 10].
Eyes:
Complication in eye following MG exposure was categorized as follow: Mild category was defined as burning, itching, tearing, redness, foreign body sensation, blurred vision, conjunctivitis, subconjunctival hemorrhage, photophobia, or conjunctival vascular changes (telangiectasia and vascular changes, edema of the eyelid, papillary changes) of 1 or both eyes. The moderate category was defined as the same changes in the mild category along with conjunctival ischemia in the limbus area, as well as unilateral or bilateral closure or cauterization of punctum. The severe category was defined as those observations in the moderate category along with symptoms of corneal involvement, epithelial and subepithelial opacity and anterior stroma in the cornea, keratopathy, pannus, hyperpigmentation around the nerve, iron deposition in the cornea and corneal vascularization, stenosis, and involvement of less than half unilaterally or bilaterally in both eyes [6, 11].
2.3. Statistical Analysis:
After collecting data, all statistical analyzes were performed in Stata software version 15. At first, the frequency and relative frequency index was used to express the basic characteristics of the subjects.
In order to describe the mortality in the cohort, we used cumulative mortality rate and mortality rate index. To calculate mortality rate index we first calculated the person-years at risk in these individuals from the time of exposure to mustard gas to the occurrence of death in the case of subject who died, or at the end of follow up for other subject. Then all-cause mortality rate and the cause-specific mortality rate were calculated.
In order to compare the mortality in the cohort with the general population (Iran), we used standardized mortality ratio (SMR) index. SMR is the number of observed deaths divided by the expected death. To calculate expected death, first, mortality rate of the general population according to sex, age (5 year age group) and calendar year were obtained using the information published by GBD 2019 Demographics Collaborators [12]. Then expected death was determined by multiplying the mortality rate of the general population and the person-years of the cohort. Finally we calculated all-cause SMR and cause-specific SMR. Also, all-cause SMR was determined according to age, sex and severity of late complication. For each SMR, a 95% confidence interval was calculated according to the assumption that the number of deaths observed follows the Poisson distribution.
Eventually, cox regression model was used to determine the relationship between the independent variables (age at time of exposure, sex, severity of late complication, evacuation and hospital admission Status) and risk of death. To do this, first, the univariate Cox regression model was used and the variables that had a significance level of less than 0.1 in this model were included in the multivariate Cox regression model. Finally, variables that had a significance level of less than 0.05 in the multiple models were considered as a factor affecting the risk of death in the cohort.