MiR-133a in Human Circulating Monocytes: A Potential Biomarker Associated with Postmenopausal Osteoporosis

Abstract

Background

Osteoporosis mainly occurs in postmenopausal women, which is characterized by low bone mineral density (BMD) due to unbalanced bone resorption by osteoclasts and formation by osteoblasts. Circulating monocytes play important roles in osteoclastogenesis by acting as osteoclast precursors and secreting osteoclastogenic factors, such as IL-1, IL-6 and TNF-α. MicroRNAs (miRNAs) have been implicated as important biomarkers in various diseases. The present study aimed to find significant miRNA biomarkers in human circulating monocytes underlying postmenopausal osteoporosis.

Methodology/Principal Findings

We used ABI TaqMan® miRNA array followed by qRT-PCR validation in circulating monocytes to identify miRNA biomarkers in 10 high and 10 low BMD postmenopausal Caucasian women. MiR-133a was upregulated (P=0.007) in the low compared with the high BMD groups in the array analyses, which was also validated by qRT-PCR (P=0.044). We performed bioinformatic target gene analysis and found three potential osteoclast-related target genes, CXCL11, CXCR3 and SLC39A1. In addition, we performed Pearson correlation analyses between the expression levels of miR-133a and the three potential target genes in the 20 postmenopausal women. We did find negative correlations between miR-133a and all the three genes though not significant.

Conclusions/Significance

This is the first in vivo miRNA expression analysis in human circulating monocytes to identify novel miRNA biomarkers underlying postmenopausal osteoporosis. Our results suggest that miR-133a in circulating monocytes is a potential biomarker for postmenopausal osteoporosis.

About the Authors:

Yang Wang

Affiliation: Osteoporosis Research Center, School of Medicine, Creighton University, Omaha, Nebraska, United States of America

Ling Li

Affiliation: Osteoporosis Research Center, School of Medicine, Creighton University, Omaha, Nebraska, United States of America

Benjamin T. Moore

Affiliation: Osteoporosis Research Center, School of Medicine, Creighton University, Omaha, Nebraska, United States of America

Xian-Hao Peng

Affiliation: Osteoporosis Research Center, School of Medicine, Creighton University, Omaha, Nebraska, United States of America

Xiang Fang

Affiliation: Biostatistical Core, Office of Research and Compliance, Creighton University, Omaha, Nebraska, United States of America

Joan M. Lappe

Affiliation: Osteoporosis Research Center, School of Medicine, Creighton University, Omaha, Nebraska, United States of America

Robert R. Recker

Affiliation: Osteoporosis Research Center, School of Medicine, Creighton University, Omaha, Nebraska, United States of America

Peng Xiao

* E-mail: [email protected]

Affiliation: Osteoporosis Research Center, School of Medicine, Creighton University, Omaha, Nebraska, United States of America

Introduction

Tissue-specific expressed microRNAs (miRNAs) are short non-coding RNA molecules that regulate gene expression, generally by destabilizing mRNAs or suppressing translation. MiRNAs have been identified as important biomarkers and regulators in various human diseases such as cancer [1], [2], diabetes [3], [4] and myocardial disease [5]. In the bone area, many miRNAs regulate osteoblastogenesis [6]–[16]. However, very few miRNAs have been related to osteoclastogenesis [17], [18]. MiRNA miR-223 played an essential role in osteoclastogenesis in a mouse osteoclast precursor cell line [18]. MiR-146a inhibited osteoclastogenesis from human circulating mononuclear cells [17].

Circulating monocytes are important cells that participate in osteoclatogenesis by acting as osteoclast precursors [19]–[22] and secreting osteoclastogenesis-related factors, such as IL-1 (interleukin-1), IL-6 and TNF-α (tumor necrosis factor-alpha) [23]–[25]. In addition, human studies have found associations of gene expression levels in circulating monocytes and osteoporosis, such as ANXA2 (annexin A2) [26], STAT1 (signal transducer and activator of transcription 1) [27], CCR3 [chemokine (C-C motif) receptor 3], HDC (histidine decarboxylase), and GCR (glucocorticoid receptor) [28].

However, no study has been conducted to identify miRNA biomarkers in circulating monocytes associated with human osteoporosis in vivo. Our present study aimed to identify differentially expressed miRNAs in circulating monocytes isolated from postmenopausal Caucasian women with discordant bone mineral density (BMD) using ABI miRNA array technology followed by qRT-PCR (quantitative RT-PCR). We found the significance of miR-133a in human circulating monocytes associated with postmenopausal osteoporosis. Further bioinformatic analysis of miR-133a identified its potential target genes that may be important in osteoclastogenesis.

Materials and Methods

Human subjects and characteristics

The study was approved by the Institutional Review Board at Creighton University, and all the subjects signed informed-consent documents before entering the project. All the subjects were Caucasians of European origin recruited from the vicinity of Creighton University in Omaha, NE. The exclusion criteria were detailed in our previous mRNA expression profiling study on B cells isolated from postmenopausal Caucasians for different BMD status [29]. The information such as age, ethnicity, menstrual status, medication history, and disease history was obtained via questionnaire. We recruited 20 unrelated postmenopausal Caucasian women, 10 with high BMD (spine or hip Z-score>0.84) and 10 with low BMD (spine or hip Z-score<−0.84). The high and low BMD groups are the top and bottom 20% BMD distributions of the age-, sex- and ethnicity-matched population. BMD (g/cm²) for the lumbar spine (L1-4) and total hip (femoral neck, trochanter, and intertrochanteric region) were measured by Hologic 4500A dual energy X-ray absorptiometry (DXA) scanners (Hologic Inc., Bedford, MA). The machine was calibrated daily. The measurement precision as reflected by the coefficient of variation (CV) was 0.9% and 1.4% for spine and hip BMD, respectively. Postmenopausal status was defined as the date of the last menses followed by at least 12 months of no menses. All the study subjects were aged 57–68. The detailed characteristics of the study subjects are summarized in Table 1.

[Figure omitted. See PDF.]

Table 1. Characteristics of the study subjects.

