A Smartphone-Fluidic Digital Imaging Analysis System for Pancreatic Islet Mass Quantification

Introduction

Type 1 diabetes mellitus (T1DM) is characterized by the autoimmune destruction of insulin-producing beta-cells within pancreatic islets. Transplantation of pancreatic islets isolated from donated, cadaveric organs can restore normal blood glucose homeostasis. Next-generation approaches with stem-cell derived islets are currently being tested in clinical trials. Islets for therapeutic application are defined as a “biological drug” by the Food and Drug Administration (FDA) and, as such, must meet product release criteria such as purity, viability, mass, and functionality (Hering et al., 2016; Ricordi et al., 2016).

Islet mass is an important parameter that determines transplant outcomes (Shapiro et al., 2000; Gangemi et al., 2008; Wang et al., 2013; Qi et al., 2014). Islet mass quantification is often performed manually under a microscope after dithizone (DTZ) staining of zinc ions in beta-cells, which helps to differentiate islet and acinar tissues. Islet mass is then calculated in terms of islet equivalents (IEq). One IEq is defined as an islet with a diameter of 150 µm. Islets are classified algorithmically into groups by size using 50 µm diameter increments. A final IEq is calculated for each group using relative conversion factors (Buchwald et al., 2009). The manual method has been widely accepted for practice; however, it is highly subjective with significant inter-and intra-operator variability (Navarro-Fontestad et al., 2012).

Several computer-assisted digital imaging analysis (DIA) systems have been implemented which show reduced variability as compared to manual counting (Niclauss et al., 2008; Gmyr et al., 2015; Buchwald et al., 2016; Wang and Kaufman, 2016). However, several shortcomings of DIA still exist. For example, islets often aggregate with or overlay each other during sample preparation, resulting in difficulties in isolating individual islets. Additionally, static 2D-based DIA methods provide limited information on purity, morphology, and cell volume. Importantly, these methods still require operator intervention, leading to human errors and biases that may not meet FDA requirements of Current Good Manufacturing Practice (cGMP).

As an alternative to existing methods of islet quantification, we may look to unique, innovative solutions involving advanced technologies. For instance, smartphones equipped with high resolution complementary metal–oxide–semiconductor (CMOS) cameras, high-performance processing units, and tailored software can be used as analytical devices for biological research and clinical diagnostics (Hutchison et al., 2015). Assisted by advancements in molecular analysis, biosensors, mathematical algorithms, microfabrication, 3D-printing, and microfluidics, smartphones have been used as portable, versatile, and highly-connected read-out platforms capable of capturing the microscopic world ranging from tissues and cells to individual DNA molecules (Lee et al., 2011; Gmyr et al., 2015; McCracken and Yoon, 2016). The imaging capability of smartphones can be extended to the function of a microscope by adding an external lens, allowing for the capture of digital images and photographs with high resolution. Furthermore, with advances in wireless communication and with the help of cloud computing, complex analytical processing can be performed after data acquisition to generate diagnostic results (Im et al., 2015). Compared to conventional laboratory microscopes, smartphone technology provides a more rapid, portable, user-friendly, and cost-effective way to allow even minimally trained users to operate the system in the field.

Another helpful technology is microfluidics, which involves the precise manipulation of fluids at a microscopic scale. Microfluidics offers unique advantages for studying pancreatic islets by closely mimicking the physiological microenvironment. It also allows for the control of stimulation cues and integration of various analytical tools. Since the early 2000s, the development of in vitro microfluidic-based tools for diabetes research has drawn significant attention. Such microfluidic devices possess many advantages over more established in vitro assays, including reduced reagent consumption requirements, low cost of manufacturing and maintenance, multiplexing capabilities, increased assay sensitivity, increased accuracy, and higher spatiotemporal resolution (Kappler et al., 2004; Adewola et al., 2010; Nourmohammadzadeh et al., 2016; Xing et al., 2016; Lu et al., 2018). An exciting development in recent years is the adaptation of smartphone technology to microfluidic platforms for point-of-care (POC) chemical and biological detection, as well as for particle counting and analysis (Long et al., 2017; Huang X. et al., 2018; Liu et al., 2019; Hassan and Zhang, 2020).

In this paper, we describe a smartphone-based microfluidic system for human islet mass quantification that can provide more comprehensive characterization of isolated human islets. This system may be adopted in the future as a product release method for diabetes cell therapies.

Materials and Methods

Design and Fabrication of Smartphone-Fluidics-Based Flow Cytometry

The Smartphone-Fluidic Digital Imaging Analysis (SFDIA) system consists primarily of two components: a microfluidic device and a smartphone system. The smartphone system includes a 3D printed frame, a magnifying lens (f = 4.51 mm), a light source (Adafruit Industries, New York, NY), and a smartphone (Google Pixel 3, Foxconn, New Taipei, Taiwan) (Figure 1). The microchannel (500 μm in diameter and 500 μm in height) is made of one layer of PDMS (Polydimethylsiloxane, Fisher Scientific, Ontario, CA) using soft-photolithography (Adewola et al., 2010; Nourmohammadzadeh et al., 2016; Xing et al., 2016). This microchannel is bonded to a glass slide (Fisher Scientific, Ontario, CA) (Figure 1A), which is divided into an islet loading area with a repeated loop channel and a straight channel as a viewing area. The repeated loop channel has a total length of 700 mm and a total liquid volume of 250 μl. It is designed to preload islets, to act as a buffer zone for focusing islets in the middle of the channel, and to separate islets in distance. The straight channel is also 700 mm in length and contains a viewing area of 2 mm² for smartphone imaging and video recording. In the Field of View (FOV), there are two embedded scale markers of 1 and 0.1 mm (Figure 1A). These markers serve as a conversion scale calibrated to equivalent pixel density.

