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

Mono, bi- and tri-exponential diffusion MRI modelling for renal solid masses and comparison with histopathological findings

Abstract

Purpose

To compare diffusion tensor imaging (DTI), intravoxel incoherent motion (IVIM), and tri-exponential models of the diffusion magnetic resonance imaging (MRI) signal for the characterization of renal lesions in relationship to histopathological findings.

Methods

Sixteen patients planned to undergo nephrectomy for kidney tumour were scanned before surgery at 3 T magnetic resonance imaging (MRI), with T2-weighted imaging, DTI and diffusion weighted imaging (DWI) using ten b-values. DTI parameters (mean diffusivity [MD] and fractional anisotropy [FA]) were obtained by iterative weighted linear least squared fitting of the DTI data and bi-, and tri-exponential fit parameters (Dbi, fstar,and Dtri, ffast,finterm) using a nonlinear fit of the multiple b-value DWI data. Average parameters were calculated for regions of interest, selecting the lesions and healthy kidney tissue. Tumour type and specificities were determined after surgery by histological examination. Mean parameter values of healthy tissue and solid lesions were compared using a Wilcoxon-signed ranked test and MANOVA.

Results

Thirteen solid lesions (nine clear cell carcinomas, two papillary renal cell carcinoma, one haemangioma and one oncocytoma) and four cysts were included. The mean MD of solid lesions are significantly (p < 0.05) lower than healthy cortex and medulla, (1.94 ± 0.32*10− 3 mm2/s versus 2.16 ± 0.12*10− 3 mm2/s and 2.21 ± 0.14*10− 3 mm2/s, respectively) whereas ffast is significantly higher (7.30 ± 3.29% versus 4.14 ± 1.92% and 4.57 ± 1.74%) and finterm is significantly lower (18.7 ± 5.02% versus 28.8 ± 5.09% and 26.4 ± 6.65%). Diffusion coefficients were high (≥2.0*10− 3 mm2/s for MD, 1.90*10− 3 mm2/s for Dbi and 1.6*10− 3 mm2/s for Dtri) in cc-RCCs with cystic structures and/or haemorrhaging and low (≤1.80*10− 3 mm2/s for MD, 1.40*10− 3 mm2/s for Dbi and 1.05*10− 3 mm2/s for Dtri) in tumours with necrosis or sarcomatoid differentiation.

Conclusion

Parameters derived from a two- or three-component fit of the diffusion signal are sensitive to histopathological features of kidney lesions.

Background

As a result of the increased use of abdominal imaging, more (asymptomatic) small (≤ 4 cm) renal masses are incidentally discovered. In a series of 173 patients only 58% of kidney tumours < 4 cm were malignant, whereas all kidney tumours > 7 cm were [1]. Hence, a substantial amount of incidentally discovered renal masses is not malignant [2,3,4]. The management of renal lesions includes radical or partial nephrectomy, minimal invasive ablative techniques or active surveillance. Because of concern for chronic kidney disease, nephron sparing surgery is preferred [5, 6] but more importantly unnecessary surgery should be avoided. One way to realize this is by distinguishing between lesion types and reliably diagnosing benign tumour types, such as oncocytoma, prior to treatment [7]. However, with currently available clinical imaging modalities, benign renal masses are indistinguishable from malignant renal masses [4, 8].

Many magnetic resonance imaging (MRI) techniques have been explored as methods to differentiate between benign and malignant renal lesions or between renal cell carcinoma (RCC) subtypes [4, 9,10,11]. One promising technique is diffusion-weighted imaging (DWI), which allows quantification of water motion in tissues without administration of exogenous contrast materials [12,13,14,15,16]. The apparent diffusion coefficient (ADC), derived from a mono-exponential model, is believed to reflect tissue cellularity as a higher tissue density will amount to more restricted diffusion, hence a lower diffusion value. However, ADC values for different subtypes may overlap, making determination of cut-off values to distinguish between benign and malignant solid renal masses problematic [17].

More complex models of diffusion, such as the diffusion tensor model (DTI) and the intravoxel incoherent motion model (IVIM), allow deriving additional information. DTI-derived parameters fractional anisotropy (FA) and mean diffusivity (MD) have been correlated with histological parameters such as cell density and nuclear grade [18]. The IVIM model is a bi-exponential model that includes molecular diffusion and microcirculation of blood in the capillary network (‘pseudodiffusion’) [19]. A combination of pseudodiffusion fraction fbi and the perfusion-free diffusion coefficient Dbi from IVIM model is able to differentiate between renal tumour types [20, 21]. Recently, the IVIM model was expanded to a three-component model by adding an additional component that accounts for intermediately fast water motion in the kidney [22, 23]. The aim of this study is to compare parameters obtained from DTI, intravoxel incoherent motion (IVIM), and tri-exponential models of the diffusion signal of kidney lesions, for the characterization of renal lesions. Because tumours are usually not uniform and may consist of several areas with different structural patterns, we compare diffusion parameters with histopathological results.

