On the lorentzian versus Gaussian character of time-domain spin-echo signals from the brain as sampled by means of gradient-echoes: Implications for quantitative transverse relaxation studies

Myelin water imaging at 0.55 T using a multigradient-echo sequence

PURPOSE

To investigate the prospects of a multigradient-echo (mGRE) acquisition for in vivo myelin water imaging at 0.55 T.

METHODS

Scans were performed on the brain of four healthy volunteers at 0.55 and 3 T, using a 3D mGRE sequence. The myelin water fraction (MWF) was calculated for both field strengths using a nonnegative least squares (NNLS) algorithm, implemented in the qMRLab suite. The quality of these maps as well as single-voxel fits were compared visually for 0.55 and 3 T.

RESULTS

The obtained MWF values at 0.55 T are consistent with previously reported ones at higher field strengths. The MWF maps are a considerable improvement over the ones at 3 T. Example fits show that 0.55 T data is better described by an exponential model than 3 T data, making the assumed multi-exponential model of the NNLS algorithm more accurate.

CONCLUSION

This first assessment shows that mGRE myelin water imaging at 0.55 T is feasible and has the potential to yield better results than at higher fields.

Information theoretic evaluation of Lorentzian, Gaussian, Voigt, and symmetric alpha-stable models of reversible transverse relaxation in cervical cancer in vivo at 3 T

OBJECTS

To better characterize cervical cancer at 3 T. MRI transverse relaxation patterns hold valuable biophysical information about cellular scale microstructure. Lorentzian modeling is typically used to represent intravoxel frequency distributions, resulting in mono-exponential decay of reversible transverse relaxation. However, deviations from mono-exponential decay are expected theoretically and observed experimentally.

MATERIALS AND METHODS

We compared the information content of four models of signal attenuation with reversible transverse relaxation. Biological phantoms and six women with cervical squamous cell carcinoma were imaged using a gradient-echo sampling of the spin-echo (GESSE) sequence. Lorentzian, Gaussian, Voigt, and Symmetric α-Stable (SAS) models were ranked using Akaike's Information Criterion (AIC), and the model retaining the highest information content was identified at each voxel as the best model.

RESULTS

The Lorentzian model resulted in information loss in large fractions of the phantoms and cervix. Gaussian and SAS models frequently had higher information content than the Lorentzian in much of the areas of interest. The Voigt model rarely surpassed the three other models in terms of information content.

DISCUSSION

Gaussian and SAS models provide better fitting of data in much of the human cervix at 3 T. Minimizing information loss through improved tissue modeling may have important implications for identifying reliable biomarkers of tumor hypoxia and iron deposition.

A relaxometry method that emphasizes practicality and availability

PURPOSE

Although many methods have been proposed to quantitatively map the main MRI parameters (e.g., T₁ , T₂ , C × M₀ ), these methods often involve special sequences not readily available on clinical scanners and/or may require long scan times. In contrast, the proposed method can readily run on most scanners, offer flexible tradeoffs between scan time and image quality, and map MRI parameters jointly to ensure spatial alignment.

METHODS

The approach is based on the multi-shot spin-echo (SE) EPI sequence. The corresponding signal equation was derived and strategies for solving it were developed. As usual with multi-shot EPI, scan time can readily be traded-off against image quality by adjusting the echo train length. Validation was performed against reference relaxometry methods, in gel phantoms with varying concentrations of gadobutrol and gadoterate meglumine contrast agents. In vivo examples are further presented, from 3 neuroradiology patients.

RESULTS

Bland-Altman analysis was performed: for T₂ , as compared to 2D SE, bias was 0.29 ms and the 95% limits of agreement ranged from -1.15 to +1.73 ms. For T₁ , compared to inversion-recovery SE (and MOLLI), bias was -20.2 ms (and -14.5 ms) and the limits of agreement ranged from -62.4 to +22.0 ms (and -53.8 to +24.9 ms). The mean relative T₁ error between the proposed method and each of the 2 reference methods was similar to that of the reference methods among themselves.

CONCLUSION

In the constellation of existing relaxometry methods, the proposed method is meant to stand out in terms of its practicality and availability.

In vivo irreversible and reversible transverse relaxation rates in human cerebral cortex via line scans at 7 T with 250 micron resolution perpendicular to the cortical surface

