Exact analytical results for ADC with oscillating diffusion sensitizing gradients

Frequency-dependent diffusion kurtosis imaging in the human brain using an oscillating gradient spin echo sequence and a high-performance head-only gradient

Measuring the time/frequency dependence of diffusion MRI is a promising approach to distinguish between the effects of different tissue microenvironments, such as membrane restriction, tissue heterogeneity, and compartmental water exchange. In this study, we measure the frequency dependence of diffusivity (D) and kurtosis (K) with oscillating gradient diffusion encoding waveforms and a diffusion kurtosis imaging (DKI) model in human brains using a high-performance, head-only MAGNUS gradient system, with a combination of b-values, oscillating frequencies (f), and echo time that has not been achieved in human studies before. Frequency dependence of diffusivity and kurtosis are observed in both global and local white matter (WM) and gray matter (GM) regions and characterized with a power-law model ∼Λ*fθ. The frequency dependences of diffusivity and kurtosis (including changes between fmin and fmax, Λ, and θ) vary over different WM and GM regions, indicating potential microstructural differences between regions. A trend of decreasing kurtosis over frequency in the short-time limit is successfully captured for in vivo human brains. The effects of gradient nonlinearity (GNL) on frequency-dependent diffusivity and kurtosis measurements are investigated and corrected. Our results show that the GNL has prominent scaling effects on the measured diffusivity values (3.5∼5.5% difference in the global WM and 6∼8% difference in the global cortex) and subsequently affects the corresponding power-law parameters (Λ, θ) while having a marginal influence on the measured kurtosis values (<0.05% difference) and power-law parameters (Λ, θ). This study expands previous OGSE studies and further demonstrates the translatability of frequency-dependent diffusivity and kurtosis measurements to human brains, which may provide new opportunities to probe human brain microstructure in health and disease.

Rat Brain Global Ischemia-Induced Diffusion Changes Revisited: Biophysical Modeling of the Water and NAA MR "Diffusion Signal"

PURPOSE

To assess changes in intracellular diffusion as a mechanism for the reduction in water ADC that accompanies brain injury. Using NAA as a marker of neuronal cytoplasmic diffusion, NAA diffusion was measured before and after global ischemia (immediately postmortem) in the female Sprague-Dawley rat.

METHODS

Diffusion-weighted PRESS spectra, with diffusion encoding in a single direction, were acquired from large voxels of rat brain gray matter in vivo and postischemia employing either pairs of pulsed half-sine-shaped gradients (in vivo and postischemia, b_max  = 19 ms/μm² ) or sinusoidal oscillating gradients (in vivo only) with frequencies of 99.2-250 Hz. A 2D randomly oriented cylinder (neurite) model gave estimates of longitudinal and transverse diffusivities (D_L and D_T , respectively). In this model, D_L represents the "free" diffusivity of NAA, whereas D_T reflects highly restricted diffusion. Using oscillating gradients, the frequency dependence of D_T [D_T (ω)] gave estimates of the cylinder (axon/dendrite) radius.

RESULTS

A 10% decrease in D_L,NAA followed global ischemia, dropping from 0.391 ± 0.012 μm² /ms to 0.350 ± 0.009 μm² /ms. Modeling D_T,NAA (ω) provided an estimate of the neurite radius of 1.0 ± 0.6 μm.

CONCLUSION

Whereas the increase in apparent intraneuronal viscosity suggested by changes in D_L,NAA may contribute to the overall reduction in water ADC associated with brain injury, it is not sufficient to be the sole explanation. Estimates of neurite radius based on D_T (ω) were consistent with literature values.

Oscillating-gradient spin-echo diffusion-weighted imaging (OGSE-DWI) with a limited number of oscillations: I. Signal equation

The oscillating-gradient spin-echo (OGSE) sequence has been promoted as a promising diffusion-weighted magnetic resonance imaging (DWI) technique for probing in vivo tissue microstructure in the frequency domain. However, due to practical restrictions on the duration and number of oscillations that a motion-probing gradient can have, the technique has limited spectral resolution and range. This work re-examines the OGSE-DWI method to clarify how these limitations are reflected in the signal model. There are several aspects of the revised framework that distinguish it from the conventional description employed for OGSE DWI. In particular, while the conventional OGSE signal model implies that the spectral density of molecular diffusion may be directly sampled in experiments, in practice information about the spectral density can be only indirectly obtained.

Probing neural tissues at small scales: Recent progress of oscillating gradient spin echo (OGSE) neuroimaging in humans

