Structural and regulatory diversity shape HLA-C protein expression levels

Minor histocompatibility antigens to predict, monitor or manipulate GvL and GvHD after allogeneic hematopoietic cell transplantation

Allogeneic hematopoietic cell transplantation (alloHCT) provides a potential curative treatment for haematological malignancies. The therapeutic Graft-versus-Leukaemia (GvL) effect is induced by donor T cells attacking patient hematopoietic (malignant) cells. However, if healthy non-hematopoietic tissues are targeted, Graft-versus-Disease (GvHD) may develop. After HLA-matched alloHCT, GvL and GvHD are induced by donor T cells recognizing polymorphic peptides presented by HLA on patient cells, so-called minor histocompatibility antigens (MiHAs). The balance between GvL and GvHD depends on the tissue distribution of MiHAs and T-cell frequencies targeting these MiHAs. T cells against broadly expressed MiHAs induce GvL and GvHD, whereas those targeting MiHAs with hematopoietic-restricted expression induce GvL without GvHD. Recently, the MiHA repertoire identified in natural immune responses after alloHCT was expanded to 159 total HLA-I-restricted MiHAs, including 14 hematopoietic-restricted MiHAs. This review explores their potential relevance to predict, monitor, and manipulate GvL and GvHD for improving clinical outcome after HLA-matched alloHCT.

The most frequent HLA alleles around the world: A fundamental synopsis

A comprehensive knowledge of human leukocyte antigen (HLA) molecular variation worldwide is essential in human population genetics research and disease association studies and is also indispensable for clinical applications such as allogeneic hematopoietic cell transplantation, where ensuring HLA compatibility between donors and recipients is paramount. Enormous progress has been made in this field thanks to several decades of HLA population studies allowing the development of helpful databases and bioinformatics tools. However, it is still difficult to appraise the global HLA population diversity in a synthetic way. We thus introduce here a novel approach, based on approximately 2000 data sets, to assess this complexity by providing a fundamental synopsis of the most frequent HLA alleles observed in different regions of the world. This new knowledge will be useful not only as a fundamental reference for basic research, but also as an efficient guide for clinicians working in the field of transplantation.

Stimulating T cell responses against patient-derived breast cancer cells with neoantigen peptide-loaded peripheral blood mononuclear cells

Breast cancer stands as a formidable global health challenge for women. While neoantigens exhibit efficacy in activating T cells specific to cancer and instigating anti-tumor immune responses, the accuracy of neoantigen prediction remains suboptimal. In this study, we identified neoantigens from the patient-derived breast cancer cells, PC-B-142CA and PC-B-148CA cells, utilizing whole-genome and RNA sequencing. The pVAC-Seq pipeline was employed, with minor modification incorporating criteria (1) binding affinity of mutant (MT) peptide with HLA (IC50 MT) ≤ 500 nm in 3 of 5 algorithms and (2) IC50 wild type (WT)/MT > 1. Sequencing results unveiled 2513 and 3490 somatic mutations, and 646 and 652 non-synonymous mutations in PC-B-142CA and PC-B-148CA, respectively. We selected the top 3 neoantigens to perform molecular dynamic simulation and synthesized 9-12 amino acid neoantigen peptides, which were then pulsed onto healthy donor peripheral blood mononuclear cells (PBMCs). Results demonstrated that T cells activated by ADGRL1E274K, PARP1E619K, and SEC14L2R43Q peptides identified from PC-B-142CA exhibited significantly increased production of interferon-gamma (IFN-γ), while PARP1E619K and SEC14L2R43Q peptides induced the expression of CD107a on T cells. The % tumor cell lysis was notably enhanced by T cells activated with MT peptides across all three healthy donors. Moreover, ALKBH6V83M and GAAI823T peptides from PC-B-148CA remarkably stimulated IFN-γ- and CD107a-positive T cells, displaying high cell-killing activity against target cancer cells. In summary, our findings underscore the successful identification of neoantigens with anti-tumor T cell functions and highlight the potential of personalized neoantigens as a promising avenue for breast cancer treatment.

Impact of HLA Class I Antigen, Killer Inhibitory Receptor, and FCGR3A Genotypes on Breast Cancer Susceptibility and Tumor Stage

BACKGROUND

The identification in breast cancer (BC) of novel genetic biomarkers regulating natural killer (NK) cell function, including the HLA, KIR, and CD16A (FCGR3A), may be still a challenge.

OBJECTIVE

We aimed to evaluate whether the combined effect of these polymorphisms has an impact on BC susceptibility and progression.

METHODS

47 BC Italian patients and healthy individuals (39 females and 66 males/ females) were genotyped by Sanger sequencing (HLA-C exon 2-4 and FCGR3A- 158V/F, 48L/R/H) and PCR-SSP typing (KIR genes).

RESULTS

HLA-C gene allele analysis showed the group C1, with HLA-C*07:02:01 allele, to be significantly associated with tumor progression (16.7% vs. 4.0%, p=0.04, OR=4.867), and instead, group C2, with HLA-C*05:01:01, was protective against disease susceptibility (0.0% vs. 7.2%, p=0.019, OR=0.087). In addition, we highlighted a significant reduction of the KIR2DS4ins in BC patients (pcorr.=0.022) and an increased combined presence of KIR2DL1 and KIR2DS1 genes in advanced BC patients compared to earlier stages (66.7% vs. 19.2%, p=0.002). The concurrent lack of KIR2DL2 and KIR2DS4 genes in the presence of HLA-C2 alleles was significantly associated with increased susceptibility to BC (p=0.012, OR=5.020) or with lymph node involvement (p=0.008, OR=6.375). Lastly, we identified different combinations of the FCGR3A-48/158 variants and KIR genes in BC patients compared to controls.

CONCLUSION

Our findings suggest that in the development of BC probably exists a disorder of the NK innate immunity influenced by KIR/HLA-C gene content and FCGR3A-158 polymorphisms and that the combined analysis of these biomarkers might help predict genetic risk scores for tailored screening of BC patients in therapy.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Proteogenomics Reveal the Overexpression of HLA-I in Cancer

An accurate quantification of HLA class I gene expression is important in understanding the interplay with the tumor microenvironment of antitumor cytotoxic T cell activities. Because HLA-I sequences are highly variable, standard RNAseq and mass spectrometry-based quantification workflows using common genome and protein sequence references do not provide HLA-I allele specific quantifications. Here, we used personalized HLA-I nucleotide and protein reference sequences based on the subjects' HLA-I genotypes and surveyed tumor and adjacent normal samples from patients across nine cancer types. Mass spectrometry using data dependent acquisition data was validated to be sufficient to estimate HLA-A protein expression at the allele level. We found that HLA-I proteins were present in significantly higher levels in tumors compared to adjacent normal tissues from 41 to 63% of head and neck squamous cell carcinoma, uterine corpus endometrial carcinoma, and clear cell renal cell carcinoma patients, and this was driven by increased levels of HLA-I gene transcripts. Most immune cell types are universally enriched in HLA-I high tumors, while endothelial and neuronal cells showed divergent relationships with HLA-I. Pathway analysis revealed that tumor senescence and autophagy activity influence the level of HLA-I proteins in glioblastoma. Genes correlated to HLA-I protein expression are mostly the ones directly involved in HLA-I function in immune response and cell death, while glycosylation genes are exclusively co-expressed with HLA-I at the protein level.

