Am J Physiol Endocrinol Metab 306: E1460 –E1467, 2014. First published April 22, 2014; doi:10.1152/ajpendo.00591.2013.

Innovative Methodology

Automated quantification of pancreatic ␤-cell mass Maria L. Golson,1 William S. Bush,2 and Marcela Brissova1 1

Division of Diabetes, Endocrinology, and Metabolism, Department of Medicine, Vanderbilt University Medical Center, Nashville, Tennessee; and 2Center for Human Genetics Research, Vanderbilt University Medical Center, Nashville, Tennessee Submitted 25 October 2013; accepted in final form 19 April 2014

Golson ML, Bush WS, Brissova M. Automated quantification of pancreatic ␤-cell mass. Am J Physiol Endocrinol Metab 306: E1460–E1467, 2014. First published April 22, 2014; doi:10.1152/ajpendo.00591.2013.—␤Cell mass is a parameter commonly measured in studies of islet biology and diabetes. However, the rigorous quantification of pancreatic ␤-cell mass using conventional histological methods is a timeconsuming process. Rapidly evolving virtual slide technology with high-resolution slide scanners and newly developed image analysis tools has the potential to transform ␤-cell mass measurement. To test the effectiveness and accuracy of this new approach, we assessed pancreata from normal C57Bl/6J mice and from mouse models of ␤-cell ablation (streptozotocin-treated mice) and ␤-cell hyperplasia (leptin-deficient mice), using a standardized systematic sampling of pancreatic specimens. Our data indicate that automated analysis of virtual pancreatic slides is highly reliable and yields results consistent with those obtained by conventional morphometric analysis. This new methodology will allow investigators to dramatically reduce the time required for ␤-cell mass measurement by automating high-resolution image capture and analysis of entire pancreatic sections. pancreatic islets; ␤-cell mass; automated analysis; diabetes PANCREATIC ␤-CELL MASS closely correlates with ␤-cell function, and it has become a parameter commonly examined in studies of islet biology and diabetes (14, 18, 19, 24, 26, 27). To accurately measure ␤-cell mass using histological sections, ⬃1.5% of the total pancreas volume needs to be systematically sampled (5). This represents a time-consuming and laborious task, since each pancreatic specimen needs to be sectioned, stained, imaged, and analyzed. Reliable measurement of ␤-cell mass in histological specimens depends on two fundamental steps: 1) the capability to capture an image of an entire pancreatic section, and 2) the analysis of digitalized sections. Conventionally, investigators imaged pancreatic sections either by low-resolution slide scanning or with a microscope at low magnification (⫻4 –5) (2, 11, 20, 25). Alternatively, smaller areas of pancreas were imaged at higher magnification with a microscope-based setup and then analyzed individually or tiled before analysis (10, 13, 25). The first approach leads to errors in identifying small islets or individual ␤-cells due to low-resolution image capture. Although microscope-based image acquisition can be performed at higher magnifications, this process is relatively slow. In addition, tiling of individual images may cause problems in subsequent image analysis steps because the edges of each image are often darker than the center. While vignetting at the margins of each tiled image can be corrected by blank-field correction, the actual image processing to achieve this goal is not intuitive and often vignetting cannot be completely elimi-

Address for reprint requests and other correspondence: M. Brissova, Div. of Diabetes, Endocrinology, and Metabolism, Dept. of Medicine, Vanderbilt Univ. Medical Center, Nashville, TN 37232 (e-mail: [email protected]). E1460

nated due to variation in image intensity between blank and tissue images (17). Image analysis for quantification of ␤-cell mass was initially performed by point counting morphometry, pioneered by Susan Bonner-Weir (9, 15, 24, 27); more recently, however, investigators have been using several lowthroughput software platforms such as MetaMorph (Molecular Devices, Sunnyvale, CA), ImageJ (National Institutes of Health, Bethesda, MD), and Volocity (PerkinElmer, Santa Clara, CA) that measure ␤-cell area relative to the area of exocrine or total pancreatic tissue (2, 9, 11, 16, 20 –22). The advent of virtual slide technology together with the development of new analytical tools for virtual slides holds great promise to transform histological analysis of pancreas, which historically has been challenging because of a disproportionately small size of endocrine vs. exocrine compartment and nonuniform distribution of islets throughout the gland. Here, we took advantage of high-resolution slide scanning to assess its suitability for ␤-cell mass measurement. This approach used a web-based digital image information management interface, allowing remote virtual slide viewing, conferencing, archiving, and image analysis. To determine whether the slide scanning system can accurately measure ␤-cell mass, we trained an image analysis algorithm within its software suite and then used trained macros to examine pancreata from normal C57Bl/6J mice and from mouse models of ␤-cell ablation [streptozotocin (STZ)-treated mice] and ␤-cell hyperplasia (leptin-deficient Lepob/ob mice). The automated analysis of histological sections sampled systematically throughout the pancreas yielded results that were highly consistent with the outcomes of the conventional morphometric approach. Both automated and conventional analyses showed that ␤-cell mass in diabetic STZ-treated mice was reduced to 10% of the control, whereas nondiabetic Lepob/ob mice had ␤-cell mass increased nearly threefold, which was similar to previously reported data (1, 23, 28). Our data indicate that the automated analysis of virtual slides is appropriate for quantitative assessment of ␤-cell mass. With all the above taken together, we present a new methodology that will allow investigators to dramatically reduce the time required for ␤-cell mass measurement by automating highresolution image acquisition and analysis of entire pancreatic sections. This new approach can be adopted for the morphometric assessment of other pancreatic cell types and applied to the analysis of human pancreas specimens. MATERIALS AND METHODS

