Association between genetic variants in CD1A and CD1D genes and pulmonary tuberculosis in an Iranian population
- Authors:
- Published online on: March 15, 2019 https://doi.org/10.3892/br.2019.1201
- Pages: 259-265
Abstract
Introduction
Mycobacterium tuberculosis (MTB) is respiratory tract infectious disease but can affect nearly all of the body, including the brain, the kidneys and the bones (1). Tuberculosis is still a major health problem and a leading cause of mortalities from infectious disease throughout the world (2). Based on a World Health Organization report 10.4 million new TB cases and 1.3 million fatalities occurred in 2016 globally (3). The prevalence of TB in Sistan and Baluchistan province is higher than other provinces of Iran. Whereas one-third of the world's population is infected with MTB, 5-10% of infected cases go on to develop active TB, which indicates that in addition to the environment host genetics factors may influence the risk of the disease (2,4,5).
Cluster of differentiation 1 (CD1) genes, mapped to human chromosome 1 (1q23.1), encoding five isoforms of the CD1 molecule, namely CD1a, CD1b, CD1c, CD1d and CD1e (6-9). CD1 is a family of glycoproteins which are expressed on antigen presenting cells (APCs) (10). They are belongs to major histocompatibility complex class Ι (MHC) and participate in presentation of lipid antigens to T cells (11) but unlike MHCs their expression is nearly identical in all humans (11). The CD1 family can be divided into three groups. Group 1 consists of CD1a, CD1b and CD1c that present lipid antigens to T-helper cells. CD1D belongs to group 2 and presents lipid antigens to natural killer T cells and group 3 that composed of CD1E (12).
It is clear that CD1 molecules participate in cell-mediated immune responses against MTB (13). Indeed, several kinds of mycobacterial lipids function as antigens in this system (11). CD1-restricted T cells identify MTB and react to mycobacterium cell wall lipid antigens (14). CD1 molecules, unlike classical MHC molecules, exhibit limited polymorphism (10,15) but, these few polymorphisms are demonstrated to be associated with susceptibility or resistance to certain diseases (12,14,16,17). Single nucleotide polymorphism (SNPs), the most common type of genetic variants can affect the expression and function of the gene (18,19). There is little data regarding the impact of CD1 gene polymorphisms and tuberculosis susceptibility (14). Seshadri et al (14) investigated the association between rs366316, rs2269714, rs411089 and rs389293 polymorphisms of CD1A gene and risk of tuberculosis. They reported that rs411089 and rs366316 variants were significantly correlated with the development of tuberculosis. To the best of our knowledge, no study investigated the impact of CDID polymorphisms on tuberculosis susceptibility. The present study aimed to examine the association between CD1A (rs411089 and rs366316) and CD1D (rs973742, rs859009 and rs859010) polymorphisms and the risk of PTB in a sample of the southeast Iranian population.
Materials and methods
Patients
A total of 352 subjects including 172 clinically diagnosed patients with pulmonary tuberculosis (PTB) between 12 and 86 years of age (69 males and 103 females), and 180 unrelated healthy individuals ranging between 20 and 85 years of age (74 males and 106 females) were recruited in the study from May 2014 to March 2018. The cases were selected from PTB patients admitted to the University Affiliated Hospital (Bouali Hospital, Zahedan, Iran; referral center for TB). The diagnosis of PTB depended on clinical findings, radiological evidence, positive smear for MTB, culture and response to antituberculosis chemotherapy as described previously (20-23). The control individuals did not have a history of TB and inflammatory disease or other of chronic infectious disease; they were of the same ethnicity as patients and living in the same area as the patients with PTB (Southeast Iran). The local Ethics Committee of the Zahedan University of Medical Sciences (Zahedan, Iran) approved the project and written informed consent was taken from all participants. Blood samples (5 ml) were collected in Na-EDTA tubes from patients and healthy controls and stored at -20 until DNA extraction. Genomic DNA was extracted from whole blood by salting-out method as described previously (24).
Genotyping
Genotyping of CD1A and CD1D was achieved by polymerase chain reaction amplification-restriction fragment length polymorphism assays. The primers were synthesized by Metabion International AG (Erding, Germany) and the restriction enzyme were purchased from Fermentas; Thermo Fisher Scientific, Inc. (Waltham, MA, USA). Primer sequences restriction enzymes and length of the fragments are summarized in Table I. Amplification was achieved in a final volume of 20 µl containing 1 µl genomic DNA (~100 ng/ml), 1 µl each primer (10 µM) and 10 µl of 2X Prime Taq Premix (Genet Bio, Daejeon, Korea), and 7 µl ddH2O. The PCR conditions were set as follows: 95˚C for 5 min, 30 cycles of 95˚C for 30 sec, 62˚C for CD1A rs411089, CD1D rs859009 and CD1D rs859010, 60˚C for CD1D rs973742, 65˚C for CD1A rs366316 and 72˚C for 30 sec and a final extension step of 72˚C for 5 min. Then, 10 µl of PCR product was digested with suitable restriction enzyme (Table II) and resolved on 2.5% agarose gel containing 0.5 µg/ml ethidium bromide and visualized under UV light. A total of ~10% of the random samples were regenotyped and the finding were confirmed as 100% concordant.
Table IPrimer sequences of polymerase chain reaction-restriction fragment length polymorphism for detection of CD1A and CD1D polymorphisms. |
Table IIThe genotypes and allele frequency distribution of CD1A polymorphism in pulmonary tuberculosis patients and control groups. |
Statistical analysis
Analysis of the data was achieved by the statistical package SPSS 20 software (IBM Corp., Armonk, NY, USA). Association among variants and PTB risk was estimated by calculation the odds ratio (OR) and 95% confidence intervals (CI) from logistic regression analysis. The Bonferroni correction test was performed by multiplying P-values of logistic regression analysis by the number of SNPs examined for each gene.
Haplotype analysis was executed using SNPStats: A web tool for the analysis of association studies (25). P<0.05 was considered to indicate a statistically significant difference.
The observed genotype frequencies in the controls were examined for Hardy-Weinberg equilibrium (HWE) using the χ2-test. The Bonferroni correction test was performed by multiplying P-values by the number of SNPs examined for each gene.
Results
Patients
The study participants include 172 patients with PTB (69 males and 103 females; mean age, 50.01±20.57 years) and 180 healthy individuals (74 males and 106 females; mean age, 49.92±15.10 years). There was no statistically significant difference between the groups regarding sex and age (P=0.486 and 0.572, respectively).
Association between the CD1A variants and PTB risk
The genotype and allelic frequencies of CD1A rs411089 and CD1A rs366316 gene polymorphisms in cases and controls are presented in Table II. The results demonstrated that the CC genotype of CD1A rs411089 polymorphism in recessive inheritance model increased the risk of PTB in comparison with TT+TC genotypes (OR=2.71, 95% CI=1.36-5.40, P=0.004). As presented in Table II there was no significant difference in genotype and allelic distribution of CD1A rs366316 gene polymorphisms between PTB patients and controls (P>0.05).
Association between the CD1D variants and PTB risk
Table III summarized the genotype and allele frequencies of CD1D rs973742, CD1D rs859009 and CD1D rs859010 polymorphisms in PTB patients and the control groups. The GC genotype as well as the GC+CC genotype of CD1D rs859009 polymorphism decreased the risk of PTB (OR=0.49, 95% CI=0.29-0.85, P=0.01 and OR=0.53, 95% CI=0.31-0.89, P=0.01, respectively). The allele frequency of CD1D rs859009 Polymorphism was not significantly different between the groups (P=0.22). The genotype and allele frequencies of CD1D rs973742 and CD1D rs859010 were not significantly different between two groups (χ2=2.206, P=0.332 and χ2=2.252, P=0.324, respectively). The results demonstrated that CD1D rs973742 and CD1D rs859010 variants were not associated with risk of/protection from PTB in codominant, dominant and recessive tested inheritance models (Table II). The results demonstrated that the rs366316, rs859010 polymorphisms but not rs411089, rs973742 and rs859009 polymorphisms were in present in HWE, Tables II and III).
Table IIIThe genotypes and allele distribution of CD1D polymorphisms in pulmonary tuberculosis patients and control groups. |
Haplotypes analysis
Haplotypes analysis was achieved (Table IV) and the results demonstrated that haplotypes TGACA, TAACA, and TAAGA significantly decreased the risk of PTB in comparison with CAGGA (rs411089C/rs366316A/rs973742G/rs859009G/rs859010A).
Table IVHaplotype association of CD1A (rs411089 and rs366316) and CD1D (rs973742, rs859009 and rs859010) variants with pulmonary tuberculosis risk. |
Discussion
CD1 is a family of antigen presenting molecules that belongs to the MHC class I molecules but, unlike those that exist on all nucleated cells, the CD1 groups are expressed only on APCs, cortical thymocytes and dendritic cells (DC) (10,26). The human CD1 proteins present lipid antigens of MTB to T cells (27). Mounting evidence revealed a link between CD1 polymorphisms and certain autoimmune and inflammatory diseases including inflammatory neuropathies (12), tuberculosis (14), cancer (28), multiple sclerosis (16,29), and Guillain-Barre syndromes (30). In the current study the possible association between CD1 and the risk pf PTB in a sample of a southeast Iranian population was investigated. The results of the present study indicated that CD1A rs411089 as well as CD1D rs859009 variant was significantly associated with the risk of PTB. The CC genotype of CD1A rs411089 polymorphism in recessive model increased the risk of PTB, while the GC genotype as well as GC+CC genotype of CD1D rs859009 polymorphism significantly decreased the risk of PTB.
To the best of the authors' knowledge, only one study has been performed on the impact of CD1 gene variants on the susceptibility to TB (14). Seshadri et al (14) in a case-population study genotyped rs366316, rs2269714, rs411089 and rs389293 in a discovery cohort of 352 cases and 382 controls. They identified that the rs366316 and rs411089 polymorphisms were significantly associated with the development of TB. In a next step they genotyped rs366316 and rs411089 in a validation cohort of 339 cases and 376 controls and another time discovered a significant association between rs411089 and TB, however they did not find any association between rs366316 and TB (14). In agreement with the results of the present study they demonstrated that the minor homozygous genotype of rs411089 was associated with an increased risk of TB in a recessive model (14). It seems that the minor homozygous genotype of rs411089 accompanied with a functional CD1A-deficiency and this polymorphism cause a low CD1A expression level. Sieling et al (31) demonstrated an increased expression of CD1 proteins (CD1A, CD1B and CD1C) at the site of tuberculoid leprosy lesions. Buettner et al (32) detected DCs in bronchoalveolar lavage specimens of tuberculosis patients with strongly upregulated of CD1A, CD83 and CCR7. Roura-Mir et al (33) have identified that MTB or its cell wall lipids caused increased expression level of CD1A, CD1B and CD1C but not CD1D. CD1A-deficient cells failed to present mycobacterium peptides to T cells. Behar et al (34) demonstrated that CD1D-deficient mice were not significantly different in their susceptibility to infection compared with control mice. This information demonstrates the use of a population genetics study to understanding the role of CD1 antigen-presenting molecules in human disease.
It has been reported that the rs366316 variant which is located in the in the promoter region is strongly associated with low CD1a surface expression (18). The rs411089, an intronic variant, is strongly linked to rs858998 and may regulate transcription independently of rs366316(14). rs973742 is located 4-kb downstream from the CD1D gene and the functional role of this variant is unknown (35). The impact of intronic rs859009 and rs859010 polymorphisms on CD1D expression has not been recognized.
There is no clear explanation for the departure from HWE in the study population regarding rs411089, rs973742 and rs859009 polymorphisms. It may be owing to a small sample size, genetic drift or consanguineous marriages, which is common in this region of the country.
A limitation of this study is that there is no data regarding smoking, diabetes mellitus, socioeconomic status and other risk factor for PTB.
In conclusion, the results of the present study revealed that CD1A (rs411089) and CD1D (rs859009) polymorphisms protect against PTB in a sample of the Iranian population. More well-designed studies in diverse ethnicities are necessary to verify the present results.
Acknowledgements
We would like to thank all individuals who participated in the study.
Funding
The present study was supported by the Deputy for research, Zahedan University of Medical Sciences, Zahedan, Iran (grant no. 7033).
Availability of data and materials
All data generated or analyzed during the present study are included in this published article.
Authors' contributions
MT, MH conceived and designed the study; HD, FB and GB performed the experiments and analyzed the data; MN performed data collection. All authors contributed to the writing of the manuscript and reviewed and approved the final manuscript.
Ethics approval and consent to participate
The local Ethics Committee of the Zahedan University of Medical Sciences approved the project and written informed consent was taken from all participants.
Patient consent for publication
Not applicable.
Conflict of interest
The authors declare no conflict of interest. The abstract of the present study was previously presented at a conference (https://www.nature.com/articles/s41431-018-0248-6).
References
Smith I: Mycobacterium tuberculosis pathogenesis and molecular determinants of virulence. Clin Microbiol Rev. 16:463–496. 2003.PubMed/NCBI View Article : Google Scholar | |
Harishankar M, Selvaraj P and Bethunaickan R: Influence of Genetic Polymorphism Towards Pulmonary Tuberculosis Susceptibility. Front Med (Lausanne). 5(213)2018.PubMed/NCBI View Article : Google Scholar | |
WHO: Global Tuberculosis Report 2017. World Health Organization, 2017. | |
Wei Z, Wenhao S, Yuanyuan M, Yang L, Daming Z, Jiangchun X and Jijun J: A single nucleotide polymorphism in the interferon-γ gene (IFNG +874 T/A) is associated with susceptibility to tuberculosis. Oncotarget. 8:50415–50429. 2017.PubMed/NCBI View Article : Google Scholar | |
Naderi M, Hashemi M and Ansari H: Macrophage migration inhibitory factor -173G>C polymorphism and risk of tuberculosis: A meta-analysis. EXCLI J. 16:313–320. 2017.PubMed/NCBI View Article : Google Scholar | |
Barral DC and Brenner MB: CD1 antigen presentation: How it works. Nat Rev Immunol. 7:929–941. 2007.PubMed/NCBI View Article : Google Scholar | |
Brigl M and Brenner MB: CD1: Antigen presentation and T cell function. Annu Rev Immunol. 22:817–890. 2004.PubMed/NCBI View Article : Google Scholar | |
Aureli A, Oumhani K, Del Beato T, Di Rocco M, Tessitore A, El Aouad R and Piancatelli D: Increased CD1D polymorphism: Identification of two novel alleles, CD1D*03 and *04, in individuals from Morocco. Int J Immunogenet. 42:287–291. 2015.PubMed/NCBI View Article : Google Scholar | |
Chancellor A, Gadola SD and Mansour S: The versatility of the CD1 lipid antigen presentation pathway. Immunology. 154:196–203. 2018.PubMed/NCBI View Article : Google Scholar | |
Siddiqui S, Visvabharathy L and Wang CR: Role of group 1 CD1-restricted T cells in infectious disease. Front Immunol. 6(337)2015.PubMed/NCBI View Article : Google Scholar | |
Van Rhijn I and Moody DB: CD1 and mycobacterial lipids activate human T cells. Immunol Rev. 264:138–153. 2015.PubMed/NCBI View Article : Google Scholar | |
Uncini A, Notturno F, Pace M and Caporale CM: Polymorphism of CD1 and SH2D2A genes in inflammatory neuropathies. J Peripher Nerv Syst. 16 (Suppl 1):48–51. 2011.PubMed/NCBI View Article : Google Scholar | |
Huang S: Targeting Innate-Like T Cells in Tuberculosis. Front Immunol. 7(594)2016.PubMed/NCBI View Article : Google Scholar | |
Seshadri C, Thuong NT, Yen NT, Bang ND, Chau TT, Thwaites GE, Dunstan SJ and Hawn TR: A polymorphism in human CD1A is associated with susceptibility to tuberculosis. Genes Immun. 15:195–198. 2014.PubMed/NCBI View Article : Google Scholar | |
DeWitt WS, Yu KKQ, Wilburn DB, Sherwood A, Vignali M, Day CL, Scriba TJ, Robins HS, Swanson WJ, Emerson RO, et al: A diverse lipid antigen-specific TCR repertoire is clonally expanded during active tuberculosis. J Immunol. 201:888–896. 2018.PubMed/NCBI View Article : Google Scholar | |
Caporale CM, Notturno F, Pace M, Aureli A, Di Tommaso V, De Luca G, Farina D, Giovannini A and Uncini A: CD1A and CD1E gene polymorphisms are associated with susceptibility to multiple sclerosis. Int J Immunopathol Pharmacol. 24:175–183. 2011.PubMed/NCBI View Article : Google Scholar | |
Kuijf ML, Geleijns K, Ennaji N, van Rijs W, van Doorn PA and Jacobs BC: Susceptibility to Guillain-Barré syndrome is not associated with CD1A and CD1E gene polymorphisms. J Neuroimmunol. 205:110–112. 2008.PubMed/NCBI View Article : Google Scholar | |
Seshadri C, Shenoy M, Wells RD, Hensley-McBain T, Andersen-Nissen E, McElrath MJ, Cheng TY, Moody DB and Hawn TR: Human CD1a deficiency is common and genetically regulated. J Immunol. 191:1586–1593. 2013.PubMed/NCBI View Article : Google Scholar | |
Lundtoft C, Awuah AA, Güler A, Harling K, Schaal H, Mayatepek E, Phillips RO, Nausch N, Owusu-Dabo E and Jacobsen M: An IL7RA exon 5 polymorphism is associated with impaired IL-7Rα splicing and protection against tuberculosis in Ghana. Genes Immun: Oct 31. 2018.(Epub ahead of print). PubMed/NCBI View Article : Google Scholar | |
Hashemi M, Eskandari-Nasab E, Moazeni-Roodi A, Naderi M, Sharifi-Mood B and Taheri M: Association of CTSZ rs34069356 and MC3R rs6127698 gene polymorphisms with pulmonary tuberculosis. Int J Tuberc Lung Dis. 17:1224–1228. 2013.PubMed/NCBI View Article : Google Scholar | |
Naderi M, Hashemi M, Safdari A, Bahari G and Taheri M: Association of genetic polymorphisms of CISH with the risk of pulmonary tuberculosis in Zahedan, Southeast Iran. Braz J Infect Dis. 20:379–383. 2016.PubMed/NCBI View Article : Google Scholar | |
Naderi M, Hashemi M and Amininia S: Association of TAP1 and TAP2 gene polymorphisms with susceptibility to pulmonary tuberculosis. Iran J Allergy Asthma Immunol. 15:62–68. 2016.PubMed/NCBI | |
Naderi M, Hashemi M, Pourmontaseri Z, Eskandari-Nasab E, Bahari G and Taheri M: TIRAP rs8177374 gene polymorphism increased the risk of pulmonary tuberculosis in Zahedan, southeast Iran. Asian Pac J Trop Med. 7:451–455. 2014.PubMed/NCBI View Article : Google Scholar | |
Hashemi M, Hanafi Bojd H, Eskandari Nasab E, Bahari A, Hashemzehi NA, Shafieipour S, Narouie B, Taheri M and Ghavami S: Association of adiponectin rs1501299 and rs266729 gene polymorphisms with nonalcoholic fatty liver disease. Hepat Mon. 13(e9527)2013.PubMed/NCBI View Article : Google Scholar | |
Solé X, Guinó E, Valls J, Iniesta R and Moreno V: SNPStats: A web tool for the analysis of association studies. Bioinformatics. 22:1928–1929. 2006.PubMed/NCBI View Article : Google Scholar | |
Ulrichs T, Moody DB, Grant E, Kaufmann SH and Porcelli SA: T-cell responses to CD1-presented lipid antigens in humans with Mycobacterium tuberculosis infection. Infect Immun. 71:3076–3087. 2003.PubMed/NCBI View Article : Google Scholar | |
Seshadri C, Thuong NT, Mai NT, Bang ND, Chau TT, Lewinsohn DM, Thwaites GE, Dunstan SJ and Hawn TR: A polymorphism in human MR1 is associated with mRNA expression and susceptibility to tuberculosis. Genes Immun. 18:8–14. 2017.PubMed/NCBI View Article : Google Scholar | |
Golmoghaddam H, Pezeshki AM, Ghaderi A and Doroudchi M: CD1a and CD1d genes polymorphisms in breast, colorectal and lung cancers. Pathol Oncol Res. 17:669–675. 2011.PubMed/NCBI View Article : Google Scholar | |
Jamshidian A, Nikseresht AR, Vessal M and Kamali-Sarvestani E: Association of CD1A +622 T/C, +737 G/C and CD1E +6129 A/G genes polymorphisms with multiple sclerosis. Immunol Invest. 39:874–889. 2010.PubMed/NCBI View Article : Google Scholar | |
Liu H, Xing Y, Guo Y, Liu P, Zhang H, Xue B, Shou J, Qian J, Peng J, Wang R, et al: Polymorphisms in exon 2 of CD1 genes are associated with susceptibility to Guillain-Barré syndrome. J Neurol Sci. 369:39–42. 2016.PubMed/NCBI View Article : Google Scholar | |
Sieling PA, Jullien D, Dahlem M, Tedder TF, Rea TH, Modlin RL and Porcelli SA: CD1 expression by dendritic cells in human leprosy lesions: Correlation with effective host immunity. J Immunol. 162:1851–1858. 1999.PubMed/NCBI | |
Buettner M, Meinken C, Bastian M, Bhat R, Stössel E, Faller G, Cianciolo G, Ficker J, Wagner M, Röllinghoff M, et al: Inverse correlation of maturity and antibacterial activity in human dendritic cells. J Immunol. 174:4203–4209. 2005.PubMed/NCBI View Article : Google Scholar | |
Roura-Mir C, Wang L, Cheng TY, Matsunaga I, Dascher CC, Peng SL, Fenton MJ, Kirschning C and Moody DB: Mycobacterium tuberculosis regulates CD1 antigen presentation pathways through TLR-2. J Immunol. 175:1758–1766. 2005.PubMed/NCBI View Article : Google Scholar | |
Behar SM, Dascher CC, Grusby MJ, Wang CR and Brenner MB: Susceptibility of mice deficient in CD1D or TAP1 to infection with Mycobacterium tuberculosis. J Exp Med. 189:1973–1980. 1999.PubMed/NCBI View Article : Google Scholar | |
Barbier M, Sabbagh A, Kasper E, Asheuer M, Ahouansou O, Pribill I, Forss-Petter S, Vidaud M, Berger J and Aubourg P: CD1 gene polymorphisms and phenotypic variability in X-linked adrenoleukodystrophy. PLoS One. 7(e29872)2012.PubMed/NCBI View Article : Google Scholar |