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

Combination of IVIM with DCE-MRI for diagnostic and prognostic evaluation of breast cancer

PURPOSE

To identify the most effective combination of DCE-MRI (Ktrans,Kep) and IVIM (D,f) and analyze the correlations of these parameters with prognostic indicators (ER, PR, and HER2, Ki-67 index, axillary lymph node (ALN) and tumor size) to improve the diagnostic and prognostic efficiency in breast cancer.

METHODS

This is a prospective study. We performed T1WI, T2WI, IVIM, DCE-MRI at 3 T MRI examinations on benign and malignant breast lesions that met the inclusion criteria. We also collected pathological results of corresponding lesions, including ER, PR, and HER2, Ki-67 index, axillary lymph node (ALN) and tumor size. The diagnostic efficacy of DCE-MRI, IVIM imaging, and their combination for benign and malignant breast lesions was assessed. Correlations between the DCE-MRI and IVIM parameters and prognostic indicators were assessed.

RESULTS

Overall,59 female patients with 62 lesions (22 benign lesions and 40 malignant lesions) were included in this study. The malignant group showed significantly lower D values (p < 0.05) and significantly higher Ktrans, Kep, and f values (p < 0.05). The AUC values of DCE, IVIM, DCE + IVIM were 0.828, 0.882, 0.901. Ktrans, Kep, D and f values were correlated with the pathological grade (p < 0.05); Ktrans was negatively correlated with ER expression (r = -0.519, p < 0.05); Kep was correlated with PR expression and the Ki-67 index (r = -0.489, 0.330, p < 0.05); the DCE and IVIM parameters showed no significant correlations with the HER2 and ALN (p > 0.05). Tumor diameter was correlated with the Kep, D and f values (r = 0.246, -0.278, 0.293; p < 0.05).

CONCLUSION

IVIM and DCE-MRI allowed differential diagnosis of benign and malignant breast lesions, and their combination showed significantly better diagnostic efficiency. DCE- and IVIM-derived parameters showed correlations with some prognostic factors for breast cancer.

Is There an Added Value of Quantitative DCE-MRI by Magnetic Resonance Dispersion Imaging for Prostate Cancer Diagnosis?

In this multicenter, retrospective study, we evaluated the added value of magnetic resonance dispersion imaging (MRDI) to standard multiparametric MRI (mpMRI) for PCa detection. The study included 76 patients, including 51 with clinically significant prostate cancer (csPCa), who underwent radical prostatectomy and had an mpMRI including dynamic contrast-enhanced MRI. Two radiologists performed three separate randomized scorings based on mpMRI, MRDI and mpMRI+MRDI. Radical prostatectomy histopathology was used as the reference standard. Imaging and histopathology were both scored according to the Prostate Imaging-Reporting and Data System V2.0 sector map. Sensitivity and specificity for PCa detection were evaluated for mpMRI, MRDI and mpMRI+MRDI. Inter- and intra-observer variability for both radiologists was evaluated using Cohen's Kappa. On a per-patient level, sensitivity for csPCa for radiologist 1 (R1) for mpMRI, MRDI and mpMRI+MRDI was 0.94, 0.82 and 0.94, respectively. For the second radiologist (R2), these were 0.78, 0.94 and 0.96. R1 detected 4% additional csPCa cases using MRDI compared to mpMRI, and R2 detected 20% extra csPCa cases using MRDI. Inter-observer agreement was significant only for MRDI (Cohen's Kappa = 0.4250, p = 0.004). The results of this study show the potential of MRDI to improve inter-observer variability and the detection of csPCa.

Dynamic Contrast-Enhanced (DCE) MRI

The non-invasive dynamic contrast-enhanced MRI (DCE-MRI) method provides valuable insights into tissue perfusion and vascularity. Primarily used in oncology, DCE-MRI is typically utilized to assess morphology and contrast agent (CA) kinetics in the tissue of interest. Interpretation of the temporal signatures of DCE-MRI data includes qualitative, semi-quantitative, and quantitative approaches. Recent advances in MRI technology allow simultaneous high spatial and temporal resolutions in DCE-MRI data acquisition on most vendor platforms, enabling the more desirable approach of quantitative data analysis using pharmacokinetic (PK) modeling. Many technical factors, including signal-to-noise ratio, temporal resolution, quantifications of arterial input function and native tissue T1, and PK model selection, need to be carefully considered when performing quantitative DCE-MRI. Standardization in data acquisition and analysis is especially important in multi-center studies.

Quantitative dynamic contrast-enhanced magnetic resonance imaging in head and neck cancer: A systematic comparison of different modelling approaches

Background and purpose

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) describes tissue microvasculature and has prognostic and predictive potential in radiotherapy for head and neck cancer (HNC). However, lack in standardization of DCE-MRI hinders comparison of studies and clinical implementation. This study investigated the accuracy and robustness of the population arterial input function (AIF), correlations between pharmacokinetic parameters and their association to T stage and human papillomavirus (HPV) status for HNC.

Materials and methods

DCE-MRI was acquired for 44 HNC patients. Population AIFs were calculated with six different approaches. DCE-MRI was analysed in primary and lymph node tumours using Tofts model (TM) with population AIFs and individual AIFs, extended TM (ETM) with individual AIFs, Brix model (BM), and areas under the curve (AUCs). Intraclass correlation, concordance correlation, Pearson correlation and Whitney Mann U test helped examining the robustness and accuracy of population AIF, correlations between DCE-MRI parameters and their association to T stage and HPV status, respectively.

Results

The population AIF was robust but differed from individual AIFs. There was significant correlation between KtransTM/ETM and ve, TM/ETM, and KtransTM/ETM and Kep, TM/ETM. ABrix and AUCs correlated for lymph nodes. Kep, Brix correlated with ABrix, KtransTM/ETM and Kep, TM/ETM for primary tumours. Kep, TM significantly decreased with increasing T stage. Both the correlations and the parameters' association to T stage were stronger for HPV negative lesions.

Conclusions

Individual AIF was preferred for accurate pharmacokinetic modelling of DCE-MRI. DCE-MRI parameters and their correlations were affected by the lesion type, HPV status and T staging.

Dynamic Contrast-Enhanced Study in the mpMRI of the Prostate-Unnecessary or Underutilised? A Narrative Review

The aim of this review is to summarise recent scientific literature regarding the clinical use of DCE-MRI as a component of multiparametric resonance imaging of the prostate. This review presents the principles of DCE-MRI acquisition and analysis, the current role of DCE-MRI in clinical practice with special regard to its role in presently available categorisation systems, and an overview of the advantages and disadvantages of DCE-MRI described in the current literature. DCE-MRI is an important functional sequence that requires intravenous administration of a gadolinium-based contrast agent and gives information regarding the vascularity and capillary permeability of the lesion. Although numerous studies have confirmed that DCE-MRI has great potential in the diagnosis and monitoring of prostate cancer, its role is still inadequate in the PI-RADS categorisation. Moreover, there have been numerous scientific discussions about abandoning the intravenous application of gadolinium-based contrast as a routine part of MRI examination of the prostate. In this review, we summarised the recent literature on the advantages and disadvantages of DCE-MRI, focusing on an overview of currently available data on bpMRI and mpMRI, as well as on studies providing information on the potential better usability of DCE-MRI in improving the sensitivity and specificity of mpMRI examinations of the prostate.

