The impact of next-generation sequencing technologies on HLA research

HLA class-I polymorphisms among the Punjabi population of Pakistan: A comparative analysis with country's other ethnic groups

The Punjabi population, constituting over 45 % of the country's total population, holds the highest prevalence in Pakistan. To understand their HLA genetics, we genotyped 389 Punjabi subjects for major Class-I loci using the PCR-SSO Luminex® method. Our study identified a total of 162 alleles, including 41 different HLA-A, 72 HLA-B, and 49 HLA-C alleles. The most common alleles included A*11:01 (14.6 %), A*01:01 (11.8 %), A*24:02 (11.3 %); B*40:06 (13.3 %), B*08:01 (10.9 %), B*51:01 (8.7 %); C*15:02 (15.5 %), C*07:02 (15.3 %), and C*04:01 (10.8 %). However, only locus B showed a significant deviation from HWE. The dominant Class I haplotype was A*24:02-B*40:06-C*15:02, followed by A*11:01-B*40:06-C*15:02, while significant LD was observed between all pairs of HLA loci. A distinct genetic makeup was observed in the Pakistani Punjabis as compared to Indian Punjabis, emphasizing the impact of the Indo-Pak partition and religious choices for marriage. In comparison to country's other ethnic groups, the Pakistani population exhibited 76 different alleles at a low field-resolution, with the Punjabi population having highest polymorphism. Phylogenetic analysis revealed that the Punjabi population is most closely related to the Sindhi population, while both populations sharing ancient connections with the Burusho population. These findings have significant implications for transplantation procedures, personalized medicine, disease susceptibility, and evolutionary studies.

Comparative Assessment of Risk and Turn-Around Time between Sequence-Based Typing and Next-Generation Sequencing for HLA Typing

This study compared laboratory risk and turn-around time (TAT) between sequence-based typing (SBT) and next-generation sequencing (NGS) for human leukocyte antigen (HLA) typing. For risk assessment, we utilized the risk priority number (RPN) score based on failure mode and effect analysis (FMEA) and a risk acceptability matrix (RAM) according to the Clinical Laboratory Standards Institute (CLSI) guidelines (EP23-A). Total TAT was documented for the analytical phase, and hands-on time was defined as manual processes conducted by medical technicians. NGS showed a significantly higher total RPN score than SBT (1169 vs. 465). NGS indicated a higher mean RPN score, indicating elevated severity and detectability scores in comparison to SBT (RPN 23 vs. 12, p = 0.001; severity 5 vs. 3, p = 0.005; detectability 5 vs. 4, p < 0.001, respectively). NGS required a greater number of steps than SBT (44 vs. 25 steps), all of which were acceptable for the RAM. NGS showed a longer total TAT, total hands-on time, and hands-on time per step than SBT (26:47:20 vs. 12:32:06, 03:59:35 vs. 00:47:39, 00:05:13 vs. 00:01:54 hh:mm:ss, respectively). Transitioning from SBT to NGS for HLA typing involves increased risk and an extended TAT. This study underscored the importance of evaluating these factors to optimize laboratory efficiency in HLA typing.

Identification of the new allele HLA-C*04:01:01:186 in a Greek individual using next generation sequencing

The newly discovered HLA-C*04:01:01:186 allele differs from HLA-C*04:01:01:01 by a single nucleotide substitution in intron 3.

Identification of the new allele HLA-A*24:632 in a Greek individual using next generation sequencing

The HLA-Α*24:632 allele differs from HLA-A*24:02:01:01 by a single nucleotide substitution in exon 2.

Frequency distribution of HLA class I and II alleles in Greek population and their significance in orchestrating the National Donor Registry Program

Human leukocyte antigens (HLA) represent one of the most polymorphic systems in humans, responsible for the identification of foreign antigens and the presentation of immune responses. Therefore, HLA is considered to play a major role in human disorders, donor-recipient matching and transplantation outcomes. This study aimed to determine the HLA class I and II alleles and haplotypes in the Greek population. Moreover, a comparative analysis of HLA alleles and haplotype frequencies found in Greek and pooled European populations was also performed to acquire a better knowledge about the HLA alleles distribution. A total number of 1896 healthy individuals were typed for their HLA alleles in the National Tissue Typing Center of Greece. High-resolution HLA typing for the HLA-A, -B, -C and -DR, -DQ, -DP with the use of the next-generation sequencing analysis was performed, followed by data analysis for establishing the HLA allele and haplotype differences. The results of this study showed that the most frequent alleles for the HLA-A were the A*02:01:01 (27.1%), *24:02:01 (14.4%), *01:01:01 (9.3%), for the HLA-B were the B*51:01:01 (15.3%), *18:01:01 (9.7%), *35:01:01 (6.8%) and for the HLA-C were the C*04:01:01 (15.4%), *07:01:01 (13.1%), *12:03:01 (9.6%). For the HLA class II, the most frequent alleles for the HLA-DRB1 were the DRB1*11:04:01 (16.4%), *16:01:01 (11.3%), *11:01:01 (9.5%), for the HLA-DQB1 were the DQB1*03:01:01 (30.5%), *05:02:01 (15.1%), *05:01:01 (10.6%) and for the HLA-DPB1 were the DPB1*04:01:01 (34.8%), *02:01:01 (11.6%), *04:02:01 (7.3%). Additionally, the most frequent haplotypes were the A*02:01:01∼C*07:01:01-B*18:01:01∼DRB1*11:04:01 (2.3%), followed by the A*01:01:01∼C*07:01:01∼B*08:01:01∼DRB1*03:01:01 (2.2%), A*24:02:01∼C*04:01:01∼B*35:02:01∼DRB1*11:04:01 (1.4%) and A*02:01:01∼C*04:01:01∼B*35:01:01-DRB1*14:01:01 (1.2%). The results herein were comparable to those obtained from the pooled European populations. Moreover, these results can be used for the improvement of the donor-recipient matching procedure and to understand better the disease association in Greece.

