Covariation between forelimb muscle anatomy and bone shape in an Australian scratch-digging marsupial: Comparison of morphometric methods

The close association between muscle and bone is broadly intuitive; however, details of the covariation between the two has not been comprehensively studied. Without quantitative understanding of how muscle anatomy influences bone shape, it is difficult to draw conclusions of the significance of many morphological traits of the skeleton. In this study, we investigated these relationships in the Quenda (Isoodon fusciventer), a scratch-digging marsupial. We quantified the relationships between forelimb muscle anatomy and bone shape for animals representing a range of body masses (124–1,952 g) using two-block partial least square analyses. Muscle anatomy was quantified as muscle mass and physiological cross-sectional area (PCSA), and we used two morphometric methods to characterize bone shape: seven indices of linear bone proportions, and landmarks analysis. Bone shape was significantly correlated with body mass, reflecting allometric bone growth. Of the seven bone indices, only shoulder moment index (SMI) and ulna robustness index (URI) showed a significant covariation with muscle anatomy. Stronger relationships between muscle anatomy and forelimb bone shape were found using the landmark coordinates: muscle mass and PCSA were correlated with the geometric shape of the scapula, humerus, and third metacarpal, but to a lesser extent with shape of the ulna. Overall, our data show that landmark coordinates are more sensitive than bone indices to capturing shape changes evident throughout ontogeny, and is therefore a more appropriate method to investigate covariation with forelimb muscle anatomy. Single-species studies investigating ontogeny require refined methods to accurately develop understanding of the important relationships between muscle force generation and bone shape remodeling. Landmark analyses provide such a method.     

Three-dimensional topology optimization model to simulate the external shapes of bone

Elucidation of the mechanism by which the shape of bones is formed is essential for understanding vertebrate development. Bones support the body of vertebrates by withstanding external loads, such as those imposed by gravity and muscle tension. Many studies have reported that bone formation varies in response to external loads. An increased external load induces bone synthesis, whereas a decreased external load induces bone resorption. This relationship led to the hypothesis that bone shape adapts to external load. In fact, by simulating this relationship through topology optimization, the internal trabecular structure of bones can be successfully reproduced, thereby facilitating the study of bone diseases. In contrast, there have been few attempts to simulate the external structure of bones, which determines vertebrate morphology. However, the external shape of bones may be reproduced through topology optimization because cells of the same type form both the internal and external structures of bones. Here, we constructed a three-dimensional topology optimization model to attempt the reproduction of the external shape of teleost vertebrae. In teleosts, the internal structure of the vertebral bodies is invariable, exhibiting an hourglass shape, whereas the lateral structure supporting the internal structure differs among species. Based on the anatomical observations, we applied different external loads to the hourglass-shaped part. The simulations produced a variety of three-dimensional structures, some of which exhibited several structural features similar to those of actual teleost vertebrae. In addition, by adjusting the geometric parameters, such as the width of the hourglass shape, we reproduced the variation in the teleost vertebrae shapes. These results suggest that a simulation using topology optimization can successfully reproduce the external shapes of teleost vertebrae. By applying our topology optimization model to various bones of vertebrates, we can understand how the external shape of bones adapts to external loads.     

Relationship between fluctuating asymmetry,morphological modality and heterozygosity in an elderly Israeli population

We have been checking the following working hypotheses: 1) There is a negative correlation between genetic heterozygosity and fluctuating asymmetry (FA); 2) FA is a measure of developmental stability/instability of the whole organism, i.e. we expect negative correlation between FA and morphological proximity of a set of mass-size variables of an individual to a population centroid; and 3) FA is a measure of character-specific stability in a population, i.e. we expect correlation between magnitude of FA and deviation of an individual from the population centroid of the bilateral characters themselves. For this purpose each individual in a sample of about 200 elderly individuals was assessed for 11 polymorphic blood systems (14 genetic loci) as well as for a set of 26 anthropometric traits: 1) a set of ten mass size variables; and 2) a set of eight pairs of bilateral measurements. Four multivariate measures of morphological centrality were computed, two measures for size and two measures of shape distances from the i^th individual to the population centroid for mass-size variables and also for the bilateral variables. A multivariate measure of FA for 8 bilateral pairs was also computed. No relationship was detected between FA and heterozygosity, or between FA and any of the four multivariate deviations. Thus, we concluded that our data do not support the listed hypotheses.     

