Association of Chemokine Genes CXCL9 and CXCL10 Polymorphisms with Tuberculosis in Pakistani Population

Saleha Tabassum1, Muneeza Zafar1,2, Munazza Raza Mirza1*, Muhammad Iqbal Choudhary1,3, Ikram ul Haq4, Imran Riaz Malik2, Ambreen Ansari1 and Mehtab Alam5

1Dr. Panjwani Center for Molecular Medicine and Drug Research, International Center for Chemical and Biological Sciences, University of Karachi, 75270, Pakistan.

2University of Sargodha, Sargodha, Pakistan.

3H.E.J. Research Institute of Chemistry, International Center for Chemical and Biological Sciences, University of Karachi, 75270, Pakistan.

4Institute of Biotechnology and Genetic Engineering, University of Sindh, Jamshoro, Pakistan.

5Department of Biochemistry, Dow University of Medical and Health Sciences, Ojha Campus. Karachi, Pakistan.

Saleha Tabassum and Muneeza Zafar equally contributed and are co-first authors.

ABSTRACT

The CXCR3 ligand family, including CXCL9, CXCL10 and CXCL11, play important role in T-cell, natural killer (NK) cell, NKT cell and dendritic cell trafficking. These chemokines were suggested for their use as alternative markers for TB diagnosis, and their possible involvement in disease pathophysiology. In the present study, we investigated the genetic association and frequency of CXCL9 and CXCL10 genes polymorphisms in TB patients of Pakistani population. A total of 260 study participants, including 122 with tuberculosis disease (TB) and 138 healthy control (HC) were included. Out of these 122 TB disease subjects (67%) were diagnosed with PTB (n = 82) and 33% with EPTB (n = 40). Tri and Tetra- ARMS-PCR techniques were used to genotype rs2276886 of CXCL9 and rs5606198 of CXCL10 for all the samples. The results showed significant association (OR, 2.1 (0.99 to 4.30), p = 0.04) of A allele of rs2276886 with EPTB group as compared to the control. In genetic model analysis, we found strong association (OR, 0.4 (0.181 to 0.94), p = 0.03) of GG genotype with male with EPTB as compared to female in gender stratified dominant model. Furthermore, in context of different ethnic groups of Pakistani subpopulation, we found Punjabi speaking’s to be significantly at higher risk (OR, 2.16 (95%CI, 1.0910 to 4.2901), p = 0.02) of developing TB as compared to other ethnic categories. Strong association of CXCL9 polymorphism (rs2276886 G>A) with TB disease in Pakistani population was concluded especially, Punjabi speaking participants were found to have higher risk of developing TB.

Article Information

Received 10 August 2021

Revised 26 September 2021

Accepted 01 October 2021

Available online 08 February 2022

(early access)

Published 29 August 2022

Authors’ Contribution

ST, MZ, MIC, MRM and AA conceptualized the work. ST and MZ performed experiments, analyzed data, drafted the original manuscript. MRM performed supervision and funding acquisition. MIC provided lab facilities and reviewed the manuscript. MRM, MA, IH, MIR, ST and MZ reviewed and edited the manuscript.

Key words

Tuberculosis, Cytokines, Polymorphism, Genotype, Inflammation

DOI: https://dx.doi.org/10.17582/journal.pjz/20210810130817

* Corresponding author: [email protected]

0030-9923/2022/0006-2871 $ 9.00/0

Copyright 2022 by the authors. Licensee Zoological Society of Pakistan.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Abbreviations

CXCL10, CXC ligand 10; MTB, Mycobacterium tuberculosis; TB, Tuberculosis; SNP, single nucleotide gene polymorphism; AFB, acid fast bacillus, PT (PTB), pulmonary tuberculosis disease; EPT (EPTB), extra pulmonary tuberculosis disease

INTRODUCTION

Tuberculosis (TB) is a serious, progressive bacterial infection caused by Mycobacterium tuberculosis (MTB) (Loddenkemper et al., 2015). The disease is endemic in half of the world with a significant burden reported from the developing countries (Doherty, 2005). About two billion or nearly one-third of the world’s population is infected with MTB. The global prevalence of the disease was estimated to be around 16-20 million, with 8 million new TB cases emerge every year (WHO, 2005; Corbett et al., 2003).Pakistan ranks sixth among the twenty-two countries with the highest burden (55%) of TB in the world (Metzger et al., 2010). Frequency of MTB varies among different ethnic groups, races, and genetic predisposition factor is also involved in the susceptibility of the disease (Nonghanphithak et al., 2015).

Lungs have been identified as the primary site of the MTB infection. Nonetheless, the M. bacillus is able to disseminate into other organs following septicemia (Loddenkemper et al., 2015). The histopathological examination of MTB patients suggests localized, proliferative and productive inflammatory responses over the course of the infection. The exudative immunity is often followed by the cell-mediated immunity, proceeds by lymphocytes including T cells and B cells. The generation of the cytotoxic CD4+ T cell response serves as the critical factor in defining the immunity (Winslow et al., 2008; Flynn and Chan, 2001). The chemokines are the primary cytokines responsible for immunological response, reported to play crucial role in the early inflammation of M. tuberculosis infection (Cooper et al., 2002; Flynn and Chan, 2001).

