Michael Doube

BoneJ: Free and extensible bone image analysis in ImageJ

Bone geometry is commonly measured on computed tomographic (CT) and X-ray microtomographic (μCT) images. We obtained hundreds of CT, μCT and synchrotron μCT images of bones from diverse species that needed to be analysed remote from scanning hardware, but found that available software solutions were expensive, inflexible or methodologically opaque. We implemented standard bone measurements in a novel ImageJ plugin, BoneJ, with which we analysed trabecular bone, whole bones and osteocyte lacunae. BoneJ is open source and free for anyone to download, use, modify and distribute.

In Vivo Time- Resolved Microtomography Reveals the Mechanics of the Blowfly Flight Motor

Time-resolved X-ray microtomography permits a real-time view of the blowfly in flight at a previously unprecedented level of detail, revealing how the tiny steering muscles work.

Sclerostin and the regulation of bone formation: Effects in hip osteoarthritis and femoral neck fracture

Remodeling imbalance in the elderly femoral neck can result in thin cortices and porosity predisposing to hip fracture. Hip osteoarthritis protects against intracapsular hip fracture. By secreting sclerostin, osteocytes may inhibit Wnt signaling and reduce bone formation by osteoblasts. We hypothesised that differences in osteocytic sclerostin expression might account for differences in osteonal bone‐formation activity between controls and subjects with hip fracture or hip osteoarthritis. Using specific antibody staining, we determined the osteocytic expression of sclerostin within osteons of the femoral neck cortex in bone removed from subjects undergoing surgery for hip osteoarthritis (hOA: 5 males, 5 females, 49 to 92 years of age) or hip fracture fixation (FNF: 5 males, 5 females, 73 to 87 years of age) and controls (C: 5 males, 6 females, 61 to 90 years of age). Sclerostin expression and distances of each osteocyte to the canal surface and cement line were assessed for all osteonal osteocytes in 636 unremodeled osteons chosen from fields (∼0.5 mm in diameter) with at least one canal staining for alkaline phosphatase (ALP), a marker of bone formation. In adjacent sections, ALP staining was used to classify basic multicellular unit (BMUs) as quiescent or actively forming bone (ALP+). The areal densities of scl− and scl+ osteocytes (number of cells per unit area) in the BMU were inversely correlated and were strong determinants of ALP status in the BMU. In controls and hip fracture patients only, sclerostin‐negative osteocytes were closer to osteonal surfaces than positively stained cells. Osteon maturity (progress to closure) was strongly associated with the proportion of osteonal osteocytes expressing sclerostin, and sclerostin expression was the chief determinant of ALP status. hOA patients had 18% fewer osteocytes per unit bone area than controls, fewer osteocytes expressed sclerostin on average than in controls, but wide variation was seen between subjects. Thus, in most hOA patients, there was increased osteonal ALP staining and reduced sclerostin staining of osteocytes. In FNF patients, newly forming osteons were similar in this respect to hOA osteons, but with closure, there was a much sharper reduction in ALP staining that was only partly accounted for by the increased proportions of osteonal osteocytes staining positive for sclerostin. There was no evidence for a greater effect on ALP expression by osteocytes near the osteonal canal. In line with data from blocking antibody experiments, osteonal sclerostin appears to be a strong determinant of whether osteoblasts actively produce bone. In hOA, reduced sclerostin expression likely mediates increased osteoblastic activity in the intracapsular cortex. In FNF, full osteonal closure is postponed, with increased porosity, in part because the proportion of osteocytes expressing sclerostin increases sharply with osteonal maturation.

Three-Dimensional Geometric Analysis of Felid Limb Bone Allometry

Background Studies of bone allometry typically use simple measurements taken in a small number of locations per bone; often the midshaft diameter or joint surface area is compared to body mass or bone length. However, bones must fulfil multiple roles simultaneously with minimum cost to the animal while meeting the structural requirements imposed by behaviour and locomotion, and not exceeding its capacity for adaptation and repair. We use entire bone volumes from the forelimbs and hindlimbs of Felidae (cats) to investigate regional complexities in bone allometry. Method/Principal Findings Computed tomographic (CT) images (16435 slices in 116 stacks) were made of 9 limb bones from each of 13 individuals of 9 feline species ranging in size from domestic cat (Felis catus) to tiger (Panthera tigris). Eleven geometric parameters were calculated for every CT slice and scaling exponents calculated at 5% increments along the entire length of each bone. Three-dimensional moments of inertia were calculated for each bone volume, and spherical radii were measured in the glenoid cavity, humeral head and femoral head. Allometry of the midshaft, moments of inertia and joint radii were determined. Allometry was highly variable and related to local bone function, with joint surfaces and muscle attachment sites generally showing stronger positive allometry than the midshaft. Conclusions/Significance Examining whole bones revealed that bone allometry is strongly affected by regional variations in bone function, presumably through mechanical effects on bone modelling. Bones phenotypic plasticity may be an advantage during rapid evolutionary divergence by allowing exploitation of the full size range that a morphotype can occupy. Felids show bone allometry rather than postural change across their size range, unlike similar-sized animals.

