Relationships between matrix mineralization, oxidative metabolism, and mitochondrial structure during ATDC5 murine chondroprogenitor cell line differentiation.

The mechanistic relationships between the progression of growth chondrocyte differentiation, matrix mineralization, oxidative metabolism, and mitochondria content and structure were examined in the ATDC5 murine chondroprogenitor cell line. The progression of chondrocyte differentiation was associated with a statistically significant (p ≤ 0.05) ~2-fold increase in oxidative phosphorylation. However, as matrix mineralization progressed, oxidative metabolism decreased. In the absence of mineralization, cartilage extracellular matrix mRNA expression for Col2a1, Aggrecan, and Col10a1 were statistically (p ≤ 0.05) ~2-3-fold greater than observed in mineralizing cultures. In contrast, BSP and Phex that are associated with promoting matrix mineralization showed statistically (p ≤ 0.05) higher ~2-4 expression, while FGF23 phosphate regulatory factor was significantly lower (~50%) in mineralizing cultures. Cultures induced to differentiate under both nonmineralizing and mineralizing media conditions showed statistically greater basal oxidative metabolism and ATP production. Maximal respiration and spare oxidative capacity were significantly elevated (p ≤ 0.05) in differentiated nonmineralizing cultures compared to those that mineralized. Increased oxidative metabolism was associated with both an increase in mitochondria volume per cell and mitochondria fusion, while mineralization diminished mitochondrial volume and appeared to be associated with fission. Undifferentiated and mineralized cells showed increased mitochondrial co-localization with the actin cytoskeletal. Examination of proteins associated with mitochondria fission and apoptosis and mitophagy, respectively, showed levels of immunological expression consistent with the increasing fission and apoptosis in mineralizing cultures. These results suggest that chondrocyte differentiation is associated with intracellular structural reorganization, promoting increased mitochondria content and fusion that enables increased oxidative metabolism. Mineralization, however, does not need energy derived from oxidative metabolism; rather, during mineralization, mitochondria appear to undergo fission and mitophagy. In summary, these studies show that as chondrocytes underwent hypertrophic differentiation, they increased oxidative metabolism, but as mineralization proceeds, metabolism decreased. Mitochondria structure also underwent a structural reorganization that was further supportive of their oxidative capacity as the chondrocytes progressed through their differentiation. Thus, the mitochondria first underwent fusion to support increased oxidative metabolism, then underwent fission during mineralization, facilitating their programed death.

Oxidative metabolism is impaired by phosphate deficiency during fracture healing and is mechanistically related to BMP induced chondrocyte differentiation.

Prior studies of acute phosphate restriction during the endochondral phase of fracture healing showed delayed chondrocyte differentiation was mechanistically linked to decreased bone morphogenetic protein signaling. In the present study, transcriptomic analysis of fracture callus gene expression in three strains of mice was used to identify differentially expressed (FDR = q ≤ 0.05) genes in response to phosphate (Pi) restriction. Ontology and pathway analysis of these genes showed that independent of genetic background, a Pi-deficient diet downregulated (p = 3.16 × 10-23) genes associated with mitochondrial oxidative phosphorylation pathways as well as multiple other pathways of intermediate metabolism. Temporal clustering was used to identify co-regulation of these specific pathways. This analysis showed that specific Ox/Phos, tricarboxylic acid cycle, pyruvate dehydrogenase. Arginine, proline metabolism genes, and prolyl 4-hydroxylase were all co-regulated in response to dietary Pi restriction. The murine C3H10T½ mesenchymal stem cell line was used to assess the functional relationships between BMP2-induced chondrogenic differentiation, oxidative metabolism and extracellular matrix formation. BMP2-induced chondrogenic differentiation of C3H10T½ was carried out in culture media in the absence or presence of ascorbic acid, the necessary co-factor for proly hydroxylation, and in media with normal and 25 % phosphate levels. BMP2 treatment led to decreased proliferation, increased protein accumulation and increased collagen and aggrecan gene expression. Across all conditions, BMP2 increased total oxidative activity and ATP synthesis. Under all conditions, the presence of ascorbate further increased total protein accumulation, proly-hydroxylation and aggrecan gene expression, oxidative capacity and ATP production. Lower phosphate levels only diminished aggrecan gene expression with no other effects of metabolic activity being observed. These data suggest that dietary phosphate restriction controls endochondral growth in vivo indirectly through BMP signaling, which upregulates oxidative activity that is linked to overall protein production and collagen hydroxylation.

