Influence of cell cycle phase on apparent diffusion coefficient in synchronized cells detected using temporal diffusion spectroscopy

Multiparametric MRI for characterization of the tumour microenvironment

Our understanding of tumour biology has evolved over the past decades and cancer is now viewed as a complex ecosystem with interactions between various cellular and non-cellular components within the tumour microenvironment (TME) at multiple scales. However, morphological imaging remains the mainstay of tumour staging and assessment of response to therapy, and the characterization of the TME with non-invasive imaging has not yet entered routine clinical practice. By combining multiple MRI sequences, each providing different but complementary information about the TME, multiparametric MRI (mpMRI) enables non-invasive assessment of molecular and cellular features within the TME, including their spatial and temporal heterogeneity. With an increasing number of advanced MRI techniques bridging the gap between preclinical and clinical applications, mpMRI could ultimately guide the selection of treatment approaches, precisely tailored to each individual patient, tumour and therapeutic modality. In this Review, we describe the evolving role of mpMRI in the non-invasive characterization of the TME, outline its applications for cancer detection, staging and assessment of response to therapy, and discuss considerations and challenges for its use in future medical applications, including personalized integrated diagnostics.

Diffusion time dependency of extracellular diffusion

PURPOSE

To quantify the variations of the power-law dependences on diffusion time t or gradient frequency f $$ f $$ of extracellular water diffusion measured by diffusion MRI (dMRI).

METHODS

Model cellular systems containing only extracellular water were used to investigate the t / f $$ t/f $$ dependence of D ex $$ {D}_{ex} $$ , the extracellular diffusion coefficient. Computer simulations used a randomly packed tissue model with realistic intracellular volume fractions and cell sizes. DMRI measurements were performed on samples consisting of liposomes containing heavy water(D₂ O, deuterium oxide) dispersed in regular water (H₂ O). D ex $$ {D}_{ex} $$ was obtained over a broad t $$ t $$ range (∼1-1000 ms) and then fit power-law equations D ex ( t ) = D const + const · t - ϑ t $$ {D}_{ex}(t)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {t}^{-{\vartheta}_t} $$ and D ex ( f ) = D const + const · f ϑ f $$ {D}_{ex}(f)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {f}^{\vartheta_f} $$ .

RESULTS

Both simulated and experimental results suggest that no single power-law adequately describes the behavior of D ex $$ {D}_{ex} $$ over the range of diffusion times of most interest in practical dMRI. Previous theoretical predictions are accurate over only limited t $$ t $$ ranges; for example, θ t = θ f = - 1 2 $$ {\theta}_t={\theta}_f=-\frac{1}{2} $$ is valid only for short times, whereas θ t = 1 $$ {\theta}_t=1 $$ or θ f = 3 2 $$ {\theta}_f=\frac{3}{2} $$ is valid only for long times but cannot describe other ranges simultaneously. For the specific t $$ t $$ range of 5-70 ms used in typical human dMRI measurements, θ t = θ f = 1 $$ {\theta}_t={\theta}_f=1 $$ matches the data well empirically.

CONCLUSION

The optimal power-law fit of extracellular diffusion varies with diffusion time. The dependency obtained at short or long t $$ t $$ limits cannot be applied to typical dMRI measurements in human cancer or liver. It is essential to determine the appropriate diffusion time range when modeling extracellular diffusion in dMRI-based quantitative microstructural imaging.

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.

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.

Comparison of pulsed and oscillating gradient diffusion-weighted MRI for characterizing hepatocellular nodules in liver cirrhosis: ex vivo study in a rat model

BACKGROUND

In contrast to classical pulsed gradient diffusion-weighted MRI, oscillating gradient diffusion-weighted MR imaging (DWI) is sensitive to short distance diffusion changes at the intracellular level.

PURPOSE

To compare the diagnostic performance of pulsed and oscillating DWI for characterizing hepatocellular nodules in a rat model of hepatic cirrhosis.

STUDY TYPE

Prospective, experimental study.

ANIMAL MODEL

Cirrhosis was induced by weekly intraperitoneal injection of diethylnitrosamine in Wistar rats.

FIELD STRENGTH/SEQUENCE

