Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Can the low and high b-value distribution influence the pseudodiffusion parameter derived from IVIM DWI in normal brain?

BMC Medical Imaging volume 20, Article number: 14 (2020) Cite this article

Abstract

Background

Our study aims to reveal whether the low b-values distribution, high b-values upper limit, and the number of excitation (NEX) influence the accuracy of the intravoxel incoherent motion (IVIM) parameter derived from multi-b-value diffusion-weighted imaging (DWI) in the brain.

Methods

This prospective study was approved by the local Ethics Committee and informed consent was obtained from each participant. The five consecutive multi-b DWI with different b-value protocols (0–3500 s/mm2) were performed in 22 male healthy volunteers on a 3.0-T MRI system. The IVIM parameters from normal white matter (WM) and gray matter (GM) including slow diffusion coefficient (D), fast perfusion coefficient (D*) and perfusion fraction (f) were compared for differences among defined groups with different IVIM protocols by one-way ANOVA.

Results

The D* and f value of WM or GM in groups with less low b-values distribution (less than or equal to 5 b-values) were significantly lower than ones in any other group with more low b-values distribution (all P <  0.05), but no significant differences among groups with more low b-values distribution (P > 0.05). In addition, no significant differences in the D, D* and f value of WM or GM were found between group with one and more NEX of low b-values distribution (all P > 0.05). IVIM parameters in normal WM and GM strongly depended on the choice of the high b-value upper limit.

Conclusions

Metrics of IVIM parameters can be affected by low and high b value distribution. Eight low b-values distribution with high b-value upper limit of 800–1000 s/mm2 may be the relatively proper set when performing brain IVIM studies.

Peer Review reports

Background

As an advanced magnetic resonance imaging (MRI) technique, diffusion-weighted imaging (DWI) is considered the one of the most sensitive pulse sequences for early pathological changes. However, previous studies have shown that perfusion can substantially confound diffusion measurements because of the incoherent motion of blood in pseudorandom capillary networks at the macroscopic level, and have proposed intravoxel incoherent motion (IVIM) imaging to measure microvascular perfusion [1]. For its ability of utilizing biexponential model to extract the perfusion-related information from a diffusion sequence [1], IVIM DWI has been successfully applied to clinical research with the great improvement of gradient coil systems for diffusion MR imaging over the recent years, especially in the body or head and neck [2,3,4,5,6,7,8,9,10,11]. Perfusion has been proven an important surrogate marker of many physiologic or pathologic processes as well as a predictor of recovery with reperfusion in patients of acute ischemic stroke [12]. Also, perfusion measurement has demonstrated improved sensitivity and predictive value for tumor grading and prognosis [13,14,15,16,17,18,19].

The IVIM model is a two-compartment model and includes terms for the fraction of received signal attributed to moving blood (fractional perfusion, f), the diffusion caused by moving blood (pseudodiffusion, D*), and a diffusion component free of perfusion effects (true molecular diffusion, D). The model requires the collection of both low and high b-values as the perfusion effect becomes largely negligible as the b-value is increased beyond approximately 200 s/mm2 [1, 20, 21].

Previous studies have supported that IVIM DWI could be potentially useful in measuring brain perfusion [22,23,24,25,26,27,28,29,30]. However, the calculated perfusion-related IVIM maps usually suffer from low image quality and high parametric variance [31]. Many factors might influence the accuracy of pseudodiffusion parameter derived from IVIM DWI, such as the low b-values distribution, number of excitation (NEX), repetition time (TR), echo time (TE), field strength, and diffusion preparation pulses [32,33,34,35,36]. To date, the b value distributions are chosen heuristically and vary greatly among researchers.

In spite of an increasing number of applications of IVIM DWI, there is no clear consensus regarding the optimal protocol that should be used. Concerning the used acquisition parameters, the amount and distribution of b values varied greatly. For example, IVIM studies in the glioma grading have used different b-value distributions and obtained different trend values for IVIM parameters [22, 37,38,39,40]. Specifically, using 13 b values (0–1000 s/mm2), Togao and his colleagues found significantly lower D and higher f in high grade glioma (HGG) than in low grade glioma (LGG), while no statistically difference of D* in two groups [37]. In another research, in which 20 b values were used (0–3500 s/mm2), reported a significant higher D* was found in HGG group [38]. Recent studies have shown pseudodiffusion in the liver tended to be underestimated when too few low b-values (0 < b < 50 s/mm2) were included in the distribution [41]. However, to our knowledge, the optimal b-value distribution for IVIM DWI of the brain is scarce.

Ideally speaking, dense sampling of b-values would certainly guarantee quantification data quality and repeatability. However, more b-value samples involve longer scan time, making the measurement less practical and vulnerable to subject motion. We hypothesize that the amount of b-values could be reduced while still enabling correct IVIM parameter estimation, without affecting the reproducibility of the technique.

Our study aims to reveal whether the low b-values distribution, high b-values upper limit, and NEX for low b-values influence the accuracy of IVIM parameter in the brain, and to determine the minimal amount and optimal b-value distribution necessary for reproducible brain IVIM parameter quantification.

Methods

Subjects

Between November and December 2014, 22 male healthy volunteers were recruited (mean age, 28 years; range, 25–34) with inclusion standard as follows: no potential vascular risk factors or diseases affecting brain microvascular perfusion including hypertension or cerebral vascular diseases, systemic metabolic disease, such as diabetes, obesity or cirrhosis; no infection or fever, and use of corticosteroid drugs; no MRI contraindication and appearing normal on conventional brain MR image (as shown in Fig. 1).

Fig. 1
figure 1

