The miRNome of canine invasive urothelial carcinoma

Abstract

Urothelial carcinoma (UC) comprises up to 2% of all naturally occurring neoplasia in dogs and can be challenging to diagnose. MicroRNAs (miRNAs) have been reported to be dysregulated in numerous diseases, including neoplasia. MiRNA expression has been evaluated in human UC, but there is limited information regarding the miRNA transcriptome of UC in dogs. Our study aimed to evaluate differential miRNA expression in bladder tissue collected from normal canine urothelium and canine invasive UC (iUC) to elucidate the dysregulated pathways in canine UC. Next-Generation RNA sequencing (RNA-Seq) was performed for dogs with UC (n = 29) and normal canine urothelium (n = 4). Raw RNA data were subjected to normalization, and pairwise comparison was performed using EdgeR with Benjamini-Hochberg FDR multiple testing correction (p < 0.05; >2-fold change) comparing tissue samples of normal urothelium to canine iUC samples. Principal component analysis and hierarchical cluster analysis were performed. MiRNA of FFPE tissue samples of separate iUC (n = 5) and normal urothelium (n = 5) were used to evaluate five miRNAs using RT-qPCR. Pathway analysis was performed utilizing miRWalk, STRING database, and Metascape utilizing KEGG pathways and GO terms databases. Twenty-eight miRNAs were differentially expressed (DE) by RNA-Seq. RT-qPCR confirmed that four miRNAs are significantly downregulated in UC compared to healthy urothelial samples (miR-105a, miR-143, miR-181a, and miR-214). Principal component analysis and hierarchical cluster analysis showed separation between miRNAs in iUC and the control group. The DE miRNAs are most often associated with gene silencing by miRNA, miRNAs in cancer, and miRNAs involved in DNA damage responses. Proteins involved include HRAS, KRAS, ARAF, RAF1, MAPK1, MAP2K1, MAPK3, FGFR3, EGFR, HBEGF, RASSF1, E2F2, E2F3, ERBB2, SRC, MMP1, and UP3KA. The differential expression of miRNAs in canine iUC compared to normal canine urothelial tissue indicates that these markers should be further evaluated for their potential role as diagnostic and therapeutic targets.

Introduction

Urothelial carcinoma (UC) comprises approximately 2% of all naturally occurring canine cancers (1), making it the most common urogenital cancer in dogs (2–4). Most UC are of the high-grade invasive form (iUC), with 78% of dogs with UC having tumors invading the bladder wall and 20% of neighboring organs (1, 2). At the time of diagnosis, approximately 20% of dogs have evidence of metastasis (3), with distant metastasis being associated with a worse prognosis (1, 2, 4). In one study, 58% of dogs were found to have distant metastasis at the time of death (1). When not controlled, the primary tumor may obstruct the urinary tract and cause death (1, 3). Canine iUC has been identified as a naturally occurring model for human muscle-invasive UC (1, 3, 4). As such, the molecular environment has been shown to be comparable between the two species (1, 3–5). There are distinct subtypes of UC identified, basal and luminal (5). Similar to humans, in dogs, the luminal subtype comprises 62% of tumors and commonly has FGFR3, ERBB2, and ERBB3 activating mutations and is generally associated with a better prognosis (5). The basal subtype represents 38% of tumors, is enriched with EGFR and HIF-1 expression, and is often metastatic at presentation (5). The basal form can have squamous and sarcomatoid histologic features and biomarkers associated with epithelial to mesenchymal transition, such as MMP9, SERPINE2, CAV1, KRT14, and RASA3 (5).

Many canine patients with UC present with nonspecific urinary tract clinical signs, such as stranguria, pollakiuria, or hematuria (2). Finding abnormal epithelial cells in the urine along with a mass lesion or thickened bladder wall may increase the suspicion of UC; however, this in itself is not diagnostic, as other conditions may present similarly (2). Additionally, concurrent lower urinary tract disease commonly results in inconclusive results by cytological evaluation of the urine. Histologic evaluation is used for the definitive diagnosis of UC (2). The histologic evaluation of the tumor biopsy requires tissue samples collected by cystoscopy, catheter biopsy, or surgery. Surgery is an invasive procedure. Surgery and cystoscopy require general anesthesia. With cystoscopy, visual examination of the bladder is possible, along with the ability to obtain biopsies for histologic evaluation; however, cystoscopy requires specialized equipment and training that are not universally available. Obtaining an adequate tissue sample size may also be challenging since it is limited by the size of the tools inserted through the urethra. Immunohistochemistry for uroplakin II and III can be used to help distinguish UC from other poorly differentiated carcinomas (1, 2). The most common marker used is uroplakin III; however, while it is a highly specific marker, it has a relatively low sensitivity (6). Currently, CADET^Ⓡ BRAF mutation testing (Antech Diagnosis, Inc.) is commercially available for urine samples and is marketed to aid in UC diagnosis. The mutation of the BRAF gene results in the activation of the MAP kinase pathway, which drives aberrant cell growth, and proliferation (7, 8). Although the BRAF mutation detection test has the advantage of being a urine-based test, it has limitations. First, BRAF mutation is present in cancers other than UC, such as prostatic carcinoma (9), which can also be detected within the urinary tract. Moreover, detecting BRAF mutations does not equate to disease development (9). Most importantly, the absence of BRAF mutations cannot rule out UC as up to 20% of the cases do not have the mutation. Thus, a negative BRAF mutation test cannot rule out UC, nor can a positive test confirm UC specifically (9). While an adjunct test, CADET^Ⓡ BRAF-PLUS, has been marketed as a companion diagnostic test used to potentially identify patients that have UC that do not harbor the BRAF mutation, with reported improved sensitivity up to 95% for UC detection (10), these tests require further study, and there is still a need for identification of additional non-invasive molecular biomarkers for the diagnosis of this disease.

MicroRNAs (miRNAs) have become established biomarkers of cancer (11, 12) because of their stability and abundance in body fluids and reported alterations in miRNA expression in neoplastic diseases (13). MiRNAs are a regulatory class of RNAs first discovered in 1993 (9); they are short (18–24 nucleotides in length) non-coding RNAs responsible for post-transcriptional regulation of gene expression. It is thought that more than half of miRNA genes are located in cancer-associated genomic regions (14). Microarray expression data has shown that aberrant miRNA expression is present in cancer, including bladder, breast, colon, gastric, lung, prostate, and thyroid cancers (14). MiRNA may exhibit oncogenic function in some tumors but act as a tumor suppressor in other types of cancer (14). While aberrant miRNA expression has been documented in UC in dogs (15–17), a complete miRNA expression profile of canine iUC has not yet been identified.

