Multicomponent water proton transverse relaxation and T2-discriminated water diffusion in myelinated and nonmyelinated nerve

Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding

Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.

Fiber-orientation independent component of R2* obtained from single-orientation MRI measurements in simulations and a post-mortem human optic chiasm

The effective transverse relaxation rate (R2*) is sensitive to the microstructure of the human brain like the g-ratio which characterises the relative myelination of axons. However, the fibre-orientation dependence of R2* degrades its reproducibility and any microstructural derivative measure. To estimate its orientation-independent part (R2,iso*) from single multi-echo gradient-recalled-echo (meGRE) measurements at arbitrary orientations, a second-order polynomial in time model (hereafter M2) can be used. Its linear time-dependent parameter, β1, can be biophysically related to R2,iso* when neglecting the myelin water (MW) signal in the hollow cylinder fibre model (HCFM). Here, we examined the performance of M2 using experimental and simulated data with variable g-ratio and fibre dispersion. We found that the fitted β1 can estimate R2,iso* using meGRE with long maximum-echo time (TEmax ≈ 54 ms), but not accurately captures its microscopic dependence on the g-ratio (error 84%). We proposed a new heuristic expression for β1 that reduced the error to 12% for ex vivo compartmental R2 values. Using the new expression, we could estimate an MW fraction of 0.14 for fibres with negligible dispersion in a fixed human optic chiasm for the ex vivo compartmental R2 values but not for the in vivo values. M2 and the HCFM-based simulations failed to explain the measured R2*-orientation-dependence around the magic angle for a typical in vivo meGRE protocol (with TEmax ≈ 18 ms). In conclusion, further validation and the development of movement-robust in vivo meGRE protocols with TEmax ≈ 54 ms are required before M2 can be used to estimate R2,iso* in subjects.

A Method for In-Vivo Mapping of Axonal Diameter Distributions in the Human Brain Using Diffusion-Based Axonal Spectrum Imaging (AxSI)

In this paper we demonstrate a generalized and simplified pipeline called axonal spectrum imaging (AxSI) for in-vivo estimation of axonal characteristics in the human brain. Whole-brain estimation of the axon diameter, in-vivo and non-invasively, across all fiber systems will allow exploring uncharted aspects of brain structure and function relations with emphasis on connectivity and connectome analysis. While axon diameter mapping is important in and of itself, its correlation with conduction velocity will allow, for the first time, the explorations of information transfer mechanisms within the brain. We demonstrate various well-known aspects of axonal morphometry (e.g., the corpus callosum axon diameter variation) as well as other aspects that are less explored (e.g., axon diameter-based separation of the superior longitudinal fasciculus into segments). Moreover, we have created an MNI based mean axon diameter map over the entire brain for a large cohort of subjects providing the reference basis for future studies exploring relation between axon properties, its connectome representation, and other functional and behavioral aspects of the brain.

Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.

Measuring compartmental T2-orientational dependence in human brain white matter using a tiltable RF coil and diffusion-T2 correlation MRI

The anisotropy of brain white matter microstructure manifests itself in orientational-dependence of various MRI contrasts, and can result in significant quantification biases if ignored. Understanding the origins of this orientation-dependence could enhance the interpretation of MRI signal changes in development, ageing and disease and ultimately improve clinical diagnosis. Using a novel experimental setup, this work studies the contributions of the intra- and extra-axonal water to the orientation-dependence of one of the most clinically-studied parameters, apparent transverse relaxation T2. Specifically, a tiltable receive coil is interfaced with an ultra-strong gradient MRI scanner to acquire multidimensional MRI data with an unprecedented range of acquisition parameters. Using this setup, compartmental T2 can be disentangled based on differences in diffusional-anisotropy, and its orientation-dependence further elucidated by re-orienting the head with respect to the main magnetic field B→0. A dependence of (compartmental) T2 on the fibre orientation w.r.t. B→0 was observed, and further quantified using characteristic representations for susceptibility- and magic angle effects. Across white matter, anisotropy effects were dominated by the extra-axonal water signal, while the intra-axonal water signal decay varied less with fibre-orientation. Moreover, the results suggest that the stronger extra-axonal T2 orientation-dependence is dominated by magnetic susceptibility effects (presumably from the myelin sheath) while the weaker intra-axonal T2 orientation-dependence may be driven by a combination of microstructural effects. Even though the current design of the tiltable coil only offers a modest range of angles, the results demonstrate an overall effect of tilt and serve as a proof-of-concept motivating further hardware development to facilitate experiments that explore orientational anisotropy. These observations have the potential to lead to white matter microstructural models with increased compartmental sensitivity to disease, and can have direct consequences for longitudinal and group-wise T2- and diffusion-MRI data analysis, where the effect of head-orientation in the scanner is commonly ignored.

Identification of water compartments in spinal cords by 2 H double quantum filtered NMR

In ² H double quantum filtered (DQF) NMR, the various water compartments are characterized by their different residual quadrupolar interactions. The spectral separation between the different signals enables the measurement of the relaxation of each compartment and the magnetization transfer (MT) between them. In the current study, five water compartments were identified in the ² H DQF spectra of porcine spinal cord. The most prominent signal was the pair of satellites with a quadrupolar splitting of about 550 Hz. ² H DQF MRI optimized for the 550 Hz quadrupolar splitting indicated that this signal originated mainly from the white matter and it was assigned to the myelin water. This splitting does not change upon changing the orientation of the spinal cord relative to the magnetic field, indicating a liquid crystalline nature. Another site exhibiting splitting of about 1500 Hz was assigned to collagenous connective tissue. The narrow central peak was assigned to a combination of intra- and inter-axonal water. The assignment of the other two sites is not certain and requires further study. The rates of MT between the various sites were recorded.

Multidimensional correlation MRI

Multidimensional correlation spectroscopy is emerging as a novel MRI modality that is well suited for microstructure and microdynamic imaging studies, especially of biological specimens. Conventional MRI methods only provide voxel-averaged and mostly macroscopically averaged information; these methods cannot disentangle intra-voxel heterogeneity on the basis of both water mobility and local chemical interactions. By correlating multiple MR contrast mechanisms and processing the data in an integrated manner, correlation spectroscopy is able to resolve the distribution of water populations according to their chemical and physical interactions with the environment. The use of a non-parametric, phenomenological representation of the multidimensional MR signal makes no assumptions about tissue structure, thereby allowing the study of microscopic structure and composition of complex heterogeneous biological systems. However, until recently, vast data requirements have confined these types of measurement to non-localized NMR applications and prevented them from being widely and successfully used in conjunction with imaging. Recent groundbreaking advancements have allowed this powerful NMR methodology to be migrated to MRI, initiating its emergence as a promising imaging approach. This review is not intended to cover the entire field of multidimensional MR; instead, it focuses on pioneering imaging applications and the challenges involved. In addition, the background and motivation that have led to multidimensional correlation MR development are discussed, along with the basic underlying mathematical concepts. The goal of the present work is to provide the reader with a fundamental understanding of the techniques developed and their potential benefits, and to provide guidance to help refine future applications of this technology.

