Very high resolution single pass HLA genotyping using amplicon sequencing on the 454 next generation DNA sequencers: Comparison with Sanger sequencing.

Compared to Sanger sequencing, next-generation sequencing offers advantages for high resolution HLA genotyping including increased throughput, lower cost, and reduced genotype ambiguity. Here we describe an enhancement of the Roche 454 GS GType HLA genotyping assay to provide very high resolution (VHR) typing, by the addition of 8 primer pairs to the original 14, to genotype 11 HLA loci. These additional amplicons help resolve common and well-documented alleles and exclude commonly found null alleles in genotype ambiguity strings. Simplification of workflow to reduce the initial preparation effort using early pooling of amplicons or the Fluidigm Access Array™ is also described. Performance of the VHR assay was evaluated on 28 well characterized cell lines using Conexio Assign MPS software which uses genomic, rather than cDNA, reference sequence. Concordance was 98.4%; 1.6% had no genotype assignment. Of concordant calls, 53% were unambiguous. To further assess the assay, 59 clinical samples were genotyped and results compared to unambiguous allele assignments obtained by prior sequence-based typing supplemented with SSO and/or SSP. Concordance was 98.7% with 58.2% as unambiguous calls; 1.3% could not be assigned. Our results show that the amplicon-based VHR assay is robust and can replace current Sanger methodology. Together with software enhancements, it has the potential to provide even higher resolution HLA typing.

A study of genotyping for management of human papillomavirus-positive, cytology-negative cervical screening results.

The effective management of women with human papillomavirus (HPV)-positive, cytology-negative results is critical to the introduction of HPV testing into cervical screening. HPV typing has been recommended for colposcopy triage, but it is not clear which combinations of high-risk HPV types provide clinically useful information. This study included 18,810 women with Hybrid Capture 2 (HC2)-positive, cytology-negative results and who were age ≥30 years from Kaiser Permanente Northern California. The median follow-up was 475 days (interquartile range [IQR], 0 to 1,077 days; maximum, 2,217 days). The baseline specimens from 482 cases of cervical intraepithelial neoplasia grade 3 or cancer (CIN3+) and 3,517 random HC2-positive noncases were genotyped using 2 PCR-based methods. Using the case-control sampling fractions, the 3-year cumulative risks of CIN3+ were calculated for each individual high-risk HPV type. The 3-year cumulative risk of CIN3+ among all women with HC2-positive, cytology-negative results was 4.6%. HPV16 status conferred the greatest type-specific risk stratification; women with HC2-positive/HPV16-positive results had a 10.6% risk of CIN3+, while women with HC-2 positive/HPV16-negative results had a much lower risk of 2.4%. The next most informative HPV types and their risks in HPV-positive women were HPV33 (5.9%) and HPV18 (5.9%). With regard to the etiologic fraction, 20 of 71 cases of cervical adenocarcinoma in situ (AIS) and adenocarcinoma in the cohort were positive for HPV18. HPV16 genotyping provides risk stratification useful for guiding clinical management; the risk among HPV16-positive women clearly exceeds the U.S. consensus risk threshold for immediate colposcopy referral. HPV18 is of particular interest because of its association with difficult-to-detect glandular lesions. There is a less clear clinical value of distinguishing the other high-risk HPV types.

Next-generation sequencing can reveal in vitro-generated PCR crossover products: some artifactual sequences correspond to HLA alleles in the IMGT/HLA database.

The high-resolution human leukocyte antigen (HLA) genotyping assay that we developed using 454 sequencing and Conexio software uses generic polymerase chain reaction (PCR) primers for DRB exon 2. Occasionally, we observed low abundance DRB amplicon sequences that resulted from in vitro PCR 'crossing over' between DRB1 and DRB3/4/5. These hybrid sequences, revealed by the clonal sequencing property of the 454 system, were generally observed at a read depth of 5%-10% of the true alleles. They usually contained at least one mismatch with the IMGT/HLA database, and consequently, were easily recognizable and did not cause a problem for HLA genotyping. Sometimes, however, these artifactual sequences matched a rare allele and the automatic genotype assignment was incorrect. These observations raised two issues: (1) could PCR conditions be modified to reduce such artifacts? and (2) could some of the rare alleles listed in the IMGT/HLA database be artifacts rather than true alleles? Because PCR crossing over occurs during late cycles of PCR, we compared DRB genotypes resulting from 28 and (our standard) 35 cycles of PCR. For all 21 cell line DNAs amplified for 35 cycles, crossover products were detected. In 33% of the cases, these hybrid sequences corresponded to named alleles. With amplification for only 28 cycles, these artifactual sequences were not detectable. To investigate whether some rare alleles in the IMGT/HLA database might be due to PCR artifacts, we analyzed four samples obtained from the investigators who submitted the sequences. In three cases, the sequences were generated from true alleles. In one case, our 454 sequencing revealed an error in the previously submitted sequence.

