Effects of arterial input function selection on kinetic parameters in brain dynamic contrast-enhanced MRI

Arterial input function estimation compensating for inflow and partial voluming in dynamic contrast-enhanced MRI

Both inflow and the partial volume effect (PVE) are sources of error when measuring the arterial input function (AIF) in dynamic contrast-enhanced (DCE) MRI. This is relevant, as errors in the AIF can propagate into pharmacokinetic parameter estimations from the DCE data. A method was introduced for flow correction by estimating and compensating the number of the perceived pulse of spins during inflow. We hypothesized that the PVE has an impact on concentration-time curves similar to inflow. Therefore, we aimed to study the efficiency of this method to compensate for both effects simultaneously. We first simulated an AIF with different levels of inflow and PVE contamination. The peak, full width at half-maximum (FWHM), and area under curve (AUC) of the reconstructed AIFs were compared with the true (simulated) AIF. In clinical data, the PVE was included in AIFs artificially by averaging the signal in voxels surrounding a manually selected point in an artery. Subsequently, the artificial partial volume AIFs were corrected and compared with the AIF from the selected point. Additionally, corrected AIFs from the internal carotid artery (ICA), the middle cerebral artery (MCA), and the venous output function (VOF) estimated from the superior sagittal sinus (SSS) were compared. As such, we aimed to investigate the effectiveness of the correction method with different levels of inflow and PVE in clinical data. The simulation data demonstrated that the corrected AIFs had only marginal bias in peak value, FWHM, and AUC. Also, the algorithm yielded highly correlated reconstructed curves over increasingly larger neighbourhoods surrounding selected arterial points in clinical data. Furthermore, AIFs measured from the ICA and MCA produced similar peak height and FWHM, whereas a significantly larger peak and lower FWHM was found compared with the VOF. Our findings indicate that the proposed method has high potential to compensate for PVE and inflow simultaneously. The corrected AIFs could thereby provide a stable input source for DCE analysis.

Measurement variability of blood-brain barrier permeability using dynamic contrast-enhanced magnetic resonance imaging

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is used to quantify the blood-brain barrier (BBB) permeability-surface area product. Serial measurements can indicate changes in BBB health, of interest to the study of normal physiology, neurological disease, and the effect of therapeutics. We performed a scan-rescan study to inform both sample size calculation for future studies and an appropriate reference change value for patient care. The final dataset included 28 healthy individuals (mean age 53.0 years, 82% female) scanned twice with mean interval 9.9 weeks. DCE-MRI was performed at 3T using a 3D gradient echo sequence with whole brain coverage, T1 mapping using variable flip angles, and a 16-min dynamic sequence with a 3.2-s time resolution. Segmentation of white and grey matter (WM/GM) was performed using a 3D magnetization-prepared gradient echo image. The influx constant Ki was calculated using the Patlak method. The primary outcome was the within-subject coefficient of variation (CV) of Ki in both WM and GM. Ki values followed biological expectations in relation to known GM/WM differences in cerebral blood volume (CBV) and consequently vascular surface area. Subject-derived arterial input functions showed marked within-subject variability which were significantly reduced by using a venous input function (CV of area under the curve 46 vs. 12%, p < 0.001). Use of the venous input function significantly improved the CV of Ki in both WM (30 vs. 59%, p < 0.001) and GM (21 vs. 53%, p < 0.001). Further improvement was obtained using motion correction, scaling the venous input function by the artery, and using the median rather than the mean of individual voxel data. The final method gave CV of 27% and 17% in WM and GM, respectively. No further improvement was obtained by replacing the subject-derived input function by one standard population input function. CV of Ki was shown to be highly sensitive to dynamic sequence duration, with shorter measurement periods giving marked deterioration especially in WM. In conclusion, measurement variability of 3D brain DCE-MRI is sensitive to analysis method and a large precision improvement is obtained using a venous input function.

Unpaired deep learning for pharmacokinetic parameter estimation from dynamic contrast-enhanced MRI without AIF measurements

DCE-MRI provides information about vascular permeability and tissue perfusion through the acquisition of pharmacokinetic parameters. However, traditional methods for estimating these pharmacokinetic parameters involve fitting tracer kinetic models, which often suffer from computational complexity and low accuracy due to noisy arterial input function (AIF) measurements. Although some deep learning approaches have been proposed to tackle these challenges, most existing methods rely on supervised learning that requires paired input DCE-MRI and labeled pharmacokinetic parameter maps. This dependency on labeled data introduces significant time and resource constraints and potential noise in the labels, making supervised learning methods often impractical. To address these limitations, we present a novel unpaired deep learning method for estimating pharmacokinetic parameters and the AIF using a physics-driven CycleGAN approach. Our proposed CycleGAN framework is designed based on the underlying physics model, resulting in a simpler architecture with a single generator and discriminator pair. Crucially, our experimental results indicate that our method does not necessitate separate AIF measurements and produces more reliable pharmacokinetic parameters than other techniques.

Implementation and evaluation of a dynamic contrast-enhanced MR perfusion protocol for glioblastoma using a 0.35 T MRI-Linac system

PURPOSE

MRI-linear accelerator (MRI-Linac) systems allow for daily tracking of MRI changes during radiotherapy (RT). Since one common MRI-Linac operates at 0.35 T, there are efforts towards developing protocols at that field strength. In this study we demonstrate the implementation of a post-contrast 3DT1-weighted (3D-T1w) and dynamic contrast-enhancement (DCE) protocol to assess glioblastoma response to RT using a 0.35 T MRI-Linac.

METHODS AND MATERIALS

The protocol implemented was used to acquire 3D-T1w and DCE data from a flow phantom and two patients with glioblastoma (a responder and a non-responder) who underwent RT on a 0.35 T MRI-Linac. The detection of post-contrast-enhanced volumes was evaluated by comparing the 3DT1w images from the 0.35 T MRI-Linac to images obtained using a 3 T scanner. The DCE data were tested temporally and spatially using data from a flow phantom and patients. Ktrans maps were derived from DCE at three time points (a week before treatment-Pre RT, four weeks through treatment-Mid RT, and three weeks after treatment-Post RT) and were validated with patients' treatment outcomes.

RESULTS

The 3D-T1w contrast-enhancement volumes were visually and volumetrically similar between 0.35 T MRI-Linac and 3 T. DCE images showed temporal stability, and associated Ktrans maps were consistent with patient response to treatment. On average, Ktrans values showed a 54 % decrease and 8.6 % increase for a responder and non-responder respectively when Pre RT and Mid RT images were compared.