FIGURE 1

FIGURE 1. (A) Design of the microfluidic device. (i) The 1,000 µm marker used to calculate scaling factor $\text{δ}$. (ii) An example of video captured by the smartphone camera. (B) Hardware design for the smartphone-fluidics-based flow cytometry system. The 3D frame is designed to hold the smartphone at the top, the microfluidic device in the middle, and the LED light source at the bottom.

The smartphone frame is 3D-printed from polylactide resin (MakerBot® PLA resin, MakerBot® Industries, New York, NY) using a MakerBot® 3D printer (MakerBot® Industries, New York, NY). The frame is designed to hold the smartphone, the microfluidic device, and the optics system (magnifying lens and LED light source) (Figure 1B). The frame is 80 mm in width and 199.3 mm in length. There is a groove channel (20.50 mm in width and 3.21 mm in height) for holding fluids. Between the microfluidic channel and smartphone camera, there is a magnifying lens (A230, Thorlabs, Newton, NJ) held in a socket (6.00 mm in diameter). The lens has a diameter of 6.34 mm with a focal length (f) of 4.51 mm and a numerical aperture (NA) of 0.55. As illustrated in Figure 1B, an LED powered by a 3.0 V lithium battery serves as the illumination source and is placed at the bottom. The LED provides a consistent diffused light to optimize video acquisition. Videos of flowing islets are recorded at a rate of 30 fps by a CMOS camera on the Google Pixel 3 (12.2 MP, f/2.0). The camera is placed directly on top of the magnifying lens with an effective FOV of 2 mm².

Video Processing Algorithm

The first step in video processing is individual frame analysis. This step mainly consists of two parts: object detection and cell segmentation, as outlined in Figure 2. External environmental factors (lighting conditions, the specific microfluidic device used, etc.) are consistent when videos are taken. As such, background subtraction is a good way to separate background from foreground.

FIGURE 2

FIGURE 2. Flow chart of video processing algorithm. (A) Flow chart of the individual frame analysis. This can be further divided into object detection and cell segmentation. (B) Flow chart of the cell matching algorithm.

Figure 3 shows the typical processing sequence of a sample video frame. First, a background image is taken. The original islet image is captured directly by the smartphone camera (Figure 3A). Both the actual input video and the background image are then smoothed by a low pass filter. The foreground is obtained by subtracting the background from frames in the input video (Figure 3B). Finally, the foreground is converted into a black and white image using a modified version of Otsu thresholding.

FIGURE 3

FIGURE 3. Demonstration of cell segmentation. (A) Original image directly obtained from the smartphone camera. (B) Image created by subtracting background from the original image. (C) Image after applying distance transform. Since distance transform calculates the distance of each pixel from its nearest boundary, the central pixels have higher values. As a result, when applying the proper threshold to the transformed image, small-sized cells as well as peripheral areas of large cells (mainly connections between cells) can be removed, leaving only the central area of large cells. (D) Image created by watershed algorithm. Scale bar = 200 µm.

The black and white image is further processed with a distance transform, converting binary foreground pixel values into distances between themselves and the nearest background pixel. As a result, the pixels around the center of each blob contain higher values, while the pixels at the edge of the blob contain lower values. Another threshold is then applied to the transformed image followed by morphological opening. These remove unwanted connections between cells small particles which are considered to be fragments of tissues instead of islets. The resulting image shows only the major, central area of the cells, with unwanted particles removed (Figure 3C). By changing the threshold, the size of the small particles to be filtered out can be controlled. Finally, Watershed algorithm was used to group the non-central foreground pixels based on the cell center clusters made in the previous step. The result is shown in Figure 3D, with cells marked out with different colors.

After analysis of each individual frame, islet matching is performed to ensure that every islet is counted only once. As is outlined in Figure 2B, islet tracking is done using a Euclidean-distance-based approach (Chowdhury et al., 2010). A feature vector is formed for each detected islet. The vector includes the x-coordinate, y-coordinate, and size of the islet. Between two consecutive frames, the weighted Euclidean distance (known as $L_2$ distance) $J$ between each cell pair is calculated using the following equation:

$J=1.0\times \frac{{\left|\right|x_1-x_2\left|\right|}_{L_2}}{min\left(r_1, r_2\right)}+0.5\times \frac{{\left|\right|y_1-y_2\left|\right|}_{L_2}}{min\left(r_1, r_2\right)}+1.0\times \frac{{\left|\right|a_1-a_2\left|\right|}_{L_2}}{min\left(a_1, a_2\right)}(1)$

where $x_1$ and $x_2$ are the x-coordinates, $y_1$ and $y_2$ are the y-coordinates, $r_1$ and $r_2$ are the radii of the two cells, and $a_1$ and $a_2$ are the areas. The weights (1, 0.5, 1) are determined based on the fluid flow rate (100 μl/min) and video frame rate (30 fps). Individual islet clusters captured from two consecutive frames are matched based on shortest distance. Via further testing, we have shown that these weights can distinguish between matched and unmatched cells.

Quantification of Islet Mass

Size, Volume, and Islet Equivalents Quantification

