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

Comprehensive characterization of tumor therapeutic response with simultaneous mapping cell size, density, and transcytolemmal water exchange

Early assessment of tumor therapeutic response is an important topic in precision medicine to optimize personalized treatment regimens and reduce unnecessary toxicity, cost, and delay. Although diffusion MRI (dMRI) has shown potential to address this need, its predictive accuracy is limited, likely due to its unspecific sensitivity to overall pathological changes. In this work, we propose a new quantitative dMRI-based method dubbed EXCHANGE (MRI of water Exchange, Confined and Hindered diffusion under Arbitrary Gradient waveform Encodings) for simultaneous mapping of cell size, cell density, and transcytolemmal water exchange. Such rich microstructural information comprehensively evaluates tumor pathologies at the cellular level. Validations using numerical simulations and in vitro cell experiments confirmed that the EXCHANGE method can accurately estimate mean cell size, density, and water exchange rate constants. The results from in vivo animal experiments show the potential of EXCHANGE for monitoring tumor treatment response. Finally, the EXCHANGE method was implemented in breast cancer patients with neoadjuvant chemotherapy, demonstrating its feasibility in assessing tumor therapeutic response in clinics. In summary, a new, quantitative dMRI-based EXCHANGE method was proposed to comprehensively characterize tumor microstructural properties at the cellular level, suggesting a unique means to monitor tumor treatment response in clinical practice.

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.

Diffusion time dependency of extracellular diffusion

PURPOSE

To quantify the variations of the power-law dependences on diffusion time t or gradient frequencyf $$ f $$ of extracellular water diffusion measured by diffusion MRI (dMRI).

METHODS

Model cellular systems containing only extracellular water were used to investigate thet / f $$ t/f $$ dependence ofD ex $$ {D}_{ex} $$ , the extracellular diffusion coefficient. Computer simulations used a randomly packed tissue model with realistic intracellular volume fractions and cell sizes. DMRI measurements were performed on samples consisting of liposomes containing heavy water(D2 O, deuterium oxide) dispersed in regular water (H2 O).D ex $$ {D}_{ex} $$ was obtained over a broadt $$ t $$ range (∼1-1000 ms) and then fit power-law equationsD ex ( t ) = D const + const · t - ϑ t $$ {D}_{ex}(t)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {t}^{-{\vartheta}_t} $$ andD ex ( f ) = D const + const · f ϑ f $$ {D}_{ex}(f)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {f}^{\vartheta_f} $$ .

RESULTS

Both simulated and experimental results suggest that no single power-law adequately describes the behavior ofD ex $$ {D}_{ex} $$ over the range of diffusion times of most interest in practical dMRI. Previous theoretical predictions are accurate over only limitedt $$ t $$ ranges; for example,θ t = θ f = - 1 2 $$ {\theta}_t={\theta}_f=-\frac{1}{2} $$ is valid only for short times, whereasθ t = 1 $$ {\theta}_t=1 $$ orθ f = 3 2 $$ {\theta}_f=\frac{3}{2} $$ is valid only for long times but cannot describe other ranges simultaneously. For the specifict $$ t $$ range of 5-70 ms used in typical human dMRI measurements,θ t = θ f = 1 $$ {\theta}_t={\theta}_f=1 $$ matches the data well empirically.

CONCLUSION

The optimal power-law fit of extracellular diffusion varies with diffusion time. The dependency obtained at short or longt $$ t $$ limits cannot be applied to typical dMRI measurements in human cancer or liver. It is essential to determine the appropriate diffusion time range when modeling extracellular diffusion in dMRI-based quantitative microstructural imaging.

The Value of Various Post-Processing Modalities of Diffusion Weighted Imaging in the Detection of Multiple Sclerosis

Diffusion tensor imaging (DTI) showed its adequacy in evaluating the normal-appearing white matter (NAWM) and lesions in the brain that are difficult to evaluate with routine clinical magnetic resonance imaging (MRI) in multiple sclerosis (MS). Recently, MRI systems have been developed with regard to software and hardware, leading to different proposed diffusion analysis methods such as diffusion tensor imaging, q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, and axonal diameter measurement. These methods have the ability to better detect in vivo microstructural changes in the brain than DTI. These different analysis modalities could provide supplementary inputs for MS disease characterization and help in monitoring the disease’s progression as well as treatment efficacy. This paper reviews some of the recent diffusion MRI methods used for the assessment of MS in vivo.

Linear multi-scale modeling of diffusion MRI data: A framework for characterization of oriented structures across length scales

Diffusion-weighted magnetic resonance imaging (DW-MRI) has evolved to provide increasingly sophisticated investigations of the human brain's structural connectome in vivo. Restriction spectrum imaging (RSI) is a method that reconstructs the orientation distribution of diffusion within tissues over a range of length scales. In its original formulation, RSI represented the signal as consisting of a spectrum of Gaussian diffusion response functions. Recent technological advances have enabled the use of ultra-high b-values on human MRI scanners, providing higher sensitivity to intracellular water diffusion in the living human brain. To capture the complex diffusion time dependence of the signal within restricted water compartments, we expand upon the RSI approach to represent restricted water compartments with non-Gaussian response functions, in an extended analysis framework called linear multi-scale modeling (LMM). The LMM approach is designed to resolve length scale and orientation-specific information with greater specificity to tissue microstructure in the restricted and hindered compartments, while retaining the advantages of the RSI approach in its implementation as a linear inverse problem. Using multi-shell, multi-diffusion time DW-MRI data acquired with a state-of-the-art 3 T MRI scanner equipped with 300 mT/m gradients, we demonstrate the ability of the LMM approach to distinguish different anatomical structures in the human brain and the potential to advance mapping of the human connectome through joint estimation of the fiber orientation distributions and compartment size characteristics.