The primary cytokines responsible for the immunological response include tumor necrosis-alpha (TNF-α), interleukins, and chemokines (Cooper et al., 2002; Flynn and Chan, 2001). Chemokines are small chemotactic peptides, secreted by macrophages and are involved in the attraction of a wide range of immune cells, including lymphocytes, monocytes, dendritic cells, granulocytes and neutrophils to the site of infection. CXC ligand 10 (CXCL10) and CXC ligand 9 (CXCL9) are important CXC chemokines, which have been extensively studied due to their role in different pathophysiological conditions such as M. tuberculosis (MTB) infection, Hhimoto’s thyroiditis (HT), chronic rheumatic heart disease, thyroiditis and hypothyroidism, etc. (Fae et al., 2013; Antonelli et al., 2011; Akahane et al., 2016). CXCL10 is a 10 kDa protein which is produced in response to interferon-gamma (IFN-γ) (Seiler et al., 2003). It regulates immune responses including inflammation, angiogenesis, and chemotaxis (Esche et al., 2005; Turner et al., 2014). Several studies showed the association of the potential CXCL10 promoter region SNPs including at 135 bp (135G→A) upstream of the transcription start site with MTB disease (Tang et al., 2009; Sheikh et al., 2015). Bhattacharyya et al. (2018) have reported significant overexpression of CXCL10 in MTB patients as compared to the control (p = 0.0002). Furthermore, the SNP in CXCL10 rs4508917 (-1447 A > G) was associated with reduced expression of CXCL10, suggesting the essential role of CXCL10 in the active phase of TB infection (Bhattacharyya et al., 2018). A study conducted on the Polish population showed the higher levels of CXCL10, which is suggested to be used as a biomarker to screen the patients with latent MTB infection (Strzelak et al., 2012).

The CXCL9, also known as MIG (monokine-induced by interferon-gamma), is another IFN-γ inducible chemokine, stimulated by mononuclear cells to attract activated T cells and monocytes through the chemokine receptor CXCR3 (Ward and Westwick, 1998). CXCL9 expression regulates early granuloma formation induced by neutrophils through CXCL9 signaling pathway. On the other hand, wild type mice presenting anti-CXCL9 were found to be associated with impaired granuloma formation (Seiler, 2003). Following MTB infection, CXCL9 stimulates the expression of other chemokines as well as mediate migration of T- and natural killer cells migration (Lande et al., 2003). Most recently, SNP (rs2276886) within the CXCL9 (MIG) gene was associated with hypothyroidism susceptibility in Japanese populations (Akahane et al., 2016). Chemokines, such as CXCL9 are a class of cytokine proteins that act as signaling molecules, regulating immune and inflammatory responses and modulating cell migration properties and localization of target cells such as leucocytes. The chemokine CXCL9 produced by local cells in inflammatory lesions, plays important roles in several autoimmune diseases (Tokunaga et al., 2018).

In Pakistani population, there is little or no data available regarding the possible association of CXCL10 and CXCL9 polymorphisms with TB (Taheri et al., 2013). Based on the studies conducted on MTB patients from different ethnic groups, we hypothesize that genetic variation in chemokine genes CXCL10 and CXCL9 are associated with TB susceptibility. Therefore, this study was designed to investigate the risk assessment of CXCL9 (rs2276886) and CXCL10 (rs5606198) gene polymorphisms with TB in Pakistani population and also found their genetic association with different ethnic groups or subpopulations.

MATERAILS AND METHODS

Selection of study participants

All procedures performed in this study, involving human participants, were approved by the Independent Ethics Committee (IEC), International Center for Chemical and Biological Sciences, University of Karachi (ERC No: ICCBS/IEC-014-BS-2016/PROTOCOL/1.0). A total of 260 study participants, including 122 clinically diagnosed TB patients and 138 healthy control from Ojha Institute of Chest Diseases, Dow University of Health Sciences (DUHS) Karachi, Pakistan were collected. The blood samples from each subject was taken after written informed consent letter. TB patients were diagnosed and assessed by expert clinicians using their clinical history, radiological evidence, sputum acid fast bacillus (AFB) smear positivity and culture studies according to the WHO Tuberculosis disease and prevention guidelines (WHO Guidelines, 2019). All the collected samples were further categorized as: Group-1 PTB (pulmonary tuberculosis), Group-2 EPTB (extra pulmonary tuberculosis) and Group-3 HC (healthy control).

PTB patients were assessed through sputum microscopy positivity for AFB or with radiological evidence of lung involvement. EPTB patients were categorized on the basis of CT scans for abdominal, skeletal and splenic involvement and MRI for meningeal involvement. The participants with any immunosuppressive conditions such as HIV, diabetes, previous history of anti tuberculosis treatment (ATT) or any other inflammatory conditions were excluded. Healthy volunteers who have any close contact with TB patients were also excluded from the study.

Blood collection and DNA isolation

