Magnetic field dependence of the distribution of NMR relaxation times in the living human brain

T2* quantification using multi-echo gradient echo sequences: a comparative study of different readout gradients

To quantify T2*, multiple echoes are typically acquired with a multi-echo gradient echo sequence using either monopolar or bipolar readout gradients. The use of bipolar readout gradients achieves a shorter echo spacing time, enabling the acquisition of a larger number of echoes in the same scan time. However, despite their relative time efficiency and the potential for more accurate quantification, a comparative investigation of these readout gradients has not yet been addressed. This work aims to compare the performance of monopolar and bipolar readout gradients for T2* quantification. The differences in readout gradients were theoretically investigated with a Cramér-Rao lower bound and validated with computer simulations with respect to the various imaging parameters (e.g., flip angle, TR, TE, TE range, and BW). The readout gradients were then compared at 3 T using phantom and in vivo experiments. The bipolar readout gradients provided higher precision than monopolar readout gradients in both computer simulations and experimental results. The difference between the two readout gradients increased for a lower SNR and smaller TE range, consistent with the prediction made using Cramér-Rao lower bound. The use of bipolar readout gradients is advantageous for regions or situations where a lower SNR is expected or a shorter acquisition time is required.

Whole-Brain Water Content Mapping Using Super-Resolution Reconstruction with MRI Acquisition in 3 Orthogonal Orientations

PURPOSE

Brain water content provides rich tissue contrast comparable to that of longitudinal relaxation time T₁ , but mapping is usually performed at modest resolution. In particular, the slice thickness in 2D mapping methods is limited. Here, we combine super-resolution reconstruction techniques with a fast water content mapping method to acquire high and isotropic resolution (0.75 mm) water content maps at 3 Tesla.

METHODS

A high-resolution multi-echo gradient echo image is super-resolution-reconstructed from 3 low-resolution, orthogonal multi-echo gradient echo image acquisitions, followed by water content mapping. The mapping accuracy and SNR of the proposed method are assessed using numerical simulations, phantom studies, and in vivo data acquired from 6 healthy volunteers at 3 Tesla. A high-resolution acquisition with an established mapping method is used as a reference.

RESULTS

Whole-brain water content maps with 0.75 mm isotropic resolution are demonstrated. No bias in the water content values was seen following super-resolution reconstruction. In the in vivo experiments, a lower SD of the mean water content values was observed with the proposed method compared to the reference method.

CONCLUSIONS

Super-resolution reconstruction of multi-echo gradient echo data is demonstrated, enabling whole-brain water content mapping with high and isotropic resolution. The accuracy of the proposed method is shown using phantoms and 6 healthy volunteers and was found to be unchanged compared to the conventional acquisition. The proposed method could increase the sensitivity of water content mapping sufficiently to enable the detection of very small lesions, such as cortical lesions in multiple sclerosis.

Development of a novel 10-echo multi-contrast sequence based on EPIK to deliver simultaneous quantification of T2 and T2 * with application to oxygen extraction fraction

PURPOSE

The simultaneous quantification of T₂ and T₂ * maps based on fast sequences for combined GE and SE acquisition has rekindled research and clinical interest by offering a wide range of attractive applications, e.g., dynamic tracking of oxygen extraction fraction (OEF). However, previously published methods based on EPI-readouts have been hindered by resolution and the number of acquired echoes.

METHODS

This work presents a novel 10-echo GE-SE EPIK (EPI with keyhole) sequence for the rapid quantification of T₂ '. T₂ /T₂ * maps from the GE-SE EPIK sequence were validated using three phantoms and 15 volunteers at 3T. The incorporation of a sliding window approach, combined with the full sampling of the k-space center inherent to EPIK, enables a high effective temporal resolution. That is, for an eight-slice breath-hold experiment, a temporal sampling rate of eight reconstructed slices per 1.1 s.

RESULTS

In comparison with repeated single-echo SE, multi-echo GE, and spectroscopy methods, the GE-SE EPIK sequence shows good agreement in quantifying T₂ /T₂ * values, while the gray matter/white matter separation yielded the expected contrast differentiation. The OEF was calculated with a view to an initial application with clinical relevance, producing results comparable to those in the literature and with good sensitivity in breath-hold experiments.

CONCLUSIONS

GE-SE EPIK provides increased resolution and more echoes, including two SEs, than comparable sequences. Moreover, GE-SE EPIK achieves this within an acquisition time of 57 s for 20 slices (matrix size = 128×128; FOV = 24 cm) and with a reasonably short TE for the final echo (114 ms). The sequence can dynamically track OEF changes in a breath-hold experiment.

The Compressed Sensing MP2RAGE as a Surrogate to the MPRAGE for Neuroimaging at 3 T

OBJECTIVES

The magnetization-prepared 2 rapid acquisition gradient echo (MP2RAGE) sequence provides quantitative T1 maps in addition to high-contrast morphological images. Advanced acceleration techniques such as compressed sensing (CS) allow its acquisition time to be compatible with clinical applications. To consider its routine use in future neuroimaging protocols, the repeatability of the segmented brain structures was evaluated and compared with the standard morphological sequence (magnetization-prepared rapid gradient echo [MPRAGE]). The repeatability of the T1 measurements was also assessed.

MATERIALS AND METHODS

Thirteen healthy volunteers were scanned either 3 or 4 times at several days of interval, on a 3 T clinical scanner, with the 2 sequences (CS-MP2RAGE and MPRAGE), set with the same spatial resolution (0.8-mm isotropic) and scan duration (6 minutes 21 seconds). The reconstruction time of the CS-MP2RAGE outputs (including the 2 echo images, the MP2RAGE image, and the T1 map) was 3 minutes 33 seconds, using an open-source in-house algorithm implemented in the Gadgetron framework.Both precision and variability of volume measurements obtained from CAT12 and VolBrain were assessed. The T1 accuracy and repeatability were measured on phantoms and on humans and were compared with literature.Volumes obtained from the CS-MP2RAGE and the MPRAGE images were compared using Student t tests (P < 0.05 was considered significant).

RESULTS

The CS-MP2RAGE acquisition provided morphological images of the same quality and higher contrasts than the standard MPRAGE images. Similar intravolunteer variabilities were obtained with the CS-MP2RAGE and the MPRAGE segmentations. In addition, high-resolution T1 maps were obtained from the CS-MP2RAGE. T1 times of white and gray matters and several deep gray nuclei are consistent with the literature and show very low variability (<1%).

CONCLUSIONS

The CS-MP2RAGE can be used in future protocols to rapidly obtain morphological images and quantitative T1 maps in 3-dimensions while maintaining high repeatability in volumetry and relaxation times.

In vivo T1 and T2 relaxation time maps of brain tissue, skeletal muscle, and lipid measured in healthy volunteers at 50 mT

PURPOSE

Low-field (B₀ < 0.1 T) MRI has generated much interest as a means of increased accessibility via reduced cost and improved portability compared to conventional clinical systems (B₀ ≥ 1.5 Tesla). Here we measure MR relaxation times at 50 mT and compare results with commonly used models based on both in vivo and ex vivo measurements.

METHODS

Using 3D turbo spin echo readouts, T₁ and T₂ maps of the human brain and lower leg were acquired on a custom-built 50 mT MRI scanner using inversion-recovery and multi-echo-based sequences, respectively. Image segmentation was performed based on a histogram analysis of the relaxation times.

RESULTS

The average T₁ times of gray matter, white matter, and cerebrospinal fluid (CSF) were 327 ± 10 ms, 275 ± 5 ms, and 3695 ± 287 ms, respectively. Corresponding values of T₂ were 102 ± 6 ms, 102 ± 6 ms, and 1584 ± 124 ms. T₁ times in the calf muscle were measured to be 171 ± 11 ms and were 130 ± 5 ms in subcutaneous and bone marrow lipid. Corresponding T₂ times were 39 ± 2 ms in muscle and 90 ± 13 ms in lipid.