Clinical-Genomic Risk Group Classification of Suspicious Lesions on Prostate Multiparametric-MRI

The utilization of multi-parametric MRI (mpMRI) in clinical decisions regarding prostate cancer patients' management has recently increased. After biopsy, clinicians can assess risk using National Comprehensive Cancer Network (NCCN) risk stratification schema and commercially available genomic classifiers, such as Decipher. We built radiomics-based models to predict lesions/patients at low risk prior to biopsy based on an established three-tier clinical-genomic classification system. Radiomic features were extracted from regions of positive biopsies and Normally Appearing Tissues (NAT) on T2-weighted and Diffusion-weighted Imaging. Using only clinical information available prior to biopsy, five models for predicting low-risk lesions/patients were evaluated, based on: 1: Clinical variables; 2: Lesion-based radiomic features; 3: Lesion and NAT radiomics; 4: Clinical and lesion-based radiomics; and 5: Clinical, lesion and NAT radiomic features. Eighty-three mpMRI exams from 78 men were analyzed. Models 1 and 2 performed similarly (Area under the receiver operating characteristic curve were 0.835 and 0.838, respectively), but radiomics significantly improved the lesion-based performance of the model in a subset analysis of patients with a negative Digital Rectal Exam (DRE). Adding normal tissue radiomics significantly improved the performance in all cases. Similar patterns were observed on patient-level models. To the best of our knowledge, this is the first study to demonstrate that machine learning radiomics-based models can predict patients' risk using combined clinical-genomic classification.

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.

Improving measurement of blood-brain barrier permeability with reduced scan time using deep-learning-derived capillary input function

PURPOSE

In Dynamic contrast-enhanced MRI (DCE-MRI), Arterial Input Function (AIF) has been shown to be a significant contributor to uncertainty in the estimation of kinetic parameters. This study is to assess the feasibility of using a deep learning network to estimate local Capillary Input Function (CIF) to estimate blood-brain barrier (BBB) permeability, while reducing the required scan time.

MATERIALS AND METHOD

A total of 13 healthy subjects (younger (<40 y/o): 8, older (> 67 y/o): 5) were recruited and underwent 25-min DCE-MRI scans. The 25 min data were retrospectively truncated to 10 min to simulate a reduced scan time of 10 min. A deep learning network was trained to predict the CIF using simulated tissue contrast dynamics with two vascular transport models. The BBB permeability (PS) was measured using 3 methods: (i) Ca-25min, using DCE-MRI data of 25 min with individually sampled AIF (Ca); (ii) Ca-10min, using truncated 10min data with AIF (Ca); and (iii) Cp-10min, using truncated 10 min data with CIF (Cp). The PS estimates from the Ca-25min method were used as reference standard values to assess the accuracy of the Ca-10min and Cp-10min methods in estimating the PS values.

RESULTS

When compared to the reference method(Ca-25min), the Ca-10min and Cp-10min methods resulted in an overestimation of PS by 217 ± 241 % and 48.0 ± 30.2 %, respectively. The Bland Altman analysis showed that the mean difference from the reference was 8.85 ± 1.78 (x10-4 min-1) with the Ca-10min, while it was reduced to 1.63 ± 2.25 (x10-4 min-1) with the Cp-10min, resulting in an average reduction of 81%. The limits of agreement also reduced by up to 39.2% with the Cp-10min. We found a 75% increase of BBB permeability in the gray matter and a 35% increase in the white matter, when comparing the older group to the younger group.

CONCLUSIONS

We demonstrated the feasibility of estimating the capillary-level input functions using a deep learning network. We also showed that this method can be used to estimate subtle age-related changes in BBB permeability with reduced scan time, without compromising accuracy. Moreover, the trained deep learning network can automatically select CIF, reducing the potential uncertainty resulting from manual user-intervention.

Stability of Multi-Parametric Prostate MRI Radiomic Features to Variations in Segmentation

Stability analysis remains a fundamental step in developing a successful imaging biomarker to personalize oncological strategies. This study proposes an in silico contour generation method for simulating segmentation variations to identify stable radiomic features. Ground-truth annotation provided for the whole prostate gland on the multi-parametric MRI sequences (T2w, ADC, and SUB-DCE) were perturbed to mimic segmentation differences observed among human annotators. In total, we generated 15 synthetic contours for a given image-segmentation pair. One thousand two hundred twenty-four unfiltered/filtered radiomic features were extracted applying Pyradiomics, followed by stability assessment using ICC(1,1). Stable features identified in the internal population were then compared with an external population to discover and report robust features. Finally, we also investigated the impact of a wide range of filtering strategies on the stability of features. The percentage of unfiltered (filtered) features that remained robust subjected to segmentation variations were T2w-36% (81%), ADC-36% (94%), and SUB-43% (93%). Our findings suggest that segmentation variations can significantly impact radiomic feature stability but can be mitigated by including pre-filtering strategies as part of the feature extraction pipeline.

Impact of Arterial Input Function and Pharmacokinetic Models on DCE-MRI Biomarkers for Detection of Vascular Effect Induced by Stroma-Directed Drug in an Orthotopic Mouse Model of Pancreatic Cancer

PURPOSE

We demonstrated earlier in mouse models of pancreatic ductal adenocarcinoma (PDA) that K^trans derived from dynamic contrast-enhanced (DCE) MRI detected microvascular effect induced by PEGPH20, a hyaluronidase which removes stromal hyaluronan, leading to reduced interstitial fluid pressure in the tumor (Clinical Cancer Res (2019) 25: 2314-2322). How the choice of pharmacokinetic (PK) model and arterial input function (AIF) may impact DCE-derived markers for detecting such an effect is not known.

PROCEDURES

Retrospective analyses of the DCE-MRI of the orthotopic PDA model are performed to examine the impact of individual versus group AIF combined with Tofts model (TM), extended-Tofts model (ETM), or shutter-speed model (SSM) on the ability to detect the microvascular changes induced by PEGPH20 treatment.

RESULTS

Individual AIF exhibit a marked difference in peak gadolinium concentration. However, across all three PK models, k_ep values show a significant correlation between individual versus group-AIF (p < 0.01). Regardless individual or group AIF, when k_ep is obtained from fitting the DCE-MRI data using the SSM, k_ep shows a significant increase after PEGPH20 treatment (p < 0.05 compared to the baseline); %change of k_ep from baseline to post-treatment is also significantly different between PEGPH20 versus vehicle group (p < 0.05). In comparison, when k_ep is derived from the TM, only the use of individual AIF leads to a significant increase of k_ep after PEGPH20 treatment, whereas the %change of k_ep is not different between PEGPH20 versus vehicle group. Group AIF but not individual AIF allows detection of a significant increase of V_p (derived from the ETM) in PEGPH20 versus vehicle group (p < 0.05). Increase of V_p is consistent with a large increase of mean capillary lumen area estimated from immunostaining.

CONCLUSION

Our results suggest that k_ep derived from SSM and V_p from ETM, both using group AIF, are optimal for the detection of microvascular changes induced by stroma-directed drug PEGPH20. These analyses provide insights in the choice of PK model and AIF for optimal DCE protocol design in mouse pancreatic cancer models.

Semi-quantitative and quantitative dynamic contrast-enhanced (DCE) MRI parameters as prostate cancer imaging biomarkers for biologically targeted radiation therapy

