MR cell size imaging with temporal diffusion spectroscopy

Joint estimation of compartment-specific T2 relaxation and tumor microstructure using multi-TE IMPULSED MRI

PURPOSE

This study aims to assess how T2 heterogeneity biases IMPULSED-derived metrics of tissue microstructure in solid tumors and evaluate the potential of estimating multi-compartmental T2 and microstructural parameters simultaneously.

METHODS

This study quantifies the impact of T2 relaxation on IMPULSED-derived microstructural parameters using computer simulations and in vivo multi-TE IMPULSED MRI in five tumor models, including brain, breast, prostate, melanoma, and colon cancer. A comprehensive T2 + IMPULSED method was developed to fit multi-compartmental T2 and microstructural parameters simultaneously. A Bayesian model selection approach was carried out voxel-wisely to determine if the T2 heterogeneity needs to be included in IMPULSED MRI in cancer.

RESULTS

Simulations suggest that T2 heterogeneity has a minor effect on the estimation of d in tissues with intermediate or high cell density, but significantly biases the estimation ofv in $$ {v}_{in} $$ with low cell density. For the in vivo animal experiments, all IMPULSED metrics exceptv in $$ {v}_{in} $$ are statistically independent on TE. For B16 tumors, the IMPULSED-derivedv in $$ {v}_{in} $$ exhibited a notable increase with longer TEs. For MDA-MB-231 tumors, IMPULSED-derivedv in $$ {v}_{in} $$ showed a significant increase with increasing TEs. The T2 + IMPULSED-derivedT 2 in $$ {T}_2^{in} $$ of all five tumor models are consistently smaller thanT 2 ex $$ {T}_2^{ex} $$ .

CONCLUSIONS

The findings from this study highlight two key observations: (i) TE has a negligible impact on IMPULSED-derived cell sizes, and (ii) the TE-dependence of IMPULSED-derived intracellular volume fractions used in T2 + IMPULSED modeling to estimateT 2 in $$ {T}_2^{in} $$ andT 2 ex $$ {T}_2^{ex} $$ . These insights contribute to the ongoing development and refinement of non-invasive MRI techniques for measuring cell sizes.

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.

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.

Diffusion time dependency of extracellular diffusion

PURPOSE

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

METHODS

Model cellular systems containing only extracellular water were used to investigate thet / f $$ t/f $$ dependence ofD 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(D2 O, deuterium oxide) dispersed in regular water (H2 O).D ex $$ {D}_{ex} $$ was obtained over a broadt $$ t $$ range (∼1-1000 ms) and then fit power-law equationsD ex ( t ) = D const + const · t - ϑ t $$ {D}_{ex}(t)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {t}^{-{\vartheta}_t} $$ andD 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 ofD ex $$ {D}_{ex} $$ over the range of diffusion times of most interest in practical dMRI. Previous theoretical predictions are accurate over only limitedt $$ 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 specifict $$ 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 longt $$ 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.

Diffusion Weighted Imaging of the Abdomen and Pelvis: Recent Technical Advances and Clinical Applications

Diffusion weighted imaging (DWI) serves as one of the most important functional magnetic resonance imaging techniques in abdominal and pelvic imaging. It is designed to reflect the diffusion of water molecules and is particularly sensitive to the malignancies. Yet, the limitations of image distortion and artifacts in single-shot DWI may hamper its widespread use in clinical practice. With recent technical advances in DWI, such as simultaneous multi-slice excitation, computed or reduced field-of-view techniques, as well as advanced shimming methods, it is possible to achieve shorter acquisition time, better image quality, and higher robustness in abdominopelvic DWI. This review discussed the recent advances of each DWI approach, and highlighted its future perspectives in abdominal and pelvic imaging, hoping to familiarize physicians and radiologists with the technical improvements in this field and provide future research directions.

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.

Selective Cell Size MRI Differentiates Brain Tumors from Radiation Necrosis

