In vivo investigation of restricted diffusion in the human brain with optimized oscillating diffusion gradient encoding

Single-shot diffusion trace-weighted MR spectroscopy: Comparison with unipolar and bipolar diffusion-weighted point-resolved spectroscopy

This study demonstrates the feasibility and performance of the point-resolved spectroscopy (PRESS)-based, single-shot diffusion trace-weighted sequence in quantifying the trace apparent diffusion coefficient (ADC) in phantom and in vivo using a 3-T MRI/MRS scanner. The single-shot diffusion trace-weighted PRESS sequence was implemented and compared with conventional diffusion-weighted (DW)-PRESS variants using bipolar and unipolar diffusion-sensitizing gradients. Nine phantom datasets were acquired using each sequence, and seven volunteers were scanned in three different brain regions to determine the range and variability of trace ADC values, and to allow a comparison of trace ADCs among the sequences. This sequence results in a comparatively stable range of trace ADC values that are statistically significantly higher than those produced from unipolar and bipolar DW-PRESS sequences. Only total n-acetylaspartate, total creatine, and total choline were reliably estimated in all sequences with Cramér-Rao lower bounds of, at most, 20%. The larger trace ADCs from the single-shot sequences are probably attributable to the shorter diffusion time relative to the other sequences. Overall, this study presents the first demonstration of the single-shot diffusion trace-weighted sequence in a clinical scanner at 3 T. The results show excellent agreement of phantom trace ADCs computed with all sequences, and in vivo ADCs agree well with the expected differences between gray and white matter. The diffusion trace-weighted sequence could provide an estimate of the trace ADC in a shorter scan time (by nearly a factor of 3) compared with conventional DW-PRESS approaches that require three separate orthogonal directions.

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.

Differentiating brain metastasis from glioblastoma by time-dependent diffusion MRI

BACKGROUND

This study was designed to investigate the use of time-dependent diffusion magnetic resonance imaging (MRI) parameters in distinguishing between glioblastomas and brain metastases.

METHODS

A retrospective study was conducted involving 65 patients with glioblastomas and 27 patients with metastases using a diffusion-weighted imaging sequence with oscillating gradient spin-echo (OGSE, 50 Hz) and a conventional pulsed gradient spin-echo (PGSE, 0 Hz) sequence. In addition to apparent diffusion coefficient (ADC) maps from two sequences (ADC50Hz and ADC0Hz), we generated maps of the ADC change (cADC): ADC50Hz - ADC0Hz and the relative ADC change (rcADC): (ADC50Hz - ADC0Hz)/ ADC0Hz × 100 (%).

RESULTS

The mean and the fifth and 95th percentile values of each parameter in enhancing and peritumoral regions were compared between glioblastomas and metastases. The area under the receiver operating characteristic curve (AUC) values of the best discriminating indices were compared. In enhancing regions, none of the indices of ADC0Hz and ADC50Hz showed significant differences between metastases and glioblastomas. The mean cADC and rcADC values of metastases were significantly higher than those of glioblastomas (0.24 ± 0.12 × 10-3mm2/s vs. 0.14 ± 0.03 × 10-3mm2/s and 23.3 ± 9.4% vs. 14.0 ± 4.7%; all p < 0.01). In peritumoral regions, no significant difference in all ADC indices was observed between metastases and glioblastomas. The AUC values for the mean cADC (0.877) and rcADC (0.819) values in enhancing regions were significantly higher than those for ADC0Hz5th (0.595; all p < 0.001).

CONCLUSIONS

The time-dependent diffusion MRI parameters may be useful for differentiating brain metastases from glioblastomas.

Fourier transform temporal diffusion spectroscopy

Temporal diffusion spectroscopy (TDS) currently uses the oscillating gradient spin echo (OGSE) experiment to measure the spectral density of translational velocity autocorrelation at single frequencies. Due to timing restrictions imposed by the transverse relaxation, the frequency selectivity and the sampling density of OGSE are limited, especially at low frequencies. We propose to overcome this problem by adopting the principles of Fourier transform spectroscopy. The new method of Fourier transform TDS (FTDS) uses two broadband gradient waveforms with different relative delays to make the spin echo attenuation sensitive to a broad range of diffusion frequencies with different harmonic modulations and calculates the spectrum by discrete Fourier transform. The method was validated by a measurement of diffusion spectra in highly restrictive tissues of a celery stalk and provided results consistent with OGSE, however, on a denser frequency grid.

Investigation of time-dependent diffusion in extra-axial brain tumors using oscillating-gradient spin-echo

Oscillating gradient spin-echo (OGSE) sequences provide access to short diffusion times and may provide insight into micro-scale internal structures of pathologic lesions based on an analysis of changes in diffusivity with differing diffusion times. We hypothesized that changes in diffusivity acquired with a shorter diffusion time may permit elucidation of properties related to the internal structure of extra-axial brain tumors. This study aimed to investigate the utility of changes in diffusivity between short and long diffusion times for characterizing extra-axial brain tumors. In total, 12 patients with meningothelial meningiomas, 13 patients with acoustic neuromas, and 11 patients with pituitary adenomas were scanned with a 3 T magnetic resonance imaging (MRI) scanner with diffusion-weighted imaging (DWI) using OGSE and pulsed gradient spin-echo (PGSE) (effective diffusion times [Δ_eff]: 6.5 ms and 35.2 ms) with b-values of 0 and 1000 s/mm². Relative percentage changes between shorter and longer diffusion times were calculated using region-of-interest (ROI) analysis of brain tumors on λ₁, λ₂, λ₃, and mean diffusivity (MD) maps. The diffusivities of PGSE, OGSE, and relative percentage changes were compared among each tumor type using a multiple comparisons Steel-Dwass test. The mean (standard deviation) MD at Δ_eff of 6.5 ms was 1.07 ± 0.23 10^-3 mm²/s, 1.19 ± 0.18 10^-3 mm²/s, 1.19 ± 0.21 10^-3 mm²/s for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. The mean (standard deviation) MD at Δ_eff of 35.2 ms was 0.93 ± 0.22 10^-3 mm²/s, 1.07 ± 0.19 10^-3 mm²/s, 0.82 ± 0.21 10^-3 mm²/s for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. The mean (standard deviation) of the relative percentage change was 15.7 ± 4.4%, 12.4 ± 8.2%, 46.8 ± 11.3% for meningothelial meningiomas, acoustic neuromas, and pituitary adenomas, respectively. Compared to meningiomas and acoustic neuromas, pituitary adenoma exhibited stronger diffusion time-dependence with diffusion times between 6.5 ms and 35.2 ms (P < 0.05). In conclusion, differences in diffusion time-dependence may be attributed to differences in the internal structures of brain tumors. DWI with a short diffusion time may provide additional information on the microstructure of each tumor and contribute to tumor diagnosis.

