Rapid high-throughput human leukocyte antigen typing by massively parallel pyrosequencing for high-resolution allele identification

Advancements in HLA Typing Techniques and Their Impact on Transplantation Medicine

HLA typing serves as a standard practice in hematopoietic stem cell transplantation to ensure compatibility between donors and recipients, preventing the occurrence of allograft rejection and graft-versus-host disease. Conventional laboratory methods that have been widely employed in the past few years, including sequence-specific primer PCR and sequencing-based typing (SBT), currently face the risk of becoming obsolete. This risk stems not only from the extensive diversity within HLA genes but also from the rapid advancement of next-generation sequencing and third-generation sequencing technologies. Third-generation sequencing systems like single-molecule real-time (SMRT) sequencing and Oxford Nanopore (ONT) sequencing have the capability to analyze long-read sequences that span entire intronic-exonic regions of HLA genes, effectively addressing challenges related to HLA ambiguity and the phasing of multiple short-read fragments. The growing dominance of these advanced sequencers in HLA typing is expected to solidify further through ongoing refinements, cost reduction, and error rate minimization. This review focuses on hematopoietic stem cell transplantation (HSCT) and explores prospective advancements and application of HLA DNA typing techniques. It explores how the adoption of third-generation sequencing technologies can revolutionize the field by offering improved accuracy, reduced ambiguity, and enhanced assessment of compatibility in HSCT. Embracing these cutting-edge technologies is essential to advancing the success rates and outcomes of hematopoietic stem cell transplantation. This review underscores the importance of staying at the forefront of HLA typing techniques to ensure the best possible outcomes for patients undergoing HSCT.

Comprehensive HLA Typing from a Current Allele Database Using Next-Generation Sequencing Data

HLA allele information is essential for a variety of medical applications, such as genomic studies of multifactorial diseases, including immune system and inflammation-related disorders, and donor selection in organ transplantation and regenerative medicine. To obtain this information, an accurate HLA typing method that is applicable for any allele registered in HLA allele databases is needed. Here we describe a method-called HLA-HD-for determining alleles from a current HLA database using next-generation sequencing (NGS) results.

A novel swab storage gel is superior to dry swab DNA collection, and enables long-range high resolution NGS HLA typing from buccal cell samples

HLA sequence-based DNA typing (SBT) by long-range PCR amplification (LR PCR) and next-generation sequencing (NGS) is a high-throughput DNA sequencing method (LR-NGS-SBT) for the efficient and sensitive detection of novel and null HLA alleles to the field-4 level of allelic resolution without phase ambiguity. However, the accuracy and reliability of the HLA typing results using buccal cells (BCs) and saliva as genetic source materials for the LR-NGS-SBT method are dependent largely on the quality of the extracted genomic DNA (gDNA) because a large degree of gDNA fragmentation can result in insufficient PCR amplification with the incorrect assignment of HLA alleles due to allele dropouts. In this study, we developed a new cost-efficient swab storage gel (SSG) for wet swab collection of BCs (BC-SSG) and evaluated its usefulness by performing different DNA analytical parameters including LR-NGS-SBT to compare the quality and quantity of gDNA extracted from BCs (in SSG or air dried), blood and saliva of 30 subjects. The BC-SSG samples after 5 days of storage revealed qualitative and quantitative DNA values equivalent to that of blood and/or saliva and better than swabs that were only air-dried (BC-nSSG). Moreover, all the gDNA extracted from blood, saliva and BC-SSG samples were HLA-typed successfully to an equivalent total of 408 alleles for each sample type. Therefore, the application of BC-SSG collection media for LR-NGS-SBT has benefits over BC dried samples (dry swabs) such as reducing retesting and the number of untestable BC samples due to insufficient DNA amplification. This article is protected by copyright. All rights reserved.

High-resolution HLA typing by long reads from the R10.3 Oxford nanopore flow cells

Nanopore sequencing has been investigated as a rapid and cost-efficient option for HLA typing in recent years. Despite the lower raw read accuracy, encouraging typing accuracy has been reported, and long reads from the platform offer additional benefits of the improved phasing of distant variants. The newly released R10.3 flow cells are expected to provide higher read-level accuracy than previous chemistries. We examined the performance of R10.3 flow cells on the MinION device in HLA typing after enrichment of target genes by multiplexed PCR. We also aimed to mimic a 1-day workflow with 8-24 samples per sequencing run. A diverse collection of 102 unique samples were typed for HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1, -DRB3/4/5 loci. The concordance rates at 2-field and 3-field resolutions were 99.5% (1836 alleles) and 99.3% (1710 alleles). We also report important quality metrics from these sequencing runs. Continued research and independent validations are warranted to increase the robustness of nanopore-based HLA typing for broad clinical application.

Distributions of HLA-A, -B, and -DRB1 alleles typed by amplicon-based next generation sequencing in Korean volunteer donors for unrelated hematopoietic stem cell transplantation

HLA genes play a pivotal role for successful hematopoietic stem cell transplantation (HSCT). There is an increasing need for sophisticated screening of donor HLA genotypes for unrelated HSCT. Next generation sequencing (NGS) has emerged as an alternative for classical Sanger sequence for HLA typing. In this study, HLA-A, -B, and -DRB1 alleles were genotyped at the allelic (6-digit) level using MiSeqDx in 26,202 volunteers from the Korean Network for Organ Sharing. Exon 2 and 3 of HLA-A and -B and exon 2 of HLA-DRB1 were amplified by polymerase chain reaction (PCR) and each allele was determined by matching the targeted exons and the reference sequence consisting of the IPD-IMGT/HLA Database. Seventy alleles of HLA-A, 102 alleles of HLA-B, and 69 alleles of HLA-DRB1 were identified. According to common and well-documented catalogs, 34 alleles in HLA-A, 61 in HLA-B, and 45 in HLA-DRB1 locus were common alleles, and 12, 14, and 11 kinds, were well-documented alleles, respectively. Thirteen novel alleles including 3 alleles in HLA-A, 8 alleles in HLA-B, and 2 alleles in HLA-DRB1 loci were found. Ten haplotypes with a frequency of more than 1.0% accounted for 22.4% of the total haplotype frequencies. Cis/trans ambiguities of HLA-A and -B loci by combination of exons 2 and 3 were analyzed to be 0.17% of 3 and 3.95% of 22 genotypes, respectively. This information on rare and novel alleles found by accurate HLA typing with NGS may be helpful for unrelated HSCT among Koreans.