Ex vivo liver MRI was performed at 7T with T₁ -weighted, T₂ -weighted, pulsed, and oscillating gradient diffusion-weighted sequences.

ASSESSMENT

Apparent diffusion coefficient from pulsed (ADC_pulsed ) and oscillating gradient (ADC_oscillating ) sequences was calculated in 82 nodules identified on the T₁ /T₂ -weighted images and on pathological examination. Two pathologists classified the nodules in three categories: benign (regenerative and low-grade dysplastic nodules), with intermediate malignancy (high-grade dysplastic nodules and early hepatocellular carcinomas) and overtly malignant (progressed hepatocellular carcinomas).

STATISTICAL TESTS

Differences between groups were assessed with Kruskal-Wallis and Mann-Whitney tests.

RESULTS

ADC, mainly ADC_oscillating , increased in the group of nodules with intermediate malignancy (ADC_pulsed : 0.75 ± 0.25 × 10^-3 mm² /s vs. 0.64 ± 0.07 × 10^-3 mm² /s in benign nodules, P = 0.025; ADC_oscillating : 0.81 ± 0.20 × 10^-3 mm² /s vs. 0.65 ± 0.13 × 10^-3 mm² /s, P = 0.0008) and ADC_pulsed decreased in the group of progressed hepatocellular carcinomas (ADC_pulsed : 0.60 ± 0.08 × 10^-3 mm² /s, P = 0.042; ADC_oscillating : 0.68 ± 0.08 × 10^-3 mm² /s, P = 0.1).

DATA CONCLUSION

ADC during hepatocarcinogenesis in rats increased in nodules with intermediate malignancy and decreased in progressed hepatocellular carcinomas. Our results suggest that oscillating gradient DWI is more sensitive to the early steps of hepatocarcinogenesis and might be useful for differentiating between high-grade dysplastic nodules / early hepatocellular carcinomas and regenerating nodules / low-grade dysplastic nodules.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2020;51:1065-1074.

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.

Oscillating gradient diffusion kurtosis imaging of normal and injured mouse brains

Recent advances in diffusion MRI employ multiple diffusion encoding schemes with varying diffusion direction, weighting, and diffusion time to investigate specific microstructural properties in biological tissues. In this study, we examined time-dependent diffusion kurtosis contrast in adult mouse brains and in neonatal mouse brains after hypoxic-ischemic (HI) injury. In vivo diffusion kurtosis maps were acquired with a short diffusion time using an oscillating gradient spin echo (OGSE) sequence at 100 Hz and with a relatively long diffusion time (20 ms) using a pulsed gradient spin echo (PGSE) sequence. In the adult mouse brain, we found that the cortex and hippocampus showed larger differences between OGSE kurtosis and PGSE kurtosis than major white matter tracts. In neonatal mouse brains with unilateral HI injury, the OGSE kurtosis map overall provided stronger edema contrast than the PGSE kurtosis map, and the differences between OGSE and PGSE kurtosis measurements in the edema region reflected heterogeneity of injury. This is the first in vivo study that has demonstrated multi-direction OGSE kurtosis contrasts in the mouse brain. Comparing PGSE and OGSE kurtosis measures may provide additional information on microstructural changes after ischemic stroke.

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.

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.

Timed, sequential administration of paclitaxel improves its cytotoxic effectiveness in a cell culture model

Paclitaxel (taxol) is a chemotherapeutic agent frequently used in combination with other anti-neoplastic drugs. It is most effective during the M phase of the cell-cycle and tends to cause synchronization in malignant cells lines. In this study, we investigated whether timed, sequential treatment based on the cell-cycle characteristics could be exploited to enhance the cytotoxic effect of paclitaxel. We characterized the cell-cycle properties of a rapidly multiplying cell line (Sp2, mouse myeloma cells) by propidium-iodide DNA staining such as the lengths of various cell cycle phases and population duplication time. Based on this we designed a paclitaxel treatment protocol that comprised a primary and a secondary, timed treatment. We found that the first paclitaxel treatment synchronized the cells at the G2/M phase but releasing the block by stopping the treatment allowed a large number of cells to enter the next cell-cycle by a synchronized manner. The second treatment was most effective during the time when these cells approached the next G2/M phase and was least effective when it occurred after the peak time of this next G2/M phase. Moreover, we found that after mixing Sp2 cells with another, significantly slower multiplying cell type (Jurkat human T-cell leukemia) at an initial ratio of 1:1, the ratio of the two different cell types could be influenced by timed sequential paclitaxel treatment at will. Our results demonstrate that knowledge of the cell-cycle parameters of a specific malignant cell type could improve the effectivity of the chemotherapy. Implementing timed chemotherapeutic treatments could increase the cytotoxicity on the malignant cells but also decrease the side-effects since other, non-malignant cell types will have different cell-cycle characteristic and be out of synch during the treatment.

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.

