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

A novel framework for human leukocyte antigen (HLA) genotyping using probe capture-based targeted next-generation sequencing and computational analysis

Human leukocyte antigen (HLA) genes play pivotal roles in numerous immunological applications. Given the immense number of polymorphisms, achieving accurate high-throughput HLA typing remains challenging. This study aimed to harness the human pan-genome reference consortium (HPRC) resources as a potential benchmark for HLA reference materials. We meticulously annotated specific four field-resolution alleles for 11 HLA genes (HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1, -DRB3, -DRB4 and -DRB5) from 44 high-quality HPRC personal genome assemblies. For sequencing, we crafted HLA-specific probes and conducted capture-based targeted sequencing of the genomic DNA of the HPRC cohort, ensuring focused and comprehensive coverage of the HLA region of interest. We used publicly available short-read whole-genome sequencing (WGS) data from identical samples to offer a comparative perspective. To decipher the vast amount of sequencing data, we employed seven distinct software tools: OptiType, HLA-VBseq, HISAT genotype, SpecHLA, T1K, QzType, and DRAGEN. Each tool offers unique capabilities and algorithms for HLA genotyping, allowing comprehensive analysis and validation of the results. We then compared these results with benchmarks derived from personal genome assemblies. Our findings present a comprehensive four-field-resolution HLA allele annotation for 44 HPRC samples. Significantly, our innovative targeted next-generation sequencing (NGS) approach for HLA genes showed superior accuracy compared with conventional short-read WGS. An integrated analysis involving QzType, T1K, and DRAGEN was developed, achieving 100% accuracy for all 11 HLA genes. In conclusion, our study highlighted the combination of targeted short-read sequencing and astute computational analysis as a robust approach for HLA genotyping. Furthermore, the HPRC cohort has emerged as a valuable assembly-based reference in this realm.

Development of an ARMS multiplex real-time PCR assay for the detection of HLA-B*13:01 genotype by detecting highly specific SNPs

OBJECTIVES

HLA-B*13:01 was strongly associated with Dapsone Hypersensitivity Syndrome (DHS). This study aimed to develop and validate a rapid and economical method for HLA-B*13:01 genotyping.

METHODS

Two tubes multiplex real-time PCR detection system comprising amplification refractory mutation system primers and TaqMan probes was established for HLA-B*13:01 genotyping. Sequence-based typing was applied to validate the accuracy of the assay.

RESULTS

The accuracy of the assay was 100% for HLA-B*13:01 genotyping. The detection limit of the new method was 0.025 ng DNA. The positive rate of HLA-B*13:01 in the Bouyei (20%, n = 50) populations was significantly higher than that in the Uighur population (4%, n = 100), Han (4.5%, n = 200), and Tibetan (1%, n = 100) ( P  < 0.05).

CONCLUSION

The proposed method is rapid and reliable for HLA-B*13:01 screening in a clinical setting.

A single-tube multiplex real-time PCR for HLA-B*38:02 genotype by detecting highly specific SNPs

Purpose: HLA-B*38:02 is closely related to carbimazole/methimazole-induced agranulocytosis. This study aimed to develop and validate a rapid and economical method for HLA-B*38:02 genotyping. Methods: A single-tube multiplex real-time PCR detection system comprising amplification refractory mutation system primers and TaqMan probes was established for HLA-B*38:02 genotyping. Sequence-based typing was applied to validate the accuracy of the assay. Results: The accuracy of the assay was 100% for HLA-B*38:02 genotyping. The detection limit of the new method was 0.05 ng DNA. The positive rate of HLA-B*38:02 in the Han (8%, n = 100), Bouyei (17.8%, n = 90) and Tibetan (12.7%, n = 110) populations was significantly higher than that in the Uighur population (1%, n = 100) (p < 0.05). Conclusion: The proposed method is rapid and reliable for HLA-B*38:02 screening in a clinical setting.

rs144012689 is a highly specific representative marker of HLA-B*15:02 in the Chinese population

Objective: To develop an accurate and rapid real-time PCR technique for HLA-B*15:02 genotyping and investigate HLA-B*15:02 allele frequency in four ethnic populations in China. Materials & methods: Based on the highly specific representative markers, a real-time PCR assay was developed for HLA-B*15:02 genotyping, and HLA-B*15:02 allele frequencies were screened in four ethnic populations of China. Sequence-based typing was used to validate the accuracy of the assay. Results: The sensitivity and specificity of the assay were 100%, and the detection limit was 0.2 ng. The frequency of HLA-B*15:02 alleles distributed in the Bouyei population was significantly higher than in the Han group (p < 0.01). Neither the Tibetan nor the Uyghur population carried the HLA-B*15:02 haplotype. Conclusion: The authors developed an accurate HLA-B*15:02 genotyping method for evaluating the risk of adverse drug reactions induced by carbamazepine in various ethnic populations in China.

Recent Advances of Human Leukocyte Antigen (HLA) Typing Technology Based on High-Throughput Sequencing

The major histocompatibility complex (MHC) in humans is a genetic region consisting of cell surface proteins located on the short arm of chromosome 6. This is also known as the human leukocyte antigen (HLA) region. The HLA region consists of genes that exhibit complex genetic polymorphisms, and are extensively involved in immune responses. Each individual has a unique set of HLAs. Donor-recipient HLA allele matching is an important factor for organ transplantation. Therefore, an established rapid and accurate HLA typing technology is instrumental to preventing graft-verses-host disease (GVHD) in organ recipients. As of recent, high-throughput sequencing has allowed for an increase read length and higher accuracy and throughput, thus achieving complete and high-resolution full-length typing. With more advanced nanotechnology used in high-throughput sequencing, HLA typing is more widely used in third-generation single-molecule sequencing. This review article summarizes some of the most widely used sequencing typing platforms and evaluates the latest developments in HLA typing kits and their clinical applications.

SARS-CoV-2 Host Immunogenetic Biomarkers