Flow diagram shows subjects selection process and multi-b-value DWI protocol. Note: Each set of low b-values was in addition to a high b-value distribution of 300, 500, 800, 1000, 1500, 2000, 3000 and 3500 s/mm2, with NEX of 2,2,2,2,2,4,4 and 6, respectively

Full size image

Brain MR image data acquisition

The whole brain MRI examinations were performed on a 3.0-T MRI system (Discovery MR750, GE Healthcare, Milwaukee, WI, USA) with a 40-mT/m maximum gradient capability and an eight-channel head coil (GE Medical Systems). Conventional MRI sequence and five consecutive multiple b-value DWI sequences of brain were performed during one examination.

Conventional MRI sequences included T1-weighted spin echo in the transverse plane (TR/TE, 1750 ms/24 ms; matrix size, 256 × 256; field of view (FOV), 24 cm × 24 cm; NEX, 1; slice thickness, 5 mm; gap, 1.5 mm), T2-weighted fast spin echo in the transverse planes (TR/TE, 4247 ms/93 ms; matrix size, 512 × 512; FOV, 24 cm × 24 cm; NEX, 1; slice thickness, 5 mm; gap, 1.5 mm) and sagittal planes (TR/TE, 10,639 ms/96 ms; matrix size, 384 × 384; FOV, 24 cm × 24 cm; NEX, 2; slice thickness, 5 mm; gap, 1.0 mm), and fluid-attenuated inversion recovery (FLAIR) in the transverse plane (TR/TE, 8000 ms/165 ms; matrix size, 256 × 256; FOV, 24 cm × 24 cm; NEX, 1; slice thickness, 5 mm; gap, 1.5 mm) [22].

Five consecutive multi-b-value DWI sequences of brain were performed after the conventional MRI. The details of five different parameter sets of low b-value distribution and NEX were showed in Fig. 1. Different parameter sets of b values were applied with a single-shot diffusion-weighted spin-echo echo-planar sequence. Parallel imaging was used with an acceleration factor of 2. A local shim box covering the whole brain was applied to minimize susceptibility artifacts. In total, 20 axial slices covering the entire brain were obtained with a 24 cm × 24 cm FOV, 5 mm slice thickness, 1.5 mm slice gap, 3000 ms TR, Minimum TE, 128 × 128 matrix and one diffusion preparation pulses [22].

IVIM DWI data processing

All data were analyzed and processed on a GE ADW4.6 workstation (MR750, GE Healthcare, Milwaukee, WI, USA), using the MADC program of Functool software. The mean IVIM parameters were measured independently by one experienced radiologist (Y.-C.,H, with 12 years of experience in radiology).

According to IVIM theory, the relationship between signal intensity and b values can be expressed based on an Eq. (1) [22]:

$$ {\mathrm{S}}_{\mathrm{b}}/{\mathrm{S}}_0=\mathrm{f}\ \exp \left(-\mathrm{b}\ \mathrm{D}\ast \right)+\left(1-\mathrm{f}\right)\exp \left(-\mathrm{b}\ \mathrm{D}\right) $$
(1)

where S0 = signal intensity at the b value of 0 s/mm2; Sb = signal intensity at the b value denoted by the subscript; D is the slow diffusion component that reflects random motion of intra- and intercellular water molecules; f is the fraction of the diffusion linked to microcirculation, and D* is the fast diffusion component representing incoherent microcirculation within the voxel. D* is usually expected to be at least one order of magnitude higher than D [2, 22, 42], the influence of D* on signal decay can be neglected for b values greater than 200 s/mm2. Eq. (1) can then be simplified, and the estimation of D can be obtained by using only b values greater than 200 s/mm2 with a simple linear fit eq. (2) [22, 43]:

$$ {\mathrm{S}}_{\mathrm{b}}/{\mathrm{S}}_0=\left(1-\mathrm{f}\right)\exp \left(-\mathrm{b}\ \mathrm{D}\right) $$
(2)

The D, D* and f values were calculated according to the bi-exponential model with eq. (1) and (2). First, for high b values (b > 200 s/mm2), Sb was first fitted to eq. (2) using a linear model, and the D value was calculated. In a second step, the f and D* values were calculated by using a nonlinear regression algorithm according to eq. (1), while keeping D constant. Freehand region of interests (ROIs) were placed in right frontal white matter (WM) and gray matter (GM) (as shown in Fig. 2). The mean ROI area was range from 25 to 35 mm2. According to the bi-exponential fitting of the diffusion signal decay over different sets of b values (as shown in Fig. 2 and Table 1), the IVIM parameter maps were generated based on the D, D* and f values derived from Eq. (1) and (2) (as shown in Fig. 3), and the mean D, D*, and f values in the corresponding ROIs were obtained, respectively.

Fig. 2
figure 2

Comparison of the bi- and mono-exponential fitting of the diffusion signal decay over a wide-range of b-values (up to 1000 s/mm2) in normal brain, with more obvious pseudodiffusion effect in GM than in WM. Note: Axial diffusion-weighted trace image (b = 1000 s/mm2) shows ROIs placed in right frontal white matter (WM) and gray matter (GM)

Full size image
Table 1 Defined groups based on IVIM DWI protocol with different b-value distribution in brain
Full size table
Fig. 3
figure 3

