Renal artery assessment with nonenhanced steady-state free precession versus contrast-enhanced MR angiography

PURPOSE: To prospectively assess the diagnostic accuracy of nonenhanced three-dimensional (3D) steady-state free precession (SSFP) magnetic resonance (MR) angiography for detection of renal artery stenosis (RAS), with breath-hold contrast material-enhanced MR angiography performed as the reference standard. MATERIALS AND METHODS: The study was local ethics committee approved; all patients gave written informed consent. Fifty-three patients (30 male, 23 female; mean age, 58 years) with arterial hypertension and suspected of having RAS were examined with 1.5-T 3D SSFP renal MR angiography. Stenosis grade, maximal visible vessel length, and subjective image quality were compared. Sensitivity, specificity, accuracy, and negative predictive value (NPV) were calculated on artery-by-artery and patient-by-patient bases. The significance of the results was assessed with the paired two-sided t test for continuous variables and with the marginal homogeneity test for categorical variables. Cohen kappa statistics were used to estimate interobserver agreement. RESULTS: One hundred eight renal arteries with 20 significant (>or=50%) stenoses were detected with contrast-enhanced MR angiography. At artery-by-artery analysis, sensitivity, specificity, accuracy, and NPV of nonenhanced SSFP MR angiography for RAS detection were 100%, 93%, 94%, and 100%, respectively, for observer 1 and 95%, 95%, 95%, and 99%, respectively, for observer 2. Corresponding patient-by-patient values were 100%, 92%, 94%, and 100%, respectively, for observer 1 and 100%, 95%, 96%, and 100%, respectively, for observer 2. Overestimation of stenosis grade with SSFP MR angiography resulted in six and four false-positive findings for readers 1 and 2, respectively. Mean maximal visible lengths of the renal arteries were 69.9 mm at contrast-enhanced MR angiography and 61.1 mm at SSFP MR angiography (P<.001). Both techniques yielded good to excellent image quality. CONCLUSION: Slab-selective inversion-prepared 3D SSFP MR angiography had high sensitivity, specificity, accuracy, and NPV for RAS detection, without the need for contrast material. However, RAS severity was overestimated in some patients.

Delayed gadolinium enhanced MRI of cartilage (dGEMRIC): Molekulare MRT-Bildgebung des Hüftgelenkknorpels [Delayed gadolinium enhanced MRI of cartilage (dGEMRIC): molecular MRI of hip joint cartilage]

Factors such as instability and impingement lead to early cartilage damage and osteoarthritis of the hip joint. The surgical outcome of joint-preserving surgery about the hip joint depends on the preoperative quality of joint cartilage.For in vivo evaluation of cartilage quality, different biochemically sensitive magnetic resonance imaging (MRI) procedures have been tested, some of which have the potential of inducing a paradigm shift in the evaluation and treatment of cartilage damage and early osteoarthritis.Instead of reacting to late sequelae in a palliative way, physicians could assess cartilage damage early on, and the treatment intensity could be adequate and based on the disease stage. Furthermore, the efficiency of different therapeutic interventions could be evaluated and monitored.This article reviews the recent application of delayed gadolinium-enhanced MRI of cartilage (dGEMRIC) and discusses its use for assessing cartilage quality in the hip joint. dGEMRIC is more sensitive to early cartilage changes in osteoarthritis than are radiographic measures and might be a helpful tool for assessing cartilage quality.

Assessment of osteoarthritis in hips with femoroacetabular impingement using delayed gadolinium enhanced MRI of cartilage

PURPOSE: To determine the feasibility of assessing early osteoarthritis (OA) in hips with femoroacetabular impingement (FAI) using delayed Gadolinium enhanced MRI of Cartilage (dGEMRIC). MATERIALS AND METHODS: Thirty-seven hips in 30 patients who had a dGEMRIC scan and radiographic evidence of FAI were identified. Clinical symptoms were assessed. Radiographic measurements were performed to determine acetabular and femoral morphology. The severity of radiographic OA was determined using Tönnis grade and minimum joint space width (JSW). On MRI, the alpha angle was measured on the sagittal oblique slices. Correlations between dGEMRIC index, patient symptoms, morphologic measurements, radiographic OA, and age were determined. RESULTS: Significant correlations were observed between dGEMRIC index, pain (P < 0.05), and alpha angle (P < 0.05). The correlation of dGEMRIC with alpha angle suggests that hips with more femoral deformity show signs of early OA. CONCLUSION: The results of osteoplasty for FAI depend on the amount of pre-existing OA in the joint. dGEMRIC may be a useful technique for diagnosis and staging of early osteoarthritis in hips with impingement.

Effect of vascular targeting agent in rat tumor model: dynamic contrast-enhanced versus diffusion-weighted MR imaging

