Radiation-induced parotid damage is one of the most common complications in patients with nasopharyngeal carcinoma (NPC) undergoing radiotherapy (RT). Intravoxel incoherent motion (IVIM) magnetic resonance (MR) imaging has been reported for evaluating irradiated parotid damage. However, the changes of IVIM perfusion-related parameters in irradiated parotid glands have not been confirmed by conventional perfusion measurements obtained from dynamic contrast-enhanced (DCE) MR imaging. The purposes of this study were to monitor radiation-induced parotid damage using IVIM and DCE MR imaging and to investigate the correlations between changes of these MR parameters.
Eighteen NPC patients underwent bilateral parotid T1-weighted, IVIM and DCE MR imaging pre-RT (2 weeks before RT) and post-RT (4 weeks after RT). Parotid volume; IVIM MR parameters, including apparent diffusion coefficient (ADC), pure diffusion coefficient (D), pseudo-diffusion coefficient (D*), and perfusion fraction (f); and DCE MR parameters, including maximum relative enhancement (MRE), time to peak (TTP), Wash in Rate, and the degree of xerostomia were recorded. Correlations of parotid MR parameters with mean radiation dose, atrophy rate and xerostomia degree, as well as the relationships between IVIM and DCE MR parameters, were investigated.
From pre-RT to post-RT, all of the IVIM and DCE MR parameters increased significantly (p < 0.001 for ADC, D, f, MRE, Wash in Rate; p = 0.024 for D*; p = 0.037 for TTP). Change rates of ADC, f and MRE were negatively correlated with atrophy rate significantly (all p < 0.05). Significant correlations were observed between the change rates of D* and MRE (r = 0.371, p = 0.026) and between the change rates of D* and TTP (r = 0.396, p = 0.017). The intra- and interobserver reproducibility of IVIM and DCE MR parameters was good to excellent (intraclass correlation coefficient, 0.633–0.983).
Early radiation-induced changes of parotid glands could be evaluated by IVIM and DCE MR imaging. Certain IVIM and DCE MR parameters were correlated significantly.
Radiotherapy (RT) is the main treatment modality for patients with nasopharyngeal carcinoma (NPC). Radiation-induced parotid damage is one of the most common complications, causing xerostomia, dysphagia and increased risk of dental caries  and severely reduces the life quality of these patients. Over the past few years, intensity-modulated radiation therapy (IMRT) has been introduced for the treatment of NPC to preserve parotid function . However, the parotid glands are sensitive to radiation , and the radiation-induced parotid damage cannot be completely avoided even with IMRT . Therefore, it is important to evaluate radiation-induced parotid damage in a timely manner preferably to preserve the function of parotid glands.
The severity of xerostomia can be evaluated based on the radiation morbidity scoring criteria proposed by the Radiation Therapy Oncology Group (RTOG) . However, this evaluation is subjective and cannot depict the morphological and pathophysiological changes in the irradiated parotid glands. Histological examination, which is the gold standard for the evaluation of radiation-induced parotid damage, is not suitable for routine clinical use due to its invasiveness. Scintigraphy, which reveals functional changes in irradiated parotid glands , involves additional radiation exposure. Magnetic resonance (MR) sialography can noninvasively depict irradiated parotid ductal damage , however without any parenchymal information about the irradiated parotid glands.
To investigate the structural and pathophysiological changes, dynamic contrast-enhanced (DCE) and diffusion-weighted (DW) MR imaging have been used for the evaluation of irradiated parotid glands [7–10]. Changes of the vascular permeability and extra-vascular extra-cellular space (EES) in the irradiated parotid glands can be successfully evaluated using DCE MR parameters (such as the transfer coefficient Ktrans and extra-vascular extra-cellular space ve). However, DCE MR imaging involves intravenous injection of gadolinium-based contrast agents, which cause additional expenditure and incur risks of nephrogenic systematic fibrosis or gadolinium deposition in the brain. DW MR imaging generates an apparent diffusion coefficient (ADC), which is affected by both water molecular diffusion and microvascular perfusion simultaneously.
Intravoxel incoherent motion (IVIM) MR imaging was initially proposed by Le Bihan et al. . The perfusion and diffusion information can be separately extracted using IVIM MR imaging with a number of b values. The signal decay at low b values is primarily attributed to perfusion, while data obtained at high b values are mainly dominated by diffusion . Microcirculation changes can be evaluated by IVIM perfusion-related parameters (perfusion fraction f and pseudo-diffusion coefficient D*), and the pure molecular diffusion (D) can reflect the Brownian movement of water molecules. In recent years, IVIM MR imaging has been widely used for the assessment of organic damage, differential diagnosis of tumours and monitoring of cancer therapy [13–16]. Marzi et al. found that IVIM MR parameters (ADC, ADClow, D and f) of parotid glands significantly changed during RT, and the changes in D values were significantly correlated with mean radiation dose . This pilot study indicated the potential of IVIM MR imaging for the evaluation of radiation-induced damage to the parotid glands. However, the changes of IVIM perfusion-related parameters in this pilot study have not been confirmed by conventional perfusion measurements obtained from DCE MR imaging.
Additionally, the precise correlation between IVIM perfusion-related and DCE MR parameters has not been confirmed. Jia et al. reported a significant correlation between f and DCE MR parameters (including Enhancement Amplitude and Maximum Slope of Increase) in NPC . However, Yuan et al. found no correlations between IVIM and DCE MR parameters in lung neoplasms . To our knowledge, the correlations between parameters derived from IVIM and DCE MR imaging in irradiated parotid glands have never been reported.
Therefore, the purposes of this study were to observe the changes in and the relationships between IVIM and DCE MR parameters of irradiated parotid glands after RT and then to correlate the change rates in parotid IVIM and DCE MR parameters with mean radiation dose, atrophy rate and xerostomia degree.
This study was approved by the institutional review board of our hospital, and all of the patients provided written informed consent. From August 2015 to April 2016, 18 patients (male, 14; female, 4; age, 24–70 years old; mean age, 50.1 ± 10.3 years) with an initial diagnosis of poorly differentiated NPC were prospectively enrolled. The inclusion criteria were: (1) a pathological diagnosis of poorly differentiated NPC through biopsy and readiness to receive IMRT with concurrent chemotherapy in our hospital; (2) scheduled to undergo MR evaluation and follow-up in our hospital; and (3) no any history of allergy to gadolinium contrast agents, with a glomerular filtration rate greater than 30 mL/min to accommodate the injection of contrast agents. The exclusion criteria included: (1) absolute MR examination contraindications, such as cardiac pacemakers, aneurysm clips, artificial cochlea implantation, etc.; (2) a history of parotid disorders, such as parotitis, parotid tumours, etc.; and (3) having received RT to the head and neck region in the past. The patients’ characteristics and the flowchart of this study are shown in Table 1 and Fig. 1, respectively.
All of the patients were treated with IMRT to the nasopharyngeal lesions and the neck lymphatic drainage areas, combined with concurrent chemotherapy (three cycles; nedaplatin 60 mg for each cycle). The total accumulated radiation dose within the tumour region was 70Gy, which was divided into two courses. In the first course, the RT field covered the nasopharyngeal lesions and neck lymphatic drainage areas (25 fractions; 2 Gy for each fraction). In the second course, the RT field was reduced to the tumour area according to the CT reset condition (10 fractions; 2 Gy for each fraction). The RT was administered as one fraction for 1 day, five fractions for 1 week and a total of 35 fractions for 7 weeks. Because the bilateral parotid and submandibular glands were quite close to the planning target volume (PTV), some attempt was undertaken to reduce their radiation doses on the premise of meeting the tumour exposure dose. The mean accumulated radiation dose to the parotid glands was calculated from the treatment planning system of Pinnacle3 (Philips Medical Systems, Fitchburg, WI, USA) and TomoTherapy HiArt (TomoTherapy, Madison, WI, USA). Because the bilateral parotid glands received different radiation doses during the course of RT, the bilateral parotid glands of each patient were analysed separately. The mean total accumulated radiation dose to the parotid glands was 28.4 ± 2.4 Gy after RT, which was less than our hospital limit for parotid radiation dose of 30–35 Gy for 50 % volume. All of the patients underwent two MR examinations within 2 weeks before RT (pre-RT) and 4 weeks after RT (post-RT), and the MR scan protocol remained identical during the course. All of the patients successfully underwent the whole therapy and follow-up MR examinations.
Clinical assessment of xerostomia
The degree of xerostomia in NPC patients was assessed 1 hour prior to each MR examination by a radiation oncologist (X.X., with 10 years of clinical experience in head and neck RT), according to the acute radiation morbidity scoring criteria proposed by the Radiation Therapy Oncology Group (RTOG) : grade 0 is no change over baseline; grade 1 indicates mild mouth dryness or slightly thickened saliva but without alterations in the baseline feeding behaviour, such as the increased use of liquids with meals; grade 2 represents moderate dryness or sticky saliva; grade 3 indicates complete dryness; and grade 4 characterizes acute salivary gland necrosis. The grade of xerostomia in each patient at each time was recorded.
All of the patients were asked to fast for 2 h before each MR examination. A full digital 3.0 T MR scanner (Ingenia, Philips Medical Systems, Best, the Netherlands) was used for the MR examinations, with a 16-channel head/neck phased array coil. Patients were placed in the supine position with the head first. The scanning sequences included: transverse T1-weighted (T1W) imaging, intravoxel incoherent motion (IVIM) MR imaging and dynamic contrast-enhanced (DCE) MR imaging. The total duration of the MR examination was approximately 10 min, 9 s.
T1W imaging was obtained with a turbo spin-echo (TSE) sequence. The other parameters were as follows: repetition time/echo time: 400–675 msec/18 msec; TSE factor: 8; matrix: 276 × 215; slice thickness: 5 mm; slice gap: default; slices: 38; field of view: 22 cm; voxel size: 0.8 mm × 0.92 mm; and number of signals averaged: 2. The duration of T1W imaging was approximately 2 min, 27 s.
IVIM MR imaging was obtained with a single-shot echo-planar imaging (SS-EPI) sequence with spectral presaturation attenuated inversion recovery (SPAIR) fat suppression before the injection of the gadolinium contrast agent. A volume shim covering the region of bilateral parotid glands was used to minimize susceptibility artefacts. The other parameters were as follows: repetition time/echo time: 6000 msec/shortest; matrix: 84 × 104; slice thickness: 4 mm; slice gap: 0.4 mm; slices: 22; field of view:22 cm; voxel size: 2.5 mm × 2.08 mm; and number of signals averaged: 2. A total of 9 b values (0, 25, 50, 75, 100, 150, 200, 500, and 800 s/mm2) were applied in the IVIM MR imaging. The duration of the IVIM MR imaging sequence was approximately 5 min, 6 s.
DCE MR imaging was obtained with a three-dimensional (3D) T1-fast field echo (FFE) sequence. Intravenous bolus injection of gadodiamide (0.2 mL/kg bodyweight, GE Healthcare Ireland, Shanghai, China) was administered at a rate of 3.0 mL/s followed by a 15 mL saline flush using an automatic power injection (Medrad Spectris Solaris EP MR Injector System; One Medrad Drive Indianola, PA, USA). The other parameters were as follows: repetition time/echo time: shortest/shortest; flip angle: 10°; matrix: 232 × 232; slice thickness: 2 mm; slice gap: default; slices: 68; field of view:30 cm; voxel size: 1.3 mm × 1.29 mm; and number of signals averaged: 1. A total of 15 dynamics at an interval of 10.4 s were obtained for each patient. The duration of the DCE MR sequence was approximately 2 min, 36 s.
All of the patients completed all of the MR examinations successfully without any discomfort or side effects.
All the MR images were analysed and measured independently by 2 radiologists (X.X., X.X.X.) with 3 and 11 years of experience, respectively, in head and neck radiology, who were blinded to the clinical information of all of the patients. Averaged values of the two radiologists’ measurements were treated as the final results for the parotid glands.
T1W images were transferred into a workstation (Extended MR WorkSpace 220.127.116.11, Philips Medical Systems, Best, the Netherlands), which was used to calculate the parotid volume owing to its perfect soft tissue contrast. The outline of each parotid gland was drawn slice by slice on T1W images to obtain the area of each slice. The volume of each parotid gland was calculated by the following equation: V = ∑Si × (ST + SG), where V represents the volume of the parotid gland, Si represents the area of the ith slice in the parotid gland, ST represents the slice thickness, and SG represents the slice gap. The atrophy rate of each parotid gland from pre-RT to post-RT was calculated by the following equation: RV = (Vpre–Vpost) / Vpre × 100 %, where RV is the atrophy rate of the parotid gland, and Vpre and Vpost are the parotid volume pre-RT and post-RT, respectively.
IVIM MR images were post-processed using DWI-Tool, developed by Philips in IDL 6.3 (ITT Visual Information Solutions, Boulder, CO, USA). D, D* and f maps were generated using the bi-exponential fit equation : Sb/S0 = (1–f) ∙ exp (−bD) + f ∙ exp [−b (D + D*)], where Sb represents the mean signal intensity at different b values of 25, 50, 75, 100, 150, 200, 500, and 800 s/mm2, S0 represents the mean signal intensity at a b value of 0 s/mm2, f represents the fraction of diffusion linked to microcirculation, exp is exponential, D represents the slow component of diffusion, and D* represents the fast component of diffusion. The ADC map was generated using the mono-exponential fit equation : ln (Sb) = ln (S0)−bADC, where ADC represents the microscopic translational motions, including the pure molecular diffusion and perfusion-related diffusion. A region of interest (ROI) was delineated manually on the largest slice of the IVIM MR images for each parotid gland to include as much as parotid parenchyma without obvious vessels. The ROI was automatically transferred between mono-exponential and bi-exponential models, and the corresponding D, D*, f and ADC values of the parotid glands were obtained. The change rates of D, D*, f and ADC values were calculated by the following equation: RIVIM-PARs = (IVIM-PARpost−IVIM-PARpre) / IVIM-PARpre × 100 %, where RIVIM-PARs is the change rate of D, D*, f and ADC values from pre-RT to post-RT, and IVIM-PARpost and IVIM-PARpre are the D, D*, f and ADC values post-RT and pre-RT, respectively.
DCE MR images were post-processed using the “Basic T1 Perfusion” function on the aforementioned workstation. The time-intensity curve was depicted automatically, and the DCE MR parameters, including the maximum relative enhancement (MRE), time to peak (TTP) and Wash in Rate, were calculated. MRE is the maximal signal enhancement of a pixel of certain dynamic relative to that same pixel with the pre-contrast dynamic. TTP is the time between the time of initial intensity and the time of peak intensity. Wash in Rate is the maximum slope between the time of initial intensity and the time of peak intensity. An ROI was drawn manually on the largest slice of the parotid gland to include as much as parotid parenchyma with visible vessels excluded, and the corresponding MRE, TTP and Wash in Rate were automatically obtained. The change rates of the MRE, TTP and Wash in Rate were calculated by the following equation: RDCE-PARs = (DCE-PARpost−DCE-PARpre) / DCE-PARpre × 100 %, where RDCE-PARs is the change rate of MRE, TTP and Wash in Rate from pre-RT to post-RT, and DCE-PARpost and DCE-PARpre are the MRE, TTP and Wash in Rate post-RT and pre-RT, respectively.
The IVIM and DCE MR parameters were repeatedly measured by the second observer with an interval of 4 weeks between measurements to evaluate intraobserver reproducibility.
Continuous numeric data with normal distribution are reported as the means ± standard deviations (SD). The paired sample t test was used to identify any significant changes of IVIM and DCE MR parameters from pre-RT to post-RT. The differences of averaged bilateral parotid RIVIM-PARs and RDCE-PARs between grade 1 and grade 2 of the post-RT xerostomia degree were analysed using the independent-samples t test. Pearson’s correlation test was used to detect the correlations between the change rates of parotid IVIM and DCE MR parameters and mean radiation dose or atrophy rates, as well as between the change rates of IVIM and DCE MR parameters. The intra- and interobserver reproducibility of IVIM and DCE MR parameters was evaluated by calculating intraclass correlation coefficient (ICC) values. The ICC was between 0 and 1, and the interpretation of ICC was as follows: < 0.20, poor; 0.21–0.40, fair; 0.41–0.60, moderate; 0.61–0.80, good; and > 0.80, excellent . Statistical analysis was performed using SPSS software, version 16.0 (SPSS Inc., Chicago, IL, USA). A two-tailed p values < 0.05 was considered statistically significant.
From pre-RT to post-RT, the volume of bilateral parotid glands significantly decreased from 26.1 ± 5.5 cm3 to 18.9 ± 3.9 cm3, with an atrophy rate of 26.5 ± 10.3 % (p < 0.001). The pre- and post-RT IVIM and DCE MR images of the bilateral parotid glands in one representative NPC patient are shown in Fig. 2.
Changes of IVIM and DCE MR parameters from pre-RT to post-RT
As shown in Table 2, all of the IVIM and DCE MR parameters increased from pre-RT to post-RT significantly (all p < 0.05).
Correlations between changes of IVIM or DCE MR parameters and atrophy rate (and mean radiation dose)
As shown in Table 3 and Fig. 3, the change rates of parotid ADC, f and MRE were negatively correlated with the atrophy rate significantly from pre-RT to post-RT (all p < 0.05). There was no significant correlation between mean radiation dose and any change rate of parotid IVIM or any DCE MR parameter (all p > 0.05).
Relationships between IVIM or DCE MR parameters and xerostomia degree
The average change rates of bilateral parotid IVIM and DCE MR parameters in patients with grade 1 xerostomia degree did not significantly differ from that in patients with grade 2 (all p > 0.05).
Correlations between IVIM and DCE MR parameters
As shown in Table 4, the change rate of D* was correlated with that of MRE (r = 0.371 and p = 0.026) and TTP (r = 0.396 and p = 0.017) significantly from pre-RT to post-RT, and there were no significant correlations between the change rates of other IVIM and DCE MR parameters.
Reproducibility of IVIM and DCE MR parameters
As shown in Table 5, the measurements of most parotid IVIM and DCE MR parameters showed excellent intra- and interobserver agreement (ICC, 0.911–0.983), except it was good for f and D* (ICC, 0.633–0.793).
Xerostomia, which is caused by irradiated parotid damage, is a common complication in NPC patients receiving RT. Morphological and microstructural changes in irradiated parotid glands can be noninvasively evaluated by MR imaging [6, 7, 10]. Tissue perfusion (D*, f) and water molecular diffusion (D) features can be quantitatively characterized by IVIM MR imaging with bi-exponential algorithms , and tissue perfusion information about the microcirculation can be described with semiquantitative DCE MR imaging [21, 22]. In this study, the changes in irradiated parotid glands from pre-RT to post-RT were successfully monitored by IVIM and DCE MR parameters, and correlations between the change rates of IVIM and DCE MR parameters were confirmed.
All of the averaged bilateral parotid IVIM MR parameters (including ADC, D, D* and f) increased significantly from pre-RT to post-RT in this study. Marzi et al. reported significant increases of parotid ADC, D, ADClow, and f values in patients with head and neck cancer from baseline to the completion of RT , consistent with our results. The significant increase of ADC and D values might result from the widespread necrosis of acinar cells induced by RT , which caused lower cell density and an augmentation of water molecular diffusion. Houweling et al. reported a significant increase in ve due to cell loss at 6 weeks after RT in oropharyngeal cancer patients , in accordance to our hypothesis. Marzi et al. attributed the increases in parotid ADClow and f on the same day of the completion of RT to radiation-induced vascular oedema, which caused vasodilation and an increase in blood volume . We speculated that the increases of D* and f values in our study shared the same pathophysiologic mechanism. Although Xu et al. reported that parotid microvascular density decreased at 4 h after RT , we considered the increase of blood volume secondary to vascular oedema to be the main effect of RT in the early phase of radiation-induced parotid damage. Furthermore, Lee et al. documented a significant increase in parotid vascular plasma volume (vp) at 3 months after RT in patients with head and neck cancer, explained by vasodilatation and increased blood volume induced by inflammation ; we share the same opinion as them.
All of the averaged bilateral parotid DCE MR parameters (including MRE, TTP and Wash in Rate) increased significantly from pre-RT to post-RT in this study. The increase of MRE and Wash in Rate might share the same mechanism that caused the increase of D* and f values, that is, vascular oedema and increased blood volume secondary to inflammation after RT. In addition, the augmentation of EES due to cell loss also promoted the accumulation of contrast agent in parotid tissue. Juan et al. found a significantly higher TTP value in irradiated parotid glands compared with non-irradiated glands , which was in agreement with our observation.
From pre-RT to post-RT, significant correlations were found between parotid atrophy rate and the change rates of ADC, f and MRE, while Marzi et al. also reported a significant correlation between the change rate of ADC values and parotid atrophy rate . However, the correlations were negative in our study and were positive in Marzi et al.’s, probably due to the discrepancy in follow-up time points between the two studies. The follow-up time point in our study was 4 weeks after RT, while Marzi et al. chose the same day as the completion of RT. Four weeks after RT, interstitial fibrosis was observed in irradiated parotid glands , which might have contributed to the negative correlation between the change rate of ADC value and the parotid atrophy rate at 4 weeks after RT. Interstitial fibrosis could not fully compensate for the increase in EES secondary to radiation-induced parotid cell loss, but it reduced the augmentation of the ADC value induced by the increased EES. Moreover, greater parotid atrophy, indicating more severe damage, was usually accompanied by more interstitial fibrosis. Therefore, a greater parotid atrophy rate might induce a smaller change rate in ADC value due to the increase in interstitial fibrosis, although Marzi et al. also reported a significant, positive correlation between the change rate of f and the parotid atrophy rate on the same day as the completion of RT , in contradiction of our observations as well. This difference might also have resulted from the different MR examination time points between us. Xu et al. reported that the parotid microvascular density in miniature pigs decreased by approximately 20 and 40 % at 4 h and 2 weeks after irradiation (25 Gy), respectively . According to this finding, we could speculate that a greater parotid atrophy rate might be accompanied by more severe microvascular damage at 4 weeks after RT. In our opinion, radiation-induced vascular oedema was the main effect in the early phase after RT, and the decreased microvascular density could not adequately compensate for and merely reduced the increase in blood volume induced by vascular oedema. Therefore, a greater parotid atrophy rate was accompanied by smaller change rates of f and MRE from pre-RT to post-RT.
Good to excellent intra- and interobserver agreement of parotid IVIM and DCE parameters was confirmed in this study. Excellent intra- and interobserver agreement of parotid IVIM MR parameters (D, D*, f) was reported by Su et al. and Xu et al. in healthy volunteers and in patients with Sjögren’s Syndrome, respectively [20, 27]. However, Patel et al. reported good to excellent reproducibility for ADC, D and f and fair reproducibility for D* on liver IVIM MR imaging . The poor reproducibility of D* might have been due to the underlying breathing motion artefacts on liver MR imaging, which were absent on parotid MR imaging. In contrast, IVIM perfusion-related parameters are mainly affected by low b values (b ≤ 100 s/mm2) . Five low b values (0, 25, 50, 75, 100 s/mm2) were adopted in our study, while Patel et al. used only three low b values (0, 50, 100 s/mm2) for liver IVIM MR imaging, which might be another reason for the poor reproducibility of D* measurement in their study.
Based on our findings in this study, the radiation-induced changes in the parotid microstructure could be reflected by IVIM diffusion-related parameters (ADC, D) noninvasively and quantitatively, and the radiation-induced microvascular changes could be assessed by IVIM perfusion-related parameters (f, D*), while parotid IVIM perfusion-related parameters shared similar change patterns with DCE MR parameters from pre-RT to post-RT, which indicated that IVIM perfusion-related parameters could serve as an alternative to DCE MR parameters in the evaluation of irradiated parotid microvascular damage, especially for those patients with renal insufficiency. Both IVIM and DCE MR imaging could serve as objective modalities for evaluating irradiated parotid damage, rather than subjective evaluation of the degree of xerostomia. The change rates of parotid ADC, f and MRE values were negatively correlated with the atrophy rate of the parotid gland from pre-RT to post-RT. It was reported that the decreased parotid gland volume was significantly correlated with decreased saliva production in patients with head-and-neck cancer undergoing RT . The correlations between IVIM or DCE MR parameters and parotid function indices (such as saliva production) will be investigated in our future study, and long-term follow-up will be performed to explore the prognostic values of IVIM and DCE MR parameters in evaluating the damage of irradiated parotid glands. Based on the findings of IVIM or DCE MR imaging, radiation oncologists could perform early interventions to avoid long-term irreversible damage to the parotid glands in NPC patients undergoing RT.
Our study had several limitations. Firstly, the sample size was relatively small although larger than in some other studies of MR imaging in the evaluation of radiation-induced parotid damage [6, 10], and a larger cohort of patients for the correlation of the change rates of IVIM and DCE MR parameters might have been more reliable. Secondly, the range of the mean parotid radiation dose was relatively limited, which might have obscured the relationships between the mean radiation dose and parotid IVIM or DCE MR parameters. We aim to enrol more patients with other head and neck tumours to enlarge the range of the mean parotid radiation dose in the future. Thirdly, the pathological information of irradiated parotid glands in NPC patients was absent in this study due to its invasiveness. Hence, we should perform animal experiments to confirm our hypotheses.
IVIM and DCE MR imaging can noninvasively evaluate the pathophysiologic changes in irradiated parotid glands, including necrosis of acinar cells, vascular oedema and interstitial fibrosis. The IVIM and DCE MR parameters in irradiated parotid glands shared the same change patterns, and significant correlations were found between the change rates of D* and MRE or TTP. IVIM and DCE MR imaging could serve as objective, quantitative and noninvasive modalities for evaluating irradiated parotid damage.
Apparent diffusion coefficient
Pure molecular diffusion
Extra-vascular extra-cellular space
Intraclass correlation coefficient
Intensity-modulated radiation therapy
Intravoxel incoherent motion
Maximum relative enhancement
4 weeks after radiotherapy
2 weeks before radiotherapy
Change rates of MRE, TTP and Wash in Rate from pre-RT to post-RT
Change rates of D, D*, f and ADC values from pre-RT to post-RT
Region of interest
Time to peak
De Castro GJ, Federico MH. Evaluation, prevention and management of radiotherapy-induced xerostomia in head and neck cancer patients. Curr Opin Oncol. 2006;18:266–70.
Kam MK, Leung SF, Zee B, Chau RM, Suen JJ, Mo R, et al. Prospective randomized study of intensity-modulated radiotherapy on salivary gland function in early-stage nasopharyngeal carcinoma patients. J Clin Oncol. 2007;25:4873–9.
Cox JD, Stetz J, Pajak TF. Toxicity criteria of the Radiation Therapy Oncology Group (RTOG) and the European Organization for Research and Treatment of Cancer (EORTC). Int J Radiat Oncol Biol Phys. 1995;31:1341–6.
Van Acker F, Flamen P, Lambin P, Maes A, Kutcher GJ, Weltens C, et al. The utility of SPECT in determining the relationship between radiation dose and salivary gland dysfunction after radiotherapy. Nucl Med Commun. 2001;22:225–31.
Juan CJ, Chen CY, Jen YM, Liu HS, Liu YJ, Hsueh CJ, et al. Perfusion characteristics of late radiation injury of parotid glands: quantitative evaluation with dynamic contrast-enhanced MRI. Eur Radiol. 2009;19:94–102.
Juan CJ, Cheng CC, Chiu SC, Jen YM, Liu YJ, Chiu HC, et al. Temporal evolution of parotid volume and parotid apparent diffusion coefficient in nasopharyngeal carcinoma patients treated by intensity-modulated radiotherapy investigated by magnetic resonance imaging: a pilot study. PLoS ONE. 2015;10, e137073.
Zhu L, Zhu LJ, Shi H, Wang HH, Yan J, Liu BR, et al. Evaluating early response of cervical cancer under concurrent chemo-radiotherapy by intravoxel incoherent motion MR imaging. BMC Cancer. 2016;16:79.
Jia QJ, Zhang SX, Chen WB, Liang L, Zhou ZG, Qiu QH, et al. Initial experience of correlating parameters of intravoxel incoherent motion and dynamic contrast-enhanced magnetic resonance imaging at 3.0 T in nasopharyngeal carcinoma. Eur Radiol. 2014;24:3076–87.
Yuan M, Zhang YD, Zhu C, Yu TF, Shi HB, Shi ZF, et al. Comparison of intravoxel incoherent motion diffusion-weighted MR imaging with dynamic contrast-enhanced MRI for differentiating lung cancer from benign solitary pulmonary lesions. J Magn Reson Imaging. 2016;43:669–79.
Su GY, Xu XQ, Wang YY, Hu H, Shen J, Hong XN, et al. Feasibility study of using intravoxel incoherent motion MRI to detect parotid gland abnormalities in early-stage Sjögren syndrome patients. J Magn Reson Imaging. 2015;43:1455–61.
Thomassin-Naggara I, Bazot M, Darai E, Callard P, Thomassin J, Cuenod CA. Epithelial ovarian tumors: value of dynamic contrast-enhanced MR imaging and correlation with tumor angiogenesis. Radiology. 2008;248:148–59.
Huang B, Wong CS, Whitcher B, Kwong DL, Lai V, Chan Q, et al. Dynamic contrast-enhanced magnetic resonance imaging for characterising nasopharyngeal carcinoma: comparison of semiquantitative and quantitative parameters and correlation with tumour stage. Eur Radiol. 2013;23:1495–502.
Lee FK, King AD, Kam MK, Ma BB, Yeung DK. Radiation injury of the parotid glands during treatment for head and neck cancer: assessment using dynamic contrast-enhanced MR imaging. Radiat Res. 2011;175:291–6.
Radfar L, Sirois DA. Structural and functional injury in minipig salivary glands following fractionated exposure to 70 Gy of ionizing radiation: an animal model for human radiation-induced salivary gland injury. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2003;96:267–74.
Xu XQ, Su GY, Liu J, Hu H, Hong XN, Shi HB, et al. Intravoxel incoherent motion MR imaging measurements of the bilateral parotid glands at 3.0-T MR: effect of age, gender and laterality in healthy adults. Br J Radiol. 2015;88:20150646.
Teshima K, Murakami R, Tomitaka E, Nomura T, Toya R, Hiraki A, et al. Radiation-induced parotid gland changes in oral cancer patients: correlation between parotid volume and saliva production. Jpn J Clin Oncol. 2009;40:42–6.
We would like to thank the Department of Radiology, Nanjing Drum Tower Hospital, The Affiliated Hospital of Nanjing University Medical School for providing the MR examination.
This work was supported by National Natural Science Foundation of China (81371516, 81501441), Foundation of National Health and Family Planning Commission of China (W201306), Social Development Foundation of Jiangsu Province (BE2015605), Natural Science Foundation of Jiangsu Province (the Youth Foundation, BK20150109 and BK20150102), Research Project of Health and Family Planning Commission of Jiangsu Province (Q201508), Six Talent Peaks Project of Jiangsu Province (2015-WSN-079) and Key Project supported by Medical Science and technology development Foundation, Nanjing Department of Health (YKK15068).
Availability of data and materials
The datasets supporting the conclusions of the article are included within the article.
NZ, CC and XD made substantial contributions to data analysis and drafting the manuscript; ML, SL and TTG had significant roles in the acquisition data and interpretation of data; LJZ, JY and BRL provided all oncological support; WBC carried out the quality control of MR examination and data analysis; JH and ZYZ made substantial contributions to conception and design. XFY and TL had significant roles in revising the manuscript. All authors have read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Written informed consent for publication of their clinical details was obtained from the patients.
Ethics approval and consent to participate
This study was approved by the institutional review board of Nanjing Drum Tower hospital, and all of the patients offered the written informed consents.
Authors and Affiliations
Department of Radiology, Nanjing Drum Tower Hospital, The Affiliated Hospital of Nanjing University Medical School, Nanjing, 210008, China
Nan Zhou, Chen Chu, Xin Dou, Ming Li, Song Liu, Jian He & Zhengyang Zhou
The Comprehensive Cancer Centre of Drum Tower Hospital, Medical School of Nanjing University & Clinical Cancer Institute of Nanjing University, Nanjing, 210008, China
Lijing Zhu, Baorui Liu & Jing Yan
Department of Radiology, Nanjing Drum Tower Hospital Clinical College of Traditional Chinese and Western Medicine, Nanjing University of Chinese Medicine, Nanjing, 210008, China
Philips Healthcare, Shanghai, 200233, China
Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA, 30322, USA
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.