- Research
- Open access
- Published:
Determination of the miRNA profile of extracellular vesicles from equine mesenchymal stem cells after different treatments
Stem Cell Research & Therapy volume 16, Article number: 162 (2025)
Abstract
Background
Osteoarthritis (OA) is a common and incurable disease in humans and animals. To gain a better understanding of the pathogenesis and identify potential treatments, miRNAs will be extracted and analysed from extracellular vesicles (EVs) of equine adipose derived mesenchymal stem cells (AdMSCs).
Methods
For this purpose we cultivated and pretreated AdMSCs under different conditions: interleukin 1β, shock wave, chondrogenic differentiation, chondrogenic differentiation under hypoxia, or after senescence. After treatment, EVs were harvested from the cell culture supernatants. Next-generation sequencing (NGS) was used to sequence the miRNAs from the EVs.
Results
A total of 89 miRNAs whose expression was significantly altered compared with that of an untreated negative control were identified. On average, 53 miRNAs were upregulated and 6 miRNAs were downregulated. Among others, the miRNAs eca-miR-101, eca-miR-143, eca-miR-145, eca-miR-146a, eca-miR-27a, eca-miR-29b, eca-miR-93, eca-miR-98, and eca-miR-221 were significantly increased after the stimulations, which, as known anti-inflammatory miRNAs, could be candidates for therapeutic use in the treatment of OA.
Conclusion
These results lay the foundation for further research into the significance and efficacy of these miRNAs so that this knowledge can be improved in further experiments and, ideally, translated into therapeutic use.
Background
Osteoarthrosis (OA), also known as osteoarthritis, is a disease that occurs more frequently, particularly in older individuals. This applies both in human medicine [1] and in veterinary medicine [2, 3]. It is considered an incurable, painful degenerative disease of the joint [4]. It is a very heterogeneous disease, accompanied by profound structural changes in the synovial membrane, the cartilage, and the subchondral bone of the joint [5]. For a long period of time, the most important part of treatment was providing patients with adequate analgesic therapy with nonsteroidal anti-inflammatory drugs [6, 7]. These are often combined with intra-articular corticosteroids, hyaluronic acid administration [8], or other adjuvant measures such as a change in diet or weight reduction [7]. Finally, the insertion of a joint prosthesis is often the method of choice to restore the function of the joint with the least possible pain, especially in human medicine [9]. Mesenchymal stem cells constitute another promising treatment method in this context [10], which in today’s context should be called mesenchymal stromal cells or medicinal signaling cells [11]. It has since been recognized that the function of MSCs is not to self-differentiate into appropriate tissues but rather a reaction in their local area, including a healing process through the coordination and stimulation of immune and tissue cells [11, 12].
In recent years, extracellular vesicles (EVs) in particular have emerged from the secretome of MSCs as special mediators [12, 13]. The use of MSC-derived EVs has several advantages over the transplantation of MSCs themselves. In addition to their significantly faster, standardizable production compared with the cultivation of autologous MSCs, they can also maintain consistent quality and reduce the risk of side effects when living cells are used [14]. Extracellular vesicles are vesicles released by nearly all cells. They are coated with a lipid double membrane and have a size in the nanometer range. EVs are released by cells into the extracellular space and can contain various substances, including proteins, lipids, and nucleic acids such as microRNA (miRNA) [15]. Supportive effects for the treatment of OA in general have already been demonstrated for EVs [16] as well as specific associations with miRNAs contained in EVs [17, 18].
This led us to investigate miRNAs in EVs, both from undifferentiated MSCs, as they could be used for the treatment of OA, and from chondrogenically differentiated MSCs, as a possible damage model for the lack of repair in joints affected by OA. Specifically, we want to find out what changes occur in the content of EVs when they are not kept under comparatively normal (cell culture) conditions, but instead are stimulated in different ways. This will involve either a damaging stimulation, such as treatment with interleukin 1β, or a more pro-cartinogenic stimulation in the form of chondrogenic differentiation of MSCs. The comparison of these two approaches should give us an indication of which changes in EVs, here regarding the content of miRNAs, result from these treatments and which potential miRNAs can be found that should be further investigated for a potential treatment of OA. We are using the equine model for this, both because we see a benefit for the care of this species, and because it is a suitable model for human OA.
Methods
AdMSC collection
The donors required for the experiments were already available in the cryobank at the required passages. Specifically, they came from horses in which adipose tissue was taken during a noninflammatory surgical procedure or directly after slaughter at the abattoir. The sampling procedure has already been published [19]. The samples were collected from three horses, a stallion, a gelding, and a mare. The average age was 13.3 (3–19) years. The authorization for this was obtained from the local authorities and is registered under the number V 54 − 19 c 20 15 h 02 GI 18/1 kTV 1/2018.
AdMSC cultivation and differentiation
Equine adipose-derived MSCs were incubated at 37 °C and 5% CO2 in a humidified atmosphere. For this purpose, cryopreserved portions of AdMSCs at passage two were rapidly thawed in a 37 °C water bath and then incubated with standard medium consisting of DMEM low glucose (DMEM-LG, order no. 11564446, Thermo Fisher Scientific, Germany), 10% fetal calf serum (FCS, batch no 201004, PAN, Germany) and 1% antibiotics (Penicillin-Streptomycin, order no. 11548876, Fisher Scientific, Germany) by centrifugation at 300 × g for 5 min, after which the cryomedium was discarded. The cells were then placed in 75 cm² cell culture flasks (order no. 734–2066, VWR, Germany) at a concentration of 0.5 million cells per bottle with 12 ml of standard medium. When the cells reached approximately 80% confluence, they were passaged again for the experiment, so that they were at passage 3.
The detection of mesenchymal stem cells is routinely performed in our laboratory following the International Society for Cellular Therapy (ISCT) standards [20] and is not described in detail here. Briefly, cells obtained from various donors, were examined for the following characteristics of stem cells: plastic adherence, differentiation into three typical lineages (osteogenic, adipogenic, and chondrogenic), and flow cytometric validation of stem cell-specific markers, such as CD90 (positive marker, monoclonal antibody clone 5E10), CD44 (positive marker, monoclonal antibody clone IM7), MHC II (negative marker, monoclonal antibody clone CVS20) and CD 45 (negative marker, monoclonal antibody clone UCHL1) [21, 22, 23].
The differentiations into the osteogenic and adipogenic lineages were carried out in 24-well plates (order no. 734–2325, VWR, Germany) for 2 weeks in the respective osteogenic media (DMEM-LG, 5% FCS, 1% antibiotics, 0.1 µM dexamethasone, 10 mM beta-glycerolphosphate, 60 µM ascorbic acid) or adipogenic media (DMEM low glucose, 5% FCS, 1% antibiotics, 1% insulin-transferrin-selenite solution (ITS, order no. ITS-H, Capricorn, Germany), 1 µM dexamethasone, 2.5 µM rosiglitazone). Chondrogenic differentiation was carried out in an antiadhesive 96-well plate with a U-shaped bottom (order no. 734–2782, VWR, Germany) in a pellet culture for 2 weeks. A medium of DMEM-LG, 1% antibiotics, 1% ITS, 0.1 µM dexamethasone, 0.9 mM sodium pyruvate, 0.17 mM ascorbic acid, 0.35 mM proline, and 10 ng/ml TGFb was used to induce chondrogenic differentiation.
After the incubation period, alizarin-red staining was used to detect osteogenic differentiation, oil red O staining was used to detect adipogenic differentiation and alcian-blue staining was used to detect chondrogenic differentiation after fixation and sectioning of the pellets [23].
AdMSC treatment
Chondrogenic differentiation
Chondrogenic differentiation was performed in a monolayer on cell culture plates with d = 10 cm (order no. 734–2321, VWR, Germany). The cells were incubated for 21 days with the chondrogenic differentiation medium (DMEM-LG, 1% antibiotics, 1% ITS, 0.1 µM dexamethasone, 0.9 mM sodium pyruvate, 0.17 mM ascorbic acid, 0.35 mM proline, 10 ng/ml TGFb).
Chondrogenic differentiation under hypoxia
In addition, chondrogenic differentiation was performed in a monolayer under hypoxia at 3% O2 as previously described.
Interleukin 1β treatment
For treatment with interleukin 1β (IL-1β, order no. 200-01B, Thermo Fisher Scientific, Germany), AdMSCs at approximately 80% confluence were incubated in a 75 cm² cell culture flask in DMEM LG, 1% P/S and 1% ITS for 3 days. For the actual treatment, 10 µg/ml IL-1β (manufacturer) was separately added to the cell culture medium every day.
Shockwave treatment
Treatment with shock waves was performed with an extracorporeal shock wave device (Piezovet 100 plus, Richard Wolf GmbH, Germany) and the FB12 G5 transducer (Richard Wolf GmbH, Germany). In a focal volume of 59 mm², an energy flux density of 1.1 mJ/mm² and a pressure of 122 MPa is applied, which corresponds to the maximum treatment level of “20”. The experimental setup is shown in the Figs. 1 and 2.
Illustration of the experimental setup for the shock wave treatment of AdMSCs. The silicone cup with the sample tube (marked with A, see details in Fig. 2), the transducer (marked with B) and the shock wave generator (marked with C)
Detailed enlargement of Fig. 1, A. The silicone cup, which is filled with water for sound wave transmission, is colored gray. The sample tube with the cell suspension (red liquid) was immersed in it. The tube is fixed in place with a green foam holder. For better visualization, the tube was pulled slightly upward from the water
Since the cells for this treatment must be present in a cell suspension to concentrate them in the focus volume, the cells obtained after passage 2 were first adjusted to a concentration of 1 × 106 cells/ml. Then, 9 ml of this suspension (containing 9 million cells) was dispensed into three 5 ml polypropylene tubes (order no. 55.526.006, Sarstedt, Germany) with 3 ml per tube. The tubes were placed in a custom silicone cup filled with water at the focus of the shock waves. The tubes were then sequentially exposed to 1000 pulses at a frequency of 8 Hz at the maximum level of the device (level 20). The AdMSCs were then recultured together in a 75 cm² cell culture flask in standard medium for 8 h. Afterwards, the FCS-containing medium was removed by washing with a PBS solution, and a further incubation period of two days followed by standard medium without FCS, supplemented with 1% ITS instead.
Induction of senescence
For the induction of senescence, the AdMSCs were further cultured until passage 20 when they reached their Hayflick limit. At passage 20, the AdMSCs were switched to FCS-free medium supplemented with 1% ITS at approximately 60% confluence and cultured for another 2 days.
In addition, senescence in a 6-well plate (order no. 734–2323, VWR, Germany) with confluent cells was confirmed after staining for beta-galactosidase. For this purpose, AdMSCs were fixed in 6-well plates with 4% paraformaldehyde, washed with PBS and then incubated overnight at 37 °C with beta-galactosidase staining solution. This staining solution consisted of a PBS solution with 150 mM NaCl, 1 mg/ml x-GAL, 5 mM K-hexacyanoferrate II, 5 mM K-hexacyanoferrate III, 2 mM MgCl2, and 40 mM citric acid. After staining, the excess dye was again rinsed off with PBS.
Isolation and collection of EVs
After the AdMSCs were cultured in serum-free medium, EVs were harvested from the cell culture supernatant. For this purpose, the EV-containing supernatant was removed from the cell culture flasks, which corresponded to a volume of 15 ml per 75 cm² flask. Since three 75 cm² flasks were used to obtain a minimum amount of EVs, a total of 45 ml of medium was obtained for each type of pretreatment (standard medium, hypoxia, IL-1β, shock wave, chondrogenic differentiation under normoxia and hypoxia) and transferred to a 50 ml centrifuge tube. These samples were centrifuged for 5 min at 2,700 × g to remove cell debris, and the supernatant was removed again and filtered through a 0.2 μm syringe filter (order no. 83.1826.001, Sarstedt, Germany).
Nanoparticle analysis and statistics
For quality assurance, 3 ml of this filtrate was then measured via a nanoparticle tracking analysis (NTA) device to ensure that the measured samples actually corresponded to the EVs. This measurement was performed via a NanoSight LM10 (Malvern Instruments Ltd., UK) and NTA 3.3 software (Malvern Instruments Ltd., UK) and is described in detail elsewhere [24]. Characterization by transmission electron microscopy (TEM) in conjunction with immunogold labeling for the tetraspanins CD9 (monoclonal antibody clone HI9a) and CD81 (monoclonal antibody clone sc-166,029) as well as western blotting of the tetraspanin CD9, was performed previously [24]. The differences in the size and number of nanoparticles obtained were analyzed for significance via one-way analysis of variance (ANOVA) with repeated measurements. If there were significant group differences, a Tukey test was then carried out to identify the significantly altered condition.
NGS analysis
A 20 ml volume was taken from the remaining supernatant and stored separately at -80 °C until the samples were collected and sent to a commercial provider for next-generation sequencing to detect the contained miRNA. For the detection of miRNAs a “TrueQuant SmallRNA Seq Kit” was used according to the manufacturer’s instructions (GenXPro GmbH, Germany). The method employs a gel-free, single-tube protocol for high-yield small RNA NGS library generation. Sequencing was performed on an Illumina NextSeq500 system via 1 × 75 sequencing cycles.
NGS data processing
Unprocessed sequencing reads were adapter-trimmed and quality-trimmed via Cutadapt (version 4.6 [25]), with the arguments “-e 0.1 -O 3 -q 20 -m 14 -n 8”. FastQC (0.11.9 [26]), was used to assess the quality of the sequencing reads. Processed sequencing reads were mapped via Bowtie2 (2.4.4 [27]), on trna (with arguments: “--sensitive --local”) and mirna (with arguments: “--local --ma 1 --score-min L,0,0.9 --mp 1,1 --rdg 2,1 --rfg 2,1”) and ENSEMBL_cdna (with arguments: “--sensitive --local”) and ENSEMBL_ncrna (with arguments: “--sensitive --local”) of EquCab3.0 (Equus caballus). The mapping was performed iteratively, meaning that only those reads not mapping to the previous reference were mapped to the next one. The Quantification of the reads mapped to each transcript was performed via HTSeq (version 2.0.2 [28]), trna (with arguments: “-i transcript_id -r name -a 0 -m union” and strandedness “no”) and mirna (with arguments: “-I transcript_id -r name -a 0 -m union” and strandedness “no”) and ENSEMBL_cdna (with arguments: “-i transcript_id -r name -a 0 -m union” and strandedness “no”) and ENSEMBL_ncrna (with arguments: “-i transcript_id -r name -a 0 -m union” and strandedness “no”). MultiQC (version 1.25.1 [29]), was used to create a single report that visualizes output from multiple tools across many samples, enabling global trends and biases to be quickly identified.
Differential expression analysis (DEA), statistics
The DEA was performed with DESeq2 (version 1.38 [30]), . Only entries that had a raw value of at least 5 in at least 2 samples were used in the DEA. The Log2FoldChange values were shrunk with “ashr” [31]. The DEA results with an FDR-adjusted p-value less than or equal to 0.05 and an absolute log2FoldChange greater than or equal to 1 were labeled as significant via the Wald test. Correction of p values after multiple testing was performed via the Benjamini-Hochberg method.
Heatmap generation
A heatmap was created via an online tool [30] to provide a clear representation of the changes in miRNA expression. The average linkage method was selected as the clustering method and Euclid was selected as the distance measurement method.
Results
Cultivation and characterization of equine adipose derived MSCs
As defined for MSCs, the AdMSCs used, exhibited clear plastic adherence to the cell culture flasks (Fig. 3). Differentiation into the adipogenic, osteogenic and chondrogenic lineages was also demonstrated (Fig. 4). According to general MSC characteristics flow cytometric analysis revealed the expression of surface markers such as CD44 and CD90 and the absence of a CD45 and MHC II expression (Fig. 5).
The three directions of AdMSC differentiation are shown. A) Intracellular staining of fat vacuoles with Oil Red O. B) Histological section through a 3D chondro-pellet after chondrogenic differentiation. The sulfated proteoglycans were stained with Alcian blue. C) Distinct nodule formation after osteogenic differentiation with calcium accumulation clearly stained red by Alizarin staining. The scale bar corresponds to 100 μm
Proof of senescence
Compared with those at lower passages, positive beta-galactosidase detection in AdMSCs at passage 20, as detected by clear bluish coloration, demonstrated the senescence characteristic of stem cells (Fig. 6).
NTA analysis of the media supernatant
After the cell culture supernatants enriched with EVs were obtained, the size and number of nanoparticles contained in the supernatants were measured via NTA. Since no clear distinction can be made between EVs and other cellular or other components, the nanoparticles contained in the medium, which also contain the EVs, must be measured at this point. The average size under all conditions (negative control, IL-1β stimulation, chondrogenic differentiation with and without hypoxia, shock wave treatment, and senescence) was relatively the same with a size of 146.01 ± 7.67 nm (x̅ ± SD, Fig. 7). There were no significant differences between the different treatment options (p = 0.075). On the other hand, there were differences in the number of nanoparticles obtained (Fig. 8). Here, the average concentration was 9.66 ± 5.39 × 108 particles/ml (x̅ ± SD). A lower concentration of nanoparticles after chondrogenic differentiation with 2.35 ± 1.19 particles/ml (x̅ ± SD) was particularly evident, compared to a significantly increased concentration after chondrogenic differentiation under hypoxia. Here, 20.23 ± 2.36 particles/ml (x̅ ± SD) were present in the cell culture supernatant. A significant difference of p = 0.029 was found between these two concentrations. In general, however, there were no significant differences between the various treatments and the negative control, i.e. untreated AdMSC, at a concentration of 8.70 ± 1.67 × 108 particles/ml (x̅ ± SD).
Particle size of EVs in the supernatant of cultured equine AdMSCs after the different treatments. The different conditions previously described are as follows: negative control without prior treatment (control), chondrogenic differentiation (chondro diff), chondrogenic differentiation under hypoxia (hypoxia), treatment with IL-1β, shockwave treatment (shockwave) and senescent AdMSCs (senescence)
EV concentrations in the supernatants of cultured equine AdMSCs after different treatments. The different conditions previously described are as follows: negative control without prior treatment (control), chondrogenic differentiation (chondro diff), chondrogenic differentiation under hypoxia (hypoxia), treatment with IL-1β, shockwave treatment (shockwave) and senescent AdMSCs (senescence). The asterisk indicates a significant difference (p = 0.029)
Compared with the respective treatment methods, the size of the respective nanoparticles was no longer extremely homogeneous. This is shown in Fig. 9.
Particle size distributions are visualized after analysis via a nanoparticle tracking analysis (NTA) instrument. Measurements were carried out after different pretreatments: negative control without prior treatment (A), chondrogenic differentiation (B), chondrogenic differentiation under hypoxia (C), treatment with IL-1β (D), shockwave treatment (E) and senescent AdMSCs (F)
NGS anaylsis
After next-generation sequencing and raw data processing, the average sequencing depth per sample was 4,872,746 reads. On average, 51.04% of the reads mapped to trna, 0.13% of the reads mapped to mirna, 2.22% of the reads mapped to ENSEMBL_cdna and 11.16% of the reads mapped to ENSEMBL_ncrna. Figure 10 shows the heatmap with the corresponding significantly altered expression. A total of 89 significantly different miRNAs were found. Among these, one to 10 miRNAs were increased, and 44 to 68 miRNAs were decreased, depending on the treatment of the AdMSCs. In particular, chondrogenic differentiation under hypoxia showed notable results, where only one miRNA was upregulated (eca-miR-1307) and 68 miRNAs were expressed at low levels (Fig. 11). An overview of the significant changes in expression after the various treatment options compared with untreated AdMSCs is shown in Tables 1, 2, 3, 4 and 5 for chondrogenic differentiation (Table 1), chondrogenic differentiation under hypoxia (Table 2), treatment with IL-1β (Table 3), shock wave treatment (Table 4) and senescence (Table 5).
Heatmap of significantly altered miRNA expressions after different treatments of equine AdMSCs. These were chondrogenic differentiation (chondro diff), stimulation with interleukin 1β (IL-1ß), shockwave treatment (shockwave), chondrogenic differentiation under hypoxia (hypoxia) and senescence in AdMSCs (senescence)
Representation of the number of significantly changed miRNA expressions after different treatments of equine AdMSCs. The sums of the significantly changed miRNA expressions compared with the negative control (NC) are shown. The various conditions previously described are as follows: negative control without prior treatment (control), chondrogenic differentiation (chondro diff), chondrogenic differentiation under hypoxia (hypoxia), treatment with IL-1β (IL 1ß), shockwave treatment (shockwave) and senescent AdMSCs (senescence)
Discussion
Since MSCs are considered to act in part through the secretion of EVs, we focused on the characterization of miRNAs within these EVs. For this purpose, equine AdMSCs were stimulated or damaged in different ways in vitro. The change in the resulting gene expression in the EVs should provide information on their supportive action toward a healing process by finding potential candidates for miRNAs via overall analysis (next generation sequencing) of the miRNAs contained in the EVs.
The equine AdMSCs used for these studies have already been obtained and characterized in advance, which is why the methods were not described in detail. Nevertheless, the donors used in this study provided the evidence required by the International Society for Cellular Therapy (ISCT) standards [20]. In addition, proof of senescence was provided for more passages via beta-galactosidase staining. In stem cells as primary cell cultures, a passage above P3 is already considered senescent [33, 34], although there may be strong donor-dependent variations [35, 36].
The various treatment options for AdMSCs before harvesting EVs should, if possible, correspond to potential typical treatment modalities associated with OA as carried out in the equine clinic. In addition to a control of chondrogenic differentiation, equine AdMSCs were subjected to IL-1β stimulation, shock wave treatment and senescence. The aim was to reproduce different scenarios of the actual state or treatment:
Chondrogenically differentiated AdMSCs under normoxia and hypoxia for the treatment of OA
-
1.
IL-1β as a damage model for OA.
-
2.
Shock wave as a treatment for OA.
-
3.
Senescence as a possible cause or additional problem of the use of MSCs in the treatment of OA.
For chondrogenic differentiation, a monolayer culture was used instead of a 3D pellet culture. This was because no EVs could be obtained in significant quantities in the media supernatants of the pellet culture normally used for chondrogenic differentiation (data not shown). This ruled out this method for subsequent NGS analysis. Since hypoxic conditions prevail in the center of the pellets, instead of culturing the cells under normoxia, chondrogenic differentiation was performed in the monolayer under hypoxic culture conditions via an incubator with O2 regulation at an oxygen concentration of 3%. This was intended to replace the conditions of the pellet culture as much as possible.
After the cell culture supernatants of the treated equine AdMSCs were collected, the number and size distribution of the contained particles (including EVs) were measured via NTA for quality control. The size of the particles contained was found to be consistent with previous studies, regardless of pretreatment. The average EV size of 146 nm was within the expected range for EVs [24, 31, 32]. In contrast, the particle concentration was significantly more heterogeneous with 9.66 ± 5.39 × 108 (x̅ ± SD) particles/ml. The significantly higher number of particles after chondrogenic differentiation under hypoxia is a described phenomenon that can be attributed to hypoxia [33, 34].
Finally, the media supernatants were sent to a service provider for NGS. A total of 89 significantly altered miRNAs were found after different treatments. These were determined from different groups after treatment of AdMSCs. As seen previously in particle counts, most miRNAs were also found after chondrogenic differentiation with hypoxia, which was significantly different from the results of the negative control. We found a total of 69 significantly altered miRNAs, of which only miR-1307 was downregulated, and all other miRNAs were upregulated. The increased expression of miRNAs stimulated by hypoxia has been described previously, and here, we identified a specific set of miRNAs, the so-called hypoxemic, which are particularly highly expressed in cells stressed by hypoxia [35]. Many miRNAs were found to be up- or downregulated under hypoxia. However, compared with our results, some of the miRNAs expected to be downregulated, namely the miRNAs eca-let-7a/c/d/f, eca-miR-101, eca-miR-16, eca-miR-186, eca-miR-199a(-5p), eca-miR-20a and eca-miR-374(a), were actually upregulated. Among the more highly expressed miRNAs is the miRNA eca-miR-424. This miRNA is known to be a key miRNA in hypoxia, as it stabilizes hypoxia-inducible factor 1 (HIF-1) and thus supports the associated cell response to hypoxia [36].
Of greatest interest, however, are certainly those miRNAs that are thought to have a protective or rather deleterious effect in OA. A central idea was that AdMSCs are stimulated by pretreatment to secrete EVs with miRNA, which has a positive effect on damaged cells such as chondrocytes. However, it should be noted that the results reported in the literature are from a variety of fluids such as synovia [43, 44] and plasma [44, 45] or from different cells in vitro, such as from primary chondrocyte cultures [46] or synovial cell cultures with MSCs [47, 48]. A list of different protective and destructive miRNAs has already been published by several authors [49,49,50,51,53]. Notably, not only can an isolated miRNA be important, but protective and destructive miRNAs are also balanced in the healthy state [52].
Thus, within the category of anti-inflammatory miRNAs, miR-101 [37], miR-143 [38], mir-145 [18], miR-146a [39], miR-27a [40] and miR-93 [41] were identified. A comparison of these miRNAs revealed that eca-miR-101 was significantly increased under all conditions, and eca-miR-146a was significantly increased only in AdMSCs stimulated with IL-1β, a miRNA that is thought to attenuate IL-1β responses [39]. In addition to anti-inflammatory miRNAs, other protective miRNAs have identified, such as miR-98, which is supposed to prevent apoptosis of chondrocytes [42], and miR-29b(-5p), which induces chondrogenesis by preventing the senescence of chondrocytes [43]. The significantly elevated expression of miR-221 has also been shown to play a supportive role in chondrogenesis [18]. In a comprehensive analysis of blood serum and synovia from horse joints with induced OA, Anderson et al. 2022 were able to postulate further miRNAs important for OA with the help of an Ingenuity Pathway Analysis, of which we were able to find a total of 23 miRNAs in our investigations, one of which has already been described here (eca-miR-93), and the following [44]: eca-let-7a, eca-let-7c, eca-miR-103, eca-miR-107a, eca-miR-10a, eca-miR-10b, eca-miR-128, eca-miR-16, eca-miR-199a-3p, eca-miR-199b-3p, eca-miR-19b, eca-miR-21, eca-miR-23a, eca-miR-25, eca-miR-26a, eca-miR-28-5p, eca-miR-30d, eca-miR-423-3p, eca-miR-532-5p, eca-miR-744, eca-miR-92a and eca-miR-99b.
In addition to the numerous miRNAs whose expression increased after the various AdMSC treatments, one miRNA, namely eca-miR-1307, was also significantly reduced in all the samples. This miRNA has often been described as a tumor marker that supports the growth of tumors and can also induce resistance to chemotherapeutic agents [45, 46]. This is thought to occur by blocking apoptosis and improving DNA repair [47]. Why this miRNA in particular showed reduced gene expression in all samples cannot be interpreted without further investigations.
In summary, many miRNAs, whose effects are still unknown or have not been sufficiently investigated, have been identified by our studies specifically and exclusively in horses. At this point, we would also like to mention the small number of animals examined for NGS. It was not possible to include more donors in this study due to the high effort involved. We tried to compensate for this as much as possible by selecting three heterogeneous donors. Due to the comparatively high number of significantly altered miRNA values, we assume that the results provide a good overview for the equine species. The applicability of the results from the equine to other species, such as other domestic mammals or primates including humans, cannot be conclusively guaranteed. However, horses can be considered model animals for musculoskeletal diseases [48, 49]. In addition, 140 equine miRNAs that are 100% identical to human miRNAs have already been identified [50]. Among these homologous miRNAs, we were able to find a total of 41 different miRNAs in our studies (Table 6).
Thus, with respect to the altered gene expression of miRNAs, it is not possible to define in general terms which miRNAs would be potentially beneficial in an OA treatment without further in vivo studies. It can be argued that unstressed AdMSCs may already produce a potent reservoir of secreted EVs, but after these treatments, their potential seems to be greatly increased. To take this into account, great attention should be given to the expression profiles after different treatments, as these profiles provide information about the conditions that may prevail in vivo in a disease. Therefore, it would be conceivable to eliminate unwanted miRNAs specifically at the site of action. This could be accomplished either with suitable, possibly synthetically produced miRNA carriers, which are protective miRNAs, or directly with anti-miRNAs [51].
Conclusion
With our studies on the EVs of equine AdMSCs after different types of stimulation, we were able to find significant changes in the expression of several miRNAs, which should be further investigated in the near future. These findings provide a valuable basis for the search for corresponding miRNAs for therapeutic purposes in humans and animals or for further deciphering the function of AdMSCs.
Data availability
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Abbreviations
- AdMSC:
-
Adipose derived mesenchymal stem cells
- ISCT:
-
International Society for Cellular Therapy
- ITS:
-
Insulin-transferrin-selenite
- NGS:
-
Next Generation Sequencing
- OA:
-
Osteoarthritis
- SD:
-
Standard derivation
References
Arden N, Nevitt MC. Osteoarthritis: epidemiology. Best Pract Res Clin Rheumatol. 2006;20:3–25. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.berh.2005.09.007.
Anderson KL, O’Neill DG, Brodbelt DC, Church DB, Meeson RL, Sargan D, et al. Prevalence, duration and risk factors for appendicular osteoarthritis in a UK dog population under primary veterinary care. Sci Rep. 2018;8:5641. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-018-23940-z.
Brommer H, van Weeren PR, Brama PAJ, Barneveld A. Quantification and age-related distribution of articular cartilage degeneration in the equine fetlock joint. Equine Vet J. 2003;35:697–701. https://doiorg.publicaciones.saludcastillayleon.es/10.2746/042516403775696357.
Yunus MHM, Nordin A, Kamal H. Pathophysiological perspective of osteoarthritis. Med (Kaunas). 2020. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/medicina56110614.
Martel-Pelletier J, Barr AJ, Cicuttini FM, Conaghan PG, Cooper C, Goldring MB, et al. Osteoarthr Nat Rev Dis Primers. 2016;2:16072. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nrdp.2016.72.
Jones AC, Doherty M. The treatment of osteoarthritis. Br J Clin Pharmacol. 1992;33:357–63. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1365-2125.1992.tb04052.x.
Rychel JK. Diagnosis and treatment of osteoarthritis. Top Companion Anim Med. 2010;25:20–5. https://doiorg.publicaciones.saludcastillayleon.es/10.1053/j.tcam.2009.10.005.
Nelson AE, Allen KD, Golightly YM, Goode AP, Jordan JM. A systematic review of recommendations and guidelines for the management of osteoarthritis: the chronic osteoarthritis management initiative of the U.S. Bone and joint initiative. Semin Arthritis Rheum. 2014;43:701–12. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.semarthrit.2013.11.012.
Ferket BS, Feldman Z, Zhou J, Oei EH, Bierma-Zeinstra SMA, Mazumdar M. Impact of total knee replacement practice: cost effectiveness analysis of data from the osteoarthritis initiative. BMJ. 2017;356:j1131. https://doiorg.publicaciones.saludcastillayleon.es/10.1136/bmj.j1131.
Goldberg A, Mitchell K, Soans J, Kim L, Zaidi R. The use of mesenchymal stem cells for cartilage repair and regeneration: a systematic review. J Orthop Surg Res. 2017;12:39. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13018-017-0534-y.
Caplan AI. Mesenchymal stem cells: time to change the name! Stem Cells Transl Med. 2017;6:1445–51. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/sctm.17-0051.
Rani S, Ryan AE, Griffin MD, Ritter T. Mesenchymal stem Cell-derived extracellular vesicles: toward Cell-free therapeutic applications. Mol Therapy: J Am Soc Gene Therapy. 2015;23:812–23. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/mt.2015.44.
Record M, Carayon K, Poirot M, Silvente-Poirot S. Exosomes as new vesicular lipid transporters involved in cell-cell communication and various pathophysiologies. Biochim Biophys Acta. 2014;1841:108–20. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.bbalip.2013.10.004.
Flamant S, Loinard C, Tamarat R. MSC beneficial effects and limitations, and MSC-derived extracellular vesicles as a new cell-free therapy for tissue regeneration in irradiated condition. Environ Adv. 2023;13:100408. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.envadv.2023.100408.
Zaborowski MP, Balaj L, Breakefield XO, Lai CP. Extracellular vesicles: composition, biological relevance, and methods of study. Bioscience. 2015;65:783–97. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/biosci/biv084.
Cavallo C, Merli G, Borzì RM, Zini N, D’Adamo S, Guescini M, et al. Small extracellular vesicles from adipose derived stromal cells significantly attenuate in vitro the NF-κB dependent inflammatory/catabolic environment of osteoarthritis. Sci Rep. 2021;11:1053. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-020-80032-7.
Wu X, Wang Y, Xiao Y, Crawford R, Mao X, Prasadam I. Extracellular vesicles: potential role in osteoarthritis regenerative medicine. J Orthop Translat. 2020;21:73–80. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jot.2019.10.012.
Zhao C, Chen J-Y, Peng W-M, Yuan B, Bi Q, Xu Y-J. Exosomes from adipose–derived stem cells promote chondrogenesis and suppress inflammation by upregulating miR–145 and miR–221. Mol Med Rep. 2020;21:1881–9. https://doiorg.publicaciones.saludcastillayleon.es/10.3892/mmr.2020.10982.
Raabe O, Reich C, Wenisch S, Hild A, Burg-Roderfeld M, Siebert H-C, Arnhold S. Hydrolyzed fish collagen induced chondrogenic differentiation of equine adipose tissue-derived stromal cells. Histochem Cell Biol. 2010;134:545–54. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00418-010-0760-4.
Dominici M, Le Blanc K, Mueller I, Slaper-Cortenbach I, Marini F, Krause D, et al. Minimal criteria for defining multipotent mesenchymal stromal cells. The international society for cellular therapy position statement. Cytotherapy. 2006;8:315–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/14653240600855905.
Leisengang S, Heilen LB, Klymiuk MC, Nürnberger F, Ott D, Wolf-Hofmann K, et al. Neuroinflammation in primary cultures of the rat spinal dorsal Horn is attenuated in the presence of adipose Tissue-Derived medicinal signalling cells (AdMSCs) in a Co-cultivation model. Mol Neurobiol. 2022;59:475–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s12035-021-02601-9.
Raabe O, Shell K, Fietz D, Freitag C, Ohrndorf A, Christ HJ, et al. Tenogenic differentiation of equine adipose-tissue-derived stem cells under the influence of tensile strain, growth differentiation factors and various oxygen tensions. Cell Tissue Res. 2013;352:509–21. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00441-013-1574-1.
Arnhold S, Elashry MI, Klymiuk MC, Geburek F. Investigation of stemness and multipotency of equine adipose-derived mesenchymal stem cells (ASCs) from different fat sources in comparison with lipoma. Stem Cell Res Ther. 2019;10:309. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13287-019-1429-0.
Klymiuk MC, Balz N, Elashry MI, Heimann M, Wenisch S, Arnhold S. Exosomes isolation and identification from equine mesenchymal stem cells. BMC Vet Res. 2019;15:42. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12917-019-1789-9.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10. https://doiorg.publicaciones.saludcastillayleon.es/10.14806/ej.17.1.200.
Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nmeth.1923.
Anders S, Pyl PT, Huber W. HTSeq—a python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/bioinformatics/btu638.
Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32:3047–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/bioinformatics/btw354.
Babicki S, Arndt D, Marcu A, Liang Y, Grant JR, Maciejewski A, Wishart DS. Heatmapper: web-enabled heat mapping for all. Nucleic Acids Res. 2016;44:W147–53. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkw419.
Welsh JA, Goberdhan DCI, O’Driscoll L, Buzas EI, Blenkiron C, Bussolati B, et al. Minimal information for studies of extracellular vesicles (MISEV2023): from basic to advanced approaches. J Extracell Vesicles. 2024;13:e12404. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/jev2.12404.
Doyle LM, Wang MZ. Overview of Extracellular Vesicles, Their Origin, Composition, Purpose, and Methods for Exosome Isolation and Analysis. Cells. 2019. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/cells8070727
Borges FT, Melo SA, Özdemir BC, Kato N, Revuelta I, Miller CA, et al. TGF-β1-containing exosomes from injured epithelial cells activate fibroblasts to initiate tissue regenerative responses and fibrosis. JASN. 2013;24:385–92. https://doiorg.publicaciones.saludcastillayleon.es/10.1681/ASN.2012101031.
Panigrahi GK, Praharaj PP, Peak TC, Long J, Singh R, Rhim JS, et al. Hypoxia-induced exosome secretion promotes survival of African-American and Caucasian prostate cancer cells. Sci Rep. 2018;8:3853. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-018-22068-4.
Nallamshetty S, Chan SY, Loscalzo J. Hypoxia: a master regulator of MicroRNA biogenesis and activity. Free Radic Biol Med. 2013;64:20–30. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.freeradbiomed.2013.05.022.
Loscalzo J. The cellular response to hypoxia: tuning the system with MicroRNAs. J Clin Invest. 2010;120:3815–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1172/JCI45105.
Dai L, Zhang X, Hu X, Liu Q, Man Z, Huang H, et al. Silencing of miR-101 prevents cartilage degradation by regulating extracellular Matrix-related genes in a rat model of osteoarthritis. Mol Ther. 2015;23:1331–40. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/mt.2015.61.
Qiu B, Xu X, Yi P, Hao Y. Curcumin reinforces MSC-derived exosomes in attenuating osteoarthritis via modulating the miR-124/NF-kB and miR-143/ROCK1/TLR9 signalling pathways. J Cell Mol Med. 2020;24:10855–65. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/jcmm.15714.
Yang C-R, Shih K-S, Liou J-P, Wu Y-W, Hsieh I-N, Lee H-Y, et al. Denbinobin upregulates miR-146a expression and attenuates IL-1β-induced upregulation of ICAM-1 and VCAM-1 expressions in osteoarthritis fibroblast-like synoviocytes. J Mol Med. 2014;92:1147–58. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00109-014-1192-8.
Qiu W-J, Xu M-Z, Zhu X-D, Ji Y-H. MicroRNA-27a alleviates IL-1β-induced inflammatory response and articular cartilage degradation via TLR4/NF-κB signaling pathway in articular chondrocytes. Int Immunopharmacol. 2019;76:105839. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.intimp.2019.105839.
Jin Z, Ren J, Qi S, RETRACTED. Human bone mesenchymal stem cells-derived exosomes overexpressing microRNA-26a-5p alleviate osteoarthritis via down-regulation of PTGS2. Int Immunopharmacol. 2020;78:105946. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.intimp.2019.105946.
Zheng W-D, Zhou F-L, Lin N, Liu J. Investigation for the role of CTX-III and microRNA-98 in diagnosis and treatment of osteoarthritis. Eur Rev Med Pharmacol Sci. 2018;22:5424–8. https://doiorg.publicaciones.saludcastillayleon.es/10.26355/eurrev_201809_15801.
Zhu J, Yang S, Qi Y, Gong Z, Zhang H, Liang K, et al. Stem cell-homing hydrogel-based miR-29b-5p delivery promotes cartilage regeneration by suppressing senescence in an osteoarthritis rat model. Sci Adv. 2022;8:eabk0011. https://doiorg.publicaciones.saludcastillayleon.es/10.1126/sciadv.abk0011.
Anderson JR, Jacobsen S, Walters M, Bundgaard L, Diendorfer A, Hackl M, et al. Small non-coding RNA landscape of extracellular vesicles from a post-traumatic model of equine osteoarthritis. Front Vet Sci. 2022;9:901269. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fvets.2022.901269.
Hu Y, Qi C, Liu X, Zhang C, Gao J, Wu Y, et al. Malignant ascites-derived exosomes promote peritoneal tumor cell dissemination and reveal a distinct MiRNA signature in advanced gastric cancer. Cancer Lett. 2019;457:142–50. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.canlet.2019.04.034.
Saberianpour S, Abkhooie L. MiR-1307: A comprehensive review of its role in various cancer. Gene Rep. 2021;25:101392. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.genrep.2021.101392.
Carotenuto P, Amato F, Lampis A, Rae C, Hedayat S, Previdi MC, et al. Modulation of pancreatic cancer cell sensitivity to FOLFIRINOX through microRNA-mediated regulation of DNA damage. Nat Commun. 2021;12:6738. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41467-021-27099-6.
Gregory MH, Capito N, Kuroki K, Stoker AM, Cook JL, Sherman SL. A review of translational animal models for knee osteoarthritis. Arthritis. 2012;2012:764621. https://doiorg.publicaciones.saludcastillayleon.es/10.1155/2012/764621.
McIlwraith CW, Frisbie DD, Kawcak CE. The horse as a model of naturally occurring osteoarthritis. Bone Joint Res. 2012;1:297–309. https://doiorg.publicaciones.saludcastillayleon.es/10.1302/2046-3758.111.2000132.
Myrzabekova MO, Labeit SB, Nyiazova R. In Silico identification of equus Caballus MiRNAs with the potential to affect human gene expression. EB. 2022;90:60–72. https://doiorg.publicaciones.saludcastillayleon.es/10.26577/eb.2022.v90.i1.06.
Yin H, Xiong G, Guo S, Xu C, Xu R, Guo P, Shu D. Delivery of Anti-miRNA for Triple-Negative breast cancer therapy using RNA nanoparticles targeting stem cell marker CD133. Mol Ther. 2019;27:1252–61. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ymthe.2019.04.018.
Szala D, Kopańska M, Trojniak J, Jabłoński J, Hanf-Osetek D, Snela S, Zawlik I. The role of MicroRNAs in the pathophysiology of osteoarthritis. Int J Mol Sci 2024. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms25126352
Acknowledgements
We gratefully thank the German Research Foundation (DFG) for the financial and personnel support of the research project.
Funding
Open Access funding enabled and organized by Projekt DEAL.
This study was funded by the German Research Foundation (DFG, grant No. AR 333/11 − 1 and WE 4319/4 − 1).
Author information
Authors and Affiliations
Contributions
MCK: Manuscript writing, JS and IDM: Performing the experiments; MIE: Manuscript revision, statistical analysis; SW and SA: Manuscript revision and study design. All the authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The authorization for the tissue collection was obtained at the 16th of May 2018 from the local authorities (Regierungspräsidium Gießen, Wetzlar, Germany) and is registered under the number V 54 − 19 c 20 15 h 02 GI 18/1 kTV 1/2018. There is no special project name for this permission.
Consent for publication
Not applicable.
Artificial intelligence (AI)
The authors declare that they have not use AI-generated work in this manuscript.
Competing interests
Not applicable.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Klymiuk, M.C., Speer, J., Marco, I.D. et al. Determination of the miRNA profile of extracellular vesicles from equine mesenchymal stem cells after different treatments. Stem Cell Res Ther 16, 162 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13287-025-04287-5
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13287-025-04287-5