Prx1 cell subpopulations identified in various tissues with diverse quiescence and activation ability following fracture and BMP2 stimulation.

The expression of Prx1 has been used as a marker to define the skeletal stem cells (SSCs) populations found within the bone marrow and periosteum that contribute to bone regeneration. However, Prx1 expressing SSCs (Prx1-SSCs) are not restricted to the bone compartments, but are also located within the muscle and able to contribute to ectopic bone formation. Little is known however, about the mechanism(s) regulating Prx1-SSCs that reside in muscle and how they participate in bone regeneration. This study compared both the intrinsic and extrinsic factors of the periosteum and muscle derived Prx1-SSCs and analyzed their regulatory mechanisms of activation, proliferation, and skeletal differentiation. There was considerable transcriptomic heterogeneity in the Prx1-SSCs found in muscle or the periosteum however in vitro cells from both tissues showed tri-lineage (adipose, cartilage and bone) differentiation. At homeostasis, periosteal-derived Prx1 cells were proliferative and low levels of BMP2 were able to promote their differentiation, while the muscle-derived Prx1 cells were quiescent and refractory to comparable levels of BMP2 that promoted periosteal cell differentiation. The transplantation of Prx1-SCC from muscle and periosteum into either the same site from which they were isolated, or their reciprocal sites showed that periosteal cell transplanted onto the surface of bone tissues differentiated into bone and cartilage cells but was incapable of similar differentiation when transplanted into muscle. Prx1-SSCs from the muscle showed no ability to differentiate at either site of transplantation. Both fracture and ten times the BMP2 dose was needed to promote muscle-derived cells to rapidly enter the cell cycle as well as undergo skeletal cell differentiation. This study elucidates the diversity of the Prx1-SSC population showing that cells within different tissue sites are intrinsically different. While muscle tissue must have factors that promote Prx1-SSC to remain quiescent, either bone injury or high levels of BMP2 can activate these cells to both proliferate and undergo skeletal cell differentiation. Finally, these studies raise the possibility that muscle SSCs are potential target for skeletal repair and bone diseases.

Transcriptome-wide association study and eQTL colocalization identify potentially causal genes responsible for human bone mineral density GWAS associations.

Genome-wide association studies (GWASs) for bone mineral density (BMD) in humans have identified over 1100 associations to date. However, identifying causal genes implicated by such studies has been challenging. Recent advances in the development of transcriptome reference datasets and computational approaches such as transcriptome-wide association studies (TWASs) and expression quantitative trait loci (eQTL) colocalization have proven to be informative in identifying putatively causal genes underlying GWAS associations. Here, we used TWAS/eQTL colocalization in conjunction with transcriptomic data from the Genotype-Tissue Expression (GTEx) project to identify potentially causal genes for the largest BMD GWAS performed to date. Using this approach, we identified 512 genes as significant using both TWAS and eQTL colocalization. This set of genes was enriched for regulators of BMD and members of bone relevant biological processes. To investigate the significance of our findings, we selected PPP6R3, the gene with the strongest support from our analysis which was not previously implicated in the regulation of BMD, for further investigation. We observed that Ppp6r3 deletion in mice decreased BMD. In this work, we provide an updated resource of putatively causal BMD genes and demonstrate that PPP6R3 is a putatively causal BMD GWAS gene. These data increase our understanding of the genetics of BMD and provide further evidence for the utility of combined TWAS/colocalization approaches in untangling the genetics of complex traits.

Inhibition of Osteoblast Differentiation by JAK2V617F Megakaryocytes Derived From Male Mice With Primary Myelofibrosis.

