FEELnc : a tool for long non-coding RNA annotation and its application to the dog transcriptome

Whole transcriptome sequencing (RNA-seq) has become a standard for cataloguing and monitoring RNA populations. One of the main bottlenecks, however, is to correctly identify the different classes of RNAs among the plethora of reconstructed transcripts, particularly those that will be translated (mRNAs) from the class of long non-coding RNAs (lncRNAs). Here, we present FEELnc (FlExible Extraction of LncRNAs), an alignment-free program that accurately annotates lncRNAs based on a Random Forest model trained with general features such as multi k-mer frequencies and relaxed open reading frames. Benchmarking versus five state-of-the-art tools shows that FEELnc achieves similar or better classification performance on GENCODE and NONCODE data sets. The program also provides specific modules that enable the user to fine-tune classification accuracy, to formalize the annotation of lncRNA classes and to identify lncRNAs even in the absence of a training set of non-coding RNAs. We used FEELnc on a real data set comprising 20 canine RNA-seq samples produced by the European LUPA consortium to substantially expand the canine genome annotation to include 10 374 novel lncRNAs and 58 640 mRNA transcripts. FEELnc moves beyond conventional coding potential classifiers by providing a standardized and complete solution for annotating lncRNAs and is freely available at https://github.com/tderrien/FEELnc.Peer reviewe

Differentiation of PC12 cells expressing estrogen receptor alpha: a new bioassay for endocrine-disrupting chemicals evaluation

International audienceXeno-estrogens, a class of endocrine disrupting chemicals (EDCs), can disturb estrogen receptor-dependent pathways involved in differentiation, proliferation or protection. Multiple methods have been developed to characterize the disturbances induced by EDCs in different cells or organs. In this study we have developed a new tool for the assessment of estrogenic compounds on differentiation. For this purpose we used the global model of NGF-induced neurite outgrowth of a pseudoneuronal PC12 cell line stably transfected with estrogen receptor alpha (PC12 ER). This new test evidences a new selectivity in which estradiol, genistein and 4-hydroxytamoxifen increased the NGF-induced neurite outgrowth of PC12 ER cells in a dose-dependent manner. In contrast, the strong estrogen agonist 17α-ethynylestradiol, the strong antagonist raloxifene and the agonist bisphenol A were unable to modify the neuritogenesis of PC12 ER cells. Therefore, the analysis of neuritogenesis in PC12 ER cells constitutes a complementary tool for the characterization of xeno-estrogen activity and also serves as a basis for further studies focusing on the mechanisms of EDCs in a neuronal context. Moreover, this test constitutes an alternative to animal testing

The tree that hides the forest: identification of common predisposing loci in several hematopoietic cancers and several dog breeds

Histiocytic sarcoma (HS) is a rare but aggressive cancer in humans and dogs. The spontaneous canine model, with the clinical, epidemiological and histological similarities with human HS and specific breed predispositions, is a unique model/opportunity to unravel the genetic bases 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 combining several breeds (Bernese mountain dog, Rottweiler, flat coated retriever and golden retriever) and three haematopoietic cancers (HS, lymphoma and mast cell tumor). Results showed that we not only refined the previously identified HS risk CDKN2A locus but we identified new loci on canine chromosomes 2, 5, 12, 14, 20, 26 and X. Capture and targeted sequencing of specific loci pointed towards the existence of regulatory variants in non coding regions and/or methylation mechanisms linked to risk haplotypes, leading to strong cancer predispositions in specific dog breeds. Our results showed that these canine cancer predisposing loci appear to be due to the additive effect of several risk haplotype involved also in other haematopoietic cancers such lymphoma or mast cell tumor, illustrating the pleiotropic nature of these canine cancer loci as observed in human oncology, thus reinforcing the interest of predisposed dog breeds to study cancer initiation and progression

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

International audienceHistiocytic 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. Author summary Because of specific breed structures and artificial selection, pet dogs suffer from numerous genetic diseases, including cancers and represent a unique spontaneous model of human cancers. Here, we focused on histiocytic sarcoma (HS), a rare and highly aggressive cancer in humans. In this study, we have used spontaneous affected and unaffected dogs from three predisposed dog breeds to identify loci involved in HS predisposition. Through genetic analyses, we showed that these canine cancer predispositions are due to the additive effect of several risk haplotypes also involved in the predisposition of other hematopoietic cancers. The corresponding chromosomal regions in humans are involved in the predisposition of several cancers and are also associated with immune traits. This study demonstrates the pleiotropic nature of these canine cancer loci as observed in human oncology, thereby reinforcing the interest of predisposed dog breeds to study mechanisms involved in cancer initiation

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

Characterisation and functional predictions of canine long non-coding RNAs

