Targeting miRNAs in Osteoblast Differentiation under Malnutrition Conditions

Aims: Previous studies reported that reduced bone formation was identified in fasting adult female mice compared with the ad libitum control group. An increasing number of studies have shown that miRNAs contribute to bone homeostasis. Unfortunately, there are minor concerns about the underlying mechanisms in osteoblastic differentiation under malnutrition conditions. Methods: We investigated microRNAs (miRNAs) in osteoblastic differentiation under malnutrition conditions using high-throughput bioinformatics approaches. To screen for targeted microRNAs, sequences were quantified by aligning reads to miRbase using miRDeep2 software. Unadjusted p-values were calculated using the Student’s t-test. Genes with a p-value of <0.05 and log 2FC (fold change) ≥ 1 were considered differentially expressed genes (DEGs). DEGs were submitted to Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses, respectively. Results: They were mainly enriched in biological process terms type and biological pathways, respectively. Particularly, we evaluated seven microRNAs, mir-494 3p, mir-466, mir-455, mir-708, mir-298, mir-92 and mir-224, which likely play roles in osteoblastogenesis in fasting adult mice. Conclusion: To our knowledge, this is the first study on the expression pattern of miRNA in osteoblasts of malnourished adult mice. These targeting miRNAs may provide a potential therapeutic approach to treat osteoporosis.


Introduction
The maintenance of total bone mass under normal conditions is a balance of os-How to cite this paper: Wang, Y.Y., He, R. and Zhong, L.J. (2018) Targeting miRNAs in Osteoblast Differentiation under Malnu-teogenesis by osteoblasts and resorption by osteoclasts.The loss of bone mass under pathological conditions is primarily caused by the failure of bone turnover, which includes decreased bone formation and enhanced bone resorption.
Energy intake influences bone turnover in anorexia nervosa patients, which are commonly young exercising women or elderly underweight women [1] [2] [3].
Interestingly, bone formation and non-bone resorption were significantly changed in adult female mice between the ad libitum and the fasting groups [6].
The osteoblast lineage arises from mesenchymal stromal cells, which are regulated by genetic and epigenetic mechanisms.Emerging evidence indicates that microRNAs (miRNAs) play crucial roles in this process, from embryonic skeletal development to the maintenance of adult bone tissue [7] [8] [9].The small non-coding RNAs, miRNAs, post-transcriptionally regulate gene expression.
Unfortunately, there is no information about the regulatory mechanisms of miRNAs in osteoblastic differentiation under malnutrition conditions.To gain further insight into these molecular mechanisms, we carried out miRNA-seq using osteoblast-enriched cells derived from lower limbs.Differentially expressed genes (DEGs) in the AL and FA groups were screened.Subsequently, the underlying functions of these two sets of DEGs were investigated by gene ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses.In total, we identified 103 DEGs, including 68 upand 35 down-regulated genes.Among them, 6 had high potential for functioning in osteoblastogenesis based on previous studies.Here, we report that these miRNAs may play important roles in bone formation.

Animals
Before the study, all experiments were reviewed and approved by the Animal Care and Use Committee of Hangzhou Normal University, Zhejiang, China.
Twelve week-old C57BL/6J female mice were housed in the Hangzhou Normal University Animal Center.Mice were maintained on a 12/12-hour light/dark cycle under controlled temperature and humidity conditions.Mice were randomly divided into two groups (n = 6/group): the food ad libitum control group and the experimental group that underwent complete fasting with water ad libitum.Fasting duration (3 days) was chosen based on the experiment conducted in previous researches [4] [10].At the start of the study, mice in the fasting group were transferred to a clean cage without access to food at Zeitgeber time 20 (ZT20).These mice were subsequently sacrificed at ZT8 on the third day of fasting [10].