Advancements in HLA Typing Techniques and Their Impact on Transplantation Medicine

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

Characterization of seven novel HLA-DPA1*01:03:01 non-coding variants by next-generation sequencing

Seven different single nucleotide substitutions in non-coding regions gave rise to novel HLA-DPA1*01:03:01 variants.

Discovery of three HLA-DQB1*03:01:01 variants: HLA-DQB1*03:01:01:54, -DQB1*03:01:01:56 and -DQB1*03:01:01:58

Three nucleotide substitutions in intronic regions give rise to the novel alleles: HLA-DQB1*03:01:01:54, -DQB1*03:01:01:56, -DQB1*03:01:01:58.

Integrative HLA typing of tumor and adjacent normal tissue can reveal insights into the tumor immune response

BACKGROUND

The HLA complex is the most polymorphic region of the human genome, and its improved characterization can help us understand the genetics of human disease as well as the interplay between cancer and the immune system. The main function of HLA genes is to recognize "non-self" antigens and to present them on the cell surface to T cells, which instigate an immune response toward infected or transformed cells. While sequence variation in the antigen-binding groove of HLA may modulate the repertoire of immunogenic antigens presented to T cells, alterations in HLA expression can significantly influence the immune response to pathogens and cancer.

METHODS

RNA sequencing was used here to accurately genotype the HLA region and quantify and compare the level of allele-specific HLA expression in tumors and patient-matched adjacent normal tissue. The computational approach utilized in the study types classical and non-classical Class I and Class II HLA alleles from RNA-seq while simultaneously quantifying allele-specific or personalized HLA expression. The strategy also uses RNA-seq data to infer immune cell infiltration into tumors and the corresponding immune cell composition of matched normal tissue, to reveal potential insights related to T cell and NK cell interactions with tumor HLA alleles.

RESULTS

The genotyping method outperforms existing RNA-seq-based HLA typing tools for Class II HLA genotyping. Further, we demonstrate its potential for studying tumor-immune interactions by applying the method to tumor samples from two different subtypes of breast cancer and their matched normal breast tissue controls.

CONCLUSIONS

The integrative RNA-seq-based HLA typing approach described in the study, coupled with HLA expression analysis, neoantigen prediction and immune cell infiltration, may help increase our understanding of the interplay between a patient's tumor and immune system; and provide further insights into the immune mechanisms that determine a positive or negative outcome following treatment with immunotherapy such as checkpoint blockade.

HLA Genetics for the Human Diseases

Highly polymorphic human leukocyte antigen (HLA) molecules (alleles) expressed by different classical HLA class I and class II genes have crucial roles in the regulation of innate and adaptive immune responses, transplant rejection and in the pathogenesis of numerous infectious and autoimmune diseases. To date, over 35,000 HLA alleles have been published from the IPD-IMGT/HLA database, and specific HLA alleles and HLA haplotypes have been reported to be associated with more than 100 different diseases and phenotypes. Next generation sequencing (NGS) technology developed in recent years has provided breakthroughs in various HLA genomic/gene studies and transplant medicine. In this chapter, we review the current information on the HLA genomic structure and polymorphisms, as well as the genetic context in which numerous disease associations have been identified in this region.

Novel alleles in the era of next-generation sequencing-based HLA typing calls for standardization and policy

Next-Generation Sequencing (NGS) has transformed clinical histocompatibility laboratories through its capacity to provide accurate, high-throughput, high-resolution typing of Human Leukocyte Antigen (HLA) genes, which is critical for transplant safety and success. As this technology becomes widely used for clinical genotyping, histocompatibility laboratories now have an increased capability to identify novel HLA alleles that previously would not be detected using traditional genotyping methods. Standard guidelines for the clinical verification and reporting of novelties in the era of NGS are greatly needed. Here, we describe the experience of a clinical histocompatibility laboratory's use of NGS for HLA genotyping and its management of novel alleles detected in an ethnically-diverse population of British Columbia, Canada. Over a period of 18 months, 3,450 clinical samples collected for the purpose of solid organ or hematopoietic stem cell transplantation were sequenced using NGS. Overall, 29 unique novel alleles were identified at a rate of ∼1.6 per month. The majority of novelties (52%) were detected in the alpha chains of class II (HLA-DQA1 and -DPA1). Novelties were found in all 11 HLA classical genes except for HLA-DRB3, -DRB4, and -DQB1. All novelties were single nucleotide polymorphisms, where more than half led to an amino acid change, and one resulted in a premature stop codon. Missense mutations were evaluated for changes in their amino acid properties to assess the potential effect on the novel HLA protein. All novelties identified were confirmed independently at another accredited HLA laboratory using a different NGS assay and platform to ensure validity in the reporting of novelties. The novel alleles were submitted to the Immuno Polymorphism Database-Immunogenetics/HLA (IPD-IMGT/HLA) for official allele name designation and inclusion in future database releases. A nationwide survey involving all Canadian HLA laboratories confirmed the common occurrence of novel allele detection but identified a wide variability in the assessment and reporting of novelties. In summary, a considerable proportion of novel alleles were identified in routine clinical testing. We propose a framework for the standardization of policies on the clinical management of novel alleles and inclusion in proficiency testing programs in the era of NGS-based HLA genotyping.

Simultaneous detection of DNA variation and methylation at HLA class II locus and immune gene promoters using targeted SureSelect Methyl-Sequencing

