Abstract

The relationships between IVIM and DCE-MRI parameters in AS are not clear. We explore the correlation between intravoxel incoherent motion (IVIM) diffusion weighted imaging (DWI) and dynamic contrast-enhanced (DCE) parameters obtained on MR images in patients with ankylosing spondylitis (AS). Forty-four patients with AS were prospectively examined using a 1.5-T MR system. IVIM DWI was performed with 11 values (range, 0–800 s/mm2) for all patients. The correlation coefficients between IVIM and DCE-MRI parameters were analyzed using Spearman’s method. Our results showed that intra- and interobserver reproducibility were excellent to relatively good (ICC = 0.804–0.981; narrow width of 95% limits of agreement). Moderate positive correlations were observed between pure molecular diffusion (Ds) and maximum enhancement (ME) and relative enhancement (RE) (, ; , , resp.). Perfusion-related diffusion (Df) showed negative moderate correlation with ME (, ). However, no correlation was observed between perfusion fraction () and any parameters of ME, RE, TTP, and BE (r = −0.093–0.213; ). In conclusion, the IVIM parameters, especially f, might play a critical role in detecting the progression of AS, because it can provide more perfusion information compared with DCE-MRI; besides the IVIM MRI is a noninvasive method.

1. Introduction

The prevalence of ankylosing spondylitis (AS) is 0.1–2% of the general population. The AS is a chronic disease with active and inactive stage and it is very important to detect the active stage of AS [1]. Sacroiliitis is the most common manifestation of AS with an upward trend to involve vertebral facet joints and other auxiliary structures. Thanks to recent introductions of dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) and Diffusion weighted imaging (DWI), AS could be diagnosed early, and the active or inactive stages of AS could be detected via observation on sacroiliitis. Moreover, previous studies have shown that DCE-MRI and DWI could predict response to the therapy by antitumour necrosis factor-α (TNF- α) in AS [2, 3].

DCE-MRI allows for estimation of tissue perfusion and permeability [3]. Specifically, using pharmacokinetic models from a DCE-MRI acquisition, semiquantitative and quantitative hemodynamic parameters, including relative enhancement (RE), maximum enhancement (ME), time to peak (TTP), and brevity of enhancement (BE), have been processed and modeled [4, 5]. Moreover, Bane et al. have reported that those parameters of the active group were significantly higher than those of the inactive group with AS and showed that DCE-MRI was valuable in assessing treatment of AS [6]. However, parameters of DCE-MRI are influenced by the injection duration time, circulation time, and contrast dosage. It should be noted that the contrast agents might increase the risk of renal fibrosis [7, 8].

DWI method enables showing the microscopic biologic structure and diffusion of water protons in tissue, without requiring an intravenous contrast agent [9]. DWI, now part of most AS MR protocols, is an effective tool for diagnosing the activity of AS and is used to monitor AS [3, 10]. A novel DWI, intravoxel incoherent motion (IVIM) is increasingly being used in AS diagnosis, with the ability to separately reveal the pure diffusion coefficient () and perfusion-related incoherent microcirculation (, Df) at the same time [11]. It is based on the theory that using a series of multiple values of , DWI signals attributes to perfusion in the capillary network in the low -values range ( < 100–150 s/mm2), and in the high -values range, the diffusion of water protons could be explored in tissue [12, 13].

Many studies have already demonstrated correlations of perfusion-related IVIM parameters and DCE-MRI perfusion pharmacokinetic parameters in brain, hepar, breast, and so on [1417]. Some studies supported the fact that the perfusion-related IVIM parameters were significantly correlated with changes in conventional perfusion techniques in cervical cancer, squamous cell carcinoma, nasopharyngeal carcinoma, and so on [14, 18, 19]. On the contrary, other studies showed that there was no correlation between aforementioned parameters in liver cirrhosis, lung cancer, and so on [20, 21]. However, for AS, a valid theory of the relationships between IVIM and DCE-MRI parameters is still lacking. Therefore, the aim of our study was to evaluate the relationship between IVIM parameters obtained using bicomponent analysis and the perfusion obtained using DCE MR imaging in AS.

