Impact of transcytolemmal water exchange on estimates of tissue microstructural properties derived from diffusion MRI

Restriction-induced time-dependent transcytolemmal water exchange: Revisiting the Kӓrger exchange model

The Kӓrger model and its derivatives have been widely used to incorporate transcytolemmal water exchange rate, an essential characteristic of living cells, into analyses of diffusion MRI (dMRI) signals from tissues. The Kӓrger model consists of two homogeneous exchanging components coupled by an exchange rate constant and assumes measurements are made with sufficiently long diffusion time and slow water exchange. Despite successful applications, it remains unclear whether these assumptions are generally valid for practical dMRI sequences and biological tissues. In particular, barrier-induced restrictions to diffusion produce inhomogeneous magnetization distributions in relatively large-sized compartments such as cancer cells, violating the above assumptions. The effects of this inhomogeneity are usually overlooked. We performed computer simulations to quantify how restriction effects, which in images produce edge enhancements at compartment boundaries, influence different variants of the Kӓrger-model. The results show that the edge enhancement effect will produce larger, time-dependent estimates of exchange rates in e.g., tumors with relatively large cell sizes (>10 μm), resulting in overestimations of water exchange as previously reported. Moreover, stronger diffusion gradients, longer diffusion gradient durations, and larger cell sizes, all cause more pronounced edge enhancement effects. This helps us to better understand the feasibility of the Kärger model in estimating water exchange in different tissue types and provides useful guidance on signal acquisition methods that may mitigate the edge enhancement effect. This work also indicates the need to correct the overestimated transcytolemmal water exchange rates obtained assuming the Kärger-model.

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

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

Exploring the potential of water channels for developing genetically encoded reporters and biosensors for diffusion-weighted MRI

Genetically encoded reporters for magnetic resonance imaging (MRI) offer a valuable technology for making molecular-scale measurements of biological processes within living organisms with high anatomical resolution and whole-organ coverage without relying on ionizing radiation. However, most MRI reporters rely on synthetic contrast agents, typically paramagnetic metals and metal complexes, which often need to be supplemented exogenously to create optimal contrast. To eliminate the need for synthetic contrast agents, we previously introduced aquaporin-1, a mammalian water channel, as a new reporter gene for the fully autonomous detection of genetically labeled cells using diffusion-weighted MRI. In this study, we aimed to expand the toolbox of diffusion-based genetic reporters by modulating aquaporin membrane trafficking and harnessing the evolutionary diversity of water channels across species. We identified a number of new water channels that functioned as diffusion-weighted reporter genes. In addition, we show that loss-of-function variants of yeast and human aquaporins can be leveraged to design first-in-class diffusion-based sensors for detecting the activity of a model protease within living cells.

Destabilized reporters for background-subtracted, chemically-gated, and multiplexed deep-tissue imaging

Tracking gene expression in deep tissues requires genetic reporters that can be unambiguously detected using tissue penetrant techniques. Magnetic resonance imaging (MRI) is uniquely suited for this purpose; however, there is a dearth of reporters that can be reliably linked to gene expression with minimal interference from background tissue signals. Here, we present a conceptually new method for generating background-subtracted, drug-gated, multiplex images of gene expression using MRI. Specifically, we engineered chemically erasable reporters consisting of a water channel, aquaporin-1, fused to destabilizing domains, which are stabilized by binding to cell-permeable small-molecule ligands. We showed that this approach allows for highly specific detection of gene expression through differential imaging. In addition, by engineering destabilized aquaporin-1 variants with orthogonal ligand requirements, it is possible to distinguish distinct subpopulations of cells in mixed cultures. Finally, we demonstrated this approach in a mouse tumor model through differential imaging of gene expression with minimal background.

A Magnetic Resonance Imaging-Chemical Exchange Saturation Transfer (MRI-CEST) Method for the Detection of Water Cycling across Cellular Membranes

Water cycling across the membrane transporters is considered a hallmark of cellular metabolism and it could be of high diagnostic relevance in the characterization of tumors and other diseases. The method relies on the response of intracellular proton exchanging molecules to the presence of extracellular Gd-based contrast agents (GBCAs). Paramagnetic GBCAs enhances the relaxation rate of water molecules in the extracellular compartment and, through membrane exchange, the relaxation enhancement is transferred to intracellular molecules. The effect is detected at the MRI-CEST (Magnetic Resonance Imaging - Chemical Exchange Saturation Transfer) signal of intracellular proton exchanging molecules. The magnitude of the change in the CEST response reports on water cycling across the membrane. The method has been tested on Red Blood Cells and on orthotopic murine models of breast cancer with different degree of malignancy (4T1, TS/A and 168FARN). The distribution of voxels reporting on membrane permeability fits well with the cells' aggressiveness and acts as an early reporter to monitor therapeutic treatments.

Detection of Treatment Response in Triple-Negative Breast Tumors to Paclitaxel Using MRI Cell Size Imaging

BACKGROUND

