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.

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.

A review on investigation of the basic contrast mechanism underlying multidimensional diffusion MRI in assessment of neurological disorders

INTRODUCTION

Multidimensional diffusion MRI (MDD MRI) is a novel diffusion technique that uses advanced gradient waveforms for microstructural tissue characterization to provide information about average rate, anisotropy and orientation of the diffusion and to disentangle the signal fraction from specific cell types i.e., elongated cells, isotropic cells and free water.

AIM

To review the diagnostic potential of MDD MRI in the clinical setting for microstructural tissue characterization in patients with neurological disorders to aid in patient care and treatment.

METHOD

A scoping review on the clinical applications of MDD MRI was conducted from original articles published in PubMed and Scopus from 2015 to 2021 using the keywords "Multidimensional diffusion MRI" OR "diffusion tensor distribution" OR "Tensor-Valued Diffusion" OR "b-tensor encoding" OR "microscopic diffusion anisotropy" OR "microscopic anisotropy" OR "microscopic fractional anisotropy" OR "double diffusion encoding" OR "triple diffusion encoding" OR "double pulsed field gradients" OR "double wave vector" OR "correlation tensor imaging" AND "brain" OR "axons".

RESULTS

Initially 145 articles were screened and after applying inclusion and exclusion criteria, nine articles were included in the final analysis. In most of these studies, microscopic diffusion anisotropy within the lesion showed deviation from the normal-appearing tissue.

CONCLUSION

Multidimensional diffusion MRI can provide better quantification and visualization of tissue microstructure than conventional diffusion MRI and can be used in the clinical setting for diagnosis of neurological disorders.

Glioma grading, molecular feature classification, and microstructural characterization using MR diffusional variance decomposition (DIVIDE) imaging

OBJECTIVE

To evaluate the potential of diffusional variance decomposition (DIVIDE) for grading, molecular feature classification, and microstructural characterization of gliomas.

MATERIALS AND METHODS

Participants with suspected gliomas underwent DIVIDE imaging, yielding parameter maps of fractional anisotropy (FA), mean diffusivity (MD), anisotropic mean kurtosis (MK_A), isotropic mean kurtosis (MK_I), total mean kurtosis (MK_T), MK_A/MK_T, and microscopic fractional anisotropy (μFA). Tumor type and grade, isocitrate dehydrogenase (IDH) 1/2 mutant status, and the Ki-67 labeling index (Ki-67 LI) were determined after surgery. Statistical analysis included 33 high-grade gliomas (HGG) and 17 low-grade gliomas (LGG). Tumor diffusion metrics were compared between HGG and LGG, among grades, and between wild and mutated IDH types using appropriate tests according to normality assessment results. Receiver operating characteristic and Spearman correlation analysis were also used for statistical evaluations.

RESULTS

FA, MD, MK_A, MK_I, MK_T, μFA, and MK_A/MK_T differed between HGG and LGG (FA: p = 0.047; MD: p = 0.037, others p < 0.001), and among glioma grade II, III, and IV (FA: p = 0.048; MD: p = 0.038, others p < 0.001). All diffusion metrics differed between wild-type and mutated IDH tumors (MK_I: p = 0.003; others: p < 0.001). The metrics that best discriminated between HGG and LGGs and between wild-type and mutated IDH tumors were MK_T and FA respectively (area under the curve 0.866 and 0.881). All diffusion metrics except FA showed significant correlation with Ki-67 LI, and MK_I had the highest correlation coefficient (r_s = 0.618).

CONCLUSION

DIVIDE is a promising technique for glioma characterization and diagnosis.

KEY POINTS

• DIVIDE metrics MK_I is related to cell density heterogeneity while MK_A and μFA are related to cell eccentricity. • DIVIDE metrics can effectively differentiate LGG from HGG and IDH mutation from wild-type tumor, and showed significant correlation with the Ki-67 labeling index. • MK_I was larger than MK_A which indicates predominant cell density heterogeneity in gliomas. • MK_A and MK_I increased with grade or degree of malignancy, however with a relatively larger increase in the cell eccentricity metric MK_A in relation to the cell density heterogeneity metric MK_I.

Disentangling white-matter damage from physiological fibre orientation dispersion in multiple sclerosis

