Effects of inflow and radiofrequency spoiling on the arterial input function in dynamic contrast-enhanced MRI: A combined phantom and simulation study: Effects of Inflow and RF Spoiling on the AIF

Improved reliability of perfusion estimation in dynamic susceptibility contrast MRI by using the arterial input function from dynamic contrast enhanced MRI

The arterial input function (AIF) plays a crucial role in estimating quantitative perfusion properties from dynamic susceptibility contrast (DSC) MRI. An important issue, however, is that measuring the AIF in absolute contrast-agent concentrations is challenging, due to uncertainty in relation to the measuredR 2 ∗ -weighted signal, signal depletion at high concentration, and partial-volume effects. A potential solution could be to derive the AIF from separately acquired dynamic contrast enhanced (DCE) MRI data. We aim to compare the AIF determined from DCE MRI with the AIF from DSC MRI, and estimated perfusion coefficients derived from DSC data using a DCE-driven AIF with perfusion coefficients determined using a DSC-based AIF. AIFs were manually selected in branches of the middle cerebral artery (MCA) in both DCE and DSC data in each patient. In addition, a semi-automatic AIF-selection algorithm was applied to the DSC data. The amplitude and full width at half-maximum of the AIFs were compared statistically using the Wilcoxon rank-sum test, applying a 0.05 significance level. Cerebral blood flow (CBF) was derived with different AIF approaches and compared further. The results showed that the AIFs extracted from DSC scans yielded highly variable peaks across arteries within the same patient. The semi-automatic DSC-AIF had significantly narrower width compared with the manual AIFs, and a significantly larger peak than the manual DSC-AIF. Additionally, the DCE-based AIF provided a more stable measurement of relative CBF and absolute CBF values estimated with DCE-AIFs that were compatible with previously reported values. In conclusion, DCE-based AIFs were reproduced significantly better across vessels, showed more realistic profiles, and delivered more stable and reasonable CBF measurements. The DCE-AIF can, therefore, be considered as an alternative AIF source for quantitative perfusion estimations in DSC MRI.

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.

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.

Robust arterial input function surrogate measurement from the superior sagittal sinus complex signal for fast dynamic contrast-enhanced MRI in the brain

PURPOSE

Accurately estimating the arterial input function for dynamic contrast-enhanced MRI is challenging. An arterial input function is typically determined from signal magnitude changes related to a contrast agent, often leading to underestimation of peak concentrations. Alternatively, signal phase recovers the accurate peak concentration for straight vessels but suffers from high noise. A recent method proposed to fit the signal in the complex plane by combining the advantages of the previous 2 methods. The purpose of this work is to refine this complex-based method to determine the venous output function (VOF), an arterial input function surrogate, from the superior sagittal sinus.

METHODS

We propose a state-of-the-art complex-based method that includes direct compensation for blood inflow and signal phase correction accounting for the curvature of the superior sagittal sinus, generally assumed collinear with B₀ . We compared the magnitude-, phase-, and complex-based VOF determination methods against various simulated biases as well as for 29 brain metastases patients.

RESULTS

Angulation of the superior sagittal sinus relative to B₀ varied widely within patients, and its effect on the signal phase caused an underestimation of peak concentrations of up to 65%. Correction significantly increased the VOF peak concentration for the phase- and complex-based VOFs in the cohort. The phase-based method recovered accurate peak concentrations but lacked precision in the tail of the VOF. Our complex-based VOF completely recovered the effect of inflow and resulted in a high-peak concentration with limited noise.

CONCLUSION

The new complex-based method resulted in high-quality VOF robust against superior sagittal sinus curvature and variations in patient positioning.

Quantifying MRI T1 relaxation in flowing blood: implications for arterial input function measurement in DCE-MRI

OBJECTIVES

