Quantification of intrapancreatic fat in type 2 diabetes by MRI

Abstract

Objectives

Accumulation of intrapancreatic fat may be important in type 2 diabetes, but widely varying data have been reported. The standard quantification by MRI in vivo is time consuming and dependent upon a high level of experience. We aimed to develop a new method which would minimise inter-observer variation and to compare this against previously published datasets.

Methods

A technique of ‘biopsying’ the image to minimise inclusion of non-parenchymal tissues was developed. Additionally, thresholding was applied to exclude both pancreatic ducts and intrusions of visceral fat, with pixels of fat values of <1% or >20% being excluded. The new MR image ‘biopsy’ (MR-opsy) was compared to the standard method by 6 independent observers with wide experience of image analysis but no experience of pancreas imaging. The effect of the new method was examined on datasets from two studies of weight loss in type 2 diabetes.

Results

At low levels of intrapancreatic fat neither the result nor the inter-observer CV was changed by MR-opsy, thresholding or a combination of the methods. However, at higher levels the conventional method exhibited poor inter-observer agreement (coefficient of variation 26.9%) and the new combined method improved the CV to 4.3% (p<0.03). Using either MR-opsy alone or with thresholding, the new methods indicated a closer relationship between decrease in intrapancreatic fat and fall in blood glucose.

Conclusion

The inter-observer variation for quantifying intrapancreatic fat was substantially improved by the new method when pancreas fat levels were moderately high. The method will improve comparability of pancreas fat measurement between research groups.

About the Authors:

Ahmad Al-Mrabeh

Affiliation: Magnetic Resonance Centre, Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, United Kingdom

ORCID http://orcid.org/0000-0003-3689-1266

Kieren G. Hollingsworth

Affiliation: Magnetic Resonance Centre, Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, United Kingdom

Sarah Steven

Affiliation: Magnetic Resonance Centre, Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, United Kingdom

Dina Tiniakos

Affiliation: Department of Cellular Pathology, Royal Victoria Infirmary, Newcastle upon Tyne Hospitals and Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, United Kingdom

Roy Taylor

* E-mail: [email protected]

Affiliation: Magnetic Resonance Centre, Institute of Cellular Medicine, Newcastle University, Newcastle upon Tyne, United Kingdom

Introduction

Type 2 diabetes never develops without substantial impairment of intrapancreatic insulin secretory capacity [1, 2]. This appears to be closely linked to increased pancreatic parenchymal fat with reduction of this depot during acute weight loss being specific to type 2 diabetes [3–7]. However, these studies are based on a time consuming, expert-dependent method of analysing magnetic resonance data, and other groups have reported a much wider range of values for pancreatic fat content in type 2 diabetes [8–16]. Controversy about results has followed [16, 17]. There is a pressing need for an easily reproducible method for precise measurement which will allow comparability between research centres. Adipose tissue expansion and intra organ fat accumulation are correlated with high levels of lipid inflammatory markers that cause oxidative stress in obesity and type 2 diabetes [18], and calorie restriction was reported to be associated with decrease in oxidative stress caused by lipid peroxidation [19]. The recent demonstration that exposure of pancreatic islets to increased fatty acids causes beta cell de-differentiation, and that this is the likely underlying mechanism for type 2 diabetes, further emphasizes the importance of precise quantification of pancreatic fat content [20, 21].

Magnetic resonance techniques allow non-invasive in vivo quantification of pancreatic fat [7, 13, 14, 22], other techniques having lower sensitivity [13, 23–26]. Values of greater than 20% have been reported [8, 9, 11], well in excess of histological estimation [24, 27]. Quantification of the percentage of fat within the parenchyma of the pancreas is challenging as it depends upon delineating tissue entirely within the organ. Our recent description of the involuted nature of the pancreas in type 2 diabetes offers some insight into the reported variability [28, 29]. The border of the normal pancreas is irregular, but far more so in type 2 diabetes. As the irregularity has been reported to be directly proportional to the fat content of the pancreas [28], it is possible that any inter-lobular intrusion of visceral fat might be interpreted as intra-pancreatic fat and hence overestimate the true value. Additionally, the volume of the pancreas is decreased by 30% in type 2 diabetes of recent onset and by 50% in type 2 diabetes of duration greater than 10 years [28, 29], making inclusion of visceral fat in any defined volume much more likely in type 2 diabetes compared with a non-diabetic group. This would obscure small differences in true parenchymal fat content and could explain at least part of the variability of intra-pancreatic fat content reported.

The original method of minimising inclusion of visceral fat intruding into the pancreas involved freehand drawing round a region to be sampled. This demands considerable experience and is extremely time consuming [5–7]. It is also prone to variation between observers. A simple rapid method which would yield consistent values between different observers would be of great benefit to further investigation of the role of intrapancreatic fat in the pathogenesis of type 2 diabetes. We have developed a new simplified method of analysing intrapancreatic fat content by MRI which minimises the extent of inclusion of extrinsic tissues. Additionally, we have applied the new method to previously published data sets [5, 6] to investigate whether its use would change the pathophysiological conclusions of previous work.

Methods

Pancreatic fat quantification

We collectively refer to intra-lobular and interlobular fat as “intrapancreatic fat” or “pancreas fat” as MRI cannot distinguish between these two fat compartments. This terminology is used throughout this paper.