Length of intact plasma membrane determines the diffusion properties of cellular water

Molecular diffusion in a boundary-free medium depends only on the molecular size, the temperature, and medium viscosity. However, the critical determinant of the molecular diffusion property in inhomogeneous biological tissues has not been identified. Here, using an in vitro system and a high-resolution MR imaging technique, we show that the length of the intact plasma membrane is a major determinant of water diffusion in a controlled cellular environment and that the cell perimeter length (CPL) is sufficient to estimate the apparent diffusion coefficient (ADC) of water in any cellular environment in our experimental system (ADC = -0.21 × CPL + 1.10). We used this finding to further explain the different diffusion kinetics of cells that are dying via apoptotic or non-apoptotic cell death pathways exhibiting characteristic changes in size, nuclear and cytoplasmic architectures, and membrane integrity. These results suggest that the ADC value can be used as a potential biomarker for cell death.

Symmetry of the gradient profile as second experimental dimension in the short-time expansion of the apparent diffusion coefficient as measured with NMR diffusometry

The time-dependent apparent diffusion coefficient as measured by pulsed gradient NMR can be used to estimate parameters of porous structures including the surface-to-volume ratio and the mean curvature of pores. In this work, the short-time diffusion limit and in particular the influence of the temporal profile of diffusion gradients on the expansion as proposed by Mitra et al. (1993) is investigated. It is shown that flow-compensated waveforms, i.e. those whose first moment is zero, are blind to the term linear in observation time, which is the term that is proportional to mean curvature and surface permeability. A gradient waveform that smoothly interpolates between flow-compensated and bipolar waveform is proposed and the degree of flow-compensation is used as a second experimental dimension. This two-dimensional ansatz is shown to yield an improved precision when characterizing the confining domain. This technique is demonstrated with simulations and in experiments performed with cylindrical capillaries of 100 μm radius.

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.

Surface-to-volume ratio mapping of tumor microstructure using oscillating gradient diffusion weighted imaging

PURPOSE

To disentangle the free diffusivity (D0 ) and cellular membrane restrictions, by means of their surface-to-volume ratio (S/V), using the frequency-dependence of the diffusion coefficient D(ω), measured in brain tumors in the short diffusion-time regime using oscillating gradients (OGSE).

METHODS

In vivo and ex vivo OGSE experiments were performed on mice bearing the GL261 murine glioma model (n = 10) to identify the relevant time/frequency (t/ω) domain where D(ω) linearly decreases with ω(-1/2) . Parametric maps (S/V, D0 ) are compared with conventional DWI metrics. The impact of frequency range and temperature (20°C versus 37°C) on S/V and D0 is investigated ex vivo.

RESULTS

The validity of the short diffusion-time regime is demonstrated in vivo and ex vivo. Ex vivo measurements confirm that the purely geometric restrictions embodied in S/V are independent from temperature and frequency range, while the temperature dependence of the free diffusivity D0 is similar to that of pure water.

CONCLUSION

Our results suggest that D(ω) in the short diffusion-time regime can be used to uncouple the purely geometric restriction effect, such as S/V, from the intrinsic medium diffusivity properties, and provides a nonempirical and objective way to interpret frequency/time-dependent diffusion changes in tumors in terms of objective biophysical tissue parameters. Magn Reson Med 76:237-247, 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.

Exact analytical results for ADC with oscillating diffusion sensitizing gradients