Utilizing nanopore sequencing technology for the rapid and comprehensive characterization of eleven HLA loci; addressing the need for deceased donor expedited HLA typing

The comprehensive characterization of human leukocyte antigen (HLA) genomic sequences remains a challenging problem. Despite the significant advantages of next-generation sequencing (NGS) in the field of Immunogenetics, there has yet to be a single solution for unambiguous, accurate, simple, cost-effective, and timely genotyping necessary for all clinical applications. This report demonstrates the benefits of nanopore sequencing introduced by Oxford Nanopore Technologies (ONT) for HLA genotyping. Samples (n = 120) previously characterized at high-resolution three-field (HR-3F) for 11 loci were assessed using ONT sequencing paired to a single-plex PCR protocol (Holotype) and to two multiplex protocols OmniType (Omixon) and NGSgo^®-MX6-1 (GenDx). The results demonstrate the potential of nanopore sequencing for delivering accurate HR-3F typing with a simple, rapid, and cost-effective protocol. The protocol is applicable to time-sensitive applications, such as deceased donor typings, enabling better assessments of compatibility and epitope analysis. The technology also allows significantly shorter turnaround time for multiple samples at a lower cost. Overall, the nanopore technology appears to offer a significant advancement over current next-generation sequencing platforms as a single solution for all HLA genotyping needs.

A long road/read to rapid high-resolution HLA typing: The nanopore perspective

Next-generation sequencing (NGS) has been widely adopted for clinical HLA typing and advanced immunogenetics researches. Current methodologies still face challenges in resolving cis-trans ambiguity involving distant variant positions, and the turnaround time is affected by testing volume and batching. Nanopore sequencing may become a promising addition to the existing options for HLA typing. The technology delivered by the MinION sequencer of Oxford Nanopore Technologies (ONT) can record the ionic current changes during the translocation of DNA/RNA strands through transmembrane pores and translate the signals to sequence reads. It features simple and flexible library preparations, long sequencing reads, portable and affordable sequencing devices, and rapid, real-time sequencing. However, the error rate of the sequencing reads is high and remains a hurdle for its broad application. This review article will provide a brief overview of this technology and then focus on the opportunities and challenges of using nanopore sequencing for high-resolution HLA typing and immunogenetics research.

Bioinformatics Strategies, Challenges, and Opportunities for Next Generation Sequencing-Based HLA Genotyping

The advent of next generation sequencing (NGS) has altered the face of genotyping the human leukocyte antigen (HLA) system in clinical, stem cell donor registry, and research contexts. NGS has led to a dramatically increased sequencing throughput at high accuracy, while being more time and cost efficient than precursor technologies. This has led to a broader and deeper profiling of the key genes in the human immunogenetic make-up. The rapid evolution of sequencing technologies is evidenced by the development of varied short-read sequencing platforms with differing read lengths and sequencing capacities to long-read sequencing platforms capable of profiling full genes without fragmentation. Concomitantly, there has been development of a diverse set of computational analyses and software tools developed to deal with the various strengths and limitations of the sequencing data generated by the different sequencing platforms. This review surveys the different modalities involved in generating NGS HLA profiling sequence data. It systematically describes various computational approaches that have been developed to achieve HLA genotyping to different degrees of resolution. At each stage, this review enumerates the drawbacks and advantages of each of the platforms and analysis approaches, thus providing a comprehensive picture of the current state of HLA genotyping technologies.

Single molecule real-time (SMRT) sequencing comes of age: applications and utilities for medical diagnostics

Short read massive parallel sequencing has emerged as a standard diagnostic tool in the medical setting. However, short read technologies have inherent limitations such as GC bias, difficulties mapping to repetitive elements, trouble discriminating paralogous sequences, and difficulties in phasing alleles. Long read single molecule sequencers resolve these obstacles. Moreover, they offer higher consensus accuracies and can detect epigenetic modifications from native DNA. The first commercially available long read single molecule platform was the RS system based on PacBio's single molecule real-time (SMRT) sequencing technology, which has since evolved into their RSII and Sequel systems. Here we capsulize how SMRT sequencing is revolutionizing constitutional, reproductive, cancer, microbial and viral genetic testing.

Application of high-throughput next-generation sequencing for HLA typing of DNA extracted from postprocessing cord blood units

Cord blood has become an acceptable source of hematopoietic stem cells for transplantation. HLA plays a major role in hematopoietic stem cell transplantation (HSCT). Typing of cord blood samples for HLA alleles has been performed based on the serological and molecular methods. However, with the advent of next-generation sequencing technology, HLA typing becomes more accurate and unambiguous (upto intron level). Contamination of cord blood cells with erythropoietic cells poses a challenge in DNA extraction and downstream application. In the present study, DNA extracted from buffy coat of cord blood samples was typed for HLA-A, -B, -C, DRB1, and DQB1 alleles by Illumina miniseq and the sequences were aligned, phased, and mapped by MIA FORA software algorithms. Most frequent alleles found were HLA A*01:01:01 (17%), A*24:02:01 (15.1%), A*11:01:01 (13.6%), B*40:06:01 (10.7%), C*06:02:01 (17.7%), C*04:01:01 (14.2%), C*15:02:01 (11.4%), C*07:02:01 (10.7%), DRB1*07:01:01 (15.9%), DRB1*10:01:01 (10.2%), DQB1*06:01:01 (17.4%), DQB1*05:01:01 (12.4%), and DQB1*05:03:01 (10.4%). One null allele (A*24:11N), two novel alleles in B loci and three rare alleles (B*40:06:04, B*51:01:05, and C*01:44) were also identified in the present study. This study shows that high-throughput, unambiguous (third-field resolution) HLA typing can be performed on cord blood samples. In order to preserve the precious sample for future use, minimal amount of cord blood samples (postprocessing) could be used for HLA typing purpose.

