Very long haplotype tracts characterized at high resolution from HLA homozygous cell lines

The 18th International HLA & Immunogenetics workshop project report: Creating fully representative MHC reference haplotypes

A fundamental endeavor of the International Histocompatibility and Immunogenetics Workshop (IHIW) was assembling a collection of DNA samples homozygous through the MHC genomic region. This collection proved invaluable for assay development in the histocompatibility and immunogenetics field, for generating the human reference genome, and furthered our understanding of MHC diversity. Defined by their HLA-A, -B, -C and -DRB1 alleles, the combined frequency of the haplotypes from these individuals is ~20% in Europe. Thus, a significant proportion of MHC haplotypes, both common and rare throughout the world, and including many associated with disease, are not yet represented. In this workshop component, we are collecting the next generation of MHC -homozygous samples, to expand, diversify and modernize this critical community resource that has been foundational to the field. We asked laboratories worldwide to identify samples homozygous through all HLA class I and/or HLA class II genes, or through whole-genome SNP genotyping or sequencing, to have extensive homozygosity tracts within the MHC region. The focus is non-Europeans or those having HLA haplotypes less common in Europeans. Through this effort, we have obtained samples from 537 individuals representing 294 distinct haplotypes, as determined by their HLA class I and II alleles, and an additional 50 haplotypes distinct in HLA class I or II alleles. Although we have expanded the diversity, many populations remain underrepresented, particularly from Africa, and we encourage further participation. The data will serve as a resource for investigators seeking to characterize variation across the MHC genomic region for disease and population studies.

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.

Next-generation HLA typing of 382 International Histocompatibility Working Group reference B-lymphoblastoid cell lines: Report from the 17th International HLA and Immunogenetics Workshop

Extended molecular characterization of HLA genes in the IHWG reference B-lymphoblastoid cell lines (B-LCLs) was one of the major goals for the 17th International HLA and Immunogenetics Workshop (IHIW). Although reference B-LCLs have been examined extensively in previous workshops complete high-resolution typing was not completed for all the classical class I and class II HLA genes. To address this, we conducted a single-blind study where select panels of B-LCL genomic DNA samples were distributed to multiple laboratories for HLA genotyping by next-generation sequencing methods. Identical cell panels comprised of 24 and 346 samples were distributed and typed by at least four laboratories in order to derive accurate consensus HLA genotypes. Overall concordance rates calculated at both 2- and 4-field allele-level resolutions ranged from 90.4% to 100%. Concordance for the class I genes ranged from 91.7 to 100%, whereas concordance for class II genes was variable; the lowest observed at HLA-DRB3 (84.2%). At the maximum allele-resolution 78 B-LCLs were defined as homozygous for all 11 loci. We identified 11 novel exon polymorphisms in the entire cell panel. A comparison of the B-LCLs NGS HLA genotypes with the HLA genotypes catalogued in the IPD-IMGT/HLA Database Cell Repository, revealed an overall allele match at 68.4%. Typing discrepancies between the two datasets were mostly due to the lower-resolution historical typing methods resulting in incomplete HLA genotypes for some samples listed in the IPD-IMGT/HLA Database Cell Repository. Our approach of multiple-laboratory NGS HLA typing of the B-LCLs has provided accurate genotyping data. The data generated by the tremendous collaborative efforts of the 17th IHIW participants is useful for updating the current cell and sequence databases and will be a valuable resource for future studies.

Introduction: MHC/KIR and governance of specificity

The MHC controls specificity, to ensure that appropriate immune responses are mounted to invading pathogens whilst maintaining tolerance to the host. It encodes molecules that act as sentinels, providing a snapshot of the health of the interior and exterior of the cell for immune surveillance. To maintain the ability to respond appropriately to any disease requires a delicate balance of expression and function, and many subtleties of the system have been described at the gene, individual and population level. The main players are the highly polymorphic classical MHC class I and class II molecules, as well as some non-classical loci of both types. Transporter associated with antigen processing (TAP) peptide transporters, proteasome components and Tapasin, encoded within the MHC, are also involved in selection of peptide for presentation. The plethora of mechanisms microorganisms use to subvert immune recognition, through blocking these antigen processing and presentation pathways, attests to the importance of HLA in resistance to infection. There is continued interest in MHC genetics in its own right, as well as in relation to KIR, to transplantation, infection, autoimmunity and reproduction. Also of topical interest, cancer immunotherapy through checkpoint inhibition depends on highly specific recognition of cancer peptide antigen and continued expression of HLA molecules. Here, we briefly introduce some background to the MHC/KIR axis in man. This special issue of immunogenetics expands on these topics, in humans and other model species.

Role of major histocompatibility complex variation in graft-versus-host disease after hematopoietic cell transplantation

