Imaging Local Diffusive Dynamics Using Diffusion Exchange Spectroscopy MRI

Novel pore size-controlled, susceptibility matched, 3D-printed MRI phantoms

PURPOSE

We report the design concept and fabrication of MRI phantoms, containing blocks of aligned microcapillaires that can be stacked into larger arrays to construct diameter distribution phantoms or fractured, to create a "powder-averaged" emulsion of randomly oriented blocks for vetting or calibrating advanced MRI methods, that is, diffusion tensor imaging, AxCaliber MRI, MAP-MRI, and multiple pulsed field gradient or double diffusion-encoded microstructure imaging methods. The goal was to create a susceptibility-matched microscopically anisotropic but macroscopically isotropic phantom with a ground truth diameter that could be used to vet advanced diffusion methods for diameter determination in fibrous tissues.

METHODS

Two-photon polymerization, a novel three-dimensional printing method is used to fabricate blocks of capillaries. Double diffusion encoding methods were employed and analyzed to estimate the expected MRI diameter.

RESULTS

Susceptibility-matched microcapillary blocks or modules that can be assembled into large-scale MRI phantoms have been fabricated and measured using advanced diffusion methods, resulting in microscopic anisotropy and random orientation.

CONCLUSION

This phantom can vet and calibrate various advanced MRI methods and multiple pulsed field gradient or diffusion-encoded microstructure imaging methods. We demonstrated that two double diffusion encoding methods underestimated the ground truth diameter.

Nanostructure of hydrogenated amorphous silicon (a-Si:H) films studied by nuclear magnetic resonance

The aim of this work is to investigate the nanostructures of nanoporous materials by studying the anisotropy of the nuclear spin-spin and spin-lattice relaxations of the guest molecules trapped in the pores. The nuclear magnetic resonance (NMR) data are analyzed in the framework of the theory of the nuclear relaxation dominated by the dipole-dipole interactions in gas or liquid species contained in nanopores. A distinctive feature of this theory is the establishment of a relationship between the degree of orientation ordering of nanopores in the host matrix and their characteristic volume and the anisotropy of the NMR relaxation times. In this work the complex experimental and theoretical approach was applied to study the nanostructure of hydrogenated amorphous silicon (a-Si:H) films. A feature of this study is the simultaneous investigation of the three (T₁, T_1ρ, and T₂) NMR relaxation times, for the same sample. This allows us to determine not only the degree of orientation ordering of nanopores but also to estimate their size (∼1 nm) and correlation times of the nanopore fluctuations. The obtained results demonstrate that the developed approach is effective in studying details of nanostructure of different nanoporous materials.

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

Span of regularization for solution of inverse problems with application to magnetic resonance relaxometry of the brain

We present a new regularization method for the solution of the Fredholm integral equation (FIE) of the first kind, in which we incorporate solutions corresponding to a range of Tikhonov regularizers into the end result. This method identifies solutions within a much larger function space, spanned by this set of regularized solutions, than is available to conventional regularization methods. An additional key development is the use of dictionary functions derived from noise-corrupted inversion of the discretized FIE. In effect, we combine the stability of solutions with greater degrees of regularization with the resolution of those that are less regularized. The span of regularizations (SpanReg) method may be widely applicable throughout the field of inverse problems.

Stabilization of parameter estimates from multiexponential decay through extension into higher dimensions

Analysis of multiexponential decay has remained a topic of active research for over 200 years. This attests to the widespread importance of this problem and to the profound difficulties in characterizing the underlying monoexponential decays. Here, we demonstrate the fundamental improvement in stability and conditioning of this classic problem through extension to a second dimension; we present statistical analysis, Monte-Carlo simulations, and experimental magnetic resonance relaxometry data to support this remarkable fact. Our results are readily generalizable to higher dimensions and provide a potential means of circumventing conventional limits on multiexponential parameter estimation.

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.

Ultrafast methods for relaxation and diffusion