Super High Resolution for Single Molecule-Sequence-Based Typing of Classical HLA Loci Using Ion Torrent PGM

Super high resolution-single molecule-sequence-based typing (SS-SBT) is an HLA DNA typing method to the field 4 level of allelic resolution (formerly known as 8-digit typing) to efficiently detect novel and null alleles without phase ambiguity by combination of long ranged PCR amplification and next-generation sequencing (NGS) technologies. In this chapter, we describe three basic steps, long ranged PCR, NGS, and allele assignment.

"Worldwide Network for Blood & Marrow Transplantation (WBMT) special article, challenges facing emerging alternate donor registries"

Hematopoietic cell transplantation (HCT) activity is increasing at an unprecedented pace with > 50,000 allogeneic transplants occurring annually worldwide. Establishing a functional HCT donor registry can be very challenging with respect to ethnicities, financial, technical, and geopolitical issues. Extensive planning steps are essential to overcome the expected challenges while establishing the registry, and to maintain its functionality. A few strategies can help move past those challenges and push the development of such registries forward. Authorities involved in HCT donor registry establishment will have to balance the advantages and costs of such a project and accommodate the emerging alternatives such as cord blood or related haploidentical transplants. Miscalculations and incomplete understanding of the various aspects of the process can have tremendous impact on the optimization of a HCT donor registry especially in developing countries. Herein we present some challenges in establishing such a registry and present potential solutions.

The MHC in the era of next-generation sequencing: Implications for bridging structure with function

The MHC continues to have the most disease-associations compared to other regions of the human genome, even in the genome-wide association study (GWAS) and single nucleotide polymorphism (SNP) era. Analysis of non-coding variation and their impact on the level of expression of HLA allotypes has shed new light on the potential mechanisms underlying HLA disease associations and alloreactivity in transplantation. Next-generation sequencing (NGS) technology has the capability of delineating the phase of variants in the HLA antigen-recognition site (ARS) with non-coding regulatory polymorphisms. These relationships are critical for understanding the qualitative and quantitative implications of HLA gene diversity. This article summarizes current understanding of non-coding region variation of HLA loci, the consequences of regulatory variation on HLA expression, the role for evolution in shaping lineage-specific expression, and the impact of HLA expression on disease susceptibility and transplantation outcomes. A role for phased sequencing methods for the MHC, and perspectives for future directions in basic and applied immunogenetic studies of the MHC are presented.

Revisiting the potential power of human leukocyte antigen (HLA) genes on relationship testing by massively parallel sequencing-based HLA typing in an extended family

The human leukocyte antigen (HLA) genes are the most polymorphic genes in the human genome and have great power in forensic applications, especially in relationship testing and personal identification. However, the extreme polymorphism of HLA has made unambiguous genotyping of these genes very challenging and resulted in the limited application in relationship testing. Fortunately, massively parallel sequencing (MPS) technology offers the promise of unambiguous and high-throughput HLA typing. In this study, 11 HLA genes were typed in one extended family residing in North China and encompassing six generations. Phase-resolved genotypes for HLA genes were generated and HLA haplotype structure was defined. The paternity/kinship index, or in other words, likelihood ratio (LR) was calculated. A total of 88 alleles were identified, of which eight alleles were newly discovered. The inheritance of HLA alleles followed Mendelian law. With the discovery of new HLA alleles and three recombination events, a total of eleven new HLA haplotypes were identified in this population. LR distribution showed that, when HLA alleles were applied, the Log₁₀LR for a single locus could reach very high and the median average Log₁₀LRs of HLA genes were much higher than that of short tandem repeat loci. The result showed that high-throughput HLA genotyping could be achieved rapidly by MPS, and the contribution of HLA genes on system performance could be high, which may be applied as a supplement in forensic genetics studies. This study was also valuable in demonstrating the genetic mechanisms governing the generation of polymorphisms of the HLA genes.

Using Exome and Amplicon-Based Sequencing Data for High-Resolution HLA Typing with ATHLATES

ATHLATES (accurate typing of human leukocyte antigen through exome sequencing) was originally developed to analyze whole-exome sequencing (exome-seq) data from the Illumina platform and to predict the HLA genotype at 2-field or higher resolution. HLA locus-specific reads are first collected by stringent read mapping to the IMGT/HLA database. ATHLATES then performs read assembly, candidate allele identification, and genotype inference. Here, we describe the protocol of using ATHLATES for the above purpose and expand the application to analyze targeted sequencing data using amplicons of full HLA genes.

Comprehensive HLA Typing from a Current Allele Database Using Next-Generation Sequencing Data

HLA allele information is essential for a variety of medical applications, such as genomic studies of multifactorial diseases, including immune system and inflammation-related disorders, and donor selection in organ transplantation and regenerative medicine. To obtain this information, an accurate HLA typing method that is applicable for any allele registered in HLA allele databases is needed. Here, we describe a method for determining alleles from a current HLA database using next-generation sequencing (NGS) results.

HLA genotyping by next-generation sequencing of complementary DNA

Background

Genotyping of the human leucocyte antigen (HLA) is indispensable for various medical treatments. However, unambiguous genotyping is technically challenging due to high polymorphism of the corresponding genomic region. Next-generation sequencing is changing the landscape of genotyping. In addition to high throughput of data, its additional advantage is that DNA templates are derived from single molecules, which is a strong merit for the phasing problem. Although most currently developed technologies use genomic DNA, use of cDNA could enable genotyping with reduced costs in data production and analysis. We thus developed an HLA genotyping system based on next-generation sequencing of cDNA.

