Polymorphisms at the G72/G30 gene locus,on 13q33, are associated with bipolar disorder in two independent pedigree series

Linkage evidence suggests that chromosome 13 (13q32-33) contains susceptibility genes for both bipolar disorder and schizophrenia. Recently, genes called "G72" and "G30" were identified, and polymorphisms of these overlapping genes were reported to be associated with schizophrenia. We studied two series of pedigrees with bipolar disorder: the Clinical Neurogenetics (CNG) pedigrees (in which linkage to illness had been previously reported at 13q32-33), with 83 samples from 22 multiplex families, and the National Institute of Mental Health (NIMH) Genetics Initiative pedigrees, with 474 samples from 152 families. Sixteen single-nucleotide polymorphisms (SNPs) were genotyped at and around the G72/G30 locus, which covered a 157-kb region encompassing the entire complementary DNA sequences of G72 and G30. We performed transmission/disequilibrium testing (TDT) and haplotype analysis, since a linkage-disequilibrium block was present at this gene locus. In the CNG and NIMH data sets, the results of global TDT of the entire haplotype set were significant and consistent (P=.0004 and P=.008, respectively). In the CNG series, the associated genotypes divided the families into those with linkage and those without linkage (partitioned by the linkage evidence). Analysis of the decay of haplotype sharing gave a location estimate that included G72/G30 in its 95% confidence interval. Although statistically significant association was not detected for individual SNPs in the NIMH data set, the same haplotype was consistently overtransmitted in both series. These data suggest that a susceptibility variant for bipolar illness exists in the vicinity of the G72/G30 genes. Taken together with the earlier report, this is the first demonstration of a novel gene(s), discovered through a positional approach, independently associated with both bipolar illness and schizophrenia.     

Identification of sialyltransferase 8B as a generalized susceptibility gene for psychotic and mood disorders on chromosome 15q25-26

We previously identified a significant bipolar spectrum disorder linkage peak on 15q25-26 using 35 extended families with a broad clinical phenotype, including bipolar disorder (types I and II), recurrent unipolar depression and schizoaffective disorder. However, the specific gene(s) contributing to this signal had not been identified. By a fine mapping association study in an Australian case-control cohort (n?=?385), we find that the sialyltransferase 8B (ST8SIA2) gene, coding for an enzyme that glycosylates proteins involved in neuronal plasticity which has previously shown association to both schizophrenia and autism, is associated with increased risk to bipolar spectrum disorder. Nominal single point association was observed with SNPs in ST8SIA2 (rs4586379, P?=?0.0043; rs2168351, P?=?0.0045), and a specific risk haplotype was identified (frequency: bipolar vs controls?=?0.41 vs 0.31; χ(2)?=?6.46, P?=?0.011, OR?=?1.47). Over-representation of the specific risk haplotype was also observed in an Australian schizophrenia case-control cohort (n?=?256) (χ(2)?=?8.41, P?=?0.004, OR?=?1.82). Using GWAS data from the NIMH bipolar disorder (n?=?2055) and NIMH schizophrenia (n?=?2550) cohorts, the equivalent haplotype was significantly over-represented in bipolar disorder (χ(2)?=?5.91, P?=?0.015, OR?=?1.29), with the same direction of effect in schizophrenia, albeit non-significant (χ(2)?=?2.3, P?=?0.129, OR?=?1.09). We demonstrate marked down-regulation of ST8SIA2 gene expression across human brain development and show a significant haplotype×diagnosis effect on ST8SIA2 mRNA levels in adult cortex (ANOVA: F(1,87)?=?6.031, P?=?0.016). These findings suggest that variation the ST8SIA2 gene is associated with increased risk to mental illness, acting to restrict neuronal plasticity and disrupt early neuronal network formation, rendering the developing and adult brain more vulnerable to secondary genetic or environmental insults.     

A study of a genetic association between the PTGS2/PLA2G4A locus and schizophrenia

Six single nucleotide polymorphisms (SNPs) present in the PTGS2/PLA2G4A locus were detected among 118 British family trios of schizophrenia patients. The transmission disequilibrium test showed that SNP4 located in the 5'-flanking region of the PLA2G4A gene was associated with schizophrenia and that the haplotype analysis also showed a genetic association between the PTGS2 gene and schizophrenia. Because these two genes are arranged in a head-to-head configuration and separated by just about 149kb of DNA, they may have a combined effect on susceptibility to schizophrenia in some cases.     

