Skip to main content
  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

  • Register
  • Log in
  • My Cart
Advertisement

Main menu

  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
    • Reviewing
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Rapid Impact Archive
    • Meeting Abstracts
    • Collections
      • Metabolism Collection
      • Editors' Picks
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • OnlineFirst
    • Editors' Picks
    • Citation
    • Author/Keyword
  • News
    • Cancer Discovery News
  • AACR Publications
    • Blood Cancer Discovery
    • Cancer Discovery
    • Cancer Epidemiology, Biomarkers & Prevention
    • Cancer Immunology Research
    • Cancer Prevention Research
    • Cancer Research
    • Clinical Cancer Research
    • Molecular Cancer Research
    • Molecular Cancer Therapeutics

User menu

  • Register
  • Log in
  • My Cart

Search

  • Advanced search
Molecular Cancer Research
Molecular Cancer Research

Advanced Search

  • Home
  • About
    • The Journal
    • AACR Journals
    • Subscriptions
    • Permissions and Reprints
    • Reviewing
  • Articles
    • OnlineFirst
    • Current Issue
    • Past Issues
    • Rapid Impact Archive
    • Meeting Abstracts
    • Collections
      • Metabolism Collection
      • Editors' Picks
  • For Authors
    • Information for Authors
    • Author Services
    • Best of: Author Profiles
    • Submit
  • Alerts
    • Table of Contents
    • OnlineFirst
    • Editors' Picks
    • Citation
    • Author/Keyword
  • News
    • Cancer Discovery News
Cancer “-omics”

Comparative Genomics Reveals Shared Mutational Landscape in Canine Hemangiosarcoma and Human Angiosarcoma

Kate Megquier, Jason Turner-Maier, Ross Swofford, Jong-Hyuk Kim, Aaron L. Sarver, Chao Wang, Sharadha Sakthikumar, Jeremy Johnson, Michele Koltookian, Mitzi Lewellen, Milcah C. Scott, Ashley J. Schulte, Luke Borst, Noriko Tonomura, Jessica Alfoldi, Corrie Painter, Rachael Thomas, Elinor K. Karlsson, Matthew Breen, Jaime F. Modiano, Ingegerd Elvers and Kerstin Lindblad-Toh
Kate Megquier
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kate Megquier
  • For correspondence: kmegq@broadinstitute.org kersli@broadinstitute.org
Jason Turner-Maier
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ross Swofford
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jong-Hyuk Kim
Department of Veterinary Clinical Sciences, College of Veterinary Medicine, University of Minnesota, St. Paul, Minnesota.Animal Cancer Care and Research Program, University of Minnesota, St. Paul, Minnesota.Masonic Cancer Center, University of Minnesota, Minneapolis, Minnesota.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jong-Hyuk Kim
Aaron L. Sarver
Animal Cancer Care and Research Program, University of Minnesota, St. Paul, Minnesota.Masonic Cancer Center, University of Minnesota, Minneapolis, Minnesota.Institute for Health Informatics, University of Minnesota, Minneapolis, Minnesota.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chao Wang
Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sharadha Sakthikumar
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sharadha Sakthikumar
Jeremy Johnson
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michele Koltookian
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mitzi Lewellen
Department of Veterinary Clinical Sciences, College of Veterinary Medicine, University of Minnesota, St. Paul, Minnesota.Animal Cancer Care and Research Program, University of Minnesota, St. Paul, Minnesota.Masonic Cancer Center, University of Minnesota, Minneapolis, Minnesota.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Milcah C. Scott
Department of Veterinary Clinical Sciences, College of Veterinary Medicine, University of Minnesota, St. Paul, Minnesota.Animal Cancer Care and Research Program, University of Minnesota, St. Paul, Minnesota.Masonic Cancer Center, University of Minnesota, Minneapolis, Minnesota.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ashley J. Schulte
Department of Veterinary Clinical Sciences, College of Veterinary Medicine, University of Minnesota, St. Paul, Minnesota.Animal Cancer Care and Research Program, University of Minnesota, St. Paul, Minnesota.Masonic Cancer Center, University of Minnesota, Minneapolis, Minnesota.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Luke Borst
Department of Clinical Sciences, North Carolina State College of Veterinary Medicine, Raleigh, North Carolina.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Noriko Tonomura
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.Tufts Cummings School of Veterinary Medicine, North Grafton, Massachusetts.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jessica Alfoldi
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Corrie Painter
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.Count Me In, Cambridge, Massachusetts.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rachael Thomas
Department of Molecular Biomedical Sciences, North Carolina State University College of Veterinary Medicine, and Comparative Medicine Institute, Raleigh, North Carolina.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Rachael Thomas
Elinor K. Karlsson
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, Massachusetts.Program in Molecular Medicine, University of Massachusetts Medical School, Worcester, Massachusetts.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Elinor K. Karlsson
Matthew Breen
Department of Molecular Biomedical Sciences, North Carolina State University College of Veterinary Medicine, and Comparative Medicine Institute, Raleigh, North Carolina.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jaime F. Modiano
Department of Veterinary Clinical Sciences, College of Veterinary Medicine, University of Minnesota, St. Paul, Minnesota.Animal Cancer Care and Research Program, University of Minnesota, St. Paul, Minnesota.Masonic Cancer Center, University of Minnesota, Minneapolis, Minnesota.Center for Immunology, University of Minnesota, Minneapolis, Minneapolis.Stem Cell Institute, University of Minnesota, Minneapolis, Minnesota.Institute for Engineering in Medicine, University of Minnesota, Minneapolis, Minnesota.Department of Laboratory Medicine and Pathology, School of Medicine, University of Minnesota, Minneapolis, Minnesota.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jaime F. Modiano
Ingegerd Elvers
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kerstin Lindblad-Toh
Broad Institute of Harvard and MIT, Cambridge, Massachusetts.Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: kmegq@broadinstitute.org kersli@broadinstitute.org
DOI: 10.1158/1541-7786.MCR-19-0221 Published December 2019
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Angiosarcoma is a highly aggressive cancer of blood vessel–forming cells with few effective treatment options and high patient mortality. It is both rare and heterogenous, making large, well-powered genomic studies nearly impossible. Dogs commonly suffer from a similar cancer, called hemangiosarcoma, with breeds like the golden retriever carrying heritable genetic factors that put them at high risk. If the clinical similarity of canine hemangiosarcoma and human angiosarcoma reflects shared genomic etiology, dogs could be a critically needed model for advancing angiosarcoma research. We assessed the genomic landscape of canine hemangiosarcoma via whole-exome sequencing (47 golden retriever hemangiosarcomas) and RNA sequencing (74 hemangiosarcomas from multiple breeds). Somatic coding mutations occurred most frequently in the tumor suppressor TP53 (59.6% of cases) as well as two genes in the PI3K pathway: the oncogene PIK3CA (29.8%) and its regulatory subunit PIK3R1 (8.5%). The predominant mutational signature was the age-associated deamination of cytosine to thymine. As reported in human angiosarcoma, CDKN2A/B was recurrently deleted and VEGFA, KDR, and KIT recurrently gained. We compared the canine data to human data recently released by The Angiosarcoma Project, and found many of the same genes and pathways significantly enriched for somatic mutations, particularly in breast and visceral angiosarcomas. Canine hemangiosarcoma closely models the genomic landscape of human angiosarcoma of the breast and viscera, and is a powerful tool for investigating the pathogenesis of this devastating disease.

Implications: We characterize the genomic landscape of canine hemangiosarcoma and demonstrate its similarity to human angiosarcoma.

Introduction

Angiosarcoma is an aggressive cancer of blood vessel–forming cells, associated with poor survival times (1–3). There is an unmet need for new diagnostics and therapies for these patients. However, the rarity of this cancer (approximately 0.01% of all cancers; refs. 4, 5) has limited large-scale genomic studies so far. Canine hemangiosarcoma is a relevant clinical model for understanding the pathophysiology of human angiosarcoma. The human and canine diseases share many clinical similarities, and hemangiosarcoma is common in dogs, occurring in some breeds (notably the golden retriever) with a frequency up to 20% (6). This means that using canine hemangiosarcoma as a model for human disease would yield sample cohorts of a magnitude inaccessible using human data alone. However, for dogs to be an effective model of this disease in the era of precision medicine, detailed genomic characterization of canine hemangiosarcoma must be undertaken, and the results directly compared with existing and emerging genomic data from human angiosarcoma.