Relaxation and diffusion NMR measurements offer an approach to studying rotational and translational motion of molecules non-invasively, and they also provide chemical resolution complementary to NMR spectra. Multidimensional experiments enable the correlation of relaxation and diffusion parameters as well as the observation of molecular exchange phenomena through relaxation or diffusion contrast. This review describes how to accelerate multidimensional relaxation and diffusion measurements significantly through spatial encoding. This so-called ultrafast Laplace NMR approach shortens the experiment time to a fraction and makes even single-scan experiments possible. Single-scan experiments, in turn, significantly facilitate the use of nuclear spin hyperpolarization methods to boost sensitivity. The ultrafast Laplace NMR method is also applicable with low-field, mobile NMR instruments, and it can be exploited in many disciplines. For example, it has been used in studies of the dynamics of fluids in porous materials, identification of intra- and extracellular metabolites in cancer cells, and elucidation of aggregation phenomena in atmospheric surfactant solutions.

Data-Driven multi-Contrast spectral microstructure imaging with InSpect: INtegrated SPECTral component estimation and mapping

•

Unsupervised learning technique for spectroscopic analysis of quantitative MRI.

•

Shares information across voxels to improve estimation of multi-dimensional or single-dimensional spectra.

•

Spectral maps are dramatically improved compared to existing approaches.

•

Can potentially identify and map tissue environments; in placental diffusion-relaxometry MRI we demonstrate that it identifies components that correspond to distinct tissue types.

We introduce and demonstrate an unsupervised machine learning technique for spectroscopic analysis of quantitative MRI experiments. Our algorithm supports estimation of one-dimensional spectra from single-contrast data, and multidimensional correlation spectra from simultaneous multi-contrast data. These spectrum-based approaches allow model-free investigation of tissue properties, but require regularised inversion of a Laplace transform or Fredholm integral, which is an ill-posed calculation. Here we present a method that addresses this limitation in a data-driven way. The algorithm simultaneously estimates a canonical basis of spectral components and voxelwise maps of their weightings, thereby pooling information across whole images to regularise the ill-posed problem. We show in simulations that our algorithm substantially outperforms current voxelwise spectral approaches. We demonstrate the method on multi-contrast diffusion-relaxometry placental MRI scans, revealing anatomically-relevant sub-structures, and identifying dysfunctional placentas. Our algorithm vastly reduces the data required to reliably estimate spectra, opening up the possibility of quantitative MRI spectroscopy in a wide range of new applications. Our InSpect code is available at github.com/paddyslator/inspect.

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.

Multidimensional correlation MRI

Multidimensional correlation spectroscopy is emerging as a novel MRI modality that is well suited for microstructure and microdynamic imaging studies, especially of biological specimens. Conventional MRI methods only provide voxel-averaged and mostly macroscopically averaged information; these methods cannot disentangle intra-voxel heterogeneity on the basis of both water mobility and local chemical interactions. By correlating multiple MR contrast mechanisms and processing the data in an integrated manner, correlation spectroscopy is able to resolve the distribution of water populations according to their chemical and physical interactions with the environment. The use of a non-parametric, phenomenological representation of the multidimensional MR signal makes no assumptions about tissue structure, thereby allowing the study of microscopic structure and composition of complex heterogeneous biological systems. However, until recently, vast data requirements have confined these types of measurement to non-localized NMR applications and prevented them from being widely and successfully used in conjunction with imaging. Recent groundbreaking advancements have allowed this powerful NMR methodology to be migrated to MRI, initiating its emergence as a promising imaging approach. This review is not intended to cover the entire field of multidimensional MR; instead, it focuses on pioneering imaging applications and the challenges involved. In addition, the background and motivation that have led to multidimensional correlation MR development are discussed, along with the basic underlying mathematical concepts. The goal of the present work is to provide the reader with a fundamental understanding of the techniques developed and their potential benefits, and to provide guidance to help refine future applications of this technology.

A Tutorial Introduction to Inverse Problems in Magnetic Resonance

