Review of treatment assessment using DCE-MRI in breast cancer radiation therapy

Meta-Analysis of MRI in Predicting Early Response to Radiotherapy and Chemotherapy in Esophageal Cancer

RATIONALE AND OBJECTIVES

At present, the application of magnetic resonance imaging (MRI) in the prediction of response to neoadjuvant therapy and concurrent chemoradiotherapy for the treatment of esophageal cancer still needs to be further explored, and its early differential value remains controversial, thus we carried out this systematic review with a meta-analysis. In the application, different MRI sequences and corresponding parameters are used for the differential diagnosis of the response to neoadjuvant therapy and concurrent chemoradiotherapy.

METHODS

All relevant studies evaluated the efficacy and response to MRI in neoadjuvant therapy or concurrent chemoradiotherapy for esophageal cancer on Pubmed, Embase, Cohrane Library, and Web of Science databases published before October 10, 2023 (inclusive) were systematically searched. A revised tool was used to assess the quality of diagnostic accuracy studies (QUADAS-2) to assess the risk of bias in the included original studies. A subgroup analysis of MRI sequences diffusion weighted imaging (DWI), dynamic contrast enhanced (DCE) and their corresponding different parameters, as well as the acquisition timepoints (before and after treatment) for different parameters, was performed during the meta-analysis. The bivariate mixed-effects model was used for meta-analysis.

RESULTS

21 studies were finally included, involving 1128 patients with esophageal cancer. The sensitivity, specificity, and area under receiver operating characteristic curve (ROC curve) of DWI sequence for identifying response to concurrent chemoradiotherapy were 0.82 (95% CI: 0.74-0.87), 0.81 (95% CI: 0.72-0.87) and 0.88 (95% CI: 0.56-0.98), respectively. The sensitivity, specificity, and area under ROC curve of DCE sequence for identifying response to concurrent chemoradiotherapy were 0.78 (95% CI: 0.70-0.84), 0.65 (95% CI: 0.59-0.70) and 0.73 (95% CI: 0.50-0.88), respectively. In patients with esophageal cancer, the sensitivity, specificity, and area under the ROC curve of DWI sequences for identifying response to neoadjuvant therapy were 0.80 (95% CI: 0.69 - 0.88), 0.81 (95% CI: 0.69 - 0.89), and 0.88 (95% CI: 0.34 - 0.99), respectively; the sensitivity, specificity, and area under the ROC curve of DCE sequences for identifying response to neoadjuvant therapy were 0.84 (95% CI: 0.76 - 0.90), 0.61 (95% CI: 0.53 - 0.68), and 0.70 (95% CI: 0.27 - 0.94), respectively.

CONCLUSIONS

Based on the available evidence, MRI had a very good value in the early identification of response to neoadjuvant therapy and concurrent chemoradiotherapy for esophageal cancer, especially DWI. Apparent diffusion coefficient (ADC) value changes before and after treatment could be used as predictors of pathological response. Also, ADC value changes before and after treatment could be used as a tool to guide clinical decision-making.

Pharmacokinetic Analysis of Enhancement-Constrained Acceleration (ECA) reconstruction-based high temporal resolution breast DCE-MRI

The high spatial and temporal resolution of dynamic contrast-enhanced MRI (DCE-MRI) can improve the diagnostic accuracy of breast cancer screening in patients who have dense breasts or are at high risk of breast cancer. However, the spatiotemporal resolution of DCE-MRI is limited by technical issues in clinical practice. Our earlier work demonstrated the use of image reconstruction with enhancement-constrained acceleration (ECA) to increase temporal resolution. ECA exploits the correlation in k-space between successive image acquisitions. Because of this correlation, and due to the very sparse enhancement at early times after contrast media injection, we can reconstruct images from highly under-sampled k-space data. Our previous results showed that ECA reconstruction at 0.25 seconds per image (4 Hz) can estimate bolus arrival time (BAT) and initial enhancement slope (iSlope) more accurately than a standard inverse fast Fourier transform (IFFT) when k-space data is sampled following a Cartesian based sampling trajectory with adequate signal-to-noise ratio (SNR). In this follow-up study, we investigated the effect of different Cartesian based sampling trajectories, SNRs and acceleration rates on the performance of ECA reconstruction in estimating contrast media kinetics in lesions (BAT, iSlope and Ktrans) and in arteries (Peak signal intensity of first pass, time to peak, and BAT). We further validated ECA reconstruction with a flow phantom experiment. Our results show that ECA reconstruction of k-space data acquired with 'Under-sampling with Repeated Advancing Phase' (UnWRAP) trajectories with an acceleration factor of 14, and temporal resolution of 0.5 s/image and high SNR (SNR ≥ 30 dB, noise standard deviation (std) < 3%) ensures minor errors (5% or 1 s error) in lesion kinetics. Medium SNR (SNR ≥ 20 dB, noise std ≤ 10%) was needed to accurately measure arterial enhancement kinetics. Our results also suggest that accelerated temporal resolution with ECA with 0.5 s/image is practical.

Comparison of DCE-MRI parametric mapping using MP2RAGE and variable flip angle T1 mapping

Quantitative dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) measures the rate of transfer of contrast agent from the vascular space to the tissue space by fitting signal-time data to pharmacokinetic models. However, these models are very sensitive to errors in T₁ mapping. Accurate T₁ mapping is necessary for high quality quantitative DCE-MRI studies. This study compares magnetization prepared rapid (two) gradient echo sequence (MP2RAGE) T₁-mapping accuracy to the conventional variable flip angle (VFA) approach, and also determines the effect of the new T₁-mapping method on the K^trans parameter. VFA and MP2RAGE T₁ values were compared to the gold standard inverse recovery (IR) method in phantom over manually drawn ROIs. In vivo, ROIs were manually drawn over prostate and prostatic lesions. Average T₁ values over ROIs were compared and K^trans maps for each method were calculated via the extended Tofts model. VFA-T₁ maps overestimated T₁ values by up to 50% compared to gold standard IR T₁ values in phantom. MP2RAGE differed by up to 9%. MP2RAGE-T₁ and K^trans values were significantly different from VFA values over prostatic lesions (p < 0.05). K^trans was consistently underestimated using VFA compared to MP2RAGE (p < 0.05). MP2RAGE T₁ maps are shown to be more accurate, leading to more reliable pharmacokinetic modeling. This can potentially lead to better lesion characterization and improve clinical outcomes.

The Influence of Data-Driven Compressed Sensing Reconstruction on Quantitative Pharmacokinetic Analysis in Breast DCE MRI

Radial acquisition with MOCCO reconstruction has been previously proposed for high spatial and temporal resolution breast DCE imaging. In this work, we characterize MOCCO across a wide range of temporal contrast enhancement in a digital reference object (DRO). Time-resolved radial data was simulated using a DRO with lesions in different PK parameters. The under sampled data were reconstructed at 5 s temporal resolution using the data-driven low-rank temporal model for MOCCO, compressed sensing with temporal total variation (CS-TV) and more conventional low-rank reconstruction (PCB). Our results demonstrated that MOCCO was able to recover curves with K^trans values ranging from 0.01 to 0.8 min⁻¹ and fixed V_e = 0.3, where the fitted results are within a 10% bias error range. MOCCO reconstruction showed less impact on the selection of different temporal models than conventional low-rank reconstruction and the greater error was observed with PCB. CS-TV showed overall underestimation in both K^trans and V_e. For the Monte-Carlo simulations, MOCCO was found to provide the most accurate reconstruction results for curves with intermediate lesion kinetics in the presence of noise. Initial in vivo experiences are reported in one patient volunteer. Overall, MOCCO was able to provide reconstructed time-series data that resulted in a more accurate measurement of PK parameters than PCB and CS-TV.