Tuned bipolar oscillating gradients for mapping frequency dispersion of diffusion kurtosis in the human brain

PURPOSE

Oscillating gradient spin-echo (OGSE) sequences have demonstrated an ability to probe time-dependent microstructural features, although they often suffer from low SNR due to increased TEs. In this work we introduce frequency-tuned bipolar (FTB) gradients as a variation of oscillating gradients with reduced TE and demonstrate their utility by mapping the frequency dispersion of kurtosis in human subjects.

METHODS

An FTB oscillating gradient waveform is presented that provides encoding of 1.5 net oscillation periods, thereby reducing the TE of the acquisition. Simulations were performed to determine an optimal protocol based on the SNR of kurtosis frequency dispersion-defined as the difference in kurtosis between pulsed and oscillating gradient acquisitions. Healthy human subjects were scanned at 7T using pulsed gradient and an optimized 23 Hz FTB protocol, which featured a maximum b-value of 2500 s/mm² . In addition, to directly compare existing methods, measurements using traditional cosine OGSE were also acquired.

RESULTS

FTB oscillating gradients demonstrated equivalent frequency-dependent diffusion measurements compared with cosine-modulated OGSE while enabling a significant reduction in TE. Optimization and in vivo results suggest that FTB gradients provide increased SNR of kurtosis dispersion maps compared with traditional cosine OGSE. The optimized FTB gradient protocol demonstrated consistent reductions in apparent kurtosis values and increased diffusivity in generated frequency dispersion maps.

CONCLUSIONS

This work presents an alternative to traditional cosine OGSE sequences, enabling more time-efficient acquisitions of frequency-dependent diffusion quantities as demonstrated through in vivo kurtosis frequency dispersion maps.

Selective Cell Size MRI Differentiates Brain Tumors from Radiation Necrosis

Brain metastasis is a common characteristic of late-stage lung cancers. High doses of targeted radiotherapy can control tumor growth in the brain but can also result in radiotherapy-induced necrosis. Current methods are limited for distinguishing whether new parenchymal lesions following radiotherapy are recurrent tumors or radiotherapy-induced necrosis, but the clinical management of these two classes of lesions differs significantly. Here, we developed, validated, and evaluated a new MRI technique termed selective size imaging using filters via diffusion times (SSIFT) to differentiate brain tumors from radiotherapy necrosis in the brain. This approach generates a signal filter that leverages diffusion time dependence to establish a cell size-weighted map. Computer simulations in silico, cultured cancer cells in vitro, and animals with brain tumors in vivo were used to comprehensively validate the specificity of SSIFT for detecting typical large cancer cells and the ability to differentiate brain tumors from radiotherapy necrosis. SSIFT was also implemented in patients with metastatic brain cancer and radiotherapy necrosis. SSIFT showed high correlation with mean cell sizes in the relevant range of less than 20 μm. The specificity of SSIFT for brain tumors and reduced contrast in other brain etiologies allowed SSIFT to differentiate brain tumors from peritumoral edema and radiotherapy necrosis. In conclusion, this new, cell size-based MRI method provides a unique contrast to differentiate brain tumors from other pathologies in the brain.

SIGNIFICANCE

