Magnetic resonance imaging of mean cell size in human breast tumors

Revealing membrane integrity and cell size from diffusion kurtosis time dependence

PURPOSE

The nonmonotonic dependence of diffusion kurtosis on diffusion time has been observed in biological tissues, yet its relation to membrane integrity and cellular geometry remains to be clarified. Here we establish and explain the characteristic asymmetric shape of the kurtosis peak. We also derive the relation between the peak timet peak $$ {t}_{\mathrm{peak}} $$ , when kurtosis reaches its maximum, and tissue parameters.

METHODS

The peak shape and its positiont peak $$ {t}_{\mathrm{peak}} $$ qualitatively follow from the adiabatic extension of the Kärger model onto the case of intra-cellular diffusivity time-dependence. This intuition is corroborated by the effective medium theory-based calculation, as well as by Monte Carlo simulations of diffusion and exchange in randomly and densely packed spheres for various values of permeability, cell fractions and sizes, and intrinsic diffusivity.

RESULTS

We establish thatt peak $$ {t}_{\mathrm{peak}} $$ is proportional to the geometric mean of two characteristic time scales: extra-cellular correlation time (determined by cell size) and intra-cellular residence time (determined by membrane permeability). When exchange is barrier-limited, the peak shape approaches a universal scaling form determined by the ratiot / t peak $$ t/{t}_{\mathrm{peak}} $$ .

CONCLUSION

Numerical simulations and theory provide an interpretation of a specific feature of kurtosis time-dependence, offering a potential biomarker for in vivo evaluation of pathology by disentangling the functional (permeability) and structural (cell size) integrity in tissues. This is relevant as the time-dependent diffusion cumulants are sensitive to pathological changes in membrane integrity and cellular structure in diseases, such as ischemic stroke, tumors, and Alzheimer's disease.

Improving Microstructural Estimation in Time-Dependent Diffusion MRI With a Bayesian Method

BACKGROUND

Accurately fitting diffusion-time-dependent diffusion MRI (td-dMRI) models poses challenges due to complex and nonlinear formulas, signal noise, and limited clinical data acquisition.

PURPOSE

Introduce a Bayesian methodology to refine microstructural fitting within the IMPULSED (Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion) model and optimize the prior distribution within the Bayesian framework.

STUDY TYPE

Retrospective.

POPULATION

Involving 69 pediatric patients (median age 6 years, interquartile range [IQR] 3-9 years, 61% male) with 41 low-grade and 28 high-grade gliomas, of which 76.8% were identified within the brainstem or cerebellum.

FIELD STRENGTH/SEQUENCE

3 T, oscillating gradient spin-echo (OGSE) and pulsed gradient spin-echo (PGSE).

ASSESSMENT

The Bayesian method's performance in fitting cell diameter ( d ), intracellular volume fraction (f in ), and extracellular diffusion coefficient (D ex ) was compared against the NLLS method, considering simulated and experimental data. The tumor region-of-interest (ROI) were manually delineated on the b0 images. The diagnostic performance in distinguishing high- and low-grade gliomas was assessed, and fitting accuracy was validated against H&E-stained pathology.

STATISTICAL TESTS

T-test, receiver operating curve (ROC), area under the curve (AUC) and DeLong's test were conducted. Significance considered at P < 0.05.

RESULTS

Bayesian methodology manifested increased accuracy with robust estimates in simulation (RMSE decreased by 29.6%, 40.9%, 13.6%, and STD decreased by 29.2%, 43.5%, and 24.0%, respectively for d ,f in , andD ex compared to NLLS), indicating fewer outliers and reduced error. Diagnostic performance for tumor grade was similar in both methods, however, Bayesian method generated smoother microstructural maps (outliers ratio decreased by 45.3% ± 19.4%) and a marginal enhancement in correlation with H&E staining result (r = 0.721 forf in compared to r = 0.698 using NLLS, P = 0.5764).

DATA CONCLUSION

The proposed Bayesian method substantially enhances the accuracy and robustness of IMPULSED model estimation, suggesting its potential clinical utility in characterizing cellular microstructure.

EVIDENCE LEVEL

3 TECHNICAL EFFICACY: Stage 1.

Influence of Multiband Technique on Temporal Diffusion Spectroscopy and Its Diagnostic Value in Breast Tumors

BACKGROUND

Temporal diffusion spectroscopy (TDS) is a noninvasive diffusion imaging technique used to characterizing cellular microstructures. The influence of multiband (MB) on TDS, particularly in breast tumor imaging remain unknown.

PURPOSE

To investigate the influence of MB on TDS in terms of scanning time, image quality, and quantitative parameters and to assess the diagnostic value of TDS with MB in breast tumors.

STUDY TYPE

Prospective.

POPULATION

Seventy-one women with 71 confirmed lesions.

FIELD STRENGTH/SEQUENCE

3.0 T; oscillating gradient spin-echo (OGSE), OGSE with MB, and pulsed gradient spin-echo, and routine magnetic resonance imaging squences.

ASSESSMENT

TDS with MB was used to assess diagnostic efficacy in differentiating benign and malignant tumors. A comparison of scanning time and image quality was performed in 17 patients. Imaging parameters were analyzed using limited spectrally edited diffusion (IMPULSED) and apparent diffusion coefficient (ADC) values were compared between MB and non-MB protocols. The cell diameter from TDS was compared with histopathological measurements in 21 patients.

STATISTICAL TESTS

Bland-Altman plot, paired t test, Mann-Whitney U test, kappa test, DeLong's test, intraclass correlation coefficient agreement, receiver operating characteristic curve, area under the curve (AUC), and simple linear regression, with statistical significance set at P < 0.05.

RESULTS

The TDS with MB protocol had a shorter average scanning time than that without MB protocol (7 minutes 22 seconds vs. 12 minutes 28 seconds); image quality was improved by reducing image artifacts. Most IMPULSED parameters and ADC values did not significantly differ between the MB and non-MB protocols (P = 0.23, P = 0.17). The IMPULSED parameters of cellularity and intracellular volume fraction achieved the highest AUC values for distinguishing breast tumors (0.865 and 0.821, respectively), surpassing the diagnostic efficiency of conventional ADC-1000 (0.776). The correlation between IMPULSED parameters and microscopic cell size was strong (r = 0.842).

DATA CONCLUSION

The MB technique improved the TDS protocol's efficiency and reduced the image artifacts. TDS parameters correlated with pathological findings and showed good performance in differentiating benign from malignant breast tumors.

PLAIN LANGUAGE SUMMARY

We explored the impact of simultaneous multislice acquisition technology on temporal diffusion spectroscopy (TDS) and whether combining this method could help distinguish benign from malignant breast tumors. Our findings showed that simultaneous multislice acquisition technology shortened the scanning time and improved image quality by reducing motion-related issues. Additionally, measurements of cell size using simultaneous multislice acquisition technology matched well with results from pathology tests. Overall, our study suggests that simultaneous multislice acquisition enhanced TDS could make breast cancer diagnosis more accurate and efficient, offering good advantages compared to conventional imaging methods.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Time-dependent diffusion magnetic resonance imaging for the analysis of parotid gland tumors

BACKGROUND

Different parotid tumors differ in terms of treatment strategies due to their distinct biological behaviors. Time-dependent diffusion magnetic resonance imaging (td-dMRI) can characterize and quantify the cytological indexes, and then aid the differential diagnosis of various tumors. However, the value of td-dMRI in the analysis of parotid gland tumors remains unclear.

PURPOSE

To investigate the value of quantitative parameters derived from td-dMRI in the analysis of parotid gland tumors.

MATERIAL AND METHODS

In total, 39 patients with parotid gland tumors were prospectively enrolled, including 24 patients with polymorphic adenomas (PAs), eight with Warthin's tumors (WTs), and seven with malignant tumors (MTs). Td-dMRI was performed for preoperative evaluation. Intracellular volume fraction (Vin), mean cell size (d), extracellular diffusion coefficient (Dex), and cellularity were obtained based on the Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion model, and compared among the three groups. One-way ANOVA, Kruskal-Wallis test, and receiver operating characteristic (ROC) curve analysis were performed for further statistical analysis as appropriate.

RESULTS

