Population biology of antigen presentation by MHC class I molecules

Comparison of human leukocyte antigen in patients with paroxysmal nocturnal hemoglobinuria of different clone sizes

Glycosylphosphatidylinositol-anchored protein-deficient hematopoietic stem and progenitor cell development caused by PIGA mutations cannot fully explain the pathogenesis of paroxysmal nocturnal hemoglobinuria (PNH). Herein, patients newly diagnosed with PNH at our hospital between April 2019 and April 2021 were recruited. The human leukocyte antigen (HLA) class I and II loci were analyzed, and patients were stratified by PNH clone sizes: small (< 50%) and large (≥ 50%). In 40 patients (29 males; 72.5%), the median PNH clone size was 72%. Thirteen (32.5%) and twenty-seven (67.5%) patients harbored small and large PNH clones, respectively. DRB1*15:01 and DQB1*06:02 had higher frequencies in patients with PNH than in healthy controls (adjusted P-value = 4.10 × 10-4 and 4.10 × 10-4, respectively). Whole HLA class I and II allele contributions differed (P = 0.046 and 0.065, not significant difference) when comparing patients with small and large PNH clones. B*13:01 and C*04:01 allelic frequencies were significantly higher in patients with small clones (P = 0.032 and P = 0.032, respectively). Patients with small clones had higher class II HLA evolutionary divergence (HED) (P = 0.041) and global class I and II HED (P = 0.019). In the entire cohort, 17 HLA aberrations were found in 11 (27.5%) patients. No significant differences in HLA aberrations were found between patients with small or large clones. In conclusion, patients with small clones tended to have a higher frequency of immune attack-associated alleles. A higher HED in patients with small clones may reflect a propensity for T cell-mediated autoimmunity. HLA aberrations were similar between patients with small and large clones.

In silico design of high-affinity antigenic peptides for HLA-B44

Cancer cell-killing by CD8+ T cells demands effective tumor antigen presentation by human leukocyte antigen class I (HLA-I) molecules. Screening and designing highly immunogenic neoantigens require quantitative computations to reliably predict HLA-peptide binding affinities. Here, with all-atom molecular dynamics (MD) simulations and free energy perturbation (FEP) methods, we design a collection of antigenic peptide candidates through in silico mutagenesis studies on immunogenic neoantigens, yielding enhanced binding affinities to HLA-B*44:02. In-depth structural dissection shows that introducing positively charged residues such as arginine to position 6 or lysine to position 7 of the candidates triggers conformational shifts in both peptides and the antigen-binding groove of the HLA, following the "induced-fit" mechanism. Enhancement in binding affinities compared to the wild-type was found in three out of five mutated candidates. The HLA pocket, capable of accommodating positively charged residues in positions from 5 to 7, is designated as the "dynamic pocket". Taken together, we showcase an effective structure-based binding affinity optimization framework for antigenic peptides of HLA-B*44:02 and underscore the importance of dynamic nature of the antigen-binding groove in concert with the anchoring motifs. This work provides structural insights for rational design of favorable HLA-peptide bindings and future developments in neoantigen-based therapeutics.

Resolving unknown nucleotides in the IPD-IMGT/HLA database by extended and full-length sequencing of HLA class I and II alleles

In the past, identification of HLA alleles was limited to sequencing the region of the gene coding for the peptide binding groove, resulting in a lack of sequence information in the HLA database, challenging HLA allele assignment software programs. We investigated full-length sequences of 19 HLA class I and 7 HLA class II alleles, and we extended another 47 HLA class I alleles with sequences of 5' and 3' UTR regions that were all not yet available in the IPD-IMGT/HLA database. We resolved 8638 unknown nucleotides in the coding sequence of HLA class I and 2139 of HLA class II. Furthermore, with full-length sequencing of the 26 alleles, more than 90 kb of sequence information was added to the non-coding sequences, whereas extension of the 47 alleles resulted in the addition of 5.5 kb unknown nucleotides to the 5' UTR and > 31.7 kb to the 3' UTR region. With this information, some interesting features were observed, like possible recombination events and lineage evolutionary origins. The continuing increase in the availability of full-length sequences in the HLA database will enable the identification of the evolutionary origin and will help the community to improve the alignment and assignment accuracy of HLA alleles.

An immunogenetic basis for lung cancer risk

Cancer risk is influenced by inherited mutations, DNA replication errors, and environmental factors. However, the influence of genetic variation in immunosurveillance on cancer risk is not well understood. Leveraging population-level data from the UK Biobank and FinnGen, we show that heterozygosity at the human leukocyte antigen (HLA)-II loci is associated with reduced lung cancer risk in smokers. Fine-mapping implicated amino acid heterozygosity in the HLA-II peptide binding groove in reduced lung cancer risk, and single-cell analyses showed that smoking drives enrichment of proinflammatory lung macrophages and HLA-II+ epithelial cells. In lung cancer, widespread loss of HLA-II heterozygosity (LOH) favored loss of alleles with larger neopeptide repertoires. Thus, our findings nominate genetic variation in immunosurveillance as a critical risk factor for lung cancer.

The effect of human leukocyte antigen genotype on survival in advanced prostate cancer treated with primary androgen deprivation therapy: the KYUCOG-1401-A study

BACKGROUND

Immune editing, in which human leukocyte antigens (HLA) have critical roles, has been suggested to shape the landscape of human cancer. This study prospectively investigated whether HLA gene zygosity is associated with the prognosis of primary androgen deprivation therapy in advanced prostate cancer.

METHODS

KYUCOG-1401-A was conducted in conjunction with a prospective clinical trial (KYUCOG-1401). Among the patients enrolled in KYUCOG-1401 and treated with primary androgen deprivation therapy, only Japanese patients were included. HLA genotypes of HLA-A, B, C, DRB1, DQB1, and DPB1 were determined. The effect of divergence of HLA genotypes on time to progression, prostate cancer-specific survival, and overall survival was evaluated.

RESULTS

Among 127 patients, homozygosity for HLA-DRB1 (HR, 95% CI; 4.05, 1.54-10.7, P = 0.0047) and HLA-DQB1 (HR, 95% CI; 3.75, 1.47-9.58, P = 0.0058) was associated with an increased risk of prostate cancer-specific mortality. Patients with higher HLA evolutionary divergence scores at HLA-DQB1 (HR, 95% CI; 0.90, 0.82-0.97, P = 0.0093) had lower risks of prostate cancer-specific mortality. Androgen-responsive gene sets were upregulated in CD4low and CD8low tumors in the prostate cancer cohort, but not in the bladder and kidney cancer cohorts.

CONCLUSIONS

This study suggested that the diversity of HLA-II loci including HLA-DRB1 and HLA-DQB1 plays an important role in advanced prostate cancer survival, contributing to improved risk stratification in advanced prostate cancer. Moreover, it was shown that CD4+ T cells play an important role in androgen deprivation therapy, suggesting that immunotherapy targeting CD4+ T cells is promising for prostate cancer.

Assessment of human leukocyte antigen-based neoantigen presentation to determine pan-cancer response to immunotherapy

Despite the central role of human leukocyte antigen class I (HLA-I) in tumor neoantigen presentation, quantitative determination of presentation capacity remains elusive. Based on a pooled pan-cancer genomic dataset of 885 patients treated with immune checkpoint inhibitors (ICIs), we developed a score integrating the binding affinity of neoantigens to HLA-I, as well as HLA-I allele divergence, termed the HLA tumor-Antigen Presentation Score (HAPS). Patients with a high HAPS were more likely to experience survival benefit following ICI treatment. Analysis of the tumor microenvironment indicated that the antigen presentation pathway was enriched in patients with a high HAPS. Finally, we built a neural network incorporating factors associated with neoantigen production, presentation, and recognition, which exhibited potential for differentiating cancer patients likely to benefit from ICIs. Our findings highlight the clinical utility of evaluating HLA-I tumor antigen presentation capacity and describe how ICI response may depend on HLA-mediated immunity.

Large scale transcriptional analysis of MHC class I haplotype diversity in sheep

Domestic sheep (Ovis aries) have been an important component of livestock agricultural production for thousands of years. Preserving genetic diversity within livestock populations maintains a capacity to respond to changing environments and rapidly evolving pathogens. MHC genetic diversity can influence immune functionality at individual and population levels. Here, we focus on defining functional MHC class I haplotype diversity in a large cohort of Scottish Blackface sheep pre-selected for high levels of MHC class II DRB1 diversity. Using high-throughput amplicon sequencing with three independent sets of barcoded primers we identified 134 MHC class I transcripts within 38 haplotypes. Haplotypes were identified with between two and six MHC class I genes, plus variable numbers of conserved sequences with very low read frequencies. One or two highly transcribed transcripts dominate each haplotype indicative of two highly polymorphic, classical MHC class I genes. Additional clusters of medium, low, and very low expressed transcripts are described, indicative of lower transcribed classical, non-classical and genes whose function remains to be determined.

