Effects of intracellular organelles on the apparent diffusion coefficient of water molecules in cultured human embryonic kidney cells

Coregistered histology sections with diffusion tensor imaging data at 200 µm resolution in meningioma tumors

A significant problem in diffusion MRI (dMRI) is the lack of understanding regarding which microstructural features account for the variability in the diffusion tensor imaging (DTI) parameters observed in meningioma tumors. A common assumption is that mean diffusivity (MD) and fractional anisotropy (FA) from DTI are inversely proportional to cell density and proportional to tissue anisotropy, respectively. Although these associations have been established across a wide range of tumors, they have been challenged for interpreting within-tumor variations where several additional microstructural features have been suggested as contributing to MD and FA. To facilitate the investigation of the biological underpinnings of DTI parameters, we performed ex-vivo DTI at 200 µm isotropic resolution on sixteen excised meningioma tumor samples. The samples exhibit a variety of microstructural features because the dataset includes meningiomas of six different meningioma types and two different grades. Diffusion-weighted signal (DWI) maps, DWI maps averaged over all directions for given b-value, signal intensities without diffusion encoding (S₀) as well as DTI parameters: MD, FA, in-plane FA (FA_IP), axial diffusivity (AD) and radial diffusivity (RD), were coregistered to Hematoxylin & Eosin- (H&E) and Elastica van Gieson-stained (EVG) histological sections by a non-linear landmark-based approach. Here, we provide DWI signal and DTI maps coregistered to histology sections and describe the pipeline for processing the raw DTI data and the coregistration. The raw, processed, and coregistered data are hosted by Analytic Imaging Diagnostics Arena (AIDA) data hub registry, and software tools for processing are provided via GitHub. We hope that data can be used in research and education concerning the link between the meningioma microstructure and parameters obtained by DTI.

Meningioma microstructure assessed by diffusion MRI: An investigation of the source of mean diffusivity and fractional anisotropy by quantitative histology

BACKGROUND

Mean diffusivity (MD) and fractional anisotropy (FA) from diffusion MRI (dMRI) have been associated with cell density and tissue anisotropy across tumors, but it is unknown whether these associations persist at the microscopic level.

PURPOSE

To quantify the degree to which cell density and anisotropy, as determined from histology, account for the intra-tumor variability of MD and FA in meningioma tumors. Furthermore, to clarify whether other histological features account for additional intra-tumor variability of dMRI parameters.

MATERIALS AND METHODS

We performed ex-vivo dMRI at 200 μm isotropic resolution and histological imaging of 16 excised meningioma tumor samples. Diffusion tensor imaging (DTI) was used to map MD and FA, as well as the in-plane FA (FA_IP). Histology images were analyzed in terms of cell nuclei density (CD) and structure anisotropy (SA; obtained from structure tensor analysis) and were used separately in a regression analysis to predict MD and FA_IP, respectively. A convolutional neural network (CNN) was also trained to predict the dMRI parameters from histology patches. The association between MRI and histology was analyzed in terms of out-of-sample (R²_OS) on the intra-tumor level and within-sample R² across tumors. Regions where the dMRI parameters were poorly predicted from histology were analyzed to identify features apart from CD and SA that could influence MD and FA_IP, respectively.

RESULTS

Cell density assessed by histology poorly explained intra-tumor variability of MD at the mesoscopic level (200 μm), as median R²_OS = 0.04 (interquartile range 0.01-0.26). Structure anisotropy explained more of the variation in FA_IP (median R²_OS = 0.31, 0.20-0.42). Samples with low R²_OS for FA_IP exhibited low variations throughout the samples and thus low explainable variability, however, this was not the case for MD. Across tumors, CD and SA were clearly associated with MD (R² = 0.60) and FA_IP (R² = 0.81), respectively. In 37% of the samples (6 out of 16), cell density did not explain intra-tumor variability of MD when compared to the degree explained by the CNN. Tumor vascularization, psammoma bodies, microcysts, and tissue cohesivity were associated with bias in MD prediction based solely on CD. Our results support that FA_IP is high in the presence of elongated and aligned cell structures, but low otherwise.

