Studying cyto and myeloarchitecture of the human cortex at ultra-high field with quantitative imaging: R1, R2* and magnetic susceptibility.

In this manuscript, the use of quantitative imaging at ultra-high field is evaluated as a mean to study cyto and myelo-architecture of the cortex. The quantitative contrasts used are the longitudinal relaxation rate (R1), apparent transverse relaxation rate (R2*) and quantitative susceptibility mapping (QSM). The quantitative contrasts were computed using high resolution in-vivo (0.65mm isotropic) brain data acquired at 7T. The performance of the different quantitative approaches was evaluated by visualizing the contrast between known highly myelinated primary sensory cortex regions and the neighbouring cortex. The transition from the inner layers to the outer layers (from white matter to the pial surface) of the human cortex, which is known to have varying cyto- and myelo architecture, was evaluated. The across cortex and through depth behaviour observed for the different quantitative maps was in good agreement between the different subjects, clearly allowing the differentiation between different Brodmann regions, suggesting these features could be used for individual cortical brain parcellation. While both R1 and R2* maps decrease monotonically from the white matter to the pial surface due to the decrease of myelin and iron between these regions, magnetic susceptibility maps have a more complex behaviour reflecting its opposing sensitivity to myelin and iron concentration.

An illustrated comparison of processing methods for MR phase imaging and QSM: combining array coil signals and phase unwrapping.

Phase imaging benefits from strong susceptibility effects at very high field and the high signal-to-noise ratio (SNR) afforded by multi-channel coils. Combining the information from coils is not trivial, however, as the phase that originates in local field effects (the source of interesting contrast) is modified by the inhomogeneous sensitivity of each coil. This has historically been addressed by referencing individual coil sensitivities to that of a volume coil, but alternative approaches are required for ultra-high field systems in which no such coil is available. An additional challenge in phase imaging is that the phase that develops up to the echo time is "wrapped" into a range of 2π radians. Phase wraps need to be removed in order to reveal the underlying phase distribution of interest. Beginning with a coil combination using a homogeneous reference volume coil - the Roemer approach - which can be applied at 3 T and lower field strengths, we review alternative methods for combining single-echo and multi-echo phase images where no such reference coil is available. These are applied to high-resolution data acquired at 7 T and their effectiveness assessed via an index of agreement between phase values over channels and the contrast-to-noise ratio in combined images. The virtual receiver coil and COMPOSER approaches were both found to be computationally efficient and effective. The main features of spatial and temporal phase unwrapping methods are reviewed, placing particular emphasis on recent developments in temporal phase unwrapping and Laplacian approaches. The features and performance of these are illustrated in application to simulated and high-resolution in vivo data. Temporal unwrapping was the fastest of the methods tested and the Laplacian the most robust in images with low SNR. © 2016 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

An in vivo study of the orientation-dependent and independent components of transverse relaxation rates in white matter.

Diffusion-weighted imaging (DWI) provides information that allows the estimation of white-matter (WM) fibre orientation and distribution, but it does not provide information about myelin density, fibre concentration or fibre size within each voxel. On the other hand, quantitative relaxation contrasts (like the apparent transverse relaxation, R2∗) offer iron and myelin-related contrast, but their dependence on the orientation of microstructure with respect to the applied magnetic field, B0 , is often neglected. The aim of this work was to combine the fibre orientation information retrieved from the DWI acquisition and the sensitivity to microstructural information from quantitative relaxation parameters. The in vivo measured quantitative transverse relaxation maps (R2 and R2∗) were decomposed into their orientation-dependent and independent components, using the DWI fibre orientation information as prior knowledge. The analysis focused on major WM fibre bundles such as the forceps major (FMj), forceps minor (FMn), cingulum (CG) and corticospinal tracts (CST). The orientation-dependent R2 parameters, despite their small size (0-1.5 Hz), showed higher variability across different fibre populations, while those derived from R2∗, although larger (3.1-4.5 Hz), were mostly bundle-independent. With this article, we have, for the first time, attempted the in vivo characterization of the orientation-(in)dependent components of the transverse relaxation rates and demonstrated that the orientation of WM fibres influences both R2 and R2∗ contrasts.

A modulated closed form solution for quantitative susceptibility mapping--a thorough evaluation and comparison to iterative methods based on edge prior knowledge.

The aim of this study is to perform a thorough comparison of quantitative susceptibility mapping (QSM) techniques and their dependence on the assumptions made. The compared methodologies were: two iterative single orientation methodologies minimizing the l2, l1TV norm of the prior knowledge of the edges of the object, one over-determined multiple orientation method (COSMOS) and a newly proposed modulated closed-form solution (MCF). The performance of these methods was compared using a numerical phantom and in-vivo high resolution (0.65 mm isotropic) brain data acquired at 7 T using a new coil combination method. For all QSM methods, the relevant regularization and prior-knowledge parameters were systematically changed in order to evaluate the optimal reconstruction in the presence and absence of a ground truth. Additionally, the QSM contrast was compared to conventional gradient recalled echo (GRE) magnitude and R2* maps obtained from the same dataset. The QSM reconstruction results of the single orientation methods show comparable performance. The MCF method has the highest correlation (corr MCF=0.95, r(2)MCF=0.97) with the state of the art method (COSMOS) with additional advantage of extreme fast computation time. The L-curve method gave the visually most satisfactory balance between reduction of streaking artifacts and over-regularization with the latter being overemphasized when the using the COSMOS susceptibility maps as ground-truth. R2* and susceptibility maps, when calculated from the same datasets, although based on distinct features of the data, have a comparable ability to distinguish deep gray matter structures.