Free circulating versus extracellular vesicle-associated microRNA expression in canine T-cell lymphoma.

Introduction: Canine lymphoma (cL) is one of the most frequent cancers in dogs. The T-cell lymphoma (TcL) is not the most common phenotype but presents an aggressive behavior. MicroRNAs (miRNAs), are small, single-stranded, non-coding RNA molecules which can circulate freely in blood or be associated with extracellular vesicles (EVs). The dysregulation of certain miRNAs has been identified in numerous types of human cancers and they have been largely investigated as possible tumors biomarkers in human medicine, while research in veterinary oncology is still scarce. The aim of this study was to compare the expression patterns of free circulating and EV-associated miRNAs in dogs with T-cell lymhoma (TcL) and healthy dogs. Methods: Eight dogs with TcL were selected as the lymphoma group (LG) and eight dogs were included as controls (Ctrl). Plasma samples were collected at the time of the diagnosis and EVs isolated with ultracentrifugation. miRNAs were extracted from both the circulating EVs and the plasma supernatant, obtaining EV-associated and free-miRNAs. Quantitative real-time PCR was performed to analyze the expression of 88 target miRNAs. Results: Ten and seven differentially expressed miRNAs between LG and Ctrl were detected in EV-associated and free-miRNAs, respectively. Among EV-associated and free-miRNAs, only has-miR-222-3p was overexpressed in both conditions. Discussion: All the differentially expressed miRNAs detected in this study, have been already described as dysregulated in other human or canine cancers. The EV-associated miRNAs, which appear to be more stable and better conserved than free-miRNAs, could be investigated in further larger studies to better assess their use as possible biomarkers for TcL.

Integrated analysis of transcriptome, methylome and copy number aberrations data of marginal zone lymphoma and follicular lymphoma in dog.

Marginal zone lymphoma (MZL) and follicular lymphoma (FL) are classified as indolent B-cell lymphomas in dogs. Aside from the clinical and histopathological similarities with the human counterpart, the molecular pathogenesis remains unclear. We integrated transcriptome, genome-wide DNA methylation and copy number aberration analysis to provide insights on the pathogenesis of canine MZL (n = 5) and FL (n = 7), also comparing them with diffuse large B-cell lymphoma (DLBCL). Transcriptome profiling highlighted the presence of similar biological processes affecting both histotypes, including BCR and TLR signalling pathways. However, FLs showed an enrichment of E2F targets, whereas MZLs were characterized by MYC-driven transcriptional activation signatures. FLs showed a distinctive loss on chr1 containing CEACAM23 and 24, conversely MZLs presented multiple recurrent gains on chr13, where MYC is located. The distribution of methylation peaks was similar between the two histotypes. Integrating data from the three omics, FLs resulted clearly separated from MZLs and DLBCL dataset. MZLs showed the enrichment of FoxM1 network and TLR associated TICAM1-dependent IRFs activation pathway. However, no specific signatures differentiated MZLs from DLBCLs. In conclusion, our study presents the first comprehensive analysis of molecular and epigenetic pathogenesis of canine FL and MZL.

Cyclin D1 immunohistochemical expression and somatic mutations in canine oral melanoma.

Canine oral melanoma (COM) is the most frequent tumour with oral localization in dogs. Copy number gains and amplifications of CCND1, a gene coding for Cyclin D1, are the most frequent chromosomal aberrations described in human non-UV induced melanomas. Twenty-eight cases of COM were retrieved from paraffin-blocks archives. A total of 4 µm thick sections were immunostained with an antibody against human Cyclin D1 and Ki-67. Cyclin D1 and Ki-67 expressions were scored through two counting methods. DNA was extracted from 20 µm thick sections of formalin-fixed paraffin-embedded blocks. Pathological and surrounding healthy tissue was extracted independently. Cyclin D1 immunolabelling was detected in 69% (18/26) while Ki-67 was present in 88.5% (23/26) of cases. Statistical analysis revealed correlation between two counting methods for Cyclin D1 (r = 0.54; P = .004) and Ki-67 (r = 0.56; P = .003). The correlation found between Ki-67 and Cyclin D1 indexes in 16/26 cases labelled by both antibodies (r = 0.7947; P = .0002) suggests a possible use of Cyclin D1 index as prognostic marker. Polymerase chain reaction analysis on CCND1 coding sequence revealed the presence of nine somatic mutations in seven samples producing synonymous, missense and stop codons. Since none of the single-nucleotide polymorphisms was found to be recurrent, it is suggested that overexpression of Cyclin D1 may be the consequence of alterations of CCND1 upstream regions or other genetic aberrations not detectable with the methodology used in this study. Future studies are needed to verify the potential use of Cyclin D1 index as prognostic indicator and to highlight the molecular events responsible for Cyclin D1 overexpression in COMs.

Long Non-Coding RNAs as Molecular Signatures for Canine B-Cell Lymphoma Characterization.

