ISSN: 1899-0967
Polish Journal of Radiology
Established by prof. Zygmunt Grudziński in 1926 Sun
Current issue Archive Manuscripts accepted About the journal Editorial board Abstracting and indexing Contact Instructions for authors Ethical standards and procedures
Editorial System
Submit your Manuscript
SCImago Journal & Country Rank
1/2020
vol. 85
 
Share:
Share:
Technology and contrast media
Review paper

Intravoxel incoherent motion magnetic resonance imaging: basic principles and clinical applications

Aleksandra E. Szubert-Franczak
1
,
Martyna Naduk-Ostrowska
1
,
Katarzyna Pasicz
2
,
Joanna Podgórska
1, 3
,
Witold Skrzyński
2
,
Andrzej Cieszanowski
1

1.
Department of Radiology I, Maria Sklodowska-Curie National Research Institute of Oncology, Warsaw, Poland
2.
Department of Medical Physics, Maria Sklodowska-Curie National Research Institute of Oncology, Warsaw, Poland
3.
2nd Department of Clinical Radiology, Medical University of Warsaw, Warsaw, Poland
Pol J Radiol 2020; 85: e624-e635
Online publish date: 2020/11/18
Article file
- Intravoxel.pdf  [0.21 MB]
Get citation
 
PlumX metrics:
 

Introduction

The intravoxel incoherent motion (IVIM) diffusion-weighted (DW) model as a possible imaging technique, using multiple b values and bi-exponential fitting for the concurrent estimation of the pure molecular water diffusion and microcirculation of blood water in randomly oriented capillaries (perfusion) was first introduced in the late 1980s by Le Bihan et al. [1]. The idea to use diffusion and IVIM magnetic resonance imaging (MRI) to acquire perfusion parameter maps was considered revolutionary but technically difficult, and, as a result, it was more than 20 years before the method started being used in clinical practice [2]. IVIM reflects the random microscopic motion of water molecules that occurs in each voxel on MR images not only in intra- or extracellular space but also in microcirculation of blood [3]. According to IVIM theory, diffusion and perfusion are affected by several tissue characteristics, including the presence of restrictive barriers within tissue, the viscosity of the fluid in which the spins are diffusing, and the velocity and fractional volume of perfusing spins [4]. Formerly, due to degradation of images caused by cardiac, respiratory, and other motion artifacts, IVIM imaging was restricted to neuroradiologic applications. Nowadays, it is pursued to apply IVIM MRI to evaluate almost entire human body. Over the last few years there has been a revival of interest in IVIM MRI and its applications in many fields, particularly in oncology [2].

Basic principles and acquisition of intravoxel incoherent motion perfusion magnetic resonance imaging

DW magnetic resonance imaging is based on the random movement of individual water protons, called Brownian motion. The Brownian motion can be described by a series of molecular jumps [5].

According to Einstein’s equation the diffusion coefficient D depends on the averages of individual displacements between two collisions and the mean molecular velocity. In the presence of a magnetic field, random displacement of individual molecules between two gradient pulses results in signal attenuation. Initially, the signal attenuation increase was described by a monoexponential function:

ss0=exp(-b×D)

where: S is the measured signal intensity, S0 is the signal intensity without the influence of diffusion, D (mm2/s) is the diffusion coefficient of water, and the b value (s/mm2) can be calculated as follows: γ – gyromagnetic ratio (MHz/T), G – amplitude of the two diffusion gradient pulses (mT/m), δ – duration of the pulses (ms), Δ – time between the two pulses (ms) [1].

The assumptions of the IVIM model, described for the first time by Le Bihan et al. in the late 1980s, are based on the translation movements at voxel levels. Microcirculation of the blood in the capillary network would mimic a pseudo-diffusion process [1]. The IVIM effect appears at low b values (i.e. b ≤ 200 s/mm2) as a deviation of tissue diffusion signal decay (Figure 1). Because the pseudo-diffusion coefficient (D*) is one order of magnitude higher than the diffusion coefficient D, the exponential decay with the pseudo-diffusion coefficient disappears faster. At low b-values, the perfusion effect predominantly contributes to the overall diffusion signal. At higher values of b, the exponential input with D modifies the signal. MRI signal attenuation is the sum of the tissue and blood component, taking the shape of biexponential decay:

Figure 1

The intravoxel incoherent motion effect

/f/fulltexts/PJR/42604/PJR-85-42604-g001_min.jpg
ss0=(1f)×(b×D)+f×exp(b×(D*+D))

where: S – the signal intensity in the pixel with diffusion gradient b, S0 – the signal intensity in the pixel without diffusion gradient, D – the true diffusion as reflected by pure molecular diffusion (mm2/s), f – the fractional perfusion related to microcirculation, and D* – pseudo-diffusion coefficient representing perfusion related diffusion or incoherent microcirculation (mm2/s). D* is sometimes referred to as ADCfast, while D may be also called ADCslow.

The IVIM sequence consists of several scans: the first without the use of a coding gradient, the next with different values of the amplitude and duration of the gradient (different b-values). An important part of the IVIM acquisition protocol is the choice of b-values (in terms of their number and distribution), resulting in a different degree of diffusion weighting in the acquired images. The b-values used in IVIM imaging are usually in the range of 0-900 s/mm2. The number of b-values ranges from 4 to 16. Number of b-values in the range of 0-200 s/mm2, where the perfusion effect is dominant, is between 2 and 12 (Table 1).

Table 1

The b-value used in clinical trials

PaperYearb-values (s/mm2)Number of b-values
Guo [9]20160, 10, 20, 30, 50, 70, 100, 150, 200, 400, 800, 100012
Liu [10]20130, 10, 20, 30, 50, 70, 100, 150, 200, 400, 800, 100012
Ichikawa [11]20130, 10, 20, 30, 40, 50, 80, 100, 200, 400, 800, 10008
Rheinheimer [12]2012[0; 50], [0; 100], [0; 150], [0; 200], [0; 300], [0; 400], [0; 600]; [0; 800]8
Cui [13]201516 b-values from 0 to 80016
Dyvorne [14]20130, 15, 30, 45, 60, 75, 90, 105, 120, 135, 150, 175, 200, 400, 600, 80016
Federau [15]20140, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 90015
Federau [16]20140, 10, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 90016
Döpfert [17]20110, 50, 500, 800
Kakite [18]20160, 15, 30, 45, 60, 75, 90, 105, 120, 135, 150, 175, 200, 400, 600, 80016
Shim [19]20150, 10, 20, 40, 60, 80, 100, 120, 140, 160, 180, 200, 300, 500, 700, 90016
Conklin [20]20160, 10, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 90016
Federau [21]20160, 10, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 90016
Xu [22]20160, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 400, 80014
Boss [23]20160, 10, 30, 60, 120, 200, 300, 430, 600, 80010
Valerio [24]20160, 10, 20, 30, 40, 50, 80, 100, 200, 400, 80011
Pang [25]20130, 188, 375, 563, 7505
Bane [26]20160, 15, 30, 45, 60, 75, 90, 105, 120, 135, 150, 175, 200, 400, 600, 80016

