Skip to main content

Identification and evaluation of circulating exosomal miRNAs for the diagnosis of postmenopausal osteoporosis

Abstract

Background

Postmenopausal osteoporosis (PMOP) is a common condition that leads to a loss of bone density and an increased risk of fractures in women. Recent evidence suggests that exosomal miRNAs are involved in regulating bone development and osteogenesis. However, exosomal miRNAs as biomarkers for PMOP diagnosis have not been systematically evaluated. In this study, we aim to identify PMOP-associated circulating exosomal miRNAs and evaluate their diagnostic performance.

Methods

We performed next-generation sequencing and bioinformatics analysis of plasma exosomal miRNAs from 12 PMOP patients and 12 non-osteoporosis controls to identify PMOP-associated exosomal miRNAs, and then validated them in an independent natural community cohort with 26 PMOP patients and 21 non-osteoporosis controls. Exosomes were isolated with the size exclusion chromatography method from the plasma of elder postmenopausal women. The plasma exosomal miRNA profiles were characterized in PMOP paired with controls with next-generation sequencing. Potential plasma exosomal miRNAs were validated by qRT-PCR in the validation cohort, and their performance in diagnosing PMOP was systematically evaluated with the receiver operating characteristic curve.

Results

Twenty-seven miRNAs were identified as differentially expressed in PMOP versus controls in sequencing data, of which six exosomal miRNAs (miR-196-5p, miR-224-5p, miR320d, miR-34a-5p, miR-9-5p, and miR-98-5p) were confirmed to be differentially expressed in PMOP patients by qRT-PCR in the validation cohort. The three miRNAs combination (miR-34a-5p + miR-9-5p + miR-98-5p) demonstrated the best diagnostic performance, with an AUC = 0.734. In addition, the number of pregnancies was found to be an independent risk factor that can improve the performance of exosomal miRNAs in diagnosing PMOP.

Conclusions

These results suggested that the plasma exosomal miRNAs had the potential to serve as noninvasive diagnostic biomarkers for PMOP.

Introduction

Osteoporosis accounts for more than 1.5 million fractures in the USA and 2.33 million fractures in China annually. Postmenopausal osteoporosis (PMOP) is the most common type of osteoporosis, characterized by decreased bone mineral density (BMD), the disintegration of bone microstructures, increased bone fragility, and easy fractures; it occurs in women after menopause as a result of increasing age and hormonal decrease [1]. As the global population ages, the incidence of PMOP rises dramatically, with about 50% of postmenopausal women worldwide affected by osteoporosis, and the prevalence of fractures in osteoporosis patients as high as 40% [2]. In 1994, a Working Group of the World Health Organization (WHO) established an operational definition of postmenopausal osteoporosis with BMD T score ≤ − 2.5 as the criteria [3]. Early detection and treatment of postmenopausal osteoporosis are essential to prevent fractures and maintain quality of life. The dual-energy X-ray absorptiometry (DXA) scan is the gold standard for diagnosing osteoporosis based on the measurement of bone mineral density (BMD) [4]. But the BMD usually show significant differences on an annual basis [5]. It is urgent to establish a more sensitive diagnosis method for osteoporosis.

Exosomes are small vesicles that are released by cells and contain a variety of biomolecules, including proteins, lipids, and RNA [6]. It participates in intercellular communication by delivering its inclusions [7]. Although different types of exosome components have been reported as markers for disease diagnosis, exosome miRNAs are the most widely reported as promising biomarkers for a range of diseases, including cancer, cardiovascular disease, and neurodegenerative diseases [8,9,10,11,12]. Compared with other components, exosomal miRNA is more abundant, and its detection method is more sensitive than protein [13]; compared to free miRNAs in the blood, exosomes can protect the miRNAs in them from degradation [10].

Previous in vitro and animal studies have shown that bone-marrow-derived mesenchymal stem cell (BMSC) derived exosomal miR-27a-3p, miR-196-5p, miR-29a, and miR-186 are involved in the mediation of bone angiogenesis and osteogenesis [14,15,16]. On the other hand, M1 macrophage-derived exosomes have been reported to aggravate bone loss in PMOP via a microRNA-98/DUSP1/JNK axis [17], and Mastocytosis-derived extracellular vesicles deliver miR-23a and miR-30a into pre-osteoblasts and prevent osteoblastogenesis and bone formation in the mouse model [18]. In addition, five plasma exosomal miRNAs, miR-224-3p, miR-25-5p, miR-302a-3p, miR-642a-3p, and miR-766-5p were identified as associated with PMOP by real-time PCR in the previous study [19]; mir-324-3p, mir-766-3p, mir-1247-5p, mir-330-5p, and mir-3124-5p have been reported to be associated with BMD in serum exosomes and may serve as potential diagnostic markers for PMOP [20]. These studies suggest that circulating exosomal miRNAs may be a promising diagnostic marker for PMOP; however, the systematic evaluation of circulating exosomal miRNAs as diagnostic markers of PMOP has not been fully investigated.

It is worth noting that in a study comparing the incidence of osteoporosis in North and South Korean women, the overall prevalence of osteoporosis was 48% (25/52) and 17% (27/156) among North Korean and South participants, while the North Korean subjects had a higher number of pregnancies[21]. Another study showed that breastfeeding was not associated with PMOP, but that the number of pregnancies increased the risk of PMOP in Chinese women [22]. This evidence suggests that the number of pregnancies may be a risk factor for PMOP.