The complex HLA-E-nonapeptide in Behçet disease

Introduction

The knowledge of the aetiology of Behçet disease (BD), an immune-mediated vasculitis, is limited. HLA-B, mainly HLA-B51, and HLA-A molecules are associated with disease, but the ultimate cause of this association remains obscure. There is evidence that NK cells participate in the etiopathology of BD. NK cells have activator and inhibitor surface receptors, like the KIR and the NKG2 families. Classical HLA-class I molecules (A, B and C) are keys in the activity control of the NK because they are KIR ligands. Most NKG2 receptors bind HLA-E, which presents only nonapeptides derived from the signal peptide of other class-I molecules.

Objective

This study investigates the contribution of the pair HLA-E and ligand, nonapeptide derived from the 3-11 sequence of the signal peptides of class I classical molecules, to the susceptibility to BD.

Methods

We analyzed the frequency of the HLA-derivated nonapeptide forms in 466 BD patients and 444 controls and an HLA-E functional dimorphism in a subgroup of patients and controls. Results: In B51 negative patients, the frequency of VMAPRTLLL was lower (70.4% versus 80.0% in controls; P=0.006, Pc=0.04, OR=0.60, 95%CI 0.41-0.86), and the frequency of VMAPRTLVL was higher (81.6% versus 71.4% in controls; P=0.004, Pc=0.03, OR=1.78, 95%CI 1.20-2.63). In homozygosity, VMAPRTLLL is protective, and VMAPRTLVL confers risk. The heterozygous condition is neutral. There were no significant differences in the distribution of the HLA-E dimorphism.

Discussion

Our results explain the association of BD with diverse HLA-A molecules, reinforce the hypothesis of the involvement of the NK cells in the disease and do not suggest a significant contribution of the HLA-E polymorphism to disease susceptibility.

"HLA-C: evolution, epigenetics, and pathological implications in the major histocompatibility complex"

HLA-C, a gene located within the major histocompatibility complex, has emerged as a prominent target in biomedical research due to its involvement in various diseases, including cancer and autoimmune disorders; even though its recent addition to the MHC, the interaction between HLA-C and KIR is crucial for immune responses, particularly in viral infections. This review provides an overview of the structure, origin, function, and pathological implications of HLA-C in the major histocompatibility complex. In the last decade, we systematically reviewed original publications from Pubmed, ScienceDirect, Scopus, and Google Scholar. Our findings reveal that genetic variations in HLA-C can determine susceptibility or resistance to certain diseases. However, the first four exons of HLA-C are particularly susceptible to epigenetic modifications, which can lead to gene silencing and alterations in immune function. These alterations can manifest in diseases such as alopecia areata and psoriasis and can also impact susceptibility to cancer and the effectiveness of cancer treatments. By comprehending the intricate interplay between genetic and epigenetic factors that regulate HLA-C expression, researchers may develop novel strategies for preventing and treating diseases associated with HLA-C dysregulation.

Precision Neoantigen Discovery Using Large-Scale Immunopeptidomes and Composite Modeling of MHC Peptide Presentation

Major histocompatibility complex (MHC)-bound peptides that originate from tumor-specific genetic alterations, known as neoantigens, are an important class of anticancer therapeutic targets. Accurately predicting peptide presentation by MHC complexes is a key aspect of discovering therapeutically relevant neoantigens. Technological improvements in mass spectrometry-based immunopeptidomics and advanced modeling techniques have vastly improved MHC presentation prediction over the past 2 decades. However, improvement in the accuracy of prediction algorithms is needed for clinical applications like the development of personalized cancer vaccines, the discovery of biomarkers for response to immunotherapies, and the quantification of autoimmune risk in gene therapies. Toward this end, we generated allele-specific immunopeptidomics data using 25 monoallelic cell lines and created Systematic Human Leukocyte Antigen (HLA) Epitope Ranking Pan Algorithm (SHERPA), a pan-allelic MHC-peptide algorithm for predicting MHC-peptide binding and presentation. In contrast to previously published large-scale monoallelic data, we used an HLA-null K562 parental cell line and a stable transfection of HLA allele to better emulate native presentation. Our dataset includes five previously unprofiled alleles that expand MHC diversity in the training data and extend allelic coverage in underprofiled populations. To improve generalizability, SHERPA systematically integrates 128 monoallelic and 384 multiallelic samples with publicly available immunoproteomics data and binding assay data. Using this dataset, we developed two features that empirically estimate the propensities of genes and specific regions within gene bodies to engender immunopeptides to represent antigen processing. Using a composite model constructed with gradient boosting decision trees, multiallelic deconvolution, and 2.15 million peptides encompassing 167 alleles, we achieved a 1.44-fold improvement of positive predictive value compared with existing tools when evaluated on independent monoallelic datasets and a 1.17-fold improvement when evaluating on tumor samples. With a high degree of accuracy, SHERPA has the potential to enable precision neoantigen discovery for future clinical applications.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation, and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here, we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues, using personalized reference genomes to mitigate technical confounding. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B, and T cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.

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.

Association between HLA-C alleles and COVID-19 severity in a pilot study with a Spanish Mediterranean Caucasian cohort

The clinical presentations of COVID-19 may range from an asymptomatic or mild infection to a critical or fatal disease. Several host factors such as elderly age, male gender, and previous comorbidities seem to be involved in the most severe outcomes, but also an impaired immune response that causes a hyperinflammatory state but is unable to clear the infection. In order to get further understanding about this impaired immune response, we aimed to determine the association of specific HLA alleles with different clinical presentations of COVID-19. Therefore, we analyzed HLA Class I and II, as well as KIR gene sequences, in 72 individuals with Spanish Mediterranean Caucasian ethnicity who presented mild, severe, or critical COVID-19, according to their clinical characteristics and management. This cohort was recruited in Madrid (Spain) during the first and second pandemic waves between April and October 2020. There were no significant differences in HLA-A or HLA-B alleles among groups. However, despite the small sample size, we found that HLA-C alleles from group C1 HLA-C*08:02, -C*12:03, or -C*16:01 were more frequently associated in individuals with mild COVID-19 (43.8%) than in individuals with severe (8.3%; p = 0.0030; pc = 0.033) and critical (16.1%; p = 0.0014; pc = 0.0154) disease. C1 alleles are supposed to be highly efficient to present peptides to T cells, and HLA-C*12:03 may present a high number of verified epitopes from abundant SARS-CoV-2 proteins M, N, and S, thereby being allegedly able to trigger an efficient antiviral response. On the contrary, C2 alleles are usually poorly expressed on the cell surface due to low association with β2-microglobulin (β2M) and peptides, which may impede the adequate formation of stable HLA-C/β2M/peptide heterotrimers. Consequently, this pilot study described significant differences in the presence of specific HLA-C1 alleles in individuals with different clinical presentations of COVID-19, thereby suggesting that HLA haplotyping could be valuable to get further understanding in the underlying mechanisms of the impaired immune response during critical COVID-19.

