Next Article in Journal
A Systematic Review of the Role of BIM in Building Sustainability Assessment Methods
Next Article in Special Issue
Subpopulations of Organoid-Forming Cells Have Different Motility
Previous Article in Journal
Concepts and Terminologies in 4D Printing
Previous Article in Special Issue
Quantitative Phase Dynamics of Cancer Cell Populations Affected by Blue Light
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Label-Free Classification of Apoptosis, Ferroptosis and Necroptosis Using Digital Holographic Cytometry

by
Kendra L. Barker
1,
Kenneth M. Boucher
1,2 and
Robert L. Judson-Torres
1,3,*
1
Huntsman Cancer Institute, University of Utah, Salt Lake City, UT 84112, USA
2
Department of Medicine, Division of Epidemiology, University of Utah School of Medicine, Salt Lake City, UT 84112, USA
3
Department of Dermatology, University of Utah School of Medicine, Salt Lake City, UT 84112, USA
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(13), 4439; https://doi.org/10.3390/app10134439
Submission received: 26 May 2020 / Revised: 23 June 2020 / Accepted: 25 June 2020 / Published: 27 June 2020
(This article belongs to the Special Issue Applications of Digital Holography in Biomedical Engineering)

Abstract

:

Featured Application

Digital holographic cytometry for classification of programmed cell death.

Abstract

Apoptosis, ferroptosis and necroptosis are three distinct forms of programmed cell death. Each of these pathways can be exploited to terminate cancer cells. One promising therapeutic strategy is to activate alternative programmed cell death pathways subsequent to cancer cells evolving mechanisms to evade apoptosis. However, the interplay between distinct programmed cell death pathways and cancer progression is complex and can paradoxically promote the disease. There is a need for high-throughput assays for real-time classification of programmed cell death, both to further investigate these important biologic processes and to assess the case-by-case efficacy of targeting each pathway in patient-derived tumor cells. Here, we sought to develop a label-free, live-imaging-based assay for classifying forms of programmed cell death with single cell resolution. We used digital holographic cytometry (DHC) to monitor human melanoma cells undergoing apoptosis, ferroptosis, and necroptosis. We developed and validated models that used DHC-derived features to classify each form of cell death with 91–93% accuracy in the test sets. We conclude that high-accuracy, high-throughput, label-free classification of apoptosis, ferroptosis and necroptosis can be achieved with DHC.

1. Introduction

The ability to evade programmed cell death is one of the hallmarks of cancer, leading to therapeutic resistance or failure [1]. Although apoptosis is perhaps the most famous, other forms of programed cell death have been well-characterized [2]. Previous studies have demonstrated that cancer cells that are resistant to one type of cell death may be more vulnerable to another. Ferroptosis and necroptosis, for example, can be exploited to induce death in cancer cells that have successfully evaded apoptosis [3,4,5,6,7,8]. However, the interplay between each of these programed cell death pathways and their involvement in cancer progression is complex. For example, pan-caspase inhibitors have been shown to inhibit apoptosis while simultaneously inducing necroptosis [9]. Induction of necroptosis can be antitumorigenic in some settings or promote tumor progression in others [10,11,12,13,14]. In another example, vulnerability of melanoma cells to induction of ferroptosis is dependent on the phenotypic subtype [15]. Such observations highlight the need for further investigation and suggest that the therapeutic potential of targeting distinct programmed cell death pathways might require case-by-case assessment. Methods that efficiently and simultaneously quantify apoptosis, ferroptosis, and necroptosis would be valuable for further exploring the clinical utility of agents that regulate cell death pathways.
Apoptosis, ferroptosis and necroptosis are all active forms of cell death, meaning they comprise defined molecular pathways that drive organized termination. The key regulatory molecules that activate, inhibit, or orchestrate each process permit molecular classification of cell death pathways using a variety of methods [2]. For example, changes in expression, post-translational processing, or subcellular localization of key proteins can be monitored using immunological methods such as fluorescence microscopy, immunoblotting, flow analysis, or enzyme-linked immunoassays (ELISA). Other useful methodologies detect outer membrane permeabilization or mitochondria function and health—both of which change during the process of apoptosis [16]. While each of these approaches has distinct advantages and limitations, none permit both high-throughput classification of all three forms of cell death in real time and with single cell resolution.
Digital Holographic Cytometry (DHC) is a label-free technique that allows for live cell imaging and real time morphologic assessment [17]. DHC uses a laser that is split into two beams—a reference beam and an object beam. The reference beam is undisturbed, while the object beam is directed through the specimen. When the two beams are merged, an image can be generated through interpretation of the resulting interference pattern, with pixel intensity representing optical thickness—a function of both physical thickness and optical density. Cells are identified using standard image segmentation techniques, and individual cell features (optical volume, optical thickness, perimeter length, etc.) are derived. Quantitative phase images can be captured over a period of time, permitting tracking of DHC-derived features as cells undergo change, for example, due to exposure to a therapeutic agent. DHC has previously been used to distinguish different cell states based upon morphology, including identification of pre-apoptotic and dying cells [18,19,20,21,22].
In addition to molecular changes, a distinct set of morphological changes accompanies each form of programmed cell death. Apoptosis is associated with an increase in cell thickness, plasma membrane blebbing, decrease in cell volume, and chromatin condensation [23,24]. Ferroptosis is associated with iron-dependent membrane lipid oxidation and damage to the plasma membrane [24,25]. Necroptosis is associated with cell swelling, an increase in cell thickness, decrease in cell area, and eventual rupture of the plasma membrane [4,23,24]. We therefore reasoned that analysis of cellular morphologic change with DHC might permit label-free, live, and high-throughput classification of apoptosis, ferroptosis and necroptosis with single cell resolution. The goal of this study was to train and test a classifier to separate three types of cell death using DHC-derived cell features.