In this study, we performed next-generation sequencing and bioinformatics analysis of plasma exosomal miRNAs from PMOP and non-osteoporosis controls, and systematically evaluated the diagnostic performance of PMOP-associated exosomal miRNAs in an independent cohort. Exosomes were isolated with the size exclusion chromatography (SEC) method from the plasma of elder postmenopausal women, whose average age was over 75 years. The plasma exosomal miRNA profiles were characterized in PMOP paired with controls with next-generation sequencing. Alterations of plasma exosomal miRNA were identified and compared with previous studies. Potential plasma exosomal miRNA markers were validated in a new independent validation cohort and their performance in diagnosing PMOP was systematically evaluated. In addition, we assessed the reliability of the number of pregnancies as a risk factor for PMOP. The purpose of this study was to identify plasma-derived exosomal miRNA biomarkers for PMOP diagnosis and evaluate their diagnostic performance.

Materials and methods

Plasma sample collection and exosome isolation

In the sequencing cohort, patients were older than 65 years with the BMD T score < − 2.5 and concurrent fractures. Healthy volunteers were over 65 years of age with the BMD T score > − 1 and no bone-related diseases (HIV and HBV infections were excluded). The validation cohort was composed of elderly women over 65 years of age from the natural community of Lu'an City, China (26 PMOP patients and 21 non-osteoporosis controls). This study was approved by the Ethics Committee of Lu’an People’s Hospital on December 28, 2020, with the approval number 2020LL-KY042. All patients and volunteers signed informed consent.

Peripheral blood samples from fasting individuals were collected in EDTA tubes following a regular venipuncture procedure. After centrifugation at 3000×g for 15 min at 4 °C, the plasma was aspirated and stored at − 80 °C before use. Dual-energy X-ray absorptiometry (GE Healthcare, USA) was performed to evaluate the bone mineral density (BMD) for subjects, and T score < − 2.5 reflected osteoporosis according to World Health Organization criteria [4].

The exosomes were isolated by SEC (size exclusion chromatography) methods as described previously with minor modifications [23]. Each 1-mL plasma was filtered by 0.8-μm PES material filter membrane (produced by Jinteng) and was diluted with 1.5-fold PBS, and it was further purified using Exosupur® columns (Echobiotech, China) according to the manufacturer’s instructions. Then, the exosomal fractions were concentrated to 200 μL by 100 kDa molecular weight cut-off Amicon® Ultra spin filters (Merck, Germany).

Nano-flow cytometer

The size distribution and particle concentration of EVs were measured by using the nano-flow cytometer (N30E Nanoflow Analyzer, NanoFCM). Briefly, the side scatters intensity (SSI) was detected by the loading of the standard polystyrene nanoparticles (250 nm) to the nano-flow cytometer. Next, the isolated exosomes sample diluted with PBS (according to BCA Protein Assay results diluted the exosomes to 1–10 ng/μL) was loaded into the nano-flow to measure the SSI. Finally, the concentration of exosomes was calculated according to the ratio of SSI to particle concentration in the standard polystyrene nanoparticles. The size distribution was calculated according to the standard curve, which was generated by standard silica nanoparticles with mixed size (68 nm, 91 nm, 113 nm, 155 nm).

Transmission electron microscopy (TEM)

Exosomes solution with 10 µL was placed on a copper mesh at room temperature for 10 min. After washing with sterile distilled water, the exosome was contrasted by uranyl acetate solution for 1 min. The sample was then dried for 2 min under incandescent light. The copper mesh was observed and photographed under a transmission electron microscope (H-7650, Hitachi Ltd., Tokyo, Japan).

Western blot analysis

The exosomes supernatant was denatured in 5 × sodium dodecyl sulfonate (SDS) buffer and subjected to western blot analysis (10% SDS–polyacrylamide gel electrophoresis; 50 µg protein/lane) using rabbit polyclonal antibody CD9 (60,232-I-Ig, Proteintech, Rosemont, IL, 1.50 mg/mL), HSP70 (ab181606, Abcam, USA, 1.12 mg/mL), Alix (ab186429, Abcam, USA, 1.25 mg/mL), and Calnexin (ab22595, Abcam, USA, 1 mg/mL). The PVDF membranes were used and blocked with 5% skim milk powder. The antibodies and secondary antibodies were diluted with NCM Universal Antibody Diluent (Cat.No: WB500D) at 1:1000 and 1:5000, respectively. The secondary antibodies were the rabbit second antibody from Proteintech (Cat.No. SA00001-2).

Exosomal RNA extraction, library preparation, and Sequencing

Total RNA was extracted from plasma exosomes using miRNeasy Serum/Plasma Advanced Kit (Qiagen, cat. No. 217204) according to the kit instruction. Exosomes (200 µL) are mixed with 60 µL Buffer RPL to release and stabilize RNA. The sample is then mixed with 20 µL Buffer RPP and centrifuged to precipitate proteins. Isopropanol is added to the supernatant to provide the appropriate conditions for RNA molecules to bind to the silica membrane. The sample is then applied to the RNeasy UCP MinElute spin column, where RNA binds to the membrane and other contaminants are washed away in subsequent wash steps. In the final step, the total RNA is eluted using RNase-free water.

Small RNA libraries were generated using the QIAseq miRNA Library Kit (Qiagen, Frederick, MD) following the manufacturer’s recommendations. In an unbiased reaction, adapters are ligated sequentially to the 3′ and 5′ ends of 5 µL total RNA on ice; the ligated RNA was reverse-transcribed into cDNA using a reverse transcription (RT) primer with Unique Molecular Index (UMI). Subsequently, the cDNA was cleaned with QIAseq Beads, and Library amplification was performed with a universal forward primer and indexing reverse primers. Then, the library quality was assessed on the Agilent Bioanalyzer 2100. Finally, the libraries were sequenced on an Illumina NovaSeq 6000 platform with 2 × 150 bp paired-end sequencing in an S4 flow cell at EchoBiotech Co. Ltd., Beijing, China.