This work introduces and provides preclinical validation of a new diffusion MRI method that exploits intrinsic differences in cell sizes to distinguish brain tumors and radiotherapy necrosis.

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 (K_iso), variance due to diffusion anisotropy (K_aniso), and microscopic kurtosis (μK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, μK is typically ignored in diffusion MRI signal modeling 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 gray 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.

Time-Dependent Diffusion MRI for Quantitative Microstructural Mapping of Prostate Cancer

Background Recently developed time-dependent diffusion MRI has potential in characterizing cellular tissue microstructures; however, its value in imaging prostate cancer (PCa) remains unknown. Purpose To investigate the feasibility of time-dependent diffusion MRI-based microstructural mapping for noninvasively characterizing cellular properties of PCa and for discriminating between clinically significant PCa and clinically insignificant disease. Materials and Methods Men with a clinical suspicion of PCa were enrolled prospectively between October 2019 and August 2020. Time-dependent diffusion MRI data were acquired with pulsed and oscillating gradient diffusion MRI sequences at an equivalent diffusion time of 7.5-30 msec on a 3.0-T scanner. Time-dependent diffusion MRI-based microstructural parameters, including cell diameter, intracellular volume fraction, cellularity, and diffusivities, were estimated with a two-compartment model. These were compared for different International Society of Urological Pathology grade groups (GGs), and their performance in discriminating clinically significant PCa (GG >1) from clinically insignificant disease (benign and GG 1) was determined with a linear discriminant analysis. The fitted microstructural parameters were validated by means of correlation with histopathologic measurements. Results In the 48 enrolled men, the time-dependent diffusion MRI measurements showed that higher GG was correlated with higher intracellular volume fraction and higher cellularity (intracellular volume fraction = 0.22, 0.36, 0.34, 0.37, and 0.40 in GGs 1-5, respectively; P < .001 at one-way analysis of variance), while lower cell diameter was found at higher GGs (diameter = 23.4, 18.3, 19.2, 17.9, and 18.5 μm in GGs 1-5, respectively; P = .002). Among all measurements derived from time-dependent diffusion MRI, cellularity achieved the highest diagnostic performance, with an accuracy of 92% (44 of 48 participants) and area under the receiver operating characteristic curve of 0.96 (95% CI: 0.87, 0.99) in discriminating clinically significant PCa from clinically insignificant disease. Microstructural mapping was supported by positive correlations between time-dependent diffusion MRI-based and pathologic examination-based intracellular volume fraction (r = 0.83; P < .001). Conclusion Time-dependent diffusion MRI-based microstructural mapping correlates with pathologic findings and demonstrates promise for characterizing prostate cancer. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Chatterjee and Oto in this issue.

Evaluating diffusion dispersion across an extended range of b-values and frequencies: Exploiting gap-filled OGSE shapes, strong gradients, and spiral readouts

Purpose

To address the long echo times and relatively weak diffusion sensitization that typically limit oscillating gradient spin‐echo (OGSE) experiments, an OGSE implementation combining spiral readouts, gap‐filled oscillating gradient shapes providing stronger diffusion encoding, and a high‐performance gradient system is developed here and utilized to investigate the tradeoff between b‐value and maximum OGSE frequency in measurements of diffusion dispersion (i.e., the frequency dependence of diffusivity) in the in vivo human brain. In addition, to assess the effects of the marginal flow sensitivity introduced by these OGSE waveforms, flow‐compensated variants are devised for experimental comparison.

Methods

Using DTI sequences, OGSE acquisitions were performed on three volunteers at b‐values of 300, 500, and 1000 s/mm² and frequencies up to 125, 100, and 75 Hz, respectively; scans were performed for gap‐filled oscillating gradient shapes with and without flow sensitivity. Pulsed gradient spin‐echo DTI acquisitions were also performed at each b‐value. Upon reconstruction, mean diffusivity (MD) maps and maps of the diffusion dispersion rate were computed.

Results

The power law diffusion dispersion model was found to fit best to MD measurements acquired at b = 1000 s/mm² despite the associated reduction of the spectral range; this observation was consistent with Monte Carlo simulations. Furthermore, diffusion dispersion rates without flow sensitivity were slightly higher than flow‐sensitive measurements.

Conclusion

The presented OGSE implementation provided an improved depiction of diffusion dispersion and demonstrated the advantages of measuring dispersion at higher b‐values rather than higher frequencies within the regimes employed in this study.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.

Characterization and correction of time-varying eddy currents for diffusion MRI

PURPOSE

To develop and test a method for reducing artifacts due to time-varying eddy currents in oscillating gradient spin-echo (OGSE) diffusion images.

METHODS

An in-house algorithm (TVEDDY), that for the first time retrospectively models eddy current decay, was tested on pulsed gradient spin echo and OGSE brain images acquired at 7 T. Image pairs were acquired using opposite polarity diffusion gradients. A three-parameter exponential decay model (two amplitudes and a time constant) was used to characterize and correct eddy current distortions by minimizing the intensity difference between image pairs. Correction performance was compared with conventional correction methods by evaluating the mean squared error (MSE) between diffusion-weighted images acquired with opposite polarity diffusion gradients. As a ground-truth comparison, images were corrected using field dynamics up to third order in space, measured using a field monitoring system.

RESULTS

Time-varying eddy currents were observed for OGSE, which introduced blurring that was not reduced using the traditional approach but was diminished considerably with TVEDDY and field monitoring-informed model-based reconstruction. No MSE difference was observed between the conventional approach and TVEDDY for pulsed gradient spin echo, but for OGSE TVEDDY resulted in significantly lower MSE than the conventional approach. The field-monitoring reconstruction had the lowest MSE for both pulsed gradient spin echo and OGSE.

CONCLUSION

This work establishes that it is possible to estimate time-varying eddy currents from the actual diffusion data, which provides substantial image-quality improvements for gradient-intensive diffusion MRI acquisitions like OGSE.

Time-dependent Diffusion in Brain Abscesses Investigated with Oscillating-gradient Spin-echo

Oscillating-gradient spin-echo sequences enable the measurement of diffusion weighting with a short diffusion time and can provide indications of internal structures. We report two cases of brain abscess in which the apparent diffusion coefficient (ADC) values appear higher at short diffusion times in comparison with those at long diffusion times. Diffusion time dependence of the ADC in brain abscesses suggests not only substrate viscosity but also restricted diffusion due to the structure within the lesions.

Feasibility of oscillating and pulsed gradient diffusion MRI to assess neonatal hypoxia-ischemia on clinical systems

Diffusion-time- (t_d) dependent diffusion MRI (dMRI) extends our ability to characterize brain microstructure by measuring dMRI signals at varying t_d. The use of oscillating gradient (OG) is essential for accessing short t_d but is technically challenging on clinical MRI systems. This study aims to investigate the clinical feasibility and value of t_d-dependent dMRI in neonatal hypoxic-ischemic encephalopathy (HIE). Eighteen HIE neonates and six normal term-born neonates were scanned on a 3 T scanner, with OG-dMRI at an oscillating frequency of 33 Hz (equivalent t_d ≈ 7.5 ms) and pulsed gradient (PG)-dMRI at a t_d of 82.8 ms and b-value of 700 s/mm². The t_d-dependence, as quantified by the difference in apparent diffusivity coefficients between OG- and PG-dMRI (ΔADC), was observed in the normal neonatal brains, and the ΔADC was higher in the subcortical white matter than the deep grey matter. In HIE neonates with severe and moderate injury, ΔADC significantly increased in the basal ganglia (BG) compared to the controls (23.7% and 10.6%, respectively). In contrast, the conventional PG-ADC showed a 12.6% reduction only in the severe HIE group. White matter edema regions also demonstrated increased ΔADC, where PG-ADC did not show apparent changes. Our result demonstrated that t_d-dependent dMRI provided high sensitivity in detecting moderate-to-severe HIE.

Time-dependent diffusion MRI probes cerebellar microstructural alterations in a mouse model of Down syndrome

The cerebellum is a complex system with distinct cortical laminar organization. Alterations in cerebellar microstructure are common and associated with many factors such as genetics, cancer and ageing. Diffusion MRI (dMRI) provides a non-invasive tool to map the brain structural organization, and the recently proposed diffusion-time (td )-dependent dMRI further improves its capability to probe the cellular and axonal/dendritic microstructures by measuring water diffusion at multiple spatial scales. The td -dependent diffusion profile in the cerebellum and its utility in detecting cerebellar disorders, however, are not yet elucidated. Here, we first deciphered the spatial correspondence between dMRI contrast and cerebellar layers, based on which the cerebellar layer-specific td -dependent dMRI patterns were characterized in both euploid and Ts65Dn mice, a mouse model of Down syndrome. Using oscillating gradient dMRI, which accesses diffusion at short td 's by modulating the oscillating frequency, we detected subtle changes in the apparent diffusivity coefficient of the cerebellar internal granular layer and Purkinje cell layer of Ts65Dn mice that were not detectable by conventional pulsed gradient dMRI. The detection sensitivity of oscillating gradient dMRI increased with the oscillating frequency at both the neonatal and adult stages. The td -dependence, quantified by ΔADC map, was reduced in Ts65Dn mice, likely associated with the reduced granule cell density and abnormal dendritic arborization of Purkinje cells as revealed from histological evidence. Our study demonstrates superior sensitivity of short-td diffusion using oscillating gradient dMRI to detect cerebellar microstructural changes in Down syndrome, suggesting the potential application of this technique in cerebellar disorders.

MR cell size imaging with temporal diffusion spectroscopy

Cytological features such as cell size and intracellular morphology provide fundamental information on cell status and hence may provide specific information on changes that arise within biological tissues. Such information is usually obtained by invasive biopsy in current clinical practice, which suffers several well-known disadvantages. Recently, novel MRI methods such as IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) have been developed for direct measurements of mean cell size non-invasively. The IMPULSED protocol is based on using temporal diffusion spectroscopy (TDS) to combine measurements of water diffusion over a wide range of diffusion times to probe cellular microstructure over varying length scales. IMPULSED has been shown to provide rapid, robust, and reliable mapping of mean cell size and is suitable for clinical imaging. More recently, cell size distributions have also been derived by appropriate analyses of data acquired with IMPULSED or similar sequences, which thus provides MRI-cytometry. This review summarizes the basic principles, practical implementations, validations, and example applications of MR cell size imaging based on TDS and demonstrates how cytometric information can be used in various applications. In addition, the limitations and potential future directions of MR cytometry are identified including the diagnosis of nonalcoholic steatohepatitis of the liver and the assessment of treatment response of cancers.