Breast cancer treatment response evaluation using the response evaluation criteria in solid tumors (RECIST) guidelines, based on tumor volume changes, has limitations, prompting interest in novel imaging markers for accurate therapeutic effect determination.

PURPOSE

To use MRI-measured cell size as a new imaging biomarker for assessing chemotherapy response in breast cancer.

STUDY TYPE

Longitudinal; animal model.

STUDY POPULATION

Triple-negative human breast cancer cell (MDA-MB-231) pellets (4 groups, n = 7) treated with dimethyl sulfoxide (DMSO) or 10 nM of paclitaxel for 24, 48, and 96 hours, and 29 mice with MDA-MB-231 tumors in right hind limbs treated with paclitaxel (n = 16) or DMSO (n = 13) twice weekly for 3 weeks.

FIELD STRENGTH/SEQUENCE

Oscillating gradient spin echo and pulsed gradient spin echo sequences at 4.7 T.

ASSESSMENT

MDA-MB-231 cells were analyzed using flowcytometry and light microscopy to assess cell cycle phases and cell size distribution. MDA-MB-231 cell pellets were MR imaged. Mice were imaged weekly, with 9, 6, and 14 being sacrificed for histology after MRI at weeks 1, 2, and 3, respectively. Microstructural parameters of tumors/cell pellets were derived by fitting diffusion MRI data to a biophysical model.

STATISTICAL TESTS

One-way ANOVA compared cell sizes and MR-derived parameters between treated and control samples. Repeated measures 2-way ANOVA with Bonferroni post-tests compared temporal changes in MR-derived parameters. A P-value <0.05 was considered statistically significant.

RESULTS

In vitro experiments showed that the mean MR-derived cell sizes of paclitaxel-treated cells increased significantly with a 24-hours treatment and decreased (P = 0.06) with a 96-hour treatment. For in vivo xenograft experiments, the paclitaxel-treated tumors showed significant decreases in cell size at later weeks. MRI observations were supported by flowcytometry, light microscopy, and histology.

DATA CONCLUSIONS

MR-derived cell size may characterize the cell shrinkage during treatment-induced apoptosis, and may potentially provide new insights into the assessment of therapeutic response.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 4.

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.

Engineering ligand stabilized aquaporin reporters for magnetic resonance imaging

Imaging transgene expression in live tissues requires reporters that are detectable with deeply penetrant modalities, such as magnetic resonance imaging (MRI). Here, we show that LSAqp1, a water channel engineered from aquaporin-1, can be used to create background-free, drug-gated, and multiplex images of gene expression using MRI. LSAqp1 is a fusion protein composed of aquaporin-1 and a degradation tag that is sensitive to a cell-permeable ligand, which allows for dynamic small molecule modulation of MRI signals. LSAqp1 improves specificity for imaging gene expression by allowing reporter signals to be conditionally activated and distinguished from the tissue background by difference imaging. In addition, by engineering destabilized aquaporin-1 variants with different ligand requirements, it is possible to image distinct cell types simultaneously. Finally, we expressed LSAqp1 in a tumor model and showed successful in vivo imaging of gene expression without background activity. LSAqp1 provides a conceptually unique approach to accurately measure gene expression in living organisms by combining the physics of water diffusion and biotechnology tools to control protein stability.

Optimized three-dimensional ultrashort echo time: Magnetic resonance fingerprinting for myelin tissue fraction mapping

Quantitative assessment of brain myelination has gained attention for both research and diagnosis of neurological diseases. However, conventional pulse sequences cannot directly acquire the myelin-proton signals due to its extremely short T2 and T2* values. To obtain the myelin-proton signals, dedicated short T2 acquisition techniques, such as ultrashort echo time (UTE) imaging, have been introduced. However, it remains challenging to isolate the myelin-proton signals from tissues with longer T2. In this article, we extended our previous two-dimensional ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) with dual-echo acquisition to three dimensional (3D). Given a relatively low proton density (PD) of myelin-proton, we utilized Cramér-Rao Lower Bound to encode myelin-proton with the maximal SNR efficiency for optimizing the MR fingerprinting design, in order to improve the sensitivity of the sequence to myelin-proton. In addition, with a second echo of approximately 3 ms, myelin-water component can be also captured. A myelin-tissue (myelin-proton and myelin-water) fraction mapping can be thus calculated. The optimized 3D UTE-MRF with dual-echo acquisition is tested in simulations, physical phantom and in vivo studies of both healthy subjects and multiple sclerosis patients. The results suggest that the rapidly decayed myelin-proton and myelin-water signal can be depicted with UTE signals of our method at clinically relevant resolution (1.8 mm isotropic) in 15 min. With its good sensitivity to myelin loss in multiple sclerosis patients demonstrated, our method for the whole brain myelin-tissue fraction mapping in clinical friendly scan time has the potential for routine clinical imaging.

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

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

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 &gt; 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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Diffusion time dependence, power-law scaling, and exchange in gray matter

Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins' interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b-1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.

Diffusion MRI signal cumulants and hepatocyte microstructure at fixed diffusion time: Insights from simulations, 9.4T imaging, and histology

