A replication study confirms the association of dendritic cell immunoreceptor (DCIR) polymorphisms with ACPA - negative RA in a large Asian cohort

Objectives

Dendritic cell immunoreceptor (DCIR) has been implicated in development of autoimmune disorders in rodent and DCIR polymorphisms were associated with anti-citrullinated proteins antibodies (ACPA)-negative rheumatoid arthritis (RA) in Swedish Caucasians. This study was undertaken to further investigate whether DCIR polymorphisms are also risk factors for the development of RA in four Asian populations originated from China and Malaysia.

Methods

We genotyped two DCIR SNPs rs2377422 and rs10840759 in Han Chinese population (1,193 cases, 1,278 controls), to assess their association with RA. Subsequently, rs2377422 was further genotyped in three independent cohorts of Malaysian-Chinese subjects (MY_Chinese, 254 cases, 206 controls), Malay subjects (MY_ Malay, 515 cases, 986 controls), and Malaysian-Indian subjects (MY_Indian, 378 cases, 285 controls), to seek confirmation of association in various ethnic groups. Meta-analysis was preformed to evaluate the contribution of rs2377422 polymorphisms to the development of ACPA-negative RA in distinct ethnic groups. Finally, we carried out association analysis of rs2377422 polymorphisms with DCIR mRNA expression levels.

Results

DCIR rs2377422 was found to be significantly associated with ACPA -negative RA in Han Chinese (OR 1.92, 95% CI 1.27–2.90, P = 0.0020). Meta-analysis confirms DCIR rs2377422 as a risk factor for ACPA-negative RA across distinct ethnic groups (OR_overall = 1.17, 95% CI 1.06–1.30, P = 0.003). The SNP rs2377422 polymorphism showed significant association with DCIR mRNA expression level, i.e. RA-risk CC genotype exhibit a significant increase in the expression of DCIR (P = 0.0023, Kruskal–Wallis).

Conclusions

Our data provide evidence for association between DCIR rs2377422 and RA in non-Caucasian populations and confirm the influence of DCIR polymorphisms on RA susceptibility, especially on ACPA-negative RA.     

Replication of putative susceptibility loci from genome-wide association studies associated with coronary atherosclerosis in Chinese Han population

Background

Coronary atherosclerosis, the main cause of cardiovascular disease, is a progressive disease. Recent Genome Wide Association Studies (GWASs) discovered several novel loci associated with coronary artery disease (CAD) or its main complication myocardial infarction (MI). In this study, we investigated the associations between previously reported CAD- and MI-associated variants and coronary atherosclerosis in Chinese Han population.

Methodology/Principal Findings

We performed a case-control association study with 2,335 coronary atherosclerosis patients and 1,078 controls undergoing coronary angiography of Chinese Han from China. Fourteen single nucleotide polymorphisms (SNPs), located at 1p13.3, 1q41, 2q36.3, 6q25.1, 9p21.3, 10q11.21 and 15q22.33, were genotyped in our sample collection. Six SNPs at 9p21 were associated with coronary atherosclerosis susceptibility (P_trend<0.05) and rs10757274 showed the most significant association (P = 2.38×10⁻⁰⁸, OR = 1.34). These associations remained significant after adjustment for multiple comparisons. Rs17465637 at 1q41 (P_trend = 6.83×10⁻⁰³, OR = 0.86) also showed significant association with coronary atherosclerosis, but the association was not significant after multiple comparisons. Additionally, rs501120 (P = 8.36×10⁻⁰³, OR = 0.80) at 10q11.21 was associated with coronary atherosclerosis in females, but did not show association in males and all participants. Variants at 1p13.3, 2q36.3, 6q25.1 and 15q22.33 showed no associations with coronary atherosclerosis and main cardiovascular risk factors in our data.

Conclusions/Significance

Our findings indicated variants at 9p21 were significantly associated with coronary atherosclerosis in Han Chinese. Variants at 1q41 showed suggestive evidence of association and variants at 10q11.21 showed suggestive evidence of association in females, which warrant further study in a larger sample.     

Association of new loci identified in European genome-wide association studies with susceptibility to type 2 diabetes in the Japanese

Background