The Human Leukocyte Antigen (HLA) locus associates with a variety of complex diseases, particularly autoimmune and inflammatory conditions. The HLA-DR15 haplotype, for example, confers the major risk for developing Multiple Sclerosis in Caucasians, pinpointing an important role in the etiology of this chronic inflammatory disease of the central nervous system. In addition to the protein-coding variants that shape the functional HLA-antigen-T cell interaction, recent studies suggest that the levels of HLA molecule expression, that are epigenetically controlled, also play a role in disease development. However, deciphering the exact molecular mechanisms of the HLA association has been hampered by the tremendous genetic complexity of the locus and a lack of robust approaches to investigate it. Here, we developed a method to specifically enrich the genomic DNA from the HLA class II locus (chr6:32,426,802-34,167,129) and proximal promoters of 2,157 immune-relevant genes, utilizing the Agilent RNA-based SureSelect Methyl-Seq Capture related method, followed by sequencing to detect genetic and epigenetic variation. We demonstrated successful simultaneous detection of the genetic variation and quantification of DNA methylation levels in HLA locus. Moreover, by the detection of differentially methylated positions in promoters of immune-related genes, we identified relevant pathways following stimulation of cells. Taken together, we present a method that can be utilized to study the interplay between genetic variance and epigenetic regulation in the HLA class II region, potentially, in a wide disease context.

The novel HLA-DQB1*03:02:01:14 allele was possibly generated by a recombination event

The novel HLA-DQB1*03:02:01:14 was likely generated by a recombination event between DQB1*03:02:01:01 and DQB1*03:03:02:01.

Detection and characterization of the novel HLA-DPA1*02:66:02N allele, with a premature stop codon in exon 2

The failure to identify HLA null alleles in bone marrow transplantation could be life-threatening because this could result in an HLA mismatch with the ability to trigger the graft-vs-host disease (GVHD) and to reduce patient's survival. In this report we describe the identification and characterization of the novel HLA-DPA1*02:66:02N allele with a non-sense codon in exon 2. This new allele was discovered in two unrelated bone marrow donors during routine HLA-typing using next-generation sequencing (NGS). DPA1*02:66:02N is homologous to DPA1*02:01:01:03 with a single nucleotide difference in exon 2, codon 50, where the replacement of C located at genomic position 3825 by T, causes the formation of a premature stop codon (TGA), resulting in a null allele. This description illustrates the benefits of HLA typing by NGS since it permits to reduce ambiguities, identify new alleles, analyze multiple HLA loci and improve transplantation outcome.

HLA allele-specific expression: Methods, disease associations, and relevance in hematopoietic stem cell transplantation

Varying HLA allele-specific expression levels are associated with human diseases, such as graft versus host disease (GvHD) in hematopoietic stem cell transplantation (HSCT), cytotoxic T cell response and viral load in HIV infection, and the risk of Crohn’s disease. Only recently, RNA-based next generation sequencing (NGS) methodologies with accompanying bioinformatics tools have emerged to quantify HLA allele-specific expression replacing the quantitative PCR (qPCR) -based methods. These novel NGS approaches enable the systematic analysis of the HLA allele-specific expression changes between individuals and between normal and disease phenotypes. Additionally, analyzing HLA allele-specific expression and allele-specific expression loss provide important information for predicting efficacies of novel immune cell therapies. Here, we review available RNA sequencing-based approaches and computational tools for NGS to quantify HLA allele-specific expression. Moreover, we explore recent studies reporting disease associations with differential HLA expression. Finally, we discuss the role of allele-specific expression in HSCT and how considering the expression quantification in recipient-donor matching could improve the outcome of HSCT.

New insights on the role of human leukocyte antigen complex in primary biliary cholangitis

Primary Biliary Cholangitis (PBC) is a rare autoimmune cholangiopathy. Genetic studies have shown that the strongest statistical association with PBC has been mapped in the human leukocyte antigen (HLA) locus, a highly polymorphic area that mostly contribute to the genetic variance of the disease. Furthermore, PBC presents high variability throughout different population groups, which may explain the different geoepidemiology of the disease. A major role in defining HLA genetic contribution has been given by genome-wide association studies (GWAS) studies; more recently, new technologies have been developed to allow a deeper understanding. The study of the altered peptides transcribed by genetic alterations also allowed the development of novel therapeutic strategies in the context of immunotolerance. This review summarizes what is known about the immunogenetics of PBC with a focus on the HLA locus, the different distribution of HLA alleles worldwide, and how HLA modifications are associated with the pathogenesis of PBC. Novel therapeutic strategies are also outlined.

Analysis of the Genomic Sequence of ABO Allele Using Next-Generation Sequencing Method

Background

Although many molecular diagnostic methods have been used for ABO genotyping, there are few reports on the full-length genomic sequence analysis of the ABO gene. Recently, next-generation sequencing (NGS) has been shown to provide fast and high-throughput results and is widely used in the clinical laboratory. Here, we established an NGS method for analyzing the sequence of the start codon to the stop codon in the ABO gene.

Study Design and Methods

Two pairs of primers covering the partial 5’-untranslated region (UTR) to 3’-UTR of the ABO gene were designed. The sequences covering from the start codon to the stop codon of the ABO gene were amplified using these primers, and an NGS method based on the overlap amplicon was developed. A total of 110 individuals, including 88 blood donors with normal phenotypes and 22 ABO subtypes, were recruited and analyzed. All these specimens were first detected by serological tests and then determined by polymerase chain reaction sequence-based typing (PCR-SBT) and NGS. The sequences, including all the intron regions for the specimens, were analyzed by bioinformatics software.

Results

Among the 88 blood donors with a normal phenotype, 48 homozygous individuals, 39 heterozygous individuals, and one individual with a novel O allele were found according to the results of the PCR-SBT method. Some single-nucleotide variants (SNV) in intronic regions were found to be specific for different ABO alleles from 48 homozygous individuals using the NGS method. Sequences in the coding region of all specimens using the NGS method were the same as those of the PCR-SBT method. Three intronic SNVs were found to be associated with the ABO subtypes, including one novel intronic SNV (c.28+5956T&gt;A). Moreover, six specimens were found to exhibit DNA recombination.

Conclusion