CONCLUSIONS

For tissues except for CSF, the measured T₁ times are much shorter than reported at higher fields and generally lie within the range of different models in the literature. As expected, T₂ times are similar to those seen at typical clinical field strengths. Analysis of the relaxation maps indicates that segmentation of white and gray matter based purely on T₁ or T₂ will be quite challenging at low field given the relatively small difference in relaxation times.

Development and validation of 3D MP-SSFP to enable MRI in inhomogeneous magnetic fields

PURPOSE

We demonstrate the feasibility of MRI with missing-pulse steady-state free precession (MP-SSFP) in a 4T magnet with artificially degraded homogeneity.

METHODS

T₁ , T₂ , and diffusion contrast of MP-SSFP was simulated with constant and alternate radiofrequency (RF) phase using an extended phase graph. To validate MP-SSFP performance in human brain imaging, MP-SSFP was tested with two types of artificially introduced inhomogeneous magnetic fields: (1) a pure linear gradient field, and (2) a pseudo-linear gradient field introduced by mounting a head-gradient set at 36 cm from the magnet isocenter. Image distortion induced by the nonlinear inhomogeneous field was corrected using B₀ mapping measured with MP-SSFP.

RESULTS

The maximum flip angle in MP-SSFP was limited to ≤10° because of the large range of resonance frequencies in the inhomogeneous magnetic fields tested in this study. Under this flip-angle limitation, MP-SSFP with constant RF phase provided advantages of higher signal-to-noise ratio and insensitivity to B₁ ⁺ field inhomogeneity as compared with an alternate RF phase. In diffusion simulation, the steady-state magnetization in constant RF phase MP-SSFP increased with an increase of static field gradient up to 8 to 21 mT/m depending on simulation parameters. Experimental results at 4T validated these findings. In human brain imaging, MP-SSFP preserved sufficient signal intensities, but images showed severe image distortion from the pseudo-linear inhomogeneous field. However, following distortion correction, good-quality brain images were achieved.

CONCLUSION

MP-SSFP appears to be a feasible MRI technique for brain imaging in an inhomogeneous magnetic field.

Design of a fast field-cycling magnetic resonance imaging system, characterization and methods for relaxation dispersion measurements around 1.5 T

The dependence of the nuclear magnetic resonance relaxation rate on the magnetic field has been widely studied, in particular, in biomedical areas with the objectives to better understand the underlying microscopic mechanisms in tissues and provide biomarkers of diseases. By combining fast-field cycling (FFC) and magnetic resonance imaging (MRI), it is possible to provide localized relaxation dispersion measurements in heterogeneous systems with recent demonstrations in solutions, biological samples, human beings, and small animals. We report here the developments and performances of a device designed for small animal FFC-MRI comprising a resistive insert technology operating inside a 1.5 T MRI system. Specific measurement methods were developed to characterize the system efficiency, response time, homogeneity, stability, and compensation. By adding a non-linear element in the system and using a dual amplifier strategy, it is shown that large field offsets can be produced during relaxation periods while maintaining precise field control during detection periods. The measurement of longitudinal nuclear magnetic relaxation dispersion (NMRD) profiles in the range of 1.08 T-1.92 T is reported, essentially displaying a linear variation in this range for common MRI contrast agents. The slopes of both the longitudinal and transverse relaxation dispersion profiles at 1.5 T are measured and validated, extending the capabilities of previous approaches. The performances of a longitudinal relaxation dispersion mapping method are finally reported, opening the way to quantitative preclinical dispersion imaging studies at a high FFC-MRI field.

A Single-Scan, Rapid Whole-Brain Protocol for Quantitative Water Content Mapping With Neurobiological Implications

Water concentration is tightly regulated in the healthy human brain and changes only slightly with age and gender in healthy subjects. Consequently, changes in water content are important for the characterization of disease. MRI can be used to measure changes in brain water content, but as these changes are usually in the low percentage range, highly accurate and precise methods are required for detection. The method proposed here is based on a long-TR (10 s) multiple-echo gradient-echo measurement with an acquisition time of 7:21 min. Using such a long TR ensures that there is no T₁ weighting, meaning that the image intensity at zero echo time is only proportional to the water content, the transmit field, and to the receive field. The receive and transmit corrections, which are increasingly large at higher field strengths and for highly segmented coil arrays, are multiplicative and can be approached heuristically using a bias field correction. The method was tested on 21 healthy volunteers at 3T field strength. Calibration using cerebral-spinal fluid values (~100% water content) resulted in mean values and standard deviations of the water content distribution in white matter and gray matter of 69.1% (1.7%) and 83.7% (1.2%), respectively. Measured distributions were coil-independent, as seen by using either a 12-channel receiver coil or a 32-channel receiver coil. In a test-retest investigation using 12 scans on one volunteer, the variation in the mean value of water content for different tissue types was ~0.3% and the mean voxel variability was ~1%. Robustness against reduced SNR was assessed by comparing results for 5 additional volunteers at 1.5T and 3T. Furthermore, water content distribution in gray matter is investigated and regional contrast reported for the first time. Clinical applicability is illustrated with data from one stroke patient and one brain tumor patient. It is anticipated that this fast, stable, easy-to-use, high-quality mapping method will facilitate routine quantitative MR imaging of water content.

Multi-Exponential Relaxometry Using l1 -Regularized Iterative NNLS (MERLIN) With Application to Myelin Water Fraction Imaging

A new parameter estimation algorithm, MERLIN, is presented for accurate and robust multi-exponential relaxometry using magnetic resonance imaging, a tool that can provide valuable insight into the tissue microstructure of the brain. Multi-exponential relaxometry is used to analyze the myelin water fraction and can help to detect related diseases. However, the underlying problem is ill-conditioned, and as such, is extremely sensitive to noise and measurement imperfections, which can lead to less precise and more biased parameter estimates. MERLIN is a fully automated, multi-voxel approach that incorporates state-of-the-art l₁ -regularization to enforce sparsity and spatial consistency of the estimated distributions. The proposed method is validated in simulations and in vivo experiments, using a multi-echo gradient-echo (MEGE) sequence at 3 T. MERLIN is compared to the conventional single-voxel l₂ -regularized NNLS (rNNLS) and a multi-voxel extension with spatial priors (rNNLS + SP), where it consistently showed lower root mean squared errors of up to 70 percent for all parameters of interest in these simulations.

An MR technique for simultaneous quantitative imaging of water content, conductivity and susceptibility, with application to brain tumours using a 3T hybrid MR-PET scanner

Approaches for the quantitative mapping of water content, electrical conductivity and susceptibility have been developed independently. The purpose of this study is to develop a method for simultaneously acquiring quantitative water content, electrical conductivity and susceptibility maps based on a 2D multi-echo gradient echo sequence. Another purpose is to investigate the changes in these properties caused by brain tumours. This was done using a 3T hybrid magnetic resonance imaging and positron emission tomography (MR-PET) scanner. Water content maps were derived after performing T₂* and transmit-receive field bias corrections to magnitude images essentially reflecting only the H₂O content contrast. Phase evolution during the multi-echo train was used to generate field maps and derive quantitative susceptibility, while the conductivity maps were retrieved from the phase value at zero echo time. Performance of the method is demonstrated on phantoms and two healthy volunteers. In addition, the method was applied to three patients with brain tumours and a comparison to maps obtained from PET using O-(2-[18 F]fluoroethyl)-L-tyrosine and clinical MR images is presented. The combined information of the water content, conductivity and susceptibility may provide additional information about the tissue viability. Future studies can benefit from the evaluation of these contrasts with shortened acquisition times.

Physical and numerical phantoms for the validation of brain microstructural MRI: A cookbook

