In vivo behavior of calcium phosphate scaffolds with four different pore sizes

The goal of the present study was to assess the effect of macropore size on the in vivo behavior of ceramic scaffolds. For that purpose, beta-tricalcium phosphate (beta-TCP) cylinders with four different macropore sizes (150, 260, 510, and 1220 microm) were implanted into drill hole defects in cancellous bone of sheep and their resorption behavior was followed for 6, 12 and 24 weeks. The scaffolds were evaluated for biocompatibility, and new bone formation was observed macroscopically, histologically and histomorphometrically. Histomorphometrical measurements were performed for the whole defect area and for the area subdivided into three concentric rings (outer, medial, and inner ring). All implants were tolerated very well as evidenced by the low amount of inflammatory cells and the absence of macroscopic signs of inflammation. Resorption proceeded fast since less than 5% ceramic remained at 24-week implantation. Hardly any effect of macropore size was observed on the in vivo response. Samples with an intermediate macropore size (510 microm) were resorbed significantly faster than samples with smaller macropore sizes (150 and 260 microm). However, this fast resorption was associated with a lower bone content and a higher soft tissue content. At 12 and 24 weeks, the latter differences had disappeared. Bone was more abundant in the outer ring than in the rest of the blocks at 6 weeks, and in the outer and medial ring compared to the inner ring at 12 weeks.

Contribution of in vivo structural measurements and load/strength ratios to the determination of forearm fracture risk in postmenopausal women

Bone structure, strength and load-strength ratios contribute to forearm fracture risk independently of areal BMD.

INTRODUCTION

Technological and conceptual advances provide new opportunities for evaluating the contribution of bone density, structure, and strength to the pathogenesis of distal forearm fractures.

MATERIALS AND METHODS

From an age-sratified random sample of Rochester, MN, women, we compared 18 with a distal forearm fracture (cases) to 18 age-matched women with no osteoporotic fracture (controls). High-resolution pQCT was used to assess volumetric BMD (vBMD), geometry, and microstructure at the ultradistal radius, the site of Colles' fractures. Failure loads in the radius were estimated from microfinite element (microFE) models derived from pQCT. Differences between case and control women were assessed, and the risk of fracture associated with each variable was estimated by logistic regression analysis.

RESULTS

Given similar heights, estimated loading in a fall on the outstretched arm was the same in cases and control. However, women with forearm fractures had inferior vBMD, geometry, microstructure, and estimated bone strength. Relative risks for the strongest determinant of fracture in each of the five main variable categories were as follows: BMD (total vBMD: OR per SD change, 4.2; 95% CI, 1.4-12), geometry (cortical thickness: OR, 4.0; 95% CI, 1.4-11), microstructure (trabecular number: OR, 2.3; 95% CI, 1.02-5.1), and strength (axial rigidity: OR, 3.8; 95% CI, 1.4-10); the factor-of-risk (fall load/microFE failure load) was 24 % greater (worse) in cases (OR, 3.0; 95% CI, 1.2-7.5). Areas under ROC curves ranged from 0.72 to 0.82 for these parameters.

CONCLUSIONS

Bone geometry, microstructure, and strength contribute to forearm fractures, as does BMD, and these additional determinants of risk promise greater insights into fracture pathogenesis.

Aging, Osteocytes, and Mechanotransduction

PURPOSE OF REVIEW

The bone is able to adapt its structure to mechanical signals via the bone remodeling process governed by mechanosensitive osteocytes. With aging, an imbalance in bone remodeling results in osteoporosis. In this review, we hypothesized that changes in lacunar morphology underlie the decreased bone mechanoresponsiveness to mechanical loading with aging.

RECENT FINDINGS

Several studies have reported considerable variations in the shape of osteocytes and their lacunae with aging. Since osteocytes can sense matrix strain directly via their cell bodies, the variations in osteocyte morphology may cause changes in osteocyte mechanosensitivity. As a consequence, the load-adaptive response of osteocytes may change with aging, even when mechanical loading would remain unchanged. Though extensive quantitative data is lacking, evidence exists that the osteocyte lacunae are becoming smaller and more spherical with aging. Future dedicated studies might reveal whether these changes would affect osteocyte mechanosensation and the subsequent biological response, and whether this is (one of) the pathways involved in age-related bone loss.