BACKGROUND: Diffuse large B-cell lymphoma (DLBCL), marginal zone lymphoma (MZL) and follicular lymphoma (FL) are the most common B-cell lymphomas (BCL) in dogs. Recent investigations have demonstrated overlaps of these histotypes with the human counterparts, including clinical presentation, biologic behavior, tumor genetics, and treatment response. The molecular mechanisms that underlie canine BCL are still unknown and new studies to improve diagnosis, therapy, and the utilization of canine species as spontaneous animal tumor models are undeniably needed. Recent work using human DLBCL transcriptomes has suggested that long non-coding RNAs (lncRNAs) play a key role in lymphoma pathogenesis and pinpointed a restricted number of lncRNAs as potential targets for further studies. RESULTS: To expand the knowledge of non-coding molecules involved in canine BCL, we used transcriptomes obtained from a cohort of 62 dogs with newly-diagnosed multicentric DLBCL, MZL and FL that had undergone complete staging work-up and were treated with chemotherapy or chemo-immunotherapy. We developed a customized R pipeline performing a transcriptome assembly by multiple algorithms to uncover novel lncRNAs, and delineate genome-wide expression of unannotated and annotated lncRNAs. Our pipeline also included a new package for high performance system biology analysis, which detects high-scoring network biological neighborhoods to identify functional modules. Moreover, our customized pipeline quantified the expression of novel and annotated lncRNAs, allowing us to subtype DLBCLs into two main groups. The DLBCL subtypes showed statistically different survivals, indicating the potential use of lncRNAs as prognostic biomarkers in future studies. CONCLUSIONS: In this manuscript, we describe the methodology used to identify lncRNAs that differentiate B-cell lymphoma subtypes and we interpreted the biological and clinical values of the results. We inferred the potential functions of lncRNAs to obtain a comprehensive and integrative insight that highlights their impact in this neoplasm.

Mutational landscape of canine B-cell lymphoma profiled at single nucleotide resolution by RNA-seq.

The genomic landscape in human B-cell lymphoma has revealed several somatic mutations and potentially relevant germline alterations affecting therapy and prognosis. Also, mutations originally described as somatic aberrations have been shown to confer cancer predisposition when occurring in the germline. The relevance of mutations in canine B-cell lymphoma is scarcely known and gene expression profiling has shown similar molecular signatures among different B-cell histotypes, suggesting other biological mechanisms underlining differences. Here, we present a highly accurate approach to identify single nucleotide variants (SNVs) in RNA-seq data obtained from 62 completely staged canine B-cell lymphomas and 11 normal B-cells used as controls. A customized variant discovery pipeline was applied and SNVs were found in tumors and differentiated for histotype. A number of known and not previously identified SNVs were significantly associated to MAPK signaling pathway, negative regulation of apoptotic process and cell death, B-cell activation, NF-kB and JAK-STAT signaling. Interestingly, no significant genetic fingerprints were found separating diffuse large B-cell lymphoma from indolent lymphomas suggesting that differences of genetic landscape are not the pivotal causative factor of indolent behavior. We also detected several variants in expressed regions of canine B-cell lymphoma and identified SNVs having a direct impact on genes. Using this brand-new approach the consequence of a gene variant is directly associated to expression. Further investigations are in progress to deeply elucidate the mechanisms by which altered genes pathways may drive lymphomagenesis and a higher number of cases is also demanded to confirm this evidence.

Longitudinal transcriptomic and genetic landscape of radiotherapy response in canine melanoma.

Canine malignant melanoma (MM) is a highly aggressive tumour with a low survival rate and represents an ideal spontaneous model for the human counterpart. Considerable progress has been recently obtained, but the therapeutic success for canine melanoma is still challenging. Little is known about the mechanisms beyond pathogenesis and melanoma development, and the molecular response to radiotherapy has never been explored before. A faster and deeper understanding of cancer mutational processes and developing mechanisms are now possible through next generation sequencing technologies. In this study, we matched whole exome and transcriptome sequencing in four dogs affected by MM at diagnosis and at disease progression to identify possible genetic mechanisms associated with therapy failure. According to previous studies, a genetic similarity between canine MM and its human counterpart was observed. Several somatic mutations were functionally related to MAPK, PI3K/AKT and p53 signalling pathways, but located in genes other than BRAF, RAS and KIT. At disease progression, several mutations were related to therapy effects. Natural killer cell-mediated cytotoxicity and several immune-system-related pathways resulted activated opening a new scenario on the microenvironment in this tumour. In conclusion, this study suggests a potential role of the immune system associated to radiotherapy in canine melanoma, but a larger sample size associated with functional studies are needed.

Array Comparative Genomic Hybridization Analysis Reveals Significantly Enriched Pathways in Canine Oral Melanoma.