An Anthropomorphic Digital Reference Object (DRO) for Simulation and Analysis of Breast DCE MRI Techniques

Advances in accelerated magnetic resonance imaging (MRI) continue to push the bounds on achievable spatial and temporal resolution while maintaining a clinically acceptable image quality. Validation tools, including numerical simulations, are needed to characterize the repeatability and reproducibility of such methods for use in quantitative imaging applications. We describe the development of a simulation framework for analyzing and optimizing accelerated MRI acquisition and reconstruction techniques used in dynamic contrast enhanced (DCE) breast imaging. The simulation framework, in the form of a digital reference object (DRO), consists of four modules that control different aspects of the simulation, including the appearance and physiological behavior of the breast tissue as well as the MRI acquisition settings, to produce simulated k-space data for a DCE breast exam. The DRO design and functionality are described along with simulation examples provided to show potential applications of the DRO. The included simulation results demonstrate the ability of the DRO to simulate a variety of effects including the creation of simulated lesions, tissue enhancement modeled by the generalized kinetic model, T1-relaxation, fat signal precession and saturation, acquisition SNR, and changes in temporal resolution.

Dynamic Contrast-enhanced and Diffusion-weighted Magnetic Resonance Imaging for Response Evaluation After Single-Dose Ablative Neoadjuvant Partial Breast Irradiation

Purpose

We aimed to evaluate changes in dynamic contrast-enhanced (DCE) and diffusion-weighted (DW) magnetic resonance imaging (MRI) scans acquired before and after single-dose ablative neoadjuvant partial breast irradiation (NA-PBI), and explore the relation between semiquantitative MRI parameters and radiologic and pathologic responses.

Methods and Materials

We analyzed 3.0T DCE and DW-MRI of 36 patients with low-risk breast cancer who were treated with single-dose NA-PBI, followed by breast-conserving surgery 6 or 8 months later. MRI was acquired before NA-PBI and 1 week, 2, 4, and 6 months after NA-PBI. Breast radiologists assessed the radiologic response and breast pathologists scored the pathologic response after surgery. Patients were grouped as either pathologic responders or nonresponders (<10% vs ≥10% residual tumor cells). The semiquantitative MRI parameters evaluated were time to enhancement (TTE), 1-minute relative enhancement (RE_1min), percentage of enhancing voxels (%EV), distribution of washout curve types, and apparent diffusion coefficient (ADC).

Results

In general, the enhancement increased 1 week after NA-PBI (baseline vs 1 week median – TTE: 15s vs 10s; RE_1min: 161% vs 197%; %EV: 47% vs 67%) and decreased from 2 months onward (6 months median – TTE: 25s; RE_1min: 86%; %EV: 12%). Median ADC increased from 0.83 × 10⁻³ mm²/s at baseline to 1.28 × 10⁻³ mm²/s at 6 months. TTE, RE_1min, and %EV showed the most potential to differentiate between radiologic responses, and TTE, RE_1min, and ADC between pathologic responses.

Conclusions

Semiquantitative analyses of DCE and DW-MRI showed changes in relative enhancement and ADC 1 week after NA-PBI, indicating acute inflammation, followed by changes indicating tumor regression from 2 to 6 months after radiation therapy. A relation between the MRI parameters and radiologic and pathologic responses could not be proven in this exploratory study.

Vascularity and Dynamic Contrast-Enhanced Breast Magnetic Resonance Imaging

Angiogenesis is a key step in the initiation and progression of an invasive breast cancer. High microvessel density by morphological characterization predicts metastasis and poor survival in women with invasive breast cancers. However, morphologic characterization is subject to variability and only can evaluate a limited portion of an invasive breast cancer. Consequently, breast Magnetic Resonance Imaging (MRI) is currently being evaluated to assess vascularity. Recently, through the new field of radiomics, dynamic contrast enhanced (DCE)-MRI is being used to evaluate vascular density, vascular morphology, and detection of aggressive breast cancer biology. While DCE-MRI is a highly sensitive tool, there are specific features that limit computational evaluation of blood vessels. These include (1) DCE-MRI evaluates gadolinium contrast and does not directly evaluate biology, (2) the resolution of DCE-MRI is insufficient for imaging small blood vessels, and (3) DCE-MRI images are very difficult to co-register. Here we review computational approaches for detection and analysis of blood vessels in DCE-MRI images and present some of the strategies we have developed for co-registry of DCE-MRI images and early detection of vascularization.

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.

Texture Analysis of DCE-MRI Intratumoral Subregions to Identify Benign and Malignant Breast Tumors

Purpose

To evaluate the potential of the texture features extracted from dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) intratumoral subregions to distinguish benign from malignant breast tumors.

Materials and Methods

A total of 299 patients with pathologically verified breast tumors who underwent breast DCE-MRI examination were enrolled in this study, including 124 benign cases and 175 malignant cases. The whole tumor area was semi-automatically segmented on the basis of subtraction images of DCE-MRI in Matlab 2018b. According to the time to peak of the contrast agent, the whole tumor area was partitioned into three subregions: early, moderate, and late. A total of 467 texture features were extracted from the whole tumor area and the three subregions, respectively. Patients were divided into training (n = 209) and validation (n = 90) cohorts by different MRI scanners. The least absolute shrinkage and selection operator (LASSO) method was used to select the optimal feature subset in the training cohort. The Kolmogorov-Smirnov test was first performed on texture features selected by LASSO to test whether the samples followed a normal distribution. Two machine learning methods, decision tree (DT) and support vector machine (SVM), were used to establish classification models with a 10-fold cross-validation method. The performance of the classification models was evaluated with receiver operating characteristic (ROC) curves.

Results

In the training cohort, the areas under the ROC curve (AUCs) for the DT_Whole model and SVM_Whole model were 0.744 and 0.806, respectively. In contrast, the AUCs of the DT_Early model (P = 0.004), DT_Late model (P = 0.015), SVM_Early model (P = 0.002), and SVM_Late model (P = 0.002) were significantly higher: 0.863 (95% CI, 0.808-0.906), 0.860 (95% CI, 0.806-0.904), 0.934 (95% CI, 0.891-0.963), and 0.921 (95% CI, 0.876-0.954), respectively. The SVM_Early model and SVM_Late model achieved better performance than the DT_Early model and DT_Late model (P = 0.003, 0.034, 0.008, and 0.026, respectively). In the validation cohort, the AUCs for the DT_Whole model and SVM_Whole model were 0.670 and 0.708, respectively. In comparison, the AUCs of the DT_Early model (P = 0.006), DT_Late model (P = 0.043), SVM_Early model (P = 0.001), and SVM_Late model (P = 0.007) were significantly higher: 0.839 (95% CI, 0.747-0.908), 0.784 (95% CI, 0.601-0.798), 0.890 (95% CI, 0.806-0.946), and 0.865 (95% CI, 0.777-0.928), respectively.

Conclusion

The texture features from intratumoral subregions of breast DCE-MRI showed potential in identifying benign and malignant breast tumors.

Five-dimensional quantitative low-dose Multitasking dynamic contrast- enhanced MRI: Preliminary study on breast cancer

PURPOSE

To develop a low-dose Multitasking DCE technique (LD-MT-DCE) for breast imaging, enabling dynamic T₁ mapping-based quantitative characterization of tumor blood flow and vascular properties with whole-breast coverage, a spatial resolution of 0.9 × 0.9 × 1.1 mm³ , and a temporal resolution of 1.4 seconds using a 20% gadolinium dose (0.02 mmol/kg).