Local delivery of bisphosphonate from coated orthopedic implants increases implants mechanical stability in osteoporotic rats

Patients with osteoporosis and joint disabilities represent a constant growing and challenging population to be treated in the musculoskeletal clinical field. Especially in the case of total hip arthroplasty, new solutions should be developed to compensate for the double negative factors, peri-implant osteolysis, and osteoporotic bone loss, affecting the quality of implant outcome. The goal of this study was then to establish a proof of concept for orthopedic implant used as Zoledronate delivery in osteoporotic rats, and in particular, to verify if this approach could increase the initial implant stability. Twenty-five female 6-month-old Wistar rats were ovariectomized 6 weeks before the implantation to induce osteoporosis. The animals were randomly separated in five groups representing the different Zoledronate concentrations in the HA coating: 0, 0.2, 2.1, 8.5, and 16 microg/implant. Histomorphometric measures and peri-implant bone volume fraction were assessed and mechanical stability tests were performed. Bone volume fraction and biomechanical results clearly illustrate the positive effect of Zoledronate coated implants in the osteoporotic rats. A remarkable result was to show the existence of a window of Zoledronate content (0.2 to 8.5 microg/implant) in which the mechanical fixation of the implant increased. We were able to establish the proof of concept for orthopedic implants used as a drug delivery system in osteoporotic rats. The local bisphosphonate delivery from a calcium phosphate coating allowed increase of the mechanical fixation of an orthopedic implant. This study shows that orthopedic implants containing bisphosphonates could be beneficial for osteoporotic patients in need of a total joint replacement.

Nondestructive micro-computed tomography for biological imaging and quantification of scaffold–bone interaction in vivo

Scaffolds, also called bioscaffolds, are needed in all tissue engineering applications as carriers for cells and biochemical factors, as constructs providing appropriate mechanical conditions, or as a combination of the two. The aim of this paper is to present recent developments in micro-computed tomography (microCT) analyses of scaffolds. The focus will be on imaging and quantification aspects in bone research, and will deal with the assessment of scaffold architecture and how it interacts with bone tissue. We show that micro-architectural imaging is a nondestructive and noninvasive procedure that allows a precise three-dimensional (3D) measurement of scaffold architecture. Direct microCT-based image analysis allows to accurately quantify scaffold porosity, surface area, and 3D measures such as pore size, pore distribution, and strut thickness; furthermore, it allows for a precise measurement of bone growth into the scaffold and onto its surface. This methodology is useful for quality control of scaffold fabrication processes, to assess scaffold degradation kinetics, and to assess bone tissue response. Even more so, in combination with bioreactors or in vivo animal models, microCT allows to qualitatively and quantitatively assess the spatial and temporal mineralization of bone tissue formation in scaffolds; such longitudinal studies improve the assessment of bone response due to scaffold architecture. Computational models will be helpful in further analyses of these data in order to improve our understanding of mechanical and biochemical stimuli on bone formation, and are likely to provide valuable knowledge to optimize scaffold design.

Computational finite element bone mechanics accurately predicts mechanical competence in the human radius of an elderly population

High-resolution peripheral quantitative computed tomography (HR-pQCT) is clinically available today and provides a non-invasive measure of 3D bone geometry and micro-architecture with unprecedented detail. In combination with microarchitectural finite element (μFE) models it can be used to determine bone strength using a strain-based failure criterion. Yet, images from only a relatively small part of the radius are acquired and it is not known whether the region recommended for clinical measurements does predict forearm fracture load best. Furthermore, it is questionable whether the currently used failure criterion is optimal because of improvements in image resolution, changes in the clinically measured volume of interest, and because the failure criterion depends on the amount of bone present. Hence, we hypothesized that bone strength estimates would improve by measuring a region closer to the subchondral plate, and by defining a failure criterion that would be independent of the measured volume of interest. To answer our hypotheses, 20% of the distal forearm length from 100 cadaveric but intact human forearms was measured using HR-pQCT. μFE bone strength was analyzed for different subvolumes, as well as for the entire 20% of the distal radius length. Specifically, failure criteria were developed that provided accurate estimates of bone strength as assessed experimentally. It was shown that distal volumes were better in predicting bone strength than more proximal ones. Clinically speaking, this would argue to move the volume of interest for the HR-pQCT measurements even more distally than currently recommended by the manufacturer. Furthermore, new parameter settings using the strain-based failure criterion are presented providing better accuracy for bone strength estimates.

