HLA-VBSeq: accurate HLA typing at full resolution from whole-genome sequencing data

Artificial intelligence and neoantigens: paving the path for precision cancer immunotherapy

Cancer immunotherapy has witnessed rapid advancement in recent years, with a particular focus on neoantigens as promising targets for personalized treatments. The convergence of immunogenomics, bioinformatics, and artificial intelligence (AI) has propelled the development of innovative neoantigen discovery tools and pipelines. These tools have revolutionized our ability to identify tumor-specific antigens, providing the foundation for precision cancer immunotherapy. AI-driven algorithms can process extensive amounts of data, identify patterns, and make predictions that were once challenging to achieve. However, the integration of AI comes with its own set of challenges, leaving space for further research. With particular focus on the computational approaches, in this article we have explored the current landscape of neoantigen prediction, the fundamental concepts behind, the challenges and their potential solutions providing a comprehensive overview of this rapidly evolving field.

Individual and population-level variability in HLA-DR associated immunogenicity risk of biologics used for the treatment of rheumatoid arthritis

Hypothesis

While conventional in silico immunogenicity risk assessments focus on measuring immunogenicity based on the potential of therapeutic proteins to be processed and presented by a global population-wide set of human leukocyte antigen (HLA) alleles to T cells, future refinements might adjust for HLA allele frequencies in different geographic regions or populations, as well for as individuals in those populations. Adjustment by HLA allele distribution may reveal risk patterns that are specific to population groups or individuals, which current methods that rely on global-population HLA prevalence may obscure.

Key findings

This analysis uses HLA frequency-weighted binding predictions to define immunogenicity risk for global and sub-global populations. A comparison of assessments tuned for North American/European versus Japanese/Asian populations suggests that the potential for anti-therapeutic responses (anti-therapeutic antibodies or ATA) for several commonly prescribed Rheumatoid Arthritis (RA) therapeutic biologics may differ, significantly, between the Caucasian and Japanese populations. This appears to align with reports of differing product-related immunogenicity that is observed in different populations.

Relevance to clinical practice

Further definition of population-level (regional) and individual patient-specific immunogenic risk profiles may enable prescription of the RA therapeutic with the highest probability of success to each patient, depending on their population of origin and/or their individual HLA background. Furthermore, HLA-specific immunogenicity outcomes data are limited, thus there is a need to expand HLA-association studies that examine the relationship between HLA haplotype and ATA in the clinic.

Large-scale integrative analysis of juvenile idiopathic arthritis for new insight into its pathogenesis

BACKGROUND

Juvenile idiopathic arthritis (JIA) is one of the most prevalent rheumatic disorders in children and is classified as an autoimmune disease (AID). While a robust genetic contribution to JIA etiology has been established, the exact pathogenesis remains unclear.

METHODS

To prioritize biologically interpretable susceptibility genes and proteins for JIA, we conducted transcriptome-wide and proteome-wide association studies (TWAS/PWAS). Then, to understand the genetic architecture of JIA, we systematically analyzed single-nucleotide polymorphism (SNP)-based heritability, a signature of natural selection, and polygenicity. Next, we conducted HLA typing using multi-ethnicity RNA sequencing data. Additionally, we examined the T cell receptor (TCR) repertoire at a single-cell level to explore the potential links between immunity and JIA risk.

RESULTS

We have identified 19 TWAS genes and two PWAS proteins associated with JIA risks. Furthermore, we observe that the heritability and cell type enrichment analysis of JIA are enriched in T lymphocytes and HLA regions and that JIA shows higher polygenicity compared to other AIDs. In multi-ancestry HLA typing, B*45:01 is more prevalent in African JIA patients than in European JIA patients, whereas DQA1*01:01, DQA1*03:01, and DRB1*04:01 exhibit a higher frequency in European JIA patients. Using single-cell immune repertoire analysis, we identify clonally expanded T cell subpopulations in JIA patients, including CXCL13+BHLHE40+ TH cells which are significantly associated with JIA risks.

CONCLUSION

Our findings shed new light on the pathogenesis of JIA and provide a strong foundation for future mechanistic studies aimed at uncovering the molecular drivers of JIA.

Advancements and hurdles in the development of a vaccine for triple-negative breast cancer: A comprehensive review of multi-omics and immunomics strategies

Triple-Negative Breast Cancer (TNBC) presents a significant challenge in oncology due to its aggressive behavior and limited therapeutic options. This review explores the potential of immunotherapy, particularly vaccine-based approaches, in addressing TNBC. It delves into the role of immunoinformatics in creating effective vaccines against TNBC. The review first underscores the distinct attributes of TNBC and the importance of tumor antigens in vaccine development. It then elaborates on antigen detection techniques such as exome sequencing, HLA typing, and RNA sequencing, which are instrumental in identifying TNBC-specific antigens and selecting vaccine candidates. The discussion then shifts to the in-silico vaccine development process, encompassing antigen selection, epitope prediction, and rational vaccine design. This process merges computational simulations with immunological insights. The role of Artificial Intelligence (AI) in expediting the prediction of antigens and epitopes is also emphasized. The review concludes by encapsulating how Immunoinformatics can augment the design of TNBC vaccines, integrating tumor antigens, advanced detection methods, in-silico strategies, and AI-driven insights to advance TNBC immunotherapy. This could potentially pave the way for more targeted and efficacious treatments.

The HLA-DRB1*09:01-DQB1*03:03 haplotype is associated with the risk for late-onset Alzheimer's disease in APOE [Formula: see text]4-negative Japanese adults

Late-onset Alzheimer's disease (LOAD) is the most common cause of dementia among those older than 65 years. The onset of LOAD is influenced by neuroinflammation. The human leukocyte antigen (HLA) system is involved in regulating inflammatory responses. Numerous HLA alleles and their haplotypes have shown varying associations with LOAD in diverse populations, yet their impact on the Japanese population remains to be elucidated. Here, we conducted a comprehensive investigation into the associations between LOAD and HLA alleles within the Japanese population. Using whole-genome sequencing (WGS) data from 303 LOAD patients and 1717 cognitively normal (CN) controls, we identified four-digit HLA class I alleles (A, B, and C) and class II alleles (DRB1, DQB1, and DPB1). We found a significant association between the HLA-DRB1*09:01-DQB1*03:03 haplotype and LOAD risk in APOE [Formula: see text]4-negative samples (odds ratio = 1.81, 95% confidence interval = 1.38-2.38, P = 2.03[Formula: see text]). These alleles not only showed distinctive frequencies specific to East Asians but demonstrated a high degree of linkage disequilibrium in APOE [Formula: see text]4-negative samples (r2 = 0.88). Because HLA class II molecules interact with T-cell receptors (TCRs), we explored potential disparities in the diversities of TCR α chain (TRA) and β chain (TRB) repertoires between APOE [Formula: see text]4-negative LOAD and CN samples. Lower diversity of TRA repertoires was associated with LOAD in APOE [Formula: see text]4-negative samples, irrespective of the HLA DRB1*09:01-DQB1*03:03 haplotype. Our study enhances the understanding of the etiology of LOAD in the Japanese population and provides new insights into the underlying mechanisms of its pathogenesis.

SpecHLA enables full-resolution HLA typing from sequencing data

