Apparent exchange rate mapping with diffusion MRI

Compartmental anisotropy of filtered exchange imaging (FEXI) in human white matter: What is happening in FEXI?

PURPOSE

To investigate the effects of compartmental anisotropy on filtered exchange imaging (FEXI) in white matter (WM).

THEORY AND METHODS

FEXI signals were measured using multiple combinations of diffusion filter and detection directions in five healthy volunteers. Additional filters, including a trace-weighted diffusion filter with trapezoidal gradients, a spherical b-tensor encoded diffusion filter, and a T2 filter, were tested with trace-weighted diffusion detection.

RESULTS

A large range of apparent exchange rates (AXR) and both positive and negative filter efficiencies (σ) were found depending on the mutual orientation of the filter and detection gradients relative to WM fiber orientation. The data demonstrated that the fast-diffusion compartment suppressed by diffusional filtering is not exclusively extra-cellular, but also intra-cellular. While not comprehensive, a simple two-compartment diffusion tensor model with water exchange was able to account qualitatively for the trends in positive and negative filtering efficiencies, while standard model imaging (SMI) without exchange could not. This two-compartment diffusion tensor model also demonstrated smaller AXR variances across subjects. When employing trace-weighted diffusion detection, AXR values were on the order of the R1 (=1/T1) of water at 3T for crossing fibers, while being less than R1 for parallel fibers.

CONCLUSION

Orientation-dependent AXR and σ values were observed when using multi-orientation filter and detection gradients in FEXI, indicating that WM FEXI models need to account for compartmental anisotropy. When using trace-weighted detection, AXR values were on the order of or less than R1, complicating the interpretation of FEXI results in WM in terms of biological exchange properties. These findings may contribute toward better understanding of FEXI results in WM.

Limitations and generalizations of the first order kinetics reaction expression for modeling diffusion-driven exchange: Implications on NMR exchange measurements

The study and modeling of water exchange in complex media using different applications of diffusion and relaxation magnetic resonance (MR) have been of interest in recent years. Most models attempt to describe this process using a first order kinetics expression, which is appropriate to describe chemical exchange; however, it may not be suitable to describe diffusion-driven exchange since it has no direct relationship to diffusion dynamics of water molecules. In this paper, these limitations are addressed through a more general exchange expression that does consider such important properties. This exchange fraction expression features a multi-exponential recovery at short times and a mono-exponential decay at long times, both of which are not captured by the first order kinetics expression. Furthermore, simplified exchange expressions containing partial information of the analyzed system's diffusion and relaxation processes and geometry are proposed, which can potentially be employed in already established estimation protocols. Finally, exchange fractions estimated from simulated MR data and derived here were compared, showing qualitative similarities but quantitative differences, suggesting that the features of the derived exchange fraction in this paper can be partially recovered by employing an existing estimation framework.

A translational MRI approach to validate acute axonal damage detection as an early event in multiple sclerosis

Axonal degeneration is a central pathological feature of multiple sclerosis and is closely associated with irreversible clinical disability. Current noninvasive methods to detect axonal damage in vivo are limited in their specificity and clinical applicability, and by the lack of proper validation. We aimed to validate an MRI framework based on multicompartment modeling of the diffusion signal (AxCaliber) in rats in the presence of axonal pathology, achieved through injection of a neurotoxin damaging the neuronal terminal of axons. We then applied the same MRI protocol to map axonal integrity in the brain of multiple sclerosis relapsing-remitting patients and age-matched healthy controls. AxCaliber is sensitive to acute axonal damage in rats, as demonstrated by a significant increase in the mean axonal caliber along the targeted tract, which correlated with neurofilament staining. Electron microscopy confirmed that increased mean axonal diameter is associated with acute axonal pathology. In humans with multiple sclerosis, we uncovered a diffuse increase in mean axonal caliber in most areas of the normal-appearing white matter, preferentially affecting patients with short disease duration. Our results demonstrate that MRI-based axonal diameter mapping is a sensitive and specific imaging biomarker that links noninvasive imaging contrasts with the underlying biological substrate, uncovering generalized axonal damage in multiple sclerosis as an early event.

Repeatability and reproducibility of apparent exchange rate measurements in yeast cell phantoms using filter-exchange imaging

OBJECTIVES

Development of a protocol for validation and quality assurance of filter-exchange imaging (FEXI) pulse sequences with well-defined and reproducible phantoms.

MATERIALS AND METHODS

A FEXI pulse sequence was implemented on a 7 T preclinical MRI scanner. Six experiments in three different test categories were established for sequence validation, demonstration of the reproducibility of phantoms and the measurement of induced changes in the apparent exchange rate (AXR). First, an ice-water phantom was used to investigate the consistency of apparent diffusion coefficient (ADC) measurements with different diffusion filters. Second, yeast cell phantoms were utilized to validate the determination of the AXR in terms of repeatability (same phantom and session), reproducibility (separate but comparable phantoms in different sessions) and directionality of diffusion encodings. Third, the yeast cell phantoms were, furthermore, used to assess potential AXR bias because of altered cell density and temperature. In addition, a treatment experiment with aquaporin inhibitors was performed to evaluate the influence of these compounds on the cell membrane permeability in yeast cells.

RESULTS

FEXI-based ADC measurements of an ice-water phantom were performed for three different filter strengths, showed good agreement with the literature value of 1.099 × 10-3 mm2/s and had a maximum coefficient of variation (CV) of 0.55% within the individual filter strengths. AXR estimation in a single yeast cell phantom and imaging session with five repetitions resulted in an overall mean value of (1.49 ± 0.05) s-1 and a CV of 3.4% between the chosen regions of interest. For three separately prepared phantoms, AXR measurements resulted in a mean value of (1.50 ± 0.04) s-1 and a CV of 2.7% across the three phantoms, demonstrating high reproducibility. Across three orthogonal diffusion directions, a mean value of (1.57 ± 0.03) s-1 with a CV of 1.9% was detected, consistent with isotropy of AXR in yeast cells. Temperature and AXR were linearly correlated (R2 = 0.99) and an activation energy EA of 37.7 kJ/mol was determined by Arrhenius plot. Furthermore, a negative correlation was found between cell density (as determined by the reference ADC/fe) and AXR (R2 = 0.95). The treatment experiment resulted in significantly decreased AXR values at different temperatures in the treated sample compared to the untreated control indicating an inhibiting effect.

CONCLUSIONS

Using ice-water and yeast cell-based phantoms, a protocol for the validation of FEXI pulse sequences was established for the assessment of stability, repeatability, reproducibility and directionality. In addition, a strong dependence of AXR on cell density and temperature was shown. As AXR is an emerging novel imaging biomarker, the suggested protocol will be useful for quality assurance of AXR measurements within a study and potentially across multiple sites.

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

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

Extra-axonal contribution to double diffusion encoding-based pore size estimates in the corticospinal tract

OBJECTIVE

To study the origin of compartment size overestimation in double diffusion encoding MRI (DDE) in vivo experiments in the human corticospinal tract. Here, the extracellular space is hypothesized to be the origin of the DDE signal. By exploiting the DDE sensitivity to pore shape, it could be possible to identify the origin of the measured signal. The signal difference between parallel and perpendicular diffusion gradient orientation can indicate if a compartment is regular or eccentric in shape. As extracellular space can be considered an eccentric compartment, a positive difference would mean a high contribution to the compartment size estimates.

MATERIALS AND METHODS