A simulation study of cell size and volume fraction mapping for tissue with two underlying cell populations using diffusion-weighted MRI

PURPOSE

To propose a method for voxel-wise estimation of cell radii and volume fractions of two cell populations when they coexist in the same MR voxel using the combination of diffusion-weighted MRI and microstructural modeling.

METHOD

Microstructure models were investigated using diffusion data simulated with the matrix method for a range of microstructures mimicking tumor tissue with two cell populations, using acquisition parameters available on preclinical scanners. The effect of noise was investigated for a subset of these microstructures. The accuracy and precision of the estimated radii and volume fractions for large and small cells R l , R s , v i n , l , v i n , s were evaluated by comparing the estimates to their true values. The stability of model fitting was characterized by the percentage of accepted fits.

RESULTS

The estimation accuracy and precision, and thus the ability to robustly distinguish the two cell populations, depended on the microstructural properties and SNR. For a SNR of 50, a minimum difference of 3 μm between the radius of the large and small cell populations was required for differentiation. Proposed modifications to the two cell population microstructure model, including constrained fits, improved the stability of fits.

CONCLUSIONS

This proof-of-concept study proposed a diffusion MRI-based method for voxel-wise estimation of cell radii and volume fractions of two cell populations when they coexist in the same MR voxel. The ability to reliably characterize tissue with two cell populations opens exciting avenues of potential applications in both tumor diagnosis and treatment monitoring.

Diffusion Tensor Imaging of Skeletal Muscle Contraction Using Oscillating Gradient Spin Echo

Diffusion tensor imaging (DTI) measures water diffusion in skeletal muscle tissue and allows for muscle assessment in a broad range of neuromuscular diseases. However, current DTI measurements, typically performed using pulsed gradient spin echo (PGSE) diffusion encoding, are limited to the assessment of non-contracted musculature, therefore providing limited insight into muscle contraction mechanisms and contraction abnormalities. In this study, we propose the use of an oscillating gradient spin echo (OGSE) diffusion encoding strategy for DTI measurements to mitigate the effect of signal voids in contracted muscle and to obtain reliable diffusivity values. Two OGSE sequences with encoding frequencies of 25 and 50 Hz were tested in the lower leg of five healthy volunteers with relaxed musculature and during active dorsiflexion and plantarflexion, and compared with a conventional PGSE approach. A significant reduction of areas of signal voids using OGSE compared with PGSE was observed in the tibialis anterior for the scans obtained in active dorsiflexion and in the soleus during active plantarflexion. The use of PGSE sequences led to unrealistically elevated axial diffusivity values in the tibialis anterior during dorsiflexion and in the soleus during plantarflexion, while the corresponding values obtained using the OGSE sequences were significantly reduced. Similar findings were seen for radial diffusivity, with significantly higher diffusivity measured in plantarflexion in the soleus muscle using the PGSE sequence. Our preliminary results indicate that DTI with OGSE diffusion encoding is feasible in human musculature and allows to quantitatively assess diffusion properties in actively contracting skeletal muscle. OGSE holds great potential to assess microstructural changes occurring in the skeletal muscle during contraction, and for non-invasive assessment of contraction abnormalities in patients with muscle diseases.

Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding

Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.

Magic DIAMOND: Multi-fascicle diffusion compartment imaging with tensor distribution modeling and tensor-valued diffusion encoding

Diffusion tensor imaging provides increased sensitivity to microstructural tissue changes compared to conventional anatomical imaging but also presents limited specificity. To tackle this problem, the DIAMOND model subdivides the voxel content into diffusion compartments and draws from diffusion-weighted data to estimate compartmental non-central matrix-variate Gamma distributions of diffusion tensors. It models each sub-voxel fascicle separately, resolving crossing white-matter pathways and allowing for a fascicle-element (fixel) based analysis of microstructural features. Alternatively, specific features of the intra-voxel diffusion tensor distribution can be selectively measured using tensor-valued diffusion-weighted acquisition schemes. However, the impact of such schemes on estimating brain microstructural features has only been studied in a handful of parametric single-fascicle models. In this work, we derive a general Laplace transform for the non-central matrix-variate Gamma distribution, which enables the extension of DIAMOND to tensor-valued encoded data. We then evaluate this "Magic DIAMOND" model in silico and in vivo on various combinations of tensor-valued encoded data. Assessing uncertainty on parameter estimation via stratified bootstrap, we investigate both voxel-based and fixel-based metrics by carrying out multi-peak tractography. We demonstrate using in silico evaluations that tensor-valued diffusion encoding significantly improves Magic DIAMOND's accuracy. Most importantly, we show in vivo that our estimated metrics can be robustly mapped along tracks across regions of fiber crossing, which opens new perspectives for tractometry and microstructure mapping along specific white-matter tracts.

Improved gradient waveforms for oscillating gradient spin-echo (OGSE) diffusion tensor imaging

The dependence of the diffusion tensor on frequency is of great interest in studies of tissue microstructure because it reveals restrictions to the Brownian motion of water molecules caused by cell membranes. Oscillating gradient spin-echo (OGSE) sequences can sample this dependence with gradient shapes for which the power spectrum of the zeroth moment is focused at a target frequency. In order to maintain the total spectral power (ie the b-value), oscillating gradient amplitudes must grow with the frequency squared. For this reason, OGSE applications on clinical MRI scanners are limited to low frequencies, for which it is difficult to obtain a narrow frequency bandwidth of the gradient moment in a useful echo time. In particular, the commonly used pair of single-period trapezoidal-cosine pulses separated by a half-period produces significant side lobes away from the target frequency. To mitigate this effect, improved OGSE waveforms are proposed, which reduce the gap between the two gradient pulses to the minimum duration required for the refocusing RF pulse. Additionally, a slight deviation from the periodicity of the waveforms is proposed in order to permit using the maximum slew rate of the gradient system for all lobes of trapezoidal waveforms while maintaining advantageous spectral properties, which is not the case for the currently used OGSE sequences. Numerical calculations validate these changes, showing that both modifications significantly narrow the gradient moment power spectrum and increase the contribution of its main lobe to the b-value, thus improving the specificity of the measurement. The utility of the new shapes is demonstrated by diffusion tensor measurements of human white matter in vivo over the range of 30-75 Hz with a b-value of nearly 1000 s/mm² , using a high-performance gradient insert. However, the improvement should increase the sampling precision of OGSE experiments for all gradient systems.

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.

