DCE-MRI model selection for investigating disruption of microvascular function in livers with metastatic disease

Dynamic contrast enhanced MRI demonstrate altered placental perfusion in the STOX1A preeclampsia mouse model

INTRODUCTION

Preeclampsia, a hypertensive disorder of pregnancy triggered by placental dysfunction, is reproduced in the murine STOX1A model, with hypertension, proteinuria, and abnormalities in umbilical and uterine Dopplers. Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is an innovative technique that provides insights into tissue perfusion. The present study aims at analyzing placental perfusion using DCE-MRI to further characterize placental defects in the STOX1A model.

METHODS

Two study groups were formed: the "TgSTOX13 pregnancy group" from mating TgSTOX13 genotype males with wild-type females, and the "wild-type pregnancy group" from mating wild-type males with wild-type females. Blood pressure, urinary albumin to creatinine ratio, and fetal weights were measured and compared between the groups, while perfusion parameters were analyzed using both conventional compartmental (1C) and free-time point-Hermite (FTPH) models in the DCE analysis.

RESULTS

Seventeen pregnant mice in the "TgSTOX13 pregnancy group" and thirteen in the "wild-type pregnant group" were included in the analysis. During late gestation, the TgSTOX13 pregnancy group exhibited higher blood pressure, elevated albumin/creatinine ratio, and decreased fetal weights compared to the wild-type pregnancy group. In the DCE analysis utilizing the 1C model, blood flow (Fb) was significantly reduced by approximately 31.8 % in the TgSTOX13 pregnancy group compared to the wild-type pregnancy group (p < 0.01), a finding corroborated by the FTPH model with a reduction estimated at 31.5 % (p < 0.01).

DISCUSSION

Our investigation successfully utilized DCE MRI to assess placental perfusion in a mouse model of preeclampsia, revealing a significant reduction of approximately 30 % in the preeclamptic mice, mirroring human pathophysiology.

Identifiability of spatiotemporal tissue perfusion models

Objective.Standard models for perfusion quantification in DCE-MRI produce a bias by treating voxels as isolated systems. Spatiotemporal models can remove this bias, but it is unknown whether they are fundamentally identifiable. The aim of this study is to investigate this question in silico using one-dimensional toy systems with a one-compartment blood flow model and a two-compartment perfusion model.Approach.For each of the two models, identifiability is explored theoretically and in-silico for three systems. Concentrations over space and time are simulated by forward propagation. Different levels of noise and temporal undersampling are added to investigate sensitivity to measurement error. Model parameters are fitted using a standard gradient descent algorithm, applied iteratively with a stepwise increasing time window. Model fitting is repeated with different initial values to probe uniqueness of the solution. Reconstruction accuracy is quantified for each parameter by comparison to the ground truth.Main results.Theoretical analysis shows that flows and volume fractions are only identifiable up to a constant, and that this degeneracy can be removed by proper choice of parameters. Simulations show that in all cases, the tissue concentrations can be reconstructed accurately. The one-compartment model shows accurate reconstruction of blood velocities and arterial input functions, independent of the initial values and robust to measurement error. The two-compartmental perfusion model was not fully identifiable, showing good reconstruction of arterial velocities and input functions, but multiple valid solutions for the perfusion parameters and venous velocities, and a strong sensitivity to measurement error in these parameters.Significance.These results support the use of one-compartment spatiotemporal flow models, but two-compartment perfusion models were not sufficiently identifiable. Future studies should investigate whether this degeneracy is resolved in more realistic 2D and 3D systems, by adding physically justified constraints, or by optimizing experimental parameters such as injection duration or temporal resolution.

A Quantitative Multiparametric MRI Analysis Platform for Estimation of Robust Imaging Biomarkers in Clinical Oncology

