The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection.

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

A combination of HLA-DP α and ß chain polymorphisms paired with a SNP in the DPB1 3' UTR region, denoting expression levels, are associated with atopic dermatitis.

Introduction: Components of the immune response have previously been associated with the pathophysiology of atopic dermatitis (AD), specifically the Human Leukocyte Antigen (HLA) Class II region via genome-wide association studies, however the exact elements have not been identified. Methods: This study examines the genetic variation of HLA Class II genes using next generation sequencing (NGS) and evaluates the resultant amino acids, with particular attention on binding site residues, for associations with AD. The Genetics of AD cohort was used to evaluate HLA Class II allelic variation on 464 subjects with AD and 384 controls. Results: Statistically significant associations with HLA-DP α and ß alleles and specific amino acids were found, some conferring susceptibility to AD and others with a protective effect. Evaluation of polymorphic residues in DP binding pockets revealed the critical role of P1 and P6 (P1: α31M + (ß84G or ß84V) [protection]; α31Q + ß84D [susceptibility] and P6: α11A + ß11G [protection]) and were replicated with a national cohort of children consisting of 424 AD subjects. Independently, AD susceptibility-associated residues were associated with the G polymorphism of SNP rs9277534 in the 3' UTR of the HLA-DPB1 gene, denoting higher expression of these HLA-DP alleles, while protection-associated residues were associated with the A polymorphism, denoting lower expression. Discussion: These findings lay the foundation for evaluating non-self-antigens suspected to be associated with AD as they potentially interact with particular HLA Class II subcomponents, forming a complex involved in the pathophysiology of AD. It is possible that a combination of structural HLA-DP components and levels of expression of these components contribute to AD pathophysiology.

Genomic characterization of HLA class I and class II genes in ethnically diverse sub-Saharan African populations: A report on novel HLA alleles.

HLA allelic variation has been well studied and documented in many parts of the world. However, African populations have been relatively under-represented in studies of HLA variation. We have characterized HLA variation from 489 individuals belonging to 13 ethnically diverse populations from rural communities from the African countries of Botswana, Cameroon, Ethiopia, and Tanzania, known to practice traditional subsistence lifestyles using next generation sequencing (Illumina) and long-reads from Oxford Nanopore Technologies. We identified 342 distinct alleles among the 11 HLA targeted genes: HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQA1, -DQB1, -DPA1, and -DPB1, with 140 of those alleles containing novel sequences that were submitted to the IPD-IMGT/HLA database. Sixteen of the 140 alleles contained novel content within the exonic regions of the genes, while 110 alleles contained novel intronic variants. Four alleles were found to be recombinants of already described HLA alleles and 10 alleles extended the sequence content of already described alleles. All 140 alleles include complete allelic sequence from the 5' UTR to the 3' UTR that are inclusive of all exons and introns. This report characterizes the HLA allelic variation from these individuals and describes the novel allelic variation present within these specific African populations.

Identification of HLA-DPA1*01:03:01:57 and HLA-DPA1*02:01:01:29 from a case-control study of atopic dermatitis.

DPA1*01:03:01:57 and DPA1*02:01:01:29 differ by a single nucleotide from their closest references, DPA1*01:03:01:02 and DPA1*02:01:01:06.

Pathogenicity and impact of HLA class I alleles in aplastic anemia patients of different ethnicities.