Several novel susceptibility loci for type 2 diabetes have been identified through genome-wide association studies (GWAS) for type 2 diabetes or quantitative traits related to glucose metabolism in European populations. To investigate the association of the 13 new European GWAS-derived susceptibility loci with type 2 diabetes in the Japanese population, we conducted a replication study using 3 independent Japanese case-control studies.

Methodology/Principal Findings

We examined the association of single nucleotide polymorphisms (SNPs) within 13 loci (MTNR1B, GCK, IRS1, PROX1, BCL11A, ZBED3, KLF14, TP53INP1, KCNQ1, CENTD2, HMGA2, ZFAND6 and PRC1) with type 2 diabetes using 4,964 participants (2,839 cases and 2,125 controls) from 3 independent Japanese samples. The association of each SNP with type 2 diabetes was analyzed by logistic regression analysis. Further, we performed combined meta-analyses for the 3 studies and previously performed Japanese GWAS data (4,470 cases vs. 3,071 controls). The meta-analysis revealed that rs2943641 in the IRS1 locus was significantly associated with type 2 diabetes, (P = 0.0034, OR = 1.15 95% confidence interval; 1.05–1.26) and 3 SNPs, rs10930963 in the MTNR1B locus, rs972283 in the KLF14 locus, and rs231362 in the KCNQ1 locus, had nominal association with type 2 diabetes in the present Japanese samples (P<0.05).

Conclusions

These results indicate that IRS1 locus may be common locus for type 2 diabetes across different ethnicities.     

Identification of six novel susceptibility loci for dyslipidemia using longitudinal exome-wide association studies in a Japanese population

Recent genome-wide association studies have identified various dyslipidemia-related genetic variants. However, most studies were conducted in a cross-sectional manner. We thus performed longitudinal exome-wide association studies of dyslipidemia in a Japanese population. We used ~244,000 genetic variants and clinical data of 6022 Japanese individuals who had undergone annual health checkups for several years. After quality control, the association of dyslipidemia-related phenotypes with 24,691 single nucleotide polymorphisms (SNPs) was tested using the generalized estimating equation model. In total, 82 SNPs were significantly (P < 2.03 × 10^?6) associated with dyslipidemia phenotypes. Of these SNPs, four (rs74416240 of TCHP, rs925368 of GIT2, rs7969300 of ATXN2, and rs12231744 of NAA25) and two (rs34902660 of SLC17A3 and rs1042127 of CDSN) were identified as novel genetic determinants of hypo-HDL- and hyper-LDL-cholesterolemia, respectively. A replication study using the cross-sectional data of 8310 Japanese individuals showed the association of the six identified SNPs with dyslipidemia-related traits.     

Identification of novel susceptibility Loci for kawasaki disease in a Han chinese population by a genome-wide association study

Kawasaki disease (KD) is an acute systemic vasculitis syndrome that primarily affects infants and young children. Its etiology is unknown; however, epidemiological findings suggest that genetic predisposition underlies disease susceptibility. Taiwan has the third-highest incidence of KD in the world, after Japan and Korea. To investigate novel mechanisms that might predispose individuals to KD, we conducted a genome-wide association study (GWAS) in 250 KD patients and 446 controls in a Han Chinese population residing in Taiwan, and further validated our findings in an independent Han Chinese cohort of 208 cases and 366 controls. The most strongly associated single-nucleotide polymorphisms (SNPs) detected in the joint analysis corresponded to three novel loci. Among these KD-associated SNPs three were close to the COPB2 (coatomer protein complex beta-2 subunit) gene: rs1873668 (p = 9.52×10⁻⁵), rs4243399 (p = 9.93×10⁻⁵), and rs16849083 (p = 9.93×10⁻⁵). We also identified a SNP in the intronic region of the ERAP1 (endoplasmic reticulum amino peptidase 1) gene (rs149481, p_best = 4.61×10⁻⁵). Six SNPs (rs17113284, rs8005468, rs10129255, rs2007467, rs10150241, and rs12590667) clustered in an area containing immunoglobulin heavy chain variable regions genes, with p_best-values between 2.08×10⁻⁵ and 8.93×10⁻⁶, were also identified. This is the first KD GWAS performed in a Han Chinese population. The novel KD candidates we identified have been implicated in T cell receptor signaling, regulation of proinflammatory cytokines, as well as antibody-mediated immune responses. These findings may lead to a better understanding of the underlying molecular pathogenesis of KD.     

Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10⁻⁶⁴) and lower levels of eicosapentaenoic acid (EPA, p = 5×10⁻⁵⁸) and docosapentaenoic acid (DPA, p = 4×10⁻¹⁵⁴). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10⁻¹²) and DPA (p = 1×10⁻⁴³) and lower docosahexaenoic acid (DHA, p = 1×10⁻¹⁵). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10⁻⁸). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.     

