Bone tissue and porous media: common features and differences studied by NMR relaxation

Despite significant differences between bone tissues and other porous media such as oilfield rocks, there are common features as well as differences in the response of NMR relaxation measurements to the internal structures of the materials. Internal surfaces contribute to both transverse (T2) and longitudinal (T1) relaxation of pore fluids, and in both cases the effects depend on, among other things, local surface-to-volume ratio (S/V). In both cases variations in local S/V can lead to distributions of relaxation times, sometimes over decades. As in rocks, it is useful to take bone data under different conditions of cleaning, saturation, and desaturation. T1 and T2 distributions are computed using UPEN. In trabecular bone it is easy to see differences in dimensions of intertrabecular spaces in samples that have been de-fatted and saturated with water, with longer T1 and T2 for larger pores. Both T1 and T2 distributions for these water-saturated samples are bimodal, separating or partly separating inter- and intratrabecular water. The T1 peak times have a ratio of from 10 to 30, depending on pore size, but for the smaller separations the distributions may not have deep minima. The T2 peak times have ratios of over 1000, with intratrabecular water represented by large peaks at a fraction of a ms, which we can observe only by single spin echoes. CPMG data show peaks at about a second, tapering down to small amplitudes by a ms. In all samples the free induction decay (FID) from an inversion-recovery (IR) T1 measurement shows an approximately Gaussian (solid-like) component, exp[-1/2 (T/TGC), with TGC approximately 11.7+/-0.7 micros (GC for "Gaussian Component"), and a liquid-like component (LLC) with initially simple-exponential decay at the rate-average time T(2-FID) for the first 100 micros. Averaging and smoothing procedures are adopted to derive T(2-FID) as a function of IR time and to get T1 distributions for both the GC and the LLC. It appears that contact with the GC, which is presumed to be 1H on collagen, leads to the T2 reduction of at least part of the LLC, which is presumed to be water. Progressive drying of the cleaned and water-saturated samples confirms that the long T1 and T2 components were in the large intertrabecular spaces, since the corresponding peaks are lost. Further drying leads to further shortening of T2 for the remaining water but eventually leads to lengthening of T1 for both the collagen and the water. After the intertrabecular water is lost by drying, T1 is the same for GC and LLC. T(2-FID) is found to be roughly 320/alpha micros, where alpha is the ratio of the extrapolated GC to LLC, appearing to indicate a time tau of about 320 micros for 1H transverse magnetization in GC to exchange with that of LLC. This holds for all samples and under all conditions investigated. The role of the collagen in relaxation is confirmed by treatment to remove the mineral component, observing that the GC remains and has the same TGC and has the same effect on the relaxation times of the associated water. Measurements on cortical bone show the same collagen-related effects but do not have the long T1 and T2 components.     

Continuous distribution analysis of marrow 1H magnetic resonance relaxation in bone

Longitudinal and transverse NMR relaxation of 1H nuclei were studied in vitro on fresh animal femur samples. A large number of data points were taken, starting at 100 micros for T(1) by inversion-recovery, at 200 micros for T(2) by single-echo sequences, and at 600 micros for T(2) by CPMG echo-trains. Quasi-continuous distributions of relaxation times were computed, giving wide distributions for all samples. Bulk marrow removed from the medullary cavity showed T(2) distributions from about 20 ms to 600 ms and T(1) distributions from about 40 ms to 2 s. The 1H nuclei in trabecular bone samples, where marrow is confined, may show long tails for T(2) at relaxation times down to 250 micros, the origin of which is still not known. These tails are absent in bulk marrow from the medullary cavity. The differences observed in T(1) distributions among trabecular bone samples are in accordance with the different marrow compositions. Discrete exponential fits were computed also, and in most cases four discrete exponential components were required to fit the experimental data adequately. However, the discrete components do not seem to correspond to any physically distinguishable separate compartments.     

Two-exponential analysis of spin-spin proton relaxation times in MR imaging using surface coils

Proton relaxation time measurements were performed on a standard whole body MR imager operating at 1.5 T using a conventional surface coil of the manufacturer. A combined CP/CPMG multiecho, multislice sequence was used for the T1 and T2 relaxation time measurements. Two repetition times of 2000 ms (30 echoes) and 600 ms (2 echoes) with 180 degrees-pulse intervals of 2 tau = 22 ms were interleaved in this sequence. A two-exponential T2 analysis of each pixel of the spin-echo images was computed in a case of an acoustic neurinoma. The two-exponential images show a "short" component (T2S) due to white and gray matter and a "long" component (T2S) due to the cerebrospinal fluid. In the fatty tissue two components with T2S = 35 +/- 3 ms and T2L = 164 +/- 7 ms were measured. Comparing with Gd-DTPA imaging the relaxation time images show a clear differentiation of vital tumor tissue and cerebrospinal fluid.     

Electron spin relaxation of triarylmethyl radicals in fluid solution

Electron spin relaxation times of a Nycomed triarylmethyl radical (sym-trityl) in water, 1:1 water:glycerol, and 1:9 water:glycerol were measured at L-band, S-band, and X-band by pulsed EPR methods. In H(2)O solution, T(1) is 17+/-1 micros at X-band at ambient temperature, is nearly independent of microwave frequency, and exhibits little dependence on viscosity. The temperature dependence of T(1) in 1:1 water:glycerol is characteristic of domination by a Raman process between 20 and 80 K. The increased spin-lattice relaxation rates at higher temperatures, including room temperature, are attributed to a local vibrational mode that modulates spin-orbit coupling. In H(2)O solution, T(2) is 11+/-1 micros at X-band, increasing to 13+/-1 micros at L-band. For more viscous solvent mixtures, T(2) is much shorter than T(1) and weakly frequency dependent, which indicates that incomplete motional averaging of hyperfine anisotropy makes a significant contribution to T(2). In water and 1:1 water:glycerol solutions continuous wave EPR linewidths are not relaxation determined, but become relaxation determined in the higher viscosity 1:9 water:glycerol solutions. The Lorentzian component of the 250-MHz linewidths as a function of viscosity is in good agreement with T(2)-determined contributions to the linewidths at higher frequencies.     

Continuous wave MRI of heterogeneous materials

A prototype continuous wave MRI system operating at 7T has been used successfully to study a variety of heterogeneous materials exhibiting T2 relaxation values ranging from 10 micros to 50 ms. Two-dimensional images of a poly(methly methacrylate) (PMMA) resolution phantom (T2=38 micros) exhibited a spatial resolution of approximately 1mm at a magnetic field gradient strength of 200 mT/m. The technique was used to study the hydration, drying, and subsequent water penetration properties of cement samples made from ordinary Portland cement, and revealed inhomogeneities arising from the cure conditions. Sandstone samples from an oil reservoir in the North Sea were also studied; structure within these materials, arising from the sedimentary bed layering in the reservoir, was found to have an effect on their water transport properties. A section from a confectionery bar (T2* approximately 50-60 ms) was also imaged, and its internal structure could be clearly discerned.     

Proton magnetic relaxation in bone marrow related to age and bone mineral density: low-resolution in vitro studies.

Detailed analysis of proton spin-spin and spin-lattice relaxation behaviors of the bone marrow in the presence of trabecular bone network was performed at low-resolution (B(0) = 0.496T) on rat vertebrae specimens deprived of spinal cord. Two groups of samples, from young and old healthy animals, were investigated before cellular necrosis had started. BMD measurements were carried out to quantify the expected age-related modifications of the trabecular bone network. 1H-MR measurements were also performed on the same samples, deprived of marrow and saturated with water, in order to control the validity of a possible interpretation of the marrow 1H-MR characteristics, in terms of marrow components, and to investigate the possible employment of these samples to study the trabecular bone network properties. We pointed out that: 1) a bimodal distribution of T(2i) and T(1i) values (distinguishing "fast" and "slow" relaxations) describes satisfactorily all the 1H-MR experimental decays; 2) age-related modifications of the trabecular bone network are marked by correlate variations of the BMD value and of the proton spin-spin relaxation rates in water saturated samples; 3) age-related modifications of marrow are underlined by variations of the average value of the "fast" T(2i) and of the "slow" T(1i) relaxation time distributions, which could be attributed to the marrow components different from the fat granules of the adipose cells.Our results suggest that studies in vitro on bone tissue, by 1H-MR techniques at low-resolution, may contribute to a better bone function characterization and, therefore, to a better clinical utilization of MRI techniques.     