Angiosarcoma can form anywhere in the vasculature. In human patients, it most commonly occurs in the skin of the head, neck, and scalp, the breast, the extremities, and less frequently in the liver, right auricle of the heart, bone, and spleen (7). Prognosis is poor, with metastatic disease occurring in approximately 50% of cases (8), and a median overall survival time of approximately 50 months for local disease and 10 months in metastatic cases (9). Treatment involves surgical resection with wide margins, plus or minus radiotherapy, as well as adjuvant chemotherapy in the metastatic disease setting (3). Many angiosarcomas are initially sensitive to doxorubicin, paclitaxel, or targeted agents, but resistance to these therapies is virtually inevitable (3).

While most cases of angiosarcoma in humans occur without known cause, there are several known risk factors. These tumors can arise secondary to radiotherapy for other cancers or chronic lymphedema (3). Other known risk factors include UV irradiation, given the typical locations of cutaneous angiosarcomas on the head and neck (10–12), as well as occupational exposure to vinyl chloride (13), arsenic exposure (14), and use of anabolic steroids (15). Genetically, angiosarcoma is associated with familial syndromes including Li-Fraumeni syndrome (TP53 mutations; ref. 16) and Klippel–Trenaunay syndrome (PIK3CA mutations; ref. 17). However, these syndromes do not solely cause angiosarcoma (7), and angiosarcomas are not present in the majority of cases.

Canine hemangiosarcoma is the histopathologic equivalent of human angiosarcoma (18), and follows a similar, aggressive clinical course. In dogs, the most common tumor locations are the spleen, right auricle of the heart, liver, and skin or subcutaneous tissue (19, 20). The different anatomic distribution seen in the human and canine diseases is likely due, at least in part, to the lack of secondary cases in dogs. Treatment protocols for dogs with hemangiosarcoma similarly involve wide surgical resection, followed by adjuvant chemotherapy. Survival times are short—a median of 4–6 months after surgical resection with adjuvant chemotherapy, with a 1-year survival rate of approximately 10% (21–23). Biological risk factors have not yet been identified in dogs. Genetically, dog breeds display differential predisposition to specific cancers, indicating that there are heritable risk factors that have become common as a result of inbreeding based on selection or drift. In a previous genome-wide association study, we identified several loci significantly associated with the risk of hemangiosarcoma in the golden retriever (24).

Recent targeted next-generation sequencing of human angiosarcomas has begun to reveal the somatic mutational spectrum of this disease. It has so far proven to be fairly heterogeneous—no pathognomonic mutations or copy number aberrations occur in all cases, and tumors from different primary locations or with different underlying etiologies have genomic differences. TP53 and genes in the MAPK pathway are frequently mutated (25), and mutations in PLCG1 and PTPRB are common, particularly in secondary angiosarcomas (25, 26). Although the PI3K pathway is activated in some human angiosarcomas (27, 28), PIK3CA mutations have so far not been commonly reported in human angiosarcoma studies (26, 29). Genes in the VEGF pathway are frequently gained or amplified, including VEGFA and KDR (25, 28). The tumor suppressor CDKN2A is frequently deleted, while the MYC oncogene is frequently amplified, most commonly in radiation-associated tumors (25, 28).

In dogs, a recent whole-exome sequencing (WES) study of a small cohort of 20 hemangiosarcomas showed that the top recurrently mutated genes were PIK3CA and TP53 (30). Earlier candidate gene studies of canine hemangiosarcomas reported mutations in TP53 (31), PTEN (32), and PDGFRA and PDGFRB (33). An analysis of somatic copy number aberrations in visceral hemangiosarcomas from five breeds found that VEGFA was recurrently gained, while CDKN2A was frequently deleted (34). MYC copy number gain was infrequent, likely reflective of the fact that secondary hemangiosarcomas are not seen in canine patients (34). While these earlier studies provide clues as to the genetic features of canine hemangiosarcoma, there is very little of the genome-wide data needed for comprehensive comparison with human angiosarcoma.

To assess the potential utility of hemangiosarcoma as a model for angiosarcoma at the molecular level, we performed the largest exome-sequencing study of hemangiosarcoma to date, and complemented our exome data with oligonucleotide array comparative genomic hybridization (oaCGH) copy number data and RNA-sequencing (RNA-seq) data in partially overlapping cohorts of hemangiosarcoma cases. We then performed comparative analyses of our results with those released by The Angiosarcoma Project direct-to-patient initiative. In this way, we have created a detailed genomic profile of this cancer in the golden retriever breed, and begun vetting this canine cancer as a comparative model for human angiosarcoma, and potentially other tumors.

Materials and Methods

Additional Materials and Methods can be found in the Supplementary Materials and Methods, and a workflow in Supplementary Fig. S1.

Canine hemangiosarcoma sample collection for exome sequencing

Samples were obtained as part of necessary diagnostic procedures with owner consent. DNA from tumor tissue and whole blood was collected from 47 golden retrievers with visceral hemangiosarcoma. Samples were collected from dogs referred for treatment at the University of Minnesota (UMN) Veterinary Medical Center and samples submitted to the Modiano laboratory for diagnostic assessment and/or use in research (n = 34), from cases seen at the North Carolina State University (NCSU) Veterinary Hospital (n = 8), or from diagnostic samples sent to Antech Diagnostics (n = 5). Cases where date of diagnosis was known (n = 36), were diagnosed between 2000 and 2014. Tumor samples were either frozen (n = 17), or formalin-fixed, paraffin-embedded (n = 30; Supplementary Table S1). Procedures involving animal use were approved by the Institutional Animal Care and Use Committees at the Broad Institute (Cambridge, MA), UMN (Minneapolis, MN), or NCSU (Raleigh, NC).

Sample preparation

DNA was extracted and sequencing libraries prepared using the Kapa Hyper Prep Kit (Supplementary Table S2). For 66 samples, additional cycles of PCR were required to obtain sufficient DNA for exome capture, while libraries from 28 samples were reconstructed from source DNA using the standard protocol (see Supplementary Information).

Exome capture

The Roche Nimblegen SeqCap-EZ capture canine exome (120705_CF3_Uppsala_Broad_EZ_HX1) was used for hybrid exome capture, following the manufacturer's protocol.

Sequencing and read alignments

The barcoded exome-captured libraries were multiplexed in pools of eight, and sequenced on the Illumina HiSeq 2500 to a target depth of 60× in the tumor and 30× in the normal, reaching a mean depth of 78× in the tumor and 63× in the normal. Reads were aligned to the CanFam3.1 reference genome using BWA. PCR duplicate reads were flagged using the Picard tool MarkDuplicates. Following the Genome Analysis Toolkit (GATK) Best Practices, we then performed Base Quality Score Recalibration using a set of approximately 19 million known canine germline variant locations.

Somatic variant calling

Somatic mutations were called using MuTect2, using the default settings with the addition of the dontUseSoftClippedBases option, which was added to avoid calling a large number of artifactual indels in our FFPE-preserved samples (see Supplementary Methods). To further refine our set of somatic variant calls, and to avoid artifacts, we also called variants using the GATK4 version of Mutect2, and kept only the consensus of calls which passed in both the GATK3 and GATK4 versions. Both sets of calls were filtered using a “panel of normals” created using all 47 normal (germline) samples, as well as germline samples from previous studies of canine lymphoma, osteosarcoma, and melanoma. Calls were further filtered to exclude oxidation artifacts, and were removed if they overlapped locations with known germline variants. Variants were further filtered if the position had low coverage (defined as a read depth less than 20 in the normal, less than 40 in the tumor, or less than four alternate allele reads in the tumor) or had excessive read depth (greater or equal to mean read depth + 5× SD), to filter out potential alignment errors. Finally, variants with a median read position < 10 were filtered to remove potential artifacts.

Variant annotation

Variants were annotated using the variant effect prediction program SnpEff v4.2. Where multiple effects were predicted for a single variant, the most damaging predicted effect was selected. Coding variants were analyzed using the smg function in Genome MuSiC 0.4 to determine which genes were significantly mutated above the background rate. The merge-concurrent-muts option was applied to count multiple mutations in the same gene within a sample as a single mutation, and a FDR threshold of 0.1 using the convolution test (CT) was applied.

Canine mutational signature discovery

Mutational signatures, considering point mutations and their genomic context, were extracted using a Bayesian nonnegative matrix factorization (NMF) algorithm. The discovered signatures were compared with known signatures in COSMIC (35), as well as those reported in the literature for human angiosarcoma (36). The overall landscape of mutations was plotted for different groupings of samples using the SomaticSignatures package.

RNA sequencing of canine hemangiosarcomas

Seventy-four snap-frozen samples were obtained from 73 dogs with hemangiosarcoma (41 golden retrievers and 32 dogs from 13 other breeds or from mixed breeding; Kim and colleagues, manuscript in preparation). RNA-seq data from 51 of the hemangiosarcoma tissues had been published previously (37, 38). RNA-seq libraries were generated as described previously (37).