Past studies described interactions between normal megakaryocytes, the platelet precursors, and bone cell precursors in the bone marrow. This relationship has also been studied in context of various mutations associated with increased number of megakaryocytes. The current study is the first to examine the effects of megakaryocytes from transgenic mice carrying the most common mutation that causes primary myelofibrosis (PMF) in humans (JAK2V617F) on bone cell differentiation. Organ level assessments of mice using micro-computed tomography showed decreased bone volume in JAK2V617F males, compared to matching controls. Tissue level histology revealed increased deposition of osteoid (bone matrix prior mineralization) in these mutated mice, suggesting an effect on osteoblast differentiation. Mechanistic studies using a megakaryocyte-osteoblast co-culture system, showed that both wild type or JAK2V617F megakaryocytes derived from male mice inhibited osteoblast differentiation, but JAK2V617F cells exerted a more significant inhibitory effect. A mouse mRNA osteogenesis array showed increased expression of Noggin, Chordin, Alpha-2-HS-glycoprotein, Collagen type IV alpha 1 and Collagen type XIV alpha 1 (mostly known to inhibit bone differentiation), and decreased expression of alkaline phosphatase, Vascular cell adhesion molecule 1, Sclerostin, Distal-less homeobox 5 and Collagen type III alpha 1 (associated with osteogenesis) in JAK2V617F megakaryocytes, compared to controls. This suggested that the mutation re-programs megakaryocytes to express a cluster of genes, which together could orchestrate greater suppression of osteogenesis in male mice. These findings provide mechanistic insight into the effect of JAK2V617F mutation on bone, encouraging future examination of patients with this or other PMF-inducing mutations.

Identification of Known and Novel Long Noncoding RNAs Potentially Responsible for the Effects of Bone Mineral Density (BMD) Genomewide Association Study (GWAS) Loci.

Osteoporosis, characterized by low bone mineral density (BMD), is the most common complex disease affecting bone and constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 associations influencing BMD. It has been shown that perturbations to long noncoding RNAs (lncRNAs) influence BMD and the activities of bone cells; however, the extent to which lncRNAs are involved in the genetic regulation of BMD is unknown. Here, we combined the analysis of allelic imbalance (AI) in human acetabular bone fragments with a transcriptome-wide association study (TWAS) and expression quantitative trait loci (eQTL) colocalization analysis using data from the Genotype-Tissue Expression (GTEx) project to identify lncRNAs potentially responsible for GWAS associations. We identified 27 lncRNAs in bone that are located in proximity to a BMD GWAS association and harbor single-nucleotide polymorphisms (SNPs) demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (false discovery rate [FDR] correction < 0.05) with BMD through TWAS and had a colocalizing eQTL (regional colocalization probability [RCP] > 0.1). The 58 lncRNAs are located in 43 BMD associations. To further support a causal role for the identified lncRNAs, we show that 23 of the 58 lncRNAs are differentially expressed as a function of osteoblast differentiation. Our approach identifies lncRNAs that are potentially responsible for BMD GWAS associations and suggest that lncRNAs play a role in the genetics of osteoporosis. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Temporal and Quantitative Transcriptomic Differences Define Sexual Dimorphism in Murine Postnatal Bone Aging.

Time is a central element of the sexual dimorphic patterns of development, pathology, and aging of the skeleton. Because the transcriptome is a representation of the phenome, we hypothesized that both sex and sex-specific temporal, transcriptomic differences in bone tissues over an 18-month period would be informative to the underlying molecular processes that lead to postnatal sexual dimorphism. Regardless of age, sex-associated changes of the whole bone transcriptomes were primarily associated not only with bone but also vascular and connective tissue ontologies. A pattern-based approach used to screen the entire Gene Expression Omnibus (GEO) database against those that were sex-specific in bone identified two coordinately regulated gene sets: one related to high phosphate-induced aortic calcification and one induced by mechanical stimulation in bone. Temporal clustering of the transcriptome identified two skeletal tissue-associated, sex-specific patterns of gene expression. One set of genes, associated with skeletal patterning and morphology, showed peak expression earlier in females. The second set of genes, associated with coupled remodeling, had quantitatively higher expression in females and exhibited a broad peak between 3 to 12 months, concurrent with the animals' reproductive period. Results of phenome-level structural assessments of the tibia and vertebrae, and in vivo and in vitro analysis of cells having osteogenic potential, were consistent with the existence of functionally unique, skeletogenic cell populations that are separately responsible for appositional growth and intramedullary functions. These data suggest that skeletal sexual dimorphism arises through sex-specific, temporally different processes controlling morphometric growth and later coupled remodeling of the skeleton during the reproductive period of the animal. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Post natal expression of Prx1 labels appendicular restricted progenitor cell populations of multiple tissues.

