# Association of circulating hsa-miRNAs with sarcopenia: the SarcoPhAge study

Marjorie Millet, Maxime Auroux, Charlotte Beaudart, Céline Demonceau, Aurélie Ladang, Etienne Cavalier, Jean-Yves Reginster, Olivier Bruyère, Roland Chapurlat, Jean-Charles Rousseau

► To cite this version:

Marjorie Millet, Maxime Auroux, Charlotte Beaudart, Céline Demonceau, Aurélie Ladang, et al. Association of circulating hsa-miRNAs with sarcopenia: the SarcoPhAge study. Aging Clinical and Experimental Research, 2024, 36 (1), pp.70. 10.1007/s40520-024-02711-z. inserm-05028892

HAL Id: inserm-05028892

https://inserm.hal.science/inserm-05028892v1

Submitted on 10 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Association of circulating hsa-miRNAs with sarcopenia: the SarcoPhAge study

Marjorie Millet<fnref n="1" />,<fnref n="2" /> · Maxime Auroux<fnref n="1" />,<fnref n="3" /> · Charlotte Beaudart<fnref n="4" />,<fnref n="5" /> · Céline Demonceau<fnref n="5" /> · Aurélie Ladang<fnref n="6" /> ·

Etienne Cavalier<fnref n="6" /> · Jean-Yves Reginster<fnref n="5" /> · Olivier Bruyère<fnref n="5" /> · Roland Chapurlat<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="7" /> · Jean-Charles Rousseau<fnref n="1" />,<fnref n="2" />

Received: 10 August 2023 / Accepted: 23 January 2024

© The Author(s) 2024

## Abstract

**Objective** To identify a microRNA signature associated to sarcopenia in community-dwelling older adults from the SarcoPhAge cohort.

**Methods** In a screening phase by next generation sequencing (NGS), we compared the hsa-miRome expression of 18 subjects with sarcopenia (79.6 ± 6.8 years, 9 men) and 19 healthy subjects without sarcopenia (77.1 ± 6 years, 9 men) at baseline. Thereafter, we have selected eight candidate hsa-miRNAs according to the NGS results and after a critical assessment of previous literature. In a validation phase and by real-time qPCR, we then analyzed the expression levels of these 8 hsa-miRNAs at baseline selecting 92 healthy subjects (74.2 ± 10 years) and 92 subjects with sarcopenia (75.3 ± 6.8 years). For both steps, the groups were matched for age and sex.

**Results** In the validation phase, serum has-miRNA-133a-3p and has-miRNA-200a-3p were significantly decreased in the group with sarcopenia vs controls [RQ: relative quantification; median (interquartile range)]: -0.16 (-1.26/+0.90) vs +0.34 (-0.73/+1.33) (p < 0.01) and -0.26 (-1.07/+0.68) vs +0.27 (-0.55/+1.10) (p < 0.01) respectively. Has-miRNA-744-5p was decreased and has-miRNA-151a-3p was increased in the group with sarcopenia vs controls, but this barely reached significance: +0.16 (-1.34/+0.79) vs +0.44 (-0.31/+1.00) (p = 0.050) and +0.35 (-0.22/+0.90) vs +0.03 (-0.68/+0.75) (p = 0.054).

**Conclusion** In subjects with sarcopenia, serum hsa-miRNA-133a-3p and hsa-miRNA-200a-3p expression were downregulated, consistent with their potential targets inhibiting muscle cells proliferation and differentiation.

**Keywords** Sarcopenia · Circulating hsa-miRNAs · Next generation sequencing · TaqMan qPCR

<figure><img src="image_2.png" /></figure>

## Introduction

Sarcopenia is a disease defined by progressive decrease in skeletal muscle mass leading to loss of strength and function [1]. Consequently, the individuals with sarcopenia are at increased risk of falls, functional decline, frailty, and mortality [2] Moreover, the amount of metabolically active tissue is reduced by sarcopenia; thus, the risk for metabolic diseases is increased including cardiovascular disease, diabetes, hypertension and hyperlipidemia [3, 4]. The overall prevalence of sarcopenia is estimated between 6 and 22% in adults aged 65 years and increases with age, making sarcopenia a major public health problem [5].

The complex pathogenesis of sarcopenia, however, remains incompletely understood. Moreover, because of several definitions and inadequate screening tools, sarcopenia remains underdiagnosed. Consequently, there is a

<footnote>Jean-Charles Rousseau
jean-charles.rousseau@inserm.fr</footnote>

<footnote n="1">1 INSERM 1033, Lyon, France</footnote>

<footnote n="2">2 PMO, Lyon, France</footnote>

<footnote n="3">3 Hôpital E. Herriot, Hospices Civils de Lyon, Lyon, France</footnote>

<footnote n="4">4 Clinical Pharmacology and Toxicology Research Unit
(URPC), NARILIS, Department of Biomedical Sciences,
Faculty of Medicine, University of Namur, Namur, Belgium</footnote>

<footnote n="5">5 WHO Collaborating Center for Epidemiology
of Musculoskeletal Health and Aging, Division of Public
Health, Epidemiology and Health Economics, University
of Liège, Liege, Belgium</footnote>

<footnote n="6">6 Department of Clinical Chemistry, CHU de Liège, University
of Liège, Liege, Belgium</footnote>

<footnote n="7">7 University of Lyon, Lyon, France</footnote>
considerable interest in the identification of specific biologic
markers that reflect quantitative and dynamic variations
in muscle tissue remodeling [6]. Potential targets include
microRNAs (miRNAs). Previous studies have shown that
hsa-miRNAs play a pivotal role in the regulation of the
major age-associated alterations in muscle pathophysiology
including skeletal cell proliferation, differentiation, apopto-
sis and senescence [2]. These small non-coding RNAs of
22–28 nucleotides in length can silence gene expression by
binding to complementary sequence on target messenger
RNA transcripts resulting in translational repression or tar-
get degradation [7]. One miRNA can inhibit a large number
of mRNAs, and one gene can be targeted by multiple miR-
NAs. Circulating miRNAs are easily accessible and quantifi-
able with high degree of sensitivity and specificity. Finally,
the remarkable miRNA stability in biofluids suggests they
could become non-invasive disease biomarkers.

So, the aim of our study was to identify a microRNA sig-
nature associated with sarcopenia. For this purpose, we have
studied the SarcoPhAge cohort, a prospective longitudinal
study following community-dwelling older subjects includ-
ing 534 participants aged 65 years or older [8], in an initial
screening phase and a subsequent validation phase.

Methods

Study population

The group for the screening and the validation phases included Belgian participants from the SarcoPhAge study (Sarcopenia and Physical impairments with advancing Age). The main methodology of the SarcoPhAge study has already been described in detail elsewhere [8]. Briefly, the SarcoPhAge cohort is a population-based cohort developed in 2013 in Liège (Belgium). The participants (n = 534) were community-dwelling, aged 65 years and older and recruited from press advertisements and general, geriatric, osteoporosis, rehabilitation and rheumatology departments from an outpatient clinic in Liège, Belgium. No specific exclusion criteria related to health or demographic characteristics were applied, except the exclusion of individuals with an amputated limb or with a body mass index (BMI) above 50 kg/m², which was required for X-ray absorptiometry. The study was approved by the Ethics Committee of the Teaching Hospital of the University of Liège (reference 2012/277) with two amendments in 2015 and 2018. All volunteers gave their written informed consent. The participants were followed up annually and a clinical research assistant performed physical examinations and health questionnaires to gather sociodemographic and anamnestic data. The diagnosis of sarcopenia was performed according to the initial definition of 2010