Purpose

Relationships between diffusion‐weighted MRI signals and hepatocyte microstructure were investigated to inform liver diffusion MRI modeling, focusing on the following question: Can cell size and diffusivity be estimated at fixed diffusion time, realistic SNR, and negligible contribution from extracellular/extravascular water and exchange?

Methods

Monte Carlo simulations were performed within synthetic hepatocytes for varying cell size/diffusivity L/D0, and clinical protocols (single diffusion encoding; maximum b‐value: {1000, 1500, 2000} s/mm²; 5 unique gradient duration/separation pairs; SNR = {∞, 100, 80, 40, 20}), accounting for heterogeneity in (D0,L) and perfusion contamination. Diffusion (D) and kurtosis (K) coefficients were calculated, and relationships between (D0,L) and (D,K) were visualized. Functions mapping (D,K) to (D0,L) were computed to predict unseen (D0,L) values, tested for their ability to classify discrete cell‐size contrasts, and deployed on 9.4T ex vivo MRI‐histology data of fixed mouse livers

Results

Relationships between (D,K) and (D0,L) are complex and depend on the diffusion encoding. Functions mapping D,K to (D0,L) captures salient characteristics of D0(D,K) and L(D,K) dependencies. Mappings are not always accurate, but they enable just under 70% accuracy in a three‐class cell‐size classification task (for SNR = 20, bmax = 1500 s/mm², δ = 20 ms, and Δ = 75 ms). MRI detects cell‐size contrasts in the mouse livers that are confirmed by histology, but overestimates the largest cell sizes.

Conclusion

Salient information about liver cell size and diffusivity may be retrieved from minimal diffusion encodings at fixed diffusion time, in experimental conditions and pathological scenarios for which extracellular, extravascular water and exchange are negligible.

Measurement of cellular-interstitial water exchange time in tumors based on diffusion-time-dependent diffusional kurtosis imaging

PURPOSE

To assess the feasibility of using diffusion-time-dependent diffusional kurtosis imaging (tDKI) to measure cellular-interstitial water exchange time (τex ) in tumors, both in animals and in humans.

METHODS

Preclinical tDKI studies at 7 T were performed with the GL261 glioma model and the 4T1 mammary tumor model injected into the mouse brain. Clinical studies were performed at 3 T with women who had biopsy-proven invasive ductal carcinoma. tDKI measurement was conducted using a diffusion-weighted STEAM pulse sequence with multiple diffusion times (20-800 ms) at a fixed echo time, while keeping the b-values the same (0-3000 s/mm2 ) by adjusting the diffusion gradient strength. The tDKI data at each diffusion time t were used for a weighted linear least-squares fit method to estimate the diffusion-time-dependent diffusivity, D(t), and diffusional kurtosis, K(t).

RESULTS

Both preclinical and clinical studies showed that, when diffusion time t ≥ 200 ms, D(t) did not have a noticeable change while K(t) decreased monotonically with increasing diffusion time in tumors and t ≥ 100 ms for the cortical ribbon of the mouse brain. The estimated τex averaged median and interquartile range (IQR) of GL261 and 4T1 tumors were 93 (IQR = 89) ms and 68 (78) ms, respectively. For the cortical ribbon, the estimated τex averaged median and IQR were 41 (34) ms for C57BL/6 and 30 (17) ms for BALB/c. For invasive ductal carcinoma, the estimated τex median and IQR of the two breast cancers were 70 (94) and 106 (92) ms.

CONCLUSION

The results of this proof-of-concept study substantiate the feasibility of using tDKI to measure cellular-interstitial water exchange time without using an exogenous contrast agent.

Microstructure Modeling of High b-Value Diffusion-Weighted Images in Glioblastoma

Apparent diffusion coefficient has limits to differentiate solid tumor from normal tissue or edema in glioblastoma (GBM). This study investigated a microstructure model (MSM) in GBM using a clinically available diffusion imaging technique. The MSM was modified to integrate with bi-polar diffusion gradient waveforms, and applied to 30 patients with newly diagnosed GBM. Diffusion-weighted (DW) images acquired on a 3 T scanner with b-values from 0 to 2500 s/mm² were fitted in volumes of interest (VOIs) of solid tumor to obtain the apparent restriction size of intracellular water (ARS), the fractional volume of intracellular water (V_in), and extracellular (D_ex) water diffusivity. The parameters in solid tumor were compared with those of other tissue types by Students’ t test. For comparison, DW images were fitted by conventional mono-exponential and bi-exponential models. ARS, D_ex, and V_in from the MSM in tumor VOIs were significantly greater than those in WM, GM, and edema (P values of .01–.001). ARS values in solid tumors (from 21.6 to 34.5 um) had absolutely no overlap with those in all other tissue types (from 0.9 to 3.5 um). V_in values showed a descending order from solid tumor (from 0.32 to 0.52) to WM, GM, and edema (from 0.05 to 0.25), consisting with the descending cellularity in these tissue types. The parameters from mono-exponential and bi-exponential models could not significantly differentiate solid tumor from all other tissue types, particularly from edema. Further development and histopathological validation of the MSM will warrant its role in clinical management of GBM.