BACKGROUND

Biologically targeted radiation therapy treatment planning requires voxel-wise characterisation of tumours. Dynamic contrast enhanced (DCE) DCE MRI has shown promise in defining voxel-level biological characteristics. In this study we consider the relative value of qualitative, semi-quantitative and quantitative assessment of DCE MRI compared with diffusion weighted imaging (DWI) and T2-weighted (T2w) imaging to detect prostate cancer at the voxel level.

METHODS

Seventy prostate cancer patients had multiparametric MRI prior to radical prostatectomy, including T2w, DWI and DCE MRI. Apparent Diffusion Coefficient (ADC) maps were computed from DWI, and semi-quantitative and quantitative parameters computed from DCE MRI. Tumour location and grade were validated with co-registered whole mount histology. Kolmogorov-Smirnov tests were applied to determine whether MRI parameters in tumour and benign voxels were significantly different. Cohen's d was computed to quantify the most promising biomarkers. The Parker and Weinmann Arterial Input Functions (AIF) were compared for their ability to best discriminate between tumour and benign tissue. Classifier models were used to determine whether DCE MRI parameters improved tumour detection versus ADC and T2w alone.

RESULTS

All MRI parameters had significantly different data distributions in tumour and benign voxels. For low grade tumours, semi-quantitative DCE MRI parameter time-to-peak (TTP) was the most discriminating and outperformed ADC. For high grade tumours, ADC was the most discriminating followed by DCE MRI parameters Ktrans, the initial rate of enhancement (IRE), then TTP. Quantitative parameters utilising the Parker AIF better distinguished tumour and benign voxel values than the Weinmann AIF. Classifier models including DCE parameters versus T2w and ADC alone, gave detection accuracies of 78% versus 58% for low grade tumours and 85% versus 72% for high grade tumours.

CONCLUSIONS

Incorporating DCE MRI parameters with DWI and T2w gives improved accuracy for tumour detection at a voxel level. DCE MRI parameters should be used to spatially characterise tumour biology for biologically targeted radiation therapy treatment planning.

Effectiveness of Dynamic Contrast Enhanced MRI with a Split Dose of Gadoterate Meglumine for Detection of Prostate Cancer

RATIONALE AND OBJECTIVES

To evaluate whether dynamic contrast enhanced (DCE) MRI with a split injection of 30% followed by 70% of a standard dose (30PSD and 70PSD) of gadoterate meglumine (DOTAREM) can improve diagnosis of prostate cancer (PCa).

MATERIALS AND METHODS

MRI for twenty patients was performed on a Philips Ingenia 3T scanner without an endorectal coil followed by subsequent radical prostatectomy. DCE 3D T1-FFE data were acquired with injection of 0.03 mmol/kg followed after 2 minutes by 0.07 mmol/kg of DOTAREM. Regions-of-interest on histologically verified PCa and normal tissue in different prostate zones and the iliac artery were drawn. Average signal intensity as function of time was calculated for each ROI and fitted by using the signal intensity form of the Tofts (SI-Tofts) model to extract physiological parameters (K^trans and v_e). In addition, the scaled arterial input function (AIF) obtained from 30PSD data was used to analyze 70PSD data.

RESULTS

The AIF obtained from 30PSD data showed both first and second passes clearly and had much higher peak magnitude than AIFs from 70PSD data. K^trans was significantly (p < 0.05) larger in PCa than in normal tissue in peripheral zone (PZ) and central zone (CZ) for both 70PSD and 70PSD data analyzed with a scaled AIF. K^trans in cancer overlapped with that of normal tissue in the transition zone (TZ). There was no statistical difference in v_e between cancer and normal tissue. Receiver operating characteristic analysis showed that use of the AIF from 30PSD data to analyze 70PSD data increased the diagnostic efficacy of K^trans in the PZ and CZ.

CONCLUSION

The split dose protocol for injection of Dotarem increased diagnostic accuracy of quantitative analysis with the SI-Tofts model.

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.

Estimation of the capillary level input function for dynamic contrast-enhanced MRI of the breast using a deep learning approach

PURPOSE

To develop a deep learning approach to estimate the local capillary-level input function (CIF) for pharmacokinetic model analysis of DCE-MRI.

METHODS

A deep convolutional network was trained with numerically simulated data to estimate the CIF. The trained network was tested using simulated lesion data and used to estimate voxel-wise CIF for pharmacokinetic model analysis of breast DCE-MRI data using an abbreviated protocol from women with malignant (n = 25) and benign (n = 28) lesions. The estimated parameters were used to build a logistic regression model to detect the malignancy.

RESULT

The pharmacokinetic parameters estimated using the network-predicted CIF from our breast DCE data showed significant differences between the malignant and benign groups for all parameters. Testing the diagnostic performance with the estimated parameters, the conventional approach with arterial input function (AIF) showed an area under the curve (AUC) between 0.76 and 0.87, and the proposed approach with CIF demonstrated similar performance with an AUC between 0.79 and 0.81.

CONCLUSION

This study shows the feasibility of estimating voxel-wise CIF using a deep neural network. The proposed approach could eliminate the need to measure AIF manually without compromising the diagnostic performance to detect the malignancy in the clinical setting.

Estimation of Contrast Agent Concentration in DCE-MRI Using 2 Flip Angles

PURPOSE

The aim of this study was to investigate the feasibility of using 2 flip angles (FAs) with an ultrashort echo time during dynamic contrast-enhanced (DCE)-magnetic resonance imaging (MRI) for estimation of plasma gadolinium (Gd) concentration without using a precontrast longitudinal relaxation time T1 (T10) measurement.

METHODS

T1-weighted DCE-MRI experiments were carried out with C57BL/6J mice using the scan protocol with 2 FAs over 3 sequential segments during 1 scan. The data with 2 FAs were used to estimate T10 (T1T) during conversion of a time-intensity curve to the time-concentration curve. Three dosages of gadolinium-based contrast agent were used to achieve a wide range of variability in Gd concentrations when measured at 10 minutes postinjection: 0.05 mmol/kg (n = 6), 0.1 mmol/kg (n = 11), and 0.15 mmol/kg (n = 7). For comparison, the signal-to-concentration conversion was also conducted using the T10 measured from the precontrast scan (T1M) as well as a constant T10 (2.1 seconds) from the literature (T1C). The Gd concentrations ([Gd]) estimated using DCE-MRI data for the time of retro-orbital blood collection ([Gd]T1T, [Gd]T1M, and [Gd]T1C, respectively) were compared against the [Gd] of the blood samples measured by inductively coupled plasma mass spectrometry ([Gd]MS). In addition, contrast kinetic model analysis was conducted on mice with GL261 brain tumors (n = 5) using the 3 different methods for T10.

RESULTS

T1T strongly correlated with T1M (r = 0.81). [Gd]T1M and [Gd]T1T were significantly different from [Gd]T1C. [Gd]T1M and [Gd]T1T were in good agreement with [Gd]MS with strong correlations (mean percentage error ± standard deviation) of r = 0.70 (16% ± 56%) and r = 0.85 (15% ± 44%), respectively. In contrast, [Gd]T1C had a weak correlation of r = 0.52 with larger errors of 33% ± 24%. The contrast kinetic model parameters of GL261 brain tumors using T1T were not significantly different from those using T1M.

