Canine Oral Melanoma Genomic and Transcriptomic Study Defines Two Molecular Subgroups with Different Therapeutical Targets.

Mucosal melanoma (MM) is a rare, aggressive clinical cancer. Despite recent advances in genetics and treatment, the prognosis of MM remains poor. Canine MM offers a relevant spontaneous and immunocompetent model to decipher the genetic bases and explore treatments for MM. We performed an integrative genomic and transcriptomic analysis of 32 canine MM samples, which identified two molecular subgroups with a different microenvironment and structural variant (SV) content. The overexpression of genes related to the microenvironment and T-cell response was associated with tumors harboring a lower content of SVs, whereas the overexpression of pigmentation-related pathways and oncogenes, such as TERT, was associated with a high SV burden. Using whole-genome sequencing, we showed that focal amplifications characterized complex chromosomal rearrangements targeting oncogenes, such as MDM2 or CDK4, and a recurrently amplified region on canine chromosome 30. We also demonstrated that the genes TRPM7, GABPB1, and SPPL2A, located in this CFA30 region, play a role in cell proliferation, and thus, may be considered as new candidate oncogenes for human MM. Our findings suggest the existence of two MM molecular subgroups that may benefit from dedicated therapies, such as immune checkpoint inhibitors or targeted therapies, for both human and veterinary medicine.

Identification of common predisposing loci to hematopoietic cancers in four dog breeds.

Histiocytic sarcoma (HS) is a rare but aggressive cancer in both humans and dogs. The spontaneous canine model, which has clinical, epidemiological, and histological similarities with human HS and specific breed predispositions, provides a unique opportunity to unravel the genetic basis of this cancer. In this study, we aimed to identify germline risk factors associated with the development of HS in canine-predisposed breeds. We used a methodology that combined several genome-wide association studies in a multi-breed and multi-cancer approach as well as targeted next-generation sequencing, and imputation We combined several dog breeds (Bernese mountain dogs, Rottweilers, flat-coated retrievers, and golden retrievers), and three hematopoietic cancers (HS, lymphoma, and mast cell tumor). Results showed that we not only refined the previously identified HS risk CDKN2A locus, but also identified new loci on canine chromosomes 2, 5, 14, and 20. Capture and targeted sequencing of specific loci suggested the existence of regulatory variants in non-coding regions and methylation mechanisms linked to risk haplotypes, which lead to strong cancer predisposition in specific dog breeds. We also showed that these canine cancer predisposing loci appeared to be due to the additive effect of several risk haplotypes involved in other hematopoietic cancers such as lymphoma or mast cell tumors as well. This illustrates the pleiotropic nature of these canine cancer loci as observed in human oncology, thereby reinforcing the interest of predisposed dog breeds to study cancer initiation and progression.

Circulating tumor DNA is detectable in canine histiocytic sarcoma, oral malignant melanoma, and multicentric lymphoma.

Circulating tumor DNA (ctDNA) has become an attractive biomarker in human oncology, and its use may be informative in canine cancer. Thus, we used droplet digital PCR or PCR for antigen receptor rearrangement, to explore tumor-specific point mutations, copy number alterations, and chromosomal rearrangements in the plasma of cancer-affected dogs. We detected ctDNA in 21/23 (91.3%) of histiocytic sarcoma (HS), 2/8 (25%) of oral melanoma, and 12/13 (92.3%) of lymphoma cases. The utility of ctDNA in diagnosing HS was explored in 133 dogs, including 49 with HS, and the screening of recurrent PTPN11 mutations in plasma had a specificity of 98.8% and a sensitivity between 42.8 and 77% according to the clinical presentation of HS. Sensitivity was greater in visceral forms and especially related to pulmonary location. Follow-up of four dogs by targeting lymphoma-specific antigen receptor rearrangement in plasma showed that minimal residual disease detection was concordant with clinical evaluation and treatment response. Thus, our study shows that ctDNA is detectable in the plasma of cancer-affected dogs and is a promising biomarker for diagnosis and clinical follow-up. ctDNA detection appears to be useful in comparative oncology research due to growing interest in the study of natural canine tumors and exploration of new therapies.

Cross-species analysis of enhancer logic using deep learning.

