Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Assessing White Matter Microstructure in Brain Regions with Different Myelin Architecture Using MRI

  • Samuel Groeschel ,

    Contributed equally to this work with: Samuel Groeschel, Gisela E. Hagberg

    Samuel.groeschel@med.uni-tuebingen.de

    Affiliation University Children's Hospital Tübingen, Germany

  • Gisela E. Hagberg ,

    Contributed equally to this work with: Samuel Groeschel, Gisela E. Hagberg

    Affiliations High Field Magnetic Resonance, Max-Planck Institute for Biological Cybernetics, Tübingen, Germany, Biomedical Magnetic Resonance, University Hospital Tübingen, Germany

  • Thomas Schultz,

    Affiliation Institute of Computer Science, University of Bonn, Germany

  • Dávid Z. Balla,

    Affiliation Department Physiology of Cognitive Processes, Max Planck Institute for Biological Cybernetics, Tübingen, Germany

  • Uwe Klose,

    Affiliation Department of Diagnostic and Interventional Neuroradiology, University Hospital, Tübingen, Germany

  • Till-Karsten Hauser,

    Affiliation Department of Diagnostic and Interventional Neuroradiology, University Hospital, Tübingen, Germany

  • Thomas Nägele,

    Affiliation Department of Diagnostic and Interventional Neuroradiology, University Hospital, Tübingen, Germany

  • Oliver Bieri,

    Affiliation Radiological Physics, University of Basel, Basel, Switzerland

  • Thomas Prasloski,

    Affiliation University of British Columbia, Vancouver, Canada

  • Alex L. MacKay,

    Affiliation University of British Columbia, Vancouver, Canada

  • Ingeborg Krägeloh-Mann,

    Affiliation University Children's Hospital Tübingen, Germany

  • Klaus Scheffler

    Affiliations High Field Magnetic Resonance, Max-Planck Institute for Biological Cybernetics, Tübingen, Germany, Biomedical Magnetic Resonance, University Hospital Tübingen, Germany

Abstract

Objective

We investigate how known differences in myelin architecture between regions along the cortico-spinal tract and frontal white matter (WM) in 19 healthy adolescents are reflected in several quantitative MRI parameters that have been proposed to non-invasively probe WM microstructure. In a clinically feasible scan time, both conventional imaging sequences as well as microstructural MRI parameters were assessed in order to quantitatively characterise WM regions that are known to differ in the thickness of their myelin sheaths, and in the presence of crossing or parallel fibre organisation.

Results

We found that diffusion imaging, MR spectroscopy (MRS), myelin water fraction (MWF), Magnetization Transfer Imaging, and Quantitative Susceptibility Mapping were myelin-sensitive in different ways, giving complementary information for characterising WM microstructure with different underlying fibre architecture. From the diffusion parameters, neurite density (NODDI) was found to be more sensitive than fractional anisotropy (FA), underlining the limitation of FA in WM crossing fibre regions. In terms of sensitivity to different myelin content, we found that MWF, the mean diffusivity and chemical-shift imaging based MRS yielded the best discrimination between areas.

Conclusion

Multimodal assessment of WM microstructure was possible within clinically feasible scan times using a broad combination of quantitative microstructural MRI sequences. By assessing new microstructural WM parameters we were able to provide normative data and discuss their interpretation in regions with different myelin architecture, as well as their possible application as biomarker for WM disorders.

Introduction

The assessment of white matter (WM) microstructure generally aims at detecting, characterising and quantifying suspected myelin pathology. These changes in the myelination of WM pathways have been shown to directly affect brain connectivity resulting in disease conditions and behavioural variation [1,2]. In a clinical context, quantitative information about WM microstructure can not only help in identifying the aetiology of a disease but can also provide biomarkers for evaluating the dynamics of the disease course and/or for evaluating treatment effects [35].

Advances in magnetic resonance imaging (MRI) techniques make the quantification of putative WM biomarkers increasingly robust and standardised, thus enabling direct comparisons between different scanners and study sites [68]. Moreover, it is likely that such advances will enable the characterisation of histology-like features of the white matter in the near future [912].

There are several established MR methods available that are sensitive to different aspects of WM microstructure, and hence of interest for the study of WM pathology. Most methods are based on the notion that magnetic properties of water surrounded by densely packed myelin are profoundly different from water moving and exchanging more freely inside the cytosol of neuronal and glial cell bodies or water in the extra-cellular space. MR-derived parameters such as the longitudinal (T1) and (effective) transverse (T2*) T2 relaxation times, the proton density, the magnetic susceptibility, the magnetization transfer (MT) effect, as well as different measures related to water diffusion are sensitive to microstructure at sub-voxel length scales. Great care was taken to select the best candidate methods that are sensitive to different aspects of the microstructure, yet can be combined in a measurement protocol of limited total duration for future use in our Children’s hospital.

A fundamental MR parameter is the longitudinal relaxation time, T1, which determines the time that a perturbed nuclear spin distribution needs to get back to equilibrium. It was shown early that in brain tissue this process is presumably driven by the lipids in the myelin sheaths [13].A more recent study demonstrated that variations in brain T1-maps can be predicted by a linear combination of two other MR parameters, magnetization transfer and effective transverse relaxation [14]. Indeed, the main contributor to T1-contrastin the brain seems to be the chemical exchange of bound and free water pools, which was parameterised by exploiting the magnetization transfer (MT) effect, and we therefore, chose to directly assess this MR parameter, rather than measure T1 in white matter tracts.

In Magnetization Transfer (MT) imaging, the effect of saturating the signal from extremely short-lived water pools, for instance originating in the lipid bilayer of myelin, is indirectly observed by comparing the MR-visible water with and without saturation transfer pulses applied on the bound pool. Indeed, the bound water pool has a very short T2 and therefore gives rise to a resonance line which is much broader than the MR-visible resonance and can therefore easily be targeted by saturation RF pulses [15]. For a full quantification of this biophysical effect, complex modelling of different water pools and quantification of the different T1 times of the pools is necessary. This endeavour can be facilitated by the dynamic analysis of increasing levels of saturation achieved by repetitive MT-pulses [16]. As an alternative to a full quantification of the multi-compartmental processes of MT, the assessment of the magnetization transfer ratio (MTR) has been shown to be clinically useful, e.g. in WM disorders [4,17]. In the present study, the fast bSSFP MT method is used, which has been shown to be a robust, sensitive and very fast measure of MTR [7,18,19].

Other methods capable of directly capturing the MR signal originating in water pools trapped between myelin sheaths are based on the transverse relaxation time, T2. The indirect effect of myelin on these water pools provokes a loss in MR-signal coherence, which is more rapid than for water pools located within the axons or in the cytoplasm. Measurements of the T2 difference between pools have proven highly relevant for characterising white matter. Pioneered by the group of MacKay [20], the separation of multiple water components by T2 relaxometry is a relatively straightforward and well-established measure of myelin water that has been shown to correlate well with histology and to remain unaffected by inflammation or unspecific changes in the overall water content [21]. Basically, by applying a train of 180° refocusing pulses (the so-called Carr-Purcell-Meibom-Gill experiment or CPMG, after its inventors) the method separates components with a long T2 time arising from the mobile water pool and the myelin trapped water with short T2 times, being below 30ms at 3T. A few caveats that have been dealt with over the years relate to the presence of stimulated echoes arising due to imperfect 180° pulses across the slice thickness and magnetization transfer effects arising from the excitation of spins in adjacent slices in case of multi-slice 2D imaging. The solution to this has been to take into account the presence of stimulated echoes in the modelling of the MR signal decay [22], on one hand, and to perform single-slice measurements with increased bandwidth for the refocusing pulses, on the other hand. These improvements have rendered the method more robust. In our study, we used the 2D CPMG approach in order to focus our investigation on a single slice with the same spatial coverage as the MR spectroscopy sequence.

Proton MR spectroscopy (MRS) offers an insight into the brain chemistry and therefore information about the composition of certain cellular components. However, the metabolic spectrum of myelin is complex with many overlapping functional groups [23] as well as compounds with short T2-relaxation times that hampers their detection [24]. Still, certain metabolite peaks have been found to be enriched in neurons or myelin membranes. Therefore, MRS gives additional valuable information for the characterisation of WM, although the technique most likely is inadequate for directly quantifying myelination [4]. The metabolite N-acetylaspartate (NAA) has been described to be enriched in neuronal structures [25,26], whereas Choline-containing compounds (Cho) were found to be associated with myelin membranes [27,28]. These metabolites show some regional variation throughout the brain, probably due to the different underlying tissue composition [29,30].

A MR parameter which is highly sensitive to the presence of myelin is the magnetic susceptibility, which basically is a measure of how well the static magnetic field can penetrate the tissue. In tissues that are more paramagnetic than the surrounding tissue, the local field experienced will be slightly higher, while the reverse happens in more diamagnetic areas. The phase of the MR image can be seen as a ‘finger print’ of the local magnetic field distribution, and regional variations in susceptibility can thus be probed through this parameter yielding higher phase values in paramagnetic areas (although the definition may differ between scanners, [31]). The white and grey matter show a strong difference in magnetic susceptibility leading to contrast differences in phase images, first observed at 7T [32]. To overcome limits posed by the orientation dependence and non-local components of the phase effect, methods that allow conversion of the phase maps into quantitative susceptibility maps (QSM) have been developed [33,34]. QSM gives a measure of the average magnetic susceptibility distribution in every unit volume (voxel)and, hence, can be used to map regional changes in tissue microstructure [35]. For clinical use, QSM is an emerging method, and for this reason we included it in a small cohort of subjects in our study. Since a multi-echo gradient echo sequence was used for this purpose, we also evaluated the presence of R2* effects.