CONCLUSION

Cell density and structure anisotropy account for variability in MD and FA_IP across tumors but cell density does not explain MD variations within the tumor, which means that low or high values of MD locally may not always reflect high or low tumor cell density. Features beyond cell density need to be considered when interpreting MD.

Utility of magnetic resonance imaging for differentiating malignant mesenchymal tumors of the uterus from T2-weighted hyperintense leiomyomas

PURPOSE

To generate a new discrimination method to distinguish between malignant mesenchymal tumors of the uterus and T2-weighted hyperintense leiomyoma based on magnetic resonance imaging findings and clinical features.

MATERIALS AND METHODS

Data from 32 tumors of 32 patients with malignant mesenchymal tumors of the uterus and from 34 tumors of 30 patients with T2-weighted hyperintense leiomyoma were analyzed. Clinical parameters, qualitative magnetic resonance imaging features, including computed diffusion-weighted imaging, and quantitative characteristics of magnetic resonance imaging of these two tumor types were compared. Predictive values for malignant mesenchymal tumors of the uterus were calculated using variant discriminant analysis.

RESULTS

The T1 bright area on qualitative assessment and mean apparent diffusion coefficient value on quantitative assessment yielded the most independent magnetic resonance imaging differentiators of malignant mesenchymal tumors of the uterus and T2-weighted hyperintense leiomyoma. The classification accuracy of the variant discriminant analysis based on three selected findings, i.e., a T1 bright area, computed diffusion-weighted imaging with a b-value of 2000s/mm² (cDWI₂₀₀₀), and T2-hypointense bands, was 84.8% (56/66), indicating high accuracy.

CONCLUSIONS

Variant discriminant analysis using the T1 bright area, cDWI₂₀₀₀, and T2-hypointense bands yielded high accuracy for differentiating between malignant mesenchymal tumors of the uterus and T2-weighted hyperintense leiomyoma.

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.

Probing neural tissues at small scales: Recent progress of oscillating gradient spin echo (OGSE) neuroimaging in humans

The detection sensitivity of diffusion MRI (dMRI) is dependent on diffusion times. A shorter diffusion time can increase the sensitivity to smaller length scales. However, the conventional dMRI uses the pulse gradient spin echo (PGSE) sequence that probes relatively long diffusion times only. To overcome this, the oscillating gradient spin echo (OGSE) sequence has been developed to probe much shorter diffusion times with hardware limitations on preclinical and clinical MRI systems. The OGSE sequence has been previously used on preclinical animal MRI systems. Recently, several studies have translated the OGSE sequence to humans on clinical MRI systems and achieved new information that is invisible using conventional PGSE sequence. This paper provides an overview of the recent progress of the OGSE neuroimaging in humans, including the technical improvements in the translation of the OGSE sequence to human imaging and various applications in different neurological disorders and stroke. Some possible future directions of the OGSE sequence are also discussed.

The apparent diffusion coefficient is a useful biomarker in predicting treatment response in patients with locally advanced rectal cancer

Background

Apparent diffusion coefficient (ADC) values achieve promising results in treatment response prediction in patients with several types of cancers.

Purpose

To determine whether ADC values predict neoadjuvant chemoradiation treatment (nCRT) response in patients with locally advanced rectal cancer (LARC).

Material and Methods

Forty-four patients with LARC who underwent magnetic resonance imaging scans before and after nCRT followed by delayed surgery were enrolled retrospectively. The sample was distributed as follows: responders (R), n = 8; and non-responders (Non-R), n = 36. Three markers of treatment response were considered: post-nCRT measures; ΔADC; and Δ%ADC. Statistical analysis included a Wilcoxon test, a Mann–Whitney U test, and a receiver operating characteristic (ROC) analysis in order to evaluate the diagnostic accuracy for each ADC value marker to differentiate between R and Non-R.

Results