CONCLUSION

Our findings support the feasibility of obtaining post-contrast 3D-T1w and DCE data from patients with glioblastoma using a 0.35 T MRI-Linac system.

Influence of scan duration on dynamic contrast -enhanced magnetic resonance imaging pharmacokinetic parameters for brain lesions

OBJECTIVE

Gadolinium-based contrast agent needs time to leak into the extravascular-extracellular space, leak back into the vascular space, and reach an equilibrium state. For this reason, acquisition times of <10 min may cause inaccurate estimation of pharmacokinetic parameters. Since no studies have been conducted on the influence of long scan times on DCE-MRI parameters in brain tumors, the aim of this study is to investigate the variation of DCE-MRI-derived kinetic parameters as a function of acquisition time, from 5 to 10 min in brain tumors.

MATERIALS AND METHODS

Fifty-two patients with histologically confirmed brain tumors were enrolled in this retrospective study, and examination at 3 T, DCE-MRI, with scan duration of 10 min, was used for retrospective generation of 6 sets of quantitative DCE-MRI maps (Ktrans, Ve and Kep) from 5 to 10 min. Features were extracted from the DCE-MRI maps in contrast enhancement (CE) volumes. Kruskal-Wallis with post-hoc correction and coefficient of variation (CoV) were used as statistical test to compare DCE-MRI maps obtained from 6 data sets.

SIGNIFICANCE

p < 0.05.

RESULTS

No differences in Ktrans features in CE volumes between different scan durations. Ve, Kep features in CE volumes were influenced by different data length. The highest number of significantly different Ve and Kep features in CE volumes were between 5 min and 10 min (p < 0.013), 5 min and 9 min (p < 0.044), 6 min and 10 min (p < 0.040). CoV of Kep was reduced from 5 min to 10 min, going from highly variable (CoV = 0.70) to mildly variable (CoV = 0.42).

CONCLUSION

Kep and Ve were time-dependent in brain tumors, so a longer scan time is needed to obtain reliable parameter values. Ktrans was found to be time-independent, as it remains the same in all 6 acquisition times and is the only reliable parameter with short acquisition times.

A Novel Multi-Model High Spatial Resolution Method for Analysis of DCE MRI Data: Insights from Vestibular Schwannoma Responses to Antiangiogenic Therapy in Type II Neurofibromatosis

This study aimed to develop and evaluate a new DCE-MRI processing technique that combines LEGATOS, a dual-temporal resolution DCE-MRI technique, with multi-kinetic models. This technique enables high spatial resolution interrogation of flow and permeability effects, which is currently challenging to achieve. Twelve patients with neurofibromatosis type II-related vestibular schwannoma (20 tumours) undergoing bevacizumab therapy were imaged at 1.5 T both before and at 90 days following treatment. Using the new technique, whole-brain, high spatial resolution images of the contrast transfer coefficient (Ktrans), vascular fraction (vp), extravascular extracellular fraction (ve), capillary plasma flow (Fp), and the capillary permeability-surface area product (PS) could be obtained, and their predictive value was examined. Of the five microvascular parameters derived using the new method, baseline PS exhibited the strongest correlation with the baseline tumour volume (p = 0.03). Baseline ve showed the strongest correlation with the change in tumour volume, particularly the percentage tumour volume change at 90 days after treatment (p < 0.001), and PS demonstrated a larger reduction at 90 days after treatment (p = 0.0001) when compared to Ktrans or Fp alone. Both the capillary permeability-surface area product (PS) and the extravascular extracellular fraction (ve) significantly differentiated the 'responder' and 'non-responder' tumour groups at 90 days (p < 0.05 and p < 0.001, respectively). These results highlight that this novel DCE-MRI analysis approach can be used to evaluate tumour microvascular changes during treatment and the need for future larger clinical studies investigating its role in predicting antiangiogenic therapy response.

The arterial input function: Spatial dependence within the imaging volume and its influence on 3D quantitative dynamic contrast-enhanced MRI for head and neck cancer

PURPOSE

To evaluate the dependence of the arterial input function (AIF) on the imaging z-axis and its effect on 3D DCE MRI pharmacokinetic parameters as mediated by the SPGR signal equation and Extended Tofts-Kermode model.

THEORY

For SPGR-based 3D DCE MRI acquisition of the head and neck, inflow effects within vessels violate the assumptions underlying the SPGR signal model. Errors in the SPGR-based AIF estimate propagate through the Extended Tofts-Kermode model to affect the output pharmacokinetic parameters.

MATERIALS AND METHODS

3D DCE-MRI data were acquired for six newly diagnosed HNC patients in a prospective single arm cohort study. AIF were selected within the carotid arteries at each z-axis location. A region of interest (ROI) was placed in normal paravertebral muscle and the Extended Tofts-Kermode model solved for each pixel within the ROI for each AIF. Results were compared to those obtained with a published population average AIF.

RESULTS

Due to inflow effect, the AIF showed extreme variation in their temporal shapes. Ktrans was most sensitive to the initial bolus concentration and showed more variation over the muscle ROI with AIF taken from the upstream portion of the carotid. kep was less sensitive to the peak bolus concentration and showed less variation for AIF taken from the upstream portion of the carotid.

CONCLUSION

Inflow effects may introduce an unknown bias to SPGR-based 3D DCE pharmacokinetic parameters. Variation in the computed parameters depends on the selected AIF location. In the context of high flow, measurements may be limited to relative rather than absolute quantitative parameters.

Reproducibility and Optimal Arterial Input Function Selection in Dynamic Contrast-Enhanced Perfusion MRI in the Healthy Brain

BACKGROUND

Dynamic contrast-enhanced MRI (DCE-MRI) has seen increasing use for quantification of low level of blood-brain barrier (BBB) leakage in various pathological disease states and correlations with clinical outcomes. However, currently there exists limited studies on reproducibility in healthy controls, which is important for the establishment of a normality threshold for future research.

PURPOSE

To investigate the reproducibility of DCE-MRI and to evaluate the effect of arterial input function (AIF) selection and manual region of interests (ROI) delineation vs. automated global segmentation.

STUDY TYPE

Prospective.

POPULATION

A total of 16 healthy controls; 11 females; mean age 28.7 years (SD 10.1).

FIELD STRENGTH/SEQUENCE

A 3T; GE DCE; 3D TFE T1WI. 2D TSE T2.

ASSESSMENT

