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.

Recent origin of HLA-DRB1 alleles and implications for human evolution.

The HLA class I and class II loci are the most highly polymorphic coding regions in the human genome. Based on the similarity of the coding sequences of alleles between species, it has been claimed that the HLA polymorphism is ancient and predates the separation of human (Homo) and chimpanzee (Pan), 4-7.4 Myr ago. Analysis of intron sequences, however, provides support for a more recent origin and for rapid generation of alleles at the HLA class II DRB1 locus. The human DRB1 alleles can be divided into groups (allelic lineages); most of these lineages have diverged from each other before the separation of Homo and Pan. Alleles within such a lineage, however, appear to be, on average, 250,000 years old, implying that the vast majority (greater than 90%) of the more than 135 contemporary human DRB1 alleles have been generated after the separation of Homo and Pan. The coalescence time of alleles within allelic lineages indicates that the effective population size (Ne) for early hominids (over the last 1 Myr) was approximately 10(4) individuals, similar to estimates based on other nuclear loci and mitochondrial DNA. With a single exception, the genetic mechanisms (gene conversion and point mutation) that have diversified the exon-2 sequences do not appear to extend into the adjacent intron sequences. The part of exon 2 encoding the beta-sheet evolves in concert with the surrounding introns, while the alpha-helix appears to have been subjected to gene conversion-like events, suggesting that such exchange events are highly localised and occur over extremely short sequence tracts.

HLA DR-DQ associations with cervical carcinoma show papillomavirus-type specificity.

Cervical carcinoma is now known to be associated with human papillomaviruses (HPV), but the evidence for a link with specific HLA loci is controversial. The role of genetic variation at the HLA class II loci and among HPV types in cervical carcinoma was investigated by PCR DNA amplification and oligonucleotide probe typing of paraffin-embedded invasive cervical cancer tissue from Hispanic patients and of cervical swabs from Hispanic controls. Certain HLA class II haplotypes (such as DRB1*1501-DQB1*0602) were associated significantly, while DR13 haplotypes were negatively associated with cervical carcinoma. These associations are HPV16-type specific. These results suggest that specific HLA class II haplotypes may influence the immune response to specific HPV-encoded epitopes and affect the risk of cervical neoplasia.

HLA class II alleles and susceptibility and resistance to insulin dependent diabetes mellitus in Mexican-American families.

The role of HLA class II alleles in genetic predisposition to insulin dependent diabetes mellitus (IDDM) was examined by PCR/oligonucleotide probe typing of 42 Mexican-American IDDM families derived from Hispanic Caucasians and Native Americans. All high risk haplotypes (HLA-DR3 and DR4) were of European origin while the most strongly protective haplotype (DRB1*1402) was Native American. Of the 16 DR-DQ DR4 haplotypes identified, only those bearing DQB1*0302 conferred risk; the DRB1 allele, however, also markedly influenced IDDM risk. The general pattern of neutral and protective haplotypes indicates that the presence of Asp-57 in the HLA-DQ beta chain does not confer IDDM protection per se and indicates that both DRB1 and DQB1 influence IDDM susceptibility as well as protection.

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.

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.

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.

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.

Molecular genetic analysis of 178 I-Abm12-reactive T cells.

We have studied the genetic diversity of the TCR repertoire to the murine alloantigen I-Abm12 by generating a panel of 178 C57BL/10-derived I-Abm12-reactive T cell hybridomas. The expression of V alpha and V beta gene families was examined in this panel and the frequency of expression of V beta, but not ofV alpha, gene families differed significantly from that observed in a companion panel of random C57BL/10-derived hybridomas. The V beta 5 gene family was expressed significantly less frequently while the V beta 14, V beta 15, and V beta 16 genes were expressed significantly more frequently in the panel of I-Abm12-reactive than in the panel of random hybridomas. The junctional regions (VJ alpha and VDJ beta) of TCR V alpha and V beta genes from selected I-Abm12-specific hybridomas were amplified using the polymerase chain reaction, and directly sequenced. Surprisingly, no conserved J alpha, D beta, J beta, or N region-encoded sequences among these selected I-Abm12-reactive TCRs were identified. Thus, the T cell response to an I-A alloantigen that differs by only three amino acid residues from the I-A molecule of the responding strain is genetically complex but nonrandom. We have estimated that the repertoire to this alloantigen is comprised of at least 37 different TCRs.

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.

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.

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.

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.

Recent advances in the polymerase chain reaction.

The polymerase chain reaction (PCR) has dramatically altered how molecular studies are conducted as well as what questions can be asked. In addition to simplifying molecular tasks typically carried out with the use of recombinant DNA technology, PCR has allowed a spectrum of advances ranging from the identification of novel genes and pathogens to the quantitation of characterized nucleotide sequences. PCR can provide insights into the intricacies of single cells as well as the evolution of species. Some recent developments in instrumentation, methodology, and applications of the PCR are presented in this review.

Direct cloning and sequence analysis of enzymatically amplified genomic sequences.

A method is described for directly cloning enzymatically amplified segments of genomic DNA into an M13 vector for sequence analysis. A 110-base pair fragment of the human beta-globin gene and a 242-base pair fragment of the human leukocyte antigen DQ alpha locus were amplified by the polymerase chain reaction method, a procedure based on repeated cycles of denaturation, primer annealing, and extension by DNA polymerase I. Oligonucleotide primers with restriction endonuclease sites added to their 5' ends were used to facilitate the cloning of the amplified DNA. The analysis of cloned products allowed the quantitative evaluation of the amplification method's specificity and fidelity. Given the low frequency of sequence errors observed, this approach promises to be a rapid method for obtaining reliable genomic sequences from nanogram amounts of DNA.

Segregation and mapping analysis of polymorphic HLA class I restriction fragments: detection of a novel fragment.

An HLA-B7 complementary DNA clone was used as a hybridization probe to analyze the segregation pattern of polymorphic class I restriction fragments in several families whose HLA types had been determined by serological techniques. In one family in which a crossover in the HLA region had occurred, a specific genomic fragment was mapped with respect to the crossover. In another family, a novel genomic fragment present in one child and absent in all other family members was observed. With the exception of this novel fragment, all polymorphic class I fragments observed in this study segregated with a serologically defined parental haplotype, a result consistent with HLA linkage.

Enzymatic amplification of beta-globin genomic sequences and restriction site analysis for diagnosis of sickle cell anemia.

Two new methods were used to establish a rapid and highly sensitive prenatal diagnostic test for sickle cell anemia. The first involves the primer-mediated enzymatic amplification of specific beta-globin target sequences in genomic DNA, resulting in the exponential increase (220,000 times) of target DNA copies. In the second technique, the presence of the beta A and beta S alleles is determined by restriction endonuclease digestion of an end-labeled oligonucleotide probe hybridized in solution to the amplified beta-globin sequences. The beta-globin genotype can be determined in less than 1 day on samples containing significantly less than 1 microgram of genomic DNA.