There is a need to develop user-friendly imaging tools estimating robust quantitative biomarkers (QIBs) from multiparametric (mp)MRI for clinical applications in oncology. Quantitative metrics derived from (mp)MRI can monitor and predict early responses to treatment, often prior to anatomical changes. We have developed a vendor-agnostic, flexible, and user-friendly MATLAB-based toolkit, MRI-Quantitative Analysis and Multiparametric Evaluation Routines ("MRI-QAMPER", current release v3.0), for the estimation of quantitative metrics from dynamic contrast-enhanced (DCE) and multi-b value diffusion-weighted (DW) MR and MR relaxometry. MRI-QAMPER's functionality includes generating numerical parametric maps from these methods reflecting tumor permeability, cellularity, and tissue morphology. MRI-QAMPER routines were validated using digital reference objects (DROs) for DCE and DW MRI, serving as initial approval stages in the National Cancer Institute Quantitative Imaging Network (NCI/QIN) software benchmark. MRI-QAMPER has participated in DCE and DW MRI Collaborative Challenge Projects (CCPs), which are key technical stages in the NCI/QIN benchmark. In a DCE CCP, QAMPER presented the best repeatability coefficient (RC = 0.56) across test-retest brain metastasis data, out of ten participating DCE software packages. In a DW CCP, QAMPER ranked among the top five (out of fourteen) tools with the highest area under the curve (AUC) for prostate cancer detection. This platform can seamlessly process mpMRI data from brain, head and neck, thyroid, prostate, pancreas, and bladder cancer. MRI-QAMPER prospectively analyzes dose de-escalation trial data for oropharyngeal cancer, which has earned it advanced NCI/QIN approval for expanded usage and applications in wider clinical trials.

Cell type-independent profiling of interactions between intracellular pathogens and the human phosphoproteome

Interactions between proteins from intracellular pathogens and host proteins in an infected cell are often mediated by post-translational modifications encoded in the host proteome. Identifying protein modifications, such as phosphorylation, that dictate these interactions remains a defining challenge in unraveling the molecular mechanisms of pathogenesis. We have developed a platform in engineered bacteria that displays over 110,000 phosphorylated human proteins coupled to a fluorescent reporter system capable of identifying the host-pathogen interactome of phosphoproteins (H-PIP). This resource broadly enables cell-type independent interrogation and discovery of proteins from intracellular pathogens capable of binding phosphorylated human proteins. As an example of the H-PIP platform, we generated a unique, high-resolution SARS-CoV-2 interaction network which expanded our knowledge of viral protein function and identified understudied areas of host pathology.

Multimodality Imaging Assessment of Desmoid Tumors: The Great Mime in the Era of Multidisciplinary Teams

Desmoid tumors (DTs), also known as desmoid fibromatosis or aggressive fibromatosis, are rare, locally invasive, non-metastatic soft tissue tumors. Although histological results represent the gold standard diagnosis, imaging represents the fundamental tool for the diagnosis of these tumors. Although histological analysis represents the gold standard for diagnosis, imaging represents the fundamental tool for the diagnosis of these tumors. DTs represent a challenge for the radiologist, being able to mimic different pathological conditions. A proper diagnosis is required to establish an adequate therapeutic approach. Multimodality imaging, including ultrasound (US), computed tomography (CT) and Magnetic Resonance Imaging (MRI), should be preferred. Different imaging techniques can also guide minimally invasive treatments and monitor their effectiveness. The purpose of this review is to describe the state-of-the-art multidisciplinary imaging of DTs; and its role in patient management.

Conventional, functional and radiomics assessment for intrahepatic cholangiocarcinoma

Background

This paper offers an assessment of diagnostic tools in the evaluation of Intrahepatic Cholangiocarcinoma (ICC).

Methods

Several electronic datasets were analysed to search papers on morphological and functional evaluation in ICC patients. Papers published in English language has been scheduled from January 2010 to December 2021.

Results

We found that 88 clinical studies satisfied our research criteria. Several functional parameters and morphological elements allow a truthful ICC diagnosis. The contrast medium evaluation, during the different phases of contrast studies, support the recognition of several distinctive features of ICC. The imaging tool to employed and the type of contrast medium in magnetic resonance imaging, extracellular or hepatobiliary, should change considering patient, departement, and regional features. Also, Radiomics is an emerging area in the evaluation of ICCs.

Post treatment studies are required to evaluate the efficacy and the safety of therapies so as the patient surveillance.

Conclusions

Several morphological and functional data obtained during Imaging studies allow a truthful ICC diagnosis.

A model selection framework to quantify microvascular liver function in gadoxetate-enhanced MRI: Application to healthy liver, diseased tissue, and hepatocellular carcinoma

PURPOSE

We introduce a novel, generalized tracer kinetic model selection framework to quantify microvascular characteristics of liver and tumor tissue in gadoxetate-enhanced dynamic contrast-enhanced MRI (DCE-MRI).

METHODS

Our framework includes a hierarchy of nested models, from which physiological parameters are derived in 2 regimes, corresponding to the active transport and free diffusion of gadoxetate. We use simulations to show the sensitivity of model selection and parameter estimation to temporal resolution, time-series duration, and noise. We apply the framework in 8 healthy volunteers (time-series duration up to 24 minutes) and 10 patients with hepatocellular carcinoma (6 minutes).

