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.

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 vp and Ktrans ) obtained with the three AIFs were compared. In the Crohn's patients Ktrans 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 Ktrans to CDEIS. RESULTS: The median vp 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 vp were significantly larger with the ROI approach (P < 0.05 for all comparisons). In the Crohn's disease study, the estimated Ktrans 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 Ktrans parameter with CDEIS, compared to existing approaches. LEVEL OF EVIDENCE: 3 Technical Efficacy Stage 1 J. Magn. Reson. Imaging 2018;47:1197-1204.

Feasibility of a fast method for B1-inhomogeneity correction for FSPGR sequences.

The Fast Spoiled Gradient Echo (FSPGR) sequence is often used in MRI to create T1-weighted images. The signal intensity generated by this sequence depends on the applied flip angle. Knowing the correct flip angle is essential for the determination of T1-maps by means of an FSPGR based Variable Flip Angle (VFA) approach. Also, quantitatively determining the concentration of contrast agent in case of Dynamic Contrast Enhanced MRI (DCE-MRI) requires knowledge of the applied flip angle. In both cases, the B1-field (in)homogeneity significantly affects the results. In this paper, we present a new method to obtain both the T1-map and B1-inhomogeneity map using scans that can each be acquired within a breath-hold. We combine two short sequences for T1 quantification: Variable Flip Angle and Look-Locker (LL). The T1-maps obtained from the LL data were used to estimate the B1-inhomogeneity inherently present in the VFA data, which was then used to correct for the VFA method's inaccurate flip angles. This way, a reliable T1-map could be computed, which was validated using both in vitro and in vivo scans. The in vitro results show that the procedure yields a substantially smaller mean deviation in T1 from the T1 measurement's gold standard (the Inversion Recovery method), while the in vivo results show both a more accurate estimation of T1 and a reduction of the influence of the B1-inhomogeneity on the signal intensity.