The Impact of Arterial Input Function Determination Variations on Prostate Dynamic Contrast-Enhanced Magnetic Resonance Imaging Pharmacokinetic Modeling: A Multicenter Data Analysis Challenge, Part II.

This multicenter study evaluated the effect of variations in arterial input function (AIF) determination on pharmacokinetic (PK) analysis of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) data using the shutter-speed model (SSM). Data acquired from eleven prostate cancer patients were shared among nine centers. Each center used a site-specific method to measure the individual AIF from each data set and submitted the results to the managing center. These AIFs, their reference tissue-adjusted variants, and a literature population-averaged AIF, were used by the managing center to perform SSM PK analysis to estimate Ktrans (volume transfer rate constant), ve (extravascular, extracellular volume fraction), kep (efflux rate constant), and τi (mean intracellular water lifetime). All other variables, including the definition of the tumor region of interest and precontrast T1 values, were kept the same to evaluate parameter variations caused by variations in only the AIF. Considerable PK parameter variations were observed with within-subject coefficient of variation (wCV) values of 0.58, 0.27, 0.42, and 0.24 for Ktrans, ve, kep, and τi, respectively, using the unadjusted AIFs. Use of the reference tissue-adjusted AIFs reduced variations in Ktrans and ve (wCV = 0.50 and 0.10, respectively), but had smaller effects on kep and τi (wCV = 0.39 and 0.22, respectively). kep is less sensitive to AIF variation than Ktrans, suggesting it may be a more robust imaging biomarker of prostate microvasculature. With low sensitivity to AIF uncertainty, the SSM-unique τi parameter may have advantages over the conventional PK parameters in a longitudinal study.

DCE-MRI of the prostate using shutter-speed vs. Tofts model for tumor characterization and assessment of aggressiveness.

PURPOSE: To quantify Tofts model (TM) and shutter-speed model (SSM) perfusion parameters in prostate cancer (PCa) and noncancerous peripheral zone (PZ) and to compare the diagnostic performance of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) to Prostate Imaging and Reporting and Data System (PI-RADS) classification for the assessment of PCa aggressiveness. MATERIALS AND METHODS: Fifty PCa patients (mean age 60 years old) who underwent MRI at 3.0T followed by prostatectomy were included in this Institutional Review Board-approved retrospective study. DCE-MRI parameters (Ktrans , ve , kep [TM&SSM] and intracellular water molecule lifetime τi [SSM]) were determined in PCa and PZ. Differences in DCE-MRI parameters between PCa and PZ, and between models were assessed using Wilcoxon signed-rank tests. Receiver operating characteristic (ROC) analysis for differentiation between PCa and PZ was performed for individual and combined DCE-MRI parameters. Diagnostic performance of DCE-MRI parameters for identification of aggressive PCa (Gleason ≥8, grade group [GG] ≥3 or pathology stage pT3) was assessed using ROC analysis and compared with PI-RADSv2 scores. RESULTS: DCE-MRI parameters were significantly different between TM and SSM and between PZ and PCa (P < 0.037). Diagnostic performances of TM and SSM for differentiation of PCa from PZ were similar (highest AUC TM: Ktrans +kep 0.76, SSM: τi +kep 0.80). PI-RADS outperformed TM and SSM DCE-MRI for identification of Gleason ≥8 lesions (AUC PI-RADS: 0.91, highest AUC DCE-MRI: Ktrans +τi SSM 0.61, P = 0.002). The diagnostic performance of PI-RADS and DCE-MRI for identification of GG ≥3 and pT3 PCa was not significantly different (P > 0.213). CONCLUSION: SSM DCE-MRI did not increase the diagnostic performance of DCE-MRI for PCa characterization. PI-RADS outperformed both TM and SSM DCE-MRI for identification of aggressive cancer. LEVEL OF EVIDENCE: 3 Technical Efficacy: Stage 2 J. MAGN. RESON. IMAGING 2017;46:837-849.

Integrative analysis of diffusion-weighted MRI and genomic data to inform treatment of glioblastoma.