PURPOSE: To compare dynamic contrast material-enhanced magnetic resonance (MR) imaging and diffusion-weighted MR imaging for noninvasive evaluation of early and late effects of a vascular targeting agent in a rat tumor model. MATERIALS AND METHODS: The study protocol was approved by the local ethics committee for animal care and use. Thirteen rats with one rhabdomyosarcoma in each flank (26 tumors) underwent dynamic contrast-enhanced imaging and diffusion-weighted echo-planar imaging in a 1.5-T MR unit before intraperitoneal injection of combretastatin A4 phosphate and at early (1 and 6 hours) and later (2 and 9 days) follow-up examinations after the injection. Histopathologic examination was performed at each time point. The apparent diffusion coefficient (ADC) of each tumor was calculated separately on the basis of diffusion-weighted images obtained with low b gradient values (ADC(low); b = 0, 50, and 100 sec/mm(2)) and high b gradient values (ADC(high); b = 500, 750, and 1000 sec/mm(2)). The difference between ADC(low) and ADC(high) was used as a surrogate measure of tissue perfusion (ADC(low) - ADC(high) = ADC(perf)). From the dynamic contrast-enhanced MR images, the volume transfer constant k and the initial slope of the contrast enhancement-time curve were calculated. For statistical analyses, a paired two-tailed Student t test and linear regression analysis were used. RESULTS: Early after administration of combretastatin, all perfusion-related parameters (k, initial slope, and ADC(perf)) decreased significantly (P < .001); at 9 days after combretastatin administration, they increased significantly (P < .001). Changes in ADC(perf) were correlated with changes in k (R(2) = 0.46, P < .001) and the initial slope (R(2) = 0.67, P < .001). CONCLUSION: Both dynamic contrast-enhanced MR imaging and diffusion-weighted MR imaging allow monitoring of perfusion changes induced by vascular targeting agents in tumors. Diffusion-weighted imaging provides additional information about intratumoral cell viability versus necrosis after administration of combretastatin.

Contrast-enhanced MR myelography in spontaneous intracranial hypotension: description of an artefact imitating CSF leakage

In contrast-enhanced (CE) MR myelography, hyperintense signal outside the intrathecal space in T1-weighted sequences with spectral presaturation inversion recovery (SPIR) is usually considered to be due to CSF leakage. We retrospectively investigated a hyperintense signal at the apex of the lung appearing in this sequence in patients with SIH (n = 5), CSF rhinorrhoea (n = 2), lumbar spine surgery (n = 1) and in control subjects (n = 6). Intrathecal application of contrast agent was performed in all patients before MR examination, but not in the control group. The reproducible signal increase was investigated with other fat suppression techniques and MR spectroscopy. All patients and controls showed strongly hyperintense signal at the apex of the lungs imitating CSF leakage into paraspinal tissue. This signal increase was identified as an artefact, caused by spectroscopically proven shift and broadening of water and lipid resonances (1-2 ppm) in this anatomical region. Only patients with SIH showed additional focal enhancement along the spinal nerve roots and/or in the spinal epidural space. In conclusion CE MR myelography with spectral selective fat suppression shows a reproducible cervicothoracic artefact, imitating CSF leakage. Selective water excitation technique as well as periradicular and epidural contrast collections may be helpful to discriminate between real pathological findings and artefacts.

5-ALA complete resections go beyond MR contrast enhancement: shift corrected volumetric analysis of the extent of resection in surgery for glioblastoma

BACKGROUND The technique of 5-aminolevulinic acid (5-ALA) tumor fluorescence is increasingly used to improve visualization of tumor tissue and thereby to increase the rate of patients with gross total resections. In this study, we measured the resection volumes in patients who underwent 5-ALA-guided surgery for non-eloquent glioblastoma and compared them with the preoperative tumor volume. METHODS We selected 13 patients who had received a complete resection according to intraoperative 5-ALA induced fluorescence and CRET according to post-operative T1 contrast-enhanced MRI. The volumes of pre-operative contrast enhancing tissue, post-operative resection cavity and resected tissue were determined through shift-corrected volumetric analysis. RESULTS The mean resection cavity (29 cm(3)) was marginally smaller than the pre-operative contrast-enhancing tumor (39 cm(3), p = 0.32). However, the mean overall resection volume (84 cm(3)) was significantly larger than the pre-operative contrast-enhancing tumor (39 cm(3), p = 0.0087). This yields a mean volume of resected 5-ALA positive, but radiological non-enhancing tissue of 45 cm(3). The mean calculated rim of resected tissue surpassed pre-operative tumor diameter by 6 mm (range 0-10 mm). CONCLUSIONS Results of the current study imply that (i) the resection cavity underestimates the volume of resected tissue and (ii) 5-ALA complete resections go significantly beyond the volume of pre-operative contrast-enhancing tumor bulk on MRI, indicating that 5-ALA also stains MRI non-enhancing tumor tissue. Use of 5-ALA may thus enable extension of coalescent tumor resection beyond radiologically evident tumor. The impact of this more extended resection method on time to progression and overall survival has not been determined, and potentially puts adjacent and functionally intact tissue at risk.