Diffusion-weighted imaging (DWI)sequences are sensitized to detect the motion of water molecules, and hence tissue properties, in the range of micrometers [36,37]. The most commonly applied model in order to quantify diffusion in the brain is the tensor model [38], forming the basis of diffusion tensor imaging (DTI). By estimating the eigen-values from this tensor, different measures that characterises diffusion along (= axial diffusivity) and perpendicular (= radial diffusivity) to the fibre pathways can be assessed. More commonly, all eigen-values are combined into the dimensionless metric fractional anisotropy (FA);while the trace of the tensor is proportional to the mean diffusivity (MD) per voxel [39,40]. While DTI metrics are useful for characterizing coherent fibre bundles, the technique is of limited use for WM areas containing crossing fibres [39,4144]. As the diffusion tensor can only detect one major diffusion direction per voxel, radial and axial diffusivities can be interpreted as diffusivities across and along the bundle only in voxels with one (parallel) fibre population, which has been shown to be true only for 10% of WM voxels [45], limiting the interpretation of DTI-parameter considerably [41]. As a valuable extension of DTI, diffusion kurtosis imaging (DKI) has been introduced [46]. By measuring the non-gaussian part of the diffusion signal, it increases the sensitivity to pathological changes [47]. Still, besides the lack of specificity to cellular components [48], DKI is also affected by the problem of crossing fibres, where the axial and radial kurtosis is difficult to interpret [49]. To overcome this major limitation a scan protocol that includes multiple diffusion-weightings can be used to derive additional DWI parameters separating the highly restricted water diffusion inside the axon from the less hindered diffusion in the extra-axonal space [50]. Such approaches have been shown to reveal the underlying fibre architecture and provide microstructural parameters with improved sensitivity to pathological changes and better specificity to histological tissue properties than the more classical DTI-based metrics [5153]. While it has been argued that a strict compartmentalization of the tissue might be an oversimplification of modelling the diffusion signal [54], several promising models have been introduced and validated by histology [9,10,55,56], although some of them require very long scan times. One of the recently introduced methods that can be achieved within clinically feasible scan times, is neurite orientation dispersion and density imaging (NODDI) [57]. Although used with certain assumptions, it attempts to disentangle and quantify intra-axonal water diffusion and fibre dispersion, making it valuable for the purpose of our study.

White matter microstructure is known to vary across the brain. Especially fibres of the cortico-spinal tract (CST) are known to have different tissue properties, compared to e.g. more frontal WM pathways, containing larger pyramidal fibres with increased axonal diameter and larger myelin sheaths [58]. Fig 1 illustrates the expected underlying differences in fibre architecture in three selected WM regions, according to histological data [5861]. In the posterior limb of the internal capsule (PLIC), tightly packed and highly myelinated parallel fibre bundles of the CST with larger axon diameters run through and proceed towards the primary motor area. They pass the level of the centrum semiovale (CS), where these large-diameter axons are still highly myelinated, but run through WM regions with crossing fibres from commissural (corpus callosum) and association pathways (superior longitudinal fasciculus) [62]. On the other hand, the frontal WM region is characterised by a mixture of crossing fibres that interconnect between neocortical areas and connect to other brain areas further downstream, showing thinner myelin sheaths and smaller axon diameters as well as larger extracellular spaces (Fig 1).