Gene expression profiling from glioblastoma (GBM) patients enables characterization of cancer into subtypes that can be predictive of response to therapy. An integrative analysis of imaging and gene expression data can potentially be used to obtain novel biomarkers that are closely associated with the genetic subtype and gene signatures and thus provide a noninvasive approach to stratify GBM patients. In this retrospective study, we analyzed the expression of 12,042 genes for 558 patients from The Cancer Genome Atlas (TCGA). Among these patients, 50 patients had magnetic resonance imaging (MRI) studies including diffusion weighted (DW) MRI in The Cancer Imaging Archive (TCIA). We identified the contrast enhancing region of the tumors using the pre- and post-contrast T1-weighted MRI images and computed the apparent diffusion coefficient (ADC) histograms from the DW-MRI images. Using the gene expression data, we classified patients into four molecular subtypes, determined the number and composition of genes modules using the gap statistic, and computed gene signature scores. We used logistic regression to find significant predictors of GBM subtypes. We compared the predictors for different subtypes using Mann-Whitney U tests. We assessed detection power using area under the receiver operating characteristic (ROC) analysis. We computed Spearman correlations to determine the associations between ADC and each of the gene signatures. We performed gene enrichment analysis using Ingenuity Pathway Analysis (IPA). We adjusted all p values using the Benjamini and Hochberg method. The mean ADC was a significant predictor for the neural subtype. Neural tumors had a significantly lower mean ADC compared to non-neural tumors ([Formula: see text]), with mean ADC of [Formula: see text] and [Formula: see text] for neural and non-neural tumors, respectively. Mean ADC showed an area under the ROC of 0.75 for detecting neural tumors. We found eight gene modules in the GBM cohort. The mean ADC was significantly correlated with the gene signature related with dendritic cell maturation ([Formula: see text], [Formula: see text]). Mean ADC could be used as a biomarker of a gene signature associated with dendritic cell maturation and to assist in identifying patients with neural GBMs, known to be resistant to aggressive standard of care.

The Impact of Arterial Input Function Determination Variations on Prostate Dynamic Contrast-Enhanced Magnetic Resonance Imaging Pharmacokinetic Modeling: A Multicenter Data Analysis Challenge.

Dynamic contrast-enhanced MRI (DCE-MRI) has been widely used in tumor detection and therapy response evaluation. Pharmacokinetic analysis of DCE-MRI time-course data allows estimation of quantitative imaging biomarkers such as Ktrans(rate constant for plasma/interstitium contrast reagent (CR) transfer) and ve (extravascular and extracellular volume fraction). However, the use of quantitative DCE-MRI in clinical prostate imaging islimited, with uncertainty in arterial input function (AIF, i.e., the time rate of change of the concentration of CR in the blood plasma) determination being one of the primary reasons. In this multicenter data analysis challenge to assess the effects of variations in AIF quantification on estimation of DCE-MRI parameters, prostate DCE-MRI data acquired at one center from 11 prostate cancer patients were shared among nine centers. Each center used its site-specific method to determine the individual AIF from each data set and submitted the results to the managing center. Along with a literature population averaged AIF, these AIFs and their reference-tissue-adjusted variants were used by the managing center to perform pharmacokinetic analysis of the DCE-MRI data sets using the Tofts model (TM). All other variables including tumor region of interest (ROI) definition and pre-contrast T1 were kept the same to evaluate parameter variations caused by AIF variations only. Considerable pharmacokinetic parameter variations were observed with the within-subject coefficient of variation (wCV) of Ktrans obtained with unadjusted AIFs as high as 0.74. AIF-caused variations were larger in Ktrans than ve and both were reduced when reference-tissue-adjusted AIFs were used. The parameter variations were largely systematic, resulting in nearly unchanged parametric map patterns. The CR intravasation rate constant, kep (= Ktrans/ve), was less sensitive to AIF variation than Ktrans (wCV for unadjusted AIFs: 0.45 for kepvs. 0.74 for Ktrans), suggesting that it might be a more robust imaging biomarker of prostate microvasculature than Ktrans.

Neuroendocrine liver metastases: Value of apparent diffusion coefficient and enhancement ratios for characterization of histopathologic grade.

PURPOSE: To assess the value of apparent diffusion coefficient (ADC) measured with diffusion-weighted imaging (DWI) and enhancement ratios (ER) measured with contrast-enhanced T1-weighted imaging (CE-T1WI) for the characterization of histopathologic tumor grade of neuroendocrine tumor liver metastases (NETLM). MATERIALS AND METHODS: Twenty-two patients with pathology-proven NETLM and pretreatment 1.5 Tesla (T) and 3T MRI including DWI were included in this Institutional Review Board-approved retrospective study. ADC histogram parameters, including mean, minimum (min), skewness, and kurtosis as well as ER, were computed for all lesions. Tumor grading was based on the World Health Organization 2010 classification. Kruskal-Wallis and Mann-Whitney test were used to assess for differences in ADC and ER between different tumor grades. MRI parameters were correlated with pathologic findings using Spearman correlation test. Receiver operating characteristic analysis was performed to determine optimum thresholds for predicting tumor grade. RESULTS: Forty-eight NETLM (mean size 3.5 cm) were analyzed with the following grade distribution: G1 (n = 25), G2 (n = 16), and G3 (n = 7). ADC-mean (×10-3 mm2 /s) of G3 tumors (0.87 ± 0.43) was significantly lower than that of G1 (1.47 ± 0.63) and G2 (1.27 ± 0.63; P = 0.042). A weak significant negative correlation was observed between ADC and tumor grade (ADC-mean: r = -0.33, P = 0.02; ADC-min: r = -0.37, P = 0.01) and Ki-67 (ADC-mean: r = -0.31, P = 0.03; ADC-min: r = -0.39, P = 0.007). AUROC, sensitivity and specificity of ADC-mean/ADC-min/ER (measured at the early arterial phase) for differentiation of G3 versus G1-G2 were 0.80/0.76/0.67, 100%/50%/70%, and 68.4%/84.2%/66.6%, respectively. CONCLUSION: ADC is a promising marker for characterization of histopathologic grade of NETLM. These results should be confirmed in a prospective study. J. Magn. Reson. Imaging 2016;44:1432-1441.

