Are complex DCE-MRI models supported by clinical data?

Dynamic Contrast-Enhanced (DCE) MRI

The non-invasive dynamic contrast-enhanced MRI (DCE-MRI) method provides valuable insights into tissue perfusion and vascularity. Primarily used in oncology, DCE-MRI is typically utilized to assess morphology and contrast agent (CA) kinetics in the tissue of interest. Interpretation of the temporal signatures of DCE-MRI data includes qualitative, semi-quantitative, and quantitative approaches. Recent advances in MRI technology allow simultaneous high spatial and temporal resolutions in DCE-MRI data acquisition on most vendor platforms, enabling the more desirable approach of quantitative data analysis using pharmacokinetic (PK) modeling. Many technical factors, including signal-to-noise ratio, temporal resolution, quantifications of arterial input function and native tissue T1, and PK model selection, need to be carefully considered when performing quantitative DCE-MRI. Standardization in data acquisition and analysis is especially important in multi-center studies.

Parametric maps of spatial two-tissue compartment model for prostate dynamic contrast enhanced MRI - comparison with the standard tofts model in the diagnosis of prostate cancer

The spatial two-tissue compartment model (2TCM) was used to analyze prostate dynamic contrast enhanced (DCE) MRI data and compared with the standard Tofts model. A total of 29 patients with biopsy-confirmed prostate cancer were included in this IRB-approved study. MRI data were acquired on a Philips Achieva 3T-TX scanner. After T2-weighted and diffusion-weighted imaging, DCE data using 3D T1-FFE mDIXON sequence were acquired pre- and post-contrast media injection (0.1 mmol/kg Multihance) for 60 dynamic scans with temporal resolution of 8.3 s/image. The 2TCM has one fast ([Formula: see text] and [Formula: see text]) and one slow ([Formula: see text] and [Formula: see text]) exchanging compartment, compared with the standard Tofts model parameters (Ktrans and kep). On average, prostate cancer had significantly higher values (p < 0.01) than normal prostate tissue for all calculated parameters. There was a strong correlation (r = 0.94, p < 0.001) between Ktrans and [Formula: see text] for cancer, but weak correlation (r = 0.28, p < 0.05) between kep and [Formula: see text]. Average root-mean-square error (RMSE) in fits from the 2TCM was significantly smaller (p < 0.001) than the RMSE in fits from the Tofts model. Receiver operating characteristic (ROC) analysis showed that fast [Formula: see text] had the highest area under the curve (AUC) than any other individual parameter. The combined four parameters from the 2TCM had a considerably higher AUC value than the combined two parameters from the Tofts model. The 2TCM is useful for quantitative analysis of prostate DCE-MRI data and provides new information in the diagnosis of prostate cancer.

Advanced MR Techniques for Preoperative Glioma Characterization: Part 1

Preoperative clinical magnetic resonance imaging (MRI) protocols for gliomas, brain tumors with dismal outcomes due to their infiltrative properties, still rely on conventional structural MRI, which does not deliver information on tumor genotype and is limited in the delineation of diffuse gliomas. The GliMR COST action wants to raise awareness about the state of the art of advanced MRI techniques in gliomas and their possible clinical translation or lack thereof. This review describes current methods, limits, and applications of advanced MRI for the preoperative assessment of glioma, summarizing the level of clinical validation of different techniques. In this first part, we discuss dynamic susceptibility contrast and dynamic contrast-enhanced MRI, arterial spin labeling, diffusion-weighted MRI, vessel imaging, and magnetic resonance fingerprinting. The second part of this review addresses magnetic resonance spectroscopy, chemical exchange saturation transfer, susceptibility-weighted imaging, MRI-PET, MR elastography, and MR-based radiomics applications. Evidence Level: 3 Technical Efficacy: Stage 2.

Impact of Arterial Input Function and Pharmacokinetic Models on DCE-MRI Biomarkers for Detection of Vascular Effect Induced by Stroma-Directed Drug in an Orthotopic Mouse Model of Pancreatic Cancer

PURPOSE

We demonstrated earlier in mouse models of pancreatic ductal adenocarcinoma (PDA) that K^trans derived from dynamic contrast-enhanced (DCE) MRI detected microvascular effect induced by PEGPH20, a hyaluronidase which removes stromal hyaluronan, leading to reduced interstitial fluid pressure in the tumor (Clinical Cancer Res (2019) 25: 2314-2322). How the choice of pharmacokinetic (PK) model and arterial input function (AIF) may impact DCE-derived markers for detecting such an effect is not known.

PROCEDURES

Retrospective analyses of the DCE-MRI of the orthotopic PDA model are performed to examine the impact of individual versus group AIF combined with Tofts model (TM), extended-Tofts model (ETM), or shutter-speed model (SSM) on the ability to detect the microvascular changes induced by PEGPH20 treatment.

RESULTS

Individual AIF exhibit a marked difference in peak gadolinium concentration. However, across all three PK models, k_ep values show a significant correlation between individual versus group-AIF (p < 0.01). Regardless individual or group AIF, when k_ep is obtained from fitting the DCE-MRI data using the SSM, k_ep shows a significant increase after PEGPH20 treatment (p < 0.05 compared to the baseline); %change of k_ep from baseline to post-treatment is also significantly different between PEGPH20 versus vehicle group (p < 0.05). In comparison, when k_ep is derived from the TM, only the use of individual AIF leads to a significant increase of k_ep after PEGPH20 treatment, whereas the %change of k_ep is not different between PEGPH20 versus vehicle group. Group AIF but not individual AIF allows detection of a significant increase of V_p (derived from the ETM) in PEGPH20 versus vehicle group (p < 0.05). Increase of V_p is consistent with a large increase of mean capillary lumen area estimated from immunostaining.

CONCLUSION