Multiple sclerosis leads to diffuse damage of the central nervous system, affecting also the normal-appearing white matter. Demyelination and axonal degeneration reduce regional fractional anisotropy in normal-appearing white matter, which can be routinely mapped with diffusion tensor imaging. However, the standard fractional anisotropy metric is also sensitive to physiological variations in orientation dispersion of white matter fibres. This complicates the detection of disease-related damage in large parts of cerebral white matter where microstructure physiologically displays a high degree of fibre dispersion. To resolve this ambiguity, we employed a novel tensor-valued encoding method for diffusion MRI, which yields a microscopic fractional anisotropy metric that is unaffected by regional variations in orientation dispersion. In 26 patients with relapsing-remitting multiple sclerosis, 14 patients with primary-progressive multiple sclerosis and 27 age-matched healthy controls, we compared standard fractional anisotropy mapping with the novel microscopic fractional anisotropy mapping method, focusing on normal-appearing white matter. Mean microscopic fractional anisotropy and standard fractional anisotropy of normal-appearing white matter were significantly reduced in both patient groups relative to healthy controls, but microscopic fractional anisotropy yielded a better reflection of disease-related white-matter alterations. The reduction in mean microscopic fractional anisotropy showed a significant positive linear relationship with physical disability, as reflected by the expanded disability status scale. Mean reduction of microscopic fractional anisotropy in normal-appearing white matter also scaled positively with individual cognitive dysfunction, as measured with the symbol digit modality test. Mean microscopic fractional anisotropy reduction in normal-appearing white matter also showed a positive relationship with total white-matter lesion load as well as lesion load in specific tract systems. None of these relationships between normal-appearing white-matter microstructure and clinical, cognitive or structural measures emerged when using mean fractional anisotropy. Together, the results provide converging evidence that microscopic fractional anisotropy mapping substantially advances the assessment of cerebral white matter in multiple sclerosis by disentangling microstructure damage from variations in physiological fibre orientation dispersion at the stage of data acquisition. Since tensor-valued encoding can be implemented in routine diffusion MRI, microscopic fractional anisotropy mapping bears considerable potential for the future assessment of disease progression in normal-appearing white matter in both relapsing-remitting and progressive forms of multiple sclerosis as well as other white-matter-related brain diseases.

Time-dependent diffusion in undulating thin fibers: Impact on axon diameter estimation

Diffusion MRI may enable non-invasive mapping of axonal microstructure. Most approaches infer axon diameters from effects of time-dependent diffusion on the diffusion-weighted MR signal by modeling axons as straight cylinders. Axons do not, however, propagate in straight trajectories, and so far the impact of the axonal trajectory on diameter estimation has been insufficiently investigated. Here, we employ a toy model of axons, which we refer to as the undulating thin fiber model, to analyze the impact of undulating trajectories on the time dependence of diffusion. We study time-dependent diffusion in the frequency domain and characterize the diffusion spectrum by its height, width, and low-frequency behavior (power law exponent). Results show that microscopic orientation dispersion of the thin fibers is the main parameter that determines the characteristics of the diffusion spectra. At lower frequencies (longer diffusion times), straight cylinders and undulating thin fibers can have virtually identical spectra. If the straight-cylinder assumption is used to interpret data from undulating thin axons, the diameter is overestimated by an amount proportional to the undulation amplitude and microscopic orientation dispersion of the fibers. At higher frequencies (shorter diffusion times), spectra from cylinders and undulating thin fibers differ. The low-frequency behavior of the spectra from the undulating thin fibers may also differ from that of cylinders, because the power law exponent of undulating fibers can reach values below 2 for experimentally relevant frequency ranges. In conclusion, we argue that the non-straight nature of axonal trajectories should not be overlooked when analyzing and interpreting diffusion MRI data.

Motion-compensated b-tensor encoding for in vivo cardiac diffusion-weighted imaging