MR cell size imaging with temporal diffusion spectroscopy

Cytological features such as cell size and intracellular morphology provide fundamental information on cell status and hence may provide specific information on changes that arise within biological tissues. Such information is usually obtained by invasive biopsy in current clinical practice, which suffers several well-known disadvantages. Recently, novel MRI methods such as IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) have been developed for direct measurements of mean cell size non-invasively. The IMPULSED protocol is based on using temporal diffusion spectroscopy (TDS) to combine measurements of water diffusion over a wide range of diffusion times to probe cellular microstructure over varying length scales. IMPULSED has been shown to provide rapid, robust, and reliable mapping of mean cell size and is suitable for clinical imaging. More recently, cell size distributions have also been derived by appropriate analyses of data acquired with IMPULSED or similar sequences, which thus provides MRI-cytometry. This review summarizes the basic principles, practical implementations, validations, and example applications of MR cell size imaging based on TDS and demonstrates how cytometric information can be used in various applications. In addition, the limitations and potential future directions of MR cytometry are identified including the diagnosis of nonalcoholic steatohepatitis of the liver and the assessment of treatment response of cancers.

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.

A simulation study of cell size and volume fraction mapping for tissue with two underlying cell populations using diffusion-weighted MRI

PURPOSE

To propose a method for voxel-wise estimation of cell radii and volume fractions of two cell populations when they coexist in the same MR voxel using the combination of diffusion-weighted MRI and microstructural modeling.

METHOD

Microstructure models were investigated using diffusion data simulated with the matrix method for a range of microstructures mimicking tumor tissue with two cell populations, using acquisition parameters available on preclinical scanners. The effect of noise was investigated for a subset of these microstructures. The accuracy and precision of the estimated radii and volume fractions for large and small cells R l , R s , v i n , l , v i n , s were evaluated by comparing the estimates to their true values. The stability of model fitting was characterized by the percentage of accepted fits.

RESULTS

The estimation accuracy and precision, and thus the ability to robustly distinguish the two cell populations, depended on the microstructural properties and SNR. For a SNR of 50, a minimum difference of 3 μm between the radius of the large and small cell populations was required for differentiation. Proposed modifications to the two cell population microstructure model, including constrained fits, improved the stability of fits.

CONCLUSIONS

This proof-of-concept study proposed a diffusion MRI-based method for voxel-wise estimation of cell radii and volume fractions of two cell populations when they coexist in the same MR voxel. The ability to reliably characterize tissue with two cell populations opens exciting avenues of potential applications in both tumor diagnosis and treatment monitoring.

MRI-cytometry: Mapping nonparametric cell size distributions using diffusion MRI

PURPOSE

This report introduces and validates a new diffusion MRI-based method, termed MRI-cytometry, which can noninvasively map intravoxel, nonparametric cell size distributions in tissues.

METHODS

MRI was used to acquire diffusion MRI signals with a range of diffusion times and gradient factors, and a model was fit to these data to derive estimates of cell size distributions. We implemented a 2-step fitting method to avoid noise-induced artificial peaks and provide reliable estimates of tumor cell size distributions. Computer simulations in silico, experimental measurements on cultured cells in vitro, and animal xenografts in vivo were used to validate the accuracy and precision of the method. Tumors in 7 patients with breast cancer were also imaged and analyzed using this MRI-cytometry approach on a clinical 3 Tesla MRI scanner.

RESULTS

Simulations and experimental results confirm that MRI-cytometry can reliably map intravoxel, nonparametric cell size distributions and has the potential to discriminate smaller and larger cells. The application in breast cancer patients demonstrates the feasibility of direct translation of MRI-cytometry to clinical applications.

CONCLUSION

The proposed MRI-cytometry method can characterize nonparametric cell size distributions in human tumors, which potentially provides a practical imaging approach to derive specific histopathological information on biological tissues.

Challenges for biophysical modeling of microstructure

The biophysical modeling efforts in diffusion MRI have grown considerably over the past 25 years. In this review, we dwell on the various challenges along the journey of bringing a biophysical model from initial design to clinical implementation, identifying both hurdles that have been already overcome and outstanding issues. First, we describe the critical initial task of selecting which features of tissue microstructure can be estimated using a model and which acquisition protocol needs to be implemented to make the estimation possible. The model performance should necessarily be tested in realistic numerical simulations and in experimental data - adapting the fitting strategy accordingly, and parameter estimates should be validated against complementary techniques, when/if available. Secondly, the model performance and validity should be explored in pathological conditions, and, if appropriate, dedicated models for pathology should be developed. We build on examples from tumors, ischemia and demyelinating diseases. We then discuss the challenges associated with clinical translation and added value. Finally, we single out four major unresolved challenges that are related to: the availability of a microstructural ground truth, the validation of model parameters which cannot be accessed with complementary techniques, the development of a generalized standard model for any brain region and pathology, and the seamless communication between different parties involved in the development and application of biophysical models of diffusion.

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.