Specificity of HLA monoclonal antibodies and their use to determine HLA expression on lymphocytes and peripheral blood stem cells

HLA Class I and II expression are known to differ locus-to-locus, however, HLA expression on the cell-surface is frequently reported as the total amount of HLA Class I or II antigens. This is despite evidence that indicates the differential expression of HLA can influence patient outcomes post-transplantation. Although numerous commercially available HLA monoclonal antibodies (mAbs) exist to characterize HLA expression, there is currently a lack of detailed information regarding their reactivities to HLA specificities. The specificities of locus-specific HLA mAbs (nine Class I and four Class II mAbs) were evaluated by two solid-phase Luminex single antigen bead assays. The reactivity patterns of these mAbs were then confirmed by flow cytometry using lymphocytes and PBSCs (peripheral blood stem cells). Out of the 13 HLA mAbs tested, only four (one Class I and three Class II mAbs) displayed intra-locus reactivity without also reacting to inter-locus specificities. Epitope analysis revealed the presence of shared epitopes across numerous HLA loci, explaining much of the observed inter-locus reactivity. The specificity of the HLA mAbs seen in solid-phase assays was confirmed against PBSCs and lymphocytes by flow cytometry. Using this method, we observed differences in the cell surface expression of HLA-C, HLA-DR, HLA-DQ, and HLA-DP between PBSCs and lymphocytes. Our results emphasize the need to characterize the reactivity patterns of HLA mAbs using solid-phase assays before their use on cells. Through understanding the reactivity of these HLA mAbs, the cellular expression of HLA can be more accurately assessed in downstream assays.

Master Transcription Factor Reprogramming Unleashes Selective Translation Promoting Castration Resistance and Immune Evasion in Lethal Prostate Cancer

Signaling rewiring allows tumors to survive therapy. Here we show that the decrease of the master regulator microphthalmia transcription factor (MITF) in lethal prostate cancer unleashes eukaryotic initiation factor 3B (eIF3B)-dependent translation reprogramming of key mRNAs conferring resistance to androgen deprivation therapy (ADT) and promoting immune evasion. Mechanistically, MITF represses through direct promoter binding eIF3B, which in turn regulates the translation of specific mRNAs. Genome-wide eIF3B enhanced cross-linking immunoprecipitation sequencing (eCLIP-seq) showed specialized binding to a UC-rich motif present in subsets of 5' untranslated regions. Indeed, translation of the androgen receptor and major histocompatibility complex I (MHC-I) through this motif is sensitive to eIF3B amount. Notably, pharmacologic targeting of eIF3B-dependent translation in preclinical models sensitizes prostate cancer to ADT and anti-PD-1 therapy. These findings uncover a hidden connection between transcriptional and translational rewiring promoting therapy-refractory lethal prostate cancer and provide a druggable mechanism that may transcend into effective combined therapeutic strategies.

SIGNIFICANCE

Our study shows that specialized eIF3B-dependent translation of specific mRNAs released upon downregulation of the master transcription factor MITF confers castration resistance and immune evasion in lethal prostate cancer. Pharmacologic targeting of this mechanism delays castration resistance and increases immune-checkpoint efficacy. This article is featured in Selected Articles from This Issue, p. 2489.

The KIR repertoire of a West African chimpanzee population is characterized by limited gene, allele, and haplotype variation

Introduction

The killer cell immunoglobulin-like receptors (KIR) play a pivotal role in modulating the NK cell responses, for instance, through interaction with major histocompatibility complex (MHC) class I molecules. Both gene systems map to different chromosomes but co-evolved during evolution. The human KIR gene family is characterized by abundant allelic polymorphism and copy number variation. In contrast, our knowledge of the KIR repertoire in chimpanzees is limited to 39 reported alleles, with no available population data. Only three genomic KIR region configurations have been mapped, and seventeen additional ones were deduced by genotyping.

Methods

Previously, we documented that the chimpanzee MHC class I repertoire has been skewed due to an ancient selective sweep. To understand the depth of the sweep, we set out to determine the full-length KIR transcriptome - in our MHC characterized pedigreed West African chimpanzee cohort - using SMRT sequencing (PacBio). In addition, the genomic organization of 14 KIR haplotypes was characterized by applying a Cas9-mediated enrichment approach in concert with long-read sequencing by Oxford Nanopore Technologies.

Results

In the cohort, we discovered 35 undescribed and 15 already recorded Patr-KIR alleles, and a novel hybrid KIR gene. Some KIR transcripts are subject to evolutionary conserved alternative splicing events. A detailed insight on the KIR region dynamics (location and order of genes) was obtained, however, only five new KIR region configurations were detected. The population data allowed to investigate the distribution of the MHC-C1 and C2-epitope specificity of the inhibitory lineage III KIR repertoire, and appears to be skewed towards C2.

Discussion

Although the KIR region is known to evolve fast, as observed in other primate species, our overall conclusion is that the genomic architecture and repertoire in West African chimpanzees exhibit only limited to moderate levels of variation. Hence, the ancient selective sweep that affected the chimpanzee MHC class I region may also have impacted the KIR system.

Correlative CD4 and CD8 T-cell immunodominance in humans and mice: Implications for preclinical testing

Antigen-specific T-cell recognition is restricted by Major Histocompatibility Complex (MHC) molecules, and differences between CD4 and CD8 immunogenicity in humans and animal species used in preclinical vaccine testing are yet to be fully understood. In this study, we addressed this matter by analyzing experimentally identified epitopes based on published data curated in the Immune Epitopes DataBase (IEDB) database. We first analyzed SARS-CoV-2 spike (S) and nucleoprotein (N), which are two common targets of the immune response and well studied in both human and mouse systems. We observed a weak but statistically significant correlation between human and H-2b mouse T-cell responses (CD8 S specific (r = 0.206, p = 1.37 × 10-13); CD4 S specific (r = 0.118, p = 2.63 × 10-5) and N specific (r = 0.179, p = 2.55 × 10-4)). Due to intrinsic differences in MHC molecules across species, we also investigated the association between the immunodominance of common Human Leukocyte Antigen (HLA) alleles for which HLA transgenic mice are available, namely, A*02:01, B*07:02, DRB1*01:01, and DRB1*04:01, and found higher significant correlations for both CD8 and CD4 (maximum r = 0.702, p = 1.36 × 10-31 and r = 0.594, p = 3.04-122, respectively). Our results further indicated that some regions are commonly immunogenic between humans and mice (either H-2b or HLA transgenic) but that others are human specific. Finally, we noted a significant correlation between CD8 and CD4 S- (r = 0.258, p = 7.33 × 1021) and N-specific (r = 0.369, p = 2.43 × 1014) responses, suggesting that discrete protein subregions can be simultaneously recognized by T cells. These findings were confirmed in other viral systems, providing general guidance for the use of murine models to test T-cell immunogenicity of viral antigens destined for human use.

The Prognostic Significance of Selected HLA Alleles on Prostate Cancer Outcome

Recently, we have shown that HLA-A*02:01 and HLA-A*24:02 in de novo metastatic prostate cancer (MPCa) have an important role in disease progression. Since de novo MPCa represents a small group among patients diagnosed with prostate cancer (PCa), it was obvious to try to extend the validity of our results to larger cohorts of PCa patients. Herein, we analyzed patients irrespective of their disease status at diagnosis to include, besides patients with MPCa, those with localized PCa (LPCa). Our goal was to specify the prognostic value of HLA-A*02:01 and HLA-A*24:02 for overall survival (OS) prospectively and for early biochemical recurrence (BCR) and castrate resistance (CR) as additional clinical endpoints in a prospective/retrospective manner, to improve clinical decisions for patients covering all stages of PCa. On univariate analysis, HLA-A alleles were significantly associated as prognostic biomarkers with early BCR (p = 0.028; HR = 1.822), OS (p = 0.013; HR = 1.547) and showed a trend for CR (p = 0.150; HR = 1.239). On multivariate analysis, HLA-A alleles proved to be independent prognosticators for early BCR (p = 0.017; HR = 2.008), CR (p = 0.005; HR = 1.615), and OS (p = 0.002; HR = 2.063). Kaplan-Meier analyses revealed that patients belonging to the HLA-A*02:01+HLA-A*24:02- group progressed much faster to BCR and CR and had also shorter OS compared to HLA-A*24:02+ patients. Patients being HLA-A*02:01-HLA-A*24:02- exhibited varying clinical outcomes, pointing to the presence of additional HLA-A alleles with potential prognostic value. Our data underline the HLA-A alleles as valuable prognostic biomarkers for PCa that may assist with the appropriate treatment and follow-up schedule based on the risk for disease progression to avoid over-diagnosis and over-treatment.