Motion is a major confound in diffusion-weighted imaging (DWI) in the body, and it is a common cause of image artefacts. The effects are particularly severe in cardiac applications, due to the nonrigid cyclical deformation of the myocardium. Spin echo-based DWI commonly employs gradient moment-nulling techniques to desensitise the acquisition to velocity and acceleration, ie, nulling gradient moments up to the 2nd order (M2-nulled). However, current M2-nulled DWI scans are limited to encode diffusion along a single direction at a time. We propose a method for designing b-tensors of arbitrary shapes, including planar, spherical, prolate and oblate tensors, while nulling gradient moments up to the 2nd order and beyond. The design strategy comprises initialising the diffusion encoding gradients in two encoding blocks about the refocusing pulse, followed by appropriate scaling and rotation, which further enables nulling undesired effects of concomitant gradients. Proof-of-concept assessment of in vivo mean diffusivity (MD) was performed using linear and spherical tensor encoding (LTE and STE, respectively) in the hearts of five healthy volunteers. The results of the M2-nulled STE showed that (a) the sequence was robust to cardiac motion, and (b) MD was higher than that acquired using standard M2-nulled LTE, where diffusion-weighting was applied in three orthogonal directions, which may be attributed to the presence of restricted diffusion and microscopic diffusion anisotropy. Provided adequate signal-to-noise ratio, STE could significantly shorten estimation of MD compared with the conventional LTE approach. Importantly, our theoretical analysis and the proposed gradient waveform design may be useful in microstructure imaging beyond diffusion tensor imaging where the effects of motion must be suppressed.

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.

Cumulant expansions for measuring water exchange using diffusion MRI

The rate of water exchange across cell membranes is a parameter of biological interest and can be measured by diffusion magnetic resonance imaging (dMRI). In this work, we investigate a stochastic model for the diffusion-and-exchange of water molecules. This model provides a general solution for the temporal evolution of dMRI signal using any type of gradient waveform, thereby generalizing the signal expressions for the Kärger model. Moreover, we also derive a general nth order cumulant expansion of the dMRI signal accounting for water exchange, which has not been explored in earlier studies. Based on this analytical expression, we compute the cumulant expansion for dMRI signals for the special case of single diffusion encoding (SDE) and double diffusion encoding (DDE) sequences. Our results provide a theoretical guideline on optimizing experimental parameters for SDE and DDE sequences, respectively. Moreover, we show that DDE signals are more sensitive to water exchange at short-time scale but provide less attenuation at long-time scale than SDE signals. Our theoretical analysis is also validated using Monte Carlo simulations on synthetic structures.

Quantification of the Intracellular Life Time of Water Molecules to Measure Transport Rates of Human Aquaglyceroporins

Orthodox aquaporins are transmembrane channel proteins that facilitate rapid diffusion of water, while aquaglyceroporins facilitate the diffusion of small uncharged molecules such as glycerol and arsenic trioxide. Aquaglyceroporins play important roles in human physiology, in particular for glycerol metabolism and arsenic detoxification. We have developed a unique system applying the strain of the yeast Pichia pastoris, where the endogenous aquaporins/aquaglyceroporins have been removed and human aquaglyceroporins AQP3, AQP7, and AQP9 are recombinantly expressed enabling comparative permeability measurements between the expressed proteins. Using a newly established Nuclear Magnetic Resonance approach based on measurement of the intracellular life time of water, we propose that human aquaglyceroporins are poor facilitators of water and that the water transport efficiency is similar to that of passive diffusion across native cell membranes. This is distinctly different from glycerol and arsenic trioxide, where high glycerol transport efficiency was recorded.

Resolution limit of cylinder diameter estimation by diffusion MRI: The impact of gradient waveform and orientation dispersion

Diffusion MRI has been proposed as a non-invasive technique for axonal diameter mapping. However, accurate estimation of small diameters requires strong gradients, which is a challenge for the transition of the technique from preclinical to clinical MRI scanners, since these have weaker gradients. In this work, we develop a framework to estimate the lower bound for accurate diameter estimation, which we refer to as the resolution limit. We analyse only the contribution from the intra-axonal space and assume that axons can be represented by impermeable cylinders. To address the growing interest in using techniques for diffusion encoding that go beyond the conventional single diffusion encoding (SDE) sequence, we present a generalised analysis capable of predicting the resolution limit regardless of the gradient waveform. Using this framework, waveforms were optimised to minimise the resolution limit. The results show that, for parallel cylinders, the SDE experiment is optimal in terms of yielding the lowest possible resolution limit. In the presence of orientation dispersion, diffusion encoding sequences with square-wave oscillating gradients were optimal. The resolution limit for standard clinical MRI scanners (maximum gradient strength 60-80 mT/m) was found to be between 4 and 8 μm, depending on the noise levels and the level of orientation dispersion. For scanners with a maximum gradient strength of 300 mT/m, the limit was reduced to between 2 and 5 μm.

NMR quantification of diffusional exchange in cell suspensions with relaxation rate differences between intra and extracellular compartments