METHODS

Magnetic resonance Multitasking was used to reconstruct 5D images with three spatial dimensions, one T₁ recovery dimension for dynamic T₁ quantification, and one DCE dimension for contrast kinetics. Kinetic parameters F p , v p , K trans , and v e were estimated from dynamic T₁ maps using the two-compartment exchange model. The LD-MT-DCE repeatability and agreement against standard-dose MT-DCE were evaluated in 20 healthy subjects. In 7 patients with triple-negative breast cancer, LD-MT-DCE image quality and diagnostic results were compared with that of standard-dose clinical DCE in the same imaging session. One-way unbalanced analysis of variance with Tukey test was performed to evaluate the statistical significance of the kinetic parameters between control and patient groups.

RESULTS

The LD-MT-DCE technique was repeatable, agreed with standard-dose MT-DCE, and showed excellent image quality. The diagnosis using LD-MT-DCE matched well with clinical results. The values of F p , v p , and K trans were significantly different between malignant tumors and normal breast tissue (P < .001, < .001, and < .001, respectively), and between malignant and benign tumors (P = .020, .003, and < .001, respectively).

CONCLUSION

The LD-MT-DCE technique was repeatable and showed excellent image quality and equivalent diagnosis compared with standard-dose clinical DCE. The estimated kinetic parameters were capable of differentiating between normal breast tissue and benign and malignant tumors.

The Effect of Registration on Voxel-Wise Tofts Model Parameters and Uncertainties from DCE-MRI of Early-Stage Breast Cancer Patients Using 3DSlicer

We quantitatively investigate the influence of image registration, using open-source software (3DSlicer), on kinetic analysis (Tofts model) of dynamic contrast enhanced MRI of early-stage breast cancer patients. We also show that registration computation time can be reduced by reducing the percent sampling (PS) of voxels used for estimation of the cost function. DCE-MRI breast images were acquired on a 3T-PET/MRI system in 13 patients with early-stage breast cancer who were scanned in a prone radiotherapy position. Images were registered using a BSpline transformation with a 2 cm isotropic grid at 100, 20, 5, 1, and 0.5PS (BRAINSFit in 3DSlicer). Signal enhancement curves were analyzed voxel-by-voxel using the Tofts kinetic model. Comparing unregistered with registered groups, we found a significant change in the 90th percentile of the voxel-wise distribution of K^trans. We also found a significant reduction in the following: (1) in the standard error (uncertainty) of the parameter value estimation, (2) the number of voxel fits providing unphysical values for the extracellular-extravascular volume fraction (v_e > 1), and (3) goodness of fit. We found no significant differences in the median of parameter value distributions (K^trans, v_e) between unregistered and registered images. Differences between parameters and uncertainties obtained using 100PS versus 20PS were small and statistically insignificant. As such, computation time can be reduced by a factor of 2, on average, by using 20PS while not affecting the kinetic fit. The methods outlined here are important for studies including a large number of post-contrast images or number of patient images.

Comparison of postoperative CT- and preoperative MRI-based breast tumor bed contours in prone position for radiotherapy after breast-conserving surgery

Objectives

To compare the target volume of tumor bed defined by postoperative computed tomography (post-CT) in prone position registered with or without preoperative magnetic resonance imaging (pre-MRI).

Methods

A total of 22 patients were included with early-stage breast invasive ductal cancer, who have undergone breast-conservative surgery and received the pre-MRI and post-CT in prone position. The MRI sequences (T1W, T2W, T2W-SPAIR, DWI, dyn-eTHRIVE, sdyn-eTHRIVE) were delineated and manually registered to CT, respectively. The clinical target volumes (CTVs) and planning target volumes (PTVs) were contoured on CT and different MRI sequences, respectively. Differences were measured in terms of consistence index (CI), dice coefficient (DC), geographical miss index (GMI), and normal tissue index (NTI).

Results

The differences of delineation volumes among CT and MRIs were significant, both in the CTVs (p = 0.035) and PTVs (p < 0.001). The values of CI and DC for sdyn-eTHRIVE registration to CT were the largest among all MRI sequences, but GMI and NTI were the smallest. No obvious linear correlation (p > 0.05) between the CI derived from the registration of CT and sdyn-eTHRIVE of CTV with the breast volume, the cavity visualization score (CVS) of CT, time interval from surgery to CT simulation, the maximum diameter of the intraoperative mass, and the number of titanium clips, respectively.

Conclusions

The CTVs and PTVs in MRI sequences were all smaller than those in CT. The pre-MRI, especially the sdyn-eTHRIVE, could be used to optimize the post-CT-based target delineation of breast cancer.

Key Points

• Registered pre-MRI to post-CT in order to improve the accuracy of target volume delineation of breast cancer.

• The CTVs and PTVs in MRI sequences were all smaller than those in CT.

• The sdyn-eTHRIVE of pre-MRIs may be a better choice to improve the delineation of CT-based CTV and PTV.

Electronic supplementary material

The online version of this article (10.1007/s00330-020-07085-0) contains supplementary material, which is available to authorized users.

Quantitative ultrasound radiomics for therapy response monitoring in patients with locally advanced breast cancer: Multi-institutional study results

BACKGROUND

Neoadjuvant chemotherapy (NAC) is the standard of care for patients with locally advanced breast cancer (LABC). The study was conducted to investigate the utility of quantitative ultrasound (QUS) carried out during NAC to predict the final tumour response in a multi-institutional setting.

METHODS

Fifty-nine patients with LABC were enrolled from three institutions in North America (Sunnybrook Health Sciences Centre (Toronto, Canada), MD Anderson Cancer Centre (Texas, USA), and Princess Margaret Cancer Centre (Toronto, Canada)). QUS data were collected before starting NAC and subsequently at weeks 1 and 4 during chemotherapy. Spectral tumour parametric maps were generated, and textural features determined using grey-level co-occurrence matrices. Patients were divided into two groups based on their pathological outcomes following surgery: responders and non-responders. Machine learning algorithms using Fisher's linear discriminant (FLD), K-nearest neighbour (K-NN), and support vector machine (SVM-RBF) were used to generate response classification models.

RESULTS

Thirty-six patients were classified as responders and twenty-three as non-responders. Among all the models, SVM-RBF had the highest accuracy of 81% at both weeks 1 and week 4 with area under curve (AUC) values of 0.87 each. The inclusion of week 1 and 4 features led to an improvement of the classifier models, with the accuracy and AUC from baseline features only being 76% and 0.68, respectively.

CONCLUSION

QUS data obtained during NAC reflect the ongoing treatment-related changes during chemotherapy and can lead to better classifier performances in predicting the ultimate pathologic response to treatment compared to baseline features alone.

Combination of DCE-MRI and DWI in Predicting the Treatment Effect of Concurrent Chemoradiotherapy in Esophageal Carcinoma

Background

Concurrent chemoradiotherapy (CCRT) is the main treatment for esophageal cancer, but the response to treatment varies from individual to individual. MR imaging methods, such as diffusion-weighted (DW) MRI and the use of dynamic contrast-enhanced (DCE) MRI, have the potential to provide additional biomarkers that could evaluate the effect of CCRT in patients with esophageal carcinoma.

Materials and Methods

