Document Detail

Mechanical implications of estrogen supplementation in early postmenopausal women.
Jump to Full Text
MedLine Citation:
PMID:  20200948     Owner:  NLM     Status:  MEDLINE    
Whereas the structural implications of drug intervention are well established, there are few data on the possible mechanical consequences of treatment. In this work we examined the changes in elastic and shear moduli (EM and SM) in a region of trabecular bone in the distal radius and distal tibia of early postmenopausal women on the basis of MRI-based micro-finite-element (microFE) analysis. Whole-section axial stiffness (AS) encompassing both trabecular and cortical compartments was evaluated as well. The study was conducted on previously acquired high-resolution images at the two anatomic sites. Images were processed to yield a 3D voxel array of bone-volume fraction (BVF), which was converted to a microFE model of hexahedral elements in which tissue modulus was set proportional to voxel BVF. The study comprised 65 early postmenopausal women (age range 45 to 55 years), of whom 32 had chosen estrogen supplementation (estradiol group); the remainder had not (control group). Subjects had been scanned at baseline and 12 and 24 months thereafter. At the distal tibia, EM and SM were reduced by 2.9% to 5.5% in the control group (p < .05 to <.005), but there was no change in the estradiol subjects. AS decreased 3.9% (4.0%) in controls (p < .005) and increased by 5.8% (6.2%) in estradiol group subjects (p < .05) at 12 (24) months. At the distal radius, EM and SM changes from baseline were not significant, but at both time points AS was increased in estradiol group subjects and decreased in controls (p < .005 to <.05), albeit by a smaller margin than at the tibia. EM and SM were strongly correlated with BV/TV (r(2) = 0.44 to 0.92) as well as with topologic parameters expressing the ratio of plates to rods (r(2) = 0.45 to 0.82), jointly explaining up to 96% of the variation in the mechanical parameters. Finally, baseline AS was strongly correlated between the two anatomic sites (r(2) = 0.58), suggesting that intersubject variations in the bone's mechanical competence follows similar mechanisms. In conclusion, the results demonstrate that micro-MRI-based microFE models are suited for the study of the mechanical implications of antiresorptive treatment. The data further highlight the anabolic effect of short-term estrogen supplementation.
Felix W Wehrli; Chamith S Rajapakse; Jeremy F Magland; Peter J Snyder
Related Documents :
20528568 - Evaluation of ugt1a1 and sult1a1 polymorphisms with lipid levels in women with differen...
9217668 - Postmenopausal estrogen replacement and tooth retention.
17493618 - Safety of testosterone treatment in postmenopausal women.
3812438 - Attributable risk ratio estimation from matched-pairs case-control data.
22791398 - Cardiovascular disease and rheumatoid arthritis: an update.
25356308 - Flexibility in men's sexual practices in response to iatrogenic erectile dysfunction af...
Publication Detail:
Type:  Journal Article; Research Support, N.I.H., Extramural; Research Support, Non-U.S. Gov't    
Journal Detail:
Title:  Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research     Volume:  25     ISSN:  1523-4681     ISO Abbreviation:  J. Bone Miner. Res.     Publication Date:  2010 Jun 
Date Detail:
Created Date:  2010-06-02     Completed Date:  2010-09-01     Revised Date:  2013-05-30    
Medline Journal Info:
Nlm Unique ID:  8610640     Medline TA:  J Bone Miner Res     Country:  United States    
Other Details:
Languages:  eng     Pagination:  1406-14     Citation Subset:  IM    
Copyright Information:
(c) 2010 American Society for Bone and Mineral Research.
Department of Radiology, University of Pennsylvania, Philadelphia, PA 19104, USA.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Biomechanics / drug effects
Dietary Supplements*
Elastic Modulus / drug effects
Estrogens / administration & dosage*,  pharmacology*
Finite Element Analysis
Image Processing, Computer-Assisted
Middle Aged
Postmenopause / drug effects*,  physiology*
Grant Support
Reg. No./Substance:

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): J Bone Miner Res
Journal ID (publisher-id): jbmr
ISSN: 0884-0431
ISSN: 1523-4681
Publisher: Wiley Subscription Services, Inc., A Wiley Company
Article Information
Download PDF
Copyright © 2010 American Society for Bone and Mineral Research
Received Day: 23 Month: 10 Year: 2009
Revision Received Day: 04 Month: 12 Year: 2009
Accepted Day: 11 Month: 1 Year: 2010
Print publication date: Month: 6 Year: 2010
Electronic publication date: Day: 29 Month: 1 Year: 2010
Volume: 25 Issue: 6
First Page: 1406 Last Page: 1414
ID: 3153138
PubMed Id: 20200948
DOI: 10.1002/jbmr.33

Mechanical Implications of Estrogen Supplementation in Early Postmenopausal Women
Felix W Wehrli1
Chamith S Rajapakse1
Jeremy F Magland1
Peter J Snyder2
1Department of Radiology, Laboratory for Structural NMR Imaging, University of PennsylvaniaPhiladelphia, PA, USA
2Department of Medicine, Division of Endocrinology, Diabetes & Metabolism, University of PennsylvaniaPhiladelphia, PA, USA
Correspondence: Address correspondence to: Felix W Wehrli, PhD, Department of Radiology, University of Pennsylvania Medical Center, Laboratory for Structural NMR Imaging, 3400 Spruce Street, 1 Founders Pavilion, Philadelphia, PA 19104, USA. E-mail:


The notion that besides bone mineral density, parameters related to architecture at the macro- and microstructural levels,(1) as well as mineral and matrix chemistry,(2) collectively summarized under the term bone quality,(3) determine overall bone strength is now generally accepted. There is also consensus that changes in many of the parameters occurring during aging and disease progression (or regression in response to treatment) can affect the bone's mechanical behavior in a manner independent of the measured changes in material density.(4)