Methods

Each HLA gene was divided into 3 or 4 target regions subjected to PCR amplification and subsequent sequencing with Ion Torrent PGM. The sequence data were then subjected to an automated analysis. The principle of the analysis was to construct candidate sequences generated from all possible combinations of variable bases and arrange them in decreasing order of the number of reads. Upon collecting candidate sequences from all target regions, 2 haplotypes were usually assigned. Cases not assigned 2 haplotypes were forwarded to 4 additional processes: selection of candidate sequences applying more stringent criteria, removal of artificial haplotypes, selection of candidate sequences with a relaxed threshold for sequence matching, and countermeasure for incomplete sequences in the HLA database.

Results

The genotyping system was evaluated using 30 samples; the overall accuracy was 97.0% at the field 3 level and 98.3% at the G group level. With one sample, genotyping of DPB1 was not completed due to short read size. We then developed a method for complete sequencing of individual molecules of the DPB1 gene, using the molecular barcode technology.

Conclusion

The performance of the automatic genotyping system was comparable to that of systems developed in previous studies. Thus, next-generation sequencing of cDNA is a viable option for HLA genotyping.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4300-7) contains supplementary material, which is available to authorized users.

HLA-HD: An accurate HLA typing algorithm for next-generation sequencing data

The accurate typing of human leukocyte antigen (HLA) alleles is critical for a variety of medical applications, such as genomic studies of multifactorial diseases, including immune system and inflammation-related disorders, and donor selection in organ transplantation and regenerative medicine. Here, we developed a new algorithm for determining HLA alleles using next-generation sequencing (NGS) results. The method consists of constructing an extensive dictionary of HLA alleles, precise mapping of the NGS reads, and calculating a score based on weighted read counts to select the most suitable pair of alleles. The developed algorithm compares the score of all allele pairs, taking into account variation not only in the domain for antigen presentation (G-DOMAIN), but also outside this domain. Using this method, HLA alleles could be determined with 6-digit precision. We showed that our method was more accurate than other NGS-based methods and revealed limitations of the conventional HLA typing technologies. Furthermore, we determined the complete genomic sequence of an HLA-A-like-pseudogene when we assembled NGS reads that had caused arguable typing, and found its identity with HLA-Y*02:01. The accuracy of the HLA-A allele typing was improved after the HLA-Y*02:01 sequence was included in the HLA allele dictionary.

Application of High-Throughput Next-Generation Sequencing for HLA Typing on Buccal Extracted DNA: Results from over 10,000 Donor Recruitment Samples

Background

Unambiguous HLA typing is important in hematopoietic stem cell transplantation (HSCT), HLA disease association studies, and solid organ transplantation. However, current molecular typing methods only interrogate the antigen recognition site (ARS) of HLA genes, resulting in many cis-trans ambiguities that require additional typing methods to resolve. Here we report high-resolution HLA typing of 10,063 National Marrow Donor Program (NMDP) registry donors using long-range PCR by next generation sequencing (NGS) approach on buccal swab DNA.

Methods

Multiplex long-range PCR primers amplified the full-length of HLA class I genes (A, B, C) from promotor to 3’ UTR. Class II genes (DRB1, DQB1) were amplified from exon 2 through part of exon 4. PCR amplicons were pooled and sheared using Covaris fragmentation. Library preparation was performed using the Illumina TruSeq Nano kit on the Beckman FX automated platform. Each sample was tagged with a unique barcode, followed by 2×250 bp paired-end sequencing on the Illumina MiSeq. HLA typing was assigned using Omixon Twin software that combines two independent computational algorithms to ensure high confidence in allele calling. Consensus sequence and typing results were reported in Histoimmunogenetics Markup Language (HML) format. All homozygous alleles were confirmed by Luminex SSO typing and exon novelties were confirmed by Sanger sequencing.

Results

Using this automated workflow, over 10,063 NMDP registry donors were successfully typed under high-resolution by NGS. Despite known challenges of nucleic acid degradation and low DNA concentration commonly associated with buccal-based specimens, 97.8% of samples were successfully amplified using long-range PCR. Among these, 98.2% were successfully reported by NGS, with an accuracy rate of 99.84% in an independent blind Quality Control audit performed by the NDMP. In this study, NGS-HLA typing identified 23 null alleles (0.023%), 92 rare alleles (0.091%) and 42 exon novelties (0.042%).

Conclusion

Long-range, unambiguous HLA genotyping is achievable on clinical buccal swab-extracted DNA. Importantly, full-length gene sequencing and the ability to curate full sequence data will permit future interrogation of the impact of introns, expanded exons, and other gene regulatory sequences on clinical outcomes in transplantation.

Determining performance characteristics of an NGS-based HLA typing method for clinical applications

This study presents performance specifications of an in-house developed human leukocyte antigen (HLA) typing assay using next-generation sequencing (NGS) on the Illumina MiSeq platform. A total of 253 samples, previously characterized for HLA-A, -B, -C, -DRB1 and -DQB1 were included in this study, which were typed at high-resolution using a combination of Sanger sequencing, sequence-specific primer (SSP) and sequence-specific oligonucleotide probe (SSOP) technologies and recorded at the two-field level. Samples were selected with alleles that cover a high percentage of HLA specificities in each of five different race/ethnic groups: European, African-American, Asian Pacific Islander, Hispanic and Native American. Sequencing data were analyzed by two software programs, Omixon's target and GenDx's NGSengine. A number of metrics including allele balance, sensitivity, specificity, precision, accuracy and remaining ambiguity were assessed. Data analyzed by the two software systems are shown independently. The majority of alleles were identical in the exonic sequences (third field) with both programs for HLA-A, -B, -C and -DQB1 in 97.7% of allele determinations. Among the remaining discrepant genotype calls at least one of the analysis programs agreed with the reference typing. Upon additional manual analysis 100% of the 2530 alleles were concordant with the reference HLA genotypes; the remaining ambiguities did not exceed 0.8%. The results demonstrate the feasibility and significant benefit of HLA typing by NGS as this technology is highly accurate, eliminates virtually all ambiguities, provides complete sequencing information for the length of the HLA gene and forms the basis for utilizing a single methodology for HLA typing in the immunogenetics labs.