Mice and STZ administration. Wild-type and Lepob/ob mice (6, 7, 12) were on a C57Bl/6J background (Jackson Laboratories, Bar Harbor, ME). For STZ treatment, STZ (Sigma, St. Louis, MO) was freshly resuspended in cold 100 mM citrate buffer (pH 4.5) at 52.5 mg/ml and injected at 175 mg/kg body wt. Diabetic mice were euthanized 8 days after STZ injection (fed blood glucose 527 ⫾ 20 http://www.ajpendo.org

Innovative Methodology ANALYSIS OF ␤-CELL MASS

mg/dl). All analyses were performed on 10-wk-old mice. Some slides used for the training of the image analysis algorithm were obtained from B6D2 mice at postnatal day 8 (P8). Animal studies were approved by the Institutional Animal Care and Use Committee at Vanderbilt University Medical Center. Tissue fixation, immunohistochemistry, and imaging. After dissection, pancreata were weighed, laid flat in cassettes, fixed for 4 h in 4% paraformaldehyde, dehydrated, and embedded in paraffin (Fig. 1A). Longitudinal pancreatic sections were cut at a 5-␮m thickness, collected at 250-␮m intervals, and plated on glass slides. This resulted in the

E1461

collection of sections from six to eight depths per pancreas. Sections were immunolabeled with guinea pig anti-insulin (Dako, Carpenteria, CA; 1:1,000) and HRP-conjugated donkey anti-guinea pig antibodies (Jackson ImmunoResearch, West Grove, PA; 1:400). Insulin immunolabeling was visualized with a DAB kit (Vector Laboratories, Burlingame, CA), and sections were counterstained with Eosin-Y (Fisher Healthcare, Houston, TX). A schematic of pancreas processing is shown in Fig. 1A. All images were acquired on a ScanScope CS slide scanner (Aperio, Vista, CA) at ⫻20 magnification (resolution of 0.5 ␮m/pixel) and uploaded into a Spectrum digital slide interface (Aperio).

Fig. 1. Workflow for automated ␤-cell mass analysis. A: schematic of pancreas processing and sectioning. B: schematic of image analysis. Building a macro for accurate image segmentation involves repeated rounds of tissue selection, training, and testing. Steps that should be omitted when tuning a macro are depicted in gray; optional steps are shown in green; required steps are shown in beige.

AJP-Endocrinol Metab • doi:10.1152/ajpendo.00591.2013 • www.ajpendo.org

Innovative Methodology E1462

ANALYSIS OF ␤-CELL MASS

Building and testing of image analysis macros. A workflow schematic of image analysis is depicted in Fig. 1B. The Aperio software suite consists of a web-based interface (Spectrum v. 11.1.0, recently renamed to eSlide Manager), in which images are stored and batch data processing occurs, and a local workstation interface (ImageScope v. 11.2.0.780) where macro training takes place. Aperio image analysis has several algorithms, including the Genetic Imagery Exploitation toolset, with histology pattern recognition capability (Genie, Los Alamos National Laboratory, Los Alamos, NM), which we used to analyze brightfield images of pancreatic sections. To train the algorithm, images were chosen to reflect differences in sectioning and staining quality. Images of interest were copied within Spectrum and placed within a Genie project. We developed two image analysis macros. Macro 1 was set to perform analysis at ⫻20 magnification. This high magnification requires more computational time but can better select individual ␤-cells and irregularly shaped islets such as those seen in postnatal pancreas and STZ-treated mice. Macro 2 analyzed images at ⫻10 magnification, which accelerated analysis fourfold. Within the Genie project, three “Classes” of objects were defined, i.e., Glass, Eosin, and DAB (representing insulin labeling). Images within the Genie project were opened with ImageScope and regions of interest assigned to each Class. It is important to select a broad range of features for each Class to ensure that Genie correctly identifies objects of interest in the tissue section. A subset of these required features is depicted in Fig. 2. To fully encompass the variety in immonolabeling/counterstaining intensity and tissue morphological features necessary for Genie training, on average seven slides were used to generate each image analysis macro. For Glass, we selected clean areas as well as areas with some bubbles and dust. Areas including exocrine tissue, large blood vessels, connective tissue, and occasional pancreas-associated lymph nodes were classified as Eosin. For Eosin and DAB Classes, we selected areas with and without section folds (including areas in and out of focus) and areas that varied in immunolabeling/counterstaining intensity. When assigning areas with section folds to Classes, selecting folds with surrounding unfolded tissue resulted in more accurate object classification. In addition, regions assigned to the DAB Class had to be selected well within islet boundaries to avoid a haloing effect around islets. Despite high DAB background within connective tissue in some samples, the careful selection of connective tissue regions allowed the macros to distinguish between the connective tissue and islets. We also selected Eosin and DAB areas adjacent to each other to teach the algorithm the spatial relationship between Eosin and DAB Classes. After regions of interest were selected, Spectrum was used to create a Genie training set montage (referred to as “Montage” for short; Fig. 2). The Montage was then opened in ImageScope, and a new Training Macro was created and saved by using the “Analysis” menu and selecting the “Genie Training v1” algorithm. We adjusted the number of iterations in training to 500 before saving the Training Macro. Next, the Training Macro was run within Spectrum on the newly created Montage. Once the training job was completed, the Montage was reopened in ImageScope and viewed with image annotations and analysis color markup to evaluate the sensitivity and specificity of the Training Macro for the defined Classes of objects such as Glass, Eosin, and DAB (Fig. 3). Low sensitivity indicated that the Training Macro was not detecting the tissue of interest although it was assigned correctly in the Montage. The Training Macro had a low specificity for a given Class if it misclassified regions of interest in a slide. We aimed for a ⬎97% training success rate to achieve a highly accurate image segmentation. We continued with the macro training until no further improvement was observed in either sensitivity or specificity. If ⬎97% sensitivity and specificity were not reached, we examined what area of the slide was classified erroneously. More examples of the misclassified area were then added to the Montage, and macro training was repeated.