Equivalent spherical diameter (ESD) has been widely used for size quantification of irregularly shaped particles (Jennings and Parslow, 1988). In this study, the equivalent spherical diameter of an individual islet is defined by the following equations:

$D_{pixel}=2\sqrt{\frac{A_{pixel}}{\text{π}}}(2)$

$D_{\mu m}=D_{pixel}\times \delta (3)$

where $D_{pixel}$ represents the islet’s equivalent diameter in pixels and $A_{pixel}$ represents the islet’s area in pixels. $D_{pixel}$ is converted to diameter in micrometers ($D_{\mu m}$) using a conversion factor $\delta$, which is calculated by measuring a built-in marker (1,000 µm) as is shown in Figure 1A. The calculation of $\delta$ is shown in Eq. 4. The length of the $1000 \mu m$ marker was measured to be 460 pixels.

$\delta =\frac{1000\mu m}{Lengt{h}_{pixel of 1000\mu m markers}}(4)$

IEq assessment is done, as previously established, by classifying $D_{pixel}$ into eight size ranges in μm (50–100, 101–150, 151–200, 201–250, 251–300, 301–350, 351–400, and >400). Then, IEq can be calculated by multiplying the number of islets in each group with corresponding multipliers (0.167, 0.667, 1.685, 3.499, 6.315, 10.352, 15.833, and 22.750, respectively) (Ricordi et al., 1990).

Under most circumstances, islets have irregular shapes. As a result, in 2-D images, an ellipse generally better represents the shape of an islet than a perfect circle (Girman et al., 2003). Thus, volume estimation is performed based on a 3-D Ellipsoid-Fitting-based algorithm. As demonstrated in Figure 4A, each islet is first fitted to an ellipse using least-squared approximation (Fitzgibbon et al., 1999) using a function provided by OpenCV™. Lengths of the major and minor axes were obtained. The 2-D ellipse can be converted to a 3-D ellipse by revolving the ellipse around the minor axis and using the major axis of the 2-D ellipse as the third axis of the ellipsoid as shown in Figure 4Ai (Girman et al., 2003). The volumes of islets in each frame are estimated using Eq. 5:

$V_{pixel}=\frac{4}{3}\pi M{a}^{2}Mi(5)$

where $Ma$ and $Mi$ represent the long axis and short axis, respectively. The final volume is the average volume across each frame as is defined in Eq. 6:

$\overline{V}=\frac{{{\displaystyle \sum \mathit{V}}}_{pixel}}{n}(6)$

$V_{pixel}$ represents the estimated volume of a cell in a single frame using pixels, and ${\overline{V}}_{pixel}$represents the average estimated volume of a cell using pixels. ${\overline{V}}_{pixel}$ can be then converted to ${\overline{V}}_{\mu m^3}$ using the scaling factor $\delta$ as shown in Eq. 7.

$\overline{V_{\mu m^3}}=\overline{V_{pixel}}\times {\delta }^3(7)$

The Ellipsoid-Fitting-Based Volume (EFV) was compared to the volume obtained based on IEq assessment (IEqV). According to the definition of IEq, IEqV can be calculated using the following equation:

$IEqV=IEq\times V_{150\mu m}(8)$

where $IEq$ is the Islet Equivalent obtained from the previous step and $V_{{150}_{\mu m}}$ is the volume of islets with a diameter of 150 μm, which is calculated to be 1.767 × 10⁶ μm³ using the sphere volume equation:

$V=\frac{4}{3}\pi r^3(9)$

FIGURE 4

FIGURE 4. (A) Demonstration of cell volume estimation. The two figures on the left show the process of ellipse fitting. The two figures on the right show the two different methods of estimating volume after ellipse fitting: (i) the calculation of ellipsoid volume when the third axis is the major axis; (ii) the same calculation when the third axis is the minor axis. (B) Demonstration of circularity and solidity assessment. The green lines in the two images represent the convex hull for the two islets. Solidity is defined as the ratio between actual cell size (white area) and the area surrounded by the convex hull. Scale bar = 200 µm.

When estimating ellipsoid volume based on 2-D ellipse information, this work initially used Eq. 5 as discussed earlier. Yet there is another way of such estimation, which is to revolve a 2-D ellipse around its major axis, and use the minor axis as the third axis of the ellipsoid as shown in Figure 4Aii (Niclauss et al., 2008) resulting in Eq. 10:

$V_{pixel}=\frac{4}{3}\pi MaM{i}^2(10)$

In this paper, both methods were implemented in order to see whether the choice of a different ellipsoid calculation method would have a significant impact on overall islet volume assessment.

Circularity and Solidity Quantification

A spherical model is generally used in the 3D representation of islets. Roundness (circularity) is used to measure the level of an islet’s shape regularity (Olehnik et al., 2017; Huang H.-H. et al., 2018). Circularity is a number ranging from 0 to 1.0. Higher circularity means that the shape of the given islet is closer to a circle. In this work, average circularity was calculated using the following equations:

$\overline{C}=\frac{\overline{A_{actual}}}{\overline{A_{estimated}}}(11)$

$\overline{{\mathit{A}}_{\mathit{estimated}}}=\frac{{\displaystyle {\sum }_0^{\mathit{n}}{\mathit{P}}^2}{}_{\mathit{actual}}/\left(4\pi \right)}{\mathit{n}}(12)$

where $\overline{C}$stands for average circularity, $\overline{A_{actual}}$ stands for actual islet area, $\overline{A_{estimated}}$ stands for average estimated particle area, and $P_{actual}$ stands for actual particle perimeter. $P_{actual}$ is obtained by counting the number of pixels on the islet’s contour, which is then converted to micrometers.