Understanding how and why MR signals and their associated relaxation rates vary with cortical depth could ultimately enable the noninvasive investigation of the laminar architecture of cerebral cortex in the living human brain. However, cortical gray matter is typically only a few millimeters thick, making it challenging to sample many cortical depths with the voxel sizes commonly used in MRI studies. Line-scan techniques provide a way to overcome this challenge and here we implemented a novel line-scan GESSE pulse sequence that allowed us to measure irreversible and reversible transverse relaxation rates-R₂ and R₂´, respectively-with extremely high resolution (250 μm) in the radial direction, perpendicular to the cortical surface. Eight healthy human subjects were scanned at 7 T using this sequence, with primary visual cortex (V1) targeted in three subjects and primary motor (M1) and somatosensory cortex (S1) targeted in the other five. In all three cortical areas, a peak in R₂ values near the central depths was seen consistently across subjects-an observation that has not been made before, to our knowledge. On the other hand, no consistent pattern was apparent for R₂´ values as a function of cortical depth. The intracortical R₂ peak reported here is unlikely to be explained by myelin content or by deoxyhemoglobin in the microvasculature; however, this peak is in accord with the laminar distribution of non-heme iron in these cortical areas, known from prior histology studies. Obtaining information about tissue microstructure via measurements of transverse relaxation (and other quantitative MR contrast mechanisms) at the extremely high radial resolutions achievable through the use of line-scan techniques could therefore bring us closer to being able to perform "in vivo histology" of the cerebral cortex.

Synthetic MRI: Technologies and Applications in Neuroradiology

Synthetic MRI is a technique that synthesizes contrast-weighted images from multicontrast MRI data. There have been advances in synthetic MRI since the technique was introduced. Although a number of synthetic MRI methods have been developed for quantifying one or more relaxometric parameters and for generating multiple contrast-weighted images, this review focuses on several methods that quantify all three relaxometric parameters (T₁ , T₂ , and proton density) and produce multiple contrast-weighted images. Acquisition, quantification, and image synthesis techniques are discussed for each method. We discuss the image quality and diagnostic accuracy of synthetic MRI methods and their clinical applications in neuroradiology. Based on this analysis, we highlight areas that need to be addressed for synthetic MRI to be widely implemented in the clinic. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 1.

A new concept for improved quantitative analysis of reversible transverse relaxation in tissues with variable microscopic field distribution

PURPOSE

The intravoxel distribution of the magnetic field strongly influences signal dephasing after RF excitation and the resulting signal decay in gradient echo-based MRI. In this work, several different field distribution models were applied and tested for analysis of microscopic field characteristics within pixels.

THEORY

A flexible model for improved pixel-wise characterization of the underlying field distribution is introduced. The proposed symmetric alpha-stable (SαS) distribution covers Lorentzian, Gaussian, and intermediate field distributions in a continuous way using a two-parametric (width and shape) function.

METHODS

The new model was applied on human brain, potatoes (homogeneous isotropic tissue), and stems of pineapple (anisotropic fibrous tissue). Effects of microscopic structure and background gradients on the shape and the widths of the microscopic field distribution were analyzed using gradient echo sampling of the spin echo and multigradient-echo sequences. Effects of non-Lorentzian shapes of microscopic field distributions on the results of common T 2 ∗ measurements with mono-exponential fitting of signal values were tested.

RESULTS

Many pixels of the examined objects showed field characteristics in between Lorentzian and Gaussian shapes. Microscopic field inhomogeneities caused by microscopic susceptibility effects and background gradients sometimes led to rather Gaussian than Lorentzian field distribution. In cases with nearly Gaussian field distribution, mono-exponential fitting of the signal decay resulted in different T 2 ∗ values, depending on the sampling points.

CONCLUSIONS

Using the concept of more flexible distributions for characterization of microscopic susceptibility effects in tissue provides better fitting of data and nearly sampling point-independent results than common T 2 ∗ measurements with mono-exponential fitting.

Spectrally-selective measurements of reversible and irreversible transverse relaxation rates from single spin-echo PRESS acquisitions in muscle

