







Abstract Use of high-resolution micro–computed tomography (µCT) imaging to assess trabecular and cortical bone morphology has grown immensely. There are several commercially available µCT systems, each with different approaches to image acquisition, evaluation, and reporting of outcomes. This lack of consistency makes it difficult to interpret reported results and to compare findings across different studies. This article addresses this critical need for standardized terminology and consistent reporting of parameters related to image acquisition and analysis, and key outcome assessments, particularly with respect to ex vivo analysis of rodent specimens. Thus the guidelines herein provide recommendations regarding (1) standardized terminology and units, (2) information to be included in describing the methods for a given experiment, and (3) a minimal set of outcome variables that should be reported. Whereas the specific research objective will determine the experimental design, these guidelines are intended to ensure accurate and consistent reporting of µCT-derived bone morphometry and density measurements. In particular, the methods section for papers that present µCT-based outcomes must include details of the following scan aspects: (1) image acquisition, including the scanning medium, X-ray tube potential, and voxel size, as well as clear descriptions of the size and location of the volume of interest and the method used to delineate trabecular and cortical bone regions, and (2) image processing, including the algorithms used for image filtration and the approach used for image segmentation. Morphometric analyses should be based on 3D algorithms that do not rely on assumptions about the underlying structure whenever possible. When reporting µCT results, the minimal set of variables that should be used to describe trabecular bone morphometry includes bone volume fraction and trabecular number, thickness, and separation. The minimal set of variables that should be used to describe cortical bone morphometry includes total cross-sectional area, cortical bone area, cortical bone area fraction, and cortical thickness. Other variables also may be appropriate depending on the research question and technical quality of the scan. Standard nomenclature, outlined in this article, should be followed for reporting of results. © 2010 American Society for Bone and Mineral Research
Integrated Bone Tissue Physiology: Anatomy and Physiology
Standardized nomenclature, symbols, and units for bone histomorphometry: A 2012 update of the report of the ASBMR Histomorphometry Nomenclature Committee
Before publication of the original version of this report in 1987, practitioners of bone histomorphometry communicated with each other in a variety of arca

Ten‐Year Simulation of the Effects of Denosumab on Bone Remodeling in Human Biopsies
ABSTRACT Postmenopausal osteoporosis is a disease manifesting in degradation of bone mass and microarchitecture, leading to weakening and increased risk of fracture. Clinical trials are an essential tool for evaluating new treatments and may provide further mechanistic understanding of their effects in vivo. However, the histomorphometry from clinical trials is limited to 2D images and reflects single time points. Biochemical markers of bone turnover give global insight into a drug's action, but not the local dynamics of the bone remodeling process and the cells involved. Additionally, comparative trials necessitate separate treatment groups, meaning only aggregated measures can be compared. In this study, in silico modeling based on histomorphometry and pharmacokinetic data was used to assess the effects of treatment versus control on μCT scans of the same biopsy samples over time, matching the changes in bone volume fraction observed in biopsies from denosumab and placebo groups through year 10 of the FREEDOM Extension trial. In the simulation, treatment decreased osteoclast number, which led to a modest increase in trabecular thickness and osteocyte stress shielding. Long‐term bone turnover suppression led to increased RANKL production, followed by a small increase in osteoclast number at the end of the 6‐month–dosing interval, especially at the end of the Extension study. Lack of treatment led to a significant loss of bone mass and structure. The study's results show how in silico models can generate predictions of denosumab cellular action over a 10‐year period, matching static and dynamic morphometric measures assessed in clinical biopsies. The use of in silico models with clinical trial data can be a method to gain further insight into fundamental bone biology and how treatments can perturb this. With rigorous validation, such models could be used for informing the design of clinical trials, such that the number of participants could be reduced to a minimum to show efficacy. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.
A multiscale mechanobiological model of bone remodelling predicts site-specific bone loss in the femur during osteoporosis and mechanical disuse
We propose a multiscale mechanobiological model of bone remodelling to investigate the site-specific evolution of bone volume fraction across the midshaft of a femur. The model includes hormonal regulation and biochemical coupling of bone cell populations, the influence of the microstructure on bone turnover rate, and mechanical adaptation of the tissue. Both microscopic and tissue-scale stress/strain states of the tissue are calculated from macroscopic loads by a combination of beam theory and micromechanical homogenisation. This model is applied to simulate the spatio-temporal evolution of a human midshaft femur scan subjected to two deregulating circumstances: (i) osteoporosis and (ii) mechanical disuse. Both simulated deregulations led to endocortical bone loss, cortical wall thinning and expansion of the medullary cavity, in accordance with experimental findings. Our model suggests that these observations are attributable to a large extent to the influence of the microstructure on bone turnover rate. Mechanical adaptation is found to help preserve intracortical bone matrix near the periosteum. Moreover, it leads to non-uniform cortical wall thickness due to the asymmetry of macroscopic loads introduced by the bending moment. The effect of mechanical adaptation near the endosteum can be greatly affected by whether the mechanical stimulus includes stress concentration effects or not.

