Bolus arrival time and its effect on tissue characterization with dynamic contrast-enhanced magnetic resonance imaging

Intravoxel Incoherent Motion Imaging on Sacroiliitis in Patients With Axial Spondyloarthritis: Correlation With Perfusion Characteristics Based on Dynamic Contrast-Enhanced Magnetic Resonance Imaging

BACKGROUND

To prospectively explore the relationship between intravoxel incoherent motion (IVIM) diffusion-weighted imaging (DWI) and dynamic contrast-enhanced MRI (DCE-MRI) parameters of sacroiliitis in patients with axial spondyloarthritis (axSpA).

METHODS

Patients with initially diagnosed axSpA prospectively underwent on 3.0 T MRI of sacroiliac joint (SIJ). The IVIM parameters (D, f, D ^*) were calculated using biexponential analysis. K ^trans, K _ep, V _e, and V _p from DCE-MRI were obtained in SIJ. The uni-variable and multi-variable linear regression analyses were used to evaluate the correlation between the parameters from these two imaging methods after controlling confounders, such as bone marrow edema (BME), age, agenda, scopes, and localization of lesions, and course of the disease. Then, their correlations were measured by calculating the Pearson's correlation coefficient (r).

RESULTS

The study eventually enrolled 234 patients (178 men, 56 women; mean age, 28.51 ± 9.50 years) with axSpA. With controlling confounders, D was independently related to K ^trans (regression coefficient [b] = 27.593, p < 0.001), K _ep (b = -6.707, p = 0.021), and V _e (b = 131.074, p = 0.003), whereas f and D ^* had no independent correlation with the parameters from DCE MRI. The correlations above were exhibited with Pearson's correlation coefficients (r) (r = 0.662, -0.408, and 0.396, respectively, all p < 0.001).

CONCLUSION

There were independent correlations between D derived from IVIM DWI and K ^trans, K _ep, and V _e derived from DCE-MRI. The factors which affect their correlations mainly included BME, gender, and scopes of lesions.

Quantifying lumbar vertebral perfusion by a Tofts model on DCE-MRI using segmental versus aortic arterial input function

The purpose of this study was to investigate the influence of arterial input function (AIF) selection on the quantification of vertebral perfusion using axial dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). In this study, axial DCE-MRI was performed on 2 vertebrae in each of eight healthy volunteers (mean age, 36.9 years; 5 men) using a 1.5-T scanner. The pharmacokinetic parameters K^trans, v_e, and v_p, derived using a Tofts model on axial DCE-MRI of the lumbar vertebrae, were evaluated using various AIFs: the population-based aortic AIF (AIF_PA), a patient-specific aortic AIF (AIF_A) and a patient-specific segmental arterial AIF (AIF_SA). Additionally, peaks and delay times were changed to simulate the effects of various AIFs on the calculation of perfusion parameters. Nonparametric analyses including the Wilcoxon signed rank test and the Kruskal–Wallis test with a Dunn–Bonferroni post hoc analysis were performed. In simulation, K^trans and v_e increased as the peak in the AIF decreased, but v_p increased when delay time in the AIF increased. In humans, the estimated K^trans and v_e were significantly smaller using AIF_A compared to AIF_SA no matter the computation style (pixel-wise or region-of-interest based). Both these perfusion parameters were significantly greater using AIF_SA compared to AIF_A.

Ultrafast dynamic contrast-enhanced breast MRI may generate prognostic imaging markers of breast cancer

Background

Ultrafast dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI)-derived kinetic parameters have demonstrated at least equivalent accuracy to standard DCE-MRI in differentiating malignant from benign breast lesions. However, it is unclear if they have any efficacy as prognostic imaging markers. The aim of this study was to investigate the relationship between ultrafast DCE-MRI-derived kinetic parameters and breast cancer characteristics.

Methods

Consecutive breast MRI examinations between February 2017 and January 2018 were retrospectively reviewed to determine those examinations that meet the following inclusion criteria: (1) BI-RADS 4–6 MRI performed on a 3T scanner with a 16-channel breast coil and (2) a hybrid clinical protocol with 15 phases of ultrafast DCE-MRI (temporal resolution of 2.7–4.6 s) followed by early and delayed phases of standard DCE-MRI. The study included 125 examinations with 142 biopsy-proven breast cancer lesions. Ultrafast DCE-MRI-derived kinetic parameters (maximum slope [MS] and bolus arrival time [BAT]) were calculated for the entire volume of each lesion. Comparisons of these parameters between different cancer characteristics were made using generalized estimating equations, accounting for the presence of multiple lesions per patient. All comparisons were exploratory and adjustment for multiple comparisons was not performed; P values < 0.05 were considered statistically significant.

Results

Significantly larger MS and shorter BAT were observed for invasive carcinoma than ductal carcinoma in situ (DCIS) (P < 0.001 and P = 0.008, respectively). Significantly shorter BAT was observed for invasive carcinomas with more aggressive characteristics than those with less aggressive characteristics: grade 3 vs. grades 1–2 (P = 0.025), invasive ductal carcinoma vs. invasive lobular carcinoma (P = 0.002), and triple negative or HER2 type vs. luminal type (P < 0.001).