Deciphering the genomic regulatory code of enhancers is a key challenge in biology because this code underlies cellular identity. A better understanding of how enhancers work will improve the interpretation of noncoding genome variation and empower the generation of cell type-specific drivers for gene therapy. Here, we explore the combination of deep learning and cross-species chromatin accessibility profiling to build explainable enhancer models. We apply this strategy to decipher the enhancer code in melanoma, a relevant case study owing to the presence of distinct melanoma cell states. We trained and validated a deep learning model, called DeepMEL, using chromatin accessibility data of 26 melanoma samples across six different species. We show the accuracy of DeepMEL predictions on the CAGI5 challenge, where it significantly outperforms existing models on the melanoma enhancer of IRF4 Next, we exploit DeepMEL to analyze enhancer architectures and identify accurate transcription factor binding sites for the core regulatory complexes in the two different melanoma states, with distinct roles for each transcription factor, in terms of nucleosome displacement or enhancer activation. Finally, DeepMEL identifies orthologous enhancers across distantly related species, where sequence alignment fails, and the model highlights specific nucleotide substitutions that underlie enhancer turnover. DeepMEL can be used from the Kipoi database to predict and optimize candidate enhancers and to prioritize enhancer mutations. In addition, our computational strategy can be applied to other cancer or normal cell types.

PTPN11 mutations in canine and human disseminated histiocytic sarcoma.

In humans, histiocytic sarcoma (HS) is an aggressive cancer involving histiocytes. Its rarity and heterogeneity explain that treatment remains a challenge. Sharing high clinical and histopathological similarities with human HS, the canine HS is conversely frequent in specific breeds and thus constitutes a unique spontaneous model for human HS to decipher the genetic bases and to explore therapeutic options. We identified sequence alterations in the MAPK pathway in at least 63.9% (71/111) of HS cases with mutually exclusive BRAF (0.9%; 1/111), KRAS (7.2%; 8/111) and PTPN11 (56.75%; 63/111) mutations concentrated at hotspots common to human cancers. Recurrent PTPN11 mutations are associated to visceral disseminated HS subtype in dogs, the most aggressive clinical presentation. We then identified PTPN11 mutations in 3/19 (15.7%) human HS patients. Thus, we propose PTPN11 mutations as key events for a specific subset of human and canine HS: the visceral disseminated form. Finally, by testing drugs targeting the MAPK pathway in eight canine HS cell lines, we identified a better anti-proliferation activity of MEK inhibitors than PTPN11 inhibitors in canine HS neoplastic cells. In combination, these results illustrate the relevance of naturally affected dogs in deciphering genetic mechanisms and selecting efficient targeted therapies for such rare and aggressive cancers in humans.

Prognostic value of somatic focal amplifications on chromosome 30 in canine oral melanoma.

Canine oral melanoma is the first malignancy of the oral cavity in dogs and is characterized by a local invasiveness and a high metastatic propensity. A better knowledge of genetic alterations is expected to improve management of this tumour. Copy number alterations are known characteristics of mucosal melanomas both in dogs and humans. The goal of this study was to explore the prognostic value of somatic focal amplifications on chromosomes (Canis Familiaris [CFA]) 10 and 30 in canine oral melanoma. The cohort included 73 dogs with oral melanoma confirmed by histology, removed surgically without adjuvant therapy and with a minimal follow-up of 6 months. Epidemiological, clinical and histological data were collected and quantitative-PCR were performed on formalin-fixed paraffin-embedded (FFPE) samples to identify specific focal amplifications. The 73 dogs included in the study had a median survival time of 220 days. Focal amplifications on CFA 10 and 30 were recurrent (49.3% and 50.7% of cases, respectively) and CFA 30 amplification was significantly associated with the amelanotic phenotype (P = .046) and high mitotic index (MI; P = .0039). CFA 30 amplification was also linked to poor prognosis (P = .0005). Other negative prognostic factors included gingiva location (P = .003), lymphadenomegaly (P = .026), tumour ulceration at diagnosis (P = .003), MI superior to 6 mitoses over 10 fields (P = .001) and amelanotic tumour (P = .029). In multivariate analyses using Cox proportional hazards regression, CFA 30 amplification (Hazard ratio [HR] = 2.08; P = .011), tumour location (HR = 2.20; P = .005) and histological pigmentation (HR = 1.87; P = .036) were significantly associated with shorter survival time. Focal amplification of CFA 30 is linked to an aggressive subset and constitutes a new prognostic factor.

Genome-Wide Analysis of Long Non-Coding RNA Profiles in Canine Oral Melanomas.

