Next Article in Journal
The Supplementation of FloraMax-B11 Did Not Affect the Bile Acid Neosynthesis and the Enterohepatic Circulation in Broiler Chickens
Next Article in Special Issue
Epithelial-to-Mesenchymal Transition and Phenotypic Marker Evaluation in Human, Canine, and Feline Mammary Gland Tumors
Previous Article in Journal
Phylogenetic Relationships of the Strongyloid Nematodes of Australasian Marsupials Based on Mitochondrial Protein Sequences
Previous Article in Special Issue
The Implication of Serum Autoantibodies in Prognosis of Canine Mammary Tumors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

RT-qPCR Expression Profiles of Selected Oncogenic and Oncosuppressor miRNAs in Formalin-Fixed, Paraffin-Embedded Canine Mammary Tumors

by
Jessica Maria Abbate
1,
Alessia Giannetto
2,
Francesca Arfuso
1,
Barbara Brunetti
3,* and
Giovanni Lanteri
2
1
Department of Veterinary Sciences, University of Messina, Polo Universitario Annunziata, 98168 Messina, Italy
2
Department of Chemical, Biological, Pharmaceutical and Environmental Sciences, University of Messina, Polo Universitario Papardo, 98166 Messina, Italy
3
Department of Veterinary Medical Sciences, Alma Mater Studiorum University of Bologna, via Tolara di Sopra, Ozzano Emilia, 40064 Bologna, Italy
*
Author to whom correspondence should be addressed.
Animals 2022, 12(21), 2898; https://doi.org/10.3390/ani12212898
Submission received: 22 September 2022 / Revised: 17 October 2022 / Accepted: 21 October 2022 / Published: 22 October 2022

Abstract

:

Simple Summary

Early diagnosis of mammary gland cancer would allow intervention in the early stage of the disease, with a better prognosis for canine patients. MiRNAs are short non-transcribed RNA molecules that can be involved in cancer as molecular regulators of tumor development and progression and recognized as early, cancer-specific biomarkers with diagnostic and prognostic value and promising targeted for cancer therapy. Our results demonstrate that MiRNAs are differentially expressed in canine mammary tumors (CMTs) compared with normal mammary gland. In particular, oncogenic miR-18a, miR-18b and miR-21 were significantly overexpressed in malignant tumors and their involvement in receptor-mediated carcinogenesis and proteoglycan remodeling, making them candidate biomarkers with a prognostic value in CMTs. We also found downregulation of the oncosuppressor miR-146b in both benign and malignant CMTs compared with the normal mammary gland. MiR-146b may regulate the production of pro-inflammatory cytokines with a crucial role in cancer development and is predicted to target genes involved in Toll-like receptor and MAPK-signaling pathways. By virtue of their possible involvement in neoplastic transformation and progression, investigated miRNAs may represent candidate biomarkers with prognostic relevance in dogs with mammary tumors.

Abstract

MicroRNAs (miRNAs) can act as oncogenes or oncosuppressor genes, and their involvement in nearly all cancer-associated processes makes these small molecules promising diagnostic and prognostic biomarkers in cancer, as well as specific targets for cancer therapy. This study aimed to investigate the expression of 7 miRNAs (miR-18a, miR-18b, miR-22, miR-124, miR-145, miR-21, miR-146b) in formalin-fixed, paraffin-embedded canine mammary tumors (CMTs) by quantitative reverse transcription polymerase chain reaction (RT-qPCR). Twenty-six mammary samples were selected, including 22 CMTs (7 benign; 15 malignant) and 4 control samples (3 normal mammary gland and 1 case of lobular hyperplasia). Oncogenic miR-18a, miR-18b and miR-21 were significantly upregulated in malignant tumors compared with control tissues (p < 0.05). Conversely, oncosuppressor miR-146b was significantly downregulated in benign and malignant mammary tumors compared with control samples (p < 0.05) while, no group-related differences in the expression levels of miR-22, miR-124 and miR-145 were found (p > 0.05). Upregulated miRNAs found here, may regulate genes involved in receptor-mediated carcinogenesis and proteoglycan remodeling in cancer; while miRNA with reduced expression can regulate genes involved in Toll-like receptor and MAPK signaling pathways. According to the results obtained in the current study, the oncogenic and oncosuppressor miRNAs analyzed here are dysregulated in CMTs and the dysregulation of miRNA targets may lead to specific altered cellular processes and key pathways involved in carcinogenesis. Of note, since oncogenic miRNAs predicted to regulate neoplastic cell proliferation and hormonal activities, they may play an active role in neoplastic transformation and/or progression, having mechanistic and prognostic relevance in CMTs.

1. Introduction