2. Patients and Methods

2.1. Patient’s Selection

Our prospective study was reviewed by the ethics committee in our hospital, and the approval with written informed consent was obtained from each participant. From December 2015 to November 2016, all patients who were enrolled our study were undergoing IVIM DW imaging and DCE imaging of the sacroiliac joint. Meanwhile, they also were assayed for blood sedimentation rate (ESR), C-reactive protein (CRP), and BASDAI.

2.2. MR Imaging Techniques

All MR imaging was performed using a 1.5-T MR (Achieva 1.5 T, Philips Healthcare, Best, the Netherlands) and utilizing a 4-channel SENSE-body coil with high-gradient performance (amplitude 33 mT/m, slew rate 80 mT/m/ms). Four standard sequences MR imaging was performed: (A) -weighted turbo-spin-echo (TSE) with fat-saturated [echo time (TE)/repetition time (TR), 20 ms/60 ms], (B) Short tau inversion recovery (STIR) (TE/TR, 60 ms/1437 ms; echo train length, 15), (C) Scan parameters of three-dimensional volumetric interpolated fat-suppressed sequence: TE = 2.2 ms, TR = 4.7 ms, the TR is the time between each echo specific for this readout; two different flip angles (5° for pre- contrast scan, and 12° for dynamic scan); acquisition matrix, 112 × 110; FOV, 400 mm × 250 mm; slice thickness, 4 mm; number of slices, 40; SENSE factor, 0; A precontrast 3D T1 FFE scan was followed by a DCE-MRI; 0.1 mmol/kg body weight of gadopentetate dimeglumine contrast agent (Magnevist; Bayer Healthcare, Berlin, Germany) was administered intravenously at a rate of 3 ml/s, followed by a 20-ml saline flush with high-pressure injector after the acquisition of four baseline dynamic scans; each 40-slice set was collected at 25 time points for approximately 8 minutes of scanning. The spatial resolution of DCE-MRI sequence was compromised for a good temporal resolution with 10.6 s, and (D) diffusion weighted MR imaging. A total of 11 -values were used: 0, 10, 20, 30, 50, 80, 100, 120, 200, 700, and 800 s/mm2. Axial MR images with bilateral sacroiliac joints coverage were collected (Table 1).

2.3. Image Analysis
2.3.1. Analysis of IVIM Parameters

The DWI signal follows the biexponential model to calculate the signal attenuation IVIM, as [22] where and are the signal intensities in the pixels with different -values of and 0. Ds could be obtained by (2) using -values greater than 200 s/mm2 [23], and then, using all -values, the perfusion fraction of tissues () and the pseudo-diffusion coefficient (Df) were calculated by substituting Ds into (1).

All curve-fitting algorithms were performed using a home developed program based on MATLAB (MathWorks, Natick, MA) and developed by the Guangdong Provincial Key Laboratory of Medical Imaging Processing, School of Biological Engineering, Southern Medical University, which has been used to perform the IVIM data of liver [24]. The parametric maps of Ds, , and Df were abstracted by biexponential fitting.

2.3.2. Analysis of DCE MR Parameters

For DCE MR parameters analysis, all data were transferred to a computer equipped with manufacturer-supplied software (The Netherlands Easy Vision, release 4.4) and workstation (Philips Medical Systems, Best, The Netherlands).

For DCE-MRI, ME, RE, TTP, and BE were automatically calculated by manually drawing four different regions of interest (ROIs) from time-signal intensity curves (TSIC) in our study, as yet reported by Kambadakone and Sahani [25]. ME is the difference between the maximum () and the precontrast signal intensity (), which is related to the volume of the extracellular space. RE is the percentage of signal intensity increases between the and the , which reflects the plasma volume (Vp) and is highly correlated with microvessel density (MVD) and vascular endothelial growth factor (VEGF). TTP reflects the time from to within a region of interest, which is directly related to permeability surface area [26]. BE is the time between wash in and wash out of contrast agents [27].

