Upregulation of miR181a/miR212 Improves Myogenic Commitment in Murine Fusion-Negative Rhabdomyosarcoma

Fusion-negative rhabdomyosarcoma (FN-RMS) is the most common soft tissue sarcoma of childhood arising from undifferentiated skeletal muscle cells from uncertain origin. Currently used therapies are poorly tumor-specific and fail to tackle the molecular machinery underlying the tumorigenicity and uncontrolled proliferation of FN-RMS. We and other groups recently found that microRNAs (miRNA) network contributes to myogenic epigenetic memory and can influence pluripotent stem cell commitments. Here, we used the previously identified promyogenic miRNAs and tailored it to the murine FN-RMS. Subsequently, we addressed the effects of miRNAs in vivo by performing syngeneic transplant of pre-treated FN-RMS cell line in C57Bl/6 mice. miRNA pre-treatment affects murine FN-RMS cell proliferation in vivo as showed by bioluminescence imaging analysis, resulting in better muscle performances as highlighted by treadmill exhaustion tests. In conclusion, in our study we identified a novel miRNA combination tackling the anti-myogenic features of FN-RMS by reducing proliferation and described novel antitumorigenic therapeutic targets that can be further explored for future pre-clinical applications.

Introduction

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma occurring in the pediatric population, accounting for ∼5% of all pediatric malignancies, and the overall event-free survival is less than 20% in patients with metastatic disease (Chen et al., 2019). RMS can be classified into two subtypes based on the histopathologic characteristics of the tumor, i.e., embryonal RMS (ERMS) and alveolar RMS (ARMS). Approximately 80% of ARMS present the chimeric proteins encoded from the fusion of PAX3 or PAX7 with FOXO1, while 20% are similar to ERMS in terms of clinical outcome (Skapek et al., 2019). ERMS, on the other hand, presents a number of causative mutations including p53, Ras, Rb1, and Pten (Rubin et al., 2011; Seki et al., 2015). It is now accepted that the fusion status better reflects the prognosis and tumor evolution with respect to tumor histology, thus fusion-positive (FP-RMS) and fusion-negative (FN-RMS) have replaced the current classification in order to identify the two RMS subtypes (Hibbitts et al., 2019; Rudzinski et al., 2021).

Given the expression of myogenic markers, it is hypothesized that RMS arises from myoblasts being restrained from fusion with syncytial muscle, thus leading to an invasive neoplastic phenotype (Kashi et al., 2015). Current treatment regimens include multimodality therapy involving surgery, high-dose chemotherapy and radiotherapy to resect the tumor, achieve remission and guarantee local control (Koscielniak et al., 1999). However, these therapeutic approaches fail to restore the myogenic propensity of FN-RMS and the persistence of chemoresistant cancer cells in the tumor are pushing research for the development of new therapies including the use of microRNAs (miRNAs). MiRNAs are small non-coding RNA molecules that bind to partially complementary mRNA sequences, resulting in target degradation or translation inhibition (Lim et al., 2005). Growing evidence is suggesting that miRNA loss-of-function can cause the progression of cancer, thus miRNA replacement therapies have emerged as promising treatment strategies for malignant neoplasms (Bonci et al., 2008; Esposito et al., 2014).

Recently, we have identified a pro-myogenic miRNA cocktail (PMC) by means of transcriptional and miRNA profile comparison between the mesodermal induced pluripotent stem cell-derived progenitors (MiPs) derived from human fibroblast and the skeletal muscle mesoangioblast (MAB)-derived MiPs (Giacomazzi et al., 2017). The identified miRNAs that were upregulated in MAB-derived MiPs were shown to rescue the myogenic intrinsic property of fibroblast-MiPs, thus paving the way for the potential contributions of the PMC in enhancing skeletal muscle differentiation in RMS, whereby the muscle differentiation is impaired.

In this study, we expand on the previously characterized murine model of FN-RMS by performing bulk RNA-seq on the cell lines derived from the model. Next, we hypothesize the use of PMC in murine FN-RMS aimed at restoring the myogenic propensity of the tumor. In order to do so, we first create a selection of PMC (s-PMC) and assess the s-PMC effectiveness on targeting proliferation reduction in murine RMS cells. Subsequently, we demonstrate its efficiency by means of syngeneic transplant of pre-treated murine FN-RMS in C57/Bl6 mice, showing cancer proliferation reduction upon treatment and functional amelioration during physical exercises.

Although deeper investigation is needed to tackle the quiescent tumor cells and metastatic niche and to rescue the tumorigenic phenotype, we identified a novel miRNA combination tackling the tumorigenic features of RMS by reducing proliferation and pinpointed novel antitumorigenic therapeutic targets that can be further explored for future pre-clinical applications.

Materials and Methods

Cell Culture

Human RD18 and murine KMR46 FN-RMS cell lines were cultured in high glucose Dulbecco’s Modified Eagles Medium (DMEM-HG) supplemented with 10% fetal bovine serum (FBS) and 1% Penicillin/Streptomycin (P/S). Cells were maintained under standard incubator conditions humidified atmosphere (95% air, 5% CO2, 37°C) and passaged twice weekly with 0.25% Trypsin-EDTA.

Co-cultures

For co-culture experiments, cells were seeded at 1:20 KMR46/human MABs ratio on collagen-coated vessels in IMDM 15% FBS medium, 1% L-glutamine, 1% non-essential amino acids, 1% sodium pyruvate, 1% ITS, 1:10,000 basic FGF, 0.2% beta-mercaptoethanol. After 24 h, medium was changed and OptiMEM with miRNA-lipofectamine or lipofectamine was added to the wells. A total of 24 h after miRNA addition, cells were differentiated in DMEM 2% horse serum medium for 96–120 h in 5% O₂/5% CO₂ at 37°C.

Reporter Cell Line