About 3 mL of blood was collected from each participant (patients and healthy controls) through venepuncture procedure in anticoagulant (EDTA) tubes for experimental analyses. DNA was isolated from each collected blood sample using standard CTAB (Cetyl trimethylammonium bromide) method (Paireder et al., 2013). Briefly, 1 mL of CTAB buffer (100 mM Tris-HCL, 20 mM EDTA, 1.4 mM NaCl, 2% CTAB, pH was adjusted as 8.0) with 5µL of β-mercaptoethanol, 20 µL (20 mg/ mL) of proteinase K enzyme and 200 µL of whole blood were added into a vial and placed in mixing block MB102 (BIOER) at 60 ºC for 1 h. Afterwards, the tubes were centrifuged for 10 min at 12,000 rpm. Following centrifugation, supernatant was transferred in to a fresh tube and 0.8 mL of chloroform/ isoamylalcohol (24:1) solution was added. The tubes were then centrifuged at 12,000 rpm for 10 min. The supernatant was transferred to the fresh tube and cold isopropanol was added to precipitate the DNA, and tubes were placed at -20ºC for 30 min. After that, tubes were centrifuged at -4°C (LABOGENE Scan speed 1730R) and supernatant was discarded. The pellet was washed with 70% ethanol, followed by centrifugation at 12,000 rpm for 5 min. Then, washing with 70% ethanol was repeated. Pellet was dried at room temperature and 100 µL of Tris-EDTA buffer was added in the tubes and stored at -20°C for further analysis. Spectrophotometry (UV-1700, Shimadzu Pharma Spec, Japan) was used to determine the yield and purity of DNA.

Genotyping

Tetra-ARMS PCR (tetra-primer amplification refractory mutation system-polymerase chain reaction) was used for rs2276886 (G/A) polymorphism of CXCL9 gene by using two sets of primers in a single reaction mixture tube. While for rs5606198 (G/A) of CXCL10 gene, Tri-ARMS PCR was applied (Taheri et al., 2013). The sequence of tetra- ARMS PCR primer sets were as following:

Pair-1 sense: 5´- ATATGCCATACATTGTGTAGC-3´

antisense: 3´- GTTCCAAGTCACTCCTGTAT-5´,

Pair-2 sense: 5´- CACTGAGAAGCTTTTATGACTAAC-3´

antisense: 3´- TCCAATATTACTCTTCAATGCAG- 5´

The resultant three genotypes of rs2276886 (G/A) were GG (408bp), AA(232bp) and for GA both bands with control band of 594bp. For tri -ARMS PCR, three primers were used as:

Forward-1 sense: 5´-TTCCTTACCTTGAATGCCACTT-3´

Forward-2sense:5´-GGAGGCTACAATAAATAATACCTTCG-3´,

Reverse antisense:3´-GGAGGCTACAATAAATAATACCTTCA-5´

The resultant genotype of rs5606198 (G/A) were GG, GA and AA with 295bp. The PCR amplification was performed in 96-well plates on a thermocycler (T100, Bio-Rad). The amplification conditions were set as the initial denaturation at 94°C for 5 minutes, followed by 35 cycles of denaturation at 94°C for 30 seconds, annealing for 30 seconds at 54°C for rs2276886 (G/A) while 62°C for rs5606198 (G/A), and extension at 72°C for 1 min. The final elongation was performed at 72°C for 5 min. The PCR products were resolved in 2% agarose gel, and the genotype patterns for each subject were recorded using the visual inspection method of the gel under the ultraviolet light by using gel documentation system (Protein Simple Fluor Chem FC3).

Statistical analysis

Allele and genotype frequencies of each polymorphism were determined by direct gene counting method. Hardy-Weinberg Equilibrium (HWE) was determined by applying the equation (p2+2pq+q2) of control population for both genetic variants (p > 0.05). Genotype distribution between case control group and TB patients were analyzed by using Chi-square analysis and Fisher’s exact tests, when the criteria for the chi-squared analysis was not fulfilled. Odd ratio (OR) was also determined by using logistic regression with their corresponding 95% CI to examine associations at the level of significance p ≤ 0.05. All the statistical analysis of the data was performed by using SPSS version 21 (SPSS Inc. Chicago, IL, USA).

RESULTS

Genotyping of the study groups

A total of 260 study participants including 122 with TB and 138 healthy control (HC) were included in present study. Out of these 122 TB disease subjects 67% were having PTB (n = 82) and 33% were with EPTB (n = 40). The mean age of the healthy controls was 30.7 ± 7 as compare to TB patients with 33.3 ± 2.1 (Table I). Genotyping was performed by T-ARMS and tri-ARMS PCR (Fig. 1A, B). By use of gene counting method, both genotypic and allelic frequencies of rs2276886 (G/A) and rs5606198 (G/A) were calculated. In HC group 85% (n=117) were having GG genotype, 14% (n = 19) were with GA genotype and 1% (n= 2) were with AA genotypes. In PTB group, 78% (n = 62) were found to have GG genotype, 20% (n= 16) were with GA genotype and 2% (n= 1) were with AA genotype. The genotypic distribution in EPTB group was 71% (n = 29) with GG genotype, 27% (n =11) with GA and 2% (n= 1) with AA genotypes, respectively (Table II).

 

