HLA genotyping in the clinical laboratory: comparison of next-generation sequencing methods

HLA-DQ genotype distribution and risk evaluation of celiac disease in Northwest China

BACKGROUND

Celiac disease (CD) is an autoimmune small bowel disease. Genetic susceptibility for CD is mainly determined by the human leukocyte antigen (HLA)-DQ haplotypes. The risk of CD conferred by HLA genotypes varies geographically and across populations, however, this has not yet been documented in Chinese patients with CD.

AIMS

To investigate the distribution of HLA-DQ and the related risks of CD development in Northwest China.

METHODS

A total of 75 CD patients and 300 healthy individuals were genotyped for HLA-DQ using the Illumina NextSeq, and the relative risks of the different genotypes were evaluated.

RESULTS

In total, 68.00% of CD patients and 21.00% of controls carried HLA-DQ2.5 heterodimers (p < 0.001). We identified four CD risk gradients. Individuals carrying a double dose of DQB1*02 had the highest risk of developing CD (1:16); however, with heterozygosis (DQB1*02:02/DQB1*02:01) having the highest risk (1:9). HLA-DQ2.5 individuals with a single copy of HLA-DQB1*02, in either the cis or trans configuration, were at a medium risk (1:38). Non-DQ2.5 carriers of DQ8 or DQ2.2 were at low risk, while only carriers of DQ7.5 or DQX.5 were at very low risk. Patients with the HLA-DQ2.5 genotype had more severe mucosal damage compared with the HLA-DQ2.5 genotype negative CD patients (70.59% vs. 41.67%, p = 0.016).

CONCLUSION

Genetic susceptibility to CD is highly prevalent in the Northwest Chinese population and the highest risk of developing CD was associated with the DQ2.5/DQ2.2 genotype. The DQ2.5 allele is involved in the severity of mucosal injury.

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.

Application of NGS Technology in Understanding the Pathology of Autoimmune Diseases

NGS technologies have transformed clinical diagnostics and broadly used from neonatal emergencies to adult conditions where the diagnosis cannot be made based on clinical symptoms. Autoimmune diseases reveal complicate molecular background and traditional methods could not fully capture them. Certainly, NGS technologies meet the needs of modern exploratory research, diagnostic and pharmacotherapy. Therefore, the main purpose of this review was to briefly present the application of NGS technology used in recent years in the understanding of autoimmune diseases paying particular attention to autoimmune connective tissue diseases. The main issues are presented in four parts: (a) panels, whole-genome and -exome sequencing (WGS and WES) in diagnostic, (b) Human leukocyte antigens (HLA) as a diagnostic tool, (c) RNAseq, (d) microRNA and (f) microbiome. Although all these areas of research are extensive, it seems that epigenetic impact on the development of systemic autoimmune diseases will set trends for future studies on this area.

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.

A new set of reagents and related software used for NGS based classical and non-classical HLA typing showing evidence for a greater HLA haplotype diversity

To evaluate the HLA typing performance of a new Long-Range PCR NGS set of reagents and its dedicated software, a panel of 41 reference homozygous cell lines from the International Histocompatibility Working Group (IHWG) and a panel of 376 volunteer bone marrow donors were analyzed for classical and non-classical HLA class I and class II genes. All results, except HLA-DPB1, were obtained without any ambiguities at the 3rd field level. Based on the high resolution performance of the reagents, a number of new alleles have been described not only for classical but also for non-classical HLA class I genes, leading to a more accurate haplotype definition. Linkage disequilibrium between HLA-A and HLA-G genes has been defined at 4th field level of resolution. Moreover, for the first time, HLA-DQA2 and DQB2 polymorphisms and their linkage disequilibrium with DQB1 were described.

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.

Long-read sequencing in deciphering human genetics to a greater depth

Through four decades' development, DNA sequencing has inched into the era of single-molecule sequencing (SMS), or the third-generation sequencing (TGS), as represented by two distinct technical approaches developed independently by Pacific Bioscience (PacBio) and Oxford Nanopore Technologies (ONT). Historically, each generation of sequencing technologies was marked by innovative technological achievements and novel applications. Long reads (LRs) are considered as the most advantageous feature of SMS shared by both PacBio and ONT to distinguish SMS from next-generation sequencing (NGS, or the second-generation sequencing) and Sanger sequencing (the first-generation sequencing). Long reads overcome the limitations of NGS and drastically improves the quality of genome assembly. Besides, ONT also contributes several unique features including ultra-long reads (ULRs) with read length above 300 kb and some close to 1 million bp, direct RNA sequencing and superior portability as made possible by pocket-sized MinION sequencer. Here, we review the history of DNA sequencing technologies and associated applications, with a special focus on the advantages as well as the limitations of ULR sequencing in genome assembly.