There are, however, several important limitations that have so far prevented the use of IVIM in routine practice. The most important is the lack of a uniform data analysis method, which is particularly important because IVIM imaging is based on the quantitative analysis of f,D, and D* parameters [1,5,6].

Le Bihan proposed acquisition with only three values of b that are theoretically sufficient to obtain IVIM results. However, more points are needed, especially for the brain due to “noise contamination” [5]. For accurately estimating IVIM parameters at least 10 of the optimal distributed b-values should be used [7].

In addition, IVIM parameters, similarly to ADC values, are possibly field-strength dependent (1.5T vs. 3.0T) [8].

Strengths and limitations of intravoxel incoherent motion magnetic resonance imaging

Understanding of the strengths and limitations of IVIM concept may help to appreciate a clinical benefit of this method (Table 2).

Table 2

Advantages and disadvantages of intravoxel incoherent motion (IVIM) magnetic resonance imaging

AdvantagesDisadvantages
Non-invasive techniqueCardiac and respiratory motion artifacts
No use of contrast agentsBulk flow phenomenon artifacts
No need of radioisotopes injectionNo standardization of acquisition protocols and models
No ionizing radiationVariance in calculated parameters among studies
Differential sensitivity to vessel size
Lack of standardization of calculation of IVIM parametres
No established values for normal organs

The most important advantage of IVIM diffusion MR imaging is that as a non-contrast perfusion imaging modality it can be used in situations in which intravenous administration of contrast agents is not clinically justified, and it may serve as an interesting alternative to contrast-enhanced perfusion MR imaging in some patients with contraindications to contrast agents, such as severely compromised renal function. Furthermore, IVIM diffusion MR imaging does not involve ionising radiation or injection of radioisotopes.

The most important disadvantages include artifacts related to cardiac and respiratory motion, artifacts from other bulk flow phenomena, no standardisation of acquisition parameters, and various algorithms used for quantitative image analysis. Moreover, the sensitivity of IVIM MRI not only varies according to vessel size but also depends on the utilised number and distribution of applied b values. Thus, IVIM-derived parameters and maps are highly variable and inconsistent. Due to the lack of standardisation of the IVIM technique, significant variance in calculated parameters among studies has been observed and no values for normal organs have been established.

Clinical applications in oncology

Abdominal and pelvic organs

Liver

IVIM was applied to study liver diseases such as liver fibrosis, nonalcoholic liver fatty disease, and focal liver lesions and to estimate possible treatment response [3,30]. With IVIM DWI, Luciani et al. suggested that the diffusion component related to the molecular displacement (D) does not differ significantly between cirrhotic and healthy livers [27].

Yamada et al. found that D was significantly lower than ADC, thus suggesting that differences reported in ADC between patients with cirrhosis and healthy patients were mainly related to the perfusion component of liver diffusion. The fact that the pure molecular diffusion coefficient, D, was similar in both the healthy liver group and the cirrhotic liver group further indicates that changes in liver architecture may be of less importance than changes in liver perfusion [28] (Figure 2).

Figure 2

Sample quantitative maps of intravoxel incoherent motion (IVIM) parameters (D and f, respectively) created using own plugin for ImageJ software for breast cancer liver metastases (arrows). The increase in blood flow/volume occurring in liver metastases can be detected and mapped with IVIM magnetic resonance imaging

/f/fulltexts/PJR/42604/PJR-85-42604-g002_min.jpg

However, the large overlap of D,f, and D* for different stages of fibrosis suggests that this technique cannot be used to accurately characterise stages of liver fibrosis [27].

In a large number of studies the IVIM technique was used for the assessment of focal liver lesions. Liver nodules including hepatocellular carcinoma (HCC), haemangioma, focal nodular hyperplasias (FNH), and liver metastasis (MET) showed overlap in terms of IVIM parameters. However, other researchers, like Yamada et al., showed the potential of IVIM MRI to differentiate HCC from haemangiomas and cysts [28].

Summarising, IVIM MRI for assessing hepatic lesions imposes some difficulties, mainly due to the possible overlap between normal and diseased liver tissue, as well as between benign and malignant lesions [10].

Pancreas

IVIM imaging has recently gained interest as a method to properly recognise and characterise pancreatic lesions.

Perfusion-related IVIM-derived parameters seem especially auspicious for differentiation between normal pancreatic parenchyma and solid pancreatic lesions with high sensitivity and specificity [29].

According to Re et al. there is a compelling difference in the average IVIM perfusion fraction of pancreatic adenocarcinoma compared to healthy pancreatic tissue [31]. IVIM MRI emerges as a promising method for characterisation of pancreatic lesions and should therefore be further investigated.

Kidney

Over the past two decades, the number of properly diagnosed renal tumours has escalated due to increased utilisation of different imaging techniques. Unfortunately, preoperative classification of solid renal masses remains challenging and defective. In light of the above, in the past few years there have been attempts to use IVIM parameters to distinguish renal tumour subtypes. Chandarana et al. [32] established the differentiation of subtypes of renal cell carcinomas using perfusion fraction (f) and tissue diffusi-vity (D*). Interestingly, the authors demonstrated a significantly higher f and lower D* in enhancing renal masses compared with markedly and poorly enhancing renal masses, differentiating these two groups with a higher accuracy than standard ADC values [33]. The combination of f and D* allowed the diagnosis of papillary renal cell carcinoma (RCC) and cystic RCC with 100% accuracy, and clear cell RCC and chromophobe RCC with 86.5% accuracy. The f was shown to have a good correlation under the curve of gadolinium concentration at 60 seconds, which is a measure of the gadolinium concentration obtained from the area under the gadolinium-enhancement signal intensity (SI) curve over the first 60 seconds (CIAUC60) after intravenous injection of contrast, which potentially enable the assessment of tumour vascularity without the need for intravenous administration of gadolinium contrast agent [32,34].

