1. Introduction
TP53 wild-type human breast tumors rarely undergo pathological complete response following anticancer therapy [
1,
2], resulting in poor survival [
3]. Their unresponsiveness to treatment has been attributed to the proneness of tumor cells with normal p53 to avoid cell death by undergoing senescence [
4], a dormant but metabolically active state. Senescent cells can persist for a long time and secrete inflammatory factors that can drive relapse by promoting survival, proliferation, and stemness [
4,
5]. Thus, understanding how senescent cancer cells can survive for long periods after treatment is paramount for developing more effective anticancer therapies. Recently, Tonnessen-Murray and colleagues have shown that
TP53 wild-type breast cancer senescent cells engulf and break down vital, neighboring ones. This process confers a survival advantage to engulfing cells, as it provides nutrients and building blocks necessary for survival, explaining in part how they can persist and maintain high metabolism rates in such unfavorable conditions [
6]. Notably, cell-in-cell structures derived from cell engulfment can be observed in patient tissues of several tumor types after chemotherapy and are associated with negative prognosis [
7,
8,
9]. As this phenomenon is not yet fully understood, the efforts made towards the development of effective strategies to eliminate cancer engulfing cells have yielded unsatisfying results so far. A peculiarity of this senescence-associated engulfing process is the formation of a big vacuole within engulfing cells. This structure is necessary for the degradation of the internalized cells [
9], but it persists over the span of several hours after their clearance, and it appears empty both in brightfield mode and following staining with several fluorescent markers. While small, transient vacuolar organelles such as lysosomes are crucial for cell survival, extreme vacuolization and large vacuoles are normally associated with cell death in mammalians [
10,
11], thus it is uncanny how cancer senescent cells can persevere in hostile conditions with a seemingly empty vacuole occupying most of their cytoplasm. Yet, despite the great deal of interest in cell-in-cell structures in the recent cancer literature, very little attention has been placed on vacuoles. As they constitute a major feature of tumor senescence-related engulfing, their description, including information about their chemical and morphological composition, is essential for the development of effective methods for the prevention and/or elimination of these aggressive senescent cells. In turn this would foster more efficient, targeted anticancer strategies for those (breast) tumors that are resistant to the conventional treatment.
Cell-in-cell and vacuolated structures are three-dimensional (3D), architecturally complex systems. The 3D reconstruction of a biological specimen and the visualization of the molecular species that compose it are often obtained with fluorescent confocal microscopy [
12], but in order to attain satisfying, reliable results with this technique, one needs to know all the potential, labelable targets and the specimen should efficiently absorb and retain the set of selected fluorescent dyes [
13]. As the vacuolar content after clearance of the engulfed cells is still undefined and the vacuolar membrane may represent a further obstacle to dye penetration, in addition to the surrounding engulfing cell one, neither of these conditions is satisfied. Moreover, the use of dyes causes the perturbation of the biological model, which often leads to artifacts and signals that do not proportionally and/or uniquely relate to the targeted feature [
13]. To overcome these problems, we employed optical diffraction tomography (ODT) and confocal Raman microscopy for label-free investigation of morphological and chemical information, respectively.
ODT is an interferometric technique that measures the refractive index (RI) of transparent samples, such as cells and subcellular structures [
14]. The RI value of an item is linearly proportional to its mass density [
15], making ODT particularly suited to probe the hollowness of engulfing-derived vacuoles, and, as the measure is label-free, it avoids typical dye-related issues [
12,
14]. Moreover, ODT measures the 3D RI distribution of a specimen, delivering 3D reconstructions of complex biological structures featuring sub-micrometer spatial resolution [
12]. However, ODT lacks the molecular specificity necessary to interpret the biochemical content of vacuoles [
12]. The solution to this problem comes from Raman spectroscopy. Raman spectroscopy uses the vibrational motions of chemical bonds as a contrast mechanism, revealing the biomolecular composition of the sample in a label-free manner, with sub-cellular spatial resolution (
i.e., micrometer scale resolution). This provides an endogenous signature that can be used as a fingerprint for the unique characterization of a specimen chemical content in terms of lipids, proteins, nucleic acids, carbohydrate concentration, and many others [
16]. Our custom-build high-speed confocal Raman microscope has been successfully used for various applications including malaria diagnosis [
17], monitoring the cellular uptake of single walled carbon nanotubes [
18], monitoring cell-drug interaction on multiple myeloma cells [
19], and monitoring UV-induced cell damage [
20]. By applying a combination of ODT and confocal Raman microscopy to MCF7
TP53 wild-type, doxorubicin (Doxo)-induced senescent engulfing cells, we were able to visualize the 3D structure of live engulfing cells, including vacuoles, and investigate vacuolar content in a completely label-free fashion, phenotyping these treacherous cells in unperturbed samples.
4. Discussion
In this study, we employed ODT and confocal Raman microscopy to gain insight into the nature and content of vacuoles that generate inside therapy-resistant cancer senescent cells after engulfing and degradation of other vital, neighboring ones. This process has been observed in several tumors
in vivo [
8] and is correlated with poor prognosis as it provides nutrients and building blocks to the engulfing cells, allowing for the persistence of deleterious cancer cells even in hostile environments [
6,
7]. By measuring the RI of vacuole-holding cells we found that vacuole density was lower than that of the cytoplasm and comparable to the PBS or culture medium. Conversely, the chemical map of vacuoles obtained through Raman microscopy differed appreciably from that of the aqueous medium (PBS or RPMI). The Raman profile of vacuoles was qualitatively identical to that of the corresponding surrounding cell, but peak intensities in vacuoles were markedly lower. As peak intensity scales directly with the concentration of the corresponding chemical species, this result indicates that biomolecules were in small concentration within the vacuoles, consistent with the low mass density in vacuoles measured through ODT. Taken together, our results indicate that engulfing-originated vacuoles were not comparable to the medium despite their appearance, instead they contained diluted biomolecules characteristics of cells, most likely derived from previously broken-down engulfed cells.
The extremely low RI of vacuoles may suggest that the breakdown of the engulfed cells within the vacuoles is not executed through enzymatic activity, as enzymatic degradation of different biomolecules requires the presence of several macromolecules/protein complexes, hence it is reasonable to assume that vacuolar density and the intensity of Raman peaks of proteins would be higher. In the last decade, it was proposed that internalized cells as found in cell-in-cell structures that originated both after senescence induction [
6] and from another, similar type of engulfment called entosis are killed by a mechanism dependent on acidified lysosomes [
35,
36]. In a recent study, Su and colleagues [
9] showed how a remarkable decline in vacuolar pH determines the death of internalized cells following entosis in several breast cancer cell lines, and described vacuoles as ‘’huge lysosomes’’. It is reasonable to hypothesize that a similar mechanism applies to our model. Notably, senescent cells have an exceptionally expanded lysosomal compartment [
37] and, while the inside of vacuoles was negative to several common fluorescent markers, we found that the vacuolar membrane was positive to the lysosome-targeting LysoTracker probe (
Supplementary Materials, Figure S2). This indicates that vacuoles really are big lysosome-like structures inside which internalized cells get broken down, and implies that the vacuole membranes effectively segregate the acid environment from the cytoplasm of the engulfing cell, to prevent it from getting damaged. This, in turn, would explain in part why vacuoles appear empty even following different staining, as the vacuolar membrane might block the dyes out.
By statistically comparing the spectra of vacuoles with those of the cells surrounding them, we found that the two groups differed solely on peak intensities, and not on the presence or absence of specific Raman bands. The Raman intensity profile of vacuoles was half that of the cells, but with the notable exception of the Raman region between 800 and 900 cm
-1 and the peak at 965, which intensities were equally great or greater in vacuoles than in cells. Raman peaks at 837, 845, 855 and 965 cm
-1 can be attributed to tyrosine [
33], while the peak at 874 cm
-1 is associated with tryptophan [
34]. Interestingly, tyrosine and tryptophan are two of the three aromatic amino acids (AAA; the third being phenylalanine), and elevated levels of AAA have been found in breast cancers [
33,
38]. AAA are critical intermediates that connect nucleotide, glucose, and lipid metabolism, and represent key precursor to many biological compounds involved in the regulation of many metabolic pathways and necessary for cellular functioning [
33,
38,
39]. As the content of vacuoles is supposedly transferred to the surrounding senescent engulfing cells, the relatively predominant presence of tyrosine and tryptophan within the vacuoles could reflect an elevated demand for these amino acids from engulfing cells to support their exceptionally elevated metabolism [
40,
41], leading to an inflated degradation of AAA-containing peptides to meet that demand. On the other hand, AAA are involved in the production of reactive oxygen species, with deleterious effects on cells [
39]. Thus, it is possible that tyrosine and tryptophan are kept within the confined space of vacuoles as long as possible, instead of being transferred quickly to the cell like the other biomolecules, to prevent them for damaging the cell itself, which would result in their relative accumulation in vacuoles. Finally, tyrosine and tryptophan are the biggest amino acids (181 and 204 Da, respectively), with great steric hindrance. The simplest explanation for their relative buildup inside the vacuoles might be that their size obstacles their transit through the vacuole membrane and into the engulfing cell cytoplasm.
In conclusion, in this work we showed a feasibility of using ODT and Raman for monitoring cell engulfment process from a limited number of cells. We demonstrated how this approach can be greatly beneficial for describing and monitoring this and other complex biological processes in a non-invasive, label-free fashion, crucial for more effective biomedical research.
Further studies with extended cell dataset will follow to confirm the mechanism underlying the degradation of internalized cells and to investigate the process through which the biomolecules derived from the breakdown of engulfed cells are transferred to the engulfing ones to be exploited as nutrients/building blocks. Although our existing custom-built Raman microscope allows high-speed Raman imaging of live cells [
21,
42], it is still based on the point-scanning measurement. Obtaining data from a large number of cells or tissue slices is challenging. We are currently building a new Raman imaging system to overcome this limitation exploiting the coherent Raman mechanism [
43]. This new system will allow collecting Raman data from large number of live cells in a short period of time so that we can extend the study with statistical analysis. Future studies will also include time-lapse monitoring of the overall process, instead of discrete snapshots as shown here.
Author Contributions
Conceptualization, S.G., J.W.K., D.P. and I.B.; methodology, S.G., J.W.K., A.B. and K.J.K.K.; software, S.G., A.B. and K.J.K.K.; formal analysis, S.G. and A.M.; writing—original draft preparation, S.G.; writing—review and editing, S.G., J.W.K., A.B., D.P, A.M. and I.B.; visualization, S.G. and I.B.; supervision, J.W.K., P.T.C.S. and I.B.; funding acquisition, I.B., P.T.C.S. and D.P. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Doxorubicin-treated cancer senescent cells engulf and degrade neighboring ones. (A) Cartoon illustrating the engulfing process: a senescent cell engulfing a neighboring one (1, 2), originating a cell-in-cell structure (3). The engulfed cell is broken down in a vacuole within the engulfing one (4). After the degradation of the engulfing cell, the vacuole appears empty (5) and shrinks over time (6) until disappearing (7). Created with BioRender.com. (B) Time-lapse image sequence of a doxorubicin-treated MCF7 cell (1) engulfing a neighboring one (2), leading to the creation of a cell-in-cell structure. Cells were treated with doxorubicin 48 hours before starting the time-lapse experiment. Scale bar: 25 μm. (C) Time-lapse image sequence of a vacuole shrinking over time inside an MCF7 doxorubicin-treated cell. Cells were treated with doxorubicin 90 hours before starting the time lapse experiment. Scale bar: 25 μm. (D) Control and doxorubicin-treated MCF7 cells stained for β-galactosidase activity (green) and nuclei (blue). In the last image, a cell-in-cell structure derived by the engulfing process. Scale bar: 25μm.
Figure 1.
Doxorubicin-treated cancer senescent cells engulf and degrade neighboring ones. (A) Cartoon illustrating the engulfing process: a senescent cell engulfing a neighboring one (1, 2), originating a cell-in-cell structure (3). The engulfed cell is broken down in a vacuole within the engulfing one (4). After the degradation of the engulfing cell, the vacuole appears empty (5) and shrinks over time (6) until disappearing (7). Created with BioRender.com. (B) Time-lapse image sequence of a doxorubicin-treated MCF7 cell (1) engulfing a neighboring one (2), leading to the creation of a cell-in-cell structure. Cells were treated with doxorubicin 48 hours before starting the time-lapse experiment. Scale bar: 25 μm. (C) Time-lapse image sequence of a vacuole shrinking over time inside an MCF7 doxorubicin-treated cell. Cells were treated with doxorubicin 90 hours before starting the time lapse experiment. Scale bar: 25 μm. (D) Control and doxorubicin-treated MCF7 cells stained for β-galactosidase activity (green) and nuclei (blue). In the last image, a cell-in-cell structure derived by the engulfing process. Scale bar: 25μm.
Figure 2.
Representative images of PFA fixed MCF7 cells by label-free cell tomography showing the main phases of the senescence-associated engulfment process. Each senescent engulfing cell undergoes these three stages in this order; two cells for each stage are shown. The yellow lines indicate the RI-based segmentation. (A) The engulfed cell (overlined with a discontinued red line) was clearly visible and took up all the space inside the cytoplasm of the engulfing cell, pressing its nucleus against the cell membrane. This caused the nucleus of the engulfing cell to acquire a half-moon shape, indicated with a white discontinued line. All the cellular components presented a higher RI than the culture medium. (B) In this phase, the engulfed cell was being degraded within a vacuole in the engulfing one; it did not take up all the space in the engulfing cell cytoplasm anymore. Some regions of the vacuole presented the same RI as the culture medium. (C) The engulfed cell was not visible anymore, meaning it had been completely degraded. In the place of the engulfed cell was a hollow space, with a RI comparable to that of the aqueous medium. (D) Representative brightfield, ODT and RI-based 3D reconstruction of a cell in the empty stage of the engulfing process. On the brightfield image, the white discontinued line outlines the nucleus of the engulfing cell, the black discontinued one indicates the vacuole, and the continue black line delineates the engulfing cell.
Figure 2.
Representative images of PFA fixed MCF7 cells by label-free cell tomography showing the main phases of the senescence-associated engulfment process. Each senescent engulfing cell undergoes these three stages in this order; two cells for each stage are shown. The yellow lines indicate the RI-based segmentation. (A) The engulfed cell (overlined with a discontinued red line) was clearly visible and took up all the space inside the cytoplasm of the engulfing cell, pressing its nucleus against the cell membrane. This caused the nucleus of the engulfing cell to acquire a half-moon shape, indicated with a white discontinued line. All the cellular components presented a higher RI than the culture medium. (B) In this phase, the engulfed cell was being degraded within a vacuole in the engulfing one; it did not take up all the space in the engulfing cell cytoplasm anymore. Some regions of the vacuole presented the same RI as the culture medium. (C) The engulfed cell was not visible anymore, meaning it had been completely degraded. In the place of the engulfed cell was a hollow space, with a RI comparable to that of the aqueous medium. (D) Representative brightfield, ODT and RI-based 3D reconstruction of a cell in the empty stage of the engulfing process. On the brightfield image, the white discontinued line outlines the nucleus of the engulfing cell, the black discontinued one indicates the vacuole, and the continue black line delineates the engulfing cell.
Figure 4.
Spontaneous Raman can discriminate between cells, vacuoles, and culture medium. (A) Representative Raman average spectra of the cell around the vacuole (red), the vacuole (green) and the (culture) medium of a vital MCF7 engulfing cell. The main Raman features are reported on the spectra. Next to spectra, the brightfield picture of the cell. (B) PCA space showing Raman spectra of medium (black), vacuoles (green), and cells surrounding vacuoles (red). (C) Scatter matrix obtained by plotting the most significant PC scores (MLR coefficients featuring p<0.01) for the separation of the three groups (cell around vacuole, vacuole, and background): PC1, PC6 and PC17. (D) Loading spectra of PC1 and PC6. The main Raman features are reported on the spectra (DNA P bond = DNA phosphodiester bond; Phe = phenylalanin). (E) Representative images of an engulfing cell imaged in brightfield mode (a) and in spontaneous Raman mode (b, c) by selecting the full spectral range (600-1800 cm-1, b) or selecting the peak 855 (c).
Figure 4.
Spontaneous Raman can discriminate between cells, vacuoles, and culture medium. (A) Representative Raman average spectra of the cell around the vacuole (red), the vacuole (green) and the (culture) medium of a vital MCF7 engulfing cell. The main Raman features are reported on the spectra. Next to spectra, the brightfield picture of the cell. (B) PCA space showing Raman spectra of medium (black), vacuoles (green), and cells surrounding vacuoles (red). (C) Scatter matrix obtained by plotting the most significant PC scores (MLR coefficients featuring p<0.01) for the separation of the three groups (cell around vacuole, vacuole, and background): PC1, PC6 and PC17. (D) Loading spectra of PC1 and PC6. The main Raman features are reported on the spectra (DNA P bond = DNA phosphodiester bond; Phe = phenylalanin). (E) Representative images of an engulfing cell imaged in brightfield mode (a) and in spontaneous Raman mode (b, c) by selecting the full spectral range (600-1800 cm-1, b) or selecting the peak 855 (c).
Figure 5.
Spontaneous Raman to discriminate between vacuoles in different stages of the engulfing process. (A) Illustration (first row) and brightfield vacuole pictures (second row) showing their classification into the three groups (i) ’’engulfed cell visible’’, (ii) ‘’big’’, and (iii) ‘’small’’. Created with BioRender.com. (B) PCA of Raman spectra of engulfed cell visible (pink), big (blue), and small (green) vacuoles. (C) Average spectra of engulfed cell visible (pink), big (blue), and small (green) vacuoles. Shaded regions display values within one standard deviation of the mean.
Figure 5.
Spontaneous Raman to discriminate between vacuoles in different stages of the engulfing process. (A) Illustration (first row) and brightfield vacuole pictures (second row) showing their classification into the three groups (i) ’’engulfed cell visible’’, (ii) ‘’big’’, and (iii) ‘’small’’. Created with BioRender.com. (B) PCA of Raman spectra of engulfed cell visible (pink), big (blue), and small (green) vacuoles. (C) Average spectra of engulfed cell visible (pink), big (blue), and small (green) vacuoles. Shaded regions display values within one standard deviation of the mean.