While initially confined to measurements in small specimens, advances in imaging technology, both computed tomography (CT) and magnetic resonance imaging (MRI), now allow detailed models to be reconstructed from high-resolution 3D images of both trabecular and cortical microarchitecture acquired in vivo in patients. These and related imaging technologies and their merits and current limitations have been reviewed extensively (see, for example, refs. (57)). On the basis of such imaging data, it has become possible to derive a host of structural measures expressing properties of the trabecular network in terms of scale, topology, and orientation, known from specimen studies to be relevant to the bone's mechanical properties.(8,9) Similarly, macro- and microstructural properties can be assessed for cortical bone, such as cortical thickness and cross-sectional area, via direct quantification of these structural measures on the basis of MRI(6) or CT.(10) Microstructural parameters such as porosity have been shown to be obtainable indirectly, for example, via MRI-based quantification of bone water.(11)

There are a number of reports showing associations between structural measures at peripheral sites such as the radius or calcaneus and vertebral or femur fracture status.(1215) It also has been demonstrated recently that modern imaging technologies, notably micro-magnetic resonance imaging (µMRI), whole-body multislice CT, and high-resolution peripheral quantitative computed tomography (HR-pQCT), are able to detect age-related loss of trabecular network integrity,(16) as well as structural changes in response to drug intervention,(1720)

Whereas there is evidence that certain structural features make trabecular bone more prone to failure (eg, rodlike versus platelike(21) or increased degree of structural anisotropy(22)), only recently, capabilities have become available to identify such structural properties on the basis of in vivo images and quantify changes serially over time in human subjects.(17,18,20) However, even though empirically established associations are known between biomechanical and structural parameters, the latter largely represent surrogates for strength. Since mechanical testing in vivo is not possible, alternate approaches have been pursued to estimate strength on the basis of image-derived finite-element meshes (see, for example, ref. (23)). Micro-finite element (µFE) models established from high-resolution µCT images yielded mechanical constants in excellent agreement with those obtained by mechanical testing.(24) Further, strong correlations with actual measures of elastic moduli and image-derived values have been shown recently to be obtainable even after downsampling the µCT images to lower resolution and superposition of noise to simulate in vivo imaging conditions.(25) These and other data instill significant confidence that meaningful biomechanical data are obtainable from in vivo image-based FE models.

Early work by Muller and Ruegegger, based on simulated atrophy of images acquired at in vivo resolution by peripheral CT in a cadaver specimen, subjected to µFE computation of elastic moduli, showed that a relatively small amount of bone loss can have disproportionately large mechanical consequences.(26) The work suggested the feasibility of evaluating disease and aging-related effects on the bone's mechanical behavior. Van Rietbergen and colleagues examined the potential of µMRI to quantify the effects of drug intervention in patients treated with a selective estrogen receptor modulator.(27) A recent article showed in 10 hypogonadal subjects that 2-year treatment with testosterone resulted in significant increases in some elastic and shear moduli,(28) computed from a region of trabecular bone in the distal tibial metaphysis, commensurate with changes in topology of the trabecular network reported earlier.(18)

The preceding study is of limited scope in terms of both the number of subjects studied and the size of the imaging volume used for construction of the µFE model. In the present work, we have performed µFE analysis on MR data in the distal tibia and radius from a recent structural imaging study in early postmenopausal women who elected estrogen supplementation compared with those who did not during a 2-year period. The image data were reprocessed yielding grayscale FE models of both central trabecular subregions used for multidirectional estimation of Young's and shear moduli, as well as whole-section simulated compressive testing, including both trabecular and cortical bone compartments, at three time points.

Materials and Methods
Subjects and image acquisition

This work is based on imaging data from a recently completed study in the authors' laboratory(20) involving 65 early postmenopausal women, 45 to 55 years of age, who had been recruited for the study. Thirty-two of the women had chosen estrogen supplementation (estradiol group); the remainder had not (control group). µMRI in the metaphysis of the distal radius and distal tibia had been performed at 137 × 137 × 410 µm3 voxel size at baseline and 12 and 24 months to determine the temporal changes in the trabecular architecture. The imaging volume was 70 × 40 × 13 mm3 (scan time 12 minutes) and 70 × 50 × 13 mm3 (scan time 16 minutes) for the distal radius and distal tibia, respectively, with the third dimension representing the axial direction. Full details of the subject characteristics and image acquisition protocol are given in ref. (20).

Image processing

Images from the prior study were reprocessed using more advanced algorithms. The raw images were corrected for translational displacements during the scan owing to involuntary subject motion with the aid of navigator echo data acquired simultaneously with the image data.(29) Motion-corrected k-space data then were Fourier-transformed to yield 32 image slices. Some examinations did not provide images of sufficient quality for analysis because of artifacts caused by subject motion during scanning that could not be removed via navigator correction (eg, rotational motion) or wrap-around artifacts resulting from field-of-view restrictions for some subjects with particularly large wrists and ankles. These images were excluded from analysis by following a procedure described previously.(20,30) Further, if the common image volume between the three time points was less than 2.5 mm in the axial direction (owing to registration errors), those images also were excluded from the analysis. In summary, the number of subjects who provided distal tibia data with adequate image quality for the baseline to 12-month (24) interval was 28 (23) for the control group and 23 (19) for the estradiol group. Similarly, the corresponding numbers for the distal radius were 23 (14) and 13 (12), respectively.

