Multi-components of T2 relaxation in ex vivo cartilage and tendon

The multi-components of T₂ relaxation in cartilage and tendon were investigated by microscopic MRI (μMRI) at 13 and 26 μm transverse resolutions. Two imaging protocols were used to quantify T₂ relaxation in the specimens, a 5-point sampling and a 60-point sampling. Both multi-exponential and non-negative-least-square (NNLS) fitting methods were used to analyze the μMRI signal. When the imaging voxel size was 6.76 × 10⁻⁴ mm³ and within the limit of practical signal-to-noise ratio (SNR) in microscopic imaging experiments, we found that (1) canine tendon has multiple T₂ components; (2) bovine nasal cartilage has a single T₂ component; and (3) canine articular cartilage has a single T₂ component. The T₂ profiles from both 5-point and 60-point methods were found to be consistent in articular cartilage. In addition, the depletion of the glycosaminoglycan component in cartilage by the trypsin digestion method was found to result in a 9.81–20.52% increase in T₂ relaxation in articular cartilage, depending upon the angle at which the tissue specimen was oriented in the magnetic field.     

The zonal architecture of human articular cartilage described by T2 relaxation time in the presence of Gd-DTPA2-

The depth-wise variation of T(2) relaxation time is known to reflect the collagen network architecture in cartilage, while the delayed Gadolinium Enhanced MRI of Cartilage (dGEMRIC) technique is sensitive to tissue proteoglycan (PG) concentration. As the cartilage PG content varies along the tissue depth, the depth-dependent accumulation of the contrast agent may affect the inherent T(2) of cartilage in a nonconstant manner. Therefore, T(2) and dGEMRIC are typically measured in separate MRI sessions. In the present in vitro MRI study at 9.4 T, depth-wise T(2) profiles and collagenous zone thicknesses as determined from T(2) maps in the absence and presence of Gd-DTPA(2-) (T(2) and T(2Gd), respectively) were compared in samples of intact human articular cartilage (n=65). These T(2) measures were further correlated with birefringence (BF) of polarized light microscopy (PLM) to quantify the ability of MRI to predict the properties of the collagen fibril network. The reproducibility of the T(2) measurement in the current setup was also studied. Typical tri-laminar collagen network architecture was observed both with and without Gd-DTPA(2-). The inverse of BF (1/BF) correlated significantly with both T(2) and T(2Gd) (r=0.91, slope=0.56 and r=0.90, slope=0.63), respectively. The statistically significant linear correlations between zone thicknesses as determined from T(2) and T(2Gd) were r=0.55 (slope=0.49), r=0.74 (slope=0.71) and r=0.95 (slope=0.94) for superficial, middle and deep tissue zones, respectively. Reproducibility of the T(2) measurement was worst for superficial cartilage. Consistent with PLM, T(2) and T(2Gd) measurements reveal highly similar depth-dependent information on collagen network in intact human cartilage. Thus, dGEMRIC and T(2) measurements in one MRI session are feasible for intact articular cartilage in vitro.     

Hydration water dynamics in biopolymers from NMR relaxation in the rotating frame

Assuming dipole-dipole interaction as the dominant relaxation mechanism of protons of water molecules adsorbed onto macromolecule (biopolymer) surfaces we have been able to model the dependences of relaxation rates on temperature and frequency. For adsorbed water molecules the correlation times are of the order of 10(-5)s, for which the dispersion region of spin-lattice relaxation rates in the rotating frame R(1)(ρ)=1/T(1)(ρ) appears over a range of easily accessible B(1) values. Measurements of T(1)(ρ) at constant temperature and different B(1) values then give the "dispersion profiles" for biopolymers. Fitting a theoretical relaxation model to these profiles allows for the estimation of correlation times. This way of obtaining the correlation time is easier and faster than approaches involving measurements of the temperature dependence of R(1)=1/T(1). The T(1)(ρ) dispersion approach, as a tool for molecular dynamics study, has been demonstrated for several hydrated biopolymer systems including crystalline cellulose, starch of different origins (potato, corn, oat, wheat), paper (modern, old) and lyophilized proteins (albumin, lysozyme).     

Assessment of interstitial water content of articular cartilage with T1 relaxation

The interstitial water content typically increases in the early degeneration of articular cartilage. Previously, T₂ relaxation has been related to water content, yet it is known to be strongly affected by the collagen orientation. Articular cartilage plugs from the bovine patella, femur and tibia (N=20) were mapped for T₁ and T₂ at 9.4 T to test the ability of T₁ relaxation to reflect cartilage water content. As a reference, water and proteoglycan (PG) contents were determined. Significant (P<.01) linear associations were demonstrated between the relaxation rates and tissue water content (R₁: r=−.81, R₂: r=−.60) and PG content (R₁: r=.75). After adjustment for the tissue water content, partial correlation analysis did not show significant associations between the relaxation rates and tissue PG content. After the effect of PGs was removed, significant (P<.05) linear correlation between the relaxation rates and tissue water content (R₁: r=−.48, R₂: r=−.50) was observed. Thus, the spin-lattice relaxation rate is proposed to provide a biomarker for water content in articular cartilage.     

Regional variations in MR relaxation of hip joint cartilage in subjects with and without femoralacetabular impingement

The objective of this study was to analyze regional variations of magnetic resonance (MR) relaxation times (T_1ρ and T₂) in hip joint cartilage of healthy volunteers and subjects with femoral acetabular impingement (FAI). Morphological and quantitative images of the hip joints of 12 healthy volunteers and 9 FAI patients were obtained using a 3 T MR scanner. Both femoral and acetabular cartilage layers in each joint were semi-automatically segmented on sagittal 3D high-resolution spoiled gradient echo (SPGR) images. These segmented regions of interest (ROIs) were automatically divided radially into twelve equal sub-regions (30⁰ intervals) based on the fitted center of the femur head. The mean value of T_1ρ/T₂ was calculated in each sub-region after superimposing the divided cartilage contours on the MR relaxation (T_1ρ/T₂) maps to quantify the relaxation times. T_1ρ and T₂ relaxation times of the femoral cartilage were significantly higher in FAI subjects compared to healthy controls (39.9 ± 3.3 msec in FAI vs. 35.4 ± 2.3 msec in controls for T_1ρ (P = 0.0020); 33.9 ± 3.1 msec in FAI vs. 31.1 ± 1.7 msec in controls for T₂ (P = 0.0160)). Sub-regional analysis showed significantly different T_1ρ and T₂ relaxation times in the anterior-superior region (R9) of the hip joint cartilage between subjects with FAI and healthy subjects, suggesting possible regional differences in cartilage matrix composition between these two groups. Receiver operating characteristic (ROC) analysis showed that sub-regional analysis in femoral cartilage was more sensitive in discriminating FAI joint cartilage from that of healthy joints than global analysis of the whole region (T_1ρ: area under the curve (AUC) = 0.981, P = 0.0001 for R9 sub-region; AUC = 0.901, P = 0.002 for whole region; T₂: AUC = 0.976, P = 0.0005 for R9 sub-region; AUC = 0.808, P = 0.0124 for whole region). The results of this study demonstrated regional variations in hip cartilage composition using MR relaxation times (T_1ρ and T₂) and suggested that analysis based on local regions was more sensitive than global measures in subjects with and without FAI.     

Single-scan longitudinal relaxation measurements in high-resolution NMR spectroscopy

Several single-scan experiments for the measurement of the longitudinal relaxation time (T1) are proposed. These experiments result in fast and accurate determinations of the relaxation rate, are relatively robust to pulse imperfections, and preserve information about the chemical shift. The method used in these experiments is to first encode the T1 values as a spatial variation of the magnetization and then to read out this variation either by applying a weak gradient during acquisition or by sequentially observing different slices of the sample. As a result, it is possible to reduce the time necessary to determine the T1 values by one or two orders of magnitude. This time saving comes at the expense of the signal-to-noise level of the resulting spectrum and some chemical shift resolution.     

Improvement of duty-cycle heating compensation in NMR spin relaxation experiments

To reliably measure NMR relaxation properties of macromolecules is a prerequisite for precise experiments that identify subtle variations in relaxation rates, as required for the determination of rotational diffusion anisotropy, CSA tensor determination, advanced motional modeling or entropy difference estimations. An underlying problem with current NMR relaxation measurement protocols is maintaining constant sample temperature throughout the execution of the relaxation series especially when rapid data acquisition is required. Here, it is proposed to use a combination of a heating compensation and a proton saturation sequence at the beginning of the NMR relaxation pulse scheme. This simple extension allows reproducible, robust and rapid acquisition of NMR spin relaxation data sets. The method is verified with (15)N spin relaxation measurements for human ubiquitin.     

2H T2rho relaxation dynamics and double-quantum filtered NMR studies

In this study 2H T2rho DQF NMR spectra of water in MCM-41 were measured. The T2rho double-quantum filtered (DQF) NMR signal is generated by applying a radio frequency (RF) field for various durations and then observed after a monitor RF pulse. It was found that the transfer between different quantum coherences by the couplings during long-duration RF fields (i.e., soft pulses) and that residual quadrupolar interaction dominates the signal decay. Knowledge of coherence transfer during long-RF pulses has special significance for the development of sophisticated multi-quantum NMR experiments especially multi-quantum MRI applications.     

On the measurement of multi-component T2 relaxation in cartilage by MR spectroscopy and imaging

The multicomponent T₂ relaxation in bovine nasal cartilage (BNC) was investigated by nuclear magnetic resonance spectroscopy using the Carr-Purcell-Meiboom-Gill (CPMG) sequence and microscopic magnetic resonance imaging (μMRI) method using a CPMG-SE imaging sequence. All experimental data were analyzed by the non-negative least square (NNLS) procedure. Only one T₂ component was found in BNC by both experimental methods (about 113 and 170 ms before and after being enzymatically digested by trypsin). Several experimental and specimen-related factors were investigated in this study, and it was found that some of them could produce artificial multi-component T₂, including the use of the standard MSME imaging sequence at certain imaging gradients.     

Distributions of transverse relaxation times for soft-solids measured in strongly inhomogeneous magnetic fields

The single-sided NMR-MOUSE sensor that operates in highly inhomogeneous magnetic fields is used to record a CPMG ¹H transverse relaxation decay by CPMG echo trains for a series of cross-linked natural rubber samples. Effective transverse relaxation rates 1/T_2,short and 1/T_2,long were determined by a bi-exponential fit. A linear dependence of transverse relaxation rates on cross-link density is observed for medium to large values of cross-link density. As an alternative to multi-exponential fits the possibility to analyze the dynamics of soft polymer network in terms of multi-exponential decays via the inverse Laplace transformation was studied. The transient regime and the effect of the T₁/T₂ ratio in inhomogeneous static and radiofrequency magnetic fields on the CPMG decays were studied numerically using a dedicated C++ program to simulate the temporal and spatial dependence of the CPMG response. A correction factor T₂/T_2,eff is derived as a function of the T₁/T₂ ratio from numerical simulations and compared with earlier results from two different well logging devices. High-resolution T₁–T₂ correlations maps are obtained by two-dimensional Laplace inversion of CPMG detected saturation recovery curves. The T₁–T₂ experimental correlations maps were corrected for the T₁/T₂ effect using the derived T₂/T_2,eff correction factor.     

Two-dimensional correlation relaxation studies of cement pastes

Two-dimensional nuclear magnetic resonance relaxation correlation studies of cement pastes have been performed on a unilateral magnet, the Surface GARField. Through these measurements, the hydration process can be observed by monitoring the evolution of porosity. Characteristic relaxation time distributions have been observed in different cement pastes: fresh white cement, prehydrated white cement and ordinary Portland cement. The observed T(1)/T(2) ratio in these cements has been shown to agree with expectations based on high field values.     

Spin-lattice relaxation and a fast T1-map acquisition method in MRI with transient-state magnetization

The magnetization under the spin-lattice relaxation and the nuclear magnetic resonance radiofrequency (RF) pulses is calculated for a signal RF pulse train and for a sequence of multiple RF pulse-trains. It is assumed that the transverse magnetization is zero when each RF pulse is applied. The result expressions can be grouped into two terms: a decay term, which is proportional to the initial magnetization M0, and a recovery term, which has no M0 dependence but strongly depends on the spin-lattice relaxation and the equilibrium magnetization Meq. In magnetic resonance pulse sequences using magnetization in transient state, the recovery term produces artifacts and can seriously degrade the function of the preparation sequence for slice selection, contrast weighting, phase encoding, etc. This work shows that the detrimental effect can be removed by signal averaging in an eliminative fashion. A novel fast data acquisition method for constructing the spin-lattice relaxation (T1) map is introduced. The method has two features: (i) By using eliminative averaging, the curve to fit the T1 value is a decay exponential function rather than a recovery one as in conventional techniques; therefore, the measurement of Meq is not required and the result is less susceptible to the accuracy of the inversion RF pulse. (ii) The decay exponential curve is sampled by using a sequence of multiple pulse-trains. An image is reconstructed from each train and represents a sample point of the curve. Hence a single imaging sequence can yield multiple sample points needed for fitting the T1 value in contrast to conventional techniques that require repeating the imaging sequence for various delay values but obtain only one sample point from each repetition.     

Oxygen-induced changes in longitudinal relaxation times in skeletal muscle

The objective was to measure the effect of 100% oxygen inhalation on T1 relaxation times in skeletal muscle. Healthy volunteers were scanned using three different MRI protocols while breathing medical air and 100% oxygen. Measurements of T1 were made from regions of interest (ROIs) within various skeletal muscle groups. Dynamic data of subjects breathing a sequence of air-oxygen-air allowed the calculation of characteristic wash-in and -out times for dissolved oxygen in muscle. Contrary to previous findings, a statistically significant decrease in T1 in skeletal muscle was observed due to oxygen inhalation. We report approximate baseline characteristic values for the response of skeletal muscle to oxygen inhalation. This measurement may provide new biomarkers for evaluation of oxygen delivery and consumption in normal and diseased skeletal muscle.     

Multivariate Image Analysis of Magnetic Resonance Images with the Direct Exponential Curve Resolution Algorithm (DECRA): Part 1: Algorithm and Model Study

Antalek and Windig recently presented a fast method to resolve a series of NMR mixture spectra, where the contribution of the components varies with a decaying exponential [B. Antalek and W. Windig,J. Am. Chem. Soc.118, 10,331–10,332 (1996); W. Windig and B. Antalek,Chemom. Intell. Lab. Syst.37, 241–254 (1997)]. The method was called DECRA (direct exponential curve resolution algorithm). In this paper DECRA will be applied to two series of magnetic resonance images. The signal of one series is based uponT₂relaxation, and the other is based uponT₁relaxation. In order to evaluate the technique, the magnetic resonance images of a phantom where used. A transformation is introduced to enable the application of DECRA to aT₁series of magnetic resonance images. A separate paper in this issue will describe the application of the techniques to magnetic resonance images of the human brain.     

Formation of p-cresol:piperazine complex in solution monitored by spin-lattice relaxation times and pulsed field gradient NMR diffusion measurements

A study of the nature of the anthelmintic p-cresol:piperazine complex in chloroform solution has been conducted using different NMR techniques: self-diffusion coefficients using DOSY; NOE, NULL, and double-selective T1 measurements to determine inter-molecular distances; and selective and non-selective T1 measurements to determine correlation times. The experimental results in solution and CP-MAS were compared to literature X-ray diffraction data using molecular modeling. It was shown that the p-cresol:piperazine complex exists in solution in a very similar manner as it does in the solid state, with one p-cresol molecule hydrogen bonded through the hydroxyl hydrogen to each nitrogen atom of piperazine. The close correspondence between the X-ray diffraction data and the inter-proton distances obtained by NULL and double selective excitation techniques indicate that those methodologies can be used to determine inter-molecular distances in solution.     

A comparative study at 3 T of sequence dependence of T2 quantitation in the knee

Objective

T₂ mapping has been used widely in detecting cartilage degeneration in osteoarthritis. Several scanning sequences have been developed in the determination of T₂ relaxation times of tissues. However, the derivation of these times may vary from sequence to sequence. This study seeks to evaluate the sequence-dependent differences in T₂ quantitation of cartilage, muscle, fat and bone marrow in the knee joint at 3 T.

Methods

Three commercial phantoms and 10 healthy volunteers were studied using 3 T MR. T₂ relaxation times of the phantoms, cartilage, muscle, subcutaneous fat and marrow were derived using spin echo (SE), multiecho SE (MESE), fast SE (FSE) with varying echo train length (ETL), spiral and spoiler gradient (SPGR) sequences. The differences between these times were then evaluated using Student's t test. In addition, the signal-to-noise ratio (SNR) efficiency and coefficient of variation of T₂ from each sequence were calculated.

Results

The average T₂ relaxation time was 36.38±5.76 ms in cartilage and 34.08±6.55 ms in muscle, ranging from 27 to 45 ms in both tissues. The times for subcutaneous fat and marrow were longer and more varying, ranging from 41 to 143 ms and from 42 to 160 ms, respectively. In FSE acquisition, relaxation time significantly increases as ETL increases (P<.05). In cartilage, the SE acquisition yields the lowest T₂ values (27.52±3.10 ms), which is significantly lower than those obtained from other sequences (P<.002). T₂ values obtained from spiral acquisition (38.27±6.45 ms) were higher than those obtained from MESE (34.35±5.62 ms) and SPGR acquisition (31.64±4.53 ms). These differences, however, were not significant (P>.05).

Conclusion

T₂ quantification can be a valuable tool for the diagnosis of degenerative disease. Several different sequences exist to quantify the relaxation times of tissues. Sequences range in scan time, SNR efficiency, reproducibility and two- or three-dimensional mapping. However, when choosing a sequence for quantitation, it is important to realize that several factors affect the measured T₂ relaxation time.     