Clinical applications of next-generation sequencing in histocompatibility and transplantation

PURPOSE OF REVIEW

Next-generation sequencing (NGS) can overcome traditional methodological barriers to facilitate detailed studies of large genomes. Here, we summarize recent NGS-based developments in histocompatibility and transplantation, and highlight the dynamic range of clinical applications achievable on this platform.

RECENT FINDINGS

Multiple NGS-based protocols have been established to achieve unambiguous human leukocyte antigen genotyping. These methods are presently engaged to serve the high-throughput demand of large bone marrow registries; however, the scalable nature of NGS makes it an equally attractive technology for select applications within solid organ transplantation. Recently, the exquisite sensitivity of NGS has been leveraged to perform noninvasive allograft monitoring by tracking the dynamics of donor-derived cell-free DNA. Further, NGS-based T-cell receptor and immunoglobulin heavy chain repertoire profiling appear to be useful in clarifying disease-specific diagnoses in certain complex allograft pathology; detecting/quantifying minimal residual disease following allogeneic stem cell transplantation; and tracking donor-reactive T cells to understand the mechanism of tolerance in kidney transplant recipients.

SUMMARY

NGS is superior to classical Sanger sequencing in its throughput, sensitivity, and the ability to provide phase-defined sequence data. These unique properties allow its broad application to diverse areas in clinical transplantation.

HLA-genotyping of clinical specimens using Ion Torrent-based NGS

We have evaluated and validated the NXType™ workflow (One Lambda, Inc.) and the accompanying TypeStream™ software on the Ion Torrent Next Generation Sequencing (NGS) platform using a comprehensive testing panel. The panel consisted of 285 genomic DNA (gDNA) samples derived from four major ethnic populations and contained 59 PT samples and 226 clinical specimens. The total number of alleles from the six loci interrogated by NGS was 3420. This validation panel provided a wide range of HLA sequence variations including many rare alleles, new variants and homozygous alleles. The NXType™ system (reagents and software) was able to correctly genotype the vast majority of these specimens. The concordance rate between SBT-derived genotypes and those generated by TypeStream™ auto-analysis ranged from 99.5% to 99.8% for the HLA-A, B, C, DRB1 and DQB1 loci, and was 98.9% for HLA-DPB1. A strategy for data review was developed that would allow correction of most of the few remaining typing errors. The entire NGS workflow from gDNA amplification to genotype assignment could be completed within 3 working days. Through this validation study, the limitations and shortcomings of the platform, specific assay system, and software algorithm were also revealed for further evaluation and improvement.

The impact of next-generation sequencing technologies on HLA research

In the past decade, the development of next-generation sequencing (NGS) has paved the way for whole-genome analysis in individuals. Research on the human leukocyte antigen (HLA), an extensively studied molecule involved in immunity, has benefitted from NGS technologies. The HLA region, a 3.6-Mb segment of the human genome at 6p21, has been associated with more than 100 different diseases, primarily autoimmune diseases. Recently, the HLA region has received much attention because severe adverse effects of various drugs are associated with particular HLA alleles. Owing to the complex nature of the HLA genes, classical direct sequencing methods cannot comprehensively elucidate the genomic makeup of HLA genes. Thus far, several high-throughput HLA-typing methods using NGS have been developed. In HLA research, NGS facilitates complete HLA sequencing and is expected to improve our understanding of the mechanisms through which HLA genes are modulated, including transcription, regulation of gene expression and epigenetics. Most importantly, NGS may also permit the analysis of HLA-omics. In this review, we summarize the impact of NGS on HLA research, with a focus on the potential for clinical applications.

Novel Approaches and Technologies in Molecular HLA Typing

The invention of the Polymerase Chain Reaction (PCR) has revolutionized molecular biology enabling gene isolation and characterization in hours rather than days. Scientists working in transplant diagnostics have proven to be pioneers in adapting this molecular technique to the clinical needs of histocompatibility testing. This chapter describes a number of novel genotyping technologies which have been used to address the challenges posed by genetic diversity seen in the extensive polymorphism in HLA genes. These novel approaches include single-stranded and duplex conformational analyses, real-time PCR, microarray hybridization, RNA-based sequencing, and the present day Next Generation Sequencing. The chapter concludes with a brief look at a possible next, Next Generation Sequencing system.

HLA typing by next-generation sequencing - getting closer to reality

Next generation sequencing (NGS) denotes novel sequencing technologies that enable the generation of a large number of clonal sequences in a single sequencing run. NGS was initially introduced for whole genome sequencing and for quantitation of viral variants or genetic mutations in tumor tissues; more recently, the potential for high resolution HLA typing and high throughput analyses has been explored. It became clear that the complexity of the HLA system implicates new challenges, especially for bioinformatics. From an economical point of view, NGS is becoming increasingly attractive for HLA typing laboratories currently relying on Sanger based sequencing. Realizing the full potential of NGS will require the development of specifically adapted typing strategies and software algorithms. In the present review, three laboratories that were among the first to perform HLA-typing using different NGS platforms, the Roche 454, the Illumina Miseq and the Ion Torrent system, respectively, give an overview of these applications and point out advantages and limitations.

Extended blood group molecular typing and next-generation sequencing