To investigate the feasibility of accurately quantifying the concentration of MRI contrast agent in flowing blood by measuring its T1 in a large vessel. Such measures are often used to obtain patient-specific arterial input functions for the accurate fitting of pharmacokinetic models to dynamic contrast enhanced MRI data. Flow is known to produce errors with this technique, but these have so far been poorly quantified and characterised in the context of pulsatile flow with a rapidly changing T1 as would be expected in vivo.

METHODS

A phantom was developed which used a mechanical pump to pass fluid at physiologically relevant rates. Measurements of T₁ were made using high temporal resolution gradient recalled sequences suitable for DCE-MRI of both constant and pulsatile flow. These measures were used to validate a virtual phantom that was then used to simulate the expected errors in the measurement of an AIF in vivo.

RESULTS

The relationship between measured T1 values and flow velocity was found to be non-linear. The subsequent error in quantification of contrast agent concentration in a measured AIF was shown.

CONCLUSIONS

The T1 measurement of flowing blood using standard DCE- MRI sequences are subject to large measurement errors which are non-linear in relation to flow velocity.

ADVANCES IN KNOWLEDGE

This work qualitatively and quantitatively demonstrates the difficulties of accurately measuring the T1 of flowing blood using DCE-MRI over a wide range of physiologically realistic flow velocities and pulsatilities. Sources of error are identified and proposals made to reduce these.

Phantom Validation of DCE-MRI Magnitude and Phase-Based Vascular Input Function Measurements

Accurate, patient-specific measurement of arterial input functions (AIF) may improve model-based analysis of vascular permeability. This study investigated factors affecting AIF measurements from magnetic resonance imaging (MRI) magnitude (AIF_MAGN) and phase (AIF_PHA) signals, and compared them against computed tomography (CT) (AIF_CT), under controlled conditions relevant to clinical protocols using a multimodality flow phantom. The flow phantom was applied at flip angles of 20° and 30°, flow rates (3-7.5 mL/s), and peak bolus concentrations (0.5-10 mM), for in-plane and through-plane flow. Spatial 3D-FLASH signal and variable flip angle T1 profiles were measured to investigate in-flow and radiofrequency-related biases, and magnitude- and phase-derived Gd-DTPA concentrations were compared. MRI AIF performance was tested against AIF_CT via Pearson correlation analysis. AIF_MAGN was sensitive to imaging orientation, spatial location, flip angle, and flow rate, and it grossly underestimated AIF_CT peak concentrations. Conversion to Gd-DTPA concentration using T1 taken at the same orientation and flow rate as the dynamic contrast-enhanced acquisition improved AIF_MAGN accuracy; yet, AIF_MAGN metrics remained variable and significantly reduced from AIF_CT at concentrations above 2.5 mM. AIF_PHA performed equivalently within 1 mM to AIF_CT across all tested conditions. AIF_PHA, but not AIF_MAGN, reported equivalent measurements to AIF_CT across the range of tested conditions. AIF_PHA showed superior robustness.

Reference region extraction by clustering for the pharmacokinetic analysis of dynamic contrast-enhanced MRI in prostate cancer

PURPOSE

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) measures changes in the concentration of an administered contrast agent to quantitatively evaluate blood circulation in a tumor or normal tissues. This method uses a pharmacokinetic analysis based on the time course of a reference region, such as muscle, rather than arterial input function. However, it is difficult to manually define a homogeneous reference region. In the present study, we developed a method for automatic extraction of the reference region using a clustering algorithm based on a time course pattern for DCE-MRI studies of patients with prostate cancer.

METHODS

Two feature values related to the shape of the time course were extracted from the time course of all voxels in the DCE-MRI images. Each voxel value of T1-weighted images acquired before administration were also added as anatomical data. Using this three-dimensional feature vector, all voxels were segmented into five clusters by the Gaussian mixture model, and one of these clusters that included the gluteus muscle was selected as the reference region.

RESULTS

Each region of arterial vessel, muscle, and fat was segmented as a different cluster from the tumor and normal tissues in the prostate. In the extracted reference region, other tissue elements including scattered fat and blood vessels were removed from the muscle region.

