Water diffusion compartmentation and anisotropy at high b values in the human brain

The clinical and pathological evaluation of patients with immunoglobulin A nephropathy by diffusion tensor imaging and intravoxel incoherent motion diffusion-weighted imaging

OBJECTIVES

To explore the efficacy of diffuse magnetic resonance imaging (MRI) for identifying clinicopathological changes in immunoglobulin A nephropathy (IgAN) patients.

METHODS

The study enrolled IgAN patients and healthy volunteers. IgAN patients were divided into Group 1 [estimated glomerular filtration rate (eGFR) ≥ 90 mL/min/1.73 m2], Group 2 (60 ≤ eGFR < 90 mL/min/1.73 m2), and Group 3 (eGFR < 60 mL/min/1.73 m2). Intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) and diffusion tensor imaging (DTI) were performed via 3.0 T magnetic resonance. Diffuse MRI, clinical, and pathological indicators were collected and analysed. P < .05 was considered statistically significant.

RESULTS

Forty-six IgAN patients and twenty-seven volunteers were enrolled. The apparent diffusion coefficient, diffusion coefficient (D), perfusion fraction (f), and fractional anisotropy (FA) were significantly different among IgAN subgroups and controls. These parameters were positively correlated with eGFR and negatively with creatinine, and inversely correlated with glomerular sclerosis, interstitial fibrosis, and tubular atrophy (all P < .05). They had significantly high area under the curve (AUC) for distinguishing IgAN patients from controls, while FA had the highest AUC in identifying Group 1 IgAN patients from volunteers.

CONCLUSIONS

DTI and IVIM-DWI had the advantage of evaluating clinical and pathological changes in IgAN patients. DTI was superior at distinguishing early IgAN patients and might be a noninvasive marker for screening early IgAN patients from healthy individuals.

ADVANCES IN KNOWLEDGE

DTI and IVIM-DWI could evaluate clinical and pathological changes and correlated with Oxford classification in IgAN patients. They could also identify IgAN patients from healthy populations, while DTI had superiority in differentiating early IgAN patients.

Diffusion tensor imaging metrics as natural markers of multiple sclerosis-induced brain disorders with a low Expanded Disability Status Scale score

Non-invasive and effective differentiation along with determining the degree of deviations compared to the healthy cohort is important in the case of various brain disorders, including multiple sclerosis (MS). Evaluation of the effectiveness of diffusion tensor metrics (DTM) in 3T DTI for recording MS-related deviations was performed using a time-acceptable MRI protocol with unique comprehensive detection of systematic errors related to spatial heterogeneity of magnetic field gradients. In a clinical study, DTMs were acquired in segmented regions of interest (ROIs) for 50 randomly selected healthy controls (HC) and 50 multiple sclerosis patients. Identical phantom imaging was performed for each clinical measurement to estimate and remove the influence of systematic errors using the b-matrix spatial distribution in the DTI (BSD-DTI) technique. In the absence of statistically significant differences due to age in healthy volunteers and patients with multiple sclerosis, the existence of significant differences between groups was proven using DTM. Moreover, a statistically significant impact of spatial systematic errors occurs for all ROIs and DTMs in the phantom and for approximately 90 % in the HC and MS groups. In the case of a single patient measurement, this appears for all the examined ROIs and DTMs. The obtained DTMs effectively discriminate healthy volunteers from multiple sclerosis patients with a low mean score on the Expanded Disability Status Scale. The magnitude of the group differences is typically significant, with an effect size of approximately 0.5, and similar in both the standard approach and after elimination of systematic errors. Differences were also observed between metrics obtained using these two approaches. Despite a small alterations in mean DTMs values for groups and ROIs (1-3 %), these differences were characterized by a huge effect (effect size ∼0.8 or more). These findings indicate the importance of determining the spatial distribution of systematic errors specific to each MR scanner and DTI acquisition protocol in order to assess their impact on DTM in the ROIs examined. This is crucial to establish accurate DTM values for both individual patients and mean values for a healthy population as a reference. This approach allows for an initial reliable diagnosis based on DTI metrics.

Romer-EPTI: rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale dMRI and microstructure imaging

Purpose

To overcome the major challenges in dMRI acquisition, including low SNR, distortion/blurring, and motion vulnerability.

Methods

A novel Romer-EPTI technique is developed to provide distortion-free dMRI with significant SNR gain, high motion-robustness, sharp spatial resolution, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free EPTI encoding. Romer enhances SNR by a simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness through a motion-aware super-resolution reconstruction, which also incorporates slice-profile and real-value diffusion, to resolve high-isotropic-resolution volumes. The in-plane encoding is performed using distortion/blurring-free EPTI, which further improves effective spatial resolution and motion robustness by preventing not only T2/T2*-blurring but also additional blurring resulting from combining encoded volumes with inconsistent geometries caused by dynamic distortions. Self-navigation was incorporated to enable efficient phase correction. Additional developments include strategies to address slab-boundary artifacts, achieve minimal TE for SNR gain at 7T, and achieve high robustness to strong phase variations at high b-values.

Results

Using Romer-EPTI, we demonstrate distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm-iso) and 7T (485-μm-iso) for the first time, with high SNR efficiency (e.g., 25 × ), and high image quality free from distortion and slab-boundary artifacts with minimal blurring. Motion experiments demonstrate Romer-EPTI's high motion-robustness and ability to recover sharp images in the presence of motion. Romer-EPTI also demonstrates significant SNR gain and robustness in high b-value (b=5000s/mm2) and time-dependent dMRI.

Conclusion

Romer-EPTI significantly improves SNR, motion-robustness, and image quality, providing a highly efficient acquisition for high-resolution dMRI and microstructure imaging.

Probing brain tissue microstructure with MRI: principles, challenges, and the role of multidimensional diffusion-relaxation encoding

Diffusion MRI uses the random displacement of water molecules to sensitize the signal to brain microstructure and to properties such as the density and shape of cells. Microstructure modeling techniques aim to estimate these properties from acquired data by separating the signal between virtual tissue 'compartments' such as the intra-neurite and the extra-cellular space. A key challenge is that the diffusion MRI signal is relatively featureless compared with the complexity of brain tissue. Another challenge is that the tissue microstructure is wildly different within the gray and white matter of the brain. In this review, we use results from multidimensional diffusion encoding techniques to discuss these challenges and their tentative solutions. Multidimensional encoding increases the information content of the data by varying not only the b-value and the encoding direction but also additional experimental parameters such as the shape of the b-tensor and the echo time. Three main insights have emerged from such encoding. First, multidimensional data contradict common model assumptions on diffusion and T2 relaxation, and illustrates how the use of these assumptions cause erroneous interpretations in both healthy brain and pathology. Second, many model assumptions can be dispensed with if data are acquired with multidimensional encoding. The necessary data can be easily acquired in vivo using protocols optimized to minimize Cramér-Rao lower bounds. Third, microscopic diffusion anisotropy reflects the presence of axons but not dendrites. This insight stands in contrast to current 'neurite models' of brain tissue, which assume that axons in white matter and dendrites in gray matter feature highly similar diffusion. Nevertheless, as an axon-based contrast, microscopic anisotropy can differentiate gray and white matter when myelin alterations confound conventional MRI contrasts.

Monte Carlo simulation of water diffusion through cardiac tissue models

Monte Carlo diffusion simulations are commonly used to establish a reliable ground truth of tissue microstructure, including for the validation of diffusion-weighted MRI. However, selecting simulation parameters is challenging and affects validity and reproducibility. We conducted experiments to investigate critical conditions in Monte Carlo simulations, such as tissue representation complexity, simulated molecules, update duration, and compartment size. Results show significant changes in microstructure characteristics when parameters are altered, emphasizing the importance of careful control for a reliable ground truth.

In Vivo Myelin Water Quantification Using Diffusion-Relaxation Correlation MRI: A Comparison of 1D and 2D Methods