Phantoms, both numerical (software) and physical (hardware), can serve as a gold standard for the validation of MRI methods probing the brain microstructure. This review aims to provide guidelines on how to build, implement, or choose the right phantom for a particular application, along with an overview of the current state-of-the-art of phantoms dedicated to study brain microstructure with MRI. For physical phantoms, we discuss the essential requirements and relevant characteristics of both the (NMR visible) liquid and (NMR invisible) phantom materials that induce relevant microstructural features detectable via MRI, based on diffusion, intra-voxel incoherent motion, magnetization transfer or magnetic susceptibility weighted contrast. In particular, for diffusion MRI, many useful phantoms have been proposed, ranging from simple liquids to advanced biomimetic phantoms consisting of hollow or plain microfibers and capillaries. For numerical phantoms, the focus is on Monte Carlo simulations of random walk, for which the basic principles, along with useful criteria to check and potential pitfalls are reviewed, in addition to a literature overview highlighting recent advances. While many phantoms exist already, the current review aims to stimulate further research in the field and to address remaining needs.

A 3D two-point method for whole-brain water content and relaxation time mapping: Comparison with gold standard methods

Quantitative imaging of the human brain is of great interest in clinical research as it enables the identification of a range of MR biomarkers useful in diagnosis, treatment and prognosis of a wide spectrum of diseases. Here, a 3D two-point method for water content and relaxation time mapping is presented and compared to established gold standard methods. The method determines free water content, H2O, and the longitudinal relaxation time, T1, quantitatively from a two-point fit to the signal equation including corrections of the transmit and receive fields. In addition, the effective transverse relaxation time, T2*, is obtained from an exponential fit to the multi-echo signal train and its influence on H2O values is corrected. The phantom results obtained with the proposed method show good agreement for H2O and T1 values with known and spectroscopically measured values, respectively. The method is compared in vivo to already established gold standard quantitative methods. For H2O and T2* mapping, the 3D two-point results were compared to a measurement conducted with a multiple-echo GRE with long TR and T1 is compared to results from a Look-Locker method, TAPIR. In vivo results show good overall agreement between the methods, but some systematic deviations are present. Besides an expected dependence of T2* on voxel size, T1 values are systematically larger in the 3D approach than those obtained with the gold standard method. This behaviour might be due to imperfect spoiling, influencing each method differently. Results for H2O differ due to differences in the saturation of cerebrospinal fluid and partial volume effects. In addition, ground truth values of in vivo studies are unknown, even when comparing to in vivo gold standard methods. A detailed region-of-interest analysis for H2O and T1 matches well published literature values.

Deducing magnetic resonance neuroimages based on knowledge from samples

PURPOSE

Because individual variance always exists, using the same set of predetermined parameters for magnetic resonance imaging (MRI) may not be exactly suitable for each participant. We propose a knowledge-based method that can repair MRI data of undesired contrast as if a new scan were acquired using imaging parameters that had been individually optimized.

METHODS

The method employed a strategy called analogical reasoning to deduce voxel-wise relaxation properties using morphological and biological similarity. The proposed framework involves steps of intensity normalization, tissue segmentation, relaxation time deducing, and image deducing.

RESULTS

This approach has been preliminarily validated using conventional MRI data at 3T from several examples, including 5 normal and 9 clinical datasets. It can effectively improve the contrast of real MRI data by deducing imaging data using optimized imaging parameters based on deduced relaxation properties. The statistics of deduced images shows a high correlation with real data that were actually collected using the same set of imaging parameters.

CONCLUSION

The proposed method of deducing MRI data using knowledge of relaxation times alternatively provides a way of repairing MRI data of less optimal contrast. The method is also capable of optimizing an MRI protocol for individual participants, thereby realizing personalized MR imaging.

Quantitative in vivo T2 mapping using fast spin echo techniques - A linear correction procedure

A method is presented for correcting the effects of stimulated and indirect echoes on quantitative T2 mapping data acquired with multiple spin echo techniques, such as turbo spin echo. In contrast to similar correction techniques proposed in the literature, the method does not require a priori knowledge of the radio frequency (RF) pulse profiles. In a first step, for the T2 mapping protocol under investigation, signal decay curves S(TE) are simulated for a range of different RF pulse profiles. The actual signal decay S(TE) is then measured on a phantom with known T2, so the approximate RF pulse profiles can be derived via comparison with the simulated decay curves. In a second step, with the RF pulses obtained from step one, signal decay curves S(TE) are simulated for different T2 values and fitted mono-exponentially, thus allowing to deduce the relationship between true T2 and the apparent T2 (T2app) values. Results show that this relationship is approximately linear, allowing for a direct correction of T2app maps. If the amplitude of the transmitted RF field (B1) does not exceed the nominal value by more than 10%, it is shown that a B1-independent correction of T2app maps yields sufficiently accurate results for T2. A B1-dependent version is also presented. The method is tested in vitro on a phantom with different T2 values and in vivo on healthy subjects.

Segmenting and validating brain tissue definitions in the presence of varying tissue contrast

We propose a method for segmenting brain tissue as either gray matter or white matter in the presence of varying tissue contrast, which can derive from either differential changes in tissue water content or increasing myelin content of white matter. Our method models the spatial distribution of intensities as a Markov Random Field (MRF) and estimates the parameters for the MRF model using a maximum likelihood approach. Although previously described methods have used similar models to segment brain tissue, accurate model of the conditional probabilities of tissue intensities and adaptive estimates of tissue properties to local intensities generates tissue definitions that are accurate and robust to variations in tissue contrast with age and across illnesses. Robustness to variations in tissue contrast is important to understand normal brain development and to identify the brain bases of neurological and psychiatric illnesses. We used simulated brains of varying tissue contrast to compare both visually and quantitatively the performance of our method with the performance of prior methods. We assessed validity of the cortical definitions by associating cortical thickness with various demographic features, clinical measures, and medication use in our three large cohorts of participants who were either healthy or who had Bipolar Disorder (BD), Autism Spectrum Disorder (ASD), or familial risk for Major Depressive Disorder (MDD). We assessed validity of the tissue definitions using synthetic brains and data for three large cohort of individuals with various neuropsychiatric disorders. Visual inspection and quantitative analyses showed that our method accurately and robustly defined the cortical mantle in brain images with varying contrast. Furthermore, associating the thickness with various demographic and clinical measures generated findings that were novel and supported by histological analyses or were supported by previous MRI studies, thereby validating the cortical definitions generated by the proposed method: (1) Although cortical thickness decreased with age in adolescents, in adults cortical thickness did not correlate significantly with age. Our synthetic data showed that the previously reported thinning of cortex in adults is likely due to decease in tissue contrast, thereby suggesting that the method generated cortical definitions in adults that were invariant to tissue contrast. In adolescents, cortical thinning with age was preserved likely due to widespread dendritic and synaptic pruning, even though the effects of decreasing tissue contrast were minimized. (3) The method generated novel finding of both localized increases and decreases in thickness of males compared to females after controlling for the differing brain sizes, which are supported by the histological analyses of brain tissue in males and females. (4) The proposed method, unlike prior methods, defined thicker cortex in BD individuals using lithium. The novel finding is supported by the studies that showed lithium treatment increased dendritic arborization and neurogenesis, thereby leading to thickening of cortex. (5) In both BD and ASD participants, associations of more severe symptoms with thinner cortex showed that correcting for the effects of tissue contrast preserved the biological consequences of illnesses. Therefore, consistency of the findings across the three large cohorts of participants, in images acquired on either 1.5T or 3T MRI scanners, and with findings from prior histological analyses provides strong evidence that the proposed method generated valid and accurate definitions of the cortex while controlling for the effects of tissue contrast.

Multislice T1 -prepared 2D single-shot EPI: analysis of a clinical T1 mapping method unbiased by B0 or B1 inhomogeneity