Association of AKT1 gene variants and protein expression in both schizophrenia and bipolar disorder

The AKT1 gene has been associated with the genetic aetiology of schizophrenia. Following the overlap model of bipolar disorder and schizophrenia, we aimed to investigate AKT1 genetic variants and protein expression in both diseases. A total of 679 subjects with European ancestry were included: 384 with schizophrenia, 130 with bipolar disorder and 165 controls. Six single nucleotide polymorphisms (SNPs) were investigated for association with the diseases using single‐ and multi‐locus analyses. AKT1 and AKT2 protein levels were measured in post‐mortem brain tissues from ante‐mortem diagnosed schizophrenia (n = 30) and bipolar disorder subjects (n = 12) and matched controls. The analysis identified a significant global distortion in schizophrenia (P = 0.0026) and a weak association in bipolar disorder (P = 0.046). A sliding window procedure showed a five‐SNP haplotype (TCGAG) to be associated with schizophrenia (P = 1.22 × 10^?4) and bipolar disorder (P = 0.0041) and a four‐SNP haplotype (TCGA) with the combined sample (1.73 × 10^?5). On the basis of selected genotypes, a significant difference in protein expression emerged between subjects (P < 0.02). In conclusion, our findings, by showing the involvement of the AKT1 gene in both schizophrenia and bipolar disorder, support the role of AKT1 in the genetics of both disorders and add support to the view that there is some genetic overlap between them.     

The interaction between TPH2 promoter haplotypes and clinical-demographic risk factors in suicide victims with major psychoses

Tryptophan hydroxylase isoform 2 (TPH2) is a rate-limiting enzyme in the biosynthesis of serotonin (5-HT) and is predominantly localized in the brain. Previous studies have suggested that there is an association between serotonergic dysfunction in the brain and suicidality. This study was designed to examine whether the -473T > A and -8396G > C polymorphisms of the TPH2 gene may be associated with completed suicide in subjects with major psychoses from the Stanley Foundation Brain Bank sample. TPH2 genotypes were determined in 69 subjects with a diagnosis of schizophrenia or bipolar disorder, among which 22 died by suicide. Genomic DNA was amplified by polymerase chain reaction and typed by automated methods. Both markers were found to be in Hardy-Weinberg equilibrium and in strong linkage disequilibrium. No association with history of suicide was found for either polymorphism. Haplotype analysis with EHAP showed no association between completed suicide and haplotype distribution (chi2 = 1.877; 3 df; P = 0.598). Nor was there any association between suicide and these genetic markers even when clinical-demographic factors were considered as covariates in the haplotype analysis. These findings suggest that these 5' marker haplotypes in the TPH2 gene do not influence suicidal behaviour.     

Association study of eight circadian genes with bipolar I disorder, schizoaffective disorder and schizophrenia

We hypothesize that circadian dysfunction could underlie, at least partially, the liability for bipolar 1 disorder (BD1). Our hypothesis motivated tests for the association between the polymorphisms of genes that mediate circadian function and liability for BD1. The US Caucasian patients with BD1 (DSM-IV criteria) and available parents were recruited from Pittsburgh and surrounding areas (n = 138 cases, 196 parents) and also selected from the NIMH Genetics Collaborative Initiative (n = 96 cases, 192 parents). We assayed 44 informative single-nucleotide polymorphisms (SNPs) from eight circadian genes in the BD1 samples. A population-based sample, specifically cord blood samples from local live births, served as community-based controls (n = 180). It was used as a contrast for genotype and haplotype distributions with those of patients. US patients with schizophrenia/schizoaffective disorder (SZ/SZA, n = 331) and available parents from Pittsburgh (n = 344) were assayed for a smaller set of SNPs based on the results from the BD1 samples. Modest associations with SNPs at ARNTL (BmaL1) and TIMELESS genes were observed in the BD1 samples. The associations were detected using family-based and case-control analyses, albeit with different SNPs. Associations with TIMELESS and PERIOD3 were also detected in the Pittsburgh SZ/SZA group. Thus far, evidence for association between specific SNPs at the circadian gene loci and BD1 is tentative. Additional studies using larger samples are required to evaluate the associations reported here.     

Genomic View of Bipolar Disorder Revealed by Whole Genome Sequencing in a Genetic Isolate