from the European Working Group on Sarcopenia in Older
People (EWGSOP) [9]. Muscle mass [Skeletal muscle
index measured with daily-calibrated dual energy X-ray
absorptiometry (DXA); cut-offs of 7.26 kg/m² for men
and 5.5 kg/m² for women] [10], muscle strength (hand-
grip strength, measured with daily-calibrated JAMAR
dynamometer; cut-offs of 30 kg for men and 20 kg for
women) [11] and physical performance [short physical
performance battery (SPPB); cut-off for the diagnosis of
sarcopenia 8 points or less out of a maximum of 12 points]
[12] were measured in all the subjects (Table 1).

The hsa-miRNA analysis by next generation sequencing (NGS) : the screening phase

Selection of patients

The expression levels of hsa-miRNAs in serum were meas-
ured in a subset of participants, 19 healthy subjects without
sarcopenia and 18 subjects with sarcopenia. The number of
samples used in this screening phase is comparable to stud-
ies previously carried out in our laboratory [13–15] and by
other authors [16–21]. Moreover, we have selected a repre-
sentative sub-group of our cohort in order to detect by NGS
the miRs with expression significantly different between
controls and sarcopenic subjects. Consequently, we do not
know which miRs will be of interest, nor the magnitude of
its expression between the two groups. As a result, we can-
can not calculate a sample size. Both groups were matched for
age and sex without any other selection criteria. The flow
chart for the design of the experimental study is reported
in Fig. 1.

NGS analysis

Total RNA was extracted from 400 µl of serum with the hsa-miRCURY Biofluids extraction kit (Exiqon®, Danemark) and analyzed by small RNA sequencing on a NextSeq500 sequencing instrument (Illumina platform). After RNA conversion into hsa-miRNA NGS libraries using NEBNEXT library generation kit, cDNA was pre-amplified prior to library purification and quantification (for detailed procedures of extraction, library preparation and purification, normalization and quality controls, see supplementary data 1). Measurements were expressed as tags per million (TPM) in which the number of reads for a particular hsa-miRNA is divided by the total number of mapped reads in a sample and multiplied by 10⁶. The miRNA level was further normalized by the trimmed mean of *m* values (TMM) method and compared between groups by the TMM ratio expressed as relative data as Log₂ (Fold change) (LogFC).

<figure><img src="image_3.png" /></figure>
<table><caption>**Table 1** Baseline characteristics of individuals in the screening phase (NGS analysis) (A) and the validation phase (qPCR analysis) (B)</caption><thead><tr><th></th><th>Sarcopenia, n=18</th><th>No sarcopenia, n=19</th><th>p value</th></tr></thead><tbody><tr><td>(A) Screening</td><td></td><td></td><td></td></tr><tr><td>Age (years)</td><td>79.6±6.8</td><td>76.4±6.3</td><td>0.14</td></tr><tr><td>Sex (women)</td><td>9 (50%)</td><td>10 (53%)</td><td>0.87</td></tr><tr><td>SPPB cut-off</td><td></td><td></td><td></td></tr><tr><td>SPPB > 8 points</td><td>6 (33.3%)</td><td>18 (94.7%)</td><td>&lt;0.0001</td></tr><tr><td>SPPB ≤ 8 points</td><td>12 (66.7%)</td><td>1 (5.3%)</td><td></td></tr><tr><td>Grip strength maximum (kg)</td><td></td><td></td><td></td></tr><tr><td>Women</td><td>20.0±4.4</td><td>35.2±11.1</td><td>0.007<sup>a</sup></td></tr><tr><td>Men</td><td>22.7±4.6</td><td>43.6±7.6</td><td>&lt;0.0001</td></tr><tr><td>Grip strength cut-off</td><td></td><td></td><td></td></tr><tr><td>Grip strength > 20 or 30 kg</td><td>0 (0.0%)</td><td>18 (94.7%)</td><td>&lt;0.0001<sup>b</sup></td></tr><tr><td>Grip strength < 20 or 30 kg</td><td>18 (100.0%)</td><td>1 (5.3%)</td><td></td></tr><tr><td>SMI (kg/m<sup>2</sup>)</td><td></td><td></td><td></td></tr><tr><td>Women</td><td>4.9±0.3</td><td>7.1±0.7</td><td>&lt;0.0001</td></tr><tr><td>Men</td><td>6.4±0.6</td><td>8.7±1.1</td><td>0.0002</td></tr><tr><td>SMI cut-off</td><td></td><td></td><td></td></tr><tr><td>SMI > 7.26 or 5.5 kg/m<sup>2</sup></td><td>0 (0.0%)</td><td>19 (9.0%)</td><td>&lt;0.0001<sup>b</sup></td></tr><tr><td>SMI < 7.26 or 5.5 kg/m<sup>2</sup></td><td>18 (100%)</td><td>0 (0.0%)</td><td></td></tr><tr><td>(B) validation</td><td></td><td></td><td></td></tr><tr><td>Age (years)</td><td>75.8±6.6</td><td>75.4±6.2</td><td>0.71</td></tr><tr><td>Sex (women)</td><td>54 (59%)</td><td>55 (60%)</td><td>0.88</td></tr><tr><td>SPPB cut-off</td><td></td><td></td><td></td></tr><tr><td>SPPB > 8 points</td><td>39 (42.4%)</td><td>72 (78.3%)</td><td>&lt;0.0001</td></tr><tr><td>SPPB ≤ 8 points</td><td>53 (42.4%)</td><td>20 (21.7%)</td><td></td></tr><tr><td>Grip strength maximum (kg)</td><td></td><td></td><td></td></tr><tr><td>Women</td><td>17.3±3.1</td><td>22.1±6.1</td><td>&lt;0.0001<sup>a</sup></td></tr><tr><td>Men</td><td>29.0±7.2</td><td>38.2±8.9</td><td>&lt;0.0001</td></tr><tr><td>Grip strength cut-off</td><td></td><td></td><td></td></tr><tr><td>Grip strength > 20 or 30 kg</td><td>14 (15.2%)</td><td>66 (71.7%)</td><td>&lt;0.0001</td></tr><tr><td>Grip strength < 20 or 30 kg</td><td>78 (84.8%)</td><td>26 (28.3%)</td><td></td></tr><tr><td>SMI kg/m<sup>2</sup></td><td></td><td></td><td></td></tr><tr><td>Women</td><td>5.1±0.3</td><td>6.4±0.9</td><td>&lt;0.0001<sup>a</sup></td></tr><tr><td>Men</td><td>6.6±0.5</td><td>8.0±0.9</td><td>&lt;0.0001</td></tr><tr><td>SMI cut-off</td><td></td><td></td><td></td></tr><tr><td>SMI > 7.26 or 5.5 kg/m<sup>2</sup></td><td>0 (0.0%)</td><td>80 (87.0%)</td><td>&lt;0.0001</td></tr><tr><td>SMI < 7.26 or 5.5 kg/m<sup>2</sup></td><td>92 (100%)</td><td>12 (13.0%)</td><td></td></tr></tbody></table>

SPPB: For the physical performance, we used the short physical performance battery (SPPB) test. This test is composed of three separate tests: balance, 4-m gait speed and chair stand tests. A score between 0 and 4 is assigned for each test, and the 3 tests are weighted equally. As suggested, the cut-off for the diagnosis of sarcopenia for this test, is a score of 8 points or less out of a maximum of 12 points [12]. Gait speed is also proposed by the EWGSOP [9] as a tool to assess physical performance. Given the fact that gait speed is one of the components of the SPPB test, we chose to use the SPPB test as diagnostic tool

