HLA class II gene associations in African American type 1 diabetes reveal a protective HLA-DRB1*03 haplotype.

AIMS: Owing to strong linkage disequilibrium between markers, pinpointing disease associations within genetic regions is difficult in European ancestral populations, most notably the very strong association of the HLA-DRB1*03-DQA1*05:01-DQB1*02:01 haplotype with Type 1 diabetes risk, which is assumed to be because of a combination of HLA-DRB1 and HLA-DQB1. In contrast, populations of African ancestry have greater haplotype diversity, offering the possibility of narrowing down regions and strengthening support for a particular gene in a region being causal. We aimed to study the human leukocyte antigen (HLA) region in African American Type 1 diabetes. METHODS: Two hundred and twenty-seven African American patients with Type 1 diabetes and 471 African American control subjects were tested for association at the HLA class II genes, HLA-DRB1, HLA-DQA1, HLA-DQB1 and 5147 single nucleotide polymorphisms across the major histocompatibility complex region using logistic regression models. Population admixture was accounted for with principal components analysis. RESULTS: Single nucleotide polymorphism marker associations were explained by the HLA associations, with the major peak over the class II loci. The HLA association overall was extremely strong, as expected for Type 1 diabetes, even in African Americans in whom diabetes diagnosis is heterogeneous. In addition, there were unique features: the HLA-DRB1*03 haplotype was split into HLA-DRB1*03:01, which confers greatest susceptibility in these samples (odds ratio 3.17, 95% CI 1.72-5.83) and HLA-DRB1*03:02, an allele rarely observed in Europeans, which confers the greatest protection in these African American samples (odds ratio 0.22, 95% CI 0.09-0.55). CONCLUSIONS: The unique diversity of the African HLA region we have uncovered supports a specific and major role for HLA-DRB1 in HLA-DRB1*03 haplotype-associated Type 1 diabetes risk.

Genetic association of zinc transporter 8 (ZnT8) autoantibodies in type 1 diabetes cases.

AIMS/HYPOTHESIS: Autoantibodies to zinc transporter 8 (ZnT8A) are associated with risk of type 1 diabetes. Apart from the SLC30A8 gene itself, little is known about the genetic basis of ZnT8A. We hypothesise that other loci in addition to SLC30A8 are associated with ZnT8A. METHODS: The levels of ZnT8A were measured in 2,239 British type 1 diabetic individuals diagnosed before age 17 years, with a median duration of diabetes of 4 years. Cases were tested at over 775,000 loci genome wide (including 53 type 1 diabetes associated regions) for association with positivity for ZnT8A. ZnT8A were also measured in an independent dataset of 855 family members with type 1 diabetes. RESULTS: Only FCRL3 on chromosome 1q23.1 and the HLA class I region were associated with positivity for ZnT8A. rs7522061T>C was the most associated single nucleotide polymorphism (SNP) in the FCRL3 region (p = 1.13 × 10(-16)). The association was confirmed in the family dataset (p ≤ 9.20 × 10(-4)). rs9258750A>G was the most associated variant in the HLA region (p = 2.06 × 10(-9) and p = 0.0014 in family cases). The presence of ZnT8A was not associated with HLA-DRB1, HLA-DQB1, HLA-A, HLA-B or HLA-C (p > 0.05). Unexpectedly, the two loci associated with the presence of ZnT8A did not alter risk of having type 1 diabetes, and the 53 type 1 diabetes risk loci did not influence positivity for ZnT8A, despite them being disease specific. CONCLUSIONS/INTERPRETATION: ZnT8A are not primary pathogenic factors in type 1 diabetes. Nevertheless, ZnT8A testing in combination with other autoantibodies facilitates disease prediction, despite the biomarker not being under the same genetic control as the disease.

Confirmation of novel type 1 diabetes risk loci in families.