Altered lacunar and vascular porosity in osteogenesis imperfecta mouse bone as revealed by synchrotron tomography contributes to bone fragility

Osteogenesis imperfecta (brittle bone disease) is caused by mutations in the collagen genes and results in skeletal fragility. Changes in bone porosity at the tissue level indicate changes in bone metabolism and alter bone mechanical integrity. We investigated the cortical bone tissue porosity of a mouse model of the disease, oim, in comparison to a wild type (WT-C57BL/6), and examined the influence of canal architecture on bone mechanical performance. High-resolution 3D representations of the posterior tibial and the lateral humeral mid-diaphysis of the bones were acquired for both mouse groups using synchrotron radiation-based computed tomography at a nominal resolution of 700nm. Volumetric morphometric indices were determined for cortical bone, canal network and osteocyte lacunae. The influence of canal porosity architecture on bone mechanics was investigated using microarchitectural finite element (μFE) models of the cortical bone. Bright-field microscopy of stained sections was used to determine if canals were vascular. Although total cortical porosity was comparable between oim and WT bone, oim bone had more numerous and more branched canals (p<0.001), and more osteocyte lacunae per unit volume compared to WT (p<0.001). Lacunae in oim were more spherical in shape compared to the ellipsoidal WT lacunae (p<0.001). Histology revealed blood vessels in all WT and oim canals. μFE models of cortical bone revealed that small and branched canals, typical of oim bone, increase the risk of bone failure. These results portray a state of compromised bone quality in oim bone at the tissue level, which contributes to its deficient mechanical properties.

Four-dimensional in vivo X-ray microscopy with projection-guided gating

Visualizing fast micrometer scale internal movements of small animals is a key challenge for functional anatomy, physiology and biomechanics. We combine phase contrast tomographic microscopy (down to 3.3 μm voxel size) with retrospective, projection-based gating (in the order of hundreds of microseconds) to improve the spatiotemporal resolution by an order of magnitude over previous studies. We demonstrate our method by visualizing 20 three-dimensional snapshots through the 150 Hz oscillations of the blowfly flight motor.

Whole-bone scaling of the avian pelvic limb

Birds form the largest extant group of bipedal animals and occupy a broad range of body masses, from grams to hundreds of kilograms. Additionally, birds occupy distinct niches of locomotor behaviour, from totally flightless strong runners such as the ratites (moa, kiwi, ostrich) to birds that may walk, dabble on water or fly. We apply a whole‐bone approach to investigate allometric scaling trends in the pelvic limb bones (femur, tibiotarsus, tarsometatarsus) from extant and recently extinct birds of greatly different size, and compare scaling between birds in four locomotor groups; flightless, burst‐flying, dabbling and flying. We also compare scaling of birds’ femoral cross‐sectional properties to data previously collected from cats. Scaling exponents were not significantly different between the different locomotor style groups, but elevations of the scaling relationships revealed that dabblers (ducks, geese, swans) have particularly short and slender femora compared with other birds of similar body mass. In common with cats, but less pronounced in birds, the proximal and distal extrema of the bones scaled more strongly than the diaphysis, and in larger birds the diaphysis occupied a smaller proportion of bone length than in smaller birds. Cats and birds have similar femoral cross‐sectional area (CSA) for the same body mass, yet birds’ bone material is located further from the bone’s long axis, leading to higher second and polar moments of area and a greater inferred resistance to bending and twisting. The discrepancy in the relationship between outer diameter to CSA may underlie birds’ reputation for having ‘light’ bones.

structure model index does not measure rods and plates in trabecular bone

