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.

Donskoy cats as a new model of oculocutaneous albinism with the identification of a splice-site variant in Hermansky-Pudlak Syndrome 5 gene.

In the feline Donskoy breed, a phenotype that breeders call "pink-eye," with associated light-brown skin, yellow irises and red-eye effect, has been described. Genealogical data indicated an autosomal recessive inheritance pattern. A single candidate region was identified by genome-wide association study and SNP-based homozygosity mapping. Within that region, we further identified HPS5 (HPS5 Biogenesis Of Lysosomal Organelles Complex 2 Subunit 2) as a strong candidate gene, since HPS5 variants have been identified in humans and animals with Hermansky-Pudlak syndrome 5 or oculocutaneous albinism. A homozygous c.2571-1G>A acceptor splice-site variant located in intron 16 of HPS5 was identified in pink-eye cats. Segregation of the variant was 100% consistent with the inheritance pattern. Genotyping of 170 cats from 19 breeds failed to identify a single carrier in non-Donskoy cats. The c.2571-1G>A variant leads to HPS5 exon-16 splicing that is predicted to produce a 52 amino acids in-frame deletion in the protein. These results support an association of the pink-eye phenotype with the c.2571-1G>A variant. The pink-eye Donskoy cat extends the panel of reported HPS5 variants and offers an opportunity for in-depth exploration of the phenotypic consequences of a new HPS5 variant.

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.

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.

A radiation hybrid map of the European sea bass (Dicentrarchus labrax) based on 1581 markers: Synteny analysis with model fish genomes.

The selective breeding of fish for aquaculture purposes requires the understanding of the genetic basis of traits such as growth, behaviour, resistance to pathogens and sex determinism. Access to well-developed genomic resources is a prerequisite to improve the knowledge of these traits. Having this aim in mind, a radiation hybrid (RH) panel of European sea bass (Dicentrarchus labrax) was constructed from splenocytes irradiated at 3000 rad, allowing the construction of a 1581 marker RH map. A total of 1440 gene markers providing ~4400 anchors with the genomes of three-spined stickleback, medaka, pufferfish and zebrafish, helped establish synteny relationships with these model species. The identification of Conserved Segments Ordered (CSO) between sea bass and model species allows the anticipation of the position of any sea bass gene from its location in model genomes. Synteny relationships between sea bass and gilthead seabream were addressed by mapping 37 orthologous markers. The sea bass genetic linkage map was integrated in the RH map through the mapping of 141 microsatellites. We are thus able to present the first complete gene map of sea bass. It will facilitate linkage studies and the identification of candidate genes and Quantitative Trait Loci (QTL). The RH map further positions sea bass as a genetic and evolutionary model of Perciformes and supports their ongoing aquaculture expansion.

Drosophila translational elongation factor-1gamma is modified in response to DOA kinase activity and is essential for cellular viability.

Drosophila translational elongation factor-1gamma (EF1gamma) interacts in the yeast two-hybrid system with DOA, the LAMMER protein kinase of Drosophila. Analysis of mutant EF1gamma alleles reveals that the locus encodes a structurally conserved protein essential for both organismal and cellular survival. Although no genetic interactions were detected in combinations with mutations in EF1alpha, an EF1gamma allele enhanced mutant phenotypes of Doa alleles. A predicted LAMMER kinase phosphorylation site conserved near the C terminus of all EF1gamma orthologs is a phosphorylation site in vitro for both Drosophila DOA and tobacco PK12 LAMMER kinases. EF1gamma protein derived from Doa mutant flies migrates with altered mobility on SDS gels, consistent with it being an in vivo substrate of DOA kinase. However, the aberrant mobility appears to be due to a secondary protein modification, since the mobility of EF1gamma protein obtained from wild-type Drosophila is unaltered following treatment with several nonspecific phosphatases. Expression of a construct expressing a serine-to-alanine substitution in the LAMMER kinase phosphorylation site into the fly germline rescued null EF1gamma alleles but at reduced efficiency compared to a wild-type construct. Our data suggest that EF1gamma functions in vital cellular processes in addition to translational elongation and is a LAMMER kinase substrate in vivo.

Survey sequencing and radiation hybrid mapping to construct comparative maps.