MRI data were acquired using a 3.0 Tesla Philips Achieva scanner (Philips, Best, The Netherlands) with a 6 channel cardiac array for signal detection. The protocol consisted of matched breath-held acquisitions of (i) a 3 point Dixon acquisition to quantify the intrapancreatic triglyceride and (ii) a balanced turbo field echo image to aid anatomical delineation of the pancreas[29]. Another 3-point Dixon acquisition was prescribed at the level of the L4-L5 intervertebral space to estimate subcutaneous and visceral fat areas in this slice. The 3 point Dixon method [30]acquires three gradient-echo scans during one breath-hold with adjacent out-of-phase and in-phase echoes (repetition time/echo times/averages/flip angle = 50ms/3.45, 4.60, 5.75ms/1/5°, bandwidth 435Hz/pixel). Field-of-view was set according to patient size (400-480x300mm), zero filled to give a resolution of 1.39x1.40mm. 12 sections of 5mm thickness were used to image the pancreas during two 17-second breath-holds, while one section was acquired at L4-L5 Custom MATLAB software was used to model the fat and water contributions to the gradient echo signals using a spectral model of fat with 6 peaks based on [31] and a single R2* component. Proton density fat fraction maps (the fat signal expressed as a percentage of the total signal) were constructed taking account of noise bias[32]. The anatomical delineation was performed on a matched balanced turbo field echo (BTFE) image. BTFE images contain a mix of T1 and T2 contrast, which distinguishes high signal intensity from vessels with visceral fat with lower intensity signals from the pancreas. It can therefore be used to clearly delineate the boundaries of the pancreas from adjacent structures, including the surrounding visceral fat, the splenic vein, the superior mesenteric vessels the inferior vena cava and duodenum. Twelve axial sections of 5mm thickness were imaged during an eight second breath-hold (repetition time/echo time/flip angle = 3.1ms/1.6ms/40°, turbo factor 95, parallel imaging factor 2, bandwidth 1156Hz per pixel). The field of view and zero filled resolution were matched to the 3 point Dixon imaging.The conventional method of freehand drawing round an area to be within the substance of the pancreas and a newly developed MR image ‘biopsy’ method (MR-opsy) were compared. For both methods, the regions of interest were selected to be within the parenchymal tissues and avoiding areas of visceral fat, main blood vessels.

For the conventional method, the ImageJ Polygon tool was used to select a region of interest in the parenchymal tissue of the pancreas head, body and tail. The region was selected to be as large as possible whilst being clear of the pancreas borders to avoid any possible contamination of surrounding visceral fat (Fig 1A).For MR-opsy, the Oval tool of ImageJ was used to select three regions of interest (~100 mm2 each) to represent equally the pancreas head, body and tail, the size of selection was chosen after pilot studies to permit easy placement entirely within the pancreas considering the irregularity in pancreas morphology (Fig 1A and 1B) [16]. In view of potential uneven distribution of parenchymal fat between different regions of the pancreas observed in some [33–37] but not all studies [12, 13, 38–40], sampling regions were placed equally throughout the pancreas to avoid possible bias. Analysis of both study datasets using the conventional methodology as originally published was carried out by experts experienced in pancreas anatomy. This was performed blinded to glucose tolerance and all clinical and metabolic markers both in the original studies and the present comparative study. Visceral and subcutaneous fat areas at L4-L5 were calculated from the L4-L5 proton density fat fraction map by thresholding and watershed analysis [41].

[Figure omitted. See PDF.]

Fig 1. Illustration of sampling methods for intrapancreatic quantification.

A: Representative MRI (3-point Dixon) slice of the pancreas was selected (upper panels). An anatomical scan was also acquired in parallel to the Dixon scan for localization purpose and for differentiation pancreas parenchymal tissue from main vessels and other adjacent abdomen tissues (lower panels). Regions of interest were carefully positioned away from pancreas borders to avoid contamination from visceral fat and away from main vessels. Conventional ROI: Polygon tool of ImageJ was used to select single ROI on the head, body and tail of pancreas away from visceral fat and main vessels. MR-opsy: Three ROIs (100mm2 each) were placed uniformly to represent different parts of the pancreas using ImageJ Oval tool away from visceral fat and main vessels. B: Magnified region of the pancreas to illustrate the size of biopsy selection (100mm2) relative to the size of an individual pixel (1.93mm2 = 1 pixel).The software reshape the oval selection (b, right) to take the nearest pixel shape (b, left).

https://doi.org/10.1371/journal.pone.0174660.g001

A step-by-step description of the process is presented in the Supplementary Methods section.

Two representative slices were selected to be assessed by each method and pancreatic fat content was calculated as the average pancreatic fat fraction of both slices.

Thresholding of fat measurement

Each image slice through the pancreas is 5mm thick to permit an adequate signal to noise ratio in the fat fraction images. In order to eliminate potential contribution of non- parenchymal tissue (visceral fat, pancreatic duct or blood vessel) within the selected region, a threshold was applied to both methods by collecting the histogram data within the area of selection and computing the resulted data to exclude pixels values outside the threshold limits which would otherwise contribute to the mean value (see step-by-step description in S1 Methods). Anonymised histological sections of pancreatic parenchymal tissues from people undergoing pancreatic surgery taken from various locations in the pancreas showed adipocyte distribution similar to the upper limit of 20% reported by Pinnick et al [27]. Hence, the maximum number of adipocytes clustered within a single voxel of pancreatic parenchymal tissues is estimated to be approximately 4000 (Fig 2), and any MRI fat signal above 20% is likely to be due to contamination by visceral fat tissue.

[Figure omitted. See PDF.]

Fig 2. An illustration of adipocytes distribution within single MRI voxel of parenchymal tissues of pancreas.

Histological section of background normal pancreatic tissue of a 48-year female undergoing pancreatectomy for a neuroendocrine tumour. The average size of single adipocyte is approximately 100μm, but adipocytes can occur in clusters. Based on average adipocyte size of 100μm, the maximum number of adipocytes likely to be present in one voxel is ~400x50 = 20000 adipocytes. The upper threshold of 20% assumes that the maximum number to be 4000 adipocytes within a single voxel of pancreas.

https://doi.org/10.1371/journal.pone.0174660.g002