Cells were transfected with lentiviral particles expressing EF1a-eGFP-P2A-fLuc and selected with 1 μg/mL puromycin for 2 weeks. GFP was used as a standard gene expression tracer in vitro. Renilla luciferase (Rluc) was used as an optical reporter gene, upon coelenterazine administration, to detect the cell line engraftment in vivo via IVIS bioluminescence imaging (BLI).

MiRNA Transfection

The miR-targeting cocktails were composed as follows: s-PMC, hsa-miR-181a-5p/-212-3p/-146b-5p/-132-5p; s-AMC, anti-hsa-miR-181a-5p/-212-3p/-146b-5p/-132-5p (all MISSION^® microRNA Mimics, Sigma-Aldrich). Before transfection, cells were seeded at 12–15.000 or 50.000 cell number in 24-wells or 6-wells, respectively. One day after plating, cells were transfected using 3 μl/mL Lipofectamine 2000 and 1.6 μl/mL of at 10 μM concentration, resuspended in OptiMEM-medium. After 60–65 h of transfection, cells were trypsinized, counted and either frozen or employed for subsequent experiments.

Immunofluorescence

Following our adjusted protocol (Ronzoni et al., 2021), fixed cells were permeabilized with 1% BSA, 0.2% TritonX-100 in PBS for 5 min and then blocked in 10% donkey serum for 1 h. After 1.5 h of incubation with primary antibodies, diluted in 1% donkey serum, samples were washed three times with PBS, incubated with anti-mouse or anti-rabbit secondary antibody (1:1000) conjugated with 488 or 594 AlexaFluor fluorochromes (Invitrogen Milan, Italy) and nuclei were counterstained with Hoechst 33258 at 1 μg/ml (Sigma, Italy). Here follows the list of primary antibodies and relative dilutions: mouse anti-MyHC (DSHB #MF20), 1:3; rabbit anti-lamin A/C (Epitomics #2966-1), 1:600; mouse anti-Myogenin (Invitrogen #AB_10977211), 1:100; rabbit anti-MyoD (Thermo Fisher Scientific #MA5-12902), 1:100. Imaging was performed at Eclipse Ti microscope (Nikon) by means of Image-Pro Plus 6.0 software (Nikon).

RNA Isolation and Quantitative Real-Time PCR

Total RNA of each sample was extracted using PureLink RNA Mini Kit (Ambion) and DNase I treatment was performed using the DNA-free kit (Ambion). Reverse transcription was performed using the Superscript III First-Strand Synthesis SuperMix (Invitrogen). A 384-well plate was prepared using Platinum SYBR Green qPCRSuperMix-UDG (Invitrogen) as SYBR Green on 1:5 diluted cDNA obtained from 1 μg total RNA (SybrGreen mix, SSIII cDNA production kit and RNA extraction kit from Thermo Fisher Scientific), using Viia7 384-plate reader (Thermo Fisher Scientific; final primer concentration, 100 nM; final volume, 10 μl; PGK, internal reference; thermal profile, 95°C 15 s, 60°C 60 s, 40×). Gene expression was analyzed in triplicate and normalized to the CT mean of Rpl13a and Gapdh as housekeeping genes. For the primer sequences, see Table 1.

TABLE 1

Table 1. List of primer sequences used for gene expression analysis.

RNA-Seq

RNA (>10 μg) samples were verified and processed by the Genomics Core (KU Leuven – UZ Leuven). RNA-sequencing libraries were constructed with the Lexogen library. Samples were indexed with unique adapters and pooled for single read (50 bp) sequencing in Illumina HiSeq2000. RNA-seq reads were aligned with TopHat v2.0.2 to the mouse genome mm10. Transcripts were assessed and quantities were determined by Cufflinks38. Differential expression levels were assessed using DESeq2. Gene Ontology Biological Process (GO:BP) pathways were identified using g:Profiler (Raudvere et al., 2019). Data has been deposited in GEO under accession code GSE175816.

Transwell Migration Assay

To examine the migration potential following miR pretreatment, following our previous titration (Camps et al., 2020) 50,000 cells were seeded in the top chambers of a transwell in serum-free medium (#3422, Corning) as previously described (Balli et al., 2020). As chemoattractant, the lower compartment contained medium supplemented with 10% FBS. Cells were incubated for 24 h and cells that did not migrate through the pores were removed by a cotton swab. Filters were fixed in 4% formaldehyde solution and stained with 0.1% crystal violet. Five pictures per filter were taken and cell number was counted using ImageJ.

Bioluminescence Imaging for in vivo Tumor Engraftment and Growth

Housing and all experimental animal procedures were approved by the Institutional Animal Care and Research Advisory Committee of KU Leuven (ECD #P089/2018). C57/Bl6 mice were injected in the femoral artery with either 1 × 10⁵ KMR46 Fluc+ cells (untreated group) or with 1 × 10⁵ KMR46 Fluc+ cells pretreated for 3 days with s-PMC (treated group). For the injection, cells were suspended in 50 μl saline water. Afterward, the mice were monitored through BLI every other day starting from day 7 after tumor injection for 10 days. For in vivo BLI scans, mice were placed in the flow chamber of IVIS^® Spectrum. Subsequently, 126 mg/kg of D-luciferin was injected subcutaneously (Breuls et al., 2021). Hence, consecutive frames were acquired until the maximum signal intensity was reached. Pulse/sec intensities were calculated by comparing the same ROI for all the animals, after subtracting the background signal coming from not injected mice.

Treadmill Exhaustion Test

A group of five C57/Bl6 mice injected in the femoral artery with 1 × 10⁵ KMR46 Fluc+ cells (untreated group) and five C57/Bl6 mice injected with 1 × 10⁵ KMR46 Fluc+ cells pretreated for 3 days with s-PMC (treated group) underwent functional tests by the treadmill exhaustion test. The test was performed at day 7, 9, 11, and 13 after the beginning of the experiment (day 0). The electric shock frequency and intensity were pulses of 200 ms/pulse of electric current with 2 pulse/s repetition rate (3 Hz) and intensity (1.22 mA), as indicated by Giarratana et al. (2020). The mice were introduced to the treadmill belt and an adaptation time of 5 min was given before the recordings (motor speed set to zero, for 5 min). A training time of 2 min at 4 m/min was set. Later on, the motor speed was set to 7 m/min, with a 1 m/min increase and a constant uphill inclination of 20°, until exhaustion and >10 s stop. The mice were weighted right after every run. Speed (m/min), distance (m), and time (min and s) were registered and used for calculating the work of each run in J. The formula here applied was: Work (J) = body mass (kg) × gravity (9.81 m = s²) × vertical speed (m/s × angle) × time (s).

Statistical Analysis

Sample size for in vitro/in vivo experiments was calculated by means of Sample Size Calculator (parameters: power, 0.80; alpha, 0.05).¹ When applicable, sample size analysis was based on average values obtained from preliminary optimization/validation trials. Two-tailed unpaired Student’s T-tests analysis was performed. All statistical analyses were conducted using Prism v9.1.0 (GraphPad).

Results

RNA-Seq of Murine RMS Model Unveils Key Signature Compared to Satellite Cells

A novel murine model of FN-RMS has been recently described which showed high molecular homology between the genetically engineered murine model (GEMM) and human RMS (Nakahata et al., 2021). From this GEMM, the cell line KMR46 was derived (from here referred to as murine FN-RMS) which was demonstrated to preserve histological resemblance and metastatic potential with respect to its primary GEMM tumors (Nakahata et al., 2021).

In order to understand the key genes that distinguish the tumor from the normally developing skeletal muscle, we performed RNA-seq on the murine FN-RMS cell line and analyzed the differentially expressed genes with respect to the murine satellite cell dataset from Machado et al. (2017) using DESeq2 (Figure 1). We employed the T0 quiescent satellite cells for this comparison as it is largely accepted that a common progenitor of satellite cells could be the cell that transforms and becomes FN-RMS (Tiffin et al., 2003). PC (Figure 1A) and distance matrix (Figure 1B) analyses RNA-seq of all samples revealed similar clustering according to the cell of origin. We identified 15,276 downregulated genes and 5638 upregulated genes (Figure 1C). The discrepancy with respect to the number of genes shown to be downregulated compared to the upregulated ones can be sought in the hypermethylation status that characterizes FN-RMS, as previously described (Mahoney et al., 2012), which leads to the chromatin closure and subsequent reduced expression of several genes. Among the genes we identified, Zbtb20, Hmcn2, Tshz2, and Pax7 were shown to be downregulated in murine FN-RMS compared to satellite cells. Zbtb20 was previously described by Alonso-Martin et al. (2016) to be one of the key genes induced in muscle stem cells and involved in myogenic progression. The downregulation of Pax7 and Hmcn2, on the other hand, points at both the deranged myogenic capacity as well as the poor cross-interactions with other cells, as Hmcn2 is a key gene coding for ECM proteins which allows for efficient Pax7+ cell homing in skeletal muscle (Faralli et al., 2011; Wurmser et al., 2020). In the upregulated gene pool, Pgam1 promotes cancer cell migration (Zhang et al., 2017), while Ppia has been shown to be correlated with poor prognosis in patients with hepatocellular carcinoma (Wang and Yu, 2019). Spp1 (also known as osteopontin) has a key role in the physiology of skeletal muscle growth, and thus its upregulation may hint at its key role in the pathogenesis of FN-RMS (Pagel et al., 2014). Finally, Nefl has a key role in intracellular transport to axons and dendrites, and has been shown to be upregulated in several cancers including head and neck squamous cell carcinoma (Huang et al., 2014).

FIGURE 1

Figure 1. Identification of differentially expressed genes in murine FN-RMS cell lines. (A) Principal component of murine satellite cells (n = 3) and FN-RMS cell lines (n = 2), with PC1 accounting for 98% of the variance observed in the analysis. (B) Sample distance matrix showing closed proximity of the replicates. (C) Volcano plot of the up- and downregulated genes in murine FN-RMS compared to satellite cells. (D) Selection of Gene Ontology Biological Process pathways up- and down-regulated in murine FN-RMS compared to satellite cells.

We then performed pathway analysis using g:Profiler (Figure 1D). Here, we identified pathways regarding metabolic process, cell cycle and chromosome organization as being upregulated in murine FN-RMS compared to satellite cells. Conversely, components of myotube formation such as T-tubule and sarcolemma were pathways found to be downregulated in murine FN-RMS. Intriguingly, cilium was one of the most downregulated pathways in murine FN-RMS, a structure that has been previously described to be present in quiescent satellite cells (Jaafar Marican et al., 2016) and its deregulation in myoblasts was shown to induce the formation of RMS (Fu et al., 2014).

s-PMC Promotes Myogenic Commitment of Murine RMS in vitro

We sought to explore the possible use of previously identified PMC as a novel miRNA cocktail to promote myogenic commitment in murine RMS. We first performed in silico analysis by cross-checking the genes targeted by PMC in MiPs with the differentially expressed genes in murine FN-RMS (Figure 2A). We found an upregulation of Anxa3, Gdf15, and Myl6 and downregulation of Tnnt1, Myh1, Hand1, Notch3, Zeb2, Bmp6, and Myo1d in murine FN-RMS compared to satellite cells. This indeed hints at the hypothesis of the use of PMC to promote myogenic commitment of murine FN-RMS and its subsequent loss of tumorigenic phenotype.

FIGURE 2

Figure 2. Selection of previously identified PMC improves myogenic potential in murine FN-RMS. Heatmaps of a selection of genes previously shown to be differentially expressed in MiPs upon PMC treatment (A) analyzed by RNA-seq and (B) confirmed in murine FN-RMS by qRT-PCR using s-PMC. (C) Myogenin expression in green and (D) quantification upon s-PMC treatment in FN-RMS cell lines compared to vehicle. Nuclei are counterstained in blue with HOECHST. s-PMC, Selected Promyogenic Cocktail. Scale bar, 100 μm. *p < 0.05. (E) Reductionist approach to identify the minimal combination of miRNAs from the s-PMC required to have the promyogenic effect. s-AMC, Selected Antimyogenic Cocktail. **p < 0.01; ***p < 0.001; ****p < 0.0001.

For our experiments, we decided to revisit the previously published PMC as it was specifically tailored for human MiPs. We therefore removed the antagomiR use, which reflected the miRNAs found to be upregulated in f-MiPs. Then, we checked for the murine sequences of the human miRNA mimics used in the PMC and found 4 out of 5 miRNAs (mmu-miR-181a-5p, mmu-miR-212-3p, mmu-miR-146b-5p, and mmu-miR-132-5p). We thus excluded miR-424 since the murine counterpart does not exist, and we named this newly identified PMC as selected PMC (s-PMC) (Supplementary Table 1). We then cross-checked the predicted target transcripts of the miRNAs with the upregulated genes in murine FN-RMS. Intriguingly, we found that 318 targets of miR-181a, 156 targets of miR-212, 138 targets of miR-146b and 12 targets of miR-132 are upregulated in murine FN-RMS compared to normal satellite cells (Table 2 and Supplementary Table 2). In order to address s-PMC effects on murine FN-RMS myogenic committment, we checked the previously identified miRNA targets and their expression after perturbation (Figure 2B). We then assessed myogenin localization, one of the key transcription factors involved in myogenic differentiation, after s-PMC treatment of murine FN-RMS (Figures 2C,D and Supplementary Figure 1). An increased percentage of myogenin was expressed in s-PMC-treated compared to vehicle-treated cells, which corroborated our hypothesis on the increased myogenic commitment of murine FN-RMS after s-PMC use. In order to identify the key miRNAs from the s-PMC responsible for the increased myogenic commitment, we performed in silico analysis using DIANA software to assess the different pathways targeted by the s-PMC (Vlachos et al., 2015; Supplementary Figure 2A). We then used a reductionist approach to identify the minimal selection of s-PMC components required to have a pro-myogenic effect (Figure 2E). To do so, we checked for the molecular expression of key myogenic genes (Acta1, Acta2, Myf5, Tnni2, Myog, Tnni1, Myod, Myh1, and Mrf4) at baseline levels, after s-PMC and after anti-myogenic cocktail. We then removed one and two miRNAs from the s-PMC and checked for the myogenic genes’ expression. Thus, we identified miR-181a/212 as the key combination required to have the pro-myogenic effect observed with the s-PMC, which was further confirmed by the in silico analysis (Supplementary Figure 2B).

TABLE 2

Table 2. List of miRNAs in the s-PMC and the key targeted genes upregulated in FN-RMS.

miR-181a/212 Improve Differentiation While Targeting Proliferation and Migration

In order to identify the effects of miR-181a/212 on murine FN-RMS, we first performed in silico analysis to identify the pathway upregulated in RMS targeted by the miRNA combination (Figure 3A). We first checked the predicted gene targets of miR-181a/212 using miRDB (Chen and Wang, 2020) and identified the pathways upregulated in murine FN-RMS and targeted by miR-181a/212 using g:Profiler. We found migration, neurogenesis, and regulation of cell differentiation to be targeted by the miRNAs (Figure 3B). We therefore aimed at characterizing these effects in vitro. We first checked the regulation of cell differentiation by adding miR-181a/212 and differentiating the murine FN-RMS. We observed an increase of fusion index of MyHC+ myotubes in miRNA-treated compared to vehicle-treated FN-RMS cells (Figure 3C). We then assessed the migration potential of murine FN-RMS after miRNA treatment by means of transwell assay and indeed we observed a reduced migration potential after miRNA treatment (Figure 3D). Finally, as we wanted to assess the specificity of this cocktail in targeting the proliferation of murine FN-RMS, we co-cultured tumor cells with human MABs and quantified the percentage of cell populations after miR-181a/212 using human specific lamin A/C antibody (Figure 3E). Indeed, we observed that healthy muscle cells were not negatively affected by the miRNA in terms of cell number and cell differentiation.

FIGURE 3

Figure 3. miR-181a/212 target differentiation, proliferation, and migration in murine FN-RMS. (A) Overlap and (B) selection of Gene Ontology Biological Process pathways upregulated in murine FN-RMS and targeted by miR-181a/212. (C) miR-181a/212 increases the myogenic commitment of murine FN-RMS as seen by the quantification of fusion index of MyHC+ myotubes. Nuclei are counterstained in blue with HOECHST. Scale bar, 100 μm. (D) Reduced migration capacity of murine FN-RMS upon miR-181a/212 treatment compared to vehicle. *p < 0.05. (E) Qualitative and quantitative effect of miR-181a/212 on murine FN-RMS proliferation in co-culture setting with human MABs compared to vehicle. Human nuclei were identified with anti-lamin A/C Abs in green. Nuclei are counterstained with HOECHST in blue. Scale bar, 100 μm. **p < 0.01; ***p < 0.001.

miR-181a/212 Pre-treatment in FN-RMS Cells Reduces Tumor Size and Improve Running Performance in Tumor Bearing-Mice

Next, we assessed the effect of miR-181a/212 treatment in FN-RMS in terms of tumor progression in vivo and possible implication in impaired physical performance. Thus, we transduced murine FN-RMS cell line with lentiviral vectors carrying fLuc to monitor its activity in vivo (Supplementary Figure 3), and pretreated the cell line with the miRNAs before syngeneic transplants. Twelve C57/Bl6 mice were introduced to the treadmill belt in training sessions to make them familiar with the treadmill and tasks. Then the mice were injected in the femoral artery with saline (control) or 1 × 10⁵ KMR46 Fluc+ cells (untreated group) or 1 × 10⁵ KMR46 Fluc+ cells pretreated for 3 days with miR-181a/212 (treated group). Bioluminescent images revealed that untreated and treated groups showed already different amount of cells at day 7 from transplantation (Figure 4A) while no signal was found in mice injected in the femoral artery with saline (control). From day 11 the amount of cells increased dramatically in the untreated group where up to 6.8 × 10¹⁰ cells were detected at day 13 (Figure 4B). Functional assays were performed at day 7, 9, 11, and 13 after cell injections in trained mice. The mice injected with miRNA-pretreated cells showed an increment in functional outcome and they were able to run for significantly longer period and distance when compared with mice injected with untreated cells (Figure 4C). The better performance was even more pronounced at day 11 and 13 from injections in concomitance with the increased of tumor cell numbers in tumor-bearing mice injected with untreated cells (tumor masses were visible only in this group of mice). Other parameters of treadmill exhaustion test such as work and power showed similar statistically significant differences among the groups and in mice injected with miRNA-pretreated cells at day 11 and 13 those values were similar to sham-operated group (Figure 4C).

FIGURE 4

Figure 4. miR181a/miR212 effects on murine FN-RMS proliferative capacity in vivo. (A) Representative bioluminescent images of C57/Bl6 mice injected in the femoral artery with saline (control) or 1 × 10⁵ KMR46 Fluc+ cells (untreated group) or 1 × 10⁵ KMR46 Fluc+ cells pretreated for 3 days with miR-181a/212 (treated group). Images were taken 7-, 9-, 11-, and 13-days post-injection. (B) Quantification of the BLI signal (n = 3). (C) Functional test by treadmill exhaustion test at 7-, 9-, 11-, and 15-days post-injection. From left to right and top to bottom, graphs of work (J), power (W), distance (m), and time of run (min). (n = 3). **p < 0.01; ***p < 0.001; ****p < 0.0001.

Discussion

The majority of FN-RMSs exhibit embryonal tumor histology differently from FP-RMS characterized by the PAX3-FOXO1 fusion protein and alveolar histology (Tiffin et al., 2003; Linardic, 2008; Davicioni et al., 2009). Emerging literature identified miRNA dysregulation in FN-RMS cells compared to healthy muscle cells involved in the inhibition of adhesion independent growth and migration. For example, miR-22 and miR-378 family members are implicated in controlling proliferation and migration of FN-RMS cells through AKT/mTOR axis (Megiorni et al., 2014; Bersani et al., 2016). Thus, miRNAs have a crucial role in both suppressing the tumorigenicity of RMS and regulating paracrine pathway signaling responsible for tumor invasion and metastatic potential. Indeed, it is generally accepted that exosomal miRNA contents reflect the characteristic miRNA of the cells releasing the extracellular vesicles. However, several reports showed that in case of RMS, a modulation of exosomal miRNA exists compared with tumor cell content. In these circumstances, it could be possible that immune and inflammatory cells, which occur coincidentally at tumor lesions, could contribute to these exosomal miRNAs enrichment (Ghayad et al., 2016).

We and other groups showed that the potential of miR-based orchestration of cell fate is evident along the striated muscle lineages and miRNA mixtures have been described as part of the epigenetic signature retained after cell reprogramming in induced pluripotent stem cells. Thus, since FN-RMS cells are sharing some of gene expression profile with the immature mesodermal progenitors, we sought to test the PMC on FN-RMS cells to rewire the myogenic epigenetic memory similarly observed in treated pluripotent stem cell derivatives. Subsequently, a s-PMC reductionist approach reveals that the exposures to miR-181a and miR-212 only are sufficient to exert the FN-RMS tumor suppressive functions. miR-181a family is highly conserved and regulates many crucial biological processes including proliferation, mitochondrial function, apoptosis, and autophagy (Ouyang et al., 2012). In addition, some authors reported that they are dysregulated also in Alzheimer’s and Parkinson’s diseases (Indrieri et al., 2020). Albeit our s-PMC cocktail was tailored to murine FN-RMS, we confirmed that it is still able to upregulate Tnnt1 and Myl6, and downregulate Bmp6 as previously shown for PMC cocktail by Giacomazzi et al. (2017), while other direct target genes (Gdf15, Anxa3, Hand1, and Myo1d) are slightly affected. Although additional studies are required to pick up the most effective therapeutic miRNA cocktail, our data points on those targeted genes able to slow down RMS progression preserving muscle performance in mice injected with miR-181a and miR-212 pre-treated FN-RMS cells.

Similar to miR-181a, the abnormal expression of miR-212 may cause a series of neurodegenerative diseases, such as Alzheimer’s disease, epilepsy and schizophrenia. More recent studies show that miR-212 is involved in various forms of cancers and can act as a diagnostic or predictive marker of many tumors, including hepatocellular carcinoma ovarian cancer and gastric cancer (Chen et al., 2020). It is not surprising that as in non-small cell lung cancer, ovarian prostate also in FN-RMS act as tumor suppressors and play a substantial role in inhibiting tumorigenesis and FN-RMS progression. This last part is further strengthened by the fact that a remarkable reduction in BLI signal was observed in mice injected with FN-RMS Fluc+ cells pretreated for three days with miR-181a/212. In addition, the evaluation of muscle performance obtained by treadmill exhaustion test at up 13 days post-injection supported the beneficial effect of s-PMC pretreatments in injected mice. We can hypothesize that this beneficial effect would not last for long period since the tumor mass was starting to grow at day 13 even after the miRNA exposure. This is indeed the main limitation of our study that did not explore the effects of multiple miRNA treatment or the stable overexpression of the miRNAs in FN-RMS cells. Nevertheless, we can conclude that as observed in other form of tumors (Chou et al., 2013; Otoukesh et al., 2020) the perturbation of miRNA signature in FN-RMS cells hampered its proliferation and migration ability in vitro and in vivo with consequence in tumor progressions. Both histone modification and DNA hypermethylation can be responsible for the dysregulation of miR-181a/212 in different cancers and their modulation can be used as important biomarker for the diagnosis or prognosis of cancer since they can be easily detected in serum and plasma. The controversial findings regarding their participation in chemoresistance and radioresistance require further studies (Zhang et al., 2010; Taketo et al., 2018). In future experiments it would be necessary to investigate which of the predicted miR-181a and miR-212 gene targets are implicated in the beneficial effect observed in this study and eventually their involvement in the resistance to chemotherapy or radiotherapy will become an important focus in the field of cancer therapy. High-throughput studies are also crucial to identify the precise signature of FN-RMS at single cell level in order to identify dysregulated miRNAs in dormant tumor cells responsible for cancer relapse and drug resistance. It is widely accepted that patient-derived 3D tumoroids better approximate the in vivo physiology and microenvironment, thus providing remarkable tool for miRNA combination screening that allows the assessment of potential therapeutic agents and mechanistic studies.

Data Availability Statement

The datasets presented in this study can be found in online repositories. Data has been deposited in GEO under accession code GSE175816.

Ethics Statement

The animal study was reviewed and approved by the Ethical Committee KU Leuven for Animal Experimentation. Biomedische Wetenschappen, Animalium Gasthuisberg O&N 1 – Herestraat 49 – bus 501, B-3000 Leuven. ECD #P089/2018.

Author Contributions

EP, NG, and MS designed the study and wrote the manuscript. EP, NG, ME, GS, C-CW, VM, and FR performed the experiments. EP and TK performed the bioinformatics analysis. EP and NG interpreted the data. JY and AU supervised the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by INTERREG – Euregio Meuse-Rhine (GYM – Generate your muscle 2020-EMR116), FWO (G0D4517N), C1-KUL3DMUSYC (C14/17/111), and Rondoufonds voor Duchenne Onderzoek (EQQ-FODUCH-O2010). EP was a recipient of FWO (1S90718N).

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.

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.

Acknowledgments

We thank Sylvia Sauvage for her technical support and Vicky Raets and Christina Vochten for their valuable assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.701354/full#supplementary-material

Supplementary Figure 1 | s-PMC and s-AMC effect on Myog and Myod expression. Immunofluorescence of murine FN-RMS after s-PMC or s-AMC treatment compared to vehicle, showing reduced or increased Myogenin (red) and Myod (green) expression after AMC or PMC, respectively. Nuclei are counterstained in blue with HOECHST. Scale bar, 100 μm.

Supplementary Figure 2 | In silico prediction of s-PMC pathways and selection of miR-181a/212. (A) DIANA-miRPath prediction of the GO pathways targeted by the single s-PMC miRNAs, and (B) by miR-181a/212 selection.

Supplementary Figure 3 | Creation of KMR46-EF1a-GFP-fLuc cell line. In vitro validation of the insertion of the cassette by checking (A) GFP expression and (B) fLuc activity upon Renilla addition.

Supplementary Table 1 | Complete list of s-PMC targets predicted by miRDB.

Supplementary Table 2 | Complete list of s-PMC targets predicted by miRDB upregulated in FN-RMS.

Footnotes

1. ^ http://www.stat.ubc.ca/~rollin/stats/ssize/index.html

References

Alonso-Martin, S., Rochat, A., Mademtzoglou, D., Morais, J., de Reynies, A., Aurade, F., et al. (2016). Gene expression profiling of muscle stem cells identifies novel regulators of postnatal myogenesis. Front. Cell Dev. Biol. 4:58. doi: 10.3389/fcell.2016.00058

PubMed Abstract | CrossRef Full Text | Google Scholar

Balli, M., Vitali, F., Janiszewski, A., Caluwé, E., Cortés-Calabuig, A., Carpentier, S., Duelen, R., et al. (2020). Autologous micrograft accelerates endogenous wound healing response through ERK-induced cell migration. Cell Death Differ 27, 1520–1538.

Google Scholar

Bersani, F., Lingua, M. F., Morena, D., Foglizzo, V., Miretti, S., Lanzetti, L., et al. (2016). Deep sequencing reveals a novel miR-22 Regulatory network with therapeutic potential in rhabdomyosarcoma. Cancer Res. 76, 6095–6106. doi: 10.1158/0008-5472.can-16-0709

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonci, D., Coppola, V., Musumeci, M., Addario, A., Giuffrida, R., Memeo, L., et al. (2008). The miR-15a-miR-16-1 cluster controls prostate cancer by targeting multiple oncogenic activities. Nat. Med. 14, 1271–1277. doi: 10.1038/nm.1880

PubMed Abstract | CrossRef Full Text | Google Scholar

Breuls, N., Giarratana, N., Yedigaryan, L., Garrido, G. M., Carai, P., Heymans, S., et al. (2021) Valproic acid stimulates myogenesis in pluripotent stem cell-derived mesodermal progenitors in a NOTCH-dependent manner. Cell Death Dis. 12:677. doi: 10.1038/s41419-021-03936-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Camps, J., Breuls, N., Sifrim, A., Giarratana, N., Corvelyn, M., Danti, L., et al. (2020). Interstitial cell remodeling promotes aberrant adipogenesis in dystrophic muscles. Cell Rep. 31:107597. doi: 10.1016/j.celrep.2020.107597

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Dorado Garcia, H., Scheer, M., and Henssen, A. G. (2019). Current and Future treatment strategies for rhabdomyosarcoma. Front. Oncol. 9:1458. doi: 10.3389/fonc.2019.01458

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Song, J., Bian, H., Yang, X., Xie, X., Zhu, Q., et al. (2020). The functions and targets of miR-212 as a potential biomarker of cancer diagnosis and therapy. J. Cell. Mol. Med. 24, 2392–2401. doi: 10.1111/jcmm.14966

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., and Wang, X. (2020). miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 48, D127–D131.

Google Scholar

Chou, J., Shahi, P., and Werb, Z. (2013). microRNA-mediated regulation of the tumor microenvironment. Cell Cycle 12, 3262–3271. doi: 10.4161/cc.26087

PubMed Abstract | CrossRef Full Text | Google Scholar

Davicioni, E., Anderson, M. J., Finckenstein, F. G., Lynch, J. C., Qualman, S. J., Shimada, H., et al. (2009). Molecular classification of rhabdomyosarcoma–genotypic and phenotypic determinants of diagnosis: a report from the children’s oncology group. Am. J. Pathol. 174, 550–564. doi: 10.2353/ajpath.2009.080631

PubMed Abstract | CrossRef Full Text | Google Scholar

Esposito, C. L., Cerchia, L., Catuogno, S., De Vita, G., Dassie, J. P., Santamaria, G., et al. (2014). Multifunctional aptamer-miRNA conjugates for targeted cancer therapy. Mol. Ther. 22, 1151–1163. doi: 10.1038/mt.2014.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Faralli, H., Martin, E., Core, N., Liu, Q. C., Filippi, P., Dilworth, F. J., et al. (2011). Teashirt-3, a novel regulator of muscle differentiation, associates with BRG1-associated factor 57 (BAF57) to inhibit myogenin gene expression. J. Biol. Chem. 286, 23498–23510. doi: 10.1074/jbc.m110.206003

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, W., Asp, P., Canter, B., and Dynlacht, B. D. (2014). Primary cilia control hedgehog signaling during muscle differentiation and are deregulated in rhabdomyosarcoma. Proc. Natl. Acad. Sci. U. S. A. 111, 9151–9156. doi: 10.1073/pnas.1323265111

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghayad, S. E., Rammal, G., Ghamloush, F., Basma, H., Nasr, R., Diab-Assaf, M., et al. (2016). Exosomes derived from embryonal and alveolar rhabdomyosarcoma carry differential miRNA cargo and promote invasion of recipient fibroblasts. Sci. Rep. 6:37088.

Google Scholar

Giacomazzi, G., Holvoet, B., Trenson, S., Caluwe, E., Kravic, B., Grosemans, H., et al. (2017). MicroRNAs promote skeletal muscle differentiation of mesodermal iPSC-derived progenitors. Nat. Commun. 8:1249.

Google Scholar

Giarratana, N., Conti, F., La Rovere, R., Gijsbers, R., Carai, P., Duelen, R., et al. (2020). MICAL2 is essential for myogenic lineage commitment. Cell Death Dis. 11:654.

Google Scholar

Hibbitts, E., Chi, Y., Hawkins, D. S., Barr, F. G., Bradley, J. A., Dasgupta, R., et al. (2019). Refinement of risk stratification for childhood rhabdomyosarcoma using FOXO1 fusion status in addition to established clinical outcome predictors: a report from the children’s oncology group. Cancer Med. 8, 6437–6448. doi: 10.1002/cam4.2504

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Z., Zhuo, Y., Shen, Z., Wang, Y., Wang, L., Li, H., et al. (2014). The role of NEFL in cell growth and invasion in head and neck squamous cell carcinoma cell lines. J. Oral Pathol. Med. 43, 191–198. doi: 10.1111/jop.12109

PubMed Abstract | CrossRef Full Text | Google Scholar

Indrieri, A., Carrella, S., Carotenuto, P., Banfi, S., and Franco, B. (2020). The pervasive role of the miR-181 family in development, neurodegeneration, and cancer. Int. J. Mol. Sci. 21:2092. doi: 10.3390/ijms21062092

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaafar Marican, N. H., Cruz-Migoni, S. B., and Borycki, A.-G. (2016). Asymmetric distribution of primary cilia allocates satellite cells for self-renewal. Stem Cell Rep. 6, 798–805. doi: 10.1016/j.stemcr.2016.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Kashi, V. P., Hatley, M. E., and Galindo, R. L. (2015). Probing for a deeper understanding of rhabdomyosarcoma: insights from complementary model systems. Nat. Rev. Cancer 15, 426–439. doi: 10.1038/nrc3961

PubMed Abstract | CrossRef Full Text | Google Scholar

Koscielniak, E., Harms, D., Henze, G., Jurgens, H., Gadner, H., Herbst, M., et al. (1999). Results of treatment for soft tissue sarcoma in childhood and adolescence: a final report of the German cooperative soft tissue sarcoma study CWS-86. J. Clin. Oncol. 17, 3706–3719. doi: 10.1200/jco.1999.17.12.3706

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, L. P., Lau, N. C., Garrett-Engele, P., Grimson, A., Schelter, J. M., Castle, J., et al. (2005). Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature 433, 769–773.

Google Scholar

Linardic, C. M. (2008). PAX3-FOXO1 fusion gene in rhabdomyosarcoma. Cancer Lett. 270, 10–18. doi: 10.1016/j.canlet.2008.03.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Machado, L., Esteves de Lima, J., Fabre, O., Proux, C., Legendre, R., Szegedi, A., et al. (2017). In situ fixation redefines quiescence and early activation of skeletal muscle stem cells. Cell Rep. 21, 1982–1993. doi: 10.1016/j.celrep.2017.10.080

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahoney, S. E., Yao, Z., Keyes, C. C., Tapscott, S. J., and Diede, S. J. (2012). Genome-wide DNA methylation studies suggest distinct DNA methylation patterns in pediatric embryonal and alveolar rhabdomyosarcomas. Epigenetics 7, 400–408. doi: 10.4161/epi.19463

PubMed Abstract | CrossRef Full Text | Google Scholar

Megiorni, F., Cialfi, S., McDowell, H. P., Felsani, A., Camero, S., Guffanti, A., et al. (2014). Deep sequencing the microRNA profile in rhabdomyosarcoma reveals down-regulation of miR-378 family members. BMC Cancer 14:880. doi: 10.1186/1471-2407-14-880

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakahata, K., Simons, B. W., Pozzo, E., Shuck, R., Kurenbekova, L., Coarfa, C., et al. (2021). MyoD-Cre driven alterations in K-Ras and p53 lead to a mouse model with histological and molecular characteristics of human rhabdomyosarcoma with direct translational applications. bioRxiv[Preprint]. doi: 10.1101/2021.06.15.448607

CrossRef Full Text | Google Scholar

Otoukesh, B., Abbasi, M., Gorgani, H. O., Farahini, H., Moghtadaei, M., Boddouhi, B., et al. (2020). MicroRNAs signatures, bioinformatics analysis of miRNAs, miRNA mimics and antagonists, and miRNA therapeutics in osteosarcoma. Cancer Cell Int. 20:254.

Google Scholar

Ouyang, Y. B., Lu, Y., Yue, S., and Giffard, R. G. (2012). miR-181 targets multiple Bcl-2 family members and influences apoptosis and mitochondrial function in astrocytes. Mitochondrion 12, 213–219. doi: 10.1016/j.mito.2011.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Pagel, C. N., Wasgewatte Wijesinghe, D. K., Taghavi Esfandouni, N., and Mackie, E. J. (2014). Osteopontin, inflammation and myogenesis: influencing regeneration, fibrosis and size of skeletal muscle. J. Cell Commun. Signal. 8, 95–103. doi: 10.1007/s12079-013-0217-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Raudvere, U., Kolberg, L., Kuzmin, I., Arak, T., Adler, P., Peterson, H., et al. (2019). g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 47, W191–W198.

Google Scholar

Ronzoni, F. L., Giarratana, N., Crippa, S., Quattrocelli, M., Cassano, M., Ceccarelli, G., et al. (2021). Guide cells support muscle regeneration and affect neuro-muscular junction organization. Int. J. Mol. Sci. 22:1939. doi: 10.3390/ijms22041939

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubin, B. P., Nishijo, K., Chen, H. I. H., Yi, X., Scheutze, D. P., Pal, R., et al. (2011). Evidence for an unanticipated relationship between undifferentiated pleomorphic sarcoma and embryonal rhabdomyosarcoma. Cancer Cell 19, 177–191. doi: 10.1016/j.ccr.2010.12.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Rudzinski, E. R., Kelsey, A., Vokuhl, C., Linardic, C. M., Shipley, J., Hettmer, S., et al. (2021). Pathology of childhood rhabdomyosarcoma: a consensus opinion document from the children’s oncology group, European paediatric soft tissue sarcoma study group, and the cooperative Weichteilsarkom Studiengruppe. Pediatr. Blood Cancer 68:e28798.

Google Scholar

Seki, M., Nishimura, R., Yoshida, K., Shimamura, T., Shiraishi, Y., Sato, Y., et al. (2015). Integrated genetic and epigenetic analysis defines novel molecular subgroups in rhabdomyosarcoma. Nat. Commun. 6:7557.

Google Scholar

Skapek, S. X., Ferrari, A., Gupta, A. A., Lupo, P. J., Butler, E., Shipley, J., et al. (2019). Rhabdomyosarcoma. Nat. Rev. Dis. Primers 5:1.

Google Scholar

Taketo, K., Konno, M., Asai, A., Koseki, J., Toratani, M., Satoh, T., et al. (2018). The epitranscriptome m6A writer METTL3 promotes chemo- and radioresistance in pancreatic cancer cells. Int. J. Oncol. 52, 621–629.

Google Scholar

Tiffin, N., Williams, R. D., Shipley, J., and Pritchard-Jones, K. (2003). PAX7 expression in embryonal rhabdomyosarcoma suggests an origin in muscle satellite cells. Br. J. Cancer 89, 327–332. doi: 10.1038/sj.bjc.6601040

PubMed Abstract | CrossRef Full Text | Google Scholar

Vlachos, I. S., Zagganas, K., Paraskevopoulou, M. D., Georgakilas, G., Karagkouni, D., Vergoulis, T., et al. (2015). DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. 43, W460–W466.

Google Scholar

Wang, S., and Yu, J. (2019). Long non-coding RNA transcribed from pseudogene PPIAP43 is associated with radiation sensitivity of small cell lung cancer cells. Oncol. Lett. 18, 4583–4592.

Google Scholar

Wurmser, M., Chaverot, N., Madani, R., Sakai, H., Negroni, E., Demignon, J., et al. (2020). SIX1 and SIX4 homeoproteins regulate PAX7+ progenitor cell properties during fetal epaxial myogenesis. Development 147: dev185975.

Google Scholar

Zhang, D., Jin, N., Sun, W., Li, X., Liu, B., Xie, Z., et al. (2017). Phosphoglycerate mutase 1 promotes cancer cell migration independent of its metabolic activity. Oncogene 36, 2900–2909. doi: 10.1038/onc.2016.446

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, P., Singh, A., Yegnasubramanian, S., Esopi, D., Kombairaju, P., Bodas, M., et al. (2010). Loss of Kelch-like ECH-associated protein 1 function in prostate cancer cells causes chemoresistance and radioresistance and promotes tumor growth. Mol. Cancer Ther. 9, 336–346. doi: 10.1158/1535-7163.mct-09-0589

PubMed Abstract | CrossRef Full Text | Google Scholar