Quantification and differential expression analysis of exosomal miRNA

The Bowtie software was used to align Clean Reads with the Silva database, GtRNAdb database, Rfam database, and Repbase database, respectively. Then, the repeats and ncRNAs, such as ribosomal RNA (rRNA), transfer RNA (tRNA), small nuclear RNA (snRNA), and small nucleolar RNA (snoRNA), were filtered.

The remaining reads were used to detect miRNAs by comparing them with known miRNAs from miRbase and Human Genome (GRCh38), respectively. The expression matrix of quantified UMI counts of miRNAs was normalized to transcript per million (TPM). Differential expression analysis of two groups was performed using the edgeR package (3.12.1) [24]. p < 0.05 was used as the parameter to select differentially expressed miRNAs (DEMs).

GO and KEGG pathway enrichment analysis

The target genes of DEMs were predicted by miRanda and RNAhybrid [25, 26]. Then, the target genes were aligned to Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database and confirmed their potential biological functions. KEGG pathway enrichment was analyzed using a Python program KOBAS [27]. GOseq R package based on Wallenius non-central hypergeometric distribution was used for GO enrichment analysis.

Verification of candidate miRNAs with qRT-PCR

The method of exosomal RNA extraction is consistent with the above section. The total RNA was then reverse-transcribed to synthesize cDNA using PrimeScript™ RT reagent Kit (Perfect Real Time) (TAKARA, RR037A). The abundance of target gene expression was detected by the TaqMan® probe using real-time qPCR (TaKaRa Ex Taq Hot Start Version, Takara, RR006A)). The 2 µL of cDNA was used as the template for each PCR reaction. Three technical replicates were performed for each qPCR reaction. The sequence of primers and probes are shown in Additional file 1. U6 was used as an internal reference, and the relative expression level of miRNA was calculated with 2−ΔΔct.

Statistical analysis

All statistical analyses were calculated by GraphPad Prism version 9.0. For pair-wise comparisons, parametric t test and nonparametric Mann–Whitney test were applied based on the data distribution (the number of pregnancies was tested by a nonparametric test); p < 0.05 was considered statistically significant. The receiver operating characteristic (ROC) curve was used to evaluate the diagnostic efficiencies of the exosomal miRNA, and logistic regression was applied to determine the diagnostic performance of the combinations. The Pearson correlation coefficient was calculated for correlation analysis.

Results

Characterization of plasma exosomes

In this study, the sequencing cohort included 12 patients with osteoporosis and accompanying fractures, while the control group consisted of 12 healthy individuals of the same age with normal bone density. The validation cohort was composed of elderly women over 65 years of age from the natural community, which consisted of 26 PMOP patients and 21 non-osteoporosis controls. The clinical information on age, height, BMD, weight, number of pregnancies, and menopause age is summarized in Table 1. We found significant differences in the number of pregnancies in both cohorts (p = 0.0013 in the sequencing set, p = 0.016 in the validation set, (Fig. 1A), and the detailed data were listed in Additional file 2), suggesting that the number of pregnancies may be a risk factor for PMOP.

Table 1 Clinical characteristics of enrolled samples
Fig. 1
figure 1

Identification of Plasma exosomes by SEC. A The comparison of the number of pregnancies between PMOP and control group in the sequencing set (p = 0.0013) and validation set (p = 0.0433), the number of pregnancies was tested by nonparametric Mann–Whitney test. B TEM image of the plasma exosomes which were cup-shaped vesicles, Scale bar on the left represents 200 nm. C NanoFCM results showed that plasma exosomes purified by SEC were enriched at 50–100 nm. D Western Blot for plasma exosomes, exosome-specific marker CD9, HSP70, and Alix were detected in plasma exosomes and the negative marker Calnexin was absent in the exosomes sample. The cell lysate here is derived from HEK293 cells, which was taken as the positive control

The exosomes were isolated from plasma by size exclusion chromatography, and the morphology and size distribution of exosomes were characterized by TEM and NanoFCM as recommended in MISEV2018 guidelines [6]. TEM and NanoFCM results showed that the plasma exosomes were cup-shaped with a size range between 50 and 150 nm (Fig. 1B, C). Exosome-specific protein markers CD9, HSP70, and Alix were detected in plasma exosomes by western blot, while the negative marker Calnexin was not present in the exosomes sample (Fig. 1D). These characterization data indicate that our isolated plasma exosomes are typical, laying the foundation for their further research.

miRNA sequencing and differential miRNA identification in PMOP plasma exosomes

To study the miRNA profiles of PMOP plasma exosomes, we generated miRNA sequencing of the sequencing cohorts (12 PMOP patients, 12 control). After removing the low-quality reads and the sequences smaller than 18 nt or longer than 32 nt, 24 samples had an average of 9.15 M clean reads with each sample from 6.21 M to 14.82 M (Additional file 3: Fig. S1), which were higher than the minimum saturated reads (2.50 M) of plasma exosome miRNA sequencing reported in the previous study [28]. There was no significant difference in clean reads between PMOP and the control group (p = 0.7179). A total of 984 known miRNAs were identified in all samples. Among them, the average number of miRNAs identified in the PMOP group was 376, and the control group was 412 with no significant difference between the two groups (p = 0.1841) (Fig. 2A).

Fig. 2
figure 2