Currently, there is no consensus whether there is a single or multiple postnatal stem cell population(s) that contribute to skeletal homeostasis and postnatal bone formation. A known population of cells that express Prx1 contributes to postnatal bone formation. Prx1 expression also connotes calvaria and appendicular tissues during embryonic development. A transgenic tamoxifen inducible Prx1 reporter mouse was used for lineage tracking, to characterize the postnatal contribution of Prx1 expressing cells in skeletal homeostasis and bone formation. Under homeostatic conditions Prx1 labeling gave rise to a transient yet rapid turnover cell population at the periosteal and endosteal surfaces, along muscle fibers, and within the medial layers of vessels both within the muscle and marrow compartments of the appendicular skeleton. Fracture and ectopic bone formation of both fore and hind limbs showed recruitment and expansion of Prx1-derived cells in newly formed bone tissues. Prx1 labeled cells were limited or absent at axial skeletal sites during both homeostasis and after induction of bone formation. Last, Prx1-derived cells differentiated into multiple cell lineages including vascular smooth muscle, adipose, cartilage, and bone cells. These results show that Prx1 expression retained its embryonic tissue specification and connotes a stem/progenitor cell populations of mesenchymal tissue progenitors.

Spatial assessment of femoral neck bone density and microstructure in hip osteoarthritis.

Osteoarthritis (OA) is known to involve profound changes in bone density and microstructure near to, and even distal to, the joint. Critically, however, a full, spatial picture of these abnormalities has not been well documented in a quantitative fashion in hip OA. Here, micro-computed tomography (44.8 µm/voxel) and data-driven computational anatomy were used to generate 3-D maps of the distribution of bone density and microstructure in human femoral neck samples with early (6F/4M, mean age = 51.3 years), moderate (14F/8M, mean age = 60 years), and severe (16F/6M, mean age = 63.3 years) radiographic OA. With increasing severity of radiographic OA, there was decreased cortical bone mineral density (BMD) (p=0.003), increased cortical thickness (p=0.001), increased cortical porosity (p=0.0028), and increased cortical cross-sectional area (p=0.0012, due to an increase in periosteal radius (p=0.018)), with no differences detected in the total femoral neck or trabecular compartment measures. No OA-related region-specific differences were detected through Statistical Parametric Mapping, but there were trends towards decreased tissue mineral density (TMD) in the inferior femoral neck with increasing OA severity (0.050 < p ≤ 0.091), possibly due to osteophytes. Overall, the lack of differences in cortical TMD among radiographic OA groups indicated that the decrease in cortical BMD with increasing OA severity was largely due to the increased cortical porosity rather than decreased tissue mineralization. As porosity is inversely associated with stiffness and strength in cortical bone, increased porosity may offset the effect that increased cortical cross-sectional area would be expected to have on reducing stresses within the femoral neck. The use of high-resolution imaging and quantitative spatial assessment in this study provide insight into the heterogeneous and multi-faceted changes in density and microstructure in hip OA, which have implications for OA progression and fracture risk.

Maternal GNAS Contributes to the Extra-Large G Protein α-Subunit (XLαs) Expression in a Cell Type-Specific Manner.

GNAS encodes the stimulatory G protein alpha-subunit (Gsα) and its large variant XLαs. Studies have suggested that XLαs is expressed exclusively paternally. Thus, XLαs deficiency is considered to be responsible for certain findings in patients with paternal GNAS mutations, such as pseudo-pseudohypoparathyroidism, and the phenotypes associated with maternal uniparental disomy of chromosome 20, which comprises GNAS. However, a study of bone marrow stromal cells (BMSC) suggested that XLαs could be biallelically expressed. Aberrant BMSC differentiation due to constitutively activating GNAS mutations affecting both Gsα and XLαs is the underlying pathology in fibrous dysplasia of bone. To investigate allelic XLαs expression, we employed next-generation sequencing and a polymorphism common to XLαs and Gsα, as well as A/B, another paternally expressed GNAS transcript. In mouse BMSCs, Gsα transcripts were 48.4 ± 0.3% paternal, while A/B was 99.8 ± 0.2% paternal. In contrast, XLαs expression varied among different samples, paternal contribution ranging from 43.0 to 99.9%. Sample-to-sample variation in paternal XLαs expression was also detected in bone (83.7-99.6%) and cerebellum (83.8 to 100%) but not in cultured calvarial osteoblasts (99.1 ± 0.1%). Osteoblastic differentiation of BMSCs shifted the paternal XLαs expression from 83.9 ± 1.5% at baseline to 97.2 ± 1.1%. In two human BMSC samples grown under osteoinductive conditions, XLαs expression was also predominantly monoallelic (91.3 or 99.6%). Thus, the maternal GNAS contributes significantly to XLαs expression in BMSCs but not osteoblasts. Altered XLαs activity may thus occur in certain cell types irrespective of the parental origin of a GNAS defect.