Multidimensional Magnetic Resonance Imaging (MRI) is a versatile tool for microstructure mapping. We use a diffusion weighted inversion recovery spin echo (DW-IR-SE) sequence with spiral readouts at ultra-strong gradients to acquire a rich diffusion-relaxation data set with sensitivity to myelin water. We reconstruct 1D and 2D spectra with a two-step convex optimization approach and investigate a variety of multidimensional MRI methods, including 1D multi-component relaxometry, 1D multi-component diffusometry, 2D relaxation correlation imaging, and 2D diffusion-relaxation correlation spectroscopic imaging (DR-CSI), in terms of their potential to quantify tissue microstructure, including the myelin water fraction (MWF). We observe a distinct spectral peak that we attribute to myelin water in multi-component T1 relaxometry, T1-T2 correlation, T1-D correlation, and T2-D correlation imaging. Due to lower achievable echo times compared to diffusometry, MWF maps from relaxometry have higher quality. Whilst 1D multi-component T1 data allows much faster myelin mapping, 2D approaches could offer unique insights into tissue microstructure and especially myelin diffusion.

Anesthetic modulation of water diffusion: Insights from a diffusion tensor imaging study

Diffusion tensor imaging (DTI) in animal models are essential for translational neuroscience studies. A critical step in animal studies is the use of anesthetics. Understanding the influence of specific anesthesia regimes on DTI-derived parameters, such as fractional anisotropy (FA) and mean diffusivity (MD), is imperative when comparing results between animal studies using different anesthetics. Here, the quantification of FA and MD under different anesthetic regimes, alpha-chloralose and isoflurane, is discussed. We also used a range of b-values to determine whether the anesthetic effect was b-value dependent. The first group of rats (n = 6) was anesthetized with alpha-chloralose (80 mg/kg), whereas the second group of rats (n = 7) was anesthetized with isoflurane (1.5%). DTI was performed with b-values of 500, 1500, and 1500s/mm2, and the MD and FA were assessed individually. Anesthesia-specific differences in MD were apparent, as manifested by the higher estimated MD under isoflurane anesthesia than that under alpha-chloralose anesthesia (P < 0.001). MD values increased with decreasing b-value in all regions studied, and the degree of increase when rats were anesthetized with isoflurane was more pronounced than that associated with alpha-chloralose (P < 0.05). FA quantitation was also influenced by anesthesia regimens to varying extents, depending on the brain regions and b-values. In conclusion, both scanning parameters and the anesthesia regimens significantly impacted the quantification of DTI indices.

Negative emotion differentiation and white matter microstructure

BACKGROUND

Deficits in the differentiation of negative emotions - the ability to specifically identify one's negative emotions - are associated with poorer mental health outcomes. However, the processes that lead to individual differences in negative emotion differentiation are not well understood, hampering our understanding of why this process is related to poor mental health outcomes. Given that disruptions in some affective processes are associated with white matter microstructure, identifying the circuitry associated with different affective processes can inform our understanding of how disturbances in these networks may lead to psychopathology. Thus, examination of how white matter microstructure relates to individual differences in negative emotion differentiation (NED) may provide insights into (i) its component processes and (ii) its relationship to brain structure.

METHOD

The relationship between white matter microstructure and NED was examined.

RESULTS

NED was related to white matter microstructure in right anterior thalamic radiation and inferior fronto-occipital fasciculus and left peri-genual cingulum.

LIMITATIONS

Although participants self-reported psychiatric diagnoses and previous psychological treatment, psychopathology was not directly targeted, and thus, the extent to which microstructure related to NED could be examined in relation to maladaptive outcomes is limited.

CONCLUSIONS

Results indicate that NED is related to white matter microstructure and suggest that pathways subserving processes that facilitate memory, semantics, and affective experience are important for NED. Our findings provide insights into the mechanisms by which individual differences in NED arise, suggesting intervention targets that may disrupt the relationship between poor differentiation and psychopathology.

A meta-analysis of the role of diffusion tensor imaging in cervical spinal cord compression

BACKGROUND AND PURPOSE

At present, the role of diffusion tensor imaging (DTI) remains controversial. This study aimed to confirm the role of DTI by comparing the differences in fractional anisotropy (FA) values between patients with cervical spinal cord compression (CSCC) and healthy individuals.

METHODS

A systematic and comprehensive literature search was conducted using the Web of Science, Embase, PubMed, and Cochrane Library databases to compare the mean FA values of patients with CSCC and healthy controls across all compression levels in the cervical spinal cord. Essential data from the literature, such as demographic information, imaging parameters, and DTI analysis method, were extracted. Fixed- or random-effect models based on I² heterogeneity were applied to the pooled and subgroup analyses.

RESULTS

Ten studies containing 445 patients and 197 healthy volunteers were eligible. The pooled results demonstrated a decrease in mean FA values across all compression levels in the experiment group compared to those in healthy controls (standardized mean difference = -1.54; 95% confidence interval = [-1.95, -1.14]; p < .001). Meta-regression revealed that the scanner field strength and DTI analysis method had a significant effect on heterogeneity.

CONCLUSIONS

Our results show that FA values in the spinal cord decline in patients with CSCC, thus confirming the crucial role of DTI in CSCC.

Applications of advanced diffusion MRI in early brain development: a comprehensive review

Brain development follows a protracted developmental timeline with foundational processes of neurodevelopment occurring from the third trimester of gestation into the first decade of life. Defining structural maturational patterns of early brain development is a critical step in detecting divergent developmental trajectories associated with neurodevelopmental and psychiatric disorders that arise later in life. While considerable advancements have already been made in diffusion magnetic resonance imaging (dMRI) for pediatric research over the past three decades, the field of neurodevelopment is still in its infancy with remarkable scientific and clinical potential. This comprehensive review evaluates the application, findings, and limitations of advanced dMRI methods beyond diffusion tensor imaging, including diffusion kurtosis imaging (DKI), constrained spherical deconvolution (CSD), neurite orientation dispersion and density imaging (NODDI) and composite hindered and restricted model of diffusion (CHARMED) to quantify the rapid and dynamic changes supporting the underlying microstructural architectural foundations of the brain in early life.

Loss of aquaporin-4 results in glymphatic system dysfunction via brain-wide interstitial fluid stagnation

The glymphatic system is a fluid transport network of cerebrospinal fluid (CSF) entering the brain along arterial perivascular spaces, exchanging with interstitial fluid (ISF), ultimately establishing directional clearance of interstitial solutes. CSF transport is facilitated by the expression of aquaporin-4 (AQP4) water channels on the perivascular endfeet of astrocytes. Mice with genetic deletion of AQP4 (AQP4 KO) exhibit abnormalities in the brain structure and molecular water transport. Yet, no studies have systematically examined how these abnormalities in structure and water transport correlate with glymphatic function. Here, we used high-resolution 3D magnetic resonance (MR) non-contrast cisternography, diffusion-weighted MR imaging (MR-DWI) along with intravoxel-incoherent motion (IVIM) DWI, while evaluating glymphatic function using a standard dynamic contrast-enhanced MR imaging to better understand how water transport and glymphatic function is disrupted after genetic deletion of AQP4. AQP4 KO mice had larger interstitial spaces and total brain volumes resulting in higher water content and reduced CSF space volumes, despite similar CSF production rates and vascular density compared to wildtype mice. The larger interstitial fluid volume likely resulted in increased slow but not fast MR diffusion measures and coincided with reduced glymphatic influx. This markedly altered brain fluid transport in AQP4 KO mice may result from a reduction in glymphatic clearance, leading to enlargement and stagnation of fluid in the interstitial space. Overall, diffusion MR is a useful tool to evaluate glymphatic function and may serve as valuable translational biomarker to study glymphatics in human disease.

Metabolic activity diffusion imaging (MADI): I. Metabolic, cytometric modeling and simulations