Quantitative MR imaging is as sensitive in detecting lesions as qualitative imaging, but it is potentially more specific in differentiating disease. T1 mapping in particular might help to assess acute ischemic stroke, multiple sclerosis, epilepsy and Alzheimer's disease better. Thus, a rapid and robust clinical technique is vital. In 1990, Ordidge and colleagues developed the multislice T1 -prepared two-dimensional (2D) single-shot echo planar imaging technique. Subsequent studies demonstrated its clinical viability, but none performed an in-depth analysis of the strengths and advantages of this T1 mapping method. Herein, theoretical and experimental evidence shows that the technique accounts for 2D slice profile effects and is unbiased by B0 or B1 inhomogeneity. This is verified explicitly by varying the linear shims, the T1 preparation flip angle and the excitation flip angle. Furthermore, it is shown that the repetition time (and hence scan time) can be reduced without a loss of T1 accuracy. Copyright © 2016 John Wiley & Sons, Ltd.

Diffusion-weighted DESS protocol optimization for simultaneous mapping of the mean diffusivity, proton density and relaxation times at 3 Tesla

PURPOSE

To design a general framework for the optimization of an MRI protocol based on the the diffusion-weighted dual-echo steady-state (DW-DESS) sequence, enabling quantitative and simultaneous mapping of proton density (PD), relaxation times T1 and T2 and diffusion coefficient D.

METHODS

A parameterization of the DW-DESS sequence minimizing the Cramér-Rao lower bound of each parameter estimate was proposed and tested in a phantom experiment. An extension of the protocol was implemented for brain imaging to return the rotationally invariant mean diffusivity (MD).

RESULTS

In an NiCl₂ -doped agar gel phantom wherein T1/T2=920/65 ms, the parameter estimation errors were below 3% for PD and T1 and below 7% for T2 and D while the measured signal-to-noise ratio always exceeded 20. In the human brain, the in vivo parametric maps obtained were overall in reasonable agreement with gold standard measurements, despite a broadening of the distributions due to physiological motion.

CONCLUSION

Within the optimization framework presented here, DW-DESS images can be quantitatively interpreted to yield four intrinsic parameters of the tissue. Currently, the method is limited by the sensitivity of the DW-DESS sequence in terms of physiological motion. Magn Reson Med 78:130-141, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Evaluation of brain ageing: a quantitative longitudinal MRI study over 7 years

OBJECTIVES

T1 relaxometry is a promising tool for the assessment of microstructural changes during brain ageing. Previous cross-sectional studies demonstrated increasing T1 values in white and decreasing T1 values in grey matter over the lifetime. However, these findings have not yet been confirmed on the basis of a longitudinal study. In this longitudinal study over 7 years, T1 relaxometry was used to investigate the dynamics of age-related microstructural changes in older healthy subjects.

METHODS

T1 mapping was performed in 17 healthy subjects (range 51-77 years) at baseline and after 7 years. Advanced cortical and white matter segmentation was used to determine mean T1 values in the cortex and white matter.

RESULTS

The analysis revealed a decrease of mean cortical T1 values over 7 years, the rate of T1 reduction being more prominent in subjects with higher age. T1 decreases were predominantly localized in the lateral frontal, parietal and temporal cortex. In contrast, mean white matter T1 values remained stable.

CONCLUSIONS

T1 mapping is shown to be sensitive to age-related microstructural changes in healthy ageing subjects in a longitudinal setting. Data of a cohort in late adulthood and the senescence period demonstrate a decrease of cortical T1 values over 7 years, most likely reflecting decreasing water content and increased iron concentrations.

KEY POINTS

• T1 mapping is sensitive to age-related microstructural changes in a longitudinal setting. • T1 decreases were predominantly localized in the lateral frontal, parietal and temporal cortex. • The rate of T1 reduction was more prominent in subjects with higher age. • These changes most likely reflect decreasing cortical water and increasing iron concentrations.

Analytical corrections of banding artifacts in driven equilibrium single pulse observation of T2 (DESPOT2)

PURPOSE

DESPOT2 is a single-component T₂ mapping technique based on bSSFP imaging. It has seen limited application because of banding artifacts and magnetization transfer (MT) effects. In this work, acquisitions are optimized to minimize MT effects, while exact and approximate analytical equations enable automatic correction of banding artifacts within the T₂ maps in mere seconds.

THEORY AND METHODS

The technique was verified on an agar phantom at 3 tesla. The T₂ resulting from four different data combination techniques was compared with the T₂ from CPMG. Two comparable DESPOT2 scan protocols (short vs. long TR/T_RF ) designed to minimize MT effects, were tested both in the phantom and in vivo. A third protocol was tested in the brain of 8 volunteers and analytical correction schemes were compared with DESPOT2-FM.

RESULTS

The T₂ measurements in agar agree with CPMG within ∼7% and in vivo results agree with values reported in the literature. The approximate analytical solutions provide increased robustness to hardware imperfections and higher T₂ -to-noise ratio than the exact solutions.

CONCLUSION

New analytical solutions enable fast and accurate whole-brain T₂ mapping from previously measured T₁ and B₁ maps, and bSSFP images with at least two phase offsets and two flip angles (=4 datasets, 8 min scan). Magn Reson Med 76:1790-1804, 2016. © 2015 International Society for Magnetic Resonance in Medicine.

A review of optimization and quantification techniques for chemical exchange saturation transfer MRI toward sensitive in vivo imaging

Chemical exchange saturation transfer (CEST) MRI is a versatile imaging method that probes the chemical exchange between bulk water and exchangeable protons. CEST imaging indirectly detects dilute labile protons via bulk water signal changes following selective saturation of exchangeable protons, which offers substantial sensitivity enhancement and has sparked numerous biomedical applications. Over the past decade, CEST imaging techniques have rapidly evolved owing to contributions from multiple domains, including the development of CEST mathematical models, innovative contrast agent designs, sensitive data acquisition schemes, efficient field inhomogeneity correction algorithms, and quantitative CEST (qCEST) analysis. The CEST system that underlies the apparent CEST-weighted effect, however, is complex. The experimentally measurable CEST effect depends not only on parameters such as CEST agent concentration, pH and temperature, but also on relaxation rate, magnetic field strength and more importantly, experimental parameters including repetition time, RF irradiation amplitude and scheme, and image readout. Thorough understanding of the underlying CEST system using qCEST analysis may augment the diagnostic capability of conventional imaging. In this review, we provide a concise explanation of CEST acquisition methods and processing algorithms, including their advantages and limitations, for optimization and quantification of CEST MRI experiments.

IVIM analysis of brain tumors: an investigation of the relaxation effects of CSF, blood, and tumor tissue on the estimated perfusion fraction

OBJECT

We sought to investigate the dependence of intravoxel incoherent motion (IVIM)-related perfusion fraction (f) estimates on the transverse relaxation of brain tissue, blood, and cerebrospinal fluid (CSF), attempting to overcome the influence of CSF on conventional f maps.

MATERIALS AND METHODS

Eighteen patients with gliomas underwent DWI with 14 b-values (0-1,300 s/mm(2)) and two distinct echo times (TEs). Regions of interest representing tumour and normal brain tissue were analysed by calculating the f values for both TEs. A mask for pixels with relevant CSF partial volume was subsequently created. The f values were tested for significant differences.

RESULTS

We found statistically significant differences between the two TEs in the f values for cortical and juxtacortical structures and non-enhancing areas of the tumour /oedema. Normal white matter and gadolinium-enhancing tumour tissue appeared insensitive to TE variation. In all tissue types examined, the masking of voxels with considerable CSF content was able to overcome issues of erroneous f estimation and calculation of f values insensitive to TE changes was feasible.

CONCLUSION

Due to the complex interaction in the relaxation rates of CSF, blood, and tumour tissue, the estimation of f values is affected by the choice of TE. Only f values in normal white matter and tumour tissue--which largely comprise blood voxels with minor CSF partial volume--may be clinically applicable in the present form of IVIM-based DWI analysis. Going a step further, and after removing voxels with heavily TE-susceptible f values, we were able to obtain accurate and TE-independent f values in contrast-enhancing tumour tissue, white matter, and grey matter, which were essentially consistent with those reported in the literature.