This study evaluated the whole miRNA expression profile (miRNome) of canine UC compared with normal canine urothelial tissue. Our work which includes RNA sequencing (RNA-Seq) to obtain a complete list of miRNAs expressed in normal canine bladder and iUC tumor tissue, will build on earlier reports of miRNA expression in canine UC (16, 17). In addition, we have validated differentially expressed (DE) miRNAs from the RNA-Seq data by quantitative reverse transcription (RT-qPCR) in a different cohort of tissue samples.

Materials and methods

RNA-sequencing

Ethics statement

Owner consent was obtained for patients undergoing cystoscopy for diagnostic evaluation. A tissue sample from this procedure was collected and saved for RNA sequencing analysis. This study was approved by the Purdue Animal Care and Use Committee (Approval Number 1111000169).

Tissue collection

Treatment naïve canine iUC tissue samples (n = 29) were collected by cytoscopic biopsy (5). Patients presenting for cytoscopic biopsy were selected with no exclusion based on age, breed, or sex (5). Diagnosis of iUC was confirmed by histologic evaluation by a board-certified pathologist (JARV); all tumors were high grade with >95% epithelial cellularity (5, 18) (Supplementary Table 1). Normal bladder tissue samples (n = 4) were obtained from dogs that underwent euthanasia for other reasons than bladder-related diseases (5). The urothelial layers were dissected out of the tissues and stored in Trizol^™ reagent (Invitrogen, Carlsbad, CA, USA) at−80°C until processing (5).

RNA sequencing

RNA sequencing procedures were run as previously described (5). Total RNA isolation and purification were achieved using the RNAeasy kit (QIAGEN, Valencia, CA, USA) following the manufacturer's instructions (5). Nanodrop spectrophotometer (NanoDrop Technologies, Wilmington, DE) and Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA) were used to evaluate the RNA quality, and an RNA integrity number of seven or greater was considered acceptable. Ribo-Zero rRNA removal kit (Illumina) was used to remove Ribosomal-RNA, and mRNA libraries were constructed (ScriptSeq v2 RNA-Seq library preparation kit, Epicenter Biotech, Madison, WI). Di- tagged cDNA was amplified by limit-cycle PCR and purified using AMPure XP System (Beck- man Coulter) (19, 20). Illumina HiSeq 4000 platform was used to generate paired-end 150 bp sequence reads, obtaining an average of 50 million reads/sample. Raw reads were cleaned for PCR artifacts, adapter trimmed, and then aligned to the canine genome (CanFam 3.1 reference genome) using COBWeb to obtain expression levels for annotated genes and isoforms (Strand NGS, v3.1, Build 235027, Agilent Technologies, Santa Clara, CA). RNA sequencing was performed by the Biomedical Genomics Core at Nationwide Children's Hospital, Columbus, OH, USA (5).

Data validation

Tissue collection

Archived formalin-fixed, paraffin-embedded (FFPE) tissues were utilized. Histologically confirmed UC (n = 5) and normal urinary bladder mucosa (n = 5) from adult dogs were selected; BRAF status was unknown in these patients. All slides were reviewed by a board-certified pathologist (JARV).

Biomarker selection

Five DE miRNAs by RNA-Seq with a fold change greater than two or with a relevant role in UC, i.e., members of the miR-143/miR-145 cluster (21–25), were chosen for validation by RT-qPCR: miR-32, miR-105a, miR-143, miR-181a-1, and miR-214.

MiRNA extraction

An Illinois sternal bone marrow aspiration needle was utilized to acquire tissue cores. The needles were cleaned before coring the samples by submerging them in bleach, 80% ethanol, and RNase-free water. H&E stained tissue sections were evaluated, then the desired area of the tumor or normal mucosa was marked on the corresponding area of the FFPE tissue block. The Illinois sternal needle was used to take a core punch out of the marked area. The core was then placed in a 2 mL microcentrifuge tube for further processing. RNA was extracted from the tissue cores using the miRNeasy FFPE kit (QIAGEN), following the protocol starting with step four with one section per sample. The control cel-miR-39 (miRNeasy- Serum/Plasma Spike-in, QIAGEN, 3.5 μL/sample) was added to each sample after adding Buffer PKD.

Reverse transcription, quantitative polymerase chain reaction

Reverse transcription to cDNA was performed with miScript Reverse Transcriptase Mix (QIAGEN) using 50 ng of RNA for each reaction. Samples were stored at−20°C until qPCR. The miScript SYBR green kit (QIAGEN) was used to run qPCR on a QuantStudio 3 Real-Time PCR System (ThermoFisher Scientific, Waltham, MA). Primers assays for selected miRNAs, miR-32, miR-105a, miR-143, miR-181a-1, and miR-214, were acquired (QIAGEN). Melting curves were evaluated for each run, ensuring all reactions had a single peak on the melting curve and verifying the amplification of a unique PCR product.

The small nuclear RNU6B and three additional miRNAs, miR-152, miR-872, and miR-1842, were evaluated as endogenous reference genes for the normalization of miRNA expression by RT-qPCR. These markers were chosen based on RNA-Seq evaluation for consistent expression throughout control and iUC samples and literature review (26–29). The exogenous control miR-39 was also used.

Statistical analysis

RNA sequencing

Statistical analysis was evaluated using Strand NGS (20), as previously described (5). Briefly, raw RNA-Seq data were normalized using TMM and DESeq (30), separately, filtered by read metrics, and quantified. For the TMM normalized data, EdgeR (31) was used for pairwise comparison, and DEseq2 was used on DESeq normalized data with Benjamini-Hochberg FDR multiple testing correction (32) (P < 0.05; >2-fold) to compare normal mucosa/urothelial layer samples vs. canine iUC samples (Supplementary Table 2). DE miRNAs identified (Table 1) were pooled and filtered by valid expression values within the standardized read count data, defined as a value above the limit of detection. Only miRNAs with at most one invalid value across all samples were retained. Missing values were imputed by the classification and regression tree (CART) method prior to downstream analysis (33, 34). Agglomerative hierarchical clustering was performed with Pearson correlation distance and average linkage (34). Principal component analysis (PCA) was performed with data scaled to unit variance (34). PCA and hierarchical clustering were performed on R (v4.1.2) with RStudio (v21.09.0).

Reference gene identification

The miRNAs expressed in the RNA-Seq data were evaluated for a standard deviation of <0.05 (n = 40). Within that group, an F-test was performed with the hypothesis that the control and test groups did not show significant variation between the control and the iUC samples. Four miRNAs had a significant (<0.05) p-value and were excluded. A two-tailed Student's t-test was performed on the included 36 miRNAs, and 19 of these miRNAs were found to have a statistically significant difference between the control group and the iUC samples and were also excluded. A literature review was performed, and 15 miRNAs with known association with neoplasia were excluded, and two miRNAs previously used as reference genes, miR-152 (26) and RNUB6 (27–29), were included. NormFinder software was used to assess the stability of the potential reference genes based on quantification cycle (CQ) values (Table 2). This algorithm ranks the candidate reference genes according to their expression stability using a mathematical model of gene expression described by Lindbjerg et al. (35).