Human Mucosal Melanoma (hMM) is an aggressive neoplasm of neuroectodermal origin with distinctive features from the more common cutaneous form of malignant melanoma (cMM). At the molecular level, hMMs are characterized by large chromosomal aberrations rather than single-nucleotide mutations, typically observed in cMM. Given the scarcity of available cases, there have been many attempts to establish a reliable animal model. In pet dogs, Canine Oral Melanoma (COM) is the most common malignant tumor of the oral cavity, sharing clinical and histological aspects with hMM. To improve the knowledge about COM's genomic DNA alterations, in the present work, formalin-fixed, paraffin-embedded (FFPE) samples of COM from different European archives were collected to set up an array Comparative Genomic Hybridization (aCGH) analysis to estimate recurrent Copy Number Aberrations (CNAs). DNA was extracted in parallel from tumor and healthy fractions and 19 specimens were successfully submitted to labeling and competitive hybridization. Data were statistically analyzed through GISTIC2.0 and a pathway-enrichment analysis was performed with ClueGO. Recurrent gained regions were detected, affecting chromosomes CFA 10, 13 and 30, while lost regions involved chromosomes CFA 10, 11, 22, and 30. In particular, CFA 13 showed a whole-chromosome gain in 37% of the samples, while CFA 22 showed a whole-chromosome loss in 25%. A distinctive sigmoidal trend was observed in CFA 10 and 30 in 25 and 30% of the samples, respectively. Comparative analysis revealed that COM and hMM share common chromosomal changes in 32 regions. MAPK- and PI3K-related genes were the most frequently involved, while pathway analysis revealed statistically significant perturbation of cancer-related biological processes such as immune response, drug metabolism, melanocytes homeostasis, and neo-angiogenesis. The latter is a new evidence of a significant involvement of neovascularization-related pathways in COMs and can provide the rationale for future application in anti-cancer targeted therapies.

Validation of epigenetic mechanisms regulating gene expression in canine B-cell lymphoma: An in vitro and in vivo approach.

Despite canine B-cell Lymphoma (BCL) representing the most common haematological tumour, epigenetic events driving development and progression are scarcely known. Recently, canine Diffuse Large BCL (DLBCL) DNA methylome by genome-wide CpG microarray has identified genes and pathways associated to pathogenesis. To validate data previously obtained by array analysis, the CLBL-1 cell line was used and the HOXD10, FGFR2, ITIH5 and RASAL3 genes were selected. CLBL-1 cells were treated with two hypomethylating drugs (HDs; IC50, 50% inhibitory concentration), i.e. azacytidine and decitabine (DEC), either alone or in combination with three histone deacetylase inhibitors (HDACis; IC20), i.e. valproic acid, trichostatin and vorinostat. Following the incubation with both HDs, an overall decrease of promoter methylation was highlighted, thus confirming target genes hypermethylation. The highest mRNA restoration was observed following the exposure to HDs combined with HDACis, and mostly with valproic acid. Contrasting results were only obtained for RASAL3. An in vivo confirmation was finally attempted treating Nod-Scid mice engrafted with CLBL-1 cells with DEC. Although DEC did not arrest tumour growth, target genes promoter methylation was significantly reduced in DEC-treated mice vs controls. Overall, this work demonstrates that CLBL-1 cell line represents a reliable in vitro model to validate the methylation-dependent silencing of key genes for BCL; moreover, it may be useful for xenograft models in mice, despite its aggressive behaviour. In future, functional studies will be performed to deepen the role of selected genes on BCL pathogenesis and progression, and their methylation-dependent mechanism of regulation.

DNA methylation profiling reveals common signatures of tumorigenesis and defines epigenetic prognostic subtypes of canine Diffuse Large B-cell Lymphoma.

Epigenetic deregulation is a hallmark of cancer characterized by frequent acquisition of new DNA methylation in CpG islands. To gain insight into the methylation changes of canine DLBCL, we investigated the DNA methylome in primary DLBCLs in comparison with control lymph nodes by genome-wide CpG microarray. We identified 1,194 target loci showing different methylation levels in tumors compared with controls. The hypermethylated CpG loci included promoter, 5'-UTRs, upstream and exonic regions. Interestingly, targets of polycomb repressive complex in stem cells were mostly affected suggesting that DLBCL shares a stem cell-like epigenetic pattern. Functional analysis highlighted biological processes strongly related to embryonic development, tissue morphogenesis and cellular differentiation, including HOX, BMP and WNT. In addition, the analysis of epigenetic patterns and genome-wide methylation variability identified cDLBCL subgroups. Some of these epigenetic subtypes showed a concordance with the clinical outcome supporting the hypothesis that the accumulation of aberrant epigenetic changes results in a more aggressive behavior of the tumor. Collectively, our results suggest an important role of DNA methylation in DLBCL where aberrancies in transcription factors were frequently observed, suggesting an involvement during tumorigenesis. These findings warrant further investigation to improve cDLBCL prognostic classification and provide new insights on tumor aggressiveness.