In addition to circularity, we applied the concept of solidity as an indication of islet shape regularity. The solidity is defined as the ratio of the area of the particle’s convex hull to the area of the actual particle:

$Solidity=\frac{\overline{A_{convexHull}}}{\overline{A_{actual}}}(13)$

The convex hull was obtained using a function provided by OpenCV™. Like circularity, solidity is a number ranging from 0 to 1.0. It is demonstrated in Figure 4B, where the green lines represent the convex hulls of the two islets. Since both islets have an elliptical shape, they have very similar circularities; however, the islet on the left is a bit more fragmented than the one on the right, and this difference is reflected in solidity.

Trapped Islet Percentage Estimation

Trapped islet percentage is estimated using the positive area of DTZ staining of zinc in the islets and computed as:

$\%DT{Z}^{+}ratio=\frac{Total dithizone positive area}{Total islet area}(14)$

The positive area normally produces a red color after DTZ staining, which can be identified based on Hue in the HSV (Hue, Saturation, Value) color space in digital image processing. In this work, pixels with a Hue value between 310 and 360 (Hue value is set in a range of 0–360) were identified as zinc positive areas.

Figure 5 includes example images of trapped islets (Figure 5A) and free islets (Figure 5B). The images on the left were taken directly from the smartphone camera. For better visualization, two heatmaps of HSV color space were generated. In the heatmaps, the DTZ stained areas appeared dark red/black, while the no-DTZ areas were bright red. 50% of the area of the two trapped islets in Figure 5A was stained with DTZ, while the free islets in Figure 5B were almost completely stained with DTZ. The DTZ positive ratios of the two trapped islets were calculated to be 67.31 and 66.73%, while the ratio of the free islets was calculated to be 99.05%.

FIGURE 5

FIGURE 5. Demonstration of trapped islet percentage. (A) Example of two trapped islets. (B) Example of a free islet. The two images on the left are raw images taken from the smartphone camera. The two images on the right are heatmaps of HSV color space. Scale bar = 200 µm.

Human Islet Preparation and Manual Mass Quantification

Human islets were isolated according to the published protocol (Qi et al., 2009a; b). In brief, a cadaver pancreas was obtained from organ procurement organizations (OPO) and islets were isolated at the University of Virginia with donors’ consent for research. The pancreas was trimmed and distended with Liberase HI (Roche Applied Science, IN) and digested in a Ricordi chamber. The digested tissues were then purified by a continuous UIC-UB gradient protocol on a cell separator (Cobe 2991, Cobe Inc., CO) and then cultured in CMRL 1066 media with 5% human albumin (CSL Behring, King of Prussia, PA) at 37°C under 5% CO₂.

Dithizone solution (DTZ.,1.5 mM, Sigma-Aldrich, MO) was used to differentiate islets and acinar tissues by staining zinc ions, which are highly concentrated in beta-cells (Ricordi et al., 1990). Tissues were incubated at 37°C for 10 min and then washed twice with PBS. Manual assessment of human islets was done under a microscope in a Petri dish labeled with measurement scales. Diameters of islets were estimated according to these scales. These diameters were used to classify islets into eight size ranges and then converted into IEq with the corresponding multipliers as described earlier in Size, Volume, and Islet Equivalents Quantification.

Human Islet Loading and Video Recording on-Chip

The counting samples of isolated human islets were picked up manually and randomly via polyethylene tubing (Intramedic, PE160) connected to a 1 ml syringe. The tubing was then connected to the inlet of the microfluidic channel. The islets were injected into the device via PHD ULTRA™ Syringe Pump (Harvard Apparatus, MA) at a flow rate of 20–30 μl/min.

The captured videos (1920 × 1,080; 1080p, 60 fps) were recorded in mp4 file format and then transferred to a computer for further video processing. The video processing algorithm was developed using python 3.8.1 within Spyder IDE (version 3.36). Video analysis was assembled for five functions: object detection, cell segmentation, cell tracking, feature extraction, and report generation. Figure 2 summarizes the procedures of object detection, cell segmentation, and cell tracking.

Statistical Analysis

Data are expressed as mean ± SD. For the processing of flowing islet videos, n = 5 experiments were conducted for each islet sample. Both manual and smartphone-based assessments were performed. Statistical significance was calculated using t-tests (p-value < 0.05).

Results and Discussion

Human Islet Quantification

In the current practice of human islet isolation, about 500 µl of cell suspension, which normally contains 0.05–0.1% of the total islet population, is picked up manually, randomly, and repeatedly for islet mass quantification (Qi et al., 2009a; b). The randomness of sample sizes often causes variation in assessment results. In this study, human islet samples with different islet masses were quantified by both manual counting and the SFDIA system for comparison. Six groups of samples were tested. Islet numbers within each sample ranged from 20 to 200. Experimental results were plotted as shown in Figure 6. Across the six groups, the average difference between the SFDIA system and manual counting was 2.91 ± 1.50%. The p-values for all six groups were greater than 0.05, indicating that there is no significant difference between the result generated by manual counting and the SFDIA system. In other words, the SFDIA system yields accurate quantification assessments close to those obtained by manual counting, regardless of islet sample sizes. However, manual counting resulted in a relatively high variation (SD = 3.1–5.0), while the SFDIA system gave consistently low variation (SD < 1). In general, for traditional manual counting, a smaller sample size creates a larger error in the results. The consistency of the SFDIA system successfully reduced human error and rendered the assessment results more reliable than those obtained from the traditional manual counting method.