Significant differences were found in all td-dMRI-derived indexes among PAs, WTs, and MTs (all P < 0.05). Vin was the sole parameter with significant differences for all sub-group comparisons (PAs vs. WTs, P < 0.001; PAs vs. MTs, P = 0.031; WTs vs. MTs, P = 0.047). With Vin values of 0.267, 0.231, and 0.260 as threshold, respectively, optimal performance levels were obtained for differentiating PAs from WTs (area under the ROC curve [AUC]=0.932, sensitivity=0.917, and specificity=0.875), PAs from MTs (AUC=0.744, sensitivity=0.833, and specificity=0.714), and WTs from MTs (AUC=0.750, sensitivity=0.875, and specificity=0.714).

CONCLUSION

Microstructural parameters derived from td-dMRI, especially Vin, might be promising imaging biomarkers for characterizing parotid gland tumors.

Accelerated MR cell size imaging through parallel acquisition technique (PAT) and simultaneous multi-slice (SMS) with local principal component analysis (LPCA) enhancement

Microstructural parameters are essential in tumor research, aiding in the understanding tumor pathogenesis, grading, and therapeutic efficacy. The imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED) model is the most widely used MR cell size imaging technique, demonstrating success in measuring microstructural parameters of solid tumors in vivo. However, its clinical application is limited by the longer scan times required for both pulsed gradient spin-echo (PGSE) and multiple oscillating gradient spin-echo (OGSE) acquisitions across a range of b-values, which can be burdensome for patients and disrupt clinical workflows. In this work, we propose and evaluate an accelerating method that integrates parallel acquisition technique (PAT) and simultaneous multi-slice (SMS) with local principal component analysis (LPCA) denoising to reduce scan times while maintaining image quality in MR cell size imaging. PGSE and OGSE (25 Hz, 50 Hz) images were acquired using P2S2 (PAT2-SMS2), P2S3 (PAT2-SMS3), and P3S3 (PAT3-SMS3) configurations, incorporating LPCA denoising, and compared to standard P2 (PAT2-SMS1) in healthy volunteers and brain tumor patients at 3 T. Additionally, clinical feasibility was further assessed through qualitative and quantitative evaluations. Qualitative assessment, conducted by two radiologists using a 5-point Likert scale, and quantitative analysis, including noise estimation, apparent diffusion coefficient (ADC) calculation, and estimation of microstructural parameters-cell diameter (Dmean), intracellular volume fraction (Vin), and extracellular diffusivity (Dex), were performed. Overall, the integration of PAT and SMS techniques reduces acquisition time by approximately 60 % compared to standard P2 acceleration, while maintaining comparable image quality and structural fidelity with LPCA denoising.

Time-dependent diffusion MRI and kinetic heterogeneity as potential imaging biomarkers for diagnosing suspicious breast lesions with 3.0-T breast MRI

PURPOSE

This study aimed to evaluate the diagnostic efficacy of time-dependent diffusion magnetic resonance imaging (td-dMRI) and dynamic contrast-enhanced MRI (DCE-MRI)-based kinetic heterogeneity in differentiating suspicious breast lesions (categorised as Breast Imaging Reporting and Data System 4 or 5).

METHODS

This prospective study included 51 females with suspicious breast lesions who underwent preoperative breast MRI, including DCE-MRI and td-dMRI. Six kinetic parameters, namely peak, persistent, plateau, washout component, predominant curve type, and heterogeneity, were extracted from the DCE series using MATLAB and SPM software. The td-dMRI data were analysed using the JOINT model to obtain five microstructural parameters and apparent diffusion coefficient at 50 ms (ADC50ms). Chi-square or Fisher's exact test and the Mann-Whitney U test were used to compare these parameters between benign and malignant breast lesions. Univariate and multivariate logistic regression analyses with forward stepwise covariate selection were performed to identify significant clinical and radiologic variables. Differential diagnostic performance was evaluated using receiver operating characteristic curves and logistic regression analyses.

RESULTS

For td-dMRI-derived parameters, the values of fin and cellularity were significantly higher in malignant breast lesions compared to benign lesions (P = 0.001 and P<0.001, respectively), while ADC50ms was significantly lower in malignant lesions (P = 0.001). In the kinetic heterogeneity analysis, the washout component was higher in malignant lesions compared to benign lesions (P = 0.003). When combining significant td-dMRI and kinetic heterogeneity parameters, the area under the curve (AUC) value was 0.875, with an accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) of 82.69 %, 86.11 %, 75.00 %, 88.57 %, and 70.59 %, respectively. Notably, margin and kinetic pattern emerged as independent predictors of malignant breast lesions (P = 0.019 and 0.006, respectively). Furthermore, incorporating these two clinical-radiologic characteristics further enhanced diagnostic accuracy, yielding an AUC of 0.969, with accuracy, sensitivity, specificity, PPV, and NPV improving to 90.38 %, 86.11 %, 100 %, 100 %, and 76.19 %, respectively.

CONCLUSIONS

Kinetic heterogeneity- and td-dMRI-derived parameters are potentially non-invasive biomarkers for distinguishing suspicious breast lesions.

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.

Diffusion time effects over the adult lifespan indicates persistent zone-specific microstructural alterations in the human prostate with aging

PURPOSE

The purpose of this study was to investigate microstructural changes in the aging adult prostate by comparing the effects of varying diffusion times using diffusion MRI, and to provide an age-related benchmark for future prostate cancer studies.

METHODS

The prostates of normal male volunteers (n = 70, 19-69 years) were scanned at 3 T with an oscillating gradient spin echo (OGSE: 6 ms), pulsed gradient spin echo (PGSE: 40 ms) and pulsed gradient stimulated echo (PGSTE: 100 ms), and anatomical T2-weighted image. Volume and mean diffusivity (MD) were measured in the peripheral (PZ) and transition zones (TZ), which were assessed versus age.

RESULTS

PZ and TZ showed quadratic age trajectories for all diffusion scans, with MD decreasing from 19 years to a minimum ˜30-40 years followed by a greater increase at older ages. Short (OGSE) and medium (PGSE) diffusion time MD had similar age trajectories, whereas long diffusion time (PGSTE) MD was significantly lower, particularly in PZ (22%). MD difference (∆MD) of OGSE-PGSTE and PGSE-PGSTE showed significant positive linear correlations with age for both PZ (larger slope) and TZ, resulting in ˜3.3x (PZ) and 1.8x (TZ) greater ∆MD from 19 to 69 years. MD and ∆MD versus age relationships differed from volume, which conversely had greater proportional growth in TZ than PZ.

CONCLUSION

The diffusion time effects suggest age-related microstructural changes consistent with development of persistently larger cell dimensions mainly in the prostate peripheral zone over the adult lifespan. This normative data can be used for comparison to prostate cancer factoring in age.

From Brownian motion to virtual biopsy: a historical perspective from 40 years of diffusion MRI

Diffusion MRI was introduced in 1985, showing how the diffusive motion of molecules, especially water, could be spatially encoded with MRI to produce images revealing the underlying structure of biologic tissues at a microscopic scale. Diffusion is one of several Intravoxel Incoherent Motions (IVIM) accessible to MRI together with blood microcirculation. Diffusion imaging first revolutionized the management of acute cerebral ischemia by allowing diagnosis at an acute stage when therapies can still work, saving the outcomes of many patients. Since then, the field of diffusion imaging has expanded to the whole body, with broad applications in both clinical and research settings, providing insights into tissue integrity, structural and functional abnormalities from the hindered diffusive movement of water molecules in tissues. Diffusion imaging is particularly used to manage many neurologic disorders and in oncology for detecting and classifying cancer lesions, as well as monitoring treatment response at an early stage. The second major impact of diffusion imaging concerns the wiring of the brain (Diffusion Tensor Imaging, DTI), allowing to obtain from the anisotropic movement of water molecules in the brain white-matter images in 3 dimensions of the brain connections making up the Connectome. DTI has opened up new avenues of clinical diagnosis and research to investigate brain diseases, neurogenesis and aging, with a rapidly extending field of application in psychiatry, revealing how mental illnesses could be seen as Connectome spacetime disorders. Adding that water diffusion is closely associated to neuronal activity, as shown from diffusion fMRI, one may consider that diffusion MRI is ideally suited to investigate both brain structure and function. This article retraces the early days and milestones of diffusion MRI which spawned over 40 years, showing how diffusion MRI emerged and expanded in the research and clinical fields, up to become a pillar of modern clinical imaging.