thumbnail
Fig 1. Illustration of different underlying fibre architecture in the three measurement regions as found by histology (e.g. [5861]: Cortico-spinal tract (CST) fibres include those with larger axon diameter and thicker myelin sheaths, which show more parallel arrangements at the level of the posterior limb of internal capsule (PLIC) and more crossing fibres at the level of the centrum semiovale (CS).

On the other hand, the frontal white matter (WM) region also contains crossing fibres, however, with thinner myelin sheaths and lesser axon diameter.

https://doi.org/10.1371/journal.pone.0167274.g001

Such microstructural differences have been shown to be detectable bya series of different MRI sequences [6365]. Furthermore, even along the CST, the regional variation in fibre architecture (crossing/bending/single fibre populations) influence microstructural MRI parameters, resulting e.g. in different DTI-derived values at the level of the internal capsule and the centrum semiovale [62,65].

These areas were thus considered appropriate to probe WM microstructure in regions with different underlying fibre architecture, with our combined set of techniques. Due to its functional relevance and known unique microstructural features [58,6365], the CST was assessed both at the level of the PLIC and the centrum semiovale and compared to a WM region in the frontal lobe for the whole brain techniques. The described MRI parameters were thus investigated, probing their ability to detect the known differences in microstructure in the selected regions. Furthermore, we aimed to provide and evaluate this myelin-sensitive MRI protocol to non-invasively probe WM microstructure in a clinical context, thereby providing normative data and interpretation of new microstructural WM parameters aimed for studies of juvenile and adult forms of leukodystrophies [66]. Therefore a matching age range was selected, covering the adolescence and young adulthood.

Methods

Subjects

Nineteen typically developing healthy adolescents (median age 14.9 years, range 9–32 years, 10 females) were recruited as controls for an ongoing study about juvenile leukodystrophies, covering the age range of the expected patients. The age distribution between the two genders were kept similar with males having a median age of 14.9 years (range 9–32 years) and females having a median age of 14.8 years (range of 10–24 years), therefore reducing the impact of a gender bias for this study. A standardised gross motor function test was performed in all subjects (GMFM 88) and demonstrated normal gross motor function in all subjects (100%). In addition, the Wechsler Intelligence Scale for children or adults was performed. The total IQ values in this group ranged from 103 to 130 (mean 114.6, SD 6.89). The study was approved by the local ethics committee. All volunteers and/or their legal guardians gave their written informed consent for participation in the study.

MRI acquisition

An MRI sequence protocol was established on a 3Tesla wide-bore scanner (“Skyra”, Siemens Healthcare GmbH, Erlangen, Germany) equipped with a 2-channel body transmit coil and a 32-channel head receiver coil, to allow data acquisition within a total (net) examination time of 42 minutes. Conventional MRI sequences consisted of 4 acquisitions, two T2-weighted turbo-spin echo scans acquired with axial (40 slices, acquisition time (TA) = 1:59min) and sagittal (30 slices, TA = 1:59min) image orientations (with the following sequence parameters: repetition time(TR)/echo time (TE) = 10810/84ms, field of view (FoV) = 190x192mm2, matrix size = 384x384, slice thickness = 3mm, reconstructed to a voxel size of 0.5x0.5x3mm3, 2 averages, grappa factor of 2); an axial T2-weighted fluid attenuated inversion recovery (FLAIR) sequence (TA = 1:30min, TR/TE = 8800/85ms, inversion time (TI) = 2480ms, flip angle = 150°, FoV = 172×230 mm2, matrix = 134×256, slice thickness 3mm, 40 slices, voxel-size: 0.9x0.9x3mm3, grappa factor of 3); and a T1-weighted Magnetization-Prepared Rapid Acquisition Gradient Echo (MPRAGE) sequence (TA 3:27min, TE/TI/TR = 4.11/900/2300ms;FA = 9°, FoV = 256×256 mm2,matrix size 256×256, 176 slices, voxel-size = 1 mm isotropic, grappa factor of 3). All axial images were acquired in the AC-PC orientation.

In addition to the conventional sequences, we also acquired data for the processing of several quantitative parameter maps. The diffusion-weighted (DW) images were acquired using a high angular resolution twice-refocused echo planar imaging (EPI) read-out with two shells (TA = 15:10min, 64 directions with b = 2000s/mm2; 30 directions with b = 700s/mm2, TE/TR = 89/9100ms, FOV = 192×192mm2; 96×96 matrix, slice thickness = 2mm, 50 interleaved contiguous slices, voxel size = 2×2×2mm³, bandwidth 1796 Hz/Px, grappa factor of 2). For EPI distortion correction, an additional b = 0 image (TA = 0:09min) was acquired with the same resolution, FoV and readout bandwidth, but with reversal of the direction of acquisition along the phase-encode (PE) axis [67].

For MR spectroscopy, two approaches were used: a chemical-shift imaging (CSI) and a single voxel (SV) sequence in order to evaluate the sensitivity and robustness of the sequences. In nine subjects, a CSI sequence with the same axial orientation as the axial T2-weighted image positioned above the lateral ventricles(TA = 3:57min, TR/TE = 1600/135ms, FoV = 160x160mm2; matrix size 12x12, interpolated to 32x32; slice thickness = 12mm, voxel size 5x5x15mm³). In another ten subjects, a SV PRESS-sequence was performed in two locations within the level of the CSI slice: one in the left frontal WM and one in the parietal WM at the central region using the central sulcus as a landmark (where the CST was expected to be located) (TA = 1:49min, TR/TE: 1600/135ms, 64 averages, voxel size 15x15x15mm3, bandwidth 1000 Hz, vector size 1024, 4 preparation scans, water suppression with bandwidth 35Hz, flip angle 90°).

Using the same localization as the CSI sequence, a Carr-Purcell-Meiboom-Gill (CPMG) sequence (TA = 4:17min, TR = 3s, TE = 10ms to 320ms in steps of 10ms, FoV = 192x192mm2; matrix size 128x128, slice thickness = 5mm, 1 slice, voxel size = 0.75×0.75×5mm3) was acquired, modified to have shorter RF pulse durations and increased bandwidth of the 180° RF pulse than the standard vendor sequence. B1-mapping for each transmit channel was performed and combined using the manufacturer’s standard approach and a voxel size that matched the CPMG sequence (TA = 0:10min, TR/TE 5000/1.72ms, nominal flip angle 8°, FoV = 192x192mm2; matrix size 128x128, slice thickness = 5mm, 1 slice, voxel size = 0.75×0.75×5mm3).

Two sets of balanced steady state free precession(balanced SSFP) 3D images were acquired using slab selective RF pulses (flip angle of 20°, FOV = 256x256mm2, matrix size = 192×192, 144 slices, isotropic voxel size1.3mm3, reconstructed to0.65x0.65x1.3mm³). One 3D data set was acquired with a long duration of the excitation pulse:1500μs (TA = 3:05min, TE/TR = 2.12/4.23ms) to minimize Magnetization Transfer (MT) effects while the second data set had a RF-pulse duration of200μs (TE/TR = 1.47/2.93ms) for the MT-weighted scans.

Images for quantitative susceptibility mapping (QSM) and R2* mapping were acquired with a 3D flow-compensated gradient-echo sequence at eight (4.5/10/15/21/27/33/39/46ms) different echo times using a TR of 50ms, flip angle of 15° and a voxel size of 1x1x2mm3 (N = 5).

Image processing

Motion and eddy-current correction was applied to the DWI data set as described before [68]. In addition, correction for EPI distortions was performed for improved alignment with anatomical images prior to calculation of the diffusion tensor [69,70]. The data was corrected for variation in signal intensity during scanning using the b0 images [71]. Maps of mean diffusivity (MD), fractional anisotropy (FA), as well as parallel (or axial, λ1) and perpendicular (or radial, λ(2+3)*0.5) eigenvalues were calculated using the standard log-linear least squares fit of the tensor model [38,72]. For the calculation of the DTI parameters, only the low b-value shell was used (30 directions with b = 700s/mm2) for better comparability to “classic” DTI studies.

The diffusional kurtosis model was fitted using quadratic cone programming. Similar to Tabesh et al. [73], we performed a linear least squares fit with constraints on diffusivities and kurtosis values that improve robustness to noise. Our implementation is based on the cvxopt package (cvxopt.org) and imposes a semi-definiteness constraint on diffusivities (d≥0.1 μm2/ms), since lower values are implausible and would render Tabesh’s upper bound on kurtosis ineffective. Maps of mean kurtosis (MK), radial kurtosis (RK), and axial kurtosis (AK) were computed based on the equations provided by Tabesh et al. [73].

Furthermore, NODDI-derived parameters, including intracellular volume fraction (ICVF, neurite density), orientation dispersion (ODI), and free water fraction (ISO) were calculated based on the publicly available implementation of the NODDI three—compartment model [57]. Even though the ICVF parameter has been referred to as a “measurable water fraction” rather than a volume fraction, as a reminder of the fact that water inside the myelin does not contribute to the DWI signal [52], we will keep the more widely used terminology.

MRI Spectra were analysed quantitatively using the LCModel [74]. Absolute concentrations of N-acetylaspartate and N-acetylaspartylglutamate (NAA), creatine and phosphocreatine (Cr), as well as choline containing compounds (Cho) were calculated. For this study, only ratios of NAA/Cr and Cho/Cr were further analysed.

MWF was fitted voxel-wise [22], in 2 steps. First, all parameters including the flip angle were fitted. In the second step, the flip angle was fixed to the voxel-specific value measured by the B1 scaled to match the median value found in the first step. Number of fitted T2 points: 32; Chi2 regularization of 1.01.

The MT ratio (MTR) map was calculated in percentage units from the co-registered non-MT and MT images using the following equation: (non-MT–MT) / non-MT.

For QSM the frequency maps were generated by using the phase unwrapping and multi-echo fitting tool from the MEDI-Toolbox [75] and the BET masking tool from FSL [76]. Background dipole field modulations were removed with RESHARP [77] and the local dipole inversion process was performed by using the fast magnitude-weighted L1-regularization technique with total variation penalty [78]. No rescaling of the QSM values was undertaken.R2* mapping was performed by assuming a mono-exponential signal decay. Non-linear least square fitting based on the Levenberg-Marquard algorithm was made pixel-wise for the square of the MRI signal decay, to minimize the influence of pixels and time-points with low signal.

For each subject, conventional images, MWF, MTR, QSM, and R2*were all co-registered (rigid body transformation using normalized mutual information) and interpolated to the subject’s FA map, therefore it was made sure that all image modalities per subject were transformed into the same space and image dimensions for further analysis.Images were processed and analysed using MRtrix (version 0.3.12, www.github.com/MRtrix3/mrtrix3 [79]), Matlab (version R2014b, www.mathworks.com) and FSL tools (version 5, http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FSL, [80]).

Data analysis

For data analysis, a region of interest (ROI) approach was chosen focusing on the three regions illustrated in Fig 1:Two regions along the cortico-spinal tract (CST) were selected, one at the level of the internal capsule, CST-PLIC, and another one located in the centrum semiovale, CST-CS, at the level where the CSI slice was placed, just above the ventricles (see Fig 2). The CST was tracked in each individual using manual delineation of the pre-central gyrus and the descending pathways at the level of the ponsas a seed as described previously [62]. Probabilistic tractography was performed using constrained spherical deconvolution using MRtrix [79,81]. A third ROI was defined in the frontal WM at the level of the MWF and CSI slice, anterior to the CST ROI (Fig 2). In order to avoid any bias of the ROI definition related to the local MR signal intensity in the modalities of interest, ROIs were drawn onto the co-registered conventional T1-weighted images with the CST overlaid, where signal intensity throughout the WM is homogeneous.

thumbnail
Fig 2. MRI protocol for this study.

T1-weighted MPRAGE image (T1, left) in two different axial planes of a 17y-old typically developing subject indicating the measurement ROIs guided by the cortico-spinal tract fibre s outlined in red: frontal WM (in blue), cortico-spinal tract at the level of the centrum semiovale (in green) and at the level of the posterior limb of the internal capsule (in yellow). Other image parameters included in the protocol were T2-weighted axial TSE (T2), mean diffusivity (MD), fractional anisotropy (FA), intracellular volume fraction (ICVF), orientation dispersion (ODI), Magnetization Transfer Ratio (MTR), Myelin Water Fraction (MWF), Chemical Shift Imaging (CSI)–MR Spectroscopy, and quantitative susceptibility mapping (QSM, the white matter fibers have been highlighted by rendering diamagnetic effects bright, and paramagnetic dark).

https://doi.org/10.1371/journal.pone.0167274.g002

For each ROI mean values were calculated across the two hemispheres, since the objective of this study was to compare between ROIs, and not to explore the hemispheric lateralization of the parameters. Boxplots of quantitative image parameters were displayed and (paired t-test) statistical analysis was conducted by paired t-test comparing the results between the ROI locations. As this study with a relatively small sample size was considered exploratory, the p-values for these analyses were not corrected for multiple comparisons and can therefore be regarded as being descriptive.

Results

Fig 2 shows the MRI results obtained by different modalities from one 17 year-old male subject. In Table 1, the results for the ROI measurements are listed and reproduced in boxplots (Figs 3 and 4). Most microstructural MRI parameter showed significant differences between the investigated WM regions.

thumbnail
Fig 3.

Boxplots of diffusion-weighted imaging parameters(Units of diffusivities are given in μm2/ms) included a) mean diffusivity (MD), b) fractional anisotropy (FA), c) orientation dispersion (ODI), d) intracellular volume fraction (ICVF),e) mean kurtosis (MK), f) radial kurtosis (RK), and g) axial kurtosis (AK). Parameters were measured in frontal white matter (WM) and two parameters within the cortico-spinal tract (CST): at the level of the posterior limb of internal capsule (PLIC) and at the level of the centrum semiovale (CS), see also Fig 1.

https://doi.org/10.1371/journal.pone.0167274.g003

thumbnail
Fig 4.

Boxplots of quantitative parameters included a) ratio of N-acetylaspartate and N-acetylaspartylglutamate (NAA) to creatine and phosphocreatine (Cr) both for chemical shift imaging (CSI) and single voxel (SV) measurements, b) ratio of choline containing compounds (Cho) to Cr both for CSI and SV, c) myelin water fraction (MWF), d) magnetization transfer ratio (MTR), e) quantitative susceptibility mapping (QSM), and f) R2*. Parameters were measured in frontal white matter (WM) and two parameters within the cortico-spinal tract (CST): at the level of the posterior limb of internal capsule (PLIC) and at the level of the centrum semiovale (CS), see also Fig 1.

https://doi.org/10.1371/journal.pone.0167274.g004

thumbnail
Table 1. Measurement results in the three regions of interest: the frontal WM, the CST at the level of the centrum semiovale (CST-CS), and at the level of the posterior limb of the internal capsule (CST-PLIC).

https://doi.org/10.1371/journal.pone.0167274.t001

Diffusion imaging revealed higher MD in the frontal WM ROI than in the CST-CS (p<0.01) and CST-PLIC (p<0.01) ROIs without differences in MD between CST-CS and PLIC (see Fig 3A). DTI-derived parameters showed significantly higher parallel and lower radial eigenvalues in the CST-PLIC region, resulting in a higher FA (p<0.001) (Fig 3B). Similarly, smaller but still significant differences with higher radial (p<0.01)and higher parallel (p = 0.066) eigenvalues were found in the frontal WM than in the CST-CS regions, resulting in a lower FA in the frontal WM region compared to CST-CS (p = 0.003).

DKI parameters revealed highest values for MK and RK in the CST-PLIC (p<0.001), clearly lower values in CST-CS, and even lower ones in the frontal WM (p<0.001). Differences in AK between the ROIs were less pronounced but still significant with lowest AK in CST-PLIC and highest in CST-CS (p<0.001).

The NODDI parameters also showed significant regional differences. The neurite density (ICVF, Fig 3D), was highest in the CST-PLIC, clearly lower in the CST-CS (p<0.001) and lowest in the frontal WM (p<0.001). In contrast, the ODI parameter (Fig 3C) was similar in the frontal WM and the CST-CS, but was significantly lower in the CST-PLIC (p<0.001).