A longitudinal microstructural MRI dataset in healthy C57Bl/6 mice at 9.4 Tesla

Multimodal microstructural MRI has shown increased sensitivity and specificity to changes in various brain disease and injury models in the preclinical setting. Here, we present an in vivo longitudinal dataset, including a subset of ex vivo data, acquired as control data and to investigate microstructural changes in the healthy mouse brain. The dataset consists of structural T2-weighted imaging, magnetization transfer ratio and saturation imaging, and advanced quantitative diffusion MRI (dMRI) methods. The dMRI methods include oscillating gradient spin echo (OGSE) dMRI and microscopic anisotropy (μA) dMRI, which provide additional insight by increasing sensitivity to smaller spatial scales and disentangling fiber orientation dispersion from true microstructural changes, respectively. The technical skills required to analyze microstructural MRI data are complex and include MRI sequence development, acquisition, and computational neuroimaging expertise. Here, we share unprocessed and preprocessed data, and scalar maps of quantitative MRI metrics. We envision utility of this dataset in the microstructural MRI field to develop and test biophysical models, methods that model temporal brain dynamics, and registration and preprocessing pipelines.

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.

IMPULSED model based cytological feature estimation with U-Net: Application to human brain tumor at 3T

PURPOSE

This work introduces and validates a deep-learning-based fitting method, which can rapidly provide accurate and robust estimation of cytological features of brain tumor based on the IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) model fitting with diffusion-weighted MRI data.

METHODS

The U-Net was applied to rapidly quantify extracellular diffusion coefficient (D_ex ), cell size (d), and intracellular volume fraction (v_in ) of brain tumor. At the training stage, the image-based training data, synthesized by randomizing quantifiable microstructural parameters within specific ranges, was used to train U-Net. At the test stage, the pre-trained U-Net was applied to estimate the microstructural parameters from simulated data and the in vivo data acquired on patients at 3T. The U-Net was compared with conventional non-linear least-squares (NLLS) fitting in simulations in terms of estimation accuracy and precision.

RESULTS

Our results confirm that the proposed method yields better fidelity in simulations and is more robust to noise than the NLLS fitting. For in vivo data, the U-Net yields obvious quality improvement in parameter maps, and the estimations of all parameters are in good agreement with the NLLS fitting. Moreover, our method is several orders of magnitude faster than the NLLS fitting (from about 5 min to <1 s).

CONCLUSION

The image-based training scheme proposed herein helps to improve the quality of the estimated parameters. Our deep-learning-based fitting method can estimate the cell microstructural parameters fast and accurately.

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

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

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.

Advanced Diffusion MR Imaging for Multiple Sclerosis in the Brain and Spinal Cord

Diffusion tensor imaging (DTI) has been established its usefulness in evaluating normal-appearing white matter (NAWM) and other lesions that are difficult to evaluate with routine clinical MRI in the evaluation of the brain and spinal cord lesions in multiple sclerosis (MS), a demyelinating disease. With the recent advances in the software and hardware of MRI systems, increasingly complex and sophisticated MRI and analysis methods, such as q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, white matter tract integrity, and multiple diffusion encoding, referred to as advanced diffusion MRI, have been proposed. These are capable of capturing in vivo microstructural changes in the brain and spinal cord in normal and pathological states in greater detail than DTI.

This paper reviews the current status of recent advanced diffusion MRI for assessing MS in vivo as part of an issue celebrating two decades of magnetic resonance in medical sciences (MRMS), an official journal of the Japanese Society of Magnetic Resonance in Medicine.

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.

Test-retest reproducibility of in vivo oscillating gradient and microscopic anisotropy diffusion MRI in mice at 9.4 Tesla

BACKGROUND AND PURPOSE

Microstructure imaging with advanced diffusion MRI (dMRI) techniques have shown increased sensitivity and specificity to microstructural changes in various disease and injury models. Oscillating gradient spin echo (OGSE) dMRI, implemented by varying the oscillating gradient frequency, and microscopic anisotropy (μA) dMRI, implemented via tensor valued diffusion encoding, may provide additional insight by increasing sensitivity to smaller spatial scales and disentangling fiber orientation dispersion from true microstructural changes, respectively. The aims of this study were to characterize the test-retest reproducibility of in vivo OGSE and μA dMRI metrics in the mouse brain at 9.4 Tesla and provide estimates of required sample sizes for future investigations.

METHODS

Twelve adult C57Bl/6 mice were scanned twice (5 days apart). Each imaging session consisted of multifrequency OGSE and μA dMRI protocols. Metrics investigated included μA, linear diffusion kurtosis, isotropic diffusion kurtosis, and the diffusion dispersion rate (Λ), which explores the power-law frequency dependence of mean diffusivity. The dMRI metric maps were analyzed with mean region-of-interest (ROI) and whole brain voxel-wise analysis. Bland-Altman plots and coefficients of variation (CV) were used to assess the reproducibility of OGSE and μA metrics. Furthermore, we estimated sample sizes required to detect a variety of effect sizes.

RESULTS

Bland-Altman plots showed negligible biases between test and retest sessions. ROI-based CVs revealed high reproducibility for most metrics (CVs < 15%). Voxel-wise CV maps revealed high reproducibility for μA (CVs ~ 10%), but low reproducibility for OGSE metrics (CVs ~ 50%).

CONCLUSION

Most of the μA dMRI metrics are reproducible in both ROI-based and voxel-wise analysis, while the OGSE dMRI metrics are only reproducible in ROI-based analysis. Given feasible sample sizes (10-15), μA metrics and OGSE metrics may provide sensitivity to subtle microstructural changes (4-8%) and moderate changes (> 6%), respectively.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

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.

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.