There has been a tremendous increase in applications of the inverse problem framework to parameter estimation in magnetic resonance. Attempting to capture both the basics of this formalism and modern developments would require an article of inordinate length. Therefore, in the following, we provide basic material as a practical introduction to the topic and an entree to the literature. First, we describe the formulation of linear and nonlinear inverse problems, with an emphasis on signal equations arising in magnetic resonance. We then describe the Fredholm equation of the first kind as a paradigm for these problems. This is followed by much more detailed considerations for determining solutions in the linear case, including central concepts such as condition number, regularization, and stability. Solution methods for nonlinear inverse problems are described next, followed by a treatment of their stability and regularization. Finally, we provide an introduction to compressed sensing, with signal reconstruction formulated as the solution to an inverse problem, making use of much of the previous material. Throughout, the emphasis is on outlines of the theory and on numerical examples, rather than on mathematical rigor and completeness.

The present and the future of microstructure MRI: From a paradigm shift to normal science

The aspiration of imaging tissue microstructure with MRI is to uncover micrometer-scale tissue features within millimeter-scale imaging voxels, in vivo. This kind of super-resolution has fueled a paradigm shift within the biomedical imaging community. However, what feels like an ongoing revolution in MRI, has been conceptually experienced in physics decades ago; from this point of view, our current developments can be seen as Thomas Kuhn's "normal science" stage of progress. While the concept of model-based quantification below the nominal imaging resolution is not new, its possibilities in neuroscience and neuroradiology are only beginning to be widely appreciated. This disconnect calls for communicating the progress of tissue microstructure MR imaging to its potential users. Here, a number of recent research developments are outlined in terms of the overarching concept of coarse-graining the tissue structure over an increasing diffusion length. A variety of diffusion models and phenomena are summarized on the phase diagram of diffusion MRI, with the unresolved problems and future directions corresponding to its unexplored domains.

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.

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.

Direct and specific assessment of axonal injury and spinal cord microenvironments using diffusion correlation imaging

We describe a practical two-dimensional (2D) diffusion MRI framework to deliver specificity and improve sensitivity to axonal injury in the spinal cord. This approach provides intravoxel distributions of correlations of water mobilities in orthogonal directions, revealing sub-voxel diffusion components. Here we use it to investigate water diffusivities along axial and radial orientations within spinal cord specimens with confirmed, tract-specific axonal injury. First, we show using transmission electron microscopy and immunohistochemistry that tract-specific axonal beading occurs following Wallerian degeneration in the cortico-spinal tract as direct sequelae to closed head injury. We demonstrate that although some voxel-averaged diffusion tensor imaging (DTI) metrics are sensitive to this axonal injury, they are non-specific, i.e., they do not reveal an underlying biophysical mechanism of injury. Then we employ 2D diffusion correlation imaging (DCI) to improve discrimination of different water microenvironments by measuring and mapping the joint water mobility distributions perpendicular and parallel to the spinal cord axis. We determine six distinct diffusion spectral components that differ according to their microscopic anisotropy and mobility. We show that at the injury site a highly anisotropic diffusion component completely disappears and instead becomes more isotropic. Based on these findings, an injury-specific MR image of the spinal cord was generated, and a radiological-pathological correlation with histological silver staining % area was performed. The resulting strong and significant correlation (r = 0.70, p < 0.0001) indicates the high specificity with which DCI detects injury-induced tissue alterations. We predict that the ability to selectively image microstructural changes following axonal injury in the spinal cord can be useful in clinical and research applications by enabling specific detection and increased sensitivity to injury-induced microstructural alterations. These results also encourage us to translate DCI to higher spatial dimensions to enable assessment of traumatic axonal injury, and possibly other diseases and disorders in the brain.

Ultrafast diffusion exchange nuclear magnetic resonance

The exchange of molecules between different physical or chemical environments due to diffusion or chemical transformations has a crucial role in a plethora of fundamental processes such as breathing, protein folding, chemical reactions and catalysis. Here, we introduce a method for a single-scan, ultrafast NMR analysis of molecular exchange based on the diffusion coefficient contrast. The method shortens the experiment time by one to four orders of magnitude. Consequently, it opens the way for high sensitivity quantification of important transient physical and chemical exchange processes such as in cellular metabolism. As a proof of principle, we demonstrate that the method reveals the structure of aggregates formed by surfactants relevant to aerosol research.