Initial experience in implementing quantitative DCE-MRI to predict breast cancer therapy response in a multi-center and multi-vendor platform setting

Quantitative dynamic contrast-enhanced (DCE) MRI as a promising method for the prediction of breast cancer response to neoadjuvant chemotherapy (NAC) has been demonstrated mostly in single-center and single-vendor platform studies. This preliminary study reports the initial experience in implementing quantitative breast DCE-MRI in multi-center (MC) and multi-vendor platform (MP) settings to predict NAC response. MRI data, including B1 mapping, variable flip angle (VFA) measurements of native tissue R1 (R1,0), and DCE-MRI, were acquired during NAC at three sites using 3T systems with Siemens, Philips, and GE platforms, respectively. High spatiotemporal resolution DCE-MRI was performed using similar vendor product sequences with k-space undersampling during acquisition and view sharing during reconstruction. A breast phantom was used for quality assurance/quality control (QA/QC) across sites. The Tofts model (TM) and shutter-speed model (SSM) were used for pharmacokinetic (PK) analysis of the DCE data. Additionally, tumor region of interest (ROI)- vs. voxel-based analyses in combination with the use of VFA-measured R1,0 vs. fixed, literature-reported R1,0 were investigated to determine the optimal analysis approach. Results from 15 patients who completed the study are reported. Voxel-based PK analysis using fixed R1,0 was deemed the optimal approach, which allowed the inclusion of data from one vendor platform where VFA measurements produced ≥100% overestimation of R1,0. The semi-quantitative signal enhancement ratio (SER) and quantitative PK parameters outperformed the tumor longest diameter (LD) in the prediction of pathologic complete response (pCR) vs. non-pCR after the first NAC cycle, whereas Ktrans consistently provided more accurate predictions than both SER and LD after the first NAC cycle and at the NAC midpoint. Both TM and SSM Ktrans and kep were excellent predictors of response at the NAC midpoint with ROC AUC >0.90, while the SSM parameters (AUC ≥0.80) performed better than their TM counterparts (AUC <0.80) after the first NAC cycle. The initial experience of this ongoing study indicates the importance of QA/QC using a phantom and suggests that deploying voxel-based PK analysis using a fixed R1,0 may mitigate random errors from R1,0 measurements across platforms and potentially eliminate the need for B1 and VFA acquisitions in MC and MP trials.

Advanced breast diffusion-weighted imaging: what are the next steps? A proposal from the EUSOBI International Breast Diffusion-weighted Imaging working group

OBJECTIVES

This study by the EUSOBI International Breast Diffusion-weighted Imaging (DWI) working group aimed to evaluate the current and future applications of advanced DWI in breast imaging.

METHODS

A literature search and a comprehensive survey of EUSOBI members to explore the clinical use and potential of advanced DWI techniques and a literature search were involved. Advanced DWI approaches such as intravoxel incoherent motion (IVIM), diffusion kurtosis imaging (DKI), and diffusion tensor imaging (DTI) were assessed for their current status and challenges in clinical implementation.

RESULTS

Although a literature search revealed an increasing number of publications and growing academic interest in advanced DWI, the survey revealed limited adoption of advanced DWI techniques among EUSOBI members, with 32% using IVIM models, 17% using non-Gaussian diffusion techniques for kurtosis analysis, and only 8% using DTI. A variety of DWI techniques are used, with IVIM being the most popular, but less than half use it, suggesting that the study identified a gap between the potential benefits of advanced DWI and its actual use in clinical practice.

CONCLUSION

The findings highlight the need for further research, standardization and simplification to transition advanced DWI from a research tool to regular practice in breast imaging. The study concludes with guidelines and recommendations for future research directions and clinical implementation, emphasizing the importance of interdisciplinary collaboration in this field to improve breast cancer diagnosis and treatment.

CLINICAL RELEVANCE STATEMENT

Advanced DWI in breast imaging, while currently in limited clinical use, offers promising improvements in diagnosis, staging, and treatment monitoring, highlighting the need for standardized protocols, accessible software, and collaborative approaches to promote its broader integration into routine clinical practice.

KEY POINTS

Increasing number of publications on advanced DWI over the last decade indicates growing research interest. EUSOBI survey shows that advanced DWI is used primarily in research, not extensively in clinical practice. More research and standardization are needed to integrate advanced DWI into routine breast imaging practice.

Advanced Diffusion-Weighted MRI for Cancer Microstructure Assessment in Body Imaging, and Its Relationship With Histology

Diffusion-weighted magnetic resonance imaging (DW-MRI) aims to disentangle multiple biological signal sources in each imaging voxel, enabling the computation of innovative maps of tissue microstructure. DW-MRI model development has been dominated by brain applications. More recently, advanced methods with high fidelity to histology are gaining momentum in other contexts, for example, in oncological applications of body imaging, where new biomarkers are urgently needed. The objective of this article is to review the state-of-the-art of DW-MRI in body imaging (ie, not including the nervous system) in oncology, and to analyze its value as compared to reference colocalized histology measurements, given that demonstrating the histological validity of any new DW-MRI method is essential. In this article, we review the current landscape of DW-MRI techniques that extend standard apparent diffusion coefficient (ADC), describing their acquisition protocols, signal models, fitting settings, microstructural parameters, and relationship with histology. Preclinical, clinical, and in/ex vivo studies were included. The most used techniques were intravoxel incoherent motion (IVIM; 36.3% of used techniques), diffusion kurtosis imaging (DKI; 16.7%), vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT; 13.3%), and imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED; 11.7%). Another notable category of techniques relates to innovative b-tensor diffusion encoding or joint diffusion-relaxometry. The reviewed approaches provide histologically meaningful indices of cancer microstructure (eg, vascularization/cellularity) which, while not necessarily accurate numerically, may still provide useful sensitivity to microscopic pathological processes. Future work of the community should focus on improving the inter-/intra-scanner robustness, and on assessing histological validity in broader contexts. LEVEL OF EVIDENCE: NA TECHNICAL EFFICACY: Stage 2.

Time-Dependent Diffusion MRI Helps Predict Molecular Subtypes and Treatment Response to Neoadjuvant Chemotherapy in Breast Cancer

Background Time-dependent diffusion MRI has the potential to help characterize tumor cell properties; however, to the knowledge of the authors, its usefulness for breast cancer diagnosis and prognostic evaluation is unknown. Purpose To investigate the clinical value of time-dependent diffusion MRI-based microstructural mapping for noninvasive prediction of molecular subtypes and pathologic complete response (pCR) in participants with breast cancer. Materials and Methods Participants with invasive breast cancer who underwent pretreatment with time-dependent diffusion MRI between February 2021 and May 2023 were prospectively enrolled. Four microstructural parameters were estimated using the IMPULSED method (a form of time-dependent diffusion MRI), along with three apparent diffusion coefficient (ADC) measurements and a relative ADC diffusion-weighted imaging parameter. Multivariable logistic regression analysis was used to identify parameters associated with each molecular subtype and pCR. A predictive model based on associated parameters was constructed, and its performance was assessed using the area under the receiver operating characteristic curve (AUC) and compared by using the DeLong test. The time-dependent diffusion MRI parameters were validated based on correlation with pathologic measurements. Results The analysis included 408 participants with breast cancer (mean age, 51.9 years ± 9.1 [SD]). Of these, 221 participants were administered neoadjuvant chemotherapy and 54 (24.4%) achieved pCR. The time-dependent diffusion MRI parameters showed reasonable performance in helping to identify luminal A (AUC, 0.70), luminal B (AUC, 0.78), and triple-negative breast cancer (AUC, 0.72) subtypes and high performance for human epidermal growth factor receptor 2 (HER2)-enriched breast cancer (AUC, 0.85), outperforming ADC measurements (all P < .05). Progesterone receptor status (odds ratio [OR], 0.08; P = .02), HER2 status (OR, 3.36; P = .009), and the cellularity index (OR, 0.01; P = .02) were independently associated with the odds of achieving pCR. The combined model showed high performance for predicting pCR (AUC, 0.88), outperforming ADC measurements and the clinical-pathologic model (AUC, 0.73 and 0.79, respectively; P < .001). The time-dependent diffusion MRI-estimated parameters correlated well with the pathologic measurements (n = 100; r = 0.67-0.81; P < .001). Conclusion Time-dependent diffusion MRI-based microstructural mapping was an effective method for helping to predict molecular subtypes and pCR to neoadjuvant chemotherapy in participants with breast cancer. © RSNA, 2024 Supplemental material is available for this article. See also the editorial by Partridge and Xu in this issue.

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

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