High throughput HLA genotyping using 454 sequencing and the Fluidigm Access Array™ System for simplified amplicon library preparation.

The human leukocyte antigen (HLA) class I and class II loci are the most polymorphic genes in the human genome; distinguishing the thousands of HLA alleles is challenging. Next generation sequencing of exonic amplicons with the 454 genome sequence (GS) FLX System and Conexio Assign ATF 454 software provides high resolution, high throughput HLA genotyping for eight class I and class II loci. HLA typing of potential donors for unrelated bone marrow donor registries typically uses a subset of these loci at high sample throughput and low cost per sample. The Fluidigm Access Array System enables the incorporation of 48 different multiplex identifiers (MIDs) corresponding to 48 genomic DNA samples with up to 48 different primer pairs in a microfluidic device generating 2304 parallel polymerase chain reactions (PCRs). Minimal volumes of reagents are used. During genomic PCR, in this 4-primer system, the outer set of primers containing the MID and the 454 adaptor sequences are incorporated into an amplicon generated by the inner HLA target-specific primers each containing a common sequence tag at the 5' end of the forward and reverse primers. Pools of the resulting amplicons are used for emulsion PCR and clonal sequencing on the 454 Life Sciences GS FLX System, followed by genotyping with Conexio software. We have genotyped 192 samples with 100% concordance to known genotypes using 8 primer pairs (covering exons 2 and 3 of HLA-A, B and C, and exon 2 of DRB1, 3/4/5 and DQB1) and 96 MIDs in a single GS FLX run. An average of 166 reads per amplicon was obtained. We have also genotyped 96 samples at high resolution (14 primer pairs covering exons 2, 3, and 4 of the class I loci and exons 2 of DRB1, 3/4/5, DQA1, DQB1, DPB1, and exon 3 of DQB1), recovering an average of 173 sequence reads per amplicon.

16(th) IHIW : review of HLA typing by NGS.

Human leucocyte antigen (HLA) genes play an important role in the success of organ transplantation and are associated with autoimmune and infectious diseases. Current DNA-based genotyping methods, including Sanger sequence-based typing (SSBT), have identified a high degree of polymorphism. This level of polymorphism makes high-resolution HLA genotyping challenging, resulting in ambiguous typing results due to an inability to resolve phase and/or defining polymorphisms lying outside the region amplified. Next-generation sequencing (NGS) may resolve the issue through the combination of clonal amplification, which provides phase information, and the ability to sequence larger regions of genes, including introns, without the additional effort or cost associated with current methods. The NGS HLA sequencing project of the 16IHIW aimed to discuss the different approaches to (i) template preparation including short- and long-range PCR amplicons, exome capture and whole genome; (ii) sequencing platforms, including GS 454 FLX, Ion Torrent PGM, Illumina MiSeq/HiSeq and Pacific Biosciences SMRT; (iii) data analysis, specifically allele-calling software. The pilot studies presented at the workshop demonstrated that although individual sequencers have very different performance characteristics, all produced sequence data suitable for the resolution of HLA genotyping ambiguities. The developments presented at this workshop clearly highlight the potential benefits of NGS in the HLA laboratory.

Investigation of the vitamin D receptor gene (VDR) and its interaction with protein tyrosine phosphatase, non-receptor type 2 gene (PTPN2) on risk of islet autoimmunity and type 1 diabetes: the Diabetes Autoimmunity Study in the Young (DAISY).