Fifty-six patients with esophageal carcinoma, verified by histopathology, underwent MRI examination before and at midtreatment (4^th week, radiotherapy 30-40 Gy) using the Siemens 3.0 T MR System. Parameter maps of apparent diffusion coefficient (ADC), and DCE maps of volume transfer constant (K ^rans), rate contrast (k _ep), and extracellular fluid space (v _e), were computed using a Siemens Company Multimodality Workplace (MMWP) model. Comparison of histogram parameters and their diagnostic performance was determined using the Mann-Whitney U test and receiver operating characteristic (ROC) analysis.

Results

56 patient MRI scans were available for analysis at baseline and at the third week, respectively. Pretreatment K ^rans, pretreatment k _ep, pretreatment ADC (P < 0.05), and during-treatment K ^rans (P < 0.05) and ΔK ^rans and ΔADC (P < 0.05) were significantly different after CCRT. Based on the binary logistic model, the ROC analysis demonstrated that the combined predictors demonstrated a high diagnostic performance with an AUC of 0.939. The sensitivity and specificity were 98.6% and 73.8%, respectively.

Conclusion

The combination of DCE and DWI can be used as an early biomarker in the prediction of the effect of CCRT three weeks after treatment in esophageal carcinoma.

Evaluating the Relationship Between Dynamic Contrast-Enhanced MRI (DCE-MRI) Parameters and Pathological Characteristics in Breast Cancer

BACKGROUND

Dynamic contrast-enhanced MRI (DCE-MRI) is used to evaluate tumor microvasculature. However, studies demonstrating an association between perfusion parameters derived from DCE-MRI and histopathologic characteristics are limited to a small set of histopathologic factors, and the results are inconsistent.

PURPOSE

To evaluate the relationship between DCE-MRI perfusion parameters and common histopathologic tumor characteristics used to predict angiogenesis and determine prognosis in breast cancer.

STUDY TYPE

Retrospective.

POPULATION

In all, 105 breast cancer patients with invasive ductal carcinoma (122 lesions).

FIELD STRENGTH/SEQUENCE

3.0T, turbo spin-echo (TSE) T₁ -weighted, fat-suppressed T₂ -weighted, TSE T₂ -weighted, and dynamic unenhanced and contrast-enhanced 3D T₁ high-resolution isotropic volume examination.

ASSESSMENT

One reviewer obtained perfusion parameters (K^trans , k_ep , v_e , and v_p ) of each breast cancer from DCE MRI using the extended Tofts model with a fixed baseline T₁ value and a population-based arterial input function. The relationship between DCE-MRI perfusion parameters and histopathologic tumor characteristics used to predict angiogenesis and determine prognosis was evaluated.

STATISTICAL TESTS

Student's t-test, Mann-Whitney U-test, analysis of variance (ANOVA), and Kruskal-Wallis test were used.

RESULTS

Triple-negative breast cancers exhibited higher K^trans and k_ep than luminal cancers (P < 0.05). Estrogen receptor (ER)-negative tumors showed higher K^trans than ER-positive tumors (P < 0.05). Progesterone receptor (PR)-negative tumors presented higher v_e than PR-positive tumors (P < 0.05). Tumors with higher Ki-67 showed higher k_ep than tumors with lower Ki-67 (P < 0.05). P53-positive tumors exhibited higher K^trans and k_ep than p53-negative tumors (P < 0.05). Higher histologic grade tumors (grade II/III) presented higher K^trans , k_ep , v_p (P < 0.05) than grade I tumors. Tumors with LVSI presented higher K^trans and k_ep than tumors without LVSI (P < 0.05).

DATA CONCLUSION

Breast cancer presenting higher K^trans and k_ep on DCE-MRI was associated with poor prognostic histopathologic factors. Therefore, pretreatment DCE-MRI perfusion parameters may be useful imaging biomarkers for the evaluation of tumor prognosis and angiogenesis.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY STAGE: 2.

DCE-MRI assessment of response to neoadjuvant SABR in early stage breast cancer: Comparisons of single versus three fraction schemes and two different imaging time delays post-SABR

PURPOSE

To determine the effect of dose fractionation and time delay post-neoadjuvant stereotactic ablative radiotherapy (SABR) on dynamic contrast-enhanced (DCE)-MRI parameters in early stage breast cancer patients.

MATERIALS AND METHODS

DCE-MRI was acquired in 17 patients pre- and post-SABR. Five patients were imaged 6-7 days post-21 Gy/1fraction (group 1), six 16-19 days post-21 Gy/1fraction (group 2), and six 16-18 days post-30 Gy/3 fractions every other day (group 3). DCE-MRI scans were performed using half the clinical dose of contrast agent. Changes in the surrounding tissue were quantified using a signal-enhancement threshold metric that characterizes changes in signal-enhancement volume (SEV). Tumour response was quantified using Ktrans and ve (Tofts model) pre- and post-SABR. Significance was assessed using a Wilcoxin signed-rank test.

RESULTS

All group 1 and 4/6 group 2 patients' SEV increased post-SABR. All group 3 patients' SEV decreased. The mean Ktrans increased for group 1 by 76% (p = 0.043) while group 2 and 3 decreased 15% (p = 0.028) and 34% (p = 0.028), respectively. For ve, there was no significant change in Group 1 (p = 0.35). Groups 2 showed an increase of 24% (p = 0.043), and Group 3 trended toward an increase (23%, p = 0.08).

CONCLUSION

Kinetic parameters measured 2.5 weeks post-SABR in both single fraction and three fraction groups were indicative of response but only the single fraction protocol led to enhancement in the surrounding tissue. Our results also suggest that DCE-MRI one-week post-SABR may be too early for response assessment, at least for single fraction SABR, whereas 2.5 weeks appears sufficiently long to minimize confounding acute effects.

Impact of the Number of Iterations in Compressed Sensing Reconstruction on Ultrafast Dynamic Contrast-enhanced Breast MR Imaging

PURPOSE

To assess the impact of the number of iterations of compressed sensing (CS) reconstruction on the kinetic parameters and image quality in dynamic contrast-enhanced (DCE)-MRI of the breast, with prospectively undersampled CS-accelerated scans.

MATERIALS AND METHODS

Breast examinations including ultrafast DCE-MRI using CS were conducted for 21 patients. Images were reconstructed with different numbers of iterations. The peak enhancement ratio of the aorta and wash-in slope, initial area under the curve, and Ktrans of the breast lesions were measured. The root mean square error and structural similarity between the images using 50 iterations and images with a lower number of iterations were evaluated as criterion for quantitative image evaluation.

RESULTS

Using an insufficient number of iterations, the contrast-enhanced effect was highly underestimated. In all semi-quantitative parameters, the number of iterations that stabilized the parameters in malignant lesions was higher than that in benign lesions. At least 15 iterations were needed for semi-quantitative parameters. For Ktrans, there were no significant differences between 10 and 50 iterations in both malignant and benign lesions.

CONCLUSION

The kinetic parameters using ultrafast DCE-MRI with CS are affected by the number of iterations, especially in malignant lesions. However, if the images are reconstructed with an adequate number of iterations, ultrafast DCE-MRI with CS can be a powerful technique having high temporal and spatial resolution.

Differentiating the lymph node metastasis of breast cancer through dynamic contrast-enhanced magnetic resonance imaging

Objective:

Lymph node metastasis is an important trait of breast cancer, and tumors with different lymph node statuses require various clinical treatments. This study was designed to evaluate the lymph node metastasis of breast cancer through pharmacokinetic and histogram analysis via dynamic contrast-enhanced (DCE) MRI.

Methods and materials:

