Monocytic lineage cells (monocytes, macrophages, and microglia) survey tissue for pathological threats, and act as hubs of the innate immune system. They exercise innate immune functions in numerous ways: by secreting cytokines and chemokines, acting through complement signaling pathways, and phagocytosing pathogens and debris. The extent of macrophage activation can tip the balance of immune response within the body during health and disease. In diseases such as cancer, the activation of monocytic cells can result in either pro or anti tumoral responses, which trigger different disease outcomes. Similarly, monocytic cells modulate the progression and remission of autoimmune and inflammatory diseases. It is therefore critical to characterize these cells and their responses, to identify steps which may be exploited therapeutically.
In vitro systems interrogate the functions of these innate immune cells in response to stimuli via specific assays. The phagocytosis assay is an essential tool in characterizing the magnitude of monocytic cell activation. Several variations, such as pulsechase assay, particle loading, and/or bead loading, exist. Upon stimulation of the monocytic cell population with drugs, cytokines, or preconditioned cell media, the results of the assay report to what extent the cells have increased or decreased their ability to engulf particles from the environment, thus giving insight into how such stimuli may affect monocytic cells during health and disease in vivo.
Here we discuss some of the technical challenges associated with deriving highquality data from standard phagocytosis assays. The variability of phagocytic activity of cells within a single sample makes this a particularly difficult function to assess. Most current studies count phagocytosed particles per cell, however, individual particles can be difficult to discern, and this measure can be misleading in cells that have higher baseline phagocytosis. Other studies simply analyze how many cells are phagocytosing (versus not), which gives a very broad and less quantitative picture of phagocytosis. Here, we have developed and propose the use of a new technical pipeline to collect and more sensitively analyze phagocytosis data. We then apply it to two different experiments for which macrophage phagocytosis is often assayed.
Materials and Reagents

Glass coverslips (Fisher Scientific, catalog number: 12545100 )

40 μm cell strainer (Fisher, catalog number: 22363547 )

Adult C57Bl/6 mice

L929 cells [or 100 U/ml recombinant MCSF (Prospec, catalog number: CYT439 )]

GL261 murine glioma cells

2.5% avertin

Hank’s Buffered Saline Solution (HBSS)

Fluorescent latex beads (carboxylatemodified polystyrene, fluorescent red; 0.5μm mean particle size, Sigma, catalog number: L3280 )

Fluorescenceconjugated AlexaFluor 488 goat antirabbit (Life Technologies, catalog number: A11008 )

Hank’s Buffered Saline Solution (HBSS) (Corning, catalog number: MT20023CV )

DMEM (Corning, catalog number: 10017CM )

FBS (VWR, catalog number: 10221866 )

AntibioticAntimycotic (Fisher, catalog number: 10221866 )

Sodium pyruvate (VWR, catalog number: 13115E )

PFA (4% paraformaldehyde) (Acros, catalog number: AC169650010 )

Triton X100 (Sigma, catalog number: X100 )

Normal goat serum (Southern Biotech, catalog number: 006001 )

BSA (VWR, catalog number: RLBSA50 )

DAPI Fluormount (Southern Biotech, catalog number: 010020 )

Iba1 antibody (Wako USA, catalog number: 1919741 )

LPS (Sigma, catalog number: L6511 )

Mouse IFNγ (Roche, catalog number: 11276905001 )
Equipment

Leica SP8X confocal imaging system
Beads (λ_{ex} ~575 nm; λ_{em} ~610 nm)
AlexaFluor 488 (λ_{ex} 490 nm; λ_{em} ~525 nm)

80 °C freezer
Software

Fiji (Schindelin et al., 2012)
Procedure

Cell Culture
Note: The cells, unless otherwise noted, are grown in 10% FBS with 1% AntibioticAntimycotic at 5% CO_{2} and 37 °C.

Generate glioma conditioned media (GCM) using GL261 murine glioma cells which are grown to 8090% confluency and serum starved for 24 h. The media are collected and stored at 80 °C until use.

Generate L cell media by collecting it from L929 cells grown to 90100% confluency and maintained for 10 days. These media are rich in colony stimulating factor 1 (CSF1), a factor used to drive proliferation and differentiation of monocytic cells into macrophages.

Monocytes are cultured from adult C57Bl/6 mice as previously described (Zhang et al., 2008): Briefly, animals are deeply anesthetized with 2.5% avertin and euthanized by cervical dislocation. Femurs and tibiae are dissected, and the internal bone cavity is flushed with 1% serum in Hank’s Buffered Saline Solution (HBSS). The collected bone marrow is filtered through a 40μm cell strainer and plated in complete media containing DMEM, 10% FBS, 1% AntibioticAntimycotic, 1% sodium pyruvate, and 20% L cell media (or 100 U/ml recombinant MCSF).