miRNA sequencing and differential miRNA identification in PMOP plasma exosomes. A Identified miRNA number in PMOP and control group, the comparisons of the two groups were based on an unpaired Student’s t test, ns indicates that there is no significant difference between the two groups (p > 0.05). B Differentially expressed miRNA volcano plot, where purple dots represent miRNAs upregulated in PMOP over control and blue dots represent down-regulated miRNAs. C Venn plots for comparing upregulated miRNAs to three other published datasets. D Venn plots for comparing down-regulated miRNAs to three other published datasets

The results of differential expression miRNA (DEMs) analysis showed that 17 miRNAs were significantly upregulated, and ten miRNAs were significantly down-regulated in PMOP patients than the control group (Fig. 2B, Table 2). Three published DEM datasets were incorporated for cross-validation (Additional file 4), which examined the exosomal miRNAs expression differences in osteoporosis patients by miRNA sequencing and used different sample sources or isolation methods: 1. Serum exosomes purified by exoEasy [29]; 2. Plasma exosomes isolated by exoEasy [19]; 3. Serum exosomes were separated by UC [20]. Five of the 17 upregulated miRNAs appeared in the other three datasets, including miR-34a-5p, miR-149-5p, miR-320d, miR-455-5p, and miR-7706 (Fig. 2C). And one of the ten downregulated miRNAs intersects with the three datasets, which is miR-451a (Fig. 2D). Additionally, eight upregulated miRNAs and one down-regulated miRNA appeared twice in the other three datasets, and one upregulated miRNA (miR-9-5p) appeared in all three datasets. Approximately, 22% (6/27) of the DEMs were reported in previous studies, indicating the reliability of our study. It also clarifies that more new PMOP-related exosomal miRNAs have been discovered in this study.

Table 2 Differential expression exosomal miRNA between PMOP with control group

Based on bioinformatics database analysis, these 27 DEMs may target a total of 919 coding genes (Additional file 5). The top 10 significant pathways in the KEGG enrichment analysis of target genes are found in Fig. 3A, and the most relevant pathways were protein digestion and absorption, dilated cardiomyopathy, endocrine resistance, and insulin secretion. The biological process (BP), molecular function (MF), and cellular component (CC) in the GO enrichment analysis are shown in Fig. 3B–D. Cellular response to stimulus, extracellular matrix organization and regulation of cell communication in BP, ion binding, anion binding, and protein-containing complex binding in MF, collagen trimer, extracellular matrix component, and cytoskeletal part in CC were the most significant enrichment GO terms.

Fig. 3
figure 3

KEGG and GO enrichment analysis of target genes for differential exosomal miRNAs in sequencing data. A A bubble plot of the top 10 enriched KEGG pathways. B A bubble plot of the top 10 enriched GOs (Biological Process). C A bubble plot of the top 10 enriched GOs (Molecular Function). D A bubble plot of the top 10 enriched GOs (Cellular Component)

Validation of PMOP-related exosomal miRNA

Our 27 DEMs and ten miRNAs that intersect in the other three datasets in the above section made a total of 37 miRNAs, which were then filtered with the median TPM > 50 both in the PMOP and control group, leaving 11 candidate miRNAs for subsequent validation. Considering that miR-9-5p was present in three datasets, it was also included in the final validation (Table 3).

Table 3 Candidate exosomal miRNAs for RT-qPCR verification

To confirm these 12 PMOP-associated exosomal miRNAs, the expression of these miRNAs was measured by RT-qPCR in an independent cohort of 47 individuals (26 PMOP samples and 21 control samples). MiR-196-5p, miR-224-5p, miR320d, miR-34a-5p, miR-9-5p, and miR-98-5p were shown to be significantly upregulated in PMOP than control groups, and all of them were consistent with our sequencing data or previously reported (miR-9-5p) (Fig. 4, Table 3). There were no significant differences of miR-152-3p, miR-200b-3p, miR320c, miR-425-5p, miR-451a, and miR-484 in the validation cohorts (Table 3, Additional file 3: Fig. S2).

Fig. 4
figure 4

Relative expression levels (calculated by 2−ΔΔct) of plasma exosomal miRNA candidates in the validation set. The comparisons of the two groups were based on an unpaired Student’s t-test, with p < 0.05 as the significance threshold. A miR-196a-5p, B miR-224-5p, (C) miR-320d, D miR-34-5p, E miR-9-5p, F miR-98-5p. *means p < 0.05

Evaluation of the diagnostic efficacy of the validated exosomal miRNA for PMOP

We assessed the diagnostic efficacy of each validated exosomal miRNA by simulating the AUC curve. The results were described in (Fig. 5A, Table 4). Among the six miRNAs, miR-9-5p and miR-320d showed the highest diagnostic performance in discriminating PMOP patients and control person, with AUC at 0.695 (95% CI 0.5443–0.8458, sensitivity = 73.08%, specificity = 61.90%) and AUC at 0.695 (95% CI 0.5426–0.8476, sensitivity = 76.92%, specificity = 57.14%), followed by miR-224-5p at AUC = 0.685 (95% CI 0.5302–0.8397, sensitivity = 84.62%, specificity = 47.62%), miR-34a-5p at AUC = 0.681 (95% CI 0.5245–0.8381, sensitivity = 80.77%, specificity = 52.38%) and miR-98-5p at AUC = 0.680 (95% CI 0.5245–0.8345, sensitivity = 100.00%, specificity = 28.57%).

Fig. 5
figure 5

Receiver operating characteristic (ROC) curves for the validated exosomal miRNAs and their combinations. A ROC curve of every single exosomal miRNA. B ROC curve for combinations of exosomal miRNAs. C ROC curve for combinations of the number of pregnancies and exosomal miRNAs. D The Pearson correlation coefficient between the number of pregnancies and the relative level of exosomal miRNA expression