Radiation hybrid (RH) mapping has become one of the most well-established techniques for economically and efficiently navigating genomes of interest. The success of the technique relies on random chromosome breakage of a target genome, which is then captured by recipient cells missing a preselected marker. Selection for hybrid cells that have DNA fragments bearing the marker of choice, plus a random set of DNA fragments from the initial irradiation, generates a set of cell lines that recapitulates the genome of the target organism several-fold. Markers or genes of interest are analyzed by PCR using DNA isolated from each cell line. Statistical tools are applied to determine both the linear order of markers on each chromosome, and the confidence of each placement. The resolution of the resulting map relies on many factors, most notably the degree of breakage from the initial radiation as well as the number of hybrid clones and mean retention value.A high-resolution RH map of a genome derived from low pass or survey sequencing (coverage from 1 to 2 times) can provide essentially the same comparative data on gene order that is derived from high-coverage (greater than x7) genome sequencing. When combined with fluorescence in situ hybridization, RH maps are complete and ordered blueprints for each chromosome. They give information about the relative order and spacing of genes and markers, and allow investigators to move between target and reference genomes, such as those of mouse or human, with ease although the approach is not limited to mammal genomes.

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.

The first radiation hybrid map of a perch-like fish: the gilthead seabream (Sparus aurata L).

Among Teleosts, Perciformes are the largest order of fishes and include numerous species of commercial importance. Perciformes also comprise species of primary interest for evolutionary studies and analysis of the sex determination systems and sex chromosome plasticity. Unfortunately, genomics tools and resources for Perciformes remain to be developed. Here, we report the production of a seabream whole-genome radiation hybrid (RH) panel in which quality was ascertained by the construction of a 2-Mb-resolution RH map. The map encompasses 440 markers (288 microsatellites, 82 gene-based markers, and 70 STS) suitable for linkage analysis and comparative mapping studies. Achievement of a RH panel and a whole-genome RH map should contribute to establishing seabream as a fish model among the Perciformes and should be of importance in aquaculture for marker-assisted selection, improvement of growth performance, and disease management. Development of RH maps in a cost-effective manner for other fishes with the described methodology will offer a powerful approach in aquaculture and will provide extended capabilities for comparing vertebrate genome evolution.

Construction of a 2-Mb resolution BAC microarray for CGH analysis of canine tumors.

Recognition of the domestic dog as a model for the comparative study of human genetic traits has led to major advances in canine genomics. The pathophysiological similarities shared between many human and dog diseases extend to a range of cancers. Human tumors frequently display recurrent chromosome aberrations, many of which are hallmarks of particular tumor subtypes. Using a range of molecular cytogenetic techniques we have generated evidence indicating that this is also true of canine tumors. Detailed knowledge of these genomic abnormalities has the potential to aid diagnosis, prognosis, and the selection of appropriate therapy in both species. We recently improved the efficiency and resolution of canine cancer cytogenetics studies by developing a small-scale genomic microarray comprising a panel of canine BAC clones representing subgenomic regions of particular interest. We have now extended these studies to generate a comprehensive canine comparative genomic hybridization (CGH) array that comprises 1158 canine BAC clones ordered throughout the genome with an average interval of 2 Mb. Most of the clones (84.3%) have been assigned to a precise cytogenetic location by fluorescence in situ hybridization (FISH), and 98.5% are also directly anchored within the current canine genome assembly, permitting direct translation from cytogenetic aberration to DNA sequence. We are now using this resource routinely for high-throughput array CGH and single-locus probe analysis of a range of canine cancers. Here we provide examples of the varied applications of this resource to tumor cytogenetics, in combination with other molecular cytogenetic techniques.

Dynamics of mammalian chromosome evolution inferred from multispecies comparative maps.

The genome organizations of eight phylogenetically distinct species from five mammalian orders were compared in order to address fundamental questions relating to mammalian chromosomal evolution. Rates of chromosome evolution within mammalian orders were found to increase since the Cretaceous-Tertiary boundary. Nearly 20% of chromosome breakpoint regions were reused during mammalian evolution; these reuse sites are also enriched for centromeres. Analysis of gene content in and around evolutionary breakpoint regions revealed increased gene density relative to the genome-wide average. We found that segmental duplications populate the majority of primate-specific breakpoints and often flank inverted chromosome segments, implicating their role in chromosomal rearrangement.

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 high-resolution comparative map of canine Chromosome 5q14.3-q33 constructed utilizing the 1.5x canine genome sequence.