Diffusion model comparison identifies distinct tumor sub-regions and tracks treatment response

PURPOSE

MRI biomarkers of tumor response to treatment are typically obtained from parameters derived from a model applied to pre-treatment and post-treatment data. However, as tumors are spatially and temporally heterogeneous, different models may be necessary in different tumor regions, and model suitability may change over time. This work evaluates how the suitability of two diffusion-weighted (DW) MRI models varies spatially within tumors at the voxel level and in response to radiotherapy, potentially allowing inference of qualitatively different tumor microenvironments.

METHODS

DW-MRI data were acquired in CT26 subcutaneous allografts before and after radiotherapy. Restricted and time-independent diffusion models were compared, with regions well-described by the former hypothesized to reflect cellular tissue, and those well-described by the latter expected to reflect necrosis or oedema. Technical and biological validation of the percentage of tissue described by the restricted diffusion microstructural model (termed %MM) was performed through simulations and histological comparison.

RESULTS

Spatial and radiotherapy-related variation in model suitability was observed. %MM decreased from a mean of 64% at baseline to 44% 6 days post-radiotherapy in the treated group. %MM correlated negatively with the percentage of necrosis from histology, but overestimated it due to noise. Within MM regions, microstructural parameters were sensitive to radiotherapy-induced changes.

CONCLUSIONS

There is spatial and radiotherapy-related variation in different models' suitability for describing diffusion in tumor tissue, suggesting the presence of different and changing tumor sub-regions. The biological and technical validation of the proposed %MM cancer imaging biomarker suggests it correlates with, but overestimates, the percentage of necrosis.

Magnetic resonance imaging of mean cell size in human breast tumors

PURPOSE

Cell size is a fundamental characteristic of all tissues, and changes in cell size in cancer reflect tumor status and response to treatments, such as apoptosis and cell-cycle arrest. Unfortunately, cell size can currently be obtained only by pathological evaluation of tumor tissue samples obtained invasively. Previous imaging approaches are limited to preclinical MRI scanners or require relatively long acquisition times that are impractical for clinical imaging. There is a need to develop cell-size imaging for clinical applications.

METHODS

We propose a clinically feasible IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) approach that can characterize mean cell sizes in solid tumors. We report the use of a combination of pulse sequences, using different gradient waveforms implemented on clinical MRI scanners and analytical equations based on these waveforms to analyze diffusion-weighted MRI signals and derive specific microstructural parameters such as cell size. We also describe comprehensive validations of this approach using computer simulations, cell experiments in vitro, and animal experiments in vivo and demonstrate applications in preoperative breast cancer patients.

RESULTS

With fast acquisitions (~7 minutes), IMPULSED can provide high-resolution (1.3 mm in-plane) mapping of mean cell size of human tumors in vivo on clinical 3T MRI scanners. All validations suggest that IMPULSED provides accurate and reliable measurements of mean cell size.

CONCLUSION

The proposed IMPULSED method can assess cell-size variations in tumors of breast cancer patients, which may have the potential to assess early response to neoadjuvant therapy.

Mapping hepatocyte size in vivo using temporal diffusion spectroscopy MRI

PURPOSE

The goal of this study is to implement a noninvasive method for in vivo mapping of hepatocyte size. This method will have a broad range of clinical and preclinical applications, as pathological changes in hepatocyte sizes are relevant for the accurate diagnosis and assessments of treatment response of liver diseases.

METHODS

Building on the concepts of temporal diffusion spectroscopy in MRI, a clinically feasible imaging protocol named IMPULSED (Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion) has been developed, which is able to report measurements of cell sizes noninvasively. This protocol acquires a selected set of diffusion imaging data and fits them to a model of water compartments in tissues to derive robust estimates of the cellular structures that restrict free diffusion. Here, we adapt and further develop this approach to measure hepatocyte sizes in vivo. We validated IMPULSED in livers of mice and rats and implemented it to image healthy human subjects using a clinical 3T MRI scanner.

RESULTS

The IMPULSED-derived mean hepatocyte sizes for rats and mice are about 15-20 µm and agree well with histological findings. Maps of mean hepatocyte size for humans can be achieved in less than 15 minutes, a clinically feasible scan time.

CONCLUSION

Our results suggest that this method has potential to overcome major limitations of liver biopsy and provide noninvasive mapping of hepatocyte sizes in clinical applications.

SANDI: A compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI

This work introduces a compartment-based model for apparent cell body (namely soma) and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging through DW-MRI presents water diffusion in white (WM) and gray (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b ​≤ ​3,000 ​s/mm² (or 3 ​ms/μm²), it has been also shown to fail in GM at high b values (b≫3,000 ​s/mm² or 3 ​ms/μm²). Here we hypothesise that the unmodelled soma compartment (i.e. cell body of any brain cell type: from neuroglia to neurons) may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma of any brain cell type is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (b_max ​= ​40,000 ​s/mm², or 40 ​ms/μm²) and human (b_max ​= ​10,000 ​s/mm², or 10 ​ms/μm²) brain. We show that SANDI defines new contrasts representing complementary information on the brain cyto- and myelo-architecture. Indeed, we show maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto- and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.

Higher-order diffusion MRI characterization of mesorectal lymph nodes in rectal cancer

PURPOSE

Mesorectal lymph node staging plays an important role in treatment decision making. Here, we explore the benefit of higher-order diffusion MRI models accounting for non-Gaussian diffusion effects to classify mesorectal lymph nodes both 1) ex vivo at ultrahigh field correlated with histology and 2) in vivo in a clinical scanner upon patient staging.

