MicroRNA Expression Changes and Integrated Pathways Associated With Poor Outcome in Canine Osteosarcoma
Deanna D. Dailey1,2,3 
Ann M. Hess
Gerrit J. Bouma
Dawn L. Duval- 1Department of Clinical Sciences, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, United States
- 2Flint Animal Cancer Center, Colorado State University, Fort Collins, CO, United States
- 3Cell and Molecular Biology Graduate Program, Colorado State University, Fort Collins, CO, United States
- 4Department of Statistics, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, United States
- 5Department of Biomedical Sciences, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins, CO, United States
- 6Tumor-Host Interactions Research Program, University of Colorado Cancer Center, Anschutz Medical Campus, Aurora, CO, United States
MicroRNAs (miRNA) are small non-coding RNA molecules involved in post-transcriptional gene regulation. Deregulation of miRNA expression occurs in cancer, and miRNA expression profiles have been associated with diagnosis and prognosis in many cancers. Osteosarcoma (OS), an aggressive primary tumor of bone, affects ~10,000 dogs each year. Though survival has improved with the addition of chemotherapy, up to 80% of canine patients will succumb to metastatic disease. Reliable prognostic markers are lacking for this disease. miRNAs are attractive targets of biomarker discovery efforts due to their increased stability in easily obtained body fluids as well as within fixed tissue. Previous studies in our laboratory demonstrated that dysregulation of genes in aggressive canine OS tumors that participate in miRNA regulatory networks is reportedly disrupted in OS or other cancers. We utilized RT-qPCR in a 384-well-plate system to measure the relative expression of 190 miRNAs in 14 canine tumors from two cohorts of dogs with good or poor outcome (disease-free interval >300 or <100 days, respectively). Differential expression analysis in this subset guided the selection of candidate miRNAs in tumors and serum samples from larger groups of dogs. We ultimately identified a tumor-based three-miR Cox proportional hazards regression model and a serum-based two-miR model, each being able to distinguish patients with good and poor prognosis via Kaplan–Meier analysis with log rank test. Additionally, we integrated miRNA and gene expression data to identify potentially important miRNA–mRNA interactions that are disrupted in canine OS. Integrated analyses of miRNAs in the three-miR predictive model and disrupted genes from previous expression studies suggest the contribution of the primary tumor microenvironment to the metastatic phenotype of aggressive tumors.
Introduction
Despite increased survival in osteosarcoma (OS) patients resulting from the addition of chemotherapy to standard treatment protocols, only about one-fourth of canine OS patients will survive longer than a year (1). New treatment strategies are needed to manage this disease and will likely include integration of targeted therapies with standard chemotherapeutics in an individualized medicine setting. To facilitate this effort, biomarkers of disease progression and response to treatment are needed to optimize the stratification of patients into groups most likely to benefit from various treatments and identify targets for development of novel therapeutics.
Previous gene expression studies in our laboratory identified the activation of the Notch signaling pathway in OS but suggested that Notch-independent changes in HES1expression resulted in low HES1 expression in the most aggressive tumors. We also identified upregulation of insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1), an oncofetal protein and known target of the let-7 tumor suppressor family of miRNAs that has been implicated in various cancers (2–4). We hypothesized that the disconnect between the HES1 and Notch pathway activation, as well as the escape of IGF2BP1 from inhibitory mechanisms present in normal adult cells, likely involved the disruption of post-transcriptional regulation by miRNAs.
miRNAs are small non-coding RNAs involved in the regulation of gene expression, providing fine-tuning of multiple cellular processes involved in the development and maintenance of homeostasis. In general, miRNAs suppress the expression of their target genes, and it is estimated that half of mammalian genes are subject to miRNA regulation via 3′ UTR binding sites (5, 6). Since a 2002 report from the Croce laboratory, the involvement of miRNA dysregulation in cancer has been well-established (7). Molecular genomic techniques such as cDNA microarrays and next-generation sequencing have been adapted to facilitate miRNA expression biomarker and novel target discovery efforts (8, 9).
A growing body of literature exploring the significance of miRNA expression changes in OS exists. Several comprehensive reviews have been written to summarize the involvement of miRNAs in OS (10–18). Major findings in OS miRNA studies include suggested or experimentally demonstrated oncogenic- or metastasis-promoting roles for miR-17-92 cluster (19–21), miR-181 family (22–24), miR-27a (23), and miR-21 (25, 26) as well as tumor-suppressive roles for miR-15/16 family members (23) and miR-34 (20, 27, 28). The roles of other miRNAs are less clear, such as the miR-29 family with reports of both elevated and decreased expression in osteosarcoma cell lines and tumors compared to “normal” controls for each sample type (20, 21, 23, 24). Significant bodies of work have explored the association of miRNAs in OS with prognosis. Loss of miRNAs located in the 14q32 locus has been associated with poor patient outcome in both human and canine OS, with the findings in human OS confirmed by multiple groups (29–31). In addition to confirming the oncogene- and tumor-suppressive roles of mir-27a and mir-16, respectively, both in vitro and in vivo, Jones et al. (23) identified tumor-based signatures associated with “osteosarcomagenesis,” metastasis, and response to chemotherapy. Several reports have included functional experiments confirming interactions between miRNAs of interest and genes previously identified as dysregulated in OS, such as loss of 14q32 miRNAs and miR-135 with upregulation of c-MYC, miR-34 with RUNX2, and miR-20a of the miR-17-92 cluster and Fas (19, 28, 31, 32).
miRNAs are attractive molecules for biomarker discovery efforts due primarily to increased stability in biologic fluids and in formalin-fixed tissues compared to other RNA molecules (33–35). These features exemplify the clinical utility of miRNA, particularly in healthcare settings where stringent sample collection and storage requirements necessary for the analysis of mRNA are not always possible. A handful of studies focusing on human OS have identified associations between miRNA expression and outcome, including studies utilizing paraffin-embedded, formalin-fixed tissues, and blood fluids (23, 29, 36). Consequently, we explored the hypothesis that cancer-associated miRNAs would be measurable in tumor and serum and associated with outcome.
Our first objective was to identify candidate biomarker miRNAs differentially expressed in tumors from different outcome groups and in all tumors relative to normal bone. Candidate miRNAs were measured in a larger group of tumors and similarly sized set of serum samples to determine associations between miRNA expression changes and patient outcome. Finally, pathway and miRNA target prediction analyses were used to integrate miRNA and gene expression data to identify potential miRNA–gene regulatory networks important for OS progression.
Materials and Methods
Patient and Tissue/Fluid Selection
Tumors in disease-free interval (DFI) cohorts from dogs with DFI >300 or <100 days treated with limb amputation followed by doxorubicin or platinum-based chemotherapy were collected as previously described (2). Normal bone was obtained from dogs with osteosarcoma from limbs post-amputation and harvested so that “normal” bone included in the study was distant from the tumor site and separated from the tumor by a joint (e.g., a femoral tumor would have matched a distal tibia bone collected). A 1–2-cm section of normal bone was collected for each sample; marrow and medullary fat were removed at collection. Supplementary Table 1 shows patient data for these groups of tumors. Thirty-three additional tumors were selected from the Colorado State University Flint Animal Cancer Center's tissue archive with post-treatment data to document disease progression and matched serum or plasma samples available for miRNA extraction and expression analysis (Supplementary Table 2, COS33). Dogs from both cohorts were confirmed to be free of metastatic lung disease at diagnosis and surgery. Following RNA extraction and quality checks of the serum or plasma samples from the second cohort, 24 of these patients were included in circulating miRNA expression analysis.
Total RNA Isolation, Quantification, and Quality Assessment (Tissues)
RNA was extracted from frozen samples using a freeze fracture device, followed by homogenization and separation of RNA from DNA and protein fractions using TRIzol® Reagent (Life Technologies, Grand Island, NY). The freeze fracture device and the samples were placed in liquid nitrogen to chill for 15–20 min. Approximately 1 cm3 of tumor tissue and up to 4 cm3 of normal metaphyseal bone were used for RNA extraction. Pulverized tissue was transferred into 2 ml/cm3 of tissue of TRIzol in 15-ml conical tubes. The tissue/TRIzol mixture was then homogenized at medium to high speed for 1 min. Homogenized samples were gently shaken, centrifuged for 1 min at 2,000 RPM, and then incubated for 5 min at room temperature. The supernatant was collected into two 1.5-ml tubes and carried forward using the TRIzol reagent manufacturer's protocol for RNA extraction.
After resuspension of the extracted RNA pellet in nuclease-free water, the mirVana™ miRNA extraction kit (Life Technologies, Grand Island, NY) was used for additional RNA purification. RNA was eluted in 50 μl nuclease-free water and treated in 20-μl batches with DNAse (2 μl 10 × DNAse buffer and 2 μl DNAse-I; DNA-free™ kit, Life Technologies) to eliminate genomic DNA contamination. RNA concentration and purity were determined using the NanoDrop 1000 spectrophotometer (NanoDrop Products, Thermo Scientific, Wilmington, DE). The quality of isolated total RNA was determined by RNA integrity number using a Bioanalyser 2100 (Agilent Technologies, Santa Clara, CA) with a RNA 6000 Nano chip. Only samples with RNA integrity number >6 were used. All samples were stored at −80°C.
Total RNA Isolation (Serum)
Archived serum samples stored at −80°C were thawed at room temperature, transferred to RNAse/DNAse Free 2-ml microcentrifuge tubes, and centrifuged for 5 min at 4°C and 16,000 × g. Exactly 200 μl of the supernatant was moved to a fresh 2-ml tube for extraction of RNA using the miRNeasy Serum/Plasma Kit (Qiagen, Valencia, CA) following the manufacturer's directions. Synthetic ce-miR-39 mimic (1.6 × 108 copies) was spiked in to each sample prior to addition of chloroform. Strict preset volumes of reagents and sample RNA were used following the manufacturer's recommendations.
Real-Time Reverse Transcriptase Quantitative PCR
cDNA synthesis of small non-coding RNAs was performed using the miScript Reverse Transcription kit (Qiagen, Valencia, CA) following the manufacturer's instructions. Briefly, reverse transcription (RT) was performed in 20-μl reactions containing 1 μg total RNA in nuclease-free water, 5 × miScript RT Buffer (Mg, dNTPs, and oligo-dTprimers), and 1 μl miScript Reverse Transcriptase Mix [poly(A)polymerase and reverse transcriptase]. Generated cDNAs were stored at −20°C until analysis. Quantitative PCR measurements were performed in 384-well PCR plates in a 6-μl reaction containing 2 × Quantitect SYBR Green master mix (Qiagen, Valencia, CA), 10 μM miRNA-specific forward primer (MWG Biotech), 10 × Universal Reverse Primer (Qiagen, Valencia, CA), 2 ng equivalent cDNA, and nuclease-free water. miRNA-specific primers were designed based on sequences of mature miRNA from MirBase (Supplementary Table 8). Samples were run in duplicate with non-template and reverse transcriptase-free (no RT) controls.
Modifications to this protocol for measurement of serum miRNA were as follows: cDNA synthesis was carried out in 10-μl reactions containing 2 μl 5 × HiSpec Buffer, 1 μl Nucleics mix, 1 μl nuclease-free water, and 5 μl total serum RNA. The serum cDNA was diluted 1:10 in nuclease-free water, and a consistent volume (0.15 μl), rather than a consistent concentration, was included in each 6-μl RT-qPCR reaction.
Data Analysis
For analysis of RT-qPCR data from tumor samples, both GeNorm (37) and NormFinder (38) were used to identify the best candidate reference miRNAs from 10 options, and data was normalized to the geometric mean of miR-30a, miR-27b, and miR-185. The 2−ΔΔCt method was used for differential expression analysis in the initial set of 14 tumors. Statistical analysis of survival data was performed using normalized and transformed expression data from 19 miRNAs in the test set of 33 tumors and 13 miRNAs in 31 serum samples using a combination of Prism and the coxph and survfit functions from the survival package in R.
Statistical Analysis
Associations between miRNA expression levels and DFI were evaluated using Cox proportional hazards linear regression. Multivariable Cox regression was then performed on a subset of candidate miRNAs (p < 0.25 from univariate analysis), utilizing both forward and backward stepwise models based on the Akaike information criterion (AIC). A risk score was calculated for each sample based on the best multivariate model, and Kaplan–Meier method was used to determine median DFI for low- and high-risk groups based on the median risk score. Comparison between groups was made with the log rank analysis, and a p-value of < 0.05 was considered significant.
This analysis pipeline was modified slightly for serum samples. Raw Ct values were first adjusted based on the expression of synthetic cel-miR-39 (33, 39). Then, two miRNAs, miR-16 and miR-21, were selected as reference miRNAs and normalized using a variation of mean centering, termed concordance correlation restricted, as described in Wylie et al. (40). This method was found to be well-suited for biofluid samples.
Results
Differentially Expressed miRNAs in Tumors From Dogs With Poor Response Compared to Those With Good Response
Expression of 188 miRNAs was measured in 14 tumors—seven tumors from dogs with DFI >300 days (good responders) and seven tumors from dogs with DFI <100 days (poor responders)—using RT-qPCR. Four miRNAs were differentially expressed in tumors from poor responders relative to those from good responders using a cutoff of p < 0.05 for significance. Nineteen miRNAs were selected based on p < 0.1, fold change >2.0, or biological interest based on human OS studies for additional expression analyses in a larger set of tumors (Supplementary Table 3, bold).
Cox proportional hazard univariable regression analysis of expression of 19 miRNAs in 33 tumors from patients with DFI ranging from 20 to 937 days identified miRNAs associated with patient outcome (Table 1). The goal of multivariate Cox proportional hazard analysis in this study was to identify the best combination of candidate miRNAs whose expression explained a significant proportion of the variability of patient outcome in this group of tumors and which would be likely to predict outcome in an independent set of canine tumors. Thus, expression values for seven miRNAs with p < 0.25 based on the univariate analysis were included in both forward and backward step-wise multivariate Cox proportional hazard regression analysis. A three-miRNA model was selected as the best model based on AIC, a measurement of model selection that takes into account the goodness-of-fit of the model with penalties for increased complexity (Table 2).
Table 1. Results of univariate Cox proportional hazard regression analysis for expression of miRNAs in canine osteosarcoma tumors (n = 33; disease-free interval range, 20–937 days).
Table 2. Three-miRNA model with lowest Akaike information criterion via both forward and backward step-wise Cox proportional hazard regression (R2 = 0.413, concordance = 0.73).
Three-miRNA Signature for Patient Outcome (DFI)
The Cox proportional hazard multivariate model with three miRNAs—miR-223-3p, miR-130a-3p, and let-7b-5p—was used to calculate the risk scores for each sample. The median risk score was used as a cutoff to discriminate samples considered high or low risk. Kaplan–Meier survival analysis with the log rank test using the three-miRNA model-based risk score distinguished patients with high risk and low risk with respective median DFIs of 123.5 and 392 days (Figure 1A, p = 0.0002, hazard ratio 3.2, 95% confidence interval 2.5–12.9). Relative expression of each miRNA in the signature (Figure 1B) indicated that Let-7b was not significantly elevated in the low-risk group, while miR-103a was significantly elevated in the low-risk group (p = 0.008), and miR-223 was significantly reduced in the low risk group (p = 0.003). Additionally, if samples were separated into cohorts of good and poor responders based on mean DFI, the three-miRNA model signature had an accuracy, based on area under the curve of 0.86 (Figure 1C).
Figure 1. Three-miRNA tumor-based predictive model. (A) Kaplan–Meier survival curve with log rank test (cutoff is median risk score: 0.8897). (B) Relative expression (2−ΔCt) of individual miRNAs in low- and high-risk groups (Mann–Whitney test). (C) Receiver operator characteristic curve for three-miRNA Cox proportional hazard-based risk score dividing outcome groups based on mean disease-free interval.
Pathway Analysis of Dysregulated miRNAs and Genes Suggests Roles for Tumor Microenvironment and IGF2BP1 Regulatory Network in Aggressive OS
We used the mirPath tool from the Diana Tools website (41, 42) with species set to human to identify the top pathways enriched for genes that are targets of the three miRNAs in our Cox proportional hazards model. The top 20 significant pathways using the microT-CDS database and the genes union function include: FoxO signaling, ECM–receptor interaction, signaling pathways regulating pluripotency of stem cells, TGF-beta signaling, cytokine–cytokine receptor interaction, and p53 signaling (Table 3). The let-7 family, being among the earliest miRNAs discovered and more widely studied, shows 682 targets in this tool, while miR-223 and miR-130a list only 367 and 552 genes, respectively. Since loss of let-7b-5p and mir-130a-3p was associated with a shorter DFI, we explored the pathways that they regulate separately from those pathways regulated by mir-223-3p which was elevated in the tumors with higher risk of metastasis. Mir-223-3p is specifically involved in transcriptional misregulation in cancer and cytokine–cytokine receptor interaction. Taken together, this suggests a role for these miRNAs in the regulation of the extracellular environment, immune system, and developmental pathways.
Table 3. Top pathways (p < 0.05) enriched for genes targeted by let-7b-5p, miR-223-3p, and/or miR-130a-3p.
We next used multiMiR, a miRNA–target interaction R package and database out of the Theodorescu lab (43), to identify either experimentally validated or predicted miRNA–mRNA interactions based on data from this study and previous studies in our laboratory (Supplementary Table 6). MultiMiR predicted potential interactions between miR-223, over-expressed in tumors from dogs with shorter DFI, and both dystonin (DST) and catenin (cadherin-associated protein), Alpha 2 (CTNNA2). Both are adhesion proteins interacting with the cytoskeleton, potentially implicating disruption of the tumor microenvironment in the aggressiveness of OS. Interactions between let-7b and six other under-expressed miRNAs and IGF2BP1 confirm that miRNA expression changes likely play a role in the high expression of this gene in tumors from dogs with poor outcome. Relative expression of let-7b and IGF2BP1 in OS tumors with a short disease-free interval via RT-qPCR confirms a statistically significant correlation, suggesting that this interaction occurs in canine OS and may contribute to outcome (Figure 2).
Figure 2. Correlation between low let-7b expression and high expression of IGF2BP1 in eight osteosarcoma tumors as determined by RT-qPCR.
Differentially Expressed miRNAs in OS Tumors Relative to Normal Bone Support Dysregulation of the Notch Pathway in OS
Expression of 188 miRNAs was also measured via RT-qPCR in seven normal bone samples. As has been our experience with gene expression, more differentially expressed miRNAs were identified with higher statistical significance and larger fold changes. Forty differentially expressed miRNAs were identified using cutoffs of p < 0.05 for significance and fold change >2; 21 miRNAs had a lower expression in tumors than normal bone, while 19 miRNAs were over-expressed in tumors (Supplementary Tables 4, 5).
Based on our previous work associating the Notch signaling pathway with OS and again using multiMiR, we sought validated interactions between 21 downregulated miRNAs and 30 upregulated Notch/HES1-associated genes as well as between 19 upregulated miRNAs and 14 downregulated Notch/HES1-associated genes. The pool of Notch/HES1-associated genes was a subset of the genes previously published (44). MultiMiR identified experimental, protein-based evidence for interactions between 21 of 41 miRNAs and 17 of 44 genes or roughly half of the miRNAs and genes entered into the analysis (Figure 3). This data supports the hypothesis that dysregulation of the Notch signaling pathway contributes to the pathogenesis of OS and likely involves disruption of miRNA regulation of Notch pathway-associated genes.
Figure 3. Notch/HES1-associated miRNA–mRNA interactions. Dysregulated genes are shown as ovals or polygons, dysregulated miRNAs are shown in text boxes. In both cases, red indicates expression that is higher in tumors than in normal bone, blue indicates expression that is lower in tumors, and purple indicates that one probe in the Affymetrix array showed NFKB1 as upregulated and another as downregulated. Genes on the left are ligands or inhibitors of Notch; genes on the right are downstream targets of the Notch signaling pathway and/or specifically interact with HES1.
Serum miRNA Changes Associated With OS Patient Outcome
Expression of 13 miRNAs in 31 serum samples from patients with DFI ranging from 20 to 772 days was analyzed using a similar Cox proportional hazard regression pipeline described for tumor miRNA expression data. The 13 miRNAs evaluated comprised a combination of 10 miRNAs selected from our analysis of tumor-derived miRNA expression and three miRNAs commonly highly expressed in human serum samples. Forward and backward stepwise Cox multivariable proportional hazard regression analysis identified a two-miRNA model (miR-23a-3p and miR-30c-5p) with the best fit based on AIC (Table 4, Figures 4A–C). The risk score based on this model separated the samples into groups, with mean DFI of 272 days for the low-risk group and 123.5 days for the high-risk group (p = 0.004, hazard ratio 2.6, 95% confidence interval 1.6–8.5).
Table 4. Two miRNA models after step-wise Cox proportional hazard regression (R2 = 0.278, concordance = 0.69).
Figure 4. Two-miRNA serum-based predictive model. (A) Kaplan–Meier survival curve with log rank test (cutoff is median risk score: 1.0372). (B) Scatter plot of risk scores in two outcome groups based on mean disease-free interval (DFI) for all 33 samples (*p = 0.014, Mann–Whitney test). (C) Receiver operator characteristic curve for serum two-miRNA Cox proportional hazard-based risk score dividing outcome groups based on mean DFI.
Tumor-Based miRNA Signature Compared to Clinical Predictors
One measure of the value of a new prognostic biomarker is its usefulness compared to other predictive markers including clinical parameters (45). For OS, the most consistent clinical indicators of outcome are proximal humerus location, weight, and serum ALP (1, 46–48). We had access to an expanded set of curated, quality-checked clinical data for a subset of our tumors (n = 24) that were included in a large retrospective study by Selmic et al. (46). Multivariate Cox proportional hazard regression of the three miRNA expression-based risk score and other clinical parameters (p < 0.25 on univariate analysis) showed that, when adjusting for these indicators, the miRNA expression based risk score remains a significant predictor of outcome (Table 5). This suggests that incorporation of miRNA expression signatures would improve the estimation of prognosis for canine patients.
Table 5. Results of univariate/multivariate analysis of factors associated with clinical outcome, including a three-miRNA expression-based risk score (tumor-derived miRNA expression).
Discussion
Aberrant miRNA expression patterns have been associated with patient outcome for a variety of human tumors. Combined with their stability in fixed tissues and less invasively obtained body fluids, miRNAs make attractive candidates for biomarker discovery efforts. In this study, we identified miRNA expression signatures from both canine OS tumor and patient serum samples that associated significantly with outcome following surgical amputation of the affected limb and standard-of-care chemotherapy. Pathway and miRNA–gene interaction analyses focused on tumor-derived miRNAs associated with poor outcome, suggesting that the interaction between OS cells and the primary tumor microenvironment may be a major determinant in the ultimate metastatic capabilities of OS tumor cells. Additional miRNA–gene interaction analyses combining expression changes identified in this study with gene expression changes from earlier studies suggest that miRNA dysregulation contributes to both (1) disruption of the Notch pathway in OS compared to normal bone and (2) deregulation of the growth-promoting oncofetal protein IGF2BP1 in the most aggressive OS tumors. Finally, we demonstrated that the tumor-based three-miRNA signature remains an independent predictor of outcome when we control for possible effects of other clinical parameters such as tumor location, patient weight, and age at diagnosis.
Although previous studies have established grading systems for canine OS (49, 50), a limitation of the current study is the lack of grading for the tumors in this data set. Meta-analysis and direct comparisons have shown limited utility for these grading systems in prognosis with simplified high- or low-grade models suggested (48, 51). Variability within tumors as well as the complexity of the criteria in the proposed grading schemes may contribute to high subjectivity. Furthermore, more than 80% of tumors will fall into high-grade histologic categories, within which variable patient outcomes may be achieved. For OS, the most consistent clinical indicators of outcome are proximal humerus location, weight, and serum ALP (1, 46–48). Among these clinical parameters (scoring p < 0.25 on univariate analysis), only location of the tumor in the proximal humerus and the three-miRNA risk score were significant predictors of disease outcome. This suggests that incorporation of miRNA expression signatures would improve the estimation of prognosis for canine patients.
Our first goal of this study was to identify miRNAs associated with progression of OS despite standard-of-care treatment including surgical amputation and doxorubicin and/or platinum-based chemotherapy. We identified a three-miRNA expression signature that separated patients into two distinct outcome groups. Within this signature, elevated expression of miR-223 and decreased expression of let-7b and miR-130a were associated with increased risk and ultimately shorter median DFI. Of these three, miR-223 was the most significantly associated with DFI based on p-value in both the univariate and multivariate regression analyses. Interestingly, the expression of miR-223 is nearly 20 times lower in OS tumors compared to normal bone (Supplementary Table 5), which is consistent with two reports in human OS (21, 23). miRNA expression analyses performed in canine cancer cell lines conducted in our laboratory showed that miR-223 expression is similarly uniformly low across canine osteosarcoma cell lines. We performed pathway and miRNA–gene regulatory analyses to identify pathways potentially affected by expression changes in miR-223. These analyses suggest that the significant increase in miR-223 expression in canine OS tumors may be either originating from or influenced by interactions with the tumor microenvironment.
Pathways enriched for both miR-223 and miR-130a included hematopoietic cell development and osteoclast differentiation. Several lines of evidence support a role for miR-223 as an important regulator of the immune response inhibiting the differentiation of classically activated (M1) macrophages and promoting anti-inflammatory and pro-tumor (M2) polarization (52–54). Notch signaling, a pathway we have found to be significantly dysregulated in aggressive canine OS, is also important for pro-inflammatory M1 polarization (55). Normal differentiation and function of osteoclasts, which are derived from bone marrow monocyte precursors, are also reliant on the expression of miR-223 (56). Given that miR-223 is highly expressed by both M2 macrophages and osteoclasts, it is possible that the increased expression of miR-223 in tumors from poor responders is originating from or induced by an interaction with the increased numbers of these cells in the tumor microenvironment. For example, Yang et al. (57) demonstrated that M2-polarized macrophages can shuttle miR-223 via exosomal transport to breast cancer cells, increasing their invasive ability. In addition, miR-223 may suppress the maturation and immunogenicity of dendritic cells to promote a tolerogenic environment (58, 59).
The role of both osteoclasts and macrophages in OS remains controversial due to a variety of factors, including the potentially different behaviors of these cells depending on the level of differentiation, polarization, and response to external stimuli (60–63). Despite this uncertainty, macrophage-activating agents (promoting pro-inflammatory M1 polarization), such as muramyl tripeptide phosphatidylethanolamine, have consistently shown promise for treatment of OS (64, 65).
In further support of the influence of the tumor microenvironment on miRNA expression changes in OS, miRNA–gene interaction analysis identified potential interactions between miR-223 and adhesion proteins DST and CTNNA2. Both are involved in actin cytoskeletal remodeling, a pathway commonly associated with metastasis (66) and identified as enriched for dysregulated genes in our previous gene expression studies (2). Changes in actin cytoskeletal remodeling are commonly triggered by cell–cell interactions, including those that may occur between tumor cells and supporting stromal cells. The decreased expression of CTNNA2 in tumors from our poor responders supports a pro-metastatic role for miR-223 as CTNNA2 acts as a tumor suppressor in both endometrial and laryngeal carcinomas (67, 68). Additional evidence for an association between miR-223 and metastasis or chemotherapy resistance has been demonstrated in recurrent ovarian tumors, renal cell metastases, and gastric cancer (69–71).
In contrast, recent reports support a potential tumor suppressor role for miR-223 in vitro (72–74). Low miR-223 expression combined with elevated expression of its target gene, epithelial cell transforming sequence 2, in OS tissues is associated with poor outcome (72, 73). It is worth noting that only one of these studies utilized patient tissues, and this included a mix of pediatric and adult tumors (age range 8–66 years) (73). Canine osteosarcoma most closely resembles the pediatric disease, while human adult OS is frequently associated with Paget's disease and may thus involve different underlying molecular mechanisms of progression (1, 75).
Another goal of this study was to integrate miRNA and gene expression data to identify key aberrant pathways contributing to pathogenesis and progression of OS. miRNA–gene expression analysis revealed seven miRNAs with low expression in aggressive tumors predicted or known to target IGF2BP1, an oncogene of interest to our lab. IGF2BP1 has a 3′ UTR that is thousands of kilobases long and with multiple well-conserved binding sites for various miRNAs. This extended 3′ UTR contains multiple polyadenylation sites, with alternate use of polyadenylation sites to produce a shortened 3′ UTR serving as a mechanism by which the gene may avoid miRNA regulation, including at least four sites for the miRNA let-7 (3). We found a statistically significant correlation between low expression of let-7b and increased expression of IGF2BP1 via RT-qPCR in eight OS tumors (Figure 2). Identification of potential miRNA regulators of this protein will facilitate additional functional studies. In addition to IGF2BP1, let-7b targets a variety of oncogenes and has been proposed to act as a tumor suppressor in human osteosarcoma by targeting insulin-like growth factor-1 receptor (76). Despite challenges with effective delivery to target tissues, restoration of tumor-suppressor miRNAs remains a rapidly growing area of research. Studies such as ours may identify new therapeutic miRNAs.
We did not identify any potential interactions between our most dysregulated miRNAs and HES1 nor was Notch signaling identified in our pathway analyses involving miRNAs aberrantly expressed between our DFI cohort tumors. This is consistent with our findings and those of Poos et al. (26) that Notch activation likely contributes to the proliferative response but does not appear to drive metastasis. To further explore the role of miRNAs in Notch activation in OS, we utilized miRNA expression changes identified by comparing tumors to normal bone. We found experimental evidence for interactions between nearly half of the dysregulated miRNAs and one or more Notch/HES1 associated. A handful of these pathways are targetable via small molecule inhibitors including Notch, Hedgehog, HER2/ERRB and PARP. Several of these have been or are under investigation for potential use in the treatment of OS (77, 78). Expression studies like ours might identify biomarkers to help stratify patients for optimal therapeutic benefit or monitor therapeutic response.
The last aim of this study was to identify expression changes of presumed tumor-associated miRNAs in the serum associated with patient outcome. Reliable, repeatable RT-qPCR results for measurement of serum are challenging due to typically low miRNA yield, which inhibits efforts to identify and control for poor-quality samples. This has contributed to inconsistencies between circulating biomarker studies and remains a considerable roadblock to the clinical utility and reliability of such screens (79).
Despite these challenges, we set out to identify a data analysis pipeline utilizing the same relatively affordable SYBR green RT-qPCR platform to measure the relative expression of serum miRNAs. We were able to ultimately identify a two-miRNA signature which successfully stratified patients into distinct outcome groups. The most significantly altered miRNA in this signature was miR-30c, which shows a progressively decreased expression from normal bone to tumors and from tumors from dogs with good outcome to tumors from dogs with poor outcome (Supplementary Tables 3–5). While promising, we acknowledge that, for all of our miRNA signatures, predictive capability in an independent tumor set remains to be established.
Conclusions
In conclusion, we successfully identified miRNA expression changes associated with patient outcome in both OS tumor and patient serum samples. miRNA–gene interactions of the disrupted miRNAs in tumors with genes identified as aberrantly expressed by previous studies (2) can be used to identify targetable pathways disrupted in OS. These studies support the value of miRNA expression studies in biomarker/target discovery efforts for OS.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/ (GSE24251).
Ethics Statement
The animal study was reviewed and approved by Colorado State University Animal Care and Use Committee. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author Contributions
DDa isolated RNA, conducted RT-qPCR and data analysis, prepared the figures, and wrote the manuscript. AH provided assistance in the statistical analysis of data. GB directed the RT-qPCR analysis of miRNAs. DDu directed the study design and sample acquisition and helped with data analysis and writing and editing of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by Morris Animal Foundation D13CA-058 (DDu), Morris Animal Foundation D08CA-053(DDu), and NIH 9T32OD010437-11(DDa).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to acknowledge Dr. Susan Lana and Irene Mok in the Flint Animal Cancer Center Biorepository for collecting and curating osteosarcoma samples.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.637622/full#supplementary-material
References
1. Withrow SJ, Vail DM, Page RL. Withrow and MacEwen's Small Animal Clinical Oncology. Saunders, PA; London: Elsevier (2013).
2. O'Donoghue LE, Ptitsyn AA, Kamstock DA, Siebert J, Thomas RS, Duval LD. Expression profiling in canine osteosarcoma: identification of biomarkers and pathways associated with outcome. BMC Cancer. (2010) 10:506. doi: 10.1186/1471-2407-10-506
3. Mayr C, Bartel PD. Widespread shortening of 3′UTRs by alternative cleavage and polyadenylation activates oncogenes in cancer cells. Cell. (2009) 138:673–84. doi: 10.1016/j.cell.2009.06.016
4. Boyerinas B, Park SM, Shomron N, Hedegaard MM, Vinther J, Andersen JS, et al. Identification of let-7-regulated oncofetal genes. Cancer Res. (2008) 68:2587–91. doi: 10.1158/0008-5472.CAN-08-0264
5. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. (2009) 136:215–33. doi: 10.1016/j.cell.2009.01.002
6. Ebert MS, Sharp PA. Roles for microRNAs in conferring robustness to biological processes. Cell. (2012) 149:515–24. doi: 10.1016/j.cell.2012.04.005
7. Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, et al. Frequent deletions and down-regulation of micro- RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. (2002) 99:15524–9. doi: 10.1073/pnas.242606799
8. Iorio MV, Croce CM. MicroRNA dysregulation in cancer: diagnostics, monitoring and therapeutics. A comprehensive review. EMBO Mol Med. (2012) 4:143–59. doi: 10.1002/emmm.201100209
9. Git A, Dvinge H, Salmon-Divon M, Osborne M, Kutter C, Hadfield J, et al. Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression. RNA. (2010) 16:991–1006. doi: 10.1261/rna.1947110
10. Ell B, Kang Y. MicroRNAs as regulators of bone homeostasis and bone metastasis. Bonekey Rep. (2014) 3:549. doi: 10.1038/bonekey.2014.44
11. Kobayashi E, Hornicek FJ, Duan Z. MicroRNA involvement in osteosarcoma. Sarcoma. (2012) 2012:359739. doi: 10.1155/2012/359739
12. Li B, Ye Z. Epigenetic alterations in osteosarcoma: promising targets. Mol Biol Rep. (2014) 41:3303–15. doi: 10.1007/s11033-014-3193-7
13. Kafchinski LA, Jones KB. MicroRNAs in osteosarcomagenesis. Adv Exp Med Biol. (2014) 804:119–27. doi: 10.1007/978-3-319-04843-7_6
14. Zhang J, Yan YG, Wang C, Zhang SJ, Yu XH, Wang WJ. MicroRNAs in osteosarcoma. Clin Chim Acta. (2015) 444:9–17. doi: 10.1016/j.cca.2015.01.025
15. Ram Kumar RM, Boro A, Fuchs B. Involvement and clinical aspects of microRNA in osteosarcoma. Int J Mol Sci. (2016) 17:877. doi: 10.3390/ijms17060877
16. Botti G, Giordano A, Feroce F, De Chiara AR, Cantile M. Noncoding RNAs as circulating biomarkers in osteosarcoma patients. J Cell Physiol. (209) 234:19249–55. doi: 10.1002/jcp.28744
17. Viera GM, Salamao KB, de Sousa GR, Baroni M, Delsin LEA, Pezuk JA, et al. miRNA signatures in childhood sarcomas and their clinical implications. Clin Transl Oncol. (2019) 21:1583–623. doi: 10.1007/s12094-019-02104-z
18. Guo S-S, Wang Y-J, Zhang G-X, Zhang W-T. Potential diagnostic value of miRNAs in peripheral blood for osteosarcoma: a meta-analysis. J Bone Oncol. (2020) 23:100307. doi: 10.1016/j.jbo.2020.100307
19. Huang G, Nishimoto K, Zhou Z, Hughes D, Kleinerman ES. miR-20a encoded by the miR-17-92 cluster increases the metastatic potential of osteosarcoma cells by regulating Fas expression. Cancer Res. (2012) 72:908–16. doi: 10.1158/0008-5472.CAN-11-1460
20. Baumhoer D, Zillmer S, Unger K, Rosemann M, Atkinson MJ, Irmler M, et al. MicroRNA profiling with correlation to gene expression revealed the oncogenic miR-17-92 cluster to be up-regulated in osteosarcoma. Cancer Genet. (2012) 205:212–9. doi: 10.1016/j.cancergen.2012.03.001
21. Namlos HM, Meza-Zepeda LA, Baroy T, Ostensen IH, Kresse SH, Kuijjer ML, et al. Modulation of the osteosarcoma expression phenotype by microRNAs. PLoS One. (2012) 7:e48086. doi: 10.1371/journal.pone.0048086
22. Hu H, Zhang Y, Cai XH, Huang JF, Cai L. Changes in microRNA expression in the MG-63 osteosarcoma cell line compared with osteoblasts. Oncol Lett. (2012) 4:1037–42. doi: 10.3892/ol.2012.866
23. Jones KB, Salah Z, Del Mare S, Galasso M, Gaudio E, Nuovo GJ, et al. miRNA signatures associate with pathogenesis and progression of osteosarcoma. Cancer Res. (2012) 72:1865–77. doi: 10.1158/0008-5472.CAN-11-2663
24. Andersen GB, Knudsen A, Hager H, Hansen LL, Tost J. miRNA profiling identifies deregulated miRNAs associated with osteosarcoma development and time to metastasis in two large cohorts. Mol Oncol. (2018) 12:114–31. doi: 10.1002/1878-0261.12154
25. Ziyan W, Shuhua Y, Xiufang W, Xiaoyun L. MicroRNA-21 is involved in osteosarcoma cell invasion and migration. Med Oncol. (2011) 28:1469–74. doi: 10.1007/s12032-010-9563-7
26. Poos K, Smida J, Nathrath M, Maugg D, Baumhoer D, Korsching E. How microRNA and transcription factor co-regulatory networks affect osteosarcoma cell proliferation. PLoS Comput Biol. (2013) 9:e1003210. doi: 10.1371/journal.pcbi.1003210
27. Tian Y, Zhang YZ, Chen W. MicroRNA-199a-3p and microRNA-34a regulate apoptosis in human osteosarcoma cells. Biosci Rep. (2014) 34:e00132. doi: 10.1042/BSR20140084
28. van der Deen M, Taipaleenmaki H, Zhang Y, Teplyuk NM, Gupta A, Cinghu S, et al. MicroRNA-34c inversely couples the biological functions of the runt-related transcription factor RUNX2 and the tumor suppressor p53 in osteosarcoma. J Biol Chem. (2013) 288:21307–19. doi: 10.1074/jbc.M112.445890
29. Kelly AD, Haibe-Kains B, Janeway KA, Hill KE, Howe E, Goldsmith J, et al. MicroRNA paraffin-based studies in osteosarcoma reveal reproducible independent prognostic profiles at 14q32. Genome Med. (2013) 5:2. doi: 10.1186/gm406
30. Sarver AL, Thayanithy V, Scott MC, Cleton-Jansen AM, Hogendoorn PC, Modiano JF, et al. MicroRNAs at the human 14q32 locus have prognostic significance in osteosarcoma. Orphanet J Rare Dis. (2013) 8:7. doi: 10.1186/1750-1172-8-7
31. Thayanithy V, Sarver AL, Kartha RV, Li L, Angstadt AY, Breen M, et al. Perturbation of 14q32 miRNAs-cMYC gene network in osteosarcoma. Bone. (2012) 50:171–81. doi: 10.1016/j.bone.2011.10.012
32. Liu Z, Zhang G, Li J, Liu J, Lv P. The tumor-suppressive microRNA-135b targets c-myc in osteoscarcoma. PLoS One. (2014) 9:e102621. doi: 10.1371/journal.pone.0102621
33. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A. (2008) 105:10513–8. doi: 10.1073/pnas.0804549105
34. Doleshal M, Magotra AA, Choudhury B, Cannon BD, Labourier E, Szafranska AE. Evaluation and validation of total RNA extraction methods for microRNA expression analyses in formalin-fixed, paraffin-embedded tissues. J Mol Diagn. (2008) 10:203–11. doi: 10.2353/jmoldx.2008.070153
35. Szafranska AE, Davison TS, Shingara J, Doleshal M, Riggenbach JA, Morrison CD, et al. Accurate molecular characterization of formalin-fixed, paraffin-embedded tissues by microRNA expression profiling. J Mol Diagn. (2008) 10:415–23. doi: 10.2353/jmoldx.2008.080018
36. Zhang C, Yao C, Li H, Wang G, He X. Serum levels of microRNA-133b and microRNA-206 expression predict prognosis in patients with osteosarcoma. Int J Clin Exp Pathol. (2014) 7:4194–203.
37. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. (2002) 3:RESEARCH0034. doi: 10.1186/gb-2002-3-7-research0034
38. Andersen CL, Jensen JL, Orntoft TF. 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. doi: 10.1158/0008-5472.CAN-04-0496
39. Kroh EM, Parkin RK, Mitchell PS, Tewari M. Analysis of circulating microRNA biomarkers in plasma and serum using quantitative reverse transcription-PCR (qRT-PCR). Methods. (2010) 50:298–301. doi: 10.1016/j.ymeth.2010.01.032
40. Wylie D, Shelton J, Choudhary A, Adai AT. A novel mean-centering method for normalizing microRNA expression from high-throughput RT-qPCR data. BMC Res Notes. (2011) 4:555. doi: 10.1186/1756-0500-4-555
41. Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M. DIANA miRPath v.2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. (2012) 40:W498–504. doi: 10.1093/nar/gks494
42. Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. (2015) 43:W460–W6. doi: 10.1093/nar/gkv403
43. Ru Y, Kechris KJ, Tabakoff B, Hoffman P, Radcliffe RA, Bowler R, et al. The multiMiR R package and database: integration of microRNA-target interactions along with their disease and drug associations. Nucleic Acids Res. (2014) 42:e133. doi: 10.1093/nar/gku631
44. Dailey DD, Anfinsen KP, Pfaff LE, Ehrhart E, Charles JB, Bonsdorff TB, et al. HES1, a target of Notch signaling, is elevated in canine osteosarcoma, but reduced in the most aggressive tumors. BMC Vet Res. (2013) 9:130. doi: 10.1186/1746-6148-9-130
45. Zoon CK, Starker EQ, Wilson AM, Emmert-Buck MR, Libutti SK, Tangrea MA. Current molecular diagnostics of breast cancer and the potential incorporation of microRNA. Expert Rev Mol Diagn. (2009) 9:455–67. doi: 10.1586/erm.09.25
46. Selmic LE, Burton JH, Thamm DH, Withrow SJ, Lana SE. Comparison of carboplatin and doxorubicin-based chemotherapy protocols in 470 dogs after amputation for treatment of appendicular osteosarcoma. J Vet Intern Med. (2014) 28:554–63. doi: 10.1111/jvim.12313
47. Ehrhart N, Dernell WS, Hoffmann WE, Weigel RM, Powers BE, Withrow SJ. Prognostic importance of alkaline phosphatase activity in serum from dogs with appendicular osteosarcoma: 75 cases (1990-1996). J Am Vet Med Assoc. (1998) 213:1002–6.
48. Boerman I, Selvarajah GT, Nielen M, Kirpensteijn J. Prognostic factors in canine appendicular osteosarcoma - a meta-analysis. BMC Vet Res. (2012) 8:56. doi: 10.1186/1746-6148-8-56
49. Kirpensteijn J, Kik M, Rutteman GR, Teske E. Prognostic significance of a new histologic grading system for canine osteosarcoma. Vet Pathol. (2002) 39:240–6. doi: 10.1354/vp.39-2-240
50. Loukopoulos P, Robinson WF. Clinicopathological relevance of tumour grading in canine osteosarcoma. J Comp Pathol. (2007) 136:65–73. doi: 10.1016/j.jcpa.2006.11.005
51. Schott CR, Tatiersky LJ, Foster RA, Wood GA. Histologic grade does not predict outcome in dogs with appendicular osteosarcoma receiving the standard of care. Vet Pathol. (2018) 55:202–11. doi: 10.1177/0300985817747329
52. Zhuang G, Meng C, Guo X, Cheruku PS, Shi L, Xu H, et al. A novel regulator of macrophage activation: miR-223 in obesity-associated adipose tissue inflammation. Circulation. (2012) 125:2892–903. doi: 10.1161/CIRCULATIONAHA.111.087817
53. Johnnidis JB, Harris MH, Wheeler RT, Stehling-Sun S, Lam MH, Kirak O, et al. Regulation of progenitor cell proliferation and granulocyte function by microRNA-223. Nature. (2008) 451:1125–9. doi: 10.1038/nature06607
54. Taibi F, Metzinger-Le Meuth V, Massy ZA, Metzinger L. miR-223: an inflammatory oncomiR enters the cardiovascular field. Biochim Biophys Acta. (2014) 1842:1001–9. doi: 10.1016/j.bbadis.2014.03.005
55. Wang YC, He F, Feng F, Liu XW, Dong GY, Qin HY, et al. Notch signaling determines the M1 versus M2 polarization of macrophages in antitumor immune responses. Cancer Res. (2010) 70:4840–9. doi: 10.1158/0008-5472.CAN-10-0269
56. Sugatani T, Hruska KA. MicroRNA-223 is a key factor in osteoclast differentiation. J Cell Biochem. (2007) 101:996–9. doi: 10.1002/jcb.21335
57. Yang M, Chen J, Su F, Yu B, Su F, Lin L, et al. Microvesicles secreted by macrophages shuttle invasion-potentiating microRNAs into breast cancer cells. Mol Cancer. (2011) 10:117. doi: 10.1186/1476-4598-10-117
58. Bros M, Youns M, Kollek V, Buchmüller D, Bollmann F, Seo EJ, et al. Differentially tolerized mouse antigen presenting cells share a common miRNA signature including enhanced mmu-miR-223-3p expression which is sufficient to imprint a protolerogenic state. Front Pharmacol. (2018) 9:915. doi: 10.3389/fphar.2018.00915
59. Tang H-C, Lai Y-Y, Zheng J, Jiang H-Y, Xu G. miR-223-3p Inhibits antigen endocytosis and presentation and promotes the tolerogenic potential of dendritic cells through targeting mannose receptor signaling and Rhob. J Immunol Res. (2020) 2020:1379458. doi: 10.1155/2020/1379458
60. Bao S, Ouyang G, Bai X, Huang Z, Ma C, Liu M, et al. Periostin potently promotes metastatic growth of colon cancer by augmenting cell survival via the Akt/PKB pathway. Cancer Cell. (2004) 5:329–39. doi: 10.1016/S1535-6108(04)00081-9
61. Segaliny AI, Mohamadi A, Dizier B, Lokajczyk A, Brion R, Lanel R, et al. Interleukin-34 promotes tumor progression and metastatic process in osteosarcoma through induction of angiogenesis and macrophage recruitment. Int J Cancer. (2015) 137:73–85. doi: 10.1002/ijc.29376
62. Munoz M, Berger M, Rosso M, Gonzalez-Ortega A, Carranza A, Covenas R. Antitumor activity of neurokinin-1 receptor antagonists in MG-63 human osteosarcoma xenografts. Int J Oncol. (2014) 44:137–46. doi: 10.3892/ijo.2013.2164
63. Buddingh EP, Kuijjer ML, Duim RA, Burger H, Agelopoulos K, Myklebost O, et al. Tumor-infiltrating macrophages are associated with metastasis suppression in high-grade osteosarcoma: a rationale for treatment with macrophage activating agents. Clin Cancer Res. (2011) 17:2110–9. doi: 10.1158/1078-0432.CCR-10-2047
64. Panni RZ, Linehan DC, DeNardo DG. Targeting tumor-infiltrating macrophages to combat cancer. Immunotherapy. (2013) 5:1075–87. doi: 10.2217/imt.13.102
65. Meyers PA, Chou AJ. Muramyl tripeptide-phosphatidyl ethanolamine encapsulated in liposomes (L-MTP-PE) in the treatment of osteosarcoma. Adv Exp Med Biol. (2014) 804:307–21. doi: 10.1007/978-3-319-04843-7_17
66. Aseervatham J. Cytoskeletal remodeling in cancer. Biology. (2020) 9:385. doi: 10.3390/biology9110385
67. Carico E, Atlante M, Giarnieri E, Raffa S, Bucci B, Giovagnoli MR, et al. E-cadherin and alpha-catenin expression in normal, hyperplastic and neoplastic endometrium. Anticancer Res. (2010) 30:4993–7.
68. Fanjul-Fernandez M, Quesada V, Cabanillas R, Cadinanos J, Fontanil T, Obaya A, et al. Cell-cell adhesion genes CTNNA2 and CTNNA3 are tumour suppressors frequently mutated in laryngeal carcinomas. Nat Commun. (2013) 4:2531. doi: 10.1038/ncomms3531
69. Laios A, O'Toole S, Flavin R, Martin C, Kelly L, Ring M, et al. Potential role of miR-9 and miR-223 in recurrent ovarian cancer. Mol Cancer. (2008) 7:35. doi: 10.1186/1476-4598-7-35
70. Wotschofsky Z, Liep J, Meyer HA, Jung M, Wagner I, Disch AC, et al. Identification of metastamirs as metastasis-associated microRNAs in clear cell renal cell carcinomas. Int J Biol Sci. (2012) 8:1363–74. doi: 10.7150/ijbs.5106
71. Zhou X, Jin W, Jia H, Yan J, Zhang G. MiR-223 promotes the cisplatin resistance of human gastric cancer cells via regulating cell cycle by targeting FBXW7. J Exp Clin Cancer Res. (2015) 34:28. doi: 10.1186/s13046-015-0145-6
72. Xu J, Yao Q, Hou Y, Xu M, Liu S, Yang L, et al. MiR-223/Ect2/p21 signaling regulates osteosarcoma cell cycle progression and proliferation. Biomed Pharmacother. (2013) 67:381–6. doi: 10.1016/j.biopha.2013.03.013
73. Zhang H, Yin Z, Ning K, Wang L, Guo R, Ji Z. Prognostic value of microRNA-223/epithelial cell transforming sequence 2 signaling in patients with osteosarcoma. Hum Pathol. (2014) 45:1430–6. doi: 10.1016/j.humpath.2014.02.018
74. Li G, Cai M, Fu D, Chen K, Sun M, Cai Z, et al. Heat shock protein 90B1 plays an oncogenic role and is a target of microRNA-223 in human osteosarcoma. Cell Physiol Biochem. (2012) 30:1481–90. doi: 10.1159/000343336
75. Ottaviani G, Jaffe N. The epidemiology of osteosarcoma. Cancer Treat Res. (2009) 152:3–13. doi: 10.1007/978-1-4419-0284-9_1
76. Zhang K, Wang W, Liu Y, Guo A, Yang D. Let-7b acts as a tumor suppressor in osteosarcoma via targeting IGF1R. Oncol Lett. (2019) 17:1646–54. doi: 10.3892/ol.2018.9793
77. Gill J, Ahluwalia MK, Geller D, Gorlick R. New targets and approaches in osteosarcoma. Pharmacol Ther. (2013) 137:89–99. doi: 10.1016/j.pharmthera.2012.09.003
78. Ma Y, Ren Y, Han EQ, Li H, Chen D, Jacobs JJ, et al. Inhibition of the Wnt-beta-catenin and Notch signaling pathways sensitizes osteosarcoma cells to chemotherapy. Biochem Biophys Res Commun. (2013) 431:274–9. doi: 10.1016/j.bbrc.2012.12.118
Keywords: osteosarcoma, bone cancer, prognosis, miRNA, microRNA, predictive signature, canine (dog)
Citation: Dailey DD, Hess AM, Bouma GJ and Duval DL (2021) MicroRNA Expression Changes and Integrated Pathways Associated With Poor Outcome in Canine Osteosarcoma. Front. Vet. Sci. 8:637622. doi: 10.3389/fvets.2021.637622
Received: 03 December 2020; Accepted: 18 February 2021;
Published: 15 April 2021.
Edited by:
Emanuela Maria Morello, University of Turin, ItalyReviewed by:
Raffaella De Maria, University of Turin, ItalyLola Llobat, Universidad CEU Cardenal Herrera, Spain
Copyright © 2021 Dailey, Hess, Bouma and Duval. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dawn L. Duval, Dawn.Duval@colostate.edu