Bipolar disorder is a common, heritable mental illness characterized by recurrent episodes of mania and depression. Despite considerable effort to elucidate the genetic underpinnings of bipolar disorder, causative genetic risk factors remain elusive. We conducted a comprehensive genomic analysis of bipolar disorder in a large Old Order Amish pedigree. Microsatellite genotypes and high-density SNP-array genotypes of 388 family members were combined with whole genome sequence data for 50 of these subjects, comprising 18 parent-child trios. This study design permitted evaluation of candidate variants within the context of haplotype structure by resolving the phase in sequenced parent-child trios and by imputation of variants into multiple unsequenced siblings. Non-parametric and parametric linkage analysis of the entire pedigree as well as on smaller clusters of families identified several nominally significant linkage peaks, each of which included dozens of predicted deleterious variants. Close inspection of exonic and regulatory variants in genes under the linkage peaks using family-based association tests revealed additional credible candidate genes for functional studies and further replication in population-based cohorts. However, despite the in-depth genomic characterization of this unique, large and multigenerational pedigree from a genetic isolate, there was no convergence of evidence implicating a particular set of risk loci or common pathways. The striking haplotype and locus heterogeneity we observed has profound implications for the design of studies of bipolar and other related disorders.     

Heat shock protein 70 gene polymorphisms are associated with paranoid schizophrenia in the Polish population

HSP70 genes have been considered as promising schizophrenia candidate genes based on their protective role in the central nervous system under stress conditions. In this study, we analyzed the potential implication of HSPA1A +190G/C, HSPA1B +1267A/G, and HSPA1L +2437T/C polymorphisms in the susceptibility to paranoid schizophrenia in a homogenous Caucasian Polish population. In addition, we investigated the association of the polymorphisms with the clinical variables of the disease. Two hundred and three patients with paranoid schizophrenia and 243 healthy controls were enrolled in the study. Polymorphisms of HSPA1A, -1B, and -1L genes were genotyped using the PCR-RFLP technique. Analyses were conducted in entire groups and in subgroups that were stratified according to gender. There were significant differences in the genotype and allele frequencies of HSPA1A polymorphism between the patients and controls. The +190CC genotype and +190C allele were over-represented in the patients and significantly increased the risk for developing schizophrenia (OR = 3.45 and OR = 1.61, respectively). Interestingly, such a risk was higher for females with the +190CC genotype than for males with the +190CC genotype (OR = 5.78 vs. OR = 2.76). We also identified the CGT haplotype as a risk haplotype for schizophrenia and demonstrated the effects of HSPA1A and HSPA1B genotypes on the psychopathology and age of onset. Our study provided the first evidence that the HSPA1A polymorphism may potentially increase the risk of developing paranoid schizophrenia. Further independent analyses in different populations to evaluate the role of gender are needed to replicate these results.     

Lack of association of EGR2 variants with bipolar disorder in Japanese population

The early growth response gene 2 (EGR2) has been recently reported to be associated with bipolar disorder in the Korean population. However replication studies in independent cohorts of same and different ethnicities are essential for establishing the credibility of a genotype–phenotype association. With this notion, in the present study we have performed a replication study of the reported association of SNPs in EGR2 in a case–control study comprising of 867 unrelated Japanese bipolar disorder patients and 895 age-, sex- and ethnicity-matched controls. Results showed no significant differences in allele and genotype frequencies of EGR2 SNPs between bipolar disorder patients and controls and also in a sex-stratified genetic analysis. The haplotype and meta-analyses also showed no significant association with bipolar disorder. In conclusion, this is the first replication study of the previously reported association of EGR2 with bipolar disorder in a larger sample set and the results showed that the EGR2 gene is unlikely to contribute to the susceptibility of bipolar disorder in a Japanese cohort.     

A Genetic Deconstruction of Neurocognitive Traits in Schizophrenia and Bipolar Disorder

Background

Impairments in cognitive functions are common in patients suffering from psychiatric disorders, such as schizophrenia and bipolar disorder. Cognitive traits have been proposed as useful for understanding the biological and genetic mechanisms implicated in cognitive function in healthy individuals and in the dysfunction observed in psychiatric disorders.

Methods

