MicroRNAs Regulate Osteogenesis and Chondrogenesis of Mouse Bone Marrow Stromal Cells

Abstract

MicroRNAs (miRNAs) are non-coding RNAs that bind to target mRNA leading to translational arrest or mRNA degradation. To study miRNA-mediated regulation of osteogenesis and chondrogenesis, we compared the expression of 35 miRNAs in osteoblasts and chondroblasts derived from mouse marrow stromal cells (MSCs). Differentiation of MSCs resulted in up- or downregulation of several miRNAs, with miR-199a expression being over 10-fold higher in chondroblasts than in undifferentiated MSCs. In addition, miR-124a was strongly upregulated during chondrogenesis while the expression of miR-96 was substantially suppressed. A systems biological analysis of the potential miRNA target genes and their interaction networks was combined with promoter analysis. These studies link the differentially expressed miRNAs to collagen synthesis and hypoxia, key pathways related to bone and cartilage physiology. The global regulatory networks described here suggest for the first time how miRNAs and transcription factors are capable of fine-tuning the osteogenic and chondrogenic differentiation of mouse MSCs.

Keywords

microRNA marrow stromal cells chondrogenesis osteogenesis

Introduction

MicroRNAs (miRNAs) are small non-coding RNA-molecules that bind to the 3’ untranslated region of mRNAs and, depending on their degree of complementarity with the target genes, induce translational repression or mRNA degradation.¹ Since the first identification of miRNAs in 1993,² hundreds of miRNAs have been identified from plants, animals and viruses.³ Consecutively, miRNAs have proven to play essential roles in diverse biological processes including early development,⁴, ⁵ cell proliferation and cell death,⁶ fat metabolism,⁷ cell differentiation,⁸, ⁹ and brain development.¹⁰ The sequences coding for miRNAs are spread around the genome, including exons, introns, 3‘-UTRs and genomic repeat-areas, and are situated either in the sense or antisense orientation with respect to the overlapping protein-coding gene. Studies carried out with embryonic stem (ES) cells indicate that miRNAs expression profiles in stem cells are different from other tissues,¹¹ suggesting that miRNAs may play an important role in stem cell self-renewal and differentiation. The expression of Argonaute genes is restricted to specific anatomical sites in mouse embryos, suggesting that short regulatory RNAs may have physiological functions during organogenesis.⁵ In addition, miRNAs are expressed in haematopoietic tissue where they participate in the regulation of haematopoietic stem cell (HSC) differentiation.⁹, ¹² However, there is only limited amount of information about the expression of miRNAs in mesenchymal stromal cells (MSCs) or their potential role in osteo- or chondrogenesis.¹³

Mesenchymal stromal cells are multipotent cells that have the potential to differentiate to various lineages of mesenchymal tissues, including bone, cartilage, adipose, tendon, and muscle.¹⁴ Compared to haematopoietic stem cells (HSCs), MSCs are rare in bone marrow, representing ~1 in 10,000 nucleated cells. Where HSCs are a well-characterized population of self-renewing cells that give rise to all mature blood cell lineages,¹⁵, ¹⁶ MSCs are less defined due to the limited understanding of MSC properties. As a result, terms such as “marrow stromal cells”, “mesenchymal progenitor cells”, “nonhaematopoietic mesenchymal stem cells” and “adult nonhaematopoietic stem cells” are used to define this cell population. Isolation and characterisation of stem cells from bone marrow rely on their immunophenotypic or functional aspects. Haematopoietic stem cells have been shown to express surface markers CD14, CD34 and CD45.¹⁷ MSCs, on the contrary, lack clearly defined surface markers and thus the isolation and characterisation of MSCs is still based on the properties described by Friedenstein already in 1970s: their adherence to plastic, spindle-shaped morphology and ability to form colonies.¹⁸, ¹⁹ In addition, MSCs are commonly characterised for their differentiation capacity and for the negative expression of haematopoietic surface markers.

In calcified tissue, MSCs are needed for bone and cartilage formation. During embryogenesis, bone formation begins with mesenchymal stem cell condensation. Membranous bone (craniofacial bones and the clavicle) is derived from MSCs that differentiate in situ into bone-forming osteoblasts and produce matrix rich in Type I collagen. Endochondral bone, which is the principal type of bone in the body, is formed by MSCs that first differentiate into chondrocytes to form a cartilagenous template for the bone. Chondrocytes secrete a matrix rich in Type II collagen and Aggrecan, and go through a genetic program driven by Sox9²⁰ leading to cartilage enlargement. In the centre of the cartilage anlage, chondrocytes become hypertrophic and start to synthesise Type X collagen that is later degraded and replaced by bone. Although transcription factors such as Sox9 and Runx2, and signalling molecules such as Indian hedgehog (Ihh), Parathyroid hormone-related protein (PTHrP), Fibroblast growth factors (FGF), and Bone morphogenetic proteins (BMPs) are involved in the regulation of endochondral bone formation,²¹ the molecular mechanisms leading to bone formation are still poorly understood. Thus, understanding the regulatory networks that control the lineage commitment and differentiation of MSCs is an important challenge.