An NGS method was established to analyze the sequence from the start codon to the stop codon of the ABO gene. Two novel ABO alleles were identified, and DNA recombination was found to exist in the ABO alleles.

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.

Immune Dysregulation in Autism Spectrum Disorder: What Do We Know about It?

Autism spectrum disorder (ASD) is a group of complex multifactorial neurodevelopmental disorders characterized by a wide and variable set of neuropsychiatric symptoms, including deficits in social communication, narrow and restricted interests, and repetitive behavior. The immune hypothesis is considered to be a major factor contributing to autism pathogenesis, as well as a way to explain the differences of the clinical phenotypes and comorbidities influencing disease course and severity. Evidence highlights a link between immune dysfunction and behavioral traits in autism from several types of evidence found in both cerebrospinal fluid and peripheral blood and their utility to identify autistic subgroups with specific immunophenotypes; underlying behavioral symptoms are also shown. This review summarizes current insights into immune dysfunction in ASD, with particular reference to the impact of immunological factors related to the maternal influence of autism development; comorbidities influencing autism disease course and severity; and others factors with particular relevance, including obesity. Finally, we described main elements of similarities between immunopathology overlapping neurodevelopmental and neurodegenerative disorders, taking as examples autism and Parkinson Disease, respectively.

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.

Personalized HLA typing leads to the discovery of novel HLA alleles and tumor‐specific HLA variants

Accurate and full‐length typing of the HLA region is important in many clinical and research settings. With the advent of next generation sequencing (NGS), several HLA typing algorithms have been developed, including many that are applicable to whole exome sequencing (WES). However, most of these solutions operate by providing the closest‐matched HLA allele among the known alleles in IPD‐IMGT/HLA Database. These database‐matching approaches have demonstrated very high performance when typing well characterized HLA alleles. However, as they rely on the completeness of the HLA database, they are not optimal for detecting novel or less well characterized alleles. Furthermore, the database‐matching approaches are also not adequate in the context of cancer, where a comprehensive characterization of somatic HLA variation and expression patterns of a tumor's HLA locus may guide therapy and clinical outcome, because of the pivotal role HLA alleles play in tumor antigen recognition and immune escape. Here, we describe a personalized HLA typing approach applied to WES data that leverages the strengths of database‐matching approaches while simultaneously allowing for the discovery of novel HLA alleles and tumor‐specific HLA variants, through the systematic integration of germline and somatic variant calling. We applied this approach on WES from 10 metastatic melanoma patients and validated the HLA typing results using HLA targeted NGS sequencing from patients where at least one HLA germline candidate was detected on Class I HLA. Targeted NGS sequencing confirmed 100% performance for the 1st and 2nd fields. In total, five out of the six detected HLA germline variants were because of Class I ambiguities at the third or fourth fields, and their detection recovered the correct HLA allele genotype. The sixth germline variant let to the formal discovery of a novel Class I allele. Finally, we demonstrated a substantially improved somatic variant detection accuracy in HLA alleles with a 91% of success rate in simulated experiments. The approach described here may allow the field to genotype more accurately using WES data, leading to the discovery of novel HLA alleles and help characterize the relationship between somatic variation in the HLA region and immunosurveillance.

The QChip1 knowledgebase and microarray for precision medicine in Qatar

Risk genes for Mendelian (single-gene) disorders (SGDs) are consistent across populations, but pathogenic risk variants that cause SGDs are typically population-private. The goal was to develop "QChip1," an inexpensive genotyping microarray to comprehensively screen newborns, couples, and patients for SGD risk variants in Qatar, a small nation on the Arabian Peninsula with a high degree of consanguinity. Over 108 variants in 8445 Qatari were identified for inclusion in a genotyping array containing 165,695 probes for 83,542 known and potentially pathogenic variants in 3438 SGDs. QChip1 had a concordance with whole-genome sequencing of 99.1%. Testing of QChip1 with 2707 Qatari genomes identified 32,674 risk variants, an average of 134 pathogenic alleles per Qatari genome. The most common pathogenic variants were those causing homocystinuria (1.12% risk allele frequency), and Stargardt disease (2.07%). The majority (85%) of Qatari SGD pathogenic variants were not present in Western populations such as European American, South Asian American, and African American in New York City and European and Afro-Caribbean in Puerto Rico; and only 50% were observed in a broad collection of data across the Greater Middle East including Kuwait, Iran, and United Arab Emirates. This study demonstrates the feasibility of developing accurate screening tools to identify SGD risk variants in understudied populations, and the need for ancestry-specific SGD screening tools.

Using synthetic chromosome controls to evaluate the sequencing of difficult regions within the human genome

BACKGROUND

Next-generation sequencing (NGS) can identify mutations in the human genome that cause disease and has been widely adopted in clinical diagnosis. However, the human genome contains many polymorphic, low-complexity, and repetitive regions that are difficult to sequence and analyze. Despite their difficulty, these regions include many clinically important sequences that can inform the treatment of human diseases and improve the diagnostic yield of NGS.

RESULTS

To evaluate the accuracy by which these difficult regions are analyzed with NGS, we built an in silico decoy chromosome, along with corresponding synthetic DNA reference controls, that encode difficult and clinically important human genome regions, including repeats, microsatellites, HLA genes, and immune receptors. These controls provide a known ground-truth reference against which to measure the performance of diverse sequencing technologies, reagents, and bioinformatic tools. Using this approach, we provide a comprehensive evaluation of short- and long-read sequencing instruments, library preparation methods, and software tools and identify the errors and systematic bias that confound our resolution of these remaining difficult regions.

CONCLUSIONS

This study provides an analytical validation of diagnosis using NGS in difficult regions of the human genome and highlights the challenges that remain to resolve these difficult regions.

HLA imputation and its application to genetic and molecular fine-mapping of the MHC region in autoimmune diseases