Analysis of exchange processes is time consuming by two-dimensional exchange NMR spectroscopy. Here the authors demonstrate a single-scan ultrafast Laplace NMR approach based on spatial encoding to measure molecular diffusion, with an increase by a factor six in the sensitivity per unit time.

Retaining information from multidimensional correlation MRI using a spectral regions of interest generator

Multidimensional correlation magnetic resonance imaging (MRI) is an emerging imaging modality that is capable of disentangling highly heterogeneous and opaque systems according to chemical and physical interactions of water within them. Using this approach, the conventional three dimensional MR scalar images are replaced with spatially resolved multidimensional spectra. The ensuing abundance in microstructural and chemical information is a blessing that incorporates a real challenge: how does one distill and refine it into images while retaining its significant components? In this paper we introduce a general framework that preserves the spectral information from spatially resolved multidimensional data. Equal weight is given to significant spectral components at the single voxel level, resulting in a summarized image spectrum. This spectrum is then used to define spectral regions of interest that are utilized to reconstruct images of sub-voxel components. Using numerical simulations we first show that, contrary to the conventional approach, the proposed framework preserves spectral resolution, and in turn, sensitivity and specificity of the reconstructed images. The retained spectral resolution allows, for the first time, to observe an array of distinct [Formula: see text]-[Formula: see text]-[Formula: see text] components images of the human brain. The robustly generated images of sub-voxel components overcome the limited spatial resolution of MRI, thus advancing multidimensional correlation MRI to fulfilling its full potential.

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.

Spatially-resolved 1H NMR relaxation-exchange measurements in heterogeneous media

In the last decades, the ¹H NMR T₂-T₂ relaxation-exchange (REXSY) technique has become an essential tool for the molecular investigation of simple and complex fluids in heterogeneous porous solids and soft matter, where the mixing-time-evolution of cross-correlated T₂-T₂ peaks enables a quantitative study of diffusive exchange kinetics in multi-component systems. Here, we present a spatially-resolved implementation of the T₂-T₂ correlation technique, named z-T₂-T₂, based on one-dimensional spatial mapping along z using a rapid frequency-encode imaging scheme. Compared to other phase-encoding methods, the adopted MRI technique has two distinct advantages: (i) is has the same experimental duration of a standard (bulk) T₂-T₂ measurement, and (ii) it provides a high spatial resolution. The proposed z-T₂-T₂ method is first validated against bulk T₂-T₂ measurements on homogeneous phantom consisting of cyclohexane uniformly imbibed in finely-sized α-Al₂O₃ particles at a spatial resolution of 0.47 mm; thereafter, its performance is demonstrated, on a layered bed of multi-sized α-Al₂O₃ particles, for revealing spatially-dependent molecular exchange kinetics properties of intra- and inter-particle cyclohexane as a function of particle size. It is found that localised z-T₂-T₂ spectra provide well resolved cross peaks whilst such resolution is lost in standard bulk T₂-T₂ data. Future prospective applications of the method lie, in particular, in the local characterisation of mass transport phenomena in multi-component porous media, such as rock cores and heterogeneous catalysts.

Improved stimulated echo in diffusion magnetic resonance imaging: introducing a π pulse for SNR enhancement

PURPOSE

Anomalous diffusion in biological tissues can be examined by diffusion MRI for various applications, including tumor diagnosis and measurement of brain fiber pathways. However, the measurement of anomalous diffusion requires high b-values for the diffusion gradient in MRI, and current MRI methods cannot provide a high SNR. This study aimed to improve on the standard stimulated echo (STE) to enhance the SNR in diffusion MRI with high b-values.

METHODS