Fig. 2. Selected tissue features for training montages with assigned Class color markup. A: clean glass. B: lightly stained acinar tissue. C: sinuous ␤-cell cluster adjacent to eosin-stained connective tissue in D. E: dirty glass. F: connective tissue with background staining. G: connective tissue. H: folded, intensely-stained acinar tissue. I: light DAB labeling within islet. J: folded, intensely-labeled DAB within an islet. Color markup for analyzed images: Glass (green), Eosin (yellow), and DAB (red or blue). All images are displayed at ⫻10 magnification; C and D are from the Macro 1 Montage; A and B and E–J are from the Macro 2 Montage. Glass, greeen class markup; Eeosin, yellow color markup; DAB, blue in C or red color markup in I and J. Yellow boxes in A–J denote a slide/tissue elements added to a montage.

Once the Training Macro reached the optimal sensitivity and specificity in our Genie Training Set Montage, we created a Classifier Macro in ImageScope by using the “Analysis” menu and selecting “Genie Classifier v1”. In the image analysis window under “Classifier,” we selected the optimized Training Macro and saved it as the Classifier Macro. The performance of the Classifier Macro was then assessed on small regions of images that were not part of the Genie Training Set Montage. If the Classifier Macro aberrantly assigned tissue types, we added regions with poorly classified tissue to the Genie Training Set Montage and repeated the process above (see Fig. 1B). Analysis of digital slides with the Classifier Macro can be carried out on a local workstation, but we used server-based batch processing to increase analysis efficiency. The batch processing allowed us to select and analyze slides on a specimen basis using a Classifier Macro optimized for each specimen type (islet phenotype). After batch processing, data for each specimen (mouse) were directly exported into an Excel spreadsheet for further statistical analyses and calculation of ␤-cell mass. Any given macro will not be universal for every investigator. Differences in tissue processing and staining between laboratories can change the intensity and staining quality, and as a result of that, the macro needs to be retrained to assure high specificity and selectivity in tissue classification. Depending on the magnification used for

AJP-Endocrinol Metab • doi:10.1152/ajpendo.00591.2013 • www.ajpendo.org

Innovative Methodology ANALYSIS OF ␤-CELL MASS

analysis and the number of tissue features included in the training montage, one training session (includes 500 iterations) generally takes between 30 and 90 min with current software and server hardware. In our experience, a working macro can be generated in less than 24 h.

E1463

Since the training takes place on the server, the investigator is involved only in reviewing the macro after each training session and in determining whether the training requires adding more tissue features or is final. We predict that the training time will greatly decrease in the foreseeable future with ongoing improvements in server hardware (e.g., increasing number of processing cores and speed) and operating systems. Analysis of ␤-cell mass. Images of whole pancreatic sections acquired on a ScanScope CS slide scanner at ⫻20 magnification were analyzed by both conventional morphometry and a macro-based automated approach. Six to eight slides were analyzed per mouse by each method. Conventional morphometry was performed with MetaMorph v. 7.7. Because MataMorph software does not have the capacity to analyze images of an entire pancreatic section captured at ⫻20 magnification, all images were physically cropped into smaller regions prior to analysis. Islets were manually traced using the MetaMorph region tool, and then islet and exocrine tissue areas were analyzed as described previously (11). Analysis of ␤-cell mass using this conventional approach required ⬃10 h per mouse. Automated analysis of ␤-cell mass ran in the server environment without any investigator supervision. Image analysis at ⫻10 magnification (defined in Classifier Macro) took between 10 and 25 min and increased to 20 – 60 min per slide at ⫻20 magnification depending on section size. The output of the Genie Classifier Macros was then used to calculate ␤-cell mass. Since the glass portion of the slide is included in image analysis, the Genie Classifier Macros return DAB%, Eosin%, Glass%, and total analysis area for each slide. To calculate DAB% in a pancreatic tissue section the following equation was used: tissue DAB% ⫽ 关共slide DAB%兲/共slide DAB% (1) ⫹ slide Eosin%兲兴 ⫻ 100 Since each slide had a different amount of total tissue, the following equations were used to calculate the total Tissue, DAB, and Eosin area per slide: slide tissue area ⫽ scanned slide area ⫻ 关共slide DAB% (2) ⫹ slide Eosin%兲 ⫻ 100兴 slide DAB area ⫽ slide tissue area ⫻ (tissue DAB%/100)

