Temporal Diffusion Ratio (TDR) for imaging restricted diffusion: Optimisation and pre-clinical demonstration.

Temporal Diffusion Ratio (TDR) is a recently proposed dMRI technique (Dell'Acqua et al., proc. ISMRM 2019) which provides contrast between areas with restricted diffusion and areas either without restricted diffusion or with length scales too small for characterisation. Hence, it has a potential for informing on pore sizes, in particular the presence of large axon diameters or other cellular structures. TDR employs the signal from two dMRI acquisitions obtained with the same, large, b-value but with different diffusion gradient waveforms. TDR is advantageous as it employs standard acquisition sequences, does not make any assumptions on the underlying tissue structure and does not require any model fitting, avoiding issues related to model degeneracy. This work for the first time introduces and optimises the TDR method in simulation for a range of different tissues and scanner constraints and validates it in a pre-clinical demonstration. We consider both substrates containing cylinders and spherical structures, representing cell soma in tissue. Our results show that contrasting an acquisition with short gradient duration, short diffusion time and high gradient strength with an acquisition with long gradient duration, long diffusion time and low gradient strength, maximises the TDR contrast for a wide range of pore configurations. Additionally, in the presence of Rician noise, computing TDR from a subset (50% or fewer) of the acquired diffusion gradients rather than the entire shell as proposed originally further improves the contrast. In the last part of the work the results are demonstrated experimentally on rat spinal cord. In line with simulations, the experimental data shows that optimised TDR improves the contrast compared to non-optimised TDR. Furthermore, we find a strong correlation between TDR and histology measurements of axon diameter. In conclusion, we find that TDR has great potential and is a very promising alternative (or potentially complement) to model-based approaches for informing on pore sizes and restricted diffusion in general.

Tensor denoising of multidimensional MRI data.

PURPOSE: To develop a denoising strategy leveraging redundancy in high-dimensional data. THEORY AND METHODS: The SNR fundamentally limits the information accessible by MRI. This limitation has been addressed by a host of denoising techniques, recently including the so-called MPPCA: principal component analysis of the signal followed by automated rank estimation, exploiting the Marchenko-Pastur distribution of noise singular values. Operating on matrices comprised of data patches, this popular approach objectively identifies noise components and, ideally, allows noise to be removed without introducing artifacts such as image blurring, or nonlocal averaging. The MPPCA rank estimation, however, relies on a large number of noise singular values relative to the number of signal components to avoid such ill effects. This condition is unlikely to be met when data patches and therefore matrices are small, for example due to spatially varying noise. Here, we introduce tensor MPPCA (tMPPCA) for the purpose of denoising multidimensional data, such as from multicontrast acquisitions. Rather than combining dimensions in matrices, tMPPCA uses each dimension of the multidimensional data's inherent tensor-structure to better characterize noise, and to recursively estimate signal components. RESULTS: Relative to matrix-based MPPCA, tMPPCA requires no additional assumptions, and comparing the two in a numerical phantom and a multi-TE diffusion MRI data set, tMPPCA dramatically improves denoising performance. This is particularly true for small data patches, suggesting that tMPPCA can be especially beneficial in such cases. CONCLUSIONS: The MPPCA denoising technique can be extended to high-dimensional data with improved performance for smaller patch sizes.

Advanced Diffusion-Weighted MRI for Cancer Microstructure Assessment in Body Imaging, and Its Relationship With Histology.

Diffusion-weighted magnetic resonance imaging (DW-MRI) aims to disentangle multiple biological signal sources in each imaging voxel, enabling the computation of innovative maps of tissue microstructure. DW-MRI model development has been dominated by brain applications. More recently, advanced methods with high fidelity to histology are gaining momentum in other contexts, for example, in oncological applications of body imaging, where new biomarkers are urgently needed. The objective of this article is to review the state-of-the-art of DW-MRI in body imaging (ie, not including the nervous system) in oncology, and to analyze its value as compared to reference colocalized histology measurements, given that demonstrating the histological validity of any new DW-MRI method is essential. In this article, we review the current landscape of DW-MRI techniques that extend standard apparent diffusion coefficient (ADC), describing their acquisition protocols, signal models, fitting settings, microstructural parameters, and relationship with histology. Preclinical, clinical, and in/ex vivo studies were included. The most used techniques were intravoxel incoherent motion (IVIM; 36.3% of used techniques), diffusion kurtosis imaging (DKI; 16.7%), vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT; 13.3%), and imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED; 11.7%). Another notable category of techniques relates to innovative b-tensor diffusion encoding or joint diffusion-relaxometry. The reviewed approaches provide histologically meaningful indices of cancer microstructure (eg, vascularization/cellularity) which, while not necessarily accurate numerically, may still provide useful sensitivity to microscopic pathological processes. Future work of the community should focus on improving the inter-/intra-scanner robustness, and on assessing histological validity in broader contexts. LEVEL OF EVIDENCE: NA TECHNICAL EFFICACY: Stage 2.