According to Ding et al. [35] IVIM-derived parameters, especially the perfusion-related parameters (D* and f), showed greater diagnostic accuracy than that of ADC values in differentiating non-ccRCCs from lipid-poor AMLs.

Prostate

Recent studies have reflected an interest in multiparametric MR imaging including diffusion MRI in the detection, staging, and post-treatment follow-up of prostate cancer [36].

The results of published papers utilising the IVIM model are inconsistent in terms of prostate cancer detection, staging, and treatment response. According to some studies, the diagnostic value of this technique compared to ADC is rather indigent [17] whereas others report showed significant differences between the IVIM parameters of benign and malignant tissue [24].

Concluding, there is a visible fluctuation between IVIM parameters in cancer and normal prostatic within the literature [37], and there is still a need for validation in the application of IVIM in the diagnosis of prostate cancer.

Head and neck

Over the last two decades IVIM imaging became a very promising technique for differentiation of tumours of the head and neck area, due to distinctive f and D.

IVIM imaging has been successfully applied for differentiation not only between benign and malignant salivary gland lesions [38] but also for discrimination of other types of tumours in this region. The perfusion-related parameter values have been significantly different among different types of head and neck tumours, including squamous cell carcinoma, benign (pleomorphic adenomas, Warthin tumour) and malignant salivary gland tumours, lymphomas, and schwannomas [39].

Additionally, low pre-treatment D and f values favour good response to treatment, whereas an increase in D during treatment confirms good response to therapy [40]. However, the role of IVIM MRI in characterisation of lymph nodes remains controversial because IVIM parameters have been shown to be inappropriate for differentiation between metastatic and non-metastatic head and neck lymph nodes [41]. On the other hand, some studies indicate that IVIM DWI is feasible in the diagnosis of lymph node metastasis, which demonstrate decreased D (due to increased nuclear-to-cytoplasmic ratio) and increased D* (due to new blood vessel formation and increased parenchymal perfusion) [42].

Furthermore, pre-treatment and mid-treatment IVIM- DWI led to the conclusion that there is a potential for prediction of chemo-radiotherapy response of cervical lymph nodes in head and neck squamous cell carcinoma patients [43].

What is more, in patients with head and neck squamous cell carcinoma (HNSCC) due to locoregional failure, high initial f in lymph nodes may be considered as an indicator of poor treatment response [44].

IVIM MRI (perfusion-related coefficients) appears promising for noninvasive evaluation of early changes of irradiated major salivary glands during radiotherapy [45]. Moreover, it is reasonable and advantageous for predicting and assessing initial responses of human papilloma virus and oropharyngeal squamous cell carcinoma SCC to chemoradiotherapy [46].

Reported artifacts, specific for this area and related to the presence of bones and air-filled structures, are additional problem in the evaluation of the head and neck with MRI. Moreover, specific artifacts related to jaw movements, swallowing, speaking, coughing, and respiration might additionally degrade acquired images.

Breast

In recent years, with the improved MR hardware, a renewed interest in IVIM has been shown and the applications of IVIM to breast lesions have been reported [47].

According to some authors, IVIM has been recognised as a very promising method for the diagnosis of breast cancer [48,51,52].

Compared to DWI with monoexponential fit, IVIM provides separate quantitative measurement for cellularity and vascularity. IVIM can be used to differentiate benign and malignant lesions with high specificity. This method can increase the diagnostic sensitivity and might play a role in screening breast MRI in high-risk women [47].

What is more, quantification of perfusion fraction, tissue diffusivity, and pseudodiffusivity provides noninvasive sensitivity to microenvironment properties without need of contrast agent. These applications have the potential to improve the specificity of breast MRI.

According to Sigmund et al., further work in a larger patient population is needed to validate the ability of IVIM for the diagnosis, differentiation, monitoring, and management of cancerous lesions in the breast and verify its role as a possible surrogate marker of the biological properties of the tumour [49].

Other lymph nodes

Some recent studies have shown that particular IVIM parameters, such as the ADC and diffusion coefficient (D), are significantly different between patients with metastatic and non-metastatic lung cancer that has spread to the lymph nodes, therefore possibly facilitating discrimination of benign and malignant mediastinal lymph nodes [53].

Furthermore, an IVIM sequence may also be helpful in diagnosing metastatic lymph nodes of rectal carcinoma. According to Qui, average D and ADC values are more sensitive than f and D* values for this purpose [54]. As reported by Yu et al. in patients with rectal cancer, metastatic mesorectal lymph nodes exhibit lower D and D* values, compared to non-metastatic lymph nodes. The authors concluded that IVIM DWI may be helpful in identifying mesorectal nodal involvement in rectal cancer patients [55].

Clinical applications in neuroradiology

There is little doubt that IVIM imaging has potential for evaluation of brain ischaemia. Wirestam et al., who utilised IVIM for the assessment of the stroke by IVIM, reported that the perfusion fraction was reduced in areas affected by ischaemia in comparison to the respective contralateral region [56].

In acute stroke, a significant decrease in f was measured in the infarct core compared with the contralateral hemisphere in three independent studies [57-59].

Furthermore, it has been found that decreased IVIM perfusion fraction f and blood flow-related parameter D* might correspond with proximal artery vasospasm development after cerebral aneurysm rupture and delayed cerebral ischaemia [60].

What is more, IVIM may be a valid and favourable method for significant measurement of brain perfusion without intravenous administration of contrast material. According to some studies, IVIM perfusion parameters are highly reactive to hyperoxygenation-induced vasoconstriction and hypercapnia-induced vasodilatation [61].

Another study showed that IVIM can be used for the assessment of cerebral small vessel disease and its complications including lacunar stroke, leukoaraiosis, or even vascular dementia [62].

Brain tumours

The most common malignant brain tumours are gliomas, which originate from glial cells (astrocytes, oligodendrocytes and ependymocytes). Due to the high mortality rate, early diagnosis and correct assessment of the grade of a glioma are crucial for the planning of a therapeutic strategy and establishing a prognosis [62,63].