The goal of this study was to test a new formalism for extracting reversible and irreversible transverse relaxation rates from resonances within typical proton muscle spectra using only a single spin echo as acquired with routine single-voxel, point-resolved echo spectroscopy (PRESS) acquisitions. Single-voxel, non-water-suppressed PRESS acquisitions within the calf muscles of four healthy subjects were performed at 1.5 T using six echo times ranging from 30 to 576 ms. Novel transverse relaxation analyses of water, choline, creatine, and lipid resonances were performed based upon the disparate relaxation sensitivities of the left versus the right sides of spectroscopically sampled spin echoes. Irreversible and reversible transverse relaxation rates R₂ and R₂ ' were extracted for water, metabolites, and lipids using echo times of 288 ms and longer. The R₂ values so obtained were compared with more conventional "gold standard" Hahn values, R_2Hahn , evaluated from the echo-time dependence of spectral peak areas generated from right-side sampling alone. Water resonances displayed biexponential Hahn signal decays, consistent with observations of decreasing R₂ values with increasing echo time via the new approach. Choline and creatine resonances displayed monoexponential echo-time decays, with R_2Hahn values in reasonable agreement with R₂ values obtained from the single-echo analyses at the longer echo times. Lipid methylene and methyl R₂ values extracted from the new approach were also in reasonable accord with R_2Hahn values. Further validation of the technique was provided through PRESS acquisitions on a water phantom to which various levels of gadolinium were added in order to manipulate transverse relaxation rates, yielding excellent agreement between water-resonance R_2Hahn and single-echo R₂ values. In summary, this work demonstrates the feasibility of measuring reversible and irreversible transverse relaxation rates for individual spectral peaks from single-echo PRESS acquisitions, enabling a reduction in overall scan time relative to the use of multiple acquisitions with varying echo time.

In vivo measurements of irreversible and reversible transverse relaxation rates in human basal ganglia at 7 T: making inferences about the microscopic and mesoscopic structure of iron and calcification deposits

The goal of this study was to measure irreversible and reversible transverse relaxation rates in the globus pallidus and putamen at 7 T, and to use these rates to make inferences about the sub-voxel structure of iron and calcification deposits. Gradient Echo Sampling of a Spin Echo (GESSE) data were acquired at 7 T on eighteen volunteers spanning a large range of ages (23-85 years), with calcifications in the globus pallidus incidentally observed in one volunteer. Maps of transverse relaxation rates were derived from the GESSE data, and the mean value of these rates in globus pallidus and putamen was estimated for each volunteer. Both irreversible and reversible transverse relaxation rates increased with the expected age-dependent iron content in these structures, except for the individual with calcifications for whom extremely large reversible relaxation rates but normal irreversible relaxation rates were found in the globus pallidus. Given the sensitivity of irreversible and reversible transverse relaxation rates to microscopic and mesoscopic field variations, respectively, our findings suggest that joint consideration of these rates may yield information not only about the amount of iron and calcification deposited in the brain, but also about the sub-voxel structure of these deposits, perhaps revealing certain aspects of their geometry and cellular distribution.

Magnetic resonance field fingerprinting

PURPOSE

To develop and evaluate the magnetic resonance field fingerprinting method that simultaneously generates T₁ , T₂ , B₀ , and B 1 + maps from a single continuous measurement.

METHODS

An encoding pattern was designed to integrate true fast imaging with steady-state precession (TrueFISP), fast imaging with steady-state precession (FISP), and fast low-angle shot (FLASH) sequence segments with varying flip angles, radio frequency (RF) phases, TEs, and gradient moments in a continuous acquisition. A multistep matching process was introduced that includes steps for integrated spiral deblurring and the correction of intravoxel phase dispersion. The method was evaluated in phantoms as well as in vivo studies in brain and lower abdomen.

RESULTS

Simultaneous measurement of T₁ , T₂ , B₀ , and B 1 + is achieved with T₁ and T₂ subsequently being less afflicted by B₀ and B 1 + variations. Phantom results demonstrate the stability of generated parameter maps. Higher undersampling factors and spatial resolution can be achieved with the proposed method as compared with solely FISP-based magnetic resonance fingerprinting. High-resolution B₀ maps can potentially be further used as diagnostic information.

CONCLUSION

The proposed magnetic resonance field fingerprinting method can estimate T₁ , T₂ , B₀ , and B 1 + maps accurately in phantoms, in the brain, and in the lower abdomen.

Design of spectral-spatial phase prewinding pulses and their use in small-tip fast recovery steady-state imaging

PURPOSE

Spectrally selective "prewinding" radiofrequency pulses can counteract B0 inhomogeneity in steady-state sequences, but can only prephase a limited range of off-resonance. We propose spectral-spatial small-tip angle prewinding pulses that increase the off-resonance bandwidth that can be successfully prephased by incorporating spatially tailored excitation patterns.

THEORY AND METHODS