Brain pseudodiffusion parametric maps from one representative subject. The top row depicts D* maps (a1-d1) calculated by 4 sets of multiple b-value DWI with different high b-value upper limit (500, 800, 1000, 1500 s/mm2, respectively) and the same low b-value distribution (0, 30, 60, 90, 120, 150, 180, 200 s/mm2), meanwhile the bottom row shows corresponding f maps (a2-d2) derived from IVIM DWI model in normal brain. As the increase of high b-value upper limit, the red color area representing fast diffusion effect in brain tissue was reduced on pseudodiffusion images, reflecting increased fast diffusion effect. At the same time, much red color area was displayed at 500 s/mm2 high b-value upper limit (a1 or a2), and relative absent red area on pseudodiffusion images when high b-value upper limit reached 1500 s/mm2 (d1 or d2)

Full size image

The monoexponential fitting of the diffusion signal decay in WM and GM ROIs (as shown in Fig. 2) was conducted by using the following eq. (3) [44]:

$$ {\mathrm{S}}_{\mathrm{b}}/{\mathrm{S}}_0=\exp \left(-\mathrm{b}\times \mathrm{ADC}\right) $$
(3)

where Sb is the signal intensity in the pixel with diffusion gradient b (all acquired b values), and S0 is the signal intensity in the pixel without diffusion gradient.

Group definition

Table 1 summarizes the defined groups based on IVIM DWI protocol with different b-value distribution. Firstly, we determine the L1 to L4 groups to explore the effect of low b-value distribution on IVIM values, by keeping the high b values (300, 500, 800, and 1000 s/mm2) and changing the amount of low b values with same NEX of 2 for each small b value [45]: group L1 (low b value: 0, 50, 150 and 200 s/mm2), L2 (low b value: 0, 30, 50, 100 and 200 s/mm2), L3 (low b value: 0, 30, 60, 90,120,150,180 and 200 s/mm2), and L4 (low b value: 0, 20, 40, 60, 80, 100, 120, 140, 160, 180 and 200 s/mm2). Secondly, for revealing the effect of NEX on the IVIM parameters, we determine the group L3–1 (low b-value distribution are same with group L3, but with a small NEX of 1 for each low b-values). In order to validate the consistency of IVIM processing, we define the group L3-L (b-value distribution and NEX are completely same with group L3) by placing the ROI in left frontal lobe during data processing. In addition, we try to explore the effect of high b-values upper limit on IVIM metrics, and define group H1 (high b value: 300, 500 and 800 s/mm2), H2 (high b value: 300, 500, 800 and 1000 s/mm2), H3 (high b value: 300, 500, 800, 1000 and 1500 s/mm2), H4 (high b value: 300, 500, 800, 1000, 1500 and 2000 s/mm2), H5 (high b value: 300, 500, 800, 1000, 1500, 2000 and 3000 s/mm2) and H6 (high b value: 300, 500, 800, 1000, 1500, 2000, 3000 and 3500 s/mm2) based on the different high b-value distribution, with a same set of low b value (0, 30, 60, 90, 120, 150, 180 and 200 s/mm2).

Statistical analysis

All statistical analyses were performed with IBM SPSS 20.0 software (IBM Corp, Chicago, IL, USA). Numerical variables were denoted as the mean and standard deviation. D, D* and f in WM or GM were tested for differences among group L1, L2, L3, L3–1, L4 and L3-L, and among group H1, H2, H3, H4, H5 and H6 by one-way ANOVA, and further post hoc multiple comparisons were performed with Bonferroni test (equal variances assumed) and Dunett’s T3 test (equal variances not assumed). D* and f values were compared for the differences between the WM and GM by using independent sample t test. P <  0.05 was considered statistically significant.

Results

Pseudodiffusion parameter of WM or GM might be underestimated using less low b-values distribution

The descriptive statistics of the D, D* and f value of WM or GM among groups of low b-value distribution are shown in Table 2. The D* and f value of WM or GM in group L1 and L2 were significantly lower than ones in any other group L3, L4 and L3-L (all P <  0.05), but no difference exists between group L1 and L2, among L3, L4, or L3-L (all P > 0.05) (D*, 5.05, 4.93, 6.11, 5.84 and 6.24 × 10− 3 mm2/s in WM, and 8.18, 8.38, 9.49, 9.79 and 9.49 × 10− 3 mm2/s in GM; f, 6.56, 6.51, 7.48, 7.64 and 7.53% in WM, and 9.39, 9.26, 11.07, 11.03 and 11.31% in GM, in group L1, L2, L3, L4 and L3-L, respectively, as shown in Fig. 4 and Table 2). There were no differences of D value in WM or GM among each group (all P > 0.05) (as shown in Table 2).

Table 2 Compared the IVIM parameter derived from multiple b-value DWI among groups with different low b-value distribution and NEX in brain
Full size table
Fig. 4
figure 4

Box plots for pseudodiffusion parameter of D* (a) and f (b) derived from IVIM DWI among groups with different distribution of low b-value and NEX in brain white matter (WM) and gray matter (GM)

Full size image

NEX for low b value did not affect metrics of IVIM-DWI in brain

There was no difference in the D, D* and f value of WM or GM between group L3 and L3–1 (D, 0.716 and 0.718 × 10− 3 mm2/s; D*, 6.105 and 5.817 × 10− 3 mm2/s; f, 7.484 and 7.550% in group L3 and L3–1, respectively, P > 0.05) (as shown in Fig. 4 and Table 2).

Parameters comparisons among different groups of high b-value distribution

