1
|
Dashti M, Malik MZ, Nizam R, Jacob S, Al-Mulla F, Thanaraj TA. Evaluation of HLA typing content of next-generation sequencing datasets from family trios and individuals of arab ethnicity. Front Genet 2024; 15:1407285. [PMID: 38859936 PMCID: PMC11163123 DOI: 10.3389/fgene.2024.1407285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Accepted: 05/07/2024] [Indexed: 06/12/2024] Open
Abstract
Introduction: HLA typing is a critical tool in both clinical and research applications at the individual and population levels. Benchmarking studies have indicated HLA-HD as the preferred tool for accurate and comprehensive HLA allele calling. The advent of next-generation sequencing (NGS) has revolutionized genetic analysis by providing high-throughput sequencing data. This study aims to evaluate, using the HLA-HD tool, the HLA typing content of whole exome, whole genome, and HLA-targeted panel sequence data from the consanguineous population of Arab ethnicity, which has been underrepresented in prior benchmarking studies. Methods: We utilized sequence data from family trios and individuals, sequenced on one or more of the whole exome, whole genome, and HLA-targeted panel sequencing technologies. The performance and resolution across various HLA genes were evaluated. We incorporated a comparative quality control analysis, assessing the results obtained from HLA-HD by comparing them with those from the HLA-Twin tool to authenticate the accuracy of the findings. Results: Our analysis found that alleles across 29 HLA loci can be successfully and consistently typed from NGS datasets. Clinical-grade whole exome sequencing datasets achieved the highest consistency rate at three-field resolution, followed by targeted HLA panel, research-grade whole exome, and whole genome datasets. Discussion: The study catalogues HLA typing consistency across NGS datasets for a large array of HLA genes and highlights assessments regarding the feasibility of utilizing available NGS datasets in HLA allele studies. These findings underscore the reliability of HLA-HD for HLA typing in underrepresented populations and demonstrate the utility of various NGS technologies in achieving accurate HLA allele calling.
Collapse
Affiliation(s)
| | | | | | | | - Fahd Al-Mulla
- Department of Genetics and Bioinformatics, Dasman Diabetes Institute, Kuwait City, Kuwait
| | | |
Collapse
|
2
|
Sverchkova A, Burkholz S, Rubsamen R, Stratford R, Clancy T. Integrative HLA typing of tumor and adjacent normal tissue can reveal insights into the tumor immune response. BMC Med Genomics 2024; 17:37. [PMID: 38281021 PMCID: PMC10821267 DOI: 10.1186/s12920-024-01808-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Accepted: 01/12/2024] [Indexed: 01/29/2024] Open
Abstract
BACKGROUND The HLA complex is the most polymorphic region of the human genome, and its improved characterization can help us understand the genetics of human disease as well as the interplay between cancer and the immune system. The main function of HLA genes is to recognize "non-self" antigens and to present them on the cell surface to T cells, which instigate an immune response toward infected or transformed cells. While sequence variation in the antigen-binding groove of HLA may modulate the repertoire of immunogenic antigens presented to T cells, alterations in HLA expression can significantly influence the immune response to pathogens and cancer. METHODS RNA sequencing was used here to accurately genotype the HLA region and quantify and compare the level of allele-specific HLA expression in tumors and patient-matched adjacent normal tissue. The computational approach utilized in the study types classical and non-classical Class I and Class II HLA alleles from RNA-seq while simultaneously quantifying allele-specific or personalized HLA expression. The strategy also uses RNA-seq data to infer immune cell infiltration into tumors and the corresponding immune cell composition of matched normal tissue, to reveal potential insights related to T cell and NK cell interactions with tumor HLA alleles. RESULTS The genotyping method outperforms existing RNA-seq-based HLA typing tools for Class II HLA genotyping. Further, we demonstrate its potential for studying tumor-immune interactions by applying the method to tumor samples from two different subtypes of breast cancer and their matched normal breast tissue controls. CONCLUSIONS The integrative RNA-seq-based HLA typing approach described in the study, coupled with HLA expression analysis, neoantigen prediction and immune cell infiltration, may help increase our understanding of the interplay between a patient's tumor and immune system; and provide further insights into the immune mechanisms that determine a positive or negative outcome following treatment with immunotherapy such as checkpoint blockade.
Collapse
Affiliation(s)
- Angelina Sverchkova
- NEC OncoImmunity, Oslo Cancer Cluster, Innovation Park, Oslo, Norway
- Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Scott Burkholz
- Flow Pharma, Inc, Warrensville Heights, Galaxy Parkway, OH, 4829, USA
| | - Reid Rubsamen
- Flow Pharma, Inc, Warrensville Heights, Galaxy Parkway, OH, 4829, USA
- University Hospitals, Cleveland Medical Center, Cleveland, OH, USA
- Case Western Reserve School of Medicine, Cleveland, OH, USA
| | - Richard Stratford
- NEC OncoImmunity, Oslo Cancer Cluster, Innovation Park, Oslo, Norway
| | - Trevor Clancy
- NEC OncoImmunity, Oslo Cancer Cluster, Innovation Park, Oslo, Norway.
| |
Collapse
|
3
|
Ricker CA, Meli K, Van Allen EM. Historical perspective and future directions: computational science in immuno-oncology. J Immunother Cancer 2024; 12:e008306. [PMID: 38191244 PMCID: PMC10826578 DOI: 10.1136/jitc-2023-008306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/07/2023] [Indexed: 01/10/2024] Open
Abstract
Immuno-oncology holds promise for transforming patient care having achieved durable clinical response rates across a variety of advanced and metastatic cancers. Despite these achievements, only a minority of patients respond to immunotherapy, underscoring the importance of elucidating molecular mechanisms responsible for response and resistance to inform the development and selection of treatments. Breakthroughs in molecular sequencing technologies have led to the generation of an immense amount of genomic and transcriptomic sequencing data that can be mined to uncover complex tumor-immune interactions using computational tools. In this review, we discuss existing and emerging computational methods that contextualize the composition and functional state of the tumor microenvironment, infer the reactivity and clonal dynamics from reconstructed immune cell receptor repertoires, and predict the antigenic landscape for immune cell recognition. We further describe the advantage of multi-omics analyses for capturing multidimensional relationships and artificial intelligence techniques for integrating omics data with histopathological and radiological images to encapsulate patterns of treatment response and tumor-immune biology. Finally, we discuss key challenges impeding their widespread use and clinical application and conclude with future perspectives. We are hopeful that this review will both serve as a guide for prospective researchers seeking to use existing tools for scientific discoveries and inspire the optimization or development of novel tools to enhance precision, ultimately expediting advancements in immunotherapy that improve patient survival and quality of life.
Collapse
Affiliation(s)
- Cora A Ricker
- Medical Oncology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
| | - Kevin Meli
- Medical Oncology, Dana-Farber Cancer Institute, Boston, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| | | |
Collapse
|
4
|
Claeys A, Merseburger P, Staut J, Marchal K, Van den Eynden J. Benchmark of tools for in silico prediction of MHC class I and class II genotypes from NGS data. BMC Genomics 2023; 24:247. [PMID: 37161318 PMCID: PMC10170851 DOI: 10.1186/s12864-023-09351-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2023] [Accepted: 04/30/2023] [Indexed: 05/11/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Arne Claeys
- Department of Human Structure and Repair, Ghent University, Ghent, Belgium
- Cancer Research Institute Ghent, Ghent, Belgium
| | - Peter Merseburger
- Department of Human Structure and Repair, Ghent University, Ghent, Belgium
- Cancer Research Institute Ghent, Ghent, Belgium
| | - Jasper Staut
- Department of Human Structure and Repair, Ghent University, Ghent, Belgium
| | - Kathleen Marchal
- Cancer Research Institute Ghent, Ghent, Belgium
- Department of Information Technology, Ghent University, IDLab, Ghent, Belgium
- Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
| | - Jimmy Van den Eynden
- Department of Human Structure and Repair, Ghent University, Ghent, Belgium.
- Cancer Research Institute Ghent, Ghent, Belgium.
| |
Collapse
|
5
|
Solomon BD, Zheng H, Dillon LW, Goldman J, Hourigan CS, Heath J, Khatri P. Prediction of HLA genotypes from single-cell transcriptome data. Front Immunol 2023; 14:1146826. [PMID: 37180102 PMCID: PMC10167300 DOI: 10.3389/fimmu.2023.1146826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 04/04/2023] [Indexed: 05/15/2023] Open
Abstract
The human leukocyte antigen (HLA) locus plays a central role in adaptive immune function and has significant clinical implications for tissue transplant compatibility and allelic disease associations. Studies using bulk-cell RNA sequencing have demonstrated that HLA transcription may be regulated in an allele-specific manner and single-cell RNA sequencing (scRNA-seq) has the potential to better characterize these expression patterns. However, quantification of allele-specific expression (ASE) for HLA loci requires sample-specific reference genotyping due to extensive polymorphism. While genotype prediction from bulk RNA sequencing is well described, the feasibility of predicting HLA genotypes directly from single-cell data is unknown. Here we evaluate and expand upon several computational HLA genotyping tools by comparing predictions from human single-cell data to gold-standard, molecular genotyping. The highest 2-field accuracy averaged across all loci was 76% by arcasHLA and increased to 86% using a composite model of multiple genotyping tools. We also developed a highly accurate model (AUC 0.93) for predicting HLA-DRB345 copy number in order to improve genotyping accuracy of the HLA-DRB locus. Genotyping accuracy improved with read depth and was reproducible at repeat sampling. Using a metanalytic approach, we also show that HLA genotypes from PHLAT and OptiType can generate ASE ratios that are highly correlated (R2 = 0.8 and 0.94, respectively) with those derived from gold-standard genotyping.
Collapse
Affiliation(s)
| | - Hong Zheng
- Institute for Immunity, Transplantation and Infection, School of Medicine, Stanford University, Stanford, CA, United States
- Center for Biomedical Informatics Research, Department of Medicine, School of Medicine, Stanford University, Stanford, CA, United States
| | - Laura W. Dillon
- Laboratory of Myeloid Malignancies, National Heart Lung and Blood Institute, Bethesda, MD, United States
| | - Jason D. Goldman
- Swedish Center for Research and Innovation, Swedish Medical Center, Seattle, WA, United States
- Providence St. Joseph Health, Renton, WA, United States
- Division of Allergy & Infectious Diseases, University of Washington, Seattle, WA, United States
| | - Christopher S. Hourigan
- Laboratory of Myeloid Malignancies, National Heart Lung and Blood Institute, Bethesda, MD, United States
| | - James R. Heath
- Institute for Systems Biology, Seattle, WA, United States
- Department of Bioengineering, University of Washington, Seattle, WA, United States
| | - Purvesh Khatri
- Institute for Immunity, Transplantation and Infection, School of Medicine, Stanford University, Stanford, CA, United States
- Center for Biomedical Informatics Research, Department of Medicine, School of Medicine, Stanford University, Stanford, CA, United States
| |
Collapse
|
6
|
Lybaert L, Lefever S, Fant B, Smits E, De Geest B, Breckpot K, Dirix L, Feldman SA, van Criekinge W, Thielemans K, van der Burg SH, Ott PA, Bogaert C. Challenges in neoantigen-directed therapeutics. Cancer Cell 2023; 41:15-40. [PMID: 36368320 DOI: 10.1016/j.ccell.2022.10.013] [Citation(s) in RCA: 31] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/07/2022] [Revised: 08/19/2022] [Accepted: 10/11/2022] [Indexed: 11/11/2022]
Abstract
A fundamental prerequisite for the efficacy of cancer immunotherapy is the presence of functional, antigen-specific T cells within the tumor. Neoantigen-directed therapy is a promising strategy that aims at targeting the host's immune response against tumor-specific antigens, thereby eradicating cancer cells. Initial forays have been made in clinical environments utilizing vaccines and adoptive cell therapy; however, many challenges lie ahead. We provide an in-depth overview of the current state of the field with an emphasis on in silico neoantigen discovery and the clinical aspects that need to be addressed to unlock the full potential of this therapy.
Collapse
Affiliation(s)
| | | | | | - Evelien Smits
- Center for Oncological Research, University of Antwerp, 2610 Wilrijk, Belgium
| | - Bruno De Geest
- Department of Pharmaceutics, Ghent University, 9000 Ghent, Belgium
| | - Karine Breckpot
- Laboratory of Molecular and Cellular Therapy, Department of Biomedical Sciences, Vrije Universiteit Brussel, Brussels, Belgium
| | - Luc Dirix
- Translational Cancer Research Unit, Center for Oncological Research, Faculty of Medicine and Health Sciences, University of Antwerp, Antwerp, Belgium
| | - Steven A Feldman
- Center for Cancer Cell Therapy, Stanford University School of Medicine, Stanford, CA, USA
| | - Wim van Criekinge
- Department of Data Analysis and Mathematical Modelling, Ghent University, Ghent, Belgium
| | - Kris Thielemans
- Laboratory of Molecular and Cellular Therapy, Department of Biomedical Sciences, Vrije Universiteit Brussel, Brussels, Belgium
| | - Sjoerd H van der Burg
- Medical Oncology, Oncode Institute, Leiden University Medical Center, Leiden, the Netherlands
| | - Patrick A Ott
- Department of Medical Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA
| | | |
Collapse
|
7
|
Thuesen NH, Klausen MS, Gopalakrishnan S, Trolle T, Renaud G. Benchmarking freely available HLA typing algorithms across varying genes, coverages and typing resolutions. Front Immunol 2022; 13:987655. [PMID: 36426357 PMCID: PMC9679531 DOI: 10.3389/fimmu.2022.987655] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Accepted: 10/10/2022] [Indexed: 11/02/2023] Open
Abstract
Identifying the specific human leukocyte antigen (HLA) allele combination of an individual is crucial in organ donation, risk assessment of autoimmune and infectious diseases and cancer immunotherapy. However, due to the high genetic polymorphism in this region, HLA typing requires specialized methods. We investigated the performance of five next-generation sequencing (NGS) based HLA typing tools with a non-restricted license namely HLA*LA, Optitype, HISAT-genotype, Kourami and STC-Seq. This evaluation was done for the five HLA loci, HLA-A, -B, -C, -DRB1 and -DQB1 using whole-exome sequencing (WES) samples from 829 individuals. The robustness of the tools to lower depth of coverage (DOC) was evaluated by subsampling and HLA typing 230 WES samples at DOC ranging from 1X to 100X. The HLA typing accuracy was measured across four typing resolutions. Among these, we present two clinically-relevant typing resolutions (P group and pseudo-sequence), which specifically focus on the peptide binding region. On average, across the five HLA loci examined, HLA*LA was found to have the highest typing accuracy. For the individual loci, HLA-A, -B and -C, Optitype's typing accuracy was the highest and HLA*LA had the highest typing accuracy for HLA-DRB1 and -DQB1. The tools' robustness to lower DOC data varied widely and further depended on the specific HLA locus. For all Class I loci, Optitype had a typing accuracy above 95% (according to the modification of the amino acids in the functionally relevant portion of the HLA molecule) at 50X, but increasing the DOC beyond even 100X could still improve the typing accuracy of HISAT-genotype, Kourami, and STC-seq across all five HLA loci as well as HLA*LA's typing accuracy for HLA-DQB1. HLA typing is also used in studies of ancient DNA (aDNA), which is often based on sequencing data with lower quality and DOC. Interestingly, we found that Optitype's typing accuracy is not notably impaired by short read length or by DNA damage, which is typical of aDNA, as long as the DOC is sufficiently high.
Collapse
Affiliation(s)
- Nikolas Hallberg Thuesen
- Evaxion Biotech, Copenhagen, Denmark
- Department of Health Technology, Section for Bioinformatics, Technical University of Denmark, Lyngby, Denmark
| | | | - Shyam Gopalakrishnan
- Section for Hologenomics, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | | | - Gabriel Renaud
- Department of Health Technology, Section for Bioinformatics, Technical University of Denmark, Lyngby, Denmark
| |
Collapse
|