Bayesian modeling of factorial time-course data with applications to a bone aging gene expression study.

Many scientific studies, especially in the biomedical sciences, generate data measured simultaneously over a multitude of units, over a period of time, and under different conditions or combinations of factors. Often, an important question of interest asked relates to which units behave similarly under different conditions, but measuring the variation over time complicates the analysis significantly. In this article we address such a problem arising from a gene expression study relating to bone aging, and develop a Bayesian statistical method that can simultaneously detect and uncover signals on three levels within such data: factorial, longitudinal, and transcriptional. Our model framework considers both cluster and time-point-specific parameters and these parameters uniquely determine the shapes of the temporal gene expression profiles, allowing the discovery and characterization of latent gene clusters based on similar underlying biological mechanisms. Our methodology was successfully applied to discover transcriptional networks in a microarray data set comparing the transcriptomic changes that occurred during bone aging in male and female mice expressing one or both copies of the bromodomain (Brd2) gene, a transcriptional regulator which exhibits an age-dependent sex-linked bone loss phenotype.

Global transcriptomic analysis of a murine osteocytic cell line subjected to spaceflight.

Bone loss is a major health concern for astronauts during long-term spaceflight and for patients during prolonged bed rest or paralysis. Growing evidence suggests that osteocytes, the most abundant cells in the mineralized bone matrix, play a key role in sensing mechanical forces applied to the skeleton and integrating the orchestrated response into subcellular biochemical signals to modulate bone homeostasis. However, the precise molecular mechanisms underlying both mechanosensation and mechanotransduction in late-osteoblast-to-osteocyte cells under microgravity (µG) have yet to be elucidated. To unravel the mechanisms by which late osteoblasts and osteocytes sense and respond to mechanical unloading, we exposed the osteocytic cell line, Ocy454, to 2, 4, or 6 days of µG on the SpaceX Dragon-6 resupply mission to the International Space Station. Our results showed that µG impairs the differentiation of osteocytes, consistent with prior osteoblast spaceflight experiments, which resulted in the downregulation of key osteocytic genes. Importantly, we demonstrate the modulation of critical glycolysis pathways in osteocytes subjected to microgravity and discovered a set of mechanical sensitive genes that are consistently regulated in multiple cell types exposed to microgravity suggesting a common, yet to be fully elucidated, genome-wide response to microgravity. Ground-based simulated microgravity experiments utilizing the NASA rotating-wall-vessel were unable to adequately replicate the changes in microgravity exposure highlighting the importance of spaceflight missions to understand the unique environmental stress that microgravity presents to diverse cell types. In summary, our findings demonstrate that osteocytes respond to µG with an increase in glucose metabolism and oxygen consumption.

Overview of Skeletal Repair (Fracture Healing and Its Assessment).

The study of postnatal skeletal repair is of immense clinical interest. Optimal repair of skeletal tissue is necessary in all varieties of elective and reparative orthopedic surgical treatments. However, the repair of fractures is unique in this context in that fractures are one of the most common traumas that humans experience and are the end-point manifestation of osteoporosis, the most common chronic disease of aging. In the first part of this introduction the basic biology of fracture healing is presented. The second part discusses the primary methodological approaches that are used to examine repair of skeletal hard tissue and specific considerations for choosing among and implementing these approaches.

Generation of Closed Transverse Fractures in Small Animals.

The most common procedure that has been developed for use in rats and mice to model fracture healing is described. The nature of the regenerative processes that may be assessed and the types of research questions that may be addressed with this model are briefly outlined. The detailed surgical protocol to generate closed simple transverse fractures is presented and general considerations when setting up an experiment using this model are described.

LOXL2 promotes aggrecan and gender-specific anabolic differences to TMJ cartilage.