SARS-CoV-2 causes generally mild symptoms, with approximately 10-20% of cases progressing to severe disease. The pathophysiologic mechanisms by which SARS-CoV-2 causes severe disease are largely unknown. Data have indicated the involvement of different immunogenetic markers such as HLA, T, and B cells, to be associated with disease outcome. This has led to interest in these genes as potential biomarkers of SARS-CoV-2 susceptibility and for predicting prognosis and response to vaccines and other therapeutic strategies. In this chapter, we discussed outline protocols for characterizing these potential biomarkers and methods for identifying SARS-CoV-2 biomarkers using the Luminex^® 100/200 technology and next-generation sequencing.

Next-generation sequencing and clinical histocompatibility testing

Histocompatibility testing is essential for donor identification and risk assessment in solid organ and hematopoietic stem cell transplant. Additionally, it is useful for identifying donor specific alleles for monitoring donor specific antibodies in post-transplant patients. Next-generation sequence (NGS) based human leukocyte antigen (HLA) typing has improved many aspects of histocompatibility testing in hematopoietic stem cell and solid organ transplant. HLA disease association testing and research has also benefited from the advent of NGS technologies. In this review we discuss the current impact and future applications of NGS typing on clinical histocompatibility testing for transplant and non-transplant purposes.

Identification of 8-Digit HLA-A, -B, -C, and -DRB1 Allele and Haplotype Frequencies in Koreans Using the One Lambda AllType Next-Generation Sequencing Kit

Background

Recent studies have successfully implemented next-generation sequencing (NGS) in HLA typing. We performed HLA NGS in a Korean population to estimate HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies up to an 8-digit resolution, which might be useful for an extended application of HLA results.

Methods

A total of 128 samples collected from healthy unrelated Korean adults, previously subjected to Sanger sequencing for 6-digit HLA analysis, were used. NGS was performed for HLA-A, -B, -C, and -DRB1 using the AllType NGS kit (One Lambda, West Hills, CA, USA), Ion Torrent S5 platform (Thermo Fisher Scientific, Waltham, MA, USA), and Type Steam Visual NGS analysis software (One Lambda).

Results

Eight HLA alleles showed frequencies of ≥10% in the Korean population, namely, A*24:02:01:01 (19.5%), A*33:03:01 (15.6%), A*02:01:01:01 (14.5%), A*11:01:01:01 (13.3%), B*15:01:01:01 (10.2%), C*01:02:01 (19.9%), C*03:04:01:02 (11.3%), and DRB1*09:01:02 (10.2%). Nine previous 6-digit HLA alleles were further identified as two or more 8-digit HLA alleles. Of these, eight alleles (A*24:02:01, B*35:01:01, B*40:01:02, B*55:02:01, B*58:01:01, C*03:02:02, C*07:02:01, and DRB1*07:01:01) were identified as two 8-digit HLA alleles, and one allele (B*51:01:01) was identified as three 8-digit HLA alleles. The most frequent four-loci haplotype was HLA-A*33:03:01-B*44:03:01:01-C*14:03-DRB1*13:02:01.

Conclusions

We identified 8-digit HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies in a healthy Korean population using NGS. These new data can be used as a representative Korean data for further disease-related HLA type analysis.

A high-affinity human TCR-like antibody detects celiac disease gluten peptide-MHC complexes and inhibits T cell activation

Antibodies specific for peptides bound to human leukocyte antigen (HLA) molecules are valuable tools for studies of antigen presentation and may have therapeutic potential. Here, we generated human T cell receptor (TCR)-like antibodies toward the immunodominant signature gluten epitope DQ2.5-glia-α2 in celiac disease (CeD). Phage display selection combined with secondary targeted engineering was used to obtain highly specific antibodies with picomolar affinity. The crystal structure of a Fab fragment of the lead antibody 3.C11 in complex with HLA-DQ2.5:DQ2.5-glia-α2 revealed a binding geometry and interaction mode highly similar to prototypic TCRs specific for the same complex. Assessment of CeD biopsy material confirmed disease specificity and reinforced the notion that abundant plasma cells present antigen in the inflamed CeD gut. Furthermore, 3.C11 specifically inhibited activation and proliferation of gluten-specific CD4⁺ T cells in vitro and in HLA-DQ2.5 humanized mice, suggesting a potential for targeted intervention without compromising systemic immunity.

Allele and haplotype frequencies of human leukocyte antigen-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 by next generation sequencing-based typing in Koreans in South Korea

Allele frequencies and haplotype frequencies of HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 have been rarely reported in South Koreans using unambiguous, phase-resolved next generation DNA sequencing. In this study, HLA typing of 11 loci in 173 healthy South Koreans were performed using next generation DNA sequencing with long-range PCR, TruSight^® HLA v2 kit, Illumina MiSeqDx platform system, and Assign^™ for TruSight^™ HLA software. Haplotype frequencies were calculated using the PyPop software. Direct counting methods were used to investigate the association with DRB1 for samples with only one copy of a particular secondary DRB locus. We compared these allele types with the ambiguous allele combinations of the IPD-IMGT/HLA database. We identified 20, 40, 26, 31, 19, 16, 4, and 16 alleles of HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1, respectively. The number of HLA-DRB3/4/5 alleles was 4, 5, and 3, respectively. The haplotype frequencies of most common haplotypes were as follows: A*33:03:01-B*44:03:01-C*14:03-DRB1*13:02:01-DQB1*06:04:01-DPB1*04:01:01 (2.89%), A*33:03:01-B*44:03:01-C*14:03 (4.91%), DRB1*08:03:02-DQA1*01:03:01-DQB1*06:01:01-DPA1*02:02:02-DPB1*05:01:01 (5.41%), DRB1*04:05:01-DRB4*01:03:01 (12.72%), DQA1*01:03:01-DQB1*06:01:01 (13.01%), and DPA1*02:02:02-DPB1*05:01:01 (30.83%). In samples with only one copy of a specific secondary DRB locus, we examined its association with DRB1. We, thus, resolved 10 allele ambiguities in HLA-B, -C (each exon 2+3), -DRB1, -DQB1, -DQA1, and -DPB1 (each exon 2) of the IPD-IMGT/HLA database. Korean population was geographically close to Japanese and Han Chinese populations in the genetic distances by multidimensional scaling (MDS) plots. The information obtained by HLA typing of the 11 extended loci by next generation sequencing may be useful for more exact diagnostic tests on various transplantations and the genetic population relationship studies in South Koreans.