Methods

Subjects

Approval of our institution’s ethical committee was obtained for this prospective study and all subjects provided written informed consent. From March 2016 to May 2017, sixteen patients (11 male, age 65 (range 50–76) years old, 5 female, age 60 (range 48–72), total group: age 64 (range 48–76) years old) who had suspected kidney tumours and were planned to undergo radical or partial nephrectomy based on standard clinical diagnostic criteria were included. After including the first five consecutive patients, patients were also selected on tumour size (≤ 4 cm on radiologic examination) in order to increase chance of including benign solid lesions. After surgical resection of the tumour, kidney tumour type was determined according to the WHO classification of tumours of the urinary system [24] by histopathological examination of 2-μm-thick sections of formalin-fixed and paraffin-embedded tumour tissue blocks using haematoxylin-eosin (HE) staining.

Scans

A T2 weighted sequence was performed for anatomical reference, followed by a DTI sequence (b = 0, 100 and 300 s/mm2 in 15 gradient direction) and a DWI sequence including ten b-values (b = 0, 10, 25, 40, 75, 100, 200, 300, 500 and 700 s/mm2 in six gradient directions) on a 3 T MRI clinical scanner (Philips, Ingenia, Philips Healthcare, Best, The Netherlands), see Table 1 for MRI acquisition details.

Table 1 MRI acquisition details of the T2-TSE, DTI and IVIM protocols

Data processing