Brain metastasis is a common characteristic of late-stage lung cancers. High doses of targeted radiotherapy can control tumor growth in the brain but can also result in radiotherapy-induced necrosis. Current methods are limited for distinguishing whether new parenchymal lesions following radiotherapy are recurrent tumors or radiotherapy-induced necrosis, but the clinical management of these two classes of lesions differs significantly. Here, we developed, validated, and evaluated a new MRI technique termed selective size imaging using filters via diffusion times (SSIFT) to differentiate brain tumors from radiotherapy necrosis in the brain. This approach generates a signal filter that leverages diffusion time dependence to establish a cell size-weighted map. Computer simulations in silico, cultured cancer cells in vitro, and animals with brain tumors in vivo were used to comprehensively validate the specificity of SSIFT for detecting typical large cancer cells and the ability to differentiate brain tumors from radiotherapy necrosis. SSIFT was also implemented in patients with metastatic brain cancer and radiotherapy necrosis. SSIFT showed high correlation with mean cell sizes in the relevant range of less than 20 μm. The specificity of SSIFT for brain tumors and reduced contrast in other brain etiologies allowed SSIFT to differentiate brain tumors from peritumoral edema and radiotherapy necrosis. In conclusion, this new, cell size-based MRI method provides a unique contrast to differentiate brain tumors from other pathologies in the brain.

SIGNIFICANCE

This work introduces and provides preclinical validation of a new diffusion MRI method that exploits intrinsic differences in cell sizes to distinguish brain tumors and radiotherapy necrosis.

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 pore imaging in the presence of extraporal water

Diffusion-weighted imaging (DWI) is a powerful, non-invasive tool which is widely used in clinical routine. Mostly, apparent diffusion coefficient maps are acquired, which cannot be related directly to cellular structure. More recently it was shown that DWI is able to reconstruct pore shapes using a specialized magnetic field gradient scheme so that cell size distributions may be obtained. So far, artificial systems have been used for experimental demonstration without extraporal signal components and relatively low gradient amplitudes. The aim of this study was to investigate the feasibility of diffusion pore imaging in the presence of extraporal fluids and to develop correction methods for the effects arising from extraporal signal contributions. Monte Carlo simulations and validation experiments on a 14.1 T NMR spectrometer equipped with a dedicated diffusion probe head were performed. Both by using a filter gradient approach suppressing extraporal signal components as well as by using post-processing methods relying on the Gaussian phase approximation, it was possible to reconstruct pore space functions in the presence of extraporal fluids with little to no deviations from the expectations. These results may be a significant step towards application of diffusion pore imaging to biological samples.

Survival Prediction Analysis in Glioblastoma With Diffusion Kurtosis Imaging

SIMPLE SUMMARY

Glioblastoma (GBM) is the most common and aggressive primary brain tumor. Diffusion kurtosis imaging (DKI) has characterized non-Gaussian diffusion behaviors in brain normal tissue and gliomas, but there are very limited efforts in investigating treatment responses of kurtosis in GBM. This study aimed to investigate whether any parameter derived from the DKI is a significant predictor of overall survival (OS). We found that the large mean, 80 and 90 percentile kurtosis values in the contrast enhanced gross tumor volume (Gd-GTV) on post-Gd T1-weighted images pre-RT were significantly associated with reduced OS. In the multivariate Cox model, the mean kurtosis Gd-GTV pre-RT after considering effects of age, extent of surgery, and methylation were significant predictors of OS. In addition, the 80 and 90 percentile kurtosis values in Gd-GTV post RT were significantly associated with progression free survival (PFS). The DKI model demonstrates the potential to predict outcomes in the patients with GBM.

PURPOSE

Non-Gaussian diffusion behaviors in gliomas have been characterized by diffusion kurtosis imaging (DKI). But there are very limited efforts in investigating the kurtosis in glioblastoma (GBM) and its prognostic and predictive values. This study aimed to investigate whether any of the diffusion kurtosis parameters derived from DKI is a significant predictor of overall survival.

METHODS AND MATERIALS

Thirty-three patients with GBM had pre-radiation therapy (RT) and mid-RT diffusion weighted (DW) images. Kurtosis and diffusion coefficient (DC) values in the contrast enhanced gross tumor volume (Gd-GTV) on post-Gd T1 weighted images pre-RT and mid-RT were calculated. Univariate and multivariate Cox models were used to evaluate the DKI parameters and clinical factors for prediction of OS and PFS.

RESULTS

The large mean kurtosis values in the Gd-GTV pre-RT were significantly associated with reduced OS (p = 0.02), but the values at mid-RT were not (p > 0.8). In the multivariate Cox model, the mean kurtosis in the Gd-GTV pre-RT (p = 0.009) was still a significant predictor of OS after adjusting effects of age, O6-Methylguanine-DNA Methyl transferase (MGMT) methylation and extent of resection. In Gd-GTV post-RT, 80 and 90 percentile kurtosis values were significant predictors (p ≤ 0.05) for progression free survival (PFS).

CONCLUSION

The DKI model demonstrates the potential to predict OS and PFS in the patients with GBM. Further development and histopathological validation of the DKI model will warrant its role in clinical management of GBM.