The descriptive statistics of the D, D* and f value of WM or GM among groups of high b-value distribution are shown in Table 3. The D and D* value both in WM and GM were significantly decreased with the increase of high b-value distribution except D and D* value between group H1 and H2 both in WM and GM, and D* value between group H5 and H6 in GM (D, 0.740, 0.718, 0.647, 0.579, 0.481 and 0.432 × 10− 3 mm2/s in WM, and 0.868, 0.839, 0.758, 0.712, 0.651 and 0.620 × 10− 3 mm2/s in GM; D*, 6.199, 5.450, 4.150, 3.379, 2.533 and 2.365 × 10− 3 mm2/s in WM, and 11.027, 9.264, 6.312, 4.986, 3.869 and 3.377 × 10− 3 mm2/s in GM, in group H1, H2, H3, H4, H5 and H6, respectively, as shown in Fig. 5a and Table 3).

Table 3 Compared the IVIM DWI parameter among groups with different high b-value upper limit in brain
Full size table
Fig. 5
figure 5

Box plots for pseudodiffusion parameter of D* (a) and f (b) derived from IVIM DWI among groups with different high b-value upper limit in brain white matter (WM) and gray matter (GM)

Full size image

The f value was significantly increased with the increase of high b-value distribution, but no significant difference was found between group H1 and H2 both in WM and GM (f, 6.595, 7.098, 12.381, 18.055, 27.827 and 32.241% in WM, and 9.781, 10.225, 13.200, 15.664, 22.445 and 25.277% in GM, in group H1, H2, H3, H4, H5 and H6, respectively, as shown in Fig. 5b and Table 3). In addition, f value in GM were higher than ones in WM both in group H1 and H2, while no significant difference was found in group H3, and inversed results obtained in group H4, H5 and H6 (as shown in Fig. 5b and Table 4).

Table 4 Compared the pseudodiffusion parameter derived from IVIM DWI between normal WM and GM among groups with different high b-value distribution
Full size table

In addition, as shown in Fig. 3, the scattered red color area representing high pseudodiffusion effect in brain tissue was reduced with the increase of high b-value upper limit on pseudodiffusion images, indicating the attenuation of fast diffusion effect.

Discussion

In this study, the influence of low b-values distribution, high b-values upper limit, and NEX for low b-values on IVIM parameters was investigated by in vivo brain multi-b-value DWI. The results of this study indicated that the D* and f value derived from IVIM DWI could be affected when there were less b-values in low b-value distribution (0 < b < 200 s/mm2). We also identified that the NEX for low b-values might not affect the IVIM metrics in brain. In addition, this study demonstrated significant differences in D, D* and f values among different distributions of high b-value derived from brain IVIM DWI of healthy volunteers.

IVIM DWI was found to be a valid and promising method to quantify water molecule diffusion and microvascular perfusion of living tissue perfusion, experiencing a remarkable revival for applications throughout the body, especially for oncologic patients, allowing earlier detection, diagnosis, staging, and monitoring of disease progression or response to therapy [24, 29, 45, 46]. Previously published studies of IVIM imaging have shown that the amount of low b-values significantly influence the pseudodiffusion metrics [41, 45]. Since the IVIM effect is usually small, more images or b-values are often acquired at low b-values than for diffusion at high b values [29]. In this study, D* and f value of WM or GM in more low b-values (more than eight b-values) were significantly larger than that in less low b-values distribution (less than five b-values). An in vivo liver MRI results in the normal controls revealed that the measured pseudodiffusion value was lower when calculated with less low b values (no b values between 0 and 50 s/mm2) compared to when calculated with more low b value (extra b = 10 and 25 s/mm2 were included) in the low b-value distribution [41], which were consistent with current research findings. Moreover, consistent with a recent study [47], there was no significant difference of the mean D* and f value among groups using eight or eleven low b-values distribution. The results indicated that insufficient numbers of low b-values might underestimate pseudodiffusion metrics. Compared with eleven b values, eight low b values (0, 30, 60, 90, 120, 150, 180 and 200 s/mm2) can significantly reduce scan time and may be sufficient for IVIM measurements.

The higher NEX was very important for high b-values to obtain higher signal to noise ratio (SNR) DW images, but it is still controversial whether the higher NEX for low b-values is necessary for calculating pseudodiffusion parameters more accurately in brain [21, 32, 41]. A main body of previous literature suggested that more images of the low b-values less than 200 s/mm2 were required for the calculation of pseudodiffusion, so higher NEX for very low b-values than for the medium b-values should be adopted in IVIM studies [31, 41]. However, we noticed that there was no significant difference of D* and f value derived from IVIM sequence using the same b-values and different NEX in current study. Although it has been confirmed that separation of perfusion from diffusion requires good SNR [21], the NEX for low b value may be less important since the SNR is higher at low b-values in brain. In this study, echo planar imaging (EPI) sequence with multi-channel receive coil and array spacial sensitivity encoding technique (ASSET) was adopted, so the current SNR of low b value DW image with low NEX might enough to satisfy the need of pseudodiffusion measurement [32], and the effect of NEX at low b values on IVIM parameters can be considered as negligible.