Complete RHD next-generation sequencing: establishment of reference RHD alleles

The Rh blood group system (ISBT004) is the second most important blood group after ABO and is the most polymorphic one, with 55 antigens encoded by 2 genes, RHD and RHCE This research uses next-generation sequencing (NGS) to sequence the complete RHD gene by amplifying the whole gene using overlapping long-range polymerase chain reaction (LR-PCR) amplicons. The aim was to study different RHD alleles present in the population to establish reference RHD allele sequences by using the analysis of intronic single-nucleotide polymorphisms (SNPs) and their correlation to a specific Rh haplotype. Genomic DNA samples (n = 69) from blood donors of different serologically predicted genotypes including R₁R₁ (DCe/DCe), R₂R₂ (DcE/DcE), R₁R₂ (DCe/DcE), R₂R_Z (DcE/DCE), R₁r (DCe/dce), R₂r (DcE/dce), and R₀r (Dce/dce) were sequenced and data were then mapped to the human genome reference sequence hg38. We focused on the analysis of hemizygous samples, as these by definition will only have a single copy of RHD For the 69 samples sequenced, different exonic SNPs were detected that correlate with known variants. Multiple intronic SNPs were found in all samples: 21 intronic SNPs were present in all samples indicating their specificity to the RHD*DAU0 (RHD*10.00) haplotype which the hg38 reference sequence encodes. Twenty-three intronic SNPs were found to be R₂ haplotype specific, and 15 were linked to R₁, R₀, and R_Z haplotypes. In conclusion, intronic SNPs may represent a novel diagnostic approach to investigate known and novel variants of the RHD and RHCE genes, while being a useful approach to establish reference RHD allele sequences.

Suitability of dried DNA for long-range PCR amplification and HLA typing by next-generation sequencing

Storage and stable shipment of genomic DNA are of great concern to laboratories that may need to perform testing off archived samples. There are some dry-state storage methods that are available that have the potential to provide a way to store samples at room temperature for long periods of time as well as offer a means to ship DNA to other facilities without the same safety concerns that come with shipping liquid samples. The recovered DNA should be of sufficient integrity such that downstream applications can be performed without concern of the sample quality. This work describes sample properties between two methods of DNA storage, dried (room temperature) and traditional (-80 °C). DNA was evaluated for purity, fragment length, and the ability to generate HLA typing using next-generation sequencing.

High-Throughput Contiguous Full-Length Next-Generation Sequencing of HLA Class I and II Genes from 96 Donors in a Single MiSeq Run

The human leukocyte antigen (HLA) genes regulate and drive the immune system, and are among the most polymorphic loci in the human genome. HLA diversity is known to play an important role in transplantation and disease association studies. There are multiple approaches to DNA-based HLA genotyping and recent advances in next-generation sequencing (NGS) technologies have facilitated the development of whole gene sequencing methods. We describe an accurate, efficient, scalable, and cost-effective approach to contiguously amplify and sequence full-length genes of six HLA class I and II loci from 96 individuals on a single Illumina MiSeq run.

Different DRB1*03:01-DQB1*02:01 haplotypes confer different risk for celiac disease