RNA-Seq validation

A modified version of the ΔΔCQ method that uses two reference genes was utilized to evaluate the relative expression of each miRNA (36). Fold change was calculated, and a Student's t-test was used to compare gene expression between controls and iUC samples (Table 3).

Pathway analysis

Pathway analysis was performed using the Metascape software (37), miRWalk (38), and STRING database (39) with respect to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (40) and Gene Ontology (GO) terms. All database searches were performed using default filtering parameters. The Metascape analysis program identified statistically enriched terms (GO/KEGG terms), then accumulative hypergeometric p-values and enrichment factors were calculated and used to filter the terms. The significant terms were hierarchically clustered into a tree based on Kappa-statistical similarities among their genes. Then 0.3 kappa score was applied as the threshold to cast the tree into term clusters (37). A Gene Set Enrichment Analysis (GSEA) was performed on the Canis lupus familiaris bladder cancer dataset (Supplementary Table 3).

Results

RNA sequencing

Two hundred and eight miRNAs were found to be expressed in the canine bladder (Supplementary Table 4), with 41 miRNAs minimally expressed (<0.001 from baseline). One hundred and thirteen miRNAs were expressed in both normal control and iUC samples, 44 were found only in normal control urothelial samples, and ten were found only in iUC samples. Twenty-eight miRNAs are found to be DE in iUC samples when compared to control samples (P < 0.05) (Table 1).

Hierarchical clustering based on filtered miRNA subset demonstrates mutual similarity of the normal tissue samples relative to the iUC samples (Figure 1). A heatmap of the Pearson's correlation distance measure compared to miRNA expression highlighted increased expression of miR-214, miR-223, miR-568, miR-24-1, miR-190a, and miR-32 (Figure 2). Principal component analysis based on filtered miRNA subset showed that normal tissue and iUC samples are clearly separated from the first principal component, with low overall variance (14.9%) (Figure 3). No distinction was found when differentiating tumor subtype (luminal or basal), immune score (score 1-5 from lowest non-T-cell inflamed to highest T-cell inflamed) (41), the presence of the BRAF mutation, or tumor grade (tumor grades 2-4 included in this study) (41) (Supplementary Figures 1, 2).

Figure 1

Figure 2

Figure 3

Reverse-transcription quantitative PCR

Reference gene identification

Both miR-1842 and RNU6B were acceptable for use as reference genes for miRNA expression using FFPE tissues, according to NormFinder (Table 2), and used in combination to increase the accuracy of the normalization (42).

RNA-Seq validation

Four of the five miRNAs evaluated (miR-105a, miR-143, miR-181a, and miR-214) were significantly downregulated (P < 0.05, FC < 2) in canine iUC compared with normal control urothelium (Table 3). Although there was a trend toward downregulation, miR-32 was not statistically DE.

Pathway analysis

Metascape software was used to evaluate the 208 expressed miRNA (Supplementary Table 5) and the 28 DE miRNAs (Supplementary Table 6). Metascape uses humans as the target organism for this evaluation due to the similarity of the molecular behavior of canine and human iUC (1, 3–5); an assumption was made that similar pathways would be enriched in dogs as are found in humans. When evaluating all expressed miRNAs, the top three out of 100 associated terms were: gene silencing by miRNA, miRNAs in cancer, and miRNA-mediated gene silencing by inhibiting translation (Figure 4A). The evaluation of only the DE miRNAs showed only the terms: gene silencing by miRNA, miRNAs in cancer, and miRNAs involved in DNA damage response (Figure 4B). The network pathway of these terms shows six distinct networks in all expressed miRNA (Figure 5A) and three distinct networks in the DE miRNAs (Figure 5B). The network nodes were also classified by p-value in the expressed miRNAs (Figure 5C) and the DE miRNAs (Figure 5D).

Figure 4

Figure 5

Twenty-nine pathways were found enriched: EGF, MDM2, ERBB2, MAP2K1, MAPK1, DAPK1, MAPK3, NRAS, HRAS, CDH1, TP53, EGFR, RAF1, RASSF1, DAPK3, E2F1, SRC, KRAS, THBS1, DAPK2, E2F3, ARAF, HBEGF, MMP2, E2F2, FGFR3, MMP1, RPS6KA5, and UPK3A. Network analysis for these pathways was further performed through STRING for the Canis lupus familiaris database (Supplementary Table 7). The DE miRNAs were evaluated for protein pathway associations using the miRWalk database (Figure 6).

Figure 6

Thirteen DE miRNAs (miR-582, miR-991-2, miR-544, miR-568, miR-450a, miR-301a, miR-876, miR-32, miR-7-3, miR-190a, miR-551b, miR-223, and miR-374b) were not previously reported nor had predicted pathways associated with iUC in dogs. Of these, five (miR-582, miR-99a-2, miR-544, miR-568, miR-450a) had no GO term association. Metascape term enrichment analysis found an association with gene silencing by miRNA (GO: 0035195), post-transcriptional gene silencing by RNA (GO: 0035194), and post-transcriptional gene silencing (GO: 0016441) for five of these miRNAs (miR-876, miR-7-3, miR-901, miR-551b, and miR-374b) in the human database. Additionally, Metascape enrichment analysis within the human literature predicted that miR-301a was associated with positive regulation of vascular associated smooth muscle cell proliferation (GO:1904707), regulation of vascular associated smooth muscle cell proliferation (GO:1904705), and positive regulation of smooth muscle cell proliferation (GO:0048661); miR-32 was associated with Wnt signaling pathway (GO:0016055), cell-cell signaling by Wnt (GO:0198738), and cell surface receptor signaling pathway involved in cell-cell signaling (GO:1905114); miR-223 was associated with negative regulation of granulocyte chemotaxis (GO:0071623), negative regulation of neutrophil chemotaxis (GO:0090024), negative regulation of neutrophil migration (GO:1902623).

Discussion

In this study, we have described the miRNome of iUC using RNA-Seq and validated the data by RT-qPCR. The RNA-Seq analysis revealed that 208 miRNAs were expressed in the urothelial tissue with only 28 DE miRNAs (differentially expressed between normal mucosa and iUC). Of the 28 DE miRNAs, 15 were predicted to have associations with pathways known to be altered in canine UC. The validation by RT-qPCR was performed in a different cohort of archived FFPE samples. Four miRNAs were confirmed to be downregulated in iUC compared to normal urothelium from the five miRNAs evaluated. Although a trend toward downregulation was seen, miR-32 did not show significant differential expression, likely due to the limited sample size. Nonetheless, independent validation combined with the pathway analysis shows that the data generated by RNA-Seq is biologically significant. The most robust pathways identified by STRING analysis included the RAS (from “rat sarcoma virus”) oncogenic pathway, the regulation of growth factors including epithelial growth factor receptor (EGFR) and fibroblast growth factor receptors (FGFRs), the danger activated protein kinase (DAPK) pathway, uroplakin protein IIIa (UPK3A), and the E2 factor (E2F) family of transcription factors pathways.