Computer simulations using MISST and in vivo experiments in eight healthy volunteers were performed. DDE experiments using a double spin-echo preparation with eight perpendicular directions were measured in vivo. The difference between parallel and perpendicular gradient orientations was analyzed using a Wilcoxon signed-rank test and a Mann-Whitney U test.

RESULTS

Simulations and MR experiments showed a statistically significant difference between parallel and perpendicular diffusion gradient orientation signals ([Formula: see text]).

CONCLUSION

The results suggest that the DDE-based size estimate may be considerably influenced by the extra-axonal compartment. However, the experimental results are also consistent with purely intra-axonal contributions in combination with a large fiber orientation dispersion.

Blood-brain barrier water exchange measurements using FEXI: Impact of modeling paradigm and relaxation time effects

PURPOSE

To evaluate potential modeling paradigms and the impact of relaxation time effects on human blood-brain barrier (BBB) water exchange measurements using FEXI (BBB-FEXI), and to quantify the accuracy, precision, and repeatability of BBB-FEXI exchange rate estimates at 3 T $$ \mathrm{T} $$ .

METHODS

Three modeling paradigms were evaluated: (i) the apparent exchange rate (AXR) model; (ii) a two-compartment model (2 CM $$ 2\mathrm{CM} $$ ) explicitly representing intra- and extravascular signal components, and (iii) a two-compartment model additionally accounting for finite compartmentalT 1 $$ {\mathrm{T}}_1 $$ andT 2 $$ {\mathrm{T}}_2 $$ relaxation times (2 CM r $$ 2{\mathrm{CM}}_r $$ ). Each model had three free parameters. Simulations quantified biases introduced by the assumption of infinite relaxation times in the AXR and2 CM $$ 2\mathrm{CM} $$ models, as well as the accuracy and precision of all three models. The scan-rescan repeatability of all paradigms was quantified for the first time in vivo in 10 healthy volunteers (age range 23-52 years; five female).

RESULTS

The assumption of infinite relaxation times yielded exchange rate errors in simulations up to 42%/14% in the AXR/2 CM $$ 2\mathrm{CM} $$ models, respectively. Accuracy was highest in the compartmental models; precision was best in the AXR model. Scan-rescan repeatability in vivo was good for all models, with negligible bias and repeatability coefficients in grey matter ofRC AXR = 0 . 43 $$ {\mathrm{RC}}_{\mathrm{AXR}}=0.43 $$  s - 1 $$ {\mathrm{s}}^{-1} $$ ,RC 2 CM = 0 . 51 $$ {\mathrm{RC}}_{2\mathrm{CM}}=0.51 $$  s - 1 $$ {\mathrm{s}}^{-1} $$ , andRC 2 CM r = 0 . 61 $$ {\mathrm{RC}}_{2{\mathrm{CM}}_r}=0.61 $$  s - 1 $$ {\mathrm{s}}^{-1} $$ .

CONCLUSION

Compartmental modelling of BBB-FEXI signals can provide accurate and repeatable measurements of BBB water exchange; however, relaxation time and partial volume effects may cause model-dependent biases.

A Method for In-Vivo Mapping of Axonal Diameter Distributions in the Human Brain Using Diffusion-Based Axonal Spectrum Imaging (AxSI)

In this paper we demonstrate a generalized and simplified pipeline called axonal spectrum imaging (AxSI) for in-vivo estimation of axonal characteristics in the human brain. Whole-brain estimation of the axon diameter, in-vivo and non-invasively, across all fiber systems will allow exploring uncharted aspects of brain structure and function relations with emphasis on connectivity and connectome analysis. While axon diameter mapping is important in and of itself, its correlation with conduction velocity will allow, for the first time, the explorations of information transfer mechanisms within the brain. We demonstrate various well-known aspects of axonal morphometry (e.g., the corpus callosum axon diameter variation) as well as other aspects that are less explored (e.g., axon diameter-based separation of the superior longitudinal fasciculus into segments). Moreover, we have created an MNI based mean axon diameter map over the entire brain for a large cohort of subjects providing the reference basis for future studies exploring relation between axon properties, its connectome representation, and other functional and behavioral aspects of the brain.

Vascular-water-exchange MRI (VEXI) enables the detection of subtle AXR alterations in Alzheimer's disease without MRI contrast agent, which may relate to BBB integrity

Blood-brain barrier (BBB) impairment is an important pathophysiological process in Alzheimer's disease (AD) and a potential biomarker for early diagnosis of AD. However, most current neuroimaging methods assessing BBB function need the injection of exogenous contrast agents (or tracers), which limits the application of these methods in a large population. In this study, we aim to explore the feasibility of vascular water exchange MRI (VEXI), a diffusion-MRI-based method proposed to assess the BBB permeability to water molecules without using a contrast agent, in the detection of the BBB breakdown in AD. We tested VEXI on a 3T MRI scanner on three groups: AD patients (AD group), mild cognitive impairment (MCI) patients due to AD (MCI group), and the age-matched normal cognition subjects (NC group). Interestingly, we find that the apparent water exchange across the BBB (AXR_BBB) measured by VEXI shows higher values in MCI compared with NC, and this higher AXR_BBB happens specifically in the hippocampus. This increase in AXR_BBB value gets larger and extends to more brain regions (medial orbital frontal cortex and thalamus) from MCI group to the AD group. Furthermore, we find that the AXR_BBB values of these three regions is correlated significantly with the impairment of respective cognitive domains independent of age, sex and education. These results suggest VEXI is a promising method to assess the BBB breakdown in AD.

Filter exchange imaging with crusher gradient modelling detects increased blood-brain barrier water permeability in response to mild lung infection

Blood-brain barrier (BBB) dysfunction occurs in many brain diseases, and there is increasing evidence to suggest that it is an early process in dementia which may be exacerbated by peripheral infection. Filter-exchange imaging (FEXI) is an MRI technique for measuring trans-membrane water exchange. FEXI data is typically analysed using the apparent exchange rate (AXR) model, yielding estimates of the AXR. Crusher gradients are commonly used to remove unwanted coherence pathways arising from longitudinal storage pulses during the mixing period. We first demonstrate that when using thin slices, as is needed for imaging the rodent brain, crusher gradients result in underestimation of the AXR. To address this, we propose an extended crusher-compensated exchange rate (CCXR) model to account for diffusion-weighting introduced by the crusher gradients, which is able to recover ground truth values of BBB water exchange (k_in) in simulated data. When applied to the rat brain, k_in estimates obtained using the CCXR model were 3.10 s^-1 and 3.49 s^-1 compared to AXR estimates of 1.24 s^-1 and 0.49 s^-1 for slice thicknesses of 4.0 mm and 2.5 mm respectively. We then validated our approach using a clinically relevant Streptococcus pneumoniae lung infection. We observed a significant 70 ± 10% increase in BBB water exchange in rats during active infection (k_in = 3.78 ± 0.42 s^-1) compared to before infection (k_in = 2.72 ± 0.30 s^-1; p = 0.02). The BBB water exchange rate during infection was associated with higher levels of plasma von Willebrand factor (VWF), a marker of acute vascular inflammation. We also observed 42% higher expression of perivascular aquaporin-4 (AQP4) in infected animals compared to non-infected controls, while levels of tight junction proteins remain consistent between groups. In summary, we propose a modelling approach for FEXI data which removes the bias in estimated water-exchange rates associated with the use of crusher gradients. Using this approach, we demonstrate the impact of peripheral infection on BBB water exchange, which appears to be mediated by endothelial dysfunction and associated with an increase in perivascular AQP4.