Similarly, pixels almost devoid of fat (<1%) are likely to represent major pancreatic ducts or blood vessels, and these cannot be discriminated on the BTFE image. The main pancreatic duct network is suggested though not segmentable on T2-weighted images, such as Fig 3A1 (a T2 weighted fast spin echo, TR/TE = 946ms/70ms with spectrally selective adiabatic inversion recovery, SPAIR, for fat saturation). A thresholding range of 1–20% was therefore applied to the original data, and the performance of both methods was compared before and after thresholding. The rationale is explained in Fig 4.

[Figure omitted. See PDF.]

Fig 3. Example of ductal system architecture within the pancreas.

A: different MRI axial acquisitions of the pancreas in T2DM subject (a1: T2-SPAIR, a2: BTFE, a3: 3-point Dixon). Pancreas of T2-SPAIR (a1) sequence was segmented and volume rendered in Drishti as described before [28], volume rendered image was colour tagged then opacity level was manipulated to show the distribution of pancreatic ductal system in white colour (Drishti version 2.6.3).

https://doi.org/10.1371/journal.pone.0174660.g003

[Figure omitted. See PDF.]

Fig 4. Colour map of pancreatic fat distribution in type 2 diabetes.

The colour map shows the wide range of fat distribution within the sampling area. This underlies the rationale for thresholding to exclude non-parenchymal tissues. The colour bar on the right shows fat level from 0% (dark blue) to 25% (red). Threshold levels were set to exclude areas of fat content less than 1% (possible blood vessels or main duct) or above 20% (visceral fat contamination). Parenchymal fat was considered to range between 1–20%. a1-a6 represent areas of varied fat content within the single MR-opsy selection (a1: < 1%, a2: 1–5%, a3:6–10%, a4:11–15%, a5:16–20%, a6: > 20%).

https://doi.org/10.1371/journal.pone.0174660.g004

Reproducibility of fat quantification

Type 2 diabetes participants with low (3.3%) and high (6.5%) pancreas fat levels were analysed to test the reproducibility of methods. MR scans acquired prior to weight loss were examined by 6 independent observers using the methods in random order. The observers have wide experience of image analysis but no previous experience of assessing pancreas. Each was asked to follow instructions and quantify intrapancreatic fat by both methods. Coefficient of variation (CV) was calculated for the 6 independent measurements and compared by both methods for the two participants.

Intervention studies design

The Counterbalance study tested the durability of type 2 diabetes reversal after a very low calorie diet in a group of 30 people with 0.5–23 years of diabetes duration [5]. Intrapancreatic fat was quantified at baseline, after 8 weeks of very low calorie diet (VLCD) and following a 6 months weight maintenance programme. Participants were considered as responders if fasting plasma glucose level <7mmol/l following VLCD and return to normal diet [5].

The bariatric surgery study evaluated change in intrapancreatic fat after weight loss at 8 weeks post-surgery in groups with (n = 18) or without (n = 9) type 2 diabetes [6].

Statistical analysis

Minitab 17 (Minitab Inc, State College, PA, USA) was used for statistical analysis. All data are presented as mean ±SEM, and p value <0.05 was considered statistically significant. Student’s paired t-test was used to measure significance.

Results

Inter-observer agreement

At low levels of intrapancreatic fat, neither the result nor the inter-observer CV (coefficient of variation) was changed by MR-opsy, thresholding or a combination of the methods (Fig 5A). At high levels of pancreatic fat, the conventional method used by non-expert observers exhibited poor inter-observer precision (CV 26.9%; Fig 5B). Application of the MR-opsy method improved the CV to 4.3% (p<0.03; Fig 5B). The components of the improvement were separately assessed. MR-opsy alone improved the precision (CV 3.5%; p = 0.02) as did application of thresholding but to a lesser extent (CV 15.2%; p<0.05). Areas of selections by the observers were investigated in order to understand the difference in CV between the methods despite similarity in mean intrapancreatic fat percentage. It was found that some observers were more or less conservative in their perception of the boundary of the pancreas and the visceral fat or areas of blood vessels that could lead to both over-estimation or under-estimation of pancreatic fat content, respectively (Fig 6).

[Figure omitted. See PDF.]

Fig 5. Reproducibility of fat quantification methods.

The inter-observer variation for each method is shown for low level (3%) of pancreatic fat (A), and high level (6%) of pancreatic fat (B). Data for both methods are shown with and without 1–20% thresholding. * p<0.05 Conventional vs MR-opsy before thresholding. ‡ p<0.05 Conventional vs MR-opsy after thresholding. † p<0.05 Conventional without thresholding vs with thresholding.

https://doi.org/10.1371/journal.pone.0174660.g005

[Figure omitted. See PDF.]

Fig 6. Example of areas selected by two observers using conventional and MR-opsy methods.

ROIs of participants with the high level of pancreatic fat (6%) were shown using conventional (a,b), and MR-opsy (c,d) methods for observers 1 and 2. Two regions with potential contribution to wrong estimation of pancreatic fat content were selected: Region 1 represents a region of focal fat which appears bright on the Dixon scan (a,c), and dark on the anatomical BTFE scan (b,d). Region 2 represents blood vessel and appears dark on the Dixon scan (a,c), and bright on the anatomical scan (b,d). It is clear that observer 1 included both areas within the selection using conventional method whereas biopsy method avoided those regions by the same observer.

https://doi.org/10.1371/journal.pone.0174660.g006

Effect of fat quantification method on counterbalance study data

At baseline, using the conventional method in expert hands, there was no significant difference in intrapancreatic fat between those who subsequently were or were not able to reverse their T2DM by weight loss (5.3±0.4% vs. 5.9±0.7%; Table 1). Application of MR-opsy, with or without thresholding did not change this (Table 1).

[Figure omitted. See PDF.]

Table 1. Counterbalance study: Pancreas fat change in responders and non-responders before and after weight loss.

https://doi.org/10.1371/journal.pone.0174660.t001

After the 8 week weight loss period using the conventional method there was a significant fall in intrapancreatic fat in both responder and non-responder groups. Use of MR-opsy, with or without thresholding did not change the significance of the decrease in intrapancreatic fat in the responders. In the non-responders the significant fall reported by the conventional method was not observed by any of the new methods (Table 1).