In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner.

Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-Gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (Kiso), variance due to diffusion anisotropy (Kaniso), and microscopic kurtosis (µK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, µK is typically ignored in diffusion MRI signal modelling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible µK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of µK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that µK significantly contributes to total diffusional kurtosis both in grey and white matter tissue but, as expected, not in the ventricles. The first µK maps of the human brain are presented, revealing that the spatial distribution of µK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring µK and assuming the multiple Gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.

Soma and Neurite Density MRI (SANDI) of the in-vivo mouse brain and comparison with the Allen Brain Atlas.

Diffusion MRI (dMRI) provides unique insights into the neural tissue milieu by probing interactions between diffusing molecules and tissue microstructure. Most dMRI techniques focus on white matter (WM) tissues, nevertheless, interest in gray matter characterizations is growing. The Soma and Neurite Density MRI (SANDI) methodology harnesses a model incorporating water diffusion in spherical objects (assumed to be associated with cell bodies) and in impermeable "sticks" (assumed to represent neurites), which potentially enables the characterization of cellular and neurite densities. Recognising the importance of rodents in animal models of development, aging, plasticity, and disease, we here employ SANDI for in-vivo preclinical imaging and provide a first validation of the methodology by comparing SANDI metrics with cellular density reflected by the Allen mouse brain atlas. SANDI was implemented on a 9.4T scanner equipped with a cryogenic coil, and in-vivo experiments were carried out on N = 6 mice. Pixelwise, ROI-based, and atlas comparisons were performed, magnitude vs. real-valued analyses were compared, and shorter acquisitions with reduced the number of b-value shells were investigated. Our findings reveal good reproducibility of the SANDI parameters, including the sphere and stick fractions, as well as sphere size (CoV < 7%, 12% and 3%, respectively). Additionally, we find a very good rank correlation between SANDI-driven sphere fraction and Allen mouse brain atlas contrast that represents cellular density. We conclude that SANDI is a viable preclinical MRI technique that can greatly contribute to research on brain tissue microstructure.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

SANDI: A compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI.

This work introduces a compartment-based model for apparent cell body (namely soma) and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging through DW-MRI presents water diffusion in white (WM) and gray (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b â€‹≤ â€‹3,000 â€‹s/mm2 (or 3 â€‹ms/µm2), it has been also shown to fail in GM at high b values (b≫3,000 â€‹s/mm2 or 3 â€‹ms/µm2). Here we hypothesise that the unmodelled soma compartment (i.e. cell body of any brain cell type: from neuroglia to neurons) may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma of any brain cell type is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (bmax â€‹= â€‹40,000 â€‹s/mm2, or 40 â€‹ms/µm2) and human (bmax â€‹= â€‹10,000 â€‹s/mm2, or 10 â€‹ms/µm2) brain. We show that SANDI defines new contrasts representing complementary information on the brain cyto- and myelo-architecture. Indeed, we show maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto- and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.

Higher-order diffusion MRI characterization of mesorectal lymph nodes in rectal cancer.

PURPOSE: Mesorectal lymph node staging plays an important role in treatment decision making. Here, we explore the benefit of higher-order diffusion MRI models accounting for non-Gaussian diffusion effects to classify mesorectal lymph nodes both 1) ex vivo at ultrahigh field correlated with histology and 2) in vivo in a clinical scanner upon patient staging. METHODS: The preclinical investigation included 54 mesorectal lymph nodes, which were scanned at 16.4 T with an extensive diffusion MRI acquisition. Eight diffusion models were compared in terms of goodness of fit, lymph node classification ability, and histology correlation. In the clinical part of this study, 10 rectal cancer patients were scanned with diffusion MRI at 1.5 T, and 72 lymph nodes were analyzed with Apparent Diffusion Coefficient (ADC), Intravoxel Incoherent Motion (IVIM), Kurtosis, and IVIM-Kurtosis. RESULTS: Compartment models including restricted and anisotropic diffusion improved the preclinical data fit, as well as the lymph node classification, compared to standard ADC. The comparison with histology revealed only moderate correlations, and the highest values were observed between diffusion anisotropy metrics and cell area fraction. In the clinical study, the diffusivity from IVIM-Kurtosis was the only metric showing significant differences between benign (0.80 ± 0.30 µm2 /ms) and malignant (1.02 ± 0.41 µm2 /ms, P = .03) nodes. IVIM-Kurtosis also yielded the largest area under the receiver operating characteristic curve (0.73) and significantly improved the node differentiation when added to the standard visual analysis by experts based on T2 -weighted imaging. CONCLUSION: Higher-order diffusion MRI models perform better than standard ADC and may be of added value for mesorectal lymph node classification in rectal cancer patients.