Mammary tumors are the second most commonly diagnosed type of cancer in dogs, accounting for up to 70% of all tumors in hormonally intact females [1,2], thus representing a significant clinical problem. Almost half of canine mammary tumors (CMTs) are malignant [1,3] and 50% of carcinomas metastasizes to local lymph nodes, lungs, and occasionally invades bones [4,5], with high mortality rate (i.e., 20–55%) if dogs are not treated on time [6,7,8]. Histological examination still remains the gold standard method for accurate classification and grading of CMTs [9], and so far, invasive surgery is considered the only therapeutic treatment with a high success rate. Indeed, chemotherapy has not been shown to be as effective in treating CMTs as in human breast cancer [10,11], and standard chemotherapeutics do not significantly improve overall survival time [12]. Therefore, due to aggressive biological behaviour and limited therapeutic opportunities, early diagnosis of CMTs would allow intervention in the early stage of the disease, improving the overall survival time of canine patients.
MicroRNAs are short (18–22 nucleotides), non-transcribed, RNA molecules in cells that bind specific sequences in messenger RNA (mRNA) strands and direct gene silencing, contributing to a complex called “RNA-silencing complex” [13,14]. MiRNAs can regulate up to 30% of genes via post-transcriptional control in target cells and, of note, a single miRNA can target hundreds of different mRNAs, affecting a multitude of transcripts [15]. Noteworthy, miRNAs regulate virtually all cellular processes in both health and disease and play a key role in cancer-associated cellular processes [13]. Dysregulation of miRNA expression profiles is a common occurrence in cancer due to epigenetic modifications, transcription alteration or as a result of DNA copy number aberrations [16]. Consequently, the genetic instability of tumors leads to changes in miRNA-controlled gene regulatory networks [16,17,18] and promotes cancer-related signaling pathways [15]. MiRNAs can act as oncogenes (oncomiRs) and as tumor suppressor genes, and each group exerts its effect depending on its level of expression [16]. Commonly, miRNAs that negatively regulate the oncoproteins are downregulated in cancer, while miRNAs that negatively regulate tumor suppressor genes are upregulated [19]. Due to their role in tumorigenesis and their abundance in tumor microenvironment [20,21,22,23], miRNAs have been extensively studied as potential cellular and molecular biomarkers for prediction, diagnosis, progression monitoring, prognosis determination and as therapeutic targets for various cancers [24,25]. They act as tissue biomarkers, and free miRNAs are also available in most biological fluids, packaged into extracellular vesicles (i.e., exosomes), enwrapped with microvesicles during their biogenesis, and/or binding to proteins, protecting them from endogenous RNases in the extracellular environment [26], also making their sampling minimally invasive.
Several non-exosomal and exosome-derived miRNAs are already known to be dysregulated in canine mammary gland cancer, and noteworthy, high similarities in terms of miRNA dynamics and functions have been found between canine mammary and human breast cancer [19,27]. MiRNAs were investigated in primary and metastatic CMT tissues to find candidate biomarkers in distinguishing between different stages of CMTs [28]. Notably, miR-21 was found significantly upregulated in mammary carcinomas and, miR-29b, miR-101, miR-125a, miR-143 and miR-145 differed significantly in metastatic cells compared with primary cancer cells, making them potential biomarkers of metastasis [28]. More recently, a total of 124 miRNAs were found to differ significantly between metastatic and non-metastatic CMTs by RT-qPCR and RNA hybridization arrays and included in a miRNA signature that predicts shorter survival [29]. MiR-141 has been shown to be upregulated in several CMT cell-lines by RT-qPCR and experimentally validated to downregulate tumor suppressor p16/INK4a [30]. Finally, upregulated miR-19b, miR-18a and miR-124 have been found in serum from dogs with mammary carcinoma and suggested as candidate biomarkers for prognosis [31,32].
Dysregulation of miRNAs expression and their involvement in tumorigenesis have made them interesting biomarkers in cancer research and, miRNAs evaluation would represent a tool to find accurate biomarkers with diagnostic and prognostic value in CMTs and to discover miRNA-based drugs as promising targeted therapies to be used in the near future. Nowadays, investigations focusing on miRNA-based biomarkers have been conducted in canine mammary cancer using surgically removed CMT tissues [19,28,29,33,34], tumor cell lines [30,35,36] and serum/plasma samples [31,37], and data obtained are sometimes conflicting depending on the RNA sources and miRNA profiling and normalization techniques. Furthermore, bioinformatic analysis to identify miRNA-controlled target genes and miRNA-regulated cancer signaling pathways are performed only occasionally. In view of the above considerations, this study aims to investigate the expression levels of selected oncogenic and oncosuppressor miRNAs (miR-18a; miR-18b; miR-22; miR-124; miR-145; miR-21, miR-146b) supposed to be involved in the pathogenesis of canine mammary cancer by quantitative reverse transcription polymerase chain reaction (RT-qPCR), and to identify target genes and specific pathways regulated by differently expressed miRNAs. Furthermore, to the Authors’ knowledge, this is the first study that reports miRNAs expression in formalin-fixed, paraffin-embedded (FFPE) canine mammary tumors.

2. Materials and Methods

2.1. Mammary Samples

Canine mammary tumors were selected among the samples referred for diagnostic purposes between 2021 and 2022, identified in the archive of the AniCura Veterinary Hospital “I Portoni Rossi,” Bologna, Italy. Histological examination was performed by a board-certified pathologist (B.B.) on three-μm thick tissue sections stained with haematoxylin-eosin (HE). Classification of CMTs was performed according to a recent classification published by the Davis-Thompson DVM Foundation [38]. Histological grading was performed according to Peña et al. [6] and based on histological features of tubules formation, nuclear pleomorphism and mitotic count. Normal mammary gland tissue samples were collected from dogs undergoing postmortem examination and used as negative controls for molecular investigations.

2.2. RNA Extraction and Reverse Transcription Reaction

Two sections, 10–15-μm thick, of a 2-mm diameter selected area for each sample block, representative of tumors and normal mammary gland and identified on HE-stained histological slides, were sampled using a disposable biopsy punch and collected in a sterile 1.5 mL microcentrifuge tube. Purification of total RNA from FFPE tissue sections was performed using miRNeasy® FFPE Kit (QIAGEN, Milan, Italy; cat. no. 217504) and Deparaffinization Solution (QIAGEN, Milan, Italy; cat. no. 19093) according to the manufacturer’s instructions.
Reverse transcription (RT) reaction was performed using miRCURY® LNA® RT Kit (QIAGEN, Milan, Italy; cat. no. 339340) and provided UniSp6 RNA Spike-in was used as an internal quality control for the cDNA synthesis. Reverse transcription reaction was set up in 10-μL total reaction volume containing 2 μL of 5× miRCURY RT Reaction Buffer, 4.5 μL of RNase-free water, 1 μL of 10× miRCURY RT Enzyme Mix, 0.5 μL of synthetic RNA spike-in. Finally, 20 ng of extracted RNA in a 2 μL total volume was added to each RT reaction tube. Obtained cDNA was stored at −20 °C for further analysis. Concentration of nucleic acids (RNA; cDNA) were measured using the Nanodrop Spectrophotometer (NanoPhotometer N50, IMPLEN, USA).

2.3. Selected miRNAs and Real-Time Quantitative PCR (RT-qPCR)

The following miRCURY® LNA® miRNA PCR Assays (QIAGEN, Milan, Italy; Cat. no. 339306) were selected: gga-miR-18a-5p (assay ID-YP02100185); gga-miR-18b-5p (assay ID-YP02100265); hsa-miR-22-3p (assay ID–YP00204606); cbr-miR-124 (assay ID–YP02103368); has-miR-145-5p (assay ID–YP00204483); has-miR-21-5p (assay ID–YP00204230); has-miR-146b-5p (assay ID–YP02119310). Assays used as reference controls for miRNA normalization included: hsa-miR-16-5p (assay ID–YP00205702) and hsa-let-7a-5p (assay ID–YP00205727) [29,36].
RT-qPCR was performed using miRCURY LNA SYBR Green PCR Kit (QIAGEN, cat. no. 339345, 339346) following the manufacturer’s instructions. Quantitative reaction was performed in duplicate per each sample in a 20 μL total reaction volume, containing 10 μL of 2× miRCURY SYBR Green Master Mix, 2 μL of specific miRNA assay, 2 μL of RNase-free water and 6 μL of cDNA (60× diluted), and performed using a Rotor-Gene® Real-Time PCR system (QIAGEN). The cycling conditions were set as follows: 95 °C for 2 min, 40 cycles of 95 °C for 10 s and 56 °C for 60 s. MiRNA expression levels were presented in terms of fold change normalized to endogenous controls using the formula Δ Δ C q .