RESULTS

The active transport regime is preferred in 98.6% of voxels in volunteers, 82.1% of patients' non-tumorous liver, and 32.2% of tumor voxels. Interpatient variations correspond to known co-morbidities. Simulations suggest both datasets have sufficient temporal resolution and signal-to-noise ratio, while patient data would be improved by using a time-series duration of at least 12 minutes.

CONCLUSIONS

In patient data, gadoxetate exhibits different kinetics: (a) between liver and tumor regions and (b) within regions due to liver disease and/or tumor heterogeneity. Our generalized framework selects a physiological interpretation at each voxel, without preselecting a model for each region or duplicating time-consuming optimizations for models with identical functional forms.

Diagnostic evaluation and ablation treatments assessment in hepatocellular carcinoma

This article provides an overview of diagnostic evaluation and ablation treatment assessment in Hepatocellular Carcinoma (HCC). Only studies, in the English language from January 2010 to January 202, evaluating the diagnostic tools and assessment of ablative therapies in HCC patients were included. We found 173 clinical studies that satisfied the inclusion criteria.HCC may be noninvasively diagnosed by imaging findings. Multiphase contrast-enhanced imaging is necessary to assess HCC. Intravenous extracellular contrast agents are used for CT, while the agents used for MRI may be extracellular or hepatobiliary. Both gadoxetate disodium and gadobenate dimeglumine may be used in hepatobiliary phase imaging. For treatment-naive patients undergoing CT, unenhanced imaging is optional; however, it is required in the post treatment setting for CT and all MRI studies. Late arterial phase is strongly preferred over early arterial phase. The choice of modality (CT, US/CEUS or MRI) and MRI contrast agent (extracelllar or hepatobiliary) depends on patient, institutional, and regional factors. MRI allows to link morfological and functional data in the HCC evaluation. Also, Radiomics is an emerging field in the assessment of HCC patients.Postablation imaging is necessary to assess the treatment results, to monitor evolution of the ablated tissue over time, and to evaluate for complications. Post- thermal treatments, imaging should be performed at regularly scheduled intervals to assess treatment response and to evaluate for new lesions and potential complications.

Correlation of radiomic features on dynamic contrast-enhanced magnetic resonance with microvessel density in hepatocellular carcinoma based on different models

Objective

To explore the correlations of radiomic features of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) with microvessel density (MVD) in patients with hepatocellular carcinoma (HCC), based on single-input and dual-input two-compartment extended Tofts (SITET and DITET) models.

Methods

We compared the quantitative parameters of SITET and DITET models for DCE-MRI in 30 patients with HCC using paired sample t-tests. The correlations of SITET and DITET model parameters with CD31-MVD and CD34-MVD were analyzed using Pearson’s correlation analysis. A diagnostic model of CD34-MVD was established and the diagnostic abilities of models for MVD were analyzed using receiver operating characteristic curve (ROC) analysis.

Results

There were significant differences between the quantitative parameters in the two kinds of models. Compared with SITET, DITET parameters showed better correlations with CD31-MVD and CD34-MVD. The K^trans and Ve radiomics features of the DITET model showed high efficiency for predicting the level of CD34-MVD according to ROC analysis, with areas under curves of 0.83 and 0.94, respectively.

Conclusion

Compared with SITET, the DITET model provides a better indication of the microcirculation of HCC and is thus more suitable for examining patients with HCC.

Diffusion model comparison identifies distinct tumor sub-regions and tracks treatment response

PURPOSE

MRI biomarkers of tumor response to treatment are typically obtained from parameters derived from a model applied to pre-treatment and post-treatment data. However, as tumors are spatially and temporally heterogeneous, different models may be necessary in different tumor regions, and model suitability may change over time. This work evaluates how the suitability of two diffusion-weighted (DW) MRI models varies spatially within tumors at the voxel level and in response to radiotherapy, potentially allowing inference of qualitatively different tumor microenvironments.

METHODS

DW-MRI data were acquired in CT26 subcutaneous allografts before and after radiotherapy. Restricted and time-independent diffusion models were compared, with regions well-described by the former hypothesized to reflect cellular tissue, and those well-described by the latter expected to reflect necrosis or oedema. Technical and biological validation of the percentage of tissue described by the restricted diffusion microstructural model (termed %MM) was performed through simulations and histological comparison.