AIMS/HYPOTHESIS: Over 50 regions of the genome have been associated with type 1 diabetes risk, mainly using large case/control collections. In a recent genome-wide association (GWA) study, 18 novel susceptibility loci were identified and replicated, including replication evidence from 2,319 families. Here, we, the Type 1 Diabetes Genetics Consortium (T1DGC), aimed to exclude the possibility that any of the 18 loci were false-positives due to population stratification by significantly increasing the statistical power of our family study. METHODS: We genotyped the most disease-predicting single-nucleotide polymorphisms at the 18 susceptibility loci in 3,108 families and used existing genotype data for 2,319 families from the original study, providing 7,013 parent-child trios for analysis. We tested for association using the transmission disequilibrium test. RESULTS: Seventeen of the 18 susceptibility loci reached nominal levels of significance (p < 0.05) in the expanded family collection, with 14q24.1 just falling short (p = 0.055). When we allowed for multiple testing, ten of the 17 nominally significant loci reached the required level of significance (p < 2.8 × 10(-3)). All susceptibility loci had consistent direction of effects with the original study. CONCLUSIONS/INTERPRETATION: The results for the novel GWA study-identified loci are genuine and not due to population stratification. The next step, namely correlation of the most disease-associated genotypes with phenotypes, such as RNA and protein expression analyses for the candidate genes within or near each of the susceptibility regions, can now proceed.

Analysis of 19 genes for association with type I diabetes in the Type I Diabetes Genetics Consortium families.

In recent years the pace of discovery of genetic associations with type I diabetes (T1D) has accelerated, with the total number of confirmed loci, including the major histocompatibility complex (MHC) region, reaching 43. However, much of the deciphering of the associations at these, and the established T1D loci, has yet to be performed in sufficient numbers of samples or with sufficient markers. Here, 257 single-nucleotide polymorphisms (SNPs) have been genotyped in 19 candidate genes (INS, PTPN22, IL2RA, CTLA4, IFIH1, SUMO4, VDR, PAX4, OAS1, IRS1, IL4, IL4R, IL13, IL12B, CEACAM21, CAPSL, Q7Z4c4(5Q), FOXP3, EFHB) in 2300 affected sib-pair families and tested for association with T1D as part of the Type I Diabetes Genetics Consortium's candidate gene study. The study had approximately 80% power at alpha=0.002 and a minor allele frequency of 0.2 to detect an effect with a relative risk (RR) of 1.20, which drops to just 40% power for a RR of 1.15. At the INS gene, rs689 (-23 HphI) was the most associated SNP (P=3.8 x 10(-31)), with the estimated RR=0.57 (95% confidence interval, 0.52-0.63). In addition, rs689 was associated with age-at-diagnosis of T1D (P=0.001), with homozygosity for the T1D protective T allele, delaying the onset of T1D by approximately 2 years in these families. At PTPN22, rs2476601 (R620W), in agreement with previous reports, was the most significantly associated SNP (P=6.9 x 10(-17)), with RR=1.55 (1.40-1.72). Evidence for association with T1D was observed for the IFIH1 SNP, rs1990760 (P=7.0 x 10(-4)), with RR=0.88 (0.82-0.95) and the CTLA4 SNP rs1427676 (P=0.0005), with RR=1.14 (1.06-1.23). In contrast, no convincing evidence of association was obtained for SUMO4, VDR, PAX4, OAS1, IRS1, IL4, IL4R, IL13, IL12B, CEACAM21 or CAPSL gene regions (http://www.T1DBase.org).

Analysis of 17 autoimmune disease-associated variants in type 1 diabetes identifies 6q23/TNFAIP3 as a susceptibility locus.

As a result of genome-wide association studies in larger sample sets, there has been an increase in identifying genes that influence susceptibility to individual immune-mediated diseases, as well as evidence that some genes are associated with more than one disease. In this study, we tested 17 single nucleotide polymorphisms (SNP) from 16 gene regions that have been reported in several autoimmune diseases including rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), multiple sclerosis (MS), ankylosing spondylitis (AS) and Crohn's disease (CD) to determine whether the variants are also associated with type 1 diabetes (T1D). In up to 8010 cases and 9733 controls we found some evidence for an association with T1D in the regions containing genes: 2q32/STAT4, 17q21/STAT3, 5p15/ERAP1 (ARTS1), 6q23/TNFAIP3 and 12q13/KIF5A/PIP4K2C with allelic P-values ranging from 3.70 x 10(-3) to 3.20 x 10(-5). These findings extend our knowledge of susceptibility locus sharing across different autoimmune diseases, and provide convincing evidence that the RA/SLE locus 6q23/TNFAIP3 is a newly identified T1D locus.

CD226 Gly307Ser association with multiple autoimmune diseases.