Reconstructing diploid sequences of human leukocyte antigen (HLA) genes, i.e., full-resolution HLA typing, from sequencing data is challenging. The high homogeneity across HLA genes and the high heterogeneity within HLA alleles complicate the identification of genomic source loci for sequencing reads. Here, we present SpecHLA, which utilizes fine-tuned reads binning and local assembly to achieve accurate full-resolution HLA typing. SpecHLA accepts sequencing data from paired-end, 10×-linked-reads, high-throughput chromosome conformation capture (Hi-C), Pacific Biosciences (PacBio), and Oxford Nanopore Technology (ONT). It can also incorporate pedigree data and genotype frequency to refine typing. In 32 Human Genome Structural Variation Consortium, Phase 2 (HGSVC2) samples, SpecHLA achieved 98.6% accuracy for G-group-resolution HLA typing, inferring entire HLA alleles with an average of three mismatches fewer, ten gaps fewer, and 590 bp less edit distance than HISAT-genotype per allele. Additionally, SpecHLA exhibited a 2-field typing accuracy of 98.6% in 875 real samples. Finally, SpecHLA detected HLA loss of heterozygosity with 99.7% specificity and 96.8% sensitivity in simulated samples of cancer cell lines.

Frequency changes in HLA-I alleles: A marker to guide immunotherapy in lung adenocarcinoma patients and its relationship with tumor mutational burden and PD-L1 expression

BACKGROUND

The aim of the study was to investigate differences in HLA-I alleles between lung adenocarcinoma patients and healthy controls and determine their association with PD-L1 expression and tumor mutational burden (TMB) to understand the mechanism underlying lung adenocarcinoma susceptibility.

METHODS

Differences in HLA allele frequencies between the two groups were analyzed in a case-control study. PD-L1 expression and TMB in lung adenocarcinoma patients were determined and their relationships with HLA-I were analyzed.

RESULTS

The lung adenocarcinoma group showed significantly higher HLA-A*30:01 (p = 0.0067, odds ratio [OR], 1.834; 95% confidence interval [CI]: 1.176-2.860), B*13:02 (p = 0.0050, OR, 1.855; 95% CI: 1.217-2.829), and C*06:02 (p = 0.0260, OR, 1.478; 95% CI: 1.060-2.060) and significantly lower B*51:01 (p = 0.0290, OR, 0.6019; 95% CI: 0.3827-0.9467), and C*14:02 (p = 0.0255, OR, 0.5089; 95% CI: 0.2781-0.9312) than the control group. Haplotype analysis results showed that HLA-A*30:01-B*13:02 (p = 0.0100, OR, 1.909; 95% CI: 1.182-3.085), A*11:01-C*01:02 (p = 0.0056, OR, 1.909; 95% CI: 1.182-3.085), A*30:01-C*06:02 (p = 0.0111, OR, 1.846; 95% CI: 1.147-2.969), and B*13:02-C*06:02 (p = 0.0067, OR, 1.846; 95% CI: 1.147-2.969) frequencies significantly increased and B*51:01-C*14:02 (p = 0.0219, OR, 0.490; 95% CI: 0.263-0.914) frequency significantly decreased in lung adenocarcinoma patients. Three-locus haplotype analysis showed that HLA-A*30:01-B*13:02-C*06:02 frequency (p = 0.0100, OR, 1.909; 95% CI: 1.182-3.085) significantly increased in patients.

CONCLUSION

HLA-A*30:01, B*13:02, and C*06:02 may be the susceptibility genes and HLA-B*51:01 and C*14:01 act as the resistance genes of lung adenocarcinoma. The changes in HLA-I allele frequencies had no association with PD-L1 expression and TMB among these patients.

Benchmark of tools for in silico prediction of MHC class I and class II genotypes from NGS data

BACKGROUND

The Human Leukocyte Antigen (HLA) genes are a group of highly polymorphic genes that are located in the Major Histocompatibility Complex (MHC) region on chromosome 6. The HLA genotype affects the presentability of tumour antigens to the immune system. While knowledge of these genotypes is of utmost importance to study differences in immune responses between cancer patients, gold standard, PCR-derived genotypes are rarely available in large Next Generation Sequencing (NGS) datasets. Therefore, a variety of methods for in silico NGS-based HLA genotyping have been developed, bypassing the need to determine these genotypes with separate experiments. However, there is currently no consensus on the best performing tool.

RESULTS

We evaluated 13 MHC class I and/or class II HLA callers that are currently available for free academic use and run on either Whole Exome Sequencing (WES) or RNA sequencing data. Computational resource requirements were highly variable between these tools. Three orthogonal approaches were used to evaluate the accuracy on several large publicly available datasets: a direct benchmark using PCR-derived gold standard HLA calls, a correlation analysis with population-based allele frequencies and an analysis of the concordance between the different tools. The highest MHC-I calling accuracies were found for Optitype (98.0%) and arcasHLA (99.4%) on WES and RNA sequencing data respectively, while for MHC-II HLA-HD was the most accurate tool for both data types (96.2% and 99.4% on WES and RNA data respectively).

CONCLUSION

The optimal strategy for HLA genotyping from NGS data depends on the availability of either WES or RNA data, the size of the dataset and the available computational resources. If sufficient resources are available, we recommend Optitype and HLA-HD for MHC-I and MHC-II genotype calling respectively.

Prediction Model with HLA-A*33:03 Reveals Number of Days to Develop Liver Cancer from Blood Test

The development of liver cancer in patients with hepatitis B is a major problem, and several models have been reported to predict the development of liver cancer. However, no predictive model involving human genetic factors has been reported to date. For the items incorporated in the prediction model reported so far, we selected items that were significant in predicting liver carcinogenesis in Japanese patients with hepatitis B and constructed a prediction model of liver carcinogenesis by the Cox proportional hazard model with the addition of Human Leukocyte Antigen (HLA) genotypes. The model, which included four items-sex, age at the time of examination, alpha-fetoprotein level (log₁₀AFP) and presence or absence of HLA-A*33:03-revealed an area under the receiver operating characteristic curve (AUROC) of 0.862 for HCC prediction within 1 year and an AUROC of 0.863 within 3 years. A 1000 repeated validation test resulted in a C-index of 0.75 or higher, or sensitivity of 0.70 or higher, indicating that this predictive model can distinguish those at high risk of developing liver cancer within a few years with high accuracy. The prediction model constructed in this study, which can distinguish between chronic hepatitis B patients who develop hepatocellular carcinoma (HCC) early and those who develop HCC late or not, is clinically meaningful.

Advances of mRNA vaccine in tumor: a maze of opportunities and challenges

High-frequency mutations in tumor genomes could be exploited as an asset for developing tumor vaccines. In recent years, with the tremendous breakthrough in genomics, intelligence algorithm, and in-depth insight of tumor immunology, it has become possible to rapidly target genomic alterations in tumor cell and rationally select vaccine targets. Among a variety of candidate vaccine platforms, the early application of mRNA was limited by instability low efficiency and excessive immunogenicity until the successful development of mRNA vaccines against SARS-COV-2 broken of technical bottleneck in vaccine preparation, allowing tumor mRNA vaccines to be prepared rapidly in an economical way with good performance of stability and efficiency. In this review, we systematically summarized the classification and characteristics of tumor antigens, the general process and methods for screening neoantigens, the strategies of vaccine preparations and advances in clinical trials, as well as presented the main challenges in the current mRNA tumor vaccine development.

The genomic and immune landscape of long-term survivors of high-grade serous ovarian cancer

Fewer than half of all patients with advanced-stage high-grade serous ovarian cancers (HGSCs) survive more than five years after diagnosis, but those who have an exceptionally long survival could provide insights into tumor biology and therapeutic approaches. We analyzed 60 patients with advanced-stage HGSC who survived more than 10 years after diagnosis using whole-genome sequencing, transcriptome and methylome profiling of their primary tumor samples, comparing this data to 66 short- or moderate-term survivors. Tumors of long-term survivors were more likely to have multiple alterations in genes associated with DNA repair and more frequent somatic variants resulting in an increased predicted neoantigen load. Patients clustered into survival groups based on genomic and immune cell signatures, including three subsets of patients with BRCA1 alterations with distinctly different outcomes. Specific combinations of germline and somatic gene alterations, tumor cell phenotypes and differential immune responses appear to contribute to long-term survival in HGSC.