Mouse fetal growth restriction through parental and fetal immune gene variation and intercellular communications cascade

Fetal growth restriction (FGR) affects 5-10% of pregnancies, and can have serious consequences for both mother and child. Prevention and treatment are limited because FGR pathogenesis is poorly understood. Genetic studies implicate KIR and HLA genes in FGR, however, linkage disequilibrium, genetic influence from both parents, and challenges with investigating human pregnancies make the risk alleles and their functional effects difficult to map. Here, we demonstrate that the interaction between the maternal KIR2DL1, expressed on uterine natural killer (NK) cells, and the paternally inherited HLA-C*0501, expressed on fetal trophoblast cells, leads to FGR in a humanized mouse model. We show that the KIR2DL1 and C*0501 interaction leads to pathogenic uterine arterial remodeling and modulation of uterine NK cell function. This initial effect cascades to altered transcriptional expression and intercellular communication at the maternal-fetal interface. These findings provide mechanistic insight into specific FGR risk alleles, and provide avenues of prevention and treatment.

Host KIR/HLA-C Genotypes Determine HIV-Mediated Changes of the NK Cell Repertoire and Are Associated With Vpu Sequence Variations Impacting Downmodulation of HLA-C

NK cells play a pivotal role in viral immunity, utilizing a large array of activating and inhibitory receptors to identify and eliminate virus-infected cells. Killer-cell immunoglobulin-like receptors (KIRs) represent a highly polymorphic receptor family, regulating NK cell activity and determining the ability to recognize target cells. Human leukocyte antigen (HLA) class I molecules serve as the primary ligand for KIRs. Herein, HLA-C stands out as being the dominant ligand for the majority of KIRs. Accumulating evidence indicated that interactions between HLA-C and its inhibitory KIR2DL receptors (KIR2DL1/L2/L3) can drive HIV-1-mediated immune evasion and thus may contribute to the intrinsic control of HIV-1 infection. Of particular interest in this context is the recent observation that HIV-1 is able to adapt to host HLA-C genotypes through Vpu-mediated downmodulation of HLA-C. However, our understanding of the complex interplay between KIR/HLA immunogenetics, NK cell-mediated immune pressure and HIV-1 immune escape is still limited. Therefore, we investigated the impact of specific KIR/HLA-C combinations on the NK cell receptor repertoire and HIV-1 Vpu protein sequence variations of 122 viremic, untreated HIV-1⁺ individuals. Compared to 60 HIV-1^- controls, HIV-1 infection was associated with significant changes within the NK cell receptor repertoire, including reduced percentages of NK cells expressing NKG2A, CD8, and KIR2DS4. In contrast, the NKG2C⁺ and KIR3DL2⁺ NK cell sub-populations from HIV-1⁺ individuals was enlarged compared to HIV-1^- controls. Stratification along KIR/HLA-C genotypes revealed a genotype-dependent expansion of KIR2DL1⁺ NK cells that was ultimately associated with increased binding affinities between KIR2DL1 and HLA-C allotypes. Lastly, our data hinted to a preferential selection of Vpu sequence variants that were associated with HLA-C downmodulation in individuals with high KIR2DL/HLA-C binding affinities. Altogether, our study provides evidence that HIV-1-associated changes in the KIR repertoire of NK cells are to some extent predetermined by host KIR2DL/HLA-C genotypes. Furthermore, analysis of Vpu sequence polymorphisms indicates that differential KIR2DL/HLA-C binding affinities may serve as an additional mechanism how host genetics impact immune evasion by HIV-1.

Fine-mapping studies distinguish genetic risks for childhood- and adult-onset asthma in the HLA region

BACKGROUND

Genome-wide association studies of asthma have revealed robust associations with variation across the human leukocyte antigen (HLA) complex with independent associations in the HLA class I and class II regions for both childhood-onset asthma (COA) and adult-onset asthma (AOA). However, the specific variants and genes contributing to risk are unknown.

METHODS

We used Bayesian approaches to perform genetic fine-mapping for COA and AOA (n=9432 and 21,556, respectively; n=318,167 shared controls) in White British individuals from the UK Biobank and to perform expression quantitative trait locus (eQTL) fine-mapping in immune (lymphoblastoid cell lines, n=398; peripheral blood mononuclear cells, n=132) and airway (nasal epithelial cells, n=188) cells from ethnically diverse individuals. We also examined putatively causal protein coding variation from protein crystal structures and conducted replication studies in independent multi-ethnic cohorts from the UK Biobank (COA n=1686; AOA n=3666; controls n=56,063).

RESULTS

Genetic fine-mapping revealed both shared and distinct causal variation between COA and AOA in the class I region but only distinct causal variation in the class II region. Both gene expression levels and amino acid variation contributed to risk. Our results from eQTL fine-mapping and amino acid visualization suggested that the HLA-DQA1*03:01 allele and variation associated with expression of the nonclassical HLA-DQA2 and HLA-DQB2 genes accounted entirely for the most significant association with AOA in GWAS. Our studies also suggested a potentially prominent role for HLA-C protein coding variation in the class I region in COA. We replicated putatively causal variant associations in a multi-ethnic cohort.

CONCLUSIONS

We highlight roles for both gene expression and protein coding variation in asthma risk and identified putatively causal variation and genes in the HLA region. A convergence of genomic, transcriptional, and protein coding evidence implicates the HLA-DQA2 and HLA-DQB2 genes and HLA-DQA1*03:01 allele in AOA.

Recent Advances in Influenza, HIV and SARS-CoV-2 Infection Prevention and Drug Treatment—The Need for Precision Medicine

Viruses, and in particular, RNA viruses, dominate the WHO’s current list of ten global health threats. Of these, we review the widespread and most common HIV, influenza virus, and SARS-CoV-2 infections, as well as their possible prevention by vaccination and treatments by pharmacotherapeutic approaches. Beyond the vaccination, we discuss the virus-targeting and host-targeting drugs approved in the last five years, in the case of SARS-CoV-2 in the last one year, as well as new drug candidates and lead molecules that have been published in the same periods. We share our views on vaccination and pharmacotherapy, their mutually reinforcing strategic significance in combating pandemics, and the pros and cons of host and virus-targeted drug therapy. The COVID-19 pandemic has provided evidence of our limited armamentarium to fight emerging viral diseases. Novel broad-spectrum vaccines as well as drugs that could even be applied as prophylactic treatments or in early phases of the viremia, possibly through oral administration, are needed in all three areas. To meet these needs, the use of multi-data-based precision medicine in the practice and innovation of vaccination and drug therapy is inevitable.

Late antibody-mediated rejection in a kidney transplant recipient: COVID 19 induced?

Background

Antibody-mediated rejection (AMR) was described in kidney transplant patients after viral infections, such as the cytomegalovirus. Very few cases were recently reported after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, probably in the context of lowering of immunosuppressive therapy. To date, no direct immunological link was proved to explain a connection between the coronavirus disease 19 (COVID-19) infection and antibody-mediated rejection (AMR) if it exists.

Case presentation

Here we try to find this association by presenting the case of a low immunological risk patient who presented, six years post-transplant, with c4d negative antibody mediated rejection due to an anti-HLA-C17 de novo donor specific antibody (DSA) after contracting the coronavirus disease 19.