2.3.3. Regions of Interest (ROI) Analysis

All regions of ROIs were positioned on DCE and DW images with -values of 0 by two radiologists, each with 10 years of experience in reading musculoskeletal system MR images with the double-blind method. The specific methods used were as follows.

Bone marrow edema (BME) of sacroiliitis in conventional MR images, DW images ( s/mm2), and DCE-MRI were used as references to determine lesion areas on IVIM and DCE-MRI parametric maps. Regions of interest (ROI) were manually positioned inside each lesion at its maximum transverse level, being along approximately 1 mm inward its contour. Meanwhile, the ROI was chosen to be as large as possible with minimum contaminations from surrounding unintended tissues, such as blood vessels, necrosis, and cystic area. However, if there were no BME in sacroiliitis, the ROI was manually placed on each side of the juxta-articular bone marrow, being adjacent to bone cortical or cartilage. At last, 4 ROIs were obtained with a size of ~2 mm2 on IVIM and DCE-MRI of each patient, and the average of the measurement was taken.

2.4. Statistical Analysis

All analyses were performed using GraphPad Prism 7.03 for Windows (GraphPad Prism, USA). Normality and homogeneity of variance were, respectively, performed by the Shapiro-Wilk test and the Levene test. IVIM and MRI-DCE parameters are reported as the mean ± standard deviation of all measurements. The three sets of measurements were obtained by two radiologists (radiologist 2 had 2 sets). In order to evaluate intra- and interobserver variability, the intraclass correlation coefficient (ICC) and coefficient of variation (CV) were calculated for the three sets of measurements. Intra- and interobserver agreement was evaluated using ICC and Bland-Altman analysis. The differences among the three sets of data and their average were assessed using one-factor Analysis of Variance (one-way ANOVA). Spearman’s rank correlation coefficient was calculated to measure the association between the IVIM and MRI-DCE parameters. Analyses of Spearman’s correlation coefficient between IVIM and DCE-MRI parameters were performed. For all tests, two-tailed P values < 0.05 were considered statistically significant.

3. Results

3.1. Patients

According to criteria of European Spondyloarthropathy Study Group in 1991 [28], a total of 44 patients finally were included (32 male, 12 female; mean age, years, range, 12–59 years; duration, years; CRP,  mm/L; ESR,  mm/H), and 19 patients were excluded because of unconfirmed AS, spondyloarthritis (SpA), osteoarthritis (OA), rheumatoid arthritis (RA), and scoliosis.

3.2. DCE and IVIM Parameters

The Shapiro-Wilk test showed that all DCE (RE, ME, TTP, and BE) and IVIM (Ds, , and Df) parameters had normal distributions (all ), and the Levene test revealed that all the variances were homogeneous (all ). Table 2 summarizes the mean and standard deviation values of the IVIM and DCE parameters. There was no significant difference among three observations of parameters (RE, ME, TTP, BE, Ds, , and Df), determined by observer 1 and 2 ( = 0.825–0.963). Table 3 presented the fact that the data of intra- and interclass coefficient correlation was 95% CI, ranging from 0.804 to 0.981. Figure 1 showed that there was a significant difference between monoexponential and biexponential curves of the IVIM diffusion in sacroiliac joints. Figure 2 illustrates the methods used to measure the parameters in the AS study in one example of the MRI-DCE images.

3.3. Intra- and Interobserver Reproducibility

Table 4 revealed that measurements of the DCE (RE, ME, TTP, and BE) and IVIM (Ds, , and Df) parameters were very reliable, possessing excellent intra- and interobserver reproducibility with the 95% CL ranging from 0.804 to 0.981. The coefficient of variation (CV) ranged from 2.6 to 11.07% (Table 3). Intra- and interobserver variations were found to be very low (Table 4), except for the ME value with −66.44–66.77 and −71.47–40.48%, intra-and interobserver variations, respectively.