The preprocessed images then were subjected to a cascade of processing steps before being used as input to the mechanical and topologic analysis (Fig. 1). First, image intensity variations across the volume produced by inhomogeneous sensitivity of the MR receive coil were corrected using a local thresholding algorithm.(31) The grayscale voxel values of the intensity-corrected images then were scaled linearly to cover the range from 0% to 100%, with pure marrow and pure bone having minimum and maximum values, respectively. We refer to the resulting 3D array as the bone-volume fraction (BVF) map, with individual voxel values representing the fraction of the voxel occupied by bone (ie, BV/TV). Subsequently, the BVF maps were subjected to two different processing paths in preparation for the full cross-sectional and subregional analysis. To enable full cross-sectional analysis, the soft tissue was first segmented out from bone by delineating the periosteal boundary using an operator-guided region-selection program developed in-house.(32) Finally, common 5-mm sections along the bone's axial direction were extracted from the images at the three time points for each subject. To prepare the images for subregional analysis, BVF maps were first sinc-interpolated by a factor of 3 in all three image coordinate axes to improve apparent resolution, yielding 45.7 × 45.7 × 137 µm3 voxels. The 12- and 24-month repeat images of each subject then were registered in 3D(33) to match the corresponding baseline volume using trilinear interpolation. Finally, 10 × 10 × 5 mm3 cuboid subvolumes were extracted from the center of the trabecular bone (TB) region for analysis.

Grayscale-based µFE modeling

To create a finite-element mesh, each voxel in the BVF map was converted directly to a hexahedral (brick) finite element with dimensions corresponding to the voxel size,(34) The bone tissue material properties were chosen as isotropic and linearly elastic, with each element's Young's modulus (YM) set to be linearly proportional to the voxel's BV/TV value such that YM = 15 GPa × BVF, while Poisson's ratio was kept constant at 0.3 for all elements. A partial threshold was applied by setting finite elements with BVF < 20% to zero to remove the contributions from image noise.

Subregional µFE analysis

To estimate the mechanical properties of the cuboid subregion from the generated µFE model, Young's and shear moduli, scalar quantities that express the relationship between stress and strain, were determined by simulating compression and shear tests on the µFE model by imposing different boundary conditions on the eight surfaces of the modeled bone region. The Young's moduli, representing the stiffness along each of the three principal axes (E11, E22, and E33) of the coordinate system (x, y, and z), were found by performing three simulated “compressive” tests along each of the three orthogonal axes.(34,35) Since the patients had been imaged in the supine position with the tibia's long axis parallel to the magnet bore, the scanner coordinate system was essentially collinear with the anatomic coordinate frame. Thus the subscripts 1, 2, and 3 correspond to anteroposterior, mediolateral, and inferosuperior directions, respectively. The value of E11, for example, was obtained by applying a small simulated displacement (compressive strain of approximately 1%) along the x direction to all nodes on one of the boundary surfaces perpendicular to the x axis in such a manner that the nodes on the opposite surface were restrained along the x direction while lateral displacements remained unrestrained. The unknown displacements at each node then were found by minimizing the total strain energy of the system. Subsequently, the resulting stress along the x direction was calculated by summing the x components for the reaction force and dividing by the surface area. Finally, E11 was estimated as the ratio of calculated stress to the applied strain. E22 and E33 also were obtained analogously. Using a similar approach, shear moduli (G23, G31, and G12), expressing the ratio of shear stress to shear strain of the modeled region, were estimated by simulating shear displacements applied to nodes on opposite surfaces along antiparallel directions.

Full cross-sectional µFE analysis

Axial stiffness of the distal radius and distal tibia was estimated via µFE analysis of the full FE mesh comprising both trabecular and cortical bone compartments generated as described earlier. First, simulated compression was applied along the bone's longitudinal axis by applying a constant displacement (∼1% strain) to all FE nodes in the proximal face of the FE mesh while keeping those in the distal face constrained.(36) The µFE system then was solved as described earlier, yielding a 3D strain map. Finally, the axial full cross-sectional modulus was obtained as the ratio of the resulting stress on the proximal face to the applied strain.

Digital-topologic analysis

In order to evaluate parameters of scale and topology of the trabecular bone network, digital-topologic analysis (DTA)(37) was performed on the same cuboid subvolumes of trabecular bone used for µFE analysis. In brief, the sinc-interpolated BVF maps obtained as described earlier were binarized such that bone and background voxels were assigned values of 1 and 0, respectively, using an empirically optimized threshold level of 20% of pure bone intensity.(38) The binary volumes then were subjected to a topology-preserving thinning algorithm(39) yielding the 3D skeleton of the trabecular bone network. Subsequently, topologic classification was performed by identifying each voxel as belonging to a surface, curve, or mutual junctions between these entities.(37) Finally, the composite topologic parameters surface-to-curve ratio (S/C) and erosion index (EI) were computed from the DTA maps of the skeleton.(37) BV/TV was computed as described in ref. (31).

Statistical analysis

The statistical analysis was analogous to the one outlined in the parent article,(20) based on multivariate analysis of variance (MANOVA). Following determination that the data were normally distributed, MANOVA was performed with time (repeated) and treatment (nonrepeated) factors along with inclusion of the treatment × time cross-product to assess second-order interaction. When a response variable was associated with significant treatment × time interaction, simple effects over time were assessed within each treatment groups using two-sided paired t tests. Between-group differences in the changes from baseline were evaluated by unpaired t tests. Intersite correlation of full cross-sectional stiffness, as well as correlations between mechanical and topologic parameters of trabecular bone, were evaluated by least-squares regression in terms of the Pearson correlation coefficient. All statistical analyses were performed using JMP Discovery Software, Version 7.02, (SAS Institute, Inc., Cary, NC), with p < .05 indicating statistical significance.

Subregional µFE analysis

The mean values of Young's and shear moduli, as well as BV/TV, all computed from matching trabecular bone subvolumes in the distal tibia at all three time points, as well as the relative changes from baseline, are summarized in Tables 1 and 2. The relative temporal changes in parameters are charted in Fig. 2. The subregional mechanical parameters did not show any significant change from baseline to the two follow-up time points in the estradiol group, suggesting that estrogen supplementation preserves the integrity of the trabecular network, as shown in some of our prior work.(20) On the other hand, in the control subjects, three of the mechanical parameters suggest significant reduction already after 12 months (3.0% to 5.5%, p = .01 to .002), whereas the change in BV/TV was not significant. At 24 months, all six mechanical parameters were reduced significantly in the control group relative to baseline values (2.9% to 5.1%, p = .03 to .002), at which point the reduction in BV/TV also became significant (1.8%, p < .05). None of the subregional between-group differences of changes from baseline were significant.