Variations of human leukocyte antigen (HLA) genes in the major histocompatibility complex region (MHC) significantly affect the risk of various diseases, especially autoimmune diseases. Fine-mapping of causal variants in this region was challenging due to the difficulty in sequencing and its inapplicability to large cohorts. Thus, HLA imputation, a method to infer HLA types from regional single nucleotide polymorphisms, has been developed and has successfully contributed to MHC fine-mapping of various diseases. Different HLA imputation methods have been developed, each with its own advantages, and recent methods have been improved in terms of accuracy and computational performance. Additionally, advances in HLA reference panels by next-generation sequencing technologies have enabled higher resolution and a more reliable imputation, allowing a finer-grained evaluation of the association between sequence variations and disease risk. Risk-associated variants in the MHC region would affect disease susceptibility through complicated mechanisms including alterations in peripheral responses and central thymic selection of T cells. The cooperation of reliable HLA imputation methods, informative fine-mapping, and experimental validation of the functional significance of MHC variations would be essential for further understanding of the role of the MHC in the immunopathology of autoimmune diseases.

The importance of establishing genetic phase in clinical medicine

Haplotyping or determination of genetic phase has always played a pivotal role in MHC (HLA studies) both in helping to understand inheritance patterns in diseases such as type 1 diabetes (T1D) and in ensuring better matching in transplantation scenarios such as haematopoietic stem cell transplantation (HSCT), using donors genetically related to the patient. In recent years the need to establish genetic phase in a number of clinical scenarios has become apparent. These include: Genetic phasing for hematopoietic stem cell transplants using unrelated donors, where the HLA haplotypes are not known but where haplotype-matched recipients fare better clinically than allele matched, but haplotype mismatched patients. The use of checkpoint inhibitors is one of the most innovative and exciting developments in cancer treatment in years. An example is the use of the monoclonal ipilimumab to block the CTLA-4 receptor which is known to contain polymorphic sites. Until the phase of these polymorphisms is known it will not be possible to determine how effectively this monoclonal will perform in individual patients. The role of miRNA single strand molecules and their effect on gene expression. Thousands of non-coding genes have been identified and have been shown to be polymorphic, as have their target genes. Genetic phasing of polymorphism both in the miRNA source genes and their targets is clearly a fertile area of research In areas such a drug metabolism where the polymorphic family of CYP genes is responsible for the metabolism of the majority of prescription drugs, determining phase of SNPs is critical to understanding drug metabolism and efficacy. In multigenic disease studies combinations of single nucleotide polymorphisms (SNPs) in participating genes require accurate phasing in order to fully appreciate their role in the disease process. In addition, the level of expression of genes (point 3) is also important in understanding disease processes at the functional level. This review outlines the techniques that are currently available for approximating phase and discusses the clinical relevance of establishing genetic phase in areas of clinical medicine outlined in points 1-3.

Approaching Genetics Through the MHC Lens: Tools and Methods for HLA Research

The current SARS-CoV-2 pandemic era launched an immediate and broad response of the research community with studies both about the virus and host genetics. Research in genetics investigated HLA association with COVID-19 based on in silico, population, and individual data. However, they were conducted with variable scale and success; convincing results were mostly obtained with broader whole-genome association studies. Here, we propose a technical review of HLA analysis, including basic HLA knowledge as well as available tools and advice. We notably describe recent algorithms to infer and call HLA genotypes from GWAS SNPs and NGS data, respectively, which opens the possibility to investigate HLA from large datasets without a specific initial focus on this region. We thus hope this overview will empower geneticists who were unfamiliar with HLA to run MHC-focused analyses following the footsteps of the Covid-19|HLA & Immunogenetics Consortium.

PharmVIP: A Web-Based Tool for Pharmacogenomic Variant Analysis and Interpretation

The increasing availability of next generation sequencing (NGS) for personal genomics could promote pharmacogenomics (PGx) discovery and application. However, current tools for analysis and interpretation of pharmacogenomic variants from NGS data are inadequate, as none offer comprehensive analytic functions in a simple, web-based platform. In addition, no tools exist to analyze human leukocyte antigen (HLA) genes for determining potential risks of immune-mediated adverse drug reaction (IM-ADR). We describe PharmVIP, a web-based PGx tool, for one-stop comprehensive analysis and interpretation of genome-wide variants obtained from NGS platforms. PharmVIP comprises three main interpretation modules covering analyses of pharmacogenes involved in pharmacokinetics, pharmacodynamics and IM-ADR. The Guideline module provides Clinical Pharmacogenetics Implementation Consortium (CPIC) drug guideline recommendations based on the translation of genotypic data in genes having guidelines. The HLA module reports HLA genotypes, potential adverse drug reactions, and the relevant drug guidelines. The Pharmacogenes module is employed for prioritizing variants according to variant effect on gene function. Detailed, customizable reports are provided as exportable files and as an interactive web version. PharmVIP is a new integrated NGS workflow for the PGx community to facilitate discovery and clinical application.

A New Human Leukocyte Antigen Typing Algorithm Combined With Currently Available Genotyping Tools Based on Next-Generation Sequencing Data and Guidelines to Select the Most Likely Human Leukocyte Antigen Genotype

Background

High-precision human leukocyte antigen (HLA) genotyping is crucial for anti-cancer immunotherapy, but existing tools predicting HLA genotypes using next-generation sequencing (NGS) data are insufficiently accurate.

Materials and Methods

We compared availability, accuracy, correction score, and complementary ratio of eight HLA genotyping tools (OptiType, HLA-HD, PHLAT, seq2HLA, arcasHLA, HLAscan, HLA*LA, and Kourami) using 1,005 cases from the 1000 Genomes Project data. We created a new HLA-genotyping algorithm combining tools based on the precision and the accuracy of tools’ combinations. Then, we assessed the new algorithm’s performance in 39 in-house samples with normal whole-exome sequencing (WES) data and polymerase chain reaction–sequencing-based typing (PCR-SBT) results.

Results