The influx constant K_i , a measure of BBB permeability, and V_p , the blood plasma volume, was calculated using the Patlak model. Cerebral blood flow (CBF) was calculated using Tikhonov model free deconvolution. Manual tissue ROIs, drawn by H.J.S. (30+ years of experience), were compared to automatic tissue segmentation.

STATISTICAL TESTS

Intraclass correlation coefficient (ICC) and repeatability coefficient (RC) was used to assess reproducibility. Bland-Altman plots were used to evaluate agreement between measurements day 1 vs. day 2, and manual vs. segmentation method.

RESULTS

K_i showed excellent reproducibility in both white and gray matter with an ICC between 0.79 and 0.82 and excellent agreement between manual ROI and automatic segmentation, with an ICC of 0.89 for K_i in WM. Furthermore, K_i values in gray and white matter conforms with histological tissue characteristics, where gray matter generally has a 2-fold higher vessel density. The highest reproducibility measures of K_i (ICC = 0.83), CBF (ICC = 0.77) and V_d (ICC = 0.83) was obtained with the AIF sampled in the internal carotid artery (ICA).

DATA CONCLUSION

DCE-MRI shows excellent reproducibility of pharmacokinetic variables derived from healthy controls.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Arterial input functions in dynamic susceptibility contrast MRI (DSC-MRI) in longitudinal evaluation of brain metastases

BACKGROUND

Dynamic susceptibility contrast magnetic resonance imaging (DSC-MRI) could be helpful to separate true disease progression from pseudo-progression in brain metastases when assessing the need for retreatment. However, the selection of arterial input functions (AIFs) is not standardized for analysis, limiting its use for this application.

PURPOSE

To compare population-based AIFs, AIFs specific to each patient, and AIFs specific to every visit in the longitudinal follow-up of brain metastases.

MATERIAL AND METHODS

Longitudinal data were collected from eight patients before treatment (6 of 8 patients) and after treatment (6-17 visits). Imaging was performed using a 1.5-T MRI system. Lesions were segmented by subtracting precontrast images from postcontrast images. Cerebral blood volume (rCBV) and cerebral blood flow (rCBF) were computed, and Pearson's product moment correlation coefficients were calculated to evaluate similarity of DSC parameters dependent on various AIF choices across time. AIF shape characteristics were compared. Parameter differences between white matter (WM) and gray matter (GM) were obtained to determine which AIF choice maximizes tissue differentiation.

RESULTS

Although DSC parameters follow similar patterns in time, the various AIF selections cause large parameter variations with relative standard deviations of up to ±60%. AIFs sampled in one patient across sessions more similar in shape than AIFs sampled across patients. Estimates of rCBV based on scan-specific AIFs differentiated better between perfusion in WM and GM than patient-specific or population-based AIFs (P ≤ 0.02).

CONCLUSION

Results indicate that scan-specific AIFs are the best choice for DSC-MRI parameter estimations in the longitudinal follow-up of brain metastases.

Fluid Mechanics Approach to Perfusion Quantification: Vasculature Computational Fluid Dynamics Simulation, Quantitative Transport Mapping (QTM) Analysis of Dynamics Contrast Enhanced MRI, and Application in Nonalcoholic Fatty Liver Disease Classification

OBJECTIVE

We quantify liver perfusion using quantitative transport mapping (QTM) method that is free of arterial input function (AIF). QTM method is validated in a vasculature computational fluid dynamics (CFD) simulation and is applied for processing dynamic contrast enhanced (DCE) MRI images in differentiating liver with nonalcoholic fatty liver disease (NAFLD) from healthy controls using pathology reference in a preclinical rabbit model.

METHODS

QTM method was validated on a liver perfusion simulation based on fluid dynamics using a rat liver vasculature model and the mass transport equation. In the NAFLD grading task, DCE MRI images of 7 adult rabbits with methionine choline-deficient diet-induced nonalcoholic steatohepatitis (NASH), 8 adult rabbits with simple steatosis (SS) were acquired and processed using QTM method and dual-input two compartment Kety's method respectively. Statistical analysis was performed on six perfusion parameters: velocity magnitude | u | derived from QTM, liver arterial blood flow LBFa, liver venous blood flow LBFv, permeability Ktrans, blood volume Vp and extravascular space volume Ve averaged in liver ROI.

RESULTS

In the simulation, QTM method successfully reconstructed blood flow, reduced error by 48% compared to Kety's method. In the preclinical study, only QTM |u| showed significant difference between high grade NAFLD group and low grade NAFLD group.

CONCLUSION

QTM postprocesses DCE-MRI automatically through deconvolution in space and time to solve the inverse problem of the transport equation. Comparing with Kety's method, QTM method showed higher accuracy and better differentiation in NAFLD classification task.

SIGNIFICANCE

We propose to apply QTM method in liver DCE MRI perfusion quantification.

Evaluation of blood-brain barrier integrity by the analysis of dynamic contrast-enhanced MRI - a comparison of quantitative and semi-quantitative methods

Disruption of the blood-brain barrier (BBB) is a key feature of various brain disorders. To assess its integrity a parametrization of dynamic magnetic resonance imaging (DCE MRI) with a contrast agent (CA) is broadly used. Parametrization can be done quantitatively or semi-quantitatively. Quantitative methods directly describe BBB permeability but exhibit several drawbacks such as high computation demands, reproducibility issues, or low robustness. Semi-quantitative methods are fast to compute, simply mathematically described, and robust, however, they do not describe the status of BBB directly but only as a variation of CA concentration in measured tissue. Our goal was to elucidate differences between five semi-quantitative parameters: maximal intensity (Imax), normalized permeability index (NPI), and difference in DCE values between three timepoints: baseline, 5 min, and 15 min (delta5-0, delta15-0, delta15-5) and two quantitative parameters: transfer constant (Ktrans) and an extravascular fraction (Ve). For the purpose of comparison, we analyzed DCE data of four patients 12-15 days after the stroke with visible CA enhancement. Calculated parameters showed abnormalities spatially corresponding with the ischemic lesion, however, findings in individual parameters morphometrically differed. Ktrans and Ve were highly correlated. Delta5-0 and delta15-0 were prominent in regions with rapid CA enhancement and highly correlated with Ktrans. Abnormalities in delta15-5 and NPI were more homogenous with less variable values, smoother borders, and less detail than Ktrans. Moreover, only delta15-5 and NPI were able to distinguish vessels from extravascular space. Our comparison provides important knowledge for understanding and interpreting parameters derived from DCE MRI by both quantitative and semi-quantitative methods.

Influence of Different Measurement Methods of Arterial Input Function on Quantitative Dynamic Contrast-Enhanced MRI Parameters in Head and Neck Cancer