Both minimum and mean ADC values were significantly higher after nCRT in the R group, while non-significant differences between basal and control ADC values were found in the non-R group. In addition, ΔADC and Δ%ADC exhibited increased values after nCRT in R when compared with non-R. ROC analysis revealed the following diagnostic performance parameters: post-nCRT: ADC_min = 1.05 × 10⁻³ mm²/s (sensitivity 61.1% and specificity 66.7%), ADC_mean = 1.50 × 10⁻³ mm²/s (sensitivity 72.2% and specificity 83.3%), ΔADC: ADC_min = 0.35 (sensitivity 66.7% and specificity 83.3%), ADC_mean = 0.50 (sensitivity 72% and specificity 83%); and Δ%ADC: ADC_min = 44% (sensitivity 66.7% and specificity 83.3%) and ADC_mean = 60% (sensitivity 83% and specificity 99%).

Conclusion

Our findings suggest that post-treatment rectal tumor ADC values, as well changes between pre- and post-treatment values, may be biomarkers for predicting treatment response in patients with LARC who underwent nCRT.

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.

Crowd Control: Effects of Physical Crowding on Cargo Movement in Healthy and Diseased Neurons

High concentration of cytoskeletal filaments, organelles, and proteins along with the space constraints due to the axon's narrow geometry lead inevitably to intracellular physical crowding along the axon of a neuron. Local cargo movement is essential for maintaining steady cargo transport in the axon, and this may be impeded by physical crowding. Molecular motors that mediate active transport share movement mechanisms that allow them to bypass physical crowding present on microtubule tracks. Many neurodegenerative diseases, irrespective of how they are initiated, show increased physical crowding owing to the greater number of stalled organelles and structural changes associated with the cytoskeleton. Increased physical crowding may be a significant factor in slowing cargo transport to synapses, contributing to disease progression and culminating in the dying back of the neuronal process. This review explores the idea that physical crowding can impede cargo movement along the neuronal process. We examine the sources of physical crowding and strategies used by molecular motors that might enable cargo to circumvent physically crowded locations. Finally, we describe sub-cellular changes in neurodegenerative diseases that may alter physical crowding and discuss the implications of such changes on cargo movement.

Oscillating-gradient diffusion magnetic resonance imaging detects acute subcellular structural changes in the mouse forebrain after neonatal hypoxia-ischemia

The recently developed oscillating-gradient diffusion MRI (OG-dMRI) technique extends our ability to examine brain structures at different spatial scales. In this study, we investigated the sensitivity of OG-dMRI in detecting cellular and subcellular structural changes in a mouse model of neonatal hypoxia ischemia (HI). Neonatal mice received unilateral HI injury or sham injury at postnatal day 10, followed by in vivo T₂-weighted and diffusion MRI of the brains at 3-6 h and 24 h after HI. Apparent diffusion coefficient (ADC) maps were acquired using conventional pulsed-gradient dMRI (PG-dMRI) and OG-dMRI with oscillating frequencies from 50 to 200 Hz. Pathology at cellular and subcellular levels was evaluated using neuronal, glial, and mitochondrial markers. We found significantly higher rates of ADC increase with oscillating frequencies (Δ_fADC) in the ipsilateral edema region, compared to the contralateral side, starting as early as 3 h after HI. Even in injured regions that showed no apparent change in PG-ADC or pseudo-normalized PG-ADC measurements, Δ_fADC remained significantly elevated. Histopathology showed swelling of sub-cellular structures in these regions with no apparent whole-cell level change. These results suggest that OG-dMRI is sensitive to subcellular structural changes in the brain after HI and is less susceptible to pseudo-normalization than PG-dMRI.

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.

Magnetic resonance temporal diffusion tensor spectroscopy of disordered anisotropic tissue