Because of hardware limitations and human safety considerations, prolonging the diffusion time (Δ) is 1 of the few methods available to realize high b-values. Here, we propose a new echo mechanism for diffusion MRI to enhance SNRs under long Δ. By introducing a π pulse at the midpoint between 2nd and 3rd π/2 pulses of STE, we refocus the magnetic moment vectors in the longitudinal plane before the third π/2 pulse is applied, which preserves the full echo signals. This sequence was compared with STE and spin echo (SE). Nine Δs were tested in a phantom. Multi b-values with 2 Δs were tested in a mouse liver, brain, and tumor.

RESULTS

Compared with STE and SE, the proposed improved STE (ISTE) exhibited an improved SNR in the phantom experiment and improved performance in the in vivo experiments.

CONCLUSION

By using the proposed echo mechanism in diffusion MRI, we enhanced the SNR of the images, which enables us to investigate diffusion behavior at higher b-values and further facilitates the development of quantitative diffusion MRI and radiomics.

Generalized Mean Apparent Propagator MRI to Measure and Image Advective and Dispersive Flows in Medicine and Biology

Water transport in biological systems spans different regimes with distinct physical behaviors: diffusion, advection, and dispersion. Identifying these regimes is of paramount importance in many in vivo applications, among them, measuring microcirculation of blood in capillary networks and cerebrospinal fluid transport in the glymphatic system. Diffusion magnetic resonance imaging (dMRI) can be used to encode water displacements, and a Fourier transform of the acquired signal furnishes a displacement probability density function known as the propagator. This transformation normally requires the use of a fast Fourier transform (FFT), which presents major feasibility challenges when scanning in vivo, mainly because of dense signal sampling, resulting in long acquisition times. A second approach to reconstruct the propagator is by using analytical representation of the signal, overcoming many of the FFT's limitations. In all analytical implementations of dMRI to date, the translational motion of water has been assumed to be exclusively diffusive, which is the case only in the absence of flow. However, retaining the phase information from the diffusion signal provides the ability to measure both mean coherent velocity and random diffusion from a single experiment. We implement and extend an analytical framework, mean apparent propagator (MAP), which can account for non-zero flow conditions. We call this method generalized MAP or GMAP. We describe a numerical optimization scheme and implement it on data from an MRI flow phantom constructed from a pack of 10- [Formula: see text] beads. The advantages of GMAP over the FFT-based method in the context of sampling density and low-flow detection were demonstrated, and analytically derived propagator moments were shown to agree with theoretical values even after data subsampling. GMAP would enable the detection of microflow in vivo that could help elucidate many important biological processes.

Water mobility spectral imaging of the spinal cord: Parametrization of model-free Laplace MRI

Diffusion magnetic resonance imaging (dMRI) of biological systems most often results in non-monoexponential signal, due to their complexity and heterogeneity. One approach to interpreting dMRI data without imposing tissue microstructural models is to fit the signal to a multiexponential function, which is sometimes referred to as an inverse Laplace transformation, and to display the coefficients as a distribution of the diffusivities, or water mobility spectra. Until recently, this method has not been used in a voxelwise manner, mainly because of heavy data requirements. With recent advancements in processing and experimental design, voxelwise Laplace MRI approaches are becoming feasible and attractive. The rich spectral information, combined with a three-dimensional image, presents a challenge because it tremendously increases the dimensionality of the data and requires a robust method for interpretation and analysis. In this work, we suggest parameterizing the empirically measured water mobility spectra using a bimodal lognormal function. This approach allows for a compact representation of the spectrum, and it also resolves overlapping spectral peaks, which allows for a robust extraction of their signal fraction. We apply the method on a fixed spinal cord sample and use it to generate robust intensity images of slow- and fast-diffusion components. Using the parametric variables, we create novel image contrasts, among them the information entropy of the water mobility spectrum, which pack unique features of the individual diffusion regimes in the investigated system.

Rapid detection of the presence of diffusion exchange