Regardless of the type of tool, the calls presented by more than six tools concordantly showed high accuracy and precision. The accuracy of the group with at least six concordant calls was 100% (97/97) in HLA-A, 98.2% (112/114) in HLA-B, 97.3% (142/146) in HLA-C. The precision of the group with at least six concordant calls was over 98% in HLA-ABC. We additionally calculated the accuracy of the combination tools considering the complementary ratio of each tool and the accuracy of each tool, and the accuracy was over 98% in all groups with six or more concordant calls. We created a new algorithm that matches the above results. It was to select the HLA type if more than six out of eight tools presented a matched type. Otherwise, determine the HLA type experimentally through PCR-SBT. When we applied the new algorithm to 39 in-house cases, there were more than six matching calls in all HLA-A, B, and C, and the accuracy of these concordant calls was 100%.

Conclusions

HLA genotyping accuracy using NGS data could be increased by combining the current HLA genotyping tools. This new algorithm could also be useful for preliminary screening to decide whether to perform an additional PCR-based experimental method instead of using tools with NGS data.

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.

High throughput analysis of MHC-I and MHC-DR diversity of Brazilian cattle populations

The major histocompatibility complex (MHC) contains many genes that play key roles in initiating and regulating immune responses. This includes the polymorphic MHCI and MHCII genes that present epitopes to CD8+ and CD4+ T-cells, respectively. Consequently, the characterisation of the repertoire of MHC genes is an important component of improving our understanding of the genetic variation that determines the outcomes of immune responses. In cattle, MHC (BoLA) research has predominantly focused on Holstein-Friesian animals (as the most economically important breed globally), although the development of high-throughput approaches has allowed the BoLA-DRB3 repertoire to be studied in a greater variety of breeds. In a previous study we reported on the development of a MiSeq-based method to enable high-throughput and high-resolution analysis of bovine MHCI repertoires. Herein, we report on the expansion of this methodology to incorporate analysis of the BoLA-DRB3 and its application to analyse MHC diversity in a large cohort of cattle from Brazil (>500 animals), including representatives from the three major Bos indicus breeds present in Brazil - Guzerat, Gir and Nelore. This large-scale description of paired MHCI-DRB3 repertoires in Bos indicus cattle has identified a small number of novel DRB3 alleles, a large number of novel MHCI alleles and haplotypes, and provided novel insights into MHCI-MHCII association - further expanding our knowledge of bovine MHC diversity.

Combined impact of HLA-allele matching and the CD34-positive cell dose on optimal unit selection for single-unit cord blood transplantation in adults

The combined effects of HLA-allele matching at six-loci (HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1) and CD34⁺ cell dose on clinical outcomes were analyzed in 1,226 adult cases with single-unit unrelated cord blood transplantation. In the six-loci analysis, low HLA-allele matches did not significantly increase the overall mortality compared to higher matches, whereas in the five-loci analysis excluding HLA-DPB1, they caused a higher overall mortality (HR 1.42, p = .002), possibly due to the graft-versus-leukemia effect of HLA-DPB1 mismatches. A lower CD34⁺ cell dose (<.50 × 10⁵/kg) resulted in higher mortality and lower engraftment; these inferior outcomes were offset by high HLA-allele matches (7-10/10 match), while the inferior outcomes of low HLA-allele matches were improved by increasing the CD34⁺ cell dose. Consideration of the combined effects of the CD34⁺ cell dose and HLA matching may expand the options for transplantable units when HLA matching or the CD34⁺ cell dose is inadequate.

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.

Towards population-specific pharmacogenomics in the era of next-generation sequencing

Pharmacogenomics (PGx) has essential roles in identifying optimal drug responders, optimizing dosage regimens and avoiding adverse events. Population-specific therapeutic interventions that tackle the genetic root causes of clinical outcomes are an important precision medicine strategy. In this perspective, we discuss next-generation sequencing genotyping and its significance for population-specific PGx applications. We emphasize the potential of NGS for preemptive pharmacogenotyping, which is crucial to population-specific clinical studies and patient care. We also provide examples that use publicly available population-based genomics data for population-specific PGx studies. Last, we discuss the remaining challenges and regulatory efforts towards improvements in this field.

Challenges in the application of NGS in the clinical laboratory

Next-generation sequencing (NGS), also known as massively parallel sequencing, has revolutionized genomic research. The current advances in NGS technology make it possible to provide high resolution, high throughput HLA typing in clinical laboratories. The focus of this review is on the recent development and implementation of NGS in clinical laboratories. Here, we examine the critical role of NGS technologies in clinical immunology for HLA genotyping. Two major NGS platforms (Illumina and Ion Torrent) are characterized including NGS library preparation, data analysis, and validation. Challenges of NGS implementation in the clinical laboratory are also discussed, including sequencing error rate, bioinformatics, result interpretation, analytic sensitivity, as well as large data storage. This review aims to promote the broader applications of NGS technology in clinical laboratories and advocate for the novel applications of NGS to drive future research.

HLA RNA Sequencing With Unique Molecular Identifiers Reveals High Allele-Specific Variability in mRNA Expression

The HLA gene complex is the most important single genetic factor in susceptibility to most diseases with autoimmune or autoinflammatory origin and in transplantation matching. Most studies have focused on the vast allelic variation in these genes; only a few studies have explored differences in the expression levels of HLA alleles. In this study, we quantified mRNA expression levels of HLA class I and II genes from peripheral blood samples of 50 healthy individuals. The gene- and allele-specific mRNA expression was assessed using unique molecular identifiers, which enabled PCR bias removal and calculation of the number of original mRNA transcripts. We identified differences in mRNA expression between different HLA genes and alleles. Our results suggest that HLA alleles are differentially expressed and these differences in expression levels are quantifiable using RNA sequencing technology. Our method provides novel insights into HLA research, and it can be applied to quantify expression differences of HLA alleles in various tissues and to evaluate the role of this type of variation in transplantation matching and susceptibility to autoimmune diseases.