Several high-throughput multiplex blood group molecular typing platforms have been developed to predict blood group antigen phenotypes. These molecular systems support extended donor/patient matching by detecting commonly encountered blood group polymorphisms as well as rare alleles that determine the expression of blood group antigens. Extended molecular typing of a large number of blood donors by high-throughput platforms can increase the likelihood of identifying donor red blood cells that match those of recipients. This is especially important in the management of multiply-transfused patients who may have developed several alloantibodies. Nevertheless, current molecular techniques have limitations. For example, they detect only predefined genetic variants. In contrast, target enrichment next-generation sequencing (NGS) is an emerging technology that provides comprehensive sequence information, focusing on specified genomic regions. Target enrichment NGS is able to assess genetic variations that cannot be achieved by traditional Sanger sequencing or other genotyping platforms. Target enrichment NGS has been used to detect both known and de novo genetic polymorphisms, including single-nucleotide polymorphisms, indels (insertions/deletions), and structural variations. This review discusses the methodology, advantages, and limitations of the current blood group genotyping techniques and describes various target enrichment NGS approaches that can be used to develop an extended blood group genotyping assay system.

OptiType: precision HLA typing from next-generation sequencing data

Motivation: The human leukocyte antigen (HLA) gene cluster plays a crucial role in adaptive immunity and is thus relevant in many biomedical applications. While next-generation sequencing data are often available for a patient, deducing the HLA genotype is difficult because of substantial sequence similarity within the cluster and exceptionally high variability of the loci. Established approaches, therefore, rely on specific HLA enrichment and sequencing techniques, coming at an additional cost and extra turnaround time.

Result: We present OptiType, a novel HLA genotyping algorithm based on integer linear programming, capable of producing accurate predictions from NGS data not specifically enriched for the HLA cluster. We also present a comprehensive benchmark dataset consisting of RNA, exome and whole-genome sequencing data. OptiType significantly outperformed previously published in silico approaches with an overall accuracy of 97% enabling its use in a broad range of applications.

Contact:szolek@informatik.uni-tuebingen.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Resolution of ambiguous HLA genotyping in korean by multi-group-specific sequence-based typing

PURPOSE

To evaluate a multi-group-specific sequence-based typing (SBT) method for resolving ambiguous results from human leukocyte antigen (HLA) genotyping.

MATERIALS AND METHODS

A total of 50 samples that showed ambiguous genotypes for at least two HLA loci from HLA-A, -B, -C and -DRB1 by the conventional SBT assay were evaluated using a new SBT test, the AVITA plus assay. The most likely HLA genotypes for the respective samples considering allele frequencies in Korean were concordant between the AVITA and conventional SBT assays.

RESULTS

An average of 3.3 loci among the HLA-A, -B, -C and -DRB1 loci per sample gave results with two or more possible allele combinations with the conventional SBT, and 48 (96.0%) out of 50 showed reduced numbers of possible genotypes for at least one HLA locus with the AVITA. A total of 41, 43, 42, and 38 cases among the 50 samples showed ambiguous results for HLA-A, -B, -C, and -DRB1 typing by the conventional SBT, respectively. The average numbers of possible allele combinations for the respective four HLA loci were 8.2, 6.7, 5.9, and 3.2, and they were reduced to 1.5, 2.2, 4.4, and 1.8, respectively, by the AVITA. Ambiguity was resolved by the AVITA in 33 (80.5%), 31 (72.1%), 17 (40.5%) and 28 (73.7%) samples among the ambiguous cases from the conventional SBT for HLA-A, -B, -C, and -DRB1 typing, respectively.

CONCLUSION

The multi-group-specific SBT method considerably reduced the number of ambiguous results, and thus may be useful for accurate HLA typing in clinical laboratories.

Allo-reactivity of mesenchymal stem cells in rhesus macaques is dose and haplotype dependent and limits durable cell engraftment in vivo

The emerging paradigm that MSCs are immune privileged has fostered the use of “off-the-shelf” allogeneic MSC-based therapies in human clinical trials. However, this approach ignores studies in experimental animals wherein transplantation of MSCs across MHC boundaries elicits measurable allo-immune responses. To determine if MSCs are hypo-immunogeneic, we characterized the immune response in rhesus macaques following intracranial administration of allogeneic vs. autologous MSCs. This analysis revealed unambiguous evidence of productive allo-recognition based on expansion of NK, B and T cell subsets in peripheral blood and detection of allo-specific antibodies in animals administered allogeneic but not autologous MSCs. Moreover, the degree of MHC class I and II mismatch between the MSC donor and recipient significantly influenced the magnitude and nature of the allo-immune response. Consistent with these findings, real-time PCR analysis of brain tissue from female recipients administered varying doses of male, allogeneic MSCs revealed a significant inverse correlation between MSC engraftment levels and cell dose. Changes in post-transplant neutrophil and lymphocyte counts also correlated with dose and were predictive of overall MSC engraftment levels. However, secondary antigen challenge failed to elicit a measurable immune response in allogeneic recipients. Finally, extensive behavior testing of animals revealed no main effect of cell dose on motor skills, social development, or temperament. Collectively, these data indicate that allogeneic MSCs are weakly immunogenic when transplanted across MHC boundaries in rhesus macaques and this negatively impacts durable engraftment levels. Therefore the use of unrelated donor MSCs should be carefully evaluated in human patients.

Cost-efficient high-throughput HLA typing by MiSeq amplicon sequencing

Background

A close match of the HLA alleles between donor and recipient is an important prerequisite for successful unrelated hematopoietic stem cell transplantation. To increase the chances of finding an unrelated donor, registries recruit many hundred thousands of volunteers each year. Many registries with limited resources have had to find a trade-off between cost and resolution and extent of typing for newly recruited donors in the past. Therefore, we have taken advantage of recent improvements in NGS to develop a workflow for low-cost, high-resolution HLA typing.

Results