Antigen presentation in cancer - mechanisms and clinical implications for immunotherapy

Over the past decade, the emergence of effective immunotherapies has revolutionized the clinical management of many types of cancers. However, long-term durable tumour control is only achieved in a fraction of patients who receive these therapies. Understanding the mechanisms underlying clinical response and resistance to treatment is therefore essential to expanding the level of clinical benefit obtained from immunotherapies. In this Review, we describe the molecular mechanisms of antigen processing and presentation in tumours and their clinical consequences. We examine how various aspects of the antigen-presentation machinery (APM) shape tumour immunity. In particular, we discuss genomic variants in HLA alleles and other APM components, highlighting their influence on the immunopeptidomes of both malignant cells and immune cells. Understanding the APM, how it is regulated and how it changes in tumour cells is crucial for determining which patients will respond to immunotherapy and why some patients develop resistance. We focus on recently discovered molecular and genomic alterations that drive the clinical outcomes of patients receiving immune-checkpoint inhibitors. An improved understanding of how these variables mediate tumour-immune interactions is expected to guide the more precise administration of immunotherapies and reveal potentially promising directions for the development of new immunotherapeutic approaches.

Complete sequences of six major histocompatibility complex haplotypes, including all the major MHC class II structures

Accurate and comprehensive immunogenetic reference panels are key to the successful implementation of population-scale immunogenomics. The 5Mbp Major Histocompatibility Complex (MHC) is the most polymorphic region of the human genome and associated with multiple immune-mediated diseases, transplant matching and therapy responses. Analysis of MHC genetic variation is severely complicated by complex patterns of sequence variation, linkage disequilibrium and a lack of fully resolved MHC reference haplotypes, increasing the risk of spurious findings on analyzing this medically important region. Integrating Illumina, ultra-long Nanopore, and PacBio HiFi sequencing as well as bespoke bioinformatics, we completed five of the alternative MHC reference haplotypes of the current (GRCh38/hg38) build of the human reference genome and added one other. The six assembled MHC haplotypes encompass the DR1 and DR4 haplotype structures in addition to the previously completed DR2 and DR3, as well as six distinct classes of the structurally variable C4 region. Analysis of the assembled haplotypes showed that MHC class II sequence structures, including repeat element positions, are generally conserved within the DR haplotype supergroups, and that sequence diversity peaks in three regions around HLA-A, HLA-B+C, and the HLA class II genes. Demonstrating the potential for improved short-read analysis, the number of proper read pairs recruited to the MHC was found to be increased by 0.06%-0.49% in a 1000 Genomes Project read remapping experiment with seven diverse samples. Furthermore, the assembled haplotypes can serve as references for the community and provide the basis of a structurally accurate genotyping graph of the complete MHC region.

Wild Boar (Sus scrofa)—Fascioloides magna Interaction from the Perspective of the MHC Genes

Fascioloidosis is a parasitic disease caused by a trematode Fascioloides magna. Since major histocompatibility complex (MHC) genes play an important role in the immune response, the aim of this study was to compare the potential differences in MHC class II SLA-DRB1 exon 2 genes between wild boar populations from infected (cases) and non-infected areas (controls). During the winter of 2021, a total of 136 wild boar tissue samples were collected, 39 cases and 97 controls. DNA was extracted and sequenced using the Illumina platform. Differences in distributions of allele combinations were calculated using the Chi-Square test for homogeneity and between proportions using the large-sample test and Fisher–Irwin test. Analysis revealed 19 previously described swine leucocyte antigen (SLA) alleles. The number of polymorphic sites was 79 (29.6%), with 99 mutations in total. Nucleotide diversity π was estimated at 0.11. Proportions of the alleles SLA-DRB1*12:05 (p = 0.0008379) and SLA-DRB1*0101 (p = 0.0002825) were statistically significantly higher in controls, and proportions of the SLA-DRB1*0602 (p = 0.006059) and SLA-DRB1*0901 (p = 0.0006601) in cases. Alleles SLA-DRB1*04:09, SLA-DRB1*0501, SLA-DRB1*11:09, and SLA-DRB1*1301 were detected only in cases, while SLA-DRB1*0404, SLA-DRB1*0701, SLA-DRB1*02:10, and SLA-DRB1*04:08 were present only in controls. We did not confirm the existence of specific alleles that could be linked to F. magna infection. Detected high variability of the MHC class II SLA-DRB1 exon 2 genes indicate high resistance potential against various pathogens.

A Review of Suggested Mechanisms of MHC Odor Signaling

Simple Summary

Immune genes of the vertebrate MHC vary among individuals. Each individual collection is optimally diverse to provide resistance against some infectious diseases but not too diverse to cause autoimmune diseases. MHC-dependent mate choice aims for optimally complementary MHC alleles. Each potential partner signals through body odor his/her MHC alleles. Identifying the signal molecules was a long-lasting puzzle solved only recently after many deviations as described. Commensal microbiota which are controlled by the individual MHC genes differ among individuals. They were suspected repeatedly to provide the signal, though mice raised germ-free could still smell MHC genes. Carrier hypotheses came in various versions, centered around the specificity of each MHC molecule for binding peptides from diseases, shown to T lymphocytes to induce the immune response. Volatiles of various origins were suggested to fill the place of the peptide and thus reflect the identity of the MHC molecule. Finally, the bound peptides themselves were identified as the sought info-chemicals. Synthesized peptides affect mate choice as predicted. Specific olfactory neurons were shown to react to these peptides but only to the anchors that define the binding specificity. Even eggs choose sperm to produce offspring with optimal MHC, though the signaling pathway needs further research.

Abstract

Although an individual’s mix of MHC immune genes determines its resistance, finding MHC-dependent mate choice occurred by accident in inbred mice. Inbred mice prefer MHC dissimilar mates, even when the choice was restricted to urine. It took decades to find the info-chemicals, which have to be as polymorphic as the MHC. Microbiota were suggested repeatedly as the origin of the odor signal though germ-free mice maintained normal preference. Different versions of the ‘carrier hypothesis’ suggested MHC molecules carry volatiles after the bound peptide is released. Theory predicted an optimal individual MHC diversity to maximize resistance. The optimally complementary mate should be and is preferred as several studies show. Thus, the odor signal needs to transmit the exact information of the sender’s MHC alleles, as do MHC ligand peptides but not microbiota. The ‘MHC peptide hypothesis’ assumes that olfactory perception of the peptide ligand provides information about the MHC protein in a key-lock fashion. Olfactory neurons react only to the anchors of synthesized MHC peptides, which reflect the binding MHC molecule’s identity. Synthesized peptides supplemented to a male’s signal affect choice in the predicted way, however, not when anchors are mutated. Also, the human brain detects smelled synthesized self-peptides as such. After mate choice, the lottery of meiosis of randomly paired oocyte and sperm haplotypes would often produce MHC non-optimal offspring. In sticklebacks, eggs select MHC-compatible sperm, thus prefer the best combination close to the population optimum.

Immunogenomics of Killer Cell Immunoglobulin-Like Receptor (KIR) and HLA Class I: Coevolution and Consequences for Human Health

Interactions of killer cell immunoglobin-like receptors (KIR) with human leukocyte antigens (HLA) class I regulate effector functions of key cytotoxic cells of innate and adaptive immunity. The extreme diversity of this interaction is genetically determined, having evolved in the ever-changing environment of pathogen exposure. Diversity of KIR and HLA genes is further facilitated by their independent segregation on separate chromosomes. That fetal implantation relies on many of the same types of immune cells as infection control places certain constraints on the evolution of KIR interactions with HLA. Consequently, specific inherited combinations of receptors and ligands may predispose to specific immune-mediated diseases, including autoimmunity. Combinatorial diversity of KIR and HLA class I can also differentiate success rates of immunotherapy directed to these diseases. Progress toward both etiopathology and predicting response to therapy is being achieved through detailed characterization of the extent and consequences of the combinatorial diversity of KIR and HLA. Achieving these goals is more tractable with the development of integrated analyses of molecular evolution, function, and pathology that will establish guidelines for understanding and managing risks. Here, we present what is known about the coevolution of KIR with HLA class I and the impact of their complexity on immune function and homeostasis.

Protective association of HLA-DPB1*04:01:01 with acute encephalopathy with biphasic seizures and late reduced diffusion identified by HLA imputation