We present a feasibility study to compare spectral and spectral-spatial prewinding pulses. These pulses add a prephasing term to the target magnetization pattern that aims to recover an assigned off-resonance bandwidth at the echo time. For spectral-spatial pulses, the design bandwidth is centered at the off-resonance frequency for each spatial location in a field map. We use these pulses in the small-tip fast recovery steady-state sequence, which is similar to balanced steady-state free precession. We investigate improvement of spectral-spatial pulses over spectral pulses using simulations and small-tip fast recovery scans of a gel phantom and human brain.

RESULTS

In simulation, spectral-spatial pulses yielded lower normalized root mean squared excitation error than spectral pulses. For both experiments, the spectral-spatial pulse images are also qualitatively better (more uniform, less signal loss) than the spectral pulse images.

CONCLUSION

Spectral-spatial prewinding pulses can prephase over a larger range of off-resonance than their purely spectral counterparts. Magn Reson Med 79:1377-1386, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Visualization of calcium phosphate cement in teeth by zero echo time 1 H MRI at high field

¹ H magnetic resonance imaging (MRI) by a zero echo time (ZTE) sequence is an excellent method to image teeth. Calcium phosphate cement (CPC) materials are applied in the restoration of tooth lesions, but it has not yet been investigated whether they can be detected by computed tomography (CT) or MRI. The aim of this study was to optimize high-field ZTE imaging to enable the visualization of a new CPC formulation implanted in teeth and to apply this in the assessment of its decomposition in vivo. CPC was implanted in three human and three goat teeth ex vivo and in three goat teeth in vivo. An ultrashort echo time (UTE) sequence with multiple flip angles and echo times was applied at 11.7 T to measure T₁ and T₂ * values of CPC, enamel and dentin. Teeth with CPC were imaged with an optimized ZTE sequence. Goat teeth implanted with CPC in vivo were imaged after 7 weeks ex vivo. T₂ * relaxation of implanted CPC, dentin and enamel was better fitted by a model assuming a Gaussian rather than a Lorentzian distribution. For CPC and human enamel and dentin, the average T₂ * values were 273 ± 19, 562 ± 221 and 476 ± 147 μs, respectively, the average T₂ values were 1234 ± 27, 963 ± 151 and 577 ± 41 μs, respectively, and the average T₁ values were 1065 ± 45, 972 ± 40 and 903 ± 7 ms, respectively. In ZTE images, CPC had a higher signal-to-noise-ratio than dentin and enamel because of the higher water content. Seven weeks after in vivo implantation, the CPC-filled lesions showed less homogeneous structures, a lower T₁ value and T₂ * separated into two components. MRI by ZTE provides excellent contrast for CPC in teeth and allows its decomposition to be followed.

Characterization of gradient echo signal decays in healthy and cancerous prostate at 3T improves with a Gaussian augmentation of the mono-exponential (GAME) model

A biomarker of cancer aggressiveness, such as hypoxia, could substantially impact treatment decisions in the prostate, especially radiation therapy, by balancing treatment morbidity (urinary incontinence, erectile dysfunction, etc.) against mortality. R2 (*) mapping with Mono-Exponential (ME) decay modeling has shown potential for identifying areas of prostate cancer hypoxia at 1.5T. However, Gaussian deviations from ME decay have been observed in other tissues at 3T. The purpose of this study is to assess whether gradient-echo signal decays are better characterized by a standard ME decay model, or a Gaussian Augmentation of the Mono-Exponential (GAME) decay model, in the prostate at 3T. Multi-gradient-echo signals were acquired on 20 consecutive patients with a clinical suspicion of prostate cancer undergoing MR-guided prostate biopsies. Data were fitted with both ME and GAME models. The information contents of these models were compared using Akaike's information criterion (second order, AICC ), in skeletal muscle, the prostate central gland (CG), and peripheral zone (PZ) regions of interest (ROIs). The GAME model had higher information content in 30% of the prostate on average (across all patients and ROIs), covering up to 67% of cancerous PZ ROIs, and up to 100% of cancerous CG ROIs (in individual patients). The higher information content of GAME became more prominent in regions that would be assumed hypoxic using ME alone, reaching 50% of the PZ and 70% of the CG as ME R2 (*) approached 40 s(-1) . R2 (*) mapping may have important applications in MRI; however, information lost due to modeling could mask differences in parameters due to underlying tissue anatomy or physiology. The GAME model improves characterization of signal behavior in the prostate at 3T, and may increase the potential for determining correlates of fit parameters with biomarkers, for example of oxygenation status.

