ABSTRACT
Neoantigens are tumor-specific peptide sequences resulting from sources such as somatic DNA mutations. Upon loading onto major histocompatibility complex (MHC) molecules, they can trigger recognition by T cells. Accurate neoantigen identification is thus critical for both designing cancer vaccines and predicting response to immunotherapies. Neoantigen identification and prioritization relies on correctly predicting whether the presenting peptide sequence can successfully induce an immune response. Because most somatic mutations are single-nucleotide variants, changes between wild-type and mutated peptides are typically subtle and require cautious interpretation. A potentially underappreciated variable in neoantigen prediction pipelines is the mutation position within the peptide relative to its anchor positions for the patient's specific MHC molecules. Whereas a subset of peptide positions are presented to the T cell receptor for recognition, others are responsible for anchoring to the MHC, making these positional considerations critical for predicting T cell responses. We computationally predicted anchor positions for different peptide lengths for 328 common HLA alleles and identified unique anchoring patterns among them. Analysis of 923 tumor samples shows that 6 to 38% of neoantigen candidates are potentially misclassified and can be rescued using allele-specific knowledge of anchor positions. A subset of anchor results were orthogonally validated using protein crystallography structures. Representative anchor trends were experimentally validated using peptide-MHC stability assays and competition binding assays. By incorporating our anchor prediction results into neoantigen prediction pipelines, we hope to formalize, streamline, and improve the identification process for relevant clinical studies.
Subject(s)
Antigens, Neoplasm , Neoplasms , Humans , Antigens, Neoplasm/genetics , T-Lymphocytes , Mutation , Peptides/geneticsABSTRACT
Somatic mutations within non-coding regions and even exons may have unidentified regulatory consequences that are often overlooked in analysis workflows. Here we present RegTools ( www.regtools.org ), a computationally efficient, free, and open-source software package designed to integrate somatic variants from genomic data with splice junctions from bulk or single cell transcriptomic data to identify variants that may cause aberrant splicing. We apply RegTools to over 9000 tumor samples with both tumor DNA and RNA sequence data. RegTools discovers 235,778 events where a splice-associated variant significantly increases the splicing of a particular junction, across 158,200 unique variants and 131,212 unique junctions. To characterize these somatic variants and their associated splice isoforms, we annotate them with the Variant Effect Predictor, SpliceAI, and Genotype-Tissue Expression junction counts and compare our results to other tools that integrate genomic and transcriptomic data. While many events are corroborated by the aforementioned tools, the flexibility of RegTools also allows us to identify splice-associated variants in known cancer drivers, such as TP53, CDKN2A, and B2M, and other genes.
Subject(s)
Neoplasms , Transcriptome , Humans , Transcriptome/genetics , Genomics , RNA Splicing/genetics , Genome , Neoplasms/genetics , Alternative Splicing/geneticsABSTRACT
Patients with multiple myeloma (MM) who are treated with lenalidomide rarely develop a secondary B-cell acute lymphoblastic leukemia (B-ALL). The clonal and biological relationship between these sequential malignancies is not yet clear. We identified 17 patients with MM treated with lenalidomide, who subsequently developed B-ALL. Patient samples were evaluated through sequencing, cytogenetics/fluorescence in situ hybridization (FISH), immunohistochemical (IHC) staining, and immunoglobulin heavy chain (IgH) clonality assessment. Samples were assessed for shared mutations and recurrently mutated genes. Through whole exome sequencing and cytogenetics/FISH analysis of 7 paired samples (MM vs matched B-ALL), no mutational overlap between samples was observed. Unique dominant IgH clonotypes between the tumors were observed in 5 paired MM/B-ALL samples. Across all 17 B-ALL samples, 14 (83%) had a TP53 variant detected. Three MM samples with sufficient sequencing depth (>500×) revealed rare cells (average of 0.6% variant allele frequency, or 1.2% of cells) with the same TP53 variant identified in the subsequent B-ALL sample. A lack of mutational overlap between MM and B-ALL samples shows that B-ALL developed as a second malignancy arising from a founding population of cells that likely represented unrelated clonal hematopoiesis caused by a TP53 mutation. The recurrent variants in TP53 in the B-ALL samples suggest a common path for malignant transformation that may be similar to that of TP53-mutant, treatment-related acute myeloid leukemia. The presence of rare cells containing TP53 variants in bone marrow at the initiation of lenalidomide treatment suggests that cellular populations containing TP53 variants expand in the presence of lenalidomide to increase the likelihood of B-ALL development.
Subject(s)
Burkitt Lymphoma , Lenalidomide , Multiple Myeloma , Precursor B-Cell Lymphoblastic Leukemia-Lymphoma , Humans , Bone Marrow/pathology , Burkitt Lymphoma/pathology , Immunoglobulin Heavy Chains/genetics , In Situ Hybridization, Fluorescence , Lenalidomide/adverse effects , Lenalidomide/therapeutic use , Multiple Myeloma/drug therapy , Mutation , Precursor B-Cell Lymphoblastic Leukemia-Lymphoma/pathologyABSTRACT
CIViC (Clinical Interpretation of Variants in Cancer; civicdb.org) is a crowd-sourced, public domain knowledgebase composed of literature-derived evidence characterizing the clinical utility of cancer variants. As clinical sequencing becomes more prevalent in cancer management, the need for cancer variant interpretation has grown beyond the capability of any single institution. CIViC contains peer-reviewed, published literature curated and expertly-moderated into structured data units (Evidence Items) that can be accessed globally and in real time, reducing barriers to clinical variant knowledge sharing. We have extended CIViC's functionality to support emergent variant interpretation guidelines, increase interoperability with other variant resources, and promote widespread dissemination of structured curated data. To support the full breadth of variant interpretation from basic to translational, including integration of somatic and germline variant knowledge and inference of drug response, we have enabled curation of three new Evidence Types (Predisposing, Oncogenic and Functional). The growing CIViC knowledgebase has over 300 contributors and distributes clinically-relevant cancer variant data currently representing >3200 variants in >470 genes from >3100 publications.
Subject(s)
Genetic Variation , Neoplasms , Humans , Neoplasms/genetics , Knowledge Bases , High-Throughput Nucleotide SequencingABSTRACT
Circulating tumor DNA (ctDNA) in peripheral blood has been used to predict prognosis and therapeutic response for triple-negative breast cancer (TNBC) patients. However, previous approaches typically use large comprehensive panels of genes commonly mutated across all breast cancers. Given the reduction in sequencing costs and decreased turnaround times associated with panel generation, the objective of this study was to assess the use of custom micro-panels for tracking disease and predicting clinical outcomes for patients with TNBC. Paired tumor-normal samples from patients with TNBC were obtained at diagnosis (T0) and whole exome sequencing (WES) was performed to identify somatic variants associated with individual tumors. Custom micro-panels of 4-6 variants were created for each individual enrolled in the study. Peripheral blood was obtained at baseline, during Cycle 1 Day 3, at time of surgery, and in 3-6 month intervals after surgery to assess variant allele fraction (VAF) at different timepoints during disease course. The VAF was compared to clinical outcomes to evaluate the ability of custom micro-panels to predict pathological response, disease-free intervals, and patient relapse. A cohort of 50 individuals were evaluated for up to 48 months post-diagnosis of TNBC. In total, there were 33 patients who did not achieve pathological complete response (pCR) and seven patients developed clinical relapse. For all patients who developed clinical relapse and had peripheral blood obtained ≤ 6 months prior to relapse (n = 4), the custom ctDNA micro-panels identified molecular relapse at an average of 4.3 months prior to clinical relapse. The custom ctDNA panel results were moderately associated with pCR such that during disease monitoring, only 11% of patients with pCR had a molecular relapse, whereas 47% of patients without pCR had a molecular relapse (Chi-Square; p-value = 0.10). In this study, we show that a custom micro-panel of 4-6 markers can be effectively used to predict outcomes and monitor remission for patients with TNBC. These custom micro-panels show high sensitivity for detecting molecular relapse in advance of clinical relapse. The use of these panels could improve patient outcomes through early detection of relapse with preemptive intervention prior to symptom onset.
Subject(s)
Circulating Tumor DNA , Triple Negative Breast Neoplasms , Humans , Circulating Tumor DNA/genetics , Triple Negative Breast Neoplasms/diagnosis , Triple Negative Breast Neoplasms/genetics , Biomarkers, Tumor/genetics , Neoplasm Recurrence, Local/pathology , PrognosisABSTRACT
The creation of visualizations to interpret genomics data remains an important aspect of data science within computational biology. The GenVisR Bioconductor package was created to lower the entry point for publication-quality graphics and has remained a popular suite of tools within this domain. GenVisR supports visualizations covering a breadth of topics including functions to produce visual summaries of copy-number alterations, somatic variants, sequence quality metrics, and more. Recently, the GenVisR package has undergone significant updates to increase performance and functionality. To demonstrate the utility of GenVisR, we present protocols for use of the updated Waterfall() function to create a customizable Oncoprint-style plot of the mutational landscape of a tumor cohort. We explain the basics of installation, data import, configuration, plotting, clinical annotation, and customization. A companion online workshop describing the GenVisR library, Waterfall() function, and other genomic visualization tools is available at genviz.org. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Generating a Waterfall() plot from original mutation data Basic Protocol 2: Adding clinical data to a Waterfall() plot Basic Protocol 3: Customizing mutation burden in Waterfall() plots Basic Protocol 4: Brief exploration of customizable options Support Protocol 1: Installing GenVisR.
Subject(s)
Neoplasms , Software , Computational Biology , Genome , Genomics , Humans , Neoplasms/geneticsABSTRACT
Allogeneic cancer vaccines are designed to induce antitumor immune responses with the goal of impacting tumor growth. Typical allogeneic cancer vaccines are produced by expansion of established cancer cell lines, transfection with vectors encoding immunostimulatory cytokines, and lethal irradiation. More than 100 clinical trials have investigated the clinical benefit of allogeneic cancer vaccines in various cancer types. Results show limited therapeutic benefit in clinical trials and currently there are no FDA approved allogeneic cancer vaccines. We used recently developed bioinformatics tools including the pVAC-seq suite of software tools to analyze DNA/RNA sequencing data from the TCGA to examine the repertoire of antigens presented by a typical allogeneic cancer vaccine, and to simulate allogeneic cancer vaccine clinical trials. Specifically, for each simulated clinical trial we modeled the repertoire of antigens presented by allogeneic cancer vaccines consisting of three hypothetical cancer cell lines to 30 patients with the same cancer type. Simulations were repeated ten times for each cancer type. Each tumor sample in the vaccine and the vaccine recipient was subjected to HLA typing, differential expression analyses for tumor associated antigens (TAAs), germline variant calling, and neoantigen prediction. These analyses provided a robust, quantitative comparison between potentially beneficial TAAs and neoantigens versus distracting antigens present in the allogeneic cancer vaccines. We observe that distracting antigens greatly outnumber shared TAAs and neoantigens, providing one potential explanation for the lack of observed responses to allogeneic cancer vaccines. This analysis provides additional rationale for the redirection of efforts towards a personalized cancer vaccine approach.
Subject(s)
Cancer Vaccines , Hematopoietic Stem Cell Transplantation , Neoplasms , Humans , Epitopes , Neoplasms/therapy , Antigens, Neoplasm/geneticsABSTRACT
The Drug-Gene Interaction Database (DGIdb, www.dgidb.org) is a web resource that provides information on drug-gene interactions and druggable genes from publications, databases, and other web-based sources. Drug, gene, and interaction data are normalized and merged into conceptual groups. The information contained in this resource is available to users through a straightforward search interface, an application programming interface (API), and TSV data downloads. DGIdb 4.0 is the latest major version release of this database. A primary focus of this update was integration with crowdsourced efforts, leveraging the Drug Target Commons for community-contributed interaction data, Wikidata to facilitate term normalization, and export to NDEx for drug-gene interaction network representations. Seven new sources have been added since the last major version release, bringing the total number of sources included to 41. Of the previously aggregated sources, 15 have been updated. DGIdb 4.0 also includes improvements to the process of drug normalization and grouping of imported sources. Other notable updates include the introduction of a more sophisticated Query Score for interaction search results, an updated Interaction Score, the inclusion of interaction directionality, and several additional improvements to search features, data releases, licensing documentation and the application framework.
Subject(s)
Crowdsourcing , Databases, Factual , Databases, Genetic , Drugs, Investigational/pharmacology , Genome, Human/drug effects , Prescription Drugs/pharmacology , Databases, Chemical , Drugs, Investigational/chemistry , Genotype , Humans , Internet , Knowledge Bases , Phenotype , Prescription Drugs/chemistry , SoftwareABSTRACT
PURPOSE: Precision oncology depends on the matching of tumor variants to relevant knowledge describing the clinical significance of those variants. We recently developed the Clinical Interpretations for Variants in Cancer (CIViC; civicdb.org) crowd-sourced, expert-moderated, and open-access knowledgebase. CIViC provides a structured framework for evaluating genomic variants of various types (eg, fusions, single-nucleotide variants) for their therapeutic, prognostic, predisposing, diagnostic, or functional utility. CIViC has a documented application programming interface for accessing CIViC records: assertions, evidence, variants, and genes. Third-party tools that analyze or access the contents of this knowledgebase programmatically must leverage this application programming interface, often reimplementing redundant functionality in the pursuit of common analysis tasks that are beyond the scope of the CIViC Web application. METHODS: To address this limitation, we developed CIViCpy (civicpy.org), a software development kit for extracting and analyzing the contents of the CIViC knowledgebase. CIViCpy enables users to query CIViC content as dynamic objects in Python. We assess the viability of CIViCpy as a tool for advancing individualized patient care by using it to systematically match CIViC evidence to observed variants in patient cancer samples. RESULTS: We used CIViCpy to evaluate variants from 59,437 sequenced tumors of the American Association for Cancer Research Project GENIE data set. We demonstrate that CIViCpy enables annotation of > 1,200 variants per second, resulting in precise variant matches to CIViC level A (professional guideline) or B (clinical trial) evidence for 38.6% of tumors. CONCLUSION: The clinical interpretation of genomic variants in cancers requires high-throughput tools for interoperability and analysis of variant interpretation knowledge. These needs are met by CIViCpy, a software development kit for downstream applications and rapid analysis. CIViCpy is fully documented, open-source, and available free online.
Subject(s)
Data Mining/methods , High-Throughput Nucleotide Sequencing/methods , Mutation , Neoplasm Proteins/genetics , Neoplasms/genetics , Software , Databases, Genetic/standards , Humans , Knowledge Bases , Neoplasms/diagnosis , Neoplasms/therapy , Precision Medicine/standards , User-Computer InterfaceABSTRACT
PURPOSE: Clinical targeted sequencing panels are important for identifying actionable variants for patients with cancer; however, existing approaches do not provide transparent and rationally designed clinical panels to accommodate the rapidly growing knowledge within oncology. MATERIALS AND METHODS: We used the Clinical Interpretations of Variants in Cancer (CIViC) database to develop an Open-Sourced CIViC Annotation Pipeline (OpenCAP). OpenCAP provides methods to identify variants within the CIViC database, build probes for variant capture, use probes on prospective samples, and link somatic variants to CIViC clinical relevance statements. OpenCAP was tested using a single-molecule molecular inversion probe (smMIP) capture design on 27 cancer samples from 5 tumor types. In total, 2,027 smMIPs were designed to target 111 eligible CIViC variants (61.5 kb of genomic space). RESULTS: When compared with orthogonal sequencing, CIViC smMIP sequencing demonstrated a 95% sensitivity for variant detection (n = 61 of 64 variants). Variant allele frequencies for variants identified on both sequencing platforms were highly concordant (Pearson's r = 0.885; n = 61 variants). Moreover, for individuals with paired tumor and normal samples (n = 12), 182 clinically relevant variants missed by orthogonal sequencing were discovered by CIViC smMIP sequencing. CONCLUSION: The OpenCAP design paradigm demonstrates the utility of an open-source and open-access database built on attendant community contributions with peer-reviewed interpretations. Use of a public repository for variant identification, probe development, and variant interpretation provides a transparent approach to build dynamic next-generation sequencing-based oncology panels.
Subject(s)
Biomarkers, Tumor/genetics , Computational Biology/methods , DNA Probes/genetics , High-Throughput Nucleotide Sequencing/standards , Molecular Sequence Annotation/methods , Neoplasms/genetics , DNA Mutational Analysis/methods , Databases, Genetic , Genetic Variation , High-Throughput Nucleotide Sequencing/methods , Humans , Molecular Sequence Annotation/standards , Molecular Targeted Therapy , Neoplasms/diagnosis , ROC Curve , Software DesignABSTRACT
PURPOSE: Following automated variant calling, manual review of aligned read sequences is required to identify a high-quality list of somatic variants. Despite widespread use in analyzing sequence data, methods to standardize manual review have not been described, resulting in high inter- and intralab variability. METHODS: This manual review standard operating procedure (SOP) consists of methods to annotate variants with four different calls and 19 tags. The calls indicate a reviewer's confidence in each variant and the tags indicate commonly observed sequencing patterns and artifacts that inform the manual review call. Four individuals were asked to classify variants prior to, and after, reading the SOP and accuracy was assessed by comparing reviewer calls with orthogonal validation sequencing. RESULTS: After reading the SOP, average accuracy in somatic variant identification increased by 16.7% (p value = 0.0298) and average interreviewer agreement increased by 12.7% (p value < 0.001). Manual review conducted after reading the SOP did not significantly increase reviewer time. CONCLUSION: This SOP supports and enhances manual somatic variant detection by improving reviewer accuracy while reducing the interreviewer variability for variant calling and annotation.
Subject(s)
High-Throughput Nucleotide Sequencing/standards , Mutation/genetics , Neoplasms/genetics , Software , Algorithms , Humans , Neoplasms/pathology , Polymorphism, Single Nucleotide/genetics , Sequence AlignmentABSTRACT
The drug-gene interaction database (DGIdb, www.dgidb.org) consolidates, organizes and presents drug-gene interactions and gene druggability information from papers, databases and web resources. DGIdb normalizes content from 30 disparate sources and allows for user-friendly advanced browsing, searching and filtering for ease of access through an intuitive web user interface, application programming interface (API) and public cloud-based server image. DGIdb v3.0 represents a major update of the database. Nine of the previously included 24 sources were updated. Six new resources were added, bringing the total number of sources to 30. These updates and additions of sources have cumulatively resulted in 56 309 interaction claims. This has also substantially expanded the comprehensive catalogue of druggable genes and anti-neoplastic drug-gene interactions included in the DGIdb. Along with these content updates, v3.0 has received a major overhaul of its codebase, including an updated user interface, preset interaction search filters, consolidation of interaction information into interaction groups, greatly improved search response times and upgrading the underlying web application framework. In addition, the expanded API features new endpoints which allow users to extract more detailed information about queried drugs, genes and drug-gene interactions, including listings of PubMed IDs, interaction type and other interaction metadata.
Subject(s)
Databases, Pharmaceutical , Genes/drug effects , Antineoplastic Agents , User-Computer InterfaceABSTRACT
The Open Regulatory Annotation database (ORegAnno) is a resource for curated regulatory annotation. It contains information about regulatory regions, transcription factor binding sites, RNA binding sites, regulatory variants, haplotypes, and other regulatory elements. ORegAnno differentiates itself from other regulatory resources by facilitating crowd-sourced interpretation and annotation of regulatory observations from the literature and highly curated resources. It contains a comprehensive annotation scheme that aims to describe both the elements and outcomes of regulatory events. Moreover, ORegAnno assembles these disparate data sources and annotations into a single, high quality catalogue of curated regulatory information. The current release is an update of the database previously featured in the NAR Database Issue, and now contains 1 948 307 records, across 18 species, with a combined coverage of 334 215 080 bp. Complete records, annotation, and other associated data are available for browsing and download at http://www.oreganno.org/.