FIGURE 6

FIGURE 6. Comparison between manual counting and smartphone-based counting. Six human islet samples of varied sample sizes were assessed by both manual counting and the SFDIA system. Results were plotted and expressed in mean +/− SD. NS: Not Significant (p-value > 0.05), N = 5.

Islet Size and Volume Estimation

In addition to cell numbers, we acquired and analyzed multiple parameters pertaining to flowing islets in the microfluidic channel, as discussed in Quantification of Islet Mass. Data acquisition and analysis were performed automatically during video processing.

Since islets do not have a perfect regular shape, their diameters were estimated as their equivalent diameters. The equivalent diameter is calculated from islet areas as seen in each video frame. Distribution of human islet diameters are shown in Figure 7A. The islet samples we tested had an average diameter of 99.88 ± 53.91 µm and a maximum diameter of 337.67 µm. The majority (>60%) had diameters between 50 and 150 µm. While the diameter obtained by the SFDIA system is the equivalent diameter as described above, manual assessments often use the length of the major axis as the diameter. As a result, the SFDIA diameter can be slightly smaller but more reliable than that estimated by manual assessment.

FIGURE 7

FIGURE 7. Islet diameter and volume assessment. (A) Distribution of islet diameter. The mean diameter of the islet sample is 99.87 ± 53.91 µm. N = 628. (B) Comparison between IEqV and EFV. N = 4. (C) Distribution curves of the volume calculated using the two different methods. N = 628.

We next compared IEqV and EFV. As discussed previously, IEqV is the volume calculated using the conventional IEq estimation method. EFV is calculated using the ellipsoid volume equation described by Eq. 5. As shown in Figure 7B, in samples 2 through 4, IEqV was higher than EFV, while in sample 1, the two had similar values. It is worth noting that sample 2 and 3 had very similar IEqV’s while their EFV’s differed dramatically. This mainly results from the IEqV calculation: islets are first grouped by diameter, and then the number of islets in each group is multiplied by the IEq conversion factor. As a result, IEqV tends to ignore small size variation and leads to greater error.

Two methods of islet volume estimation were compared: V_major uses the major axis as the third axis of the ellipse, while V_minor uses the minor axis. As shown in Figure 7C, our results indicated that the V_major curve was slightly to the right of that of V_minor as expected. However, the two curves largely overlapped, indicating that the two methods generated similar results in terms of islet volume estimation. As such, the choice of the third axis in a 2D-ellipse fitting model does not significantly affect islet volume estimation. Our system’s islet volume estimation better represents actual islet volume than traditional IEq-derived islet volume. Thus, our system provides more accurate information to clinicians preparing for islet transplants.

Islet Fragmentation: Circularity and Solidity Estimation

Islet fragmentation is often caused by chemical and mechanical stress during islet isolation and post-isolation culture (Kin et al., 2008). Fragmentation is a key parameter in the evaluation of islet quality. It has also been considered an important influencer of post-transplant graft survival rates. Therefore, the accurate assessment of human islet fragmentation prior to transplantation is important for improving clinical outcomes.

Circularity has been widely used to measure islet fragmentation in islet assessments (Huang H.-H. et al., 2018). An ellipse/ellipsoid model is considered to be a better estimation of islet size and volume than a circular model, since circularity alone may not be enough to quantify shape regularity. In this study, both circularity and solidity were calculated to describe islet morphology. Our human islet samples had an average circularity of 0.591 ± 0.133. As shown in Figure 8A, circularity of the majority of islets (∼83%) ranged between 0.4 and 0.7. As in Figure 8B, the solidity of the samples averaged 0.853 ± 0.107. The majority of islets (∼70%) had solidity values between 0.8 and 1. The high solidity values demonstrate a low rate of fragmentation in our preparation. The low circularity values indicate that the shapes of most islets in our samples are better represented as ellipses than perfect circles. This implies that islet solidity can be used as an improved quantification method for describing islet morphology (roundness, shape, and fragmentation).

FIGURE 8

FIGURE 8. Circularity and solidity assessment. (A) Islet circularity distribution. The islet sample has an average circularity of 0.591 ± 0.127. N = 628. (B) Islet solidity distribution. The islet sample has an average solidity of 0.853 ± 0.107. N = 628.

Islet Purity: Trapped Islet Percentage

Despite recent progress in the field of human islet transplantation, low purity of islet products still majorly affects its success rate. Islet cell percentage determines competition for oxygen and other nutrients between islet and non-islet tissues within transplanted allografts. The conventional method to estimate the purity of isolated islets involves DTZ staining and manually counting DTZ-positive islets under a light microscope. The operator-dependent nature of this process and the use of islet count as a unit can lead to significant overestimation of purity (Pisania et al., 2010; Kitzmann et al., 2014). Our software automatically analyzes the trapped islet percentage based on area data acquired from image processing. This provides a more accurate estimation of islet product purity.

