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 (
link)]. 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 (
link)]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 (
link)] 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 (
link)]. The anatomical delineation was performed on a matched balanced turbo field echo (BTFE) image. BTFE images contain a mix of T
1 and T
2 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 mm
2 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 (
link)]. In view of potential uneven distribution of parenchymal fat between different regions of the pancreas observed in some [33 (
link)–37 (
link)] but not all studies [12 (
link), 13 (
link), 38 (
link)–40 (
link)], 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 (
link)].
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.
Al-Mrabeh A., Hollingsworth K.G., Steven S., Tiniakos D, & Taylor R. (2017). Quantification of intrapancreatic fat in type 2 diabetes by MRI. PLoS ONE, 12(4), e0174660.