Hepatocellular carcinoma: IVIM diffusion quantification for prediction of tumor necrosis compared to enhancement ratios.

PURPOSE: To correlate intra voxel incoherent motion (IVIM) diffusion parameters of liver parenchyma and hepatocellular carcinoma (HCC) with degree of liver/tumor enhancement and necrosis; and to assess the diagnostic performance of diffusion parameters vs. enhancement ratios (ER) for prediction of complete tumor necrosis. PATIENTS AND METHODS: In this IRB approved HIPAA compliant study, we included 46 patients with HCC who underwent IVIM diffusion-weighted (DW) MRI in addition to routine sequences at 3.0 T. True diffusion coefficient (D), pseudo-diffusion coefficient (D*), perfusion fraction (PF) and apparent diffusion coefficient (ADC) were quantified in tumors and liver parenchyma. Tumor ER were calculated using contrast-enhanced imaging, and degree of tumor necrosis was assessed using post-contrast image subtraction. IVIM parameters and ER were compared between HCC and background liver and between necrotic and viable tumor components. ROC analysis for prediction of complete tumor necrosis was performed. RESULTS: 79 HCCs were assessed (mean size 2.5 cm). D, PF and ADC were significantly higher in HCC vs. liver (p < 0.0001). There were weak significant negative/positive correlations between D/PF and ER, and significant correlations between D/PF/ADC and tumor necrosis (for D, r 0.452, p < 0.001). Among diffusion parameters, D had the highest area under the curve (AUC 0.811) for predicting complete tumor necrosis. ER outperformed diffusion parameters for prediction of complete tumor necrosis (AUC > 0.95, p < 0.002). CONCLUSION: D has a reasonable diagnostic performance for predicting complete tumor necrosis, however lower than that of contrast-enhanced imaging.

Comparison of gadoxetic acid to gadobenate dimeglumine for assessment of biliary anatomy of potential liver donors.

PURPOSE: To compare MRI using gadobenate dimeglumine (Gd-BOPTA) vs. gadoxetic acid disodium (Gd-EOB-DTPA) for the assessment of biliary anatomy of potential liver donors. METHODS: 76 potential liver donors (39 M/37 F, mean 38 years) who underwent 1.5T MRI using Gd-BOPTA (n = 37) or Gd-EOB-DTPA (n = 39) were retrospectively evaluated. T2 cholangiogram (T2 MRC) and delayed hepatobiliary phase (HBP) T1 cholangiogram (T1 MRC) (performed during HBP 20 min after injection of Gd-EOB-DTPA and 1-2 h after Gd-BOPTA injection) were obtained in addition to MR angiogram/venogram. Two independent observers evaluated image quality (IQ) and conspicuity scores (CS) of the biliary system. Biliary anatomy was assessed in 3 reading sessions (T2 MRC, T1 MRC, and combined T2/T1 MRC). Reference standard consisted of consensus reading of two separate observers of all image sets, clinical/surgical information and intraoperative cholangiogram when available. Datasets were compared using the Mann-Whitney U test or Chi-squared test. RESULTS: There was no difference in IQ for T1 MRC using either contrast agent or T2 MRC vs. T1 MRC for both observers (all p values >0.07). There was superior CS for T2 MRC vs. Gd-BOPTA T1 MRC for both observers and T2 MRC vs. Gd-EOB for one observer (p < 0.001). No difference was found for biliary variant detection for T1 MRC (with either contrast agent) vs. T2 MRC. Combined T2/T1 MRC demonstrated improved sensitivity for biliary variant detection using Gd-BOPTA for both observers (p < 0.004) and Gd-EOB-DTPA for one observer (p < 0.001). CONCLUSION: Equivalent image quality was found for T1 MRC obtained with Gd-BOPTA or Gd-EOB-DTPA and T2 MRC. T1 MRC is equivalent to T2 MRC for detection of variant biliary anatomy, and the combination of sequences may have added value.