Acute encephalopathy with biphasic seizures and late reduced diffusion (AESD) is a severe syndrome of acute encephalopathy that affects infants and young children. AESD is a polygenic disorder preceded by common viral infections with high fever. We conducted an association study of human leukocyte antigen (HLA) regions with AESD using HLA imputation. SNP genotyping was performed on 254 Japanese patients with AESD and 799 healthy controls. We conducted 3-field HLA imputation for 14 HLA genes based on Japanese-specific references using data from our previous genome-wide association study. After quality control, 208 patients and 737 controls were included in the analysis of HLA alleles. We then compared the carrier frequencies of HLA alleles and haplotypes between the patients and controls. HLA-DPB1*04:01:01 showed a significant association with AESD, exerting a protective effect against the disease (p = 0.0053, p_corrected = 0.042, odds ratio = 0.43, 95% confidence interval = 0.21-0.80). The allele frequency of HLA-DPB1*04:01:01 was lower in East Asians than in Caucasians, which may partially account for the higher incidence of AESD in the Japanese population. The present results demonstrate the importance of fine-mapping of the HLA region to investigate disease susceptibilities and elucidate the pathogenesis of AESD.

Predictive Markers for Immune Checkpoint Inhibitors in Non-Small Cell Lung Cancer

Immune checkpoint inhibitors (ICIs) have dramatically improved the outcomes of non-small cell lung cancer patients and have increased the possibility of long-term survival. However, few patients benefit from ICIs, and no predictive biomarkers other than tumor programmed cell death ligand 1 (PD-L1) expression have been established. Hence, the identification of biomarkers is an urgent issue. This review outlines the current understanding of predictive markers for the efficacy of ICIs, including PD-L1, tumor mutation burden, DNA mismatch repair deficiency, microsatellite instability, CD8⁺ tumor-infiltrating lymphocytes, human leukocyte antigen class I, tumor/specific genotype, and blood biomarkers such as peripheral T-cell phenotype, neutrophil-to-lymphocyte ratio, interferon-gamma, and interleukin-8. A tremendous number of biomarkers are in development, but individual biomarkers are insufficient. Tissue biomarkers have issues in reproducibility and accuracy because of intratumoral heterogeneity and biopsy invasiveness. Furthermore, blood biomarkers have difficulty in reflecting the tumor microenvironment and therefore tend to be less predictive for the efficacy of ICIs than tissue samples. In addition to individual biomarkers, the development of composite markers, including novel technologies such as machine learning and high-throughput analysis, may make it easier to comprehensively analyze multiple biomarkers.

Characterization of a novel HLA-A*11:335 allele resulting from a rare interlocus recombination involving HLA-A*11:01:01:01/126 and HLA-H*02:07/14/18 alleles with nanopore sequencing, in a volunteer from the China Marrow Donor Program

BACKGROUND

The major histocompatibility complex (MHC) in humans includes three classical class I loci (A, B, and C), which are important biomarkers for the transplantation of organs and hematopoietic stem cells. In the MHC, polymorphism is known to be extremely high while interlocus recombination is rare. We report a rare interlocus recombination between HLA-A and HLA-H, which was analyzed using next generation sequencing and nanopore sequencing.

METHODS

In the sample, the genotypes of HLA-A, B, C, DRB1, and DQB1 were firstly determined using the methods of sequence-specific primer, sequence-specific oligonucleotide, Sanger's sequencing, and NGS; however, HLA-A could not be phased. Nanopore sequencing was finally utilized to distinguish the sequence of the novel allele.

RESULTS

Finally, the novel HLA-A*11:335 allele was identified as an interlocus recombination involving HLA-A*11:01:01:01/126 and HLA-H*02:07/14/18 alleles; this was mainly achieved by nanopore sequencing.

CONCLUSIONS

The identification of the interlocus recombination indicated that nanopore sequencing can be helpful in the characterization of novel alleles with complex rearrangements. Interlocus recombination has been identified as one of the mechanisms involved in the generation of novel HLA alleles.

Class I/Class II HLA Evolutionary Divergence Ratio Is an Independent Marker Associated With Disease-Free and Overall Survival After Allogeneic Hematopoietic Stem Cell Transplantation for Acute Myeloid Leukemia

Class I Human Leukocyte Antigen (HLA) evolutionary divergence (HED) is a metric which reflects immunopeptidome diversity and has been associated with immune checkpoint inhibitor responses in solid tumors. Its impact and interest in allogeneic hematopoietic stem cell transplantation (HCT) have not yet been thoroughly studied. This study analyzed the clinical and immune impact of class I and II HED in 492 acute myeloid leukemia (AML) recipients undergoing HCT. The overall cohort was divided into a training (n=338) and a testing (n=132) set. Univariate cox screening found a positive impact of a high class I HED and a negative impact of a high class II HED on both disease-free (DFS) and overall survival (OS). These results were combined in a unique marker, class I/class II HED ratio, and assessed in the testing cohort. The final multivariate cox model confirmed the positive impact of a high versus low class I/class II HED ratio on both DFS (Hazard Ratio (HR) 0.41 [95% CI 0.2-0.83]; p=0.01) and OS (HR 0.34 [0.19-0.59]; p<0.001), independently of HLA matching and other HCT parameters. No significant association was found between the ratio and graft-versus-host disease (GvHD) nor with neutrophil and platelet recovery. A high class I HED was associated with a tendency for an increase in NK, CD8 T-cell, and B cell recovery at 12 months. These results introduce HED as an original and independent prognosis marker reflecting immunopeptidome diversity and alloreactivity after HCT.

High Resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1, and ∼DQB1 Diversity in South African Populations

Background: Lack of HLA data in southern African populations hampers disease association studies and our understanding of genetic diversity in these populations. We aimed to determine HLA diversity in South African populations using high resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1 and ∼DQB1 data, from 3005 previously typed individuals.

Methods: We determined allele and haplotype frequencies, deviations from Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD) and neutrality test. South African HLA class I data was additionally compared to other global populations using non-metrical multidimensional scaling (NMDS), genetic distances and principal component analysis (PCA).

Results: All loci strongly (p < 0.0001) deviated from HWE, coupled with excessive heterozygosity in most loci. Two of the three most frequent alleles, HLA ∼DQA1*05:02 (0.2584) and HLA ∼C*17:01 (0.1488) were previously reported in South African populations at lower frequencies. NMDS showed genetic distinctness of South African populations. Phylogenetic analysis and PCA clustered our current dataset with previous South African studies. Additionally, South Africans seem to be related to other sub-Saharan populations using HLA class I allele frequencies.

Discussion and Conclusion: Despite the retrospective nature of the study, data missingness, the imbalance of sample sizes for each locus and haplotype pairs, and induced methodological difficulties, this study provides a unique and large HLA dataset of South Africans, which might be a useful resource to support anthropological studies, disease association studies, population based vaccine development and donor recruitment programs. We additionally provide simulated high resolution HLA class I data to augment the mixed resolution typing results generated from this study.

hDirect-MAP: projection-free single-cell modeling of response to checkpoint immunotherapy

There is a lack of robust generalizable predictive biomarkers of response to immune checkpoint blockade in multiple types of cancer. We develop hDirect-MAP, an algorithm that maps T cells into a shared high-dimensional (HD) expression space of diverse T cell functional signatures in which cells group by the common T cell phenotypes rather than dimensional reduced features or a distorted view of these features. Using projection-free single-cell modeling, hDirect-MAP first removed a large group of cells that did not contribute to response and then clearly distinguished T cells into response-specific subpopulations that were defined by critical T cell functional markers of strong differential expression patterns. We found that these grouped cells cannot be distinguished by dimensional-reduction algorithms but are blended by diluted expression patterns. Moreover, these identified response-specific T cell subpopulations enabled a generalizable prediction by their HD metrics. Tested using five single-cell RNA-seq or mass cytometry datasets from basal cell carcinoma, squamous cell carcinoma and melanoma, hDirect-MAP demonstrated common response-specific T cell phenotypes that defined a generalizable and accurate predictive biomarker.

HLA polymorphisms and risk of glioblastoma in Koreans

Purpose

Immune responses for cancer cells can be altered according to genetic variation of human leukocyte antigen (HLA). Association of HLA polymorphism with risk of various cancer types is well known. However, the association between HLA and glioblastoma (GBM) remains uncertain. We sought to evaluate the association of HLA polymorphism with risk of GBM development in Koreans.

Materials and methods

A case-control study was performed to identify the odds ratios (OR) of HLA class I and II genes for GBM. The control group consisted of 142 healthy Korean volunteers, and the GBM group was 80 patients with newly diagnosed GBM at our institution. HLA class I (-A, -B, and–C) and class II (-DR, -DQ, and–DP) genotyping was performed by high-resolution polymerase chain reaction (PCR)-sequence-based typing (PCR-SBT) methods.

Results

There were significantly decreased frequencies of HLA-A*26:02 (OR 0.22 CI 0.05–0.98), HLA-C*08:01 (OR 0.29 CI 0.10–0.87), and HLA-DRB1*08:03 (OR 0.32 CI 0.11–0.98), while there was significantly increased frequency of HLA-C*04:01 (OR 2.29 CI 1.05–4.97). In analysis of haplotypes, the frequency of DRB1*14:05-DQB1*05:03 was significantly decreased (OR 0.22 CI 0.05–0.98).