Abstract Long non-coding RNAs (lncRNAs) are a family of heterogeneous RNAs that play major roles in multiple biological processes. We recently identified an extended repertoire of more than 10,000 lncRNAs of the domestic dog however, predicting their biological functionality remains challenging. In this study, we have characterised the expression profiles of 10,444 canine lncRNAs in 26 distinct tissue types, representing various anatomical systems. We showed that lncRNA expressions are mainly clustered by tissue type and we highlighted that 44% of canine lncRNAs are expressed in a tissue-specific manner. We further demonstrated that tissue-specificity correlates with specific families of canine transposable elements. In addition, we identified more than 900 conserved dog-human lncRNAs for which we show their overall reproducible expression patterns between dog and human through comparative transcriptomics. Finally, co-expression analyses of lncRNA and neighbouring protein-coding genes identified more than 3,400 canine lncRNAs, suggesting that functional roles of these lncRNAs act as regulatory elements. Altogether, this genomic and transcriptomic integrative study of lncRNAs constitutes a major resource to investigate genotype to phenotype relationships and biomedical research in the dog species

Identification of a Missense Variant in MFSD12 Involved in Dilution of Phaeomelanin Leading to White or Cream Coat Color in Dogs.

White coat color in mammals has been selected several times during the domestication process. Numerous dog breeds are fixed for one form of white coat color that involves darkly pigmented skin. The genetic basis of this color, due to the absence of pigment in the hairs, was suggested to correspond to extreme dilution of the phaeomelanin, by both the expression of only phaeomelanin (locus E) and its extreme dilution (locus I). To go further, we performed genome-wide association studies (GWAS) using a multiple breed approach. The first GWAS, using 34 white dogs and 128 non-white dogs, including White Shepherds, Poodles, Cotons de Tulear and Bichons allowed us to identify two significantly associated loci on the locus E and a novel locus on chromosome 20. A second GWAS using 15 other breeds presenting extreme phaeomelanin dilution confirmed the position of locus I on the chromosome 20 (position 55 Mb = 6 × 10). Using whole-genome sequencing, we identified a missense variant in the first exon of MFSD12, a gene recently identified to be involved in human, mouse and horse pigmentation. We confirmed the role of this variant in phaeomelanin dilution of numerous canine breeds, and the conserved role of MFSD12 in mammalian pigmentation

An exome sequencing based approach for genome-wide association studies in the dog

Abstract Genome-wide association studies (GWAS) are widely used to identify loci associated with phenotypic traits in the domestic dog that has emerged as a model for Mendelian and complex traits. However, a disadvantage of GWAS is that it always requires subsequent fine-mapping or sequencing to pinpoint causal mutations. Here, we performed whole exome sequencing (WES) and canine high-density (cHD) SNP genotyping of 28 dogs from 3 breeds to compare the SNP and linkage disequilibrium characteristics together with the power and mapping precision of exome-guided GWAS (EG-GWAS) versus cHD-based GWAS. Using simulated phenotypes, we showed that EG-GWAS has a higher power than cHD to detect associations within target regions and less power outside target regions, with power being influenced further by sample size and SNP density. We analyzed two real phenotypes (hair length and furnishing), that are fixed in certain breeds to characterize mapping precision of the known causal mutations. EG-GWAS identified the associated exonic and 3′UTR variants within the FGF5 and RSPO2 genes, respectively, with only a few samples per breed. In conclusion, we demonstrated that EG-GWAS can identify loci associated with Mendelian phenotypes both within and across breeds

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

International audienceCanine 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

Interest of circulating tumor DNA as a biomarker for canine cancers: illustration in histiocytic sarcoma, oral malignant melanoma and multicentric lymphoma

Circulating tumor DNA (ctDNA) has become an attractive biomarker in human oncology and may be informative in cancer-affected dogs. By performing ddPCR or PARR methods, we detected tumor-specific point mutations, copy number alterations and chromosomal rearrangements in the plasma of cancer-affected dogs. It allowed the detection of ctDNA in 2/8 (25%) oral malignant melanoma cases, 12/13 (92.3%) lymphoma cases and 21/23 (91.3%) histiocytic sarcoma (HS) cases. The value of ctDNA to diagnose HS was explored in 133 dogs including 49 with HS. In this cohort, screening recurrent PTPN11 mutations in plasma had a specificity of 98.8%, and a sensitivity between 42.8-77% according to HS clinical presentation, being higher in internal forms, especially with pulmonary location. Regarding lymphoma, the follow-up of four dogs showed that the minimal residual disease detection by targeting lymphoma-specific antigen receptor rearrangement in the plasma was concordant with the clinical evaluation. Moreover, ctDNA analysis appeared interesting to assess treatment response and to predict relapse.This study shows that ctDNA is detectable in the plasma of cancer-affected dogs and is a relevant biomarker for diagnosis and clinical follow-up. With a growing interest in integrating natural canine tumors to explore new therapies, this biomarker appears promising in comparative oncology research