Grip strength: as recommended by the EWGSOP [9], we measured subjects' handgrip strength to determine their muscle strength. Therefore, we used a hydraulic dynamometer (acquired from Saehan Corporation, MSD Europe Bvba, Belgium), calibrated at the beginning of the study for 10, 40 and 90 kg, that subjects had to squeeze as hard as possible three times with each hand (dominant and nondominant). We used the highest result out of the six measurements recorded in our analysis [11]. For the diagnosis of sarcopenia, we used the recommended cut-offs of 30 kg for men and 20 kg for women

SMI: Skeletal muscle index calculated by dividing the appendicular skeletal muscle mass (ASM) by the height squared. The proposed cut-off of 7.26 kg/m² for men and 5.5 kg/m² for women were used [10]

<footnote><sup>a</sup>Mann-Whitney U test performed due to the non-normality</footnote>

<footnote><sup>b</sup>Fisher's exact test performed because N<5</footnote>
<figure><img src="image_4.png" /><figcaption>Fig. 1 Flow diagram describing the methodology</figcaption></figure>

# The hsa-miRNA analysis by real-time quantitative polymerase chain reaction: the validation phase

## Selection of patients

Expression levels of serum hsa-miRNA were measured in 92 subjects with sarcopenia corresponding to the total number of individuals with sarcopenia in the SarcoPhAge cohort and in 92 healthy subjects matched for age and sex.

## Preselection of hsa-miRNAs

In the validation phase, we have planned to validate the hsa-miRNAs significantly differentially expressed in the screening phase between patients with sarcopenia and controls. In this screening phase, the prespecified threshold of significance *p* < 0.05 was used after false discovery rate correction using the Benjamini–Hochberg method. Then, we selected 6 additional hsa-miRNA by critical assessment of literature to highlight hsa-miRNAs that are differentially expressed in biofluids of patients with sarcopenia compared with controls [2, 22–28] respectively: 21-5p, 133a-5p, 146a-5p, 151a-3p, 486-5p and 744-5p (for miRNA sequences, see supplementary Table 1).

# Hsa-miRNAs quantification by real-time quantitative polymerase chain reaction (RT-qPCR)

Total RNA was extracted from 100 µl serum with the Nucleospin miRNA plasma kit (Macherey–Nagel®) according to

manufacturer's protocol. A lysis buffer solution containing 1 µl of glycogen as RNA carrier and a synthetic spike-in control RNA (cel-miRNA-39-3p, 0.1 ng) as exogenous control was added to serum. After protein precipitation and washes, RNA was eluted from the microcolumn and stored at –80 °C. The cDNA was synthesized from serum total RNA using a TaqMan Advanced miRNA cDNA synthesis kit (Applied Biosystems, ThermoFisher Scientific). The qPCR amplification was performed on cDNA obtained by miRNA-Amp reaction of reverse transcription reaction, using Fast Advanced Master Mix and TaqMan Advanced miRNA Assays (Applied Biosystems, ThermoFisher Scientific). Hsa-miRNAs quantification was performed by the TaqMan® Advanced miRNA technology (Applied Biosystems, ThermoFisher Scientific) in duplicate of 8 hsa-miRNAs target, 3 hsa-miRNAs endogeneous controles and 1 miRNA as exogenous control by real-time PCR reaction on a Quantstudio 7 Flex (Applied Biosystems) according to the manufacturer's protocol. We used the software Expression Suite (Applied Biosystems) to express the hsa-miRNA level as relative quantification (RQ). The threshold cycle (Ct) value of each hsa-miRNA were normalized with the mean of expression level of three endogenous controls. RQ was calculated as 2⁻²⁻ΔΔCT, with ΔCT = (CT hsa-miRNA – CT mean of endogenous controls) and ΔΔCT = (ΔCT of the hsa-miRNA –ΔCT mean of the hsa-miRNAs through all samples) and converted as FC = Log2(2⁻²⁻ΔΔCT). Data were normalized with the mean of expression level of three endogenous hsa-miRNAs: hsa-miR-191-5p, hsa-miR-222-3p and hsa-miR-361-5p, that are known to be ubiquitously
expressed and without reported impact on sarcopenia disease or bone [13, 14]. The exogenous spike cel-miRNA-39-3p was used as a qPCR quality control (for protocol details, see supplementary data 1 and for data normalization, see supplementary Fig.1).

## Statistical analysis

### NGS (screening phase) and qPCR (validation phase) analysis

In the NGS step, analysis of hsa-miRNAs was performed by Qiagen using the EdgeR software package. For RT-qPCR, RQ was calculated with Expression-Suite Software. A hsa-miRNA with a *p* value ≤ 0.05 and a false discovery rate of 5% (Benjamini-Hochberg FDR correction for NGS approach) was considered as differentially expressed. Wilcoxon tests were used to compare hsa-miRNA levels between subjects with and without sarcopenia because of the skewness of the data.

### Baseline characteristics of subjects included in this study

Normality of the variables was assessed using a Kolmogorov-Smirnov test, observation of the difference between the mean and the median, histogram and QQ plot. Qualitative variables are presented as numbers (*N*) and percentages (%) and continuous variables are expressed as medians and quartiles. We compared the characteristics between healthy subjects and individuals with sarcopenia using the Mann-Whitney *U* test for quantitative variables. Statistical significance was defined as a two-tailed *p* value < 0.05.

The analyses, graphics and figures were performed using GraphPad Prism version 8.4.3 for Windows, GraphPad Software, San Diego, California, USA, http://www graph pad.com.

## Results

Baseline characteristics of participants are displayed in Table 1. There was no sizeable difference between these age- and sex-matched individuals.

### Screening: serum hsa-miRNA profiling of subjects with sarcopenia and control individuals

Expression levels of serum hsa-miRNA were measured in 19 healthy subjects without sarcopenia (77.1 ± 6 years, 9 men) and in 18 subjects with sarcopenia (79.6 ± 6.8 years, 9 men). We identified 383 has-miRs with an expression level ≥ 1 TPM and 196 with an expression level ≥ 10 TPM in all samples. When we compared both groups, 43 hsa-miRNAs showed differential expression (*p* < 0.05)

between controls and individuals with sarcopenia, 31 up-regulated and 12 down-regulated in the group with sarcopenia vs controls. After the Benjamini-Hochberg FDR correction, hsa-miRNA-200a-3p and hsa-miRNA-668-3p remained significantly downregulated in the group with sarcopenia compared to the control group (*p* < 0.05, FDR at 5%) (Table 2).

### Validation: differential expression of candidate hsa-miRNAs in serum of subjects with sarcopenia and control individuals