(3)

slide Eosin area ⫽ slide tissue area ⫺ slide DAB area

(4)

To calculate DAB% per pancreas the following equation was used: cumulative pancreas DAB% n n ⫽ ( i⫽1 i ⫽ slide DAB arean) ⁄ [( i⫽1 slide DAB arean) (5) n ⫹ ( i⫽1 slide Eosin arean)] ⫻ 100







where n ⫽ number of slides/pancreas/mouse. ␤-Cell mass was calculated by the following equation: ␤-cell mass ⫽ pancreas weight ⫻ 共cumulative pancreas DAB%/100兲 (6) Statistical analysis. Mouse weights, pancreas weights, and ␤-cell mass were analyzed using Student’s t-test with two tails and one-way

Fig. 3. Pattern recognition algorithm allows for automated detection of islets in pancreatic sections. Pancreatic sections from wild-type (A and B=), STZ-treated (C and D=), and Lepob//ob mice (E and F=) on a C57Bl6/J background were labeled for insulin and counterstained with eosin. Slides were scanned at ⫻20 magnification. A–F: original images; insulin (brown), eosin (pink). A=–F=: color markup of tissue Classes identified in A–F by image analysis with Macro 1 (A=–D=) or Macro 2 (E= and F=). Glass (green), Eosin (yellow), DAB (blue in A=–D= or red in E= and F=). Scale bar in A, C, and E is 1 mm; scale bar in B, D, and F is 100 ␮m. AJP-Endocrinol Metab • doi:10.1152/ajpendo.00591.2013 • www.ajpendo.org

Innovative Methodology E1464

ANALYSIS OF ␤-CELL MASS

ANOVA using Prism v. 6 (GraphPad Software, La Jolla, CA). Two-way repeated-measures ANOVA was performed using the anova function of STATA 12.1 (StataCorp, College Station, TX) with slide as the repeated measure. A standard two-way ANOVA was used to compare ␤-cell mass measured by automated vs. conventional method across three mouse strains. Post hoc t-tests were performed to examine pairwise differences among mouse strains. P values ⬍ 0.05 were considered significant. RESULTS