Cells are kept at 5% CO_{2} and 37 °C. Macrophage differentiation is achieved after 79 days of culture.

Phagocytosis assay
Experiment 1: Macrophages are plated on glass coverslips and incubated with either complete media, as described above, or complete media plus 50% GCM for 24 h.
Experiment 2: Macrophages, as above, are incubated with saline (control) or both 100 ng/ml LPS (Sigma) and 100 U/ml IFNγ (Roche) for 24 h.

Fluorescent latex beads (carboxylatemodified polystyrene, fluorescent red; 0.5μm mean particle size, Sigma, L3280 ) are added to a final concentration of 0.1 μg/ml, and incubated for 2 h. After rinsing with HBSS to remove excess beads, cells are fixed in 4% PFA at room temperature for 20 min.

Fixed cells are rinsed with PBS and blocked for 1 h in 0.1% Triton X100 and 5% normal goat serum. They are incubated with antiIba1 (1:500, Wako) in 0.01% Triton X100 and 30 mg/ml BSA overnight at 4 °C.

The cells are then rinsed with PBS, and incubated with fluorescenceconjugated AlexaFluor 488 goat antirabbit for 1 h.

Cells are rinsed again with PBS and coverslipped using DAPI Fluormount (Southern Biotech).

Confocal Imaging
Images are acquired on a Leica SP8X confocal imaging system. 10 µm zstacks are taken at 40x at a resolution of 512 x 512 pixels. Confocal images are processed using Fiji (Schindelin et al., 2012). For this procedure, five fields are imaged, and at least 100 cells are quantified in total per treatment. For evaluation of phagocytosis, total number of cells, and number of cells with internalized beads are counted. The images are further analyzed by counting the number of beads per cell, and measuring cell area and the integrated density of sumprojected pixels from the fluorescent bead channel. Integrated density values are analyzed for outliers before downstream analysis. The stepbystep image analysis protocol is provided in Supplementary document 1.
Data analysis
Data analysis is performed in R for Windows. The method and interactive UI can be accessed via GitHub for users to upload and analyze data (RShinyApp; GitHub Repo).
Results
A high throughput method of data collection
In the first experiment, we quantified macrophage phagocytosis in response to GCM. Tumor associated macrophages are known to be highly activated and protumorigenic. Their phagocytic activity may be related to their tumor supporting behavior by mechanisms that promote wound healing and antiinflammatory phenotypes. Confocal imaging allows complete 3D capture of macrophages engulfing fluorescent beads (Figure 1A). Standard methods attempt to count the number of beads that each cell has engulfed. However, some cells exhibit low, basal level of phagocytic activity, while others become stimulated and intake more beads. This nonstandardized phagocytic activity can be problematic for quantification and analysis. The beads overlap within and across zstack layers, making distinguishing and counting individual beads difficult. Further, counting the number of beads per cell in a single field of view can take more than 15 min, making data acquisition timeconsuming for large numbers of technical replicates. This highlights the technical challenges associated with this method of quantification.
To address this, others have quantified the number of cells positive or negative for bead signal, thereby creating a ratio of cells which are phagocytic in response to treatment (Lian et al., 2016; Roy et al., 2018). While more timeefficient, this method does not capture the changes in relative amount of phagocytosis occurring, i.e., phagocytic activity of each cell. As seen in Figure1A, cells engulf variable amounts of beads. Other methods to access phagocytosis, such as fluorescentlyactivated cell sorting (FACS), have also been used (Pul et al., 2013), however they require more preparation time, experimental controls (negative fluorescent controls, compensation techniques etc.), and several costly reagents and conjugated antibodies.
We developed a technical pipeline which takes advantage of the discrete fluorescent signal produced by the latex beads. By summing the fluorescent pixels across the zstack of the bead channel, we create a single metric. By using region of interest (ROI) functionality, we can capture this metric for each cell within a field of view. This method of data collection allows the analysis of all cells within a field in under 2 min, appreciably faster than counting individual beads/cell. It allows hundreds of cells to be analyzed in a single experiment.
To validate this fluorescent metric, we compared the signal produced by this method with the counted number of beads per cell, as previously performed (Figure 1B). Simple linear regression analysis was used. The model yielded high significance (P < 0.001) with a slope coefficient of 0.53 and an adjusted r^{2} = 0.363. This relationship was also confirmed with a Spearman’s rho value of 0.683. Thus, generating a fluorescent signal by summing the pixels throughout the zstack is an accurate and sensitive measurement, reflecting the number of beads each cell has engulfed. It provides a high throughput alternative to existing quantification methods.
Figure 1. Confocal imaging and analysis of macrophages following incubation with fluorescent latex beads. A. Macrophages stained for Iba1 (green) can be seen engulfing fluorescent latex beads (red) in both treatment conditions. Orthogonal views are also shown, demonstrating complete uptake of the beads by the cells. Note the varying levels of bead signal among cells. DAPI marks nuclei. B. Scatter plot comparing the fluorescence signal derived from beads versus the actual number of beads counted per cell. y = 22.40619 + 0.52782x. P < 2.2e^{16}. Adjusted r^{2} = 0.3634. Spearman’s = 0.683. C. Histogram and density plots of hypothetical datasets. For both populations dist1 and dist2, n = 100, mean = 50, median = 50. D. Histogram and density plots of experimentally derived data from panel A, Experiment 1. Complete descriptive statistics are in Table S1. Scale bar for panel A is 20μ.
Utilizing empirical cumulative distribution functions as a measure of phagocytic activity
The data rendered by this method create a continuous distribution of fluorescent signal across all cells within an experimental group (Figure 1D). Standard statistical methods would use a parametric or nonparametric statistical test to compare the means or medians of each distribution, respectively. However, let us consider two hypothetical distributions, for which the probability distribution functions (PDF) are noncongruent, but the mean and median between them are equal (Figure 1C, Table S1). Both parametric and nonparametric tests reject the alternative hypothesis and infer no difference between the two distributions, yet it is clear visually that the distribution of events is vastly different between the two samples. This suggests that for nonparametric, multimodal distributions the standard ttest or MannWhitney U test is insufficient, and can break down under certain conditions. Interpreted biologically, these tests may lack the ability to detect the ‘shifts’ in phagocytic influx that occur when macrophages are challenged, for example by treatment with GCM.
If we consider the empirical cumulative distribution function (ECDF) of these distributions, it more aptly represents the changes that exist as the PDF expands (Figure 2A). To compare these distributions, we applied principals from the field of optimal mass transport, otherwise known as Earth mover’s theory. These methods seek to compare data based on their underlying ECDFs. To compute differences between the ECDF derived from our hypothetical and experimental data, we next performed the KolmogorovSmirnov test (KS test), where the D statistic given by:
Where F_{1,n} and F_{2,m} are the ECDF from two 1 dimensional datasets, and sup is the supremum function, computes the absolute maximum distance between the ECDFs (Figure 2A). To reject the null hypothesis Ho, which states that both samples come from the same population, at a given significance levelα, D_{n,m} > D_{n,m,α}, where D_{n,m,α} is the critical value given by
For the hypothetical distributions dist1 and dist2 the two sample KS test yields D = 0.5 and P = 2.778e^{11}. This clearly indicates the advantage of using this method for multimodal distributions over standard ttest or MannWhitney U test. When applied to our experimental data, KS results in D = 0.24 and P = 2.2e^{4} (Figure 2B), indicating that the probability that a macrophage has increased phagocytic activity is more accurately reflected in the collected data by the shift in its PDF, which is captured by analyzing the ECDF. These results also demonstrate that exposure to GCM significantly increases macrophage phagocytic activity.
Comparing KS test to reported methods
To further assess KS test sensitivity, we performed a second experiment exposing macrophages to the canonical proinflammatory/M1 activating agents LPS and IFNγ. We calculated the percentage of phagocytosing cells by dividing the number of cells positive for engulfed beads by the total number of cells within the field of view, as reported (Lian et al., 2016; Roy et al., 2018). The results yielded a trend towards increase in phagocytosis in cells exposed to LPS and IFNγ, however the high variability and lack of technical replicates fail to reject the null hypothesis using standard two sample parametric and nonparametric tests (Figure 2C). We then employed the KS test to generate PDF and ECDF of the two treatment groups (Figure 2D). This resulted in D = 0.15 and P = 0.04, thus detecting a small but significant shift in phagocytic activity, which can be seen by the rightward shift in the histogram, density curve, and ECDF. While not as large of an increase relative to that observed after GCM treatment, this result agrees with published literature (Wu et al., 2009; Scheraga et al., 2016; Caponegro et al., 2019) and validates the approach, highlighting its high throughput and sensitivity.
Figure 2. Comparison of quantification of phagocytosis using traditional metrics vs the proposed new method of analysis. A. Empirical cumulative distribution functions (ECDFs) of nonparametric, multimodal hypothetical datasets dist1 and dist2. The ECDF more accurately portrays the changes that occur in the underlying PDF. B. ECDFs of data experimentally derived from Figure 1A, Experiment 1. The D statistic is plotted, showing the maximal distance between the curves. C. Using traditional metrics to quantify phagocytosis in vitro, cells that were positive for bead engulfment were divided by the total number of cells per field of view. Data are derived from Experiment 2. D. Top panel: Using the new methodology histogram and density curves of data derived from Experiment 2. A rightshift in the PDF is denoted by an arrow. Bottom panel: ECDFs of data derived from Experiment 2. The D statistic is plotted, showing the maximal distance between the curves.
Discussion
There is a lack of standardization in collecting and analyzing phagocytosis assay data. While several techniques offer solutions to quantification of actively phagocytosing cells, the increase of phagocytic activity, i.e., the relative amount or magnitude of phagocytic increases/decreases, is not investigated. Here we propose using the KS test following transformation of the continuous data into an ECDF. The KS test is supported by the theory of optimal mass transport, which ultimately describes two distributions, and seeks to calculate the distances between them. Many principles of optimal mass transport theory, such as those that utilize the Wasserstein distance metric or Earth Mover’s distance, have been used in various areas of research from engineering, to machine learning, to cell biology (Kolouri et al., 2018; Schiebinger et al., 2019). Such analyses of CDFs thus deeply describe the higher order PDF and realworld data from which it is derived. This is a sensitive and accurate technical pipeline to collect and analyze data from monocytic cell phagocytosis assays. The finer detail of data derived reflects the magnitude of changes these innate immune cells undergo. While we show that the signal captured by our method correlates directly with number of beads per cell, investigators should be aware of experimental caveats, such as background from nonengulfed beads, cellcell engulfment, as well as incomplete capture of the axial plane.
Acknowledgments
We would like to thank members of the Tsirka lab and Dr Rome Sandhu for helpful discussions and edits. This work was partially supported by K12GM102778 (KKT), R25 CA214272 (MT) and SBU seed funding (SET).
Competing interests
The authors declare no conflicts of interest or competing interests.
References