The RAS pathway regulates many cellular processes such as differentiation, proliferation, motility, and transformation (43); It is one of the first described oncogenic pathways and the most altered in human cancers. RAS then activates the mitogen-activated protein kinase (MAPK) pathway in the BRAF/MEK/MAPK signaling axis, resulting in the upregulation of transcription factors that support cell proliferation and survival (44). These pathways are regulated by many post-transcriptional factors, including miRNAs (44). Herein, we identified the downregulation of miR-145, miR-214, and miR-532 in iUC, which are associated with RAS. MiR-145 and miR-214 are experimentally confirmed to target RAS expression; miR-145 is a tumor suppressor found downregulated in several cancers and associated with cancer progression by targeting RAS (45). MiR-145 is also part of the p53 and c-Myc regulatory network, and its downregulation is associated with cancer initiation and development through this pathway. Notably, miR-532 is downregulated in renal cancer through its role in the MAPK pathway (46), and its downregulation results in the proliferation and invasion of bladder cancer in humans through the activation of the Wnt/β-catenin signaling pathway (47). Additionally, our study identified three downregulated miRNAs associated with MAPK regulation by STRING analysis: miR-214, miR-652, and miR-874. MiR-214 induces the MAPK signaling pathway activity increasing cell proliferation and inhibiting apoptosis (48). While miR-652 is downregulated herein and in many cancers; however, it is found upregulated in human bladder cancer (49), adding to the studies reporting different molecular mechanisms of urothelial carcinoma in dogs and humans (50). Similarly, miR-874 was not DE in a study of the human miRNome and was used as a reference gene in that study (29). Comparative studies are warranted to elucidate miRNA regulation of UC in dogs and humans, as many research groups suggest using the dog as a natural animal model of disease (4). Nevertheless, our study suggests that the downregulation of miRNAs associated with the RAS/MAPK pathways may play a role in canine iUC tumorigenesis (23, 51–53).

It is noteworthy that BRAF is part of the RAS pathway as an activator of the MAPK (43). The BRAF mutation, identified in ~80% of canine UC cases (8, 10), results in inappropriate activation of the MAPK pathway with consequent abnormal proliferation and differentiation of cancers, including human and canine UC (43, 54). Interestingly, experimental data in human cancers strongly suggest that dysregulated miRNA expression can activate RAS independently of their oncogenic mutations (54). However, it is not clear if there is a direct association between the altered miRNAs associated with those pathways and the presence of the mutation in this study, as no clear sample separation was evidenced by PCA analysis when analyzing the presence of BRAF mutation in the RNA-Seq data, likely because most iUC samples (25/29) had the BRAF mutation. Nevertheless, additional experimental studies must be conducted to explore those associations.

Growth factors have been extensively researched in cancers, including bladder cancer, because of their role in cancer biology and clinical relevance as prognostic and therapeutic targets, and EGFR is a well-characterized initiator of the RAS signaling pathway upon activation by EGF (30). Both EGFs and FGFRs have reported prognostic and therapeutic value in human UC. Abnormal expression of the EGFR pathway was associated with muscle-invasive bladder cancer and indicated a poor prognosis, and FGFR3 mutations or overexpression indicated a more favorable prognosis (55, 56). Moreover, anti-EGFR treatment has shown promise in treating metastatic urothelial carcinoma in humans (57), as well as antitumor activity in canine iUC (58). Herein we found downregulation of miRNAs predicted to regulate EGF (miR-105a, miR-216a, miR-532, miR-874) and FGFR (miR-329b). These miRNAs have been shown to have numerous associations with various other cancers (59–65). Furthermore, another miRNA (miR-764) predicted to regulate the growth factor ERBB (erythroblastic leukemia viral oncogene homolog), also known as HER, was downregulated in iUC compared to controls. ERBB belongs to a family of proteins containing receptor tyrosine kinases related to the EGFR, and this interaction promotes the development of UC. ERBB functions as a receptor for a tyrosine kinase that can directly affect the RAS signaling pathway (66). ERBB has been associated with other neoplasms; of most note is its association with breast cancer (66). In canine UC, ERBB has been found to be expressed most frequently in the luminal subtype of UC, and the luminal subtype has been associated with a better prognosis than the basal subtype (5, 41). Although, in this study, no distinguishable separation was found between miRNAs in the iUC subtypes on PCA (Supplementary Figure 2). Downregulation of these regulatory miRNAs will result in overexpression of their target genes. Given the role of these growth factors in cancer development and consequent use as prognostic and therapeutic markers, it is plausible to suggest that their correspondent regulatory miRNAs are also potential targets and should be further investigated. Several studies in animal models have shown that downregulated miRNAs that act as tumor-suppressors can be restored and potentially aid in cancer therapy (67).

MiRNA-361 is reported to regulate epithelial to mesenchymal transition through the extracellular signaling-related kinases (ERK), signaled by RAS(68), and acts as a tumor suppressor in prostate cancer by targeting signal transducer and activator of transcription-6 (STAT6), a pro-apoptotic signaler to BCL-xL (69). This study reports two downregulated miRNAs, miR-532 and miR-361, which were predicted to regulate UPK3A expression. Uroplakins are membrane proteins that terminally differentiate to create specialized areas of urothelial membrane called urothelial plaques, which give the bladder its flexibility and the ability to be impenetrable (70). MiRNAs have been previously reported to affect the development of uroplakins by targeting the tight junction-related proteins during the differentiation of urothelial cells (71). Interestingly, overexpression of UPK3A is well-documented in urothelial carcinomas, and UPK3A is widely used as a marker of UC by immunohistochemistry in humans and dogs (6, 72).

Death-associated protein kinases (DAPK) are a family of proteins that regulate apoptosis and non-apoptotic cell death (73). Overexpression of members of the DAPK family (DAPK, Zipk, and DRP-1) is known to cause cell death (73). DAPK has been shown to encode calcium/calmodulin-regulated serine/threonine kinase to induce apoptosis and suppress tumor growth (74). The influx of calcium induces phosphorylation and activation of DAPK, and activation of Beclin-1, inducing phagophore formation (75). This method of tumor suppression has been documented in human bladder cancer (74, 76), and methylation of DAPK has been suggested as a diagnostic marker for bladder cancer in humans (73, 74, 76, 77). In addition, DAPK is one of the pro-apoptotic genes activated by p53 (75). Many DE downregulated miRNAs in this study were found to affect the DAPK family of proteins: miR-145, miR-181a, miR-329b, miR-764, miR-874, miR-652, and miR-143. Interestingly, upon literature review, only miR-103 and miR-107 have been shown to target DAPK, promoting the development of colorectal cancer (75). The findings of this study suggest that there are numerous previously unreported miRNAs that may suppress DAPK and contribute to canine iUC development.