To enable accurate parameter fitting all scans were corrected for (breathing) motion before further processing. Due to differences in motion between the right and left kidneys, they were cropped and processed as separate data sets, as described previously [22]. All pre-processing was performed using diffusion imaging analysis package DTItools [github.com/mfroeling/DTITools] [25] and image registration toolbox Elastix [http://elastix.isi.uu.nl/] [26]. First, T2 scans were processed to correct for slice by slice misalignment due to acquisition in multiple breath-holds using a rigid 2D registration algorithm after being resampled to 2 mm isotropic using a single interpolation method. Finally, all DWI data was corrected for breathing motion, by registering them to the unweighted volume using a rigid 2D b-spline registration algorithm after which the DWI data was registered to the reference T2 scan using a 3D affine registration algorithm [22].

Parameter maps

From the DTI data the FA and MD were calculated using an iterative weighted linear least squares algorithm with outlier rejection using ExploreDTI [27]. From the IVIM data, bi- and tri-exponential diffusion decay parameters, i.e. the mean diffusion (Dbi for bi-exponential and Dtri for tri-exponential fitting), and the signal fraction attributed to pseudo-diffusion (fstar for bi-exponential and ffast and finterm for tri-exponential fitting), were obtained by fitting a two and three-component model to the multiple b-value DWI data, as described previously [22, 23]. To make a comparison between the DTI and IVIM data, the mean diffusion from a mono-exponential fit, Dmono was also obtained.

Regions of interest (ROIs) to segment the tumour volumes were manually defined on the combined T2 and DWI data by the principal researcher (S.v.B., 4 years of experience) in agreement with an experienced radiologists (C.K., 12 years of experience) using image segmentation toolbox ITK snap [28]. ROIs were placed inside the tumour, rather than following the contour, to limit the contribution of the signal from other tissue types due to partial volume effect or imprecise image registration. For comparison of tumour tissue with healthy kidney parenchyma, the cortex and medulla in the healthy contralateral kidneys were segmented using an automated algorithm as in [22]. The mean and standard deviation of the diffusion parameters were obtained for healthy cortex and medulla and lesion ROIs. Parameter maps MD, Dbi, Dtri, fstar, ffast and finterm were obtained for visual comparison with T2 and (if available) photographs of the gross appearance of the resected kidney tumours before histological examination.

Statistical analysis

All statistical tests were performed using SPSS (version 23.0. Armonk, NY: IBM Corp.). Healthy cortex and medulla were compared with all solid lesions using a Wilcoxon Signed Ranks test. The means of parameters MD, FA, Dbi, Dtri, fstar, ffast and finterm in healthy cortex and medulla, different types of RCCs, cysts and benign solid lesions were compared using multivariate analysis of variance (MANOVA). Bonferroni correction was applied, and a p-value < 0.05 was considered significant for all statistical tests.

Results

Subjects and scans

All patients were successfully scanned. One scan was removed before processing due to artefacts resulting from a lower-back implantation. After visual inspection following processing two other scans were removed, due to poor motion correction results and an error in the data. In the remaining thirteen scans, thirteen solid lesions (average size of maximum diameter, determined by histopathological examination 3.85 cm, range 0.8–7.5 cm) and five fluid-filled cysts were found. Examples of the raw acquired data and the data after motion correction and image registration together with the ROI placement in one lesion are shown in Fig. 1.

Fig. 1
figure 1

Anatomical reference T2 image before (a) and after (b) motion correction and image registration diffusion weighted images before (d) and after (e) processing, ROI selecting tumor tissue in T2 (c) and same ROI projected on DWI (f)

Histological examination

Of the solid lesions, eleven were considered to be malignant (nine clear cell RCCs (cc-RCCs), two papillary cell RCCs (p-RCCs) and two were considered to be benign (one haemangioma of the kidney capsule and one oncocytoma). Of the nine cc-RCCs, only one had a homogeneous microstructural pattern consisting mainly of clear cells. Others had considerable amounts of necrosis, sarcomatoid differentiation, haemorrhaging or deviating growth patterns such as papillary, tubular or cystic growth. In Fig. 2, histopathological features of several kidney tumour tissues are displayed.

Fig. 2
figure 2

Histopathological features of kidney (tumor) tissue types, sections are stained with haematoxylin-eosin. a Normal kidney tissue, with kidney cortex including glomeruli on the right and medulla including tubular structures to the left. Magnification: 100x, scale bar represents 100 μm, (b) Papillary RCC which presents with papillae and the presence of macrophages in nuclei. Magnification: 200x, (c) Clear cell RCC with to the right a necrotic area. Magnification: 400x, (d) Clear cell RCC with to the right micro cystic structures. Magnification: 200x,, (e) Sarcomatoid differentiation in RCC, displaying elongated, spindle-shaped cells, high cellularity and cellular atypia. Magnification: 400x, (f) Oncocytoma, displaying granular eosinophilic cytoplasm. Magnification: 400x scale bar represents 50 μm in (b-f)

Parameter maps

Figure 3 shows MRI data (T2 and DWI b = 0) parameters maps and (where available) photographs of the gross appearance of the resected kidney tumours for: a contralateral unaffected kidney (A-D), a case of cc-RCC (E-H) and the oncocytoma (I-L). In the Additional file 1: Figure S1 the gross appearance, MRI data and parameter maps of all tumour types are shown. The fractions of the diffusion components are shown as merged f-maps where ffast, finterm and fslow are colour coded red, blue and green, respectively. For the unaffected kidney, ffast (red) was high in those areas with a high blood flow (e.g. large blood vessels) whereas finterm (blue) was high in areas with free water (e.g. the pyelum). The diffusion coefficient (Dtri) was homogeneous throughout healthy kidney parenchyma.

Fig. 3
figure 3

Diffusion-derived parameter maps of an unaffected kidney (a-d), a kidney with renal cell carcinoma (e-h), and a kidney with oncocytoma (i-l). First row (e, j): gross appearance of the whole kidney or tumour after nephrectomy, (a, f, k): anatomical reference (after processing) (T2-TSE), second row (b, g, l): the unweighted image of the diffusion scan after motion correction and image registration (DWI-b0), third row (c, h, m): a merge of the fraction maps from the tri-exponential fit red = ffast, blue = finterm, green = fslow (1- finterm - ffast), fourth row (d, i, n): diffusion coefficient from the tri-exponential fit, (Dtri)

Upon visual examination, in the maps obtained in a cc-RCC, Dtri was lower throughout the tumour, and fslow had a high contribution. For the oncocytoma, the lesion did not seem much different from normal kidney parenchyma, although fslow and Dtri appeared higher.

In the p-RCC, the merged f-maps showed a small contribution from ffast and a larger contribution from fslow. The Dtri map indicated a low diffusion coefficient. In the cc-RCC with sarcomatoid differentiation and the cc-RCC with papillary growth the photographs and the diffusion parameter maps showed a more heterogeneous make-up, indicating a more complex tumour. In the cysts, fslow had a high contribution and Dtri was high. For the haemangioma, fslow and Dtri appeared higher.

Parameter analysis

The mean and standard deviations of parameter values of grouped lesions and measurements of healthy cortex and medulla are given in Table 2. In Additional file 2: Figure S2 Dmono is plotted against MD for each lesion. In Additional file 3: Table S1, the values for Dmono are displayed together with the other diffusion coefficients MD, Dbi, and Dtri. MD and Dmono have a similar order of magnitude for each (group of) lesion, whereas Dbi and Dtri are structurally lower.

Table 2 Mean (standard differentiation) of the diffusion parameters obtained from DTI, and the bi- and tri-exponential models for each of the segmented tissue types

Solid lesions had a significantly lower MD (1.94 ± 0.32 10− 3 mm2/s) than healthy tissue (2.16 ± 0.12 10− 3 mm2/s for cortex, p = 0.019 and 2.21 ± 0.14 10− 3 mm2/s for medulla, p = 0.009) and a significantly lower (p = 0.023) Dbi (2.02 ± 0.11 10− 3 mm2/s) than healthy medulla (1.71 ± 0.43 10− 3 mm2/s), but no significant difference in Dtri. MD, Dbi and Dtri were all significantly higher (p < 0.002) in cysts than healthy tissue and other lesions (see Table 2).

Compared to healthy tissue (0.38 ± 0.09 for cortex and 0.39 ± 0.08 for medulla), FA was higher in all solid lesions (0.47 ± 0.11), RCCs (0.46 ± 0.10) and benign solid lesions (0.48 ± 0.19) but similar in cysts (0.37 ± 0.11), however these differences were not significant.

fstar did not show significant differences between healthy tissue (10.1 ± 2.58% for cortex, 9.69 ± 2.90% % for medulla), and solid lesions (11.6 ± 3.88%). finterm was significantly higher in healthy tissue (28.8 ± 5.09% for cortex and 26.4 ± 6.65% for medulla) than solid lesions (18.7 ± 5.02%, p = 0.003 for cortex, p = 0.011 for medulla) and RCCs (18.3 ± 5.35%, p = 0.000 for cortex, p = 0.009 for medulla). ffast was a significantly lower for healthy tissue (4.14 ± 1.92% for cortex, p = 0.009, and 4.57 ± 1.74% for medulla, p = 0.033) than all solid lesions (7.30 ± 3.29%), and healthy cortex had a significantly lower ffast than RCCs (7.21 ± 2.88%, p = 0.034). Cysts had a significantly lower fstar (1.88 ± 1.60%) than healthy tissue (p = 0.000 for cortex, 0.001 for medulla), RCCs (11.6 ± 3.45%, p = 0.000) and benign lesions (11.6 ± .83%, p = 0.009), a significantly lower finterm (8.31 ± 2.02%, p = 0.034) than RCCs and a significantly lower ffast (1.18 ± 1.70%) than RCCs (p = 0.001) and benign solid lesions (7.76 ± 6.75%, p = 0.030).

Parameter values separated according to tumour type are given in Table 3, together with values for individual heterogeneous RCCs. There was a large variation in all parameters among the cc-RCCs. Diffusion coefficients were high (≥ 2.0*10− 3 mm2/s for MD, 1.90*10− 3 mm2/s for Dbi and 1.6*10− 3 mm2/s for Dtri) in cc-RCCs with cystic structures (see Fig. 2d for microscopic photograph) and/or haemorrhaging and low (≤ 1.80*10− 3 mm2/s for MD, 1.40*10− 3 mm2/s for Dbi and 1.05*10− 3 mm2/s for Dtri) in tumours with necrosis (Fig. 2c) or sarcomatoid differentiation (Fig. 2e). fstar and ffast is high (16.62 and 12.43% respectively) in a tumour with extensive haemorrhaging, and low (≤ 11% for fstar and ≤ 6.5% for ffast) in tumours with cystic structures and/or necrosis. finterm is particularly high in a cc-RCC with micro-cystic structures and haemorrhaging.

Table 3 Mean (standard differentiation) of the diffusion parameters obtained from DTI, and the bi- and tri-exponential models for each type of solid lesion

The haemangioma had the highest diffusivity in all models (2.38*10− 3 mm2/s for MD, 2.45*10− 3 mm2/s for Dbi and 1.99*10− 3 mm2/s for Dtri). Additionally, the haemangioma had the highest FA (0.61). The oncocytoma (Fig. 2f) had the highest fstar (17.1%) and ffast (12.5%), although the cc-RCC with papillary growth had a comparable fstar (16.5%). The cysts and haemangioma had a low fstar (7.95 and 6.06%, respectively) and ffast (2.83 and 2.99% respectively). finterm was in the range of 19–23% except for p-RCC (16%) and cysts (32%).

In Fig. 4, each lesion is represented by a dot in the scatter plot, plotting DTI parameters MD versus FA (Fig. 3a), IVIM parameters Dbi versus fstar (Fig. 3b), three-component parameters Dtri versus finterm (Fig. 3c) and Dtri versus ffast (Fig. 3d). Cysts were recognizable by a high MD, Dbi and Dtri, but they were more grouped in Dtri versus finterm/ffast due to a consistently low finterm and ffast. In Dbi versus fstar and Dtri versus ffast several general groupings were identifiable; the cysts are located in the lower right corner, the cc-RCCs in the middle and p-RCC (Fig. 2b for microscopic photograph) in the middle to the left. The oncocytoma is located above the cc-RCCs and the haemangioma between cc-RCCs and the cysts. However, one of the cc-RCCs is closely located to the oncocytoma due to a high value of fstar and ffast. The p-RCCs and cc-RCC with sarcomatoid differentiation seemed to have a slightly lower Dbi and Dtri than most cc-RCCs, but one cc-RCC has a much lower diffusivity than all other lesions.

Fig. 4
figure 4

scatter plots, each point represents a single lesion. a: FA vs. MD (DTI fit), b: Dbi vs. fstar (bi-exponential fit). c: Dtri vs. finterm (tri-exponential fit). d: Dtri vs. ffast (tri-exponential fit). RCC = renal cell carcinoma, cc-RCC = clear cell renal cell carcinoma, p-RCC = papillary renal cell carcinoma

Figures 5 and 6 show box plots for MD (Fig. 5a), Dbi (Fig. 5b), Dtri (Fig. 5c), FA (Fig. 6a), fstar, (Fig. 6b), ffast (Fig. 6c) and finterm (Fig. 6d). The finterm and ffast show more pronounced differences and less overlap between individual lesions than fstar. and compared to MD, both Dbi and Dtri gave more pronounced differences and less overlap between tissue types.

Fig. 5
figure 5

Boxplots for diffusion coefficients, (a) MD, (b) Dbi, and (c) Dtri, for individual lesions, cysts, RCCs and healthy cortex and medulla. Points represent individual lesions, error bars groups of lesions. RCC = renal cell carcinoma, cc-RCC = clear cell renal cell carcinoma, p-RCC = papillary renal cell carcinoma

Fig. 6
figure 6

Boxplots for (a) FA, (b) fstar, (c) ffast, and (d) and finterm or individual lesions, cysts, RCCs and healthy cortex and medulla. Points represent individual lesions, error bars groups of lesions. RCC = renal cell carcinoma, cc-RCC = clear cell renal cell carcinoma, p-RCC = papillary renal cell carcinoma

Discussion

In this study, we have compared the DTI, IVIM and three-compartment models of the diffusion signal for the characterization of renal lesions. In the study population of sixteen patients who received radical or partial nephrectomy, two lesions were found benign upon histological examination. This indicates that resection was unnecessary and could have been prevented if these lesions were identified as non-malignant prior to surgery. Based on our results the haemangioma could potentially have been identified using diffusion-derived parameters, notably, a high Dbi and Dtri and a low fstar and ffast.

Previous studies applying DWI to evaluate kidney tumour type concluded that kidney lesions generally have a lower ADC than normal kidney tissue [13, 17, 29] and some studies also found a higher mean ADC value in benign lesions than in malignant lesions [15, 30, 31]. Studies typically find a higher ADC in oncocytomas than in RCCs [12,13,14,15, 17, 30,31,32] and therefore, ADC is likely to be a valuable parameter in evaluating tumour type. However, in itself it is not sufficient as a clinical index due to overlap in values between tumour tissue types [16, 17, 33].

In line with previous studies MD in this study was significantly lower in solid lesions than in healthy tissue, and RCCs had a lower MD than benign lesions and cysts. Differences in diffusivity between tissue types are usually attributed to differences in tissue cellularity, a higher cellularity will result in more restricted diffusion and therefore, more aggressive lesions are expected to present with a lower diffusivity [15, 21]. Our study also showed similar results: cysts had the highest diffusion coefficients whereas lesions with a high degree of necrosis had low diffusivity, which we have assigned to an increase in with a diffusion-restricting elements, such as macromolecules, and disorganisation.

In the two-component IVIM model, a fast moving component is separated from the diffusion signal, resulting in a lower diffusion coefficient Dbi. The fraction fstar of the diffusion signal that is attributed to fast moving water has been correlated to renal tumour vascularity [21]. As in this study, a lower MD and Dbi for p-RCC than cc-RCC and similar values for fstar in p-RCCs were found previously [21]. A tri-exponential fit of the diffusion signal in the kidney was previously shown to provide additional information on structures associated with pseudodiffusion by separating the fast from intermediate pseudodiffusion resulting in signal fractions finterm associated with a diffusion rate in the order of magnitude of free water and ffast associated with perfusion [22, 23], and comparable values to this study for parameters derived from DTI, two- and three-compartment fitting for healthy cortex and medulla were reported [22].

MD was the only diffusion coefficient to be significantly different in healthy tissue from RCCs but within the tumour type groups the MD range was wider whereas the differences between tumour types were more pronounced in Dbi and Dtri,. Because diffusion coefficients Dbi and Dtri exclude the fraction of the diffusion signal that is attributed to fast water movement, they were lower but more precise than MD. Therefore, Dbi and Dtri better reflect tissue diffusion, making both these parameters more specific for tissue cellularity. However, none of the diffusion coefficients could be used to reliably distinguish between lesion types, as diffusion coefficients overlap. For a better contrast between lesion types, tissue cellularity (MD, Dbi or Dtri) can be combined with a measure for vascularisation (fstar or ffast). For example, a combination of Dbi and fstar was previously shown to discriminate between renal tumour subtypes [21].

There was a large variation in diffusion parameter values between individual cc-RCC tumours due to tumour heterogeneity and differences in microstructural make-up. In addition, drawing statistical inferences from this study is limited due to the small number of cases and the sensitivity of the analysis to outliers. Therefore we have also analysed the results from individual lesions, associating histopathological characteristics with the interpretations of diffusion parameters outlined above. In this analysis, low diffusion coefficients were found in cc-RCCs with a high tissue density (due to extensive necrosis or sarcomatoid differentiation) whereas high diffusion coefficients were found in cc-RCCs with cystic structures. Additionally, tumours with a high perfusion rate are characterized by a high value of fstar or ffast whereas tumours with a low perfusion rate, such as the haemangioma, cysts and the cc-RCCs with cystic structures have lower fstar or ffast. Hence, these diffusion parameters seem to be indicative of histopathological features of kidney tumours.

Although both fstar and ffast seem to correlate to vascularisation, only the three-component parameters ffast and finterm show significant differences between different tissues whereas two-component parameter fstar does not, showing that the tri-exponential model provides additional information over the bi-exponential model.

Because of the limited amount of cases it is impossible to formulate conclusions regarding the characterisation of kidney tumour type. Therefore, we have also analysed individual lesions relating our findings to histopathological details. The initial results from this analysis indicate that diffusion parameters are sensitive to histopathological features of kidney lesions, which is a first step towards non-invasive characterisation of these lesions prior to treatment. An improvement to this study would be to spatially correlate histopathology to parameter measurements and maps [34]. This would result in more specific validation of diffusion parameters, confirming the correlation between diffusion parameters to histopathological features of tissues, such as cellularity, perfusion and cystic structures. In addition, to establish what parameter values should be used to confidently distinguish between benign and malignant lesions and draw statistically significant conclusions, a larger study population should be included. However, since kidney tumour type is unknown before histopathological evaluation, researchers have no control over which tumour types are included. To increase the amount of benign kidney tumours, only small (≤ 4 cm) renal masses (about 40% benign [1]) can be included. Additionally, this study shows that parameters derived from the DTI sequence (MD and FA) do not provide additional information over parameters derived from a multiple b-value sequence (from two- and three- component fits). Hence, the DTI sequence can be omitted, decreasing total scanner time with one third, to about 30 min.

Conclusion

In conclusion, parameters derived from a two- or three-component fit of the diffusion signal are sensitive to histopathological features of kidney lesions.

Abbreviations

ADC:

Apparent diffusion coefficient

cc-RCC:

Clear cell renal cell carcinoma

DTI:

Diffusion tensor imaging

DWI:

Diffusion weighted imaging

FA:

Fractional anisotropy

HE:

Haematoxylin-eosin

IVIM:

Intravoxel incoherent motion

MD:

Mean diffusivity

p-RCC:

Papillary cell renal cell carcinoma

RCC:

Renal cell carcinoma

ROI:

Region of interest

WHO:

World Health Organisation

References

  1. Duchene DA, Lotan Y, Cadeddu JA, Sagalowsky AI, Koeneman KS. Histopathology of surgically managed renal tumors: analysis of a contemporary series. Urology. 2003;62(5):827–30.

    Article  Google Scholar 

  2. Thompson RH, Kurta JM, Kaag M, Tickoo SK, Kundu S, Katz D, Nogueira L, Reuter VE, Russo P. Tumor size is associated with malignant potential in renal cell carcinoma cases. J Urol. 2009;181(5):2033–6.

    Article  Google Scholar 

  3. Rothman J, Egleston B, Wong YN, Iffrig K, Lebovitch S, Uzzo RG. Histopathological characteristics of localized renal cell carcinoma correlate with tumor size: a SEER analysis. J Urol. 2009;181(1):29–33 discussion 33-4.

    Article  Google Scholar 

  4. Kang SK, Huang WC, Pandharipande PV, Chandarana H. Solid renal masses: what the numbers tell us. AJR Am J Roentgenol. 2014;202(6):1196–206.

    Article  Google Scholar 

  5. Li L, Lau WL, Rhee CM, Harley K, Kovesdy CP, Sim JJ, Jacobsen S, Chang A, Landman J, Kalantar-Zadeh K. Risk of chronic kidney disease after cancer nephrectomy. Nat Rev Nephr. 2014;10:135–45.

    Article  Google Scholar 

  6. Van Poppel H, Becker F, Cadeddu JA, Gill IS, Janetschek G, Jewett MA, Laguna MP, Marberger M, Montorsi F, Polascik TJ, Ukimura O, Zhu G. Treatment of localised renal cell carcinoma. Eur Urol. 2011;60(4):662–72.

    Article  Google Scholar 

  7. Kutikov A, Fossett LK, Ramchandani P, Tomaszewski JE, Siegelman ES, Banner MP, Van Arsdalen KN, Wein AJ, Malkowicz SB. Incidence of benign pathologic findings at partial nephrectomy for solitary renal mass presumed to be renal cell carcinoma on preoperative imaging. Urology. 2006;68(4):737–40.

    Article  Google Scholar 

  8. Prasad SR, Surabhi VR, Menias CO, Raut AA, Chintapalli KN. Benign renal neoplasms in adults: cross-sectional imaging findings. AJR Am J Roentgenol. 2008;190(1):158–64.

    Article  Google Scholar 

  9. Ramamurthy NK, Moosavi B, McInnes MD, Flood TA, Schieda N. Multiparametric MRI of solid renal masses: pearls and pitfalls. Clin Radiol. 2015;70(3):304–16.

    Article  CAS  Google Scholar 

  10. Pedrosa I, Alsop DC, Rofsky NM. Magnetic resonance imaging as a biomarker in renal cell carcinoma. Cancer. 2009;115(10 Suppl):2334–45.

    Article  Google Scholar 

  11. Gilet AG, Kang SK, Kim D, Chandarana H. Advanced renal mass imaging: diffusion and perfusion MRI. Curr Urol Rep. 2012;13(1):93–8.

    Article  Google Scholar 

  12. Agnello F, Roy C, Bazille G, Galia M, Midiri M, Charles T, Lang H. Small solid renal masses: characterization by diffusion-weighted MRI at 3 T. Clin Radiol. 2013;68(6):e301–8.

    Article  CAS  Google Scholar 

  13. Mytsyk Y, Dutka I, Borys Y, Komnatska I, Shatynska-Mytsyk I, Farooqi AA, Gazdikova K, Caprnda M, Rodrigo L, Kruzliak P. Renal cell carcinoma: applicability of the apparent coefficient of the diffusion-weighted estimated by MRI for improving their differential diagnosis, histologic subtyping, and differentiation grade. Int Urol Nephrol. 2017;49:215–24.

    Article  CAS  Google Scholar 

  14. Inci E, Hocaoglu E, Aydin S, Cimilli T. Diffusion-weighted magnetic resonance imaging in evaluation of primary solid and cystic renal masses using the Bosniak classification. Eur J Radiol. 2012;81(5):815–20.

    Article  Google Scholar 

  15. Taouli B, Thakur RK, Mannelli L, Babb JS, Kim S, Hecht EM, Lee VS, Israel GM. Renal lesions: characterization with diffusion-weighted imaging versus contrast-enhanced MR imaging. Radiology. 2009;251(2):398–407.

    Article  Google Scholar 

  16. Lassel EA, Rao R, Schwenke C, Schoenberg SO, Michaely HJ. Diffusion-weighted imaging of focal renal lesions: a meta-analysis. Eur Radiol. 2014;24(1):241–9.

    Article  CAS  Google Scholar 

  17. Sevcenco S, Heinz-Peer G, Ponhold L, Javor D, Kuehhas FE, Klingler HC, Remzi M, Weibl P, Shariat SF, Baltzer PA. Utility and limitations of 3-tesla diffusion-weighted magnetic resonance imaging for differentiation of renal tumors. Eur J Radiol. 2014;83(6):909–13.

    Article  CAS  Google Scholar 

  18. Feng Q, Fang W, Sun XP, Sun SH, Zhang RM, Ma ZJ. Renal clear cell carcinoma: diffusion tensor imaging diagnostic accuracy and correlations with clinical and histopathological factors. Clin Radiol. 2017;72(7):560–4.

    Article  CAS  Google Scholar 

  19. Le Bihan D, Breton E, Lallemand D, Aubin M-L, Vignaud J, Laval-Jeantet M. Separation of diffusion and perfusion in Intravoxel incoherent Moting MR imaging. Radiology. 1988;168:9.

    Article  Google Scholar 

  20. Rheinheimer S, Stieltjes B, Schneider F, Simon D, Pahernik S, Kauczor HU, Hallscheidt P. Investigation of renal lesions by diffusion-weighted magnetic resonance imaging applying intravoxel incoherent motion-derived parameters--initial experience. Eur J Radiol. 2012;81(3):e310–6.

    Article  CAS  Google Scholar 

  21. Chandarana H, Kang SK, Wong S, Rusinek H, Zhang JL, Arizono S, Huang WC, Melamed J, Babb JS, Suan EF, Lee VS, Sigmund EE. Diffusion-weighted Intravoxel incoherent motion imaging of renal tumors with histopathologic correlation. Investig Radiol. 2012;47(12):8.

    Article  Google Scholar 

  22. van Baalen S, Leemans A, Dik P, Lilien MR, Ten Haken B, Froeling M. Intravoxel incoherent motion modeling in the kidneys: comparison of mono-, bi-and triexponential fit. J Magn Reson Imaging. 2017;46(1):228–39.

    Article  Google Scholar 

  23. van der Bel R, Gurney-Champion OJ, Froeling M, Stroes ESG, Nederveen AJ, Krediet CTP. A tri-exponential model for intravoxel incoherent motion analysis of the human kidney: in silico and during pharmacological renal perfusion modulation. Eur J Radiol. 2017;91:168–74.

    Article  Google Scholar 

  24. H. Moch, A.L. Cubilla, P.A. Humphrey, V.E. Reuter, T.M. Ulbright, WHO classification of tumours of the urinary system and male genital organs, International Agency for Research on Cancer 2016.

  25. Froeling M, Nederveen AJ, Heijtel DFR, Lataster A, Bos C, Nicolay K, Maas M, Drost MR, Strijkers GJ. Diffusion-tensor MRI reveals the complex muscle architecture of the human forearm. J Magn Reson Imaging. 2012;36(1):237–48.

    Article  Google Scholar 

  26. Klein S, Staring MM, Murphy K, Viergever MA, Pluim JPW. Elastix: a toolbox for intensity-based medical image registration. IEEE Trans Med Imaging. 2010;29, 9(1).

    Article  Google Scholar 

  27. Leemans A, Jeurissen B, Sijbers J, Jones DK. ExploreDTI: a graphical toolbox for processing, analyzing and visualizing diffusion MRI data. In: Proc 17th Annual Meeting ISMRM, Honolulu; 2009.

  28. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, Gerig G. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. NeuroImage. 2006;31(3):1116–28.

    Article  Google Scholar 

  29. Cova M, Squillaci E, Stacul F, Manenti G, Gava S, Simonetti G, Pozzi-Mucelli R. Diffusion-weighted MRI in the evaluation of renal lesions: preliminary results. Br J Radiol. 2004;77(922):851–7.

    Article  CAS  Google Scholar 

  30. Doganay S, Kocakoc E, Cicekci M, Aglamis S, Akpolat N, Orhan I. Ability and utility of diffusion-weighted MRI with different b values in the evaluation of benign and malignant renal lesions. Clin Radiol. 2011;66(5):420–5.

    Article  CAS  Google Scholar 

  31. Sobh DM, El Hawary GESM, Abou El Ghar M, El-Diasty TAEM, El-Sayed Settein M, ElShaer S, Tantawy MSE. Role of diffusion weighted MR imaging in characterization of focal kidney and upper urinary tract lesions. Egypt J Radiol Nucl Med. 2016;47(4):1689–700.

    Article  Google Scholar 

  32. Galmiche C, Bernhard JC, Yacoub M, Ravaud A, Grenier N, Cornelis F, Multiparametric I. MRI useful for differentiating Oncocytomas from Chromophobe renal cell carcinomas? AJR Am J Roentgenol. 2017;208(2):343–50.

    Article  Google Scholar 

  33. Liu JH, Tian SF, Ju Y, Li Y, Chen AL, Chen LH, Liu AL. Apparent diffusion coefficient measurement by diffusion weighted magnetic resonance imaging is a useful tool in differentiating renal tumors. BMC Cancer. 2015;15:292.

    Article  Google Scholar 

  34. Littooij AS. Apparent diffusion coefficient as it relates to histopathology findings in post-chemotherapy nephroblastoma: a feasibility study. Pediatric Radiology Epub ahead of print. 2017;7.

Download references

Acknowledgements

We are grateful to all the patients who have participated in this study. We thank Casper Jansen for histopathological support, Aletta Goolkate-Geerlig for technical support.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not for-profit sectors.

Availability of data and materials

The approval received by our local ethics institution does not allow for making data and material available.

Author information

Authors and Affiliations

Authors

Contributions

All authors listed have contributed sufficiently to the project to be included as authors, and all those who are qualified to be authors are listed in the author byline. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Sophie van Baalen.

Ethics declarations

Ethics approval and consent to participate

This study was approved by our local institutional review board and written informed consent was given by all subjects.

Consent for publication

All authors consent to publish this research.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Figure S1. Diffusion-derived parameter maps of each tumor type: an unaffected kidney (A-D), RCC (E-I), clear cell renal cell carcinoma (cc-RCC) with sarcomatoid differentiation (J-N), papillary cell clear cell carcinoma (O-S), cc-RCC with papillary growth (T-X), hemangioma (Y-Ba), simple cyst (Ca-Fa), RCC with micro cysts (Ga-Ka), oncocytoma (La-Pa). First row: gross appearance of the whole kidney or tumor after nephrectomy, second row: anatomical reference (after processing) which is used to manually draw a mask of the whole kidney and tumor (T2-TSE), third row: the unweighted image of the diffusion scan after processing and masking (DWI-b0), fourth row: a merge of the fraction maps from the tri-exponential fit, red = ffast, blue = finterm, green = fslow (1- finterm - ffast), fifth row: diffusion coefficient from the tri-exponential fit (Dtri). (TIF 36317 kb)

Additional file 2:

Figure S2. Dmono plotted against other diffusion coefficients MD (A), Dbi (B) and Dmono (C) for each lesions. MD versus Dmono displays good correlation, whereas Dbi and Dtri are structurally lower. Cor-heal = healthy cortex, med-heal = healthy medulla, ccRCC = clear cell renal cell carcinoma, pRCC = papillary cell Rhema = haemangioma, onco = oncocytoma, RCC, ccRCC-pRCC = ccRCC with papillary growth, ccRCC-sarc = ccRCC with sarcomatoid differentiation, ccRCC-cyst is ccRCC with micro-cystic structures. (TIF 1288 kb)

Additional file 3:

Table S1. Comparison of diffusion coefficients MD, Dmono, Dbi, Dtri. (PDF 682 kb)

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

van Baalen, S., Froeling, M., Asselman, M. et al. Mono, bi- and tri-exponential diffusion MRI modelling for renal solid masses and comparison with histopathological findings. Cancer Imaging 18, 44 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s40644-018-0178-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s40644-018-0178-0

Keywords