Importance of individual rods and plates in the assessment of bone quality and their contribution to bone stiffness

Local morphometry based on the assessment of individual rods and plates was applied to 42 human vertebral trabecular bone samples. Results showed that multiple linear regression models based on local morphometry as a measure for bone microstructure helped improving our understanding of the role of local structural changes in the determination of bone stiffness as assessed from direct and computational biomechanics.

INTRODUCTION

In a recent study, we proposed a method for local morphometry of trabecular bone, i.e., morphometry as applied to individual rods and plates. In this study, we used this method to study the relative importance of local morphometry in the assessment of bone architecture and its relative contribution to the stiffness of human vertebral bone.

MATERIALS AND METHODS

We extracted 42 human trabecular bone autopsies from nine intact spinal columns. The cylindrical samples were imaged with muCT to assess bone microstructure. From these images, global and local morphometric indices were derived and related to Young's modulus as assessed by experimental uniaxial compression testing (Emeas) and computational finite element analysis (EFE).

RESULTS

We found the best single predictor for Young's modulus to be apparent bone volume density (BV/TV), which explained 89% of the variance in EFE when fitted with a power law. A multiple linear regression model combining mean trabecular spacing (Tb.Sp), mean slenderness of the rods (<Ro.Sl>), and the relative amount of rod volume to total bone volume (Ro.BV/BV) was able to explain 90% of the variance in EFE. This model could not be improved by adding BV/TV as an independent variable. Furthermore, we found that mean trabecular thickness of the rods was significantly related to EFE (r2 = 0.42), whereas mean trabecular thickness of plates had no correlation to Young's modulus. Because the globally determined trabecular thickness does not discriminate between rods and plates, this index had only a poor predictive power for EFE (r2 = 0.09), showing the importance of local analysis of individual rods and plates.

CONCLUSIONS

From these results, we conclude that models based on local morphometry help improving our understanding of the relative importance of local structural changes in the determination of the stiffness of bone. Separate analysis of individual rods and plates may help to better predict age and disease-related fractures as well as to shed new light on the effect of pharmaceutical intervention in the prevention of such fractures beyond BMD.

Differential regulation of bone and body composition in male mice with combined inactivation of androgen and estrogen receptor-alpha

Osteoporosis and muscle frailty are important health problems in elderly men and may be partly related to biological androgen activity. This androgen action can be mediated directly through stimulation of the androgen receptor (AR) or indirectly through stimulation of estrogen receptor-alpha (ERalpha) following aromatization of androgens into estrogens. To assess the differential action of AR and ERalpha pathways on bone and body composition, AR-ERalpha double-knockout mice were generated and characterized. AR disruption decreased trabecular bone mass, whereas ERalpha disruption had no additional effect on the AR-dependent trabecular bone loss. In contrast, combined AR and ERalpha inactivation additionally reduced cortical bone and muscle mass compared with either AR or ERalpha disruption alone. ERalpha inactivation--in the presence or absence of AR--increased fat mass. We demonstrate that AR activation is solely responsible for the development and maintenance of male trabecular bone mass. Both AR and ERalpha activation, however, are needed to optimize the acquisition of cortical bone and muscle mass. ERalpha activation alone is sufficient for the regulation of fat mass. Our findings clearly define the relative importance of AR and ERalpha signaling on trabecular and cortical bone mass as well as body composition in male mice.

Implant stability is affected by local bone microstructural quality