The behavior of adaptive bone-remodeling simulation models

The process of adaptive bone remodeling can be described mathematically and simulated in a computer model, integrated with the finite element method. In the model discussed here, cortical and trabecular bone are described as continuous materials with variable density. The remodeling rule applied to simulate the remodeling process in each element individually is, in fact, an objective function for an optimization process, relative to the external load. Its purpose is to obtain a constant, preset value for the strain energy per unit bone mass, by adapting the density. If an element in the structure cannot achieve that, it either turns to its maximal density (cortical bone) or resorbs completely.

It is found that the solution obtained in generally a discontinuous patchwork. For a two-dimensional proximal femur model this patchwork shows a good resemblance with the density distribution of a real proximal femur.

It is shown that the discontinuous end configuration is dictated by the nature of the differential equations describing the remodeling process. This process can be considered as a nonlinear dynamical system with many degrees of freedom, which behaves divergent relative to the objective, leading to many possible solutions. The precise solution is dependent on the parameters in the remodeling rule, the load and the initial conditions. The feedback mechanism in the process is self-enhancing; denser bone attracts more strain energy, whereby the bone becomes even more dense. It is suggested that this positive feedback of the attractor state (the strain energy field) creates order in the end configuration. In addition, the process ensures that the discontinuous end configuration is a structure with a relatively low mass, perhaps a minimal-mass structure, although this is no explicit objective in the optimization process.

It is hypothesized that trabecular bone is a chaotically ordered structure which can be considered as a fractal with characteristics of optimal mechanical resistance and minimal mass, of which the actual morphology depends on the local (internal) loading characteristics, the sensor-cell density and the degree of mineralization.     

Anisotropic properties of human cortical bone with osteogenesis imperfecta

The heterogeneity of bone shape and size variation is modulated by genetic, mechanical, nutritional, and hormonal patterning throughout its lifetime. Microstructural changes across cross sections are a result of mechanistic optimization that results over the years of evolution while being based on universal, time-invariant ingredients and patterns. Here we report changes across anatomical sections of bone with osteogenesis imperfecta (OI) that undermines the work of evolution through genetic mutation. This work examines the microstructure and molecular composition of different anatomical positions (anterior, medial, posterior, and lateral regions) in the diaphysis of an OI human tibia. The study shows that although there is no significant microstructural difference, molecular changes are observed using FTIR revealing differences in molecular composition of the four anatomical positions. In addition, the nanomechanical properties of anterior section of OI bone seem more heterogeneous. The nanomechanical properties of interstitial lamellae in all these bone samples are consistently greater than those of osteonal lamellae. The nanomechanical properties of bone depend on its anatomical section and on the measurement direction as well. Variations in molecular structure with anatomical positions and also corresponding differences in nanomechanical properties are reported. These are compared to those observed typically in healthy bone illustrating the unique influence of OI on bone multiscale behavior which results from an evolutionary process lasting for many years.     

Comparison of different hip prosthesis shapes considering micro-level bone remodeling and stress-shielding criteria using three-dimensional design space topology optimization

