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

Simultaneous quantification of perfusion, permeability, and leakage effects in brain gliomas using dynamic spin-and-gradient-echo echoplanar imaging MRI

OBJECTIVE

To determine the feasibility and biologic correlations of dynamic susceptibility contrast (DSC), dynamic contrast enhanced (DCE), and quantitative maps derived from contrast leakage effects obtained simultaneously in gliomas using dynamic spin-and-gradient-echo echoplanar imaging (dynamic SAGE-EPI) during a single contrast injection.

MATERIALS AND METHODS

Thirty-eight patients with enhancing brain gliomas were prospectively imaged with dynamic SAGE-EPI, which was processed to compute traditional DSC metrics (normalized relative cerebral blood flow [nrCBV], percentage of signal recovery [PSR]), DCE metrics (volume transfer constant [Ktrans], extravascular compartment [ve]), and leakage effect metrics: ΔR2,ss* (reflecting T2*-leakage effects), ΔR1,ss (reflecting T1-leakage effects), and the transverse relaxivity at tracer equilibrium (TRATE, reflecting the balance between ΔR2,ss* and ΔR1,ss). These metrics were compared between patient subgroups (treatment-naïve [TN] vs recurrent [R]) and biological features (IDH status, Ki67 expression).

RESULTS

In IDH wild-type gliomas (IDHwt-i.e., glioblastomas), previous exposure to treatment determined lower TRATE (p = 0.002), as well as higher PSR (p = 0.006), Ktrans (p = 0.17), ΔR1,ss (p = 0.035), ve (p = 0.006), and ADC (p = 0.016). In IDH-mutant gliomas (IDHm), previous treatment determined higher Ktrans and ΔR1,ss (p = 0.026). In TN-gliomas, dynamic SAGE-EPI metrics tended to be influenced by IDH status (p ranging 0.09-0.14). TRATE values above 142 mM-1s-1 were exclusively seen in TN-IDHwt, and, in TN-gliomas, this cutoff had 89% sensitivity and 80% specificity as a predictor of Ki67 > 10%.

CONCLUSIONS

Dynamic SAGE-EPI enables simultaneous quantification of brain tumor perfusion and permeability, as well as mapping of novel metrics related to cytoarchitecture (TRATE) and blood-brain barrier disruption (ΔR1,ss), with a single contrast injection.

CLINICAL RELEVANCE STATEMENT

Simultaneous DSC and DCE analysis with dynamic SAGE-EPI reduces scanning time and contrast dose, respectively alleviating concerns about imaging protocol length and gadolinium adverse effects and accumulation, while providing novel leakage effect metrics reflecting blood-brain barrier disruption and tumor tissue cytoarchitecture.

KEY POINTS

• Traditionally, perfusion and permeability imaging for brain tumors requires two separate contrast injections and acquisitions. • Dynamic spin-and-gradient-echo echoplanar imaging enables simultaneous perfusion and permeability imaging. • Dynamic spin-and-gradient-echo echoplanar imaging provides new image contrasts reflecting blood-brain barrier disruption and cytoarchitecture characteristics.

Multiparametric analysis from dynamic susceptibility contrast-enhanced perfusion MRI to evaluate malignant brain tumors

BACKGROUND AND PURPOSE

Dynamic susceptibility contrast-enhanced (DSC) MR perfusion is a valuable technique for distinguishing brain tumors. Diagnostic potential of measurable parameters derived from preload leakage-corrected-DSC-MRI remains somewhat underexplored. This study aimed to evaluate these parameters for differentiating primary CNS lymphoma (PCNSL), glioblastoma, and metastasis.

METHODS

Thirty-nine patients with pathologically proven PCNSL (n = 14), glioblastoma (n = 14), and metastasis (n = 11) were analyzed. Five DSC parameters-relative CBV (rCBV), percentage of signal recovery (PSR), downward slope (DS), upward slope (US), and first-pass slope ratio-were derived from tumor-enhancing areas. Diagnostic performance was assessed using receiver operating characteristic curve analysis.

RESULTS

RCBV was higher in metastasis (4.58; interquartile range [IQR]: 2.54) and glioblastoma (3.98; IQR: 1.87), compared with PCNSL (1.46; IQR: 0.29; p = .00006 for both). rCBV better distinguished metastasis and glioblastoma from PCNSL, with an area under the curve (AUC) of 0.97 and 0.99, respectively. PSR was higher in PCNSL (88.11; IQR: 21.21) than metastases (58.30; IQR: 22.28; p = .0002), while glioblastoma (74.54; IQR: 21.23) presented almost significant trend-level differences compared to the others (p≈.05). AUCs were 0.79 (PCNSL vs. glioblastoma), 0.91 (PCNSL vs. metastasis), and 0.78 (glioblastoma vs. metastasis). DS and US parameters were statistically significant between glioblastoma (-109.92; IQR: 152.71 and 59.06; IQR: 52.87) and PCNSL (-47.36; IQR: 44.30 and 21.68; IQR: 16.85), presenting AUCs of 0.86 and 0.87.

