Footprintsofgeneticsusceptibilitytopulmonarytuberculosis: CytokinegenevariantsinnorthIndians Abhimanyu,MridulaBose,PankajJha * &IndianGenomeVariationConsortium Department of Microbiology, Vallabhbhai Patel Chest Institute, University of Delhi & * Genomics & Molecular Medicine, CSIR-Institute of Genomics & Integrative Biology, Delhi, India ReceivedMarch1,2011 Background & objectives: Tuberculosis is (TB) responsible for high morbidity and mortality worldwide. Cytokines play a major role in defense against Mycobacterium tuberculosis infection. Polymorphisms in the genes encoding the various pro- and anti-inammatory cytokines have been associated with tuberculosis susceptibility. In this study we examined association of 25 sequence polymorphisms in six candidate cytokine genes namely IFNG, TNFB, IL4, IL1RA, IL1B and IL12 and their related haplotypes with risk of developing pulmonary tuberculosis (PTB) among north Indians. Methods: Pulmonary TB (n=110) patients and 215 healthy controls (HC) from north India were genotyped. Puried multiplex PCR products were subjected to mass spectrometry using Sequenom MassARRAY platform to generate the genotypes in a population-based case-control study. Results: Using multiple corrections, signicant overall risk against PTB was observed at seven loci which included variants in IFNG at rs1861493 and rs1861494; IL1RA at rs4252019, IL4 variant rs2070874, IL12 variants rs3212220, rs2853694 and TNFB variant rs1041981. Analysis of gene structure revealed two haplotype blocks formed by IFNG variants rs1861493 and rs1861494. The TA haplotype was signicantly over-represented (P=0.011) in the cases showing a two-fold risk in the current population (Odds ratio=1.59 CI=1.101 to 2.297) and TNFB variants at rs2229094 and rs1041981 contributed to two haplotypes which were in strong linkage disequilibrium (LD) with AT haplotype showing a three-fold risk (P=0.0011, Odds ratio=3, CI=0.1939 to 0.7445) of developing PTB in north Indians. Interpretation & conclusions: Our study showed six novel associations of cytokine gene variants with susceptibility to PTB in north Indians. Variants of IFNG and TNFB emerged as factors imposing a signicant risk of developing PTB in north Indians apart from risk indicated by IL1RA, IL4 and IL12. Key wordsCytokinegenevariant-haplotype-Mycobacterium tuberculosis-pulmonarytuberculosis-singlenucleotidepolymorphisms 763 IndianJMedRes135,May2012,pp763-770 Tuberculosis (TB) causes signicant morbidity and mortality throughout the world 1 . The vast majorityof individualsinfected withMycobacterium tuberculosis (up to 95%) remain healthy, probably because of mounting an effective immune response against M. tuberculosis. In 1949, Haldane proposed that the maintenance of multiple genes that confer relative susceptibilities on the host to infectious diseases wouldbe favoured byevolution. Insupport of this hypothesis, certain populations appear to be 764 INDIANJMEDRES,MAY 2012 at risk for both increased susceptibility to infection 2 andprogressiveclinicaldiseaseduetomycobacteria 3 . Severalcase-controlstudieshaveidentiedassociation betweenTBandcandidategenespotentiallyinvolved in immune response to TB 4,5 . A growing body of evidencesupports arole ofhost geneticcomponents inthedevelopmentoftuberculosis.Theobservationof familialclusteringofdiseasewithhigherconcordance oftuberculosisdiseaseinmonozygoticversusdizygotic twins 6 , the ethnic clustering of tuberculosis disease withahigherprevalenceoftuberculosisinindividuals ofrecentAfricandescent 2 ,aswellasthedemonstration of both common polymorphisms and rare mutations which confer susceptibility to mycobacterial species inhumans 7 pointsignicantlyinthisdirection.These studies suggest that unique environment and natural selectivefactorsmayberesponsibleforthedevelopment ofethnic-specichostgeneticfactorsassociatedwith TB. The rst step in innate host defense is cellular uptake of M. tuberculosis, which involves different cellularreceptorsandhumoralfactors.Thesubsequent inammatoryresponseisregulatedbytheproductionof pro-andanti-inammatorycytokinesandchemokines. Interferon-gamma (IFN-γ one of the most important cytokines involved in macrophage activation, stimulatinganti-tumourandanti-microbicidalactivities as well as expression of MHC-II 8,9 . Interleukin-4 (IL-4), an anti-inammatory cytokine has been implicated to downregulate IFN-γ, and thus has a deleterious effect on TB patients 10 . It also promotes the induction of Th2 cells 11 . IL-12, a heterodimeric pro-inammatory cytokine produced by activated macrophages,monocytes,β-lymphocytesanddendritic cellsistheprincipalTh1responseinducingcytokine 11 . Thiscytokine isimportant forsustaining asufcient number of memory/effector Th1 cells to mediate long-term protection to intracellular pathogen. Like tumournecrosisfactor-alpha(TNF-a),IL-1bismainly produced by monocytes, macrophages, anddendritic cells 12 . Intuberculosispatients,IL-1bisexpressedin excess 13 andatthesiteofdisease 14 .Implicatedmainly intuberculosis pleurisy, ausually self-resolvingtype ofprimarytuberculosis,onemayhypothesizethatan increasedIL-1b/IL-1Raratioprotectsagainstamore severeformoftuberculosis. TNF-b or lymhotoxin-alpha (LTa) is considered to be a proinammatory cytokine and it is shown that secreted LTa is essential for the control of an intracellularbacterialinfection 15 .RecentlyAllieet al 16 suggestedthatLTαmightnothaveacriticalroleinhost defensetoacutemycobacterialinfection,independent ofTNF,butcertainlyacontributionofLTαinthecontrol ofchronic M. tuberculosis infectionisobserved 17 . Association studies from north India probing multiplelociacrossthespectrumofcandidatecytokine genesarescanty.Thepresentstudy,therefore,wasaimed to bring in focus certain unexplored polymorphisms in the context of tuberculosis susceptibility in north Indianpopulation.Theroleandimportanceofgenetic backgroundintuberculosishasnowbecomeunivocal withethnicityplayingacrucialrole.Probingnewloci relating to tuberculosis susceptibility could suggest novel approach in pharmacogenomics and therapy to combat this pathogen. Also it could provide an insightinto predictingindividual’sgeneticproneness totuberculosisandofbeingfuturediagnostictoolfor preventivetherapyagainsttuberculosis. Material & Methods Study population: PTB patients above 18 yr of age (n=110)wereenrolledrandomlyinthestudybetween 2010-11 from Rajan Babu Institute of Pulmonary Medicine and Tuberculosis (RBIPMT), Kingsway Camp,NewDelhi(India).Thestudywascarriedout in Department of Microbiology, V.P. Chest Institute, University of Delhi, Delhi. Enrolled patients were category I cases, clinically and radiologically (chest X-ray) diagnosed for pulmonary tuberculosis and conrmed by sputum microscopy and culture for Mycobacterium following the guidelines of Revised NationalTBControlProgramme(RNCTP),Ministry of Health and FamilyWelfare, Government of India (http://www.tbcindia.nic.in). All patients were given free anti-tuberculosis drugs under DOTS (Directly Observed Treatment, short course) regimen of the GovernmentofIndia.ThemeanageofPTBcaseswas 31.89 ± 2.6yr while the ratioof male : femalewas 47:53. Patients having any immunosuppressive presentation such as diabetes mellitus or HIV co- infection which areconsidered to be risk factorsfor tuberculosis development, and patients suspected to have extra-pulmonary tuberculosis along with pulmonarytuberculosiswereexcludedfromthestudy. Structured questionnaires were usedto document all otherrelevantinformationsuchasage,sex,ethnicity, socio-economic status, BCG vaccinations, and previous family history of tuberculosis. The healthy control(HC)groupconsistedof215randomlychosen nonconsanguineous BCG vaccinated students and laboratorypersonnelfromthevariousdepartmentsof theUniversityofDelhiwhowerewillingtoparticipate in the study with no signs, symptoms or history of previous mycobacterial infection. For HC mean age was29.31±.82yrandtheratioofmale:femalewas 43:57. Analysis of population stratication: Serious effort was made to avoid any false-positives arising as a result of population stratication. The self reported ethnicity of each subject and his/her parents was carefully considered. In addition, the genotype data weresubjectedtoEIGENSTRATprincipalcomponent analysis for population stratication correction as illustratedbyPriceet al 18 . Allindividualswerebriefedaboutthestudyanda signedinformedconsentwasobtainedfromthepatient or his or her guardians before sample collection. The study was approved by the ethics committee of VallabhbhaiPatelChestInstitute,UniversityofDelhi, India. DNA extraction: Three ml of venous blood was collected in BD vacutainers containing ethylene diaminetetraaceticacid(EDTA)asanticoagulantand kept frozen until use. Genomic DNA was extracted from frozen whole blood using QiaAMP DNA kit (Qiagen,Germany).ExtractedDNAwasquantiedby spectrophotometery, checkedfor purityand storedat -20 o Cuntilfurtheranalyses. SNP selection and genotyping:Sixcandidatecytokine genes namely IFNG, TNFB, IL4, IL1RA, IL1B and IL12B, were selected owing to their suggested role in tuberculosis pathogenesis. All single nucleotide polymorphisms (SNPs) selected for genotyping were accessed from the public dbSNP (http://www. ncbi.nih.gov) and the HapMap (http://www.hapmap. org/).MostoftheselectedSNPsarefromtheintronic regionsofthecorrespondinggenes.Wereasonedthat notonlythechangesinthepromoterbutalsoofother unexploredregionsofthegenemayhamperitsnormal functioningleading todisease.The parameterstaken intoaccountwhileSNPselectionwerethefrequency of<0.01indbSNP,reportedallelefrequencyofatleast 20percentintwoworldpopulations(fromHapmap), average spacing 1 kb but in closely spaced minor allelefrequencywascarefullyconsidered.Inaddition, reportedheterozygositywasconsideredinaneffortto minimizeselectionofhomozygousloci. All SNPs were genotyped using the matrix- assisted laser desorption/ionization time-of-ight (MALDI-TOF) mass spectrometry (Sequenom Inc., USA). Assays for all SNPs were designed using SpectroDESIGNER software (Sequenom Inc., USA) All SNPs were genotyped using the iPLEX assays (www.sequenom.com/iplex). Briey, as template, 5 ng of genomic DNA was used in a multiplex PCR reaction.ThePCRproductwasfurtherpuriedbefore theprimerextensionreactiontogenerateallele-specic baseextensionproducts.Thebase-extensionproducts weredetectedintheMALDI-TOFmassspectrometer to determinegenotypes. Genetic and statistical analyses: Hardy-Weinberg equilibriumwascalculatedinbothPTBcasesandHC separatelytoensurethatthesampleswerewithinallelic populationequilibriumbyusingHaploviewv4.2(http:// www.broad.mit.edu/mpg/haploview/).A stringent cut- offofferedbytheHaploviewv4.2wasusedforfurther analysis (minimum genotype =75% and minimum minor allele frequency 0.0010). The samples and variationsfailingthistestwerenotselectedforfurther analysis.PLINKv1.07(http://pngu.mgh.harvard.edu/ purcell/plink/)wasusedtotestformultiplecomparison andPvalueafterBonferronicorrectionswasconsidered signicant.Haplotypeblockgenerationwasperformed using the algorithm by Gabriel et al 19 implemented in the Haploview software which was also used for initial association testing. The statistical signicance ofPvalueofhaplotypeswasassessedbypermutation analysis(N=10,000)withHaploviewv4.2. Geneticassociationtestingwasdoneusinga2x2 contingencytable.Oddsratio,twotailedPvaluewas calculatedforalleles.2x2Computationsweredone using GraphPad Prism (version 5.00 for Windows, Graph Pad Software, San Diego California, USA; www.graphpad.com). Two-tailed P<0.05 was consideredstatisticallysignicant. Results Table Ishows thelocation and characteristicsof theSNPsincludedinthestudyandTableIIshowsthe associationsaftermultiplecorrectionscarriedoutusing PLINK (http://pngu.mgh.harvard.edu/purcell/plink/) whichwerefoundtobeassociatedwithsusceptibility toPTBinnorthIndiansinthisstudy. Population stratication correction: To access any underlyingstructureinthestudypopulationthatcould ABHIMANYUet al:CYTOKINEGENEVARIANTSINPTB 765 Table II.Allelicassociationsinafteradjustmentformultipletesting Gene dbSNP a rsID Case (n=110), control (n=215)frequencies Oddsratio (95%CI) Chisquare Pvalue * P bonferroni # IFNG rs1861493 0.962,0.869 3.8(1.7-8.6) 12.089 5.00E-04 0.00659 rs1861494 0.946,0.859 3.0(1.5-5.6) 10.466 0.0012 0.01581 IL4 rs2070874 0.387,0.255 1.8(1.3-2.6) 10.708 0.0011 0.01387 TNFB rs1041981 0.356,0.238 1.7(1.2-2.6) 8.649 0.0033 0.03618 IL12 rs2853694 0.607,0.478 1.6(1.2-2.4) 8.854 0.0029 0.0399 rs3212220 0.419,0.263 2.0(1.4-2.9) 14.572 1.00E-04 0.00175 IL1RA rs4252019 1.000,0.935 14.0(1.8-103.5) 13.643 2.00E-04 0.00287 SNP,singlenucleotidepolymorphism; * unadjustedP-value; # Pvalueafterbonferronimultipletestingcorrection; a dbSNP,theSNP database(http://www.ncbi.nlm.nih.gov/projects/SNP);P<0.05wasconsideredsignicant Table I. Locationandbase-pairpositionsofsingleneucleotidepolymorphisms(SNPs)ofvariouscytokinegenespassingtheexclusion criteriaandminorallelefrequency(MAF)incontrols Gene name dbSNP a rsID Basechange Chromosome position Location MAF controls References IFNG rs1861493 A/G 68551196 Intron4 0.13 New;thisstudy rs1861494 C/T 68551409 Intron4 0.14 New;thisstudy IL4 rs2070874 C/T 132009710 5’-UTR 0.25 Molleret al,2010 20 TNFB rs1041981 A/C 31540784 Exon4 0.26 New;thisstudy IL12 rs2853694 A/C 158749088 Intron4 0.5 New;thisstudy rs3212220 G/T 158754195 Intron4 0.28 Molleret al,2010 20 IL1RA rs4252019 C/T 113889119 Intron5 0.05 New;thisstudy a dbSNP,theSNPdatabase(http://www.ncbi.nlm.nih.gov/projects/SNP) confoundtheapparentgeneticassociationpopulation straticationcorrectionwascarriedoutusingEigenstrat PrincipalComponentanalysismethodasillustratedby Priceet al 18 .The methodmodels ancestrydifference between cases and controls and any other compared groupbasedonthesuppliedgenotypedata.Ourcases andcontrolsformedahomogenousgroupdevoidofany stratication.According to Indian Genome Variation Consortium (IGVC) 20 north Indians fall into Indo- European lineage. Our cases and controls matched withsuppliedmarkerdataofIndo-Europeanancestry therebyrulingoutcompletelyanyunderlyingstructure inthepopulation. Allelic association of cytokine SNPs and the risk of pulmonary tuberculosis:Amongthe25studiedSNPs, from six candidate cytokine genes the variants of IFNG, IL1RA, IL4, IL12 andTNFBwerefoundtobe associatedwithsusceptibilitytoPTBinnorthIndians. Allstudiedvariantspassingtheexclusioncriteriawere in Hardy-Weinberg equilibrium in both cases and controls.Allelic associationwhen probed invariants passingtheexclusioncriteriayieldedsixlocishowing highriskforPTBsusceptibility. IFNG polymorphism and PTB susceptibility: After adjusting for multiple testing corrections the IFNG intronic variants at rs1861493 [χ 2 =12.089, P bonferroni =0.006593,oddsratio(95%CI)=3.8(1.7-8.6)] andrs1861494(χ 2 =10.466,P bonferroni =0.01581,odds ratio(95%CI)=3.0(1.5-5.6)]showedasignicantrisk ofdevelopingpulmonarytuberculosisinnorthIndians with over-representation of the associated A and T allelesamongPTBpatients,respectively.Investigation of the gene structure and linkage disequilibrium patternshowed haplotypesformed byIFNG variants rs1861493andrs1861494whichwereinhighlinkage disequilibrium (LD) (Fig.). Three combinations of haplotypewereseennamelyTC,CCandTA,ofwhich TA haplotype was over-represented in the cases and imposed a two-fold risk of developing pulmonary tuberculosisinnorthIndians(TableIII). 766 INDIANJMEDRES,MAY 2012 Table III.Heplotypeblocksandfrequencies Blocks Haplotype frequency Case(n=110), control(n=215) frequencies Chisquare Permutations Pvalue # Oddsratio (95%CI) Block1 TC 0.41 0.36,0.43 2.82 0.093 0.75 (0.53-1.1) CC 0.32 0.29,0.33 0.86 0.353 0.85 (0.59-1.2) TA 0.27 0.34,0.24 6.46 0.04 1.59 (1.1-2.3) Block2 AT 0.88 0.95,0.86 10.66 0.005 2.9 (1.5-5.6) GC 0.10 0.05,0.13 8.85 0.017 0.38 (0.2-0.7) # Pvalueafterperformingpermutation(n=10,000);P<0.05wasconsideredsignicant Fig. Linkagedisequilibrium(LD)plotandhaplotypestructureof cytokinegenevariantsinPTBcases.D’valuesaredisplayedwithin eachdiamond,missingvaluesindicateD’=100%.Colourscheme gradientindicatesr 2 values.Lengthofeachblock,inkilobases(kb), isshowninbrackets. IL4 polymorphism and PTB susceptibility:IL4variant rs2070874[x 2 =10.708,P bonferroni =0.01387,oddsratio (95%CI)=1.8(1.3-2.6)]showedatwo-foldriskbyT alleleinnorthIndians.TheotherstudiedIL-4variant rs2243270passingtheexclusioncriteriadidnotshow any association towards susceptibility to pulmonary tuberculosisinthispopulation. IL1RA polymorphism and PTB susceptibility: The signicantlyassociatedlocusofIL1RAincludedintronic variantatrs4252019[χ 2 =13.643,P bonferroni =0.00287, Oddsratio(95%CI)=14.0(1.8-103.5)]showinga14- foldrisk.Othervariantsuchasrs315919andrs380092 didnotshowanyassociationtowardssusceptibilityto pulmonarytuberculosisinthispopulation. IL12 polymorphism and PTB susceptibility: IL12 variants rs3212220 [χ 2 =14.572, P bonferroni = 0.00175, Oddsratio(95%CI)=2.0(1.4-2.9)]andrs2853694 [χ 2 =8.854,P bonferroni = 0.0399,odds ratio(95%CI) = 1.6(1.2-2.4)]showedatwo-foldriskassociatedwith TandAalleles,respectively. IL1B polymorphism and PTB susceptibility: The selectedIL1Bvariantsdidnotshowanydirectinuence onPTBsusceptibilityinnorthIndians. TNFB polymorphism and PTB susceptibility: TNFB variantsatrs1041981[χ 2 =8.649,P bonferroni =0.03618, Oddsratio(95%CI)=1.7(1.2-2.6)]asynonymous changeshowedatwo-foldriskofassociationforPTB innorthIndians. Interestinglyrs1041981 contributed to a haplotype block with rs2229094 conrming the importanceofthislocusinriskofdevelopingPTBin northIndians.ThetwohaplotypesobservedwereAT and GC of which AT was over-represented in PTB casesandimposedathree-foldriskofdevelopingPTB innorthIndians. Discussion Thehostgeneticbiascontributingtosusceptibility and progression of pulmonary tuberculosis might ABHIMANYUet al:CYTOKINEGENEVARIANTSINPTB 767 involve interactions between multiple alleles located on different genes and chromosomes 21 . In order to overcome this drawback we planned selection of differentcytokinegeneandmultiplelocitocoverawide spectrumofimmuneresponseassociatedcytokines. Case-control studies involving carefully chosen locus across ethnicities are valiant means of identifying novel associations pertaining to disease susceptibility.Associationthatarisesmaybe aresult ofthepolymorphisminquestionbeingfunctionalorit beinginlinkagedisequilibriumwithanotherfunctional allele or a result of confounding association due to population stratication. To overcome such false positives, we carefully considered the self reported ethnicityofthestudygroupsandfurthercheckedfor any genetic heterogeneity in our data by Eigenstrat principal component analysis illustrated by Price et al 18 andfoundthatthepresentdatawerefreefromany underlying population structure. Thus, this uniform datarepresentnorthIndianpopulationforassociation analysis. The IFN-γ being a crucial cytokine in immunopathogenesisofTBhasbeensubjecttoseveral polymorphisms studies for pulmonary tuberculosis susceptibility. The locus probed here namely rs1861494 has not been studied in susceptibility to PTB but extensively studied in many other diseases suchasleprosy 22 andasthma 23 .Kumaret al 24 foundan associationofthislocuswithsusceptibilitytoasthma inIndians andcould identifya haplotype.They also showed that alleles of rs1861494 A/G have differential afnitytobindtoputativenuclearfactor.Inthepresentstudy, wefoundsignicantriskforthelocusinsusceptibility toPTB.The otherprobed locusrs1861493 hasbeen studied in idiopathic inammatory myopathy 24 and asthma 23 but notin pulmonary tuberculosis.We also identied a risk haplotypecontributed by rs1861493 and rs1861494 emphasizing the importance of the above mentioned loci as risk factors for developing pulmonarytuberculosisinnorthIndians. IL4locusrs2070874hasbeenanimportantlocus ofinvestigationinvariousdiseasesincluding asthma andrheumatoidarthritis 25 .ItsroleinTBwasreported nottobesignicantinIranianpulmonaryTBpatients 26 and recently in South Africans TB patients also the locus did not show any association 27 . In the present study this locusshowed a two-fold riskin the north Indianpopulation. IL1RAlocusrs4252019hasshownsignicantrisk ofdevelopmentofpulmonaryTBinnorthIndians.The variantrs4252019hasbeenshowntobeassociatedwith prostatecancerrisk 28 butnotpulmonarytuberculosis. Interestingly, the variant showed a 14-fold risk of developing PTB in the population studied here and emerged as a major locus to look out for in further studies. IL12 variantsrs3212220andrs2853694showeda signicantriskassociatedwithdevelopmentofPTBin northIndians.Thevariantrs321220hasbeenshownto contributetoahaplotypebyMolleret al 20 .Wehavealso predicteditsimportanceinourpreviousstudy 29 .Based ontheanalysisofserumIL-12level,wedemonstrated that forIL12 variant rs3212220TTgenotype among active PTB cases showed signicantly higher serum IL-12levelwhencomparedtoeitherGTorGG.The present study revealed T allele to be a risk allele in the present population. Similarly, rs2853694 a novel variantinthecontextofdevelopingtuberculosis 29 was predictedtobeofimportanceandwasvalidatedinthe presentstudy.Forrs2853694amongactivePTBcases AA genotype showed a trend towards higher serum IL-12levelincontrasttoareversetrendobservedin HCwhereAAaccountedforlowserumIL-12 29 .The presentstudyshowedAalleleatrs2853694tobearisk allele for the north Indian population in the context of PTB susceptibility. An interesting observation was that both the higher serum cytokine producers i.e.TTgenotypeforrs3212220andAAgenotypefor rs2853694 emerged as respective risk alleles T and A for thispopulation, indicating that overproduction of IL-12 by these individuals might be interfering withthe cytokinehomeostasis andthus affectingthe immunefunctionofthecytokineintheseindividuals makingthempronetoinfection.Ourobservationwas furthersupportedbytheworkofLeandroet al 30 ,who indicatedthatroleofIL-12aspotentinducerofIFN-γ liedinitsefcacyatlowconcentrations.Inthepresent studyitis observedthat thePTB patientswith IL12 riskallelegenotypesarenotefcientinducersofIFN-γ whichinturninterfereswiththeprotectiveimmunity intheseindividuals,whereasalowproleofIL-12in HCelicitsaneffectiveandoptimalimmuneresponse renderingtheseindividualshealthy. TNFB though not usually considered for PTB associationstudies,wastakenupinthecurrentstudy becauseofitsroleincontrolofintracellularbacterial infection 15 . The variant rs1041981 emerged as a 768 INDIANJMEDRES,MAY 2012 signicantrisk locusfor PTBsusceptibility innorth Indians. The variant also contributed to a haplotype with rs2229094 and reinstated the role of TNFB polymorphismsinPTB. Overall, ve of the loci namely rs1861493 and rs1861494 (IFNG), rs4252019 (IL1RA) rs1041981 (TNFB) and rs2853694 (IL12) studied in patients of pulmonary tuberculosis showed a signicant risk towards susceptibility to pulmonary tuberculosis in northIndians.Wealsoreportherethesignicantrisk imposedbyIL4 variantrs2070874 intheactivePTB patients.Sixnewassociationsandthreenewassociated haplotypescontributingtothespectrumofcytokinegene polymorphismsandriskofdevelopingtuberculosisin generalandnorthIndiansinparticular,weredetected. Acknowledgment Theauthorsthankallpatientsandvolunteersforparticipating inthisstudy.ThesupportoftheMedicalSuperintendentandstaff atRajanBabuInstituteofPulmonaryMedicineandTuberculosis (RBIPMT), Kingsway Camp,New Delhi(India) forthe helpin sample collection is acknowledged. Authors acknowledge the CouncilofScienticandIndustrialResearch(CSIR),NewDelhi, for nancial support.The rst author was the Junior Research Fellow(JRF)intheCSIRproject. References 1. RaviglioneMC,SniderJrDE,KochiA.Globalepidemiology of tuberculosis: morbidity and mortality of a worldwide epidemic.JAMA1995;273:220-6. 2. Stead WW, Senner JW, Reddick WT, Lofgren JP. Racial differences in susceptibility to infection by Mycobacterium tuberculosis.N Engl J Med 1990;322:422–7. 3. Bellamy RJ, HillAV. Host genetic susceptibility to human tuberculosis. Novartis Found Symp1998;217:3-13. 4. NewportMJ,HuxleyCM,HustonS, HawrylowiczCM, Oostra BA, Williamson R, et al. A mutation in the interferon- γ-receptorgeneandsusceptibilitytomycobacterialinfection. N Engl J Med 1996;335 :1941-5. 5. Altare F, DurandyA,Lammas D, Emile JF, Lamhamedi S, LeDeistF,et al.Impairmentofmycobacterialimmunityin humaninterleukin-12receptordeciency.Science1998;280: 1432-5 6. Comstock GW. Tuberculosis in twins: a re-analysis of the Prophitsurvey.Am Rev Respir Dis1978;117:621-4. 7. Dofnger R, Dupuis S, Picard C, Fieschi C, Feinberg J, Barcenas-MoralesG,et al.InheriteddisordersofIL-12and IFNgamma-mediatedimmunity:amoleculargeneticsupdate. Mol Immunol2002;38:903-9. 8. Pace JL, Russell SW, Torres BA, Johnson HM, Gray PW. Recombinantmouseγ-interferoninducestheprimingstepin macrophageactivationfortumorcellkilling.J Immunol1983; 130:2011-3. 9. NathanCF,MurrayHW,WiebeME,RubinBY.Identication of interferon-γ as the lymphokine that activates human macrophageoxidativemetabolismandanti-microbialactivity. J Exp Med1983;58:670-89. 10. PowrieF,CoffmanRL.Inhibitionofcell-mediatedimmunity byIL4andIL10.Res Immunol 1993;144:639-43. 11. Abbas AK, Murphy KM, Sher A. Functional diversity of helperTlymphocytes.Nature 1996;383:787-93. 12. RoachTI,BartonCH,ChatterjeeD,BlackwellJM.Macrophage activation: lipoarabinomannan from avirulent and virulent strainsofMycobacterium tuberculosis differentiallyinduces theearlygenesc-fos,KC,JE,andtumornecrosisfactor-alpha. J Immunol 1993;150:1886-96. 13. Schauf V, Rom WN, Smith KA, Sampaio EP, Meyn PA, TramontanaJM,et al.Cytokinegeneactivationandmodied responsiveness to interleukin-2 in theblood of tuberculosis patients.J Infect Dis 1993;168:1056-9. 14. Law K,Weiden M, HarkinT,Tchou-WongK, Chi C, Rom WN.Increasedreleaseofinterleukin-1beta,interleukin-6,and tumornecrosisfactor-alphabybronchoalveolarcellslavaged frominvolvedsitesinpulmonarytuberculosis.Am J Respir Crit Care Med 1996;153:799-804. 15. RoachDR,BriscoeH,SaundersB,FranceMP,RimintonS, BrittonWJ.Secretedlymphotoxin-aisessentialforthecontrol ofanintracellularbacterialinfection.J Exp Med 2001;193: 239-46. 16. AllieN,KeetonR,CourtN,AbelB, FickL,VasseurV,et al. Limitedroleforlymphotoxinαinthehostimmuneresponse to Mycobacterium tuberculosis. J Immunol 2010;185:4292- 301. 17. JacobsM,TogbeD,FremondC,SamarinaA,AllieN,Botha T,et al.Tumornecrosisfactoriscriticaltocontroltuberculosis infection.Microbes Infect2007;5:623-8. 18. Price AL, Patterson J, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for straticationingenome-wideassociationstudies.Nat Genet 2006;38:904. 19. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, BlumenstielB, et al.Thestructureofhaplotypeblocksin thehumangenome.Science 2002;296:2225-9. 20. Indian GenomeVariation Consortium. Genetic landscapeof thepeopleofIndia:acanvasfordiseasegeneexploration.J Genet2008;87:3-20. 21. HillAV.Theimmunogeneticsofhumaninfectiousdiseases. Annu Rev Immunol 1998;16:593-617. 22. CardosoCC, PereiraAC, Brito-de-SouzaVN, Dias-Baptista IM, ManieroVC, VenturiniJ, et al.IFNG+874T>Asingle nucleotide polymorphism is associated with leprosy among Brazilians. Hum Genet2010;128:481-90. 23. KumarA, GhoshB.Asinglenucleotidepolymorphism(A > G)inintron3ofIFNgammageneisassociatedwithasthma. Genes Immun2008;9:294-301. 24. ChinoyH,SalwayF,JohnS,FertigN,TaitBD,OddisCV,et al.Interferon-gammaandInterleukin-4genepolymorphisms inUKCaucasianidiopathicinammatorymyopathypatients. Ann Rheum Dis2007;66:970-3. ABHIMANYUet al:CYTOKINEGENEVARIANTSINPTB 769 Reprint requests:Prof.MridulaBose,DepartmentofMicrobiology,VallabhbhaiPatelChestInstitute, UniversityofDelhi,Delhi110007,India e-mail:mridulabose@hotmail.com 25. BalsaA,DelAmoJ,BlancoF,CalizR,SilvaL,SanmartiR, et al. Prediction offunctional impairmentand remissionin rheumatoid arthritis patients by biochemical variables and geneticpolymorphisms.Rheumatology2010;49 :458-66. 26. AmirzargarAA, RezaeiN, JabbariH, DaneshAA, Khosravi F, Hajabdolbaghi M, et al. Cytokine single nucleotide polymorphisms in Iranian patients with pulmonary tuberculosis.Eur Cytokine Netw 2006;17 :84-9. 27. Möller M, Nebel A, Helden PDV, Schreiber S, Hoal EG. Analysis of eight genes modulating interferon gamma and human genetic susceptibility to tuberculosis: a case-control associationstudy.BMC Infect Dis 2010;10:154. 28. TindallEA,HayesVM,PetersenDC.Inammatorygenetic markersofprostatecancerrisk.Cancers 2010;2 :1198-220. 29. Abhimanyu, Mangangcha IR, Jha P, Arora K, Mukerji M, BanavalikerJN, et al.Differentialserumcytokinelevels are associated with cytokine gene polymorphisms in north Indianswithactivepulmonarytuberculosis.Infect Genet Evol 2011;11 :1015-22. 30. LeandroAC, RochaMA, CardosoCS, Bonecini-AlmeidaMG. GeneticpolymorphismsinvitaminDreceptor,vitaminD-binding protein, Toll-like receptor 2, nitric oxide synthase 2, and interferon-gammagenesanditsassociationwithsusceptibility totuberculosis.Braz J Med Biol Res 2009;42 :312-22. 770 INDIANJMEDRES,MAY 2012 . CI=0.1939 to 0.7445) of developing PTB in north Indians. Interpretation & conclusions: Our study showed six novel associations of cytokine gene variants with susceptibility to PTB in north Indians. . Footprints of genetic susceptibility to pulmonary tuberculosis: Cytokine gene variants in north Indians Abhimanyu,MridulaBose,PankajJha * &IndianGenomeVariationConsortium Department. different cellularreceptorsandhumoralfactors.Thesubsequent in ammatoryresponseisregulatedbytheproduction of pro-andanti -in ammatorycytokinesandchemokines. Interferon-gamma (IFN-γ one of the most important cytokines involved in