Next-generation sequencing reveals new information about HLA allele and haplotype diversity in a large European American population

Divergent HLA variations and heterogeneous expression but recurrent HLA loss-of- heterozygosity and common HLA-B and TAP transcriptional silencing across advanced pediatric solid cancers

The human leukocyte antigen (HLA) system is a major factor controlling cancer immunosurveillance and response to immunotherapy, yet its status in pediatric cancers remains fragmentary. We determined high-confidence HLA genotypes in 576 children, adolescents and young adults with recurrent/refractory solid tumors from the MOSCATO-01 and MAPPYACTS trials, using normal and tumor whole exome and RNA sequencing data and benchmarked algorithms. There was no evidence for narrowed HLA allelic diversity but discordant homozygosity and allele frequencies across tumor types and subtypes, such as in embryonal and alveolar rhabdomyosarcoma, neuroblastoma MYCN and 11q subtypes, and high-grade glioma, and several alleles may represent protective or susceptibility factors to specific pediatric solid cancers. There was a paucity of somatic mutations in HLA and antigen processing and presentation (APP) genes in most tumors, except in cases with mismatch repair deficiency or genetic instability. The prevalence of loss-of-heterozygosity (LOH) ranged from 5.9 to 7.7% in HLA class I and 8.0 to 16.7% in HLA class II genes, but was widely increased in osteosarcoma and glioblastoma (~15-25%), and for DRB1-DQA1-DQB1 in Ewing sarcoma (~23-28%) and low-grade glioma (~33-50%). HLA class I and HLA-DR antigen expression was assessed in 194 tumors and 44 patient-derived xenografts (PDXs) by immunochemistry, and class I and APP transcript levels quantified in PDXs by RT-qPCR. We confirmed that HLA class I antigen expression is heterogeneous in advanced pediatric solid tumors, with class I loss commonly associated with the transcriptional downregulation of HLA-B and transporter associated with antigen processing (TAP) genes, whereas class II antigen expression is scarce on tumor cells and occurs on immune infiltrating cells. Patients with tumors expressing sufficient HLA class I and TAP levels such as some glioma, osteosarcoma, Ewing sarcoma and non-rhabdomyosarcoma soft-tissue sarcoma cases may more likely benefit from T cell-based approaches, whereas strategies to upregulate HLA expression, to expand the immunopeptidome, and to target TAP-independent epitopes or possibly LOH might provide novel therapeutic opportunities in others. The consequences of HLA class II expression by immune cells remain to be established. Immunogenetic profiling should be implemented in routine to inform immunotherapy trials for precision medicine of pediatric cancers.

Report from the extended HLA-DPA1 ~ promoter ~ HLA-DPB1 haplotype of the 18th international HLA and immunogenetics workshop

The primary goal of the HLA-DPA1 ~ promoter ~ HLA-DPB1 haplotype component of the 18th IHIWS was to characterise the extended haplotypes within the HLA-DP region and survey the extent of genetic diversity in this region across human populations. In this report, we analysed single-nucleotide polymorphisms (SNPs) in 255 subjects from 6 different cohorts. The results from the HLA-DP haplotype component have validated findings from the initial pilot study. SNPs in this region were inherited in strong linkage, particularly HLA-DPA1, SNP-linked promoter haplotypes and motifs in exon 2 of HLA-DPB1. We reported 17 SNP-linked haplotypes in the promoter region. Together with HLA-DPA1 and HLA-DPB1 alleles, they formed 74 distinct extended HLA-DP haplotypes in 438 sequences. We also observed the presence of region-specific alleles and promoter haplotypes. Our approach involved phasing extended SNPs including promoter SNPs, HLA-DPA1 and HLA-DPB1 alleles, in a 22 kb region, GRCh38/hg38 (chr6:33,064,111-33,086,679), followed by clustering of these SNPs as one extended haplotype. This hierarchical clustering revealed four major clades, suggesting that haplotypes within each clade may have diverged from a common ancestral haplotype and undergone similar evolutionary processes. The correlation between HLA-DPA1 and the promoter region raises questions about the role of HLA-DPA1 antigen in the heterodimer. This finding requires validation on a larger sample size specifically designed for anthropological analysis. Nevertheless, the results from this study highlight the clinical potential of selecting better-matched donors for patients awaiting haematopoietic stem cell transplants from genetically overlapping groups that share common ancestral haplotypes.

Impact of Preformed Donor-Specific Anti-HLA-Cw and Anti-HLA-DP Antibodies on Acute Antibody-Mediated Rejection in Kidney Transplantation

Given the risk of rejection, the presence of preformed donor specific antibodies (DSA) contraindicates transplantation in most allocation systems. However, HLA-Cw and -DP DSA escape this censorship. We performed a multicentric observational study, in which the objective was to determinate risk factors of acute antibody-mediated rejection (aABMR) in recipients transplanted with preformed isolated Cw- or DP-DSA. Between 2010 and 2019, 183 patients were transplanted with a preformed isolated Cw- or DP-DSA (92 Cw-DSA; 91 DP-DSA). At 2 years, the incidence of aABMR was 12% in the Cw-DSA group, versus 28% in the DP-DSA group. Using multivariable Cox regression model, the presence of a preformed DP-DSA was associated with an increased risk of aABMR (HR = 2.32 [1.21-4.45 (p = 0.001)]) compared with Cw-DSA. We also observed a significant association between the DSA's MFI on the day of transplant and the risk of aABMR (HR = 1.09 [1.08-1.18], p = 0.032), whatever the DSA was. Interaction term analysis found an increased risk of aABMR in the DP-DSA group compared with Cw-DSA, but only for MFI below 3,000. These results may plead for taking these antibodies into account in the allocation algorithms, in the same way as other DSA.