A retrospective analysis was conducted to quantitatively evaluate the lymph node statuses of patients with breast cancer. A total of 75 patients, i.e. 34 patients with lymph node metastasis and 41 patients without lymph node metastasis, were involved in this research. Of the patients with lymph node metastases, 19 had sentinel lymph node metastasis, and 15 had axillary lymph node metastasis. MRI was conducted using a 3.0 T imaging device. Segmentation was carried out on the regions of interest (ROIs) in breast tumors under DCE-MRI, and pharmacokinetic and histogram parameters were calculated from the same ROIs. Mann–Whitney U test was performed, and receiver operating characteristic curves for the parameters of the two groups were constructed to determine their diagnostic values.

Results:

Pharmacokinetic parameters, including Ktrans, Kep, area under the curve of time–concentration, and time to peak, which were derived from the extended Tofts linear model for DCE-MRI, could highlight the tumor areas in the breast and reveal the increased perfusion. Conversely, the pharmacokinetic parameters showed no significant difference between the patients with and without lymph node metastases. By contrast, the parameters from the histogram analysis yielded promising results. The entropy of the ROIs exhibited the best diagnostic ability between patients with and without lymph node metastases (p < 0.01, area under the curve of receiver operating characteristic = 0.765, specificity = 0.706, sensitivity = 0.780).

Conclusion:

In comparison with the pharmacokinetic parameters, the histogram analysis of the MR images could reveal the differences between patients with and without lymph node metastases. The entropy from the histogram indicated that the diagnostic ability was highly sensitive and specific.

Advances in knowledge:

This research gave out a promising result on the differentiating lymph node metastases through histogram analysis on tumors in DCE-MR images. Histogram could reveal the tumors heterogenicity between patients with different lymph node status.

DCE-MRI-Derived Volume Transfer Constant (Ktrans) and DWI Apparent Diffusion Coefficient as Predictive Markers of Short- and Long-Term Efficacy of Chemoradiotherapy in Patients With Esophageal Cancer

This study aimed to evaluate both the short- and long-term efficacies of chemoradiotherapy in relation to the treatment of esophageal cancer . This was achieved through the use of dynamic contrast-enhanced magnetic resonance imaging–derived volume transfer constant and diffusion weighted imaging–derived apparent diffusion coefficient . Patients with esophageal cancer were assigned into the sensitive and resistant groups based on respective efficacies in chemoradiotherapy. Dynamic contrast-enhanced magnetic resonance imaging and diffusion weighted imaging were used to measure volume transfer constant and apparent diffusion coefficient, while computed tomography was used to calculate tumor size reduction rate. Pearson correlation analyses were conducted to analyze correlation between volume transfer constant, apparent diffusion coefficient, and the tumor size reduction rate. Receiver operating characteristic curve was constructed to analyze the short-term efficacy of volume transfer constant and apparent diffusion coefficient, while Kaplan-Meier curve was employed for survival rate analysis. Cox proportional hazard model was used for the risk factors for prognosis of patients with esophageal cancer. Our results indicated reduced levels of volume transfer constant, while increased levels were observed in ADC_min, ADC_mean, and ADC_max following chemoradiotherapy. A negative correlation was determined between ADC_min, ADC_mean, and ADC_max, as well as in the tumor size reduction rate prior to chemoradiotherapy, whereas a positive correlation was uncovered postchemoradiotherapy. Volume transfer constant was positively correlated with tumor size reduction rate both before and after chemoradiotherapy. The 5-year survival rate of patients with esophageal cancer having high ADC_min, ADC_mean, and ADC_max and volume transfer constant before chemoradiotherapy was greater than those with respectively lower values. According to the Cox proportional hazard model, ADC_mean, clinical stage, degree of differentiation, and tumor stage were all confirmed as being independent risk factors in regard to the prognosis of patients with EC. The findings of this study provide evidence suggesting that volume transfer constant and apparent diffusion coefficient as being tools allowing for the evaluation of both the short- and long-term efficacies of chemoradiotherapy esophageal cancer treatment.

Dynamic contrast-enhanced MRI for advanced esophageal cancer response assessment after concurrent chemoradiotherapy

PURPOSE

We aimed to evaluate the treatment response of patients with esophageal cancer after concurrent chemoradiation therapy (CRT) using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

METHODS

This retrospective study included 59 patients with histologically confirmed esophageal squamous cell carcinoma. The patients underwent DCE-MRI before and 4 weeks after CRT. Patients with complete response were defined as the CR group; partial response, stable disease, and progressive disease patients were defined as the non-CR group. DCE-MRI parameters (Ktrans, Ve, and Kep) were measured and compared between pre- and post-CRT in the CR and non-CR groups, respectively. Pre-CRT and post-CRT parameters were used to calculate the absolute change and the ratio of change. DCE-MRI parameters were compared between the CR and non-CR groups. Receiver operating characteristic (ROC) curves were used to verify diagnostic performance.

RESULTS

Patients with higher T-stage esophageal cancer might present with poorer response. After CRT, the Ktrans and Kep values significantly decreased in the CR group, whereas only Kep value decreased in the non-CR group. The post-Ktrans and post-Kep values were observed to be significantly lower in the CR group than in the non-CR group. The absolute change and ratio of change of both Ktrans and Kep were higher in the CR group than in the non-CR group. Based on ROC analysis, the ratio of change in Ktrans was the best parameter to assess treatment response (AUC= 0.840).

CONCLUSION

DCE-MRI parameters are valuable in predicting and assessing concurrent CRT response for advanced esophageal cancer.

Role and future of MRI in radiation oncology

Technical innovations and developments in areas such as disease localization, dose calculation algorithms, motion management and dose delivery technologies have revolutionized radiation therapy resulting in improved patient care with superior outcomes. A consequence of the ability to design and accurately deliver complex radiation fields is the need for improved target visualization through imaging. While CT imaging has been the standard of care for more than three decades, the superior soft tissue contrast afforded by MR has resulted in the adoption of this technology in radiation therapy. With the development of real time MR imaging techniques, the problem of real time motion management is enticing. Currently, the integration of an MR imaging and megavoltage radiation therapy treatment delivery system (MR-linac or MRL) is a reality that has the potential to provide improved target localization and real time motion management during treatment. Higher magnetic field strengths provide improved image quality potentially providing the backbone for future work related to image texture analysis-a field known as Radiomics-thereby providing meaningful information on the selection of future patients for radiation dose escalation, motion-managed treatment techniques and ultimately better patient care. On-going advances in MRL technologies promise improved real time soft tissue visualization, treatment margin reductions, beam optimization, inhomogeneity corrected dose calculation, fast multileaf collimators and volumetric arc radiation therapy. This review article provides rationale, advantages and disadvantages as well as ideas for future research in MRI related to radiation therapy mainly in adoption of MRL.

Hierarchical segmentation using equivalence test (HiSET): Application to DCE image sequences

Dynamical contrast enhanced (DCE) imaging allows non invasive access to tissue micro-vascularization. It appears as a promising tool to build imaging biomarkers for diagnostic, prognosis or anti-angiogenesis treatment monitoring of cancer. However, quantitative analysis of DCE image sequences suffers from low signal to noise ratio (SNR). SNR may be improved by averaging functional information in a large region of interest when it is functionally homogeneous. We propose a novel method for automatic segmentation of DCE image sequences into functionally homogeneous regions, called DCE-HiSET. Using an observation model which depends on one parameter a and is justified a posteriori, DCE-HiSET is a hierarchical clustering algorithm. It uses the p-value of a multiple equivalence test as dissimilarity measure and consists of two steps. The first exploits the spatial neighborhood structure to reduce complexity and takes advantage of the regularity of anatomical features, while the second recovers (spatially) disconnected homogeneous structures at a larger (global) scale. Given a minimal expected homogeneity discrepancy for the multiple equivalence test, both steps stop automatically by controlling the Type I error. This provides an adaptive choice for the number of clusters. Assuming that the DCE image sequence is functionally piecewise constant with signals on each piece sufficiently separated, we prove that DCE-HiSET will retrieve the exact partition with high probability as soon as the number of images in the sequence is large enough. The minimal expected homogeneity discrepancy appears as the tuning parameter controlling the size of the segmentation. DCE-HiSET has been implemented in C++ for 2D and 3D image sequences with competitive speed.