CONCLUSIONS

Our proposed method can automatically extract the reference region using the clustering algorithm with three types of features based on the time course pattern and anatomical data. This method may be useful for evaluating tumor circulatory function in DCE-MRI studies.

Pharmacokinetic modeling of dynamic contrast-enhanced MRI using a reference region and input function tail

PURPOSE

Quantitative analysis of dynamic contrast-enhanced MRI (DCE-MRI) requires an arterial input function (AIF) which is difficult to measure. We propose the reference region and input function tail (RRIFT) approach which uses a reference tissue and the washout portion of the AIF.

METHODS

RRIFT was evaluated in simulations with 100 parameter combinations at various temporal resolutions (5-30 s) and noise levels (σ = 0.01-0.05 mM). RRIFT was compared against the extended Tofts model (ETM) in 8 studies from patients with glioblastoma multiforme. Two versions of RRIFT were evaluated: one using measured patient-specific AIF tails, and another assuming a literature-based AIF tail.

RESULTS

RRIFT estimated the transfer constant K trans and interstitial volume v e with median errors within 20% across all simulations. RRIFT was more accurate and precise than the ETM at temporal resolutions slower than 10 s. The percentage error of K trans had a median and interquartile range of -9 ± 45% with the ETM and -2 ± 17% with RRIFT at a temporal resolution of 30 s under noiseless conditions. RRIFT was in excellent agreement with the ETM in vivo, with concordance correlation coefficients (CCC) of 0.95 for K trans , 0.96 for v e , and 0.73 for the plasma volume v p using a measured AIF tail. With the literature-based AIF tail, the CCC was 0.89 for K trans , 0.93 for v e and 0.78 for v p .

CONCLUSIONS

Quantitative DCE-MRI analysis using the input function tail and a reference tissue yields absolute kinetic parameters with the RRIFT method. This approach was viable in simulation and in vivo for temporal resolutions as low as 30 s.

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.

A functional form for a representative individual arterial input function measured from a population using high temporal resolution DCE MRI

PURPOSE

To measure the arterial input function (AIF), an essential component of tracer kinetic analysis, in a population of patients using an optimized dynamic contrast-enhanced (DCE) imaging sequence and to estimate inter- and intrapatient variability. From these data, a representative AIF that may be used for realistic simulation studies can be extracted.

METHODS

Thirty-nine female patients were imaged on multiple visits before and during a course of neoadjuvant chemotherapy for breast cancer. A total of 97 T₁ -weighted DCE studies were analyzed including bookend estimates of T₁ and model-fitting to each individual AIF. Area under the curve and cardiac output were estimated from each first pass peak, and these data were used to assess inter- and intrapatient variability of the AIF.

RESULTS

Interpatient variability exceeded intrapatient variability of the AIF. There was no change in cardiac output as a function of MR visit (mean value 5.6 ± 1.1 L/min) but baseline blood T₁ increased significantly following the start of chemotherapy (which was accompanied by a decrease in hematocrit).

CONCLUSION

The AIF in an individual patient can be measured reproducibly but the variability of AIFs between patients suggests that use of a population AIF will decrease the precision of tracer kinetic analysis performed in cross-patient comparison studies. A representative AIF is presented that is typical of the population but retains the characteristics of an individually measured AIF.

An extended reference region model for DCE-MRI that accounts for plasma volume