The natural history of de novo donor-specific HLA antibodies after kidney transplantation

Background

De novo donor-specific HLA antibodies (dnDSA) are key factors in the diagnosis of antibody-mediated rejection (ABMR) and related to graft loss.

Methods

This retrospective study was designed to evaluate the natural course of dnDSA in graft function and kidney allograft survival and to assess the impact of mean fluorescence intensity (MFI) evolution as detected by annual Luminex^® screening. All 400 kidney transplant recipients with 731 dnDSA against the last graft (01/03/2000-31/05/2021) were included.

Results

During 8.3 years of follow-up, ABMR occurred in 24.8% and graft loss in 33.3% of the cases, especially in patients with class I and II dnDSA, and those with multiple dnDSA. We observed frequent changes in MFI with 5-year allograft survivals post-dnDSA of 74.0% in patients with MFI reduction ≥ 50%, 62.4% with fluctuating MFI (MFI reduction ≥ 50% and doubling), and 52.7% with doubling MFI (log-rank p < 0.001). Interestingly, dnDSA in 168 (24.3%) cases became negative at some point during follow-up, and 38/400 (9.5%) patients became stable negative, which was associated with better graft survival. Multivariable analysis revealed the importance of MFI evolution and rejection, while class and number of dnDSA were not contributors in this model.

Conclusion

In summary, we provide an in-depth analysis of the natural course of dnDSA after kidney transplantation, first evidence for the impact of MFI evolution on graft outcomes, and describe a relevant number of patients with a stable disappearance of dnDSA, related to better allograft survival.

Distribution of HLA-B Alleles and Haplotypes in Qatari: Recommendation for Establishing Pharmacogenomic Markers Screening for Drug Hypersensitivity

Human leukocyte antigen (HLA) proteins are present at the cellular surface of antigen-presenting cells and play a crucial role in the adaptive immune response. Class I genes, specifically certain HLA-B alleles, are associated with adverse drug reactions (ADRs) and are used as pharmacogenetic markers. Although ADRs are a common causes of hospitalization and mortality, the data on the prevalence of HLA-B pharmacogenetics markers in Arab countries are scarce. In this study, we investigated the frequencies of major HLA-B pharmacogenomics markers in the Qatari population. Next-generation sequencing data from 1,098 Qatari individuals were employed for HLA-B typing using HLA-HD version 1.4.0 and IPD-IMGT/HLA database. In addition, HLA-B pharmacogenetics markers were obtained from the HLA Adverse Drug Reaction Database. In total, 469 major HLA-B pharmacogenetic markers were identified, with HLA-B*51:01 being the most frequent pharmacogenetic marker (26.67%) in the Qatari population. Moreover, HLA-B*51:01 is associated with phenytoin- and clindamycin-induced ADRs. The second most frequent pharmacogenetic marker was the HLA-B*58:01 allele (6.56%), which is associated with allopurinol-induced ADRs. The third most frequent pharmacogenetic marker was the HLA-B*44:03 allele, which is associated with phenytoin-induced ADRs. The establishment of a pharmacogenetics screening program in Qatar for cost effective interventions aimed at preventing drug-induced hypersensitivity can be aided by the highly prevalent HLA-B pharmacogenetic markers detected here.

Allele and Haplotype Frequencies of HLA-A, -B, -C, and -DRB1 Genes in 3,750 Cord Blood Units From a Kinh Vietnamese Population

The frequencies and diversities of human leukocyte antigen (HLA) alleles and haplotypes are representative of ethnicities. Matching HLA alleles is essential for many clinical applications, including blood transfusion, stem cell transplantation, and tissue/organ transplantation. To date, the information about the frequencies and distributions of HLA alleles and haplotypes among the Kinh Vietnamese population is limited because of the small sample size. In this study, more than 3,750 cord blood units from individuals belonging to the Kinh Vietnamese population were genotyped using PCR sequence-specific oligonucleotide (PCR-SSO) for HLA testing. The results of the study demonstrated that the most frequently occurring HLA-A, -B, -C, and -DRB1 alleles were A*11:01 (25%), A*24:02 (12.3%), A*02:01 (11.2); A*03:03 (8.95%), A*02:03 (7.81%), A*29:01 (7.03%); B*15:02 (15.1%), B*46:01 (10.7%), B*58:01 (7.65%), B*38:02 (7.29%); C*08:01 (17.2), C*07:02 (16.2%), C*01:02 (15.2), C*03:02 (8.3%), C*15:05 (6.13); DRB1*12:02 (31.0%), DRB1*09:01 (10.47%), DRB1*15:02 (7.54%); DRB1*07:01 (6.68%), DRB1*10:01 (6.63%), respectively, with the highest allele diversity level observed in locus B (93 alleles). The most frequent haplotypes of two-locus combinations of HLA-A–B, HLA-A–C, HLA-A–DRB1, HLA-B–C, HLA-B–DRB1, and HLA-C–DRB1 haplotypes were A*11:01–B*15:02 (7.63%), A*11:01–C*08:01 (7.98%), A*11:01–DRB1*12:02 (10.56%), B*15:02–C*08:01 (14.0%), B*15:02–DRB1*12:02 (10.47%), and C*08:01–DRB1*12:02 (11.38%), respectively. In addition, the most frequent haplotypes of three- and four-locus sets of HLA-A–B–C, HLA-A–B–DRB1, HLA-A–C–DRB1, HLA-B–C–DRB1, and HLA-A–B–C–DRB1 were A*11:01–B*15:02–C*08:01 (7.57%), A*11:01–B*15:02–DRB1*12:02 (5.39%), A*11:01–C*08:01–DRB1*12:02 (5.54%), B*15:02–C*08:01–DRB1*12:02 (10.21%), and A*11:01–B*15:02–C*08:01–DRB1*12:02 (5.45%), respectively. This study provides critical information on the frequencies and distributions of HLA alleles and haplotypes in the Kinh Vietnamese population, accounting for more than 85% of Vietnamese citizens. It paves the way to establish an umbilical cord blood bank for cord blood transplantation programs in Vietnam.

