Mechanical loading exerts a profound influence on bone density and architecture, but the exact mechanism is unknown. Our study shows that expression of the neurological transcriptional factor zinc finger of the cerebellum 1 (ZIC1) is markedly increased in trabecular bone biopsies in the lumbar spine compared with the iliac crest, skeletal sites of high and low mechanical stress, respectively. Human trabecular bone transcriptome analyses revealed a strong association between ZIC1 mRNA levels and gene transcripts characteristically associated with osteoblasts, osteocytes and osteoclasts. This supposition is supported by higher ZIC1 expression in iliac bone biopsies from postmenopausal women with osteoporosis compared with age-matched control subjects, as well as strongly significant inverse correlation between ZIC1 mRNA levels and BMI-adjusted bone mineral density (BMD) (Z-score). ZIC1 promoter methylation was decreased in mechanically loaded vertebral bone compared to unloaded normal iliac bone, and its mRNA levels correlated inversely with ZIC1 promoter methylation, thus linking mechanical stress to epigenetic control of gene expression. The findings were corroborated in cultures of rat osteoblast progenitors and osteoblast-like cells. This study demonstrates for the first time how skeletal epigenetic changes that are affected by mechanical forces give rise to marked alteration in bone cell transcriptional activity and translate to human bone pathophysiology.

This study explores how select microRNAs (miRNAs) influence bone structure in humans and in transgenic mice. In trabecular bone biopsies from 84 postmenopausal women (healthy, osteopenic, and osteoporotic), we demonstrate that DLEU2 (deleted in lymphocytic leukemia 2)-encoded miR-15a-5p is strongly positively associated with bone mineral density (BMD) at different skeletal sites. In bone transcriptome analyses, miR-15a-5p levels correlated positively with the osteocyte characteristic transcripts SOST (encoding sclerostin) and MEPE (Matrix Extracellular Phosphoglycoprotein), while the related miR-15b-5p showed a negative association with BMD and osteoblast markers. The data imply that these miRNAs have opposite roles in bone remodeling with distinct actions on bone cells. Expression quantitative trait loci (eQTL) variants confirmed earlier DLEU2 associations. Furthermore, a novel variant (rs12585295) showed high localization with transcriptionally active chromatin states in osteoblast primary cell cultures. The supposition that DLEU2-encoded miRNAs have an important regulatory role in bone remodeling was further confirmed in a transgenic mice model showing that miR-15a/16-1-deleted mice had significantly higher percentage bone volume and trabecular number than the wild type, possibly due to prenatal actions. However, the three-point mechanical break force test of mice femurs showed a positive correlation between strength and miR-15a-5p/miR-16-5p levels, indicating differential effects on cortical and trabecular bone. Moreover, these miRNAs appear to have distinct and complex actions in mice prenatally and in adult humans, impacting BMD and microstructure by regulating bone cell transcription. However, detailed interactions between these miRNAs and their downstream mechanisms in health and disease need further clarification.