<i>In silico</i> biology of bone modelling and remodelling: adaptation
Modelling and remodelling are the processes by which bone adapts its shape and internal structure to external influences. However, the cellular mechanisms triggering osteoclastic resorption and osteoblastic formation are still unknown. In order to investigate current biological theories, in silico models can be applied. In the past, most of these models were based on the continuum assumption, but some questions related to bone adaptation can be addressed better by models incorporating the trabecular microstructure. In this paper, existing simulation models are reviewed and one of the microstructural models is extended to test the hypothesis that bone adaptation can be simulated without particular knowledge of the local strain distribution in the bone. Validation using an experimental murine loading model showed that this is possible. Furthermore, the experimental model revealed that bone formation cannot be attributed only to an increase in trabecular thickness but also to structural reorganization including the growth of new trabeculae. How these new trabeculae arise is still an unresolved issue and might be better addressed by incorporating other levels of hierarchy, especially the cellular level. The cellular level sheds light on the activity and interplay between the different cell types, leading to the effective change in the whole bone. For this reason, hierarchical multi-scale simulations might help in the future to better understand the biomathematical laws behind bone adaptation.

New perspectives on body size and shape evolution in dinosaurs
ABSTRACT Diversity in the body shapes and sizes of dinosaurs was foundational to their widespread success during the Mesozoic era. The ability to quantify body size and form reliably is therefore critical to the study of dinosaur biology and evolution. Body mass estimates for any given fossil animal are, in theory, most informative when derived from volumetric models that account for the three‐dimensional shapes of the entire body. In addition to providing estimates of total body mass, volumetric approaches can be used to determine the inertial properties of specific body segments and the overall distribution of mass throughout the body, each of which are essential for the modelling and interpretation of form–function relationships and their associations with ecology. However, the determination of body volumes in fossil taxa is often subjective, and may be sensitive to varied artistic inference. This highlights the need for an approach to body mass estimation in which body segment volumes are systematically constrained by quantitative scaling relationships between the hard tissues that fossilise and the soft tissues only observable in extant taxa. To this end, we used recently published skeletal to soft tissue volumetric scaling factors derived from CT data of extant sauropsids to estimate body segment mass properties from skeletal models of 52 non‐avian dinosaurs representing the majority of major clades and body plans. The body masses estimated by this study range from less than 200 g in the tiny avialan Yixianornis to over 60 tonnes in the giant sauropod Patagotitan , which is currently the largest dinosaur known from mostly complete skeletal remains. From our models, we infer that many previous reconstructions of soft tissue envelopes may be too small, and that many dinosaurs were therefore heavier than previous estimates. Our models generally overlap with the range of body mass estimates derived from limb bone shaft dimensions, but with considerable quantitative variability among major clades. This suggests that different taxa either differed in skeletal to soft tissue volume ratios, or that their limb bone dimensions varied relative to body mass, perhaps related to differences in locomotor dynamics and postural evolution. Our models also allowed us to investigate variation in mass distribution and body proportions across different dinosaurs from a perspective grounded in extant anatomical data, framing long‐standing hypotheses about their form, function, and behaviour in a quantitative context. For example, reconstructed disparity in whole‐body centres of mass reflects a broad array of postures in different dinosaur clades, while the lack of strong positive allometry in the dimensions of the weight‐bearing limb segments relative to total body mass corroborates previous studies suggesting an overall decrease in dinosaur locomotor performance as body size increased.

