We evaluate the reproducibility of ktrans values derived from dynamic contrast MRI images in patients with newly diagnosed glioblastoma. Particular focus is put on the reproducibility across choices of T1 values and arterial input functions (AIFs) to the Tofts-Kermode model. Reproducibility is assessed across multiple pre-therapy baseline visits in a 45 patient cohort. Our model based on static population AIFs and static global T1 values had excellent reproducibility compared to other models including unique individual AIFs and T1 maps. There is relatively little concordance between all tested models, but individual AIFs led to higher mean ktrans values.
Purpose
We evaluated the repeatability of ktrans obtained from enhancement regions in dynamic contrast-enhanced MRI (DCE-MRI) images of patients with newly diagnosed glioblastoma. Each patient had two pre-therapy baseline scans in the context of a clinical trial. Ktrans data were obtained from the modified Tofts-Kermode model using different arterial input functions (AIFs) and T1 tissue mappings [1]. Four ktrans values were extracted, based on having either a population-based arterial input function (AIF) or an automatically determined individual-level AIF, and either a global assumed T1 value or an individual variable flip angle T1 map. In addition to assessing repeatability, we also assessed the correlations between all combinations of T1 and AIFs. Ktrans has previously been implicated as a predictive imaging biomarker for patients with glioblastoma and other tumors, and assessing their repeatability in vivo, both across patient visits and processing methods, is essential in determining their statistical power in clinical settings [2]. Repeatability has previously been examined for many sites and tissues [3], but few studies have assessed repeatability with a large patient cohort on a 3T scanner. Furthermore, concerns over gadolinium dosage make such an imaging set unlikely to be recreated [4].Methods
Forty-five patients with newly-diagnosed glioblastoma (age: 22-74 years, mean: 56 years) across two IRB approved studies were included in this evaluation. Two baseline scans were acquired before treatment onset, two to four days apart without any interventions or changes in therapy. DCE-MRI images were retrieved with a repetition time of 6.8ms, echo-time of 2.73ms, and a flip angle of 10 degrees all on a 3T Siemens scanner. A neuro-oncologist drew regions of interest outlining the contrast-enhancing tumor on post-contrast T1-weighted images, which were then registered to the space of the DCE-MRI images. Ktrans maps were generated by an in-house Matlab script that performs perfectly on the QIBA v6 DCE-MRI phantom using the modified Tofts-Kermode model [5]. Population-based AIFs were derived from the Parker model [6], individual AIFs were averaged between automatically selected voxels within a manual ROI via nordicICE [7], T1 maps were derived from acquired variable flip angle maps via Freesurfer’s T1 Mapping tool (angles = 2,5,10,15,30) [8], and global static T1 was set to 1500ms. T1 maps, automatic AIFs, and parametric maps were visually inspected for accuracy, and unusable AIFs and images were excluded. This left a double-baseline cohort of n=43 patients for static T1 and population AIF, n=38 for T1 mapping, n=34 for automatic AIFs, and n=26 for automatic AIFs with T1 mapping. Repeatability was assessed via the Bland-Altman repeatability coefficient (RC), a measure of expected error between repeat measurements (CI = .95%), and inter-correlation was assessed via the concordance correlation coefficient (CCC). Scaling differences were assessed by voxel-wide means calculated for voxels with ktrans ranging from 0 to 1.