Accuracy and precision of statistical descriptors obtained from multidimensional diffusion signal inversion algorithms

In biological tissues, typical MRI voxels comprise multiple microscopic environments, the local organization of which can be captured by microscopic diffusion tensors. The measured diffusion MRI signal can, therefore, be written as the multidimensional Laplace transform of an intravoxel diffusion tensor distribution (DTD). Tensor-valued diffusion encoding schemes have been designed to probe specific features of the DTD, and several algorithms have been introduced to invert such data and estimate statistical descriptors of the DTD, such as the mean diffusivity, the variance of isotropic diffusivities, and the mean squared diffusion anisotropy. However, the accuracy and precision of these estimations have not been assessed systematically and compared across methods. In this article, we perform and compare such estimations in silico for a one-dimensional Gamma fit, a generalized two-term cumulant approach, and two-dimensional and four-dimensional Monte-Carlo-based inversion techniques, using a clinically feasible tensor-valued acquisition scheme. In particular, we compare their performance at different signal-to-noise ratios (SNRs) for voxel contents varying in terms of the aforementioned statistical descriptors, orientational order, and fractions of isotropic and anisotropic components. We find that all inversion techniques share similar precision (except for a lower precision of the two-dimensional Monte Carlo inversion) but differ in terms of accuracy. While the Gamma fit exhibits infinite-SNR biases when the signal deviates strongly from monoexponentiality and is unaffected by orientational order, the generalized cumulant approach shows infinite-SNR biases when this deviation originates from the variance in isotropic diffusivities or from the low orientational order of anisotropic diffusion components. The two-dimensional Monte Carlo inversion shows remarkable accuracy in all systems studied, given that the acquisition scheme possesses enough directions to yield a rotationally invariant powder average. The four-dimensional Monte Carlo inversion presents no infinite-SNR bias, but suffers significantly from noise in the data, while preserving good contrast in most systems investigated.

Diffusion-prepared 3D gradient spin-echo sequence for improved oscillating gradient diffusion MRI

PURPOSE

Oscillating gradient (OG) enables the access of short diffusion times for time-dependent diffusion MRI (dMRI); however, it poses several technical challenges for clinical use. This study proposes a 3D oscillating gradient-prepared gradient spin-echo (OGprep-GRASE) sequence to improve SNR and shorten acquisition time for OG dMRI on clinical scanners.

METHODS

The 3D OGprep-GRASE sequence consisted of global saturation, diffusion encoding, fat saturation, and GRASE readout modules. Multiplexed sensitivity-encoding reconstruction was used to correct the phase errors between multiple shots. We compared the scan time and SNR of the proposed sequence and the conventional 2D-EPI sequence for OG dMRI at 30-90-mm slice coverage. We also examined the time-dependent diffusivity changes with OG dMRI acquired at frequencies of 50 Hz and 25 Hz and pulsed-gradient dMRI at diffusion times of 30 ms and 60 ms.

RESULTS

The OGprep-GRASE sequence reduced the scan time by a factor of 1.38, and increased the SNR by 1.74-2.27 times compared with 2D EPI for relatively thick slice coverage (60-90 mm). The SNR gain led to improved diffusion-tensor reconstruction in the multishot protocols. Image distortion in 2D-EPI images was also reduced in GRASE images. Diffusivity measurements from the pulsed-gradient dMRI and OG dMRI showed clear diffusion-time dependency in the white matter and gray matter of the human brain, using both the GRASE and EPI sequences.

CONCLUSION

The 3D OGprep-GRASE sequence improved scan time and SNR and reduced image distortion compared with the 2D multislice acquisition for OG dMRI on a 3T clinical system, which may facilitate the clinical translation of time-dependent dMRI.

Differentiation of high-grade and low-grade intra-axial brain tumors by time-dependent diffusion MRI

INTRODUCTION

Oscillating gradient spin-echo (OGSE) sequences enable acquisitions with shorter diffusion times. There is growing interest in the effect of diffusion time on apparent diffusion coefficient (ADC) values in patients with cancer. However, little evidence exists regarding its usefulness for differentiating between high-grade and low-grade brain tumors. The purpose of this study is to investigate the utility of changes in the ADC value between short and long diffusion times in distinguishing low-grade and high-grade brain tumors.

MATERIAL AND METHODS

Eleven patients with high-grade brain tumors and ten patients with low-grade brain tumors were scanned using a 3 T magnetic resonance imaging with diffusion-weighted imaging (DWI) using OGSE and PGSE (effective diffusion time [Δ_eff]: 6.5 ms and 35.2 ms) and b-values of 0 and 1000 s/mm². Using a region of interest (ROI) analysis of the brain tumors, we measured the ADC for two Δeff (ADC_Δeff) values and computed the subtraction ADC (ΔADC = ADC_{6.5 ms} - ADC_{35.2 ms}) and the relative ADC (ΔADC = (ADC_{6.5 ms} - ADC_{35.2 ms}) / ADC_{35.2 ms} × 100). The maximum values for the subtraction ADC (ΔADC_max) and the relative ADC (rADC_max) on the ROI were compared between low-grade and high-grade tumors using the Wilcoxon rank-sum test. A P-value <.05 was considered significant. The ROIs were also placed in the normal white matter of patients with high- and low-grade brain tumors, and ΔADC_max values were determined.

RESULTS

High-grade tumors had significantly higher ΔADC_max and rADC_max than low-grade tumors. The ΔADC_max values of the normal white matter were lower than the ΔADC_max of high- and low-grade brain tumors.

CONCLUSION

The dependence of ADC values on diffusion time between 6.5 ms and 35.2 ms was stronger in high-grade tumors than in low-grade tumors, suggesting differences in internal tissue structure. This finding highlights the importance of reporting diffusion times in ADC evaluations and might contribute to the grading of brain tumors using DWI.

Magnetic resonance imaging of mean cell size in human breast tumors

PURPOSE

Cell size is a fundamental characteristic of all tissues, and changes in cell size in cancer reflect tumor status and response to treatments, such as apoptosis and cell-cycle arrest. Unfortunately, cell size can currently be obtained only by pathological evaluation of tumor tissue samples obtained invasively. Previous imaging approaches are limited to preclinical MRI scanners or require relatively long acquisition times that are impractical for clinical imaging. There is a need to develop cell-size imaging for clinical applications.

METHODS

We propose a clinically feasible IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) approach that can characterize mean cell sizes in solid tumors. We report the use of a combination of pulse sequences, using different gradient waveforms implemented on clinical MRI scanners and analytical equations based on these waveforms to analyze diffusion-weighted MRI signals and derive specific microstructural parameters such as cell size. We also describe comprehensive validations of this approach using computer simulations, cell experiments in vitro, and animal experiments in vivo and demonstrate applications in preoperative breast cancer patients.

RESULTS