Technical considerations on the validity of blood oxygenation level-dependent-based MR assessment of vascular deoxygenation

A blood oxygenation level-dependent (BOLD)-based apparent relative oxygen extraction fraction (rOEF) as a semi-quantitative marker of vascular deoxygenation has recently been introduced in clinical studies of patients with glioma and stroke, yielding promising results. These rOEF measurements are based on independent quantification of the transverse relaxation times T2 and T2* and relative cerebral blood volume (rCBV). Simulations demonstrate that small errors in any of the underlying measures may result in a large deviation of the calculated rOEF. Therefore, we investigated the validity of such measurements. For this, we evaluated the quantitative measurements of T2 and T2* at 3 T in a gel phantom, in healthy subjects and in healthy tissue of patients with brain tumors. We calculated rOEF maps covering large portions of the brain from T2, T2* and rCBV [routinely measured in patients using dynamic susceptibility contrast (DSC)], and obtained rOEF values of 0.63 ± 0.16 and 0.90 ± 0.21 in healthy-appearing gray matter (GM) and white matter (WM), respectively; values of about 0.4 are usually reported. Quantitative T2 mapping using the fast, clinically feasible, multi-echo gradient spin echo (GRASE) approach yields significantly higher values than much slower multiple single spin echo (SE) experiments. Although T2* mapping is reliable in magnetically homogeneous tissues, uncorrectable macroscopic background gradients and other effects (e.g. iron deposition) shorten T2*. Cerebral blood volume (CBV) measurement using DSC and normalization to WM yields robust estimates of rCBV in healthy-appearing brain tissue; absolute quantification of the venous fraction of CBV, however, is difficult to achieve. Our study demonstrates that quantitative measurements of rOEF are currently biased by inherent difficulties in T2 and CBV quantification, but also by inadequacies of the underlying model. We argue, however, that standardized, reproducible measurements of apparent T2, T2* and rCBV may still allow the estimation of a meaningful apparent rOEF, which requires further validation in clinical studies.

Effect of magnetic field and iron content on NMR proton relaxation of liver, spleen and brain tissues

Iron accumulation is observed in liver and spleen during hemochromatosis and important neurodegenerative diseases involve iron overload in brain. Storage of iron is ensured by ferritin, which contains a magnetic core. It causes a darkening on T2 -weighted MR images. This work aims at improving the understanding of the NMR relaxation of iron-loaded human tissues, which is necessary to develop protocols of iron content measurements by MRI. Relaxation times measurements on brain, liver and spleen samples were realized at different magnetic fields. Iron content was determined by atomic emission spectroscopy. For all samples, the longitudinal relaxation rate (1/T1 ) of tissue protons decreases with the magnetic field up to 1 T, independently of iron content, while their transverse relaxation rate (1/T2 ) strongly increases with the field, either linearly or quadratically, or a combination thereof. The extent of the inter-echo time dependence of 1/T2 also varies according to the sample. A combination of theoretical models is necessary to describe the relaxation of iron-containing tissues. This can be due to the presence, inside tissues, of ferritin clusters of different sizes and densities. When considering all samples, a correlation (r(2)  = 0.6) between 1/T1 and iron concentration is observed at 7.0 T. In contrast the correlation between 1/T2 and iron content is poor, even at high field (r(2)  = 0.14 at 7.0 T). Our results show that MRI methods based on T1 or T2 measurements will easily detect an iron overloading at high magnetic field, but will not provide an accurate quantification of tissue iron content at low iron concentrations.

Analysis of proton-density bias corrections based on T1 measurement for robust quantification of water content in the brain at 3 Tesla

PURPOSE

Estimating tissue water content using high field MRI, such as 3 Tesla (T), is challenging due to the difficulty in dissociating the radio frequency inhomogeneity pattern from the signal arising from tissue intrinsic proton density (PD) variations. To overcome this problem the longitudinal relaxation time T1 can be combined with an initial guess of the PD to yield the desired PD bias correction. However, it is necessary to know whether T1 effects, i.e., any effect contributing to T1 while being independent of tissue hydration, influence the estimated correction.

METHODS

Twenty-five healthy subjects underwent a quantitative 3T MRI protocol enabling acquisition of 64 slices with 1 mm in-plane resolution and 2 mm slice thickness in 14 min. Influence of T1 effects on the estimated water content map is evaluated using a dedicated method including T1 and T2 * information and region of interest-based water content values are compared with the literature.

RESULTS

Our analysis indicates that the PD bias correction based on T1 is largely insensitive to T1 effects. Besides, water content results are in good agreement with literature values obtained at 1.5T.

CONCLUSION

This study demonstrates the applicability of a PD bias correction based on T1 to yield tissue water content at 3T.

The effects of changing water content, relaxation times, and tissue contrast on tissue segmentation and measures of cortical anatomy in MR images

Water content is the dominant chemical compound in the brain and it is the primary determinant of tissue contrast in magnetic resonance (MR) images. Water content varies greatly between individuals, and it changes dramatically over time from birth through senescence of the human life span. We hypothesize that the effects that individual- and age-related variations in water content have on contrast of the brain in MR images also have important, systematic effects on in vivo, MRI-based measures of regional brain volumes. We also hypothesize that changes in water content and tissue contrast across time may account for age-related changes in regional volumes, and that differences in water content or tissue contrast across differing neuropsychiatric diagnoses may account for differences in regional volumes across diagnostic groups. We demonstrate in several complementary ways that subtle variations in water content across age and tissue compartments alter tissue contrast, and that changing tissue contrast in turn alters measures of the thickness and volume of the cortical mantle: (1) We derive analytic relations describing how age-related changes in tissue relaxation times produce age-related changes in tissue gray-scale intensity values and tissue contrast; (2) We vary tissue contrast in computer-generated images to assess its effects on tissue segmentation and volumes of gray matter and white matter; and (3) We use real-world imaging data from adults with either Schizophrenia or Bipolar Disorder and age- and sex-matched healthy adults to assess the ways in which variations in tissue contrast across diagnoses affects group differences in tissue segmentation and associated volumes. We conclude that in vivo MRI-based morphological measures of the brain, including regional volumes and measures of cortical thickness, are a product of, or at least are confounded by, differences in tissue contrast across individuals, ages, and diagnostic groups, and that differences in tissue contrast in turn likely derive from corresponding differences in water content of the brain across individuals, ages, and diagnostic groups.

MR parameter quantification with magnetization-prepared double echo steady-state (MP-DESS)

PURPOSE

The mapping of MR relaxation times and proton density has been the subject of research in medical imaging for many years, as it offers the possibility for longitudinal investigation of disease and the correlation with related biochemical processes. The purpose of this study is to provide a fast mapping protocol, which simultaneously acquires MR relaxation times and relative proton density without compromising accuracy and precision.

METHODS

This work presents a novel magnetization-prepared double echo steady-state (MP-DESS) sequence, which was designed to be sensitive to parameter variations of interest, and insensitive to variations of confounding variables. It provides high sensitivity against variations of the MR relaxation times, high acquisition efficiency, and it is insensitive to off-resonance. Accurate phase graph modeling of the MP-DESS signal is used to obtain unbiased parameter estimates.

RESULTS

The approach is validated in phantom and in vivo measurements. A whole-brain acquisition of 1.4-mm isotropic resolution was acquired in 15 min. Comparisons to gold-standard methods suggest a mapping precision of 5% for T1 and M0 , and below 10% for T2.

CONCLUSION

A new quantitative imaging technique is introduced that allows fast and isotropic simultaneous MR parameter mapping of T1, T2, and M0.

Toward real-time temperature monitoring in fat and aqueous tissue during magnetic resonance-guided high-intensity focused ultrasound using a three-dimensional proton resonance frequency T1 method

PURPOSE