Myelin Water Imaging Demonstrates Lower Brain Myelination in Children and Adolescents With Poor Reading Ability

Magnetic resonance imaging (MRI) provides a means to non-invasively investigate the neurological links with dyslexia, a learning disability that affects one’s ability to read. Most previous brain MRI studies of dyslexia and reading skill have used structural or diffusion imaging to reveal regional brain abnormalities. However, volumetric and diffusion MRI lack specificity in their interpretation at the microstructural level. Myelin is a critical neural component for brain function and plasticity, and as such, deficits in myelin may impact reading ability. MRI can estimate myelin using myelin water fraction (MWF) imaging, which is based on evaluation of the proportion of short T2 myelin-associated water from multi-exponential T2 relaxation analysis, but has not yet been applied to the study of reading or dyslexia. In this study, MWF MRI, intelligence, and reading assessments were acquired in 20 participants aged 10–18 years with a wide range of reading ability to investigate the relationship between reading ability and myelination. Group comparisons showed markedly lower MWF by 16–69% in poor readers relative to good readers in the left and right thalamus, as well as the left posterior limb of the internal capsule, left/right anterior limb of the internal capsule, left/right centrum semiovale, and splenium of the corpus callosum. MWF over the entire group also correlated positively with three different reading scores in the bilateral thalamus as well as white matter, including the splenium of the corpus callosum, left posterior limb of the internal capsule, left anterior limb of the internal capsule, and left centrum semiovale. MWF imaging from T2 relaxation suggests that myelination, particularly in the bilateral thalamus, splenium, and left hemisphere white matter, plays a role in reading abilities. Myelin water imaging thus provides a potentially valuable in vivo imaging tool for the study of dyslexia and its remediation.

Data processing in NMR relaxometry using the matrix pencil

The matrix pencil method (MPM) is explored for stable, reproducible data processing in nuclear magnetic resonance (NMR) relaxometry. Data from one-dimensional and two-dimensional relaxometry experiments designed to measure transverse relaxation T₂, longitudinal relaxation T₁, diffusion coefficient D values, and their correlations in a standard olive oil/water mixture serve as a platform available to any NMR spectroscopist to compare the performance of the MPM to the benchmark inverse Laplace transform (ILT). The data from two practical examples, including the drying of a solvent polymer system and the enzymatic digestion of polysialic acid, were also explored with the MPM and ILT. In the cases considered here, the MPM appears to outperform the ILT in terms of resolution and stability in the determination of fundamental constants for complex materials and mixtures.

Retaining information from multidimensional correlation MRI using a spectral regions of interest generator

Multidimensional correlation magnetic resonance imaging (MRI) is an emerging imaging modality that is capable of disentangling highly heterogeneous and opaque systems according to chemical and physical interactions of water within them. Using this approach, the conventional three dimensional MR scalar images are replaced with spatially resolved multidimensional spectra. The ensuing abundance in microstructural and chemical information is a blessing that incorporates a real challenge: how does one distill and refine it into images while retaining its significant components? In this paper we introduce a general framework that preserves the spectral information from spatially resolved multidimensional data. Equal weight is given to significant spectral components at the single voxel level, resulting in a summarized image spectrum. This spectrum is then used to define spectral regions of interest that are utilized to reconstruct images of sub-voxel components. Using numerical simulations we first show that, contrary to the conventional approach, the proposed framework preserves spectral resolution, and in turn, sensitivity and specificity of the reconstructed images. The retained spectral resolution allows, for the first time, to observe an array of distinct [Formula: see text]-[Formula: see text]-[Formula: see text] components images of the human brain. The robustly generated images of sub-voxel components overcome the limited spatial resolution of MRI, thus advancing multidimensional correlation MRI to fulfilling its full potential.

MR-based age-related effects on the striatum, globus pallidus, and thalamus in healthy individuals across the adult lifespan

While numerous studies have used magnetic resonance imaging (MRI) to elucidate normative age-related trajectories in subcortical structures across the human lifespan, there exists substantial heterogeneity among different studies. Here, we investigated the normative relationships between age and morphology (i.e., volume and shape), and microstructure (using the T1-weighted/T2-weighted [T1w/T2w] signal ratio as a putative index of myelin and microstructure) of the striatum, globus pallidus, and thalamus across the adult lifespan using a dataset carefully quality controlled, yielding a final sample of 178 for the morphological analyses, and 162 for the T1w/T2w analyses from an initial dataset of 253 healthy subjects, aged 18-83. In accordance with previous cross-sectional studies of adults, we observed age-related volume decrease that followed a quadratic relationship between age and bilateral striatal and thalamic volumes, and a linear relationship in the globus pallidus. Our shape indices consistently demonstrated age-related posterior and medial areal contraction bilaterally across all three structures. Beyond morphology, we observed a quadratic inverted U-shaped relationship between T1w/T2w signal ratio and age, with a peak value occurring in middle age (at around 50 years old). After permutation testing, the Akaike information criterion determined age relationships remained significant for the bilateral globus pallidus and thalamus, for both the volumetric and T1w/T2w analyses. Our findings serve to strengthen and expand upon previous volumetric analyses by providing a normative baseline of morphology and microstructure of these structures to which future studies investigating patients with various disorders can be compared.

SyMRI detects delayed myelination in preterm neonates

OBJECTIVES

The software "SyMRI" generates different MR contrasts and characterizes tissue properties based on a single acquisition of a multi-dynamic multi-echo (MDME)-FLAIR sequence. The aim of this study was to assess the applicability of "SyMRI" in the assessment of myelination in preterm and term-born neonates. Furthermore, "SyMRI" was compared with conventional MRI.

METHODS

A total of 30 preterm and term-born neonates were examined at term-equivalent age using a standardized MRI protocol. MDME sequence (acquisition time, 5 min, 24 s)-based post-processing was performed using "SyMRI". Myelination was assessed by scoring seven brain regions on quantitative T1-/T2-maps, generated by "SyMRI" and on standard T1-/T2-weighted images, acquired separately. Analysis of covariance (ANCOVA) (covariate, gestational age (GA) at MRI (GAMRI)) was used for group comparison.

RESULTS

