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.

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.

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.

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.

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.

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.

High-resolution, high-throughput HLA genotyping by next-generation sequencing.

The human leukocyte antigen (HLA) class I and class II loci are the most polymorphic genes in the human genome. Hematopoietic stem cell transplantation requires allele-level HLA typing at multiple loci to select the best matched unrelated donors for recipient patients. In current methods for HLA typing, both alleles of a heterozygote are amplified and typed or sequenced simultaneously, often making it difficult to unambiguously determine the sequence of the two alleles. Next-generation sequencing methods clonally propagate in parallel millions of single DNA molecules, which are then also sequenced in parallel. Recently, the read lengths obtainable by one such next-generation sequencing method (454 Life Sciences, Inc.) have increased to >250 nucleotides. These clonal read lengths make possible setting the phase of the linked polymorphisms within an exon and thus the unambiguous determination of the sequence of each HLA allele. Here we demonstrate this capacity as well as show that the throughput of the system is sufficiently high to enable a complete, 7-locus HLA class I and II typing for 24 or 48 individual DNAs in a single GS FLX sequencing run. Highly multiplexed amplicon sequencing is facilitated by the use of sample-specific internal sequence tags (multiplex identification tags or MIDs) in the primers that allow pooling of samples yet maintain the ability to assign sequences to specific individuals. We have incorporated an HLA typing software application developed by Conexio Genomics (Freemantle, Australia) that assigns HLA genotypes for these 7 loci (HLA-A, -B, -C, DRB1, DQA1, DQB1, DPB1), as well as for DRB3, DRB4, and DRB5 from 454 sequence data. The potential of this HLA sequencing system to analyze chimeric mixtures is demonstrated here by the detection of a rare HLA-B allele in a mixture of two homozygous cell lines (1/100), as well as by the detection of the rare nontransmitted maternal allele present in the blood of a severe combined immunodeficiency disease syndrome (SCIDS) patient.

The identification of a new HLA-DPB1 allele (*1302) in one family and the detection of a recombination event between the DR and the DQ regions in another Caucasian T1DGC family.

The analysis of families collected by the T1DGC and typed at high resolution for the HLA class I and class II loci provided an opportunity for identifying new alleles and rare recombination events. In one American Caucasian family, a novel allele (HLA-DPB1*1302), detected as an unusual pattern of probe binding, was identified in the mother and in one child. Amplicons from both individuals were sequenced and a new variant of DPB1*1301 with an A->T mutation [TAC to TTC in codon 64, (amino acid 35); Y to F] was confirmed. In another American Caucasian family, one child inherited an unusual haplotype, DRB1*1501-DQA1*0102-DQB1*0609-DPA1*0103-DPB1*0601 resulting from a recombination between the DRB1 loci on the maternal chromosomes DRB1*1501-DQA1*0102-DQB1*0602-DPA1*0103-DPB1*0401 and DRB1*1302-DQA1*0102-DQB1*0609-DPA1*0103-DPB1*0601.

The frequent and conserved DR3-B8-A1 extended haplotype confers less diabetes risk than other DR3 haplotypes.

AIM: The goal of this study was to develop and implement methodology that would aid in the analysis of extended high-density single nucleotide polymorphism (SNP) major histocompatibility complex (MHC) haplotypes combined with human leucocyte antigen (HLA) alleles in relation to type 1 diabetes risk. METHODS: High-density SNP genotype data (2918 SNPs) across the MHC from the Type 1 Diabetes Genetics Consortium (1240 families), in addition to HLA data, were processed into haplotypes using PedCheck and Merlin, and extended DR3 haplotypes were analysed. RESULTS: With this large dense set of SNPs, the conservation of DR3-B8-A1 (8.1) haplotypes spanned the MHC (>/=99% SNP identity). Forty-seven individuals homozygous for the 8.1 haplotype also shared the same homozygous genotype at four 'sentinel' SNPs (rs2157678 'T', rs3130380 'A', rs3094628 'C' and rs3130352 'T'). Conservation extended from HLA-DQB1 to the telomeric end of the SNP panels (3.4 Mb total). In addition, we found that the 8.1 haplotype is associated with lower risk than other DR3 haplotypes by both haplotypic and genotypic analyses [haplotype: p = 0.009, odds ratio (OR) = 0.65; genotype: p = 6.3 x 10(-5), OR = 0.27]. The 8.1 haplotype (from genotypic analyses) is associated with lower risk than the high-risk DR3-B18-A30 haplotype (p = 0.01, OR = 0.23), but the DR3-B18-A30 haplotype did not differ from other non-8.1 DR3 haplotypes relative to diabetes association. CONCLUSION: The 8.1 haplotype demonstrates extreme conservation (>3.4 Mb) and is associated with significantly lower risk for type 1 diabetes than other DR3 haplotypes.

MHC loci affecting cervical cancer risk: distinguishing the effects of HLA-DQB1 and non-HLA genes TNF, LTA, TAP1 and TAP2.