Since the late 1980s, computational analysis of total hip arthroplasty (THA) prosthesis components has been completed using macro-level bone remodeling algorithms. The utilization of macro-sized elements requires apparent bone densities to predict cancellous bone strength, thereby, preventing visualization and analysis of realistic trabecular architecture. In this study, we utilized a recently developed structural optimization algorithm, design space optimization (DSO), to perform a micro-level three-dimensional finite element bone remodeling simulation on the human proximal femur pre- and post-THA. The computational simulation facilitated direct performance comparison between two commercially available prosthetic implant stems from Zimmer Inc.: the Alloclassic and the Mayo conservative. The novel micro-level approach allowed the unique ability to visualize the trabecular bone adaption post-operation and to quantify the changes in bone mineral content by region. Stress-shielding and strain energy distribution were also quantified for the immediate post-operation and the stably fixated, post-remodeling conditions. Stress-shielding was highest in the proximal region and remained unchanged post-remodeling; conversely, the mid and distal portions show large increases in stress, suggesting a distal shift in the loadpath. The Mayo design conserves bone mass, while simultaneously reducing the incidence of stress-shielding compared to the Alloclassic, revealing a key benefit of the distinctive geometry. Several important factors for stable fixation, determined in clinical evaluations from the literature, were evident in both designs: high levels of proximal bone loss and distal bone densification. The results suggest this novel computational framework can be utilized for comparative hip prosthesis shape, uniquely considering the post-operation bone remodeling as a design criterion.     

Modelling external bone adaptation using evolutionary structural optimisation

External remodelling is significant in the bone healing process, and it is essential to predict the bone external shape in the design of artificial bone grafts. This paper demonstrates the effectiveness of the evolutionary structural optimisation (ESO) method for the simulation of bone morphology. A two-dimensional ESO strategy is developed which is capable of finding the modified bone topology beginning with any geometry under any loading conditions. The morphology of bone structure is described by the quantitative bone adaptation theory, which is integrated with the finite element method. The evolutionary topology optimisation process is introduced to find the bone shape. A rectangle, which occupies a larger space than the external shape of the bone structure, is specified as a design domain; the evolutionary process iteratively eliminates and redistributes material throughout the domain to obtain an optimum arrangement of bone materials. The technique has been tested on a wide range of examples. In this paper, the formation of trabecular bone architecture around an implant is studied; as another example, the growth of the coronal section of a vertebral body is predicted. The examples support the assertion that the external shape of bone structure can be successfully predicted by the proposed ESO procedure.     

Declining tibial curvature parallels ∼6150 years of decreasing mobility in central european agriculturalists

Long bones respond to mechanical loading through functional adaptation in a suite of morphological characteristics that together ensure structural competence to in vivo loading. As such, adult bone structure is often used to make inferences about past behavior from archaeological remains. However, such biomechanical approaches often investigate change in just one aspect of morphology, typically cross‐sectional morphology or trabecular structure. The relationship between longitudinal bone curvature and mobility patterns is less well understood, particularly in the tibia, and it is unknown how tibial curvature and diaphyseal cross‐sectional geometry interact to meet the structural requirements of loading. This study examines tibial curvature and its relationship with diaphyseal cross‐sectional geometry (CSG) and body size in preindustrial Central Europeans spanning ～6150 years following the introduction of agriculture in the region. Anteroposterior centroid displacement from the proximo‐distal longitudinal axis was quantified at nine diaphyseal section locations (collectively representative of diaphyseal curvature) in 216 tibial three‐dimensional laser scans. Results documented significant and corresponding temporal declines in midshaft centroid displacement and CSG properties. Significant correlations were found between mid‐diaphyseal centroid displacement and all mobility‐related CSG properties, while the relationship weakened toward the diaphyseal ends. No significant relationship was found between centroid displacement and body size variables with the exception of the most distal section location. Results support a relationship between tibial curvature and cross‐sectional geometry among prehistoric Central European agricultural populations, and suggest that changes in mechanical loading may have influenced a suite of morphological features related to bone adaptation in the lower limb. Am J Phys Anthropol 157:260–275, 2015. © 2015 Wiley Periodicals, Inc.     

Bone density and the lightweight skeletons of birds