Identification of genetic loci associated with growth traits at weaning in yak through a genome-wide association study

A multilocus GWAS was performed to explore the genetic architecture of four growth traits in yak. In total, 354 female yaks for which measurements of body weight (BW), withers height (WH), body length (BL) and chest girth (CG) at weaning were available underwent genotyping with the Illumina BovineHD BeadChip (770K). After quality control, we retained 98 688 SNPs and 354 animals for GWAS analysis. We identified seven, 18, seven and nine SNPs (corresponding to seven, 17, seven and eight candidate genes) associated with BW, WH, BL and CG at weaning respectively. Interestingly, most of these candidate genes were reported to be involved in growth-related processes such as muscle formation, lipid deposition, feed efficiency, carcass composition and development of the central and peripheral nervous system. Our results offer novel insight into the molecular architecture underpinning yak growth traits. Further functional analyses are thus warranted to explore the molecular mechanisms whereby these genes affect these traits of interest.     

On the meta-analysis of genome-wide association studies: a robust and efficient approach to combine population and family-based studies

For the meta-analysis of genome-wide association studies, we propose a new method to adjust for the population stratification and a linear mixed approach that combines family-based and unrelated samples. The proposed approach achieves similar power levels as a standard meta-analysis which combines the different test statistics or p values across studies. However, by virtue of its design, the proposed approach is robust against population admixture and stratification, and no adjustments for population admixture and stratification, even in unrelated samples, are required. Using simulation studies, we examine the power of the proposed method and compare it to standard approaches in the meta-analysis of genome-wide association studies. The practical features of the approach are illustrated with a meta-analysis of three genome-wide association studies for Alzheimer's disease. We identify three single nucleotide polymorphisms showing significant genome-wide association with affection status. Two single nucleotide polymorphisms are novel and will be verified in other populations in our follow-up study.     

Anticitrullinated protein antibody (ACPA) in rheumatoid arthritis: influence of an interaction between HLA-DRB1 shared epitope and a deletion polymorphism in glutathione s-transferase in a cross-sectional study

Introduction

A deletion polymorphism in glutathione S-transferase Mu-1 (GSTM1-null) has previously been implicated to play a role in rheumatoid arthritis (RA) risk and progression, although no prior investigations have examined its associations with anticitrullinated protein antibody (ACPA) positivity. The purpose of this study was to examine the associations of GSTM1-null with ACPA positivity in RA and to assess for evidence of interaction between GSTM1 and HLA-DRB1 shared epitope (SE).

Methods

Associations of GSTM1-null with ACPA positivity were examined separately in two RA cohorts, the Veterans Affairs Rheumatoid Arthritis (VARA) registry (n = 703) and the Study of New-Onset RA (SONORA; n = 610). Interactions were examined by calculating an attributable proportion (AP) due to interaction.

Results

A majority of patients in the VARA registry (76%) and SONORA (69%) were positive for ACPA with a similar frequency of GSTM1-null (53% and 52%, respectively) and HLA-DRB1 SE positivity (76% and 71%, respectively). The parameter of patients who had ever smoked was more common in the VARA registry (80%) than in SONORA (65%). GSTM1-null was significantly associated with ACPA positivity in the VARA registry (odds ratio (OR), 1.45; 95% confidence interval (CI), 1.02 to 2.05), but not in SONORA (OR, 1.00; 95% CI, 0.71 to 1.42). There were significant additive interactions between GSTM1 and HLA-DRB1 SE in the VARA registry (AP, 0.49; 95% CI, 0.21 to 0.77; P < 0.001) in ACPA positivity, an interaction replicated in SONORA (AP, 0.38; 95% CI, 0.00 to 0.76; P = 0.050).

Conclusions

