Search

Supplementary Methods

Table S1. Canine cohort metadata Table S2. Canine library preparation, amplification, and complexity Table S3. COSMIC Cancer Gene Census genes Table S4. Summary of canine RNA-seq data Table S5. RNA-seq validation Table S6. RNA-seq survey Table S7. Enrichment of protein domains in hemangiosarcoma and angiosarcoma Table S8. Canine somatic copy number aberrations by oaCGH Table S9. Comparison of top SCNAs in human Angiosarcoma Project data with canine oaCGH data Figure S1 - Canine exome sequencing workflow Figure S2 - Canine SMGs called in normal vs. overamplified libraries Figure S3 - Mutational landscape of entire canine cohort Figure S4 - Mutational landscape: FFPE vs Frozen Figure S5 - Mutational landscape: Canine Tumor Location Figure S6 - FFPE vs Frozen exome SCNA data clustering

Supplemental Table 7

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.

Angiosarcoma is a highly aggressive cancer of blood vessel-forming cells with high fatality and few effective treatment options. It is both rare and heterogenous, making large, well powered genomic studies nearly impossible. In dogs, angiosarcoma is common, with breeds like the golden retriever carrying heritable genetic factors that put them at very high risk. If the clinical similarity of canine 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 angiosarcoma via whole exome sequencing (47 golden retriever angiosarcomas) and RNA sequencing (74 angiosarcomas from multiple breeds). The predominant mutational signature was the age-associated deamination of cytosine to thymine, and somatic coding mutations occurred most frequently in the tumor suppressorTP53(59.6% of cases) as well as two genes in the PI3K pathway: the oncogenePIK3CA(29.8%) and its regulatory subunitPIK3R1(8.5%). We compared the canine data to human data recently released by The Angiosarcoma Project, and found the same genes and many of the same pathways significantly enriched for somatic mutations, most notably protein kinases, glycoproteins, fibronectin Type III domains, EGF-like domains, and cell adhesion proteins such as cadherins. As in human angiosarcoma,CDKN2A/Bwas recurrently deleted andVEGFA, KDR, and KITrecurrently gained. Canine angiosarcoma closely models human angiosarcoma on a genomic level, and is a powerful tool for investigating the pathogenesis of this devastating disease.

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.

Canine degenerative myelopathy (DM) is a naturally occurring neurodegenerative disease with similarities to some forms of amyotrophic lateral sclerosis (ALS). Most dogs that develop DM are homozygous for a common superoxide dismutase 1 gene (SOD1) mutation. However, not all dogs homozygous for this mutation develop disease. We performed a genome-wide association analysis in the Pembroke Welsh Corgi (PWC) breed comparing DM-affected and -unaffected dogs homozygous for the SOD1 mutation. The analysis revealed a modifier locus on canine chromosome 25. A haplotype within the SP110 nuclear body protein (SP110) was present in 40% of affected compared with 4% of unaffected dogs (P = 1.5 × 10(-5)), and was associated with increased probability of developing DM (P = 4.8 × 10(-6)) and earlier onset of disease (P = 1.7 × 10(-5)). SP110 is a nuclear body protein involved in the regulation of gene transcription. Our findings suggest that variations in SP110-mediated gene transcription may underlie, at least in part, the variability in risk for developing DM among PWCs that are homozygous for the disease-related SOD1 mutation. Further studies are warranted to clarify the effect of this modifier across dog breeds.

Dogs, with their breed-determined limited genetic background, are great models of human disease including cancer. Canine B-cell lymphoma and hemangiosarcoma are both malignancies of the hematologic system that are clinically and histologically similar to human B-cell non-Hodgkin lymphoma and angiosarcoma, respectively. Golden retrievers in the US show significantly elevated lifetime risk for both B-cell lymphoma (6%) and hemangiosarcoma (20%). We conducted genome-wide association studies for hemangiosarcoma and B-cell lymphoma, identifying two shared predisposing loci. The two associated loci are located on chromosome 5, and together contribute ~20% of the risk of developing these cancers. Genome-wide p-values for the top SNP of each locus are 4.6×10-7 and 2.7×10-6, respectively. Whole genome resequencing of nine cases and controls followed by genotyping and detailed analysis identified three shared and one B-cell lymphoma specific risk haplotypes within the two loci, but no coding changes were associated with the risk haplotypes. Gene expression analysis of B-cell lymphoma tumors revealed that carrying the risk haplotypes at the first locus is associated with down-regulation of several nearby genes including the proximal gene TRPC6, a transient receptor Ca2+-channel involved in T-cell activation, among other functions. The shared risk haplotype in the second locus overlaps thevesicle transport and release gene STX8. Carrying the shared risk haplotype is associated with gene expression changes of 100 genes enriched for pathways involved in immune cell activation. Thus, the predisposing germ-line mutations in B-cell lymphoma and hemangiosarcoma appear to be regulatory, and affect pathways involved in T-cell mediated immune response in the tumor. This suggests that the interaction between the immune system and malignant cells plays a common role in the tumorigenesis of these relatively different cancers., Golden Retriever Foundation, Morris Animal Foundation (Grant D10CA-501), National Institutes of Health (U.S.) (P30CA077598 (MCC Core)), Land of PureGold Foundation, Inc., Uppsala University, American Kennel Club Canine Health Foundation, Inc. (Grant 422), American Kennel Club Canine Health Foundation, Inc. (Grant 615), American Kennel Club Canine Health Foundation, Inc. (Grant 2254), American Kennel Club Canine Health Foundation, Inc. (Grant 1131), American Kennel Club Canine Health Foundation, Inc. (Grant 593), American Kennel Club Canine Health Foundation, Inc. (Grant 1168), American Kennel Club Canine Health Foundation, Inc. (Grant 1169), National Institutes of Health (U.S.) (RO1CA112211), American Kennel Club Canine Health Foundation, Inc. (Grant 1147), American Kennel Club Canine Health Foundation, Inc. (Starlight Fund), National Institutes of Health (U.S.) (NIH Short-term Training Grant (T32 RR18267)), Swedish Medical Research Council, University of Minnesota, United States. Army Medical Research and Materiel Command (Training Grant (W81XWH-06-1-064)), European Research Council (ERC Young Investigator Award), European Science Foundation (European Young Investigator Award Program)