BACKGROUND

Head and neck cancer (HNC) is the sixth most prevalent cancer worldwide. Dynamic contrast-enhanced MRI (DCE-MRI) helps in diagnosis and prognosis. Quantitative DCE-MRI requires an arterial input function (AIF), which affects the values of pharmacokinetic parameters (PKP).

PURPOSE

To evaluate influence of four individual AIF measurement methods on quantitative DCE-MRI parameters values (K^trans , v_e , k_ep , and v_p ), for HNC and muscle.

STUDY TYPE

Prospective.

POPULATION

A total of 34 HNC patients (23 males, 11 females, age range 24-91) FIELD STRENGTH/SEQUENCE: A 3 T; 3D SPGR gradient echo sequence with partial saturation of inflowing spins.

ASSESSMENT

Four AIF methods were applied: automatic AIF (AIFa) with up to 50 voxels selected from the whole FOV, manual AIF (AIFm) with four voxels selected from the internal carotid artery, both conditions without (Mc-) or with (Mc+) motion correction. Comparison endpoints were peak AIF values, PKP values in tumor and muscle, and tumor/muscle PKP ratios.

STATISTICAL TESTS

Nonparametric Friedman test for multiple comparisons. Nonparametric Wilcoxon test, without and with Benjamini Hochberg correction, for pairwise comparison of AIF peak values and PKP values for tumor, muscle and tumor/muscle ratio, P value ≤ 0.05 was considered statistically significant.

RESULTS

Peak AIF values differed significantly for all AIF methods, with mean AIFmMc+ peaks being up to 66.4% higher than those for AIFaMc+. Almost all PKP values were significantly higher for AIFa in both, tumor and muscle, up to 76% for mean K^trans values. Motion correction effect was smaller. Considering tumor/muscle parameter ratios, most differences were not significant (0.068 ≤ Wilcoxon P value ≤ 0.8).

DATA CONCLUSION

We observed important differences in PKP values when using either AIFa or AIFm, consequently choice of a standardized AIF method is mandatory for DCE-MRI on HNC. From the study findings, AIFm and inflow compensation are recommended. The use of the tumor/muscle PKP ratio should be of interest for multicenter studies.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Maximum Entropy Technique and Regularization Functional for Determining the Pharmacokinetic Parameters in DCE-MRI

This paper aims to solve the arterial input function (AIF) determination in dynamic contrast-enhanced MRI (DCE-MRI), an important linear ill-posed inverse problem, using the maximum entropy technique (MET) and regularization functionals. In addition, estimating the pharmacokinetic parameters from a DCE-MR image investigations is an urgent need to obtain the precise information about the AIF-the concentration of the contrast agent on the left ventricular blood pool measured over time. For this reason, the main idea is to show how to find a unique solution of linear system of equations generally in the form of [Formula: see text] named an ill-conditioned linear system of equations after discretization of the integral equations, which appear in different tomographic image restoration and reconstruction issues. Here, a new algorithm is described to estimate an appropriate probability distribution function for AIF according to the MET and regularization functionals for the contrast agent concentration when applying Bayesian estimation approach to estimate two different pharmacokinetic parameters. Moreover, by using the proposed approach when analyzing simulated and real datasets of the breast tumors according to pharmacokinetic factors, it indicates that using Bayesian inference-that infer the uncertainties of the computed solutions, and specific knowledge of the noise and errors-combined with the regularization functional of the maximum entropy problem, improved the convergence behavior and led to more consistent morphological and functional statistics and results. Finally, in comparison to the proposed exponential distribution based on MET and Newton's method, or Weibull distribution via the MET and teaching-learning-based optimization (MET/TLBO) in the previous studies, the family of Gamma and Erlang distributions estimated by the new algorithm are more appropriate and robust AIFs.

Surrogate vascular input function measurements from the superior sagittal sinus are repeatable and provide tissue-validated kinetic parameters in brain DCE-MRI

Accurate vascular input function (VIF) derivation is essential in brain dynamic contrast-enhanced (DCE) MRI. The optimum site for VIF estimation is, however, debated. This study sought to compare VIFs extracted from the internal carotid artery (ICA) and its branches with an arrival-corrected vascular output function (VOF) derived from the superior sagittal sinus (VOF_SSS). DCE-MRI datasets from sixty-six patients with different brain tumours were retrospectively analysed and plasma gadolinium-based contrast agent (GBCA) concentration-time curves used to extract VOF/VIFs from the SSS, the ICA, and the middle cerebral artery. Semi-quantitative parameters across each first-pass VOF/VIF were compared and the relationship between these parameters and GBCA dose was evaluated. Through a test-retest study in 12 patients, the repeatability of each semiquantitative VOF/VIF parameter was evaluated; and through comparison with histopathological data the accuracy of kinetic parameter estimates derived using each VOF/VIF and the extended Tofts model was also assessed. VOF_SSS provided a superior surrogate global input function compared to arteries, with greater contrast-to-noise (p < 0.001), higher peak (p < 0.001, repeated-measures ANOVA), and a greater sensitivity to interindividual plasma GBCA concentration. The repeatability of VOF_SSS derived semi-quantitative parameters was good to excellent (ICC = 0.717-0.888) outperforming arterial based approaches. In contrast to arterial VIFs, kinetic parameters obtained using a SSS derived VOF permitted detection of intertumoural differences in both microvessel surface area and cell density within resected tissue specimens. These results support the usage of an arrival-corrected VOF_SSS as a surrogate vascular input function for kinetic parameter mapping in brain DCE-MRI.

Quantitative transport mapping (QTM) for differentiating benign and malignant breast lesion: Comparison with traditional kinetics modeling and semi-quantitative enhancement curve characteristics

PURPOSE

To test the feasibility of using quantitative transport mapping (QTM) method, which is based on the inversion of transport equation using spatial deconvolution without any arterial input function, for automatically postprocessing dynamic contrast enhanced MRI (DCE-MRI) to differentiate malignant and benign breast tumors.

MATERIALS AND METHODS

Breast DCE-MRI data with biopsy confirmed malignant (n = 13) and benign tumors (n = 13) was used to assess QTM velocity (|u|) and diffusion coefficient (D), volume transfer constant (Ktrans), volume fraction of extravascular extracellular space (Ve) from kinetics method, and traditional enhancement curve characteristics (ECC: amplitude A, wash-in rate α, wash-out rate β). A Mann-Whitney U test and receiver operating characteristic curve (ROC) analysis were performed to assess the diagnostic performance of these parameters for distinguishing between benign and malignant tumors.