Spatiotemporal features of DCE-MRI for breast cancer diagnosis

BACKGROUND AND OBJECTIVE

Breast cancer is a major cause of mortality among women if not treated in early stages. Previous works developed non-invasive diagnosis methods using imaging data, focusing on specific sets of features that can be called spatial features or temporal features. However, limited set of features carry limited information, requiring complex classification methods to diagnose the disease. For non-invasive diagnosis, different imaging modalities can be used. DCE-MRI is one of the best imaging techniques that provides temporal information about the kinetics of the contrast agent in suspicious lesions along with acceptable spatial resolution.

METHODS

We have extracted and studied a comprehensive set of features from spatiotemporal space to obtain maximum available information from the DCE-MRI data. Then, we have applied a feature fusion technique to remove common information and extract a feature set with maximum information to be used by a simple classification method. We have also implemented conventional feature selection and classification methods and compared them with our proposed approach.

RESULTS

Experimental results obtained from DCE-MRI data of 26 biopsy or short-term follow-up proven patients illustrate that the proposed method outperforms alternative methods. The proposed method achieves a classification accuracy of 99% without missing any of the malignant cases.

CONCLUSIONS

The proposed method may help physicians determine the likelihood of malignancy in breast cancer using DCE-MRI without biopsy.

Influence of B1-Inhomogeneity on Pharmacokinetic Modeling of Dynamic Contrast-Enhanced MRI: A Simulation Study

OBJECTIVE

To simulate the B₁-inhomogeneity-induced variation of pharmacokinetic parameters on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

MATERIALS AND METHODS

B₁-inhomogeneity-induced flip angle (FA) variation was estimated in a phantom study. Monte Carlo simulation was performed to assess the FA-deviation-induced measurement error of the pre-contrast R₁, contrast-enhancement ratio, Gd-concentration, and two-compartment pharmacokinetic parameters (K^trans, v_e, and v_p).

RESULTS

B₁-inhomogeneity resulted in -23-5% fluctuations (95% confidence interval [CI] of % error) of FA. The 95% CIs of FA-dependent % errors in the gray matter and blood were as follows: -16.7-61.8% and -16.7-61.8% for the pre-contrast R₁, -1.0-0.3% and -5.2-1.3% for the contrast-enhancement ratio, and -14.2-58.1% and -14.1-57.8% for the Gd-concentration, respectively. These resulted in -43.1-48.4% error for K^trans, -32.3-48.6% error for the v_e, and -43.2-48.6% error for v_p. The pre-contrast R₁ was more vulnerable to FA error than the contrast-enhancement ratio, and was therefore a significant cause of the Gd-concentration error. For example, a -10% FA error led to a 23.6% deviation in the pre-contrast R₁, -0.4% in the contrast-enhancement ratio, and 23.6% in the Gd-concentration. In a simulated condition with a 3% FA error in a target lesion and a -10% FA error in a feeding vessel, the % errors of the pharmacokinetic parameters were -23.7% for K^trans, -23.7% for v_e, and -23.7% for v_p.

CONCLUSION

Even a small degree of B₁-inhomogeneity can cause a significant error in the measurement of pharmacokinetic parameters on DCE-MRI, while the vulnerability of the pre-contrast R₁ calculations to FA deviations is a significant cause of the miscalculation.

A Pilot Study Evaluating the Use of Dynamic Contrast-Enhanced Perfusion MRI to Predict Local Recurrence After Radiosurgery on Spinal Metastases

Purpose:

Dynamic contrast-enhanced magnetic resonance imaging offers noninvasive characterization of the vascular microenvironment and hemodynamics. Stereotactic radiosurgery, or stereotactic body radiation therapy, engages a vascular component of the tumor response which may be detectable using dynamic contrast-enhanced magnetic resonance imaging. The purpose of this study is to examine whether dynamic contrast-enhanced magnetic resonance imaging can be used to predict local tumor recurrence in patients with spinal bone metastases who undergo high-dose radiotherapy with stereotactic radiosurgery.

Materials and Methods:

We conducted a study of 30 patients with spinal metastases who underwent dynamic contrast-enhanced magnetic resonance imaging before and after radiotherapy. Twenty patients received single-fraction stereotactic radiosurgery (24 Gy), while 10 received hypofractionated stereotactic radiosurgery (3-5 fractions, 27-30 Gy total). Kaplan-Meier analysis was used to estimate the actuarial local recurrence rates. Two perfusion parameters (K^trans: permeability and V_p: plasma volume) were measured for each metastasis. Percentage change in parameter values from pre- to posttreatment was calculated and compared.

Results:

At 20-month median follow-up, 5 of the 30 patients had pathological evidence of local recurrence. One- and 3-year actuarial local recurrence rates were 24% and 44% for the hypofractionated stereotactic radiosurgery cohort versus 5% and 16% for the single-fraction stereotactic radiosurgery cohort (P = .20). The average change in V_p and K^trans for patients without local recurrence versus those with local recurrence was −76% and −66% versus +28% and −14% (P < .01 for both). With a cutoff point of −20%, V_p had a sensitivity, specificity, positive predictive value, and negative predictive value of 100%, 98%, 91%, and 100%, respectively, for the detection of local recurrence following high-dose radiotherapy. Using this definition, dynamic contrast-enhanced magnetic resonance imaging identified local recurrence up to 18 months (mean [standard deviation], 6.6 [6.8] months) earlier than standard magnetic resonance imaging.

Conclusions:

We demonstrated that changes in perfusion parameters, particularly V_p, after high-dose radiotherapy to spinal bone metastases were predictive of local tumor recurrence. These changes predicted local recurrence on average >6 months earlier than standard imaging did.

Dynamic fractal signature dissimilarity analysis for therapeutic response assessment using dynamic contrast-enhanced MRI

PURPOSE

To develop a dynamic fractal signature dissimilarity (FSD) method as a novel image texture analysis technique for the quantification of tumor heterogeneity information for better therapeutic response assessment with dynamic contrast-enhanced (DCE)-MRI.

METHODS

A small animal antiangiogenesis drug treatment experiment was used to demonstrate the proposed method. Sixteen LS-174T implanted mice were randomly assigned into treatment and control groups (n = 8/group). All mice received bevacizumab (treatment) or saline (control) three times in two weeks, and one pretreatment and two post-treatment DCE-MRI scans were performed. In the proposed dynamic FSD method, a dynamic FSD curve was generated to characterize the heterogeneity evolution during the contrast agent uptake, and the area under FSD curve (AUCFSD) and the maximum enhancement (MEFSD) were selected as representative parameters. As for comparison, the pharmacokinetic parameter K(trans) map and area under MR intensity enhancement curve AUCMR map were calculated. Besides the tumor's mean value and coefficient of variation, the kurtosis, skewness, and classic Rényi dimensions d1 and d2 of K(trans) and AUCMR maps were evaluated for heterogeneity assessment for comparison. For post-treatment scans, the Mann-Whitney U-test was used to assess the differences of the investigated parameters between treatment/control groups. The support vector machine (SVM) was applied to classify treatment/control groups using the investigated parameters at each post-treatment scan day.