In order to study the role of miRNAs in osteo- and chondrogenesis, miRNA expression profiles of osteoblasts and chondroblasts derived from mouse MSCs were compared. Subsequently, target prediction studies carried out with the differentially expressed miRNAs were combined with pathway analyses to gain more insight into the cellular functions potentially regulated by these miRNAs. Bioinformatics studies have shown that the promoter regions of miRNAs seem to contain similar regulatory motifs as the promoter regions of protein coding genes.²² In order to investigate whether the studied miRNAs could form regulatory networks with transcription factors (TFs) involved in osteo- or chondrogenesis, the promoter regions of the differentially expressed miRNAs were analysed. We present here multiple lines of evidence to suggest that in addition to haematopoietic cells, miRNAs are also involved in the regulation of lineage commitment in mesenchymal cells.

Materials and Methods

Cell culture and RNA extraction

All cell culture reagents, unless otherwise stated, were purchased from Gibco Invitrogen (U.S.A.). Total RNA was extracted from cultured cells before and after osteo- or chondrogenic induction using the mirVana miRNA Isolation Kit following the manufacturer's protocol (Ambion, U.S.A.). To remove genomic DNA contamination, total RNA samples were digested with DNase I (NEB, U.S.A.). RNA concentrations were quantified using an Eppendorf Biophotometer (Eppendorf, U.S.A.).

Bone marrow cells were isolated from 8–12 week-old male C57BL × DBA mice according to a previously described method.²³ Briefly, cells were isolated from the tibiae and femora by flushing them from the bone marrow cavity using a 10 ml syringe with a 25 gauge needle and medium consisting of RPMI-1640, 12% iFCS, 100 U/ml penicillin and 100 µg/ml streptomycin. A primary culture of plastic adherent cells from mouse bone marrow is a heterogeneous population of mesenchymal and hematopoietic stem cells.²⁴ For the selection of mesenchymal stem cells, bone marrow cells were incubated 2 hours at 37 °C on a plastic culture dish containing RPMI-1640 medium described above (12% iFCS, 100 U/ml penicillin and 100 µg/ml streptomycin) to remove rapidly adherent cells.¹⁸, ¹⁹ Unattached cells were collected and cultured in cell culture flasks at the initial density of 1 × 10⁶ cells/cm². Non-adherent cells were removed 48 hours later and adherent cells were washed with phosphate-buffered saline (PBS). Cells were further cultured with a twice-weekly medium replacement (half of the medium replaced). When confluent, cells were detached using trypsin-EDTA and re-plated at the density of 10 000 cells/cm². RPMI medium has been demonstrated to inhibit the growth of hematopoietic cells in culture²⁵ and cultures were therefore maintained in RPMI-1640 for 1 to 2 weeks.²⁶ Finally, adherent cells were detached by a trypsin-EDTA treatment and expanded by plating them in DMEM medium supplemented with 12% iFCS, 100 U/ml penicillin and 100 µg/ml streptomycin at the density of 1 000 cells/cm². Cells were cultured in described medium until confluent (1 to 2 weeks), thereafter trypsinized, immunophenotypically characterised and subjected to osteoblastic or chondrogenic differentiation.

For immunophenotypic characterisation, MSCs were plated on chamber slides, cultured to confluency and then stained for surface markers Ly-6A/E stem cell antigen 1 (Sca-1) (BD Biosciences, U.S.A.), CD34 (Vision Biosystems, U.S.A.) and CD45 (Dako, Denmark). Cells were incubated for 1 hour with primary antibodies against Sca-1 (1:100), CD45 (1:100) and CD34 (1:400) diluted in DAKO ChemMate Antibody Diluent (Dako, Denmark). Cells were washed with TNT-buffer (0.1 M Tris pH 7.5, 0.15 M NaCl, 0.05% Tween20) and incubated with Post-Blocking reagent (DPVO+500Post, Immunovision Technologies Co., Netherlands) for 20 minutes. Cells were washed three times with TNT buffer and incubated with Poly-HRP anti-Mouse IgG (Immunovision Technologies Co., Netherlands) for 30 minutes. The secondary antibody was detected with DAB (Zymed Invitrogen, U.S.A.) and nuclei were counterstained with Mayer's hematoxylin.