A new biological bone remodeling in silico model combined with advanced discretization methods
Abstract Bone remodeling remains a highly researched topic investigated by many strands of science. The main purpose of this work is formulating a new computational framework for biological simulation, extending the version of the bone remodeling model previously proposed by Komarova. Thus, considering only the biological aspect of the remodeling process, the action of osteoclasts and osteoblasts is taken into account as well as its impact on bone mass. It is conducted a spatiotemporal analysis of a remodeling cycle obtaining a dynamic behavior of bone cells very similar to the biological process already described in the literature. The numerical example used is based on bone images obtained with scanning electron microscopy. During simulation, it is possible to observe the variation of bone's architecture through isomaps. These maps are obtained through the combination of biological bone remodeling models with three distinct numerical techniques—finite element method (FEM), radial point interpolation method (RPIM), and natural neighbor radial point interpolation method (NNRPIM). A study combining these numerical techniques allows to compare their performance. Ultimately, this work supports the inclusion of meshless methods due to their smoother results and its easiness to be combined with medical images from CT scans and MRI.

Control of Bone Matrix Properties by Osteocytes
Osteocytes make up 90-95% of the cellular content of bone and form a rich dendritic network with a vastly greater surface area than either osteoblasts or osteoclasts. Osteocytes are well positioned to play a role in bone homeostasis by interacting directly with the matrix; however, the ability for these cells to modify bone matrix remains incompletely understood. With techniques for examining the nano- and microstructure of bone matrix components including hydroxyapatite and type I collagen becoming more widespread, there is great potential to uncover novel roles for the osteocyte in maintaining bone quality. In this review, we begin with an overview of osteocyte biology and the lacunar-canalicular system. Next, we describe recent findings from in vitro models of osteocytes, focusing on the transitions in cellular phenotype as they mature. Finally, we describe historical and current research on matrix alteration by osteocytes in vivo, focusing on the exciting potential for osteocytes to directly form, degrade, and modify the mineral and collagen in their surrounding matrix.

Modeling the Human Bone–Tumor Niche: Reducing and Replacing the Need for Animal Data
ABSTRACT Bone is the most common site for cancer metastasis. Understanding the interactions within the complex, heterogeneous bone–tumor microenvironment is essential for the development of new therapeutics. Various animal models of tumor‐induced bone disease are routinely used to provide valuable information on the relationship between cancer cells and the skeleton. However, new model systems exist that offer an alternative approach to the use of animals and might more accurately reveal the cellular interactions occurring within the human bone–tumor niche. This review highlights replacement models that mimic the bone microenvironment and where cancer metastases and tumor growth might be assessed alongside bone turnover. Such culture models include the use of calcified regions of animal tissue and scaffolds made from bone mineral hydroxyapatite, synthetic polymers that can be manipulated during manufacture to create structures resembling trabecular bone surfaces, gel composites that can be modified for stiffness and porosity to resemble conditions in the tumor–bone microenvironment. Possibly the most accurate model system involves the use of fresh human bone samples, which can be cultured ex vivo in the presence of human tumor cells and demonstrate similar cancer cell–bone cell interactions as described in vivo. In addition, the use of mathematical modeling and computational biology approaches provide an alternative to preliminary animal testing. The use of such models offers the capacity to mimic significant elements of the human bone–tumor environment, and complement, refine, or replace the use of preclinical models. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

<i>In silico</i> models of bone remodeling from macro to nano—from organ to cell
Abstract Computational modeling is a tool through which researchers can achieve a greater understanding of the mechanisms governing biological systems. In the field of bone biology, a plethora of models exist which attempt to replicate and investigate bone's dynamic behavior at different scales. At organ level, models are continuum based and describe the variation of bone's apparent density as a function of both biological and external mechanical stimuli. At tissue level, models include bone microarchitecture and more descriptive parameters such as trabecular thickness, osteoclast resorption depth, and activation frequency. Finally, at cell level, models employ partial differential equations to describe complex cellular interactions in the temporal domain. Although informative, these models exist in isolation. Consequently, their interpretation is limited. In this review, we present an overview of the organ‐, tissue‐, and cell‐level models and assess their ability to reflect bone's metabolic processes reliably. Existing interscale synergies are then presented along with a computational framework which could be exploited to achieve a fully integrated, multiscale modeling approach. WIREs Syst Biol Med 2011 3 241–251 DOI: 10.1002/wsbm.115 This article is categorized under: Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models