Evidence mounts that the steady-state cellular water efflux (unidirectional) first-order rate constant (k_io [s^-1 ]) magnitude reflects the ongoing, cellular metabolic rate of the cytolemmal Na⁺ , K⁺ -ATPase (NKA), ^c MR_NKA (pmol [ATP consumed by NKA]/s/cell), perhaps biology's most vital enzyme. Optimal ¹ H₂ O MR k_io determinations require paramagnetic contrast agents (CAs) in model systems. However, results suggest that the homeostatic metabolic k_io biomarker magnitude in vivo is often too large to be reached with allowable or possible CA living tissue distributions. Thus, we seek a noninvasive (CA-free) method to determine k_io in vivo. Because membrane water permeability has long been considered important in tissue water diffusion, we turn to the well-known diffusion-weighted MRI (DWI) modality. To analyze the diffusion tensor magnitude, we use a parsimoniously primitive model featuring Monte Carlo simulations of water diffusion in virtual ensembles comprising water-filled and -immersed randomly sized/shaped contracted Voronoi cells. We find this requires two additional, cytometric properties: the mean cell volume (V [pL]) and the cell number density (ρ [cells/μL]), important biomarkers in their own right. We call this approach metabolic activity diffusion imaging (MADI). We simulate water molecule displacements and transverse MR signal decays covering the entirety of b-space from pure water (ρ = V = 0; k_io undefined; diffusion coefficient, D₀ ) to zero diffusion. The MADI model confirms that, in compartmented spaces with semipermeable boundaries, diffusion cannot be described as Gaussian: the nanoscopic D (D_n ) is diffusion time-dependent, a manifestation of the "diffusion dispersion". When the "well-mixed" (steady-state) condition is reached, diffusion becomes limited, mainly by the probabilities of (1) encountering (ρ, V), and (2) permeating (k_io ) cytoplasmic membranes, and less so by D_n magnitudes. Importantly, for spaces with large area/volume (A/V; claustrophobia) ratios, this can happen in less than a millisecond. The model matches literature experimental data well, with implications for DWI interpretations.

Quantitating Diffusion Enhancement in Pore Hierarchies

A microporous continuum traversed by a set of mutually perpendicular channels is considered to be a model for a hierarchically porous system of the mesoporous zeolite type. Transient profiles of molecular uptake as determined by kinetic Monte Carlo (kMC) simulation are found to be in excellent agreement with the result attained by the application of the two-region model (the Kärger model) of molecular diffusion. In particular, it is found that, in the two limiting cases referred to as fast exchange and slow exchange, there exist two simple analytical expressions for the rate of molecular uptake and hence for the quantification of transport enhancement in comparison with the purely microporous adsorbent. In the general case, transport enhancement is simply recognized by the reciprocal addition of the expressions in the two limiting cases.

Generalisation of continuous time random walk to anomalous diffusion MRI models with an age-related evaluation of human corpus callosum

Diffusion MRI measures of the human brain provide key insight into microstructural variations across individuals and into the impact of central nervous system diseases and disorders. One approach to extract information from diffusion signals has been to use biologically relevant analytical models to link millimetre scale diffusion MRI measures with microscale influences. The other approach has been to represent diffusion as an anomalous transport process and infer microstructural information from the different anomalous diffusion equation parameters. In this study, we investigated how parameters of various anomalous diffusion models vary with age in the human brain white matter, particularly focusing on the corpus callosum. We first unified several established anomalous diffusion models (the super-diffusion, sub-diffusion, quasi-diffusion and fractional Bloch-Torrey models) under the continuous time random walk modelling framework. This unification allows a consistent parameter fitting strategy to be applied from which meaningful model parameter comparisons can be made. We then provided a novel way to derive the diffusional kurtosis imaging (DKI) model, which is shown to be a degree two approximation of the sub-diffusion model. This link between the DKI and sub-diffusion models led to a new robust technique for generating maps of kurtosis and diffusivity using the sub-diffusion parameters β_SUB and D_SUB. Superior tissue contrast is achieved in kurtosis maps based on the sub-diffusion model. 7T diffusion weighted MRI data for 65 healthy participants in the age range 19-78 years was used in this study. Results revealed that anomalous diffusion model parameters α and β have shown consistent positive correlation with age in the corpus callosum, indicating α and β are sensitive to tissue microstructural changes in ageing.

Exploring arterial tissue microstructural organization using non-Gaussian diffusion magnetic resonance schemes

The purpose of this study was to characterize the alterations in microstructural organization of arterial tissue using higher-order diffusion magnetic resonance schemes. Three porcine carotid artery models namely; native, collagenase treated and decellularized, were used to estimate the contribution of collagen and smooth muscle cells (SMC) on diffusion signal attenuation using gaussian and non-gaussian schemes. The samples were imaged in a 7 T preclinical scanner. High spatial and angular resolution diffusion weighted images (DWIs) were acquired using two multi-shell (max b-value = 3000 s/mm²) acquisition protocols. The processed DWIs were fitted using monoexponential, stretched-exponential, kurtosis and bi-exponential schemes. Directionally variant and invariant microstructural parametric maps of the three artery models were obtained from the diffusion schemes. The parametric maps were used to assess the sensitivity of each diffusion scheme to collagen and SMC composition in arterial microstructural environment. The inter-model comparison showed significant differences across the considered models. The bi-exponential scheme based slow diffusion compartment (Ds) was highest in the absence of collagen, compared to native and decellularized microenvironments. In intra-model comparison, kurtosis along the radial direction was the highest. Overall, the results of this study demonstrate the efficacy of higher order dMRI schemes in mapping constituent specific alterations in arterial microstructure.

The relationship between diffusion heterogeneity and microstructural changes in high-grade gliomas using Monte Carlo simulations

PURPOSE

Diffusion-weighted imaging (DWI) may aid accurate tumor grading. Decreased diffusivity and increased diffusion heterogeneity measures have been observed in high-grade gliomas using the non-monoexponential models for DWI. However, DWI measures concerning tissue characteristics in terms of pathophysiological and structural changes are yet to be established. Thus, this study aims to investigate the relationship between the diffusion measurements and microstructural changes in the presence of high-grade gliomas using a three-dimensional Monte Carlo simulation with systematic changes of microstructural parameters.

METHODS

Water diffusion was simulated in a microenvironment along with changes associated with the presence of high-grade gliomas, including increases in cell density, nuclear volume, extracellular volume (VF_ex), and extracellular tortuosity (λ_ex), and changes in membrane permeability (P_mem). DWI signals were simulated using a pulsed gradient spin-echo sequence. The sequence parameters, including the maximum gradient strength and diffusion time, were set to be comparable to those of clinical scanners and advanced human MRI systems. The DWI signals were fitted using the gamma distribution and diffusional kurtosis models with b-values up to 6000 and 2500 s/mm², respectively.

RESULTS

The diffusivity measures (apparent diffusion coefficients (ADC), D_gamma of the gamma distribution model and D_app of the diffusional kurtosis model) decreased with increases in cell density and λ_ex, and a decrease in P_mem. These diffusivity measures increased with increases in nuclear volume and VF_ex. The diffusion heterogeneity measures (σ_gamma of the gamma distribution model and K_app of the diffusional kurtosis model) increased with increases in cell density or nuclear volume at the low P_mem, and a decrease in P_mem. Increased σ_gamma was also associated with an increase in VF_ex.

CONCLUSION

Among simulated microstructural changes, only increases in cell density at low P_mem or decreases in P_mem corresponded to both the decreased diffusivity and increased diffusion heterogeneity measures. The results suggest that increases in cell density at low P_mem or decreases in P_mem may be associated with the diffusion changes observed in high-grade gliomas.

Exploring cell membrane water exchange in aquaporin-4-deficient ischemic mouse brain using diffusion-weighted MRI

BACKGROUND

Aquaporin-4 is a membrane channel protein that is highly expressed in brain astrocytes and facilitates the transport of water molecules. It has been suggested that suppression of aquaporin-4 function may be an effective treatment for reducing cellular edema after cerebral infarction. It is therefore important to develop clinically applicable measurement systems to evaluate and better understand the effects of aquaporin-4 suppression on the living body.

METHODS

Animal models of focal cerebral ischemia were created by surgically occluding the middle cerebral artery of wild-type and aquaporin-4 knockout mice, after which multi-b-value multi-diffusion-time diffusion-weighted imaging measurements were performed. Data were analyzed with both the apparent diffusion coefficient (ADC) model and a compartmental water-exchange model.

RESULTS

ADCs were estimated for five different b value ranges. The ADC of aquaporin-4 knockout mice in the contralateral region was significantly higher than that of wild-type mice for each range. In contrast, aquaporin-4 knockout mice had significantly lower ADC than wild-type mice in ischemic tissue for each b-value range. Genotype-dependent differences in the ADC were particularly significant for the lowest ranges in normal tissue and for the highest ranges in ischemic tissue. The ADCs measured at different diffusion times were significantly different for both genotypes. Fitting of the water-exchange model to the ischemic region data found that the water-exchange time in aquaporin-4 knockout mice was approximately 2.5 times longer than that in wild-type mice.

CONCLUSIONS

Multi-b-value multi-diffusion-time diffusion-weighted imaging may be useful for in vivo research and clinical diagnosis of aquaporin-4-related diseases.