RESULTS

Between malignant and benign tumors, there was a significant difference in |u| and Ktrans, (p = 0.0066, 0.0274, respectively), but not in D, Ve, A, α and β (p = 0.1119, 0.2382, 0.4418,0.2592 and 0.9591, respectively). ROC area-under-the-curve was 0.82, 0.75 (95% confidence level 0.60-0.95, 0.51-0.90) for |u| and Ktrans, respectively.

CONCLUSION

QTM postprocesses DCE-MRI automatically through deconvolution in space and time to solve the inverse problem of the transport equation. Comparing with traditional kinetics method and ECC, QTM method showed better diagnostic accuracy in differentiating benign from malignant breast tumors in this study.

Modified Maximum Entropy Method and Estimating the AIF via DCE-MRI Data Analysis

Background: For the kinetic models used in contrast-based medical imaging, the assignment of the arterial input function named AIF is essential for the estimation of the physiological parameters of the tissue via solving an optimization problem. Objective: In the current study, we estimate the AIF relayed on the modified maximum entropy method. The effectiveness of several numerical methods to determine kinetic parameters and the AIF is evaluated—in situations where enough information about the AIF is not available. The purpose of this study is to identify an appropriate method for estimating this function. Materials and Methods: The modified algorithm is a mixture of the maximum entropy approach with an optimization method, named the teaching-learning method. In here, we applied this algorithm in a Bayesian framework to estimate the kinetic parameters when specifying the unique form of the AIF by the maximum entropy method. We assessed the proficiency of the proposed method for assigning the kinetic parameters in the dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), when determining AIF with some other parameter-estimation methods and a standard fixed AIF method. A previously analyzed dataset consisting of contrast agent concentrations in tissue and plasma was used. Results and Conclusions: We compared the accuracy of the results for the estimated parameters obtained from the MMEM with those of the empirical method, maximum likelihood method, moment matching (“method of moments”), the least-square method, the modified maximum likelihood approach, and our previous work. Since the current algorithm does not have the problem of starting point in the parameter estimation phase, it could find the best and nearest model to the empirical model of data, and therefore, the results indicated the Weibull distribution as an appropriate and robust AIF and also illustrated the power and effectiveness of the proposed method to estimate the kinetic parameters.

The blood-brain barrier disruption after syncope: a dynamic contrast-enhanced magnetic resonance imaging study: A case report

RATIONALE

Using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), we demonstrated blood-brain barrier (BBB) disruption following syncope.

PATIENT CONCERNS

A 45-year-old man experienced syncope with a chief complaint of syncope (duration: 1 minutes), 1 day before visiting a university hospital for examination. He had no history of medical problems and was not taking any medications. This episode was the first in his lifetime.

DIAGNOSES

After syncope, the patient did not have any illnesses or symptoms, such as headache, cognitive deficits, or somnolence.

INTERVENTIONS

Cardiac evaluation did not reveal any abnormal findings. In addition, in conventional brain and chest computed tomography and brain MRI, no abnormal lesions were observed.

OUTCOMES

DCE-MRI of the patient showed bright blue colored lines within the sulci throughout the cerebral cortex. The regions of interest, including bright blue colored lines, had significantly higher Ktrans values (6.86 times higher) than those in healthy control participants. These findings are indicative of BBB disruption of the vessels in the sulci.

LESSONS

Using DCE-MRI, we demonstrated BBB disruption following syncope. DCE-MRI is a useful tool for the detection of BBB disruption following syncope.

Blood-Brain Barrier Permeability in Patients With Reversible Cerebral Vasoconstriction Syndrome Assessed With Dynamic Contrast-Enhanced MRI

BACKGROUND AND OBJECTIVES

Blood-brain barrier (BBB) disruption has been proposed to be important in the pathogenesis of reversible cerebral vasoconstriction syndrome (RCVS), but not all patients present an identifiable macroscopic BBB disruption; that is, visible contrast leakage on contrast-enhanced T2 fluid-attenuated inversion recovery imaging. This study aimed to evaluate microscopic BBB permeability and its dynamic change in patients with RCVS.

METHODS

This prospective cohort implemented 3T dynamic contrast-enhanced MRI. We measured microscopic BBB permeability by determining the whole-brain and white matter hyperintensity (WMH) K^trans values and evaluated the correlation of whole-brain K^trans permeability with clinical and vascular measures in transcranial color-coded sonography.

RESULTS

In total, 176 patients (363 scans) were analyzed and separated into acute (≦30 days) and remission (≧90 days) groups based on the onset-to-examination time. Whole-brain K^trans values were similar between patients with and without macroscopic BBB disruption in either acute or remission stage. The whole-brain K^trans was significantly decreased (p < 0.001) from acute to remission stages. The WMH K^trans was significantly higher than mirror references and decreased from acute to remission stages (p < 0.001). Whole-brain K^trans correlated with mean pulsatility index (r _s = 0.5, p = 0.029), mean resistance index (r _s = 0.662, p = 0.002), and distal-to-proximal ratio of resistance index (r _s = 0.801, p < 0.001) of M1 segment of middle cerebral arteries at around 10-15 days after onset. The time-trend curve of whole-brain K^trans depicted dynamic changes during disease course, similar to temporal trends of vasoconstrictions and WMH.

DISCUSSION

Patients with RCVS presented increased microscopic brain permeability during acute stage, even without discernible macroscopic BBB disruption. The dynamic changes in BBB permeability may be related to impaired cerebral microvascular compliance and WMH formation.

How tissue T1-variability influences DCE-MRI perfusion parameters estimation of recurrent high-grade glioma after surgery followed by radiochemotherapy

BACKGROUND

Quantification of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) kinetic parameters (KPs) requires a determination of native tissue T1. Two approaches are adopted: (i) tissue T1-maps are acquired; and (ii) an a priori T1 value (fT1) is fixed for all patients (fT1-approach). Although it is more attractive, the fT1-approach might bias the results of KP calculations due to tissue T1 variability.

PURPOSE

To quantify the tissue T1 variability of recurrent high-grade glioma (HGG) and the error in KP estimation when the fT1-approach is adopted.

MATERIAL AND METHODS

We reviewed the postoperative MRI scans of 28 patients with recurrent HGG after radiochemotherapy. MRI study included T1-maps from multiple-dynamic multiple-echo imaging, DCE-MRI, and contrast enhanced T1-weighted images. KPs were calculated using T1-map and fT1-approach.

RESULTS