The skeletons of birds are universally described as lightweight as a result of selection for minimizing the energy required for flight. From a functional perspective, the weight (mass) of an animal relative to its lift-generating surfaces is a key determinant of the metabolic cost of flight. The evolution of birds has been characterized by many weight-saving adaptations that are reflected in bone shape, many of which strengthen and stiffen the skeleton. Although largely unstudied in birds, the material properties of bone tissue can also contribute to bone strength and stiffness. In this study, I calculated the density of the cranium, humerus and femur in passerine birds, rodents and bats by measuring bone mass and volume using helium displacement. I found that, on average, these bones are densest in birds, followed closely by bats. As bone density increases, so do bone stiffness and strength. Both of these optimization criteria are used in the design of strong and stiff, but lightweight, manmade airframes. By analogy, increased bone density in birds and bats may reflect adaptations for maximizing bone strength and stiffness while minimizing bone mass and volume. These data suggest that both bone shape and the material properties of bone tissue have played important roles in the evolution of flight. They also reconcile the conundrum of how bird skeletons can appear to be thin and delicate, yet contribute just as much to total body mass as do the skeletons of terrestrial mammals.     

On the optimal shape of hip implants

The success of a total hip arthroplasty is strongly related to the initial stability of the femoral component and to the stress shielding effect. In fact, for cementless stems, initial stability is essential to promote bone ingrowth into the stem coating. An inefficient primary stability is also a cause of thigh pain. In addition, the bone adaptation after the surgery can lead to an excessive bone loss and, consequently, can compromise the success of the implant. These factors depend on prosthesis design, namely on material, interface conditions and shape. Although, surgeons use stems with very different geometries, new computational tools using structural optimization methods have been used to achieve a better design in order to improve initial stability and therefore, the implant durability. In this work, a multi-criteria shape optimization process is developed to study the relationship between implants performance and geometry. The multi-criteria objective function takes into account the initial stability of the femoral stem and the effect of stress shielding on bone adaptation after the surgery. Then, the optimized stems are tested using a concurrent model for bone remodeling and osseointegration to evaluate long-term performance. Additionally, the sensitivity to misalignments is analyzed, since femoral stems are often placed in varus or valgus position. Results show that the different criteria are contradictory resulting in different characteristics for the hip stem. However, the multi-criteria formulation leads to compromise solutions, with a combination of the geometric characteristics obtained for each criterion separately.     

The variation of the cross-sectional shape in the long bones of birds and mammals

The transverse and sagittal diameters of the long bones were measured in a sample of 53 species of eutherian mammals and 36 species of birds. The scaling of the transverse and sagittal diameters of each bone to body mass was calculated. For each bone the ratio of sagittal/transverse diameter was calculated, as an expression of the cross-sectional shape of the bones. The distributions of the ratios were not significantly different from normality in all the avian bones and in the mammalian femur and tibia. In most cases, the mean of the distribution was significantly different from 1 (circular shape). The analysis shows that changes in the ratio can be caused by selective factors, considering the correlation predicted between the breaking moments and the radii, but at the same time the cross-sectional shape of mammalian and avian long bones may have a phylogenetic basis. Finally, the previous assumption of relationship between bone curvature and stress predictability, is also discussed.     

Exploring sexual dimorphism of human occipital and temporal bones through geometric morphometrics in an identified Western-European sample

Sex estimation is a paramount step of bioprofiling in both forensic anthropology and osteoarchaeology. When the pelvis is not optimally preserved, anthropologists commonly rely on the cranium to accurately estimate sex. Over the last decades, the geometric morphometric (GM) approach has been used to determine sexual dimorphism of the crania, in size and shape, overcoming some difficulties of traditional visual and metric methods. This article aims to investigate sexual dimorphism of the occipital and temporal region through GM analysis in a metapopulation of 50 Western-European identified individuals. Statistical analyses were performed to compare centroid size and shape data between sexes through the examination of distinct functional modules. Regression and Procrustes ANOVA were used to examine allometric and asymmetrical implications. Discriminant functions, combining size and shape data, were established. Significant dimorphism in size was found, with males having larger crania, confirming the major influence size has on cranial morphology. Allometric relationships were found to be statistically significant in both right and left temporal bones while shape differences between sexes were only significant on the right temporal bone. The visualization of the mean consensus demonstrated that males displayed a larger mastoid process associated with a reduced mastoid triangle and less projected occipital condyles. This exploratory study confirms that GM analysis represents an effective way to quantitatively capture shape of dimorphic structures, even on complex rounded ones such as the mastoid region. Further examination in a larger sample would be valuable to design objective visualization tools that can improve morphoscopic sex estimation methods.     