The present study investigated the association between variants in the vitamin D receptor gene (VDR) and protein tyrosine phosphatase, non-receptor type 2 gene (PTPN2), as well as an interaction between VDR and PTPN2 and the risk of islet autoimmunity (IA) and progression to type 1 diabetes (T1D). The Diabetes Autoimmunity Study in the Young (DAISY) has followed children at increased risk of T1D since 1993. Of the 1692 DAISY children genotyped for VDR rs1544410, VDR rs2228570, VDR rs11568820, PTPN2 rs1893217, and PTPN2 rs478582, 111 developed IA, defined as positivity for GAD, insulin or IA-2 autoantibodies on 2 or more consecutive visits, and 38 IA positive children progressed to T1D. Proportional hazards regression analyses were conducted. There was no association between IA development and any of the gene variants, nor was there evidence of a VDR*PTPN2 interaction. Progression to T1D in IA positive children was associated with the VDR rs2228570 GG genotype (HR: 0.49, 95% CI: 0.26-0.92) and there was an interaction between VDR rs1544410 and PTPN2 rs1893217 (p(interaction)=0.02). In children with the PTPN2 rs1893217 AA genotype, the VDR rs1544410 AA/AG genotype was associated with a decreased risk of T1D (HR: 0.24, 95% CI: 0.11-0.53, p=0.0004), while in children with the PTPN2 rs1893217 GG/GA genotype, the VDR rs1544410 AA/AG genotype was not associated with T1D (HR: 1.32, 95% CI: 0.43-4.06, p=0.62). These findings should be replicated in larger cohorts for confirmation. The interaction between VDR and PTPN2 polymorphisms in the risk of progression to T1D offers insight concerning the role of vitamin D in the etiology of T1D.

Use of class I and class II HLA loci for predicting age at onset of type 1 diabetes in multiple populations.

AIMS/HYPOTHESIS: The study aimed to assess, in multiple populations, the role of HLA alleles on early and late age at onset of type 1 diabetes. METHODS: Stepwise linear regression models were used to determine which HLA class I and class II risk alleles to include. High-resolution genotyping data for patients from the Type 1 Diabetes Genetics Consortium (T1DGC) collection (n = 2,278) and four independent cohorts from Denmark, Sardinia and the USA (Human Biological Data Interchange [HBDI] and Joslin Diabetes Center) (n = 1,324) (total n = 3,602) were used to assess the role of HLA variation on age of onset and predict early onset (age ≤ 5 years) and late onset (age ≥ 15 years) of type 1 diabetes. RESULTS: In addition to carriage of HLA class I alleles A*24:02, B*39:06, B*44:03 and B*18:01, HLA class II DRB1-DQB1 loci significantly contributed to age at onset, explaining 3.4% of its variance in the combined data. HLA genotypes, together with sex, were able to predict late onset in all cohorts studied, with AUC values ranging from 0.58 to 0.63. Similar AUC values (0.59-0.70) were obtained for early onset for most cohorts, except in the Sardinian study, in which none of the models tested had significant predictive power. CONCLUSIONS/INTERPRETATION: HLA associations with age of onset are consistent across most white populations and HLA information can predict some of the risk of early and late onset of type 1 diabetes. Considerable heterogeneity was observed between Sardinian and other populations, particularly with regard to early age of onset.

HLA DNA typing: past, present, and future.

The human leukocyte antigen (HLA) class I and class II loci are the most polymorphic genes in the human genome, with a highly clustered and patchwork pattern of sequence motifs. In the three decades since the first HLA gene was isolated by molecular cloning (a cDNA clone of HLA-B7), thousands of alleles have been identified and the names and sequences of all known alleles have been curated in the IMGT/HLA database. This extensive allelic diversity made and continues to make high-resolution HLA DNA typing very challenging. The first attempt at HLA DNA typing involved restriction fragment length polymorphism (RFLP) analysis, but this approach had many limitations. The development of PCR in 1985 allowed for the amplification of the polymorphic exons of the HLA class I and class II genes and the analysis of the polymorphic sequence motifs with sequence-specific oligonucleotide (SSO) hybridization probes. The immobilization of these probes on membranes and later on beads, along with primer sets for sequence-specific priming (SSP), gave rise to the current set of HLA typing reagents. Sanger sequencing has provided high-resolution typing but, in many cases, genotyping 'ambiguity' remains an issue. In the past few years, the introduction of next-generation sequencing, with the critical properties of massively parallel and clonal sequencing, has significantly reduced HLA genotyping ambiguity. Here, our lab's efforts to develop high-resolution and high-throughput HLA DNA typing using the 454 Sequencing System are reviewed, and the potential future developments and applications of HLA DNA typing are discussed.