Water transport across cell membranes can be measured non-invasively with diffusion NMR. We present a method to quantify the intracellular lifetime of water in cell suspensions with short transverse relaxation times, T₂, and also circumvent the confounding effect of different T₂ values in the intra- and extracellular compartments. Filter exchange spectroscopy (FEXSY) is specifically sensitive to exchange between compartments with different apparent diffusivities. Our investigation shows that FEXSY could yield significantly biased results if differences in T₂ are not accounted for. To mitigate this problem, we propose combining FEXSY with diffusion-relaxation correlation experiment, which can quantify differences in T₂ values in compartments with different diffusivities. Our analysis uses a joint constrained fitting of the two datasets and considers the effects of diffusion, relaxation and exchange in both experiments. The method is demonstrated on yeast cells with and without human aquaporins.

Quantification of microcirculatory parameters by joint analysis of flow-compensated and non-flow-compensated intravoxel incoherent motion (IVIM) data

The aim of this study was to improve the accuracy and precision of perfusion fraction and blood velocity dispersion estimates in intravoxel incoherent motion (IVIM) imaging, using joint analysis of flow-compensated and non-flow-compensated motion-encoded MRI data. A double diffusion encoding sequence capable of switching between flow-compensated and non-flow-compensated encoding modes was implemented. In vivo brain data were collected in eight healthy volunteers and processed using the joint analysis. Simulations were used to compare the performance of the proposed analysis method with conventional IVIM analysis. With flow compensation, strong rephasing was observed for the in vivo data, approximately cancelling the IVIM effect. The joint analysis yielded physiologically reasonable perfusion fraction maps. Estimated perfusion fractions were 2.43 ± 0.81% in gray matter, 1.81 ± 0.90% in deep gray matter, and 1.64 ± 0.72% in white matter (mean ± SD, n = 8). Simulations showed improved accuracy and precision when using joint analysis of flow-compensated and non-flow-compensated data, compared with conventional IVIM analysis. Double diffusion encoding with flow compensation was feasible for in vivo imaging of the perfusion fraction in the brain. The strong rephasing implied that blood flowing through the cerebral microvascular system was closer to the ballistic limit than the diffusive limit.

Apparent exchange rate for breast cancer characterization

Although diffusion MRI has shown promise for the characterization of breast cancer, it has low specificity to malignant subtypes. Higher specificity might be achieved if the effects of cell morphology and molecular exchange across cell membranes could be disentangled. The quantification of exchange might thus allow the differentiation of different types of breast cancer cells. Based on differences in diffusion rates between the intra- and extracellular compartments, filter exchange spectroscopy/imaging (FEXSY/FEXI) provides non-invasive quantification of the apparent exchange rate (AXR) of water between the two compartments. To test the feasibility of FEXSY for the differentiation of different breast cancer cells, we performed experiments on several breast epithelial cell lines in vitro. Furthermore, we performed the first in vivo FEXI measurement of water exchange in human breast. In cell suspensions, pulsed gradient spin-echo experiments with large b values and variable pulse duration allow the characterization of the intracellular compartment, whereas FEXSY provides a quantification of AXR. These experiments are very sensitive to the physiological state of cells and can be used to establish reliable protocols for the culture and harvesting of cells. Our results suggest that different breast cancer subtypes can be distinguished on the basis of their AXR values in cell suspensions. Time-resolved measurements allow the monitoring of the physiological state of cells in suspensions over the time-scale of hours, and reveal an abrupt disintegration of the intracellular compartment. In vivo, exchange can be detected in a tumor, whereas, in normal tissue, the exchange rate is outside the range experimentally accessible for FEXI. At present, low signal-to-noise ratio and limited scan time allows the quantification of AXR only in a region of interest of relatively large tumors.

NMR diffusion-encoding with axial symmetry and variable anisotropy: Distinguishing between prolate and oblate microscopic diffusion tensors with unknown orientation distribution

We introduce a nuclear magnetic resonance method for quantifying the shape of axially symmetric microscopic diffusion tensors in terms of a new diffusion anisotropy metric, DΔ, which has unique values for oblate, spherical, and prolate tensor shapes. The pulse sequence includes a series of equal-amplitude magnetic field gradient pulse pairs, the directions of which are tailored to give an axially symmetric diffusion-encoding tensor b with variable anisotropy bΔ. Averaging of data acquired for a range of orientations of the symmetry axis of the tensor b renders the method insensitive to the orientation distribution function of the microscopic diffusion tensors. Proof-of-principle experiments are performed on water in polydomain lyotropic liquid crystals with geometries that give rise to microscopic diffusion tensors with oblate, spherical, and prolate shapes. The method could be useful for characterizing the geometry of fluid-filled compartments in porous solids, soft matter, and biological tissues.

Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: applications in healthy volunteers and in brain tumors