Table I. Demographic characteristics of TB patients and healthy controls.

Groups	N (%)	Gender M:F (ratio)	Mean age	SD ±
HC	138 (53)	42/96	33.7	7.5
TBD	122 (47)	54/68	33.6	2.1
PTB	82 (67)	38/44	33.6	2.5
EPTB	40 (33)	16/24	32.75	1.2

 

HC, healthy control; TB, tuberculosis disease; PTB, pulmonary tuberculosis disease; EPTB, extra pulmonary tuberculosis.

 

No departure from Hardy Weinberg equilibrium was observed for both rs2276886 (G/A) of CXCL9 and rs5606198 (G/A) of CXCL10 (χ2 = 1.34, 1.47, p > 0.05). The genotypic frequency of rs506198 (G/A) was found to be extremely low (minor allele frequency < 0.01) in the studied group therefore, this SNP was not included in the further analysis. In rs2276886 (G/A), no significant difference was found between genotypic frequencies of all the groups, however there was a significant difference between allelic frequency of EPTB and HC group (χ2 = 3.3, p = 0.04). The minor allele A of rs2276886 (G/A) was found significantly associated with EPTB disease group (OR, 2.1 (0.99 to 4.30) as shown in Table I. Furthermore, to assess the genotypic effect on all the groups, different genetic models including dominant, recessive and codominant were constructed. A significant association of GG genotype with EPTB group was found (OR, 0.4 (0.181 to 0.94), p = 0.03) as compared to the PTB and HC groups in the dominant model (Table II). No significant association was found in the remaining other genetic models. Each genetic model was further stratified for gender in each group. The male with GG genotype were found to be significantly associated (OR, 0.2 (0.060 to 0.813),

 

 

Table II. Distribution of genotype and allele frequencies between the PTB, EPTB and HC groups.

Gene	Groups (n)	Genotypic frequencies (%)					Allelic frequencies (%)
		GG	GA	AA	x2	P-value	G	A	x2	P-value	OR (95%CI)
CXCL9	HC	117 (85)	19 (14)	2 (1)			253	23
	PTB	62 (78)	16 (20)	1 (2)	1.5	0.45	140	18	1.09	0.29	1.41 (0.73 to 2.71)
	EPTB	29 (71)	11 (27)	1 (2)	3.4	0.12	69	13	3.3	0.04*	2.1 (0.99 to 4.30)

 

Note: *Significance, (p < 0.05).

For abbreviations, see Table I.

 

p = 0.02) with EPTB disease as compared to the PTB and HC groups in the dominant model. While no significant association was found in other gender stratified genetic models (Tables III, IV, V).

 

Table III. Genetic analysis of rs2276886 (G/A) under the dominant model.

Group	GA+AA	GG	p-value	OR (95% CI)
Total Control PT-D EPT-D Male Control PT-D EPT-D Female Control PT-D EPT-D	20 (15) 17(22) 12 (29) 7 (16) 7 (26) 7 (47) 14 (15) 11 (21) 5 (19)	117 (85) 62 (78) 29 (71) 36 (84) 20 (74) 8 (53) 81 (85) 42 (79) 21 (81)	0.19 0.03* 0.32 0.02* 0.35 0.57	0.62 (0.304 to 1.275) 0.4 (0.181 to 0.94) 0.6 (0.170 to 1.811) 0.2 (0.060 to 0.813) 0.7 (0.275 to 1.580) 0.72 (0.234 to 2.243)

 

For abbreviations, see Table I.

 

Table IV. Genetic analysis of rs2276886 (G/A) under the recessive model.

Group	GG+GA	AA	p-value	OR (95% CI)
Total Control PT-D EPT-D Male Control PT-D EPT-D Female Control PT-D EPT-D	136 (99) 78 (99) 40 (98) 42 (98) 26 (96) 14 (93) 94 (99) 52 (98) 26 (100)	2 (1) 1 (1) 1 (2) 1 (2) 1 (4) 1 (7) 1 (1) 1 (2) 0	0.31 0.13 0.73 0.44 0.67 0.36	3.48 (0.31 to 1.08) 6.4 (0.57 to 1.21) 1.6 (0.09 to 1.95) 3 (0.17 to 1.19) 1.8 (0.11 to 2.50) 3.6 (0.21 to 1.79)

 

*Significance, (p < 0.05).

For abbreviations, see Table I.

 

In the current study, the genetic association of rs2276886 (G/A) with different ethnic groups of Pakistani population was also conducted. A regression analysis showed the significant association of minor allele (A) with Punjabi speaking ethnic group of TB disease (OR, 2.16 (1.0910 to 4.2901), p = 0.02) as compared to all other ethnic groups (Table VI).

 

Table V. Genetic analysis of rs2276886 (G/A) under the codominance model.

Group	GG+GA	GA+AA	p-value	OR (95% CI)
Total Control PT-D EPT-D Male Control PT-D EPT-D Female Control PT-D EPT-D	119(85) 78 (82) 40 (77) 42 (86) 26 (79) 14 (67) 94 (87) 52 (83) 26 (84)	22 (15) 17 (18) 12 (23) 7 (14) 7 (21) 7 (33) 14 (13) 11 (17) 5 (16)	0.6 0.22 0.41 0.71 0.42 0.65	0.8 (0.42 to 1.69) 0.6 (0.27 to 1.37) 0.7 (0.19 to 1.96) 3 (0.89 to 1.05) 0.4 (0.29 to 1.66) 1.29 (0.42 to 3.91)

 

*Significance, (p < 0.05).

For abbreviations, see Table I.

 

DISCUSSION

TB disease is generally infecting the pulmonary region, however the infection can also break into the extra-pulmonary nodules (Barr and Yates, 2017; Said et al., 2017). Previous data have shown that the intensity of cytokines in response to MTB antigens were affected from genetic polymorphisms (Qiu et al., 2001). Both the CXCL9 and CXCL10 gene SNPs are found to be associated with several diseases including TB (Tang et al., 2009; Slight and Khader, 2013). The current case control study was based on the rationale that these genes are located next to the interferon-gamma (IFN-γ) gene on chromosome 4 and are under the influence of this cytokine for their induction. Therefore, we investigated the influence of these SNPs at the promoter or regulatory region of CXCL9 and CXCL10 genes in TB patients. The focus was on risk of disease susceptibility due to these polymorphisms.

In the present study, we found the association of A allele of rs2276886 (G/A) with EPTB group which suggests that a person with A allele have two times more risk to have EPTB as compared to the normal healthy control. In concordance to our study, Nonghanphithak et al. (2016) have also demonstrated that polymorphism in chemokine genes including CXCL9 were associated with susceptibility to develop TB among the Asians. However, in Japanese and Chinese populations, the rs2276886 was found to be associated with Hashimoto’s thyroiditis (HT) (Akahane et al., 2016; Mo et al., 2019). In another study, Wilson et al. (2013) found the association of different polymorphism of CXCL9 (rs3733236), located within the untranslated region, which was associated with liver fibrosis in both the humans and the animal models. Beside these facts, very limited data was available

 

Table VI. Regression analysis of selected polymorphism for risk of tuberculosis diseases with different ethnic groups in Pakistani population.

Ethnic groups	CXCL9 tuberculosis group Allelic frequencies		Control group Allelic frequencies		OR (95% CI)	p- value
	G	A	G	A
Punjabi	72	12	39	3	2.16 (0.5766 to 8.1419)	0.25
Sindhi	34	8	9	3	0.75 (0.1549 to 3.2170)	0.65
Pashtun	47	9	25	5	0.91 (0.2895 to 3.1663)	0.94
Urdu speaking	176	30	165	13	2.16 (1.0910 to 4.2901)	0.02*
Others	15	5	5	3	0.6 (0.0962 to 3.2069)	0.51

 

*Significance, (p < 0.05).

 

regarding the association of CXCL9 SNPs with TB patients. It provides the strength to our study in particular to explore the effects of genetic predisposition of CXCR3 ligands (CXCL9 and CXCL10) on TB disease in Pakistani population. We also proposed the genetic variant (rs2276886) of CXCL9 as a potential genetic marker for TB patients.

In the present study, we found very low frequency of CXCL10 (rs5606198) SNP variant that was not included in further analysis. This low frequency might be due to the population variations with in different ethnicities. Furthermore, many cofounding factors including immunosuppressive conditions and false positive control group selection made difficult to study host genetic susceptibility to TB. In gender stratified dominant genetic models of rs2276886 (G/A) of CXCL9, we found a strong association of GG genotype with male (OR = 0.2 (0.060 to 0.813), p = 0.02) as compared to the female of EPTB group. This represents that wild type GG genotype of rs2276886 may provide protective effect to male from TB diseases susceptibility as compared to the female. The gender based influence on activation of several immune responses have been observed in various diseases and it might contribute to the disease susceptibility (Whitacre, 2001). Whether CXCL9 expression and function is influenced by sex has not been studied yet. In contrast to our findings, Biswas et al., (2020), reported another important chemokine CCL2-2518, to be significantly associated with GG genotype of female patients in South Indian population in comparison to male patients. No previous data was found with relevance to association of rs2276886 (G/A) of CXCL9 with TB. To validate our results, large size population based studies are recommended.

The importance of considering ethnicity can be determined by the observation of the gender-related case reports in the Eastern and Western neighbors of Pakistan, suggesting a difference in the epidemiology (Dogar et al., 2012). As of Karachi, being a highly populated metropolitan city of Pakistan, having people from diverse ethnic groups. So, the current study also focused on the difference in the frequency distribution of rs2276886 (G/A) in various ethnic groups. We found that Urdu speaking group with A allele carrier of CXCL9 variant rs2276886 were at two times higher risk to develop TB with reference to G allele (OR, 2.16 (1.0910 to 4.2901), p = 0.02*). Due to the unavailability of data including plasma concentration of CXCL9 and CXCL10 proteins, we were unable to explore the genotype-phenotype relationship with respect to TB disease. To validate our results, large size population based study with different ethnic groups is recommended.

CONCLUSION

In summary, we found strong association of CXCL9 polymorphism (rs2276886 G>A) with TB disease in Pakistani population. In subpopulation analysis. Punjabi speaking people have higher risk to develop TB as compared to other ethnic groups. Furthermore, in gender-based stratification, we found that male have low risk to develop TB as compared to the female. This data reinforces the critical importance of rs2276886 (G>A) of CXCL9 as a genetic marker for TB. In future, to validate our results, large population-based association studies are recommended to better understand the genotypic-phenotypic effects of chemokines that might pave a way to better control the TB susceptibility and to understand its deep pathobiology.

ACKNOWLEDGEMENTS

This work was supported by our Institute, Dr. Panjwani Center for Molecular Medicine and Drug Research (PCMD), ICCBS, University of Karachi. We thank the Institute for supporting this work by providing recurring grant for chemicals and consumables.

Declarations

Compliance with ethical standards

All procedures performed in this study, involving human participants, were approved by the Independent Ethics Committee (IEC), International Center for Chemical and Biological Sciences, University of Karachi (ERC No: ICCBS/IEC-014-BS-2016/PROTOCOL/1.0). Informed consent was obtained from all participants included in the study.

Consent for publication

All authors approved the manuscript and consented to its publication.

Statement of conflict of interests

The authors have declared no conflict of interest.

REFERENCES

Akahane, M., Watanabe, M., Inoue, N., Miyahara, Y., Arakawa, Y., Inoue, Y., and Iwatani, Y., 2016. Association of the polymorphisms of chemokine genes (IL8, RANTES, MIG, IP10, MCP1 and IL16) with the pathogenesis of autoimmune thyroid diseases. Autoimmunity, 49: 312-319. https://doi.org/10.3109/08916934.2015.1134507

Ansari, A., Talat, N., Jamil, B., Hasan, Z., Razzaki, T., Dawood, G., and Hussain, R., 2009. Cytokine gene polymorphisms across tuberculosis clinical spectrum in Pakistani patients. PLoS One, 4: e4778. https://doi.org/10.1371/journal.pone.0004778

Antonelli, A., Ferrari, M.S., Frascerra, S., Galetta, F. and Franzoni, F., 2011. Circulating chemokine (CXC motif) ligand (CXCL)9 is increased in aggressive chronic autoimmune thyroiditis, in association with CXCL10. Cytokine, 55: 288-293.

Bangs, C.D., and Donlon, T.A., 2005. Metaphase chromosome preparation from cultured peripheral blood cells. Curr. Protoct. Hum. Genet., 45: 4-11. https://doi.org/10.1002/0471142905.hg0401s45

Barr, D.A., and Yates, T.A., 2017. Whole-genome sequencing identifies nosocomial transmission of extra-pulmonary mycobacterium tuberculosis. Int. J. Med., 110: 611-612. https://doi.org/10.1093/qjmed/hcx115

Bhattacharyya, C, Majumder, P.P., and Pandit, B., 2018. CXCL10 is overexpressed in active tuberculosis patients compared to M. tuberculosis exposed household contacts. Tuberculosis (Edinb), 109: 8–16. https://doi.org/10.1016/j.tube.2018.01.005

Biswas, S.K., Mittal, M., Sinha, E., Singh, V., Arela, N., Bajaj, B. and Mohanty, K.K., 2020. 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. https://doi.org/10.1155/2020/1019639

Chibba, M., 2009. Financial inclusion, poverty reduction and the millennium development goals. Eur. J. Dev. Res., 21: 213-230. https://doi.org/10.1057/ejdr.2008.17

Cooper, A.M., Kipnis, A., Turner, J., Magram, J., Ferrante, J., and Orme, I.M., 2002. Mice lacking bioactive IL-12 can generate protective, antigen-specific cellular responses to mycobacterial infection only if the IL-12 p40 subunit is present. J. Immunol., 168: 1322–1327. https://doi.org/10.4049/jimmunol.168.3.1322

Corbett, E.L., Watt, C.J., Walker, N., Maher, D., Williams, B.G., Raviglione, M.C., and Dye, C., 2003. The Growing Burden of Tuberculosis: Global Trends and Interactions with the HIV Epidemic. Arch. Int. Med., 163: 1009–1021. https://doi.org/10.1001/archinte.163.9.1009

Dogar, O.F., Shah, S.K., Chughtai, A.A., and Qadeer, E., 2012. Gender disparity in tuberculosis cases in eastern and western Provinces of Pakistan. BMC. Infect. Dis., 12: 1-7.

Doherty, T.M., 2005. Real world TB vaccines: Clinical trials in TB-endemic regions. Vaccine, 23: 2109–2114. https://doi.org/10.1016/j.vaccine.2005.01.060

Domingo, Gonzalez, R., Prince, O., Cooper, A., and Khader, S.A., 2016. Cytokines and chemokines in mycobacterium tuberculosis infection. Microbiol. Spectr., 1: 33-72. https://doi.org/10.1128/9781555819569.ch2

Esche, C., Stellato, C., and Beck, L.A., 2005. Chemokines: Key players in innate and adaptive immunity. J. Invest. Dermatol., 125: 615–628. https://doi.org/10.1111/j.0022-202X.2005.23841.x

Faé C.K., Palacios, A.S., Nogueira, G.L., Oshiro, E.S., Demarchi, L.M.F. and Angelina, M. B., 2013. CXCL9/Mig mediates T cells recruitment to valvular tissue lesions of chronic rheumatic heart disease patients. Inflammation, 36: 800-811.

Flynn, J.L., and Chan, J., 2001. Immunology of tuberculosis. Annu. Rev. Immunol., 19: 93–129. https://doi.org/10.1146/annurev.immunol.19.1.93

Hoff, S.T., Salman, A.M., Ruhwald, M., Ravn, P., Brock, I., Elsheikh, N., Andersen, P., and Agger, E.M., 2015. Human B cells produce chemokine CXCL10 in the presence of mycobacterium tuberculosis specific T Cells. Tuberculosis (Edinb), 95: 40–47. https://doi.org/10.1016/j.tube.2014.10.005

Hussain, R., Talat, N., Shahid, F., and Dawood, G., 2009. Biomarker changes associated with tuberculin skin test (TST) conversion: A two-year longitudinal follow-up study in exposed household contacts. PLoS One, 4: e7444. https://doi.org/10.1371/journal.pone.0007444

Kotrych, D., Dziedziejko, V., Safranow, K., Drozdzik, M., and Pawlik, A., 2015. CXCL9 and CXCL10 gene polymorphisms in patients with rheumatoid arthritis. Rheumatol. Int., 35: 1319-1323. https://doi.org/10.1007/s00296-015-3234-0

Lande, R., Giacomini, E., Grassi, T., Remoli, M.E., Iona, E., Miettinen, M., Julkunen, I., and Coccia, E.M., 2003. IFN-alpha beta released by mycobacterium tuberculosis-infected human dendritic cells induces the expression of CXCL10: Selective recruitment of NK and activated T cells. J. Immunol., 170: 1174–1182. https://doi.org/10.4049/jimmunol.170.3.1174

Loddenkemper, R., Lipman, M., and Zumla, A., 2015. Clinical aspects of adult tuberculosis. Cold Spring Harb. Perspect. Med., 6: a017848. https://doi.org/10.1101/cshperspect.a017848

Lu, H., Wang, J., Gao, B., Chen, S., Xiang, G., Yang, S., Gong, Y., and Guan, L., 2015. The association between the CC chemokine ligand 5-28 C> G gene polymorphism and tuberculosis susceptibility. Saudi med. J., 36: 1400. https://doi.org/10.15537/smj.2015.12.12264

Magri, M.C., Alvarez, M.S., Iogi, A.A., Alves, G.M., Manchiero, C., Dantas, B.P., Prata, T.V., Nunes, A.K., and Tengan, F.M., 2021. Study of CXCL9-11 gene polymorphisms in liver fibrosis among patients with chronic hepatitis C. Pathog. Dis., 79: ftab007. https://doi.org/10.1093/femspd/ftab007

Maher, D., Chaulet, P., Spinaci, S., and Harries, A., 1997. Treatment of tuberculosis: Guidelines for national programmes. Second edition., (Ed. 2), 1-77.

Metzger, P., Baloch, N.A., Kazi, G.N., and Bile, K.M., 2010. Tuberculosis control in Pakistan: Reviewing a decade of success and challenges. East Mediterr. Hlth. J., 16(Supp.): 47-53. https://doi.org/10.26719/2010.16.Supp.47

Mo, D., Li, J., Peng, L., Liu, Z., Wang, J., and Yuan, J., 2019. Genetic polymorphisms on 4q21. 1 contributed to the risk of Hashimoto’s thyroiditis. Genet. Test. Mol. Biomark., 23: 837-842. https://doi.org/10.1089/gtmb.2019.0125

Nogueira, L.G., Santos, R.H., Ianni, BM., Fiorelli, A.I, Mairena, E.C., Benvenuti, L.A., Frade, A., Donadi, E., Dias, F., Saba, B., and Wang, H.T., 2012. Myocardial chemokine expression and intensity of myocarditis in chagas cardiomyopathy are controlled by polymorphisms in CXCL9 and CXCL10. PLoS. Negl. Trop. Dis., 6: e1867. https://doi.org/10.1371/journal.pntd.0001867

Nonghanphithak, D., Reechaipichitkul, W., Namwat, W., Lulitanond, V., Naranbhai, V., and Faksri, K. 2016. Genetic polymorphisms of CCL2 associated with susceptibility to latent tuberculous infection in Thailand. Int. J. Tuberc. Lung Dis., 20: 1242-1248. https://doi.org/10.5588/ijtld.16.0017

Nonghanphithak, D., Reechaipichitkul, W., Namwat, W., Wongwajana, S., Lulitanond, V., and Faksri, K., 2015. Comparison of CXCL9 polymorphism between pulmonary tuberculosis patients and healthy controls in Northeast Thailand. Srinagarind Med. J., 30: 432–438.

Paireder, S., Werner, B., Bailer, J., Werther, W., Schmid, E., Patzak, B., and Cichna-Markl, M., 2013. Comparison of protocols for DNA extraction from long-term preserved formalin fixed tissues. Anal. Biochem., 439: 152-160. https://doi.org/10.1016/j.ab.2013.04.006

Qiu, B., Frait, K.A., Reich, F., Komuniecki, E., and Chensue, S.W., 2001. Chemokine expression dynamics in mycobacterial (type-1) and schistosomal (type-2) antigen-elicited pulmonary granuloma formation. Am. J. Pathol., 158: 1503-1515. https://doi.org/10.1016/S0002-9440(10)64101-6

Said, A.M., Shehata, S.M., Lutfy, S.M., Elsammak, A., E.l-Azony, A., and Roshdy, M., 2017. Assessment of radiological pulmonary shadows in patients with extra pulmonary neoplasm. Egypt J. Chest Dis. Tuberc., 66: 681–686. https://doi.org/10.1016/j.ejcdt.2017.09.001

Seiler, P., Aichele, P., Bandermann, S., Hauser, A.E., Lu, B., Gerard, N.P., Gerard, C., Ehlers., S., Mollenkopf, H.J., and Kaufmann, S.H., 2003. Early granuloma formation after aerosol mycobacterium tuberculosis infection is regulated by neutrophils via CXCR3-signaling chemokines. Eur. J. Immunol., 33: 2676–2686. https://doi.org/10.1002/eji.200323956

Sheikh, G., Neela, V.S., Pydi, S.S., Suryadevara, N.C., Gaddam, R., Gaddam, S.L., Auzumeedi, S.K., and Valluri, V.L., 2015. Genetic association of interferon gamma induced protein-10 (IP-10), CXCL-10 gene polymorphisms with TB pleurisy susceptibility in south Indian population. Open J. Immunol., 5: 72–78. https://doi.org/10.4236/oji.2015.52008

Slight, S.R., and Khader, S.A., 2013. Chemokines shape the immune responses to tuberculosis. Cytokine Growth Factor Rev., 24: 105-113. https://doi.org/10.1016/j.cytogfr.2012.10.002

Strzelak, A., Komorowska, Piotrowska, A., and Ziołkowski, J., 2012. CXCL10/IP-10 as a new biomarker for Mycobacterium tuberculosis infection. Pol. Merkur. Lekarski., 33: 342–345.

Taheri, M., Kouhpayeh, H.R., Hosseinalizadeh, T., Naderi, M., Bahari, G., and Hashemi, M., 2013. A functional polymorphism in promoter of the CXCL10 Gene (-135 G/A) associated with pulmonary tuberculosis. Arch. Clin. Infect. Dis., 8: e15542. https://doi.org/10.5812/archcid.15542

Tang, N.L., Fan, H.P., Chang, K.C., Ching, J.K., Kong, K.P., Yew, W.W., Kam, K.M., Leung, C.C., Tam, C.M., Blackwell, J., and Chan, C.Y., 2009. Genetic association between a chemokine gene CXCL-10 (IP-10, interferon gamma inducible protein 10) and susceptibility to tuberculosis. Clin. Chim. Acta, 406: 98–102. https://doi.org/10.1016/j.cca.2009.06.006

Tokunaga, R., Zhang, W., Naseem, M., Puccini, A., Berger, D.M., Soni, S., Michelle, M.M., Baba, H. and Lenz, H.J., 2018. CXCL9, CXCL10, CXCL11/CXCR3 axis for immune activation - A target for novel cancer therapy. Cancer Treat. Rev., 63: 40-47.

Turner, M.D., Nedjai, B., Hurst., and Pennington, D.J., 2014. Cytokines and chemokines: At the crossroads of cell signalling and inflammatory disease. Biochim. biophys. Acta, 1843: 2563–2582. https://doi.org/10.1016/j.bbamcr.2014.05.014

Ward, S.G., and Westwick, J., 1998. Chemokines: Understanding their role in T- lymphocyte biology. Biochem. J., 333: 457–470. https://doi.org/10.1042/bj3330457

Whitacre, C.C., 2001. Sex differences in autoimmune disease. Nat. Immunol., 2: 777-780. https://doi.org/10.1038/ni0901-777

WHO guidelines on tuberculosis infection prevention and control, 2019 update, World Health Organization; Geneva. License: CC BY-NC-SA 3.0 IGO.80

Wilson, N., Driss, A., Solomon, W., Dickinson-Copeland, C., Salifu, H., Jain, V., and Stiles, J., 2013. CXCL10 gene promoter polymorphism-1447A> G correlates with plasma CXCL10 levels and is associated with male susceptibility to cerebral malaria. PLoS One, 8: e81329. https://doi.org/10.1371/journal.pone.0081329

Winslow, G.M., Cooper. A., Reiley, W., Chatterjee, M., and Woodland, D.L., 2008. Early T-cell responses in tuberculosis immunity. Immunol. Rev., 225: 284–299. https://doi.org/10.1111/j.1600-065X.2008.00693.x

WHO, 2005. Tuberculosis control in the South-East Asia and Western Pacific Regions 2005. A bi-regional report, World Health Organization.