https://doi.org/10.1371/journal.pone.0034641.t001

As shown in Table 1, both hip and spine BMD were significantly different between the high and low BMD groups. For age, weight and height traits, only height showed marginal difference between the two BMD groups. However, height only demonstrated a very small effect on quantitative BMD variations [30]. Moreover, in this study, BMD was classified as a quality trait into two categories, the low and the high BMD. Therefore, the effect of height on BMD can be ignored in this study.

Monocyte isolation

Blood mononuclear cells (MNCs) from 70 ml peripheral blood from each study subject were separated by density gradients with UNI-SEP tubes containing a solution of 5.6% polysucrose and 9.6% sodium metrizoate with a density of 1.077 g/ml (Novamed, Jerusalem, Israel). Monocytes were isolated by a negative isolation kit, Dynabeads® Untouched™ Human Monocytes (Dynal Biotech, Lake Success, NY, USA), which contains a cocktail of CD2, CD7, CD16, CD19, CD56 and CD235a antibodies to deplete T cells, B cells, natural killer cells, erythrocytes and granulocytes, leaving monocytes naive and free of the surface-bound antibody and beads. The purity of isolated monocytes was assessed by flow cytometry with fluorescence labeled antibodies CD19-PE and CD45-FITC (BD Biosciences, San Jose, CA USA), and the average purity is about 85% with 3% deviation.

Total RNA extraction

The mirVana miRNA Isolation Kit (Ambion, Austin, Texas, USA) was used to extract total RNA including miRNAs from each cell sample following the manufacturer's protocol. Total RNA concentration and integrity were evaluated by an Agilent 2100 Bioanalyzer (Agilent, Palo Alto, CA, USA). Each RNA sample has a high quality with an excellent integrity number >9.0.

MiRNA array procedures

We used TaqMan® Human MicroRNA Array v1.0 (Applied Biosystems, Foster City, CA, USA) to perform miRNA expression profiling for each RNA sample. Each array covers 365 human miRNAs and endogenous controls RNU48 and RNU44. First, TaqMan miRNA Multiplex Reverse Transcription Kit (Applied Biosystems) was used for the RT reaction. For each RNA sample, the RT reaction was performed in a 63 µl reaction system including 1.8 µl 100 mM dNTPs, 18 µl Reverse Transcriptase (50 U/ml), 9 µl 10× RT Buffer, 1.13 µl RNase Inhibitor (20 U/µl), 16 µl sample RNA, and 17.08 µl nuclease-free water. The reaction conditions were as follows: 30 min at 16°C, 30 min at 42°C, and 5 min at 85°C. After that, we mixed 450 µl diluted RT reaction product (diluted 62.5-fold) with 450 µl TaqMan Universal PCR Master Mix (ABI) and loaded 100 µl real-time PCR reaction mix into each port of the array card (8 ports/card). The real-time qRT-PCR for each array was carried out on an Applied BioSystems 7900HT Fast Real-time PCR System with the following reaction conditions: 2 min at 50°C, 10 min at 95°C, 40 cycles of 15 sec at 95°C plus 1 min at 60°C. For each array card, there was only one probe for each target miRNA.

In the miRNA array data analysis, the raw expression level was determined by the cycle number at which the reaction crossed a predetermined cycle threshold (CT) as identified for each miRNA probe. The relative quantity (RQ) of each miRNA for each sample is determined by 2^−ΔΔCT, where ΔCT=(CT_{Target miRNA}−CT_{endogenous control RNU48}) and ΔΔCT=(ΔCT−average ΔCT of all the samples). The RQ data were used for student's t test to identify differentially expressed miRNAs between the high and the low BMD groups.

qRT-PCR for miRNAs

To correct for the multiple-testing comparison and eliminate false positive results in the miRNA array analysis, we conducted qRT-PCR among the same 20 RNA samples to further validate the identified significant miRNAs in the array analysis. Two-step qRT-PCR was used to confirm the differentially expressed miRNAs. The first step is RT of cDNA and the second step is real-time quantitative PCR. All the reagents are provided by Applied Biosystems. The RT reaction was performed in a 15 µl volume, containing 1.5 µl Taqman RT Buffer (10×), 0.15 µl 100 mM dNTPs (100 mM), 1.0 µl Reverse Transcriptase, 0.19 µl RNase inhibitor (20 U/µl), 3.0 µl specific miRNA primer, 100 ng total RNA, and nuclease-free water to make the final volume 15 µl. The real-time quantitative PCR was performed in a 20 µl reaction volume using standard protocols on the Applied Biosystems 7900HT System. Briefly, 2.5 µl cDNA was mixed with 10.0 µl TaqMan universal PCR master mix (2×), 1.0 µl TaqMan miRNA assay and 6.5 µl nuclease-free water. The reaction conditions were the same as the above real-time PCR in the array experiments. For each RNA sample, the target miRNA and RNU48 reactions were run as triplicates in the same plate. The RQ of each miRNA for each sample is determined by 2^−ΔΔCT, where ΔCT=(average of triplicate CT_{Target miRNA}−average of triplicate CT_{endogenous control RNU48}) and ΔΔCT=(ΔCT−average ΔCT of all the samples). The RQ data were used for student's t test between the two groups.

Target gene prediction and verification