The HLA-Cw17 activated the antibody-dependent cell-mediated cytotoxicity via the KIR2DS1 positive NK cells.

Discussion and conclusions

This case report may prove a direct role for COVID-19 infection in AMRs in the kidney transplant recipients, leading us to closely monitor kidney transplant recipients, especially if they have “at-risk” donor antigens.

Principles of Virtual Crossmatch Testing for Kidney Transplantation

Human leukocyte antigens (HLAs) are the primary determinants of alloimmunity. A crossmatch test is a test that determines the immunologic risk of a recipient with a potential donor by ensuring that there are no transplant-relevant circulating antibodies in the recipient directed against donor antigens. Physical crossmatch (PXM) tests, such as complement-dependent cytotoxicity crossmatch (CDCXM) and flow cytometry crossmatch (FCXM), require mixing of patient serum and donor cells, are labor intensive, and are logistically challenging. Virtual crossmatch (VXM) test assesses immunologic compatibility between recipient and potential donor by analyzing the results of 2 independently done physical laboratory tests—patient anti-HLA antibody and donor HLA typing. The goal of VXM is pretransplant risk stratification—though there is no consensus on whether such risk assessment involves predicting the PXM result or the posttransplant outcome. Although the concept of VXM is not new, the advent of solid-phase assays for detecting circulating antibodies in the recipient directed against individual HLA and DNA-based methods for typing donor HLA specificities at a higher resolution makes the routine use of VXM a reality. Accordingly, VXM may be applied at different scenarios—both for sensitized and nonsensitized patients. Implementation of VXM-based approach has resulted in statistically significant reduction in cold ischemia time without an increase in hyperacute rejection episodes. Though there are considerable challenges, VXM is expected to be used more often in the future, depending on the transplant center’s tolerance of immunologic risk.

High cell surface expression and peptide binding affinity of HLA-DQA1*05:03, a susceptible allele of neuromyelitis optica spectrum disorders (NMOSD)

Neuromyelitis optica spectrum disorder (NMOSD) is a relapsing autoimmune disease characterized by the presence of pathogenic autoantibodies, anti-aquaporin 4 (AQP4) antibodies. Recently, HLA-DQA1*05:03 was shown to be significantly associated with NMOSD in a Japanese patient cohort. However, the specific mechanism by which HLA-DQA1*05:03 is associated with the development of NMOSD has yet to be elucidated. In the current study, we revealed that HLA-DQA1*05:03 exhibited significantly higher cell surface expression levels compared to other various DQA1 alleles, and that its expression strongly depended on the amino acid sequence of the α1 domain, with a preference for leucine at position 75. Moreover, in silico analysis indicated that the HLA-DQ encoded by HLA-DQA1*05:03 preferentially presents immunodominant AQP4 peptides, and that the peptide major histocompatibility complexes (pMHCs) are more energetically stable in the presence of HLA-DQA1*05:03 than other HLA-DQA1 alleles. In silico 3D structural models were also applied to investigate the validity of the energetic stability of pMHCs. Taken together, our findings indicate that HLA-DQA1*05:03 possesses a distinct property to play a pathogenic role in the development of NMOSD.

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 New Kid on the Block: HLA-C, a Key Regulator of Natural Killer Cells in Viral Immunity

The human leukocyte antigen system (HLA) is a cluster of highly polymorphic genes essential for the proper function of the immune system, and it has been associated with a wide range of diseases. HLA class I molecules present intracellular host- and pathogen-derived peptides to effector cells of the immune system, inducing immune tolerance in healthy conditions or triggering effective immune responses in pathological situations. HLA-C is the most recently evolved HLA class I molecule, only present in humans and great apes. Differentiating from its older siblings, HLA-A and HLA-B, HLA-C exhibits distinctive features in its expression and interaction partners. HLA-C serves as a natural ligand for multiple members of the killer-cell immunoglobulin-like receptor (KIR) family, which are predominately expressed by natural killer (NK) cells. NK cells are crucial for the early control of viral infections and accumulating evidence indicates that interactions between HLA-C and its respective KIR receptors determine the outcome and progression of viral infections. In this review, we focus on the unique role of HLA-C in regulating NK cell functions and its consequences in the setting of viral infections.

21-Hydroxylase-Specific CD8+ T Cells in Autoimmune Addison's Disease Are Restricted by HLA-A2 and HLA-C7 Molecules

Objectives

CD8+ T cells targeting 21-hydroxylase (21OH) are presumed to play a central role in the destruction of adrenocortical cells in autoimmune Addison’s disease (AAD). Earlier reports have suggested two immunodominant CD8+ T cell epitopes within 21OH: LLNATIAEV (21OH_342-350), restricted by HLA-A2, and EPLARLEL (21OH_431-438), restricted by HLA-B8. We aimed to characterize polyclonal CD8+ T cell responses to the proposed epitopes in a larger patient cohort with AAD.

Methods

Recombinant fluorescent HLA-peptide multimer reagents were used to quantify antigen-specific CD8+ T cells by flow cytometry. Interferon-gamma (IFNγ) Elispot and biochemical assays were used to functionally investigate the 21OH-specific T cells, and to map the exactly defined epitopes of 21OH.

Results

We found a significantly higher frequency of HLA-A2 restricted LLNATIAEV-specific cells in patients with AAD than in controls. These cells could also be expanded in vitro in an antigen specific manner and displayed a robust antigen-specific IFNγ production. In contrast, only negligible frequencies of EPLARLEL-specific T cells were detected in both patients and controls with limited IFNγ response. However, significant IFNγ production was observed in response to a longer peptide encompassing EPLARLEL, 21OH_430-447, suggesting alternative dominant epitopes. Accordingly, we discovered that the slightly offset ARLELFVVL (21OH_434-442) peptide is a novel dominant epitope restricted by HLA-C7 and not by HLA-B8 as initially postulated.

Conclusion

We have identified two dominant 21OH epitopes targeted by CD8+ T cells in AAD, restricted by HLA-A2 and HLA-C7, respectively. To our knowledge, this is the first HLA-C7 restricted epitope described for an autoimmune disease.

ERAP1 Controls the Autoimmune Response against Melanocytes in Psoriasis by Generating the Melanocyte Autoantigen and Regulating Its Amount for HLA-C*06:02 Presentation