Wnt10b Increases Postnatal Bone Formation by Enhancing Osteoblast Differentiation
Abstract Overexpression of Wnt10b from the osteocalcin promoter in transgenic mice increases postnatal bone mass. Increases in osteoblast perimeter, mineralizing surface, and bone formation rate without detectable changes in pre-osteoblast proliferation, osteoblast apoptosis, or osteoclast number and activity suggest that, in this animal model, Wnt10b primarily increases bone mass by stimulating osteoblastogenesis. Introduction: Wnt signaling regulates many aspects of development including postnatal accrual of bone. Potential mechanisms for how Wnt signaling increases bone mass include regulation of osteoblast and/or osteoclast number and activity. To help differentiate between these possibilities, we studied mice in which Wnt10b is expressed specifically in osteoblast lineage cells or in mice devoid of Wnt10b. Materials and Methods: Transgenic mice, in which mouse Wnt10b is expressed from the human osteocalcin promoter (Oc-Wnt10b), were generated in C57BL/6 mice. Transgene expression was evaluated by RNase protection assay. Quantitative assessment of bone variables was done by radiography, μCT, and static and dynamic histomorphometry. Mechanisms of bone homeostasis were evaluated with assays for BrdU, TUNEL, and TRACP5b activity, as well as serum levels of C-terminal telopeptide of type I collagen (CTX). The endogenous role of Wnt10b in bone was assessed by dynamic histomorphometry in Wnt10b−/− mice. Results: Oc-Wnt10b mice have increased mandibular bone and impaired eruption of incisors during postnatal development. Analyses of femoral distal metaphyses show significantly higher BMD, bone volume fraction, and trabecular number. Increased bone formation is caused by increases in number of osteoblasts per bone surface, rate of mineral apposition, and percent mineralizing surface. Although number of osteoclasts per bone surface is not altered, Oc-Wnt10b mice have increased total osteoclast activity because of higher bone mass. In Wnt10b−/− mice, changes in mineralizing variables and osteoblast perimeter in femoral distal metaphyses were not observed; however, bone formation rate is reduced because of decreased total bone volume and trabecular number. Conclusions: High bone mass in Oc-Wnt10b mice is primarily caused by increased osteoblastogenesis, with a minor contribution from elevated mineralizing activity of osteoblasts.

Physiological Bone Remodeling: Systemic Regulation and Growth Factor Involvement
Bone remodeling is essential for adult bone homeostasis. It comprises two phases: bone formation and resorption. The balance between the two phases is crucial for sustaining bone mass and systemic mineral homeostasis. This review highlights recent work on physiological bone remodeling and discusses our knowledge of how systemic and growth factors regulate this process.

A review of mathematical modeling of bone remodeling from a systems biology perspective
Bone remodeling is an essential, delicately balanced physiological process of coordinated activity of bone cells that remove and deposit new bone tissue in the adult skeleton. Due to the complex nature of this process, many mathematical models of bone remodeling have been developed. Each of these models has unique features, but they have underlying patterns. In this review, the authors highlight the important aspects frequently found in mathematical models for bone remodeling and discuss how and why these aspects are included when considering the physiology of the bone basic multicellular unit, which is the term used for the collection of cells responsible for bone remodeling. The review also emphasizes the view of bone remodeling from a systems biology perspective. Understanding the systemic mechanisms involved in remodeling will help provide information on bone pathology associated with aging, endocrine disorders, cancers, and inflammatory conditions and enhance systems pharmacology. Furthermore, some features of the bone remodeling cycle and interactions with other organ systems that have not yet been modeled mathematically are discussed as promising future directions in the field.

