Monitoring Human Neutrophil Differentiation by Digital Holographic Microscopy

We report on the usefulness of digital holographic microscopy (DHM) for the assessment of human neutrophil differentiation from myeloid cells. The cell and nuclear regions have been designated by image segmentation of the optical phase function, and the changes of the cell nucleus morphology in relation to the whole cell morphology have been examined during the process of granulocytic differentiation into mature neutrophils in PLB-985 cell line. Nucleus phase volume and circularity and the ratios between the nucleus and the cell projected area and volume provide a reliable set of parameters to characterize the maturation process. As control, cell differentiation has been monitored in parallel using standard nucleus staining and fluorescence imaging. From this research, it emerged that DHM can be used as a valid label-free solution, alternatively to the standard staining technology.


INTRODUCTION
Human neutrophils are central players in innate immunity, a major component of inflammatory responses, and a leading model for cell motility and chemotaxis. The neutrophils are generated from granulocyte-monocyte progenitors (GMPs) in the bone marrow, where they undergo several stages of maturation [1]. Primary, mature blood neutrophils have a short lifetime and lack proliferation and transfection capacity, thus limiting their ex vivo experimental use in the laboratory [2]. Therefore, protocols using human myeloid leukemia cell lines, such as HL-60 and PLB-985 (which provides higher differentiation efficiency), have been developed for neutrophil differentiation in vivo [3,4]. In the presence of appropriate inducing agents, these cell lines can undergo granulocyte differentiation into mature neutrophil-like granulocytes [1,2,5]. The assessment of the differentiation process is based on cellular morphological parameters, the molecular and functional properties of differentiated cells, defining a neutrophil-like phenotype which closely resembles the mature blood neutrophils [5,6].
Non-differentiated cells are considered as being at the promyelocytic stage. Activation by chemical agents transforms them along the myeloid pathway into myelocytes (Ms), metamyelocytes (MMs), band neutrophils (BNs) and, lastly, segmented neutrophils (SNs) [1,2,7]. This transformation process takes about 4-6 days. Monitoring the changes in nuclear morphology represents one of the most reliable technique to assess different stages and the success of the differentiation process [6,8,9]. Thus, cells in SN stage, presenting polymorph nucleus formed usually by two to five nuclear segments, are considered as mature neutrophil-like cells. In the first two stages [promyelocyte (PM) and M], cells have compact ellipsoidal nuclei occupying almost all of the cell body. Nuclei of the MMs begin to be distorted and present indentations. BNs have more distorted nuclei, often with a "horseshoe" shape, while in the final stage, SNs present multilobed nuclei connected between them [6,8,9].
In biological laboratories, the current practice to examine the morphology of the nucleus during neutrophil differentiation involves nuclear staining and epifluorescence microscopy [7,[10][11][12]. Since these techniques require sample preparation and provide relatively poor information, they are mostly used to evaluate the final stage of the differentiation process. Moreover, the labeling methods have a potential risk of altering cellular functions and phototoxicity from repeated excitation of the dye in the imaging process, which limits the duration of the time-lapse imaging [12]. The use of non-quantitative phase microscopy, which has become a pervasive tool for contrast enhancement and qualitative examination of cellular morphology [13], is limited by the low contrast enhancement for these cells.
Label-free quantitative phase imaging (QPI) [14][15][16], including digital holographic microscopy (DHM) [17][18][19], has been recently implemented in different setups to study morphology, structure, and dynamics of unstained living cells. QPI is based on the measurement of the optical phase (OP) shift introduced by the sample cell, which allows to calculate the optical path difference (OPD): OPD = OP · 2π/λ, where λ is the wavelength of the laser beam. The OPD combines the physical thickness, h, and the difference between the refractive index (RI) of the cell, n c , and the RI of the medium, n m : OPD = h · n, where n = (n c − n m ). The RI of the cell represents an important parameter which can be used to correlate with other cell biophysical parameters, such as dry and wet mass, protein concentration, elasticity, and conductivity, and to study certain cell metabolic activities, such as cell division and infection [20]. To isolate the RI from the thickness in OPD, several solutions are possible [21]. Assuming that the RI of the suspension medium is known and that cells in suspension are spherical, the simplest and fastest method to determine the integral RI is to approximate the local thickness [22]. This method has been used to show changes in the subcellular structure and the RI of differentiating myeloid precursor cells within 1 day of differentiation induction [23] and to demonstrate that cell nuclei have lower RI than cytoplasm [24]. Another approach is performing two OPD measurements, each of them with a different surrounding medium [25,26] or a different wavelength [27], thus obtaining two linear equations with two unknowns which enable the decoupling of the integral RI from the thickness. Other solutions for the thickness-RI coupling problem are thoroughly reviewed in Dardikman and Shaked [21].
Nonetheless, various cell morphological features can also be directly extracted from the OPD, e.g., OPD mean/median, phase volume, phase surface area, dry mass, phase sphericity, phase statistical parameters, and energy [28]. Such features have been successfully combined with machine learning algorithms in cytometric classification of cancer cells and blood cells [29], phenotyping of cell lines [30], classification of leukocytes flowing in microfluidics [31], or automatic detection of Plasmodium falciparum [32]. Combining QPI with dye exclusion cell volumetric imaging, the cell volume can be quantified independently of the RI of the cell, enabling high-throughput neutrophil differentiation from other white blood cells flowing in an optofluidic device [33].
Using QPI with a common arrangement for DHM, in this study, we investigate on the usefulness of a set of cell parameters extracted directly from OPD to supervise neutrophil differentiation from the PLB-985 cell line. Since the differentiation protocol includes five stages, namely, PMs, Ms, MMs, BNs, and SNs, and lasts for about 5-6 days, a set of OPD images were taken every day and analyzed. The reconstructed OPD functions for a set of cells in full field were first segmented to separate the cells from the background and then a second segmentation was applied at the single cell level to demarcate the central region with high phase values, designating to the nuclear region. The OPD and a set of nine morphological parameters were calculated for each cell.
Comparing the values of the parameters calculated for cells from different stages, we discuss their relevance and show that DHM can be successfully employed to distinguish between different cell stages of neutrophil differentiation.

Cell Preparation
Neutrophil differentiation from its promyelocytic progenitors was performed following the dimethyl sulfoxide (DMSO) protocol described in Tucker et al. [4]. Briefly, PLB-985 cells (DSMZ, Braunschweig, Germany) were cultured in RPMI 1640 media supplemented with 10% fetal bovine serum (FBS), 1% L-Glutamin, and 100 mg/mL and 100 U/mL (P/S) of streptomycin and penicillin, respectively, at 37 • C in a humidified atmosphere with 5% CO 2 . Cell cultures were passed three times per week, maintaining cell densities between 10 5 and 2 × 10 6 cells per mL. The cells were differentiated into a neutrophil-like state by culturing at an initial density of 2×10 5 using RPMI 1640 media supplemented with 10% FBS concentration and 1.3% DMSO for 6 days.
Cell differentiation was monitored every 24 h for 6 days beginning from day 0. For each measurement, 1 mL of cell suspension was taken from the culture flask. The suspension was then centrifuged, and the supernatant was removed to exclude debris and impurities. Then, the sample was resuspended and incubated for 5 min at 37 • C. To measure the cells, 50 µL of cell suspension was placed in a 18 mm diameter glass slide. About 2-3 min was required for the cells to sediment before imaging it. For DHM imaging, the sample was resuspended in RPMI 1640 media, while for fluorescence imaging, the sample was resuspended in a Hoechst 33342 solution of 1:1,000 in RPMI 1640 media as described in the study by Chazotte [11].

Digital Holographic Microscopy
We used a simplified version of a previously developed DHM off-axis configuration, based on a Mach-Zehnder interferometer [34,35]. The laser beam (λ = 520 nm, Thorlabs LP520-SF15) was split into object and reference beams by a fiber optic coupler (Thorlabs TW560R3F1) and recombined by a cube beam splitter (Thorlabs, BS079) to generate the hologram, which was recorded on a sCMOS camera (Thorlabs CS2100M-USB Quantalux). The power P of the laser was set to P < 1 mW, and the exposure time to 1 < t < 5 ms. The magnification of the microscope was 33×, with a lateral spatial resolution d ∼ 600 nm. An aspheric lens (Thorlabs C230TMD-A) with the numerical aperture NA = 0.55 was used as the objective lens. The size of the sCMOS sensor was 4.8 × 4.8 µm. The inter-fringe, i, of the interference pattern was adjusted to five pixels on the sensor, corresponding to 750 nm at the sample plane.
The numerical reconstruction of the optical phase was performed using a method based on the spatial filtering at the Fourier plane from the recorded hologram [36]. Thus, the complex amplitude distribution of the transmitted frequency band-pass was recovered by applying a Fourier transform over the recorded hologram while considering a circular spatial filtering mask located at one of the diffraction orders. After the filtering and centering process at the Fourier domain, the complex amplitude was Fourier transformed again to propagate it to the sample plane and retrieve the wrapped OP. The final step, phase unwrapping, is used to get the reconstructed OP function. Since these issues are of particular significance in the reconstruction process, we defined the parameters for spatial filtering, propagation, and phase unwrapping by calibration with silica microbeads of different diameters (3 − 30µm) [34]. Microbeads are the perfect sample for calibration purposes since they are microspheres of a single material (silica, polystyerene, etc.) with known value of the RI.

Segmentation of the OPD Function for Cell and Nuclear Region Designation
In the first step, the OPD map was segmented to isolate the cells from the background by applying the bimodal Otsu's method [37,38] implemented in MATLAB, using the following two functions, multitresh (L = 2 levels) and imquantize, from the Image Processing Toolbox. To remove the incongruous cell spots resulting from segmentation, the phase volume of each cell spot was calculated and those spots with the phase volume below a predetermined threshold were discarded.
The designation of the nuclear region was accomplished in a second segmentation applied to each cell, using multitresh with (L = 4). The number of segmentation levels is supported by the Histogram-based Valley Estimation Method presented in the study by Huang et al. [39], where the number of clusters for an image to be properly segmented depends on the number of valleys in the histogram of the image. After segmentation, the first level is assigned to the background pixels, the second and the third levels to the cell cytoplasm region, and the fourth level to cell nucleus. The segmentation method is described in detail in Supplementary Material.

RESULTS AND DISCUSSION
In order to assess Optical phase (OP) imaging by DHM as a freelabel technique for characterization and discrimination of the cell stages, we monitored the differentiation of PLB-985 cell line from PMs to SNs and imaged the cells at intervals of 24 h for 6 days during the neutrophil differentiation cycle.
As described in section 2.3, segmentation of the reconstructed OP for an image containing several cells was used to select the cells from the background and delineate the estimated nucleus from the other part of the cell. An example of the image taken at day 5 is illustrated in Figure 1D, showing the presence of multiple types of cells as indicated by the morphology of the cell and the designated nuclear regions, characterized by high values of the OPD function. At a first qualitative inspection, using the morphological phenotypes reported in literature and as already mentioned in the "Introduction" section, we can recognize one SN, three BNs, two MMs, one M. BNs are the most frequent, as expected for day 5. The presence of other types of cells is justified by the various biological times employed by different Frontiers in Physics | www.frontiersin.org FIGURE 3 | is indicated by blue color. The legend of the parameters is shown top right: PSA N , PSA C -nucleus PSA and cell PSA, respectively; C N , C C -nucleus circularity and cell circularity, respectively; OV N , OV C -nucleus OV and cell OV, respectively; PSA N/C , OV N/C -ratio between the cell and nucleus for PSA and OV; SNG-segmented nucleus. cells to transform [1]. This observation suggests that OPD and the segmented OPD could be used to characterize different stages during neutrophil differentiation.
Since the results reported in the literature are based on nuclear staining and fluorescence imaging, we also performed image analysis of the nucleus morphological changes using Hoechst staining and live cell fluorescence imaging (Nikon TE2003-E inverted microscope, 60X, DAPI fluorescence cube). An example of nuclear morphology changes during the five stages of the neutrophil differentiation is shown in Figure 2. The first two rows are brightfield and DIC images, which help to define the cell shape but do not allow to discriminate the nucleus. Complementarily, the fluorescence images illustrated in the third row for the same cell allow the identification of the nucleus. At the beginning of the differentiation, corresponding to PM and M cells, the cell nucleus changes round to oval, occupying more than 50% of the cell projected surface area. The oval shape begins to be preponderant to round shape for M, and the ratio between nucleus and cell area is also slightly smaller than Ms for PMs. As the cell matures, the nucleus begins to be indented and becomes asymmetric with respect to the cell center (MM). A more mature cell presents highly indented nucleus and often a "horseshoes" shape (BN), arriving to the segmentation of the nucleus at the end (SN).
Free-label DHM imaging for the five stages of the cell differentiation is exemplified in the last two rows of Figure 2. The OPD and the segmented OPD for an isolated cell from each of the five stages are represented, confirming the observations acquired from brightfield and fluorescence imaging on the cell and nucleus configuration during the differentiation cycle. Moreover, the OPD values in the cell and nuclear regions assigned by segmentation provide additional information on the cell morphological features as the OP volume.
We monitored the neutrophil differentiation recording five DHM image fields/day at intervals of 24 h starting from cell plating. Two different differentiation cycles were completed. The OPD was calculated for each field, and a qualitative inspection was performed as discussed above and showed in Figure 1, discarding the cells which did not correspond to the differentiation stage of the respective day. We created then five groups of N > 20 cells, corresponding to the five stages: PM, M, MM, BN, and SN, and calculated the parameters defined above for each cell.
Cellular morphologic parameters as projected surface area (PSA), OP volume (OV), and circularity (C) have been effectively used for the characterization of red blood cells [35,40,41] in QPI and DHM. In this study, we propose similar parameters for the nucleus, including the ratios between the values for nucleus and cell as parameters. We define the following morphological parameters: 1. Projected surface area (PSA) of the cell and nuclear regions: where N is the number of pixels within the cell nucleus region, ps = 4.8 * 4.8 = 23.04 µm 2 is the area of a single pixel of the CMOS sensor, and M = 33× is the microscope magnification. 2. Optical Phase Volume (OV) of the cell and nuclear regions: where OPDi is the optical phase difference of cell/nucleus corresponding to the pixel. 3. Circularity (C) of the cell and nuclear regions: where Pm is the perimeter of the cell/nucleus. Circularity expresses the roundness of an object, e.g., C = 1 means a perfectly round cell, as a circular disk. 4. PSA and OV nucleus to cell, N/C, aspect ratio: where N stays for nucleus, C for cell, PSA for projected surface area, and OV for optical volume.
The mean and SD values for each parameter are presented in Figure 3A for each group of cells. This study shows that PSA and OV monotonously decrease as the cells maturate, both for the cell and the nucleus region. The cell C, C C , is preserved close to C C = 1 for PM, N, and MM, while decreasing slightly for the last two stages, BN and SN (C C > 0.9). As the cells present more protrusions in these two stages, the values for C are also more dispersed. The C values for nucleus, C N , decrease more than C N > 0.8 for PM and M, C N = 0.8 ± 0.08 for MM, C N = 0.69±0.13 for BN and C N = 0.63±0.15 for SN, indicating that the nucleus undergoes substantial morphological changes during the differentiation process. The relative morphological changes of the nucleus vs. cell are outlined by the N/C ratio for area and volume, respectively, PSA N/C and OV N/C , in the last row of Figure 3A. Both ratios decrease monotonically during cell maturation, with a more pronounced change in OV N/C . We note that the values between neighboring stages overlap for some parameters, indicating that a set of parameters should be considered rather than a single parameter to characterize the differentiation stages. In order to test the weight of each parameter, we used the Mann-Whitney U test (or two-sided Wilcoxon rank sum test) for each pair of classes. The results are schematically illustrated in Figure 3B: p < 0.05 (blue box) is considered as a threshold for good distinction between stages, 0.05 < p < 0.1 (gray) as acceptable, and p > 0.1 (white) as non-confident.
It can be seen that all eight parameters are individually suitable to distinguish SN from PM. All parameters except cell PSA C and cell OV C are useful to discriminate SN from M and MM. Only three parameters (OV N/C , PSA N/C , OV N/C ) are available to distinguish between SN and BN. Note that all of them were extracted from the OPD, which is provided by DHM only. Additionally, the presence of a segmented nucleus, which is confidently detected from OPD, represents another useful parameter which distinguishes SN from all the other stages (crossed box in Figure 3B).
All eight parameters are suitable to distinguish BN from PM cells, four of the eight parameters to separate BN from M cells, and only two parameters to discriminate BN from neighboring stage, MM: the nucleus and cell circularity, C N and C C . We notice that these two parameters are complementary to the parameters which distinguish BN from the other neighboring stage, SN. Moreover, considering 0.05 < p < 0.1 for an acceptable discrimination, the nucleus area, PSA N , becomes an eligible parameter to distinguish between BN and MM or M.
As regards the MM cells, seven of the eight parameters are useful to distinguish MM from PM cells, but none is adequate using the criterium p < 0.05 to separate MM from M. However, considering p < 0.1, both NA and NC can be used to separate M from MM. Finally, five of eight parameters can be used to separate M from PM cells. These parameters are complementary to the parameters used to distinguish between BN and MM, and only two parameters are in common with the discrimination between BN and SN.
Five parameters can be used to separate M from PM cells. These parameters are complementary to the parameters used to distinguish between BN and MM, and only two parameters are in common with the discrimination between BN and SN.
According to the imposed criterium, p < 0.05, the nearest two classes seem to be MM and M cells. However, by relaxing the criterium to p < 0.1, two more parameters are indicated as suitable for MM and M cells discrimination.
The analysis of the results show that all these parameters are useful to discriminate between cells in different differentiation stages. Three of these parameters, OV N , OV C , OV N/C , which have a central role for cell characterization and discrimination between differentiation stages, are extracted from the OPD, which is specific for DHM. These preliminary results suggest that a combination of the parameters derived by means of DHM in a multiparameter testing scheme would enable the characterization and monitoring of the different stages of cell differentiation.
We note that the OPD function was used directly, without decoupling the RI from the cell height. Although, the designation of the nuclear region by segmentation cannot be claimed as precise identification of the nucleus, the method appears useful to monitor the neutrophil differentiation cycle. Considering the recent findings showing that the nuclear RI is lower than the cytoplasm [24,42], this affects the OV parameter. We calculated the OV of the value which corresponds to a physical volume bigger for the nucleus relative to the same value of the OV for the cell. Although the absolute values of the volume change, the relative differences between stages are not expected to be altered considerably because the changes are uniformly applied.
Another interesting point is the possible change of the cellular RI during the HL60 cells differentiation [23], indicating that the cells become less dense during differentiation. Since we measured only the OPD, we could not confirm that this observation also applied for PLB-985 cells cultured on substrate, but we observed a significant decrease of the PSA N/C , which was associated with the decrease in the RI in this study.

CONCLUSION
In this study, we performed cell morphological analysis during promyelocytic differentiation into SN for both cell and nuclear regions using DHM as free-label optical imaging technique. We defined a set of eight morphological parameters and investigated their contribution to characterize cells belonging to one of the five maturation stages: PMs, Ms, MMs, BNs, and SNs.
We proved that DHM provides additional information with respect to the standard fluorescence imaging of stained nucleus. This information is related to the volumes of cells and nucleus and allows a better discrimination between different stages of neutrophil differentiation. We believe that DHM has been employed for the first time for monitoring neutrophil differentiation in PLB 945 cells. The results suggest DHM is a suitable free-label technique providing better cell characterization and discrimination between differentiation stages.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
DC and ML conceived the project, analyzed and interpreted the data, and wrote the manuscript. DC designed and implemented the DHM. ML performed the experiments. All authors contributed to the article and approved the submitted version.