Somatic copy number detection in canine hemangiosarcomas

Analysis of somatic copy number aberrations (SCNA) in our canine WES data was limited by the inclusion of both frozen and formalin-fixed samples. The FFPE samples appeared to have an increased number of amplifications and deletions when compared with the frozen samples, possibly as a result of DNA fragmentation during the fixation process (39). While we were able to partially account for this using correction for guanine and cytosine content (40), we were not able to completely remove the bias between the two groups (Supplementary Fig. S6). Hence, we instead used oaCGH data to interrogate copy number changes in a cohort of 69 golden retrievers with hemangiosarcoma tumors, of which 28 were also included in the exome sequencing cohort. SCNAs were called from oligonucleotide array comparative genomic hybridization (oaCGH), as described previously (34), using a approximately 180,000-feature microarray (Agilent Technologies) with approximately 26 kb resolution throughout the dog genome.

Accessing human data from The Angiosarcoma Project for comparative analysis

We compared the somatic mutations and SCNAs present in our cohort of canine hemangiosarcomas to the results reported by The Angiosarcoma Project, a direct-to-patient sequencing project run by the Count Me In initiative (JoinCountMeIn.org). Somatic mutations and SCNAs derived from WES data from 48 samples from 36 patients were downloaded from cBioPortal, and somatic mutations classified as nonsense, missense, splice site, or splice region were kept for analysis. For our analysis, we selected only the annotated primary tumor site from each patient, including 30 tumors—17 breast, five visceral (one each of right atrium, spleen, lung, abdomen, and bladder wall), and eight head, face, neck, and scalp (HFNS) tumors.

Pathway analysis in canine and human data

Pathway analysis was performed on the canine and human mutational data using the DAVID Functional Annotation Tool. For the overall enrichment analyses, all genes with nonsynonymous splice site, or splice region somatic mutations were used. Functional annotation charts were created for Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways using the default options, using the Benjamini–Hochberg method to control the FDR. Analysis was performed independently for the canine data and the human data, mapping genes from both species to Homo sapiens reference to avoid confounding due to differences in gene annotation between the species. Canine Ensembl gene IDs were mapped to their 1:1 orthologs in the human GRCh37 reference genome annotation using Ensembl's BioMart (41).

Data availability

The WES data, as well as the RNA-seq data that has not been previously published, has been submitted to the NCBI Sequence Read Archive. WES data: BioProject PRJNA552034, BioSamples SAMN12173468 - SAMN12173561. RNA-seq data: BioProject PRJNA562916, BioSamples SAMN12659339 - SAMN12659361.

Results

Simple somatic mutations

TP53 and PIK3CA are commonly mutated in canine hemangiosarcomas.

The seven significantly mutated genes (SMGs) in the canine hemangiosarcomas contained well-known cancer genes, including TP53, as well as two genes in the PI3K pathway (Table 1; Fig. 1). Tumor suppressor TP53 was most frequently mutated (28/47 cases, 59.6%), with all 28 cases carrying at least one mutation affecting the DNA-binding domain (Fig. 2). Oncogene PIK3CA (14/47, 29.8%) and its regulatory subunit PIK3R1 (4/47, 8.5%) were both mutated. Ten of the 14 cases with PIK3CA mutations had a mutation at amino acid position 1047, a hotspot frequently mutated in many types of human cancers (ref. 42; Fig. 2). The remaining four SMGs were ORC1 (4/47, 8.5%), RASA1 (4/47, 8.5%), ARPC1A (3/47, 6.4%), and ENSCAFG00000017407 (2/47, 4.3%), a “one-to-many” ortholog of human ATP5PD with 27 paralogs in the canine genome. Overall, 59 genes that were mutated at least once in the canine dataset are annotated as having likely causal somatic mutations in the COSMIC Cancer Gene Census (43), including TP53, PIK3CA, and PIK3R1 (Supplementary Table S3).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Per-sample annotation of metadata. A, Somatic mutations in the seven significantly mutated genes. B, Somatic mutations in the PI3K gene family by tumor location. “Other PI3K” category includes mutations in PIK3CB, PIK3C2G, PIK3C3, PIK3R1, and PIK3R5). SNV, single-nucleotide variant; Indel, insertion/deletion; SCNA, somatic copy number aberration. SCNA frequencies are calculated on the basis of oaCGH data from 69 canine hemangiosarcoma samples.

View this table:
  • View inline
  • View popup
Table 1.

Significantly mutated genes in golden retriever hemangiosarcoma tumors

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Lollipop plots showing the positions of the tumor protein p53 (TP53; A) and phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA; B) mutations in the canine and human data. Canine mutational locations were lifted over to the human hg19 reference genome using the LiftOver tool. Green, missense mutation; Black, truncating mutation.

Variant allele fraction and effect predictions.

For further insight into whether the four SMGs not annotated in COSMIC are likely to be driver mutations, we investigated the variant allele fraction (VAF) of each and whether missense mutations were predicted to be deleterious or tolerated using the SIFT (44) score via Ensembl's Variant Effect Predictor (45). The median VAF for the four RASA1 mutations was 0.12 (range 0.08–0.23). Three mutations were nonsense mutations, with the fourth a missense mutation predicted to be deleterious (SIFT score 0), supporting a causal role for these mutations. Similarly, the median VAF of the three ARPC1A mutations was 0.16 (range 0.12–0.17), with one nonsense mutation, and two missense mutations at the same position, predicted to be deleterious (SIFT score 0). The median VAF for the four ORC1 mutations (all at the same position) was lower, at 0.08, and the change was predicted to be tolerated (SIFT score 0.06). The two mutations in ENSCAFG00000017407 also had a low median VAF (0.05), and were predicted to be tolerated (SIFT score 1).

Comparison with top mutated genes in human angiosarcoma.

The most commonly mutated genes in the human data were TP53 (8/30 tumors, 29%), followed by KDR, LRP2, RYR2, and ABCA13 with mutations in 7 of 30 tumors each, and PIK3CA, FLG, ASXL3, MYH14, and UNC13C, mutated in 6 of 30 tumors each. However, the distribution of these mutations varied by tumor location. In the human data, breast tumors were the only location to carry PIK3CA mutations (n = 6 patients), and were the most common location for KDR mutations. TP53 mutations occurred in all locations, but were more common in HFNS and visceral tumors. KDR was much less frequently mutated in the canine cohort, with only one case having a nonsynonymous mutation, predicted to be tolerated. Two human tumors also had mutations in the canine SMG RASA1. One was a nonsense mutation, the other predicted to be deleterious.

Differential patterns of mutations found in some tumor locations.

The pattern of mutations in the PI3K gene family varied by tumor location in the canine data, as well. In the canine cases, a slightly higher proportion of heart tumors (11/16, 68.8%) compared with splenic tumors (9/22, 40.9%) had PI3K alterations, although this difference was not significant (Pχ2 = 0.09). Liver tumors had slightly fewer overall alterations in the PI3K gene family (3/8, 37.5%), and were less likely to have PIK3CA mutations (1/8, 12.5%), although this was not statistically significant.

Comparison of mutational burden by tumor location in canine and human data.

We found fewer nonsynonymous and splice site/region mutations in the canine tumors than reported in the human tumors overall, while the number of mutations per sample varied significantly by tumor location in human cases. In the canine cohort, there were a median of 23 nonsynonymous coding mutations per sample (range 1–152; Supplementary Table S1), which was closest to the human breast (median = 32, n = 17), and visceral (median = 40, n = 5) tumors. In the dogs, there was a small difference in the total number of nonsynonymous and splice site/region mutations between heart (median = 31, n = 15, outlier with 152 mutations removed) and splenic (n = 22, median = 17.5) tumors (PTukeyHSD = 0.04). In the human data, there was a significant difference in the median number of mutations by tumor location, with head, face, neck, and scalp (HFNS, n = 8) tumors having a higher mutational burden (median = 596.5) than breast or visceral tumors (PANOVA = 3.6 × 10−5; ref. 46).

Comparative pathway analysis highlights similarity between canine and human tumors.