After the 6 month weight maintenance period, using the conventional method there was a significant fall in intrapancreatic fat in both responder and non-responder groups. In the responders, all methods observed the previously reported significant decrease in intrapancreatic fat content. In the non-responders, using MR-opsy with or without thresholding no significant change was observed (Table 1). Use of thresholding alone in the non-responders decreased the apparent extent of fall in pancreas fat (Table 1).

Effect of fat quantification method on bariatric surgery study

At baseline, use of any of the methods showed intrapancreatic fat to be significantly higher in the group with type 2 diabetes compared with the normal glucose tolerant group (conventional: 6.6±0.5% vs. 5.1±0.2%; MR-opsy with thresholding: 6.4±0.3 vs. 5.1±0.6%; Table 2). The fall in intrapancreatic fat during weight loss in the type 2 diabetic group remained significant using all methods. Conversely, weight loss brought about no change within the NGT group between baseline and 8 weeks after surgery using conventional or new methods (Table 2).

[Figure omitted. See PDF.]

Table 2. Bariatric surgery study: Pancreas fat change in type 2 diabetes (T2DM) and normal glucose tolerance (NGT) participants before and after weight loss.

https://doi.org/10.1371/journal.pone.0174660.t002

Variability of fat distribution within the pancreas areas

Fat distribution varied significantly between the head and other parts of the pancreas for the Counterbalance study (Table 3). At the baseline of the study, fat percentage was higher in the head of the pancreas compared with the body or the tail using conventional, MR-opsy alone or MR-opsy with thresholding methods indicating heterogeneity among different pancreatic tissues in fat distribution (MR-opsy plus thresholding: head vs. body, p = 0.006; head vs. tail, p = 0.01).

[Figure omitted. See PDF.]

Table 3. Fat% in different regions of the pancreas before and after intervention studies.

https://doi.org/10.1371/journal.pone.0174660.t003

Correlation between pancreas fat and some body characteristics

No correlation between body weight and BMI with pancreatic fat was found in the Counterbalance study using any quantification method. In the bariatric surgery study, there was correlation between pancreatic fat and weight using the conventional method only(r = 0. 5, p = 0.04). There was no correlation observed between fat content and age in both studies. Interestingly, we found significant correlation between pancreatic fat and diabetes duration within the Counterbalance study (r = 0.48, p = 0.008).

Discussion

Reproducible quantification of intra-parenchymal pancreas fat is important to allow comparisons between data from different research groups, and this is especially important as absolute differences in pancreas fat between type 2 diabetes and normal are modest [5–7]. We demonstrate that higher inter-observer agreement can be achieved using MR-opsy compared with the conventional region of interest method when intrapancreatic fat levels are higher and pancreas volume is lower [28]. As intrapancreatic fat increases and pancreas volume decreases with increasing disease duration [28, 29], the data are of particular relevance to this disease state. Re-analysis using the new method of previously published intervention studies of type 2 diabetes, which used conventional methodology applied by experts, did not change the previously reported pathophysiological implications.

Several studies have demonstrated the association between increased intrapancreatic fat and type 2 diabetes. In diabetes-prone rodent models of type 2 diabetes overfeeding brings about impairment of beta cell function, and this susceptibility to lipid availability is reflected in studies on isolated islets [3, 4, 42–44]. In humans predisposed to develop type 2 diabetes, prolonged Intralipid infusion severely impairs beta-cell function [45]. Conversely, removal of excess lipid from the environment of the pancreatic islet allows return of normal insulin secretion in early type 2 diabetes [5, 7]. This has also been observed in isolated islets [3]. The apparent relationship of this lipid depot to the pathophysiology of type 2 diabetes emphasizes the importance of methodology for precise measurement.

Homogeneity of fat distribution within the pancreas is a topic of great debate [12, 13, 33–40]. The series of studies on people with type 2 diabetes showed a degree of variability in fat content between head, body and tail. Given that the biological relevance of this work is to investigate any effect of fat upon overall beta cell function and that these are distributed throughout the pancreas, inclusion of data from each region in a mean to represent the whole pancreas is justified in order to represent fat distribution in the whole pancreas. Although selection of one region could be sufficient under certain conditions of homogeneous pancreas fat distribution such as in the study of very obese people, use of the MR-opsy method is still appropriate. A potential disadvantage of the method could arise if there was marked heterogeneity of fat content between regions of the pancreas, but the present observations and those of others suggest that this is rare. The pancreas in type 2 diabetes is 30–50% smaller than normal [28, 29]. The decrease in volume as diabetes duration increases is accompanied by notable increase in irregularity in the pancreas borders. This implies greater likelihood of inclusion of the extra-pancreatic fat which exists between lobules [28]. The contribution of pancreatic ducts or blood vessels which cannot be identified in the image, has previously been overlooked. It is notable that the mean level of pancreas fat increased as a result of 1–20% thresholding (Tables 1 and 2). The conventional method of pancreas fat quantification using magnetic resonance imaging has resulted in a wide range of reported pancreas fat content [7, 13, 46] and the present data suggest that this would be minimised by use of MR-opsy. In the present study the observers, who were experienced in image analysis but not in studying the pancreas, reported that placement of the 100 mm2 MR-opsies was not challenging, and was also rapid (approximately 5 minutes vs. up to 30 minutes for conventional drawing round a region of interest).

Short duration type 2 diabetes can be reversed after weight loss with restoration of normal beta cell function and this has been reported to be associated with a fall in intrapancreatic fat content [5–7]. Application of the new method resulted in identification of no change in intrapancreatic fat in the longer duration group (which did not respond to weight loss by normalising plasma glucose). These subjects had smaller, more irregular pancreases than the responders, and the new method is more likely to reflect true intra-pancreatic fat levels. In the bariatric study, the type 2 diabetes participants exhibited a good return to normal glucose control [6], and quantitation of intrapancreatic fat by either convention method in expert hands or by the new method showed a significant decrease.