The tissue T1 variability of recurrent HGG was relevant. The absolute error in KP estimation, as a function of the deviation of fT1 from the true value, was 8% every 100 ms. The difference between the KPs obtained with fT1-approach from fT1 values of 1300, 1390, and 1500 ms and their reference values were mostly within the 95% confidence interval (± 1.96 standard deviation). Conversely, using fT1 values of 900, 1200, 1600, and 1900 ms causes a significant error in KP estimation (P<0.05).

CONCLUSION

Recurrent HGG is characterized by a substantial T1 variability. Although the fT1-approach does not account for this variability, it results in a minor effect on the KP estimations provided the fT1 value is in the range of 1300-1500 ms.

A comprehensive evaluation and impact of normalization of generalized tracer kinetic model parameters to characterize blood-brain-barrier permeability in normal-appearing and tumor tissue regions of patients with glioma

RATIONALE AND OBJECTIVES

To comprehensively evaluate robustness and variations of DCE-MRI derived generalized-tracer-kinetic-model (GTKM) parameters in healthy and tumor tissues and impact of normalization in mitigating these variations on application to glioma.

MATERIALS (PATIENTS) AND METHODS

A retrospective study included pre-operative 31 high-grade-glioma(HGG), 22 low-grade-glioma(LGG) and 33 follow-up data from 10 patients a prospective study with 4 HGG subjects. Voxel-wise GTKM was fitted to DCE-MRI data to estimate K^trans, v_e, v_b. Simulations were used to evaluate noise sensitivity. Variation of parameters with-respect-to arterial-input-function (AIF) variation and data length were studied. Normalization of parameters with-respect-to mean values in gray-matter (GM) and white-matter (WM) regions (GM-Type-2, WM-Type-2) and mean curves (GM-Type-1, WM-Type-1) were also evaluated. Co-efficient-of-variation(CoV), relative-percentage-error (RPE), Box-Whisker plots, bar graphs and t-test were used for comparison.

RESULTS

GTKM was fitted well in all tissue regions. K^trans and v_e in contrast-enhancing (CE) has shown improved noise sensitivity in longer data. v_b was reliable in all tissues. Mean AIF and C(t) peaks showed ~38% and ~35% variations. During simulation, normalizations have mitigated variations due to changes in AIF amplitude in K^trans and v_b.. v_e was less sensitive to normalizations. CoV of K^trans and v_b has reduced ~70% after GM-Type-1 normalization and ~80% after GM-Type-2 normalization, respectively. GM-Type-1 (p = 0.003) and GM-Type-2 (p = 0.006) normalizations have significantly improved differentiation of HGG and LGG using K^trans.

CONCLUSION

K^trans and v_b can be reliably estimated in normal-appearing brain tissues and can be used for normalization of corresponding parameters in tumor tissues for mitigating inter-subject variability due to errors in AIF. Normalized K^trans and v_b provided improved differentiation of HGG and LGG.

Robust arterial input function surrogate measurement from the superior sagittal sinus complex signal for fast dynamic contrast-enhanced MRI in the brain

PURPOSE

Accurately estimating the arterial input function for dynamic contrast-enhanced MRI is challenging. An arterial input function is typically determined from signal magnitude changes related to a contrast agent, often leading to underestimation of peak concentrations. Alternatively, signal phase recovers the accurate peak concentration for straight vessels but suffers from high noise. A recent method proposed to fit the signal in the complex plane by combining the advantages of the previous 2 methods. The purpose of this work is to refine this complex-based method to determine the venous output function (VOF), an arterial input function surrogate, from the superior sagittal sinus.

METHODS

We propose a state-of-the-art complex-based method that includes direct compensation for blood inflow and signal phase correction accounting for the curvature of the superior sagittal sinus, generally assumed collinear with B₀ . We compared the magnitude-, phase-, and complex-based VOF determination methods against various simulated biases as well as for 29 brain metastases patients.

RESULTS

Angulation of the superior sagittal sinus relative to B₀ varied widely within patients, and its effect on the signal phase caused an underestimation of peak concentrations of up to 65%. Correction significantly increased the VOF peak concentration for the phase- and complex-based VOFs in the cohort. The phase-based method recovered accurate peak concentrations but lacked precision in the tail of the VOF. Our complex-based VOF completely recovered the effect of inflow and resulted in a high-peak concentration with limited noise.

CONCLUSION

The new complex-based method resulted in high-quality VOF robust against superior sagittal sinus curvature and variations in patient positioning.

Operator dependency of arterial input function in dynamic contrast-enhanced MRI

Objective

To investigate the effect of inter-operator variability in arterial input function (AIF) definition on kinetic parameter estimates (KPEs) from dynamic contrast-enhanced (DCE) MRI in patients with high-grade gliomas.

Methods

The study included 118 DCE series from 23 patients. AIFs were measured by three domain experts (DEs), and a population AIF (pop-AIF) was constructed from the measured AIFs. The DE-AIFs, pop-AIF and AUC-normalized DE-AIFs were used for pharmacokinetic analysis with the extended Tofts model. AIF-dependence of KPEs was assessed by intraclass correlation coefficient (ICC) analysis, and the impact on relative longitudinal change in K^trans was assessed by Fleiss’ kappa (κ).

Results

There was a moderate to substantial agreement (ICC 0.51–0.76) between KPEs when using DE-AIFs, while AUC-normalized AIFs yielded ICC 0.77–0.95 for K^trans, k_ep and v_e and ICC 0.70 for v_p. Inclusion of the pop-AIF did not reduce agreement. Agreement in relative longitudinal change in K^trans was moderate (κ = 0.591) using DE-AIFs, while AUC-normalized AIFs gave substantial (κ = 0.809) agreement.

Discussion

AUC-normalized AIFs can reduce the variation in kinetic parameter results originating from operator input. The pop-AIF presented in this work may be applied in absence of a satisfactory measurement.

7T dynamic contrast-enhanced MRI for the detection of subtle blood-brain barrier leakage

Background and Purpose

Dynamic contrast‐enhanced MRI (DCE‐MRI) can be employed to assess the blood–brain barrier (BBB) integrity. Detection of BBB leakage at lower field strengths (≤3T) is cumbersome as the signal is noisy, while leakage can be subtle. Utilizing the increased signal‐to‐noise ratio at higher field strengths, we explored the application of 7T DCE‐MRI for assessing BBB leakage.

Methods