3.4. Correlation between DCE and IVIM Parameters

The correlations between the DCE-MRI parameters (RE, ME, TTP, and BE) and IVIM (Ds, , and Df) parameters were evaluated using Spearman’s values and shown in Table 5. Moderate and positive correlations were, respectively, discovered between Ds and some of DCE parameters, including ME (, ) and RE (, ). The Ds to TTP (, ) and Ds to BE (, ) correlations were poor. The Df to ME correlation was moderate and negative (, ), while the Df to RE and TTP correlations were poor and negative (, ; , , respectively). However, there were no correlations found between and any parameters of ME, RE, TTP, and BE (Figure 3).

4. Discussion

In our study, signal-to-noise ratios (SNRs) in the raw EPI images were exceeded over 13.5, and 8 equally spaced -values under 200 s/mm2 were performed. Previous studies show that low SNRs affect the measurement at low values in IVIM [29], and that few low -values under 100 s/mm2 (0, 50, and 100 s/mm2) might decrease the accuracy of IVIM measurements [15, 20]. Therefore, the reliable Df and were obtained in sacroiliitis with AS.

In general, parameters of IVIM can be interpreted as separate intra- and extravascular compartments, with the parameter Ds referring to extravascular compartments, and with the parameters and Df referring to the intravascular compartment [20]. Moreover, Ds, , and Df are not influenced by each other [30, 31]. The perfusion parameters of IVIM pertained to different characteristics; namely, Df values reflect endovascular blood flow velocity and the mean capillary segment length, while values reflect the volume of vascular and extracellular space [30]. Both of them correlate with the amount of normal angiogenesis with intact vessels in terms of basement membrane thickness and pericyte coverage [32]. They could detect the increasing number of capillaries and blood vessels involved in angiogenesis and new bone formation in sacroiliitis with AS [11, 33].

Previous studies demonstrated that DCE-MRI necessarily involves not only exchange dynamics of intra- and extravascular but also other factors, such as blood vessel density, the pattern of blood delivery, vascular permeability, and distribution of contrast agent in lesions [34]. Parameters of DCE-MRI are related to microvascular volume or flow that can reflect some composite encompassing total tracer transit. Among parameters of DCE-MRI, ME is related to the volume of the extracellular space, RE reflects Vp and MVD, (3) TTP is directly related to the constant rate between extravascular space and blood Vp (permeability surface area) [26], which is generally affected by the injection and circulation time [7], and (4) BE also depends on various factors, like contrast dosage, MVD, blood flow, capillary wall permeability, composition of the extracellular space, and so on [8].

We found moderate and positive correlations between Ds in IVIM and ME and RE from DCE-MRI within sacroiliitis with AS. Our results were consistent with results reported by Bourillon et al. with the correlation between Ds and ME [35]. This might be explained by considering restricted diffusion within BME in AS because of inflammatory cells increase, the increase in the volume of extracellular space and MVD associated with angioedema [10, 35]. Poor and positive correlations between Ds in IVIM, and TTP and BE from DCE-MRI also were observed in sacroiliitis with AS. It has been proved that TTP increased due to the delayed enhancement within sacroiliitis in AS [6], and it might be accompanying higher BE.

We also found a moderate and negative correlation between Df and ME, which was not consistent with the result reported by Bane et al. in renal parenchyma with the poor and negative correlation between Df and ME [6]. This different result might be caused by the reliable Df obtained in our study. The correlation between Df and ME can imply that the amount of blood and serum of capillary transported to the marrow cavity was increased [11], and the volume of the extracellular space was decreased, in sacroiliitis with AS.

Poor and negative correlations between Df in IVIM and RE, TTP, and BE from DCE-MRI were obtained. François et al. have reported that the change of edema, myxoid, fibrosis, bone sclerosis, and the number of hematopoietic and fat cells contribute to comprehending of the MRI findings in sacroiliitis with AS [36]. These histopathologies of AS could explain partially the poor and negative correlations between the Df in IVIM and RE, TTP, and BE from DCE-MRI. Other influential factors might exist, like MVD, blood flow, capillary wall permeability, composition of the extracellular space, the injection time, contrast dosage, and so on [7, 8]. Hence, this poor correlation between Df, and TTP and BE might be meaningless.