Normal Brain and Brain Tumor ADC: Changes Resulting From Variation of Diffusion Time and/or Echo Time in Pulsed-Gradient Spin Echo Diffusion Imaging

OBJECTIVES

Increasing gradient performance on modern magnetic resonance imaging scanners has profoundly reduced the attainable diffusion and echo times for clinically available pulsed-gradient spin echo (PGSE) sequences. This study investigated how this may impact the measured apparent diffusion coefficient (ADC), which is considered an important diagnostic marker for differentiation between normal and abnormal brain tissue and for therapeutic follow-up.

MATERIALS AND METHODS

Diffusion time and echo time dependence of the ADC were evaluated on a high-performance 3 T magnetic resonance imaging scanner. Diffusion PGSE brain scans were performed in 10 healthy volunteers and in 10 brain tumor patients using diffusion times of 16, 40, and 70 ms, echo times of 60, 75, and 104 ms at 3 b-values (0, 100, and 1000 s/mm 2 ), and a maximum gradient amplitude of 68 mT/m. A low gradient performance system was also emulated by reducing the diffusion encoding gradient amplitude to 19 mT/m. In healthy subjects, the ADC was measured in 6 deep gray matter regions and in 6 white matter regions. In patients, the ADC was measured in the solid part of the tumor.

RESULTS

With increasing diffusion time, a small but significant ADC increase of up to 2.5% was observed for 6 aggregate deep gray matter structures. With increasing echo time or reduced gradient performance, a small but significant ADC decrease of up to 2.6% was observed for 6 aggregate white matter structures. In tumors, diffusion time-related ADC changes were inconsistent without clear trend. For tumors with diffusivity above 1.0 μm 2 /ms, with prolonged echo time, there was a pronounced ADC increase of up to 12%. Meanwhile, for tumors with diffusivity at or below 1.0 μm 2 /ms, no change or a reduction was observed. Similar results were observed for gradient performance reduction, with an increase of up to 21%. The coefficient of variation determined in repeat experiments was 2.4%.

CONCLUSIONS

For PGSE and the explored parameter range, normal tissue ADC changes seem negligible. Meanwhile, observed tumor ADC changes can be relevant if ADC is used as a quantitative biomarker and not merely assessed by visual inspection. This highlights the importance of reporting all pertinent timing parameters in ADC studies and of considering these effects when building scan protocols for use in multicenter investigations.

Time-dependent diffusion MRI-based microstructural mapping for differentiating high-grade serous ovarian cancer from serous borderline ovarian tumor

PURPOSE

To investigate the value of microstructural characteristics derived from time-dependent diffusion MRI in distinguishing high-grade serous ovarian cancer (HGSOC) from serous borderline ovarian tumor (SBOT) and the associations of immunohistochemical markers with microstructural features.

METHODS

Totally 34 HGSOC and 12 SBOT cases who received preoperative pelvic MRI were retrospectively included in this study. Two radiologists delineated the tumors to obtain the regions of interest (ROIs). Time-dependent diffusion MRI signals were fitted by the IMPULSED (imaging microstructural parameters using limited spectrally edited diffusion) model, to extract microstructural parameters, including fraction of the intracellular component (fin), cell diameter (d), cellularity and extracellular diffusivity (Dex). Apparent diffusion coefficient (ADC) values were obtained from standard diffusion-weighted imaging (DWI). The parameters of HGSOCs and SBOTs were compared, and the diagnostic performance was evaluated. The associations of microstructural indexes with immunopathological parameters were assessed, including Ki-67, P53, Pax-8, ER and PR.

RESULTS

In this study, fin, cellularity, Dex and ADC had good diagnostic performance levels in differentiating HGSOC from SBOT, with AUCs of 0.936, 0.909, 0.902 and 0.914, respectively. There were no significant differences in diagnostic performance among these parameters. Spearman analysis revealed in the HGSOC group, cellularity had a significant positive correlation with P53 expression (P = 0.028, r = 0.389) and Dex had a significant positive correlation with Pax-8 expression (P = 0.018, r = 0.415). ICC showed excellent agreement for all parameters.

CONCLUSION

Time-dependent diffusion MRI had value in evaluating the microstructures of HGSOC and SBOT and could discriminate between these tumors.

Nonparametric distributions of tensor-valued Lorentzian diffusion spectra for model-free data inversion in multidimensional diffusion MRI

Magnetic resonance imaging (MRI) is the method of choice for noninvasive studies of micrometer-scale structures in biological tissues via their effects on the time- and frequency-dependent (restricted) and anisotropic self-diffusion of water. While new designs of time-dependent magnetic field gradient waveforms have enabled disambiguation between different aspects of translational motion that are convolved in traditional MRI methods relying on single pairs of field gradient pulses, data analysis for complex heterogeneous materials remains a challenge. Here, we propose and demonstrate nonparametric distributions of tensor-valued Lorentzian diffusion spectra, or "D(ω) distributions," as a general representation with sufficient flexibility to describe the MRI signal response from a wide range of model systems and biological tissues investigated with modulated gradient waveforms separating and correlating the effects of restricted and anisotropic diffusion.

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.

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.

In Vivo and Post-mortem Comparisons of IVIM/Time-dependent Diffusion MR Imaging Parameters in Melanoma and Breast Cancer Xenograft Models

PURPOSE

We aimed to investigate the changes in intravoxel incoherent motion (IVIM) and diffusion parameters between in vivo and post-mortem conditions and the time dependency of these parameters using two different mouse tumor models with different vessel lumen sizes.

METHODS

Six B16 and six MDA-MB-231 xenograft mice were scanned using 7 Tesla MRI under both in vivo/post-mortem conditions. Diffusion weighted imaging with 17 b-values (0-3000 s/mm2) were obtained at two diffusion times (9 and 27.6 ms). The shifted apparent diffusion coefficient (sADC) using 2 b-values (200 and 1500 s/mm2), non-Gaussian diffusion and IVIM parameters (ADC0, K, fIVIM) were estimated at each of the diffusion times. The results were evaluated by repeated measures two-way analysis of variance and post hoc Bonferroni test.

RESULTS

In B16 tumors, fIVIM significantly decreased with post-mortem conditions (from 12.6 ± 6.5% to 5.2 ± 1.9%, P < 0.05 at long diffusion time; from 11.0 ± 2.4% to 4.6 ± 2.7%, P < 0.05 at short diffusion time). In MDA-MB-231 tumors, fIVIM also significantly decreased (from 8.8 ± 3.8% to 2.6 ± 1.1%, P < 0.05 at long; from 7.9 ± 5.4% to 2.9 ± 1.1%, P < 0.05 at short). No diffusion time dependency was observed (P = 0.59 in B16 and P = 0.77 in MDA-MB-231). The sADC and ADC0 values tended to decrease and the K value tended to increase after sacrificing and when increasing the diffusion time.

CONCLUSION

The fIVIM values dropped after sacrificing, confirming that IVIM MRI is a promising quantitative parameter to evaluate blood microcirculation. The presence of residual post-mortem fIVIM values suggested that the influence of water molecule diffusion in the blood lumen may contribute to the IVIM effect. Diffusion MRI parameter's time dependency and those changes after sacrificing could possibly provide additional insights into diffusion hindrance mechanisms.

Developing an AI-empowered head-only ultra-high-performance gradient MRI system for high spatiotemporal neuroimaging