Acquired aplastic anemia (AA) is caused by autoreactive T cell-mediated destruction of early hematopoietic cells. Somatic loss of human leukocyte antigen (HLA) class I alleles was identified as a mechanism of immune escape in surviving hematopoietic cells of some patients with AA. However, pathogenicity, structural characteristics, and clinical impact of specific HLA alleles in AA remain poorly understood. Here, we evaluated somatic HLA loss in 505 patients with AA from 2 multi-institutional cohorts. Using a combination of HLA mutation frequencies, peptide-binding structures, and association with AA in an independent cohort of 6,323 patients from the National Marrow Donor Program, we identified 19 AA risk alleles and 12 non-risk alleles and established a potentially novel AA HLA pathogenicity stratification. Our results define pathogenicity for the majority of common HLA-A/B alleles across diverse populations. Our study demonstrates that HLA alleles confer different risks of developing AA, but once AA develops, specific alleles are not associated with response to immunosuppression or transplant outcomes. However, higher pathogenicity alleles, particularly HLA-B*14:02, are associated with higher rates of clonal evolution in adult patients with AA. Our study provides insights into the immune pathogenesis of AA, opening the door to future autoantigen identification and improved understanding of clonal evolution in AA.

KIR-favorable TCR-αß/CD19-depleted haploidentical HCT in children with ALL/AML/MDS: primary analysis of the PTCTC ONC1401 trial.

We performed a prospective multicenter study of T-cell receptor αß (TCR-αß)/CD19-depleted haploidentical hematopoietic cell transplantation (HCT) in children with acute leukemia and myelodysplastic syndrome (MDS), to determine 1-year disease-free survival (DFS) and compare 2-year outcomes with recipients of other donor cell sources. Fifty-one patients aged 0.7 to 21 years were enrolled; donors were killer immunoglobulin-like receptor (KIR) favorable based on ligand mismatch and/or high B content. The 1-year DFS was 78%. Superior 2-year DFS and overall survival (OS) were noted in patients <10 years of age, those treated with reduced toxicity conditioning (RTC) rather than myeloablative conditioning, and children with minimal residual disease <0.01% before HCT. Multivariate analysis comparing the KIR-favorable haploidentical cohort with controls showed similar DFS and OS compared with other donor cell sources. Multivariate analysis also showed a marked decrease in the risk of grades 2 to 4 and 3 to 4 acute graft versus host disease (aGVHD), chronic GVHD, and transplant-related mortality vs other donor cell sources. Ethnic and racial minorities accounted for 53% of enrolled patients, and data from a large cohort of recipients/donors screened for KIR showed that >80% of recipients had a KIR-favorable donor by our definition, demonstrating that this approach is broadly applicable to groups often unable to find donors. This prospective, multicenter study showed improved outcomes using TCR-αß/CD19-depleted haploidentical donors using RTC for children with acute leukemia and MDS. Randomized trials comparing this approach with matched unrelated donors are warranted. This trial was registered at https://clinicaltrials.gov as #NCT02646839.

HLA Class I Polymorphisms Influencing Both Peptide Binding and KIR Interactions Are Associated with Remission among Children with Atopic Dermatitis: A Longitudinal Study.

Atopic dermatitis (AD) is a disease of immune dysregulation and skin barrier dysfunction with a relapsing, remitting course and has been associated with several different genetic risk variants. HLA represent a highly variable set of genes that code for cell surface protein molecules involved in the Ag-specific immune response, including the regulation or functioning of T cells, NK cells, and APCs. The purpose of this study was to evaluate associations between HLA class I polymorphisms and the progression of AD over time. We evaluated the associations of AD symptoms and HLA class I polymorphisms based on high-resolution two-field typing in a longitudinal cohort of children with AD (up to 10 y of follow-up). Seven hundred and ninety-two children were evaluated every 6 mo, resulting in 12,752 AD evaluations. Using generalized estimating equations and corrected p values, B*44:02 was found to be associated with AD remission (1.83 [1.35, 2.47]; p = 0.0015). The HLA-B residues at position 116 (d-aspartate) and 80 (T-threonine) were associated with remission (1.42 [1.13, 1.76], p = 0.003; corrected p = 0.028) and (1.45 [1.17, 1.80], p = 0.0008; corrected p = 0.0024), respectively. B80T is a killer-cell Ig-like receptor (KIR) site. Our findings reveal that two axes of immune response (T cell and NK cell) may influence disease progression. Identifying binding pocket changes in addition to other factors (e.g., allergens) that increase the risk or severity of AD can improve our understanding of the immunologic mechanisms associated with AD and may lead to personalized therapies for improving patient care.