HLA Diversity in Saudi Population: High Frequency of Homozygous HLA Alleles and Haplotypes

Human leukocyte antigens (HLA) diversity has a tremendous impact on shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host-pathogen interactions. Here, we conducted a retrospective study of HLA class I and class II homozygosity at allelic and haplotype levels in unrelated individuals genotyped from 2012 to 2016 in a tertiary hospital in the capital of Saudi Arabia. Among 5,000 individuals, 2,773 individuals meet inclusion criteria and were retrospectively analyzed for HLA-A, -B, -C-DRB1, and -DQB1 homozygosity at allelic and haplotype levels. HLA molecular typing was performed using a commercial reverse sequence-specific oligonucleotide (rSSO) kit. We were able to identify 15 HLA-A, 20 HLA-B, 11 HLA-C, 13 HLA-DRB1, and five HLA-DQB1 homozygous alleles demonstrating a very low genetic diversity in the Saudi population. The highest homozygosity in HLA class I was found in locus C followed by A and B (20.3% > 16.1% > 15.5%; p < 0.001) where the most homozygote alleles were A*02 (9.2%), B*51 and B*50 (5.7% and 3.7%), and C*07, C*06, and C*15 (7.2%, 5.48%, and 3.3%) and in HLA class II, the highest homozygosity was found in locus DQB1 compared to DRB1 (31.71% > 19.2%; p < 0.001), with the most common homozygote alleles being DRB1*07 and DRB1*04 (5.33% and 4.2%) and DQB1*02, DQB1*06, and DQB1*03 (13.55%, 7.92%, and 7.64%). The frequency of finding an individual with one homozygote allele was (24.6%), two homozygote alleles (13.5%), three homozygote alleles (4.7%), four homozygote alleles (3.4%), and five alleles were (4.8%). The most frequent homozygote haplotypes are A*23∼C*06∼B*50∼DRB1*07∼DQB1*02 and A*02∼C*06∼B*50∼DRB1*07∼DQB1*02. This study shows low diversity of both class I and II alleles and haplotypes in the Saudi population, which would have a significant impact on shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host-pathogen interactions.

HLA-A*03, the hemochromatosis ancestral haplotype, and phenotypes of referred hemochromatosis probands with HFE p.C282Y homozygosity

Background

Human leukocyte antigen (HLA)-A*03, hemochromatosis ancestral haplotype marker, was associated with greater iron overload in hemochromatosis cohorts reported before discovery of the HFE gene. We sought to learn whether an A*03-linked locus influences phenotypes in referred HFE p.C282Y homozygotes.

Methods

We tabulated these phenotypes in probands with p.C282Y homozygosity: age, transferrin saturation (TS), serum ferritin (SF), conditions related to iron overload, fibrosis-four variables (FIB-4) index and aspartate aminotransferase-to-platelet ratio index (APRI) predictors of severe hepatic fibrosis, and iron removed to achieve depletion (QFe/age). We analyzed phenotypes of men and women separately across three A*03 subgroups.

Results

There were 104 men (57.8%) and 76 women (42.2%). Mean age (SD) was 49 ± 13 y. Mean TS was 79 ± 17%. Median SF (range) was 715 µg/L (28, 6103). Related conditions included: hemochromatosis arthropathy (21.7%); type 2 diabetes (18.9%); hypogonadotropic hypogonadism (5.8% of men); cardiomyopathy (0%); and cirrhosis (10.0%). Median QFe/age was 61 mg/y (0, 714). A*03 homozygosity, heterozygosity, and no A*03 occurred in 37 (20.6%), 104 (57.8%), and 39 probands (21.7%), respectively. In men, mean TS and median SF were significantly higher in A*03 homozygotes than heterozygotes but not A*03-negative probands. In men, median APRI was significantly lower in A*03 heterozygotes than homozygotes and A*03-negative probands. No other phenotypes, including QFe/age, differed significantly across A*03 subgroups in either men or women.

Conclusions

Our results suggest that an A*03-linked locus does not influence phenotypes in referred HFE p.C282Y homozygotes. It is unlikely that heritable factors that modify phenotypes of p.C282Y homozygotes are linked to the hemochromatosis ancestral haplotype.