For the distal radius (Tables 3 and 4), the temporal changes in mean values of subregional mechanical parameters did not reach statistical significance except for E22 (at 12 months) and G23 (at 24 months) for the estradiol group, for which 6.1% ± 2.2% (p = .02) and 3.8% ± 1.5% (p = .04) increases were observed, respectively. Further, the detected change in BV/TV was nonsignificant for both groups at both time points. Again, at this site, there was no significant difference in the temporal changes between groups.

Full cross-sectional µFE analysis

Tables 1 through 4 also provide mean values and relative temporal changes in full cross-sectional axial stiffness at the two anatomic sites from baseline to the follow-up time points. Figure 3 displays the fractional changes from baseline to 12 and 24 months. Estrogen supplementation increased whole-section axial stiffness after 12 and 24 months of treatment at both anatomic sites. At the tibia, the increases were 5.8% and 6.3% (both p = .03) at 12 and 24 months, respectively. At the radius, axial stiffness increased 2.5% (p = .02) at 12 months and 3.2% at 24 months, but the latter effect did not quite reach significance (p = .065). In control subjects, the data further revealed substantial and highly significant reductions in axial stiffness (3.9% and 4.0%, p = .0004 and .001 at 12 and 24 months, respectively) at the tibia and 3.6% and 3.7% (p = .002 and .03) at the radius at the two time points.

The intergroup differences in the observed temporal changes in axial stiffness were highly significant at both follow-up time points and for both anatomic sites. However, the values were greater for the distal tibia than for the distal radius at 12 months (9.7% ± 2.5%, p = .0001, versus 6.1% ± 1.5%, p = .0003). At 24 months, the corresponding intergroup changes were 10.2% ± 2.3% (p = 0.0002) at the tibia versus 6.8% ± 2.1% (p = 0.004) at the radius. Further, the full cross-sectional axial stiffness values at tibial and radial sites were highly correlated with each other, indicating the systemic nature of these parameters (Fig. 4).

Relationship between mechanical and topologic parameters

Strong correlations were observed between the mechanical and composite topologic parameters (S/C and EI) computed from the same trabecular bone subvolumes at the distal tibia (Table 5). The strength of the correlations differed for the three directions of Young's moduli, being largest for the axial direction (E33). Further, the predictive ability of two-parameter linear regression models involving BV/TV and either S/C or EI improved relative to a single predictor. Similar trends were observed for the correlations between mechanical and composite topologic parameters computed from the distal radius. However, the correlations (r2 ranging from 0.24 to 0.89) generally were somewhat weaker than those observed at distal tibia.


The results from this work show that µMRI-based structural images from trabecular and cortical bone at the distal extremities acquired in patients are suited to derive finite-element models for estimating the mechanical consequences of estrogen loss following menopause and estrogen supplementation. This study is only one of a small number of in-vivo image-based studies compellingly showing not only the structural but also mechanical manifestations of intervention. A prior study involving a smaller number of subjects had provided evidence of the architectural(18) and, subsequently, mechanical effects of testosterone replacement in men with hypogonadism.(28) Unlike in that work, the effects expected here were smaller in consideration that the women studied were early postmenopausal (mean approximately 1 year from last menstrual period), in contrast to testosterone-deficient men who had been hypogonadal for very long times and therefore eliciting a very strong treatment response. Further, this work goes considerably beyond the prior pilot study in that we performed subregional µFE analysis of a small cuboid trabecular subregion (which, however, was larger than in ref 28)—10 × 10 × 5 mm3 versus 5 × 5 × 5 mm3—but also examined the role of the cortex by performing full cross-sectional stiffness analysis similar to that described by MacNeil and Boyd.(40) Further, we studied the hypothesized mechanical treatment effects at both the distal tibia and the radius, that is, a load-bearing and a non-load-bearing anatomic site. We further brought advanced registration methods to bear for mutual alignment of the images at the three time points and with the scanner coordinate frame.(33)

In the tibia, the subregion analysis revealed significant reductions in all six moduli in the control group (ie, women who elected not to take estrogen) ranging from 3% to 5% from baseline after 24 months (p < .05 to <.005), whereas only three of the moduli were significantly lower after 12 months. In contrast, in the estradiol group, all six mechanical constants remained unaltered at both follow-up time points. Between-group changes in the elastic moduli, in contrast to whole-section axial stiffness values, were not significant. These results are somewhat at variance with the structural changes reported for the same cohort of patients in ref. (20). This discrepancy is likely due to the constraints imposed by the mechanical analysis in the choice of the subvolume in the present work (10 × 10 × 5 mm3 parallelepiped at the centroid versus larger irregularly shaped region for prior structural analysis in ref 20). It is interesting to note, however, that despite the highly significant decrements in the mechanical constants in the control subjects, BV/TV decreased only marginally at 24 months (−1.8%, p = .03), whereas the reduction at 12 months did not reach significance (−1.0%, p = .07). The data thus suggest that the predicted changes in elastic and shear moduli are largely the result of structural changes, notably the observed reduction in plate-to-rod ratio observed previously in the data.(20) Similar topologic changes had been previously observed in paired biopsies following early menopause.(41)

Of particular relevance are the effects of intervention on whole-section axial stiffness at both anatomic sites. In the tibia, women in the control group suffered a loss in estimated stiffness of about 4% within a year (p < .005) while gaining over 6% in response to estrogen supplementation during the same period (p < .05). In the radius, whole-section axial stiffness declined by about 4% during the initial 12 months (p < .005) and increased marginally in the estradiol group (by about 2.5%, p < .05). Of significance further are the large intergroup differences in change from baseline: 10% after 12 months in the tibia (p = .0001) and 6% in the radius (p < .0005).