2. Materials and Methods

2.1. Cell Culture

The BRAFV600E-mutated, human-derived, metastatic melanoma cell line 501mel (CVCL_4633) was obtained from the University of California, San Francisco (a generous gift from Dr. Boris Bastian). It was obtained under ethical guidelines and informed consent was signed [26]. As all information regarding the donating patient has been de-identified, the use of the line for this study is not considered Human Subjects Research by the National Institutes of Health. Cells were cultured in growth media consisting of 10% Fetal Bovine Serum (VWR, Radnor, PA, USA), 1% L-Glutamine (Gibco Thermo Fisher Scientific, Gaithersburg, MD, United States), 1% Non-essential Amino Acid supplement (Gibco Thermo Fisher Scientific, Gaithersburg, MD, United States), and 1% Penicillin-Streptomycin (Gibco Thermo Fisher Scientific, Gaithersburg, MD, United States) in RPMI-1640 medium (Gibco Thermo Fisher Scientific, Gaithersburg, MD, United States). The cells were passaged twice a week at a ratio of 1 to 10. The cells were kept in a standard 37 °C humidified mammalian cell incubator with 5% CO2 (Thermo Fisher Scientific, Waltham, MA, United States).
Preliminary dose–response experiments were conducted to determine the IC50 of each compound. Briefly, the cells were seeded in 6-well plates (Sarstedt, Nümbrecht, Germany) at a density of 60k cells per well in 3 mL of growth media. The cells were incubated for approximately 24 h after being seeded. Growth media was then aspirated and 3 mL fresh growth media supplemented with indicated concentrations of dimethyl sulfoxide (DMSO, MilliporeSigma, St. Louis, MO, United States), staurosporine (Tocris Bioscience, Minneapolis, MN, United States, Cat. No. 1285, 0.001–10.0 μm in DMSO), erastin (Selleckchem, Houston, TX, United States, Cat. No. S7242, 0.1–30 μm in DMSO), or shikonin (Selleckchem, Houston, TX, United States, Cat. No. S8279, 0.1–10.0 μm in DMSO) was added. The standard plate lid was replaced with an optically clear Hololid (Phase Holographic Imaging, Lund, Sweden). After addition of the Hololid, plates were incubated for 45 min to eliminate condensation. The dose–response application in AppSuite (v2.1.2.0, Phase Holographic Imaging, Lund, Sweden) was then used to determine the IC50.
For classification experiments, cells were cultured and treated as above, except that each plate contained a single well of each compound at the calculated IC50. The 6-well plate was then loaded onto the DHC platform and imaged via AppSuite. The parameters were set to 48 h, with images taken every 1 h from 20 random fields. The training and test sets were generated using independent passages, plates of cells, and freshly generated working solutions of each compound, and were conducted at different times. In two test set experiments, the erastin-treated well was lost due to technical reasons resulting in a sample size of two instead of four for this condition.

2.2. Digital Holographic Imaging and Cytometry

The M4 HoloMonitor (Phase Holographic Imaging, Lund, Sweden) uses a form of quantitative phase imaging called digital holographic microscopy to create a digitally reconstructed image of cells [17]. Images were segmented and analyzed using the Hstudio software (v2.7.5, Phase Holographic Imaging, Lund, Sweden). HStudio was used to derive 32 quantitative cell features for individual cells. For each condition, image stacks were monitored for cell death and time points selected at which cells were undergoing morphological changes prior to cell death. In control conditions, M-phase cells were excluded from the analyses.

2.3. Data Analysis and Statistics

Pearson’s correlation coefficient was used to assess the association between all features. Blocks of features with significant correlation (r > 0.96) were refined to single features. One-way ANOVAs with Tukey’s multiple comparisons test were used to determine significant differences in unique feature means between each of the treatments. The Tukey multiple comparison procedure was chosen as the goal was to adjust for all pairwise comparisons. Features that were significantly associated (adjusted p value < 0.0001) with at least two treatments were selected for further analysis. Decision trees were constructed by first identifying individual features that best separated each pair of conditions. Thresholds were identified by optimizing for accuracy of classification. Finally, the ordering of nodes was conducted manually through a rational consideration of the feature distributions. Sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and accuracy were calculated using standard methods [27].

3. Results

3.1. Identification of Informative DHC-Derived Features