Table 4 Diagnostic efficiencies of plasma exosomal miRNA

We further explored whether multiple exosomal miRNA combinations could be used to improve the diagnostic capacity of discriminating PMOP patients. When two exosomal miRNAs were combined, miR-34a-5p + miR-9-5p (AUC = 0.725) and mi-34a-5p + miR-98-5p (AUC = 0.723) were the two best combinations. MiR-34a-5p + miR-9-5p + miR-98-5p (AUC = 0.734) and miR-320d + miR-34a-5p + miR-98-5p (AUC = 0.729) performed best when the three exosomal miRNAs were combined (Fig. 5B, Table 4). The miR-320d + miR-34a-5p + miR-98-5p combination exhibited high sensitivity (96.15%) with a specificity of 38.10%. The addition of more than three exosomal miRNAs does not significantly improve the performance of PMOP detection.

In the section above, we found that the number of pregnancies (PNs) is a risk factor for PMOP, so we combined it with validated exosomal miRNAs. Interestingly, the number of pregnancies could effectively enhance the efficacy of exosomal miRNA for PMOP diagnosis (Fig. 5C, Table 4). First, the number of pregnancies as a single predictor had an AUC of 0.699. When PNs was combined with a single exosomal miRNA, miR-320d is the best choice, with an AUC of 0.751; when it was combined with two exosomal miRNAs, miR-34a-5p + miR-9-5p (AUC = 0.762, sensitivity = 84.62%, specificity = 71.43%) and 320d + miR-9-5p (AUC = 0.760, sensitivity = 92.31%, specificity = 52.38%) show better performance. The best combination of three exosomal miRNAs and PNs is miR-34a-5p + miR-9-5p + miR-98-5p (AUC = 0.764), which was consistent with the exosomal miRNAs training results alone. Pearson correlation analysis of the PNs with exosomal miRNA expression levels showed that the PNs correlated with exosomal miRNA expression was low (< 0.6) (Fig. 5D), indicating that the number of pregnancies (PNs) is an independent predictor that can improve the performance of plasma exosomal miRNAs in diagnosing PMOP.

Discussion

Here, we identified 27 PMOP-related plasma exosomal miRNAs by next-generation sequencing, and 12 candidate exosomal miRNAs were selected for qRT-PCR validation in another independent cohort, of which six exosomal miRNAs (miR-196-5p, miR-224-5p, miR320d, miR-34a-5p, miR-9-5p, and miR-98-5p) were confirmed to be significantly associated with PMOP. At the same time, we evaluated the performance of exosomal miRNAs in the diagnosis of PMOP for the first time, and the combination of miR-34a-5p + miR-9-5p + miR-98-5p had a better predictive performance with an AUC = 0.734. In addition, the number of pregnancies was found to be an independent risk factor that can improve the performance of exosomal miRNAs in diagnosing PMOP, increasing AUC from 0.734 to 0.764.

In this study, the sequencing cohort included patients with osteoporosis and accompanying fractures, while the control group consisted of healthy individuals of the same age with normal bone density. The validation cohort was composed of elderly women over 65 years of age from the natural community, which represented a real-world situation. The SEC method was used for the isolation of plasma exosomes, which was a more reproduced approach than ultracentrifugation (UC) [30]. Previous studies comparing the applicability of SEC, UC, exoEasy, and exoQuick in plasma exosomal RNomics study showed that SEC had the minimum miRNA binding protein AGO2 and the highest amounts of exosome-specific miRNA, but the lowest nonspecific miRNA [31]. It demonstrated that SEC was a more suitable method for plasma exosomal miRNA research. Then, exosome identification is performed strictly following the guidelines for MISEV2018 [6]. At the sequencing data level, each sample exceeded the minimum number of saturated reads (2.5 M) for plasma exosome miRNA sequencing previously reported in the literature, and the average number of miRNAs identified per sample was close to the number of saturated miRNAs reported in the study of 400 [28]. These details demonstrated that this study is based on high data quality.

A comparison of the sequencing-identified differential expression exosomal miRNAs (DEMs) with the other three DEM datasets revealed that nearly 22% (6/27) of the DEMs were reported in previous studies, indicating the credibility of this study. However, there are a large number of non-intersecting DEMs in different studies. We infer that this is mainly due to different sample collection formats and exosome isolation methods. Studies have reported that when collecting serum, platelets release a large number of exosomes during blood agglutination, thereby diluting the abundance of exosomes in the blood, so plasma can better feedback the real exosomes in human blood than serum collection [32]. Another important reason was that in addition to miRNAs carried in exosomes, there were also a large number of free miRNAs and miRNAs bound to proteins such as AGO2, LDL/HDL, and different isolation methods for exosomes may retain different types of non-exosome miRNAs, thus reflecting different miRNA profile [31, 33]. Therefore, as a blood exosomal biomarker research, it was very important to strictly control the sample collection form and exosome isolation and other test methods; here, we chose plasma collection and SEC isolation exosomes as much as possible to feedback on the real situation of exosomal miRNA in the blood and reduce the non-exosomes-derived miRNA interference.