New-generation slide scanning system allows rapid, highresolution image acquisition of full pancreatic sections. To determine whether an automated slide scanning and analysis system can accurately assess ␤-cell mass, we used three mouse models with different islet phenotypes; C57Bl6/J controls and STZ-treated and leptin-deficient (Lebob/ob) mice. Body weights (15.1 ⫾ 0.7 g, control; 11.8 ⫾ 1.3g, STZ-treated; and 31.9 ⫾ 1.8g, Lepob/ob mice) and wet pancreas weights (280 ⫾ 21 mg, control; 182 ⫾ 12 mg, STZ-treated; and 270 ⫾ 12 mg, Lepob/ob mice) of all three cohorts were consistent with prior reports (1, 6, 7, 12). Pancreatic sections were immunolabeled and imaged at ⫻20 magnification using the workflow outlined in Fig. 1A. The mean and median sizes of sections were 113.6 and 113.7 mm2, respectively, and the mean and median scan times per section were 73.49 and 75.2 s, respectively. Classifier Macros automate image analysis of pancreatic sections. High-resolution images of pancreatic sections were used to design Classifier Macros for automated analysis of ␤-cell mass (Fig. 1B). Detailed instructions for building Classifier Macros are outlined in MATERIALS AND METHODS. We reasoned that using pancreatic sections with islet features encompassing both normal islets and islets from STZtreated mice would allow us to develop a generic pattern recognition macro sufficiently versatile for all three islet phenotypes. With this goal in mind, we selected P8 tissues for Genie algorithm training, since islets at this stage are small and ␤-cells are often not adjacent to each other, similar to those in STZ-treated animals. To ensure high resolution for single ␤-cell detection, this Macro was trained at ⫻20 magnification using a Genie Training Set Montage that contained 32 different slide/tissue features. The final Classifier Macro (Macro 1) correctly recognized all three classes of objects, Glass, Eosin, and DAB (Fig. 3, A–D=), and was able to detect individual ␤-cells (Fig. 3B=). We noted that, in some pancreatic sections from the STZ-treated group, small areas between ␤-cells were occasionally misclassified as DAB, and glass in between eosinstained tissue was occasionally misclassified as Eosin. However, this small variation in Macro 1 specificity was ⬍10% and did not affect ␤-cell mass quantification in the STZ-treated cohort. Initially, we thought that Macro 1 would be sufficiently sensitive and specific for all three types of tissues, but in Lepob/ob pancreas it was unable to detect the interior of some islets. This may have been due to larger islet size or the lower intensity of DAB staining in Lepob/ob islets that could be associated with decreased insulin content per ␤-cell (8). Because of this lower islet recognition by Macro 1 in Lepob/ob pancreas and memory limits imposed by the size of our Training Montage, we created Macro 2 (Fig. 3, E–F=). To improve the detection of Lepob/ob islets, we incorporated more slide/tissue features into the Genie Training Set Montage (total

of 42 features; some examples shown in Fig. 2). Some of these features included lightly stained islets that are common to Lepob/ob pancreas (Fig. 2I). However, the larger Montage required decreasing the magnification for Macro training from ⫻20 to ⫻10 to balance image processing capacity and analysis speed. Although Macro 2 accurately classified normal and Lepob/ob islets (DAB), it was not suitable for detection of individual ␤-cells in STZ-treated animals due to the lower magnification selected for training. We also noted that the Training Macro frequently displayed a “halo” effect, where some pixels around islets in the adjacent exocrine tissue were classified as DAB (data not shown). This effect was minimized in the analysis by selecting DAB⫹ areas well within islet borders and adding this feature to the Genie Training Set Montage; see Fig. 2, I and J, for examples of interior outlining of islets that results in minimal haloing in Fig. 3, E–F=. Because of the vast differences in islet morphological features between STZ-treated and Lepob/ob mice, we subsequently used Macro 1 for analysis of ␤-cell ablation and Macro 2 to measure ␤-cell hyperplasia. This approach ensured high fidelity of islet detection and analysis accuracy. In contrast, Macros 1 and 2 could be used interchangeably for analysis of ␤-cell mass in C57Bl6/J controls with ⬍10% variation in the analysis results (data not shown). Automated image analysis accurately measures ␤-cell mass. Once the design and testing of Macros 1 and 2 were successfully completed, we applied them to the measurement of tissue DAB% (Eq. 1) in control, STZ-treated, and Lepob/ob mice. Each pancreas was subjected to a systematic section sampling at 250-␮m intervals throughout the tissue block (Fig. 1A). ␤-Cell mass was then calculated using Eq. 6. We first examined the distribution of islets throughout the pancreas by plotting tissue DAB% in relationship to section depth within the pancreas. Our data indicate that tissue DAB% varied considerably with section depth in all three groups of mice (Fig. 4, A–C; note different scale on y-axis in panels A–C). A similar observation was reported by Chintinne at al. (5), who analyzed ␤-cell mass in adult rat pancreas. The average coefficient of variation was 45% in controls (range 20 –77%), 46% in STZ-treated mice (range 20 – 85%), and 32% in Lepob/ob mice (range 20 – 41%). One outlier in the Lepob/ob group had tissue DAB% similar to that of controls (Fig. 4C). Despite this variation, tissue DAB% was significantly different among control, STZ-treated, and Lepob/ob mice at every tissue depth examined (Fig. 4D). The exception was the last portion of some tissue blocks (slide 8), where sections were available from only two to three animals per group (Fig. 4D). Next, we validated our automated macro-based approach by comparing the automated and conventional morphometric analyses for the entire set of scanned images from each mouse. To test whether our macros produced significantly different pancreas DAB% measurements across all mice, we performed a two-way repeated-measures ANOVA, considering each slide/tissue section to be a repeated measure within each mouse (Fig. 5A). We modeled a between effect to compare our automated analysis vs. the conventional method. If our macros produced significantly different pancreas DAB% results across this diverse set of samples, we would see a statistically significant effect of the analysis choice; this effect was not significant (P ⫽ 0.6044). On the basis of these results, we found no

AJP-Endocrinol Metab • doi:10.1152/ajpendo.00591.2013 • www.ajpendo.org

Innovative Methodology ANALYSIS OF ␤-CELL MASS

E1465

DISCUSSION

Defining how pancreatic ␤-cell mass relates to function requires accurate measurement of ␤-cell mass. This task has been challenging due to low abundance of islets and their heterogeneous distribution throughout the pancreas. Here, we show that these challenges can be overcome with a rapidly evolving virtual slide technology that provides automated highresolution image acquisition and analysis of pancreatic sections. Histological analysis of ␤-cell mass relies on the capability to image entire pancreatic sections. New-generation slide scanning systems, such as the one used in the current study, allow automation of this process by capturing images of entire tissue sections at high magnification (⫻20 – 40) with unprecedented speed and quality. For example, the imaging time required to generate our dataset from all three animal cohorts was only 2.5 h. Completing this task with conventional instrumentation would likely take several days if not weeks. In contrast to older scanners, new slide scanning at high magnification dramatically increased image resolution, allowing detection of individual ␤-cells. This is especially important for the monitoring