To investigate whether DHC-derived features could be used to distinguish cells undergoing apoptosis, ferroptosis and necroptosis, we generated four replicate datasets (Figure 1). Human melanoma cells (501mel) were treated with compounds previously established as inducing apoptosis [28,29,30,31] (staurosporine, 0.1 μm), ferroptosis [8,25,32,33,34,35] (erastin, 2 μm), or necroptosis [36,37,38,39,40,41] (shikonin, 1.5 μm). Cells treated with only DMSO were included as controls (not treated, NT). Cultures were live-imaged and DHC-derived features were generated for each cell. Two of the experiments contained all four conditions, whereas two were performed without erastin. One of the four-condition experiments was designated the Training Set and the rest as Test Sets. To first determine whether data generated from independent experiments were comparable, we evaluated the optical thickness (average pixel intensity in object) of all cells within the NT conditions (Figure 2). Optical thickness is a fundamental measurement of DHC, from which other morphological features are derived, so we reasoned that comparable datasets should yield comparable distributions of this feature in NT conditions, where biological variability should be minimal. Indeed, optical thickness distributions were comparable across datasets, indicative of a stable platform and reproducible features. We next compared images from each condition (Figure 3). We observed notable and distinct morphological changes upon treatment with each of the three compounds. The distribution of average optical thickness demonstrated somewhat greater variability in the treated conditions as compared to the NT conditions, indicative of the expected increase in biological variability (Figure 2). We therefore investigated whether changes in morphological features were sufficiently consistent across experiments to classify each form of programed cell death despite biological variability.
Many of the morphological features generated by HStudio are highly correlated. In order to avoid over-training, we sought to limit the features considered in classification models. We calculated the correlation coefficient for each combination of features in the Training Set (Figure 4). For features with perfect or near-perfect (r > 0.96) correlations, only one was selected for further analysis. The choice of feature was either arbitrary or, where appropriate, driven by ease of biological interpretation (for example, optical thickness was selected over phaseshift). We compared the distribution of all selected features across the four conditions in the Training Set. Nine features demonstrated a significant difference in distribution (adj. P < 0.0001) between at least two conditions (Figure 5 and Supplemental Tables 1–4). These nine features were selected as the basis for two approaches to classifier generation (Figure 1).

3.2. Optical Volume and Optical Thickness Distinguish Form of Programmed Cell Death

In our first approach, we sought to train a classifier that, when presented with a dying cell, could predict the form of programmed cell death. For this reason, healthy cells were excluded from this analysis (Figure 1). Cells undergoing necroptosis presented a significant increase in thickness, as previously observed using other systems, whereas cells undergoing ferroptosis presented substantially decreased volume (Figure 6). Consequently, the form of programmed cell death could be classified with high accuracy using just these two features (Figure 7). We observed an overall accuracy of 89% of the Training Set. When the trained model was then presented with unmodified data (i.e., not normalized for batch effects) from each of the three Test Sets, accuracies remained high, ranging from 87 to 97% (Supplemental Table 5). Accuracy were high for the datasets lacking the erastin condition, indicating the performance of the model was not dependent on receiving equally-sized groups as input. When all Test Set data were combined, the model performed with an overall accuracy of 91% (Figure 7).

3.3. A four-Node Classifier Distinguishes Form of Programmed Cell Death from Healthy Cells

Although encouraging, our first approach is only useful when classifying a cell pre-identified as undergoing cell death. A classifier that could identify different forms of programmed cell death from the diverse morphologies adopted by healthy cells could be a more versatile tool. This is a more challenging demand due, in large part, to the morphological heterogeneity of both healthy cells and staurosporine-treated cells (Figure 3 and Figure 5). Although we noted several features with distributions that were significantly different across each of the four conditions (for example, Figure 5, optical thickness average (ave)), the range of values in each condition were sufficiently overlapping to prohibit accurate classification of individual cells based upon any one or two features. However, when considering four features, we succeeded in constructing a decision tree that generated overall accuracies of 92% and 93% when presented with the Training and Test Sets, respectively (Figure 8). Necroptosis and ferroptosis remained identifiable with 97% and 95% accuracies based upon optical thickness and optical volume alone. Of the remaining unclassified cells, apoptotic cells were identified with an accuracy of 88% by a combination of large fluctuation in phaseshift (potentially due to the established process of membrane blebbing) and a decrease in length. Similar to the two-node classifier, presentation of each Test Set individually resulted in a tight range of high overall accuracies (91–98%), even in the absence of erastin, indicating the performance of the model was not dependent on a single condition or Test Set (Supplemental Table 6).

4. Discussion