Non-invasive prediction of portal pressures using CT and MRI in chronic liver disease.

PURPOSE: To assess the diagnostic value of a fast scoring system based on non-invasive cross-sectional imaging to predict portal hypertension (PH) in patients with liver disease. METHODS: In this retrospective study, we included patients who underwent contrast-enhanced CT or MRI within 3 months of hepatic venous pressure gradient (HVPG) measurements. Two independent observers provided an imaging-based scoring system (max of 9): number of variceal sites, volume of ascites, and spleen size. ROC analysis was performed to predict the presence of PH (HVPG ≥ 5 mmHg) and clinically significant PH (HVPG ≥ 10 mmHg). RESULTS: Our cohort consists of 143 patients with mean HVPG of 13.1 ± 2.0 mmHg. Mean PH scores from the two observers were 3.9 ± 2.7 and 3.2 ± 2.5. There was a significant correlation between PH score and HVPG (r = 0.58, p < 0.001 for both observers) with high inter-observer agreement (kappa 0.71). AUCs of 0.78-0.76 and 0.83-0.81 were observed for diagnosing HVPG ≥ 5 mmHg and HVPG ≥ 10 mmHg, respectively, for observers 1 and 2. CONCLUSIONS: We have developed a fast PH imaging-based composite score, which could be used for non-invasive detection of clinically significant PH.

Prospective comparison of magnetic resonance imaging to transient elastography and serum markers for liver fibrosis detection.

BACKGROUND & AIMS: Establishing accurate non-invasive methods of liver fibrosis quantification remains a major unmet need. Here, we assessed the diagnostic value of a multiparametric magnetic resonance imaging (MRI) protocol including diffusion-weighted imaging (DWI), dynamic contrast-enhanced (DCE)-MRI and magnetic resonance elastography (MRE) in comparison with transient elastography (TE) and blood tests [including ELF (Enhanced Liver Fibrosis) and APRI] for liver fibrosis detection. METHODS: In this single centre cross-sectional study, we prospectively enrolled 60 subjects with liver disease who underwent multiparametric MRI (DWI, DCE-MRI and MRE), TE and blood tests. Correlation was assessed between non-invasive modalities and histopathologic findings including stage, grade and collagen content, while accounting for covariates such as age, sex, BMI, HCV status and MRI-derived fat and iron content. ROC curve analysis evaluated the performance of each technique for detection of moderate-to-advanced liver fibrosis (F2-F4) and advanced fibrosis (F3-F4). RESULTS: Magnetic resonance elastography provided the strongest correlation with fibrosis stage (r = 0.66, P < 0.001), inflammation grade (r = 0.52, P < 0.001) and collagen content (r = 0.53, P = 0.036). For detection of moderate-to-advanced fibrosis (F2-F4), AUCs were 0.78, 0.82, 0.72, 0.79, 0.71 for MRE, TE, DCE-MRI, DWI and APRI, respectively. For detection of advanced fibrosis (F3-F4), AUCs were 0.94, 0.77, 0.79, 0.79 and 0.70, respectively. CONCLUSIONS: Magnetic resonance elastography provides the highest correlation with histopathologic markers and yields high diagnostic performance for detection of advanced liver fibrosis and cirrhosis, compared to DWI, DCE-MRI, TE and serum markers.

DCE-MRI of hepatocellular carcinoma: perfusion quantification with Tofts model versus shutter-speed model--initial experience.

OBJECTIVE: To quantify hepatocellular carcinoma (HCC) perfusion and flow with the fast exchange regime-allowed Shutter-Speed model (SSM) compared to the Tofts model (TM). MATERIALS AND METHODS: In this prospective study, 25 patients with HCC underwent DCE-MRI. ROIs were placed in liver parenchyma, portal vein, aorta and HCC lesions. Signal intensities were analyzed employing dual-input TM and SSM models. ART (arterial fraction), K (trans) (contrast agent transfer rate constant from plasma to extravascular extracellular space), ve (extravascular extracellular volume fraction), kep (contrast agent intravasation rate constant), and τi (mean intracellular water molecule lifetime) were compared between liver parenchyma and HCC, and ART, K (trans), v e and k ep were compared between models using Wilcoxon tests and limits of agreement. Test-retest reproducibility was assessed in 10 patients. RESULTS: ART and v e obtained with TM; ART, ve, ke and τi obtained with SSM were significantly different between liver parenchyma and HCC (p < 0.04). Parameters showed variable reproducibility (CV range 14.7-66.5% for both models). Liver K (trans) and ve; HCC ve and kep were significantly different when estimated with the two models (p < 0.03). CONCLUSION: Our results show differences when computed between the TM and the SSM. However, these differences are smaller than parameter reproducibilities and may be of limited clinical significance.