The reference region model (RRM) for dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides pharmacokinetic parameters without requiring the arterial input function. A limitation of the RRM is that it assumes that the blood plasma volume in the tissue of interest is zero, but this is often not true in highly vascularized tissues, such as some tumours. This study proposes an extended reference region model (ERRM) to account for tissue plasma volume. Furthermore, ERRM was combined with a two-fit approach to reduce the number of fitting parameters, and this was named the constrained ERRM (CERRM). The accuracy and precision of RRM, ERRM and CERRM were evaluated in simulations covering a range of parameters, noise and temporal resolutions. These models were also compared with the extended Tofts model (ETM) on in vivo glioblastoma multiforme data. In simulations, RRM overestimated K^trans by over 10% at v_p  = 0.01 under noiseless conditions. In comparison, ERRM and CERRM were both accurate, with CERRM showing better precision when noise was included. On in vivo data, CERRM provided maps that had the highest agreement with ETM, whilst also being robust at temporal resolutions as poor as 30 s. ERRM can provide pharmacokinetic parameters without an arterial input function in tissues with non-negligible v_p where RRM provides inaccurate estimates. The two-fit approach, named CERRM, further improves on the accuracy and precision of ERRM.

Estimating the arterial input function from dynamic contrast-enhanced MRI data with compensation for flow enhancement (II): Applications in spine diagnostics and assessment of crohn's disease

BACKGROUND

Pharmacokinetic (PK) models can describe microvascular density and integrity. An essential component of PK models is the arterial input function (AIF) representing the time-dependent concentration of contrast agent (CA) in the blood plasma supplied to a tissue.

PURPOSE/HYPOTHESIS

To evaluate a novel method for subject-specific AIF estimation that takes inflow effects into account.

STUDY TYPE

Retrospective study.

SUBJECTS

Thirteen clinical patients referred for spine-related complaints; 21 patients from a study into luminal Crohn's disease with known Crohn's Disease Endoscopic Index of Severity (CDEIS).

FIELD STRENGTH/SEQUENCE

Dynamic fast spoiled gradient echo (FSPGR) at 3T.

ASSESSMENT

A population-averaged AIF, AIFs derived from distally placed regions of interest (ROIs), and the new AIF method were applied. Tofts' PK model parameters (including v_p and K^trans ) obtained with the three AIFs were compared. In the Crohn's patients K^trans was correlated to CDEIS.

STATISTICAL TESTS

The median values of the PK model parameters from the three methods were compared using a Mann-Whitney U-test. The associated variances were statistically assessed by the Brown-Forsythe test. Spearman's rank correlation coefficient was computed to test the correlation of K^trans to CDEIS.

RESULTS

The median v_p was significantly larger when using the distal ROI approach, compared to the two other methods (P < 0.05 for both comparisons, in both applications). Also, the variances in v_p were significantly larger with the ROI approach (P < 0.05 for all comparisons). In the Crohn's disease study, the estimated K^trans parameter correlated better with the CDEIS (r = 0.733, P < 0.001) when the proposed AIF was used, compared to AIFs from the distal ROI method (r = 0.429, P = 0.067) or the population-averaged AIF (r = 0.567, P = 0.011).

DATA CONCLUSION

The proposed method yielded realistic PK model parameters and improved the correlation of the K^trans parameter with CDEIS, compared to existing approaches.

LEVEL OF EVIDENCE

3 Technical Efficacy Stage 1 J. Magn. Reson. Imaging 2018;47:1197-1204.

Estimating the arterial input function from dynamic contrast-enhanced MRI data with compensation for flow enhancement (I): Theory, method, and phantom experiments

BACKGROUND

The arterial input function (AIF) represents the time-dependent arterial contrast agent (CA) concentration that is used in pharmacokinetic modeling.

PURPOSE

To develop a novel method for estimating the AIF from dynamic contrast-enhanced (DCE-) MRI data, while compensating for flow enhancement.

STUDY TYPE

Signal simulation and phantom measurements.

PHANTOM MODEL

Time-intensity curves (TICs) were simulated for different numbers of excitation pulses modeling flow effects. A phantom experiment was performed in which a solution (without CA) was passed through a straight tube, at constant flow velocity.

FIELD STRENGTH/SEQUENCE

Dynamic fast spoiled gradient echo (FSPGRs) at 3T MRI, both in the simulations and in the phantom experiment. TICs were generated for a duration of 373 seconds and sampled at intervals of 1.247 seconds (300 timepoints).

ASSESSMENT