RESULTS

Spatial and radiotherapy-related variation in model suitability was observed. %MM decreased from a mean of 64% at baseline to 44% 6 days post-radiotherapy in the treated group. %MM correlated negatively with the percentage of necrosis from histology, but overestimated it due to noise. Within MM regions, microstructural parameters were sensitive to radiotherapy-induced changes.

CONCLUSIONS

There is spatial and radiotherapy-related variation in different models' suitability for describing diffusion in tumor tissue, suggesting the presence of different and changing tumor sub-regions. The biological and technical validation of the proposed %MM cancer imaging biomarker suggests it correlates with, but overestimates, the percentage of necrosis.

Dynamic contrast-enhanced MRI to assess hepatocellular carcinoma response to Transarterial chemoembolization using LI-RADS criteria: A pilot study

PURPOSE

To identify quantitative dynamic contrast-enhanced (DCE)-MRI perfusion parameters indicating tumor response of hepatocellular carcinoma (HCC) to transarterial chemoembolization (TACE).

MATERIALS AND METHODS

This prospective pilot study was approved by our institutional review board; written and informed consent was obtained for each participant. Patients underwent DCE-MRI examinations before and after TACE. A variable flip-angle unenhanced 3D mDixon sequence was performed for T1 mapping. A dynamic 4D mDixon sequence was performed after contrast injection for assessing dynamic signal enhancement. Nonparametric analysis was conducted on the time-intensity curves. Parametric analysis was performed on the time-concentration curves using a dual-input single-compartment model. Treatment response according to Liver Reporting and Data System (LI-RADS) v2018 was used as the reference standard. The comparisons within groups (before vs. after treatment) and between groups (nonviable vs. equivocal or viable tumor) were performed using nonparametric bootstrap taking into account the clustering effect of lesions in patients.

RESULTS

Twenty-eight patients with 52 HCCs (size: 10-104 mm) were evaluated. For nonviable tumors (n = 27), time to peak increased from 62.5 ± 18.2 s before to 83.3 ± 12.8 s after treatment (P< 0.01). For equivocal or viable tumors (n = 25), time to peak and mean transit time significantly increased (from 54.4 ± 24.1 s to 69.5 ± 18.9 s, P < 0.01 and from 14.2 ± 11.8 s to 33.9 ± 36.8 s, P= 0.01, respectively) and the transfer constant from the extracellular and extravascular space to the central vein significantly decreased from 14.8 ± 14.1 to 8.1 ± 9.1 s^-1 after treatment (P= 0.01).

CONCLUSION

This prospective pilot DCE-MRI study showed that time to peak significantly changed after TACE treatment for both groups (nonviable tumors and equivocal or viable tumors). In our cohort, several perfusion parameters may provide an objective marker for differentiation of treatment response after TACE in HCC patients.

Detection of Sources in Non-Negative Blind Source Separation by Minimum Description Length Criterion

While non-negative blind source separation (nBSS) has found many successful applications in science and engineering, model order selection, determining the number of sources, remains a critical yet unresolved problem. Various model order selection methods have been proposed and applied to real-world data sets but with limited success, with both order over- and under-estimation reported. By studying existing schemes, we have found that the unsatisfactory results are mainly due to invalid assumptions, model oversimplification, subjective thresholding, and/or to assumptions made solely for mathematical convenience. Building on our earlier work that reformulated model order selection for nBSS with more realistic assumptions and models, we report a newly and formally revised model order selection criterion rooted in the minimum description length (MDL) principle. Adopting widely invoked assumptions for achieving a unique nBSS solution, we consider the mixing matrix as consisting of deterministic unknowns, with the source signals following a multivariate Dirichlet distribution. We derive a computationally efficient, stochastic algorithm to obtain approximate maximum-likelihood estimates of model parameters and apply Monte Carlo integration to determine the description length. Our modeling and estimation strategy exploits the characteristic geometry of the data simplex in nBSS. We validate our nBSS-MDL criterion through extensive simulation studies and on four real-world data sets, demonstrating its strong performance and general applicability to nBSS. The proposed nBSS-MDL criterion consistently detects the true number of sources, in all of our case studies.

Simultaneous acquisition sequence for improved hepatic pharmacokinetics quantification accuracy (SAHA) for dynamic contrast-enhanced MRI of liver

PURPOSE

To propose a simultaneous acquisition sequence for improved hepatic pharmacokinetics quantification accuracy (SAHA) method for liver dynamic contrast-enhanced MRI.