Accurate imputation of human leukocyte antigens with CookHLA

The recent development of imputation methods enabled the prediction of human leukocyte antigen (HLA) alleles from intergenic SNP data, allowing studies to fine-map HLA for immune phenotypes. Here we report an accurate HLA imputation method, CookHLA, which has superior imputation accuracy compared to previous methods. CookHLA differs from other approaches in that it locally embeds prediction markers into highly polymorphic exons to account for exonic variability, and in that it adaptively learns the genetic map within MHC from the data to facilitate imputation. Our benchmarking with real datasets shows that our method achieves high imputation accuracy in a wide range of scenarios, including situations where the reference panel is small or ethnically unmatched.

HLA-A, -B, -C, -DRB1 allele and haplotype frequencies of the Korean population and performance characteristics of HLA typing by next-generation sequencing

INTRODUCTION

Human leukocyte antigen (HLA) identification at the allelic level is important for haematopoietic stem cell transplantation (HSCT). Next-generation sequencing (NGS) resolves ambiguous alleles by determining the phase of the polymorphisms. The aim of this study was to validate the software for HLA-SBT (sequence-based typing), assess Korean allele frequency, and characterise the performance of NGS-HLA typing.

METHODS

From the 2009 to 2016 registry, 1293 unrelated healthy donors with a complete dataset of previously characterised HLA-A, -B, -C, and -DRB1 loci were selected and assessed for frequency, haplotype inference, and relative linkage disequilibrium. For performance characteristics of NGS-HLA, alleles included in 1293 cases and ambiguous or alleles assigned as new by SBT-HLA software, or unassigned alleles were included. A total of 91 and 41 quality control samples resulted in 1056 alleles (132 samples × 4 loci × 2 diploid) for analysis. The GenDx NGSgo kit was used for NGS-HLA typing using the Illumina MiSeq platform.

RESULTS

A panel of 132 samples covered 231 alleles, including 53 HLA-A, 80 HLA-B, 43 HLA-C, and 55 HLA-DRB1 by HLA-SBT typing. Comparison of SBT-HLA and NGS-HLA typing showed 99.7% (1053/1056) concordance and discrepant cases were resolved by manual evaluation. Typing by NGS resulted in 67 HLA-A, 112 HLA-B, 71 HLA-C, and 72 HLA-DRB1 alleles. A total of 132 ambiguous, 4 new, and 1 unassigned alleles by HLA-SBT were resolved by NGS-HLA typing.

CONCLUSIONS

NGS-HLA typing provided robust and conclusive results without ambiguities, and its implementation could support HSCT in clinical settings.

Pathogenesis of liver injury in Takayasu arteritis: advanced understanding leads to new horizons

Liver injury in Takayasu arteritis (TA) is a rare phenomenon. Most symptoms are nonspecific, and the exact pathogenesis remains to be elucidated. Early diagnosis and new treatment methods are important for an improved prognosis. A summary of the clinical information and mechanistic analyses may contribute to making an early diagnosis and development of new treatment methods. A PubMed search was conducted using the specific key words “Takayasu arteritis” and “liver” or “hepatitis” or “hepatic”. Symptoms and treatment of TA with an accompanying liver injury were reviewed retrospectively. Many factors are presumed to be involved in the mechanism of TA with liver injury, including the immune response, genes, infections, and gut microbiota. There are several lines of evidence indicating that immune dysfunction is the main pathogenic factor that triggers granuloma formation in TA patients. However, the role of genetics and infections has not been fully confirmed. Recently, the gut microbiota has emerged as an essential component in the process. We reviewed in detail the current concepts that support the complex pathogenesis of TA accompanied by liver injury, and we presented recent theories from the literature. Finally, we discussed future research directions of liver injury in TA.

Analysis of HLA-G long-read genomic sequences in mother-offspring pairs with preeclampsia

Preeclampsia is a pregnancy-induced disorder that is characterized by hypertension and is a leading cause of perinatal and maternal–fetal morbidity and mortality. HLA-G is thought to play important roles in maternal–fetal immune tolerance, and the associations between HLA-G gene polymorphisms and the onset of pregnancy-related diseases have been explored extensively. Because contiguous genomic sequencing is difficult, the association between the HLA-G genotype and preeclampsia onset is controversial. In this study, genomic sequences of the HLA-G region (5.2 kb) from 31 pairs of mother–offspring genomic DNA samples (18 pairs from normal pregnancies/births and 13 from preeclampsia births) were obtained by single-molecule real-time sequencing using the PacBio RS II platform. The HLA-G alleles identified in our cohort matched seven known HLA-G alleles, but we also identified two new HLA-G alleles at the fourth-field resolution and compared them with nucleotide sequences from a public database that consisted of coding sequences that cover the 3.1-kb HLA-G gene span. Intriguingly, a potential association between preeclampsia onset and the poly T stretch within the downstream region of the HLA-G*01:01:01:01 allele was found. Our study suggests that long-read sequencing of HLA-G will provide clues for characterizing HLA-G variants that are involved in the pathophysiology of preeclampsia.

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.

Novel HLA Class I Alleles Outside the Extended DR3 Haplotype Are Protective against Autoimmune Hepatitis

HLA class II allele, DRB1*03:01, is the most common genetic risk factor for autoimmune hepatitis (AIH), but other unrecognized HLA related risks exist.

HLA-B*39:01:01 is a novel risk factor for antithyroid drug-induced agranulocytosis in Japanese population