RESULTS

The tumor mean K(trans) and its heterogeneity measurements d1 and d2 values showed significant differences between treatment/control groups in the second post-treatment scan. In contrast, the relative values (in reference to the pretreatment value) of AUCFSD and MEFSD in both post-treatment scans showed significant differences between treatment/control groups. When using AUCFSD and MEFSD as SVM input for treatment/control classification, the achieved accuracies were 93.8% and 93.8% at first and second post-treatment scan days, respectively. In comparison, the classification accuracies using d1 and d2 of K(trans) map were 87.5% and 100% at first and second post-treatment scan days, respectively.

CONCLUSIONS

As quantitative metrics of tumor contrast agent uptake heterogeneity, the selected parameters from the dynamic FSD method accurately captured the therapeutic response in the experiment. The potential application of the proposed method is promising, and its addition to the existing DCE-MRI techniques could improve DCE-MRI performance in early assessment of treatment response.

Difference in the intratumoral distributions of extracellular-fluid and intravascular MR contrast agents in glioblastoma growth

Contrast enhancement by an extracellular-fluid contrast agent (CA) (Gd-DOTA) depends primarily on the blood-brain-barrier permeability (bp), and transverse-relaxation change caused by intravascular T₂ CA (superparamagnetic iron oxide nanoparticles, SPIONs) is closely associated with the blood volume (BV). Pharmacokinetic (PK) vascular characterization based on single-CA-using dynamic contrast-enhanced MRI (DCE-MRI) has shown significant measurement variation according to the molecular size of the CA. Based on this recognition, this study used a dual injection of Gd-DOTA and SPIONs for tracing the changes of bp and BV in C6 glioma growth (Days 1 and 7 after the tumor volume reached 2 mL). bp was quantified according to the non-PK parameters of Gd-DOTA-using DCE-MRI (wash-in rate, maximum enhancement ratio and initial area under the enhancement curve (IAUC)). BV was estimated by SPION-induced ΔR₂ * and ΔR₂ . With validated measurement reliability of all the parameters (coefficients of variation ≤10%), dual-contrast MRI demonstrated a different region-oriented distribution between Gd-DOTA and SPIONs within a tumor as follows: (a) the BV increased stepwise from the tumor center to the periphery; (b) the tumor periphery maintained the augmented BV to support continuous tumor expansion from Day 1 to Day 7; (c) the internal tumor area underwent significant vascular shrinkage (i.e. decreased ΔR₂ and ΔR₂ ) as the tumor increased in size; (d) the tumor center showed greater bp-indicating parameters, i.e. wash-in rate, maximum enhancement ratio and IAUC, than the periphery on both Days 1 and 7 and (e) the tumor center showed a greater increase of bp than the tumor periphery in tumor growth, as suggested to support tumor viability when there is insufficient blood supply. In the MRI-histologic correlation, a prominent BV increase in the tumor periphery seen in MRI was verified with increased fluorescein isothiocyanate-dextran signals and up-regulated immunoreactivity of CD31-VEGFR. In conclusion, the spatiotemporal alterations of BV and bp in glioblastoma growth, i.e. augmented BV in the tumor periphery and increased bp in the center, can be sufficiently evaluated by MRI with dual injection of extracellular-fluid Gd chelates and intravascular SPION.

Assessment of Treatment Response With Diffusion-Weighted MRI and Dynamic Contrast-Enhanced MRI in Patients With Early-Stage Breast Cancer Treated With Single-Dose Preoperative Radiotherapy: Initial Results

Single-dose preoperative stereotactic body radiotherapy is a novel radiotherapy technique for the early-stage breast cancer, and the treatment response pattern of this technique needs to be investigated on a quantitative basis. In this work, dynamic contrast-enhanced magnetic resonance imaging and diffusion-weighted magnetic resonance imaging were used to study the treatment response pattern in a unique cohort of patients with early-stage breast cancer treated with preoperative radiation. Fifteen female qualified patients received single-dose preoperative radiotherapy with 1 of the 3 prescription doses: 15 Gy, 18 Gy, and 21 Gy. Magnetic resonance imaging scans including both diffusion-weighted magnetic resonance imaging and dynamic contrast-enhanced magnetic resonance imaging were acquired before radiotherapy for planning and after radiotherapy but before surgical resection. In diffusion-weighted magnetic resonance imaging, the regional averaged apparent diffusion coefficient was calculated. In dynamic contrast-enhanced magnetic resonance imaging, quantitative parameters K (trans) and v e were evaluated using the standard Tofts model based on the average contrast agent concentration within the region of interest, and the semiquantitative initial area under the concentration curve (iAUC6min) was also recorded. These parameters' relative changes after radiotherapy were calculated for gross tumor volume, clinical target volume, and planning target volume. The initial results showed that after radiotherapy, initial area under the concentration curve significantly increased in planning target volume (P < .006) and clinical target volume (P < .006), and v e significantly increased in planning target volume (P < .05) and clinical target volume (P < .05). Statistical studies suggested that linear correlations between treatment dose and the observed parameter changes exist in most examined tests, and among these tests, the change in gross tumor volume regional averaged apparent diffusion coefficient (P < .012) and between treatment dose and planning target volume K (trans) (P < .029) were found to be statistically significant. Although it is still preliminary, this pilot study may be useful to provide insights for future works.

Accelerated Brain DCE-MRI Using Iterative Reconstruction With Total Generalized Variation Penalty for Quantitative Pharmacokinetic Analysis: A Feasibility Study

PURPOSE

To investigate the feasibility of using undersampled k-space data and an iterative image reconstruction method with total generalized variation penalty in the quantitative pharmacokinetic analysis for clinical brain dynamic contrast-enhanced magnetic resonance imaging.

METHODS

Eight brain dynamic contrast-enhanced magnetic resonance imaging scans were retrospectively studied. Two k-space sparse sampling strategies were designed to achieve a simulated image acquisition acceleration factor of 4. They are (1) a golden ratio-optimized 32-ray radial sampling profile and (2) a Cartesian-based random sampling profile with spatiotemporal-regularized sampling density constraints. The undersampled data were reconstructed to yield images using the investigated reconstruction technique. In quantitative pharmacokinetic analysis on a voxel-by-voxel basis, the rate constant K^trans in the extended Tofts model and blood flow F_B and blood volume V_B from the 2-compartment exchange model were analyzed. Finally, the quantitative pharmacokinetic parameters calculated from the undersampled data were compared with the corresponding calculated values from the fully sampled data. To quantify each parameter's accuracy calculated using the undersampled data, error in volume mean, total relative error, and cross-correlation were calculated.

RESULTS

The pharmacokinetic parameter maps generated from the undersampled data appeared comparable to the ones generated from the original full sampling data. Within the region of interest, most derived error in volume mean values in the region of interest was about 5% or lower, and the average error in volume mean of all parameter maps generated through either sampling strategy was about 3.54%. The average total relative error value of all parameter maps in region of interest was about 0.115, and the average cross-correlation of all parameter maps in region of interest was about 0.962. All investigated pharmacokinetic parameters had no significant differences between the result from original data and the reduced sampling data.

CONCLUSION

With sparsely sampled k-space data in simulation of accelerated acquisition by a factor of 4, the investigated dynamic contrast-enhanced magnetic resonance imaging pharmacokinetic parameters can accurately estimate the total generalized variation-based iterative image reconstruction method for reliable clinical application.