The anisotropy of water diffusion in brain tissue is affected by both disease and development. This change can be detected using diffusion MRI and is often quantified by the fractional anisotropy (FA) derived from diffusion tensor imaging (DTI). Although FA is sensitive to anisotropic cell structures, such as axons, it is also sensitive to their orientation dispersion. This is a major limitation to the use of FA as a biomarker for "tissue integrity", especially in regions of complex microarchitecture. In this work, we seek to circumvent this limitation by disentangling the effects of microscopic diffusion anisotropy from the orientation dispersion. The microscopic fractional anisotropy (μFA) and the order parameter (OP) were calculated from the contrast between signal prepared with directional and isotropic diffusion encoding, where the latter was achieved by magic angle spinning of the q-vector (qMAS). These parameters were quantified in healthy volunteers and in two patients; one patient with meningioma and one with glioblastoma. Finally, we used simulations to elucidate the relation between FA and μFA in various micro-architectures. Generally, μFA was high in the white matter and low in the gray matter. In the white matter, the largest differences between μFA and FA were found in crossing white matter and in interfaces between large white matter tracts, where μFA was high while FA was low. Both tumor types exhibited a low FA, in contrast to the μFA which was high in the meningioma and low in the glioblastoma, indicating that the meningioma contained disordered anisotropic structures, while the glioblastoma did not. This interpretation was confirmed by histological examination. We conclude that FA from DTI reflects both the amount of diffusion anisotropy and orientation dispersion. We suggest that the μFA and OP may complement FA by independently quantifying the microscopic anisotropy and the level of orientation coherence.

Noninvasive mapping of water diffusional exchange in the human brain using filter-exchange imaging

We present the first in vivo application of the filter-exchange imaging protocol for diffusion MRI. The protocol allows noninvasive mapping of the rate of water exchange between microenvironments with different self-diffusivities, such as the intracellular and extracellular spaces in tissue. Since diffusional water exchange across the cell membrane is a fundamental process in human physiology and pathophysiology, clinically feasible and noninvasive imaging of the water exchange rate would offer new means to diagnose disease and monitor treatment response in conditions such as cancer and edema. The in vivo use of filter-exchange imaging was demonstrated by studying the brain of five healthy volunteers and one intracranial tumor (meningioma). Apparent exchange rates in white matter range from 0.8±0.08 s(-1) in the internal capsule, to 1.6±0.11 s(-1) for frontal white matter, indicating that low values are associated with high myelination. Solid tumor displayed values of up to 2.9±0.8 s(-1). In white matter, the apparent exchange rate values suggest intra-axonal exchange times in the order of seconds, confirming the slow exchange assumption in the analysis of diffusion MRI data. We propose that filter-exchange imaging could be used clinically to map the water exchange rate in pathologies. Filter-exchange imaging may also be valuable for evaluating novel therapies targeting the function of aquaporins.

Apparent exchange rate mapping with diffusion MRI

Water exchange through the cell membranes is an important feature of cells and tissues. The rate of exchange is determined by factors such as membrane lipid composition and organization, as well as the type and activity of aquaporins. A method for noninvasively estimating the rate of water exchange would be useful for characterizing pathological conditions, e.g., tumors, multiple sclerosis, and ischemic stroke, expected to be associated with a change of the membrane barrier properties. This study describes the filter exchange imaging method for determining the rate of water exchange between sites having different apparent diffusion coefficients. The method is based on the filter-exchange pulsed gradient spin-echo NMR spectroscopy experiment, which is here modified to be compatible with the constraints of clinical MR scanners. The data is analyzed using a model-free approach yielding maps of the apparent exchange rate, here being introduced in analogy with the concept of the apparent diffusion coefficient. Proof-of-principle experiments are performed on microimaging and whole-body clinical scanners using yeast suspension phantoms. The limitations and appropriate experimental conditions are examined. The results demonstrate that filter exchange imaging is a fast and reliable method for characterizing exchange, and that it has the potential to become a powerful diagnostic tool.