Measuring water exchange on a preclinical MRI system using filter exchange and diffusion time dependent kurtosis imaging

PURPOSE

Filter exchange imaging (FEXI) and diffusion time (t)-dependent diffusion kurtosis imaging (DKI(t)) are both sensitive to water exchange between tissue compartments. The restrictive effects of tissue microstructure, however, introduce bias to the exchange rate obtained by these two methods, as their interpretation conventionally rely on the Kärger model of barrier limited exchange between Gaussian compartments. Here, we investigated whether FEXI and DKI(t) can provide comparable exchange rates in ex vivo mouse brains.

THEORY AND METHODS

FEXI and DKI(t) data were acquired from ex vivo mouse brains on a preclinical MRI system. Phase cycling and negative slice prewinder gradients were used to minimize the interferences from imaging gradients.

RESULTS

In the corpus callosum, apparent exchange rate (AXR) from FEXI correlated with the exchange rate (the inverse of exchange time, 1/τex ) from DKI(t) along the radial direction. In comparison, discrepancies between FEXI and DKI(t) were found in the cortex due to low filter efficiency and confounding effects from tissue microstructure.

CONCLUSION

The results suggest that FEXI and DKI(t) are sensitive to the same exchange processes in white matter when separated from restrictive effects of microstructure. The complex microstructure in gray matter, with potential exchange among multiple compartments and confounding effects of microstructure, still pose a challenge for FEXI and DKI(t).

Transmembrane water-efflux rate measured by magnetic resonance imaging as a biomarker of the expression of aquaporin-4 in gliomas

The water-selective channel protein aquaporin-4 (AQP4) contributes to the migration and proliferation of gliomas, and to their resistance to therapy. Here we show, in glioma cell cultures, in subcutaneous and orthotopic gliomas in rats, and in glioma tumours in patients, that transmembrane water-efflux rate is a sensitive biomarker of AQP4 expression and can be measured via conventional dynamic-contrast-enhanced magnetic resonance imaging. Water-efflux rates correlated with stages of glioma proliferation as well as with changes in the heterogeneity of intra-tumoural and inter-tumoural AQP4 in rodent and human gliomas following treatment with temozolomide and with the AQP4 inhibitor TGN020. Regions with low water-efflux rates contained higher fractions of stem-like slow-cycling cells and therapy-resistant cells, suggesting that maps of water-efflux rates could be used to identify gliomas that are resistant to therapies.

In vivo methods for imaging blood-brain barrier function and dysfunction