The MR Spectroscopy (Fig 4A and 4B) values revealed higher NAA/Cr ratios in the CST-CS compared to the frontal WM, regardless of acquisition methodology, with slightly lower standard deviations, and hence improved significance for the CSI approach (p<0.001) than for SV measurements (p = 0.011). On the other hand, both methods showed lower Cho/Cr ratios in the CST-CS region than in the frontal WM ROI (p<0.01).

Quantitative values obtained by Myelin water fraction imaging (MWF, Fig 4C) were found to be significantly higher in the CST-CS region compared to the frontal WM (p<0.001). The Magnetization Transfer Ration (MTR), however, was slightly higher in frontal WM compared to the CST-CS region, while significantly higher MTR values were found in the CST-PLIC than in the frontal WM and CST-CS region (p<0.001). Similarly, we found no significant differences between the CST-CS and the frontal WM in the QSMs while a clearly higher diamagnetic susceptibility was measured in the CST-PLIC (p<0.001).On the other hand, R2* was not significantly different between the three ROIs.

Discussion

In this study we were able to probe WM microstructure by MRI in brain regions with different underlying fibre architecture. In a clinically feasible scan time, various myelin-sensitive quantitative MR methods could be acquired. The results demonstrate that the underlying (normal) WM fibre architecture influences significantly and differently the behaviour of several quantitative MRI parameters.

High values were observed in all myelin-sensitive MR parameters assessed along the CST, including MTR, QSM, MWF, FA, MK, and ICVF, reflecting the unique microstructural characteristics of the CST with thicker myelin sheaths and tightly packed and partly larger fibres (Fig 1) [5861]. Aside with this finding the MD values were low, reflecting a high degree of diffusion restriction and little free extracellular water diffusion due to aligned and tightly packed (myelin) membranes. On the other hand, in the frontal WM region, we observed higher MD and Cho/Cr levels, and lower FA, NAA/Cr, ICVF, and MWF levels, when compared to the more parietal WM of the centrum semi-ovale, where the CST fibres pass through. Interestingly, we found no differences in ODI and QSM between these two areas, suggesting that these tissue parameters are less influenced by the different proportions of intra- to extra-cellular compartments. ODI and QSM are known to reflect microstructural directionality, which is similar in the two crossing fibre regions (Fig 1).

In the following we will more critically assess and interpret our findings for the different MR modalities.

DWI parameters

The cortico-spinal tract is known to contain fibres with larger axon diameter and myelin thickness due to long descending fibres mainly from the primary motor cortex [58,82]. Since intra-axonal water diffusion is restricted mainly by cell membranes and myelin sheaths, a thicker myelin in the CST leads to a decrease of measured apparent diffusivities radial to the main fibre direction, and to an increase in Fractional Anisotropy. Accordingly, in the frontal WM region water diffusion was less hindered as indicated by a higher MD and a lower FA. Also, MK was highly increased in the PLIC-CST (mainly due to increased radial kurtosis), in line with previous reports [83] and the association of diffusion kurtosis with myelination [84,85]. At the same time, higher intra-axonal water diffusion along the CST leads to higher axial diffusivity, further contributing to the high FA in the CST-PLIC. Here, axial kurtosis was lower compared to crossing fibre regions of the centrum semiovale, similar to previous observations [47].

The microstructural differences between these areas were also evident from the ICVF parameter, in line with results shown recently [11]. Clearly higher ICVF values were found in both CST-regions, compared to the frontal WM, most likely reflecting the increased intra-axonal water diffusion of the larger axons in the CST. Although it has to be noted that ICVF alone might not sufficiently be able to disentangle the influences of axonal density and myelination [52].

When interpreting diffusion parameters, the dispersion of fibre orientations isknown to be a critical factor, limiting the interpretation of FA or eigenvalues as measures for fibre integrity in crossing fibre regions [43,44,86]. As a result, parallel fibre arrangements, indicated by a low ODI, like in the PLIC, resulted in high FA, based on low radial and high parallel diffusivity as modelled by the diffusion tensor. These results are in line with previous diffusion studies of the CST, which showed high anisotropy at the level of the PLIC and low anisotropy at the level of the centrum semiovale [62,64,65]. On the other hand, in voxels with more complex fibre orientations, like crossings and fanning in the frontal WM and centrum semiovale level of the CST, corresponding FA and eigenvalues are highly influenced by the inhomogeneity of fibre orientations, making them difficult to interpret [44,62]. The same holds true for the diffusion kurtosis measures. While the mean kurtosis has been reported to be sensitive to pathological WM changes, and correlated with myelination and maturation of WM [84,85], it is difficult to interpret axial and radial kurtosis in regions of crossing fibres [49].

We found that the ICVF parameter was less influenced by the spread of fibre orientation, better reflecting the integrity of WM pathways in crossing fibre regions, as was recently shown by others and underlined by histology [10]. Interestingly, a recent study also showed ICVF to be more sensitive than FA in detecting brain maturation and myelination [87,88]. The same has been found in a comprehensive recent study [41], where a more specific axonal density parameter was found to be highly correlated with myelin content, whether in single or crossing fibre regions, different to tensor parameters (like FA), which were limited by the effect of crossing fibres.

While the assessment of conventional DTI parameters requires less scan time, the use of a high angular resolution sequence with different b-values allows a more detailed analysis of the diffusion signal per voxel [39]. Both crossing fibres and different diffusion compartments can be investigated by adding only a few more minutes of scan time. The additional data enables the calculation of other microstructural parameter maps, like in our study. Other recently developed and promising DWI parameters include apparent fibre density [89], hindrance modulated orientational anisotropy [90],axon diameter [91,92], neurite density and dispersion [57,93], permeability [94], tortuosity [48], intra-axonal water fraction [48], cellularity [95], vascularity and cell size [96], myelin density [10], per-axon diffusion coefficients [97] and others. Several of these models are clinically practicable, and a more detailed comparison is a potential topic for a more specialized work.

However, notes of caution have been raised regarding the (over)simplification of compartmentalizing the tissue by modelling the diffusion [54]. Also, see Burcaw et al [98] and Fieremans et al [99] regarding the impossibility of mapping inner diameters on today’s clinical scanners in humans. Furthermore, we note that, due to unavoidable assumptions and simplifications, even so-called microstructural parameters do not allow us to unambiguously deduce from the MRI data alone which exact differences in tissue microstructure led to the observed differences in diffusion. The NODDI model used in our study is no exception to this. In particular, in order to infer the reported parameters, NODDI assumes fixed values for the diffusivities in its three compartments. This means that values of the ICVF parameter might be biased estimates of the true measurable intra-cellular water fraction, and differences in ICVF might, in part, reflect the model’s inability to adapt to regional differences in diffusivities [52,100]. However, the main goal of our current work is to empirically identify clinically feasible MR-based parameters that allow us to reliably distinguish between regions with known differences in myelin architecture. We found ICVF to be a promising candidate, even though, as for most available alternatives, the mechanism behind it is not fully known.

Magnetic Resonsance Spectroscopy

Also, MRS parameters (NAA and Cho) were found to be different in frontal compared to the CST region within the centrum semiovale. These regional differences are well documented in the literature [29,30,101], but their interpretation is not straightforward. While NAA is supposed to be a neuronal marker [25,26], Cho is thought to reflect cell membrane turnover [27,28]. Therefore, higher NAA and lower Choline in CST-CS might indicate a higher ratio of glial to axonal cell compartments in the frontal lobe. On the other hand, early studies have shown the involvement of NAA in lipid biosynthesis [102], and may therefore be indicative of a reduced availability due to on-going myelination in frontal areas, or an increased supply of metabolites in the CST. With regard to preferred measurement methodology, the CSI acquisition can be regarded as a more informative approach due to a greater spatial coverage in less acquisition time, at least at clinical magnetic field strengths where the influence of B1 inhomogeneity and the chemical shift artefact on spectral localization is acceptably low. In addition, we observed a greater precision in terms of less variance for CSI than SV, in line with the greater efficiency of the image-based approach.

Myelin Water Fraction

CPMG based imaging of the myelin water fraction is of interest since it correlates well with histology and yields a parameter which remains unaffected by inflammation or unspecific changes in the overall water content [20,21]. We found this parameter to be increased in the CST area compared to the more frontal regions of the WM. This confirms previous results from MWF measurements of the CST [63]. Whether this is indeed due to increased myelin content, leading to a decreased exchange of water between the extra and intracellular compartments in presence of thicker myelin sheaths, remains to be elucidated. A limitation of our study is posed by the use of the 2D CPMG approach, which did not allow whole brain coverage. Therefore we have no data about the CST located within the PLIC. Nevertheless, this approach allowed us to obtain robust measurements within a relatively short acquisition time (~4min) in a slice which could be matched to the slice used for MRS. Indeed, the draw-back of whole brain, multi-spin-echo based observations of the short T2 component is the long measurement times of 15min or more. This has triggered the quest for alternative, faster, gradient echo-based methods, albeit at the cost of more complex data analysis [103]. This approach assesses the myelin water fraction by constraining fits of the measured MRI data to a three-compartment model of the tissue, including magnetization exchange between the different water pools. The potential bias introduced by the modelling per se thus somewhat hampers this approach. Indeed, it has to be shown to yield MWF values that are several times greater than those furnished by the classical CPMG approach [104].

Magnetization Transfer Ratio

MTR was found to be increased in PLIC and not different between frontal WM and the CST region located within the centrum semiovale. A previous study also showed higher MTR in PLIC compared to CST in the centrum semiovale [65]. Recently, it was shown that hypomyelinating disorders are accompanied by a pattern of MTR, that closely match histopathological findings with less MT-effects in areas lacking myelin, and increased MTR in areas with more myelin [4,17]. Nevertheless, MTR remains a relative measure, and depends on the exact implementation of the MT-pulses and how these play out in presence of B1 inhomogeneities related to the design of the transmit RF-coil, besides the MR sequence parameters (TR) and the magnetic field strength per se [105107]. These factors hence need particular attention. In this respect recent advances based on the balanced SSFP sequence hold promise. For instance this approach has been shown to be a robust, sensitive and very fast measure of MTR which generates values that are comparable across multiple measurement sites [7,18,19]. In the present study, we therefore chose to include the fast balanced SSFP MT method as an integral part of our white matter protocol.