With fast acquisitions (~7 minutes), IMPULSED can provide high-resolution (1.3 mm in-plane) mapping of mean cell size of human tumors in vivo on clinical 3T MRI scanners. All validations suggest that IMPULSED provides accurate and reliable measurements of mean cell size.

CONCLUSION

The proposed IMPULSED method can assess cell-size variations in tumors of breast cancer patients, which may have the potential to assess early response to neoadjuvant therapy.

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.

Oscillating diffusion-encoding with a high gradient-amplitude and high slew-rate head-only gradient for human brain imaging

PURPOSE

We investigate the importance of high gradient-amplitude and high slew-rate on oscillating gradient spin echo (OGSE) diffusion imaging for human brain imaging and evaluate human brain imaging with OGSE on the MAGNUS head-gradient insert (200 mT/m amplitude and 500 T/m/s slew rate).

METHODS

Simulations with cosine-modulated and trapezoidal-cosine OGSE at various gradient amplitudes and slew rates were performed. Six healthy subjects were imaged with the MAGNUS gradient at 3T with OGSE at frequencies up to 100 Hz and b = 450 s/mm² . Comparisons were made against standard pulsed gradient spin echo (PGSE) diffusion in vivo and in an isotropic diffusion phantom.

RESULTS

Simulations show that to achieve high frequency and b-value simultaneously for OGSE, high gradient amplitude, high slew rates, and high peripheral nerve stimulation limits are required. A strong linear trend for increased diffusivity (mean: 8-19%, radial: 9-27%, parallel: 8-15%) was observed in normal white matter with OGSE (20 Hz to 100 Hz) as compared to PGSE. Linear fitting to frequency provided excellent correlation, and using a short-range disorder model provided radial long-term diffusivities of D_∞,MD = 911 ± 72 µm² /s, D_∞,PD = 1519 ± 164 µm² /s, and D_∞,RD = 640 ± 111 µm² /s and correlation lengths of l_{c ,MD} = 0.802 ± 0.156 µm, l_{c ,PD} = 0.837 ± 0.172 µm, and l_{c ,RD} = 0.780 ± 0.174 µm. Diffusivity changes with OGSE frequency were negligible in the phantom, as expected.

CONCLUSION

The high gradient amplitude, high slew rate, and high peripheral nerve stimulation thresholds of the MAGNUS head-gradient enables OGSE acquisition for in vivo human brain imaging.

Diffusion time dependency along the human corpus callosum and exploration of age and sex differences as assessed by oscillating gradient spin-echo diffusion tensor imaging

Conventional diffusion imaging uses pulsed gradient spin echo (PGSE) waveforms with diffusion times of tens of milliseconds (ms) to infer differences of white matter microstructure. The combined use of these long diffusion times with short diffusion times (<10 ​ms) enabled by oscillating gradient spin echo (OGSE) waveforms can enable more sensitivity to changes of restrictive boundaries on the scale of white matter microstructure (e.g. membranes reflecting the axon diameters). Here, PGSE and OGSE images were acquired at 4.7 ​T from 20 healthy volunteers aged 20-73 years (10 males). Mean, radial, and axial diffusivity, as well as fractional anisotropy were calculated in the genu, body and splenium of the corpus callosum (CC). Monte Carlo simulations were also conducted to examine the relationship of intra- and extra-axonal radial diffusivity with diffusion time over a range of axon diameters and distributions. The results showed elevated diffusivities with OGSE relative to PGSE in the genu and splenium (but not the body) in both males and females, but the OGSE-PGSE difference was greater in the genu for males. Females showed positive correlations of OGSE-PGSE diffusivity difference with age across the CC, whereas there were no such age correlations in males. Simulations of radial diffusion demonstrated that for axon sizes in human brain both OGSE and PGSE diffusivities were dominated by extra-axonal water, but the OGSE-PGSE difference nonetheless increased with area-weighted outer-axon diameter. Therefore, the lack of OGSE-PGSE difference in the body is not entirely consistent with literature that suggests it is composed predominantly of axons with large diameter. The greater OGSE-PGSE difference in the genu of males could reflect larger axon diameters than females. The OGSE-PGSE difference correlation with age in females could reflect loss of smaller axons at older ages. The use of OGSE with short diffusion times to sample the microstructural scale of restriction implies regional differences of axon diameters along the corpus callosum with preliminary results suggesting a dependence on age and sex.

Diffusion dispersion imaging: Mapping oscillating gradient spin-echo frequency dependence in the human brain

PURPOSE

Oscillating gradient spin-echo (OGSE) diffusion MRI provides information about the microstructure of biological tissues by means of the frequency dependence of the apparent diffusion coefficient (ADC). ADC dependence on OGSE frequency has been explored in numerous rodent studies, but applications in the human brain have been limited and have suffered from low contrast between different frequencies, long scan times, and a limited exploration of the nature of the ADC dependence on frequency.

THEORY AND METHODS

Multiple frequency OGSE acquisitions were acquired in healthy subjects at 7T to explore the power-law frequency dependence of ADC, the "diffusion dispersion." Furthermore, a method for optimizing the estimation of the ADC difference between different OGSE frequencies was developed, which enabled the design of a highly efficient protocol for mapping diffusion dispersion.

RESULTS

For the first time, evidence of a linear dependence of ADC on the square root of frequency in healthy human white matter was obtained. Using the optimized protocol, high-quality, full-brain maps of apparent diffusion dispersion rate were also demonstrated at an isotropic resolution of 2 mm in a scan time of 6 min.

CONCLUSIONS

This work sheds light on the nature of diffusion dispersion in the healthy human brain and introduces full-brain diffusion dispersion mapping at clinically relevant scan times. These advances may lead to new biomarkers of pathology or improved microstructural modeling.

Brain lateralization probed by water diffusion at the atomic to micrometric scale

Combined neutron scattering and diffusion nuclear magnetic resonance experiments have been used to reveal significant interregional asymmetries (lateralization) in bovine brain hemispheres in terms of myelin arrangement and water dynamics at micron to atomic scales. Thicker myelin sheaths were found in the left hemisphere using neutron diffraction. 4.7 T dMRI and quasi-elastic neutron experiments highlighted significant differences in the properties of water dynamics in the two hemispheres. The results were interpreted in terms of hemisphere-dependent cellular composition (number of neurons, cell distribution, etc.) as well as specificity of neurological functions (such as preferential networking).

Comparison of pulsed and oscillating gradient diffusion-weighted MRI for characterizing hepatocellular nodules in liver cirrhosis: ex vivo study in a rat model

BACKGROUND

In contrast to classical pulsed gradient diffusion-weighted MRI, oscillating gradient diffusion-weighted MR imaging (DWI) is sensitive to short distance diffusion changes at the intracellular level.

PURPOSE