Fig. 4. Percentage of tissue DAB representing pancreatic section area labeled for insulin varies in the pancreas with section depth. Percentage of tissue DAB was plotted in relationship to the tissue section location within the pancreas. A: Control mice, n ⫽ 5. B: STZ-treated mice, n ⫽ 5. C: Lepob/ob mice, n ⫽ 5. Note different scale on y-axis in A, B, and C. D: percentage of tissue DAB plotted as a function of section depth was averaged for all mice per each group. *P ⬍ 0.05, **P ⬍ 0.01, ***P ⬍ 0.001.

significant difference between our automated analysis and the conventional approach (Fig. 5A). We also considered the ␤-cell mass computed across all mice and performed a two-way ANOVA. The effect of the analysis choice was not statistically significant (P ⫽ 0.1863), but, as expected, differences in ␤-cell mass across the three mouse strains were highly significant (P ⬍ 0.0001). We then performed pairwise post hoc t-tests between mouse strains for each method (Fig. 5B). Each pairwise comparison was statistically significant. Our data showed that, independently of the analysis method, diabetic STZ-treated mice had ␤-cell mass reduced by 90% compared with controls, whereas ␤-cell mass was increased approximately threefold in the leptin-deficient mice, which is consistent with previously published reports (1, 23, 28).

Fig. 5. Validation of automated analysis for ␤-cell mass measurement. A: pancreas DAB% was not significantly different when analyzed by automated or conventional method across all mice (P ⫽ 0.6044). B: automated measurements of ␤-cell mass in control, STZ-treated, and Lepob/ob mice were compared with results generated by conventional morphometry. As expected, ␤-cell mass was significantly reduced in STZ-treated mice and increased in Lepob/ob mice compared with WT mice (P ⬍ 0.0001), but no significant difference was detected between automated and conventional analysis method (P ⫽ 0.1863). Automated analysis: ****P ⬍ 0.0001, STZ-treated vs. WT mice; ###P ⬍ 0.001, STZ-treated vs. Lepob/ob mice; **P ⬍ 0.01, Lepob/ob vs. WT mice. Conventional analysis: ****P ⬍ 0.0001, STZ-treated vs. WT mice; ##P ⬍ 0.01, STZ-treated vs. Lepob/ob mice; *P ⬍ 0.05 Lepob/ob vs. WT mice.

AJP-Endocrinol Metab • doi:10.1152/ajpendo.00591.2013 • www.ajpendo.org

Innovative Methodology E1466

ANALYSIS OF ␤-CELL MASS

of ␤-cell mass loss, which is frequently associated with significantly perturbed islet architecture. In addition to new-generation slide scanners, virtual slide technology offers web-based digital slide management and, more importantly, new tools for automation of digital slide analysis. This automation is based on specialized algorithms for image segmentation and is further supported by batch processing of digital slides. Our data from three different strains of mice demonstrate that automated analysis of ␤-cell mass yields results that are highly consistent with those obtained by conventional morphometry and with previously published data (1, 23, 28). To automate ␤-cell mass measurement, we used the pattern recognition algorithm because it is able to classify image components based on a number of features including morphology, color, and texture. This is a great advantage compared with a simple color deconvolution approach, since the pattern recognition algorithm can be trained to develop macros that can easily distinguish the staining in ␤-cells from the background in non-␤-cell compartments such as red blood cells, connective tissue, or acinar tissue. We used DAB labeling for ␤-cell visualization and eosin as the tissue counterstain, but other color development systems and counterstains can be utilized. Changing the stain or counterstain color will require developing a new macro, since the color is one of the tissue features used by the algorithm for object classification. On the other hand, the variation in staining intensity can be addressed by a simple retraining of the existing classifier macro. While the pattern recognition algorithm can be trained to recognize various morphological features of interest, our data suggest that the current version of the algorithm has some limitations (mostly related to the computational power of the current hardware) when it comes to development of a single macro capable of detecting a broad range of islet phenotypes, such as those in Lepob/ob and STZtreated mice. This may also apply to ␤-cell detection, for example, in embryonic and postnatal pancreatic tissues or in human pancreas where ␤-cells are more intermingled with other endocrine cell types. Nevertheless, one can circumvent this limitation by developing an islet phenotype-specific macro with high specificity and selectivity for ␤-cell detection. Furthermore, we show that some macros can be shared between different experiments. For example, Macro 1 was applicable to ␤-cell mass analysis in normal, STZ-treated, and neonatal mice (data not shown). As we mentioned in MATERIALS AND METHODS, the tissue Classifier Macro reports data summary per each slide. However, the image analysis suite used in our study has provisions to combine tissue classifier macros with other algorithms, and this allows reporting parameters such as area, proliferation, etc. on an islet basis, which we will explore in the future. Training of Classifier Macro and image analysis run in a server environment can be performed with minimal time input from an investigator in contrast to conventional morphometry. The automation of image analysis not only saves the investigator an enormous amount of time but also provides consistency and eliminates observer bias. Our data showed that the tissue DAB% varied considerably with section depth in the pancreas, even within the same experimental group; thus, choosing only one or a few random slides within the pancreas will lead to inaccurate ␤-cell mass results (5). This variation in tissue DAB% is likely associated with inherent differences in

