Neoantigen Landscape Supports Feasibility of Personalized Cancer Vaccine for Follicular Lymphoma.

Personalized cancer vaccines designed to target neoantigens represent a promising new treatment paradigm in oncology. In contrast to classical idiotype vaccines, we hypothesized that 'polyvalent' vaccines could be engineered for the personalized treatment of follicular lymphoma (FL) using neoantigen discovery by combined whole exome sequencing (WES) and RNA sequencing (RNA-Seq). Fifty-eight tumor samples from 57 patients with FL underwent WES and RNA-Seq. Somatic and B-cell clonotype neoantigens were predicted and filtered to identify high-quality neoantigens. B-cell clonality was determined by alignment of B-cell receptor (BCR) CDR3 regions from RNA-Seq data, grouping at the protein level, and comparison to the BCR repertoire from healthy individuals using RNA-Seq data. An average of 52 somatic mutations per patient (range: 2-172) were identified, and two or more (median: 15) high-quality neoantigens were predicted for 56 of 58 FL samples. The predicted neoantigen peptides were composed of missense mutations (77%), indels (9%), gene fusions (3%), and BCR sequences (11%). Building off of these preclinical analyses, we initiated a pilot clinical trial using personalized neoantigen vaccination combined with PD-1 blockade in patients with relapsed or refractory FL (#NCT03121677). Synthetic long peptide (SLP) vaccines targeting predicted high-quality neoantigens were successfully synthesized for and administered to all four patients enrolled. Initial results demonstrate feasibility, safety, and potential immunologic and clinical responses. Our study suggests that a genomics-driven personalized cancer vaccine strategy is feasible for patients with FL, and this may overcome prior challenges in the field.

Single-cell T-cell receptor repertoire profiling in dogs.

Spontaneous cancers in companion dogs are robust models of human disease. Tracking tumor-specific immune responses in these models requires reagents to perform species-specific single cell T cell receptor sequencing (scTCRseq). scTCRseq and integration with scRNA data have not been demonstrated on companion dogs with cancer. Here, five healthy dogs, two dogs with T cell lymphoma and four dogs with melanoma are selected to demonstrate applicability of scTCRseq in a cancer immunotherapy setting. Single-cell suspensions of PBMCs or lymph node aspirates are profiled using scRNA and dog-specific scTCRseq primers. In total, 77,809 V(D)J-expressing cells are detected, with an average of 3498 (348 - 5,971) unique clonotypes identified per sample. In total, 29/34, 40/40, 22/22 and 9/9 known functional TRAV, TRAJ, TRBV and TRBJ gene segments are observed respectively. Pseudogene or otherwise defective gene segments are also detected supporting re-annotation of several as functional. Healthy dogs exhibit highly diverse repertoires, T cell lymphomas exhibit clonal repertoires, and vaccine-treated melanoma dogs are dominated by a small number of highly abundant clonotypes. scRNA libraries define large clusters of V(D)J-expressing CD8+ and CD4 + T cells. Dominant clonotypes observed in melanoma PBMCs are predominantly CD8 + T cells, with activated phenotypes, suggesting possible anti-tumor T cell populations.

Ultra-Deep Sequencing Reveals the Mutational Landscape of Classical Hodgkin Lymphoma.

The malignant Hodgkin and Reed Sternberg (HRS) cells of classical Hodgkin lymphoma (cHL) are scarce in affected lymph nodes, creating a challenge to detect driver somatic mutations. As an alternative to cell purification techniques, we hypothesized that ultra-deep exome sequencing would allow genomic study of HRS cells, thereby streamlining analysis and avoiding technical pitfalls. To test this, 31 cHL tumor/normal pairs were exome sequenced to approximately 1,000× median depth of coverage. An orthogonal error-corrected sequencing approach verified >95% of the discovered mutations. We identified mutations in genes novel to cHL including: CDH5 and PCDH7, novel stop gain mutations in IL4R, and a novel pattern of recurrent mutations in pathways regulating Hippo signaling. As a further application of our exome sequencing, we attempted to identify expressed somatic single-nucleotide variants (SNV) in single-nuclei RNA sequencing (snRNA-seq) data generated from a patient in our cohort. Our snRNA analysis identified a clear cluster of cells containing a somatic SNV identified in our deep exome data. This cluster has differentially expressed genes that are consistent with genes known to be dysregulated in HRS cells (e.g., PIM1 and PIM3). The cluster also contains cells with an expanded B-cell clonotype further supporting a malignant phenotype. This study provides proof-of-principle that ultra-deep exome sequencing can be utilized to identify recurrent mutations in HRS cells and demonstrates the feasibility of snRNA-seq in the context of cHL. These studies provide the foundation for the further analysis of genomic variants in large cohorts of patients with cHL. SIGNIFICANCE: Our data demonstrate the utility of ultra-deep exome sequencing in uncovering somatic variants in Hodgkin lymphoma, creating new opportunities to define the genes that are recurrently mutated in this disease. We also show for the first time the successful application of snRNA-seq in Hodgkin lymphoma and describe the expression profile of a putative cluster of HRS cells in a single patient.

Endothelial cells are a key target of IFN-g during response to combined PD-1/CTLA-4 ICB treatment in a mouse model of bladder cancer.

To explore mechanisms of response to combined PD-1/CTLA-4 immune checkpoint blockade (ICB) treatment in individual cell types, we generated scRNA-seq using a mouse model of invasive urothelial carcinoma with three conditions: untreated tumor, treated tumor, and tumor treated after CD4+ T cell depletion. After classifying tumor cells based on detection of somatic variants and assigning non-tumor cell types using SingleR, we performed differential expression analysis, overrepresentation analysis, and gene set enrichment analysis (GSEA) within each cell type. GSEA revealed that endothelial cells were enriched for upregulated IFN-g response genes when comparing treated cells to both untreated cells and cells treated after CD4+ T cell depletion. Functional analysis showed that knocking out IFNgR1 in endothelial cells inhibited treatment response. Together, these results indicated that IFN-g signaling in endothelial cells is a key mediator of ICB induced anti-tumor activity.

Mutations associated with progression in follicular lymphoma predict inferior outcomes at diagnosis: Alliance A151303.

Follicular lymphoma (FL) is clinically heterogeneous, with select patients tolerating extended watch-and-wait, whereas others require prompt treatment, suffer progression of disease within 24 months of treatment (POD24), and/or experience aggressive histologic transformation (t-FL). Because our understanding of the relationship between genetic alterations in FL and patient outcomes remains limited, we conducted a clinicogenomic analysis of 370 patients with FL or t-FL (from Cancer and Leukemia Group B/Alliance trials 50402/50701/50803, or real-world cohorts from Washington University School of Medicine, Cleveland Clinic, or University of Miami). FL subsets by grade, stage, watch-and-wait, or POD24 status did not differ by mutation burden, whereas mutation burden was significantly higher in relapsed/refractory (rel/ref) FL and t-FL than in newly diagnosed (dx) FL. Nonetheless, mutation burden in dx FL was not associated with frontline progression-free survival (PFS). CREBBP was the only gene more commonly mutated in FL than in t-FL yet mutated CREBBP was associated with shorter frontline PFS in FL. Mutations in 20 genes were more common in rel/ref FL or t-FL than in dx FL, including 6 significantly mutated genes (SMGs): STAT6, TP53, IGLL5, B2M, SOCS1, and MYD88. We defined a mutations associated with progression (MAP) signature as ≥2 mutations in these 7 genes (6 rel/ref FL or t-FL SMGs plus CREBBP). Patients with dx FL possessing a MAP signature had shorter frontline PFS, revealing a 7-gene set offering insight into FL progression risk potentially more generalizable than the m7-Follicular Lymphoma International Prognostic Index (m7-FLIPI), which had modest prognostic value in our cohort. Future studies are warranted to validate the poor prognosis associated with a MAP signature in dx FL, potentially facilitating novel trials specifically in this high-risk subset of patients.

Genomic landscape and gene expression profiles of feline oral squamous cell carcinoma.

Feline oral squamous cell carcinoma (FOSCC) is a cancer of the squamous cell lining in the oral cavity and represents up to 80% of all oral cancers in cats, with a poor prognosis. We have used whole exome sequencing (WES) and RNA sequencing of the tumor to discover somatic mutations and gene expression changes that may be associated with FOSCC occurrence. FOSCC offers a potential comparative model to study human head and neck squamous cell carcinoma (HNSCC) due to its similar spontaneous formation, and morphological and histological features. In this first study using WES to identify somatic mutations in feline cancer, we have identified tumor-associated gene mutations in six cats with FOSCC and found some overlap with identified recurrently mutated genes observed in HNSCC. Four samples each had mutations in TP53, a common mutation in all cancers, but each was unique. Mutations in other cellular growth control genes were also found such as KAT2B and ARID1A. Enrichment analysis of FOSCC gene expression profiles suggests a molecular similarity to human OSCC as well, including alterations in epithelial to mesenchymal transition and IL6/JAK/STAT pathways. In this preliminary study, we present exome and transcriptome results that further our understanding of FOSCC.

Endothelial cells are a key target of IFN-g during response to combined PD-1/CTLA-4 ICB treatment in a mouse model of bladder cancer.

To explore mechanisms of response to combined PD-1/CTLA-4 immune checkpoint blockade (ICB) treatment in individual cell types, we generated scRNA-seq using a mouse model of invasive urothelial carcinoma with three conditions: untreated tumor, treated tumor, and tumor treated after CD4+ T cell depletion. After classifying tumor cells based on detection of somatic variants and assigning non-tumor cell types using SingleR, we performed differential expression analysis, overrepresentation analysis, and gene set enrichment analysis (GSEA) within each cell type. GSEA revealed that endothelial cells were enriched for upregulated IFN-g response genes when comparing treated cells to both untreated cells and cells treated after CD4+ T cell depletion. Functional analysis showed that knocking out IFNgR1 in endothelial cells inhibited treatment response. Together, these results indicated that IFN-g signaling in endothelial cells is a key mediator of ICB induced anti-tumor activity.

Integrated analysis of genomic and transcriptomic data for the discovery of splice-associated variants in cancer.

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.

Distinct clonal identities of B-ALLs arising after lenolidomide therapy for multiple myeloma.

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.

CIViCdb 2022: evolution of an open-access cancer variant interpretation knowledgebase.

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.

Personalized ctDNA micro-panels can monitor and predict clinical outcomes for patients with triple-negative breast cancer.

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.

Integrative genomic analysis reveals low T-cell infiltration as the primary feature of tobacco use in HPV-positive oropharyngeal cancer.

Although tobacco use is an independent adverse prognostic feature in HPV(+) oropharyngeal squamous cell carcinoma (OPSCC), the biologic features associated with tobacco use have not been systematically investigated. We characterized genomic and immunologic features associated with tobacco use through whole exome sequencing, mRNA hybridization, and immunohistochemical staining in 47 HPV(+) OPSCC tumors. Low expression of transcripts in a T cell-inflamed gene expression profile (TGEP) was associated with tobacco use at diagnosis and lower overall and disease-free survival. Tobacco use was associated with an increased proportion of T > C substitutions and a lower proportion of expected mutational signatures, but not with increases in mutational burden or recurrent oncogenic mutations. Our findings suggest that rather than increased mutational burden, tobacco's primary and clinically relevant association in HPV(+) OPSCC is immunosuppression of the tumor immune microenvironment. Quantitative assays of T cell infiltration merit further study as prognostic markers in HPV(+) OPSCC.

Characterization of the Genomic and Immunologic Diversity of Malignant Brain Tumors through Multisector Analysis.

Despite some success in secondary brain metastases, targeted or immune-based therapies have shown limited efficacy against primary brain malignancies such as glioblastoma (GBM). Although the intratumoral heterogeneity of GBM is implicated in treatment resistance, it remains unclear whether this diversity is observed within brain metastases and to what extent cancer cell-intrinsic heterogeneity sculpts the local immune microenvironment. Here, we profiled the immunogenomic state of 93 spatially distinct regions from 30 malignant brain tumors through whole-exome, RNA, and T-cell receptor sequencing. Our analyses identified differences between primary and secondary malignancies, with gliomas displaying more spatial heterogeneity at the genomic and neoantigen levels. In addition, this spatial diversity was recapitulated in the distribution of T-cell clones in which some gliomas harbored highly expanded but spatially restricted clonotypes. This study defines the immunogenomic landscape across a cohort of malignant brain tumors and contains implications for the design of targeted and immune-based therapies against intracranial malignancies. SIGNIFICANCE: This study describes the impact of spatial heterogeneity on genomic and immunologic characteristics of gliomas and brain metastases. The results suggest that gliomas harbor significantly greater intratumoral heterogeneity of genomic alterations, neoantigens, and T-cell clones than brain metastases, indicating the importance of multisector analysis for clinical or translational studies.This article is highlighted in the In This Issue feature, p. 1.

Phase 1/dose expansion trial of brentuximab vedotin and lenalidomide in relapsed or refractory diffuse large B-cell lymphoma.

New therapies are needed for patients with relapsed/refractory (rel/ref) diffuse large B-cell lymphoma (DLBCL) who do not benefit from or are ineligible for stem cell transplant and chimeric antigen receptor therapy. The CD30-targeted, antibody-drug conjugate brentuximab vedotin (BV) and the immunomodulator lenalidomide (Len) have demonstrated promising activity as single agents in this population. We report the results of a phase 1/dose expansion trial evaluating the combination of BV/Len in rel/ref DLBCL. Thirty-seven patients received BV every 21 days, with Len administered continuously for a maximum of 16 cycles. The maximum tolerated dose of the combination was 1.2 mg/kg BV with 20 mg/d Len. BV/Len was well tolerated with a toxicity profile consistent with their use as single agents. Most patients required granulocyte colony-stimulating factor support because of neutropenia. The overall response rate was 57% (95% CI, 39.6-72.5), complete response rate, 35% (95% CI, 20.7-52.6); median duration of response, 13.1 months; median progression-free survival, 10.2 months (95% CI, 5.5-13.7); and median overall survival, 14.3 months (95% CI, 10.2-35.6). Response rates were highest in patients with CD30+ DLBCL (73%), but they did not differ according to cell of origin (P = .96). NK cell expansion and phenotypic changes in CD8+ T-cell subsets in nonresponders were identified by mass cytometry. BV/Len represents a potential treatment option for patients with rel/ref DLBCL. This combination is being further explored in a phase 3 study (registered on https://clinicaltrials.org as NCT04404283). This trial was registered on https://clinicaltrials.gov as NCT02086604.

Exploring the Genomic Landscape of Cancer Patient Cohorts with GenVisR.

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.

Immunogenomic profiling and pathological response results from a clinical trial of docetaxel and carboplatin in triple-negative breast cancer.

PURPOSE: Patients with triple-negative breast cancer (TNBC) who do not achieve pathological complete response (pCR) following neoadjuvant chemotherapy have a high risk of recurrence and death. Molecular characterization may identify patients unlikely to achieve pCR. This neoadjuvant trial was conducted to determine the pCR rate with docetaxel and carboplatin and to identify molecular alterations and/or immune gene signatures predicting pCR. EXPERIMENTAL DESIGN: Patients with clinical stages II/III TNBC received 6 cycles of docetaxel and carboplatin. The primary objective was to determine if neoadjuvant docetaxel and carboplatin would increase the pCR rate in TNBC compared to historical expectations. We performed whole-exome sequencing (WES) and immune profiling on pre-treatment tumor samples to identify alterations that may predict pCR. Thirteen matching on-treatment samples were also analyzed to assess changes in molecular profiles. RESULTS: Fifty-eight of 127 (45.7%) patients achieved pCR. There was a non-significant trend toward higher mutation burden for patients with residual cancer burden (RCB) 0/I versus RCB II/III (median 80 versus 68 variants, p 0.88). TP53 was the most frequently mutated gene, observed in 85.7% of tumors. EGFR, RB1, RAD51AP2, SDK2, L1CAM, KPRP, PCDHA1, CACNA1S, CFAP58, COL22A1, and COL4A5 mutations were observed almost exclusively in pre-treatment samples from patients who achieved pCR. Seven mutations in PCDHA1 were observed in pre-treatment samples from patients who did not achieve pCR. Several immune gene signatures including IDO1, PD-L1, interferon gamma signaling, CTLA4, cytotoxicity, tumor inflammation signature, inflammatory chemokines, cytotoxic cells, lymphoid, PD-L2, exhausted CD8, Tregs, and immunoproteasome were upregulated in pre-treatment samples from patients who achieved pCR. CONCLUSION: Neoadjuvant docetaxel and carboplatin resulted in a pCR of 45.7%. WES and immune profiling differentiated patients with and without pCR. TRIAL REGISTRATION: Clinical trial information: NCT02124902, Registered 24 April 2014 & NCT02547987, Registered 10 September 2015.

Impact of a 40-Gene Targeted Panel Test on Physician Decision Making for Patients With Acute Myeloid Leukemia.

PURPOSE: Physicians treating hematologic malignancies increasingly order targeted sequencing panels to interrogate recurrently mutated genes. The precise impact of these panels on clinical decision making is not well understood. METHODS: Here, we report our institutional experience with a targeted 40-gene panel (MyeloSeq) that is used to generate a report for both genetic variants and variant allele frequencies for the treating physician (the limit of mutation detection is approximately one AML cell in 50). RESULTS: In total, 346 sequencing reports were generated for 325 patients with suspected hematologic malignancies over an 8-month period (August 2018 to April 2019). To determine the influence of genomic data on clinical care for patients with acute myeloid leukemia (AML), we analyzed 122 consecutive reports from 109 patients diagnosed with AML and surveyed the treating physicians with a standardized questionnaire. The panel was ordered most commonly at diagnosis (61.5%), but was also used to assess response to therapy (22.9%) and to detect suspected relapse (15.6%). The panel was ordered at multiple timepoints during the disease course for 11% of patients. Physicians self-reported that 50 of 114 sequencing reports (44%) influenced clinical care decisions in 44 individual patients. Influences were often nuanced and extended beyond identifying actionable genetic variants with US Food and Drug Administration-approved drugs. CONCLUSION: This study provides insights into how physicians are currently using multigene panels capable of detecting relatively rare AML cells. The most influential way to integrate these tools into clinical practice will be to perform prospective clinical trials that assess patient outcomes in response to genomically driven interventions.

Unraveling the chaotic genomic landscape of primary and metastatic canine appendicular osteosarcoma with current sequencing technologies and bioinformatic approaches.

Osteosarcoma is a rare disease in children but is one of the most common cancers in adult large breed dogs. The mutational landscape of both the primary and pulmonary metastatic tumor in two dogs with appendicular osteosarcoma (OSA) was comprehensively evaluated using an automated whole genome sequencing, exome and RNA-seq pipeline that was adapted for this study for use in dogs. Chromosomal lesions were the most common type of mutation. The mutational landscape varied substantially between dogs but the lesions within the same patient were similar. Copy number neutral loss of heterozygosity in mutant TP53 was the most significant driver mutation and involved a large region in the middle of chromosome 5. Canine and human OSA is characterized by loss of cell cycle checkpoint integrity and DNA damage response pathways. Mutational profiling of individual patients with canine OSA would be recommended prior to targeted therapy, given the heterogeneity seen in our study and previous studies.