Our results suggest that k_ep derived from SSM and V_p from ETM, both using group AIF, are optimal for the detection of microvascular changes induced by stroma-directed drug PEGPH20. These analyses provide insights in the choice of PK model and AIF for optimal DCE protocol design in mouse pancreatic cancer models.

Surrogate vascular input function measurements from the superior sagittal sinus are repeatable and provide tissue-validated kinetic parameters in brain DCE-MRI

Accurate vascular input function (VIF) derivation is essential in brain dynamic contrast-enhanced (DCE) MRI. The optimum site for VIF estimation is, however, debated. This study sought to compare VIFs extracted from the internal carotid artery (ICA) and its branches with an arrival-corrected vascular output function (VOF) derived from the superior sagittal sinus (VOF_SSS). DCE-MRI datasets from sixty-six patients with different brain tumours were retrospectively analysed and plasma gadolinium-based contrast agent (GBCA) concentration-time curves used to extract VOF/VIFs from the SSS, the ICA, and the middle cerebral artery. Semi-quantitative parameters across each first-pass VOF/VIF were compared and the relationship between these parameters and GBCA dose was evaluated. Through a test-retest study in 12 patients, the repeatability of each semiquantitative VOF/VIF parameter was evaluated; and through comparison with histopathological data the accuracy of kinetic parameter estimates derived using each VOF/VIF and the extended Tofts model was also assessed. VOF_SSS provided a superior surrogate global input function compared to arteries, with greater contrast-to-noise (p < 0.001), higher peak (p < 0.001, repeated-measures ANOVA), and a greater sensitivity to interindividual plasma GBCA concentration. The repeatability of VOF_SSS derived semi-quantitative parameters was good to excellent (ICC = 0.717-0.888) outperforming arterial based approaches. In contrast to arterial VIFs, kinetic parameters obtained using a SSS derived VOF permitted detection of intertumoural differences in both microvessel surface area and cell density within resected tissue specimens. These results support the usage of an arrival-corrected VOF_SSS as a surrogate vascular input function for kinetic parameter mapping in brain DCE-MRI.

Rat Brain Global Ischemia-Induced Diffusion Changes Revisited: Biophysical Modeling of the Water and NAA MR "Diffusion Signal"

PURPOSE

To assess changes in intracellular diffusion as a mechanism for the reduction in water ADC that accompanies brain injury. Using NAA as a marker of neuronal cytoplasmic diffusion, NAA diffusion was measured before and after global ischemia (immediately postmortem) in the female Sprague-Dawley rat.

METHODS

Diffusion-weighted PRESS spectra, with diffusion encoding in a single direction, were acquired from large voxels of rat brain gray matter in vivo and postischemia employing either pairs of pulsed half-sine-shaped gradients (in vivo and postischemia, b_max  = 19 ms/μm² ) or sinusoidal oscillating gradients (in vivo only) with frequencies of 99.2-250 Hz. A 2D randomly oriented cylinder (neurite) model gave estimates of longitudinal and transverse diffusivities (D_L and D_T , respectively). In this model, D_L represents the "free" diffusivity of NAA, whereas D_T reflects highly restricted diffusion. Using oscillating gradients, the frequency dependence of D_T [D_T (ω)] gave estimates of the cylinder (axon/dendrite) radius.

RESULTS

A 10% decrease in D_L,NAA followed global ischemia, dropping from 0.391 ± 0.012 μm² /ms to 0.350 ± 0.009 μm² /ms. Modeling D_T,NAA (ω) provided an estimate of the neurite radius of 1.0 ± 0.6 μm.

CONCLUSION

Whereas the increase in apparent intraneuronal viscosity suggested by changes in D_L,NAA may contribute to the overall reduction in water ADC associated with brain injury, it is not sufficient to be the sole explanation. Estimates of neurite radius based on D_T (ω) were consistent with literature values.

Study on the Value of DCE-MRI in Differentiating Glossitis and Tongue Cancer and the Intratumour Heterogeneity

Objective

DCE-MRI is an imaging technique that reflects the blood perfusion status of the tissue’s microcirculation. The purpose of this article is to explore the clinical value of dynamic contrast-enhanced (DCE)-MRI in distinguishing benign and malignant tongue lesions and the internal heterogeneity of a tumour.

Methods

The patients were divided into a tongue cancer group (22 patients) and a glossitis group (7 patients) based on the pathology results. All of the patients underwent DCE-MRI examination.

Results

The results of this study showed that the volume transfer constant (Ktrans), rate constant (Kep), contrast enhancement ratio (CER) and initial area under the gadolinium contrast agent concentration time curve (IAUGG) values of the tongue cancer group were significantly higher than those of the glossitis group, and the difference was statistically significant (P < 0.05). However, the extravascular extracellular volume fraction (Ve), fractional plasma volume (fPV), maximum slope (MaxSlope), and bolus arrival time (BAT) values measured by DCE-MRI in the tongue cancer group were not significantly different from those in the glossitis group (P > 0.05). The results of this study showed that the Ktrans, Kep, and IAUGG values measured by DCE-MRI had a good ability to distinguish tongue inflammation from tumours and Ktrans threshold of 0.484 has the best discriminative ability among them. The mean Ktrans values of stage I–II lesions were significantly higher than that of stage III–IV lesion (p = 0.045).

Conclusion

DCE-MRI is effective in distinguishing between benign and malignant tongue lesions and the internal heterogeneity of the tumour; it is worth following up in a larger study.

Clinical Registration Number

Research registry 6393.

Non-Invasive Evaluation of Cerebral Microvasculature Using Pre-Clinical MRI: Principles, Advantages and Limitations