As shown in Figure 9, we compared the percentage of trapped islets between two groups: middle-layer and top-layer human islets. These two layers generally vary significantly in morphology and tissue composition. In Figure 9A, the average staining ratio of the middle-layer sample was 0.576 ± 0.223. The distribution curve was skewed to the left. ∼72% of islets had DTZ positive ratios between 0.3 and 0.7. These islets had a mean circularity of 0.534 ± 0.102 and a mean solidity of 0.793 ± 0.101. In contrast, Figure 9B shows the result of top-layer islets. The average staining ratio among the top-layer islet samples was 0.735 ± 0.213, with the majority (∼67% of total islets) above 0.7. These islets had a higher mean circularity (0.636 ± 0.139) and a higher mean solidity (0.854 ± 0.107). With our video processing program, we obtained results consistent with expected values from typical islet preparations. We also produced more detailed and reliable information on islet quality than purity estimations from standard manual quantification. Our higher-quality data may prove more useful to clinicians.

FIGURE 9

FIGURE 9. Trapped islet percentage assessment. (A) Assessment for middle layer islets. The islet sample has an average DTZ⁺ ratio of 0.576 ± 0.223, an average circularity of 0.534 ± 0.102, and an average solidity of 0.793 ± 0.102. N = 198. (B) Assessment for top layer islets. The islet sample has an average DTZ⁺ ratio of 0.735 ± 0.213, an average circularity of 0.635 ± 0.139, and an average solidity of 0.853 ± 0.107. N = 530.

Conclusion

We described a novel and dynamic smartphone-based digital imaging system integrated with microfluidic technology for pancreatic islet mass quantification. Our software program can track multiple moving islets and separate closely attached islets. SFDIA generated little variance in islet counting across multiple experiments, exhibiting improved consistency as compared to the conventional manual counting method. In addition, the system estimated size using an ellipsoid model and returned additional islet parameters such as volume, circularity, solidity, and trapped islet percentage. These parameters are impossible to directly identify with the conventional method. Using our low-cost and portable system, reliable islet parameters can be easily obtained during the preparation of islet biologics.

Several limitations of this study need be mentioned. First, due to the relatively low resolution and frame rate of a smartphone camera, video processing requires a consistent fluid flow. This consistent fluid flow assures a reliable reading of islet parameters but requires an external pump system. To overcome this issue, upgrading to a higher-quality camera is favorable. With the upgrade, it is possible to simplify the system by removing the external power source for fluid delivery. This results in a pumpless microfluidic system as we introduced previously for islet function tests (Xing et al., 2016). Another limitation of our methodology is that our imaging is captured in 2D, but islet mass parameters are truly dependent on the 3D structure of the islets. Although our system is simple and has significant advantages over traditional manual counting, a 3D scanning and modeling approach can further improve accuracy.

In conclusion, manual counting tends to overestimate the islet mass, the volume, and the purity of islet samples due to human errors or method limitations. Through advanced analysis, our SFDIA program analyzes more islet quantification parameters beyond traditional islet mass quantification. As more information becomes accessible via our system, we suggest an alternative, multi-parameter islet quality assessment method. This new method would allow operators or physicians to access these islets’ properties objectively and with higher accuracy. Future applications of this smartphone-fluidic system could include a functional potency assay for a comprehensive product release test.

Data Availability Statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author Contributions

HA, YW, and YX for aid in designing the experiments. XY developed the software and performed most of the experiments. YX designed the microfluidic device. YX and PZ designed the system frame. YH performed cell preparation for the experiments. XY and YX analyzed the data. All authors contributed to discussions and manuscript preparation. XY, YX, and YW wrote the manuscript with input from all authors. EL assisted with manuscript revision. MR, HA, YW, and JO supervised the project.

Funding

This work was supported by the NIH R01 DK122253 (HA, JO, and YW), NIH NIDDK R25 DK105924-01 (JO) and UVA Launchpad for diabetes (HA, JO, and YW).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors would like to acknowledge the help from MAE Rapid Prototyping Lab and Machine Labs in fabricating the 3D-printed frame of the system.

References

Adewola, A. F., Lee, D., Harvat, T., Mohammed, J., Eddington, D. T., Oberholzer, J., et al. (2010). Microfluidic Perifusion and Imaging Device for Multi-Parametric Islet Function Assessment. Biomed. Microdevices 12 (3), 409–417. doi:10.1007/s10544-010-9398-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Buchwald, P., Bernal, A., Echeverri, F., Tamayo-Garcia, A., Linetsky, E., and Ricordi, C. (2016). Fully Automated Islet Cell Counter (ICC) for the Assessment of Islet Mass, Purity, and Size Distribution by Digital Image Analysis. Cel Transpl. 25 (10), 1747–1761. doi:10.3727/096368916X691655

PubMed Abstract | CrossRef Full Text | Google Scholar

Buchwald, P., Wang, X., Khan, A., Bernal, A., Fraker, C., Inverardi, L., et al. (2009). Quantitative Assessment of Islet Cell Products: Estimating the Accuracy of the Existing Protocol and Accounting for Islet Size Distribution. Cel Transpl. 18 (10), 1223–1235. doi:10.3727/096368909X476968

CrossRef Full Text | Google Scholar

Chowdhury, A. S., Chatterjee, R., Ghosh, M., and Ray, N. (2010). “Cell Tracking in Video Microscopy Using Bipartite Graph Matching,” in 20th International Conference on Pattern Recognition: IEEE, Istanbul, Turkey, August 23–26, 2010, 2456–2459.

CrossRef Full Text | Google Scholar

Fitzgibbon, A., Pilu, M., and Fisher, R. B. (1999). Direct Least Square Fitting of Ellipses. IEEE Trans. Pattern Anal. Machine Intell. 21 (5), 476–480. doi:10.1109/34.765658

CrossRef Full Text | Google Scholar