The detection sensitivity of diffusion MRI (dMRI) is dependent on diffusion times. A shorter diffusion time can increase the sensitivity to smaller length scales. However, the conventional dMRI uses the pulse gradient spin echo (PGSE) sequence that probes relatively long diffusion times only. To overcome this, the oscillating gradient spin echo (OGSE) sequence has been developed to probe much shorter diffusion times with hardware limitations on preclinical and clinical MRI systems. The OGSE sequence has been previously used on preclinical animal MRI systems. Recently, several studies have translated the OGSE sequence to humans on clinical MRI systems and achieved new information that is invisible using conventional PGSE sequence. This paper provides an overview of the recent progress of the OGSE neuroimaging in humans, including the technical improvements in the translation of the OGSE sequence to human imaging and various applications in different neurological disorders and stroke. Some possible future directions of the OGSE sequence are also discussed.

Diffusion-time dependence of diffusional kurtosis in the mouse brain

PURPOSE

To investigate diffusion-time dependency of diffusional kurtosis in the mouse brain using pulsed-gradient spin-echo (PGSE) and oscillating-gradient spin-echo (OGSE) sequences.

METHODS

3D PGSE and OGSE kurtosis tensor data were acquired from ex vivo brains of adult, cuprizone-treated, and age-matched control mice with diffusion-time (t_D ) ~ 20 ms and frequency (f) = 70 Hz, respectively. Further, 2D acquisitions were performed at multiple times/frequencies ranging from f = 140 Hz to t_D = 30 ms with b-values up to 4000 s/mm² . Monte Carlo simulations were used to investigate the coupled effects of varying restriction size and permeability on time/frequency-dependence of kurtosis with both diffusion-encoding schemes. Simulations and experiments were further performed to investigate the effect of varying number of cycles in OGSE waveforms.

RESULTS

Kurtosis and diffusivity maps exhibited significant region-specific changes with diffusion time/frequency across both gray and white matter areas. PGSE- and OGSE-based kurtosis maps showed reversed contrast between gray matter regions in the cerebellar and cerebral cortex. Localized time/frequency-dependent changes in kurtosis tensor metrics were found in the splenium of the corpus callosum in cuprizone-treated mouse brains, corresponding to regional demyelination seen with histological assessment. Monte Carlo simulations showed that kurtosis estimates with pulsed- and oscillating-gradient waveforms differ in their sensitivity to exchange. Both simulations and experiments showed dependence of kurtosis on number of cycles in OGSE waveforms for non-zero permeability.

CONCLUSION

The results show significant time/frequency-dependency of diffusional kurtosis in the mouse brain, which can provide sensitivity to probe intrinsic cellular heterogeneity and pathological alterations in gray and white matter.

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.

Concise derivation of oscillating-gradient-derived ADC