Genome-wide association studies provide insight into multigenic diseases through the identification of susceptibility genes and etiological pathways. In addition, the identification of shared variants among autoimmune disorders provides insight into common disease pathways. We previously reported an association of a nonsynonymous single nucleotide polymorphism (SNP) rs763361/Gly307Ser in the immune response gene CD226 on chromosome 18q22 with type 1 diabetes (T1D) susceptibility. Here, we report efforts toward identifying the causal variant by exonic resequencing and tag SNP mapping of the 18q22 region in both T1D and multiple sclerosis (MS). In addition to the analysis of newly available samples in T1D (2088 cases and 3289 controls) and autoimmune thyroid disease (AITD) (821 cases and 1920 controls), resulting in strong support for the Ser(307) association with T1D (P=3.46 x 10(-9)) and continued potential evidence for AITD (P=0.0345), we provide evidence for association of Gly307Ser with MS (P=4.20 x 10(-4)) and rheumatoid arthritis (RA) (P=0.017). The Ser(307) allele of rs763361 in exon 7 of CD226 predisposes to T1D, MS, and possibly AITD and RA, and based on the tag SNP analysis, could be the causal variant.

No association of multiple type 2 diabetes loci with type 1 diabetes.

AIMS/HYPOTHESIS: We used recently confirmed type 2 diabetes gene regions to investigate the genetic relationship between type 1 and type 2 diabetes, in an average of 7,606 type 1 diabetic individuals and 8,218 controls, providing >80% power to detect effects as small as an OR of 1.11 at a false-positive rate of 0.003. METHODS: The single nucleotide polymorphisms (SNPs) with the most convincing evidence of association in 12 type 2 diabetes-associated gene regions, PPARG, CDKAL1, HNF1B, WFS1, SLC30A8, CDKN2A-CDKN2B, IGF2BP2, KCNJ11, TCF7L2, FTO, HHEX-IDE and THADA, were analysed in type 1 diabetes cases and controls. PPARG and HHEX-IDE were additionally tested for association in 3,851 type 1 diabetes families. Tests for interaction with HLA class II genotypes, autoantibody status, sex, and age-at-diagnosis of type 1 diabetes were performed with all 12 gene regions. RESULTS: Only PPARG and HHEX-IDE showed any evidence of association with type 1 diabetes cases and controls (p = 0.004 and p = 0.003, respectively; p > 0.05 for other SNPs). The potential association of PPARG was supported by family analyses (p = 0.003; p (combined) = 1.0 x 10(-4)). No SNPs showed evidence of interaction with any covariate (p > 0.05). CONCLUSIONS/INTERPRETATION: We found no convincing genetic link between type 1 and type 2 diabetes. An association of PPARG (rs1801282/Pro12Ala) could be consistent with its known function in inflammation. Hence, our results reinforce evidence suggesting that type 1 diabetes is a disease of the immune system, rather than being due to inherited defects in beta cell function or regeneration or insulin resistance.

Confirmation of HLA class II independent type 1 diabetes associations in the major histocompatibility complex including HLA-B and HLA-A.

AIM: Until recently, human leucocyte antigen (HLA) class II-independent associations with type 1 diabetes (T1D) in the Major Histocompatibility Complex (MHC) region were not adequately characterized owing to insufficient map coverage, inadequate statistical approaches and strong linkage disequilibrium spanning the entire MHC. Here we test for HLA class II-independent associations in the MHC using fine mapping data generated by the Type 1 Diabetes Genetics Consortium (T1DGC). METHODS: We have applied recursive partitioning to the modelling of the class II loci and used stepwise conditional logistic regression to test approximately 1534 loci between 29 and 34 Mb on chromosome 6p21, typed in 2240 affected sibpair (ASP) families. RESULTS: Preliminary analyses confirm that HLA-B (at 31.4 Mb), HLA-A (at 30.0 Mb) are associated with T1D independently of the class II genes HLA-DRB1 and HLA-DQB1 (P = 6.0 x 10(-17) and 8.8 x 10(-13), respectively). In addition, a second class II region of association containing the single-nucleotide polymorphism (SNP), rs439121, and the class II locus HLA-DPB1, was identified as a T1D susceptibility effect which is independent of HLA-DRB1, HLA-DQB1 and HLA-B (P = 9.2 x 10(-8)). A younger age-at-diagnosis of T1D was found for HLA-B*39 (P = 7.6 x 10(-6)), and HLA-B*38 was protective for T1D. CONCLUSIONS: These analyses in the T1DGC families replicate our results obtained previously in approximately 2000 cases and controls and 850 families. Taking both studies together, there is evidence for four T1D-associated regions at 30.0 Mb (HLA-A), 31.4 Mb (HLA-B), 32.5 Mb (rs9268831/HLA-DRA) and 33.2 Mb (rs439121/HLA-DPB1) that are independent of HLA-DRB1/HLA-DQB1. Neither study found evidence of independent associations at HLA-C, HLA-DQA1 loci nor in the UBD/MAS1L or ITPR3 gene regions. These studies show that to find true class II-independent effects, large, well-powered sample collections are required and be genotyped with a dense map of markers. In addition, a robust statistical methodology that fully models the class II effects is necessary. Recursive partitioning is a useful tool for modelling these multiallelic systems.