It is known that low bone quality, caused for instance by osteoporosis, not only increases the risk of fractures, but also decreases the performance of fracture implants; yet the specific mechanisms behind this phenomenon are still largely unknown. We hypothesized that especially peri-implant bone microstructure affects implant stability in trabecular bone, to a greater degree than more distant bone. To test this hypothesis we performed a computational study on implant stability in trabecular bone. Twelve humeral heads were measured using micro-computed tomography. Screws were inserted digitally into these heads at 25 positions. In addition, at each screw location, a virtual biopsy was taken. Bone structural quality was quantified by morphometric parameters. The stiffness of the 300 screw-bone constructs was quantified as a measure of implant stability. Global bone density correlated moderately with screw-bone stiffness (r2=0.52), whereas local bone density was a very good predictor (r2=0.91). The best correlation with screw-bone stiffness was found for local bone apparent Young's modulus (r2=0.97), revealing that not only bone mass but also its arrangement in the trabecular microarchitecture are important for implant stability. In conclusion, we confirmed our hypothesis that implant stability is affected by the microstructural bone quality of the trabecular bone in the direct vicinity of the implant. Local bone density was the best single morphometric predictor of implant stability. The best predictability was provided by the mechanical competence of the peri-implant bone. A clinical implication of this work is that apparently good bone stock, such as assessed by DXA, does not guarantee good local bone quality, and hence does not guarantee good implant stability. New tools that could quantify the structural or mechanical quality of the peri-implant bone may help improve the surgical intervention in reaching better clinical outcomes for screw fixation.

Evaluation of an in situ formed synthetic hydrogel as a biodegradable membrane for guided bone regeneration

The aim of the present study was to test whether or not the application of an in situ formed synthetic hydrogel made of polyethylene glycol (PEG) used as a biodegradable membrane for guided bone regeneration will result in the same amount of bone regeneration as with the use of an expanded polytetrafluoro-ethylene (ePTFE) membrane. In eight New Zealand White rabbits, four evenly distributed 6 mm diameter defects were drilled into the calvarial bone. Three treatment modalities were evenly distributed among the 32 defects: hydroxyapatite (HA)/tricalciumphosphate (TCP) granules covered at the outer and inner surface with a PEG membrane (test), HA/TCP granules covered at the outer and inner surface with an ePTFE membrane (positive control) and HA/TCP granules alone without membranes (negative control). After 4 weeks, the animals were sacrificed and the calvarial bones were removed. The area fraction of newly formed bone was determined by histomorphometrical analysis of the vertical sections from the middle of the defect and by micro-computed tomography of the entire defect. Multiple regression analysis (SAS GLM) was used to model the amount of new bone formation. The quantitative histomorphometric analysis clearly revealed higher values of newly formed bone for the two membrane groups compared with the negative control group. The average area fractions of newly formed bone measured within the former defect amounted to 20.3+/-9.5% for the PEG membrane, 18.9+/-9.9% for the ePTFE membrane, and 7.3+/-5.3% for the sites with no membrane. The micro-computed tomography also showed higher values of new bone formation for the PEG and for the ePTFE groups compared with the negative control group. The GLM revealed a highly significant effect of the treatment on the amount of bone formation (P=0.0048). The values for the negative control group were significantly lower than the ones found in the PEG membrane group (P=0.0017), whereas the ePTFE membrane group showed no significant difference from the PEG membrane group. It is concluded that the PEG membrane can be used successfully as a biodegradable barrier membrane in the treatment of non-critical-size defects in the rabbit skull, and leads to similar amounts of bone regeneration as an ePTFE membrane.

Non-invasive bone competence analysis by high-resolution pQCT: an in vitro reproducibility study on structural and mechanical properties at the human radius

Osteoporosis is defined as a skeletal disorder characterized by compromised bone strength. Bone strength depends, among others, on bone density, bone geometry and its internal architecture. With the recent introduction of a new generation high-resolution 3D peripheral quantitative computed tomography (HR-pQCT) system, direct quantification of structural bone parameters has become feasible. Furthermore, it has recently been demonstrated that bone mechanical competence can be derived from HR-pQCT based micro-finite element modeling (microFE). However, reproducibility data for HR-pQCT-derived mechanical indices is not well-known. Therefore, the aim of this study was to quantify reproducibility of HR-pQCT-derived indices. We measured 14 distal formalin-fixed cadaveric forearms three times and analyzed three different regions for each measurement. For each region cortical and trabecular parameters were determined. Reproducibility was assessed with respect to precision error (PE) and intraclass correlation coefficient (ICC). Reproducibility values were found to be best in all three regions for the full bone compartment with an average PE of 0.79%, followed by the cortical compartment (PE=1.19%) and the trabecular compartment with an average PE of 2.31%. The mechanical parameters showed similar reproducibility (PE=0.48%-2.93% for bone strength and stiffness, respectively). ICC showed a very high reproducibility of subject-specific measurements, ranging from 0.982 to 1.000, allowing secure identification of individual donors ranging from healthy to severely osteoporotic subjects. From these in vitro results we conclude that HR-pQCT derived morphometric and mechanical parameters are highly reproducible such that differences in bone structure and strength can be detected with a reproducibility error smaller than 3%; hence, the technique has a high potential to become a tool for detecting bone quality and bone competence of individual subjects.