A dual‐time resolution DCE‐MRI method was implemented at 7T and a slow injection rate (0.3 ml/s) and low dose (3 mmol) served to obtain signal changes linearly related to the gadolinium concentration, that is, minimized for T₂^* degradation effects. With the Patlak graphical approach, the leakage rate (K_i) and blood plasma volume fraction (v_p) were calculated. The method was evaluated in 10 controls, an ischemic stroke patient, and a patient with a transient ischemic attack.

Results

K_i and v_p were significantly higher in gray matter compared to white matter of all participants. These K_i values were higher in both patients compared to the control subjects. Finally, for the lesion identified in the ischemic stroke patient, higher leakage values were observed compared to normal‐appearing tissue.

Conclusion

We demonstrate how a dual‐time resolution DCE‐MRI protocol at 7T, with administration of half the clinically used contrast agent dose, can be used for assessing subtle BBB leakage. Although the feasibility of DCE‐MRI for assessing the BBB integrity at 3T is well known, we showed that a continuous sampling DCE‐MRI method tailored for 7T is also capable of assessing leakage with a high sensitivity over a range of K_i values.

Sources of systematic error in DCE-MRI estimation of low-level blood-brain barrier leakage

PURPOSE

Dynamic contrast-enhanced (DCE) -MRI with Patlak model analysis is increasingly used to quantify low-level blood-brain barrier (BBB) leakage in studies of pathophysiology. We aimed to investigate systematic errors due to physiological, experimental, and modeling factors influencing quantification of the permeability-surface area product PS and blood plasma volume v_p , and to propose modifications to reduce the errors so that subtle differences in BBB permeability can be accurately measured.

METHODS

Simulations were performed to predict the effects of potential sources of systematic error on conventional PS and v_p quantification: restricted BBB water exchange, reduced cerebral blood flow, arterial input function (AIF) delay and B 1 + error. The impact of targeted modifications to the acquisition and processing were evaluated, including: assumption of fast versus no BBB water exchange, bolus versus slow injection of contrast agent, exclusion of early data from model fitting and B 1 + correction. The optimal protocol was applied in a cohort of recent mild ischaemic stroke patients.

RESULTS

Simulation results demonstrated substantial systematic errors due to the factors investigated (absolute PS error ≤ 4.48 × 10^-4 min^-1 ). However, these were reduced (≤0.56 × 10^-4 min^-1 ) by applying modifications to the acquisition and processing pipeline. Processing modifications also had substantial effects on in-vivo normal-appearing white matter PS estimation (absolute change ≤ 0.45 × 10^-4 min^-1 ).

CONCLUSION

Measuring subtle BBB leakage with DCE-MRI presents unique challenges and is affected by several confounds that should be considered when acquiring or interpreting such data. The evaluated modifications should improve accuracy in studies of neurodegenerative diseases involving subtle BBB breakdown.

Quantitative transport mapping (QTM) of the kidney with an approximate microvascular network

PURPOSE

Proof-of-concept study of mapping renal blood flow vector field according to the inverse solution to a mass transport model of time resolved tracer-labeled MRI data.

THEORY AND METHODS

To determine tissue perfusion according to the underlying physics of spatiotemporal tracer concentration variation, the mass transport equation is integrated over a voxel with an approximate microvascular network for fitting time-resolved tracer imaging data. The inverse solution to the voxelized transport equation provides the blood flow vector field, which is referred to as quantitative transport mapping (QTM). A numerical microvascular network modeling the kidney with computational fluid dynamics reference was used to verify the accuracy of QTM and the current Kety's method that uses a global arterial input function. Multiple post-label delay arterial spin labeling (ASL) of the kidney on seven subjects was used to assess QTM in vivo feasibility.

RESULTS

Against the ground truth in the numerical model, the error in flow estimated by QTM (18.6%) was smaller than that in Kety's method (45.7%, 2.5-fold reduction). The in vivo kidney perfusion quantification by QTM (cortex: 443 ± 58 mL/100 g/min and medulla: 190 ± 90 mL/100 g/min) was in the range of that by Kety's method (482 ± 51 mL/100 g/min in the cortex and 242 ± 73 mL/100 g/min in the medulla), and QTM provided better flow homogeneity in the cortex region.

CONCLUSIONS

QTM flow velocity mapping is feasible from multi-delay ASL MRI data based on inverting the transport equation. In a numerical simulation, QTM with deconvolution in space and time provided more accurate perfusion quantification than Kety's method with deconvolution in time only.

Effects of artery input function on dynamic contrast-enhanced MRI for determining grades of gliomas

OBJECTIVE

To evaluate the effect of artery input function (AIF) derived from different arteries for pharmacokinetic modeling on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) parameters in the grading of gliomas.

METHODS

49 patients with pathologically confirmed gliomas were recruited and underwent DCE-MRI. A modified Tofts model with different AIFs derived from anterior cerebral artery (ACA), ipsilateral and contralateral middle cerebral artery (MCA) and posterior cerebral artery (PCA) was used to estimate quantitative parameters such as K^trans (volume transfer constant) and V_e (fractional extracellular-extravascular space volume) for distinguishing the low grade glioma from high grade glioma. The K^trans and V_e were compared between different arteries using Two Related Samples Tests (TRST) (i.e. Wilcoxon Signed Ranks Test). In addition, these parameters were compared between the low and high grades as well as between the grade II and III using the Mann-Whitney U-test. A p-value of less than 0.05 was regarded as statistically significant.

RESULTS

All the patients completed the DCE-MRI successfully. Sharp wash-in and wash-out phases were observed in all AIFs derived from the different arteries. The quantitative parameters (K^trans and V_e) calculated from PCA were significant higher than those from ACA and MCA for low and high grades, respectively (p < 0.05). Despite the differences of quantitative parameters derived from ACA, MCA and PCA, the K^trans and V_e from any AIFs could distinguish between low and high grade, however, only K^trans from any AIFs could distinguish grades II and III. There was no significant correlation between parameters and the distance from the artery, which the AIF was extracted, to the tumor.

CONCLUSION

Both quantitative parameters K^trans and V_e calculated using any AIF of ACA, MCA, and PCA can be used for distinguishing the low- from high-grade gliomas, however, only K^trans can distinguish grades II and III.

ADVANCES IN KNOWLEDGE

We sought to assess the effect of AIF on DCE-MRI for determining grades of gliomas. Both quantitative parameters K^trans and V_e calculated using any AIF of ACA, MCA, and PCA can be used for distinguishing the low- from high-grade gliomas.

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.

Multiparametric MRI and Coregistered Histology Identify Tumor Habitats in Breast Cancer Mouse Models

