Sequencing and association analysis of the type 1 diabetes-linked region on chromosome 10p12-q11.

BACKGROUND: In an effort to locate susceptibility genes for type 1 diabetes (T1D) several genome-wide linkage scans have been undertaken. A chromosomal region designated IDDM10 retained genome-wide significance in a combined analysis of the main linkage scans. Here, we studied sequence polymorphisms in 23 Mb on chromosome 10p12-q11, including the putative IDDM10 region, to identify genes associated with T1D. RESULTS: Initially, we resequenced the functional candidate genes, CREM and SDF1, located in this region, genotyped 13 tag single nucleotide polymorphisms (SNPs) and found no association with T1D. We then undertook analysis of the whole 23 Mb region. We constructed and sequenced a contig tile path from two bacterial artificial clone libraries. By comparison with a clone library from an unrelated person used in the Human Genome Project, we identified 12,058 SNPs. We genotyped 303 SNPs and 25 polymorphic microsatellite markers in 765 multiplex T1D families and followed up 22 associated polymorphisms in up to 2,857 families. We found nominal evidence of association in six loci (P = 0.05 - 0.0026), located near the PAPD1 gene. Therefore, we resequenced 38.8 kb in this region, found 147 SNPs and genotyped 84 of them in the T1D families. We also tested 13 polymorphisms in the PAPD1 gene and in five other loci in 1,612 T1D patients and 1,828 controls from the UK. Overall, only the D10S193 microsatellite marker located 28 kb downstream of PAPD1 showed nominal evidence of association in both T1D families and in the case-control sample (P = 0.037 and 0.03, respectively). CONCLUSION: We conclude that polymorphisms in the CREM and SDF1 genes have no major effect on T1D. The weak T1D association that we detected in the association scan near the PAPD1 gene may be either false or due to a small genuine effect, and cannot explain linkage at the IDDM10 region.

Analysis of polymorphisms in 16 genes in type 1 diabetes that have been associated with other immune-mediated diseases.

BACKGROUND: The identification of the HLA class II, insulin (INS), CTLA-4 and PTPN22 genes as determinants of type 1 diabetes (T1D) susceptibility indicates that fine tuning of the immune system is centrally involved in disease development. Some genes have been shown to affect several immune-mediated diseases. Therefore, we tested the hypothesis that alleles of susceptibility genes previously associated with other immune-mediated diseases might perturb immune homeostasis, and hence also associate with predisposition to T1D. METHODS: We resequenced and genotyped tag single nucleotide polymorphisms (SNPs) from two genes, CRP and FCER1B, and genotyped 27 disease-associated polymorphisms from thirteen gene regions, namely FCRL3, CFH, SLC9A3R1, PADI4, RUNX1, SPINK5, IL1RN, IL1RA, CARD15, IBD5-locus (including SLC22A4), LAG3, ADAM33 and NFKB1. These genes have been associated previously with susceptibility to a range of immune-mediated diseases including rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), Graves' disease (GD), psoriasis, psoriatic arthritis (PA), atopy, asthma, Crohn disease and multiple sclerosis (MS). Our T1D collections are divided into three sample subsets, consisting of set 1 families (up to 754 families), set 2 families (up to 743 families), and a case-control collection (ranging from 1,500 to 4,400 cases and 1,500 to 4,600 controls). Each SNP was genotyped in one or more of these subsets. Our study typically had approximately 80% statistical power for a minor allele frequency (MAF) >5% and odds ratios (OR) of 1.5 with the type 1 error rate, alpha = 0.05. RESULTS: We found no evidence of association with T1D at most of the loci studied 0.02

Analysis of polymorphisms of the interleukin-18 gene in type 1 diabetes and Hardy-Weinberg equilibrium testing.

Recently, the interleukin-18 cytokine gene (IL18) was reported to be associated with type 1 diabetes. In the present report, we calculated that the reported genotypes of the two 5' region/promoter single nucleotide polymorphisms (SNPs), -607 (C-->A) (rs1946518) and -137 (G-->C) (rs187238), were not in Hardy-Weinberg equilibrium (HWE). We therefore investigated the association of the -607 and -137 SNPs in a U.K. type 1 diabetic Caucasian case-control collection (1,560 case and 1,715 control subjects tested at -607 and 4,323 case and 4,610 control subjects tested at -137) as well as a type 1 diabetic Caucasian collection comprised of families of European ancestry (1,347 families tested at -137 and 1,356 families tested at -607). No evidence for association with type 1 diabetes was found, including for the -607 A/A and C/A genotypes. To evaluate whether common variation elsewhere in the gene was associated with disease susceptibility, we analyzed eight IL18 tag SNPs in a type 1 diabetic case-control collection (1,561 case and 1,721 control subjects). No evidence for association was obtained (P = 0.11). We conclude that common allelic variation in IL18 is unlikely to contribute substantially to type 1 diabetes susceptibility in the populations tested and recommend routine application of tests for HWE in population-based studies for genetic association.