Generation of individualized immunocompatible endothelial cells from HLA-I-matched human pluripotent stem cells

BACKGROUND

Endothelial cells (ECs) derived from human-induced pluripotent stem cell (iPSC) are a valuable cell resource for cardiovascular regeneration. To avoid time-consuming preparation from primary autologous cells, the allogeneic iPSC-ECs are being expected to become "off-the-shelf" cell products. However, allorejection caused by HLA mismatching is a major barrier for this strategy. Although the "hypoimmunogenic" iPSCs could be simply generated by inhibition of HLA-I expression via β-2 microglobulin knockout (B2M KO), the deletion of HLA-I expression will activate natural killer (NK) cells, which kill the HLA-I negative cells. To inhibit NK activation, we proposed to generate HLA-matched iPSCs based on patient's HLA genotyping by HLA exchanging approach to express the required HLA allele.

METHODS

To establish a prototype of HLA exchanging system, the expression of HLA-I molecules of iPSCs was inhibited by CRISPR/Cas9-mediated B2M KO, and then HLA-A*11:01 allele, as a model molecule, was introduced into B2M KO iPSCs by lentiviral gene transfer. HLA-I-modified iPSCs were tested for their pluripotency and ability to differentiate into ECs. The stimulation of iPSC-EC to allogeneic T and NK cells was detected by respective co-culture of PBMC-EC and NK-EC. Finally, the iPSC-ECs were used as the seeding cells to re-endothelialize the decellularized valves.

RESULTS

We generated the iPSCs only expressed one HLA-A allele (HLA-A *11:01) by B2M KO plus HLA gene transfer. These HLA-I-modified iPSCs maintained pluripotency and furthermore were successfully differentiated into functional ECs assessed by tube formation assay. Single HLA-A*11:01-matched iPSC-ECs significantly less induced the allogeneic response of CD8+ T cell and NK cells expressing matched HLA-A*11:01 and other HLA-A,-B and -C alleles. These cells were successfully used to re-endothelialize the decellularized valves.

CONCLUSIONS

In summary, a simple HLA-I exchanging system has been created by efficient HLA engineering of iPSCs to evade both of the alloresponse of CD8+ T cells and the activation of NK cells. This technology has been applied to generate iPSC-ECs for the engineering of cellular heart valves. Our strategy should be extremely useful if the "off-the-shelf" and "non-immunogenic" allogeneic iPSCs were created for the common HLA alleles.

The polymorphism of HLA-A, -C, -B, -DRB3/4/5, -DRB1, -DQB1 loci in Zhejiang Han population, China using NGS technology

The distributions of HLA allele and haplotype are various in the populations. Currently, the data for HLA alleles and haplotypes at three fields resolution level in Chinese Han population is rare. Here, the HLA alleles and haplotypes of the 1734 cord blood samples from Zhejiang Han population, China were reported at three fields resolution. All samples were randomly collected from the Zhejiang Cord Blood Bank, China. HLA-A, -B, -C, -DRB1, -DQB1, -DRB3/4/5 loci was genotyped using next generation sequencing method. The genotypes of the samples were assigned using the HLA TypeStream Visual Software version 1.2.0. The frequency of alleles, haplotype estimation and linkage disequilibrium analysis were performed with the Arlequin software 3.5.2.2. It was found that the top three frequent alleles of HLA-A, -B, -C, -DRB1, -DQB1, -DRB3/4/5 loci were A*11:01:01 (25.81%), A*24:02:01 (16.70%), A*02:01:01 (10.61%); B*40:01:02 (15.97%), B*46:01:01 (11.48%), B*58:01:01 (7.96%); C*07:02:01 (19.03%), C*01:02:01 (17.65%), C*03:04:01 (10.41%); DRB1*09:01:02G (17.96%), DRB1*12:02:01 (9.57%), DRB1*08:03:02 (9.54%); DQB1*03:01:01G (21.05%), DQB1*03:03:02 (19.15%), DQB1*06:01:01G (12.08%); DRB4*01:03:01 (25.72%), DRB3*02:02:01 (20.27%), DRB5*01:01:01 (10.96%), respectively. A total of 1528 distinct A∼C∼B∼DRB3/4/5∼DRB1∼DQB1 haplotypes were estimated, and the top three most common haplotypes were A*33:03:01∼C*03:02:02∼B*58:01:01∼DRB3*02:02:01∼DRB1*03:01:01∼ DQB1*02:01:01 (4.02%), A*30:01:01∼C*06:02:01∼B*13:02:01∼DRB4*01:03:01∼ DRB1*07:01:01 ∼DQB1*02:02:01 (3.11%) and A*02:07:01∼C*01:02:01∼B*46:01:01 ∼DRB4*01:03:01∼DRB1*09:01:02G∼DQB1*03:03:02 (3.05%). Some alleles of different HLA loci were shown strong linkage disequilibrium. In conclusion, the data of allele and haplotype of HLA-A, -B, -C, -DRB1, -DQB1 and -DRB3/4/5 loci at three fields resolution level were obtained in Zhejiang Han population, thus contributing to analyze the HLA ploymorphism in the populations.

Computational Design and Analysis of a Multi-epitope Against Influenza A virus