CONCLUSIONS

This study substantiates the feasibility of using the 2-FA approach during DCE-MRI scan to estimate [Gd] in the plasma without using an extra scan to perform precontrast T1 measurements.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

DCE-MRI of Brain Fluid Barriers: In Vivo Water Cycling at the Human Choroid Plexus

Metabolic deficits at brain-fluid barriers are an increasingly recognized feature of cognitive decline in older adults. At the blood-cerebrospinal fluid barrier, water is transported across the choroid plexus (CP) epithelium against large osmotic gradients via processes tightly coupled to activity of the sodium/potassium pump. Here, we quantify CP homeostatic water exchange using dynamic contrast-enhanced MRI and investigate the association of the water efflux rate constant (k_co) with cognitive dysfunction in older individuals. Temporal changes in the longitudinal relaxation rate constant (R₁) after contrast agent bolus injection were measured in a CP region of interest in 11 participants with mild cognitive dysfunction [CI; 73 ± 6 years] and 28 healthy controls [CN; 72 ± 7 years]. k_co was determined from a modified two-site pharmacokinetic exchange analysis of the R₁ time-course. K^trans, a measure of contrast agent extravasation to the interstitial space was also determined. Cognitive function was assessed by neuropsychological test performance. k_co averages 5.8 ± 2.7 s^-1 in CN individuals and is reduced by 2.4 s^-1 [ca. 40%] in CI subjects. Significant associations of k_co with global cognition and multiple cognitive domains are observed. K^trans averages 0.13 ± 0.07 min^-1 and declines with age [-0.006 ± 0.002 min^-1 yr^-1], but shows no difference between CI and CN individuals or association with cognitive performance. Our findings suggest that the CP water efflux rate constant is associated with cognitive dysfunction and shows an age-related decline in later life, consistent with the metabolic disturbances that characterize brain aging.

Multi-parametric MRI (mpMRI) for treatment response assessment of radiation therapy

Magnetic resonance imaging (MRI) plays an important role in the modern radiation therapy (RT) workflow. In comparison with computed tomography (CT) imaging, which is the dominant imaging modality in RT, MRI possesses excellent soft-tissue contrast for radiographic evaluation. Based on quantitative models, MRI can be used to assess tissue functional and physiological information. With the developments of scanner design, acquisition strategy, advanced data analysis, and modeling, multiparametric MRI (mpMRI), a combination of morphologic and functional imaging modalities, has been increasingly adopted for disease detection, localization, and characterization. Integration of mpMRI techniques into RT enriches the opportunities to individualize RT. In particular, RT response assessment using mpMRI allows for accurate characterization of both tissue anatomical and biochemical changes to support decision-making in monotherapy of radiation treatment and/or systematic cancer management. In recent years, accumulating evidence have, indeed, demonstrated the potentials of mpMRI in RT response assessment regarding patient stratification, trial benchmarking, early treatment intervention, and outcome modeling. Clinical application of mpMRI for treatment response assessment in routine radiation oncology workflow, however, is more complex than implementing an additional imaging protocol; mpMRI requires additional focus on optimal study design, practice standardization, and unified statistical reporting strategy to realize its full potential in the context of RT. In this article, the mpMRI theories, including image mechanism, protocol design, and data analysis, will be reviewed with a focus on the radiation oncology field. Representative works will be discussed to demonstrate how mpMRI can be used for RT response assessment. Additionally, issues and limits of current works, as well as challenges and potential future research directions, will also be discussed.

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.

Intrasession Reliability of Arterial Spin-Labeled MRI-Measured Noncontrast Perfusion in Glioblastoma at 3 T

Arterial spin-labeled magnetic resonance imaging can provide quantitative perfusion measurements in the brain and can be potentially used to evaluate therapy response assessment in glioblastoma (GBM). The reliability and reproducibility of this method to measure noncontrast perfusion in GBM, however, are lacking. We evaluated the intrasession reliability of brain and tumor perfusion in both healthy volunteers and patients with GBM at 3 T using pseudocontinuous labeling (pCASL) and 3D turbo spin echo (TSE) using Cartesian acquisition with spiral profile reordering (CASPR). Two healthy volunteers at a single time point and 6 newly diagnosed patients with GBM at multiple time points (before, during, and after chemoradiation) underwent scanning (total, 14 sessions). Compared with 3D GraSE, 3D TSE-CASPR generated cerebral blood flow maps with better tumor-to-normal background tissue contrast and reduced image distortions. The intraclass correlation coefficient between the 2 runs of 3D pCASL with TSE-CASPR was consistently high (≥0.90) across all normal-appearing gray matter (NAGM) regions of interest (ROIs), and was particularly high in tumors (0.98 with 95% confidence interval [CI]: 0.97-0.99). The within-subject coefficients of variation were relatively low in all normal-appearing gray matter regions of interest (3.40%-7.12%), and in tumors (4.91%). Noncontrast perfusion measured using 3D pCASL with TSE-CASPR provided robust cerebral blood flow maps in both healthy volunteers and patients with GBM with high intrasession repeatability at 3 T. This approach can be an appropriate noncontrast and noninvasive quantitative perfusion imaging method for longitudinal assessment of therapy response and management of patients with GBM.

AI applications to medical images: From machine learning to deep learning

PURPOSE

Artificial intelligence (AI) models are playing an increasing role in biomedical research and healthcare services. This review focuses on challenges points to be clarified about how to develop AI applications as clinical decision support systems in the real-world context.

METHODS

A narrative review has been performed including a critical assessment of articles published between 1989 and 2021 that guided challenging sections.

RESULTS

We first illustrate the architectural characteristics of machine learning (ML)/radiomics and deep learning (DL) approaches. For ML/radiomics, the phases of feature selection and of training, validation, and testing are described. DL models are presented as multi-layered artificial/convolutional neural networks, allowing us to directly process images. The data curation section includes technical steps such as image labelling, image annotation (with segmentation as a crucial step in radiomics), data harmonization (enabling compensation for differences in imaging protocols that typically generate noise in non-AI imaging studies) and federated learning. Thereafter, we dedicate specific sections to: sample size calculation, considering multiple testing in AI approaches; procedures for data augmentation to work with limited and unbalanced datasets; and the interpretability of AI models (the so-called black box issue). Pros and cons for choosing ML versus DL to implement AI applications to medical imaging are finally presented in a synoptic way.

CONCLUSIONS

Biomedicine and healthcare systems are one of the most important fields for AI applications and medical imaging is probably the most suitable and promising domain. Clarification of specific challenging points facilitates the development of such systems and their translation to clinical practice.

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

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

DCE-MRI Quantitative Parameters as Predictors of Treatment Response in Patients With Locally Advanced Cervical Squamous Cell Carcinoma Underwent CCRT

Purpose

To evaluate the predictive value of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) quantitative parameters in treatment response to concurrent chemoradiotherapy (CCRT) for locally advanced cervical squamous cell carcinoma (LACSC).

Methods and materials