The effect of estradiol is clearly greater at the tibia than at the radius, suggesting that the load-bearing site may be more sensitive to antiresorptive treatment. Possible synergistic effects of antiresorptive agents with mechanical stimulation have been studied widely (see, for example, ref (42)). However, most work relies on studies in animals. For example, Westerlind and colleagues(43) found that estrogen administration to ovariectomized rats reduced bone loss in the unloaded and prevented loss in the loaded limb following unilateral sciatic neurotomy. Tromp and colleagues(44) also found an interaction between estrogen and mechanical loading on bone in the ovariectomized rat model.

The highly significant between-group changes from baseline in whole-section axial stiffness at both anatomic sites emphasize the role of the cortex as being affected by estrogen depletion and supplementation. The findings are counter to the notion that short-term remodeling changes after menopause (both loss owing to reduction in estrogen levels or accrual from estrogen supplementation) affect primarily trabecular bone, which remodels faster than cortical bone. Further, the within-group increases in trabecular bone elastic and shear moduli, as well as whole-section axial stiffness in the estradiol group suggest that estrogen supplementation not only preserves but actually improves bone mechanical competence, in contrast to prior structural studies indicating mere preservation of architecture.(45) Nevertheless, the anabolic effects of estrogen replacement are well documented. Khastgir and colleagues found in older women large increases in trabecular bone volume and bone mineral density (BMD). as well as increased wall thickness, in response to estrogen supplementation.(46) Similarly, based on paired bone biopsies obtained before and after 2 years of hormone replacement in postmenopausal women, Vedi and colleagues demonstrated significantly higher cortical width.(47)

The lack of further significant gains from 12 to 24 months is somewhat puzzling but parallels the structural observations made previously.(20) Nevertheless, in the tibia, four of the elastic moduli that had not changed significantly after 12 months were all lower after 24 months (Fig. 2). In contrast, there was no further decline in whole-section stiffness after 24 months and only an insignificant further increase in the estradiol group. There are several possible reasons for these observations. First, lifestyle changes in the control group could have mitigated further declines in structural and mechanical parameters. Second, in these very early postmenopausal women, the degradation in bone quality is likely to be largest early after the onset of menopause. Third, the number of subjects proceeding to the 24-month time point was lower than those who were examined at baseline.

This study has limitations because it ignores effects of treatment on the intrinsic bone tissue properties. Paschalis and colleagues(2) found in early postmenopausal women that 2 years of estrogen/progestin supplementation increased mineral/matrix ratio, mineral crystallinity/maturity, and the relative ratio of collagen cross-links, clearly showing that antiresorptive treatment affects material properties. The same group recently has provided evidence that subjects with vertebral fractures differ in their matrix and mineral composition (ie, collagen maturity, mineral/matrix ratio, and crystallinity).(48) However, the work did not provide data on structural group differences, and BMD was matched at the proximal femur rather than at the vertebrae.

As with all CT-based µFE work using in vivo structural images to create the FE model.(49,50) our model analyzes the structural implications only of longitudinal changes on mechanical parameters. While CT can, in principle, make use of variations in mineralization density, such differences are not distinguishable from partial-volume averaging in the limited spatial resolution regime of in vivo imaging, where resolution is comparable with trabecular thickness, nor would CT be able to quantify changes in matrix chemistry. Also, with few exceptions,(36) all prior in vivo µFE work was based on thresholded images in that a fixed tissue modulus was given to the voxels assigned to bone. We found in preliminary work grayscale-based µFE analyses to be more accurate than those relying on binarized images,(51) which is what prompted us to avoid binarization, since it is well known to be fraught with error at the resolution achievable by in vivo structural imaging of trabecular bone. Finally, the study is limited to the linear regime in that ultimate strength was not evaluated, which could be obtained from the present data by invoking the Pistoia criterion,(52) which suggests that failure will occur once 2% of the bone tissue is strained beyond 0.7%.(52) Nevertheless, there is substantial evidence that yield strength and stiffness are highly correlated (r2 > 0.95; see, for example, ref. (36)).

In conclusion, MRI-based µFE analysis is able to predict the short-term mechanical implications of drug intervention of both trabecular and cortical bone. The data provide evidence of an anabolic effect of estrogen supplementation in terms of whole-section axial stiffness at both the distal tibia and radius, but the effects, in terms of both reduction in mechanical competence in the absence of estradiol administration and the observed increases in subjects receiving hormone supplementation, were greater at the tibia than at the radius.

We are indebted to Dr Andrew Cucchiara for advice on the statistical analysis and acknowledge Mr J Love's assistance in processing the images. This work was supported by NIH Grants RO1 AR41443, AR53156, and AR55647, CTRC University of Pennsylvania School of Medicine, and Novartis Pharmaceuticals.


All authors state that they have no conflicts of interest.