METHODS

The preclinical investigation included 54 mesorectal lymph nodes, which were scanned at 16.4 T with an extensive diffusion MRI acquisition. Eight diffusion models were compared in terms of goodness of fit, lymph node classification ability, and histology correlation. In the clinical part of this study, 10 rectal cancer patients were scanned with diffusion MRI at 1.5 T, and 72 lymph nodes were analyzed with Apparent Diffusion Coefficient (ADC), Intravoxel Incoherent Motion (IVIM), Kurtosis, and IVIM-Kurtosis.

RESULTS

Compartment models including restricted and anisotropic diffusion improved the preclinical data fit, as well as the lymph node classification, compared to standard ADC. The comparison with histology revealed only moderate correlations, and the highest values were observed between diffusion anisotropy metrics and cell area fraction. In the clinical study, the diffusivity from IVIM-Kurtosis was the only metric showing significant differences between benign (0.80 ± 0.30 μm² /ms) and malignant (1.02 ± 0.41 μm² /ms, P = .03) nodes. IVIM-Kurtosis also yielded the largest area under the receiver operating characteristic curve (0.73) and significantly improved the node differentiation when added to the standard visual analysis by experts based on T₂ -weighted imaging.

CONCLUSION

Higher-order diffusion MRI models perform better than standard ADC and may be of added value for mesorectal lymph node classification in rectal cancer patients.

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.

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.

INDIANA: An in-cell diffusion method to characterize the size, abundance and permeability of cells

NMR and MRI diffusion experiments contain information describing the shape, size, abundance, and membrane permeability of cells although extracting this information can be challenging. Here we present the INDIANA (IN-cell DIffusion ANAlysis) method to simultaneously and non-invasively measure cell abundance, effective radius, permeability and intrinsic relaxation rates and diffusion coefficients within the inter- and intra-cellular populations. The method couples an experimental dataset comprising stimulated-echo diffusion measurements, varying both the gradient strength and the diffusion delay, together with software to fit a model based on the Kärger equations to robustly extract the relevant parameters. A detailed error analysis is presented by comparing the results from fitting simulated data from Monte Carlo simulations, establishing its effectiveness. We note that for parameters typical of mammalian cells the approach is particularly effective, and the shape of the underlying cells does not unduly affect the results. Finally, we demonstrate the performance of the experiment on systems of suspended yeast and mammalian cells. The extracted parameters describing cell abundance, size, permeability and relaxation are independently validated.

Quantitative temporal diffusion spectroscopy as an early imaging biomarker of radiation therapeutic response in gliomas: A preclinical proof of concept

PURPOSE

This study aims to test the ability of quantitative temporal diffusion spectroscopy (qTDS) to assess cellular changes associated with radiation-induced cell death in a rat glioma model.

METHODS AND MATERIALS

Tumor response to a single fraction of 20 Gy of x-ray radiation was investigated in a rat glioma (9L) model. Tumor response was monitored longitudinally at postinoculation days 21, 23, and 25, using a specific implementation of qTDS with acronym IMPULSED (Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion), as well as conventional diffusion and high-resolution anatomic imaging. IMPULSED method combines diffusion-weighted signals acquired over a range of diffusion times that are then analyzed and interpreted using a theoretical model of water diffusion in tissues, which generates parametric maps depicting cellular and subcellular structural information on a voxel-wise basis. Results from different metrics were compared statistically.

RESULTS

A single dose of 20 Gy x-ray radiation significantly prolonged survival of 9L-bearing rats. The mean cell sizes of irradiated tumors decreased (P < .005) after radiation treatment, which we associate with cell shrinkage and the formation of small cellular bodies during apoptosis and necrosis. A combination of IMPULSED-derived parameters (mean cell size d and extracellular structural parameter β ex ) separated 90% of irradiated tumors from the nonirradiated cases at post inoculation day 23, whereas a combination of tumor growth and conventional apparent diffusion coefficient did not differentiate irradiated tumors from nonirradiated tumors.

CONCLUSIONS

This proof-of-concept study demonstrates the IMPULSED method to be a new method for deriving quantitative microstructural parameters in a preclinical tumor model. The method provides unique information based on the diffusion time dependency of diffusion magnetic resonance imaging, which cannot be obtained by conventional diffusion weighted imaging methods, and the results have a close correlation with primary biologic markers of treatment efficacy, such as cell death and survival.

A biomimetic tumor tissue phantom for validating diffusion-weighted MRI measurements

