Back to Journals » International Journal of General Medicine » Volume 15

Association of Three SNPs Loci of Kelch-Like-ECH-Associated Protein 1 (Human) with Tuberculosis in Chinese Han Population

Authors Wu S, Ding X, Yang Q, Wang M , He JQ

Received 6 May 2022

Accepted for publication 27 July 2022

Published 1 August 2022 Volume 2022:15 Pages 6365—6372

DOI https://doi.org/10.2147/IJGM.S373555

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 3

Editor who approved publication: Dr Scott Fraser



Shouquan Wu,* Xiaojuan Ding,* Qianlan Yang, Minggui Wang, Jian-Qing He

Department of Respiratory and Critical Care Medicine, West China Hospital, Sichuan University, Chengdu, People’s Republic of China

*These authors contributed equally to this work

Correspondence: Jian-Qing He, Department of Respiratory and Critical Care Medicine, West China Hospital, Sichuan University, 37, Guo Xue Alley, Chengdu, 610041, People’s Republic of China, Tel +86-18980602293, Fax +86-028-85422571, Email [email protected]

Purpose: Progression from latent tuberculosis infection (LTBI) to pulmonary TB (PTB) was associated with genetic polymorphisms, but there were limited genetic polymorphism data on LTBI and PTB. We aimed at examining the association of KEAP1 gene polymorphisms with PTB and LTBI.
Patients and Methods: PTB patients and close contacts of PTB patients were recruited from West China Hospital of Sichuan University. After obtaining the patient’s consent, we draw 2– 5mL of blood from the patient’s peripheral vein. Tag-SNPs of KEAP1 were chosen according to previous studies. The genotyping was done by improved multiplex ligase detection reaction (iMLDR). We used logistic regression to assess the association of SNPs with LTBI/PTB, with sex and age as covariates.
Results: A total of 209 PTB patients, 201 LTBI, and 204 HCS were included in the present study. Three Tag-SNPs were included in this study. Significant association was found for KEAP1 rs1048290 between LTBI and HCS. Compared with the KEAP1 rs1048290 CC genotype, genotype GC had an 38% decreased risk for development LTBI (P = 0.043, OR = 0.62, 95% CI: 0.039– 0.98). We also found that SNPs in KEAP1 were significantly related to PTB compared to LTBI. Compared with the rs11545829G allele, allele A had an 30% decreased risk for development PTB (P = 0.034, OR = 0.70, 95% CI: 0.51– 0.97). We also found the rs11668429 polymorphism was related to PTB. Compared with TT, GT had a significantly increased risk of LTBI developing into PTB (P = 0.041, OR = 1.68, 95% CI: 1.02– 2.77).
Conclusion: Our study suggested that KEAP1 polymorphisms were significantly related to susceptibility to PTB and LTBI subjects.

Keywords: latent tuberculosis infection, pulmonary tuberculosis, SNPs, association study

Introduction

Tuberculosis (TB), caused by Mycobacterium TB, is a serious infectious disease with a mortality rate second only to COVID-19.1 TB is a serious infectious disease worldwide, causes 1.3 million deaths every year. In 2020, 9.87 million TB patients were diagnosed worldwide, with an incidence rate of 1.27/100,000. In China, the estimated number of new TB cases in 2020 is 842,000 second only to India.1 However, only 5–15% of infected people progress to TB.2 It is unclear why only a small number of infected people progress to active disease.

Previous researches have shown that host genetics, specific strains of Mycobacterium TB (MTB) and environmental factors may explain why TB incidence varies across specific ethnic groups, geographic regions, genders and age groups.3,4 It was suggested that different host gene polymorphisms were closely associated with TB.5 Our previous basic research also confirmed the role of gene polymorphism in TB.6 The study of various host genes associated with TB can improve the understanding of the pathogenesis TB and the development of prevention or treatment strategies.

MTB infection can lead to oxidative stress in the intracellular environment, which can hinder the progress of TB.7 Nuclear factor erythroid 2-related factor 2 (Nrf2) protein is the main mediator of many anti-oxidant pathway genes expression. Our previous studies have shown that Nrf2 is associated with TB.6 NRF2 is commonly present in the cytoplasm and forms a complex with Kelch-like Ech-associated protein 1 (KEAP1), which leads to inhibition of Nrf2 activity temporarily.8 KEAP1 is located on the short arm of chromosome 19. The C-terminal Kelch domain of Keap1 binds Nrf2 through the evolutionarily conserved kelch repeat sequence motif. As a sensor of ROS, KEAP1 plays an important role in protecting cells from oxidative damage.9 Study has shown that Keap1 is involved in the repair process of MTB infection in macrophages.10 Keap1 gene polymorphisms have been reported to be associated with multiple diseases.11–13 However, the association between KEAP1 gene polymorphisms and TB is currently lacking.