In the United States, 5-12% of adults have at least one symptom of temporomandibular joint (TMJ) disorders, including TMJ osteoarthritis (TMJ-OA). However, there is no chondroprotective agent that is approved for clinical application. We showed that LOXL2 is elevated in the regenerative response during fracture healing in mice and has a critical role in chondrogenic differentiation. Indeed, LOXL2 is an anabolic effector that attenuates pro-inflammatory signaling in OA cartilage of the TMJ and knee joint, induces chondroprotective and regenerative responses, and attenuates NF-kB signaling. The specific goal of the study was to evaluate if adenoviral delivery of LOXL2 is anabolic to human and mouse TMJ condylar cartilage in vivo and evaluate the protective and anabolic effect on cartilage-specific factors. We employed two different models to assess TMJ-OA. In one model, clinical TMJ-OA cartilage from 5 different samples in TMJ-OA cartilage plugs were implanted subcutaneously in nude mice. Adenovirus LOXL2 -treated implants showed higher mRNA levels of LOXL2, ACAN, and other anabolic genes compared to the adenovirus-Empty-treated implants. Further characterization by RNA-seq analysis showed LOXL2 promotes proteoglycan networks and extracellular matrix in human TMJ-OA cartilage implants in vivo. In order to evaluate if LOXL2-induced functional and sex-linked differences, both male and female four-month-old chondrodysplasia (Cho/+) mice, which develop progressive TMJ-OA due to a point mutation in the Col11a1 gene, were subjected to intraperitoneal injection with Adv-RFP-LOXL2 every 2 weeks for 12 weeks. The data showed that adenovirus delivery of LOXL2 upregulated LOXL2 and aggrecan (Acan), whereas MMP13 expression was slightly downregulated. The fold change expression of Acan and Runx2 induced by Adv-RFP-LOXL2 was higher in females compared to males. Interestingly, Adv-RFP-LOXL2 injection significantly increased Rankl expression in male but there was no change in females, whereas VegfB gene expression was increased in females, but not in males, as compared to those injected with Adv-RFP-Empty in respective groups. Our findings indicate that LOXL2 can induce specifically the expression of Acan and other anabolic genes in two preclinical models in vivo. Further, LOXL2 has beneficial functions in human TMJ-OA cartilage implants and promotes gender-specific anabolic responses in Cho/+ mice with progressive TMJ-OA, suggesting its merit for further study as an anabolic therapy for TMJ-OA.

An in vivo Comparison Study Between Strontium Nanoparticles and rhBMP2.

The osteoinductive property of strontium was repeatedly proven in the last decades. Compelling in vitro data demonstrated that strontium hydroxyapatite nanoparticles exert a dual action, by promoting osteoblasts-driven matrix secretion and inhibiting osteoclasts-driven matrix resorption. Recombinant human bone morphogenetic protein 2 (rhBMP2) is a powerful osteoinductive biologic, used for the treatment of vertebral fractures and critically-sized bone defects. Although effective, the use of rhBMP2 has limitations due its recombinant morphogen nature. In this study, we examined the comparison between two osteoinductive agents: rhBMP2 and the innovative strontium-substituted hydroxyapatite nanoparticles. To test their effectiveness, we independently loaded Gelfoam sponges with the two osteoinductive agents and used the sponges as agent-carriers. Gelfoam are FDA-approved biodegradable medical devices used as delivery system for musculoskeletal defects. Their porous structure and spongy morphology make them attractive in orthopedic field. The abiotic characterization of the loaded sponges, involving ion release pattern and structure investigation, was followed by in vivo implantation onto the periosteum of healthy mice and comparison of the effects induced by each implant was performed. Abiotic analysis demonstrated that strontium was continuously released from the sponges over 28 days with a pattern similar to rhBMP2. Histological observations and gene expression analysis showed stronger endochondral ossification elicited by strontium compared to rhBMP2. Osteoclast activity was more inhibited by strontium than by rhBMP2. These results demonstrated the use of sponges loaded with strontium nanoparticles as potential bone grafts might provide better outcomes for complex fractures. Strontium nanoparticles are a novel and effective non-biologic treatment for bone injuries and can be used as novel powerful therapeutics for bone regeneration.

Epidemiology, Clinical Assessments, and Current Treatments of Nonunions.

PURPOSE OF REVIEW: The failure of bony union following a fracture, termed a fracture nonunion, has severe patient morbidity and economic consequences. This review describes current consensuses and future directions of investigation for determining why, detecting when, and effective treatment if this complication occurs. RECENT FINDINGS: Current nonunion investigation is emphasizing an expanded understanding of the biology of healing. This has led to assessments of the immune environment, multiple cytokines and morphogenetic factors, and the role of skeletogenic stem cells in the development of nonunion. Detecting biological markers and other objective diagnostic criteria is also a current objective of nonunion research. Treatment approaches in the near future will likely be dominated by the development of specific adjunct therapies to the nonunion surgical management, which will be informed by an expanded mechanistic understanding of nonunion biology. Current consensus among orthopedists is that improved diagnosis and treatment of nonunion hinges first on discoveries at the bench side with later translation to the clinic.