Here we present the application of DHC to the classification of human melanoma cells undergoing three active forms of programmed cell death—apoptosis, ferroptosis, and necroptosis. This approach to monitoring and distinguishing forms of programmed cell death has several advantages over other methods. First, digital holographic imaging is label-free and therefore does not require the manipulation of the cells through addition of dyes or introduction of exogenous genetic material. Second, the assay is compatible with live-imaging, permitting real-time monitoring of cell viability and assessment of the kinetics of each form of programmed cell death. Third, the approach has single-cell resolution, and can therefore provide information on the percentages of a population that undergo different forms of programmed cell death. Finally, DHC is a non-invasive and non-endpoint assay, such that once the analysis is complete, the cells remain unadulterated and can be used for further molecular analyses. We expect this assay will be a valuable addition to the toolbox for investigating programed cell death, and that it will be particularly useful for assessing the dynamics and interplay of different forms of cell death in the same culture.
The strengths of this assay make it particularly well suited for real-time assessment of the therapeutic response of cells directly isolated from human tumors. The potential clinical utility of assessing tumor-specific sensitivity to a collection of potential single agent or combination therapies using in vitro approaches has long been appreciated [42]. However, several demanding requirements impede the successful implementation of this approach as a widespread clinically valuable tool for personalized medicine. First, in order to provide efficient turnaround time and to most accurately capture the in situ cell state, tumor-derived cells must be minimally expanded in vitro. As a label-free approach that requires few cells per condition, the assay presented here requires little in vitro manipulation and selection. Second, an ideal assay would capture the phenotypic heterogeneity of both the naïve tumor cells and their response to therapeutic agents. Capturing heterogeneity of response is crucial for assessing the adoption of “persister states”—cells that survive an initial toxic treatment, often leading to therapeutic resistance, but also increasing susceptibility to other therapeutic agents [2]. The assay presented here could, for example, simultaneously compare not only the efficacy of two different agents, but also determine whether phenotypic subsets of tumor cells were predominately affected, and whether different forms of programmed cell death were triggered. Such information would be instructive when selecting candidate combination therapies. For these reasons, DHC platforms have already been implemented in the high-throughput real-time screening of tumor-derived cells [43,44]. The approach presented here offers additional granularity to the information regarding cell death provided by these systems.
We present here the classification of three forms of programmed cell death. However, other major forms of programmed cell death, most notably, pyroptosis, have also been described2. We limited our analyses to the forms of cell death that we could confidently and uniformly induce based upon prior work. It is possible that similar approaches could yield models that also identify pyroptosis. One limitation of this study is the inclusion of only single agents to induce each form of cell death. Each agent used in this study (staurosporin, erastin, shinkonin) is an established and classic inducer of the associated form of cell death [25,29,30,31,32,33,34,35,37,38,39,40,41]. However, it remains possible that the morphological changes measured here are specific to the agent, not the form of cell death. Similarly, an additional limitation to the work presented here is the analysis of a single human melanoma line. As we and others have applied DHC to the classification of different cellular phenotypes [18,45,46,47,48,49,50,51], one emerging concern is that single models cannot be trained that accurately classify diverse cell lines or conditions. Indeed, as the approach is based exclusively on morphology, it should be expected that cells with distinct morphologies cannot be assessed by the same model with any degree of accuracy. For these reasons, we emphasize that, as for any assay, the optimal parameters must be determined for each cellular system, experimental condition, reagent, and biological question. We do not encourage the blind application of the exact decision nodes and thresholds defined here to cell lines other than human melanoma cells. Rather, it is our intent, that by visualizing each step of model development, we have provided a workflow for generating high-accuracy classifiers using DHC-derived features that is accessible to all cell biologists interested in applying this method to their cellular systems.

5. Conclusions

We conclude that three forms of programmed cell death—apoptosis, ferroptosis and necroptosis—can be accurately classified using DHC-derived features, label-free, in real-time, and with single-cell resolution.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-3417/10/13/4439/s1, Table S1: Nine DHC-derived features from Training Set. Table S2: Nine DHC-derived features from Test Set 1. Table S3: Nine DHC-derived features from Test Set 2. Table S4: Nine DHC-derived features from Test Set 3. Table S5: Performance of two-node decision tree with each individual Test Set. Table S6: Performance of four-node decision tree with each individual Test Set.

Author Contributions

Conceptualization, R.L.J.-T.; methodology, K.L.B and R.L.J.-T.; validation, K.L.B. and K.M.B.; formal analysis, K.L.B. and K.M.B.; investigation, K.L.B.; writing—original draft preparation, K.L.B and R.L.J.-T.; visualization, K.L.B., and R.L.J.-T.; supervision, R.L.J.-T.; funding acquisition, R.L.J.-T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a University of Utah 1U4U Innovation fund awarded to Tom Zangle, Alexander Lex and Robert L Judson-Torres and the National Cancer Institute of the National Institutes of Health Award Number P30CA042014 to Huntsman Cancer Institute at the University of Utah.

Acknowledgments

Research reported in this publication utilized the Cancer Biostatistics Shared Resource at Huntsman Cancer Institute at the University of Utah.

Conflicts of Interest

The authors declare no conflict of interest