Researchers reported that both D* and f could be used to differentiate gliomas because these parameters were higher in high-grade gliomas (HGG; WHO: III-IV) in comparison to low-grade gliomas (LGG; WHO: I-II) [16,64]. In subsequent reports the authors concluded that IVIM-derived metrics characterising both perfusion (D*, f) and diffusion (D) have been assessed as promising imaging biomarkers in preoperative differentiation of glioma grades. The f times D* value demonstrated the best diagnostic performance in grading gliomas compared with other parameters and was the only parameter showing a significant difference between grade III and grade IV [65]. Additionally, according to some scientists IVIM metrics could be a potential biomarker for survival in patients with glioblastomas as it has been found out that f and D* measured in contrast-enhancing regions correlate well with response to therapy and survival [66].

It is also worth noticing that there are studies showing that IVIM may be used to differentiate between glioblastoma and primary central nervous system lymphoma (PCNSL) [19,67,68]. According to authors mean maximum f is significantly higher in the glioblastoma group than in the atypical primary central nervous system lymphoma that often mimics glioblastoma [67,68].

The most recent studies showed that the f-value can be used as a noninvasive quantitative imaging measure to directly assess the vascular volume fraction in brain tumours. Moreover, one study showed an obvious correlation with the histological vascular density in meningiomas [69].

Moreover, IVIM parameters have been used to monitor the tumours treated with antiangiogenic and vascular target agents that have been recently considered an alternative or complementary therapy to conventional cancer treatments [45,62,70-73].

According to some authors, IVIM may be used to differentiate radiation necrosis from tumour progression in brain metastases that are treated with stereotactic radiotherapy (radiosurgery). The f-value was found to be low in the region of the radiation necrosis, whereas tumour recurrence was more heterogeneous and showed higher values [74].

Other clinical applications

IVIM MR imaging might also be used as a diagnostic tool in gynaecology, orthopaedics, dermatology, and ophthalmology. Presentation of these applications is beyond the scope of this paper; however, some examples of various clinical applications according to the PubMed database are presented in Table 3.

Table 3

Clinical applications of intravoxel incoherent motion magnetic resonance imaging

NeuroradiologyOncologyOther
Cerebral infarctionBrain tumorsRenal function, e.g. renal fibrosis, early changes in diabetic kidneys
Brain perfusionPrimary and non-metastatic head and neck tumorsCardiac imaging
Pulsatility of the brain’s microvasculatureLymph nodal metastasisDermatomyositis imaging
Pancreatic lesionsHuman Achilles tendon imaging
Prostate lesionsLactating breast lesion imaging
Breast tumorsOrbital lesions imaging
Liver lesionsVertebral bone marrow imaging
Renal lesionsThyroid benign and malignant nodule differentiation
Cervical cancerTissue characterization of the uterine fibroids
Soft tissue tumorsIschemic optic neuropathy
Diagnosis of pediatric solid abdominal tumorsTemporomandibular joint disorder imaging
Monitoring treatment efficacy of chemo- or radiotherapy, effectiveness of antiangiogenic drugs and vascular targeting agentsKnee joint in children with juvenile idiopathic arthritis imaging
Invasive fungal infection in the lung treatment response imaging
Parotid gland imaging in patients with Sjögren diasease
Placental perfusion imaging

Conclusions

The IVIM method is a DW MRI technique incorporating low b-values (< 200 s/mm2) for the estimation of perfusion parameters, which offers many advantages in comparison to the conventionally used DWI technique. It is a noninvasive substitute for the measurement of perfusion that does not require intravenous injection of exogenous contrast agents, such as gadolinium, and can be a great source of information on perfusion and diffusion simultaneously [62].

Today, the most promising clinical application for IVIM perfusion MRI is oncology. The graphs demonstrate the number of articles on implementation of IVIM MRI in oncology (Figure 3), compared to the total number of articles about this technique published between 1986 and 2018 (Figure 4).

Figure 3

Number of articles about intravoxel incoherent motion magnetic resonance imaging 1987-2018 (PubMed)

/f/fulltexts/PJR/42604/PJR-85-42604-g003_min.jpg
Figure 4

Number of articles about intravoxel incoherent motion magnetic resonance imaging in oncology 1986-2018 (PubMed)

/f/fulltexts/PJR/42604/PJR-85-42604-g004_min.jpg

Recently, there has been greater interest in using the IVIM technique to examine organ function and to evaluate treatment response. There is growing evidence that the additional functional information related to tumour angiogenesis may assist diagnosis and could have prognostic value in oncology patients.

Applications of IVIM may include detection and characterisation of widespread disease, such as both primary and metastatic tumours.

Published studies indicate that, in the near future, IVIM may become an important diagnostic imaging modality to evaluate brain, head and neck, salivary gland, pancreatic, prostate, breast, liver, and renal lesions, as well as lymph node metastasis [2].

However, the main problem with implementation of IVIM for routine MR imaging is lack of standardisation of both image acquisition and analysis, resulting in significant variance in calculated parameters among studies [75,76] (Table 4). Thus, this technique still needs to be refined in terms of acquisition, analysis, and application.

Table 4

Characteristics of the included studies