The proposed method first estimates the number of pulses that spins have received, and then uses this knowledge to accurately estimate the CA concentration.

STATISTICAL TESTS

The difference between the median of the estimated number of pulses and the true value was determined, as well as the interquartile range (IQR) of the estimations. The estimated CA concentrations were evaluated in the same way. The estimated number of pulses was also used to calculate flow velocity.

RESULTS

The difference between the median estimated and reference number of pulses varied from -0.005 to -1.371 (corresponding IQRs: 0.853 and 48.377) at true values of 10 and 180 pulses, respectively. The difference between the median estimated CA concentration and the reference value varied from -0.00015 to 0.00306 mmol/L (corresponding IQRs: 0.01989 and 1.51013 mmol/L) at true values of 0.5 and 8.0 mmol/l, respectively, at an intermediate value of 100 pulses. The estimated flow velocities in the phantom were within 10% of the reference value.

DATA CONCLUSION

The proposed method accurately corrects the MRI signal affected by the inflow effect.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;47:1190-1196.

Semi-parametric arterial input functions for quantitative dynamic contrast enhanced magnetic resonance imaging in mice

OBJECTIVE

An extension of single- and multi-channel blind deconvolution is presented to improve the estimation of the arterial input function (AIF) in quantitative dynamic contrast enhanced magnetic resonance imaging (DCE-MRI).

METHODS

The Lucy-Richardson expectation-maximization algorithm is used to obtain estimates of the AIF and the tissue residue function (TRF). In the first part of the algorithm, nonparametric estimates of the AIF and TRF are obtained. In the second part, the decaying part of the AIF is approximated by three decaying exponential functions with the same delay, giving an almost noise free semi-parametric AIF. Simultaneously, the TRF is approximated using the adiabatic approximation of the Johnson-Wilson (aaJW) pharmacokinetic model.

RESULTS

In simulations and tests on real data, use of this AIF gave perfusion values close to those obtained with the corresponding previously published nonparametric AIF, and are more noise robust.

CONCLUSION

When used subsequently in voxelwise perfusion analysis, these semi-parametric AIFs should give more correct perfusion analysis maps less affected by recording noise than the corresponding nonparametric AIFs, and AIFs obtained from arteries.

SIGNIFICANCE

This paper presents a method to increase the noise robustness in the estimation of the perfusion parameter values in DCE-MRI.

Increased robustness in reference region model analysis of DCE MRI using two-step constrained approaches

PURPOSE

Reference region models (RRMs) can quantify tumor perfusion in dynamic contrast-enhanced MRI without an arterial input function. Inspection of the RRM reveals that one of the free parameters in the fit is uniquely linked to the reference region and is common to all voxels. A two-step approach is proposed that takes this constraint into account.

METHODS

Three constrained RRM (CRRM) approaches were devised and evaluated. Simulations were performed to compare their accuracy and precision over a range of noise and temporal resolutions. The CRRM was also applied on a virtual phantom that simulates different perfusion values. In vivo evaluation was performed on data from breast cancer and soft tissue sarcoma.

RESULTS

In simulations, the CRRM consistently improved precision and had better accuracy at low signal-to-noise ratio (SNR). In virtual phantom, the CRRMs were able to fit voxels that had similar kinetics to the reference tissue, whereas the unconstrained models failed to accurately fit these voxels. In the in vivo data, the constrained approaches produced parameter maps that had less variability and were in better agreement with the Tofts model.

CONCLUSION

These findings indicate that the two-step fitting approach of the CRRM can reduce the variability of perfusion estimates for quantifying perfusion with dynamic contrast-enhanced (DCE) MRI. Magn Reson Med 78:1547-1557, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Reduced respiratory motion artifacts using structural similarity in fast 2D dynamic contrast enhanced MRI of liver lesions