In the present study, we also evaluated how the upper limit for the high b-value distribution influences the in vivo IVIM measurements. We detected significant differences in parameters D, D* and f value among majority of high b-value distribution groups both in WM and GM. As the upper limit for the high b-value distribution increased, the D and D* values decreased and f value increased significantly. As a matter of fact, non-Gaussian diffusion becomes visible when high b values are used, and the degree of diffusion related signal attenuation or apparent diffusion coefficient decreases when the high b value increases [29]. It is, thus, mandatory to indicate that D value decreases with the b value increases according to the IVIM model eq. (2), and correspondingly, an accurate D* and f value can be obtained when high b-value upper limit at less than 1000 s/mm2. However, by using such high b values (b value larger than 1500 s/mm2), the IVIM model reaches its limitation, as it cannot give a proper perfusion measurement [29], thus the decreased D value may lead to a biased D* and f value. This phenomenon may explain the contradictory results in many IVIM studies. For example, the f values were significant higher in HGG than in LGG using relatively reasonable upper limit for the high b-value distribution (b value less than 1000 s/mm2) in some studies [37, 40], but other studies using higher upper limit for the high b-value distribution (b value larger than 1500 s/mm2) got inconsistent results in that the f values were significantly lower in HGG than in LGG group [22, 38, 39]. Also demonstrated in current study, the f values were significant lower in WM than in GM when using a high b-value upper limit of less than 1000 s/mm2, but converse results were demonstrated when using a high b-value upper limit higher than 1500 s/mm2. In addition, at high b values, because of the nature of the MR imaging signal (a magnitude signal that cannot be negative), there is always some background noise signal left [29], which may affect IVIM parameter metrics. In summary, DWI model when b-values larger than 1500 s/mm2 is a non-linear fitting curve (diffusion is non-Gaussian), more elements of kurtosis imaging rather than IVIM, we cannot obtain accurate pseudodiffusion metrics by using bi-exponential equation.

It is well known that blood flow in randomly oriented capillaries (at voxel level) mimics a random walk, which results in a pseudodiffusion effect in the presence of diffusion encoding gradient pulses. Indeed, the IVIM imaging has a differential sensitivity to vessel sizes, according to the range of b values that are used [29]. Similarly, in our study, the red color area representing significant pseudodiffusion effect within the brain tissue was reduced with the increase of high b-value upper limit on pseudodiffusion images, indicating the attenuation of fast diffusion effect derived from blood flow in small vessel or cerebrospinal fluid (CSF) space [48]. We noticed that the red color area or fast diffusion effect was relative absent in brain tissue on pseudodiffusion images using higher upper limit for the high b-value distribution (b value higher than 1500 s/mm2). But when a rather low b-value upper limit, for example less than 500 s/mm2, was adopted, considerable partial volume from marked IVIM effect of small vessel and CSF space could influence the pseudiffusion measurement derived from brain micro circulation [48]. Thus, from this point of view, we recommend a proper high b-value upper limit (800–1000 s/mm2) as the referred standard for brain IVIM DWI studies.

One important consideration of using the IVIM technique based on DW MRI is obtaining reliable and repeatable data from patients. In our study, to ensure data accuracy or consistency, we obtained IVIM parameters from contralateral frontal lobe in each participant. Analysis revealed no significant difference of IVIM parameters between two measurements in right and left frontal lobe, which suggested the reliability of the obtained data in our study.

Our study has some limitations. First, no simulation experiments were performed to look at the effects of low b-values on IVIM calculations before in vivo MRI, however, our designs of the b-value selection based on the previous simulation and in vivo research [36, 41, 45, 49]. Second, food intake was not controlled, food intake could potentially influence brain perfusion, and however, pseudodiffusion parameters comparisons were made between different b-value combinations within the same scans and not between subjects, these effects should be minimal in this study. Thirdly, only two-step fitting method was used in this study, another fitting model such as one-step direct fitting technique or Bayesian fitting was not assessed. Comparing with one-step direct fitting technique, the segmented IVIM fitting method was used to increase robustness under biological conditions. Bayesian modeling is capable of producing more visually pleasing IVIM parameter maps than least squares approaches, but their potential to mask certain tissue features demands caution during implementation [50]. Finally, in the current study, hand-drawn ROIs used for each group were not absolutely identical, which may lead to a sampling bias.

Conclusion

This study demonstrated that D* and f value derived from in vivo brain IVIM-DWI could be affected when there were less low b-values distribution (0 < b < 200 s/mm2). Eight low b-values with upper limit for the high b-value distribution of 800–1000 s/mm2 might be the relatively proper set when performing brain IVIM-DWI. We also identified that the NEX selection for low b-value does not influence the brain IVIM metrics.

Availability of data and materials

The datasets analyzed during the current study are not publicly available. These data could only be accessed to researchers to ensure participant confidentiality.

Abbreviations

ANOVA:

Analysis of variance

ASSET:

Array spacial sensitivity encoding technique

CSF:

Cerebrospinal fluid

D:

Slow diffusion coefficient

D*:

Fast perfusion coefficient

DWI:

Diffusion-weighted imaging

EPI:

Echo planar imaging

f:

Perfusion fraction

FLAIR:

Fluid-attenuated inversion recovery

FOV:

Field of view

GM:

Gray matter

HGG:

High grade glioma

IVIM:

Intravoxel incoherent motion

LGG:

Low grade glioma

MRI:

Magnetic resonance imaging

NEX:

Number of excitation

ROI:

Region of interest

SNR:

Signal to noise ratio

TE:

Echo time

TR:

Repetition time

WM:

White matter