Multi-locus HLA class I and II allele and haplotype associations with follicular lymphoma.

Follicular lymphoma (FL) is an indolent, sometimes, fatal disease characterized by recurrence at progressively shorter intervals and is frequently refractive to therapy. Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) in the human leukocyte antigen (HLA) region on chromosome 6p21.32-33 that are statistically significantly associated with FL risk. Low to medium resolution typing of single or multiple HLA genes has provided an incomplete picture of the total genetic risk imparted by this highly variable region. To gain further insight into the role of HLA alleles in lymphomagenesis and to investigate the independence of validated SNPs and HLA alleles with FL risk, high-resolution HLA typing was conducted using next-generation sequencing in 222 non-Hispanic White FL cases and 220 matched controls from a larger San Francisco Bay Area population-based case-control study of lymphoma. A novel protective association was found between the DPB1*03:01 allele and FL risk [odds ratio (OR) = 0.39, 95% confidence interval (CI) = 0.21-0.68]. Extended haplotypes DRB1*01:01-DQA1*01:01-DQB1*05:01 (OR = 2.01, 95% CI = 1.22-3.38) and DRB1*15-DQA1*01-DQB1*06 (OR = 0.55, 95% CI = 0.36-0.82) also influenced FL risk. Moreover, DRB1*15-DQA1*01-DQB1*06 was highly correlated with an established FL risk locus, rs2647012. These results provide further insight into the critical roles of HLA alleles and SNPs in FL pathogenesis that involve multi-locus effects across the HLA region.

Confirmation of novel type 1 diabetes risk loci in families.

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

No association of vitamin D intake or 25-hydroxyvitamin D levels in childhood with risk of islet autoimmunity and type 1 diabetes: the Diabetes Autoimmunity Study in the Young (DAISY).

AIMS/HYPOTHESIS: The aim of the study was to investigate the association between vitamin D intake and status and the risk of islet autoimmunity (IA) and subsequent type 1 diabetes in children at increased risk of type 1 diabetes. METHODS: The Diabetes Autoimmunity Study in the Young (DAISY) in Denver, CO, USA, has been following children at increased risk of diabetes since 1993. As of February 2011, 198 children developed IA during follow-up of 2,644 DAISY children. Vitamin D intake and plasma 25-hydroxyvitamin D [25(OH)D] were measured longitudinally. Proportional hazards regression analyses of time to IA, or type 1 diabetes in IA-positive children, were conducted, with vitamin D intake and 25(OH)D as time-varying covariates. HRs were calculated for a standard deviation difference in exposure, with adjustment for confounders. RESULTS: Intake of vitamin D was not associated with the risk of IA (adjusted HR 1.13; 95% CI 0.95, 1.35; p = 0.18) nor progression to diabetes in IA-positive children (adjusted HR 1.30; 95% CI 0.91, 1.86; p = 0.15). Moreover, 25(OH)D level was not associated with the risk of IA (adjusted HR 1.12; 95% CI 0.88, 1.43; p = 0.36), nor progression to diabetes in IA-positive children (adjusted HR 0.91; 95% CI 0.68, 1.22; p = 0.54). In the 128 children in whom we measured 25(OH)D at 9 months of age, 25(OH)D was not associated with risk of IA (n = 30 IA-positive children) (adjusted HR 1.02; 95% CI 0.96, 1.07; p = 0.58). CONCLUSIONS/INTERPRETATION: Neither vitamin D intake nor 25(OH)D levels throughout childhood were associated with the risk of IA or progression to type 1 diabetes in our population.

Pathway-based analysis of genetic susceptibility to cervical cancer in situ: HLA-DPB1 affects risk in Swedish women.