Concurrent typing of over 4000 samples by long-range PCR amplicon-based NGS and rSSO revealed the need to verify NGS typing for HLA allelic dropouts

Hematopoietic stem cell transplantation (HSCT) from HLA-matched donors significantly decreases the risks of graft-rejection and graft-versus-host disease. Long-range PCR- amplicon-based next-generation sequencing (NGS) is increasingly used as a standalone method in clinical laboratories to determine HLA compatibility for HSCT and solid-organ transplantation. We hypothesized that an allelic dropout is a frequent event in the long-range PCR amplicon-based NGS HLA typing method. To test the hypothesis, we typed 4,006 samples concurrently using a commercially available long-range PCR amplicon-based NGS-typing and short exon-specific amplicon-based reverse sequence-specific oligonucleotide (rSSO) methods. The concordance between the NGS and rSSO typing results was 100% at HLA-A, -B, -C, -DRB1, -DRB3, -DRB5, -DQA1, DPA1 loci. However, 4.5% of the samples (179/4006) showed allelic-dropouts at one of the other three loci: HLA-DRB4 (3.9%), HLA-DPB1 (0.4%), and HLA-DQB1*(0.15%). The allelic-dropouts are not associated with specific haplotypes, and some dropouts can be reagent lot-specific. Although DRB1-DRB3/4/5-DQB1 linkages help to diagnose these allelic-dropouts in some cases, the rSSO typing was crucial to identify the dropouts in DQB1 and DPB1 loci. These results uncover the critical limitations of using long-range PCR amplicon-based NGS as a standalone method in clinical histocompatibility laboratories and advocate the need for strategies to diagnose and resolve allelic-dropouts.

HLA class I and II associations with common enteric pathogens in the first year of life

BACKGROUND

genetic susceptibility to infection is mediated by numerous host factors, including the highly diverse, classical human leukocyte antigen (HLA) genes, which are critical genetic determinants of immunity. We systematically evaluated the effect of HLA alleles and haplotypes on susceptibility to 12 common enteric infections in children during the first year of life in an urban slum of Dhaka, Bangladesh.

METHODS

a birth cohort of 601 Bangladeshi infants was prospectively monitored for diarrhoeal disease. Each diarrhoeal stool sample was analyzed for enteric pathogens by multiplex TaqMan Array Card (TAC). High resolution genotyping of HLA class I (A and B) and II (DRB1, DQA1, and DQB1) genes was performed by next-generation sequencing. We compared the frequency of HLA alleles and haplotypes between infected and uninfected children.

FINDINGS

we identified six individual allele associations and one five-locus haplotype association. One allele was associated with protection: A*24:02 - EAEC. Five alleles were associated with increased risk: A*24:17 - typical EPEC, B*15:01 - astrovirus, B*38:02 - astrovirus, B*38:02 - Cryptosporidium and DQA1*01:01 - Cryptosporidium. A single five-locus haplotype was associated with protection: A*11:01~B*15:02~DRB1*12:02~DQA1*06:01~DQB1*03:01- adenovirus 40/41.

INTERPRETATION

our findings suggest a role for HLA in susceptibility to early enteric infection for five pathogens. Understanding the genetic contribution of HLA in susceptibility has important implications in vaccine design and understanding regional differences in incidence of enteric infection.

FUNDING

this research was supported by the National Institute of Health (NIH) and the Bill and Melinda Gates Foundation.

Next-generation sequencing technologies: An overview

Since the days of Sanger sequencing, next-generation sequencing technologies have significantly evolved to provide increased data output, efficiencies, and applications. These next generations of technologies can be categorized based on read length. This review provides an overview of these technologies as two paradigms: short-read, or "second-generation," technologies, and long-read, or "third-generation," technologies. Herein, short-read sequencing approaches are represented by the most prevalent technologies, Illumina and Ion Torrent, and long-read sequencing approaches are represented by Pacific Biosciences and Oxford Nanopore technologies. All technologies are reviewed along with reported advantages and disadvantages. Until recently, short-read sequencing was thought to provide high accuracy limited by read-length, while long-read technologies afforded much longer read-lengths at the expense of accuracy. Emerging developments for third-generation technologies hold promise for the next wave of sequencing evolution, with the co-existence of longer read lengths and high accuracy.

New challenges, new opportunities: Next generation sequencing and its place in the advancement of HLA typing

The Human Leukocyte Antigen (HLA) system has a critical role in immunorecognition, transplantation, and disease association. Early typing techniques provided the foundation for genotyping methods that revealed HLA as one of the most complex, polymorphic regions of the human genome. Next Generation Sequencing (NGS), the latest molecular technology introduced in clinical tissue typing laboratories, has demonstrated advantages over other established methods. NGS offers high-resolution sequencing of entire genes in time frames and price points considered unthinkable just a few years ago, contributing a wealth of data informing histocompatibility assessment and standards of clinical care. Although the NGS platforms share a high-throughput massively parallel processing model, differing chemistries provide specific strengths and weaknesses. Research-oriented Third Generation Sequencing and related advances in bioengineering continue to broaden the future of NGS in clinical settings. These diverse applications have demanded equally innovative strategies for data management and computational bioinformatics to support and analyze the unprecedented volume and complexity of data generated by NGS. We discuss some of the challenges and opportunities associated with NGS technologies, providing a comprehensive picture of the historical developments that paved the way for the NGS revolution, its current state and future possibilities for HLA typing.

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.

Identification of a single nucleotide deletion in the novel HLA-DQB1*06:379N allele, detected by Polymerase Chain Reaction-Sequence Based Typing but not by Next Generation Sequencing