Demonstration of nonlinearity bias in the measurement of the apparent diffusion coefficient in multicenter trials.

PURPOSE: Characterize system-specific bias across common magnetic resonance imaging (MRI) platforms for quantitative diffusion measurements in multicenter trials. METHODS: Diffusion weighted imaging (DWI) was performed on an ice-water phantom along the superior-inferior (SI) and right-left (RL) orientations spanning ± 150 mm. The same scanning protocol was implemented on 14 MRI systems at seven imaging centers. The bias was estimated as a deviation of measured from known apparent diffusion coefficient (ADC) along individual DWI directions. The relative contributions of gradient nonlinearity, shim errors, imaging gradients, and eddy currents were assessed independently. The observed bias errors were compared with numerical models. RESULTS: The measured systematic ADC errors scaled quadratically with offset from isocenter, and ranged between -55% (SI) and 25% (RL). Nonlinearity bias was dependent on system design and diffusion gradient direction. Consistent with numerical models, minor ADC errors (± 5%) due to shim, imaging and eddy currents were mitigated by double echo DWI and image coregistration of individual gradient directions. CONCLUSION: The analysis confirms gradient nonlinearity as a major source of spatial DW bias and variability in off-center ADC measurements across MRI platforms, with minor contributions from shim, imaging gradients and eddy currents. The developed protocol enables empiric description of systematic bias in multicenter quantitative DWI studies.

3D T1 relaxometry pre and post gadoxetic acid injection for the assessment of liver cirrhosis and liver function.

PURPOSE: To assess the diagnostic value of a 3D dual-flip-angle (DFA) T1 mapping technique with whole liver coverage before and after gadoxetic acid injection for assessment of cirrhosis and liver function, compared to blood tests (APRI: aspartate aminotransferase-to-platelet ratio index). MATERIALS AND METHODS: A total of 133 patients who underwent gadoxetic acid-enhanced liver MRI including a 3D FLASH DFA-T1 mapping sequence before and 20min post-contrast (hepatobiliary phase, HBP) were included in this retrospective IRB approved study. T1 values (msec) were measured on pre-contrast and during HBP in liver parenchyma, ΔT1 (%) was calculated as [(T1 pre-T1 post)/T1 pre]×100. T1 and ΔT1 values were compared between cirrhotic and non-cirrhotic patients and between patients stratified using Child-Pugh and Model for End-Stage Liver Disease (MELD) scores using Mann-Whitney U test. Diagnostic performance of T1 mapping parameters vs. APRI for diagnosing cirrhosis and for assessing degree of liver dysfunction was evaluated using ROC analysis. RESULTS: Fifty non-cirrhotic and 83 cirrhotic patients [Child-Pugh A (n=41), B (n=31) and C (n=11)] were included. There was no significant difference in pre-contrast T1 values between cirrhotic and non-cirrhotic patients. T1-HBP and ΔT1 values were significantly different in patients with cirrhosis (p<0.0001) and higher MELD scores (>17) (p=0.003). ΔT1 showed significant strong correlations with Child-Pugh and MELD scores (r=-0.7, p<0.0001; r=-0.56, p<0.001 respectively). Similar AUCs (p=0.9) for detection of liver cirrhosis were observed for T1 HBP (0.83), ΔT1 (0.86) and APRI (0.85); however APRI showed limited sensitivity (≤55%) in comparison with ΔT1 (74.7%) and T1 HBP (80.7%). CONCLUSION: 3D DFA-T1 mapping sequence used before and after gadoxetic acid injection is useful for the diagnosis of cirrhosis and for the assessment of liver function.

Simultaneous measurement of hepatic and splenic stiffness using MR elastography: preliminary experience.

PURPOSE: To compare MR elastography (MRE) using a single and a dual driver excitation for the quantification of hepatic and splenic stiffness (HS and SS), and to investigate the performance of HS and SS measured with single or dual driver excitation for the detection of liver cirrhosis in subjects with liver disease. PATIENTS AND METHODS: This prospective HIPAA compliant and IRB approved study involved 49 subjects who underwent MRE at 3.0T, comparing three different acquisition methods (single driver on the liver, single driver on the spleen and dual driver acoustic excitation). A Mann-Whitney test was used to assess changes in stiffness values. Bland-Altman analysis was used to compare single and dual driver configurations for each organ. Performance for detection of liver cirrhosis was assessed using ROC analysis. Pearson correlation was used to estimate the dependence of HS and SS on spleen size. RESULTS: There were 40 noncirrhotic and 9 cirrhotic patients. There was good agreement between stiffness values measured with a single or a dual driver (Bland-Altman limits of agreement -14.3 % to 18.9 % and -18.1 % to 29.7 %, CV 6.4 % and 9.4 %, for HS and SS. respectively). HS and SS were higher in subjects with liver cirrhosis (p < 0.001), with excellent detection performance (AUROC range 0.87-0.93). SS correlated strongly with spleen size (r = 0.69, p < 0.001), while HS showed weak correlation (r = 0.38, p = 0.006). CONCLUSION: Using a dual acoustic driver configuration, hepatic and splenic stiffness can be simultaneously estimated with good concordance with single driver measurement.

