Integrated molecular and multiparametric MRI mapping of high-grade glioma identifies regional biologic signatures.

Sampling restrictions have hindered the comprehensive study of invasive non-enhancing (NE) high-grade glioma (HGG) cell populations driving tumor progression. Here, we present an integrated multi-omic analysis of spatially matched molecular and multi-parametric magnetic resonance imaging (MRI) profiling across 313 multi-regional tumor biopsies, including 111 from the NE, across 68 HGG patients. Whole exome and RNA sequencing uncover unique genomic alterations to unresectable invasive NE tumor, including subclonal events, which inform genomic models predictive of geographic evolution. Infiltrative NE tumor is alternatively enriched with tumor cells exhibiting neuronal or glycolytic/plurimetabolic cellular states, two principal transcriptomic pathway-based glioma subtypes, which respectively demonstrate abundant private mutations or enrichment in immune cell signatures. These NE phenotypes are non-invasively identified through normalized K2 imaging signatures, which discern cell size heterogeneity on dynamic susceptibility contrast (DSC)-MRI. NE tumor populations predicted to display increased cellular proliferation by mean diffusivity (MD) MRI metrics are uniquely associated with EGFR amplification and CDKN2A homozygous deletion. The biophysical mapping of infiltrative HGG potentially enables the clinical recognition of tumor subpopulations with aggressive molecular signatures driving tumor progression, thereby informing precision medicine targeting.

Diffusion time dependency of extracellular diffusion.

PURPOSE: To quantify the variations of the power-law dependences on diffusion time t or gradient frequency f $$ f $$ of extracellular water diffusion measured by diffusion MRI (dMRI). METHODS: Model cellular systems containing only extracellular water were used to investigate the t / f $$ t/f $$ dependence of D ex $$ {D}_{ex} $$ , the extracellular diffusion coefficient. Computer simulations used a randomly packed tissue model with realistic intracellular volume fractions and cell sizes. DMRI measurements were performed on samples consisting of liposomes containing heavy water(D2 O, deuterium oxide) dispersed in regular water (H2 O). D ex $$ {D}_{ex} $$ was obtained over a broad t $$ t $$ range (∼1-1000 ms) and then fit power-law equations D ex ( t ) = D const + const · t - ϑ t $$ {D}_{ex}(t)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {t}^{-{\vartheta}_t} $$ and D ex ( f ) = D const + const · f ϑ f $$ {D}_{ex}(f)={D}_{\mathrm{const}}+\mathrm{const}\cdotp {f}^{\vartheta_f} $$ . RESULTS: Both simulated and experimental results suggest that no single power-law adequately describes the behavior of D ex $$ {D}_{ex} $$ over the range of diffusion times of most interest in practical dMRI. Previous theoretical predictions are accurate over only limited t $$ t $$ ranges; for example, θ t = θ f = - 1 2 $$ {\theta}_t={\theta}_f=-\frac{1}{2} $$ is valid only for short times, whereas θ t = 1 $$ {\theta}_t=1 $$ or θ f = 3 2 $$ {\theta}_f=\frac{3}{2} $$ is valid only for long times but cannot describe other ranges simultaneously. For the specific t $$ t $$ range of 5-70 ms used in typical human dMRI measurements, θ t = θ f = 1 $$ {\theta}_t={\theta}_f=1 $$ matches the data well empirically. CONCLUSION: The optimal power-law fit of extracellular diffusion varies with diffusion time. The dependency obtained at short or long t $$ t $$ limits cannot be applied to typical dMRI measurements in human cancer or liver. It is essential to determine the appropriate diffusion time range when modeling extracellular diffusion in dMRI-based quantitative microstructural imaging.

Longitudinal evaluation of myofiber microstructural changes in a preclinical ALS model using the transverse relaxivity at tracer equilibrium (TRATE): A preliminary study.

T2⁎ relaxivity contrast imaging may serve as a potential imaging biomarker for amyotrophic lateral sclerosis (ALS) by noninvasively quantifying the tissue microstructure. In this preliminary longitudinal study, we investigated the Transverse Relaxivity at Tracer Equilibrium (TRATE) in three muscle groups between SOD1-G93A (ALS model) rat and a control population at two different timepoints. The control group was time matched to the ALS group such that the second timepoint was the onset of disease. We observed a statistically significant decrease in TRATE over time in the gastrocnemius, tibialis, and digital flexor muscles in the SOD1-G93A model (p-value = 0.003, 0.008, 0.005; respectively), whereas TRATE did not change over time in the control group (p-value = 0.4777, 0.6837, 0.9682; respectively). Immunofluorescent staining revealed a decrease in minimum fiber area and cell density in the SOD1-G93A model when compared to the control group (p-value = 6.043E-10 and 2.265E-10, respectively). These microstructural changes observed from histology align with the theorized biophysical properties of TRATE. We demonstrate that TRATE can longitudinally differentiate disease associated atrophy from healthy muscle and has potential to serve as a biomarker for disease progression and ultimately therapy response in patients with ALS.

Systematic assessment of multi-echo dynamic susceptibility contrast MRI using a digital reference object.

