| Literature DB >> 36010568 |
Chenglong Wang1, Sabine Stöckl1, Shushan Li1, Marietta Herrmann2, Christoph Lukas1, Yvonne Reinders3, Albert Sickmann3,4,5, Susanne Grässel1.
Abstract
Osteoporosis, or steroid-induced osteonecrosis of the hip, is accompanied by increased bone marrow adipogenesis. Such a disorder of adipogenic/osteogenic differentiation, affecting bone-marrow-derived mesenchymal stem cells (BMSCs), contributes to bone loss during aging. Here, we investigated the effects of extracellular vesicles (EVs) isolated from human (h)BMSCs during different stages of osteogenic differentiation on the osteogenic and adipogenic differentiation capacity of naïve (undifferentiated) hBMSCs. We observed that all EV groups increased viability and proliferation capacity and suppressed the apoptosis of naïve hBMSCs. In particular, EVs derived from hBMSCs at late-stage osteogenic differentiation promoted the osteogenic potential of naïve hBMSCs more effectively than EVs derived from naïve hBMSCs (naïve EVs), as indicated by the increased gene expression of COL1A1 and OPN. In contrast, the adipogenic differentiation capacity of naïve hBMSCs was inhibited by treatment with EVs from osteogenic differentiated hBMSCs. Proteomic analysis revealed that osteogenic EVs and naïve EVs contained distinct protein profiles, with pro-osteogenic and anti-adipogenic proteins encapsulated in osteogenic EVs. We speculate that osteogenic EVs could serve as an intercellular communication system between bone- and bone-marrow adipose tissue, for transporting osteogenic factors and thus favoring pro-osteogenic processes. Our data may support the theory of an endocrine circuit with the skeleton functioning as a ductless gland.Entities:
Keywords: ECM remodeling; adipogenic differentiation; bone regeneration; extracellular vesicles; mesenchymal stem cells; osteogenic differentiation; osteogenic potential
Mesh:
Year: 2022 PMID: 36010568 PMCID: PMC9406723 DOI: 10.3390/cells11162491
Source DB: PubMed Journal: Cells ISSN: 2073-4409 Impact factor: 7.666
Figure 1Overview of experimental set-up; EVs = extracellular vesicles.
Primer sequences for qPCR.
| Gene | Primer Sequences (5′–3′) |
|---|---|
| GAPDH | Fwd, CTGACTTCAACAGCGACACC |
| Rev, CCCTGTTGCTGTAGCCAAAT | |
| TBP | Fwd, TTGTACCGCAGCTGCAAAAT |
| Rev, TATATTCGGCGTTTCGGGCA | |
| ALP | Fwd, CCTCCTCGGAAGACACTCTG |
| Rev, CCTCCTCGGAAGACACTCTG | |
| BGLAP | Fwd, GTGCAGAGTCCAGCAAAGGT |
| Rev, TCAGCCAACTCGTCACAGTC | |
| COL1A1 | Fwd, ACGTCCTGGTGAAGTTGGTC |
| Rev, ACCAGGGAAGCCTCTCTCTC | |
| RUNX2 | Fwd, CGGAATGCCTCTGCTGTTATG |
| Rev, GCTTCTGTCTGTGCCTTCTG | |
| OPN | Fwd, TGAAACGAGTCAGCTGGATG |
| Rev, TGAAATTCATGGCTGTGGAA | |
| C/EBPβ | Fwd, AACCTGGAGACGCAGCACAA |
| Rev, GAACAAGTTCCGCAGGGTGC | |
| PPARγ | Fwd, GACCAAAGCAAAGGCGAGGG |
| Rev, CCCTGAAAGATGCGGATGGC | |
| ADIPOQ | Fwd, AAGGAGATCCAGGTCTTATTGGTC |
| Rev, CGAATGGGCATGTTGGGGAT |
Figure 2Molecular characterization of EVs. (A–C) Particle size distribution of PBS, naïve hBMSCs-EVs and osteogenic EVs was measured by NTA analysis. n = 4. (D,E) Quantitative comparison among PBS, naïve hBMSCs-EVs and osteogenic EVs in count and size measured by NTA analysis. n = 4. * p < 0.05; two-tailed Mann–Whitney U-test was used to compare differences between different groups. (F) Uptake of EVs by hBMSCs. PKH26-labeled EVs_D28 (red) were internalized by hBMSCs and visualized with fluorescence microscopy. Cell nuclei were stained with DAPI. Scale bar: 100 μm. n = 2.
Figure 3Effects of EVs derived from different stages of BMSC osteogenesis on naïve hBMSC proliferation, viability and apoptosis. (A) For proliferation analysis, naïve hBMSCs were stimulated for 3 days with the different EV groups or PBS (no EVs) in medium containing 1% EV-depleted FCS; n = 6–7. Analysis of viability (B) and apoptosis (C) of naïve hBMSCs after stimulation for 3 days with the different EV groups or PBS (no EVs) in medium containing 1% EVs-depleted FCS; n = 6–7. Results were calculated as percentage to the control group (no EVs, shown by the dotted line); * p < 0.05; Wilcoxon signed-rank test was used when no EVs group was set to 100%, and two-tailed Mann–Whitney U-test was used to compare differences between different EV-treated groups.
Figure 4Evaluation of osteogenic differentiation of naïve hBMSCs after EV treatment. (A) Alizarin Red staining of hBMSCs after 21 days of osteogenic differentiation and treatment with the different EVs groups. Microscopic view (scale bar: 0.5 mm) and macroscopic view (scale bar: 1 cm); n = 6–8. (B) Quantification of Alizarin Red staining; n = 6–8. (C) Alkaline Phosphatase (ALP) activity of hBMSC after 14 days of osteogenic differentiation and treatment with the different EV groups; n = 7. Results were calculated as percentage to the control group (no EVs, shown by the dotted line); * p < 0.05, ** p < 0.01; Wilcoxon signed-rank test was used when no EVs group was set to 100%, and two-tailed Mann–Whitney U-test was used to compare differences between different EV-treated groups.
Figure 5Evaluation of osteogenic marker gene expression in hBMSCs after EV treatment. (A–E) Gene expression level of the osteogenic marker genes (OPN, BGLAP, ALP, RUNX2, and COL1A1) were analyzed after 14 days of osteogenic differentiation of hBMSCs and simultaneous stimulation with the different EV groups. Results were calculated as percentage to the control group (no EVs, shown by the dotted line); * p < 0,05; Wilcoxon signed-rank test was used when no EVs group was set to 100%, and two-tailed Mann–Whitney U-test was used to compare differences between different EV-treated groups; n = 5–7.
Figure 6Evaluation of adipogenic differentiation of naïve hBMSCs after EV treatment. (A) Microscopic view of Oil red O staining of hBMSCs after 21 days of adipogenic differentiation of hBMSCs and treatment with the different EVs groups; n = 8. (B) Quantification of Oil red O staining; n = 8. (C–E) Gene expression level of adipogenic markers such as ADIPOQ, C/EBPα, and PPARγ after 14 days of adipogenic differentiation and simultaneous stimulation with the different EVs groups; n = 6–8. Results were calculated as percentage to the control group (no EVs, shown by the dotted line); * p < 0.05, ** p < 0.01; Wilcoxon signed-rank test was used when no EVs group was set to 100%, and two-tailed Mann–Whitney U-test was used to compare differences between different EV-treated groups. Scale bar: 100 μm.
Figure 7Venn diagram of total proteins and heatmap of significantly different proteins. (A) Venn diagram showing the distinct profiles of total proteins in the EVs_D28-35 and EVs_D0; n = 3. (B) Heatmap showing the distinct significant protein profiles of the EVs_D28-35 compared to EVs_D0; n = 3.
The expression/level ratio of selected pro-osteogenic proteins in EVs_D28-35 compared with that in EVs_D0.
| UniProt | Protein Names | Gene Names | EVs_D28-35 |
|---|---|---|---|
| P23142 | Fibulin-1 | FBLN1 | 100 |
| P51888 | Prolargin | PRELP | 100 |
| P08493 | Matrix Gla protein | MGP | 63.867 |
| P01034 | Cystatin-C | CST3 | 7.567 |
| Q8IX30 | Signal peptide, CUB and EGF-like domain-containing protein 3 | SCUBE3 | 7.149 |
| P11362 | Fibroblast growth factor receptor 1 | FGFR1 | 6.795 |
| O76076 | WNT1-inducible-signaling pathway protein 2 | CCN5 | 6.312 |
| P02751 | IGFBP4 | FN1 | 6.29 |
| Q96CG8 | Collagen triple helix repeat-containing protein 1 | CTHRC1 | 5.782 |
| P22692 | Insulin-like growth factor-binding protein 4 | IGFBP4 | 4.895 |
| P24592 | Insulin-like growth factor-binding protein 6 | IGFBP6 | 4.756 |
| P07585 | Decorin | DCN | 4.365 |
| P05452 | Tetranectin | CLEC3B | 4.298 |
| P08253 | 72 kDa type IV collagenase | MMP2 | 4.06 |
| P01344 | Insulin-like growth factor II | IGF2 | 4.017 |
| P18065 | Insulin-like growth factor-binding protein 2 | IGFBP2 | 3.611 |
| P17936 | Insulin-like growth factor-binding protein 3 | IGFBP3 | 3.582 |
| P09429 | High mobility group protein B1 | HMGB1 | 3.223 |
| P12109 | Collagen alpha-1(VI) chain | COL6A1 | 3.163 |
| P13611 | Versican core protein | VCAN | 3.093 |
| P49747 | Cartilage oligomeric matrix protein | COMP | 2.569 |
| Q4LDE5 | Sushi, von Willebrand factor type A, EGF and pentraxin domain-containing protein 1 | SVEP1 | 2.515 |
The expression/level ratio of selected anti-adipogenic proteins in EVs_D28-35 compared with that in EVs_D0.
| UniProt | Protein Names | Gene Names | EVs_D28-35 |
|---|---|---|---|
| Q8IUX7 | Adipocyte enhancer-binding protein 1 | AEBP1 | 100 |
| Q76M96 | Coiled-coil domain-containing protein 80 | CCDC80 | 100 |
| P05090 | Apolipoprotein D | APOD | 8.375 |
| Q9NRY6 | Phospholipid scramblase 3 | PLSCR3 | 2.635 |
The expression/level ratio of selected proteins involved in ECM remodeling in EVs_D28-35 compared with that in EVs_D0.
| UniProt | Protein Names | Gene Names | EVs_D28-35 |
|---|---|---|---|
| P55268 | Laminin subunit beta-2 | LAMB2 | 6.053 |
| Q9Y6C2 | EMILIN-1 | EMILIN1 | 5.021 |
| P07942 | Laminin subunit beta-1 | LAMB1 | 4.939 |
| P98095 | Fibulin-2 | FBLN2 | 4.838 |
| Q14112 | Nidogen-2 | NID2 | 3.652 |
| Q16363 | Laminin subunit alpha-4 | LAMA4 | 3.344 |
| Q15113 | ProcollagenC-endopeptidase enhancer 1 | PCOLCE | 2.897 |
| P12111 | Collagen alpha-3(VI) chain | COL6A3 | 2.492 |
GO terms related to cell metabolism (proliferation, apoptosis, viability) in the GO biological process analysis of upregulated proteins in EVs_D28-35 compared with that in EVs_D0.
| GO-Term ID | Term Description | Observed Gene Count | False Discovery Rate ( | Candidate Gene Name |
|---|---|---|---|---|
| GO:0008284 | Positive regulation of cell population proliferation | 22 | 8.71 × 10−6 | TIMP1, MMP2, SFRP1, COMP, LAMB1, IGFBP2, LAMC1, SDCBP, PDCD6, HMGB2, NPM1, AQP1, CTHRC1, FBLN1, GAS6, HMGB1, FN1, HTRA1, CST3, IGF2, FGFR1, NACA |
| GO:0008285 | Negative regulation of cell population proliferation | 13 | 0.0125 | SFRP1, SERPINF1, TIMP2, NPM1, IGFBP6, PODN, FBLN1, CAV1, APOD, SRPX, IGFBP3, CD9, B2M |
| GO:0060548 | Negative regulation of cell death | 17 | 0.0047 | TIMP1, SFRP1, COMP, SERPINF1, HMGB2, NPM1, HSP90B1, AXL, AQP1, GAS6, CAV1, FN1, SET, WISP2, HSPG2, CST3, NACA |
| GO:0031324 | Negative regulation of cellular metabolic processes | 28 | 0.0293 | TIMP1, SFRP1, AEBP1, SERPINF1, SDCBP, SLC4A1, TIMP2, APLP2, COL6A3, HMGB2, NPM1, AQP1, NCL, FBLN1, GAS6, HIST2H2AC, CAV1, APOD, HMGB1, SET, EMILIN1, IGFBP3, PROS1, MVP, CST3, IGF2, FGFR1, NACA |
| GO:0043085 | Positive regulation of catalytic activity | 21 | 0.0070 | SFRP1, PCOLCE, CHI3L1, TIMP2, PDCD6, ACE, HMGB2, AXL, IGFBP6, FBLN1, GAS6, HSP90AA1, CAV1, HMGB1, |
| GO:0043086 | Negative regulation of catalytic activity | 15 | 0.0051 | TIMP1, SFRP1, SERPINF1, |
Figure 8Functional enrichment analysis of significantly regulated proteins. Gene ontology (GO) analysis for the top 20 terms of upregulated cellular components (A), molecular functions (B), and biological processes (C). Kyoto Encyclopedia of Genes and Genomes enrichment analysis for significant proteins were clustered, and the top 20 pathways are shown; n = 3. Note: −log10 (p = 0.05) = 1.30 (the bottom horizontal axis number). −log10 (p = 0.01) = 2. −log10 (p = 0.001) = 3.
Figure 9Summary Osteogenic primed EVs positively regulate osteogenic differentiation and negatively regulate adipogenic differentiation capacity of naïve hBMSCs.