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.

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.

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.

Isolation and characterization of two canine melanoma cell lines: new models for comparative oncology.

BACKGROUND: Metastatic melanoma is one of the most aggressive forms of cancer in humans. Among its types, mucosal melanomas represent one of the most highly metastatic and aggressive forms, with a very poor prognosis. Because they are rare in Caucasian individuals, unlike cutaneous melanomas, there has been fewer epidemiological, clinical and genetic evaluation of mucosal melanomas. Moreover, the lack of predictive models fully reproducing the pathogenesis and molecular alterations of mucosal melanoma makes its treatment challenging. Interestingly, dogs are frequently affected by melanomas of the oral cavity that are characterized, as their human counterparts, by focal infiltration, recurrence, and metastasis to regional lymph nodes, lungs and other organs. In dogs, some particular breeds are at high risk, suggesting a specific genetic background and strong genetic drivers. Altogether, the striking homologies in clinical presentation, histopathological features, and overall biology between human and canine mucosal melanomas make dogs invaluable natural models with which to investigate tumor development, including tumor ætiology, and develop tailored treatments. METHODS: We developed and characterized two canine oral melanoma cell lines from tumors isolated from dog patients with distinct clinical profiles; with and without lung metastases. The cells were characterized using immunohistochemistry, pharmacology and genetic studies. RESULTS: We have developed and immunohistochemically, genetically, and pharmacologically characterized. Two cell lines (Ocr_OCMM1X & Ocr_OCMM2X) were produced through mouse xenografts originating from two clinically contrasting melanomas of the oral cavity. Their exhaustive characterization showed two distinct biological and genetic profiles that are potentially linked to the stage of malignancy at the time of diagnosis and sample collection of each melanoma case. These cell lines thus constitute relevant tools with which to perform genetic and drug screening analyses for a better understanding of mucosal melanomas in dogs and humans. CONCLUSIONS: The aim of this study was to establish and characterize xenograft-derived canine melanoma cell lines with different morphologies, genetic features and pharmacological sensitivities that constitute good predictive models for comparative oncology. These cell lines are relevant tools to advance the use of canine mucosal melanomas as natural models for the benefit of both veterinary and human medicine.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains.

BACKGROUND: Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation. RESULTS: We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains. CONCLUSIONS: In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

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.

ERK-regulated differential expression of the Mitf 6a/b splicing isoforms in melanoma.

The master regulator of the melanocyte lineage Mitf is intimately involved in development as well as melanoma, controlling cell survival, differentiation, proliferation and metastasis/migration. Consistent with its central role, Mitf expression and Mitf post-translational modifications are tightly regulated. An additional potential level of regulation is afforded by differential splicing of Mitf exon-6 leading to the generation of two isoforms that differ by the presence of six amino-acids in the Mitf (+) isoform and which have differential effects on cell cycle progression. However, whether the ratio of the two isoforms is regulated and whether their expression correlates with melanoma progression is not known. Here, we show that the differential expression of the Mitf 6a/b isoforms is dependent on the MAPKinase signalling, being linked to the activation of MEK1-ERK2, but not to N-RAS/B-RAF mutation status. In addition, quantification of Mitf 6a/b splicing forms in 86 melanoma samples revealed substantially increased levels of the Mitf (-) form in a subset of metastatic melanomas. The results suggest that differential expression of the Mitf 6a/b isoforms may represent an additional mechanism for regulating Mitf function and melanoma biology.

Analysis of splicing patterns by pyrosequencing.

Several different mRNAs can be produced from a given pre-mRNA by regulated alternative splicing, or as the result of deregulations that may lead to pathological states. Analysing splicing patterns is therefore of importance to describe and understand developmental programs, cellular responses to internal or external cues, or human diseases. We describe here a method, Pyrosequencing Analysis of Splicing Patterns (PASP), that combines RT-PCR and pyrosequencing of PCR products. We demonstrated that: (i) Ratios of two pure RNAs mixed in various proportions were accurately measured by PASP; (ii) PASP can be adapted to virtually any splicing event, including mutually exclusive exons, complex patterns of exon skipping or inclusion, and alternative 3' terminal exons; (iii) In extracts from different organs, the proportions of RNA isoforms measured by PASP reflected those measured by other methods. The PASP method is therefore reliable for analysing splicing patterns. All steps are done in 96-wells microplates, without gel electrophoresis, opening the way to high-throughput comparisons of RNA from several sources.