Diffusion exchange spectroscopy (DEXSY) provides a detailed picture of how fluids in different microenvironments communicate with one another but requires a large amount of data. For DEXSY MRI, a simple measure of apparent exchanging fractions may suffice to characterize and differentiate materials and tissues. Reparameterizing signal intensity from a PGSE-storage-PGSE experiment as a function of the sum, b_s=b₁+b₂, and difference b_d=b₂-b₁ of the diffusion encodings separates diffusion weighting from exchange weighting. Exchange leads to upward curvature along a slice of constant b_s. Exchanging fractions can be measured rapidly by a finite difference approximation of the curvature using four data points. The method is generalized for non-steady-state and multi-site exchange. We apply the method to image exchanging fractions and calculate exchange rates of water diffusing across the bulk water interface of a glass capillary array.

A novel MRI phantom to study interstitial fluid transport in the glymphatic system

The glymphatic system is a recently discovered transport system, mediated by cerebral spinal fluid (CSF), that clears metabolic and cellular waste products in the brain. This system's function in the brain is analogous to that of the lymphatic system in the rest of the mammalian body. It is hypothesized that CSF clears harmful chemicals from the brain by flowing through interstitial spaces in the brain during sleep. While there is growing recognition of the critical role the glymphatic system plays in maintaining normal brain health and in explaining pathology, there are few noninvasive imaging methods that measure and characterize the efficacy of glymphatic transport in vivo. In this study we designed, constructed, and tested a glymphatic transport magnetic resonance imaging (MRI) flow phantom, which combines regions that mimic CSF-filled ventricles and brain interstitial space. We tested high- and low-q space diffusion MRI and diffusion tensor imaging (DTI) acquisitions to determine if they could detect, measure, and map interstitial glymphatic flows. The results suggest that, under certain flow conditions, diffusion-weighted MRI can detect the enhanced mixing that occurs during glymphatic clearance.

CEST, ASL, and magnetization transfer contrast: How similar pulse sequences detect different phenomena

Chemical exchange saturation transfer (CEST), arterial spin labeling (ASL), and magnetization transfer contrast (MTC) methods generate different contrasts for MRI. However, they share many similarities in terms of pulse sequences and mechanistic principles. They all use RF pulse preparation schemes to label the longitudinal magnetization of certain proton pools and follow the delivery and transfer of this magnetic label to a water proton pool in a tissue region of interest, where it accumulates and can be detected using any imaging sequence. Due to the versatility of MRI, differences in spectral, spatial or motional selectivity of these schemes can be exploited to achieve pool specificity, such as for arterial water protons in ASL, protons on solute molecules in CEST, and protons on semi-solid cell structures in MTC. Timing of these sequences can be used to optimize for the rate of a particular delivery and/or exchange transfer process, for instance, between different tissue compartments (ASL) or between tissue molecules (CEST/MTC). In this review, magnetic labeling strategies for ASL and the corresponding CEST and MTC pulse sequences are compared, including continuous labeling, single-pulse labeling, and multi-pulse labeling. Insight into the similarities and differences among these techniques is important not only to comprehend the mechanisms and confounds of the contrasts they generate, but also to stimulate the development of new MRI techniques to improve these contrasts or to reduce their interference. This, in turn, should benefit many possible applications in the fields of physiological and molecular imaging and spectroscopy.

Brain active transmembrane water cycling measured by MR is associated with neuronal activity

PURPOSE

fMRI is widely used to study brain activity. Unfortunately, conventional fMRI methods assess neuronal activity only indirectly, through hemodynamic coupling. Here, we show that active, steady-state transmembrane water cycling (AWC) could serve as a basis for a potential fMRI mechanism for direct neuronal activity detection.

METHODS

AWC and neuronal actitivity in rat organotypic cortical cultures were simultaneously measured with a hybrid MR-fluorescence system. Perfusion with a paramagnetic MRI contrast agent, Gadoteridol, allows NMR determination of the kinetics of transcytolemmal water exchange. Changes in intracellular calcium concentration, [Ca_i ²⁺ ] were used as a proxy of neuronal activity and were monitored by fluorescence imaging.

RESULTS

When we alter neuronal activity by titrating with extracellular [K⁺ ] near the normal value, we see an AWC response resembling Na⁺ -K⁺ -ATPase (NKA) Michaelis-Menten behavior. When we treat with the voltage-gated sodium channel inhibitor, or with an excitatory postsynaptic inhibitor cocktail, we see AWC decrease by up to 71%. AWC was found also to be positively correlated with the basal level of spontaneous activity, which varies in different cultures.