Several studies reported the robustness of MR-based fat quantification methods [22, 47, 48]. A recent phantom study evaluated the reproducibility of MRI fat quantification technique between research centres, MR scanner vendors, field strengths, and acquisition protocols [49] emphasizing the importance of standardized image analysis technique for precise comparison. However, use of a phantom does not reflect the complexities introduced by variable inclusion of visceral fat and fluid filled intra-organ ducts.

The published studies employing MR to quantify pancreas fat content used a wide range of methods for sampling size and location selected for fat quantification (Table 4). This can partially explain the discrepancy in reported pancreas fat content. Of these studies, some reported a significant relationship between diabetes or insulin resistance and pancreas fat [6, 8–12, 14, 15, 29, 40, 50, 51]. Other studies reported no significant difference in pancreas fat content between type 2 diabetes and non-diabetic controls [7, 10, 13, 37, 52, 53].

[Figure omitted. See PDF.]

Table 4. Summary of up-to-date studies employed MR for fat quantification in the pancreas.

https://doi.org/10.1371/journal.pone.0174660.t004

The remaining studies did not compare between non-diabetic and diabetic groups [6, 22, 35, 38, 39, 46, 54–60]. Whereas magnetic resonance imaging methods allow subsequent selection of the volume to analyse, magnetic resonance spectroscopy depends upon acquiring data from a volume of the body pre-selected by imaging. Consequently it is particularly susceptible to inclusion of visceral fat due to respiratory and other movement in the scanner. Such spectroscopy methods tend to report higher pancreas fat content (up to 24%). Hu et al. [22] reported that MR spectroscopy was less accurate than imaging for pancreas fat quantification due to the difficulty in voxel positioning, and this is consistent with the present data on the effect of selection of region of interest for fat quantification. By combining (a) selection of several regions to represent tissues from the whole organ; (b) restriction of size of the selected region to decrease contamination from visceral fat; and (c) applying thresholding to exclude contribution from non-parenchymal tissues, an improvement in inter-observer agreement is observed.

The major limitation of the current study is the lack of a gold standard for non-invasive quantitation of fat solely within the parenchymal tissue of the pancreas. At present, neither Dixon nor anatomical scans can differentiate between parenchymal tissue and ductal or small vascular structures in the pancreas. Optimization of image acquisition for differentiating between those small structures is demanding and currently being developed. For example the T2- SPAIR sequence used in image 3 might allow exclusion of the main pancreatic duct within the MR-opsy selection. As one of the major limiting factors for resolution of pancreas imaging is breath-hold duration, development of sparse scanning techniques which acquire data more rapidly may be expected to permit higher resolution imaging [61]. Under condition of severe pancreas fat infiltration of parenchymal tissues, the performance of the MR-opsy method alone can be limited. Nonetheless, the proposed 20% threshold to exclude areas of visceral fat invasion remains useful under such circumstances, and values close to 20% should trigger detailed examination of the pancreas anatomy when selecting regions of interest. In conclusion, quantification of fat within the pancreas by MRI is significantly affected by the method of sampling and the new MR-opsy method allows higher inter-observer agreement. Application of this standardised new method with thresholding should permit measurement of changes in true intrapancreatic fat content which can reliably be compared between different research groups.

Supporting information

[Figure omitted. See PDF.]

S1 Methods.

https://doi.org/10.1371/journal.pone.0174660.s001

(DOC)

Acknowledgments

We thank all our participants, and L Ward, T Hodgson and D Wallace, research radiographers.

Author Contributions

1. Conceptualization: RT.

2. Data curation: KGH.

3. Formal analysis: AA-M SS KGH.

4. Funding acquisition: RT KGH.

5. Investigation: AA-M SS KGH.

6. Methodology: AA-M KGH RT.

7. Project administration: RT.

8. Resources: KGH DT.

9. Software: KGH.

10. Supervision: RT.

11. Validation: AA-M KGH RT.

12. Visualization: AA-M KGH RT.

13. Writing – original draft: AA-M.

14. Writing – review & editing: AA-M KGH SS DT RT.

Citation: Al-Mrabeh A, Hollingsworth KG, Steven S, Tiniakos D, Taylor R (2017) Quantification of intrapancreatic fat in type 2 diabetes by MRI. PLoS ONE 12(4): e0174660. https://doi.org/10.1371/journal.pone.0174660

References

1. Ferrannini E, Nannipieri M, Williams K, Gonzales C, Haffner SM, Stern MP. Mode of onset of type 2 diabetes from normal or impaired glucose tolerance. Diabetes. 2004;53(1):160–5. pmid:14693710

2. Tabak AG, Jokela M, Akbaraly TN, Brunner EJ, Kivimaki M, Witte DR. Trajectories of glycaemia, insulin sensitivity, and insulin secretion before diagnosis of type 2 diabetes: an analysis from the Whitehall II study. Lancet. 2009;373(9682):2215–21. Epub 2009/06/12. PubMed Central PMCID: PMC2726723. pmid:19515410

3. Pinnick K, Neville M, Clark A, Fielding B. Reversibility of metabolic and morphological changes associated with chronic exposure of pancreatic islet beta-cells to fatty acids. Journal of cellular biochemistry. 2010;109(4):683–92. pmid:20069570

4. Lee Y, Hirose H, Ohneda M, Johnson JH, McGarry JD, Unger RH. Beta-cell lipotoxicity in the pathogenesis of non-insulin-dependent diabetes mellitus of obese rats: impairment in adipocyte-beta-cell relationships. Proceedings of the National Academy of Sciences of the United States of America. 1994;91(23):10878–82. Epub 1994/11/08. PubMed Central PMCID: PMCPmc45129. pmid:7971976