AutorLesionParameterParameter valueSensitivitySpecificityAccuracyAUCP value
Chandarana et al. [32]Clear cell RCCADC × 10-3 mm2/s2.2 ± 0.492.9%58.3%76.9%0.67
D × 10-3 mm2/s1.8 ± 0.686%62.5%75.0%0.68
f0.25 ± 0.05100%62.5%82.7%0.74
D* × 10-3 mm2/s19 ± 13.50,53
Cystic RCC and papillary RCCf< 0.16100%
Clear cell RCC and chromophobe RCCf> 0.16100%
Ding et al. [35]Papillary RCC and chromophobe RCC vs. fat poor AMLADC> 1.3939.1%100%60.0%0.6340.167
D> 0.9756.5%100%71.4%0.7570.002
D*≤ 28.0387.0%75%82.9%0.822< 0.001
f≤ 13.6143.5%100%62.9%0.783< 0.001
Sumi et al. [38]Normal parotid glandD, × 10-3 mm2/s0.765 ± 0.138
f0.148 ± 0.043
D*, × 10-3 mm2/s41.68 ± 21.91
Normal submandibular glandD, × 10-3 mm2/s0.999 ± 0.154
f0.133 ± 0.037
D*, × 10-3 mm2/s54.16 ± 23.24
Malignant salivary gland tumorD, × 10-3 mm2/s0.96 ± 0.22
f0.103 ± 0.0050
D*, × 10-3 mm2/s21.99 ± 19.01
Pleomorphic adenomaD, × 10-3 mm2/s1.38 ± 0.30
f0.066 ± 0.031
D*, × 10-3 mm2/s10.53 ± 3.48
Warthin tumorD, × 10-3 mm2/s0.61 ± 0.11
f0.156 ± 0.039
D*, × 10-3 mm2/s42.64 ± 20.17
Discriminating pleomorphic adenomas from Warthin tumorsD, × 10-3 mm2/s≥ 1.1100%100%100%
f≤ 0.183%100%90%
D*, × 10-3 mm2/s≤ 23100%100%100%
Discriminating malignant salivary gland tumors from benign onesD, × 10-3 mm2/s0.8 < D < 1.164%100%87%
D*, × 10-3 mm2/s10 ≤ D*≤ 2373%65%68%
Combined D and D*, × 10-3 mm2/s0.8 < D < 1.1, D ≤ 0,8 and D* ≤ 23, or 1.1 ≤ D < 1.4 and D* ≥ 12100%100%100%
Sumi et al. [39]LymphomasPP perfusion-related parameter0.09 ± 0.04
D, × 10-3 mm2/s0.47 ± 0.07
Malignant salivary gland tumorsPP0.22 ± 0.07
D, × 10-3 mm2/s1.03 ± 0.16
Squamous cell carcinomasPP0.15 ± 0.04
D, × 10-3 mm2/s0.82 ± 0.17
Pleomorphic adenomasPP0.13 ± 0.02
D, × 10-3 mm2/s1.44 ± 0.39
SchwannomasPP0.23 ± 0.08
D, × 10-3 mm2/s1.26 ± 0.20
Warthin tumorsPP0.19 ± 0.04
D, × 10-3 mm2/s0.73 ± 0.22
Liang L et AL. (42)Malignant lymph nodes vs benign lymph nodesD*, × 10-3 mm2/s120.89 ± 26.94 vs. 68.78 ± 17.720.0001
D, × 10-3 mm2/s0.57 ± 0.12 vs. 0.74 ± 0.210.0001
f (%)0.20 ± 0.02 vs. 0.32 ± 0.050.0001
ADC, × 10-3 mm2/s1.08 ± 0.26 vs. 1.25 ± 0.190.035
Federau et al. [57]Acute ischemic stroke - brain infarct core vs. contralateral sidef0.026 ± 0.019 vs. 0.056 ± 0.0252.2 × 10−6
D, mm2/s3.9 ± 0.79 × 10-4 vs. 7.5 ± 0.86 × 10-41.3 × 10-20
Suo et al. [58]Ischemic stroke vs contralateral hemisphereADC, × 10-3 mm2/s0.43 ± 0.10 vs. 0.73 ± 0.07< 0.001
D, × 10-3 mm2/s0.42 ± 0.10 vs. 0.72 ± 0.07< 0.001
D*, × 10-3 mm2/s10.20 ± 4.17 vs. 10.87 ± 4.750.218
f (%)4.29 ± 2.01 vs. 7.97 ± 2.03< 0.001
fD*, × 10-3 mm2/s0.49± 0.27 vs. 0.94 ± 0.42< 0.001
Togao et al. [64]Differentiating high-grade gliomas from low-grade gliomasD, × 10-3 mm2/s≤ 1.25100%56.2%0.78
ADC, × 10-3 mm2/s≤ 1.2996.6%50.0%0.73
D*, × 10-3 mm2/s≥ 8.4375.9%50.0%0.60
f (%)≥ 7.796.6%81.2%0.95
Suh et al. [67]Glioblastoma group vs. atypical primary central nervous system lymphomafmaxreader 1: 0.101 ± 0.016
vs. 0.021 ± 0.010
reader 2: 0.107 ± 0.024
vs. 0.027 ± 0.015
reader 1: 89.5%
reader 2: 84.2%
reader 1: 95.1%
reader 2: 95.1%
reader 1: < 0.001
reader 2: < 0.001
Liu et al. [10]Breast malignant tumorADC, × 10-3 mm2/s0.95 (0.83, 1.06)93%90%
D, × 10-3 mm2/s0.85 (0.77, 0.98)90%93%
f10.34 (7.68, 11.88)88%54%
D*,× 10-3 mm2/s94.71 (70.33, 113.23)85%41%
Breast benign lesionsADC, × 10-3 mm2/s1.39 (1.32, 1.50)95%90%
D, × 10-3 mm2/s1.99 (1.77, 2.03)92%96%
f6.83 (4.72, 10.33)89%56%
D*, × 10-3 mm2/s107.49 (83.20, 131.19)82%44%
Breast simple cystADC, × 10-3 mm2/s1.96 (1.73, 2.18)
D, × 10-3 mm2/s1.35 (1.26, 1.44)
f1.69 (0.70, 3.47)
D*, × 10-3 mm2/s99.33 (87.04, 155.69)
Valerio et al. [24]Prostate cancerADC, × 10-3 mm2/s0.76 ± 0.27
D, × 10-3 mm2/s0.99 ± 0.38
f9.35 ± 5.97
D*, × 10-3 mm2/s15.56 ±12.91
Yamada et al. [28]HCCADC, × 10-3 mm2/s1.10 ± 0.18
D, × 10-3 mm2/s1.02 ±10.17
f0.15 ±10.07
D*, × 10-3 mm2/s
Liver metastasisADC, × 10-3 mm2/s1.26 ± 0.25
D, × 10-3 mm2/s1.16 ± 0.18
f0.22 ± 0.09
D*, × 10-3 mm2/s
Liver hemangiomaADC, × 10-3 mm2/s1.56 ± 0.22
D, × 10-3 mm2/s1.31 ± 0.21
f0.35 ± 0.10
D*, × 10-3 mm2/s
Liver cystADC, × 10-3 mm2/s3.01 ± 0.28
D, × 10-3 mm2/s3.03 ±0.22
f0.00 ± 0.01
De Robertis et al. [29]Pancreatic cancerADC, × 10-3 mm2/s1.41 (1.02-1.73)
D, × 10-3 mm2/s1.42 (0.87-1.97)77%88%0.821
f5.82 (2.12-18.84)91%100%0.989
D*, × 10-3 mm2/s10.23 (1.79-57.85)94%96%0.952
Neuroendocrine neoplasmsADC, × 10-3 mm2/s1.28 (1.01-1.89)
D, × 10-3 mm2/s1.2 (0.54-1.56)
f27.61 (3.22–80.83)
D*, × 10-3 mm2/s28.9 (12.16-73.08)
Autoimmune pancreatitisADC, × 10-3 mm2/s1.24 (0.94-1.64)
D, × 10-3 mm2/s1.16 (0.84-1.76)
f9.87 (4.18-11.60)100%100%1.000
D*, × 10-3 mm2/s10.62 (7.6-15.15)100%98%0.988