No evidence for association of the TATA-box binding protein glutamine repeat sequence or the flanking chromosome 6q27 region with type 1 diabetes.

Susceptibility to the autoimmune disease type 1 diabetes has been linked to human chromosome 6q27 and, moreover, recently associated with one of the genes in the region, TATA box-binding protein (TBP). Using a much larger sample of T1D families than those studied by others, and by extensive re-sequencing of nine other genes in the proximity, in which we identified 279 polymorphisms, 83 of which were genotyped in up to 725 T1D multiplex and simplex families, we obtained no evidence for association of the TBP CAG/CAA (glutamine) microsatellite repeat sequence with disease, or for nine other genes, PDCD2, PSMB1, KIAA1838, DLL1, dJ894D12.4, FLJ25454, FLJ13162, FLJ11152, PHF10 and CCR6. This study also provides an exon-based tag single nucleotide polymorphism map for these 10 genes that can be used for analysis of other diseases.

Localization of a type 1 diabetes locus in the IL2RA/CD25 region by use of tag single-nucleotide polymorphisms.

As part of an ongoing search for genes associated with type 1 diabetes (T1D), a common autoimmune disease, we tested the biological candidate gene IL2RA (CD25), which encodes a subunit (IL-2R alpha) of the high-affinity interleukin-2 (IL-2) receptor complex. We employed a tag single-nucleotide polymorphism (tag SNP) approach in large T1D sample collections consisting of 7,457 cases and controls and 725 multiplex families. Tag SNPs were analyzed using a multilocus test to provide a regional test for association. We found strong statistical evidence in the case-control collection (P=6.5x10(-8)) for a T1D locus in the CD25 region of chromosome 10p15 and replicated the association in the family collection (P=7.3x10(-3); combined P=1.3x10(-10)). These results illustrate the utility of tag SNPs in a chromosome-regional test of disease association and justify future fine mapping of the causal variant in the region.

Construction and analysis of tag single nucleotide polymorphism maps for six human-mouse orthologous candidate genes in type 1 diabetes.

BACKGROUND: One strategy to help identify susceptibility genes for complex, multifactorial diseases is to map disease loci in a representative animal model of the disorder. The nonobese diabetic (NOD) mouse is a model for human type 1 diabetes. Linkage and congenic strain analyses have identified several NOD mouse Idd (insulin dependent diabetes) loci, which have been mapped to small chromosome intervals, for which the orthologous regions in the human genome can be identified. Here, we have conducted re-sequencing and association analysis of six orthologous genes identified in NOD Idd loci: NRAMP1/SLC11A1 (orthologous to Nramp1/Slc11a1 in Idd5.2), FRAP1 (orthologous to Frap1 in Idd9.2), 4-1BB/CD137/TNFRSF9 (orthologous to 4-1bb/Cd137/Tnrfrsf9 in Idd9.3), CD101/IGSF2 (orthologous to Cd101/Igsf2 in Idd10), B2M (orthologous to B2m in Idd13) and VAV3 (orthologous to Vav3 in Idd18). RESULTS: Re-sequencing of a total of 110 kb of DNA from 32 or 96 type 1 diabetes cases yielded 220 single nucleotide polymorphisms (SNPs). Sixty-five SNPs, including 54 informative tag SNPs, and a microsatellite were selected and genotyped in up to 1,632 type 1 diabetes families and 1,709 cases and 1,829 controls. CONCLUSION: None of the candidate regions showed evidence of association with type 1 diabetes (P values > 0.2), indicating that common variation in these key candidate genes does not play a major role in type 1 diabetes susceptibility in the European ancestry populations studied.

No evidence of association or interaction between the IL4RA, IL4, and IL13 genes in type 1 diabetes.