Characterisation of microvessel blood velocity and segment length in the brain using multi-diffusion-time diffusion-weighted MRI

Multi-diffusion-time diffusion-weighted MRI can probe tissue microstructure, but the method has not been widely applied to the microvasculature. At long diffusion-times, blood flow in capillaries is in the diffusive regime, and signal attenuation is dependent on blood velocity (v) and capillary segment length (l). It is described by the pseudo-diffusion coefficient (D*=vl/6) of intravoxel incoherent motion (IVIM). At shorter diffusion-times, blood flow is in the ballistic regime, and signal attenuation depends on v, and not l. In theory, l could be estimated using D* and v. In this study, we compare the accuracy and repeatability of three approaches to estimating v, and therefore l: the IVIM ballistic model, the velocity autocorrelation model, and the ballistic approximation to the velocity autocorrelation model. Twenty-nine rat datasets from two strains were acquired at 7 T, with b-values between 0 and 1000 smm^-2 and diffusion times between 11.6 and 50 ms. Five rats were scanned twice to assess scan-rescan repeatability. Measurements of l were validated using corrosion casting and micro-CT imaging. The ballistic approximation of the velocity autocorrelation model had lowest bias relative to corrosion cast estimates of l, and had highest repeatability.

On the potential for mapping apparent neural soma density via a clinically viable diffusion MRI protocol

•

B-tensor encoding enables estimation of spherical cellular structures in the brain.

•

Spherical compartments may provide markers for apparent neural soma density.

•

Model parameters can be estimated in a fast and robust way using deep learning.

•

Practical acquisition times are achievable on widely available clinical scanners.

Diffusion MRI is a valuable tool for probing tissue microstructure in the brain noninvasively. Today, model-based techniques are widely available and used for white matter characterisation where their development is relatively mature. Conversely, tissue modelling in grey matter is more challenging, and no generally accepted models exist. With advances in measurement technology and modelling efforts, a clinically viable technique that reveals salient features of grey matter microstructure, such as the density of quasi-spherical cell bodies and quasi-cylindrical cell projections, is an exciting prospect. As a step towards capturing the microscopic architecture of grey matter in clinically feasible settings, this work uses a biophysical model that is designed to disentangle the diffusion signatures of spherical and cylindrical structures in the presence of orientation heterogeneity, and takes advantage of B-tensor encoding measurements, which provide additional sensitivity compared to standard single diffusion encoding sequences. For the fast and robust estimation of microstructural parameters, we leverage recent advances in machine learning and replace conventional fitting techniques with an artificial neural network that fits complex biophysical models within seconds. Our results demonstrate apparent markers of spherical and cylindrical geometries in healthy human subjects, and in particular an increased volume fraction of spherical compartments in grey matter compared to white matter. We evaluate the extent to which spherical and cylindrical geometries may be interpreted as correlates of neural soma and neural projections, respectively, and quantify parameter estimation errors in the presence of various departures from the modelling assumptions. While further work is necessary to translate the ideas presented in this work to the clinic, we suggest that biomarkers focussing on quasi-spherical cellular geometries may be valuable for the enhanced assessment of neurodevelopmental disorders and neurodegenerative diseases.

Statistical Evaluation of Different Mathematical Models for Diffusion Weighted Imaging of Prostate Cancer Xenografts in Mice

Purpose

To evaluate fitting quality and repeatability of four mathematical models for diffusion weighted imaging (DWI) during tumor progression in mouse xenograft model of prostate cancer.

Methods

Human prostate cancer cells (PC-3) were implanted subcutaneously in right hind limbs of 11 immunodeficient mice. Tumor growth was followed by weekly DWI examinations using a 7T MR scanner. Additional DWI examination was performed after repositioning following the fourth DWI examination to evaluate short term repeatability. DWI was performed using 15 and 12 b-values in the ranges of 0-500 and 0-2000 s/mm², respectively. Corrected Akaike information criteria and F-ratio were used to evaluate fitting quality of each model (mono-exponential, stretched exponential, kurtosis, and bi-exponential).

Results

Significant changes were observed in DWI data during the tumor growth, indicated by ADC_m, ADC_s, and ADC_k. Similar results were obtained using low as well as high b-values. No marked changes in model preference were present between the weeks 1−4. The parameters of the mono-exponential, stretched exponential, and kurtosis models had smaller confidence interval and coefficient of repeatability values than the parameters of the bi-exponential model.

Conclusion

Stretched exponential and kurtosis models showed better fit to DWI data than the mono-exponential model and presented with good repeatability.

Diffusion-weighted MRI of ischemic stroke at 3T: Value of synthetic b-values

OBJECTIVES

Diffusion-weighted imaging (DWI) plays a crucial role in the diagnosis of ischemic stroke. We assessed the value of computed and acquired high b-value DWI in comparison with conventional b = 1000 s mm^-2 DWI for ischemic stroke at 3T.

METHODS

We included 36 patients with acute ischemic stroke who presented with diffusion abnormalities on DWI performed within 24 h of symptom onset. B-values of 0, 500, 1000 and 2000 s mm^-2 were acquired. Synthetic images with b-values of 1000, 1500, 2000 and 2500 s mm^-2 were computed. Two readers compared synthetic (syn) and acquired (acq) b = 2000 s mm^-2 images with acquired b = 1000 s mm^-2 images in terms of lesion detection rate, image quality, presence of uncertain hyperintensities and lesion conspicuity. Readers also selected their preferred b-value. Contrast ratio (CR) measurements were performed. Non-parametrical statistical tests and weighted Cohens' κ tests were computed.

RESULTS

Syn1000 and syn1500 matched acq1000 images in terms of lesion detection rate, image quality and presence of uncertain hyperintensities but presented with significantly improved lesion conspicuity (p < 0.01) and were frequently selected as preferred b-values. Acq2000 images exhibited a similar lesion detection rate and improved lesion conspicuity (p < 0.01) but worse image quality (p < 0.01) than acq1000 images. Syn2000 and syn2500 images performed significantly worse (p < 0.01) than acq1000 images in most or all categories. CR significantly increased with increasing b-values.

CONCLUSION

Synthetic images at b = 1000 and 1500 s mm^-2 and acquired DWI images at b = 2000 s mm^-2 may be of clinical value due to improved lesion conspicuity.

ADVANCES IN KNOWLEDGE

Synthetic b-values enable improved lesion conspicuity for DWI of ischemic stroke.

Transient changes in white matter microstructure during general anesthesia

Cognitive dysfunction after surgery under general anesthesia is a well-recognized clinical phenomenon in the elderly. Physiological effects of various anesthetic agents have been studied at length. Very little is known about potential effects of anesthesia on brain structure. In this study we used Diffusion Tensor Imaging to compare the white matter microstructure of healthy control subjects under sevoflurane anesthesia with their awake state. Fractional Anisotropy, a white mater integrity index, transiently decreases throughout the brain during sevoflurane anesthesia and then returns back to baseline. Other DTI metrics such as mean diffusivity, axial diffusivity and radial diffusivity were increased under sevoflurane anesthesia. Although DTI metrics are age dependent, the transient changes due to sevoflurane were independent of age and sex. Volumetric analysis shows various white matter volumes decreased whereas some gray matter volumes increased during sevoflurane anesthesia. These results suggest that sevoflurane anesthesia has a significant, but transient, effect on white matter microstructure. In spite of the transient effects of sevoflurane anesthesia there were no measurable effects on brain white matter as determined by the DTI metrics at 2 days and 7 days following anesthesia. The role of white matter in the loss of consciousness under anesthesia will need to be studied and MRI studies with subjects under anesthesia will need to take these results into account.

Interpolation and Averaging of Diffusion MRI Multi-Compartment Models

Multi-compartment models (MCM) are increasingly used to characterize the brain white matter microstructure from diffusion-weighted imaging (DWI). Their use in clinical studies is however limited by the inability to resample an MCM image towards a common reference frame, or to construct atlases from such brain microstructure models. We propose to solve this problem by first identifying that these two tasks amount to the same problem. We propose to tackle it by viewing it as a simplification problem, solved thanks to spectral clustering and the definition of semi-metrics between several usual compartments encountered in the MCM literature. This generic framework is evaluated for two models: the multi-tensor model where individual fibers are modeled as individual tensors and the diffusion direction imaging (DDI) model that differentiates intra- and extra-axonal components of each fiber. Results on simulated data, simulated transformations and real data show the ability of our method to well interpolate MCM images of these types. We finally present as an application an MCM template of normal controls constructed using our approach.