METHODS

The proposed SAHA simultaneously acquired high temporal-resolution 2D images for vascular input function extraction using Cartesian sampling and 3D large-coverage high spatial-resolution liver dynamic contrast-enhanced images using golden angle stack-of-stars acquisition in an interleaved way. Simulations were conducted to investigate the accuracy of SAHA in pharmacokinetic analysis. A healthy volunteer and three patients with cirrhosis or hepatocellular carcinoma were included in the study to investigate the feasibility of SAHA in vivo.

RESULTS

Simulation studies showed that SAHA can provide closer results to the true values and lower root mean square error of estimated pharmacokinetic parameters in all of the tested scenarios. The in vivo scans of subjects provided fair image quality of both 2D images for arterial input function and portal venous input function and 3D whole liver images. The in vivo fitting results showed that the perfusion parameters of healthy liver were significantly different from those of cirrhotic liver and HCC.

CONCLUSIONS

The proposed SAHA can provide improved accuracy in pharmacokinetic modeling and is feasible in human liver dynamic contrast-enhanced MRI, suggesting that SAHA is a potential tool for liver dynamic contrast-enhanced MRI. Magn Reson Med 79:2629-2641, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A semi-automatic method for the extraction of the portal venous input function in quantitative dynamic contrast-enhanced CT of the liver

OBJECTIVE

To aid the extraction of the portal venous input function (PVIF) from axial dynamic contrast-enhanced CT images of the liver, eliminating the need for full manual outlining of the vessel across time points.

METHODS

A cohort of 20 patients undergoing perfusion CT imaging of the liver was examined. Dynamic images of the liver were reformatted into contiguous thin slices. A region of interest was defined within a transverse section of the portal vein on a single contrast-enhanced image. This region of interest was then computationally projected across all thin slices for all time points to yield a semi-automated PVIF curve. This was compared against the "gold-standard" PVIF curve obtained by conventional manual outlining.

RESULTS

Bland-Altman plots of curve characteristics indicated no substantial difference between automated and manual PVIF curves [concordance correlation coefficient in the range (0.66, 0.98)]. No substantial differences were shown by Bland-Altman plots of derived pharmacokinetic parameters when a suitable kinetic model was applied in each case [concordance correlation coefficient in range (0.92, 0.95)].

CONCLUSION

This semi-automated method of extracting the PVIF performed equivalently to a "gold-standard" manual method for assessing liver function. Advances in knowledge: This technique provides a quick, simple and effective solution to the problems incurred by respiration motion and partial volume factors in the determination of the PVIF in liver dynamic contrast-enhanced CT.

Hepatic function imaging using dynamic Gd-EOB-DTPA enhanced MRI and pharmacokinetic modeling

PURPOSE

To determine whether pharmacokinetic modeling parameters with different output assumptions of dynamic contrast-enhanced MRI (DCE-MRI) using Gd-EOB-DTPA correlate with serum-based liver function tests, and compare the goodness of fit of the different output assumptions.

METHODS

A 6-min DCE-MRI protocol was performed in 38 patients. Four dual-input two-compartment models with different output assumptions and a published one-compartment model were used to calculate hepatic function parameters. The Akaike information criterion fitting error was used to evaluate the goodness of fit. Imaging-based hepatic function parameters were compared with blood chemistry using correlation with multiple comparison correction.

RESULTS

The dual-input two-compartment model assuming venous flow equals arterial flow plus portal venous flow and no bile duct output better described the liver tissue enhancement with low fitting error and high correlation with blood chemistry. The relative uptake rate Kir derived from this model was found to be significantly correlated with direct bilirubin (r = -0.52, P = 0.015), prealbumin concentration (r = 0.58, P = 0.015), and prothrombin time (r = -0.51, P = 0.026).

CONCLUSION

It is feasible to evaluate hepatic function by proper output assumptions. The relative uptake rate has the potential to serve as a biomarker of function. Magn Reson Med 78:1488-1495, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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.

The use of error-category mapping in pharmacokinetic model analysis of dynamic contrast-enhanced MRI data