The blood-brain barrier (BBB) is the interface between the central nervous system and systemic circulation. It tightly regulates what enters and is removed from the brain parenchyma and is fundamental in maintaining brain homeostasis. Increasingly, the BBB is recognised as having a significant role in numerous neurological disorders, ranging from acute disorders (traumatic brain injury, stroke, seizures) to chronic neurodegeneration (Alzheimer's disease, vascular dementia, small vessel disease). Numerous approaches have been developed to study the BBB in vitro, in vivo, and ex vivo. The complex multicellular structure and effects of disease are difficult to recreate accurately in vitro, and functional aspects of the BBB cannot be easily studied ex vivo. As such, the value of in vivo methods to study the intact BBB cannot be overstated. This review discusses the structure and function of the BBB and how these are affected in diseases. It then discusses in depth several established and novel methods for imaging the BBB in vivo, with a focus on MRI, nuclear imaging, and high-resolution intravital fluorescence microscopy.

Diffusion MRI with pulsed and free gradient waveforms: Effects of restricted diffusion and exchange

Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm and exchange rates in the simulated range of 0 to 20 s - 1 , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.

What does FEXI measure?

Filter-exchange imaging (FEXI) has already been utilized in several biomedical studies for evaluating the permeability of cell membranes. The method relies on suppressing the extracellular signal using strong diffusion weighting (the mobility filter causing a reduction in the overall diffusivity) and monitoring the subsequent diffusivity recovery. Using Monte Carlo simulations, we demonstrate that FEXI is sensitive not uniquely to the transcytolemmal exchange but also to the geometry of involved compartments: complex geometry offers locations where spins remain unaffected by the mobility filter; moving to other locations afterwards, such spins contribute to the diffusivity recovery without actually permeating any membrane. This exchange mechanism is a warning for those who aim to use FEXI in complex media such as brain gray matter and opens wide scope for investigation towards crystallizing the genuine membrane permeation and characterizing the compartment geometry.

Improving MR cell size imaging by inclusion of transcytolemmal water exchange

The goal of the current study is to include transcytolemmal water exchange in MR cell size imaging using the IMPULSED model for more accurate characterization of tissue cellular properties (e.g., apparent volume fraction of intracellular space v in ) and quantification of indicators of transcytolemmal water exchange. We propose a heuristic model that incorporates transcytolemmal water exchange into a multicompartment diffusion-based method (IMPULSED) that was developed previously to extract microstructural parameters (e.g., mean cell size d and apparent volume fraction of intracellular space v in ) assuming no water exchange. For t diff ≤ 5 ms, the water exchange can be ignored, and the signal model is the same as the IMPULSED model. For t diff ≥ 30 ms, we incorporated the modified Kärger model that includes both restricted diffusion and exchange between compartments. Using simulations and previously published in vitro cell data, we evaluated the accuracy and precision of model-derived parameters and determined how they are dependent on SNR and imaging parameters. The joint model provides more accurate d values for cell sizes ranging from 10 to 12 microns when water exchange is fast (e.g., intracellular water pre-exchange lifetime τ in ≤ 100 ms) than IMPULSED, and reduces the bias of IMPULSED-derived estimates of v in , especially when water exchange is relatively slow (e.g., τ in > 200 ms). Indicators of transcytolemmal water exchange derived from the proposed joint model are linearly correlated with ground truth τ in values and can detect changes in cell membrane permeability induced by saponin treatment in murine erythroleukemia cancer cells. Our results suggest this joint model not only improves the accuracy of IMPULSED-derived microstructural parameters, but also provides indicators of water exchange that are usually ignored in diffusion models of tissues.

Role of Transmembrane Water Exchange in Glioma Invasion/Migration: In Vivo Preclinical Study by Relaxometry at Very Low Magnetic Field

This work shows that the longitudinal relaxation differences observed at very low magnetic fields between invasion/migration and proliferation processes on glioma mouse models in vivo are related to differences in the transmembrane water exchange basically linked to the aquaporin expression changes. Three glioma mouse models were used: Glio6 and Glio96 as invasion/migration models and U87 as cell proliferation model. In vivo proton longitudinal relaxation-rate constants (R1) at very low fields were measured by fast field cycling NMR (FFC-NMR). The tumor contribution to the observed proton relaxation rate, R1tum (U87: 12.26 ± 0.64 s−1; Glio6: 3.76 ± 0.88 s−1; Glio96: 6.90 ± 0.64 s−1 at 0.01 MHz), and the intracellular water lifetime, τin (U87: 826 ± 19 ms; Glio6: 516 ± 8 ms; Glio96: 596 ± 15 ms), were found to be good diagnostic hallmarks to distinguish invasion/migration from proliferation (p < 0.01 and 0.001). Overexpression of AQP4 and AQP1 were assessed in invasion/migration models, highlighting the pathophysiological role of these two aquaporins in water exchange that, in turn, determine the lower values in the observed R1 relaxation rate constant in glioma invasion/migration. Overall, our findings demonstrate that τin and R1 (measured at very low fields) are relevant biomarkers, discriminating invasion/migration from proliferation in vivo. These results highlight the use of FFC-NMR and FFC-imaging to assess the efficiency of drugs that could modulate aquaporin functions.

Diffusion pore imaging in the presence of extraporal water

Diffusion-weighted imaging (DWI) is a powerful, non-invasive tool which is widely used in clinical routine. Mostly, apparent diffusion coefficient maps are acquired, which cannot be related directly to cellular structure. More recently it was shown that DWI is able to reconstruct pore shapes using a specialized magnetic field gradient scheme so that cell size distributions may be obtained. So far, artificial systems have been used for experimental demonstration without extraporal signal components and relatively low gradient amplitudes. The aim of this study was to investigate the feasibility of diffusion pore imaging in the presence of extraporal fluids and to develop correction methods for the effects arising from extraporal signal contributions. Monte Carlo simulations and validation experiments on a 14.1 T NMR spectrometer equipped with a dedicated diffusion probe head were performed. Both by using a filter gradient approach suppressing extraporal signal components as well as by using post-processing methods relying on the Gaussian phase approximation, it was possible to reconstruct pore space functions in the presence of extraporal fluids with little to no deviations from the expectations. These results may be a significant step towards application of diffusion pore imaging to biological samples.

Disentangling the effects of restriction and exchange with diffusion exchange spectroscopy

Diffusion exchange spectroscopy (DEXSY) is a multidimensional NMR technique that can reveal how water molecules exchange between compartments within heterogeneous media, such as biological tissue. Data from DEXSY experiments is typically processed using numerical inverse Laplace transforms (ILTs) to produce a diffusion-diffusion spectrum. A tacit assumption of this ILT approach is that the signal behavior is Gaussian - i.e., the spin echo intensity decays exponentially with the degree of diffusion weighting. The assumptions that underlie Gaussian signal behavior may be violated, however, depending on the gradient strength applied and the sample under study. We argue that non-Gaussian signal behavior due to restrictions is to be expected in the study of biological tissue using diffusion NMR. Further, we argue that this signal behavior can produce confounding features in the diffusion-diffusion spectra obtained from numerical ILTs of DEXSY data - entangling the effects of restriction and exchange. Specifically, restricted signal behavior can result in broadening of peaks and in the appearance of illusory exchanging compartments with distributed diffusivities, which pearl into multiple peaks if not highly regularized. We demonstrate these effects on simulated data. That said, we suggest the use of features in the signal acquisition domain that can be used to rapidly probe exchange without employing an ILT. We also propose a means to characterize the non-Gaussian signal behavior due to restrictions within a sample using DEXSY measurements with a near zero mixing time or storage interval. We propose a combined acquisition scheme to independently characterize restriction and exchange with various DEXSY measurements, which we term Restriction and Exchange from Equally-weighted Double and Single Diffusion Encodings (REEDS-DE). We test this method on ex vivo neonatal mouse spinal cord - a sample consisting primarily of gray matter - using a low-field, static gradient NMR system. In sum, we highlight critical shortcomings of prevailing DEXSY analysis methods that conflate the effects of restriction and exchange, and suggest a viable experimental approach to disentangle them.

The direction-dependence of apparent water exchange rate in human white matter

Transmembrane water exchange is a potential biomarker in the diagnosis and understanding of cancers, brain disorders, and other diseases. Filter-exchange imaging (FEXI), a special case of diffusion exchange spectroscopy adapted for clinical applications, has the potential to reveal different physiological water exchange processes. However, it is still controversial whether modulating the diffusion encoding gradient direction can affect the apparent exchange rate (AXR) measurements of FEXI in white matter (WM) where water diffusion shows strong anisotropy. In this study, we explored the diffusion-encoding direction dependence of FEXI in human brain white matter by performing FEXI with 20 diffusion-encoding directions on a clinical 3T scanner in-vivo. The results show that the AXR values measured when the gradients are perpendicular to the fiber orientation (0.77 ± 0.13 s ^- 1, mean ± standard deviation of all the subjects) are significantly larger than the AXR estimates when the gradients are parallel to the fiber orientation (0.33 ± 0.14 s ^- 1, p < 0.001) in WM voxels with coherently-orientated fibers. In addition, no significant correlation is found between AXRs measured along these two directions, indicating that they are measuring different water exchange processes. What's more, only the perpendicular AXR rather than the parallel AXR shows dependence on axonal diameter, indicating that the perpendicular AXR might reflect transmembrane water exchange between intra-axonal and extra-cellular spaces. Further finite difference (FD) simulations having three water compartments (intra-axonal, intra-glial, and extra-cellular spaces) to mimic WM micro-environments also suggest that the perpendicular AXR is more sensitive to the axonal water transmembrane exchange than parallel AXR. Taken together, our results show that AXR measured along different directions could be utilized to probe different water exchange processes in WM.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Spatiotemporal characterisation of ischaemic lesions in transient stroke animal models using diffusion free water elimination and mapping MRI with echo time dependence

BACKGROUND AND PURPOSE

The excess fluid as a result of vasogenic oedema and the subsequent tissue cavitation obscure the microstructural characterisation of ischaemic tissue by conventional diffusion and relaxometry MRI. They lead to a pseudo-normalisation of the water diffusivity and transverse relaxation time maps in the subacute and chronic phases of stroke. Within the context of diffusion MRI, the free water elimination and mapping method (FWE) with echo time dependence has been proposed as a promising approach to measure the amount of free fluid in brain tissue robustly and to eliminate its biasing effect on other biomarkers. In this longitudinal study of transient middle cerebral artery occlusion (MCAo) in the rat brain, we investigated the use of FWE MRI with echo time dependence for the characterisation of the tissue microstructure and explored the potential of the free water fraction as a novel biomarker of ischaemic tissue condition.

METHODS

Adult rats received a transient MCAo. Diffusion- and transverse relaxation-weighted MRI experiments were performed longitudinally, pre-occlusion and on days 1, 3, 4, 5, 6, 7 and 10 after MCAo on four rats. Histology was performed for non-stroke and 1, 3 and 10 days after MCAo on three different rats at each time point.

RESULTS

The free water fraction was homogeneously increased in the ischaemic cortex one day after stroke. Between three and ten days after stroke, the core of the ischaemic tissue showed a progressive normalisation in the amount of free water, whereas the inner and outer border zones of the ischaemic cortex depicted a large, monotonous increase with time. The specific lesions in brain sections were verified by H&E and immunostaining. The tissue-specific diffusion and relaxometry MRI metrics in the ischaemic cortex were significantly different compared to their conventional counterpart.

CONCLUSIONS

Our results demonstrate that the free water fraction in FWE MRI with echo time dependence is a valuable biomarker, sensitive to the progressive degeneration in ischaemic tissue. We showed that part of the heterogeneity previously observed in conventional parameter maps can be accounted for by a heterogeneous distribution of free water in the tissue. Our results suggest that the temporal evolution of the free fluid fraction map at the core and inner border zone can be associated with the pathological changes linked to the evolution of vasogenic oedema. Namely, the homogeneous increase in free water one day after stroke and its tendency to normalise in the core of the ischaemic cortex starting three days after stroke, followed by a progressive increase in free water at the inner border zone from three to ten days after stroke. Finally, the monotonous increase in free fluid in the outer border zone of the cortex reflects the formation of fluid-filled cysts.

Assessment of the effects of mimicking tissue microstructural properties on apparent diffusion coefficient and apparent exchange rate in diffusion MRI via a series of specially designed phantoms

PURPOSE

Diffusion MRI provides a valuable tool for imaging tissue microstructure. However, due to the lack of related experimental methods and specially designed phantoms, no experimental study has been conducted yet to quantitatively assess the effects of membrane permeability, intracellular volume fraction (IVF), and intracellular diffusivity on the apparent diffusion coefficient (ADC) obtained from diffusion weighted imaging (DWI), and the effects of membrane permeability on the apparent exchange rate (AXR) obtained from filter exchange imaging (FEXI).

METHODS

A series of phantoms with three adjustable parameters was designed to mimic tissue microstructural properties including membrane permeability, IVF, and intracellular diffusivity. Quantitative experiments were conducted to assess the effects of these properties on ADC and AXR. DWI scans were performed to obtain axial and radial ADC values. FEXI scans were performed to obtain AXR values.

RESULTS

Axial ADC values range from 1.148 μm² /ms to 2.157 μm² /ms, and radial ADC values range from 0.904 μm² /ms to 2.067 μm² /ms. Radial ADC decreased with a decrease in fiber permeability. Decreased axial and radial ADC values with increased intra-fiber volume fraction, and increased polyvinylpyrrolidone (PVP) concentration of the intra-fiber space were observed. AXR values range from 2.1 s^-1 to 4.9 s^-1 . AXR increases with fiber permeability.

CONCLUSION

The proposed phantoms can quantitatively evaluate the effects of mimicking tissue microstructural properties on ADC and AXR. This new phantom design provides a potential method for further understanding the biophysical mechanisms underlying the change in ADC and diffusion exchange.

SPHERIOUSLY? The challenges of estimating sphere radius non-invasively in the human brain from diffusion MRI

The Soma and Neurite Density Imaging (SANDI) three-compartment model was recently proposed to disentangle cylindrical and spherical geometries, attributed to neurite and soma compartments, respectively, in brain tissue. There are some recent advances in diffusion-weighted MRI signal encoding and analysis (including the use of multiple so-called 'b-tensor' encodings and analysing the signal in the frequency-domain) that have not yet been applied in the context of SANDI. In this work, using: (i) ultra-strong gradients; (ii) a combination of linear, planar, and spherical b-tensor encodings; and (iii) analysing the signal in the frequency domain, three main challenges to robust estimation of sphere size were identified: First, the Rician noise floor in magnitude-reconstructed data biases estimates of sphere properties in a non-uniform fashion. It may cause overestimation or underestimation of the spherical compartment size and density. This can be partly ameliorated by accounting for the noise floor in the estimation routine. Second, even when using the strongest diffusion-encoding gradient strengths available for human MRI, there is an empirical lower bound on the spherical signal fraction and radius that can be detected and estimated robustly. For the experimental setup used here, the lower bound on the sphere signal fraction was approximately 10%. We employed two different ways of establishing the lower bound for spherical radius estimates in white matter. The first, examining power-law relationships between the DW-signal and diffusion weighting in empirical data, yielded a lower bound of 7μm, while the second, pure Monte Carlo simulations, yielded a lower limit of 3μm and in this low radii domain, there is little differentiation in signal attenuation. Third, if there is sensitivity to the transverse intra-cellular diffusivity in cylindrical structures, e.g., axons and cellular projections, then trying to disentangle two diffusion-time-dependencies using one experimental parameter (i.e., change in frequency-content of the encoding waveform) makes spherical radii estimates particularly challenging. We conclude that due to the aforementioned challenges spherical radii estimates may be biased when the corresponding sphere signal fraction is low, which must be considered.

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study

This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI.

Time dependence in diffusion MRI predicts tissue outcome in ischemic stroke patients

Purpose:

Reperfusion therapy enables effective treatment of ischemic stroke presenting within 4–6 hours. However, tissue progression from ischemia to infarction is variable, and some patients benefit from treatment up until 24 hours. Improved imaging techniques are needed to identify these patients. Here, it was hypothesized that time dependence in diffusion MRI may predict tissue outcome in ischemic stroke.

Methods:

Diffusion MRI data were acquired with multiple diffusion times in five non-reperfused patients at 2, 9, and 100 days after stroke onset. Maps of “rate of kurtosis change” (k), mean kurtosis, ADC, and fractional anisotropy were derived. The ADC maps defined lesions, normal-appearing tissue, and the lesion tissue that would either be infarcted or remain viable by day 100. Diffusion parameters were compared (1) between lesions and normal-appearing tissue, and (2) between lesion tissue that would be infarcted or remain viable.

Results:

Positive values of k were observed within stroke lesions on day 2 (P = .001) and on day 9 (P = .023), indicating diffusional exchange. On day 100, high ADC values indicated infarction of 50 ± 20% of the lesion volumes. Tissue infarction was predicted by high k values both on day 2 (P = .026) and on day 9 (P = .046), by low mean kurtosis values on day 2 (P = .043), and by low fractional anisotropy values on day 9 (P = .029), but not by low ADC values.

Conclusions:

Diffusion time dependence predicted tissue outcome in ischemic stroke more accurately than the ADC, and may be useful for predicting reperfusion benefit.

Compartmental diffusion and microstructural properties of human brain gray and white matter studied with double diffusion encoding magnetic resonance spectroscopy of metabolites and water

Double diffusion encoding (DDE) of the water signal offers a unique ability to separate the effect of microscopic anisotropic diffusion in structural units of tissue from the overall macroscopic orientational distribution of cells. However, the specificity in detected microscopic anisotropy is limited as the signal is averaged over different cell types and across tissue compartments. Performing side-by-side water and metabolite DDE spectroscopic (DDES) experiments provides complementary measures from which intracellular and extracellular microscopic fractional anisotropies (μFA) and diffusivities can be estimated. Metabolites are largely confined to the intracellular space and therefore provide a benchmark for intracellular μFA and diffusivities of specific cell types. By contrast, water DDES measurements allow examination of the separate contributions to water μFA and diffusivity from the intra- and extracellular spaces, by using a wide range of b values to gradually eliminate the extracellular contribution. Here, we aimed to estimate tissue and compartment specific human brain microstructure by combining water and metabolites DDES experiments. We performed our DDES measurements in two brain regions that contain widely different amounts of white matter (WM) and gray matter (GM): parietal white matter (PWM) and occipital gray matter (OGM) in a total of 20 healthy volunteers at 7 Tesla. Metabolite DDES measurements were performed at b = 7199 s/mm², while water DDES measurements were performed with a range of b values from 918 to 7199 s/mm². The experimental framework we employed here resulted in a set of insights pertaining to the morphology of the intracellular and extracellular spaces in both gray and white matter. Results of the metabolite DDES experiments in both PWM and OGM suggest a highly anisotropic intracellular space within neurons and glia, with the possible exception of gray matter glia. The water μFA obtained from the DDES results at high b values in both regions converged with that of the metabolite DDES, suggesting that the signal from the extracellular space is indeed effectively suppressed at the highest b value. The μFA measured in the OGM significantly decreased at lower b values, suggesting a considerably lower anisotropy of the extracellular space in GM compared to WM. In PWM, the water μFA remained high even at the lowest b value, indicating a high degree of organization in the interstitial space in WM. Tortuosity values in the cytoplasm for water and tNAA, obtained with correlation analysis of microscopic parallel diffusivity with respect to GM/WM tissue fraction in the volume of interest, are remarkably similar for both molecules, while exhibiting a clear difference between gray and white matter, suggesting a more crowded cytoplasm and more complex cytomorphology of neuronal cell bodies and dendrites in GM than those found in long-range axons in WM.

Apparent exchange rate imaging: On its applicability and the connection to the real exchange rate

PURPOSE

Water exchange between the intracellular and extracellular space can be measured using apparent exchange rate (AXR) imaging. The aim of this study was to investigate the relationship between the measured AXR and the geometry of diffusion restrictions, membrane permeability, and the real exchange rate, as well as to explore the applicability of AXR for typical human measurement settings.

METHODS

The AXR measurements and the underlying exchange rates were simulated using the Monte Carlo method with different geometries, size distributions, packing densities, and a broad range of membrane permeabilities. Furthermore, the influence of SNR and sequence parameters was analyzed.

RESULTS

The estimated AXR values correspond to the simulated values and show the expected proportionality to membrane permeability, except for fast exchange (ie, AXR > 20 - 30 s - 1 ) and small packing densities. Moreover, it was found that the duration of the filter gradient must be shorter than 2 · AX R - 1 . In cell size and permeability distributions, AXR depends on the average surface-to-volume ratio, permeability, and the packing density. Finally, AXR can be reliably determined in the presence of orientation dispersion in axon-like structures with sufficient gradient sampling (ie, 30 gradient directions).

CONCLUSION

Currently used experimental settings for in vivo human measurements are well suited for determining AXR, with the exception of single-voxel analysis, due to limited SNR. The detection of changes in membrane permeability in diseased tissue is nonetheless challenging because of the AXR dependence on further factors, such as packing density and geometry, which cannot be disentangled without further knowledge of the underlying cell structure.

Gradient waveform design for tensor-valued encoding in diffusion MRI

Diffusion encoding along multiple spatial directions per signal acquisition can be described in terms of a b-tensor. The benefit of tensor-valued diffusion encoding is that it unlocks the 'shape of the b-tensor' as a new encoding dimension. By modulating the b-tensor shape, we can control the sensitivity to microscopic diffusion anisotropy which can be used as a contrast mechanism; a feature that is inaccessible by conventional diffusion encoding. Since imaging methods based on tensor-valued diffusion encoding are finding an increasing number of applications we are prompted to highlight the challenge of designing the optimal gradient waveforms for any given application. In this review, we first establish the basic design objectives in creating field gradient waveforms for tensor-valued diffusion MRI. We also survey additional design considerations related to limitations imposed by hardware and physiology, potential confounding effects that cannot be captured by the b-tensor, and artifacts related to the diffusion encoding waveform. Throughout, we discuss the expected compromises and tradeoffs with an aim to establish a more complete understanding of gradient waveform design and its impact on accurate measurements and interpretations of data.

The sensitivity of diffusion MRI to microstructural properties and experimental factors

Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.

Toward Quantification: Microstructure and Magnetic Resonance Fingerprinting

Quantitative magnetic resonance imaging (MRI) is a long-standing challenge. We advocate that the origin of the problem is the simplification applied in commonly used models of the MRI signal relation to the target parameters of biological tissues. Two research fields are briefly reviewed as ways to respond to the challenge of quantitative MRI, both experiencing an exponential growth right now. Microstructure MRI strives to build physiology-based models from cells to signal and, given the signal, back to the cells again. Magnetic resonance fingerprinting aims at efficient simultaneous determination of multiple signal parameters. The synergy of these yet disjoined approaches promises truly quantitative MRI with specific target-oriented diagnostic tools rather than universal imaging methods.

In vivo observation and biophysical interpretation of time-dependent diffusion in human cortical gray matter

The dependence of the diffusion MRI signal on the diffusion time t is a hallmark of tissue microstructure at the scale of the diffusion length. Here we measure the time-dependence of the mean diffusivity D(t) and mean kurtosis K(t) in cortical gray matter and in 25 ​gray matter sub-regions, in 10 healthy subjects. Significant diffusivity and kurtosis time-dependence is observed for t=21.2-100 ​ms, and is characterized by a power-law tail ∼t-ϑ with dynamical exponent ϑ. To interpret our measurements, we systematize the relevant scenarios and mechanisms for diffusion time-dependence in the brain. Using the effective medium theory formalism, we derive an exact relation between the power-law tails in D(t) and K(t). The estimated dynamical exponent ϑ≃1/2 in both D(t) and K(t) is consistent with one-dimensional diffusion in the presence of randomly positioned restrictions along neurites. We analyze the short-range disordered statistics of synapses on axon collaterals in the cortex, and perform one-dimensional Monte Carlo simulations of diffusion restricted by permeable barriers with a similar randomness in their placement, to confirm the ϑ=1/2 exponent. In contrast, the Kärger model of exchange is less consistent with the data since it does not capture the diffusivity time-dependence, and the estimated exchange time from K(t) falls below our measured t-range. Although we cannot exclude exchange as a contributing factor, we argue that structural disorder along neurites is mainly responsible for the observed time-dependence of diffusivity and kurtosis. Our observation and theoretical interpretation of the t-1/2 tail in D(t) and K(t) altogether establish the sensitivity of a macroscopic MRI signal to micrometer-scale structural heterogeneities along neurites in human gray matter in vivo.

Double diffusion encoding and applications for biomedical imaging

Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.

Machine learning based white matter models with permeability: An experimental study in cuprizone treated in-vivo mouse model of axonal demyelination

The intra-axonal water exchange time (τ_i), a parameter associated with axonal permeability, could be an important biomarker for understanding and treating demyelinating pathologies such as Multiple Sclerosis. Diffusion-Weighted MRI (DW-MRI) is sensitive to changes in permeability; however, the parameter has so far remained elusive due to the lack of general biophysical models that incorporate it. Machine learning based computational models can potentially be used to estimate such parameters. Recently, for the first time, a theoretical framework using a random forest (RF) regressor suggests that this is a promising new approach for permeability estimation. In this study, we adopt such an approach and for the first time experimentally investigate it for demyelinating pathologies through direct comparison with histology. We construct a computational model using Monte Carlo simulations and an RF regressor in order to learn a mapping between features derived from DW-MRI signals and ground truth microstructure parameters. We test our model in simulations, and find strong correlations between the predicted and ground truth parameters (intra-axonal volume fraction f: R² =0.99, τ_i: R² =0.84, intrinsic diffusivity d: R² =0.99). We then apply the model in-vivo, on a controlled cuprizone (CPZ) mouse model of demyelination, comparing the results from two cohorts of mice, CPZ (N=8) and healthy age-matched wild-type (WT, N=8). We find that the RF model estimates sensible microstructure parameters for both groups, matching values found in literature. Furthermore, we perform histology for both groups using electron microscopy (EM), measuring the thickness of the myelin sheath as a surrogate for exchange time. Histology results show that our RF model estimates are very strongly correlated with the EM measurements (ρ = 0.98 for f, ρ = 0.82 for τ_i). Finally, we find a statistically significant decrease in τ_i in all three regions of the corpus callosum (splenium/genu/body) of the CPZ cohort (<τ_i>=310ms/330ms/350ms) compared to the WT group (<τ_i>=370ms/370ms/380ms). This is in line with our expectations that τ_i is lower in regions where the myelin sheath is damaged, as axonal membranes become more permeable. Overall, these results demonstrate, for the first time experimentally and in vivo, that a computational model learned from simulations can reliably estimate microstructure parameters, including the axonal permeability .

Pulse sequences for measuring exchange rates between proton species: From unlocalised NMR spectroscopy to chemical exchange saturation transfer imaging

Within the field of NMR spectroscopy, the study of chemical exchange processes through saturation transfer techniques has a long history. In the context of MRI, chemical exchange techniques have been adapted to increase the sensitivity of imaging to small fractions of exchangeable protons, including the labile protons of amines, amides and hydroxyls. The MR contrast is generated by frequency-selective irradiation of the labile protons, which results in a reduction of the water signal associated with transfer of the labile protons' saturated magnetization to the protons of the surrounding free water. The signal intensity depends on the rate of chemical exchange and the concentration of labile protons as well as on the properties of the irradiation field. This methodology is referred to as CEST (chemical exchange saturation transfer) imaging. Applications of CEST include imaging of molecules with short transverse relaxation times and mapping of physiological parameters such as pH, temperature, buffer concentration and chemical composition due to the dependency of this chemical exchange effect on all these parameters. This article aims to describe these effects both theoretically and experimentally. In depth analysis and mathematical modelling are provided for all pulse sequences designed to date to measure the chemical exchange rate. Importantly, it has become clear that the background signal from semi-solid protons and the presence of the Nuclear Overhauser Effect (NOE), either through direct dipole-dipole mechanisms or through exchange-relayed signals, complicates the analysis of CEST effects. Therefore, advanced methods to suppress these confounding factors have been developed, and these are also reviewed. Finally, the experimental work conducted both in vitro and in vivo is discussed and the progress of CEST imaging towards clinical practice is presented.

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

INTRODUCTION

To combine numerical simulations, in vitro and in vivo experiments to evaluate the feasibility of measuring diffusion exchange across the cell membrane with diffusion exchange spectroscopy (DEXSY).

METHODS

DEXSY acquisitions were simulated over a range of permeabilities in nerve tissue and yeast substrates. In vitro measurements were performed in a yeast substrate and in vivo measurements in mouse tumor xenograft models, all at 9.4 T.

RESULTS

Diffusion exchange was observed in simulations over a physiologically relevant range of cell permeability values. In vitro and in vivo measures also provided evidence of diffusion exchange, which was quantified with the Diffusion Exchange Index (DEI).

CONCLUSIONS

Our findings provide preliminary evidence that DEXSY can be used to make in vivo measurements of diffusion exchange and cell membrane permeability.

Microstructure with diffusion MRI: what scale we are sensitive to?

Diffusion-weighted MRI is the forerunner of the rapidly developed microstructure MRI (μMRI) aimed at in vivo evaluation of the cellular tissue architecture. This brief review focuses on the spatiotemporal scales of the microstructure that are accessible using different diffusion MRI techniques and the need to weight the measurability against the interpretability of results. Diffusion phenomena and models are first classified in two-dimensional space (the q-t-plane) of the measurement with narrow gradient pulses. Three-dimensional parameter space of the Stejskal-Tanner diffusion weighting adds more phenomena to this collection. Modern measurement techniques with larger number of parameters are briefly discussed under the overarching idea of diffusion weighting matching the geometry of the targeted cell species.

Real-time measurement of diffusion exchange rate in biological tissue

Diffusion exchange spectroscopy (DEXSY) provides a means to isolate the signal attenuation associated with exchange from other sources of signal loss. With the total diffusion weighting b₁+b₂=b_s held constant, DEXSY signals acquired with b₁=0 or b₂=0 have no exchange weighting, while a DEXSY signal acquired with b₁=b₂ has maximal exchange weighting. The exchange rate can be estimated by fitting a diffusion exchange model to signals acquired with variable mixing times. Conventionally, acquired signals are normalized by a signal with b₁=0 and b₂=0 to remove the decay due to spin-lattice relaxation. Instead, division by a signal with equal b_s but b₁=0 or b₂=0 reduces spin-lattice relaxation weighting of the apparent exchange rate (AXR). Furthermore, apparent diffusion-weighted R₁ relaxation rates can be estimated from non-exchange-weighted DEXSY signals. Estimated R₁ values are utilized to remove signal decay due to spin-lattice relaxation from exchange-weighted signals, permitting a more precise estimate of AXR with less data. Data reduction methods are proposed and tested with regards to statistical accuracy and precision of AXR estimates on simulated and experimental data. Simulations show that the methods are capable of accurately measuring the ground-truth exchange rate. The methods remain accurate even when the assumption that DEXSY signals attenuate with b is violated, as occurs for restricted diffusion. Experimental data was collected from fixed neonatal mouse spinal cord samples at 25 and 7°C using the strong static magnetic field gradient produced by a single-sided permanent magnet (i.e., an NMR MOUSE). The most rapid method for exchange measurements requires only five data points (an 80 s experiment as implemented) and achieves a similar level of accuracy and precision to the baseline method using 44 data points. This represents a significant improvement in acquisition speed, overcoming a barrier which has limited the use of DEXSY on living specimen.

Dedicated diffusion phantoms for the investigation of free water elimination and mapping: insights into the influence of T2 relaxation properties

Conventional diffusion-weighted (DW) MRI suffers from free water contamination due to the finite voxel size. The most common case of free water contamination occurs with cerebrospinal fluid (CSF) in voxels located at the CSF-tissue interface, such as at the ventricles in the human brain. Another case refers to intra-tissue free water as in vasogenic oedema. In order to avoid the bias in diffusion metrics, several multi-compartment methods have been introduced, which explicitly model the presence of a free water compartment. However, fitting multi-compartment models in DW MRI represents a well known ill conditioned problem. Although during the last decade great effort has been devoted to mitigating this estimation problem, the research field remains active. The aim of this work is to introduce the design, characterise the NMR properties and demonstrate the use of two dedicated anisotropic diffusion fibre phantoms, useful for the study of free water elimination (FWE) and mapping models. In particular, we investigate the recently proposed FWE diffusion tensor imaging approach, which takes explicit account of differences in the transverse relaxation times between the free water and tissue compartments.

Constant gradient FEXSY: A time-efficient method for measuring exchange

Filter-Exchange NMR Spectroscopy (FEXSY) is a method for measurement of apparent transmembranal water exchange rates. The experiment is comprised of two co-linear sequential pulsed-field gradient (PFG) blocks, separated by a mixing period in which exchange takes place. The first block remains constant and serves as a diffusion-based filter that removes signal coming from fast-diffusing water. The mixing time and the gradient area (q-value) of the second block are varied on repeated iterations to produce a 2D data set that is analyzed using a bi-compartmental model which assumes that intra- and extra-cellular water are slow and fast diffusing, respectively. Here we suggest a variant of the FEXSY method in which measurements for different mixing times are taken at a constant gradient. This Constant Gradient FEXSY (CG-FEXSY) allows for the determination of the exchange rate by using a smaller 1D data set, so that the same information can be gathered during a considerably shorter scan time. Furthermore, in the limit of high diffusion weighting, such that the extra-cellular water signal is removed while the intra-cellular signal is retained, CG-FEXSY also allows for determination of the intra-cellular mean residence time (MRT). The theoretical results are validated on a living yeast cells sample and on a fixed porcine optic nerve, where the values obtained from the two methods are shown to be in agreement.

Magnetic resonance measurements of cellular and sub-cellular membrane structures in live and fixed neural tissue

We develop magnetic resonance (MR) methods for real-time measurement of tissue microstructure and membrane permeability of live and fixed excised neonatal mouse spinal cords. Diffusion and exchange MR measurements are performed using the strong static gradient produced by a single-sided permanent magnet. Using tissue delipidation methods, we show that water diffusion is restricted solely by lipid membranes. Most of the diffusion signal can be assigned to water in tissue which is far from membranes. The remaining 25% can be assigned to water restricted on length scales of roughly a micron or less, near or within membrane structures at the cellular, organelle, and vesicle levels. Diffusion exchange spectroscopy measures water exchanging between membrane structures and free environments at 100 s^-1.

In vivo magnetic resonance imaging and spectroscopy. Technological advances and opportunities for applications continue to abound

Over the past decades, the field of in vivo magnetic resonance (MR) has built up an impressive repertoire of data acquisition and analysis technologies for anatomical, functional, physiological, and molecular imaging, the description of which requires many book volumes. As such it is impossible for a few authors to have an authoritative overview of the field and for a brief article to be inclusive. We will therefore focus mainly on data acquisition and attempt to give some insight into the principles underlying current advanced methods in the field and the potential for further innovation. In our view, the foreseeable future is expected to show continued rapid progress, for instance in imaging of microscopic tissue properties in vivo, assessment of functional and anatomical connectivity, higher resolution physiologic and metabolic imaging, and even imaging of receptor binding. In addition, acquisition speed and information content will continue to increase due to the continuous development of approaches for parallel imaging (including simultaneous multi-slice imaging), compressed sensing, and MRI fingerprinting. Finally, artificial intelligence approaches are becoming more realistic and will have a tremendous effect on both acquisition and analysis strategies. Together, these developments will continue to provide opportunity for scientific discovery and, in combination with large data sets from other fields such as genomics, allow the ultimate realization of precision medicine in the clinic.

In vivo magnetic resonance imaging of treatment-induced apoptosis

Imaging apoptosis could provide an early and specific means to monitor tumor responses to treatment. To date, despite numerous attempts to develop molecular imaging approaches, there is still no widely-accepted and reliable method for in vivo imaging of apoptosis. We hypothesized that the distinct cellular morphologic changes associated with treatment-induced apoptosis, such as cell shrinkage, cytoplasm condensation, and DNA fragmentation, can be detected by temporal diffusion spectroscopy imaging (TDSI). Cetuximab-induced apoptosis was assessed in vitro and in vivo with cetuximab-sensitive (DiFi) and insensitive (HCT-116) human colorectal cancer cell lines by TDSI. TDSI findings were complemented by flow cytometry and immunohistochemistry. Cell cycle analysis and flow cytometry detected apoptotic cell shrinkage in cetuximab-treated DiFi cells, and significant apoptosis was confirmed by histology. TDSI-derived parameters quantified key morphological changes including cell size decreases during apoptosis in responsive tumors that occurred earlier than gross tumor volume regression. TDSI provides a unique measurement of apoptosis by identifying cellular characteristics, particularly cell shrinkage. The method will assist in understanding the underlying biology of solid tumors and predict tumor response to therapies. TDSI is free of any exogenous agent or radiation, and hence is very suitable to be incorporated into clinical applications.

Multidimensional diffusion MRI with spectrally modulated gradients reveals unprecedented microstructural detail

Characterization of porous media is essential in a wide range of biomedical and industrial applications. Microstructural features can be probed non-invasively by diffusion magnetic resonance imaging (dMRI). However, diffusion encoding in conventional dMRI may yield similar signatures for very different microstructures, which represents a significant limitation for disentangling individual microstructural features in heterogeneous materials. To solve this problem, we propose an augmented multidimensional diffusion encoding (MDE) framework, which unlocks a novel encoding dimension to assess time-dependent diffusion specific to structures with different microscopic anisotropies. Our approach relies on spectral analysis of complex but experimentally efficient MDE waveforms. Two independent contrasts to differentiate features such as cell shape and size can be generated directly by signal subtraction from only three types of measurements. Analytical calculations and simulations support our experimental observations. Proof-of-concept experiments were applied on samples with known and distinctly different microstructures. We further demonstrate substantially different contrasts in different tissue types of a post mortem brain. Our simultaneous assessment of restriction size and shape may be instrumental in studies of a wide range of porous materials, enable new insights into the microstructure of biological tissues or be of great value in diagnostics.

Diffusion MRI in the brain - Theory and concepts

Over the past two decades, diffusion MRI has become an essential tool in neuroimaging investigations. This is due to its sensitivity to the motion of water molecules as they diffuse through the microstructural environment, allowing diffusion MRI to be used as a 'probe' of tissue microstructure. Furthermore, this sensitivity is strongly direction-dependent, notably in brain white matter, due to the alignment of structures that restrict or hinder the motion of water molecules, notably axonal membranes. This provides a means of inferring the orientation of fibres in vivo, and by use of appropriate fibre-tracking algorithms, of delineating the path of white matter tracts in the brain. The ability to perform so-called tractography in humans in vivo non-invasively is unique to diffusion MRI, and is now used in applications such as neurosurgery planning and more broadly within investigations of brain connectomics. This review describes the theory and concepts of diffusion MRI and describes its most important areas of application in the brain, with a strong focus on tractography.

VERDICT MRI validation in fresh and fixed prostate specimens using patient-specific moulds for histological and MR alignment

The VERDICT framework for modelling diffusion MRI data aims to relate parameters from a biophysical model to histological features used for tumour grading in prostate cancer. Validation of the VERDICT model is necessary for clinical use. This study compared VERDICT parameters obtained ex vivo with histology in five specimens from radical prostatectomy. A patient-specific 3D-printed mould was used to investigate the effects of fixation on VERDICT parameters and to aid registration to histology. A rich diffusion data set was acquired in each ex vivo prostate before and after fixation. At both time points, data were best described by a two-compartment model: the model assumes that an anisotropic tensor compartment represents the extracellular space and a restricted sphere compartment models the intracellular space. The effect of fixation on model parameters associated with tissue microstructure was small. The patient-specific mould minimized tissue deformations and co-localized slices, so that rigid registration of MRI to histology images allowed region-based comparison with histology. The VERDICT estimate of the intracellular volume fraction corresponded to histological indicators of cellular fraction, including high values in tumour regions. The average sphere radius from VERDICT, representing the average cell size, was relatively uniform across samples. The primary diffusion direction from the extracellular compartment of the VERDICT model aligned with collagen fibre patterns in the stroma obtained by structure tensor analysis. This confirmed the biophysical relationship between ex vivo VERDICT parameters and tissue microstructure from histology.

Imaging brain microstructure with diffusion MRI: practicality and applications

This article gives an overview of microstructure imaging of the brain with diffusion MRI and reviews the state of the art. The microstructure-imaging paradigm aims to estimate and map microscopic properties of tissue using a model that links these properties to the voxel scale MR signal. Imaging techniques of this type are just starting to make the transition from the technical research domain to wide application in biomedical studies. We focus here on the practicalities of both implementing such techniques and using them in applications. Specifically, the article summarizes the relevant aspects of brain microanatomy and the range of diffusion-weighted MR measurements that provide sensitivity to them. It then reviews the evolution of mathematical and computational models that relate the diffusion MR signal to brain tissue microstructure, as well as the expanding areas of application. Next we focus on practicalities of designing a working microstructure imaging technique: model selection, experiment design, parameter estimation, validation, and the pipeline of development of this class of technique. The article concludes with some future perspectives on opportunities in this topic and expectations on how the field will evolve in the short-to-medium term.