Abdominal 4D flow MR imaging in a breath hold: combination of spiral sampling and dynamic compressed sensing for highly accelerated acquisition.

PURPOSE: To develop a highly accelerated phase-contrast cardiac-gated volume flow measurement (four-dimensional [4D] flow) magnetic resonance (MR) imaging technique based on spiral sampling and dynamic compressed sensing and to compare this technique with established phase-contrast imaging techniques for the quantification of blood flow in abdominal vessels. MATERIALS AND METHODS: This single-center prospective study was compliant with HIPAA and approved by the institutional review board. Ten subjects (nine men, one woman; mean age, 51 years; age range, 30-70 years) were enrolled. Seven patients had liver disease. Written informed consent was obtained from all participants. Two 4D flow acquisitions were performed in each subject, one with use of Cartesian sampling with respiratory tracking and the other with use of spiral sampling and a breath hold. Cartesian two-dimensional (2D) cine phase-contrast images were also acquired in the portal vein. Two observers independently assessed vessel conspicuity on phase-contrast three-dimensional angiograms. Quantitative flow parameters were measured by two independent observers in major abdominal vessels. Intertechnique concordance was quantified by using Bland-Altman and logistic regression analyses. RESULTS: There was moderate to substantial agreement in vessel conspicuity between 4D flow acquisitions in arteries and veins (κ = 0.71 and 0.61, respectively, for observer 1; κ = 0.71 and 0.44 for observer 2), whereas more artifacts were observed with spiral 4D flow (κ = 0.30 and 0.20). Quantitative measurements in abdominal vessels showed good equivalence between spiral and Cartesian 4D flow techniques (lower bound of the 95% confidence interval: 63%, 77%, 60%, and 64% for flow, area, average velocity, and peak velocity, respectively). For portal venous flow, spiral 4D flow was in better agreement with 2D cine phase-contrast flow (95% limits of agreement: -8.8 and 9.3 mL/sec, respectively) than was Cartesian 4D flow (95% limits of agreement: -10.6 and 14.6 mL/sec). CONCLUSION: The combination of highly efficient spiral sampling with dynamic compressed sensing results in major acceleration for 4D flow MR imaging, which allows comprehensive assessment of abdominal vessel hemodynamics in a single breath hold.

DCE-MRI of the liver: reconstruction of the arterial input function using a low dose pre-bolus contrast injection.

PURPOSE: To assess the quality of the arterial input function (AIF) reconstructed using a dedicated pre-bolus low-dose contrast material injection imaged with a high temporal resolution and the resulting estimated liver perfusion parameters. MATERIALS AND METHODS: In this IRB-approved prospective study, 24 DCE-MRI examinations were performed in 21 patients with liver disease (M/F 17/4, mean age 56 y). The examination consisted of 1.3 mL and 0.05 mmol/kg of gadobenate dimeglumine for pre-bolus and main bolus acquisitions, respectively. The concentration-curve of the abdominal aorta in the pre-bolus acquisition was used to reconstruct the AIF. AIF quality and shape parameters obtained with pre-bolus and main bolus acquisitions and the resulting estimated hepatic perfusion parameters obtained with a dual-input single compartment model were compared between the 2 methods. Test-retest reproducibility of perfusion parameters were assessed in three patients. RESULTS: The quality of the pre-bolus AIF curve was significantly better than that of main bolus AIF. Shape parameters peak concentration, area under the time activity curve of gadolinium contrast at 60 s and upslope of pre-bolus AIF were all significantly higher, while full width at half maximum was significantly lower than shape parameters of main bolus AIF. Improved liver perfusion parameter reproducibility was observed using pre-bolus acquisition [coefficient of variation (CV) of 4.2%-38.7% for pre-bolus vs. 12.1-71.4% for main bolus] with the exception of distribution volume (CV of 23.6% for pre-bolus vs. 15.8% for main bolus). The CVs between pre-bolus and main bolus for the perfusion parameters were lower than 14%. CONCLUSION: The AIF reconstructed with pre-bolus low dose contrast injection displays better quality and shape parameters and enables improved liver perfusion parameter reproducibility, although the resulting liver perfusion parameters demonstrated no clinically significant differences between pre-bolus and main bolus acquisitions.