Common allelic variants of exons 10, 12, and 33 of the thyroglobulin gene are not associated with autoimmune thyroid disease in the United Kingdom.

Thyroglobulin (Tg) is a major autoantigen for autoimmune thyroid disease (AITD). The Tg gene (Tg) has been mapped to chromosome 8q24, which has recently been linked in two independent studies to AITD. Association of specific alleles of microsatellite markers within Tg itself supports a role for Tg as a good candidate susceptibility locus for AITD. Resequencing of the Tg exons has led to the identification of a number of novel single nucleotide polymorphisms, four of which have been reported to be associated with AITD. Resequencing of Tg in Caucasian subjects in the United Kingdom (UK) has confirmed the presence of four single nucleotide polymorphisms in exons 10, 12, and 33. However, in the largest case-control association study to date with adequate power to detect the reported effect if present, we found no evidence for association of the Tg DNA variants with AITD in the UK. These data suggest that the recently identified single nucleotide polymorphisms do not have a causal role for AITD in the UK. At this stage, we cannot exclude the Tg region as harboring a susceptibility locus for AITD, and only large scale sequencing and fine mapping of the region, including neighboring genes, will allow us to identify any potential causal variants within this region.

A type 1 diabetes subgroup with a female bias is characterised by failure in tolerance to thyroid peroxidase at an early age and a strong association with the cytotoxic T-lymphocyte-associated antigen-4 gene.

AIMS/HYPOTHESIS: HLA haplotypes DRB1*03_DQB1*02 and DRB1*04_DQB1*0302, and allelic variation of the T cell regulatory gene cytotoxic T-lymphocyte-associated antigen-4 (CTLA4) and of the T cell activation gene protein tyrosine phosphatase, non-receptor type 22 (lymphoid) (PTPN22) have been associated with type 1 diabetes and autoimmune thyroid disease. Using thyroid peroxidase autoantibodies (TPOAbs) as an indicator of thyroid autoimmunity, we assessed whether the association of these loci is different in type 1 diabetes patients with TPOAbs than in those without. MATERIALS AND METHODS: TPOAbs were measured in 4,364 type 1 diabetic patients from across Great Britain, 67% of whom were aged under 18 years. These patients and 6,866 geographically matched control subjects were genotyped at CTLA4, PTPN22, HLA-DRB1 and HLA-DQB1. RESULTS: TPOAbs were detected in 462 (10.6%) of the type 1 diabetic patients. These patients had a stronger association with CTLA4 (odds ratio [OR] = 1.49 for the G allele of the single nucleotide polymorphism rs3087243; 95% CI = 1.29-1.72) than did the TPOAbs-negative patients (p = 0.0004; OR = 1.16; 95% CI = 1.10-1.24) or type 1 diabetes patients overall (OR = 1.20; 95% CI = 1.13-1.27). The ratio of women:men was higher (1.94:1) in this subgroup than in type 1 diabetes patients without TPOAbs (0.94:1; p = 1.86 x 10(-15)). TPOAbs status did not correlate with age at diagnosis of type 1 diabetes or with PTPN22 (Arg620Trp; rs2476601). CONCLUSIONS/INTERPRETATION: Our results identify a subgroup of type 1 diabetic patients that is sensitive to allelic variation of the negative regulatory molecule CTLA-4 and indicate that TPOAbs testing could be used to subclassify type 1 diabetes patients for inclusion in genetic, biological or clinical studies.

A systematic approach to the assessment of known TNF-alpha polymorphisms in Graves' disease.

Single-nucleotide polymorphisms (SNPs) within the tumour necrosis factor alpha (TNF-alpha) gene on chromosome 6p21.3 have been associated with many autoimmune diseases; however, results have been conflicting and accurate allele frequencies have never been established in a UK Caucasian population. The aim of this study was to assess the frequency of 22 known TNF-alpha SNPs in a UK Caucasian control population and investigate association of all polymorphisms with >5% minor allele frequency in a large case-control data set of patients with Graves' disease (GD). Eight of the 22 SNPs had minor allele frequencies >5% and were investigated further. The other 14 SNPs were present in the UK population at frequencies ranging from 0 to 4.7%. A significant increase of the A allele of the -238 SNP was seen in GD patients when compared with control subjects (9.6 vs 6.8%, respectively; P=0.003) and mirrored in the genotype distribution (P=0.009). Furthermore, association of the -238 SNP appears not to be due to linkage disequilibrium of the known HLA-DRB1(*)03 associations with GD. This study has established accurate allele frequencies of TNF-alpha SNPs in a UK population and provides preliminary evidence for association of the TNF-alpha gene with GD.