Recent advances in neuroscience requires high-resolution MRI to decipher the structural and functional details of the brain. Developing a high-performance gradient system is an ongoing effort in the field to facilitate high spatial and temporal encoding. Here, we proposed a head-only gradient system NeuroFrontier, dedicated for neuroimaging with an ultra-high gradient strength of 650 mT/m and 600 T/m/s. The proposed system features in 1) ultra-high power of 7MW achieved by running two gradient power amplifiers using a novel paralleling method; 2) a force/torque balanced gradient coil design with a two-step mechanical structure that allows high-efficiency and flexible optimization of the peripheral nerve stimulation; 3) a high-density integrated RF system that is miniaturized and customized for the head-only system; 4) an AI-empowered compressed sensing technique that enables ultra-fast acquisition of high-resolution images and AI-based acceleration in q-t space for diffusion MRI (dMRI); and 5) a prospective head motion correction technique that effectively corrects motion artifacts in real-time with 3D optical tracking. We demonstrated the potential advantages of the proposed system in imaging resolution, speed, and signal-to-noise ratio for 3D structural MRI (sMRI), functional MRI (fMRI) and dMRI in neuroscience applications of submillimeter layer-specific fMRI and dMRI. We also illustrated the unique strength of this system for dMRI-based microstructural mapping, e.g., enhanced lesion contrast at short diffusion-times or high b-values, and improved estimation accuracy for cellular microstructures using diffusion-time-dependent dMRI or for neurite microstructures using q-space approaches.

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 dependent diffusion MRI: effect of diffusion-time on microstructural mapping and prediction of prognostic features in breast cancer

OBJECTIVES

This study aimed to evaluate the effect of achievable td on the accuracy of microstructural mapping based on simulation and patient experiments, and investigate the feasibility of td-dMRI in distinguishing prognostic factors in breast cancer patients.

METHODS

Simulation was performed using different td settings. Patients with breast cancer were enrolled prospectively between November 2020 and January 2021, who underwent oscillating and pulsed gradient encoded dMRI on a 3-T scanner using short-/long-td protocol with oscillating frequency up to 50/33 Hz. Data were fitted with a two-compartment model to estimate cell diameter (d), intracellular fraction (fin), and diffusivities. Estimated microstructural markers were used to differentiate immunohistochemical receptor status and the presence of lymph node (LN), which were correlated with histopathological measurements.

RESULTS

Simulation results showed that d fitted from the short-td protocol significantly reduced estimation error than those from long-td (2.07 ± 1.51% versus 3.05 ± 1.92%, p < 0.0001) while the estimation error of fin was robust to different protocols. Among a total of 37 breast cancer patients, the estimated d was significantly higher in HER2-positive and LN-positive (p < 0.05) groups compared to their negative counterparts only using the short-td protocol. Histopathological validation in a subset of 6 patients with whole slide images showed the estimated d was highly correlated with measurements from H&E staining (r = 0.84, p = 0.03) only using the short-td protocol.

CONCLUSIONS

The results indicated the necessity of short-td for accurate microstructural mapping in breast cancer. The current td-dMRI with a total acquisition time of 4.5 min showed its potential in the diagnosis of breast cancer.

KEY POINTS

• Short td is important for accurate microstructural mapping in breast cancer using the td-dMRI technique, based on simulation and histological validation. • The 4.5-min td-dMRI protocol showed potential clinical value for breast cancer, given the difference in cell diameter between HER2/LN positive and negative groups.

Histological and molecular classifications of pediatric glioma with time-dependent diffusion MRI-based microstructural mapping

BACKGROUND

Gliomas are the most common type of central nervous system tumors in children, and the combination of histological and molecular classification is essential for prognosis and treatment. Here, we proposed a newly developed microstructural mapping technique based on diffusion-time-dependent diffusion MRI td-dMRI theory to quantify tumor cell properties and tested these microstructural markers in identifying histological grade and molecular alteration of H3K27.

METHODS

This prospective study included 69 pediatric glioma patients aged 6.14 ± 3.25 years old, who underwent td-dMRI with pulsed and oscillating gradient diffusion sequences on a 3T scanner. dMRI data acquired at varying tds were fitted into a 2-compartment microstructural model to obtain intracellular fraction (fin), cell diameter, cellularity, etc. Apparent diffusivity coefficient (ADC) and T1 and T2 relaxation times were also obtained. H&E stained histology was used to validate the estimated microstructural properties.

RESULTS

For histological classification of low- and high-grade pediatric gliomas, the cellularity index achieved the highest area under the receiver-operating-curve (AUC) of 0.911 among all markers, while ADC, T1, and T2 showed AUCs of 0.906, 0.885, and 0.886. For molecular classification of H3K27-altered glioma in 39 midline glioma patients, cell diameter showed the highest discriminant power with an AUC of 0.918, and the combination of cell diameter and extracellular diffusivity further improved AUC to 0.929. The td-dMRI estimated fin correlated well with the histological ground truth with r = 0.7.

CONCLUSIONS

The td-dMRI-based microstructural properties outperformed routine MRI measurements in diagnosing pediatric gliomas, and the different microstructural features showed complementary strength in histological and molecular classifications.

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.

Evaluation of a continuous-time random-walk diffusion model for the differentiation of malignant and benign breast lesions and its association with Ki-67 expression

The purpose of the current study was to evaluate the performance of a continuous-time random-walk (CTRW) diffusion model for differentiating malignant and benign breast lesions and to consider the potential association between CTRW parameters and the Ki-67 expression. Sixty-four patients (46.2 ± 11.4 years) with breast lesions (29 malignant and 35 benign) were evaluated with the CTRW model, intravoxel incoherent motion model, and diffusion-weighted imaging. Echo planar diffusion-weighted imaging was conducted using 13 b-values (0-3000 s/mm² ). Three CTRW model parameters, including an anomalous diffusion coefficient D_m , and two parameters related to temporal and spatial diffusion heterogeneity, α and β, respectively, were obtained, and had MRI b-values of 0-3000 s/mm² . Receiver operating characteristic (ROC) analysis was conducted to determine the sensitivity, specificity, and diagnostic accuracy of CTRW parameters for differentiating malignant from benign breast lesions. In malignant breast lesions, the CTRW parameters D_m , α, and β were significantly lower than the corresponding parameters of benign breast lesions. In the malignant breast lesion group, the CTRW parameter D_m was significantly lower in high Ki-67 expression than in low Ki-67 expression. In ROC analysis, the combination of CTRW parameters (D_m , α, β) demonstrated the highest area under the curve value (0.985) and diagnostic accuracy (94.23%) in differentiating malignant and benign breast lesions. The CTRW model effectively differentiated malignant from benign breast lesions. The CTRW diffusion model offers a new way for noninvasive assessment of breast malignancy and better understanding of the proliferation of malignant lesions.

Profiling specific cell populations within the inflammatory tumor microenvironment by oscillating-gradient diffusion-weighted MRI

BACKGROUND

The inflammatory tumor microenvironment (TME) is formed by various immune cells, being closely associated with tumorigenesis. Especially, the interaction between tumor-infiltrating T-cells and macrophages has a crucial impact on tumor progression and metastatic spread. The purpose of this study was to investigate whether oscillating-gradient diffusion-weighted MRI (OGSE-DWI) enables a cell size-based discrimination between different cell populations of the TME.

METHODS

Sine-shaped OGSE-DWI was combined with the Imaging Microstructural Parameters Using Limited Spectrally Edited Diffusion (IMPULSED) approach to measure microscale diffusion distances, here relating to cell sizes. The accuracy of IMPULSED-derived cell radii was evaluated using in vitro spheroid models, consisting of either pure cancer cells, macrophages, or T-cells. Subsequently, in vivo experiments aimed to assess changes within the TME and its specific immune cell composition in syngeneic murine breast cancer models with divergent degrees of malignancy (4T1, 67NR) during tumor progression, clodronate liposome-mediated depletion of macrophages, and immune checkpoint inhibitor (ICI) treatment. Ex vivo analysis of IMPULSED-derived cell radii was conducted by immunohistochemical wheat germ agglutinin staining of cell membranes, while intratumoral immune cell composition was analyzed by CD3 and F4/80 co-staining.

RESULTS

OGSE-DWI detected mean cell radii of 8.8±1.3 µm for 4T1, 8.2±1.4 µm for 67NR, 13.0±1.7 for macrophage, and 3.8±1.8 µm for T-cell spheroids. While T-cell infiltration during progression of 4T1 tumors was observed by decreasing mean cell radii from 9.7±1.0 to 5.0±1.5 µm, increasing amount of intratumoral macrophages during progression of 67NR tumors resulted in increasing mean cell radii from 8.9±1.2 to 12.5±1.1 µm. After macrophage depletion, mean cell radii decreased from 6.3±1.7 to 4.4±0.5 µm. T-cell infiltration after ICI treatment was captured by decreasing mean cell radii in both tumor models, with more pronounced effects in the 67NR tumor model.