To present a three-dimensional (3D) segmented echoplanar imaging (EPI) pulse sequence implementation that provides simultaneously the proton resonance frequency shift temperature of aqueous tissue and the longitudinal relaxation time (T1 ) of fat during thermal ablation.

METHODS

The hybrid sequence was implemented by combining a 3D segmented flyback EPI sequence, the extended two-point Dixon fat and water separation, and the double flip angle T1 mapping techniques. High-intensity focused ultrasound (HIFU) heating experiments were performed at three different acoustic powers on excised human breast fat embedded in ex vivo porcine muscle. Furthermore, T1 calibrations with temperature in four different excised breast fat samples were performed, yielding an estimate of the average and variation of dT1 /dT across subjects.

RESULTS

The water only images were used to mask the complex original data before computing the proton resonance frequency shift. T1 values were calculated from the fat-only images. The relative temperature coefficients were found in five fat tissue samples from different patients and ranged from 1.2% to 2.6%/°C.

CONCLUSION

The results demonstrate the capability of real-time simultaneous temperature mapping in aqueous tissue and T1 mapping in fat during HIFU ablation, providing a potential tool for treatment monitoring in organs with large fat content, such as the breast.

Imaging white matter in human brainstem

The human brainstem is critical for the control of many life-sustaining functions, such as consciousness, respiration, sleep, and transfer of sensory and motor information between the brain and the spinal cord. Most of our knowledge about structure and organization of white and gray matter within the brainstem is derived from ex vivo dissection and histology studies. However, these methods cannot be applied to study structural architecture in live human participants. Tractography from diffusion-weighted magnetic resonance imaging (MRI) may provide valuable insights about white matter organization within the brainstem in vivo. However, this method presents technical challenges in vivo due to susceptibility artifacts, functionally dense anatomy, as well as pulsatile and respiratory motion. To investigate the limits of MR tractography, we present results from high angular resolution diffusion imaging of an intact excised human brainstem performed at 11.1 T using isotropic resolution of 0.333, 1, and 2 mm, with the latter reflecting resolution currently used clinically. At the highest resolution, the dense fiber architecture of the brainstem is evident, but the definition of structures degrades as resolution decreases. In particular, the inferred corticopontine/corticospinal tracts (CPT/CST), superior (SCP) and middle cerebellar peduncle (MCP), and medial lemniscus (ML) pathways are clearly discernable and follow known anatomical trajectories at the highest spatial resolution. At lower resolutions, the CST/CPT, SCP, and MCP pathways are artificially enlarged due to inclusion of collinear and crossing fibers not inherent to these three pathways. The inferred ML pathways appear smaller at lower resolutions, indicating insufficient spatial information to successfully resolve smaller fiber pathways. Our results suggest that white matter tractography maps derived from the excised brainstem can be used to guide the study of the brainstem architecture using diffusion MRI in vivo.

Quantitative multi-parameter mapping of R1, PD(*), MT, and R2(*) at 3T: a multi-center validation

Multi-center studies using magnetic resonance imaging facilitate studying small effect sizes, global population variance and rare diseases. The reliability and sensitivity of these multi-center studies crucially depend on the comparability of the data generated at different sites and time points. The level of inter-site comparability is still controversial for conventional anatomical T1-weighted MRI data. Quantitative multi-parameter mapping (MPM) was designed to provide MR parameter measures that are comparable across sites and time points, i.e., 1 mm high-resolution maps of the longitudinal relaxation rate (R1 = 1/T1), effective proton density (PD(*)), magnetization transfer saturation (MT) and effective transverse relaxation rate (R2(*) = 1/T2(*)). MPM was validated at 3T for use in multi-center studies by scanning five volunteers at three different sites. We determined the inter-site bias, inter-site and intra-site coefficient of variation (CoV) for typical morphometric measures [i.e., gray matter (GM) probability maps used in voxel-based morphometry] and the four quantitative parameters. The inter-site bias and CoV were smaller than 3.1 and 8%, respectively, except for the inter-site CoV of R2(*) (<20%). The GM probability maps based on the MT parameter maps had a 14% higher inter-site reproducibility than maps based on conventional T1-weighted images. The low inter-site bias and variance in the parameters and derived GM probability maps confirm the high comparability of the quantitative maps across sites and time points. The reliability, short acquisition time, high resolution and the detailed insights into the brain microstructure provided by MPM makes it an efficient tool for multi-center imaging studies.

Trends and properties of human cerebral cortex: correlations with cortical myelin content

"In vivo Brodmann mapping" or non-invasive cortical parcellation using MRI, especially by measuring cortical myelination, has recently become a popular research topic, though myeloarchitectonic cortical parcellation in humans previously languished in favor of cytoarchitecture. We review recent in vivo myelin mapping studies and discuss some of the different methods for estimating myelin content. We discuss some ways in which myelin maps may improve surface registration and be useful for cross-modal and cross-species comparisons, including some preliminary cross-species results. Next, we consider neurobiological aspects of why some parts of cortex are more myelinated than others. Myelin content is inversely correlated with intracortical circuit complexity - in general, more myelin content means simpler and perhaps less dynamic intracortical circuits. Using existing PET data and functional network parcellations, we examine metabolic differences in the differently myelinated cortical functional networks. Lightly myelinated cognitive association networks tend to have higher aerobic glycolysis than heavily myelinated early sensory-motor ones, perhaps reflecting greater ongoing dynamic anabolic cortical processes. This finding is consistent with the hypothesis that intracortical myelination may stabilize intracortical circuits and inhibit synaptic plasticity. Finally, we discuss the future of the in vivo myeloarchitectural field and cortical parcellation--"in vivo Brodmann mapping"--in general.

In vivo functional and myeloarchitectonic mapping of human primary auditory areas

In contrast to vision, where retinotopic mapping alone can define areal borders, primary auditory areas such as A1 are best delineated by combining in vivo tonotopic mapping with postmortem cyto- or myeloarchitectonics from the same individual. We combined high-resolution (800 μm) quantitative T(1) mapping with phase-encoded tonotopic methods to map primary auditory areas (A1 and R) within the "auditory core" of human volunteers. We first quantitatively characterize the highly myelinated auditory core in terms of shape, area, cortical depth profile, and position, with our data showing considerable correspondence to postmortem myeloarchitectonic studies, both in cross-participant averages and in individuals. The core region contains two "mirror-image" tonotopic maps oriented along the same axis as observed in macaque and owl monkey. We suggest that these two maps within the core are the human analogs of primate auditory areas A1 and R. The core occupies a much smaller portion of tonotopically organized cortex on the superior temporal plane and gyrus than is generally supposed. The multimodal approach to defining the auditory core will facilitate investigations of structure-function relationships, comparative neuroanatomical studies, and promises new biomarkers for diagnosis and clinical studies.

In vivo measurement of transverse relaxation time in the mouse brain at 17.6 T

PURPOSE

To establish regional T1 and T2 values of the healthy mouse brain at ultra-high magnetic field strength of 17.6 T and to follow regional brain T1 and T2 changes with age.

METHODS

In vivo T1 and T2 values in the C57BL/6J mouse brain were followed with age using multislice-multiecho sequence and multiple spin echo saturation recovery with variable repetition time sequence, respectively, at 9.4 and 17.6 T. Gadolinium-tetra-azacyclo-dodecane-tetra-acetic acid phantoms were used to validate in vivo T2 measurements. Student's t-test was used to compare mean relaxation values.

RESULTS

A field-dependent decrease in T2 is shown and validated with phantom measurements. T2 values at 17.6 T typically increased with age in multiple brain regions except in the hypothalamus and the caudate-putamen, where a slight decrease was observed. Furthermore, T1 values in various brain regions of young and old mice are presented at 17.6 T. A large gain in signal-to-noise ratio was observed at 17.6 T.

CONCLUSIONS