Constrained modeling for spectroscopic measurement of bi-exponential spin-lattice relaxation of water in vivo

The (1)H NMR water signal from spectroscopic voxels localized in gray matter contains contributions from tissue and cerebral spinal fluid (CSF). A typically weak CSF signal at short echo times makes separating the tissue and CSF spin-lattice relaxation times (T(1)) difficult, often yielding poor precision in a bi-exponential relaxation model. Simulations show that reducing the variables in the T(1) model by using known signal intensity values significantly improves the precision of the T(1) measurement. The method was validated on studies on eight healthy subjects (four males and four females, mean age 21 +/- 2 years) through a total of twenty-four spectroscopic relaxation studies. Each study included both T(1) and spin-spin relaxation (T(2)) experiments. All volumes were localized along the Sylvian fissure using a stimulated echo localization technique with a mixing time of 10 ms. The T(2) experiment consisted of 16 stimulated echo acquisitions ranging from a minimum echo time (TE) of 20 ms to a maximum of 1000 ms, with a repetition time of 12 s. All T(1) experiments consisted of 16 stimulated echo acquisition, using a homospoil saturation recovery technique with a minimum recovery time of 50 ms and a maximum 12 s. The results of the T(2) measurements provided the signal intensity values used in the bi-exponential T(1) model. The mean T(1) values when the signal intensities were constrained by the T(2) results were 1055.4 ms +/- 7.4% for tissue and 5393.5 ms +/- 59% for CSF. When the signal intensities remained free variables in the model, the mean T(1) values were 1085 ms +/- 19.4% and 5038.8 ms +/- 113.0% for tissue and CSF, respectively. The resulting improvement in precision allows the water tissue T(1) value to be included in the spectroscopic characterization of brain tissue.     