Contrast-enhanced ultrasound for myocardial perfusion imaging.

Ultrasound contrast agents are gas-filled microbubbles that enhance visualization of cardiac structures, function and blood flow during contrast-enhanced ultrasound (CEUS). An interesting cardiovascular application of CEUS is myocardial contrast echocardiography, which allows real-time myocardial perfusion imaging. The intraoperative use of this technically challenging imaging method is limited at present, although several studies have examined its clinical utility during cardiac surgery in the past. In the present review we provide general information on the basic principles of CEUS and discuss the methodology and technical aspects of myocardial perfusion imaging.

3D morphometry using automated aortic segmentation in native MR angiography: an alternative to contrast enhanced MRA?

INTRODUCTION Native-MR angiography (N-MRA) is considered an imaging alternative to contrast enhanced MR angiography (CE-MRA) for patients with renal insufficiency. Lower intraluminal contrast in N-MRA often leads to failure of the segmentation process in commercial algorithms. This study introduces an in-house 3D model-based segmentation approach used to compare both sequences by automatic 3D lumen segmentation, allowing for evaluation of differences of aortic lumen diameters as well as differences in length comparing both acquisition techniques at every possible location. METHODS AND MATERIALS Sixteen healthy volunteers underwent 1.5-T-MR Angiography (MRA). For each volunteer, two different MR sequences were performed, CE-MRA: gradient echo Turbo FLASH sequence and N-MRA: respiratory-and-cardiac-gated, T2-weighted 3D SSFP. Datasets were segmented using a 3D model-based ellipse-fitting approach with a single seed point placed manually above the celiac trunk. The segmented volumes were manually cropped from left subclavian artery to celiac trunk to avoid error due to side branches. Diameters, volumes and centerline length were computed for intraindividual comparison. For statistical analysis the Wilcoxon-Signed-Ranked-Test was used. RESULTS Average centerline length obtained based on N-MRA was 239.0±23.4 mm compared to 238.6±23.5 mm for CE-MRA without significant difference (P=0.877). Average maximum diameter obtained based on N-MRA was 25.7±3.3 mm compared to 24.1±3.2 mm for CE-MRA (P<0.001). In agreement with the difference in diameters, volumes obtained based on N-MRA (100.1±35.4 cm(3)) were consistently and significantly larger compared to CE-MRA (89.2±30.0 cm(3)) (P<0.001). CONCLUSIONS 3D morphometry shows highly similar centerline lengths for N-MRA and CE-MRA, but systematically higher diameters and volumes for N-MRA.

MR Angiography at 3 Tesla to Assess Proximal Internal Carotid Artery Stenoses: Contrast-Enhanced or 3D Time-of-Flight MR Angiography?

PURPOSE The aim of this study was to compare the diagnostic accuracy of 3D time-of-flight (TOF-MRA) and contrast-enhanced (CE-MRA) magnetic resonance angiography at 3 T for detection and quantification of proximal high-grade stenosis using multidetector computed tomography angiography (MDCTA) as reference standard. METHODS The institutional ethics committee approved this prospective study. A total of 41 patients suspected of having internal carotid artery (ICA) stenosis underwent both MDCTA and MRA. CE-MRA and TOF-MRA were performed using a 3.0-T imager with a dedicated eight-element cervical coil. ICA stenoses were measured according to the North American Symptomatic Carotid Endarterectomy Trial criteria and categorized as 0-25 % (minimal), 25-50 % (mild), 50-69 % (moderate), 70-99 % (high grade), and 100 % (occlusion). Sensitivity and specificity for the detection of high-grade ICA stenoses (70-99 %) and ICA occlusions were determined. In addition, intermodality agreement was assessed with κ-statistics for detection of high-grade ICA stenoses (70-99 %) and ICA occlusions. RESULTS A total of 80 carotid arteries of 41 patients were reviewed. Two previously stented ICAs were excluded from analysis. On MDCTA, 7 ICAs were occluded, 12 ICAs presented with and 63 without a high-grade ICA stenosis (70-99 %). For detecting 70-99 % stenosis, both 3D TOF-MRA and CE-MRA were 91.7 % sensitive and 98.5 % specific, respectively. Both MRA techniques were highly sensitive (100 %), and specific (CE-MRA, 100 %; TOF-MRA, 98.7 %) for the detection of ICA occlusion. However, TOF-MRA misclassified one high-grade stenosis as occlusion. Intermodality agreement for detection of 70-99 % ICA stenoses was excellent between TOF-MRA and CE-MRA [κ = 0.902, 95 % confidence interval (CI) = 0.769-1.000], TOF-MRA and MDCTA (κ = 0.902, 95 % CI = 0.769-1.000), and CE-MRA and MDCTA (κ = 0.902, 95 % CI = 0.769-1.000). CONCLUSION Both 3D TOF-MRA and CE-MRA at 3 T are reliable tools for detecting high-grade proximal ICA stenoses (70-99 %). 3D TOF-MRA might misclassify pseudo-occlusions as complete occlusions. If there are no contraindications for CE-MRA, CE-MRA is recommended as primary MR imaging modality.