CONCLUSION

Metastasis and glioblastoma can be better differentiated from PCNSL through rCBV. PSR demonstrated higher differential performance compared to the other parameters and seemed useful, allowing a proper distinction among all, particularly between metastasis and glioblastoma, where rCBV failed. Finally, DS and US were only helpful in differentiating glioblastoma from PCNSL.

Image-localized biopsy mapping of brain tumor heterogeneity: A single-center study protocol

Brain cancers pose a novel set of difficulties due to the limited accessibility of human brain tumor tissue. For this reason, clinical decision-making relies heavily on MR imaging interpretation, yet the mapping between MRI features and underlying biology remains ambiguous. Standard (clinical) tissue sampling fails to capture the full heterogeneity of the disease. Biopsies are required to obtain a pathological diagnosis and are predominantly taken from the tumor core, which often has different traits to the surrounding invasive tumor that typically leads to recurrent disease. One approach to solving this issue is to characterize the spatial heterogeneity of molecular, genetic, and cellular features of glioma through the intraoperative collection of multiple image-localized biopsy samples paired with multi-parametric MRIs. We have adopted this approach and are currently actively enrolling patients for our 'Image-Based Mapping of Brain Tumors' study. Patients are eligible for this research study (IRB #16-002424) if they are 18 years or older and undergoing surgical intervention for a brain lesion. Once identified, candidate patients receive dynamic susceptibility contrast (DSC) perfusion MRI and diffusion tensor imaging (DTI), in addition to standard sequences (T1, T1Gd, T2, T2-FLAIR) at their presurgical scan. During surgery, sample anatomical locations are tracked using neuronavigation. The collected specimens from this research study are used to capture the intra-tumoral heterogeneity across brain tumors including quantification of genetic aberrations through whole-exome and RNA sequencing as well as other tissue analysis techniques. To date, these data (made available through a public portal) have been used to generate, test, and validate predictive regional maps of the spatial distribution of tumor cell density and/or treatment-related key genetic marker status to identify biopsy and/or treatment targets based on insight from the entire tumor makeup. This type of methodology, when delivered within clinically feasible time frames, has the potential to further inform medical decision-making by improving surgical intervention, radiation, and targeted drug therapy for patients with glioma.

Standardized brain tumor imaging protocols for clinical trials: current recommendations and tips for integration

Standardized MRI acquisition protocols are crucial for reducing the measurement and interpretation variability associated with response assessment in brain tumor clinical trials. The main challenge is that standardized protocols should ensure high image quality while maximizing the number of institutions meeting the acquisition requirements. In recent years, extensive effort has been made by consensus groups to propose different "ideal" and "minimum requirements" brain tumor imaging protocols (BTIPs) for gliomas, brain metastases (BM), and primary central nervous system lymphomas (PCSNL). In clinical practice, BTIPs for clinical trials can be easily integrated with additional MRI sequences that may be desired for clinical patient management at individual sites. In this review, we summarize the general concepts behind the choice and timing of sequences included in the current recommended BTIPs, we provide a comparative overview, and discuss tips and caveats to integrate additional clinical or research sequences while preserving the recommended BTIPs. Finally, we also reflect on potential future directions for brain tumor imaging in clinical trials.

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.

Progress of radiological‑pathological workflows in the differential diagnosis between primary central nervous system lymphoma and high‑grade glioma (Review)

Primary central nervous system lymphoma (PCNSL) and high‑grade glioma (HGG) are distinct entities of the CNS with completely distinct treatments. The treatment of PCNSL is chemotherapy‑based, while surgery is the first choice for HGG. However, the clinical features of the two entities often overlap, and a clear pathological diagnosis is important for subsequent management, especially for the management of PCNSL. Stereotactic biopsy is recognized as one of the minimally invasive alternatives for evaluating the involvement of the CNS. However, in the case of limited tissue materials, the differential diagnosis between the two entities is still difficult. In addition, some patients are too ill to tolerate a needle biopsy. Therefore, combining imaging, histopathology and laboratory examinations is essential in order to make a clear diagnosis as soon as possible. The present study reviews the progress of comparative research on both imaging and laboratory tests based on the pathophysiological changes of the two entities, and proposes an integrative and optimized diagnostic process, with the purpose of building a better understanding for neurologists, hematologists, radiologists and pathologists.

Evaluation of Amyotrophic Lateral Sclerosis-Induced Muscle Degeneration Using Magnetic Resonance-Based Relaxivity Contrast Imaging (RCI)

(1) Background: This work characterizes the sensitivity of magnetic resonance-based Relaxivity Contrast Imaging (RCI) to Amyotrophic Lateral Sclerosis (ALS)-induced changes in myofiber microstructure. Transverse Relaxivity at Tracer Equilibrium (TRATE), an RCI-based parameter, was evaluated in the lower extremities of ALS patients and healthy subjects. (2) Methods: In this IRB-approved study, 23 subjects (12 ALS patients and 11 healthy controls) were scanned at 3T (Philips, The Netherlands). RCI data were obtained during injection of a gadolinium-based contrast agent. TRATE, fat fraction and T₂ measures, were compared in five muscle groups of the calf muscle, between ALS and control populations. TRATE was also evaluated longitudinally (baseline and 6 months) and was compared to clinical measures, namely ALS Functional Rating Scale (ALSFRS-R) and Hand-Held Dynamometry (HHD), in a subset of the ALS population. (3) Results: TRATE was significantly lower (p < 0.001) in ALS-affected muscle than in healthy muscle in all muscle groups. Fat fraction differences between ALS and healthy muscle were statistically significant for the tibialis anterior (p = 0.01), tibialis posterior (p = 0.004), and peroneus longus (p = 0.02) muscle groups but were not statistically significant for the medial (p = 0.07) and lateral gastrocnemius (p = 0.06) muscles. T₂ differences between ALS and healthy muscle were statistically significant for the tibialis anterior (p = 0.004), peroneus longus (p = 0.004) and lateral gastrocnemius (p = 0.03) muscle groups but were not statistically significant for the tibialis posterior (p = 0.06) and medial gastrocnemius (p = 0.07) muscles. Longitudinally, TRATE, averaged over all patients, decreased by 28 ± 16% in the tibialis anterior, 47 ± 18% in the peroneus longus, 25 ± 19% in the tibialis posterior, 29 ± 14% in the medial gastrocnemius and 35 ± 18% in the lateral gastrocnemius muscles between two timepoints. ALSFRS-R scores were stable in two of four ALS patients. HHD scores decreased in three of four ALS patients. (4) Conclusion: RCI-based TRATE was shown to consistently differentiate ALS-affected muscle from healthy muscle and also provide a quantitative measure of longitudinal muscle degeneration.

Uncertainty quantification in the radiogenomics modeling of EGFR amplification in glioblastoma

Radiogenomics uses machine-learning (ML) to directly connect the morphologic and physiological appearance of tumors on clinical imaging with underlying genomic features. Despite extensive growth in the area of radiogenomics across many cancers, and its potential role in advancing clinical decision making, no published studies have directly addressed uncertainty in these model predictions. We developed a radiogenomics ML model to quantify uncertainty using transductive Gaussian Processes (GP) and a unique dataset of 95 image-localized biopsies with spatially matched MRI from 25 untreated Glioblastoma (GBM) patients. The model generated predictions for regional EGFR amplification status (a common and important target in GBM) to resolve the intratumoral genetic heterogeneity across each individual tumor-a key factor for future personalized therapeutic paradigms. The model used probability distributions for each sample prediction to quantify uncertainty, and used transductive learning to reduce the overall uncertainty. We compared predictive accuracy and uncertainty of the transductive learning GP model against a standard GP model using leave-one-patient-out cross validation. Additionally, we used a separate dataset containing 24 image-localized biopsies from 7 high-grade glioma patients to validate the model. Predictive uncertainty informed the likelihood of achieving an accurate sample prediction. When stratifying predictions based on uncertainty, we observed substantially higher performance in the group cohort (75% accuracy, n = 95) and amongst sample predictions with the lowest uncertainty (83% accuracy, n = 72) compared to predictions with higher uncertainty (48% accuracy, n = 23), due largely to data interpolation (rather than extrapolation). On the separate validation set, our model achieved 78% accuracy amongst the sample predictions with lowest uncertainty. We present a novel approach to quantify radiogenomics uncertainty to enhance model performance and clinical interpretability. This should help integrate more reliable radiogenomics models for improved medical decision-making.

Consensus recommendations for a dynamic susceptibility contrast MRI protocol for use in high-grade gliomas

Despite the widespread clinical use of dynamic susceptibility contrast (DSC) MRI, DSC-MRI methodology has not been standardized, hindering its utilization for response assessment in multicenter trials. Recently, the DSC-MRI Standardization Subcommittee of the Jumpstarting Brain Tumor Drug Development Coalition issued an updated consensus DSC-MRI protocol compatible with the standardized brain tumor imaging protocol (BTIP) for high-grade gliomas that is increasingly used in the clinical setting and is the default MRI protocol for the National Clinical Trials Network. After reviewing the basis for controversy over DSC-MRI protocols, this paper provides evidence-based best practices for clinical DSC-MRI as determined by the Committee, including pulse sequence (gradient echo vs spin echo), BTIP-compliant contrast agent dosing (preload and bolus), flip angle (FA), echo time (TE), and post-processing leakage correction. In summary, full-dose preload, full-dose bolus dosing using intermediate (60°) FA and field strength-dependent TE (40-50 ms at 1.5 T, 20-35 ms at 3 T) provides overall best accuracy and precision for cerebral blood volume estimates. When single-dose contrast agent usage is desired, no-preload, full-dose bolus dosing using low FA (30°) and field strength-dependent TE provides excellent performance, with reduced contrast agent usage and elimination of potential systematic errors introduced by variations in preload dose and incubation time.

Validation of combined use of DWI and percentage signal recovery-optimized protocol of DSC-MRI in differentiation of high-grade glioma, metastasis, and lymphoma

PURPOSE

With conventional MRI, it is often difficult to effectively differentiate between contrast-enhancing brain tumors, including primary central nervous system lymphoma (PCNSL), high-grade glioma (HGG), and metastasis. This study aimed to assess the discrimination ability of the parameters obtained from DWI and the percentage signal recovery- (PSR-) optimized protocol of DSC-MRI between these three tumor types at an initial step.

METHODS

DSC-MRI using a PSR-optimized protocol (TR/TE = 1500/30 ms, flip angle = 90°, no preload) and DWI of 99 solitary enhancing tumors (60 HGGs, 24 metastases, 15 PCNSLs) were retrospectively assessed before treatment. rCBV, PSR, ADC in the tumor core and rCBV, and ADC in peritumoral edema were measured. The differences were evaluated using one-way ANOVA, and the diagnostic performance was evaluated using ROC curve analysis.

RESULTS

PSR in the tumor core showed the best discriminating performance in differentiating these three tumor types with AUC values of 0.979 for PCNSL vs. others and 0.947 for HGG vs. metastasis. The ADC was only helpful in the tumor core and distinguishing PCNSLs from others (AUC = 0.897).

CONCLUSION

Different from CBV-optimized protocols (preload, intermediate FA), PSR derived from the PSR-optimized protocol seems to be the most important parameter in the differentiation of HGGs, metastases, and PCNSLs at initial diagnosis. This property makes PSR remarkable and carries the need for comprehensive DSC-MRI protocols, which provides PSR sensitivity and CBV accuracy together, such as the preload use of the PSR-optimized protocol before the CBV-optimized protocol.

R2* Relaxation Affects Pharmacokinetic Analysis of Dynamic Contrast-Enhanced MRI in Cancer and Underestimates Treatment Response at 7 T

Effective transverse relaxivity of gadolinium-based contrast agents is often neglected in dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). Here, we assess time and tissue dependence of R2* enhancement and its impact on pharmacokinetic parameter quantification and treatment monitoring. Multiecho DCE-MRI was performed at 7 T on mice bearing subcutaneous TOV-21G human ovarian cancer xenografts (n = 8) and on the transgenic adenocarcinoma of the mouse prostate (TRAMP) model (n = 7). Subsequently, the TOV-21G tumor-bearing mice were treated with bevacizumab and rescanned 2 days later. Pharmacokinetic analysis (extended Tofts model) was performed using either the first echo signal only (standard single-echo DCE-MRI) or the estimated signal at TE = 0 derived from exponential fitting of R2* relaxation (R2*-corrected). Neglecting R2* enhancement causes underestimation of Gd-DOTA concentration (peak enhancement underestimated by 9.4%-16% in TOV-21G tumors and 13%-20% in TRAMP prostates). Median K^trans and v_e were underestimated in every mouse (TOV-21G K^trans: 11%-19%, TOV-21G v_e: 5.3%-8.9%; TRAMP K^trans: 8.6%-19%, TRAMP v_e: 12%-21%). Bevacizumab treatment reduced K^trans in all TOV-21G tumors after 48 hours. Treatment effect was significantly greater in all tumors after R2* correction (median change of -0.050 min^-1 in R2*-corrected K^trans vs. -0.037 min^-1 in uncorrected K^trans). R2* enhancement in DCE-MRI is both time- and tissue-dependent and may not be negligible at 7 T in tissue with high K^trans. This has consequences for the use of K^trans and other DCE-MRI parameters as biomarkers, because treatment effect size can be underestimated when R2* enhancement is neglected.