Therefore, in this study, we hypothesized that KEAP1 genetic polymorphisms might play an important role in TB susceptibility. To test this hypothesis, 3 tag single-nucleotide polymorphisms (Tag-SNPs) in KEAP1 were analyzed to determine the role of Tag-SNPs in pulmonary TB (PTB) and latent TB infection (LTBI) in Chinese population.

Materials and Methods

Inclusion of Objects

A total of 209 PTB patients and 405 close contacts of PTB patients were recruited from West China Hospital of Sichuan University. Sample size calculations were from our previous studies.14 All participants were genetically unrelated Chinese Han people, and they were over 16 years old. The diagnosis of PTB was trough sputum culture and/or lung lavage by bronchoscopy and/or lung percutaneous biopsy. Patients with complications of diabetes, hepatitis B and/or C, chronic obstructive pulmonary disease, immune-mediated disease and HIV were excluded.

The close contacts of PTB in this study included accompanying persons, colleagues, spouses, and staff of the TB ward. We followed up for one year, excluding those who progressed to TB. Finally, based on interferon gamma release assay results, chest imaging, and symptoms, the 405 close contacts were finally divided into LTBI and healthy control (HCS) at enrolment and at one-year follow-up. Definition of close contacts is based on our previously published research.14

After obtaining the patient’s consent, we draw 2–5mL of blood from the patient’s peripheral vein. The blood collection tubes contained ethylene diamine tetra acetic acid (EDTA) and the whole blood was subsequently stored in a −80°C freezer. DNA was extracted from whole blood using a genomic DNA purification kit (Axygen Scientific, Inc., Union City, CA, USA) according to the instructions. The DNA samples were stored at −80°C for further investigation. All participants signed an informed consent form. Parents or legal guardians of participants under the age of 18 signed the informed consent form. This study has passed ethical review. This study was approved by the ethical committee of the West China Hospital Institutional Review Board (ethics number: 2019;761). Our study complies with the Declaration of Helsinki.

Tag-SNPs (SNP with high linkage disequilibrium in a genome and can represent multiple SNPs) of KEAP1 were chosen according to previous studies.15

Genotyping

The genotyping was done by improved multiplex ligase detection reaction (iMLDR), and the Shanghai Genesky Biotechnology company provides technical support. 5% of the samples were repeated for verification. The genotyping steps are briefly described below. The iMLDR process begins with multiplex PCR amplification of the mutant region. The amplification products are purified by exonuclease and shrimp alkalase. The purified product is used in a double ligase reaction. During the reaction, each locus contains a 5` terminal allele-specific probe and an immediately following 3` terminal site-specific probe. A specific ligating sequence is added to the 5` termina of each allele-specific primer and ligated to the specifically fluorescently labeled primer to achieve labeling. The tested allelic site ligation products are obtained in 2 ligation reactions, including the fluorescent labeling reaction and the allele-specific ligation reaction. The allele-specific ligation products of each locus are distinguished by labeling with different fluorescence, while the different loci are distinguished by 3` terminal probes and different length ligation sequences.

Statistical Analysis