Our study also identified miRNAs involved in regulating the transcription factors E2F2 and E2F3. The E2F family of transcription factors has been extensively studied in the context of cancer development; they play a crucial role in cell cycle progression preferentially through the retinoblastoma protein p53 pathways. It is also well-established that E2F transcription factors are regulated by miRNAs. On the other hand, these factors can regulate miRNA expression through negative feedback in a complex regulatory network (78). Altered expression of these factors has been reported in UC and prostatic carcinoma (79). Three miRNAs with predicted E2F2 association were DE: miR-490 was upregulated and miR-181a and miR-532 were downregulated, whereas two miRNAs, miR-329b and miR-181a, associated with E2F3 were downregulated. Increased expression of E2F3 has been found with advanced tumor stage of human bladder cancer (79). Experimental studies linking these miRNAs with the E2F family are lacking; our study serves as base information for further investigations on these interactions.

Thirteen miRNAs with no reported canine UC pathway association were DE in this study and represent an opportunity for further evaluation. While in the human literature, five of these miRNAs (miR-876, miR-7-3, miR-901, miR-551b, and miR-374b) are involved in gene silencing, known to be a function of miRNAs and a contributing factor to cancer development, three of these have associations with vascular development, cell signaling, and negative chemotaxis regulation (miR-301a, miR-32, and miR-223, respectively) (80–82). MiR-301a is also reported to downregulate p63, resulting in decreased expression of E-cadherin, which is associated with increased potential for prostate cancer reoccurrence (83). MiR-301a is also associated with increased cell proliferation and invasion in colorectal cancer by targeting RUNX family transcription factor 3 (84) and targeting of suppressor of cytokine signaling 6 (85). In pancreatic cancer, miR-301a acts as an activator of the nuclease factor kappa beta (86), associated with numerous cell processes such as cell tumorigenesis, apoptosis, immunity, and inflammation. MiR-32 targets PTEN, a tumor suppressor that inhibits angiogenesis and cell proliferation (87, 88). MiR-223 has numerous previously reported targets and associations with several cancers and inflammatory conditions (89). Our findings indicate that numerous miRNAs and pathways are understudied that may contribute to the development of canine iUC.

A limitation of this study was that different breeds were used for both iUC and control samples. While this may be representative of the general canine population presenting with iUC and the histologic appearance and behavior of iUC in dogs are similar regardless of the breed (4, 5, 41, 90), the effects of the difference in the genetic source are not evaluated in this study. Additionally, although the genes discussed in this study, EGFR, MDM2, ERBB2, MAPK1, HRAS, CDH1, TP53, RASSF1, THBS1, HBEGF, FGFR3, MMP1, UPK3A, and MAP2K1, were also dysregulated in the transcriptome (5), a direct connection of miRNA and mRNA regulation cannot be made as miRNAs are known to have multiple gene targets (13). Data integration and miRNA functional studies in canine iUC are warranted.

In conclusion, this study has characterized the miRNome in canine iUC. We describe 15 DE miRNAs associated with 29 known protein pathways altered in canine iUC, 13 DE miRNAs not associated with known canine iUC pathways, and three associated with pathways not predicted to be altered in iUC. Although the potential of miRNA use as diagnostic markers, prognostic markers, and therapeutic targets has been recognized (91), the paucity of studies on canine diseases and samples indicates a need for additional investigation. This study highlights numerous miRNAs and miRNA targets, supporting previously reported BRAF/MEK/MAPK signaling pathway and newly identified DAPK signaling, UPK3A, and E2F transcription factors that advance the understanding of miRNA interactions and canine iUC development. The discovery of multiple altered miRNAs in different pathways is noteworthy since they can be further investigated as diagnostic markers of iUC, especially for cases that do not harbor common mutations.

Funding

Funding for RNA-Seq evaluation was provided by the National Institutes of Health/National Cancer Institute grant P30CA023168 Supp (DK). The funders had no role in study design, data collection, or analysis. The Purdue University Andrews Fellowship provided salary support for MV. Additional funding was provided by AS' laboratory.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1.

  MutsaersAJWidmetWRKnappDW. Canine transitional cell carcinoma. J Vet Intern Med. (2003) 17:136–44. 10.1111/j.1939-1676.2003.tb02424.x

  • CrossRef
  • Google Scholar