Cervical cancer has been associated with specific human leukocyte antigen (HLA) haplotypes/alleles and with polymorphisms at the nearby non-HLA loci TNF, LTA, TAP1 and TAP2. Distinguishing effects of individual loci in the major histocompatibility complex (MHC) region are difficult due to the complex linkage disequilibrium (LD) pattern characterized by high LD, punctuated by recombination hot spots. We have evaluated the association of polymorphism at HLA class II DQB1 and the TNF, LTA, TAP1 and TAP2 genes with cervical cancer risk, using 1306 familial cases and 288 controls. DQB1 was strongly associated; alleles *0301, *0402 and (*)0602 increased cancer susceptibility, whereas *0501 and *0603 decreased susceptibility. Among the non-HLA loci, association was only detected for the TAP2 665 polymorphism, and interallelic disequilibrium analysis indicated that this could be due to LD with DQB1. As the TAP2 665 association was seen predominantly in non-carriers of DQB1 susceptibility alleles, we hypothesized that TAP2 665 may have an effect not attributable to LD with DQB1. However, a logistic regression analysis suggested that TAP2 665 was strongly influenced by LD with DQB1. Our results emphasize the importance of large sample sizes and underscore the necessity of examining both HLA and non-HLA loci in the MHC to assign association to the correct locus.

Differential effects of DRB1*0301 and DQA1*0501-DQB1*0201 on the activation and progression of islet cell autoimmunity.

Autoimmune diabetes shows extreme variation in age of onset and clinical presentation, although most studies have been done in children with the most severe subtype. Disease risk is strongly associated with HLA-DRB1*0301-DQA1*0501-DQB1*0201 (DR3-DQ2), but it has not been possible to separate the effects of the DR and DQ alleles. We have identified a large Bedouin kindred in which a high prevalence of islet autoimmunity is associated with two different DR3 haplotypes, one carrying the usual DQ2 and the other carrying DQA1*0102-DQB1*0502 (DQ5). Results of prospective follow-up studies indicate that DR3 is associated with the initial activation of islet autoimmunity whereas DQ2 is associated with early-onset and severe clinical disease. The association signals map to a 350-kb interval, thus implicating primary effects for DR3 and DQ2. Overall, our results emphasize the importance of prospective genetic studies that examine the full range of variation in the initiation, progression and expression of autoimmune disease.

Relative predispositional effects of HLA class II DRB1-DQB1 haplotypes and genotypes on type 1 diabetes: a meta-analysis.

The direct involvement of the human leukocyte antigen class II DR-DQ genes in type 1 diabetes (T1D) is well established, and these genes display a complex hierarchy of risk effects at the genotype and haplotype levels. We investigated, using data from 38 studies, whether the DR-DQ haplotypes and genotypes show the same relative predispositional effects across populations and ethnic groups. Significant differences in risk within a population were considered, as well as comparisons across populations using the patient/control (P/C) ratio. Within a population, the ratio of the P/C ratios for two different genotypes or haplotypes is a function only of the absolute penetrance values, allowing ranking of risk effects. Categories of consistent predisposing, intermediate ('neutral'), and protective haplotypes were identified and found to correlate with disease prevalence and the marked ethnic differences in DRB1-DQB1 frequencies. Specific effects were identified, for example for predisposing haplotypes, there was a statistically significant and consistent hierarchy for DR4 DQB1*0302s: DRB1*0405 =*0401 =*0402 > *0404 > *0403, with DRB1*0301 DQB1*0200 (DR3) being significantly less predisposing than DRB1*0402 and more than DRB1*0404. The predisposing DRB1*0401 DQB1*0302 haplotype was relatively increased compared with the protective haplotype DRB1*0401 DQB1*0301 in heterozygotes with DR3 compared with heterozygotes with DRB1*0101 DQB1*0501 (DR1). Our results show that meta-analyses and use of the P/C ratio and rankings thereof can be valuable in determining T1D risk factors at the haplotype and amino acid residue levels.

14th International HLA and Immunogenetics Workshop: report of progress in methodology, data collection, and analyses.

The Biostatistics Component of the 13th International Histocompatibility Workshop (IHWS) developed the PyPop (Python for Population Genomics) software framework for high-throughput analysis and quality control (QC) assessments of highly polymorphic genotype data. Since its initial release, the software has had several new analysis modules added to it. These additions, combined with improved data filtering and QC modules, facilitate analyses of data at different levels (allele, haplotype, amino acid sequence, and nucleotide sequence). Since the 13th IHWS, much of the human leukocyte antigen (HLA) data from the workshop, QCed via PyPop and other methods, have been made publicly available through the Major Histocompatibility Complex database web site at the National Center for Biotechnology Information (http://ncbi.nih.gov/mhc/). The Anthropology/Human Genetic Diversity component (AHGDC) data have been used in a variety of studies. Prugnolle et al. used this data to corroborate a model of pathogen-driven selection as a factor related to high levels of diversity at HLA loci. Using a comparative genomics approach contrasting results for HLA and non-HLA markers, Meyer et al. analyzed a subset of the 13th IHWS AHGDC data and showed that HLA loci show detectable signs of both natural selection and the demographic history of populations.