This study introduces the use of 'error-category mapping' in the interpretation of pharmacokinetic (PK) model parameter results derived from dynamic contrast-enhanced (DCE-) MRI data. Eleven patients with metastatic renal cell carcinoma were enrolled in a multiparametric study of the treatment effects of bevacizumab. For the purposes of the present analysis, DCE-MRI data from two identical pre-treatment examinations were analysed by application of the extended Tofts model (eTM), using in turn a model arterial input function (AIF), an individually-measured AIF and a sample-average AIF. PK model parameter maps were calculated. Errors in the signal-to-gadolinium concentration ([Gd]) conversion process and the model-fitting process itself were assigned to category codes on a voxel-by-voxel basis, thereby forming a colour-coded 'error-category map' for each imaged slice. These maps were found to be repeatable between patient visits and showed that the eTM converged adequately in the majority of voxels in all the tumours studied. However, the maps also clearly indicated sub-regions of low Gd uptake and of non-convergence of the model in nearly all tumours. The non-physical condition ve ≥ 1 was the most frequently indicated error category and appeared sensitive to the form of AIF used. This simple method for visualisation of errors in DCE-MRI could be used as a routine quality-control technique and also has the potential to reveal otherwise hidden patterns of failure in PK model applications.

Dynamic contrast-enhanced magnetic resonance imaging of the liver: Applications in treatment of hepatic malignancies with vascular targeting agents

Dynamic contrast-enhanced magnetic resonance (DCE-MR) imaging of the liver as a trendy technique can be applied in various kinds of liver diseases to evaluate perfusion and vascular characteristics of liver tissue and tumor. It has been proved that DCE-MR imaging plays an important role in the treatment of liver malignancies with vascular targeting agents. This review aims to give an overview of DCE-MR imaging of the liver in terms of semi-quantitative analysis methods, common quantitative analysis models and contrast agents and discuss its application value in the treatment of liver malignancies with vascular targeting agents.

Direct parametric reconstruction from undersampled (k, t)-space data in dynamic contrast enhanced MRI

The Magnetic Resonance Imaging (MRI) signal can be made sensitive to functional parameters that provide information about tissues. In dynamic contrast enhanced (DCE) MRI these functional parameters are related to the microvasculature environment and the concentration changes that occur rapidly after the injection of a contrast agent. Typically DCE images are reconstructed individually and kinetic parameters are estimated by fitting a pharmacokinetic model to the time-enhancement response; these methods can be denoted as "indirect". If undersampling is present to accelerate the acquisition, techniques such as kt-FOCUSS can be employed in the reconstruction step to avoid image degradation. This paper suggests a Bayesian inference framework to estimate functional parameters directly from the measurements at high temporal resolution. The current implementation estimates pharmacokinetic parameters (related to the extended Tofts model) from undersampled (k, t)-space DCE MRI. The proposed scheme is evaluated on a simulated abdominal DCE phantom and prostate DCE data, for fully sampled, 4 and 8-fold undersampled (k, t)-space data. Direct kinetic parameters demonstrate better correspondence (up to 70% higher mutual information) to the ground truth kinetic parameters (of the simulated abdominal DCE phantom) than the ones derived from the indirect methods. For the prostate DCE data, direct kinetic parameters depict the morphology of the tumour better. To examine the impact on cancer diagnosis, a peripheral zone prostate cancer diagnostic model was employed to calculate a probability map for each method.

An investigation into the effects of temporal resolution on hepatic dynamic contrast-enhanced MRI in volunteers and in patients with hepatocellular carcinoma

This study investigated the effect of temporal resolution on the dual-input pharmacokinetic (PK) modelling of dynamic contrast-enhanced MRI (DCE-MRI) data from normal volunteer livers and from patients with hepatocellular carcinoma. Eleven volunteers and five patients were examined at 3 T. Two sections, one optimized for the vascular input functions (VIF) and one for the tissue, were imaged within a single heart-beat (HB) using a saturation-recovery fast gradient echo sequence. The data was analysed using a dual-input single-compartment PK model. The VIFs and/or uptake curves were then temporally sub-sampled (at interval ▵t = [2-20] s) before being subject to the same PK analysis. Statistical comparisons of tumour and normal tissue PK parameter values using a 5% significance level gave rise to the same study results when temporally sub-sampling the VIFs to HB < ▵t <4 s. However, sub-sampling to ▵t > 4 s did adversely affect the statistical comparisons. Temporal sub-sampling of just the liver/tumour tissue uptake curves at ▵t ≤ 20 s, whilst using high temporal resolution VIFs, did not substantially affect PK parameter statistical comparisons. In conclusion, there is no practical advantage to be gained from acquiring very high temporal resolution hepatic DCE-MRI data. Instead the high temporal resolution could be usefully traded for increased spatial resolution or SNR.