[i] ADC - apparent diffusion coefficient; D - diffusion coefficient; D* - pseudodiffusion coefficient; f - perfusion fraction; AUC - area under curve; RCC - renal cell carcinoma; AML- angiomyolipoma; HCC - hepatocellular carcinoma; PP - perfusion-related parameter.

Conflict of interest

The authors report no conflict of interest.

References

1 

Le Bihan D, Breton E, Lallemand D, et al. Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology 1988; 168: 497-505.

2 

Iima M, Le Bihan D. Clinical intravoxel incoherent motion and diffusion MR imaging: past, present, and future. Radiology 2016; 278: 13-32.

3 

Li YT, Cercueil JP, Yuan J, et al. Liver intravoxel incoherent motion (IVIM) magnetic resonance imaging: a comprehensive review of published data on normal values and applications for fibrosis and tumor evaluation. Quant Imaging Med Surg 2017; 7: 59-78.

4 

Le Bihan D, Johansen-Berg H. Diffusion MRI at 25: exploring brain tissue structure and function. Neuroimage 2012; 61: 324-341.

5 

Le Bihan D. What can we see with IVIM MRI? Neuroimage 2019; 187: 56-67.

6 

Park HJ, Sung YS, Lee SS, et al. Intravoxel incoherent motion diffusion-weighted MRI of the abdomen: the effect of fitting algorithms on the accuracy and reliability of the parameters. J Magn Reson Imaging 2017; 45: 1637-1647.

7 

Lemke A, Stieltjes B, Schad L, Laun F. Toward an optimal distribution of b values for intravoxel incoherent motion imaging. Magn Reson Imaging 2011; 29: 766-776.

8 

Barbieri S, Donati OF, Froehlich JM, et al. Comparison of intravoxel incoherent motion parameters across MR imagers and field strengths: evaluation in upper abdominal organs. Radiology 2016; 279: 784-794.

9 

Guo W, Luo D, Lin M, et al. Pretreatment intra-voxel incoherent motion diffusion-weighted imaging (IVIM-DWI) in predicting induction chemotherapy response in locally advanced hypopharyngeal carcinoma. Medicine (Baltimore) 2016; 95: e3039.

10 

Liu C, Liang C, Liu Z, et al. Intravoxel incoherent motion (IVIM) in evaluation of breast lesions: comparison with conventional DWI. Eur J Radiol 2013; 82: e782-789.

11 

Ichikawa S, Motosugi U, Ichikawa T, et al. Intravoxel incoherent motion imaging of the kidney: alterations in diffusion and perfusion in patients with renal dysfunction. Magn Reson Imaging 2013; 31: 414-417.

12 

Rheinheimer S, Stieltjes B, Schneider F, et al. Investigation of renal lesions by diffusion-weighted magnetic resonance imaging applying intravoxel incoherent motion-derived parameters–initial experience. Eur J Radiol 2012; 81: e310-316.

13 

Cui Y, Dyvorne H, Besa C, et al. IVIM diffusion-weighted imaging of the liver at 3.0T: comparison with 1.5T. Eur J Radiol Open 2015; 2: 123-128.

14 

Dyvorne HA, Galea N, Nevers T, et al. Diffusion-weighted imaging of the liver with multiple b values: effect of diffusion gradient polarity and breathing acquisition on image quality and intravoxel incoherent motion parameters–a pilot study. Radiology 2013; 266: 920-929.

15 

Federau C, O’Brien K, Meuli R, et al. Measuring brain perfusion with intravoxel incoherent motion (IVIM): initial clinical experience. J Magn Reson Imaging 2014; 39: 624-632.

16 

Federau C, Meuli R, O’Brien K, et al. Perfusion measurement in brain gliomas with intravoxel incoherent motion MRI. AJNR Am J Neuroradiol 2014; 35: 256-262.

17 

Döpfert J, Lemke A, Weidner A, Schad LR. Investigation of prostate cancer using diffusion-weighted intravoxel incoherent motion imaging. Magn Reson Imaging 2011; 29: 1053-1058.

18 

Kakite S, Dyvorne HA, Lee KM, et al. Hepatocellular carcinoma: IVIM diffusion quantification for prediction of tumor necrosis compared to enhancement ratios. Eur J Radiol Open 2015; 3: 1-7.

19 

Shim WH, Kim HS, Choi CG, Kim SJ. Comparison of apparent diffusion coefficient and intravoxel incoherent motion for differentiating among glioblastoma, metastasis, and lymphoma focusing on diffusion-related parameter. PLoS One 2015; 10: e0134761.

20 

Conklin J, Heyn C, Roux M, et al. A simplified model for intravoxel incoherent motion perfusion imaging of the brain. AJNR Am J Neuroradiol 2016; 37: 2251-2257.

21 

Federau C, Cerny M, Roux M, et al. IVIM perfusion fraction is prognostic for survival in brain glioma. Clin Neuroradiol 2017; 27: 485-492.

22 

Xu XQ, Choi YJ, Sung YS, et al. Intravoxel incoherent motion MR imaging in the head and neck: correlation with dynamic contrast-enhanced MR imaging and diffusion-weighted imaging. Korean J Radiol 2016; 17: 641-649.

23 

Boss A, Barth B, Filli L, et al. Simultaneous multi-slice echo planar diffusion weighted imaging of the liver and the pancreas: optimization of signal-to-noise ratio and acquisition time and application to intravoxel incoherent motion analysis. Eur J Radiol 2016; 85: 1948-1955.

24 

Valerio M, Zini C, Fierro D, et al. 3T multiparametric MRI of the prostate: does intravoxel incoherent motion diffusion imaging have a role in the detection and stratification of prostate cancer in the peripheral zone? Eur J Radiol 2016; 85: 790-794.