Quantitative Susceptibility Mapping and R2*

The QSM maps showed results that were similar to previous studies [108,109], with very high diamagnetic values in the PLIC, caused potentially by tightly packed and myelinated tracts leading to some form of anisotropy of susceptibility. Indeed, the residual orientation dependence of the major white matter fibre tracts in QSM images is puzzling. Most likely, these variations originate in a myelin-driven anisotropy of the susceptibility, with a greater (diamagnetic) susceptibility variation across the myelinated axons than along them [110112]. Such phenomena make the QSM calculations more intricate, and may necessitate the introduction of a susceptibility tensor [113][112].and corrections based on the Generalized Lorentzian Tensor Approach (GLTA) rather than the widely used sphere of Lorentz concept [112]. These recent developments and biophysical insights call for the development of validated software tools for the scientific community in the near future. Although we have not developed and used such tools for our data, we did look into the possibility that the obtained data show orientation dependence. For this purpose we investigated the dependence of the frequency maps on the fibre orientation obtained from the DTI data. According to the GLTA theory, cylindrical cavities should not perturb the field when oriented along it, while in case of a perpendicular orientation additional corrections are called for. In view of this, we investigated the voxel-wise variation of the frequency shift maps against the square of the sine of the angle between the external field and the orientation of the diffusion vector. In our limited set of subjects we found significant variation of the frequency shift with fibre orientation consistent with the GLTA in three out of five subjects (data not shown). This kind of analysis could not explain the high QSM values observed in the PLIC area, and therefore other explanations in terms of the g-ratio, fibre density or averaging of effects in fibres with different orientation residing inside the voxels on the susceptibility must be investigated in the future. Such approaches, although more exact, would likely require that images are acquired for different head orientations, and are not (yet) applicable in the clinic, especially if additional MRI modalities are to be acquired. The simpler approach is based on MRI images acquired for a single head orientation, and the final results rely on the selected method for quantification. Several different flavours for QSM quantification based on the sphere of Lorentz approach have been proposed, and these differ mainly by how unwanted background signals are dealt with and which type of mathematical constraints are put on the solution to the inversion problem [114118].

In contrast to the QSM data, we did find significant orientation dependent effects of R2*, that varied as the fourth power of the sine of the angle between the main DTI orientation and the external field, in accordance with previous reports [111,119,120]. Again the PLIC area stood out in terms of greater than expected values, suggesting consistent effects of a thicker myelin sheet across MRI modalities (QSM, MTR and R2*). On the other hand, in ANT and CST that have similar characteristics in terms of crossing fibers, no significant differences in R2* were found. In future studies, more refined models that incorporate a greater number of microcompartments that give rise to multi-exponential decay may be required[121].

Towards combining quantitative sequences for the assessment of WM pathology

This study shows that different quantitative MRI parameters might add complementary information for the assessment of WM microstructure. Depending on the underlying tissue architecture of the region of interest, the research question and the suspected type and distribution of WM pathology, a combination of WM MRI parameters might be optimal for the quantification and characterisation of WM microstructure [12].

Furthermore, when studying WM pathology, it is important to have appropriate control data, as MR parameters vary throughout the brain. Not only the cortico-spinal tract [6365] as demonstrated in our study, but also callosal fibres [11] and other major fibre bundles in the brain have been shown by MR parameters to vary in their microstructural properties [10,122], confirmed by the findings from anatomical ex-vivo histology studies in humans [60,123]. Several of the MR parameters could be validated with histopathology [10,21,124,125].

Moving towards in vivo histology, studies attempt to assess axon diameter, myelin content or measure the g-ratio [10,11,122], which is the ratio between axon and total fibre diameter [126].Advances in this field of MRI based in vivo histology are promising [12]. For example, over the years several studies have shown that even intra-cortical myelin can be probed through the T1-shortening caused by myelin, giving rise to MR signal differences along the cortical rim [127129]. In white matter fibres, T1 undergoes profound changes in development [130]. Recently, it was shown in adults that variations in T1 can be explained in terms of a combination of other quantitative MR parameters, regardless whether grey or white matter areas are considered [14]. The main contributor to T1 is the magnetization transfer effect, MT, followed by water content and iron, which was assessed in terms of the effective transverse relaxation rate. Therefore, we arguably chose to directly assess these MR parameters, rather than measure T1 in white matter tracts. Since a substantial gain in sensitivity towards microstructure can be obtained by separating the short from the long T1 components [131,132] future studies may prove that this parameter is actually a sensitive, if not specific, measure for assessing myelin.

A critical and limiting factor for MRI protocols used in a clinical setting is the total duration of the exam. Often, trade-offs are made in terms of lower spatial resolution, especially when acquiring a combination of sequences or when subjects do not tolerate long scan times due to their age or brain dysfunction. Choosing carefully the best combination of several sequences allowing a variety of MRI parameters to be analysed [8], and new hardware [133] or acquisition strategies [134,135] will certainly help to overcome this current limitation.

Conclusion

For the first time, a large number of modern microstructural MRI sequences have been systematically compared on a single well-defined cohort. Within a clinically feasible scan time we were able to assess and characterise WM microstructure using a combination of microstructural MRI sequences. By assessing new quantitative WM parameters we were able to provide normative data and discuss their interpretation in regions with different myelin architecture. The results demonstrate that the underlying normal WM fibre architecture influences significantly and differently the behaviour of several quantitative MRI parameters. In tightly packed parallel fibre bundles high values were observed in all myelin-sensitive parameters, including MTR, QSM, FA, and ICVF, and lower MD and ODI, reflecting more diffusion restriction and less free extracellular water diffusion due to aligned and tightly packed (myelin) membranes. On the other hand, in the frontal WM with underlying crossing fibres architecture, higher MD and Cho, and lower FA, NAA, ICVF, and MWF were observed compared to the more parietal WM, where CST fibres with larger axon diameters and thicker myelin sheaths pass through. These normative data are essential when using MR parameters as biomarker for WM disorders.

Acknowledgments

This study was funded by the Medical Faculty of the University of Tübingen, Germany. We also acknowledge support by Deutsche Forschungsgemeinschaft (DFG) and Open Access Publishing Fund of the University of Tübingen, Germany. T.S. was funded by DFG grant SCHU 3040/1-1; G.E.H. by the ministry of Science, Research and the Arts of Baden-Württemberg (Az: 32-771-8-1504.12/1/1).

Author Contributions

  1. Conceptualization: SG GEH TS IKM.
  2. Formal analysis: SG GEH TS DZB.
  3. Funding acquisition: SG IKM KS.
  4. Investigation: SG GEH TS UK TKH.
  5. Methodology: SG GEH TS DZB UK OB TP ALM.
  6. Project administration: SG TN IKM KS.
  7. Software: SG GEH TS DZB UK OB TP ALM.
  8. Supervision: KS IKM.
  9. Writing – original draft: SG GEH TS.
  10. Writing – review & editing: SG GEH TS DZB UK TKH TN OB TP ALM IKM KS.