Regional, age and gender differences in architectural measures of bone quality and their correlation to bone mechanical competence in the human radius of an elderly population

An accurate prediction of bone strength in the human radius is of major interest because distal radius fractures are amongst the most common in humans. The objective of this study was to determine gender and age-related changes in bone morphometry at the radius and how these relate to bone strength. Specifically, our aims were to (i) analyze gender differences to get an insight into different bone quantities and qualities between women and men, (ii) to determine which microarchitectural bone parameters would best correlate with strength, (iii) to find the region of interest for the best assessment of bone strength, and (iv) to determine how loss of bone quality depends on age. Intact right forearms of 164 formalin-fixed cadavers from a high-risk elderly population were imaged with a new generation high-resolution pQCT scanner (HR-pQCT). Morphometric indices were derived for six different regions and were related to failure load as assessed by experimental uniaxial compression testing. Significant gender differences in bone quantity and quality were found that correlated well with measured failure load. The most relevant region to determine failure load based on morphometric indices assessed in this study was located just below the proximal end of the subchondral plate; this region differed from the one measured clinically today. Trends in bone changes with increasing age were found, even though for all morphometric indices the variation between subjects was large in comparison to the observed age-related changes. We conclude that HR-pQCT systems can determine how gender and age-related changes in morphometric parameters relate to bone strength, and that HR-pQCT is a promising tool for the assessment of bone quality in patient populations.

Tissue modulus calculated from beam theory is biased by bone size and geometry: implications for the use of three-point bending tests to determine bone tissue modulus

Current practice to determine bone tissue modulus of murine cortical bone is to estimate it from three-point bending tests, using Euler-Bernoulli beam theory. However, murine femora are not perfect beams; hence, results can be inaccurate. Our aim was to assess the accuracy of beam theory, which we tested for two commonly used inbred strains of mice, C57BL/6 (B6) and C3H/He (C3H). We measured the three-dimensional structure of male and female B6 and C3H femora (N=20/group) by means of micro-computed tomography. For each femur five micro-finite element (micro-FE) models were created that simulated three-point bending tests with varying distances between the supports. Tissue modulus was calculated from beam theory using micro-FE results. The accuracy of beam theory was assessed by comparing the beam theory-derived moduli with the modulus as used in the micro-FE analyses. An additional set of fresh-frozen femora (10 B6 and 12 C3H) was biomechanically tested and subjected to the same micro-FE analyses. These combined experimental-computational analyses enabled an unbiased assessment of specimen-specific tissue modulus. We found that by using beam theory, tissue modulus was underestimated for all femora. Femoral geometry and size had strong effects on beam theory-derived tissue moduli. Owing to their relatively thin cortex, underestimation was markedly higher for B6 than for C3H. Underestimation was dependent on support width in a strain-specific manner. From our combined experimental-computational approach we calculated tissue moduli of 12.0+/-1.3 GPa and 13.4+/-2.1 GPa for B6 and C3H, respectively. We conclude that tissue moduli in murine femora are strongly underestimated when calculated from beam theory. Using image-based micro-FE analyses we could precisely quantify this underestimation. We showed that previously reported murine inbred strain-specific differences in tissue modulus are largely an effect of geometric differences, not accounted for by beam theory. We suggest a re-evaluation of the tissue properties obtained from three-point bending tests, especially in mouse genetics.

Abnormal bone microarchitecture and evidence of osteoblast dysfunction in premenopausal women with idiopathic osteoporosis

CONTEXT

Idiopathic osteoporosis (IOP) in premenopausal women is an uncommon disorder of uncertain pathogenesis in which fragility fractures occur in otherwise healthy women with intact gonadal function. It is unclear whether women with idiopathic low bone mineral density and no history of fragility fractures have osteoporosis.

OBJECTIVE

The objective of the study was to elucidate the microarchitectural and remodeling features of premenopausal women with IOP.

DESIGN

