PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Magn Reson Imaging. Author manuscript; available in PMC Aug 1, 2012.
Published in final edited form as:
PMCID: PMC3144507
NIHMSID: NIHMS284312
MRI of trabecular bone using a DDIF contrast imaging sequence
Dionyssios Mintzopoulos,1 Jerome L. Ackerman,1 and Yi-Qiao Song1,2
1 Athinoula A Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129, USA
2 Schlumberger-Doll Research, 1 Hampshire Street, Cambridge, Massachusetts 02139, USA
Address correspondence to: Dionyssios Mintzopoulos, PhD, Athinoula A Martinos Center for Biomedical Imaging, MGH Building 149, 13th Street/5th Ave, Charlestown, MA 02129, USA, Tel: (617) 785-4566, dionyssi/at/nmr.mgh.harvard.edu
Purpose
To characterize the DDIF (Decay due to Diffusion in the Internal Field) method using intact animal trabecular bone specimens of varying trabecular structure and porosity, under ex vivo conditions closely resembling in vivo physiological conditions. The DDIF method provides a diffusion contrast which is related to the surface-to-volume ratio of the porous structure of bones. DDIF has previously been used successfully to study marrow-free trabecular bone, but the DDIF contrast hitherto had not been tested in intact specimens containing marrow and surrounded by soft tissue.
Materials and Methods
DDIF imaging was implemented on a 4.7 T small-bore, horizontal, animal scanner. Ex vivo results on fresh bone specimens containing marrow were obtained at body temperature. Control measurements were carried out in surrounding tissue and saline.
Results
Significant DDIF effect was observed for trabecular bone samples, while it was considerably smaller for soft tissue outside the bone and for lipids. Additionally, significant differences were observed between specimens of different trabecular structure.
Conclusion
The DDIF contrast is feasible despite the reduction of the diffusion constant and of T1 in such conditions, increasing our confidence that DDIF imaging in vivo may be clinically viable for bone characterization.
Keywords: trabecular bone, imaging, DDIF, diffusion, MRI contrast
Fragility fracture is a serious and costly public health issue. With the growth of the elderly segment of the United States (and global) population, there is increasing incidence of metabolic bone diseases. In the U.S. alone, approximately 1.5 million fractures occur annually while the incidence of osteoporosis (low bone mass) is estimated at 14 million patients and the incidence of osteopenia (reduced bone mass) at 30 million, with a concomitant annual projected economic impact of $20 billion (1,2). One-third of patients admitted to U.S. hospitals for hip (femoral neck) fracture die within one year, primarily due to the rapid decline in mobility and general quality of life (3). Improved methods for screening, diagnosis and treatment monitoring of metabolic bone disease are of great importance.
Bone strength is a key predictor of fracture risk (4). Direct mechanical testing of bone strength is destructive and invasive. In clinical in vivo applications, bone strength cannot be measured directly. Instead, it is estimated by the use of biomarkers. The widely accepted clinical biomarker standard for bone strength is bone mineral density (BMD), estimated by Dual-Energy X-ray Absorptiometry (DXA). In fact, the World Health Organization defines osteoporosis in terms of the DXA BMD score (5) which is used in fracture risk assessment tools such as FRAX (6). DXA-derived scores are predictive risk factors for fractures (7). DXA measurements provide two-dimensional measures of bone density (the signal is integrated along the third dimension, depth) and, for that reason, cannot distinguish trabecular architecture which is important to the overall biomechanical function of bone (8,9) and in bone disease (e.g. osteoporosis) (10). Currently, much is still not understood about bone strength and, therefore, noninvasive methods for characterizing bone mechanical strength are needed (4,11). Both trabecular and cortical bones contribute significantly to the overall ability of bones to withstand loading (12,13). Trabecular bone has a higher surface-to-volume ratio than cortical bone, which translates in higher metabolic activity and more rapid turnover (14,15). The classic technique to assess trabecular microarchitecture in patients is histomorphometry of biopsy specimens by light microscopy, most typically of the trabecular bone of the iliac crest (16), while μCT is rapidly becoming the current standard. However, biopsy is generally a method of last resort for both patients and physicians because it is a minor surgical procedure. Furthermore, biopsy may not be statistically representative of the skeleton as a whole, and cannot be performed in load-bearing bones (which are the bones of most concern in fracture risk).
The porous geometry of trabecular bone makes it amenable to study by proton magnetic resonance methods such as spectroscopy (NMR, MRS) and imaging (MRI). Diffusible water, which is present in the interstitial (pore) space of trabecular bone, can be studied using magnetic resonance methods to probe the pore space of trabecular structure. As a result, there is strong interest in applying non-invasive, clinically applicable magnetic resonance methods for the in vivo study of trabecular bone structure (17,18).
Other noninvasive imaging modalities, based on high-resolution MRI (1825), X-ray Computed Tomography (CT) and Ultrasound (US) (26,27) seek to improve fracture risk assessment by characterizing the microarchitecture of trabecular bone and estimating bone mechanical properties by post-processing of digitized images (2830). Considerable effort has been expended in the development of high-resolution micro-MRI to spatially resolve trabeculae (18,3135), as well as in the development of software to derive accurate imaging-based metrics that potentially correlate with biomechanical properties of healthy, aging, or diseased bone. These methods could yield statistical parameters to describe the microarchitecture, building on classic histomorphometry results (19,21,25,28,3645). Furthermore, solid state MRI is being developed to measure solid bone composition, another factor in bone strength and metabolic state (46). Additionally, much work has been done on the measurement of relaxation times, more specifically the reversible component equation M1 of the transverse relaxation time equation M2 (4749). Relaxation times may be influenced by static microscopic magnetic susceptibility distributions (resulting from the bone structure and geometry), but also by a number of other factors such as molecular diffusion through the resulting magnetic field gradients chemical exchange, and the presence of paramagnetic ions. Therefore, relaxation times may in some cases indirectly reflect biomechanical properties such as Young’s modulus (50); however, voxel-wise values of equation M3 (or equation M4) are dependent on the voxel size and they reflect the average orientation of trabecular elements relative to B0 (51). Lastly, other MR methods of trabecular bone characterization include manipulation of multiple quantum coherences to modulate the dipolar field, but the interpretation of these results is difficult (5254).
Here we focus on an alternative method, DDIF (Decay due to Diffusion in the Internal Field). The DDIF contrast is generated from differences in magnetic susceptibility in the tissue. Susceptibility imaging in itself is an active area of research (18,48,50,52,55), but the DDIF contrast is not the same as susceptibility imaging. DDIF was originally used in studies of porous media. It was developed as a spectroscopic tool for characterizing porous media such as the rock of a petroleum reservoir, by harnessing the spin dephasing which occurs when molecules diffuse through magnetic field gradients arising from spatial variations in magnetic susceptibility among the various structural phases of the material (5658). In porous media the susceptibility differences at the tissue-matrix interface result in internal field gradients that influence the proton signal due to diffusion in the porous space (59). The DDIF method is based on a stimulated echo sequence (56). A particle (for example a proton) diffusing in a porous medium from location x1 to x2, in the absence of external gradients accumulates a phase factor increment Δ[var phi] given from the relation
equation M5
[1]
where the symbol equation M6 denotes the local internal magnetic field at location x and and time t, and γ is the gyromagnetic ratio of the proton. The echo signal is proportional to the integral of that phase factor over a measure which is the diffusion propagator in pore space. One may gain a qualitative understanding of the DDIF signal in terms of diffusion propagator eigenfunctions (57,60). In order to gain an intuitive understanding, it helps to visualize the case of a simple two-component material, composed of a water-like fluid and of a solid matrix. The internal field induced within the fluid in the pore space has a complex non-uniform spatial structure. Its gradient is often largest near the tissue interface, i.e. near the solid matrix surface. As a proton diffuses, performing a random walk through the pore space, it accumulates phase in proportion to the local instantaneous internal field (Equation 1). The total phase accumulated depends on the surface-to-volume ratio of the porous material which, in turn, influences Bi. The MR signal is an integral of the phase factor over time, resulting in signal loss (DDIF weighting) as a function of diffusion time. When the diffusion or mixing time t is sufficiently long, each proton effectively samples all paths and the magnetization decays with a T1-like (slower) decay. In the context of a biomedical application (bone imaging) the DDIF contrast is the result of susceptibility-induced diffusion decay from the porous structure of the trabecular bone. The DDIF contrast does not require high-resolution MRI in order to resolve the trabeculae, and it provides a diffusion contrast that is related to the geometrical structure of the trabecular bone.
The DDIF method has been previously applied ex vivo to bovine trabecular specimens whose marrow had been removed and were immersed in saline (61,62). These ex vivo studies demonstrated that DDIF correlated with the mechanical compressive strength of ex vivo bovine bone specimens, that the measured DDIF decay constants correlated with simulated DDIF decay constants based on the trabecular bone surface geometry, and that DDIF measurements agreed well with the susceptibility differences induced by trabecular pores (approximately up to 200 μm). These ex vivo specimens approximate well the two-component idealization, but real bone is not a simple two-component material. Bone marrow is a complex material which to first approximation can be viewed as a two-component liquid itself, comprised of slowly diffusing large lipid molecules (adipocytic, yellow, marrow) and fast diffusing watery haematopoietic red marrow containing iron-carrying hemoglobin. These components differ from each other both in diffusion and in relaxation time constants (6366). In both the watery and lipid components, water is by far the most diffusible molecule, and is the dominant source of DDIF contrast. However, the reduction of water diffusion constant and of T1 in marrow compared to bulk water may reduce the DDIF contrast and make it harder to measure experimentally. The goal of this study is to examine whether DDIF is feasible in realistic bone specimens.
This study had two major aims. First, we implemented and demonstrated a DDIF imaging sequence with reduced sensitivity to imaging gradients. Second, we examined the performance of DDIF on fresh bone specimens containing marrow and studied at physiological (body) temperature. While this is an ex vivo study, we aimed to realistically approximate the conditions of an in vivo DDIF human study in a clinical context.
Animal Bone Samples
Fresh animal bone specimens were acquired from a local meatpacker (Superior Farms Boston) and from local supermarkets (Shaw’s Supermarkets Inc.). The bone specimens had been refrigerated to approximately 4°C prior to purchase. Upon purchase, the bone specimens were cut into smaller pieces able to fit in a standard 50 mL Falcon tube (Becton, Dickinson, Franklin Lakes, NJ, USA), using a hand saw (hacksaw) in order to minimize local damage from the heat generated from friction. Subsequently, the specimens were stored at −15°C. Seven specimens were examined for the purposes of this paper: four BV specimens, two PV, and one BR.
Some soft tissue was intentionally left in the specimens to enable comparison of the DDIF response between soft tissue and bone. Muscle tissue was carefully excised from the specimens. Prior to scanning, each specimen was placed in a 50-mL centrifuge tube and the remaining space in the tube was filled with normal saline (0.9% NaCl). The tube was then wrapped in a continuously heated water blanket and maintained at 34°C. Temperature equilibration was achieved in typically 40 min. Following temperature equilibration the specimen was placed in the coil and the coil was inserted into the magnet with the heating blanket still wrapped on one side, in order to maintain the temperature during the experiment.
Imaging
Imaging was performed in an Oxford Instruments (Oxford, UK) 4.7 T 33 cm horizontal bore magnet equipped with a Bruker BioSpin (Karlsruhe, Germany) Avance console, a Bruker gradient system capable of 40 G/cm, and a Bruker volume coil of 7 cm inner diameter and 10 cm active length. Shimming was performed using two repetitions of the automated shimming procedure, followed by manual shimming when necessary (typical water linewidth, ~15–20 Hz).
The DDIF imaging pulse sequence is schematically shown in Figure 1. It produces an image of the stimulated echo. All other possible signals are dephased. The algebraic conditions for coherence pathway selection in the stimulated echo using gradient lobes are well known (for example, chapter 10 of (67)). The gradient lobes between the second and third pulses serve to dephase unwanted coherences (gradient lobes labeled “2”). The slice-selecting gradient (gradient lobe “3”) is placed along the third pulse, as well. An additional gradient (“1”) follows after slice selection and before phase encoding, and it dephases the FID produced by the third pulse. In order to balance the coherence pathway for the stimulated echo, the same gradient is repeated (“1a”) between the first and second RF pulses. This is an external gradient that biases the DDIF by inducing a gradient-dependent decay rate, and so it must be minimized. In general, we can write the following equation for the measured DDIF decay rate:
equation M7
[2]
where equation M8 is the DDIF decay rate due to the internal field and equation M9 is the extra decay rate induced by the sequence (i.e. by externally applied gradients), defined from the expression
equation M10
[3]
where γ= 2π × 42.58MHz/T is the gyromagnetic ratio of the proton. The gradient was minimized empirically, finding the smallest possible balanced gradients (labeled “1” and “1a” in Figure 1) such that the FID from the third pulse was dephased and the stimulated echo image was formed. That optimization was carried out on a saline phantom. The gradient vector (1.8 G/cm, 1.8 G/cm, 1.3 G/cm) was played for one ms, which corresponds to a bias ( equation M11) for water (D ≈ 2.5 × 10−5 cm2/s). A minimal phase-cycling scheme was used. Two images were averaged with relative phase π between them (ie, [0 0 0], [π π π]). The DDIF encoding time is the period TE (echo time). The DDIF mixing, or diffusion, time is the period TM. For DDIF imaging, the echo time TE is set to a fixed value and a collection of images is acquired with variable TM. Typically, in our experiment TE = 10 ms and seven to nine TM values ranging from 15 ms to 800 ms were acquired at fixed TE. The overall repetition time TR was kept constant (2000 ms). Frequency-selective imaging was achieved by making the 2nd RF pulse frequency-selective (ten-lobed sinc pulse, 600 Hz pulse bandwidth, 10.35 ms). For lipid excitation, the carrier frequency of the sinc pulse was placed at −660 Hz (−3.3 ppm) relative to the water signal. Typical image acquisition times for a 64 × 64 matrix were approximately 4 minutes for the minimum number of averages (Navg = 2 according to the phase-cycling scheme), and field of view 40 mm (0.625 mm × 0.625 mm × 1.5 mm pixel dimensions).
Image Analysis
DDIF images were acquired for fixed encoding time TE and a set of variable mixing times TM. The set of all these images was used to produce a DDIF decay curve. Regions of interest (ROI) were drawn by hand on the images to delineate areas containing trabecular bone, saline in the 50 mL Falcon tube, and fat or other remaining tissue outside the cortical bone (MRIcroN, http://www.cabiatl.com/mricro/). Subsequent analysis was carried out using custom-designed code written in MATLAB. For each of the images in a given DDIF series, image intensities were estimated voxel-wise for each of the various ROI. The ROI-averaged intensities were used to produce the DDIF decay curve (for each specimen and for each ROI) which was used to estimate ROI-wise decay constants.
Estimation of T1 and DDIF decay constants
T1 was estimated using a multi-point inversion-recovery spin-echo (IRSE) pulse sequence. The repetition time (TR) was kept fixed for all images, while the inversion time, TI, was varied. T1 was estimated from parametric fits to the image-intensity curves as a function of TI. To account for the fact that TR may not be regarded as infinite relative to T1 and TI, the IRSE data (the ROI-wise averaged signal intensities S) were fitted to the following parametric equation for three unknown parameters, the background constant c, intensity a, and T1 relaxation.
equation M12
[4]
Equation (4) is valid in the limit T E[double less-than sign] TI, TR (for example see Ref (68)). The DDIF decay curves were fitted to a monoexponential decay plus constant noise (three-parameter model)
equation M13
[5]
Fitting was performed in MATLAB using lsqnonlin for parameter estimation. The MATLAB algorithm nlparci was used to produce confidence intervals from the parameter estimates, the residuals, and the Jacobian matrix.
Statistical Analysis
Between-groups comparison was performed using repeated-measures ANOVA with Bonferroni correction for multiple comparisons. Criterion for significance was P = 0.05 corrected. Analysis was performed using SPSS (SPSS v12, SPSS Inc., Chicago, IL). Effect sizes in Table 3 were calculated as Cohen’s d.
Table 3
Table 3
Statistical comparisons among trabecular DDIF time constants.
Computation of weighted averages
Weighted mean and weighted variance were calculated following standard methods (69). Given a set of N measurements xi with associated uncertainties σxi, the weighted mean is given by
equation M14
[6]
The weights are calculated from the uncertainties,
equation M15
[7]
The pooled (weighted) variance of the mean is given by
equation M16
[8]
Testing of the fitting algorithm
The fitting functions were tested on synthetic a priori data that were created with the three-component formula
equation M17
[9]
where η ~ N(0,σ2). The variance σ2 was calibrated so as to realize a priori (synthetic) χ2 similar to the experimental χ2. In order to mimic the experimental number of degrees of freedom, the time points t for the synthetic data were in the same range and had the same spacing as the MR time variable (mixing time TM), typically, 6–8 timepoints from TM = 10 ms to TM = 800 or 900 ms. The synthetic data were fitted with four models (Table 1): (a) a monoexponential decay (two-parameter model: one amplitude and one decay constant); (b) a monoexponential decay plus constant noise (three-parameter model); (c) a biexponential decay (four-parameter model: two amplitudes and two decay constant); and (d) a triexponential decay (six-parameter model: three amplitudes and three decay constants. A typical result of these numerical experiments (for a particular realization of the vector η) is shown in Table 1. Fitting was performed in MATLAB using lsqnonlin for parameter estimation. The MATLAB algorithm nlparci was used to produce confidence intervals from the parameter estimates, the residuals, and the Jacobian matrix). In general, it is well known that the inverse problem of how to estimate parameters by fitting sums of exponentials to a dataset (here, fitting to the DDIF decay curves) is notoriously ill-conditioned (for example see Ref (70)). The results in Table 1 quantify this for typical fits to synthetic data that approximate the experimental DDIF decay curves. However, the true, experimental, data follow a multi-exponential decay with an unknown number of terms. Therefore, experimental data were fitted with model (b).
Table 1
Table 1
Fitting algorithm: Typical performance on synthetic data
Effect Size
Effect size was measured with Cohen’s d (71).
Typical DDIF images are shown in Figures 2 and and3.3. Figure 2 shows DDIF images of a phantom and of a bovine vertebral specimen. The selective-excitation experiment (Figure 2, Left) was carried out on a phantom composed of two tubes attached to each other, a 3 mL syringe tube filled with saline and a second, narrower, tube filled with olive oil. These images were acquired with three different types of RF excitation: with hard RF pulses; with water-selective RF excitation; and with lipid-selective RF excitation. The residual image, defined as the difference between the hard-pulse image and the sum of the water- and lipid-selective images, (oil/water phantom, Figure 2d, Left; bovine specimen, Figure 2d, Right) is within the noise level in each case, verifying that the two selective-excitation images sum to the hard-pulse excitation image. The trabecular area of the bovine vertebral (BV) specimens had very little lipid content. The image using hard-pulse excitation (Figure 2a, Right) is similar to that using the water-selective excitation (Figure 2b, Right); and the DDIF decay curves from the corresponding experiments are also similar (Figure 4 and Figure 6). The situation was different for the porcine vertebral (PV) specimen (Figure 3), where the trabecular ROI exhibits significant lipid content (Figure 3c). The bovine rib (BR) specimen also had a significant lipid marrow contribution (image not shown).
Typical DDIF decay curves from the trabecular ROI are shown in Figure 4. Figure 4(a) shows decays from two bovine vertebral (BV) specimens. Both BV specimens were measured with both the hard-pulse and the water-selective sequences. These two specimens exhibit low lipid content. A replicate water-selective DDIF set of data was acquired three days later on the second BV specimen. The result was consistent with the first experiment, demonstrating the repeatability of DDIF imaging. Figure 4(b) shows all the water (red marrow) trabecular decay curves from all the specimens measured for this work (BV, BR, PV). The bovine specimens fall in a range with the rib specimen being on one end of the range (the lowest curve of the group) and the vertebral specimens having slightly slower decays and being above it. The porcine vertebral specimens resulted in significantly different curves compared to the bovine specimens.
Typical experimental data from the multi-point inversion-recovery spin-echo pulse sequence used to determine T1 are shown in Figure 5. The experimental data are from specimens BV4 and PV2. They show the inversion recovery curves for the trabecular, for the saline, and for the lipid ROI, together with their parametric fits. Parametric fits were estimated as previously outlined in the Methods section.
DDIF decay curves of a typical bovine (specimen BV4) and of a typical porcine specimen (specimen PV1) are plotted in Figures 6 and and7,7, respectively. The best-fitting lines through the data are parametric fits to the monoexponential three-parameter model, as outlined in the Methods section. The hard-pulse (non-selective, NS) and water-selective (S) excitations result in identical decay curves. The saline DDIF curve is clearly very different from the trabecular DDIF curve. Similar data are shown in Figure 7 for porcine specimen PV1. Specifically, shown in Figure 7 are, the water-selective trabecular DDIF decay of PV1 (faster decaying curve) as well as the lipid-selective DDIF curves of both the fatty marrow inside the trabecular ROI, and of the lipid tissue outside the bone. The (intensity-normalized) lipid DDIF curves are identical to each other, demonstrating that there is no DDIF effect for the lipids.
The performance of the algorithms used to solve the inverse problem, and to estimate DDIF decay constants from data, is illustrated in Table 1. Table 1 lists typical numerical results on synthetic a priori data that realistically mimic the experimental DDIF data. Two standard algorithms, functions nls (R) and lsqnonlin (MATLAB) were used. Both functions implement a nonlinear least squares optimization employing the Levenberg-Marquardt algorithm (72,73). The bias in the parameter estimation is inherent to the procedure, as can be seen from the right-most columns of Table 1 (fit to data on noiseless data).
Estimated DDIF decay constants are summarized in Table 2. Data were fitted with a three-parameter monoexponential fit as described in Methods. and estimated T1 relaxation constants, arranged for tissue and ROI. The saline T1 estimate from the PV specimens is very similar to that of the BV specimens, as expected. Statistical comparisons among trabecular DDIF constants across specimens of different trabecular structure (i.e. over different specimens BV, PV, and BR) are shown in Table 3. Effect size was computed as Cohen’s d.
Table 2
Table 2
Summary values of DDIF decay time constant and of estimated T1 values.
In summary, significant DDIF effect was observed for trabecular bone samples. The DDIF effect was considerably smaller for soft tissue outside the bone, and for lipids. Further, significant differences were observed between specimens of different trabecular structure.
Previous work on trabecular bone specimens whose marrow was removed (62) has shown a significant correlation between DDIF measurement and metrics of structural parameters such as the surface-to-volume ratio. However, such clean bone specimens are sufficiently different from in vivo bone tissue in many ways. First, the T1 values of marrow water and fat can be significantly shorter than those of pure water. Second, the water diffusion constant in marrow is lower than the self diffusion constant in pure water. Third, diffusion in tissue is further reduced due to the presence of cells and other structures. Both a shorter T1 and a reduced diffusion may reduce the DDIF contrast in vivo. Thus, it is an important step to examine the bone specimens with marrow in order to ascertain the feasibility of DDIF measurements for future in vivo applications.
Clear differences between the DDIF signal in the trabecular areas of the bone specimens and the DDIF signal in saline, were indeed observed (Figure 6). Generally, DDIF decay constants in saline are comparable to the T1 relaxation constants, whereas the water signal in the trabecular ROI decays faster, demonstrating trabecular DDIF contrast, in agreement with theoretical considerations about the DDIF (56) and with prior literature on DDIF measurements in porous media and in cleaned bone (56,58,62). In contrast, there was no difference between the DDIF decay in fatty marrow within trabecular bone and the DDIF decay in fat outside the bone (Figure 7, lipid images acquired with lipid-selective RF pulses). This is consistent with the fact that fat molecules are much larger than water and their diffusion constant is significantly smaller, and thus the DDIF effect is much smaller, too. In summary, these results demonstrate the observation of the DDIF contrast for water, and the lack thereof for lipids, in agreement with previous literature on the diffusion of water-like red marrow and of fatty marrow (lipids diffusing much more slowly than water) (6366).
In this work, specimens of different trabecular structure (bovine vs porcine) were imaged in order to test the feasibility of DDIF on different bones. Indeed, both the DDIF decay time constant and T1 are significantly different between PV and BV and the effect size was correspondingly large (see Table 3). Typically, effect sizes of 0.2–0.3 are considered to be small; 0.5, medium; and above 0.8–1.0 large (71). Effect sizes in Table 3 are clearly large, ranging from 1.4 to 5.4.
In conclusion, an imaging DDIF sequence was applied in an ex vivo study of fresh trabecular specimens. Significantly different DDIF signals among trabecular bone marrow, soft tissue outside the bone, and water were observed. Additionally, significant differences between specimens of different trabecular structure (bovine and porcine trabecular vertebral specimens) were observed. Together, these results indicate that DDIF imaging is possible in spite of the reduction of T1 and diffusion coefficient in bone marrow, and suggest that the application of DDIF in vivo is possible for improving bone characterization.
Acknowledgments
This work was supported in part by Schlumberger-Doll Research, Cambridge, MA, National Institutes of Health grants S10-RR16811 and P41-RR14075 (National Institute of Research Resources), the Athinoula A. Martinos Center for Biomedical Imaging, and the MIND Institute.
1. Rousculp MD, Long SR, Wang S, Schoenfeld MJ, Meadows ES. Economic burden of osteoporosis-related fractures in Medicaid. Value Health. 2007;10(2):144–152. [PubMed]
2. Becker DJ, Kilgore ML, Morrisey MA. The societal burden of osteoporosis. Curr Rheumatol Rep. 2010;12(3):186–191. [PubMed]
3. Roche JJ, Wenn RT, Sahota O, Moran CG. Effect of comorbidities and postoperative complications on mortality after hip fracture in elderly people: prospective observational cohort study. BMJ. 2005;331(7529):1374. [PMC free article] [PubMed]
4. Bouxsein ML, Seeman E. Quantifying the material and structural determinants of bone strength. Best Pract Res Clin Rheumatol. 2009;23(6):741–753. [PubMed]
5. Genant HK, Cooper C, Poor G, et al. Interim report and recommendations of the World Health Organization Task-Force for Osteoporosis. Osteoporos Int. 1999;10(4):259–264. [PubMed]
6. Blake GM, Fogelman I. An update on dual-energy x-ray absorptiometry. Semin Nucl Med. 2010;40(1):62–73. [PubMed]
7. Johnell O, Kanis JA, Oden A, et al. Predictive value of BMD for hip and other fractures. J Bone Miner Res. 2005;20(7):1185–1194. [PubMed]
8. Parfitt AM, Mathews CH, Villanueva AR, Kleerekoper M, Frame B, Rao DS. Relationships between surface, volume, and thickness of iliac trabecular bone in aging and in osteoporosis. Implications for the microanatomic and cellular mechanisms of bone loss. J Clin Invest. 1983;72(4):1396–1409. [PMC free article] [PubMed]
9. Kleerekoper M, Villanueva AR, Stanciu J, Rao DS, Parfitt AM. The role of three-dimensional trabecular microstructure in the pathogenesis of vertebral compression fractures. Calcif Tissue Int. 1985;37(6):594–597. [PubMed]
10. Aaron JE, Johnson DR, Paxton S, Kanis JA. Secondary osteoporosis and the microanatomy of trabecular bone. Clin Rheumatol. 1989;8 (Suppl 2):84–88. [PubMed]
11. Melton LJ, 3rd, Riggs BL, Keaveny TM, et al. Structural determinants of vertebral fracture risk. J Bone Miner Res. 2007;22(12):1885–1892. [PubMed]
12. Oden ZM, Selvitelli DM, Hayes WC, Myers ER. The effect of trabecular structure on DXA-based predictions of bovine bone failure. Calcif Tissue Int. 1998;63(1):67–73. [PubMed]
13. Currey J. Bones: Structure and Mechanics. Princeton University Press; 2006.
14. Hernandez CJ, Gupta A, Keaveny TM. A biomechanical analysis of the effects of resorption cavities on cancellous bone strength. J Bone Miner Res. 2006;21(8):1248–1255. [PMC free article] [PubMed]
15. Hernandez CJ. How can bone turnover modify bone strength independent of bone mass? Bone. 2008;42(6):1014–1020. [PMC free article] [PubMed]
16. Kumarasinghe DD, Perilli E, Tsangari H, et al. Critical molecular regulators, histomorphometric indices and their correlations in the trabecular bone in primary hip osteoarthritis. Osteoarthritis Cartilage. 2010 [PubMed]
17. Carballido-Gamio J, Majumdar S. Clinical utility of microarchitecture measurements of trabecular bone. Curr Osteoporos Rep. 2006;4(2):64–70. [PubMed]
18. Wehrli FW. Structural and functional assessment of trabecular and cortical bone by micro magnetic resonance imaging. J Magn Reson Imaging. 2007;25(2):390–409. [PubMed]
19. Majumdar S, Kothari M, Augat P, et al. High-resolution magnetic resonance imaging: three-dimensional trabecular bone architecture and biomechanical properties. Bone. 1998;22(5):445–454. [PubMed]
20. Genant HK, Gordon C, Jiang Y, et al. Advanced imaging of the macrostructure and microstructure of bone. Horm Res. 2000;54 (Suppl 1):24–30. [PubMed]
21. Wehrli FW, Saha PK, Gomberg BR, et al. Role of magnetic resonance for assessing structure and function of trabecular bone. Top Magn Reson Imaging. 2002;13(5):335–355. [PubMed]
22. Link TM, Majumdar S. Current diagnostic techniques in the evaluation of bone architecture. Curr Osteoporos Rep. 2004;2(2):47–52. [PubMed]
23. Krug R, Banerjee S, Han ET, Newitt DC, Link TM, Majumdar S. Feasibility of in vivo structural analysis of high-resolution magnetic resonance images of the proximal femur. Osteoporos Int. 2005;16(11):1307–1314. [PubMed]
24. Wehrli FW, Song HK, Saha PK, Wright AC. Quantitative MRI for the assessment of bone structure and function. NMR Biomed. 2006;19(7):731–764. [PubMed]
25. Magland JF, Wehrli FW. Trabecular bone structure analysis in the limited spatial resolution regime of in vivo MRI. Acad Radiol. 2008;15(12):1482–1493. [PMC free article] [PubMed]
26. Faulkner KG, Gluer CC, Majumdar S, Lang P, Engelke K, Genant HK. Noninvasive measurements of bone mass, structure, and strength: current methods and experimental techniques. AJR Am J Roentgenol. 1991;157(6):1229–1237. [PubMed]
27. Genant HK, Lang TF, Engelke K, et al. Advances in the Noninvasive Assessment of Bone Density, Quality, and Structure. Calcif Tissue Int. 1996;59(7):10–15. [PubMed]
28. Folkesson J, Carballido-Gamio J, Eckstein F, Link TM, Majumdar S. Local bone enhancement fuzzy clustering for segmentation of MR trabecular bone images. Med Phys. 2010;37(1):295–302. [PubMed]
29. Issever AS, Link TM, Kentenich M, et al. Trabecular bone structure analysis in the osteoporotic spine using a clinical in vivo setup for 64-slice MDCT imaging: comparison to microCT imaging and microFE modeling. J Bone Miner Res. 2009;24(9):1628–1637. [PubMed]
30. Huber MB, Carballido-Gamio J, Fritscher K, et al. Development and testing of texture discriminators for the analysis of trabecular bone in proximal femur radiographs. Med Phys. 2009;36(11):5089–5098. [PubMed]
31. Genant HK, Majumdar S. High-resolution magnetic resonance imaging of trabecular bone structure. Osteoporos Int. 1997;7 (Suppl 3):S135–139. [PubMed]
32. Techawiboonwong A, Song HK, Magland JF, Saha PK, Wehrli FW. Implications of pulse sequence in structural imaging of trabecular bone. J Magn Reson Imaging. 2005;22(5):647–655. [PubMed]
33. Banerjee S, Choudhury S, Han ET, et al. Autocalibrating parallel imaging of in vivo trabecular bone microarchitecture at 3 Tesla. Magn Reson Med. 2006;56(5):1075–1084. [PubMed]
34. Magland J, Vasilic B, Wehrli FW. Fast low-angle dual spin-echo (FLADE): a new robust pulse sequence for structural imaging of trabecular bone. Magn Reson Med. 2006;55(3):465–471. [PubMed]
35. Krug R, Carballido-Gamio J, Burghardt AJ, et al. Assessment of trabecular bone structure comparing magnetic resonance imaging at 3 Tesla with high-resolution peripheral quantitative computed tomography ex vivo and in vivo. Osteoporos Int. 2008;19(5):653–661. [PubMed]
36. Siffert RS, Luo GM, Cowin SC, Kaufman JJ. Dynamic relationships of trabecular bone density, architecture, and strength in a computational model of osteopenia. Bone. 1996;18(2):197–206. [PubMed]
37. Hwang SN, Wehrli FW, Williams JL. Probability-based structural parameters from three-dimensional nuclear magnetic resonance images as predictors of trabecular bone strength. Med Phys. 1997;24(8):1255–1261. [PubMed]
38. Ouyang X, Selby K, Lang P, et al. High resolution magnetic resonance imaging of the calcaneus: age-related changes in trabecular structure and comparison with dual X-ray absorptiometry measurements. Calcif Tissue Int. 1997;60(2):139–147. [PubMed]
39. McCalden RW, McGeough JA, Court-Brown CM. Age-related changes in the compressive strength of cancellous bone. The relative importance of changes in density and trabecular architecture. J Bone Joint Surg Am. 1997;79(3):421–427. [PubMed]
40. Wehrli FW, Hwang SN, Ma J, Song HK, Ford JC, Haddad JG. Cancellous bone volume and structure in the forearm: noninvasive assessment with MR microimaging and image processing. Radiology. 1998;206(2):347–357. [PubMed]
41. Link TM, Majumdar S, Grampp S, et al. Imaging of trabecular bone structure in osteoporosis. Eur Radiol. 1999;9(9):1781–1788. [PubMed]
42. Pothuaud L, Porion P, Lespessailles E, Benhamou CL, Levitz P. A new method for three-dimensional skeleton graph analysis of porous media: application to trabecular bone microarchitecture. J Microsc. 2000;199(Pt 2):149–161. [PubMed]
43. 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 osteoporosis. J Bone Miner Res. 2001;16(8):1520–1531. [PubMed]
44. Laib A, Newitt DC, Lu Y, Majumdar S. New model-independent measures of trabecular bone structure applied to in vivo high-resolution MR images. Osteoporos Int. 2002;13(2):130–136. [PubMed]
45. Jinnai H, Watashiba H, Kajihara T, Nishikawa Y, Takahashi M, Ito M. Surface curvatures of trabecular bone microarchitecture. Bone. 2002;30(1):191–194. [PubMed]
46. Cao H, Nazarian A, Ackerman JL, et al. Quantitative 31P NMR spectroscopy and 1H MRI measurements of bone mineral and matrix density differentiate metabolic bone diseases in rat models. Bone. 2010;46(6):1582–1590. [PMC free article] [PubMed]
47. Ma J, Wehrli FW. Method for image-based measurement of the reversible and irreversible contribution to the transverse-relaxation rate. J Magn Reson B. 1996;111(1):61–69. [PubMed]
48. Song HK, Wehrli FW, Ma J. Field strength and angle dependence of trabecular bone marrow transverse relaxation in the calcaneus. J Magn Reson Imaging. 1997;7(2):382–388. [PubMed]
49. Wehrli FW, Ma J, Hopkins JA, Song HK. Measurement of R’2 in the presence of multiple spectral components using reference spectrum deconvolution. J Magn Reson. 1998;131(1):61–68. [PubMed]
50. Chung H, Wehrli FW, Williams JL, Kugelmass SD. Relationship between NMR transverse relaxation, trabecular bone architecture, and strength. Proc Natl Acad Sci U S A. 1993;90(21):10250–10254. [PubMed]
51. Chung HW, Hwang SN, Yeung HN, Wehrli FW. Mapping of the magnetic-field distribution in cancellous bone. J Magn Reson B. 1996;113(2):172–176. [PubMed]
52. Bouchard LS, Wehrli FW, Chin CL, Warren WS. Structural anisotropy and internal magnetic fields in trabecular bone: coupling solution and solid dipolar interactions. J Magn Reson. 2005;176(1):27–36. [PubMed]
53. Capuani S, Curzi F, Alessandri FM, Maraviglia B, Bifone A. Characterization of trabecular bone by dipolar demagnetizing field MRI. Magn Reson Med. 2001;46(4):683–689. [PubMed]
54. Capuani S, Alessandri FM, Bifone A, Maraviglia B. Multiple spin echoes for the evaluation of trabecular bone quality. Magma. 2002;14(1):3–9. [PubMed]
55. Ma J, Wehrli FW, Song HK. Fast 3D large-angle spin-echo imaging (3D FLASE) Magn Reson Med. 1996;35(6):903–910. [PubMed]
56. Song YQ, Ryu S, Sen PN. Determining multiple length scales in rocks. Nature. 2000;406(6792):178–181. [PubMed]
57. Song YQ. Pore sizes and pore connectivity in rocks using the effect of internal field. Magn Reson Imaging. 2001;19(3–4):417–421. [PubMed]
58. Song YQ. Using Internal Magnetic Fields to Obtain Pore Size Distributions of Porous Media. Concepts in Magnetic Resonance Part A. 2003;18A(2):97–110.
59. Hurlimann MD. Effective Gradients in Porous Media Due to Susceptibility Differences. J Magn Reson. 1998;131(2):232–240. [PubMed]
60. Lisitza NV, Song YQ. The behavior of diffusion eigenmodes in the presence of internal magnetic field in porous media. J Chem Phys. 2001;114(20):9120–9124.
61. Sigmund EE, Cho H, Chen P, et al. Diffusion-based MR methods for bone structure and evolution. Magn Reson Med. 2008;59(1):28–39. [PubMed]
62. Sigmund EE, Cho H, Song YQ. High-resolution MRI of internal field diffusion-weighting in trabecular bone. NMR Biomed. 2009;22(4):436–448. [PubMed]
63. Jaramillo D. Whole-body MR imaging, bone diffusion imaging: how and why? Pediatr Radiol. 2010;40(6):978–984. [PubMed]
64. Ababneh ZQ, Beloeil H, Berde CB, Ababneh AM, Maier SE, Mulkern RV. In vivo lipid diffusion coefficient measurements in rat bone marrow. Magn Reson Imaging. 2009;27(6):859–864. [PubMed]
65. Lehnert A, Machann J, Helms G, Claussen CD, Schick F. Diffusion characteristics of large molecules assessed by proton MRS on a whole-body MR system. Magn Reson Imaging. 2004;22(1):39–46. [PubMed]
66. Yeung DK, Wong SY, Griffith JF, Lau EM. Bone marrow diffusion in osteoporosis: evaluation with quantitative MR diffusion imaging. J Magn Reson Imaging. 2004;19(2):222–228. [PubMed]
67. Bernstein MA, King KF, Zhou XJ. Handbook of MRI Pulse Sequences. Elsevier Science; 2004.
68. Haacke ME, Brown RW, Thompson MR, Ramesh V. Magnetic Resonance Imaging: Physical Principles and Sequence Design. New York: Wiley; 1999.
69. Bevington PR, Robinson DK. Data Reduction and Error Analysis for the Physical Sciences. New York: Mc Graw Hill Higher Education; 2002.
70. Transtrum MK, Machta BB, Sethna JP. Why are nonlinear fits to data so challenging? Phys Rev Lett. 2010;104(6):060201. [PubMed]
71. Cohen J. The statistical power of abnormal-social psychological research: a review. J Abnorm Soc Psychol. 1962;65:145–153. [PubMed]
72. Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes: The Art of Scientific Computing. Cambridge University Press; 2007.
73. Nocedal J, Wright S. Numerical Optimization. New York: Springer; 2006.