In this report we describe the characterization of a novel null HLA-DQB1 allele, detected by polymerase chain reaction-sequence-based typing but not by next-generation sequencing (NGS). The new allele, HLA-DQB1*06:379N, was discovered in an Italian patient with acute myeloid leukemia and also in one of her healthy siblings. The new allele is largely homologous to DQB1*06:02:01:01 with a T deletion in exon 2, in codon 11, which causes a frameshift and the formation of an unusual and premature TGA STOP in codon 27. This case report underlines the importance of not removing Sanger method sequencing from routine use for high-resolution HLA typing, but to maintain it together with NGS technology.

Next-generation sequencing of 11 HLA loci in a large dengue vaccine cohort from the Philippines

HLA genotyping by next-generation sequencing (NGS) has evolved with significant advancements in the last decade. Here we describe full-length HLA genotyping of 11 loci in 612 individuals comprising a dengue vaccine cohort from Cebu province in the Philippines. The multi-locus individual tagging NGS (MIT-NGS) method that we developed initially for genotyping 4-6 loci in one MiSeq run was expanded to 11 loci including HLA-A, B, C, DPA1, DPB1, DQA1, DQB1, DRB1, and DRB3/4/5. This change did not affect the overall coverage or depth of the sequencing reads. HLA alleles with frequencies greater than 10% were A*11:01:01, A*24:02:01, A*24:07:01, A*34:01:01, B*38:02:01, B*15:35, B*35:05:01, C*07:02:01, C*04:01:01, DPA1*02:02:02, DPB1*05:01:01, DPB1*01:01:01, DQA1*01:02:01, DQA1*06:01:01, DQB1*05:02:01, DQB1*03:01:01, DRB1*15:02:01, DRB1*12:02:01, DRB3*03:01:03, DRB4*01:03:01, and DRB5*01:01:01. Improvements in sequencing library preparation provide uniform and even coverage across all exons and introns. This has led to a marked reduction in allele imbalance and dropout. Furthermore, including more loci, such as DRB3/4/5, decreases cross-mapping and incorrect allele assignment at the DRB1 locus. The increased number of loci sequenced for each sample does not reduce the number of samples that can be multiplexed on a single MiSeq run and is therefore more cost-efficient. We believe that such improvements will help HLA genotyping by NGS to gain momentum over other conventional methods by increasing confidence in the calls.

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.

Concordance between predicted HLA type using next generation sequencing data generated for non-HLA purposes and clinical HLA type

We explored the feasibility of obtaining accurate HLA type using pre-existing NGS data not generated for HLA purposes. 83 exomes and 500 targeted NGS pharmacogenomic panels were analyzed using Omixon HLA Explore, OptiType, and/or HLA-Genotyper software. Results were compared against clinical HLA genotyping. 765 (94.2%) Omixon and 769 (94.7%) HLA-Genotyper of 812 germline allele calls across class I/II loci and 402 (99.5%) of 404 OptiType class I calls were concordant to the second field (i.e. HLA-A*02:01). An additional 19 (2.3%) Omixon, 39 (4.8%) HLA-Genotyper, and 2 (0.5%) OptiType allele calls were first field concordant (i.e. HLA-A*02). Using Omixon, four alleles (0.4%) were discordant and 24 (3.0%) failed to call, while 4 alleles (0.4%) were discordant using HLA-Genotyper. Tumor exomes were also evaluated and were 85.4%, 91.6%, and 100% concordant (Omixon and HLA-Genotyper with 96 alleles tested, and Optitype with 48 class I alleles, respectively). The 15 exomes and 500 pharmacogenomic panels were 100% concordant for each pharmacogenomic allele tested. This work has broad implications spanning future clinical care (pharmacogenomics, tumor response to immunotherapy, autoimmunity, etc.) and research applications.

High-Resolution HLA Typing of HLA-A, -B, -C, -DRB1, and -DQB1 in Kinh Vietnamese by Using Next-Generation Sequencing

Human leukocyte antigen (HLA) genotyping displays the particular characteristics of HLA alleles and haplotype frequencies in each population. Although it is considered the current gold standard for HLA typing, high-resolution sequence-based HLA typing is currently unavailable in Kinh Vietnamese populations. In this study, high-resolution sequence-based HLA typing (3-field) was performed using an amplicon-based next-generation sequencing platform to identify the HLA-A, -B, -C, -DRB1, and -DQB1 alleles of 101 unrelated healthy Kinh Vietnamese individuals from southern Vietnam. A total of 28 HLA-A, 41 HLA-B, 21 HLA-C, 26 HLA-DRB1, and 25 HLA-DQB1 alleles were identified. The most frequently occurring HLA alleles were A^∗11:01:01, B^∗15:02:01, C^∗07:02:01, DRB1^∗12:02:01, and DQB1^∗03:01:01. Haplotype calculation showed that A^∗29:01:01∼B^∗07:05:01, DRB1^∗12:02:01∼DQB1^∗3:01:01, A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01, A^∗33:03:01∼B^∗58:01:01∼DRB1^∗03:01:01, and A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01∼DRB1^∗10:01:01∼DQB1^∗05:01:01 were the most common haplotypes in the southern Kinh Vietnamese population. Allele distribution and haplotype analyses demonstrated that the Vietnamese population shares HLA features with South-East Asians but retains unique characteristics. Data from this study will be potentially applicable in medicine and anthropology.

Isolation and characterization of new human carrier peptides from two important vaccine immunogens