PURPOSE: Brain tumor dynamic susceptibility contrast (DSC) MRI is adversely impacted by T1 and T2∗ contrast agent leakage effects that result in inaccurate hemodynamic metrics. While multi-echo acquisitions remove T1 leakage effects, there is no consensus on the optimal set of acquisition parameters. Using a computational approach, we systematically evaluated a wide range of acquisition strategies to determine the optimal multi-echo DSC-MRI perfusion protocol. METHODS: Using a population-based DSC-MRI digital reference object (DRO), we assessed the influence of preload dosing (no preload and full dose preload), field strength (1.5 and 3T), pulse sequence parameters (echo time, repetition time, and flip angle), and leakage correction on relative cerebral blood volume (rCBV) and flow (rCBF) accuracy. We also compared multi-echo DSC-MRI protocols with standard single-echo protocols. RESULTS: Multi-echo DSC-MRI is highly consistent across all protocols, and multi-echo rCBV (with or without use of a preload dose) had higher accuracy than single-echo rCBV. Regression analysis showed that choice of repetition time and flip angle had minimal impact on multi-echo rCBV and rCBV, indicating the potential for significant flexibility in acquisition parameters. The echo time combination had minimal impact on rCBV, though longer echo times should be avoided, particularly at higher field strengths. Leakage correction improved rCBV accuracy in all cases. Multi-echo rCBF was less biased than single-echo rCBF, although rCBF accuracy was reduced overall relative to rCBV. CONCLUSIONS: Multi-echo acquisitions were more robust than single-echo, essentially decoupling both repetition time and flip angle from rCBV accuracy. Multi-echo acquisitions obviate the need for preload dosing, although leakage correction to remove residual T2∗ leakage effects remains compulsory for high rCBV accuracy.

A Population-Based Digital Reference Object (DRO) for Optimizing Dynamic Susceptibility Contrast (DSC)-MRI Methods for Clinical Trials.

The standardization and broad-scale integration of dynamic susceptibility contrast (DSC)-magnetic resonance imaging (MRI) have been confounded by a lack of consensus on DSC-MRI methodology for preventing potential relative cerebral blood volume inaccuracies, including the choice of acquisition protocols and postprocessing algorithms. Therefore, we developed a digital reference object (DRO), using physiological and kinetic parameters derived from in vivo data, unique voxel-wise 3-dimensional tissue structures, and a validated MRI signal computational approach, aimed at validating image acquisition and analysis methods for accurately measuring relative cerebral blood volume in glioblastomas. To achieve DSC-MRI signals representative of the temporal characteristics, magnitude, and distribution of contrast agent-induced T1 and changes observed across multiple glioblastomas, the DRO's input parameters were trained using DSC-MRI data from 23 glioblastomas (>40 000 voxels). The DRO's ability to produce reliable signals for combinations of pulse sequence parameters and contrast agent dosing schemes unlike those in the training data set was validated by comparison with in vivo dual-echo DSC-MRI data acquired in a separate cohort of patients with glioblastomas. Representative applications of the DRO are presented, including the selection of DSC-MRI acquisition and postprocessing methods that optimize CBV accuracy, determination of the impact of DSC-MRI methodology choices on sample size requirements, and the assessment of treatment response in clinical glioblastoma trials.

Assessing tumor cytoarchitecture using multiecho DSC-MRI derived measures of the transverse relaxivity at tracer equilibrium (TRATE).

PURPOSE: In brain tumor dynamic susceptibility contrast (DSC)-MRI studies, multiecho acquisition methods are used to quantify the dynamic changes in T1 and T2 * that occur when contrast agent (CA) extravasates. Such methods also enable the estimation of the effective tissue CA transverse relaxivity. The goal of this study was to evaluate the sensitivity of the transverse relaxivity at tracer equilibrium (TRATE) to tumor cytoarchitecture. METHODS: Computational and in vitro studies were used to evaluate the biophysical basis of TRATE. In 9L, C6, and human brain tumors, TRATE, the apparent diffusion coefficient (ADC), the CA transfer constant (K(trans) ), the extravascular extracellular volume fraction (ve ), and histological data were compared. RESULTS: Simulations and in vitro results indicate that TRATE is highly sensitive to variations in cellular properties such as cell size and density. The histologic cell density and TRATE values were significantly higher in 9L tumors as compared to C6 tumors. In animal and human tumors, a voxel-wise comparison of TRATE with ADC, ve , and K(trans) maps showed low spatial correlation. CONCLUSION: The assessment of TRATE is clinically feasible and its sensitivity to tissue cytoarchitectural features not present in other imaging methods indicate that it could potentially serve as a unique structural signature or "trait" of cancer.

An efficient computational approach to characterize DSC-MRI signals arising from three-dimensional heterogeneous tissue structures.

The systematic investigation of susceptibility-induced contrast in MRI is important to better interpret the influence of microvascular and microcellular morphology on DSC-MRI derived perfusion data. Recently, a novel computational approach called the Finite Perturber Method (FPM), which enables the study of susceptibility-induced contrast in MRI arising from arbitrary microvascular morphologies in 3D has been developed. However, the FPM has lower efficiency in simulating water diffusion especially for complex tissues. In this work, an improved computational approach that combines the FPM with a matrix-based finite difference method (FDM), which we call the Finite Perturber the Finite Difference Method (FPFDM), has been developed in order to efficiently investigate the influence of vascular and extravascular morphological features on susceptibility-induced transverse relaxation. The current work provides a framework for better interpreting how DSC-MRI data depend on various phenomena, including contrast agent leakage in cancerous tissues and water diffusion rates. In addition, we illustrate using simulated and micro-CT extracted tissue structures the improved FPFDM along with its potential applications and limitations.