Measuring diffusion exchange across the cell membrane with DEXSY (Diffusion Exchange Spectroscopy).

INTRODUCTION: To combine numerical simulations, in vitro and in vivo experiments to evaluate the feasibility of measuring diffusion exchange across the cell membrane with diffusion exchange spectroscopy (DEXSY). METHODS: DEXSY acquisitions were simulated over a range of permeabilities in nerve tissue and yeast substrates. In vitro measurements were performed in a yeast substrate and in vivo measurements in mouse tumor xenograft models, all at 9.4 T. RESULTS: Diffusion exchange was observed in simulations over a physiologically relevant range of cell permeability values. In vitro and in vivo measures also provided evidence of diffusion exchange, which was quantified with the Diffusion Exchange Index (DEI). CONCLUSIONS: Our findings provide preliminary evidence that DEXSY can be used to make in vivo measurements of diffusion exchange and cell membrane permeability.

Layer-specific connectivity revealed by diffusion-weighted functional MRI in the rat thalamocortical pathway.

Investigating neural activity from a global brain perspective in-vivo has been in the domain of functional Magnetic Resonance Imaging (fMRI) over the past few decades. The intricate neurovascular couplings that govern fMRI's blood-oxygenation-level-dependent (BOLD) functional contrast are invaluable in mapping active brain regions, but they also entail significant limitations, such as non-specificity of the signal to active foci. Diffusion-weighted functional MRI (dfMRI) with relatively high diffusion-weighting strives to ameliorate this shortcoming as it offers functional contrasts more intimately linked with the underlying activity. Insofar, apart from somewhat smaller activation foci, dfMRI's contrasts have not been convincingly shown to offer significant advantages over BOLD-driven fMRI, and its activation maps relied on significant modelling. Here, we study whether dfMRI could offer a better representation of neural activity in the thalamocortical pathway compared to its (spin-echo (SE)) BOLD counterpart. Using high-end forepaw stimulation experiments in the rat at 9.4 T, and with significant sensitivity enhancements due to the use of cryocoils, we show for the first time that dfMRI signals exhibit layer specificity, and, additionally, display signals in areas devoid of SE-BOLD responses. We find that dfMRI signals in the thalamocortical pathway cohere with each other, namely, dfMRI signals in the ventral posterolateral (VPL) thalamic nucleus cohere specifically with layers IV and V in the somatosensory cortex. These activity patterns are much better correlated (compared with SE-BOLD signals) with literature-based electrophysiological recordings in the cortex as well as thalamus. All these findings suggest that dfMRI signals better represent the underlying neural activity in the pathway. In turn, these advanatages may have significant implications towards a much more specific and accurate mapping of neural activity in the global brain in-vivo.