In the preparation of glycoconjugate vaccines in clinical practice, two highly immunogenic carrier proteins, CRM₁₉₇ and tetanus toxoid (TT), are predominantly conjugated with the capsular polysaccharides (CPSs) of bacterial pathogens. In addition, TT has long been used as an effective vaccine to prevent tetanus. While these carrier proteins play an important role in immunogenicity and vaccine design alike, their defined human major histocompatibility complex class II (MHCII) T cell epitopes are inadequately characterized. In this current work, we use mass spectrometry to identify the peptides from these carrier proteins that are naturally processed and presented by human B cells via MHCII pathway. The MHCII-presented peptides are screened for their T cell stimulation using primary CD4+ T cells from four healthy adult donors. These combined methods reveal a subset of eleven CD4+ T cell epitopes that proliferate and stimulate human T cells with diverse MHCII allelic repertoire. Six of these peptides stand out as potential immunodominant epitopes by responding in three or more donors. Additionally, we provide evidence of these natural epitopes eliciting more significant T cell responses in donors than previously published TT peptides selected from T cell epitope screening. This study serves toward understanding carrier protein immune responses and thus enables the use of these peptides in developing novel knowledge-based vaccines to combat persisting problems in glycoconjugate vaccine design.

Human Leukocyte Antigen alleles as an aid to STR in complex forensic DNA samples

Human biological samples with multiple contributors remain one of the most challenging aspects of DNA typing within a forensic science context. With the increasing sensitivity of commercially available kits allowing detection of low template DNA, complex mixtures are now a standard component of forensic DNA evidence. Over the years, various methods and techniques have been developed to try to resolve the issue of mixed profiles. However, forensic DNA analysis has relied on the same markers to generate DNA profiles for the past 30 years causing considerable challenges in the deconvolution of complex mixed samples. The future of resolving complicated DNA mixtures may rely on utilising markers that have been previously applied to gene typing of non-forensic relevance. With Massively Parallel Sequencing (MPS), techniques becoming more popular and accessible even epigenetic markers have become a source of interest for forensic scientists. The aim of this review is to consider the potential of alleles from the Human Leukocyte Antigen (HLA) complex as effective forensic markers. While Massively Parallel Sequencing of HLA is routinely used in clinical laboratories in fields such as transplantation, pharmacology or population studies, there have not been any studies testing its suitability for forensic casework samples.

Adapterama II: universal amplicon sequencing on Illumina platforms (TaggiMatrix)

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

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.

The Role of Major Histocompatibility Complex in Organ Transplantation- Donor Specific Anti-Major Histocompatibility Complex Antibodies Analysis Goes to the Next Stage

Organ transplantation has progressed with the comprehension of the major histocompatibility complex (MHC). It is true that the outcome of organ transplantation largely relies on how well rejection is managed. It is no exaggeration to say that to be well acquainted with MHC is a shortcut to control rejection. In human beings, MHC is generally recognized as human leukocyte antigens (HLA). Under the current circumstances, the number of alleles is still increasing, but the function is not completely understood. Their roles in organ transplantation are of vital importance, because mismatches of HLA alleles possibly evoke both cellular and antibody-mediated rejection. Even though the control of cellular rejection has improved by recent advances of immunosuppressants, there is no doubt that antibody-mediated rejection (AMR), which is strongly correlated with donor-specific anti-HLA antibodies (DSA), brings a poor outcome. Thus, to diagnose and treat AMR correctly is a clear proposition. In this review, we would like to focus on the detection of intra-graft DSA as a recent trend. Overall, here we will review the current knowledge regarding MHC, especially with intra-graft DSA, and future perspectives: HLA epitope matching; eplet risk stratification; predicted indirectly recognizable HLA epitopes etc. in the context of organ transplantation.

Development of a rapid and reliable single-tube multiplex real-time PCR method for HLA-A*24:02 genotyping

Aim: HLA-A*24:02 is significantly associated with cutaneous adverse drug reactions caused by aromatic antiepileptic drugs. Here, we aimed to establish a fast and reliable detection method for HLA-A*24:02 genotyping. Methods: A single-tube multiplex quantitative real-time polymerase chain reaction (qPCR) assay for HLA-A*24:02 genotyping was established by combining allele-specific primers with TaqMan probes. Results: A 100% concordance was observed between qPCR and SBT result in 106 Han subjects. The detection limit of the new method was 0.05 ng DNA. The positive rate of HLA-A*24:02 in Tibetans (55.6%, n = 81) was significantly higher than those in Han (34%, n = 106), Uighur (27.5%, n = 102), Bouyei (25.9%, n = 116) and Miao populations (26.5%, n = 113). Conclusion: The newly established qPCR assay was reliable for HLA-A*24:02 screening in clinical applications.

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.

Deciphering the fine nucleotide diversity of full HLA class I and class II genes in a well-documented population from sub-Saharan Africa

With the aim to understand how next‐generation sequencing (NGS) improves both our assessment of genetic variation within populations and our knowledge on HLA molecular evolution, we sequenced and analysed 8 HLA loci in a well‐documented population from sub‐Saharan Africa (Mandenka). The results of full‐gene NGS‐MiSeq sequencing compared with those obtained by traditional typing techniques or limited sequencing strategies showed that segregating sites located outside exon 2 are crucial to describe not only class I but also class II population diversity. A comprehensive analysis of exons 2, 3, 4 and 5 nucleotide diversity at the 8 HLA loci revealed remarkable differences among these gene regions, notably a greater variation concentrated in the antigen recognition sites of class I exons 3 and some class II exons 2, likely associated with their peptide‐presentation function, a lower diversity of HLA‐C exon 3, possibly related to its role as a KIR ligand, and a peculiar molecular diversity of HLA‐A exon 2, revealing demographic signals. Based on full‐length HLA sequences, we also propose that the most frequent DRB1 allele in the studied population, DRB1*13:04, emerged from an allelic conversion involving 3 potential alleles as donors and DRB1*11:02:01 as recipient. Finally, our analysis revealed a high occurrence of the DRB1*13:04‐DQA1*05:05:01‐DQB1*03:19 haplotype, possibly resulting from a selective sweep due to protection to Onchorcerca volvulus, a prevalent pathogen in West Africa. This study unveils highly relevant information on the molecular evolution of HLA genes in relation to their immune function, calling for similar analyses in other populations living in contrasting environments.

High-Resolution Full-Length HLA Typing Method Using Third Generation (Pac-Bio SMRT) Sequencing Technology