When we selected miRNAs for validation, we not only considered the results of our sequencing screening, but also included miRNAs that intersected in at least two datasets. At last, six miRNAs were validated by qPCR to prove to be PMOP-associated plasma exosomal miRNAs, including miR-196-5p, miR-224-5p, miR-320d, miR-34a-5p, miR-9-5p, and miR-98-5p. Among them, the upregulation of miR-225-5p has been shown to inhibit osteoblast differentiation by increasing the expression of Pai-1 in the lumbar spine of a rat model of congenital kyphoscoliosis [34]. Muscle-derived miR-34a-5p increased with age in circulating extracellular vesicles and induced senescence of bone marrow stem cells in mice [35]. And the overexpression of miR-9-5p inhibited osteogenic differentiation of bone marrow mesenchymal stem cells by targeting DDX17 [36]. miR-98-5p had been reported to prevent bone regeneration by targeting HMGA2 [37], and the inhibition of miR-98-5p expression could activate PI3K/AKT/GSK3β signaling to promote preosteoblast viability and differentiation by targeting BMP2 [38]. These mechanistic reports support the potential of miR-225-5p, miR-34a-5p, miR-9-5p, and miR-98-5p as diagnostic biomarkers for osteoporosis. The association of miR-196-5p and miR-320d with osteoporosis had not been elaborated, but the MIR196A2 gene expressing miR-196-5p had been reported to be associated with bone mineral density and fracture risk classification in GWAS studies [39]. Even though some of those miRNAs were reported to be associated with bone development and osteoblast differentiation, we were the first to explore and evaluate them as circulating exosomal biomarkers for diagnosing PMOP.

The diagnostic performance evaluation showed that miR-9-5p and miR-320d had the best performance in the single exosomal miRNA model with an AUC = 0.695, miR-34a-5p + miR-9-5p performed best when the two exosomal miRNAs were combined (AUC = 0.725), mi andR-34a-5p + miR-9-5p + miR-98-5p (AUC = 0.734) and miR-320d + miR-34a-5p + miR-98-5p (AUC = 0.729) were the best combinations of the three markers. The miR-320d + miR-34a-5p + miR-98-5p combination exhibited high sensitivity (96.15%) with a specificity of 38.10%. Notably, this study was the first report to evaluate the performance of exosomal miRNAs in diagnosing PMOP. The number of pregnancies was found to be a risk factor for PMOP. Several previous studies have supported the number of pregnancies as a risk factor for PMOP, the conclusion supported by multiple previous studies [21, 22]. However, we were the first study to reveal that the number of pregnancies combined with exosomal miRNA could improve the predictive performance of exosomal miRNA for PMOP.

There were still many limitations to this study that needed to be clarified. Although the size exclusion chromatography method could remove most of the free protein-bound miRNAs in plasma, some miRNAs bound to LDLs would be co-isolated because some LDLs were similar in size to exosomes [40]; therefore, the interference of free miRNAs in plasma on the results could not be completely excluded. The second limitation was that the patients in the sequencing cohort were all patients with osteoporotic fractures, while the validation cohort came from the natural population of the community, which may result in some of the differential exosomal miRNAs screened by sequencing not being validated in the validation cohort; at the same time, the sample size of the validation cohort was relatively small, which may lead to the overfitting of performance evaluation. Finally, although there were differences in AUC values between different diagnostic models, the differences between different models did not pass the Delong test (Additional file 6) [41], which did not indicate significant differences in the diagnostic performance of different models, which might require further validation of larger samples.

Conclusion

In this study, six PMOP-associated plasma exosomal miRNAs (miR-196-5p, miR-224-5p, miR-320d, miR-34a-5p, miR-9-5p, and miR-98-5p) were identified and validated, and their performance in identifying PMOP was initially evaluated. It was also revealed that the number of pregnancies as an independent risk factor for PMOP could enhance the performance of these exosomal miRNAs in diagnosing PMOP. Further validation of large cohorts is required to confirm the performance of these exosomal miRNAs for the diagnosis of PMOP. What tissues and cells release these exosomal miRNAs also need to be further revealed.

Availability of data and materials

The sequencing data of 24 samples had been uploaded to the SRA database (BioProject: PRJNA951192).