Relevance of time-dependence for clinically viable diffusion imaging of the spinal cord.

PURPOSE: Time-dependence is a key feature of the diffusion-weighted (DW) signal, knowledge of which informs biophysical modelling. Here, we study time-dependence in the human spinal cord, as its axonal structure is specific and different from the brain. METHODS: We run Monte Carlo simulations using a synthetic model of spinal cord white matter (WM) (large axons), and of brain WM (smaller axons). Furthermore, we study clinically feasible multi-shell DW scans of the cervical spinal cord (b = 0; b = 711 s mm-2 ; b = 2855 s mm-2 ), obtained using three diffusion times (Δ of 29, 52 and 76 ms) from three volunteers. RESULTS: Both intra-/extra-axonal perpendicular diffusivities and kurtosis excess show time-dependence in our synthetic spinal cord model. This time-dependence is reflected mostly in the intra-axonal perpendicular DW signal, which also exhibits strong decay, unlike our brain model. Time-dependence of the total DW signal appears detectable in the presence of noise in our synthetic spinal cord model, but not in the brain. In WM in vivo, we observe time-dependent macroscopic and microscopic diffusivities and diffusion kurtosis, NODDI and two-compartment SMT metrics. Accounting for large axon calibers improves fitting of multi-compartment models to a minor extent. CONCLUSIONS: Time-dependence of clinically viable DW MRI metrics can be detected in vivo in spinal cord WM, thus providing new opportunities for the non-invasive estimation of microstructural properties. The time-dependence of the perpendicular DW signal may feature strong intra-axonal contributions due to large spinal axon caliber. Hence, a popular model known as "stick" (zero-radius cylinder) may be sub-optimal to describe signals from the largest spinal axons.

Accurate estimation of microscopic diffusion anisotropy and its time dependence in the mouse brain.

Microscopic diffusion anisotropy (µA) has been recently gaining increasing attention for its ability to decouple the average compartment anisotropy from orientation dispersion. Advanced diffusion MRI sequences, such as double diffusion encoding (DDE) and double oscillating diffusion encoding (DODE) have been used for mapping µA, usually using measurements from a single b shell. However, the accuracy of µA estimation vis-à-vis different b-values was not assessed. Moreover, the time-dependence of this metric, which could offer additional insights into tissue microstructure, has not been studied so far. Here, we investigate both these concepts using theory, simulation, and experiments performed at 16.4T in the mouse brain, ex-vivo. In the first part, simulations and experimental results show that the conventional estimation of microscopic anisotropy from the difference of D(O)DE sequences with parallel and orthogonal gradient directions yields values that highly depend on the choice of b-value. To mitigate this undesirable bias, we propose a multi-shell approach that harnesses a polynomial fit of the signal difference up to third order terms in b-value. In simulations, this approach yields more accurate µA metrics, which are similar to the ground-truth values. The second part of this work uses the proposed multi-shell method to estimate the time/frequency dependence of µA. The data shows either an increase or no change in µA with frequency depending on the region of interest, both in white and gray matter. When comparing the experimental results with simulations, it emerges that simple geometric models such as infinite cylinders with either negligible or finite radii cannot replicate the measured trend, and more complex models, which, for example, incorporate structure along the fibre direction are required. Thus, measuring the time dependence of microscopic anisotropy can provide valuable information for characterizing tissue microstructure.

Low frequency oscillating gradient spin-echo sequences improve sensitivity to axon diameter: An experimental study in viable nerve tissue.