To compare the diagnostic performance of pulsed and oscillating DWI for characterizing hepatocellular nodules in a rat model of hepatic cirrhosis.

STUDY TYPE

Prospective, experimental study.

ANIMAL MODEL

Cirrhosis was induced by weekly intraperitoneal injection of diethylnitrosamine in Wistar rats.

FIELD STRENGTH/SEQUENCE

Ex vivo liver MRI was performed at 7T with T₁ -weighted, T₂ -weighted, pulsed, and oscillating gradient diffusion-weighted sequences.

ASSESSMENT

Apparent diffusion coefficient from pulsed (ADC_pulsed ) and oscillating gradient (ADC_oscillating ) sequences was calculated in 82 nodules identified on the T₁ /T₂ -weighted images and on pathological examination. Two pathologists classified the nodules in three categories: benign (regenerative and low-grade dysplastic nodules), with intermediate malignancy (high-grade dysplastic nodules and early hepatocellular carcinomas) and overtly malignant (progressed hepatocellular carcinomas).

STATISTICAL TESTS

Differences between groups were assessed with Kruskal-Wallis and Mann-Whitney tests.

RESULTS

ADC, mainly ADC_oscillating , increased in the group of nodules with intermediate malignancy (ADC_pulsed : 0.75 ± 0.25 × 10^-3 mm² /s vs. 0.64 ± 0.07 × 10^-3 mm² /s in benign nodules, P = 0.025; ADC_oscillating : 0.81 ± 0.20 × 10^-3 mm² /s vs. 0.65 ± 0.13 × 10^-3 mm² /s, P = 0.0008) and ADC_pulsed decreased in the group of progressed hepatocellular carcinomas (ADC_pulsed : 0.60 ± 0.08 × 10^-3 mm² /s, P = 0.042; ADC_oscillating : 0.68 ± 0.08 × 10^-3 mm² /s, P = 0.1).

DATA CONCLUSION

ADC during hepatocarcinogenesis in rats increased in nodules with intermediate malignancy and decreased in progressed hepatocellular carcinomas. Our results suggest that oscillating gradient DWI is more sensitive to the early steps of hepatocarcinogenesis and might be useful for differentiating between high-grade dysplastic nodules / early hepatocellular carcinomas and regenerating nodules / low-grade dysplastic nodules.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2020;51:1065-1074.

Evidence of the diffusion time dependence of intravoxel incoherent motion in the brain

PURPOSE

To investigate the diffusion time (T_D ) dependence of intravoxel incoherent motion (IVIM) signals in the brain.

METHODS

A 3-compartment IVIM model was proposed to characterize 2 types of microcirculatory flows in addition to tissue water in the brain: flows that cross multiple vascular segments (pseudo-diffusive) and flows that stay in 1 segment (ballistic) within T_D . The model was first evaluated using simulated flow signals. Experimentally, flow-compensated (FC) pulsed-gradient spin-echo (PGSE) and oscillating-gradient spin-echo (OGSE) sequences were tested using a flow phantom and then used to examine IVIM signals in the mouse brain with T_D ranging from ~2.5 ms to 40 ms on an 11.7T scanner.

RESULTS

By fitting the model to simulated flow signals, we demonstrated the T_D dependency of the estimated fraction of pseudo-diffusive flow and the pseudo-diffusion coefficient (D*), which were dictated by the characteristic timescale of microcirculatory flow (τ). Flow phantom experiments validated that the OGSE and FC-PGSE sequences were not susceptible to the change in flow velocity. In vivo mouse brain data showed that both the estimated fraction of pseudo-diffusive flow and D* increased significantly as T_D increased.

CONCLUSION

We demonstrated that IVIM signals measured in the brain are T_D -dependent, potentially because more microcirculatory flows approach the pseudo-diffusive limit as T_D increases with respect to τ. Measuring the T_D dependency of IVIM signals may provide additional information on microvascular flows in the brain.

Oscillating-gradient diffusion magnetic resonance imaging detects acute subcellular structural changes in the mouse forebrain after neonatal hypoxia-ischemia

The recently developed oscillating-gradient diffusion MRI (OG-dMRI) technique extends our ability to examine brain structures at different spatial scales. In this study, we investigated the sensitivity of OG-dMRI in detecting cellular and subcellular structural changes in a mouse model of neonatal hypoxia ischemia (HI). Neonatal mice received unilateral HI injury or sham injury at postnatal day 10, followed by in vivo T₂-weighted and diffusion MRI of the brains at 3-6 h and 24 h after HI. Apparent diffusion coefficient (ADC) maps were acquired using conventional pulsed-gradient dMRI (PG-dMRI) and OG-dMRI with oscillating frequencies from 50 to 200 Hz. Pathology at cellular and subcellular levels was evaluated using neuronal, glial, and mitochondrial markers. We found significantly higher rates of ADC increase with oscillating frequencies (Δ_fADC) in the ipsilateral edema region, compared to the contralateral side, starting as early as 3 h after HI. Even in injured regions that showed no apparent change in PG-ADC or pseudo-normalized PG-ADC measurements, Δ_fADC remained significantly elevated. Histopathology showed swelling of sub-cellular structures in these regions with no apparent whole-cell level change. These results suggest that OG-dMRI is sensitive to subcellular structural changes in the brain after HI and is less susceptible to pseudo-normalization than PG-dMRI.

Diffusion tensor distribution imaging

Conventional diffusion MRI yields voxel-averaged parameters that suffer from ambiguities for heterogeneous anisotropic materials such as brain tissue. Using principles from solid-state NMR spectroscopy, we have previously introduced the shape of the diffusion encoding tensor as a separate acquisition dimension that disentangles isotropic and anisotropic contributions to the observed diffusivities, thereby allowing for unconstrained data inversion into diffusion tensor distributions with "size," "shape," and orientation dimensions. Here we combine our recent non-parametric data inversion algorithm and data acquisition protocol with an imaging pulse sequence to demonstrate spatial mapping of diffusion tensor distributions using a previously developed composite phantom with multiple isotropic and anisotropic components. We propose a compact format for visualizing two-dimensional arrays of the distributions, new scalar parameters quantifying intra-voxel heterogeneity, and a binning procedure giving maps of all relevant parameters for each of the components resolved in the multidimensional distribution space.

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.

Effects of nongaussian diffusion on "isotropic diffusion" measurements: An ex-vivo microimaging and simulation study