The apparent diffusion coefficient (ADC) is analyzed for the case of oscillating diffusion sensitizing gradients. Exact analytical expressions are obtained in the high-frequency expansion of the ADC for an arbitrary number of oscillations N. These expressions are universal and valid for arbitrary system geometry. The validity conditions of the high-frequency expansion of ADC are obtained in the framework of a simple 1D model of restricted diffusion. These conditions are shown to be substantially different for cos- and sin-type gradients: for the cos-type gradients, the high-frequency expansion is valid when the period of a single oscillation is smaller than the characteristic diffusion time, the frequency dependence of ADC being practically the same for any N. In contrast, for the sin-type gradients, the high-frequency regime can be achieved only when the total diffusion time is smaller than the characteristic diffusion time, the frequency dependence of ADC being different for different N.

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.

Oscillating and pulsed gradient diffusion magnetic resonance microscopy over an extended b-value range: implications for the characterization of tissue microstructure

Oscillating gradient spin-echo (OGSE) pulse sequences have been proposed for acquiring diffusion data with very short diffusion times, which probe tissue structure at the subcellular scale. OGSE sequences are an alternative to pulsed gradient spin echo measurements, which typically probe longer diffusion times due to gradient limitations. In this investigation, a high-strength (6600 G/cm) gradient designed for small-sample microscopy was used to acquire OGSE and pulsed gradient spin echo data in a rat hippocampal specimen at microscopic resolution. Measurements covered a broad range of diffusion times (TDeff = 1.2-15.0 ms), frequencies (ω = 67-1000 Hz), and b-values (b = 0-3.2 ms/μm2). Variations in apparent diffusion coefficient with frequency and diffusion time provided microstructural information at a scale much smaller than the imaging resolution. For a more direct comparison of the techniques, OGSE and pulsed gradient spin echo data were acquired with similar effective diffusion times. Measurements with similar TDeff were consistent at low b-value (b < 1 ms/μm(2) ), but diverged at higher b-values. Experimental observations suggest that the effective diffusion time can be helpful in the interpretation of low b-value OGSE data. However, caution is required at higher b, where enhanced sensitivity to restriction and exchange render the effective diffusion time an unsuitable representation. Oscillating and pulsed gradient diffusion techniques offer unique, complementary information. In combination, the two methods provide a powerful tool for characterizing complex diffusion within biological tissues.

Surface-to-volume ratio with oscillating gradients

Restrictions to diffusion result in the dispersion of the bulk diffusion coefficient. We derive the exact universal high-frequency behavior of the diffusion coefficient in terms of the surface-to-volume ratio of the restrictions. This frequency dependence can be applied to quantify structure of complex samples with NMR using oscillating field gradients and static-gradient CPMG. We also demonstrate the inter-relations between different equivalent diffusion metrics, and describe how to calculate the effect of restrictions for arbitrary gradient waveforms.

Dependence of temporal diffusion spectra on microstructural properties of biological tissues

The apparent diffusion coefficient (ADC) measured using magnetic resonance imaging methods provides information on microstructural properties of biological tissues, and thus has found applications as a useful biomarker for assessing changes such as those that occur in ischemic stroke and cancer. Conventional pulsed gradient spin echo methods are in widespread use and provide information on, for example, variations in cell density. The oscillating gradient spin echo (OGSE) method has the additional ability to probe diffusion behaviors more readily at short diffusion times, and the temporal diffusion spectrum obtained by the OGSE method provides a unique tool for characterizing tissues over different length scales, including structural features of intracellular spaces. It has previously been reported that several tissue properties can affect ADC measurements significantly, and the precise biophysical mechanisms that account for ADC changes in different situations are still unclear. Those factors may vary in importance depending on the time and length scale over which measurements are made. In the present work, a comprehensive numerical simulation is used to investigate the dependence of the temporal diffusion spectra measured by OGSE methods on different microstructural properties of biological tissues, including cell size, cell membrane permeability, intracellular volume fraction, intranucleus and intracytoplasm diffusion coefficients, nuclear size and T(2) relaxation times. Some unique characteristics of the OGSE method at relatively high frequencies are revealed. The results presented in the paper offer a framework for better understanding possible causes of diffusion changes and may be useful to assist the interpretation of diffusion data from OGSE measurements.