Continuous variables were calculated using the Student’s t-test and displayed as means. Categorical variables and Weinberg equilibrium (HWE) were assessed using the X2-test and displayed as counts and percentages. We used logistic regression to assess the association of SNPs with LTBI/PTB, with sex and age as covariates. The relationship of haplotypes between each group was analyzed using online software (http://analysis.bio-x.cn). P values smaller than 0.05 were considered statistically significant. SPSS software (SPSS, Inc., Chicago, IL, USA) was used for all statistics analysis.

Results

Demographics of the Participants and results of Quality Control

A total of 209 PTB patients, 201 LTBI, and 204 HCS between the age group of 16 and 87 years were included in the present study. The male to female ratio of the three groups did not differ significantly; However, the mean age of the groups was observed to be significantly higher in PTB and HCS (P=0.027) (Table 1).

Table 1 Baseline Information of the Participants

A total of 3 Tag-SNPs were included in this study. According to the 5% samples repeat genotyping, the genotyping call rate of the three SNPs was 99.78%, and the accuracy rate was 100%. All 3 SNPs did not deviate from HWE in the control subjects. All SNPs information is shown in Table 2.

Table 2 Basic Information of All SNPs in Our Study

Polymorphisms of the KEAP1 in the Three Groups

Three SNPs in the KEAP1 gene was evaluated for associations with LTBI or PTB. We compared genotypes of LTBI with PTB and HCS, respectively. The distribution of genotypes among different groups is shown in Table 3. Significant association was found for KEAP1 rs1048290 between LTBI and HCS. Compared with the KEAP1 rs1048290 CC genotype, genotype GC had an 38% decreased risk for development LTBI (P = 0.043, OR = 0.62, 95% CI: 0.039–0.98). We also found that SNPs in KEAP1 were significantly related to PTB compared to LTBI. Compared with the rs11545829G allele, allele A had an 30% decreased risk for development PTB (P = 0.034, OR = 0.70, 95% CI: 0.51–0.97). We also found the rs11668429 polymorphism was related to PTB. Compared with TT, GT had a significantly increased risk of LTBI developing into PTB (P = 0.041, OR = 1.68, 95% CI: 1.02–2.77). Finally, we compared the differences in alleles and genotypes of KEAP1 between PTB and HCS. We found that rs11545829 AA (P = 0.017, OR = 0.84, 95% CI: 0.73–0.97) genotype and A (P = 0.031, OR = 0.71, 95% CI: 0.52–0.97) allele were associated with PTB.

Table 3 Association Between KEAP1 Genotypic/Allelic Frequencies and LTBI/PTB

LD Patterns and Haplotype Analysis

Pair-wise linkage disequilibrium (r2) analysis was conducted among the 3 SNPs of KEAP1. Haplotype blocks were built based on the r2 values (>0.7) (Figure 1). rs11545829 and rs11668429 were found in strong LD in the PTB and LTBI groups. No haplotypes were significantly associated with PTB or LTBI (Table 4).

Table 4 Haplotypes of the KEAP1 Genes and Their Distributions in the Three Groups

Figure 1 LD of KEAP1 gene polymorphisms in the both PTB vs LTBI (left) and PTBI vs HCS (right) populations LD r2 values (range from 0 to 1) for all pairs of SNPs are presented as percentages. Shading from white to black indicates LD measured as r2 (range from 0 to 1).

Discussion

Most previous genetic studies on TB have used associations between TB and genetic polymorphisms, with controls including LTBI and uninfected individuals. However, few scholars have studied the association between LTBI and gene polymorphism. We designed three groups, including PTB, LTBI, and HCS, to identify candidate genes for PTB and LTBI to find genetic markers at different stages of TB development. This study is the first to find that KEAP1 gene polymorphisms were associated with PTB and LTBI.

The past studies based on immunocompromised population, twin comparisons, candidate genes, and the genome-wide association study show that host genetic factors affect the progress of TB.16–19 Oxidative stress plays an important role in the susceptibility to TB according to literature.7 Some key factors related to TB have been identified, but further study of the fine-tuning of oxidative stress is necessary so as to further understand the role of antioxidant balance in TB treatment. The oxidative stress signaling pathway involved in KEAP1 was associated with a variety of diseases, such as insulin resistance,20 tumor,21 and neuronal apoptosis.22 Some studies have also shown that it is involved in respiratory diseases, like respiratory infections.23 Meanwhile, MTB is able to reduce the KEAP1 gene expression level.24 It has also been shown that rs1048290 and rs11545829 gene polymorphisms are associated with KEAP1 gene expression.25 Therefore, based on the association of KEAP1 with multiple respiratory diseases and the effect of MTB on KEAP1 expression, we hypothesized that KEAP1 polymorphims plays a role in PTB.

The association between KEAP1 polymorphisms and TB has not been studied, and this study was the first to explore the relationship between KEAP1 polymorphisms and PTB. Previously published studies have shown that KEAP1 gene polymorphisms were associated with the risk of various respiratory diseases. For example, KEAP1 rs11085735 and rs1048287 were associated with chronic obstructive pulmonary disease (COPD) and asthma.26 They also found that KEAP1 rs9676881 was closely related to COPD. Ungvári et al demonstrated that KEAP1 polymorphisms are associated with asthma.26 However, the association of KEAP1 gene polymorphisms with TB has not been clearly established.

In a cohort of study in a Russian population, the rs1048290 polymorphism was reported to be associated with COPD, with a particularly significant association in the additive model. In their study, the C allele frequency was lower in the COPD group.12 Another study showed that the rs1048290 gene polymorphism was a susceptibility factor for rectal cancer. Among them, the CC genotype and the C allele increase the risk of colorectal cancer.25 Meanwhile, rs1048290 was associated with KEAP1 mRNA expression level, and KEAP1 expression level was relatively decreased in CC genotype. In our study, rs1048290 GC genotype compared with CC was a protective factor for LTBI compare with HCS. rs11545829 was in the fifth exon in KEAP1. There were few studies on the rs11545829 gene polymorphisms, and a study in Chinese population showed that the rs11545829 AA genotype has a protective effect on type 2 diabetes mellitus.27 Another study showed that the rs1048290-rs11545829 GT haplotype was a low risk factor for colorectal cancer.25 In the current study, we found that the rs11545829 A allele was a protective factor for PTB. There is only one study on the rs11668429 polymorphism so far.26 Ungvári et al explored the association of rs11668429 with childhood asthma, but they found no association between them. In the present study, the rs11668429 polymorphism was associated with PTB, and GT genotype was a risk factor for PTB. More researches are needed to validate our results.

There are still some deficiencies in this study. First, the included SNPs loci lack functional validation. Therefore, the relationship between the KEAP1 gene and PTB is still unknown. Second, the results of this study have not been validated in an independent population and therefore may lead to type I errors. Finally, the association of clinical information and KEAP1 polymorphisms was not analyzed due to the incomplete collection of clinical data.

Conclusion

In conclusion, our study suggested that KEAP1 polymorphisms were significantly related to susceptibility to PTB and LTBI subjects. This study may reflect the role of gene polymorphisms in different stages of TB progression, potential to help detect and treat TB early, but more studies are needed to verify the results.

Data Sharing Statement

Please contact corresponding author for original data.

Funding

The National Natural Science Foundation of China [grant number 81870015] supported the research. The Medical Science and Technology Project of Sichuan Provincial Health Commission [grant number 21PJ025] supported the research. Post-Doctor Research Project, West China Hospital, Sichuan University [grant number 2021HXBH075] supported the research.

Disclosure

The authors report no conflicts of interest in this work.

References

1. World Health Organization. Global Tuberculosis Report 2021. Geneva: World Health Organization; 2021.

2. Vynnycky E, Fine PE. Lifetime risks, incubation period, and serial interval of tuberculosis. Am J Epidemiol. 2000;152(3):247–263. doi:10.1093/aje/152.3.247

3. Hill AV. The genomics and genetics of human infectious disease susceptibility. Annu Rev Genomics Hum Genet. 2001;2(1):373–400. doi:10.1146/annurev.genom.2.1.373

4. Casanova J-L, Abel L. Genetic dissection of immunity to mycobacteria: the human model. Annu Rev Immunol. 2002;20(1):581–620. doi:10.1146/annurev.immunol.20.081501.125851

5. Biswas SK, Mittal M, Sinha E, et al. Exploring the role of CC motif chemokine ligand-2 single nucleotide polymorphism in pulmonary tuberculosis: a genetic association study from North India. J Immunol Res. 2020;2020. doi:10.1155/2020/1019639

6. Ji G, Zhang M, Liu Q, et al. Functional polymorphism in the NFE2L2 gene associated with tuberculosis susceptibility. Front Immunol. 2021;12:1714. doi:10.3389/fimmu.2021.660384

7. Shastri MD, Shukla SD, Chong WC, et al. Role of oxidative stress in the pathology and management of human tuberculosis. Oxid Med Cell Longev. 2018;2018:1–10. doi:10.1155/2018/7695364

8. Niture SK, Khatri R, Jaiswal AK. Regulation of Nrf2—an update. Free Radic Biol Med. 2014;66:36–44. doi:10.1016/j.freeradbiomed.2013.02.008

9. Taguchi K, Motohashi H, Yamamoto M. Molecular mechanisms of the Keap1–Nrf2 pathway in stress response and cancer evolution. Genes Cells. 2011;16(2):123–140. doi:10.1111/j.1365-2443.2010.01473.x

10. Sun Q, Shen X, Ma J, Lou H, Zhang Q. Activation of Nrf2 signaling by oltipraz inhibits death of human macrophages with Mycobacterium tuberculosis infection. Biochem Biophys Res Commun. 2020;531(3):312–319. doi:10.1016/j.bbrc.2020.07.026

11. Hartikainen JM, Tengström M, Winqvist R, et al. KEAP1 genetic polymorphisms associate with breast cancer risk and survival outcomesKEAP1 genetic polymorphism predicts outcome in breast cancer. Clin Cancer Res. 2015;21(7):1591–1601. doi:10.1158/1078-0432.CCR-14-1887

12. Korytina GF, Akhmadishina LZ, Aznabaeva YG, et al. Associations of the NRF2/KEAP1 pathway and antioxidant defense gene polymorphisms with chronic obstructive pulmonary disease. Gene. 2019;692:102–112. doi:10.1016/j.gene.2018.12.061

13. Wong TF, Yoshinaga K, Monma Y, et al. Association of keap1 and nrf2 genetic mutations and polymorphisms with endometrioid endometrial adenocarcinoma survival. Int J Gynecol Canc. 2011;21:8. doi:10.1097/IGC.0b013e31822d0eb2

14. Wu S, Huang W, Wang D, et al. Evaluation of TLR 2, TLR 4, and TOLLIP polymorphisms for their role in tuberculosis susceptibility. Apmis. 2018;126(6):501–508. doi:10.1111/apm.12855

15. Wu S, Wang Y, Zhang M, Shrestha SS, Wang M, He J-Q. Genetic polymorphisms of IL1B, IL6, and TNFα in a Chinese Han population with pulmonary tuberculosis. In: BioMed Research International; 2018:2018.

16. Abel L, El-Baghdadi J, Bousfiha AA, Casanova J-L, Schurr E. Human genetics of tuberculosis: a long and winding road. Philos Trans R Soc B. 2014;369(1645):20130428. doi:10.1098/rstb.2013.0428

17. Sester M, Van Leth F, Bruchfeld J, et al. Risk assessment of tuberculosis in immunocompromised patients. A TBNET study. Am J Respir Crit Care Med. 2014;190(10):1168–1176. doi:10.1164/rccm.201405-0967OC

18. Thye T, Owusu-Dabo E, Vannberg FO, et al. Common variants at 11p13 are associated with susceptibility to tuberculosis. Nat Genet. 2012;44(3):257–259. doi:10.1038/ng.1080

19. Chimusa ER, Zaitlen N, Daya M, et al. Genome-wide association study of ancestry-specific TB risk in the South African Coloured population. Hum Mol Genet. 2014;23(3):796–809. doi:10.1093/hmg/ddt462

20. Ding X, Jian T, Wu Y, et al. Ellagic acid ameliorates oxidative stress and insulin resistance in high glucose-treated HepG2 cells via miR-223/keap1-Nrf2 pathway. Biomed Pharmacother. 2019;110:85–94. doi:10.1016/j.biopha.2018.11.018

21. Sajadimajd S, Khazaei M. Oxidative stress and cancer: the role of Nrf2. Curr Cancer Drug Targets. 2018;18(6):538–557. doi:10.2174/1568009617666171002144228

22. Li D, Bai X, Jiang Y, Cheng Y. Butyrate alleviates PTZ-induced mitochondrial dysfunction, oxidative stress and neuron apoptosis in mice via Keap1/Nrf2/HO-1 pathway. Brain Res Bull. 2021;168:25–35. doi:10.1016/j.brainresbull.2020.12.009

23. Liu Q, Gao Y, Ci X. Role of Nrf2 and its activators in respiratory diseases. Oxid Med Cell Longev. 2019;2019. doi:10.1155/2019/7090534

24. Sonavane A. Mycobacterium tuberculosis Rv3034c regulates mTORC1 and PPAR-γ dependant pexophagy mechanism to control redox levels in macrophages. Cell Microbiol. 2020;22:e13214.

25. Gao L, Yuan F, Che G, et al. Epigenetic modifications but not genetic polymorphisms regulate KEAP1 expression in colorectal cancer. J Cell Biochem. 2019;120(8):12311–12320. doi:10.1002/jcb.28495

26. Ungvári I, Hadadi É, Virág V, et al. Relationship between air pollution, NFE2L2 gene polymorphisms and childhood asthma in a Hungarian population. J Community Genet. 2012;3(1):25–33. doi:10.1007/s12687-011-0075-8

27. Fan C, Zhan Z, Zhang X, et al. Research for type 2 diabetes mellitus in endemic arsenism areas in central China: role of low level of arsenic exposure and KEAP1 rs11545829 polymorphism. Arch Toxicol. 2022;96:1–11. doi:10.1007/s00204-021-03184-z

Creative Commons License © 2022 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.