CONCLUSIONS

OGSE-DWI provides a versatile tool for non-invasive profiling of the inflammatory TME by assessing the dominating cell type T-cells or macrophages.

Inversion-Recovery-Prepared Oscillating Gradient Sequence Improves Diffusion-Time Dependency Measurements in the Human Brain

BACKGROUND

Oscillating gradient diffusion MRI (dMRI) enables measurements at a short diffusion-time (t_d ), but it is challenging for clinical systems. Particularly, the low b-value and low resolution may give rise to cerebrospinal fluid (CSF) contamination.

PURPOSE

To assess the effect of CSF partial volume on t_d -dMRI measurements and efficacy of inversion-recovery (IR) prepared oscillating and pulsed gradient dMRI sequence to improve t_d -dMRI measurements in the human brain.

STUDY TYPE

Prospective.

SUBJECTS

Ten normal volunteers and six glioma patients.

FIELD STRENGTH/SEQUENCE

A 3 T; three-dimensional (3D) IR-prepared oscillating gradient-prepared gradient spin-echo (GRASE) and two-dimensional (2D) IR-prepared oscillating gradient echo-planar imaging (EPI) sequences.

ASSESSMENT

We assessed the t_d -dependent patterns of apparent diffusion coefficient (ADC) in several gray and white matter structures, including the hippocampal subfields (head, body, and tail), cortical gray matter, thalamus, and posterior white matter in normal volunteers. Pulsed gradient (0 Hz) and oscillating gradients at frequencies of 20 Hz, 40 Hz, and 60 Hz dMRI were acquired with GRASE and EPI sequences with or without the IR module. We also tested the t_d -dependency patterns in glioma patients using the EPI sequence with or without the IR module.

STATISTICAL TESTS

The differences in ADC across the different t_d s were compared by one-way ANOVA followed by post hoc pairwise t-tests with Bonferroni correction.

RESULTS

In the healthy subjects, brain regions that were possibly contaminated by CSF signals, such as the hippocampus (head, body, and tail) and cortical gray matter, t_d -dependent ADC changes were only significant with the IR-prepared 2D and 3D sequences but not with the non-IR sequences. In brain glioblastomas patients, significantly higher t_d -dependence was observed in the tumor region with the IR module than that without IR (slope = 0.0196 μm² /msec² vs. 0.0034 μm² /msec² ).

CONCLUSION

The IR-prepared sequence effectively suppressed the CSF partial volume effect and significantly improved the t_d -dependent measurements in the human brain.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 1.

Tuned bipolar oscillating gradients for mapping frequency dispersion of diffusion kurtosis in the human brain

PURPOSE

Oscillating gradient spin-echo (OGSE) sequences have demonstrated an ability to probe time-dependent microstructural features, although they often suffer from low SNR due to increased TEs. In this work we introduce frequency-tuned bipolar (FTB) gradients as a variation of oscillating gradients with reduced TE and demonstrate their utility by mapping the frequency dispersion of kurtosis in human subjects.

METHODS

An FTB oscillating gradient waveform is presented that provides encoding of 1.5 net oscillation periods, thereby reducing the TE of the acquisition. Simulations were performed to determine an optimal protocol based on the SNR of kurtosis frequency dispersion-defined as the difference in kurtosis between pulsed and oscillating gradient acquisitions. Healthy human subjects were scanned at 7T using pulsed gradient and an optimized 23 Hz FTB protocol, which featured a maximum b-value of 2500 s/mm² . In addition, to directly compare existing methods, measurements using traditional cosine OGSE were also acquired.

RESULTS

FTB oscillating gradients demonstrated equivalent frequency-dependent diffusion measurements compared with cosine-modulated OGSE while enabling a significant reduction in TE. Optimization and in vivo results suggest that FTB gradients provide increased SNR of kurtosis dispersion maps compared with traditional cosine OGSE. The optimized FTB gradient protocol demonstrated consistent reductions in apparent kurtosis values and increased diffusivity in generated frequency dispersion maps.

CONCLUSIONS

This work presents an alternative to traditional cosine OGSE sequences, enabling more time-efficient acquisitions of frequency-dependent diffusion quantities as demonstrated through in vivo kurtosis frequency dispersion maps.

Metabolic activity diffusion imaging (MADI): II. Noninvasive, high-resolution human brain mapping of sodium pump flux and cell metrics

We introduce a new ¹ H₂ O magnetic resonance approach: metabolic activity diffusion imaging (MADI). Numerical diffusion-weighted imaging decay simulations characterized by the mean cellular water efflux (unidirectional) rate constant (k_io ), mean cell volume (V), and cell number density (ρ) are produced from Monte Carlo random walks in virtual stochastically sized/shaped cell ensembles. Because of active steady-state trans-membrane water cycling (AWC), k_io reflects the cytolemmal Na⁺ , K⁺ ATPase (NKA) homeostatic cellular metabolic rate (^c MR_NKA ). A digital 3D "library" contains thousands of simulated single diffusion-encoded (SDE) decays. Library entries match well with disparate, animal, and human experimental SDE decays. The V and ρ values are consistent with estimates from pertinent in vitro cytometric and ex vivo histopathological literature: in vivo V and ρ values were previously unavailable. The library allows noniterative pixel-by-pixel experimental SDE decay library matchings that can be used to advantage. They yield proof-of-concept MADI parametric mappings of the awake, resting human brain. These reflect the tissue morphology seen in conventional MRI. While V is larger in gray matter (GM) than in white matter (WM), the reverse is true for ρ. Many brain structures have k_io values too large for current, invasive methods. For example, the median WM k_io is 22s^-1 ; likely reflecting mostly exchange within myelin. The k_io •V product map displays brain tissue ^c MR_NKA variation. The GM activity correlates, quantitatively and qualitatively, with the analogous resting-state brain ¹⁸ FDG-PET tissue glucose consumption rate (^t MR_glucose ) map; but noninvasively, with higher spatial resolution, and no pharmacokinetic requirement. The cortex, thalamus, putamen, and caudate exhibit elevated metabolic activity. MADI accuracy and precision are assessed. The results are contextualized with literature overall homeostatic brain glucose consumption and ATP production/consumption measures. The MADI/PET results suggest different GM and WM metabolic pathways. Preliminary human prostate results are also presented.

Towards differentiation of brain tumor from radiation necrosis using multi-parametric MRI: Preliminary results at 4.7 T using rodent models

BACKGROUND

It remains a clinical challenge to differentiate brain tumors from radiation-induced necrosis in the brain. Despite significant improvements, no single MRI method has been validated adequately in the clinical setting.

METHODS

Multi-parametric MRI (mpMRI) was performed to differentiate 9L gliosarcoma from radiation necrosis in animal models. Five types of MRI methods probed complementary information on different scales i.e., T2 (relaxation), CEST based APT (probing mobile proteins/peptides) and rNOE (mobile macromolecules), qMT (macromolecules), diffusion based ADC (cell density) and SSIFT iAUC (cell size), and perfusion based DSC (blood volume and flow).

RESULTS

For single MRI parameters, iAUC and ADC provide the best discrimination of radiation necrosis and brain tumor. For mpMRI, a combination of iAUC, ADC, and APT shows the best classification performance based on a two-step analysis with the Lasso and Ridge regressions.

CONCLUSION

A general mpMRI approach is introduced to choosing candidate multiple MRI methods, identifying the most effective parameters from all the mpMRI parameters, and finding the appropriate combination of chosen parameters to maximize the classification performance to differentiate tumors from radiation necrosis.

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 Breast MRI: Current Standard and Emerging Techniques

The role of diffusion weighted imaging (DWI) as a biomarker has been the subject of active investigation in the field of breast radiology. By quantifying the random motion of water within a voxel of tissue, DWI provides indirect metrics that reveal cellularity and architectural features. Studies show that data obtained from DWI may provide information related to the characterization, prognosis, and treatment response of breast cancer. The incorporation of DWI in breast imaging demonstrates its potential to serve as a non-invasive tool to help guide diagnosis and treatment. In this review, current technical literature of diffusion-weighted breast imaging will be discussed, in addition to clinical applications, advanced techniques, and emerging use in the field of radiomics.