References

  1. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [Green Version]
  2. Green, D.R. The Coming Decade of Cell Death Research: Five Riddles. Cell 2019, 177, 1094–1107. [Google Scholar] [CrossRef]
  3. Gong, Y.; Fan, Z.; Luo, G.; Yang, C.; Huang, Q.; Fan, K.; Cheng, H.; Jin, K.; Ni, Q.; Yu, X.; et al. The role of necroptosis in cancer biology and therapy. Mol Cancer 2019, 18. [Google Scholar] [CrossRef] [Green Version]
  4. Choi, M.; Price, D.; Ryter, S.W.; Choi, A. Necroptosis: A crucial pathogenic mediator of human disease. JCI Insight 2019, 4. [Google Scholar] [CrossRef] [PubMed]
  5. Christofferson, D.E.; Yuan, J. Necroptosis as an alternative form of programmed cell death. Curr Opin Cell Biol. 2010, 22, 263–268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bebber, C.M.; Müller, F.; Clemente, L.P.; Weber, J.; Karstedt, S.V. Ferroptosis in Cancer Cell Biology. Cancers (Basel) 2020, 12, 164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Mou, Y.; Wang, J.; Wu, J.; He, D.; Zhang, C.; Duan, C.; Li, B. Ferroptosis, a new form of cell death: Opportunities and challenges in cancer. J Hematol Oncol. 2019, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Dixon, S.J.; Lemberg, K.M.; Lamprecht, M.R.; Skouta, R.; Zaitsev, E.M.; Gleason, C.E.; Patel, D.N.; Bauer, A.J.; Gantley, A.M.; Yang, W. Ferroptosis: An iron-dependent form of nonapoptotic cell death. Cell 2012, 149, 1060–1072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Werthmöller, N.; Frey, B.; Wunderlich, R.; Fietkau, R.; Gaipl, U.S. Modulation of radiochemoimmunotherapy-induced B16 melanoma cell death by the pan-caspase inhibitor zVAD-fmk induces anti-tumor immunity in a HMGB1-, nucleotide- and T-cell-dependent manner. Cell Death Dis. 2015, 6. [Google Scholar] [CrossRef] [Green Version]
  10. Seifert, L.; Werba, G.; Tiwari, S.; Ly, N.N.G.; Alothman, S.; Alqunaibit, D.; Avanzi, A.; Barilla, R.M.; Daley, D.; Greco, S.h.; et al. The necrosome promotes pancreatic oncogenesis via CXCL1 and Mincle-induced immune suppression. Nature 2016, 532, 245–249. [Google Scholar] [CrossRef] [Green Version]
  11. Strilic, B.; Yang, L.; Juárez, J.A.; Wachsmuth, L.; Han, K.; Müller, U.; Pasparakis, M.; Offermanns, S. Tumour-cell-induced endothelial cell necroptosis via death receptor 6 promotes metastasis. Nature 2016, 536, 215–218. [Google Scholar] [CrossRef] [PubMed]
  12. McCormick, K.D.; Ghosh, A.; Trivedi, S.; Wang, L.; Coyne, C.B.; Ferris, R.L.; Sarkar, S.N. Innate immune signaling through differential RIPK1 expression promote tumor progression in head and neck squamous cell carcinoma. Carcinogenesis 2016, 37, 522–529. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Geserick, P.; Wang, J.; Schilling, R.; Horn, S.; Harris, P.A.; Bertin, J.; Gough, P.J.; Feoktistova, M.; Leverkus, M. Absence of RIPK3 predicts necroptosis resistance in malignant melanoma. Cell Death Dis. 2015, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ke, H.; Augustine, C.K.; Gandham, V.D.; Jin, J.Y.; Tyler, D.S.; Akiyama, S.K.; Hall, R.; Zhang, J. CYLD inhibits melanoma growth and progression through suppression of the jnk/ap-1 and β1-integrin signaling pathways. J Invest Dermatol. 2013, 133, 221–229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tsoi, J.; Robert, L.; Paraiso, K.; Galvan, C.; Sheu, K.; Lay, J.; Wong, D.; Atefi, M.; Shirazi, R.; Wang, X.; et al. Multi-stage Differentiation Defines Melanoma Subtypes with Differential Vulnerability to Drug-Induced Iron-Dependent Oxidative Stress. Cancer Cell 2018, 33, 890–904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Galluzzi, L.; Zamzami, N.; Rouge, T.; Lemaire, C.; Brenner, C.; Kroemer, G. Methods for the assessment of mitochondrial membrane permeabilization in apoptosis. Apoptosis 2007, 12, 803–813. [Google Scholar] [CrossRef] [Green Version]
  17. Alm, K.; El-Schich, Z.; Miniotis, M.F.; Wingren, A.G.; Janicke, B.; Oredsson, S.M.; Mihaylova, E. Holography—Basic Principles and Contemporary Applications; IntechOpen: London, UK, 2013; Cells and Holograms—Holograms and Digital Holographic Microscopy as a Tool to Study the Morphology of Living Cells. [Google Scholar]
  18. Hejna, M.; Jorapur, A.; Song, J.S.; Judson, R.L. High accuracy label-free classification of single-cell kinetic states from holographic cytometry of human melanoma cells. Sci Rep. 2017, 7. [Google Scholar] [CrossRef]
  19. El-Schich, Z.; Mölder, A.; Tassidis, H.; Härkönen, P.; Miniotis, M.F.; Wingren, A.G. Induction of morphological changes in death-induced cancer cells monitored by holographic microscopy. J. Struct. Biol. 2015, 189, 207–212. [Google Scholar] [CrossRef]
  20. Nicolas, P.; Kühn, J.; Corinne, M.; Pascal, J.; Christian, D.; Magistretti, P.J.; Marquet, P. Early cell death detection with digital holographic microscopy. PLoS ONE 2012, 7. [Google Scholar] [CrossRef]
  21. Huang, D.; Roy, I.J.; Murray, G.F.; JasonReed, Z.T.A.; Teitell, M.A. Identifying fates of cancer cells exposed to mitotic inhibitors by quantitative phase imaging. Analyst. 2020, 145, 97–106. [Google Scholar] [CrossRef]
  22. Khmaladze, A.; Matz, R.L.; Epstein, T.; Jasensky, J.; Holl, M.M.B.; Chen, Z. Cell volume changes during apoptosis monitored in real time using digital holographic microscopy. J. Struct. Biol. 2012, 178, 270–278. [Google Scholar] [CrossRef] [PubMed]
  23. Kroemer, G.; Dallaporta, B.; Resche-Rigon, M. The Mitochondrial Death/Life Regulator in Apoptosis and Necrosis. Annu Rev Physiol. 1998, 60, 619–642. [Google Scholar] [CrossRef] [PubMed]
  24. Kroemer, G.; Galluzzi, L.; Vandenabeele, P.; Abrams, J. Classification of cell death: Recommendations of the Nomenclature Committee on Cell Death 2009. Cell Death Differ. 2009, 16, 3–11. [Google Scholar] [CrossRef] [PubMed]
  25. Xie, Y.; Hou, W.; Song, X.; Yu, Y.; Huang, J.; Sun, X.; Kang, R.; Tang, D. Ferroptosis: Process and function. Cell Death Differ. 2016, 23, 369–379. [Google Scholar] [CrossRef] [Green Version]
  26. Rose, A.E.; Poliseno, L.; Wang, J.; Clark, M.; Pearlman, A.; Wang, G.; Vega-Saenz, d.M.E.; Medicherla, R.; Christos, P.; Shapiro, R.; et al. Integrative genomics identifies molecular alterations that challenge the linear model of melanoma progression. Cancer Res. 2011, 71, 2561–2571. [Google Scholar] [CrossRef] [Green Version]
  27. Parikh, R.; Mathai, A.; Parikh, S.; Sekhar, G.; Thomas, R. Understanding and using sensitivity, specificity and predictive values. Indian J. Ophthalmol. 2008, 56, 45–50. [Google Scholar] [CrossRef]
  28. Zhang, X.D.; Gillespie, S.K.; Hersey, P. Staurosporine induces apoptosis of melanoma by both caspase-dependent and -independent apoptotic pathways. Mol. Cancer Ther. 2004, 3, 187–197. [Google Scholar]
  29. Antonsson, A.; Persson, J.L. Induction of apoptosis by staurosporine involves the inhibition of expression of the major cell cycle proteins at the G2/M checkpoint accompanied by alterations in Erk and Akt kinase activities. Anticancer Res. 2009, 29, 2893–2898. [Google Scholar]
  30. Malsy, M.; Bitzinger, D.; Graf, B.; Bundscherer, A. Staurosporine induces apoptosis in pancreatic carcinoma cells PaTu 8988t and Panc-1 via the intrinsic signaling pathway. Eur. J. Med. Res. 2019, 24, 5. [Google Scholar] [CrossRef]
  31. Chae, H.-J.; Kang, J.; Byun, J.; Han, K.; Kim, D.; Oh, S.; Kim, H.; Chae, S.-W.; Kim, H.R. Molecular mechanism of staurosporine-induced apoptosis in osteoblasts. Pharmacol Res. 2000, 42, 373–381. [Google Scholar] [CrossRef]
  32. Shibata, Y.; Yasui, H.; Higashikawa, K.; Miyamoto, N.; Kuge, Y. Erastin, a ferroptosis-inducing agent, sensitized cancer cells to X-ray irradiation via glutathione starvation in vitro and in vivo. PLoS ONE 2019, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Wang, H.; Liu, C.; Zhao, Y.; Zhang, W.; Xu, K.; Li, D.; Zhou, Y.; Li, H.; Xiao, G.; Lu, B.; et al. Inhibition of LONP1 protects against erastin-induced ferroptosis in Pancreatic ductal adenocarcinoma PANC1 cells. Biochem Biophys Res Commun. 2020, 522, 1063–1068. [Google Scholar] [CrossRef] [PubMed]
  34. Yagoda, N.; Rechenberg, M.v.; Zaganjor, E.; Bauer, A.; Yang, W.; Fridman, D.; Wolpaw, A.; Smukste, I.; Peltier, J.; Boniface, J.; et al. RAS-RAF-MEK-dependent oxidative cell death involving voltage-dependent anion channels. Nature 2007, 447, 864–868. [Google Scholar] [CrossRef] [Green Version]
  35. Dolma, S.; Lessnick, S.L.; Hahn, W.C.; Stockwell, B. Identification of genotype-selective antitumor agents using synthetic lethal chemical screening in engineered human tumor cells. Cancer Cell 2003, 3, 285–296. [Google Scholar] [CrossRef] [Green Version]
  36. Han, W.; Li, L.; Qiu, S.; Lu, Q.; Pan, Q.; Gu, Y.; Luo, J.; Hu, X. Shikonin circumvents cancer drug resistance by induction of a necroptotic death. Mol Cancer Ther. 2007, 6, 1641–1649. [Google Scholar] [CrossRef] [Green Version]
  37. Fu, Z.; Deng, B.; Liao, Y.; Shan, L.; Cai, Z. The anti-tumor effect of shikonin on osteosarcoma by inducing RIP1 and RIP3 dependent necroptosis. BMC Cancer. 2013, 13, 580. [Google Scholar] [CrossRef] [Green Version]
  38. Wada, N.; Kawano, Y.; Fujiwara, S.; Kikukawa, Y.; Okuno, Y.; Tasaki, M.; Ueda, M.; Ando, Y.; Yoshinaga, K.; Ri, M.; et al. Shikonin, dually functions as a proteasome inhibitor and a necroptosis inducer in multiple myeloma cells. Int J. Oncol. 2015, 46, 963–972. [Google Scholar] [CrossRef] [Green Version]
  39. Huang, C.; Luo, Y.; Zhao, J.; Yang, F.; Zhao, H.; Fan, W.; Ge, P. Shikonin Kills Glioma Cells through Necroptosis Mediated by RIP-1. Sobol RW, ed. PLoS ONE 2013, 8, e66326. [Google Scholar] [CrossRef] [Green Version]
  40. Shahsavari, Z.; Karami-Tehrani, F.; Salami, S. Shikonin induced necroptosis via reactive oxygen species in the T-47D breast cancer cell line. Asian Pacific J. Cancer Prev. 2015, 16, 7261–7266. [Google Scholar] [CrossRef] [Green Version]
  41. Lu, B.; Gong, X.; Wang, Z.; Ding, Y.; Wang, C.; Luo, T.; Piao, M.; Meng, F.; Chi, G.; Luo, Y. Shikonin induces glioma cell necroptosis in vitro by ROS overproduction and promoting RIP1/RIP3 necrosome formation. Acta Pharmacol Sin. 2017, 38, 1543–1553. [Google Scholar] [CrossRef] [Green Version]
  42. Shrivastav, S.; Bonar, R.A.; Stone, K.R.; Paulson, D.F. An in Vitro Assay Procedure to Test Chemotherapeutic Drugs on Cells from Human Solid Tumors. Cancer Res. 1981, 40, 4438–4442. [Google Scholar]
  43. Huang, D.; Leslie, K.; Guest, D.; Yeshcheulova, O.; Roy, I.; Piva, M.; Moriceau, G.; Zangle, T.; Lo, R.; Teitell, M.; et al. High-Speed Live-Cell Interferometry: A New Method for Quantifying Tumor Drug Resistance and Heterogeneity. Anal. Chem. 2018, 90, 3299–3306. [Google Scholar] [CrossRef] [PubMed]
  44. Reed, J.; Chun, J.; Zangle, T.; Kalim, S.; Hong, J.; Pefley, S.; Zheng, X.; Gimzewski, J.; Teitell, M. Rapid, massively parallel single-cell drug response measurements via live cell interferometry. Biophys J. 2011, 101, 1025–1031. [Google Scholar] [CrossRef] [Green Version]
  45. Lam, V.; Nguyen, T.; Bui, V.; Chung, B.M.; Chang, L.; Nehmetallah, G.; Raub, C. Quantitative scoring of epithelial and mesenchymal qualities of cancer cells using machine learning and quantitative phase imaging. J. Biomed. Opt. 2020, 25. [Google Scholar] [CrossRef] [PubMed]
  46. Calin, V.; Mihailescu, M.; Scarlat, E.; Baluta, A.; Calin, D.; Kovacs, E.; Savopol, T.; Moisescu, M. Evaluation of the metastatic potential of malignant cells by image processing of digital holographic microscopy data. FEBS Open Bio 2017, 7, 1527–1538. [Google Scholar] [CrossRef] [Green Version]
  47. Mangal, J.; Monga, R.; Mathur, S.; Dinda, A.; Joseph, J.; Ahlawat, S.; Khare, K. Unsupervised organization of cervical cells using bright-field and single-shot digital holographic microscopy. J. Biophotonics 2019, 12. [Google Scholar] [CrossRef]
  48. Sebesta, M.; Egelberg, P.J.; Langberg, A.; Lindskov, J.-H.; Alm, K.; Janicke, B. HoloMonitor M4: holographic imaging cytometer for real-time kinetic label-free live-cell analysis of adherent cells. In Proceedings of the SPIE, San Francisco, CA, USA, 9 March 2016; 2016; p. 971813. [Google Scholar]
  49. Roitshtain, D.; Wolbromsky, L.; Bal, E.; Greenspan, H.; Satterwhite, L.; Shaked, N. Quantitative phase microscopy spatial signatures of cancer cells. Cytom. Part A. 2017, 91, 482–493. [Google Scholar] [CrossRef] [Green Version]
  50. Mölder, A.L.; Persson, J.; El-Schich, Z.; Czanner, S.; Gjörloff-Wingren, A. Supervised classification of etoposide-treated in vitro adherent cells based on noninvasive imaging morphology. J. Med. Imaging. 2017, 4. [Google Scholar] [CrossRef] [Green Version]
  51. Péter, B.; Nador, J.; Juhász, K.; Dobos, A.; Kőrösi, L.; Székács, I.; Patko, D.; Horvath, R. Incubator proof miniaturized Holomonitor to in situ monitor cancer cells exposed to green tea polyphenol and preosteoblast cells adhering on nanostructured titanate surfaces: validity of the measured parameters and their corrections. J. Biomed. Opt. 2015, 20. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic of study design. In four independent experiments, 501mel cells were treated with DMSO (not treated, N.T.), staurosporin (Sta.), erastin (Era.) or shikonin (Shi.). One experiment was designated the Training Set and used to generate two classification models—one to distinguish the form of programmed cell death among dying cells and one to distinguish the form of programmed cell death and healthy cells. The reproducibility of both models was tested with the remaining datasets.