Mapping axon diameters within the central and peripheral nervous system could play an important role in our understanding of nerve pathways, and help diagnose and monitor an array of neurological disorders. Numerous diffusion MRI methods have been proposed for imaging axon diameters, most of which use conventional single diffusion encoding (SDE) spin echo sequences. However, a growing number of studies show that oscillating gradient spin echo (OGSE) sequences can provide additional advantages over conventional SDE sequences. Recent theoretical results suggest that this is especially the case in realistic scenarios, such as when fibres have unknown or dispersed orientation. In the present study, we adopt the ActiveAx approach to experimentally investigate the extent of these advantages by comparing the performances of SDE and trapezoidal OGSE in viable nerve tissue. We optimise SDE and OGSE ActiveAx protocols for a rat peripheral nerve tissue and test their performance using Monte Carlo simulations and a 800 mT/m gradient strength pre-clinical imaging experiment. The imaging experiment uses excised sciatic nerve from a rat's leg placed in a MRI compatible viable isolated tissue (VIT) maintenance chamber, which keeps the tissue in a viable physiological state that preserves the structural complexity of the nerve and enables lengthy scan times. We compare model estimates to histology, which we perform on the nerve post scanning. Optimisation produces a three-shell SDE and OGSE ActiveAx protocol, with the OGSE protocol consisting of one SDE sequence and two low-frequency oscillating gradient waveform sequences. Both simulation and imaging results show that the OGSE ActiveAx estimates of the axon diameter index have a higher accuracy and a higher precision compared to those from SDE. Histology estimates of the axon diameter index in our nerve tissue samples are 4-5.8 µm and these are excellently matched with the OGSE estimates 4.2-6.5 µm, while SDE overestimates at 5.2-8 µm for the same sample. We found OGSE estimates to be more precise with on average a 0.5 µm standard deviation compared to the SDE estimates which have a 2 µm standard deviation. When testing the robustness of the estimates when the number of the diffusion gradient directions reduces, we found that both OGSE and SDE estimates are affected, however OGSE is more robust to these changes than the SDE. Overall, these results suggest, quantitatively and in in vivo conditions, that low-frequency OGSE sequences may provide improved accuracy of axon diameter mapping compared to standard SDE sequences.

Incomplete initial nutation diffusion imaging: An ultrafast, single-scan approach for diffusion mapping.

PURPOSE: Diffusion MRI is confounded by the need to acquire at least two images separated by a repetition time, thereby thwarting the detection of rapid dynamic microstructural changes. The issue is exacerbated when diffusivity variations are accompanied by rapid changes in T2 . The purpose of the present study is to accelerate diffusion MRI acquisitions such that both reference and diffusion-weighted images necessary for quantitative diffusivity mapping are acquired in a single-shot experiment. METHODS: A general methodology termed incomplete initial nutation diffusion imaging (INDI), capturing two diffusion contrasts in a single shot, is presented. This methodology creates a longitudinal magnetization reservoir that facilitates the successive acquisition of two images separated by only a few milliseconds. The theory behind INDI is presented, followed by proof-of-concept studies in water phantom, ex vivo, and in vivo experiments at 16.4 and 9.4 T. RESULTS: Mean diffusivities extracted from INDI were comparable with diffusion tensor imaging and the two-shot isotropic diffusion encoding in the water phantom. In ex vivo mouse brain tissues, as well as in the in vivo mouse brain, mean diffusivities extracted from conventional isotropic diffusion encoding and INDI were in excellent agreement. Simulations for signal-to-noise considerations identified the regimes in which INDI is most beneficial. CONCLUSIONS: The INDI method accelerates diffusion MRI acquisition to single-shot mode, which can be of great importance for mapping dynamic microstructural properties in vivo without T2 bias. Magn Reson Med 79:2198-2204, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

An optimized framework for quantitative magnetization transfer imaging of the cervical spinal cord in vivo.