A high-density map of the region of canine Chromosome 5 (CFA5) surrounding the evolutionary breakpoint between human Chromosomes 1p32 and 17pll was constructed by integrating a radiation hybrid map including 41 microsatellites, 10 BACs, and 59 genes and a linkage map including 18 markers. A collection of canine genomic survey sequences providing 1.5x coverage was used to identify dog orthologs of human genes, proving instrumental in the development of this map. Of particular interest is the canine BHD gene, within which we have previously described a single nucleotide polymorphism associated with Hereditary Multifocal Renal Cystadenocarcinoma and Nodular Dermatofibrosis (RCND) in German Shepherd dogs. The corresponding region of the human genome is particularly gene rich, containing genes involved in development, metabolism, and cancer that are likely to be of interest in future mapping studies. This current mapping effort on CFA5 expands the degree to which initial findings of linkage in canine families can be followed by successful positional cloning efforts and increases the value of the human genome sequence for defining candidate genes. Moreover, this study demonstrates the utility of genomic survey sequences when combined with accurate genome maps for rapid mapping of disease susceptibility loci.

CRH_Server: an online comparative and radiation hybrid mapping server for the canine genome.

UNLABELLED: CRH_Server is an on line Comparative and Radiation Hybrid mapping Server dedicated to canine genomics. CRH_Server allows users to compute their own RH data using the current canine RH map, and allows comparative dog/human mapping analyses. Finally, it suggests multiple options for storage and queries of the dog RH database. AVAILABILITY: http://idefix.univ-rennes1.fr:8080/Dogs/rh-server.html. SUPPLEMENTARY INFORMATION: All information is available at http://idefix.univ-rennes1.fr:8080/Dogs/help_rh-server.html.

A mutation in the canine BHD gene is associated with hereditary multifocal renal cystadenocarcinoma and nodular dermatofibrosis in the German Shepherd dog.

Hereditary multifocal renal cystadenocarcinoma and nodular dermatofibrosis (RCND) is a naturally occurring canine kidney cancer syndrome that was originally described in German Shepherd dogs. The disease is characterized by bilateral, multifocal tumors in the kidneys, uterine leiomyomas and nodules in the skin consisting of dense collagen fibers. We previously mapped RCND to canine chromosome 5 (CFA5) with a highly significant LOD score of 16.7 (theta=0.016). We have since narrowed the RCND interval following selection and RH mapping of canine genes from the 1.3 x canine genome sequence. These sequences also allowed for the isolation of gene-associated BACs and the characterization of new microsatellite markers. Ordering of newly defined markers and genes with regard to recombinants localizes RCND to a small chromosomal region that overlaps the human Birt-Hogg-Dubé locus, suggesting the same gene may be responsible for both the dog and the phenotypically similar human disease. We herein describe a disease-associated mutation in exon 7 of canine BHD that leads to the mutation of a highly conserved amino acid of the encoded protein. The absence of recombinants between the disease locus and the mutation in US and Norwegian dogs separated by several generations is consistent with this mutation being the disease-causing mutation. Strong evidence is provided that the RCND mutation may have a homozygous lethal effect (P<0.01).

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.

Characterization and radiation hybrid mapping of expressed sequence tags from the canine brain.

Maps of the canine genome are now developing rapidly. Most of the markers on the current integrated canine radiation hybrid/genetic linkage/cytogenetic map are highly polymorphic microsatellite (type II) markers that are very useful for mapping disease loci. However, there is still an urgent need for the mapping of gene-based (type I) markers that are required for comparative mapping, as well as identifying candidate genes for disease loci that have been genetically mapped. We constructed an adult brain cDNA library as a resource to increase the number of gene-based markers on the canine genome map. Eighty-one percent of the 2700 sequenced expressed sequence tags (ESTs) represented unique sequences. The canine brain ESTs were compared with sequences in public databases to identify putative canine orthologs of human genes. One hundred nine of the canine ESTs were mapped on the latest canine radiation hybrid (RH) panel to determine the location of the respective canine gene. The addition of these new gene-based markers revealed three conserved segments (CS) between human and canine genomes previously detected by fluorescence in situ hybridization (FISH), but not by RH mapping. In addition, five new CS between dog and human were identified that had not been detected previously by RH mapping or FISH. This work has increased the number of gene-based markers on the canine RH map by approximately 30% and indicates the benefit to be gained by increasing the gene content of the current canine comparative map.