Utility of Percentage Signal Recovery and Baseline Signal in DSC-MRI Optimized for Relative CBV Measurement for Differentiating Glioblastoma, Lymphoma, Metastasis, and Meningioma

BACKGROUND AND PURPOSE

The percentage signal recovery in non-leakage-corrected (no preload, high flip angle, intermediate TE) DSC-MR imaging is known to differ significantly for glioblastoma, metastasis, and primary CNS lymphoma. Because the percentage signal recovery is influenced by preload and pulse sequence parameters, we investigated whether the percentage signal recovery can still differentiate these common contrast-enhancing neoplasms using a DSC-MR imaging protocol designed for relative CBV accuracy (preload, intermediate flip angle, low TE).

MATERIALS AND METHODS

We retrospectively analyzed DSC-MR imaging of treatment-naïve, pathology-proved glioblastomas (n = 14), primary central nervous system lymphomas (n = 7), metastases (n = 20), and meningiomas (n = 13) using a protocol designed for relative CBV accuracy (a one-quarter-dose preload and single-dose bolus of gadobutrol, TR/TE = 1290/40 ms, flip angle = 60° at 1.5T). Mean percentage signal recovery, relative CBV, and normalized baseline signal intensity were compared within contrast-enhancing lesion volumes. Classification accuracy was determined by receiver operating characteristic analysis.

RESULTS

Relative CBV best differentiated meningioma from glioblastoma and from metastasis with areas under the curve of 0.84 and 0.82, respectively. The percentage signal recovery best differentiated primary central nervous system lymphoma from metastasis with an area under the curve of 0.81. Relative CBV and percentage signal recovery were similar in differentiating primary central nervous system lymphoma from glioblastoma and from meningioma. Although neither relative CBV nor percentage signal recovery differentiated glioblastoma from metastasis, mean normalized baseline signal intensity achieved 86% sensitivity and 50% specificity.

CONCLUSIONS

Similar to results for non-preload-based DSC-MR imaging, percentage signal recovery for one-quarter-dose preload-based, intermediate flip angle DSC-MR imaging differentiates most pair-wise comparisons of glioblastoma, metastasis, primary central nervous system lymphoma, and meningioma, except for glioblastoma versus metastasis. Differences in normalized post-preload baseline signal for glioblastoma and metastasis, reflecting a snapshot of dynamic contrast enhancement, may motivate the use of single-dose multiecho protocols permitting simultaneous quantification of DSC-MR imaging and dynamic contrast-enhanced MR imaging parameters.

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 T₁ and T 2 ∗ contrast agent leakage effects that result in inaccurate hemodynamic metrics. While multi-echo acquisitions remove T₁ 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 T 2 ∗ leakage effects remains compulsory for high rCBV accuracy.

Accurate Patient-Specific Machine Learning Models of Glioblastoma Invasion Using Transfer Learning

BACKGROUND AND PURPOSE

MR imaging-based modeling of tumor cell density can substantially improve targeted treatment of glioblastoma. Unfortunately, interpatient variability limits the predictive ability of many modeling approaches. We present a transfer learning method that generates individualized patient models, grounded in the wealth of population data, while also detecting and adjusting for interpatient variabilities based on each patient's own histologic data.

MATERIALS AND METHODS

We recruited patients with primary glioblastoma undergoing image-guided biopsies and preoperative imaging, including contrast-enhanced MR imaging, dynamic susceptibility contrast MR imaging, and diffusion tensor imaging. We calculated relative cerebral blood volume from DSC-MR imaging and mean diffusivity and fractional anisotropy from DTI. Following image coregistration, we assessed tumor cell density for each biopsy and identified corresponding localized MR imaging measurements. We then explored a range of univariate and multivariate predictive models of tumor cell density based on MR imaging measurements in a generalized one-model-fits-all approach. We then implemented both univariate and multivariate individualized transfer learning predictive models, which harness the available population-level data but allow individual variability in their predictions. Finally, we compared Pearson correlation coefficients and mean absolute error between the individualized transfer learning and generalized one-model-fits-all models.

RESULTS