Genetic susceptibility to visceral leishmaniasis in The Sudan: linkage and association with IL4 and IFNGR1.

Longitudinal studies in Sudan show ethnic differences in incidence and clinical phenotypes associated with Leishmania donovani. Immunologically, bias in type 1 vs type 2 cytokine responses is important. To determine whether polymorphisms at IL4/IL9 or IFNGR1 contribute to susceptibility, we examined 59 multicase families of visceral leishmaniasis (VL) with/without post Kala-azar dermal leishmaniasis (PKDL). Multipoint nonparametric analysis (Allegro) linked IL4/IL9 to VL per se (P=0.002). Transmission disequilibrium testing with robust variance estimates confirmed association in the presence of linkage between VL per se and IL4 (P=0.008) but not IL9. Stepwise logistic regression analysis showed both IL4RP2 and IL4RP1 markers contributed significantly to the association, suggesting a common disease-associated haplotype. In contrast, IFNGR1 was linked (P=0.031) and associated (P=0.007) to PKDL but not VL or VL per se. Hence, polymorphism in a type 2 cytokine gene influences underlying susceptibility to VL, whereas IFNGR1 is specifically related to susceptibility to PKDL.

Evidence for a cluster of genes on chromosome 17q11-q21 controlling susceptibility to tuberculosis and leprosy in Brazilians.

The region of conserved synteny on mouse chromosome 11/human 17q11-q21 is known to carry a susceptibility gene(s) for intramacrophage pathogens. The region is rich in candidates including NOS2A, CCL2/MCP-1, CCL3/MIP-1alpha, CCL4/MIP-1beta, CCL5/RANTES, CCR7, STAT3 and STAT5A/5B. To examine the region in man, we studied 92 multicase tuberculosis (627 individuals) and 72 multicase leprosy (372 individuals) families from Brazil. Multipoint nonparametric analysis (ALLEGRO) using 16 microsatellites shows two peaks of linkage for leprosy at D17S250 (Z(lr) score 2.34; P=0.01) and D17S1795 (Z(lr) 2.67; P=0.004) and a single peak for tuberculosis at D17S250 (Z(lr) 2.04; P=0.02). Combined analysis shows significant linkage (peak Z(lr) 3.38) at D17S250, equivalent to an allele sharing LOD score 2.48 (P=0.0004). To determine whether one or multiple genes contribute, 49 informative single nucleotide polymorphisms were typed in candidate genes. Family-based allelic association testing that was robust to family clustering demonstrated significant associations with tuberculosis susceptibility at four loci separated by intervals (NOS2A-8.4 Mb-CCL18-32.3 kb-CCL4-6.04 Mb-STAT5B) up to several Mb. Stepwise conditional logistic regression analysis using a case/pseudo-control data set showed that the four genes contributed separate main effects, consistent with a cluster of susceptibility genes across 17q11.2.

Testing the possible negative association of type 1 diabetes and atopic disease by analysis of the interleukin 4 receptor gene.

Variations in the interleukin 4 receptor A (IL4RA) gene have been reported to be associated with atopy, asthma, and allergy, which may occur less frequently in subjects with type 1 diabetes (T1D). Since atopy shows a humoral immune reactivity pattern, and T1D results from a cellular (T lymphocyte) response, we hypothesised that alleles predisposing to atopy could be protective for T1D and transmitted less often than the expected 50% from heterozygous parents to offspring with T1D. We genotyped seven exonic single nucleotide polymorphisms (SNPs) and the -3223 C>T SNP in the putative promoter region of IL4RA in up to 3475 T1D families, including 1244 Finnish T1D families. Only the -3223 C>T SNP showed evidence of negative association (P=0.014). There was some evidence for an interaction between -3233 C>T and the T1D locus IDDM2 in the insulin gene region (P=0.001 in the combined and P=0.02 in the Finnish data set). We, therefore, cannot rule out a genetic effect of IL4RA in T1D, but it is not a major one.