LACSC patients underwent CCRT had DCE-MRI before (e0) and after 3 days of treatment (e3). Extended Tofts Linear model with a user arterial input function was adopted to generate quantitative measurements. Endothelial transfer constant (K^trans), reflux rate (K_ep), fractional extravascular extracellular space volume (V_e), and fractional plasma volume (V_p) were calculated, and percentage changes ΔK^trans, ΔK_ep, ΔV_e, and ΔV_p were computed. The correlations of these measurements with the tumor regression rate were analyzed. The predictive value of these parameters on treatment outcome was generated by the receiver operating characteristic (ROC) curve. Univariate and multivariate logistic regression analyses were conducted to find the independent variables.

Results

K^trans-e0, K_ep -e0, ΔK^trans, and ΔV_e were positively correlated with the tumor regression rate. Mean values of K^trans-e0, K^trans-e3, ΔK^trans, and ΔV_e were higher in the non-residual tumor group than residual tumor group and were independent prognostic factors for predicting residual tumor occurrence. K^trans-e3 showed the highest area under the curve (AUC) for treatment response prediction.

Conclusions

Quantitative parameters at e0 and e3 from DCE-MRI could be used as potential indicators for predicting treatment response of LACSC.

Dose-response assessment by quantitative MRI in a phase 1 clinical study of the anti-cancer vascular disrupting agent crolibulin

The vascular disrupting agent crolibulin binds to the colchicine binding site and produces anti-vascular and apoptotic effects. In a multisite phase 1 clinical study of crolibulin (NCT00423410), we measured treatment-induced changes in tumor perfusion and water diffusivity (ADC) using dynamic contrast-enhanced MRI (DCE-MRI) and diffusion-weighted MRI (DW-MRI), and computed correlates of crolibulin pharmacokinetics. 11 subjects with advanced solid tumors were imaged by MRI at baseline and 2–3 days post-crolibulin (13–24 mg/m²). ADC maps were computed from DW-MRI. Pre-contrast T₁ maps were computed, co-registered with the DCE-MRI series, and maps of area-under-the-gadolinium-concentration-curve-at-90 s (AUC_90s) and the Extended Tofts Model parameters k^trans, v_e, and v_p were calculated. There was a strong correlation between higher plasma drug \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}^{max}$$\end{document}Cmax and a linear combination of (1) reduction in tumor fraction with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${AUC}_{90s}>15.8$$\end{document}AUC90s>15.8 mM s, and, (2) increase in tumor fraction with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${v}_{e}<0.3$$\end{document}ve<0.3. A higher plasma drug AUC was correlated with a linear combination of (1) increase in tumor fraction with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{ADC}} < 1.1 \times 10^{ - 3} \;{\text{mm}}^{2} /{\text{s}}$$\end{document}ADC<1.1×10-3mm2/s, and, (2) increase in tumor fraction with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_{e}<0.3$$\end{document}ve<0.3. These findings are suggestive of cell swelling and decreased tumor perfusion 2–3 days post-treatment with crolibulin. The multivariable linear regression models reported here can inform crolibulin dosing in future clinical studies of crolibulin combined with cytotoxic or immune-oncology agents.

Efficient DCE-MRI Parameter and Uncertainty Estimation Using a Neural Network

Quantitative DCE-MRI provides voxel-wise estimates of tracer-kinetic parameters that are valuable in the assessment of health and disease. These maps suffer from many known sources of variability. This variability is expensive to compute using current methods, and is typically not reported. Here, we demonstrate a novel approach for simultaneous estimation of tracer-kinetic parameters and their uncertainty due to intrinsic characteristics of the tracer-kinetic model, with very low computation time. We train and use a neural network to estimate the approximate joint posterior distribution of tracer-kinetic parameters. Uncertainties are estimated for each voxel and are specific to the patient, exam, and lesion. We demonstrate the methods' ability to produce accurate tracer-kinetic maps. We compare predicted parameter ranges with uncertainties introduced by noise and by differences in post-processing in a digital reference object. The predicted parameter ranges correlate well with tracer-kinetic parameter ranges observed across different noise realizations and regression algorithms. We also demonstrate the value of this approach to differentiate significant from insignificant changes in brain tumor pharmacokinetics over time. This is achieved by enforcing consistency in resolving model singularities in the applied tracer-kinetic model.

BRAINSTORM: A Multi-Institutional Phase 1/2 Study of RRx-001 in Combination With Whole Brain Radiation Therapy for Patients With Brain Metastases

PURPOSE

To determine the recommended phase 2 dose of RRx-001, a radiosensitizer with vascular normalizing properties, when used with whole-brain radiation therapy (WBRT) for brain metastases and to assess whether quantitative changes in perfusion magnetic resonance imaging (MRI) after RRx-001 correlate with response.

METHODS AND MATERIALS

Five centers participated in this phase 1/2 trial of RRx-001 given once pre-WBRT and then twice weekly during WBRT. Four dose levels were planned (5 mg/m², 8.4 mg/m², 16.5 mg/m², 27.5 mg/m²). Dose escalation was managed by the time-to-event continual reassessment method algorithm. Linear mixed models were used to correlate change in 24-hour T1, Ktrans (capillary permeability), and fractional plasma volume with change in tumor volume.

RESULTS

Between 2015 and 2017, 31 patients were enrolled. Two patients dropped out before any therapy. Median age was 60 years (range, 30-76), and 12 were male. The most common tumor types were melanoma (59%) and non-small cell lung cancer (18%). No dose limiting toxicities were observed. The most common severe adverse event was grade 3 asthenia (6.9%, 2 of 29). The median intracranial response rate was 46% (95% confidence interval, 24-68) and median overall survival was 5.2 months (95% confidence interval, 4.5-9.4). No neurologic deaths occurred. Among 10 patients undergoing dynamic contrast-enhanced MRI, a reduction in Vp 24 hours after RRx-001 was associated with reduced tumor volume at 1 and 4 months (P ≤ .01).

CONCLUSIONS

The addition of RRx-001 to WBRT is well tolerated with favorable intracranial response rates. Because activity was observed across all dose levels, the recommended phase 2 dose is 10 mg twice weekly. A reduction in fractional plasma volume on dynamic contrast-enhanced MRI 24 hours after RRx-001 suggests antiangiogenic activity associated with longer-term tumor response.

Multiparameter MRI Predictors of Long-Term Survival in Glioblastoma Multiforme

Standard-of-care multiparameter magnetic resonance imaging (MRI) scans of the brain were used to objectively subdivide glioblastoma multiforme (GBM) tumors into regions that correspond to variations in blood flow, interstitial edema, and cellular density. We hypothesized that the distribution of these distinct tumor ecological "habitats" at the time of presentation will impact the course of the disease. We retrospectively analyzed initial MRI scans in 2 groups of patients diagnosed with GBM, a long-term survival group comprising subjects who survived >36 month postdiagnosis, and a short-term survival group comprising subjects who survived ≤19 month postdiagnosis. The single-institution discovery cohort contained 22 subjects in each group, while the multi-institution validation cohort contained 15 subjects per group. MRI voxel intensities were calibrated, and tumor voxels clustered on contrast-enhanced T1-weighted and fluid-attenuated inversion-recovery (FLAIR) images into 6 distinct "habitats" based on low- to medium- to high-contrast enhancement and low-high signal on FLAIR scans. Habitat 6 (high signal on calibrated contrast-enhanced T1-weighted and FLAIR sequences) comprised a significantly higher volume fraction of tumors in the long-term survival group (discovery cohort, 35% ± 6.5%; validation cohort, 34% ± 4.8%) compared with tumors in the short-term survival group (discovery cohort, 17% ± 4.5%, P < .03; validation cohort, 16 ± 4.0%, P < .007). Of the 6 distinct MRI-defined habitats, the fractional tumor volume of habitat 6 at diagnosis was significantly predictive of long- or short-term survival. We discuss a possible mechanistic basis for this association and implications for habitat-driven adaptive therapy of GBM.