islet location within individual animals, tissue embedding, and the tissue sectioning process. Our data suggest that sections sampled systematically from four to five depths while avoiding the tissue periphery are sufficient to accurately assess differences between mice with robust variances in ␤-cell mass, such as those in control, STZ-treated, and Lepob/ob mice. However, more subtle differences in ␤-cell mass may require increasing the frequency for pancreatic section sampling, as for example those reported recently by Rankin et al. (23). We propose a standardized approach for quantifying mouse ␤-cell mass (outlined in Fig. 1) that includes 1) measuring pancreas wet weight, 2) embedding pancreas in a head-to-tail orientation, 3) systematic sampling of longitudinal sections throughout the pancreas, 4) automated imaging of whole pancreatic sections, and 5) automating image analysis. The systematic sampling of human pancreas and measurement of human ␤-cell mass are more complex due to pancreas size, but the protocols developed recently by Campbell-Thompson and colleagues (3, 4) should allow this procedure to be standardized. In conclusion, having a fast, reliable, and automated method for histological measurement of ␤-cell mass will greatly benefit analysis of islet/pancreas phenotypes and validation of new approaches and biomarkers for imaging of ␤-cell mass in vivo. Moreover, we envision that this methodology for automation of ␤-cell morphometry can be applied to the analysis of other islet endocrine cell types and human pancreatic specimens. ACKNOWLEDGMENTS We thank Greg Poffenberger, Nora Kayton, and Anastasia Coldren for technical assistance with the mice and Drs. Maureen Gannon and Alvin C. Powers (Vanderbilt University) for their critical reading of the manuscript and helpful suggestions. GRANTS This work was partially supported by grants to Dr. Alvin C. Powers from the National Institute of Diabetes and Digestive and Kidney Diseases (DK69603, DK-68764, DK-89572, DK-66636, DK-63439, DK-72473, DK-89538) and the Juvenile Diabetes Research Foundation International and a Merit Review Award from the Veterans Affairs Research Service (BX000666). Imaging and analysis of pancreatic sections were performed in the Vanderbilt Islet Procurement and Analysis Core supported by the Vanderbilt Diabetes Research and Training Center (DK-20593). DISCLOSURES No conflicts of interest, financial or otherwise, are declared by the author(s). AUTHOR CONTRIBUTIONS Author contributions: M.L.G. and M.B. conception and design of research; M.L.G. performed experiments; M.L.G., W.S.B., and M.B. analyzed data; M.L.G., W.S.B., and M.B. interpreted results of experiments; M.L.G. and W.S.B. prepared figures; M.L.G. and M.B. drafted manuscript; M.L.G., W.S.B., and M.B. edited and revised manuscript; M.L.G., W.S.B., and M.B. approved final version of manuscript. REFERENCES 1. Bock T, Pakkenberg B, Buschard K. Increased islet volume but unchanged islet number in ob/ob mice. Diabetes 52: 1716 –1722, 2003. 2. Brissova M, Shostak A, Shiota M, Wiebe PO, Poffenberger G, Kantz J, Chen Z, Carr C, Jerome WG, Chen J, Baldwin HS, Nicholson W, Bader DM, Jetton T, Gannon M, Powers AC. Pancreatic islet production of vascular endothelial growth factor-A is essential for islet vascularization, revascularization, and function. Diabetes 55: 2974 –2985, 2006. 3. Campbell-Thompson ML, Heiple T, Montgomery E, Zhang L, Schneider L. Staining protocols for human pancreatic islets. J Vis Exp: e4068, 2012.

AJP-Endocrinol Metab • doi:10.1152/ajpendo.00591.2013 • www.ajpendo.org