Tumor cell density significantly correlated with relative CBV (r = 0.33, P < .001), and T1-weighted postcontrast (r = 0.36, P < .001) on univariate analysis after correcting for multiple comparisons. With single-variable modeling (using relative CBV), transfer learning increased predictive performance (r = 0.53, mean absolute error = 15.19%) compared with one-model-fits-all (r = 0.27, mean absolute error = 17.79%). With multivariate modeling, transfer learning further improved performance (r = 0.88, mean absolute error = 5.66%) compared with one-model-fits-all (r = 0.39, mean absolute error = 16.55%).

CONCLUSIONS

Transfer learning significantly improves predictive modeling performance for quantifying tumor cell density in glioblastoma.

Imaging vascular and hemodynamic features of the brain using dynamic susceptibility contrast and dynamic contrast enhanced MRI

In the context of neurologic disorders, dynamic susceptibility contrast (DSC) and dynamic contrast enhanced (DCE) MRI provide valuable insights into cerebral vascular function, integrity, and architecture. Even after two decades of use, these modalities continue to evolve as their biophysical and kinetic basis is better understood, with improvements in pulse sequences and accelerated imaging techniques and through application of more robust and automated data analysis strategies. Here, we systematically review each of these elements, with a focus on how their integration improves kinetic parameter accuracy and the development of new hemodynamic biomarkers that provide sub-voxel sensitivity (e.g., capillary transit time and flow heterogeneity). Regarding contrast mechanisms, we discuss the dipole-dipole interactions and susceptibility effects that give rise to simultaneous T₁, T₂ and T₂^∗ relaxation effects, including their quantification, influence on pulse sequence parameter optimization, and use in methods such as vessel size and vessel architectural imaging. The application of technologic advancements, such as parallel imaging, simultaneous multi-slice, undersampled k-space acquisitions, and sliding window strategies, enables improved spatial and/or temporal resolution of DSC and DCE acquisitions. Such acceleration techniques have also enabled the implementation of, clinically feasible, simultaneous multi-echo spin- and gradient echo acquisitions, providing more comprehensive and quantitative interrogation of T₁, T₂ and T₂^∗ changes. Characterizing these relaxation rate changes through different post-processing options allows for the quantification of hemodynamics and vascular permeability. The application of different biophysical models provides insight into traditional hemodynamic parameters (e.g., cerebral blood volume) and more advanced parameters (e.g., capillary transit time heterogeneity). We provide insight into the appropriate selection of biophysical models and the necessary post-processing steps to ensure reliable measurements while minimizing potential sources of error. We show representative examples of advanced DSC- and DCE-MRI methods applied to pathologic conditions affecting the cerebral microcirculation, including brain tumors, stroke, aging, and multiple sclerosis. The maturation and standardization of conventional DSC- and DCE-MRI techniques has enabled their increased integration into clinical practice and use in clinical trials, which has, in turn, spurred renewed interest in their technological and biophysical development, paving the way towards a more comprehensive assessment of cerebral hemodynamics.

Optimization of DSC MRI Echo Times for CBV Measurements Using Error Analysis in a Pilot Study of High-Grade Gliomas

BACKGROUND AND PURPOSE

The optimal TE must be calculated to minimize the variance in CBV measurements made with DSC MR imaging. Simulations can be used to determine the influence of the TE on CBV, but they may not adequately recapitulate the in vivo heterogeneity of precontrast T2*, contrast agent kinetics, and the biophysical basis of contrast agent-induced T2* changes. The purpose of this study was to combine quantitative multiecho DSC MRI T2* time curves with error analysis in order to compute the optimal TE for a traditional single-echo acquisition.

MATERIALS AND METHODS

Eleven subjects with high-grade gliomas were scanned at 3T with a dual-echo DSC MR imaging sequence to quantify contrast agent-induced T2* changes in this retrospective study. Optimized TEs were calculated with propagation of error analysis for high-grade glial tumors, normal-appearing white matter, and arterial input function estimation.

RESULTS

The optimal TE is a weighted average of the T2* values that occur as a contrast agent bolus transverses a voxel. The mean optimal TEs were 30.0 ± 7.4 ms for high-grade glial tumors, 36.3 ± 4.6 ms for normal-appearing white matter, and 11.8 ± 1.4 ms for arterial input function estimation (repeated-measures ANOVA, P < .001).

CONCLUSIONS

Greater heterogeneity was observed in the optimal TE values for high-grade gliomas, and mean values of all 3 ROIs were statistically significant. The optimal TE for the arterial input function estimation is much shorter; this finding implies that quantitative DSC MR imaging acquisitions would benefit from multiecho acquisitions. In the case of a single-echo acquisition, the optimal TE prescribed should be 30-35 ms (without a preload) and 20-30 ms (with a standard full-dose preload).

Characterizing the Influence of Preload Dosing on Percent Signal Recovery (PSR) and Cerebral Blood Volume (CBV) Measurements in a Patient Population With High-Grade Glioma Using Dynamic Susceptibility Contrast MRI

With DSC-MRI, contrast agent leakage effects in brain tumors can either be leveraged for percent signal recovery (PSR) measurements or be adequately resolved for accurate relative cerebral blood volume (rCBV) measurements. Leakage effects can be dimished by administration of a preload dose before imaging and/or specific postprocessing steps. This study compares the consistency of both PSR and rCBV measurements as a function of varying preload doses in a retrospective analysis of 14 subjects with high-grade gliomas. The scans consisted of 6 DSC-MRI scans during 6 sequential bolus injections (0.05 mmol/kg). Mean PSR was calculated for tumor and normal-appearing white matter regions of interest. DSC-MRI data were corrected for leakage effects before computing mean tumor rCBV. Statistical differences were seen across varying preloads for tumor PSR (P value = 4.57E-24). Tumor rCBV values did not exhibit statistically significant differences across preloads (P value = .14) and were found to be highly consistent for clinically relevant preloads (intraclass correlation coefficient = 0.93). For a 0.05 mmol/kg injection bolus and pulse sequence parameters used, the highest PSR contrast between normal-appearing white matter and tumor occurs when no preload is used. This suggests that studies using PSR as a biomarker should acquire DSC-MRI data without preload. The finding that leakage-corrected rCBV values do not depend on the presence or dose of preload contradicts that of previous studies with dissimilar acquisition protocols. This further confirms the sensitivity of rCBV to preload dosing schemes and pulse sequence parameters and highlights the importance of standardization efforts for achieving multisite rCBV consistency.

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 T₁ and T₂* 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.

Radiogenomics to characterize regional genetic heterogeneity in glioblastoma

BACKGROUND

Glioblastoma (GBM) exhibits profound intratumoral genetic heterogeneity. Each tumor comprises multiple genetically distinct clonal populations with different therapeutic sensitivities. This has implications for targeted therapy and genetically informed paradigms. Contrast-enhanced (CE)-MRI and conventional sampling techniques have failed to resolve this heterogeneity, particularly for nonenhancing tumor populations. This study explores the feasibility of using multiparametric MRI and texture analysis to characterize regional genetic heterogeneity throughout MRI-enhancing and nonenhancing tumor segments.

METHODS

We collected multiple image-guided biopsies from primary GBM patients throughout regions of enhancement (ENH) and nonenhancing parenchyma (so called brain-around-tumor, [BAT]). For each biopsy, we analyzed DNA copy number variants for core GBM driver genes reported by The Cancer Genome Atlas. We co-registered biopsy locations with MRI and texture maps to correlate regional genetic status with spatially matched imaging measurements. We also built multivariate predictive decision-tree models for each GBM driver gene and validated accuracies using leave-one-out-cross-validation (LOOCV).

RESULTS

We collected 48 biopsies (13 tumors) and identified significant imaging correlations (univariate analysis) for 6 driver genes: EGFR, PDGFRA, PTEN, CDKN2A, RB1, and TP53. Predictive model accuracies (on LOOCV) varied by driver gene of interest. Highest accuracies were observed for PDGFRA (77.1%), EGFR (75%), CDKN2A (87.5%), and RB1 (87.5%), while lowest accuracy was observed in TP53 (37.5%). Models for 4 driver genes (EGFR, RB1, CDKN2A, and PTEN) showed higher accuracy in BAT samples (n = 16) compared with those from ENH segments (n = 32).

CONCLUSION

MRI and texture analysis can help characterize regional genetic heterogeneity, which offers potential diagnostic value under the paradigm of individualized oncology.

On the Use of DSC-MRI for Measuring Vascular Permeability

BACKGROUND AND PURPOSE

Contrast agent extravasation has been shown to confound brain tumor perfusion measurements with DSC-MR imaging, necessitating the use of correction techniques (eg, Weisskoff, Bjornerud). Leakage parameters (K2 and K(a)) postulated to reflect vessel permeability can be extracted from these correction methods; however, the biophysical interpretation of these parameters and their relationship to commonly used MR imaging measures of vascular permeability (eg, contrast agent volume transfer constant, [K(trans)]) remain unclear. Given that vascular density, as assessed by blood volume, and vascular permeability, as reflected by K(trans) (and potentially K2 or K(a)), report on unique and clinically informative vascular characteristics, there is a compelling interest to simultaneously assess these features.

MATERIALS AND METHODS