5. Steven S, Hollingsworth KG, Al-Mrabeh A, Avery L, Aribisala B, Caslake M, et al. Very low calorie diet and 6 months of weight stability in type 2 diabetes: Pathophysiologic changes in responders and non-responders. Diabetes care. 2016a;

6. Steven S, Hollingsworth KG, Small PK, Woodcock SA, Pucci A, Aribisala B, et al. Weight Loss Decreases Excess Pancreatic Triacylglycerol Specifically in Type 2 Diabetes. Diabetes care. 2016b;39(1):158–65. Epub 2015/12/03.

7. Lim EL, Hollingsworth KG, Aribisala BS, Chen MJ, Mathers JC, Taylor R. Reversal of type 2 diabetes: normalisation of beta cell function in association with decreased pancreas and liver triacylglycerol. Diabetologia. 2011;54(10):2506–14. Epub 2011/06/10. PubMed Central PMCID: PMCPmc3168743. pmid:21656330

8. van der Zijl NJ, Goossens GH, Moors CC, van Raalte DH, Muskiet MH, Pouwels PJ, et al. Ectopic fat storage in the pancreas, liver, and abdominal fat depots: impact on beta-cell function in individuals with impaired glucose metabolism. The Journal of clinical endocrinology and metabolism. 2011;96(2):459–67. Epub 2010/11/19. pmid:21084401

9. Gaborit B, Abdesselam I, Kober F, Jacquier A, Ronsin O, Emungania O, et al. Ectopic fat storage in the pancreas using 1H-MRS: importance of diabetic status and modulation with bariatric surgery-induced weight loss. International journal of obesity (2005). 2015;39(3):480–7. Epub 2014/07/22. pmid:25042860

10. Begovatz P, Koliaki C, Weber K, Strassburger K, Nowotny B, Nowotny P, et al. Pancreatic adipose tissue infiltration, parenchymal steatosis and beta cell function in humans. Diabetologia. 2015;58(7):1646–55. Epub 2015/03/06. pmid:25740696

11. Tushuizen ME, Bunck MC, Pouwels PJ, Bontemps S, van Waesberghe JH, Schindhelm RK, et al. Pancreatic fat content and beta-cell function in men with and without type 2 diabetes. Diabetes care. 2007;30(11):2916–21. Epub 2007/08/02. pmid:17666465

12. Idilman IS, Tuzun A, Savas B, Elhan AH, Celik A, Idilman R, et al. Quantification of liver, pancreas, kidney, and vertebral body MRI-PDFF in non-alcoholic fatty liver disease. Abdominal imaging. 2015;40(6):1512–9. Epub 2015/02/27. pmid:25715922

13. Patel NS, Peterson MR, Brenner DA, Heba E, Sirlin C, Loomba R. Association between novel MRI-estimated pancreatic fat and liver histology-determined steatosis and fibrosis in non-alcoholic fatty liver disease. Alimentary pharmacology & therapeutics. 2013a;37(6):630–9. Epub 2013/02/07. PubMed Central PMCID: PMCPmc4136524.

14. Lingvay I, Esser V, Legendre JL, Price AL, Wertz KM, Adams-Huet B, et al. Noninvasive quantification of pancreatic fat in humans. The Journal of clinical endocrinology and metabolism. 2009;94(10):4070–6. Epub 2009/09/24. PubMed Central PMCID: PMCPmc2758729. pmid:19773401

15. Ma J, Song Z, Yan F. Detection of hepatic and pancreatic fat infiltration in type II diabetes mellitus patients with IDEAL-Quant using 3.0T MR: comparison with single-voxel proton spectroscopy. Chinese medical journal. 2014;127(20):3548–52. Epub 2014/10/16. pmid:25316227

16. Hollingsworth KG, Al-Mrabeh A, Steven S, Taylor R. Pancreatic triacylglycerol distribution in type 2 diabetes. Diabetologia. 2015;58(11):2676–8. Epub 2015/08/02. pmid:26232098

17. Begovatz P, Bierwagen A, Lundbom J, Roden M. Pancreatic triacylglycerol distribution in type 2 diabetes. Reply to Hollingsworth K. G., Al Mrabeh A., Steven S. et al [letter]. Diabetologia. 2015;58(11):2679–81. Epub 2015/09/25. pmid:26399402

18. Furukawa S, Fujita T, Shimabukuro M, Iwaki M, Yamada Y, Nakajima Y, et al. Increased oxidative stress in obesity and its impact on metabolic syndrome. The Journal of clinical investigation. 2004;114(12):1752–61. Epub 2004/12/16. PubMed Central PMCID: PMCPMC535065. pmid:15599400

19. Buchowski MS, Hongu N, Acra S, Wang L, Warolin J, Roberts LJ 2nd. Effect of modest caloric restriction on oxidative stress in women, a randomized trial. PloS one. 2012;7(10):e47079. Epub 2012/10/17. PubMed Central PMCID: PMCPMC3465282. pmid:23071718

20. Talchai C, Xuan S, Lin HV, Sussel L, Accili D. Pancreatic beta cell dedifferentiation as a mechanism of diabetic beta cell failure. Cell. 2012;150(6):1223–34. Epub 2012/09/18. PubMed Central PMCID: PMCPmc3445031. pmid:22980982

21. White MG, Shaw JAM, Taylor R. Type 2 Diabetes: The Pathologic Basis of Reversible β-Cell Dysfunction. Diabetes care. 2016;39(11):2080–8. pmid:27926891

22. Hu HH, Kim HW, Nayak KS, Goran MI. Comparison of fat-water MRI and single-voxel MRS in the assessment of hepatic and pancreatic fat fractions in humans. Obesity (Silver Spring, Md). 2010;18(4):841–7. Epub 2009/10/17. PubMed Central PMCID: PMCPmc2847037.

