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

Systematic discovery of neoepitope-HLA pairs for neoantigens shared among patients and tumor types

The broad application of precision cancer immunotherapies is limited by the number of validated neoepitopes that are common among patients or tumor types. To expand the known repertoire of shared neoantigen-human leukocyte antigen (HLA) complexes, we developed a high-throughput platform that coupled an in vitro peptide-HLA binding assay with engineered cellular models expressing individual HLA alleles in combination with a concatenated transgene harboring 47 common cancer neoantigens. From more than 24,000 possible neoepitope-HLA combinations, biochemical and computational assessment yielded 844 unique candidates, of which 86 were verified after immunoprecipitation mass spectrometry analyses of engineered, monoallelic cell lines. To evaluate the potential for immunogenicity, we identified T cell receptors that recognized select neoepitope-HLA pairs and elicited a response after introduction into human T cells. These cellular systems and our data on therapeutically relevant neoepitopes in their HLA contexts will aid researchers studying antigen processing as well as neoepitope targeting therapies.

Utilizing Massively Parallel Sequencing (MPS) of Human Leukocyte Antigen (HLA) Gene Polymorphism to Assess Relatedness in Deficiency Parentage Testing

In the realm of DNA testing with legal implications, the reliability and precision of genetic markers play a pivotal role in confirming or negating paternity claims. This study aimed to assess the potential utility of human leukocyte antigen (HLA) gene polymorphism through massively parallel sequencing (MPS) technology as robust forensic markers for parentage testing involving genetic deficiencies. It sought to redefine the significance of HLA genes in this context. Data on autosomal short tandem repeat (aSTR) mutational events across 18 paternity cases involving 16 commonly employed microsatellite loci were presented. In instances where traditional aSTR analysis failed to establish statistical certainty, kinship determination was pursued via HLA genotyping, encompassing the amplification of 17 linked HLA loci. Within the framework of this investigation, phase-resolved genotypes for HLA genes were meticulously generated, resulting in the definition of 34 inherited HLA haplotypes. An impressive total of 274 unique HLA alleles, which were classified at either the field 3 or 4 level, were identified, including the discovery of four novel HLA alleles. Likelihood ratio (LR) values, which indicated the likelihood of the observed data under a true biological relationship versus no relationship, were subsequently calculated. The analysis of the LR values demonstrated that the HLA genes significantly enhanced kinship determination compared with the aSTR analysis. Combining LR values from aSTR markers and HLA loci yielded conclusive outcomes in duo paternity cases, showcasing the potential of HLA genes and MPS technology for deeper insights and diversity in genetic testing. Comprehensive reference databases and high-resolution HLA typing across diverse populations are essential. Reintegrating HLA alleles into forensic identification complements existing markers, creating a potent method for future forensic analysis.

86 novel HLA-E alleles discovered through full-gene sequencing of 6227 hematopoietic cell transplant patients and unrelated donors

Until recently the number of alleles of the nonclassical HLA class I gene HLA-E documented in the IPD-IMGT/HLA Database was small and as a result, the gene was often not considered to be notably polymorphic. Here, we describe our work in identifying and submitting 86 novel HLA-E alleles after full-gene single-molecule real-time (SMRT) DNA sequencing of 6227 DNA samples. These samples were comprised of 2468 patients undergoing hematopoietic cell transplantation and 3759 unrelated potential donors. A total of 111 unique HLA-E alleles were detected in this cohort. The majority of novel alleles (79.1%) contained polymorphisms in intronic regions, highlighting the significant undiscovered variation present in the noncoding regions of the HLA-E gene.

Development of an Immortalized Porcine Fibroblast Cell Panel With Different Swine Leukocyte Antigen Genotypes

Immortalized cell lines are valuable resources to expand the molecular characterization of major histocompatibility complex genes and their presented antigens. We generated a panel of immortalized cell lines by transfecting human telomerase reverse transcriptase (hTERT) into primary fibroblast cells prepared from ear, fetal, and lung tissues of 10 pigs from five breeds and successfully cultured them for 30-45 passages. The cell growth characteristic of the immortalized fibroblasts was similar to that of primary fibroblast, which was unable to form colonies on soft agar. The genotypes of major swine leukocyte antigen (SLA) genes, including three classical class I (SLA-1, -2, and -3) and three class II genes (DQB1, DRB1, and DQA), were determined using high-resolution typing. A total of 58 alleles, including a novel allele for SLA-2, were identified. Each cell line was unique. A cell line derived from a National Institutes of Health miniature pig was homozygous across the six major SLA genes. The expression levels of SLA classical class I genes varied among the cell lines and were slightly upregulated in the immortalized compared to the primary cells based on semiquantitative reverse transcription polymerase chain reaction. The immortalized porcine fibroblast cell lines with diverse SLA haplotypes that were developed in this study have potential to be applied in studies regarding the molecular characteristics and genetic structure of SLA genes and epitope-major histocompatibility complex interactions in pigs.