This study establishes for the first time the normative T1 and T2 values at 17.6 T over different mouse brain regions with age. The estimates of in vivo T1 and T2 will be useful to optimize pulse sequences for optimal image contrast at 17.6 T and will serve as baseline values against which disease-related relaxation changes can be assessed in mice.

Automated detection of white matter signal abnormality using T2 relaxometry: application to brain segmentation on term MRI in very preterm infants

Hyperintense white matter signal abnormalities, also called diffuse excessive high signal intensity (DEHSI), are observed in up to 80% of very preterm infants on T2-weighted MRI scans at term-equivalent age. DEHSI may represent a developmental stage or diffuse microstructural white matter abnormalities. Automated quantitative assessment of DEHSI severity may help resolve this debate and improve neonatal brain tissue segmentation. For T2-weighted sequence without fluid attenuation, the signal intensity distribution of DEHSI greatly overlaps with that of cerebrospinal fluid (CSF) making its detection difficult. Furthermore, signal intensities of T2-weighted images are susceptible to magnetic field inhomogeneity. Increased signal intensities caused by field inhomogeneity may be confused with DEHSI. To overcome these challenges, we propose an algorithm to detect DEHSI using T2 relaxometry, whose reflection of the rapid changes in free water content provides improved distinction between CSF and DEHSI over that of conventional T2-weighted imaging. Moreover, the parametric transverse relaxation time T2 is invulnerable to magnetic field inhomogeneity. We conducted computer simulations to select an optimal detection parameter and to validate the proposed method. We also demonstrated that brain tissue segmentation is further enhanced by incorporating DEHSI detection for both simulated preterm infant brain images and in vivo in very preterm infants imaged at term-equivalent age.

Investigation of the spatial correlation in human white matter and the influence of age using 3-dimensional variography applied to MP-RAGE data

A novel method for the quantification of heterogeneity and spatial correlation in 3D MP-RAGE images of white matter is presented. The technique is based on the variogram, a tool commonly used in geosciences for the analysis of spatial data, and was tailored to the special requirements of MR image analysis. Influences from intensity non-uniformities, noise and arbitrary greyscale were quantified and considered in the calculations. The obtained variograms were fitted with spherical model functions to infer parameters that quantify heterogeneity and size of the correlation structures of the tissue. Numerically generated samples with well-defined correlation properties were employed to validate the estimation process and to provide an interpretation of the parameters obtained. It is shown that the method gives reliable results in an interval of correlation structures sized between 2mm and 20mm. The method was applied to 24 MP-RAGE datasets of healthy female volunteers ranging in age from 19 to 73 years. White matter was found to have two prominent correlation structures with sizes of approximately 3mm and 23 mm. The heterogeneity of the smaller structure increases significantly with age (r=0.83, p<10(-6)).

Brain iron deposition and sequence characteristics in Parkinsonism: comparison of SWI, T₂* maps, T₂-weighted-, and FLAIR-SPACE

OBJECTIVES

To compare quantitatively T2- and T2*-based magnetic resonance imaging sequences in patients with symptoms of Parkinson disease and to evaluate the information content of those sequences regarding brain iron concentration.

MATERIALS AND METHODS

We imaged 51 patients with symptoms of Parkinson disease on 3-T magnetic resonance imaging with T2-weighted sampling perfection with application optimized contrasts using different flip-angle evolution (SPACE), fluid attenuation inversion recovery (FLAIR)-SPACE, susceptibility-weighted imaging (SWI), and parametric T2* sequence (MapIt). Signal analysis was performed in 22 regions of interest in the brain.

RESULTS

Correlations (r2 = 0.82...0.96) with brain iron concentration were excellent. Contrast and tissue separability ratios were best in the T2* maps and FLAIR-SPACE, respectively. Good correlations of contrast were reached between SWI and both T2-weighted SPACE and FLAIR-SPACE. Their relation to quantitative T2* values was reminiscent of a quadratic curve shape. However, separation into gray and white matter revealed a linear positive and negative correlation, respectively.

CONCLUSIONS

SWI showed potential in differentiating illnesses characterized by brain iron deposition. Closely similar information was given by T2-weighted SPACE and FLAIR-SPACE, whereas other sequence comparisons revealed dispersion from intersequence agreement.

In vivo quantification of T₂ anisotropy in white matter fibers in marmoset monkeys

T₂-weighted MRI at high field is a promising approach for studying noninvasively the tissue structure and composition of the brain. However, the biophysical origin of T₂ contrast, especially in white matter, remains poorly understood. Recent work has shown that R₂ (=1/T₂) may depend on the tissue's orientation relative to the static magnetic field (B(0)) and suggested that this dependence could be attributed to local anisotropy in the magnetic properties of brain tissue. In the present work, we analyzed high-resolution, multi-gradient-echo images of in vivo marmoset brains at 7T, and compared them with ex vivo diffusion tensor images, to show that R₂ relaxation in white matter is highly sensitive to the fiber orientation relative to the main field. We directly demonstrate this orientation dependence by performing in vivo multi-gradient-echo experiments in two orthogonal brain positions, uncovering a nearly 50% change in the R₂ relaxation rate constant of the optic radiations. We attribute this substantial R₂ anisotropy to local subvoxel susceptibility effects arising from the highly ordered and anisotropic structure of the myelin sheath.

Exploiting spatial information to estimate metabolite levels in two-dimensional MRSI of heterogeneous brain lesions

MRSI provides MR spectra from multiple adjacent voxels within a body volume represented as a two- or three-dimensional matrix, allowing the measurement of the distribution of metabolites over this volume. The spectra of these voxels are usually analyzed one by one, without exploiting their spatial context. In this article, we present an advanced metabolite quantification method for MRSI data, in which the available spatial information is considered. A nonlinear least-squares algorithm is proposed in which prior knowledge is included in the form of proximity constraints on the spectral parameters within a grid and optimized starting values. A penalty term that promotes a spatially smooth spectral parameter map is added to the fitting algorithm. This method is adaptive, in the sense that several sweeps through the grid are performed and each solution may tune some hyperparameters at run-time. Simulation studies of MRSI data showed significantly improved metabolite estimates after the inclusion of spatial information. Improved metabolite maps were also demonstrated by applying the method to in vivo MRSI data. Overlapping peaks or peaks of compounds present at low concentration can be better quantified with the proposed method than with single-voxel approaches. The new approach compares favorably against the multivoxel approach embedded in the well-known quantification software LCModel.

Measuring the absolute water content of the brain using quantitative MRI

Methods for quantitative imaging of the brain are presented and compared. Highly precise and accurate mapping of the absolute water content and distribution, as presented here, requires a significant number of corrections and also involves mapping of other MR parameters. Here, either T(1) and T(2)(*) or T(2) is mapped, and several corrections involving the measurement of temperature, transmit and receive B(1) inhomogeneities and signal extrapolation to zero TE are applied. Information about the water content of the whole brain can be acquired in clinically acceptable measurement times (10 or 20 min). Since water content is highly regulated in the healthy brain, pathological changes can be easily identified and their evolution or correlation with other manifestations of the disease investigated. In addition to voxel-based total water content, information about the different environments of water can be gleaned from qMRI. The myelin water fraction can be extracted from the fit of very high-SNR multiple-echo T(2) decay curves with a superposition of a large number of exponentials. Diseases involving de- or dysmyelination can be investigated and lead to novel observations regarding the water compartmentalisation in tissue, despite the limited spatial coverage. In conclusion, quantitative MRI is emerging as an unparalleled tool for the study of the normal and diseased brain, replacing the customary time-space environment of the sequential mixed-contrast MRI with a multi-NMR-parametric space in which tissue microscopy is increasingly revealed.

[Relaxation times T1, T2, and T2* of apples, pears, citrus fruits, and potatoes with a comparison to human tissues]