Gangemi, A., Salehi, P., Hatipoglu, B., Martellotto, J., Barbaro, B., Kuechle, J. B., et al. (2008). Islet Transplantation for Brittle Type 1 Diabetes: the UIC Protocol. Am. J. Transpl. 8 (6), 1250–1261. doi:10.1111/j.1600-6143.2008.02234.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Girman, P., Kříž, J., Friedmanský, J., and Saudek, F. (2003). Digital Imaging as a Possible Approach in Evaluation of Islet Yield. Cel Transpl. 12 (2), 129–133. doi:10.3727/000000003108746713

CrossRef Full Text | Google Scholar

Gmyr, V., Bonner, C., Lukowiak, B., Pawlowski, V., Dellaleau, N., Belaich, S., et al. (2015). Automated Digital Image Analysis of Islet Cell Mass Using Nikon's Inverted Eclipse Ti Microscope and Software to Improve Engraftment May Help to advance the Therapeutic Efficacy and Accessibility of Islet Transplantation across Centers. Cel Transpl. 24 (1), 1–9. doi:10.3727/096368913X667493

CrossRef Full Text | Google Scholar

Hassan, S.-u., and Zhang, X. (2020). Design and Fabrication of Capillary-Driven Flow Device for Point-Of-Care Diagnostics. Biosensors 10 (4), 39. doi:10.3390/bios10040039

PubMed Abstract | CrossRef Full Text | Google Scholar

Hering, B. J., Clarke, W. R., Bridges, N. D., Eggerman, T. L., Alejandro, R., Bellin, M. D., et al. (2016). Phase 3 Trial of Transplantation of Human Islets in Type 1 Diabetes Complicated by Severe Hypoglycemia. Dia Care 39 (7), 1230–1240. doi:10.2337/dc15-1988

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, H.-H., Harrington, S., and Stehno-Bittel, L. (2018a). The Flaws and Future of Islet Volume Measurements. Cel Transpl. 27 (7), 1017–1026. doi:10.1177/0963689718779898

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X., Xu, D., Chen, J., Liu, J., Li, Y., Song, J., et al. (2018b). Smartphone-based Analytical Biosensors. Analyst 143 (22), 5339–5351. doi:10.1039/c8an01269e

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutchison, J. R., Erikson, R. L., Sheen, A. M., Ozanich, R. M., and Kelly, R. T. (2015). Reagent-free and Portable Detection of Bacillus Anthracis Spores Using a Microfluidic Incubator and Smartphone Microscope. Analyst 140 (18), 6269–6276. doi:10.1039/c5an01304f

PubMed Abstract | CrossRef Full Text | Google Scholar

Im, H., Castro, C. M., Shao, H., Liong, M., Song, J., Pathania, D., et al. (2015). Digital Diffraction Analysis Enables Low-Cost Molecular Diagnostics on a Smartphone. Proc. Natl. Acad. Sci. USA 112 (18), 5613–5618. doi:10.1073/pnas.1501815112

PubMed Abstract | CrossRef Full Text | Google Scholar

Jennings, B., and Parslow, K. (1988). Particle Size Measurement: the Equivalent Spherical Diameter. Proc. R. Soc. Lond. A. Math. Phys. Sci. 419 (1856), 137–149.

CrossRef Full Text | Google Scholar

Kappler, J. A., Starr, C. J., Chan, D. K., Kollmar, R., and Hudspeth, A. J. (2004). A Nonsense Mutation in the Gene Encoding a Zebrafish Myosin VI Isoform Causes Defects in Hair-Cell Mechanotransduction. Proc. Natl. Acad. Sci. 101 (35), 13056–13061. doi:10.1073/pnas.0405224101

PubMed Abstract | CrossRef Full Text | Google Scholar

Kin, T., Senior, P., O’Gorman, D., Richer, B., Salam, A., and Shapiro, A. M. J. (2008). Risk Factors for Islet Loss during Culture Prior to Transplantation. Transpl. Int. 21 (11), 1029–1035. doi:10.1111/j.1432-2277.2008.00719.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kitzmann, J. P., Karatzas, T., Mueller, K. R., Avgoustiniatos, E. S., Gruessner, A. C., Balamurugan, A. N., et al. (2014). Islet Preparation Purity Is Overestimated, and Less Pure Fractions Have Lower post-culture Viability before Clinical Allotransplantation. Transplant. Proc. 46 (6), 1953–1955. doi:10.1016/j.transproceed.2014.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. A., Leitao, R., Zheng, G., Yang, S., Rodriguez, A., and Yang, C. (2011). Color Capable Sub-pixel Resolving Optofluidic Microscope and its Application to Blood Cell Imaging for Malaria Diagnosis. PLoS One 6 (10), e26127. doi:10.1371/journal.pone.0026127

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Geng, Z., Fan, Z., Liu, J., and Chen, H. (2019). Point-of-care Testing Based on Smartphone: The Current State-Of-The-Art (2017-2018). Biosens. Bioelectron. 132, 17–37. doi:10.1016/j.bios.2019.01.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, K. D., Woodburn, E. V., Le, H. M., Shah, U. K., Lumetta, S. S., and Cunningham, B. T. (2017). Multimode Smartphone Biosensing: the Transmission, Reflection, and Intensity Spectral (TRI)-analyzer. Lab. Chip 17 (19), 3246–3257. doi:10.1039/c7lc00633k

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, S., Dugan, C. E., and Kennedy, R. T. (2018). Microfluidic Chip with Integrated Electrophoretic Immunoassay for Investigating Cell-Cell Interactions. Anal. Chem. 90 (8), 5171–5178. doi:10.1021/acs.analchem.7b05304