The human HLA genes are among the most polymorphic genes in the human genome. Therefore, it is very difficult to find two unrelated individuals with identical HLA molecules. As a result, HLA Class I and Class II genes are routinely sequenced or serotyped for organ transplantation, autoimmune disease-association studies, drug hypersensitivity research, and other applications. However, these methods were able to give two or four digit data, which was not sufficient enough to understand the completeness of haplotypes of HLA genes. To overcome these limitations, we here described end-to-end workflow for sequencing of HLA class I and class II genes using third generation sequencing, SMRT technology. This method produces fully-phased, unambiguous, allele-level information on the PacBio System.

Clinical features, biochemistry and HLA-DRB1 status in youth-onset type 1 diabetes in Pakistan

Published information on diabetes in Pakistani youth is limited. We aimed to investigate the demographic, clinical, and biochemical features, and HLA-DRB1 alleles in new cases of diabetes affecting children and adolescents <22 years of age. The study was conducted at Baqai Institute of Diabetology and Endocrinology in Karachi from June 2013-December 2015. One hundred subjects aged <22 years at diagnosis were enrolled. Demographic characteristics, clinical information, biochemical parameters (blood glucose, HbA1c, C-peptide, glutamic acid decarboxylase 65 (GAD65) and islet antigen 2 (IA-2) autoantibodies) were measured. DNA from 100 subjects and 200 controls was extracted and genotyped for HLA-DRB1 using high-resolution genotyping technology. Ninety-nine subjects were clinically diagnosed as type 1 diabetes (T1D) and one as type 2 diabetes (T2D). Of the 99 with T1D, 57 (57.6%) were males and 42 (42.4%) females, with mean age at diagnosis 11.0 ± 5.2 years (range 1.6-21.7 years) and peaks at six and fifteen years. Fifty-seven subjects were assessed within one month of diagnosis and all within eleven months. For the subjects diagnosed as T1D, mean C-peptide was 0.63 ± 0.51 nmol/L (1.91 ± 1.53 ng/mL), with 16 (16.2%) IA2 positive, 53 (53.5%) GAD-65 positive, and 10 (10.1%) positive for both autoantibodies. In T1D patients, the allele DRB1*03:01 demonstrated highly significant T1D association (p < 10^-16), with no apparent risk conferred by DRB1*04:xx alleles. CONCLUSIONS: Heterogeneous forms of T1D appear more common in children and youth in Pakistan than in European populations. Individual understanding of such cases could enable improved management strategies and healthier outcomes.

"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.

HLA typing using next-generation sequencing for Chinese juvenile- and adult-onset myasthenia gravis patients

To compare HLA typing between juvenile- and adult-onset myasthenia gravis (MG), we enrolled 101 children (age ≤12 years) and 168 adults (age ≥20 years) with MG. We excluded patients with histories of thymoma, thyroid disease, or other autoimmune disease. We selected 41 seronegative juvenile-onset patients with ocular symptoms only, and 41 seropositive adult-onset patients with generalized symptoms. We used next-generation sequencing for typing and analysis of HLA genes (Loci: A, B, C, DPA1, DPB1, DQA1, DQB1 and DRB1). Haplotypes HLA-A∗02:07:01-B∗46:01:01-C∗01:02:01-DQA1∗01:01:01-DQB1∗03:03:02-DRB1∗09:01:02, HLA-A∗11:01:01, HLA-A∗24:02:01, and HLA-DPA1∗02:02:02 were found to be related to juvenile-onset MG and HLA-A∗01:01:01, HLA-A∗02:03:01, HLA-C∗03:04:01, and HLA-DQB1∗06:02:01 to adult-onset MG. Therefore, our findings suggested that HLA typing might determine the heterogeneity between AChR-Ab negative juvenile-onset and AChR-Ab positive adult-onset Chinese MG patients.

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.

Epidemiology of childhood-onset type 1 diabetes in Azerbaijan: Incidence, clinical features, biochemistry, and HLA-DRB1 status

AIMS

Determine the incidence and typology of diabetes in children in Azerbaijan.

METHODS

Clinical features, C-peptide, autoantibodies (glutamic acid decarboxylase 65 (GAD65) and islet antigen 2 (IA-2)), and HLA-DRB1 status were studied in 106 subjects <18 years of age who were recently diagnosed. 104 cases were consecutive. Incidence was determined for Baku and Absheron regions, where ascertainment is estimated to be essentially 100%.

RESULTS

104 of the 106 (98%) were diagnosed with type 1 diabetes, one with type 2 diabetes and one with atypical diabetes. Type 1 diabetes incidence in Baku City and Absheron was 7.05 per 100,000 population <15 years per year. Peak age of onset was 10 years. There was a slight male preponderance (male:female 1.17:1), and no temporal association with seasons. Almost all type 1 diabetes subjects presented with classic symptoms including a high incidence (58%) of diabetic ketoacidosis. 86% presented with low C-peptide values (<0.13 nmol/L, <0.40 ng/mL) and 74% were positive for at least one type 1 diabetes-related autoantibody.

CONCLUSIONS

Azerbaijan has a moderate type 1 diabetes incidence and clinical, biochemical and genetic features similar to that in European populations.

Full-length next-generation sequencing of HLA class I and II genes in a cohort from Thailand