Combined MR-relaxation and MR-cryoporometry in the study of bone microstructure.

MR-Relaxation (MRR) of 1H nuclei and MR-Cryoporometry (MRC) are combined to assess their feasibility and their potential in the study of bone microstructure. In principle, both techniques are able to give information on the structure of the pore space confining the fluids. Cow femur samples were carefully cored and cleaned in order to remove the natural fluids inside. For MRR analysis quasi-continuous distributions of T(1) and T(2) were obtained on samples fully saturated with water. Cyclohexane was used as a saturating fluid for MRC analysis. All T(1) and T(2) quasi-continuous distributions of water confined in bone samples are more than three decades wide, showing sufficient details to differentiate the samples. Pore size distributions obtained by MRC also differentiate the samples showing different characteristics of the pore space structure in the range of the highest sensitivity of the method (typically 3 to 100 nm, mesopore range). In particular, in samples where MRR shows a large fraction of signal with relaxation times below 10(2) ms, MRC indicates a large fraction of pore volume with pore sizes in the mesopore range.     

A feasibility study of in vivo T1rho imaging of the intervertebral disc

PURPOSE: Recent studies have proposed that magnetic resonance (MR) T1rho relaxation time is associated with loss of macromolecules. The depletion of macromolecules in the matrix of the intervertebral disc may be an initiating factor in degenerative disc disease. The purpose of this study was to test the feasibility of quantifying T1rho relaxation time in phantoms and intervertebral discs of healthy volunteers using in vivo MR imaging at 3 T. MATERIALS AND METHODS: A multislice T1rho spiral sequence was used to quantify T1rho relaxation time in phantoms with different agarose concentrations and in the intervertebral discs of 11 healthy volunteers (mean age=31.3 years; age range=23-60 years; gender: 5 females, 6 males). RESULTS: The phantom studies demonstrated the feasibility of using spiral imaging at 3 T. The in vivo results indicate that the median T1rho value of the nucleus (116.6+/-21.4 ms) is significantly greater (P<0.05) than that of the annulus (84.1+/-11.7 ms). The correlations between the age of the volunteers and T1rho relaxation time in the nucleus (r2=-0.82; P=0.0001) and the annulus (r2=-0.37; P=0.04) were significant. A trend of decreasing T1rho values from L3-4 to L4-5 to L5-S1 was evident. CONCLUSION: The results of this study suggest that in vivo T1rho quantification is feasible and may potentially be a clinical tool in identifying early degenerative changes in the intervertebral disc.     

From porous media to trabecular bone relaxation analysis: spatial variation of marrow 1H relaxation time distributions detected in vitro by quasi-continuous distribution analysis.

Quasi-continuous distributions of T(1) and T(2) of 1H nuclei were analyzed in vitro at 20MHz on some twenty fresh bone samples of pig femur. Large numbers of data points allowed a detailed investigation. Relaxation data were inverted by UPEN (Uniform PENalty inversion). In all samples the widths of the distributions, covering more than two decades, are not even close to being compatible with single exponential components. Moreover, the T(1) and T(2) distributions show enough character to distinguish the samples. We observe a spatial variation of these characteristics and in particular a second peak centered at 500-600 ms appearing in some proximal femur samples. The quasi-continuous distribution allows one to correlate the water content of the sample with parts of the distributions in specific time ranges. The signal fraction with T(1) values longer than a cutoff time of about 170 ms is in very good agreement with the water content of the samples and is significantly larger in the group of samples cored from proximal femur. Also T(2) distributions differentiate the samples, and the signal fraction with T(2) shorter than about 30 ms is significantly larger in the group of distal femur samples.