On the contrary, there are no significant correlations between and the DCE-MRI parameters in our results. Zhao et al. have shown that had no statistical difference between active and inactive stage of AS [11]. Conversely, Bane et al. had reported statistical differences for the parameters of DCE-MRI between the active and inactive stage of sacroiliitis [6]. Both of results support our conclusion that was not correlated with the parameters based on DCE-MRI. Our results might demonstrate that IVIM provides the different perfusion information with DCE-MRI due to exogenous contrast agents and the different biophysical sensitivities, although both of them could reflect vascular volume and extracellular space [30]. Therefore, we speculated that might explain the recurrence of sacroiliitis without contrast-enhancing lesion. Further studies should be performed to reveal the reason why in sacroiliitis increases in the inactive stage of AS.

There were some limitations in our study. First, both DCE-MRI and DWI require longer scan times. The increased likelihood of patient movement increases susceptibility to motion artifacts and could introduce errors for the curve fit in the postprocessing phase biexponential mapping. Movement of the subjects was minimized by careful body fixation, and the images were coregistered in postprocessing. It is difficult to accurately match ROI of IVIM and DCE-MRI in juxta-articular bone marrow with sacroiliitis of AS, and partial volume effects might result in a little error of parameters between IVIM and DCE-MRI in the lesion. Second, in our study, due to the impact of SNRs and spatial resolution, small lesions (<5 mm) in BME were easily missed. Third, because only semiquantitative parameters were derived from DCE, the correlation between perfusion fraction and any DCE parameters should be studied through quantitative parameters, like mean transit time, arterial/blood vessel fraction, total flow, and so on. Finally, some heterogeneous flow preconditions, such as the rate of injection, impaired cardiac output, and proximal arterial stenosis, were not considered in parameters of DCE-MRI, which might affect perfusion parameters of DCE-MRI in lesions with AS, while parameters from IVIM could be measured repeatedly without concerning effects of exogenous contrast agents on renal functions. However, the SNR with 1.5 T and 4 channel body coils is relatively very low. Maybe a unilateral imaging will make the SNR higher in the future. The imaging time for DWI is quite long even with 6 mm slice thickness. A better resolution with maybe fewer values needs to be tested to make it more clinically useful. In the future, we will standardize MR imaging protocols and generation of parametric imaging, and we will perform model-based dynamic analytical parameters, and we will calculate mean transit time, arterial/blood vessel fraction, total flow, and so on. Whether correlates with the total flow/MTT is particularly interesting and important.

Accompanying AS treatment, a long-term observation on IVIM should be carried out, exploring the role of in the recurrence of active symptoms from an inactive stage of AS.

5. Conclusions

In this study we suggest that the IVIM perfusion imaging using 1.5-T MRI can be used in diagnosis and detection of the progression of AS. Ds in IVIM showed moderated positive correlations with ME and RE from DCE-MRI, and Df showed a moderate negative correlation with ME. However, was not related to any parameters of ME, RE, TTP, and BE. The IVIM MRI has be highlighted as a novel and effective technique to detect the progression of AS, because it can provide more perfusion information compared with DCE-MRI; besides the IVIM MRI is a noninvasive method.

Conflicts of Interest

The authors have no conflicts of interest to report.

Authors’ Contributions

Yinghua Zhao and Qun Zhang contributed equally to this article, and Shaolin Li and Zhiming Xiang contributed equally to this article too.

Acknowledgments

The authors thank Shujian Fan and Saeed Jerban for proofreading the manuscript. This work is supported in part by grants from the National Scientific Foundation of China (no. 81471810) and the Science and Technology Plan Projects of Guangdong Province, China (nos. 2014A020211018, 2014A020212399).