Caponegro, M. D., Torres, L. F., Rastegar, C., Rath, N., Anderson, M. E., Robinson, J. K. and Tsirka, S. E. (2019). Pifithrinmu modulates microglial activation and promotes histological recovery following spinal cord injury. CNS Neurosci Ther 25(2): 200214.

Kolouri, S., Pope, P. E., Martin, C. E. and Rohde, G. K. (2018) Slicedwasserstein autoencoder: an embarrassingly simple generative model. arXiv eprints.

Lian, H., Roy, E. and Zheng, H. (2016). Microglial phagocytosis assay. Bioprotocol 6(21): e1988.

Pul, R., Chittappen, K. P. and Stangel, M. (2013). Quantification of microglial phagocytosis by a flow cytometerbased assay. Methods Mol Biol 1041: 121127.

Roy, S., Bag, A. K., Dutta, S., Polavaram, N. S., Islam, R., Schellenburg, S., Banwait, J., Guda, C., Ran, S., Hollingsworth, M. A., Singh, R. K., Talmadge, J. E., Muders, M. H., Batra, S. K. and Datta, K. (2018). Macrophagederived neuropilin2 exhibits novel tumorpromoting functions. Cancer Res 78(19): 56005617.

Scheraga, R. G., Abraham, S., Niese, K. A., Southern, B. D., Grove, L. M., Hite, R. D., McDonald, C., Hamilton, T. A. and Olman, M. A. (2016). TRPV4 mechanosensitive ion channel regulates lipopolysaccharidestimulated macrophage phagocytosis. J Immunol 196(1): 428436.