Innovative Methodology ANALYSIS OF ␤-CELL MASS 4. Campbell-Thompson ML, Montgomery EL, Foss RM, Kolheffer KM, Phipps G, Schneider L, Atkinson MA. Collection protocol for human pancreas. J Vis Exp: e4039, 2012. 5. Chintinne M, Stangé G, Denys B, In’t Veld P, Hellemans K, PipeleersMarichal M, Ling Z, Pipeleers D. Contribution of postnatally formed small beta cell aggregates to functional beta cell mass in adult rat pancreas. Diabetologia 53: 2380 –2388, 2010. 6. Coleman DL, Hummel KP. The influence of genetic background on the expression of the obese (Ob) gene in the mouse. Diabetologia 9: 287–293, 1973. 7. Coleman DL. Obese and diabetes: two mutant genes causing diabetesobesity syndromes in mice. Diabetologia 14: 141–148, 1978. 8. Dai C, Brissova M, Reinert RB, Nyman L, Liu EH, Thompson C, Shostak A, Shiota M, Takahashi T, Powers AC. Pancreatic islet vasculature adapts to insulin resistance through dilation and not angiogenesis. Diabetes 62: 4144 –4153, 2013. 9. Drott CJ, Olerud J, Emanuelsson H, Christoffersson G, Carlsson PO. Sustained beta-cell dysfunction but normalized islet mass in aged thrombospondin-1 deficient mice. PLoS ONE 7: e47451, 2012. 10. Golson ML, Le Lay J, Gao N, Brämswig N, Loomes KM, Oakey R, May CL, White P, Kaestner KH. Jagged1 is a competitive inhibitor of Notch signaling in the embryonic pancreas. Mech Dev 126: 687–699, 2009. 11. Golson ML, Misfeldt AA. High fat diet regulation of-cell proliferation and-cell mass. Open Endocrinol 4: 66 –77, 2010. 12. Ingalls AM, Dickie MM, Shell GD. Obese, a new mutation in the house mouse. J Hered 41: 317–318, 1950. 13. Jetton TL, Lausier J, Larock K, Trotman WE, Larmie B, Habibovic A, Peshavaria M, Leahy JL. Mechanisms of compensatory beta-cell growth in insulin-resistant rats: roles of Akt kinase. Diabetes 54: 2294 – 2304, 2005. 14. Kahn SE, Carr DB, Faulenbach MV, Utzschneider KM. An examination of beta-cell function measures and their potential use for estimating beta-cell mass. Diabetes Obes Metab 10, Suppl 4: 63–76, 2008. 15. Kulkarni RN, Holzenberger M, Shih DQ, Ozcan U, Stoffel M, Magnuson MA, Kahn CR. beta-cell-specific deletion of the Igf1 receptor leads to hyperinsulinemia and glucose intolerance but does not alter beta-cell mass. Nat Genet 31: 111–115, 2002. 16. Liu Z, Tanabe K, Bernal-Mizrachi E, Permutt MA. Mice with beta cell overexpression of glycogen synthase kinase-3beta have reduced beta cell mass and proliferation. Diabetologia 51: 623–631, 2008.

E1467

17. Marty GD. Blank-field correction for achieving a uniform white background in brightfield digital photomicrographs. BioTechniques 42: 716 – 718 –720, 2007. 18. Matveyenko AV, Butler PC. Relationship between beta-cell mass and diabetes onset. Diabetes Obes Metab 10, Suppl 4: 23–31, 2008. 19. McCulloch DK, Koerker DJ, Kahn SE, Bonner-Weir S, Palmer JP. Correlations of in vivo beta-cell function tests with beta-cell mass and pancreatic insulin content in streptozocin-administered baboons. Diabetes 40: 673–679, 1991. 20. Peshavaria M, Larmie BL, Lausier J, Satish B, Habibovic A, Roskens V, Larock K, Everill B, Leahy JL, Jetton TL. Regulation of pancreatic beta-cell regeneration in the normoglycemic 60% partial-pancreatectomy mouse. Diabetes 55: 3289 –3298, 2006. 21. Rachdi L, Balcazar N, Osorio-Duque F, Elghazi L, Weiss A, Gould A, Chang-Chen KJ, Gambello MJ, Bernal-Mizrachi E. Disruption of Tsc2 in pancreatic beta cells induces beta cell mass expansion and improved glucose tolerance in a TORC1-dependent manner. Proc Natl Acad Sci USA 105: 9250 –9255, 2008. 22. Rankin MM, Kushner JA. Adaptive beta-cell proliferation is severely restricted with advanced age. Diabetes 58: 1365–1372, 2009. 23. Rankin MM, Wilbur CJ, Rak K, Shields EJ, Granger A, Kushner JA. ␤-Cells are not generated in pancreatic duct ligation-induced injury in adult mice. Diabetes 62: 1634 –1645, 2013. 24. Steil GM, Trivedi N, Jonas JC, Hasenkamp WM, Sharma A, BonnerWeir S, Weir GC. Adaptation of ␤-cell mass to substrate oversupply: enhanced function with normal gene expression. Am J Physiol Endocrinol Metab 280: E788 –E796, 2001. 25. Teta M, Rankin MM, Long SY, Stein GM, Kushner JA. Growth and regeneration of adult beta cells does not involve specialized progenitors. Dev Cell 12: 817–826, 2007. 26. Weir GC, Bonner-Weir S. Islet ␤ cell mass in diabetes and how it relates to function, birth, and death. Ann NY Acad Sci 1281: 92–105, 2013. 27. Xu G, Stoffers DA, Habener JF, Bonner-Weir S. Exendin-4 stimulates both beta-cell replication and neogenesis, resulting in increased beta-cell mass and improved glucose tolerance in diabetic rats. Diabetes 48: 2270 –2276, 1999. 28. Yin D, Tao J, Lee DD, Shen J, Hara M, Lopez J, Kuznetsov A, Philipson LH, Chong AS. Recovery of islet beta-cell function in streptozotocin- induced diabetic mice: an indirect role for the spleen. Diabetes 55: 3256 –3263, 2006.

AJP-Endocrinol Metab • doi:10.1152/ajpendo.00591.2013 • www.ajpendo.org

Copyright of American Journal of Physiology: Endocrinology & Metabolism is the property of American Physiological Society and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.

Automated quantification of pancreatic β-cell mass.

β-Cell mass is a parameter commonly measured in studies of islet biology and diabetes. However, the rigorous quantification of pancreatic β-cell mass ...
3MB Sizes 2 Downloads 4 Views