Autoimmune diseases develop when autoantigens activate previously quiescent self-reactive lymphocytes. Gene-gene interaction between certain HLA class I risk alleles and variants of the endoplasmic reticulum aminopeptidase ERAP1 controls the risk for common immune-mediated diseases, including psoriasis, ankylosing spondylitis, and Behçet disease. The functional mechanisms underlying this statistical association are unknown. In psoriasis, HLA-C*06:02 mediates an autoimmune response against melanocytes by autoantigen presentation. Using various genetically modified cell lines together with an autoreactive psoriatic TCR in a TCR activation assay, we demonstrate in this study that in psoriasis, ERAP1 generates the causative melanocyte autoantigen through trimming N-terminal elongated peptide precursors to the appropriate length for presentation by HLA-C*06:02. An ERAP1 risk haplotype for psoriasis produced the autoantigen much more efficiently and increased HLA-C expression and stimulation of the psoriatic TCR by melanocytes significantly more than a protective haplotype. Compared with the overall HLA class I molecules, cell surface expression of HLA-C decreased significantly more upon ERAP1 knockout. The combined upregulation of ERAP1 and HLA-C on melanocytes in psoriasis lesions emphasizes the pathogenic relevance of their interaction in patients. We conclude that in psoriasis pathogenesis, the increased generation of an ERAP1-dependent autoantigen by an ERAP1 risk haplotype enhances the likelihood that autoantigen presentation by HLA-C*06:02 will exceed the threshold for activation of potentially autoreactive T cells, thereby triggering CD8+ T cell-mediated autoimmune disease. These data identify ERAP1 function as a central checkpoint and promising therapeutic target in psoriasis and possibly other HLA class I-associated diseases with a similar genetic predisposition.

Inside the pocket: Critical elements of HLA-mediated susceptibility to cervical precancerous lesions

Human papillomavirus (HPV) infection is a necessary cause for cervical cancer (CC), but it also depends on genetic factors, such as HLA polymorphism. However, few reports addressed the role of amino acids residues at the HLA peptide-binding cleft in HPV-related cervical disease. Therefore, we aimed to investigate the association between HLA-B, HLA-C, and HLA-DRB1 polymorphism and amino acid residues composing the pockets of the peptide-binding cleft of the respective polypeptide chains with cervical intraepithelial neoplasia (CIN II/III). HLA typing was performed by PCR-SSOP in 184 women with CIN II/III and 174 controls from South Brazil. Associations were estimated by multivariate logistic regression. FDR test was performed to correct the p-value for multiple comparisons. HLA-DRB1*13:01 was associated with protection against CIN II/III, while HLA-C*03:04 was associated with susceptibility. The amino acid residues isoleucine, tyrosine, and leucine at positions 95, 116, and 163 of HLA-C, respectively, were associated with CIN II/III susceptibility. In contrast, serine at positions 11 and 13 of HLA-DRB1 was associated with protection against the disease. Our results confirm previously reported associations between HLA and cervical diseases caused by HPV and suggest a role for amino acid residues at different positions of HLA-C and HLA-DRB1 in CIN II/III. This finding may be further explored to better understand the genetic risk and the influence of immune response to CC development.

Negative trade-off between neoantigen repertoire breadth and the specificity of HLA-I molecules shapes antitumor immunity

Human leukocyte antigen class I (HLA-I) genes shape our immune response against pathogens and cancer. Certain HLA-I variants can bind a wider range of peptides than others, a feature that could be favorable against a range of viral diseases. However, the implications of this phenomenon on cancer immune response are unknown. Here we quantified peptide repertoire breadth (or promiscuity) of a representative set of HLA-I alleles and found that patients with cancer who were carrying HLA-I alleles with high peptide-binding promiscuity have significantly worse prognosis after immune checkpoint inhibition. This can be explained by a reduced capacity of the immune system to discriminate tumor neopeptides from self-peptides when patients carry highly promiscuous HLA-I variants, shifting the regulation of tumor-infiltrating T cells from activation to tolerance. In summary, HLA-I peptide-binding specificity shapes neopeptide immunogenicity and the self-immunopeptidome repertoire in an antagonistic manner, and could underlie a negative trade-off between antitumor immunity and genetic susceptibility to viral infections.

An immunogenetic view of COVID-19

Meeting the challenges brought by the COVID-19 pandemic requires an interdisciplinary approach. In this context, integrating knowledge of immune function with an understanding of how genetic variation influences the nature of immunity is a key challenge. Immunogenetics can help explain the heterogeneity of susceptibility and protection to the viral infection and disease progression. Here, we review the knowledge developed so far, discussing fundamental genes for triggering the innate and adaptive immune responses associated with a viral infection, especially with the SARS-CoV-2 mechanisms. We emphasize the role of the HLA and KIR genes, discussing what has been uncovered about their role in COVID-19 and addressing methodological challenges of studying these genes. Finally, we comment on questions that arise when studying admixed populations, highlighting the case of Brazil. We argue that the interplay between immunology and an understanding of genetic associations can provide an important contribution to our knowledge of COVID-19.

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.

Withdrawn: Precision Neoantigen Discovery Using Large-scale Immunopeptidomes and Composite Modeling of MHC Peptide Presentation

This article has been withdrawn by the authors. A publication of the manuscript with the correct figures and tables has been approved and the authors state the conclusions of the manuscript remain unaffected. Specifically, errors are in Figure 6A, Supplementary Figure 10B, Supplementary Figure 10C, and Supplementary Table 5. The details of the errors are as follows: the HLA types for one sample were incorrectly assigned because of a tumor/normal mislabeling from the biobank vendor. Due to the differing HLA types between the tumor and normal sample, the sequence analysis established that the HLA alleles for this patient had been deleted (HLA LOH). The authors conclude that this was an artifact caused by the normal sample mislabeling. The corrected version can be accessed (Pyke, R.M., Mellacheruvu, D., Dea, S., Abbott, C.W., Zhang, S.V., Philips, N.A., Harris, J., Bartha, G., Desai, S., McClory, R., West, J., Snyder, M,P., Chen, R., Boyle, S.M. (2022) Precision Neoantigen Discovery Using Large-Scale Immunopeptidomics and Composite Modeling of MHC Peptide Presentation. Mol. Cell. Proteomics 22, 100506

HLA Class II Presentation Is Specifically Altered at Elevated Temperatures in the B-Lymphoblastic Cell Line JY

Human leukocyte antigen (HLA) molecules play critical roles in our adaptive immune system by signaling a cell's health status to the immune system, through presentation of small peptides. Understanding HLA biology is important because of its prominent role in autoimmune diseases and cancer immunotherapy. Although both the HLA class I and class II antigen processing and presentation pathways have been studied extensively, the fundamental rules in HLA class II antigen presentation still remain less understood. To clarify the mechanistic and adaptive differences between the HLA systems, we challenged a B lymphoblastic cell line (JY), widely used as model system in studying antigen presentation, with a high temperature treatment to mimic a "fever-like state", representing one of the most common physiological responses to infection. In the absence of real invading pathogenic peptides to present, we could focus on delineating the intrinsic HLA pathway adaptations in response to high temperature in this particular cell line. Following a three-pronged approach, we performed quantitative analyses of the proteome, the HLA class I ligandome, as well as the HLA class II ligandome. The data reveals that elevated temperature may already prepare these cells for an immune-like response through increased HLA class II presentation capacity and specific release of, from the invariant chain originating, CLIP peptides. Interestingly, at high temperature, prominent changes in the composition of the CLIP repertoire were observed, with enrichment of peptides containing C-terminal extensions beyond the CLIP-core region. Collectively, these illustrate intriguing temperature sensitive adaptations in this B cell line.

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.

Anthem: a user customised tool for fast and accurate prediction of binding between peptides and HLA class I molecules

Neopeptide-based immunotherapy has been recognised as a promising approach for the treatment of cancers. For neopeptides to be recognised by CD8+ T cells and induce an immune response, their binding to human leukocyte antigen class I (HLA-I) molecules is a necessary first step. Most epitope prediction tools thus rely on the prediction of such binding. With the use of mass spectrometry, the scale of naturally presented HLA ligands that could be used to develop such predictors has been expanded. However, there are rarely efforts that focus on the integration of these experimental data with computational algorithms to efficiently develop up-to-date predictors. Here, we present Anthem for accurate HLA-I binding prediction. In particular, we have developed a user-friendly framework to support the development of customisable HLA-I binding prediction models to meet challenges associated with the rapidly increasing availability of large amounts of immunopeptidomic data. Our extensive evaluation, using both independent and experimental datasets shows that Anthem achieves an overall similar or higher area under curve value compared with other contemporary tools. It is anticipated that Anthem will provide a unique opportunity for the non-expert user to analyse and interpret their own in-house or publicly deposited datasets.

Comprehensive analysis of circular RNA in oral leukoplakia: upregulated circHLA-C as a potential biomarker for diagnosis and prognosis

Background

Emerging evidence indicates that circular RNAs (circRNAs) play an indispensable role in a variety of tumors, yet the function of circRNAs in premalignant lesions is still obscure. Oral leukoplakia (OLK) is one of the most common premalignant lesions of the oral mucosa. Our study aimed to comprehensively investigate whether circRNAs contribute to the occurrence and development of OLK.

Methods

We obtained six pairs of OLK and normal oral mucosal (NOM) tissue samples and subjected them to high-throughput sequencing to detect the expression of circRNA. In total, 26 pairs of NOM and OLK tissues were used for validation. Key circRNAs were selected and further validated by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), ribonuclease (RNase) R digestion, and Sanger sequencing. Visualization analysis of circular human leukocyte antigen-C (circHLA-C) was performed in the UCSC Genome Browser (genome.ucsc.edu). Functional analysis of differentially expressed (DE) circRNAs were processed by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Furthermore, TargetScan (www.targetscan.org) was applied to predict targeted micro RNAs (miRNAs) and messenger RNAs (mRNAs) of circRNAs and a competing endogenous RNA (ceRNA) network related with identified circRNAs was constructed in Cytoscape (v2.8.0).