Hepatobiliary Tumor Organoids Reveal HLA Class I Neoantigen Landscape and Antitumoral Activity of Neoantigen Peptide Enhanced with Immune Checkpoint Inhibitors

Neoantigen-directed therapy lacks preclinical models recapitulating neoantigen characteristics of original tumors. It is urgent to develop a platform to assess T cell response for neoantigen screening. Here, immunogenic potential of neoantigen-peptides of tumor tissues and matched organoids (n = 27 pairs) are analyzed by Score tools with whole genome sequencing (WGS)-based human leukocyte antigen (HLA)-class-I algorithms. The comparisons between 9203 predicted neoantigen-peptides from 2449 mutations of tumor tissues and 9991 ones from 2637 mutations of matched organoids demonstrate that organoids preserved majority of genetic features, HLA alleles, and similar neoantigen landscape of original tumors. Higher neoantigen load is observed in tumors with early stage. Multiomics analysis combining WGS, RNA-seq, single-cell RNA-seq, mass spectrometry filters out 93 candidate neoantigen-peptides with strong immunogenic potential for functional validation in five organoids. Immunogenic peptides are defined by inducing increased CD107aCD137IFN-γ expressions and IFN-γ secretion of CD8 cells in flow cytometry and enzyme-linked immunosorbent assay assays. Nine immunogenic peptides shared by at least two individuals are validated, including peptide from TP53^R90S . Organoid killing assay confirms the antitumor activity of validated immunogenic peptide-reactive CD8 cells, which is further enhanced in the presence of immune checkpoint inhibitors. The study characterizes HLA-class-I neoantigen landscape in hepatobiliary tumor, providing practical strategy with tumor organoid model for neoantigen-peptide identification in personalized immunotherapy.

Environmental risk factors and genetic markers of Kaposi's sarcoma associated herpesvirus infection among Uygur population in Xinjiang, China

OBJECTIVES

Kaposi sarcoma-associated herpesvirus (KSHV) is endemic in Xinjiang, China. Determinants of KSHV seropositivity among high-risk groups are not well understood. We seek to identify genetic and environmental predisposing factors for KSHV infection among Uygurs in this endemic region.

METHODS

A cross-sectional study was performed among Uygur population in Xinjiang, China. KSHV-antibodies were detected using immunofluorescence assay (IFA) and human leukocyte antigen (HLA) alleles were genotyped. Univariate and multivariate logistic regression analysis were applied to explore the environmental and genetic risk factors of KSHV seropositivity.

RESULTS

A total of 721 participants were included. The seroprevalence of KSHV was 24.1% among this population. Sweet-food preference (OR 1.85, 95%CI 1.03-3.34), and coronary heart disease (OR 1.91, 95 %CI 1.24-2.94) were statistically correlated with KSHV infection. HLA-DQB1*06:09 were found to significantly increase the risk of KSHV infection under all 3 models (OR_Allelic =4.06; OR_Dominant =3.27; OR_Recessive =8.06). Six SNPs (SNP0260, SNP0361, SNP0797, SNP0852, SNP1159, SNP1375) in the DQB1 and DRB1 region and haploid type GTCTAACTAATC in block 17 were statistically associated with KSHV infection.

CONCLUSION

We demonstrated that genetic variations in HLA-DQB1/DRB1 and environmental risk factors were strongly associated with KSHV infection among this population. This article is protected by copyright. All rights reserved.

Genomic predictors of response to PD-1 inhibition in children with germline DNA replication repair deficiency

Cancers arising from germline DNA mismatch repair deficiency or polymerase proofreading deficiency (MMRD and PPD) in children harbour the highest mutational and microsatellite insertion–deletion (MS-indel) burden in humans. MMRD and PPD cancers are commonly lethal due to the inherent resistance to chemo-irradiation. Although immune checkpoint inhibitors (ICIs) have failed to benefit children in previous studies, we hypothesized that hypermutation caused by MMRD and PPD will improve outcomes following ICI treatment in these patients. Using an international consortium registry study, we report on the ICI treatment of 45 progressive or recurrent tumors from 38 patients. Durable objective responses were observed in most patients, culminating in a 3 year survival of 41.4%. High mutation burden predicted response for ultra-hypermutant cancers (>100 mutations per Mb) enriched for combined MMRD + PPD, while MS-indels predicted response in MMRD tumors with lower mutation burden (10–100 mutations per Mb). Furthermore, both mechanisms were associated with increased immune infiltration even in ‘immunologically cold’ tumors such as gliomas, contributing to the favorable response. Pseudo-progression (flare) was common and was associated with immune activation in the tumor microenvironment and systemically. Furthermore, patients with flare who continued ICI treatment achieved durable responses. This study demonstrates improved survival for patients with tumors not previously known to respond to ICI treatment, including central nervous system and synchronous cancers, and identifies the dual roles of mutation burden and MS-indels in predicting sustained response to immunotherapy.

Hypermutation and microsatellite burden determine responses and long-term survival following PD-1 blockade in children and young adults with refractory cancers resulting from germline DNA replication repair deficiency.

In silico tools and databases for designing cancer immunotherapy