Influenza A viruses are among the most studied viruses, however no effective prevention against influenza infection has been developed. So, designing an effective vaccine against Influenza A virus is a critical issue in the field of medical biotechnology. For this reason, to combat this disease, we have designed a novel multi-epitope vaccine candidate based on the several conserved and potential linear B-cell and T-cell binding epitopes by using the in silico approach. This vaccine consists of an ER signal conserved sequence, the PADRE conserved epitope and two conserved epitopes of Influenza matrix protein 2. T-cell binding epitopes from Matrix protein 2 were predicted by in silico tools of epitope prediction. The selected epitopes were joined by flexible linkers and physicochemical properties, toxicity, and allergenecity were investigated. The designed vaccine was antigenic, immunogenic, and non-allergenic with suitable physicochemical properties and has higher solubility. The final multi-epitope construct was modeled, confirmed by different programs and the molecular interactions with immune receptors were considered. The molecular docking assay indicated the interactions with immune-stimulatory toll-like receptor 3 (TLR3) and major histocompatibility complex class I (MHCI). The HADDOCK and H DOCK servers were used to make docking analysis, respectively. The docking analysis indicated a strong and stable binding interaction between the vaccine construct with major histocompatibility complex (MHC) class I and toll-like receptor 3. Overall, the findings suggest that the current vaccine may be a promising vaccine to prevent Influenza infection.

Challenges for the standardized reporting of NGS HLA genotyping: Surveying gaps between clinical and research laboratories

Next generation sequencing (NGS) is being applied for HLA typing in research and clinical settings. NGS HLA typing has made it feasible to sequence exons, introns and untranslated regions simultaneously, with significantly reduced labor and reagent cost per sample, rapid turnaround time, and improved HLA genotype accuracy. NGS technologies bring challenges for cost-effective computation, data processing and exchange of NGS-based HLA data. To address these challenges, guidelines and specifications such as Genotype List (GL) String, Minimum Information for Reporting Immunogenomic NGS Genotyping (MIRING), and Histoimmunogenetics Markup Language (HML) were proposed to streamline and standardize reporting of HLA genotypes. As part of the 17th International HLA and Immunogenetics Workshop (IHIW), we implemented standards and systems for HLA genotype reporting that included GL String, MIRING and HML, and found that misunderstanding or misinterpretations of these standards led to inconsistencies in the reporting of NGS HLA genotyping results. This may be due in part to a historical lack of centralized data reporting standards in the histocompatibility and immunogenetics community. We have worked with software and database developers, clinicians and scientists to address these issues in a collaborative fashion as part of the Data Standard Hackathons (DaSH) for NGS. Here we report several categories of challenges to the consistent exchange of NGS HLA genotyping data we have observed. We hope to address these challenges in future DaSH for NGS efforts.

Genome-wide association study reveals an association between the HLA-DPB1∗02:01:02 allele and wheat-dependent exercise-induced anaphylaxis

Wheat-dependent exercise-induced anaphylaxis (WDEIA) is a life-threatening food allergy triggered by wheat in combination with the second factor such as exercise. The identification of potential genetic risk factors for this allergy might help high-risk individuals before consuming wheat-containing food. We aimed to identify genetic variants associated with WDEIA. A genome-wide association study was conducted in a discovery set of 77 individuals with WDEIA and 924 control subjects via three genetic models. The associations were confirmed in a replication set of 91 affected individuals and 435 control individuals. Summary statistics from the combined set were analyzed by meta-analysis with a random-effect model. In the discovery set, a locus on chromosome 6, rs9277630, was associated with WDEIA in the dominant model (OR = 3.95 [95% CI, 2.31-6.73], p = 7.87 × 10-8). The HLA-DPB1∗02:01:02 allele displayed the most significant association with WDEIA (OR = 4.51 [95% CI, 2.66-7.63], p = 2.28 × 10-9), as determined via HLA imputation following targeted sequencing. The association of the allele with WDEIA was confirmed in replication samples (OR = 3.82 [95% CI, 2.33-6.26], p = 3.03 × 10-8). A meta-analysis performed in the combined set revealed that the HLA-DPB1∗02:01:02 allele was significantly associated with an increased risk of WDEIA (OR = 4.13 [95% CI, 2.89-5.93], p = 1.06 × 10-14). Individuals carrying the HLA-DPB1∗02:01:02 allele have a significantly increased risk of WDEIA. Further validation of these findings in independent multiethnic cohorts is needed.

Relevance of Polymorphic KIR and HLA Class I Genes in NK-Cell-Based Immunotherapies for Adult Leukemic Patients

Simple Summary

Immunotherapies are promising approaches to curing different acute leukemias. Natural killer (NK) cells are lymphocytes that are efficient in the elimination of leukemic cells. NK-cell-based immunotherapies are particularly attractive, but the landscape of the heterogeneity of NK cells must be deciphered. This review provides an overview of the polymorphic KIR and HLA class I genes that modulate the NK cell repertoire and how these markers can improve the outcomes of patients with acute leukemia. A better knowledge of these genetic markers that are linked to NK cell subsets that are efficient against hematological diseases will optimize hematopoietic stem-cell donor selection and NK immunotherapy design.

Abstract