25 

Pang Y, Turkbey B, Bernardo M, et al. Intravoxel incoherent motion MR imaging for prostate cancer: an evaluation of perfusion fraction and diffusion coefficient derived from different b-value combinations. Magn Reson Med 2013; 69: 553-562.

26 

Bane O, Wagner M, Zhang JL, et al. Assessment of renal function using intravoxel incoherent motion diffusion-weighted imaging and dynamic contrast-enhanced MRI. J Magn Reson Imaging 2016; 44: 317-326.

27 

Luciani A, Vignaud A, Cavet M, et al. Liver cirrhosis: intravoxel incoherent motion MR imaging–pilot study. Radiology 2008; 249: 891-899.

28 

Yamada I, Aung W, Himeno Y, et al. Diffusion coefficients in abdominal organs and hepatic lesions: evaluation with intravoxel incoherent motion echo-planar MR imaging. Radiology 1999; 210: 617-623.

29 

De Robertis R, Cardobi N, Ortolani S, et al. Intravoxel incoherent motion diffusion-weighted MR imaging of solid pancreatic masses: reliability and usefulness for characterization. Abdom Radiol (NY) 2019; 44: 131-139.

30 

Chiaradia M, Baranes L, Van Nhieu JT, et al. Intravoxel incoherent motion (IVIM) MR imaging of colorectal liver metastases: are we only looking at tumor necrosis? J Magn Reson Imaging 2014; 39: 317-325.

31 

Re TJ, Lemke A, Klauss M, et al. Enhancing pancreatic adenocarcinoma delineation in diffusion derived intravoxel incoherent motion f-maps through automatic vessel and duct segmentation. Magn Reson Med 2011; 66: 1327-1332.

32 

Chandarana H, Kang SK, Wong S, et al. Diffusion-weighted intravoxel incoherent motion imaging of renal tumors with histopathologic correlation. Invest Radiol 2012; 47: 688-696.

33 

Campbell N, Rosenkrantz AB, Pedrosa I. MRI phenotype in renal cancer: is it clinically relevant? Top Magn Reson Imaging 2014; 23: 95-115.

34 

Wu Y, Kwon YS, Labib M, et al. Magnetic resonance imaging as a biomarker for renal cell carcinoma. Dis Markers 2015; 2015: 648495.

35 

Ding Y, Zeng M, Rao S, et al. Comparison of biexponential and monoexponential model of diffusion-weighted imaging for distinguishing between common renal cell carcinoma and fat poor angiomyolipoma. Korean J Radiol 2016; 17: 853-863.

36 

Johnson LM, Turkbey B, Figg WD, Choyke PL. Multiparametric MRI in prostate cancer management. Nat Rev Clin Oncol 2014; 11: 346-353.

37 

Shinmoto H, Tamura C, Soga S, et al. An intravoxel incoherent motion diffusion-weighted imaging study of prostate cancer. AJR Am J Roentgenol 2012; 199: W496-500.

38 

Sumi M, Van Cauteren M, Sumi T, et al. Salivary gland tumors: use of intravoxel incoherent motion MR imaging for assessment of diffusion and perfusion for the differentiation of benign from malignant tumors. Radiology 2012; 263: 770-777.

39 

Sumi M, Nakamura T. Head and neck tumors: assessment of perfusion-related parameters and diffusion coefficients based on the intravoxel incoherent motion model. AJNR Am J Neuroradiol 2013; 34: 410-416.

40 

Noij DP, Martens RM, Marcus JT, et al. Intravoxel incoherent motion magnetic resonance imaging in head and neck cancer: A systematic review of the diagnostic and prognostic value. Oral Oncol 2017; 68: 81-91.

41 

Hejduk B, Bobek-Billewicz B, Rutkowski T, et al. Application of intravoxel incoherent motion (IVIM) model for differentiation between metastatic and non-metastatic head and neck lymph nodes. Pol J Radiol 2017; 82: 506-510.

42 

Liang L, Luo X, Lian Z, et al. Lymph node metastasis in head and neck squamous carcinoma: efficacy of intravoxel incoherent motion magnetic resonance imaging for the differential diagnosis. Eur J Radiol 2017; 90: 159-165.

43 

Marzi S, Piludu F, Sanguineti G, et al. The prediction of the treatment response of cervical nodes using intravoxel incoherent motion diffusion-weighted imaging. Eur J Radiol 2017; 92: 93-102.

44 

Hauser T, Essig M, Jensen A, et al. Prediction of treatment response in head and neck carcinomas using IVIM-DWI: evaluation of lymph node metastasis. Eur J Radiol 2014; 83: 783-787.

45 

Marzi S, Forina C, Marucci L, et al. Early radiation-induced changes evaluated by intravoxel incoherent motion in the major salivary glands. J Magn Reson Imaging 2015; 41: 974-982.

46 

Ding Y, Hazle JD, Mohamed AS, et al. Intravoxel incoherent motion imaging kinetics during chemoradiotherapy for human papillomavirus-associated squamous cell carcinoma of the oropharynx: preliminary results from a prospective pilot study. NMR Biomed 2015; 28: 1645-1654.

47 

Liu C, Liang C, Liu Z, et al. Intravoxel incoherent motion (IVIM) in evaluation of breast lesions: comparison with conventional DWI. Eur J Radiol 2013; 82: e782-e789.

48 

Iima M, Yano K, Kataoka M, et al. Quantitative non-Gaussian diffusion and intravoxel incoherent motion magnetic resonance imaging: differentiation of malignant and benign breast lesions. Invest Radiol 2015; 50: 205-211.

49 

Sigmund EE, Cho GY, Kim S, et al. Intravoxel incoherent motion imaging of tumor microenvironment in locally advanced breast cancer. Magn Reson Med 2011; 65: 1437-1447.

50 

Iima M, Le Bihan D, Okumura R, et al. Apparent diffusion coefficient as an MR imaging biomarker of low-risk ductal carcinoma in situ: a pilot study. Radiology 2011; 260: 364-372.

51 

Kamitani T, Hatakenaka M, Yabuuchi H, et al. Detection of axillary node metastasis using diffusion-weighted MRI in breast cancer. Clin Imaging 2013; 37: 56-61.

52 

