Association between in vivo bone formation and ex vivo migratory capacity of human bone marrow stromal cells
Stem Cell Research & Therapy volume 6, Article number: 196 (2015)
There is a clinical need for developing systemic transplantation protocols for use of human skeletal stem cells (also known bone marrow stromal stem cells) (hBMSC) in tissue regeneration. In systemic transplantation studies, only a limited number of hBMSC home to injured tissues suggesting that only a subpopulation of hBMSC possesses “homing” capacity. Thus, we tested the hypothesis that a subpopulation of hBMSC defined by ability to form heterotopic bone in vivo, is capable of homing to injured bone.
We tested ex vivo and in vivo homing capacity of a number of clonal cell populations derived from telomerized hBMSC (hBMSC-TERT) with variable ability to form heterotopic bone when implanted subcutaneously in immune deficient mice. In vitro transwell migration assay was used and the in vivo homing ability of transplanted hBMSC to bone fractures in mice was visualized by bioluminescence imaging (BLI). In order to identify the molecular phenotype associated with enhanced migration, we carried out comparative DNA microarray analysis of gene expression of hBMSC-derived high bone forming (HBF) clones versus low bone forming (LBF) clones.
HBF clones were exhibited higher ex vivo transwell migration and following intravenous injection, better in vivo homing ability to bone fracture when compared to LBF clones. Comparative microarray analysis of HBF versus LBF clones identified enrichment of gene categories of chemo-attraction, adhesion and migration associated genes. Among these, platelet-derived growth factor receptor (PDGFR) α and β were highly expressed in HBF clones. Follow up studies showed that the chemoattractant effects of PDGF in vitro was more enhanced in HBF compared to LBF clones and this effect was reduced in presence of a PDGFRβ-specific inhibitor: SU-16 f. Also, PDGF exerted greater chemoattractant effect on PDGFRβ+ cells sorted from LBF clones compared to PDGFRβ- cells.
Our data demonstrate phenotypic and molecular association between in vivo bone forming ability and migratory capacity of hBMSC. PDGFRβ can be used as a potential marker for the prospective selection of hBMSC populations with high migration and bone formation capacities suitable for clinical trials for enhancing bone regeneration.
Human skeletal stem cells (also known as human bone marrow-derived stromal cells (hBMSC)) are adult multipotent stem cells located in the bone marrow perivascular niche and are recruited to bone formation sites during bone remodeling . During recent years, hBMSC have been tested in a number of clinical trials for their ability to enhance tissue repair including tissue regeneration where hBMSC were injected locally at the sites of tissue injury; for example, bone fracture [2–4] or ischemic myocardium [5–8]. However, systemic intravenous infusion is more suitable for clinical cell transplantation and is employed for hematopoietic stem cell (HSC) transplantation with success and where HSCs, following homing from systemic circulation to bone marrow, engraft and initiate hematopoiesis .
Several studies have demonstrated that systemically injected bone marrow-derived stromal cells (BMSC) can home to damaged tissues in animal models of brain injury , skeletal disorders [11–13], and acute radiation syndrome [14, 15]. However, the number of BMSC that home and engraft in injured tissues is usually small and most of the infused BMSC get entrapped in the lungs [16, 17]. The explanation for these phenomena is still missing because the mechanisms governing migration of BMSC to injured tissues are poorly understood .
Cultured hBMSC are a heterogeneous population of cells that when analyzed at a clonal level exhibit variations in cell morphology, proliferation, and differentiation capacity [19, 20]. Recently, we have also demonstrated that clonal heterogeneity of the hBMSC population reflects functional heterogeneity with respect to cell capacity for osteoblast adipocyte differentiation or immune functions [21, 22].
Here we hypothesized the existence of clonal heterogeneity in the ability of hBMSC to home to injured tissues (e.g., bone fractures) and that hBMSC with good bone-forming capacity will be more efficient at homing to bone fracture sites. To test this hypothesis, we examined the ex vivo and in vivo migratory capacity of a number of clonal cell populations isolated from telomerized hBMSC that exhibit variation in their ability to form heterotopic bone when implanted in vivo . Our results demonstrate that there is phenotypic association between the in vivo bone formation and migratory capacity to bone fracture sites, and furthermore identified platelet-derived growth factor receptor (PDGFR)α and PDGFRβ as potential markers for the hBMSC population with enhanced migratory function.
Human mesenchymal stem cell culture
As a model for primary hBMSC, we employed our well-characterized telomerized hBMSC-TERT cell line, established by ectopic expression of the catalytic subunit of human telomerase as described previously . The hBMSC-TERT cells exhibit a stable cellular and molecular phenotype during in vitro culture similar to that of primary hBMSC . The derivation and characterization of hBMSC-TERT+Bone, hBMSC-TERT–Bone, and high bone-forming (HBF) and low bone-forming (LBF) single-cell clones have been described previously by our group . In brief, the hBMSC-TERT+Bone subpopulation was derived from early-passage hBMSC-TERT cells (population doubling level 77), and showed high capacity for in vivo heterotopic bone formation, while the hBMSC-TERT–Bone subpopulation was derived from hBMSC-TERT cells (population doubling level 233), and showed LBF capacity. Both the HBF and LBF single-cell clones were derived from hBMSC-TERT+Bone cells by the limiting dilution method.
The cells were cultured in a standard growth medium containing minimal essential medium (MEM) (Gibco, Invitrogen, Herlev, Denmark) supplemented with 10 % fetal calf serum (FCS) (Biochrom, Berlin, Germany) and 1 % penicillin/streptomycin (Gibco, Invitrogen, Herlev, Denmark) at 37 °C in a humidified atmosphere containing 5 % CO2. The medium was changed every third day until cells became 90 % confluent.
Studying cell spreading and focal adhesion formations
Cell culture plates were coated with fibronectin (10 μg/ml) in phosphate-buffered saline (PBS) for 2 hours at 37 °C, rinsed twice with PBS, and blocked with 1 % bovine serum albumin (BSA) for 1 hour. Cells were trypsinized, washed twice with MEM, and resuspended in serum-free medium for 1 hour at room temperature. Cell were then replated onto fibronectin-coated plates in standard culture medium supplemented with platelet-derived growth factor (PDGF)-BB (100 ng/ml) for 30 minutes at 37 °C. Cells were fixed in 4 % paraformaldehyde for 10 minutes, washed with PBS, and stained for F-actin with Phalloidin–fluorescein isothiocyanate (FITC) (Sigma, Brøndby, Denmark) and for focal adhesion with mouse monoclonal anti-Vinculin antibody (Sigma, Brøndby, Denmark) using Alexa Fluor® 488-conjugated rabbit anti-mouse IgG (H + L) as secondary antibody (Cell Signaling). Fluorescent staining was visualized by the Operetta® High Content Imaging system (Perkin Elmer, Rodgau, Germany) at 20× magnification. Fluorescent images were analyzed using Harmony® High Content Imaging Analysis Software (Perkin Elmer, Rodgau, Germany).
Microarray analysis was performed on hBMSC-TERT+Bone, hBMSC-TERT–Bone, and hBMSC-TERT-derived single-cell clones: HBF (three clones; DD8, AD10, BB10) and LBF (three clones; CF1, CB4, CB8). To perform global gene expression analysis, 500 ng quality-checked total RNA per sample in triplicate were amplified and labeled according to the manufacturer’s protocol (Illumina TotalPrep RNA Amplification Kit; Ambion, Austin, TX, USA ). The resulting biotinylated cRNA was purified and hybridized to Illumina HumanRef-8 v3 Expression BeadChips (Illumina, San Diego, CA, USA ) on the Illumina Beadstation 500 platform. This was followed by washing, blocking, staining with streptavidin-Cy3, and quantitative detection of the fluorescent image of the array as specified by the manufacturer.
Raw data were processed using the Gene Expression Module version 1.8.0 provided with the GenomeStudio software (Illumina). This included background subtraction and normalization according to the “rank invariant” algorithm. Genes were considered “expressed” if the corresponding “Detection P-Value” given by the GenomeStudio software was p det <0.01. Differential gene expression analysis was computed using the Illumina Custom Model. p values of differentially expressed genes (“Diff P-Values”) were modified by the Benjamini and Hochberg false discovery rate (FDR) correction algorithm . Genes were considered differentially expressed genes if: they exhibited twofold upregulation or downregulation in average signal intensity; the corresponding FDR-adjusted “Diff P-Value” was p adj <0.05; or the gene of interest was expressed in at least one of the samples under consideration.
Differentially upregulated genes in group 1 (hBMSC-TERT+Bone and three HBF clones) when compared with group 2 (hBMSC-TERT–Bone and three LBF clones) are presented in Table 1.
The microarray data have been deposited in the public repository [GEO:GSE69358].
Total RNA isolation and real-time PCR
RNA was isolated from cultured cells using Trizol® (Bio-Rad, Herlev, Denmark) according to the manufacturer’s protocol. cDNA was synthesized using a commercial revertAid H minus first strand cDNA synthesis kit (Fermentas, Copenhagen, Denmark) according to the instruction manual. Real-time PCR was performed in the iCycler IQ detection system (Bio-Rad, Herlev, Denmark) using SYBR® Green I as described previously . Primer sequences are indicated in Table S1 in Additional file 1. The expression level for each target gene was calculated using the comparative Ct formula 1/(2∆Ct) and data were presented as relative expression to the reference genes (RG) (HPRT1 and UBC). Data were analyzed using Microsoft Excel 2000 to generate relative expression values.
Flow cytometry analysis
Cells were trypsinized and washed with PBS containing 0.5 % BSA (Fluorescence Activated Cell Sorting (FACS) buffer) and 50 μg/ml trypsin inhibitor. Cells were then incubated with primary antibody for PDGFRβ (0.025 mg/ml, PR292a; R&D Systems (Abingdon, UK) diluted 1:50 in FACS buffer for 30 minutes on ice. The cells were washed twice with ice-cold FACS buffer and incubated with FITC-conjugated goat anti-mouse antibody (R&D Systems) for 30 minutes on ice. The cells were analyzed by flow cytometer Cell Lab Quanta TMSC (Beckman Coulter using Kaluza version1.2 (Beckman Coulter, Copenhagen, Denmark).
Magnetic activated cell sorting of hBMSC-TERT–Bone
Cells were rinsed with PBS, lifted with trypsin, and diluted in Magnetic activated cell sorting (MACS) buffer (PBS with 0.5 % (w/v) BSA and 2 mM ethylenediamine tetraacetic acid (EDTA)) before incubation with PDGFRβ antibody (R&D Systems) for 20 minutes at 4 °C. Cells were washed with MACS buffer, incubated with anti-mouse-beads (Miltenyi Biotec, Lund, Sweden) for 15 minutes at 4 °C, and run through a MS column (Miltenyi Biotec, Lund, Sweden) according to the manufacturer’s instructions.
Generating the Luciferase-overexpressing hBMSC-TERT+Bone and hBMSC-TERT–Bone
The hBMSC-TERT+Bone and hBMSC-TERT–Bone cells were transduced with a retroviral vector containing the firefly luciferase gene (Luc2) producing luciferase-containing hBMSC-TERT (hBMSC-TERT-Luc). The luciferase reporter gene luc2 (Photinus pyralis) was subcloned into the retroviral vector pBABE-puro (addgene Plasmid 1764). Retroviral vector was transfected into the Phoenix-AMPHO packaging cells (CRL-3213™; ATCC), and cultured in Dulbecco’s modified Eagle’s medium (DMEM) medium (Gibco, Invitrogen, Herlev, Denmark) supplemented with 10 % FCS until 70–80 % confluent, using the FuGENE 6 (Roche, Hvidovre, Denmark) reagent according to the manufacturer’s instruction. For cell transduction, the supernatants collected from Phoenix packaging cells containing virus particles were filtered with a 0.45 μm filter and added to hBMSC-TERT cell lines in the presence of Polybrene (Sigma, Brøndby, Denmark). Stably transduced cells were selected on the antibiotic puromycine.
In vitro differentiation studies and cytochemical staining
Cells were plated at 2 × 104 cells/cm2 in six-well plates in modified Dulbecco medium (MEM) (GIBCO, Herlev, Denmark) containing 10 % fetal bovine serum (FBS) (GIBCO), 100 U/ml penicillin (GIBCO), and 100 μg/ml streptomycin (GIBCO). One day after plating, media were replaced with osteogenic media including culture media supplemented with 10 mM β-glycerophosphate (Sigma, Brøndby, Denmark), 10 nM dexamethasone (Sigma, Brøndby, Denmark), 10 nM Calcitriol (1,25-dihydroxyvitamin D3; Leo, Ballerup, Denmark), and 50 μg/ml vitamin C (Sigma, Brøndby, Denmark). The osteogenic medium was changed every other day for 2 weeks.
For adipocyte differentiation, cells were plated at 3 × 104 cells/cm2 in six-well plates in culture media. The day after, culture media was removed and adipogenic media was added including MEM (GIBCO), 10 % FBS (GIBCO), 10 % horse serum (GIBCO), 100 U/ml penicillin (GIBCO), 100 μg/ml streptomycin (GIBCO), 100 nM dexamethasone (Sigma, Brøndby, Denmark), 0.25 mM 3-isobutyl-1-methylxanthine (IBMX, Gibco, Herlev, Denmark), 1 μM BRL (Rosiglitazone; Sigma, Brøndby, Denmark), and 500 nM insulin (Sigma, Brøndby, Denmark,). The adipogenic medium was changed every other day for 20 days.
Alkaline phosphatase staining
At 7 and 14 days during osteoblast differentiation, cells were fixed with acetone/citrate (10 mM) buffer, pH 4.2 (1.5:1), for 5 minutes at room temperature and stained with the 1:1 mixture of 0.2 mg/ml Napthol-AS-TR-phosphate (Sigma) and 1.3 mg/ml Fast Red TR solution (Sigma) in 0.1 M Tris buffer, pH 9.0, for 1 hour at room temperature.
Alizarin Red staining
Alizarin Red staining was performed on days 7 and 14 of osteoblast differentiation to assess matrix mineralization of osteoblastic cultures. Cells were fixed with 70 % ice-cold ethanol for 1 hour at −20 °C, and stained with 40 mM Alizarin Red S (AR-S; Sigma, Brøndby, Denmark), pH 4.2, for 10 minutes at room temperature. The amount of mineralized matrix (bound stain) was quantified by eluting the Alizarin Red stain, using 20-minute incubation of the cultures in 10 % (w/vol) cetylpyridinium chloride solution. The absorbance of the eluted stain was measured at 570 nm using FLUOstar Omega multimode microplate reader (Leica, Ballerup, Denmark).
Oil Red O staining
Adipocytic cultures from days 10 and 20 were fixed with 4 % paraformaldehyde for 10 minutes at room temperature, rinsed with 3 % isopropanol solution, and stained with Oil Red O (Sigma, USA) solution for 1 hour at room temperature. Oil Red O solution consist of 25 mg Oil Red O dye, 5 ml of 60 % isopropanol, and 3.35 ml H2O. The bound dye was eluted by 100 % isopropanol and absorbance measured at 490 nm using FLUOstar Omega multimode microplate reader.
Images were taken using a Leica® DM4500B microscope equipped with a motorized stage (Märzhäuser Wetzlar GmbH, Wetzlar, Germany) and a Leica® DFC300 FX camera. Images were acquired with the Surveyor® software (Objective Imaging, Leica, Ballerup, Denmark).
Boyden chamber transwell migration assay
For migration experiments, 90 % confluent cells were starved for 24 hours in low-glucose Dulbecco’s minimal essential medium (LG-DMEM; GIBCO, Herlev, Denmark) supplemented with 0.1 % (w/v) BSA (Sigma, Brøndby, Denmark), harvested, and assayed using the AP48 Migration Assay system and 8 μm pore size polycarbonate membranes (NeuroProbe, Inc., USA). Membranes were coated for 1 hour at 37 °C with 5 μg/ml fibronectin (human placenta; Sigma) and 10 μg/ml rat tail collagen type I (Sigma, Brøndby, Denmark) in migration media consisting of LG-DMEM (GIBCO) supplemented with 0.2 % FBS (GIBCO). hBMSC-TERT cells in LG-DMEM with 0.1 % BSA (control medium) with or without pretreatment (see later) were added to the wells. Control media with or without PDGFs were added to the lower chamber. The assembled assay was incubated at 37 °C for 4 hours. After removal of cells on the upper side of the membrane with a rubber scraper, cell migration through the membrane was quantitated by staining cells with the Hemacolor staining kit (Merck), and cells were counted from one picture/well taken at 10× magnification by a computer-assisted cell-counting program (Visiopharm, Hørsholm, Denmark). The following concentrations of factors were used: PDGF-BB (PeproTech GmbH Hamburg, Germany), 10–100 ng/ml; PDGF-AA (PeproTech GmbH Hamburg, Germany), 100 ng/ml; PDGF-AB (PeproTech GmbH Hamburg, Germany), 100 ng/ml; and PDGFRβ inhibitor (SU16f; Tocris), 10–200 nM. Then 1 μl/ml vehicle (dimethyl sulfoxide (DMSO)) were added to the cells just prior to assay assembly.
In vivo homing studies
Ethical permission for animal study was granted from the Danish National Authority on animal experiments (2012-15-2934-00559). Immunodeficient NOD scid (NOD.CB17-Prkdc scid /J) mice were used and operations were performed under intraperitoneal injection of ketamine (100 mg/kg) and xylazine (10 mg/kg). Temgesic (100 mg/kg) was given subcutaneously every 8 hours to relieve postoperative pain for at least 3 days post procedure.
Closed femur fracture
The closed femur fracture procedure was carried out as described by Bonnarens and Einhorn . In brief, an incision was made just medial to the patella, and then the longitudinal fibers of the quadriceps were divided and the patella dislocated laterally exposing the femur condyles. A 0.5 mm needle was inserted between the condyles and drilled in a retrograde fashion into the medullary cavity of the femur (functioning as a nail fixating the fracture). Afterwards, a standardized reproducible closed fracture was induced using a guillotine fracture apparatus . A radiograph was acquired following fracture induction (using Faxitron MX-20; Faxitron X-Ray , Lincolnshire, USA) to confirm the fracture induction and needle position.
Transplantation studies and bioluminescence imaging
Immunodeficient NOD/SCID (NOD/LtSz/Prkdc scid /J) 12-week-old female mice (n = 24) were randomly divided into the two main fracture study groups. The hBMSC-TERT+Bone (n = 12) or hBMSC-TERT–Bone cells (n = 12) were injected intravenously via tail vein at a dose of 1 × 106 cells/mouse suspended in 200 μl PBS per injection. For noninvasive imaging analysis, mice were injected with d-Luciferin (150 mg/kg, intraperitoneal injection; Caliber Life Sciences), and after 10 minutes the mice were placed in a supine position for image acquisition with 5 minutes of acquisition time using the IVIS Imaging System, 200 Series (Xenogen Corp., Rodgau, Germany). Bioluminescence color images were overlaid on gray-scale photographic images of the mice to allow for localization of the light source using the Livingimage™ software (V. 2.11) overlay (Xenogen Corp., Rodgau, Germany). To equalize comparisons across animals and between groups, the scales were fixed. Regions of interest (ROI) were manually selected, and signal intensity was expressed in terms of number of photons/cm2/second.
Data are presented as mean ± standard error of the mean SEM. Statistical analysis was performed using unpaired, two-tailed Student’s t test for comparison between two groups and one-way analysis of variance for comparison between more than two groups using Graphpad Prism version 5.0 software (Graphpad Prism Software, Inc., La Jolla, CA, USA). For all experiments, p <0.05 compared with control was considered significant.
hBMSC-TERT+Bone cells demonstrate high ex vivo directed migratory function
To study whether the ability to form heterotopic bone in vivo is associated with an enhancement of the migratory capacity of hBMSC, we compared the ex vivo transwell migration ability of two populations derived from hBMSC-TERT cells: one capable of in vivo bone formation (hBMSC-TERT+Bone) versus one that lacks this ability (hBMSC-TERT–Bone) . As shown in Fig. 1a, the transwell migration ability of hBMSC-TERT+Bone toward PDGF-BB was significantly higher compared with hBMSC-TERT–Bone.
hBMSC-TERT+Bone cells exhibit distinct cytoskeletal properties
Since cell migration includes a cascade of cellular events of cell spreading and focal adhesion formation , we compared cytoskeletal and adhesion properties of hBMSC-TERT+Bone versus hBMSC-TERT–Bone when cultured on fibronectin-coated plates for 2 hours and stained for F-actin and vinculin. As shown in Fig. 1b, hBMSC-TERT+Bone exhibited more cells spreading with formation of cytoskeletal projections, increased mean cell area, and reduced mean cell roundness. In addition, the focal adhesion number per cell area, as visualized by vinculin staining, was increased in hBMSC-TERT+Bone (Fig. 1b).
hBMSC-TERT+Bone cells exhibit enhanced in vivo homing to fractured bone site
We further compared the in vivo homing ability of hBMSC-TERT+Bone versus hBMSC-TERT–Bone to bone fracture sites in mice. hBMSC-TERT+Bone and hBMSC-TERT–Bone were stably transduced with a retroviral vector expressing the luciferase gene 2 (Luc2). Luc-overexpressing cell lines exhibited high levels of luciferase enzyme activity and maintained their differentiation capacity similar to the parental cell lines (Figure S1a, b in Additional file 1). hBMSC-TERT-Luc+Bone and hBMSC-TERT-Luc–Bone were administrated intravenously via tail vein injection in SCID/NOD mice with a stabilized femur bone fracture immediately following fracture induction. Noninvasive bioluminescence imaging using the IVIS Imaging System showed trapping of both cell lines in the lungs after 1 day (Fig. 1c). hBMSC-TERT+Bone cells and not hBMSC-TERT–Bone, were able to migrate to the fracture sites by day 6 (Fig. 1c).
Molecular phenotype of hBMSC single-cell clones with enhanced migratory function
Since hBMSC-TERT+Bone and hBMSC-TERT–Bone are nonclonal populations, we further studied the migratory capacity of three hBMSC-TERT-derived single-cell clones with either HBF capacity (HBF1, HBF2, and HBF3) or three cell clones with LBF capacity (LBF1, LBF2, and LBF3) . HBF clones showed higher migratory capacity toward PDGF-BB as compared with LBF clones, supporting our hypothesis of the existence of an association between the in vivo bone formation capacity and ex vivo migratory capacity of hBMSC (Fig. 2a).
To identify the molecular signature underlying enhanced hBMSC migratory function, we performed a global DNA microarray analysis comparing the transcriptome profile of hBMSC-TERT+Bone, hBMSC-TERT–Bone, the three HBF clones (HBF1, HBF2, and HBF3), and the three LBF clones (LBF1, LBF2, and LBF3). Cluster analysis revealed molecular clustering of HBF clones to hBMSC-TERT+Bone and LBF clones to hBMSC-TERT–Bone (Fig. 2b). We identified 427 genes (cutoff level ≥2-fold, p 0.05) that are upregulated in hBMSC-TERT+Bone and HBF clones compared with hBMSC-TERT–Bone and LBF clones (Table S2 in Additional file 1). Gene annotation analysis based on molecular functions revealed significant enrichment of gene categories of cell migration (12 %), cell adhesion (17 %), cytoskeletal genes (12 %), and bone-related genes (16 %) in hBMSC-TERT+Bone and HBF clones (Fig. 2c, Table 1). Some of the significantly upregulated genes by HBF clones were confirmed by real-time RT-PCR (Fig. 2d).
PDGFRα and PDGFRβ are associated with enhanced hBMSC migration
To identify surface marker(s) that can be used prospectively to isolate the hBMSC population with enhanced migratory functions, we examined differentially expressed genes in HBF clones. Both PDGFRα (CD140A) and PDGFRβ (CD140B) were found to be significantly upregulated by hBMSC-TERT+Bone and HBF clones versus hBMSC-TERT–Bone and LBF clones by 4.16-fold and 2.14-fold respectively (Table 1). These data have also been confirmed by quantitative RT-PCR analysis (Fig. 3a, b).
To test for the specific contribution of PDGFRs to migratory functions, we tested the effect of PDGF and other growth factors known to exert chemoattractant function: insulin-like growth factor 1 (IGF1), stromal cell-derived factor 1 (SDF1), and tumor necrosis factor alpha (TNFα). The ex vivo transwell migratory capacity of hBMSC-TERT+Bone was significantly higher toward PDGF-BB compared with other chemoattractants (Fig. 4a). We next confirmed the specificity of chemotactic effects of PDGF isoforms on hBMSC-TERT+Bone. As shown in Fig. 4b, hBMSC-TERT+Bone exhibited higher migration toward heterodimeric PDGF-AB or homodimeric PDGF-AA and PDGF-BB proteins and their migration was decreased by adding the same concentration of PDGF isoform to both lower and upper chambers in the transwell migration assay (Fig. 4b). Furthermore, we demonstrated a dose-dependent stimulatory effect of PDGF-BB on increasing the migration of hBMSC-TERT+Bone versus hBMSC-TERT–Bone (Fig. 4c). In addition, selective inhibition of PDGFRβ signaling in hBMSC-TERT+Bone using the chemical inhibitor SU-16f suppressed the migration of hBMSC-TERT+Bone toward PDGF-BB in a dose-dependent manner (Fig. 4d).
Identifying PDGFRβ as a marker for prospective isolation of a population of hBMSC with migratory function
To study the possible use of PDGFRβ as a candidate marker for isolating hBMSC with increased migratory capacity, MACS-enriched PDGFRβ+ cells from hBMSC-TERT–Bone (Fig. 5a) were obtained and the PDGFRβ+ cells were examined for their response to PDGF isoforms in the transwell migration assay. Sorted PDGFRβ+ cells exhibited enhanced chemotactic migration toward different PDGF isoforms as compared with the PDGFRβ-negative population (Fig. 5b).
In the present study, we demonstrate that a subpopulation of hBMSC with known ability to form bone in vivo exhibits enhanced ability to migrate in vitro and to home to bone fractures in vivo. In addition, we identified a molecular signature that defines this population and one surface marker, PDGFRβ (CD 140b), that can be employed for prospective isolation of the cells with enhanced migratory capacity from heterogeneous hBMSC cultures.
hBMSC are capable of multilineage differentiation into various mesoderm-type cells including osteoblastic cells , which is the basis for their use in cell-based therapies for bone tissue regeneration; for example, in repair of bone defects and nonhealed fractures [34–36]. However, cultured hBMSC represent a heterogeneous population with respect to differentiation capacity (so-called progenitor functions) as well as nonprogenitor functions. For example, not all hBMSC maintained under standard culture conditions are capable of bone formation . Similarly, when hBMSC are infused intravenously, only a small population of cells home to injured tissue . The efficient use of hBMSC in therapy thus requires defining an ex vivo molecular phenotype predictive for the in vivo behavior . In the present study, we have investigated the correlation between two characteristics: cellular migratory behavior and ability to form heterotopic bone, as well as homing capacity to bone fractures in vivo. For this purpose, we employed previously characterized hBMSC clones with various capacities for in vivo heterotopic bone formation. Our data suggest that in vivo heterotopic bone formation capacity is associated with enhanced in vivo migratory abilities to bone fractures.
By comparing the transcriptome profile of hBMSC-TERT+Bone and hBMSC-TERT–Bone and clones with different heterotopic bone-forming capacities, we identified a molecular signature linking homing and in vivo bone formation potential. This molecular signature included expression of high levels of known genes important for cell migration and homing. Among these were chemokine (C-X-C motif) receptor 7 (CXCR7) and chemokine (C-X-C motif) ligand 12 (CXCL12) , extracellular matrix, and adhesion molecules (e.g., collagen, type XVI, alpha 1 (COL16A1) and vascular cell adhesion molecule 1 (VCAM-1)) .
We have demonstrated that hBMSC with bone-forming capacity can home to bone fracture, suggesting that these migratory cells are recruited from the circulation under physiological conditions of bone formation and bone regeneration. However, the contribution of homed hBMSC to bone regeneration at bone fracture sites has been a contentious issue . Some studies have demonstrated that circulating hBMSC can home and contribute to fracture healing [36, 41, 42]. Also, patients with recent fractures have a large number of circulating osteoprogenitor (i.e., hBMSC-like) cells in their peripheral circulation, providing circumstantial evidence of recruiting bone-forming cells to bone fracture sites . On the other hand, some studies have demonstrated that the contribution of homed hBMSC to bone fracture healing is limited and that fracture healing is mediated by locally recruited cells possibly from periosteum or bone marrow [44, 45]. Our study suggests that there exists a population of hBMSC with enhanced bone-forming and migratory capacity. Further prospective studies are needed to quantitate the contribution of this cell population to fracture healing and bone regeneration.
Some of the molecules identified in the molecular signature defining high bone formation and migratory capacity do not seem to contribute to the migratory phenotype. For example, CXCR7 expression was highly upregulated in cells with HBF capacity, but we were unable to detect an increased ex vivo chemotactic migration toward its ligand SDF-1. It is possible that other noncanonical receptors are involved in regulating the migratory function of hBMSC. It is also plausible that hBMSC differ from leucocytes regarding the molecules employed for transendothelial migration. In support of this notion, we did not detect increased homing of hBMSC to bone fracture in vivo in hBMSC overexpressing CXCR4 (data not shown), a known factor enhancing leukocyte and HSC homing .
We identified high expression levels of PDGFRα and PDGFRβ to be associated with enhanced chemotactic migration of hBMSC toward PDGF isoforms. Therefore, PDGFRβ can be used as a possible marker to enrich for hBMSC with high ex vivo migratory capacity. PDGFs, via binding to their cognate receptors PDGFRα and PDGFRβ, have been reported previously as potent chemoattractants for multiple cell types including hBMSC [47–50]. Our study corroborates and extends previous findings by showing the possible prospective isolation of hBMSC with enhanced migratory capacity based on PDGFRβ expression.
Our study has some limitations. First, we employed a telomerized cell line (hBMSC-TERT) and not primary hBMSC for isolating subpopulations of BMSC. hBMSC-TERT has the advantage of having a stable phenotype during long-term culture and our extensive experience with hBMSC-TERT demonstrate that these cells maintain a phenotype similar to primary hBMSC in vitro and in vivo  and has comparable molecular signature to primary hBMSC . While hMSC-TERT has been employed in a clinical trial as a vehicle for drug delivery , the cells are not suitable for clinical trials for bone regeneration due to their genetic modification and the risk of transformation. Future studies on the relevance of our findings to in vivo bone regeneration should thus be conducted using primary human cells in preclinical animal models . Second, the identified PDGFRβ can be both a marker and functional protein contributing to enhance migration. Our studies do not distinguish between these two possibilities. Third, prospective isolation of PDGFRβ+ cells from heterogeneous hBMSC cultures that can both home and contribute to fracture healing needs to be performed.
This study demonstrates that cultured hBMSC, in addition to cellular heterogeneity, exhibit functional heterogeneity and that bone formation and homing abilities are closely linked characteristics that can be defined by a common molecular signature. Further studies are needed to determine the role of PDGFRβ in hBMSC homing and to isolate hBMSC populations based on their molecular signature to be tested in preclinical animal models for enhancing fracture healing and bone regeneration.
Bone marrow-derived stromal cells
Bovine serum albumin
Ethylenediamine tetraacetic acid
Fetal bovine serum
Fetal calf serum
False discovery rate
High bone forming
Human bone marrow-derived stromal stem cells
Hematopoietic stem cell
Insulin-like growth factor 1
Low bone forming
Low-glucose Dulbecco’s minimal essential medium
Minimal essential medium
Platelet-derived growth factor
Platelet-derived growth factor receptor
Regions of interest
Stromal cell derived factor 1
Standard error of the mean
Tumor necrosis factor alpha
Jafari A, Harkness L, Zaher W, Kassem M. Adult stromal (skeletal, mesenchymal) stem cells: advances towards clinical applications. In: Turksen K, editor. Adult stem cells. New York: Springer; 2014. p. 359–73.
Nejadnik H, Hui JH, Feng Choong EP, Tai BC, Lee EH. Autologous bone marrow-derived mesenchymal stem cells versus autologous chondrocyte implantation: an observational cohort study. Am J Sports Med. 2010;38:1110–6.
Quarto R, Mastrogiacomo M, Cancedda R, Kutepov SM, Mukhachev V, Lavroukov A, et al. Repair of large bone defects with the use of autologous bone marrow stromal cells. N Engl J Med. 2001;344:385–6.
Hernigou P, Poignard A, Beaujean F, Rouard H. Percutaneous autologous bone-marrow grafting for nonunions. Influence of the number and concentration of progenitor cells. J Bone Joint Surg Am. 2005;87:1430–7.
Yelda T, Berrin U, Murat S, Aytac O, Sevgi B, Yasemin S, et al. Intracoronary stem cell infusion in heart transplant candidates. Tohoku J Exp Med. 2007;213:113–20.
Suzuki K, Murtuza B, Suzuki N, Smolenski RT, Yacoub MH. Intracoronary infusion of skeletal myoblasts improves cardiac function in doxorubicin-induced heart failure. Circulation. 2001;104:I213–7.
Yang Q, Zhang J, Jiang J. Intracoronary transplantation of genetically modified mesenchymal stem cells, a novel method to close muscular ventricular septal defects. Med Hypotheses. 2011;77:505–7.
Chen SL, Fang WW, Ye F, Liu YH, Qian J, Shan SJ, et al. Effect on left ventricular function of intracoronary transplantation of autologous bone marrow mesenchymal stem cell in patients with acute myocardial infarction. Am J Cardiol. 2004;94:92–5.
Magnon C, Frenette PS. Hematopoietic stem cell trafficking. In: StemBook. Cambridge, MA: Harvard Stem Cell Institute; 2008.
Ji JF, He BP, Dheen ST, Tay SS. Interactions of chemokines and chemokine receptors mediate the migration of mesenchymal stem cells to the impaired site in the brain after hypoglossal nerve injury. Stem Cells. 2004;22:415–27.
Shi M, Li J, Liao L, Chen B, Li B, Chen L, et al. Regulation of CXCR4 expression in human mesenchymal stem cells by cytokine treatment: role in homing efficiency in NOD/SCID mice. Haematologica. 2007;92:897–904.
Devine SM, Bartholomew AM, Mahmud N, Nelson M, Patil S, Hardy W, et al. Mesenchymal stem cells are capable of homing to the bone marrow of non-human primates following systemic infusion. Exp Hematol. 2001;29:244–55.
Horwitz EM, Gordon PL, Koo WKK, Marx JC, Neel MD, McNall RY, et al. Isolated allogeneic bone marrow-derived mesenchymal cells engraft and stimulate growth in children with osteogenesis imperfecta: implications for cell therapy of bone. Proc Natl Acad Sci U S A. 2002;99:8932–7.
Yang X, Balakrishnan I, Torok-Storb B, Pillai MM. Marrow stromal cell infusion rescues hematopoiesis in lethally irradiated mice despite rapid clearance after infusion. Adv Hematol. 2012;2012:142530.
Lange C, Brunswig-Spickenheier B, Cappallo-Obermann H, Eggert K, Gehling UM, Rudolph C, et al. Radiation rescue: mesenchymal stromal cells protect from lethal irradiation. PLoS One. 2011;6, e14486.
Eggenhofer E, Benseler V, Kroemer A, Popp FC, Geissler EK, Schlitt HJ, et al. Mesenchymal stem cells are short-lived and do not migrate beyond the lungs after intravenous infusion. Front Immunol. 2012;3:297.
Schrepfer S, Deuse T, Reichenspurner H, Fischbein MP, Robbins RC, Pelletier MP. Stem cell transplantation: the lung barrier. Transplant Proc. 2007;39:573–6.
Karp JM, Leng Teo GS. Mesenchymal stem cell homing: the devil is in the details. Cell Stem Cell. 2009;4:206–16.
Kuznetsov SA, Krebsbach PH, Satomura K, Kerr J, Riminucci M, Benayahu D, et al. Single-colony derived strains of human marrow stromal fibroblasts form bone after transplantation in vivo. J Bone Miner Res. 1997;12:1335–47.
Muraglia A, Cancedda R, Quarto R. Clonal mesenchymal progenitors from human bone marrow differentiate in vitro according to a hierarchical model. J Cell Sci. 2000;113:1161–6.
Larsen KH, Frederiksen CM, Burns JS, Abdallah BM, Kassem M. Identifying a molecular phenotype for bone marrow stromal cells with in vivo bone forming capacity. J Bone Miner Res. 2009;25:796.
Post S, Abdallah BM, Bentzon JF, Kassem M. Demonstration of the presence of independent pre-osteoblastic and pre-adipocytic cell populations in bone marrow-derived mesenchymal stem cells. Bone. 2008;43:32–9.
Simonsen JL, Rosada C, Serakinci N, Justesen J, Stenderup K, Rattan SI, et al. Telomerase expression extends the proliferative life-span and maintains the osteogenic potential of human bone marrow stromal cells. Nat Biotechnol. 2002;20:592–6.
Abdallah BM, Haack-Sorensen M, Burns JS, Elsnab B, Jakob F, Hokland P, et al. Maintenance of differentiation potential of human bone marrow mesenchymal stem cells immortalized by human telomerase reverse transcriptase gene despite [corrected] extensive proliferation. Biochem Biophys Res Commun. 2005;326:527–38.
Ambion. Illumina TotalPrep RNA Amplification Kit. www.ambion.com.
Illumina HumanRef-8 v3 Expression BeadChips. www.illumina.com.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B (Stat Method) 1995:289–300.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4:3.
Bonnarens F, Einhorn TA. Production of a standard closed fracture in laboratory animal bone. J Orthop Res. 1984;2:97–101.
Parsons JT, Horwitz AR, Schwartz MA. Cell adhesion: integrating cytoskeletal dynamics and cellular tension. Nat Rev Mol Cell Biol. 2010;11:633–43.
Abdallah BM, Jafari A, Zaher W, Qiu W, Kassem M. Skeletal (stromal) stem cells: an update on intracellular signaling pathways controlling osteoblast differentiation. Bone 2014
Tseng SS, Lee MA, Reddi AH. Nonunions and the potential of stem cells in fracture-healing. J Bone Joint Surg Am. 2008;90:92–8.
Granero-Molto F, Weis JA, Miga MI, Landis B, Myers TJ, O'Rear L, et al. Regenerative effects of transplanted mesenchymal stem cells in fracture healing. Stem Cells. 2009;27:1887–98.
Devine MJ, Mierisch CM, Jang E, Anderson PC, Balian G. Transplanted bone marrow cells localize to fracture callus in a mouse model. J Orthop Res. 2002;20:1232–9.
Wynn RF, Hart CA, Corradi-Perini C, O'Neill L, Evans CA, Wraith JE, et al. A small proportion of mesenchymal stem cells strongly expresses functionally active CXCR4 receptor capable of promoting migration to bone marrow. Blood. 2004;104:2643–5.
Dar A, Kollet O, Lapidot T. Mutual, reciprocal SDF-1/CXCR4 interactions between hematopoietic and bone marrow stromal cells regulate human stem cell migration and development in NOD/SCID chimeric mice. Exp Hematol. 2006;34:967–75.
Aldridge V, Garg A, Davies N, Bartlett DC, Youster J, Beard H, et al. Human mesenchymal stem cells are recruited to injured liver in a beta1 integrin and CD44 dependent manner. Hepatology. 2012;56:1063.
Bielby R, Jones E, McGonagle D. The role of mesenchymal stem cells in maintenance and repair of bone. Injury. 2007;38:S26–32.
Shirley D, Marsh D, Jordan G, McQuaid S, Li G. Systemic recruitment of osteoblastic cells in fracture healing. J Orthop Res. 2005;23:1013–21.
Tu Q, Valverde P, Li S, Zhang J, Yang P, Chen J. Osterix overexpression in mesenchymal stem cells stimulates healing of critical-sized defects in murine calvarial bone. Tissue Eng. 2007;13:2431–40.
Eghbali-Fatourechi GZ, Lamsam J, Fraser D, Nagel D, Riggs BL, Khosla S. Circulating osteoblast-lineage cells in humans. N Engl J Med. 2005;352:1959–66.
Arnsdorf EJ, Jones LM, Carter DR, Jacobs CR. The periosteum as a cellular source for functional tissue engineering. Tissue Eng A. 2009;15:2637–42.
Einhorn TA. The cell and molecular biology of fracture healing. Clin Orthop Relat Res 1998;355 Suppl:S7–21.
Lapidot T, Kollet O. The essential roles of the chemokine SDF-1 and its receptor CXCR4 in human stem cell homing and repopulation of transplanted immune-deficient NOD/SCID and NOD/SCID/B2m(null) mice. Leukemia. 2002;16:1992–2003.
Facchiano A, De Marchis F, Turchetti E, Facchiano F, Guglielmi M, Denaro A, et al. The chemotactic and mitogenic effects of platelet-derived growth factor-BB on rat aorta smooth muscle cells are inhibited by basic fibroblast growth factor. J Cell Sci. 2000;113:2855–63.
Seppa H, Grotendorst G, Seppa S, Schiffmann E, Martin GR. Platelet-derived growth factor in chemotactic for fibroblasts. J Cell Biol. 1982;92:584–8.
Ozaki Y, Nishimura M, Sekiya K, Suehiro F, Kanawa M, Nikawa H, et al. Comprehensive analysis of chemotactic factors for bone marrow mesenchymal stem cells. Stem Cells Dev. 2007;16:119–29.
Phipps MC, Xu Y, Bellis SL. Delivery of platelet-derived growth factor as a chemotactic factor for mesenchymal stem cells by bone-mimetic electrospun scaffolds. PLoS One. 2012;7:e40831.
Al-Nbaheen M, Vishnubalaji R, Ali D, Bouslimi A, Al-Jassir F, Megges M, et al. Human stromal (mesenchymal) stem cells from bone marrow, adipose tissue and skin exhibit differences in molecular phenotype and differentiation potential. Stem Cell Rev. 2012.
Serakinci N, Guldberg P, Burns JS, Abdallah B, Schrodder H, Jensen T, et al. Adult human mesenchymal stem cell as a target for neoplastic transformation. Oncogene. 2004;23:5095–8.
This study was supported by a grant from the Danish Innovation fond (“Efficient Homing into Injured Tissues: Systemic Cell-based Therapy using Mesenchymal Stem Cells (MSC)”, 09-067112), Odense University Hospital, Odense, Denmark. WZ has received a fellowship from the Saudi Ministry of Higher Education and support from the College of Medicine Research Center, Deanship of Scientific Research, King Saud University.
The authors declare that they have no competing interests.
RKA and WZ carried out in vitro migration studies and in vivo homing experiments, analyzed data, and revised the manuscript. KHL and ND performed molecular genetics and in vivo studies, analyzed data, and revised the manuscript. KD, WW, and JA carried out and analyzed the microarray study and revised the manuscript. BMA performed the cytoskeletal study, contributed to the study design, and drafted and revised the manuscript. MK conceived of the study, and participated in the study design and revised the manuscript. All authors read and approved the manuscript.
Rikke K. Andersen and Walid Zaher contributed equally to this work.
Figure S1. Showing generation and differentiation of hBMSC-TERT-Luc cell lines, Table S1 presenting a list of primers used for RT-PCR, and Table S2 presenting genes upregulated by group 1 (hMSC-TERT+Bone and three HBF clones) over group 2 (hMSC-TERT–Bone and three LBF clones).(PDF 317 kb)
About this article
Cite this article
Andersen, R.K., Zaher, W., Larsen, K.H. et al. Association between in vivo bone formation and ex vivo migratory capacity of human bone marrow stromal cells. Stem Cell Res Ther 6, 196 (2015). https://doi.org/10.1186/s13287-015-0188-9