Intravoxel incoherent motion diffusion imaging of the liver: optimal b-value subsampling and impact on parameter precision and reproducibility.

PURPOSE: To increase diffusion sampling efficiency in intravoxel incoherent motion (IVIM) diffusion-weighted imaging (DWI) of the liver by reducing the number of diffusion weightings (b-values). MATERIALS AND METHODS: In this IRB approved HIPAA compliant prospective study, 53 subjects (M/F 38/15, mean age 52 ± 13 y) underwent IVIM DWI at 1.5T using 16 b-values (0-800s/mm(2)), with 14 subjects having repeat exams to assess IVIM parameter reproducibility. A biexponential diffusion model was used to quantify IVIM hepatic parameters (PF: perfusion fraction, D: true diffusion and D*: pseudo diffusion). All possible subsets of the 16 b-values were probed, with number of b values ranging from 4 to 15, and corresponding parameters were quantified for each subset. For each b-value subset, global parameter estimation error was computed against the parameters obtained with all 16 b-values and the subsets providing the lowest error were selected. Interscan estimation error was also evaluated between repeat exams to assess reproducibility of the IVIM technique in the liver. The optimal b-values distribution was selected such that the number of b-values was minimal while keeping parameter estimation error below interscan reproducibility error. RESULTS: As the number of b-values decreased, the estimation error increased for all parameters, reflecting decreased precision of IVIM metrics. Using an optimal set of 4 b-values (0, 15, 150 and 800s/mm(2)), the errors were 6.5, 22.8 and 66.1% for D, PF and D* respectively. These values lie within the range of test-retest reproducibility for the corresponding parameters, with errors of 12.0, 32.3 and 193.8% for D, PF and D* respectively. CONCLUSION: A set of 4 optimized b-values can be used to estimate IVIM parameters in the liver with significantly shorter acquisition time (up to 75%), without substantial degradation of IVIM parameter precision and reproducibility compared to the 16 b-value acquisition used as the reference.

DCE-MRI of the liver: effect of linear and nonlinear conversions on hepatic perfusion quantification and reproducibility.

PURPOSE: To evaluate the effect of different methods to convert magnetic resonance (MR) signal intensity (SI) to gadolinium concentration ([Gd]) on estimation and reproducibility of model-free and modeled hepatic perfusion parameters measured with dynamic contrast-enhanced (DCE)-MRI. MATERIALS AND METHODS: In this Institutional Review Board (IRB)-approved prospective study, 23 DCE-MRI examinations of the liver were performed on 17 patients. SI was converted to [Gd] using linearity vs. nonlinearity assumptions (using spoiled gradient recalled echo [SPGR] signal equations). The [Gd] vs. time curves were analyzed using model-free parameters and a dual-input single compartment model. Perfusion parameters obtained with the two conversion methods were compared using paired Wilcoxon test. Test-retest and interobserver reproducibility of perfusion parameters were assessed in six patients. RESULTS: There were significant differences between the two conversion methods for the following parameters: AUC60 (area under the curve at 60 s, P < 0.001), peak gadolinium concentration (Cpeak, P < 0.001), upslope (P < 0.001), Fp (portal flow, P = 0.04), total hepatic flow (Ft, P = 0.007), and MTT (mean transit time, P < 0.001). Our preliminary results showed acceptable to good reproducibility for all model-free parameters for both methods (mean coefficient of variation [CV] range, 11.87-23.7%), except for upslope (CV = 37%). Among modeled parameters, DV (distribution volume) had CV <22% with both methods, PV and MTT showed CV <21% and <29% using SPGR equations, respectively. Other modeled parameters had CV >30% with both methods. CONCLUSION: Linearity assumption is acceptable for quantification of model-free hepatic perfusion parameters while the use of SPGR equations and T1 mapping may be recommended for the quantification of modeled hepatic perfusion parameters.

Quantitative liver MRI combining phase contrast imaging, elastography, and DWI: assessment of reproducibility and postprandial effect at 3.0 T.