23. Saisho Y, Butler AE, Meier JJ, Monchamp T, Allen-Auerbach M, Rizza RA, et al. Pancreas volumes in humans from birth to age one hundred taking into account sex, obesity, and presence of type-2 diabetes. Clinical anatomy (New York, NY). 2007;20(8):933–42. Epub 2007/09/20. PubMed Central PMCID: PMCPmc2680737.

24. Kim SY, Kim H, Cho JY, Lim S, Cha K, Lee KH, et al. Quantitative assessment of pancreatic fat by using unenhanced CT: pathologic correlation and clinical implications. Radiology. 2014;271(1):104–12. Epub 2014/01/31. pmid:24475851

25. Ou HY, Wang CY, Yang YC, Chen MF, Chang CJ. The association between nonalcoholic fatty pancreas disease and diabetes. PloS one. 2013;8(5):e62561. Epub 2013/05/15. PubMed Central PMCID: PMCPmc3643962. pmid:23671610

26. Della Corte C, Mosca A, Majo F, Lucidi V, Panera N, Giglioni E, et al. Nonalcoholic fatty pancreas disease and Nonalcoholic fatty liver disease: more than ectopic fat. Clinical endocrinology. 2015;83(5):656–62. Epub 2015/07/24. pmid:26201937

27. Pinnick KE, Collins SC, Londos C, Gauguier D, Clark A, Fielding BA. Pancreatic ectopic fat is characterized by adipocyte infiltration and altered lipid composition. Obesity (Silver Spring, Md). 2008;16(3):522–30. Epub 2008/02/02.

28. Al-Mrabeh A, Hollingsworth K, Steven S, Taylor R. Morphology of the pancreas in type 2 diabetes: effect of weight loss with or without normalisation of insulin secretory capacity. Diabetologia. 2016;

29. Macauley M, Percival K, Thelwall PE, Hollingsworth KG, Taylor R. Altered Volume, Morphology and Composition of the Pancreas in Type 2 Diabetes. PloS one. 2015;10(5):e0126825. pmid:25950180

30. Glover GH, Schneider E. Three-point Dixon technique for true water/fat decomposition with B0 inhomogeneity correction. Magnetic resonance in medicine. 1991;18(2):371–83. Epub 1991/04/01. pmid:2046518

31. Hamilton G, Yokoo T, Bydder M, Cruite I, Schroeder ME, Sirlin CB, et al. In vivo characterization of the liver fat (1)H MR spectrum. NMR in biomedicine. 2011;24(7):784–90. Epub 2011/08/13. PubMed Central PMCID: PMCPMC3860876. pmid:21834002

32. Liu CY, McKenzie CA, Yu H, Brittain JH, Reeder SB. Fat quantification with IDEAL gradient echo imaging: correction of bias from T(1) and noise. Magnetic resonance in medicine. 2007;58(2):354–64. Epub 2007/07/27. pmid:17654578

33. Matsumoto S, Mori H, Miyake H, Takaki H, Maeda T, Yamada Y, et al. Uneven fatty replacement of the pancreas: evaluation with CT. Radiology. 1995;194(2):453–8. Epub 1995/02/01. pmid:7824726

34. Isserow JA, Siegelman ES, Mammone J. Focal fatty infiltration of the pancreas: MR characterization with chemical shift imaging. AJR American journal of roentgenology. 1999;173(5):1263–5. Epub 1999/11/30. pmid:10541101

35. Kim HJ, Byun JH, Park SH, Shin YM, Kim PN, Ha HK, et al. Focal fatty replacement of the pancreas: usefulness of chemical shift MRI. AJR American journal of roentgenology. 2007;188(2):429–32. Epub 2007/01/24. pmid:17242252

36. Pezzilli R, Calculli L. Pancreatic steatosis: Is it related to either obesity or diabetes mellitus? World journal of diabetes. 2014;5(4):415–9. Epub 2014/08/16. PubMed Central PMCID: PMCPmc4127578. pmid:25126389

37. Kühn JP, Berthold F, Mayerle J, Volzke H, Reeder SB, Rathmann W, et al. Pancreatic Steatosis Demonstrated at MR Imaging in the General Population: Clinical Relevance. Radiology. 2015;276(1):129–36. Epub 2015/02/07. PubMed Central PMCID: PMCPmc4554208. pmid:25658037

38. Li J, Xie Y, Yuan F, Song B, Tang C. Noninvasive quantification of pancreatic fat in healthy male population using chemical shift magnetic resonance imaging: effect of aging on pancreatic fat content. Pancreas. 2011;40(2):295–9. Epub 2010/12/24. pmid:21178651

39. Livingstone RS, Begovatz P, Kahl S, Nowotny B, Strassburger K, Giani G, et al. Initial clinical application of modified Dixon with flexible echo times: hepatic and pancreatic fat assessments in comparison with (1)H MRS. Magma (New York, NY). 2014;27(5):397–405. Epub 2013/12/07.

40. Chai J, Liu P, Jin E, Su T, Zhang J, Shi K, et al. MRI chemical shift imaging of the fat content of the pancreas and liver of patients with type 2 diabetes mellitus. Experimental and Therapeutic Medicine. 2016;11(2):476–80. pmid:26893633

41. Hallsworth K, Fattakhova G, Hollingsworth KG, Thoma C, Moore S, Taylor R, et al. Resistance exercise reduces liver fat and its mediators in non-alcoholic fatty liver disease independent of weight loss. Gut. 2011;60(9):1278–83. Epub 2011/06/29. PubMed Central PMCID: PMCPMC3152868. pmid:21708823

42. Elks ML. Chronic perifusion of rat islets with palmitate suppresses glucose-stimulated insulin release. Endocrinology. 1993;133(1):208–14. pmid:8319569

43. Lalloyer F, Vandewalle B, Percevault F, Torpier G, Kerr-Conte J, Oosterveer M, et al. Peroxisome proliferator-activated receptor alpha improves pancreatic adaptation to insulin resistance in obese mice and reduces lipotoxicity in human islets. Diabetes. 2006;55(6):1605–13. Epub 2006/05/30. pmid:16731822