In 25/30 patients (83.3%) (18 preterm and seven term-born neonates), "SyMRI" acquisitions were successfully performed. "SyMRI"-based myelination scores were significantly lower in preterm compared with term-born neonates (ANCOVA: T1: F(1, 22) = 7.420, p = 0.012; T2: F(1, 22) = 5.658, p = 0.026). "SyMRI"-based myelination scores positively correlated with GAMRI (T1: r = 0.662, n = 25, p ≤ 0.001; T2: r = 0.676, n = 25, p ≤ 0.001). The myelination scores based on standard MRI did not correlate with the GAMRI. No significant differences between preterm and term-born neonates were detectable.

CONCLUSIONS

"SyMRI" is a highly promising MR technique for neonatal brain imaging. "SyMRI" is superior to conventional MR sequences in the visual detection of delayed myelination in preterm neonates.

KEY POINTS

• By providing multiple MR contrasts, "SyMRI" is a time-saving method in neonatal brain imaging. • Differences concerning the myelination in term-born and preterm infants are visually detectable on T1-/T2-weighted maps generated by "SyMRI". • "SyMRI" allows a faster and more sensitive assessment of myelination compared with standard MR sequences.

Searching for the neurite density with diffusion MRI: Challenges for biophysical modeling

In vivo mapping of the neurite density with diffusion MRI (dMRI) is a high but challenging aim. First, it is unknown whether all neurites exhibit completely anisotropic ("stick-like") diffusion. Second, the "density" of tissue components may be confounded by non-diffusion properties such as T2 relaxation. Third, the domain of validity for the estimated parameters to serve as indices of neurite density is incompletely explored. We investigated these challenges by acquiring data with "b-tensor encoding" and multiple echo times in brain regions with low orientation coherence and in white matter lesions. Results showed that microscopic anisotropy from b-tensor data is associated with myelinated axons but not with dendrites. Furthermore, b-tensor data together with data acquired for multiple echo times showed that unbiased density estimates in white matter lesions require data-driven estimates of compartment-specific T2 values. Finally, the "stick" fractions of different biophysical models could generally not serve as neurite density indices across the healthy brain and white matter lesions, where outcomes of comparisons depended on the choice of constraints. In particular, constraining compartment-specific T2 values was ambiguous in the healthy brain and had a large impact on estimated values. In summary, estimating neurite density generally requires accounting for different diffusion and/or T2 properties between axons and dendrites. Constrained "index" parameters could be valid within limited domains that should be delineated by future studies.

Myelin Water Fraction Imaging of the Brain in Children with Prenatal Alcohol Exposure

BACKGROUND

Prenatal alcohol exposure (PAE) is linked to alterations of cerebral white matter, including volume and nonspecific diffusion magnetic resonance imaging (MRI) indices of microstructure in humans. Some animal models of PAE have demonstrated myelination deficiencies, but myelin levels have not yet been evaluated in individuals with PAE. Multiecho T₂ MRI offers a quantitative method to estimate myelin water fraction (MWF; related to myelin content) noninvasively, which was used here to evaluate brain myelination in children with PAE.

METHODS

Participants with PAE (n = 10, 6 females, mean age 13.9 years, range 7 to 18 years) and controls (n = 14, 11 females, mean age 13.2 years, range 9 to 16 years) underwent 3T MRI of the brain. T₂ images (15 minutes acquisition for 32 echoes) were used to create MWF maps from which mean MWF was measured in 12 regions of interest (ROIs) including 8 in white matter and 4 in deep gray matter.

RESULTS

As expected, across the combined sample, MWF was highest for major white matter tracts such as the internal capsule and genu/splenium of the corpus callosum (10 to 18%) while the caudate and putamen had MWF less than 5%. Mean MWF was similar across 11/12 brain white and gray matter regions for the PAE and control groups (L/R internal capsule, major forceps, putamen, caudate nucleus, L minor forceps, genu and splenium of corpus callosum). In the PAE group, MWF was positively correlated with age in the genu of corpus callosum and right minor forceps, notably 2 frontal tracts.

CONCLUSIONS

Given comparable MRI-derived myelination fraction measures in PAE relative to controls, white matter alterations shown in other imaging studies, such as diffusion tensor imaging, may reflect microstructural anomalies related to axon caliber and density.

Measuring intra-axonal T2 in white matter with direction-averaged diffusion MRI

PURPOSE

To demonstrate how the T₂ relaxation time of intra-axonal water (T_2a ) in white matter can be measured with direction-averaged diffusion MRI.

METHODS

For b-values larger than about 4000 s/mm² , the direction-averaged diffusion MRI signal from white matter is dominated by the contribution from water within axons, which enables T_2a to be estimated by acquiring data for multiple TE values and fitting a mono-exponential decay curve. If given a value of the intra-axonal diffusivity, an extension of the method allows the extra-axonal relaxation time (T_2e ) to be calculated also. This approach was applied to estimate T_2a in white matter for 3 healthy subjects at 3 T, as well as T_2e for a selected set of assumed intra-axonal diffusivities.

RESULTS

The estimated T_2a values ranged from about 50 ms to 110 ms, with considerable variation among white matter regions. For white matter tracts with primarily collinear fibers, T_2a was found to depend on the angle of the tract relative to the main magnetic field, which is consistent with T_2a being affected by magnetic field inhomogeneities arising from spatial differences in magnetic susceptibility. The T_2e values were significantly smaller than the T_2a values across white matter regions for several plausible choices of the intra-axonal diffusivity.

CONCLUSION

The relaxation time for intra-axonal water in white matter can be determined in a straightforward manner by measuring the direction-averaged diffusion MRI signal with a large b-value for multiple TEs. In healthy brain, T_2a is greater than T_2e and varies considerably with anatomical region.

TE dependent Diffusion Imaging (TEdDI) distinguishes between compartmental T2 relaxation times

Biophysical modeling of macroscopic diffusion-weighted MRI signal in terms of microscopic cellular parameters holds the promise of quantifying the integrity of white matter. Unfortunately, even fairly simple multi-compartment models of proton diffusion in the white matter do not provide a unique, biophysically plausible solution. Here we report a nontrivial diffusion MRI signal dependence on echo time (TE) in human white matter in vivo. We demonstrate that such TE dependence originates from compartment-specific T2 values and that it is a promising "orthogonal measure" able to break the degeneracy in parameter estimation, and to yield important relaxation metrics robustly. We thereby enable the precise estimation of the intra- and extra-axonal water T2 relaxation times, which is precluded by a limited signal-to-noise ratio when using multi-echo relaxometry alone.

MRI-based assessment of function and dysfunction in myelinated axons