Spatial constrains and information content of sub-genomic regions of the human genome.

Complexity metrics and machine learning (ML) models have been utilized to analyze the lengths of segmental genomic entities of DNA sequences (exonic, intronic, intergenic, repeat, unique) with the purpose to ask questions regarding the segmental organization of the human genome within the size distribution of these sequences. For this we developed an integrated methodology that is based upon the reconstructed phase space theorem, the non-extensive statistical theory of Tsallis, ML techniques, and a technical index, integrating the generated information, which we introduce and named complexity factor (COFA). Our analysis revealed that the size distribution of the genomic regions within chromosomes are not random but follow patterns with characteristic features that have been seen through its complexity character, and it is part of the dynamics of the whole genome. Finally, this picture of dynamics in DNA is recognized using ML tools for clustering, classification, and prediction with high accuracy.

Complex Linkage Disequilibrium Effects in HLA-DPB1 Expression and Molecular Mismatch Analyses of Transplantation Outcomes.

BACKGROUND: HLA molecular mismatch (MM) is a risk factor for de novo donor-specific antibody (dnDSA) development in solid organ transplantation. HLA expression differences have also been associated with adverse outcomes in hematopoietic cell transplantation. We sought to study both MM and expression in assessing dnDSA risk. METHODS: One hundred three HLA-DP-mismatched solid organ transplantation pairs were retrospectively analyzed. MM was computed using amino acids (aa), eplets, and, supplementarily, Grantham/Epstein scores. DPB1 alleles were classified as rs9277534-A (low-expression) or rs9277534-G (high-expression) linked. To determine the associations between risk factors and dnDSA, logistic regression, linkage disequilibrium (LD), and population-based analyses were performed. RESULTS: A high-risk AA:GX (recipient:donor) expression combination (X = A or G) demonstrated strong association with HLA-DP dnDSA (P = 0.001). MM was also associated with HLA-DP dnDSA when evaluated by itself (eplet P = 0.007, aa P = 0.003, Grantham P = 0.005, Epstein P = 0.004). When attempting to determine the relative individual effects of the risk factors in multivariable analysis, only AA:GX expression status retained a strong association (relative risk = 18.6, P = 0.007 with eplet; relative risk = 15.8, P = 0.02 with aa), while MM was no longer significant (eplet P = 0.56, aa P = 0.51). Importantly, these risk factors are correlated, due to LD between the expression-tagging single-nucleotide polymorphism and polymorphisms along HLA-DPB1. CONCLUSIONS: The MM and expression risk factors each appear to be strong predictors of HLA-DP dnDSA and to possess clinical utility; however, these two risk factors are closely correlated. These metrics may represent distinct ways of characterizing a common overlapping dnDSA risk profile, but they are not independent. Further, we demonstrate the importance and detailed implications of LD effects in dnDSA risk assessment and possibly transplantation overall.

Genomic characterization of MICA gene using multiple next generation sequencing platforms: A validation study.

We have developed a protocol regarding the genomic characterization of the MICA gene by next generation sequencing (NGS). The amplicon includes the full length of the gene and is about 13 kb. A total of 156 samples were included in the study. Ninety-seven of these samples were previously characterized at MICA by legacy methods (Sanger or sequence specific oligonucleotide) and were used to evaluate the accuracy, precision, specificity, and sensitivity of the assay. An additional 59 DNA samples of unknown ethnicity volunteers from the United States were only genotyped by NGS. Samples were chosen to contain a diverse set of alleles. Our NGS approach included a first round of sequencing on the Illumina MiSeq platform and a second round of sequencing on the MinION platform by Oxford Nanopore Technology (ONT), on selected samples for the purpose of either characterizing new alleles or setting phase among multiple polymorphisms to resolve ambiguities or generate complete sequence for alleles that were only partially reported in the IMGT/HLA database. Complete consensus sequences were generated for every allele sequenced with ONT, extending from the 5' untranslated region (UTR) to the 3' UTR of the MICA gene. Thirty-two MICA sequences were submitted to the IMGT/HLA database including either new alleles or filling up the gaps (exonic, intronic and/or UTRs) of already reported alleles. Some of the challenges associated with the characterization of these samples are discussed.

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.

