Association between TAP gene polymorphisms and tuberculosis susceptibility in a Han Chinese population in Guangdong.

Tuberculosis (TB) is an important public health problem. Studies indicated that TAP plays a key role in the presentation and transport of antigenic peptides during anti-M.tb infection. Given the important biological role of the TAP gene involved in anti-M.tb infection, a family-based case-control study including 133 tuberculosis patients, 107 healthy household contacts, and 173 healthy controls was conducted to assess the association between TAP gene polymorphisms and TB susceptibility. The basic information of subjects and their blood samples were collected. Four SNPs including rs1135216, rs1057141, rs241447, and rs3819721 were genotyped by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Our results suggested that BMI, residence, bedroom crowding, indoor humidity, fitness activities, history of smoking, and TB exposure history were associated with the occurrence of tuberculosis (P < 0.05). A significant association was observed between the TAP1 rs1135216 CT/CC genotype and increased TB risk, and the ORs were 2.56 (95% CI 1.31-4.99) and 6.73 (95% CI 1.33-34.02), respectively. TAP2 rs3819721 GG genotype carriers also showed an increased risk of TB when compared TB patients to healthy household contacts. Haplotype analysis revealed that the haplotype CT at rs1057141 and rs1135216 (OR = 11.34, 95% CI 1.49-86.56; OR = 7.45, 95% CI 1.43-38.76), as well as TA at rs241447 and rs3819721 (OR = 2.20, 95% CI 1.07-4.56) had a significantly increased risk of TB. The genetic risk scores (GRS) analysis of the four loci indicated that the risk of tuberculosis increased with increasing GRS scores in TB vs HHC (Ptrend = 0.010) and in TB vs HC (Ptrend = 0.001). In conclusion, our findings suggested that the SNPs of rs1135216 and rs3819721 were associated with TB susceptibility among the tuberculosis-prone families in the Chinese Han population and the risk of developing tuberculosis increases with the number of risk alleles, which could help identify high-risk groups in time and take scientific preventive measures. Further cohort studies with large samples are needed to validate the role of TAP gene variants on TB susceptibility.

Introduction

Tuberculosis (TB), a chronic infectious disease caused by Mycobacterium tuberculosis (M.tb), is a serious threat to human health and a major public health problem in developing countries (Zumla et al. 2013). The World Health Organization (WHO) reported that there were 10 million new cases of TB in 2020 and China has the second highest TB burden in the world, accounting for 9% of the world’s TB cases (Chakaya et al. 2021). The Coronavirus Disease 2019 (COVID-19) ravages the world, and extraordinary health demands during the pandemic have diminished efforts to diagnose and treat people with tuberculosis (Kuehn 2020). While one-third of the world’s population were infected with M.tb according to the statistics, only 10% of infected individuals eventually develop active TB, the other 90% may never develop symptomatic disease (Tang et al. 2011; Moonan et al. 2020). This suggests that a large proportion of these individuals may be naturally resistant to TB infection and that host genetic susceptibility plays a non-negligible role in the development of TB. However, the specific role played by specific genetic factors in the development of tuberculosis is unknown, so further research into host genetic factors associated with M.tb infection could help us uncover new molecular markers and provide potential targets for tuberculosis prevention and treatment research.

Several studies support the influence of host genetics on TB susceptibility, including case observations, twin studies, candidate gene association and genome-wide association studies (Thye et al. 2010; Vannberg et al. 2011). Moreover, in recent years, there have been a large number of previous studies proving the evidence of association between host genetic polymorphisms and TB susceptibility (Selvaraj et al. 2015; Yuan et al. 2016). These genetic variants are mainly involved in the body’s immune response to M.tb and their polymorphisms may lead to increased genetic susceptibility to TB. Therefore, it is essential to determine the relationships between the genetic variations and the susceptibility to TB to comprehensively understand the pathogenesis of mycobacterial infection (Zhang et al. 2020).

TB is considered to be an immune-related disease. The incidence and development of TB are associated with host immune status, which is mostly decided by host immune regulation-related genes (Verway et al. 2013). Transporter-associated antigen processing (TAP) is responsible for transferring the endogenous antigen peptides from the cytoplasm to the endoplasmic reticulum (ER) and then handling and presenting them to immune cells (Harriff et al. 2013). TAP complex is a heterodimer consisted of two subunits encoded by TAP1 and TAP2 genes, a cluster of genes located within the HLA class II region of the HLA between the HLA-DP and -DQ loci on the short arm of human chromosome 6 (6p21.3) (Ramos et al. 2009). Single-nucleotide polymorphisms in the TAP1 and TAP2 genes may result in structural changes that prevent TAP heterodimer formation and alter antigen recognition and presentation, as evidenced in several experimental models, human case studies and malignancies (Zhang et al. 2002; Tang et al. 2001). These polymorphisms can serve as markers for human disease. Several studies have supported the association of TAP genes with various multifactorial diseases, such as immune-mediated, viral, and bacterial diseases and various types of cancers (Ozbas-Gerceker et al. 2013; Abitew et al. 2020; Qian et al. 2017; Thu et al. 2016). Studies were also performed on TAP polymorphisms and TB in various populations, including Iran, Koreans, Chinese, etc. (Naderi et al. 2016; Roh et al. 2015; Zhang et al. 2021). However, the research conclusions are not always consistent. Notably, the research in Han Chinese population into the relationship between TAP gene polymorphisms and TB is limited.