Attempts to identify susceptibility loci that, on their own, have marginal main effects by use of gene-gene interaction tests have increased in popularity. The results obtained from analyses of epistasis are, however, difficult to interpret. Gene-gene interaction, albeit only marginally significant, has recently been reported for the interleukin-4 and interleukin-13 genes (IL4 and IL13) with the interleukin-4 receptor A gene (IL4RA), contributing to the susceptibility of type 1 diabetes (T1D). We aimed to replicate these findings by genotyping both large family and case-control data sets and by using previously published data. Gene-gene interaction tests were performed using linear regression models in cases only. We did not find any single-locus associations with T1D and did not obtain evidence of gene-gene interaction. Additional support from independent samples will be even more important in the study of gene-gene interactions and other subgroup analyses.

Development of an integrated genome informatics, data management and workflow infrastructure: a toolbox for the study of complex disease genetics.

The genetic dissection of complex disease remains a significant challenge. Sample-tracking and the recording, processing and storage of high-throughput laboratory data with public domain data, require integration of databases, genome informatics and genetic analyses in an easily updated and scaleable format. To find genes involved in multifactorial diseases such as type 1 diabetes (T1D), chromosome regions are defined based on functional candidate gene content, linkage information from humans and animal model mapping information. For each region, genomic information is extracted from Ensembl, converted and loaded into ACeDB for manual gene annotation. Homology information is examined using ACeDB tools and the gene structure verified. Manually curated genes are extracted from ACeDB and read into the feature database, which holds relevant local genomic feature data and an audit trail of laboratory investigations. Public domain information, manually curated genes, polymorphisms, primers, linkage and association analyses, with links to our genotyping database, are shown in Gbrowse. This system scales to include genetic, statistical, quality control (QC) and biological data such as expression analyses of RNA or protein, all linked from a genomics integrative display. Our system is applicable to any genetic study of complex disease, of either large or small scale.

Replication of an association between the lymphoid tyrosine phosphatase locus (LYP/PTPN22) with type 1 diabetes, and evidence for its role as a general autoimmunity locus.

In the genetic analysis of common, multifactorial diseases, such as type 1 diabetes, true positive irrefutable linkage and association results have been rare to date. Recently, it has been reported that a single nucleotide polymorphism (SNP), 1858C>T, in the gene PTPN22, encoding Arg620Trp in the lymphoid protein tyrosine phosphatase (LYP), which has been shown to be a negative regulator of T-cell activation, is associated with an increased risk of type 1 diabetes. Here, we have replicated these findings in 1,388 type 1 diabetic families and in a collection of 1,599 case and 1,718 control subjects, confirming the association of the PTPN22 locus with type 1 diabetes (family-based relative risk (RR) 1.67 [95% CI 1.46-1.91], and case-control odds ratio (OR) 1.78 [95% CI 1.54-2.06]; overall P = 6.02 x 10(-27)). We also report evidence for an association of Trp(620) with another autoimmune disorder, Graves' disease, in 1,734 case and control subjects (P = 6.24 x 10(-4); OR 1.43 [95% CI 1.17-1.76]). Taken together, these results indicate a more general association of the PTPN22 locus with autoimmune disease.

Association analysis of the lymphocyte-specific protein tyrosine kinase (LCK) gene in type 1 diabetes.

Prior data associating the expression of lymphocyte-specific protein tyrosine kinase (LCK) with type 1 diabetes, its critical function in lymphocytes, and the linkage of the region to diabetes in the nonobese diabetic (NOD) mouse model make LCK a premier candidate for a susceptibility gene. Resequencing of LCK in 32 individuals detected seven single nucleotide polymorphisms (SNPs) with allele frequencies >3%, including four common SNPs previously reported. These and six other SNPs from dbSNP were genotyped in a two-stage strategy using 2,430 families and were all shown not to be significantly associated with type 1 diabetes. We conclude that a major role for the common LCK polymorphisms in type 1 diabetes is unlikely. However, we cannot rule out the possibility of there being a causal variant outside the exonic, intronic, and untranslated regions studied.

Investigating the utility of combining phi29 whole genome amplification and highly multiplexed single nucleotide polymorphism BeadArray genotyping.

