Evidence for Susceptibility Loci from Four Genome-Wide Linkage Scans in 1,435 Multiplex Families
http://www.100md.com
糖尿病学杂志 2005年第10期
1 Benaroya Research Institute, Seattle, Washington
2 Roche Molecular Systems, Alameda, California
3 Pasteur Institute, Paris, France
4 The Western Australian Institute of Medical Research, Perth, Australia
5 Steno Diabetes Center, Gentofte, Denmark
6 Juvenile Diabetes Research Foundation/Wellcome Trust Diabetes Inflammation Laboratory, Cambridge University, Cambridge, U.K
7 Wake Forest University School of Medicine, Winston-Salem, North Carolina
ABSTRACT
Type 1 diabetes is a common, multifactorial disease with strong familial clustering (genetic risk ratio [S] 15). Approximately 40% of the familial aggregation of type 1 diabetes can be attributed to allelic variation of HLA loci in the major histocompatibility complex on chromosome 6p21 (locus-specific S 3). Three other disease susceptibility loci have been clearly demonstrated based on their direct effect on risk, INS (chromosome 11p15, allelic odds ratio [OR] 1.9), CTLA4 (chromosome 2q33, allelic OR 1.2), and PTPN22 (chromosome 1p13, allelic OR 1.7). However, a large proportion of type 1 diabetes clustering remains unexplained. We report here on a combined linkage analysis of four datasets, three previously published genome scans, and one new genome scan of 254 families, which were consolidated through an international consortium for type 1 diabetes genetic studies (www.t1dgc.org) and provided a total sample of 1,435 families with 1,636 affected sibpairs. In addition to the HLA region (nominal P = 2.0 x 10eC52), nine noneCHLA-linked regions showed some evidence of linkage to type 1 diabetes (nominal P < 0.01), including three at (or near) genome-wide significance (P < 0.05): 2q31-q33, 10p14-q11, and 16q22-q24. In addition, after taking into account the linkage at the 6p21 (HLA) region, there was evidence supporting linkage for the 6q21 region (empiric P < 10eC4). More than 80% of the genome could be excluded as harboring type 1 diabetes susceptibility genes of modest effect (S 1.3) that could be detected by linkage. This study represents one of the largest linkage studies ever performed for any common disease. The results demonstrate some consistency emerging for the existence of susceptibility loci on chromosomes 2q31-q33, 6q21, 10p14-q11, and 16q22-q24 but diminished support for some previously reported locations.
Type 1 diabetes is the third most prevalent chronic disease of childhood, affecting up to 0.4% of children in some populations by age 30 years, with an overall lifetime risk of nearly 1% (1,2). It is believed that a large proportion of cases of type 1 diabetes result from the autoimmune destruction of the pancreatic cells, leading to complete dependence on exogenous insulin to regulate blood glucose levels (3). The etiology of type 1 diabetes is only partially characterized, but it is recognized that both genetic and environmental determinants are important in defining disease risk. Type 1 diabetes clusters in families, based on population-based twin and family studies (4) but does not segregate with a known mode of inheritance (5). The incidence and the age at onset of type 1 diabetes in some populations have changed dramatically since 1950 (6eC8). These data, coupled with the incomplete concordance for the phenotype in monozygotic twins (30eC70%), suggest that the penetrance of type 1 diabetes susceptibility alleles is strongly influenced by environmental factors (4).
Type 1 diabetes is strongly clustered in families with an overall genetic risk ratio (S) of 15 (9). At least one locus that contributes strongly to this familial clustering resides within the major histocompatibility complex (MHC) on chromosome 6p21. Genetic, functional, structural, and model studies all suggest that the HLA class II genes (HLA-DRB1 and -DQB1) likely represent the primary determinants of IDDM1. The frequency of HLA class II susceptibility alleles also correlates well with the population incidence of type 1 diabetes (10). These studies suggest that the MHC (IDDM1) may account for nearly 40% of the observed familial clustering of type 1 diabetes, with a locus-specific S of 3 (11).
Given that HLA alone cannot explain the familial clustering of type 1 diabetes, several, perhaps many, genes remain to be identified. First, in the general population, individuals who carry the high-risk haplotypic combination DRB104-DQB10302/DRB103-DQB10201 have 5% absolute risk of type 1 diabetes. However, within affected sibpair families, this genotype has 20% risk (5,12). Second, three non-HLA loci have been identified based on genetic association studies: IDDM2 [INS, 11p15 (13eC16)], IDDM12 [CTLA4, 2q33 (17)], and LYP/PTPN22 [1p13 (18eC20)]. Finally, the observed risk of type 1 diabetes in first- and second-degree relatives declines in a pattern consistent with multiplicative effects of multiple loci (11).
The first type 1 diabetes genome-wide scans for linkage, using fewer than 100 affected sibpair families, identified chromosome 6p21 (IDDM1) as the major type 1 diabetes risk locus (21,22). Subsequent studies supported non-HLA loci on chromosomes 11q13 (IDDM4) and 6q25 (IDDM5) in families from the U.K., the U.S., and France and on chromosome 15q26 (IDDM3) in families from Canada. Using both linkage and association approaches, other putative type 1 diabetes susceptibility loci were identified on chromosomes 18q12-q21 (IDDM6), 2q33 (IDDM7), 6q27 (IDDM8), 3q22-q25 (IDDM9), 10p11-q11 (IDDM10), 14q24-q31 (IDDM11), 2q31-q33 (IDDM12), 2q34-q35 (IDDM13), 6q21 (IDDM15), 14q32 (IDDM16), 10q25 (IDDM17), 5q33 (IDDM18), 7p15-p13 (GCK), 1q42, 16q22-q24, Xp11 (conditional on HLA-DR genotype), and 8q22-q24 (8eC10,12eC17). Even though statistical evidence supporting linkage for some of these regions was strong in the initial reports, most regions have not been clearly established in multiple populations (9,23).
A major barrier to type 1 diabetes gene identification, given the likely small locus-specific contribution (low S) for non-HLA genes, is the relatively small number of newly ascertained affected sibpair families with type 1 diabetes. In addition, previous studies have used available samples from a variety of collections, making the compilation of linkage results difficult because of apparent overlap in families analyzed and uncertainty in the equivalence of allele coding. To facilitate the genetic analysis of type 1 diabetes, results from two previous genome scans of type 1 diabetes were merged, comprising 767 families and 831 affected sibpairs from the U.K. and U.S. (24). The combined analyses supported linkage to at least six non-HLA regions to type 1 diabetes, including IDDM2 (INS, nominal P = 6.5 x 10eC4), 2q31-q33 (P = 5.1 x 10eC4), and 10p11 (P = 3.2 x 10eC4). A third genome scan of 424 type 1 diabetic families with at least two affected relative pairs (464 affected pairs) from Scandinavia (25) found no evidence for linkage at these latter three loci but did support linkage on chromosomes 5q11.2 (nominal P = 8.1 x 10eC4) and 16p13 (P = 1.6 x 10eC4). The IDDM15 region on chromosome 6q21 supported linkage (P = 7.0 x 10eC7) when HLA was taken into consideration in a combined analysis of U.S., French, and Scandinavian families.
We present here a joint analysis of data from these three prior genome-wide scans (U.S., U.K., and Scandinavia) as well as 254 new families collected for this study, a total of 1,435 multiplex families, for linkage to type 1 diabetes. With an average map information content of 67% (from 400 polymorphic microsatellite markers in each scan), this family collection provides 95% power to detect a locus with locus-specific S 1.3 and P = 10eC4. In addition to HLA, there was nominal evidence for linkage of type 1 diabetes to 10 other chromosome regions, including 6q21 (IDDM15) and 3 that reached genome-wide levels of significance, 2q31-q33 (IDDM12 and IDDM7), 10p11-q14 (IDDM10), and 16p12-q24. These data support the existence of non-HLA susceptibility loci for type 1 diabetes and strengthen support for a subset of loci previously proposed to contribute to type 1 diabetes risk.
RESEARCH DESIGN AND METHODS
Four sets of Caucasian families provided genome scan data for the combined analyses. Three sets of families have been previously published—U.K. (21,26), U.S. (24,26), and Scandinavia (25)—and one set of 254 families that were newly assembled for this study. The new collection of DNA samples from 254 families was obtained from several sources. DNA samples from 47 U.K. families were identified from the Diabetes U.K. Warren repository (28) that had not been genotyped previously. Families not previously used in published genome scans from the U.S. were contributed by the Joslin Diabetes Center (121 families) and from the Human Biological Data Interchange (76 families). Ten families from Australia were collected by investigators at the Walter and Eliza Hall Institute as previously described (29). In total, there were 1,435 families containing 6,899 individuals (6,358 with genome scan data). A total of 3,109 individuals were affected (with type 1 diabetes), and of these, 3,072 had genotype data (for details of the samples, see table in online appendix [available at http://diabetes.diabetesjournals.org]).
Genotyping.
Microsatellite marker genotyping technologies and allele scoring conventions varied between the different laboratories providing the data for previously published type 1 diabetic families. Details of genotyping of the U.K. (21), U.S. (27), and Scandinavian (25) families have been previously described. The 254 new families were genotyped by the Center for the Inheritance of Disease Research (http://www.cidr.jhmi.edu/) using a panel of 405 microsatellite markers. Because direct merging of genotypes (by standardized allele size) was not possible, within-family recoded genotype data from all four sources were merged into a single database. Family-naming conventions between the samples were normalized, and individual marker names were modified to indicate the laboratory of origin for the genotyping. After elimination of marker inconsistencies (see below), genetic markers were selected to form the analysis panel. Multiple independent genotyping occurred for some markers on a subset of individuals by different laboratories. In these cases, only one marker was used for the current analysis. Unless a marker showed inconsistencies in identity-by-descent (IBD) sharing, the marker that was included for analysis was the one scored in the most samples. A total of 1,190 markers were included in the combined analyses. Markers that previously had been added to maps because of association with type 1 diabetes were excluded to avoid bias in the multipoint linkage results. The excluded markers were located on chromosome 2 (alpha4, ND1, D2S152, CTLA4, and IGFBP5) and chromosome 11 (INS, TH).
Statistical analysis.
Before integration of the genetic data, marker error detection and pedigree structure within each dataset were made using PREST software (30). This method uses the genome scan data to determine the likelihood of each specified relationship given the genetic data. Unlikely marker genotypes were resolved by recoding the specific genotype to "unknown." Occurrences of nonpaternity were resolved by changing the pedigree structure to that which was most likely and then repeating the analysis to confirm appropriate relationships. An integrated marker map was developed by using public databases (Mammalian Genotyping Service, http://research.marshfieldclinic.org/genetics/Genotyping_Service/mgsver2.htm; Southampton, http://cedar.genetics.soton.ac.uk/public_html/; Cooperative Human Linkage Center, http://gai.nci.nih.gov/CHLC/; deCODE, http://www.decode.is) as well as physical map and genome sequence information from the University of California at Santa Cruz (http://genome.ucsc.edu/) using primer sequences in BLAST searches against the genome sequence. All analyses were based on this "consensus" map. Single and multipoint linkage analyses (based on the consensus map order and distances) were performed using GeneHunter-plus [S(pairs) option (31eC33)]. Examination of double recombinants was performed using Merlin software (34). Information content was estimated using Allegro (35).
Estimation of IBD statistics and resulting likelihoods under the null and alternative hypotheses were computed within each dataset. These dataset-specific likelihoods were then combined for the combined linkage analyses. Multipoint linkage analyses were performed, and maximized logarithm of odds (LOD) scores were calculated under an exponential model with constrained between 0 and 2 (32). Genome-wide empirical P values were determined by simulating Mendelian transmission with families maintaining the patterns of missing data observed in the sample (36). The fraction of LOD scores observed greater than the nominal value, based on simulations of 10,000 replicates across the genome, provided the estimated genome-wide P value. Exclusion mapping was performed using the MapMaker/SIBS program (37).
Previous studies have identified a potential type 1 diabetes susceptibility locus near (but not within) the HLA complex (IDDM15) (25,27,38). In the presence of strong support for linkage due to IDDM1, determining the support for IDDM15 is complex because of the positive correlation among the IBD proportions for linked loci. Without accounting for this positive correlation, statistical tests are biased toward inferring an epistatic relationship. A simple approach when IBD is known is to compare the observed correlation, r, between the IBD estimates at two loci (i.e., IDDM1 and IDDM15) with the theoretical correlation, = (1eC2)2, for two loci separated by recombination fraction, . The statistic t = (z eC )neC 3, where z = ln () and = ln () allows for the test of "interaction" that contrasts the observed IBD at IDDM15 based on the distance of IDDM15 from IDDM1 and the expected IBD given that distance. Using a sex-averaged map, IBD estimates for each sibpair in the data were computed, and a single IBD estimate was selected from each pedigree. Thus, each of these IBD estimates is independent. Under the null hypothesis of no interaction between these loci, t approximately has a standard normal distribution. Observed correlations greater than reflect increased sharing at IDDM15 over that expected from IBD at IDDM1 and the hypothesized genetic distance of IDDM15 from IDDM1. A series of 10,000 simulations were performed as described above (36), with the correlation in IBD computed between IDDM1 and IDDM15. The empirical P value for significance of IDDM15 was based on the number of simulated correlations greater than that observed in the original data.
RESULTS
Linkage analysis.
Analysis of 1,190 genetic markers in 1,435 families revealed that the strongest evidence for linkage to type 1 diabetes was on chromosome 6p21 (nominal P = 2.0 x 10eC52) in the MHC (Fig. 1).There were nine noneCHLA-linked regions with nominal evidence supporting linkage to type 1 diabetes (P < 0.01), including 2q31-q33 (IDDM7 and IDDM12; P = 9.0 x 10eC5; genome-wide P = 0.016), 10p14-q11 (IDDM10, P = 1.2 x 10eC4; genome-wide P = 0.021), and 16q22-q24 (P = 4.9 x 10eC4; genome-wide P = 0.075). The estimates of genetic relative risk varied by region, with IDDM1 having the largest locus-specific effect (S 3.3) and the other sites having low locus-specific effects (S 1.1eC1.2). Individual linkage plots by chromosome are provided in the online appendix (supplementary figures, upper curves).
Exclusion mapping.
Each of the four datasets had the equivalent of an 9-centiMorgan (cM) map, providing an average marker information content of 66%. The range of information from the markers was 62% for chromosome 15 to 73% for chromosome 16. In the combined data, 82% of the genome could be excluded at LOD < eC2 for loci of effect size S 1.3 (supplementary figures, lower curves), and >95% could be excluded for S 1.5. Several entire chromosomes could be excluded (chromosomes 7, 8, 18, 20, 21, and 22). The majority of chromosome 6 (due to the strongly linked 6p21/MHC region) and <50% of chromosomes 16, 19, and X could be excluded for S 1.3. For effects of S 1.1, only 6% of the genome could be excluded. The extent of exclusion for each chromosome is shown in Table 1.
IDDM15.
The t-statistic was computed to determine the increase in sharing at IDDM15, beyond that expected from sharing at IDDM1 and expected decay in sharing due to genetic distance. In the current data, the test statistic was computed for IDDM15 at a position 37 cM from IDDM1 (at 47 cM on chromosome 6), using a sex-averaged map. The test statistic was also computed for a range of map positions within 5 cM (32eC42 cM) with similar results. At the 37-cM distance from IDDM1 and IDDM15, the expected correlation coefficient between the IBD estimates under the null hypothesis was estimated as = 0.2276. The observed Pearson’s correlation coefficient between the IBD estimates at IDDM1 and IDDM15, using 1,401 informative pedigrees, was r = 0.3132 (empirical P < 1.0 x 10eC4). These results support an HLA-independent effect in the IDDM15 region.
DISCUSSION
In a previous combined analysis of U.K. and U.S. families (24), it was concluded that an effort to merge and jointly analyze existing families would be required to clarify the role of noneCHLA-linked loci in type 1 diabetes. In the present study, this effort has been achieved under the auspices of the Type 1 Diabetes Genetics Consortium (http://www.t1dgc.org). We have assembled families and merged data from three large genome scans and added new data from 254 families not previously scanned. This increased sample size has allowed the exclusion of >80% of the human genome for locus-specific, but population-independent, effects of S 1.3. In addition to continued support for type 1 diabetes susceptibility related to the MHC (IDDM1) and INS (IDDM2), we identified eight regions that supported noneCHLA-linked susceptibility. Furthermore, we identified three chromosomes that contained extensive areas for which linkage could not be excluded—chromosomes 16, 19, and X—that could benefit from further genotyping to increase information content and the analysis of additional families.
Three noneCHLA-linked regions provided support for linkage at, or near, the genome-wide level of significance (P < 0.05): chromosome 2q31-q33 (IDDM7 and IDDM12), 10p14-q11 (IDDM10), and 16q22-q24. These locations were unlikely to have occurred by chance, based on our simulations. Together with the six other noneCHLA-linked regions that exhibited evidence of linkage at nominal P < 0.01 (Table 2) and the 10-cM map, the data indicate a strong non-HLA genetic effect for type 1 diabetes (39). Furthermore, none of the three most strongly linked regions exhibited support for linkage in the 408 families from Scandinavia (25).
Strong support for linkage (nominal P = 7.0 x 10eC7), after taking into account linkage to HLA) to the IDDM15 locus (6q21), was observed previously in a combined analysis of French, U.S., and Scandinavian families (25,38). In the present study, we have obtained support for IDDM15 (empirical P < 1.0 x 10eC4). The ability to further define the effects of this locus will be facilitated by increased information content in the HLA region and in the region surrounding IDDM15 to better estimate the observed IBD sharing and model the residual linkage to the HLA region, including taking into account sex-specific genetic map differences.
The IDDM12 locus lies within the 2q31-q33 region and has been attributed to single nucleotide polymorphisms (SNPs) in the 3'-untranslated region of CTLA4 (15); however, the modest S value predicted from the odds ratios (ORs) (1.1eC1.2) of the disease-associated SNPs at CTLA4 (S 1.01) in type 1 diabetes seems unlikely to fully account for the magnitude of the observed evidence for linkage (regional S 1.19). This result suggests the presence of other loci in the 2q31-q33 region, if this linkage is confirmed in other future studies. Originally, IDDM7 at chromosome 2q33 was assigned on the basis of evidence of allelic association of the D2S152 microsatellite marker, but this association has not been substantiated. The evidence supporting linkage in the current study does not include the putative IDDM13 locus at chromosome 2q34-q35 (29). Linkage of type 1 diabetes to 10p14-q13 (IDDM10) is well supported by the current and past studies (24,26); however, there has been little follow-up other than association analyses of the functional candidate gene GAD2, suggesting that this gene is not a type 1 diabetes susceptibility locus (40,41). The observed locus-specific effects for 2q31-q33 (S 1.19) and 10p14-q11 (S 1.12) suggest that a single common susceptibility allele would have an allelic association (OR) 3. This effect should be identifiable in a fine-mapping association study using dense SNP maps across the regions.
Support for a type 1 diabetes susceptibility locus on chromosome 16p12-q11.1, which was observed independently in both the combined U.K. and U.S. families (nominal P = 4.5 x 10eC3) and in the Scandinavian families (P = 2 x 10eC4), remained in the present study (P = 3.3 x 10eC3). A recent analysis of four rheumatoid arthritis genome scans (42) reported evidence for linkage at chromosomes 6p21 (HLA; P = 2 x 10eC5) and 16p-cen (P = 0.004). Because rheumatoid arthritis, antithyroid autoimmune disease, and type 1 diabetes cluster in families more often than expected by chance (43), evidence for linkage for any one of these autoimmune diseases could be informative for others. Evidence for linkage in U.K. families with early-onset rheumatoid arthritis (44) to chromosome 16p has previously been demonstrated (P = 3.2 x 10eC4). The comparison between linkage scan results for type 1 diabetes and rheumatoid arthritis provides other interesting similarities. The largest, single, combined scan of rheumatoid arthritis families (45) reported significant linkage of rheumatoid arthritis to chromosome 6p21 (HLA; P = 5 x 10eC12), and some evidence (P < 0.005) of rheumatoid arthritis linked to six other regions (1q43, 6q21, 10q21, 12q12, 17p13, and 18q21). This overlap with potential type 1 diabetes susceptibility at 6q21, 12q12, and 16p-cen may not be coincidental in the etiology of these autoimmune diseases.
Recently, evidence for association of type 1 diabetes with alleles in the PTPN22 locus (chromosome 1p13) has been reported (18), and this association has been confirmed (19,20). PTPN22 encodes a lymphoid-specific tyrosine phosphatase (LYP) and is also associated with autoimmune thyroid disease, rheumatoid arthritis, and SLE (46). The absence of evidence supporting linkage of type 1 diabetes to chromosome 1p13 (D1S206) in the 1,435 families studied here is not surprising given the magnitude of the PTPN22 association with type 1 diabetes. The OR of PTPN22 is large (1.7) but the S is 1.05. Thus, to detect linkage at P < 0.001 with 50% power, a sample of >8,000 affected sibpair families would be required, using a fully informative genetic map. Assuming a multiplicative model, the contribution of PTPN22 to type 1 diabetes (based on the observed OR) is 2%, much lower than HLA (40eC50%). Nevertheless, the knowledge that PTPN22 may be involved in risk to type 1 diabetes, as well as in other autoimmune diseases, is significant and could provide insight into modulating T-cell activity for disease prevention.
Several previously supported regions of linkage have diminished support in the current analyses. Type 1 diabetes susceptibility locus on chromosome 1q42 was strongly supported (nominal P = 9.8 x 10eC5) in a study of 679 U.K. and U.S. families (24,47) but exhibited decreasing linkage support in a follow-up analysis of 616 families using a denser map in the region (P = 4.0 x 10eC4). The evaluation of previously supported regions is difficult, even in the present study with over 1,600 affected sibpairs, particularly for regions with low s (e.g., s 1.1). For example, the region on 1q42 was originally supported with MLS = 3.31 and s 1.5 (27). Although the current sample excludes this region at the reported s 1.5, support for linkage in this region has decreased with increasing sample size from a LOD = 2.20 (24), to the current LOD = 0.87 (nominal P = 1.4 x 10eC2) with s 1.05. At the current estimated magnitude of genetic effect, the 1q42 region could not be excluded.
In a combined analysis of animal model and human linkage data from a number of autoimmune diseases, chromosome 18q12-q21 demonstrated evidence of linkage, which has now been supported by the analysis of congenic strains in NOD mice (48). There was no support for loci on chromosome 18 at P < 0.05 in the current study. Three independent studies of type 1 diabetes have reported linkage to chromosome 8q (49eC51), but there was no support for 8q in this study. Additional previously reported loci with relatively little support for linkage in the current study include IDDM4 (11q13), IDDM6 (18q12-q21), IDDM9 (3q22-q25), IDDM11 (14q24-q31), IDDM16 (14q32), IDDM17 (10q25), and IDDM18 (5q33). These data suggest that these putative type 1 diabetes susceptibility loci represent either false-positive results or have very small effects that may be more readily detected in certain populations because of variation in allele frequencies or other factors, including the possibility of population-specific genetic or environmental effects.
Linkage and fine mapping studies in mouse models of SLE (52) and of type 1 diabetes (53) have demonstrated that a single linkage peak may be composed of several susceptibility loci. In human populations, a linkage signal may be observed by the chance clustering of several disease loci, each with relatively weak locus-specific effects. The presence of multiple susceptibility loci may also account, in part, for broad linkage peaks often observed in studies of complex, common diseases. This underlying complexity would also increase the difficulty to obtain convincing results in future fine-mapping association studies. Both development of novel analytical approaches and increased sample size will be necessary to resolve this apparent complexity (54,55). Through the efforts of consortia (such as our international effort), it will be possible to increase the number of families for type 1 diabetes (http://www.t1dgc.org), which would increase power and allow exclusion of loci s 1.2, as well as provide standardized samples and reagents for future fine-mapping studies.
These results suggest two parallel tracks for the identification of type 1 diabetes susceptibility loci. First, systematic fine mapping of all the variants responsible for the HLA linkage to type 1 diabetes is justified, especially within the 4-Mb HLA region. Second, further exploration of potential non-HLA regions described here is now justified, including chromosomes 2q31-q33, 6q21, 10p14-q11, and 16q22-q24. Because sample sizes in linkage and association studies have historically been small and few genes (out of 25,000) have been studied in depth, future collaborative efforts and establishment of accessible resources for study should increase the yield of true disease susceptibility loci.
Appendix
Members of the Type 1 Diabetes Genetics Consortium
Asia-Pacific Network: Eesh Bhatia, Francois Bonnici, Pik To Cheung, Peter Colman, Andrew Cotterill, Jenny Couper, Ric Cutfield, Elizabeth Davis, Tim Davis, Paul Dixon, Kim Donaghue, Paul Drury, Mark Harris, Len Harrison, Tim Jones, Uma Kanga, Alok Kanungo, Betty Kek, Jeremy Krebs, Yann-Jinn Lee, Margaret Lloyd, Amanda Loth, Narinder Mehra, Grant Morahan, C.B. Sanjeevi, Brian Tait, Mike Varney, Jinny Willis, Loke Kah Yin, Goh Siok Ying.
European Network: Francisco J. Ampudia-Blasco, Jesus Argente, Magdalena Avbelj, Gulja Babadjanova, Klaus Badenhoop, Lubomir Barak, Christos Bartsocas, Tadej Battelino, Emilia Belda, Polly Bingley, Bernhard O. Boehm, Eulalia Brugues, Raffaella Buzzetti, Joyce Carlson, Luis Castano, Anna Casu, Ondrej Cinek, Raquel Corripio, Alberto de Leiva, V.A. Ruiz Esquide, Ana Fagulha, Marta Hernandez Garcia, Cristian Guja, Alona Hamou, Erifili Hatziagelaki, Simon Heath, Kaire Heilman, Nora Hosszufalusi, Constantin Ionescu-Tirgoviste, Cecile Julier, Ida Kinalska, Ingrid Kockum, Kalinka Koprivarova, Adam Kretowski, Dora Krikovszky, Nebojsa Lalic, Nicole Lambracht, Merce Lara, Mark Lathrop, Katharina Laubner, Claire Levy-Marchal, Johnny Ludvigsson, Laszlo Madacsy, Mara Marga, Didac Mauricio, Jorn Nerup, Antanas Norkus, Anna Okruszko, Xavier Palomer, Teresa Pedro, Moshe Phillip, Valdis Pirags, Flemming Pociot, Galina Popova, Paolo Pozzilli, Bart O. Roep, Silke Rosinger, Ana Maria Varela Sande, Ilhan Satman, Edith Schober, Jochen Seufert, Jan Skrha, Gyula Soltesz, Giatgen Spinas, Juraj Stanik, Tanya Szendeffy, Vallo Tillmann, Dag Undlien, Vaidotas Urbanavicius, Luciana Valente, Bart Van der Auwera, Andriani Vazeou-Gerasimidi, Dzilda Velickiene, Ana Wagner, Markus Walter, Alistair Williams, Miroslav Wurzburger, Anette-G Ziegler.
North American Network: Alan Aldrich, Mark Anderson, Christophe Benoist, Noureddine Berka, Patrick Concannon, Mark Daly, Jayne Danska, Larry Dolan, David Donaldson, Alessandro Doria, Janice Dorman, George Eisenbarth, Henry Erlich, Pamela Fain, Rosanna Fiallo-Scharer, Kenneth Gabbay, Daniel Geraghty, Soumitra Ghosh, Steven Gitelman, Nat Goodman, Gregory Goodwin, Carla Greenbaum, William Hagopian, John Hansen, Joel Hirschhorn, Leroy Hood, Kevin Kaiserman, Jean Lawrence, Victoria Magnuson, Jennifer Marks, John Mayberry, Elizabeth Mayer-Davis, Richard McIndoe, Brad McNeney, Eric Mickelson, Antoinette Moran, Gerald Nepom, Janelle Noble, Jill Norris, Tihamer Orban, David Owerbach, Andrew Paterson, Catherine Pihoker, Constantin Polychronakos, Alberto Pugliese, Philip Raskin, Marian Rewers, Henry Rodriquez, Jerry Rotter, Monique Roy, Desmond Schatz, Gary Schoch, Jin-Xiong She, Richard Spielman, Andrea Steck, Kent Taylor, Jay Tischfield, Ellen Toth, Diane Wherrett, Stephen Willi, Darrell Wilson, Lue Ping Zhao.
U.K. Network: Francesco Cucca, David Dunger, Simon Howell, Sarah Nutland, Helen Rance, Luc Smink, John Todd, Neil Walker, Barry Widmer, Heather Withers.
Coordinating Center: Don Babcock, Stephanie Beck, Mark Brown, Cralen Davis, Mark Espeland, Clara Gorodezky, Jason Griffin, Mark Hall, Teresa Harnish, Laura Hemrick, John Hepler, Joan Hilner, Letitia Howard, Ethan Lange, Carl Langefeld, Josyf Mychaleckyj, June Pierce, David Reboussin, Stephen Rich, Scott Rushing, Michele Sale, Elizabeth Sides, Beverly Snively, Michael Steffes, Augy Thiel, Lynne Wagenknecht, Dustin Williams, Jianzhao Xu, Belinda Youngdahl.
ACKNOWLEDGMENTS
P.C. has received support from the National Institute of Diabetes and Digestive and Kidney Diseases and Juvenile Diabetes Research Foundation. S.S.R. has received support from the National Institute of Diabetes and Digestive and Kidney Diseases and the JDRF. J.A.T. has received support from the JDRF and the Wellcome Trust.
We acknowledge the assistance of Mark Brown, Mathew Barber, Heather Cordell, Neil Walker, Carl Langefeld and Nancy Cox.
FOOTNOTES
A complete list of the members of the Type 1 Diabetes Genetics Consortium can be found in the APPENDIX.
Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org.
IBD, identity-by-descent; LOD, logarithm of odds; MHC, major histocompatibility complex; SNP, single nucleotide polymorphism
REFERENCES
Rewers M, LaPorte RE, King H, Tuomilehto J: Trends in the prevalence and incidence of diabetes: insulin-dependent diabetes mellitus in childhood. World Health Stat Q41 :179 eC189,1988
Rewers MJ: The changing face of the epidemiology of insulin-dependent diabetes mellitus (IDDM): research designs and models of disease causation. Ann Med23 :419 eC426,1991
Genuth S, Alberti KG, Bennett P, Buse J, Defronzo R, Kahn R, Kitzmiller J, Knowler WC, Lebovitz H, Lernmark A, Nathan D, Palmer J, Rizza R, Saudek C, Shaw J, Steffes M, Stern M, Tuomilehto J, Zimmet P; Expert Committee on the Diagnosis and Classification of Diabetes Mellitus: Follow-up report on the diagnosis of diabetes mellitus. Diabetes Care31 :3160 eC3167,2003
Hirshhorn JN: Genetic epidemiology of type 1 diabetes. Pediatr Diabetes4 :87 eC100,2003
Noble JA, Valdes AM, Cook M, Klitz W, Thomson G, Erlich HA: The role of HLA class II genes in insulin-dependent diabetes mellitus: molecular analysis of 180 Caucasian, multiplex families. Am J Hum Genet59 :1134 eC1148,1996
Hyttinen V, Kaprio J, Kinnunen L, Koskenvuo M, Tuomilehto J: Genetic liability of type 1 diabetes and the onset age among 22,650 young Finnish twin pairs: a nationwide follow-up study. Diabetes52 :1052 eC1055,2003
Svensson J, Carstensen B, Molbak A, Christau B, Mortensen HB, Nerup J, Borch-Johnsen K: Increased risk of childhood type 1 diabetes in children born after 1985. Diabetes Care25 :2197 eC2201,2002
Onkamo P, Vaananen S, Karvonen M, Tuomilehto J: Worldwide increase in incidence of type 1 diabetes: the analysis of the data on published incidence trends. Diabetologia42 :1395 eC1403,1999
Pociot F, McDermott MF: Genetics of type 1 diabetes mellitus (Review). Genes Immun3 :235 eC249,2002
Ronningen KS, Keiding N, Green A, EURODIAB ACE Study Group: Europe and diabetes: correlations between the incidence of childhood-onset type I diabetes in Europe and HLA genotypes. Diabetologia44 (Suppl. 3) :B51 eCB59,2001
Rich SS: Mapping genes in diabetes: a genetic epidemiological perspective. Diabetes39 :1315 eC1319,1990
Valdes AM, Noble JA, Genin E, Clerget-Darpoux F, Erlich HA, Thomson G: Modeling of HLA class II susceptibility to type 1 diabetes reveals an effect associated with DPB1. Genet Epidemiol21 :212 eC223,2001
Bell GI, Horita S, Karam JH: A polymorphic locus near the human insulin gene is associated with insulin-dependent diabetes mellitus. Diabetes33 :176 eC183,1984
Julier C, Hyer RN, Daviews J, Merlin F, Soularue P, Briant L, Cathelineau G, Deschamps I, Rotter JI, Froguel P: Insulin-IGF2 region on chromosome 11p encodes a gene implicated in HLA-DR4-dependent diabetes susceptibility. Nature354 :155 eC159,1991
Julier C, Lucassen A, Villedieu P, Delepine M, Levy-Marchal C, Danze PM, Bianchi F, Boitard C, Froguel P, Bell J: Multiple DNA variant association analysis: application to the insulin gene region in type I diabetes. Am J Hum Genet55 :1247 eC1254,1994
Barratt BJ, Payne F, Lowe CE, Hermann R, Healy BC, Harold D, Concannon P, Gharani N, McCarthy MI, Olavesen MG, McCormack R, Guja C, Ionescu-Tirgoviste C, Undlien DE, Ronningen KS, Gillespie KM, Tuomilehto-Wolf E, Tuomilehto J, Bennett ST, Clayton DG, Cordell HJ, Todd JA: Remapping the insulin gene/IDDM2 locus in type 1 diabetes. Diabetes53 :1884 eC1889,2004
Ueda H, Howson JM, Esposito L, Heward J, Snook H, Chamberlain G, Rainbow DB, Hunter KM, Smith AN, Di Genova G, Herr MH, Dahlman I, Payne F, Smyth D, Lowe C, Twells RC, Howlett S, Healy B, Nutland S, Rance HE, Everett V, Smink LJ, Lam AC, Cordell HJ, Walker NM, Bordin C, Hulme J, Motzo C, Cucca F, Hess JF, Metzker ML, Rogers J, Gregory S, Allahabadia A, Nithiyananthan R, Tuomilehto-Wolf E, Tuomilehto J, Bingley P, Gillespie KM, Undlien DE, Ronningen KS, Guja C, Ionescu-Tirgoviste C, Savage DA, Maxwell AP, Carson DJ, Patterson CC, Franklyn JA, Clayton DG, Peterson LB, Wicker LS, Todd JA, Gough SC: Association of the T-cell regulatory gene CTLA4 with susceptibility to autoimmune disease. Nature423 :506 eC511,2003
Bottini N, Musumeci L, Alonso A, Rahmouni S, Nika K, Rostamkhani M, MacMurray J, Meloni GF, Lucarelli P, Pellecchia M, Eisenbarth GS, Comings D, Mustelin T: A functional variant of lymphoid tyrosine phosphatase is associated with type I diabetes. Nat Genet36 :337 eC338,2004
Smyth DJ, Cooper JD, Collins JE, Heward JM, Franklyn JA, Howson JMM, Walker N, Vella A, Nutland S, Rance HE, Maier LM, Barratt BJ, Guja C, Ionescu-Tirgoviste CA, Savage D, Dunger DB, Widmer B, Strachan DP, Ring SM, Clayton DG, Twells RCJ, Gough SCL, Todd JA: 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. Diabetes53 :3020 eC3023,2004
Onengut-Gumuscu S, Ewens KG, Spielman RS, Concannon P: A functional polymorphism (1858C/T) in the PTPN22 gene is linked and associated with type 1 diabetes in multiplex families. Genes Immun5 :678 eC680,2004
Davies JL, Kawaguchi Y, Bennett ST, Copeman JB, Cordell HJ, Pritchard LE, Reed PW, Gough SCL, Jenkins SC, Palmer SM, Balfour KM, Row B, Farrall M, Barnett AH, Bain SC, Todd JA: A genome-wide search for human type 1 diabetes susceptibility genes. Nature371 :130 eC136,1994
Hashimoto L, Habita C, Beressi J, Delepine M, Besse C, Cambon-Thomsen A, Deschamps I, Rotter J, Djoulah S, James M, Froguel P, Weissenbach J, Lathrop GM, Julier C: Genetic mapping of a susceptibility locus for insulin-dependent diabetes mellitus on chromosome 11q. Nature371 :161 eC164,1994
Melanitou E, Fain P, Eisenbarth GS: Genetics of type 1A (immune mediated) diabetes. J Autoimmun21 :93 eC98,2003
Cox NJ, Wapelhorst B, Morrison VA, Johnson L, Pinchuk L, Spielman RS, Todd JA, Concannon P: Seven regions of the genome show evidence of linkage to type 1 diabetes in a consensus analysis of 767 multiplex families. Am J Hum Genet69 :820 eC830,2001
Mein CA, Esposito L, Dunn MG, Johnson GC, Timms AE, Goy JV, Smith AN, Sebag-Montefiore L, Merriman ME, Wilson AJ, Pritchard LE, Cucca F, Barnett AH, Bain SC, Todd JA: A search for type 1 diabetes susceptibility genes in families from the United Kingdom. Nat Genet19 :297 eC300,1998
Concannon P, Gogolin-Ewens KJ, Hinds D, Wapelhorst B, Morrison VA, Stirling B, Mitra M, Farmer J, Williams SR, Cox NJ, Bell GI, Risch N, Spielman RS: A second-generation screen of the human genome for susceptibility to insulin-dependent diabetes mellitus (IDDM). Nat Genet19 :292 eC296,1998
European Consortium for IDDM Genome Studies: A genomewide scan for type 1-diabetes susceptibility in Scandinavian families: identification of new loci with evidence of interactions. Am J Hum Genet69 :1301 eC1313,2001
Bain SC, Todd JA, Barnett AH: The British Diabetic Association: Warren repository. Autoimmunity7 :83 eC85,1990
Morahan G, Huang D, Tait B, Colman PG, Harrison LC: Markers on distal chromosome 2q linked to insulin-dependent diabetes mellitus. Science272 :1811 eC1813,1996
McPeek MS, Sun L: Statistical tests for detection of misspecified relationships by use of genome screen data. Am J Hum Genet66 :1076 eC1094,2000
Kruglyak L, Daly MJ, Reeve-Daly MP, Lander ES: Parametric and non-parametric linkage analysis: a unified multipoint approach. Am J Hum Genet58 :1347 eC1363,1996
Kong A, Cox NJ: Allele-sharing models: LOD scores and accurate linkage tests. Am J Hum Genet61 :1179 eC1188,1997
Langefeld CD, Boehnke M: Multiple trait nonparametric linkage regression. Am J Hum Genet65 :A258 ,1999
Abecasis GR, Cherny SS, Cookson WOC, Cardon LR: MERLIN: rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet30 :97 eC101,2002
Gudbjartsson DF, Jonasson K, Frigge M, Kong A: Allegro, a new computer program for multipoint linkage analysis. Nat Genet25 :12 eC13,2000
Hauser ER, Boehnke M: Genetic linkage analysis of complex genetic traits by using affected sibling pairs. Biometrics54 :1238 eC1246,1998
Kruglyak L, Lander ES: Complete multipoint sib-pair analysis of qualitative and quantitative traits. Am J Hum Genet57 :439 eC454,1995
Delepine M, Pociot F, Habita C, Hashimoto L, Froguel P, Rotter J, Combon-Thomsen A, Deschamps I, Djoulah S, Weissenbach J, Nerup J, Lathrop M, Julier C: Evidence of a non-MHC susceptibility locus in type 1 diabetes linked to HLA on chromosome 6. Am J Hum Genet60 :174 eC187,1997
Cucca F, Goy JV, Kawaguchi Y, Esposito L, Merriman ME, Wilson AJ, Cordell HJ, Bain SC, Todd JA: A male-female bias in type 1 diabetes and linkage to chromosome Xp in MHC HLA-DR3-positive patients. Nat Genet3 :301 eC302,1998
Wiltshire S, Cardon LR, McCarthy MI: Evaluating the results of genomewide linkage scans of complex traits by locus counting. Am J Hum Genet71 :1175 eC1182,2002
Wapelhorst B, Bell GI, Risch N, Spielman RS, Concannon P: Linkage and association studies in insulin-dependent diabetes with a new dinucleotide repeat polymorphism at the GAD65 locus. Autoimmunity21 :127 eC130,1995
Johnson GC, Payne F, Nutland S, Stevens H, Tuomilehto-Wolf E, Tuomilehto J, Todd JA: A comprehensive, statistically powered analysis of GAD2 in type 1 diabetes. Diabetes51 :2866 eC2870,2002
Fisher SA, Lanchbury JS, Lewis CM: Meta-analysis of four rheumatoid arthritis genome-wide linkage studies: confirmation of a susceptibility locus on chromosome 16. Arthritis Rheum48 :1200 eC1206,2003
Tait KF, Marshall T, Berman J, Carr-Smith J, Rowe B, Todd JA, Bain SC, Barnett AH, Gough SC: Clustering of autoimmune disease in parents of siblings from the type 1 diabetes Warren repository. Diabet Med21 :358 eC362,2004
MacKay K, Eyre S, Myerscough A, Millicic A, Barton A, Laval S, Barrett J, Lee D, While S, Jhon S, Brown MA, Bell J, Silman A, Ollier W, Wordsworth P, Worthington J: Whole-genome linkage analysis of rheumatoid arthritis susceptibility loci in 252 affected sibling pairs in the United Kingdom. Arthritis Rheum46 :632 eC639,2002
Jawaheer D, Seldin MF, Amos CI, Chen WV, Shigeta R, Etzel C, Damle A, Xiao X, Chen D, Lum RF, Monteiro J, Kern M, Criswell LA, Albani S, Nelson JL, Clegg DO, Pope R, Schroeder HW Jr, Bridges SL Jr, Pisetsky DS, Ward R, Kastner DL, Wilder RL, Pincus T, Callahan LF, Flemming D, Wener MH, Gregersen PK, North American Rheumatoid Arthritis Consortium: Screening the genome for rheumatoid arthritis susceptibility genes: a replication study and combined analysis of 512 multicase families. Arthritis Rheum48 :906 eC916,2003
Begovich AB, Carlton VE, Honigberg LA, Schrodi SJ, Chokkalingam AP, Alexander HC, Ardlie KG, Huang Q, Smith AM, Spoerke JM, Conn MT, Chang M, Chang SY, Saiki RK, Catanese JJ, Leong DU, Garcia VE, McAllister LB, Jeffery DA, Lee AT, Batliwalla F, Remmers E, Criswell LA, Seldin MF, Kastner DL, Amos CI, Sninsky JJ, Gregersen PK: A missense single-nucleotide polymorphism in a gene encoding a protein pyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am J Hum Genet75 :330 eC337,2004
Ewens KG, Johnson LN, Wapelhorst B, O’Brien K, Gutin S, Morrison VA, Street C, Gregory SG, Spielman RS, Concannon P: Linkage and association with type 1 diabetes on chromosome 1q42. Diabetes51 :3318 eC3325,2002
Merriman T, Twells R, Merriman M, Eaves I, Cox R, Cucca F, McKinney P, Shield J, Baum D, Bosi E, Pozzilli P, Nistico L, Buzzetti R, Joner G, Ronningen K, Thorsby E, Undlien D, Pociot F, Nerup J, Bain S, Barnett A, Todd J: Evidence by allelic association-dependent methods for a type 1 diabetes polygene (IDDM6) on chromosome 18q21. Hum Mol Genet6 :1003 eC1010,1997
Luo DF, Buzzetti R, Rotter JI, Maclaren NK, Raffel LJ, Nistico L, Govannini C, Pozzilli P, Thomson G, She JX: Confirmation of three susceptibility genes to insulin-dependent diabetes mellitus: IDDM4, IDDM5 and IDDM8. Hum Mol Genet5 :693 eC698,1996
Sale MM, FitzGerald LM, Charlesworth JC, Bowden DW, Rich SS: Evidence for a novel type 1 diabetes susceptibility locus on chromosome 8. Diabetes51 (Suppl. 3) :S316 eCS319,2002
Vaessen N, Heutink P, Houwing-Duistermaat JJ, Snijders PJ, Rademaker T, Testers L, Batstra MR, Sandkuijl LA, van Duijn CM, Oostra BA: A genome-wide search for linkage disequilibrium with type 1 diabetes in a recent genetically isolated population from the Netherlands. Diabetes51 :856 eC859,2002
Nguyen C, Limaye N, Wakeland EK: Susceptibility genes in the pathogenesis of murine lupus. Arthritis Res4 (Suppl. 3) :S255 eCS263,2002
Todd JA, Wicker LS: Genetic protection from the inflammatory disease type 1 diabetes in humans and animal models. Immunity15 :387 eC395,2001
Rich SS, Concannon P: Challenges and strategies for investigating the genetic complexity of common human diseases. Diabetes51 (Suppl. 3) :S288 eCS294,2002
Pociot F, Karlsen AE, Pedersen CB, Aalund M, Nerup J for the European Consortium for IDDM Genome Studies: Novel analytical methods applied to type 1 diabetes genome-scan data. Am J Hum Genet74 :647 eC660,2004(Patrick Concannon, Henry )
2 Roche Molecular Systems, Alameda, California
3 Pasteur Institute, Paris, France
4 The Western Australian Institute of Medical Research, Perth, Australia
5 Steno Diabetes Center, Gentofte, Denmark
6 Juvenile Diabetes Research Foundation/Wellcome Trust Diabetes Inflammation Laboratory, Cambridge University, Cambridge, U.K
7 Wake Forest University School of Medicine, Winston-Salem, North Carolina
ABSTRACT
Type 1 diabetes is a common, multifactorial disease with strong familial clustering (genetic risk ratio [S] 15). Approximately 40% of the familial aggregation of type 1 diabetes can be attributed to allelic variation of HLA loci in the major histocompatibility complex on chromosome 6p21 (locus-specific S 3). Three other disease susceptibility loci have been clearly demonstrated based on their direct effect on risk, INS (chromosome 11p15, allelic odds ratio [OR] 1.9), CTLA4 (chromosome 2q33, allelic OR 1.2), and PTPN22 (chromosome 1p13, allelic OR 1.7). However, a large proportion of type 1 diabetes clustering remains unexplained. We report here on a combined linkage analysis of four datasets, three previously published genome scans, and one new genome scan of 254 families, which were consolidated through an international consortium for type 1 diabetes genetic studies (www.t1dgc.org) and provided a total sample of 1,435 families with 1,636 affected sibpairs. In addition to the HLA region (nominal P = 2.0 x 10eC52), nine noneCHLA-linked regions showed some evidence of linkage to type 1 diabetes (nominal P < 0.01), including three at (or near) genome-wide significance (P < 0.05): 2q31-q33, 10p14-q11, and 16q22-q24. In addition, after taking into account the linkage at the 6p21 (HLA) region, there was evidence supporting linkage for the 6q21 region (empiric P < 10eC4). More than 80% of the genome could be excluded as harboring type 1 diabetes susceptibility genes of modest effect (S 1.3) that could be detected by linkage. This study represents one of the largest linkage studies ever performed for any common disease. The results demonstrate some consistency emerging for the existence of susceptibility loci on chromosomes 2q31-q33, 6q21, 10p14-q11, and 16q22-q24 but diminished support for some previously reported locations.
Type 1 diabetes is the third most prevalent chronic disease of childhood, affecting up to 0.4% of children in some populations by age 30 years, with an overall lifetime risk of nearly 1% (1,2). It is believed that a large proportion of cases of type 1 diabetes result from the autoimmune destruction of the pancreatic cells, leading to complete dependence on exogenous insulin to regulate blood glucose levels (3). The etiology of type 1 diabetes is only partially characterized, but it is recognized that both genetic and environmental determinants are important in defining disease risk. Type 1 diabetes clusters in families, based on population-based twin and family studies (4) but does not segregate with a known mode of inheritance (5). The incidence and the age at onset of type 1 diabetes in some populations have changed dramatically since 1950 (6eC8). These data, coupled with the incomplete concordance for the phenotype in monozygotic twins (30eC70%), suggest that the penetrance of type 1 diabetes susceptibility alleles is strongly influenced by environmental factors (4).
Type 1 diabetes is strongly clustered in families with an overall genetic risk ratio (S) of 15 (9). At least one locus that contributes strongly to this familial clustering resides within the major histocompatibility complex (MHC) on chromosome 6p21. Genetic, functional, structural, and model studies all suggest that the HLA class II genes (HLA-DRB1 and -DQB1) likely represent the primary determinants of IDDM1. The frequency of HLA class II susceptibility alleles also correlates well with the population incidence of type 1 diabetes (10). These studies suggest that the MHC (IDDM1) may account for nearly 40% of the observed familial clustering of type 1 diabetes, with a locus-specific S of 3 (11).
Given that HLA alone cannot explain the familial clustering of type 1 diabetes, several, perhaps many, genes remain to be identified. First, in the general population, individuals who carry the high-risk haplotypic combination DRB104-DQB10302/DRB103-DQB10201 have 5% absolute risk of type 1 diabetes. However, within affected sibpair families, this genotype has 20% risk (5,12). Second, three non-HLA loci have been identified based on genetic association studies: IDDM2 [INS, 11p15 (13eC16)], IDDM12 [CTLA4, 2q33 (17)], and LYP/PTPN22 [1p13 (18eC20)]. Finally, the observed risk of type 1 diabetes in first- and second-degree relatives declines in a pattern consistent with multiplicative effects of multiple loci (11).
The first type 1 diabetes genome-wide scans for linkage, using fewer than 100 affected sibpair families, identified chromosome 6p21 (IDDM1) as the major type 1 diabetes risk locus (21,22). Subsequent studies supported non-HLA loci on chromosomes 11q13 (IDDM4) and 6q25 (IDDM5) in families from the U.K., the U.S., and France and on chromosome 15q26 (IDDM3) in families from Canada. Using both linkage and association approaches, other putative type 1 diabetes susceptibility loci were identified on chromosomes 18q12-q21 (IDDM6), 2q33 (IDDM7), 6q27 (IDDM8), 3q22-q25 (IDDM9), 10p11-q11 (IDDM10), 14q24-q31 (IDDM11), 2q31-q33 (IDDM12), 2q34-q35 (IDDM13), 6q21 (IDDM15), 14q32 (IDDM16), 10q25 (IDDM17), 5q33 (IDDM18), 7p15-p13 (GCK), 1q42, 16q22-q24, Xp11 (conditional on HLA-DR genotype), and 8q22-q24 (8eC10,12eC17). Even though statistical evidence supporting linkage for some of these regions was strong in the initial reports, most regions have not been clearly established in multiple populations (9,23).
A major barrier to type 1 diabetes gene identification, given the likely small locus-specific contribution (low S) for non-HLA genes, is the relatively small number of newly ascertained affected sibpair families with type 1 diabetes. In addition, previous studies have used available samples from a variety of collections, making the compilation of linkage results difficult because of apparent overlap in families analyzed and uncertainty in the equivalence of allele coding. To facilitate the genetic analysis of type 1 diabetes, results from two previous genome scans of type 1 diabetes were merged, comprising 767 families and 831 affected sibpairs from the U.K. and U.S. (24). The combined analyses supported linkage to at least six non-HLA regions to type 1 diabetes, including IDDM2 (INS, nominal P = 6.5 x 10eC4), 2q31-q33 (P = 5.1 x 10eC4), and 10p11 (P = 3.2 x 10eC4). A third genome scan of 424 type 1 diabetic families with at least two affected relative pairs (464 affected pairs) from Scandinavia (25) found no evidence for linkage at these latter three loci but did support linkage on chromosomes 5q11.2 (nominal P = 8.1 x 10eC4) and 16p13 (P = 1.6 x 10eC4). The IDDM15 region on chromosome 6q21 supported linkage (P = 7.0 x 10eC7) when HLA was taken into consideration in a combined analysis of U.S., French, and Scandinavian families.
We present here a joint analysis of data from these three prior genome-wide scans (U.S., U.K., and Scandinavia) as well as 254 new families collected for this study, a total of 1,435 multiplex families, for linkage to type 1 diabetes. With an average map information content of 67% (from 400 polymorphic microsatellite markers in each scan), this family collection provides 95% power to detect a locus with locus-specific S 1.3 and P = 10eC4. In addition to HLA, there was nominal evidence for linkage of type 1 diabetes to 10 other chromosome regions, including 6q21 (IDDM15) and 3 that reached genome-wide levels of significance, 2q31-q33 (IDDM12 and IDDM7), 10p11-q14 (IDDM10), and 16p12-q24. These data support the existence of non-HLA susceptibility loci for type 1 diabetes and strengthen support for a subset of loci previously proposed to contribute to type 1 diabetes risk.
RESEARCH DESIGN AND METHODS
Four sets of Caucasian families provided genome scan data for the combined analyses. Three sets of families have been previously published—U.K. (21,26), U.S. (24,26), and Scandinavia (25)—and one set of 254 families that were newly assembled for this study. The new collection of DNA samples from 254 families was obtained from several sources. DNA samples from 47 U.K. families were identified from the Diabetes U.K. Warren repository (28) that had not been genotyped previously. Families not previously used in published genome scans from the U.S. were contributed by the Joslin Diabetes Center (121 families) and from the Human Biological Data Interchange (76 families). Ten families from Australia were collected by investigators at the Walter and Eliza Hall Institute as previously described (29). In total, there were 1,435 families containing 6,899 individuals (6,358 with genome scan data). A total of 3,109 individuals were affected (with type 1 diabetes), and of these, 3,072 had genotype data (for details of the samples, see table in online appendix [available at http://diabetes.diabetesjournals.org]).
Genotyping.
Microsatellite marker genotyping technologies and allele scoring conventions varied between the different laboratories providing the data for previously published type 1 diabetic families. Details of genotyping of the U.K. (21), U.S. (27), and Scandinavian (25) families have been previously described. The 254 new families were genotyped by the Center for the Inheritance of Disease Research (http://www.cidr.jhmi.edu/) using a panel of 405 microsatellite markers. Because direct merging of genotypes (by standardized allele size) was not possible, within-family recoded genotype data from all four sources were merged into a single database. Family-naming conventions between the samples were normalized, and individual marker names were modified to indicate the laboratory of origin for the genotyping. After elimination of marker inconsistencies (see below), genetic markers were selected to form the analysis panel. Multiple independent genotyping occurred for some markers on a subset of individuals by different laboratories. In these cases, only one marker was used for the current analysis. Unless a marker showed inconsistencies in identity-by-descent (IBD) sharing, the marker that was included for analysis was the one scored in the most samples. A total of 1,190 markers were included in the combined analyses. Markers that previously had been added to maps because of association with type 1 diabetes were excluded to avoid bias in the multipoint linkage results. The excluded markers were located on chromosome 2 (alpha4, ND1, D2S152, CTLA4, and IGFBP5) and chromosome 11 (INS, TH).
Statistical analysis.
Before integration of the genetic data, marker error detection and pedigree structure within each dataset were made using PREST software (30). This method uses the genome scan data to determine the likelihood of each specified relationship given the genetic data. Unlikely marker genotypes were resolved by recoding the specific genotype to "unknown." Occurrences of nonpaternity were resolved by changing the pedigree structure to that which was most likely and then repeating the analysis to confirm appropriate relationships. An integrated marker map was developed by using public databases (Mammalian Genotyping Service, http://research.marshfieldclinic.org/genetics/Genotyping_Service/mgsver2.htm; Southampton, http://cedar.genetics.soton.ac.uk/public_html/; Cooperative Human Linkage Center, http://gai.nci.nih.gov/CHLC/; deCODE, http://www.decode.is) as well as physical map and genome sequence information from the University of California at Santa Cruz (http://genome.ucsc.edu/) using primer sequences in BLAST searches against the genome sequence. All analyses were based on this "consensus" map. Single and multipoint linkage analyses (based on the consensus map order and distances) were performed using GeneHunter-plus [S(pairs) option (31eC33)]. Examination of double recombinants was performed using Merlin software (34). Information content was estimated using Allegro (35).
Estimation of IBD statistics and resulting likelihoods under the null and alternative hypotheses were computed within each dataset. These dataset-specific likelihoods were then combined for the combined linkage analyses. Multipoint linkage analyses were performed, and maximized logarithm of odds (LOD) scores were calculated under an exponential model with constrained between 0 and 2 (32). Genome-wide empirical P values were determined by simulating Mendelian transmission with families maintaining the patterns of missing data observed in the sample (36). The fraction of LOD scores observed greater than the nominal value, based on simulations of 10,000 replicates across the genome, provided the estimated genome-wide P value. Exclusion mapping was performed using the MapMaker/SIBS program (37).
Previous studies have identified a potential type 1 diabetes susceptibility locus near (but not within) the HLA complex (IDDM15) (25,27,38). In the presence of strong support for linkage due to IDDM1, determining the support for IDDM15 is complex because of the positive correlation among the IBD proportions for linked loci. Without accounting for this positive correlation, statistical tests are biased toward inferring an epistatic relationship. A simple approach when IBD is known is to compare the observed correlation, r, between the IBD estimates at two loci (i.e., IDDM1 and IDDM15) with the theoretical correlation, = (1eC2)2, for two loci separated by recombination fraction, . The statistic t = (z eC )neC 3, where z = ln () and = ln () allows for the test of "interaction" that contrasts the observed IBD at IDDM15 based on the distance of IDDM15 from IDDM1 and the expected IBD given that distance. Using a sex-averaged map, IBD estimates for each sibpair in the data were computed, and a single IBD estimate was selected from each pedigree. Thus, each of these IBD estimates is independent. Under the null hypothesis of no interaction between these loci, t approximately has a standard normal distribution. Observed correlations greater than reflect increased sharing at IDDM15 over that expected from IBD at IDDM1 and the hypothesized genetic distance of IDDM15 from IDDM1. A series of 10,000 simulations were performed as described above (36), with the correlation in IBD computed between IDDM1 and IDDM15. The empirical P value for significance of IDDM15 was based on the number of simulated correlations greater than that observed in the original data.
RESULTS
Linkage analysis.
Analysis of 1,190 genetic markers in 1,435 families revealed that the strongest evidence for linkage to type 1 diabetes was on chromosome 6p21 (nominal P = 2.0 x 10eC52) in the MHC (Fig. 1).There were nine noneCHLA-linked regions with nominal evidence supporting linkage to type 1 diabetes (P < 0.01), including 2q31-q33 (IDDM7 and IDDM12; P = 9.0 x 10eC5; genome-wide P = 0.016), 10p14-q11 (IDDM10, P = 1.2 x 10eC4; genome-wide P = 0.021), and 16q22-q24 (P = 4.9 x 10eC4; genome-wide P = 0.075). The estimates of genetic relative risk varied by region, with IDDM1 having the largest locus-specific effect (S 3.3) and the other sites having low locus-specific effects (S 1.1eC1.2). Individual linkage plots by chromosome are provided in the online appendix (supplementary figures, upper curves).
Exclusion mapping.
Each of the four datasets had the equivalent of an 9-centiMorgan (cM) map, providing an average marker information content of 66%. The range of information from the markers was 62% for chromosome 15 to 73% for chromosome 16. In the combined data, 82% of the genome could be excluded at LOD < eC2 for loci of effect size S 1.3 (supplementary figures, lower curves), and >95% could be excluded for S 1.5. Several entire chromosomes could be excluded (chromosomes 7, 8, 18, 20, 21, and 22). The majority of chromosome 6 (due to the strongly linked 6p21/MHC region) and <50% of chromosomes 16, 19, and X could be excluded for S 1.3. For effects of S 1.1, only 6% of the genome could be excluded. The extent of exclusion for each chromosome is shown in Table 1.
IDDM15.
The t-statistic was computed to determine the increase in sharing at IDDM15, beyond that expected from sharing at IDDM1 and expected decay in sharing due to genetic distance. In the current data, the test statistic was computed for IDDM15 at a position 37 cM from IDDM1 (at 47 cM on chromosome 6), using a sex-averaged map. The test statistic was also computed for a range of map positions within 5 cM (32eC42 cM) with similar results. At the 37-cM distance from IDDM1 and IDDM15, the expected correlation coefficient between the IBD estimates under the null hypothesis was estimated as = 0.2276. The observed Pearson’s correlation coefficient between the IBD estimates at IDDM1 and IDDM15, using 1,401 informative pedigrees, was r = 0.3132 (empirical P < 1.0 x 10eC4). These results support an HLA-independent effect in the IDDM15 region.
DISCUSSION
In a previous combined analysis of U.K. and U.S. families (24), it was concluded that an effort to merge and jointly analyze existing families would be required to clarify the role of noneCHLA-linked loci in type 1 diabetes. In the present study, this effort has been achieved under the auspices of the Type 1 Diabetes Genetics Consortium (http://www.t1dgc.org). We have assembled families and merged data from three large genome scans and added new data from 254 families not previously scanned. This increased sample size has allowed the exclusion of >80% of the human genome for locus-specific, but population-independent, effects of S 1.3. In addition to continued support for type 1 diabetes susceptibility related to the MHC (IDDM1) and INS (IDDM2), we identified eight regions that supported noneCHLA-linked susceptibility. Furthermore, we identified three chromosomes that contained extensive areas for which linkage could not be excluded—chromosomes 16, 19, and X—that could benefit from further genotyping to increase information content and the analysis of additional families.
Three noneCHLA-linked regions provided support for linkage at, or near, the genome-wide level of significance (P < 0.05): chromosome 2q31-q33 (IDDM7 and IDDM12), 10p14-q11 (IDDM10), and 16q22-q24. These locations were unlikely to have occurred by chance, based on our simulations. Together with the six other noneCHLA-linked regions that exhibited evidence of linkage at nominal P < 0.01 (Table 2) and the 10-cM map, the data indicate a strong non-HLA genetic effect for type 1 diabetes (39). Furthermore, none of the three most strongly linked regions exhibited support for linkage in the 408 families from Scandinavia (25).
Strong support for linkage (nominal P = 7.0 x 10eC7), after taking into account linkage to HLA) to the IDDM15 locus (6q21), was observed previously in a combined analysis of French, U.S., and Scandinavian families (25,38). In the present study, we have obtained support for IDDM15 (empirical P < 1.0 x 10eC4). The ability to further define the effects of this locus will be facilitated by increased information content in the HLA region and in the region surrounding IDDM15 to better estimate the observed IBD sharing and model the residual linkage to the HLA region, including taking into account sex-specific genetic map differences.
The IDDM12 locus lies within the 2q31-q33 region and has been attributed to single nucleotide polymorphisms (SNPs) in the 3'-untranslated region of CTLA4 (15); however, the modest S value predicted from the odds ratios (ORs) (1.1eC1.2) of the disease-associated SNPs at CTLA4 (S 1.01) in type 1 diabetes seems unlikely to fully account for the magnitude of the observed evidence for linkage (regional S 1.19). This result suggests the presence of other loci in the 2q31-q33 region, if this linkage is confirmed in other future studies. Originally, IDDM7 at chromosome 2q33 was assigned on the basis of evidence of allelic association of the D2S152 microsatellite marker, but this association has not been substantiated. The evidence supporting linkage in the current study does not include the putative IDDM13 locus at chromosome 2q34-q35 (29). Linkage of type 1 diabetes to 10p14-q13 (IDDM10) is well supported by the current and past studies (24,26); however, there has been little follow-up other than association analyses of the functional candidate gene GAD2, suggesting that this gene is not a type 1 diabetes susceptibility locus (40,41). The observed locus-specific effects for 2q31-q33 (S 1.19) and 10p14-q11 (S 1.12) suggest that a single common susceptibility allele would have an allelic association (OR) 3. This effect should be identifiable in a fine-mapping association study using dense SNP maps across the regions.
Support for a type 1 diabetes susceptibility locus on chromosome 16p12-q11.1, which was observed independently in both the combined U.K. and U.S. families (nominal P = 4.5 x 10eC3) and in the Scandinavian families (P = 2 x 10eC4), remained in the present study (P = 3.3 x 10eC3). A recent analysis of four rheumatoid arthritis genome scans (42) reported evidence for linkage at chromosomes 6p21 (HLA; P = 2 x 10eC5) and 16p-cen (P = 0.004). Because rheumatoid arthritis, antithyroid autoimmune disease, and type 1 diabetes cluster in families more often than expected by chance (43), evidence for linkage for any one of these autoimmune diseases could be informative for others. Evidence for linkage in U.K. families with early-onset rheumatoid arthritis (44) to chromosome 16p has previously been demonstrated (P = 3.2 x 10eC4). The comparison between linkage scan results for type 1 diabetes and rheumatoid arthritis provides other interesting similarities. The largest, single, combined scan of rheumatoid arthritis families (45) reported significant linkage of rheumatoid arthritis to chromosome 6p21 (HLA; P = 5 x 10eC12), and some evidence (P < 0.005) of rheumatoid arthritis linked to six other regions (1q43, 6q21, 10q21, 12q12, 17p13, and 18q21). This overlap with potential type 1 diabetes susceptibility at 6q21, 12q12, and 16p-cen may not be coincidental in the etiology of these autoimmune diseases.
Recently, evidence for association of type 1 diabetes with alleles in the PTPN22 locus (chromosome 1p13) has been reported (18), and this association has been confirmed (19,20). PTPN22 encodes a lymphoid-specific tyrosine phosphatase (LYP) and is also associated with autoimmune thyroid disease, rheumatoid arthritis, and SLE (46). The absence of evidence supporting linkage of type 1 diabetes to chromosome 1p13 (D1S206) in the 1,435 families studied here is not surprising given the magnitude of the PTPN22 association with type 1 diabetes. The OR of PTPN22 is large (1.7) but the S is 1.05. Thus, to detect linkage at P < 0.001 with 50% power, a sample of >8,000 affected sibpair families would be required, using a fully informative genetic map. Assuming a multiplicative model, the contribution of PTPN22 to type 1 diabetes (based on the observed OR) is 2%, much lower than HLA (40eC50%). Nevertheless, the knowledge that PTPN22 may be involved in risk to type 1 diabetes, as well as in other autoimmune diseases, is significant and could provide insight into modulating T-cell activity for disease prevention.
Several previously supported regions of linkage have diminished support in the current analyses. Type 1 diabetes susceptibility locus on chromosome 1q42 was strongly supported (nominal P = 9.8 x 10eC5) in a study of 679 U.K. and U.S. families (24,47) but exhibited decreasing linkage support in a follow-up analysis of 616 families using a denser map in the region (P = 4.0 x 10eC4). The evaluation of previously supported regions is difficult, even in the present study with over 1,600 affected sibpairs, particularly for regions with low s (e.g., s 1.1). For example, the region on 1q42 was originally supported with MLS = 3.31 and s 1.5 (27). Although the current sample excludes this region at the reported s 1.5, support for linkage in this region has decreased with increasing sample size from a LOD = 2.20 (24), to the current LOD = 0.87 (nominal P = 1.4 x 10eC2) with s 1.05. At the current estimated magnitude of genetic effect, the 1q42 region could not be excluded.
In a combined analysis of animal model and human linkage data from a number of autoimmune diseases, chromosome 18q12-q21 demonstrated evidence of linkage, which has now been supported by the analysis of congenic strains in NOD mice (48). There was no support for loci on chromosome 18 at P < 0.05 in the current study. Three independent studies of type 1 diabetes have reported linkage to chromosome 8q (49eC51), but there was no support for 8q in this study. Additional previously reported loci with relatively little support for linkage in the current study include IDDM4 (11q13), IDDM6 (18q12-q21), IDDM9 (3q22-q25), IDDM11 (14q24-q31), IDDM16 (14q32), IDDM17 (10q25), and IDDM18 (5q33). These data suggest that these putative type 1 diabetes susceptibility loci represent either false-positive results or have very small effects that may be more readily detected in certain populations because of variation in allele frequencies or other factors, including the possibility of population-specific genetic or environmental effects.
Linkage and fine mapping studies in mouse models of SLE (52) and of type 1 diabetes (53) have demonstrated that a single linkage peak may be composed of several susceptibility loci. In human populations, a linkage signal may be observed by the chance clustering of several disease loci, each with relatively weak locus-specific effects. The presence of multiple susceptibility loci may also account, in part, for broad linkage peaks often observed in studies of complex, common diseases. This underlying complexity would also increase the difficulty to obtain convincing results in future fine-mapping association studies. Both development of novel analytical approaches and increased sample size will be necessary to resolve this apparent complexity (54,55). Through the efforts of consortia (such as our international effort), it will be possible to increase the number of families for type 1 diabetes (http://www.t1dgc.org), which would increase power and allow exclusion of loci s 1.2, as well as provide standardized samples and reagents for future fine-mapping studies.
These results suggest two parallel tracks for the identification of type 1 diabetes susceptibility loci. First, systematic fine mapping of all the variants responsible for the HLA linkage to type 1 diabetes is justified, especially within the 4-Mb HLA region. Second, further exploration of potential non-HLA regions described here is now justified, including chromosomes 2q31-q33, 6q21, 10p14-q11, and 16q22-q24. Because sample sizes in linkage and association studies have historically been small and few genes (out of 25,000) have been studied in depth, future collaborative efforts and establishment of accessible resources for study should increase the yield of true disease susceptibility loci.
Appendix
Members of the Type 1 Diabetes Genetics Consortium
Asia-Pacific Network: Eesh Bhatia, Francois Bonnici, Pik To Cheung, Peter Colman, Andrew Cotterill, Jenny Couper, Ric Cutfield, Elizabeth Davis, Tim Davis, Paul Dixon, Kim Donaghue, Paul Drury, Mark Harris, Len Harrison, Tim Jones, Uma Kanga, Alok Kanungo, Betty Kek, Jeremy Krebs, Yann-Jinn Lee, Margaret Lloyd, Amanda Loth, Narinder Mehra, Grant Morahan, C.B. Sanjeevi, Brian Tait, Mike Varney, Jinny Willis, Loke Kah Yin, Goh Siok Ying.
European Network: Francisco J. Ampudia-Blasco, Jesus Argente, Magdalena Avbelj, Gulja Babadjanova, Klaus Badenhoop, Lubomir Barak, Christos Bartsocas, Tadej Battelino, Emilia Belda, Polly Bingley, Bernhard O. Boehm, Eulalia Brugues, Raffaella Buzzetti, Joyce Carlson, Luis Castano, Anna Casu, Ondrej Cinek, Raquel Corripio, Alberto de Leiva, V.A. Ruiz Esquide, Ana Fagulha, Marta Hernandez Garcia, Cristian Guja, Alona Hamou, Erifili Hatziagelaki, Simon Heath, Kaire Heilman, Nora Hosszufalusi, Constantin Ionescu-Tirgoviste, Cecile Julier, Ida Kinalska, Ingrid Kockum, Kalinka Koprivarova, Adam Kretowski, Dora Krikovszky, Nebojsa Lalic, Nicole Lambracht, Merce Lara, Mark Lathrop, Katharina Laubner, Claire Levy-Marchal, Johnny Ludvigsson, Laszlo Madacsy, Mara Marga, Didac Mauricio, Jorn Nerup, Antanas Norkus, Anna Okruszko, Xavier Palomer, Teresa Pedro, Moshe Phillip, Valdis Pirags, Flemming Pociot, Galina Popova, Paolo Pozzilli, Bart O. Roep, Silke Rosinger, Ana Maria Varela Sande, Ilhan Satman, Edith Schober, Jochen Seufert, Jan Skrha, Gyula Soltesz, Giatgen Spinas, Juraj Stanik, Tanya Szendeffy, Vallo Tillmann, Dag Undlien, Vaidotas Urbanavicius, Luciana Valente, Bart Van der Auwera, Andriani Vazeou-Gerasimidi, Dzilda Velickiene, Ana Wagner, Markus Walter, Alistair Williams, Miroslav Wurzburger, Anette-G Ziegler.
North American Network: Alan Aldrich, Mark Anderson, Christophe Benoist, Noureddine Berka, Patrick Concannon, Mark Daly, Jayne Danska, Larry Dolan, David Donaldson, Alessandro Doria, Janice Dorman, George Eisenbarth, Henry Erlich, Pamela Fain, Rosanna Fiallo-Scharer, Kenneth Gabbay, Daniel Geraghty, Soumitra Ghosh, Steven Gitelman, Nat Goodman, Gregory Goodwin, Carla Greenbaum, William Hagopian, John Hansen, Joel Hirschhorn, Leroy Hood, Kevin Kaiserman, Jean Lawrence, Victoria Magnuson, Jennifer Marks, John Mayberry, Elizabeth Mayer-Davis, Richard McIndoe, Brad McNeney, Eric Mickelson, Antoinette Moran, Gerald Nepom, Janelle Noble, Jill Norris, Tihamer Orban, David Owerbach, Andrew Paterson, Catherine Pihoker, Constantin Polychronakos, Alberto Pugliese, Philip Raskin, Marian Rewers, Henry Rodriquez, Jerry Rotter, Monique Roy, Desmond Schatz, Gary Schoch, Jin-Xiong She, Richard Spielman, Andrea Steck, Kent Taylor, Jay Tischfield, Ellen Toth, Diane Wherrett, Stephen Willi, Darrell Wilson, Lue Ping Zhao.
U.K. Network: Francesco Cucca, David Dunger, Simon Howell, Sarah Nutland, Helen Rance, Luc Smink, John Todd, Neil Walker, Barry Widmer, Heather Withers.
Coordinating Center: Don Babcock, Stephanie Beck, Mark Brown, Cralen Davis, Mark Espeland, Clara Gorodezky, Jason Griffin, Mark Hall, Teresa Harnish, Laura Hemrick, John Hepler, Joan Hilner, Letitia Howard, Ethan Lange, Carl Langefeld, Josyf Mychaleckyj, June Pierce, David Reboussin, Stephen Rich, Scott Rushing, Michele Sale, Elizabeth Sides, Beverly Snively, Michael Steffes, Augy Thiel, Lynne Wagenknecht, Dustin Williams, Jianzhao Xu, Belinda Youngdahl.
ACKNOWLEDGMENTS
P.C. has received support from the National Institute of Diabetes and Digestive and Kidney Diseases and Juvenile Diabetes Research Foundation. S.S.R. has received support from the National Institute of Diabetes and Digestive and Kidney Diseases and the JDRF. J.A.T. has received support from the JDRF and the Wellcome Trust.
We acknowledge the assistance of Mark Brown, Mathew Barber, Heather Cordell, Neil Walker, Carl Langefeld and Nancy Cox.
FOOTNOTES
A complete list of the members of the Type 1 Diabetes Genetics Consortium can be found in the APPENDIX.
Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org.
IBD, identity-by-descent; LOD, logarithm of odds; MHC, major histocompatibility complex; SNP, single nucleotide polymorphism
REFERENCES
Rewers M, LaPorte RE, King H, Tuomilehto J: Trends in the prevalence and incidence of diabetes: insulin-dependent diabetes mellitus in childhood. World Health Stat Q41 :179 eC189,1988
Rewers MJ: The changing face of the epidemiology of insulin-dependent diabetes mellitus (IDDM): research designs and models of disease causation. Ann Med23 :419 eC426,1991
Genuth S, Alberti KG, Bennett P, Buse J, Defronzo R, Kahn R, Kitzmiller J, Knowler WC, Lebovitz H, Lernmark A, Nathan D, Palmer J, Rizza R, Saudek C, Shaw J, Steffes M, Stern M, Tuomilehto J, Zimmet P; Expert Committee on the Diagnosis and Classification of Diabetes Mellitus: Follow-up report on the diagnosis of diabetes mellitus. Diabetes Care31 :3160 eC3167,2003
Hirshhorn JN: Genetic epidemiology of type 1 diabetes. Pediatr Diabetes4 :87 eC100,2003
Noble JA, Valdes AM, Cook M, Klitz W, Thomson G, Erlich HA: The role of HLA class II genes in insulin-dependent diabetes mellitus: molecular analysis of 180 Caucasian, multiplex families. Am J Hum Genet59 :1134 eC1148,1996
Hyttinen V, Kaprio J, Kinnunen L, Koskenvuo M, Tuomilehto J: Genetic liability of type 1 diabetes and the onset age among 22,650 young Finnish twin pairs: a nationwide follow-up study. Diabetes52 :1052 eC1055,2003
Svensson J, Carstensen B, Molbak A, Christau B, Mortensen HB, Nerup J, Borch-Johnsen K: Increased risk of childhood type 1 diabetes in children born after 1985. Diabetes Care25 :2197 eC2201,2002
Onkamo P, Vaananen S, Karvonen M, Tuomilehto J: Worldwide increase in incidence of type 1 diabetes: the analysis of the data on published incidence trends. Diabetologia42 :1395 eC1403,1999
Pociot F, McDermott MF: Genetics of type 1 diabetes mellitus (Review). Genes Immun3 :235 eC249,2002
Ronningen KS, Keiding N, Green A, EURODIAB ACE Study Group: Europe and diabetes: correlations between the incidence of childhood-onset type I diabetes in Europe and HLA genotypes. Diabetologia44 (Suppl. 3) :B51 eCB59,2001
Rich SS: Mapping genes in diabetes: a genetic epidemiological perspective. Diabetes39 :1315 eC1319,1990
Valdes AM, Noble JA, Genin E, Clerget-Darpoux F, Erlich HA, Thomson G: Modeling of HLA class II susceptibility to type 1 diabetes reveals an effect associated with DPB1. Genet Epidemiol21 :212 eC223,2001
Bell GI, Horita S, Karam JH: A polymorphic locus near the human insulin gene is associated with insulin-dependent diabetes mellitus. Diabetes33 :176 eC183,1984
Julier C, Hyer RN, Daviews J, Merlin F, Soularue P, Briant L, Cathelineau G, Deschamps I, Rotter JI, Froguel P: Insulin-IGF2 region on chromosome 11p encodes a gene implicated in HLA-DR4-dependent diabetes susceptibility. Nature354 :155 eC159,1991
Julier C, Lucassen A, Villedieu P, Delepine M, Levy-Marchal C, Danze PM, Bianchi F, Boitard C, Froguel P, Bell J: Multiple DNA variant association analysis: application to the insulin gene region in type I diabetes. Am J Hum Genet55 :1247 eC1254,1994
Barratt BJ, Payne F, Lowe CE, Hermann R, Healy BC, Harold D, Concannon P, Gharani N, McCarthy MI, Olavesen MG, McCormack R, Guja C, Ionescu-Tirgoviste C, Undlien DE, Ronningen KS, Gillespie KM, Tuomilehto-Wolf E, Tuomilehto J, Bennett ST, Clayton DG, Cordell HJ, Todd JA: Remapping the insulin gene/IDDM2 locus in type 1 diabetes. Diabetes53 :1884 eC1889,2004
Ueda H, Howson JM, Esposito L, Heward J, Snook H, Chamberlain G, Rainbow DB, Hunter KM, Smith AN, Di Genova G, Herr MH, Dahlman I, Payne F, Smyth D, Lowe C, Twells RC, Howlett S, Healy B, Nutland S, Rance HE, Everett V, Smink LJ, Lam AC, Cordell HJ, Walker NM, Bordin C, Hulme J, Motzo C, Cucca F, Hess JF, Metzker ML, Rogers J, Gregory S, Allahabadia A, Nithiyananthan R, Tuomilehto-Wolf E, Tuomilehto J, Bingley P, Gillespie KM, Undlien DE, Ronningen KS, Guja C, Ionescu-Tirgoviste C, Savage DA, Maxwell AP, Carson DJ, Patterson CC, Franklyn JA, Clayton DG, Peterson LB, Wicker LS, Todd JA, Gough SC: Association of the T-cell regulatory gene CTLA4 with susceptibility to autoimmune disease. Nature423 :506 eC511,2003
Bottini N, Musumeci L, Alonso A, Rahmouni S, Nika K, Rostamkhani M, MacMurray J, Meloni GF, Lucarelli P, Pellecchia M, Eisenbarth GS, Comings D, Mustelin T: A functional variant of lymphoid tyrosine phosphatase is associated with type I diabetes. Nat Genet36 :337 eC338,2004
Smyth DJ, Cooper JD, Collins JE, Heward JM, Franklyn JA, Howson JMM, Walker N, Vella A, Nutland S, Rance HE, Maier LM, Barratt BJ, Guja C, Ionescu-Tirgoviste CA, Savage D, Dunger DB, Widmer B, Strachan DP, Ring SM, Clayton DG, Twells RCJ, Gough SCL, Todd JA: 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. Diabetes53 :3020 eC3023,2004
Onengut-Gumuscu S, Ewens KG, Spielman RS, Concannon P: A functional polymorphism (1858C/T) in the PTPN22 gene is linked and associated with type 1 diabetes in multiplex families. Genes Immun5 :678 eC680,2004
Davies JL, Kawaguchi Y, Bennett ST, Copeman JB, Cordell HJ, Pritchard LE, Reed PW, Gough SCL, Jenkins SC, Palmer SM, Balfour KM, Row B, Farrall M, Barnett AH, Bain SC, Todd JA: A genome-wide search for human type 1 diabetes susceptibility genes. Nature371 :130 eC136,1994
Hashimoto L, Habita C, Beressi J, Delepine M, Besse C, Cambon-Thomsen A, Deschamps I, Rotter J, Djoulah S, James M, Froguel P, Weissenbach J, Lathrop GM, Julier C: Genetic mapping of a susceptibility locus for insulin-dependent diabetes mellitus on chromosome 11q. Nature371 :161 eC164,1994
Melanitou E, Fain P, Eisenbarth GS: Genetics of type 1A (immune mediated) diabetes. J Autoimmun21 :93 eC98,2003
Cox NJ, Wapelhorst B, Morrison VA, Johnson L, Pinchuk L, Spielman RS, Todd JA, Concannon P: Seven regions of the genome show evidence of linkage to type 1 diabetes in a consensus analysis of 767 multiplex families. Am J Hum Genet69 :820 eC830,2001
Mein CA, Esposito L, Dunn MG, Johnson GC, Timms AE, Goy JV, Smith AN, Sebag-Montefiore L, Merriman ME, Wilson AJ, Pritchard LE, Cucca F, Barnett AH, Bain SC, Todd JA: A search for type 1 diabetes susceptibility genes in families from the United Kingdom. Nat Genet19 :297 eC300,1998
Concannon P, Gogolin-Ewens KJ, Hinds D, Wapelhorst B, Morrison VA, Stirling B, Mitra M, Farmer J, Williams SR, Cox NJ, Bell GI, Risch N, Spielman RS: A second-generation screen of the human genome for susceptibility to insulin-dependent diabetes mellitus (IDDM). Nat Genet19 :292 eC296,1998
European Consortium for IDDM Genome Studies: A genomewide scan for type 1-diabetes susceptibility in Scandinavian families: identification of new loci with evidence of interactions. Am J Hum Genet69 :1301 eC1313,2001
Bain SC, Todd JA, Barnett AH: The British Diabetic Association: Warren repository. Autoimmunity7 :83 eC85,1990
Morahan G, Huang D, Tait B, Colman PG, Harrison LC: Markers on distal chromosome 2q linked to insulin-dependent diabetes mellitus. Science272 :1811 eC1813,1996
McPeek MS, Sun L: Statistical tests for detection of misspecified relationships by use of genome screen data. Am J Hum Genet66 :1076 eC1094,2000
Kruglyak L, Daly MJ, Reeve-Daly MP, Lander ES: Parametric and non-parametric linkage analysis: a unified multipoint approach. Am J Hum Genet58 :1347 eC1363,1996
Kong A, Cox NJ: Allele-sharing models: LOD scores and accurate linkage tests. Am J Hum Genet61 :1179 eC1188,1997
Langefeld CD, Boehnke M: Multiple trait nonparametric linkage regression. Am J Hum Genet65 :A258 ,1999
Abecasis GR, Cherny SS, Cookson WOC, Cardon LR: MERLIN: rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet30 :97 eC101,2002
Gudbjartsson DF, Jonasson K, Frigge M, Kong A: Allegro, a new computer program for multipoint linkage analysis. Nat Genet25 :12 eC13,2000
Hauser ER, Boehnke M: Genetic linkage analysis of complex genetic traits by using affected sibling pairs. Biometrics54 :1238 eC1246,1998
Kruglyak L, Lander ES: Complete multipoint sib-pair analysis of qualitative and quantitative traits. Am J Hum Genet57 :439 eC454,1995
Delepine M, Pociot F, Habita C, Hashimoto L, Froguel P, Rotter J, Combon-Thomsen A, Deschamps I, Djoulah S, Weissenbach J, Nerup J, Lathrop M, Julier C: Evidence of a non-MHC susceptibility locus in type 1 diabetes linked to HLA on chromosome 6. Am J Hum Genet60 :174 eC187,1997
Cucca F, Goy JV, Kawaguchi Y, Esposito L, Merriman ME, Wilson AJ, Cordell HJ, Bain SC, Todd JA: A male-female bias in type 1 diabetes and linkage to chromosome Xp in MHC HLA-DR3-positive patients. Nat Genet3 :301 eC302,1998
Wiltshire S, Cardon LR, McCarthy MI: Evaluating the results of genomewide linkage scans of complex traits by locus counting. Am J Hum Genet71 :1175 eC1182,2002
Wapelhorst B, Bell GI, Risch N, Spielman RS, Concannon P: Linkage and association studies in insulin-dependent diabetes with a new dinucleotide repeat polymorphism at the GAD65 locus. Autoimmunity21 :127 eC130,1995
Johnson GC, Payne F, Nutland S, Stevens H, Tuomilehto-Wolf E, Tuomilehto J, Todd JA: A comprehensive, statistically powered analysis of GAD2 in type 1 diabetes. Diabetes51 :2866 eC2870,2002
Fisher SA, Lanchbury JS, Lewis CM: Meta-analysis of four rheumatoid arthritis genome-wide linkage studies: confirmation of a susceptibility locus on chromosome 16. Arthritis Rheum48 :1200 eC1206,2003
Tait KF, Marshall T, Berman J, Carr-Smith J, Rowe B, Todd JA, Bain SC, Barnett AH, Gough SC: Clustering of autoimmune disease in parents of siblings from the type 1 diabetes Warren repository. Diabet Med21 :358 eC362,2004
MacKay K, Eyre S, Myerscough A, Millicic A, Barton A, Laval S, Barrett J, Lee D, While S, Jhon S, Brown MA, Bell J, Silman A, Ollier W, Wordsworth P, Worthington J: Whole-genome linkage analysis of rheumatoid arthritis susceptibility loci in 252 affected sibling pairs in the United Kingdom. Arthritis Rheum46 :632 eC639,2002
Jawaheer D, Seldin MF, Amos CI, Chen WV, Shigeta R, Etzel C, Damle A, Xiao X, Chen D, Lum RF, Monteiro J, Kern M, Criswell LA, Albani S, Nelson JL, Clegg DO, Pope R, Schroeder HW Jr, Bridges SL Jr, Pisetsky DS, Ward R, Kastner DL, Wilder RL, Pincus T, Callahan LF, Flemming D, Wener MH, Gregersen PK, North American Rheumatoid Arthritis Consortium: Screening the genome for rheumatoid arthritis susceptibility genes: a replication study and combined analysis of 512 multicase families. Arthritis Rheum48 :906 eC916,2003
Begovich AB, Carlton VE, Honigberg LA, Schrodi SJ, Chokkalingam AP, Alexander HC, Ardlie KG, Huang Q, Smith AM, Spoerke JM, Conn MT, Chang M, Chang SY, Saiki RK, Catanese JJ, Leong DU, Garcia VE, McAllister LB, Jeffery DA, Lee AT, Batliwalla F, Remmers E, Criswell LA, Seldin MF, Kastner DL, Amos CI, Sninsky JJ, Gregersen PK: A missense single-nucleotide polymorphism in a gene encoding a protein pyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am J Hum Genet75 :330 eC337,2004
Ewens KG, Johnson LN, Wapelhorst B, O’Brien K, Gutin S, Morrison VA, Street C, Gregory SG, Spielman RS, Concannon P: Linkage and association with type 1 diabetes on chromosome 1q42. Diabetes51 :3318 eC3325,2002
Merriman T, Twells R, Merriman M, Eaves I, Cox R, Cucca F, McKinney P, Shield J, Baum D, Bosi E, Pozzilli P, Nistico L, Buzzetti R, Joner G, Ronningen K, Thorsby E, Undlien D, Pociot F, Nerup J, Bain S, Barnett A, Todd J: Evidence by allelic association-dependent methods for a type 1 diabetes polygene (IDDM6) on chromosome 18q21. Hum Mol Genet6 :1003 eC1010,1997
Luo DF, Buzzetti R, Rotter JI, Maclaren NK, Raffel LJ, Nistico L, Govannini C, Pozzilli P, Thomson G, She JX: Confirmation of three susceptibility genes to insulin-dependent diabetes mellitus: IDDM4, IDDM5 and IDDM8. Hum Mol Genet5 :693 eC698,1996
Sale MM, FitzGerald LM, Charlesworth JC, Bowden DW, Rich SS: Evidence for a novel type 1 diabetes susceptibility locus on chromosome 8. Diabetes51 (Suppl. 3) :S316 eCS319,2002
Vaessen N, Heutink P, Houwing-Duistermaat JJ, Snijders PJ, Rademaker T, Testers L, Batstra MR, Sandkuijl LA, van Duijn CM, Oostra BA: A genome-wide search for linkage disequilibrium with type 1 diabetes in a recent genetically isolated population from the Netherlands. Diabetes51 :856 eC859,2002
Nguyen C, Limaye N, Wakeland EK: Susceptibility genes in the pathogenesis of murine lupus. Arthritis Res4 (Suppl. 3) :S255 eCS263,2002
Todd JA, Wicker LS: Genetic protection from the inflammatory disease type 1 diabetes in humans and animal models. Immunity15 :387 eC395,2001
Rich SS, Concannon P: Challenges and strategies for investigating the genetic complexity of common human diseases. Diabetes51 (Suppl. 3) :S288 eCS294,2002
Pociot F, Karlsen AE, Pedersen CB, Aalund M, Nerup J for the European Consortium for IDDM Genome Studies: Novel analytical methods applied to type 1 diabetes genome-scan data. Am J Hum Genet74 :647 eC660,2004(Patrick Concannon, Henry )