The localization regime in a nutshell

High diffusion-sensitizing magnetic field gradients have been more and more often applied nowadays to achieve a better characterization of the microstructure. As the resulting spin-echo signal significantly deviates from the conventional Gaussian form, various models have been employed to interpret these deviations and to relate them with the microstructural properties of a sample. In this paper, we argue that the non-Gaussian behavior of the signal is a generic universal feature of the Bloch-Torrey equation. We provide a simple yet rigorous description of the localization regime emerging at high extended gradients and identify its origin as a symmetry breaking at the reflecting boundary. We compare the consequent non-Gaussian signal decay to other diffusion NMR regimes such as slow-diffusion, motional-narrowing and diffusion-diffraction regimes. We emphasize limitations of conventional perturbative techniques and advocate for non-perturbative approaches which may pave a way to new imaging modalities in this field.

Higher b-values improve the correlation between diffusion MRI and the cortical microarchitecture

PURPOSE

In diffusion MRI (dMRI), it remains unclear to know how much increase of b-value is conveying additional biological meaning. We tested the correlations between cortical microarchitecture and diffusion metrics computed from standard (1000 s/mm²), high (3000 s/mm²), to very high (5000 s/mm²) b-value dMRI.

METHODS

Healthy volunteers were scanned with a dMRI pulse sequence that was first optimized together with a T1-WI and T2-WI. Averaged cortical surface map of estimated myelin (T1-WI/T2-WI) was compared with surface maps of mean diffusivity (MD) computed from each b-value (MD₁₀₀₀, MD₃₀₀₀, and MD₅₀₀₀) and to surface map of mean kurtosis (MK computed from the 0-, 1000-, to 3000-s/mm² shells) in 360 cortical parcels using Spearman correlations, multiple linear regressions, and Akaike information criteria (AIC).

RESULTS

Surface map from MD₁₀₀₀ showed variations not related to myelin but the MD₃₀₀₀ and MD₅₀₀₀ maps inversely mirrored estimated myelin map; lower MD values being observed in more myelinated cortical areas. MK mirrored myelinated cortical areas. Quantitatively, Spearman correlations between myelin and MD became more and more negative as long as b-values increased while the correlation was positive between myelin and MK. Multiple regression models confirmed negative associations between myelin and MD that were significantly better from MD₁₀₀₀ to MD₃₀₀₀ and MD₅₀₀₀ (R² = 0.33, p < 0.001; R² = 0.43, p < 0.001; and R² = 0.50, p < 0.001) and positive association between myelin and MK (R² = 0.53, p < 0.001). Comparisons of the 3 statistical models showed the best performances with MK and MD₅₀₀₀ (AIC_MK < AIC_MD5000 < AIC_MD3000 < AIC_MD1000).

CONCLUSION

Higher b-values are more closely related to subtle cellular variations of the cortical microarchitecture.

Quantitative analysis of intervertebral disc degeneration using Q-space imaging in a rat model

The degree of intervertebral disc (IVD) degeneration is qualitatively evaluated on T2-weighted imaging (T2WI). However, it is difficult to assess subtle changes in IVD degeneration using T2WI. Q-space imaging (QSI) is a quantitative diffusion-weighted imaging modality used to detect subtle changes in microenvironments. This study aimed to evaluate whether QSI can detect the inhibitory effects of the antioxidant N-acetylcysteine (NAC) in IVD degeneration. We classified female Wistar rats into control, puncture, and NAC groups (n = 5 per group). In the puncture and NAC groups, IVDs were punctured using a needle. The antioxidant NAC, which suppresses the progression of IVD degeneration, was orally administered in the NAC group 1 week prior to puncture. The progression and inhibitory effect of NAC in IVD degeneration were assessed using magnetic resonance imaging (MRI): IVD height, T2 mapping, apparent diffusion coefficient (ADC), and QSI. MRI was performed using a 7-Tesla system with a conventional probe (20 IVDs in each group). QSI parameters that were assessed included Kurtosis, the probability at zero displacement (ZDP), and full width at half maximum (FWHM). IVD degeneration by puncture was confirmed by histology, IVD height, T2 mapping, ADC, and all QSI parameters (P < .001); however, the inhibitory effect of NAC was confirmed only by QSI parameters (Kurtosis and ZDP: both P < .001; FWHM: P < .01). Kurtosis had the largest effect size (Kurtosis: 1.13, ZDP: 1.06, and FWHM: 1.02) when puncture and NAC groups were compared. QSI has a higher sensitivity than conventional quantitative methods for detecting the progressive change and inhibitory effect of NAC in IVD degeneration.

Efficiency of High and Standard b Value Diffusion-Weighted Magnetic Resonance Imaging in Grading of Gliomas

Background

Glioma is the most common fatal malignant tumor of the CNS. Early detection of glioma grades based on diffusion-weighted imaging (DWI) properties is considered one of the most recent noninvasive promising tools in the assessment of glioma grade and could be helpful in monitoring patient prognosis and response to therapy.

Aim

This study aimed to investigate the accuracy of DWI at both standard and high b values (b = 1000 s/mm² and b = 3000 s/mm²) to distinguish high-grade glioma (HGG) from low-grade glioma (LGG) in clinical practice based on histopathological results.

Materials and Methods

Twenty-three patients with glioma had DWI at l.5 T MR using two different b values (b = 1000 s/mm² and b = 3000 s/mm²) at Al-Shifa Medical Complex after obtaining ethical and administrative approvals, and data were collected from March 2019 to March 2020. Minimum, maximum, and mean of apparent diffusion coefficient (ADC) values were measured through drawing region of interest (ROI) on a solid part at ADC maps. Data were analyzed by using the MedCalc analysis program, version 19.0.4, receiver operating characteristic (ROC) curve analysis was done, and optimal cutoff values for grading gliomas were determined. Sensitivity and specificity were also calculated.

Results

The obtained results showed the ADC_mean, ADC_ratio, ADC_max, and ADC_min were performed to differentiate between LGG and HGG at both standard and high b values. Moreover, ADC values were inversely proportional to glioma grade, and these differences are more obvious at high b value. Minimum ADC values using standard b value were 1.13 ± 0.17 × 10⁻³ mm²/s, 0.89 ± 0.85 × 10⁻³ mm²/s, and 0.82 ± 0.17 × 10⁻³ mm²/s for grades II, III, and IV, respectively. Concerning high b value, ADC_min values were 0.76 ± 0.07 × 10⁻³ mm²/s, 0.61 ± 0.01 × 10⁻³ mm²/s, and 0.48 ± 0.07 × 10⁻³ mm²/s for grades II, III, and IV, respectively. ADC values were inversely correlated with results of glioma grades, and the correlation was stronger at ADC₃₀₀₀ (r = −0.722, P ≤ 0.001). The ADC₃₀₀₀ achieved the highest diagnostic accuracy with an area under the curve (AUC) of 0.618, 100% sensitivity, 85.7% specificity, and 85.7% accuracy for glioma grading at a cutoff point of ≤0.618 × 10⁻³ mm²/s. The high b value showed stronger agreement with histopathology compared with standard b value results (k = 0.89 and 0.79), respectively.

Conclusion

The ADC values decrease with an increase in tumor cellularity. Meanwhile, high b value provides better tissue contrast by reflecting more tissue diffusivity. Therefore, ADC-derived parameters at high b value are more useful in the grading of glioma than those obtained at standard b value. They might be a better surrogate imaging sequence in the preoperative evaluation of gliomas.

Analysis of Apparent Diffusion Coefficients of the Brain in Healthy Controls: A Comparison Study between Single-Shot Echo-Planar Imaging and Read-out-Segmented Echo-Planar Imaging

OBJECTIVE

To compare apparent diffusion coefficients (ADCs) of brain segments by using two diffusion-weighted imaging acquisition modes, single-shot echo-planar imaging (ss-EPI) and read-out-segmented echo-planar imaging (rs-EPI), and to assess their correlation and agreement in healthy controls.

MATERIALS AND METHODS