References

  1. Le Bihan D, Breton E, Lallemand D, Aubin ML, Vignaud J, Laval-Jeantet M. Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology. 1988;168(2):497–505.

    Article  PubMed  Google Scholar 

  2. Luciani A, Vignaud A, Cavet M, Nhieu JT, Mallat A, Ruel L, Laurent A, Deux JF, Brugieres P, Rahmouni A. Liver cirrhosis: intravoxel incoherent motion MR imaging--pilot study. Radiology. 2008;249(3):891–9.

    Article  PubMed  Google Scholar 

  3. Lemke A, Laun FB, Klauss M, Re TJ, Simon D, Delorme S, Schad LR, Stieltjes B. Differentiation of pancreas carcinoma from healthy pancreatic tissue using multiple b-values: comparison of apparent diffusion coefficient and intravoxel incoherent motion derived parameters. Investig Radiol. 2009;44(12):769–75.

    Article  Google Scholar 

  4. Sigmund EE, Cho GY, Kim S, Finn M, Moccaldi M, Jensen JH, Sodickson DK, Goldberg JD, Formenti S, Moy L. Intravoxel incoherent motion imaging of tumor microenvironment in locally advanced breast cancer. Magn Reson Med. 2011;65(5):1437–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Chandarana H, Kang SK, Wong S, Rusinek H, Zhang JL, Arizono S, Huang WC, Melamed J, Babb JS, Suan EF, et al. Diffusion-weighted intravoxel incoherent motion imaging of renal tumors with histopathologic correlation. Investig Radiol. 2012;47(12):688–96.

    Article  Google Scholar 

  6. Sumi M, Van Cauteren M, Sumi T, Obara M, Ichikawa Y, Nakamura T. 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(3):770–7.

    Article  PubMed  Google Scholar 

  7. Hauser T, Essig M, Jensen A, Laun FB, Munter M, Maier-Hein KH, Stieltjes B. Prediction of treatment response in head and neck carcinomas using IVIM-DWI: evaluation of lymph node metastasis. Eur J Radiol. 2014;83(5):783–7.

    Article  PubMed  Google Scholar 

  8. Kuru TH, Roethke MC, Stieltjes B, Maier-Hein K, Schlemmer HP, Hadaschik BA, Fenchel M. Intravoxel incoherent motion (IVIM) diffusion imaging in prostate cancer - what does it add? J Comput Assist Tomogr. 2014;38(4):558–64.

    Article  PubMed  Google Scholar 

  9. Sumi M, Nakamura T. Head and neck tumours: combined MRI assessment based on IVIM and TIC analyses for the differentiation of tumors of different histological types. Eur Radiol. 2014;24(1):223–31.

    Article  PubMed  Google Scholar 

  10. Woo S, Lee JM, Yoon JH, Joo I, Han JK, Choi BI. Intravoxel incoherent motion diffusion-weighted MR imaging of hepatocellular carcinoma: correlation with enhancement degree and histologic grade. Radiology. 2014;270(3):758–67.

    Article  PubMed  Google Scholar 

  11. Filli L, Wurnig MC, Luechinger R, Eberhardt C, Guggenberger R, Boss A. Whole-body intravoxel incoherent motion imaging. Eur Radiol. 2015;25(7):2049–58.

    Article  PubMed  Google Scholar 

  12. Albers GW, Thijs VN, Wechsler L, Kemp S, Schlaug G, Skalabrin E, Bammer R, Kakuda W, Lansberg MG, Shuaib A, et al. Magnetic resonance imaging profiles predict clinical response to early reperfusion: the diffusion and perfusion imaging evaluation for understanding stroke evolution (DEFUSE) study. Ann Neurol. 2006;60(5):508–17.

    Article  PubMed  Google Scholar 

  13. Law M, Young RJ, Babb JS, Peccerelli N, Chheang S, Gruber ML, Miller DC, Golfinos JG, Zagzag D, Johnson G. Gliomas: predicting time to progression or survival with cerebral blood volume measurements at dynamic susceptibility-weighted contrast-enhanced perfusion MR imaging. Radiology. 2008;247(2):490–8.

    Article  PubMed  Google Scholar 

  14. Baidya Kayal E, Kandasamy D, Khare K, Bakhshi S, Sharma R, Mehndiratta A. Intravoxel incoherent motion (IVIM) for response assessment in patients with osteosarcoma undergoing neoadjuvant chemotherapy. Eur J Radiol. 2019;119:108635.

    Article  PubMed  Google Scholar 

  15. Wang YC, Hu DY, Hu XM, Shen YQ, Meng XY, Tang H, Li Z. Assessing the early response of advanced cervical Cancer to Neoadjuvant chemotherapy using Intravoxel incoherent motion diffusion-weighted magnetic resonance imaging: a pilot study. Chin Med J. 2016;129(6):665–71.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Pieper CC, Sprinkart AM, Meyer C, Konig R, Schild HH, Kukuk GM, Murtz P. Evaluation of a simplified Intravoxel incoherent motion (IVIM) analysis of diffusion-weighted imaging for prediction of tumor size changes and imaging response in breast Cancer liver metastases undergoing Radioembolization: a retrospective single center analysis. Medicine (Baltimore). 2016;95(14):e3275.

    Article  Google Scholar 

  17. Valerio M, Zini C, Fierro D, Giura F, Colarieti A, Giuliani A, Laghi A, Catalano C, Panebianco V. 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(4):790–4.

    Article  PubMed  Google Scholar 

  18. Ye X, Chen S, Tian Y, You B, Zhang W, Zhao Y, Jiang T, Hu B, Li H. 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(4):1073–80.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Xu Q, Xu Y, Sun H, Chan Q, Shi K, Song A, Wang W. Quantitative intravoxel incoherent motion parameters derived from whole-tumor volume for assessing pathological complete response to neoadjuvant chemotherapy in locally advanced rectal cancer. J Magn Reson Imaging. 2018;48(1):248–58.

    Article  PubMed  Google Scholar 

  20. Le Bihan D, Turner R, MacFall JR. Effects of intravoxel incoherent motions (IVIM) in steady-state free precession (SSFP) imaging: application to molecular diffusion imaging. Magn Reson Med. 1989;10(3):324–37.

    Article  PubMed  Google Scholar 

  21. Le Bihan D. Intravoxel incoherent motion perfusion MR imaging: a wake-up call. Radiology. 2008;249(3):748–52.

    Article  PubMed  Google Scholar 

  22. Hu YC, Yan LF, Wu L, Du P, Chen BY, Wang L, Wang SM, Han Y, Tian Q, Yu Y, et al. Intravoxel incoherent motion diffusion-weighted MR imaging of gliomas: efficacy in preoperative grading. Sci Rep. 2014;4:7208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Federau C, Sumer S, Becce F, Maeder P, O'Brien K, Meuli R, Wintermark M. Intravoxel incoherent motion perfusion imaging in acute stroke: initial clinical experience. Neuroradiology. 2014;56(8):629–35.

    Article  CAS  PubMed  Google Scholar 

  24. Federau C, O'Brien K, Meuli R, Hagmann P, Maeder P. Measuring brain perfusion with intravoxel incoherent motion (IVIM): initial clinical experience. J Magn Reson Imaging. 2014;39(3):624–32.

    Article  PubMed  Google Scholar 

  25. Bisdas S, Braun C, Skardelly M, Schittenhelm J, Teo TH, Thng CH, Klose U, Koh TS. Correlative assessment of tumor microcirculation using contrast-enhanced perfusion MRI and intravoxel incoherent motion diffusion-weighted MRI: is there a link between them? NMR Biomed. 2014;27(10):1184–91.

    Article  PubMed  Google Scholar 

  26. Suo S, Cao M, Zhu W, Li L, Li J, Shen F, Zu J, Zhou Z, Zhuang Z, Qu J, et al. Stroke assessment with intravoxel incoherent motion diffusion-weighted MRI. NMR Biomed. 2016;29(3):320–8.

    Article  PubMed  Google Scholar 

  27. Sun J, Yu X, Jiaerken Y, Song R, Huang P, Wang C, Yuan L, Mao Y, Guo Y, Yu H, et al. The relationship between microvasculature in white matter hyperintensities and cognitive function. Brain Imaging Behav. 2017;11(2):503–11.

    Article  PubMed  Google Scholar 

  28. Lee JH, Cheong H, Lee SS, Lee CK, Sung YS, Huh JW, Song JA, Choe H. Perfusion assessment using Intravoxel incoherent motion-based analysis of diffusion-weighted magnetic resonance imaging: validation through phantom experiments. Investig Radiol. 2016;51(8):520–8.

    Article  Google Scholar 

  29. Iima M, Le Bihan D. Clinical Intravoxel incoherent motion and diffusion MR imaging: past, present, and future. Radiology. 2016;278(1):13–32.

    Article  PubMed  Google Scholar 

  30. van Bussel FC, Backes WH, Hofman PA, van Oostenbrugge RJ, Kessels AG, van Boxtel MP, Schram MT, Stehouwer CD, Wildberger JE, Jansen JF. On the interplay of microvasculature, parenchyma, and memory in type 2 diabetes. Diabetes Care. 2015;38(5):876–82.

    Article  PubMed  CAS  Google Scholar 

  31. Lemke A, Stieltjes B, Schad LR, Laun FB. Toward an optimal distribution of b values for intravoxel incoherent motion imaging. Magn Reson Imaging. 2011;29(6):766–76.

    Article  PubMed  Google Scholar 

  32. Celik A. Effect of imaging parameters on the accuracy of apparent diffusion coefficient and optimization strategies. Diagn Interv Radiol. 2016;22(1):101–7.

    Article  PubMed  Google Scholar 

  33. Rydhog AS, van Osch MJ, Lindgren E, Nilsson M, Latt J, Stahlberg F, Wirestam R, Knutsson L. Intravoxel incoherent motion (IVIM) imaging at different magnetic field strengths: what is feasible? Magn Reson Imaging. 2014;32(10):1247–58.

    Article  PubMed  Google Scholar 

  34. Cui Y, Dyvorne H, Besa C, Cooper N, Taouli B. IVIM Diffusion-weighted Imaging of the Liver at 3.0T: Comparison with 1.5T. Eur J Radiol Open. 2015;2:123–8.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Wetter A, Nensa F, Lipponer C, Guberina N, Olbricht T, Schenck M, Schlosser TW, Gratz M, Lauenstein TC. High and ultra-high b-value diffusion-weighted imaging in prostate cancer: a quantitative analysis. Acta Radiol. 2015;56(8):1009–15.

    Article  PubMed  Google Scholar 

  36. Jambor I, Merisaari H, Aronen HJ, Jarvinen J, Saunavaara J, Kauko T, Borra R, Pesola M. Optimization of b-value distribution for biexponential diffusion-weighted MR imaging of normal prostate. J Magn Reson Imaging. 2014;39(5):1213–22.

    Article  PubMed  Google Scholar 

  37. Togao O, Hiwatashi A, Yamashita K, Kikuchi K, Mizoguchi M, Yoshimoto K, Suzuki SO, Iwaki T, Obara M, Van Cauteren M, et al. Differentiation of high-grade and low-grade diffuse gliomas by intravoxel incoherent motion MR imaging. Neuro-Oncology. 2016;18(1):132–41.

    Article  PubMed  Google Scholar 

  38. Lin Y, Li J, Zhang Z, Xu Q, Zhou Z, Zhang Z, Zhang Y, Zhang Z. Comparison of Intravoxel incoherent motion diffusion-weighted MR imaging and arterial spin labeling MR imaging in Gliomas. Biomed Res Int. 2015;2015:234245.

    PubMed  PubMed Central  Google Scholar 

  39. Bai Y, Lin Y, Tian J, Shi D, Cheng J, Haacke EM, Hong X, Ma B, Zhou J, Wang M. Grading of Gliomas by using Monoexponential, Biexponential, and stretched exponential diffusion-weighted MR imaging and diffusion kurtosis MR imaging. Radiology. 2015;278(2):496–504.

    Article  PubMed  Google Scholar 

  40. Federau C, Meuli R, O'Brien K, Maeder P, Hagmann P. Perfusion measurement in brain gliomas with intravoxel incoherent motion MRI. Am J Neuroradiol. 2014;35(2):256–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Cohen AD, Schieke MC, Hohenwalter MD, Schmainda KM, Schmainda KM. The effect of low b-values on the intravoxel incoherent motion derived pseudodiffusion parameter in liver. Magn Reson Med. 2015;73(1):306–11.

    Article  PubMed  Google Scholar 

  42. Jalnefjord O, Andersson M, Montelius M, Starck G, Elf AK, Johanson V, Svensson J, Ljungberg M. Comparison of methods for estimation of the intravoxel incoherent motion (IVIM) diffusion coefficient (D) and perfusion fraction (f). MAGMA. 2018;31(6):715–23.

    Article  CAS  PubMed  Google Scholar 

  43. Li GF, Duan SJ, Yan LF, Wang W, Jing Y, Yan WQ, Sun Q, Wang SM, Nan HY, Xu TY, et al. Intravoxel incoherent motion diffusion-weighted MR imaging parameters predict pathological classification in thymic epithelial tumors. Oncotarget. 2017;8(27):44579–92.

    PubMed  PubMed Central  Google Scholar 

  44. Wei Y, Gao F, Wang M, Huang Z, Tang H, Li J, Wang Y, Zhang T, Wei X, Zheng D, et al. Intravoxel incoherent motion diffusion-weighted imaging for assessment of histologic grade of hepatocellular carcinoma: comparison of three methods for positioning region of interest. Eur Radiol. 2019;29(2):535–44.

    Article  PubMed  Google Scholar 

  45. 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(3):874–81.

    Article  PubMed  Google Scholar 

  46. Le Bihan D. What can we see with IVIM MRI? NeuroImage. 2019;187:56–67.

    Article  PubMed  Google Scholar 

  47. Malagi AV, Das CJ, Khare K, Calamante F, Mehndiratta A. Effect of combination and number of b values in IVIM analysis with post-processing methodology: simulation and clinical study. MAGMA. 2019;32(5):519–27.

    Article  PubMed  Google Scholar 

  48. Bisdas S, Klose U. IVIM analysis of brain tumors: an investigation of the relaxation effects of CSF, blood, and tumor tissue on the estimated perfusion fraction. MAGMA. 2015;28(4):377–83.

    Article  PubMed  Google Scholar 

  49. Dyvorne H, Jajamovich G, Kakite S, Kuehn B, Taouli B. Intravoxel incoherent motion diffusion imaging of the liver: optimal b-value subsampling and impact on parameter precision and reproducibility. Eur J Radiol. 2014;83(12):2109–13.

    Article  PubMed  PubMed Central  Google Scholar 

  50. While PT. A comparative simulation study of bayesian fitting approaches to intravoxel incoherent motion modeling in diffusion-weighted MRI. Magn Reson Med. 2017;78(6):2373–87.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Dr. Tian-Yong Xu in GE Healthcare China for providing technical support regarding the appropriate applied multiple b-value DWI sequence.