The purpose of this work was to improve dynamic contrast enhanced MRI (DCE-MRI) of liver lesions by removing motion corrupted images as identified by a structural similarity (SSIM) algorithm, and to assess the effect of this correction on the pharmacokinetic parameter K^trans using automatically determined arterial input functions (AIFs). Fifteen patients with colorectal liver metastases were measured twice with a T₁ weighted multislice 2D FLASH sequence for DCE-MRI (time resolution 1.2 s). AIFs were automatically derived from contrast inflow in the aorta of each patient. Thereafter, SSIM identified motion corrupted images of the liver were removed from the DCE dataset. From this corrected data set K^trans and its reproducibility were determined. Using the SSIM algorithm a median fraction of 46% (range 37-50%) of the liver images in DCE time series was labeled as motion distorted. Rejection of these images resulted in a significantly lower median K^trans (p < 0.05) and lower coefficient of repeatability of K^trans in liver metastases compared with an analysis without correction. SSIM correction improves the reproducibility of the DCE-MRI parameter K^trans in liver metastasis and reduces contamination of K^trans values of lesions by that of surrounding normal liver tissue.

Phase-based vascular input function: Improved quantitative DCE-MRI of atherosclerotic plaques

PURPOSE

Quantitative pharmacokinetic modeling of dynamic contrast-enhanced (DCE)-MRI can be used to assess atherosclerotic plaque microvasculature, which is an important marker of plaque vulnerability. Purpose of the present study was (1) to compare magnitude- versus phase-based vascular input functions (m-VIF vs ph-VIF) used in pharmacokinetic modeling and (2) to perform model calculations and flow phantom experiments to gain more insight into the differences between m-VIF and ph-VIF.

METHODS

Population averaged m-VIF and ph-VIFs were acquired from 11 patients with carotid plaques and used for pharmacokinetic analysis in another 17 patients. Simulations, using the Bloch equations and the MRI scan geometry, and flow phantom experiments were performed to determine the effect of local blood velocity on the magnitude and phase signal enhancement.

RESULTS

Simulations and flow phantom experiments revealed that flow within the lumen can lead to severe underestimation of m-VIF, while this is not the case for the ph-VIF. In line, the peak concentration of the m-VIF is significantly lower than ph-VIF (p < 0.001), in vivo. Quantitative model parameters for m- and ph-VIF differed in absolute values but were moderate to strongly correlated with each other [K(trans) Spearman's ρ > 0.93 (p < 0.001) and vp Spearman's ρ > 0.58 (p < 0.05)].

CONCLUSIONS

m-VIF is strongly influenced by local blood velocity, which leads to underestimation of the contrast medium concentration. Therefore, it is advised to use ph-VIF for DCE-MRI analysis of carotid plaques for accurate quantification.

Improved accuracy and precision of tracer kinetic parameters by joint fitting to variable flip angle and dynamic contrast enhanced MRI data

PURPOSE

To improve the accuracy and precision of tracer kinetic model parameter estimates for use in dynamic contrast enhanced (DCE) MRI studies of solid tumors.

THEORY

Quantitative DCE-MRI requires an estimate of precontrast T1 , which is obtained prior to fitting a tracer kinetic model. As T1 mapping and tracer kinetic signal models are both a function of precontrast T1 it was hypothesized that its joint estimation would improve the accuracy and precision of both precontrast T1 and tracer kinetic model parameters.

METHODS

Accuracy and/or precision of two-compartment exchange model (2CXM) parameters were evaluated for standard and joint fitting methods in well-controlled synthetic data and for 36 bladder cancer patients. Methods were compared under a number of experimental conditions.

RESULTS

In synthetic data, joint estimation led to statistically significant improvements in the accuracy of estimated parameters in 30 of 42 conditions (improvements between 1.8% and 49%). Reduced accuracy was observed in 7 of the remaining 12 conditions. Significant improvements in precision were observed in 35 of 42 conditions (between 4.7% and 50%). In clinical data, significant improvements in precision were observed in 18 of 21 conditions (between 4.6% and 38%).

CONCLUSION