The apparent diffusion coefficient (ADC) is analyzed for the case of oscillating diffusion-sensitizing gradients in the high-frequency regime. We provide a concise derivation of the analytical expression for the ADC for an arbitrary number of gradient oscillations N and initial phase φ. It is demonstrated that an ultimate goal - to determine the surface-to-volume ratio (S/V) from MR measurements by using oscillating gradients - can be achieved with cosine-type gradients (φ = 0) for an arbitrary N. However, to determine S/V employing gradients with φ ≠ 0 (including the sine-type gradients) and arbitrary N additionally requires prior knowledge of the time-dependent diffusion coefficient D(t). The latter is rarely known a priori but can be estimated under certain limiting conditions: (i) in the short time regime, when the total diffusion time of the measurements, t, is smaller than the characteristic diffusion time of the microstructural system of interest, an analytical expression for D(t) is available (Mitra's expression) and this allows S/V to be determined in the short time regime with sine-type gradients; (ii) in the important case of purely restricted diffusion, D(t) → 0 at sufficiently long time, the signal becomes independent of φ and behaves as for the cosine-type gradients, thus, allowing determination of S/V.

Inferring diameters of spheres and cylinders using interstitial water

OBJECT

Most early methods to infer axon diameter distributions using magnetic resonance imaging (MRI) used single diffusion encoding sequences such as pulsed gradient spin echo (SE) and are thus sensitive to axons of diameters > 5 μm. We previously simulated oscillating gradient (OG) SE sequences for diffusion spectroscopy to study smaller axons including the majority constituting cortical connections. That study suggested the model of constant extra-axonal diffusion breaks down at OG accessible frequencies. In this study we present data from phantoms to test a time-varying interstitial apparent diffusion coefficient.

MATERIALS AND METHODS

Diffusion spectra were measured in four samples from water packed around beads of diameters 3, 6 and 10 μm; and 151 μm diameter tubes. Surface-to-volume ratios, and diameters were inferred.

RESULTS

The bead pore radii estimates were 0.60±0.08 μm, 0.54±0.06 μm and 1.0±0.1 μm corresponding to bead diameters ranging from 2.9±0.4 μm to 5.3±0.7 μm, 2.6±0.3 μm to 4.8±0.6 μm, and 4.9±0.7 μm to 9±1 μm. The tube surface-to-volume ratio estimate was 0.06±0.02 μm^-1 corresponding to a tube diameter of 180±70 μm.

CONCLUSION

Interstitial models with OG inferred 3-10 μm bead diameters from 0.54±0.06 μm to 1.0±0.1 μm pore radii and 151 μm tube diameters from 0.06±0.02 μm^-1 surface-to-volume ratios.

Validation of surface-to-volume ratio measurements derived from oscillating gradient spin echo on a clinical scanner using anisotropic fiber phantoms

A diffusion measurement in the short-time surface-to-volume ratio (S/V) limit (Mitra et al., Phys Rev Lett. 1992;68:3555) can disentangle the free diffusion coefficient from geometric restrictions to diffusion. Biophysical parameters, such as the S/V of tissue membranes, can be used to estimate microscopic length scales non-invasively. However, due to gradient strength limitations on clinical MRI scanners, pulsed gradient spin echo (PGSE) measurements are impractical for probing the S/V limit. To achieve this limit on clinical systems, an oscillating gradient spin echo (OGSE) sequence was developed. Two phantoms containing 10 fiber bundles, each consisting of impermeable aligned fibers with different packing densities, were constructed to achieve a range of S/V values. The frequency-dependent diffusion coefficient, D(ω), was measured in each fiber bundle using OGSE with different gradient waveforms (cosine, stretched cosine, and trapezoidal), while D(t) was measured from PGSE and stimulated-echo measurements. The S/V values derived from the universal high-frequency behavior of D(ω) were compared against those derived from quantitative proton density measurements using single spin echo (SE) with varying echo times, and from magnetic resonance fingerprinting (MRF). S/V estimates derived from different OGSE waveforms were similar and demonstrated excellent correlation with both SE- and MRF-derived S/V measures (ρ  ≥  0.99). Furthermore, there was a smoother transition between OGSE frequency f and PGSE diffusion time when using teffS/V=9/64f, rather than the commonly used t_eff  = 1/(4f), validating the specific frequency/diffusion time conversion for this regime. Our well-characterized fiber phantom can be used for the calibration of OGSE and diffusion modeling techniques, as the S/V ratio can be measured independently using other MR modalities. Moreover, our calibration experiment offers an exciting perspective of mapping tissue S/V on clinical systems.

Pulsed and oscillating gradient MRI for assessment of cell size and extracellular space (POMACE) in mouse gliomas

Solid tumor microstructure is related to the aggressiveness of the tumor, interstitial pressure and drug delivery pathways, which are closely associated with treatment response, metastatic spread and prognosis. In this study, we introduce a novel diffusion MRI data analysis framework, pulsed and oscillating gradient MRI for assessment of cell size and extracellular space (POMACE), and demonstrate its feasibility in a mouse tumor model. In vivo and ex vivo POMACE experiments were performed on mice bearing the GL261 murine glioma model (n = 8). Since the complete diffusion time dependence is in general non-analytical, the tumor microstructure was modeled in an appropriate time/frequency regime by impermeable spheres (radius Rcell , intracellular diffusivity Dics ) surrounded by extracellular space (ECS) (approximated by constant apparent diffusivity Decs in volume fraction ECS). POMACE parametric maps (ECS, Rcell , Dics , Decs ) were compared with conventional diffusion-weighted imaging metrics, electron microscopy (EM), alternative ECS determination based on effective medium theory (EMT), and optical microscopy performed on the same samples. It was shown that Decs can be approximated by its long time tortuosity limit in the range [1/(88 Hz)-31 ms]. ECS estimations (44 ± 7% in vivo and 54 ± 11% ex vivo) were in agreement with EMT-based ECS and literature on brain gliomas. Ex vivo, ECS maps correlated well with optical microscopy. Cell sizes (Rcell  = 4.8 ± 1.3 in vivo and 4.3 ± 1.4 µm ex vivo) were consistent with EM measurements (4.7 ± 1.8 µm). In conclusion, Rcell and ECS can be quantified and mapped in vivo and ex vivo in brain tumors using the proposed POMACE method. Our experimental results support the view that POMACE provides a way to interpret the frequency or time dependence of the diffusion coefficient in tumors in terms of objective biophysical parameters of neuronal tissue, which can be used for non-invasive monitoring of preclinical cancer studies and treatment efficacy. Copyright © 2016 John Wiley & Sons, Ltd.

Symmetry of the gradient profile as second experimental dimension in the short-time expansion of the apparent diffusion coefficient as measured with NMR diffusometry

The time-dependent apparent diffusion coefficient as measured by pulsed gradient NMR can be used to estimate parameters of porous structures including the surface-to-volume ratio and the mean curvature of pores. In this work, the short-time diffusion limit and in particular the influence of the temporal profile of diffusion gradients on the expansion as proposed by Mitra et al. (1993) is investigated. It is shown that flow-compensated waveforms, i.e. those whose first moment is zero, are blind to the term linear in observation time, which is the term that is proportional to mean curvature and surface permeability. A gradient waveform that smoothly interpolates between flow-compensated and bipolar waveform is proposed and the degree of flow-compensation is used as a second experimental dimension. This two-dimensional ansatz is shown to yield an improved precision when characterizing the confining domain. This technique is demonstrated with simulations and in experiments performed with cylindrical capillaries of 100 μm radius.

Surface-to-volume ratio mapping of tumor microstructure using oscillating gradient diffusion weighted imaging

PURPOSE

To disentangle the free diffusivity (D0 ) and cellular membrane restrictions, by means of their surface-to-volume ratio (S/V), using the frequency-dependence of the diffusion coefficient D(ω), measured in brain tumors in the short diffusion-time regime using oscillating gradients (OGSE).

METHODS

In vivo and ex vivo OGSE experiments were performed on mice bearing the GL261 murine glioma model (n = 10) to identify the relevant time/frequency (t/ω) domain where D(ω) linearly decreases with ω(-1/2) . Parametric maps (S/V, D0 ) are compared with conventional DWI metrics. The impact of frequency range and temperature (20°C versus 37°C) on S/V and D0 is investigated ex vivo.

RESULTS

The validity of the short diffusion-time regime is demonstrated in vivo and ex vivo. Ex vivo measurements confirm that the purely geometric restrictions embodied in S/V are independent from temperature and frequency range, while the temperature dependence of the free diffusivity D0 is similar to that of pure water.

CONCLUSION

Our results suggest that D(ω) in the short diffusion-time regime can be used to uncouple the purely geometric restriction effect, such as S/V, from the intrinsic medium diffusivity properties, and provides a nonempirical and objective way to interpret frequency/time-dependent diffusion changes in tumors in terms of objective biophysical tissue parameters. Magn Reson Med 76:237-247, 2016. © 2015 Wiley Periodicals, Inc.

Time-Dependent Influence of Cell Membrane Permeability on MR Diffusion Measurements

PURPOSE

To investigate the influence of cell membrane permeability on diffusion measurements over a broad range of diffusion times.

METHODS

Human myelogenous leukemia K562 cells were cultured and treated with saponin to selectively alter cell membrane permeability, resulting in a broad physiologically relevant range of 0.011-0.044 μm/ms. Apparent diffusion coefficient (ADC) values were acquired with the effective diffusion time (Δeff ) ranging from 0.42 to 3000 ms. Cosine-modulated oscillating gradient spin echo (OGSE) measurements were performed to achieve short Δeff from 0.42 to 5 ms, while stimulated echo acquisitions were used to achieve long Δeff from 11 to 2999 ms. Computer simulations were also performed to support the experimental results.

RESULTS

Both computer simulations and experiments in vitro showed that the influence of membrane permeability on diffusion MR measurements is highly dependent on the choice of diffusion time, and it is negligible only when the diffusion time is at least one order of magnitude smaller than the intracellular exchange lifetime.

CONCLUSION

The influence of cell membrane permeability on the measured ADCs is negligible in OGSE measurements at moderately high frequencies. By contrast, cell membrane permeability has a significant influence on ADC and quantitative diffusion measurements at low frequencies such as those sampled using conventional pulsed gradient methods.

Fast and robust measurement of microstructural dimensions using temporal diffusion spectroscopy

Mapping axon sizes non-invasively is of interest for neuroscientists and may have significant clinical potential because nerve conduction velocity is directly dependent on axon size. Current approaches to measuring axon sizes using diffusion-weighted MRI, e.g. q-space imaging with pulsed gradient spin echo (PGSE) sequences usually require long scan times and high q-values to detect small axons (diameter <2μm). The oscillating gradient spin echo (OGSE) method has been shown to be able to achieve very short diffusion times and hence may be able to detect smaller axons with high sensitivity. In the current study, OGSE experiments were performed to measure the inner diameters of hollow microcapillaries with a range of sizes (∼1.5-19.3μm) that mimic axons in the human central nervous system. The results suggest that OGSE measurements, even with only moderately high frequencies, are highly sensitive to compartment sizes, and a minimum of two ADC values with different frequencies may be sufficient to extract the microcapillary size accurately. This suggests that the OGSE method may serve as a fast and robust measurement method for mapping axon sizes non-invasively.