The aim of the project was a systematic assessment of relaxation times of different fruits and vegetables and a comparison to values of human tissues. Results provide an improved basis for selection of plant phantoms for development of new MR techniques and sequences. Vessels filled with agar gel are mostly used for this purpose, preparation of which is effortful and time-consuming. In the presented study apples, (malus, 8 species), pears, (pyrus, 2 species), citrus fruits (citrus, 5 species) and uncooked potatoes (solanum tuberosum, 8 species) from the supermarket were examined which are easily available nearly all-the-year. T1, T2 and T2* relaxation times of these nature products were measured on a 1.5 Tesla MR system with adapted examination protocols and mono-exponential fitting, and compared to literature data of human parenchyma tissues, fatty tissue and body fluid (cerebrospinal fluid). Resulting values were as follows: apples: T1: 1486-1874 ms, T2: 163-281 ms, T2*: 2.3-3.2 ms; pears: T1: 1631-1969 ms, T2: 119-133 ms, T2* : 10.1-10.6 ms, citrus fruits (pulp) T1: 2055-2632 ms, T2: 497-998 ms, T2* : 151-182 ms; citrus fruits (skin) T1: 561-1669 ms, T2: 93-119 ms; potatoes: T1: 1011-1459 ms, T2: 166 - 210 ms, T2* : 20 - 30 ms. All T1-values of the examined objects (except for potatoes and skins of citrus fruits) were longer than T1 values of human tissues. Also T2 values (except for pears and skins of citrus fruits) of the fruits and the potatoes tended to be longer. T2* values of apples, pears and potatoes were shorter than in healthy human tissue. Results show relaxation values of many fruits to be not exactly fitting to human tissue, but with suitable selection of the fruits and optionally with an adaption of measurement parameters one can achieve suitable contrast and signal characteristics for some purposes.

Unified segmentation based correction of R1 brain maps for RF transmit field inhomogeneities (UNICORT)

Quantitative mapping of the longitudinal relaxation rate (R1 = 1/T1) in the human brain enables the investigation of tissue microstructure and macroscopic morphology which are becoming increasingly important for clinical and neuroimaging applications. R1 maps are now commonly estimated from two fast high-resolution 3D FLASH acquisitions with variable excitation flip angles, because this approach is fast and does not rely on special acquisition techniques. However, these R1 maps need to be corrected for bias due to RF transmit field (B1⁺) inhomogeneities, requiring additional B1⁺ mapping which is usually time consuming and difficult to implement. We propose a technique that simultaneously estimates the B1⁺ inhomogeneities and R1 values from the uncorrected R1 maps in the human brain without need for B1⁺ mapping. It employs a probabilistic framework for unified segmentation based correction of R1 maps for B1⁺ inhomogeneities (UNICORT). The framework incorporates a physically informed generative model of smooth B1⁺ inhomogeneities and their multiplicative effect on R1 estimates. Extensive cross-validation with the established standard using measured B1⁺ maps shows that UNICORT yields accurate B1⁺ and R1 maps with a mean deviation from the standard of less than 4.3% and 5%, respectively. The results of different groups of subjects with a wide age range and different levels of atypical brain anatomy further suggest that the method is robust and generalizes well to wider populations. UNICORT is easy to apply, as it is computationally efficient and its basic framework is implemented as part of the tissue segmentation in SPM8.

Automated brain tissue segmentation based on fractional signal mapping from inversion recovery Look-Locker acquisition

Most current automated segmentation methods are performed on T(1)- or T(2)-weighted MR images, relying on relative image intensity that is dependent on other MR parameters and sensitive to B(1) magnetic field inhomogeneity. Here, we propose an image segmentation method based on quantitative longitudinal magnetization relaxation time (T(1)) of brain tissues. Considering the partial volume effect, fractional volume maps of brain tissues (white matter, gray matter, and cerebrospinal fluid) were obtained by fitting the observed signal in an inversion recovery procedure to a linear combination of three exponential functions, which represents the relaxations of each of the tissue types. A Look-Locker acquisition was employed to accelerate the acquisition process. The feasibility and efficacy of this proposed method were evaluated using simulations and experiments. The potential applications of this method in the study of neurological disease as well as normal brain development and aging are discussed.

Accuracy of real-time MR temperature mapping in the brain: a comparison of fast sequences

PURPOSE

To compare magnetic resonance (MR) thermometry based on the proton resonance frequency (PRF) method using a single shot echoplanar imaging (ss EPI) sequence to both of the standard sequences, gradient echo (GRE) and segmented echoplanar imaging (seg EPI) in the in vivo human brain, at 1.5T and 3T.

MATERIAL AND METHODS

Repetitive MR thermometry was performed on the brain of six volunteers using GRE, seg EPI, and ss EPI sequences on whole-body 1.5T and 3T clinical systems using comparable acquisition parameters. Phase stability and temperature data precision in the human head were determined over 12 min for the three sequences at both field strengths. An ex-vivo swine skeletal muscle model was used to evaluate temperature accuracy of the ss EPI sequence during heating by high intensity focused ultrasound (HIFU).

RESULTS

In-vivo examinations of brain revealed an average temperature precision of 0.37 °C/0.39 °C/0.16 °C at 3T for the GRE/seg EPI/ss EPI sequences. At 1.5T, a precision of 0.58 °C/0.63 °C/0.21 °C was achieved. In the ex-vivo swine model, a strong correlation of temperature data derived using ss EPI and GRE sequences was found with a temperature deviation <1 °C.

CONCLUSION

The ss EPI sequence was the fastest and the most precise sequence for MR thermometry, with significantly higher accuracy compared to GRE.

Quantification of ethanol methyl (1)H magnetic resonance signal intensity following intravenous ethanol administration in primate brain

In vivo(1)H magnetic resonance spectroscopy (MRS) can be used to directly monitor brain ethanol. Previously, studies of human subjects have lead to the suggestion that the ethanol methyl (1)H MRS signal intensity relates to tolerance to ethanol's intoxicating effects. More recently, the ethanol (1)H MRS signal intensity has been recognized to vary between brain gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) due to differences in T(2) within these environments. The methods presented here extend ethanol MRS techniques to non-human primate subjects. Twelve monkeys were administered ethanol while sedated and positioned within a 3T MRI system. Chemical shift imaging (CSI) measurements were performed following intravenous infusion of 1g/kg ethanol. Magnetic resonance imaging (MRI) data were also recorded for each monkey to provide volume fractions of GM, WM, and CSF for each CSI spectrum. To estimate co-variance of ethanol MRS intensity with GM, WM, and CSF volume fractions, the relative contribution of each tissue subtype was determined following corrections for radiofrequency pulse profile non-uniformity, chemical shift artifacts, and differences between the point spread function in the CSI data and the imaging data. The ethanol MRS intensity per unit blood ethanol concentration was found to differ between GM, WM, and CSF. Individual differences in MRS intensity were larger in GM than WM. This methodology demonstrates the feasibility of ethanol MRS experiments and analysis in non-human primate subjects, and suggests GM may be a site of significant variation in ethanol MRS intensity between individuals.

Effects of echo time on diffusion quantification of brain white matter at 1.5 T and 3.0 T

The aim was to investigate the effects of echo time (TE) on diffusion quantification of brain white matter. Seven rhesus monkeys (all males; age, 4-6 years; weight, 5-7 kg) underwent diffusion tensor imaging (DTI) with a series of TEs in 1.5 T and 3.0 T MR scanners. The mean diffusivity (MD), fractional anisotropy (FA), primary (lambda(1)), and transverse eigenvalues (lambda(23)) were measured in a region of interest at the bilateral internal capsule. Pearson correlation showed that the FA and lambda(1) increased and lambda(23) decreased with TE both at 1.5 T and 3.0 T except for the MD. Repeated measurement analysis of variance (ANOVA) also showed significantly higher FA and lower MD and lambda(23) at 3.0 T than those at 1.5 T (P<0.01), but no statistical differences were found in lambda(1) between these two field strengths (P=0.709). These findings implied that TE and field strength might influence diffusion quantification in brain white matter.