Conclusion

This study suggests that genetic variations of HLA may affect GBM development in Koreans. Further investigations with larger sample sizes are needed to delineate any potential role of the HLA polymorphisms in the pathogenesis of GBM development.

HLA Expression in Relation to HLA Type in Classic Hodgkin Lymphoma Patients

Simple Summary

Classic Hodgkin lymphoma (cHL) is a B-cell malignancy with involvement of Epstein–Barr virus (EBV) in about 30% of the European population. The risk to develop cHL is strongly linked to genetic variants in the human leukocyte antigen (HLA) genomic region and to certain HLA alleles. This may be caused by the function of HLA alleles, or by genetic linkage to non-HLA genes. HLA can present EBV-derived and tumour-cell specific antigens and this may lead to anti-tumour immune responses. However, the tumour cells downregulate HLA expression in a proportion of the cases, which may result in immune escape. In this study, we tested whether the loss of HLA expression is related to the presence of certain protective HLA alleles. We found that loss and retention of HLA expression is indeed associated with presence of known susceptibility HLA alleles. These findings suggest that HLA itself is involved in development of cHL.

Abstract

Several human leukocyte antigen (HLA) alleles are strongly associated with susceptibility to classic Hodgkin lymphoma (cHL), also in subgroups stratified for presence of the Epstein–Barr virus (EBV). We tested the hypothesis that the pressure on cHL tumour cells to lose HLA expression is associated with HLA susceptibility alleles. A meta-analysis was carried out to identify consistent protective and risk HLA alleles in a combined cohort of 839 cHL patients from the Netherlands and the United Kingdom. Tumour cell HLA expression was studied in 338 cHL cases from these two cohorts and correlated to the presence of specific susceptibility HLA alleles. Carriers of the HLA-DRB1*07 protective allele frequently lost HLA class II expression in cHL overall. Patients carrying the HLA-DRB1*15/16 (DR2) risk allele retained HLA class II expression in EBV− cHL and patients with the HLA-B*37 risk allele retained HLA class I expression more frequently than non-carriers in EBV+ cHL. The other susceptibility alleles showed no significant differences in expression. Thus, HLA expression by tumour cells is associated with a subset of the protective and risk alleles. This strongly suggests that HLA associations in cHL are related to peptide binding capacities of specific HLA alleles.

The Genomic Organization of the LILR Region Remained Largely Conserved Throughout Primate Evolution: Implications for Health And Disease

The genes of the leukocyte immunoglobulin-like receptor (LILR) family map to the leukocyte receptor complex (LRC) on chromosome 19, and consist of both activating and inhibiting entities. These receptors are often involved in regulating immune responses, and are considered to play a role in health and disease. The human LILR region and evolutionary equivalents in some rodent and bird species have been thoroughly characterized. In non-human primates, the LILR region is annotated, but a thorough comparison between humans and non-human primates has not yet been documented. Therefore, it was decided to undertake a comprehensive comparison of the human and non-human primate LILR region at the genomic level. During primate evolution the organization of the LILR region remained largely conserved. One major exception, however, is provided by the common marmoset, a New World monkey species, which seems to feature a substantial contraction of the number of LILR genes in both the centromeric and the telomeric region. Furthermore, genomic analysis revealed that the killer-cell immunoglobulin-like receptor gene KIR3DX1, which maps in the LILR region, features one copy in humans and great ape species. A second copy, which might have been introduced by a duplication event, was observed in the lesser apes, and in Old and New World monkey species. The highly conserved gene organization allowed us to standardize the LILR gene nomenclature for non-human primate species, and implies that most of the receptors encoded by these genes likely fulfill highly preserved functions.

NSCLC Biomarkers to Predict Response to Immunotherapy with Checkpoint Inhibitors (ICI): From the Cells to In Vivo Images

Lung cancer remains the leading cause of cancer-related death, and it is usually diagnosed in advanced stages (stage III or IV). Recently, the availability of targeted strategies and of immunotherapy with checkpoint inhibitors (ICI) has favorably changed patient prognosis. Treatment outcome is closely related to tumor biology and interaction with the tumor immune microenvironment (TME). While the response in molecular targeted therapies relies on the presence of specific genetic alterations in tumor cells, accurate ICI biomarkers of response are lacking, and clinical outcome likely depends on multiple factors that are both host and tumor-related. This paper is an overview of the ongoing research on predictive factors both from in vitro/ex vivo analysis (ranging from conventional pathology to molecular biology) and in vivo analysis, where molecular imaging is showing an exponential growth and use due to technological advancements and to the new bioinformatics approaches applied to image analyses that allow the recovery of specific features in specific tumor subclones.

Effect of HLA genotype on intravesical recurrence after bacillus Calmette-Guérin therapy for non-muscle-invasive bladder cancer

The intravesical administration of bacillus Calmette-Guérin (BCG) is widely used to control the intravesical recurrence of non-muscle-invasive bladder cancer (NMIBC). This study aimed to reveal the effects of zygosity on human leukocyte antigen (HLA) genes and individual HLA genotypes on intravesical recurrence after intravesical BCG therapy for NMIBC. This study included Japanese patients who had received intravesical BCG for NMIBC. HLA genotyping of HLA-A, B, C, and DRB1 was performed. The effect of HLA zygosity and HLA genotype on intravesical recurrence was evaluated. Among 195 patients, those homozygous for the HLA-B supertype were more likely than those heterozygous for the HLA-B supertype to experience intravesical recurrence by univariate analysis (hazard ratio [HR], 95% confidence interval [CI]; 1.87, 1.14-3.05, P = 0.012) and multivariate analysis (HR, 95% CI; 2.26, 1.02-5.01, P = 0.045). Patients with B07 or B44 had a decreased risk of intravesical recurrence by univariate analysis (HR, 95% CI; 0.43, 0.24-0.78, P = 0.0056) and multivariate analysis (HR, 95% CI; 0.36, 0.16-0.82, P = 0.016). This study suggests the importance of the diversity and specificity of HLA-B loci in the antitumor effect of BCG immunotherapy for NMIBC. These findings may contribute to the delineation of risk strata for BCG therapy and improve the medical management of NMIBC.

Antigen Processing, Presentation, and Tolerance: Role in Autoimmune Skin Diseases

Autoreactive T cells pose a constant risk for the emergence of autoimmune skin diseases in genetically predisposed individuals carrying certain HLA risk alleles. Immune tolerance mechanisms are opposed by broad HLA-presented self-immunopeptidomes, a predefined repertoire of polyspecific TCRs, the continuous generation of new antibody specificities by somatic recombination of Ig genes in B cells, and heightened proinflammatory reactivity. Increased autoantigen presentation by HLA molecules, cross-activation of pathogen-induced T cells against autologous structures, altered metabolism of self-proteins, and excessive production of proinflammatory signals may all contribute to the breakdown of immune tolerance and the development of autoimmune skin diseases.

The prevalence of HLA-I LOH in Chinese pan-cancer patients and genomic features of patients harboring HLA-I LOH

HLA-I LOH may facilitate immune evasion. However, large population studies on the prevalence of HLA-I LOH across different cancer types and in relation to mutational profiles are lacking, in particular, in the Chinese population. In this study, analysis was performed in 1504 advanced pan-cancer patients and 134 early-stage non-small-cell lung cancer patients using a 1021-gene panel. The consistency between the 1021-gene panel and whole-exome sequencing was evaluated in 45 samples, where concordant results were obtained in 95.6% (43/45) of the samples. Analytical results revealed that the prevalence of HLA-I LOH in tumor tissue presents considerable differences across cancer types. HLA-I LOH was relevant to genomic instability, reflected in higher tumor mutation burden level. HLA-I LOH occurs more frequently in MSS samples than in MSI-H samples. The alteration frequencies of p53 pathway, RTK/RAS pathway, Notch pathway, Hippo pathway, and Nrf2 pathway in HLA-I LOH group were significantly higher than that in HLA-I stable group (p < .0001, p < .0001, p = .032, p = .013, p = .003, respectively). In DNA damage response pathways, alterations in the checkpoint factor pathway and Fanconi anemia pathway are enriched in HLA-I LOH group (p < .0001, p = .023, respectively). Besides, HLA-I LOH was accompanied by higher mutation rates of several tumor suppressors, including TP53 and LRP1B. These results may shed light on follow-up tumor immunology research.

Potential implications of six American Indian patients with myopathy, statin exposure and anti-HMGCR antibodies

OBJECTIVES

Statin-associated autoimmune myopathy is a rare condition associated with the formation of autoantibodies to 3-hydroxy-3-methylglutaryl-coenzyme A reductase. Underlying environmental and genetic risk factors remain poorly understood. American Indians have high rates of cardiovascular disease and associated co-morbidities that require lipid-lowering therapies. We observed this autoimmune myopathy in a series of American Indian statin users in rural Arizona.