Real-Time Quantitative Assessment of Accuracy and Precision of Blood Volume Derived from DCE-MRI in Individual Patients During a Clinical Trial

Accuracy and precision of quantitative imaging (QI) metrics should be assessed in real time in each patient during a clinical trial to support QI-based decision-making. We developed a framework for real-time quantitative assessment of QI metrics and evaluated accuracy and precision of dynamic contrast-enhanced (DCE)-magnetic resonance imaging (MRI)–derived blood volume (BV) in a clinical trial for head and neck cancers. Patients underwent DCE-MRI before and after 2 weeks of radiation therapy (2wkRT). A mean as a reference value and a repeatability coefficient (RC) of BV values established from n patients in cerebellum volumes of interest (VOIs), which were normal and affected little by therapy, served as accuracy and precision measurements. The BV maps of a new patient were called accurate and precise if the values in cerebellum VOIs and the difference between the 2 scans agreed with the respective mean and RC with 95% confidence. The new data could be used to update reference values. Otherwise, the data were flagged for further evaluation before use in the trial. BV maps from 62 patients enrolled on the trial were evaluated. Mean BV values were 2.21 (±0.14) mL/100 g pre-RT and 2.22 (±0.17) mL/100 g at 2wkRT; relative RC was 15.9%. The BV maps from 3 patients were identified to be inaccurate and imprecise before use in the clinical trial. Our framework of real-time quantitative assessment of QI metrics during a clinical trial can be translated to different QI metrics and organ-sites for supporting QI-based decision-making that warrants success of a clinical trial.

Developing a Pipeline for Multiparametric MRI-Guided Radiation Therapy: Initial Results from a Phase II Clinical Trial in Newly Diagnosed Glioblastoma

Quantitative mapping of hyperperfused and hypercellular regions of glioblastoma has been proposed to improve definition of tumor regions at risk for local recurrence following conventional radiation therapy. As the processing of the multiparametric dynamic contrast-enhanced (DCE-) and diffusion-weighted (DW-) magnetic resonance imaging (MRI) data for delineation of these subvolumes requires additional steps that go beyond the standard practices of target definition, we sought to devise a workflow to support the timely planning and treatment of patients. A phase II study implementing a multiparametric imaging biomarker for tumor hyperperfusion and hypercellularity consisting of DCE-MRI and high b-value DW-MRI to guide intensified (75 Gy/30 fractions) radiation therapy (RT) in patients with newly diagnosed glioblastoma was launched. In this report, the workflow and the initial imaging outcomes of the first 12 patients are described. Among all the first 12 patients, treatment was initiated within 6 weeks of surgery and within 2 weeks of simulation. On average, the combined hypercellular volume and high cerebral blood volume/tumor perfusion volume were 1.8 times smaller than the T1 gadolinium abnormality and 10 times smaller than the FLAIR abnormality. Hypercellular volume and high cerebral blood volume/tumor perfusion volume each identified largely distinct regions and showed 57% overlap with the enhancing abnormality, and minimal-to-no extension outside of the FLAIR. These results show the feasibility of implementing a workflow for multiparametric magnetic resonance-guided radiation therapy into clinical trials with a coordinated multidisciplinary team, and the unique and complementary tumor subregions identified by the combination of high b-value DW-MRI and DCE-MRI.

QIN Benchmarks for Clinical Translation of Quantitative Imaging Tools

The Quantitative Imaging Network of the National Cancer Institute is in its 10th year of operation, and research teams within the network are developing and validating clinical decision support software tools to measure or predict the response of cancers to various therapies. As projects progress from development activities to validation of quantitative imaging tools and methods, it is important to evaluate the performance and clinical readiness of the tools before committing them to prospective clinical trials. A variety of tests, including special challenges and tool benchmarking, have been instituted within the network to prepare the quantitative imaging tools for service in clinical trials. This article highlights the benchmarking process and provides a current evaluation of several tools in their transition from development to validation.

Advances in Diffusion and Perfusion MRI for Quantitative Cancer Imaging

PURPOSE OF REVIEW

This article is to review recent technical developments and their clinical applications in cancer imaging quantitative measurement of cellular and vascular properties of the tumors.

RECENT FINDINGS

Rapid development of fast Magnetic Resonance Imaging (MRI) technologies over last decade brought new opportunities in quantitative MRI methods to measure both cellular and vascular properties of tumors simultaneously.

SUMMARY

Diffusion MRI (dMRI) and dynamic contrast enhanced (DCE)-MRI have become widely used to assess the tissue structural and vascular properties, respectively. However, the ultimate potential of these advanced imaging modalities has not been fully exploited. The dependency of dMRI on the diffusion weighting gradient strength and diffusion time can be utilized to measure tumor perfusion, cellular structure, and cellular membrane permeability. Similarly, DCE-MRI can be used to measure vascular and cellular membrane permeability along with cellular compartment volume fractions. To facilitate the understanding of these potentially important methods for quantitative cancer imaging, we discuss the basic concepts and recent developments, as well as future directions for further development.

Impact of (k,t) sampling on DCE MRI tracer kinetic parameter estimation in digital reference objects

Purpose

To evaluate the impact of (k,t) data sampling on the variance of tracer‐kinetic parameter (TK) estimation in high‐resolution whole‐brain dynamic contrast enhanced magnetic resonance imaging (DCE‐MRI) using digital reference objects. We study this in the context of TK model constraints, and in the absence of other constraints.

Methods

Three anatomically and physiologically realistic brain‐tumor digital reference objects were generated. Data sampling strategies included uniform and variable density; zone‐based, lattice, pseudo‐random, and pseudo‐radial; with 50‐time frames and 4‐fold to 25‐fold undersampling. In all cases, we assume a fully sampled first time frame, and prior knowledge of the arterial input function. TK parameters were estimated by indirect estimation (i.e., image‐time‐series reconstruction followed by model fitting), and direct estimation from the under‐sampled data. We evaluated methods based on the Cramér‐Rao bound and Monte‐Carlo simulations, over the range of signal‐to‐noise ratio (SNR) seen in clinical brain DCE‐MRI.

Results

Lattice‐based sampling provided the lowest SDs, followed by pseudo‐random, pseudo‐radial, and zone‐based. This ranking was consistent for the Patlak and extended Tofts model. Pseudo‐random sampling resulted in 19% higher averaged SD compared to lattice‐based sampling. Zone‐based sampling resulted in substantially higher SD at undersampling factors above 10. CRB analysis showed only a small difference between uniform and variable density for both lattice‐based and pseudo‐random sampling up to undersampling factors of 25.

Conclusion

Lattice sampling provided the lowest SDs, although the differences between sampling schemes were not substantial at low undersampling factors. The differences between lattice‐based and pseudo‐random sampling strategies with both uniform and variable density were within the range of error induced by other sources, at up to 25‐fold undersampling.

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.

Pharmacokinetic modeling of dynamic contrast-enhanced MRI using a reference region and input function tail

PURPOSE