We have established a straightforward three-step workflow for high-throughput HLA typing: Exons 2 and 3 of HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1 are amplified by PCR on Fluidigm Access Array microfluidic chips. Illumina sequencing adapters and sample specific tags are directly incorporated during PCR. Upon pooling and cleanup, 384 samples are sequenced in a single Illumina MiSeq run. We developed “neXtype” for streamlined data analysis and HLA allele assignment. The workflow was validated with 1140 samples typed at 6 loci. All neXtype results were concordant with the Sanger sequences, demonstrating error-free typing of more than 6000 HLA loci. Current capacity in routine operation is 12,000 samples per week.

Conclusions

The workflow presented proved to be a cost-efficient alternative to Sanger sequencing for high-throughput HLA typing. Despite the focus on cost efficiency, resolution exceeds the current standards of Sanger typing for donor registration.

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.

Major histocompatibility complex genomics and human disease

Over several decades, various forms of genomic analysis of the human major histocompatibility complex (MHC) have been extremely successful in picking up many disease associations. This is to be expected, as the MHC region is one of the most gene-dense and polymorphic stretches of human DNA. It also encodes proteins critical to immunity, including several controlling antigen processing and presentation. Single-nucleotide polymorphism genotyping and human leukocyte antigen (HLA) imputation now permit the screening of large sample sets, a technique further facilitated by high-throughput sequencing. These methods promise to yield more precise contributions of MHC variants to disease. However, interpretation of MHC-disease associations in terms of the functions of variants has been problematic. Most studies confirm the paramount importance of class I and class II molecules, which are key to resistance to infection. Infection is likely driving the extreme variation of these genes across the human population, but this has been difficult to demonstrate. In contrast, many associations with autoimmune conditions have been shown to be specific to certain class I and class II alleles. Interestingly, conditions other than infections and autoimmunity are also associated with the MHC, including some cancers and neuropathies. These associations could be indirect, owing, for example, to the infectious history of a particular individual and selective pressures operating at the population level.

HLA typing from RNA-seq data using hierarchical read weighting [corrected]

Correctly matching the HLA haplotypes of donor and recipient is essential to the success of allogenic hematopoietic stem cell transplantation. Current HLA typing methods rely on targeted testing of recognized antigens or sequences. Despite advances in Next Generation Sequencing, general high throughput transcriptome sequencing is currently underutilized for HLA haplotyping due to the central difficulty in aligning sequences within this highly variable region. Here we present the method, HLAforest, that can accurately predict HLA haplotype by hierarchically weighting reads and using an iterative, greedy, top down pruning technique. HLAforest correctly predicts >99% of allele group level (2 digit) haplotypes and 93% of peptide-level (4 digit) haplotypes of the most diverse HLA genes in simulations with read lengths and error rates modeling currently available sequencing technology. The method is very robust to sequencing error and can predict 99% of allele-group level haplotypes with substitution rates as high as 8.8%. When applied to data generated from a trio of cell lines, HLAforest corroborated PCR-based HLA haplotyping methods and accurately predicted 16/18 (89%) major class I genes for a daughter-father-mother trio at the peptide level. Major class II genes were predicted with 100% concordance between the daughter-father-mother trio. In fifty HapMap samples with paired end reads just 37 nucleotides long, HLAforest predicted 96.5% of allele group level HLA haplotypes correctly and 83% of peptide level haplotypes correctly. In sixteen RNAseq samples with limited coverage across HLA genes, HLAforest predicted 97.7% of allele group level haplotypes and 85% of peptide level haplotypes correctly.

ATHLATES: accurate typing of human leukocyte antigen through exome sequencing

Human leukocyte antigen (HLA) typing at the allelic level can in theory be achieved using whole exome sequencing (exome-seq) data with no added cost but has been hindered by its computational challenge. We developed ATHLATES, a program that applies assembly, allele identification and allelic pair inference to short read sequences, and applied it to data from Illumina platforms. In 15 data sets with adequate coverage for HLA-A, -B, -C, -DRB1 and -DQB1 genes, ATHLATES correctly reported 74 out of 75 allelic pairs with an overall concordance rate of 99% compared with conventional typing. This novel approach should be broadly applicable to research and clinical laboratories.

Routine performance and errors of 454 HLA exon sequencing in diagnostics

Background

Next-generation sequencing (NGS) has changed genomics significantly. More and more applications strive for sequencing with different platforms. Now, in 2012, after a decade of development and evolution, NGS has been accepted for a variety of research fields. Determination of sequencing errors is essential in order to follow next-generation sequencing beyond research use only. This study describes the overall 454 system performance of using multiple GS Junior runs with an in-house established and validated diagnostic assay for human leukocyte antigen (HLA) exon sequencing. Based on this data, we extracted, evaluated and characterized errors and variants of 60 HLA loci per run with respect to their adjacencies.

Results

We determined an overall error rate of 0.18% in a total of 118,484,408 bases. 31.3% of all reads analyzed (n=349,503) contain one or more errors. The largest group are deletions that account for 50% of the errors. Incorrect bases are not distributed equally along sequences and tend to be more frequent at sequence ends. Certain sequence positions in the middle or at the beginning of the read accumulate errors. Typically, the corresponding quality score at the actual error position is lower than the adjacent scores.

Conclusions

Here we present the first error assessment in a human next-generation sequencing diagnostics assay in an amplicon sequencing approach. Improvements of sequence quality and error rate that have been made over the years are evident and it is shown that both have now reached a level where diagnostic applications become feasible. Our presented data are better than previously published error rates and we can confirm and quantify the often described relation of homopolymers and errors. Nevertheless, a certain depth of coverage is needed, in particular with challenging areas of the sequencing target. Furthermore, the usage of error correcting tools is not essential but might contribute towards the capacity and efficiency of a sequencing run.

Diagnostic applications of next generation sequencing in immunogenetics and molecular oncology

SUMMARY