44. Zhou YP, Grill VE. Long-term exposure of rat pancreatic islets to fatty acids inhibits glucose-induced insulin secretion and biosynthesis through a glucose fatty acid cycle. Journal of Clinical Investigation. 1994;93(2):870–6. pmid:8113418

45. Storgaard H, Jensen CB, Vaag AA, Volund A, Madsbad S. Insulin secretion after short- and long-term low-grade free fatty acid infusion in men with increased risk of developing type 2 diabetes. Metabolism: clinical and experimental. 2003;52(7):885–94.

46. Schwenzer NF, Machann J, Martirosian P, Stefan N, Schraml C, Fritsche A, et al. Quantification of pancreatic lipomatosis and liver steatosis by MRI: comparison of in/opposed-phase and spectral-spatial excitation techniques. Investigative radiology. 2008;43(5):330–7. Epub 2008/04/22. pmid:18424954

47. Reeder SB, McKenzie CA, Pineda AR, Yu H, Shimakawa A, Brau AC, et al. Water-fat separation with IDEAL gradient-echo imaging. Journal of magnetic resonance imaging: JMRI. 2007;25(3):644–52. Epub 2007/02/28. pmid:17326087

48. Sengupta S, Smith DS, Gifford A, Welch EB. Whole-body continuously moving table fat-water MRI with dynamic B0 shimming at 3 Tesla. Magnetic resonance in medicine. 2016;76(1):183–90. Epub 2015/07/23. PubMed Central PMCID: PMCPMC4775426. pmid:26198380

49. Hernando D, Sharma SD, Aliyari Ghasabeh M, Alvis BD, Arora SS, Hamilton G, et al. Multisite, multivendor validation of the accuracy and reproducibility of proton-density fat-fraction quantification at 1.5T and 3T using a fat-water phantom. Magnetic resonance in medicine. 2016. Epub 2016/04/16. PubMed Central PMCID: PMCPMC4835219.

50. Pacifico L, Chiesa C. Pancreatic fat and hepato-metabolic features in obese children with nonalcoholic fatty liver disease. Clinical endocrinology. 2015. Epub 2015/08/08.

51. Patel NS, Peterson MR, Lin GY, Feldstein A, Schnabl B, Bettencourt R, et al. Insulin Resistance Increases MRI-Estimated Pancreatic Fat in Nonalcoholic Fatty Liver Disease and Normal Controls. Gastroenterology research and practice. 2013b;2013:498296. Epub 2013/12/19. PubMed Central PMCID: PMCPmc3855930.

52. Heni M, Machann J, Staiger H, Schwenzer NF, Peter A, Schick F, et al. Pancreatic fat is negatively associated with insulin secretion in individuals with impaired fasting glucose and/or impaired glucose tolerance: a nuclear magnetic resonance study. Diabetes/metabolism research and reviews. 2010;26(3):200–5. Epub 2010/03/13. pmid:20225188

53. Wicklow BA, Griffith AT, Dumontet JN, Venugopal N, Ryner LN, McGavock JM. Pancreatic Lipid Content Is Not Associated with Beta Cell Dysfunction in Youth-Onset Type 2 Diabetes. Canadian journal of diabetes. 2015;39(5):398–404. Epub 2015/06/24. pmid:26099932

54. Sijens PE, Edens MA, Bakker SJ, Stolk RP. MRI-determined fat content of human liver, pancreas and kidney. World journal of gastroenterology: WJG. 2010;16(16):1993–8. Epub 2010/04/27. PubMed Central PMCID: PMCPmc2860076. pmid:20419836

55. Le KA, Ventura EE, Fisher JQ, Davis JN, Weigensberg MJ, Punyanitya M, et al. Ethnic differences in pancreatic fat accumulation and its relationship with other fat depots and inflammatory markers. Diabetes care. 2011;34(2):485–90. Epub 2011/01/29. PubMed Central PMCID: PMCPmc3024373. pmid:21270204

56. Szczepaniak LS, Victor RG, Mathur R, Nelson MD, Szczepaniak EW, Tyer N, et al. Pancreatic steatosis and its relationship to beta-cell dysfunction in humans: racial and ethnic variations. Diabetes care. 2012;35(11):2377–83. Epub 2012/09/13. PubMed Central PMCID: PMCPmc3476895. pmid:22968187

57. Targher G, Rossi AP, Zamboni GA, Fantin F, Antonioli A, Corzato F, et al. Pancreatic fat accumulation and its relationship with liver fat content and other fat depots in obese individuals. Journal of endocrinological investigation. 2012;35(8):748–53. Epub 2011/10/08. pmid:21979274

58. Kovanlikaya A, Mittelman SD, Ward A, Geffner ME, Dorey F, Gilsanz V. Obesity and fat quantification in lean tissues using three-point Dixon MR imaging. Pediatric radiology. 2005;35(6):601–7. Epub 2005/03/24. pmid:15785930

59. Wong VW, Wong GL, Yeung DK, Abrigo JM, Kong AP, Chan RS, et al. Fatty pancreas, insulin resistance, and beta-cell function: a population study using fat-water magnetic resonance imaging. The American journal of gastroenterology. 2014;109(4):589–97. Epub 2014/02/05. pmid:24492753

60. Cohen M, Syme C, Deforest M, Wells G, Detzler G, Cheng HL, et al. Ectopic fat in youth: the contribution of hepatic and pancreatic fat to metabolic disturbances. Obesity (Silver Spring, Md). 2014;22(5):1280–6. Epub 2014/01/10.

61. Mann LW, Higgins DM, Peters CN, Cassidy S, Hodson KK, Coombs A, et al. Accelerating MR Imaging Liver Steatosis Measurement Using Combined Compressed Sensing and Parallel Imaging: A Quantitative Evaluation. Radiology. 2016;278(1):247–56. Epub 2015/07/29. pmid:26218662