This study is the first to show that the GSTM1-null genotype, a common genetic variant, exerts significant additive interaction with HLA-DRB1 SE on the risk of ACPA positivity in RA. Since GSTM1 has known antioxidant functions, these data suggest that oxidative stress may be important in the development of RA-specific autoimmunity in genetically susceptible individuals.     

Genetic determinants of cardiovascular events among women with migraine: a genome-wide association study

Background

Migraine is associated with an increased risk for cardiovascular disease (CVD). Both migraine and CVD are highly heritable. However, the genetic liability for CVD among migraineurs is unclear.

Methods

We performed a genome-wide association study for incident CVD events during 12 years of follow-up among 5,122 migraineurs participating in the population-based Women''s Genome Health Study. Migraine was self-reported and CVD events were confirmed after medical records review. We calculated odds ratios (OR) and 95% confidence intervals (CI) and considered a genome-wide p-value <5×10⁻⁸ as significant.

Results

Among the 5,122 women with migraine 164 incident CVD events occurred during follow-up. No SNP was associated with major CVD, ischemic stroke, myocardial infarction, or CVD death at the genome-wide level; however, five SNPs showed association with p<5×10⁻⁶. Among migraineurs with aura rs7698623 in MEPE (OR = 6.37; 95% CI 3.15–12.90; p = 2.7×10⁻⁷) and rs4975709 in IRX4 (OR = 5.06; 95% CI 2.66–9.62; p = 7.7×10⁻⁷) appeared to be associated with ischemic stroke, rs2143678 located close to MDF1 with major CVD (OR = 3.05; 95% CI 1.98–4.69; p = 4.3×10⁻⁷), and the intergenic rs1406961 with CVD death (OR = 12.33; 95% CI 4.62–32.87; p = 5.2×10⁻⁷). Further, rs1047964 in BACE1 appeared to be associated with CVD death among women with any migraine (OR = 4.67; 95% CI 2.53–8.62; p = 8.0×10⁻⁷).

Conclusion

Our results provide some suggestion for an association of five SNPs with CVD events among women with migraine; none of the results was genome-wide significant. Four associations appeared among migraineurs with aura, two of those with ischemic stroke. Although our population is among the largest with migraine and incident CVD information, these results must be treated with caution, given the limited number of CVD events among women with migraine and the low minor allele frequencies for three of the SNPs. Our results await independent replication and should be considered hypothesis generating for future research.     

Mapping of quantitative trait loci for carcass traits in a Japanese Black (Wagyu) cattle population

To detect quantitative trait loci (QTL) that influence economically important traits in a purebred Japanese Black cattle population, we performed a preliminary genome-wide scan using 187 microsatellite markers across a paternal half-sib family composed of 258 offspring. We located six QTL at the 1% chromosome-wise level on bovine chromosomes (BTA) 4, 6, 13, 14 and 21. A second screen of these six QTL regions using 138 additional paternal offspring half-sib from the same sire, provided further support for five QTL: carcass weight on BTA14 (22-39 cM), one for rib thickness on BTA6 (27-58 cM) and three for beef marbling score (BMS) on BTA4 (59-67 cM), BTA6 (68-89 cM) and BTA21 (75-84 cM). The location of QTL for subcutaneous fat thickness on BTA13 was not supported by the second screen (P > 0.05). We determined that the combined contribution of the three QTLs for BMS was 10.1% of the total variance. The combined phenotypic average of these three Q was significantly different (P < 0.001) from those of other allele combinations. Analysis of additional half-sib families will be necessary to confirm these QTL.     

Association of TNFAIP3 interacting protein 1, TNIP1 with systemic lupus erythematosus in a Japanese population: a case-control association study

Introduction

TNFAIP3 interacting protein 1, TNIP1 (ABIN-1) is involved in inhibition of nuclear factor-κB (NF-κB) activation by interacting with TNF alpha-induced protein 3, A20 (TNFAIP3), an established susceptibility gene to systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). Recent genome-wide association studies revealed association of TNIP1 with SLE in the Caucasian and Chinese populations. In this study, we investigated whether the association of TNIP1 with SLE was replicated in a Japanese population. In addition, association of TNIP1 with RA was also examined.

Methods