Quantitative analysis of dynamic contrast-enhanced MRI (DCE-MRI) requires an arterial input function (AIF) which is difficult to measure. We propose the reference region and input function tail (RRIFT) approach which uses a reference tissue and the washout portion of the AIF.

METHODS

RRIFT was evaluated in simulations with 100 parameter combinations at various temporal resolutions (5-30 s) and noise levels (σ = 0.01-0.05 mM). RRIFT was compared against the extended Tofts model (ETM) in 8 studies from patients with glioblastoma multiforme. Two versions of RRIFT were evaluated: one using measured patient-specific AIF tails, and another assuming a literature-based AIF tail.

RESULTS

RRIFT estimated the transfer constant K trans and interstitial volume v e with median errors within 20% across all simulations. RRIFT was more accurate and precise than the ETM at temporal resolutions slower than 10 s. The percentage error of K trans had a median and interquartile range of -9 ± 45% with the ETM and -2 ± 17% with RRIFT at a temporal resolution of 30 s under noiseless conditions. RRIFT was in excellent agreement with the ETM in vivo, with concordance correlation coefficients (CCC) of 0.95 for K trans , 0.96 for v e , and 0.73 for the plasma volume v p using a measured AIF tail. With the literature-based AIF tail, the CCC was 0.89 for K trans , 0.93 for v e and 0.78 for v p .

CONCLUSIONS

Quantitative DCE-MRI analysis using the input function tail and a reference tissue yields absolute kinetic parameters with the RRIFT method. This approach was viable in simulation and in vivo for temporal resolutions as low as 30 s.

Optimal mass transport kinetic modeling for head and neck DCE-MRI: Initial analysis

PURPOSE

Current state-of-the-art models for estimating the pharmacokinetic parameters do not account for intervoxel movement of the contrast agent (CA). We introduce an optimal mass transport (OMT) formulation that naturally handles intervoxel CA movement and distinguishes between advective and diffusive flows.

METHOD

Ten patients with head and neck squamous cell carcinoma (HNSCC) were enrolled in the study between June 2014 and October 2015 and underwent DCE MRI imaging prior to beginning treatment. The CA tissue concentration information was taken as the input in the data-driven OMT model. The OMT approach was tested on HNSCC DCE data that provides quantitative information for forward flux ( Φ F ) and backward flux ( Φ B ). OMT-derived Φ F was compared with the volume transfer constant for CA, K trans , derived from the Extended Tofts Model (ETM).

RESULTS

The OMT-derived flows showed a consistent jump in the CA diffusive behavior across the images in accordance with the known CA dynamics. The mean forward flux was 0.0082 ± 0.0091 ( min - 1 ) whereas the mean advective component was 0.0052 ± 0.0086 ( min - 1 ) in the HNSCC patients. The diffusive percentages in forward and backward flux ranged from 8.67% to 18.76% and 12.76% to 30.36%, respectively. The OMT model accounts for intervoxel CA movement and results show that the forward flux ( Φ F ) is comparable with the ETM-derived K trans .

CONCLUSIONS

This is a novel data-driven study based on optimal mass transport principles applied to patient DCE imaging to analyze CA flow in HNSCC.

Repeatability of Multiparametric Prostate MRI Radiomics Features

In this study we assessed the repeatability of radiomics features on small prostate tumors using test-retest Multiparametric Magnetic Resonance Imaging (mpMRI). The premise of radiomics is that quantitative image-based features can serve as biomarkers for detecting and characterizing disease. For such biomarkers to be useful, repeatability is a basic requirement, meaning its value must remain stable between two scans, if the conditions remain stable. We investigated repeatability of radiomics features under various preprocessing and extraction configurations including various image normalization schemes, different image pre-filtering, and different bin widths for image discretization. Although we found many radiomics features and preprocessing combinations with high repeatability (Intraclass Correlation Coefficient > 0.85), our results indicate that overall the repeatability is highly sensitive to the processing parameters. Neither image normalization, using a variety of approaches, nor the use of pre-filtering options resulted in consistent improvements in repeatability. We urge caution when interpreting radiomics features and advise paying close attention to the processing configuration details of reported results. Furthermore, we advocate reporting all processing details in radiomics studies and strongly recommend the use of open source implementations.

Use of Indicator Dilution Principle to Evaluate Accuracy of Arterial Input Function Measured With Low-Dose Ultrafast Prostate Dynamic Contrast-Enhanced MRI

Accurately measuring arterial input function (AIF) is essential for quantitative analysis of dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI). We used the indicator dilution principle to evaluate the accuracy of AIF measured directly from an artery following a low-dose contrast media ultrafast DCE-MRI. In total, 15 patients with biopsy-confirmed localized prostate cancers were recruited. Cardiac MRI (CMRI) and ultrafast DCE-MRI were acquired on a Philips 3 T Ingenia scanner. The AIF was measured at iliac arties following injection of a low-dose (0.015 mmol/kg) gadolinium (Gd) contrast media. The cardiac output (CO) from CMRI (CO_CMRI) was calculated from the difference in ventricular volume at diastole and systole measured on the short axis of heart. The CO from DCE-MRI (CO_DCE) was also calculated from the AIF and dose of the contrast media used. A correlation test and Bland–Altman plot were used to compare CO_CMRI and CO_DCE. The average (±standard deviation [SD]) area under the curve measured directly from local AIF was 0.219 ± 0.07 mM·min. The average (±SD) CO_CMRI and CO_DCE were 6.52 ± 1.47 L/min and 6.88 ± 1.64 L/min, respectively. There was a strong positive correlation (r = 0.82, P < .01) and good agreement between CO_CMRI and CO_DCE. The CO_DCE is consistent with the reference standard CO_CMRI. This indicates that the AIF can be measured accurately from an artery with ultrafast DCE-MRI following injection of a low-dose contrast media.

Dynamic Contrast-enhanced MRI Detects Responses to Stroma-directed Therapy in Mouse Models of Pancreatic Ductal Adenocarcinoma

PURPOSE

The dense stroma underlies the drug resistance of pancreatic ductal adenocarcinoma (PDA) and has motivated the development of stroma-directed drugs. Our objective is to test the concept that dynamic contrast-enhanced (DCE) MRI using FDA-approved contrast media, an imaging method sensitive to the tumor microenvironment, can detect early responses to stroma-directed drug.

EXPERIMENTAL DESIGN

Imaging studies were performed in three mouse models exhibiting high desmoplastic reactions: the autochthonous PDA in genetically engineered mice (KPC), an orthotopic model in syngeneic mice, and a xenograft model of human PDA in athymic mice. An investigational drug, PEGPH20 (pegvorhyaluronidase alfa), which degrades hyaluronan (HA) in the stroma of PDA, was injected alone or in combination with gemcitabine.

RESULTS

At 24 hours after a single injection of PEGPH20, Ktrans , a DCE-MRI-derived marker that measures how fast a unit volume of contrast media is transferred from capillaries to interstitial space, increased 56% and 50% from baseline in the orthotopic and xenograft tumors, respectively, compared with a 4% and 6% decrease in vehicle groups (both P < 0.05). Similarly, after three combined treatments, Ktrans in KPC mice increased 54%, whereas it decreased 4% in controls treated with gemcitabine alone (P < 0.05). Consistently, after a single injection of PEGPH20, tumor HA content assessed by IHC was reduced substantially in all three models while drug delivery (measured by paclitaxel accumulation in tumor) was increased by 2.6-fold.