Mucosal melanomas (MM) are rare aggressive cancers in humans, and one of the most common forms of oral cancers in dogs. Similar biological and histological features are shared between MM in both species, making dogs a powerful model for comparative oncology studies of melanomas. Although exome sequencing recently identified recurrent coding mutations in canine MM, little is known about changes in non-coding gene expression, and more particularly, in canine long non-coding RNAs (lncRNAs), which are commonly dysregulated in human cancers. Here, we sampled a large cohort (n = 52) of canine normal/tumor oral MM from three predisposed breeds (poodles, Labrador retrievers, and golden retrievers), and used deep transcriptome sequencing to identify more than 400 differentially expressed (DE) lncRNAs. We further prioritized candidate lncRNAs by comparative genomic analysis to pinpoint 26 dog-human conserved DE lncRNAs, including SOX21-AS, ZEB2-AS, and CASC15 lncRNAs. Using unsupervised co-expression network analysis with coding genes, we inferred the potential functions of the DE lncRNAs, suggesting associations with cancer-related genes, cell cycle, and carbohydrate metabolism Gene Ontology (GO) terms. Finally, we exploited our multi-breed design to identify DE lncRNAs within breeds. This study provides a unique transcriptomic resource for studying oral melanoma in dogs, and highlights lncRNAs that may potentially be diagnostic or therapeutic targets for human and veterinary medicine.

Discovery of Human-Similar Gene Fusions in Canine Cancers.

Canine cancers represent a tremendous natural resource due to their incidence and striking similarities to human cancers, sharing similar clinical and pathologic features as well as oncogenic events, including identical somatic mutations. Considering the importance of gene fusions as driver alterations, we explored their relevance in canine cancers. We focused on three distinct human-comparable canine cancers representing different tissues and embryonic origins. Through RNA-Seq, we discovered similar gene fusions as those found in their human counterparts: IGK-CCND3 in B-cell lymphoma, MPB-BRAF in glioma, and COL3A1-PDGFB in dermatofibrosarcoma protuberans-like. We showed not only similar partner genes but also identical breakpoints leading to oncogene overexpression. This study demonstrates similar gene fusion partners and mechanisms in human-dog corresponding tumors and allows for selection of targeted therapies in preclinical and clinical trials with pet dogs prior to human trials, within the framework of personalized medicine. Cancer Res; 77(21); 5721-7. ©2017 AACR.

Naturally occurring melanomas in dogs as models for non-UV pathways of human melanomas.

Spontaneously occurring melanomas are frequent in dogs. They appear at the same localizations as in humans, i.e. skin, mucosal sites, nail matrix and eyes. They display variable behaviors: tumors at oral localizations are more frequent and aggressive than at other anatomical sites. Interestingly, dog melanomas are associated with strong breed predispositions and overrepresentation of black-coated dogs. Epidemiological analysis of 2350 affected dogs showed that poodles are at high risk of developing oral melanoma, while schnauzers or Beauce shepherds mostly developped cutaneous melanoma. Clinical and histopathological analyses were performed on a cohort of 153 cases with a 4-yr follow-up. Histopathological characterization showed that most canine tumors are intradermal and homologous to human rare morphological melanomas types - 'nevocytoid type' and 'animal type'-. Tumor cDNA sequencing data, obtained from 95 dogs for six genes, relevant to human melanoma classification, detected somatic mutations in oral melanoma, in NRAS and PTEN genes, at human hotspot sites, but not in BRAF. Altogether, these findings support the relevance of the dog model for comparative oncology of melanomas, especially for the elucidation of non-UV induced pathways.

The MTAP-CDKN2A locus confers susceptibility to a naturally occurring canine cancer.

BACKGROUND: Advantages offered by canine population substructure, combined with clinical presentations similar to human disorders, makes the dog an attractive system for studies of cancer genetics. Cancers that have been difficult to study in human families or populations are of particular interest. Histiocytic sarcoma is a rare and poorly understood neoplasm in humans that occurs in 15% to 25% of Bernese Mountain Dogs (BMD). METHODS: Genomic DNA was collected from affected and unaffected BMD in North America and Europe. Both independent and combined genome-wide association studies (GWAS) were used to identify cancer-associated loci. Fine mapping and sequencing narrowed the primary locus to a single gene region. RESULTS: Both populations shared the same primary locus, which features a single haplotype spanning MTAP and part of CDKN2A and is present in 96% of affected BMD. The haplotype is within the region homologous to human chromosome 9p21, which has been implicated in several types of cancer. CONCLUSIONS: We present the first GWAS for histiocytic sarcoma in any species. The data identify an associated haplotype in the highly cited tumor suppressor locus near CDKN2A. These data show the power of studying distinctive malignancies in highly predisposed dog breeds. IMPACT: Here, we establish a naturally occurring model of cancer susceptibility due to CDKN2 dysregulation, thus providing insight about this cancer-associated, complex, and poorly understood genomic region.

Epidemiology, pathology, and genetics of histiocytic sarcoma in the Bernese mountain dog breed.