Optimal Mass Distribution Prediction for Human Proximal Femur with Bi-modulus Property

Simulation of the mass distribution in a human proximal femur is important to provide a reasonable therapy scheme for a patient with osteoporosis. An algorithm is developed for prediction of optimal mass distribution in a human proximal femur under a given loading environment. In this algorithm, the bone material is assumed to be bi-modulus, i.e., the tension modulus is not identical to the compression modulus in the same direction. With this bi-modulus bone material, a topology optimization method, i.e., modified SIMP approach, is employed to determine the optimal mass distribution in a proximal femur. The effects of the difference between two moduli on the final material distribution are numerically investigated. Numerical results obtained show that the mass distribution in bi-modular bone materials is different from that in traditional isotropic material. As the tension modulus is less than the compression modulus for bone tissues, the amount of mass required to support tension loads is greater than that required by isotropic material for the same daily activities including one-leg stance, abduction and adduction.     

Mapping ligament insertion sites onto bone surfaces in knee by co-registration of CT and digitization data

This paper describes a new methodology that enables mapping of the ligament insertion sites onto bone surfaces in the knee joint by co-registration of the data acquired using digitization and computed tomography (CT). Local coordinate systems on the distal femur and proximal tibia were established by three spherical fiducial markers rigidly affixed to each bone. The fiducial marker centroid locations were identified by a least-squares sphere-fitting algorithm. An optimization correction procedure was proposed to mitigate the effect of the target registration error (TRE) on the alignment of coordinate systems for co-registration. A test with four cadaveric specimens demonstrated successful mapping of insertion sites for five ligaments. Fiducial registration error (FRE) as measured by the differences in inter-marker distances between the two modalities was smaller than 2%. The optimization procedure corrected the insertion site invisibility or partial visibility problem and improved the overall mapping quality, as indicated by substantial reduction of the mean and dispersion of distances from digitized insertion site points to the bone surfaces.     

Continuing periosteal apposition. II: The significance of peak bone mass, strain equilibrium, and age-related activity differentials for mechanical compensation in human tubular bones

It is generally presumed that compensation for the reduction of bone strength by progressive endosteal bone loss in adults is provided by continuing periosteal apposition (CPA) of new lamellar bone. However, the appropriate magnitude of compensatory bone growth, and the parameters that operate to determine that magnitude, are unknown. This paper examines the mechanical compensation hypothesis in a series of right-circular tubular bone analogues. Under this hypothesis, the stated objective of CPA is maintenance of the cross-sectional geometric properties of the element. These include the second and polar moments of area, as well as the cortical area of the section (I, J, and CA, respectively). This study assumes that, as resorption and apposition proceed, geometric change is isometric (shape preserving). The analysis suggests that for a given rate of endosteal bone loss (the stimulus), the magnitude of periosteal growth (the response) required to maintain geometric strength is determined by the maximum ratio (CT0) of the radial distances from the section centroid to the endosteal and periosteal surfaces (i.e., cortical thickness prior to the onset of progressive endosteal bone loss, or peak bone mass). The analysis also indicates that, for any given individual, the amount of compensatory periosteal gain required may be very small. This is particularly true for individuals having a large CT0 and for whom the magnitude of dynamic loading imparted to the skeleton declines with advancing age. This finding is illustrated in a model that relates concepts of bone surface remodeling equilibria and age-related activity differentials.     