Repetitive electrical activity produces microstructural alteration in myelinated axons, which may afford the opportunity to noninvasively monitor function of myelinated fibers in peripheral nervous system (PNS)/CNS pathways. Microstructural changes were assessed via two different magnetic-resonance-based approaches: diffusion fMRI and dynamic T2 spectroscopy in the ex vivo perfused bullfrog sciatic nerves. Using this robust, classical model as a platform for testing, we demonstrate that noninvasive diffusion fMRI, based on standard diffusion tensor imaging (DTI), can clearly localize the sites of axonal conduction blockage as might be encountered in neurotrauma or other lesion types. It is also shown that the diffusion fMRI response is graded in proportion to the total number of electrical impulses carried through a given locus. Dynamic T2 spectroscopy of the perfused frog nerves point to an electrical-activity-induced redistribution of tissue water and myelin structural changes. Diffusion basis spectrum imaging (DBSI) reveals a reversible shift of tissue water into a restricted isotropic diffusion signal component. Submyelinic vacuoles are observed in electron-microscopy images of tissue fixed during electrical stimulation. A slowing of the compound action potential conduction velocity accompanies repetitive electrical activity. Correlations between electrophysiology and MRI parameters during and immediately after stimulation are presented. Potential mechanisms and interpretations of these results are discussed.

Delineation of early brain development from fetuses to infants with diffusion MRI and beyond

Dynamic macrostructural and microstructural changes take place from the mid-fetal stage to 2 years after birth. Delineating structural changes of the brain during early development provides new insights into the complicated processes of both typical development and the pathological mechanisms underlying various psychiatric and neurological disorders including autism, attention deficit hyperactivity disorder and schizophrenia. Decades of histological studies have identified strong spatial and functional maturation gradients in human brain gray and white matter. The recent improvements in magnetic resonance imaging (MRI) techniques, especially diffusion MRI (dMRI), relaxometry imaging, and magnetization transfer imaging (MTI) have provided unprecedented opportunities to non-invasively quantify and map the early developmental changes at whole brain and regional levels. Here, we review the recent advances in understanding early brain structural development during the second half of gestation and the first two postnatal years using modern MR techniques. Specifically, we review studies that delineate the emergence and microstructural maturation of white matter tracts, as well as dynamic mapping of inhomogeneous cortical microstructural organization unique to fetuses and infants. These imaging studies converge into maturational curves of MRI measurements that are distinctive across different white matter tracts and cortical regions. Furthermore, contemporary models offering biophysical interpretations of the dMRI-derived measurements are illustrated to infer the underlying microstructural changes. Collectively, this review summarizes findings that contribute to charting spatiotemporally heterogeneous gray and white matter structural development, offering MRI-based biomarkers of typical brain development and setting the stage for understanding aberrant brain development in neurodevelopmental disorders.

Inferring brain tissue composition and microstructure via MR relaxometry

MRI relaxometry is sensitive to a variety of tissue characteristics in a complex manner, which makes it both attractive and challenging for characterizing tissue. This article reviews the most common water proton relaxometry measures, T₁, T₂, and T₂^*, and reports on their development and current potential to probe the composition and microstructure of brain tissue. The development of these relaxometry measures is challenged by the need for suitably accurate tissue models, as well as robust acquisition and analysis methodologies. MRI relaxometry has been established as a tool for characterizing neural tissue, particular with respect to myelination, and the potential for further development exists.

Cerebral White Matter Maturation Patterns in Preterm Infants: An MRI T2 Relaxation Anisotropy and Diffusion Tensor Imaging Study

BACKGROUND AND PURPOSE

Preterm birth is associated with worse neurodevelopmental outcome, but brain maturation in preterm infants is poorly characterized with standard methods. We evaluated white matter (WM) of infant brains at term-equivalent age, as a function of gestational age at birth, using multimodal magnetic resonance imaging (MRI).

METHODS

Infants born very preterm (<32 weeks gestation) and late preterm (33-36 weeks gestation) were scanned at 3 T at term-equivalent age using diffusion tensor imaging (DTI) and T2 relaxometry. MRI data were analyzed using tract-based spatial statistics, and anisotropy of T2 relaxation was also determined. Principal component analysis and linear discriminant analysis were applied to seek the variables best distinguishing very preterm and late preterm groups.

RESULTS

Across widespread regions of WM, T2 is longer in very preterm infants than in late preterm ones. These effects are more prevalent in regions of WM that myelinate earlier and faster. Similar effects are obtained from DTI, showing that fractional anisotropy (FA) is lower and radial diffusivity higher in the very preterm group, with a bias toward earlier myelinating regions. Discriminant analysis shows high sensitivity and specificity of combined T2 relaxometry and DTI for the detection of a distinct WM development pathway in very preterm infants. T2 relaxation is anisotropic, depending on the angle between WM fiber and magnetic field, and this effect is modulated by FA.

CONCLUSIONS

Combined T2 relaxometry and DTI characterizes specific patterns of retarded WM maturation, at term equivalent age, in infants born very preterm relative to late preterm.

Magnetic resonance microdynamic imaging reveals distinct tissue microenvironments

Magnetic resonance imaging (MRI) provides a powerful set of tools with which to investigate biological tissues noninvasively and in vivo. Tissues are heterogeneous in nature; an imaging voxel contains an ensemble of different cells and extracellular matrix components. A long-standing challenge has been to infer the content of and interactions among these microscopic tissue components within a macroscopic imaging voxel. Spatially resolved multidimensional relaxation-diffusion correlation (REDCO) spectroscopy holds the potential to deliver such microdynamic information. However, to date, vast data requirements have mostly relegated these type of measurements to nuclear magnetic resonance applications and prevented them from being widely and successfully used in conjunction with imaging. By using a novel data acquisition and processing strategy in this study, spatially resolved REDCO could be performed in reasonable scanning times with excellent prospects for clinical applications. This new MR imaging framework-which we term "magnetic resonance microdynamic imaging (MRMI)"-permits the simultaneous noninvasive and model-free quantification of multiple subcellular, cellular, and interstitial tissue microenvironments within a voxel. MRMI is demonstrated with a fixed spinal cord specimen, enabling the quantification of microscopic tissue components with unprecedented specificity. Tissue components, such as axons, neuronal and glial soma, and myelin were identified on the basis of their multispectral signature within individual imaging voxels. These tissue elements could then be composed into images and be correlated with immunohistochemistry findings. MRMI provides novel image contrasts of tissue components and a new family of microdynamic biomarkers that could lead to new diagnostic imaging approaches to probe biological tissue alterations accompanied by pathological or developmental changes.