Resolving MiSeq-Generated Ambiguities in HLA-DPB1 Typing by Using the Oxford Nanopore Technology.

The technical limitations of current next-generation sequencing technologies, combined with an ever-increasing number of human leukocyte antigen (HLA) alleles, form the basis for the additional ambiguities encountered at an increasing rate in clinical practice. HLA-DPB1 characterization, particularly, generates a significant percentage of ambiguities (25.5%), posing a challenge for accurate and unambiguous HLA-DPB1 genotyping. Phasing of exonic heterozygous positions between exon 2 and all other downstream exons has been the major cause of ambiguities. In this study, the Oxford Nanopore MinION, a third-generation sequencing technology, was used to resolve the phasing. The accurate MiSeq sequencing data, combined with the long reads obtained from the MinION platform, allow for the resolution of the tested ambiguities.

Development of hemolytic paroxysmal nocturnal hemoglobinuria without graft loss following hematopoietic stem cell transplantation for acquired aplastic anemia.

PNH is the most common clonal hematopoietic disorder arising in patients with aAA. PNH is caused by mutations in PIGA, a gene that encodes the catalytic subunit of an enzyme involved in the biosynthesis of GPI anchors, transmembrane glycolipids required for cell surface expression of many proteins. PNH clones likely arise as immune escape mechanisms in aAA by preventing CD1D-restricted T-cell recognition of GPI anchors and GPI-linked autoantigens. Though many patients with aAA treated with IST will develop subclinical PNH clones, only a subset will develop PNH disease, characterized by increased thrombosis, intravascular hemolysis, and potential for severe organ dysfunction. In contrast to IST, allogeneic HSCT for patients with aAA is thought to cure bone marrow aplasia and prevent hematopoietic clonal evolution to PNH. Herein, we present a phenomenon of host-derived PNH disease arising in a patient with aAA many years following MSD-BMT, highlighting the importance of monitoring for this clonal disease in aAA patients with stable mixed donor/recipient chimerism after HSCT. We also provide a literature review for similar occurrences of PNH arising after HSCT.

Next-generation HLA typing of 382 International Histocompatibility Working Group reference B-lymphoblastoid cell lines: Report from the 17th International HLA and Immunogenetics Workshop.

Extended molecular characterization of HLA genes in the IHWG reference B-lymphoblastoid cell lines (B-LCLs) was one of the major goals for the 17th International HLA and Immunogenetics Workshop (IHIW). Although reference B-LCLs have been examined extensively in previous workshops complete high-resolution typing was not completed for all the classical class I and class II HLA genes. To address this, we conducted a single-blind study where select panels of B-LCL genomic DNA samples were distributed to multiple laboratories for HLA genotyping by next-generation sequencing methods. Identical cell panels comprised of 24 and 346 samples were distributed and typed by at least four laboratories in order to derive accurate consensus HLA genotypes. Overall concordance rates calculated at both 2- and 4-field allele-level resolutions ranged from 90.4% to 100%. Concordance for the class I genes ranged from 91.7 to 100%, whereas concordance for class II genes was variable; the lowest observed at HLA-DRB3 (84.2%). At the maximum allele-resolution 78 B-LCLs were defined as homozygous for all 11 loci. We identified 11 novel exon polymorphisms in the entire cell panel. A comparison of the B-LCLs NGS HLA genotypes with the HLA genotypes catalogued in the IPD-IMGT/HLA Database Cell Repository, revealed an overall allele match at 68.4%. Typing discrepancies between the two datasets were mostly due to the lower-resolution historical typing methods resulting in incomplete HLA genotypes for some samples listed in the IPD-IMGT/HLA Database Cell Repository. Our approach of multiple-laboratory NGS HLA typing of the B-LCLs has provided accurate genotyping data. The data generated by the tremendous collaborative efforts of the 17th IHIW participants is useful for updating the current cell and sequence databases and will be a valuable resource for future studies.