Personalized HLA typing leads to the discovery of novel HLA alleles and tumor‐specific HLA variants

Accurate and full‐length typing of the HLA region is important in many clinical and research settings. With the advent of next generation sequencing (NGS), several HLA typing algorithms have been developed, including many that are applicable to whole exome sequencing (WES). However, most of these solutions operate by providing the closest‐matched HLA allele among the known alleles in IPD‐IMGT/HLA Database. These database‐matching approaches have demonstrated very high performance when typing well characterized HLA alleles. However, as they rely on the completeness of the HLA database, they are not optimal for detecting novel or less well characterized alleles. Furthermore, the database‐matching approaches are also not adequate in the context of cancer, where a comprehensive characterization of somatic HLA variation and expression patterns of a tumor's HLA locus may guide therapy and clinical outcome, because of the pivotal role HLA alleles play in tumor antigen recognition and immune escape. Here, we describe a personalized HLA typing approach applied to WES data that leverages the strengths of database‐matching approaches while simultaneously allowing for the discovery of novel HLA alleles and tumor‐specific HLA variants, through the systematic integration of germline and somatic variant calling. We applied this approach on WES from 10 metastatic melanoma patients and validated the HLA typing results using HLA targeted NGS sequencing from patients where at least one HLA germline candidate was detected on Class I HLA. Targeted NGS sequencing confirmed 100% performance for the 1st and 2nd fields. In total, five out of the six detected HLA germline variants were because of Class I ambiguities at the third or fourth fields, and their detection recovered the correct HLA allele genotype. The sixth germline variant let to the formal discovery of a novel Class I allele. Finally, we demonstrated a substantially improved somatic variant detection accuracy in HLA alleles with a 91% of success rate in simulated experiments. The approach described here may allow the field to genotype more accurately using WES data, leading to the discovery of novel HLA alleles and help characterize the relationship between somatic variation in the HLA region and immunosurveillance.

HLA class II immunopeptidomics reveals that co-inherited HLA-allotypes within an extended haplotype can improve proteome coverage for immunosurveillance

Human leucocyte antigen (HLA) class II molecules in humans are encoded by three different loci, HLA-DR, -DQ, and -DP. These molecules share approximately 70% sequence similarity and all present peptide ligands to circulating T cells. While the peptide repertoires of numerous HLA-DR, -DQ, and -DP allotypes have been examined, there have been few reports on the combined repertoire of these co-inherited molecules expressed in a single cell as an extended HLA haplotype. Here we describe the endogenous peptide repertoire of a human B lymphoblastoid cell line (C1R) expressing the class II haplotype HLA-DR12/DQ7/DP4. We have identified 71350 unique naturally processed peptides presented collectively by HLA-DR12, HLA-DQ7, or HLA-DP4. The resulting "haplodome" is complemented by the cellular proteome defined by standard LC-MS/MS approaches. This large dataset has shed light on properties of these class II ligands especially the preference for membrane and extracellular source proteins. Our data also provides insights into the co-evolution of these conserved haplotypes of closely linked and co-inherited HLA molecules; which together increase sequence coverage of cellular proteins for immune surveillance with minimal overlap between each co-inherited HLA-class II allomorph.

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

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

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

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

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

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

Recombinant HA-based vaccine outperforms split and subunit vaccines in elicitation of influenza-specific CD4 T cells and CD4 T cell-dependent antibody responses in humans

Although traditional egg-based inactivated influenza vaccines can protect against infection, there have been significant efforts to develop improved formats to overcome disadvantages of this platform. Here, we have assessed human CD4 T cell responses to a traditional egg-based influenza vaccine with recently available cell-derived vaccines and recombinant baculovirus-derived vaccines. Adults were administered either egg-derived Fluzone^®, mammalian cell-derived Flucelvax^® or recombinant HA (Flublok^®). CD4 T cell responses to each HA protein were assessed by cytokine EliSpot and intracellular staining assays. The specificity and magnitude of antibody responses were quantified by ELISA and HAI assays. By all criteria, Flublok vaccine exhibited superior performance in eliciting both CD4 T cell responses and HA-specific antibody responses, whether measured by mean response magnitude or percent of responders. Although the mechanism(s) underlying this advantage is not yet clear, it is likely that both qualitative and quantitative features of the vaccines impact the response.