The risk of M.tb infection is much higher in household contacts of tuberculosis patients, ranging from 30 to 80% (Marks et al. 2000). In the present study, we intend to use family-based case–control study to avoid the false-positive results that might occur in traditional case–control studies due to population stratification. The aim of this study was to investigate the associations between polymorphisms of TAP gene and susceptibility to TB in the Han population in Guangdong province. Concurrently, we also attempted to determine whether the haplotypes covering these SNPs were linked to the development of TB and elucidate the immunopathological mechanisms of tuberculosis and provide a scientific basis for the prevention and treatment of tuberculosis.

Materials and methods

Ethics statement

All individuals voluntarily participated in the study, signed the informed consent before conducting the epidemiological investigation and collecting peripheral blood samples. All experiments were performed in accordance with relevant named guidelines and regulations and approved by the audit of the Ethics Committee of the First Affiliated Hospital of Guangdong Pharmaceutical University.

Subjects

Subjects were recruited from January 2017 to January 2020 among TB-prone families and healthy control families in the Han population in Guangdong province. A total of 413 people were enrolled in the study, including 133 TB patients, 107 genetically related healthy household controls contacts (HHCs), and 173 healthy controls (HCs).

According to one of the questions in Guangdong Province TB Drug Resistance Surveillance Questionnaire: ‘does any resident of your household suffer from tuberculosis’, we documented and verified familial clustering (families with two or more patients) from 2017 to 2020. Then we identified the probands (the first patient) and their blood family members and recruited them into the TB-prone family group.

The TB diagnosis was based on the Guidelines for The Implementation of Tuberculosis Control Program in China: (1) clinical symptoms suggestive of TB; (2) evidence of active TB in the upper lobes of the chest radiograph; (3) sputum smears positive with or without chest radiographic evidence of active TB; (4) excluding the individuals with diseases such as lung carcinoma, pneumonia, diabetes, HIV/AIDS and those under immunosuppressive treatment.

The healthy control families were from the general population with the same socioeconomic status, family structure and ethnic background as that of the TB-prone families, and without blood ties. Definition: (1) HHC (healthy household contacts): individuals who have a kinship with the proband, residing in the same house. (2) HC (healthy controls): individuals without exposure to M.tb. Considering that age is a potential confounding variable, the proband of each control family was matched to the proband of a tuberculosis-prone family, the matching criterion was an age difference no higher or smaller than five years. Other demographic characteristics are not matched to avoid information loss due to over-matching. All subjects had no kinship with the members in the TB-prone families and did not suffer from tuberculosis, HIV infection and diabetes.

Measurement of risk factors

There are numerous potential environmental risk factors for tuberculosis. Some of the most studied factors include the history of smoking, alcohol consumption, BCG vaccination, crowded housing and exposure to TB patients. BMI, fitness activities, socioeconomic conditions, indoor humidity, lighting and hygiene may also influence the development of TB, and studies were relatively limited. This study was based on the literature to comprehensively examine the environmental factors that may influence the development of tuberculosis and then analyses the true link between the TAP gene and tuberculosis while controlling for environmental factors. Survey data was collected using epidemiological questionnaires and face-to-face field interviews. The questionnaire included general demographic characteristics, such as gender, age, occupation, education level, Bacillus Calmette Guerin (BCG) vaccination history, and personal behavior, such as smoking status, alcohol consumption and fitness situation, as well as living environmental factors, such as lighting, humidity, ventilation and bedroom congestion. The definition and judgment criteria of related indicators as follows: smoking was defined as having at least one cigarette a day for 3 months, and drinking was defined as having at least three times a week for 1 year. At least three times a week, each time more than 30 min of fitness activities and lasting more than half a year is considered to have fitness, as a measure of nutritional status, a BMI of less than 18.5 is defined as low weight malnutrition according to the WHO.

DNA extraction and genotyping

Approximately, 5.0 mL of peripheral blood was collected from participants into a tube containing EDTA (ethylenediaminetetraacetic acid) and then stored in a low temperature refrigerator at − 80 ℃. Genomic DNA was isolated from peripheral blood cells using the E.Z.N.A. SE Blood DNA Kit D3471-02 (Omega Bio-Tek) according to the manufacturers’ protocol. DNA concentrations were determined using the NanoDrop-2000 spectrophotometer (Thermo Scientific).

Genotyping of SNPs at rs1057141 and rs1135216 of TAP1 gene and rs241447 and rs3819721 of TAP2 was analyzed by polymerase chain reaction-restriction fragment length polymorphism (PCR–RFLP). Each 5 µL PCR sample contained 0.25 µL TaqMan probe (Thermo Fisher Scientific), 2.5 µL Mix (Takara Biotechnology), 0.5 µL DNA, and 1.75 µL ultra-pure water. The mixture was added to a 384-well plate and placed in a PCR apparatus (ABI Prism® 7900HT) for automatic genotyping. The PCR cycling parameters were set as follows: pre-denaturation at 95 ℃ for 30 s, denaturation at 95 ℃ for 5 s, annealing extension at 60 ℃ for 30 s, and 40 cycles were repeated.

Blank control and 5% duplicate control were set on each reaction plate in the experiment, and the consistent typing rate of repeated samples reached 100%.

Statistical analysis