Alterations to the cerebral microcirculation have been recognized to play a crucial role in the development of neurodegenerative disorders. However, the exact role of the microvascular alterations in the pathophysiological mechanisms often remains poorly understood. The early detection of changes in microcirculation and cerebral blood flow (CBF) can be used to get a better understanding of underlying disease mechanisms. This could be an important step towards the development of new treatment approaches. Animal models allow for the study of the disease mechanism at several stages of development, before the onset of clinical symptoms, and the verification with invasive imaging techniques. Specifically, pre-clinical magnetic resonance imaging (MRI) is an important tool for the development and validation of MRI sequences under clinically relevant conditions. This article reviews MRI strategies providing indirect non-invasive measurements of microvascular changes in the rodent brain that can be used for early detection and characterization of neurodegenerative disorders. The perfusion MRI techniques: Dynamic Contrast Enhanced (DCE), Dynamic Susceptibility Contrast Enhanced (DSC) and Arterial Spin Labeling (ASL), will be discussed, followed by less established imaging strategies used to analyze the cerebral microcirculation: Intravoxel Incoherent Motion (IVIM), Vascular Space Occupancy (VASO), Steady-State Susceptibility Contrast (SSC), Vessel size imaging, SAGE-based DSC, Phase Contrast Flow (PC) Quantitative Susceptibility Mapping (QSM) and quantitative Blood-Oxygenation-Level-Dependent (qBOLD). We will emphasize the advantages and limitations of each strategy, in particular on applications for high-field MRI in the rodent's brain.

A free time point model for dynamic contrast enhanced exploration

Dynamic-Contrast-Enhanced (DCE) Imaging has been widely studied to characterize microcirculatory disorders associated with various diseases. Although numerous studies have demonstrated its diagnostic interest, the physiological interpretation using pharmacokinetic models often remains debatable. Indeed, to be interpretable, a model must provide, at first instance, an accurate description of the DCE data. However, the evaluation and optimization of this accuracy remain rather limited in DCE. Here we established a non-linear Free-Time-Point-Hermite (FTPH) data-description model designed to fit DCE data accurately. Its performance was evaluated on data generated using two contrasting pharmacokinetic microcirculatory hypotheses (MH). The accuracy of data description of the models was evaluated by calculating the mean squared error (QE) from initial and assessed tissue impulse responses. Then, FTPH assessments were provided to blinded observers to evaluate if these assessments allowed observers to identify MH in their data. Regardless of the initial pharmacokinetic model used for data generation, QE was lower than 3% for the noise-free datasets and increased up to 10% for a signal-to-noise-ratio (SNR) of 20. Under SNR = 20, the sensitivity and specificity of the MH identification were over 80%. The performance of the FTPH model was higher than that of the B-Spline model used as a reference. The accuracy of the FTPH model regardless of the initial MH provided an opportunity to have a reference to check the accuracy of other pharmacokinetic models.

Machine Learning-Based Radiomics Nomogram for Detecting Extramural Venous Invasion in Rectal Cancer

Objective

To establish and validate a radiomics nomogram based on the features of the primary tumor for predicting preoperative pathological extramural venous invasion (EMVI) in rectal cancer using machine learning.

Methods

The clinical and imaging data of 281 patients with primary rectal cancer from April 2012 to May 2018 were retrospectively analyzed. All the patients were divided into a training set (n = 198) and a test set (n = 83) respectively. The radiomics features of the primary tumor were extracted from the enhanced computed tomography (CT), the T2-weighted imaging (T2WI) and the gadolinium contrast-enhanced T1-weighted imaging (CE-TIWI) of each patient. One optimal radiomics signature extracted from each modal image was generated by receiver operating characteristic (ROC) curve analysis after dimensionality reduction. Three kinds of models were constructed based on training set, including the clinical model (the optimal radiomics signature combining with the clinical features), the magnetic resonance imaging model (the optimal radiomics signature combining with the mrEMVI status) and the integrated model (the optimal radiomics signature combining with both the clinical features and the mrEMVI status). Finally, the optimal model was selected to create a radiomics nomogram. The performance of the nomogram to evaluate clinical efficacy was verified by ROC curves and decision curve analysis curves.

Results

The radiomics signature constructed based on T2WI showed the best performance, with an AUC value of 0.717, a sensitivity of 0.742 and a specificity of 0.621. The radiomics nomogram had the highest prediction efficiency, of which the AUC was 0.863, the sensitivity was 0.774 and the specificity was 0.801.

Conclusion

The radiomics nomogram had the highest efficiency in predicting EMVI. This may help patients choose the best treatment strategy and may strengthen personalized treatment methods to further optimize the treatment effect.

Imaging-based characterization of convective tissue properties

Convective transport is an important phenomenon for nanomedicine delivery. We present an imaging-based approach to recover tissue properties that are significant in the accumulation of nanoparticles delivered via systemic methods. The classical pharmacokinetic analysis develops governing equations for the particle transport from a first principle mass balance. Fundamentally, the governing equations for compartmental mass balance represent a spatially invariant mass transport between compartments and do not capture spatially variant convection phenomena. Further, the parameters recovered from this approach do not necessarily have direct meaning with respect to the governing equations for convective transport. In our approach, a framework is presented for directly measuring permeability in the sense of Darcy flow through porous tissue. Measurements from our approach are compared to an extended Tofts model as a control. We demonstrate that a pixel-wise iterative clustering algorithm may be applied to reduce the parameter space of the measurements. We show that measurements obtained from our approach are correlated with measurements obtained from the extended Tofts model control. These correlations demonstrate that the proposed approach contains similar information to an established compartmental model and may be useful in providing an alternative theoretical framework for parameterizing mathematical models for treatment planning and diagnostic studies involving nanomedicine where convection dominated effects are important.

Analysis Protocol for Dynamic Contrast Enhanced (DCE) MRI of Renal Perfusion and Filtration

Here we present an analysis protocol for dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) data of the kidneys. It covers comprehensive steps to facilitate signal to contrast agent concentration mapping via T1 mapping and the calculation of renal perfusion and filtration parametric maps using model-free approaches, model free analysis using deconvolution, the Toft's model and a Bayesian approach.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This analysis protocol chapter is complemented by two separate chapters describing the basic concept and experimental procedure.

Dynamic Contrast Enhancement (DCE) MRI-Derived Renal Perfusion and Filtration: Basic Concepts

Dynamic contrast-enhanced (DCE) MRI monitors the transit of contrast agents, typically gadolinium chelates, through the intrarenal regions, the renal cortex, the medulla, and the collecting system. In this way, DCE-MRI reveals the renal uptake and excretion of the contrast agent. An optimal DCE-MRI acquisition protocol involves finding a good compromise between whole-kidney coverage (i.e., 3D imaging), spatial and temporal resolution, and contrast resolution. By analyzing the enhancement of the renal tissues as a function of time, one can determine indirect measures of clinically important single-kidney parameters as the renal blood flow, glomerular filtration rate, and intrarenal blood volumes. Gadolinium-containing contrast agents may be nephrotoxic in patients suffering from severe renal dysfunction, but otherwise DCE-MRI is clearly useful for diagnosis of renal functions and for assessing treatment response and posttransplant rejection.Here we introduce the concept of renal DCE-MRI, describe the existing methods, and provide an overview of preclinical DCE-MRI applications to illustrate the utility of this technique to measure renal perfusion and glomerular filtration rate in animal models.This publication is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This introduction is complemented by two separate publications describing the experimental procedure and data analysis.

Exploring the Inter-voxel Information in Pharmacokinetic Maps for Cervical Carcinoma Prediction

Physiological parameters can be estimated from dynamic contrast enhanced magnetic resonance imaging (DCEMRI) data using pharmacokinetic models. This work evaluates the performance of various pharmacokinetic models through a retrospective study on cervix cancer, including two generalized kinetic models and three 2-compartment exchange models (2CXMs). In the current clinical practice, region of interest (ROI) is treated as a whole and the models are assessed by their top pharmacokinetic parameters. We explore the intervoxel relationship in the pharmacokinetic parameter maps and demonstrate that, for those insignificant parameters, texture descriptors can largely improve their discriminative power. Multi-parametric classifiers are developed to fuse the information carried by physiological parameters and the descriptors. Assessed merely by the top parameter, the DP (distributed parameter) model is the best one with an area under the ROC (receiver operating characteristic) curve (AUC) of 0.80; by combining multiple pharmacokinetic parameters, the ExTofts model is the winner with an AUC of 0.837. Finally, the classifier of the AATH (adiabatic approximation to the tissue homogeneity) model build on combined features achieves an AUC of 0.92.Clinical Relevance - Using data from 36 cervical cancer patients and 17 normal subjects, this work quantitatively compared the various pharmacokinetic models and provided recommendations for model selection in cervical cancer diagnosis.

Machine Learning-Based Radiomics Nomogram Using Magnetic Resonance Images for Prediction of Neoadjuvant Chemotherapy Efficacy in Breast Cancer Patients

Purpose: The construction and validation of a radiomics nomogram based on machine learning using magnetic resonance image (MRI) for predicting the efficacy of neoadjuvant chemotherapy (NACT) in patients with breast cancer (BCa).

Methods: This retrospective investigation consisted of 158 patients who were diagnosed with BCa and underwent MRI before NACT, of which 33 patients experienced pathological complete response (pCR) by the postoperative pathological examination. The patients with BCa were divided into the training set (n = 110) and test set (n = 48) randomly. The features were selected by the maximum relevance minimum redundancy (mRMR) and absolute shrinkage and selection operator (LASSO) algorithm in the training set. In return, the radiomics signature was established using machine learning. The predictive score of each patient was calculated using the radiomics signature formula. Finally, the predictive scores and clinical factors were used to perform the multivariate logistic regression and construct the nomogram. Receiver operating characteristics (ROC) analyses were used to assess and validate the diagnostic accuracy of the nomogram in the test set. Lastly, the usefulness of the nomogram was confirmed via decision curve analysis (DCA).

Results: The radiomics signature was well-discriminated in the training set [AUC 0.835, specificity 71.32%, and sensitivity 82.61%], and test set (AUC 0.834, specificity 73.21%, and sensitivity 80%). Containing the radiomics signature and hormone status, the radiomics nomogram showed good calibration and discrimination in the training set [AUC 0.888, specificity 79.31%, and sensitivity 86.96%] and test set (AUC 0.879, specificity 82.19%, and sensitivity 83.57%). The decision curve indicated the clinical usefulness of our nomogram.

Conclusion: Our radiomics nomogram showed good discrimination in patients with BCa who experience pCR after NACT. The model may aid physicians in predicting how specific patients may respond to BCa treatments in the future.

18F-FDG PET and DCE kinetic modeling and their correlations in primary NSCLC: first voxel-wise correlative analysis of human simultaneous [18F]FDG PET-MRI data

OBJECTIVES

To decipher the correlations between PET and DCE kinetic parameters in non-small-cell lung cancer (NSCLC), by using voxel-wise analysis of dynamic simultaneous [18F]FDG PET-MRI.

MATERIAL AND METHODS

Fourteen treatment-naïve patients with biopsy-proven NSCLC prospectively underwent a 1-h dynamic [18F]FDG thoracic PET-MRI scan including DCE. The PET and DCE data were normalized to their corresponding T1-weighted MR morphological space, and tumors were masked semi-automatically. Voxel-wise parametric maps of PET and DCE kinetic parameters were computed by fitting the dynamic PET and DCE tumor data to the Sokoloff and Extended Tofts models respectively, by using in-house developed procedures. Curve-fitting errors were assessed by computing the relative root mean square error (rRMSE) of the estimated PET and DCE signals at the voxel level. For each tumor, Spearman correlation coefficients (rs) between all the pairs of PET and DCE kinetic parameters were estimated on a voxel-wise basis, along with their respective bootstrapped 95% confidence intervals (n = 1000 iterations).

RESULTS

Curve-fitting metrics provided fit errors under 20% for almost 90% of the PET voxels (median rRMSE = 10.3, interquartile ranges IQR = 8.1; 14.3), whereas 73.3% of the DCE voxels showed fit errors under 45% (median rRMSE = 31.8%, IQR = 22.4; 46.6). The PET-PET, DCE-DCE, and PET-DCE voxel-wise correlations varied according to individual tumor behaviors. Beyond this wide variability, the PET-PET and DCE-DCE correlations were mainly high (absolute rs values > 0.7), whereas the PET-DCE correlations were mainly low to moderate (absolute rs values < 0.7). Half the tumors showed a hypometabolism with low perfused/vascularized profile, a hallmark of hypoxia, and tumor aggressiveness.

CONCLUSION

A dynamic "one-stop shop" procedure applied to NSCLC is technically feasible in clinical practice. PET and DCE kinetic parameters assessed simultaneously are not highly correlated in NSCLC, and these correlations showed a wide variability among tumors and patients. These results tend to suggest that PET and DCE kinetic parameters might provide complementary information. In the future, this might make PET-MRI a unique tool to characterize the individual tumor biological behavior in NSCLC.

Efficient DCE-MRI Parameter and Uncertainty Estimation Using a Neural Network

Quantitative DCE-MRI provides voxel-wise estimates of tracer-kinetic parameters that are valuable in the assessment of health and disease. These maps suffer from many known sources of variability. This variability is expensive to compute using current methods, and is typically not reported. Here, we demonstrate a novel approach for simultaneous estimation of tracer-kinetic parameters and their uncertainty due to intrinsic characteristics of the tracer-kinetic model, with very low computation time. We train and use a neural network to estimate the approximate joint posterior distribution of tracer-kinetic parameters. Uncertainties are estimated for each voxel and are specific to the patient, exam, and lesion. We demonstrate the methods' ability to produce accurate tracer-kinetic maps. We compare predicted parameter ranges with uncertainties introduced by noise and by differences in post-processing in a digital reference object. The predicted parameter ranges correlate well with tracer-kinetic parameter ranges observed across different noise realizations and regression algorithms. We also demonstrate the value of this approach to differentiate significant from insignificant changes in brain tumor pharmacokinetics over time. This is achieved by enforcing consistency in resolving model singularities in the applied tracer-kinetic model.

Clinical Implementation of a Free-Breathing, Motion-Robust Dynamic Contrast-Enhanced MRI Protocol to Evaluate Pleural Tumors

OBJECTIVE. The purpose of this study was to develop a motion insensitive clinical dynamic contrast-enhanced MRI (DCE-MRI) protocol to assess the response of pleural tumors in clinical trials. MATERIALS AND METHODS. Thirty-two patients with pleura-based lesions were administered contrast material and imaged with gradient-recalled echo DCE-MRI sequence variants: either a traditional cartesian k-space acquisition (FLASH), a time-resolved imaging with stochastic trajectories acquisition (TWIST), or a radial stack-of-stars acquisition (radial) sequence in addition to other standard-of-care imaging sequences. Each image acquisition's sensitivity to motion was evaluated by comparing the motion of the thoracic border in 3D throughout the acquisition. One-way ANOVA was used to compare the image quality between different acquisitions. The 95% CIs were calculated for mean thoracic border displacement. The effects of motion on kinetic parameter estimation were explored with simulations according to clinically acquired data. RESULTS. Radial was the most motion-robust sequence with subvoxel mean displacement in the superior-inferior direction (0.4 ± 1.2 [SD] mm). FLASH showed intermediate displacement (4.6 ± 2.0 mm), whereas TWIST was most sensitive to motion (6.4 ± 3.4 mm). Signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) of the images acquired with the radial sequence were on par or better than the FLASH and TWIST sequences when reconstructed with an improved density compensation algorithm. Simulations showed that motion on scans showing pleural-based lesions can lead to markedly inaccurate kinetic parameter estimation and inappropriate kinetic model convergence within a nested model analysis. CONCLUSION. A practical radial k-space trajectory sequence that provides motion-insensitive pharmacokinetic parameters was incorporated as part of the DCE-MRI protocol of pleural tumors. Validation and usefulness in clinical trials assessing response to therapy is needed.

Toward noninvasive quantification of adipose tissue oxygenation with MRI

BACKGROUND

Molecular oxygen (O₂) plays a key role in normal and pathological adipose tissue function, yet technologies to measure its role in adipose tissue function are limited. O₂ is paramagnetic and, in principle, directly influences the magnetic resonance (MR) ¹H longitudinal relaxation rate constant of lipids, R₁; thus, we hypothesize that MR imaging (MRI) can directly measure adipose O₂ via a simple measure of R₁.

METHODS

R₁ was measured in a 4.7T preclinical MRI system at discrete oxygen partial pressure (pO₂) levels. These measures were made in vitro in an idealized system and in vivo in subcutaneous and visceral white adipose of rodents. pO₂ was determined using an invasive fiber-optic oxygen monitor. From the MRI and fiber optic data we determined the "relaxivity" of O₂ in lipid, a critical parameter in converting the MRI-based R₁ measurement into pO₂. We used breathing gas challenge to estimate the changes in lipid pO₂ (ΔpO₂).

RESULTS

The relaxivity of O₂ in lipid was determined to be 1.7·10^-3 ± 4·10^-4 mmHg^-1s^-1 at 4.7T and 37 °C, and was consistent between in vitro and in vivo adipose tissue. There was a strong, significant correlation between MRI- and gold standard OxyLite-based measurements of lipid ΔpO₂ for in vivo visceral and subcutaneous fat depots in rodents.

CONCLUSION

This study lays the foundation for a direct, noninvasive measure of adipose pO₂ using MRI and will allow for noninvasive measurement of O₂ flux in adipose tissue. The proposed approach would be of particular importance in the interrogation of the pathogenesis of type 2 diabetes, where it has been suggested that adipose tissue hypoxia is an independent driver of insulin resistance pathway.

DCE-MRI pharmacokinetic parameter maps for cervical carcinoma prediction

Pharmacokinetic parameters estimated from dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) time course data enable the physio-biological interpretation of tissue angiogenesis. This study aims to develop machine learning approaches for cervical carcinoma prediction based on pharmacokinetic parameters. The performance of individual parameters was assessed in terms of their efficacy in differentiating cancerous tissue from normal cervix tissue. The effect of combining parameters was evaluated using the following two approaches: the first approach was based on support vector machines (SVMs) to combine the parameters from one pharmacokinetic model or across several models; the second approach was based on a novel method called APITL (artificial pharmacokinetic images for transfer learning), which was designed to fully utilize the comprehensive pharmacokinetic information acquired from DCE-MRI data. A "winner-takes-all" strategy was employed to consolidate the slice-wise prediction into subject-wise prediction. Experiments were carried out with a dataset comprising 36 patients with cervical cancer and 17 healthy subjects. The results demonstrated that parameter V_e, representing volume fraction of the extracellular extravascular space (EES), attained high discriminative power regardless of the pharmacokinetic model used for estimation. An approximately 10% improvement in the accuracy was achieved with the SVM approach. The APITL method further outperformed SVM and attained a subject-wise prediction accuracy of 94.3%. Our experiment demonstrated that APITL could predict cervical carcinoma with high accuracy and had potential in clinical applications.

Dynamic contrast-enhanced MRI model selection for predicting tumor aggressiveness in papillary thyroid cancers

The purpose of this study was to identify the optimal tracer kinetic model from T1 -weighted dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) data and evaluate whether parameters estimated from the optimal model predict tumor aggressiveness determined from histopathology in patients with papillary thyroid carcinoma (PTC) prior to surgery. In this prospective study, 18 PTC patients underwent pretreatment DCE-MRI on a 3 T MR scanner prior to thyroidectomy. This study was approved by the institutional review board and informed consent was obtained from all patients. The two-compartment exchange model, compartmental tissue uptake model, extended Tofts model (ETM) and standard Tofts model were compared on a voxel-wise basis to determine the optimal model using the corrected Akaike information criterion (AICc) for PTC. The optimal model is the one with the lowest AICc. Statistical analysis included paired and unpaired t-tests and a one-way analysis of variance. Bonferroni correction was applied for multiple comparisons. Receiver operating characteristic (ROC) curves were generated from the optimal model parameters to differentiate PTC with and without aggressive features, and AUCs were compared. ETM performed best with the lowest AICc and the highest Akaike weight (0.44) among the four models. ETM was preferred in 44% of all 3419 voxels. The ETM estimates of Ktrans in PTCs with the aggressive feature extrathyroidal extension (ETE) were significantly higher than those without ETE (0.78 ± 0.29 vs. 0.34 ± 0.18 min-1 , P = 0.005). From ROC analysis, cut-off values of Ktrans , ve and vp , which discriminated between PTCs with and without ETE, were determined at 0.45 min-1 , 0.28 and 0.014 respectively. The sensitivities and specificities were 86 and 82% (Ktrans ), 71 and 82% (ve ), and 86 and 55% (vp ), respectively. Their respective AUCs were 0.90, 0.71 and 0.71. We conclude that ETM Ktrans has shown potential to classify tumors with and without aggressive ETE in patients with PTC.

Bayesian pharmacokinetic modeling of dynamic contrast-enhanced magnetic resonance imaging: validation and application

Tracer-kinetic analysis of dynamic contrast-enhanced magnetic resonance imaging data is commonly performed with the well-known Tofts model and nonlinear least squares (NLLS) regression. This approach yields point estimates of model parameters, uncertainty of these estimates can be assessed e.g. by an additional bootstrapping analysis. Here, we present a Bayesian probabilistic modeling approach for tracer-kinetic analysis with a Tofts model, which yields posterior probability distributions of perfusion parameters and therefore promises a robust and information-enriched alternative based on a framework of probability distributions. In this manuscript, we use the quantitative imaging biomarkers alliance (QIBA) Tofts phantom to evaluate the Bayesian tofts model (BTM) against a bootstrapped NLLS approach. Furthermore, we demonstrate how Bayesian posterior probability distributions can be employed to assess treatment response in a breast cancer DCE-MRI dataset using Cohen's d. Accuracy and precision of the BTM posterior distributions were validated and found to be in good agreement with the NLLS approaches, and assessment of therapy response with respect to uncertainty in parameter estimates was found to be excellent. In conclusion, the Bayesian modeling approach provides an elegant means to determine uncertainty via posterior distributions within a single step and provides honest information about changes in parameter estimates.

Reliability of dynamic contrast-enhanced magnetic resonance imaging data in primary brain tumours: a comparison of Tofts and shutter speed models

Purpose

The purpose of this study is to investigate the robustness of pharmacokinetic modelling of DCE-MRI brain tumour data and to ascertain reliable perfusion parameters through a model selection process and a stability test.

Methods

DCE-MRI data of 14 patients with primary brain tumours were analysed using the Tofts model (TM), the extended Tofts model (ETM), the shutter speed model (SSM) and the extended shutter speed model (ESSM). A no-effect model (NEM) was implemented to assess overfitting of data by the other models. For each lesion, the Akaike Information Criteria (AIC) was used to build a 3D model selection map. The variability of each pharmacokinetic parameter extracted from this map was assessed with a noise propagation procedure, resulting in voxel-wise distributions of the coefficient of variation (CV).

Results

The model selection map over all patients showed NEM had the best fit in 35.5% of voxels, followed by ETM (32%), TM (28.2%), SSM (4.3%) and ESSM (< 0.1%). In analysing the reliability of K^trans, when considering regions with a CV < 20%, ≈ 25% of voxels were found to be stable across all patients. The remaining 75% of voxels were considered unreliable.

Conclusions

The majority of studies quantifying DCE-MRI data in brain tumours only consider a single model and whole tumour statistics for the output parameters. Appropriate model selection, considering tissue biology and its effects on blood brain barrier permeability and exchange conditions, together with an analysis on the reliability and stability of the calculated parameters, is critical in processing robust brain tumour DCE-MRI data.

Improved repeatability of dynamic contrast-enhanced MRI using the complex MRI signal to derive arterial input functions: a test-retest study in prostate cancer patients

Purpose

The arterial input function (AIF) is a major source of uncertainty in tracer kinetic (TK) analysis of dynamic contrast‐enhanced (DCE)‐MRI data. The aim of this study was to investigate the repeatability of AIFs extracted from the complex signal and of the resulting TK parameters in prostate cancer patients.

Methods

Twenty‐two patients with biopsy‐proven prostate cancer underwent a 3T MRI exam twice. DCE‐MRI data were acquired with a 3D spoiled gradient echo sequence. AIFs were extracted from the magnitude of the signal (AIF_MAGN), phase (AIF_PHASE), and complex signal (AIF_COMPLEX). The Tofts model was applied to extract K^trans, k_ep and v_e. Repeatability of AIF curve characteristics and TK parameters was assessed with the within‐subject coefficient of variation (wCV).

Results

The wCV for peak height and full width at half maximum for AIF_COMPLEX (7% and 8%) indicated an improved repeatability compared to AIF_MAGN (12% and 12%) and AIF_PHASE (12% and 7%). This translated in lower wCV values for K^trans (11%) with AIF_COMPLEX in comparison to AIF_MAGN (24%) and AIF_PHASE (15%). For k_ep, the wCV was 16% with AIF_MAGN, 13% with AIF_PHASE, and 13% with AIF_COMPLEX.

Conclusion

Repeatability of AIF_PHASE and AIF_COMPLEX is higher than for AIF_MAGN, resulting in a better repeatability of TK parameters. Thus, use of either AIF_PHASE or AIF_COMPLEX improves the robustness of quantitative analysis of DCE‐MRI in prostate cancer.

Robust estimation of hemo-dynamic parameters in traditional DCE-MRI models

PURPOSE

In dynamic contrast enhanced (DCE) MRI, separation of signal contributions from perfusion and leakage requires robust estimation of parameters in a pharmacokinetic model. We present and quantify the performance of a method to compute tissue hemodynamic parameters from DCE data using established pharmacokinetic models.

METHODS

We propose a Bayesian scheme to obtain perfusion metrics from DCE MRI data. Initial performance is assessed through digital phantoms of the extended Tofts model (ETM) and the two-compartment exchange model (2CXM), comparing the Bayesian scheme to the standard Levenberg-Marquardt (LM) algorithm. Digital phantoms are also invoked to identify limitations in the pharmacokinetic models related to measurement conditions. Using computed maps of the extra vascular volume (ve) from 19 glioma patients, we analyze differences in the number of un-physiological high-intensity ve values for both ETM and 2CXM, using a one-tailed paired t-test assuming un-equal variance.

RESULTS

The Bayesian parameter estimation scheme demonstrated superior performance over the LM technique in the digital phantom simulations. In addition, we identified limitations in parameter reliability in relation to scan duration for the 2CXM. DCE data for glioma and cervical cancer patients was analyzed with both algorithms and demonstrated improvement in image readability for the Bayesian method. The Bayesian method demonstrated significantly fewer non-physiological high-intensity ve values for the ETM (p<0.0001) and the 2CXM (p<0.0001).

CONCLUSION

We have demonstrated substantial improvement of the perceptive quality of pharmacokinetic parameters from advanced compartment models using the Bayesian parameter estimation scheme as compared to the LM technique.

Impact of contrast agent injection duration on dynamic contrast-enhanced MRI quantification in prostate cancer

The volume transfer constant K^trans , which describes the leakage of contrast agent (CA) from vasculature into tissue, is the most commonly reported quantitative parameter for dynamic contrast-enhanced (DCE-) MRI. However, the variation in reported K^trans values between studies from different institutes is large. One of the primary sources of uncertainty is quantification of the arterial input function (AIF). The aim of this study is to determine the influence of the CA injection duration on the AIF and tracer kinetic analysis (TKA) parameters (i.e. K^trans , k_ep and v_e ). Thirty-one patients with prostate cancer received two DCE-MRI examinations with an injection duration of 5 s in the first examination and a prolonged injection duration in the second examination, varying between 7.5 s and 30 s. The DCE examination was carried out on a 3.0 T MRI scanner using a transversal T₁ -weighted 3D spoiled gradient echo sequence (300 s duration, dynamic scan time of 2.5 s). Data of 29 of the 31 were further analysed. AIFs were determined from the phase signal in the left and right femoral arteries. K^trans , k_ep and v_e were estimated with the standard Tofts model for regions of healthy peripheral zone and tumour tissue. We observed a significantly smaller peak height and increased width in the AIF for injection durations of 15 s and longer. However, we did not find significant differences in K^trans , k_ep or v_e for the studied injection durations. The study demonstrates that the TKA parameters K^trans , k_ep and v_e , measured in the prostate, do not show a significant change as a function of injection duration.

Pancreatic imaging: Current status of clinical practices and small animal studies