We acquired multiecho DSC-MR imaging data, allowing the simultaneous computation and voxelwise comparison of single- and dual-echo derived measures of K2, K(a) and K(trans) in patients with glioma. This acquisition enabled the investigation of competing T1 and T2* leakage effects and TE dependency on these parameters.

RESULTS

K2 and K(a) displayed nonsignificant (P = .150 and P = .060, respectively) voxelwise linear correlations with K(trans), while a significant (P < .001) inverse relationship was observed between K2 and Ka (coefficient of determination [r(2)] = 0.466-0.984). Significantly different (P < .005) mean estimates were found between voxels exhibiting predominately T1 and T2* effects for K2 and K(a). K(trans), however, was observed to be similar between these voxels (0.109 versus 0.092 minutes(-1)). Significant differences (P < .001) in extracellular-extravascular volume fraction (v(e)) (0.285 versus 0.167) were also observed between cohorts. Additionally, K2 and K(a) were found to have a significant quadratic relationship (P = .031 and P = .005, respectively) with v(e).

CONCLUSIONS

Estimates of vascular permeability in brain tumors may be simultaneously acquired from multiple-echo DSC-MR imaging via K(trans); however, caution should be used in assuming a similar relationship for K2 and K(a).

Multi-Parametric MRI and Texture Analysis to Visualize Spatial Histologic Heterogeneity and Tumor Extent in Glioblastoma

Background

Genetic profiling represents the future of neuro-oncology but suffers from inadequate biopsies in heterogeneous tumors like Glioblastoma (GBM). Contrast-enhanced MRI (CE-MRI) targets enhancing core (ENH) but yields adequate tumor in only ~60% of cases. Further, CE-MRI poorly localizes infiltrative tumor within surrounding non-enhancing parenchyma, or brain-around-tumor (BAT), despite the importance of characterizing this tumor segment, which universally recurs. In this study, we use multiple texture analysis and machine learning (ML) algorithms to analyze multi-parametric MRI, and produce new images indicating tumor-rich targets in GBM.

Methods

We recruited primary GBM patients undergoing image-guided biopsies and acquired pre-operative MRI: CE-MRI, Dynamic-Susceptibility-weighted-Contrast-enhanced-MRI, and Diffusion Tensor Imaging. Following image coregistration and region of interest placement at biopsy locations, we compared MRI metrics and regional texture with histologic diagnoses of high- vs low-tumor content (≥80% vs <80% tumor nuclei) for corresponding samples. In a training set, we used three texture analysis algorithms and three ML methods to identify MRI-texture features that optimized model accuracy to distinguish tumor content. We confirmed model accuracy in a separate validation set.

Results

We collected 82 biopsies from 18 GBMs throughout ENH and BAT. The MRI-based model achieved 85% cross-validated accuracy to diagnose high- vs low-tumor in the training set (60 biopsies, 11 patients). The model achieved 81.8% accuracy in the validation set (22 biopsies, 7 patients).

Conclusion

Multi-parametric MRI and texture analysis can help characterize and visualize GBM’s spatial histologic heterogeneity to identify regional tumor-rich biopsy targets.

Validation of a T1 and T2* leakage correction method based on multiecho dynamic susceptibility contrast MRI using MION as a reference standard

PURPOSE

A combined biophysical- and pharmacokinetic-based method is proposed to separate, quantify, and correct for both T1 and T2* leakage effects using dual-echo dynamic susceptibility contrast (DSC) acquisitions to provide more accurate hemodynamic measures, as validated by a reference intravascular contrast agent (CA).

THEORY AND METHODS

Dual-echo DSC-MRI data were acquired in two rodent glioma models. The T1 leakage effects were removed and also quantified to subsequently correct for the remaining T2* leakage effects. Pharmacokinetic, biophysical, and combined biophysical and pharmacokinetic models were used to obtain corrected cerebral blood volume (CBV) and cerebral blood flow (CBF), and these were compared with CBV and CBF from an intravascular CA.

RESULTS

T1 -corrected CBV was significantly overestimated compared with MION CBV, while T1 + T2*-correction yielded CBV values closer to the reference values. The pharmacokinetic and simplified biophysical methods showed similar results and underestimated CBV in tumors exhibiting strong T2* leakage effects. The combined method was effective for correcting T1 and T2* leakage effects across tumor types.

CONCLUSION

Correcting for both T1 and T2* leakage effects yielded more accurate measures of CBV. The combined correction method yields more reliable CBV measures than either correction method alone, but for certain brain tumor types (e.g., gliomas), the simplified biophysical method may provide a robust and computationally efficient alternative. Magn Reson Med 76:613-625, 2016. © 2015 Wiley Periodicals, Inc.