The human leukocyte antigen (HLA) genes are highly variable and are known to play an important role in disease outcomes, including infectious diseases. Prior knowledge of HLA polymorphisms in a population usually forms the basis for an effective case-control study design. As a prelude to future disease association analyses, we report HLA class I and II diversity in 334 unrelated donors from a Dengue vaccine efficacy trial conducted in Thailand. Long-range PCR amplification of six HLA loci was performed on DNA extracted from saliva samples. HLA-A, -B, -C, -DPB1, -DQB1 and -DRB1 were genotyped using a next-generation sequencing method presented at the 17th International HLA and Immunogenetics Workshop. In total, we identified 201 HLA alleles, including 35 HLA-A, 57 HLA-B, 28 HLA-C, 24 HLA-DPB1, 21 HLA-DQB1 and 36 HLA-DRB1 alleles. Very common HLA alleles with frequencies greater than 10 percent were A∗11:01:01, A∗33:03:01, A∗24:02:01, B∗46:01:01, C∗07:02:01, C∗01:02:01, C∗08:01:01, DPB1∗05:01:01, DPB1∗13:01:01, DPB1∗04:01:01, DPB1∗02:01:02, DQB1∗03:01:01, DQB1∗05:02:01, DQB1∗03:03:02, DRB1∗12:02:01, DRB1∗09:01:02, and DRB1∗15:02:01. A novel HLA allele, B∗15:450, had a non-synonymous substitution and occurred in more than one donor. Population-based full-length NGS HLA typing is more conclusive and provides a sound foundation for exploring disease association in a given population.

New DQA1 allele specific antibody against epitope 2D (an exon 1 encoded amino acid). Considerations for alleles under the same P-group designation

The majority of polymorphisms of the Human Leukocyte Antigen (HLA) proteins are clustered at the peptide binding domain (PBD), defined as the area coded by exon 2 and 3 for class I and exon 2 for class II. HLA alleles with the same amino acid (AA) sequence at the PBD are considered functionally equivalent and can be grouped under the same P-group designation. Here we present a case of a kidney recipient, typed as DQA1*01:04, 01:05 and DQB1*05:01P, 05:03, who developed antibodies against all DQ antigens on our Luminex Single Antigen (LSA) panel. Our LSA panel does not include DQA1*01:05 or 01:04, but both alleles belong to the DQA1*01:01P group and beads carrying DQA1*01:01 tested positive. Mature protein sequence alignment demonstrated a single AA mismatch between DQA1*01:04/01:05 and DQA1*01:01 located at position 2 (G vs D), which is encoded by exon 1. Luminex assay by another manufacturer which include a bead carrying patient's own DQA1* type and crossmatch studies with surrogate donors confirmed the presence of an antibody against mismatched epitope 2D. This case illustrates that alleles included in the same P-group may have polymorphisms able to trigger immunological responses and brings attention to the fact that some mature HLA proteins express AA encoded by exon 1, which is structurally part of the PBD.

Multiplexed Nanopore Sequencing of HLA-B Locus in Māori and Pacific Island Samples

The human leukocyte antigen (HLA) system encodes the human major histocompatibility complex (MHC). HLA-B is the most polymorphic gene in the MHC class I region and many HLA-B alleles have been associated with adverse drug reactions (ADRs) and disease susceptibility. The frequency of such HLA-B alleles varies by ethnicity, and therefore it is important to understand the prevalence of such alleles in different population groups. Research into HLA involvement in ADRs would be facilitated by improved methods for genotyping key HLA-B alleles. Here, we describe an approach to HLA-B typing using next generation sequencing (NGS) on the MinION™ nanopore sequencer, combined with data analysis with the SeqNext-HLA software package. The nanopore sequencer offers the advantages of long-read capability and single molecule reads, which can facilitate effective haplotyping. We developed this method using reference samples as well as individuals of New Zealand Māori or Pacific Island descent, because HLA-B diversity in these populations is not well understood. We demonstrate here that nanopore sequencing of barcoded, pooled, 943 bp polymerase chain reaction (PCR) amplicons of 49 DNA samples generated ample read depth for all samples. HLA-B alleles were assigned to all samples at high-resolution with very little ambiguity. Our method is a scaleable and efficient approach for genotyping HLA-B and potentially any other HLA locus. Finally, we report our findings on HLA-B genotypes of this cohort, which adds to our understanding of HLA-B allele frequencies among Māori and Pacific Island people.

Construction of full-length Japanese reference panel of class I HLA genes with single-molecule, real-time sequencing

Human leukocyte antigen (HLA) is a gene complex known for its exceptional diversity across populations, importance in organ and blood stem cell transplantation, and associations of specific alleles with various diseases. We constructed a Japanese reference panel of class I HLA genes (ToMMo HLA panel), comprising a distinct set of HLA-A, HLA-B, HLA-C, and HLA-H alleles, by single-molecule, real-time (SMRT) sequencing of 208 individuals included in the 1070 whole-genome Japanese reference panel (1KJPN). For high-quality allele reconstruction, we developed a novel pipeline, Primer-Separation Assembly and Refinement Pipeline (PSARP), in which the SMRT sequencing and additional short-read data were used. The panel consisted of 139 alleles, which were all extended from known IPD-IMGT/HLA sequences, contained 40 with novel variants, and captured more than 96.5% of allelic diversity in 1KJPN. These newly available sequences would be important resources for research and clinical applications including high-resolution HLA typing, genetic association studies, and analyzes of cis-regulatory elements.

High resolution HLA analysis reveals independent class I haplotypes and amino-acid motifs protective for multiple sclerosis

We investigated association between HLA class I and class II alleles and haplotypes, and KIR loci and their HLA class I ligands, with multiple sclerosis (MS) in 412 European-American MS patients and 419 ethnically-matched controls, using next generation sequencing. The DRB1*15:01~DQB1*06:02 haplotype was highly predisposing (odds ratio (OR) = 3.98; 95% confidence interval (CI) = 3−5.31; p-value (p) = 2.22E−16), as was DRB1*03:01~DQB1*02:01 (OR = 1.63; CI = 1.19–2.24; p = 1.41E−03). Hardy-Weinberg (HW) analysis in MS patients revealed a significant DRB1*03:01~DQB1*02:01 homozyote excess (15 observed, 8.6 expected; p = 0.016). The OR for this genotype (5.27; CI = 1.47–28.52; p = 0.0036) suggests a recessive MS risk model. Controls displayed no HW deviations. The C*03:04~B*40:01 haplotype (OR = 0.27; CI = 0.14–0.51; p = 6.76E−06) was highly protective for MS, especially in haplotypes with A*02:01 (OR = 0.15; CI = 0.04–0.45; p = 6.51E−05). By itself, A*02:01 is moderately protective, (OR = 0.69; CI = 0.54–0.87; p = 1.46E−03), and haplotypes of A*02:01 with the HLA-B Thr80 Bw4 variant (Bw4T) more so (OR = 0.53; CI = 0.35–0.78; p = 7.55E−04). Protective associations with the Bw4 KIR ligand resulted from linkage disequilibrium (LD) with DRB1*15:01, but the Bw4T variant was protective (OR = 0.64; CI = 0.49–0.82; p = 3.37E−04) independent of LD with DRB1*15:01. The Bw4I variant was not associated with MS. Overall, we find specific class I HLA polymorphisms to be protective for MS, independent of the strong predisposition conferred by DRB1*15:01.