1. Seeman E,Delmas PD. Bone quality–the material and structural basis of bone strength and fragilityN Engl J MedYear: 200635422506116723616
2. Paschalis EP,Boskey AL,Kassem M,Eriksen EF. Effect of hormone replacement therapy on bone quality in early postmenopausal womenJ Bone Miner ResYear: 200318955912817747
3. Recker RR,Barger-Lux MJ. The elusive concept of bone qualityCurr Osteoporos RepYear: 200429710016036089
4. Sarkar S,Mitlak BH,Wong M,Stock JL,Black DM,Harper KD. Relationships between bone mineral density and incident vertebral fracture risk with raloxifene therapyJ Bone Miner ResYear: 20021711011771654
5. Genant HK,Jiang Y. Advanced imaging assessment of bone qualityAnn N Y Acad SciYear: 200610684102816831940
6. Wehrli FW. Structural and functional assessment of trabecular and cortical bone by micro magnetic resonance imagingJ Magn Reson ImagingYear: 20072539040917260403
7. Bouxsein ML. Technology insight: noninvasive assessment of bone strength in osteoporosisNat Clin Pract RheumatolYear: 20084310818431371
8. Hwang SN,Wehrli FW,Williams JL. Probability-based structural parameters from 3D NMR images as predictors of trabecular bone strengthMedical PhysicsYear: 199724125512619284249
9. Ulrich D,van Rietbergen B,Laib A,Ruegsegger P. The ability of three-dimensional structural indices to reflect mechanical aspects of trabecular boneBoneYear: 199925556010423022
10. Genant HK,Engelke K,Prevrhal S. Advanced CT bone imaging in osteoporosisRheumatology (Oxford)Year: 200847Suppl 4iv91618556648
11. Techawiboonwong A,Song HK,Leonard MB,Wehrli FW. Cortical bone water: in vivo quantification with ultrashort echo-time MR imagingRadiologyYear: 20082488243318632530
12. Link TM,Majumdar S,Augat P,et al. In vivo high resolution MRI of the calcaneus: differences in trabecular structure in osteoporosis patientsJournal of Bone and Mineral ResearchYear: 1998131175829661082
13. Wehrli FW,Gomberg BR,Saha PK,Song HK,Hwang SN,Snyder PJ. Digital topological analysis of in vivo magnetic resonance microimages of trabecular bone reveals structural implications of osteoporosisJ Bone Miner ResYear: 20011615203111499875
14. Boutry N,Cortet B,Dubois P,Marchandise X,Cotten A. Trabecular bone structure of the calcaneus: preliminary in vivo MR imaging assessment in men with osteoporosisRadiologyYear: 20032277081712676974
15. Boutroy S,Bouxsein ML,Munoz F,Delmas PD. In Vivo Assessment of Trabecular Bone Microarchitecture by High-Resolution Peripheral Quantitative Computed TomographyJ Clin Endocrinol MetabYear: 200590680515
16. Khosla S,Riggs BL,Atkinson EJ,et al. Effects of sex and age on bone microstructure at the ultradistal radius: a population-based noninvasive in vivo assessmentJ Bone Miner ResYear: 2006211243116355281
17. Chesnut CH 3rd,Majumdar S,Newitt DC,et al. Effects of salmon calcitonin on trabecular microarchitecture as determined by magnetic resonance imaging: results from the QUEST studyJ Bone Miner ResYear: 20052015486116059627
18. Benito M,Vasilic B,Wehrli FW,et al. Effect of testosterone replacement on bone architecture in hypogonadal menJ Bone Miner ResYear: 20052017859116160736
19. Graeff C,Timm W,Nickelsen TN,et al. Monitoring Teriparatide associated Changes in Vertebral Microstructure by High-Resolution Computed Tomography in vivo: Results from the EUROFORS studyJ Bone Miner ResYear: 20072214263317547537
20. Wehrli FW,Ladinsky GA,Jones C,et al. In vivo magnetic resonance detects rapid remodeling changes in the topology of the trabecular bone network after menopause and the protective effect of estradiolJ Bone Miner ResYear: 2008237304018251704
21. Liu XS,Sajda P,Saha PK,Wehrli FW,Guo XE. Quantification of the roles of trabecular microarchitecture and trabecular type in determining the elastic modulus of human trabecular boneJ Bone Miner ResYear: 20062116081716995816
22. Ciarelli TE,Fyhrie DP,Schaffler MB,Goldstein SA. Variations in three-dimensional cancellous bone architecture of the proximal femur in female hip fractures and in controlsJ Bone Miner ResYear: 200015324010646112
23. Keaveny TM,Morgan EF,Niebur GL,Yeh OC. Biomechanics of trabecular boneAnnu Rev Biomed EngYear: 200133073311447066
24. Kabel J,van Rietbergen B,Dalstra M,Odgaard A,Huiskes R. The role of an effective isotropic tissue modulus in the elastic properties of cancellous boneJ BiomechYear: 1999326738010400354
25. Rajapakse CS,Magland J,Zhang XH,et al. Implications of noise and resolution on mechanical properties of trabecular bone estimated by image-based finite-element analysisJ Orthop Res.Year: 2009
26. Muller R,Ruegsegger P. Analysis of mechanical properties of cancellous bone under conditions of simulated bone atrophyJ BiomechYear: 1996291053608817372
27. van Rietbergen B,Majumdar S,Newitt D,MacDonald B. High-resolution MRI and micro-FE for the evaluation of changes in bone mechanical properties during longitudinal clinical trials: application to calcaneal bone in postmenopausal women after one year of idoxifene treatmentClin Biomech (Bristol, Avon)Year: 200217818
28. Zhang XH,Liu XS,Vasilic B,et al. In vivo microMRI-based finite element and morphological analyses of tibial trabecular bone in eugonadal and hypogonadal men before and after testosterone treatmentJ Bone Miner ResYear: 20082314263418410234
29. Song HK,Wehrli FW. In vivo micro-imaging using alternating navigator echoes with applications to cancellous bone structural analysisMagnetic Resonance in MedicineYear: 19994194795310332878
30. Ladinsky GA,Vasilic B,Popescu AM,et al. Trabecular Structure Quantified With the MRI-Based Virtual Bone Biopsy in Postmenopausal Women Contributes to Vertebral Deformity Burden Independent of Areal Vertebral BMDJ Bone Miner ResYear: 200823647417784842
31. Vasilic B,Wehrli FW. A novel local thresholding algorithm for trabecular bone volume fraction mapping in the limited spatial resolution regime of in-vivo MRIIEEE Trans Med ImagingYear: 2005241574158516353372
32. Rajapakse CS,Magland JF,Wald MJ,et al. Estimation of Relative Stiffness Contributions of Cortical and Trabecular Compartments by MRI-based Finite Element AnalysisJournal of Bone and Mineral ResearchYear: 200823S366S366
33. Magland JF,Jones CE,Leonard MB,Wehrli FW. Retrospective 3D registration of trabecular bone MR images for longitudinal studiesJ Magn Reson ImagingYear: 2009291182619097098
34. Van Rietbergen B,Odgaard A,Kabel J,Huiskes R. Direct mechanics assessment of elastic symmetries and properties of trabecular bone architectureJournal of BiomechanicsYear: 199629165316578945668
35. Hollister SJ,Brennan JM,Kikuchi N. A homogenization sampling procedure for calculating trabecular bone effective stiffness and tissue level stressJ BiomechYear: 199427433448188724
36. Macneil JA,Boyd SK. Bone strength at the distal radius can be estimated from high-resolution peripheral quantitative computed tomography and the finite element methodBoneYear: 20084212031318358799
37. Gomberg BG,Saha PK,Song HK,Hwang SN,Wehrli FW. Application of topological analysis to magnetic resonance images of human trabecular boneIEEE Transactions on Medical ImagingYear: 20001916617410875701
38. Magland JF,Wehrli FW. Trabecular bone structure analysis in the limited spatial resolution regime of in vivo MRIAcad RadiolYear: 20081514829319000865
39. Manzanera A,Bernard TM,Preteux P,Longuet B. N-dimensional skeletonization: a unified mathematical frameworkJournal of Electronic ImagingYear: 2002112537
40. MacNeil JA,Boyd SK. Load distribution and the predictive power of morphological indices in the distal radius and tibia by high resolution peripheral quantitative computed tomographyBoneYear: 2007411293717442649
41. Akhter MP,Lappe JM,Davies KM,Recker RR. Transmenopausal changes in the trabecular bone structureBoneYear: 200741111617499038
42. Jee WS,Tian XY. The benefit of combining non-mechanical agents with mechanical loading: a perspective based on the Utah Paradigm of Skeletal PhysiologyJ Musculoskelet Neuronal InteractYear: 20055110815951626
43. Westerlind KC,Wronski TJ,Ritman EL,et al. Estrogen regulates the rate of bone turnover but bone balance in ovariectomized rats is modulated by prevailing mechanical strainProc Natl Acad Sci U S AYear: 19979441992049108129
44. Tromp AM,Bravenboer N,Tanck E,et al. Additional weight bearing during exercise and estrogen in the rat: the effect on bone mass, turnover, and structureCalcif Tissue IntYear: 2006794041517160577
45. Vedi S,Croucher PI,Garrahan NJ,Compston JE. Effects of hormone replacement therapy on cancellous bone microstructure in postmenopausal womenBoneYear: 19961969728830991
46. Khastgir G,Studd J,Holland N,Alaghband-Zadeh J,Fox S,Chow J. Anabolic effect of estrogen replacement on bone in postmenopausal women with osteoporosis: histomorphometric evidence in a longitudinal studyJ Clin Endocrinol MetabYear: 2001862899511232014
47. Vedi S,Bell KL,Loveridge N,Garrahan N,Purdie DW,Compston JE. The effects of hormone replacement therapy on cortical bone in postmenopausal women. A histomorphometric studyBoneYear: 200333330413678773
48. Gourion-Arsiquaud S,Faibish D,Myers E,et al. Use of FTIR spectroscopic imaging to identify parameters associated with fragility fractureJ Bone Miner ResYear: 20092415657119419303
49. Boutroy S,Van Rietbergen B,Sornay-Rendu E,Munoz F,Bouxsein ML,Delmas PD. Finite element analysis based on in vivo HR-pQCT images of the distal radius is associated with wrist fracture in postmenopausal womenJ Bone Miner ResYear: 200823392917997712
50. Macneil JA,Boyd SK. Accuracy of high-resolution peripheral quantitative computed tomography for measurement of bone qualityMed Eng Phys.Year: 2007
51. Magland JF,Rajapakse CS,Wald MJ,et al. Grayscale MR image based finite element mechanical modeling of trabecular bone at in vivo resolution 30th Annual Meeting of the American Society for Bone and Mineral ResearchYear: 2008Montréal, Quebéc; CanadaASBMRS310
52. Pistoia W,van Rietbergen B,Lochmuller EM,Lill CA,Eckstein F,Ruegsegger P. Estimation of distal radius failure load with micro-finite element analysis models based on three-dimensional peripheral quantitative computed tomography imagesBoneYear: 200230842812052451


