The IPD Project: a centralised resource for the study of polymorphism in genes of the immune system

Genetic Diversity and Ethnic Tapestry of Kazakhstan as Inferred from HLA Polymorphism and Population Dynamics: A Comprehensive Review

Background: The human leukocyte antigen (HLA) system represents the most polymorphic segment within human DNA sequences and constitutes a core component of immune defense responses and in understanding population genetics. This research investigates the distribution of HLA class I and II polymorphisms across different ethnic groups in Kazakhstan, offering valuable insights into the genetic diversity and demographic evolution within this region. Methods: We performed an in-depth examination of HLA class I and II polymorphisms across diverse ethnic communities living in Kazakhstan, including Kazakhs, Russians, Uzbeks, Ukrainians, Germans, Tatars, and Koreans. Utilizing data from high-resolution HLA typing studies allowed us to assess allele frequencies alongside haplotype distributions while analyzing genetic interrelations between these populations. Additionally, we performed comparative assessments with global HLA databases to determine the genetic affiliations between these groups and their relationships with neighboring and more distant populations. Results: Our study revealed over 200 HLA alleles within the analyzed populations, and significant variations were observed in their allele and haplotype frequencies. Notably, the Kazakh group exhibited strong genetic ties to Asian and Siberian demographics; conversely, other ethnicities showed associations reflective of their historical roots. Notable alleles included HLA-A*02:01, B*07:02, C*07:02, DRB1*07:01, and DQB1*03:01, commonly observed across various groups. Linkage disequilibrium analysis revealed the presence of population-specific haplotypes, highlighting distinct genetic structures within these communities. Conclusions: The findings highlight the significant genetic diversity in Kazakhstan, influenced by its geographical location at the crossroads of Europe and Asia. These results are pertinent to immunogenetics, transplantation medicine, and personalized healthcare within Kazakhstan and adjacent regions. Future research should expand the sample size and explore disease associations to enhance our comprehension of HLA genetics across Central Asia.

The 2024 IPD-MHC database update: a comprehensive resource for major histocompatibility complex studies