Figure 1. Schematic of study design. In four independent experiments, 501mel cells were treated with DMSO (not treated, N.T.), staurosporin (Sta.), erastin (Era.) or shikonin (Shi.). One experiment was designated the Training Set and used to generate two classification models—one to distinguish the form of programmed cell death among dying cells and one to distinguish the form of programmed cell death and healthy cells. The reproducibility of both models was tested with the remaining datasets.
Applsci 10 04439 g001
Figure 2. Distribution of average (ave) optical thickness per cell from indicated conditions across four experiments.
Figure 2. Distribution of average (ave) optical thickness per cell from indicated conditions across four experiments.
Applsci 10 04439 g002
Figure 3. Digital holographic cytometry (DHC) images of indicated conditions. Expected form of programmed cell death indicated in parentheses.
Figure 3. Digital holographic cytometry (DHC) images of indicated conditions. Expected form of programmed cell death indicated in parentheses.
Applsci 10 04439 g003
Figure 4. Matrix demonstrating correlation coefficients between DHC-derived cellular features.
Figure 4. Matrix demonstrating correlation coefficients between DHC-derived cellular features.
Applsci 10 04439 g004
Figure 5. Distributions of indicated DHC-derived single cell features across four conditions. Data represent Training Set only. Ordinary one-way ANOVA with Tukey’s multiple comparisons test; Adjusted p value > 0.05 (n.s.), < 0.05 (*), < 0.01 (**), < 0.001 (***), or < 0.0001(****).
Figure 5. Distributions of indicated DHC-derived single cell features across four conditions. Data represent Training Set only. Ordinary one-way ANOVA with Tukey’s multiple comparisons test; Adjusted p value > 0.05 (n.s.), < 0.05 (*), < 0.01 (**), < 0.001 (***), or < 0.0001(****).
Applsci 10 04439 g005
Figure 6. Average optical thickness plotted against optical volume for cells in indicated conditions. Dotted lines indicate threshold used for 2-node model presented in Figure 7.
Figure 6. Average optical thickness plotted against optical volume for cells in indicated conditions. Dotted lines indicate threshold used for 2-node model presented in Figure 7.
Applsci 10 04439 g006
Figure 7. Two-node decision tree for classifying form of programmed cell death. Graphs depict all datapoints from Training Set not yet classified at each node. Dotted lines and “x” indicate threshold used to classify. Tables provide raw results and calculated metrics of accuracy for Training Set and combined Test Sets.
Figure 7. Two-node decision tree for classifying form of programmed cell death. Graphs depict all datapoints from Training Set not yet classified at each node. Dotted lines and “x” indicate threshold used to classify. Tables provide raw results and calculated metrics of accuracy for Training Set and combined Test Sets.
Applsci 10 04439 g007
Figure 8. Four-node decision tree for classifying form of programmed cell death when healthy cells are included. Graphs depict all datapoints from Training Set not yet classified at each node. Dotted lines and “x” indicate threshold used to classify. Tables provide raw results and calculated metrics of accuracy for Training Set and combined Test Sets.
Figure 8. Four-node decision tree for classifying form of programmed cell death when healthy cells are included. Graphs depict all datapoints from Training Set not yet classified at each node. Dotted lines and “x” indicate threshold used to classify. Tables provide raw results and calculated metrics of accuracy for Training Set and combined Test Sets.
Applsci 10 04439 g008

Share and Cite

MDPI and ACS Style

Barker, K.L.; Boucher, K.M.; Judson-Torres, R.L. Label-Free Classification of Apoptosis, Ferroptosis and Necroptosis Using Digital Holographic Cytometry. Appl. Sci. 2020, 10, 4439. https://doi.org/10.3390/app10134439

AMA Style

Barker KL, Boucher KM, Judson-Torres RL. Label-Free Classification of Apoptosis, Ferroptosis and Necroptosis Using Digital Holographic Cytometry. Applied Sciences. 2020; 10(13):4439. https://doi.org/10.3390/app10134439

Chicago/Turabian Style

Barker, Kendra L., Kenneth M. Boucher, and Robert L. Judson-Torres. 2020. "Label-Free Classification of Apoptosis, Ferroptosis and Necroptosis Using Digital Holographic Cytometry" Applied Sciences 10, no. 13: 4439. https://doi.org/10.3390/app10134439

APA Style

Barker, K. L., Boucher, K. M., & Judson-Torres, R. L. (2020). Label-Free Classification of Apoptosis, Ferroptosis and Necroptosis Using Digital Holographic Cytometry. Applied Sciences, 10(13), 4439. https://doi.org/10.3390/app10134439

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop