Simulating the effect of input errors on the accuracy of Tofts' pharmacokinetic model parameters

Unified Bayesian network for uncertainty quantification of physiological parameters in dynamic contrast enhanced (DCE) MRI of the liver

Objective. Physiological parameter estimation is affected by intrinsic ambiguity in the data such as noise and model inaccuracies. The aim of this work is to provide a deep learning framework for accurate parameter and uncertainty estimates for DCE-MRI in the liver.Approach. Concentration time curves are simulated to train a Bayesian neural network (BNN). Training of the BNN involves minimization of a loss function that jointly minimizes the aleatoric and epistemic uncertainties. Uncertainty estimation is evaluated for different noise levels and for different out of distribution (OD) cases, i.e. where the data during inference differs strongly to the data during training. The accuracy of parameter estimates are compared to a nonlinear least squares (NLLS) fitting in numerical simulations andin vivodata of a patient suffering from hepatic tumor lesions.Main results. BNN achieved lower root-mean-squared-errors (RMSE) than the NLLS for the simulated data. RMSE of BNN was on overage of all noise levels lower by 33% ± 1.9% forktrans, 22% ± 6% forveand 89% ± 5% forvpthan the NLLS. The aleatoric uncertainties of the parameters increased with increasing noise level, whereas the epistemic uncertainty increased when a BNN was evaluated with OD data. For thein vivodata, more robust parameter estimations were obtained by the BNN than the NLLS fit. In addition, the differences between estimated parameters for healthy and tumor regions-of-interest were significant (p< 0.0001).Significance. The proposed framework allowed for accurate parameter estimates for quantitative DCE-MRI. In addition, the BNN provided uncertainty estimates which highlighted cases of high noise and in which the training data did not match the data during inference. This is important for clinical application because it would indicate cases in which the trained model is inadequate and additional training with an adapted training data set is required.

The arterial input function: Spatial dependence within the imaging volume and its influence on 3D quantitative dynamic contrast-enhanced MRI for head and neck cancer

PURPOSE

To evaluate the dependence of the arterial input function (AIF) on the imaging z-axis and its effect on 3D DCE MRI pharmacokinetic parameters as mediated by the SPGR signal equation and Extended Tofts-Kermode model.

THEORY

For SPGR-based 3D DCE MRI acquisition of the head and neck, inflow effects within vessels violate the assumptions underlying the SPGR signal model. Errors in the SPGR-based AIF estimate propagate through the Extended Tofts-Kermode model to affect the output pharmacokinetic parameters.

MATERIALS AND METHODS

3D DCE-MRI data were acquired for six newly diagnosed HNC patients in a prospective single arm cohort study. AIF were selected within the carotid arteries at each z-axis location. A region of interest (ROI) was placed in normal paravertebral muscle and the Extended Tofts-Kermode model solved for each pixel within the ROI for each AIF. Results were compared to those obtained with a published population average AIF.

RESULTS

Due to inflow effect, the AIF showed extreme variation in their temporal shapes. K^trans was most sensitive to the initial bolus concentration and showed more variation over the muscle ROI with AIF taken from the upstream portion of the carotid. k_ep was less sensitive to the peak bolus concentration and showed less variation for AIF taken from the upstream portion of the carotid.

CONCLUSION

Inflow effects may introduce an unknown bias to SPGR-based 3D DCE pharmacokinetic parameters. Variation in the computed parameters depends on the selected AIF location. In the context of high flow, measurements may be limited to relative rather than absolute quantitative parameters.

Characterizing the Bone Marrow Environment in Advanced-Stage Myelofibrosis during Ruxolitinib Treatment Using PET/CT and MRI: A Pilot Study

Current diagnostic criteria for myelofibrosis are largely based on bone marrow (BM) biopsy results. However, these have several limitations, including sampling errors. Explorative studies have indicated that imaging might form an alternative for the evaluation of disease activity, but the heterogeneity in BM abnormalities complicates the choice for the optimal technique. In our prospective diagnostic pilot study, we aimed to visualize all BM abnormalities in myelofibrosis before and during ruxolitinib treatment using both PET/CT and MRI. A random sample of patients was scheduled for examinations at baseline and after 6 and 18 months of treatment, including clinical and laboratory examinations, BM biopsies, MRI (T1-weighted, Dixon, dynamic contrast-enhanced (DCE)) and PET/CT ([¹⁵O]water, [¹⁸F]NaF)). At baseline, all patients showed low BM fat content (indicated by T1-weighted MRI and Dixon), increased BM blood flow (as measured by [¹⁵O]water PET/CT), and increased osteoblastic activity (reflected by increased skeletal [¹⁸F]NaF uptake). One patient died after the baseline evaluation. In the others, BM fat content increased to various degrees during treatment. Normalization of BM blood flow (as reflected by [¹⁵O]water PET/CT and DCE-MRI) occurred in one patient, who also showed the fastest clinical response. Vertebral [¹⁸F]NaF uptake remained stable in all patients. In evaluable cases, histopathological parameters were not accurately reflected by imaging results. A case of sampling error was suspected. We conclude that imaging results can provide information on functional processes and disease distribution throughout the BM. Differences in early treatment responses were especially reflected by T1-weighted MRI. Limitations in the gold standard hampered the evaluation of diagnostic accuracy.

Influence of Different Measurement Methods of Arterial Input Function on Quantitative Dynamic Contrast-Enhanced MRI Parameters in Head and Neck Cancer

BACKGROUND

Head and neck cancer (HNC) is the sixth most prevalent cancer worldwide. Dynamic contrast-enhanced MRI (DCE-MRI) helps in diagnosis and prognosis. Quantitative DCE-MRI requires an arterial input function (AIF), which affects the values of pharmacokinetic parameters (PKP).

PURPOSE

To evaluate influence of four individual AIF measurement methods on quantitative DCE-MRI parameters values (K^trans , v_e , k_ep , and v_p ), for HNC and muscle.

STUDY TYPE

Prospective.

POPULATION

A total of 34 HNC patients (23 males, 11 females, age range 24-91) FIELD STRENGTH/SEQUENCE: A 3 T; 3D SPGR gradient echo sequence with partial saturation of inflowing spins.

ASSESSMENT

Four AIF methods were applied: automatic AIF (AIFa) with up to 50 voxels selected from the whole FOV, manual AIF (AIFm) with four voxels selected from the internal carotid artery, both conditions without (Mc-) or with (Mc+) motion correction. Comparison endpoints were peak AIF values, PKP values in tumor and muscle, and tumor/muscle PKP ratios.

STATISTICAL TESTS

Nonparametric Friedman test for multiple comparisons. Nonparametric Wilcoxon test, without and with Benjamini Hochberg correction, for pairwise comparison of AIF peak values and PKP values for tumor, muscle and tumor/muscle ratio, P value ≤ 0.05 was considered statistically significant.

RESULTS

Peak AIF values differed significantly for all AIF methods, with mean AIFmMc+ peaks being up to 66.4% higher than those for AIFaMc+. Almost all PKP values were significantly higher for AIFa in both, tumor and muscle, up to 76% for mean K^trans values. Motion correction effect was smaller. Considering tumor/muscle parameter ratios, most differences were not significant (0.068 ≤ Wilcoxon P value ≤ 0.8).

DATA CONCLUSION

We observed important differences in PKP values when using either AIFa or AIFm, consequently choice of a standardized AIF method is mandatory for DCE-MRI on HNC. From the study findings, AIFm and inflow compensation are recommended. The use of the tumor/muscle PKP ratio should be of interest for multicenter studies.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Maximum Entropy Technique and Regularization Functional for Determining the Pharmacokinetic Parameters in DCE-MRI

This paper aims to solve the arterial input function (AIF) determination in dynamic contrast-enhanced MRI (DCE-MRI), an important linear ill-posed inverse problem, using the maximum entropy technique (MET) and regularization functionals. In addition, estimating the pharmacokinetic parameters from a DCE-MR image investigations is an urgent need to obtain the precise information about the AIF-the concentration of the contrast agent on the left ventricular blood pool measured over time. For this reason, the main idea is to show how to find a unique solution of linear system of equations generally in the form of [Formula: see text] named an ill-conditioned linear system of equations after discretization of the integral equations, which appear in different tomographic image restoration and reconstruction issues. Here, a new algorithm is described to estimate an appropriate probability distribution function for AIF according to the MET and regularization functionals for the contrast agent concentration when applying Bayesian estimation approach to estimate two different pharmacokinetic parameters. Moreover, by using the proposed approach when analyzing simulated and real datasets of the breast tumors according to pharmacokinetic factors, it indicates that using Bayesian inference-that infer the uncertainties of the computed solutions, and specific knowledge of the noise and errors-combined with the regularization functional of the maximum entropy problem, improved the convergence behavior and led to more consistent morphological and functional statistics and results. Finally, in comparison to the proposed exponential distribution based on MET and Newton's method, or Weibull distribution via the MET and teaching-learning-based optimization (MET/TLBO) in the previous studies, the family of Gamma and Erlang distributions estimated by the new algorithm are more appropriate and robust AIFs.

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.

Modified Maximum Entropy Method and Estimating the AIF via DCE-MRI Data Analysis

Background: For the kinetic models used in contrast-based medical imaging, the assignment of the arterial input function named AIF is essential for the estimation of the physiological parameters of the tissue via solving an optimization problem. Objective: In the current study, we estimate the AIF relayed on the modified maximum entropy method. The effectiveness of several numerical methods to determine kinetic parameters and the AIF is evaluated—in situations where enough information about the AIF is not available. The purpose of this study is to identify an appropriate method for estimating this function. Materials and Methods: The modified algorithm is a mixture of the maximum entropy approach with an optimization method, named the teaching-learning method. In here, we applied this algorithm in a Bayesian framework to estimate the kinetic parameters when specifying the unique form of the AIF by the maximum entropy method. We assessed the proficiency of the proposed method for assigning the kinetic parameters in the dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), when determining AIF with some other parameter-estimation methods and a standard fixed AIF method. A previously analyzed dataset consisting of contrast agent concentrations in tissue and plasma was used. Results and Conclusions: We compared the accuracy of the results for the estimated parameters obtained from the MMEM with those of the empirical method, maximum likelihood method, moment matching (“method of moments”), the least-square method, the modified maximum likelihood approach, and our previous work. Since the current algorithm does not have the problem of starting point in the parameter estimation phase, it could find the best and nearest model to the empirical model of data, and therefore, the results indicated the Weibull distribution as an appropriate and robust AIF and also illustrated the power and effectiveness of the proposed method to estimate the kinetic parameters.

Pseudo Test-Retest Evaluation of Millimeter-Resolution Whole-Brain Dynamic Contrast-enhanced MRI in Patients with High-Grade Glioma

Background

Advances in sub-Nyquist–sampled dynamic contrast-enhanced (DCE) MRI enable monitoring of brain tumors with millimeter resolution and whole-brain coverage. Such undersampled quantitative methods need careful characterization regarding achievable test-retest reproducibility.

Purpose

To demonstrate a fully automated high-resolution whole-brain DCE MRI pipeline with 30-fold sparse undersampling and estimate its reproducibility on the basis of reference regions of stable tissue types during multiple posttreatment time points by using longitudinal clinical images of high-grade glioma.

Materials and Methods

Two methods for sub-Nyquist–sampled DCE MRI were extended with automatic estimation of vascular input functions. Continuously acquired three-dimensional k-space data with ramped-up flip angles were partitioned to yield high-resolution, whole-brain tracer kinetic parameter maps with matched precontrast-agent T1 and M₀ maps. Reproducibility was estimated in a retrospective study in participants with high-grade glioma, who underwent three consecutive standard-of-care examinations between December 2016 and April 2019. Coefficients of variation and reproducibility coefficients were reported for histogram statistics of the tracer kinetic parameters plasma volume fraction and volume transfer constant (K^trans) on five healthy tissue types.

Results

The images from 13 participants (mean age ± standard deviation, 61 years ± 10; nine women) with high-grade glioma were evaluated. In healthy tissues, the protocol achieved a coefficient of variation less than 57% for median K^trans, if K^trans was estimated consecutively. The maximum reproducibility coefficient for median K^trans was estimated to be at 0.06 min^–1 for large or low-enhancing tissues and to be as high as 0.48 min^–1 in smaller or strongly enhancing tissues.

Conclusion

A fully automated, sparsely sampled DCE MRI reconstruction with patient-specific vascular input function offered high spatial and temporal resolution and whole-brain coverage; in healthy tissues, the protocol estimated median volume transfer constant with maximum reproducibility coefficient of 0.06 min^–1 in large, low-enhancing tissue regions and maximum reproducibility coefficient of less than 0.48 min^–1 in smaller or more strongly enhancing tissue regions.

Published under a CC BY 4.0 license.

Online supplemental material is available for this article.

See also the editorial by Lenkinski in this issue.

Artificial intelligence and imaging biomarkers for prostate radiation therapy during and after treatment

Magnetic resonance imaging (MRI) is increasingly used in the management of prostate cancer (PCa). Quantitative MRI (qMRI) parameters, derived from multi-parametric MRI, provide indirect measures of tumour characteristics such as cellularity, angiogenesis and hypoxia. Using Artificial Intelligence (AI), relevant information and patterns can be efficiently identified in these complex data to develop quantitative imaging biomarkers (QIBs) of tumour function and biology. Such QIBs have already demonstrated potential in the diagnosis and staging of PCa. In this review, we explore the role of these QIBs in monitoring treatment response during and after PCa radiotherapy (RT). Recurrence of PCa after RT is not uncommon, and early detection prior to development of metastases provides an opportunity for salvage treatments with curative intent. However, the current method of monitoring treatment response using prostate-specific antigen levels lacks specificity. QIBs, derived from qMRI and developed using AI techniques, can be used to monitor biological changes post-RT providing the potential for accurate and early diagnosis of recurrent disease.

Quantitative MRI: Defining repeatability, reproducibility and accuracy for prostate cancer imaging biomarker development

INTRODUCTION

Quantitative MRI (qMRI) parameters have been increasingly used to develop predictive models to accurately monitor treatment response in prostate cancer after radiotherapy. To reliably detect changes in signal due to treatment response, predictive models require qMRI parameters with high repeatability and reproducibility. The purpose of this study was to measure qMRI parameter uncertainties in both commercial and in-house developed phantoms to guide the development of robust predictive models for monitoring treatment response.

MATERIALS AND METHODS

ADC, T1, and R2* values were acquired across three 3 T scanners with a prostate-specific qMRI protocol using the NIST/ISMRM system phantom, RSNA/NIST diffusion phantom, and an in-house phantom. A B1 field map was acquired to correct for flip angle inhomogeneity in T1 maps. All sequences were repeated in each scan to assess within-session repeatability. Weekly scans were acquired on one scanner for three months with the in-house phantom. Between-session repeatability was measured with test-retest scans 6-months apart on all scanners with all phantoms. Accuracy, defined as percentage deviation from reference value for ADC and T1, was evaluated using the system and diffusion phantoms. Repeatability and reproducibility coefficients of variation (%CV) were calculated for all qMRI parameters on all phantoms.

RESULTS

Overall, repeatability CV of ADC was <2.40%, reproducibility CV was <3.98%, and accuracy ranged between -8.0% to 2.7% across all scanners. Applying B1 correction on T1 measurements significantly improved the repeatability and reproducibility (p<0.05) but increased error in accuracy (p<0.001). Repeatability and reproducibility of R2* was <4.5% and <7.3% respectively in the system phantom across all scanners.

CONCLUSION

Repeatability, reproducibility, and accuracy in qMRI parameters from a prostate-specific protocol was estimated using both commercial and in-house phantoms. Results from this work will be used to identify robust qMRI parameters for use in the development of predictive models to longitudinally monitor treatment response for prostate cancer in current and future clinical trials.

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.

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.

MRI protocol optimization for quantitative DCE-MRI of the spine

PURPOSE

In this study we systematically investigated different Dynamic Contrast Enhancement (DCE)-MRI protocols in the spine, with the goal of finding an optimal protocol that provides data suitable for quantitative pharmacokinetic modelling (PKM).

MATERIALS AND METHODS

In 13 patients referred for MRI of the spine, DCE-MRI of the spine was performed with 2D and 3D MRI protocols on a 3T Philips Ingenuity MR system. A standard bolus of contrast agent (Dotarem - 0.2ml/kg body weight) was injected intravenously at a speed of 3ml/s. Different techniques for acceleration and motion compensation were tested: parallel imaging, partial-Fourier imaging and flow compensation. The quality of the DCE MRI images was scored on the basis of SNR, motion artefacts due to flow and respiration, signal enhancement, quality of the T₁ map and of the arterial input function, and quality of pharmacokinetic model fitting to the extended Tofts model.

RESULTS

Sagittal 3D sequences are to be preferred for PKM of the spine. Acceleration techniques were unsuccessful due to increased flow or motion artefacts. Motion compensating gradients failed to improve the DCE scans due to the longer echo time and the T₂* decay which becomes more dominant and leads to signal loss, especially in the aorta. The quality scoring revealed that the best method was a conventional 3D gradient-echo acquisition without any acceleration or motion compensation technique. The priority in the choice of sequence parameters should be given to reducing echo time and keeping the dynamic temporal resolution below 5s. Increasing the number of acquisition, when possible, helps towards reducing flow artefacts. In our setting we achieved this with a sagittal 3D slab with 5 slices with a thickness of 4.5mm and two acquisitions.

CONCLUSION

The proposed DCE protocol, encompassing the spine and the descending aorta, produces a realistic arterial input function and dynamic data suitable for PKM.

Assessment of DCE-MRI parameters for brain tumors through implementation of physiologically-based pharmacokinetic model approaches for Gd-DOTA

Dynamic-contrast enhanced magnetic resonance imaging (DCE-MRI) is used for detailed characterization of pathology of lesions sites, such as brain tumors, by quantitative analysis of tracer's data through the use of pharmacokinetic (PK) models. A key component for PK models in DCE-MRI is the estimation of the concentration-time profile of the tracer in a nearby vessel, referred as Arterial Input Function (AIF). The aim of this work was to assess through full body physiologically-based pharmacokinetic (PBPK) model approaches the PK profile of gadoteric acid (Gd-DOTA) and explore potential application for parameter estimation in DCE-MRI based on PBPK-derived AIFs. The PBPK simulations were generated through Simcyp(®) platform and the predicted PK parameters for Gd-DOTA were compared with available clinical data regarding healthy volunteers and renal impairment patients. The assessment of DCE-MRI parameters was implemented by utilizing similar virtual profiles based on gender, age and weight to clinical profiles of patients diagnosed with glioblastoma multiforme. The PBPK-derived AIFs were then used to compute DCE-MRI parameters through the Extended Tofts Model and compared with the corresponding ones derived from image-based AIF computation. The comparison involved: (i) image measured AIF of patients vs AIF of in silico profile, and, (ii) population average AIF vs in silico mean AIFs. The results indicate that PBPK-derived AIFs allowed the estimation of comparable imaging biomarkers with those calculated from typical DCE-MRI image analysis. The incorporation of PBPK models and potential utilization of in silico profiles to real patient data, can provide new perspectives in DCE-MRI parameter estimation and data analysis.

Dynamic Contrast-Enhanced Magnetic Resonance Imaging Using Pharmacokinetic Modeling: Initial Experience in Patients With Early Arthritis

OBJECTIVE

Analysis of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) using pharmacokinetic modeling (PKM) provides quantitative measures that mirror microvessel integrity and can be used as an objective marker of the level of synovial inflammation. The aim of this study was to investigate the PKM parameters K(trans) , kep , and ve in a prospective cohort of disease-modifying antirheumatic drug (DMARD)-naive patients with early arthritis, and to validate the results by assessing their correlation with the number of synovial endothelial cells (ECs).

METHODS

Forty-seven patients with early arthritis (arthritis duration <1 year, DMARD naive; comprising 14 patients with rheumatoid arthritis, 22 with unclassified arthritis, 6 with spondyloarthritis [SpA], and 5 with other arthritides) were included. At baseline, DCE-MRI was performed on an inflamed knee joint of each patient. These images were used to calculate the K(trans) (volume transfer constant between the plasma and extracellular extravascular space [EES]), the kep (transfer constant between the EES and plasma), and the ve (fractional volume of the EES). Second, markers of disease activity were collected. Finally, vascularity was evaluated by immunohistochemical analysis of synovial tissue samples obtained from the inflamed knee joints, using antibodies to detect von Willebrand factor (vWF), a marker of ECs.

RESULTS

The 3 PKM parameters differed significantly between diagnostic groups at baseline, with the highest K(trans) value being observed in patients with SpA (median 0.050/minute, interquartile range [IQR] 0.041- 0.069). Furthermore, the K(trans) , kep , and ve values correlated significantly with markers of disease activity. Finally, the PKM parameters K(trans) and kep , but not ve , correlated significantly with synovial expression of vWF (r = 0.647, P = 0.004 for K(trans) ; r = 0.614, P = 0.007 for kep ; r = 0.398, P = 0.102 for ve ).

CONCLUSION

These results suggest that the K(trans) , kep , and ve can be used to detect synovial inflammation in patients with early arthritis, and these PKM parameters may be helpful in differential diagnosis. This approach may also be useful in translational research analyzing tissue microcirculation and angiogenesis.

Mathematical Models of Contrast Transport Kinetics for Cancer Diagnostic Imaging: A Review

Angiogenesis plays a fundamental role in cancer growth and the formation of metastasis. Novel cancer therapies aimed at inhibiting angiogenic processes and/or disrupting angiogenic tumor vasculature are currently being developed and clinically tested. The need for earlier and improved cancer diagnosis, and for early evaluation and monitoring of therapeutic response to angiogenic treatment, have led to the development of several imaging methods for in vivo noninvasive assessment of angiogenesis. The combination of dynamic contrast-enhanced imaging with mathematical modeling of the contrast agent kinetics enables quantitative assessment of the structural and functional changes in the microvasculature that are associated with tumor angiogenesis. In this paper, we review quantitative imaging of angiogenesis with dynamic contrast-enhanced magnetic resonance imaging, computed tomography, and ultrasound.

Comparison of actual with default hematocrit value in dynamic contrast enhanced MR perfusion quantification in grading of human glioma

PURPOSE

Dynamic contrast enhanced (DCE) MRI is used to grade and to monitor the progression of glioma while on treatment. Usually, a fixed hematocrit (Hct) value for adults is assumed to be ~45%; however, it is actually known for individual variations. Purpose of this study was to investigate the effect of measured Hct values in glioma grading using DCE-MRI.

MATERIALS AND METHODS

Fifty glioma patients were included in this study. Kinetic and hemodynamic parameters were estimated for each patient using assumed as well as measured Hct values. To look the changes in Hct value over time, Hct was measured multiple times from 10 of these glioma patients who were on treatment. Simulation was done to look for the effect of extreme variations of Hct values on perfusion metrics. The data was compared to look for significant differences in the perfusion metrics derived from assumed and measured Hct values.

RESULTS

The measured Hct value in patients was found to be (40.4±4.28)%. The sensitivity and specificity of DCE-MRI parameters in glioma grading were not significantly influenced by using measured vis-a-vis assumed Hct values. The serial Hct values from 10 patients who were on treatment showed a fluctuation of 15-20% over time. The simulated data showed linear influence of Hct values on kinetic parameters. The tumor grading was altered on altering the Hct values in borderline cases.

CONCLUSION