Lysyl Oxidase-Like 2 Protects against Progressive and Aging Related Knee Joint Osteoarthritis in Mice.

BACKGROUND: The goal of this study was to determine if adenovirus-delivered LOXL2 protects against progressive knee osteoarthritis (OA), assess its specific mechanism of action; and determine if the overexpression of LOXL2 in transgenic mice can protect against the development of OA-related cartilage damage and joint disability. METHODS: Four-month-old Cho/+ male and female mice were intraperitoneally injected with either Adv-RFP-LOXL2 or an empty vector twice a month for four months. The proteoglycan levels and the expression of anabolic and catabolic genes were examined by immunostaining and qRT-PCR. The effect of LOXL2 expression on signaling was tested via the pro-inflammatory cytokine IL1ß in the cartilage cell line ATDC5. Finally; the OA by monosodium iodoacetate (MIA) injection was also induced in transgenic mice with systemic overexpression of LOXL2 and examined gene expression and joint function by treadmill tests and assessment of allodynia. RESULTS: The adenovirus treatment upregulated LOXL2; Sox9; Acan and Runx2 expression in both males and females. The Adv-RFP-LOXL2 injection; but not the empty vector injection increased proteoglycan staining and aggrecan expression but reduced MMP13 expression. LOXL2 attenuated IL-1ß-induced phospho-NF-κB/p65 and rescued chondrogenic lineage-related genes in ATDC5 cells; demonstrating one potential protective mechanism. LOXL2 attenuated phospho-NF-κB independent of its enzymatic activity. Finally; LOXL2-overexpressing transgenic mice were protected from MIA-induced OA-related functional changes; including the time and distance traveled on the treadmill and allodynia. CONCLUSION: Our study demonstrates that systemic LOXL2 adenovirus or LOXL2 genetic overexpression in mice can protect against OA. These findings demonstrate the potential for LOXL2 gene therapy for knee-OA clinical treatment in the future.

Clustering of temporal gene expression data with mixtures of mixed effects models with a penalized likelihood.

MOTIVATION: Clustering algorithms like K-Means and standard Gaussian mixture models (GMM) fail to account for the structure of variability of replicated data or repeated measures over time. Additionally, a priori cluster number assumptions add an additional complexity to the process. Current methods to optimize cluster labels and number can be inaccurate or computationally intensive for temporal gene expression data with this additional variability. RESULTS: An extension to a model-based clustering algorithm is proposed using mixtures of mixed effects polynomial regression models and the EM algorithm with an entropy penalized log-likelihood function (EPEM). The EPEM is used to cluster temporal gene expression data with this additional variability. The addition of random effects in our model decreased the misclassification error when compared to mixtures of fixed effects models or other methods such as K-Means and GMM. Applying our method to microarray data from a fracture healing study revealed distinct temporal patterns of gene expression. AVAILABILITY AND IMPLEMENTATION: https://github.com/darlenelu72/EPEM-GMM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Hypophosphatemia Regulates Molecular Mechanisms of Circadian Rhythm.

Transcriptomic analysis showed that the central circadian pathway genes had significantly altered expression in fracture calluses from mice fed a low phosphate diet. This led us to hypothesize that phosphate deficiency altered the circadian cycle in peripheral tissues. Analysis of the expression of the central clock genes over a 24-36 hour period in multiple peripheral tissues including fracture callus, proximal tibia growth plate and cardiac tissues after 12 days on a low phosphate diet showed higher levels of gene expression in the hypophosphatemia groups (p < 0.001) and a 3 to 6 hour elongation of the circadian cycle. A comparative analysis of the callus tissue transcriptome genes that were differentially regulated by hypophosphatemia with published data for the genes in bone that are diurnally regulated identified 1879 genes with overlapping differential regulation, which were shown by ontology assessment to be associated with oxidative metabolism and apoptosis. Network analysis of the central circadian pathway genes linked their expression to the up regulated expression of the histone methyltransferase gene EZH2, a gene that when mutated in both humans and mice controls overall skeletal growth. These data suggest that phosphate is an essential metabolite that controls circadian function in both skeletal and non skeletal peripheral tissues and associates its levels with the overall oxidative metabolism and skeletal growth of animals.