Different causative factors acting on the pancreas can result in diseases such as pancreatitis, diabetes and pancreatic tumors. The high incidence and mortality of pancreatic diseases have placed diagnostic imaging in a crucial position in daily clinical practice. In this mini-review article different pancreatic imaging techniques are discussed, from the standard clinical imaging modalities and state of the art clinical magnetic resonance imaging techniques to current situations in pre-clinical pancreatic imaging studies. In particular, the challenges of pre-clinical rodent pancreatic imaging are addressed, with both the image acquisition techniques and the post-processing methods for rodent pancreatic imaging elaborated.

Reproducibility and relative stability in magnetic resonance imaging indices of tumor vascular physiology over a period of 24h in a rat 9L gliosarcoma model

PURPOSE

The objective was to study temporal changes in tumor vascular physiological indices in a period of 24h in a 9L gliosarcoma rat model.

METHODS

Fischer-344 rats (N=14) were orthotopically implanted with 9L cells. At 2weeks post-implantation, they were imaged twice in a 24h interval using dynamic contrast enhanced magnetic resonance imaging (DCE-MRI). Data-driven model-selection-based analysis was used to segment tumor regions with varying vascular permeability characteristics. The region with the maximum number of estimable parameters of vascular kinetics was chosen for comparison across the two time points. It provided estimates of three parameters for an MR contrast agent (MRCA): i) plasma volume (v_p), ii) forward volumetric transfer constant (K^trans) and interstitial volume fraction (v_e, ratio of K^trans to reverse transfer constant, k_ep). In addition, MRCA extracellular distribution volume (V_D) was estimated in the tumor and its borders, along with tumor blood flow (TBF) and peritumoral MRCA flux. Descriptors of parametric distributions were compared between the two times. Tumor extent was examined by hematoxylin and eosin (H&E) staining. Picrosirus red staining of secreted collagen was performed as an additional index for 9L cells.

RESULTS

Test-retest differences between population summaries for any parameter were not significant (paired t and Wilcoxon signed rank tests). Bland-Altman plots showed no apparent trends between the differences and averages of the test-retest measures for all indices. The intraclass correlation coefficients showed moderate to almost perfect reproducibility for all of the parameters, except v_p. H&E staining showed tumor infiltration in parenchyma, perivascular space and white matter tracts. Collagen staining was observed along the outer edges of main tumor mass.

CONCLUSION

The data suggest the relative stability of these MR indices of tumor microenvironment over a 24h duration in this gliosarcoma model.

Modeling Dynamic Contrast-Enhanced MRI Data with a Constrained Local AIF

PURPOSE

This study aims to develop a constrained local arterial input function (cL-AIF) to improve quantitative analysis of dynamic contrast-enhanced (DCE)-magnetic resonance imaging (MRI) data by accounting for the contrast-agent bolus amplitude error in the voxel-specific AIF.

PROCEDURES

Bayesian probability theory-based parameter estimation and model selection were used to compare tracer kinetic modeling employing either the measured remote-AIF (R-AIF, i.e., the traditional approach) or an inferred cL-AIF against both in silico DCE-MRI data and clinical, cervical cancer DCE-MRI data.

RESULTS

When the data model included the cL-AIF, tracer kinetic parameters were correctly estimated from in silico data under contrast-to-noise conditions typical of clinical DCE-MRI experiments. Considering the clinical cervical cancer data, Bayesian model selection was performed for all tumor voxels of the 16 patients (35,602 voxels in total). Among those voxels, a tracer kinetic model that employed the voxel-specific cL-AIF was preferred (i.e., had a higher posterior probability) in 80 % of the voxels compared to the direct use of a single R-AIF. Maps of spatial variation in voxel-specific AIF bolus amplitude and arrival time for heterogeneous tissues, such as cervical cancer, are accessible with the cL-AIF approach.

CONCLUSIONS

The cL-AIF method, which estimates unique local-AIF amplitude and arrival time for each voxel within the tissue of interest, provides better modeling of DCE-MRI data than the use of a single, measured R-AIF. The Bayesian-based data analysis described herein affords estimates of uncertainties for each model parameter, via posterior probability density functions, and voxel-wise comparison across methods/models, via model selection in data modeling.

Estimating breast tumor blood flow during neoadjuvant chemotherapy using interleaved high temporal and high spatial resolution MRI

PURPOSE

To evaluate an interleaved MRI sampling strategy that acquires both high temporal resolution (HTR) dynamic contrast-enhanced (DCE) data for quantifying breast tumor blood flow (TBF) and high spatial resolution (HSR) DCE data for clinical reporting, following a single standard injection of contrast agent.

METHODS

A simulation study was used to evaluate the performance of the interleaved technique under different conditions. In a prospective clinical study, 18 patients with primary breast cancer, who were due to undergo neoadjuvant chemotherapy (NACT), were examined using interleaved HTR and HSR DCE-MRI at 1.5 Tesla. Tumor regions of interest were analyzed with a two-compartment tracer kinetic model. Paired parameters (n = 10) from the data acquired before and post-cycle 2 of NACT were compared using the nonparametric Wilcoxon signed-rank test.

RESULTS

Simulations demonstrated that TBF was reliably estimated using the proposed strategy. The region of interest analysis revealed significant changes in TBF (0.81-0.43 mL/min/mL; P = 0.002) following two cycles of NACT. The HSR data were reported in the normal way and enabled the assessment of tumor volume, which decreased by 53% following NACT (P = 0.065).

CONCLUSIONS

TBF can be measured reliably using the proposed strategy without compromising a standard clinical protocol. Furthermore, in our feasibility study, TBF decreased significantly following NACT, whereas capillary permeability surface-area product did not. Magn Reson Med 79:317-326, 2018. © 2017 International Society for Magnetic Resonance in Medicine.