Repeatability of dynamic contrast enhanced vp parameter in healthy subjects and patients with brain tumors

DCE-MRI quantitative analysis and MRI-based radiomics for predicting the early efficacy of microwave ablation in lung cancers

OBJECTIVES

To evaluate the feasibility and value of dynamic contrast-enhanced MRI (DCE-MRI) quantitative analysis and MRI-based radiomics in predicting the efficacy of microwave ablation (MWA) in lung cancers (LCs).

METHODS

Forty-three patients with LCs who underwent DCE-MRI within 24 h of receiving MWA were enrolled in the study and divided into two groups according to the modified response evaluation criteria in solid tumors (m-RECIST) criteria: the effective treatment (complete response + partial response + stable disease, n = 28) and the ineffective treatment (progressive disease, n = 15). DCE-MRI datasets were processed by Omni. Kinetics software, using the extended tofts model (ETM) and exchange model (ECM) to yield pharmacokinetic parameters and their histogram features. Changes in quantitative perfusion parameters were compared between the two groups. Scientific research platform ( https://medresearch.shukun.net/ ) was used for radiomics analysis. A total of 1874 radiomic features were extracted for each tumor by manually segmentation of T1WI and Contrast-enhanced of T1WI (Ce-T1WI) fat inhibition sequence. The performances of radiomics models were evaluated by the receiver operating characteristic curve. Based on radiomics features, survival curves were generated by Kaplan-Meier survival analysis to evaluate patient outcomes. P < 0.05 was set for the significance threshold.

RESULTS

The Vp value of ECM was significantly higher in the ineffective group compared to the effective group (p = 0.027). Additionally, the skewness, and kurtosis of Vp (p = 0.020 and 0.013, respectively) derived from ETM and Fp (p = 0.027 and 0.030, respectively) from ECM as well as the quantiles were higher in the ineffective group than in the effective group. Significant statistical differences were observed in the energy and entropy of Ve (p = 0.044 and 0.025, respectively) and Vp (p = 0.025 and 0.026, respectively) between the two groups. In the process of model construction, seven features from T1WI, five features from Ce-T1WI, and ten features from combined sequences were ultimately selected. The area under the curve (AUC) values for the T1WI model, Ce-T1WI model, and combined model were 0.910, 0.890, 0.965 in the training group, and 0.850, 0.700, 0.850 in the test group, respectively.

CONCLUSIONS

DCE-MRI quantitative analysis and MRI-based radiomics may be helpful in assessing the early response to MWA in patients with LCs.

Greybox: A hybrid algorithm for direct estimation of tracer kinetic parameters from undersampled DCE-MRI data

BACKGROUND

A variety of deep learning-based and iterative approaches are available to predict Tracer Kinetic (TK) parameters from fully sampled or undersampled dynamic contrast-enhanced (DCE) MRI data. However, both the methods offer distinct benefits and drawbacks.

PURPOSE

To propose a hybrid algorithm (named as 'Greybox'), using both model- as well as DL-based, for solving a multi-parametric non-linear inverse problem of directly estimating TK parameters from undersampled DCE MRI data, which is invariant to undersampling rate.

METHODS

The proposed algorithm was inspired by plug-and-play algorithms used for solving linear inverse imaging problems. This technique was tested for its effectiveness in solving the nonlinear ill-posed inverse problem of generating 3D TK parameter maps from four-dimensional (4D; Spatial + Temporal) retrospectively undersampled k-space data. The algorithm learns a deep learning-based prior using UNET to estimate theK trans $\mathbf {K_{trans}}$ andV p $\mathbf {V_{p}}$ parameters based on the Patlak pharmacokinetic model, and this trained prior was utilized to estimate the TK parameter maps using an iterative gradient-based optimization scheme. Unlike the existing DL models, this network is invariant to the undersampling rate of the input data. The proposed method was compared with the total variation-based direct reconstruction technique on brain, breast, and prostate DCE-MRI datasets for various undersampling rates using the Radial Golden Angle (RGA) scheme. For the breast dataset, an indirect estimation using the Fast Composite Splitting algorithm was utilized for comparison. Undersampling rates of 8 × $\times$ , 12 × $\times$ and 20 × $\times$ were used for the experiments, and the results were compared using the PSNR and SSIM as metrics. For the breast dataset of 10 patients, data from four patients were utilized for training (1032 samples), two for validation (752 samples), and the entire volume of four patients for testing. Similarly, for the prostate dataset of 18 patients, 10 patients were utilized for training (720 samples), five for validation (216 samples), and the whole volume of three patients for testing. For the brain dataset of nineteen patients, ten patients were used for training (3152 samples), five for validation (1168 samples), and the whole volume of four patients for testing. Statistical tests were also conducted to assess the significance of the improvement in performance.

RESULTS

The experiments showed that the proposed Greybox performs significantly better than other direct reconstruction methods. The proposed algorithm improved the estimatedK trans $\mathbf {K_{trans}}$ andV p $\mathbf {V_{p}}$ in terms of the peak signal-to-noise ratio by up to 3 dB compared to other standard reconstruction methods.

CONCLUSION

The proposed hybrid reconstruction algorithm, Greybox, can provide state-of-the-art performance in solving the nonlinear inverse problem of DCE-MRI. This is also the first of its kind to utilize convolutional neural network-based encodings as part of the plug-and-play priors to improve the performance of the reconstruction algorithm.

Mathematical modelling of the dynamics of image-informed tumor habitats in a murine model of glioma

Tumors exhibit high molecular, phenotypic, and physiological heterogeneity. In this effort, we employ quantitative magnetic resonance imaging (MRI) data to capture this heterogeneity through imaging-based subregions or "habitats" in a murine model of glioma. We then demonstrate the ability to model and predict the growth of the habitats using coupled ordinary differential equations (ODEs) in the presence and absence of radiotherapy. Female Wistar rats (N = 21) were inoculated intracranially with 10⁶ C6 glioma cells, a subset of which received 20 Gy (N = 5) or 40 Gy (N = 8) of radiation. All rats underwent diffusion-weighted and dynamic contrast-enhanced MRI at up to seven time points. All MRI data at each visit were subsequently clustered using k-means to identify physiological tumor habitats. A family of four models consisting of three coupled ODEs were developed and calibrated to the habitat time series of control and treated rats and evaluated for predictive capability. The Akaike Information Criterion was used for model selection, and the normalized sum-of-square-error (SSE) was used to evaluate goodness-of-fit in model calibration and prediction. Three tumor habitats with significantly different imaging data characteristics (p < 0.05) were identified: high-vascularity high-cellularity, low-vascularity high-cellularity, and low-vascularity low-cellularity. Model selection resulted in a five-parameter model whose predictions of habitat dynamics yielded SSEs that were similar to the SSEs from the calibrated model. It is thus feasible to mathematically describe habitat dynamics in a preclinical model of glioma using biology-based ODEs, showing promise for forecasting heterogeneous tumor behavior.

High-Grade Glioma Treatment Response Monitoring Biomarkers: A Position Statement on the Evidence Supporting the Use of Advanced MRI Techniques in the Clinic, and the Latest Bench-to-Bedside Developments. Part 1: Perfusion and Diffusion Techniques

Objective

Summarize evidence for use of advanced MRI techniques as monitoring biomarkers in the clinic, and highlight the latest bench-to-bedside developments.

Methods

Experts in advanced MRI techniques applied to high-grade glioma treatment response assessment convened through a European framework. Current evidence regarding the potential for monitoring biomarkers in adult high-grade glioma is reviewed, and individual modalities of perfusion, permeability, and microstructure imaging are discussed (in Part 1 of two). In Part 2, we discuss modalities related to metabolism and/or chemical composition, appraise the clinic readiness of the individual modalities, and consider post-processing methodologies involving the combination of MRI approaches (multiparametric imaging) or machine learning (radiomics).

Results

High-grade glioma vasculature exhibits increased perfusion, blood volume, and permeability compared with normal brain tissue. Measures of cerebral blood volume derived from dynamic susceptibility contrast-enhanced MRI have consistently provided information about brain tumor growth and response to treatment; it is the most clinically validated advanced technique. Clinical studies have proven the potential of dynamic contrast-enhanced MRI for distinguishing post-treatment related effects from recurrence, but the optimal acquisition protocol, mode of analysis, parameter of highest diagnostic value, and optimal cut-off points remain to be established. Arterial spin labeling techniques do not require the injection of a contrast agent, and repeated measurements of cerebral blood flow can be performed. The absence of potential gadolinium deposition effects allows widespread use in pediatric patients and those with impaired renal function. More data are necessary to establish clinical validity as monitoring biomarkers. Diffusion-weighted imaging, apparent diffusion coefficient analysis, diffusion tensor or kurtosis imaging, intravoxel incoherent motion, and other microstructural modeling approaches also allow treatment response assessment; more robust data are required to validate these alone or when applied to post-processing methodologies.

Conclusion

Considerable progress has been made in the development of these monitoring biomarkers. Many techniques are in their infancy, whereas others have generated a larger body of evidence for clinical application.

High-Grade Glioma Treatment Response Monitoring Biomarkers: A Position Statement on the Evidence Supporting the Use of Advanced MRI Techniques in the Clinic, and the Latest Bench-to-Bedside Developments. Part 2: Spectroscopy, Chemical Exchange Saturation, Multiparametric Imaging, and Radiomics

Objective

To summarize evidence for use of advanced MRI techniques as monitoring biomarkers in the clinic, and to highlight the latest bench-to-bedside developments.

Methods

The current evidence regarding the potential for monitoring biomarkers was reviewed and individual modalities of metabolism and/or chemical composition imaging discussed. Perfusion, permeability, and microstructure imaging were similarly analyzed in Part 1 of this two-part review article and are valuable reading as background to this article. We appraise the clinic readiness of all the individual modalities and consider methodologies involving machine learning (radiomics) and the combination of MRI approaches (multiparametric imaging).

Results

The biochemical composition of high-grade gliomas is markedly different from healthy brain tissue. Magnetic resonance spectroscopy allows the simultaneous acquisition of an array of metabolic alterations, with choline-based ratios appearing to be consistently discriminatory in treatment response assessment, although challenges remain despite this being a mature technique. Promising directions relate to ultra-high field strengths, 2-hydroxyglutarate analysis, and the use of non-proton nuclei. Labile protons on endogenous proteins can be selectively targeted with chemical exchange saturation transfer to give high resolution images. The body of evidence for clinical application of amide proton transfer imaging has been building for a decade, but more evidence is required to confirm chemical exchange saturation transfer use as a monitoring biomarker. Multiparametric methodologies, including the incorporation of nuclear medicine techniques, combine probes measuring different tumor properties. Although potentially synergistic, the limitations of each individual modality also can be compounded, particularly in the absence of standardization. Machine learning requires large datasets with high-quality annotation; there is currently low-level evidence for monitoring biomarker clinical application.

Conclusion

Advanced MRI techniques show huge promise in treatment response assessment. The clinical readiness analysis highlights that most monitoring biomarkers require standardized international consensus guidelines, with more facilitation regarding technique implementation and reporting in the clinic.

A systematic review on the use of quantitative imaging to detect cancer therapy adverse effects in normal-appearing brain tissue

Cancer therapy for both central nervous system (CNS) and non-CNS tumors has been previously associated with transient and long-term cognitive deterioration, commonly referred to as 'chemo fog'. This therapy-related damage to otherwise normal-appearing brain tissue is reported using post-mortem neuropathological analysis. Although the literature on monitoring therapy effects on structural magnetic resonance imaging (MRI) is well established, such macroscopic structural changes appear relatively late and irreversible. Early quantitative MRI biomarkers of therapy-induced damage would potentially permit taking these treatment side effects into account, paving the way towards a more personalized treatment planning.This systematic review (PROSPERO number 224196) provides an overview of quantitative tomographic imaging methods, potentially identifying the adverse side effects of cancer therapy in normal-appearing brain tissue. Seventy studies were obtained from the MEDLINE and Web of Science databases. Studies reporting changes in normal-appearing brain tissue using MRI, PET, or SPECT quantitative biomarkers, related to radio-, chemo-, immuno-, or hormone therapy for any kind of solid, cystic, or liquid tumor were included. The main findings of the reviewed studies were summarized, providing also the risk of bias of each study assessed using a modified QUADAS-2 tool. For each imaging method, this review provides the methodological background, and the benefits and shortcomings of each method from the imaging perspective. Finally, a set of recommendations is proposed to support future research.

Repeatability of tumor perfusion kinetics from dynamic contrast-enhanced MRI in glioblastoma

Background

Dynamic contrast-enhanced MRI (DCE-MRI) parameters have been shown to be biomarkers for treatment response in glioblastoma (GBM). However, variations in analysis and measurement methodology complicate determination of biological changes measured via DCE. The aim of this study is to quantify DCE-MRI variations attributable to analysis methodology and image quality in GBM patients.

Methods

The Extended Tofts model (eTM) and Leaky Tracer Kinetic Model (LTKM), with manually and automatically segmented vascular input functions (VIFs), were used to calculate perfusion kinetic parameters from 29 GBM patients with double-baseline DCE-MRI data. DCE-MRI images were acquired 2-5 days apart with no change in treatment. Repeatability of kinetic parameters was quantified with Bland-Altman and percent repeatability coefficient (%RC) analysis.

Results

The perfusion parameter with the least RC was the plasma volume fraction (v _p ), with a %RC of 53%. The extra-cellular extra-vascular volume fraction (v _e ) %RC was 82% and 81%, for extended Tofts-Kety Model (eTM) and LTKM respectively. The %RC of the volume transfer rate constant (K ^trans ) was 72% for the eTM, and 82% for the LTKM, respectively. Using an automatic VIF resulted in smaller %RCs for all model parameters, as compared to manual VIF.

Conclusions

As much as 72% change in K ^trans (eTM, autoVIF) can be attributable to non-biological changes in the 2-5 days between double-baseline imaging. Poor K ^trans repeatability may result from inferior temporal resolution and short image acquisition time. This variation suggests DCE-MRI repeatability studies should be performed institutionally, using an automatic VIF method and following quantitative imaging biomarkers alliance guidelines.

Quantitative magnetic resonance imaging towards clinical application in multiple sclerosis

Quantitative MRI provides biophysical measures of the microstructural integrity of the CNS, which can be compared across CNS regions, patients, and centres. In patients with multiple sclerosis, quantitative MRI techniques such as relaxometry, myelin imaging, magnetization transfer, diffusion MRI, quantitative susceptibility mapping, and perfusion MRI, complement conventional MRI techniques by providing insight into disease mechanisms. These include: (i) presence and extent of diffuse damage in CNS tissue outside lesions (normal-appearing tissue); (ii) heterogeneity of damage and repair in focal lesions; and (iii) specific damage to CNS tissue components. This review summarizes recent technical advances in quantitative MRI, existing pathological validation of quantitative MRI techniques, and emerging applications of quantitative MRI to patients with multiple sclerosis in both research and clinical settings. The current level of clinical maturity of each quantitative MRI technique, especially regarding its integration into clinical routine, is discussed. We aim to provide a better understanding of how quantitative MRI may help clinical practice by improving stratification of patients with multiple sclerosis, and assessment of disease progression, and evaluation of treatment response.

DCE-MRI in Glioma, Infiltration Zone and Healthy Brain to Assess Angiogenesis: A Biopsy Study

Purpose

To explore the focal predictability of vascular growth factor expression and neovascularization using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) in glioma.

Methods

120 brain biopsies were taken in vital tumor, infiltration zone and normal brain tissue of 30 glioma patients: 17 IDH(isocitrate dehydrogenase)-wildtype glioblastoma (GBM), 1 IDH-wildtype astrocytoma °III (together prognostic group 1), 3 IDH-mutated GBM (group 2), 3 anaplastic astrocytomas IDH-mutated (group 3), 4 anaplastic oligodendrogliomas and 2 low-grade oligodendrogliomas (together prognostic group 4). A mixed linear model evaluated the predictabilities of microvessel density (MVD), vascular area ratio (VAR), mean vessel size (MVS), vascular endothelial growth factor and receptors (VEGF-A, VEGFR‑2) and vascular endothelial-protein tyrosine phosphatase (VE-PTP) expression from Tofts model kinetic and model-free curve parameters.

Results

All kinetic parameters were associated with VEGF‑A (all p < 0.001) expression. K^trans, k_ep and v_e were associated with VAR (p = 0.006, 0.004 and 0.01, respectively) and MVS (p = 0.0001, 0.02 and 0.003, respectively) but not MVD (p = 0.84, 0.74 and 0.73, respectively). Prognostic groups differed in K^trans (p = 0.007) and v_e (p = 0.004) values measured in the infiltration zone. Despite significant differences of VAR, MVS, VEGF‑A, VEGFR‑2, and VE-PTP in vital tumor tissue and the infiltration zone (p = 0.0001 for all), there was no significant difference between kinetic parameters measured in these zones.

Conclusion

The DCE-MRI kinetic parameters show correlations with microvascular parameters in vital tissue and also reveal blood-brain barrier abnormalities in the infiltration zones adequate to differentiate glioma prognostic groups.

Supplementary Information

The online version of this article (10.1007/s00062-021-01015-3) contains supplementary material, which is available to authorized users.