PURPOSE

To develop a biomimetic tumor tissue phantom which more closely reflects water diffusion in biological tissue than previously used phantoms, and to evaluate the stability of the phantom and its potential as a tool for validating diffusion-weighted (DW) MRI measurements.

METHODS

Coaxial-electrospraying was used to generate micron-sized hollow polymer spheres, which mimic cells. The bulk structure was immersed in water, providing a DW-MRI phantom whose apparent diffusion coefficient (ADC) and microstructural properties were evaluated over a period of 10 months. Independent characterization of the phantom's microstructure was performed using scanning electron microscopy (SEM). The repeatability of the construction process was investigated by generating a second phantom, which underwent high resolution synchrotron-CT as well as SEM and MR scans.

RESULTS

ADC values were stable (coefficients of variation (CoVs) < 5%), and varied with diffusion time, with average values of 1.44 ± 0.03 µm2 /ms (Δ = 12 ms) and 1.20 ± 0.05 µm2 /ms (Δ = 45 ms). Microstructural parameters showed greater variability (CoVs up to 13%), with evidence of bias in sphere size estimates. Similar trends were observed in the second phantom.

CONCLUSION

A novel biomimetic phantom has been developed and shown to be stable over 10 months. It is envisaged that such phantoms will be used for further investigation of microstructural models relevant to characterizing tumor tissue, and may also find application in evaluating acquisition protocols and comparing DW-MRI-derived biomarkers obtained from different scanners at different sites. Magn Reson Med 80:147-158, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Microstructure Characterization of Bone Metastases from Prostate Cancer with Diffusion MRI: Preliminary Findings

PURPOSE

To examine the usefulness of rich diffusion protocols with high b-values and varying diffusion time for probing microstructure in bone metastases. Analysis techniques including biophysical and mathematical models were compared with the clinical apparent diffusion coefficient (ADC).

METHODS

Four patients were scanned using 13 b-values up to 3,000 s/mm2 and diffusion times ranging 18-52 ms. Data were fitted to mono-exponential ADC, intravoxel incoherent motion (IVIM), Kurtosis and Vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT) models. Parameters from the models were compared using correlation plots.

RESULTS

ADC and IVIM did not fit the data well, failing to capture the signal at high b-values. The Kurtosis model best explained the data in many voxels, but in voxels exhibiting a more time-dependent signal, the VERDICT model explained the data best. The ADC correlated significantly (p < 0.004) with the intracellular diffusion coefficient (r = 0.48), intracellular volume fraction (r = -0.21), and perfusion fraction (r = 0.46) parameters from VERDICT, suggesting that these factors all contribute to ADC contrast. The mean kurtosis correlated with the intracellular volume fraction parameter (r = 0.26) from VERDICT, consistent with the hypothesis that kurtosis relates to cellularity, but also correlated weakly with the intracellular diffusion coefficient (r = 0.18) and cell radius (r = 0.16) parameters, suggesting that it may be difficult to attribute physical meaning to kurtosis.

CONCLUSION

Both Kurtosis and VERDICT explained the diffusion signal better than ADC and IVIM, primarily due to poor fitting at high b-values in the latter two models. The Kurtosis and VERDICT models captured information at high b using parameters (Kurtosis or intracellular volume fraction and radius) that do not have a simple relationship with ADC and that may provide additional microstructural information in bone metastases.

Intracellular water preexchange lifetime in neurons and astrocytes

PURPOSE

To determine the intracellular water preexchange lifetime, τ_i , the "average residence time" of water, in the intracellular milieu of neurons and astrocytes. The preexchange lifetime is important for modeling a variety of MR data sets, including relaxation, diffusion-sensitive, and dynamic contrast-enhanced data sets.

METHODS

Herein, τ_i in neurons and astrocytes is determined in a microbead-adherent, cultured cell system. In concert with thin-slice selection, rapid flow of extracellular media suppresses extracellular signal, allowing determination of the transcytolemmal-exchange-dominated, intracellular T₁ . With this knowledge, and that of the intracellular T₁ in the absence of exchange, τ_i can be derived.

RESULTS

Under normal culture conditions, τ_i for neurons is 0.75 ± 0.05 s versus 0.57 ± 0.03 s for astrocytes. Both neuronal and astrocytic τ_i s decrease within 30 min after the onset of oxygen-glucose deprivation, with the astrocytic τ_i showing a substantially greater decrease than the neuronal τ_i .

CONCLUSIONS

Given an approximate intra- to extracellular volume ratio of 4:1 in the brain, these data imply that, under normal physiological conditions, an MR experimental characteristic time of less than 0.012 s is required for a nonexchanging, two-compartment (intra- and extracellular) model to be valid for MR studies. This characteristic time shortens significantly (i.e., 0.004 s) under injury conditions. Magn Reson Med 79:1616-1627, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

In vivo imaging of cancer cell size and cellularity using temporal diffusion spectroscopy

PURPOSE

A temporal diffusion MRI spectroscopy based approach has been developed to quantify cancer cell size and density in vivo.

METHODS