References

  1. 1. Catani M, Ffytche DH. The rises and falls of disconnection syndromes. Brain. 2005;128: 2224–2239. pmid:16141282
  2. 2. Johansen-Berg H. Behavioural relevance of variation in white matter microstructure. Curr Opin Neurol. 2010;23: 351–358. pmid:20581685
  3. 3. Barkhof F, Calabresi PA, Miller DH, Reingold SC. Imaging outcomes for neuroprotection and repair in multiple sclerosis trials. NatRevNeurol. 2009;5: 256–266.
  4. 4. Pouwels PJW, Vanderver A, Bernard G, Wolf NI, Dreha-Kulczewksi SF, Deoni SCL, et al. Hypomyelinating leukodystrophies: translational research progress and prospects. Ann Neurol. 2014;76: 5–19. pmid:24916848
  5. 5. Horsfield MA, Jones DK. Applications of diffusion-weighted and diffusion tensor MRI to white matter diseases—a review. NMR Biomed. 2002;15: 570–577. pmid:12489103
  6. 6. Deoni SC, Williams SCR, Jezzard P, Suckling J, Murphy DGM, Jones DK. Standardized structural magnetic resonance imaging in multicentre studies using quantitative T1 and T2 imaging at 1.5 T. Neuroimage. 2008;40: 662–671. pmid:18221894
  7. 7. Gloor M, Scheffler K, Bieri O. Intrascanner and interscanner variability of magnetization transfer-sensitized balanced steady-state free precession imaging. Magn Reson Med. 2011;65: 1112–1117. pmid:21413076
  8. 8. Weiskopf N, Suckling J, Williams G, Correia MM, Inkster B, Tait R, et al. Quantitative multi-parameter mapping of R1, PD(*), MT, and R2(*) at 3T: a multi-center validation. Front Neurosci. 2013;7: 95. pmid:23772204
  9. 9. Jelescu IO, Zurek M, Winters KV, Veraart J, Rajaratnam A, Kim NS, et al. In vivo quantification of demyelination and recovery using compartment-specific diffusion MRI metrics validated by electron microscopy. NeuroImage. 2016;132: 104–114. pmid:26876473
  10. 10. Sepehrband F, Clark KA, Ullmann JFP, Kurniawan ND, Leanage G, Reutens DC, et al. Brain tissue compartment density estimated using diffusion-weighted MRI yields tissue parameters consistent with histology. Hum Brain Mapp. 2015;36: 3687–3702. pmid:26096639
  11. 11. Stikov N, Campbell JSW, Stroh T, Lavelée M, Frey S, Novek J, et al. In vivo histology of the myelin g-ratio with magnetic resonance imaging. NeuroImage. 2015;118: 397–405. pmid:26004502
  12. 12. Weiskopf N, Mohammadi S, Lutti A, Callaghan MF. Advances in MRI-based computational neuroanatomy: from morphometry to in-vivo histology. Curr Opin Neurol. 2015;28: 313–322. pmid:26132532
  13. 13. Koenig SH, Brown RD, Spiller M, Lundbom N. Relaxometry of brain: why white matter appears bright in MRI. Magn Reson Med. 1990;14: 482–495. pmid:2355830
  14. 14. Callaghan MF, Helms G, Lutti A, Mohammadi S, Weiskopf N. A general linear relaxometry model of R1 using imaging data. Magn Reson Med. 2015;73: 1309–1314. pmid:24700606
  15. 15. Forsén S, Hoffman RA. Study of Moderately Rapid Chemical Exchange Reactions by Means of Nuclear Magnetic Double Resonance. The Journal of Chemical Physics. 1963;39: 2892.
  16. 16. Helms G, Hagberg GE. In vivo quantification of the bound pool T1 in human white matter using the binary spin-bath model of progressive magnetization transfer saturation. Phys Med Biol. 2009;54: N529–540. pmid:19904029
  17. 17. Dreha-Kulaczewski SF, Brockmann K, Henneke M, Dechent P, Wilken B, Gärtner J, et al. Assessment of myelination in hypomyelinating disorders by quantitative MRI. J Magn Reson Imaging. 2012;36: 1329–1338. pmid:22911904
  18. 18. Garcia M, Gloor M, Bieri O, Wetzel SG, Radue E-W, Scheffler K. MTR variations in normal adult brain structures using balanced steady-state free precession. Neuroradiology. 2011;53: 159–167. pmid:20480154
  19. 19. Gloor M, Scheffler K, Bieri O. Quantitative magnetization transfer imaging using balanced SSFP. Magn Reson Med. 2008;60: 691–700. pmid:18727085
  20. 20. MacKay A, Whittall K, Adler J, Li D, Paty D, Graeb D. In vivo visualization of myelin water in brain by magnetic resonance. Magn Reson Med. 1994;31: 673–677. pmid:8057820
  21. 21. Laule C, Leung E, Lis DKB, Traboulsee AL, Paty DW, MacKay AL, et al. Myelin water imaging in multiple sclerosis: quantitative correlations with histopathology. Mult Scler. 2006;12: 747–753. pmid:17263002
  22. 22. Prasloski T, Mädler B, Xiang Q-S, MacKay A, Jones C. Applications of stimulated echo correction to multicomponent T2 analysis. Magn Reson Med. 2012;67: 1803–1814. pmid:22012743
  23. 23. Wilhelm MJ, Ong HH, Wehrli SL, Li C, Tsai P-H, Hackney DB, et al. Direct magnetic resonance detection of myelin and prospects for quantitative imaging of myelin density. Proc Natl Acad Sci USA. 2012;109: 9605–9610. pmid:22628562
  24. 24. Kreis R, Ernst T, Ross BD. Absolute Quantitation of Water and Metabolites in the Human Brain. II. Metabolite Concentrations. Journal of Magnetic Resonance, Series B. 1993;102: 9–19.
  25. 25. Birken DL, Oldendorf WH. N-acetyl-L-aspartic acid: a literature review of a compound prominent in 1H-NMR spectroscopic studies of brain. Neurosci Biobehav Rev. 1989;13: 23–31. pmid:2671831
  26. 26. Moffett JR, Namboodiri AMA. Preface: a brief review of N-acetylaspartate. Adv Exp Med Biol. 2006;576: vii–xiii. pmid:16802700
  27. 27. Brenner RE, Munro PM, Williams SC, Bell JD, Barker GJ, Hawkins CP, et al. The proton NMR spectrum in acute EAE: the significance of the change in the Cho:Cr ratio. Magn Reson Med. 1993;29: 737–745. pmid:8350716
  28. 28. Frahm J, Hanefeld F. Localized proton magnetic resonance spectroscopy of cerebral metabolites. Neuropediatrics. 1996;27: 64–69. pmid:8737820
  29. 29. Horská A, Kaufmann WE, Brant LJ, Naidu S, Harris JC, Barker PB. In vivo quantitative proton MRSI study of brain development from childhood to adolescence. J Magn Reson Imaging. 2002;15: 137–143. pmid:11836768
  30. 30. Pouwels PJ, Frahm J. Regional metabolite concentrations in human brain as determined by quantitative localized proton MRS. Magn Reson Med. 1998;39: 53–60. pmid:9438437
  31. 31. Hagberg GE, Welch EB, Greiser A. The sign convention for phase values on different vendor systems: definition and implications for susceptibility-weighted imaging. Magn Reson Imaging. 2010;28: 297–300. pmid:19699598
  32. 32. Duyn JH, van Gelderen P, Li T-Q, de Zwart JA, Koretsky AP, Fukunaga M. High-field MRI of brain cortical substructure based on signal phase. Proc Natl Acad Sci USA. 2007;104: 11796–11801. pmid:17586684
  33. 33. Schweser F, Deistung A, Reichenbach JR. Foundations of MRI phase imaging and processing for Quantitative Susceptibility Mapping (QSM). Z Med Phys. 2016;26: 6–34. pmid:26702760
  34. 34. Schäfer A, Wharton S, Gowland P, Bowtell R. Using magnetic field simulation to study susceptibility-related phase contrast in gradient echo MRI. NeuroImage. 2009;48: 126–137. pmid:19520176
  35. 35. He X, Yablonskiy DA. Biophysical mechanisms of phase contrast in gradient echo MRI. Proc Natl Acad Sci USA. 2009;106: 13558–13563. pmid:19628691
  36. 36. Concha L. A macroscopic view of microstructure: Using diffusion-weighted images to infer damage, repair, and plasticity of white matter. Neuroscience. 2013;
  37. 37. Beaulieu C. The basis of anisotropic water diffusion in the nervous system—a technical review. NMR Biomed. 2002;15: 435–455. pmid:12489094
  38. 38. Basser PJ, Mattiello J, LeBihan D. Estimation of the effective self-diffusion tensor from the NMR spin echo. J Magn Reson B. 1994;103: 247–254. pmid:8019776
  39. 39. Tournier J-D, Mori S, Leemans A. Diffusion tensor imaging and beyond. Magn Reson Med. 2011;65: 1532–1556. pmid:21469191
  40. 40. Basser PJ, Jones DK. Diffusion-tensor MRI: theory, experimental design and data analysis—a technical review. NMR Biomed. 2002;15: 456–467. pmid:12489095
  41. 41. De Santis S, Drakesmith M, Bells S, Assaf Y, Jones DK. Why diffusion tensor MRI does well only some of the time: variance and covariance of white matter tissue microstructure attributes in the living human brain. Neuroimage. 2014;89: 35–44. pmid:24342225
  42. 42. Farquharson S, Tournier J-D, Calamante F, Fabinyi G, Schneider-Kolsky M, Jackson GD, et al. White matter fiber tractography: why we need to move beyond DTI. J Neurosurg. 2013;118: 1367–1377. pmid:23540269
  43. 43. Jones DK. Challenges and limitations of quantifying brain connectivity in vivo with diffusion MRI. Imaging Med. 2010;2: 341–355.
  44. 44. Wheeler-Kingshott CA, Cercignani M. About “axial” and “radial” diffusivities. Magn ResonMed. 2009;61: 1255–1260.
  45. 45. Jeurissen B, Leemans A, Tournier J-D, Jones DK, Sijbers J. Investigating the prevalence of complex fiber configurations in white matter tissue with diffusion magnetic resonance imaging. Hum Brain Mapp. 2013;34: 2747–2766. pmid:22611035
  46. 46. Jensen JH, Helpern JA, Ramani A, Lu H, Kaczynski K. Diffusional kurtosis imaging: the quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn Reson Med. 2005;53: 1432–1440. pmid:15906300
  47. 47. Cheung MM, Hui ES, Chan KC, Helpern JA, Qi L, Wu EX. Does diffusion kurtosis imaging lead to better neural tissue characterization? A rodent brain maturation study. NeuroImage. 2009;45: 386–392. pmid:19150655
  48. 48. Fieremans E, Jensen JH, Helpern JA. White matter characterization with diffusional kurtosis imaging. Neuroimage. 2011;58: 177–188. pmid:21699989
  49. 49. Ankele M, Schultz T. Quantifying Microstructure in Fiber Crossings with Diffusional Kurtosis. In: Navab N, Hornegger J, Wells WM, Frangi AF, editors. Medical Image Computing and Computer-Assisted Intervention—MICCAI 2015. Cham: Springer International Publishing; 2015. pp. 150–157. Available: http://link.springer.com/10.1007/978-3-319-24553-9_19
  50. 50. Assaf Y, Basser PJ. Composite hindered and restricted model of diffusion (CHARMED) MR imaging of the human brain. NeuroImage. 2005;27: 48–58. pmid:15979342
  51. 51. Adluru G, Gur Y, Anderson JS, Richards LG, Adluru N, DiBella EVR. Assessment of white matter microstructure in stroke patients using NODDI. Conf Proc IEEE Eng Med Biol Soc. 2014;2014: 742–745. pmid:25570065
  52. 52. Jelescu IO, Veraart J, Adisetiyo V, Milla SS, Novikov DS, Fieremans E. One diffusion acquisition and different white matter models: How does microstructure change in human early development based on WMTI and NODDI? Neuroimage. 2015;107: 242–256. pmid:25498427
  53. 53. Billiet T, Mädler B, D’Arco F, Peeters R, Deprez S, Plasschaert E, et al. Characterizing the microstructural basis of “unidentified bright objects” in neurofibromatosis type 1: A combined in vivo multicomponent T2 relaxation and multi-shell diffusion MRI analysis. Neuroimage Clin. 2014;4: 649–658. pmid:24936416
  54. 54. Kiselev VG, Il’yasov KA. Is the “biexponential diffusion” biexponential? Magnetic Resonance in Medicine. 2007;57: 464–469. pmid:17326171
  55. 55. Kaden E, Kelm ND, Carson RP, Does MD, Alexander DC. Multi-compartment microscopic diffusion imaging. NeuroImage. 2016;139: 346–359. pmid:27282476
  56. 56. Panagiotaki E, Schneider T, Siow B, Hall MG, Lythgoe MF, Alexander DC. Compartment models of the diffusion MR signal in brain white matter: A taxonomy and comparison. NeuroImage. 2012;59: 2241–2254. pmid:22001791
  57. 57. Zhang H, Schneider T, Wheeler-Kingshott CA, Alexander DC. NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage. 2012;61: 1000–1016. pmid:22484410
  58. 58. Yagishita A, Nakano I, Oda M, Hirano A. Location of the corticospinal tract in the internal capsule at MR imaging. Radiology. 1994;191: 455–460. pmid:8153321
  59. 59. Caminiti R, Ghaziri H, Galuske R, Hof PR, Innocenti GM. Evolution amplified processing with temporally dispersed slow neuronal connectivity in primates. Proceedings of the National Academy of Sciences. 2009;106: 19551–19556.
  60. 60. Liewald D, Miller R, Logothetis N, Wagner H-J, Schüz A. Distribution of axon diameters in cortical white matter: an electron-microscopic study on three human brains and a macaque. Biol Cybern. 2014;108: 541–557. pmid:25142940
  61. 61. Perge JA, Niven JE, Mugnaini E, Balasubramanian V, Sterling P. Why do axons differ in caliber? J Neurosci. 2012;32: 626–638. pmid:22238098
  62. 62. Groeschel S, Tournier J-D, Northam GB, Baldeweg T, Wyatt J, Vollmer B, et al. Identification and interpretation of microstructural abnormalities in motor pathways in adolescents born preterm. Neuroimage. 2014;87: 209–219. pmid:24185027
  63. 63. Russell-Schulz B, Laule C, Li DKB, MacKay AL. What causes the hyperintense T2-weighting and increased short T2 signal in the corticospinal tract? Magn Reson Imaging. 2013;31: 329–335. pmid:22921941
  64. 64. Hervé P-Y, Cox EF, Lotfipour AK, Mougin OE, Bowtell RW, Gowland PA, et al. Structural properties of the corticospinal tract in the human brain: a magnetic resonance imaging study at 7 Tesla. Brain Struct Funct. 2011;216: 255–262. pmid:21409416
  65. 65. Reich DS, Smith SA, Jones CK, Zackowski KM, van Zijl PC, Calabresi PA, et al. Quantitative characterization of the corticospinal tract at 3T. AJNR Am J Neuroradiol. 2006;27: 2168–2178. pmid:17110689
  66. 66. van der Knaap MS, Valk J. Magnetic resonance in myelin, myelination, and myelin disorders. Berlin Heidelberg: Springer; 2005.
  67. 67. Andersson JLR, Skare S, Ashburner J. How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging. Neuroimage. 2003;20: 870–888. pmid:14568458
  68. 68. Raffelt, D., Tournier, J.-D., Salvado, O., Connelly, A. Mask-Based Motion and Eddy-Current Correction of High b-value Diffusion-Weighted Images. Proceedings of the ISMRM. 2012; 3555.
  69. 69. Holland D, Kuperman JM, Dale AM. Efficient correction of inhomogeneous static magnetic field-induced distortion in Echo Planar Imaging. Neuroimage. 2010;50: 175–183. pmid:19944768
  70. 70. Smith RE, Tournier J-D, Calamante F, Connelly A. Anatomically-constrained tractography: improved diffusion MRI streamlines tractography through effective use of anatomical information. Neuroimage. 2012;62: 1924–1938. pmid:22705374
  71. 71. Vos SB, Tax CMW, Luijten PR, Ourselin S, Leemans A, Froeling M. The importance of correcting for signal drift in diffusion MRI. Magn Reson Med. 2016;
  72. 72. Basser PJ, Pierpaoli C. Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. J Magn Reson B. 1996;111: 209–219. pmid:8661285
  73. 73. Tabesh A, Jensen JH, Ardekani BA, Helpern JA. Estimation of tensors and tensor-derived measures in diffusional kurtosis imaging. Magnetic Resonance in Medicine. 2011;65: 823–836. pmid:21337412
  74. 74. Provencher SW. Estimation of metabolite concentrations from localized in vivo proton NMR spectra. Magn Reson Med. 1993;30: 672–679. pmid:8139448
  75. 75. Liu T, Liu J, de Rochefort L, Spincemaille P, Khalidov I, Ledoux JR, et al. Morphology enabled dipole inversion (MEDI) from a single-angle acquisition: Comparison with COSMOS in human brain imaging. Magnetic Resonance in Medicine. 2011;66: 777–783. pmid:21465541
  76. 76. Smith SM. Fast robust automated brain extraction. Hum Brain Mapp. 2002;17: 143–155. pmid:12391568
  77. 77. Sun H, Wilman AH. Background field removal using spherical mean value filtering and Tikhonov regularization. Magn Reson Med. 2014;71: 1151–1157. pmid:23666788
  78. 78. Bilgic B, Fan AP, Polimeni JR, Cauley SF, Bianciardi M, Adalsteinsson E, et al. Fast quantitative susceptibility mapping with L1-regularization and automatic parameter selection: Fast QSM with L1-Regularization. Magnetic Resonance in Medicine. 2014;72: 1444–1459. pmid:24259479
  79. 79. Tournier J-D, Calamante F, Connelly A. MRtrix: Diffusion tractography in crossing fiber regions. International Journal of Imaging Systems and Technology. 2012;22: 53–66.
  80. 80. Jenkinson M, Beckmann CF, Behrens TEJ, Woolrich MW, Smith SM. FSL. Neuroimage. 2012;62: 782–790. pmid:21979382
  81. 81. Tournier J-D, Calamante F, Connelly A. Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage. 2007;35: 1459–1472. pmid:17379540
  82. 82. Hirayama K, Tsubaki T, Toyokura Y, Okinaka S. The representation of the pyramidal tract in the internal capsule and basis pedunculi. A study based on three cases of amyotrophic lateral sclerosis. Neurology. 1962;12: 337–342. pmid:13907643
  83. 83. Minati L, Aquino D, Rampoldi S, Papa S, Grisoli M, Bruzzone MG, et al. Biexponential and diffusional kurtosis imaging, and generalised diffusion-tensor imaging (GDTI) with rank-4 tensors: a study in a group of healthy subjects. Magnetic Resonance Materials in Physics, Biology and Medicine. 2007;20: 241–253.
  84. 84. Falangola MF, Jensen JH, Babb JS, Hu C, Castellanos FX, Di Martino A, et al. Age-related non-Gaussian diffusion patterns in the prefrontal brain. Journal of Magnetic Resonance Imaging. 2008;28: 1345–1350. pmid:19025941
  85. 85. Grinberg F, Maximov II, Farrher E, Neuner I, Amort L, Thönneßen H, et al. Diffusion kurtosis metrics as biomarkers of microstructural development: a comparative study of a group of children and a group of adults. NeuroImage. 2016;
  86. 86. Jones DK, Knösche TR, Turner R. White matter integrity, fiber count, and other fallacies: the do’s and don’ts of diffusion MRI. Neuroimage. 2013;73: 239–254. pmid:22846632
  87. 87. Chang YS, Owen JP, Pojman NJ, Thieu T, Bukshpun P, Wakahiro MLJ, et al. White Matter Changes of Neurite Density and Fiber Orientation Dispersion during Human Brain Maturation. PLoS ONE. 2015;10: e0123656. pmid:26115451
  88. 88. Kunz N, Zhang H, Vasung L, O’Brien KR, Assaf Y, Lazeyras F, et al. Assessing white matter microstructure of the newborn with multi-shell diffusion MRI and biophysical compartment models. Neuroimage. 2014;
  89. 89. Raffelt D, Tournier J-D, Rose S, Ridgway GR, Henderson R, Crozier S, et al. Apparent Fibre Density: a novel measure for the analysis of diffusion-weighted magnetic resonance images. Neuroimage. 2012;59: 3976–3994. pmid:22036682
  90. 90. Dell’Acqua F, Simmons A, Williams SCR, Catani M. Can spherical deconvolution provide more information than fiber orientations? Hindrance modulated orientational anisotropy, a true-tract specific index to characterize white matter diffusion. Hum Brain Mapp. 2013;34: 2464–2483. pmid:22488973
  91. 91. Alexander DC, Hubbard PL, Hall MG, Moore EA, Ptito M, Parker GJM, et al. Orientationally invariant indices of axon diameter and density from diffusion MRI. Neuroimage. 2010;52: 1374–1389. pmid:20580932
  92. 92. Assaf Y, Blumenfeld-Katzir T, Yovel Y, Basser PJ. AxCaliber: a method for measuring axon diameter distribution from diffusion MRI. Magn Reson Med. 2008;59: 1347–1354. pmid:18506799
  93. 93. Jespersen SN, Bjarkam CR, Nyengaard JR, Chakravarty MM, Hansen B, Vosegaard T, et al. Neurite density from magnetic resonance diffusion measurements at ultrahigh field: comparison with light microscopy and electron microscopy. Neuroimage. 2010;49: 205–216. pmid:19732836
  94. 94. Sukstanskii AL, Yablonskiy DA, Ackerman JJH. Effects of permeable boundaries on the diffusion-attenuated MR signal: insights from a one-dimensional model. J Magn Reson. 2004;170: 56–66. pmid:15324758
  95. 95. Wang Y, Wang Q, Haldar JP, Yeh F-C, Xie M, Sun P, et al. Quantification of increased cellularity during inflammatory demyelination. Brain. 2011;134: 3590–3601. pmid:22171354
  96. 96. Panagiotaki E, Walker-Samuel S, Siow B, Johnson SP, Rajkumar V, Pedley RB, et al. Noninvasive quantification of solid tumor microstructure using VERDICT MRI. Cancer Res. 2014;74: 1902–1912. pmid:24491802
  97. 97. Kaden E, Kruggel F, Alexander DC. Quantitative mapping of the per-axon diffusion coefficients in brain white matter. Magn Reson Med. 2016;75: 1752–1763. pmid:25974332
  98. 98. Burcaw LM, Fieremans E, Novikov DS. Mesoscopic structure of neuronal tracts from time-dependent diffusion. NeuroImage. 2015;114: 18–37. pmid:25837598
  99. 99. Fieremans E, Burcaw LM, Lee H-H, Lemberskiy G, Veraart J, Novikov DS. In vivo observation and biophysical interpretation of time-dependent diffusion in human white matter. Neuroimage. 2016;129: 414–427. pmid:26804782
  100. 100. Jelescu IO, Veraart J, Fieremans E, Novikov DS. Degeneracy in model parameter estimation for multi-compartmental diffusion in neuronal tissue: Degeneracy in Model Parameter Estimation of Diffusion in Neural Tissue. NMR in Biomedicine. 2016;29: 33–47. pmid:26615981
  101. 101. Maudsley AA, Domenig C, Govind V, Darkazanli A, Studholme C, Arheart K, et al. Mapping of brain metabolite distributions by volumetric proton MR spectroscopic imaging (MRSI). Magnetic Resonance in Medicine. 2009;61: 548–559. pmid:19111009
  102. 102. D’Adamo AF, Yatsu FM. Acetate metabolism in the nervous system. N-acetyl-L-aspartic acid and the biosynthesis of brain lipids. J Neurochem. 1966;13: 961–965. pmid:5927765
  103. 103. Deoni SC, Rutt BK, Arun T, Pierpaoli C, Jones DK. Gleaning multicomponent T1 and T2 information from steady-state imaging data. Magn Reson Med. 2008;60: 1372–1387. pmid:19025904
  104. 104. Zhang J, Kolind SH, Laule C, MacKay AL. Comparison of myelin water fraction from multiecho T2 decay curve and steady-state methods. Magn Reson Med. 2015;73: 223–232. pmid:24515972
  105. 105. Cercignani M, Symms MR, Ron M, Barker GJ. 3D MTR measurement: from 1.5 T to 3.0 T. Neuroimage. 2006;31: 181–186. pmid:16413794
  106. 106. Helms G, Dathe H, Dechent P. Modeling the influence of TR and excitation flip angle on the magnetization transfer ratio (MTR) in human brain obtained from 3D spoiled gradient echo MRI. Magn Reson Med. 2010;64: 177–185. pmid:20572139
  107. 107. Tofts PS, Steens SCA, Cercignani M, Admiraal-Behloul F, Hofman P a. M, van Osch MJP, et al. Sources of variation in multi-centre brain MTR histogram studies: body-coil transmission eliminates inter-centre differences. MAGMA. 2006;19: 209–222. pmid:16957936
  108. 108. Liu C, Li W, Tong KA, Yeom KW, Kuzminski S. Susceptibility-weighted imaging and quantitative susceptibility mapping in the brain. J Magn Reson Imaging. 2014;
  109. 109. Li W, Wu B, Liu C. Quantitative susceptibility mapping of human brain reflects spatial variation in tissue composition. Neuroimage. 2011;55: 1645–1656. pmid:21224002
  110. 110. Wharton S, Bowtell R. Effects of white matter microstructure on phase and susceptibility maps. Magn Reson Med. 2015;73: 1258–1269. pmid:24619643
  111. 111. Wharton S, Bowtell R. Gradient echo based fiber orientation mapping using R2* and frequency difference measurements. Neuroimage. 2013;83: 1011–1023. pmid:23906549
  112. 112. Yablonskiy DA, Sukstanskii AL. Generalized Lorentzian Tensor Approach (GLTA) as a biophysical background for quantitative susceptibility mapping. Magn Reson Med. 2015;73: 757–764. pmid:25426775
  113. 113. Liu C, Li W, Wu B, Jiang Y, Johnson GA. 3D fiber tractography with susceptibility tensor imaging. Neuroimage. 2012;59: 1290–1298. pmid:21867759
  114. 114. Schweser F, Deistung A, Lehr BW, Reichenbach JR. Differentiation between diamagnetic and paramagnetic cerebral lesions based on magnetic susceptibility mapping. Med Phys. 2010;37: 5165–5178. pmid:21089750
  115. 115. Langkammer C, Bredies K, Poser BA, Barth M, Reishofer G, Fan AP, et al. Fast quantitative susceptibility mapping using 3D EPI and total generalized variation. Neuroimage. 2015;111: 622–630. pmid:25731991
  116. 116. de Rochefort L, Brown R, Prince MR, Wang Y. Quantitative MR susceptibility mapping using piece-wise constant regularized inversion of the magnetic field. Magn Reson Med. 2008;60: 1003–1009. pmid:18816834
  117. 117. Bilgic B, Pfefferbaum A, Rohlfing T, Sullivan EV, Adalsteinsson E. MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping. Neuroimage. 2012;59: 2625–2635. pmid:21925274
  118. 118. Wharton S, Schäfer A, Bowtell R. Susceptibility mapping in the human brain using threshold-based k-space division. Magn Reson Med. 2010;63: 1292–1304. pmid:20432300
  119. 119. Bender B, Klose U. The in vivo influence of white matter fiber orientation towards B(0) on T2* in the human brain. NMR Biomed. 2010;23: 1071–1076. pmid:20665897
  120. 120. Cherubini A, Péran P, Hagberg GE, Varsi AE, Luccichenti G, Caltagirone C, et al. Characterization of white matter fiber bundles with T2* relaxometry and diffusion tensor imaging. Magn Reson Med. 2009;61: 1066–1072. pmid:19253372
  121. 121. Sati P, van Gelderen P, Silva AC, Reich DS, Merkle H, de Zwart JA, et al. Micro-compartment specific T2* relaxation in the brain. NeuroImage. 2013;77: 268–278. pmid:23528924
  122. 122. Mohammadi S, Carey D, Dick F, Diedrichsen J, Sereno MI, Reisert M, et al. Whole-Brain In-vivo Measurements of the Axonal G-Ratio in a Group of 37 Healthy Volunteers. Front Neurosci. 2015;9: 441. pmid:26640427
  123. 123. Aboitiz F, Scheibel AB, Fisher RS, Zaidel E. Fiber composition of the human corpus callosum. Brain Res. 1992;598: 143–153. pmid:1486477
  124. 124. MacKay A, Laule C, Vavasour I, Bjarnason T, Kolind S, Mädler B. Insights into brain microstructure from the T2 distribution. Magn Reson Imaging. 2006;24: 515–525. pmid:16677958
  125. 125. Schmierer K, Scaravilli F, Altmann DR, Barker GJ, Miller DH. Magnetization transfer ratio and myelin in postmortem multiple sclerosis brain. Ann Neurol. 2004;56: 407–415. pmid:15349868
  126. 126. Rushton W a. H. A theory of the effects of fibre size in medullated nerve. J Physiol (Lond). 1951;115: 101–122.
  127. 127. Gelman N, Ewing JR, Gorell JM, Spickler EM, Solomon EG. Interregional variation of longitudinal relaxation rates in human brain at 3.0 T: relation to estimated iron and water contents. Magn Reson Med. 2001;45: 71–79. pmid:11146488
  128. 128. Bock NA, Kocharyan A, Liu JV, Silva AC. Visualizing the entire cortical myelination pattern in marmosets with magnetic resonance imaging. J Neurosci Methods. 2009;185: 15–22. pmid:19737577
  129. 129. Stüber C, Morawski M, Schäfer A, Labadie C, Wähnert M, Leuze C, et al. Myelin and iron concentration in the human brain: a quantitative study of MRI contrast. Neuroimage. 2014;93 Pt 1: 95–106.
  130. 130. Schneider J, Kober T, Bickle Graz M, Meuli R, Hüppi PS, Hagmann P, et al. Evolution of T1 Relaxation, ADC, and Fractional Anisotropy during Early Brain Maturation: A Serial Imaging Study on Preterm Infants. AJNR Am J Neuroradiol. 2015;
  131. 131. Labadie C, Lee J-H, Rooney WD, Jarchow S, Aubert-Frécon M, Springer CS, et al. Myelin water mapping by spatially regularized longitudinal relaxographic imaging at high magnetic fields. Magn Reson Med. 2014;71: 375–387. pmid:23468414
  132. 132. Rioux JA, Levesque IR, Rutt BK. Biexponential longitudinal relaxation in white matter: Characterization and impact on T1 mapping with IR-FSE and MP2RAGE. Magn Reson Med. 2015;
  133. 133. McNab JA, Edlow BL, Witzel T, Huang SY, Bhat H, Heberlein K, et al. The Human Connectome Project and beyond: initial applications of 300 mT/m gradients. Neuroimage. 2013;80: 234–245. pmid:23711537
  134. 134. Barth M, Breuer F, Koopmans PJ, Norris DG, Poser BA. Simultaneous multislice (SMS) imaging techniques. Magn Reson Med. 2015;
  135. 135. Feinberg DA, Moeller S, Smith SM, Auerbach E, Ramanna S, Gunther M, et al. Multiplexed echo planar imaging for sub-second whole brain FMRI and fast diffusion imaging. PLoS ONE. 2010;5: e15710. pmid:21187930