• 2.

  FulkersonCMKnappDW. Management of transitional cell carcinoma of the urinary bladder in dogs: a review. Vet J. (2015) 205:217-25. 10.1016/j.tvjl.2015.01.017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3.

  ShapiroSGRaghunathSWilliamsCMotsinger-ReifAACullenJMLiuTet al. Canine urothelial carcinoma: genomically aberrant and comparatively relevant. Chromosom Res. (2015) 23:311–31. 10.1007/s10577-015-9471-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4.

  KnappDWRamos-VaraJAMooreGEDhawanDBonneyPLYoungKE. Urinary bladder cancer in dogs, a naturally occurring model for cancer biology and drug development. ILAR J. (2014) 55:100–18. 10.1093/ilar/ilu018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5.

  DhawanDHahnNMRamos-VaraJAKnappDW. Naturally-occurring canine invasive urothelial carcinoma harbors luminal and basal transcriptional subtypes found in human muscle invasive bladder cancer. PLoS Genet. (2018) 14:1–15. 10.1371/journal.pgen.1007571

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  Ramos-VaraJAMillerMABoucherMRoudabushAJohnsonGC. Immunohistochemical detection of uroplakin III, cytokeratin 7, and cytokeratin 20 in canine urothelial tumors. Vet Pathol. (2003) 40:55–62. 10.1354/vp.40-1-55

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7.

  RossmanPZabkaTSRupleATuerckDRamos-VaraJALiuLet al. Phase I/II trial of vemurafenib in dogs with naturally occurring, BRAF-mutated urothelial carcinoma. Mol Cancer Ther. (2021) 20:2177–88. 10.1158/1535-7163.MCT-20-0893

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8.

  OstranderEDeckerBParkerHGDhawanDKwonEMK. Homologous mutation to human BRAF V600E is common in naturally occurring canine bladder cancer—evidence for a relevant model system and urine-based diagnostic test. Mol Cancer Res. (2012) 17:1310–4. 10.1201/b11826

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9.

  MochizukiHShapiroSGBreenM. Detection of BRAF mutation in urine DNA as a molecular diagnostic for canine urothelial and prostatic carcinoma. PLoS ONE. (2015) 10:1–12. 10.1371/journal.pone.0144170

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10.

  WileyCWiseCFBreenM. Novel noninvasive diagnostics. Vet Clin North Am Small Anim Pract. (2019) 49:781–91. 10.1016/j.cvsm.2019.05.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11.

  LawrieCHGalSDunlopHMPushkaranBLigginsAPPulfordKet al. Detection of elevated levels of tumour-associated microRNAs in serum of patients with diffuse large B-cell lymphoma. Br J Haematol. (2008) 141:672–5. 10.1111/j.1365-2141.2008.07077.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12.

  MitchellPSParkinRKKrohEMFritzBRWymanSKPogosova-AgadjanyanELet al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A. (2008) 105:10513–8. 10.1073/pnas.0804549105

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13.

  HammondSM. An overview of microRNAs. Adv Drug Deliv Rev. (2015) 87:3–14. 10.1016/j.addr.2015.05.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14.

  ReddyKB. MicroRNA (miRNA) in cancer. Cancer Cell Int. (2015) 15:4–9. 10.1186/s12935-015-0185-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15.

  SahabiKSelvarajahGTAbdullahRCheahYKTanGC. Comparative aspects of microRNA expression in canine and human cancers. J Vet Sci. (2018) 19:162–71. 10.4142/jvs.2018.19.2.162

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16.

  VinallRLKentMSdeVere WhiteRW. Expression of microRNAs in urinary bladder samples obtained from dogs with grossly normal bladders, inflammatory bladder disease, or transitional cell carcinoma. Am J Vet Res. (2012) 73:1626–33. 10.2460/ajvr.73.10.1626

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17.

  KentMSZwingenbergerAWestroppJLBarrettLEDurbin-JohnsonBPGhoshPet al. MicroRNA profiling of dogs with transitional cell carcinoma of the bladder using blood and urine samples. BMC Vet Res. (2017) 13:1–13. 10.1186/s12917-017-1259-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18.

  ChengLMacLennanGTLopez-BeltranA. Histologic grading of urothelial carcinoma: a reappraisal. Hum Pathol. (2012) 43:2097–108. 10.1016/j.humpath.2012.01.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19.

  MortazaviAWilliamsBAMcCueKSchaefferLWoldB. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. (2008) 5:621–8. 10.1038/nmeth.1226

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20.

  WangZGersteinMSnyderM. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. (2009) 10:57–63. 10.1038/nrg2484

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21.

  ItesakoTSekiNYoshinoHChiyomaruTYamasakiTHidakaHet al. The MicroRNA expression signature of bladder cancer by deep sequencing: the functional significance of the miR-195/497 cluster. PLoS ONE. (2014) 9:19–22. 10.1371/journal.pone.0084311

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22.

  SongTXiaWShaoNZhangXWangCWuYet al. Differential miRNA expression profiles in bladder urothelial carcinomas. Asian Pac J Cancer Prev. (2010) 11:905–11.

  • Pubmed Abstract
  • Google Scholar

• 23.

  NoguchiSYasuiYIwasakiJKumazakiMYamadaNNaitoSet al. Replacement treatment with microRNA-143 and−145 induces synergistic inhibition of the growth of human bladder cancer cells by regulating PI3K/Akt and MAPK signaling pathways. Cancer Lett. (2013) 328:353–61. 10.1016/j.canlet.2012.10.017

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24.

  VilladsenSBBramsenJBOstenfeldMSWiklundEDFristrupNGaoSet al. The miR-143/-145 cluster regulates plasminogen activator inhibitor-1 in bladder cancer. Br J Cancer. (2012) 106:366–74. 10.1038/bjc.2011.520

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25.

  AvgerisMMavridisKTokasTStravodimosKFragoulisGEScorilasA. Uncovering the clinical utility of miR-143, miR-145 and miR-224 for predicting the survival of bladder cancer patients following treatment. Carcinogenesis. (2015) 36:528–37. 10.1093/carcin/bgv024

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26.

  LambaVGhodke-PuranikYGuanWLambaJK. Identification of suitable reference genes for hepatic microRNA quantitation. BMC Res Notes. (2014) 7:129. 10.1186/1756-0500-7-129

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27.

  HanYChenJZhaoXLiangCWangYSunLet al. MicroRNA expression signatures of bladder cancer revealed by deep sequencing. PLoS ONE. (2011) 6:1–6. 10.1371/journal.pone.0018286

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28.

  CattoJWFMiahSOwenHCBryantHMyersKDudziecEet al. Distinct MicroRNA alterations characterize high- and low-grade bladder cancer. Cancer Res. (2009) 69:8472–81. 10.1158/0008-5472.CAN-09-0744

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29.

  RatertNMeyerHAJungMMollenkopfHJWagnerIMillerKet al. Reference miRNAs for miRNAome analysis of urothelial carcinomas. PLoS ONE. (2012) 7:1–10. 10.1371/journal.pone.0039309

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30.

  LoveMAndersSHuberW. Analyzing RNA-seq data with DESeq2. Bioconductor. (2017) 2:1–63.

  • Google Scholar

• 31.

  RobinsonMDMcCarthyDJSmythGK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. (2009) 26:139–40. 10.1093/bioinformatics/btp616

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32.

  HaynesW. Encyclopedia of Systems Biology. 1st ed Encyclopedia of Systems Biology.New York, NY:Springer New York : Imprint: Springer (2013).

  • Google Scholar

• 33.

  BuurenS. Flexible Imputation of Missing Data2nd ed. New York, NY: Chapman & Hall/CRC (2012).

  • Google Scholar

