-
The basic clinical and exposure data of the study participants (all male) are shown in Table 1. The statistical findings showed that there is no significant difference in age, BMI, daily exposure time, job history, and smoking history between the case and control groups. Regarding the case group, they had a mean age of 34.26 years (with a range of 21–58 years) and a mean exposure period of 11.8 hours (range of 8–16 h) over a job span of 7.38 years (range of 2–20 years). Moreover, their average BMI was 25.85 kg/m2 (range of 16.48–35.83 kg/m2) and 5.3 % of them had a history of smoking. The above information for the control group was as follows; a mean age of 34.47 years (with a range of 24–51 years), a mean exposure period of 11.66 hours (range of 4–15 h) over a job span of 7.84 years (range of 2–24), a mean BMI of 25.85 kg/m2 (range of 18.6–33.9 kg/m2), and 5.8% of them were smokers.
Table 1. Characteristics of the BTEX and other variables in two groups of the study participants.
Variables Case (n=138) Control (n=145) P-values Age (years), mean ± SD 34.26 ± 6.50 34.47 ± 5.36 0.67 BMI (kg/m2), mean ± SD 25.85 ± 3.38 25.40 ± 4.78 0.84 Length of exposure in a day (hours), mean ± SD 11.80 ± 1.44 11.66 ± 1.45 0.45 Job history (year), mean ± SD 7.38 ± 3.75 7.84 ± 4.16 0.47 Smoking, n (%) 0.87 Yes 7 (5.3) 8 (5.8) NO 125 (94.7) 131 (94.2) BTEX-exposure measurements (ppm), mean ± SD Benzene 0.57 ± 1.51 0.71 ± 1.59 0.99 Toluene 1.02 ± 2.03 1.11 ± 2.17 0.88 Ethylbenzene 1.44 ± 2.52 1.24 ± 2.70 0.27 O-xylene 0.437 ± 1.33 0.414 ± 1.64 0.27 P-xylene 2.47 ± 3.06 2.28 ± 3.04 0.49 M xylene 1.70 ± 2.46 1.79 ± 2.64 0.89 Note. Student’s t-test was used to determine the statistical significance of the differences in age, BMI, daily exposure time, and job history, and the chi-square test was used for smoking history; P-value < 0.05 was considered statistically significant (in bold). BMI, body mass index; BTEX, benzene toluene ethylbenzene xylene. *Some data were missing due to lack of information. -
As mentioned in the previously, in this study we measured the concentration of each of the BTEX components in the breathing air environment of the study participants to evaluate the exposure. The results of this assessment showed that the mean of arithmetic concentrations of benzene, toluene, and ethylbenzene in the air environment of the population was 0.70 (standard deviation = 1.67 ppm), 1.10 (standard deviation = 2.14 ppm), and 1.34 (standard deviation = 2.6 ppm), respectively. Moreover, the composition of xylene is considered as a set mixture of ortho, para, and meta isomers whose average concentrations were 0.44 (standard deviation = 1.48 ppm), 2.36 (standard deviation = 3.04), and 1.75 (standard deviation = 2.54 ppm), respectively. The findings also showed that there was no significant difference in the concentrations of various BTEX compounds in the breathing environments of both patients and healthy individuals (Table 1).
-
As shown in Table 2, a significant difference was observed in the levels of some blood indices between the two study groups. The mean RDW and lymphocyte count were significantly higher in patients than in healthy controls (P = 0.018 and P = 0.001, respectively). The mean HB, MCV, MCH, MCHC, and granulocyte counts revealed a significant decrease in the case group compared to the control group (P = 0.001, P = 0.004, P = 0.001, P = 0.002, and P = 0.001, respectively). Other blood indices, including HCT, WBC, RBC, PLT count, and monocyte count, were not significantly different between the two groups. Additionally, regarding the RBC-related indices, including RBC count, HB, MCV, HCT, MCH, mean MCHC, and red RDW, it was found that their abnormalities in patients was 15.21%, 22.46%, 13.76%, 21.01%, 25.36%, 23.18%, and 9.42%, respectively. Moreover, 7.20% of the case group showed abnormal WBC index. The abnormality rates for lymphocytes, monocytes, and granulocytes was 62.31%, 10.86%, and 8.69%, respectively. For the platelet index, we observed a 10.14% aberration (Supplementary Tables, available in www.besjournal.com).
Table 2. Data on hematological indices in case and control subjects
Hematological indices Case (n = 138) Control (n=145) P-values WBC (×109/L) 6.97 ± 1.86 7.03 ± 1.51 0.80 RBC (×109/L) 5.04 ± 0.559 5. 0 ± 0.360 0.47 HB (g/dL) 14.77 ± 1.50 15.32 ± 0.840 0.001 HCT (%) 43.52 ± 3.65 44.22 ± 2.24 0.051 MCV (fL) 86.20 ± 9.30 88.67 ± 3.82 0.004 MCH (Pg) 29.64 ± 3.67 30.75 ± 1.71 0.001 MCHC (g/dL) 34.10 ± 1.89 34.66 ± 0.911 0.002 RDW (%) 12.41 ± 0.968 12.13 ± 0.562 0.018 Platelets (×109/L) 212 ± 46.64 210.8 ± 38.52 0.908 Lymphocytes (%) 44.72 ± 10.87 37.96 ± 5.72 0.001 Monocytes (%) 2.92 ± 1.59 3.08 ± 1.79 0.936 Granulocytes (%) 51.07 ± 10.68 57.97 ± 6.38 0.001 Note. Values are given as the mean ± standard deviation; Student t-test for comparing the difference of hematological indices in the case and control groups was used; P-value <0.05 was considered statistically significant (in bold). WBC, white blood cell count; RBC, red blood cell count; HB, hemoglobin; HCT, hematocrit; MCV, mean corpuscular volume; MCH, mean corpuscular hemoglobin; MCHC, mean corpuscular hemoglobin concentration; RDW, red cell distribution width. -
As described in Table 3, in the current population, the genotype frequencies of miR-100-rs543412 C>T in the cases and controls were consistent with HWE (P > 0.05), whereas the genotypes of miR-506-rs5905008 A>G were not in equilibrium (P > 0.05). Regarding the comparison of allelic and genotype frequencies of miR-100-rs543412 C>T polymorphism between the two groups, we observed that patients differ from healthy individuals in both allelic (P = 0.005) and genotype (P = 0.013) frequencies, in which the patient group carried significantly fewer T minor allele (39.5% vs. 51.4 %, P = 0.005) and TT minor genotype (13.8% vs. 24.1%, P = 0.013) compared with the healthy group.
Table 3. Genotype and allele distribution of the selected polymorphisms in study population groups.
SNPs Cases (n=138),
N (%)Controls (n=145)
N (%)P-Value P-HWE miR-100-rs543412c>T Genotypes
CC
CT
TT
48(34.8)
71(51.4)
19(13.8)
31(21.4)
79(54.5)
35(24.1)
0.013
0.24Alleles
C
T
167(60.5)
109(39.5)
141(48.6)
149(51.4)0.005 miR-506-rs5905008a>G
Genotypes
AA
AG
GG
21(15.2)
17(12.3)
100(72.5)
23(15.9)
23(15.9)
99(68.3)
0.66
0.001Alleles
A
G
217(78.6)
59(21.4)
221(76.2)
69(23.8)
0.49Note. Person χ2 test used for difference in distributions between the case and control groups. A goodness-of-fit chi-squared test was used to evaluate the Hardy-Weinberg equilibrium in the study. population; P-value < 0.05 was considered statistically significant (in bold). SNPs, single nucleotide polymorphisms; P-HWE, P value for Hardy-Weinberg equilibrium. We also used multiple genetic models of inheritance by regression analysis to determine the risk or protective relationship of each allele and genotype of this SNP with abnormalities in blood indices in the exposed subjects (Table 4). It was observed that TT and CT genotypes of the rs543412 C > T was significantly associated with the decreased risk of abnormal blood indices under genetic models of codominant (CT vs. CC, OR: 0.546, 95% CI: 0.312–0.957, P = 0.034 & TT vs. CC, OR: 0.335, 95% CI: 0.158–0.711, P = 0.004), dominant (CT+TT vs. CC, OR: 0.486, 95% CI: 0.285–0.830, P = 0.008), and recessive (TT vs. CC+CT, OR: 0.507, 95% CI: 0.267–0.962, P = 0.038). It should be noted that this finding was maintained in both the adjusted (for variables of age, BMI, smoking, and length of exposure to BTEXs) and unadjusted states, except for the CT heterozygous genotype in the codominant model, which was not significant, with a borderline value (P = 0.054). Finally, we conducted further stratification analyses to evaluate the independent effect of each genotype on each hematological parameter in BTEX-exposed workers, and the findings were not significant for any subgroup (Supplementary Tables).
Table 4. Association results between miR-SNP genotypes and patterns of hematological parameters in multiple inheritance models.
SNPs Genotype models Genotypes Population groups Unadjusted OR (95%CI) P Adjusted OR (95%CI) P Case (n = 138) Control (n = 145) miR-100-RS543412C>T Codominant C/CC/TT/T 48 (34.8%) 31 (21.4%) 1 1 71 (51.4%) 79 (54.5%) 0.580 (0.334-1.010) 0.054 0.546 (0.312-0.957) 0.034 19 (13.8%) 35 (24.1%) 0.351 (0.171-0.719) 0.004 0.335 (0.158-0.711) 0.004 Dominant C/CC/T+T/T 48 (34.8%) 31 (21.4%) 1 1 90 (65.2%) 114 (78.6%) 0.510 (0.300-0.866) 0.013 0.486 (0.285-0.830) 0.008 Recessive C/C+C/TT/T 119 (86.2%) 110 (75. 9%) 1 1 19 (13.8%) 35 (24.1%) 0.502 (0.271-0.929) 0.028 0.507 (0.267-0.962) 0.038 Overdominant C/C+T/TC/T 67 (48.6%) 66 (45.5%) 1 1 71 (51.4%) 79 (54.4%) 0.885 (0.555-1.413) 0.609 0.822 (0.510 -1.325) 0.421 miR-506-RS5905008A>G Codominant A/AA/GG/G 21 (15.2%) 23 (15.9%) 1 1 17 (12.3%) 23 (15.9%) 0.732 (0.369 -1.453) 0.37 0.684 (0.339 -1.379) 0.28 100 (72.5%) 99 (68.3%) 0.904 (0.470 -1.738) 0.76 0.893 (0.457 -1.745) 0.74 Dominant A/AA/G+G/G 21 (15.2%) 23 (15.9%) 1 1 117 (84.8%) 122 (84.1%) 0.952 (0.500 -1.812) 0.88 0.933 (0.482 -1.804) 0. 83 Recessive A/A+A/GG/G 38 (27.5%) 46 (31.7%) 1 100 (72.5%) 99 (68.3%) 0.818 (0.490 -1.364) 0.44 0.780 (0.463 -1.315) 0. 35 Overdominant A/A+G/GA/G 121 (87.7%) 122 (84.1%) 1 1 17 (12.3%) 23 (15.9%) 0.745 (0.379 -1.464) 0.39 0.702 (0.352 -1.399) 0.31 Note. OR, Odds ratio; 95% CI, 95% confidence interval; P < 0.05, considered statistically significant; logistic regression model, OR adjusted for age, BMI, smoking and length of exposure to benzene. No significant differences were observed in the genotype and allele frequencies of the miR-506-rs5905008 A>G SNP between subjects with abnormal hematological indices and healthy controls. In addition, we did not find significant statistical values for the regression analysis data under different genetic models or for the results of its relationship with each subgroup of blood parameters. Nevertheless, when the combined genotypes of the miR-100-rs543412 C>T and miR-506-rs5905008 A>G polymorphisms were compared between the two groups using combined genotype analysis, it was revealed that the combined heterozygote genotype of the two SNPs was significantly different (P = 0.028) between the two groups and was associated with a reduced risk of blood index abnormalities (OR: 0.089, 95% CI: 0.009–0.857, P = 0.03) (Table 5).
Table 5. The results of the combined genotype analysis of the two polymorphisms under study.
Combined genotypes N (%) Cases Controls P-value OR (95%CI) P-value CCAA 17 (6.0) 9 (6.5) 8 (5.5)
0/028ref − CCAG 10 (3.5) 6 (4.3) 4 (2.8) 1.33 (0.274-6.49) 0.99 CCGG 52 (18.4) 33 (23.9) 19 (13.1) 1.544 (0.510- 4.670) 0.44 CTAA 16 (5.7) 11 (8.0) 5 (3.4) 1.956 (0.471- 8.11) 0.35 CTAG 25 (8.8) 9 (6.5) 16 (11.0) 0.089 (0.009- 0.857) 0.03 CTGG 109 (38.5) 51 (37.0) 58 (40.0) 0.500 (0.143- 1.753) 0.27 TTAA 11 (3.9) 1 (0.7) 10 (6.9) 0.782 (0.281- 2.176) 0.63 TTAG 5 (1.8) 2 (1.4) 3 (2.1) 0.593 (0.078- 4.498) 0.61 TTGG 38 (13.4) 16 (11.6) 22 (15.2) 0.646 (0.205- 2.041) 0.45 Note. Logistic regression model, OR, odds ratio, 95% CI 95% confidence interval P < 0.05 considered statistically significant.
doi: 10.3967/bes2024.110
A miR-100 Polymorphism Signature is Protectively Associated with Hematological Abnormalities in Individuals Exposed to Benzene, Toluene, Ethylbenzene, and Xylene
-
Abstract:
Objective The DNA double-strand break (DSB) repair system plays a key role in eliminating DNA damage in hematopoietic cells caused by benzene, toluene, ethylbenzene, and xylene (BTEX) compounds, and the polymorphisms of genes controlling this DNA repair system are linked to the different genetic susceptibilities of individuals to respond to the effects of BTEXs. In addition, there is some evidence that BTEX can induce hematological abnormalities through changing the function of micro RNAs regulating the DSB system, but there is no report on the association of their polymorphisms with BTEX-related disorders. This study aimed to explore the relationship between important polymorphisms in miRNA-100 and miRNA-506, two key regulators of the hematopoietic DSB system, including miRNA-100 RS543412 C>T and miRNA-506 RS5905008 A>G, and hematological abnormalities in BTEX-exposed workers. Methods This study included 138 patients who were exposed to BTEX compounds and had one or more abnormalities in their hematological indices and 145 healthy individuals with the same exposure conditions but without any abnormal hematological defects. Polymorphism genotyping was performed using polymerase chain reaction restriction fragment length polymorphism. Results The results revealed that frequency of the T allele as well as the TT and CT genotypes of the miR-100-RS543412C>T single nucleotide polymorphism was significantly lower in patients than in healthy participants (allelic P = 0.005; genotype P = 0.013). Moreover, individuals with TT and CT genotypes in the codominant (CT vs. CC, OR: 0.546, P = 0.034; TT vs. CC, OR: 0.335, P = 0.004), dominant (CT+TT vs. CC, OR: 0.486, P = 0.008), and recessive (TT vs. CC+CT, OR: 0.507, P = 0.038) models had a lower risk of hematological abnormalities. We also observed that the combined heterozygous genotype of the two polymorphisms was significantly different (P = 0.028) between the two groups and was associated with a reduced risk of abnormalities in blood indices (OR: 0.089, P = 0.03). Conclusion These findings suggest that this miR-single nucleotide polymorphism may be a protective non-coding signature underlying the risk of BTEX exposure-related hematological abnormalities. However, this hypothesis requires further investigation. -
Key words:
- BTEX compounds /
- Genetic susceptibility /
- miR-100 SNP /
- rs543412C>T /
- Hematological abnormalities
&These authors contributed equally to this work.
注释:1) CONFLICTS OF INTEREST: -
Table 1. Characteristics of the BTEX and other variables in two groups of the study participants.
Variables Case (n=138) Control (n=145) P-values Age (years), mean ± SD 34.26 ± 6.50 34.47 ± 5.36 0.67 BMI (kg/m2), mean ± SD 25.85 ± 3.38 25.40 ± 4.78 0.84 Length of exposure in a day (hours), mean ± SD 11.80 ± 1.44 11.66 ± 1.45 0.45 Job history (year), mean ± SD 7.38 ± 3.75 7.84 ± 4.16 0.47 Smoking, n (%) 0.87 Yes 7 (5.3) 8 (5.8) NO 125 (94.7) 131 (94.2) BTEX-exposure measurements (ppm), mean ± SD Benzene 0.57 ± 1.51 0.71 ± 1.59 0.99 Toluene 1.02 ± 2.03 1.11 ± 2.17 0.88 Ethylbenzene 1.44 ± 2.52 1.24 ± 2.70 0.27 O-xylene 0.437 ± 1.33 0.414 ± 1.64 0.27 P-xylene 2.47 ± 3.06 2.28 ± 3.04 0.49 M xylene 1.70 ± 2.46 1.79 ± 2.64 0.89 Note. Student’s t-test was used to determine the statistical significance of the differences in age, BMI, daily exposure time, and job history, and the chi-square test was used for smoking history; P-value < 0.05 was considered statistically significant (in bold). BMI, body mass index; BTEX, benzene toluene ethylbenzene xylene. *Some data were missing due to lack of information. Table 2. Data on hematological indices in case and control subjects
Hematological indices Case (n = 138) Control (n=145) P-values WBC (×109/L) 6.97 ± 1.86 7.03 ± 1.51 0.80 RBC (×109/L) 5.04 ± 0.559 5. 0 ± 0.360 0.47 HB (g/dL) 14.77 ± 1.50 15.32 ± 0.840 0.001 HCT (%) 43.52 ± 3.65 44.22 ± 2.24 0.051 MCV (fL) 86.20 ± 9.30 88.67 ± 3.82 0.004 MCH (Pg) 29.64 ± 3.67 30.75 ± 1.71 0.001 MCHC (g/dL) 34.10 ± 1.89 34.66 ± 0.911 0.002 RDW (%) 12.41 ± 0.968 12.13 ± 0.562 0.018 Platelets (×109/L) 212 ± 46.64 210.8 ± 38.52 0.908 Lymphocytes (%) 44.72 ± 10.87 37.96 ± 5.72 0.001 Monocytes (%) 2.92 ± 1.59 3.08 ± 1.79 0.936 Granulocytes (%) 51.07 ± 10.68 57.97 ± 6.38 0.001 Note. Values are given as the mean ± standard deviation; Student t-test for comparing the difference of hematological indices in the case and control groups was used; P-value <0.05 was considered statistically significant (in bold). WBC, white blood cell count; RBC, red blood cell count; HB, hemoglobin; HCT, hematocrit; MCV, mean corpuscular volume; MCH, mean corpuscular hemoglobin; MCHC, mean corpuscular hemoglobin concentration; RDW, red cell distribution width. Table 3. Genotype and allele distribution of the selected polymorphisms in study population groups.
SNPs Cases (n=138),
N (%)Controls (n=145)
N (%)P-Value P-HWE miR-100-rs543412c>T Genotypes
CC
CT
TT
48(34.8)
71(51.4)
19(13.8)
31(21.4)
79(54.5)
35(24.1)
0.013
0.24Alleles
C
T
167(60.5)
109(39.5)
141(48.6)
149(51.4)0.005 miR-506-rs5905008a>G
Genotypes
AA
AG
GG
21(15.2)
17(12.3)
100(72.5)
23(15.9)
23(15.9)
99(68.3)
0.66
0.001Alleles
A
G
217(78.6)
59(21.4)
221(76.2)
69(23.8)
0.49Note. Person χ2 test used for difference in distributions between the case and control groups. A goodness-of-fit chi-squared test was used to evaluate the Hardy-Weinberg equilibrium in the study. population; P-value < 0.05 was considered statistically significant (in bold). SNPs, single nucleotide polymorphisms; P-HWE, P value for Hardy-Weinberg equilibrium. Table 4. Association results between miR-SNP genotypes and patterns of hematological parameters in multiple inheritance models.
SNPs Genotype models Genotypes Population groups Unadjusted OR (95%CI) P Adjusted OR (95%CI) P Case (n = 138) Control (n = 145) miR-100-RS543412C>T Codominant C/CC/TT/T 48 (34.8%) 31 (21.4%) 1 1 71 (51.4%) 79 (54.5%) 0.580 (0.334-1.010) 0.054 0.546 (0.312-0.957) 0.034 19 (13.8%) 35 (24.1%) 0.351 (0.171-0.719) 0.004 0.335 (0.158-0.711) 0.004 Dominant C/CC/T+T/T 48 (34.8%) 31 (21.4%) 1 1 90 (65.2%) 114 (78.6%) 0.510 (0.300-0.866) 0.013 0.486 (0.285-0.830) 0.008 Recessive C/C+C/TT/T 119 (86.2%) 110 (75. 9%) 1 1 19 (13.8%) 35 (24.1%) 0.502 (0.271-0.929) 0.028 0.507 (0.267-0.962) 0.038 Overdominant C/C+T/TC/T 67 (48.6%) 66 (45.5%) 1 1 71 (51.4%) 79 (54.4%) 0.885 (0.555-1.413) 0.609 0.822 (0.510 -1.325) 0.421 miR-506-RS5905008A>G Codominant A/AA/GG/G 21 (15.2%) 23 (15.9%) 1 1 17 (12.3%) 23 (15.9%) 0.732 (0.369 -1.453) 0.37 0.684 (0.339 -1.379) 0.28 100 (72.5%) 99 (68.3%) 0.904 (0.470 -1.738) 0.76 0.893 (0.457 -1.745) 0.74 Dominant A/AA/G+G/G 21 (15.2%) 23 (15.9%) 1 1 117 (84.8%) 122 (84.1%) 0.952 (0.500 -1.812) 0.88 0.933 (0.482 -1.804) 0. 83 Recessive A/A+A/GG/G 38 (27.5%) 46 (31.7%) 1 100 (72.5%) 99 (68.3%) 0.818 (0.490 -1.364) 0.44 0.780 (0.463 -1.315) 0. 35 Overdominant A/A+G/GA/G 121 (87.7%) 122 (84.1%) 1 1 17 (12.3%) 23 (15.9%) 0.745 (0.379 -1.464) 0.39 0.702 (0.352 -1.399) 0.31 Note. OR, Odds ratio; 95% CI, 95% confidence interval; P < 0.05, considered statistically significant; logistic regression model, OR adjusted for age, BMI, smoking and length of exposure to benzene. Table 5. The results of the combined genotype analysis of the two polymorphisms under study.
Combined genotypes N (%) Cases Controls P-value OR (95%CI) P-value CCAA 17 (6.0) 9 (6.5) 8 (5.5)
0/028ref − CCAG 10 (3.5) 6 (4.3) 4 (2.8) 1.33 (0.274-6.49) 0.99 CCGG 52 (18.4) 33 (23.9) 19 (13.1) 1.544 (0.510- 4.670) 0.44 CTAA 16 (5.7) 11 (8.0) 5 (3.4) 1.956 (0.471- 8.11) 0.35 CTAG 25 (8.8) 9 (6.5) 16 (11.0) 0.089 (0.009- 0.857) 0.03 CTGG 109 (38.5) 51 (37.0) 58 (40.0) 0.500 (0.143- 1.753) 0.27 TTAA 11 (3.9) 1 (0.7) 10 (6.9) 0.782 (0.281- 2.176) 0.63 TTAG 5 (1.8) 2 (1.4) 3 (2.1) 0.593 (0.078- 4.498) 0.61 TTGG 38 (13.4) 16 (11.6) 22 (15.2) 0.646 (0.205- 2.041) 0.45 Note. Logistic regression model, OR, odds ratio, 95% CI 95% confidence interval P < 0.05 considered statistically significant. -
[1] Baberi Z, Azhdarpoor A, Hoseini M, et al. Monitoring benzene, toluene, ethylbenzene, and xylene (BTEX) levels in mixed-use residential-commercial buildings in Shiraz, Iran: assessing the carcinogenicity and non-carcinogenicity risk of their inhabitants. Int J Environ Res Public Health, 2022; 19, 723. doi: 10.3390/ijerph19020723 [2] Bolden AL, Kwiatkowski CF, Colborn T. New look at BTEX: are ambient levels a problem? Environ Sci Technol, 2015; 49, 5261-76. [3] Rafiee A, Delgado-Saborit JM, Sly PD, et al. Lifestyle and occupational factors affecting exposure to BTEX in municipal solid waste composting facility workers. Sci Total Environ, 2019; 656, 540−6. doi: 10.1016/j.scitotenv.2018.11.398 [4] Sagbo F, Lawin HB, Atindehou M, et al. Effects of BTEX exposure on hematological and c-reactive-protein in professional and non professional motorcycle drivers in cotonou/benin. J Environ Pollut Hum Health, 2020; 8, 1−5. [5] Chen Q, Sun H, Zhang JY, et al. The hematologic effects of BTEX exposure among elderly residents in Nanjing: a cross-sectional study. Environ Sci Pollut Res, 2019; 26, 10552−61. doi: 10.1007/s11356-019-04492-9 [6] Jafari Roshan S, Mansoori Y, Hosseini SR, et al. Genetic variations in ATM and H2AX loci contribute to risk of hematological abnormalities in individuals exposed to BTEX chemicals. J Clin Lab Anal, 2022; 36, e24321. doi: 10.1002/jcla.24321 [7] Lim SK, Shin HS, Yoon KS, et al. Risk assessment of volatile organic compounds benzene, toluene, ethylbenzene, and xylene (BTEX) in consumer products. J Toxicol Environ Health Part A, 2014; 77, 1502−21. doi: 10.1080/15287394.2014.955905 [8] Isinkaralar K. Theoretical removal study of gas BTEX onto activated carbon produced from Digitalis purpurea L. biomass. Biomass Convers Biorefin, 2022; 12, 4171−81. doi: 10.1007/s13399-022-02558-2 [9] Allahabady A, Yousefi Z, Ali Mohammadpour Tahamtan R, et al. Measurement of BTEX (benzene, toluene, ethylbenzene and xylene) concentration at gas stations. Environ Health Eng Manag, 2022; 9, 23−31. doi: 10.34172/EHEM.2022.04 [10] Adebambo TH, Fox DT, Otitoloju AA. Toxicological study and genetic basis of BTEX susceptibility in Drosophila melanogaster. Front Genet, 2020; 11, 594179. doi: 10.3389/fgene.2020.594179 [11] Xiong F, Li Q, Zhou B, et al. Oxidative stress and genotoxicity of long-term occupational exposure to low levels of BTEX in gas station workers. Int J Environ Res Public Health, 2016; 13, 1212. doi: 10.3390/ijerph13121212 [12] Cordiano R, Papa V, Cicero N, et al. Effects of benzene: hematological and hypersensitivity manifestations in resident living in oil refinery areas. Toxics, 2022; 10, 678. doi: 10.3390/toxics10110678 [13] Shen M, Lan Q, Zhang L P, et al. Polymorphisms in genes involved in DNA double-strand break repair pathway and susceptibility to benzene-induced hematotoxicity. Carcinogenesis, 2006; 27, 2083−9. doi: 10.1093/carcin/bgl061 [14] Li N, Chen HZ, Wang J. DNA damage and repair in the hematopoietic system. Acta Biochim Biophys Sin, 2022; 54, 847−57. doi: 10.3724/abbs.2022053 [15] Valikhani M, Shokuhian M, Gholami MS, et al. DNA double-strand break repair and hematologic malignancy, a review. SID, https://sid.ir/paper/946596/en. [2017]. [16] Zhao BL, Rothenberg E, Ramsden DA, et al. The molecular basis and disease relevance of non-homologous DNA end joining. Nat Rev Mol Cell Biol, 2020; 21, 765−81. doi: 10.1038/s41580-020-00297-8 [17] Peraza-Vega RI, Valverde M, Rojas E. Interactions between miRNAs and double-strand breaks DNA repair genes, pursuing a fine-tuning of repair. Int J Mol Sci, 2022; 23, 3231. doi: 10.3390/ijms23063231 [18] Szatkowska M, Krupa R. Regulation of DNA damage response and homologous recombination repair by microRNA in human cells exposed to ionizing radiation. Cancers, 2020; 12, 1838. doi: 10.3390/cancers12071838 [19] Wang YD, Cheng J, Li DC, et al. Modulation of DNA repair capacity by ataxia telangiectasia mutated gene polymorphisms among polycyclic aromatic hydrocarbons-exposed workers. Toxicol Sci, 2011; 124, 99−108. doi: 10.1093/toxsci/kfr216 [20] Lan Q, Zhang LP, Shen M, et al. Polymorphisms in cytokine and cellular adhesion molecule genes and susceptibility to hematotoxicity among workers exposed to benzene. Cancer Res, 2005; 65, 9574-81. [21] Takagi M. DNA damage response and hematological malignancy. Int J Hematol, 2017; 106, 345−56. doi: 10.1007/s12185-017-2226-0 [22] Zhang GH, Ren JC, Luo MK, et al. Association of BER and NER pathway polymorphism haplotypes and micronucleus frequencies with global DNA methylation in benzene-exposed workers of China: effects of DNA repair genes polymorphisms on genetic damage. Mutat Res Genet Toxicol Environ Mutagen, 2019; 839, 13−20. doi: 10.1016/j.mrgentox.2019.01.006 [23] SiamiGorji S, Jorjani I, Tahamtan A, et al. Effects of microRNAs polymorphism in cancer progression. Med J Islam Repub Iran, 2020; 34, 3. [24] León-Mejía G, Quintana-Sosa M, de Moya Hernandez Y, et al. DNA repair and metabolic gene polymorphisms affect genetic damage due to diesel engine exhaust exposure. Environ Sci Pollut Res Int, 2020; 27, 20516−26. doi: 10.1007/s11356-020-08533-6 [25] Sisto R, Cavallo D, Ursini CL, et al. Direct and oxidative dna damage in a group of painters exposed to VOCs: dose – response relationship. Front Public Health, 2020; 8, 445. doi: 10.3389/fpubh.2020.00445 [26] Wu HZ, Li SQ, Hu XY, et al. Associations of mRNA expression of DNA repair genes and genetic polymorphisms with cancer risk: a bioinformatics analysis and meta-analysis. J Cancer, 2019; 10, 3593−607. doi: 10.7150/jca.30975 [27] Vodicka P, Stetina R, Polakova V, et al. Association of DNA repair polymorphisms with DNA repair functional outcomes in healthy human subjects. Carcinogenesis, 2007; 28, 657−64. [28] Mousavi N, Tafvizi F, Mansoori Y. Genetic polymorphisms of base excision repair gene XRCC1 and susceptibility to benzene among employees of chemical industries. Gene Rep, 2021; 23, 101081. doi: 10.1016/j.genrep.2021.101081 [29] Li YX, Tong Y, Liu JQ, et al. The role of MicroRNA in DNA damage response. Front Genet, 2022; 13, 850038. doi: 10.3389/fgene.2022.850038 [30] Wei HY, Zhang J, Tan KH, et al. Benzene-induced aberrant miRNA expression profile in hematopoietic progenitor cells in C57BL/6 mice. Int J Mol Sci, 2015; 16, 27058−71. doi: 10.3390/ijms161126001 [31] Ma XJ, Zhang X, Luo J, et al. MiR-486-5p-directed MAGI1/Rap1/RASSF5 signaling pathway contributes to hydroquinone-induced inhibition of erythroid differentiation in K562 cells. Toxicol in Vitro, 2020; 66, 104830. doi: 10.1016/j.tiv.2020.104830 [32] Yu CH, Yang SQ, Li L, et al. Identification of potential pathways and microRNA-mRNA networks associated with benzene metabolite hydroquinone-induced hematotoxicity in human leukemia K562 cells. BMC Pharmacol Toxicol, 2022; 23, 20. doi: 10.1186/s40360-022-00556-8 [33] Hu DL, Peng XW, Liu YG, et al. Overexpression of miR-221 in peripheral blood lymphocytes in petrol station attendants: a population based cross-sectional study in southern China. Chemosphere, 2016; 149, 8−13. doi: 10.1016/j.chemosphere.2016.01.083 [34] Yokoi T, Nakajima M. microRNAs as mediators of drug toxicity. Annu Rev Pharmacol Toxicol, 2013; 53, 377−400. doi: 10.1146/annurev-pharmtox-011112-140250 [35] Mishra PJ, Bertino JR. MicroRNA polymorphisms: the future of pharmacogenomics, molecular epidemiology and individualized medicine. Pharmacogenomics, 2009; 10, 399−416. doi: 10.2217/14622416.10.3.399 [36] Konoshenko MY, Bryzgunova OE, Laktionov PP. miRNAs and androgen deprivation therapy for prostate cancer. Biochim Biophys Acta Rev Cancer, 2021; 1876, 188625. doi: 10.1016/j.bbcan.2021.188625 [37] Sun Y, Wang HX, Luo CB. MiR-100 regulates cell viability and apoptosis by targeting ATM in pediatric acute myeloid leukemia. Biochem Biophys Res Commun, 2020; 522, 855−61. doi: 10.1016/j.bbrc.2019.11.156 [38] Said F, Tantawy M, Sayed A, et al. Clinical significance of MicroRNA-29a and MicroRNA-100 gene expression in pediatric acute myeloid leukemia. J Pediatr Hematol/Oncol, 2022; 44, e391−5. doi: 10.1097/MPH.0000000000002168 [39] Zheng YS, Zhang H, Zhang XJ, et al. MiR-100 regulates cell differentiation and survival by targeting RBSP3, a phosphatase-like tumor suppressor in acute myeloid leukemia. Oncogene, 2012; 31, 80−92. doi: 10.1038/onc.2011.208 [40] Bao XH, Liu XJ, Li F, et al. Limited MOMP, ATM, and their roles in carcinogenesis and cancer treatment. Cell Biosci, 2020; 10, 81. doi: 10.1186/s13578-020-00442-y [41] Xu JF, Bradley N, He Y. Structure and function of the apical PIKKs in double-strand break repair. Curr Opin Struct Biol, 2023; 82, 102651. doi: 10.1016/j.sbi.2023.102651 [42] Tan JP, Sun XY, Zhao HL, et al. Double‐strand DNA break repair: molecular mechanisms and therapeutic targets. MedComm, 2023; 4, e388. doi: 10.1002/mco2.388 [43] Roden C, Gaillard J, Kanoria S, et al. Novel determinants of mammalian primary microRNA processing revealed by systematic evaluation of hairpin-containing transcripts and human genetic variation. Genome Res, 2017; 27, 374−84. doi: 10.1101/gr.208900.116 [44] Xue Y, Yang XY, Hu SY, et al. A genetic variant in miR‐100 is a protective factor of childhood acute lymphoblastic leukemia. Cancer Med, 2019; 8, 2553−60. doi: 10.1002/cam4.2082 [45] Liu GY, Yang D, Rupaimoole R, et al. Augmentation of response to chemotherapy by microRNA-506 through regulation of RAD51 in serous ovarian cancers. J Natl Cancer Inst, 2015; 107, djv108. [46] Liu GY, Xue FX, Zhang W. miR-506: a regulator of chemo-sensitivity through suppression of the RAD51-homologous recombination axis. Chin J Cancer, 2015; 34, 1−3. [47] Chen NF, Meng Z, Song JJ, et al. miR-506 in patients with chronic myeloid leukemia and its effect on apoptosis of K562 cells. Am J Transl Res, 2021; 13, 9413−20. [48] Zhu XW, Wang J, Zhu MX, et al. MicroRNA-506 inhibits the proliferation and invasion of mantle cell lymphoma cells by targeting B7-H3. Biochem Biophys Res Commun, 2019; 508, 1067−73. doi: 10.1016/j.bbrc.2018.12.055 [49] Health. Division of Physical Sciences. NIOSH manual of analytical methods. 4th ed. US Department of Health and Human Services, Public Health Service, Centers for Disease Control and Prevention, National Institute for Occupational Safety and Health, Division of Physical Sciences and Engineering, 1994. [50] Minatel BC, Sage AP, Anderson C, et al. Environmental arsenic exposure: from genetic susceptibility to pathogenesis. Environ Int, 2018; 112, 183−97. doi: 10.1016/j.envint.2017.12.017 [51] Lash LH. Environmental and genetic factors influencing kidney toxicity. Semin Nephrol, 2019; 39, 132−40. doi: 10.1016/j.semnephrol.2018.12.003 [52] Abubakar MB, Sanusi KO. Influence of GSTM1 and GSTT1 genetic polymorphisms on petrol-induced toxicities: a systematic review. Meta Gene, 2020; 26, 100796. doi: 10.1016/j.mgene.2020.100796 [53] Chiarella P, Capone P, Sisto R. The role of genetic polymorphisms in the occupational exposure. In: Çalışkan M, Erol O, Öz G C. The Recent Topics in Genetic Polymorphisms. IntechOpen. 2019. [54] Jalilian S, Sabzalipour S, Mohammadi Rouzbahani M, et al. Assessing the effect of BTEX on blood and spirometry parameters staff in a petroleum refinery. Front Public Health, 2022; 10, 1037413. doi: 10.3389/fpubh.2022.1037413 [55] Mozzoni P, Poli D, Pinelli S, et al. Benzene exposure and MicroRNAs expression: in vitro, in vivo and human findings. Int J Environ Res Public Health, 2023; 20, 1920. doi: 10.3390/ijerph20031920 [56] Dai K, Wang C, Yao W, et al. Expression level and function analysis of serum miRNAs in workers with occupational exposure to benzene series. Chemosphere, 2023; 313, 137460. doi: 10.1016/j.chemosphere.2022.137460 [57] Sisto R, Capone P, Cerini L, et al. Occupational exposure to volatile organic compounds affects microRNA profiling: Towards the identification of novel biomarkers. Toxicol Rep, 2020; 7, 700−10. doi: 10.1016/j.toxrep.2020.05.006 [58] Tessitore A, Cicciarelli G, Del Vecchio F, et al. MicroRNAs in the DNA damage/repair network and cancer. Int J Genomics, 2014; 2014, 820248. [59] Santos JSD, Zunta GL, Negrini AB, et al. The association of a single-nucleotide variant in the microRNA-146a with advanced colorectal cancer prognosis. Tumour Biol, 2020; 42, 1010428320923856. [60] Machowska M, Galka-Marciniak P, Kozlowski P. Consequences of genetic variants in miRNA genes. Comput Struct Biotechnol J, 2022; 20, 6443−57. doi: 10.1016/j.csbj.2022.11.036 [61] Teklu G, Negash M, Asefaw T, et al. Effect of gasoline exposure on hematological parameters of gas station workers in Mekelle city, Tigray Region, Northern Ethiopia. J Blood Med, 2021; 839-47. [62] Abdulkareem AQ, Mohammed BMA. Occupational risk assessments of hematological parameters alterations of kar oil refinery workers in Erbil Province, Kurdistan Region, Iraq. J Duhok Univ, 2023; 26, 106−17. doi: 10.26682/sjuod.2023.26.2.10 [63] Doherty BT, Kwok RK, Curry MD, et al. Associations between blood BTEXS concentrations and hematologic parameters among adult residents of the U. S. Gulf States. Environ Res, 2017; 156, 579−87. doi: 10.1016/j.envres.2017.03.048 [64] Zhang ZR, Liu X, Guo CF, et al. Hematological effects and benchmark doses of long-term co-exposure to benzene, toluene, and xylenes in a follow-up study on petrochemical workers. Toxics, 2022; 10, 502. doi: 10.3390/toxics10090502 [65] Kasemy ZA, Kamel GM, Abdel-Rasoul GM, et al. Environmental and health effects of benzene exposure among Egyptian taxi drivers. J Environ Public Health, 2019; 2019; 7078024. [66] Stubbins RJ, Korotev S, Godley LA. Germline CHEK2 and ATM variants in myeloid and other hematopoietic malignancies. Curr Hematol Malig Rep, 2022; 17, 94−104. doi: 10.1007/s11899-022-00663-7 [67] Parvez A, Mahjabeen I, Mehmood A, et al. Expression variations of DNA damage response genes ATM and ATR in blood cancer patients. Mol Genet Genomics, 2023; 298, 1173−83. doi: 10.1007/s00438-023-02043-z [68] Fortin J, Bassi C, Ramachandran P, et al. Concerted roles of PTEN and ATM in controlling hematopoietic stem cell fitness and dormancy. J Clin Invest, 2021; 131, e131698. doi: 10.1172/JCI131698 [69] Kyriakidis I, Kyriakidis K, Tsezou A. Micrornas and the diagnosis of childhood acute lymphoblastic leukemia: systematic review, meta-analysis and re-analysis with novel small RNA-Seq tools. Cancers, 2022; 14, 3976. doi: 10.3390/cancers14163976 [70] Kargar M, Farsani MA, Garavand J, et al. Engaging of the mTOR signaling pathway by miR100 and miR101 in de novo acute myeloid leukemia. Res Square, 2022, doi: 10.21203/rs.3.rs-2248376/v1. [71] Guo YH, Deng XD, Dai K, et al. Benchmark dose estimation based on oxidative damage in Chinese workers exposed to benzene series compounds. Environ Toxicol Pharmacol, 2023; 100, 104150. doi: 10.1016/j.etap.2023.104150 [72] Montero-Montoya RD, López-Vargas R, Méndez-Serrano A, et al. Increased micronucleus frequencies in reticulocytes of children exposed to industrial pollution: oxidative stress and the OGG1 S326C polymorphism. Mutat Res Genet Toxicol Environ Mutagen, 2020; 853, 503170. doi: 10.1016/j.mrgentox.2020.503170 [73] Chen CS, Hseu YC, Liang SH, et al. Assessment of genotoxicity of methyl-tert-butyl ether, benzene, toluene, ethylbenzene, and xylene to human lymphocytes using comet assay. J Hazard Mater, 2008; 153, 351−6. doi: 10.1016/j.jhazmat.2007.08.053 [74] Zhang LG, Zhang M, Wang H, et al. Comprehensive review of genetic association studies and meta-analysis on polymorphisms in microRNAs and urological neoplasms risk. Sci Rep, 2018; 8, 3776. doi: 10.1038/s41598-018-21749-4 [75] Gago-Fuentes R, Oksenych V. Non-homologous end joining factors XLF, PAXX and DNA-PKcs maintain the neural stem and progenitor cell population. Biomolecules, 2021; 11, 20. [76] Nourozi MA, Neghab M, Bazzaz JT, et al. Association between polymorphism of GSTP1, GSTT1, GSTM1 and CYP2E1 genes and susceptibility to benzene-induced hematotoxicity. Arch Toxicol, 2018; 92, 1983−90. doi: 10.1007/s00204-017-2104-9 [77] Xing CH, Chen Q, Li GL, et al. Microsomal epoxide hydrolase (EPHX1) polymorphisms are associated with aberrant promoter methylation of ERCC3 and hematotoxicity in benzene‐exposed workers. Environ Mol Mutagen, 2013; 54, 397−405. doi: 10.1002/em.21786