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
(). 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 (map.
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( />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 ( />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( />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( />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:
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