CONCLUSIONS

These data demonstrated a DCE-MRI marker, Ktrans , can detect early responses to stroma-directed drug and reveal the sustained effect of combination treatment (PEGPH20+ gemcitabine).

Accuracy and Performance of Functional Parameter Estimation Using a Novel Numerical Optimization Approach for GPU-Based Kinetic Compartmental Modeling

Quantitative kinetic parameters derived from dynamic contrast-enhanced (DCE) data are dependent on signal measurement quality and choice of pharmacokinetic model. However, the fundamental optimization analysis method is equally important and its impact on pharmacokinetic parameters has been mostly overlooked. We examine the effects of those choices on accuracy and performance of parameter estimation using both computer processing unit and graphical processing unit (GPU) numerical optimization implementations and evaluate the improvements offered by a novel optimization approach. A test framework was developed where experimentally derived population-average arterial input function and randomly sampled parameter sets {K_trans, K_ep, V_b, τ} were used to generate known tissue curves. Five numerical optimization algorithms were evaluated: sequential quadratic programming, downhill simplex (Nelder–Mead), pattern search, simulated annealing, and differential evolution. This was combined with various objective function implementation details: delay approximation, discretization and varying sampling rates. Then, impact of noise and CPU/GPU implementation was tested for speed and accuracy. Finally, the optimal method was compared to conventional implementation as applied to clinical DCE computed tomography. Nelder–Mead, differential evolution and sequential quadratic programming produced good results on clean and noisy input data outperforming simulated annealing and pattern search in terms of speed and accuracy in the respective order of 10⁻⁸%, 10⁻⁷%, and ×10⁻⁶%). A novel approach for DCE numerical optimization (infinite impulse response with fractional delay approximation) was implemented on GPU for speed increase of at least 2 orders of magnitude. Applied to clinical data, the magnitude of overall parameter error was <10%.

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.

Modified MR dispersion imaging in prostate dynamic contrast-enhanced MRI

BACKGROUND

An estimation of an intravascular dispersion parameter was previously proposed to improve the overall accuracy and precision of the model parameters, but the high computation complexity can limit its practical usability in prostate dynamic contrast-enhanced MRI (DCE-MRI).

PURPOSE

To compare and evaluate the model fitting uncertainty and error in the model parameter estimation using different DCE-MRI analysis models and to evaluate the ability of the intravascular dispersion parameter to delineate between noncancerous and cancerous prostate tissue in the transition and peripheral zones.

STUDY TYPE

Retrospective.

POPULATION

Fifty-three patients who underwent radical prostatectomy.

FIELD STRENGTH/SEQUENCE

3 T/3D RF-spoiled gradient echo sequence.

ASSESSMENT

The coefficient of variation was used to assess the model fitting uncertainty by adding random noise to the time-concentration curves, and the Akaike information criterion was used to assess the model fitting error. The parametric maps derived from four DCE-MRI analysis models were evaluated by evaluating the delineation between noncancerous tissue and prostate cancer or clinically significant prostate cancer.

STATISTICAL TESTS

The receiver operating curve analysis was performed to compare the ability to delineate between noncancerous and prostate cancer tissue in the transition and peripheral zones.

RESULTS

Both MR dispersion imaging (MRDI) and Weinmann analysis models had the maximum coefficient of variation in different tissue types, while the model fitting uncertainty of modified (m)MRDI was similar to the standard Toft model. In mMRDI, the model fitting error was minimum, and the delineation between noncancerous and clinically significant prostate cancer tissue was improved in both transition (area under the curve [AUC] = 0.92) and peripheral zones (AUC = 0.92), in comparison with MRDI (AUC = 0.89 and AUC = 0.85, respectively).

DATA CONCLUSION

The mMRDI showed promising results in detecting prostate cancer while maintaining a similar model fitting uncertainty.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:1307-1317.

The Use of Quantitative Imaging in Radiation Oncology: A Quantitative Imaging Network (QIN) Perspective

Modern radiation therapy is delivered with great precision, in part by relying on high-resolution multidimensional anatomic imaging to define targets in space and time. The development of quantitative imaging (QI) modalities capable of monitoring biologic parameters could provide deeper insight into tumor biology and facilitate more personalized clinical decision-making. The Quantitative Imaging Network (QIN) was established by the National Cancer Institute to advance and validate these QI modalities in the context of oncology clinical trials. In particular, the QIN has significant interest in the application of QI to widen the therapeutic window of radiation therapy. QI modalities have great promise in radiation oncology and will help address significant clinical needs, including finer prognostication, more specific target delineation, reduction of normal tissue toxicity, identification of radioresistant disease, and clearer interpretation of treatment response. Patient-specific QI is being incorporated into radiation treatment design in ways such as dose escalation and adaptive replanning, with the intent of improving outcomes while lessening treatment morbidities. This review discusses the current vision of the QIN, current areas of investigation, and how the QIN hopes to enhance the integration of QI into the practice of radiation oncology.

Imaging vascular and hemodynamic features of the brain using dynamic susceptibility contrast and dynamic contrast enhanced MRI

In the context of neurologic disorders, dynamic susceptibility contrast (DSC) and dynamic contrast enhanced (DCE) MRI provide valuable insights into cerebral vascular function, integrity, and architecture. Even after two decades of use, these modalities continue to evolve as their biophysical and kinetic basis is better understood, with improvements in pulse sequences and accelerated imaging techniques and through application of more robust and automated data analysis strategies. Here, we systematically review each of these elements, with a focus on how their integration improves kinetic parameter accuracy and the development of new hemodynamic biomarkers that provide sub-voxel sensitivity (e.g., capillary transit time and flow heterogeneity). Regarding contrast mechanisms, we discuss the dipole-dipole interactions and susceptibility effects that give rise to simultaneous T₁, T₂ and T₂^∗ relaxation effects, including their quantification, influence on pulse sequence parameter optimization, and use in methods such as vessel size and vessel architectural imaging. The application of technologic advancements, such as parallel imaging, simultaneous multi-slice, undersampled k-space acquisitions, and sliding window strategies, enables improved spatial and/or temporal resolution of DSC and DCE acquisitions. Such acceleration techniques have also enabled the implementation of, clinically feasible, simultaneous multi-echo spin- and gradient echo acquisitions, providing more comprehensive and quantitative interrogation of T₁, T₂ and T₂^∗ changes. Characterizing these relaxation rate changes through different post-processing options allows for the quantification of hemodynamics and vascular permeability. The application of different biophysical models provides insight into traditional hemodynamic parameters (e.g., cerebral blood volume) and more advanced parameters (e.g., capillary transit time heterogeneity). We provide insight into the appropriate selection of biophysical models and the necessary post-processing steps to ensure reliable measurements while minimizing potential sources of error. We show representative examples of advanced DSC- and DCE-MRI methods applied to pathologic conditions affecting the cerebral microcirculation, including brain tumors, stroke, aging, and multiple sclerosis. The maturation and standardization of conventional DSC- and DCE-MRI techniques has enabled their increased integration into clinical practice and use in clinical trials, which has, in turn, spurred renewed interest in their technological and biophysical development, paving the way towards a more comprehensive assessment of cerebral hemodynamics.