Luo N, Su D, Jin G, et al. Apparent diffusion coefficient ratio between axillary lymph node with primary tumor to detect nodal metastasis in breast cancer patients. J Magn Reson Imaging 2013; 38: 824-828.

53 

Ye X, Chen S, Tian Y, et al. A preliminary exploration of the intravoxel incoherent motion applied in the preoperative evaluation of mediastinal lymph node metastasis of lung cancer. J Thorac Dis 2017; 9: 1073-1080.

54 

Qiu L, Liu XL, Liu SR, et al. Role of quantitative intravoxel incoherent motion parameters in the preoperative diagnosis of nodal metastasis in patients with rectal carcinoma. J Magn Reson Imaging 2016; 44: 1031-1039.

55 

Yu XP, Wen L, Hou J, et al. Discrimination between metastatic and nonmetastatic mesorectal lymph nodes in rectal cancer using intravoxel incoherent motion diffusion-weighted magnetic resonance imaging. Acad Radiol 2016; 23: 479-485.

56 

Wirestam R, Brockstedt S, Lindgren A, et al. The perfusion fraction in volunteers and in patients with ischaemic stroke. Acta Radiol 1997; 38: 961-964.

57 

Federau C, Sumer S, Becce F, et al. Intravoxel incoherent motion perfusion imaging in acute stroke: initial clinical experience. Neuroradiology 2014; 56: 629-635.

58 

Suo S, Cao M, Zhu W, et al. Stroke assessment with intravoxel incoherent motion diffusion-weighted MRI. NMR Biomed 2016; 29: 320-328.

59 

Yao Y, Zhang S, Tang X, et al. Intravoxel incoherent motion diffusion-weighted imaging in stroke patients: initial clinical experience. Clin Radiol 2016; 71: 938.e11-16.

60 

Heit JJ, Wintermark M, Martin BW, et al. Reduced intravoxel incoherent motion microvascular perfusion predicts delayed cerebral ischemia and vasospasm after aneurysm rupture. Stroke 2018; 49: 741-745.

61 

Federau C, Maeder P, O’Brien K, Browaeys P, Meuli R, Hagmann P. Quantitative measurement of brain perfusion with intravoxel incoherent motion MR imaging. Radiology 2012; 265: 874-881.

62 

Paschoal AM, Leoni RF, Dos Santos AC, Paiva FF. Intravoxel incoherent motion MRI in neurological and cerebrovascular diseases. Neuroimage Clin 2018; 20: 705-714.

63 

Louis DN, Perry A, Reifenberger G, et al. The 2016 World Health Organization classification of tumors of the central nervous system: a summary. Acta Neuropathol 2016; 131: 803-820.

64 

Togao O, Hiwatashi A, Yamashita K, et al. Differentiation of high-grade and low-grade diffuse gliomas by intravoxel incoherent motion MR imaging. Neuro Oncol 2016; 18: 132-141.

65 

Shen N, Zhao L, Jiang J, et al. Intravoxel incoherent motion diffusion-weighted imaging analysis of diffusion and microperfusion in grading gliomas and comparison with arterial spin labeling for evaluation of tumor perfusion. J Magn Reson Imaging 2016; 44: 620-632.

66 

Puig J, Sanchez-Gonzalez J, Blasco G, et al. Intravoxel incoherent motion metrics as potential biomarkers for survival in glioblastoma. PLoS One 2016; 11: e0158887.

67 

Suh CH, Kim HS, Lee SS, et al. Atypical imaging features of primary central nervous system lymphoma that mimics glioblastoma: utility of intravoxel incoherent motion MR imaging. Radiology 2014; 272: 504-513.

68 

Yamashita K, Hiwatashi A, Togao O, et al. Diagnostic utility of intravoxel incoherent motion mr imaging in differentiating primary central nervous system lymphoma from glioblastoma multiforme. J Magn Reson Imaging 2016; 44: 1256-1261.

69 

Togao O, Hiwatashi A, Yamashita K, et al. Measurement of the perfusion fraction in brain tumors with intravoxel incoherent motion MR imaging: validation with histopathological vascular density in meningiomas. Br J Radiol 2018; 91: 20170912.

70 

Cui Y, Zhang C, Li X, et al. Intravoxel incoherent motion diffusion-weighted magnetic resonance imaging for monitoring the early response to ZD6474 from nasopharyngeal carcinoma in nude mouse. Sci Rep 2015; 5: 16389.

71 

Joo I, Lee JM, Grimm R, et al. Monitoring vascular disrupting therapy in a rabbit liver tumor model: relationship between tumor perfusion parameters at IVIM diffusion-weighted MR imaging and those at dynamic contrast-enhanced MR imaging. Radiology 2016; 278: 104-113.

72 

Joo I, Lee JM, Han JK, Choi BI. Intravoxel incoherent motion diffusion-weighted MR imaging for monitoring the therapeutic efficacy of the vascular disrupting agent CKD-516 in rabbit VX2 liver tumors. Radiology 2014; 272: 417-426.

73 

Yang SH, Lin J, Lu F, et al. Evaluation of antiangiogenic and antiproliferative effects of sorafenib by sequential histology and intravoxel incoherent motion diffusion-weighted imaging in an orthotopic hepatocellular carcinoma xenograft model. J Magn Reson Imaging 2017; 45: 270-280.

74 

Detsky JS, Keith J, Conklin J, et al. Differentiating radiation necrosis from tumor progression in brain metastases treated with stereotactic radiotherapy: utility of intravoxel incoherent motion perfusion MRI and correlation with histopathology. J Neurooncol 2017; 134: 433-441.

75 

Yablonskiy DA, Sukstanskii AL. Theoretical models of the diffusion weighted MR signal. NMR Biomed 2010; 23: 661-681.

76 

Cieszanowski A, Pasicz K, Podgórska J, et al. Reproducibility of intravoxel incoherent motion of liver on a 3.0T scanner: free-breathing and respiratory-triggered sequences acquired with different numbers of excitations. Pol J Radiol 2018; 83: e437-e445.

Copyright: © Polish Medical Society of Radiology This is an Open Access article distributed under the terms of the Creative Commons Attribution-Noncommercial-No Derivatives 4.0 International (CC BY-NC-ND 4.0). License allowing third parties to download articles and share them with others as long as they credit the authors and the publisher, but without permission to change them in any way or use them commercially.
 
Quick links
© 2024 Termedia Sp. z o.o.
Developed by Bentus.