With the introduction of the next generation sequencing (NGS) technologies, remarkable new diagnostic applications have been established in daily routine. Implementation of NGS is challenging in clinical diagnostics, but definite advantages and new diagnostic possibilities make the switch to the technology inevitable. In addition to the higher sequencing capacity, clonal sequencing of single molecules, multiplexing of samples, higher diagnostic sensitivity, workflow miniaturization, and cost benefits are some of the valuable features of the technology. After the recent advances, NGS emerged as a proven alternative for classical Sanger sequencing in the typing of human leukocyte antigens (HLA). By virtue of the clonal amplification of single DNA molecules ambiguous typing results can be avoided. Simultaneously, a higher sample throughput can be achieved by tagging of DNA molecules with multiplex identifiers and pooling of PCR products before sequencing. In our experience, up to 380 samples can be typed for HLA-A, -B, and -DRB1 in high-resolution during every sequencing run. In molecular oncology, NGS shows a markedly increased sensitivity in comparison to the conventional Sanger sequencing and is developing to the standard diagnostic tool in detection of somatic mutations in cancer cells with great impact on personalized treatment of patients.

Rapid, scalable and highly automated HLA genotyping using next-generation sequencing: a transition from research to diagnostics

BACKGROUND

Human leukocyte antigen matching at allelic resolution is proven clinically significant in hematopoietic stem cell transplantation, lowering the risk of graft-versus-host disease and mortality. However, due to the ever growing HLA allele database, tissue typing laboratories face substantial challenges. In light of the complexity and the high degree of allelic diversity, it has become increasingly difficult to define the classical transplantation antigens at high-resolution by using well-tried methods. Thus, next-generation sequencing is entering into diagnostic laboratories at the perfect time and serving as a promising tool to overcome intrinsic HLA typing problems. Therefore, we have developed and validated a scalable automated HLA class I and class II typing approach suitable for diagnostic use.

RESULTS

A validation panel of 173 clinical and proficiency testing samples was analysed, demonstrating 100% concordance to the reference method. From a total of 1,273 loci we were able to generate 1,241 (97.3%) initial successful typings. The mean ambiguity reduction for the analysed loci was 93.5%. Allele assignment including intronic sequences showed an improved resolution (99.2%) of non-expressed HLA alleles.

CONCLUSION

We provide a powerful HLA typing protocol offering a short turnaround time of only two days, a fully integrated workflow and most importantly a high degree of typing reliability. The presented automated assay is flexible and can be scaled by specific primer compilations and the use of different 454 sequencing systems. The workflow was successfully validated according to the policies of the European Federation for Immunogenetics. Next-generation sequencing seems to become one of the new methods in the field of Histocompatibility.

Clustering and alignment of polymorphic sequences for HLA-DRB1 genotyping

Located on Chromosome 6p21, classical human leukocyte antigen genes are highly polymorphic. HLA alleles associate with a variety of phenotypes, such as narcolepsy, autoimmunity, as well as immunologic response to infectious disease. Moreover, high resolution genotyping of these loci is critical to achieving long-term survival of allogeneic transplants. Development of methods to obtain high resolution analysis of HLA genotypes will lead to improved understanding of how select alleles contribute to human health and disease risk. Genomic DNAs were obtained from a cohort of n = 383 subjects recruited as part of an Ulcerative Colitis study and analyzed for HLA-DRB1. HLA genotypes were determined using sequence specific oligonucleotide probes and by next-generation sequencing using the Roche/454 GSFLX instrument. The Clustering and Alignment of Polymorphic Sequences (CAPSeq) software application was developed to analyze next-generation sequencing data. The application generates HLA sequence specific 6-digit genotype information from next-generation sequencing data using MUMmer to align sequences and the R package diffusionMap to classify sequences into their respective allelic groups. The incorporation of Bootstrap Aggregating, Bagging to aid in sorting of sequences into allele classes resulted in improved genotyping accuracy. Using Bagging iterations equal to 60, the genotyping results obtained using CAPSeq when compared with sequence specific oligonucleotide probe characterized 4-digit genotypes exhibited high rates of concordance, matching at 759 out of 766 (99.1%) alleles.

New pharmacogenetic test for detecting an HLA-A*31: 01 allele using the InvaderPlus assay

BACKGROUND & AIMS

Carbamazepine (CBZ) is widely used for the treatment of epilepsy and other neurological disorders. However, 3-5% of CBZ-treated individuals suffer from cutaneous adverse drug reactions (cADRs). Recently, in a genome-wide association study, HLA-A*31:01 has been reported to be a strong genetic marker for CBZ-induced cADRs in both Japanese and European populations. As most of the available methods for HLA genotyping are laborious, the development of a simple and rapid genotyping method for HLA-A*31:01 is desirable from the viewpoint of a clinical pharmacogenetic test.

METHODS

More than 1700 sequences for HLA-A alleles were obtained from the MHC database of the National Center for Biotechnology Information (dbMHC). Several HLA-A*31:01-discriminating single-nucleotide polymorphisms were selected. These SNPs were used for sequence-specific primer PCR (SSP-PCR) and for the target site of the Invader reaction. By combining SSP-PCR with a target-specific Invader reaction, we designed two sets of primers/probes for HLA-A*31:01 allele detection. The performance of both sets was evaluated using 90 Asian HapMap samples. Further evaluation was carried out using another 376 Japanese samples and 90 CEU (European) and 90 YRI (African) HapMap samples.

RESULTS

Our assay specifically detected an HLA-A*31:01 allele in a total of 466 individuals of the Asian population. Furthermore, the assay correctly identified HLA-A*31:01-positive carriers from the CEU and the YRI population, respectively, implying that the assay has potential for application to other ethnic groups.

CONCLUSION

We developed a new HLA-A*31:01-detecting method by a combination of SSP-PCR with target-specific InvaderPlus technology. As our assay is rapid and accurate, it is hoped that this method will be used in a pharmacogenetic test in a clinical setting to avoid CBZ-induced cADRs.