Exploring the origins of echo-time-dependent quantitative susceptibility mapping (QSM) measurements in healthy tissue and cerebral microbleeds

Quantitative susceptibility mapping (QSM) is increasingly used to measure variation in tissue composition both in the brain and in other areas of the body in a range of disease pathologies. Although QSM measurements were originally believed to be independent of the echo time (TE) used in the gradient-recalled echo (GRE) acquisition from which they are derived; recent literature (Sood et al., 2016) has shown that these measurements can be highly TE-dependent in a number of brain regions. In this work we systematically investigate possible causes of this effect through analysis of apparent frequency and QSM measurements derived from data acquired at multiple TEs in vivo in healthy brain regions and in cerebral microbleeds (CMBs); QSM data acquired in a gadolinium-doped phantom; and in QSM data derived from idealized simulated phase data. Apparent frequency measurements in the optic radiations (OR) and central corpus callosum (CC) were compared to those predicted by a 3-pool white matter model, however the model failed to fully explain contrasting frequency profiles measured in the OR and CC. Our results show that TE-dependent QSM measurements can be caused by a failure of phase unwrapping algorithms in and around strong susceptibility sources such as CMBs; however, in healthy brain regions this behavior appears to result from intrinsic non-linear phase evolution in the MR signal. From these results we conclude that care must be taken when deriving frequency and QSM measurements in strong susceptibility sources due to the inherent limitations in phase unwrapping; and that while signal compartmentalization due to tissue microstructure and content is a plausible cause of TE-dependent frequency and QSM measurements in healthy brain regions, better sampling of the MR signal and more complex models of tissue are needed to fully exploit this relationship.

Degeneracy in model parameter estimation for multi-compartmental diffusion in neuronal tissue

The ultimate promise of diffusion MRI (dMRI) models is specificity to neuronal microstructure, which may lead to distinct clinical biomarkers using noninvasive imaging. While multi-compartment models are a common approach to interpret water diffusion in the brain in vivo, the estimation of their parameters from the dMRI signal remains an unresolved problem. Practically, even when q space is highly oversampled, nonlinear fit outputs suffer from heavy bias and poor precision. So far, this has been alleviated by fixing some of the model parameters to a priori values, for improved precision at the expense of accuracy. Here we use a representative two-compartment model to show that fitting fails to determine the five model parameters from over 60 measurement points. For the first time, we identify the reasons for this poor performance. The first reason is the existence of two local minima in the parameter space for the objective function of the fitting procedure. These minima correspond to qualitatively different sets of parameters, yet they both lie within biophysically plausible ranges. We show that, at realistic signal-to-noise ratio values, choosing between the two minima based on the associated objective function values is essentially impossible. Second, there is an ensemble of very low objective function values around each of these minima in the form of a pipe. The existence of such a direction in parameter space, along which the objective function profile is very flat, explains the bias and large uncertainty in parameter estimation, and the spurious parameter correlations: in the presence of noise, the minimum can be randomly displaced by a very large amount along each pipe. Our results suggest that the biophysical interpretation of dMRI model parameters crucially depends on establishing which of the minima is closer to the biophysical reality and the size of the uncertainty associated with each parameter.

In vivo mapping of human spinal cord microstructure at 300mT/m

The ability to characterize white matter microstructure non-invasively has important applications for the diagnosis and follow-up of several neurological diseases. There exists a family of diffusion MRI techniques, such as AxCaliber, that provide indices of axon microstructure, such as axon diameter and density. However, to obtain accurate measurements of axons with small diameters (<5μm), these techniques require strong gradients, i.e. an order of magnitude higher than the 40-80mT/m currently available in clinical systems. In this study we acquired AxCaliber diffusion data at a variety of different q-values and diffusion times in the spinal cord of five healthy subjects using a 300mT/m whole body gradient system. Acquisition and processing were optimized using state-of-the-art methods (e.g., 64-channel coil, template-based analysis). Results consistently show an average axon diameter of 4.5+/-1.1μm in the spinal cord white matter. Diameters ranged from 3.0μm (gracilis) to 5.9μm (spinocerebellar tracts). Values were similar across laterality (left-right), but statistically different across spinal cord pathways (p<10(-5)). The observed trends are similar to those observed in animal histology. This study shows, for the first time, in vivo mapping of axon diameter in the spinal cord at 300mT/m, thus creating opportunities for applications in spinal cord diseases.

White-matter relaxation time and myelin water fraction differences in young adults with autism

BACKGROUND

Increasing evidence suggests that autism is associated with abnormal white-matter (WM) anatomy and impaired brain 'connectivity'. While myelin plays a critical role in synchronized brain communication, its aetiological role in autistic symptoms has only been indirectly addressed by WM volumetric, relaxometry and diffusion tensor imaging studies. A potentially more specific measure of myelin content, termed myelin water fraction (MWF), could provide improved sensitivity to myelin alteration in autism.

METHOD

We performed a cross-sectional imaging study that compared 14 individuals with autism and 14 age- and IQ-matched controls. T 1 relaxation times (T 1), T 2 relaxation times (T 2) and MWF values were compared between autistic subjects, diagnosed using the Autism Diagnostic Interview - Revised (ADI-R), with current symptoms assessed using the Autism Diagnostic Observation Schedule (ADOS) and typical healthy controls. Correlations between T 1, T 2 and MWF values with clinical measures [ADI-R, ADOS, and the Autism Quotient (AQ)] were also assessed.

RESULTS

Individuals with autism showed widespread WM T 1 and MWF differences compared to typical controls. Within autistic individuals, worse current social interaction skill as measured by the ADOS was related to reduced MWF although not T 1. No significant differences or correlations with symptoms were observed with respect to T 2.

CONCLUSIONS

Autistic individuals have significantly lower global MWF and higher T 1, suggesting widespread alteration in tissue microstructure and biochemistry. Areas of difference, including thalamic projections, cerebellum and cingulum, have previously been implicated in the disorder; however, this is the first study to specifically indicate myelin alteration in these regions.

Whole brain myelin mapping using T1- and T2-weighted MR imaging data