References

  1. Camacho PM, Petak SM, Binkley N, Diab DL, Eldeiry LS, Farooki A, Harris ST, Hurley DL, Kelly J, Lewiecki EM, et al. American association of clinical endocrinologists/American college of endocrinology clinical practice guidelines for the diagnosis and treatment of postmenopausal osteoporosis-2020 update. Endocr Pract. 2020;26(Suppl 1):1–46.

    Article  PubMed  Google Scholar 

  2. Rachner TD, Khosla S, Hofbauer LC. Osteoporosis: now and the future. Lancet. 2011;377(9773):1276–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Kanis JA, Melton LJ 3rd, Christiansen C, Johnston CC, Khaltaev N. The diagnosis of osteoporosis. J Bone Miner Res. 1994;9(8):1137–41.

    Article  CAS  PubMed  Google Scholar 

  4. Celi M, Rao C, Scialdoni A, Tempesta V, Gasbarra E, Pistillo P, Tarantino U. Bone mineral density evaluation in osteoporosis: why yes and why not? Aging Clin Exp Res. 2013;25(Suppl 1):S47-49.

    Article  PubMed  Google Scholar 

  5. Warming L, Hassager C, Christiansen C. Changes in bone mineral density with age in men and women: a longitudinal study. Osteoporos Int. 2002;13(2):105–12.

    Article  CAS  PubMed  Google Scholar 

  6. Thery C, Witwer KW, Aikawa E, Alcaraz MJ, Anderson JD, Andriantsitohaina R, Antoniou A, Arab T, Archer F, Atkin-Smith GK, et al. Minimal information for studies of extracellular vesicles 2018 (MISEV2018): a position statement of the international society for extracellular vesicles and update of the MISEV2014 guidelines. J Extracell Vesicles. 2018;7(1):1535750.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, Lotvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007;9(6):654–9.

    Article  CAS  PubMed  Google Scholar 

  8. Chandrasekera D, Katare R. Exosomal microRNAs in diabetic heart disease. Cardiovasc Diabetol. 2022;21(1):122.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Manier S, Liu CJ, Avet-Loiseau H, Park J, Shi J, Campigotto F, Salem KZ, Huynh D, Glavey SV, Rivotto B, et al. Prognostic role of circulating exosomal miRNAs in multiple myeloma. Blood. 2017;129(17):2429–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Min L, Zhu S, Chen L, Liu X, Wei R, Zhao L, Yang Y, Zhang Z, Kong G, Li P, et al. Evaluation of circulating small extracellular vesicles derived miRNAs as biomarkers of early colon cancer: a comparison with plasma total miRNAs. J Extracell Vesicles. 2019;8(1):1643670.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Melo SA, Luecke LB, Kahlert C, Fernandez AF, Gammon ST, Kaye J, LeBleu VS, Mittendorf EA, Weitz J, Rahbari N, et al. Glypican-1 identifies cancer exosomes and detects early pancreatic cancer. Nature. 2015;523(7559):177–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Rastogi S, Sharma V, Bharti PS, Rani K, Modi GP, Nikolajeff F, Kumar S. The evolving landscape of exosomes in neurodegenerative diseases: exosomes characteristics and a promising role in early diagnosis. Int J Mol Sci. 2021;22(1):440.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wei Z, Batagov AO, Schinelli S, Wang J, Wang Y, El Fatimy R, Rabinovsky R, Balaj L, Chen CC, Hochberg F, et al. Coding and noncoding landscape of extracellular RNA released by human glioma stem cells. Nat Commun. 2017;8(1):1145.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Lai G, Zhao R, Zhuang W, Hou Z, Yang Z, He P, Wu J, Sang H. BMSC-derived exosomal miR-27a-3p and miR-196b-5p regulate bone remodeling in ovariectomized rats. PeerJ. 2022;10:e13744.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Li L, Zhou X, Zhang JT, Liu AF, Zhang C, Han JC, Zhang XQ, Wu S, Zhang XY, Lv FQ. Exosomal miR-186 derived from BMSCs promote osteogenesis through hippo signaling pathway in postmenopausal osteoporosis. J Orthop Surg Res. 2021;16(1):23.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Lu GD, Cheng P, Liu T, Wang Z. BMSC-derived exosomal miR-29a promotes angiogenesis and osteogenesis. Front Cell Dev Biol. 2020;8:608521.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Yu L, Hu M, Cui X, Bao D, Luo Z, Li D, Li L, Liu N, Wu Y, Luo X, et al. M1 macrophage-derived exosomes aggravate bone loss in postmenopausal osteoporosis via a microRNA-98/DUSP1/JNK axis. Cell Biol Int. 2021;45(12):2452–63.

    Article  CAS  PubMed  Google Scholar 

  18. Kim DK, Bandara G, Cho YE, Komarow HD, Donahue DR, Karim B, Baek MC, Kim HM, Metcalfe DD, Olivera A. Mastocytosis-derived extracellular vesicles deliver miR-23a and miR-30a into pre-osteoblasts and prevent osteoblastogenesis and bone formation. Nat Commun. 2021;12(1):2527.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Kong D, Chen T, Zheng X, Yang T, Zhang Y, Shao J. Comparative profile of exosomal microRNAs in postmenopausal women with various bone mineral densities by small RNA sequencing. Genomics. 2021;113(3):1514–21.

    Article  CAS  PubMed  Google Scholar 

  20. Shi H, Jiang X, Xu C, Cheng Q. MicroRNAs in serum exosomes as circulating biomarkers for postmenopausal osteoporosis. Front Endocrinol (Lausanne). 2022;13:819056.

    Article  PubMed  Google Scholar 

  21. Kim KJ, An JH, Kim KJ, Yu JH, Kim NH, Yoo HJ, Kim HY, Seo JA, Kim NH, Choi KM, et al. Prevalence of osteoporosis among North Korean women refugees living in South Korea: a comparative cross-sectional study. BMJ Open. 2020;10(6):e036230.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Yan G, Huang Y, Cao H, Wu J, Jiang N, Cao X. Association of breastfeeding and postmenopausal osteoporosis in Chinese women: a community-based retrospective study. BMC Womens Health. 2019;19(1):110.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Boing AN, van der Pol E, Grootemaat AE, Coumans FA, Sturk A, Nieuwland R. Single-step isolation of extracellular vesicles by size-exclusion chromatography. J Extracell Vesicles. 2014;3:23430.

    Article  Google Scholar 

  24. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    Article  CAS  PubMed  Google Scholar 

  25. Betel D, Wilson M, Gabow A, Marks DS, Sander C. The microRN.Aorg resource: targets and expression. Nucleic Acids Res. 2008;36(1):149–53.

    Google Scholar 

  26. Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R. Fast and effective prediction of microRNA/target duplexes. RNA. 2004;10(10):1507–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21(19):3787–93.

    Article  CAS  PubMed  Google Scholar 

  28. Srinivasan S, Yeri A, Cheah PS, Chung A, Danielson K, De Hoff P, Filant J, Laurent CD, Laurent LD, Magee R, et al. Small RNA sequencing across diverse biofluids identifies optimal methods for exRNA isolation. Cell. 2019;177(2):446-462 e416.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Shao JL, Li H, Zhang XR, Zhang X, Li ZZ, Jiao GL, Sun GD. Identification of serum exosomal MicroRNA expression profiling in menopausal females with osteoporosis by high-throughput sequencing. Curr Med Sci. 2020;40(6):1161–9.

    Article  CAS  PubMed  Google Scholar 

  30. Koh YQ, Almughlliq FB, Vaswani K, Peiris HN, Mitchell MD. Exosome enrichment by ultracentrifugation and size exclusion chromatography. Front Biosci (Landmark Ed). 2018;23(5):865–74.

    CAS  PubMed  Google Scholar 

  31. Yang Y, Wang Y, Wei S, Zhou C, Yu J, Wang G, Wang W, Zhao L. Extracellular vesicles isolated by size-exclusion chromatography present suitability for RNomics analysis in plasma. J Transl Med. 2021;19(1):104.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Coumans FAW, Brisson AR, Buzas EI, Dignat-George F, Drees EEE, El-Andaloussi S, Emanueli C, Gasecka A, Hendrix A, Hill AF, et al. Methodological guidelines to study extracellular vesicles. Circ Res. 2017;120(10):1632–48.

    Article  CAS  PubMed  Google Scholar 

  33. Van Deun J, Mestdagh P, Sormunen R, Cocquyt V, Vermaelen K, Vandesompele J, Bracke M, De Wever O, Hendrix A. The impact of disparate isolation methods for extracellular vesicles on downstream RNA profiling. J Extracell Vesicles. 2014;3:24858.

    Article  Google Scholar 

  34. Ishiwata S, Iizuka H, Sonoda H, Tsunoda D, Tajika Y, Chikuda H, Koibuchi N, Shimokawa N. Upregulated miR-224-5p suppresses osteoblast differentiation by increasing the expression of Pai-1 in the lumbar spine of a rat model of congenital kyphoscoliosis. Mol Cell Biochem. 2020;475(1–2):53–62.

    Article  CAS  PubMed  Google Scholar 

  35. Fulzele S, Mendhe B, Khayrullin A, Johnson M, Kaiser H, Liu Y, Isales CM, Hamrick MW. Muscle-derived miR-34a increases with age in circulating extracellular vesicles and induces senescence of bone marrow stem cells. Aging (Albany NY). 2019;11(6):1791–803.

    Article  CAS  PubMed  Google Scholar 

  36. He C, Liu M, Ding Q, Yang F, Xu T. Upregulated miR-9-5p inhibits osteogenic differentiation of bone marrow mesenchymal stem cells under high glucose treatment. J Bone Miner Metab. 2022;40(2):208–19.

    Article  CAS  PubMed  Google Scholar 

  37. Zheng F, Wang F, Xu Z. MicroRNA-98-5p prevents bone regeneration by targeting high mobility group AT-Hook 2. Exp Ther Med. 2019;18(4):2660–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Zheng F, Zhang F, Wang F. Inhibition of miR-98-5p promotes high glucose-induced suppression of preosteoblast proliferation and differentiation via the activation of the PI3K/AKT/GSK3beta signaling pathway by targeting BMP2. Mol Med Rep. 2022;26(3):1–9.

    Article  Google Scholar 

  39. Styrkarsdottir U, Stefansson OA, Gunnarsdottir K, Thorleifsson G, Lund SH, Stefansdottir L, Juliusson K, Agustsdottir AB, Zink F, Halldorsson GH, et al. GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures. Nat Commun. 2019;10(1):2054.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011;13(4):423–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44(3):837–45.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank Prof Liangye Sun and Guanyi Kong for providing guidance.