Molecular diffusion measured with diffusion weighted MRI (DWI) offers a probe for tissue microstructure. However, inferring microstructural properties from conventional DWI data is a complex inverse problem and has to account for heterogeneity in sizes, shapes and orientations of the tissue compartments contained within an imaging voxel. Alternative experimental means for disentangling the signal signatures of such features could provide a stronger link between the data and its interpretation. Double diffusion encoding (DDE) offers the possibility to factor out variation in compartment shapes from orientational dispersion of anisotropic domains by measuring the correlation between diffusivity in multiple directions. Time dependence of the diffusion is another effect reflecting the dimensions and distributions of barriers. In this paper we extend on DDE with a modified version of the oscillating gradient spin echo (OGSE) experiment, giving a basic contrast mechanism closely linked to both the temporal diffusion spectrum and the compartment anisotropy. We demonstrate our new method on post mortem brain tissue and show that we retrieve the correct temporal diffusion tensor spectrum in synthetic data from Monte Carlo simulations of random walks in a range of disordered geometries of different sizes and shapes.

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.

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

PURPOSE

To investigate the influence of transcytolemmal water exchange on estimates of tissue microstructural parameters derived from diffusion MRI using conventional PGSE and IMPULSED methods.

METHODS

Computer simulations were performed to incorporate a broad range of intracellular water life times τin (50-∞ ms), cell diameters d (5-15 μm), and intrinsic diffusion coefficient Din (0.6-2 μm2 /ms) for different values of signal-to-noise ratio (SNR) (10 to 50). For experiments, murine erythroleukemia (MEL) cancer cells were cultured and treated with saponin to selectively change cell membrane permeability. All fitted microstructural parameters from simulations and experiments in vitro were compared with ground-truth values.

RESULTS

Simulations showed that, for both PGSE and IMPULSED methods, cell diameter d can be reliably fit with sufficient SNR (≥ 50), whereas intracellular volume fraction fin is intrinsically underestimated due to transcytolemmal water exchange. Din can be reliably fit only with sufficient SNR and using the IMPULSED method with short diffusion times. These results were confirmed with those obtained in the cell culture experiments in vitro.

CONCLUSION

For the sequences and models considered in this study, transcytolemmal water exchange has minor effects on the fittings of d and Din with physiologically relevant membrane permeabilities if the SNR is sufficient (> 50), but fin is intrinsically underestimated. Magn Reson Med 77:2239-2249, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

A comparative assessment of preclinical chemotherapeutic response of tumors using quantitative non-Gaussian diffusion MRI

BACKGROUND

Diffusion-weighted MRI (DWI) signal attenuation is often not mono-exponential (i.e. non-Gaussian diffusion) with stronger diffusion weighting. Several non-Gaussian diffusion models have been developed and may provide new information or higher sensitivity compared with the conventional apparent diffusion coefficient (ADC) method. However the relative merits of these models to detect tumor therapeutic response is not fully clear.

METHODS

Conventional ADC, and three widely-used non-Gaussian models, (bi-exponential, stretched exponential, and statistical model), were implemented and compared for assessing SW620 human colon cancer xenografts responding to barasertib, an agent known to induce apoptosis via polyploidy. Bayesian Information Criterion (BIC) was used for model selection among all three non-Gaussian models.

RESULTS

All of tumor volume, histology, conventional ADC, and three non-Gaussian DWI models could show significant differences between control and treatment groups after four days of treatment. However, only the non-Gaussian models detected significant changes after two days of treatment. For any treatment or control group, over 65.7% of tumor voxels indicate the bi-exponential model is strongly or very strongly preferred.

CONCLUSION

Non-Gaussian DWI model-derived biomarkers are capable of detecting tumor earlier chemotherapeutic response of tumors compared with conventional ADC and tumor volume. The bi-exponential model provides better fitting compared with statistical and stretched exponential models for the tumor and treatment models used in the current work.

Fast and simplified mapping of mean axon diameter using temporal diffusion spectroscopy