Despite recent advancements in MR imaging, non-invasive mapping of myelin in the brain still remains an open issue. Here we attempted to provide a potential solution. Specifically, we developed a processing workflow based on T1-w and T2-w MR data to generate an optimized myelin enhanced contrast image. The workflow allows whole brain mapping using the T1-w/T2-w technique, which was originally introduced as a non-invasive method for assessing cortical myelin content. The hallmark of our approach is a retrospective calibration algorithm, applied to bias-corrected T1-w and T2-w images, that relies on image intensities outside the brain. This permits standardizing the intensity histogram of the ratio image, thereby allowing for across-subject statistical analyses. Quantitative comparisons of image histograms within and across different datasets confirmed the effectiveness of our normalization procedure. Not only did the calibrated T1-w/T2-w images exhibit a comparable intensity range, but also the shape of the intensity histograms was largely corresponding. We also assessed the reliability and specificity of the ratio image compared to other MR-based techniques, such as magnetization transfer ratio (MTR), fractional anisotropy (FA), and fluid-attenuated inversion recovery (FLAIR). With respect to these other techniques, T1-w/T2-w had consistently high values, as well as low inter-subject variability, in brain structures where myelin is most abundant. Overall, our results suggested that the T1-w/T2-w technique may be a valid tool supporting the non-invasive mapping of myelin in the brain. Therefore, it might find important applications in the study of brain development, aging and disease.

Advances in multiple sclerosis and its variants: conventional and newer imaging techniques

Multiple sclerosis (MS) and its variants are inflammatory as well as neurodegenerative diseases that diffusely affect the central nervous system (CNS). There is a poor correlation between traditional imaging findings and symptoms in patients with MS. Current research in conventional magnetic resonance (MR) imaging of MS and related diseases includes optimization of hardware and pulse sequences and the development of automated and semiautomated techniques to measure and quantify disease burden. Advanced nonconventional MR techniques such as diffusion tensor and functional MR imaging probe the changes found in the CNS, and correlate these findings with clinical measures of disease.

Prediction of cartilage compressive modulus using multiexponential analysis of T(2) relaxation data and support vector regression

Evaluation of mechanical characteristics of cartilage by magnetic resonance imaging would provide a noninvasive measure of tissue quality both for tissue engineering and when monitoring clinical response to therapeutic interventions for cartilage degradation. We use results from multiexponential transverse relaxation analysis to predict equilibrium and dynamic stiffness of control and degraded bovine nasal cartilage, a biochemical model for articular cartilage. Sulfated glycosaminoglycan concentration/wet weight (ww) and equilibrium and dynamic stiffness decreased with degradation from 103.6 ± 37.0 µg/mg ww, 1.71 ± 1.10 MPa and 15.3 ± 6.7 MPa in controls to 8.25 ± 2.4 µg/mg ww, 0.015 ± 0.006 MPa and 0.89 ± 0.25MPa, respectively, in severely degraded explants. Magnetic resonance measurements were performed on cartilage explants at 4 °C in a 9.4 T wide-bore NMR spectrometer using a Carr-Purcell-Meiboom-Gill sequence. Multiexponential T2 analysis revealed four water compartments with T2 values of approximately 0.14, 3, 40 and 150 ms, with corresponding weight fractions of approximately 3, 2, 4 and 91%. Correlations between weight fractions and stiffness based on conventional univariate and multiple linear regressions exhibited a maximum r(2) of 0.65, while those based on support vector regression (SVR) had a maximum r(2) value of 0.90. These results indicate that (i) compartment weight fractions derived from multiexponential analysis reflect cartilage stiffness and (ii) SVR-based multivariate regression exhibits greatly improved accuracy in predicting mechanical properties as compared with conventional regression.

Brain development in preterm infants assessed using advanced MRI techniques

Infants who are born preterm have a high incidence of neurocognitive and neurobehavioral abnormalities, which may be associated with impaired brain development. Advanced magnetic resonance imaging (MRI) approaches, such as diffusion MRI (d-MRI) and functional MRI (fMRI), provide objective and reproducible measures of brain development. Indices derived from d-MRI can be used to provide quantitative measures of preterm brain injury. Although fMRI of the neonatal brain is currently a research tool, future studies combining d-MRI and fMRI have the potential to assess the structural and functional properties of the developing brain and its response to injury.

Investigation of outer cortical magnetisation transfer ratio abnormalities in multiple sclerosis clinical subgroups

BACKGROUND

Pathological abnormalities including demyelination and neuronal loss are reported in the outer cortex in multiple sclerosis (MS).

OBJECTIVE

We investigated for in vivo evidence of outer cortical abnormalities by measuring the magnetisation transfer ratio (MTR) in MS patients of different subgroups.

METHODS

Forty-four relapsing-remitting (RR) (mean age 41.9 years, median Expanded Disability Status Scale (EDSS) 2.0), 25 secondary progressive (SP) (54.1 years, EDSS 6.5) and 19 primary progressive (PP) (53.1 years, EDSS 6.0) MS patients and 35 healthy control subjects (mean age 39.2 years) were studied. Three-dimensional (3D) 1×1×1mm(3) T1-weighted images and MTR data were acquired. The cortex was segmented, then subdivided into outer and inner bands, and MTR values were calculated for each band.

RESULTS

In a pairwise analysis, mean outer cortical MTR was lower than mean inner cortical MTR in all MS groups and controls (p<0.001). Compared with controls, outer cortical MTR was decreased in SPMS (p<0.001) and RRMS (p<0.01), but not PPMS. Outer cortical MTR was lower in SPMS than PPMS (p<0.01) and RRMS (p<0.01).

CONCLUSIONS

Lower outer than inner cortical MTR in healthy controls may reflect differences in myelin content. The lowest outer cortical MTR was seen in SPMS and is consistent with more extensive outer cortical (including subpial) pathology, such as demyelination and neuronal loss, as observed in post-mortem studies of SPMS patients.

Reprint of "Quantitative evaluation of brain development using anatomical MRI and diffusion tensor imaging"

The development of the brain is structure-specific, and the growth rate of each structure differs depending on the age of the subject. Magnetic resonance imaging (MRI) is often used to evaluate brain development because of the high spatial resolution and contrast that enable the observation of structure-specific developmental status. Currently, most clinical MRIs are evaluated qualitatively to assist in the clinical decision-making and diagnosis. The clinical MRI report usually does not provide quantitative values that can be used to monitor developmental status. Recently, the importance of image quantification to detect and evaluate mild-to-moderate anatomical abnormalities has been emphasized because these alterations are possibly related to several psychiatric disorders and learning disabilities. In the research arena, structural MRI and diffusion tensor imaging (DTI) have been widely applied to quantify brain development of the pediatric population. To interpret the values from these MR modalities, a "growth percentile chart," which describes the mean and standard deviation of the normal developmental curve for each anatomical structure, is required. Although efforts have been made to create such a growth percentile chart based on MRI and DTI, one of the greatest challenges is to standardize the anatomical boundaries of the measured anatomical structures. To avoid inter- and intra-reader variability about the anatomical boundary definition, and hence, to increase the precision of quantitative measurements, an automated structure parcellation method, customized for the neonatal and pediatric population, has been developed. This method enables quantification of multiple MR modalities using a common analytic framework. In this paper, the attempt to create an MRI- and a DTI-based growth percentile chart, followed by an application to investigate developmental abnormalities related to cerebral palsy, Williams syndrome, and Rett syndrome, have been introduced. Future directions include multimodal image analysis and personalization for clinical application.

