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

Table S1. Demographic and pathological features of human patients with angiosarcoma Table S2. Demographic and pathological features of dogs with HSA Table S3. Primary antibodies for immunohistochemical evaluation Table S4. Gene ontology and biological functions of fusion partner genes in human angiosarcomas Table S5. Gene ontology and biological functions of fusion partner genes in canine HSAs Table S6. Detection of fusion genes in canine HSAs by RT-PCR Table S7. Biological functions associated with 490 upregulated genes in human angiosarcomas Table S8. Association of frequency between fusion genes, TP53 and PIK3CA mutations Table S9. Association of fusion genes with dog signalment and tumor histology Table S10. Association of frequency between fusion genes, TP53 and PIK3CA mutations Table S11. 41 differentially expressed genes in PF (n = 4) compared to P (n = 9) group Table S12. 212 differentially expressed genes in PF (n = 4) compared to PK (n = 6) group Table S13. 377 differentially expressed genes in PF (n = 4) compared to PFK (n = 4) group Table S14. Immunohistochemical features of human angiosarcomas Table S15. Immunohistochemical features of canine HSAs

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.

Sporadic angiosarcomas are aggressive vascular sarcomas whose rarity and genomic complexity present significant obstacles in deciphering the pathogenic significance of individual genetic alterations. Numerous fusion genes have been identified across multiple types of cancers, but their existence and significance remain unclear in sporadic angiosarcomas. In this study, we leveraged RNA-sequencing data from 13 human angiosarcomas and 76 spontaneous canine hemangiosarcomas to identify fusion genes associated with spontaneous vascular malignancies. Ten novel protein-coding fusion genes, including TEX2-PECAM1 and ATP8A2-FLT1, were identified in seven of the 13 human tumors, with two tumors showing mutations of TP53. HRAS and NRAS mutations were found in angiosarcomas without fusions or TP53 mutations. We found 15 novel protein-coding fusion genes including MYO16-PTK2, GABRA3-FLT1, and AKT3-XPNPEP1 in 11 of the 76 canine hemangiosarcomas; these fusion genes were seen exclusively in tumors of the angiogenic molecular subtype that contained recurrent mutations in TP53, PIK3CA, PIK3R1, and NRAS. In particular, fusion genes and mutations of TP53 cooccurred in tumors with higher frequency than expected by random chance, and they enriched gene signatures predicting activation of angiogenic pathways. Comparative transcriptomic analysis of human angiosarcomas and canine hemangiosarcomas identified shared molecular signatures associated with activation of PI3K/AKT/mTOR pathways. Our data suggest that genome instability induced by TP53 mutations might create a predisposition for fusion events that may contribute to tumor progression by promoting selection and/or enhancing fitness through activation of convergent angiogenic pathways in this vascular malignancy.Implications:This study shows that, while drive events of malignant vasoformative tumors of humans and dogs include diverse mutations and stochastic rearrangements that create novel fusion genes, convergent transcriptional programs govern the highly conserved morphologic organization and biological behavior of these tumors in both species.

S1. Representative putative inter-chromosomal fusion transcripts in human AS and canine HSA. SCLT1-NIPBL fusion in human angiosarcoma. S2. Fusion genes and their association with tumor content and sequencing depth in canine HSA tissues. S3. A, Primers designed for amplification of SCLT1-NIPBL fusion transcripts identified in human angiosarcoma. B, SCLT1-NIPBL fusion gene amplification by quantitative RT-PCR. Expression value represents 2(-â^†Ct)(Fusion gene/GAPDH). C, Primers designed for amplification of four fusion transcripts identified in canine hemangiosarcomas. D, Fusion gene amplification was confirmed in canine tumors by quantitative RT-PCR. E, Electrophoresis (1% agarose gel) of DNA product by RTPCR was done. S4. Genetic alterations of fusion partner genes identified in human AS. S5. Cell type enrichment analysis in ASs and TCGA sarcomas. S6. Three distinct molecular groups of canine HSAs. S7. Kaplan-Meier probability of survival for dogs with angiogenic and inflammatory HSA. S8. Association of the fusion events with demographic and histopathologic features in canine HSA. S9. Gene pathways associated with commonly enriched genes in tumors with fusion genes and TP53 mutations in canine HSA. S10. Violin plots display IHC scores of p53, p-p53 (Ser15), and p-p53 (Ser20) in human AS and canine HSA. S11. Violin plots present IHC scores of AKT and p-AKT (Thr308) in human AS and canine HSA. S12. IHC scores of mTOR and p-mTOR (Ser2448) in human AS.

Osteosarcoma is a debilitating bone cancer that affects humans, especially children and adolescents. A homologous form of osteosarcoma spontaneously occurs in dogs, and its differential incidence observed across breeds allows for the investigation of tumor mutations in the context of multiple genetic backgrounds. Using whole-exome sequencing and dogs from three susceptible breeds (22 golden retrievers, 21 Rottweilers, and 23 greyhounds), we found that osteosarcoma tumors show a high frequency of somatic copy-number alterations (SCNA), affecting key oncogenes and tumor-suppressor genes. The across-breed results are similar to what has been observed for human osteosarcoma, but the disease frequency and somatic mutation counts vary in the three breeds. For all breeds, three mutational signatures (one of which has not been previously reported) and 11 significantly mutated genes were identified. TP53 was the most frequently altered gene (83% of dogs have either mutations or SCNA in TP53), recapitulating observations in human osteosarcoma. The second most frequently mutated gene, histone methyltransferase SETD2, has known roles in multiple cancers, but has not previously been strongly implicated in osteosarcoma. This study points to the likely importance of histone modifications in osteosarcoma and highlights the strong genetic similarities between human and dog osteosarcoma, suggesting that canine osteosarcoma may serve as an excellent model for developing treatment strategies in both species.Significance: Canine osteosarcoma genomics identify SETD2 as a possible oncogenic driver of osteosarcoma, and findings establish the canine model as a useful comparative model for the corresponding human disease. Cancer Res; 78(13); 3421–31. ©2018 AACR.