Association between dynamic contrast enhanced MRI imaging features and WHO histopathological grade in patients with invasive ductal breast cancer

The present study aimed to investigate the dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) and World Health Organization (WHO) histopathological grade in patients with invasive ductal breast cancer. A retrospective analysis on the results of DCE-MRI of 92 patients, who were diagnosed with invasive ductal breast cancer following surgery or biopsy, and these results were correlated with WHO histopathological grade. The statistical analysis demonstrated that the tumor size, shape and characteristics of early enhancement were associated with the WHO histopathological grade: The larger the lesion's long diameter, the higher the WHO histopathological grade; the WHO histopathological grades of round and oval masses were relatively lower, while those of lobulated and irregular masses were higher; and tumors with heterogeneous and ring-like enhancement exhibited higher WHO histopathological grades, while those of homogeneous enhancement were lower. The lesion's margin shape was not associated with the WHO histopathological grade. The present study demonstrates that features of DCE-MRI and WHO histopathological grade in patients with invasive ductal breast cancer are correlated, and these MRI features could be used to evaluate the biological behavior and prognosis of lesions.

Dynamic contrast-enhanced MRI for oncology drug development

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is a promising tool for evaluating tumor vascularity, as it can provide vasculature-derived, functional, and quantitative parameters. To implement DCE-MRI parameters as biomarkers for monitoring the effect of antiangiogenic or vascular-disrupting treatment, two crucial elements of surrogate endpoint, ie, validation and qualification, should be satisfied. Although early studies have shown the accuracy and reliability of DCE-MRI parameters for evaluating treatment-driven vascular alterations, there have been an increasing number of studies demonstrating the limitations of DCE-MRI parameters as surrogate endpoints. Therefore, in order to improve the application of DCE-MRI parameters in drug development, it is necessary to establish a standardized evaluation method and to determine the correct therapeutics-oriented meaning of individual DCE-MRI parameter. In this regard, this article describes the biophysical background and data acquisition/analysis techniques of DCE-MRI while focusing on the validation and qualification issues. Specifically, the causes of disagreement and confusion encountered in the preclinical and clinical trials using DCE-MRI are presented in detail. Finally, considering these limitations, we present potential strategies to optimize implementation of DCE-MRI. J. Magn. Reson. Imaging 2016;44:251-264.

A Comparison of Lumpectomy Cavity Delineations Between Use of Magnetic Resonance Imaging and Computed Tomography Acquired With Patient in Prone Position for Radiation Therapy Planning of Breast Cancer

PURPOSE

To compare lumpectomy cavity (LC) and planning target volume (PTV) delineated with the use of magnetic resonance imaging (MRI) and computed tomography (CT) and to examine the possibility of replacing CT with MRI for radiation therapy (RT) planning for breast cancer.

METHODS AND MATERIALS

MRI and CT data were acquired for 15 patients with early-stage breast cancer undergoing lumpectomy during RT simulation in prone positions, the same as their RT treatment positions. The LCs were delineated manually on both CT (LC-CT) and MRI acquired with 4 sequences: T1, T2, STIR, and DCE. Various PTVs were created by expanding a 15-mm margin from the corresponding LCs and from the union of the LCs for the 4 MRI sequences (PTV-MRI). Differences were measured in terms of cavity visualization score (CVS) and dice coefficient (DC).

RESULTS

The mean CVSs for T1, T2, STIR, DCE, and CT defined LCs were 3.47, 3.47, 3.87, 3.50. and 2.60, respectively, implying that the LC is mostly visible with a STIR sequence. The mean reductions of LCs from those for CT were 22%, 43%, 36%, and 17% for T1, T2, STIR, and DCE, respectively. In 14 of 15 cases, MRI (union of T1, T2, STIR, and DCE) defined LC included extra regions that would not be visible from CT. The DCs between CT and MRI (union of T1, T2, STIR, and DCE) defined volumes were 0.65 ± 0.20 for LCs and 0.85 ± 0.06 for PTVs. There was no obvious difference between the volumes of PTV-MRI and PTV-CT, and the average PTV-STIR/PTV-CT volume ratio was 0.83 ± 0.23.

CONCLUSIONS

The use of MRI improves the visibility of LC in comparison with CT. The volumes of LC and PTV generated based on a MRI sequence are substantially smaller than those based on CT, and the PTV-MRI volumes, defined by the union of T1, T2, STIR, and DCE, were comparable with those of PTV-CT for most of the cases studied.

An efficient calculation method for pharmacokinetic parameters in brain permeability study using dynamic contrast-enhanced MRI

PURPOSE

To develop an efficient method for calculating pharmacokinetic (PK) parameters in brain DCE-MRI permeability studies.

METHODS

A linear least-squares fitting algorithm based on a derivative expression of the two-compartment PK model was proposed to analytically solve for the PK parameters. Noise in the expression was minimized through low-pass filtering. Simulation studies were conducted in which the proposed method was compared with two existing methods in terms of accuracy and efficiency. Five in vivo brain studies were demonstrated for potential clinical application.

RESULTS

In the simulation studies using chosen parameter values, the calculated percent difference of K(trans) by the proposed method was <5.0% with a temporal resolution (Δt) < 5 s, and the accuracies of all parameter results were better or comparable to existing methods. When analyzed within certain parameter intensity ranges, the proposed method was more accurate than the existing methods and improved the efficiency by a factor of up to 458 for a Δt = 1 s and up to 38 for a Δt = 5 s. In the in vivo study, the calculated parameters using the proposed method were comparable to those using the existing methods with improved efficiencies.

CONCLUSIONS

An efficient method was developed for the accurate and efficient calculation of parameters in brain DCE-MRI permeability studies.

GPU-accelerated compartmental modeling analysis of DCE-MRI data from glioblastoma patients treated with bevacizumab

The compartment model analysis using medical imaging data is the well-established but extremely time consuming technique for quantifying the changes in microvascular physiology of targeted organs in clinical patients after antivascular therapies. In this paper, we present a first graphics processing unit-accelerated method for compartmental modeling of medical imaging data. Using this approach, we performed the analysis of dynamic contrast-enhanced magnetic resonance imaging data from bevacizumab-treated glioblastoma patients in less than one minute per slice without losing accuracy. This approach reduced the computation time by more than 120-fold comparing to a central processing unit-based method that performed the analogous analysis steps in serial and more than 17-fold comparing to the algorithm that optimized for central processing unit computation. The method developed in this study could be of significant utility in reducing the computational times required to assess tumor physiology from dynamic contrast-enhanced magnetic resonance imaging data in preclinical and clinical development of antivascular therapies and related fields.

Treatment assessment of radiotherapy using MR functional quantitative imaging

Recent developments in magnetic resonance (MR) functional quantitative imaging have made it a potentially powerful tool to assess treatment response in radiation therapy. With its abilities to capture functional information on underlying tissue characteristics, MR functional quantitative imaging can be valuable in assessing treatment response and as such to optimize therapeutic outcome. Various MR quantitative imaging techniques, including diffusion weighted imaging, diffusion tensor imaging, MR spectroscopy and dynamic contrast enhanced imaging, have been investigated and found useful for assessment of radiotherapy. However, various aspects including data reproducibility, interpretation of biomarkers, image quality and data analysis impose challenges on applications of MR functional quantitative imaging in radiotherapy assessment. All of these challenging issues shall be addressed to help us understand whether MR functional quantitative imaging is truly beneficial and contributes to future development of radiotherapy. It is evident that individualized therapy is the future direction of patient care. MR functional quantitative imaging might serves as an indispensable tool towards this promising direction.