Funding

This work was supported by the Science and Technology Innovation Development Foundation of Tangdu Hospital (No. 2015JCYJ010 and 2017LCYJ004) in the data collection and manuscript revision.

Author information

Author notes

  1. Yu-Chuan Hu and Lin-Feng Yan contributed equally to this work.

Authors and Affiliations

  1. Department of Radiology and Functional and Molecular Imaging Key Lab of Shaanxi Province, Tangdu Hospital, Air Force Medical University (Fourth Military Medical University), Xi’an, 710038, Shaanxi, People’s Republic of China

    Yu-Chuan Hu, Lin-Feng Yan, Yu Han, Shi-Jun Duan, Qian Sun, Gang-Feng Li, Wen Wang & Guang-Bin Cui

  2. MR Research China, GE Healthcare China, Beijing, 100176, China

    Xiao-Cheng Wei & Dan-Dan Zheng

Authors
  1. Yu-Chuan Hu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lin-Feng Yan
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yu Han
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Shi-Jun Duan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Qian Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Gang-Feng Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Wen Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xiao-Cheng Wei
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Dan-Dan Zheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Guang-Bin Cui
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

GBC conceived the study. YCH, LFY, XCW and DDZ participated in the study design. YCH, LFY, YH, SJD, QS, GFL and WW performed the data acquisition. YCH, LFY and WW participated in the statistical analyses. All authors participated in the data interpretation. YCH drafted the first version of the report. All authors revised and approved the final draft of the report.

Corresponding author

Correspondence to Guang-Bin Cui.

Ethics declarations

Ethics approval and consent to participate

This study was conducted in accordance with the ethical guidelines of the Declaration of Helsinki (version 2002). The brain MR imaging of healthy volunteers was approved by the Medical Ethics Committee of Tangdu Hospital of the Fourth Military Medical University, and informed written consent were obtained from all participants.

Consent for publication

The manuscript is approved by all participants for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hu, YC., Yan, LF., Han, Y. et al. Can the low and high b-value distribution influence the pseudodiffusion parameter derived from IVIM DWI in normal brain?. BMC Med Imaging 20, 14 (2020). https://doi.org/10.1186/s12880-020-0419-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12880-020-0419-0

Keywords

BMC Medical Imaging

ISSN: 1471-2342

Contact us