Schiebinger, G., Shu, J., Tabaka, M., Cleary, B., Subramanian, V., Solomon, A., Gould, J., Liu, S., Lin, S., Berube, P., Lee, L., Chen, J., Brumbaugh, J., Rigollet, P., Hochedlinger, K., Jaenisch, R., Regev, A. and Lander, E. S. (2019). Optimaltransport analysis of singlecell gene expression identifies developmental trajectories in reprogramming. Cell 176(4): 928943 e922.

Schindelin, J., ArgandaCarreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., Preibisch, S., Rueden, C., Saalfeld, S., Schmid, B., Tinevez, J. Y., White, D. J., Hartenstein, V., Eliceiri, K., Tomancak, P. and Cardona, A. (2012). Fiji: an opensource platform for biologicalimage analysis. Nat Methods 9(7): 676682.

Thompson, K.K., and Tsirka, S.E. (2020). Guanabenz modulates microglia and macrophages during demyelination. Scientific reports 10, 19333.

Wu, T. T., Chen, T. L. and Chen, R. M. (2009). Lipopolysaccharide triggers macrophage activation of inflammatory cytokine expression, chemotaxis, phagocytosis, and oxidative ability via a tolllike receptor 4dependent pathway: validated by RNA interference. Toxicol Lett 191(23): 195202.

Zhang, X., Goncalves, R. and Mosser, D. M. (2008). The isolation and characterization of murine macrophages. Curr Protoc Immunol Chapter 14: Unit 14 11.