Mapping axon diameter is of interest for the potential diagnosis and monitoring of various neuronal pathologies. Advanced diffusion-weighted MRI methods have been developed to measure mean axon diameters non-invasively, but suffer major drawbacks that prevent their direct translation into clinical practice, such as complex non-linear data fitting and, more importantly, long scanning times that are usually not tolerable for most human subjects. In the current study, temporal diffusion spectroscopy using oscillating diffusion gradients was used to measure mean axon diameters with high sensitivity to small axons in the central nervous system. Axon diameters have been found to be correlated with a novel metric, DDR⊥ (the rate of dispersion of the perpendicular diffusion coefficient with gradient frequency), which is a model-free quantity that does not require complex data analyses and can be obtained from two diffusion coefficient measurements in clinically relevant times with conventional MRI machines. A comprehensive investigation including computer simulations and animal experiments ex vivo showed that measurements of DDR⊥ agree closely with histological data. In humans in vivo, DDR⊥ was also found to correlate well with reported mean axon diameters in human corpus callosum, and the total scan time was only about 8 min. In conclusion, DDR⊥ may have potential to serve as a fast, simple and model-free approach to map the mean axon diameter of white matter in clinics for assessing axon diameter changes.

Influence of water compartmentation and heterogeneous relaxation on quantitative magnetization transfer imaging in rodent brain tumors

PURPOSE

The goal of this study was to investigate the influence of water compartmentation and heterogeneous relaxation properties on quantitative magnetization transfer (qMT) imaging in tissues, and in particular whether a two-pool model is sufficient to describe qMT data in brain tumors.

METHODS

Computer simulations and in vivo experiments with a series of qMT measurements before and after injection of Gd-DTPA were performed. Both off-resonance pulsed saturation (pulsed) and on-resonance selective inversion recovery (SIR) qMT methods were used, and all data were fit with a two-pool model only.

RESULTS

Simulations indicated that a two-pool fitting of four-pool data yielded accurate measures of pool size ratio (PSR) of macromolecular versus free water protons when there were fast transcytolemmal exchange and slow R1 recovery. The fitted in vivo PSR of both pulsed and SIR qMT methods showed no dependence on R1 variations caused by different concentrations of Gd-DTPA during wash-out, whereas the fitted kex (magnetization transfer exchange rate) changed significantly with R1 .

CONCLUSION

A two-pool model provides reproducible estimates of PSR in brain tumors independent of relaxation properties in the presence of relatively fast transcytolemmal exchange, whereas estimates of kex are biased by relaxation variations. In addition, estimates of PSR in brain tumors using the pulsed and SIR qMT methods agree well with one another. Magn Reson Med 76:635-644, 2016. © 2015 Wiley Periodicals, Inc.

Structural information revealed by the dispersion of ADC with frequency

Diffusion MRI provides a non-invasive means to characterize tissue microstructure at varying length scales. Temporal diffusion spectra reveal how the apparent diffusion coefficient (ADC) varies with frequency. When measured using oscillating gradient spin echo sequences, the manner in which ADC disperses with gradient frequency (which is related to the reciprocal of diffusion time) provides information on the characteristic dimensions of restricting structures within the medium. For example, the dispersion of ADC with oscillating gradient frequency (ΔfADC) has been shown to correlate with axon sizes in white matter and provide novel tissue contrast in images of mouse hippocampus and cerebellum. However, despite increasing interest in applying frequency-dependent ADC to derive novel information on tissue, the interpretations of ADC spectra are not always clear. In this study, the relation between ADC spectra and restricting dimensions are further elucidated and used to derive novel image contrast related to the sizes of intrinsic microstructures.

Time-Dependent Influence of Cell Membrane Permeability on MR Diffusion Measurements

PURPOSE

To investigate the influence of cell membrane permeability on diffusion measurements over a broad range of diffusion times.

METHODS

Human myelogenous leukemia K562 cells were cultured and treated with saponin to selectively alter cell membrane permeability, resulting in a broad physiologically relevant range of 0.011-0.044 μm/ms. Apparent diffusion coefficient (ADC) values were acquired with the effective diffusion time (Δeff ) ranging from 0.42 to 3000 ms. Cosine-modulated oscillating gradient spin echo (OGSE) measurements were performed to achieve short Δeff from 0.42 to 5 ms, while stimulated echo acquisitions were used to achieve long Δeff from 11 to 2999 ms. Computer simulations were also performed to support the experimental results.