PubMed Abstract | CrossRef Full Text | Google Scholar

McCracken, K. E., and Yoon, J.-Y. (2016). Recent Approaches for Optical Smartphone Sensing in Resource-Limited Settings: a Brief Review. Anal. Methods 8 (36), 6591–6601. doi:10.1039/c6ay01575a

CrossRef Full Text | Google Scholar

Navarro-Fontestad, C., Mangas-Sanjuán, V., González-Álvarez, I., García-Arieta, A., Fernández-Teruel, C., Casabó-Alós, V. G., et al. (2012). Computer Simulations as a Tool for Optimizing Bioequivalence Trials. Readings Adv. Pharmacokinet. - Theor. Methods Appl. 17, 15. doi:10.5772/34611

CrossRef Full Text | Google Scholar

Niclauss, N., Sgroi, A., Morel, P., Baertschiger, R., Armanet, M., Wojtusciszyn, A., et al. (2008). Computer-assisted Digital Image Analysis to Quantify the Mass and Purity of Isolated Human Islets before Transplantation. Transplantation 86 (11), 1603–1609. doi:10.1097/TP.0b013e31818f671a

PubMed Abstract | CrossRef Full Text | Google Scholar

Nourmohammadzadeh, M., Xing, Y., Lee, J. W., Bochenek, M. A., Mendoza-Elias, J. E., McGarrigle, J. J., et al. (2016). A Microfluidic Array for Real-Time Live-Cell Imaging of Human and Rodent Pancreatic Islets. Lab. Chip 16 (8), 1466–1472. doi:10.1039/c5lc01173f

PubMed Abstract | CrossRef Full Text | Google Scholar

Olehnik, S. K., Fowler, J. L., Avramovich, G., and Hara, M. (2017). Quantitative Analysis of Intra- and Inter-individual Variability of Human Beta-Cell Mass. Sci. Rep. 7 (1), 16398. doi:10.1038/s41598-017-16300-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Pisania, A., Weir, G. C., O'Neil, J. J., Omer, A., Tchipashvili, V., Lei, J., et al. (2010). Quantitative Analysis of Cell Composition and Purity of Human Pancreatic Islet Preparations. Lab. Invest. 90 (11), 1661–1675. doi:10.1038/labinvest.2010.124

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, M., Barbaro, B., Wang, S., Wang, Y., Hansen, M., and Oberholzer, J. (2009a). Human Pancreatic Islet Isolation: Part I: Digestion and Collection of Pancreatic Tissue. JoVE 27, 1125. doi:10.3791/1125

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, M., Barbaro, B., Wang, S., Wang, Y., Hansen, M., and Oberholzer, J. (2009b). Human Pancreatic Islet Isolation: Part II: Purification and Culture of Human Islets. JoVE 27, 1343. doi:10.3791/1343

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, M., Kinzer, K., Danielson, K. K., Martellotto, J., Barbaro, B., Wang, Y., et al. (2014). Five-year Follow-Up of Patients with Type 1 Diabetes Transplanted with Allogeneic Islets: the UIC Experience. Acta Diabetol. 51 (5), 833–843. doi:10.1007/s00592-014-0627-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ricordi, C., Goldstein, J. S., Balamurugan, A. N., Szot, G. L., Kin, T., Liu, C., et al. (2016). National Institutes of Health-Sponsored Clinical Islet Transplantation Consortium Phase 3 Trial: Manufacture of a Complex Cellular Product at Eight Processing Facilities. Diabetes 65 (11), 3418–3428. doi:10.2337/db16-0234

PubMed Abstract | CrossRef Full Text | Google Scholar

Ricordi, C., Gray, D. W. R., Hering, B. J., Kaufman, D. B., Warnock, G. L., Kneteman, N. M., et al. (1990). Islet Isolation Assessment in Man and Large Animals. Acta Diabet. Lat 27 (3), 185–195. doi:10.1007/BF02581331

PubMed Abstract | CrossRef Full Text | Google Scholar

Shapiro, A. M. J., Lakey, J. R. T., Ryan, E. A., Korbutt, G. S., Toth, E., Warnock, G. L., et al. (2000). Islet Transplantation in Seven Patients with Type 1 Diabetes Mellitus Using a Glucocorticoid-free Immunosuppressive Regimen. N. Engl. J. Med. 343 (4), 230–238. doi:10.1056/NEJM200007273430401

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L.-J., and Kaufman, D. B. (2016). Digital Image Analysis to Assess Quantity and Morphological Quality of Isolated Pancreatic Islets. Cel Transpl. 25 (7), 1219–1225. doi:10.3727/096368915X689947

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Danielson, K. K., Ropski, A., Harvat, T., Barbaro, B., Paushter, D., et al. (2013). Systematic Analysis of Donor and Isolation Factor's Impact on Human Islet Yield and Size Distribution. Cel Transpl. 22 (12), 2323–2333. doi:10.3727/096368912X662417

CrossRef Full Text | Google Scholar

Xing, Y., Nourmohammadzadeh, M., Elias, J. E. M., Chan, M., Chen, Z., McGarrigle, J. J., et al. (2016). A Pumpless Microfluidic Device Driven by Surface Tension for Pancreatic Islet Analysis. Biomed. Microdevices 18 (5), 80. doi:10.1007/s10544-016-0109-4

PubMed Abstract | CrossRef Full Text | Google Scholar