The IPD-MHC Database project (http://www.ebi.ac.uk/ipd/mhc/) serves as a comprehensive and expertly curated repository for major histocompatibility complex (MHC) sequences from non-human species, providing the necessary infrastructure and tools to study the function and evolution of this highly polymorphic genomic region. In its latest version, the IPD-MHC database has expanded both in content and in the tools for data visualization and comparison. The database now hosts over 18 000 MHC alleles from 125 species, organized into eleven taxonomic groups, all manually curated and named by the Comparative MHC Nomenclature Committee. A cetacean section has recently been included, offering researchers valuable data to study the immune system of whales, dolphins, and porpoises, as well establishing the official nomenclature platform for the Cetacea Leukocyte Antigens (CeLA). In response to user demand and reflecting broader trends in bioinformatics and immunogenetics, IPD-MHC now includes the predicted tertiary structure of over 8000 alleles and allows comparison and visualisation of allele variation within and between species at single residue resolution. These latest developments maintain the critically important link between official nomenclature of curated alleles and the ability to analyse this complex polymorphism using the most up to date methods within a single repository.

High KIR diversity in Uganda and Botswana children living with HIV

Killer-cell immunoglobulin-like receptors (KIRs) are essential components of the innate immune system found on the surfaces of natural killer (NK) cells. The KIRs encoding genes are located on chromosome 19q13.4 and are genetically diverse across populations. KIRs are associated with various disease states including HIV progression, and are linked to transplantation rejection and reproductive success. However, there is limited knowledge on the diversity of KIRs from Uganda and Botswana HIV-infected paediatric cohorts, with high endemic HIV rates. We used next-generation sequencing technologies on 312 (246 Uganda, 66 Botswana) samples to generate KIR allele data and employed customised bioinformatics techniques for allelic, allotype and disease association analysis. We show that these sample sets from Botswana and Uganda have different KIRs of different diversities. In Uganda, we observed 147 vs 111 alleles in the Botswana cohort, which had a more than 1 % frequency. We also found significant deviation towards homozygosity for the KIR3DL2 gene for both rapid (RPs) and long-term non-progressors (LTNPs)in the Ugandan cohort. The frequency of the bw4-80I ligand was also significantly higher among the LTNPs than RPs (8.9 % Vs 2.0%, P-value: 0.032). In the Ugandan cohort, KIR2DS4*001 (OR: 0.671, 95 % CI: 0.481-0.937, FDR adjusted Pc=0.142) and KIR2DS4*006 (OR: 2.519, 95 % CI: 1.085-5.851, FDR adjusted Pc=0.142) were not associated with HIV disease progression after adjustment for multiple testing. Our study results provide additional knowledge of the genetic diversity of KIRs in African populations and provide evidence that will inform future immunogenetics studies concerning human disease susceptibility, evolution and host immune responses.

A genome assembly and transcriptome atlas of the inbred Babraham pig to illuminate porcine immunogenetic variation

The inbred Babraham pig serves as a valuable biomedical model for research due to its high level of homozygosity, including in the major histocompatibility complex (MHC) loci and likely other important immune-related gene complexes, which are generally highly diverse in outbred populations. As the ability to control for this diversity using inbred organisms is of great utility, we sought to improve this resource by generating a long-read whole genome assembly and transcriptome atlas of a Babraham pig. The genome was de novo assembled using PacBio long reads and error-corrected using Illumina short reads. Assembled contigs were then mapped to the porcine reference assembly, Sscrofa11.1, to generate chromosome-level scaffolds. The resulting TPI_Babraham_pig_v1 assembly is nearly as contiguous as Sscrofa11.1 with a contig N50 of 34.95 Mb and contig L50 of 23. The remaining sequence gaps are generally the result of poor assembly across large and highly repetitive regions such as the centromeres and tandemly duplicated gene families, including immune-related gene complexes, that often vary in gene content between haplotypes. We also further confirm homozygosity across the Babraham MHC and characterize the allele content and tissue expression of several other immune-related gene complexes, including the antibody and T cell receptor loci, the natural killer complex, and the leukocyte receptor complex. The Babraham pig genome assembly provides an alternate highly contiguous porcine genome assembly as a resource for the livestock genomics community. The assembly will also aid biomedical and veterinary research that utilizes this animal model such as when controlling for genetic variation is critical.

Genetic complexity of killer-cell immunoglobulin-like receptor genes in human pangenome assemblies

The killer-cell immunoglobulin-like receptor (KIR) gene complex, a highly polymorphic region of the human genome that encodes proteins involved in immune responses, poses strong challenges in genotyping owing to its remarkable genetic diversity and structural intricacy. Accurate analysis of KIR alleles, including their structural variations, is crucial for understanding their roles in various immune responses. Leveraging the high-quality genome assemblies from the Human Pangenome Reference Consortium (HPRC), we present a novel bioinformatic tool, the structural KIR annoTator (SKIRT), to investigate gene diversity and facilitate precise KIR allele analysis. In 47 HPRC-phased assemblies, SKIRT identifies a recurrent novel KIR2DS4/3DL1 fusion gene in the paternal haplotype of HG02630 and maternal haplotype of NA19240. Additionally, SKIRT accurately identifies eight structural variants and 15 novel nonsynonymous alleles, all of which are independently validated using short-read data or quantitative polymerase chain reaction. Our study has discovered a total of 570 novel alleles, among which eight haplotypes harbor at least one KIR gene duplication, six haplotypes have lost at least one framework gene, and 75 out of 94 haplotypes (79.8%) carry at least five novel alleles, thus confirming KIR genetic diversity. These findings are pivotal in providing insights into KIR gene diversity and serve as a solid foundation for understanding the functional consequences of KIR structural variations. High-resolution genome assemblies offer unprecedented opportunities to explore polymorphic regions that are challenging to investigate using short-read sequencing methods. The SKIRT pipeline emerges as a highly efficient tool, enabling the comprehensive detection of the complete spectrum of KIR alleles within human genome assemblies.

The Primate Major Histocompatibility Complex: An Illustrative Example of Gene Family Evolution

Gene families are groups of evolutionarily-related genes. One large gene family that has experienced rapid evolution is the Major Histocompatibility Complex (MHC), whose proteins serve critical roles in innate and adaptive immunity. Across the ~60 million year history of the primates, some MHC genes have turned over completely, some have changed function, some have converged in function, and others have remained essentially unchanged. Past work has typically focused on identifying MHC alleles within particular species or comparing gene content, but more work is needed to understand the overall evolution of the gene family across species. Thus, despite the immunologic importance of the MHC and its peculiar evolutionary history, we lack a complete picture of MHC evolution in the primates. We readdress this question using sequences from dozens of MHC genes and pseudogenes spanning the entire primate order, building a comprehensive set of gene and allele trees with modern methods. Overall, we find that the Class I gene subfamily is evolving much more quickly than the Class II gene subfamily, with the exception of the Class II MHC-DRB genes. We also pay special attention to the often-ignored pseudogenes, which we use to reconstruct different events in the evolution of the Class I region. We find that despite the shared function of the MHC across species, different species employ different genes, haplotypes, and patterns of variation to achieve a successful immune response. Our trees and extensive literature review represent the most comprehensive look into MHC evolution to date.

A pan-cetacean MHC amplicon sequencing panel developed and evaluated in combination with genome assemblies

The major histocompatibility complex (MHC) is a highly polymorphic gene family that is crucial in immunity, and its diversity can be effectively used as a fitness marker for populations. Despite this, MHC remains poorly characterised in non-model species (e.g., cetaceans: whales, dolphins and porpoises) as high gene copy number variation, especially in the fast-evolving class I region, makes analyses of genomic sequences difficult. To date, only small sections of class I and IIa genes have been used to assess functional diversity in cetacean populations. Here, we undertook a systematic characterisation of the MHC class I and IIa regions in available cetacean genomes. We extracted full-length gene sequences to design pan-cetacean primers that amplified the complete exon 2 from MHC class I and IIa genes in one combined sequencing panel. We validated this panel in 19 cetacean species and described 354 alleles for both classes. Furthermore, we identified likely assembly artefacts for many MHC class I assemblies based on the presence of class I genes in the amplicon data compared to missing genes from genomes. Finally, we investigated MHC diversity using the panel in 25 humpback and 30 southern right whales, including four paternity trios for humpback whales. This revealed copy-number variable class I haplotypes in humpback whales, which is likely a common phenomenon across cetaceans. These MHC alleles will form the basis for a cetacean branch of the Immuno-Polymorphism Database (IPD-MHC), a curated resource intended to aid in the systematic compilation of MHC alleles across several species, to support conservation initiatives.

Donor KIR2DL1 Allelic Polymorphism Influences Posthematopoietic Progenitor Cell Transplantation Outcomes in the T Cell Depleted and Reduced Intensity Conditioning Setting

The majority of established KIR clinical assessment algorithms used for donor selection for hematopoietic progenitor cell transplantation (HPCT) evaluate gene content (presence/absence) of the KIR gene complex. In comparison, relatively little is known about the impact of KIR allelic polymorphism. By analyzing donors of T cell depleted (TcD) reduced intensity conditioning (RIC) HPCT, this study investigated the influence on post-transplant outcome of 2 polymorphic residues of the inhibitory KIR2DL1. The aim of this study was to expand upon existing research into the influence of KIR2DL1 allelic polymorphism upon post-transplant outcome. The effects of allele groups upon transplant outcomes were investigated within a patient cohort using a defined treatment protocol of RIC with TcD. Using phylogenetic data, KIR2DL1 allelic polymorphism was categorized into groups on the basis of variation within codons 114 and 245 (positive or negative for the following groups: KIR2DL1*002/001g, KIR2DL1*003, KIR2DL1*004g) and the identification of null alleles. The influence of these KIR2DL1 allele groups in hematopoietic progenitor cell transplantation (HPCT) donors was assessed in the post-transplant data of 86 acute myelogenous leukemia patients receiving RIC TcD HPCT at a single center. KIR2DL1 allele groups in the donor significantly impacted upon 5-year post-transplant outcomes in RIC TcD HPCT. Donor KIR2DL1*003 presented the greatest influence upon post-transplant outcomes, with KIR2DL1*003 positive donors severely reducing 5-year post-transplant overall survival (OS) compared to those receiving a transplant from a KIR2DL1*003 negative donor (KIR2DL1*003 pos versus neg: 27.0% versus 60.0%, P = .008, pc = 0.024) and disease-free survival (DFS) (KIR2DL1*003 pos versus neg: 23.5% versus 60.0%, P = .004, pc = 0.012), and increasing 5-year relapse incidence (KIR2DL1*003 pos versus neg: 63.9% versus 27.2%, P = .009, pc = 0.027). KIR2DL1*003 homozygous and KIR2DL1*003 heterozygous grafts did not present significantly different post-transplant outcomes. Donors possessing the KIR2DL1*002/001 allele group were found to significantly improve post-transplant outcomes, with donors positive for the KIR2DL1*004 allele group presenting a trend towards improvement. KIR2DL1*002/001 allele group (KIR2DL1*002/001g) positive donors improved 5-year OS (KIR2DL1*002/001g pos versus neg: 56.4% versus 27.2%, P = .009, pc = 0.024) and DFS (KIR2DL1*002/001g pos versus neg: 53.8% versus 25.5%, P = .018, pc = 0.036). KIR2DL1*004 allele group (KIR2DL1*004g) positive donors trended towards improving 5-year OS (KIR2DL1*004g pos versus neg: 53.3% versus 35.5%, P = .097, pc = 0.097) and DFS (KIR2DL1*004g pos versus neg: 50.0% versus 33.9%, P = .121, pc = 0.121), and reducing relapse incidence (KIR2DL1*004g pos versus neg: 33.1% versus 54.0%, P = .079, pc = 0.152). The presented findings suggest donor selection algorithms for TcD RIC HPCT should consider avoiding KIR2DL1*003 positive donors, where possible, and contributes to the mounting evidence that KIR assessment in donor selection algorithms should reflect the conditioning regime protocol used.

Exceptional diversity of KIR and HLA class I in Egypt

Genetically determined variation of killer cell immunoglobulin like receptors (KIR) and their HLA class I ligands affects multiple aspects of human health. Their extreme diversity is generated through complex interplay of natural selection for pathogen resistance and reproductive health, combined with demographic structure and dispersal. Despite significant importance to multiple health conditions of differential effect across populations, the nature and extent of immunogenetic diversity is under-studied for many geographic regions. Here, we describe the first high-resolution analysis of KIR and HLA class I combinatorial diversity in Northern Africa. Analysis of 125 healthy unrelated individuals from Cairo in Egypt yielded 186 KIR alleles arranged in 146 distinct centromeric and 79 distinct telomeric haplotypes. The most frequent haplotypes observed were KIR-A, encoding two inhibitory receptors specific for HLA-C, two that are specific for HLA-A and -B, and no activating receptors. Together with 141 alleles of HLA class I, 75 of which encode a KIR ligand, we identified a mean of six distinct interacting pairs of inhibitory KIR and HLA allotypes per individual. We additionally characterize 16 KIR alleles newly identified in the study population. Our findings place Egyptians as one of the most highly diverse populations worldwide, with important implications for transplant matching and studies of immune-mediated diseases.

Using a novel structure/function approach to select diverse swine major histocompatibility complex 1 alleles to predict epitopes for vaccine development

MOTIVATION

Swine leukocyte antigens (SLAs) (i.e. swine major histocompatibility complex proteins) conduct a fundamental role in swine immunity. To generate a protective vaccine across an outbred species, such as pigs, it is critical that epitopes that bind to diverse SLA alleles are used in the vaccine development process. We introduced a new strategy for epitope prediction.

RESULTS

We employed molecular dynamics simulation to identify key amino acids for interactions with epitopes. We developed an algorithm wherein each SLA-1 is compared to a crystalized reference allele with unique weighting for non-conserved amino acids based on R group and position. We then performed homology modeling and electrostatic contact mapping to visualize how relatively small changes in sequences impacted the charge distribution in the binding site. We selected eight diverse SLA-1 alleles and performed homology modeling followed, by protein-peptide docking and binding affinity analyses, to identify porcine reproductive and respiratory syndrome virus matrix protein epitopes that bind with high affinity to these alleles. We also performed docking analysis on the epitopes identified as strong binders using NetMHCpan 4.1. Epitopes predicted to bind to our eight SLA-1 alleles had equivalent or higher energetic interactions than those predicted to bind to the NetMHCpan 4.1 allele repertoire. This approach of selecting diverse SLA-1 alleles, followed by homology modeling, and docking simulations, can be used as a novel strategy for epitope prediction that complements other available tools and is especially useful when available tools do not offer a prediction for SLAs/major histocompatibility complex.

AVAILABILITY AND IMPLEMENTATION

The data underlying this article are available in the online Supplementary Material.

Allelic variation of KIR and HLA tunes the cytolytic payload and determines functional hierarchy of NK cell repertoires

The functionality of natural killer (NK) cells is tuned during education and is associated with remodeling of the lysosomal compartment. We hypothesized that genetic variation in killer cell immunoglobulin-like receptor (KIR) and HLA, which is known to influence the functional strength of NK cells, fine-tunes the payload of effector molecules stored in secretory lysosomes. To address this possibility, we performed a high-resolution analysis of KIR and HLA class I genes in 365 blood donors and linked genotypes to granzyme B loading and functional phenotypes. We found that granzyme B levels varied across individuals but were stable over time in each individual and genetically determined by allelic variation in HLA class I genes. A broad mapping of surface receptors and lysosomal effector molecules revealed that DNAM-1 and granzyme B levels served as robust metric of the functional state in NK cells. Variation in granzyme B levels at rest was tightly linked to the lytic hit and downstream killing of major histocompatibility complex-deficient target cells. Together, these data provide insights into how variation in genetically hardwired receptor pairs tunes the releasable granzyme B pool in NK cells, resulting in predictable hierarchies in global NK cell function.

Innate receptors modulating adaptive T cell responses: KIR-HLA interactions and T cell-mediated control of chronic viral infections

Killer-cell immunoglobulin-like receptors (KIRs) are mainly expressed on natural killer (NK) cells and are key regulators of innate immune responses. NK cells are the first responders in the face of infection and help promote placentation during pregnancy; the importance of KIRs in these NK-mediated processes is well-established. However, mounting evidence suggests that KIRs also have a prominent and long-lasting effect on the adaptive immune system. Here, we review the evidence for the impact of KIRs on T cell responses with a focus on the clinical significance of this interaction.

The VGNC: expanding standardized vertebrate gene nomenclature

The Vertebrate Gene Nomenclature Committee (VGNC) was established in 2016 as a sister project to the HUGO Gene Nomenclature Committee, to approve gene nomenclature in vertebrate species without an existing dedicated nomenclature committee. The VGNC aims to harmonize gene nomenclature across selected vertebrate species in line with human gene nomenclature, with orthologs assigned the same nomenclature where possible. This article presents an overview of the VGNC project and discussion of key findings resulting from this work to date. VGNC-approved nomenclature is accessible at https://vertebrate.genenames.org and is additionally displayed by the NCBI, Ensembl, and UniProt databases.

Identification of a DRB3*011:01-restricted CD4+ T cell response against bovine respiratory syncytial virus fusion protein

Although Human Respiratory Syncytial Virus (HRSV) is a significant cause of severe respiratory disease with high morbidity and mortality in pediatric and elderly populations worldwide there is no licensed vaccine. Bovine Respiratory Syncytial Virus (BRSV) is a closely related orthopneumovirus with similar genome structure and high homology between structural and nonstructural proteins. Like HRSV in children, BRSV is highly prevalent in dairy and beef calves and known to be involved in the etiology of bovine respiratory disease, in addition to being considered an excellent model for HRSV. Commercial vaccines are currently available for BRSV, though improvements in efficacy are needed. The aims of this study were to identify CD4+ T cell epitopes present in the fusion glycoprotein of BRSV, an immunogenic surface glycoprotein that mediates membrane fusion and a major target of neutralizing antibodies. Overlapping peptides representing three regions of the BRSV F protein were used to stimulate autologous CD4+ T cells in ELISpot assays. T cell activation was observed only in cells from cattle with the DRB3*011:01 allele by peptides from AA249-296 of the BRSV F protein. Antigen presentation studies with C-terminal truncated peptides further defined the minimum peptide recognized by the DRB3*011:01 allele. Computationally predicted peptides presented by artificial antigen presenting cells further confirmed the amino acid sequence of a DRB3*011:01 restricted class II epitope on the BRSV F protein. These studies are the first to identify the minimum peptide length of a BoLA-DRB3 class II-restricted epitope in BRSV F protein.

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.

Enhanced Extracellular Transfer of HLA-DQ Activates CD3+ Lymphocytes towards Compromised Treg Induction in Celiac Disease

Celiac disease (CeD) manifests with autoimmune intestinal inflammation from gluten and genetic predisposition linked to human leukocyte antigen class-II (HLA-II) gene variants. Antigen-presenting cells facilitate gluten exposition through the interaction of their surface major histocompatibility complex (MHC) with the T cell receptor (TCR) on T lymphocytes. This fundamental mechanism of adaptive immunity has broadened upon recognition of extracellular exosomal MHC, raising awareness of an alternative means for antigen presentation. This study demonstrates that conditioned growth media (CGM) previously exposed to monocyte-derived dendritic cells from CeD significantly downregulates the CD3⁺ lineage marker of control T cells. Such increased activation was reflected in their elevated IL-2 secretion. Exosome localization motif identification and quantification within HLA-DQA1 and HLA-DQB1 transcripts highlighted their significant prevalence within HLA-DQB1 alleles associated with CeD susceptibility. Flow cytometry revealed the strong correlation between HLA-DQ and the CD63 exosomal marker in T cells exposed to CGM from MoDCs sourced from CeD patients. This resulted in lower concentrations of CD25⁺ CD127⁻ T cells, suggestive of their compromised induction to T-regulatory cells associated with CeD homeostasis. This foremost comparative study deciphered the genomic basis and extracellular exosomal effects of HLA transfer on T lymphocytes in the context of CeD, offering greater insight into this auto-immune disease.

Highly-contiguous bovine genomes underpin accurate functional analyses and updated nomenclature of MHC class I

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3' end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3' coding sequence, we identified distinct subgroups of BoLA-3 and BoLA-6 that represent distinct genetic loci. We further examined allele-specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle.

Rapid Characterization of Complex Killer Cell Immunoglobulin-Like Receptor (KIR) Regions Using Cas9 Enrichment and Nanopore Sequencing

Long-read sequencing approaches have considerably improved the quality and contiguity of genome assemblies. Such platforms bear the potential to resolve even extremely complex regions, such as multigenic immune families and repetitive stretches of DNA. Deep sequencing coverage, however, is required to overcome low nucleotide accuracy, especially in regions with high homopolymer density, copy number variation, and sequence similarity, such as the MHC and KIR gene clusters of the immune system. Therefore, we have adapted a targeted enrichment protocol in combination with long-read sequencing to efficiently annotate complex KIR gene regions. Using Cas9 endonuclease activity, segments of the KIR gene cluster were enriched and sequenced on an Oxford Nanopore Technologies platform. This provided sufficient coverage to accurately resolve and phase highly complex KIR haplotypes. Our strategy eliminates PCR-induced amplification errors, facilitates rapid characterization of large and complex multigenic regions, including its epigenetic footprint, and is applicable in multiple species, even in the absence of a reference genome.

DR2S: an integrated algorithm providing reference-grade haplotype sequences from heterozygous samples

Background

High resolution HLA genotyping of donors and recipients is a crucially important prerequisite for haematopoetic stem-cell transplantation and relies heavily on the quality and completeness of immunogenetic reference sequence databases of allelic variation.

Results

Here, we report on DR2S, an R package that leverages the strengths of two sequencing technologies—the accuracy of next-generation sequencing with the read length of third-generation sequencing technologies like PacBio’s SMRT sequencing or ONT’s nanopore sequencing—to reconstruct fully-phased high-quality full-length haplotype sequences. Although optimised for HLA and KIR genes, DR2S is applicable to all loci with known reference sequences provided that full-length sequencing data is available for analysis. In addition, DR2S integrates supporting tools for easy visualisation and quality control of the reconstructed haplotype to ensure suitability for submission to public allele databases.

Conclusions

DR2S is a largely automated workflow designed to create high-quality fully-phased reference allele sequences for highly polymorphic gene regions such as HLA or KIR. It has been used by biologists to successfully characterise and submit more than 500 HLA alleles and more than 500 KIR alleles to the IPD-IMGT/HLA and IPD-KIR databases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04153-0.

Exploiting Tumor Neoantigens to Target Cancer Evolution: Current Challenges and Promising Therapeutic Approaches

Immunotherapeutic manipulation of the antitumor immune response offers an attractive strategy to target genomic instability in cancer. A subset of tumor-specific somatic mutations can be translated into immunogenic and HLA-bound epitopes called neoantigens, which can induce the activation of helper and cytotoxic T lymphocytes. However, cancer immunoediting and immunosuppressive mechanisms often allow tumors to evade immune recognition. Recent evidence also suggests that the tumor neoantigen landscape extends beyond epitopes originating from nonsynonymous single-nucleotide variants in the coding exome. Here we review emerging approaches for identifying, prioritizing, and immunologically targeting personalized neoantigens using polyvalent cancer vaccines and T-cell receptor gene therapy. SIGNIFICANCE: Several major challenges currently impede the clinical efficacy of neoantigen-directed immunotherapy, such as the relative infrequency of immunogenic neoantigens, suboptimal potency and priming of de novo tumor-specific T cells, and tumor cell-intrinsic and -extrinsic mechanisms of immune evasion. A deeper understanding of these biological barriers could help facilitate the development of effective and durable immunotherapy for any type of cancer, including immunologically "cold" tumors that are otherwise therapeutically resistant.

Two Human Monoclonal HLA-Reactive Antibodies Cross-React with Mamu-B*008, a Rhesus Macaque MHC Allotype Associated with Control of Simian Immunodeficiency Virus Replication

MHC class I molecules play an important role in adaptive immune responses against intracellular pathogens. These molecules are highly polymorphic, and many allotypes have been characterized. In a transplantation setting, a mismatch between MHC allotypes may initiate an alloimmune response. Rhesus macaques (Macaca mulatta, Mamu) are valuable as a preclinical model species in transplantation research as well as to evaluate the safety and efficacy of vaccine candidates. In both lines of research, the availability of nonhuman primate MHC-reactive mAbs may enable in vitro monitoring and detection of presence of particular Mamu molecules. In this study, we screened a collection of thoroughly characterized HLA class I-specific human mAbs for cross-reactivity with rhesus macaque MHC class I allotypes. Two mAbs, OK4F9 and OK4F10, recognize an epitope that is defined by isoleucine (I) at amino acid position 142 that is present on the Indian rhesus macaque Mamu-B*008:01 allotype, which is an allotype known to be associated with elite control of SIV replication. The reactive pattern of a third mAb, MUS4H4, is more complex and includes an epitope shared on Mamu-A2*05:01 and -B*001:01-encoded Ags. This is the first description, to our knowledge, of human HLA-reactive mAbs that can recognize Mamu allotypes, and these can be useful tools for in vitro monitoring the presence of the relevant allelic products. Moreover, OK4F9 and OK4F10 can be powerful mAbs for application in SIV-related research.

Haplotype Motif-Based Models for KIR-Genotype Informed Selection of Hematopoietic Cell Donors Fail to Predict Outcome of Patients With Myelodysplastic Syndromes or Secondary Acute Myeloid Leukemia

Results from registry studies suggest that harnessing Natural Killer (NK) cell reactivity mediated through Killer cell Immunoglobulin-like Receptors (KIR) could reduce the risk of relapse after allogeneic Hematopoietic Cell Transplantation (HCT). Several competing models have been developed to classify donors as KIR-advantageous or disadvantageous. Basically, these models differ by grouping donors based on distinct KIR–KIR–ligand combinations or by haplotype motif assignment. This study aimed to validate different models for unrelated donor selection for patients with Myelodysplatic Syndromes (MDS) or secondary Acute Myeloid Leukemia (sAML). In a joint retrospective study of the European Society for Blood and Marrow Transplantation (EBMT) and the Center for International Blood and Marrow Transplant Research (CIBMTR) registry data from 1704 patients with secondary AML or MDS were analysed. The cohort consisted mainly of older patients (median age 61 years) with high risk disease who had received chemotherapy-based reduced intensity conditioning and anti-thymocyte globulin prior to allogeneic HCT from well-matched unrelated stem cell donors. The impact of the predictors on Overall Survival (OS) and relapse incidence was tested in Cox regression models adjusted for patient age, a modified disease risk index, performance status, donor age, HLA-match, sex-match, CMV-match, conditioning intensity, type of T-cell depletion and graft type. KIR genes were typed using high-resolution amplicon-based next generation sequencing. In univariable and multivariable analyses none of the models predicted OS and the risk of relapse consistently. Our results do not support the hypothesis that optimizing NK-mediated alloreactivity is possible by KIR-genotype informed selection of HLA-matched unrelated donors. However, in the context of allogeneic transplantation, NK-cell biology is complex and only partly understood. KIR-genes are highly diverse and current assignment of haplotype motifs based on the presence or absence of selected KIR genes is over-simplistic. As a consequence, further research is highly warranted and should integrate cutting edge knowledge on KIR genetics, and NK-cell biology into future studies focused on homogeneous groups of patients and treatment modalities.