CONCLUSIONS

These results suggest that AWC is associated with neuronal activity and NKA activity is a major contributor in coupling AWC to neuronal activity. Although AWC comprises steady-state, homeostatic transmembrane water exchange, our analysis also yields a simultaneous measure of the average cell volume, which reports any slower net transmembrane water transport.

Improving neurophysiological biomarkers for functional myoclonic movements

Introduction

Differentiating between functional jerks (FJ) and organic myoclonus can be challenging. At present, the only advanced diagnostic biomarker to support FJ is the Bereitschaftspotential (BP). However, its sensitivity is limited and its evaluation subjective. Recently, event related desynchronisation in the broad beta range (13–45 Hz) prior to functional generalised axial (propriospinal) myoclonus was reported as a possible complementary diagnostic marker for FJ. Here we study the value of ERD together with a quantified BP in clinical practice.

Methods

Twenty-nine patients with FJ and 16 patients with cortical myoclonus (CM) were included. Jerk-locked back-averaging for determination of the ‘classical’ and quantified BP, and time-frequency decomposition for the event related desynchronisation (ERD) were performed. Diagnostic gain, sensitivity and specificity were obtained for individual and combined techniques.

Results

We detected a classical BP in 14/29, a quantitative BP in 15/29 and an ERD in 18/29 patients. At group level we demonstrate that ERD in the broad beta band preceding a jerk has significantly higher amplitude in FJ compared to CM (respectively −0.14 ± 0.13 and +0.04 ± 0.09 (p < 0.001)). Adding ERD to the classical BP achieved an additional diagnostic gain of 53%. Furthermore, when combining ERD with quantified and classical BP, an additional diagnostic gain of 71% was achieved without loss of specificity.

Conclusion

Based on the current findings we propose to the use of combined beta ERD assessment and quantitative BP analyses in patients with a clinical suspicion for all types of FJ with a negative classical BP.

•

Differentiating between functional jerks and organic myoclonus is often difficult.

•

The sensitivity of a positive bereitschaftspotential for diagnosing FJ is low.

•

FJ are preceded by event related desynchronisation (ERD) in the beta band.

•

Combining beta ERD and BP improves diagnosing FJ.

•

The ERD method is of special relevance in possible FJ patients with a negative BP.

Magnetic resonance microdynamic imaging reveals distinct tissue microenvironments

Magnetic resonance imaging (MRI) provides a powerful set of tools with which to investigate biological tissues noninvasively and in vivo. Tissues are heterogeneous in nature; an imaging voxel contains an ensemble of different cells and extracellular matrix components. A long-standing challenge has been to infer the content of and interactions among these microscopic tissue components within a macroscopic imaging voxel. Spatially resolved multidimensional relaxation-diffusion correlation (REDCO) spectroscopy holds the potential to deliver such microdynamic information. However, to date, vast data requirements have mostly relegated these type of measurements to nuclear magnetic resonance applications and prevented them from being widely and successfully used in conjunction with imaging. By using a novel data acquisition and processing strategy in this study, spatially resolved REDCO could be performed in reasonable scanning times with excellent prospects for clinical applications. This new MR imaging framework-which we term "magnetic resonance microdynamic imaging (MRMI)"-permits the simultaneous noninvasive and model-free quantification of multiple subcellular, cellular, and interstitial tissue microenvironments within a voxel. MRMI is demonstrated with a fixed spinal cord specimen, enabling the quantification of microscopic tissue components with unprecedented specificity. Tissue components, such as axons, neuronal and glial soma, and myelin were identified on the basis of their multispectral signature within individual imaging voxels. These tissue elements could then be composed into images and be correlated with immunohistochemistry findings. MRMI provides novel image contrasts of tissue components and a new family of microdynamic biomarkers that could lead to new diagnostic imaging approaches to probe biological tissue alterations accompanied by pathological or developmental changes.