We performed transiliac biopsies after tetracycline labeling in 104 women: 45 with fragility fractures (IOP), 19 with idiopathic low bone mineral density (Z score ≤-2.0) and 40 controls. Biopsies were analyzed by two-dimensional quantitative histomorphometry and three-dimensional microcomputed tomography. Bone stiffness was estimated using finite element analysis.

RESULTS

Compared with controls, affected women had thinner cortices; fewer, thinner, more widely separated, and heterogeneously distributed trabeculae; reduced stiffness; and lower osteoid width and mean wall width. All parameters were indistinguishable between women with IOP and idiopathic low bone mineral density. Although there were no group differences in dynamic histomorphometric remodeling parameters, serum calciotropic hormones, bone turnover markers, or IGF-I, subjects in the lowest tertile of bone formation rate had significantly lower osteoid and wall width, more severely disrupted microarchitecture, lower stiffness, and higher serum IGF-I than those in the upper two tertiles, suggesting that women with low turnover IOP have osteoblast dysfunction with resistance to IGF-I. Subjects with high bone turnover had significantly higher serum 1,25 dihydroxyvitamin D levels and a nonsignificant trend toward higher serum PTH and urinary calcium excretion.

CONCLUSIONS

These results suggest that the diagnosis of IOP should not require a history of fracture. Women with IOP may have high, normal or low bone turnover; those with low bone turnover have the most marked deficits in microarchitecture and stiffness. These results also suggest that the pathogenesis of idiopathic osteoporosis is heterogeneous and may differ according to remodeling activity.

Load sharing and ligament strains in balanced, overstuffed and understuffed UKA. A validated finite element analysis

The aim of this study was to quantify the effects of understuffing and overstuffing UKA on bone stresses, load distribution and ligament strains. For that purpose, a numerical knee model of a cadaveric knee was developed and was validated against experimental measurements on that same knee. Good agreement was found among the numerical and experimental results. This study showed that, even if a medial UKA is well-aligned with normal soft tissue tension and with correct thickness of the tibia component, it induces a stiffness modification in the joint that alters the load distribution between the medial and lateral compartments, the bone stress and the ligament strain potentially leading to an osteoarthritic progression.

The different contributions of cortical and trabecular bone to implant anchorage in a human vertebra

The quality of the peri-implant bone and the strength of the bone-implant interface are important factors for implant anchorage. With regard to peri-implant bone, cortical and trabecular compartments both contribute to the load transfer from the implant to the surrounding bone but their relative roles have yet to be investigated in detail. However, this knowledge is crucial for the better understanding of implant failure and for the development of new implants. This is especially true for osteoporotic bone, which is characterized by a deterioration of the trabecular architecture and a thinning of the cortical shell, leading to a higher probability of implant loosening. The aim of this study was to investigate the relative biomechanical roles of cortical and trabecular bone on implant pull-out stiffness in human vertebrae. The starting point of our investigation was a micro-computed tomography scan of an adult human vertebra. The cortical shell was identified and an implant was digitally inserted into the vertebral body. Pull-out tests were simulated with micro-finite element analysis and the apparent stiffness of the system with various degrees of shell thickness and bone volume fraction was computed. Our computational models demonstrated that cortical bone, although being very thin, plays a major role in the mechanical competence of the bone-implant construct.

Age-related changes in female mouse cortical bone microporosity

Osteocyte lacunae are small cavities within the bone matrix. Their dimensions and spatial arrangement affect bone mechanical properties. Furthermore, their size and shape affect the strain in bone tissue close to the lacunae; hence, they are supposed to affect the mechanosensory function of the osteocytes residing in the lacunae. It was the purpose of this study to quantify the morphological features of osteocyte lacunae, whether these are affected by aging and whether these vary among different anatomical location. In addition, we aimed at quantifying the vascular canals as these affect bone's microporosity too. We quantified the microporosity in the fibular midshaft of young-adult and old female C57BL/6 mice. Using micro-computed tomography (μCT), we found that advanced age was associated with a significantly decreased vascular canal number and volume density. In aged mice, the mean volume of the lacuna was significantly smaller than in young animals and they were more round. Lacuna number density close to the neutral axis of the fibula was higher in older mice than in young ones. The characterization of bone microporosity presents a first step in further unraveling their potential role in age-related reductions in bone strength.