High-Throughput microRNA Sequencing
To separate osteoblasts from the whole bone, muscle, connective tissue and periosteum were removed from femurs.Bone marrow cells in the diaphyses of fe-murs, which were cut at the metaphyses, were flushed out using phosphate-buffered saline.Osteoblast-enriched cells were collected using a micro-interdental brush (Kobayashi Pharmaceutical Co. Ltd., Osaka, Japan).
Total RNA was extracted using a Qiagen miRNeasy Mini Kit.The integrity of RNA was verified by running it on a 1% agarose gel and staining with ethidium bromide.The concentration of RNA was measured using a Nanodrop 2000 (Thermo, USA).Purified RNAs were used to construct small RNA libraries according to Illumina's recommendations.Briefly, purified sRNAs were ligated with 3' and 5' adapters (Illumina, San Diego, CA, USA).The first strand of cDNA was synthesised using reverse transcription, and the synthesised cDNAs were subjected to a PCR amplification process of 15 cycles, each consisting of the basic steps of denaturation (10 s at 98˚C), annealing (30 s at 60˚C) and extension (15 s at 72˚C).The amplified products (approximately 140 bp) were purified by 6% polyacrylamide gel electrophoresis.Following the purification of amplified cDNAs, products were sequenced on a Hiseq2000 (Illumina, San Diego, CA, USA) using the Illumina paired-end/single-end RNA-seq approach.

Screening and Functional annotation for DEGs
The sequences were quantified by aligning reads to miRbase (version 21, http://www.mirbase.org)using the miRDeep2 package [11].The expression of miRNAs was normalised by transcripts per million (read count * 10 6 /total transcript count).Analysis of the differential expression of miRNA was performed using Bioconductor's package DEGs.P-values were estimated using a negative binomial distribution model and local regression to estimate the relationship between the dispersion and the mean of each miRNA.Raw p-values were considered statistically significant for p < 0.05, were adjusted for multiple testing using the Benjamini and Hochberg procedure with a 10% false discovery rate threshold.A fold change (FC) threshold of 2 (or log 2FC > 1) was used to select modulated miRNAs.
Functional enrichment of target genes was assessed based on the biological process category in GO [12] and KEGG annotation terms [13].GO and KEGG signalling pathway analyses were performed using the GO Function package (version 1.14.0) in Bioconductor (http://www.bioconductor.org/packages/release/bioc/html/GOFunction.html)[14], which conducted the standard hypergeometric test.A p-value of <0.05 was considered to indicate a statistically significant difference.

Screening for DEGs
Between samples from the fasting and control groups, we identified 103 DEGs with a p-value < 0.05 and log 2FC ≥1, including 68 up-and 35 down-regulated genes.Some microRNAs identified were clearly osteoblastogenesis-related genes [15]- [21].Among others, seven were likely to have a role in the process of os-teoblastogenesis.We show detailed information on them in Table 1.

Functional and Pathway Analyses of DEGs
According to the GO and KEGG databases, we analysed the functions and pathways of all the predicted target genes from differentially expressed miRNAs.
Significant functions and pathways were filtered based on p-values < 0.05.These genes were classified into functional categories based on GO predictions for biological processes and molecular function.As shown in Figure 1, the main GO categories for these genes were biological regulation, binding, single-organism cellular process, regulation of biological process, cellular response to stimulus and metabolic processes.In addition, target genes were enriched in the MAPK signalling pathway, axon guidance, the Wnt signalling pathway, focal adhesion, endocytosis, the T cell receptor signalling pathway and pathways involved in cancer (Figure 2).

Discussions
A class of short noncoding single stranded RNAs, miRNAs, have emerged as essential posttranscriptional regulators of gene expression.An increasing number of studies have shown that miRNAs contribute to bone homeostasis, which suggests that targeting miRNAs may provide a therapeutic approach to treat Table 1.Typical DEGs screened between the AL and FA groups.Inhibits proliferation of MC3T3E1 cells [22]; inhibits proliferation and metastasis of osteosarcoma [23]; induces cellular senescence of oral squamous carcinoma cells [24].
GO analysis for putative target genes.A p-value of <0.05 and a false discovery rate of <0.05 were used as thresholds to select significant GO categories; lgP is the negative logarithm of the p-value.
KEGG pathway analysis for putative target genes.A p-value of <0.05 and a false discovery rate of <0.05 were used as the threshold to select significant KEGG pathways; lgP is the negative logarithm of the p-value.

KEGG pathway
osteoporosis [37].Unfortunately, there are minor concerns about the underlying mechanisms in osteoblastic differentiation under malnutrition conditions.We investigated miRNAs identified in osteoblast-enriched cells analysed using high-throughput bioinformatics approaches.
To our knowledge, this is the first study to examine miRNA expression patterns in osteoblasts of adult mice under malnutrition conditions.However, the detailed functions associated with these DEGs in osteoblastic differentiation require further investigation, and the target genes of these miRNAs need experimental validation.