Rapid, highly accurate and cost-effective open-source simultaneous complete HLA typing and phasing of class I and II alleles using nanopore sequencing

Accurate rapid genotyping of the genes within the HLA region presents many difficulties because of the complexity of this region. Here we present the results of our proof of concept nanopore-based long read polymerase chain reaction (PCR) solution for HLA genotyping. For 15 HLA anthropology-based samples and 13 NHS Blood and Transplant derived samples 40 ng of genomic DNA underwent long-range PCR for class I and II HLA alleles. Pooled PCR products were sequenced on the Oxford Nanopore MinIoON R9.4.1 flow cell. Sequenced reads had HLA genotype assigned with HLA-LA. Called genotypes were compared with reference derived from a combination of short-read next-generation sequencing, Sanger sequence and/or single-site polymorphism (SSP) typing. For concordance, accuracy was 100%, 98.4%, 97.5% and 95.1% for the first, second, third and fourth fields, respectively, to four field accuracy where it was available, otherwise three field in 28 samples for class I calls and 17 samples for class II calls. Phasing of maternal and paternal alleles, as well as phasing based identification of runs of homozygosity, was shown successfully. Time for assay run was 8 hours and the reconstruction of HLA typing data was 15 minutes. Assay cost was £55 ($80USD)/sample. We have developed a rapid and cost-effective long-range PCR and nanopore sequencing-based assay that can genotype the genes within HLA region to up to four field accuracy, identify runs of homozygosity in HLA, reconstruct maternal and paternal haplotypes and can be scaled from multi-sample runs to a single sample.

A long road/read to rapid high-resolution HLA typing: The nanopore perspective

Next-generation sequencing (NGS) has been widely adopted for clinical HLA typing and advanced immunogenetics researches. Current methodologies still face challenges in resolving cis-trans ambiguity involving distant variant positions, and the turnaround time is affected by testing volume and batching. Nanopore sequencing may become a promising addition to the existing options for HLA typing. The technology delivered by the MinION sequencer of Oxford Nanopore Technologies (ONT) can record the ionic current changes during the translocation of DNA/RNA strands through transmembrane pores and translate the signals to sequence reads. It features simple and flexible library preparations, long sequencing reads, portable and affordable sequencing devices, and rapid, real-time sequencing. However, the error rate of the sequencing reads is high and remains a hurdle for its broad application. This review article will provide a brief overview of this technology and then focus on the opportunities and challenges of using nanopore sequencing for high-resolution HLA typing and immunogenetics research.

Performance of a multiplexed amplicon-based next-generation sequencing assay for HLA typing

BACKGROUND

Next-generation sequencing (NGS) has enabled efficient high-resolution typing of human leukocyte antigen (HLA) genes with minimal ambiguity. Most commercially available assays amplify individual or subgroup of HLA genes by long-range PCR followed by library preparation and sequencing. The AllType assay simplifies the workflow by amplifying 11 transplant-relevant HLA genes in one PCR reaction. Here, we report the performance of this unique workflow evaluated using 218 genetically diverse samples.

METHODS

Five whole genes (HLA-A/B/C/DQA1/DPA1) and six near-whole genes (HLA-DRB1/DRB345/DQB1/DPB1; excluding exon 1 and part of intron 1) were amplified in a multiplexed, long-range PCR. Manual library preparation was performed per manufacturer's protocol, followed by template preparation and chip loading on the Ion Chef, and sequencing on the Ion S5 sequencer. Pre-specified rules for quality control and repeat testing were followed; technologists were blinded to the reference results. The concordance between AllType and reference results was determined at 2-field resolution. We also describe the ranges of input DNA and library concentrations, read number per sample and per locus, and key health metrics in relation to typing results.

RESULTS

The concordance rates were 98.6%, 99.8% and 99.9% at the sample (n = 218), genotype (n = 1688), and allele (n = 3376) levels, respectively. Three genotypes were discordant, all of which shared the same G group typing results with the reference. Most ambiguous genotypes (116 out of 144, 80.6%) were due to the lack of exon 1 and intron 1 coverage for HLA-DRB1/DRB345/DQB1/DPB1 genes. A broad range of input DNA concentrations and library concentrations were tolerated. Per sample read numbers were adequate for accurate genotyping. Per locus read numbers showed some inter-lot variations, and a trend toward improved inter-locus balance was observed with later lots of reagents.

CONCLUSION

The AllType assay on the Ion Chef/Ion S5 platform offers a robust and efficient workflow for clinical HLA typing at the 2-field resolution. The multiplex PCR strategy simplifies the laboratory procedure without compromising the typing accuracy.