Results

Profile data showed that 366 circRNAs were significantly altered in OLK tissues, including 65 upregulated and 301 downregulated circRNA transcripts. Compared with sequencing results, seven selected circRNAs expressed the same changing tendency. The amplest upregulated circRNA in our sequencing data, circHLA-C, was confirmed through back-splice junction sequences by Sanger sequencing after RNase R digestion. Correlation analysis demonstrated that circHLA-C correlated positively with the degree of dysplasia. Furthermore, receiver operating characteristic (ROC) curve analysis indicated that circHLA-C had potential diagnostic value with excellent accuracy and specificity.

Conclusions

According to the literature, we were the first to uncover the expression profiles of circRNAs in OLK. Our research performed a comprehensive bioinformatics analysis of DE circRNAs in OLK and identified circHLA-C as a promising diagnostic biomarker with potential as a therapeutic genetic target for OLK.

From Chickens to Humans: The Importance of Peptide Repertoires for MHC Class I Alleles

In humans, killer immunoglobulin-like receptors (KIRs), expressed on natural killer (NK) and thymus-derived (T) cells, and their ligands, primarily the classical class I molecules of the major histocompatibility complex (MHC) expressed on nearly all cells, are both polymorphic. The variation of this receptor-ligand interaction, based on which alleles have been inherited, is known to play crucial roles in resistance to infectious disease, autoimmunity, and reproduction in humans. However, not all the variation in response is inherited, since KIR binding can be affected by a portion of the peptide bound to the class I molecules, with the particular peptide presented affecting the NK response. The extent to which the large multigene family of chicken immunoglobulin-like receptors (ChIRs) is involved in functions similar to KIRs is suspected but not proven. However, much is understood about the two MHC-I molecules encoded in the chicken MHC. The BF2 molecule is expressed at a high level and is thought to be the predominant ligand of cytotoxic T lymphocytes (CTLs), while the BF1 molecule is expressed at a much lower level if at all and is thought to be primarily a ligand for NK cells. Recently, a hierarchy of BF2 alleles with a suite of correlated properties has been defined, from those expressed at a high level on the cell surface but with a narrow range of bound peptides to those expressed at a lower level on the cell surface but with a very wide repertoire of bound peptides. Interestingly, there is a similar hierarchy for human class I alleles, although the hierarchy is not as wide. It is a question whether KIRs and ChIRs recognize class I molecules with bound peptide in a similar way, and whether fastidious to promiscuous hierarchy of class I molecules affect both T and NK cell function. Such effects might be different from those predicted by the similarities of peptide-binding based on peptide motifs, as enshrined in the idea of supertypes. Since the size of peptide repertoire can be very different for alleles with similar peptide motifs from the same supertype, the relative importance of these two properties may be testable.

Thermostability profiling of MHC-bound peptides: a new dimension in immunopeptidomics and aid for immunotherapy design

The features of peptide antigens that contribute to their immunogenicity are not well understood. Although the stability of peptide-MHC (pMHC) is known to be important, current assays assess this interaction only for peptides in isolation and not in the context of natural antigen processing and presentation. Here, we present a method that provides a comprehensive and unbiased measure of pMHC stability for thousands of individual ligands detected simultaneously by mass spectrometry (MS). The method allows rapid assessment of intra-allelic and inter-allelic differences in pMHC stability and reveals profiles of stability that are broader than previously appreciated. The additional dimensionality of the data facilitated the training of a model which improves the prediction of peptide immunogenicity, specifically of cancer neoepitopes. This assay can be applied to any cells bearing MHC or MHC-like molecules, offering insight into not only the endogenous immunopeptidome, but also that of neoepitopes and pathogen-derived sequences.

Clinical relevance of donor-specific antibodies directed at HLA-C: A long road to acceptance

In solid organ transplantation (SOT), the clinical relevance of donor-specific antibodies (DSA) directed at anti-HLA-A, -B, -DR and -DQ antigens is largely recognized while it is still a matter of debate for DSA directed at HLA-C. In this review, we summarize the peculiarities of HLA-C among class I HLA antigens as well as their immunogenicity, which underlie the clinical relevance of HLA-C locus and anti-HLA-C DSA in SOT. Many factors, both intrinsic and extrinsic to the HLA-C gene and HLA-C protein, explain its lower expression in comparison with HLA-A and -B. This lower expression can explain the apparent lower immunogenicity of HLA-C leading to a lower prevalence and strength of anti-HLA-C antibodies. Nevertheless, HLA-C antigens are truly immunogenic and preformed anti-HLA-C DSA are clinically relevant. Indeed, anti-HLA-C DSA are able to bind donor cells and to activate the complement pathway both ex vivo and in vivo. In line with this, numerous clinical studies now show that preformed DSA directed at native HLA-C molecules induce poorer graft outcomes. We then plead for the inclusion of HLA-C in all transplant allocation systems and we propose a strategy to cope with anti-HLA-C DSA in SOT. Beyond SOT, anti-HLA-C antibodies generate a growing interest in the allo-HCT, transfusion and obstetrics fields, while new concepts such as the role of the "missing-self" in solid organ rejection places HLA-C as an inescapable actor in transplant tolerance.