Target gene specificity of USF-1 is directed via p38-mediated phosphorylation-dependent acetylation.

How transcription factors interpret the output from signal transduction pathways to drive distinct programs of gene expression is a key issue that underpins development and disease. The ubiquitously expressed basic-helix-loop-helix leucine zipper upstream stimulating factor-1 binds E-box regulatory elements (CANNTG) to regulate a wide number of gene networks. In particular, USF-1 is a key component of the tanning process. Following UV irradiation, USF-1 is phosphorylated by the p38 stress-activated kinase on threonine 153 and directly up-regulates expression of the POMC, MC1R, TYR, TYRP-1 and DCT genes. However, how phosphorylation on Thr-153 might affect the activity of USF-1 is unclear. Here we show that, in response to DNA damage, oxidative stress and cellular infection USF-1 is acetylated in a phospho-Thr-153-dependent fashion. Phospho-acetylated USF-1 is nuclear and interacts with DNA but displays altered gene regulatory properties. Phospho-acetylated USF-1 is thus proposed to be associated with loss of transcriptional activation properties toward several target genes implicated in pigmentation process and cell cycle regulation. The identification of this critical stress-dependent USF-1 modification gives new insights into understanding USF-1 gene expression modulation associated with cancer development.

UV-induced expression of key component of the tanning process, the POMC and MC1R genes, is dependent on the p-38-activated upstream stimulating factor-1 (USF-1).

Protection against UV-mediated DNA damage and the onset of oncogenesis is afforded by the tanning response in which UV irradiation triggers melanocytes to increase production of melanin that is then transferred to keratinocytes. A key component of the tanning process is the UV-mediated induction of the pro-opiomelanocortin (POMC) and MC1R genes encoding the alpha-melanocyte-stimulating hormone and its receptor, respectively, which play a crucial role in pigmentation by regulating the intracellular levels of cAMP. How these genes are regulated in response to UV irradiation is not known. Here we have shown that UV-induced activation of the POMC and MC1R promoters is mediated by p38 stress-activated kinase signaling to the transcription factor, upstream stimulating factor-1 (USF-1). Importantly, melanocytes derived from USF-1 -/- mice exhibit a defective UV response and fail to activate POMC and MC1R expression in response to UV irradiation. The results define USF-1 as a critical UV-responsive activator of genes implicated in protection from solar radiation.

Plasmodium falciparum glycogen synthase kinase-3: molecular model, expression, intracellular localisation and selective inhibitors.

Worldwide increasing resistance of Plasmodium falciparum to common anti-malaria agents calls for the urgent identification of new drugs. Glycogen synthase kinase-3 (GSK-3) represents a potential screening target for the identification of such new compounds. We have cloned PfGSK-3, the P. falciparum gene homologue of GSK-3 beta. It encodes a 452-amino-acid, 53-kDa protein with an unusual N-terminal extension but a well-conserved catalytic domain. A PfGSK-3 tridimensional homology model was generated on the basis of the recently crystallised human GSK-3 beta. It illustrates how the regions involved in the active site, in substrate binding (P+4 phosphate binding domain) and in activity regulation are highly conserved. Recombinant PfGSK-3 phosphorylates GS-1, a GSK-3-specific peptide substrate, glycogen synthase, recombinant axin and the microtubule-binding protein tau. Neither native nor recombinant PfGSK-3 binds to axin. Expression and intracellular localisation of PfGSK-3 were investigated in the erythrocytic stages. Although PfGSK-3 mRNA is present in similar amounts at all stages, the PfGSK-3 protein is predominantly expressed at the early trophozoite stage. Once synthesized, PfGSK-3 is rapidly transported to the erythrocyte cytoplasm where it associates with vesicle-like structures. The physiological functions of PfGSK-3 for the parasite remain to be elucidated. A series of GSK-3 beta inhibitors were tested on both PfGSK-3 and mammalian GSK-3beta. Remarkably these enzymes show a partially divergent sensitivity to the compounds, suggesting that PfGSK-3 selective compounds might be identified.