BACKGROUND: Sustainable DNA resources and reliable high-throughput genotyping methods are required for large-scale, long-term genetic association studies. In the genetic dissection of common disease it is now recognised that thousands of samples and hundreds of thousands of markers, mostly single nucleotide polymorphisms (SNPs), will have to be analysed. In order to achieve these aims, both an ability to boost quantities of archived DNA and to genotype at low costs are highly desirable. We have investigated phi29 polymerase Multiple Displacement Amplification (MDA)-generated DNA product (MDA product), in combination with highly multiplexed BeadArray genotyping technology. As part of a large-scale BeadArray genotyping experiment we made a direct comparison of genotyping data generated from MDA product with that from genomic DNA (gDNA) templates. RESULTS: Eighty-six MDA product and the corresponding 86 gDNA samples were genotyped at 345 SNPs and a concordance rate of 98.8% was achieved. The BeadArray sample exclusion rate, blind to sample type, was 10.5% for MDA product compared to 5.8% for gDNA. CONCLUSIONS: We conclude that the BeadArray technology successfully produces high quality genotyping data from MDA product. The combination of these technologies improves the feasibility and efficiency of mapping common disease susceptibility genes despite limited stocks of gDNA samples.

Lack of association of the Ala(45)Thr polymorphism and other common variants of the NeuroD gene with type 1 diabetes.

Variation in genes necessary for normal functioning and development of beta-cells, e.g., NEUROD1, which encodes a transcription factor for the insulin gene and is important in beta-cell development, causes maturity-onset diabetes of the young. Some studies have reported an association between a nonsynonymous Ala(45)Thr (+182G-->A) single nucleotide polymorphism (SNP) in NEUROD1 and type 1 diabetes, but this result has not been consistently found. To clarify this, we genotyped Ala(45)Thr in 2,434 type 1 diabetic families of European descent and Caucasian ethnicity from five different countries. Taking the allele frequency of 36% for Thr(45) and an odds ratio (OR) of 1.2, this sample provided >99% power to detect an association (P < 0.05). We could not confirm the association (P = 0.77). No evidence of population heterogeneity in the lack of association of Thr(45) with type 1 diabetes was observed. To evaluate the possibility that another NEUROD1 variant was associated with type 1 diabetes, we resequenced the gene in 32 U.K. affected individuals and identified and genotyped all common SNPs (minor allele frequency >10%; n = 5) in 786 families. We report no evidence of association of these common variants in NEUROD1 and type 1 diabetes in these samples.

Analysis of the type 2 diabetes-associated single nucleotide polymorphisms in the genes IRS1, KCNJ11, and PPARG2 in type 1 diabetes.

It has been proposed that type 1 and 2 diabetes might share common pathophysiological pathways and, to some extent, genetic background. However, to date there has been no convincing data to establish a molecular genetic link between them. We have genotyped three single nucleotide polymorphisms associated with type 2 diabetes in a large type 1 diabetic family collection of European descent: Gly972Arg in the insulin receptor substrate 1 (IRS1) gene, Glu23Lys in the potassium inwardly-rectifying channel gene (KCNJ11), and Pro12Ala in the peroxisome proliferative-activated receptor gamma2 gene (PPARG2). We were unable to confirm a recently published association of the IRS1 Gly972Arg variant with type 1 diabetes. Moreover, KCNJ11 Glu23Lys showed no association with type 1 diabetes (P > 0.05). However, the PPARG2 Pro12Ala variant showed evidence of association (RR 1.15, 95% CI 1.04-1.28, P = 0.008). Additional studies need to be conducted to confirm this result.

Haplotype tag single nucleotide polymorphism analysis of the human orthologues of the rat type 1 diabetes genes Ian4 (Lyp/Iddm1) and Cblb.

The diabetes-prone BioBreeding (BB) and Komeda diabetes-prone (KDP) rats are both spontaneous animal models of human autoimmune, T-cell-associated type 1 diabetes. Both resemble the human disease, and consequently, susceptibility genes for diabetes found in these two strains can be considered as potential candidate genes in humans. Recently, a frameshift deletion in Ian4, a member of the immune-associated nucleotide (Ian)-related gene family, has been shown to map to BB rat Iddm1. In the KDP rat, a nonsense mutation in the T-cell regulatory gene, Cblb, has been described as a major susceptibility locus. Following a strategy of examining the human orthologues of susceptibility genes identified in animal models for association with type 1 diabetes, we identified single nucleotide polymorphisms (SNPs) from each gene by resequencing PCR product from at least 32 type 1 diabetic patients. Haplotype tag SNPs (htSNPs) were selected and genotyped in 754 affected sib-pair families from the U.K. and U.S. Evaluation of disease association by a multilocus transmission/disequilibrium test (TDT) gave a P value of 0.484 for IAN4L1 and 0.692 for CBLB, suggesting that neither gene influences susceptibility to common alleles of human type 1 diabetes in these populations.