We conducted bioinformatic sequence analysis of each significant miRNA to identify potential target genes [31]. MiRNAs normally repress gene expression by base pairing at complementarity sites mainly but not exclusively in the 3′-untraslated region (3′-UTR) of the target mRNAs [32], [33]. The currently available miRNA target gene databases are all limited in the 3′-UTR analyses. We used both miRDB (http://www.miRDB.org/) and TargetScan (http://www.targetscan.org/) databases to predict target genes by searching for the presence of conserved 8-mer and 7-mer sites in their 3′-UTRs that match the seed region of each significant miRNA [34]. In addition, we also conducted qRT-PCR for the potential target genes of the significant miRNA among the same 20 RNA samples. Similar to miRNA qRT-PCR, the mRNA qRT-PCR was also composed of RT and real-time qPCR. The first step is RT of cDNA and the second step is real-time quantitative PCR. The RT and qPCR were in 100 µl and 25 µl volumes, respectively, following the company's standard protocols (Applied Biosystems). For each RNA sample, the target mRNA and internal control β-actin were run as triplicates in the same plate. We used the same calculation for RQ 2^−ΔΔCT as we did for miRNA qRT-PCR and performed student's t test between the two groups.

Results

MiRNA array analyses

Among the 365 miRNAs in the array, the expression of many miRNAs were missing among the 20 study samples, probably due to tissue-specific expression or extremely low expression. To obtain enough power, we selected miRNAs that were expressed in at least 5 samples in each BMD group for the analyses. According to this criterion, 156 qualified miRNAs (Table S1) were subject to the statistical analyses and two miRNAs, miR-133a and miR-382, showed significant upregulation in the low BMD group compared with the high BMD group (Figure 1). Specifically, miR-133a displayed a fold change of 6.48 between the low and high BMD groups as mean ± SD (4.21±2.15 vs. 0.65±0.75, P=0.007), and miR-382 showed a fold change of 3.65 between the low and high BMD groups (2.74±2.18 vs. 0.75±0.63, P=0.027).

[Figure omitted. See PDF.]

Figure 1. Expression levels (2^−ΔΔCT) of significant miRNAs measured by array analysis in circulating monocytes in the low and high BMD groups (**: P<0.01; *: P<0.05).

https://doi.org/10.1371/journal.pone.0034641.g001

qRT-PCR for miRNAs

We further performed qRT-PCR to validate the differential expression of miR-133a and miR-382. However, only the upregulation of miR-133a in monocytes in the low vs. the high BMD group (2.21±2.08 vs. 0.76±0.37) was validated by qRT-PCR (P=0.044). The difference in expression of miR-382 in monocytes in the low vs. the high BMD group (6.56±2.84 vs. 7.93±9.73) was not significant (P=0.67) (Figure 2).

[Figure omitted. See PDF.]

Figure 2. Expression levels (2^−ΔΔCT) of miR-133a and miR-382 measured by qRT-PCR analysis in circulating monocytes in the low and high BMD groups (*: P<0.05).

https://doi.org/10.1371/journal.pone.0034641.g002

Target gene prediction and verification

Our bioinformatic sequence analyses identified 226 potential target genes that are overlapped in both miRDB and TargetScan databases. By searching relevant references for all the 226 genes, we found that three of them are related to osteoclastogenesis, CXCL11 [chemokine (C-X-C motif) ligand 11], CXCR3 [chemokine (C-X-C motif) receptor 3], and SLC39A1 [solute carrier family (zinc transporter), member 1]. Table 2 demonstrates the specific putative binding sites of miR-133a in the 3′ UTRs of the three genes. We conducted qRT-PCR analyses for all the three genes among the same 20 study samples and did not find significant differential expression. In addition, we performed correlation analysis of the expression levels of miR-133a and each gene. All three genes did demonstrate negative correlation with miR-133a, although they were not significant (P>0.05) (Table 3).

[Figure omitted. See PDF.]

Table 2. Putative binding sites of miR-133a in predicted target genes in humans.

https://doi.org/10.1371/journal.pone.0034641.t002

[Figure omitted. See PDF.]

Table 3. Correlation of expression levels (2^−ΔΔCT) of miR-133a with those of three potential target genes and the ratio of expression levels of each gene in the high and low BMD groups, as measured by qRT-PCR.

https://doi.org/10.1371/journal.pone.0034641.t003

Discussion

In this study, we aimed to identify important miRNAs in human circulating monocytes associated with discordant BMD status in postmenopausal Caucasian women. We found significant upregulation of miR-133a in the low BMD group in both the array and the qRT-PCR analyses.

Human mature miR-133a is encoded by two genes: MIR133A1 for miR-133a1 at 18q11.2 (194,036,59–194,077,46 bp) and MIR133A2 for miR-133a2 at 20q13.33 (611,601,19–611,642,20 bp). Both genes encode different pre-mature miRNAs but generate the same mature miR-133a sequence. Interestingly, human genetic studies also found the association of 18q11.2 to osteoporosis-related traits [35] and linkage of 20q13 to bone phenotypes [36]–[38]. In humans, there are two types of miR-133 miRNA isoforms, miR-133a and miR-133b, with one base difference (g-a) in the last nucleotide at the 3′ end (miRBase: http://www.mirbase.org/). The ABI miRNA array used in this study includes both miR-133a and miR-133b probes. Interestingly, miR-133b is marginally unregulated in the low vs. the high BMD groups (1.51±0.67 vs. 0.95±0.63, P=0.08). In addition, we detected miR-133a expression levels in circulating B cells from the same 20 high or low BMD postmenopausal women. Circulating B cells were isolated by Dynabeads® CD19 (Pan B) (Dynal Biotech). However, miR-133a was not differentially expressed in B cells between the high and the low BMD groups (P=0.49). Therefore, miR-133a is most likely to be a monocyte specific biomarker underlying postmenopausal osteoporosis.

Many studies demonstrated that miR-133 and 133a are important in the development of muscle, such as skeletal muscle [39]–[43], and heart/cardiovascular muscle [44]–[47]. In bone, particularly, miR-133 and 133a have been found to regulate osteoblastogenesis by targeting and regulating Runx2 expression [6], [48]. A recent study also demonstrated that miR-133a was upregulated in osteoblast-like periodontal ligament stem cells treated with ibandronate, a nitrogen-containing bisphosphonate that inhibits bone resorption and is widely used to treat osteoporosis [49]. However, our study for the first time showed the association of miR-133a expression levels in circulating monocytes, the osteoclast precursors, with postmenopausal osteoporosis.

To further predict what genes are targeted and regulated by miR-133a in monocytes in bone metabolism, we used two miRNA target gene predicting databases (miRDB and TargetScan) [50]–[55]. Our bioinformatic sequence analyses and reference searching identified three potential target genes of miR-133a related to the inhibition of osteoclastogenesis, which are CXCL11, CXCR3, and SLC39A1 (Table 2). CXCL11 is a small cytokine of the CXC chemokine family. CXCL11 has been shown to inhibit osteoclast differentiation of CD14+ monocytes [56]. CXCR3 is a Gα_i protein-coupled receptor in the CXC chemokine receptor family. CXCR3 expression significantly decreased during osteoclast differentiation [57]. The SLC39A1 gene encodes zinc transporter 1 (ZIP1). Zinc deficiency has been correlated with reduction of bone growth and development of osteoporosis [58], [59]. SLC39A1 has been detected in osteoclasts and inhibited osteoclastogenesis and osteoclast function through zinc uptake [60].

All three genes did show negative correlation with miR-133a, though not significant (Table 3). Since one single miRNA normally regulates the expression of hundreds of genes, the regulatory effect of each gene may be small. Osteoporosis is a complex disease and regulated by multiple genes [61]. To investigate the combined effects, we also utilized a principle component analysis (PCA) to analyze the qRT-PCR observations from all three potential target genes. The PCA was performed using the correlation matrix. The analysis reduced the original data into two principle components (PC) that account for 58.17% and 32.78% of the total variance, respectively. A plot of the second PC against the first PC in the 10 high and 10 low BMD subjects is shown in Figure 3. We can see that both PCs are above the average levels in 6 out of 10 high BMD subjects and below the average levels in 7 out of 10 low BMD subjects, which means that the three genes are systematically upregulated in the high vs. low BMD groups. This result is consistently correlated with the downregulation of miR-133a in the high vs. low BMD groups.

[Figure omitted. See PDF.]

Figure 3. PCA of the expression levels of the three potential target genes measured by qRT-PCR in the 10 high and 10 low BMD subjects.

In each axis, “0" represents the average level of each principle component.

https://doi.org/10.1371/journal.pone.0034641.g003

There are some limitations for this study. In the miRNA array data, we performed the statistical tests for the 156 miRNAs according to expression data available in at least five samples in each group, and used raw P values but not adjusted P values for multiple tests. Actually, we did not find significant miRNAs after multiple testing adjustments by Bonferroni. However, the significant qRT-PCR P value confirmed the differential expression of miR-133a in the array analyses. The qRT-PCR validation largely solved the multiple testing problem. The sample size of 20 is relatively small. However, it is enough for an initial biomarker screening according to our previous [29], [62] and other [63], [64] gene expression profiling studies. Moreover, we are also planning to validate the miR-133a biomarker in a larger independent population. We did not find significant correlations of miR-133a and selected potential target genes. First, it may be due to the limited sample size to get enough power to detect the significant correlation. Second, many miRNAs inhibit gene expression by only suppressing protein translation instead of mRNA degradation [32], [33]. Third, there may be other unknown target genes of miR-133a in circulating monocytes that affect osteoclastogenesis. Therefore, we will test the current and potential novel target gene expression and correlation with miR-133a in a bigger population at both mRNA and protein levels in the future.

In summary, this is an initial miRNA expression study to indentify miRNA biomarkers in human circulating monocytes underlying postmenopausal osteoporosis. Our study suggested that miR-133a is a potential miRNA biomarker and/or regulatory element in circulating monocytes for postmenopausal osteoporosis.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

MiRNA array results of 156 miRNAs that were expressed in at least 5 samples in each BMD group. *: the ratio of the mean expression values between the low and high BMD groups.

https://doi.org/10.1371/journal.pone.0034641.s001

(DOC)

Acknowledgments

We thank Dr. Xiang-Ding Chen for his help with the miRNA array experiments.

Author Contributions

Conceived and designed the experiments: YW LL BTM XHP XF JML RRR PX. Performed the experiments: YW LL BTM XHP PX. Analyzed the data: YW LL BTM XF PX. Contributed reagents/materials/analysis tools: PX. Wrote the paper: YW LL BTM XHP XF JML RRR PX.

Citation: Wang Y, Li L, Moore BT, Peng X-H, Fang X, Lappe JM, et al. (2012) MiR-133a in Human Circulating Monocytes: A Potential Biomarker Associated with Postmenopausal Osteoporosis. PLoS ONE 7(4): e34641. https://doi.org/10.1371/journal.pone.0034641

References

1. Li Y, Vandenboom TG, Wang Z, Kong D, Ali S, et al. (2010) miR-146a suppresses invasion of pancreatic cancer cells. Cancer Res 70: 1486–1495.Y. LiTG VandenboomZ. WangD. KongS. Ali2010miR-146a suppresses invasion of pancreatic cancer cells.Cancer Res7014861495

2. Ma L, Teruya-Feldstein J, Weinberg RA (2007) Tumour invasion and metastasis initiated by microRNA-10b in breast cancer. Nature 449: 682–688.L. MaJ. Teruya-FeldsteinRA Weinberg2007Tumour invasion and metastasis initiated by microRNA-10b in breast cancer.Nature449682688

3. Balasubramanyam M, Aravind S, Gokulakrishnan K, Prabu P, Sathishkumar C, et al. (2011) Impaired miR-146a expression links subclinical inflammation and insulin resistance in Type 2 diabetes. Mol Cell Biochem 351: 197–205.M. BalasubramanyamS. AravindK. GokulakrishnanP. PrabuC. Sathishkumar2011Impaired miR-146a expression links subclinical inflammation and insulin resistance in Type 2 diabetes.Mol Cell Biochem351197205

4. Poy MN, Eliasson L, Krutzfeldt J, Kuwajima S, Ma X, et al. (2004) A pancreatic islet-specific microRNA regulates insulin secretion. Nature 432: 226–230.MN PoyL. EliassonJ. KrutzfeldtS. KuwajimaX. Ma2004A pancreatic islet-specific microRNA regulates insulin secretion.Nature432226230

5. Thum T, Gross C, Fiedler J, Fischer T, Kissler S, et al. (2008) MicroRNA-21 contributes to myocardial disease by stimulating MAP kinase signalling in fibroblasts. Nature 456: 980–984.T. ThumC. GrossJ. FiedlerT. FischerS. Kissler2008MicroRNA-21 contributes to myocardial disease by stimulating MAP kinase signalling in fibroblasts.Nature456980984

6. Li Z, Hassan MQ, Volinia S, van Wijnen AJ, Stein JL, et al. (2008) A microRNA signature for a BMP2-induced osteoblast lineage commitment program. Proc Natl Acad Sci U S A 105: 13906–13911.Z. LiMQ HassanS. VoliniaAJ van WijnenJL Stein2008A microRNA signature for a BMP2-induced osteoblast lineage commitment program.Proc Natl Acad Sci U S A1051390613911

7. Eskildsen T, Taipaleenmaki H, Stenvang J, Abdallah BM, Ditzel N, et al. (2011) MicroRNA-138 regulates osteogenic differentiation of human stromal (mesenchymal) stem cells in vivo. Proc Natl Acad Sci U S A 108: 6139–6144.T. EskildsenH. TaipaleenmakiJ. StenvangBM AbdallahN. Ditzel2011MicroRNA-138 regulates osteogenic differentiation of human stromal (mesenchymal) stem cells in vivo.Proc Natl Acad Sci U S A10861396144

8. Zhang J, Tu Q, Bonewald LF, He X, Stein G, et al. (2011) Effects of miR-335-5p in modulating osteogenic differentiation by specifically downregulating Wnt antagonist DKK1. J Bone Miner Res 26: 1953–1963.J. ZhangQ. TuLF BonewaldX. HeG. Stein2011Effects of miR-335-5p in modulating osteogenic differentiation by specifically downregulating Wnt antagonist DKK1.J Bone Miner Res2619531963

9. Hassan MQ, Gordon JA, Beloti MM, Croce CM, van Wijnen AJ, et al. (2010) A network connecting Runx2, SATB2, and the miR-23a∼27a∼24-2 cluster regulates the osteoblast differentiation program. Proc Natl Acad Sci U S A 107: 19879–19884.MQ HassanJA GordonMM BelotiCM CroceAJ van Wijnen2010A network connecting Runx2, SATB2, and the miR-23a∼27a∼24-2 cluster regulates the osteoblast differentiation program.Proc Natl Acad Sci U S A1071987919884

10. Kapinas K, Kessler C, Ricks T, Gronowicz G, Delany AM (2010) miR-29 modulates Wnt signaling in human osteoblasts through a positive feedback loop. J Biol Chem 285: 25221–25231.K. KapinasC. KesslerT. RicksG. GronowiczAM Delany2010miR-29 modulates Wnt signaling in human osteoblasts through a positive feedback loop.J Biol Chem2852522125231

11. Huang J, Zhao L, Xing L, Chen D (2010) MicroRNA-204 regulates Runx2 protein expression and mesenchymal progenitor cell differentiation. Stem Cells 28: 357–364.J. HuangL. ZhaoL. XingD. Chen2010MicroRNA-204 regulates Runx2 protein expression and mesenchymal progenitor cell differentiation.Stem Cells28357364

12. Mizuno Y, Yagi K, Tokuzawa Y, Kanesaki-Yatsuka Y, Suda T, et al. (2008) miR-125b inhibits osteoblastic differentiation by down-regulation of cell proliferation. Biochem Biophys Res Commun 368: 267–272.Y. MizunoK. YagiY. TokuzawaY. Kanesaki-YatsukaT. Suda2008miR-125b inhibits osteoblastic differentiation by down-regulation of cell proliferation.Biochem Biophys Res Commun368267272

13. Luzi E, Marini F, Sala SC, Tognarini I, Galli G, et al. (2008) Osteogenic differentiation of human adipose tissue-derived stem cells is modulated by the miR-26a targeting of the SMAD1 transcription factor. J Bone Miner Res 23: 287–295.E. LuziF. MariniSC SalaI. TognariniG. Galli2008Osteogenic differentiation of human adipose tissue-derived stem cells is modulated by the miR-26a targeting of the SMAD1 transcription factor.J Bone Miner Res23287295

14. Inose H, Ochi H, Kimura A, Fujita K, Xu R, et al. (2009) A microRNA regulatory mechanism of osteoblast differentiation. Proc Natl Acad Sci U S A 106: 20794–20799.H. InoseH. OchiA. KimuraK. FujitaR. Xu2009A microRNA regulatory mechanism of osteoblast differentiation.Proc Natl Acad Sci U S A1062079420799

15. Kim YJ, Bae SW, Yu SS, Bae YC, Jung JS (2009) miR-196a regulates proliferation and osteogenic differentiation in mesenchymal stem cells derived from human adipose tissue. J Bone Miner Res 24: 816–825.YJ KimSW BaeSS YuYC BaeJS Jung2009miR-196a regulates proliferation and osteogenic differentiation in mesenchymal stem cells derived from human adipose tissue.J Bone Miner Res24816825

16. Li H, Xie H, Liu W, Hu R, Huang B, et al. (2009) A novel microRNA targeting HDAC5 regulates osteoblast differentiation in mice and contributes to primary osteoporosis in humans. J Clin Invest 119: 3666–3677.H. LiH. XieW. LiuR. HuB. Huang2009A novel microRNA targeting HDAC5 regulates osteoblast differentiation in mice and contributes to primary osteoporosis in humans.J Clin Invest11936663677

17. Nakasa T, Shibuya H, Nagata Y, Niimoto T, Ochi M (2011) The inhibitory effect of microRNA-146a expression on bone destruction in collagen-induced arthritis. Arthritis Rheum 63: 1582–1590.T. NakasaH. ShibuyaY. NagataT. NiimotoM. Ochi2011The inhibitory effect of microRNA-146a expression on bone destruction in collagen-induced arthritis.Arthritis Rheum6315821590

18. Sugatani T, Hruska KA (2007) MicroRNA-223 is a key factor in osteoclast differentiation. J Cell Biochem 101: 996–999.T. SugataniKA Hruska2007MicroRNA-223 is a key factor in osteoclast differentiation.J Cell Biochem101996999

19. Fujikawa Y, Quinn JM, Sabokbar A, McGee JO, Athanasou NA (1996) The human osteoclast precursor circulates in the monocyte fraction. Endocrinology 137: 4058–4060.Y. FujikawaJM QuinnA. SabokbarJO McGeeNA Athanasou1996The human osteoclast precursor circulates in the monocyte fraction.Endocrinology13740584060

20. Higuchi S, Tabata N, Tajima M, Ito M, Tsurudome M, et al. (1998) Induction of human osteoclast-like cells by treatment of blood monocytes with anti-fusion regulatory protein-1/CD98 monoclonal antibodies. J Bone Miner Res 13: 44–49.S. HiguchiN. TabataM. TajimaM. ItoM. Tsurudome1998Induction of human osteoclast-like cells by treatment of blood monocytes with anti-fusion regulatory protein-1/CD98 monoclonal antibodies.J Bone Miner Res134449

21. Matayoshi A, Brown C, DiPersio JF, Haug J, Abu-Amer Y, et al. (1996) Human blood-mobilized hematopoietic precursors differentiate into osteoclasts in the absence of stromal cells. Proc Natl Acad Sci U S A 93: 10785–10790.A. MatayoshiC. BrownJF DiPersioJ. HaugY. Abu-Amer1996Human blood-mobilized hematopoietic precursors differentiate into osteoclasts in the absence of stromal cells.Proc Natl Acad Sci U S A931078510790

22. Purton LE, Lee MY, Torok-Storb B (1996) Normal human peripheral blood mononuclear cells mobilized with granulocyte colony-stimulating factor have increased osteoclastogenic potential compared to nonmobilized blood. Blood 87: 1802–1808.LE PurtonMY LeeB. Torok-Storb1996Normal human peripheral blood mononuclear cells mobilized with granulocyte colony-stimulating factor have increased osteoclastogenic potential compared to nonmobilized blood.Blood8718021808

23. Cohen-Solal ME, Graulet AM, Denne MA, Gueris J, Baylink D, et al. (1993) Peripheral monocyte culture supernatants of menopausal women can induce bone resorption: involvement of cytokines. J Clin Endocrinol Metab 77: 1648–1653.ME Cohen-SolalAM GrauletMA DenneJ. GuerisD. Baylink1993Peripheral monocyte culture supernatants of menopausal women can induce bone resorption: involvement of cytokines.J Clin Endocrinol Metab7716481653

24. Cohen-Solal ME, Boitte F, Bernard-Poenaru O, Denne MA, Graulet AM, et al. (1998) Increased bone resorbing activity of peripheral monocyte culture supernatants in elderly women. J Clin Endocrinol Metab 83: 1687–1690.ME Cohen-SolalF. BoitteO. Bernard-PoenaruMA DenneAM Graulet1998Increased bone resorbing activity of peripheral monocyte culture supernatants in elderly women.J Clin Endocrinol Metab8316871690

25. Pacifici R (1996) Estrogen, cytokines, and pathogenesis of postmenopausal osteoporosis. J Bone Miner Res 11: 1043–1051.R. Pacifici1996Estrogen, cytokines, and pathogenesis of postmenopausal osteoporosis.J Bone Miner Res1110431051

26. Deng FY, Lei SF, Zhang Y, Zhang YL, Zheng YP, et al. (2011) Peripheral blood monocyte-expressed ANXA2 gene is involved in pathogenesis of osteoporosis in humans. Mol Cell Proteomics 10: M111.FY DengSF LeiY. ZhangYL ZhangYP Zheng2011Peripheral blood monocyte-expressed ANXA2 gene is involved in pathogenesis of osteoporosis in humans.Mol Cell Proteomics10M111

27. Chen XD, Xiao P, Lei SF, Liu YZ, Guo YF, et al. (2010) Gene expression profiling in monocytes and SNP association suggest the importance of the STAT1 gene for osteoporosis in both Chinese and Caucasians. J Bone Miner Res 25: 339–355.XD ChenP. XiaoSF LeiYZ LiuYF Guo2010Gene expression profiling in monocytes and SNP association suggest the importance of the STAT1 gene for osteoporosis in both Chinese and Caucasians.J Bone Miner Res25339355

28. Liu YZ, Dvornyk V, Lu Y, Shen H, Lappe JM, et al. (2005) A novel pathophysiological mechanism for osteoporosis suggested by an in vivo gene expression study of circulating monocytes. J Biol Chem 280: 29011–29016.YZ LiuV. DvornykY. LuH. ShenJM Lappe2005A novel pathophysiological mechanism for osteoporosis suggested by an in vivo gene expression study of circulating monocytes.J Biol Chem2802901129016

29. Xiao P, Chen Y, Jiang H, Liu YZ, Pan F, et al. (2008) In vivo genome-wide expression study on human circulating B cells suggests a novel ESR1 and MAPK3 network for postmenopausal osteoporosis. J Bone Miner Res 23: 644–654.P. XiaoY. ChenH. JiangYZ LiuF. Pan2008In vivo genome-wide expression study on human circulating B cells suggests a novel ESR1 and MAPK3 network for postmenopausal osteoporosis.J Bone Miner Res23644654

30. Zerbini CA, Latorre MR, Jaime PC, Tanaka T, Pippa MG (2000) Bone mineral density in Brazilian men 50 years and older. Braz J Med Biol Res 33: 1429–1435.CA ZerbiniMR LatorrePC JaimeT. TanakaMG Pippa2000Bone mineral density in Brazilian men 50 years and older.Braz J Med Biol Res3314291435

31. Rajewsky N (2006) microRNA target predictions in animals. Nat Genet 38: SupplS8–13.N. Rajewsky2006microRNA target predictions in animals.Nat Genet38SupplS813

32. He L, Hannon GJ (2004) MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet 5: 522–531.L. HeGJ Hannon2004MicroRNAs: small RNAs with a big role in gene regulation.Nat Rev Genet5522531

33. Bartel DP (2004) MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116: 281–297.DP Bartel2004MicroRNAs: genomics, biogenesis, mechanism, and function.Cell116281297

34. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, et al. (2007) MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell 27: 91–105.A. GrimsonKK FarhWK JohnstonP. Garrett-EngeleLP Lim2007MicroRNA targeting specificity in mammals: determinants beyond seed pairing.Mol Cell2791105

35. Hsu YH, Zillikens MC, Wilson SG, Farber CR, Demissie S, et al. (2010) An integration of genome-wide association study and gene expression profiling to prioritize the discovery of novel susceptibility Loci for osteoporosis-related traits. PLoS Genet 6: e1000977.YH HsuMC ZillikensSG WilsonCR FarberS. Demissie2010An integration of genome-wide association study and gene expression profiling to prioritize the discovery of novel susceptibility Loci for osteoporosis-related traits.PLoS Genet6e1000977

36. Ralston SH, Galwey N, MacKay I, Albagha OM, Cardon L, et al. (2005) Loci for regulation of bone mineral density in men and women identified by genome wide linkage scan: the FAMOS study. Hum Mol Genet 14: 943–951.SH RalstonN. GalweyI. MacKayOM AlbaghaL. Cardon2005Loci for regulation of bone mineral density in men and women identified by genome wide linkage scan: the FAMOS study.Hum Mol Genet14943951

37. Deng FY, Xiao P, Lei SF, Zhang L, Yang F, et al. (2007) Bivariate whole genome linkage analysis for femoral neck geometric parameters and total body lean mass. J Bone Miner Res 22: 808–816.FY DengP. XiaoSF LeiL. ZhangF. Yang2007Bivariate whole genome linkage analysis for femoral neck geometric parameters and total body lean mass.J Bone Miner Res22808816

38. Mitchell BD, Cole SA, Bauer RL, Iturria SJ, Rodriguez EA, et al. (2000) Genes influencing variation in serum osteocalcin concentrations are linked to markers on chromosomes 16q and 20q. J Clin Endocrinol Metab 85: 1362–1366.BD MitchellSA ColeRL BauerSJ IturriaEA Rodriguez2000Genes influencing variation in serum osteocalcin concentrations are linked to markers on chromosomes 16q and 20q.J Clin Endocrinol Metab8513621366

39. Deng Z, Chen JF, Wang DZ (2011) Transgenic overexpression of miR-133a in skeletal muscle. BMC Musculoskelet Disord 12: 115.Z. DengJF ChenDZ Wang2011Transgenic overexpression of miR-133a in skeletal muscle.BMC Musculoskelet Disord12115

40. McCarthy JJ, Esser KA (2007) MicroRNA-1 and microRNA-133a expression are decreased during skeletal muscle hypertrophy. J Appl Physiol 102: 306–313.JJ McCarthyKA Esser2007MicroRNA-1 and microRNA-133a expression are decreased during skeletal muscle hypertrophy.J Appl Physiol102306313

41. Rao PK, Kumar RM, Farkhondeh M, Baskerville S, Lodish HF (2006) Myogenic factors that regulate expression of muscle-specific microRNAs. Proc Natl Acad Sci U S A 103: 8721–8726.PK RaoRM KumarM. FarkhondehS. BaskervilleHF Lodish2006Myogenic factors that regulate expression of muscle-specific microRNAs.Proc Natl Acad Sci U S A10387218726

42. Huang MB, Xu H, Xie SJ, Zhou H, Qu LH (2011) Insulin-like growth factor-1 receptor is regulated by microRNA-133 during skeletal myogenesis. PLoS One 6: e29173.MB HuangH. XuSJ XieH. ZhouLH Qu2011Insulin-like growth factor-1 receptor is regulated by microRNA-133 during skeletal myogenesis.PLoS One6e29173

43. Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, et al. (2006) The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet 38: 228–233.JF ChenEM MandelJM ThomsonQ. WuTE Callis2006The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation.Nat Genet38228233

44. Bostjancic E, Zidar N, Stajer D, Glavac D (2010) MicroRNAs miR-1, miR-133a, miR-133b and miR-208 are dysregulated in human myocardial infarction. Cardiology 115: 163–169.E. BostjancicN. ZidarD. StajerD. Glavac2010MicroRNAs miR-1, miR-133a, miR-133b and miR-208 are dysregulated in human myocardial infarction.Cardiology115163169

45. Villar AV, Merino D, Wenner M, Llano M, Cobo M, et al. (2011) Myocardial gene expression of microRNA-133a and myosin heavy and light chains, in conjunction with clinical parameters, predict regression of left ventricular hypertrophy after valve replacement in patients with aortic stenosis. Heart 97: 1132–1137.AV VillarD. MerinoM. WennerM. LlanoM. Cobo2011Myocardial gene expression of microRNA-133a and myosin heavy and light chains, in conjunction with clinical parameters, predict regression of left ventricular hypertrophy after valve replacement in patients with aortic stenosis.Heart9711321137

46. Torella D, Iaconetti C, Catalucci D, Ellison GM, Leone A, et al. (2011) MicroRNA-133 controls vascular smooth muscle cell phenotypic switch in vitro and vascular remodeling in vivo. Circ Res 109: 880–893.D. TorellaC. IaconettiD. CatalucciGM EllisonA. Leone2011MicroRNA-133 controls vascular smooth muscle cell phenotypic switch in vitro and vascular remodeling in vivo.Circ Res109880893

47. Liu N, Williams AH, Kim Y, McAnally J, Bezprozvannaya S, et al. (2007) An intragenic MEF2-dependent enhancer directs muscle-specific expression of microRNAs 1 and 133. Proc Natl Acad Sci U S A 104: 20844–20849.N. LiuAH WilliamsY. KimJ. McAnallyS. Bezprozvannaya2007An intragenic MEF2-dependent enhancer directs muscle-specific expression of microRNAs 1 and 133.Proc Natl Acad Sci U S A1042084420849

48. Zhang Y, Xie RL, Croce CM, Stein JL, Lian JB, et al. (2011) A program of microRNAs controls osteogenic lineage progression by targeting transcription factor Runx2. Proc Natl Acad Sci U S A 108: 9863–9868.Y. ZhangRL XieCM CroceJL SteinJB Lian2011A program of microRNAs controls osteogenic lineage progression by targeting transcription factor Runx2.Proc Natl Acad Sci U S A10898639868

49. Zhou Q, Zhao ZN, Cheng JT, Zhang B, Xu J, et al. (2011) Ibandronate promotes osteogenic differentiation of periodontal ligament stem cells by regulating the expression of microRNAs. Biochem Biophys Res Commun 404: 127–132.Q. ZhouZN ZhaoJT ChengB. ZhangJ. Xu2011Ibandronate promotes osteogenic differentiation of periodontal ligament stem cells by regulating the expression of microRNAs.Biochem Biophys Res Commun404127132

50. Wang X (2008) miRDB: a microRNA target prediction and functional annotation database with a wiki interface. RNA 14: 1012–1017.X. Wang2008miRDB: a microRNA target prediction and functional annotation database with a wiki interface.RNA1410121017

51. Wang X, El Naqa IM (2008) Prediction of both conserved and nonconserved microRNA targets in animals. Bioinformatics 24: 325–332.X. WangIM El Naqa2008Prediction of both conserved and nonconserved microRNA targets in animals.Bioinformatics24325332

52. Lewis BP, Burge CB, Bartel DP (2005) Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 120: 15–20.BP LewisCB BurgeDP Bartel2005Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets.Cell1201520

53. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, et al. (2007) MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell 27: 91–105.A. GrimsonKK FarhWK JohnstonP. Garrett-EngeleLP Lim2007MicroRNA targeting specificity in mammals: determinants beyond seed pairing.Mol Cell2791105

54. Garcia DM, Baek D, Shin C, Bell GW, Grimson A, et al. (2011) Weak seed-pairing stability and high target-site abundance decrease the proficiency of lsy-6 and other microRNAs. Nat Struct Mol Biol 18: 1139–1146.DM GarciaD. BaekC. ShinGW BellA. Grimson2011Weak seed-pairing stability and high target-site abundance decrease the proficiency of lsy-6 and other microRNAs.Nat Struct Mol Biol1811391146

55. Friedman RC, Farh KK, Burge CB, Bartel DP (2009) Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 19: 92–105.RC FriedmanKK FarhCB BurgeDP Bartel2009Most mammalian mRNAs are conserved targets of microRNAs.Genome Res1992105

56. Coelho LF, Magno de Freitas AG, Mennechet FJ, Blangy A, Uze G (2005) Interferon-alpha and -beta differentially regulate osteoclastogenesis: role of differential induction of chemokine CXCL11 expression. Proc Natl Acad Sci U S A 102: 11917–11922.LF CoelhoAG Magno de FreitasFJ MennechetA. BlangyG. Uze2005Interferon-alpha and -beta differentially regulate osteoclastogenesis: role of differential induction of chemokine CXCL11 expression.Proc Natl Acad Sci U S A1021191711922

57. Grassi F, Piacentini A, Cristino S, Toneguzzi S, Cavallo C, et al. (2003) Human osteoclasts express different CXC chemokines depending on cell culture substrate: molecular and immunocytochemical evidence of high levels of CXCL10 and CXCL12. Histochem Cell Biol 120: 391–400.F. GrassiA. PiacentiniS. CristinoS. ToneguzziC. Cavallo2003Human osteoclasts express different CXC chemokines depending on cell culture substrate: molecular and immunocytochemical evidence of high levels of CXCL10 and CXCL12.Histochem Cell Biol120391400

58. Eberle J, Schmidmayer S, Erben RG, Stangassinger M, Roth HP (1999) Skeletal effects of zinc deficiency in growing rats. J Trace Elem Med Biol 13: 21–26.J. EberleS. SchmidmayerRG ErbenM. StangassingerHP Roth1999Skeletal effects of zinc deficiency in growing rats.J Trace Elem Med Biol132126

59. Nishi Y (1996) Zinc and growth. J Am Coll Nutr 15: 340–344.Y. Nishi1996Zinc and growth.J Am Coll Nutr15340344

60. Khadeer MA, Sahu SN, Bai G, Abdulla S, Gupta A (2005) Expression of the zinc transporter ZIP1 in osteoclasts. Bone 37: 296–304.MA KhadeerSN SahuG. BaiS. AbdullaA. Gupta2005Expression of the zinc transporter ZIP1 in osteoclasts.Bone37296304

61. Liu YJ, Shen H, Xiao P, Xiong DH, Li LH, et al. (2006) Molecular genetic studies of gene identification for osteoporosis: a 2004 update. J Bone Miner Res 21: 1511–1535.YJ LiuH. ShenP. XiaoDH XiongLH Li2006Molecular genetic studies of gene identification for osteoporosis: a 2004 update.J Bone Miner Res2115111535

62. Chen XD, Xiao P, Lei SF, Liu YZ, Guo YF, et al. (2010) Gene expression profiling in monocytes and SNP association suggest the importance of the STAT1 gene for osteoporosis in both Chinese and Caucasians. J Bone Miner Res 25: 339–355.XD ChenP. XiaoSF LeiYZ LiuYF Guo2010Gene expression profiling in monocytes and SNP association suggest the importance of the STAT1 gene for osteoporosis in both Chinese and Caucasians.J Bone Miner Res25339355

63. Atzmon G, Yang XM, Muzumdar R, Ma XH, Gabriely I, et al. (2002) Differential gene expression between visceral and subcutaneous fat depots. Horm Metab Res 34: 622–628.G. AtzmonXM YangR. MuzumdarXH MaI. Gabriely2002Differential gene expression between visceral and subcutaneous fat depots.Horm Metab Res34622628

64. Linder K, Arner P, Flores-Morales A, Tollet-Egnell P, Norstedt G (2004) Differentially expressed genes in visceral or subcutaneous adipose tissue of obese men and women. J Lipid Res 45: 148–154.K. LinderP. ArnerA. Flores-MoralesP. Tollet-EgnellG. Norstedt2004Differentially expressed genes in visceral or subcutaneous adipose tissue of obese men and women.J Lipid Res45148154