Hct values influence the hemodynamic and kinetic metrics linearly and may affect glioma grading. However, perfusion metrics values might change significantly with large change in Hct values, especially in patients who are on chemotherapy necessitating its use in the DCE model.

Reproducibility of Dynamic Contrast-Enhanced MRI in Renal Cell Carcinoma: A Prospective Analysis on Intra- and Interobserver and Scan-Rescan Performance of Pharmacokinetic Parameters

The objective of this study was to investigate the intra- and interobserver as well as scan-rescan reproducibility of quantitative parameters of renal cell carcinomas (RCCs) with dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). A total of 21 patients with clear cell RCCs (17 men, 4 woman; age 37-69 years, mean age 54.6 years, mean size, 5.0 ± 2.2 cm) were prospectively recruited from September 2012 to November 2012. Patients underwent paired DCE-MRI studies on a 3.0 T MR system with an interval of 48 to 72 hours. The extended-Tofts model and population-based arterial input function were used to calculate kinetic parameters. Three observers defined the 2-dimensional whole-tumor region of interest at the slice with the maximum diameter of the RCC. Intraobserver and scan-rescan differences were assessed using paired t tests, whereas interobserver differences using two-way analysis of variance. Intra- and interobserver reproducibility and scan-rescan reproducibility were evaluated using within-subject coefficient of variation (wCoV) and intraclass correlation coefficient (ICC). There were no significant intra-, interobserver, or scan-rescan differences in parameters (all P > 0.05). All ICCs for intra- and interobserver agreements were >0.75 (P < 0.05), whereas the scan-rescan agreement was moderate to good; V(e) (0.764, 95% confidence interval [CI]: 0.378-0.925) and K(ep) (0.906, 95% CI: 0.710-0.972) had higher ICC than K(trans) (0.686; 95% CI: 0.212-0.898) and V(p) (0.657; 95% CI: 0.164-0.888). In intra- and interobserver variability analyses, all parameters except V(p) had low wCoV values. K(trans) and V(e) had slightly lower intraobserver wCoV (1.2% and 0.9%) compared with K(ep) (3.7%), whereas all 3 of these parameters had similar interobserver wCoV values (2.5%, 3.1%, and 2.9%, respectively). Regarding scan-rescan variability, K(trans) and K(ep) showed slightly higher variation (15.6% and 15.4%) than V(e) (10.1%). V(p) had the largest wCoV in all variability analyses (all >30%). DCE-MRI demonstrated good intra- and interobserver reproducibility and moderate to good scan-rescan performance in the assessment of RCC using K(trans), K(ep), and V(e) as parameters under noncontinuous scanning mode. V(p) showed poor reproducibility, and thus may not be suitable for this scanning protocol.

Dynamic Contrast-Enhanced Perfusion MRI of High Grade Brain Gliomas Obtained with Arterial or Venous Waveform Input Function

BACKGROUND AND PURPOSE

The aim of this study was to evaluate the differences in dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) perfusion estimates of high-grade brain gliomas (HGG) due to the use of an input function (IF) obtained respectively from arterial (AIF) and venous (VIF) approaches by two different commercially available software applications.

METHODS

This prospective study includes 20 patients with pathologically confirmed diagnosis of high-grade gliomas. The data source was processed by using two DCE dedicated commercial packages, both based on the extended Toft model, but the first customized to obtain input function from arterial measurement and the second from sagittal sinus sampling. The quantitative parametric perfusion maps estimated from the two software packages were compared by means of a region of interest (ROI) analysis. The resulting input functions from venous and arterial data were also compared.

RESULTS

No significant difference has been found between the perfusion parameters obtained with the two different software packages (P-value < .05). The comparison of the VIFs and AIFs obtained by the two packages showed no statistical differences.

CONCLUSIONS

Direct comparison of DCE-MRI measurements with IF generated by means of arterial or venous waveform led to no statistical difference in quantitative metrics for evaluating HGG. However, additional research involving DCE-MRI acquisition protocols and post-processing would be beneficial to further substantiate the effectiveness of venous approach as the IF method compared with arterial-based IF measurement.