Antithyroid drug (ATD) is a mainstay of Graves' disease (GD). About 0.1-0.5% of patients with GD treated with ATD exhibit ATD-induced agranulocytosis, which is characterized by severe reduction of circulating neutrophils. Immune-mediated responses have been proposed as a possible mechanism for the pathogenesis of ATD-induced agranulocytosis. Although it has been reported that the HLA class II allele (HLA-DRB1*08:03) was associated with ATD-induced agranulocytosis in multiple populations, the entire HLA region have not been explored in Japanese. Therefore, we performed HLA sequencing for 10 class I and 11 class II genes in 87 patients with ATD-induced agranulocytosis and 384 patients with GD who did not show ATD-induced agranulocytosis. By conducting case-control association studies at the HLA allele and haplotype levels, we replicated the association between HLA-DRB1*08:03:02 and ATD-induced agranulocytosis (P = 5.2 × 10^-7, odds ratio = 2.80), and identified HLA-B*39:01:01 as an independent risk factor (P = 1.4 × 10^-3, odds ratio = 3.35). To verify reproducibility of the novel association of HLA-B*39:01:01, we reanalyzed allele frequency data for HLA-B*39:01:01 from previous case-control association studies. The association of HLA-B*39:01:01 was significantly replicated in Chinese (P = 9.0 × 10^-3), Taiwanese (P = 1.1 × 10^-3), and European populations (P = 5.2 × 10^-4). A meta-analysis combining results from the previous and current studies reinforced evidence of the association between HLA-B*39:01:01 and ATD-induced agranulocytosis (P_meta = 1.2 × 10^-9, pooled OR = 3.66, 95% CI; 2.41-5.57). The results of this study will provide a better understanding of the pathogenesis of ATD-induced agranulocytosis in the context of HLA-mediated hypersensitivity reactions.

SOAPTyping: an open-source and cross-platform tool for sequence-based typing for HLA class I and II alleles

Background

The human leukocyte antigen (HLA) gene family plays a key role in the immune response and thus is crucial in many biomedical and clinical settings. Utilizing Sanger sequencing, the golden standard technology for HLA typing enables accurate identification of HLA alleles in high-resolution. However, only the commercial software, such as uTYPE, SBT-Assign, and SBTEngine, and very few open-source tools could be applied to perform HLA typing based on Sanger sequencing.

Results

We developed a user-friendly, cross-platform and open-source desktop application, known as SOAPTyping, for Sanger-based typing in HLA class I and II alleles. SOAPTyping can produce accurate results with a comprehensible protocol and featured functions. Moreover, SOAPTyping supports a more advanced group-specific sequencing primers (GSSP) module to solve the ambiguous typing results. We used SOAPTyping to analyze 36 samples with known HLA typing from the University of California Los Angeles (UCLA) International HLA DNA Exchange platform and 100 anonymous clinical samples, and the HLA typing results from SOAPTyping are identical to the golden results and 5.5 times faster than commercial software uTYPE, which shows the usability of SOAPTyping.

Conclusions

We introduce the SOAPTyping as the first open-source and cross-platform HLA typing software with the capability of producing high-resolution HLA typing predictions from Sanger sequence data.

Capturing Differential Allele-Level Expression and Genotypes of All Classical HLA Loci and Haplotypes by a New Capture RNA-Seq Method

The highly polymorphic human major histocompatibility complex (MHC) also known as the human leukocyte antigen (HLA) encodes class I and II genes that are the cornerstone of the adaptive immune system. Their unique diversity (>25,000 alleles) might affect the outcome of any transplant, infection, and susceptibility to autoimmune diseases. The recent rapid development of new next-generation sequencing (NGS) methods provides the opportunity to study the influence/correlation of this high level of HLA diversity on allele expression levels in health and disease. Here, we describe the NGS capture RNA-Seq method that we developed for genotyping all 12 classical HLA loci (HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, HLA-DRA, HLA-DRB1, HLA-DRB3, HLA-DRB4, and HLA-DRB5) and assessing their allelic imbalance by quantifying their allele RNA levels. This is a target enrichment method where total RNA is converted to a sequencing-ready complementary DNA (cDNA) library and hybridized to a complex pool of RNA-specific HLA biotinylated oligonucleotide capture probes, prior to NGS. This method was applied to 161 peripheral blood mononuclear cells and 48 umbilical cord blood cells of healthy donors. The differential allelic expression of 10 HLA loci (except for HLA-DRA and HLA-DPA1) showed strong significant differences (P < 2.1 × 10-15). The results were corroborated by independent methods. This newly developed NGS method could be applied to a wide range of biological and medical questions including graft rejections and HLA-related diseases.

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.

Identifying the HLA DRB1-DQB1 molecules and predicting epitopes associated with high-risk HPV infection clearance and redetection

Several determining factors are involved in HPV infection outcomes; human leukocyte antigen (HLA) polymorphisms have been described as related factors. This study has ascertained the effect of genetic variation on HLA-DRB1 and DQB1 genes on HPV-16/-18/-31/-33/-45 and -58 clearance and redetection in Colombian women. PCR and qPCR were used for viral identification and the Illumina MiSeq system was used for HLA-typing of cervical samples (n = 276). Survival models were adjusted for identifying alleles/haplotypes related to HPV clearance/redetection; L1/L2 protein-epitope binding to MHC-II molecules was also predicted. Significant associations suggested effects favouring or hampering clearance/redetection events depending on the viral type involved in infection, e.g. just DRB1*12:01:01G favoured HPV-16 (coeff: 4.8) and HPV-45 clearance (coeff: 12.65) whilst HPV-18 (coeff: 2E-15), HPV-31 (coeff: 8E-17) and HPV-58 hindered elimination (coeff: 1E-14). An effect was only observed for some alelles when configured as haplotypes, e.g. DRB1*04:07:01G (having the greatest frequency in the target population) was associated with DQB1*02:01:1G or *03:02:03. Epitope prediction identified 23 clearance-related peptides and 29 were redetection-related; eight might have been related to HPV-16/-18 and -58 persistence and one to HPV-18 elimination. HLA allele/haplotype relationship with the course of HPV infection (clearance/redetection) depended on the infecting HPV type, in line with the specific viral epitopes displayed.