Accuracy and precision of DCE-MRI parameter estimates are improved when signal models are fit jointly rather than sequentially. Magn Reson Med 76:1270-1281, 2016. © 2015 Wiley Periodicals, Inc.

Sensitivity of quantitative myocardial dynamic contrast-enhanced MRI to saturation pulse efficiency, noise and t1 measurement error: Comparison of nonlinearity correction methods

PURPOSE

To compare methods designed to minimize or correct signal nonlinearity in quantitative myocardial dynamic contrast-enhanced (DCE) MRI.

METHODS

DCE-MRI studies were simulated and data acquired in eight volunteers. Signal nonlinearity was corrected using either a dual-bolus approach or model-based correction using proton-density weighted imaging (conventional or dual-sequence acquisition) or T1 data (native or bookend). Scanning of healthy and infarcted myocardium at 3 T was simulated, including noise, saturation imperfection and T1 measurement error. Data were analyzed using model-based deconvolution with a one-compartment (mono-exponential) model.

RESULTS

Substantial variation between methods was demonstrated in volunteers. In simulations the dual-bolus method proved stable for realistic levels of saturation efficiency but demonstrated bias due to residual nonlinearity. Model-based methods performed ideally in the absence of confounding error sources and were generally robust to noise or saturation imperfection, except for native T1 based correction which was highly sensitive to the latter. All methods demonstrated large variation in accuracy above an over-saturation level where baseline signal was nulled. For the dual-sequence approach this caused substantial bias at the saturation efficiencies observed in volunteers.

CONCLUSION

The choice of nonlinearity correction method in myocardial DCE-MRI impacts on accuracy and precision of estimated parameters, particularly in the presence of nonideal saturation.

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

Pharmacokinetic modeling in Dynamic Contrast Enhanced (DCE)-MRI is an elegant and useful method that provides valuable insight into angiogenesis in cancer and inflammatory diseases. Despite its widespread use, the reliability of the model results is still questioned, as many factors hamper the calculation of the model's parameters, resulting in the poor reproducibility and accuracy of the method. Pharmacokinetic modeling relies on the knowledge of inputs such as the Arterial Input Function (AIF) and of the tissue contrast agent concentration, both of which are difficult to accurately measure. Any errors in the measurement of either of the inputs propagate into the calculated pharmacokinetic model parameters (PMPs), and the significance of the effect depends on the source of the measurement error. In this work we systematically investigate the effect of the incorrect estimation of the parameters describing the inputs of the model on the calculated PMPs when using Tofts' model. Furthermore, we analyze the dependence of these errors on the native values of the PMPs. We show that errors on the measurement of the native T1 as well as errors on the parameters describing the initial peak of the AIF have the largest impact on the calculated PMPs. The parameter whose error has the least effect is the one describing the slow decay of the AIF. The effect of input parameter (IP) errors on the calculated PMPs is found to be dependent on the native set of PMPs: this is particularly true for the errors in the flip angle, and for the errors in parameters describing the initial AIF peak. Conversely the effect of T1 and AIF scaling errors on the calculated PMPs is only slightly dependent on the native PMPs.

Comparison of dynamic contrast-enhanced MRI and quantitative SPECT in a rat glioma model

Pharmacokinetic modeling of dynamic contrast-enhanced (DCE) MRI data provides measures of the extracellular-extravascular volume fraction (v(e) ) and the volume transfer constant (K(trans) ) in a given tissue. These parameter estimates may be biased, however, by confounding issues such as contrast agent and tissue water dynamics, or assumptions of vascularization and perfusion made by the commonly used model. In contrast to MRI, radiotracer imaging with SPECT is insensitive to water dynamics. A quantitative dual-isotope SPECT technique was developed to obtain an estimate of v(e) in a rat glioma model for comparison with the corresponding estimates obtained using DCE-MRI with a vascular input function and reference region model. Both DCE-MRI methods produced consistently larger estimates of v(e) in comparison to the SPECT estimates, and several experimental sources were postulated to contribute to these differences.