A novel imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED) method selects a specific limited diffusion spectral window for an accurate quantification of cell sizes ranging from 10 to 20 μm in common solid tumors. In practice, it is achieved by a combination of a single long diffusion time pulsed gradient spin echo (PGSE) and three low-frequency oscillating gradient spin echo (OGSE) acquisitions. To validate our approach, hematoxylin and eosin staining and immunostaining of cell membranes, in concert with whole slide imaging, were used to visualize nuclei and cell boundaries, and hence, enabled accurate estimates of cell size and cellularity.

RESULTS

Based on a two compartment model (incorporating intra- and extracellular spaces), accurate estimates of cell sizes were obtained in vivo for three types of human colon cancers. The IMPULSED-derived apparent cellularities showed a stronger correlation (r = 0.81; P < 0.0001) with histology-derived cellularities than conventional ADCs (r = -0.69; P < 0.03).

CONCLUSION

The IMPULSED approach samples a specific region of temporal diffusion spectra with enhanced sensitivity to length scales of 10-20 μm, and enables measurements of cell sizes and cellularities in solid tumors in vivo. Magn Reson Med 78:156-164, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Robust and fast nonlinear optimization of diffusion MRI microstructure models

Advances in biophysical multi-compartment modeling for diffusion MRI (dMRI) have gained popularity because of greater specificity than DTI in relating the dMRI signal to underlying cellular microstructure. A large range of these diffusion microstructure models have been developed and each of the popular models comes with its own, often different, optimization algorithm, noise model and initialization strategy to estimate its parameter maps. Since data fit, accuracy and precision is hard to verify, this creates additional challenges to comparability and generalization of results from diffusion microstructure models. In addition, non-linear optimization is computationally expensive leading to very long run times, which can be prohibitive in large group or population studies. In this technical note we investigate the performance of several optimization algorithms and initialization strategies over a few of the most popular diffusion microstructure models, including NODDI and CHARMED. We evaluate whether a single well performing optimization approach exists that could be applied to many models and would equate both run time and fit aspects. All models, algorithms and strategies were implemented on the Graphics Processing Unit (GPU) to remove run time constraints, with which we achieve whole brain dataset fits in seconds to minutes. We then evaluated fit, accuracy, precision and run time for different models of differing complexity against three common optimization algorithms and three parameter initialization strategies. Variability of the achieved quality of fit in actual data was evaluated on ten subjects of each of two population studies with a different acquisition protocol. We find that optimization algorithms and multi-step optimization approaches have a considerable influence on performance and stability over subjects and over acquisition protocols. The gradient-free Powell conjugate-direction algorithm was found to outperform other common algorithms in terms of run time, fit, accuracy and precision. Parameter initialization approaches were found to be relevant especially for more complex models, such as those involving several fiber orientations per voxel. For these, a fitting cascade initializing or fixing parameter values in a later optimization step from simpler models in an earlier optimization step further improved run time, fit, accuracy and precision compared to a single step fit. This establishes and makes available standards by which robust fit and accuracy can be achieved in shorter run times. This is especially relevant for the use of diffusion microstructure modeling in large group or population studies and in combining microstructure parameter maps with tractography results.

•

Evaluate robustness of fit, accuracy, precision for diffusion microstructure models

•

Test three optimization algorithms and three parameter initialization strategies

•

GPU implementation removes run time constraints; whole brain fit within minutes

•

Powell conjugate-direction algorithm has superior fit, accuracy, precision

•

Initialization approaches are important for crossing fiber microstructure models

Evaluation and comparison of diffusion MR methods for measuring apparent transcytolemmal water exchange rate constant

Two diffusion-based approaches, CG (constant gradient) and FEXI (filtered exchange imaging) methods, have been previously proposed for measuring transcytolemmal water exchange rate constant k_in, but their accuracy and feasibility have not been comprehensively evaluated and compared. In this work, both computer simulations and cell experiments in vitro were performed to evaluate these two methods. Simulations were done with different cell diameters (5, 10, 20μm), a broad range of k_in values (0.02-30s^-1) and different SNR's, and simulated k_in's were directly compared with the ground truth values. Human leukemia K562 cells were cultured and treated with saponin to selectively change cell transmembrane permeability. The agreement between measured k_in's of both methods was also evaluated. The results suggest that, without noise, the CG method provides reasonably accurate estimation of k_in especially when it is smaller than 10s^-1, which is in the typical physiological range of many biological tissues. However, although the FEXI method overestimates k_in even with corrections for the effects of extracellular water fraction, it provides reasonable estimates with practical SNR's and more importantly, the fitted apparent exchange rate AXR showed approximately linear dependence on the ground truth k_in. In conclusion, either CG or FEXI method provides a sensitive means to characterize the variations in transcytolemmal water exchange rate constant k_in, although the accuracy and specificity is usually compromised. The non-imaging CG method provides more accurate estimation of k_in, but limited to large volume-of-interest. Although the accuracy of FEXI is compromised with extracellular volume fraction, it is capable of spatially mapping k_in in practice.