A novel mathematical model of bone remodelling cycles for trabecular bone at the cellular level
After an initial phase of growth and development, bone undergoes a continuous cycle of repair, renewal and optimisation by a process called remodelling. This paper describes a novel mathematical model of the trabecular bone remodelling cycle. It is essentially formulated to simulate a remodelling event at a fixed position in the bone, integrating bone removal by osteoclasts and formation by osteoblasts. The model is developed to construct the variation in bone thickness at a particular point during the remodelling event, derived from standard bone histomorphometric analyses. The novelties of the approach are the adoption of a predator–prey model to describe the dynamic interaction between osteoclasts and osteoblasts, using a genetic algorithm–based solution; quantitative reconstruction of the bone remodelling cycle; and the introduction of a feedback mechanism in the bone formation activity to co-regulate bone thickness. The application of the model is first demonstrated by using experimental data recorded for normal (healthy) bone remodelling to predict the temporal variation in the number of osteoblasts and osteoclasts. The simulated histomorphometric data and remodelling cycle characteristics compare well with the specified input data. Sensitivity studies then reveal how variations in the model’s parameters affect its output; it is hoped that these parameters can be linked to specific biochemical factors in the future. Two sample pathological conditions, hypothyroidism and primary hyperparathyroidism, are examined to demonstrate how the model could be applied more broadly, and, for the first time, the osteoblast and osteoclast populations are predicted for these conditions. Further data are required to fully validate the model’s predictive capacity, but this work shows it has potential, especially in the modelling of pathological conditions and the optimisation of the treatment of those conditions.
Osteonal and hemi‐osteonal remodeling: The spatial and temporal framework for signal traffic in adult human bone
Abstract The bone replacement process in the adult skeleton is known as remodeling. When bone is removed by osteoclasts, new bone is laid down by osteoblasts in the same place, because the load bearing requirement is unchanged. Bone is usually replaced because it is too old to carry out its function, which is mainly mechanical in cortical bone and mainly support for homeostasis and hematopoiesis in cancellous bone. Remodeling always begins on a quiescent bone surface, separated from the marrow by flat lining cells that are one of the two modes of terminal differentiation of osteoblasts. Lining cells are gatekeepers, able to be informed of the need for remodeling, and to either execute or mediate all four components of its activation‐selection and preparation of the site, recruitment of mononuclear preosteoclasts, budding of new capillaries, and attraction of preosteoclasts to the chosen site where they fuse into multinucleated osteoclasts. In cortical bone, osteonal remodeling is carried out by a complex and unique structure, the basic multicellular unit (BMU) that comprises a cutting cone of osteoclasts in front, a closing cone lined by osteoblasts following behind, and connective tissue, blood vessels and nerves filling the cavity. The BMU maintains its size, shape and internal organization for many months as it travels through bone in a controlled direction. Individual osteoclast nuclei are short‐lived, turning over about 8% per d, replaced by new preosteoclasts that originated in the bone marrow and travel in the circulation to the site of resorption. Refilling of bone at each successive cross‐sectional location is accomplished by a team of osteoblasts, probably originating from precursors within the local connective tissue, all assembled within a narrow window of time, at the right location, and in the right orientation to the surface. Each osteoblast team forms bone most rapidly at its onset and slows down progressively. Some of the osteoblasts are buried as osteocytes, some die, and the remainder gradually assume the shape of lining cells. Cancellous bone is more accessible to study than cortical bone, but is geometrically complex. Although remodeling conforms to the same sequence of surface activation, resorption and formation, its three‐dimensional organization is difficult to visualize from two‐dimensional histologic sections. But the average sizes of resorption sites, formation sites, and completed structural units increase progressively, as they do in cortical bone, indicating that the cancellous BMU travels across the surface digging a trench rather than a tunnel, but maintaining its size, shape and individual identity by the continuous recruitment of new cells, just as in cortical bone, a process that can be visualized as hemiosteonal remodeling. The conclusion that all remodeling is carried out by individual BMUs has important implications for bone biology, since many questions about how BMUs operate cannot be answered by studying either intact organisms or isolated cell systems. Many different steps in remodeling and many factors that influence each step have been identified, but very little is known about how the process is regulated in vivo to achieve its biologic purposes; most factors studied to date are likely permissive rather than regulatory in nature. Based on the proposed conceptual model of the BMU, much in vitro experimentation is relevant to the growth, modeling and repair of bone, but not to its remodeling in the adult skeleton. Further progress in the understanding of in vivo physiology will require the characterization of gene expression in individual cells to be related to the spatial and temporal organization of the BMU. This is likely to be possible only for osteonal remodeling in cortical bone in which, because of its geometric simplicity, individual BMUs can consistently be observed in two‐dimensional, longitudinal sections. © 1994 Wiley‐Liss, Inc.

The gut microbiota regulates bone mass in mice
Abstract The gut microbiota modulates host metabolism and development of immune status. Here we show that the gut microbiota is also a major regulator of bone mass in mice. Germ-free (GF) mice exhibit increased bone mass associated with reduced number of osteoclasts per bone surface compared with conventionally raised (CONV-R) mice. Colonization of GF mice with a normal gut microbiota normalizes bone mass. Furthermore, GF mice have decreased frequency of CD4+ T cells and CD11b+/GR 1 osteoclast precursor cells in bone marrow, which could be normalized by colonization. GF mice exhibited reduced expression of inflammatory cytokines in bone and bone marrow compared with CONV-R mice. In summary, the gut microbiota regulates bone mass in mice, and we provide evidence for a mechanism involving altered immune status in bone and thereby affected osteoclast-mediated bone resorption. Further studies are required to evaluate the gut microbiota as a novel therapeutic target for osteoporosis. © 2012 American Society for Bone and Mineral Research.