[Figure ID: fig01]
Fig. 1 

Image acquisition and processing flow to compute full cross-sectional axial stiffness and subregional mechanical and topologic parameters of trabecular bone from µMRI images of distal radius and distal tibia.

[Figure ID: fig02]
Fig. 2 

Within-group mean temporal changes (%) in subregional mechanical parameters and BV/TV of distal tibia from baseline to 12 and 24 months. The p values indicate the significance of the computed changes. Between-group changes from baseline were not significant.

[Figure ID: fig03]
Fig. 3 

Within-group mean temporal changes (%) in full cross-sectional axial stiffness of (A) distal tibia and (B) distal radius from baseline to the two follow-up time points. Error bars represent standard error of mean changes. Intergroup differences in the temporal changes were significant at both follow-up time points and for both anatomic sites (p < .005). For details, see text.

[Figure ID: fig04]
Fig. 4 

Comparison of full cross-sectional stiffness (MPa) of distal tibia and distal radius computed from µMR-based µFE analysis at baseline.

[TableWrap ID: tbl1] Table 1 

Young's Moduli (E11, E22, E33) and Shear Moduli (G12, G13, G23), Cross-Sectional Stiffness (All in MPa), and BV/TV (%) at the Distal Tibia Derived from µMRI Images in Control Subjects at Baseline and 12 and 24 Months, Including Relative Changes of Mean Values from Baseline to Both Follow-Up Time Points

Mean value Baseline to 12 mos. Baseline to 24 mos./

