The characteristics of the subjects were listed in Table 1. A total of 510 TB patients including 318 males and 192 females and 508 unrelated healthy controls involving 316 males and 192 females were employed. All participants were of Han ethnicity, and lived in or around the Xi’an city. The mean ages of TB cases and healthy controls were 41.90 ± 14.83 years old and 41.14 ± 18.42 years old, respectively. No significant differences between all individuals in age or gender proportion were observed.
Table 1
Demographic distribution of TB patients and Healthy controls
|
Cases (n = 510) No. (%)
|
Controls (n = 508) No. (%)
|
p
|
Age, year (mean ± SD)
|
41.90 ± 14.83
|
41.14 ± 18.42
|
0.469
|
Gender
|
|
|
0.961
|
Male
|
318 (62.35)
|
316 (62.20)
|
|
Female
|
192 (37.65)
|
192 (37.80)
|
|
In our study, the detail information of these SNPs (the position, minor allele frequency (MAF) and HWE p value) were summarized in Table 2. We found that rs6531668 of TLR6 gene was association with TB risk (OR = 1.28, 95% CI = 1.04–1.57, p = 0.023). Analysis of the genotype frequencies of TLR6 in different genetic models, our results suggested that rs6531668 was correlated with increased the risk of TB in dominant model (OR = 1.30, 95% CI = 1.01–1.67, p = 0.046) and log–additive model (OR = 1.26, 95% CI = 1.03–1.55, p = 0.026). However, no significant association was identified between rs3775073, rs3796508, rs12644563, rs12640631, rs12642845 and TB risk.
In stratified analysis by age, the results confirmed that rs6531668 was increased TB risk in age > 41 years under the co–dominant model (OR = 1.75, 95% CI = 1.14–2.68, p = 0.010), the dominant model (OR = 1.77, 95% CI = 1.19–2.62, p = 0.005) and the log–additive model (OR = 1.51, 95% CI = 1.11–2.06, p = 0.005). Interestingly, we observed that the G/C – C/C genotype of rs12642845 was significantly interacted with decreased TB risk in age > 41 years old under the dominant model (OR = 0.67, 95% CI = 0.46–0.99, p = 0.042) (Table 4). When gender stratification was analyzed, the results confirmed that rs3775073 G/A genotype and rs12642845 G/C genotype were decreased TB risk in males (OR = 0.70, 95% CI = 0.50–0.98, p = 0.036; OR = 0.71, 95% CI = 0.51–0.99, p = 0.042, respectively). However, no association of rs3775073 and rs12642845 and TB risk were observed in females. The results were shown in Table 5.
We further performed LD and haplotype analyses of all SNPs in our study. Strong LD analysis displayed the existence of two blocks in TLR6 SNPs (Fig. 1). LD was found between rs3775073 and rs6531688, among rs12640631 and rs12642845. Besides, we observed that the CTA and TCG haplotypes were associated with reduced risk of TB (OR = 0.74, 95% CI = 0.55–0.99, p = 0.043; OR = 0.47, 95% CI = 0.28–0.79, p = 0.004, respectively) (Table 6).
Moreover, we used Regulome DB and Haploreg v4.1 to determine the potential functional effect of these SNPs. The Regulome DB suggested that rs6531688 and rs3775073 were the binding site for the transcription factor. The Regulome DB score was 1f and 2b. No data were shown the prediction functional of rs3796508 and rs12640631 (Supplementary Table 1).
Table 2
Basic characteristics and allele frequencies among TLR6 SNPs
SNP
|
Chr
|
Allele
|
MAF
|
HWE p–Value
|
OR (95% CI)
|
pa
|
Case
|
Control
|
rs3775073
|
4
|
A/G
|
0.31
|
0.33
|
0.163
|
0.89 (0.74–1.07)
|
0.205
|
rs3796508
|
4
|
A/G
|
0.04
|
0.05
|
0.639
|
0.72 (0.47–1.10)
|
0.130
|
rs6531668
|
4
|
A/G
|
0.25
|
0.20
|
0.276
|
1.28 (1.04–1.57)
|
0.023*
|
rs12644563
|
4
|
C/T
|
0.04
|
0.05
|
0.148
|
0.72 (0.47–1.10)
|
0.130
|
rs12640631
|
4
|
C/T
|
0.08
|
0.07
|
0.340
|
1.14 (0.83–1.58)
|
0.425
|
rs12642845
|
4
|
C/G
|
0.21
|
0.24
|
0.272
|
0.87 (0.71–1.07)
|
0.188
|
Table 3
The association between six SNPs within the TLR6gene and the risk of TB.
SNP
|
Model
|
Genotype
|
Cases
|
Controls
|
OR (95%CI)
|
pa
|
rs3775073
|
Co–dominant
|
G/G
|
245
|
219
|
1.00
|
|
|
|
G/A
|
216
|
240
|
0.80 (0.62–1.04)
|
0.090
|
|
|
A/A
|
48
|
49
|
0.88 (0.57–1.36)
|
0.556
|
|
Dominant
|
G/G
|
245
|
219
|
1.00
|
|
|
|
G/A – A/A
|
264
|
289
|
0.81 (0.63–1.04)
|
0.100
|
|
Recessive
|
G/G – G/A
|
461
|
459
|
1.00
|
|
|
|
A/A
|
48
|
49
|
0.98 (0.64–1.49)
|
0.919
|
|
Log–additive
|
–
|
–
|
–
|
0.88 (0.73–1.07)
|
0.191
|
rs3796508
|
Co–dominant
|
G/G
|
471
|
457
|
1.00
|
|
|
|
G/A
|
39
|
49
|
0.78 (0.50–1.21)
|
0.262
|
|
|
A/A
|
0
|
2
|
–
|
–
|
|
Dominant
|
G/G
|
471
|
457
|
1.00
|
|
|
|
G/A – A/A
|
39
|
51
|
0.75 (0.48–1.16)
|
0.190
|
|
Recessive
|
G/G – G/A
|
510
|
506
|
1.00
|
|
|
|
A/A
|
0
|
2
|
–
|
–
|
|
Log–additive
|
–
|
–
|
–
|
0.73 (0.48–1.11)
|
0.139
|
rs6531668
|
Co–dominant
|
G/G
|
294
|
326
|
1.00
|
|
|
|
G/A
|
175
|
157
|
1.24 (0.95–1.62)
|
0.118
|
|
|
A/A
|
37
|
25
|
1.65 (0.97–2.81)
|
0.065
|
|
Dominant
|
G/G
|
294
|
326
|
1.00
|
|
|
|
G/A – A/A
|
212
|
182
|
1.30 (1.01–1.67)
|
0.046*
|
|
Recessive
|
G/G – G/A
|
469
|
483
|
1.00
|
|
|
|
A/A
|
37
|
25
|
1.53 (0.91–2.59)
|
0.110
|
|
Log–additive
|
–
|
–
|
–
|
1.26 (1.03–1.55)
|
0.026*
|
rs12644563
|
Co–dominant
|
T/T
|
471
|
458
|
1.00
|
|
|
|
T/C
|
39
|
47
|
0.81 (0.52–1.27)
|
0.363
|
|
|
C/C
|
0
|
3
|
–
|
–
|
|
Dominant
|
T/T
|
471
|
458
|
1.00
|
|
|
|
T/C–C/C
|
39
|
50
|
0.76 (0.49–1.18)
|
0.228
|
|
Recessive
|
T/T– T/C
|
510
|
505
|
1.00
|
|
|
|
C/C
|
0
|
3
|
–
|
–
|
|
Log–additive
|
–
|
–
|
–
|
0.73 (0.48–1.11)
|
0.144
|
rs12640631
|
Co–dominant
|
T/T
|
430
|
437
|
1.00
|
|
|
|
T/C
|
75
|
67
|
1.15 (0.80–1.64)
|
0.453
|
|
|
C/C
|
5
|
4
|
1.26 (0.34–4.73)
|
0.733
|
|
Dominant
|
T/T
|
430
|
437
|
1.00
|
|
|
|
T/C–C/C
|
80
|
71
|
1.15 (0.81–1.63)
|
0.423
|
|
Recessive
|
T/T– T/C
|
505
|
504
|
1.00
|
|
|
|
C/C
|
5
|
4
|
1.24 (0.33–4.63)
|
0.753
|
|
Log–additive
|
–
|
–
|
–
|
1.14 (0.83–1.57)
|
0.417
|
rs12642845
|
Co–dominant
|
G/G
|
318
|
289
|
1.00
|
|
|
|
G/C
|
165
|
195
|
1.03 (0.58–1.82)
|
0.925
|
|
|
C/C
|
27
|
24
|
0.77 (0.59–1.00)
|
0.050
|
|
Dominant
|
G/G
|
318
|
289
|
1.00
|
|
|
|
G/C–CC
|
192
|
219
|
0.80 (0.62–1.03)
|
0.078
|
|
Recessive
|
G/G–G/C
|
483
|
484
|
1.00
|
|
|
|
C/C
|
27
|
24
|
1.13 (0.64–1.99)
|
0.665
|
|
Log–additive
|
–
|
–
|
–
|
0.87 (0.71–1.07)
|
0.192
|
Table 4
Age stratification analysis of association between TLR6 and the risk of TB.
SNP
|
Model
|
Genotype
|
≤ 41
|
> 41
|
OR (95% CI)
|
Pa
|
OR (95% CI)
|
Pa
|
rs3775073
|
Co–dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A
|
0.76 (0.52–1.10)
|
0.141
|
0.80 (0.54–1.19)
|
0.275
|
|
|
A/A
|
0.82 (0.43–1.57)
|
0.547
|
0.86 (0.45–1.63)
|
0.641
|
|
Dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A – A/A
|
0.77 (0.54–1.09)
|
0.143
|
0.81 (0.56–1.18)
|
0.279
|
|
Recessive
|
G/G – G/A
|
1.00
|
|
1.00
|
|
|
|
A/A
|
0.94 (0.51–1.76)
|
0.855
|
0.95 (0.52–1.76)
|
0.880
|
|
Log–additive
|
–
|
0.84 (0.64–1.11)
|
0.221
|
0.88 (0.66–1.17)
|
0.377
|
rs3796508
|
Co–dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A
|
0.95 (0.49–1.85)
|
0.884
|
–
|
–
|
|
|
A/A
|
–
|
–
|
–
|
–
|
|
Dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A – A/A
|
0.92 (0.48–1.76)
|
0.792
|
0.60 (0.32–1.14)
|
0.118
|
|
Recessive
|
G/G – G/A
|
1.00
|
|
1.00
|
|
|
|
A/A
|
–
|
–
|
–
|
–
|
|
Log–additive
|
–
|
0.89 (0.47–1.66)
|
0.706
|
0.60 (0.32–1.14)
|
0.118
|
rs6531668
|
Co–dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A
|
1.10 (0.76–1.61)
|
0.613
|
1.75 (1.14–2.68)
|
0.010*
|
|
|
A/A
|
1.73 (0.76–3.95)
|
0.194
|
1.83 (0.87–3.85)
|
0.111
|
|
Dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A – A/A
|
1.17 (0.81–1.67)
|
0.405
|
1.77 (1.19–2.62)
|
0.005*
|
|
Recessive
|
G/G – G/A
|
1.00
|
|
1.00
|
|
|
|
A/A
|
1.67 (0.74–3.76)
|
0.219
|
1.53 (0.74–3.17)
|
0.256
|
|
Log–additive
|
–
|
1.19 (0.88–1.61)
|
0.250
|
1.51 (1.11–2.06)
|
0.009*
|
rs12644563
|
Co–dominant
|
T/T
|
1.00
|
|
1.00
|
|
|
|
T/C
|
1.00 (0.51–1.97)
|
0.998
|
–
|
–
|
|
|
C/C
|
–
|
–
|
–
|
–
|
|
Dominant
|
T/T
|
1.00
|
|
1.00
|
|
|
|
T/C–C/C
|
0.89 (0.46–1.73)
|
0.736
|
0.67 (0.36–1.26)
|
0.212
|
|
Recessive
|
T/T– T/C
|
1.00
|
|
1.00
|
|
|
|
C/C
|
–
|
–
|
–
|
–
|
|
Log–additive
|
–
|
0.82 (0.44–1.51)
|
0.517
|
0.67 (0.36–1.26)
|
0.212
|
rs12640631
|
Co–dominant
|
T/T
|
1.00
|
|
1.00
|
|
|
|
T/C
|
1.03 (0.61–1.74)
|
0.913
|
1.23 (0.73–2.07)
|
0.437
|
|
|
C/C
|
0.73 (0.13–4.08)
|
0.719
|
–
|
–
|
|
Dominant
|
T/T
|
1.00
|
|
1.00
|
|
|
|
T/C–C/C
|
1.00 (0.60–1.67)
|
0.991
|
1.29 (0.77–2.16)
|
0.335
|
|
Recessive
|
T/T– T/C
|
1.00
|
|
1.00
|
|
|
|
C/C
|
0.73 (0.13–4.06)
|
0.716
|
–
|
–
|
|
Log–additive
|
–
|
0.98 (0.62–1.55)
|
0.931
|
1.34 (0.81–2.2)
|
0.255
|
rs12642845
|
Co–dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/C
|
0.82 (0.56–1.19)
|
0.293
|
0.62 (0.42–0.93)
|
0.021*
|
|
|
C/C
|
0.92 (0.40–2.12)
|
0.847
|
1.09 (0.46–2.56)
|
0.845
|
|
Dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/C–CC
|
0.83 (0.58–1.19)
|
0.310
|
0.67 (0.46–0.99)
|
0.042*
|
|
Recessive
|
G/G–G/C
|
1.00
|
|
1.00
|
|
|
|
C/C
|
1.00 (0.44–2.27)
|
0.992
|
1.30 (0.56–3.01)
|
0.541
|
|
Log–additive
|
–
|
0.88 (0.65–1.19)
|
0.392
|
0.79 (0.58–1.09)
|
0.152
|
Table 5
Gender stratification analysis of association between TLR6 and the risk of TB.
SNP
|
Model
|
Genotype
|
Males
|
Females
|
OR (95% CI)
|
Pa
|
OR (95% CI)
|
Pa
|
rs3775073
|
Co–dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A
|
0.70 (0.50–0.98)
|
0.036*
|
0.99 (0.66–1.51)
|
0.980
|
|
|
A/A
|
0.91 (0.54–1.53)
|
0.712
|
0.80 (0.35–1.81)
|
0.589
|
|
Dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A – A/A
|
0.74 (0.54–1.01)
|
0.059
|
0.97 (0.65–1.45)
|
0.872
|
|
Recessive
|
G/G – G/A
|
1.00
|
|
1.00
|
|
|
|
A/A
|
1.08 (0.65–1.77)
|
0.776
|
0.80 (0.36–1.76)
|
0.579
|
|
Log–additive
|
–
|
0.86 (0.68–1.09)
|
0.202
|
0.94 (0.68–1.31)
|
0.716
|
rs3796508
|
Co–dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A
|
0.86 (0.5–1.49)
|
0.594
|
0.66 (0.32–1.38)
|
0.272
|
|
|
A/A
|
–
|
–
|
–
|
–
|
|
Dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A – A/A
|
0.83 (0.48–1.44)
|
0.511
|
0.63 (0.30–1.30)
|
0.213
|
|
Recessive
|
G/G – G/A
|
1.00
|
|
1.00
|
|
|
|
A/A
|
–
|
–
|
–
|
–
|
|
Log–additive
|
–
|
0.81 (0.47–1.38)
|
0.438
|
0.61 (0.30–1.24)
|
0.174
|
rs6531668
|
Co–dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A
|
1.12 (0.80–1.58)
|
0.507
|
1.41 (0.92–2.18)
|
0.116
|
|
|
A/A
|
1.57 (0.82–3.00)
|
0.173
|
1.76 (0.69–4.48)
|
0.237
|
|
Dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/A – A/A
|
1.19 (0.86–1.64)
|
0.291
|
1.45 (0.96–2.20)
|
0.077
|
|
Recessive
|
G/G – G/A
|
1.00
|
|
1.00
|
|
|
|
A/A
|
1.51 (0.80–2.85)
|
0.206
|
1.54 (0.61–3.88)
|
0.357
|
|
Log–additive
|
–
|
1.19 (0.92–1.54)
|
0.179
|
1.37 (0.97–1.94)
|
0.071
|
rs12644563
|
Co–dominant
|
T/T
|
1.00
|
|
1.00
|
|
|
|
T/C
|
0.89 (0.51–1.56)
|
0.693
|
0.70 (0.33–1.48)
|
0.351
|
|
|
C/C
|
–
|
–
|
–
|
–
|
|
Dominant
|
T/T
|
1.00
|
|
1.00
|
|
|
|
T/C–C/C
|
0.83 (0.48–1.44)
|
0.513
|
0.67 (0.32–1.39)
|
0.278
|
|
Recessive
|
T/T– T/C
|
1.00
|
|
1.00
|
|
|
|
C/C
|
–
|
–
|
–
|
–
|
|
Log–additive
|
–
|
0.79 (0.47–1.33)
|
0.373
|
0.65 (0.32–1.31)
|
0.227
|
rs12640631
|
Co–dominant
|
T/T
|
1.00
|
|
1.00
|
|
|
|
T/C
|
1.18 (0.75–1.85)
|
0.479
|
1.11 (0.62–1.98)
|
0.722
|
|
|
C/C
|
0.62 (0.10–3.78)
|
0.607
|
3.06 (0.31–2.73)
|
0.336
|
|
Dominant
|
T/T
|
1.00
|
|
1.00
|
|
|
|
T/C–C/C
|
1.14 (0.73–1.77)
|
0.566
|
1.18 (0.68–2.07)
|
0.557
|
|
Recessive
|
T/T– T/C
|
1.00
|
|
1.00
|
|
|
|
C/C
|
0.61 (0.10–3.70)
|
0.591
|
3.01 (0.31–2.25)
|
0.342
|
|
Log–additive
|
–
|
1.09 (0.72–1.64)
|
0.686
|
1.23 (0.74–2.04)
|
0.431
|
rs12642845
|
Co–dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/C
|
0.71 (0.51–0.99)
|
0.042*
|
0.90 (0.59–1.38)
|
0.627
|
|
|
C/C
|
1.34 (0.64–2.79)
|
0.439
|
0.70 (0.27–1.80)
|
0.455
|
|
Dominant
|
G/G
|
1.00
|
|
1.00
|
|
|
|
G/C–CC
|
0.77 (0.56–1.05)
|
0.102
|
0.87 (0.58–1.32)
|
0.514
|
|
Recessive
|
G/G–G/C
|
1.00
|
|
1.00
|
|
|
|
C/C
|
1.52 (0.74–3.14)
|
0.257
|
0.73 (0.28–1.85)
|
0.501
|
|
Log–additive
|
–
|
0.88 (0.68–1.15)
|
0.350
|
0.87 (0.62–1.23)
|
0.428
|
Table 6
Haplotype analysis of TLR6 SNPs in association with the risk of TB.
Haplotype
|
Freq
|
|
|
Cases
|
Controls
|
OR (95% CI)
|
Pa
|
TTA
|
0.87
|
0.88
|
0.87 (0.66–1.13)
|
0.284
|
CTA
|
0.74
|
0.91
|
0.74 (0.55–0.99)
|
0.043*
|
CCG
|
0.90
|
0.33
|
0.90 (0.74–1.09)
|
0.282
|
CTG
|
1.08
|
0.54
|
1.08 (0.90–1.28)
|
0.418
|
TTC
|
0.87
|
0.24
|
0.87 (0.71–1.07)
|
0.192
|
CCG
|
0.73
|
0.05
|
0.73 (0.48–1.11)
|
0.144
|
TCG
|
0.48
|
0.98
|
0.47 (0.28–0.79)
|
0.004*
|
TTG
|
0.93
|
0.31
|
0.93 (0.77–1.13)
|
0.479
|