A case-control association study was conducted on the TNIP1 single nucleotide polymorphism (SNP) rs7708392 in 364 Japanese SLE patients, 553 RA patients and 513 healthy controls.

Results

Association of TNIP1 rs7708392C was replicated in Japanese SLE (allele frequency in SLE: 76.5%, control: 69.9%, P = 0.0022, odds ratio [OR] 1.40, 95% confidence interval [CI] 1.13-1.74). Notably, the risk allele frequency in the healthy controls was considerably greater in Japanese (69.9%) than in Caucasians (24.3%). A tendency of stronger association was observed in the SLE patients with renal disorder (P = 0.00065, OR 1.60 [95%CI 1.22-2.10]) than in all SLE patients (P = 0.0022, OR 1.40 [95%CI 1.13-1.74]). Significant association with RA was not observed, regardless of the carriage of human leukocyte antigen DR β1 (HLA-DRB1) shared epitope. Significant gene-gene interaction between TNIP1 and TNFAIP3 was detected neither in SLE nor RA.

Conclusions

Association of TNIP1 with SLE was confirmed in a Japanese population. TNIP1 is a shared SLE susceptibility gene in the Caucasian and Asian populations, but the genetic contribution appeared to be greater in the Japanese and Chinese populations because of the higher risk allele frequency. Taken together with the association of TNFAIP3, these observations underscore the crucial role of NF-κB regulation in the pathogenesis of SLE.     

The association of rs1149048 polymorphism in Matrilin-1(MATN1) gene with adolescent idiopathic scoliosis susceptibility: a meta-analysis

Several previous studies have evaluated the association between rs1149048 polymorphism in the matrilin-1 gene (MATN1) and the risk of adolescent idiopathic scoliosis (AIS). However the results of those studies were inconsistent. We conducted this meta-analysis to assess whether rs1149048 polymorphism was involved in the risk of AIS and evaluated the associations in different ethnicities. Electronic databases, such as: PubMed, EMBASE, WANFANG databases in any languages up to Dec 2012 were searched to assess the association between rs1149048 polymorphism and AIS. Meta-analysis was performed by STATA 12.0 software to estimate the pooled odds ratio (OR) and the 95 % confidence interval (CI). Finally four papers including five studies which involved 1436 AIS patients and 1,879 controls were identified for this meta-analysis. The results showed that G allele of the rs1149048 was significantly associated with increased AIS risk [OR = 1.13, 95 % CI (1.02–1.25), P = 0.023]. As for genotype (GG vs. GA + AA), homozygous GG genotype was also found to be a risk factor of developing AIS. The subgroup meta-analysis results showed G allele and GG genotype were significantly associated with AIS in Asian group but not in Caucasian group. Neither Egger’s test nor Begg’s test found evidence of publication bias in current study (P > 0.05). In summary, this meta-analysis found an overall significant association of rs1149048 polymorphism with risk of AIS, especially in Asian population. The relationship between rs1149048 polymorphism and AIS in other ethnic population is needed to be investigated.     

Lack of association between catalase gene polymorphism (T/C exon 9) and susceptibility to vitiligo in a Turkish population

Accumulation of hydrogen peroxide (H(2)O(2)) and low catalase (CAT) activity have been demonstrated in the epidermis of vitiligo patients. We investigated a possible association between the CAT exon 9 (Asp-389) gene and vitiligo susceptibility in the Turkish population. Thirty-four patients with vitiligo and 49 gender, age and ethnic matched controls were enrolled in the study. Genotyping was done by PCR-RFLP. The CAT exon 9 (Asp-389) genotype and allele frequencies of vitiligo patients did not differ significantly from those of healthy controls. We found no association between CAT (Asp-389) gene polymorphism and vitiligo susceptibility in Turkish vitiligo patients.     

Novel association of ABO histo-blood group antigen with soluble ICAM-1: results of a genome-wide association study of 6,578 women