File contains list of GISTIC Defined Recurrent SCNA genes that overlap between Humans and Dogs.

Supplementary Figure 1 - Workflow for Somatic Point and Indel Mutation detection. Supplementary Figure 2 - Heatmap counts of samples with recurrent focal SCNAs per chromosome. Supplementary Figure 3 - oaCGH profiles penetrance plots and GISTIC analysis. Supplementary Figure 4 - Venn diagrams of shared and unique somatic point mutations and somatic indel mutations across the three breeds. Supplementary Figure 5 - Scatter plot shows age of diagnosis is positively correlated to the mutation counts across all breeds. Supplementary Figure 6 - Mutational signature distribution per breed shows enrichment for the novel signature in Golden Retriever. Supplementary Figure 7 - Overview and functional classes for the mutations in histone modifier genes across the three breeds. Supplementary Figure 8 - Overview of germ-line candidate mutation across three breeds and their functional impacts.

File contains genes that have focal somatic copy number deletions.

File contains supplementary tables and figures for the main article. Supplementary Table 1 - Source for tools and databases used in the methods for analysis. Supplementary Table 2 - Summary of Somatic Point and Indel Mutation across the three breeds. Supplementary Table 3 - Summary of GISTIC-defined recurrent large-scale SCNAs and gene frequency. Supplementary Table 4 - List of GISTIC classified amplification genes. See attached excel spreadsheet. Supplementary Table 5 - List of GISTIC classified deletion genes. See attached excel spreadsheet. Supplementary Table 6 - List of GISTIC genes and their predicted roles in cancer. Supplementary Table 7 - Significantly over-represented pathways for GISTIC- and MuSiC-defined genes. Supplementary Table 8 - Significantly mutated genes as identified by MuSiC and the proportion of samples per breed that show mutation in them. Supplementary Table 9 - Putative driver genes as detected by MutSigCV algorithm Supplementary Table 10 - Putative driver genes assessed by VEP high-impact mutations and their variant allele frequency. Supplementary Table 11 - Mutations and their effects on p53 protein. Supplementary Table 12 - Mutations and their effects on SETD2 protein. Supplementary Table 13 - Count of non-silent mutations per breed/per germline for OSA susceptibility genes Supplementary Table 14 - Overlap of GISTIC-defined recurrent SCNA genes between humans and dogs. See attached excel spreadsheet.

File contains genes that have focal somatic copy number amplifications

File contains additional information regarding methods for identifying putative driver genes with MutSigCV and high-impact mutations based analysis. Also methods for predicting impacts of TP53 and SETD2 genes. And further information about the oaCGH analysis for somatic copy number aberrations across a subset of 22 dogs.

PDF file - 843K, Unsupervised hierarchical clustering of cBCLs.

PDF file - 881K, Western blot comparing phospho-NF-κB (p65) levels in "ABC-like" and "GCB-like" canine samples.

Supplementary Tables 1-2 from Distinct B-Cell and T-Cell Lymphoproliferative Disease Prevalence among Dog Breeds Indicates Heritable Risk

XLSX file - 40K, Canine genes identified by Ingenuity Pathway Analysis as part of the NF-kB or B-cell receptor signaling pathways (Fig. 5).

PDF file - 60K

Immunophenotypes in lymphoproliferative diseases (LPD) are prognostically significant, yet causative factors for these conditions, and specifically those associated with heritable risk, remain elusive. The full spectrum of LPD seen in humans occurs in dogs, but the incidence and lifetime risk of naturally occurring LPD differs among dog breeds. Taking advantage of the limited genetic heterogeneity that exists within dog breeds, we tested the hypothesis that the prevalence of LPD immunophenotypes would differ among different breeds. The sample population included 1,263 dogs representing 87 breeds. Immunophenotype was determined by the presence of clonal rearrangements of immunoglobulin heavy chain or T-cell receptor γ chain. The probability of observing the number of B-cell or T-cell tumors in a particular breed or breed group was compared with three reference populations. Significance was computed using χ2 test, and logistic regression was used to confirm binomial predictions. The data show that, among 87 breeds tested, 15 showed significant differences from the prevalence of LPD immunophenotypes seen across the dog population as a whole. More significantly, elevated risk for T-cell LPD seems to have arisen ancestrally and is retained in related breed groups, whereas increased risk for B-cell disease may stem from different risk factors, or combinations of risk factors, arising during the process of breed derivation and selection. The data show that domestic dogs provide a unique and valuable resource to define factors that mediate risk as well as genes involved in the initiation of B-cell and T-cell LPD.

PDF file - 41K, Distribution of lymphoma samples into "Ongoing" or "Static" categories based on the proportion of IGHV subclones that are identical at the CDR3 region for each sample.

XLSX file - 116K, Canine genes that distinguish "ABC-like" and "GCB-like" subtypes (Fig. 5).