We have conducted a pathway-based analysis of genome-wide single-nucleotide polymorphism (SNP) data in order to identify genetic susceptibility factors for cervical cancer in situ. Genotypes derived from Affymetrix 500k or 5.0 arrays for 1076 cases and 1426 controls were analyzed for association, and pathways with enriched signals were identified using the SNP ratio test. The most strongly associated KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were Asthma (empirical P=0.03), Folate biosynthesis (empirical P=0.04) and Graft-versus-host disease (empirical P=0.05). Among the 11 top-ranking pathways were 6 related to the immune response with the common denominator being genes in the major histocompatibility complex (MHC) region on chromosome 6. Further investigation of the MHC revealed a clear effect of HLA-DPB1 polymorphism on disease susceptibility. At a functional level, DPB1 alleles associated with risk and protection differ in key amino-acid residues affecting peptide-binding motifs in the extracellular domains. The results illustrate the value of pathway-based analysis to mine genome-wide data, and point to the importance of the MHC region and specifically the HLA-DPB1 locus for susceptibility to cervical cancer.

A multi-site study using high-resolution HLA genotyping by next generation sequencing.

The high degree of polymorphism at human leukocyte antigen (HLA) class I and class II loci makes high-resolution HLA typing challenging. Current typing methods, including Sanger sequencing, yield ambiguous typing results because of incomplete genomic coverage and inability to set phase for HLA allele determination. The 454 Life Sciences Genome Sequencer (GS FLX) next generation sequencing system coupled with conexio atf software can provide very high-resolution HLA genotyping. High-throughput genotyping can be achieved by use of primers with multiplex identifier (MID) tags to allow pooling of the amplicons generated from different individuals prior to sequencing. We have conducted a double-blind study in which eight laboratory sites performed amplicon sequencing using GS FLX standard chemistry and genotyped the same 20 samples for HLA-A, -B, -C, DPB1, DQA1, DQB1, DRB1, DRB3, DRB4, and DRB5 (DRB3/4/5) in a single sequencing run. The average sequence read length was 250 base pairs and the average number of sequence reads per amplicon was 672, providing confidence in the allele assignments. Of the 1280 genotypes considered, assignment was possible in 95% of the cases. Failure to assign genotypes was the result of researcher procedural error or the presence of a novel allele rather than a failure of sequencing technology. Concordance with known genotypes, in cases where assignment was possible, ranged from 95.3% to 99.4% for the eight sites, with overall concordance of 97.2%. We conclude that clonal pyrosequencing using the GS FLX platform and CONEXIO ATF software allows reliable identification of HLA genotypes at high resolution.

An update to HLA nomenclature, 2010.

The WHO Nomenclature Committee for Factors of the HLA System met during the 15th International Histocompatibility and Immunogenetics Workshop in Buzios, Brazil in September 2008. This update is an extract of the main report that documents the additions and revisions to the nomenclature of human leukocyte antigen (HLA) specificities following the principles established in previous reports.

Overview of the Type I Diabetes Genetics Consortium.

The Type I Diabetes Genetics Consortium (T1DGC) is an international, multicenter research program with two primary goals. The first goal is to identify genomic regions and candidate genes whose variants modify an individual's risk of type I diabetes (T1D) and help explain the clustering of the disease in families. The second goal is to make research data available to the research community and to establish resources that can be used by, and that are fully accessible to, the research community. To facilitate the access to these resources, the T1DGC has developed a Consortium Agreement (http://www.t1dgc.org) that specifies the rights and responsibilities of investigators who participate in Consortium activities. The T1DGC has assembled a resource of affected sib-pair families, parent-child trios, and case-control collections with banks of DNA, serum, plasma, and EBV-transformed cell lines. In addition, both candidate gene and genome-wide (linkage and association) studies have been performed and displayed in T1DBase (http://www.t1dbase.org) for all researchers to use in their own investigations. In this supplement, a subset of the T1DGC collection has been used to investigate earlier published candidate genes for T1D, to confirm the results from a genome-wide association scan for T1D, and to determine associations with candidate genes for other autoimmune diseases or with type II diabetes that may be involved with beta-cell function.

Current status and the future for the genetics of type I diabetes.