METHODS

We reviewed the charts of six American Indian patients with statin-associated autoimmune myopathy. We provide an illustrative case in addition to summaries of clinical presentations and treatment courses.

RESULTS

This is the first report of statin-associated autoimmune myopathy in American Indians. These cases were all identified at the same geographically isolated hospital that exclusively serves an American Indian population with only 1800 statin users. There is relatively low migration. Each case was consistent with the previously described classical presentations for the disease. All six of our cases had diabetes and developed myopathy on high-dose atorvastatin, often with a recent change in statin type or dose.

CONCLUSION

Providers serving American Indians need to be aware of the possibility of statin-associated autoimmune myopathy and familiar with its presentation. Larger, inclusive, population-based investigations are needed to elucidate risk factors for this condition, in particular the potential interactions between predisposing HLA alleles, diabetes and specific statin exposures. This is necessary to identify effective and safe lipid-lowering medications.

Tumor microenvironment characterization in stage IV gastric cancer

Immunotherapy is remarkably affected by the immune environment of the principal tumor. Nonetheless, the immune environment's clinical relevance in stage IV gastric cancer (GC) is largely unknown. The gene expression profiles of 403 stage IV GC patients in the three cohorts: GEO (Gene Expression Omnibus, GSE84437 (n=292) and GSE62254 (n=77), and TCGA (The Cancer Genome Atlas, n=34) were used in the present study. Using four publicly available stage IV GC expression datasets, 29 immune signatures were expression profiled, and on this basis, we classified stage IV GC. The classification was conducted using the hierarchical clustering method. Three stage IV GC subtypes L, M, and H were identified representing low, medium, and high immunity, respectively. Immune correlation analysis of these three types revealed that Immune H exhibited a better prognostic outcome as well as a higher immune score compared with other subtypes. There was a noticeable difference in the three subgroups of HLA genes. Further, on comparing with other subtypes, CD86, CD80, CD274, CTLA4, PDCD1, and PDCD1LG2 had higher expression in the Immunity H subtype. In stage IV GC, potentially positive associations between immune and pathway activities were displayed, due to the enrichment of pathways including TNF signaling, Th-17 cell differentiation, and JAK-STAT signaling pathways in Immunity H vs Immunity L subtypes. External cohorts from TCGA cohort ratified these results. The identification of stage IV GC subtypes has potential clinical implications in stage IV GC treatment.

The influence of evolutionary history on human health and disease

Nearly all genetic variants that influence disease risk have human-specific origins; however, the systems they influence have ancient roots that often trace back to evolutionary events long before the origin of humans. Here, we review how advances in our understanding of the genetic architectures of diseases, recent human evolution and deep evolutionary history can help explain how and why humans in modern environments become ill. Human populations exhibit differences in the prevalence of many common and rare genetic diseases. These differences are largely the result of the diverse environmental, cultural, demographic and genetic histories of modern human populations. Synthesizing our growing knowledge of evolutionary history with genetic medicine, while accounting for environmental and social factors, will help to achieve the promise of personalized genomics and realize the potential hidden in an individual's DNA sequence to guide clinical decisions. In short, precision medicine is fundamentally evolutionary medicine, and integration of evolutionary perspectives into the clinic will support the realization of its full potential.

Using PacBio SMRT data for identification of class I MHC alleles in a wildlife species, Zalophus californianus (California sea lion)

High allelic polymorphism and association with disease susceptibility has made the genes encoding major histocompatibility complex (MHC) antigen presentation molecules in humans, domesticated animals, and wildlife species of wide interest to ecologists, evolutionary biologists, and health specialists. The often multifaceted polygenism and extreme polymorphism of this immunogenetic system have made it especially difficult to characterize in non-model species. Here we compare and contrast the workflows of traditional Sanger sequencing of plasmid-cloned amplicons to Pacific Biosciences SMRT circular consensus sequencing (CCS) in their ability to capture alleles of MHC class I in a wildlife species where characterization of these genes was absent. We assessed two California sea lions (Zalophus californianus), a species suffering from a high prevalence of an aggressive cancer associated with a sexually transmitted gamma herpesvirus. In this pilot study, SMRT CCS proved superior in identifying more alleles from each animal than the more laborious plasmid cloning/Sanger workflow (12:7, 10:7), and no alleles were identified with the cloning/Sanger approach that were not identified by SMRT CCS. We discuss the advantages and disadvantages of each approach including cost, allele rarefaction, and sequence fidelity.

Host genetics and infectious disease: new tools, insights and translational opportunities

Understanding how human genetics influence infectious disease susceptibility offers the opportunity for new insights into pathogenesis, potential drug targets, risk stratification, response to therapy and vaccination. As new infectious diseases continue to emerge, together with growing levels of antimicrobial resistance and an increasing awareness of substantial differences between populations in genetic associations, the need for such work is expanding. In this Review, we illustrate how our understanding of the host–pathogen relationship is advancing through holistic approaches, describing current strategies to investigate the role of host genetic variation in established and emerging infections, including COVID-19, the need for wider application to diverse global populations mirroring the burden of disease, the impact of pathogen and vector genetic diversity and a broad array of immune and inflammation phenotypes that can be mapped as traits in health and disease. Insights from study of inborn errors of immunity and multi-omics profiling together with developments in analytical methods are further advancing our knowledge of this important area.

Infectious diseases are an ever-present global threat. In this Review, Kwok, Mentzer and Knight discuss our latest understanding of how human genetics influence susceptibility to disease. Furthermore, they discuss emerging progress in the interplay between host and pathogen genetics, molecular responses to infection and vaccination, and opportunities to bring these aspects together for rapid responses to emerging diseases such as COVID-19.

A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system