Sets of genes associated with a range of cognitive functions often impaired in schizophrenia and bipolar disorder were generated from a genome-wide association study (GWAS) on a sample comprising 670 healthy Norwegian adults who were phenotyped for a broad battery of cognitive tests. These gene sets were then tested for enrichment of association in GWASs of schizophrenia and bipolar disorder. The GWAS data was derived from three independent single-centre schizophrenia samples, three independent single-centre bipolar disorder samples, and the multi-centre schizophrenia and bipolar disorder samples from the Psychiatric Genomics Consortium.

Results

The strongest enrichments were observed for visuospatial attention and verbal abilities sets in bipolar disorder. Delayed verbal memory was also enriched in one sample of bipolar disorder. For schizophrenia, the strongest evidence of enrichment was observed for the sets of genes associated with performance in a colour-word interference test and for sets associated with memory learning slope.

Conclusions

Our results are consistent with the increasing evidence that cognitive functions share genetic factors with schizophrenia and bipolar disorder. Our data provides evidence that genetic studies using polygenic and pleiotropic models can be used to link specific cognitive functions with psychiatric disorders.     

SIRT1 gene, schizophrenia and bipolar disorder in the Japanese population: an association study

Several lines of evidence suggest that alterations in circadian rhythms might be associated with the pathophysiology of psychiatric disorders such as schizophrenia and bipolar disorder (BP). A recent study reported that SIRT1 is a molecule that plays an important role in the circadian clock system. Therefore, to evaluate the association among the SIRT1 gene, schizophrenia and BP, we conducted a case-control study of Japanese population samples (1158 schizophrenia patients, 1008 BP patients and 2127 controls) with four tagging SNPs (rs12778366, rs2273773, rs4746720 and rs10997875) in the SIRT1 gene. Marker-trait association analysis was used to evaluate the allele and the genotype association with the χ(2) test, and haplotype association analysis was evaluated with a likelihood ratio test. We showed an association between rs4746720 in the SIRT1 gene and schizophrenia in the allele and the genotype analysis. However, the significance of these associations did not survive after Bonferroni's correction for multiple testing. On the other hand, the SIRT1 gene was associated with Japanese schizophrenia in a haplotype-wise analysis (global P(all markers) = 4.89 × 10(-15)). Also, four tagging SNPs in the SIRT1 gene were not associated with BP. In conclusion, the SIRT1 gene may play an important role in the pathophysiology of schizophrenia in the Japanese population.     

Molecular mechanism of schizophrenia with reference to disrupted-in-schizophrenia 1 (DISC1)

Disrupted-in-schizophrenia 1 (DISC1) is a gene disrupted by a (1:1) (q42.1;q14.3) translocation that segregates with major psychiatric disorders in a Scottish family. To elucidate how DISC1 confers susceptibility to psychiatric disorders, identification of the molecules, which bind to the domain close to the translocation breakpoint in the DISC1 gene, was performed and fasciculation and elongation protein zeta-1 (Fez1), a novel DISC1-interacting protein, termed DISC1-binding zinc-finger protein (DBZ) and Kendrin were identified. The DISC1-Fez1 interaction is up-regulated by nerve growth factor (NGF) and involved in neurite extension. Transient dissociation of the DISC1-DBZ interaction by pituitary adenylate cyclase-activating polypeptide (PACAP) causes neurite extension. Furthermore, single-nucleotide polymorphisms association studies in a Japanese population have shown the relation of the Fez1, PACAP and PACAP receptor (PAC1) genes to schizophrenia. In schizophrenia with DISC1 translocation carrier, the DISC1-Fez1 and DISC1-DBZ interaction is disrupted, and it is likely that neural circuit formation remains immature, suggesting that schizophrenia is a neurodevelopmental disease. On the other hand, the DISC1-Kendrin interaction is suggested to be involved in microtubule network formation and an association between single-nucleotide polymorphisms of the Kendrin gene and bipolar disease has also been suggested in a Japanese population. This demonstrates that a part of bipolar disease is also a neurodevelopmental disorder.     

Bipolar I disorder and schizophrenia: a 440-single-nucleotide polymorphism screen of 64 candidate genes among Ashkenazi Jewish case-parent trios