Since the mid-1990s, the biology and functions of natural killer (NK) cells have been deeply investigated in healthy individuals and in people with diseases. These effector cells play a particularly crucial role after allogeneic hematopoietic stem-cell transplantation (HSCT) through their graft-versus-leukemia (GvL) effect, which is mainly mediated through polymorphic killer-cell immunoglobulin-like receptors (KIRs) and their cognates, HLA class I ligands. In this review, we present how KIRs and HLA class I ligands modulate the structural formation and the functional education of NK cells. In particular, we decipher the current knowledge about the extent of KIR and HLA class I gene polymorphisms, as well as their expression, interaction, and functional impact on the KIR⁺ NK cell repertoire in a physiological context and in a leukemic context. In addition, we present the impact of NK cell alloreactivity on the outcomes of HSCT in adult patients with acute leukemia, as well as a description of genetic models of KIRs and NK cell reconstitution, with a focus on emergent T-cell-repleted haplo-identical HSCT using cyclosphosphamide post-grafting (haplo-PTCy). Then, we document how the immunogenetics of KIR/HLA and the immunobiology of NK cells could improve the relapse incidence after haplo-PTCy. Ultimately, we review the emerging NK-cell-based immunotherapies for leukemic patients in addition to HSCT.

Tri-SNP polymorphism in the intron of HLA-DRA1 affects type 1 diabetes susceptibility in the Finnish population

Genes in the HLA class II region include the most important inherited risk factors for type 1 diabetes (T1D) although also polymorphisms outside the HLA region modulate the predisposition to T1D. This study set out to confirm a recent observation in which a novel expression quantitative trait locus was formed by three single nucleotide polymorphisms (SNP) in the intron of HLA-DRA1 in DR3-DQ2 haplotypes. The SNPs significantly increased the risk for T1D in DR3-DQ2 homozygous individuals and we intended to further explore this association, in the Finnish population, by comparing two DR3-DQ2 positive genotypes. Cohorts with DR3-DQ2/DR3-DQ2 (N = 570) and DR3-DQ2/DR1-DQ5 (N = 1035) genotypes were studied using TaqMan analysis that typed for rs3135394, rs9268645 and rs3129877. The tri-SNP haplotype was significantly more common in cases than controls in the DR3-DQ2/DR3-DQ2 cohort (OR = 1.70 CI 95% = 1.15-2.51P = 0.007). However, no significant associations could be observed in the DR3-DQ2/DR1-DQ5 cohort.

Next-Generation Sequencing Identifies Extended HLA Class I and II Haplotypes Associated With Early-Onset and Late-Onset Myasthenia Gravis in Italian, Norwegian, and Swedish Populations

Genetic susceptibility to myasthenia gravis (MG) associates with specific HLA alleles and haplotypes at the class I and II regions in various populations. Previous studies have only examined alleles at a limited number of HLA loci that defined only broad serotypes or alleles defined at the protein sequence level. Consequently, genetic variants in noncoding and untranslated HLA gene segments have not been fully explored but could also be important determinants for MG. To gain further insight into the role of HLA in MG, we applied next-generation sequencing to analyze sequence variation at eleven HLA genes in early-onset (EO) and late-onset (LO) non-thymomatous MG patients positive for the acetylcholine receptor (AChR) antibodies and ethnically matched controls from Italy, Norway, and Sweden. For all three populations, alleles and haplotype blocks present on the ancestral haplotype AH8.1 were associated with risk in AChR-EOMG patients. HLA-B*08:01:01:01 was the dominant risk allele in Italians (OR = 3.28, P = 1.83E-05), Norwegians (OR = 3.52, P = 4.41E-16), and in Swedes HLA-B*08:01 was the primary risk allele (OR = 4.24, P <2.2E-16). Protective alleles and haplotype blocks were identified on the HLA-DRB7, and HLA-DRB13.1 class II haplotypes in Italians and Norwegians, whereas in Swedes HLA-DRB7 exhibited the main protective effect. For AChR-LOMG patients, the HLA-DRB15.1 haplotype and associated alleles were significantly associated with susceptibility in all groups. The HLA-DR13-HLA-DR-HLA-DQ haplotype was associated with protection in all AChR-LOMG groups. This study has confirmed and extended previous findings that the immunogenetic predisposition profiles for EOMG and LOMG are distinct. In addition, the results are consistent with a role for non-coding HLA genetic variants in the pathogenesis of MG.

High Resolution Haplotype Analyses of Classical HLA Genes in Families With Multiple Sclerosis Highlights the Role of HLA-DP Alleles in Disease Susceptibility