It is well-recognized that solid tumors are genomically, anatomically, and physiologically heterogeneous. In general, more heterogeneous tumors have poorer outcomes, likely due to the increased probability of harboring therapy-resistant cells and regions. It is hypothesized that the genomic and physiologic heterogeneity are related, because physiologically distinct regions will exert variable selection pressures leading to the outgrowth of clones with variable genomic/proteomic profiles. To investigate this, methods must be in place to interrogate and define, at the microscopic scale, the cytotypes that exist within physiologically distinct subregions ("habitats") that are present at mesoscopic scales. MRI provides a noninvasive approach to interrogate physiologically distinct local environments, due to the biophysical principles that govern MRI signal generation. Here, we interrogate different physiologic parameters, such as perfusion, cell density, and edema, using multiparametric MRI (mpMRI). Signals from six different acquisition schema were combined voxel-by-voxel into four clusters identified using a Gaussian mixture model. These were compared with histologic and IHC characterizations of sections that were coregistered using MRI-guided 3D printed tumor molds. Specifically, we identified a specific set of MRI parameters to classify viable-normoxic, viable-hypoxic, nonviable-hypoxic, and nonviable-normoxic tissue types within orthotopic 4T1 and MDA-MB-231 breast tumors. This is the first coregistered study to show that mpMRI can be used to define physiologically distinct tumor habitats within breast tumor models. SIGNIFICANCE: This study demonstrates that noninvasive imaging metrics can be used to distinguish subregions within heterogeneous tumors with histopathologic correlation.

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 K^trans (volume transfer rate constant), v_e (extravascular, extracellular volume fraction), k_ep (efflux rate constant), and τ_i (mean intracellular water lifetime). All other variables, including the definition of the tumor region of interest and precontrast T₁ 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 K^trans, v_e, k_ep, and τ_i, respectively, using the unadjusted AIFs. Use of the reference tissue-adjusted AIFs reduced variations in K^trans and v_e (wCV = 0.50 and 0.10, respectively), but had smaller effects on k_ep and τ_i (wCV = 0.39 and 0.22, respectively). k_ep is less sensitive to AIF variation than K^trans, 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.

Perfusion MRI in treatment evaluation of glioblastomas: Clinical relevance of current and future techniques

Treatment evaluation of patients with glioblastomas is important to aid in clinical decisions. Conventional MRI with contrast is currently the standard method, but unable to differentiate tumor progression from treatment-related effects. Pseudoprogression appears as new enhancement, and thus mimics tumor progression on conventional MRI. Contrarily, a decrease in enhancement or edema on conventional MRI during antiangiogenic treatment can be due to pseudoresponse and is not necessarily reflective of a favorable outcome. Neovascularization is a hallmark of tumor progression but not for posttherapeutic effects. Perfusion-weighted MRI provides a plethora of additional parameters that can help to identify this neovascularization. This review shows that perfusion MRI aids to identify tumor progression, pseudoprogression, and pseudoresponse. The review provides an overview of the most applicable perfusion MRI methods and their limitations. Finally, future developments and remaining challenges of perfusion MRI in treatment evaluation in neuro-oncology are discussed. Level of Evidence: 3 Technical Efficacy: Stage 4 J. Magn. Reson. Imaging 2019;49:11-22.

Meningioma assessment: Kinetic parameters in dynamic contrast-enhanced MRI appear independent from microvascular anatomy and VEGF expression

BACKGROUND AND PURPOSE

Kinetic parameters of T1-weighted dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) are considered to be influenced by microvessel environment. This study was performed to explore the extent of this association for meningiomas.

MATERIALS AND METHODS

DCE-MRI kinetic parameters (contrast agent transfer constants K^trans and k_ep, volume fractions v_p and v_e) were determined in pre-operative 3T MRI of meningioma patients for later biopsy sites (19 patients; 15 WHO I^o, no previous radiation, and 4 WHO III^o pre-radiated recurrent tumors). Sixty-three navigated biopsies were consecutively retrieved. Biopsies were immunohistochemically investigated with endothelial marker CD34 and VEGF antibodies, stratified in a total of 4383 analysis units and computationally assessed for VEGF expression and vascular parameters (vessel density, vessel quantity, vascular fraction within tissue [vascular area ratio], vessel wall thickness). Derivability of kinetic parameters from VEGF expression or microvascularization was determined by mixed linear regression analysis. Tissue kinetic and microvascular parameters were tested for their capacity to identify the radiation status in a subanalysis.

RESULTS

Kinetic parameters were neither significantly related to the corresponding microvascular parameters nor to tissue VEGF expression. There was no significant association between microvessel density and its presumed correlate v_p (P=0.07). The subgroup analysis of high-grade radiated meningiomas showed a significantly reduced microvascular density (AUC 0.91; P<0.0001) and smaller total vascular fraction (AUC 0.73; P=0.01).

CONCLUSIONS

In meningioma, DCE-MRI kinetic parameters neither allow for a reliable prediction of tumor microvascularization, nor for a prediction of VEGF expression. Kinetic parameters seem to be determined from different independent factors.

Robust and efficient pharmacokinetic parameter non-linear least squares estimation for dynamic contrast enhanced MRI of the prostate

PURPOSE

To describe an efficient numerical optimization technique using non-linear least squares to estimate perfusion parameters for the Tofts and extended Tofts models from dynamic contrast enhanced (DCE) MRI data and apply the technique to prostate cancer.

METHODS

Parameters were estimated by fitting the two Tofts-based perfusion models to the acquired data via non-linear least squares. We apply Variable Projection (VP) to convert the fitting problem from a multi-dimensional to a one-dimensional line search to improve computational efficiency and robustness. Using simulation and DCE-MRI studies in twenty patients with suspected prostate cancer, the VP-based solver was compared against the traditional Levenberg-Marquardt (LM) strategy for accuracy, noise amplification, robustness to converge, and computation time.

RESULTS

The simulation demonstrated that VP and LM were both accurate in that the medians closely matched assumed values across typical signal to noise ratio (SNR) levels for both Tofts models. VP and LM showed similar noise sensitivity. Studies using the patient data showed that the VP method reliably converged and matched results from LM with approximate 3× and 2× reductions in computation time for the standard (two-parameter) and extended (three-parameter) Tofts models. While LM failed to converge in 14% of the patient data, VP converged in the ideal 100%.

CONCLUSION

The VP-based method for non-linear least squares estimation of perfusion parameters for prostate MRI is equivalent in accuracy and robustness to noise, while being more reliably (100%) convergent and computationally about 3× (TM) and 2× (ETM) faster than the LM-based method.