Classic models for dynamic contrast-enhanced MRI

Dynamic contrast-enhanced MRI (DCE-MRI) is a functional MRI method where T1 -weighted MR images are acquired dynamically after bolus injection of a contrast agent. The data can be interpreted in terms of physiological tissue characteristics by applying the principles of tracer-kinetic modelling. In the brain, DCE-MRI enables measurement of cerebral blood flow (CBF), cerebral blood volume (CBV), blood-brain barrier (BBB) permeability-surface area product (PS) and the volume of the interstitium (ve ). These parameters can be combined to form others such as the volume-transfer constant K(trans) , the extraction fraction E and the contrast-agent mean transit times through the intra- and extravascular spaces. A first generation of tracer-kinetic models for DCE-MRI was developed in the early 1990s and has become a standard in many applications. Subsequent improvements in DCE-MRI data quality have driven the development of a second generation of more complex models. They are increasingly used, but it is not always clear how they relate to the models of the first generation or to the model-free deconvolution methods for tissues with intact BBB. This lack of understanding is leading to increasing confusion on when to use which model and how to interpret the parameters. The purpose of this review is to clarify the relation between models of the first and second generations and between model-based and model-free methods. All quantities are defined using a generic terminology to ensure the widest possible scope and to reveal the link between applications in the brain and in other organs.

Tracer kinetic modelling in MRI: estimating perfusion and capillary permeability

The tracer-kinetic models developed in the early 1990s for dynamic contrast-enhanced MRI (DCE-MRI) have since become a standard in numerous applications. At the same time, the development of MRI hardware has led to increases in image quality and temporal resolution that reveal the limitations of the early models. This in turn has stimulated an interest in the development and application of a second generation of modelling approaches. They are designed to overcome these limitations and produce additional and more accurate information on tissue status. In particular, models of the second generation enable separate estimates of perfusion and capillary permeability rather than a single parameter K(trans) that represents a combination of the two. A variety of such models has been proposed in the literature, and development in the field has been constrained by a lack of transparency regarding terminology, notations and physiological assumptions. In this review, we provide an overview of these models in a manner that is both physically intuitive and mathematically rigourous. All are derived from common first principles, using concepts and notations from general tracer-kinetic theory. Explicit links to their historical origins are included to allow for a transfer of experience obtained in other fields (PET, SPECT, CT). A classification is presented that reveals the links between all models, and with the models of the first generation. Detailed formulae for all solutions are provided to facilitate implementation. Our aim is to encourage the application of these tools to DCE-MRI by offering researchers a clearer understanding of their assumptions and requirements.

Phase-based arterial input functions in humans applied to dynamic contrast-enhanced MRI: potential usefulness and limitations

OBJECT

Phase-based arterial input functions (AIFs) provide a promising alternative to standard magnitude-based AIFs, for example, because inflow effects are avoided. The usefulness of phase-based AIFs in clinical dynamic contrast-enhanced MRI (DCE-MRI) was investigated, and relevant pitfalls and sources of uncertainty were identified.

MATERIALS AND METHODS

AIFs were registered from eight human subjects on, in total, 21 occasions. AIF quality was evaluated by comparing AIFs from right and left internal carotid arteries and by assessing the reliability of blood plasma volume estimates.

RESULTS

Phase-based AIFs yielded an average bolus peak of 3.9 mM and a residual concentration of 0.37 mM after 3 min, (0.033 mmol/kg contrast agent injection). The average blood plasma volume was 2.7% when using the AIF peak in the estimation, but was significantly different (p < 0.0001) and less physiologically reasonable when based on the AIF tail concentration. Motion-induced phase shifts and accumulation of contrast agent in background tissue regions were identified as main sources of uncertainty.

CONCLUSION

Phase-based AIFs are a feasible alternative to magnitude AIFs, but sources of errors exist, making quantification difficult, especially of the AIF tail. Improvement of the technique is feasible and also required for the phase-based AIF approach to reach its full potential.