Multiple sclerosis (MS) susceptibility shows strong genetic associations with HLA alleles and haplotypes. We genotyped 11 HLA genes in 477 non-Hispanic European MS patients and their 954 unaffected parents using a validated next-generation sequencing (NGS) methodology. HLA haplotypes were assigned unequivocally by tracing HLA allele transmissions. We explored HLA haplotype/allele associations with MS using the genotypic transmission disequilibrium test (gTDT) and multiallelic TDT (mTDT). We also conducted a case-control (CC) study with all patients and 2029 healthy unrelated ethnically matched controls. We performed separate analyses of 54 extended multi-case families by reviewing transmission of haplotype blocks. The haplotype fragment including DRB5*01:01:01~DRB1*15:01:01:01 was significantly associated with predisposition (gTDT: p < 2.20e-16; mTDT: p =1.61e-07; CC: p < 2.22e-16) as reported previously. A second risk allele, DPB1*104:01 (gTDT: p = 3.69e-03; mTDT: p = 2.99e-03; CC: p = 1.00e-02), independent from the haplotype bearing DRB1*15:01 was newly identified. The allele DRB1*01:01:01 showed significant protection (gTDT: p = 8.68e-06; mTDT: p = 4.50e-03; CC: p = 1.96e-06). Two DQB1 alleles, DQB1*03:01 (gTDT: p = 2.86e-03; mTDT: p = 5.56e-02; CC: p = 4.08e-05) and DQB1*03:03 (gTDT: p = 1.17e-02; mTDT: p = 1.16e-02; CC: p = 1.21e-02), defined at two-field level also showed protective effects. The HLA class I block, A*02:01:01:01~C*03:04:01:01~B*40:01:02 (gTDT: p = 5.86e-03; mTDT: p = 3.65e-02; CC: p = 9.69e-03) and the alleles B*27:05 (gTDT: p = 6.28e-04; mTDT: p = 2.15e-03; CC: p = 1.47e-02) and B*38:01 (gTDT: p = 3.20e-03; mTDT: p = 6.14e-03; CC: p = 1.70e-02) showed moderately protective effects independently from each other and from the class II associated factors. By comparing statistical significance of 11 HLA loci and 19 haplotype segments with both untruncated and two-field allele names, we precisely mapped MS candidate alleles/haplotypes while eliminating false signals resulting from ‘hitchhiking’ alleles. We assessed genetic burden for the HLA allele/haplotype identified in this study. This family-based study including the highest-resolution of HLA alleles proved to be powerful and efficient for precise identification of HLA genotypes associated with both, susceptibility and protection to development of MS.

High-resolution HLA allele and haplotype frequencies in several unrelated populations determined by next generation sequencing: 17th International HLA and Immunogenetics Workshop joint report

The primary goal of the unrelated population HLA diversity (UPHD) component of the 17th International HLA and Immunogenetics Workshop was to characterize HLA alleles at maximum allelic-resolution in worldwide populations and re-evaluate patterns of HLA diversity across populations. The UPHD project included HLA genotype and sequence data, generated by various next-generation sequencing methods, from 4,240 individuals collated from 12 different countries. Population data included well-defined large datasets from the USA and smaller samples from Europe, Australia, and Western Asia. Allele and haplotype frequencies varied across populations from distant geographical regions. HLA genetic diversity estimated at 2- and 4-field allelic resolution revealed that diversity at the majority of loci, particularly for European-descent populations, was lower at the 2-field resolution. Several common alleles with identical protein sequences differing only by intronic substitutions were found in distinct haplotypes, revealing a more detailed characterization of linkage between variants within the HLA region. The examination of coding and non-coding nucleotide variation revealed many examples in which almost complete biunivocal relations between common alleles at different loci were observed resulting in higher linkage disequilibrium. Our reference data of HLA profiles characterized at maximum resolution from many populations is useful for anthropological studies, unrelated donor searches, transplantation, and disease association studies.

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-DQB1 mismatch increase risk of severe bleeding independently in recipients of allogeneic stem cell transplant

Severe bleeding is a major cause of death in acute leukemia (AL) patients with graft-versus-host disease (GVHD) after allogene hematopoietic stem-cell transplantation (allo-HSCT). However, the prognostic value and prediction of HSCT-associated severe bleeding in GVHD patients have not been reported in cohort studies. We did a retrospective analysis of 200 AL patients with GVHD after allo-HSCT from Feb 1, 2014, to Dec 1, 2015. Multivariate analysis showed that the severe bleeding class was associated with the risk of death (HR 2.26, 95% CI 1.31-3.92, p<0.001***). In order to predict severe bleeding and figure out the solution to bleeding events, we established a multiple logistic regression model. HLA-DQB1 unmatching, megakaryocyte reconsititution failure, and III or IV GVHD were the independent risk factors for severe bleeding. Among all the variations above, OR of HLA-DQB1 was the highest (OR: 16.02, 95% CI: 11.54-48.68). Adding HLA-DQB1 to other factors improved the reclassification for predicting severe bleeding (NRI=0.195, z=2.634, p=0.008**; IDI=0.289, z=3.249, p<0.001***). Lasso regression was used to select variants. A nomogram of the logistic model was generated and displayed. Calibration curve demonstrated excellent accuracy in estimating severe bleeding (C index of 0.935). HLA-DQB1 showed excellent efficacy of predicting severe bleeding in HSCT patients.

Standard reference sequences for submission of HLA genotyping for the 18th International HLA and Immunogenetics Workshop

The International human leukocyte antigen (HLA) and Immunogenetics Workshops (IHIWs) have fostered international collaborations of researchers and experts in the fields of HLA, histocompatibility and immunology. These IHIW collaborations have comprised many projects focused on achieving a variety of specific goals. The international and collaborative nature of these projects necessitates the collection and analysis of complex data generated in multiple laboratories, often using multiple methods of acquisition. Collection and storage of these data in a consistent way adds value to IHIW projects, which can be extended to future work. DNA‐based genotyping data, especially HLA genotyping data, can be transmitted in the form of a Histoimmunogenetics Markup Language (HML) document. HML facilitates clear communication of a genotype and supporting metadata, such as, sequencing platform, laboratory assays, consensus sequence, and interpretation. Sequence information can be reported relative to known reference sequences, which add meaning and context to genotypes. Selecting the correct reference sequence for a given allele sequence is nuanced, and guidelines have emerged through collaborative community efforts such as Data Standards Hackathons. Here, we describe the guidelines established for the selection of reference sequences to be used in transmission of HLA (and MICA/MICB) genotyping data for the 18th IHIW.