While circulating levels of soluble Intercellular Adhesion Molecule 1 (sICAM-1) have been associated with diverse conditions including myocardial infarction, stroke, malaria, and diabetes, comprehensive analysis of the common genetic determinants of sICAM-1 is not available. In a genome-wide association study conducted among 6,578 participants in the Women's Genome Health Study, we find that three SNPs at the ICAM1 (19p13.2) locus (rs1799969, rs5498 and rs281437) are non-redundantly associated with plasma sICAM-1 concentrations at a genome-wide significance level (P<5×10⁻⁸), thus extending prior results from linkage and candidate gene studies. We also find that a single SNP (rs507666, P=5.1×10⁻²⁹) at the ABO (9q34.2) locus is highly correlated with sICAM-1 concentrations. The novel association at the ABO locus provides evidence for a previously unknown regulatory role of histo-blood group antigens in inflammatory adhesion processes.     

A genome-wide association study of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality

Background

Numerous quantitative trait loci (QTL) have been detected in pigs over the past 20 years using microsatellite markers. However, due to the low density of these markers, the accuracy of QTL location has generally been poor. Since 2009, the dense genome coverage provided by the Illumina PorcineSNP60 BeadChip has made it possible to more accurately map QTL using genome-wide association studies (GWAS). Our objective was to perform high-density GWAS in order to identify genomic regions and corresponding haplotypes associated with production traits in a French Large White population of pigs.

Methods

Animals (385 Large White pigs from 106 sires) were genotyped using the PorcineSNP60 BeadChip and evaluated for 19 traits related to feed intake, growth, carcass composition and meat quality. Of the 64 432 SNPs on the chip, 44 412 were used for GWAS with an animal mixed model that included a regression coefficient for the tested SNPs and a genomic kinship matrix. SNP haplotype effects in QTL regions were then tested for association with phenotypes following phase reconstruction based on the Sscrofa10.2 pig genome assembly.

Results

Twenty-three QTL regions were identified on autosomes and their effects ranged from 0.25 to 0.75 phenotypic standard deviation units for feed intake and feed efficiency (four QTL), carcass (12 QTL) and meat quality traits (seven QTL). The 10 most significant QTL regions had effects on carcass (chromosomes 7, 10, 16, 17 and 18) and meat quality traits (two regions on chromosome 1 and one region on chromosomes 8, 9 and 13). Thirteen of the 23 QTL regions had not been previously described. A haplotype block of 183 kb on chromosome 1 (six SNPs) was identified and displayed three distinct haplotypes with significant (0.0001 < P < 0.03) associations with all evaluated meat quality traits.

Conclusions

GWAS analyses with the PorcineSNP60 BeadChip enabled the detection of 23 QTL regions that affect feed consumption, carcass and meat quality traits in a LW population, of which 13 were novel QTL. The proportionally larger number of QTL found for meat quality traits suggests a specific opportunity for improving these traits in the pig by genomic selection.     

IFNG +874T/A, IL10 -1082G/A and TNF -308G/A polymorphisms in association with tuberculosis susceptibility: a meta-analysis study

Susceptibility to infectious diseases is influenced by genetic background and efficient cellular immune activation is responsible for protection. In tuberculosis (TB), interferon-gamma (IFNγ) is crucial to control intracellular growth of Mycobacterium tuberculosis while interleukin-10 (IL-10) has an antagonistic role. Tumor necrosis factor (TNF) is a central mediator of granuloma formation and control of bacilli spread synergizing with IFNγ to hamper M. tuberculosis infection. Single nucleotide polymorphisms (SNPs) located at these genes could influence cytokine levels and regulate resistance and susceptibility to TB. The aim of this study was to determine the association of the interferon-gamma gene (IFNG) +874T/A, interleukin-10 gene (IL10) -1082G/A and tumor necrosis factor gene (TNF) -308G/A SNPs with TB in several populations using meta-analysis. We searched for association studies correlating these polymorphisms and TB using pre-established keywords in Medline. Meta-analysis was conducted with random effects models to account for heterogeneity between studies. Eleven studies were included in the IFNG +874T/A meta-analysis, while eight were used for the IL10 -1082G/A, and 10 were employed for TNF -308G/A. Data were analyzed in respect to associations between alleles, genotypes and minor allele carriers. Statistically significant results were found only for IFNG. The +874T allele of IFNG showed a protective significant association (OR = 0.75; 95% CI, 0.634–0.887; P = 0.0008). Though not significant, IL10 presented a trend towards protection when only studies with pulmonary TB patients were considered. This data reinforces the critical importance of IFNG +874T/A as a genetic marker for TB resistance and this information can be used for better design of a TB vaccine. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.