T2-weighted (T2W) images, rs-EPI, and ss-EPI of 30 healthy subjects were acquired using a 3T magnetic resonance scanner. The T2W images were co-registered to the rs-EPI and ss-EPI, which were then segmented into the gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) to generate masking templates. ADC maps of rs-EPI and ss-EPI were also segmented into the GM, WM, and CSF by using the generated templates. ADCs of rs-EPI and ss-EPI were compared using Student's t tests and correlated using Pearson's correlation coefficients. Bland-Altman plots were used to assess the agreement between acquisitions.

RESULTS

ADCs of rs-EPI and ss-EPI were significantly different in the GM (p < 0.001) and WM (p < 0.001). ADCs showed high agreement and correlation in the whole brain and CSF (r > 0.988; p < 0.001). ADC of the WM showed the least correlation (r = 0.894; p < 0.001), and ADCs of the WM and GM showed poor agreement. Pearson's correlation equations for each brain segment were y = 1.1x - 59.4 (GM), y = 1.45x - 255 (WM), and y = 0.98x - 63.5 (CSF), where x and y indicated ADCs of rs-EPI and ss-EPI, respectively.

CONCLUSION

While ADCs of rs-EPI and ss-EPI showed high correlation and agreement in the whole brain and CSF, ADCs of the WM and GM showed significant differences and large variability, reflecting brain parenchymal inhomogeneity due to different regional microenvironments. ADCs of different acquisition methods should be interpreted carefully, especially in intra-individual comparisons.

Feasibility of IVIM parameters from diffusion-weighted imaging at 11.7T MRI for detecting ischemic changes in common carotid artery occlusion rats

This study aimed to investigate whether intravoxel incoherent motion (IVIM) parameters can identify ischemic changes in the rat cerebral cortex using a preclinical ultra-high-field 11.7 Tesla magnetic resonance imaging (11.7TMRI) scanner. In nine female Wistar rats (eight weeks old), diffusion-weighted imaging (DWI) for IVIM analysis was successfully performed before (Pre) and after unilateral (UCCAO) and bilateral (BCCAO) common carotid artery occlusion. From the acquired DWI signals averaged in six regions of interest (ROI) placed on the cortex, volume fraction of perfusion compartment (F), pseudo diffusion coefficient (D*), F × D* and apparent diffusion coefficient (ADC) were determined as IVIM parameters in the following three DWI signal models: the bi-exponential, kurtosis, and tri-exponential model. For a subgroup analysis, four rats that survived two weeks after BCCAO were assigned to the long survival (LS) group, whereas the non-LS group consisted of the remaining five animals. Each IVIM parameter change among three phases (Pre, UCCAO and BCCAO) was statistically examined in each ROI. Then, the change in each rat group was also examined for subgroup analysis. All three models were able to identify cerebral ischemic change and damage as IVIM parameter change among three phases. Furthermore, the kurtosis model could identify the parameter changes in more regions than the other two models. In the subgroup analysis with the kurtosis model, ADC in non-LS group significantly decreased between UCCAO and BCCAO but not in LS group. IVIM parameters at 11.7TMRI may help us to detect the subtle ischemic change; in particular, with the kurtosis model.

Multicomponent diffusion analysis reveals microstructural alterations in spinal cord of a mouse model of amyotrophic lateral sclerosis ex vivo

The microstructure changes associated with degeneration of spinal axons in amyotrophic lateral sclerosis (ALS) may be reflected in altered water diffusion properties, potentially detectable with diffusion-weighted (DW) MRI. Prior work revealed the classical mono-exponential model fails to precisely depict decay in DW signal at high b-values. In this study, we aim to investigate signal decay behaviors at ultra-high b-values for non-invasive assessment of spinal cord alterations in the transgenic SOD1^G93A mouse model of ALS. A multiexponential diffusion analysis using regularized non-negative least squares (rNNLS) algorithm was applied to a series of thirty DW MR images with b-values ranging from 0 to 858,022 s/mm² on ex vivo spinal cords of transgenic SOD1^G93A and age-matched control mice. We compared the distributions of measured diffusion coefficient fractions between the groups. The measured diffusion weighted signals in log-scale showed non-linear decay behaviors with increased b-values. Faster signal decays were observed with diffusion gradients applied parallel to the long axis of the spinal cord compared to when oriented in the transverse direction. Multiexponential analysis at the lumbar level in the spinal cord identified ten subintervals. A significant decrease of diffusion coefficient fractions was found in the ranges of [1.63×10⁻⁸,3.70×10⁻⁶] mm²/s (P = 0.0002) and of [6.01×10⁻⁶,4.20×10⁻⁵] mm²/s (P = 0.0388) in SOD1^G93A mice. Anisotropic diffusion signals persisted at ultra-high b-value DWIs of the mouse spinal cord and multiexponential diffusion analysis offers the potential to evaluate microstructural alterations of ALS-affected spinal cord non-invasively.

Direction-averaged diffusion-weighted MRI signal using different axisymmetric B-tensor encoding schemes

Purpose

It has been shown, theoretically and in vivo, that using the Stejskal‐Tanner pulsed‐gradient, or linear tensor encoding (LTE), and in tissue exhibiting a “stick‐like” diffusion geometry, the direction‐averaged diffusion‐weighted MRI signal at high b‐values ( 7000<b<10000s/mm2) follows a power‐law, decaying as 1/b. It has also been shown, theoretically, that for planar tensor encoding (PTE), the direction‐averaged diffusion‐weighted MRI signal decays as 1/b. We aimed to confirm this theoretical prediction in vivo. We then considered the direction‐averaged signal for arbitrary b‐tensor shapes and different tissue substrates to look for other conditions under which a power‐law exists.

Methods

We considered the signal decay for high b‐values for encoding geometries ranging from 2‐dimensional PTE, through isotropic or spherical tensor encoding to LTE. When a power‐law behavior was suggested, this was tested using in silico simulations and, when appropriate, in vivo using ultra‐strong (300 mT/m) gradients.

Results

Our in vivo results confirmed the predicted 1/b power law for PTE. Moreover, our analysis showed that using an axisymmetric b‐tensor a power‐law only exists under very specific conditions: (a) “stick‐like” tissue geometry and purely LTE or purely PTE waveforms; and (b) "pancake‐like" tissue geometry and a purely LTE waveform.

Conclusions

A complete analysis of the power‐law dependencies of the diffusion‐weighted signal at high b‐values has been performed. Only three specific forms of encoding result in a power‐law dependency, pure linear and pure PTE when the tissue geometry is “stick‐like” and pure LTE when the tissue geometry is "pancake‐like". The different exponents of these encodings could be used to provide independent validation of the presence of different tissue geometries in vivo.

Modeling the diffusion-weighted imaging signal for breast lesions in the b = 200 to 3000 s/mm2 range: quality of fit and classification accuracy for different representations

PURPOSE

To evaluate different non-Gaussian representations for the diffusion-weighted imaging (DWI) signal in the b-value range 200 to 3000 s/mm² in benign and malignant breast lesions.

METHODS

Forty-three patients diagnosed with benign (n = 18) or malignant (n = 25) tumors of the breast underwent DWI (b-values 200, 600, 1200, 1800, 2400, and 3000 s/mm² ). Six different representations were fit to the average signal from regions of interest (ROIs) at different b-value ranges. Quality of fit was assessed by the corrected Akaike information criterion (AICc), and the Friedman test was used for assessing representation ranks. The area under the curve (AUC) of receiver operating characteristic curves were used to evaluate the power of derived parameters to differentiate between malignant and benign lesions. The lesion ROI was divided in central and peripheral parts to assess potential effect of heterogeneity. Sensitivity to noise-floor correction was also evaluated.

RESULTS

The Padé exponent was ranked as the best based on AICc, whereas 3 models (kurtosis, fractional, and biexponential) achieved the highest AUC = 0.99 for lesion differentiation. The monoexponential model at b_max = 600 s/mm² already provides AUC = 0.96, with considerably shorter acquisition time and simpler analysis. Significant differences between central and peripheral parts of lesions were found in malignant lesions. The mono- and biexponential models were most stable against varying degrees of noise-floor correction.

CONCLUSION

Non-Gaussian representations are required for fitting of the DWI curve at high b-values in breast lesions. However, the added clinical value from the high b-value data for differentiation of benign and malignant lesions is not clear.

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).

Anomalous water dynamics in brain: a combined diffusion magnetic resonance imaging and neutron scattering investigation