We compared the gene functional annotations enriched in hemangiosarcoma to those enriched in breast and visceral angiosarcoma and those enriched in HFNS angiosarcoma (Table 2). KEGG pathways enriched in canine hemangiosarcoma overlapped with KEGG pathways enriched in both subtypes of angiosarcoma. Pathways enriched in hemangiosarcoma included a larger fraction of those enriched in visceral and breast angiosarcoma (75%, vs. 48% of pathways enriched in HFNS angiosarcoma); however, the total number of pathways shared was greater between hemangiosarcoma and HFNS angiosarcoma (10 pathways). Pathways shared between hemangiosarcoma and visceral and breast angiosarcoma included pathways in cancer (Pcanine = 3.6 × 10−4, Pvisceral/breast = 4.8 × 10−2), Rap1 signaling pathway (Pcanine = 5.8 × 10−3, Pvisceral/breast = 9.0 × 10−3) and central carbon metabolism in cancer (Pcanine = 3.0 × 10−3, Pvisceral/breast = 2.3 × 10−2). Pathways shared between hemangiosarcoma, visceral and breast angiosarcoma, and HFNS angiosarcoma included focal adhesion pathway (Pcanine = 1.4 × 10−3, Pvisceral/breast = 2.2 × 10−2, PHFNS = 8.1 × 10−4), calcium signaling pathway (Pcanine = 2.8 × 10−2, Pvisceral/breast = 5.4 × 10−3, PHFNS = 1.9 × 10−2), and long-term depression (Pcanine = 1.5 × 10−2, Pvisceral/breast = 7.1 × 10−3, PHFNS = 8.6 × 10−3).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2.

KEGG pathways enriched in hemangiosarcoma, visceral and breast angiosarcoma, and HFNS angiosarcoma

In addition, functional annotation enrichment analysis of the genes mutated in hemangiosarcoma, visceral and breast angiosarcoma, and HFNS angiosarcoma found that they were enriched in many of the same protein domains. Shared domains include Fibronectin type III (Pcanine = 9.5 × 10−9, Pvisceral/breast = 6.3 × 10−5, PHFNS = 3.0 × 10−15), Epidermal growth factor-like domain (Pcanine = 1.5 × 10−3, Pvisceral/breast = 5.6 × 10−5, PHFNS = 4.0 × 10−16), and Tyrosine protein kinase active site (Pcanine = 6.9 × 10−4, Pvisceral/breast = 3.8 × 10−3, PHFNS = 2.3 × 10−6). Both visceral and breast and HFNS angiosarcomas were enriched in cadherins (Pvisceral/breast = 4.9 × 10−4, PHFNS = 4.7 × 10−16), HFNS highly so, while hemangiosarcoma was not (Pcanine = 0.14). In addition, HFNS angiosarcoma was enriched in ABC-transporters (P = 2.8 × 10−7) while visceral and breast angiosarcoma and canine hemangiosarcoma were not (Supplementary Table S7).

Pathways and gene families previously reported in the angiosarcoma literature.

Several pathways and gene families have been previously reported to be affected in the human angiosarcoma literature. Mutations in both PLCG1 and PTPRB have been reported to be recurrent (26), particularly in secondary angiosarcoma. In the human data, 12 patients had mutations in eight phospholipase C genes, including PLCG1 (n = 5, Table 3). Seven canine cases (15%) had mutations in PLC genes, including one in PLCG1. Twenty protein tyrosine phosphatases were mutated in 12 samples in human, including 2 patients with PTPRB mutations. While none of the canine tumors had somatic mutations in PTPRB, 8 cases (17%) harbored a mutation in one of nine other protein tyrosine phosphatase genes (Table 3).

View this table:
  • View inline
  • View popup
Table 3.

Overlap of mutated genes between canine hemangiosarcoma and human angiosarcoma tumors in several gene families and pathways

The MAPK pathway has been reported to be frequently affected by mutations in angiosarcoma (25), and we found the visceral and breast angiosarcomas to be enriched in mutations in MAPK pathway genes (P = 0.02, Table 2). We noted mutations in 25 genes in the MAPK pathway in the canine cohort, including the SMGs TP53 and RASA1. In the human data, 67 genes in the MAPK pathway were mutated. Thirteen of these genes were mutated in both the canine and human data (Table 3). Although we found mutations in genes involved in histone methyltransferase/demethylase activity in both species, consistent with earlier human studies (47), we saw no significant enrichment on a pathway level in either the canine or human cohorts (Table 3). Tyrosine protein kinases were enriched in both hemangiosarcoma and angiosarcoma. In the human data, 55 genes in the GO Protein Tyrosine Kinase Activity pathway were mutated, while 22 genes in this pathway were mutated in the canine data, with 12 genes mutated in both human and canine samples (Table 3). Low-density lipoprotein receptors were enriched for mutations in the canine (P = 0.03) and human HFNS tumors (P = 0.05), and LRP2 was one of the most frequently mutated genes in the human dataset (mutated in 7/30 tumors). While mutations in the LDL receptor–related protein family have not been widely reported in human angiosarcoma, they have recently been shown to play a role in a variety of cancers (48). We found that four LDL receptor–related protein genes were mutated in both hemangiosarcoma and angiosarcoma, with an additional gene mutated in angiosarcoma only (Table 3).

Mutational signature of aging is present across all canine tumors.

Analysis of mutational signatures, which capture the mutational landscape of tumors and are shaped by both genetic factors and environmental exposure, revealed similarities between canine hemangiosarcoma and human angiosarcoma, as reported by Thibodeau and colleagues (36). We analyzed the signatures of trinucleotide mutational frequencies present in golden retriever hemangiosarcomas, and found a strong signature of mutations arising through spontaneous deamination of methylated cytosines in CpG islands, corresponding to COSMIC signature SBS1 and generally associated with aging (ref. 35; Fig. 3; Supplementary Figs. S3–S5). In addition, we see a faint signature not described in COSMIC (Fig. 3). The overall mutational landscape was consistent across all canine hemangiosarcoma tumors, regardless of tumor location or tissue preservation method (Supplementary Figs. S2–S4).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Mutational signatures called using Bayesian NMF in the entire canine hemangiosarcoma cohort, showing the count of mutations (y-axis) in each trinucleotide context (x-axis). This analysis reveals the COSMIC SBS1 aging signature (top) and a faint novel signature (bottom).

Canine RNA-seq analysis validates exome mutations and extends findings to additional breeds.

We validated our exome-sequencing somatic mutation calls using RNA-seq data from a partially overlapping set of dogs (n = 74 tumor samples from 73 dogs from 14 breeds, as well as mixed-breed dogs; Supplementary Table S4). Thirteen golden retrievers from the exome cohort with mutations in SMGs had the same tumor location included in the RNA-seq data, with a total of 25 SMG mutations among them. We excluded two frameshift deletion variants from the validation, as they showed evidence of deletion in the pileup data, but also unexpected alternate alleles. Of the remaining 23 variants, we validated 15 (65%) mutations. Of the eight sites which were not validated, one was a nonsense mutation, which may have been subject to nonsense-mediated decay, five had a variant allele fraction (VAF) less than 0.1, and three variants came from samples with tumor purity estimate (≤ 35%) on histologic examination, suggesting that tumor heterogeneity and stromal contamination may be a contributing factor to variants not being replicated (Supplementary Table S5). In the RNA-seq data, four tumor samples from three dogs with variants in the SMGs were from a different tumor site than was included in the exome discovery cohort. In this subset, we were able to detect the mutation discovered in the same individual, but a different tumor location in 4 of 8 (50%) cases. Overall, the group of variants, which we were unable to validate, came from tumor samples with a significantly lower estimated tumor purity on histology (Pttest = 0.04), but this difference was not significant based on tumor purity estimates from the program ESTIMATE (Pttest = 0.4; Supplementary Table S5).