Application of High-Throughput Next-Generation Sequencing for HLA Typing on Buccal Extracted DNA

Next-generation sequencing (NGS) is increasingly recognized for its ability to deliver high-resolution and high-throughput HLA genotyping. As a result, there is active interest in applying NGS technologies to perform high volume bone marrow donor recruitment typing. Currently, buccal-based DNA specimens are considered a noninvasive and cost-effective method for registry typing. Here, we describe the feasibility of using long-range PCR and clonal sequencing by Illumina MiSeq to deliver unambiguous HLA typing on buccal-based donor recruitment samples.

Hapl-o-Mat: open-source software for HLA haplotype frequency estimation from ambiguous and heterogeneous data

BACKGROUND

Knowledge of HLA haplotypes is helpful in many settings as disease association studies, population genetics, or hematopoietic stem cell transplantation. Regarding the recruitment of unrelated hematopoietic stem cell donors, HLA haplotype frequencies of specific populations are used to optimize both donor searches for individual patients and strategic donor registry planning. However, the estimation of haplotype frequencies from HLA genotyping data is challenged by the large amount of genotype data, the complex HLA nomenclature, and the heterogeneous and ambiguous nature of typing records.

RESULTS

To meet these challenges, we have developed the open-source software Hapl-o-Mat. It estimates haplotype frequencies from population data including an arbitrary number of loci using an expectation-maximization algorithm. Its key features are the processing of different HLA typing resolutions within a given population sample and the handling of ambiguities recorded via multiple allele codes or genotype list strings. Implemented in C++, Hapl-o-Mat facilitates efficient haplotype frequency estimation from large amounts of genotype data. We demonstrate its accuracy and performance on the basis of artificial and real genotype data.

CONCLUSIONS

Hapl-o-Mat is a versatile and efficient software for HLA haplotype frequency estimation. Its capability of processing various forms of HLA genotype data allows for a straightforward haplotype frequency estimation from typing records usually found in stem cell donor registries.

2.7 million samples genotyped for HLA by next generation sequencing: lessons learned

Background

At the DKMS Life Science Lab, Next Generation Sequencing (NGS) has been used for ultra-high-volume high-resolution genotyping of HLA loci for the last three and a half years. Here, we report on our experiences in genotyping the HLA, CCR5, ABO, RHD and KIR genes using a direct amplicon sequencing approach on Illumina MiSeq and HiSeq 2500 instruments.

Results

Between January 2013 and June 2016, 2,714,110 samples largely from German, Polish and UK-based potential stem cell donors have been processed. 98.9% of all alleles for the targeted HLA loci (HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1) were typed at high resolution or better. Initially a simple three-step workflow based on nanofluidic chips in conjunction with 4-primer amplicon tagging was used. Over time, we found that this setup results in PCR artefacts such as primer dimers and PCR-mediated recombination, which may necessitate repeat typing. Split workflows for low- and high-DNA-concentration samples helped alleviate these problems and reduced average per-locus repeat rates from 3.1 to 1.3%. Further optimisations of the workflow included the use of phosphorothioate oligos to reduce primer degradation and primer dimer formation, and employing statistical models to predict read yield from initial template DNA concentration to avoid intermediate quantification of PCR products. Finally, despite the populations typed at DKMS Life Science Lab being relatively homogenous genetically, an analysis of 1.4 million donors processed between January 2015 and May 2016 led to the discovery of 1,919 distinct novel HLA alleles.

Conclusions

Amplicon-based NGS HLA genotyping workflows have become the workhorse in high-volume tissue typing of registry donors. The optimisation of workflow practices over multiple years has led to insights and solutions that improve the efficiency and robustness of short amplicon based genotyping workflows.

Electronic supplementary material

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

Comparison of exome-based HLA class I genotyping tools: identification of platform-specific genotyping errors

Accurate human leukocyte antigen (HLA) genotyping is critical in studies involving the immune system. Several algorithms to estimate HLA genotypes from whole-exome data were developed. We compared the accuracy of seven algorithms, including Optitype, Polysolver and PHLAT, as well as investigated patterns and possible causes of miscalls using 12 clinical samples and 961 individuals from the 1000 Genomes Project. Optitype showed the highest accuracy of 97.2% for HLA class I alleles at the second field resolution, followed by 94.0% in Polysolver and 85.6% in PHLAT. In Optitype, 34 (21.1%) of 161 miscalls were across different serological types, and common miscalls were HLA-A*26:01 to HLA-A*25:01, HLA-B*45:01 to HLA-B*44:15 and HLA-C*08:02 to HLA-C*05:01 with error rates of 4.1%, 10.0% and 4.1%, respectively. In Polysolver, 193 (55.9%) of 345 miscalls occurred across different serological alleles, and a specific pattern of genotyping error from HLA-A*25:01 to HLA-A*26:01 was observed in 93.3% of HLA-A*25:01 carriers, due to dropping of HLA-A*25:01 sequence reads during the extraction process of HLA reads. In PHLAT, 147 (59.8%) of 246 miscalls in HLA-A were due to erroneous assignment of multiple alleles to either HLA-A*01:22 or HLA-A*01:81. These results suggest that careful considerations needed to be taken when using exome-based HLA class I genotyping data and applying these results in clinical settings.