PURPOSE: To develop a framework to fully characterize quantitative magnetization transfer indices in the human cervical cord in vivo within a clinically feasible time. METHODS: A dedicated spinal cord imaging protocol for quantitative magnetization transfer was developed using a reduced field-of-view approach with echo planar imaging (EPI) readout. Sequence parameters were optimized based in the Cramer-Rao-lower bound. Quantitative model parameters (i.e., bound pool fraction, free and bound pool transverse relaxation times [ T2F, T2B], and forward exchange rate [kFB ]) were estimated implementing a numerical model capable of dealing with the novelties of the sequence adopted. The framework was tested on five healthy subjects. RESULTS: Cramer-Rao-lower bound minimization produces optimal sampling schemes without requiring the establishment of a steady-state MT effect. The proposed framework allows quantitative voxel-wise estimation of model parameters at the resolution typically used for spinal cord imaging (i.e. 0.75 × 0.75 × 5 mm3 ), with a protocol duration of ∼35 min. Quantitative magnetization transfer parametric maps agree with literature values. Whole-cord mean values are: bound pool fraction = 0.11(±0.01), T2F = 46.5(±1.6) ms, T2B = 11.0(±0.2) µs, and kFB = 1.95(±0.06) Hz. Protocol optimization has a beneficial effect on reproducibility, especially for T2B and kFB . CONCLUSION: The framework developed enables robust characterization of spinal cord microstructure in vivo using qMT. Magn Reson Med 79:2576-2588, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Double oscillating diffusion encoding and sensitivity to microscopic anisotropy.

PURPOSE: To introduce a novel diffusion pulse sequence, namely double oscillating diffusion encoding (DODE), and to investigate whether it adds sensitivity to microscopic diffusion anisotropy (µA) compared to the well-established double diffusion encoding (DDE) methodology. METHODS: We simulate measurements from DODE and DDE sequences for different types of microstructures exhibiting restricted diffusion. First, we compare the effect of varying pulse sequence parameters on the DODE and DDE signal. Then, we analyse the sensitivity of the two sequences to the microstructural parameters (pore diameter and length) which determine µA. Finally, we investigate specificity of measurements to particular substrate configurations. RESULTS: Simulations show that DODE sequences exhibit similar signal dependence on the relative angle between the two gradients as DDE sequences, however, the effect of varying the mixing time is less pronounced. The sensitivity analysis shows that in substrates with elongated pores and various orientations, DODE sequences increase the sensitivity to pore diameter, while DDE sequences are more sensitive to pore length. Moreover, DDE and DODE sequence parameters can be tailored to enhance/suppress the signal from a particular range of substrates. CONCLUSIONS: A combination of DODE and DDE sequences maximize sensitivity to µA, compared to using just the DDE method. Magn Reson Med 78:550-564, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

PGSE, OGSE, and sensitivity to axon diameter in diffusion MRI: Insight from a simulation study.

PURPOSE: To identify optimal pulsed gradient spin-echo (PGSE) and oscillating gradient spin-echo (OGSE) sequence settings for maximizing sensitivity to axon diameter in idealized and practical conditions. METHODS: Simulations on a simple two-compartment white matter model (with nonpermeable cylinders) are used to investigate a wide space of clinically plausible PGSE and OGSE sequence parameters with trapezoidal diffusion gradient waveforms. Signal sensitivity is measured as a derivative of the signal with respect to axon diameter. Models of parallel and dispersed fibers are investigated separately to represent idealized and practical conditions. RESULTS: Simulations show that, for the simple case of gradients perfectly perpendicular to straight parallel fibers, PGSE always gives maximum sensitivity. However, in real-world scenarios where fibers have unknown and dispersed orientation, low-frequency OGSE provides higher sensitivity. Maximum sensitivity results show that on current clinical scanners (Gmax = 60 mT/m, signal to noise ratio (SNR) = 20) axon diameters below 6 µm are indistinguishable from zero. Scanners with stronger gradient systems such as the Massachusetts General Hospital (MGH) Connectom scanner (Gmax = 300 mT/m) can extend this sensitivity limit down to 2-3 µm, probing a much greater proportion of the underlying axon diameter distribution. CONCLUSION: Low-frequency OGSE provides additional sensitivity to PGSE in practical situations. OGSE is particularly advantageous for systems with high performance gradients.