Assessing the utilization of high-resolution 2-field HLA typing in solid organ transplantation.

HLA typing in solid organ transplantation (SOT) is necessary for determining HLA-matching status between donor-recipient pairs and assessing patients' anti-HLA antibody profiles. Histocompatibility has traditionally been evaluated based on serologically defined HLA antigens. The evolution of HLA typing and antibody identification technologies, however, has revealed many limitations with using serologic equivalents for assessing compatibility in SOT. The significant improvements to HLA typing introduced by next-generation sequencing (NGS) require an assessment of the impact of this technology on SOT. We have assessed the role of high-resolution 2-field HLA typing (HR-2F) in SOT by retrospectively evaluating NGS-typed pre- and post-SOT cases. HR-2F typing was highly instructive or necessary in 41% (156/385) of the cases. Several pre- and posttransplant scenarios were identified as being better served by HR-2F typing. Five different categories are presented with specific case examples. The experience of another center (Temple University Hospital) is also included, whereby 21% of the cases required HR-2F typing by Sanger sequencing, as supported by other legacy methods, to properly address posttransplant anti-HLA antibody issues.

Characterization of 108 Genomic DNA Reference Materials for 11 Human Leukocyte Antigen Loci: A GeT-RM Collaborative Project.

The highly polymorphic human leukocyte antigen (HLA) genes, located in the human major histocompatibility complex, encode the class I and II antigen-presenting molecules, which are centrally involved in the immune response. HLA typing is used for several clinical applications, such as transplantation, pharmacogenetics, and diagnosis of autoimmune disease. HLA typing is highly complex because of the homology of HLA genes and pseudogenes and the extensive polymorphism in the population. The Centers for Disease Control and Prevention established the Genetic Testing Reference Materials Coordination Program (GeT-RM) in partnership with the genetics community to improve the availability of genomic DNA reference materials necessary for quality assurance of genetic laboratory testing. The GeT-RM together with three clinical laboratories and the Coriell Cell Repositories have characterized genomic DNA obtained from a panel of 108 cell lines for all HLA classic polymorphic loci: HLA-A, B, C, DRB1, DRB3, DRB4, DRB5, DQA1, DQB1, DPA1, and DPB1. The goal was to develop a publicly available and renewable source of well-characterized genomic DNA reference materials to support molecular HLA typing assay development, validation, and verification, quality control, and proficiency testing. These genomic DNA samples are publicly available from the National Institutes of General Medical Science Repository at the Coriell Cell Repositories.

Somatic HLA Mutations Expose the Role of Class I-Mediated Autoimmunity in Aplastic Anemia and its Clonal Complications.

Acquired aplastic anemia (aAA) is an acquired deficiency of early hematopoietic cells, characterized by inadequate blood production, and a predisposition to myelodysplastic syndrome (MDS) and leukemia. Although its exact pathogenesis is unknown, aAA is thought to be driven by Human Leukocyte Antigen (HLA)-restricted T cell immunity, with earlier studies favoring HLA class II-mediated pathways. Using whole exome sequencing (WES), we recently identified two aAA patients with somatic mutations in HLA class I genes. We hypothesized that HLA class I mutations are pathognomonic for autoimmunity in aAA, but were previously underappreciated because the Major Histocompatibility Complex (MHC) region is notoriously difficult to analyze by WES. Using a combination of targeted deep sequencing of HLA class I genes and single nucleotide polymorphism array (SNP-A) genotyping we screened 66 aAA patients for somatic HLA class I loss. We found somatic HLA loss in eleven patients (17%), with thirteen loss-of-function mutations in HLA-A*33:03, HLA-A*68:01, HLA-B*14:02 and HLA-B*40:02 alleles. Three patients had more than one mutation targeting the same HLA allele. Interestingly, HLA-B*14:02 and HLA-B*40:02 were significantly overrepresented in aAA patients, compared to ethnicity-matched controls. Patients who inherited the targeted HLA alleles, regardless of HLA mutation status, had a more severe disease course with more frequent clonal complications as assessed by WES, SNP-A, and metaphase cytogenetics, and more frequent secondary MDS. The finding of recurrent HLA class I mutations provides compelling evidence for a predominant HLA class I-driven autoimmunity in aAA, and establishes a novel link between aAA patients' immunogenetics and clonal evolution.