Bone marrow segmentation based on a combined consideration of transverse relaxation processes and Dixon oscillations

The aim of this study was to demonstrate that gradient-echo sampling of single spin echoes can be used to isolate the signal from trabecular bone marrow, with high-quality segmentation and surface reconstructions resulting from the application of simple post-processing strategies. Theoretical expressions of the time-domain single-spin-echo signal were used to simulate signals from bone marrow, non-bone fatty deposits and muscle. These simulations were compared with and used to interpret signals obtained by the application of the gradient-echo sampling of a spin-echo sequence to image the knee and surrounding tissues at 1.5 T. Trabecular bone marrow has a much higher reversible transverse relaxation rate than surrounding non-bone fatty deposits and other musculoskeletal tissues. This observation, combined with a choice of gradient-echo spacing that accentuates Dixon-type oscillations from chemical-shift interference effects, enabled the isolation of bone marrow signal from surrounding tissues through the use of simple image subtraction and thresholding. Three-dimensional renderings of the marrow surface were then readily generated with this approach - renderings that may prove useful for bone morphology assessment, e.g. for the measurement of femoral anteversion. In conclusion, understanding the behavior of signals from bone marrow and surrounding tissue as a function of time through a spin echo facilitates the segmentation and reconstruction of bone marrow surfaces using straightforward post-processing strategies that are typically available on modern radiology workstations.

Dual-pathway multi-echo sequence for simultaneous frequency and T2 mapping

PURPOSE

To present a dual-pathway multi-echo steady state sequence and reconstruction algorithm to capture T2, T2(∗) and field map information.

METHODS

Typically, pulse sequences based on spin echoes are needed for T2 mapping while gradient echoes are needed for field mapping, making it difficult to jointly acquire both types of information. A dual-pathway multi-echo pulse sequence is employed here to generate T2 and field maps from the same acquired data. The approach might be used, for example, to obtain both thermometry and tissue damage information during thermal therapies, or susceptibility and T2 information from a same head scan, or to generate bonus T2 maps during a knee scan.

RESULTS

Quantitative T2, T2(∗) and field maps were generated in gel phantoms, ex vivo bovine muscle, and twelve volunteers. T2 results were validated against a spin-echo reference standard: A linear regression based on ROI analysis in phantoms provided close agreement (slope/R(2)=0.99/0.998). A pixel-wise in vivo Bland-Altman analysis of R2=1/T2 showed a bias of 0.034 Hz (about 0.3%), as averaged over four volunteers. Ex vivo results, with and without motion, suggested that tissue damage detection based on T2 rather than temperature-dose measurements might prove more robust to motion.

CONCLUSION

T2, T2(∗) and field maps were obtained simultaneously, from the same datasets, in thermometry, susceptibility-weighted imaging and knee-imaging contexts.

Characterizing gradient echo signal decays in gynecologic cancers at 3T using a Gaussian augmentation of the monoexponential (GAME) model

PURPOSE

To assess whether R2* mapping with a standard Monoexponential (ME) or a Gaussian Augmentation of the Monoexponential (GAME) decay model better characterizes gradient-echo signal decays in gynecological cancers after external beam radiation therapy at 3T, and evaluate implications of modeling for noninvasive identification of intratumoral hypoxia.

MATERIALS AND METHODS

Multi-gradient-echo signals were acquired on 25 consecutive patients with gynecologic cancers and three healthy participants during inhalation of different oxygen concentrations at 3T. Data were fitted with both ME and GAME models. Models were compared using F-tests in tumors and muscles in patients, muscles, cervix, and uterus in healthy participants, and across oxygenation levels.

RESULTS

GAME significantly improved fitting over ME (P < 0.05): Improvements with GAME covered 34% of tumor regions-of-interest on average, ranging from 6% (of a vaginal tumor) to 68% (of a cervical tumor) in individual tumors. Improvements with GAME were more prominent in areas that would be assumed hypoxic based on ME alone, reaching 90% as ME R2* approached 100 Hz. Gradient echo decay parameters at different oxygenation levels were not significantly different (P = 0.81).

CONCLUSION

R2* may prove sensitive to hypoxia; however, inaccurate representations of underlying data may limit the success of quantitative assessments. Although the degree to which R2 or σ values correlate with hypoxia remains unknown, improved characterization with GAME increases the potential for determining any correlates of fit parameters with biomarkers, such as oxygenation status. J. MAGN. RESON. IMAGING 2016;44:1020-1030.