3D Morphometric Analysis Reveals Similar Ecomorphs for Early Kangaroos (Macropodidae) and Fanged Kangaroos (Balbaridae) from the Riversleigh World Heritage Area,Australia

Understanding feeding ecology of extinct kangaroos is fundamental to understanding the evolution of kangaroos and the Australia paleoenvironment during the Oligo-Miocene. Comparisons with extant species have suggested that the macropodiforms of the Oligo/Miocene (kangaroos and allies) from the Riversleigh World Heritage Area, northern Australia, were predominantly folivorous browsers or fungivores, unlike the majority of extant species. To further test this hypothesis, we investigate the relationship between variation in cranial and mandibular shape of extant and extinct macropodiforms and ecological factors such as diet, locomotion, and body mass using 3D geometric morphometric analysis of 42 living species and eight extinct species from two radiations (the extinct clade of Balbaridae and some early representatives of the extant Macropodidae. Dietary class (fungivore, browser, grazer, and mixed feeder) correlated strongly with variation in cranial shape (20–25% of variance explained). There was also significant association between cranial shape, and both locomotor mode and body mass. In a principal component analysis of shape variation for crania (including the shape of the molar row), Riversleigh macropodiforms cluster with extant folivorous browsers on principal components (PC) 1 and 3, providing support for previous interpretations of these species as browsing kangaroos. However, as a group and regardless of phylogenetic association, the shape centroid of extinct species differs significantly from that of extant species. Riversleigh macropodiforms cluster with regular hoppers or arboreal tree kangaroos, but this may be a result of the correlation between diet and locomotor mode in kangaroos. Their similarity to extant browsers supports previous interpretations of rainforest and woodland environments at Riversleigh during the early and middle Miocene, respectively. Procrustes ANOVA Analysis of the full shape dataset and diet also shows that diet accounts for a significant portion of variation; however, when phylogeny is taken into account these results become nonsignificant. In analyses of dentary shape, some balbarid species cluster with extant mixed feeders, although this may reflect phylogenetic differences rather than ecological signal.

     

内蒙古宁城山嘴子辽墓契丹族颅骨的人类学特征

本文对内蒙古自治区赤峰市宁城县山嘴子辽代墓地出土的古代契丹族颅骨的人类学特征进行了研究。9例男性颅骨和3例女性颅骨中的大多数标本在种族特征上均可归入现代亚洲蒙古人种中的西伯利亚(北亚)人种范畴,少数标本中的个别特征则显示出某种程度上的东亚人种或北极(东北亚)人种的影响。在若干古代和近代对比组中,山嘴子辽代契丹族的体质特征分别与汉代鲜卑族和近代蒙古族最为接近。     

The effect of exercise and nutrition on the mechanostat

In this review, we discuss the effect of increased and decreased loading and nutrition deficiency on muscle and bone mass and strength (and bone length and architecture) independently and combined. Both exercise and nutrition are integral components of the mechanostat model but both have distinctly different roles. Mechanical strain imparted by muscle action is responsible for the development of the external size and shape of the bone and subsequently the bone strength. In contrast, immobilization during growth results in reduced growth in bone length and a loss of bone strength due to large losses in bone mass (a result of endosteal resorption in cortical bone and trabecular thinning) and changes in geometry (bone shafts do not develop their characteristic shape but rather develop a rounded default shape). The use of surrogate measures for peak muscle forces acting on bone (muscle strength, size, or mass) limits our ability to confirm a cause-and-effect relationship between peak muscle force acting on bone and changes in bone strength. However, the examples presented in this review support the notion that under adequate nutrition, exercise has the potential to increase peak muscle forces acting on bone and thus can lead to a proportional increase in bone strength. In contrast, nutrition alone does not influence muscle or bone in a dose-dependent manner. Muscle and bone are only influenced when there is nutritional deficiency--and in this case the effect is profound. Similar to immobilization, the immediate effect of malnutrition is a reduction in longitudinal growth. More specifically, protein and energy malnutrition results in massive bone loss due to endosteal resorption in cortical bone and trabecular thinning. Unlike loading however, there is indirect evidence that severe malnutrition when associated with menstrual dysfunction can shift the mechanostat set point upward, thus leading to less bone accrual for a given amount of bone strain.     

Size, structure and gender: lessons about fracture risk

The differences in age-related fracture risks among men and women must reflect gender differences in the relevant variables. We are concerned here with gender differences in structural variables that relate to the size and shape of bones. As children grow, their bones grow in diameter through periosteal modeling. Studies show that radial growth is driven by mechanical forces and is not just "genetically programmed". Moving bone mass farther from the center of the diaphysis makes it more effective in resisting bending and twisting forces, and disproportionately so in comparison to changes in bone mass. Gender differences in long bone structure appear to arise because the bone cells of males and females function in different hormonal environments which affect their responses to mechanical loading. In girls, bone formation on the metacarpal periosteal surface essentially stops at puberty, and is replaced by formation on the endosteal surface, reducing endosteal diameter until about age 20. Bone strength is 60% greater in male metacarpals than in those of females because bone is added periosteally in boys and endosteally in girls. At menopause endosteal resorption resumes, accompanied by slow periosteal apposition, weakening cortical structure. Similar phenomena occur in such critical regions as the femoral neck. Another fundamental gender difference in skeletal development is that whole body bone mineral content increases in linear proportion to lean body mass throughout skeletal maturation in boys, but in girls there is a distinct increase in the slope of this relationship at puberty, when estrogen rises. Frost's hypothesis is that this reflects an effect of estrogen on bone's mechanostat set point, and this is increasingly supported by data showing that estrogen and mechanical strain act through a common pathway in osteoblast-like cells. If Frost's hypothesis is correct, the mechanostat is set for maximal effect of mechanical loading on bone gain during the 2-3 years preceding menarche. During the childbearing years, the set point is at an intermediate level, and at menopause, it shifts again to place the skeleton into the metabolic equivalent of a disuse state. The most direct approach to resolving this problem would be to simulate the putative effect of estrogen on the set point itself.     

Osteoblast recruitment to sites of bone formation in skeletal development,homeostasis, and regeneration

During endochondral bone development, bone‐forming osteoblasts have to colonize the regions of cartilage that will be replaced by bone. In adulthood, bone remodeling and repair require osteogenic cells to reach the sites that need to be rebuilt, as a prerequisite for skeletal health. A failure of osteoblasts to reach the sites in need of bone formation may contribute to impaired fracture repair. Conversely, stimulation of osteogenic cell recruitment may be a promising osteo‐anabolic strategy to improve bone formation in low bone mass disorders such as osteoporosis and in bone regeneration applications. Yet, still relatively little is known about the cellular and molecular mechanisms controlling osteogenic cell recruitment to sites of bone formation. In vitro, several secreted growth factors have been shown to induce osteogenic cell migration. Recent studies have started to shed light on the role of such chemotactic signals in the regulation of osteoblast recruitment during bone remodeling. Moreover, trafficking of osteogenic cells during endochondral bone development and repair was visualized in vivo by lineage tracing, revealing that the capacity of osteoblast lineage cells to move into new bone centers is largely confined to undifferentiated osteoprogenitors, and coupled to angiogenic invasion of the bone‐modeling cartilage intermediate. It is well known that the presence of blood vessels is absolutely required for bone formation, and that a close spatial and temporal relationship exists between osteogenesis and angiogenesis. Studies using genetically modified mouse models have identified some of the molecular constituents of this osteogenic–angiogenic coupling. This article reviews the current knowledge on the process of osteoblast lineage cell recruitment to sites of active bone formation in skeletal development, remodeling, and repair, considering the role of chemo‐attractants for osteogenic cells and the interplay between osteogenesis and angiogenesis in the control of bone formation. Birth Defects Research (Part C) 99:170–191, 2013. © 2013 Wiley Periodicals, Inc.