RESULTS

Both computer simulations and experiments in vitro showed that the influence of membrane permeability on diffusion MR measurements is highly dependent on the choice of diffusion time, and it is negligible only when the diffusion time is at least one order of magnitude smaller than the intracellular exchange lifetime.

CONCLUSION

The influence of cell membrane permeability on the measured ADCs is negligible in OGSE measurements at moderately high frequencies. By contrast, cell membrane permeability has a significant influence on ADC and quantitative diffusion measurements at low frequencies such as those sampled using conventional pulsed gradient methods.

Quantification of cell size using temporal diffusion spectroscopy

PURPOSE

A new approach has been developed to quantify cell sizes and intracellular volume fractions using temporal diffusion spectroscopy with diffusion-weighted acquisitions.

METHODS

Temporal diffusion spectra may be used to characterize tissue microstructure by measuring the effects of restrictions over a range of diffusion times. Oscillating gradients have been used previously to probe variations on cellular and subcellular scales, but their ability to accurately measure cell sizes larger than 10 μm is limited. By combining measurements made using oscillating gradient spin echo (OGSE) and a conventional pulsed gradient spin echo (PGSE) acquisition with a single, relatively long diffusion time, we can accurately quantify cell sizes and intracellular volume fractions.

RESULTS

Based on a two compartment model (incorporating intra- and extracellular spaces), accurate estimates of cell sizes and intracellular volume fractions were obtained in vitro for (i) different cell types with sizes ranging from 10 to 20 μm, (ii) different cell densities, and (iii) before and after anticancer treatment.

CONCLUSION

Hybrid OGSE-PGSE acquisitions sample a larger region of temporal diffusion spectra and can accurately quantify cell sizes over a wide range. Moreover, the maximum gradient strength used was lower than 15 G/cm, suggesting that this approach is translatable to practical MR imaging.

Mapping mean axon diameter and axonal volume fraction by MRI using temporal diffusion spectroscopy

Mapping mean axon diameter and intra-axonal volume fraction may have significant clinical potential because nerve conduction velocity is directly dependent on axon diameter, and several neurodegenerative diseases affect axons of specific sizes and alter axon counts. Diffusion-weighted MRI methods based on the pulsed gradient spin echo (PGSE) sequence have been reported to be able to assess axon diameter and volume fraction non-invasively. However, due to the relatively long diffusion times used, e.g. >20ms, the sensitivity to small axons (diameter<2μm) is low, and the derived mean axon diameter has been reported to be overestimated. In the current study, oscillating gradient spin echo (OGSE) diffusion sequences with variable frequency gradients were used to assess rat spinal white matter tracts with relatively short effective diffusion times (1-5ms). In contrast to previous PGSE-based methods, the extra-axonal diffusion cannot be modeled as hindered (Gaussian) diffusion when short diffusion times are used. Appropriate frequency-dependent rates are therefore incorporated into our analysis and validated by histology-based computer simulation of water diffusion. OGSE data were analyzed to derive mean axon diameters and intra-axonal volume fractions of rat spinal white matter tracts (mean axon diameter of ~1.27-5.54μm). The estimated values were in good agreement with histology, including the small axon diameters (<2.5μm). This study establishes a framework for the quantification of nerve morphology using the OGSE method with high sensitivity to small axons.

Fast and robust measurement of microstructural dimensions using temporal diffusion spectroscopy

Mapping axon sizes non-invasively is of interest for neuroscientists and may have significant clinical potential because nerve conduction velocity is directly dependent on axon size. Current approaches to measuring axon sizes using diffusion-weighted MRI, e.g. q-space imaging with pulsed gradient spin echo (PGSE) sequences usually require long scan times and high q-values to detect small axons (diameter <2μm). The oscillating gradient spin echo (OGSE) method has been shown to be able to achieve very short diffusion times and hence may be able to detect smaller axons with high sensitivity. In the current study, OGSE experiments were performed to measure the inner diameters of hollow microcapillaries with a range of sizes (∼1.5-19.3μm) that mimic axons in the human central nervous system. The results suggest that OGSE measurements, even with only moderately high frequencies, are highly sensitive to compartment sizes, and a minimum of two ADC values with different frequencies may be sufficient to extract the microcapillary size accurately. This suggests that the OGSE method may serve as a fast and robust measurement method for mapping axon sizes non-invasively.