Funding

This research was supported by the fund of Anhui medical university (2019xkj222).

Author information

Authors and Affiliations

Authors

Contributions

ZS, SJ, CY, XC, JC, J Chen, and YW completed most of the experiments and analyzed the data. XC, JC, CZ, and JX collected the sample and isolated the exosomes. ZS, SJ, and CY wrote the first draft. SS and GT designed the study, proofread, and revised the submission. All authors contributed to the article and approved the submitted version.

Corresponding author

Correspondence to Song Shao.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the ethics committee of Anhui Medical University, and all participants had signed informed consent. All experiments were in line with the principles of The Code of Ethics of the World Medical Association (Declaration of Helsinki) for research experiments involving human subjects.

Competing interests

The authors declare that they have no conflict of interest. Author Jing Chen and Yuan Wang are employed by Echo Biotech Co., Ltd, Beijing, P. R. China. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. The sequence of primers and probes for candidate exosomal miRNAs.

Additional file 2

. The number of pregnancies in both cohorts.

Additional file 3

. The Supplementary Figures of this article. Fig. S1. The clean reads number of every sequencing sample, after removing the low-quality reads and the sequences smaller than 18 nt or longer than 32 nt. Fig. S2. Relative expression levels (calculated by 2-ΔΔct) of plasma exosomal miRNA candidates in the validation set. The comparisons of the two groups were based on an unpaired Student’s t-test, with p < 0.05 as the significance threshold. This figure shows miRNAs without significant differences. (A) miR-152-3p, (B) miR-220b-3p, (C) miR-320c, (D) miR-425-5p, (E) miR-451a, (F) miR-484.

Additional file 4

. The differential expression exosomal miRNAs of three published datasets.

Additional file 5

. The target gene list of differential expression exosomal miRNAs in sequencing data.

Additional file 6

. The p-value of Delong test between different diagnostic models.

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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sun, Z., Shi, J., Yang, C. et al. Identification and evaluation of circulating exosomal miRNAs for the diagnosis of postmenopausal osteoporosis. J Orthop Surg Res 18, 533 (2023). https://doi.org/10.1186/s13018-023-04020-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-023-04020-z

Keywords