Specific Class I HLA Supertypes but Not HLA Zygosity or Expression Are Associated with Outcomes following HLA-Matched Allogeneic Hematopoietic Cell Transplant: HLA Supertypes Impact Allogeneic HCT Outcomes

Maximizing the probability of antigen presentation to T cells through diversity in HLAs can enhance immune responsiveness and translate into improved clinical outcomes, as evidenced by the association of heterozygosity and supertypes at HLA class I loci with improved survival in patients with advanced solid tumors treated with immune checkpoint inhibitors. We investigated the impact of HLA heterozygosity, supertypes, and surface expression on outcomes in adult and pediatric patients with acute myeloid leukemia (AML), myelodysplastic syndrome, acute lymphoblastic leukemia, and non-Hodgkin lymphoma who underwent 8/8 HLA-matched, T cell replete, unrelated, allogeneic hematopoietic cell transplant (HCT) from 2000 to 2015 using patient data reported to the Center for International Blood and Marrow Transplant Research. HLA class I heterozygosity and HLA expression were not associated with overall survival, relapse, transplant-related mortality (TRM), disease-free survival (DFS), and acute graft-versus-host disease following HCT. The HLA-B62 supertype was associated with decreased TRM in the entire patient cohort (hazard ratio [HR], 0.79; 95% CI, 0.69 to 0.90; P = .00053). The HLA-B27 supertype was associated with worse DFS in patients with AML (HR = 1.21; 95% CI, 1.10 to 1.32; P = .00005). These findings suggest that the survival benefit of HLA heterozygosity seen in solid tumor patients receiving immune checkpoint inhibitors does not extend to patients undergoing allogeneic HCT. Certain HLA supertypes, however, are associated with TRM and DFS, suggesting that similarities in peptide presentation between supertype members play a role in these outcomes. Beyond implications for prognosis following HCT, these findings support the further investigation of these HLA supertypes and the specific immune peptides important for transplant outcomes.

Hla-C genetic diversity and evolutionary insights in two samples from Brazil and Benin

Human leukocyte antigen-C (HLA-C) is a classical HLA class I molecule that binds and presents peptides to cytotoxic T lymphocytes in the cell surface. HLA-C has a dual function because it also interacts with Killer-cell immunoglobulin-like receptors (KIR) receptors expressed in natural killer and T cells, modulating their activity. The structure and diversity of the HLA-C regulatory regions, as well as the relationship among variants along the HLA-C locus, are poorly addressed, and few population-based studies explored the HLA-C variability in the entire gene in different population samples. Here we present a molecular and bioinformatics method to evaluate the entire HLA-C diversity, including regulatory sequences. Then, we applied this method to survey the HLA-C diversity in two population samples with different demographic histories, one highly admixed from Brazil with major European contribution, and one from Benin with major African contribution. The HLA-C promoter and 3'UTR were very polymorphic with the presence of few, but highly divergent haplotypes. These segments also present conserved sequences that are shared among different primate species. Nucleotide diversity was higher in other segments rather than exons 2 and 3, particularly around exon 5 and the second half of the 3'UTR region. We detected evidence of balancing selection on the entire HLA-C locus and positive selection in the HLA-C leader peptide, for both populations. HLA-C motifs previously associated with KIR interaction and expression regulation are similar between both populations. Each allele group is associated with specific regulatory sequences, reflecting the high linkage disequilibrium along the entire HLA-C locus in both populations.

A fast and reliable method for detecting SNP rs67384697 (Hsa-miR-148a binding site) by a single run of allele-specific real-time PCR

Surface expression of human leukocyte antigen (HLA)-class I molecules is critical for modulating T/natural killer lymphocytes' effector functions. Among HLA molecules, HLA-C, the most recently evolved form of class I antigens, is subjected to both transcriptional and multiple post-transcriptional regulation mechanisms affecting its cell surface expression. Among the latter a region placed in the 3' untranslated region of HLA-C transcript contains the single nucleotide polymorphism (SNP) rs67384697 "G-ins/del" that has been found to be strictly associated with surface levels of HLA-C allomorphs because of the effect on the binding site of a microRNA (Hsa-miR-148a). Higher expression of HLA-C has been proved to influence HIV-1 infection via a better control of viremia and a slower disease progression. More importantly, the analysis of SNP rs67384697 "G-ins/del" combined with the evaluation of the HLA-Bw4/-Bw6 C1/C2 supratype, as well as the killer immunoglobulin-like receptor genetic asset, has proved to be pivotal in defining the status of Elite Controllers in the Caucasian population. Here we describe a new reliable and fast method of allele-specific real-time PCR to monitor the integrity/disruption of the binding site of the microRNA Hsa-miR-148a in a high-throughput format that can be easily applied to studies involving large cohorts of individuals.

Binding affinities of 438 HLA proteins to complete proteomes of seven pandemic viruses and distributions of strongest and weakest HLA peptide binders in populations worldwide

We report detailed peptide‐binding affinities between 438 HLA Class I and Class II proteins and complete proteomes of seven pandemic human viruses, including coronaviruses, influenza viruses and HIV‐1. We contrast these affinities with HLA allele frequencies across hundreds of human populations worldwide. Statistical modelling shows that peptide‐binding affinities classified into four distinct categories depend on the HLA locus but that the type of virus is only a weak predictor, except in the case of HIV‐1. Among the strong HLA binders (IC₅₀ ≤ 50), we uncovered 16 alleles (the top ones being A*02:02, B*15:03 and DRB1*01:02) binding more than 1% of peptides derived from all viruses, 9 (top ones including HLA‐A*68:01, B*15:25, C*03:02 and DRB1*07:01) binding all viruses except HIV‐1, and 15 (top ones A*02:01 and C*14:02) only binding coronaviruses. The frequencies of strongest and weakest HLA peptide binders differ significantly among populations from different geographic regions. In particular, Indigenous peoples of America show both higher frequencies of strongest and lower frequencies of weakest HLA binders. As many HLA proteins are found to be strong binders of peptides derived from distinct viral families, and are hence promiscuous (or generalist), we discuss this result in relation to possible signatures of natural selection on HLA promiscuous alleles due to past pathogenic infections. Our findings are highly relevant for both evolutionary genetics and the development of vaccine therapies. However they should not lead to forget that individual resistance and vulnerability to diseases go beyond the sole HLA allelic affinity and depend on multiple, complex and often unknown biological, environmental and other variables.

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-affinity oligoclonal TCRs define effective adoptive T cell therapy targeting mutant KRAS-G12D