Background Patients presenting with head and neck squamous cell carcinoma of unknown primary (HNSCCUP) remain challenging clinical scenarios as large variation exists in practices used to locate the primary. Objective To perform a systematic review of the literature and offer recommendations for oropharyngeal biopsies in HNSCCUP. Method Pubmed, Medline and Embase were searched to identify studies from inception to October 2021. The Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) guidelines were followed. Results 483 articles were included and screened, 40 studies met the inclusion criteria, including over 3400 patients from the original articles and 1575 patients from 3 meta-analyses. The primary site identification rate following random biopsies or deep tissue biopsies is less than 5% in most studies. The mean detection rate following ipsilateral tonsillectomy is 34%; two pooled analyses indicate that the mean detection rate following tongue base mucosectomy is 64%, with this figure rising when the tonsils are negative. Conclusions High level evidence is lacking, with heterogeneity in the reported studies. Published meta analyses are based on retrospective data. There is little evidence supporting the practice of random/non-directed oropharyngeal biopsies. Available evidence supports palatine tonsillectomy and tongue base mucosectomy compared to deep tissue biopsies.

Background Natural and anthropogenic disasters can have long-lasting impacts on the genetics and structure of impacted populations. The 1986 Chernobyl Nuclear Power Plant disaster led to extensive contamination of the local environment and the wildlife therein. Several ecological, environmental, and genetic studies reported various effects of this disaster on animal, insect, and plant species; however, little work has been done to investigate the genetics of the free-breeding dogs that occupy the Chernobyl Exclusion Zone (CEZ). Results We define the population genetic structure of two groups of dogs that reside within the CEZ, one around the reactor site itself and another living within Chernobyl City. We found little evidence of gene flow and a significant degree of genetic differentiation between the two populations dogs, suggesting that these are two distinct populations despite occupying areas located just 16 km apart. With an FST-based outlier analysis, we then performed a genome-wide scan for evidence of directional selection within the dog populations. We found 391 outlier loci associated with genomic regions influenced by directional selection, from which we identified 52 candidate genes. Conclusions Our genome scan highlighted outlier loci within or near genomic regions under directional selection, possibly in response to the multi-generational exposure faced. In defining the population structure and identifying candidate genes for these dog populations, we take steps towards understanding how these types of prolonged exposures have impacted these populations.

Molecular profiling studies have shown that 85% of canine urothelial carcinomas (UC) harbor an activating BRAF V595E mutation, which is orthologous to the V600E variant found in several human cancer subtypes. In dogs, this mutation provides both a powerful diagnostic marker and a potential therapeutic target; however, due to their relative infrequency, the remaining 15% of cases remain understudied at the molecular level. We performed whole exome sequencing analysis of 28 canine urine sediments exhibiting the characteristic DNA copy number signatures of canine UC, in which the BRAF V595E mutation was undetected (UDV595E specimens). Among these we identified 13 specimens (46%) harboring short in-frame deletions within either BRAF exon 12 (7/28 cases) or MAP2K1 exons 2 or 3 (6/28 cases). Orthologous variants occur in several human cancer subtypes and confer structural changes to the protein product that are predictive of response to different classes of small molecule MAPK pathway inhibitors. DNA damage response and repair genes, and chromatin modifiers were also recurrently mutated in UDV595E specimens, as were genes that are positive predictors of immunotherapy response in human cancers. Our findings suggest that short in-frame deletions within BRAF exon 12 and MAP2K1 exons 2 and 3 in UDV595E cases are alternative MAPK-pathway activating events that may have significant therapeutic implications for selecting first-line treatment for canine UC. We developed a simple, cost-effective capillary electrophoresis genotyping assay for detection of these deletions in parallel with the BRAF V595E mutation. The identification of these deletion events in dogs offers a compelling cross-species platform in which to study the relationship between somatic alteration, protein conformation, and therapeutic sensitivity.

Sporadic angiosarcomas are aggressive vascular sarcomas whose rarity and genomic complexity present significant obstacles in deciphering the pathogenic significance of individual genetic alterations. Numerous fusion genes have been identified across multiple types of cancers, but their existence and significance remain unclear in sporadic angiosarcomas. In this study, we leveraged RNA-sequencing data from 13 human angiosarcomas and 76 spontaneous canine hemangiosarcomas to identify fusion genes associated with spontaneous vascular malignancies. Ten novel protein-coding fusion genes, including TEX2-PECAM1 and ATP8A2-FLT1, were identified in seven of the 13 human tumors, with two tumors showing mutations of TP53. HRAS and NRAS mutations were found in angiosarcomas without fusions or TP53 mutations. We found 15 novel protein-coding fusion genes including MYO16-PTK2, GABRA3-FLT1, and AKT3-XPNPEP1 in 11 of the 76 canine hemangiosarcomas; these fusion genes were seen exclusively in tumors of the angiogenic molecular subtype that contained recurrent mutations in TP53, PIK3CA, PIK3R1, and NRAS. In particular, fusion genes and mutations of TP53 cooccurred in tumors with higher frequency than expected by random chance, and they enriched gene signatures predicting activation of angiogenic pathways. Comparative transcriptomic analysis of human angiosarcomas and canine hemangiosarcomas identified shared molecular signatures associated with activation of PI3K/AKT/mTOR pathways. Our data suggest that genome instability induced by TP53 mutations might create a predisposition for fusion events that may contribute to tumor progression by promoting selection and/or enhancing fitness through activation of convergent angiogenic pathways in this vascular malignancy. Implications: This study shows that, while drive events of malignant vasoformative tumors of humans and dogs include diverse mutations and stochastic rearrangements that create novel fusion genes, convergent transcriptional programs govern the highly conserved morphologic organization and biological behavior of these tumors in both species.