Technical Note: Cortical thickness and density estimation from clinical CT using a prior thickness-density relationship

PURPOSE

Cortical thickness and density are critical components in determining the strength of bony structures. Computed tomography (CT) is one possible modality for analyzing the cortex in 3D. In this paper, a model-based approach for measuring the cortical bone thickness and density from clinical CT images is proposed.

METHODS

Density variations across the cortex were modeled as a function of the cortical thickness and density, location of the cortex, density of surrounding tissues, and imaging blur. High resolution micro-CT data of cadaver proximal femurs were analyzed to determine a relationship between cortical thickness and density. This thickness-density relationship was used as prior information to be incorporated in the model to obtain accurate measurements of cortical thickness and density from clinical CT volumes. The method was validated using micro-CT scans of 23 cadaver proximal femurs. Simulated clinical CT images with different voxel sizes were generated from the micro-CT data. Cortical thickness and density were estimated from the simulated images using the proposed method and compared with measurements obtained using the micro-CT images to evaluate the effect of voxel size on the accuracy of the method. Then, 19 of the 23 specimens were imaged using a clinical CT scanner. Cortical thickness and density were estimated from the clinical CT images using the proposed method and compared with the micro-CT measurements. Finally, a case-control study including 20 patients with osteoporosis and 20 age-matched controls with normal bone density was performed to evaluate the proposed method in a clinical context.

RESULTS

Cortical thickness (density) estimation errors were 0.07 ± 0.19 mm (-18 ± 92 mg/cm(3)) using the simulated clinical CT volumes with the smallest voxel size (0.33 × 0.33 × 0.5 mm(3)), and 0.10 ± 0.24 mm (-10 ± 115 mg/cm(3)) using the volumes with the largest voxel size (1.0 × 1.0 × 3.0 mm(3)). A trend for the cortical thickness and density estimation errors to increase with voxel size was observed and was more pronounced for thin cortices. Using clinical CT data for 19 of the 23 samples, mean errors of 0.18 ± 0.24 mm for the cortical thickness and 15 ± 106 mg/cm(3) for the density were found. The case-control study showed that osteoporotic patients had a thinner cortex and a lower cortical density, with average differences of -0.8 mm and -58.6 mg/cm(3) at the proximal femur in comparison with age-matched controls (p-value < 0.001).

CONCLUSIONS

This method might be a promising approach for the quantification of cortical bone thickness and density using clinical routine imaging techniques. Future work will concentrate on investigating how this approach can improve the estimation of mechanical strength of bony structures, the prevention of fracture, and the management of osteoporosis.

Mechanical Loading Differentially Affects Osteocytes in Fibulae from Lactating Mice Compared to Osteocytes in Virgin Mice: Possible Role for Lacuna Size

Hormonal changes during lactation are associated with profound changes in bone cell biology, such as osteocytic osteolysis, resulting in larger lacunae. Larger lacuna shape theoretically enhances the transmission of mechanical signals to osteocytes. We aimed to provide experimental evidence supporting this theory by comparing the mechanoresponse of osteocytes in the bone of lactating mice, which have enlarged lacunae due to osteocytic osteolysis, with the response of osteocytes in bone from age-matched virgin mice. The osteocyte mechanoresponse was measured in excised fibulae that were cultured in hormone-free medium for 24 h and cyclically loaded for 10 min (sinusoidal compressive load, 3000 µε, 5 Hz) by quantifying loading-related changes in Sost mRNA expression (qPCR) and sclerostin and β-catenin protein expression (immunohistochemistry). Loading decreased Sost expression by ~ threefold in fibulae of lactating mice. The loading-induced decrease in sclerostin protein expression by osteocytes was larger in lactating mice (55% decrease ± 14 (± SD), n = 8) than virgin mice (33% decrease ± 15, n = 7). Mechanical loading upregulated β-catenin expression in osteocytes in lactating mice by 3.5-fold (± 0.2, n = 6) which is significantly (p < 0.01) higher than the 1.6-fold increase in β-catenin expression by osteocytes in fibulae from virgin mice (± 0.12, n = 4). These results suggest that osteocytes in fibulae from lactating mice with large lacunae may respond stronger to mechanical loading than those from virgin mice. This could indicate that osteocytes residing in larger lacuna show a stronger response to mechanical loading.

In silico 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.