Water diffusion is an optimal tool for investigating the architecture of brain tissue on which modern medical diagnostic imaging techniques rely. However, intrinsic tissue heterogeneity causes systematic deviations from pure free-water diffusion behaviour. To date, numerous theoretical and empirical approaches have been proposed to explain the non-Gaussian profile of this process. The aim of this work is to shed light on the physics piloting water diffusion in brain tissue at the micrometre-to-atomic scale. Combined diffusion magnetic resonance imaging and first pioneering neutron scattering experiments on bovine brain tissue have been performed in order to probe diffusion distances up to macromolecular separation. The coexistence of free-like and confined water populations in brain tissue extracted from a bovine right hemisphere has been revealed at the micrometre and atomic scale. The results are relevant for improving the modelling of the physics driving intra- and extracellular water diffusion in brain, with evident benefit for the diffusion magnetic resonance imaging technique, nowadays widely used to diagnose, at the micrometre scale, brain diseases such as ischemia and tumours.

Searching for the neurite density with diffusion MRI: Challenges for biophysical modeling

In vivo mapping of the neurite density with diffusion MRI (dMRI) is a high but challenging aim. First, it is unknown whether all neurites exhibit completely anisotropic ("stick-like") diffusion. Second, the "density" of tissue components may be confounded by non-diffusion properties such as T2 relaxation. Third, the domain of validity for the estimated parameters to serve as indices of neurite density is incompletely explored. We investigated these challenges by acquiring data with "b-tensor encoding" and multiple echo times in brain regions with low orientation coherence and in white matter lesions. Results showed that microscopic anisotropy from b-tensor data is associated with myelinated axons but not with dendrites. Furthermore, b-tensor data together with data acquired for multiple echo times showed that unbiased density estimates in white matter lesions require data-driven estimates of compartment-specific T2 values. Finally, the "stick" fractions of different biophysical models could generally not serve as neurite density indices across the healthy brain and white matter lesions, where outcomes of comparisons depended on the choice of constraints. In particular, constraining compartment-specific T2 values was ambiguous in the healthy brain and had a large impact on estimated values. In summary, estimating neurite density generally requires accounting for different diffusion and/or T2 properties between axons and dendrites. Constrained "index" parameters could be valid within limited domains that should be delineated by future studies.

Monoexponential, Biexponential, and stretched-exponential models using diffusion-weighted imaging: A quantitative differentiation of breast lesions at 3.0T

BACKGROUND

Diffusion-weighted imaging (DWI) plays an important role in the differentiation of malignant and benign breast lesions.

PURPOSE

To investigate the utility of various diffusion parameters obtained from monoexponential, biexponential, and stretched-exponential DWI models in the differential diagnosis of breast lesions.

STUDY TYPE

Prospective.

POPULATION

Sixty-one patients (age range: 25-68 years old; mean age: 46 years old) with 31 malignant lesions, 42 benign lesions, and 28 normal breast tissues diagnosed initially by clinical palpation, ultrasonography, or conventional mammography were enrolled in the study from January to September 2016.

FIELD STRENGTH

3.0T MR scanner, T₁ WI, T₂ WI, DWI (conventional and multi-b values), dynamic contrast-enhanced.

ASSESSMENT

The apparent diffusion coefficient (ADC) was calculated by monoexponential analysis. The diffusion coefficient (ADC_slow ), pseudodiffusion coefficient (ADC_fast ), and perfusion fraction (f) were calculated using the biexponential model. The distributed diffusion coefficient (DDC) and water molecular diffusion heterogeneity index (α) were obtained using a stretched-exponential model. All parameters were compared for malignant tumors, benign tumors, and normal breast tissues. A receiver operating characteristic curve was used to compare the ability of these parameters, in order to differentiate benign and malignant breast lesions.

STATISTICAL TESTS

All statistical analyses were performed using statistical software (SPSS).

RESULTS

ADC, ADC_slow , f, DDC, and α values were significantly lower in malignant tumors when compared with normal breast tissues and benign tumors (P < 0.05). However, ADC and f had higher area under the receiver operating characteristic curve (AUC) values (0.889 and 0.919, respectively).

DATA CONCLUSION

The parameters derived from the biexponential and stretched-exponential DWI could provide additional information for differentiating between benign and malignant breast tumors when compared with conventional diffusion parameters.

LEVEL OF EVIDENCE

4 Technical Efficacy: Stage 4 J. Magn. Reson. Imaging 2019;50:1461-1467.

Characterization of focal liver lesions using the stretched exponential model: comparison with monoexponential and biexponential diffusion-weighted magnetic resonance imaging

OBJECTIVE

To compare the stretched exponential model of diffusion-weighted imaging (DWI) with monoexponential and biexponential models in terms of the ability to characterize focal liver lesions (FLLs).

METHODS

This retrospective study included 180 patients with FLLs who underwent magnetic resonance imaging including DWI with nine b values at 3.0 T. The distributed diffusion coefficient (DDC) and intravoxel diffusion heterogeneity index (α) from a stretched exponential model; true diffusion coefficient (D_t), pseudo-diffusion coefficient (D_p), and perfusion fraction (f) from a biexponential model; and apparent diffusion coefficient (ADC) were calculated for each lesion. Diagnostic performances of the parameters were assessed through receiver operating characteristic (ROC) analysis. For 20 patients with treated hepatic metastases, the correlation between the DWI parameters and the percentage of tumor necrosis on pathology was evaluated using the Spearman correlation coefficient.

RESULTS

DDC had the highest area under the ROC curve (AUC, 0.905) for differentiating malignant from benign lesions, followed by D_t (0.903) and ADC (0.866), without significant differences among them (DDC vs. D_t, p = 0.946; DDC vs. ADC, p = 0.157). For distinguishing hypovascular from hypervascular lesions, and hepatocellular carcinoma from metastasis, f  had a significantly higher AUC than the other DWI parameters (p < 0.05). The α had the strongest correlation with the degree of tumor necrosis (ρ = 0.655, p = 0.002).

CONCLUSION

The DDC from stretched exponential model of DWI demonstrated excellent diagnostic performance for differentiating malignant from benign FLLs. The α is promising for evaluating the degree of necrosis in treated metastases.

KEY POINTS

• The stretched exponential DWI model is valuable for characterizing focal liver lesions. • The DDC from stretched exponential model shows excellent performance for differentiating malignant from benign focal liver lesions. • The α from stretched exponential model is promising for evaluating the degree of necrosis in hepatic metastases after chemotherapy.

Utility of Diffusion-weighted Imaging to Decrease Unnecessary Biopsies Prompted by Breast MRI: A Trial of the ECOG-ACRIN Cancer Research Group (A6702)

PURPOSE

Conventional breast MRI is highly sensitive for cancer detection but prompts some false positives. We performed a prospective, multicenter study to determine whether apparent diffusion coefficients (ADCs) from diffusion-weighted imaging (DWI) can decrease MRI false positives.Experimental Design: A total of 107 women with MRI-detected BI-RADS 3, 4, or 5 lesions were enrolled from March 2014 to April 2015. ADCs were measured both centrally and at participating sites. ROC analysis was employed to assess diagnostic performance of centrally measured ADCs and identify optimal ADC thresholds to reduce unnecessary biopsies. Lesion reference standard was based on either definitive biopsy result or at least 337 days of follow-up after the initial MRI procedure.

RESULTS

Of 107 women enrolled, 67 patients (median age 49, range 24-75 years) with 81 lesions with confirmed reference standard (28 malignant, 53 benign) and evaluable DWI were analyzed. Sixty-seven of 81 lesions were BI-RADS 4 (n = 63) or 5 (n = 4) and recommended for biopsy. Malignancies exhibited lower mean in centrally measured ADCs (mm²/s) than benign lesions [1.21 × 10^-3 vs.1.47 × 10^-3; P < 0.0001; area under ROC curve = 0.75; 95% confidence interval (CI) 0.65-0.84]. In centralized analysis, application of an ADC threshold (1.53 × 10^-3 mm²/s) lowered the biopsy rate by 20.9% (14/67; 95% CI, 11.2%-31.2%) without affecting sensitivity. Application of a more conservative threshold (1.68 × 10^-3 mm²/s) to site-measured ADCs reduced the biopsy rate by 26.2% (16/61) but missed three cancers.

CONCLUSIONS

DWI can reclassify a substantial fraction of suspicious breast MRI findings as benign and thereby decrease unnecessary biopsies. ADC thresholds identified in this trial should be validated in future phase III studies.