2.4. Gene Target Predictions, Gene Ontology and KEGG Pathway Enrichment

Predicted target genes for each miRNA were generated using the miRDB online resource and analysis platform (http://mirdb.org (accessed on 30 August 2022) and all targets having scores greater than 80 were selected, representing the most confident gene predictions [39]. All target genes obtained in the miRDB were uploaded into the DAVID database (https://david.ncifcrf.gov (accessed on 2 September 2022) using official gene symbols, and enriched gene ontology terms and KEGG pathways were identified [40,41,42].

2.5. Statistical Analysis

All data were tested for normal distribution using the Kolmogorov-Smirnov test. Since data resulted in abnormal distribution (p < 0.05), a nonparametric statistical analysis was applied. In particular, the Kruskal-Wallis test was performed to assess statistically significant differences among samples, followed by Dunn’s multiple comparison test. Statistical analyses were performed using GraphPad Prism version 4.00 (GraphPad Software, San Diego, CA, USA, 2003). The DAVID software provided raw and Benjamini-corrected p-values, and a crude threshold of p < 0.06 was selected for enriched KEGG pathways [43].

3. Results

3.1. Mammary Samples and Histological Classification of CMTs

A total of 26 mammary samples from dogs with a median age of 10 years (range: 4–12) were used in this study. In particular, 22 CMTs were used: 7 benign tumors (Group B) and 15 malignant tumors (Group M). Malignant tumors were: grade I, 8 samples (Group M1); grade II, 2 samples (Group M2); grade III, 5 samples (Group M3). Three samples of normal mammary gland and 1 sample of lobular hyperplasia were used as negative control cases (Group C). Histological classification and grading for CMTs is reported in Table 1.

3.2. miRNAs Expression Levels

Selected miRNAs were detected in all samples. Statistical analysis of data revealed a significant effect of group (p < 0.05) on the expression levels of miR-18a, miR-18b, miR-21 and miR-146b, whereas no group-related differences in the expression levels of miR-22, miR-124 and miR-145 were found (p > 0.05) (Figure 1 and Figure 2).
MiR-18a, miR-18b and miR-21 were upregulated in malignant tumors compared with the normal/hyperplastic mammary gland; furthermore, the expression level of both miR-18a and miR-21 was upregulated in malignant neoplasms compared with benign tumors. Conversely, miR-146b was significantly downregulated in both benign and malignant mammary tumors compared with control cases, and significant differences were also observed in grade II malignant tumors compared with both grade I and III, and benign tumors (Figure 2).

3.3. MicroRNAs Targets, Functional Annotation and Pathway Enrichment of All miRNAs

Target genes of statistically significant up- or downregulated miRNAs were predicted using miRDB online resource. The number of totals predicts identified per each miRNA ranged from 205 to 780 genes, and predicted genes having a target score greater than 80 were considered, including 45 target genes for both miR-18a and miR-18b, 69 for miR-21 and 47 for miR-146b. Detailed list of predicted target genes for significant up- and downregulated miRNAs is reported in Table S1.
Gene functional annotation and pathway enrichment analysis was performed using the DAVID bioinformatic tool. The complete list of target genes (n = 158) obtained in miRDB, was submitted in the DAVID database using official gene symbols and, 156 out of 158 uploaded targets were successfully mapped to canine genes, whereas two gene symbols (i.e., LOC483164; INSYN2) were identified as unknown. The gene functional annotation was carried out organizing the annotation types in three categories: biological process, cellular component, and molecular function.
Biological processes included genes mainly involved in the regulation of transcriptional activity and protein binding, and enriched biological term included the “Ubiquitin conjugation pathway” (p-value = 0.017) associated with genes FBXL3, TRAF6, KLHL15, PELI1, RMND5A, SOCS5. Cellular component items focused on nucleus (p-value = 0.019) with the involvement of the following genes: BBX, ELF2, ESF1, GATAD2B, HORMAD1, RORA, TRAF6, APPL1, CREBL2, ESR1, HSF5, HNRNPD, IRF2, KLHL15, MAGOH, NFIA, NFIB, PITX2, PPARA, PBRM1, UIMC1.
Molecular processes mainly involved genes implicated in DNA binding activity and regulation of transcription. Enriched molecular functions included: DNA-binding (p-value = 0.0017), associated with genes ELF2, RORA, THAP5, ESR1, HSF5, IRF2, NFIA, NFIB, PITX2, PPARA, STRBP, ZNF287; the function “Activator” (p-value = 0.027) associated with genes HIF1A, IRF2, NFIA, NFIB, PPARA.
Enriched pathways obtained from KEGG pathway enrichment analysis using DAVID Bioinformatics Resource are reported in Table 2.
The graphic representation of the enriched “Pathways in cancer” involving 10 targets genes for oncogenic and oncosuppressor miRNAs is reported in Figure 3. Pathway information is generated by KEGG using DAVID Bioinformatic tool (https://david.ncifcrf.gov (accessed on 2 September 2022).

3.4. Gene Ontology and Pathway Enrichment Differentiating between Up- and Downregulated miRNAs

A total of 112 predicts obtained for upregulated miR-18a, miR-18b and miR-21, were uploaded in the DAVID database, and 111 out of 112 were successfully mapped. The target genes controlled by the upregulated miRNAs were significantly enriched in the biological processes of “Ubiquitin conjugation pathway” (p-value = 0.021), associated with FBXL3, KLHL15, PELI1, RMND5A, SOCS5 genes; and “Transcription regulation” (p-value = 0.043) associated with genes ELF2, GATAD2B, ESR1, IRF2, PPARA, PBRM1, ZNF287. Enriched cellular components focused on the nucleus (p-value = 0.0087), associated with genes BBX, ELF2, ESF1, GATAD2B, RORA, CREBL2, ESR1, HSF5, IRF2, KLHL15, MAGOH, NFIA, NFIB, PITX2, PPARA, PBRM1, UIMC1. Enriched molecular terms included: DNA-binding (p-value = 0.00041) associated with ELF2, RORA, THAP5, ESR1, HSF5, IRF2, NFIA, NFIB, PITX2, PPARA, ZNF287; “Activator” (p-value= 0.009) associated with genes HIF1A, IRF2, NFIA, NFIB, PPARA.
Enriched pathways for upregulated miRNAs included: “Chemical carcinogenesis–receptor activation” (p-value = 0.018) associated with ESR1, FGF18, JAG1, KPNA6, PPARA genes; the pathway “Proteoglycans in cancer” (p-value = 0.089) associated with genes ESR1, FRS2, HIF1A, PDCD4.
Analysis of downregulated miR-146b was based on 46 successfully mapped predicted targets. Enriched biological terms included “Spermatogenesis- Differentiation” (p-value= 0.049) associated with genes HORMAD1 and STRBP. No enriched molecular functions and cellular components were found. Enriched pathways included: “Toll-like receptor signaling pathway” (p-value = 0.027) associated with genes CD80, TRAF6 and IRAK1; “MAPK signaling pathway” (p-value = 0.038) associated with genes MET, TRAF6, ERBB4 and IRAK1.

4. Discussion

MicroRNAs have functional roles in various diseases and in cancer, where they can act as oncogenes or tumor suppressor genes [16]. Differentially expressed miRNAs have been implicated in neoplastic conditions in dogs as molecular regulators of tumor development and progression, and circulatory miRNAs have been recognized as specific biomarkers for individual tumor type, including mammary gland tumors [28,29,30]. In the present study, the expression levels of 7 miRNAs presumably involved in the pathogenesis of CMTs were investigated and, to the best of the Authors’ knowledge, this is the first study that reports miRNAs expression level using FFPE mammary samples. The results obtained here suggest that the selected miRNAs are differentially expressed in neoplastic samples compared with the nonneoplastic mammary gland, with significantly up- and down-regulated miRNAs that could represent potential diagnostic and prognostic biomarkers, as well as promising targeted therapy in canine patients with mammary cancer. In particular, oncogenic miR-18a, miR-18b and miR-21 were significantly upregulated in mammary tumors compared with the normal mammary gland; whereas, the progressive decrease in the expression of oncosuppressor miR-146b was significantly related to malignancy.
Upregulation of both miR-18a and miR-18b in malignant tumors found here largely corroborates previous studies on miRNAs in CMTs [30,31,36]. In particular, miR-18a was found significantly overexpressed in malignant canine mammary epithelial cell lines and in their exosomes [30,36] as well as a significantly higher serum concentration of miR-18a was found in dogs with mammary carcinomas and histological evidence of lymphatic invasion, representing a possible candidate biomarker for prognosis [31]. Based on gene target prediction analysis, both miR-18a and miR-18b show a high probability of targeting to 3′ UTR of the estrogen receptor ESR1 (score 99 in miRDB) mRNA. Furthermore, miR-22 investigated in this study was also predicted to target ESR1 mRNA (miRDB score 84). However, unlike other studies [36], expression level of miR-22 in our study did not differ significantly between CMTs and nonneoplastic conditions. It has long been recognized that human and canine mammary neoplasms lose estrogen receptor (ER) expression along with increasing tumor grade [44,45], and a high target score for ESR1 mRNA for both miR-18a and miR-18b, together with their overexpression in malignant cases, indicate that specific miRNAs such as miR-18a and miR-18b may contribute to this loss of hormone receptor activity, providing a potential molecular mechanism for the loss of ER with increased grade of malignancy. Therefore, both miRNAs, may represent non-invasive biomarkers of hormonal status and phenotype with a prognostic value in canine mammary tumors.
Significant upregulation of miR-21 was found in malignant mammary cases compared with both benign tumors and control cases in our study. It has been well established that miR-21 overexpression is considered a hallmark of carcinogenesis in CMTs [19,46]. In particular, miR-21 selectively targets tumor suppressor genes and, acting as oncomiR through the inhibition of apoptosis in neoplastic cells [19,46]. Furthermore, in human breast tumors, miR-21 has been shown to simultaneously downregulate multiple metastasis-related tumor suppressor genes, including programmed cell death 4 (PDCD4), encouraging cell invasion and metastasis [46]. Although with a very low target score, gene target analysis performed in our study, included PDCD4 mRNA as a predicted target of miR-21 (score 81 in miRDB). Conversely, different from human breast cancer, miR-21 had a very high likelihood of targeting to 3′ UTR of the fibroblast growth factor 18 FGF18 (score 97 in miRDB) mRNA. In breast cancer, FGF18 promotes cell proliferation through the ERK/c-Myc signaling pathway and stimulates the production of epithelial-to-mesenchymal transition (EMT) factors, thus, promoting neoplastic cell migration and invasion [47]. Noteworthy, in the current study, the highest level of miR-21 expression was observed in the two malignant tumors with evidence of neoplastic emboli in lymphatic vessels (i.e., tubular carcinoma, grade III; inflammatory carcinoma, grade III) (Data not shown). In primary breast cancer, miR-21 has been included in a miRNA signature that predicts tumor recurrence and shorter survival [48], and, overexpression was positively correlated with metastatic carcinomas and poor prognosis, make it a good metastatic and prognostic biomarker in canine patients [28].
The Gene Ontology and pathway analysis in our study suggest that upregulated miRNAs are engaged in transcription regulation among enriched biological processes, inducing receptor activation and modulating expression patterns of proteoglycans and their biosynthetic activities. Therefore, upregulation of miRNAs target genes leads to deregulation of specific cellular and molecular processes probably involved in the pathogenesis of CMTs. Receptor-mediated carcinogenesis is a major mechanism of cancer induction, and in particular estrogen signaling pathway has been associated with mammary carcinogenesis [1,49]. Long-term exposure to the hormone estrogen and its binding to the estrogen receptor positively regulates the transcription of more proliferative genes and suppresses apoptotic genes, leading to the proliferation of cells with DNA mutations [49]. Regarding proteoglycans (PGs), they are the main biomolecules of the extracellular matrix and recent evidence shows that they represent an important target of microRNAs, and dysregulated expression of specific microRNAs results in aberrant expression patterns of proteoglycans and their biosynthetic activities [50]. Especially in malignant tumors, extensive remodeling of the tumor stroma is associated with marked alterations in proteoglycan expression and structural variability, which in turn affects tumor cell proliferation, migration and epithelial-to-mesenchymal transition [50]. These macromolecules mainly contribute to the formation of a temporary permissive matrix for tumor growth, influencing tissue organization, cell-cell and cell-matrix interactions and cell signaling, also regulating tumor stroma angiogenesis and the development of drug resistance [51].
Downregulation of miR-146b was identified in CMTs compared with nonneoplastic mammary tissues in the present study. In breast cancer cell, miR-146b seems to act as tumor suppressor genes [52], while, to the best of the Authors’ knowledge no studies investigate the role of this miRNA in canine mammary cancer, to date. In human breast cancer cells, miR-146b overexpression significantly inhibits cell proliferation and EMT by upregulating E-cadherin and downregulating vimentin expression, modulating the transcription of the EMT-related genes [52]. Furthermore, miR-146b inhibits NF-κB-mediated production of IL-6, a proinflammatory cytokine associated with cancer, STAT3 activity, and IL-6/STAT3-driven breast cancer cells migration and invasion [53]. Of note, as the gene encoding miR-146b is decreased in neoplastic cells, therapies that reintroduce or stimulate this miRNA have been suggested as beneficial to patients [53].
Target prediction analysis for miR-146b in our study suggests interleukin 1 receptor associated kinase 1 IRAK1 mRNA as the most predicted (miRDB target score 99). In triple-negative breast cancer, IRAK-1 regulates the expression of inflammatory genes by immune cells [54]. Since inflammatory cytokines interleukin 1 (IL-1) and interleukin 6 (IL-6) play an important role in the regulation of cancer development [55], the role of miR-146b in regulating cytokines production in CMTs is worthy of future investigations. Interestingly, enriched pathways for miR-146b included “Toll-like receptor signaling pathway” and Toll-like receptors (TLRs), modulate activation/ inhibition of immune and non-immune cells, primarily affecting the production of pro-inflammatory cytokines. Especially, TLR4 protein was found to be involved in breast cancer development, and its overexpression negatively correlates with prognosis [56]. Finally, abnormally activated MAPK-(mitogen-activated protein kinase) signaling pathway is a common finding in different type of cancer, including breast cancer, being involved in proliferation, migration and invasion characteristics of neoplastic cells [57,58].

5. Conclusions

Oncogenic miR-18a, miR-18b, miR-21, and oncosuppressor miR-146b are dysregulated in FFPE canine mammary tumors and, dysregulation of miRNAs targets leads to altered specific cellular processes and key pathways involved in carcinogenesis. Since investigated oncogenic miRNAs predicted to regulate neoplastic cell proliferation and hormone activities, they could play an active role in neoplastic transformation and/or progression of CMTs, having mechanistic and prognostic relevance. Among early biomarkers of cancer, miRNAs evaluation has been shown to be a very promising method for diagnosis and prognosis, as well as to discover miRNA-based drugs as targeted therapies. Further studies are needed to accurately identify targets genes and to confirm mechanistic and prognostic relevance of oncogenic miRNAs analyzed here, to be considered as valuable candidate biomarkers for canine patients ‘prognosis in the near future. Furthermore, the potential beneficial role of miR-146b also in dogs with mammary tumors is worthy of further investigations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani12212898/s1, Table S1: List of predicted target genes for dysregulated miRNAs generated using the miRDB online analysis platform.

Author Contributions

Conceptualization, J.M.A.; methodology, J.M.A., A.G.; F.A. and B.B.; software, J.M.A., A.G. and F.A.; validation, A.G. and B.B., formal analysis, J.M.A., A.G., F.A. and B.B.; data curation, J.M.A.; writing—original draft preparation, J.M.A.; writing—review and editing, A.G., F.A., B.B. and G.L.; supervision, A.G., B.B. and G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from the owners of the animals.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sorenmo, K. Canine mammary gland tumors. Vet. Clin. Am. Small Anim. Pract. 2003, 33, 573–596. [Google Scholar] [CrossRef]
  2. Moe, L. Population-based incidence of mammary tumors in some dog breeds. J. Reprod. Fertil. Suppl. 2001, 57, 439–443. [Google Scholar] [PubMed]
  3. Salas, Y.; Marquez, A.; Diaz, D.; Romero, L. Epidemiological study of mammary tumors in female dogs diagnosed during the period 2002-2012: A growing animal health problem. PLoS ONE 2015, 10, e0127381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Rosol, T.J.; Tannehill-Gregg, S.H.; LeRoy, B.E.; Mandl, S.; Contag, C.H. Animal models of bone metastasis. Cancer 2003, 97, 748–757. [Google Scholar] [CrossRef] [PubMed]
  5. Gilbertson, S.R.; Kurzman, I.D.; Zachrau, R.E.; Hurvitz, A.I.; Black, M.M. Canine Mammary Epithelial Neoplasms: Biologic Implications of Morphologic Characteristics Assessed in 232 Dogs. Vet. Pathol. 1983, 20, 127–142. [Google Scholar] [CrossRef] [Green Version]
  6. Peña, L.; De Andrés, P.J.; Clemente, M.; Cuesta, P.; Pérez-Alenza, M.D. Prognostic value of histological grading in noninflammatory canine mammary carcinomas in a prospective study with two-year follow-up: Relationship with clinical and histological characteristics. Vet. Pathol. 2013, 50, 94–105. [Google Scholar] [CrossRef] [PubMed]
  7. Philibert, J.C.; Snyder, P.W.; Glickman, N.; Glickman, L.T.; Knapp, D.W.; Waters, D.J. Influence of host factors on survival in dogs with malignant mammary gland tumors. J. Vet. Intern. Med. 2003, 17, 102–106. [Google Scholar] [CrossRef] [PubMed]
  8. Peña, L.; Nieto, A.I.; Pérez-Alenza, D.; Cuesta, P.; Castaño, M. Immunohistochemical detection of Ki-67 and PCNA in canine mammary tumors: Relationship to clinical pathological variables. J. Vet. Diagn. Investig. 1998, 10, 237–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Canadas, A.; França, M.; Pereira, C.; Vilaça, R.; Vilhena, H.; Tinoco, F.; Silva, M.J.; Ribeiro, J.; Medeiros, R.; Oliveira, P.; et al. Canine Mammary Tumors: Comparison of Classification and Grading Methods in a Survival Study. Vet. Pathol. 2019, 56, 208–219. [Google Scholar] [CrossRef] [PubMed]
  10. Levi, M.; Brunetti, B.; Sarli, G.; Benazzi, C. Immunohistochemical Expression of P-glycoprotein and Breast Cancer Resistance Protein in Canine Mammary Hyperplasia, Neoplasia and Supporting Stroma. J. Comp. Pathol. 2016, 155, 277–285. [Google Scholar] [CrossRef]
  11. Cassali, G.D.; Lavalle, G.E.; De Nardi, A.B.; Ferreira, E.; Bertagnolli, A.C.; Estrela-Lima, A.; Alessi, A.C.; Daleck, C.R.; Salgado, B.S.; Fernandes, C.G.; et al. Consensus for the Diagnosis, Prognosis and Treatment of Canine Mammary Tumors. Braz. J. Vet. Pathol. 2011, 4, 153–180. [Google Scholar]
  12. Simon, D.; Schoenrock, D.; Baumgärtner, W.; Nolfe, I. Postoperative Adjuvant Treatment of Invasive Malignant Mammary Gland Tumors in Dogs with Doxorubicin and Docetaxel. J. Vet. Intern. Med. 2006, 20, 1184–1190. [Google Scholar] [CrossRef]
  13. Cullen, J.M.; Breen, M. An Overview of Molecular Cancer Pathogenesis, Prognosis, and Diagnosis. In Tumors in Domestic Animals, 5th ed.; Meuten, D.J., Ed.; Wiley Blackwell: Ames, IA, USA, 2017; pp. 1–26. [Google Scholar]
  14. Xie, X.; Lu, J.; Kulbokas, E.J.; Golub, T.R.; Mootha, V.; Lindblad-Toh, K.; Lander, E.S.; Kellis, M. Systematic discovery of regulatory motifs in human promoters and 3’ UTRs by comparison of several mammals. Nature 2005, 434, 338–345. [Google Scholar] [CrossRef] [Green Version]
  15. Fabian, M.R.; Sonenberg, N. The mechanics of miRNA-mediated gene silencing: A look under the hood of miRISC. Nat. Struct. Mol. Bol. 2012, 19, 586–593. [Google Scholar] [CrossRef] [PubMed]
  16. Croce, C.M. Causes and consequences of microRNA dysregulation in cancer. Nat. Rev. Genet. 2009, 10, 704–714. [Google Scholar] [CrossRef]
  17. Kaszak, I.; Ruszczak, A.; Kanafa, S.; Kacprzak, K.; Król, M.; Jurka, P. Current biomarkers of canine mammary tumors. Acta Vet. Scand. 2018, 60, 66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Thayanithy, V.; Park, C.; Sarver, A.L.; Kartha, R.V.; Korpela, D.M.; Graef, A.J.; Steer, C.J.; Modiano, J.F.; Subramanian, S. Combinatorial treatment of DNA and chromatin-modifying drugs cause cell death in human and canine osteosarcoma cell lines. PLoS ONE 2012, 7, e43720. [Google Scholar] [CrossRef]
  19. Boggs, R.M.; Wright, Z.M.; Stickney, M.J.; Porter, W.W.; Murphy, K.E. MicroRNA expression in canine mammary cancer. Mamm. Genome 2008, 19, 561–569. [Google Scholar] [CrossRef] [Green Version]
  20. Uzuner, E.; Ulu, G.T.; Gürler, S.B.; Baran, Y. The role of MiRNA in Cancer: Pathogenesis, Diagnosis, and Treatment. In miRNomics. Methods in Molecular Biology; Allmer, J., Yousef, M., Eds.; Springer: New York, NY, USA, 2022; pp. 375–422. [Google Scholar]
  21. Shah, V.; Shah, J. Recent trends in targeting miRNAs for cancer therapy. J. Pharm. Pharmacol. 2020, 72, 1732–1749. [Google Scholar] [CrossRef]
  22. Tan, W.; Liu, B.; Qu, S.; Liang, G.; Luo, W.; Gong, C. MicroRNAs and cancer: Key paradigms in molecular therapy. Oncol. Lett. 2018, 15, 2735–2742. [Google Scholar]
  23. Peng, Y.; Croce, C.M. The role of MicroRNAs in human cancer. Signal Transduct. Target. Ther. 2016, 1, 15004. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Henry, N.L.; Hayes, D.F. Cancer biomarkers. Mol. Oncol. 2012, 6, 140–146. [Google Scholar] [CrossRef] [PubMed]
  25. Chatterjee, S.K.; Zetter, B.R. Cancer biomarkers: Knowing the present and predicting the future. Future Oncol. Lend. Eng. 2005, 1, 37–50. [Google Scholar] [CrossRef] [PubMed]
  26. Iorio, M.V. MicroRNA dysregulation in cancer: Diagnostics, monitoring and therapeutics. A comprehensive review. EMBO Mol. Med. 2012, 4, 143–159. [Google Scholar] [CrossRef] [PubMed]
  27. Petroušková, P.; Hudáková, N.; Maloveská, M.; Humeník, F.; Cizkova, D. Non-Exosomal and Exosome-Derived miRNAs as Promising in Canine Mammary Cancer. Life 2022, 12, 524. [Google Scholar] [CrossRef] [PubMed]
  28. Von Deetzen, M.C.; Schmeck, B.T.; Gruber, A.D.; Klopfleisch, R. Malignancy associated MicroRNA Expression Changes in Canine Mammary Cancer of Different Malignancies. ISRN Vet. Sci. 2014, 2014, 148597. [Google Scholar] [CrossRef] [Green Version]
  29. Bulkowska, M.; Rybicka, A.; Senses, K.M.; Ulewicz, K.; Witt, K.; Szymanska, J.; Taciak, B.; Klopfleisch, R.; Hellmén, E.; Dolka, I.; et al. MicroRNA expression patterns in canine mammary cancer show significant differences between metastatic and non-metastatic tumours. BMC Cancer 2017, 17, 728. [Google Scholar] [CrossRef] [Green Version]
  30. Lutful Kabir, F.M.; Delnnocentes, P.; Bird, R.C. Altered microRNA Expression Profiles and Regulation of INK4A/CDKN2A Tumor Suppressor Genes in Canine Breast Cancer Models. J. Cell. Biochem. 2015, 116, 2956–2969. [Google Scholar] [CrossRef]
  31. Fish, E.J.; Martinez-Romero, E.G.; Delnnocentes, P.; Koehler, J.W.; Prasad, N.; Smith, A.N.; Bird, R.C. Circulating microRNA as biomarkers of canine mammary carcinoma in dogs. J. Vet. Intern. Med. 2020, 34, 1282–1290. [Google Scholar] [CrossRef]
  32. Heishima, K.; Ichikawa, Y.; Yoshida, K.; Iwasaki, R.; Sakai, H.; Nakagawa, T.; Tanaka, Y.; Hoshino, Y.; Okamura, Y.; Murakami, M.; et al. Circulating microRNA-214 and -126 as potential biomarkers for canine neoplastic disease. Sci. Rep. 2017, 7, 2301. [Google Scholar] [CrossRef] [Green Version]
  33. Jeong, S.J.; Lee, K.H.; Nam, A.R.; Cho, J.Y. Genome-wide methylation profiling in canine mammary tumors reveals miRNA candidates associated with Human Breast Cancer. Cancers 2019, 11, 1466. [Google Scholar] [CrossRef] [Green Version]
  34. Losiewicz, K.; Chmielewska-Krzesinska, M.; Socha, P.; Jakimiuk, A.; Wasowicz, K. MiRNA-21, miRNA-10b, and miRNA-34a expression in canine mammary gland neoplasms. J. Vet. Res. 2014, 58, 447–451. [Google Scholar] [CrossRef]
  35. Osaki, T.; Suden, Y.; Sugiyama, A.; Azuma, K.; Murahata, Y.; Tsuka, T.; Ito, N.; Imagawa, T.; Okamoto, Y. Establishment of a canine mammary gland tmor cel line and characterization of its miRNA expression. J. Vet. Sci. 2016, 17, 385–390. [Google Scholar] [CrossRef] [PubMed]
  36. Fish, E.J.; Irizarry, K.J.; Delnnocentes, P.; Ellis, C.J.; Prasad, N.; Moss, A.G.; Bird, R.C. Malignant canine mammary epithelial cells shed exosomes containing differentially expressed microRNA that regulate oncogenic network. BMC Cancer 2018, 18, 832. [Google Scholar] [CrossRef] [PubMed]
  37. Ramadan, E.S.; Salem, N.Y.; Enam, I.A.; AbdElKader, N.A.; Farghali, H.A.; Khattab, M.S. MicroRNA-21 expression, serum tumor markers, and immunohistochemistry in canine mammary tumors. Vet. Res. Commun. 2022, 2, 377–388. [Google Scholar] [CrossRef]
  38. Zappulli, V.; Peña, L.; Rasotto, R.; Goldschmidt, M.H.; Gama, A.; Scruggs, J.L. Surgical Pathology of Tumors of Domestic Animals Volume 2: Mammary Tumors; Kiupel, M., Ed.; Davis-Thompson DVM Foundation: Washington, DC, USA, 2019; Volume 2. [Google Scholar]
  39. Wong, N.; Wang, X. miRDB: An online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015, 43, 46–52. [Google Scholar] [CrossRef] [Green Version]
  40. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022, 50, W216–W221. [Google Scholar] [CrossRef]
  41. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  42. Dennis, G.; Sherman, B.T.; Hosack, D.A.; Yang, J.; Gao, W.; Lane, H.C.; Lempicki, R.A. DAVID: Database for annotation, visualization, and integrated discovery. Genome Biol. 2003, 4, P3. [Google Scholar] [CrossRef] [Green Version]
  43. Irizarry, K.J.L.; Chan, A.; Kettle, D.; Kezian, S.; Ma, D.; Palacios, L.; Li, Q.Q.; Keeler, C.L.; Drechsler, Y. Bioinformatics Analysis of Chicken miRNAs Associated with Monocyte to Macrophage Differentiation and Subsequent IFNγ Stimulated Activation. Microrna 2017, 6, 53–70. [Google Scholar] [CrossRef]
  44. Chang, C.C.; Tsai, M.H.; Liao, J.W.; Chan, J.P.W.; Wong, M.L.; Chang, S.C. Evaluation of hormone receptor expression for use in predicting survival of female dogs with malignant mammary gland tumors. J. Am. Med. Assoc. 2009, 235, 391–396. [Google Scholar] [CrossRef] [PubMed]
  45. Chang, S.C.; Chang, C.C.; Chang, T.J.; Wong, M.L. Prognostic factors associated with survival two years after surgery in dogs with malignant mammary tumors: 79 cases (1998–2002). J. Am. Vet. Med. Assoc. 2005, 227, 1625–1629. [Google Scholar] [CrossRef] [PubMed]
  46. Zhu, S.; Wu, H.; Wu, F.; Nie, D.; Sheng, S.; Mo, Y.Y. MicroRNA-21 targets tumor suppressor genes in invasion and metastasis. Cell Res. 2008, 18, 350–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Yu, Z.; Lou, L.; Zhao, Y. Fibroblast growth factor 18 promotes the growth, migration and invasion of MDA-MB-231 cells. Oncol. Rep. 2018, 40, 704–714. [Google Scholar] [CrossRef] [Green Version]
  48. Papadaki, C.; Stratigos, M.; Markakis, G.; Spiliotaki, M.; Mastrostamatis, G.; Nikolaou, C.; Mavroudis, D.; Agelaki, S. Circulating microRNAs in the early prediction of disease recurrence in primary breast cancer. Breast Cancer Res. 2018, 20, 72. [Google Scholar] [CrossRef] [Green Version]
  49. Yue, W.; Yager, J.D.; Wang, J.P.; Jupe, E.R.; Santen, R.J. Estrogen receptor-dependent and independent mechanisms of breast cancer carcinogenesis. Steroids 2013, 78, 161–170. [Google Scholar] [CrossRef]
  50. Ibrahim, S.A.; Hassan, H.; Götte, M. MicroRNA regulation of proteoglycan function in cancer. FEBS J. 2014, 281, 5009–5022. [Google Scholar] [CrossRef]
  51. Theocharis, A.D.; Karamanos, N.K. Proteoglycans remodeling in cancer: Underlying molecular mechanisms. Matrix Biol. 2019, 75–76, 220–259. [Google Scholar] [CrossRef]
  52. Li, S.; Hao, J.; Hong, Y.; Mai, J.; Huang, W. Long Non-coding RNA NEAT1 Promotes the Proliferation, Migration, and Metastasis of Human Breast-Cancer Cells by Inhibiting miR-146b-5p Expression. Cancer Manag. Res. 2020, 12, 6091–6101. [Google Scholar] [CrossRef]
  53. Xiang, M.; Birkbak, N.J.; Vafaizadeh, V.; Walker, S.R.; Yeh, J.E.; Liu, S.; Kroll, Y.; Boldin, M.; Taganov, K.; Groner, B.; et al. STAT3 Induction of miR-146b forms a Feedback Loop to Inhibit the NF-kB to IL-6 Signaling Axis and STAT3-driven Cancer Phenotypes. Sci. Signal. 2014, 7, ra11. [Google Scholar] [CrossRef] [Green Version]
  54. Wee, Z.N.; Yatim, S.M.J.M.; Kohlbauer, V.K.; Feng, M.; Goh, J.Y.; Bao, Y.; Lee, P.L.; Zhang, S.; Wang, P.P.; Lim, E.; et al. IRAK1 is a therapeutic target that drives breast cancer metastasis and resistance to paclitaxel. Nat. Commun. 2015, 6, 8746. [Google Scholar] [CrossRef] [Green Version]
  55. Dmitrieva, O.S.; Shilovskiy, I.P.; Khaitov, M.R.; Grievennikov, S.I. Interleukins 1 and 6 as main mediators of inflammation of cancer. Biochemistry 2016, 81, 80–90. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, X.; Yu, X.; Wang, Q.; Lu, Y.; Chen, H. Expression and clinical significance of SABTB1 and TLR4 in breast cancer. Oncol. Lett. 2017, 14, 3611–3615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. MacCorkle, R.A.; Tan, T.H. Mitogen-activated protein kinases in cell-cycle control. Cell Biochem. Biophys. 2005, 43, 451–461. [Google Scholar] [CrossRef]
  58. Martin, G.S. Cell signaling and cancer. Cancer Cell 2003, 4, 167–174. [Google Scholar] [CrossRef]
Figure 1. Expression levels of miR-18a, miR-18b, miR-21, miR-22. MiR-18a, miR-18b and miR-21 were upregulated in malignant tumors (M1-3) compared with the control group (C). The expression level of miR-18a and miR-21 was significantly upregulated in malignant neoplasms compared with benign tumors (B). Significances (p < 0.05): a vs. Group C, b vs. Group B.
Figure 1. Expression levels of miR-18a, miR-18b, miR-21, miR-22. MiR-18a, miR-18b and miR-21 were upregulated in malignant tumors (M1-3) compared with the control group (C). The expression level of miR-18a and miR-21 was significantly upregulated in malignant neoplasms compared with benign tumors (B). Significances (p < 0.05): a vs. Group C, b vs. Group B.
Animals 12 02898 g001
Figure 2. Expression levels of miR-124, miR-145, miR-146b. MiR-146b expression was significantly downregulated in benign (B) and malignant (M1-3) tumors compared with control mammary gland (C). Significant upregulation for miR-146b expression was observed in malignant tumors grade II (M2) compared with grade I (M1) and III (M3) and benign tumors (B). Significances (p < 0.05): a vs. Group C; b vs. Group B, M1 and M3.
Figure 2. Expression levels of miR-124, miR-145, miR-146b. MiR-146b expression was significantly downregulated in benign (B) and malignant (M1-3) tumors compared with control mammary gland (C). Significant upregulation for miR-146b expression was observed in malignant tumors grade II (M2) compared with grade I (M1) and III (M3) and benign tumors (B). Significances (p < 0.05): a vs. Group C; b vs. Group B, M1 and M3.
Animals 12 02898 g002
Figure 3. Enriched “Pathways in cancer” in Canine Mammary Tumors. Predicted target genes for oncogenic and oncosuppressor miRNAs are marked with red star, including 5 target genes for miR-146b (MET; TRAF6; APPL1; F2; IL12A); 3 predicted genes for both miR-18a and miR-18b (ESR1; HIF1A; ITGA6) and 2 predicts for miR-21 (FGF18; JAG1).
Figure 3. Enriched “Pathways in cancer” in Canine Mammary Tumors. Predicted target genes for oncogenic and oncosuppressor miRNAs are marked with red star, including 5 target genes for miR-146b (MET; TRAF6; APPL1; F2; IL12A); 3 predicted genes for both miR-18a and miR-18b (ESR1; HIF1A; ITGA6) and 2 predicts for miR-21 (FGF18; JAG1).
Animals 12 02898 g003
Table 1. Histological Classification and Grading for Canine Mammary Tumors.
Table 1. Histological Classification and Grading for Canine Mammary Tumors.
Tumor TypeHistological ClassificationNo.Grade of MalignancyLymphatic Invasion
Benign Tumor (n = 7)Tubular Adenoma1
Tubulopapillary Adenoma1
Complex Adenoma2
Intraductal Papillary Adenoma2
Benign Mixed Tumor1
Malignant Tumor (n = 15)Tubular Carcinoma1
1
I
III
Yes
Tubulopapillary Carcinoma1I
Complex Carcinoma1
1
1
I
I
I
Intraductal Papillary Carcinoma1
1
I
I
Mixed Carcinoma1
1
I
II
Solid Carcinoma1
1
1
1
II
III
III
III
Inflammatory Carcinoma1IIIYes
Table 2. Pathway Enrichment for all statistically significant expressed miRNAs.
Table 2. Pathway Enrichment for all statistically significant expressed miRNAs.
KEGG PathwayGenesNumber of Genes (%)p-ValueFold EnrichmentCorrected p-Value
Neurotrophin signaling pathwayTRAF6, FRS2, IRAK1, MAP3K1, NTF35 (3.2)0.0145.30.82
Alcoholic liver diseaseTRAF6, CPT1C, IRAK1, IL12A, PPARA5 (3.2)0.0204.70.82
Proteoglycans in cancerMET, ERBB4, ESR1, FRS2, HIF1A, PDCD46 (3.8)0.0233.60.82
Pathways in cancerMET, TRAF6, APPL1, F2, ESR1, FGF18, HIF1A, ITGA6, IL12A, JAG110 (6.4)0.0242.30.82
MAPK signaling pathwayMET, TRAF6, ERBB4, FGF18, IRAK1, MAP3K1, NTF37 (4.5)0.0253.00.82
Toll-like receptor signaling pathwayCD80, TRAF6, IRAK1, IL12A4 (2.6)0.0415.11.0
ToxoplasmosisTRAF6, ITGA6, IRAK1, IL12A4 (2.6)0.0524.71.0
Chemical carcinogenesis–receptor activationESR1, FGF18, JAG1, KPNA6, PPARA5 (3.2)0.0703.21.0
RIG-I-like receptor signaling pathwayTRAF6, IL12A, MAP3K13 (1.9)0.0935.81.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abbate, J.M.; Giannetto, A.; Arfuso, F.; Brunetti, B.; Lanteri, G. RT-qPCR Expression Profiles of Selected Oncogenic and Oncosuppressor miRNAs in Formalin-Fixed, Paraffin-Embedded Canine Mammary Tumors. Animals 2022, 12, 2898. https://doi.org/10.3390/ani12212898

AMA Style

Abbate JM, Giannetto A, Arfuso F, Brunetti B, Lanteri G. RT-qPCR Expression Profiles of Selected Oncogenic and Oncosuppressor miRNAs in Formalin-Fixed, Paraffin-Embedded Canine Mammary Tumors. Animals. 2022; 12(21):2898. https://doi.org/10.3390/ani12212898

Chicago/Turabian Style

Abbate, Jessica Maria, Alessia Giannetto, Francesca Arfuso, Barbara Brunetti, and Giovanni Lanteri. 2022. "RT-qPCR Expression Profiles of Selected Oncogenic and Oncosuppressor miRNAs in Formalin-Fixed, Paraffin-Embedded Canine Mammary Tumors" Animals 12, no. 21: 2898. https://doi.org/10.3390/ani12212898

APA Style

Abbate, J. M., Giannetto, A., Arfuso, F., Brunetti, B., & Lanteri, G. (2022). RT-qPCR Expression Profiles of Selected Oncogenic and Oncosuppressor miRNAs in Formalin-Fixed, Paraffin-Embedded Canine Mammary Tumors. Animals, 12(21), 2898. https://doi.org/10.3390/ani12212898

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop