GESFIDE-PROPELLER Approach for Simultaneous R2 and R2* Measurements in the Abdomen

Purpose

To investigate the feasibility of combining GESFIDE with PROPELLER sampling approaches for simultaneous abdominal R2 and R2* mapping.

Materials and Methods

R2 and R2* measurements were performed in 9 healthy volunteers and phantoms using the GESFIDE-PROPELLER and the conventional Cartesian-sampling GESFIDE approaches.

Results

Images acquired with the GESFIDE-PROPELLER sequence effectively mitigated the respiratory motion artifacts, which were clearly evident in the images acquired using the conventional GESFIDE approach. There was no significant difference between GESFIDE-PROPELLER and reference MGRE R2* measurements (p = 0.162) whereas the Cartesian-sampling based GESFIDE methods significantly overestimated R2* values compared to MGRE measurements (p < 0.001).

Conclusion

The GESFIDE-PROPELLER sequence provided high quality images and accurate abdominal R2 and R2* maps while avoiding the motion artifacts common to the conventional Cartesian-sampling GESFIDE approaches.     

Optimal echo spacing for multi-echo imaging measurements of Bi-exponential T2 relaxation

Calculations, analytical solutions, and simulations were used to investigate the trade-off of echo spacing and receiver bandwidth for the characterization of bi-exponential transverse relaxation using a multi-echo imaging pulse sequence. The Cramer–Rao lower bound of the standard deviation of the four parameters of a two-pool model was computed for a wide range of component T₂ values and echo spacing. The results demonstrate that optimal echo spacing (TE_opt) is not generally the minimal available given other pulse sequence constraints. The TE_opt increases with increasing value of the short T₂ time constant and decreases as the ratio of the long and short time constant decreases. A simple model of TE_opt as a function of the two T₂ time constants and four empirically derived scalars is presented.     

Detection of glutamate and glutamine (Glx) by turbo spectroscopic imaging

Turbo spectroscopic imaging (TSI) is a spin echo spectroscopic imaging technique in which two or more echoes are acquired per excitation to reduce the acquisition time. The application of TSI has primarily been limited to the detection of uncoupled spins because the signal from coupled spins is modulated as a function of echo time. In this work we demonstrate how the TSI sequence can be modified to observe spins like the C₂ protons of Glx (≈3.75 ppm) which are involved solely in weak-coupling interactions. The technique exploits the chemical shift displacement effect by employing TSI refocusing pulses that have bandwidths which are less than the chemical shift difference between the target spins and the spins to which they are weakly coupled. The modified TSI sequence rewinds the J-evolution of the target protons in the slice of interest independently of the echo time or echo spacing, thereby removing any signal variation between successive echoes (apart from T₂ relaxation effects). In this study we tailored the narrow-bandwidth TSI sequence for observation of the C₂ Glx protons. The echo time was experimentally optimized to minimize signal contamination from myo-inositol, and the efficacy of the method was verified on phantom solutions of Glx and on brain in vivo.     

MSE-MRI sequence optimisation for measurement of bi- and tri-exponential T2 relaxation in a phantom and fruit

The transverse relaxation signal from vegetal cells can be described by multi-exponential behaviour, reflecting different water compartments. This multi-exponential relaxation is rarely measured by conventional MRI imaging protocols; mono-exponential relaxation times are measured instead, thus limiting information about of the microstructure and water status in vegetal cells. In this study, an optimised multiple spin echo (MSE) MRI sequence was evaluated for assessment of multi-exponential transverse relaxation in fruit tissues. The sequence was designed for the acquisition of a maximum of 512 echoes. Non-selective refocusing RF pulses were used in combination with balanced crusher gradients for elimination of spurious echoes. The study was performed on a bi-compartmental phantom with known T₂ values and on apple and tomato fruit. T₂ decays measured in the phantom and fruit were analysed using bi- and tri-exponential fits, respectively. The MRI results were compared with low field non-spatially resolved NMR measurements performed on the same samples.     

Microstructural analysis of foam by use of NMR R2 dispersion

The spin–spin relaxation rate R₂ (=1/T₂) in hydrogel foams measured by use of a multiple spin echo sequence is found to be dependent on the echo time spacing. This property, referred to as R₂-dispersion, originates to a large extent from molecular self-diffusion of water within internal field gradients that result from magnetic susceptibility differences between the gel and air phase. Another contribution to the R₂ relaxation rate is surface relaxation. Numerical simulations are performed to investigate the relation between the foam microstructure (the mean air bubble radius and standard deviation of the air bubble radius) and foam composition properties (such as magnetic susceptibilities, diffusion coefficient and surface relaxivity) at one hand and the R₂-dispersion at the other hand. The simulated R₂-dispersions of gel foam are in agreement with the measured R₂-dispersions. By correlating the R₂-dispersion parameters and simulated microstructure properties a semi-empirical relationship is obtained that enables the mean air bubble size to be derived from measured R₂-dispersion curves. The R₂-derived mean air bubble size of a hydrogel foam is in agreement with the bubble size measured with X-ray micro-CT. This illustrates the feasibility of using ¹H R₂-dispersion measurements to determine the size of air bubbles in hydrogel foams and of alveoli in lung tissue.     

Paramagnetic relaxation in sandstones: distinguishing T1 and T2 dependence on surface relaxation, internal gradients and dependence on echo spacing

This work provides a generalized theory of proton relaxation in inhomogeneous magnetic fields. Three asymptotic regimes of relaxation are identified depending on the shortest characteristic time scale. Numerical simulations illustrate that the relaxation characteristics in the regimes such as the T(1)/T(2) ratio and echo spacing dependence are determined by the time scales. The theoretical interpretation is validated for fluid relaxation in porous media in which field inhomogeneity is induced due to susceptibility contrast of fluids and paramagnetic sites on pore surfaces. From a set of measurements on model porous media, we conclude that when the sites are small enough, no dependence on echo spacing is observed with conventional low-field NMR spectrometers. Echo spacing dependence is observed when the paramagnetic materials become large enough or form a 'shell' around each grain such that the length scale of the region of induced magnetic gradients is large compared to the diffusion length during the time of the echo spacing. The theory can aid in interpretation of diffusion measurements in porous media as well as imaging experiments in presence of contrast agents used in MRI.     

Single-slice mapping of ultrashort T(2)

In this communication we present a method for single-slice mapping of ultrashort transverse relaxation times T(2). The RF pulse sequence consists of a spin echo preparation of the magnetization followed by slice-selective ultrashort echo time (UTE) imaging with radial k-space sampling. In order to keep the minimum echo time as small as possible, avoid out-of-slice contamination and signal contamination due to unwanted echoes, the implemented pulse sequence employs a slice-selective 180° RF refocusing pulse and a 4-step phase cycle. The slice overlap of the two slice-selective RF pulses was investigated. An acceptable Gaussian slice profile could be achieved by adjusting the strength of the two slice-selection gradients. The method was tested on a short T(2) phantom consisting of an arrangement of a roll of adhesive tape, an eraser, a piece of modeling dough made of Plasticine?, and a 10% w/w agar gel. The T(2) measurements on the phantom revealed exponential signal decays for all samples with T(2)(adhesive tape)=(0.5 ± 0.1)ms, T(2)(eraser)=(2.33 ± 0.07)ms, T(2)(Plasticine?)=(2.8 ± 0.06)ms, and T(2)(10%agar)=(9.5 ± 0.83)ms. The T(2) values obtained by the mapping method show good agreement with the T(2) values obtained by a non-selective T(2) measurement. For all samples, except the adhesive tape, the effective transverse relaxation time T(2)(?) was significantly shorter than T(2). Depending on the scanner hardware the presented method allows mapping of T(2) down to a few hundreds of microseconds. Besides investigating material samples, the presented method can be used to study the rapidly decaying MR-signal from biological tissue (e.g.: bone, cartilage, and tendon) and quadrupolar nuclei (e.g.: (23)Na, (35)Cl, and (17)O).     

Comparative study of motions in dimethylsulfone by noise excitation and solid echo spectroscopy

2H NMR spectra of dimethylsulfone were measured with noise excitation and solid echo NMR spectroscopy in the temperature range from 125 to 355 K. Besides the known fact that broad NMR spectra can be measured with both methods, in comparable times it is shown that for noise excitation, the signal loss is negligible compared to echo spectroscopy in the regime when the correlation times of the motions are of the order of magnitude of the echo pulse spacing. For simulating the dynamic NMR spectra acquired with noise excitation, only the motional process must be taken into account and relaxation can be neglected. Furthermore, the problem of restricted acquisition bandwidth in noise NMR spectroscopy is discussed.     

Bi-exponential proton transverse relaxation rate (R2) image analysis using RF field intensity-weighted spin density projection: potential for R2 measurement of iron-loaded liver

A bi-exponential proton transverse relaxation rate (R(2)) image analysis technique has been developed that enables the discrimination of dual compartment transverse relaxation behavior in systems with rapid transverse relaxation enhancement. The technique is particularly well suited to single spin-echo imaging studies where a limited number of images are available for analysis. The bi-exponential R(2) image analysis is facilitated by estimation of the initial proton spin density signal within the region of interest weighted by the RF field intensities. The RF field intensity-weighted spin density map is computed by solving a boundary value problem presented by a high spin density, long T(2) material encompassing the region for analysis. The accuracy of the bi-exponential R(2) image analysis technique is demonstrated on a simulated dual compartment manganese chloride phantom system with relaxation rates and relative population densities between the two compartments similar to the bi-exponential transverse relaxation behavior expected of iron loaded liver. Results from analysis of the phantoms illustrate the potential of bi-exponential R(2) image analysis with RF field intensity-weighted spin density projection for quantifying transverse relaxation enhancement as it occurs in liver iron overload.     

Measurement of in vivo multi-component T2 relaxation times for brain tissue using multi-slice T2 prep at 1.5 and 3 T

The objective of this study was to implement a clinically relevant multi-slice multi-echo imaging sequence in order to quantify multi-component T2 relaxation times for normal volunteers at both 1.5 and 3 T. Multi-echo data were fitted using a nonnegative least square algorithm. Twelve echo data with nonlinear echo sampling were acquired using a receive-only eight-channel phased array coil and volume head coil for phantoms and normal volunteers, and compared to 32-echo data with linear echo sampling. It was observed that the performance of the 180 degrees refocusing trains was more spatially uniform for the receive-only eight-channel phased array coil than for the head coil, particularly at 3 T. The phantom study showed that the estimated T2 relaxation times were accurate and reproducible for both single- and multi-slice acquisition from a commercial phantom with known T2 relaxation times. Short T2 components (T2 <50 ms) were mainly observed within the white matter for normal volunteers, and the fraction of short T2 water components (i.e., myelin water) was 7-12% of total water. It was observed that the calculated myelin water fraction map from the nonlinearly sampled 12-echo data was comparable with that from the linearly sampled 32-echo data. Quantification of T2 relaxation times from multi-slice images was accomplished with a clinically acceptable scan times (16 min) for normal volunteers by using a nonselective T2 prep imaging sequence. The use of the eight-channel head coil involved more accurate quantification of T2 relaxation times particularly when the number of echoes was limited.     

Evaluation of the echolocation model for range estimation of multiple closely spaced objects

Experimental evidence indicates that bats can use frequency-modulated echolocation to identify objects with an accuracy of less than 1 μs. However, when modeling this process, it is difficult to estimate the delay times of multiple closely spaced objects by analyzing the echo spectrum, because the sequence of delay separations cannot be determined without information on the temporal changes in the interference patterns of the echoes. To extract the temporal changes, Gaussian chirplets with a carrier frequency compatible with bat emission sweep rates are introduced. The delay time for object 1 (T(1)) is estimated from the echo spectrum around the onset time. The T(2) is obtained by adding the T(1) to the delay separation between objects 1 and 2. Further objects are located in sequence by this procedure. Here echoes were measured from single and multiple objects at a low signal-to-noise ratio. It was confirmed that the delay time for a single object could be estimated with an accuracy of about 1.3 μs. The range accuracy was less than 6 μs when the frequency bandwidth was less than 10 kHz. The delay time for multiple closely spaced objects could be estimated with a high range resolution by extracting the interference pattern.     

Relaxation effects in the quantification of fat using gradient echo imaging

Quantification of fat has been investigated using images acquired from multiple gradient echoes. The evolution of the signal with echo time and flip angle was measured in phantoms of known fat and water composition and in 21 research subjects with fatty liver. Data were compared to different models of the signal equation, in which each model makes different assumptions about the T1 and/or T2* relaxation effects. A range of T1, T2*, fat fraction and number of echoes was investigated to cover situations of relevance to clinical imaging. Results indicate that quantification is most accurate at low flip angles (to minimize T1 effects) with a small number of echoes (to minimize spectral broadening effects). At short echo times, the spectral broadening effects manifest as a short apparent T2 for the fat component.     

The effect of varying echo spacing within a multiecho acquisition: better characterization of long T2 components

A 48-echo pulse sequence with five different echo-spacing combinations was examined to determine how one can most effectively measure the T2 relaxation characteristics of cerebral tissue containing a long T2 component. For each scan, the first 32 echoes had an echo spacing of 10 ms, while the spacing for Echoes 33-48 (DeltaTE2) was 10, 20, 30, 40 or 50 ms. In an in vivo study using 10 normal volunteers, it was found that the resolution of T2 distribution peaks for both myelin water (approximately 20 ms) and intracellular/extracellular (IE) water (approximately 80 ms) improved as DeltaTE2 increased. The geometric mean T2 values of the main peak agreed within the error for all DeltaTE2 values. A phantom study simulated T2 relaxation distributions that are expected in the brains of patients with demyelinating diseases. For phantoms in which the T2 values of the IE and lesion (200-500 ms) water compartments were separated by at least a factor of 3, each compartment in the distribution was better resolved when DeltaTE2=40 or 50 ms. On the basis of these results, we recommend the use of extended DeltaTE2 values for imaging patients with lesions, without the risk of losing valuable short T2 information.     

High-speed spectroscopic imaging for cancellous bone marrow R(2)* mapping and lipid quantification

In this work an interleaved multiple-gradient-echo chemical shift imaging (IMGE-CSI) technique was designed, implemented and evaluated at 1.5 and 4T for high-resolution lipid quantification and R(2)* measurement in-vivo. The method is analogous to echo planar CSI but utilizes conventional gradient echoes, exploiting the principle of spectroscopic bandwidth extension by interleaving temporally offset gradient-echo trains. It is shown that IMGE-CSI is able to measure true fat volume fraction in oil/water mixtures with high accuracy, not possible with Dixon-type methods which approximate the spectrum as consisting of only two spectral components. Correlation of the CSI- derived volume fractions with volumetry afforded r(2) > 0.99 with a slope of 0.98. The method is shown to be able to quantify regional variations in bone marrow composition in vivo with a spatial resolution of 2.5 x 2.5 x 5 mm(3.) R(2)* was obtained by multi-line spectral curve fitting. For the measurement of R(2)* in cancellous bone marrow the method is shown to agree well with time-domain fitting techniques but is superior in instances where the marrow has both hematopoietic and fatty constituents. Finally, excellent inter-scan reproducibility (1% coefficient of variation for global means and medians) was achieved, yielding r(2) = 0.98 of the test-retest correlation for three scans in four test subjects. In conclusion, IMGE-CSI is found to enable highly accurate lipid quantification and measurement of cancellous bone marrow R(2)* at spatial resolutions and scan times typical of standard clinical protocols.     

Fast T2-weighted MR imaging: impact of variation in pulse sequence parameters on image quality and artifacts

The purpose of this study was to quantitatively evaluate in a phantom model the practical impact of alteration of key imaging parameters on image quality and artifacts for the most commonly used fast T(2)-weighted MR sequences. These include fast spin-echo (FSE), single shot fast spin-echo (SSFSE), and spin-echo echo-planar imaging (EPI) pulse sequences. We developed a composite phantom with different T1 and T2 values, which was evaluated while stationary as well as during periodic motion. Experiments involved controlled variations in key parameters including effective TE, TR, echo spacing (ESP), receive bandwidth (BW), echo train length (ETL), and shot number (SN). Quantitative analysis consisted of signal-to-noise ratio (SNR), image nonuniformity, full-width-at-half-maximum (i.e., blurring or geometric distortion) and ghosting ratio. Among the fast T(2)-weighted sequences, EPI was most sensitive to alterations in imaging parameters. Among imaging parameters that we tested, effective TE, ETL, and shot number most prominently affected image quality and artifacts. Short T(2) objects were more sensitive to alterations in imaging parameters in terms of image quality and artifacts. Optimal clinical application of these fast T(2)-weighted imaging pulse sequences requires careful attention to selection of imaging parameters.     

Identification of patients with hereditary haemochromatosis by magnetic resonance imaging and spectroscopic relaxation time measurements

A total of 4302 healthy blood donors were screened for elevated serum ferritin and transferrin saturation. Fifteen had increased serum ferritin at a follow-up examination. Five relatives of these donors also entered the study. Eleven patients had elevated liver iron concentrations, while five had normal liver iron concentrations. The R₂ relaxation rate in the liver was first measured with a conventional multi-spin-echo imaging sequence, and then by a volume-selective spectroscopic multi-spin-echo sequence, in order to achieve a minimum echo time of 4 msec. No correlation was found between the relaxation rate R₂ and the liver iron concentration, when R₂ was calculated from the imaging data. Multi-exponential transverse relaxation could be resolved when the spectroscopic sequence was used. A strong correlation between the initial slope of the relaxation curve and the liver iron concentration was found (r = 0.90, p < 0.001). Signal intensity ratios between liver and muscle were calculated from the first three echoes in the multi-echo imaging sequence, and from a gradient echo sequence. A strong correlation between the logarithm of the signal intensity ratios and the liver iron concentration was found. Although both spectroscopic T₂ relaxation time measurements and signal intensity ratios could be used to quantify liver iron concentration, the gradient echo imaging seemed to be the best choice. Gradient echo imaging could be performed during a single breath hold, so motion artifacts could be avoided. The accuracy of liver iron concentration estimates from signal intensity ratios in the gradient echo images was about 35%.     

Optimization of timing in the Carr-Purcell-Meiboom-Gill sequence.

The Carr-Purcell-Meiboom-Gill sequence is widely used in grossly inhomogeneous fields to characterize fluid saturated porous media. The distribution of T(2) relaxation times is a measure of the distribution of pore sizes. We present here a theoretical analysis of the spin dynamics of this sequence in strongly inhomogeneous fields. Based on this analysis, we optimize the timing of the sequence for maximal signal bandwidth. It is shown that the initial pulse spacing should be decreased by a time 2t(90)/pi. Experimental results in a strayfield set-up confirm the theoretical analysis. The optimized timing increases the measured signal bandwidth and increases the ratio of signal-to-noise by 1.2 dB without affecting the measured relaxation time.     

An investigation of the origins of contrast in NMR spin echo images of plant tissue.

The effects that the spatial distribution of water protons and their transverse relaxation times have on the image contrast of spin echo images of courgette was investigated. The T2-weighted image of courgette contains the most anatomical information. The image contrast was explained using a phenomenological theory based on the Bloch equations, which gave an insight into the morphology and microdynamics of water in the plant tissue. The perceived contrast in the spin echo images of courgette, glucose and Sephadex bead solutions can be dramatically altered by keeping all the imaging acquisition parameters constant, such as the recycle and echo time, but reducing the interpulse spacing by introducing a CPMG train of 180 degrees pulses into the middle of the sequence. These changes were interpreted by considering the microenvironment of the water. This work demonstrates that the origin of image contrast in T2-weighted images of plant tissue can be understood using the water proton transverse relaxation theory developed by Hills et al.     

Parametric multiecho proton spectroscopic imaging: Application to the rat brain in vivo

A parametric multiecho variant of proton spectroscopic imaging (SI) is presented using a multiecho SI sequence with uniform phase-encoding of all echoes within each echo train. The acquisition of SI data sets at different echo times (TE) increases the amount of information obtained within the same total measuring time as in standard SI measurements. The gain in information can be used: (a) to choose the most appropriate TE for each metabolite signal with respect to T₂, spin coupling, or problems caused by peak overlap; (b) to measure the relaxation time T₂ of metabolite signals with high spatial resolution; or (c) to improve the signal-to-noise ratio for metabolite signals with long T₂ values by adding spectra calculated from consecutive echoes. The method was tested in vivo on healthy rat brain and applied to study metabolic changes in rat brain lesions.     

Nuclear magnetic resonance in inhomogeneous magnetic fields

The response of the spin system has been investigated by numerical simulations in the case of a nuclear magnetic resonance (NMR) experiment performed in inhomogeneous static and radiofrequency fields. The particular case of the NMR-MOUSE was considered. The static field and the component of the radiofrequency field perpendicular to the static field were evaluated as well as the spatial distribution of the maximum NMR signal detected by the surface coil. The NMR response to various pulse sequences was evaluated numerically for the case of an ensemble of isolated spins (1/2). The behavior of the echo train in Carr-Purcell-like pulse sequences used for measurements of transverse relaxation and self-diffusion was simulated and compared with the experiment. The echo train is shown to behave qualitatively differently depending on the particular phase schemes used in these pulse sequences. Different echo trains are obtained, because of the different superposition of Hahn and stimulated echoes forming mixed echoes as a result of the spatial distribution of pulse flip angles. The superposition of Hahn and stimulated echoes originating from different spatial regions leads to distortions of the mixed echoes in intensity, shape, and phase. The volume selection produced by Carr-Purcell-like pulse sequences is also investigated for the NMR-MOUSE. The developed numerical simulation procedure is useful for understanding a variety of experiments performed with the NMR-MOUSE and for improving its performance. Copyright 2000 Academic Press.