Evaluation of liver parenchyma and perfusion using dynamic contrast-enhanced computed tomography and contrast-enhanced ultrasonography in captive green iguanas (Iguana iguana) under general anesthesia.

BACKGROUND Contrast-enhanced diagnostic imaging techniques are considered useful in veterinary and human medicine to evaluate liver perfusion and focal hepatic lesions. Although hepatic diseases are a common occurrence in reptile medicine, there is no reference to the use of contrast-enhanced ultrasound (CEUS) and contrast-enhanced computed tomography (CECT) to evaluate the liver in lizards. Therefore, the aim of this study was to evaluate the pattern of change in echogenicity and attenuation of the liver in green iguanas (Iguana iguana) after administration of specific contrast media. RESULTS An increase in liver echogenicity and density was evident during CEUS and CECT, respectively. In CEUS, the mean ± SD (median; range) peak enhancement was 19.9% ± 7.5 (18.3; 11.7-34.6). Time to peak enhancement was 134.0 ± 125.1 (68.4; 59.6-364.5) seconds. During CECT, first visualization of the contrast medium was at 3.6 ± 0.5 (4; 3-4) seconds in the aorta, 10.7 ± 2.2 (10.5; 7-14) seconds in the hepatic arteries, and 15 ± 4.5 (14.5; 10-24) seconds in the liver parenchyma. Time to peak was 14.1 ± 3.4 (13; 11-21) and 31 ± 9.6 (29; 23-45) seconds in the aorta and the liver parenchyma, respectively. CONCLUSION CEUS and dynamic CECT are practical means to determine liver hemodynamics in green iguanas. Distribution of contrast medium in iguana differed from mammals. Specific reference ranges of hepatic perfusion for diagnostic evaluation of the liver in iguanas are necessary since the use of mammalian references may lead the clinician to formulate incorrect diagnostic suspicions.

Evaluation of treatment response after chemoembolisation (TACE) in hepatocellular carcinoma using real time image fusion of contrast-enhanced ultrasound (CEUS) and computed tomography (CT)--preliminary results.

OBJECTIVE To evaluate treatment response of hepatocellular carcinoma (HCC) after transarterial chemoembolization (TACE) with a new real-time imaging fusion technique of contrast-enhanced ultrasound (CEUS) with multi-slice detection computed tomography (CT) in comparison to conventional post-interventional follow-up. MATERIAL AND METHODS 40 patients with HCC (26 male, ages 46-81 years) were evaluated 24 hours after TACE using CEUS with ultrasound volume navigation and image fusion with CT compared to non-enhanced CT and follow-up contrast-enhanced CT after 6-8 weeks. Reduction of tumor vascularization to less than 25% was regarded as "successful" treatment, whereas reduction to levels >25% was considered as "partial" treatment response. Homogenous lipiodol retention was regarded as successful treatment in non-enhanced CT. RESULTS Post-interventional image fusion of CEUS with CT was feasible in all 40 patients. In 24 patients (24/40), post-interventional image fusion with CEUS revealed residual tumor vascularity, that was confirmed by contrast-enhanced CT 6-8 weeks later in 24/24 patients. In 16 patients (16/40), post-interventional image fusion with CEUS demonstrated successful treatment, but follow-up CT detected residual viable tumor (6/16). Non-enhanced CT did not identify any case of treatment failure. Image fusion with CEUS assessed treatment efficacy with a specificity of 100%, sensitivity of 80% and a positive predictive value of 1 (negative predictive value 0.63). CONCLUSIONS Image fusion of CEUS with CT allows a reliable, highly specific post-interventional evaluation of embolization response with good sensitivity without any further radiation exposure. It can detect residual viable tumor at early state, resulting in a close patient monitoring or re-therapy.

New diagnostic and therapeutic approaches to treat ventricular tachycardias originating at the summit of the left ventricle role of merged hemodynamic-MRI and alternative ablation sources

The left ventricular (LV) summit is the most common site of idiopathic epicardial LV arrhythmias and frequently represents a diagnostic and a therapeutic challenge.1 We present a case of sustained monomorphic ventricular tachycardia (SMVT) originating at the LV summit that underwent failed cryosurgical epicardial ablation and was successfully treated with the aid of merged hemodynamic and contrast-enhanced MRI (CE-MRI).

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.