Expression levels of serum hsa-miRNA were measured in 92 healthy subjects without sarcopenia (74.2 ± 10 years) and in the 92 subjects with sarcopenia (75.3 ± 6.8 years). Both groups were matched for age (*p* = 0.71) and sex (*p* = 0.88). There was no significant difference between the two groups, so we made no further adjustment. Serum has-miRNA-133a-3p and has-miRNA-200a-3p were significantly decreased in the group with sarcopenia vs controls (RQ: relative quantification; median [Interquartile range]: -0.16 (-1.26/+0.90) vs +0.34 (-0.73/+1.33) (*p* < 0.01) and -0.26 (-1.07/+0.68) vs +0.27 (-0.55/+1.10) (*p* < 0.01), respectively. Has-miRNA-744-5p was decreased and has-miRNA-151a-3p was increased in the group with sarcopenia compared to the control group, but this barely reached significance: +0.16 (-1.34/+0.79) vs +0.44 (-0.31/+1.00) (*p* = 0.050) and +0.35 (-0.22/+0.90) vs +0.03 (-0.68/+0.75) (*p* = 0.054) respectively (Table 3 and Fig. 2).

## Discussion

We have evaluated by real-time qPCR whether the two differentially expressed hsa-miRNAs (200a-3p and 668-3p) in subjects with sarcopenia observed in the screening step were also aberrantly expressed in the validation cohort. In addition, we measured the expression levels of 6 additional hsa-miRNAs (21-5p, 133a-5p, 146a-5p, 151a-3p, 486-5p and 744-5p) previously highlighted in prior studies as having a dysregulated expression in subjects suffering from sarcopenia.

### Screening phase

The NGS approach revealed that 43 hsa-miRNAs had serum levels of expression significantly different between the group with sarcopenia and the healthy subjects. After FDR correction, only 2 of them remained significant downregulated, hsa-miRNA-200a-3p and 668-3p. Of the remaining 41 hsa-miRNAs, twelve were downregulated and twenty-nine
<table><caption>**Table 2** Differential expression of the miR levels in the serum from healthy subjects and subjects with sarcopenia analyzed by next generation sequencing</caption><thead><tr><th>miRNAs</th><th>Control (TMM)</th><th>Sarcopenic (TMM)</th><th>logFC</th><th>p value</th><th>FDR</th></tr></thead><tbody><tr><td>hsa-miR-668-3p</td><td>3.820</td><td>0.230</td><td>3.287</td><td>1.720 E-06</td><td>1.018 E-03</td></tr><tr><td>hsa-miR-200a-3p</td><td>146.330</td><td>50.500</td><td>1.535</td><td>5.739E-05</td><td>1.699 E-02</td></tr><tr><td>hsa-miR-6833-3p</td><td>0.250</td><td>2.850</td><td>-3.122</td><td>4.772E-04</td><td>9.416 E-02</td></tr><tr><td>hsa-miR-1246</td><td>135.040</td><td>293.180</td><td>-1.116</td><td>9.589 E-04</td><td>1.419 E-01</td></tr><tr><td>hsa-miR-4638-5p</td><td>0.370</td><td>3.450</td><td>-3.053</td><td>1.250 E-03</td><td>1.480 E-01</td></tr><tr><td>hsa-miR-509-3-5p</td><td>0.100</td><td>2.330</td><td>-3.702</td><td>1.568 E-03</td><td>1.547 E-01</td></tr><tr><td>hsa-miR-6821-3p</td><td>0.040</td><td>1.380</td><td>-3.593</td><td>1.935 E-03</td><td>1.637 E-01</td></tr><tr><td>hsa-miR-4433-3p</td><td>4.040</td><td>11.630</td><td>-1.514</td><td>2.608 E-03</td><td>1.930 E-01</td></tr><tr><td>hsa-miR-885-5p</td><td>41.430</td><td>14.990</td><td>1.468</td><td>3.516 E-03</td><td>2.313 E-01</td></tr><tr><td>hsa-miR-150-5p</td><td>526.010</td><td>1271.250</td><td>-1.273</td><td>3.965 E-03</td><td>2.348 E-01</td></tr><tr><td>hsa-miR-3591-5p</td><td>1.870</td><td>0.380</td><td>1.989</td><td>8.071 E-03</td><td>3.690 E-01</td></tr><tr><td>hsa-miR-31-5p</td><td>0.140</td><td>1.510</td><td>-2.841</td><td>8.625 E-03</td><td>3.690 E-01</td></tr><tr><td>hsa-miR-140-3p</td><td>1288.650</td><td>2243.880</td><td>-0.800</td><td>8.795 E-03</td><td>3.690 E-01</td></tr><tr><td>hsa-miR-181b-5p</td><td>248.470</td><td>377.000</td><td>-0.602</td><td>8.915 E-03</td><td>3.690 E-01</td></tr><tr><td>hsa-miR-29b-3p</td><td>15.010</td><td>38.780</td><td>-1.368</td><td>9.450 E-03</td><td>3.690 E-01</td></tr><tr><td>hsa-miR-146b-5p</td><td>1362.870</td><td>2484.160</td><td>-0.866</td><td>9.974 E-03</td><td>3.690 E-01</td></tr><tr><td>hsa-miR-15b-3p</td><td>64.730</td><td>105.750</td><td>-0.707</td><td>1.261 E-02</td><td>4.350 E-01</td></tr><tr><td>hsa-miR-4664-3p</td><td>5.360</td><td>1.410</td><td>1.837</td><td>1.460 E-02</td><td>4.350 E-01</td></tr><tr><td>hsa-miR-1285-3p</td><td>13.510</td><td>25.050</td><td>-0.901</td><td>1.477 E-02</td><td>4.350 E-01</td></tr><tr><td>hsa-miR-181a-5p</td><td>3903.140</td><td>5312.120</td><td>-0.445</td><td>1.558 E-02</td><td>4.350 E-01</td></tr><tr><td>hsa-miR-4669</td><td>1.390</td><td>8.030</td><td>-2.422</td><td>1.567 E-02</td><td>4.350 E-01</td></tr><tr><td>hsa-miR-629-5p</td><td>603.340</td><td>861.850</td><td>-0.515</td><td>1.656 E-02</td><td>4.350 E-01</td></tr><tr><td>hsa-miR-29b-2-5p</td><td>3.100</td><td>0.750</td><td>1.857</td><td>1.690 E-02</td><td>4.350 E-01</td></tr><tr><td>hsa-miR-155-5p</td><td>206.090</td><td>302.870</td><td>-0.555</td><td>1.798 E-02</td><td>4.435 E-01</td></tr><tr><td>hsa-miR-1273 g-3p</td><td>4.910</td><td>1.630</td><td>1.570</td><td>1.900 E-02</td><td>4.500 E-01</td></tr><tr><td>hsa-miR-203a</td><td>74.690</td><td>219.060</td><td>-1.554</td><td>2.141 E-02</td><td>4.874 E-01</td></tr><tr><td>hsa-miR-186-5p</td><td>1138.270</td><td>1749.830</td><td>-0.620</td><td>2.374 E-02</td><td>5.204 E-01</td></tr><tr><td>hsa-miR-142-5p</td><td>2176.960</td><td>3522.360</td><td>-0.694</td><td>2.768 E-02</td><td>5.852 E-01</td></tr><tr><td>hsa-miR-181c-3p</td><td>4.700</td><td>12.170</td><td>-1.345</td><td>3.131 E-02</td><td>6.202 E-01</td></tr><tr><td>hsa-miR-205-5p</td><td>22.210</td><td>85.530</td><td>-1.949</td><td>3.143 E-02</td><td>6.202 E-01</td></tr><tr><td>hsa-miR-1292-5p</td><td>26.360</td><td>42.390</td><td>-0.682</td><td>3.579 E-02</td><td>6.407 E-01</td></tr><tr><td>hsa-miR-101-3p</td><td>901.930</td><td>1587.000</td><td>-0.815</td><td>3.662 E-02</td><td>6.407 E-01</td></tr><tr><td>hsa-miR-543</td><td>1318.830</td><td>810.030</td><td>0.703</td><td>3.699 E-02</td><td>6.407 E-01</td></tr><tr><td>hsa-miR-211-5p</td><td>0.340</td><td>2.720</td><td>-2.737</td><td>3.742 E-02</td><td>6.407 E-01</td></tr><tr><td>hsa-miR-671-5p</td><td>1.500</td><td>4.220</td><td>-1.378</td><td>3.868 E-02</td><td>6.407 E-01</td></tr><tr><td>hsa-miR-29c-3p</td><td>10.360</td><td>20.190</td><td>-0.958</td><td>3.896 E-02</td><td>6.407 E-01</td></tr><tr><td>hsa-miR-107</td><td>126.140</td><td>191.000</td><td>-0.599</td><td>4.272 E-02</td><td>6.692 E-01</td></tr><tr><td>hsa-miR-485-3p</td><td>112.840</td><td>53.340</td><td>1.081</td><td>4.497 E-02</td><td>6.692 E-01</td></tr><tr><td>hsa-miR-485-5p</td><td>151.120</td><td>95.240</td><td>0.668</td><td>4.561 E-02</td><td>6.692 E-01</td></tr><tr><td>hsa-miR-224-5p</td><td>140.780</td><td>89.380</td><td>0.654</td><td>4.741 E-02</td><td>6.692 E-01</td></tr><tr><td>hsa-miR-378a-3p</td><td>660.340</td><td>927.940</td><td>-0.491</td><td>4.752 E-02</td><td>6.692 E-01</td></tr><tr><td>hsa-miR-342-5p</td><td>396.460</td><td>555.760</td><td>-0.487</td><td>4.761 E-02</td><td>6.692 E-01</td></tr><tr><td>hsa-miR-382-5p</td><td>216.140</td><td>141.400</td><td>0.610</td><td>4.861 E-02</td><td>6.692 E-01</td></tr></tbody></table>

When we compared the sarcopenic vs healthy groups, 43 miRs showed differential expression (p < 0.05) between healthy subjects and subjects with sarcopenia. After Benjamini-Hochberg false discovery rate (FDR) correction, hsa-miR-200a-3p and hsa-miR-668-3p remained significantly different between healthy subjects and subjects with sarcopenia (p < 0.05, FDR at 5%) (in bold)

upregulated in the group with sarcopenia compared to the healthy group (see Table 2). Most of them are known to be associated with multiple pathways of muscle metabolism.

A first group may be linked to vascular biology dysregulation. Atherosclerosis may be influenced by hsa-miRNAs-155-5p [29], hsa-miRNA-150-5p [30], hsa-miRNA-382-5p
<table><caption>**Table 3** Differential expression of the miR level in the serum from healthy subjects and subjects with sarcopenia analyzed by quantitative real-time PCR</caption><thead><tr><th rowspan="2">Target</th><th colspan="4">Sarcopenic group</th><th colspan="4">Control group</th><th rowspan="2">p value</th></tr><tr><th>Median</th><th>25% percentile</th><th>75% percentile</th><th>n</th><th>Median</th><th>25% percentile</th><th>75% percentile</th><th>n</th></tr></thead><tbody><tr><td>hsa.miR.21.5p</td><td>-0.14</td><td>-0.51</td><td>0.35</td><td>92</td><td>0.06</td><td>-0.47</td><td>0.74</td><td>92</td><td>0.085</td></tr><tr><td><strong>hsa.miR.133a.3p</strong></td><td><strong>-0.16</strong></td><td><strong>-1.26</strong></td><td><strong>0.90</strong></td><td><strong>85</strong></td><td><strong>0.34</strong></td><td><strong>-0.73</strong></td><td><strong>1.33</strong></td><td><strong>90</strong></td><td><strong>0.009</strong></td></tr><tr><td>hsa.miR.146a.5p</td><td>0.05</td><td>-0.72</td><td>0.93</td><td>92</td><td>-0.11</td><td>-0.61</td><td>0.86</td><td>92</td><td>0.572</td></tr><tr><td><strong>hsa.miR.151a.3p</strong></td><td><strong>0.35</strong></td><td><strong>-0.22</strong></td><td><strong>0.90</strong></td><td><strong>92</strong></td><td><strong>0.03</strong></td><td><strong>-0.68</strong></td><td><strong>0.75</strong></td><td><strong>92</strong></td><td><strong>0.054</strong></td></tr><tr><td><strong>hsa.miR.200a.3p</strong></td><td><strong>-0.26</strong></td><td><strong>-1.07</strong></td><td><strong>0.68</strong></td><td><strong>92</strong></td><td><strong>0.27</strong></td><td><strong>-0.55</strong></td><td><strong>1.10</strong></td><td><strong>92</strong></td><td><strong>0.009</strong></td></tr><tr><td>hsa.miR.486.5p</td><td>0.27</td><td>-0.67</td><td>1.05</td><td>92</td><td>-0.07</td><td>-0.79</td><td>0.80</td><td>92</td><td>0.108</td></tr><tr><td>hsa.miR.668.3p</td><td>0.01</td><td>-1.77</td><td>2.14</td><td>82</td><td>0.34</td><td>-2.11</td><td>2.73</td><td>82</td><td>0.710</td></tr><tr><td><strong>hsa.miR.744.5p</strong></td><td><strong>0.16</strong></td><td><strong>-1.34</strong></td><td><strong>0.79</strong></td><td><strong>90</strong></td><td><strong>0.44</strong></td><td><strong>-0.31</strong></td><td><strong>1.00</strong></td><td><strong>90</strong></td><td><strong>0.050</strong></td></tr></tbody></table>

The significant miRs are in bold

<figure><img src="image_5.png" /><figcaption>**Fig. 2** The serum expression levels of the eight miRs analyzed in the validation phase in individual groups of subjects with or without sarcopenia after normalization with the three control miRs and standard-</figcaption></figure>

ized to cel-miRNA-39. Data are shown using a box-and-whisker plot where boxes show the 25th, 50th (median) and 75th percentiles, and the ends of the whiskers indicate the maximum and minimum

[31], hsa-miRNA-205-5p [32, 33] and hsa-miRNA-211-5p [34]. Hsa-miRNA 146-5p [35] and hsa-miRNA-485-3p may affect vascular calcification [36]. Migration, proliferation and differentiation of smooth muscle cells for hsa-miRNA 101-3p, 342-5p and 1246 [37–39] and diseases related to airway smooth muscle cells, asthma and chronic obstructive pulmonary disease may relate to hsa-miRNA 101-3p, 140-3p, 155-3p, 224-5p and 509-5p [40–44].

The second group of hsa-miRNAs included those involved in heart disease as a cardioprotective mechanism for hsa-miRNA 342-5p [45], anti-cardiomyocyte apoptosis after infarction for hsa-miRNA 140-3p [46], inhibition of collagen formation after a myocardial injury for hsa-miRNA 142-5p [47] and hsa-miRNA 485-5p [48], atrial fibrillation for hsa-miRNA 671-5p [49] and chronic heart failure for hsa-miRNA 1285-3p [50].

The third group encompasses hsa-miRNAs related to myoblast metabolism, such as hsa-miRNA 378-3p

influencing vascularization in skeletal muscle [51], 885-5p for determining muscle type-specific tissue formation and maintenance [52], hsa-miRNA 203 for muscle differentiation and development [53, 54], hsa-miRNA 205-5p and has-miRNA-342-5p for muscle proliferation and migration of human aortic vascular smooth muscle cells [32, 38] and differentiation of vascular smooth muscle cells regarding hsa-miRNA 142-5p [55].

Finally, the last group includes hsa-miRNAs known to be involved in sarcopenia and exercise as hsa-miRNA 107, hsa-miRNA 29b, hsa-miRNA 181a-5p, hsa-miRNA 15b-3p and hsa-miRNA 186-5p [2]. Moreover, gene ontology (GO) enrichment analysis revealed that bone mineralization, muscle hypertrophy, vascular smooth muscle cell differentiation, regulation of IL-6 synthesis, protein kinase activity and axon guidance are among the main pathways significantly associated with differentially expressed hsa-miRNAs (see supplementary Table 2).
## Validation phase

In the validation step, we analyzed the serum expression levels of 8 hsa-miRNAs at baseline selecting 92 subjects with sarcopenia and 92 healthy subjects matched for age and sex. Among the eight selected hsa-miRNAs, four of them, 21-5p, 146a-5p, 486-5p and 668-3p were not significantly dysregulated between both groups.

In contrast, in the subjects with sarcopenia, the expression of hsa-miRNA-133a-3p and hsa-miRNA-200a-3p were significantly down regulated while the expression of hsa-miRNA-744-5p and hsa-miRNA-151a-3p were down and upregulated respectively but barely significant.

## Potential explanations for the hsa-miRNAs dysregulations in the validation phase

Recently, circulating hsa-miRNAs 146a-5p and 486-5p have been shown to be downregulated in a group of subjects with sarcopenia compared to those of normal and dynapenic groups [22]. The targets of hsa-miRNA-486 suggest that it may modulate the progressive loss of muscle mass in older adults [56, 57]. The hsa-miRNA-146a-5p is an anti-inflammatory hsa-miRNA that negatively regulates the inflammatory response [58]. Moreover, it also modulates cellular senescence and the mitochondrial metabolism [59]. However, we did not find any associations with sarcopenia for hsa-miRNAs 146a-5p and 486-5p in the SarcoPhAge cohort. It could be argued that the subjects did not belong to the same ethnical group, Chinese for Liu et al. and Caucasian in our study. Moreover, the sarcopenia definitions used were different as well as the normalization process for the qPCR quantification.

For the circulating hsa-miRNA-21-5p, four previous studies did not find any association with sarcopenia which is consistent with our findings [22, 60–62]. Comparisons are challenging, however, because the definitions of sarcopenia were not the same, conventional definition for Chen et al. and Asian Working Group for Sarcopenia (AWGS) for the three others and the quantification techniques were different, SYBGreen for He et al. and locked nucleic acid (LNA) probe for Liu et al. in plasma and serum for Chen et al. It could be noticed that hsa-miRNA-21 targets Forkhead box protein O1 (FoxO1) and phosphatase and TENsin homolog (PTEN), two genes involved in sarcopenia dysregulation [63, 64].

Finally, hsa-miRNA-688 has not been quantified in circulation before and we did not manage to quantify it when we measured its level in the largest validation cohort. This may be due to the very low level of its expression in the screening phase, under 4 TPM. However, in silico analysis suggested that hsa-miRNA-668-3p targets nuclear factor-kappa B inhibitor alpha (NFKBIA) that inhibits the transcription factor NF-kB involved in the myogenesis regulation [65]. It

also participates in the regulation of skeletal muscle growth by the inhibition of insulin-like growth factor 2 receptor (IGF2R) expression [66].

To sum up, this lack of significance for these four hsa-miRNAs can be attributed to differences between our study and the previous ones or we can speculate that the metabolic pathways dysregulated by these hsa-miRNAs do not affect sarcopenia as evaluated in the SarcoPhAge cohort.

However, four hsa-miRNAs had expression significantly different between the two groups. In subjects with sarcopenia, hsa-miRNA-133a-3p and hsa-miRNA-200a-3p expression was downregulated. Both are known to promote myoblast differentiation and proliferation. The hsa-miRNA-133a is primarily produced by skeletal muscle and it belongs to the myomiRNA family. It promotes muscle differentiation and proliferation [28]. In addition, circulating hsa-miRNA-133a was strongly associated with creatine kinase indicating that the muscle damage contributes at least in part to plasma hsa-miRNA-133a levels [25]. As for the hsa-miRNA-200a-3p, it promotes skeletal muscle satellite cells development by targeting transforming growth factor beta 2 (TGF-β2) and regulating the TFG-β2/SMAD signaling pathway [67]. Previously, Liu et al. [22, 68], did not observed a significant dysregulation between healthy subjects and subjects with sarcopenia for has-miRNA-133a in human plasma. Valaskova et al., also looking at hsa-miRNA-133a and b in human plasma [53], found that the severity of muscle performance deterioration was correlated with the decrease of hsa-miRNA-133 expression, similarly to our results.

To sum up, the hsa-miRNA-133a-3p and hsa-miRNA-200a-3p decreases are therefore consistent with their known roles suggesting that these two hsa-miRNAs could be potential circulating markers of sarcopenia.

In contrast, the variations of hsa-miRNA-744-5p and hsa-miRNA-151a-3p in subjects with sarcopenia are the opposite of what we could expect. Hsa-miRNA-744-5p is down regulated while it inhibits myoblast differentiation and myogenesis. The hsa-miRNA-744-5p sequestration by circSNX29 positively upregulates the expression of Wnt family member 5a (Wnt5a) and Ca²⁺/calmodulin-dependent protein kinase II delta (CaMKIIδ) promoting myoblast differentiation [23]. In contrast, Wei et al. showed that hsa-miRNA-151a-3p promotes myoblast proliferation but also it downregulates slow muscle gene expression by targeting ATPase2, an intracellular calcium pump located in the cell sarcoplasmic reticula. This induces the shift in muscle fiber type, upregulating myosin heavy chain beta (MHC-β) stimulating muscle contraction [27].

To sum up, these variations of hsa-miRNA-151a-3p and 744-5p may reflect the adaptation of metabolic activity of muscle cells to lower muscle mass in order to maintain the steady state of muscle quality. These variations

<figure><img src="image_6.png" /></figure>
are possibly inadequate to compensate for the muscle loss
leading to sarcopenia.

Finally, to the best of our knowledge, five studies have
been previously performed concerning circulating hsa-
miRNAs in serum [60] or plasma [22, 61, 62, 68]. How-
ever, the ethnicity of participants, the criteria for sarcope-
nia evaluation, the qPCR quantification technique and the
normalization of the results were different across studies.

Strengths and limitations of the study

We have investigated the association between the serum
hsa-miRNAs expression and sarcopenia in a well-char-
acterized population of subjects, the SarcoPhAge cohort.
We used DXA to evaluate muscle mass of subjects. We
have chosen to carry out the validation step of our study
by analysis of duplicate samples on TaqMan microRNA
arrays in order to reduce handling variability to improve
the reproducibility between the samples. In addition to
the three biologic endogenous controls, we have included
a spike in to control for efficiency in the RNA isolation
procedure and cDNA synthesis reaction. Moreover, our
study combines the validation of NGS results and the first
attempt to replicate previously published results, a process
that had been lacking in the hsa-miRNA research field.
The SarcoPhAge cohort, however, mainly includes vol-
unteers, which can limit the generalizability of our results
to all women and men in Belgium. In addition, the loss of
muscle mass and poor muscle strength or physical perfor-
mance in these older adult participants may be only par-
tially attributable to physiologic aging, and the influences
of genetic selection, lifestyle, and nutritional status or the
differences in other characteristics among the two groups
cannot be excluded. Moreover, the current experimental
results have not provided direct evidence to clarify how
circulating hsa-miRNAs regulate the sarcopenic processes
in the older adults.

In conclusion, we have shown a dysregulation of
has-miRNA-133a-3p and has-miRNA-200a-3p in older
patients with sarcopenia. Our results provide two potential
new therapeutic targets, hsa-miRNA-151a-3p and 744-5p
for the treatment and/or prevention of sarcopenia.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s40520-024-02711-z.

**Author contributions** The first draft of the manuscript was written by JCR. All authors commented on previous versions of the manuscript, and read and approved the final manuscript. MA and CD performed the statistical analysis. MM carried out data collection. JCR, RC and MM conceived the miR study. CB, EC, AL, JYR and OB conceived the SarcoPhAge study.

Declarations

**Conflict of interest** The authors have no competing interests to declare that are relevant to the content of this article.

**Informed consent and ethical approval** All patients gave informed con-
sent and the SarcoPhAge study was approved by the Ethics Committee
of the University Teaching Hospital of Liège (number 2012/277).

**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/.

References

1. Cruz-Jentoft AJ, Bahat G, Bauer J et al (2019) Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing 48:601

2. Javanmardifard Z, Shahrbanian S, Mowla SJ (2021) MicroRNAs associated with signaling pathways and exercise adaptation in sarcopenia. Life Sci 285:119926

3. Petermann-Rocha F, Balntzi V, Gray SR et al (2022) Global prevalence of sarcopenia and severe sarcopenia: a systematic review and meta-analysis. J Cachexia Sarcopenia Muscle 13:86-99

4. Yilmaz O, Bahat G (2019) Muscle mass adjustment method affects association of sarcopenia and sarcopenic obesity with metabolic syndrome. Geriatr Gerontol Int 19:272

5. Ethgen O, Beaudart C, Buckinx F et al (2017) The future prevalence of sarcopenia in Europe: a claim for public health action. Calcif Tissue Int 100:229-234

6. Ladang A, Beaudart C, Reginster JY et al (2023) Biochemical markers of musculoskeletal health and aging to be assessed in clinical trials of drugs aiming at the treatment of sarcopenia: consensus paper from an expert group meeting organized by the European Society for Clinical and Economic Aspects of Osteoporosis, Osteoarthritis and Musculoskeletal Diseases (ESCEO) and the Centre Academique de Recherche et d'Experimentation en Sante (CARES SPRL), Under the auspices of the World Health Organization Collaborating Center for the epidemiology of musculoskeletal conditions and aging. Calcif Tissue Int 112:197-217

7. Margolis LM, Rivas DA (2018) Potential role of microRNA in the anabolic capacity of skeletal muscle with aging. Exerc Sport Sci Rev 46:86-91

8. Beaudart C, Reginster JY, Petermans J et al (2015) Quality of life and physical components linked to sarcopenia: the SarcoPhAge study. Exp Gerontol 69:103-110

9. Cruz-Jentoft AJ, Baeyens JP, Bauer JM et al (2010) Sarcopenia: European consensus on definition and diagnosis: report of the European working group on sarcopenia in older people. Age Ageing 39:412-423

10. Baumgartner RN, Koehler KM, Gallagher D et al (1998) Epidemiology of sarcopenia among the elderly in New Mexico. Am J Epidemiol 147:755-763
11. Roberts HC, Denison HJ, Martin HJ et al (2011) A review of the measurement of grip strength in clinical and epidemiological studies: towards a standardised approach. Age Ageing 40:423-429

12. Guralnik JM, Ferrucci L, Pieper CF et al (2000) Lower extremity function and subsequent disability: consistency across studies, predictive models, and value of gait speed alone compared with the short physical performance battery. J Gerontol A Biol Sci Med Sci 55:M221-M231

13. Legrand MA, Millet M, Merle B et al (2020) A signature of circulating miRNAs associated with fibrous dysplasia of bone: the mirDys study. J Bone Miner Res 35:1881-1892

14. Rousseau JC, Millet M, Croset M et al (2020) Association of circulating microRNAs with prevalent and incident knee osteoarthritis in women: the OFELY study. Arthritis Res Ther 22:2

15. Mercier-Guery A, Millet M, Merle B et al (2023) Dysregulation of microRNAs in adult osteogenesis imperfecta: the miROI study. J Bone Miner Res 38:1665-1678

16. Zhang L, Yang M, Marks P et al (2012) Serum non-coding RNAs as biomarkers for osteoarthritis progression after ACL injury. Osteoar Cartil 20:1631-1637

17. Borgonio Cuadra VM, Gonzalez-Huerta NC, Romero-Cordoba S et al (2014) Altered expression of circulating microRNA in plasma of patients with primary osteoarthritis and in silico analysis of their pathways. PLoS ONE 9:e97690

18. Beyer C, Zampetaki A, Lin NY et al (2015) Signature of circulating microRNAs in osteoarthritis. Ann Rheum Dis 74:e18

19. Li YH, Tavallaee G, Tokar T et al (2016) Identification of synovial fluid microRNA signature in knee osteoarthritis: differentiating early- and late-stage knee osteoarthritis. Osteoarth Cartil 24:1577-1586

20. Kong R, Gao J, Si Y et al (2017) Combination of circulating miR-19b-3p, miR-122-5p and miR-486-5p expressions correlates with risk and disease severity of knee osteoarthritis. Am J Transl Res 9:2852-2864

21. Ntoumou E, Tzetis M, Braoudaki M et al (2017) Serum microRNA array analysis identifies miR-140-3p, miR-33b-3p and miR-671-3p as potential osteoarthritis biomarkers involved in metabolic processes. Clin Epigenetics 9:127

22. Liu HC, Han DS, Hsu CC et al (2021) Circulating microRNA-486 and microRNA-146a serve as potential biomarkers of sarcopenia in the older adults. BMC Geriatr 21:86

23. Peng S, Song C, Li H et al (2019) Circular RNA SNX29 sponges miR-744 to regulate proliferation and differentiation of myoblasts by activating the Wnt5a/Ca(2+) signaling pathway. Mol Ther Nucleic Acids 16:481-493

24. Brown DM, Goljanek-Whysall K (2015) microRNAs: modulators of the underlying pathophysiology of sarcopenia? Ageing Res Rev 24:263-273

25. Qaisar R, Karim A, Muhammad T et al (2021) Circulating microRNAs as biomarkers of accelerated sarcopenia in chronic heart failure. Glob Heart 16:56

26. Iannone F, Montesanto A, Cione E et al (2020) Expression patterns of muscle-specific miR-133b and miR-206 correlate with nutritional status and sarcopenia. Nutrients 12:297

27. Wei H, Li Z, Wang X et al (2015) MicroRNA-151-3p regulates slow muscle gene expression by targeting ATP2a2 in skeletal muscle cells. J Cell Physiol 230:1003-1012

28. Chen JF, Mandel EM, Thomson JM et al (2006) The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet 38:228-233

29. Gonzalez-Lopez P, Ares-Carral C, Lopez-Pastor AR et al (2022) Implication of miR-155-5p and miR-143-3p in the vascular insulin resistance and instability of human and experimental atherosclerotic plaque. Int J Mol Sci 23:10253

30. Bian Y, Cai W, Lu H et al (2021) miR-150-5p affects AS plaque with ASMC proliferation and migration by STAT1. Open Med (Wars) 16:1642-1652

31. Shi Z, Zhu Q, Fan J (2021) lncRNA TUG1 promotes atherosclerosis progression by targeting miR-382-5p. Int J Clin Exp Pathol 14:972-979

32. Huang P, Zhang Y, Wang F et al (2022) MiRNA-205-5p regulates the ERBB4/AKT signaling pathway to inhibit the proliferation and migration of HAVSMCs induced by ox-LDL. Pathol Res Pract 233:153858

33. Huang P, Wang F, Zhang Y et al (2023) Icariin alleviates atherosclerosis by regulating the miR-205-5p/ERBB4/AKT signaling pathway. Int Immunopharmacol 114:109611

34. Zhang Y, Wang H, Xia Y (2021) The expression of miR-211-5p in atherosclerosis and its influence on diagnosis and prognosis. BMC Cardiovasc Disord 20:371

35. Yang J, Zhou X, Lu J et al (2022) miR-146-5p restrains calcification of vascular smooth muscle cells by suppressing TRAF6. Open Med (Wars) 17:1515-1527

36. Wei X, Su Y, Li Q et al (2021) Analysis of crucial genes, pathways and construction of the molecular regulatory networks in vascular smooth muscle cell calcification. Exp Ther Med 21:589

37. Pan D, Liu G, Li B et al (2021) MicroRNA-1246 regulates proliferation, invasion, and differentiation in human vascular smooth muscle cells by targeting cystic fibrosis transmembrane conductance regulator (CFTR). Pflugers Arch 473:231-240

38. Bi S, Peng Q, Liu W et al (2020) MicroRNA-342-5p activates the Akt signaling pathway by downregulating PIK3R1 to modify the proliferation and differentiation of vascular smooth muscle cells. Exp Ther Med 20:239

39. Park N, Kang H (2020) BMP-induced microRNA-101 expression regulates vascular smooth muscle cell migration. Int J Mol Sci 21:4764

40. Fang L, Wang X, Zhang M et al (2022) MicroRNA-101-3p Suppresses mTOR and causes mitochondrial fragmentation and cell degeneration in COPD. Can Respir J 2022:5933324

41. Meng J, Zou Y, Hou L et al (2022) MiR-140-3p ameliorates the inflammatory response of airway smooth muscle cells by targeting HMGB1 to regulate the JAK2/STAT3 signaling pathway. Cell J 24:673-680

42. Wang G, Tao X, Peng L (2022) miR-155-5p regulates hypoxia-induced pulmonary artery smooth muscle cell function by targeting PYGL. Bioengineered 13:12985-12997

43. Zhuang Z, Zhou Y, Xu J et al (2022) Anti-asthmatic miR-224-5p inhibits the FHL1/MAPK pathway to repress airway smooth muscle cell proliferation in a murine model of asthma-like airway inflammation. Allergy Asthma Clin Immunol 18:88

44. Wang J, Jiang R, Tan Y et al (2022) Human pulmonary artery smooth muscle cell dysfunction is regulated by miR-509-5p in hypoxic environment. Cell Cycle 21:1212-1221

45. Hou Z, Qin X, Hu Y et al (2019) Longterm exercise-derived exosomal miR-342-5p: a novel exerkine for cardioprotection. Circ Res 124:1386-1400

46. Wei DZ, Lin C, Huang YQ et al (2017) Ellagic acid promotes ventricular remodeling after acute myocardial infarction by up-regulating miR-140-3p. Biomed Pharmacother 95:983-989

47. Wang Z, Fu M, Li Y (2020) miR-142-5p and miR-212-5p cooperatively inhibit the proliferation and collagen formation of cardiac fibroblasts by regulating c-Myc/TP53INP1. Can J Physiol Pharmacol 98:314-323

48. Zhou J, Li L, Hu H et al (2020) Circ-HIPK2 accelerates cell apoptosis and autophagy in myocardial oxidative injury by sponging miR-485-5p and targeting ATG101. J Cardiovasc Pharmacol 76:427-436

49. Xiao S, Zhou Y, Liu Q et al (2021) Identification of pivotal microRNAs and target genes associated with persistent atrial fibrillation
based on bioinformatics analysis. Comput Math Methods Med
2021:6680211

50. Li H, Fan J, Yin Z et al (2016) Identification of cardiac-related circulating microRNA profile in human chronic heart failure. Oncotarget 7:33–45

51. Krist B, Podkalicka P, Mucha O et al (2020) miR-378a influ-
ences vascularization in skeletal muscles. Cardiovasc Res
116:1386–1397

52. Muroya S, Taniguchi M, Shibata M et al (2013) Profiling of dif-
ferentially expressed microRNA and the bioinformatic target gene
analyses in bovine fast- and slow-type muscles by massively paral-
lel sequencing. J Anim Sci 91:90–103

53. Lu C, Wu J, Xiong S et al (2017) MicroRNA-203a regulates fast muscle differentiation by targeting dmrt2a in zebrafish embryos. Gene 625:49–54

54. Shen X, Cui C, Tang S et al (2022) MyoG-enhanced circGPD2 regulates chicken skeletal muscle development by targeting miR-203a. Int J Biol Macromol 222:2212–2224

55. Wu W, Shang Y, Dai S et al (2020) Downregulation of miR-142-5p inhibits human aortic smooth muscle cell proliferation and migration by targeting MKL2. Mol Med Rep 22:277–285

56. Small EM, O'Rourke JR, Moresi V et al (2010) Regulation of PI3-
kinase/Akt signaling by muscle-enriched microRNA-486. Proc
Natl Acad Sci U S A 107:4218–4223

57. Xu M, Chen X, Chen D et al (2017) FoxO1: a novel insight into its
molecular mechanisms in the regulation of skeletal muscle differ-
entiation and fiber type specification. Oncotarget 8:10662–10674

58. Park H, Huang X, Lu C et al (2015) MicroRNA-146a and micro-
RNA-146b regulate human dendritic cell apoptosis and cytokine
production by targeting TRAF6 and IRAK1 proteins. J Biol Chem
290:2831–2841

59. Jiang M, Xiang Y, Wang D et al (2012) Dysregulated expres-
sion of miR-146a contributes to age-related dysfunction of mac-
rophages. Aging Cell 11:29–40

60. Chen Z, Bemben MG, Bemben DA (2019) Bone and muscle spe-
cific circulating microRNAs in postmenopausal women based on
osteoporosis and sarcopenia status. Bone 120:271–278

61. He N, Zhang YL, Zhang Y et al (2020) Circulating microRNAs in plasma decrease in response to sarcopenia in the elderly. Front Genet 11:167

62. He N, Zhang Y, Zhang Y et al (2021) Increasing fracture risk associates with plasma circulating microRNAs in aging people's sarcopenia. Front Physiol 12:678610

63. Luo A, Yan H, Liang J et al (2017) MicroRNA-21 regu-
lates hepatic glucose metabolism by targeting FOXO1. Gene
627:194–201

64. Zhao MY, Wang LM, Liu J et al (2018) MiR-21 Suppresses anoikis through targeting PDCD4 and PTEN in human esophageal adenocarcinoma. Curr Med Sci 38:245-251

65. Bakkar N, Wang J, Ladner KJ et al (2008) IKK/NF-kappaB regu-
lates skeletal myogenesis via a signaling switch to inhibit dif-
ferentiation and promote mitochondrial biogenesis. J Cell Biol
180:787–802

66. Bella P, Farini A, Banfi S et al (2020) Blockade of IGF2R improves muscle regeneration and ameliorates Duchenne muscular dystrophy. EMBO Mol Med 12:e11019

67. Yin H, He H, Shen X et al (2020) MicroRNA profiling reveals an abundant miR-200a-3p promotes skeletal muscle satellite cell development by targeting TGF-beta2 and regulating the TGF-beta2/SMAD signaling pathway. Int J Mol Sci 21:3274

68. Valaskova S, Gazova A, Vrbova P et al (2021) The severity of muscle performance deterioration in sarcopenia correlates with circulating muscle tissue-specific miRNAs. Physiol Res 70:S91–S98

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.