Histiocytic sarcoma (HS) refers to a highly aggressive and frequently disseminated neoplastic disease belonging to the class of canine histiocytic proliferative disorders. Disseminated HS (previously called malignant histiocytosis) is highly breed specific, with Bernese mountain dogs (BMDs), rottweilers, and retrievers having a high prevalence with a frequency of approximately 25% in the BMD breed. We collected DNA samples and clinical information from 800 BMDs, of which 200 are affected by HS. To better characterize the physiopathology and epidemiology, an in-depth analysis of 89 BMD cases has been performed. The mean age of onset was 6.5 years, males and females being equally affected. The clinical features, biochemical parameters, and pathological features have been determined. The life span after diagnosis has been estimated to be 49 days. A large BMD pedigree of 327 dogs, 121 of which are affected, was assembled. Using a subset of 160 BMDs, encompassing 21 complete sibships, we now propose an oligogenic transmission mode of the disease. Whole-genome linkage scans as well as association studies using a case/control analysis, in parallel with expression profiling of neoplastic versus normal histiocytes, are all underway. Altogether, these complementary approaches are expected to localize the genes for HS in the BMD, leading to advances in our knowledge of histiocyte diseases in dogs and humans.

Canine population structure: assessment and impact of intra-breed stratification on SNP-based association studies.

BACKGROUND: In canine genetics, the impact of population structure on whole genome association studies is typically addressed by sampling approximately equal numbers of cases and controls from dogs of a single breed, usually from the same country or geographic area. However one way to increase the power of genetic studies is to sample individuals of the same breed but from different geographic areas, with the expectation that independent meiotic events will have shortened the presumed ancestral haplotype around the mutation differently. Little is known, however, about genetic variation among dogs of the same breed collected from different geographic regions. METHODOLOGY/PRINCIPAL FINDINGS: In this report, we address the magnitude and impact of genetic diversity among common breeds sampled in the U.S. and Europe. The breeds selected, including the Rottweiler, Bernese mountain dog, flat-coated retriever, and golden retriever, share susceptibility to a class of soft tissue cancers typified by malignant histiocytosis in the Bernese mountain dog. We genotyped 722 SNPs at four unlinked loci (between 95 and 271 per locus) on canine chromosome 1 (CFA1). We showed that each population is characterized by distinct genetic diversity that can be correlated with breed history. When the breed studied has a reduced intra-breed diversity, the combination of dogs from international locations does not increase the rate of false positives and potentially increases the power of association studies. However, over-sampling cases from one geographic location is more likely to lead to false positive results in breeds with significant genetic diversity. CONCLUSIONS: These data provide new guidelines for association studies using purebred dogs that take into account population structure.

Facilitating genome navigation: survey sequencing and dense radiation-hybrid gene mapping.

Accurate and comprehensive sequence coverage for large genomes has been restricted to only a few species of specific interest. Lower sequence coverage (survey sequencing) of related species can yield a wealth of information about gene content and putative regulatory elements. But survey sequences lack long-range continuity and provide only a fragmented view of a genome. Here we show the usefulness of combining survey sequencing with dense radiation-hybrid (RH) maps for extracting maximum comparative genome information from model organisms. Based on results from the canine system, we propose that from now on all low-pass sequencing projects should be accompanied by a dense, gene-based RH map-construction effort to extract maximum information from the genome with a marginal extra cost.

A 1-Mb resolution radiation hybrid map of the canine genome.

The purebred dog population consists of >300 partially inbred genetic isolates or breeds. Restriction of gene flow between breeds, together with strong selection for traits, has led to the establishment of a unique resource for dissecting the genetic basis of simple and complex mammalian traits. Toward this end, we present a comprehensive radiation hybrid map of the canine genome composed of 3,270 markers including 1,596 microsatellite-based markers, 900 cloned gene sequences and ESTs, 668 canine-specific bacterial artificial chromosome (BAC) ends, and 106 sequence-tagged sites. The map was constructed by using the RHDF5000-2 whole-genome radiation hybrid panel and computed by using MULTIMAP and TSP/CONCORDE. The 3,270 markers map to 3,021 unique positions and define an average intermarker distance corresponding to 1 Mb. We also define a minimal screening set of 325 highly informative well spaced markers, to be used in the initiation of genome-wide scans. The well defined synteny between the dog and human genomes, established in part as a function of this work by the identification of 85 conserved fragments, will allow follow-up of initial findings of linkage by selection of candidate genes from the human genome sequence. This work continues to define the canine system as the method of choice in the pursuit of the genes causing mammalian variation and disease.