To study neoplasia in tissue culture, cell lines representing the evolution of normal cells to tumor cells are needed. To produce such cells, we developed the AGMK1-9T7 cell line, established cell banks at 10-passage intervals, and characterized their biological properties. Here we examine the evolution of chromosomal DNA copy-number aberrations and miRNA expression in this cell line from passage 1 to the acquisition of a tumorigenic phenotype at passage 40. We demonstrated the use of a human microarray platform for DNA copy-number profiling of AGMK1-9T7 cells using knowledge of synteny to ‘recode’ data from human chromosome coordinates to those of the African green monkey. This approach revealed the accumulation of DNA copy-number gains and losses in AGMK1-9T7 cells from passage 3 to passage 40, which spans the period in which neoplastic transformation occurred. These alterations occurred in the sequences of genes regulating DNA copy-number imbalance of several genes that regulate endothelial cell angiogenesis, survival, migration, and proliferation. Regarding miRNA expression, 195 miRNAs were up- or down-regulated at passage 1 at levels that appear to be biologically relevant (i.e., log2 fold change >2.0 (q

Despite the strikingly high worldwide prevalence of vulvovaginal candidiasis (VVC), treatment options for recurrent VVC (RVVC) remain limited, with many women experiencing failed clinical treatment with frontline azoles. Further, the cause of onset and recurrence of disease is largely unknown, with few studies identifying potential mechanisms of treatment failure. This study aimed to assess a panel of clinical samples from healthy women and those with RVVC to investigate the influence of Candida, the vaginal microbiome, and how their interaction influences disease pathology. 16S rRNA sequencing characterized disease by a reduction in specific health-associated Lactobacillus species, such as Lactobacillus crispatus, coupled with an increase in Lactobacillus iners. In vitro analysis showed that Candida albicans clinical isolates are capable of heterogeneous biofilm formation, and we found the presence of hyphae and C. albicans aggregates in vaginal lavage fluid. Additionally, the ability of Lactobacillus to inhibit C. albicans biofilm formation and biofilm-related gene expression was demonstrated. Using RNA sequencing technology, we were able to identify a possible mechanism by which L. crispatus may contribute to re-establishing a healthy vaginal environment through amino acid acquisition from C. albicans. This study highlights the potential formation and impact of Candida biofilms in RVVC. Additionally, it suggests that RVVC is not entirely due to an arbitrary switch in C. albicans from commensal to pathogen and that understanding interactions between this yeast and vaginal Lactobacillus species may be crucial to elucidating the cause of RVVC and developing appropriate therapies. IMPORTANCE RVVC is a significant burden, both economically and for women's health, but its prevalence is poorly documented globally due to the levels of self-treatment. Identifying triggers for development and recurrence of VVC and the pathogenesis of the microbes involved could considerably improve prevention and treatment options for women with recurrent, azole-resistant cases. This study therefore aimed to examine the interkingdom dynamics from healthy women and those with RVVC using next-generation sequencing techniques and to further investigate the molecular interactions between C. albicans and L. crispatus in a relevant biofilm coculture system.

An amendment to this paper has been published and can be accessed via the original article.

Introduction Iatrogenic facial nerve palsy is distressing to the patient and clinician. The deformity is aesthetically displeasing, and can be functionality problematic for oral competence, dental lip trauma and speech. Furthermore such injuries have litigation implications. Marginal mandibular nerve (MMN) palsy causes an obvious asymmetrical smile. MMN is at particular risk during procedures such as rhytidoplasties, mandibular fracture, tumour resection and neck dissections. Cited causes for the high incidence are large anatomical variations, unreliable landmarks, an exposed course and tumour grade or nodal involvement dictating requisite nerve sacrifice. An alternative cause for post-operative asymmetry is damage to the cervical branch of the facial nerve or platysmal dysfunction. This tends to have a transient course and recovers. Distinction between MMN palsy and palsy of the cervical branch of the facial nerve should therefore be made. In 1979 Ellenbogen differentiated between MMN palsy and “Pseudo-paralysis of the mandibular branch of the facial nerve”. Despite this, there is paucity in the literature & confusion amongst clinicians in distinguishing between these palsies, and there is little regarding these post-operative sequelae and neck dissections. Method This article reflects on the surgical anatomy of the MMN and cervical nerve in relation to danger zones during lymphadenectomy. The authors review the anatomy of the smile. Finally, we utilise case studies to evaluate the differences between MMN palsy and its pseudo-palsy to allow clinical differentiation. Conclusion Here we present a simple method for clinical differentiation between these two prognostically different injuries, allowing appropriate reassurance, therapy & management.

Vulvovaginal Candidiasis (VVC) is the most prevalent Candida infection in humans affecting 75% of women at least once throughout their lifetime. In its debilitating recurrent form, RVVC is estimated to affect 140 million women annually. Despite this strikingly high prevalence, treatment options for RVVC remain limited with many women experiencing failed clinical treatment with frontline azoles. Further, the cause of onset and recurrence of disease is largely unknown with few studies identifying potential mechanisms of failed treatment. This study aimed to assess a panel of clinical samples from healthy women and those with RVVC to investigate the influence of Candida, vaginal microbiome and antagonism between Candida and Lactobacillus on disease pathology. 16S rRNA sequencing characterised disease by a reduction in specific health-associated Lactobacillus such as L. crispatus, coupled with an increase in L. iners. In vitro analysis showed Candida albicans clinical isolates are capable of heterogeneous biofilm formation and show the presence of hyphae and C. albicans aggregates in vaginal lavage. Additionally, the ability of Lactobacillus to inhibit C. albicans biofilm formation and biofilm-related gene expression was demonstrated. Using RNA sequencing technology, we were able to exploit a possible mechanism by which L. crispatus may aim to re-establish a healthy vaginal environment through amino-acid acquisition from C. albicans. This study suggests RVVC is not entirely due to an arbitrary switch in C. albicans from commensal to pathogen and understanding interactions between the yeast and vaginal Lactobacillus species may be more crucial to elucidating the cause of RVVC and developing appropriate therapies.

Mast cell tumor (MCT) is the most common skin malignancy of domestic dogs and presents with a widely variable clinical behavior. Although activating KIT mutations are present in approximately 20% of canine MCTs, molecular etiology is largely unknown for the majority of this cancer. Characterization of genomic alterations in canine MCTs may identify genomic regions and/or genes responsible for their development and progression, facilitating the discovery of new therapeutic targets and improved clinical management of this heterogeneous cancer. We performed genome-wide DNA copy number analysis of 109 primary MCTs derived from three popular canine breeds (the Boxer, Labrador Retriever, and Pug) as well as nontarget breeds using oligonucleotide array comparative genomic hybridization (oaCGH). We demonstrated a stepwise accumulation of numerical DNA copy number aberrations (CNAs) as tumor grade increases. DNA sequencing analysis revealed that KIT mutations were found less frequently in the Pug tumors and were strongly associated with high histological grade. Tumors with KIT mutations showed genome-wide aberrant copy number profiles, with frequent CNAs involving genes in the p53 and RB pathways, whereas CNAs were very limited in tumors with wild-type KIT. We evaluated the presence of four CNAs to predict aggressive tumor phenotypes. This approach predicted aggressive tumors with a sensitivity of 78-94% and specificity of 88-93%, when using oaCGH and droplet digital PCR platforms. Further investigation of genome regions identified in this study may lead to the development of a molecular tool for classification and prognosis, as well as identification of therapeutic target molecules.

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 histiocytic malignancies (HM) are rare across the general dog population, but overrepresented in certain breeds, such as Bernese mountain dog and flat-coated retriever. Accurate diagnosis relies on immunohistochemical staining to rule out histologically similar cancers with different prognoses and treatment strategies (e.g., lymphoma and hemangiosarcoma). HM are generally treatment refractory with overall survival of less than 6 months. A lack of understanding regarding the mechanisms of disease development and progression hinders development of novel therapeutics. While the study of human tumors can benefit veterinary medicine, the rarity of the suggested orthologous disease (dendritic cell sarcoma) precludes this. This study aims to improve the understanding of underlying disease mechanisms using genome-wide DNA copy number and gene expression analysis of spontaneous HM across several dog breeds. Extensive DNA copy number disruption was evident, with losses of segments of chromosomes 16 and 31 detected in 93% and 72% of tumors, respectively. Droplet digital PCR (ddPCR) evaluation of these regions in numerous cancer specimens effectively discriminated HM from other common round cell tumors, including lymphoma and hemangiosarcoma, resulting in a novel, rapid diagnostic aid for veterinary medicine. Transcriptional analysis demonstrated disruption of the spindle assembly complex, which is linked to genomic instability and reduced therapeutic impact in humans. A key signature detected was up-regulation of Matrix Metalloproteinase 9 (MMP9), supported by an immunohistochemistry-based assessment of MMP9 protein levels. Since MMP9 has been linked with rapid metastasis and tumor aggression in humans, the data in this study offer a possible mechanism of aggression in HM.

Follicular lymphomas (FLs) typically exhibit a chromosome translocation that induces constitutive expression of the anti-apoptotic bcl2 protein and accumulation of additional molecular defects. This rearrangement offers a promising therapeutic target, but its nature as a fundamental driver of FL pathogenesis remains unclear as 15% of cases lack the translocation. We performed an integrated immunohistochemical and genomic investigation of 10 naturally occurring FL cases from domestic dogs, showing that, as with human tumours, they exhibit marked heterogeneity in the frequency and intensity of bcl2 protein expression. Genomic copy number aberrations were infrequent and broadly consistent with those of other canine B-cell lymphoma subtypes. None of the canine FL specimens exhibited a rearrangement consistent with the hallmark translocation of human FL, despite their remarkable histomorphologic similarity. Parallel exploration of canine and human cases may reveal alternative tumour-initiating mechanisms other than BCL2 disruption, yielding a more complete definition of the molecular pathogenesis of FL.

Angiosarcoma is an aggressive, albeit rare cancer in humans. Angiosarcomas are vascular malignancies that can occur anywhere in the body, and their metastatic propensity is high. The cause of the vast majority of sporadic angiosarcomas is unknown, and no therapeutic targets have been identified to improve outcomes. Angiosarcomas are genomically complex, with chaotic karyotypes and massive chromosomal abnormalities. Despite the genomic complexity, angiosarcomas share a histological morphology that consists of disorganized, malignant vessel-forming cells. Hemangiosarcoma (HSA) is a common cancer of dogs; it shares clinical and morphological features, as well as aspects of its mutational landscape with human angiosarcoma. Our previous work has revealed that canine HSAs and human angiosarcomas share transcriptional signatures that establish angiogenic and inflammatory molecular subtypes. A comparative genomics approach is useful to apply knowledge from appropriately powered canine studies to inform research into human angiosarcomas. In this study, we leveraged next generation RNA-Seq data from a cohort of 76 spontaneous canine HSAs and from thirteen human angiosarcomas to identify fusion genes. Fifteen novel protein-coding fusion genes including MYO16-PTK2, GABRA3-FLT1, and AKT3-XPNPEP1 were identified in 11 of the 76 canine HSAs; these fusion genes were exclusively seen in tumors of the angiogenic molecular subtype. Mutations of TP53 and fusion genes co-occurred in tumors with higher frequency than expected by random chance. Pathway analysis revealed that co-occurring mutations of TP53 and PIK3CA were associated with gene expression signatures of chromatin remodeling and immunosuppression, and co-occurrence of fusion genes and TP53 mutation enriched a gene signature predicting activation of angiogenic pathways. In human angiosarcomas, ten novel protein-coding fusion genes, including TEX2-PECAM1 and ATP8A2-FLT1, were identified in 7 of 13 tumors, with two showing mutations of TP53. Immunohistochemical assays showed that canine HSA and human angiosarcoma activated p53, AKT, and mTOR signaling pathways independent of fusion genes and mutational conditions, suggesting that both tumors activate convergent signaling pathways to retain the ontogenetical properties. In summary, our comparative analysis identified shared molecular signatures between canine HSA and human angiosarcoma. Specifically, the data suggests that genomic instability induced by TP53 mutations might create a predisposition for fusion events that may contribute to tumor progression by promoting selection and/or enhancing fitness through activation of convergent angiogenic pathways. Our ongoing work seeks to define the key molecular programs that establish the mutational landscape, which consequently activates convergent pathways that contribute to angiosarcoma development. Citation Format: Jong Hyuk Kim, Kate Megquier, Aaron L. Sarver, Rachael Thomas, Ashley J. Schulte, Chao Wang, Ingegerd Elvers, Elinor Karlsson, Matthew Breen, Kerstin Lindblad-Toh, Jaime F. Modiano. Molecular mechanisms that activate convergent oncogenic pathway in genomically complex angiosarcoma [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 195.

The goals of this work are biomarker and drug discovery to advance cancer diagnostic and therapeutic strategies in humans through the study of naturally-occurring canine cancers. Many factors, including shared environment, intact host immunity, and greater cancer gene family homology between dogs and humans than mice and humans, make spontaneous canine cancers valuable complementary models of human cancer. Transcriptomic profiling utilizing RNA sequencing (RNA SEQ) of 5 canine cancers, including osteosarcoma, melanoma, B cell lymphoma, T cell lymphoma, and pulmonary carcinoma, have revealed five distinct gene co-expression models. From these unique module expression profiles, cancer-specific gene panels were derived. A similar analysis performed on existing RNA-SEQ data from human tumor samples produced cancer-specific human gene panels. Comparison of the canine and human gene panels found significant correlation (Spearman correlation > 0.6) which supports the translational relevance of naturally-occurring canine cancers to their human counterparts. Further, proteomic profiles derived from matched tumor tissue and peripheral blood samples mirror those of the tumor transcriptome, demonstrating these cancer-specific gene panels and their encoded proteins may represent robust canine diagnostic biomarker and/or therapeutic target candidates. Therapeutic hypotheses associated with each cancer specific gene panel were derived through matching of drugs documented to alter expression of panel genes in opposition to that exhibited by each cancer type. We identified 60 candidate drugs and screened them against a panel of canine cancer cell lines, finding 40 drugs that inhibited cell growth by 75% or more. Three or more active compounds were found for each cell line. From these 40 active compounds we then derived 30 synergistic drug combinations with the requirement that that they alter two or more panel genes in opposition to that exhibited by each cancer type. Additional studies to document drug synergism are underway and those confirmed as such will be considered good drug combination candidates for future canine clinical trials. Biomarker and drug combination candidates that perform well in canine clinical trials will then be considered for human trials. This work exemplifies the type of approach meant to further establish the comparative relevance of canine to human cancer and provide opportunities to explore hypotheses related to detection and treatment in both species. Citation Format: Amy Kathleen LeBlanc, Christina N. Mazcko, Matthew Breen, Rachael Thomas, Douglas H. Thamm. A comparative oncology approach to biomarker and drug discovery for cancer diagnosis and treatment in dogs and humans [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 6134.

The utility of the domestic cat as a model system for biomedical studies was constrained for many years by the absence of a comprehensive feline reference genome sequence assembly. While such a resource now exists, the cat continues to lag behind the domestic dog in terms of integration into the &lsquo, One Health&rsquo, era of molecular medicine. Stimulated by the advances being made within the evolving field of comparative cancer genomics, we developed a microarray platform that allows rapid and sensitive detection of DNA copy number aberrations in feline tumors using comparative genomic hybridization analysis. The microarray comprises 110,456 unique oligonucleotide probes anchored at mean intervals of 22.6 kb throughout the feline reference genome sequence assembly, providing ~350-fold higher resolution than was previously possible using this technique. We demonstrate the utility of this resource through genomic profiling of a feline injection-site sarcoma case, revealing a highly disrupted profile of DNA copy number imbalance involving several key cancer-associated genes including KIT, TP53, PTEN, FAS and RB1. These findings were supported by targeted fluorescence in-situ hybridization analysis, which identified major alterations in chromosome structure, including complex intrachromosomal reorganization events typical of those seen in aggressive soft-tissue sarcomas of other species. We then characterized a second mass that was identified at a nearby site in the same patient almost 12 months later. This mass demonstrated a remarkably conserved genomic profile consistent with a recurrence of the original tumor, however the detection of subtle differences reflected evolution of the tumor over time. These findings exemplify the diverse potential of this microarray platform to incorporate domestic cat cancers into comparative and translational research efforts in molecular oncology.

OBJECTIVE To evaluate gene expression and DNA copy number in adipose tissue-derived stromal cells (ADSCs) and in ADSC-derived neurosphere-like cell clusters (ADSC-NSCs) generated from tissues of chronically paraplegic dogs. ANIMALS 14 client-owned paraplegic dogs. PROCEDURES Dorsal subcutaneous adipose tissue (< 1 cm3) was collected under general anesthesia; ADSCs were isolated and cultured. Third-passage ADSCs were cultured in neural cell induction medium to generate ADSC-NSCs. Relative gene expression of mesenchymal cell surface marker CD90 and neural progenitor marker nestin was assessed in ADSCs and ADSC-NSCs from 3 dogs by quantitative real-time PCR assay; expression of these and various neural lineage genes was evaluated for the same dogs by reverse transcription PCR assay. Percentages of cells expressing CD90, nestin, glial fibrillary acidic protein (GFAP), and tubulin β 3 class III (TUJ1) proteins were determined by flow cytometry for all dogs. The DNA copy number stability (in samples from 6 dogs) and neural cell differentiation (14 dogs) were assessed with array-comparative genomic hybridization analysis and immunocytochemical evaluation, respectively. RESULTS ADSCs and ADSC-NSCs expressed neural cell progenitor and differentiation markers; GFAP and microtubule-associated protein 2 were expressed by ADSC-NSCs but not ADSCs. Relative gene expression of CD90 and nestin was subjectively higher in ADSC-NSCs than in ADSCs. Percentages of ADSC-NSCs expressing nestin, GFAP, and TUJ1 proteins were substantially higher than those of ADSCs. Cells expressing neuronal and glial markers were generated from ADSC-NSCs and had no DNA copy number instability detectable by the methods used. CONCLUSIONS AND CLINICAL RELEVANCE Results suggested ADSCs can potentially be a safe and clinically relevant autologous source for canine neural progenitor cells. Further research is needed to verify these findings.

BACKGROUND In 2012, over 240,000 men were diagnosed with prostate cancer and over 28,000 died from the disease. Animal models of prostate cancer are vital to understanding its pathogenesis and developing therapeutics. Canine models in particular are useful due to their similarities to late-stage, castration-resistant human disease with osteoblastic bone metastases. This study established and characterized a novel canine prostate cancer cell line that will contribute to the understanding of prostate cancer pathogenesis.

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype, and fewer than half of patients are cured with standard first-line therapy. To improve therapeutic options, better animal models that accurately mimic human DLBCL (hDLBCL) are needed. Canine DLBCL, one of the most common cancers in veterinary oncology, is morphologically similar to hDLBCL and is treated using similar chemotherapeutic protocols. With genomic technologies, it is now possible to molecularly evaluate dogs as a potential large-animal model for hDLBCL. We evaluated canine B-cell lymphomas (cBCL) using immunohistochemistry (IHC) and gene expression profiling. cBCL expression profiles were similar in many ways to hDLBCLs. For instance, a subset had increased expression of NF-κB pathway genes, mirroring human activated B-cell (ABC)–type DLBCL. Furthermore, immunoglobulin heavy chain ongoing mutation status, which is correlated with ABC/germinal center B-cell cell of origin in hDLBCL, separated cBCL into two groups with statistically different progression-free and overall survival times. In contrast with hDLBCL, cBCL rarely expressed BCL6 and MUM1/IRF4 by IHC. Collectively, these studies identify molecular similarities to hDLBCL that introduce pet dogs as a representative model of hDLBCL for future studies, including therapeutic clinical trials. Cancer Res; 73(16); 5029–39. ©2013 AACR.

Established cell lines are utilized extensively to study tumor biology and preclinical therapeutic development. However, they may not accurately recapitulate the heterogeneity of their corresponding primary disease. B-cell tumor cells are especially difficult to maintain under conventional culture conditions, limiting access to samples that faithfully represent this disease for preclinical studies. Here, we used primary canine diffuse large B-cell lymphoma to establish a culture system that reliably supports the growth of these cells. CD40 ligand, either expressed by feeder cells or provided as a soluble two-trimeric form, was sufficient to support primary lymphoma cells in vitro. The tumor cells retained their original phenotype, clonality and known karyotypic abnormalities after extended expansion in culture. Finally, we illustrate the utility of the feeder cell-free culture system for comparable assessment of cytotoxicity using dog and human B-cell malignancies. We conclude that this system has broad applications for in vitro preclinical development for B-cell malignancies.

Molecular characterization of tumour cell lines is increasingly regarded as a prerequisite for defining their validity as models of in vivo neoplasia. We present the first comprehensive catalogue of genomic and transcriptional characteristics of five widely used canine lymphoid tumour cell lines. High-resolution microarray-based comparative genomic hybridization defined their unique profiles of genomic DNA copy number imbalance. Multicolour fluorescence in situ hybridization identified aberrant gains of MYC, KIT and FLT3 and deletions of PTEN and CDKN2 in individual cell lines, and also revealed examples of extensive structural chromosome reorganization. Gene expression profiling and RT-PCR analyses defined the relationship between genomic imbalance and transcriptional dysregulation in each cell line, clarifying their relevance as models of discrete functional pathways with biological and therapeutic significance. In combination, these data provide an extensive resource of molecular data for directing the appropriate use of these cell lines as tools for studying canine lymphoid neoplasia.

Osteosarcoma (OS) is the most commonly diagnosed malignant bone tumor in humans and dogs, characterized in both species by extremely complex karyotypes exhibiting high frequencies of genomic imbalance. Evaluation of genomic signatures in human OS using array comparative genomic hybridization (aCGH) has assisted in uncovering genetic mechanisms that result in disease phenotype. Previous low-resolution (10‐20 Mb) aCGH analysis of canine OS identified a wide range of recurrent DNA copy number aberrations, indicating extensive genomic instability. In this study, we profiled 123 canine OS tumors by 1 Mb-resolution aCGH to generate a dataset for direct comparison with current data for human OS, concluding that several high frequency aberrations in canine and human OS are orthologous. To ensure complete coverage of gene annotation, we identified the human refseq genes that map to these orthologous aberrant dog regions and found several candidate genes warranting evaluation for OS involvement. Specifically, subsequenct FISH and qRT-PCR analysis of RUNX2, TUSC3, and PTEN indicated that expression levels correlated with genomic copy number status, showcasing RUNX2 as an OS associated gene and TUSC3 as a possible tumor suppressor candidate. Together these data demonstrate the ability of genomic comparative oncology to identify genetic abberations which may be important for OS progression. Large scale screening of genomic imbalance in canine OS further validates the use of the dog as a suitable model for human cancers, supporting the idea that dysregulation discovered in canine cancers will provide an avenue for complementary study in human counterparts. V C 2011 Wiley-Liss, Inc.

Identification of the genomic regions most intimately associated with non-Hodgkin's lymphoma (NHL) pathogenesis is confounded by the genetic heterogeneity of human populations. We hypothesize that the restricted genetic variation of purebred dogs, combined with the contrasting architecture of the human and canine karyotypes, will increase the penetrance of fundamental NHL-associated chromosomal aberrations in both species. We surveyed non-random aneuploidy in 150 canine NHL cases, revealing limited genomic instability compared to their human counterparts and no evidence for CDKN2A/B deletion in canine B-cell NHL. ‘Genomic recoding’ of canine NHL data into a ‘virtual human’ chromosome format showed remarkably few regions of copy number aberration (CNA) shared between both species; restricted to regions of dog chromosomes 13 and 31, and human chromosomes 8 and 21. Our data suggest that gene discovery in NHL may be enhanced through comparative studies exploiting the less complex association between CNAs and tumor pathogenesis in canine patients.

Angiosarcoma is an aggressive, albeit rare cancer in humans. The cause of the vast majority of sporadic angiosarcomas is unknown, mortality is high, and no therapeutic targets have been identified to improve outcomes. Hemangiosarcoma (HSA) is a common cancer of dogs, and it shares histopathologic features with human angiosarcoma. In our previous work, canine HSAs were classified into angiogenic, inflammatory, and adipogenic subtypes based on transcriptional profiles. However, the genetic and molecular events that regulate transcriptional subtypes in angiosarcoma are not currently understood. Our goal was to use a comparative genomics approach to apply knowledge from appropriately powered canine studies to inform our research into human sarcomas. In this study, we identified recurrent mutations in RNASeq data from 93 HSAs and 16 nonmalignant controls, based on mutations first identified in exomes from 42 paired tumor and normal samples. In addition to identifying recurrent somatic mutations we also identified translocation fusions, allowing elucidation of oncogenic mechanisms for vascular endothelial growth factor receptors (VEGFR), phosphoinositide-3 kinase (PIK3) signaling pathways, and the p53 DNA damage repair pathway in canine HSA. Significantly, mutational signatures were associated with distinct molecular subtypes of canine hemangiosarcomas, and both the angiogenic and the inflammatory subtypes were apparent in RNASeq data from human angiosarcomas (n=14), suggesting that comparable etiologic mechanisms are operative in the canine and human disease. Our ongoing work seeks to understand how the molecular mechanisms give rise to molecular subtypes of angiosarcoma by defining the association between driver mutations, signaling pathway alterations and transcriptional patterns, which should allow us to identify rational therapeutic targets. Citation Format: Jong Hyuk Kim, Kate Megquier, Aaron L. Sarver, Rachael Thomas, Chao Wang, Ingegerd Elvers, Elinor Karlsson, Matthew Breen, Kerstin Lindblad-Toh, Jaime F. Modiano. Mutational and transcriptomic profiling identify distinct angiogenic and inflammatory subtypes of angiosarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5357.

Despite their clinical significance and substantial human health burden, fungal infections remain relatively under-appreciated. The widespread overuse of antibiotics and the increasing requirement for indwelling medical devices provides an opportunistic potential for the overgrowth and colonization of pathogenic Candida species on both biological and inert substrates. Indeed, it is now widely recognized that biofilms are a highly important part of their virulence repertoire. Candida albicans is regarded as the primary fungal biofilm forming species, yet there is also increasing interest and growing body of evidence for non-Candida albicans species (NCAS) biofilms, and interkingdom biofilm interactions. C. albicans biofilms are heterogeneous structures by definition, existing as three-dimensional populations of yeast, pseudo-hyphae, and hyphae, embedded within a self-produced extracellular matrix. Classical molecular approaches, driven by extensive studies of laboratory strains and mutants, have enhanced our knowledge and understanding of how these complex communities develop, thrive, and cause host-mediated damage. Yet our clinical observations tell a different story, with differential patient responses potentially due to inherent biological heterogeneity from specific clinical isolates associated with their infections. This review explores some of the recent advances made in an attempt to explore the importance of working with clinical isolates, and what this has taught us.