Bipolar, schizophrenia, and schizoaffective disorders are common, highly heritable psychiatric disorders, for which familial coaggregation, as well as epidemiological and genetic evidence, suggests overlapping etiologies. No definitive susceptibility genes have yet been identified for any of these disorders. Genetic heterogeneity, combined with phenotypic imprecision and poor marker coverage, has contributed to the difficulty in defining risk variants. We focused on families of Ashkenazi Jewish descent, to reduce genetic heterogeneity, and, as a precursor to genomewide association studies, we undertook a single-nucleotide polymorphism (SNP) genotyping screen of 64 candidate genes (440 SNPs) chosen on the basis of previous linkage or of association and/or biological relevance. We genotyped an average of 6.9 SNPs per gene, with an average density of 1 SNP per 11.9 kb in 323 bipolar I disorder and 274 schizophrenia or schizoaffective Ashkenazi case-parent trios. Using single-SNP and haplotype-based transmission/disequilibrium tests, we ranked genes on the basis of strength of association (P<.01). Six genes (DAO, GRM3, GRM4, GRIN2B, IL2RB, and TUBA8) met this criterion for bipolar I disorder; only DAO has been previously associated with bipolar disorder. Six genes (RGS4, SCA1, GRM4, DPYSL2, NOS1, and GRID1) met this criterion for schizophrenia or schizoaffective disorder; five replicate previous associations, and one, GRID1, shows a novel association with schizophrenia. In addition, six genes (DPYSL2, DTNBP1, G30/G72, GRID1, GRM4, and NOS1) showed overlapping suggestive evidence of association in both disorders. These results may help to prioritize candidate genes for future study from among the many suspected/proposed for schizophrenia and bipolar disorders. They provide further support for shared genetic susceptibility between these two disorders that involve glutamate-signaling pathways.     

A Cytogenetic Abnormality and Rare Coding Variants Identify ABCA13 as a Candidate Gene in Schizophrenia, Bipolar Disorder, and Depression

Schizophrenia and bipolar disorder are leading causes of morbidity across all populations, with heritability estimates of ∼80% indicating a substantial genetic component. Population genetics and genome-wide association studies suggest an overlap of genetic risk factors between these illnesses but it is unclear how this genetic component is divided between common gene polymorphisms, rare genomic copy number variants, and rare gene sequence mutations. We report evidence that the lipid transporter gene ABCA13 is a susceptibility factor for both schizophrenia and bipolar disorder. After the initial discovery of its disruption by a chromosome abnormality in a person with schizophrenia, we resequenced ABCA13 exons in 100 cases with schizophrenia and 100 controls. Multiple rare coding variants were identified including one nonsense and nine missense mutations and compound heterozygosity/homozygosity in six cases. Variants were genotyped in additional schizophrenia, bipolar, depression (n > 1600), and control (n > 950) cohorts and the frequency of all rare variants combined was greater than controls in schizophrenia (OR = 1.93, p = 0.0057) and bipolar disorder (OR = 2.71, p = 0.00007). The population attributable risk of these mutations was 2.2% for schizophrenia and 4.0% for bipolar disorder. In a study of 21 families of mutation carriers, we genotyped affected and unaffected relatives and found significant linkage (LOD = 4.3) of rare variants with a phenotype including schizophrenia, bipolar disorder, and major depression. These data identify a candidate gene, highlight the genetic overlap between schizophrenia, bipolar disorder, and depression, and suggest that rare coding variants may contribute significantly to risk of these disorders.     

DISC1 gene polymorphisms and the risk of schizophrenia in an Iranian population: A preliminary study

Background: Schizophrenia, schizoaffective disorder, and bipolar illness are common psychological disorders with high heritability and variable phenotypes. The disrupted in schizophrenia 1 ( DISC1) gene, on chromosome 1q42, has an essential role in neurite outgrowth and cell signaling. The purpose of this study was to investigate the association of three single-nucleotide polymorphisms (SNPs; rs6675281, rs2255340, and rs2738864) with schizophrenia disorder. These three SNPs were chosen as they had been used in most of the previous studies. Methods: In a case-control study of Iranian population for the first time 778 blood samples were collected including, 402 schizophrenic patients and 376 healthy controls. Genomic DNA was extracted from peripheral blood using DNA extraction kit (BioFlux Co). The genotypes of rs6675281, rs2255340, and rs2738864 were detected by nested allele-specific multiplex polymersae chain reaction (PCR). Results: Our data revealed that the three SNPs are significantly associated with schizophrenia (rs2255349 C>T: confidence interval (CI), 2.115 to 3.268; P = 0.0000 OR: 2.629; rs2738864 C>T: CI, 1.538 to 2.339; P = 0.0000 OR: 1.897; rs6675281 C>T: CI, 2.788 to 4.662; P = 0.0009241 OR: 3.605). Through applying the expectation-maximization (EM) algorithm, we calculated the haplotype frequency, and finally performed haplotype analysis with Bonferroni correction and data preprocessing methods and the results showed rs66875281 to have the highest association. Discussion: Our findings primarily showed that DISC1 gene polymorphisms contribute to schizophrenia risk and have a significant association with this disorder among Iranian population. The strategy was found to be easy, rapid, specific, and consistent for the co-occurring detection of the DISC1 polymorphisms. We could finally confirm that the polymorphisms are related to schizophrenia studied in Iranian population.     