Effects of off-resonance irradiation, cross-relaxation, and chemical exchange on steady-state magnetization and effective spin-lattice relaxation times

In the presence of an off-resonance radiofrequency field, recovery of longitudinal magnetization to a steady state is not purely monoexponential. Under reasonable conditions with zero initial magnetization, recovery is nearly exponential and an effective relaxation rate constant R_1eff = 1/T_1eff can be obtained. Exact and approximate formulas for R_1eff and steady-state magnetization are derived from the Bloch equations for spins undergoing cross-relaxation and chemical exchange between two sites in the presence of an off-resonance radiofrequency field. The relaxation formulas require that the magnetization of one spin is constant, but not necessarily zero, while the other spin relaxes. Extension to three sites with one radiofrequency field is explained. The special cases of off-resonance effects alone and with cross-relaxation or chemical exchange, cross-relaxation alone, and chemical exchange alone are compared. The inaccuracy in saturation transfer measurements of exchange rate constants by published formulas is discussed for the creatine kinase reaction.     

Methods for correlating T1ρ and FID components in wideline H NMR studies of motionally heterogeneous polymer systems

We address the problem of correlating the observed FID and T_1ρ components in wideline ¹H relaxation measurements of motionally heterogeneous polymers, and show that different methods of data treatment can highlight different aspects of the correlations present. For a sample of polypropylene we find that the T_1ρ relaxation behaviour is driven by relaxation associated with the intermediate FID component, which strongly suggests a motionally inhomogeneous amorphous region in the sample.     

In vivo NMR T2 relaxation of experimental brain tumors in the cat: A multiparameter tissue characterization

Experimental gliomas (F98) were inoculated in cat brain for the systematic study of their in vivo T₂ relaxation time behavior. With a CPMG multi-echo imaging sequence, a train of 16 echoes was evaluated to obtain the transverse relaxation time and the magnetization M(0) at time t = 0. The magnetization decay curves were analyzed for biexponentiality. All tissues showed monoexponential T₂, only that of the ventricular fluid and part of the vital tumor tissue were biexponential. Based on these NMR relaxation parameters the tissues were characterized, their correct assignment being assured by comparison with histological slices. T₂ of normal grey and white matter was 74 ± 6 and 72 ± 6 msec, respectively. These two tissue types were distinguished through M(0) which for white matter was only 0.88 of the intensity of grey matter in full agreement with water content, determined from tissue specimens. At the time of maximal tumor growth and edema spread a tissue differentiation was possible in NMR relaxation parameter images. Separation of the three tissue groups of normal tissue, tumor and edema was based on T₂ with T_2(normal) < T_2(tumor) < T_2(edema). Using M(0) as a second parameter the differentiation was supported, in particular between white matter and tumor or edema. Animals were studied at 1–4 wk after tumor implantation to study tumor development. The magnetization M(0) of both tumor and peritumoral edema went through a maximum between the second and third week of tumor growth. T₂ of edema was maximal at the same time with 133 ± 4 msec, while the relaxation time of tumor continued to increase during the whole growth period, reaching values of 114 ± 12 msec at the fourth week. Thus, a complete characterization of pathological tissues with NMR relaxometry must include a detailed study of the developmental changes of these tissues to assure correct experimental conditions for the goal of optimal contrast between normal and pathological regions in the NMR images.     

Flip-back, an old trick to face highly contrasted relaxation times: application in the characterization of pharmaceutical mixtures by CPMAS NMR

The ¹³C–¹H CPMAS with flip-back pulse NMR experiment is revisited in view of applications to pharmaceutical mixtures. The analysis of the kinetics of relaxation and CP transfer with and without the flip-back pulse shows that a significant gain in ¹³C signal can be expected (thus in experimental time) from the flip-back pulse for protons with long T₁. The gain is of the order of T₁ of the protons expressed in seconds. The experiment is applied on samples with highly contrasted spin-lattice relaxation times T₁ for protons, situation encountered in pharmaceutical mixtures. The application of the flip-back increases significantly the relative signal intensity of the component with the longer T₁, making this component detectable even after using short recycle delays. Therefore, this CPMAS with flip-back experiment could be used routinely to get ¹³C CPMAS NMR spectra of mixtures in constant experimental time and signal-to-noise ratio without the need for optimization of the recycle delays, and for whatever may be the degree of crystallinity of the active principal ingredient (API) and/or excipients.