Osteogenic differentiation was induced by culturing the long-term selected MSCs in osteogenic medium consisting of phenol red-free αMEM, 12% iFCS, 10 mM Na-β-glycerophosphate (Fluka BioChemika, Switzerland), 50 µgrams/ml ascorbic acid 2-phosphate (Sigma-Aldrich, U.S.A.), 100 U/ml penicillin and 100 µg/ml streptomycin for 3 weeks in cell culture flasks and 24-well plates at the initial density of 10 000 cells/cm². During the first week, the culture medium was supplemented with 10 nM dexamethasone. The cultures were terminated by RNA extraction or fixation in 3% paraformaldehyde. To demonstrate osteoblastic differentiation, cells were stained for alkaline phosphatase (ALP) (Sigma-Aldrich, U.S.A.) and bone nodules were detected by von Kossa staining.²⁷ To induce chondrogenic differentiation, 200 000 of the cultured MSCs were placed in a 15-ml polypropylene tube and centrifuged (6 min 500 × g) to form a micromass pellet culture.²⁸ Cell pellets were cultured for 21 days in chondroinductive medium consisting of high-glucose DMEM supplemented with 10 ng/ml TGF-(33 (R&D Systems, UK), 10^–7 M dexamethasone, 50 µg/ml ascorbic acid 2-phosphate, 40 µg/ml L-proline (Sigma Aldrich, U.S.A.), 100|lg/ml sodium pyruvate (Sigma Aldrich, U.S.A.), 50 mg/ml ITS+ Premix (BD Biosciences, U.S.A.). Media were changed every 3 to 4 days. After three weeks of culture, the pellets were either lysed for RNA extraction or fixed for 2.5 hours in 4% paraformaldehyde. To evaluate chondrogenic differentiation, pellets were embedded in paraffin, cut into 5 µm sections and stained with toluidine blue for proteoglycans. Presence of type II collagen was detected by 6B3 monoclonal antibody raised against chicken type II collagen²⁹ following the method described earlier.³⁰

Osteogenic and chondrogenic gene expression

To further confirm the osteogenic and chondrogenic differentiation of the long-term selected MSCs, the total RNAs were analyzed with RT-PCR for the expression of selected osteogenic and chondrogenic transcripts (Table S1). The osteogenic genes included Type I collagen, Osteocalcin, Osterix (Sp7), and Runx2, whereas chondrogenic differentiation was evaluated based on the expression of Type II and X collagens and Sox9. One µg of DNase I treated total RNA was reverse transcribed using M-MLV reverse transcriptase (Promega Corporation, UK). cDNAs were amplified using DyNAzyme II DNA Polymerase (Finnzymes, Finland) for 35 cycles and analysed on 1,5% agarose gel. Amplification of GAPDH and L19 served as loading controls.

miRNA expression

A total of 35 miRNAs were selected for the follow-up based on their expression in haematopoietic tissues⁹, ³¹, ³² or mouse embryonic stem (ES) cells¹¹ or based on computational predictions on physiologically important genes related to bone and cartilage function (Table S2). The expression profiles of these miRNAs were detected by quantitative real-time PCR (qRT-PCR). Amplifications were performed using Taq DNA Polymerase (ABgene) and mirVana qRT-PCR miRNA Detection Kit (Ambion, U.S.A.) following the manufacturer's instructions. For each sample, a total of 37 different reactions were performed in triplicate with mirVana qRT-PCR Primer Sets. Out of these, 35 were specific for miRNAs and two, U6 snRNA and 5S rRNA, were used for normalisation. For qRT-PCR reactions, 50 ng of DNAse-treated total RNA was used and a no-template reaction was performed for each primer set. RT-reactions, 30 min at 37 °C followed by 10 min at 95 °C, were incubated in MJ Research PTC-200. Real-time PCR reactions were performed in MJ Research PTC-200 DNA Engine Cycler in optical strips. The reactions were incubated at 95 °C for 3 min, followed by 40 cycles of 95 °C for 15 sec and 60 °C for 30 sec. At the end, a dissociation analysis (melt-curve) from 56 °C to 90 °C was performed. Ct data were determined using default threshold settings. End-point reactions were analysed on a 15% neutral PAGE to discriminate between the correct amplification products (~90 bp) and potential primer dimers.

Micro-RNA expression data was normalised to U6 snRNA and 5S rRNA according to the manufacturer's recommendations. Relative quantification of miRNA expression was calculated with the 2^–ΔΔCt method³³ where undifferentiated cells were set as a calibrator sample. Standard error of the normalized expression was calculated by applying the differential equation of Gauss.³⁴ The miRNA expression values were compared between undifferentiated and differentiated cells, and between osteoblasts and chondroblasts. MicroRNAs, whose expression was changed at least 5-fold or 2-fold, respectively, were selected for target predictions and promoter analyses.

Target predictions

Target prediction tools TargetScan (http://www.targetscan.org), PicTar (http://pictar.bio.nyu.edu) and miRanda (http://www.microrna.org) were utilized in order to find out possible targets genes for the differentially expressed miRNAs. Since prediction algorithms often result in false positives,³⁵ the target gene lists of the three programs were combined and an intersection set, including only the target genes found with all three prediction algorithms, was created for each miRNA separately. Predicted target gene identifiers were first converted into a common nomenclature, and the results were then combined for each miRNA separately using data available from the latest Ensembl release 44.

Pathway analysis

In order to elucidate the physiological role of differentially expressed miRNAs, their target genes were analyzed through the use of Ingenuity Pathways Analysis version 5.0 (Ingenuity^® Systems, www.ingenuity.com). The intersection sets were uploaded into the application, and each gene identifier was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base. The genes in the intersection sets were overlaid onto a global molecular network developed from information contained in the Ingenuity Pathways Knowledge Base, and networks were then algorithmically generated based on their connectivity. The functional analysis identified the biological functions or diseases that were most significant to the intersection sets. The significance of the association between the intersection set and pathway was measured in two ways. First, a ratio of the number of genes from the intersection set that mapped to the pathway divided by the total number of genes that mapped to the pathway was calculated. Second, Fischer's exact test was used to calculate a p-value determining the probability that the association between the genes in the intersection set and the pathway is explained by chance alone. Graphical representations of the molecular relationships between genes/gene products were also produced. Genes/gene products were represented as nodes, and the biological relationship between two nodes was represented as an edge (line). All edges are supported by at least 1 reference from the literature. Nodes were displayed using various shapes that represent the functional class of the gene product.

Promoter analysis

To examine the potential transcription factors (TFs) involved in miRNA regulation, promoter analysis was performed and conserved binding sites for known TFs were predicted for the upstream regions of differentially expressed miRNA genes. As target prediction algorithms, transcription factor binding site (TFBS) prediction tools often result in false positives. In order to improve the authenticity of the predictions, a technique known as phylogenetic footprinting was applied. Thus, human and mouse orthological sequences were compared and only conserved binding sites were accepted for the analysis. Regulatory regions, 500 bp upstream and 100 bp downstream from the starting locus, of human and mouse orthological pre-miRNAs (Table S3) were retrieved from Ensemble release 44. TFBS predictions were performed by Conreal³⁶ using high quality vertebrate matrices from TRANSFAC professional 9.4.³⁷, ³⁸

Results

Cell culture

In order to compare miRNA expression in osteoblasts and chondroblasts, mesenchymal cells were isolated and enriched from bone marrow. To confirm the negative selection of HSCs and the positive selection of MSCs, cells were characterised for the expression of CD34, CD45 and Sca-1 surface markers (Fig. 1). After enrichment procedure, cells were positive for stem cell marker Sca-1³⁹ and negative for hematopoietic surface antigens CD34 and CD45.¹⁷ To compare miRNA expression in osteoblasts and chondroblasts, long-term selected MSCs were induced to differentiate into the desired cell types. For osteoblastic differentiation, cells were cultured for one to 3 weeks in osteogenic medium, followed by alkaline phosphatase- and von Kossa stainings (Fig. 2A-D). For chondrogenic differentiation, cells were cultured in chondrogenic medium as micromass pellet cultures for 3 weeks, followed by histological evaluation. Toluidine blue staining demonstrated the presence of proteoglycans and immunohistochemical staining showed deposition of type II collagen in the cell pellets (Fig. 2E–F).

Figure 1

Surface antigen expression in MSCs.

Figure 2

Evaluation of osteogenic and chondrogenic differentiation of MSCs

Osteogenic and chondrogenic gene expression

To further evaluate the phenotypes of the in vitro differentiated cell cultures, expression of osteogenic and chondrogenic genes was analysed with RT-PCR (Table S1). After a three-week culture in osteoinductive medium, long-term selected MSCs expressed osteoblast-related genes Type I collagen, Osteocalcin, Osterix and Runx2, but also Type X collagen and Sox9. Cells grown in chondroinductive conditions expressed Runx2, Type II collagen, Type X collagen and Sox9, and were negative for the expression of Type I collagen, Osteocalcin and Osterix (Fig. 2G).

miRNA expression

The mirVana qRT-PCR miRNA Detection Kit (Ambion) was used in order to compare miRNA expression in osteoblasts and chondroblasts. Undifferentiated MSCs were used as a control and data were normalized as described in Experimental procedures. Figure 3 shows the miRNA expression profiles of osteoblasts and chondroblasts derived from MSCs. Relative expression value of 1 represents the expression level of a specific miRNA in undifferentiated MSCs. Only 2 out of 35 miRNAs were differentially (>5-fold) expressed in osteoblasts compared to the undifferentiated cells whereas in chondroblasts the expression levels of 7 out of 35 miRNAs changed at least 5-fold during the differentiation (Table 1). In addition, 8 out of 35 miRNAs were differentially expressed (at least 2-fold) between osteoblasts and chondroblasts derived from long-term selected MSCs (Table 1).

Figure 3

miRNA expression in osteoblasts and chondroblasts derived from MSCs.

Table 1

Differentially expressed miRNAs in MSC-derived osteoblasts (O) and chondroblasts (C).

Comparison of osteoblasts and chondroblasts to MSCs (at least 5-fold difference)
Osteoblasts	Expression relative to MSCs	Chondroblasts	Expression relative to MSCs
miR-24	7,66	miR-18	0,44 ↓
miR-199b	5,30	miR-24	5,54
		miR-96	0,13 ↓
		miR-101	5,66
		miR-124a	8,18
		miR-199a	12,13
		miR-199b	7,46

Comparison between osteoblasts and chondroblasts (at least 2-fold difference)
O vs. C	Ratio	C vs. O	Ratio
miR-18	2,65	miR-124a	3,50
miR-31	3,05	miR-130a	2,61
miR-96	5,73	miR-142-3p	2,35
miR-130b	3,31	miR-199a	5,25

Pathway analysis

Intersection sets (Fig. 4) of potential miRNA target genes were uploaded into Ingenuity Pathways Analysis version 5.0 (Ingenuity^® Systems, www.ingenuity.com), resulting in multiple interaction networks. An exceptionally high score and low p-value, as calculated by IPA, was observed for 5 miRNAs out of 11 and the most significant biological functions related to individual miRNA target gene pathways are shown in Table 2.

Figure 4

Intersection sets of the differentially expressed miRNAs constitute complex interaction networks.

Table 2

Systemic analysis of specific biological functions related to predicted miRNA target genes.

miRNA	Specific biological functions	–log(p-value)	Focus genes
miR-96	Developmental disorder	5,0	13
	Cellular movement	4,1	10
	Cellular assembly and organisation	4,0	18
	Cell morphology	4,0	17
	Connective tissue development and function	4,0	10
	Cancer	3,2	14
	Small molecule biochemistry	3,0	19
	Lipid metabolism	3,0	11
	Genetic disorder	2,5	10
miR-124a	Cancer	3,5	11
	Cellular development	3,0	16
	Cellular movement	2,5	14
miR-130a	Cellular movement	5,2	13
	Cell signalling	3,3	18
	Gene expression	3,0	14
miR-130b	Cell cycle	4,5	12
	Cell death	3,5	15
	Cellular movement	2,8	12
miR-199a	Cell death	4,9	10
	Cancer	3,0	10

A total of 79 biological functions were evaluated for each miRNA intersection set by Ingenuity Pathways Analysis. The statistically significant biological functions with 10 or more focus genes (target genes of a given miRNA) are listed here.

Promoter analysis

Information about eukaryotic transcription factors, their genomic binding sites and DNA-binding profiles is stored in databases such as TRANSFAC³⁷ and JASPAR.⁴⁰ The upstream regions of differently expressed miRNA genes (Table S2) were analysed in order to study the potential transcription factors (TFs) involved in miRNA regulation. When the results were combined with data obtained from target predictions and IPA analyses, it could be noted that 3 transcription factors (PBX1, PPARγ and HIF1α) that were predicted as targets for the differentially expressed miRNAs had also binding sites in the upstream region of the same miRNAs (Table S4). PBX1 (pre-B-cell leukemia homeobox) is a potential target for miR-101 that was over 5-fold upregulated in chondroblasts. PPARγ (peroxisome proliferator-activated receptor γ) is potentially regulated by miR-130b that came up when comparing the miRNA expression between osteoblasts and chondroblasts. HIF1α was predicted to be regulated by miR-199a that was over 12-fold upregulated in chondroblasts.

When the interplay of miRNA target genes and TFs was analysed on the basis of published observations, a global regulatory network could be observed (Fig. 4B). miR-199a that was upregulated in chondroblasts was found to target HIF1α. miR-124a was also upregulated in chondroblasts with RFX1 as its target. miR-96 expression was strongly suppressed in chondroblasts, and it was found to target SOX5, a transcription factor that controls chondrogenesis. The previously described cartilage-specific miR-140 target HDAC4⁴¹ could also be explained by the regulatory network presented in Figure 4, yet the expression of miR-140 remained constant in our experiments. When the existence of genomic miRNA clusters were analysed for the miRNAs studied here, two miRNAs were found to be located in clusters, namely miR-199a/miR-214 and miR-96/miR-182/miR-183 . These clusters may strengthen the regulation of Sox5-Sox6 axis leading to type II collagen responses, or of HIF1α-PGF-axis leading to various hypoxia responses (Fig. 4).

Discussion

Although miRNAs have been shown to play an important role in cell differentiation, their contribution to osteo- or chondrogenesis has not been previously demonstrated. We present here multiple lines of evidence to suggest that miRNAs are an integral part of the transcription factor network regulating bone marrow stem cell differentiation and proliferation. Promoter analyses and target predictions carried out with differentially expressed miRNAs show that transcription factors may act as direct miRNA targets, but also as upstream regulators of miRNA target genes, or the miRNAs themselves. This way they constitute loops that strengthen or attenuate regulatory events.

Since there was no previous data about miRNA expression in mesenchymal stromal cells, we selected the miRNAs for this experiment based on their expression in haematopoietic tissue. We hypothesised that the haematopoietic miRNAs may also target genes involved in the differentiation of mesenchymal tissue within the bone marrow microenvironment. Haematopoietic miRNAs were supplemented with miRNAs found from mouse embryonic stem cells and with miRNAs selected based on target prediction studies carried out with genes involved in bone or cartilage function. In addition, miR-140 was included based on its suggested role in chondrogenesis.⁴¹

To gain more insight into the cellular functions possibly affected by the studied miRNAs, the predicted target genes of the differently expressed miRNAs were analysed. We hypothesised that if miRNA target genes are physiologically relevant, they should produce significant interaction networks in the pathway analysis. Extremely low p-values and extensive interactions between genes predicted as miRNA targets were observed, which cannot be explained by chance alone. Bioinformatics’ approaches have suggested that miRNA expression may be regulated by transcription factors.⁴² When the interaction networks between miRNAs and transcription factors were computationally studied, thousands of human genes were suggested to be regulated by miRNA-TF interactions.⁴³ As miRNAs are known to target many TFs, the regulatory network appears to be very complex. Composite loops⁴³ were observed in our analysis for 3 TFs: PBX1, PPARγ and HIF1α. These harboured binding sites in the upstream region of miRNAs while they were also predicted as target genes for the same miRNA. In addition, numerous TF binding sites were observed in the upstream regions of the differentially expressed miRNAs.

It became obvious from these analyses that only a fraction, if any, of the miRNA responses are such that lead to a single, easily predictable outcome in cells. Instead, accumulation of several miRNA responses eventually may lead to physiological responses (Fig. 5). An example of such a response is the downregulation of Type I collagen expression by miR-124a. The TF RFX1 (Regulatory factor X1) is a target of miR-124a, which was upregulated in chondroblasts. Another example is Sox5, which regulates Type II collagen expression during chondrogenesis but also lipase expression, raising the possibility that Sox5 has also a role in adipogenesis and lineage commitment. In MSC cultures, β-glycerophosphate and ascorbic acid induce osteoblast formation and lead to miR-96 upregulation. During osteoblastic differentiation, adipocyte or chondrocyte genes, like Sox5-regulated lipase or Type II collagen, must be turned off e.g. by miR-96 in favour of osteogenic differentiation. The complex regulatory networks described here indicate that miRNA-TF interactions are powerful modulators of bone marrow stromal cell differentiation. The specific biological functions related to the target genes of differentially expressed miRNAs suggest that they are involved in developmental pathways including cellular development and movement, cell morphology, cell signalling, cell death, and connective tissue development and function.

Figure 5

Model for miRNA-TF interactions during the differentiation of bone marrow stromal cells into bone-forming osteoblasts or chondroblasts.

Footnotes

Acknowledgments

We thank Dr K. G. Büki for the critical reading of the manuscript, Dr H. Vuorikoski for kind support with IPA analyses and Mrs T. Oivanen for technical assistance. This work was supported by TEKES, the Research Council for Health of the Academy of Finland (projects 116307, 203446, and 205581), the National Drug Discovery Graduate School and Turku Graduate School for Biomedical Sciences.

Supplementary Materials

Table S4.

Biological roles for TFs with binding sites upstream of miRNA target genes. The highlighted TFs represent composite loops where the TF binds to a miRNA but is also a target for the same miRNA. All interactions are supported by published data, with key references listed below.

miRNA(s)	TF	Downstream target gene(s)	Physiological response(s)	Reference(s)
miR-101	PBX1	MyoD, COL1A1, BMP4	Suppression of myogenic differentiation in favour of osteogenesis or chondrogenesis. PBX1 is regulated e.g. by retinoic acid.	Maves et al. (2007) Development 134:3371–3382
miR-130b	PPARγ	Multiple adipocyte genes, PPARγ→UCP1 → NRF1	Suppression of adipogenic differentiation in favour of osteogenesis or chondrogenesis. PPARγ is regulated e.g. by dexamethasone and NFκB.”	Hong et al. (2005) Science 309:1074–1078
miR-199a	HIF1α	SOX9, HDAC4, RUNX2, ARNT	Regulation of cell differentiation and proliferation in hypoxia. The major regulator of HIF1a is oxygen.	Wang et al. (2007) J Clin Invest 117:1616–1626; Schipani et al. (2001) Genes Dev 15:2865–2876
miR-124a	RFX1	Type 1 collagen	Suppression of osteogenic differentiation in favour of chondrogenesis. RFX1 is regulated e.g. by interferon gamma.	Xu et al. (2006) J Biol Chem 281:9260–9270
miR-96	SOX5	Type II collagen	Supression of chondrogenic differentiation in favour of osteogenesis. SOX5 is regulated e.g. by NFkB, retinoic acid and CDC42.	Lefebvre et al. (1998) EMBO J 17:5718–5733; Smits et al. (2001) Dev Cell 1:277–290

Abbreviations: PBX, pre-B-cell leukemia homeobox 1; PPARγ, peroxisome proliferator-activated receptor γ; HIF1α, hypoxia-inducible factor 1α; RFX1, Regulatory factor X1; SOX5, sex determining region Y-box 5.

References

Bartel

2004. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell, 116: 281–97.

Lee

, Feinbaum

, and Ambros

1993. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell, 75: 843–54.

Griffiths-Jones

, Grocock

, van Dongen

2006. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res., 34: D140–4.

Reinhart

, Slack

, Basson

2000. The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature, 403: 901–6.

Lü

, Qian

, Chen

2005. Differential expression of components of the microRNA machinery during mouse organogenesis. Biochem. Biophys. Res. Commun., 334: 319–23.

Brennecke

, Hipfner

, Stark

2003. bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell, 113: 25–36.

, Vernooy

, Guo

2003. The Drosophila microRNA Mir-14 suppresses cell death and is required for normal fat metabolism. Curr. Biol, 13: 790–5.

Dostie

, Mourelatos

, Yang

2003. Numerous microRNPs in neuronal cells containing novel microRNAs. RNA, 9: 180–6.

Chen

, Li

, Lodish

2004. MicroRNAs modulate hematopoietic lineage differentiation. Science, 303: 83–6.

10.

Krichevsky

, King

, Donahue

2003. A microRNA array reveals extensive regulation of microRNAs during brain development. RNA, 9: 1274–81.

11.

Houbaviy

, Murray

, and Sharp

2003. Embryonic stem cell-specific MicroRNAs. Dev. Cell, 5: 351–8.

12.

Garzon

, Pichiorri

, Palumbo

2006. MicroRNA fingerprints during human megakaryocytopoiesis. Proc. Natl. Acad. Sci. U.S.A., 103: 5078–83.

13.

Lakshmipathy

, and Hart

2007. Concise Review: Micro RNA Expression in Multipotent Mesenchymal Stromal Cells. Stem Cells.

14.

Pittenger

, Mackay

, Beck

1999. Multilineage potential of adult human mesenchymal stem cells. Science, 284: 143–7.

15.

Kondo

, Wagers

, Manz

2003. Biology of hematopoietic stem cells and progenitors: implications for clinical application. Annu Rev. Immunol., 21: 759–806.

16.

Spangrude

, Heimfeld

, and Weissman

1988. Purification and characterization of mouse hematopoietic stem cells. Science, 241: 58–62.

17.

Zhao

, Glesne

, and Huberman

2003. A human peripheral blood monocyte-derived subset acts as pluripotent stem cells. Proc. Natl. Acad. Sci. U.S.A., 100: 2426–31.

18.

Friedenstein

, Gorskaja

, and Kulagina

1976. Fibroblast precursors in normal and irradiated mouse hematopoietic organs. Exp. Hematol, 4: 267–74.

19.

Prockop

1997. Marrow stromal cells as stem cells for nonhematopoietic tissues. Science, 276: 71–4.

20.

Lefebvre

, Behringer

, and de Crombrugghe

2001. L-Sox5, Sox6 and Sox9 control essential steps of the chondrocyte differentiation pathway. Osteoarthritis Cartilage, 9 Suppl A: S69–75.

21.

Kronenberg

2003. Developmental regulation of the growth plate. Nature, 423: 332–6.

22.

Xie

, Lu

, Kulbokas

2005. Systematic discovery of regulatory motifs in human promoters and 3’ UTRs by comparison of several mammals. Nature, 434: 338–45.

23.

Takahashi

, Yamana

, Yoshiki

1988. Osteoclast-like cell formation and its regulation by osteotropic hormones in mouse bone marrow cultures. Endocrinology, 122: 1373–82.

24.

Phinney

, Kopen

, Isaacson

1999. Plastic adherent stromal cells from the bone marrow of commonly used strains of inbred mice: variations in yield, growth, and differentiation. J. Cell Biochem., 72: 570–85.

25.

Eaves

, and Eaves

1988. Maintenance and proliferation control of primitive hemopoietic progenitors in long-term cultures of human marrow cells. Blood Cells, 14: 355–68.

26.

Peister

, Mellad

J.A.

, Larson

B.L.

2004. Adult stem cells from bone marrow (MSCs) isolated from different strains of inbred mice vary in surface epitopes, rates of proliferation, and differentiation potential. Blood, 103: 1662–8.

27.

von Kossa

Uber die im organismus kunztlich erzeugen Verkalkungen. Vol 29: Beitr Anat; 1901: 163–202.

28.

Johnstone

, Hering

, Caplan

1998. In vitro chondrogenesis of bone marrow-derived mesenchymal progenitor cells. Exp. Cell Res., 238: 265–72.

29.

Salminen

, Vuorio

, and Säämänen

2001. Expression of Sox9 and type IIA procollagen during attempted repair of articular cartilage damage in a transgenic mouse model of osteoarthritis. Arthritis Rheum., 44: 947–55.

30.

Linsenmayer

, and Hendrix

1980. Monoclonal antibodies to connective tissue macromolecules: type II collagen. Biochem. Biophys. Res. Commun., 92: 440–6.

31.

Chen

, and Lodish

2005. MicroRNAs as regulators of mammalian hematopoiesis. Semin. Immunol., 17: 155–65.

32.

Monticelli

, Ansel

, Xiao

2005. MicroRNA profiling of the murine hematopoietic system. Genome Biol, 6: R71.

33.

Livak

, and Schmittgen

2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods, 25: 402–8.

34.

Muller

, Janovjak

, Miserez

2002. Processing of gene expression data generated by quantitative real-time RT-PCR. Biotechniques, 32: 1372–1374, 1376,1378–9.

35.

Kuhn

, Martin

, Feldman

2008. Experimental validation of miRNA targets. Methods, 44: 47–54.

36.

Berezikov

, Guryev

, Plasterk

2004. CONREAL: conserved regulatory elements anchored alignment algorithm for identification of transcription factor binding sites by phylogenetic footprinting. Genome Res., 14: 170–8.

37.

Wingender

, Dietze

, Karas

1996. TRANSFAC: a database on transcription factors and their DNA binding sites. Nucleic Acids Res., 24: 238–41.

38.

Matys

, Kel-Margoulis

, Fricke

2006. TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res., 34: D108–10.

39.

Baddoo

, Hill

, Wilkinson

2003. Characterization of mesenchymal stem cells isolated from murine bone marrow by negative selection. J. Cell Biochem., 89: 1235–49.

40.

Sandelin

, Alkema

, Engström

2004. JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res., 32: D91–94.

41.

Tuddenham

, Wheeler

, Ntounia-Fousara

2006. The cartilage specific microRNA-140 targets histone deacetylase 4 in mouse cells. FEBS Lett, 580: 4214–7.

42.

Lee

, Li

, Brower-Sinning

2007. Regulatory circuit of human microRNA biogenesis. PLoS Comput. Biol, 3: e67.

43.

Shalgi

, Lieber

, Oren

2007. Global and Local Architecture of the Mammalian microRNA-Transcription Factor Regulatory Network. PLoS Comput. Biol, 3: e131.