Association study between kynurenine 3-monooxygenase gene and schizophrenia in the Japanese population

Several lines of evidence suggest that metabolic changes in the kynurenic acid (KYNA) pathway are related to the etiology of schizophrenia. The inhibitor of kynurenine 3-monooxygenase (KMO) is known to increase KYNA levels, and the KMO gene is located in the chromosome region associated with schizophrenia, 1q42-q44. Single-marker and haplotype analyses for 6-tag single nucleotide polymorphisms (SNPs) of KMO were performed (cases = 465, controls = 440). Significant association of rs2275163 with schizophrenia was observed by single-marker comparisons (P = 0.032) and haplotype analysis including this SNP (P = 0.0049). Significant association of rs2275163 and haplotype was not replicated using a second, independent set of samples (cases = 480, controls = 448) (P = 0.706 and P = 0.689, respectively). These results suggest that the KMO is unlikely to be related to the development of schizophrenia in Japanese.     

Possible association of beta-arrestin 2 gene with methamphetamine use disorder, but not schizophrenia

Recent investigations suggest that the AKT/glycogen synthase kinase 3 (GSK3) signaling cascade may be associated with the pathophysiology of schizophrenia and methamphetamine (METH) use disorder. One important molecule related to this cascade is beta-arrestin 2 (ARRB2). We therefore conducted a genetic case-control association analysis of the gene for ARRB2 with schizophrenia and METH use disorder in a Japanese population (547 people with schizophrenia, 177 with METH use disorder and 546 controls). A possible association of 'tag single nucleotide polymorphisms (SNPs)' was found in METH use disorder (rs1045280: P(genotype) = 0.0118, P(allele) = 0.00351; rs2036657: P(allele) = 0.0431; rs4790694: P(genotype) = 0.0167, P(allele) = 0.0202), but no association was found with schizophrenia. We also evaluated the gene-gene interactions among ARRB2, AKT1, and GSK3B, which we previously reported for each of these diseases. However, no interaction was seen in our samples. This is the first association analysis of ARRB2, and our results indicate that ARRB2 may play a role in the pathophysiology of METH use disorder.     

Multiple SNPs in intron 41 of thyroglobulin gene are associated with autoimmune thyroid disease in the Japanese population

BACKGROUND: The etiology of the autoimmune thyroid diseases (AITDs), Graves' disease (GD) and Hashimoto's thyroiditis (HT), is largely unknown. However, genetic susceptibility is believed to play a major role. Two whole genome scans from Japan and from the US identified a locus on chromosome 8q24 that showed evidence for linkage with AITD and HT. Recent studies have demonstrated an association between thyroglobulin (Tg) polymorphisms and AITD in Caucasians, suggesting that Tg is a susceptibility gene on 8q24. OBJECTIVES: The objective of the study was to refine Tg association with AITD, by analyzing a panel of 25 SNPs across an extended 260 kb region of the Tg. METHODS: We studied 458 Japanese AITD patients (287 GD and 171 HT patients) and 221 matched Japanese control subjects in association studies. Case-control association studies were performed using 25 Tg single nucleotide polymorphisms (SNPs) chosen from a database of the Single Nucleotide Polymorphism Database (dbSNP). Haplotype analysis was undertaken using the computer program SNPAlyze version 7.0. PRINCIPAL FINDINGS AND CONCLUSIONS: In total, 5 SNPs revealed association with GD (P<0.05), with the strongest SNP associations at rs2256366 (P?=?0.002) and rs2687836 (P?=?0.0077), both located in intron 41 of the Tg gene. Because of the strong LD between these two strongest associated variants, we performed the haplotype analysis, and identified a major protective haplotype for GD (P?=?0.001). These results suggested that the Tg gene is involved in susceptibility for GD and AITD in the Japanese.