Designing novel diffusion-weighted pulse sequences to probe tissue microstructure beyond the conventional Stejskal-Tanner family is currently of broad interest. One such technique, multidimensional diffusion MRI, has been recently proposed to afford model-free decomposition of diffusion signal kurtosis into terms originating from either ensemble variance of isotropic diffusivity or microscopic diffusion anisotropy. This ability rests on the assumption that diffusion can be described as a sum of multiple Gaussian compartments, but this is often not strictly fulfilled. The effects of nongaussian diffusion on single shot isotropic diffusion sequences were first considered in detail by de Swiet and Mitra in 1996. They showed theoretically that anisotropic compartments lead to anisotropic time dependence of the diffusion tensors, which causes the measured isotropic diffusivity to depend on gradient frame orientation. Here we show how such deviations from the multiple Gaussian compartments assumption conflates orientation dispersion with ensemble variance in isotropic diffusivity. Second, we consider additional contributions to the apparent variance in isotropic diffusivity arising due to intracompartmental kurtosis. These will likewise depend on gradient frame orientation. We illustrate the potential importance of these confounds with analytical expressions, numerical simulations in simple model geometries, and microimaging experiments in fixed spinal cord using isotropic diffusion encoding waveforms with 7.5 ms duration and 3000 mT/m maximum amplitude.

Diffusion time dependence of microstructural parameters in fixed spinal cord

Biophysical modelling of diffusion MRI is necessary to provide specific microstructural tissue properties. However, estimating model parameters from data with limited diffusion gradient strength, such as clinical scanners, has proven unreliable due to a shallow optimization landscape. On the other hand, estimation of diffusion kurtosis (DKI) parameters is more robust, and its parameters may be connected to microstructural parameters, given an appropriate biophysical model. However, it was previously shown that this procedure still does not provide sufficient information to uniquely determine all model parameters. In particular, a parameter degeneracy related to the relative magnitude of intra-axonal and extra-axonal diffusivities remains. Here we develop a model of diffusion in white matter including axonal dispersion and demonstrate stable estimation of all model parameters from DKI in fixed pig spinal cord. By employing the recently developed fast axisymmetric DKI, we use stimulated echo acquisition mode to collect data over a two orders of magnitude diffusion time range with very narrow diffusion gradient pulses, enabling finely resolved measurements of diffusion time dependence of both net diffusion and kurtosis metrics, as well as model intra- and extra-axonal diffusivities, and axonal dispersion. Our results demonstrate substantial time dependence of all parameters except volume fractions, and the additional time dimension provides support for intra-axonal diffusivity to be larger than extra-axonal diffusivity in spinal cord white matter, although not unambiguously. We compare our findings for the time-dependent compartmental diffusivities to predictions from effective medium theory with reasonable agreement.

Oscillating gradient diffusion kurtosis imaging of normal and injured mouse brains

Recent advances in diffusion MRI employ multiple diffusion encoding schemes with varying diffusion direction, weighting, and diffusion time to investigate specific microstructural properties in biological tissues. In this study, we examined time-dependent diffusion kurtosis contrast in adult mouse brains and in neonatal mouse brains after hypoxic-ischemic (HI) injury. In vivo diffusion kurtosis maps were acquired with a short diffusion time using an oscillating gradient spin echo (OGSE) sequence at 100 Hz and with a relatively long diffusion time (20 ms) using a pulsed gradient spin echo (PGSE) sequence. In the adult mouse brain, we found that the cortex and hippocampus showed larger differences between OGSE kurtosis and PGSE kurtosis than major white matter tracts. In neonatal mouse brains with unilateral HI injury, the OGSE kurtosis map overall provided stronger edema contrast than the PGSE kurtosis map, and the differences between OGSE and PGSE kurtosis measurements in the edema region reflected heterogeneity of injury. This is the first in vivo study that has demonstrated multi-direction OGSE kurtosis contrasts in the mouse brain. Comparing PGSE and OGSE kurtosis measures may provide additional information on microstructural changes after ischemic stroke.

Magnetic resonance temporal diffusion tensor spectroscopy of disordered anisotropic tissue

Molecular diffusion measured with diffusion weighted MRI (DWI) offers a probe for tissue microstructure. However, inferring microstructural properties from conventional DWI data is a complex inverse problem and has to account for heterogeneity in sizes, shapes and orientations of the tissue compartments contained within an imaging voxel. Alternative experimental means for disentangling the signal signatures of such features could provide a stronger link between the data and its interpretation. Double diffusion encoding (DDE) offers the possibility to factor out variation in compartment shapes from orientational dispersion of anisotropic domains by measuring the correlation between diffusivity in multiple directions. Time dependence of the diffusion is another effect reflecting the dimensions and distributions of barriers. In this paper we extend on DDE with a modified version of the oscillating gradient spin echo (OGSE) experiment, giving a basic contrast mechanism closely linked to both the temporal diffusion spectrum and the compartment anisotropy. We demonstrate our new method on post mortem brain tissue and show that we retrieve the correct temporal diffusion tensor spectrum in synthetic data from Monte Carlo simulations of random walks in a range of disordered geometries of different sizes and shapes.

Brain Metabolite Diffusion from Ultra-Short to Ultra-Long Time Scales: What Do We Learn, Where Should We Go?

In vivo diffusion-weighted MR spectroscopy (DW-MRS) allows measuring diffusion properties of brain metabolites. Unlike water, most metabolites are confined within cells. Hence, their diffusion is expected to purely reflect intracellular properties, opening unique possibilities to use metabolites as specific probes to explore cellular organization and structure. However, interpretation and modeling of DW-MRS, and more generally of intracellular diffusion, remains difficult. In this perspective paper, we will focus on the study of the time-dependency of brain metabolite apparent diffusion coefficient (ADC). We will see how measuring ADC over several orders of magnitude of diffusion times, from less than 1 ms to more than 1 s, allows clarifying our understanding of brain metabolite diffusion, by firmly establishing that metabolites are neither massively transported by active mechanisms nor massively confined in subcellular compartments or cell bodies. Metabolites appear to be instead diffusing in long fibers typical of neurons and glial cells such as astrocytes. Furthermore, we will evoke modeling of ADC time-dependency to evaluate the effect of, and possibly quantify, some structural parameters at various spatial scales, departing from a simple model of hollow cylinders and introducing additional complexity, either short-ranged (such as dendritic spines) or long-ranged (such as cellular fibers ramification). Finally, we will discuss the experimental feasibility and expected benefits of extending the range of diffusion times toward even shorter and longer values.

Spatial Restriction within Intracranial Epidermoid Cysts Observed Using Short Diffusion-time Diffusion-weighted Imaging

We report two cases of pathologically proven intracranial epidermoid cysts. Both cases were scanned with diffusion-weighted imaging using pulsed gradient spin-echo (PGSE) and oscillating gradient spin-echo (OGSE; 50 Hz) prototype sequences with diffusion times of 47.3 ms and 8.5 ms, respectively. The apparent diffusion coefficient measured by OGSE was higher than that measured by PGSE, indicating the spatial restriction of water diffusion in the laminated keratin layers within the cyst as demonstrated by histopathology.