Assessing a single targeted next generation sequencing for human leukocyte antigen typing protocol for interoperability, as performed by users with variable experience.

BACKGROUND: A simplified protocol for HLA-typing -by NGS, developed for use with the Illumina MiSeq, was performed by technologists with varying NGS experience to assess accuracy and reproducibility. METHODS: Technologists from six laboratories typed the same 16 samples at HLA-A, B, C, DRB1, and DQB1. The protocol includes long range PCR, library preparation, and paired-end 250bp sequencing. Two indexing strategies were employed: locus-specific indexing whereby each locus was tagged uniquely and sample-specific indexing whereby all 5 loci for a sample were pooled prior to library preparation. Sequence analysis was performed with two software packages, Target HLA (Omixon) and NGSengine (GenDx). RESULTS: The average number of sequence reads per library was 387,813; however, analysis was limited to 40,000 reads for locus-indexed libraries and 200,000 reads for sample-indexed libraries resulting in an average depth of coverage of 1444 reads per locus. Sufficient reads for genotype analysis were obtained for 98.4% of libraries. Genotype accuracy was >97% in pooled amplicons and >99% in individual amplicons by both software analysis. Inter-laboratory reproducibility was 99.7% and only cause of discordance was cross-contamination of a single amplicon. CONCLUSIONS: This NGS HLA-typing protocol is simple, reproducible, scalable, highly accurate and amenable to clinical testing.

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.

Generation of Full-Length Class I Human Leukocyte Antigen Gene Consensus Sequences for Novel Allele Characterization.

BACKGROUND: Routine, high-resolution human leukocyte antigen (HLA) genotyping by next generation sequencing within clinical immunogenetics laboratories can now provide the full-length gene sequence characterization of fully phased HLA alleles. This powerful technique provides insights into HLA variation beyond the traditionally characterized antigen recognition domain, providing sequence annotation across the entire gene including untranslated and intronic regions and may be used to characterize novel alleles from massively parallel sequencing runs. METHODS: We evaluated the utility of the Omixon Holotype HLA assay to generate credible, fully phased full-length gene consensus sequences for 50 individuals at major histocompatibility complex, class I, A (HLA-A), HLA-B, and HLA-C loci (300 genotyped alleles in total) to identify and characterize novel class I HLA alleles using our downstream analytical pipeline. RESULTS: Our analysis revealed that 7.7% (23/300) of genotyped class I HLA alleles contain novel polymorphisms. Interestingly, all of the novel alleles identified by our analysis were found to harbor sequence variations within intronic regions of the respective locus. In total our analysis identified 17 unique novel class I HLA alleles from 23 of the 300 genotyped alleles and generated full-length gene sequence annotations for 9 previously incompletely annotated HLA class I allele sequences derived from 14 of the 300 genotyped alleles. CONCLUSIONS: The demonstrated utility of the Omixon Holotype HLA assay in combination with our downstream analytical framework to generate fully phased, full-length gene consensus sequences for the identification and characterization of novel HLA alleles, facilitates the study of HLA polymorphism beyond the antigen recognition domain in human health and disease.