Complete cancer regression occurs in a subset of patients following adoptive T cell therapy (ACT) of ex vivo expanded tumor-infiltrating lymphocytes (TILs). However, the low success rate presents a great challenge to broader clinical application. To provide insight into TIL-based immunotherapy, we studied a successful case of ACT where regression was observed against tumors carrying the hotspot mutation G12D in the KRAS oncogene. Four T cell receptors (TCRs) made up the TIL infusion and recognized two KRAS-G12D neoantigens, a nonamer and a decamer, all restricted by human leukocyte antigen (HLA) C*08:02. Three of them (TCR9a, 9b, and 9c) were nonamer-specific, while one was decamer-specific (TCR10). We show that only mutant G12D but not the wild-type peptides stabilized HLA-C*08:02 due to the formation of a critical anchor salt bridge to HLA-C. Therapeutic TCRs exhibited high affinities, ranging from nanomolar to low micromolar. Intriguingly, TCR binding affinities to HLA-C inversely correlated with their persistence in vivo, suggesting the importance of antigenic affinity in the function of therapeutic T cells. Crystal structures of TCR-HLA-C complexes revealed that TCR9a to 9c recognized G12D nonamer with multiple conserved contacts through shared CDR2β and CDR3α. This allowed CDR3β variation to confer different affinities via a variable HLA-C contact, generating an oligoclonal response. TCR10 recognized an induced and distinct G12D decamer conformation. Thus, this successful case of ACT included oligoclonal TCRs of high affinity recognizing distinct conformations of neoantigens. Our study revealed the potential of a structural approach to inform clinical efforts in targeting KRAS-G12D tumors by immunotherapy and has general implications for T cell-based immunotherapies.

Tuning of human NK cells by endogenous HLA-C expression

NK cells are primarily responsible for detecting malignant or pathogen-infected cells, and their function is influenced both by stress-associated activating signals and opposing inhibitory signals from receptors that recognize self MHC. The receptors that produce this inhibitory signal shift from the NKG2A:HLA-E system to that of KIR:HLA as the NK cells mature. This maturation is associated with an increase in lytic activity, as well as an increase in HLA-C protein levels controlled by the NK-specific HLA-C promoter, NK-Pro. We propose that modulation of the translatability of HLA-C transcripts in NK cells constitutes an evolutionary mechanism to control cis inhibitory signaling by HLA-C, which fine tunes NK cell activity. Furthermore, the high degree of variability in KIR receptor affinity for HLA alleles, as well as the variable expression levels of both KIR and HLA, suggest an evolutionary requirement for the tuning of NK lytic activity. Various data have demonstrated that mature NK cells may gain or lose lytic activity when placed in different environments. This indicates that NK cell activity may be more a function of constant tuning by inhibitory signals, rather than a static, irreversible "license to kill" granted to mature NK cells. Inhibitory signaling controls the filling of the cytolytic granule reservoir, which becomes depleted if there are insufficient inhibitory signals, leading to a hyporesponsive NK cell. We propose a novel model for the tuning of human NK cell activity via cis interactions in the context of recent findings on the mechanism of NK education.

A large peptidome dataset improves HLA class I epitope prediction across most of the human population

Prediction of HLA epitopes is important for the development of cancer immunotherapies and vaccines. However, current prediction algorithms have limited predictive power, in part because they were not trained on high-quality epitope datasets covering a broad range of HLA alleles. To enable prediction of endogenous HLA class I-associated peptides across a large fraction of the human population, we used mass spectrometry to profile >185,000 peptides eluted from 95 HLA-A, -B, -C and -G mono-allelic cell lines. We identified canonical peptide motifs per HLA allele, unique and shared binding submotifs across alleles and distinct motifs associated with different peptide lengths. By integrating these data with transcript abundance and peptide processing, we developed HLAthena, providing allele-and-length-specific and pan-allele-pan-length prediction models for endogenous peptide presentation. These models predicted endogenous HLA class I-associated ligands with 1.5-fold improvement in positive predictive value compared with existing tools and correctly identified >75% of HLA-bound peptides that were observed experimentally in 11 patient-derived tumor cell lines.

Tuning of NK-Specific HLA-C Expression by Alternative mRNA Splicing

A complex system regulating HLA-C expression in NK cells, driven by an NK-specific promoter that produces alternatively spliced variants of the 5'-UTR has been recently identified. Exon content of the NK-specific 5'-UTR varies strikingly across HLA-C alleles, with some exons being allele specific. In order to investigate the possibility that allelic variation in the 5'-UTR modulates HLA-C expression levels, cDNAs containing several distinct classes of 5'-UTR were compared. Subtle changes in 5'-UTR content had a significant effect on the expression of HLA-C ^* 03 and HLA-C ^* 12 cDNA clones, suggesting that alternative splicing can fine-tune the level of protein expression. The HLA-C ^* 06 allele was found to be highly expressed in relation to the other alleles studied. However, its increased expression was primarily associated with differences in the peptide-binding groove. Although the impact of allele-specific alternative splicing of NK-Pro transcripts on protein levels can be modest when compared with the effect of changes in peptide-loading, alternative splicing may represent an additional regulatory mechanism to fine-tune HLA-C levels within NK cells in distinct tissue environments or at different stages of maturation in order to achieve optimal levels of missing-self recognition.

Expression profile of HLA-B*38:55Q allele

The identification of null or questionably expressed HLA allelic variants is a major issue in HLA diagnostics, because the mistyping of the aberrant expression of such alleles can have a major impact on the outcome of both hematopoietic stem cell transplantation (HSCT) and solid organ transplants. It is debated how questionable (Q) alleles, because of their unknown expression profile, should be considered in an allogenic HSCT setting. The HLA-B*38:55Q allele was detected as an HLA-B blank specificity; DNA sequencing identified a single polymorphism at position 373 in exon 3 (TGC > CGC), which results in the replacement of cysteine 101 with an arginine in the HLA-B heavy chain, thus, impairing disulfide bridge formation in the alpha-2 domain, essential for the normal expression of the HLA molecules. In order to determine the RNA and protein expression profile of this allelic variant, we analyzed antigenic expression at different levels, transcriptional and transductional, using a combination of cellular methods, such as serological testing and flow cytometric analysis, polymerase chain reaction (PCR) sequence-specific primer (SSP) cDNA group-specific amplification and immunocytochemical assay, demonstrating the prevalent cytoplasmatic distribution of the HLA-B*38:55Q protein. Our findings suggest that in matching process the HLA-B*38:55Q allele needs to be considered as a low expressed allele, able to elicit an allogenic T-cell response in vivo and impair the transplant outcome.

The Human Immunopeptidome Project: A Roadmap to Predict and Treat Immune Diseases

The science that investigates the ensembles of all peptides associated to human leukocyte antigen (HLA) molecules is termed "immunopeptidomics" and is typically driven by mass spectrometry (MS) technologies. Recent advances in MS technologies, neoantigen discovery and cancer immunotherapy have catalyzed the launch of the Human Immunopeptidome Project (HIPP) with the goal of providing a complete map of the human immunopeptidome and making the technology so robust that it will be available in every clinic. Here, we provide a long-term perspective of the field and we use this framework to explore how we think the completion of the HIPP will truly impact the society in the future. In this context, we introduce the concept of immunopeptidome-wide association studies (IWAS). We highlight the importance of large cohort studies for the future and how applying quantitative immunopeptidomics at population scale may provide a new look at individual predisposition to common immune diseases as well as responsiveness to vaccines and immunotherapies. Through this vision, we aim to provide a fresh view of the field to stimulate new discussions within the community, and present what we see as the key challenges for the future for unlocking the full potential of immunopeptidomics in this era of precision medicine.