• 34.

  JaskowiakPACampelloRJGBCostaIG. Proximity measures for clustering gene expression microarray data: a validation methodology and a comparative analysis. IEEE/ACM Trans Comput Biol Bioinforma. (2013) 10:845–57. 10.1109/TCBB.2013.9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35.

  AndersenCLJensenJLØrntoftTF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. (2004) 64:5245–50. 10.1158/0008-5472.CAN-04-0496

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36.

  RiedelGRüdrichUFekete-DrimuszNMannsMPVondranFWRBockM. An extended ΔCT-method facilitating normalisation with multiple reference genes suited for quantitative RT-PCR analyses of human hepatocyte-like cells. PLoS ONE. (2014) 9:93031. 10.1371/journal.pone.0093031

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37.

  ZhouYZhouBPacheLChangMKhodabakhshiAHTanaseichukOet al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. (2019) 10:1523. 10.1038/s41467-019-09234-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38.

  StichtCDe La TorreCParveenAGretzN. Mirwalk: an online resource for prediction of microrna binding sites. PLoS ONE. (2018) 13:1–6. 10.1371/journal.pone.0206239

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39.

  SzklarczykDGableALNastouKCLyonDKirschRPyysaloSet al. The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. (2021) 49:D605–12. 10.1093/nar/gkaa1074

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40.

  KanehisaMGotoSKEGG. Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30. 10.1093/nar/28.1.27

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41.

  SommerBCDhawanDRupleARamos-VaraJAHahnNMUtturkarSMet al. Basal and luminal molecular subtypes in naturally-occurring canine urothelial carcinoma are associated with tumor immune signatures and dog breed. Bl Cancer. (2021) 7:317–33. 10.3233/BLC-201523

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42.

  YangXHatfieldJTHinzeSJMuXAndersonPJPowellBC. Bone to pick: the importance of evaluating reference genes for RT-qPCR quantification of gene expression in craniosynostosis and bone-related tissues and cells. BMC Res Notes. (2012) 5:222. 10.1186/1756-0500-5-222

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43.

  BoulalasIZaravinosADelakasDSpandidosDA. Mutational analysis of the BRAF gene in transitional cell carcinoma of the bladder. Int J Biol Markers. (2009) 24:17–21. 10.1177/172460080902400103

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44.

  GimpleRCWangX. RAS: striking at the core of the oncogenic circuitry. Front Oncol. (2019) 9:965. 10.3389/fonc.2019.00965

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45.

  YinYYanZPLuNNXuQHeJQianXet al. Downregulation of miR-145 associated with cancer progression and VEGF transcriptional activation by targeting N-RAS and IRS1. Biochim Biophys Acta Gene Regul Mech. (2013) 1829:239–47. 10.1016/j.bbagrm.2012.11.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46.

  ZhaiWMaJZhuRXuCZhangJChenYet al. MiR-532-5p suppresses renal cancer cell proliferation by disrupting the ETS1-mediated positive feedback loop with the KRAS-NAP1L1/P-ERK axis. Br J Cancer. (2018) 119:591–604. 10.1038/s41416-018-0196-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47.

  XieXPanJHanXChenW. Downregulation of microRNA-532-5p promotes the proliferation and invasion of bladder cancer cells through promotion of HMGB3/Wnt/β-catenin signaling. Chem Biol Interact. (2019) 300:73–81. 10.1016/j.cbi.2019.01.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48.

  ZhangHSunPWangYLYuXFTongJJ. MiR-214 promotes proliferation and inhibits apoptosis of oral cancer cells through MAPK/ERK signaling pathway. Eur Rev Med Pharmacol Sci. (2020) 24:3710–6. 10.26355/eurrev_202004_20834

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49.

  StevensMTSaundersBM. Targets and regulation of microRNA-652-3p in homoeostasis and disease. J Mol Med. (2021) 99:755–69. 10.1007/s00109-021-02060-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50.

  CroniseKEDasSHernandezBGReganDPDaileyDDMcGeachanRIet al. Characterizing the molecular and immune landscape of canine bladder cancer. Vet Comp Oncol. (2022) 20:69–81. 10.1111/vco.12740

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51.

  LinTDongWHuangJPanQFanXZhangCet al. MicroRNA-143 as a tumor suppressor for bladder cancer. J Urol. (2009) 181:1372–80. 10.1016/j.juro.2008.10.149

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52.

  MengualLLozanoJJIngelmo-TorresMGazquezCRibalMJAlcarazA. Using microRNA profiling in urine samples to develop a non-invasive test for bladder cancer. Int J Cancer. (2013) 133:2631–41. 10.1002/ijc.28274

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53.

  KentOAChivukulaRRMullendoreMWentzelEAFeldmannGLeeKHet al. Repression of the miR-143/145 cluster by oncogenic Ras initiates a tumor-promoting feed-forward pathway. Genes Dev. (2010) 24:2754–9. 10.1101/gad.1950610

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54.

  MochizukiHKennedyKShapiroSGBreenMB. BRAF mutations in canine cancers. PLoS ONE. (2015) 10:1–9. 10.1371/journal.pone.0129534

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55.

  BlackPCDinneyCPN. Growth factors and receptors as prognostic markers in urothelial carcinoma. Curr Urol Rep. (2008) 9:55–61. 10.1007/s11934-008-0011-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56.

  Le GouxCVacherSSchnitzlerABarry DelongchampsNZerbibMPeyromaureMet al. Assessment of prognostic implication of a panel of oncogenes in bladder cancer and identification of a 3-gene signature associated with recurrence and progression risk in non-muscle-invasive bladder cancer. Sci Rep. (2020) 10:1–11. 10.1038/s41598-020-73642-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57.

  WongYNLitwinSVaughnDCohenSPlimackERLeeJet al. Phase II trial of cetuximab with or without paclitaxel in patients with advanced urothelial tract carcinoma. J Clin Oncol. (2012) 30:3545–51. 10.1200/JCO.2012.41.9572

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58.

  JackSMadhivananKRamadesikanSSubramanianSEdwardsDFElzeyBDet al. A novel, safe, fast and efficient treatment for Her2-positive and negative bladder cancer utilizing an EGF-anthrax toxin chimera. Int J Cancer. (2020) 146:449–60. 10.1002/ijc.32719

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59.

  HoneywellDRCabritaMAZhaoHDimitroulakosJAddisonCL. miR-105 inhibits prostate tumour growth by suppressing CDK6 levels. PLoS ONE. (2013) 8:70515. 10.1371/journal.pone.0070515

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60.

  ShenGRongXZhaoJYangXLiHJiangHet al. MicroRNA-105 suppresses cell proliferation and inhibits PI3K/AKT signaling in human hepatocellular carcinoma. Carcinogenesis. (2014) 35:2748–55. 10.1093/carcin/bgu208

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61.

  YamadaYAraiTKatoMKojimaSSakamotoSKomiyaAet al. Role of pre-miR-532 (miR-532-5p and miR-532-3p) in regulation of gene expression and molecular pathogenesis in renal cell carcinoma. Am J CLin Exp Urol. (2019) 7:11–30.

  • Pubmed Abstract
  • Google Scholar