Genotyping complex multigene families in novel systems is particularly challenging. Target primers frequently amplify simultaneously multiple loci leading to high PCR and sequencing artefacts such as chimeras and allele amplification bias. Most genotyping pipelines have been validated in nonmodel systems whereby the real genotype is unknown and the generation of artefacts may be highly repeatable. Further hindering accurate genotyping, the relationship between artefacts and genotype complexity (i.e. number of alleles per genotype) within a PCR remains poorly described. Here, we investigated the latter by experimentally combining multiple known major histocompatibility complex (MHC) haplotypes of a model organism (chicken, Gallus gallus, 43 artificial genotypes with 2-13 alleles per amplicon). In addition to well-defined 'optimal' primers, we simulated a nonmodel species situation by designing 'cross-species' primers based on sequence data from closely related Galliform species. We applied a novel open-source genotyping pipeline (ACACIA; https://gitlab.com/psc_santos/ACACIA), and compared its performance with another, previously published pipeline (AmpliSAS). Allele calling accuracy was higher when using ACACIA (98.5% versus 97% and 77.8% versus 75% for the 'optimal' and 'cross-species' data sets, respectively). Systematic allele dropout of three alleles owing to primer mismatch in the 'cross-species' data set explained high allele calling repeatability (100% when using ACACIA) despite low accuracy, demonstrating that repeatability can be misleading when evaluating genotyping workflows. Genotype complexity was positively associated with nonchimeric artefacts, chimeric artefacts (nonlinearly by levelling when amplifying more than 4-6 alleles) and allele amplification bias. Our study exemplifies and demonstrates pitfalls researchers should avoid to reliably genotype complex multigene families.

Genomics-based immuno-oncology: bridging the gap between immunology and tumor biology

The first hypotheses about how the immune system affects cancers were proposed in the early 20th century. These early concepts about cancer immunosurveillance were further developed in the decades that followed, but a detailed understanding of cancer immunity remained elusive. It was only recently, through the advent of high-throughput technologies, that scientists gained the ability to profile tumors with a resolution that allowed for granular assessment of both tumor cells and the tumor microenvironment. The advent of immune checkpoint inhibitors (ICIs), which have proven to be effective cancer therapies in many malignancies, has spawned great interest in developing biomarkers for efficacy, an endeavor that highlighted the value of dissecting tumor immunity using large-scale methods. Response to ICI therapy has been shown to be a highly complex process, where the dynamics of tumor and immune cells is key to success. The need to understand the biologic mechanisms at the tumor-immune interface has given rise to the field of cancer immunogenomics, a discipline that aims to bridge the gap between cancer genomics and classical immunology. We provide a broad overview of this emerging branch of translational science, summarizing common platforms used and recent discoveries in the field, which are having direct clinical implications. Our discussion will be centered around the genetic foundations governing tumor immunity and molecular determinants associated with clinical benefit from ICI therapy. We emphasize the importance of molecular diversity as a driver of anti-tumor immunity and discuss how these factors can be probed using genomic approaches.

Concanamycin A counteracts HIV-1 Nef to enhance immune clearance of infected primary cells by cytotoxic T lymphocytes

Nef is an HIV-encoded accessory protein that enhances pathogenicity by down-regulating major histocompatibility class I (MHC-I) expression to evade killing by cytotoxic T lymphocytes (CTLs). A potent Nef inhibitor that restores MHC-I is needed to promote immune-mediated clearance of HIV-infected cells. We discovered that the plecomacrolide family of natural products restored MHC-I to the surface of Nef-expressing primary cells with variable potency. Concanamycin A (CMA) counteracted Nef at subnanomolar concentrations that did not interfere with lysosomal acidification or degradation and were nontoxic in primary cell cultures. CMA specifically reversed Nef-mediated down-regulation of MHC-I, but not CD4, and cells treated with CMA showed reduced formation of the Nef:MHC-I:AP-1 complex required for MHC-I down-regulation. CMA restored expression of diverse allotypes of MHC-I in Nef-expressing cells and inhibited Nef alleles from divergent clades of HIV and simian immunodeficiency virus, including from primary patient isolates. Lastly, we found that restoration of MHC-I in HIV-infected cells was accompanied by enhanced CTL-mediated clearance of infected cells comparable to genetic deletion of Nef. Thus, we propose CMA as a lead compound for therapeutic inhibition of Nef to enhance immune-mediated clearance of HIV-infected cells.

Blood Monocyte Subsets with Activation Markers in Relation with Macrophages in Non-Small Cell Lung Cancer

Simple Summary

This study characterized monocyte subtypes: classical, intermediate, and non-classical with the expression of surfaces markers: CD62L, CD11c, CD18, HLA-DR in non-small cell lung cancer patients (NSCLC) compared to healthy controls and correlations between monocyte subtypes and macrophages in the lung cancer microenvironment. We confirmed the presence of various monocyte subtypes in the blood with predominance of classic monocytes and a higher proportion of classical and intermediate monocytes in NSCLC patients than in healthy subjects. Our observation that intermediate monocytes with CD11c+ and HLA-DR+ expression correlation with the amount of macrophages from the lung cancer microenvironment may indicate role of these cells in cancer immunity. A high proportion of monocytes with low expression of CD62L indicates participation of monocytes in attenuation of anticancer response. The detection and monitoring of the presented monocyte subsets in the blood might be a useful test in lung cancer.

Abstract

(1) The cells from the monocyte line play an important role as regulators of cancer development and progression. Monocytes present pro- and anti-tumor immunity and differentiation into macrophages. Macrophages are predominant in the lung cancer environment and could be evaluated by bronchoalveolar lavage fluid (BALF). (2) The aim of the study was analysis of monocytes: classical, intermediate and non-classical with expression of: CD62L, CD11c, CD18, HLA-DR in non-small cell lung cancer (NSCLC) and their correlation with BALF macrophages from lungs with cancer (clBALF) and healthy lungs (hlBALF). (3) A total of 24 patients with NSCLC and 20 healthy donors were investigated. Monocyte subtyping and macrophage counts were performed by flow cytometry. (4) There are three types in peripheral blood (PB): classical monocytes (CD14++CD16-), intermediate (CD14+CD16+) and non-classical (CD14-/+CD16++). We noticed a higher proportion of classical and intermediate monocytes in lung cancer than in healthy donors (76.2 vs. 67.3, and 7.9 vs. 5.2 p < 0.05). We observed a higher proportion of macrophages in clBALF then in hlBALF. A higher CD62L expression on all monocyte subtypes in healthy donors than in study group was found. There were positive correlations between: classical CD11c+, intermediate CD11c+, intermediate HLA-DR+ monocytes in PB with macrophages in clBALF. We did not observe these correlations with macrophages from hlBALF. (5) A predominance of classical and intermediate monocytes in lung cancer and the correlation between intermediate monocytes with CD11c+ and HLA-DR+ and macrophages from the NSCLC milieu support a role of monocyte-line cells in cancer immunity. A high proportion of monocytes with low expression of CD62L indicates the participation of monocytes in attenuation of anticancer response.

Predictive biomarkers for cancer immunotherapy with immune checkpoint inhibitors

Although the clinical development of immune checkpoint inhibitors (ICIs) therapy has ushered in a new era of anti-tumor therapy, with sustained responses and significant survival advantages observed in multiple tumors, most patients do not benefit. Therefore, more and more attention has been paid to the identification and development of predictive biomarkers for the response of ICIs, and more in-depth and comprehensive understanding has been continuously explored in recent years. Predictive markers of ICIs efficacy have been gradually explored from the expression of intermolecular interactions within tumor cells to the expression of various molecules and cells in tumor microenvironment, and been extended to the exploration of circulating and host systemic markers. With the development of high-throughput sequencing and microarray technology, a variety of biomarker strategies have been deeply explored and gradually achieved the process from the identification of single marker to the development of multifactorial synergistic predictive markers. Comprehensive predictive-models developed by integrating different types of data based on different components of tumor-host interactions is the direction of future research and will have a profound impact in the field of precision immuno-oncology. In this review, we deeply analyze the exploration course and research progress of predictive biomarkers as an adjunctive tool to tumor immunotherapy in effectively identifying the efficacy of ICIs, and discuss their future directions in achieving precision immuno-oncology.

HLA Class I Allele Expression and Clinical Outcome in De Novo Metastatic Prostate Cancer

The prognostic value of human leukocyte antigen (HLA) class I molecules in prostate cancer (PCa) remains unclear. Herein, we investigated the prognostic relevance of the most frequently expressed HLA-A alleles in Greece (A*02:01 and HLA-A*24:02) in de novo metastatic hormone-sensitive PCa (mPCa), which is a rare and aggressive disease characterized by a rapid progression to castration-resistance (CR) and poor overall survival (OS), contributing to almost 50% of PCa-related deaths. We identified 56 patients who had either progressed to CR (these patients were retrospectively analyzed for the time to the progression of CR and prospectively for OS) or had at least three months’ follow-up postdiagnosis without CR progression and, thus, were prospectively analyzed for both CR and OS. Patients expressing HLA-A*02:01 showed poor clinical outcomes vs. HLA-A*02:01⁻negative patients. HLA-A*24:02⁻positive patients progressed slower to CR and had increased OS. Homozygous HLA-A*02:01 patients progressed severely to CR, with very short OS. Multivariate analyses ascribed to both HLA alleles significant prognostic values for the time to progression (TTP) to CR and OS. The presence of HLA-A*02:01 and HLA-A*24:02 alleles in de novo mPCa patients are significantly and independently associated with unfavorable or favorable clinical outcomes, respectively, suggesting their possible prognostic relevance for treatment decision-making in the context of precision medicine.

Hepatitis B infection outcome is associated with novel human leukocyte antigen variants in Ghanaian cohort

IMPACT STATEMENT

Genetic association studies can determine the effect size of gene loci on disease outcomes. In the arena of HBV infections, HLA alleles that associate with HBV outcomes can be used in clinical management decisions. This potential translational utility can shape the future management of HBV infections by identifying at-risk individuals and tailoring medical interventions accordingly. This precision medicine motif is currently only a nascent idea. However, it has stakes that may well override the current "wait and see" approach of clinical management of HBV infections. Here, we have identified HLA alleles associated with HBV outcome in a Ghanaian cohort. Our findings support the motif that HLA alleles associate with HBV outcome along geo-ethnic lines. This buttresses the need for further population pivoted studies. In the long term, our findings add to efforts towards the development of an HLA molecular-based algorithm for predicting HBV infection outcomes.

Host-parasite interactions in non-native invasive species are dependent on the levels of standing genetic variation at the immune locus

Background

Parasites may mediate the success of biological invasions through their effect on host fitness and thus, on host population growth and stability. However, a release from the pressure of parasites is strongly related to the genetic differentiation of the host. In invasive host populations, the number of available genetic variants, allowing them to ‘fight’ the infection, are likely to be influenced by founder events and genetic drift. The level standing genetic variation of invasive populations may be crucial in successfully adapting to new environments and resisting diseases. We studied invasive populations of raccoon that experienced a random reduction in genetic diversity during the establishment and evaluated the relationship between host immune genetic diversity and intestinal parasites infection.

Results

We distinguished two different genetic clusters that are characterized by different sets of functionally relevant MHC-DRB alleles. Both clusters were characterized by considerably different allele-parasite associations and different levels of parasite infection. The specific resistance MHC-DRB alleles explained the lower prevalence of Digenea parasites. An increased infection intensity was related to the presence of two MHC-DRB alleles. One of these alleles significantly decreased in frequency over time, causing a decrease of Digenea abundance in raccoons in consecutive years.

Conclusions

Our findings suggest that intestinal parasites can exert selective pressure on an invasive host with lowered levels of immune genetic diversity and contribute to promoting local adaptation over time. The random genetic drift that created the two different genetic clusters in the invasive raccoon range imposed completely different MHC-parasite associations, strongly associated with the infection status of populations. Our findings underline the role of standing genetic variation in shaping host-parasite relationships and provide empirical support that functional genetic variation may be, at least partly, responsible for differences in the success of invasive populations.

Evolution of major histocompatibility complex class I genes in the sable Martes zibellina (Carnivora, Mustelidae)

The molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immune response among vertebrates. We investigated the molecular evolution of MHC class I genes in the sable Martes zibellina. We isolated 26 MHC class I sequences, including 12 putatively functional sequences and 14 pseudogene sequences, from 24 individuals from two geographic areas of northeast China. The number of putatively functional sequences found in a single individual ranged from one to five, which might be at least 1-3 loci. We found that both balancing selection and recombination contribute to evolution of MHC class I genes in M. zibellina. In addition, we identified a candidate nonclassical MHC class I lineage in Carnivora, which may have preceded the divergence (about 52-57 Mya) of Caniformia and Feliformia. This may contribute to further understanding of the origin and evolution of nonclassical MHC class I genes. Our study provides important immune information of MHC for M. zibellina, as well as other carnivores.

Differential prognostic impact of CD8+ T cells based on human leucocyte antigen I and PD-L1 expression in microsatellite-unstable gastric cancer

BACKGROUND

The aim of the study was to determine the human leucocyte antigen class-I (HLA-I), programmed death-ligand 1 (PD-L1) expression and tumour-infiltrating lymphocytes (TILs) of microsatellite instability-high gastric cancer.

METHODS

The HLA-I expression type was determined by immunohistochemistry of HLA-A, HLA-B, HLA-C and β2-microglobulin in the centre of the tumour (CT) and in the invasive margin (IM) of samples from 293 patients (total loss vs. preserved type). PD-L1 expression and TIL density was examined immunohistochemically. HLA-I genotyping was also performed.

RESULTS

The expression loss of the HLA-I molecules was significantly associated with low TIL density. According to survival analyses, the HLA-I expression type and PD-L1 positivity were not independent prognostic factors. The TIL density had no prognostic implication when survival analysis was performed for the whole patient group; however, high CD8⁺ TIL infiltration was significantly associated with good prognosis in only HLA-I-preserved-type/PD-L1-positive group (p = 0.034). The homozygosity of the HLA-I allele was more frequently observed in the total loss type group.

CONCLUSIONS

We confirmed differential prognostic implication of CD8⁺ TILs according to the HLA-I and PD-L1 expression. Determination of the HLA-I expression could be helpful to select patients who would benefit from anti-PD-1/PD-L1 therapy.

Genomic evidence for MHC disassortative mating in humans

Although pervasive in many animal species, the evidence for major histocompatibility complex (MHC) disassortative mating in humans remains inconsistent across studies. Here, to revisit this issue, we analyse dense genotype data for 883 European and Middle Eastern couples. To distinguish MHC-specific effects from socio-cultural confounders, the pattern of relatedness between spouses in the MHC region is compared to the rest of the genome. Couples from Israel exhibit no significant pattern of relatedness across the MHC region, whereas across the genome, they are more similar than random pairs of individuals, which may reflect social homogamy and/or cousin marriages. On the other hand, couples from The Netherlands and more generally from Northern Europe are significantly more MHC-dissimilar than random pairs of individuals, and this pattern of dissimilarity is extreme when compared with the rest of the genome. Our findings support the hypothesis that the MHC influences mate choice in humans in a context-dependent way: MHC-driven preferences may exist in all populations but, in some populations, social constraints over mate choice may reduce the ability of individuals to rely on such biological cues when choosing their mates.

Association of HLA Antigen Mismatch With Risk of Developing Skin Cancer After Solid-Organ Transplant

Importance

Risk factors for the development of skin cancer after solid-organ transplant can inform clinical care, but data on these risk factors are limited.

Objective

To study the association between HLA antigen mismatch and skin cancer incidence after solid-organ transplant.

Design, Setting, and Participants

This retrospective cohort study is a secondary analysis of the multicenter Transplant Skin Cancer Network study of 10 649 adults who underwent a primary solid-organ transplant between January 1, 2003, and December 31, 2003, or between January 1, 2008, and December 31, 2008. These participants were identified through the Scientific Registry of Transplant Recipients standard analysis files, which contain data collected mostly by the Organ Procurement and Transplantation Network. Participants were matched to skin cancer outcomes by medical record review. This study was conducted from August 1, 2016, to July 31, 2017.

Main Outcomes and Measures

The primary outcome was time to diagnosis of posttransplant skin cancer, including squamous cell carcinoma, melanoma, and Merkel cell carcinoma. The HLA antigen mismatch was calculated based on the 2016 Organ Procurement and Transplantation Network guidelines. Risk of skin cancer was analyzed using a multivariate Cox proportional hazards regression model.

Results

In total, 10 649 organ transplant recipients (6776 men [63.6%], with a mean [SD] age of 51 [12] years) contributed 59 923 years of follow-up. For each additional mismatched allele, a 7% to 8% reduction in skin cancer risk was found (adjusted hazard ratio [HR], 0.93; 95% CI, 0.87-0.99; P = .01). Subgroup analysis found the protective effect of HLA antigen mismatch to be statistically significant in lung (adjusted HR, 0.70; 95% CI, 0.56-0.87; P = .001) and heart (adjusted HR, 0.75; 95% CI, 0.60-0.93; P = .008) transplant recipients but not for recipients of liver, kidney, or pancreas. The degree of HLA-DR mismatch, but not HLA-A or HLA-B mismatch, was the most statistically significant for skin cancer risk (adjusted HR, 0.85; 95% CI, 0.74-0.97; P = .01).

Conclusions and Relevance

The HLA antigen mismatch appears to be associated with reductions in the risk of skin cancer after solid-organ transplant among heart and lung transplant recipients; this finding suggests that HLA antigen mismatch activates the tumor surveillance mechanisms that protect against skin cancer in transplant recipients and that skin cancer risk may be higher in patients who received a well-matched organ.

Regulatory T cells and co-evolution of allele-specific MHC recognition by the TCR

What is the evolutionary mechanism for the TCR-MHC-conserved interaction? We extend Dembic's model (Dembic Z. In, Scand J Immunol e12806, 2019) of thymus positive selection for high-avidity anti-self-MHC Tregs among double (CD4 + CD8+)-positive (DP) developing thymocytes. This model is based on competition for self-MHC (+ Pep) complexes presented on cortical epithelium. Such T cells exit as CD4 + CD25+FoxP3 + thymic-derived Tregs (tTregs). The other positively selected DP T cells are then negatively selected on medulla epithelium removing high-avidity anti-self-MHC + Pep as T cells commit to CD4 + or CD8 + lineages. The process is likened to the competitive selection and affinity maturation in Germinal Centre for the somatic hypermutation (SHM) of rearranged immunoglobulin (Ig) variable region (V[D]Js) of centrocytes bearing antigen-specific B cell receptors (BCR). We now argue that the same direct SHM processes for TCRs occur in post-antigenic Germinal Centres, but now occurring in peripheral pTregs. This model provides a potential solution to a long-standing problem previously recognized by Cohn and others (Cohn M, Anderson CC, Dembic Z. In, Scand J Immunol e12790, 2019) of how co-evolution occurs of species-specific MHC alleles with the repertoire of their germline TCR V counterparts. We suggest this is not by 'blind', slow, and random Darwinian natural selection events, but a rapid structured somatic selection vertical transmission process. The pTregs bearing somatic TCR V mutant genes then, on arrival in reproductive tissues, can donate their TCR V sequences via soma-to-germline feedback as discussed in this journal earlier. (Steele EJ, Lindley RA. In, Scand J Immunol e12670, 2018) The high-avidity tTregs also participate in the same process to maintain a biased, high-avidity anti-self-MHC germline V repertoire.

Single haplotype admixture models using large scale HLA genotype frequencies to reproduce human admixture

The human leukocyte antigen (HLA) is the most polymorphic region in humans. Anthropologists use HLA to trace populations' migration and evolution. However, recent admixture between populations can mask the ancestral haplotype frequency distribution. We present a statistical method based on high-resolution HLA haplotype frequencies to resolve population admixture using a non-negative matrix factorization formalism and validated using haplotype frequencies from 56 world populations. The result is a minimal set of source components (SCs) decoding roughly 90% of the total variance in the studied admixtures. These SCs agree with the geographical distribution, phylogenies, and recent admixture events of the studied groups. With the growing population of multi-ethnic individuals, or individuals that do not report race/ethnic information, the HLA matching process for stem-cell and solid organ transplants is becoming more challenging. The presented algorithm provides a framework that facilitates the breakdown of highly admixed populations into SCs, which can be used to better match the rapidly growing population of multi-ethnic individuals worldwide.