Diverse responses to vascular disrupting agent combretastatin a4 phosphate: a comparative study in rats with hepatic and subcutaneous tumor allografts using MRI biomarkers, microangiography, and histopathology

OBJECTIVE

Differently located tumors of the same origin may exhibit diverse responses to the same therapeutics. To test this hypothesis, we compared the responses of rodent hepatic and subcutaneous engrafts of rhabdomyosarcoma-1 (R1) to a vascular disrupting agent Combretastatin A4 phosphate (CA4P).

METHODS

Twelve WAG/Rij rats, each bearing three R1 implanted in the right and left hepatic lobes and subcutaneously in the thoracic region, received CA4P intravenously at 5 mg/kg (n = 6) or solvent (n = 6). Therapeutic responses were compared interindividually and intraindividually among tumors of different sites till 48 hours after injection using in vivo MRI, postmortem digital microangiography, and histopathology.

RESULTS

MRI revealed that the subcutaneous tumors (STs) significantly increased in volume than hepatic tumors (HTs) 48 hours after CA4P (P < .05). Relative to vehicle controls and treated group at baseline, necrosis ratio, apparent diffusion coefficient, and enhancement ratio changed slightly with the STs but significantly with HTs (P < .05) after CA4P treatment. Vessel density derived from microangiography was significantly lower in STs compared to HTs without CA4P treatment. CA4P treatment resulted in decreased vessel density in HTs, while it did not affect vessel density in STs. MRI and microangiography outcomes were supported by histopathologic findings.

CONCLUSIONS

MRI and microangiography allowed quantitative comparison of therapeutic responses to CA4P in rats with multifocal tumors. The discovered diverse effects of the same drug on tumors of the same origin but different locations emphasize the presence of cancer heterogeneity and the importance of individualization of drug delivery.

Characterizing tumor response to chemotherapy at various length scales using temporal diffusion spectroscopy

Measurements of apparent diffusion coefficient (ADC) using magnetic resonance imaging (MRI) have been suggested as potential imaging biomarkers for monitoring tumor response to treatment. However, conventional pulsed-gradient spin echo (PGSE) methods incorporate relatively long diffusion times, and are usually sensitive to changes in cell density and necrosis. Diffusion temporal spectroscopy using the oscillating gradient spin echo (OGSE) sequence is capable of probing short length scales, and may detect significant intracellular microstructural changes independent of gross cell density changes following anti-cancer treatment. To test this hypothesis, SW620 xenografts were treated by barasertib (AZD1152), a selective inhibitor of Aurora B kinase which causes SW620 cancer cells to develop polyploidy and increase in size following treatment, ultimately leading to cell death through apoptosis. Following treatment, the ADC values obtained by both the PGSE and low frequency OGSE methods increased. However, the ADC values at high gradient frequency (i.e. short diffusion times) were significantly lower in treated tumors, consistent with increased intracellular restrictions/hindrances. This suggests that ADC values at long diffusion times are dominated by tumor microstructure at long length scales, and may not convey unambiguous information of subcellular space. While the diffusion temporal spectroscopy provides more comprehensive means to probe tumor microstructure at various length scales. This work is the first study to probe intracellular microstructural variations due to polyploidy following treatment using diffusion MRI in vivo. It is also the first observation of post-treatment ADC changes occurring in opposite directions at short and long diffusion times. The current study suggests that temporal diffusion spectroscopy potentially provides pharmacodynamic biomarkers of tumor early response which distinguish microstructural variations following treatment at both the subcellular and supracellular length scales.