Impact of tissue properties on time-dependent alterations in apparent diffusion coefficient: a phantom study using oscillating-gradient spin-echo and pulsed-gradient spin-echo sequences

Purpose

The purpose of this study was to investigate whether the changes in apparent diffusion coefficients (ADCs) due to differences in diffusion time reflect tissue properties in actual measurements of phantoms.

Materials and methods

Various n-alkane phantoms and sucrose/collagen phantoms with various collagen densities were set up with and without polyvinyl alcohol (PVA) foam with an average pore diameter of 300 μm. Thus, n-alkanes or sucrose/collagen represented substrate viscosity and the presence of PVA foam represented tissue structure with septum. Diffusion-weighted images with various diffusion times (7.71–60 ms) were acquired using pulsed-gradient spin-echo (PGSE) and oscillating-gradient spin-echo (OGSE) sequences. The ADCs of the phantoms with and without PVA foam were calculated.

Results

The ADCs of some of the phantoms without PVA decreased with diffusion times decreased. In the n-alkane phantoms, only C₈H₁₈ showed significantly different ADCs depending on the use of PVA foam in the OGSE sequence. On the other hand, sucrose/collagen phantoms showed significant differences according to diffusion time. The ADCs of the phantoms decreased as the molecular size of the n-alkanes or collagen density of the sucrose/collagen phantom increased. Compared to phantoms without PVA foam, the ADC of the phantoms with PVA foam decreased as the diffusion time increased.

Conclusion

Changes in ADCs due to differences in diffusion time reflect tissue properties in actual measurements of phantoms. These changes in ADCs can be used for tissue characterization in vivo.

Investigation of breast cancer microstructure and microvasculature from time-dependent DWI and CEST in correlation with histological biomarkers

We investigated the associations of time-dependent DWI, non-Gaussian DWI, and CEST parameters with histological biomarkers in a breast cancer xenograft model. 22 xenograft mice (7 MCF-7 and 15 MDA-MB-231) were scanned at 4 diffusion times [T_d = 2.5/5 ms with 11 b-values (0–600 s/mm²) and T_d = 9/27.6 ms with 17 b-values (0–3000 s/mm²), respectively]. The apparent diffusion coefficient (ADC) was estimated using 2 b-values in different combinations (ADC^0–600 using b = 0 and 600 s/mm² and shifted ADC [sADC^200–1500] using b = 200 and 1500 s/mm²) at each of those diffusion times. Then the change (Δ) in ADC/sADC between diffusion times was evaluated. Non-Gaussian diffusion and intravoxel incoherent motion (IVIM) parameters (ADC₀, the virtual ADC at b = 0; K, Kurtosis from non-Gaussian diffusion; f, the IVIM perfusion fraction) were estimated. CEST images were acquired and the amide proton transfer signal intensity (APT SI) were measured. The ΔsADC_9–27.6 (between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{sADC}}_{{9\,{\text{ms}}}}^{200{-}1500}$$\end{document}sADC9ms200-1500 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{sADC}}_{{27.6\,{\text{ms}}}}^{200{-}1500}$$\end{document}sADC27.6ms200-1500 and ΔADC_2.5_sADC_27.6 (between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{ADC}}_{{2.5\, {\text{ms}}}}^{0{-}600}$$\end{document}ADC2.5ms0-600 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{sADC}}_{{27.6\,{\text{ms}}}}^{200{-}1500}$$\end{document}sADC27.6ms200-1500) was significantly larger for MCF-7 groups, and ΔADC_2.5_sADC_27.6 was positively correlated with Ki67_max and APT SI. ADC₀ decreased significantly in MDA-MB-231 group and K increased significantly with T_d in MCF-7 group. APT SI and cellular area had a moderately strong positive correlation in MDA-MB-231 and MCF-7 tumors combined, and there was a positive correlation in MDA-MB-231 tumors. There was a significant negative correlation between APT SI and the Ki-67-positive ratio in MDA-MB-231 tumors and when combined with MCF-7 tumors. The associations of ΔADC_2.5_sADC_27.6 and API SI with Ki-67 parameters indicate that the T_d-dependent DW and CEST parameters are useful to predict the histological markers of breast cancers.

Time-Dependent Diffusion MRI for Quantitative Microstructural Mapping of Prostate Cancer

Background Recently developed time-dependent diffusion MRI has potential in characterizing cellular tissue microstructures; however, its value in imaging prostate cancer (PCa) remains unknown. Purpose To investigate the feasibility of time-dependent diffusion MRI-based microstructural mapping for noninvasively characterizing cellular properties of PCa and for discriminating between clinically significant PCa and clinically insignificant disease. Materials and Methods Men with a clinical suspicion of PCa were enrolled prospectively between October 2019 and August 2020. Time-dependent diffusion MRI data were acquired with pulsed and oscillating gradient diffusion MRI sequences at an equivalent diffusion time of 7.5-30 msec on a 3.0-T scanner. Time-dependent diffusion MRI-based microstructural parameters, including cell diameter, intracellular volume fraction, cellularity, and diffusivities, were estimated with a two-compartment model. These were compared for different International Society of Urological Pathology grade groups (GGs), and their performance in discriminating clinically significant PCa (GG >1) from clinically insignificant disease (benign and GG 1) was determined with a linear discriminant analysis. The fitted microstructural parameters were validated by means of correlation with histopathologic measurements. Results In the 48 enrolled men, the time-dependent diffusion MRI measurements showed that higher GG was correlated with higher intracellular volume fraction and higher cellularity (intracellular volume fraction = 0.22, 0.36, 0.34, 0.37, and 0.40 in GGs 1-5, respectively; P < .001 at one-way analysis of variance), while lower cell diameter was found at higher GGs (diameter = 23.4, 18.3, 19.2, 17.9, and 18.5 μm in GGs 1-5, respectively; P = .002). Among all measurements derived from time-dependent diffusion MRI, cellularity achieved the highest diagnostic performance, with an accuracy of 92% (44 of 48 participants) and area under the receiver operating characteristic curve of 0.96 (95% CI: 0.87, 0.99) in discriminating clinically significant PCa from clinically insignificant disease. Microstructural mapping was supported by positive correlations between time-dependent diffusion MRI-based and pathologic examination-based intracellular volume fraction (r = 0.83; P < .001). Conclusion Time-dependent diffusion MRI-based microstructural mapping correlates with pathologic findings and demonstrates promise for characterizing prostate cancer. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Chatterjee and Oto in this issue.

Clinical experience of tensor-valued diffusion encoding for microstructure imaging by diffusional variance decomposition in patients with breast cancer

BACKGROUND

Diffusion-weighted imaging plays a key role in magnetic resonance imaging (MRI) of breast tumors. However, it remains unclear how to interpret single diffusion encoding with respect to its link with tissue microstructure. The purpose of this retrospective cross-sectional study was to use tensor-valued diffusion encoding to investigate the underlying microstructure of invasive ductal carcinoma (IDC) and evaluate its potential value in a clinical setting.

METHODS

We retrospectively reviewed biopsy-proven breast cancer patients who underwent preoperative breast MRI examination from July 2020 to March 2021. We reviewed the MRI of 29 patients with 30 IDCs, including analysis by diffusional variance decomposition enabled by tensor-valued diffusion encoding. The diffusion parameters of mean diffusivity (MD), total mean kurtosis (MKT), anisotropic mean kurtosis (MKA), isotropic mean kurtosis (MKI), macroscopic fractional anisotropy (FA), and microscopic fractional anisotropy (µFA) were estimated. The parameter differences were compared between IDC and normal fibroglandular breast tissue (FGBT), as well as the association between the diffusion parameters and histopathologic items.

RESULTS

The mean value of MD in IDCs was significantly lower than that of normal FGBT (1.07±0.27 vs. 1.34±0.29, P<0.001); however, MKT, MKA, MKI, FA, and µFA were significantly higher (P<0.005). Among all the diffusion parameters, MKI was positively correlated with the tumor size on both MRI and pathological specimen (rs=0.38, P<0.05 vs. rs=0.54, P<0.01), whereas MKT had a positive correlation with the tumor size in the pathological specimen only (rs=0.47, P<0.02). In addition, the lymph node (LN) metastasis group had significantly higher MKT, MKA, and µFA compared to the metastasis negative group (P<0.05).