The Type I Diabetes Genetics Consortium (T1DGC) is an international collaboration whose primary goal is to identify genes whose variants modify an individual's risk of type I diabetes (T1D). An integral part of the T1DGC's mission is the establishment of clinical and data resources that can be used by, and that are fully accessible to, the T1D research community (http://www.t1dgc.org). The T1DGC has organized the collection and analyses of study samples and conducted several major research projects focused on T1D gene discovery: a genome-wide linkage scan, an intensive evaluation of the human major histocompatibility complex, a detailed examination of published candidate genes, and a genome-wide association scan. These studies have provided important information to the scientific community regarding the function of specific genes or chromosomal regions on T1D risk. The results are continually being updated and displayed (http://www.t1dbase.org). The T1DGC welcomes all investigators interested in using these data for scientific endeavors on T1D. The T1DGC resources provide a framework for future research projects, including examination of structural variation, re-sequencing of candidate regions in a search for T1D-associated genes and causal variants, correlation of T1D risk genotypes with biomarkers obtained from T1DGC serum and plasma samples, and in-depth bioinformatics analyses.

Association analysis of SNPs in the IL4R locus with type I diabetes.

The Type I Diabetes Genetics Consortium (T1DGC) has collected thousands of multiplex and simplex families with type I diabetes (T1D) with the goal of identifying genes involved in T1D susceptibility. These families have all been genotyped for the HLA class I and class II loci and a subset of samples has been typed for an major histocompatibility complex (MHC) single-nucleotide polymorphism (SNP) panel. In addition, the T1DGC has genotyped SNPs in candidate genes to evaluate earlier reported T1D associations. Individual SNPs and SNP haplotypes in IL4R, which encodes the alpha-chain of the IL4 and IL13 receptors, have been associated with T1D in some reports, but not in others. In this study, 38 SNPs in IL4R were genotyped using the Sequenom iPLEX Gold MassARRAY technology in 2042 multiplex families from nine cohorts. Association analyses (transmission-disequilibrium test and parental-disequilibrium test) were performed on individual SNPs and on three-SNP haplotypes. Analyses were also stratified on the high-risk HLA DR3/DR4-DQB1*0302 genotype. A modest T1D association in HBDI families (n=282) was confirmed in this larger collection of HBDI families (n=424). The variant alleles at the non-synonymous SNPs (rs1805011 (E400A), rs1805012 (C431R), and rs1801275 (Q576R)), which are in strong linkage disequilibrium, were negatively associated with T1D risk. These SNPs were more associated with T1D among non-DR3/DR4-DQB1*0302 genotypes than DR3/DR4-DQB1*0302 genotypes. This association was stronger, both in terms of odds ratio and P-values, than the initial report of the smaller collection of HBDI families. However, the IL4R SNPs and the three-SNP haplotype containing the variant alleles were not associated with T1D in the total data. Thus, in the overall families, these results do not show evidence for an association of SNPs in IL4R with T1D.

Evidence for association of the TCF7 locus with type I diabetes.

The Type I Diabetes Genetics Consortium (T1DGC) has collected thousands of multiplex and simplex families with type I diabetes (T1D) with the goal of identifying genes involved in T1D susceptibility. These families have been genotyped for the HLA class I and class II loci and, recently, for a genome-wide panel of single-nucleotide polymorphisms (SNPs). In addition, multiple SNPs in specific candidate genes have been genotyped in these families in an attempt to evaluate previously reported T1D associations, including the C883A (Pro-Thr) polymorphism in exon 2 of TCF7, a T-cell transcription factor. The TCF7 883A allele was associated with T1D in subjects with T1D not carrying the high-risk HLA genotype DR3/DR4. A panel of 11 SNPs in TCF7 was genotyped in 2092 families from 9 cohorts of the T1DGC. SNPs at two positions in TCF7 were associated with T1D. One associated SNP, C883A (rs5742913), was reported earlier to have a T1D association. A second SNP, rs17653687, represents a novel T1D susceptibility allele in TCF7. After stratification on the high T1D risk DR3/DR4 genotype, the variant (A) allele of C883A was significantly associated with T1D among non-DR3/DR4 cases (transmission=55.8%, P=0.004; OR=1.26) but was not significantly associated in the DR3/DR4 patient subgroup, replicating the earlier report. The reference A allele of intronic SNP rs17653687 was modestly associated with T1D in both DR3/DR4 strata (transmission=54.4% in DR3/DR4; P=0.03; transmission=52.9% in non-DR3/DR4; P=0.03). These results support the previously reported association of the non-synonymous Pro-Thr SNP in TCF7 with T1D, and suggest that other alleles at this locus may also confer risk.