Quantitative evaluation of brain development using anatomical MRI and diffusion tensor imaging

The development of the brain is structure-specific, and the growth rate of each structure differs depending on the age of the subject. Magnetic resonance imaging (MRI) is often used to evaluate brain development because of the high spatial resolution and contrast that enable the observation of structure-specific developmental status. Currently, most clinical MRIs are evaluated qualitatively to assist in the clinical decision-making and diagnosis. The clinical MRI report usually does not provide quantitative values that can be used to monitor developmental status. Recently, the importance of image quantification to detect and evaluate mild-to-moderate anatomical abnormalities has been emphasized because these alterations are possibly related to several psychiatric disorders and learning disabilities. In the research arena, structural MRI and diffusion tensor imaging (DTI) have been widely applied to quantify brain development of the pediatric population. To interpret the values from these MR modalities, a "growth percentile chart," which describes the mean and standard deviation of the normal developmental curve for each anatomical structure, is required. Although efforts have been made to create such a growth percentile chart based on MRI and DTI, one of the greatest challenges is to standardize the anatomical boundaries of the measured anatomical structures. To avoid inter- and intra-reader variability about the anatomical boundary definition, and hence, to increase the precision of quantitative measurements, an automated structure parcellation method, customized for the neonatal and pediatric population, has been developed. This method enables quantification of multiple MR modalities using a common analytic framework. In this paper, the attempt to create an MRI- and a DTI-based growth percentile chart, followed by an application to investigate developmental abnormalities related to cerebral palsy, Williams syndrome, and Rett syndrome, have been introduced. Future directions include multimodal image analysis and personalization for clinical application.

Suppressing magnetization exchange effects in stimulated-echo diffusion experiments

Exchange of nuclear magnetization between spin pools, either by chemical exchange or by cross-relaxation or both, has a significant influence on the signal attenuation in stimulated-echo-type pulsed field gradient experiments. Hence, in such cases the obtained molecular self-diffusion coefficients can carry a large systematic error. We propose a modified stimulated echo pulse sequence that contains T2-filters during the z-magnetization store period. We demonstrate, using a common theoretical description for chemical exchange and cross-relaxation, that these filters suppress the effects of exchange on the diffusional decay in that frequent case where one of the participating spin pools is immobile and exhibits a short T2. We demonstrate the performance of this experiment in an agarose/water gel. We posit that this new experiment has advantages over other approaches hitherto used, such as that consisting of measuring separately the magnetization exchange rate, if suitable by Goldman-Shen type experiments, and then correcting for exchange effects within the framework of a two-site exchange model. We also propose experiments based on selective decoupling and applicable in systems with no large T2 difference between the different spin pools.

Temporal phase correction of multiple echo T2 magnetic resonance images

Typically, magnetic resonance imaging (MRI) analysis is performed on magnitude data, and multiple echo T2 data consist of numerous images of the same slice taken with different echo spacing, giving voxel-wise temporal sampling of the noise as the signals decay according to T2 relaxation. Magnitude T2 decay data has Rician distributed noise which is characterized by a change in the noise distribution from Gaussian, through a transitional region, to Rayleigh as the signal to noise ratio decreases with increasing echo time. Non-Gaussian noise distributions may produce errors in the commonly applied non-negative least squares (NNLS) algorithm that is used to assess multiple echo decays for compartmentalized water environments through the creation of T2 distributions. Typically, Gaussian noise is sought by performing spatial-based phase correction on the MRI data however, these methods cannot capitalize on the temporal information available from multiple echo T2 acquisitions. Here we describe a temporal phase correction (TPC) algorithm that utilizes the temporal noise information available in multiple echo T2 acquisitions to put the relevant decay information in the Real portion of the decay data and leave only noise in the Imaginary portion. We apply this TPC algorithm to create real-valued multiple echo T2 data from human subjects measured at 1.5 T. We show that applying TPC causes changes in the T2 distribution estimates; notably the possible resolution of separate extracellular and intracellular water environments, and the disappearance of the commonly labeled cerebrospinal fluid peak, which might be an artefact observed in many previously published multiple echo T2 analyses.

The role of tissue microstructure and water exchange in biophysical modelling of diffusion in white matter

Biophysical models that describe the outcome of white matter diffusion MRI experiments have various degrees of complexity. While the simplest models assume equal-sized and parallel axons, more elaborate ones may include distributions of axon diameters and axonal orientation dispersions. These microstructural features can be inferred from diffusion-weighted signal attenuation curves by solving an inverse problem, validated in several Monte Carlo simulation studies. Model development has been paralleled by microscopy studies of the microstructure of excised and fixed nerves, confirming that axon diameter estimates from diffusion measurements agree with those from microscopy. However, results obtained in vivo are less conclusive. For example, the amount of slowly diffusing water is lower than expected, and the diffusion-encoded signal is apparently insensitive to diffusion time variations, contrary to what may be expected. Recent understandings of the resolution limit in diffusion MRI, the rate of water exchange, and the presence of microscopic axonal undulation and axonal orientation dispersions may, however, explain such apparent contradictions. Knowledge of the effects of biophysical mechanisms on water diffusion in tissue can be used to predict the outcome of diffusion tensor imaging (DTI) and of diffusion kurtosis imaging (DKI) studies. Alterations of DTI or DKI parameters found in studies of pathologies such as ischemic stroke can thus be compared with those predicted by modelling. Observations in agreement with the predictions strengthen the credibility of biophysical models; those in disagreement could provide clues of how to improve them. DKI is particularly suited for this purpose; it is performed using higher b-values than DTI, and thus carries more information about the tissue microstructure. The purpose of this review is to provide an update on the current understanding of how various properties of the tissue microstructure and the rate of water exchange between microenvironments are reflected in diffusion MRI measurements. We focus on the use of biophysical models for extracting tissue-specific parameters from data obtained with single PGSE sequences on clinical MRI scanners, but results obtained with animal MRI scanners are also considered. While modelling of white matter is the central theme, experiments on model systems that highlight important aspects of the biophysical models are also reviewed.

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.

Gradient echo plural contrast imaging--signal model and derived contrasts: T2*, T1, phase, SWI, T1f, FST2*and T2*-SWI

