Application Progress of Intravoxel Incoherent Motion Diffusion Weighted Magnetic Resonance Imaging in Pancreatic Cancer ()
1. Introduction
Pancreatic cancer is known as “the king of cancers” and is the most malignant gastrointestinal tumor with occult onset, rapid progression and poor prognosis [1]. According to the data released by the American cancer society in 2019, the incidence of pancreatic cancer among male and female malignant tumors ranks 10th and 9th respectively, the mortality rate ranks 4th, and the 5-year survival rate is only 9% [2] [3]. Early symptoms of pancreatic cancer are not typical. Once the cancer is found to be late, it is very likely, so early diagnosis is of great significance for the prognosis of patients with pancreatic cancer. Diffusion weighted imaging (DWI), which was first used in the diagnosis of brain diseases, is the only non-invasive imaging technique that can reflect the diffusion of water molecules in vivo without the need of contrast medium injection and quantify water molecules diffusion and the pathological characteristics of tumors by measuring apparent diffusion coefficient (ADC) [4]. The traditional DWI imaging technology is based on the single exponential model, and the ADC value is the combined effect of diffusion of water molecules in tissues and microcirculation perfusion, which is higher than the diffusion coefficient of pure water molecules [5] [6]. It cannot really reflect the characteristics of the microstructure of the organization. Le Bitan et al. [7] put forward the intravoxel incoherent motion imaging technology in 1986, which is based on a dual exponential model with multiple b-values for analysis, enabling to separately reflect pure tissue diffusion and information about microcirculation perfusion. IVIM-DWI can simultaneously provide quantitative analysis to more accurately assess the movement of water molecules in the tissue. In recent years, intravoxel incoherent motion has been widely used in the pancreas. This article will review the principles of IVIM imaging and its application in pancreatic diseases.
2. Basics of IVIM DWI Imaging
The traditional DWI uses the single exponential model whose formula is
(1)
to calculate ADC value, where b value is the diffusion sensitivity coefficient in the unit of s/mm2, reflecting the sensitivity degree of water molecule dispersion, the ADC is the apparent diffusion coefficient whose unit is mm2/s. S0 refers to the signal strength in voxel on the b = 0, Sb refers to the signal strength in voxel at a specific b value. This formula assumes that DWI signals show a single exponential attenuation, and at least two scans with different b values are required for linear fitting to obtain the regression slope. However, the tissue contains not only the water molecule diffusion effect but also the capillary microcirculation perfusion effect. In order to distinguish the two, the biexponential IVIM model is used to describe the relationship between the attenuation of signals in the tissue and the b-value. According to the equation proposed by Le Bitan [8]
(2)
where Sb refers to the signal intensity in voxel at a specific b value, S0 refers to the signal intensity in voxel on the b = 0, f is the perfusion fraction (%), which represents the proportion of the diffusion effect of voxel microcirculation perfusion in the overall diffusion effect, D* (mm2/s) is a false diffusion coefficient, which reflects the blood perfusion diffusion effect; D (mm2/s) is a true diffusion coefficient, which indicates the pure water molecule diffusion effect. Since D* is significantly higher than D [5], when b takes a highb value (>200 mm2/s), the microcirculation perfusion is not sensitive to the signal intensity. The formula can be simplified as
(3)
the perfusion coefficient D of pure water molecules can be obtained and then the values of D* and f can be obtained through the nonlinear regression algorithm by substituting D into Equation (2). Conversely, when b takes a low value, the attenuation of tissue signals is mainly caused by microcirculation perfusion. The closer b is to 0, the more it can be considered that the attenuation of tissue signals is caused by microcirculation perfusion [9].
3. Clinical Application of IVIM in Pancreatic Diseases
3.1. Application of IVIM in Pancreatic Cancer
Pancreatic ductal adenocarcinoma has the characteristics of late detection, easy recurrence, and resistance to chemotherapeutic drugs, with a high mortality rate. Yoon et al. [10] found that with the increase of pancreatic fibrosis, IVIM (ADC,D, f, and D*) related parameters all showed a downward trend, but only D* and f have significant difference between different stages of fibrosis. Hecht et al. [11] used IVIM-DWI related parameters to study the histopathology of pancreatic ductal adenocarcinoma (PDAC) and found that D value was negatively correlated with matrix fibrosis (pooled r = −0.46, p = 0.04) and positively correlated with f (r = 0.44, p = 0.05) suggesting that both perfusion and diffusion could lead to interstitial fibrosis hyperplasia. This is consistent with the fact that PDAC is a tumor with deficient blood supply, it is possible that pancreatic cancer cells are rich in fibrous matrix, which limits their diffusion. ADC values (derived by single index and double index fitting) were significantly reduced in tumors with dense fibrosis and could be used as a biomarker for fibrosis hyperplasia of pancreatic ductal adenocarcinoma, which was similar to the results of previous studies [10] [12] [13]. However, the variation trend of f value in Hecht et al.’s study is different from that of Yoon et al. [10]. It is possible that the study of Hecht et al. selected non-tumor tissue adjacent to the lesion rather than normal healthy pancreatic tissue. In the study of Zhang [14] et al., it was believed that diffusion parameters such as ADC and D value were valuable for the evaluation of pancreatic fibrosis stage, and D value decreased with the increase of fibrosis degree, which was consistent with the relationship between D value and fibrosis in the previous two studies. The FDA believes that gemcitabine can significantly improve the quality of life of patients with pancreatic cancer and prolong the life of patients. It can be used as the gold standard for the treatment of pancreatic cancer [15]. The mechanism is that gemcitabine can inhibit the growth of cancer cells and promote cell apoptosis. Studies [16] [17] have shown that gemcitabine has different efficacy in the treatment of pancreatic cancer, and drug resistance may occur in some patients. Therefore, early evaluation of chemotherapy effect is of great significance for the treatment of pancreatic cancer. Li et al. [18] established 24 nude mice model of pancreatic cancer and divided them into gemcitabine group (n = 12) and normal saline group (n = 12). Three b values of DWI (b = 50, 400 and 800 s/mm2) and multiple b values of IVIM-DWI (b = 0, 25, 50, 80, 100, 300, 500, 800 s/mm) were used to compare ADC, D, f and D* values of the two groups. The results showed that compared with normal saline group, ADC and D values in the tumor area of the guitar gemcitabine group were significantly increased. There was no significant difference between D* and f values, suggesting that ADC values and D values could be early imaging indicators of chemotherapy response of pancreatic cancer.
3.2. Application of IVIM in the Prediction of Pathological Grade of Pancreatic Ductal Adenocarcinoma and the Judgment of Lymph Node Metastasis
ADC derived from DWI can be used to evaluate the microstructure of early disease tissues and is widely used in the diagnosis of pancreatic diseases (PC). However, there have been different reports on the role of ADC values in the diagnosis of pancreatic cancer. Report [19] indicated that the ADC value of PC was higher than that of the normal pancreas group, while some studies showed the ADC value of PC is much lower than the normal pancreas [20] [21]. The ADC value includes not only pure water molecular diffusion but also microcirculation effects, making the results inaccurate in these studies, in addition the ADC value calculated by traditional DWI is greatly affected by b value, so different b value selection will also affect the final result. Ma et al. [22] used IVIM technology to predict the grading of pancreatic ductal adenocarcinoma with 15b values (0, 10, 20, 40, 60, 80, 100, 150, 200, 400, 800, 1000, 1200, 1500, and 2000 s/mm2). The ADC values, D values, f values, and D* values were measured for 38 pathologically confirmed PDAC. The results showed that the mean D value of well/moderately differentiated PDAC was lower than that of poorly differentiated PDAC and there was a statistical difference. However, the mean f value of well/moderately differentiated PDAC was higher than that of poorly differentiated PDAC, and the f value was higher than that of the D value (0.894 > 0.865) when differentiating the area under the subject curve of well/moderately differentiated pancreatic cancer from poorly differentiated pancreatic cancer. Matrix reaction and fibrosis are the main histopathological features in pancreatic cancer, well/moderately differentiated PDAC shows abundant fibrosis and glandular formation characterized by neoplastic tubular and duct-like structures which contain massive mucus rich in macromolecular protein restraining diffusion of water molecules [23]. Thef value represents the perfusion information of the tissue and has shown significant positive correlation with microvascular density (MVD) [24]. Study has shown that well/moderately differentiated PDAC contains higher MVD than poorly differentiated group, which may explain the higher f value of the former than the latter. It is concluded that IVIM technique can be used as a potential non-invasive tool to predict PDAC pathological grading to some extent.
In the study of PDAC lymph node metastasis, Rong et al. [25] found that IVIM parameters (ADC, D value, f value and D* value) were correlated with lymph node metastasis. This study used ROC curve analysis on 59 lymph nodes of 15 patients with pancreatic cancer. The IVIM parameters of metastatic and non-metastatic lymph nodes were measured. The average value of each parameter in the metastatic group was significantly lower than that in the non-metastatic group. The D value was the highest in identifying the metastasis of lymph nodes, followed by ADC (0.940), D* (0.867) andf (0.855), the AUC increased to 0.997 when the parameters D and f were combined, suggesting that the combination of IVIM parameters may be more helpful in diagnosing lymph node metastasis. So far, there have been few studies on the pathological grade and lymph node metastasis of pancreatic cancer with IVIM technology, so it is necessary to expand the sample size to further study the clinical significance of IVIM derived parameters.
3.3. Application of IVIM in Differentiating Pathological Classification of Pancreatic Cancer
Pancreatic ductal adenocarcinoma (PDAC) and pancreatic neuroendocrine tumor (pNET) are the first and second most common cases of pancreatic cancer [26], and the treatment and prognosis of the two tumors are completely different. Compared with PDAC, pNET has better response to chemotherapy, higher resectability and longer overall survival. In view of the high resectability and relatively good prognosis of pNET, it is significant to distinguish pNET from PDAC before surgical treatment [27].
By retrospectively analyzing the clinical data of 18 pNET patients and 32 PAC patients, Wang et al. [28] discussed the value of texture analysis of IVIM parameters in differentiating pancreatic neuroendocrine tumor (pNET) from pancreatic adenocarcinoma (PAC). The study showed that the mean f value of pNET group was significantly higher than that of PAC group (27.0% vs 19.0%, P = 0.001), while the mean D* value and mean D value had no significant difference between the two groups. The author believed that f value is beneficial to distinguish PNET and PAC, which is similar to previous research results [29] [30]. It also shows that multi-texture analysis of IVIM parameters can provide more diagnostic information than IVIM parameters and can be used as an effective tool to distinguish pNET and PDAC. However, because of the small number of subjects and the complexity of tumor analysis, it is difficult to perform in clinical practice. Ma et al. [31] studied 37 cases of PDAC and 17 cases of high-grade pNET by using IVIM-DWI model. TheD* value of pNET was significantly higher than PDAC, which was the same with the results of Kang et al. [32]. The D value is significantly lower than PDAC, which is consistent with the conclusion of Ma et al. [33]. However, it is different from Kang et al. [32] ’s conclusion thatD value has no effect in distinguishing PDAC and PNET, which may be related to the number of b value selection. Kang et al. [32] used 9b values in the study, while 15b values were used in the Ma’s study. The more b is selected, the better the image quality is, which may affect the measurement results. The combination of D value and D* value can achieve the highest diagnostic efficiency, and the specificity and sensitivity for the identification of high-level pNET and PDAC are 76.9% and 100%, respectively. Taken together, these results indicate that IVIM-DWI is likely to be a very valuable biomarker for the identification of pancreatic tumors. Kim et al. [34] performed pancreatic ductal adenocarcinoma (n = 60), neuroendocrine tumors (n = 15), and solid pseudopapillary tumors (n = 9), 20 cases of pancreatitis [13 cases of acute pancreatitis, 7 patients with AIP (autoimmune pancreatitis)] and 30 healthy volunteers (normal pancreas) underwent 10 b-value IVIW-DWI examination, and the f-values and D* values of PDAC, ACP, AIP, and SPT were significantly reduced compared with normal pancreas which is similar to previous research [35]. In terms of pNET, f-values are higher than PDAC, SPT, ACP, and AIP, but usually only pay attention to the statistical differences between pNET and PDAC. The f-value in IVIM parameters has the best diagnostic performance among all the parameters that distinguish PDAC from NP. Compared with pNET, the decrease in f-value of PDCA is consistent with the results of Klau et al. [30]. Since PDAC is pathologically a tumor with insufficient blood supply and rich in fiber matrix, and there is a good correlation between perfusion related parameterf value and microvascular density (MVD), thef value and MVD of PDAC group are significantly lower than that of pancreatic neuroendocrine tumor group. Study [36] showed that thef value of AIP and PDAC group was significantly lower than that of the normal control group, while the f value of PDAC group was lower than that of AIP group, which was consistent with the results of De Robertis et al. [37]. The above research results show that IVIM-DWI perfusion-related parameters (D*, f-values) are helpful in determining the degree of internal blood supply of the tumor and in the diagnosis and differential diagnosis of pancreatic neoplastic diseases. Some recent research results related to IVIM are summarized in Table 1.
4. Advantages and Disadvantages of IVIM
IVIM-DWI can obtain and distinguish the perfusion information of water molecules and microcirculation by using the dual-exponential model with multiple
Table 1. Different researches on IVIM derived parameters and its diagnostic performance.
a: PDAC = pancreatic ductal adenocarcinoma; PAC = pancreatic adenocarcinoma; NP = normal pancreatic; pNET = pancreatic neuroendocrine neoplasms.
different b values, to provide more accurate and comprehensive diagnostic information. Various studies have shown that the quantitative analysis of IVIM parameters can not only reflect the histopathological characteristics of the tissue, but also play an important role in the classification of pancreatic tumors. Moreover, there is no need to use contrast agent in IVIM-DWI, which has an absolute advantage for patients with renal insufficiency. However, there is no unified international standard for the size and quantity of b value, and the difference of b value may lead to slight differences in different research results, which may affect the repeatability of parameters. The image signal-to-noise ratio (SNR) will decrease while the selected b value is higher than 1000 s/mm2, which directly affects the image quality. Poor quality of image would affect the accuracy of measurement. However, with a smaller b value, the image with a high SNR can be obtained, but it is not sensitive to the detection of water molecule diffusion motion, and it is easy to miss the lesions. The more and higher b values would increase scanning time, causing motion artifacts and affecting the post-processing of the parameters. Therefore, more in-depth clinical studies are needed to explore and improve the optimal size and quantity of b value, the optimal fitting model and the accurate post-processing measurement, so as to obtain a unified scanning parameter scheme, which can better reflect the technical value of IVIM-DWI and obtain a more comprehensive understanding of the disease.