Structure model index (SMI) is widely used to measure rods and plates in trabecular bone. It exploits the change in surface curvature that occurs as a structure varies from spherical (SMI = 4), to cylindrical (SMI = 3) to planar (SMI = 0). The most important assumption underlying SMI is that the entire bone surface is convex and that the curvature differential is positive at all points on the surface. The intricate connections within the trabecular continuum suggest that a high proportion of the surface could be concave, violating the assumption of convexity and producing regions of negative differential. We implemented SMI in the BoneJ plugin and included the ability to measure the amounts of surface that increased or decreased in area after surface mesh dilation, and the ability to visualize concave and convex regions. We measured SMI and its positive (SMI+) and negative (SMI−) components, bone volume fraction (BV/TV), the fraction of the surface that is concave (CF), and mean ellipsoid factor (EF) in trabecular bone using 38 X-ray microtomography (XMT) images from a rat ovariectomy model of sex steroid rescue of bone loss, and 169 XMT images from a broad selection of 87 species’ femora (mammals, birds, and a crocodile). We simulated bone resorption by eroding an image of elephant trabeculae and recording SMI and BV/TV at each erosion step. Up to 70%, and rarely <20%, of the trabecular surface is concave (CF 0.155–0.700). SMI is unavoidably influenced by aberrations induced by SMI−, which is strongly correlated with BV/TV and CF. The plate-to-rod transition in bone loss is an erroneous observation resulting from the close and artifactual relationship between SMI and BV/TV. SMI cannot discern between the distinctive trabecular geometries typical of mammalian and avian bone, whereas EF clearly detects birds’ more plate-like trabeculae. EF is free from confounding relationships with BV/TV and CF. SMI results reported in the literature should be treated with suspicion. We propose that EF should be used instead of SMI for measurements of rods and plates in trabecular bone.

Changes in mineralised tissue at the site of origin of condylar fracture are present before athletic training in Thoroughbred horses

Abstract AIM: To show that changes are present at the site of origin of metacarpal condylar fracture in young Thoroughbred horses before they enter race training. METHOD: Bone slices, 2 mm thick, in three mediolateral planes through the centre of rotation of the metacarpo-phalangeal joint (MCPJ) of both distal third metacarpal bones (Mc3) of 12 Thoroughbred horses aged 17 months, were imaged using point-projection digital X-ray imaging (µXR). RESULTS: In some horses, linear or ovoid radiolucency was found in articular calcified cartilage (ACC) and subchondral bone of the palmaro-distal aspect of the sagittal groove, exactly at the site of more advanced stages of condylar fatigue fracture. An incidental finding was ovoid radiolucency in the apex of the dorso-distal aspect of the sagittal ridge, with or without fragmentation or disturbance of the subchondral mineralised tissue line, resembling equine osteochondrosis. CONCLUSIONS AND CLINICAL RELEVANCE: The findings imply that the aetiology of condylar fatigue fracture in young Thoroughbred horses includes abnormality in development of the bone and joint that is present before athletic activity occurs.

Osteocyte recruitment declines as the osteon fills in: Interacting effects of osteocytic sclerostin and previous hip fracture on the size of cortical canals in the femoral neck

There is little information on the distribution of osteocytes within the individual cortical osteon, but using direct 3-D imaging in a single subject, Hannah et al. found a gradient with a two-fold higher density of cells adjacent to the cement line compared to near the canal. Since a limiting factor for bone formation might be the availability of osteoblasts due to their recruitment as osteocytes, we studied distributions of osteonal osteocytes in frozen sections of the femoral neck cortex. Osteocytes were stained with an anti-sclerostin antibody and counter-stained with toluidine blue. Adjacent sections were stained for alkaline phosphatase (ALP). Each osteonal osteocyte was categorised as being sclerostin-positive (scl+) or negative (scl-). ImageJ was used to measure the perimeter and area of each osteon and canal, while special purpose routines were used to measure the minimum distances of each osteocyte from the cement line and the canal. Canal area was strongly correlated with osteon area. Osteocytes were most dense close to the cement line; and their areal density within the matrix declined up to three-fold between the cement line and the canal, depending on osteon diameter. Large and small osteons had similar densities of osteocytes close to the cement line, but fractured neck of femur cases had significantly lower densities of osteocytes close to the canal. Higher osteocyte density close to the canal was associated with ALP expression. It is concluded that entombment of osteocytes newly drawn from the osteoblast pool into the mineralising matrix is independent of preceding bone resorption depth. As osteonal infilling proceeds, osteocyte formation declines more rapidly than matrix formation, leading to a progressive reduction in osteocyte density. A shrinking supply of precursor osteoblasts due to previous osteocyte recruitment, apoptosis, or both could produce this effect. In a statistically significant contrast, sclerostin negative osteocytes adjacent to the canal had the expected effect of reducing canal size in controls but this was not seen in hip fracture. This demonstrated the failure of osteonal osteoblasts to sustain bone formation through a complete remodelling cycle in osteoporosis, perhaps due to insufficient osteoblasts remaining capable of mineralized matrix formation. The failure of osteocytic sclerostin suppression to associate with bone formation in these osteons might alternatively be explained by downstream interference with sclerostins effect on wnt signalling.