HLA diversity in the Argentinian Umbilical Cord Blood Bank: frequencies according to donor's reported ancestry and geographical distribution

Umbilical cord blood (UCB) is a suitable source for hematopoietic stem cell transplantation. The study of HLA genes by next generation sequencing is commonly used in transplants. Donor/patient HLA matching is often higher within groups of common ancestry, however "Hispanic" is a broad category that fails to represent Argentina's complex genetic admixture. Our aim is to describe HLA diversity of banked UCB units collected across the country taking into consideration donor's reported ancestral origins as well as geographic distribution. Our results showed an evenly distribution of units mainly for 2 groups: of European and of Native American descent, each associated to a defined geographic location pattern (Central vs. North regions). We observed differences in allele frequency distributions for some alleles previously described in Amerindian populations: for Class I (A*68:17, A*02:11:01G, A*02:22:01G, B*39:05:01, B*35:21, B*40:04, B*15:04:01G, B*35:04:01, B*51:13:01) and Class II (DRB1*04:11:01, DRB1*04:07:01G/03, DRB1*08:02:01, DRB1*08:07, DRB1*09:01:02G, DRB1*14:02:01, DRB1*16:02:01G). Our database expands the current knowledge of HLA diversity in Argentinian population. Although further studies are necessary to fully comprehend HLA heterogeneity, this report should prove useful to increase the possibility of finding compatible donors for successful allogeneic transplant and to improve recruitment strategies for UCB donors across the country.

Identification of the novel HLA-A*02:837 and -A*02:888 alleles by next-generation sequencing in two Chinese individuals

Compared with HLA-A*02:06:01:01, HLA-A*02:837 and HLA-A*02:888 show one single nucleotide substitution respectively.

T cell analysis in vaccination

A major problem in the analysis of vaccine candidates is the lack of any agreed upon surrogates of efficacy, which means that for diseases that depend on a strong T cell response (HIV, TB especially) the only option is to perform an efficacy trial, involving thousands of subjects, enormous costs, and years before the results are known [1]. We also know that T cell responses are an important part of most pathogen responses, and so identifying key T cell response metrics in early vaccine trials would be generally useful. Given our ignorance of what the most important variables are, what would we like to measure and how can this be accomplished, especially given the explosion of new technologies that are available? What follows is a consideration of what should be measured, with the caveat that some of these will be more important than others.

Genomic variations in EBNA3C of EBV associate with posttransplant lymphoproliferative disorder

Epstein-Barr Virus (EBV) is a ubiquitous virus linked to a variety of lymphoid and epithelial malignancies. In solid organ and hematopoietic stem cell transplant recipients, EBV is causally associated with posttransplant lymphoproliferative disorder (PTLD), a group of heterogeneous lymphoid diseases. EBV+ B cell lymphomas that develop in the context of PTLD are generally attributed to the immunosuppression required to promote graft survival, but little is known regarding the role of EBV genome diversity in the development of malignancy. We deep-sequenced the EBV genome from the peripheral blood of 18 solid organ transplant recipients, including 6 PTLD patients. Sequences from 6 EBV+ spontaneous lymphoblastoid B cell lines (SLCL) were similarly analyzed. The EBV genome from PTLD patients had a significantly greater number of variations than EBV from transplant recipients without PTLD. Importantly, there were 15 nonsynonymous variations, including 8 in the latent cycle gene EBNA3C that were associated with the development of PTLD. One of the nonsynonymous variations in EBNA3C is located within a previously defined T cell epitope. These findings suggest that variations in the EBV genome can contribute to the pathogenesis of PTLD.

Complete nucleotide sequence characterization of DRB5 alleles reveals a homogeneous allele group that is distinct from other DRB genes

Next Generation Sequencing allows for testing and typing of entire genes of the HLA region. A better and comprehensive sequence assessment can be achieved by the inclusion of full gene sequences of all the common alleles at a given locus. The common alleles of DRB5 are under-characterized with the full exon-intron sequence of two alleles available. In the present study the DRB5 genes from 18 subjects alleles were cloned and sequenced; haplotype analysis showed that 17 of them had a single copy of DRB5 and one consanguineous subject was homozygous at all HLA loci. Methodological approaches including robust and efficient long-range PCR amplification, molecular cloning, nucleotide sequencing and de novo sequence assembly were combined to characterize DRB5 alleles. DRB5 sequences covering from 5'UTR to the end of intron 5 were obtained for DRB5*01:01, 01:02 and 02:02; partial coverage including a segment spanning exon 2 to exon 6 was obtained for DRB5*01:03, 01:08N and 02:03. Phylogenetic analysis of the generated sequences showed that the DRB5 alleles group together and have distinctive differences with other DRB loci. Novel intron variants of DRB5*01:01:01, 01:02 and 02:02 were identified. The newly characterized DRB5 intron variants of each DRB5 allele were found in subjects harboring distinct associations with alleles of DRB1, B and/or ethnicity. The new information provided by this study provides reference sequences for HLA typing methodologies. Extending sequence coverage may lead to identify the disease susceptibility factors of DRB5 containing haplotypes while the unexpected intron variations may shed light on understanding of the evolution of the DRB region.