PURPOSE: To quantify short-term reproducibility (in fasting conditions) and postprandial changes after a meal in portal vein (PV) flow parameters measured with phase contrast (PC) imaging, liver diffusion parameters measured with multiple b value diffusion-weighted imaging (DWI) and liver stiffness (LS) measured with MR elastography (MRE) in healthy volunteers and patients with liver disease at 3.0 T. MATERIALS AND METHODS: In this IRB-approved prospective study, 30 subjects (11 healthy volunteers and 19 liver disease patients; 23 males, 7 females; mean age 46.5 y) were enrolled. Imaging included 2D PC imaging, multiple b value DWI and MRE. Subjects were initially scanned twice in fasting state to assess short-term parameter reproducibility, and then scanned 20 min. after a liquid meal. PV flow/velocity, LS, liver true diffusion coefficient (D), pseudodiffusion coefficient (D*), perfusion fraction (PF) and apparent diffusion coefficient (ADC) were measured in fasting and postprandial conditions. Short-term reproducibility was assessed in fasting conditions by measuring coefficients of variation (CV) and Bland-Altman limits of agreement. Differences in MR metrics before and after caloric intake and between healthy volunteers and liver disease patients were assessed. RESULTS: PV flow parameters, D, ADC and LS showed good to excellent short-term reproducibility in fasting state (CV <16%), while PF and D* showed acceptable and poor reproducibility (CV = 20.4% and 51.6%, respectively). PV flow parameters and LS were significantly higher (p<0.04) in postprandial state while liver diffusion parameters showed no significant change (p>0.2). LS was significantly higher in liver disease patients compared to healthy volunteers both in fasting and postprandial conditions (p<0.001). Changes in LS were significantly correlated with changes in PV flow (Spearman rho = 0.48, p = 0.013). CONCLUSIONS: Caloric intake had no/minimal/large impact on diffusion/stiffness/portal vein flow, respectively. PC MRI and MRE but not DWI should be performed in controlled fasting state.

Variations of dynamic contrast-enhanced magnetic resonance imaging in evaluation of breast cancer therapy response: a multicenter data analysis challenge.

Pharmacokinetic analysis of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) time-course data allows estimation of quantitative parameters such as K (trans) (rate constant for plasma/interstitium contrast agent transfer), v e (extravascular extracellular volume fraction), and v p (plasma volume fraction). A plethora of factors in DCE-MRI data acquisition and analysis can affect accuracy and precision of these parameters and, consequently, the utility of quantitative DCE-MRI for assessing therapy response. In this multicenter data analysis challenge, DCE-MRI data acquired at one center from 10 patients with breast cancer before and after the first cycle of neoadjuvant chemotherapy were shared and processed with 12 software tools based on the Tofts model (TM), extended TM, and Shutter-Speed model. Inputs of tumor region of interest definition, pre-contrast T1, and arterial input function were controlled to focus on the variations in parameter value and response prediction capability caused by differences in models and associated algorithms. Considerable parameter variations were observed with the within-subject coefficient of variation (wCV) values for K (trans) and v p being as high as 0.59 and 0.82, respectively. Parameter agreement improved when only algorithms based on the same model were compared, e.g., the K (trans) intraclass correlation coefficient increased to as high as 0.84. Agreement in parameter percentage change was much better than that in absolute parameter value, e.g., the pairwise concordance correlation coefficient improved from 0.047 (for K (trans)) to 0.92 (for K (trans) percentage change) in comparing two TM algorithms. Nearly all algorithms provided good to excellent (univariate logistic regression c-statistic value ranging from 0.8 to 1.0) early prediction of therapy response using the metrics of mean tumor K (trans) and k ep (=K (trans)/v e, intravasation rate constant) after the first therapy cycle and the corresponding percentage changes. The results suggest that the interalgorithm parameter variations are largely systematic, which are not likely to significantly affect the utility of DCE-MRI for assessment of therapy response.

Maximum-parsimony haplotype frequencies inference based on a joint constrained sparse representation of pooled DNA.

BACKGROUND: DNA pooling constitutes a cost effective alternative in genome wide association studies. In DNA pooling, equimolar amounts of DNA from different individuals are mixed into one sample and the frequency of each allele in each position is observed in a single genotype experiment. The identification of haplotype frequencies from pooled data in addition to single locus analysis is of separate interest within these studies as haplotypes could increase statistical power and provide additional insight. RESULTS: We developed a method for maximum-parsimony haplotype frequency estimation from pooled DNA data based on the sparse representation of the DNA pools in a dictionary of haplotypes. Extensions to scenarios where data is noisy or even missing are also presented. The resulting method is first applied to simulated data based on the haplotypes and their associated frequencies of the AGT gene. We further evaluate our methodology on datasets consisting of SNPs from the first 7Mb of the HapMap CEU population. Noise and missing data were further introduced in the datasets in order to test the extensions of the proposed method. Both HIPPO and HAPLOPOOL were also applied to these datasets to compare performances. CONCLUSIONS: We evaluate our methodology on scenarios where pooling is more efficient relative to individual genotyping; that is, in datasets that contain pools with a small number of individuals. We show that in such scenarios our methodology outperforms state-of-the-art methods such as HIPPO and HAPLOPOOL.