Graft-versus-host disease (GVHD) remains a significant potentially life-threatening complication of allogeneic hematopoietic cell transplantation (HCT). Since the discovery of the human leukocyte antigen (HLA) system over 50 years ago, significant advances have clarified the nature of HLA variation between transplant recipients and donors as a chief etiology of GVHD. New information on coding and non-coding gene variation and GVHD risk provides clinicians with options to consider selected mismatched donors when matched donors are not available. These advances have increased the availability of unrelated donors for patients in need of a transplant and have lowered the overall morbidity and mortality of HCT.

Sequences of 95 human MHC haplotypes reveal extreme coding variation in genes other than highly polymorphic HLA class I and II

The most polymorphic part of the human genome, the MHC, encodes over 160 proteins of diverse function. Half of them, including the HLA class I and II genes, are directly involved in immune responses. Consequently, the MHC region strongly associates with numerous diseases and clinical therapies. Notoriously, the MHC region has been intractable to high-throughput analysis at complete sequence resolution, and current reference haplotypes are inadequate for large-scale studies. To address these challenges, we developed a method that specifically captures and sequences the 4.8-Mbp MHC region from genomic DNA. For 95 MHC homozygous cell lines we assembled, de novo, a set of high-fidelity contigs and a sequence scaffold, representing a mean 98% of the target region. Included are six alternative MHC reference sequences of the human genome that we completed and refined. Characterization of the sequence and structural diversity of the MHC region shows the approach accurately determines the sequences of the highly polymorphic HLA class I and HLA class II genes and the complex structural diversity of complement factor C4A/C4B. It has also uncovered extensive and unexpected diversity in other MHC genes; an example is MUC22, which encodes a lung mucin and exhibits more coding sequence alleles than any HLA class I or II gene studied here. More than 60% of the coding sequence alleles analyzed were previously uncharacterized. We have created a substantial database of robust reference MHC haplotype sequences that will enable future population scale studies of this complicated and clinically important region of the human genome.

HLA class I variation in Iranian Lur and Kurd populations: high haplotype and allotype diversity with an abundance of KIR ligands

HLA-A, -B and -C alleles of 285 individuals, representing three Iranian Lur populations and one Iranian Kurd population were sequenced completely, yielding human leukocyte antigen (HLA) class I genotypes at high resolution and filling four fields of the official HLA nomenclature. Each population has 87-99 alleles, evenly distributed between the three HLA class I genes, 145 alleles being identified in total. These alleles were already known, named and deposited in the HLA database. The alleles form 316 different HLA A-B-C haplotypes, with each population having between 80 and 112 haplotypes. The four Iranian populations form a related group that is distinguished from other populations, including other Iranians. All four KIR ligands - the A3/11, Bw4, C1 and C2 epitopes - are well represented, particularly Bw4, which is carried by three high-frequency allotypes: HLA-A*24:02, HLA-A*32:01 and HLA-B*51:01. In the Lur and Kurd populations, between 82% and 94% of individuals have the Bw4 epitope, the ligand for KIR3DL1. HLA-B*51:01 is likely of Neandertal origin and associated with Behcet's disease, also known as the Silk Road disease. The Lordegan Lur have the highest frequency of HLA-B*51:01 in the world. This allele is present on 46 Lur and Kurd haplotypes. Present at lower frequency is HLA-B*51:08, which is also associated with Behcet's disease. In the four Iranian populations, 31 haplotypes encode both Bw4(+) HLA-A and Bw4(+) HLA-B, a dual combination of Bw4 epitopes that is relatively rare in other populations, worldwide. This study both demonstrates and emphasizes the value of studying HLA class I polymorphism at highest resolution in anthropologically well-defined populations.

Defining KIR and HLA Class I Genotypes at Highest Resolution via High-Throughput Sequencing

The physiological functions of natural killer (NK) cells in human immunity and reproduction depend upon diverse interactions between killer cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands: HLA-A, HLA-B, and HLA-C. The genomic regions containing the KIR and HLA class I genes are unlinked, structurally complex, and highly polymorphic. They are also strongly associated with a wide spectrum of diseases, including infections, autoimmune disorders, cancers, and pregnancy disorders, as well as the efficacy of transplantation and other immunotherapies. To facilitate study of these extraordinary genes, we developed a method that captures, sequences, and analyzes the 13 KIR genes and HLA-A, HLA-B, and HLA-C from genomic DNA. We also devised a bioinformatics pipeline that attributes sequencing reads to specific KIR genes, determines copy number by read depth, and calls high-resolution genotypes for each KIR gene. We validated this method by using DNA from well-characterized cell lines, comparing it to established methods of HLA and KIR genotyping, and determining KIR genotypes from 1000 Genomes sequence data. This identified 116 previously uncharacterized KIR alleles, which were all demonstrated to be authentic by sequencing from source DNA via standard methods. Analysis of just two KIR genes showed that 22% of the 1000 Genomes individuals have a previously uncharacterized allele or a structural variant. The method we describe is suited to the large-scale analyses that are needed for characterizing human populations and defining the precise HLA and KIR factors associated with disease. The methods are applicable to other highly polymorphic genes.