Haplotype structure, LD blocks, and uneven recombination within the LRP5 gene.

Patterns of linkage disequilibrium (LD) in the human genome are beginning to be characterized, with a paucity of haplotype diversity in "LD blocks," interspersed by apparent "hot spots" of recombination. Previously, we cloned and physically characterized the low-density lipoprotein-receptor-related protein 5 (LRP5) gene. Here, we have extensively analysed both LRP5 and its flanking three genes, spanning 269 kb, for single nucleotide polymorphisms (SNPs), and we present a comprehensive SNP map comprising 95 polymorphisms. Analysis revealed high levels of recombination across LRP5, including a hot-spot region from intron 1 to intron 7 of LRP5, where there are 109 recombinants/Mb (4882 meioses), in contrast to flanking regions of 14.6 recombinants/Mb. This region of high recombination could be delineated into three to four hot spots, one within a 601-bp interval. For LRP5, three haplotype blocks were identified, flanked by the hot spots. Each LD block comprised over 80% common haplotypes, concurring with a previous study of 14 genes that showed that common haplotypes account for at least 80% of all haplotypes. The identification of hot spots in between these LD blocks provides additional evidence that LD blocks are separated by areas of higher recombination.

Association of the T-cell regulatory gene CTLA4 with susceptibility to autoimmune disease.

Genes and mechanisms involved in common complex diseases, such as the autoimmune disorders that affect approximately 5% of the population, remain obscure. Here we identify polymorphisms of the cytotoxic T lymphocyte antigen 4 gene (CTLA4)--which encodes a vital negative regulatory molecule of the immune system--as candidates for primary determinants of risk of the common autoimmune disorders Graves' disease, autoimmune hypothyroidism and type 1 diabetes. In humans, disease susceptibility was mapped to a non-coding 6.1 kb 3' region of CTLA4, the common allelic variation of which was correlated with lower messenger RNA levels of the soluble alternative splice form of CTLA4. In the mouse model of type 1 diabetes, susceptibility was also associated with variation in CTLA-4 gene splicing with reduced production of a splice form encoding a molecule lacking the CD80/CD86 ligand-binding domain. Genetic mapping of variants conferring a small disease risk can identify pathways in complex disorders, as exemplified by our discovery of inherited, quantitative alterations of CTLA4 contributing to autoimmune tissue destruction.

Linkage and association mapping of the LRP5 locus on chromosome 11q13 in type 1 diabetes.

Linkage of chromosome 11q13 to type 1 diabetes (T1D) was first reported from genome scans (Davies et al. 1994; Hashimoto et al. 1994) resulting in P <2.2 x 10(-5) (Luo et al. 1996) and designated IDDM4 ( insulin dependent diabetes mellitus 4). Association mapping under the linkage peak using 12 polymorphic microsatellite markers suggested some evidence of association with a two-marker haplotype, D11S1917*03-H0570POLYA*02, which was under-transmitted to affected siblings and over-transmitted to unaffected siblings ( P=1.5 x 10(-6)) (Nakagawa et al. 1998). Others have reported evidence for T1D association of the microsatellite marker D11S987, which is approximately 100 kb proximal to D11S1917 (Eckenrode et al. 2000). We have sequenced a 400-kb interval surrounding these loci and identified four genes, including the low-density lipoprotein receptor related protein (LRP5) gene, which has been considered as a functional candidate gene for T1D (Hey et al. 1998; Twells et al. 2001). Consequently, we have developed a comprehensive SNP map of the LRP5 gene region, and identified 95 SNPs encompassing 269 kb of genomic DNA, characterised the LD in the region and haplotypes (Twells et al. 2003). Here, we present our refined linkage curve of the IDDM4 region, comprising 32 microsatellite markers and 12 SNPs, providing a peak MLS=2.58, P=5 x 10(-4), at LRP5 g.17646G>T. The disease association data, largely focused in the LRP5 region with 1,106 T1D families, provided no further evidence for disease association at LRP5 or at D11S987. A second dataset, comprising 1,569 families from Finland, failed to replicate our previous findings at LRP5. The continued search for the variants of the putative IDDM4 locus will greatly benefit from the future development of a haplotype map of the genome.