Gradient Echo Plural Contrast Imaging (GEPCI) is a post processing technique that, based on a widely available multiple gradient echo sequence, allows simultaneous generation of naturally co-registered images with various contrasts: T1 weighted, R2*=1/T2* maps and frequency (f) maps. Herein, we present results demonstrating the capability of GEPCI technique to generate image sets with additional contrast characteristics obtained by combing the information from these three basic contrast maps. Specifically, we report its ability to generate GEPCI-susceptibility weighted images (GEPCI-SWI) with improved SWI contrast that is free of T1 weighting and RF inhomogeneities; GEPCI-SWI-like images with the contrast similar to original SWI; T1f images that offer superior GM/WM matter contrast obtained by combining the GEPCI T1 and frequency map data; Fluid Suppressed T2* (FST2*) images that utilize GEPCI T1 data to suppress CSF signal in T2* maps and provide contrast similar to FLAIR T2 weighted images; and T2*-SWI images that combine SWI contrast with quantitative T2* map and offer advantages of visualizing venous structure with hyperintense T2* lesions (e.g. MS lesions). To analyze GEPCI images we use an improved algorithm for combining data from multi-channel RF coils and a method for unwrapping phase/frequency maps that takes advantage of the information on phase evolution as a function of gradient echo time in GEPCI echo train.

Diffusion MRI at 25: exploring brain tissue structure and function

Diffusion MRI (or dMRI) came into existence in the mid-1980s. During the last 25 years, diffusion MRI has been extraordinarily successful (with more than 300,000 entries on Google Scholar for diffusion MRI). Its main clinical domain of application has been neurological disorders, especially for the management of patients with acute stroke. It is also rapidly becoming a standard for white matter disorders, as diffusion tensor imaging (DTI) can reveal abnormalities in white matter fiber structure and provide outstanding maps of brain connectivity. The ability to visualize anatomical connections between different parts of the brain, non-invasively and on an individual basis, has emerged as a major breakthrough for neurosciences. The driving force of dMRI is to monitor microscopic, natural displacements of water molecules that occur in brain tissues as part of the physical diffusion process. Water molecules are thus used as a probe that can reveal microscopic details about tissue architecture, either normal or in a diseased state.

In vivo assessment of myelination by phase imaging at high magnetic field

The present study evaluated the potential of using the phase of T2* weighted MR images to characterize myelination during brain development and pathology in rodents at 9.4 T. Phase contrast correlated with myelin content assessed by histology and suggests that most contrast between white and cortical gray matter is modulated by myelin. Ex vivo experiments showed that gray-white matter phase contrast remains unchanged after iron extraction. In dysmyelinated shiverer mice, phase imaging correlated strongly with myelin staining, showing reduced contrast between white and gray matter when compared to healthy controls. We conclude that high-resolution phase images, acquired at high field, allow assessment of myelination and dysmyelination.

Relaxation time mapping in multiple sclerosis

Several relaxation mapping techniques have been proposed to quantitatively assess disease-related brain tissue changes in multiple sclerosis. Newer developments also account for the distribution of hydrogen protons in different tissue compartments, and therefore provide markers for myelin and macromolecular content. This article will cover the broad spectrum of the pulse sequences and analysis techniques related to this topic that are currently available. Various technical and practical limitations linked with specific approaches will be discussed. These include acquisition time, accuracy and precision, radiofrequency absorption and limited coverage of the brain. Finally, the application of these techniques in the context of multiple sclerosis will be reviewed.

Human brain iron mapping using atlas-based T2 relaxometry

Several in vivo quantitative MRI techniques have been proposed as surrogate measures to map iron content in the human brain. The majority of in vivo quantitative MRI iron mapping methods used the age-dependent iron content data based on postmortem data. In this work, we fused atlas-based human brain volumetry obtained on a large cohort of healthy adults using FreeSurfer with T(2) relaxation time measurements. We provide a brain atlas-based T(2) relaxation time map, which was subsequently used along with published postmortem iron content data to obtain a map of iron content in subcortical and cortical gray matter. We have also investigated the sensitivity of the linear model relating transverse relaxation rate with published iron content to the number of regions used. Our work highlights the challenges encountered on using the simple model along with postmortem data to infer iron content in several brain regions where postmortem iron data are scant (e.g., corpus callosum, amygdale).

Algebraic T2 estimation improves detection of right temporal lobe epilepsy by MR T2 relaxometry

Seizure related abnormalities may be detected with T2 relaxometry, which involves quantitative estimation of T2 values. Accounting for the partial-volume effect of cerebrospinal fluid (CSF) is important, especially for voxel-based relaxometry, VBR. With a mono-exponential decay model, this can be accomplished by including a baseline constant. An algebraic calculation, which accommodates this constant, offers improved T2 estimation speed over the commonly used non-linear fitting approach. Our objective was to compare the algebraic approach against three fitting approaches for the detection of seizure related abnormalities. We tested the performance of the four methods in the presence of noise using simulated data as well as real data acquired at 3 T with a Carr-Purcell-Meiboom-Gill sequence from 45 healthy subjects and 24 patients with confirmed right temporal lobe epilepsy. A quantitative analysis was performed on spatially normalized data by measuring T2 in various regions and with a whole brain tissue segmentation analysis. The detection rate of hippocampal T2 changes in patients was assessed by comparing the regional T2 measurements from each patient against the control data with a z-score threshold of 2.33. The algebraic method yielded high sensitivity for detection of hippocampal abnormalities in the epileptic patients in regional assessment and in follow-up single-subject VBR. This can be attributed to the relatively small variance across healthy subjects and improved precision in the presence of CSF and noise in simulation. In conclusion, the algebraic method is better than fitting based on faster calculation speed and better sensitivity for detecting seizure-related T2 changes.

Atlas-based investigation of human brain tissue microstructural spatial heterogeneity and interplay between transverse relaxation time and radial diffusivity

Microstructural metrics obtained using magnetic resonance imaging (MRI) such as transverse relaxation time and radial diffusivity have been used as in vivo markers of human brain tissue integrity. Considering the sensitivity of these parameters to some common biophysical contributors and their structural and spatial heterogeneity, we hypothesized that strong inter and intra-regional associations exist between these variables providing evidence to possible interplay between transverse relaxation time and radial diffusivity. To validate our hypothesis we obtained high resolution anatomical T1-weighted data and fused it with T2-relaxometry and diffusion tensor imaging (DTI) data on a cohort of healthy adults. The anatomical data were parcellated using FreeSurfer and then coaligned and fused with the T2 and DTI maps. Our data reveal some association between transverse relaxation and radial diffusivity that may help toward the interpretation and modeling of the biophysical contributors to the measured MRI metrics.

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.