Immunotherapy is a rapidly growing therapy for cancer which have numerous benefits over conventional treatments like surgery, chemotherapy, and radiation. Overall survival of cancer patients has improved significantly due to the use of immunotherapy. It acts as a novel pillar for treating different malignancies from their primary to the metastatic stage. Recent preferments in high-throughput sequencing and computational immunology leads to the development of targeted immunotherapy for precision oncology. In the last few decades, several computational methods and resources have been developed for designing immunotherapy against cancer. In this review, we have summarized cancer-associated genomic, transcriptomic, and mutation profile repositories. We have also enlisted in silico methods for the prediction of vaccine candidates, HLA binders, cytokines inducing peptides, and potential neoepitopes. Of note, we have incorporated the most important bioinformatics pipelines and resources for the designing of cancer immunotherapy. Moreover, to facilitate the scientific community, we have developed a web portal entitled ImmCancer (https://webs.iiitd.edu.in/raghava/immcancer/), comprises cancer immunotherapy tools and repositories.

A New Human Leukocyte Antigen Typing Algorithm Combined With Currently Available Genotyping Tools Based on Next-Generation Sequencing Data and Guidelines to Select the Most Likely Human Leukocyte Antigen Genotype

Background

High-precision human leukocyte antigen (HLA) genotyping is crucial for anti-cancer immunotherapy, but existing tools predicting HLA genotypes using next-generation sequencing (NGS) data are insufficiently accurate.

Materials and Methods

We compared availability, accuracy, correction score, and complementary ratio of eight HLA genotyping tools (OptiType, HLA-HD, PHLAT, seq2HLA, arcasHLA, HLAscan, HLA*LA, and Kourami) using 1,005 cases from the 1000 Genomes Project data. We created a new HLA-genotyping algorithm combining tools based on the precision and the accuracy of tools’ combinations. Then, we assessed the new algorithm’s performance in 39 in-house samples with normal whole-exome sequencing (WES) data and polymerase chain reaction–sequencing-based typing (PCR-SBT) results.

Results

Regardless of the type of tool, the calls presented by more than six tools concordantly showed high accuracy and precision. The accuracy of the group with at least six concordant calls was 100% (97/97) in HLA-A, 98.2% (112/114) in HLA-B, 97.3% (142/146) in HLA-C. The precision of the group with at least six concordant calls was over 98% in HLA-ABC. We additionally calculated the accuracy of the combination tools considering the complementary ratio of each tool and the accuracy of each tool, and the accuracy was over 98% in all groups with six or more concordant calls. We created a new algorithm that matches the above results. It was to select the HLA type if more than six out of eight tools presented a matched type. Otherwise, determine the HLA type experimentally through PCR-SBT. When we applied the new algorithm to 39 in-house cases, there were more than six matching calls in all HLA-A, B, and C, and the accuracy of these concordant calls was 100%.

Conclusions

HLA genotyping accuracy using NGS data could be increased by combining the current HLA genotyping tools. This new algorithm could also be useful for preliminary screening to decide whether to perform an additional PCR-based experimental method instead of using tools with NGS data.

Whole-Genome Sequencing of a 900-Year-Old Human Skeleton Supports Two Past Migration Events from the Russian Far East to Northern Japan

Recent studies on paleogenomics have reported some Paleolithic and Neolithic genomes that have provided new insights into the human population history in East and Northeast Asia. However, there remain some cases where more recent migration events need to be examined to elucidate the detailed formation process of local populations. Although the area around northern Japan is one of the regions archaeologically suggested to have been affected by migration waves after the Neolithic period, the genetic source of these migrations are still unclear. Thus, genomic data from such past migrant populations would be highly informative to clarify the detailed formation process of local populations in this region. Here, we report the genome sequence of a 900-year-old adult female (NAT002) belonging to the prehistoric Okhotsk people, who have been considered to be the past migrants to northern Japan after the Neolithic period. We found a close relationship between NAT002 and modern Lower Amur populations and past admixture events between the Amur, Jomon, and Kamchatka ancestries. The admixture dating suggested migration of Amur-related ancestry at approximately 1,600 BP, which is compatible with the archaeological evidence regarding the settlement of the Okhotsk people. Our results also imply migration of Kamchatka-related ancestry at approximately 2,000 BP. In addition, human leukocyte antigen (HLA) typing detected the HLA-B*40 allele, which is reported to increase the risk of arthritis, suggesting the genetic vulnerability of NAT002 to hyperostosis, which was observed around her chest clavicle.

An immunogenetic view of COVID-19

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

Technological advances in cancer immunity: from immunogenomics to single-cell analysis and artificial intelligence

Immunotherapies play critical roles in cancer treatment. However, given that only a few patients respond to immune checkpoint blockades and other immunotherapeutic strategies, more novel technologies are needed to decipher the complicated interplay between tumor cells and the components of the tumor immune microenvironment (TIME). Tumor immunomics refers to the integrated study of the TIME using immunogenomics, immunoproteomics, immune-bioinformatics, and other multi-omics data reflecting the immune states of tumors, which has relied on the rapid development of next-generation sequencing. High-throughput genomic and transcriptomic data may be utilized for calculating the abundance of immune cells and predicting tumor antigens, referring to immunogenomics. However, as bulk sequencing represents the average characteristics of a heterogeneous cell population, it fails to distinguish distinct cell subtypes. Single-cell-based technologies enable better dissection of the TIME through precise immune cell subpopulation and spatial architecture investigations. In addition, radiomics and digital pathology-based deep learning models largely contribute to research on cancer immunity. These artificial intelligence technologies have performed well in predicting response to immunotherapy, with profound significance in cancer therapy. In this review, we briefly summarize conventional and state-of-the-art technologies in the field of immunogenomics, single-cell and artificial intelligence, and present prospects for future research.

Association of Protective HLA-A With HLA-B∗27 Positive Ankylosing Spondylitis

Objectives

To further elucidate the role of the MHC in ankylosing spondylitis by typing 17 genes, searching for HLA-B^∗27 independent associations and assessing the impact of sex on this male biased disease.

Methods

High-confidence two-field resolution genotyping was performed on 310 cases and 2196 controls using an n-1 concordance method. Protein-coding variants were called from next-generation sequencing reads using up to four software programs and the consensus result recorded. Logistic regression tests were applied to the dataset as a whole, and also in stratified sets based on sex or HLA-B^∗27 status. The amino acids driving association were also examined.

Results

Twenty-five HLA protein-coding variants were significantly associated to disease in the population. Three novel protective associations were found in a HLA-B^∗27 positive population, HLA-A^∗24:02 (OR = 0.4, CI = 0.2–0.7), and HLA-A amino acids Leu95 and Gln156. We identified a key set of seven loci that were common to both sexes, and robust to change in sample size. Stratifying by sex uncovered three novel risk variants restricted to the female population (HLA-DQA1^∗04.01, -DQB1^∗04:02, -DRB1^∗08:01; OR = 2.4–3.1). We also uncovered a set of neutral variants in the female population, which in turn conferred strong effects in the male set, highlighting how population composition can lead to the masking of true associations.

Conclusion

Population stratification allowed for a nuanced investigation into the tightly linked MHC region, revealing novel HLA-B^∗27 signals as well as replicating previous HLA-B^∗27 dependent results. This dissection of signals may help to elucidate sex biased disease predisposition and clinical progression.

Multi-step screening of neoantigens' HLA- and TCR-interfaces improves prediction of survival

Improvement of risk stratification through prognostic biomarkers may enhance the personalization of cancer patient monitoring and treatment. We used Ancer, an immunoinformatic CD8, CD4, and regulatory T cell neoepitope screening system, to perform an advanced neoantigen analysis of genomic data derived from the urothelial cancer cohort of The Cancer Genome Atlas. Ancer demonstrated improved prognostic stratification and five-year survival prediction compared to standard analyses using tumor mutational burden or neoepitope identification using NetMHCpan and NetMHCIIpan. The superiority of Ancer, shown in both univariate and multivariate survival analyses, is attributed to the removal of neoepitopes that do not contribute to tumor immunogenicity based on their homology with self-epitopes. This analysis suggests that the presence of a higher number of unique, non-self CD8- and CD4-neoepitopes contributes to cancer survival, and that prospectively defining these neoepitopes using Ancer is a novel prognostic or predictive biomarker.

Systematic comparative study of computational methods for HLA typing from next-generation sequencing

The human leukocyte antigen (HLA) system plays an important role in hematopoietic stem cell transplantation (HSCT) and organ transplantations, immune disorders as well as oncological immunotherapy. However, HLA typing remains a challenging task due to the high level of polymorphism and homology among HLA genes. Based on the high-throughput next-generation sequencing data, new HLA typing algorithms and software tools were developed. But there is still a deficit of systematic comparative studies to assist in the selection of the optimal analytical approaches under different conditions. Here, we present a detailed comparison of eight software tools for HLA typing on different real datasets (whole-genome sequencing, whole-exome sequencing and transcriptomic sequencing data) and in-silico samples with different sequencing lengths, depths, and error rates. We figure out the algorithms with the best efficiency in different scenarios, and demonstrate the effect of different raw reads on analytical performances. Our results provide a comprehensive picture of specifications and performances of the eight existing HLA genotyping algorithms, which could assist researchers in selecting the most appropriate tool for specific raw datasets.

Benchmarking the Human Leukocyte Antigen Typing Performance of Three Assays and Seven Next-Generation Sequencing-Based Algorithms

With the great progress made recently in next generation sequencing (NGS) technology, sequencing accuracy and throughput have increased, while the cost for data has decreased. Various human leukocyte antigen (HLA) typing algorithms and assays have been developed and have begun to be used in clinical practice. In this study, we compared the HLA typing performance of three HLA assays and seven NGS-based HLA algorithms and assessed the impact of sequencing depth and length on HLA typing accuracy based on 24 benchmarked samples. The algorithms HISAT-genotype and HLA-HD showed the highest accuracy at both the first field and the second field resolution, followed by HLAscan. Our internal capture-based HLA assay showed comparable performance with whole exome sequencing (WES). We found that the minimal depth was 100X for HISAT-genotype and HLA-HD to obtain more than 90% accuracy at the third field level. The top three algorithms were quite robust to the change of read length. Thus, we recommend using HISAT-genotype and HLA-HD for NGS-based HLA genotyping because of their higher accuracy and robustness to read length. We propose that a minimal sequence depth for obtaining more than 90% HLA typing accuracy at the third field level is 100X. Besides, targeting capture-based NGS HLA typing may be more suitable than WES in clinical practice due to its lower sequencing cost and higher HLA sequencing depth.

New challenges, new opportunities: Next generation sequencing and its place in the advancement of HLA typing

The Human Leukocyte Antigen (HLA) system has a critical role in immunorecognition, transplantation, and disease association. Early typing techniques provided the foundation for genotyping methods that revealed HLA as one of the most complex, polymorphic regions of the human genome. Next Generation Sequencing (NGS), the latest molecular technology introduced in clinical tissue typing laboratories, has demonstrated advantages over other established methods. NGS offers high-resolution sequencing of entire genes in time frames and price points considered unthinkable just a few years ago, contributing a wealth of data informing histocompatibility assessment and standards of clinical care. Although the NGS platforms share a high-throughput massively parallel processing model, differing chemistries provide specific strengths and weaknesses. Research-oriented Third Generation Sequencing and related advances in bioengineering continue to broaden the future of NGS in clinical settings. These diverse applications have demanded equally innovative strategies for data management and computational bioinformatics to support and analyze the unprecedented volume and complexity of data generated by NGS. We discuss some of the challenges and opportunities associated with NGS technologies, providing a comprehensive picture of the historical developments that paved the way for the NGS revolution, its current state and future possibilities for HLA typing.

State-of-the-art genome inference in the human MHC

The Major Histocompatibility Complex (MHC) on the short arm of chromosome 6 is associated with more diseases than any other region of the genome; it encodes the antigen-presenting Human Leukocyte Antigen (HLA) proteins and is one of the key immunogenetic regions of the genome. Accurate genome inference and interpretation of MHC association signals have traditionally been hampered by the region's uniquely complex features, such as high levels of polymorphism; inter-gene sequence homologies; structural variation; and long-range haplotype structures. Recent algorithmic and technological advances have, however, significantly increased the accessibility of genetic variation in the MHC; these developments include (i) accurate SNP-based HLA type imputation; (ii) genome graph approaches for variation-aware genome inference from next-generation sequencing data; (iii) long-read-based diploid de novo assembly of the MHC; (iv) cost-effective targeted MHC sequencing methods. Applied to hundreds of thousands of samples over the last years, these technologies have already enabled significant biological discoveries, for example in the field of autoimmune disease genetics. Remaining challenges concern the development of integrated methods that leverage haplotype-resolved de novo assembly of the MHC for the development of improved MHC genotyping methods for short reads and the construction of improved reference panels for SNP-based imputation. Improved genome inference in the MHC can crucially contribute to an improved genetic and functional understanding of many immune-related phenotypes and diseases.

Benchmarking HLA genotyping and clarifying HLA impact on survival in tumor immunotherapy

Human leukocyte antigen (HLA) genotyping gains intensive attention due to its critical role in cancer immunotherapy. It is still a challenging issue to generate reliable HLA genotyping results through in silico tools. In addition, the survival impact of HLA alleles in tumor prognosis and immunotherapy remains controversial. In this study, the benchmarking of HLA genotyping on TCGA is performed and a ‘Gun‐Bullet’ model which helps to clarify the survival impact of HLA allele is presented. The performance of HLA class I genotyping is generally better than class II. POLYSOLVER, OptiType, and xHLA perform generally better at HLA class I calling with an accuracy of 0.954, 0.949, and 0.937, respectively. HLA‐HD obtained the highest accuracy of 0.904 on HLA class II alleles calling. Each HLA genotyping tool displayed specific error patterns. The ensemble HLA calling from the top‐3 tools is superior to any individual one. HLA alleles show distinct survival impact among cancers. Cytolytic activity (CYT) was proposed as the underlying mechanism to interpret the survival impact of HLA alleles in the ‘Gun‐Bullet’ model for fighting cancer. A strong HLA allele plus a high tumor mutation burden (TMB) could stimulate intensive immune CYT, leading to extended survival. We established an up to now most reliable TCGA HLA benchmark dataset, composing of concordance alleles generated from eight prevalently used HLA genotyping tools. Our findings indicate that reliable HLA genotyping should be performed based on concordance alleles integrating multiple tools and incorporating TMB background with HLA genotype, which helps to improve the survival prediction compared to HLA genotyping alone.

Evolution of Genomic and T-cell Repertoire Heterogeneity of Malignant Pleural Mesothelioma Under Dasatinib Treatment

PURPOSE

Malignant pleural mesothelioma (MPM) is considered an orphan disease with few treatment options. Despite multimodality therapy, the majority of MPMs recur and eventually become refractory to any systemic treatment. One potential mechanism underlying therapeutic resistance may be intratumor heterogeneity (ITH), making MPM challenging to eradicate. However, the ITH architecture of MPM and its clinical impact have not been well studied.

EXPERIMENTAL DESIGN

We delineated the immunogenomic ITH by multiregion whole-exome sequencing and T-cell receptor (TCR) sequencing of 69 longitudinal MPM specimens from nine patients with resectable MPM, who were treated with dasatinib.

RESULTS

The median total mutation burden before dasatinib treatment was 0.65/Mb, similar with that of post-dasatinib treatment (0.62/Mb). The median proportion of mutations shared by any given pair of two tumor regions within the same tumors was 80% prior to and 83% post-dasatinib treatment indicating a relatively homogenous genomic landscape. T-cell clonality, a parameter indicating T-cell expansion and reactivity, was significantly increased in tumors after dasatinib treatment. Furthermore, on average, 82% of T-cell clones were restricted to individual tumor regions, with merely 6% of T-cell clones shared by all regions from the same tumors indicating profound TCR heterogeneity. Interestingly, patients with higher T-cell clonality and higher portion of T cells present across all tumor regions in post-dasatinib-treated tumors had significantly longer survival.

CONCLUSIONS

Despite the homogeneous genomic landscape, the TCR repertoire is extremely heterogeneous in MPM. Dasatinib may potentially induce T-cell response leading to improved survival.

High-Throughput Sequencing to Investigate Associations Between HLA Genes and Metamizole-Induced Agranulocytosis

Background and Objective: Agranulocytosis is a rare and potentially life-threatening complication of metamizole (dipyrone) intake that is characterized by a loss of circulating neutrophil granulocytes. While the mechanism underlying this adverse drug reaction is not well understood, involvement of the immune system has been suggested. In addition, associations between genetic variants in the Human Leukocyte Antigen (HLA) region and agranulocytosis induced by other drugs have been reported. The aim of the present study was to assess whether genetic variants in classical HLA genes are associated with the susceptibility to metamizole-induced agranulocytosis (MIA) in a European population by targeted resequencing of eight HLA genes.

Design: A case-control cohort of Swiss patients with a history of neutropenia or agranulocytosis associated with metamizole exposure (n = 53), metamizole-tolerant (n = 39) and unexposed controls (n = 161) was recruited for this study. A high-throughput resequencing (HTS) and high-resolution typing method was used to sequence and analyze eight HLA loci in a discovery subset of this cohort (n = 31 cases, n = 38 controls). Identified candidate alleles were investigated in the full Swiss cohort as well as in two independent cohorts from Germany and Spain using HLA imputation from genome-wide SNP array data. In addition, variant calling based on HTS data was performed in the discovery subset for the class I genes HLA-A, -B, and -C using the HLA-specific mapper hla-mapper.

Results: Eight candidate alleles (p < 0.05) were identified in the discovery subset, of which HLA-C^∗04:01 was associated with MIA in the full Swiss cohort (p < 0.01) restricted to agranulocytosis (ANC < 0.5 × 10⁹/L) cases. However, no candidate allele showed a consistent association in the Swiss, German and Spanish cohorts. Analysis of individual sequence variants in class I genes produced consistent results with HLA typing but did not reveal additional small nucleotide variants associated with MIA.

Conclusion: Our results do not support an HLA-restricted T cell-mediated immune mechanism for MIA. However, we established an efficient high-resolution (three-field) eight-locus HTS HLA resequencing method to interrogate the HLA region and demonstrated the feasibility of its application to pharmacogenetic studies.

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.

The premise of personalized immunotherapy for cancer dormancy

Progress in cancer therapies has resulted in improved survival of patients with early stage breast cancer. However, mortality remains high in patients with distant recurrence of the disease after initially successful treatment of early stage breast cancer. To this end, tumor recurrences have been attributed to the presence of dormant tumor cells in breast cancer patients and cancer survivors. Current clinical practice guidelines recommend a “wait and watch” approach for tumor recurrence. This is because of our limited understanding of tumor dormancy. Dormant tumor cells are quiescent, and thus, do not respond to chemotherapies or radiation therapies, and they are not operable. Therefore, immunotherapy is the only option for the treatment of tumor dormancy. However, gaps in our knowledge as to dormancy-specific antigens prevent a relapse preventing vaccine design. Here, we provide a critical review of cancer immunotherapy, and discuss empirical evidence related to naturally-occurring tumor dormancy and treatment-induced tumor dormancy at the site of primary tumor and in distant organs before and after cancer therapies. Finally, we suggest that personalized vaccines targeting dormancy-associated neoantigens, which can be given to patients with early stage disease after the completion of neoadjuvant therapies and tumor resection as well as to cancer survivors could eliminate relapse causing dormant cells and offer a cure for cancer.

Using RNA Sequencing to Characterize the Tumor Microenvironment

RNA sequencing (RNA-seq) is an integral tool in immunogenomics, allowing for interrogation of the transcriptome of a tumor and its microenvironment. Analytical methods to deconstruct the genomics data can then be applied to infer gene expression patterns associated with the presence of various immunocyte populations. High quality RNA-seq is possible from formalin-fixed, paraffin-embedded (FFPE), fresh-frozen, and fresh tissue, with a wide variety of sequencing library preparation methods, sequencing platforms, and downstream bioinformatics analyses currently available. Selection of an appropriate library preparation method is largely determined by tissue type, quality of RNA, and quantity of RNA. Downstream of sequencing, many analyses can be applied to the data, including differential gene expression analysis, immune gene signature analysis, gene pathway analysis, T/B-cell receptor inference, HLA inference, and viral transcript quantification. In this chapter, we will describe our workflow for RNA-seq from bulk tissue to evaluable data, including extraction of RNA, library preparation methods, sequencing of libraries, alignment and quality assurance of data, and initial downstream analyses of RNA-seq data to extract relevant immunogenomics features. Systems biology methods that draw additional insights by integrating these features are covered further in Chapters 28 - 30 .

SweHLA: the high confidence HLA typing bio-resource drawn from 1000 Swedish genomes

There is a need to accurately call human leukocyte antigen (HLA) genes from existing short-read sequencing data, however there is no single solution that matches the gold standard of Sanger sequenced lab typing. Here we aimed to combine results from available software programs, minimizing the biases of applied algorithm and HLA reference. The result is a robust HLA population resource for the published 1000 Swedish genomes, and a framework for future HLA interrogation. HLA 2nd-field alleles were called using four imputation and inference methods for the classical eight genes (class I: HLA-A, HLA-B, HLA-C; class II: HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, HLA-DRB1). A high confidence population set (SweHLA) was determined using an n−1 concordance rule for class I (four software) and class II (three software) alleles. Results were compared across populations and individual programs benchmarked to SweHLA. Per gene, 875 to 988 of the 1000 samples were genotyped in SweHLA; 920 samples had at least seven loci called. While a small fraction of reference alleles were common to all software (class I = 1.9% and class II = 4.1%), this did not affect the overall call rate. Gene-level concordance was high compared to European populations (>0.83%), with COX and PGF the dominant SweHLA haplotypes. We noted that 15/18 discordant alleles (delta allele frequency >2) were previously reported as disease-associated. These differences could in part explain across-study genetic replication failures, reinforcing the need to use multiple software solutions. SweHLA demonstrates a way to use existing NGS data to generate a population resource agnostic to individual HLA software biases.

Systematic genetic analysis of the MHC region reveals mechanistic underpinnings of HLA type associations with disease

The MHC region is highly associated with autoimmune and infectious diseases. Here we conduct an in-depth interrogation of associations between genetic variation, gene expression and disease. We create a comprehensive map of regulatory variation in the MHC region using WGS from 419 individuals to call eight-digit HLA types and RNA-seq data from matched iPSCs. Building on this regulatory map, we explored GWAS signals for 4083 traits, detecting colocalization for 180 disease loci with eQTLs. We show that eQTL analyses taking HLA type haplotypes into account have substantially greater power compared with only using single variants. We examined the association between the 8.1 ancestral haplotype and delayed colonization in Cystic Fibrosis, postulating that downregulation of RNF5 expression is the likely causal mechanism. Our study provides insights into the genetic architecture of the MHC region and pinpoints disease associations that are due to differential expression of HLA genes and non-HLA genes.

Neoantigen prediction from genomic and transcriptomic data

Tumor cells acquire distinct genetic characteristics as a means to survive and proliferate indefinitely. Changes in the genetic code can also translate in changes at the protein level, therefore creating a distinguishable signature unique for tumor cells, and absent in normal tissues. The presence of discernable moieties in tumors is particularly attractive because it represents a therapeutic opportunity to target tumor cells with specificity, while sparing non-transformed cells. In this sense neoantigens, short peptides containing a mutated sequence, are seen attractive therapeutic targets because of their confinement within tumor cells. Neoantigens can be recognized with high affinity and specificity by tumor-targeting T cells, which consequently can initiate a potent anti-tumor immune response. While this is feasible and it has been tested in numerous cancer types including melanoma, colon and lung cancer, to mention a few, there are technical challenges in identifying immunogenic neoantigens. In this manuscript we address the topic of neoantigen identification from tumor samples, offering a technical overview of the bioinformatic methods utilized to profile the neoantigenic load of tumor samples obtained from clinical specimens. This is meant to guide readers through the steps of neoantigen identification using genomic data, by suggesting tools and methods that can provide, with a high degree of confidence, reliable results for downstream in vitro and in vivo applications.

Predicting HLA class II antigen presentation through integrated deep learning

Accurate prediction of antigen presentation by human leukocyte antigen (HLA) class II molecules would be valuable for vaccine development and cancer immunotherapies. Current computational methods trained on in vitro binding data are limited by insufficient training data and algorithmic constraints. Here we describe MARIA (major histocompatibility complex analysis with recurrent integrated architecture; https://maria.stanford.edu/ ), a multimodal recurrent neural network for predicting the likelihood of antigen presentation from a gene of interest in the context of specific HLA class II alleles. In addition to in vitro binding measurements, MARIA is trained on peptide HLA ligand sequences identified by mass spectrometry, expression levels of antigen genes and protease cleavage signatures. Because it leverages these diverse training data and our improved machine learning framework, MARIA (area under the curve = 0.89-0.92) outperformed existing methods in validation datasets. Across independent cancer neoantigen studies, peptides with high MARIA scores are more likely to elicit strong CD4+ T cell responses. MARIA allows identification of immunogenic epitopes in diverse cancers and autoimmune disease.

Bioinformatics Strategies, Challenges, and Opportunities for Next Generation Sequencing-Based HLA Genotyping

The advent of next generation sequencing (NGS) has altered the face of genotyping the human leukocyte antigen (HLA) system in clinical, stem cell donor registry, and research contexts. NGS has led to a dramatically increased sequencing throughput at high accuracy, while being more time and cost efficient than precursor technologies. This has led to a broader and deeper profiling of the key genes in the human immunogenetic make-up. The rapid evolution of sequencing technologies is evidenced by the development of varied short-read sequencing platforms with differing read lengths and sequencing capacities to long-read sequencing platforms capable of profiling full genes without fragmentation. Concomitantly, there has been development of a diverse set of computational analyses and software tools developed to deal with the various strengths and limitations of the sequencing data generated by the different sequencing platforms. This review surveys the different modalities involved in generating NGS HLA profiling sequence data. It systematically describes various computational approaches that have been developed to achieve HLA genotyping to different degrees of resolution. At each stage, this review enumerates the drawbacks and advantages of each of the platforms and analysis approaches, thus providing a comprehensive picture of the current state of HLA genotyping technologies.

Best practices for bioinformatic characterization of neoantigens for clinical utility

Neoantigens are newly formed peptides created from somatic mutations that are capable of inducing tumor-specific T cell recognition. Recently, researchers and clinicians have leveraged next generation sequencing technologies to identify neoantigens and to create personalized immunotherapies for cancer treatment. To create a personalized cancer vaccine, neoantigens must be computationally predicted from matched tumor-normal sequencing data, and then ranked according to their predicted capability in stimulating a T cell response. This candidate neoantigen prediction process involves multiple steps, including somatic mutation identification, HLA typing, peptide processing, and peptide-MHC binding prediction. The general workflow has been utilized for many preclinical and clinical trials, but there is no current consensus approach and few established best practices. In this article, we review recent discoveries, summarize the available computational tools, and provide analysis considerations for each step, including neoantigen prediction, prioritization, delivery, and validation methods. In addition to reviewing the current state of neoantigen analysis, we provide practical guidance, specific recommendations, and extensive discussion of critical concepts and points of confusion in the practice of neoantigen characterization for clinical use. Finally, we outline necessary areas of development, including the need to improve HLA class II typing accuracy, to expand software support for diverse neoantigen sources, and to incorporate clinical response data to improve neoantigen prediction algorithms. The ultimate goal of neoantigen characterization workflows is to create personalized vaccines that improve patient outcomes in diverse cancer types.

HLA-VBSeq v2: improved HLA calling accuracy with full-length Japanese class-I panel

HLA-VBSeq is an HLA calling tool developed to infer the most likely HLA types from high-throughput sequencing data. However, there is still room for improvement in specific genetic groups because of the diversity of HLA alleles in human populations. Here, we present HLA-VBSeq v2, a software application that makes use of a new Japanese HLA reference panel to enhance calling accuracy for Japanese HLA class-I genes. Our analysis showed significant improvements in calling accuracy in all HLA regions, with prediction accuracies achieving over 99.0, 97.8, and 99.8% in HLA-A, B and C, respectively.

3.5KJPNv2: an allele frequency panel of 3552 Japanese individuals including the X chromosome

The first step towards realizing personalized healthcare is to catalog the genetic variations in a population. Since the dissemination of individual-level genomic information is strictly controlled, it will be useful to construct population-level allele frequency panels with easy-to-use interfaces. In the Tohoku Medical Megabank Project, we sequenced nearly 4000 individuals from a Japanese population and constructed an allele frequency panel of 3552 individuals after removing related samples. The panel is called the 3.5KJPNv2. It was constructed by using a standard pipeline including the 1KGP and gnomAD algorithms to reduce technical biases and to allow comparisons to other populations. Our database is the first large-scale panel providing the frequencies of variants present on the X chromosome and on the mitochondria in the Japanese population. All the data are available on our original database at https://jmorp.megabank.tohoku.ac.jp.

A new database provides information on the frequency of genetic variations within 3552 Japanese individuals, and facilitates comparisons with other populations. The reference panel, constructed by Kengo Kinoshita of Tohoku University, Sendai, and colleagues in Japan is also the first large-scale database to provide genetic variation frequency information on the X chromosome and mitochondrial DNA in the Japanese population. The methods used to sequence the genetic data are similar to those used in other large databases, allowing comparisons with other populations. The population size and methods used to compile the database overcome limitations in previous Japanese reference panels. This and similar databases that catalog genetic variations within populations can improve efforts towards personalizing healthcare and contribute to the study of human population genetics. The database is publicly available online.

ALPHLARD-NT: Bayesian Method for Human Leukocyte Antigen Genotyping and Mutation Calling through Simultaneous Analysis of Normal and Tumor Whole-Genome Sequence Data

Human leukocyte antigen (HLA) genes provide useful information on the relationship between cancer and the immune system. Despite the ease of obtaining these data through next-generation sequencing methods, interpretation of these relationships remains challenging owing to the complexity of HLA genes. To resolve this issue, we developed a Bayesian method, ALPHLARD-NT, to identify HLA germline and somatic mutations as well as HLA genotypes from whole-exome sequencing (WES) and whole-genome sequencing (WGS) data. ALPHLARD-NT showed 99.2% accuracy for WGS-based HLA genotyping and detected five HLA somatic mutations in 25 colon cancer cases. In addition, ALPHLARD-NT identified 88 HLA somatic mutations, including recurrent mutations and a novel HLA-B type, from WES data of 343 colon adenocarcinoma cases. These results demonstrate the potential of ALPHLARD-NT for conducting an accurate analysis of HLA genes even from low-coverage data sets. This method can become an essential tool for comprehensive analyses of HLA genes from WES and WGS data, helping to advance understanding of immune regulation in cancer as well as providing guidance for novel immunotherapy strategies.

Construction of full-length Japanese reference panel of class I HLA genes with single-molecule, real-time sequencing

Human leukocyte antigen (HLA) is a gene complex known for its exceptional diversity across populations, importance in organ and blood stem cell transplantation, and associations of specific alleles with various diseases. We constructed a Japanese reference panel of class I HLA genes (ToMMo HLA panel), comprising a distinct set of HLA-A, HLA-B, HLA-C, and HLA-H alleles, by single-molecule, real-time (SMRT) sequencing of 208 individuals included in the 1070 whole-genome Japanese reference panel (1KJPN). For high-quality allele reconstruction, we developed a novel pipeline, Primer-Separation Assembly and Refinement Pipeline (PSARP), in which the SMRT sequencing and additional short-read data were used. The panel consisted of 139 alleles, which were all extended from known IPD-IMGT/HLA sequences, contained 40 with novel variants, and captured more than 96.5% of allelic diversity in 1KJPN. These newly available sequences would be important resources for research and clinical applications including high-resolution HLA typing, genetic association studies, and analyzes of cis-regulatory elements.

Evaluation of computational programs to predict HLA genotypes from genomic sequencing data

Motivation

Despite being essential for numerous clinical and research applications, high-resolution human leukocyte antigen (HLA) typing remains challenging and laboratory tests are also time-consuming and labour intensive. With next-generation sequencing data becoming widely accessible, on-demand in silico HLA typing offers an economical and efficient alternative.

Results

In this study we evaluate the HLA typing accuracy and efficiency of five computational HLA typing methods by comparing their predictions against a curated set of > 1000 published polymerase chain reaction-derived HLA genotypes on three different data sets (whole genome sequencing, whole exome sequencing and transcriptomic sequencing data). The highest accuracy at clinically relevant resolution (four digits) we observe is 81% on RNAseq data by PHLAT and 99% accuracy by OptiType when limited to Class I genes only. We also observed variability between the tools for resource consumption, with runtime ranging from an average of 5 h (HLAminer) to 7 min (seq2HLA) and memory from 12.8 GB (HLA-VBSeq) to 0.46 GB (HLAminer) per sample.

While a minimal coverage is required, other factors also determine prediction accuracy and the results between tools do not correlate well. Therefore, by combining tools, there is the potential to develop a highly accurate ensemble method that is able to deliver fast, economical HLA typing from existing sequencing data.

Genome analyses for the Tohoku Medical Megabank Project towards establishment of personalized healthcare

Personalized healthcare (PHC) based on an individual's genetic make-up is one of the most advanced, yet feasible, forms of medical care. The Tohoku Medical Megabank (TMM) Project aims to combine population genomics, medical genetics and prospective cohort studies to develop a critical infrastructure for the establishment of PHC. To date, a TMM CommCohort (adult general population) and a TMM BirThree Cohort (birth+three-generation families) have conducted recruitments and baseline surveys. Genome analyses as part of the TMM Project will aid in the development of a high-fidelity whole-genome Japanese reference panel, in designing custom single-nucleotide polymorphism (SNP) arrays specific to Japanese, and in estimation of the biological significance of genetic variations through linked investigations of the cohorts. Whole-genome sequencing from >3,500 unrelated Japanese and establishment of a Japanese reference genome sequence from long-read data have been done. We next aim to obtain genotype data for all TMM cohort participants (>150,000) using our custom SNP arrays. These data will help identify disease-associated genomic signatures in the Japanese population, while genomic data from TMM BirThree Cohort participants will be used to improve the reference genome panel. Follow-up of the cohort participants will allow us to test the genetic markers and, consequently, contribute to the realization of PHC.

Empirical evaluation of variant calling accuracy using ultra-deep whole-genome sequencing data

In the design of whole-genome sequencing (WGS) studies, sequencing depth is a crucial parameter to define variant calling accuracy and study cost, with no standard recommendations having been established. We empirically evaluated the variant calling accuracy of the WGS pipeline using ultra-deep WGS data (approximately 410×). We randomly sampled sequence reads and constructed a series of simulation WGS datasets with a variety of gradual depths (n = 54; from 0.05× to 410×). Next, we evaluated the genotype concordances of the WGS data with those in the SNP microarray data or the WGS data using all the sequence reads. In addition, we assessed the accuracy of HLA allele genotyping using the WGS data with multiple software tools (PHLAT, HLA-VBseq, HLA-HD, and SNP2HLA). The WGS data with higher depths showed higher concordance rates, and >13.7× depth achieved as high as >99% of concordance. Comparisons with the WGS data using all the sequence reads showed that SNVs achieved >95% of concordance at 17.6× depth, whereas indels showed only 60% concordance. For the accuracy of HLA allele genotyping using the WGS data, 13.7× depth showed sufficient accuracy while performance heterogeneity among the software tools was observed (the highest concordance of 96.9% was observed with HLA-HD). Improvement in HLA genotyping accuracy by further increasing the depths was limited. These results suggest a medium degree of the WGS depth setting (approximately 15×) to achieve both accurate SNV calling and cost-effectiveness, whereas relatively higher depths are required for accurate indel calling.

The molecular landscape of synchronous colorectal cancer reveals genetic heterogeneity

Synchronous colorectal cancers (syCRCs), which present two or more lesions at diagnosis, are rare and pose a great challenge for clinical management. Although some predisposing factors associated with syCRCs have been studied with limited accession, the full repertoire of genomic events among the lesions within an individual and the causes of syCRCs remain unclear. We performed whole-exome sequencing of 40 surgical tumour samples of paired lesions from 20 patients to characterize the genetic alterations. Lesions from same patient showed distinct landscapes of somatic aberrations and shared few mutations, which suggests that they originate and develop independently, although they shared the similar genetic background. Canonical genes, such as APC, KRAS, TP53 and PIK3CA, were frequently mutated in the syCRCs, and most of them show different mutation profile compared with solitary colorectal cancer. We identified a recurrent somatic alteration (K15fs) in RPL22 in 25% of the syCRCs. Functional analysis indicated that mutated RPL22 may suppress cell apoptosis and promote the epithelial–mesenchymal transition (EMT). Potential drug targets were identified in several signalling pathways, and they present great discrepancy between lesions from the same patient. Our data show that the syCRCs within the same patient present great genetic heterogeneity, and they may be driven by distinct molecular events and develop independently. The discrepancy of potential drug targets and mutation burden in lesions from one patient provides valuable information in clinical management for patients with syCRCs.

Local mutational diversity drives intratumoral immune heterogeneity in non-small cell lung cancer

Combining whole exome sequencing, transcriptome profiling, and T cell repertoire analysis, we investigate the spatial features of surgically-removed biopsies from multiple loci in tumor masses of 15 patients with non-small cell lung cancer (NSCLC). This revealed that the immune microenvironment has high spatial heterogeneity such that intratumoral regional variation is as large as inter-personal variation. While the local total mutational burden (TMB) is associated with local T-cell clonal expansion, local anti-tumor cytotoxicity does not directly correlate with neoantigen abundance. Together, these findings caution against that immunological signatures can be predicted solely from TMB or microenvironmental analysis from a single locus biopsy.

The role of neoantigen in immune checkpoint blockade therapy

Immune checkpoint inhibitor induces tumor rejection by activated host immune system. The anti-tumor immune response consists of capture, presentation, recognition of neoantigen, as well as subsequent killing of tumor cell. Due to the interdependence among this series of stepwise events, neoantigen profoundly influences the efficacy of anti-immune checkpoint therapy. Moreover, the neoantigen-specific T cell reactivity is the cornerstone of multiple immunotherapies. In fact, several strategies targeting neoantigen have been attempted for synergetic effect with immune checkpoint inhibitor. Increasing neoantigen presentation to immune system by oncolytic virus, radiotherapy, or cancer vaccine is feasible to enhance neoantigen-specific T cell reactivity in theory. However, some obstacles have not been overcome in practice such as dynamic variation of neoantigen landscape, identification of potential neoantigen, maintenance of high T cell titer post vaccination. In addition, adoptive T cell transfer is another approach to enhance neoantigen-specific T cell reactivity, especially for patients with severe immunosuppression. In this review, we highlighted the advancements of neoantigen and innovative explorations of utilization of neoantigen repertoire in immune checkpoint blockade therapy.

ALPHLARD: a Bayesian method for analyzing HLA genes from whole genome sequence data

Background

Although human leukocyte antigen (HLA) genotyping based on amplicon, whole exome sequence (WES), and RNA sequence data has been achieved in recent years, accurate genotyping from whole genome sequence (WGS) data remains a challenge due to the low depth. Furthermore, there is no method to identify the sequences of unknown HLA types not registered in HLA databases.

Results

We developed a Bayesian model, called ALPHLARD, that collects reads potentially generated from HLA genes and accurately determines a pair of HLA types for each of HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, and -DRB1 genes at 3rd field resolution. Furthermore, ALPHLARD can detect rare germline variants not stored in HLA databases and call somatic mutations from paired normal and tumor sequence data. We illustrate the capability of ALPHLARD using 253 WES data and 25 WGS data from Illumina platforms. By comparing the results of HLA genotyping from SBT and amplicon sequencing methods, ALPHLARD achieved 98.8% for WES data and 98.5% for WGS data at 2nd field resolution. We also detected three somatic point mutations and one case of loss of heterozygosity in the HLA genes from the WGS data.

Conclusions

ALPHLARD showed good performance for HLA genotyping even from low-coverage data. It also has a potential to detect rare germline variants and somatic mutations in HLA genes. It would help to fill in the current gaps in HLA reference databases and unveil the immunological significance of somatic mutations identified in HLA genes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5169-9) contains supplementary material, which is available to authorized users.