Conclusions

Ultrafast DCE-MRI-derived parameters showed a strong relationship with some breast cancer characteristics, especially histopathology and molecular subtype.

Quantitative Imaging Informatics for Cancer Research

PURPOSE

We summarize Quantitative Imaging Informatics for Cancer Research (QIICR; U24 CA180918), one of the first projects funded by the National Cancer Institute (NCI) Informatics Technology for Cancer Research program.

METHODS

QIICR was motivated by the 3 use cases from the NCI Quantitative Imaging Network. 3D Slicer was selected as the platform for implementation of open-source quantitative imaging (QI) tools. Digital Imaging and Communications in Medicine (DICOM) was chosen for standardization of QI analysis outputs. Support of improved integration with community repositories focused on The Cancer Imaging Archive (TCIA). Priorities included improved capabilities of the standard, toolkits and tools, reference datasets, collaborations, and training and outreach.

RESULTS

Fourteen new tools to support head and neck cancer, glioblastoma, and prostate cancer QI research were introduced and downloaded over 100,000 times. DICOM was amended, with over 40 correction proposals addressing QI needs. Reference implementations of the standard in a popular toolkit and standalone tools were introduced. Eight datasets exemplifying the application of the standard and tools were contributed. An open demonstration/connectathon was organized, attracting the participation of academic groups and commercial vendors. Integration of tools with TCIA was improved by implementing programmatic communication interface and by refining best practices for QI analysis results curation.

CONCLUSION

Tools, capabilities of the DICOM standard, and datasets we introduced found adoption and utility within the cancer imaging community. A collaborative approach is critical to addressing challenges in imaging informatics at the national and international levels. Numerous challenges remain in establishing and maintaining the infrastructure of analysis tools and standardized datasets for the imaging community. Ideas and technology developed by the QIICR project are contributing to the NCI Imaging Data Commons currently being developed.

Selection of Fitting Model and Arterial Input Function for Repeatability in Dynamic Contrast-Enhanced Prostate MRI

RATIONALE AND OBJECTIVES

Analysis of dynamic contrast-enhanced (DCE) magnetic resonance imaging is notable for the variability of calculated parameters. The purpose of this study was to evaluate the level of measurement variability and error/variability due to modeling in DCE magnetic resonance imaging parameters.

MATERIALS AND METHODS

Two prostate DCE scans were performed on 11 treatment-naïve patients with suspected or confirmed prostate peripheral zone cancer within an interval of less than two weeks. Tumor-suspicious and normal-appearing regions of interest (ROI) in the prostate peripheral zone were segmented. Different Tofts-Kety based models and different arterial input functions, with and without bolus arrival time (BAT) correction, were used to extract pharmacokinetic parameters. The percent repeatability coefficient (%RC) of fitted model parameters K_trans, v_e, and k_ep was calculated. Paired t-tests comparing parameters in tumor-suspicious ROIs and in normal-appearing tissue evaluated each parameter's sensitivity to pathology.

RESULTS

Although goodness-of-fit criteria favored the four-parameter extended Tofts-Kety model with the BAT correction included, the simplest two-parameter Tofts-Kety model overall yielded the best repeatability scores. The best %RC in the tumor-suspicious ROI was 63% for k_ep, 28% for v_e, and 83% for K_trans . The best p values for discrimination between tissues were p <10^-5 for k_ep and K_trans, and p = 0.11 for v_e. Addition of the BAT correction to the models did not improve repeatability.

CONCLUSION

The parameter k_ep, using an arterial input functions directly measured from blood signals, was more repeatable than K_trans. Both K_trans and k_ep values were highly discriminatory between healthy and diseased tissues in all cases. The parameter v_e had high repeatability but could not distinguish the two tissue types.

Bolus arrival time estimation in dynamic contrast-enhanced magnetic resonance imaging of small animals based on spline models

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is used to quantify perfusion and vascular permeability. In most cases a bolus arrival time (BAT) delay exists between the arterial input function (AIF) and the contrast agent arrival in the tissue of interest which needs to be estimated. Existing methods for BAT estimation are tailored to tissue concentration curves, which have a fast upslope to the peak as frequently observed in patient data. However, they may give poor results for curves that do not have this characteristic shape such as tissue concentration curves of small animals. In this paper, we propose a method for BAT estimation of signals that do not have a fast upslope to their peak. The model is based on splines which are able to adapt to a large variety of concentration curves. Furthermore, the method estimates BATs on a continuous time scale. All relevant model parameters are automatically determined by generalized cross validation. We use simulated concentration curves of small animal and patient settings to assess the accuracy and robustness of our approach. The proposed method outperforms a state-of-the-art method for small animal data and it gives competitive results for patient data. Finally, it is tested on in vivo acquired rat data where accuracy of BAT estimation was also improved upon the state-of-the-art method. The results indicate that the proposed method is suitable for accurate BAT estimation of DCE-MRI data, especially for small animals.

Exercise-stimulated arterial transit time in calf muscles measured by dynamic contrast-enhanced magnetic resonance imaging