Celiac disease is associated with the HLA-DR3-DQA1*05:01-DQB1*02:01 and DR4-DQA1*03:01-DQB1*03:02 haplotypes. In addition, there are currently over 40 non-HLA loci associated with celiac disease. This study extends previous analyses on different HLA haplotypes in celiac disease using next generation targeted sequencing. Included were 143 patients with celiac disease and 135 non-celiac disease controls investigated at median 9.8 years (1.4-18.3 years). PCR-based amplification of HLA and sequencing with Illumina MiSeq technology were used for extended sequencing of the HLA class II haplotypes HLA-DRB1, DRB3, DRB4, DRB5, DQA1 and DQB1, respectively. Odds ratios were computed marginally for every allele and haplotype as the ratio of allelic frequency in patients and controls as ratio of exposure rates (RR), when comparing a null reference with equal exposure rates in cases and controls. Among the extended HLA haplotypes, the strongest risk haplotype for celiac disease was shown for DRB3*01:01:02 in linkage with DQA1*05:01-DQB1*02:01 (RR = 6.34; P-value < .0001). In a subpopulation analysis, DRB3*01:01:02-DQA1*05:01-DQB1*02:01 remained the most significant in patients with Scandinavian ethnicity (RR = 4.63; P < .0001) whereas DRB1*07:01:01-DRB4*01:03:01-DQA1*02:01-DQB1*02:02:01 presented the highest risk of celiac disease among non-Scandinavians (RR = 7.94; P = .011). The data also revealed 2 distinct celiac disease risk DR3-DQA1*05:01-DQB*02:01 haplotypes distinguished by either the DRB3*01:01:02 or DRB3*02:02:01 alleles, indicating that different DRB1*03:01-DQB1*02:01 haplotypes confer different risk for celiac disease. The associated risk of celiac disease for DR3-DRB3*01:01:02-DQA1*05:01-DQB1*02:01 is predominant among patients of Scandinavian ethnicity.

Transplant genetics and genomics

Ever since the discovery of the major histocompatibility complex, scientific and clinical understanding in the field of transplantation has been advanced through genetic and genomic studies. Candidate-gene approaches and recent genome-wide association studies (GWAS) have enabled a deeper understanding of the complex interplay of the donor-recipient interactions that lead to transplant tolerance or rejection. Genetic analysis in transplantation, when linked to demographic and clinical outcomes, has the potential to drive personalized medicine by enabling individualized risk stratification and immunosuppression through the identification of variants associated with immune-mediated complications, post-transplant disease or alterations in drug-metabolizing genes.

Targeted Next-Generation Sequencing for Human Leukocyte Antigen Typing in a Clinical Laboratory: Metrics of Relevance and Considerations for Its Successful Implementation

CONTEXT

- Numerous feasibility studies to type human leukocyte antigens (HLAs) by next-generation sequencing (NGS) have led to the development of vendor-supported kits for HLA typing by NGS. Some clinical laboratories have introduced HLA-NGS, and many are investigating the introduction. Standards from accrediting agencies form the regulatory framework for introducing this test into clinical laboratories.

OBJECTIVES

- To provide an assessment of metrics and considerations relevant to the successful implementation of clinical HLA-NGS typing, and to provide as a reference a validated HLA-NGS protocol used clinically since December 2013 at the Children's Hospital of Philadelphia (Philadelphia, Pennsylvania).

DATA SOURCES

- The HLA-NGS has been performed on 2532 samples. The initial 1046 and all homozygous samples were also typed by an alternate method. The HLA-NGS demonstrated 99.7% concordance with the alternate method. Ambiguous results were most common at the DPB1 locus because of a lack of phasing between exons 2 and 3 or the unsequenced exon 1 (533 of 2954 alleles; 18.04%) and the DRB1 locus because of not sequencing exon 1 (75 of 3972 alleles; 1.89%). No ambiguities were detected among the other loci. Except for 2 false homozygous samples, all homozygous samples (1891) demonstrated concordance with the alternate method. The article is organized to address the critical elements in the preanalytic, analytic, and postanalytic phases of introducing this assay into the clinical laboratory.

CONCLUSIONS

- The results demonstrate that HLA typing by NGS is a highly accurate, reproducible, efficient method that provides more-complete sequencing information for the length of the HLA gene and can be the single methodology for HLA typing in clinical immunogenetics laboratories.

High resolution HLA haplotyping by imputation for a British population bioresource

This study aimed to establish the occurrence and frequency of HLA alleles and haplotypes for a healthy British Caucasian population bioresource from Oxfordshire. We present the results of imputation from HLA SNP genotyping data using SNP2HLA for 5553 individuals from Oxford Biobank, defining one- and two-field alleles together with amino acid polymorphisms. We show that this achieves a high level of accuracy with validation using sequence-specific primer amplification PCR. We define six- and eight-locus HLA haplotypes for this population by Bayesian methods implemented using PHASE. We determine patterns of linkage disequilibrium and recombination for these individuals involving classical HLA loci and show how analysis within a haplotype block structure may be more tractable for imputed data. Our findings contribute to knowledge of HLA diversity in healthy populations and further validate future large-scale use of HLA imputation as an informative approach in population bioresources.