Introduction
Tuberculosis (TB) is a chronic infectious disease, and it was reported that 10.0 million new cases were detected and 1.2 million deaths from TB occurred among HIV-negative people in 2018 [1]. It has been estimated that one-third of the whole population is infected with Mycobacterium tuberculosis (MTB); 10% of these will develop TB while the other 90% may never develop symptomatic disease. Previous studies revealed that TB is a complex disease that is affected by many factors, including genetic variants, nutrition status and environmental factors, as well as interactions between the host and pathogen.
As a disease that is caused by a pathogen, pulmonary tuberculosis (PTB) is partly regarded as an immune-related disease. The incidence and development of PTB are associated with host immune status, which is mostly decided by host immune regulation-related genes [Reference Verway2], especially those responsible for cellular immunity, such as Toll-like receptor (TLR) family genes (which play roles in the recognition of exogenous material), as well as immune regulation-related cytokine family genes [Reference Montoya3], inflammation response mediators, etc. Thus, host susceptibility studies mainly focused on the related genes and their interactions with the pathogen in recent years. As vital regulators for MTB recognition and immune activation, TLRs are regarded as critical elements for TB infection. Previous studies revealed that single nucleotide polymorphisms (SNPs) in TLR family genes were associated with TB risk [Reference Fouad4–Reference Sinha6], and downstream genes were implicated too, such as tumour necrosis factor [Reference Joshi7], interleukins [Reference Qi H8], etc.
Transporter-associated antigen processing (TAP) is responsible for dealing with exogenous pathogenic microorganisms and then handling and presenting them to immune cells [Reference Harriff9]. Consequently, some mutations located within key areas of TAP genes can lead to MHC class I deficiency and the reduction of CD8+ T cells. Some genetic variants within TAP genes were found to be associated with cervical intraepithelial neoplasia [Reference Natter10] as well as haematological malignancies [Reference Ozbas-Gerceker11]. In recent years, Sunder et al. found that the SNPs at position 637 of TAP1 were associated with HIV-TB co-infection [Reference Sunder12]. Other studies were performed on TAP polymorphisms and TB in various populations, including Koreans, Chinese, etc. [Reference Fouad4, Reference Wang13, Reference Gomez14]; however, the study conclusions were mostly inconsistent. Notably, the associations among the Chinese population needed to be clarified, and the present study aims to elucidate and verify the potential correlations.
Previous studies revealed that age is a vital factor of TB incidence, and ageing may lead to a continuous increase in the risk of TB [Reference Cui15]. Gender is considered to be another factor implicated in TB incidence, and TB is more common in males than females [Reference Wanli Kang16]. In order to clarify the impact of genetic factors on TB risk in different age and gender groups, the present study performed an analysis of the stratification to provide more clues for further study.
Methods
Ethics statement
The study was approved by the ethics committee of Zhejiang Provincial Center for Disease Control and Prevention, and written informed consents were obtained.
Study population
A population-based case–control study in Changshan and Jiangshan counties in Zhejiang province was conducted. In the present study, 168 newly diagnosed smear-positive PTB cases were recruited from 1 January 2016 to 30 June 2017. Meanwhile, 251 clinically diagnosed non-TB patients or healthy controls were enrolled in the study from the same hospital, and all the controls were verified to never previously have TB.
All recruited subjects were instructed to complete a structured questionnaire with demographic information as well as life style-related questions including smoking, drinking, dietary habits, etc. Health-related information of the recruited subjects was also collected.
Genotyping
Genotyping of SNPs at rs1057141and rs1135216 was analysed by polymerase chain reaction-restriction fragment length polymorphism. The primer sequences and the restriction endonuclease utilised, as well as their product sizes, are shown in Table 1.
Primers adopted for genotyping of the three SNPs
Each 10 μl PCR sample contained 30 ng of the DNA template, 5 μl of 2X Prime Taq Premix (Genet Bio, Korea), 2 mM of each primer and the appropriate amount of ddH2O. The PCR cycling parameters were set as follows: DNA denaturation at 95 °C for 5 min followed by 35 cycles of annealing at 95 °C for 30 s and extension at 72 °C for 30 s, with a final extension of 72 °C for 10 min. Three μl of PCR product was incubated at 37 °C with 3 U of each endonuclease, 4 μl of 2X buffer and the appropriate amount of ddH2O, and was then separated onto 2.5% agarose gels. The SNP of rs2228396 was detected by the MassARRAY® MALDI-TOF System. The primer sequences employed are shown in Table 1. Five per cent of the DNA samples were randomly selected and rechecked with the same method.
Statistical analysis
In the present study, SPSS software version 19.0 was employed, and the χ 2 test was used to compare genotype distribution between control and case groups, as well as the distribution of demographic factors, including age, gender, educational status, occupation, family income, etc. Non-conditional logistic regression was adopted and the odds ratios (ORs) and 95% confidence intervals (95% CIs) were computed to estimate the contributions of some alleles to PTB incidence. Haplotype analysis including two and three SNPs was performed with online software (http://analysis.bio-x.cn/myAnalysis.php), and frequencies <0.03 were disregarded during analysis. The associations of the selected haplotypes with PTB incidence were estimated with ORs and 95% CIs. A two-sided P value <0.05 was considered to be significant.
Results
Distribution of demographic factors in case and control groups
A total of 168 PTB cases and 251 controls were recruited in the study. As shown in Table 2, a significant difference was detected in the distribution of age between the case and control groups (P = 0.030). Differences in education level and family income between the case and control groups were both significant, and the P values were 0.001 and 0.014, respectively. No significant differences were detected in sex, occupations, or nationalities between case and control groups.
Demographic characteristics of recruited subjects
Associations of SNP with PTB
Adjusted associations of three loci (rs1135216, rs1057141 and rs2228396) with PTB risk for age, sex and family income are shown in Table 3. No significant association was detected in the three SNPs.
Associations of the SNPs with PTB
a Adjusted for age, sex and family income.
After stratification by age, the associations adjusting for sex and family income are shown in Table 4. A marginal significance was detected among those aged <60 years regarding the association of the rs1135216 AG/GG genotype with PTB risk (OR = 1.97, 95% CI 0.93–4.19, P = 0.077). No similar association was detected in other SNPs among the population aged <60 years.
Associations of the SNPs with PTB among different age groups
a Adjusted for sex and family income.
Among those aged ≥60 years, a significant association was detected between the rs1057141 GG genotype and increased PTB risk (OR = 0.17, 95% CI 0.04–0.79, P = 0.024); however, no significant association was observed in the SNPs at rs1135216 and rs2228396.
As shown in Table 5, no significant association was detected among males or females in any of the three SNPs after stratification by sex.
Associations of the SNPs with PTB among male and female population
a Adjusted for age and family income.
Associations of the selected haplotypes with PTB incidence
While haplotypes combined with two SNPs for analysis, the results are shown in Table 6. The results revealed that those carrying the haplotype AT at rs1057141 and rs2228396 have a significantly increased risk of PTB (OR = 2.83, 95% CI 1.30–6.14). No significant association was detected among the other haplotypes including the two SNPs.
Associations of the haplotypes including two SNPs with PTB risk
As shown in Table 7, when all three SNPs were included for analysis, the haplotype AAT at rs1057141/rs1135216/rs2228396 was significantly associated with increased PTB risk (OR = 2.89, 95% CI 1.34–6.27). No significant association was detected among the other haplotypes.
Associations of the haplotypes including three SNPs with PTB risk
Discussion
TB is a common chronic infectious disease, which is influenced by immune status. In the present study, three SNPs located within immune regulation-related genes TAP1/2 were selected, and their contributions to PTB risk were analysed. A marginal association was detected between the rs1057141 GG genotype and PTB risk. After stratification by age, we found that the rs1057141 GG genotype carriers showed a reduced risk to PTB among those aged ≥60 years, and a marginal association was observed between the rs1135216 GG/AG genotype and PTB risk among those aged <60 years. Carriers with the haplotype AT at rs1057141 and rs2228396 or the haplotype AAT at rs1057141, rs1135216 and rs2228396 have a higher susceptibility to PTB.
TAPs are a key class of proteins that translocate antigenic peptides into the endoplasmic reticulum lumen for loading onto MHC class I molecules [Reference Harriff9]. Since the TAP gene is closely related to antigen processing and immune regulation, the current research involving TAP mainly focuses on immune-related diseases, including one study that investigated the association of the SNP at rs1135216 of TAP1 with the onset of allergic rhinitis [Reference Liu, Chen and Qi17], and another study which found that the SNP at rs1057141 was associated with an increased risk of atopic disease among African populations [Reference Liu, Chen and Qi17]. However, no significant association was detected between either rs1057141 or rs1135216 and allergic rhinitis among the Xinjiang Han Chinese population [Reference Chen and Zhang18]. These studies indicate that TAP-related disease incidence might be affected by ethnic characteristics and other factors.
TB is a chronic infectious disease that is considered to be affected by demographic characteristics that might cause the decline of immune status, as well as genetic factors. Recent studies revealed that rs4148871, rs4148876 and rs2857103 in TAP2 were associated with PTB risk, as well as the SNPs at rs4148871 and rs2857103 [Reference Thu19]. Another study conducted among the Iranian population revealed that the G allele at rs1135216 in TAP1 was significantly associated with an increased risk of TB, while the SNP at rs241447 in TAP2 was associated with a reduced risk of TB [Reference Naderi and Amininia20]. Additionally, another study conducted in South Korea also confirmed that the TAP2 gene was associated with the occurrence, development and even recurrence of TB [Reference Roh21]. However, no significant association was detected in another study [Reference Gomez14].
An accumulation of previous studies revealed that age is a vital factor associated with TB incidence [Reference Fernandes22], and our present study also revealed this as an important factor, as the difference between the case and control group was significant in our study. Family income was also found to be a factor. In addition, sex was also found to be implicated in TB risk, and TB was reported to be more common in males than females, although no significance was detected in our study between the case and control groups. Thus, during statistical analysis, age, sex and family income were included as covariates in the logistic regression model to discover the real associations of genetic variants with TB.
Stratification is considered to be an effective tool to eliminate the influence of certain factors, so stratifications by sex and age were conducted in our study to find out potential associations. Our present study found that rs1135216 and rs1057141 showed no effect on PTB risk among the whole study population. However, the GG genotype at rs1057141 was associated with a reduced PTB risk among those aged 60 or more years after stratification by age. Just as previous studies reported, ageing is a key process responsible for the decline of body immunity [Reference Guzzetta and Kirschner23], and the immune system might be disabled by the dysfunction of innate inactivation [Reference Bailey24]. Thus, we speculate that the association of SNP at rs1057141 with PTB might be induced by age-related immune system alteration. Immunoregulation of the body is a complex process that involves the participation of many factors, including intracellular signalling mediators, which may be involved in age-related alterations in neuroendocrine-immune interactions [Reference Muller, Di Benedetto and Pawelec25]. As a result, a significant protective effect was detected among those with a reduced immune system due to increased age, while no significant association was detected among the younger population with better immune conditions. Wang et al. also revealed that those aged 51 years and older had a higher risk of developing TB [Reference Wang26].
Haplotype combination and analysis is a useful statistical analysis strategy that combines multi-loci together to discover genetic predictors of disease incidence risk [Reference Gandhi27, Reference Dumaidi28]. Our study found that both haplotypes, including AT at rs1057141 and rs2228396 and AAT at rs1057141, rs1135216 and rs2228396, were significantly associated with increased PTB risk. The results revealed that the combination of multi-site genetic variants was a useful and valuable tool to provide deep insight and more clues to discover disease incidence risk.
Our study identified that the SNP of rs1057141 was associated with decreased PTB risk among the older population, and the SNP was validated to be a functional genetic variant that could induce functional change. However, a previous study revealed that the effect could be affected by age, sex and family income, etc. Considering the distribution of the factors, we performed the analysis by adjusting for the potential factors, and consequently measured the association by stratification.
Notably, sample size might be a limit of the present study, so it will be necessary to validate the real effect of the SNPs with a study evaluating a larger sample size.
Data
The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.
Acknowledgement
The study was supported by the projects of Study on Molecular Mechanism and Risk Prediction of Pulmonary Tuberculosis Based on microRNA Regulation (2020keyan512) and Study on Spectrum and Epidemic Characteristics of Non-tuberculous Mycobacterium in Zhejiang Province (2017keyan284).
Conflict of interest
None.
Open access