The RNA-seq cohort also confirmed that the TP53 and PI3K pathway mutations were not breed-specific, but were present across different breeds. Excluding the 19 golden retrievers who were also included in the exome cohort, we looked to see whether somatic mutations corresponding to those observed in the exome cohort could be found in other cases, which included Golden Retrievers (n = 22), Portuguese water dogs (PWD, n = 6), German Shepherd dogs (GSD, n = 6), mixed breed dogs (n = 6), boxers (n = 2), Labrador retrievers (Labs, n = 2), Keeshonds (n = 2), and other breeds (n = 8, one each of American Staffordshire Terrier, Bernese Mountain Dog, Bichon Frise, Briard, Bullmastiff, Gordon setter, Parson's Russell Terrier, and Saluki). We found mutations at the same sites as discovered in the exome cohort in TP53 (n = 14, 26%, 5 breeds), PIK3CA (n = 11, 20%, 6 breeds), and PIK3R1 (n = 4, 7%, 4 breeds; Supplementary Table S6).

SCNAs

SCNAs in canine hemangiosarcoma recurrently affect known cancer genes.

We surveyed SCNAs in genes known to be involved in hemangiosarcoma and angiosarcoma. The genes most recurrently affected by DNA copy number aberrations in the oaCGH data were VEGFA, with copy number gain in 19% of cases, KDR, gained in 22%, KIT, gained in 17%, and the tumor suppressor CDKN2A/B, deleted in 22% (Supplementary Table S8). The MYC oncogene had copy number gain in 9% of cases. Copy number aberrations in the top significantly mutated genes were relatively rare (Fig. 1A; Supplementary Table S8).

Copy number gains in KDR and losses in AXIN1 are common in both dogs and humans.

Comparison of SCNAs in the human data and canine oaCGH data revealed recurrent copy number aberrations in known cancer genes KDR and AXIN1 in both species. Copy number gains in KDR occurred in approximately 27% of human samples, and 22% of canine samples. AXIN1 was lost in 20% of human and 22% of canine samples. In addition, the genes PTK6, ARFRP1, and RTEL1 showed copy number loss in 17% of human and canine samples; however, these genes are close together, and also showed copy number gains in a number of canine samples (PTK6, 8.7%; ARFRP1, 13%; RTEL1, 13%; Supplementary Table S9).

SCNA profiles differ among cases with and without PIK3CA mutations.

We examined the SCNA profiles of hemangiosarcoma cases with and without TP53 and PIK3CA mutations in the 28 cases with both exome sequencing and oaCGH data. There were no significant differences in the relative frequency of any given CNA between the CNA profiles of cases with and without TP53 mutations. However, significant differences were detected between cases with PIK3CA mutations and those without (Fig. 4), using a two-tailed Fisher exact test and a minimum differential threshold of 25% between the two groups. A region on chromosome 11 at approximately 22.6 Mb, near the UBE2B (ubiquitin-conjugating enzyme E2B) gene, was deleted in 4 of 11 cases with PIK3CA mutations, and 0 of 17 without (P < 0.016). Similarly, the CDKN2B gene, located distally on chromosome 11 at 41.2 Mb, was deleted in 4 of 11 cases with PIK3CA mutations, and 0 of 17 cases without (P < 0.016). A region on chromosome 24 at 21.2 Mb was gained in 7 of 17 cases without PIK3CA mutation and 0 of 11 with PIK3CA mutation (P < 0.023). This region overlaps the antiapoptotic BCL2L1 gene. Broad copy number gains along the length of chromosome 31 were more frequent in cases without PIK3CA mutations compared with those with PIK3CA mutations.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Comparison of DNA copy number aberration profiles between cases with PIK3CA mutation and those without. Blue, copy number gain; red, copy number loss. (Thomas and colleagues, manuscript in preparation; ref. 34)

Discussion

Detailed molecular profiling of canine hemangiosarcoma has revealed both similarities and differences in the genetic landscape between dogs and humans. In particular, visceral canine hemangiosarcoma showed strong similarities to human angiosarcoma of the viscera and breast. Our findings have important implications for comparative oncology, as the study of canine hemangiosarcoma has the potential to improve our understanding of the pathophysiology of both canine hemangiosarcoma and human angiosarcoma, and to improve treatment and outcomes in both species.

Tumor suppressor TP53 was the top significantly mutated gene in the canine data. The majority of mutations occurred in the DNA-binding domain, likely causing loss of function (Fig. 2). TP53 was also the only significantly mutated gene in the human Angiosarcoma Project data, and it has been frequently reported as mutated in targeted sequencing studies of angiosarcoma (25).

We found that the PI3K pathway was commonly mutated in canine hemangiosarcoma. A total of 23 canine tumors (48.9%) had at least one somatic mutation affecting this gene family (Fig. 1A). Tumors with a mutation in the PI3K family tended to have only one mutation in this family. The PI3K pathway is one of the most commonly altered pathways in cancer, playing an important role in signal transduction leading to cell proliferation, survival, differentiation, and regulation of metabolism and immunity (49, 50). PIK3CA is an oncogene (51) that has been shown to be mutated in human glioblastoma, breast, gastric, colorectal, lung, and endometrial cancers (52). Ten of the 14 PIK3CA mutations in our canine cohort occurred at amino acid position 1047, a mutational hotspot in many human cancers (42). Mutations within this domain have been shown to increase catalytic activity (53).

Of the less frequently mutated SMGs, PIK3R1 is annotated as a likely driver gene by COSMIC. In addition, variant allele fraction and SIFT scores support a potential role for RASA1 and ARPC1A as driver genes. RASA1 is a negative regulator of the RAS and MAPK pathways, and plays an important role in vascular formation (54, 55). Germline RASA1 mutations can cause capillary malformation—arteriovenous malformation syndrome (56). Somatic mutations in this gene have been found in a subset of human basal cell carcinomas, and expression has been correlated with survival in invasive ductal breast carcinomas and hepatocellular carcinomas (57, 58). ARPC1A plays an important role in regulating the actin cytoskeleton, which functions in the migration and invasiveness of pancreatic carcinoma cells (59). There is less support for a causal role for ORC1 and ENSCAFG00000017407; however, it is important to note that the SIFT score does not annotate activating mutations; for example, the hotspot mutations in PIK3CA (known to be drivers in many human cancers) are also predicted to be tolerated.

There were also potentially important differences in somatic mutations between the two species. In the human data, mutations in TP53 and PIK3CA tended to be mutually exclusive, while we did not see this pattern in the canine tumors. In addition, PIK3CA mutations were exclusively found in breast tumors in the small human dataset, while we found them to be common in cardiac and splenic tumors in dogs. Within the canine visceral tumors, we found fewer mutations in PIK3CA in liver tumors. These differences in distribution of PIK3CA mutations by tumor location may be due to genetic heterogeneity of the cancer, with tumors in different locations activating the PI3K pathway at different points, or relying on alterations in different pathways to affect an essential protein downstream. It is also possible that this difference is due to the small number of human visceral angiosarcomas currently sequenced, and the small number of liver tumors in the canine cohort.

Another potentially important difference between the two species is that, while copy number gains in KDR are common in both species, somatic mutations in this gene were seen in over 20% of human tumors, but in only one canine tumor. As the KDR receptor is upstream of the PI3K pathway, it is possible that mutations in either may lead to a similar phenotype. We also saw more copy number gains of VEGFA in our canine cohort (19% vs. 0 in the human data), and as VEGFA is upstream of KDR, this copy number gain may serve a similar role to KDR mutations in the canine tumors. In addition, whole genome sequencing will be necessary to determine whether there are common regulatory mutations affecting this gene in canine hemangiosarcoma.

In the human data, mutation rates were significantly different between different tumor locations, with HFNS angiosarcomas having a much higher mutational burden, as reported by Painter and colleagues (46). It is possible that the higher mutational burden in these tumors is due to UV exposure. It would be interesting to compare these findings to canine cutaneous and subcutaneous hemangiosarcomas, which were not included in this study. In the canine visceral hemangiosarcoma, a slight difference was found between mutational burden in heart and splenic tumors. It is possible that the golden retrievers have a lower overall somatic mutation burden than humans, as they likely have a higher germline risk burden, given the high incidence of hemangiosarcoma in the breed. Further studies in dogs will be necessary to determine whether the lower mutation rate is correlated with breed risk, or if in this case it is an artifact of how mutations were called.

Tumors in both dogs and humans were enriched for mutations in protein tyrosine kinases, which are important regulators of cellular growth and division signals and are commonly mutated in cancers. There were also recurrent mutations in the protein tyrosine phosphatase gene family in both species. This may suggest an alternate mechanism of tyrosine kinase overactivation, as protein tyrosine phosphatases deactivate tyrosine kinase signaling by dephosphorylating proteins in opposition to kinase phosphorylation (60). In addition, phospholipase C proteins play a crucial role in cellular signaling pathways by hydrolyzing phosphatidylinositol 4,5-bisphosphate (PIP2) into the second messengers DAG and IP3, passing on signals from receptor tyrosine kinases (61). The canine hemangiosarcomas were enriched for phospholipase C genes, including PLCG1, and PLCG1 mutations were recurrent in the human tumors.

The shared enriched pathways between canine hemangiosarcomas and human angiosarcomas provide insight into disease pathogenesis for both species. Tyrosine kinase inhibitors have been effective against angiosarcoma in the clinic, but tumor heterogeneity and the development of resistance have limited their long-term utility. They have also shown promise against hemangiosarcoma in vitro (62), but so far have been less promising in the veterinary clinic (63). Future investigation of the interaction between the many affected pathways will help to determine the potential for combination therapy targeting multiple of these pathways or a common downstream effector to combat resistance.

A recent study in 13 radiation-induced and three spontaneous breast angiosarcomas detected the irradiation signature, as well as the aging signature, and a unique C>T signature (36). The median age of patients in this human cohort was 74.5 years (36). Canine hemangiosarcomas looked very similar, in that we primarily saw the aging signature and low levels of a novel C>T signature. This novel signature bore some resemblance to the signature reported in the human angiosarcomas, including higher levels of C>T mutations at C nucleotides flanked by A-A or A-T; however, there were also differences, such as a high number of mutations flanked by T-G in the dogs, and a low number of these in the human cohort. A larger study will be needed to decipher whether this is a novel angiosarcoma-related signature or whether it represents noise. The lack of the irradiation signature was anticipated as our canine data did not include any tumors secondary to radiotherapy.

We examined copy number changes in canine hemangiosarcoma in genes previously reported to be affected in hemangiosarcoma and angiosarcoma. Most common were copy number gain of VEGFA and KDR, and loss of CDKN2A/B. The MYC oncogene, which has been reported amplified in human angiosarcoma and canine hemangiosarcoma, was only rarely gained in our dataset and showed no evidence of high-level amplification. This makes sense, given that it is more common in radiation-induced tumors, which were not present in the canine cohort.

Our data suggest that visceral canine hemangiosarcoma could be developed as a model for primary human angiosarcoma. Detailed molecular characterization of canine hemangiosarcoma revealed many similarities, but also some important differences, between canine hemangiosarcoma and human angiosarcoma. Future work should include analysis of larger cohorts, including recalling all currently available human angiosarcoma data using the same methods, to decipher potential molecular subtypes and to facilitate a more complete comparison between tumors in different locations. An integrated understanding of the interaction between mutations in the many enriched signaling pathways may be useful for determining treatment strategy; for example, the feasibility of combinations of targeted inhibitors or the prevention of convergent resistance. Our data suggest that clinical trials evaluating therapeutic approaches in dogs with this disease might also inform human medicine.

Disclosure of Potential Conflicts of Interest

J.F. Modiano is a partner (paid consultant) for Veterinary Research Associates, reports receiving a commercial research grant from Boston Scientific, and has ownership interest (including patents) in Half Moon Bay Biotechnology, LLC. No potential conflict of interests were disclosed by the other authors.

Authors' Contributions

Conception and design: C. Painter, M. Breen, J.F. Modiano, I. Elvers, K. Lindblad-Toh

Development of methodology: K. Megquier, J. Turner-Maier, R. Swofford, A.L. Sarver, S. Sakthikumar, K. Lindblad-Toh

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): K. Megquier, R. Swofford, J.-H. Kim, J. Johnson, M. Koltookian, M. Lewellen, M.C. Scott, A.J. Schulte, L. Borst, N. Tonomura, C. Painter, R. Thomas, M. Breen, J.F. Modiano

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): K. Megquier, J. Turner-Maier, J.-H. Kim, A.L. Sarver, C. Wang, S. Sakthikumar, L. Borst, R. Thomas, E.K. Karlsson, J.F. Modiano, I. Elvers, K. Lindblad-Toh