Detecting white matter alterations in multiple sclerosis using advanced diffusion magnetic resonance imaging

Multiple sclerosis is a neurodegenerative and inflammatory disease, a hallmark of which is demyelinating lesions in the white matter. We hypothesized that alterations in white matter microstructures can be non-invasively characterized by advanced diffusion magnetic resonance imaging. Seven diffusion metrics were extracted from hybrid diffusion imaging acquisitions via classic diffusion tensor imaging, neurite orientation dispersion and density imaging, and q-space imaging. We investigated the sensitivity of the diffusion metrics in 36 sets of regions of interest in the brain white matter of six female patients (age 52.8 ± 4.3 years) with multiple sclerosis. Each region of interest set included a conventional T2-defined lesion, a matched perilesion area, and normal-appearing white matter. Six patients with multiple sclerosis (n = 5) or clinically isolated syndrome (n = 1) at a mild to moderate disability level were recruited. The patients exhibited microstructural alterations from normal-appearing white matter transitioning to perilesion areas and lesions, consistent with decreased tissue restriction, decreased axonal density, and increased classic diffusion tensor imaging diffusivity. The findings suggest that diffusion compartment modeling and q-space analysis appeared to be sensitive for detecting subtle microstructural alterations between perilesion areas and normal-appearing white matter.

Comparison of diffusion-weighted MRI and anti-Stokes Raman scattering (CARS) measurements of the inter-compartmental exchange-time of water in expression-controlled aquaporin-4 cells

We performed multi-b and multi-diffusion-time diffusion-weighted magnetic resonance imaging on aquaporin-4-expressing (AQ) and -non-expressing (noAQ) cells, and demonstrated a clear difference between the signals from the two cell types. The data were interpreted using a two-compartment (intra and extracellular spaces) model including inter-compartmental exchange. It was also assumed that restricted diffusion of water molecules inside the cells leads to the intracellular diffusion coefficient being inversely proportional to the diffusion-time. Estimates of the water-exchange-times obtained with this model are compared to those measured using an independent optical imaging technique (coherent anti-Stokes Raman scattering imaging, CARS). For both techniques it was found that the exchange-time estimated for the noAQ cells was significantly longer than that for the AQ cells.

Optimising sampling patterns for bi-exponentially decaying signals

A recently reported method, based on the Cramér-Rao Lower Bound theory, for optimising sampling patterns for a wide range of nuclear magnetic resonance (NMR) experiments is applied to the problem of optimising sampling patterns for bi-exponentially decaying signals. Sampling patterns are optimised by minimizing the percentage error in estimating the most difficult to estimate parameter of the bi-exponential model, termed the objective function. The predictions of the method are demonstrated in application to pulsed field gradient NMR data recorded for the two-component diffusion of a binary mixture of methane/ethane in a zeolite. It is shown that the proposed method identifies an optimal sampling pattern with the predicted objective function being within 10% of that calculated from the experiment dataset. The method is used to advise on the number of sampled points and the noise level needed to resolve two-component systems characterised by a range of ratios of populations and diffusion coefficients. It is subsequently illustrated how the method can be used to reduce the experiment acquisition time while still being able to resolve a given two-component system.

The Kärger vs bi-exponential model: Theoretical insights and experimental validations

We revise three common models accounting for water exchange in pulsed-gradient spin-echo measurements: a bi-exponential model with time-dependent water fractions, the Kärger model, and a modified Kärger model designed for restricted diffusion, e.g. inside cells. The three models are compared and applied to experimental data from yeast cell suspensions. The Kärger model and the modified Kärger model yield very close results and accurately fit the data. The bi-exponential model, although less rigorous, has a natural physical interpretation and suggests a new experimental modality to estimate the water exchange time.

Diffusion MRI of cancer: From low to high b-values

Following its success in early detection of cerebral ischemia, diffusion-weighted imaging (DWI) has been increasingly used in cancer diagnosis and treatment evaluation. These applications are propelled by the rapid development of novel diffusion models to extract biologically valuable information from diffusion-weighted MR signals, and significant advances in MR hardware that has enabled image acquisition with high b-values. This article reviews recent technical developments and clinical applications in cancer imaging using DWI, with a special emphasis on high b-value diffusion models. The article is organized in four sections. First, we provide an overview of diffusion models that are relevant to cancer imaging. The model parameters are discussed in relation to three tissue properties-cellularity, vascularity, and microstructures. An emphasis is placed on characterization of microstructural heterogeneity, given its novelty and close relevance to cancer. Second, we illustrate diffusion MR clinical applications in each of the following three categories: 1) cancer detection and diagnosis; 2) cancer grading, staging, and classification; and 3) cancer treatment response prediction and evaluation. Third, we discuss several practical issues, including selection of image acquisition parameters, reproducibility and reliability, motion management, image distortion, etc., that are commonly encountered when applying DWI to cancer in clinical settings. Lastly, we highlight a few ongoing challenges and provide some possible future directions, particularly in the area of establishing standards via well-organized multicenter clinical trials to accelerate clinical translation of advanced DWI techniques to improving cancer care on a large scale. Level of Evidence: 5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:23-40.

Influence of the extracellular matrix on water mobility in subcortical gray matter

PURPOSE

Water mobility in tissues is related to the microstructure that modulates diffusion and spin relaxation. Previous work has shown that the extracellular matrix (ECM) impacts water diffusion in cartilage. To investigate if similar contributions to image contrast exist for brain, which is characterized by a substantially lower ECM content, diffusion and relaxation were studied in fixed samples from goat and human thalamus before and after enzymatic digestion of ECM compounds. Selected experiments in human corpus callosum were included for comparing subcortical gray matter and white matter.

METHODS

Digestion of matrix components was achieved by treatment with hyaluronidase. Nonlocalized pulsed field gradient measurements were performed with <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>b</mml:mi></mml:math> values between 0.6 and 18,000 s/mm² at 3T and temperatures between 0°C and 20°C, in addition to T₁ and T₂ relaxation measurements. The data were fitted to multiexponential models to account for different water compartments. After the measurements, the samples were sliced and stained for ECM-sensitive markers to verify efficient digestion.

RESULTS

Microstructural alterations associated with hyaluronan digestion did not lead to measurable effects on water diffusion or <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>T</mml:mi> <mml:mn>2</mml:mn></mml:msub> </mml:math> . However, T₁ of the main relaxographic component, attributed to intra-/extracellular water, decreased by 7%.

CONCLUSION

Investigations with very strong gradients did not reveal a detectable effect on water diffusion or <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msub><mml:mi>T</mml:mi> <mml:mn>2</mml:mn></mml:msub> </mml:math> after hyaluronan removal, indicating that the brain ECM content is too low to produce a detectable effect. The subtle alteration of T₁ upon hyaluronidase treatment might reflect a modulation of intercompartmental water exchange properties.

STROBE-A preliminary investigation of IVIM-DWI in cardiac imaging

This study aims to explore the possibility to apply intravoxel incoherent motion-magnetic resonance imaging (IVIM-MRI) in cardiac imaging.Multi-b-value diffusion-weighted imaging (DWI) sequence scanning was performed on 12 healthy volunteers. A double exponential model was adopted, and the b-value sequence was 0, 20, 60, 80, 120, 200, and 600 second/mm. The D-value, D*-value, and f-value of the anterior posterior and lateral walls of the ventricular septum were respectively measured on the short axis section of the heart, the parameters of the myocardium in different blood supply areas in each segment were recorded, and the measured data of these different segments were compared using analysis of variance.Among these 12 healthy volunteers, the D-value, D*-value, and f-value of these 72 segments were not exactly equal, the D-values of the myocardium in the 5th and 11th segment were lower than those in the 2nd, 3rd, 8th, and 9th segments, and the pairwise differences were statistically significant (P < .001). Furthermore, the difference in D-value between the 5th and 11th segments was not statistically significant (P = 1.000). The D*-value and f-value of the myocardium in the 5th and 11th segment were higher than those in the 2nd, 3rd, 8th, and 9th segments, and the pairwise differences were statistically significant (P < .001). Furthermore, the differences in D*-value and f-value between the 5th and 11th segments was not statistically significant (P = .214, .787).The intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) technique can quantitatively reflect the diffusion and blood perfusion status of the myocardium.