Parameter Baseline 12 mos. 24 mos. Relative change (%) p Relative change (%) p
E11 243 ± 75 228 ± 70 233 ± 51 −5.53 .006 −5.14 .005
E22 199 ± 38 202 ± 46 192 ± 39 1.44 .309 −3.71 .018
E33 827 ± 164 799 ± 151 805 ± 126 −3.01 .013 −3.55 .003
G12 155 ± 30 155 ± 31 151 ± 28 −0.21 .850 −2.93 .005
G13 181 ± 45 172 ± 42 175 ± 32 −4.60 .002 −4.41 .002
G23 103 ± 24 102 ± 25 101 ± 21 −1.51 .231 −3.64 .005
Stiffness 1222 ± 144 1184 ± 150 1156 ± 126 −3.89 .0004 −4.03 .0011
BV/TV 11.0 ± 1.2 10.9 ± 1.1 10.9 ± 1.0 −0.99 .071 −1.76 .026

[TableWrap ID: tbl2] Table 2 

Subregional Young's Moduli (E11, E22, E33), Shear Moduli (G12, G13, G23), Cross-Sectional Stiffness (All in MPa), and BV/TV (%) at the Distal Tibia, Derived from µMR Images in Estradiol Subjects at Baseline and 12 and 24 Months, Including Relative Changes of Mean Values from Baseline to Both Follow-Up Time Points

Mean value Baseline to 12 mos. Baseline to 24 mos.

Parameter Baseline 12 mos. 24 mos. Relative change (%) p Relative change (%) p
E11 230 ± 68 231 ± 60 241 ± 59 2.89 .49 3.79 .44
E22 193 ± 38 194 ± 36 194 ± 30 0.96 .66 −1.63 .60
E33 840 ± 149 838 ± 140 859 ± 128 0.27 .90 0.10 .97
G12 153 ± 26 153 ± 24 156 ± 19 0.45 .79 −0.58 .80
G13 179 ± 45 179 ± 39 184 ± 40 1.75 .55 1.64 .61
G23 97 ± 20 98 ± 18 99 ± 16 2.38 .32 0.37 .91
Stiffness 1158 ± 139 1218 ± 116 1177 ± 80 5.76 .03 6.22 .03
BV/TV 11.2 ± 1.1 11.1 ± 1.0 11.3 ± 1.0 −0.35 .68 −0.26 .81

[TableWrap ID: tbl3] Table 3 

Subregional Young's Moduli (E11, E22, E33) and Shear Moduli (G12, G13, G23), Cross-Sectional Stiffness (All in MPa) and BV/TV (%) at the Distal Radius Derived from µMR Images in Control Subjects at Baseline and 12 and 24 Months, Including Relative Changes of Mean Values from Baseline to Both Follow-Up Time Points

Mean value Baseline to 12 mos. Baseline to 24 mos.

Parameter Baseline 12 mos. 24 mos. Relative change (%) p Relative change (%) p
E11 302 ± 84 290 ± 67 286 ± 62 2.89 .49 3.79 .44
E22 404 ± 99 391 ± 82 373 ± 73 0.96 .66 −1.63 .60
E33 838 ± 117 813 ± 102 824 ± 127 0.27 .90 0.10 .97
G12 249 ± 47 242 ± 38 235 ± 37 0.45 .79 −0.58 .80
G13 211 ± 40 203 ± 33 201 ± 31 1.75 .55 1.64 .65
G23 170 ± 49 161 ± 40 152 ± 35 2.38 .32 0.37 .91
Stiffness 1547 ± 206 1542 ± 207 1487 ± 188 −3.57 .002 −3.65 .03
BV/TV 9.8 ± 1.2 9.8 ± 1.0 10.2 ± 1.7 −0.35 .68 −0.26 .81

[TableWrap ID: tbl4] Table 4 

Subregional Young's Moduli (E11, E22, E33) and Shear Moduli (G12, G13, G23), Cross-Sectional Stiffness (All in MPa), and BV/TV (%) at the Distal Radius Derived from µMR Images in Estradiol Subjects at Baseline and 12 and 24 months, Including Relative Changes of Mean Values from Baseline to Both Follow-Up Time Points

Mean value Baseline to 12 mos. Baseline to 24 mos.

Parameter Baseline 12 mos. 24 mos. Relative change (%) p Relative change (%) p
E11 281 ± 82 310 ± 111 327 ± 132 2.89 .490 3.79 .442
E22 376 ± 94 407 ± 114 409 ± 120 0.96 .655 −1.63 .596
E33 813 ± 88 799 ± 108 823 ± 133 0.27 .901 0.10 .966
G12 230 ± 42 239 ± 53 245 ± 60 0.45 .792 −0.58 .795
G13 197 ± 39 207 ± 53 212 ± 62 1.75 .554 1.64 .647
G23 157 ± 50 172 ± 65 177 ± 74 2.38 .322 0.37 .907
Stiffness 1382 ± 384 1477 ±163 1486 ± 145 2.52 .024 3.16 .065
BV/TV 9.4 ± 0.9 10.7 ± 3.6 12.7 ± 5.6 −0.35 .677 −0.26 .811

[TableWrap ID: tbl5] Table 5 

R2 Values for the Correlation Between Mechanical and Structural Parameters Computed from the Same 10 × 10 × 5 mm3 Trabecular-Bone Subregion of the Distal Tibia at Baseline for All Subjects

E11 E22 E33 G23 G31 G12
BV/TV 0.71 0.44 0.92 0.69 0.81 0.66
S/C 0.76 0.45 0.85 0.68 0.82 0.62
EI 0.64 0.35 0.74 0.54 0.72 0.50
BV/TV and S/C 0.80 0.47 0.96 0.74 0.88 0.69
BV/TV and EI 0.77 0.44 0.95 0.71 0.87 0.66

All correlations were significant (p < .0001).

Article Categories:
  • Original Article

Keywords: MRI, trabecular bone, cortical bone, structure, finite element, menopause, estrogen.

Previous Document:  Osteocalcin gene polymorphisms influence concentration of serum osteocalcin and enhance fracture ide...
Next Document:  Enhanced differentiation of human embryonic stem cells to mesenchymal progenitors by inhibition of T...