The primary goal of this study was to evaluate arterial transit time (ATT) in exercise-stimulated calf muscles as a promising indicator of muscle function. Following plantar flexion, ATT was measured by dynamic contrast-enhanced (DCE) MRI in young and elderly healthy subjects and patients with peripheral artery disease (PAD). In the young healthy subjects, gastrocnemius ATT decreased significantly (P < 0.01) from 4.3 ± 1.5 to 2.4 ± 0.4 sec when exercise load increased from 4 lbs to 16 lbs. For the same load of 4 lbs, gastrocnemius ATT was lower in the elderly healthy subjects (3.2 ± 1.1 sec; P = 0.08) and in the PAD patients (2.4 ± 1.2 sec; P = 0.02) than in the young healthy subjects. While the sensitivity of the exercise-stimulated ATT is diagnostically useful, it poses a challenge for arterial spin labeling (ASL), a noncontrast MRI method for measuring muscle perfusion. As a secondary goal of this study, we assessed the impact of ATT on ASL-measured perfusion with ASL data of multiple post labeling delays (PLDs) acquired from a healthy subject. Perfusion varied substantially with PLD in the activated gastrocnemius, which can be attributed to the ATT variability as verified by a simulation. In conclusion, muscle ATT is sensitive to exercise intensity, and it potentially reflects the functional impact of aging and PAD on calf muscles. For precise measurement of exercise-stimulated muscle perfusion, it is recommended that ATT be considered when quantifying muscle ASL data.

Prostate DCE-MRI with <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>B</mml:mi> <mml:mn>1</mml:mn> <mml:mo>+</mml:mo></mml:msubsup> </mml:mrow> </mml:math> correction using an approximated analytical approach

PURPOSE

To develop and evaluate a practical <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>B</mml:mi> <mml:mn>1</mml:mn> <mml:mo>+</mml:mo></mml:msubsup> </mml:mrow> </mml:math> correction method for prostate dynamic contrast-enhanced (DCE) MRI analysis.

THEORY

We proposed a simple analytical <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>B</mml:mi> <mml:mn>1</mml:mn> <mml:mo>+</mml:mo></mml:msubsup> </mml:mrow> </mml:math> correction method using a Taylor series approximation to the steady-state spoiled gradient echo signal equation. This approach only requires <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>B</mml:mi> <mml:mn>1</mml:mn> <mml:mo>+</mml:mo></mml:msubsup> </mml:mrow> </mml:math> maps and uncorrected pharmacokinetic (PK) parameters as input to estimate the corrected PK parameters.

METHODS

The proposed method was evaluated using a prostate digital reference object (DRO), and 82 in vivo prostate DCE-MRI cases. The approximated analytical correction was compared with the ground truth PK parameters in simulation, and compared with the reference numerical correction in in vivo experiments, using percentage error as the metric.

RESULTS

The prostate DRO results showed that our approximated analytical approach provided residual error less than 0.4% for both K^trans and v_e , compared to the ground truth. This noise-free residual error was smaller than the noise-induced error using the reference numerical correction, which had a minimum error of 2.1+4.3% with baseline signal-to-noise ratio of 234.5. For the 82 in vivo cases, K^trans and v_e percentage error compared to the reference numerical correction method had a mean of 0.1% (95% central range of [0.0%, 0.2%]) across the prostate volume.

CONCLUSION

The approximated analytical <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>B</mml:mi> <mml:mn>1</mml:mn> <mml:mo>+</mml:mo></mml:msubsup> </mml:mrow> </mml:math> correction method provides comparable results with less than 0.2% error within 95% central range, compared to reference numerical <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>B</mml:mi> <mml:mn>1</mml:mn> <mml:mo>+</mml:mo></mml:msubsup> </mml:mrow> </mml:math> correction. The proposed method is a practical solution for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>B</mml:mi> <mml:mn>1</mml:mn> <mml:mo>+</mml:mo></mml:msubsup> </mml:mrow> </mml:math> correction in prostate DCE-MRI because of its simple implementation.

Vastly accelerated linear least-squares fitting with numerical optimization for dual-input delay-compensated quantitative liver perfusion mapping

PURPOSE

To propose an efficient algorithm to perform dual input compartment modeling for generating perfusion maps in the liver.

METHODS

We implemented whole field-of-view linear least squares (LLS) to fit a delay-compensated dual-input single-compartment model to very high temporal resolution (four frames per second) contrast-enhanced 3D liver data, to calculate kinetic parameter maps. Using simulated data and experimental data in healthy subjects and patients, whole-field LLS was compared with the conventional voxel-wise nonlinear least-squares (NLLS) approach in terms of accuracy, performance, and computation time.

RESULTS

Simulations showed good agreement between LLS and NLLS for a range of kinetic parameters. The whole-field LLS method allowed generating liver perfusion maps approximately 160-fold faster than voxel-wise NLLS, while obtaining similar perfusion parameters.

CONCLUSIONS

Delay-compensated dual-input liver perfusion analysis using whole-field LLS allows generating perfusion maps with a considerable speedup compared with conventional voxel-wise NLLS fitting. Magn Reson Med 79:2415-2421, 2018. © 2017 International Society for Magnetic Resonance in Medicine.