Writing, review, and/or revision of the manuscript: K. Megquier, J.-H. Kim, A.L. Sarver, S. Sakthikumar, M.C. Scott, L. Borst, N. Tonomura, R. Thomas, E.K. Karlsson, M. Breen, J.F. Modiano, I. Elvers, K. Lindblad-Toh

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J. Turner-Maier, R. Swofford, J. Johnson, M.C. Scott, J. Alfoldi

Study supervision: E.K. Karlsson, M. Breen, J.F. Modiano, I. Elvers, K. Lindblad-Toh

Acknowledgments

The authors would like to thank all of the dogs and owners who participated in our research, as well as the veterinarians who collected samples. The results included here include the use of data from The Angiosarcoma Project (https://ascproject.org/), a project of Count Me In (https://joincountmein.org/), downloaded March 2019. We would also like to thank Dr. Scott Moroff, Vice-President and Chief Scientific Officer of Antech Diagnostics, for contributing tumor and paired blood samples to our research. This work was funded by American Kennel Club (AKC) Canine Health Foundation (CHF) grants #422 (to J.F. Modiano), 1131 (to J.F. Modiano), and 1889-G (to J.F. Modiano, M. Breen, K. Lindblad-Toh, E.K. Karlsson), NIH grants R03CA191713 (to J.F. Modiano), P30CA077598 (NIH Comprehensive Cancer Center Support Grant to the Masonic Cancer Center, University of Minnesota), R37CA218570 (to E.K. Karlsson, C. Painter), and R24OD018250 (to E.K. Karlsson), Cancerfonden (K. Lindblad-Toh), National Canine Cancer Foundation (NCCF) grants DM06CO-003 (to J.F. Modiano) and JHK15MN-004 (to J.-H. Kim), and Morris Animal Foundation grant D10CA-501 (to J.F. Modiano, M. Breen, K. Lindblad-Toh). K. Lindblad-Toh is the recipient of a Distinguished Professor award from the Swedish Research Council. I. Elvers is supported by a postdoctoral fellowship from the Swedish Medical Research Council, SSMF. J.F. Modiano is supported by the Alvin and June Perlman Chair in Animal Oncology at the University of Minnesota. M. Breen is supported in part by the Oscar J. Fletcher Distinguished Professorship in Comparative Oncology Genetics at NC State University. A.L. Sarver is supported by NCI R50 CA211249.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Footnotes

  • Note: Supplementary data for this article are available at Molecular Cancer Research Online (http://mcr.aacrjournals.org/).

  • Mol Cancer Res 2019;17:2410–21

  • Received February 22, 2019.
  • Revision received July 12, 2019.
  • Accepted September 25, 2019.
  • Published first September 30, 2019.
  • ©2019 American Association for Cancer Research.

References

  1. 1.↵
    1. Penel N,
    2. Marréaud S,
    3. Robin YM,
    4. Hohenberger P
    . Angiosarcoma: state of the art and perspectives. Crit Rev Oncol Hematol 2011;80:257–63.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Abraham JA,
    2. Hornicek FJ,
    3. Kaufman AM,
    4. Harmon DC,
    5. Springfield DS,
    6. Raskin KA,
    7. et al.
    Treatment and outcome of 82 patients with angiosarcoma. Ann Surg Oncol 2007;14:1953–67.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Florou V,
    2. Wilky BA
    . Current and future directions for angiosarcoma therapy. Curr Treat Options Oncol 2018;19:14.
    OpenUrl
  4. 4.↵
    1. Antonescu C
    . Malignant vascular tumors—an update. Mod Pathol 2014;27:S30.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Siegel R,
    2. Naishadham D,
    3. Jemal A
    . Cancer statistics, 2013. CA Cancer J Clin 2013;63:11–30.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Glickman L,
    2. Glickman N,
    3. Thorpe R
    . The Golden Retriever Club of America National Health Survey 1998–1999 [Internet]. Available from: https://www.grca.org/wp-content/uploads/2015/08/healthsurvey1-1.pdf.
  7. 7.↵
    1. Young RJ,
    2. Brown NJ,
    3. Reed MW,
    4. Hughes D,
    5. Woll PJ
    . Angiosarcoma. Lancet Oncol 2010;11:983–91.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Mark RJ,
    2. Poen JC,
    3. Tran LM,
    4. Fu YS,
    5. Juillard GF
    . Angiosarcoma: a report of 67 patients and a review of the literature. Cancer 1996;77:2400–6.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Lahat G,
    2. Dhuka AR,
    3. Hallevi H,
    4. Xiao L,
    5. Zou C,
    6. Smith KD,
    7. et al.
    Angiosarcoma: clinical and molecular insights. Ann Surg 2010;251:1098–106.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Stenbäck F
    . Cellular injury and cell proliferation in skin carcinogenesis by UV light. Oncology 1975;31:61–75.
    OpenUrlPubMed
  11. 11.↵
    1. Cioffi A,
    2. Reichert S,
    3. Antonescu CR,
    4. Maki RG
    . Angiosarcomas and other sarcomas of endothelial origin. Hematol Oncol Clin North Am 2013;27:975–88.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Maddox JC,
    2. Evans HL
    . Angiosarcoma of skin and soft tissue: a study of forty-four cases. Cancer 1981;48:1907–21.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Makk L,
    2. Creech JL,
    3. Whelan JG Jr.,
    4. Johnson MN
    . Liver damage and angiosarcoma in vinyl chloride workers. A systematic detection program. JAMA 1974;230:64–8.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Centeno JA,
    2. Mullick FG,
    3. Martinez L,
    4. Page NP,
    5. Gibb H,
    6. Longfellow D,
    7. et al.
    Pathology related to chronic arsenic exposure. Environ Health Perspect 2002;110:883–6.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Falk H,
    2. Thomas LB,
    3. Popper H,
    4. Ishak KG
    . Hepatic angiosarcoma associated with androgenic-anabolic steroids. Lancet 1979;2:1120–3.
    OpenUrlPubMed
  16. 16.↵
    1. Calvete O,
    2. Martinez P,
    3. Garcia-Pavia P,
    4. Benitez-Buelga C,
    5. Paumard-Hernández B,
    6. Fernandez V,
    7. et al.
    A mutation in the POT1 gene is responsible for cardiac angiosarcoma in TP53-negative Li-Fraumeni-like families. Nat Commun 2015;6:8383.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Ploegmakers MJ,
    2. Pruszczynski M,
    3. De Rooy J,
    4. Kusters B,
    5. Veth RP
    . Angiosarcoma with malignant peripheral nerve sheath tumour developing in a patient with Klippel–Trénaunay–Weber syndrome. Sarcoma 2005;9:137–40.
    OpenUrl
  18. 18.↵
    1. Kim JH,
    2. Graef AJ,
    3. Dickerson EB,
    4. Modiano JF
    . Pathobiology of hemangiosarcoma in dogs: research advances and future perspectives. Vet Sci China 2015;2:388–405.
    OpenUrl
  19. 19.↵
    1. Oksanen A
    . Haemangiosarcoma in dogs. J Comp Pathol 1978;88:585–95.
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Brown NO,
    2. Patnaik AK,
    3. MacEwen EG
    . Canine hemangiosarcoma: retrospective analysis of 104 cases. J Am Vet Med Assoc 1985;186:56–8.
    OpenUrlPubMed
  21. 21.↵
    1. Withrow SJ,
    2. Page RL
    . Withrow and MacEwen's small animal clinical oncology. Elsevier Health Sciences; 2013.
  22. 22.↵
    1. Wendelburg KM,
    2. Price LL,
    3. Burgess KE,
    4. Lyons JA,
    5. Lew FH,
    6. Berg J
    . Survival time of dogs with splenic hemangiosarcoma treated by splenectomy with or without adjuvant chemotherapy: 208 cases (2001-2012). J Am Vet Med Assoc 2015;247:393–403.
    OpenUrlPubMed
  23. 23.↵
    1. Sorenmo KU,
    2. Jeglum KA,
    3. Helfand SC
    . Chemotherapy of canine hemangiosarcoma with doxorubicin and cyclophosphamide. J Vet Intern Med 1993;7:370–6.
    OpenUrlPubMed
  24. 24.↵
    1. Tonomura N,
    2. Elvers I,
    3. Thomas R,
    4. Megquier K,
    5. Turner-Maier J,
    6. Howald C,
    7. et al.
    Genome-wide association study identifies shared risk loci common to two malignancies in golden retrievers. PLoS Genet 2015;11:e1004922.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Murali R,
    2. Chandramohan R,
    3. Möller I,
    4. Scholz SL,
    5. Berger M,
    6. Huberman K,
    7. et al.
    Targeted massively parallel sequencing of angiosarcomas reveals frequent activation of the mitogen activated protein kinase pathway. Oncotarget 2015;6:36041–52.
    OpenUrlCrossRef
  26. 26.↵
    1. Behjati S,
    2. Tarpey PS,
    3. Sheldon H,
    4. Martincorena I,
    5. Van Loo P,
    6. Gundem G,
    7. et al.
    Recurrent PTPRB and PLCG1 mutations in angiosarcoma. Nat Genet 2014;46:376–9.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Italiano A,
    2. Chen CL,
    3. Thomas R,
    4. Breen M,
    5. Bonnet F,
    6. Sevenet N,
    7. et al.
    Alterations of the p53 and PIK3CA/AKT/mTOR pathways in angiosarcomas: a pattern distinct from other sarcomas with complex genomics. Cancer 2012;118:5878–87.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Wagner MJ,
    2. Ravi V,
    3. Menter DG,
    4. Sood AK
    . Endothelial cell malignancies: new insights from the laboratory and clinic. NPJ Precision Oncol 2017;1:11.
    OpenUrl
  29. 29.↵
    1. Je EM,
    2. An CH,
    3. Yoo NJ,
    4. Lee SH
    . Mutational analysis of PIK3CA, JAK2, BRAF, FOXL2, IDH1, AKT1 and EZH2 oncogenes in sarcomas. APMIS 2012;120:635–9.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Wang G,
    2. Wu M,
    3. Maloneyhuss MA,
    4. Wojcik J,
    5. Durham AC,
    6. Mason NJ,
    7. et al.
    Actionable mutations in canine hemangiosarcoma. PLoS One 2017;12:e0188667.
    OpenUrl
  31. 31.↵
    1. Mayr B,
    2. Zwetkoff S,
    3. Schaffner G,
    4. Reifinger M
    . Tumour suppressor gene p53 mutation in a case of haemangiosarcoma of a dog. Acta Vet Hung 2002;50:157–60.
    OpenUrlPubMed
  32. 32.↵
    1. Dickerson EB,
    2. Thomas R,
    3. Fosmire SP,
    4. Lamerato-Kozicki AR,
    5. Bianco SR,
    6. Wojcieszyn JW,
    7. et al.
    Mutations of phosphatase and tensin homolog deleted from chromosome 10 in canine hemangiosarcoma. Vet Pathol 2005;42:618–32.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Abou Asa S,
    2. Mori T,
    3. Maruo K,
    4. Khater A,
    5. El-Sawak A,
    6. Abd el-Aziz E,
    7. et al.
    Analysis of genomic mutation and immunohistochemistry of platelet-derived growth factor receptors in canine vascular tumours. Vet Comp Oncol 2015;13:237–45.
    OpenUrl
  34. 34.↵
    1. Thomas R,
    2. Borst L,
    3. Rotroff D,
    4. Motsinger-Reif A,
    5. Lindblad-Toh K,
    6. Modiano JF,
    7. et al.
    Genomic profiling reveals extensive heterogeneity in somatic DNA copy number aberrations of canine hemangiosarcoma. Chromosome Res 2014;22:305–19.
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Forbes SA,
    2. Beare D,
    3. Gunasekaran P,
    4. Leung K,
    5. Bindal N,
    6. Boutselakis H,
    7. et al.
    COSMIC: exploring the world's knowledge of somatic mutations in human cancer. Nucleic Acids Res 2015;43:D805–11.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Thibodeau BJ,
    2. Lavergne V,
    3. Dekhne N,
    4. Benitez P,
    5. Amin M,
    6. Ahmed S,
    7. et al.
    Mutational landscape of radiation-associated angiosarcoma of the breast. Oncotarget 2018;9:10042–53.
    OpenUrl
  37. 37.↵
    1. Gorden BH,
    2. Kim JH,
    3. Sarver AL,
    4. Frantz AM,
    5. Breen M,
    6. Lindblad-Toh K,
    7. et al.
    Identification of three molecular and functional subtypes in canine hemangiosarcoma through gene expression profiling and progenitor cell characterization. Am J Pathol 2014;184:985–95.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Borgatti A,
    2. Koopmeiners JS,
    3. Sarver AL,
    4. Winter AL,
    5. Stuebner K,
    6. Todhunter D,
    7. et al.
    Safe and effective sarcoma therapy through bispecific targeting of EGFR and uPAR. Mol Cancer Ther 2017;16:956–65.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. Do H,
    2. Dobrovic A
    . Sequence artifacts in DNA from formalin-fixed tissues: causes and strategies for minimization. Clin Chem 2015;61:64–71.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Benjamini Y,
    2. Speed TP
    . Summarizing and correcting the GC content bias in high-throughput sequencing. Nucleic Acids Res 2012;40:e72.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Kinsella RJ,
    2. Kähäri A,
    3. Haider S,
    4. Zamora J,
    5. Proctor G,
    6. Spudich G,
    7. et al.
    Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database 2011;2011:bar030.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Lawrence MS,
    2. Stojanov P,
    3. Mermel CH,
    4. Robinson JT,
    5. Garraway LA,
    6. Golub TR,
    7. et al.
    Discovery and saturation analysis of cancer genes across 21 tumour types. Nature 2014;505:495–501.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Sondka Z,
    2. Bamford S,
    3. Cole CG,
    4. Ward SA,
    5. Dunham I,
    6. Forbes SA
    . The COSMIC Cancer Gene Census: describing genetic dysfunction across all human cancers. Nat Rev Cancer 2018;18:696–705.
    OpenUrl
  44. 44.↵
    1. Ng PC,
    2. Henikoff S
    . SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res 2003;31:3812–4.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. McLaren W,
    2. Gil L,
    3. Hunt SE,
    4. Riat HS,
    5. Ritchie GRS,
    6. Thormann A,
    7. et al.
    The ensembl variant effect predictor. Genome Biol 2016;17:122.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Painter C,
    2. Jain E,
    3. Dunphy M,
    4. Anastasio E,
    5. McGillicuddy M,
    6. Stoddard R,
    7. et al.
    High mutation burden and response to immune checkpoint inhibitors in angiosarcomas of the scalp and face. In: Proceedings of the Fourth CRI-CIMT-EATI-AACR International Cancer Immunotherapy Conference: Translating Science into Survival; Sept 30–Oct 3, 2018; New York, NY. Philadelphia (PA): AACR; 2019. Abstract nr B085.
  47. 47.↵
    1. Kunze K,
    2. Spieker T,
    3. Gamerdinger U,
    4. Nau K,
    5. Berger J,
    6. Dreyer T,
    7. et al.
    A recurrent activating PLCG1 mutation in cardiac angiosarcomas increases apoptosis resistance and invasiveness of endothelial cells. Cancer Res 2014;74:6173–83.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Gonias SL,
    2. Karimi-Mostowfi N,
    3. Murray SS,
    4. Mantuano E,
    5. Gilder AS
    . Expression of LDL receptor-related proteins (LRPs) in common solid malignancies correlates with patient survival. PLoS One 2017;12:e0186649.
    OpenUrlCrossRef
  49. 49.↵
    1. Fruman DA,
    2. Chiu H,
    3. Hopkins BD,
    4. Bagrodia S,
    5. Cantley LC,
    6. Abraham RT
    . The PI3K pathway in human disease. Cell 2017;170:605–35.
    OpenUrlCrossRefPubMed
  50. 50.↵
    1. Cully M,
    2. You H,
    3. Levine AJ,
    4. Mak TW
    . Beyond PTEN mutations: the PI3K pathway as an integrator of multiple inputs during tumorigenesis. Nat Rev Cancer 2006;6:184–92.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Chang HW,
    2. Aoki M,
    3. Fruman D,
    4. Auger KR,
    5. Bellacosa A,
    6. Tsichlis PN,
    7. et al.
    Transformation of chicken cells by the gene encoding the catalytic subunit of PI 3-kinase. Science 1997;276:1848–50.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Samuels Y,
    2. Wang Z,
    3. Bardelli A,
    4. Silliman N,
    5. Ptak J,
    6. Szabo S,
    7. et al.
    High frequency of mutations of the PIK3CA gene in human cancers. Science 2004;304:554.
    OpenUrlFREE Full Text
  53. 53.↵
    1. Hon WC,
    2. Berndt A,
    3. Williams RL
    . Regulation of lipid binding underlies the activation mechanism of class IA PI3-kinases. Oncogene 2012;31:3655–66.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Pamonsinlapatham P,
    2. Hadj-Slimane R,
    3. Lepelletier Y,
    4. Allain B,
    5. Toccafondi M,
    6. Garbay C,
    7. et al.
    p120-Ras GTPase activating protein (RasGAP): a multi-interacting protein in downstream signaling. Biochimie 2009;91:320–8.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Henkemeyer M,
    2. Rossi DJ,
    3. Holmyard DP,
    4. Puri MC,
    5. Mbamalu G,
    6. Harpal K,
    7. et al.
    Vascular system defects and neuronal apoptosis in mice lacking ras GTPase-activating protein. Nature 1995;377:695–701.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Eerola I,
    2. Boon LM,
    3. Mulliken JB,
    4. Burrows PE,
    5. Dompmartin A,
    6. Watanabe S,
    7. et al.
    Capillary malformation-arteriovenous malformation, a new clinical and genetic disorder caused by RASA1 mutations. Am J Hum Genet 2003;73:1240–9.
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. Liu Y,
    2. Liu T,
    3. Sun Q,
    4. Niu M,
    5. Jiang Y,
    6. Pang D
    . Downregulation of Ras GTPase-activating protein 1 is associated with poor survival of breast invasive ductal carcinoma patients. Oncol Rep 2015;33:119–24.
    OpenUrl
  58. 58.↵
    1. Chen YL,
    2. Huang WC,
    3. Yao HL,
    4. Chen PM,
    5. Lin PY,
    6. Feng FY,
    7. et al.
    Down-regulation of RASA1 is associated with poor prognosis in human hepatocellular carcinoma. Anticancer Res 2017;37:781–5.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    1. Laurila E,
    2. Savinainen K,
    3. Kuuselo R,
    4. Karhu R,
    5. Kallioniemi A
    . Characterization of the 7q21-q22 amplicon identifies ARPC1A, a subunit of the Arp2/3 complex, as a regulator of cell migration and invasion in pancreatic cancer. Genes Chromosomes Cancer 2009;48:330–9.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Elson A
    . Stepping out of the shadows: oncogenic and tumor-promoting protein tyrosine phosphatases. Int J Biochem Cell Biol 2018;96:135–47.
    OpenUrl
  61. 61.↵
    1. Kim MJ,
    2. Kim E,
    3. Ryu SH,
    4. Suh PG
    . The mechanism of phospholipase C-γ1 regulation. Exp Mol Med 2000;32:101.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Lyles SE,
    2. Milner RJ,
    3. Kow K,
    4. Salute ME
    . In vitro effects of the tyrosine kinase inhibitor, masitinib mesylate, on canine hemangiosarcoma cell lines. Vet Comp Oncol 2012;10:223–35.
    OpenUrlPubMed
  63. 63.↵
    1. Gardner HL,
    2. London CA,
    3. Portela RA,
    4. Nguyen S,
    5. Rosenberg MP,
    6. Klein MK,
    7. et al.
    Maintenance therapy with toceranib following doxorubicin-based chemotherapy for canine splenic hemangiosarcoma. BMC Vet Res 2015;11:131.
    OpenUrl
View Abstract
PreviousNext
Back to top
Molecular Cancer Research: 17 (12)
December 2019
Volume 17, Issue 12
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Editorial Board (PDF)

Sign up for alerts

View this article with LENS

Open full page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Molecular Cancer Research article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Comparative Genomics Reveals Shared Mutational Landscape in Canine Hemangiosarcoma and Human Angiosarcoma
(Your Name) has forwarded a page to you from Molecular Cancer Research
(Your Name) thought you would be interested in this article in Molecular Cancer Research.
Citation Tools
Comparative Genomics Reveals Shared Mutational Landscape in Canine Hemangiosarcoma and Human Angiosarcoma
Kate Megquier, Jason Turner-Maier, Ross Swofford, Jong-Hyuk Kim, Aaron L. Sarver, Chao Wang, Sharadha Sakthikumar, Jeremy Johnson, Michele Koltookian, Mitzi Lewellen, Milcah C. Scott, Ashley J. Schulte, Luke Borst, Noriko Tonomura, Jessica Alfoldi, Corrie Painter, Rachael Thomas, Elinor K. Karlsson, Matthew Breen, Jaime F. Modiano, Ingegerd Elvers and Kerstin Lindblad-Toh
Mol Cancer Res December 1 2019 (17) (12) 2410-2421; DOI: 10.1158/1541-7786.MCR-19-0221

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Comparative Genomics Reveals Shared Mutational Landscape in Canine Hemangiosarcoma and Human Angiosarcoma
Kate Megquier, Jason Turner-Maier, Ross Swofford, Jong-Hyuk Kim, Aaron L. Sarver, Chao Wang, Sharadha Sakthikumar, Jeremy Johnson, Michele Koltookian, Mitzi Lewellen, Milcah C. Scott, Ashley J. Schulte, Luke Borst, Noriko Tonomura, Jessica Alfoldi, Corrie Painter, Rachael Thomas, Elinor K. Karlsson, Matthew Breen, Jaime F. Modiano, Ingegerd Elvers and Kerstin Lindblad-Toh
Mol Cancer Res December 1 2019 (17) (12) 2410-2421; DOI: 10.1158/1541-7786.MCR-19-0221
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Disclosure of Potential Conflicts of Interest
    • Authors' Contributions
    • Acknowledgments
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF
Advertisement

Related Articles

Cited By...

More in this TOC Section

  • Outcome and Gene Expression in post-NACT HGSC
  • Genetic and Epigenetic Characterization of Somatotropinomas
Show more Cancer “-omics”
  • Home
  • Alerts
  • Feedback
  • Privacy Policy
Facebook  Twitter  LinkedIn  YouTube  RSS

Articles

  • Online First
  • Current Issue
  • Past Issues
  • Rapid Impact Archive
  • Meeting Abstracts

Information for

  • Authors
  • Subscribers
  • Advertisers
  • Librarians
  • Reviewers

About MCR

  • About the Journal
  • Editorial Board
  • Permissions
  • Submit a Manuscript
AACR logo

Copyright © 2019 by the American Association for Cancer Research.

Molecular Cancer Research
eISSN: 1557-3125
ISSN: 1541-7786

Advertisement