SPSS 26.0 software was employed for statistical analysis of the data. Continuous data are reported as mean ± standard deviation. Categorical variables are shown as number and proportion. Independent sample t test for continuous data and Chi-square test for categorical data were used. Hardy–Weinberg Equilibrium (HWE) was tested using Chi-square tests for each SNP. The associations between genotypes and TB incidence were calculated by computing the odds ratio (OR) and 95% confidence intervals (95% CI) from non-conditional logistic regression analyses. Linkage disequilibrium (LD) and haplotype analysis for four SNPs was performed with online software (http://analysis.bio-x.cn/myAnalysis.php) (Li et al. 2009), and frequencies < 0.03 were disregarded during analysis. The association of genetic risk scores (GRS) for the four loci of the TAP gene with the development of TB was analyzed by constructing simple genetic risk scores and weighted genetic risk scores. P value less than 0.05 was considered statistically significant.

Results

Demographic characteristics

The demographic characteristics of study participants are shown in Table 1. A total of 170 families were enrolled in this study, 85 tuberculosis-prone families and 85 healthy control families, which included 133 TB patients, and their 107 household contacts and 173 healthy controls, with an average age of 43.65 ± 16.06, 34.11 ± 19.71, and 41.90 ± 14.75 years, respectively.

Table 1

Demographic characteristics of TB (patients), HHC (healthy household contacts) and HC (healthy controls)

Variable	TB, n (%)	HHC, n (%)	HC, n (%)	TB vs HHC		TB vs HC
				χ2	P value	χ2	P value
Age				16.214	0.001	1.12	0.773
≤ 24	17 (12.8)	35 (32.7)	23 (13.3)
25–34	25 (18.8)	20 (18.7)	39 (22.5)
35–44	18 (13.5)	15 (14.0)	26 (15.0)
≥ 45	73 (54.9)	37(34.6)	85 (49.1)
Gender				14.56	< 0.001	19.91	< 0.001
Male	91 (68.4)	47 (43.9)	74 (42.8)
Female	42 (31.6)	60 (56.1)	99 (57.2)
Level of education				2.00	0.367	47.44	< 0.001
Junior and below	93 (69.9)	83 (77.6)	62 (35.8)
Senior high school	29 (21.8)	16 (15.0)	40 (23.1)
Undergraduate and above	11 (8.3)	8 (7.5)	71 (41.0)
Occupations				20.44	0.001	48.12	< 0.001
Students	4 (3.0)	22 (20.6)	9 (5.2)
Business	18 (13.5)	11 (10.3)	25 (14.5)
Public officials	11 (8.3)	8 (7.5)	64 (37.0)
Workers	33 (24.8)	16 (15.0)	23 (13.3)
Farmers	25 (18.8)	18 (16.8)	7 (4.0)
Others	42 (31.6)	32 (29.9)	45 (26.0)

The distribution of age, gender and occupations were significantly different between TB and HHC (P < 0.05). A significant difference was achieved in the distribution of TB and HC regarding gender and level of education and occupations (P < 0.05), while the distribution of level of education between TB and HHC, as well as the distribution of age between TB and HC, had no significant difference.

Environmental risk factors

Environmental risk factors’ distribution between TB and HHC

As shown in Table 2, a higher proportion of smokers (50.4%), non-exercisers (68.4%) and individuals with TB exposure history (90.2%) was observed among the TB compared with HHC (15.0, 52.3 and 79.4%, respectively), and the differences were statistically significant (P < 0.05). However, no significant difference was detected when comparing BMI, medical insurance, monthly household income and BCG vaccine history between the two groups. For the living environment factors, as TB and HHC are in the same household and have the same living environment, no comparison is made.

Table 2

Comparison of environmental factors between TB (patients) and HHC (healthy household contacts)

Variable	Total	TB, n (%)	HHC, n (%)	χ2	P value
BMI (kg/m2)	240			0.95	0.330
≤ 18.5	60	30 (22.6)	30 (28.0)
> 18.5	180	103 (77.4)	77 (72.0)
Medical insurance	240			2.65	0.103
No	29	11 (8.3)	16 (15.0)
Yes	211	122 (91.7)	91 (85.0)
Monthly household income	240			1.07	0.586
≤ 3000	79	46 (34.6)	33 (30.8)
3001–7999	110	57 (42.9)	53 (49.5)
≥ 8000	51	30 (22.6)	21 (19.6)
Smoking status	240			32.89	< 0.001
Yes	83	67 (50.4)	16 (15.0)
No	157	66 (49.6)	91(85.0)
Drinking status	240			1.35	0.245
Yes	65	40 (30.1)	25 (23.4)
No	175	93 (69.9)	82 (76.6)
Fitness activities	240			6.46	0.011
Yes	93	42 (31.6)	51 (47.7)
No	147	91 (68.4)	56 (52.3)
BCG vaccine	240			1.50	0.221
Yes	149	78 (58.6)	71 (66.4)
No	91	55 (41.4)	36 (33.6)
History of TB exposure	240			5.54	0.019
Yes	225	120 (90.2)	85 (79.4)
No	15	13 (9.8)	22 (20.6)

BMI body mass index, BCG vaccine Bacillus Calmette Guerin vaccine

Environmental risk factors’ distribution between TB and HC

The comparison results of environmental factors between TB and HC are summarized in Table 3. The HC reported higher proportions of individuals with a BMI > 18.5 (91.9%), high monthly household income (43.9%), fitness activities (68.2%), urban residence (72.8%), good indoor sanitation (85.0%), good indoor lighting (79.2%) and a history of BCG vaccination (64.2%) than TB (P < 0.05). Nevertheless, the proportions of smokers (50.4%), individuals in crowded housing (80.5%), indoor moist (74.4%), and TB contact messengers (90.2%) were higher among TB than in HC, and the difference was statistically significant (P < 0.05).

Table 3

Comparison of environmental factors between TB (patients) and HC (healthy controls)

Variable	Total	TB, n (%)	HC, n (%)	χ2	P value
BMI (kg/m2)	306			12.78	< 0.001
≤ 18.5	44	30 (22.6)	14 (8.1)
> 18.5	262	103 (77.4)	159 (91.9)
Medical insurance	306			2.42	0.120
No	18	11 (8.3)	7(4.0)
Yes	288	122 (91.7)	166 (96.0)
Monthly household income	306			16.37	< 0.001
≤ 3000	82	46 (34.6)	36 (20.8)
3001–7999	118	57 (42.9)	61 (35.3)
≥ 8000	106	30 (22. 6)	76 (43.9)
Smoking status	306			39.46	< 0.001
Yes	101	67 (50.4)	29 (16.8)
No	205	66 (49.6)	144 (83.2)
Drinking status	306			2.48	0.116
Yes	107	40 (30.1)	67 (38.7)
No	199	93 (69.9)	106 (61.3)
Fitness activities	306			40.44	< 0.001
Yes	160	42 (31.6)	118 (68.2)
No	146	91 (68.4)	55 (31.8)
Residence	306			44.68	< 0.001
Rural	134	87 (65.4)	47 (27.2)
Towns	172	46 (34.6)	126 (72.8)
Bedroom occupants	306			17.62	< 0.001
≤ 1 person/room	99	26 (19.5)	73 (42.2)
≥ 2 person/room	207	107 (80.5)	100 (57.8)
Indoor ventilation	306			1.70	0.192
Everyday	270	121 (91.0)	149 (86.1)
Definite time	36	12 (9.0)	24 (13.9)
Indoor humidity	306			4.51	0.034
Dry	208	34 (25.6)	64 (37.0)
Moist	98	99 (74.4)	109 (63.0)
Indoor lighting	306			14.19	< 0.001
Good	216	79 (59.4)	137 (79.2)
General	90	54 (40.6)	36 (20.8)
Indoor sanitation	306			15.99	< 0.001
Good	234	87 (65.4)	147 (85.0)
General	72	46 (34.6)	26 (15.0)
BCG vaccine	306			15.76	< 0.001
Yes	166	55 (41.4)	111 (64.2)
No	140	78 (58.6)	62 (35.8)
History of TB exposure	306			79.11	< 0.001
Yes	190	120 (90.2)	70 (17.9)
No	116	13 (9.8)	103 (82.1)

BMI body mass index, BCG vaccine Bacillus Calmette Guerin vaccine

Multivariate logistic regression analysis

Table 4 shows the results of multivariate logistic regression analysis for TB vs HHC and TB vs HC. When comparing TB to HHC, the risk of tuberculosis in smokers is 3.61 times higher than in non-smokers (P = 0.003, OR = 3.61, 95% CI 1.54–8.45). Other statistically significant difference included the history of TB exposure (P = 0.014, OR = 2.99, 95% CI 1.25–7.16). Fitness activities can effectively reduce the risk of tuberculosis (P = 0.038, OR = 0.58, 95% CI 0.28–0.96).

Table 4

Multivariate logistic regression analysis of risk factors regarding TB vs HHC and TB vs HC

Group	Variable	β	Standard deviation	χ2	P value	OR	95% CI
TB vs HHC	Age
	≤ 24					1.00
	25–34	0.14	0.55	0.07	0.796	1.15	0.39–3.38
	35–44	0.41	0.59	0.47	0.493	1.50	0.47–4.79
	≥ 45	1.01	0.52	3.68	0.055	2.73	0.98–7.63
	Gender
	Male					1.00
	Female	− 0.56	0.39	2.06	0.151	0.57	0.27–1.23
	Occupations
	Students					1.00
	Business	1.58	0.82	3.77	0.052	4.88	0.99–24.11
	Public officials	1.10	0.90	1.49	0.223	3.00	0.51–17.55
	Workers	1.45	0.75	3.72	0.054	4.27	0.98–18.67
	Farmers	1.28	0.82	2.40	0.121	3.59	0.71–18.02
	Others	1.25	0.73	2.94	0.086	3.49	0.84–14.56
	Smoking status
	No					1.00
	Yes	1.28	0.43	8.75	0.003	3.61	1.54–8.45
	Fitness activities
	No					1.00
	Yes	− 0.66	0.32	4.30	0.038	0.52	0.28–0.96
	History of TB exposure
	No					1.00
	Yes	1.09	0.45	6.02	0.014	2.99	1.25–7.16
TB vs HC	Gender
	Male					1.00
	Female	− 0.56	0.48	1.36	0.244	0.57	0.23–1.46
	Level of education
	Junior and below					1.00
	Senior high school	0.27	0.46	0.33	0.563	1.30	0.53–3.21
	Undergraduate and above	− 1.10	0.57	3.77	0.052	0.33	0.11–1.01
	Occupations
	Students					1.00
	Business	0.17	1.07	0.02	0.877	1.18	0.15–9.57
	Public officials	− 0.23	1.04	0.05	0.824	0.79	0.10–6.13
	Workers	0.68	1.01	0.46	0.498	1.98	0.28–14.23
	Farmers	1.56	1.12	1.93	0.165	4.75	0.53–42.77
	Others	0.16	1.00	0.03	0.871	1.18	0.17–8.34
	BMI (kg/m2)
	≤ 18.5					1.00
	> 18.5	− 1.20	0.56	4.62	0.032	0.30	0.10–0.90
	Monthly household income
	≤ 3000					1.00
	3001–7999	0.22	0.46	0.23	0.630	1.25	0.51–3.09
	≥ 8000	− 0.16	0.51	0.09	0.762	0.86	0.31–2.34
	Smoking status
	No					1.00
	Yes	1.00	0.49	4.11	0.043	2.72	1.03–7.15
	Fitness activities
	No					1.00
	Yes	− 0.90	0.39	5.48	0.019	0.41	0.19–0.86
	Residence
	Rural					1.00
	towns	− 0.96	0.40	5.83	0.016	0.38	0.18–0.84
	Bedroom occupants
	≤ 1 person/room					1.00
	≥ 2 person/room	1.27	0.41	9.85	0.002	3.56	1.61–7.89
	Indoor humidity
	Dry					1.00
	Moist	1.01	0.46	4.95	0.026	2.75	1.13–6.70
	Indoor lighting
	Good					1.00
	General	0.48	0.43	1.21	0.271	1.61	0.69–3.75
	Indoor sanitation
	Good					1.00
	General	0.95	0.50	3.70	0.054	2.60	0.98–6.84
	BCG vaccine
	Yes					1.00
	No	− 0.60	0.37	2.65	0.104	0.55	0.27–1.13
	History of TB exposure
	No					1.00
	Yes	2.56	0.46	31.46	< 0.001	12.93	5.29–31.61

BMI body mass index, BCG vaccine Bacillus Calmette Guerin vaccine

Multivariate analysis of TB and HC showed that BMI higher than 18.5 (P = 0.032, OR = 0.30, 95% CI 0.10–0.90), fitness activities (P = 0.019, OR = 0.41, 95% CI 0.19–0.86) and living in towns (P = 0.016, OR = 0.38, 95% CI 0.18–0.84) were the protective factors for tuberculosis. While the history of smoking (P = 0.043, OR = 2.72, 95% CI 1.03–7.15), bedroom crowding (P = 0.002, OR = 3.56, 95% CI 1.61–7.89), indoor moist (P = 0.026, OR = 2.75, 95% CI 1.13– 6.70) and history of TB exposure (P < 0.001, OR = 12.93, 95% CI 5.29–31.61) increased the risk of tuberculosis.

Association between TAP gene polymorphisms and susceptibility to tuberculosis

Hardy–Weinberg equilibrium test

The genotype distributions of TAP1 rs1135216 and rs1057141, as well as TAP2 rs241447 and rs3819721 were found to be in Hardy–Weinberg equilibrium in TB, HHC and HC (P > 0.05). The results of the HWE check and MAF are shown in Supplementary Table 5.

Table 5

Genotype and allele distribution of TAP1 rs1135216 and rs1057141 in TB, HHC and HC

SNPs	Genetic model	Genotype/allele	TB, n (%)	HHC, n (%)	HC, n (%)	TB vs HHC		TB vs HC
						OR (95% CI)a	P value	OR (95% CI)b	P value
rs1135216	Co-dominant	TT	75 (56.4)	81 (75.7)	125 (72.2)	1.00	0.002	1.00	0.092
		CT	48 (36.1)	24 (22.4)	46 (26.6)	2.56 (1.31–4.99)		1.67 (0.82–3.43)
		CC	10 (7.5)	2 (1.9)	2 (1.2)	6.73 (1.33–34.02)		5.48 (0.76–39.53)
	Dominant	TT	75 (56.4)	81 (75.7)	125 (72.2)	1.00	0.001	1.00	0.072
		CT–CC	58 (43.6)	26 (24.3)	48 (27.8)	2.92 (1.55–5.53)		1.88 (0.94–3.77)
	Recessive	TT–CT	123 (92.5)	105 (98.1)	171 (98.8)	1.00	0.025	1.00	0.094
		CC	10 (7.5)	2 (1.9)	2 (1.2)	5.14 (1.04–25.53)		4.49 (0.63–31.84)
	Super-dominance	TT–CC	85 (63.9)	83 (77.6)	127 (73.4)	1.00	0.010	1.00	0.271
		CT	48 (36.1)	24 (22.4)	46 (26.6)	2.28 (1.1.18–4.39)		1.48 (0.73–2.99)
	Alleles	T	198 (74.4)	186 (86.9)	296 (85.5)	1.00	< 0.001	1.00	0.001
		C	68 (25.6)	28 (13.1)	50 (14.5)	2.28 (1.41–3.70)		2.03 (1.35–3.06)
rs1057141	Co-dominant	TT	70 (52.6)	50 (46.7)	103 (59.5)	1.00	0.930	1.00	0.732
		CT	53 (39.9)	50 (46.7)	59 (34.1)	0.89 (0.50–1.60)		1.29 (0.64–2.58)
		CC	10 (7.5)	7 (6.5)	11 (6.4)	0.98 (0.32–2.99)		0.89 (0.27–2.97)
	Dominant	TT	70 (52.6)	50 (46.7)	103 (59.5)	1.00	0.721	1.00	0.582
		CT–CC	63 (47.4)	57 (53.3)	70 (40.5)	0.90 (0.51–1.59)		1.20 (0.62–2.31)
	Recessive	TT–CT	123 (92.5)	100 (93.5)	162 (93.6)	1.00	0.962	1.00	0.723
		CC	10 (7.5)	7 (6.5)	11 (6.4)	1.03 (0.35–3.05)		0.81 (0.25–2.59)
	Super-dominance	TT–CC	80 (60.1)	57 (53.3)	114 (65.9)	1.00	0.734	1.00	0.441
		CT	53 (39.9)	50 (46.7)	59 (34.1)	0.89 (0.51–1.58)		1.31 (0.66–2.57)
	Alleles	T	193 (72.6)	150 (70.1)	265 (76.6)	1.00	0.553	1.00	0.254
		C	73 (27.4)	64 (29.9)	81 (23.4)	0.89 (0.60–1.32)		1.24 (0.86–1.79)
rs241447	Co-dominant	TT	63 (47.4)	42 (39.2)	83 (48.0)	1.00	0.171	1.00	0.531
		CT	58 (43.6)	54 (50.5)	72 (41.6)	0.58 (0.32–1.06)		1.41 (0.72–2.78)
		CC	12 (9.0)	11 (10.3)	18 (10.4)	0.56 (0.20–1.55)		0.87 (0.26–2.86)
	Dominant	TT	63 (47.4)	42 (39.2)	83 (48.0)	1.00	0.062	1.00	0.422
		CT–CC	70 (52.6)	65 (60.8)	90 (52.0)	0.58 (0.32–1.03)		1.30 (0.68–2.48)
	Recessive	TT–CT	121 (91.0)	96 (89.7)	155 (89.6)	1.00	0.552	1.00	0.600
		CC	12 (9.0)	11 (10.3)	18 (10.4)	0.75 (0.29–1.95)		0.74 (0.23–2.33)
	Super-dominance	TT–CC	75 (56.4)	53 (49.5)	101 (58.4)	1.00	0.134	1.00	0.270
		CT	58 (43.6)	54 (50.5)	72 (41.6)	0.65 (0.37–1.15)		1.44 (0.75–2.78)
	Alleles	T	184 (69.2)	138 (64.5)	238 (68.8)	1.00	0.278	1.00	0.768
		C	82 (30.8)	76 (35.5)	108 (31.2)	0.81 (0.55–1.19)		1.05 (0.75–1.48)
rs3819721	Co-dominant	GG	61 (45.9)	71 (66.3)	160 (61.3)	1.00	0.021	1.00	0.160
		AG	58 (43.6)	31 (29.0)	57 (33.0)	2.28 (1.22–4.27)		0.51 (0.24–1.09)
		AA	14 (10.5)	5 (4.7)	10 (5.7)	2.57 (0.72–9.16)		0.46 (0.11–2.05)
	Dominant	GG	61 (45.9)	71 (66.4)	106 (61.3)	1.00	0.006	1.00	0.058
		AG–AA	72 (54.1)	36 (33.6)	67 (38.7)	2.32 (1.27–4.24)		0.50 (0.24–1.03)
	Recessive	GG–AG	119 (89.5)	102 (95.3)	163 (94.2)	1.00	0.363	1.00	0.450
		AA	14 (10.5)	5 (4.7)	10 (5.8)	1.75 (0.51–6.03)		0.57 (0.13–2.47)
	Super-dominance	GG–AA	75 (56.4)	76 (71.0)	116 (67.0)	1.00	0.022	1.00	0.110
		AG	58 (43.6)	31 (29.0)	57 (33.0)	2.04 (1.11–3.76)		0.54 (0.25–1.15)
	Alleles	G	180 (67.7)	173 (80.8)	289 (83.5)	1.00	0.001	1.00	< 0.001
		A	86 (32.3)	41 (19.2)	57 (16.5)	2.02 (1.32–3.09)		2.42 (1.65–3.55)

aAdjusted by smoking history, fitness activities and history of TB exposure

bAdjusted by BMI, smoking status, fitness activities, residence, bedroom occupants, indoor humidity and history of TB exposure

Four genetic patterns of TAP1 and TAP2 genes

As shown in Table 5, a significant association was observed between TAP1 rs1135216 and susceptibility to tuberculosis in four genetic models (CT vs TT: OR = 2.56, 95% CI 1.31–4.99; CC vs TT: OR = 6.73, 95% CI 1.33–34.02; CT–CC vs TT: OR = 2.92, 95% CI 1.55–5.53; CC vs TT–CT: OR = 5.14, 95% CI 1.04–25.53; CT vs TT–CC: OR = 2.28, 95% CI 1.18–4.39) regarding TB vs HHC. However, there was no significant association for such locus when comparing TB and HC (P > 0.05). For TAP2 rs3819721, the AG genotype in the co-dominant model (AG vs GG: OR = 2.28, 95% CI 1.22–4.27), GG genotype in the dominant model (GG vs AG–AA: OR = 2.32, 95% CI 1.27–4.24), as well as AG genotype in the super-dominance model (AG vs GG–AA: OR = 2.04, 95% CI 1.11–3.76) was found to be associated to TB when compared to HHC.

Different gene models of TAP1 rs1057141 and TAP2 rs241447 were not associated with tuberculosis susceptibility either in TB vs HHC (P > 0.05) or in TB vs HC (P > 0.05).

The C allele of TAP1 rs1135216 was found to be positively associated with the risk of tuberculosis in TB and HHC (P < 0.001, OR = 2.28, 95% CI 1.41–3.70), as much as in TB and HC (P = 0.001, OR = 2.03, 95% CI 1.35–3.06). The A allele of TAP2 rs3819721 also increased the susceptibility to tuberculosis. No significant difference in the distribution of alleles of TAP1 rs1057141 and TAP2 rs241447 in TB vs HHC (P > 0.05) and in TB vs HC (P > 0.05).

Haplotype analysis of TAP1 and TAP2 genes

Haplotype refers to the linear arrangement of a group of related multiple SNP sites on a chromosome or in a certain region (Jakobsson et al. 2008). Studies have shown that haplotypes as genetic markers have higher polymorphism and statistical efficiency than single SNPs (Davidson 2000).

As shown in Table 6, the Linkage Disequilibrium (LD) analyses indicated that polymorphisms rs1135216 and rs1057141 in TB vs HHC (D′ = 0.55, r2 = 0.19) and TB vs HC (D′ = 0.61, r2 = 0.27) were not independent from another, whereas the strength of association between rs241447 and rs3819721 was observed to be low.

Table 6

Analysis of association between haplotype and pulmonary tuberculosis susceptibility

Haplotype	TB, n (%)	HHC, n (%)	HC, n (%)	TB vs HHC				TB vs HC
				OR (95% CI)a		P value		OR (95% CI)b		P value
TAP1
TT	164 (61.8)	149 (69.5)	259 (75.0)	1.00				1.00
TC	34 (12.6)	37 (17.4)	37 (10.6)	0.59 (0.31–1.14)		0.120		0.94 (0.42–2.11)		0.890
CC	39 (14.8)	27 (12.5)	44 (12.8)	1.70 (0.92–3.15)		0.090		1.40 (0.72–2.73)		0.320
CT	29 (10.8)	1 (0.6)	6 (1.6)	11.34 (1.49–86.56)		0.020		7.45 (1.43–38.76)		0.018
				D′	0.55	r2	0.19	D′	0.61	r2	0.27
TAP2
TG	129 (48.6)	114 (53.2)	(57.5)	1.00				1.00
CG	51 (19.1)	59 (27.6)	(20.2)	0.74 (0.40–1.38)		0.350		1.17 (0.63–2.17)		0.61
TA	55 (20.6)	24 (11.3)	(11.3)	2.20 (1.07–4.56)		0.034		0.53 (0.22–1.27)		0.16
CA	31 (11.7)	17 (7.9)	(11.0)	1.13 (0.45–2.81)		0.800		0.73 (0.29–1.85)		0.51
				D'	0.29	r2	0.04	D′	0.27	r2	0.02

aAdjusted by smoking history, fitness activities and history of TB exposure

bAdjusted by BMI, smoking status, fitness activities, residence, bedroom occupants, indoor humidity and history of TB exposure

Four haplotypes were formed at the TAP1 rs1135216 and rs1057141, which were TT, TC, CC and CT. Take TT haplotype as reference, CT haplotype was associated with an increased risk of tuberculosis in TB vs HC (P = 0.020, OR = 11.34, 95% CI 1.49–86.56) and in TB vs HHC (P = 0.018, OR = 7.45, 95% CI 1.43–38.76). No other haplotypes of TAP1 were found to be statistically associated with susceptibility to tuberculosis either in TB vs HC or in TB vs HHC.

The two loci of TAP2 gene rs241447 and rs3819721 form four haplotypes of TG, CG, TA, and CA. Compared with the TG haplotype, individuals with the TA haplotype have a 1.2-fold increased risk of tuberculosis regarding TB vs HHC (P = 0.034, OR = 2.20, 95% CI 1.07–4.56).

Genetic risk score of the TAP gene

Individual genetic risk scores ranged from a minimum of 0 to a maximum of 6. The subgroups were 0–2, 3–4 and 5–6 (Table 7). In TB vs HHC, the risk of TB increased 3.4-fold (P = 0.038, OR = 4.40, 95% CI 1.09–17.76) for individuals with a GRS score of 5–6 compared to those with a GRS score of 0–2, after adjusting for smoking history, fitness activities and history of TB exposure. Similarly, when comparing TB and HC, individuals with a GRS score of 3–4 had a higher risk of TB than those with a GRS score of 0–2 (P = 0.047, OR = 2.13, 95% CI 1.01–4.47). The results of the trend test showed that the risk of tuberculosis increased with increasing GRS scores in TB vs HHC (Ptrend = 0.010) and in TB vs HC (Ptrend = 0.001).

Table 7

The relationship between the GRS of TAP gene and tuberculosis susceptibility

GRSc	TB, n (%)	HHC, n (%)	HC, n (%)	TB vs HHC		TB vs HC
				OR (95% CI)a	P value	OR (95% CI)b	P value
0–2	77 (57.9)	77 (72.0)	131 (75.7)	1.00		1.00
3–4	44 (33.1)	27 (25.2)	31 (20.8)	1.74 (0.91–3.33)	0.096	2.13 (1.01–4.47)	0.047
5–6	12 (9.0)	3 (2.8)	6 (3.5)	4.40 (1.09–17.76)	0.038	3.09 (0.84–11.41)	0.090

aAdjusted by smoking history, fitness activities and history of TB exposure

bAdjusted by BMI, smoking status, fitness activities, residence, bedroom occupants, indoor humidity and history of TB exposure

cThe risk alleles at rs1135216, rs1057141, rs241447 and rs3819721 are C, C, C, and A

Discussion

In addition to socioeconomic factors such as living conditions and poverty, host genetic susceptibility may be another important risk factor for the development of tuberculosis. In this family-based case–control study, our results demonstrated that SNPs at rs1135216 of the TAP1 gene and rs3819721 of the TAP2 gene were significantly associated with an increased risk of TB in the Chinese Han population. The haplotype of rs1135216*C/rs1057141*T and rs3819721*T/rs241447*A was found to be associated with high susceptibility to TB. Moreover, the risk of developing tuberculosis increases with the number of risk alleles. These findings are important for the timely identification of people at high risk of developing TB and for elucidating the molecular mechanisms of TB development.

In current study, we found that the rate of smoking among TB patients was relatively higher, accounting for 50.4%. Smokers had three times the risk of developing TB compared to non-smokers. Smoking was shown to be a significant risk factor for developing TB and is consistent with other studies (Velen et al. 2020; Htet et al. 2018). The study also observed that keeping a certain frequency of regular fitness could improve the body’s immunity and BMI greater than 18.5 was an effective protective factor against tuberculosis. Several studies (Lönnroth et al. 2010) had shown that BMI could have an effect on the development of tuberculosis. Moreover, multivariate analysis indicated that living in rural areas and crowded bedrooms, as well as moist indoors may be a risk factor for the occurrence of tuberculosis. This may be due to the fact that rural areas are less hygienic than urban areas and M.tb survives more easily in moist environments. Individuals with a history of TB exposure are also at increased risk of developing TB. Literature (Horsburgh and Rubin 2011; Acuña-Villaorduña et al. 2018) published had demonstrated that close contact to infectious tuberculosis cases is considered the strongest predictor of M.tb infection, more constant exposures to infectious M.tb and subsequently to more risk of progression to tuberculosis, which is in accordance with the findings of this research.

The present study sought to identify four polymorphisms in the coding regions of TAP1 and TAP2 genes and observed all the polymorphisms in the subjects enrolled in the study. Significant association between TAP1 rs1135216 and TAP2 rs3819721 genes and tuberculosis was observed when compared TB patients and healthy household contacts. Individuals containing TAP1 CC homozygote and CT heterozygote were found to be strongly associated with tuberculosis infection (OR = 6.73, 2.56, respectively). Similar to the study by Wang et al. (2012) showed that rs1135216 exhibited a significant relation to tuberculosis susceptibility. However, the study carried out in a Northwestern Colombian population also failed to detect any association of TAP1 gene with TB disease (Gomez et al. 2006). In addition, the present study also showed that the AG genotype of TAP2 rs3819721 was more common in tuberculosis patients than that in healthy household contacts, suggesting that AG genotype significantly increased risk of tuberculosis infection. In contrast, no association has been found between the risk of tuberculosis and the polymorphisms at TAP1 rs1057141 and TAP2 rs241447, which is not fully concordant with the results of studies by Naderi et al. (2016) and Zhang et al. (2021) found that polymorphisms at the TAP2 rs241447 locus can reduce the risk of tuberculosis. This may be due to the different races and regions of the subjects.

Haplotype combination and analysis combines multiple-loci together to discover genetic predictors of disease incidence risk (Gandhi et al. 2020). In this study, the rs1057141 and rs241447 polymorphism was not found to be significantly associated with TB susceptibility. However, when the rs1135216 and rs1057141, as well as rs3819721 and rs241447 SNPs were combined in haplotypes, the rs1135216*C/rs1057141*T and rs3819721*T/rs241447*A were found to be significantly associated with increased tuberculosis risk, and the OR were 11.34 (95% CI 1.49–86.56) and 2.20 (95% CI 1.07–4.56), respectively. The haplotype CT and TA were genetic markers that may predict tuberculosis risk. Consequently, it seems more important to determine the association between polymorphism and phenotype based on the configuration of haplotypes compared with alleles and genotypes. The results revealed that the combination of multiple genetic variants is a useful and valuable tool, providing deeper insights and additional clues for discovering the risk of disease incidence.

The number of genetic variants affecting the development of complex diseases is large and genetic risk scores combine the weak effects of multiple genetic variants to improve the predictability of genetic variants for disease (Li et al. 2020). Compared with those carrying 0–2 risk alleles, those carrying 3–4 risk alleles and those carrying 5–6 risk alleles increased the risk of tuberculosis by 2.13 times and 4.40 times, respectively, with increased disease risk with higher scores. The four genetic polymorphisms of the TAP1 and TAP2 genes have a certain cumulative effect in the pathogenesis of tuberculosis.

There are still some limitations in the current study that should be addressed. First, the number of families in our group was still relatively small. Second, the age matching in the study resulted in a degree of selection bias. Third, this study focused on Han Chinese subjects in Guangdong, China, and the results may not be applicable to other populations. In the future, studies in multiple and large samples are needed to accurately evaluate the true role of TAP gene polymorphisms in the pathogenesis of tuberculosis.

In summary, our research identified that the SNPs of rs1135216 in the TAP1 gene and rs3819721 in the TAP2 gene were associated with TB susceptibility among the Chinese Han tuberculosis-prone families, and the risk of developing tuberculosis increases with the number of risk alleles. These findings add to our understanding of the possible impact of TAP gene variants on the development of tuberculosis.

Below is the link to the electronic supplementary material.

Supplementary file1 (DOCX 14 kb)