• 62.

  HuSZhengQWuHWangCLiuTZhouW. miR-532 promoted gastric cancer migration and invasion by targeting NKD1. Life Sci. (2017) 177:15–9. 10.1016/j.lfs.2017.03.019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63.

  WangFChangJTHKaoCJHuangRS. High expression of MIR-532-5p, a tumor suppressor, leads to better prognosis in ovarian cancer both in vivo and in vitro. Mol Cancer Ther. (2016) 15:1123–31. 10.1158/1535-7163.MCT-15-0943

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64.

  GuCCaiJXuZZhouSYeLYanQet al. MiR-532-3p suppresses colorectal cancer progression by disrupting the ETS1/TGM2 axis-mediated Wnt/β-catenin signaling. Cell Death Dis. (2019) 10:739. 10.1038/s41419-019-1962-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65.

  ShiahSGHsiaoJRChangWMChenYWJinYTWongTYet al. Downregulated miR329 and miR410 promote the proliferation and invasion of oral squamous cell carcinoma by targeting Wnt-7b. Cancer Res. (2014) 74:7560–72. 10.1158/0008-5472.CAN-14-0978

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66.

  GlaserAPFantiniDShilatifardASchaefferEMMeeksJJ. The evolving genomic landscape of urothelial carcinoma. Nat Rev Urol. (2017) 14:215–29. 10.1038/nrurol.2017.11

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67.

  KasinskiALSlackFJ. MicroRNAs en route to the clinic: progress in validating and targeting microRNAs for cancer therapy. Nat Rev Cancer. (2011) 11:849–64. 10.1038/nrc3166

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68.

  HuJLiLChenHZhangGLiuHKongRet al. MiR-361-3p regulates ERK1/2-induced EMT via DUSP2 mRNA degradation in pancreatic ductal adenocarcinoma. Cell Death Dis. (2018) 9:807. 10.1038/s41419-018-0839-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69.

  LiuDTaoTXuBChenSLiuCZhangLet al. MiR-361-5p acts as a tumor suppressor in prostate cancer by targeting signal transducer and activator of transcription-6(STAT6). Biochem Biophys Res Commun. (2014) 445:151–6. 10.1016/j.bbrc.2014.01.140

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70.

  WuXRKongXPPellicerAKreibichGSunTT. Uroplakins in urothelial biology, function, and disease. Kidney Int. (2009) 75:1153–65. 10.1038/ki.2009.73

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71.

  ChungPJKChiLMChenCLLiangCLLinCTChangYXet al. MicroRNA-205 targets tight junction-related proteins during urothelial cellular differentiation. Mol Cell Proteomics. (2014) 13:2321–36. 10.1074/mcp.M113.033563

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72.

  KaufmannOVolmerigJDietelM. Uroplakin III is a highly specific and moderately sensitive immunohistochemical marker for primary and metastatic urothelial carcinomas. Am J Clin Pathol. (2000) 113:683–7. 10.1309/PYQC-17CB-063T-Q07J

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73.

  GozuacikDKimchiA. DAPk protein family and cancer. Autophagy. (2006) 2:74–9. 10.4161/auto.2.2.2459

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74.

  DaiLMaCZhangZZengSLiuATangSet al. DAPK promoter methylation and bladder cancer risk: a systematic review and meta-analysis. PLoS ONE. (2016) 11:1–15. 10.1371/journal.pone.0167228

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75.

  SuZYangZXuYChenYYuQ. MicroRNAs in apoptosis, autophagy and necroptosis. Oncotarget. (2015) 6:8474–90. 10.18632/oncotarget.3523

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76.

  VinciSGiannariniGSelliCKuncovaJVillariDValentFet al. Quantitative methylation analysis of BCL2, hTERT, and DAPK promoters in urine sediment for the detection of non-muscle-invasive urothelial carcinoma of the bladder: a prospective, two-center validation study. Urol Oncol Semin Orig Investig. (2011) 29:150–6. 10.1016/j.urolonc.2009.01.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77.

  LeeHHKimSH. Review of non-invasive urinary biomarkers in bladder cancer. Transl Cancer Res. (2020) 9:6554–64. 10.21037/tcr-20-1990

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78.

  EmmrichSPützerBM. Checks and balances: E2F - MicroRNA crosstalk in cancer control. Cell Cycle. (2010) 9:2555–67. 10.4161/cc.9.13.12061

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79.

  OlssonAYFeberAEdwardsSTe PoeleRGiddingsIMersonSet al. Role of E2F3 expression in modulating cellular proliferation rate in human bladder and prostate cancer cells. Oncogene. (2007) 26:1028–37. 10.1038/sj.onc.1209854

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80.

  RoussosETCondeelisJSPatsialouA. Chemotaxis in cancer. Nat Rev Cancer. (2011) 11:573–87. 10.1038/nrc3078

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81.

  GranotZJablonskaJ. Distinct functions of neutrophil in cancer and its regulation. Mediators Inflamm. (2015) 2015:701067. 10.1155/2015/701067

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82.

  LustigBBehrensJ. The Wnt signaling pathway and its role in tumor development. J Cancer Res Clin Oncol. (2003) 129:199–221. 10.1007/s00432-003-0431-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83.

  NamRKBenatarTWallisCJDAmemiyaYYangWGarbensAet al. MiR-301a regulates E-cadherin expression and is predictive of prostate cancer recurrence. Prostate. (2016) 76:869–84. 10.1002/pros.23177

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84.

  WangMLiCYuBSuLLiJJuJet al. Overexpressed miR-301a promotes cell proliferation and invasion by targeting RUNX3 in gastric cancer. J Gastroenterol. (2013) 48:1023–33. 10.1007/s00535-012-0733-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85.

  FangYSunBXiangJChenZ. MiR-301a promotes colorectal cancer cell growth and invasion by directly targeting SOCS6. Cell Physiol Biochem. (2015) 35:227–36. 10.1159/000369690

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86.

  LuZLiYTakwiALiBZhangJConklinDJet al. miR-301a as an NF-κB activator in pancreatic cancer cells. EMBO J. (2011) 30:57–67. 10.1038/emboj.2010.296

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87.

  WuWYangJFengXWangHYeSYangPet al. MicroRNA-32 (miR-32) regulates phosphatase and tensin homologue (PTEN) expression and promotes growth, migration, and invasion in colorectal carcinoma cells. Mol Cancer. (2013) 12:1–11. 10.1186/1476-4598-12-30

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88.

  YanS-YChenM-MLiG-MWangY-QFanJ-G. MiR-32 induces cell proliferation, migration, and invasion in hepatocellular carcinoma by targeting PTEN. Tumor Biol. (2015) 36:4747–55. 10.1007/s13277-015-3124-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89.

  HaneklausMGerlicMO'NeillLAJMastersSL. MiR-223: infection, inflammation and cancer. J Intern Med. (2013) 274:215–26. 10.1111/joim.12099

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90.

  DhawanDRamos-VaraJAStewartJCZhengRKnappDW. Canine invasive transitional cell carcinoma cell lines: in vitro tools to complement a relevant animal model of invasive urinary bladder cancer. Urol Oncol Semin Orig Investig. (2009) 27:284–92. 10.1016/j.urolonc.2008.02.015

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 91.

  SempereLFAzmiASMooreA. microRNA-based diagnostic and therapeutic applications in cancer medicine. Wiley Interdiscip Rev RNA. (2021) 12:1–54. 10.1002/wrna.1662

  • Pubmed Abstract
  • CrossRef
  • Google Scholar