CONCLUSIONS

Tensor-valued diffusion encoding enables a useful non-invasive method for characterizing breast cancers with information on tissue microstructures. Particularly, µFA could be a potential imaging biomarker for evaluating breast cancers prior to surgery or chemotherapy.

Advanced Diffusion MR Imaging for Multiple Sclerosis in the Brain and Spinal Cord

Diffusion tensor imaging (DTI) has been established its usefulness in evaluating normal-appearing white matter (NAWM) and other lesions that are difficult to evaluate with routine clinical MRI in the evaluation of the brain and spinal cord lesions in multiple sclerosis (MS), a demyelinating disease. With the recent advances in the software and hardware of MRI systems, increasingly complex and sophisticated MRI and analysis methods, such as q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, white matter tract integrity, and multiple diffusion encoding, referred to as advanced diffusion MRI, have been proposed. These are capable of capturing in vivo microstructural changes in the brain and spinal cord in normal and pathological states in greater detail than DTI.

This paper reviews the current status of recent advanced diffusion MRI for assessing MS in vivo as part of an issue celebrating two decades of magnetic resonance in medical sciences (MRMS), an official journal of the Japanese Society of Magnetic Resonance in Medicine.

Quantitative magnetic resonance imaging and tumor forecasting of breast cancer patients in the community setting

This protocol describes a complete data acquisition, analysis and computational forecasting pipeline for employing quantitative MRI data to predict the response of locally advanced breast cancer to neoadjuvant therapy in a community-based care setting. The methodology has previously been successfully applied to a heterogeneous patient population. The protocol details how to acquire the necessary images followed by registration, segmentation, quantitative perfusion and diffusion analysis, model calibration, and prediction. The data collection portion of the protocol requires ~25 min of scanning, postprocessing requires 2-3 h, and the model calibration and prediction components require ~10 h per patient depending on tumor size. The response of individual breast cancer patients to neoadjuvant therapy is forecast by application of a biophysical, reaction-diffusion mathematical model to these data. Successful application of the protocol results in coregistered MRI data from at least two scan visits that quantifies an individual tumor's size, cellularity and vascular properties. This enables a spatially resolved prediction of how a particular patient's tumor will respond to therapy. Expertise in image acquisition and analysis, as well as the numerical solution of partial differential equations, is required to carry out this protocol.

Non-invasive cell-tracking methods for adoptive T cell therapies

Adoptive T cell therapies (ACT) have demonstrated groundbreaking results in blood cancers and melanoma. Nevertheless, their significant cost, the occurrence of severe adverse events, and their poor performance in solid tumors are important hurdles hampering more widespread applicability. In vivo cell tracking allows instantaneous and non-invasive monitoring of the distribution, tumor homing, persistence, and redistribution to other organs of infused T cells in patients. Furthermore, cell tracking could aid in the clinical management of patients, allowing the detection of non-responders or severe adverse events at an early stage. This review provides a concise overview of the main principles and potential of cell tracking, followed by a discussion of the clinically relevant labeling strategies and their application in ACT.

Feasibility of oscillating and pulsed gradient diffusion MRI to assess neonatal hypoxia-ischemia on clinical systems

Diffusion-time- (t_d) dependent diffusion MRI (dMRI) extends our ability to characterize brain microstructure by measuring dMRI signals at varying t_d. The use of oscillating gradient (OG) is essential for accessing short t_d but is technically challenging on clinical MRI systems. This study aims to investigate the clinical feasibility and value of t_d-dependent dMRI in neonatal hypoxic-ischemic encephalopathy (HIE). Eighteen HIE neonates and six normal term-born neonates were scanned on a 3 T scanner, with OG-dMRI at an oscillating frequency of 33 Hz (equivalent t_d ≈ 7.5 ms) and pulsed gradient (PG)-dMRI at a t_d of 82.8 ms and b-value of 700 s/mm². The t_d-dependence, as quantified by the difference in apparent diffusivity coefficients between OG- and PG-dMRI (ΔADC), was observed in the normal neonatal brains, and the ΔADC was higher in the subcortical white matter than the deep grey matter. In HIE neonates with severe and moderate injury, ΔADC significantly increased in the basal ganglia (BG) compared to the controls (23.7% and 10.6%, respectively). In contrast, the conventional PG-ADC showed a 12.6% reduction only in the severe HIE group. White matter edema regions also demonstrated increased ΔADC, where PG-ADC did not show apparent changes. Our result demonstrated that t_d-dependent dMRI provided high sensitivity in detecting moderate-to-severe HIE.

Time-dependent diffusion MRI probes cerebellar microstructural alterations in a mouse model of Down syndrome

The cerebellum is a complex system with distinct cortical laminar organization. Alterations in cerebellar microstructure are common and associated with many factors such as genetics, cancer and ageing. Diffusion MRI (dMRI) provides a non-invasive tool to map the brain structural organization, and the recently proposed diffusion-time (td )-dependent dMRI further improves its capability to probe the cellular and axonal/dendritic microstructures by measuring water diffusion at multiple spatial scales. The td -dependent diffusion profile in the cerebellum and its utility in detecting cerebellar disorders, however, are not yet elucidated. Here, we first deciphered the spatial correspondence between dMRI contrast and cerebellar layers, based on which the cerebellar layer-specific td -dependent dMRI patterns were characterized in both euploid and Ts65Dn mice, a mouse model of Down syndrome. Using oscillating gradient dMRI, which accesses diffusion at short td 's by modulating the oscillating frequency, we detected subtle changes in the apparent diffusivity coefficient of the cerebellar internal granular layer and Purkinje cell layer of Ts65Dn mice that were not detectable by conventional pulsed gradient dMRI. The detection sensitivity of oscillating gradient dMRI increased with the oscillating frequency at both the neonatal and adult stages. The td -dependence, quantified by ΔADC map, was reduced in Ts65Dn mice, likely associated with the reduced granule cell density and abnormal dendritic arborization of Purkinje cells as revealed from histological evidence. Our study demonstrates superior sensitivity of short-td diffusion using oscillating gradient dMRI to detect cerebellar microstructural changes in Down syndrome, suggesting the potential application of this technique in cerebellar disorders.

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.

Multidimensional Diffusion Magnetic Resonance Imaging for Characterization of Tissue Microstructure in Breast Cancer Patients: A Prospective Pilot Study

Diffusion-weighted imaging is a non-invasive functional imaging modality for breast tumor characterization through apparent diffusion coefficients. Yet, it has so far been unable to intuitively inform on tissue microstructure. In this IRB-approved prospective study, we applied novel multidimensional diffusion (MDD) encoding across 16 patients with suspected breast cancer to evaluate its potential for tissue characterization in the clinical setting. Data acquired via custom MDD sequences was processed using an algorithm estimating non-parametric diffusion tensor distributions. The statistical descriptors of these distributions allow us to quantify tissue composition in terms of metrics informing on cell densities, shapes, and orientations. Additionally, signal fractions from specific cell types, such as elongated cells (bin1), isotropic cells (bin2), and free water (bin3), were teased apart. Histogram analysis in cancers and healthy breast tissue showed that cancers exhibited lower mean values of "size" (1.43 ± 0.54 × 10-3 mm2/s) and higher mean values of "shape" (0.47 ± 0.15) corresponding to bin1, while FGT (fibroglandular breast tissue) presented higher mean values of "size" (2.33 ± 0.22 × 10-3 mm2/s) and lower mean values of "shape" (0.27 ± 0.11) corresponding to bin3 (p < 0.001). Invasive carcinomas showed significant differences in mean signal fractions from bin1 (0.64 ± 0.13 vs. 0.4 ± 0.25) and bin3 (0.18 ± 0.08 vs. 0.42 ± 0.21) compared to ductal carcinomas in situ (DCIS) and invasive carcinomas with associated DCIS (p = 0.03). MDD enabled qualitative and quantitative evaluation of the composition of breast cancers and healthy glands.

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

PURPOSE

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

METHOD

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

RESULTS

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

CONCLUSIONS

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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