Skip to main content
Advertisement
  • Loading metrics

Phosphoinositide species and filamentous actin formation mediate engulfment by senescent tumor cells

  • Wesley D. Frey,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Supervision, Validation, Visualization, Writing – review & editing

    Affiliation Tulane School of Medicine, Department of Biochemistry and Molecular Biology, New Orleans, Louisiana, United States of America

  • Ashlyn Y. Anderson,

    Roles Formal analysis, Visualization

    Affiliation Tulane School of Medicine, Department of Biochemistry and Molecular Biology, New Orleans, Louisiana, United States of America

  • Hyemin Lee,

    Roles Formal analysis, Methodology, Validation, Visualization

    Affiliation Tulane School of Medicine, Department of Biochemistry and Molecular Biology, New Orleans, Louisiana, United States of America

  • Julie B. Nguyen,

    Roles Data curation, Methodology, Writing – review & editing

    Affiliation Tulane School of Medicine, Department of Biochemistry and Molecular Biology, New Orleans, Louisiana, United States of America

  • Emma L. Cowles,

    Roles Data curation, Visualization, Writing – review & editing

    Affiliation Tulane School of Medicine, Department of Biochemistry and Molecular Biology, New Orleans, Louisiana, United States of America

  • Hua Lu,

    Roles Supervision, Writing – review & editing

    Affiliation Tulane School of Medicine, Department of Biochemistry and Molecular Biology, New Orleans, Louisiana, United States of America

  • James G. Jackson

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing

    jjacks8@tulane.edu

    Affiliation Tulane School of Medicine, Department of Biochemistry and Molecular Biology, New Orleans, Louisiana, United States of America

Abstract

Cancer cells survive chemotherapy and cause lethal relapse by entering a senescent state that facilitates expression of many phagocytosis/macrophage-related genes that engender a novel cannibalism phenotype. We used biosensors and live-cell imaging to reveal the basic steps and mechanisms of engulfment by senescent human and mouse tumor cells. We show filamentous actin in predator cells was localized to the prey cell throughout the process of engulfment. Biosensors to various phosphoinositide (PI) species revealed increased concentration and distinct localization of predator PI(4) P and PI(4,5)P2 at the prey cell during early stages of engulfment, followed by a transient burst of PI(3) P before and following internalization. PIK3C2B, the kinase responsible for generating PI(3)P, was required for complete engulfment. Inhibition or knockdown of Clathrin, known to associate with PIK3C2B and PI(4,5)P2, severely impaired engulfment. In sum, our data reveal the most fundamental cellular processes of senescent cell engulfment, including the precise localizations and dynamics of actin and PI species throughout the entire process.

Introduction

Cancer cells can survive chemotherapy and cause lethal relapse if they avoid cell death in conditions of (1) DNA damage and/or mitotic stress caused by therapy; and (2) nutrient deprivation. The breast cancer cells most likely to survive chemotherapy are TP53 wild type [1,2], and patients with these tumors have very poor survival [3]. We have previously shown that TP53 wild-type breast cancer cells survive chemotherapy by entering a senescent state that arrests the cell cycle to prevent mitotic catastrophe [4,5] and blocks apoptosis by activity of BCL-XL and/or MCL1 [6,7].

After entering senescence, tumor cells must still support elevated metabolism and survive an environment with unpredictable access to vasculature and nutrient sources. This should be problematic for senescent cells, as they have a high metabolic burden [812] that includes the need to produce cytokines and chemokines of the senescence-associated secretory phenotype (SASP) [4,5,1319]. SASP factors drive protumor phenotypes such as metastasis, stemness, and survival [15,16,2025]. Senescent cells have increased capacity for autophagy, but this is a limited source for nutrients [26].

We recently identified a novel phenotype of chemotherapy-induced senescent cells: the ability to engulf and break down entire cells in the lysosome [27]. Engulfment occurred after exposure of all 8 cultured cancer cell lines tested to different chemotherapy drugs and in vivo in a syngeneic mouse mammary tumor model. Chemotherapy-treated, senescent cells expressed an entire program of genes related to macrophages and phagocytosis. Engulfment of surrounding cells was unrelated to entosis and occurred at a very high frequency: 15% to 40% of senescent cells in culture had engulfed a neighboring cell. Engulfment conferred a survival advantage [27]. This phenotype had never been previously described and provides an explanation for many characteristics of the senescent cells that contribute to tumor relapse [27].

What is still not known is the mechanism of whole cell engulfment by senescent cells. Various types of engulfment, such as phagocytosis, efferocytosis, or macropinocytosis, are mediated by complex programs that share some features [28,29]. Here we show the program of whole cell engulfment is driven by actin dynamics and generation and localization of phosphoinositide (PI) species that regulate different stages of the engulfment process, and, surprisingly, we found that Clathrin is required for the process.

Results

Filamentous actin in predator cells localizes to prey during all stages of engulfment

Actin reorganization is required for fundamental cell processes, including motility and phagocytosis. To monitor dynamics of filamentous (F)-actin formation during senescent cell engulfment, we expressed the biosensor LifeAct-GFP [30], a fusion protein that marks F-actin, in MCF-7 cells that were previously shown to engulf when senescent [6,27]. LifeAct-GFP cells were treated with doxorubicin (DOXO), and interactions with mCherry-expressing, not-treated (NT) cells were cocultured (“DOXO-NT”) and visualized by live-cell imaging and confocal microscopy.

Over a time-course of 6 hours covering early stages of engulfment (contact, partial overtopping), confocal microscopy revealed F-actin concentrated in the predator cell at the points of contact with the prey cell, and at the leading edge of lamellar sheets during overtopping of the prey cell (Fig 1A and S1A and S1B Video). Following apparent internalization, F-actin concentrated and dissipated repeatedly around the prey cell(s) during the degradation process of approximately 7 hours (Fig 1B and S2A and S2B Video). F-actin formed structures reminiscent of those present intermittently on lysosomes and endosomes during propulsion through the cell or those associated with phagosomes and required for their mechanical deformation [31,32].

thumbnail
Fig 1. Filamentous actin in predator cells localizes to prey cells throughout the engulfment process.

Indicated “predator” MCF-7 cell lines were treated with 250 nM doxorubicin (DOXO) for 24 hours, washed, and plated on indicated not-treated (NT) MCF-7 “prey” cells, creating a “DOXO-NT” culture. Cells were imaged over days 4–8 as noted. (A) Volume view reconstruction of time course live-cell confocal imaging of a senescent LifeAct-GFP-MCF-7 cell during contact and partial overtopping of proliferating mCherry-MCF-7 cells on day 7 post-doxorubicin. Concentrated areas of F-actin are indicated by closed arrows, scale bar = 10 μm. (B) Volume view reconstruction of time course live-cell confocal imaging of a senescent LifeAct-GFP-MCF-7 cell following engulfment and through partial digestion of proliferating mCherry-MCF-7 cells. The top of the predator cell was removed above the 51.2% axial plane to visualize the internalized mCherry prey cell. Examples of concentrated F-actin are indicated by closed arrows, scale bar = 10 μm. (C-E) Axial slices along the Z-axis and volume view reconstruction of senescent LifeAct-GFP-MCF-7 cells during engulfment of a proliferating mCherry-MCF-7 cell. (C) shows contact and early overtopping, scale bar = 20 μm; (D) shows a mostly overtopped prey cell and a volume view with upper part of predator removed as in (B) scale bar = 20 μm. (E) shows a fully overtopped, engulfed cells with F-actin concentrated in lamellipodia and filipodia above and below the prey, scale bar = 25 μm. 3D volume view is shown from the bottom. (F) Volume view reconstruction of a senescent LifeAct-GFP-MCF-7 cell after fully engulfing a proliferating mCherry-MCF-7 cell. The top and side of the predator were removed along axial and sagittal planes to visualize engulfed cell, scale bar = 25 μm. (G, H) Confocal imaging with z-stack projections shown along right and bottom for senescent LifeAct-GFP-MCF-7 cells after fully engulfing proliferating mCherry-MCF-7 cells. At right of (G), volume views to visualize localization of the F-actin ring in the predator with an engulfed mCherry prey cell. Top right is top view, bottom right is side view, scale bar = 10 μm. At bottom of (H), axial slices displaying formation of smaller rings around digested portions of prey cell. Scale bar = 25 μm. (I) Schematic of a senescent LifeAct-GFP cell with 2 engulfed mCherry prey cells. Shown are example ROIs from the cell membrane, the cytoplasm, and the “ring” around the cell being engulfed. (J) Five individual LifeAct-GFP-MCF-7 predator cells were followed during engulfment of pHrodo red stained NT prey cells, and LifeAct-GFP intensity was measured at cytoplasm, nonengulfing membrane, and ring [as depicted in (I)] and pHrodo intensity was measured (right). The left graph shows the ratio of LifeAct-GFP intensity in the ring:cytoplasm (light green) and membrane:cytoplasm over 7 hours for 5 cells; error bars represent SEM. The right graph shows pHrodo red intensity over the same time course. Underlying data can be found at S1 Data. (K) Time course images of 1 representative cell used to generate data for (J). Top: merged green/red channels with ROI used for pixel intensity measurements marked in yellow. Middle: merged green/red channels. Lower: green channel. (L) Model of predator cell F-actin localization during the 5 stages of engulfing a prey cell.

https://doi.org/10.1371/journal.pbio.3001858.g001

To more precisely localize F-actin in senescent predator cells during various stages of engulfment, we visualized LifeAct-GFP in axial cross sections, 3D volume view reconstructions, and videos. In early stages, we found F-actin concentrated across a broad area of contact with the prey cell, including advancing filapodial structures (Fig 1C and S3 Video). F-actin in the predator cell was tightly associated with the prey cell throughout the overtopping process (Fig 1D and S4 Video). At the time of apparent internalization, F-actin was highly concentrated in filopodial extensions present under and over the top of the prey cell (Fig 1E and S5 Video). Following engulfment and during early digestion as the prey cell repositioned within the predator, F-actin localized to a concentrated “belt-like” structure around the circumference of the prey cell (Fig 1F–1H).

To more quantitatively examine F-actin dynamics during maturation of the phagosome, we stained prey cells with pHrodo, a pH-sensitive dye that fluoresces brightly in the acidic environment of the lysosome. Thus, an increase in pHrodo (red) fluorescence intensity can only occur following maturation and acidification of the phagosome after the prey cell has been internalized by the predator cell. We found that predator F-actin was localized at a ring-like structure associated with the prey cell (as depicted in Fig 1I) before engulfment and increased through engulfment, as pHrodo fluorescence increased (Fig 1J). Following the increase in pHrodo intensity as the prey cell was digested, the F-actin localized at the prey cell diminished sharply but was still maintained at a higher level than in the surrounding cytoplasm (Fig 1J, Ring Intensity/Cytoplasm Intensity). The fluctuating predator F-actin levels localized to the prey cell contrasted with levels observed in nonengulfing areas of the membrane edge, where F-actin was present, but intensity never diminished over the engulfment time-course (Fig 1J, Membrane Intensity/Cytoplasm Intensity). Fig 1K shows a time course of 1 representative LifeAct-GFP-expressing predator cell engulfing a pHrodo (red) stained prey cell. The predator “ring,” “membrane,” and “cytoplasm” areas used to calculate intensity are outlined and diagrammed in Fig 1I.

Similar dynamics were observed for 4226 mammary tumor cell line (S1 Fig). Taken together, the visualization of actin dynamics suggested 5 steps of the engulfment process: predator–prey contact, overtopping of the prey cell (ranging from approximately 10% to 90%), rounding of the prey cell, internalization, and degradation of the prey cell (Fig 1L).

Phospoinositide species regulate the entire process of whole cell engulfment by senescent cells

PI species interact with and modulate the activity of proteins that reorganize the actin cytoskeleton and are thus critical for the control of actin dynamics in macrophages during phagocytosis of various large targets [3338].

We used specific biosensors to investigate involvement of different PI species in engulfment by senescent cells [39]. Live-cell time course imaging showed PLCD1-GFP [detecting PI(4,5)P2] and SidMx2-GFP [detecting PI(4)P] in the predator cell were closely associated with the prey cell during initial contact, at the overtopping stage, and through the point the predator cell became rounded and was detached from the substrate (Fig 2A and 2B). Both species dissipated in concentration as the digestion process ensued (Fig 2A and 2B). We observed in the predator cell a transient burst of highly concentrated 2xFYVE-GFP [detecting PI(3)P] enveloping the prey after overtopping that remained through rounding and internalization (Fig 2C). Imaging showed biosensors detecting PI(3,4)P2, PI(3,4,5)P3, and a mutant PLCD1-GFP biosensor that does not associate with any PI species (Negative Ctrl) were not concentrated or only diffusely associated with the engulfed cell (Fig 2D–2F), suggesting these species are not involved or are below the sensitivity of detection of their biosensors.

thumbnail
Fig 2. PI biosensors expressed in predator cells reveal PI(4,5)P2, PI(4)P, and PI(3)P, but not PI(3,4,5)P3 or PI(4)P, localize to the prey cells during engulfment.

MCF-7 cell lines expressing biosensors that detect indicated PI species were treated with 250 nM doxorubicin for 24 hours, washed, and plated on untreated mCherry-MCF-7 “prey” cells. DOXO-NT cultures were imaged over days 3–8. Time course live-cell imaging of senescent MCF-7 cells that express (A) PLCD1-GFP marking PI(4,5)P2; (B) P4M-SidMx2-GFP marking PI(4)P; (C) 2xFYVE-GFP marking PI(3)P; (D) TAPP1-GFP marking PI(3,4)P2; (E) BTK-GFP marking PI(3,4,5)P3; (F) PLCD1(R40L)-GFP mutant that does not bind PI species (Negative Ctrl), throughout the entire process of engulfing mCherry-MCF-7 cells. Scale bar = 100 μm.

https://doi.org/10.1371/journal.pbio.3001858.g002

To more precisely localize PI species in senescent predator cells during various stages of engulfment, we visualized the biosensors specific for PI(4,5)P2, PI(4)P, and PI(3)P, in axial cross sections, 3D volume view reconstructions, and videos. We found that PI(4,5)P2 was concentrated at the leading edge of the membrane as it contacted prey cells (Fig 3A, upper) and in the cell membrane as it advanced over and around the prey (Fig 3A, lower, and S6 Video). PI(4)P was similarly concentrated at contact early in engulfment (Fig 3B, upper, closed arrow), as the prey was partially overtopped (Fig 3B, upper, and S7 Video), and after complete overtopping and rounding of the prey cell (Fig 3B, lower). Fig 3C shows GFP fused to a mutant PI biosensor varies minimally in concentration within a predator and is only diffusely associated with prey cells during early through late stages of engulfment (Fig 3C). PI(4,5)P2 and PI(4)P both localized to the cell membrane and points of contact at the nonengulfing edges of the cell (Fig 3A and 3B). To determine if these PI species concentrate to a greater extent at sites of engulfment, we quantified the intensity at different regions of the engulfing cell (as depicted in Fig 1I) and compared to another membrane localizing protein that is not involved in engulfment, LYN11 [40,41]. We found the ratio of PLCD1-GFP to LYN11 was greater at the engulfing edge of the phagocytic cup than the nonengulfing edge or in the cytoplasm, suggesting PI(4,5)P2 was preferentially concentrated at sites of engulfment (Fig 3D). Similar results were observed for PI(4)P as detected by SidMx2-GFP (Fig 3E). Quantification over time lapse images showed PI(4,5)P2 dissipated at the prey cell before pHrodo fluorescence, suggesting PI(4,5)P2 was not involved in phagosome maturation post-engulfment (Fig 3F and 3G).

thumbnail
Fig 3. Confocal imaging shows precise localizations of PI(4,5)P2, PI(4)P, and PI(3)P in predator cells during engulfment.

Axial slices along the Z-axis and volume view reconstructions of senescent MCF-7 cells, expressing different PI biosensors, during different stages of engulfment in DOXO-NT cultures. (A) PLCD1-GFP localization in predator cells during partial (top) and complete (bottom) overtopping of an mCherry-MCF-7 cell (top), scale bar = 20 μm (bottom), scale bar = 10 μm; (B) P4M-SidMx2-GFP localization in a predator cell during partial (top) and complete (bottom) overtopping of an mCherry-MCF-7 cell, scale bar = 10 μm. Closed arrows indicate examples of concentrated P4M-SidMx2-GFP localized to lamellipodia advancing over the prey cell. (C) Localization of PLCD1(R40L)-GFP mutant that does not bind PI species, during early (closed arrow), mid (dashed arrow), and late stage (open arrow) engulfments of mCherry-MCF-7 cells, scale bar = 10 μm. (D) MCF-7 cells coexpressing PLCD1-GFP and LYN11-mCherry (a control that localizes to membranes but is not involved in phagocytic processes) were imaged while engulfing NIR-MCF-7 cells. Left: representative images of 3 cells and the ROI used for quantitation indicated. Right: pixel intensities for the indicated ROI were determined and the ratio of intensity for PLCD1-GFP to LYN11-mCherry was calculated for 2–3 areas of the membrane edge actively contacting the prey cell (engulfing edge), 2–3 areas of the uninvolved membrane (nonengulfing edge), and 1 measurement per cell for the whole cytoplasm. Ratio calculations for 6 individual cells are shown as violin plots with mean and SEM indicated by red dashed line and black dashed line, respectively. Underlying data can be found at S1 Data. (E) Ratios as in (D) were calculated and shown for MCF-7 cells coexpressing SidMx2-GFP and LYN11-mCherry. Underlying data can be found at S1 Data. (F) Five PLCD1-GFP-MCF-7 predator cells were followed during engulfment of pHrodo red stained prey cells, and PLCD1-GFP intensity was measured at cytoplasm, membrane, and ring (as depicted in Fig 1I) and pHrodo intensity was measured (right). The left graph shows the ratio of PLCD1-GFP intensity of the ring:cytoplasm (light green) and membrane:cytoplasm over 12 hours for 5 cells; error bars represent SEM. The right graph shows pHrodo red intensity +/− SEM over the same time course. Underlying data can be found at S1 Data. (G) Time course images of 1 representative cell used to generate data for (F). Top: merged green/red channels with ROI used for pixel intensity measurements marked in yellow. Middle: merged green/red channels. Lower: green channel. Scale bar = 100 μm. In Fig 4D and 4E, two-way ANOVA was used for analysis, p-value < 0.05 = *, p-value < 0.01 = **, p-value < 0.001 = ***.

https://doi.org/10.1371/journal.pbio.3001858.g003

Confocal imaging of 2xFYVE-GFP expressing cells revealed striking images of PI(3)P completely enveloping the prey after rounding (Fig 4A, upper and lower, S8 Video). Quantification of fluorescence intensity in different channels across the predator/prey showed predator cell PI(3)P concentrated at the prey cell far exceeding levels observed for LYN11-mCherry (Fig 4B and 4C). Time lapse imaging with pHrodo-labeled prey cells showed 2xFYVE-GFP was highly concentrated during the transient burst and remained enriched during digestion of the prey cell (Fig 4D and 4E).

thumbnail
Fig 4. Confocal imaging shows precise localization of PI(3)P in predator cells during engulfment.

(A) Axial slices along the Z-axis and volume view reconstructions of a senescent MCF-7 cell expressing 2xFYVE-GFP and an mCherry expressing MCF-7 prey cell in a DOXO-NT culture. Top panel scale bar = 10 μm; bottom panel scale bar = 1 μm. (B, C) Representative image of a 2xFYVE-GFP cell engulfing a NIR-MCF-7 cell. Pixel intensity was determined across the distance of the yellow line and plotted in (C) for 5 individual cells. Scale bar = 60 μm. Underlying data can be found at S1 Data. (D) Five 2xFYVE-GFP-MCF-7 predator cells were followed during engulfment of pHrodo red stained prey cells, and 2xFYVE-GFP intensity was measured at cytoplasm, membrane, and ring (as depicted in Fig 1I) and pHrodo intensity was measured (right). The right graph shows the ratio of 2xFYVE-GFP intensity of the ring:cytoplasm (light green) over multiple hours surrounding pHrodo intensity increase for 5 cells. Error bars represent SEM. Underlying data can be found at S1 Data. (D) Representative images of one 2xFYVE-GFP-MCF-7 predator cell engulfing pHrodo-labeled MCF-7 prey used for Fig 4D (right). ROIs used for quantification are shown in top panel. Scale bar = 100 μm. (E) Time course of a 2xFYVE-GFP-MCF-7 cell engulfing an MCF-7-mCherry prey cell. Dashed arrows indicate predator/prey contact; closed arrows indicate the overtopped prey cell with 2xFYVE concentration; open arrows indicate the prey cell no longer overtopped mCherry prey cell. Scale bar = 50 μm. (F) Axial slices along the Z-axis and volume view reconstructions of senescent 2xFYVE-GFP-MCF-7 predator cells engulfing mCherry-MCF-7 prey cells. Leftmost image shows the bottom axial slice. Closed arrows indicate 2xFYVE-GFP concentration at the prey cell. Scale bar = 10 μm.

https://doi.org/10.1371/journal.pbio.3001858.g004

In relatively rare instances, predator PI(3)P was observed concentrating at prey cells that were ultimately never engulfed (Fig 4F), suggesting PI(3)P localization was occurring before internalization. Indeed, confocal imaging showed PI(3)P tightly concentrated at prey cells that were clearly still attached to the substrate (Fig 4G). These data contrast to macrophages that localize PI(3)P only following phagocytic cup closure and during phagosome maturation [35,4245], suggesting a novel role for PI(3)P generation in engulfment by senescent cells.

We observed accumulation of PI(4,5)P2, PI(4)P, PI(3)P at multiple stages of engulfment (Figs 24). To precisely discern differences in their localization and timing during engulfment, we coexpressed biosensors for PI(4,5)P2 and either PI(4)P (Fig 5A) or PI(3)P (Fig 5B) in GFP/mCherry combinations in the same predator cell, and imaged engulfment of near-infrared fluorescent protein (miRFP-713, “NIR”) expressing prey cells. Time-lapse imaging showed near perfect overlap of predator PI(4,5)P2 and PI(4)P during contact and through overtopping of prey cells (Fig 5A, NIR prey cells are blue). Confocal imaging in volume view, of axial slices, and in separated color channels also clearly showed colocalization of PI(4,5)P2 and PI(4)P in a predator cell at different stages of engulfing multiple prey cells (Fig 5C, NIR prey cells are pink). These data suggest similar or identical roles in engulfment and that phosphorylation of the 4 position on the inositol ring is significant.

thumbnail
Fig 5. PI(4,5)P2, PI(4)P colocalize at contact through overtopping of prey, then dissipate as PI(3)P localizes tightly around prey during later stages.

(A) Time course live-cell imaging of a senescent MCF-7 cell that expresses PLCD1-mCherry and P4M-SidMx2-GFP, throughout the entire process of engulfing NIR-MCF-7 prey cells in a DOXO-NT culture. Closed arrows indicate examples PLCD1-mCherry and P4M-SidMx2-GFP colocalization during engulfment. (B) Time course live-cell imaging of a senescent MCF-7 cell that expresses PLCD1-mCherry and 2xFYVE-GFP, throughout the entire process of engulfing a NIR-MCF-7 cell. Closed arrows indicate localization of PLCD1-GFP, but not 2xFYVE-mCherry, at mid stage engulfment. Open arrows indicate 2xFYVE-mCherry localization, but not PLCD1-GFP, at late-stage engulfment. Scale bar = 100 μm for A+B. (C) Volume view reconstruction of a senescent MCF-7 cell expressing PLCD1-GFP and P4M-SidMx2-mCherry that is engulfing 3 NIR-MCF-7 cells. Middle panels show separate color channels; bottom panel shows axial planes from bottom to top. Closed arrows indicate examples PLCD1-mCherry and P4M-SidMx2-GFP colocalization during engulfment, scale bar = 10 μm. (D) Volume view reconstruction of a senescent MCF-7 cell expressing PLCD1-mCherry and 2xFYVE-GFP, which is engulfing 2 NIR-MCF-7 cells. Middle panels show separate color channels; bottom panels show axial planes from bottom to top. Closed arrows indicate localization of PLCD1-mCherry, but not 2xFYVE-GFP, at mid stage engulfment. Open arrows indicate 2xFYVE-mCherry localization, but not PLCD1-GFP, at an internalized cell, scale bar = 10 μm.

https://doi.org/10.1371/journal.pbio.3001858.g005

Further imaging showed PI(4,5)P2 and PI(3)P were localized at entirely different stages of engulfment (Fig 5B and 5D). We found in the predator cell that PI(4,5)P2 was concentrated at the prey cell first, during contact and overtopping, but then dissipated as PI(3)P began to concentrate and localize. We observed very little overlap or colocalization during the transition (Fig 5B). This pattern of localization was evident in confocal imaging in a predator cell engulfing 2 NIR expressing prey cells (pink): one after internalization and one at overtopping. Imaging shows PI(4,5)P2 was concentrated and localized at the overtopped prey cell (Fig 5D) but not the internalized prey cell, while PI(3)P was concentrated and localized to the internalized, rounded prey cell (Fig 5D, green localization, open arrow). Similar PI species localizations were observed in 4226 cell line (S2 Fig). Dynamic localization of different PI species throughout the process of engulfment is depicted in S3 Fig.

PIK3C2B facilitates internalization

Based on imaging data (Figs 25), the PI species most likely to play a nonredundant, critical role for the final, irreversible stages of engulfment is PI(3)P. This species is generated by 4 different kinases [46], and we used CRISPR-Cas9 to derive single-cell knockout clones of each (Figs 6A and S4A). Engulfment rates determined on IncuCyte (Materials and methods, S6 Fig), and viability of predator cells with knockout of PIK3C2G, PIK3C3, PIK3C2A were not reduced (S4B Fig). PIK3C2B has a well-described role in cancer biology, therapy resistance, cell motility, and cytoskeletal dynamics [4754]. We found that cells deficient in PIK3C2B were relatively normal and were not more sensitive to the chemotherapy (Fig 6B, right, and S4D Fig, right). We found multiple PIK3C2B knockout clones had reduced rates of whole cell engulfment (Fig 6B, left, and S4D Fig, left). In PIK3C2B knockout cell lines, we observed that predators formed lamellar projections, made contact with and successfully overtopped prey (S4E and S4F Fig), but were less likely to complete engulfment of adherent prey (Fig 6B).

thumbnail
Fig 6. Predator cell PIK3C2B colocalizes with PI(3)P during engulfment and is required for complete internalization.

(A) CRISPR-Cas9 single-cell clone PIK3C2B knockouts of MCF-7 were screened by immunoblot. Clones indicated by asterisk were chosen for further testing. (B) Predator cell engulfment rates (left) and confluency as a measure of viability (right) for senescent MCF-7 parental cells and 3 PIK3C2B knockout clones were determined by time course imaging in a DOXO-NT culture. Underlying data can be found at S1 Data. (C) Time course live-cell imaging of a senescent MCF-7 cell expressing 2xFYVE-mCherry and PIK3C2B-GFP throughout the entire process of engulfing an NIR-MCF-7 cell, scale bar = 100 μm. (D) Left: Merged (top) and separate (middle, lower) color channels of axial planes from bottom to top of a senescent MCF-7 cell expressing 2xFYVE-mCherry and PIK3C2B-GFP, engulfing NIR-MCF-7 cells. Right, volume view reconstruction of the same image, scale bar = 10 μm.

https://doi.org/10.1371/journal.pbio.3001858.g006

Time course live cell imaging showed PIK3C2B-GFP fusion protein began to accumulate near prey cells during overtopping and increased in concentration through rounding and internalization of the prey and was partially colocalized with PI(3)P during much of this process (Fig 6C). Confocal microscopy showed both PIK3C2B-GFP and PI(3)P were present independently in the predator and did not completely overlap in localization, but significant areas of overlap were observed at the site of an overtopped prey cell (Fig 6D). In both time-lapse and confocal imaging, we noted the cytoplasmic PIK3C2B protein was diffusely present and did not associate as tightly to prey cells as the biosensors that detected membrane bound PI species.

PI3 kinase activity and generation of PI(3,4,5)P3 plays a critical role in phagocytosis of large substrates [55, 56], and thus it was surprising to observe no localization during senescent cell engulfment (Fig 2). To determine if PI(3,4,5)P3 contributes to engulfment in our senescent cells without actually localizing to sites of engulfment, we treated cells with 2 selective PI3 kinase inhibitors, BKM-120 (targeting p110α/β/δ/γ) [57] and GDC-0941 (targeting p110α/δ) [58]. Because the entire culture of predator/prey is exposed to drug, any observed change in engulfment could be due to drug effects on the proliferating prey cell, and not the senescent predator cell. To mitigate this, we assayed engulfment of senescent predator/senescent prey cells (S5 Fig, DOXO-DOXO) in cultures as we have previously [27].

Using concentrations that were mostly nontoxic to the cells (S5A and S5B Fig), we observed no decrease in predator engulfment rates for MCF-7 or 4226 mammary tumor cells (S5C and S5D Fig). MCF-7 cells were more sensitive to the drugs (S5C Fig), but at concentrations that were beginning to cause cell death (S5C Fig, lower), more engulfment was actually observed (S5C Fig, upper). It is possible this could be the engulfment of dead/dying cells through efferocytosis by the senescent predator cells, as has been previously observed [27]. Nevertheless, we did not observe reduced engulfment rates in MCF-7 or 4226 cell lines treated with potent inhibitors of PI3 kinase. Taken together with localization data (Fig 2), these data suggest generation of PI(3,4,5)P3 is not critical for senescent cells to engulf other whole cells.

Clathrin is involved in senescent cell engulfment

After finding evident involvement of PIK3C2B in senescent cell engulfment, we investigated potential downstream mechanisms. Interestingly, PIK3C2B directly associates and colocalizes with Clathrin [59], a protein known to create cages for endocytosis and facilitate nucleation of F-actin during engulfment of substrates ranging from small to large [60,61]. Components of the Clathrin complex also associate with PI(4,5)P2, which we found highly enriched at sites of engulfment (Figs 24) [62,63]. In senescent cells not actively engulfing or contacting other cells, Clathrin was localized to its typical, perinuclear space (Fig 7A). When a prey cell was overtopped by a predator cell, Clathrin concentrated and localized around and over the prey, colocalized with PIK3C2B (Fig 7A, closed arrows). Following internalization, as the prey cells were digested, Clathrin was redistributed back to the perinuclear space (Fig 7A, open arrows). Next, we treated cells with 3 drugs that all have in common inhibition of Clathrin-mediated endocytosis (CME) [64] (Dynasore [65], Filipin3 [66], MBCD [67]). While these drugs are not highly specific to Clathrin, we did observe reduced engulfment rates at the day 8 posttreatment time point (Figs 7B and S7A). Filipin3 treatment decreased engulfment by 23.7% and 61.8% in MCF-7 and 4226 cells, respectively. Inhibition of dynamin using Dynasore decreased engulfment by 52% in MCF-7 cells and by 58.1% in 4226 cells. MBCD reduced engulfment by 74% in MCF-7 and 55.9% in 4226 (Figs 7 and S7A). Though we identified doses of drugs that were toxic to predator and prey cells (Fig 7B, red lines), engulfment was still inhibited at a wide range of well-tolerated concentrations (Figs 7B and S7A). To rule out drug effects on prey cells (rather than on predator cell engulfment mechanisms), we treated cultures of senescent predator/senescent prey (DOXO-DOXO) cells, and, in these conditions, similar reductions in engulfment rate were observed (S7B Fig). Because, in addition to Clathrin, these drugs can inhibit other proteins, we used shRNA to knock down CLTC (S7C Fig). CLTC knockdown dramatically reduced rates of engulfment (71.5% in MCF-7 and C60.6% in MPE600), supporting the results with the CME inhibitors (Figs 7B and S7D, shCLTC).

thumbnail
Fig 7. Clathrin relocalizes to sites of engulfment and is required for efficient completion of internalization.

(A) Time course live-cell imaging of a senescent MCF-7 cell expressing CLTC-mCherry and PIK3C2B-GFP, throughout the entire process of engulfing an NIR-MCF-7 prey cell in a DOXO-NT culture. Scale bar = 100 μm. (B) Predator cell engulfment rates (upper) and confluency (lower) in DOXO-NT cultures were determined for senescent GFP-MCF-7 cells that were treated with indicated drugs or made to express an shRNA construct targeting CLTC. Concentrations of Dynasore, Filipin3, and MBCD producing toxicity are displayed as red lines (Dynasore 75 uM, 100 uM and MBCD 6 mg/mL). Underlying data can be found at S1 Data.

https://doi.org/10.1371/journal.pbio.3001858.g007

Discussion

Our group previously showed breast tumor cells that enter senescence to avoid apoptosis and mitotic catastrophe after chemotherapy express many genes related to different mechanisms of engulfment [27]. These include genes related to phagocytosis, efferocytosis, macropinocytosis, and endocytosis, thus resulting in a novel phenotype that includes whole cell cannibalism. Here, we determined the fundamental mechanisms mediating the process by which a senescent tumor cells cannibalizes a healthy neighboring cell. We defined 5 basic steps in the process (Fig 1L), and we show F-actin is involved throughout the process, from contact to relatively late in digestion. Recent studies have described 2 distinct modes of actin-mediated macrophage engulfment, lamellipodial and filopodial, depending on the spatial constriction or crowding of the cells [68]. As senescent predator cells contacted and overtopped prey cells, we observed formation of both lamellar sheets and elongated filopodial extensions that were rich in F-actin (Fig 1). However, the predominant type of engulfment in senescent breast cancer cells appears to be the sheet-like overtopping shown in Fig 1A. As fully engulfed cells were processed, F-actin was highly concentrated around the circumference of the engulfed cell. As the prey cell became smaller and broken down from one large vesicle to multiple smaller ones, concentrated F-actin structures formed, disappeared, and reformed continuously (Fig 1B and S2 Video). These observations are highly reminiscent of transient actin structures formed in macrophages that mediate mechanical breakdown (described as “chewing”) of internalized bacteria or red blood cells [32,69].

PI species mediate numerous cell processes, including engulfments such as phagocytosis. During engulfment, PI(4)P and PI(4,5)P2 become concentrated and precisely localized at the prey cell. The fact these 2 species appear to have overlapping localization suggests possible redundancy that could, perhaps, extend to the kinases responsible for generating them. PI(3)P was the only species detected at the critical stages of rounding and internalization, and at least 1 kinase responsible for generating PI(3)P was required for efficient engulfment. Distinct from other forms of phagocytosis [35], we observed enriched localization of predator PI(3)P at prey cells before internalization (Fig 4D–4G). Further, while PI3 kinase activity and generation of PI(3,4,5)P3 is critical for phagocytosis of large substrates by macrophages [55,56], in senescent cells, PI(3,4,5)P3 did not localize to sites of engulfment, and phagocytosis occurred efficiently in the presence of PI3 kinase inhibitors (S5 Fig).

While our cells with knockout of PIK3C2B had reduced rates of engulfment, we cannot rule out involvement or compensation by the other kinases responsible for generating PI(3)P: PIK3C3, PIK3C2A, and PIK3C2G [46]. Knockout of these kinases, however, had no effect on engulfment, suggesting a primary role for PIK3C2B. It is also possible the kinase(s) necessary and/or sufficient for PI(3)P generation are tissue or cell type specific.

Identifying Clathrin involvement was surprising but also made sense in light of evident roles for PI(4,5)P2 and PIK3C2B, both known to interact with Clathrin [59]. Roles for Clathrin in mediating engulfment of larger targets are emerging [61], and Clathrin involvement in engulfment by senescent cells may be related to its role as an actin network organizer rather than classic, coated pit endocytosis [61]. Recently, Clathrin was also shown to be involved in the final stages of phagocytosis: resolution of phagosomes [70].

In sum, we present the first data demonstrating the mechanistic underpinnings of the remarkable cell cannibalism phenotype of breast tumor cells that enter senescence to survive chemotherapy. Understanding how senescent cells engulf could inform strategies to eliminate them, thus improving treatments in the cancer most difficult to eradicate.

Materials and methods

Cell culture

MCF-7 cells were from ATCC and cultured in complete EMEM media (ATCC). MPE600 cells were a gift of Joe Gray (Oregon Health Sciences University) and were cultured in complete DMEM media. 4226 cells, created from a spontaneous MMTV-Wnt1 tumor, were cultured in complete DMEM and have been previously described [6,27].

Vector construction

Plasmids for tracking PI species were from Addgene. Specifically, mCherry-Clathrin LC-15 (Addgene #55019) [71]; PLCD1(PH)-mCherry (Addgene #36075) [72]; PH-PLCD1-GFP (Addgene #51407), PH-PLCD1(R40L)-GFP (Addgene #51408) [73]; pLentiLifeACT-EGFP BlastR (Addgene #84383) [74]; PH-Btk-GFP (Addgene #51463), GFP-PH-TAPP1 (Addgene #161985); and pLenti-EGFP-P4M-SidMx2 (Addgene #136997), pLenti-EGFP-2xFYVE (Addgene #136996) were gifts from Ken-Ichi Takemaru. Plasmids were also constructed by Vectorbuilder: VB200917 is a lenti-mCherry-shCLTC containing U6-driven shRNA sequence 5′-CGTGTTCTTGTAACCTTTATT-3′. VB200917 is a lenti-mCherry-shLuciferase containing U6-driven shRNA sequence 5′- ATGTTTACTACACTCGGATAT-3′ and was used to create lentivirus to infect MCF7, 4226, and MPE600 cell lines. Infected cells were sorted for mCherry using FACS as above and tested for engulfment. Human PIK3C2B-GFP fusion overexpression vector is a lentiviral CMV-driven plasmid from VectorBuilder. miRFP-713 (near-infrared fluorescent protein) overexpression vector is a lentiviral EFS-driven plasmid from VectorBuilder. 2xFYVE-domain-mCherry fusion protein is lentiviral plasmid from VectorBuilder. Lentiviral EFS-driven LYN11 mCherry fusion protein construct was also produced by VectorBuilder.

Lentiviral infections

Lentivirus was produced using 293T cells and a standard Lipofectamine 2000 protocol according to the manufacturer’s instructions (Thermo Fisher Scientific). Viral supernatants were collected 48 to 72 hours after transfection and filtered (0.45 μm). MCF-7 cells were plated at 166,000 cells in a 6-well plate, then spinfected with the filtered viral supernatant at 1,100 x g for 20 minutes the following day. Positive cells were FACS sorted using a BD FACSAriaII as previously published [5].

CRISPR-mediated gene deletion

Guides for CRISPR-mediated gene editing (human target sequences: PIK3C2Bsg1 ATGGCGTATGTCATCAACG and PIK3C2Bsg2 CATGGCGTATGTCATCAAC) (Mouse targeting sequences: Pik3c2bsg2 GACCTGTGGTCGAAAGCTCC, Pik3c2bsg3 ATCAAGTTCACCAGCCCCCC) were designed to target regions of the PI3Ka domain of PIK3C2B. These guides were hybridized to their antisense oligos and ligated into BsmBI-digested pLentiCRISPRV2-mCherry [Addgene #75161] backbone. Resultant vectors were sequence verified by sanger sequencing via Genewiz. Lentivirus was created as described above and used for infection of MCF7. Infected cells were sorted for mCherry using a BD FACSAriaII cell sorter and serially diluted giving rise to single-cell clones. Individual clones were grown out and tested for PIK3C2B expression by western blot. Clones confirmed to have homozygous deletions of PIK3C2B were then used for engulfment assays. Single-cell clones with knockouts of PIk3C2A, PIK3C3, and PIK3C2G in MCF7 were created using the same backbone and targeting strategy. The target sequences used for PIK3C2A were PIK3C2Asg1 GGATCCTCTAAGTAAGCCTA, PIK3C2Asg2 CAATATAATGCTTCGAAGCA. The target sequence used for PIK3C3 was PIK3C3sg3 GGCTGAAGTTCTTGATACAG. The target sequences used for PIK3C2G were PIK3C2Gsg2 AATGGCGATGAATCTTTGCT, PIK3C2Gsg3 TGTGACGGTCACATACTCCC.

Drug treatments

Doxorubicin was from Santa Cruz Biotech (SC-200923A), resuspended in water at a concentration of 10 mg/mL, and then diluted in PBS to 50 μg/mL and used at concentrations indicated. Filipin III (Filipin3) was from Sigma (Cat# F4767). Dynasore was from Abcam (Cat# 120192). Methyl-B-Cyclodextrin (MBCD) was from Sigma-Aldrich (Cat# 332615). Filipin III, Dynasore, or MBCD were added to cultures at concentrations shown immediately after removing doxorubicin. BKM-120 was purchased from Apex Bio (Cat#A3015, Batch#2) and resuspended in DMSO for a stock concentration of 2.44 mM. GDC-0941 was purchased from Apex Bio (Cat#A8210, Batch#2) and resuspended in DMSO for a stock concentration of 1.95 mM. The drug/media were changed every 2 days until timelines were completed. pHrodo red labeling kit for phagocytosis was purchased from Sartorius (Cat#4649).

Immunoblotting

Cell lines were treated and harvested as indicated in the figures. Lysates were prepared using TNESV buffer (1% NP-40, 50 mM Tris-HCl (pH 7.5), 100 mM NaCl, 2 mM EDTA, and 10 mM NaVO4), and 20 μg of total protein was separated by 10% SDS-PAGE and then transferred to nitrocellulose membranes using Bio-Rad reagents and apparatus. Immunoblotting was performed as previously described [4] for PIK3C2B (goat polyclonal, AF7249; R&D Systems), CLTC (rabbit monoclonal, 4796S (D3C6); Cell Signaling), and actin (mouse monoclonal, MA5-15739; Invitrogen). All antibodies were used at 1:1,000 dilution.

Live-cell imaging

IncuCyte S3 using software version 2019B and IncuCyte Sx5 using software version2021B were utilized for live-cell visualization of GFP, mCherry, and miRFP-713 constructs in cells. Images were taken every 2 hours for at least 25 images per well at 20× magnification from day 3 to day 8 following doxorubicin treatment.

Engulfment quantification

IncuCyte S3 and SX5 live cell imagers were used to quantify engulfment. Images were taken every 2 hours from the 3D-7D post-doxorubicin time point. Parameters were optimized to quantify images as follows: First, thresholds were set on predator cells for fluorescence intensity (to exclude background fluorescent areas) and size (large size restriction used to exclude dead cells). Once these constraints were set, we established settings for quantifying the number of engulfed prey cells. As the prey cells are made smaller and rounded when fully internalized by the predator cells, we adjusted size and circularity parameters to account for these circumstances. Cell outline area setting was adjusted from the default (1) on both predators and prey to exclude normally occurring cell–cell contacts. The areas of predator/prey overlap were then quantified as a function of the number of predator cells, yielding the percentage of predator cells engulfing at each time point. Finally, to minimize counting overlaps due to floating cellular debris and cells that were overtopped but never internalized, we averaged the counts from 3 × 2-hour time points, yielding an average percentage engulfment for each 6-hour time point. For MCF-7: Segmentation Adjustment: 1 (default), Red MCF7 Predators: Tophat radius set to 45 μm, Threshold (OCU): 0.4, Edge sensitivity at −35, Hole fill: 60 μm2, Pixel size: 5, Area inclusion: 160 to 20,000 μm2, Eccentricity maximum: 0.99; Green MCF7 Prey: Tophat radius set to 10 μm, Threshold (GCU): 0.15, Edge sensitivity at −25, Hole fill: 500 μm2, Pixel size: −3, Area minimum 10 μm2, Eccentricity inclusion 0.4 to 0.9. For MPE600, settings were the same except for the following: Predator pixel size was increased to 5 as they had less background cell/cell contact; Prey Tophat radius was increased to 10 μm; Prey threshold (GCU) was increased to 0.75; Prey pixel size was also increased to 2, Area: maximum set to 80 μm2. For 4226 cells, Segmentation Adjustment: 2, Red 4226 Predators: Tophat radius: 75 μm, Threshold (OCU): 0.15, Edge sensitivity at −40, Hole fill: 400 μm2, Pixel size to stop edge effect −5, Area inclusion: 324 to 20,000 μm2, Green 4226 Prey: Tophat radius set to 25 μm, Green threshold 2.5 GCU, Edge sensitivity at −35, Hole fill: 400 μm2, Pixel size to stop edge effect −3, Area inclusion 18 to 160 μm2, Eccentricity inclusion 0.3 to 0.96, Green threshold 2.5 GCU, Red + Green (Overlap) Parameters: Minimum Area: 5 μm. For each engulfment experiment, as a measure of predator viability over time, surface area of the predator cells was determined.

Quantitative analysis of biosensor localization

MCF7 predator cells expressing both biosensor-GFP and LYN11-mCherry were mixed with MCF7 prey expressing miRFP713 and imaged on IncuCyte live-cell imager. Images at specific stages of engulfment were exported as separate channels (GFP-biosensor, Orange-Lyn11, NIR-prey).

The images were imported into ImageJ and converted to greyscale. ImageJ line tool was used to manually quantify PI(3)P (as in Fig 4B), and polygon tool was used for creating ROI for quantification of PI(4)P and PI(4,5)P2 (as in Figs 1I, 3D and 3E) ROI measurements from each cell were normalized to the lowest mean intensity (LMI) for the GFP channel from the same cell, producing enrichment for the GFP channel at the indicated locations (i.e., nonengulfing membrane, cytoplasm, ring forming at the prey cell, as shown in Fig 1I). Both biosensor-GFP and LYN11-mCherry intensities were independently normalized to their respective LMI within the cell. The normalized biosensor-GFP intensity was then divided by normalized LYN11-mCherry intensity, giving the overall enrichment of the biosensor species relative to LYN11. Examples of these regions are shown in Figs 1K (upper), 3D, 3E and 4C. To generate graphs for PI(3)P, 10 independent engulfing cells were quantified. For PI(4)P and PI(4,5)P2, 6 independent engulfing cells were quantified.

pHrodo time-lapse quantification

pHrodo-labeling of prey cells was preformed (in accordance with manufacturer’s instructions) on MCF-7 cells prior to incubation with senescent MCF7 predator cells. In brief, MCF-7 prey cells were plated (2.5 million cells in 10 cm dish) and treated with 0.25 uM doxorubicin. Fresh media was added 24 hours later. After 3 days, MCF-7 prey cells were trypsinized and counted, 100,000 cells were washed with 1 mL of prewarmed “IncuCyte pHrodo wash buffer,” then briefly rinsed in 1 mL of prewarmed “IncuCyte pHrodo labeling buffer.” Cells were then spun down at 400 RCF × 5 min and resuspended in 1 mL of pHrodo red labeling dye at a concentration of 100 ng/ml in IncuCyte pHrodo labeling buffer and incubated at 37°C for 1 hour. Resulting pHrodo-labeled cells were pelleted, rinsed with PBS, resuspended in growth media, and plated on senescent (day 3 post-DOXO) MCF7-GFP cells at a ratio of 1 predator cell:4 prey cells. Quantification of biosensor-GFP predator cells engulfing pHrodo-labeled prey over time was preformed using ImageJ at time points centered around the point of greatest increase in pHrodo fluorescence between 2 time points over the time course. Three regions in the cytoplasm were measured and averaged for each image. The cell membrane and predator/prey contact points were manually outlined using the polygon tool and normalized to mean cytoplasm ROI intensity. Time course quantification was performed for 5 cells. Data points plotted represent mean intensity ratios for ROI indicated (e.g., ring intensity/cytoplasm intensity).

Microscopy

Confocal images of fluorescent cells were captured on a Nikon TiE-2 Inverted Research Microscope Nikon A1rSi Laser Point Scanning Confocal, with a Plan Apo λ 60× oil objective for both cell lines and tissue sections. Image brightness and background levels were adjusted uniformly and minimally for entire images using NIS-Elements software or Photoshop when necessary for visual clarity in publication. In the case of experiments with miRFP713 (near-infrared) constructs, imaging was acquired in the Cy5 channel and spectrally unmixed in Nikon’s NIS elements software version 5.02.01 with default settings.

Statistics

In all IncuCyte experiments with cellular engulfment quantified, 25 fields of view per well were taken for 3 wells of each treatment group and 3 wells of the control group every 2 hours. Two biological replicates of each experiment were performed. The percentages from 3 consecutive time points (6 hours) were then averaged. For IncuCyte experiments containing 2 conditions, unpaired t tests with 2-stage step-up (Benjamini, Krieger and Yekutieli) model was performed. A minimum of 2,000 predator cells were counted for each 2-hour time point. Engulfment percentages were compared for the final time point of the engulfment experiments. For IncuCyte experiments containing 3 or more conditions, 2-way ANOVA with a Dunnett’s post-test to correct for multiple comparisons was performed using Prism 9 software. Error bars indicate SEM. P values (*) < 0.05, (**) < 0.002, (***) < 0.0002, (****) < 0.0001.

Supporting information

S1 Fig. Filamentous actin in senescent 4226 mammary tumor predator cells localizes to prey cells throughout the engulfment process.

Predator 4226 cell line expressing LifeAct-GFP were treated with 750 nM doxorubicin for 24 hours, washed, and plated on untreated mCherry 4226 prey cells. On day 7 post-doxorubicin, cells were fixed, DAPI stained, and imaged on the confocal microscope. (A) Axial slices of a predator cell that has partially overtopped a prey cell. Dashed arrows indicate LifeAct-GFP localization. (B) Axial slices of a predator cell that has partially overtopped a prey cell. Closed arrows indicate LifeAct-GFP localization. (C) Axial slices of a predator cell that has engulfed and broken down a prey cell. Open arrows indicate LifeAct-GFP localization. Scale bar = 10 μm.

https://doi.org/10.1371/journal.pbio.3001858.s001

(PDF)

S2 Fig. PI biosensors expressed in predator 4226 cells reveal PI(4,5)P2, PI(4)P, and PI(3)P localize to the prey cells during engulfment.

4226 mammary tumor cell lines expressing biosensors that detect indicated PI species were treated with 750 nM doxorubicin for 24 hours, washed, and plated on untreated mCherry-4226 “prey” cells. Cultures were imaged over days 3–8. Time course live-cell imaging of senescent 4226 cells that express (A) PLCD1-GFP marking PI(4,5)P2; (B) P4M-SidMx2-GFP marking PI(4)P; (C) 2xFYVE-GFP marking PI(4)P; (D) PLCD1(R40L)-GFP (Negative Ctrl) mutant that does not bind PI species, throughout the entire process of engulfing mCherry-4226 cells. Scale bar = 100 μm.

https://doi.org/10.1371/journal.pbio.3001858.s002

(PDF)

S3 Fig. Model of PI species localization during the 5 stages of engulfment.

During early stages of contact and overtopping, predator PI(4)P and PI(4,5)P2 are highly localized to prey cells and remain concentrated at rounding. PI(3)P is localized at rounding, internalization, and is present during early stages of digestion before dissipating.

https://doi.org/10.1371/journal.pbio.3001858.s003

(PDF)

S4 Fig. PI3C2B but not PIK3C2A, PIK3C3, PIK3C2G is required for engulfment.

(A) In MCF-7 cells, CRISPR-Cas9 knockout of PIK3C2A, PIK3C3, PIK3C2G single-cell clones were screened by immunoblot as indicated in the figure. Clones indicated by asterisk were chosen for further testing. (B) Predator cell engulfment rates (upper) and confluency as a measure of viability (lower) for senescent MCF-7 parental cells and indicated knockout clones were determined by time course imaging in IncuCyte. Underlying data can be found at S1 Data. (C) CRISPR-Cas9 single-cell clone Pik3c2b knockouts of 4226 cells were screened by immunoblot. (D) Predator cell engulfment rates (left) and confluency (right) of senescent 4226 parental cells and 3 Pik3c2b knockout clones were determined by time course imaging. Underlying data can be found at S1 Data. (E, F) Time course live-cell imaging of a senescent MCF-7 cell expressing mCherry (E) or MCF-7-PIK3C2B-KO knockout (F) throughout the entire process of engulfing MCF-7-GFP untreated (NT) cells. Open arrows indicate successful engulfment by a senescent MCF-7 predator cell; closed arrows show representative failed engulfment by MCF-7-PIK3C2B-KO cells. Scale bar = 100 μm.

https://doi.org/10.1371/journal.pbio.3001858.s004

(PDF)

S5 Fig. Inhibition of PI3 kinase does not reduce engulfment rates.

(A-D) MCF-7 (A) and 4226 (B) cells were made senescent by doxorubicin treatment, then treated with the indicated selective PI3K inhibitor and cell viability was determined by MTT assay. Underlying data can be found at S1 Data. (C, D) Mixed cultures of mCherry and GFP expressing cells were made senescent by doxorubicin to generate senescent predator/senescent prey cultures (DOXO-DOXO cultures), and engulfment rates (upper graphs) and confluency as a measure of viability (lower graphs) were determined for MCF-7 (C) and 4226 cells (D) treated with the PI3 kinase inhibitor as indicated in the figure. Underlying data can be found at S1 Data.

https://doi.org/10.1371/journal.pbio.3001858.s005

(PDF)

S6 Fig. Visualization of IncuCyte parameters for quantification of engulfment rates.

MCF-7 predator cells were treated with 250 nM doxorubicin for 24 hours, washed, and plated on untreated mCherry-MCF-7 “prey” cells. Cultures were imaged over days 4–8. Upper displays unmasked predator (GFP) and prey (mCherry) cells during engulfment. Lower shows a representative imaging mask used for quantification of cell-in-cell. Pink masking is used for predator cells, while blue masking is used for prey cells. Overlap of the 2 masks is indicated by white mask in the final frame. Cells are not counted as “engulfed” until prey are in a rounded and fully circumscribed state of overlap. Closed arrows indicate a cell negative for engulfment; open arrows indicate a cell counted as a positive.

https://doi.org/10.1371/journal.pbio.3001858.s006

(PDF)

S7 Fig. Clathrin is required for efficient completion of internalization.

(A) Predator cell engulfment rates (upper) and confluency (lower) were determined for senescent 4226-GFP cells that were treated with indicated drugs. Underlying data can be found at S1 Data. (B) Predator cell engulfment rates (upper) and confluency (lower) were determined for MCF-7 senescent predator/senescent prey cells (“DOXO-DOXO” cultures) that were treated with indicated drugs. Underlying data can be found at S1 Data. (C) Immunoblot verification of shRNA knockdown of CLTC in MCF-7 and MPE600 cells as indicated. Actin loading control is shown in lower panels. (D) Predator cell engulfment rates (upper) and confluency (lower) were determined for senescent GFP-MPE600 cells that were made to express an shRNA construct targeting CLTC. Underlying data can be found at S1 Data.

https://doi.org/10.1371/journal.pbio.3001858.s007

(PDF)

S1 Raw images. The uncropped immunoblots for Figs 6A, S4A, S4C and S7C.

https://doi.org/10.1371/journal.pbio.3001858.s009

(PDF)

S1 Video. Corresponds to Fig 1A.

(A) Time lapse of confocal imaging and volume view reconstruction of doxorubicin-induced senescent MCF-7-LifeAct cell engulfing multiple untreated MCF-7 mCherry cells over the course of 6.75 hours with 11 z-planes taken at a spacing of 1.6 μm per plane. Video was recorded at 4 frames per hour. (B) The same video as in (A), except the top part of the predator cell above the axial plane was removed to visualize cytoplasmic actin.

https://doi.org/10.1371/journal.pbio.3001858.s010

(AVI)

S2 Video. Corresponds to Fig 1B.

(A, B) Time lapse of confocal imaging and volume view reconstructions of 1 doxorubicin-induced senescent MCF-7-LifeAct-GFP cell that has engulfed and is breaking down multiple untreated MCF-7 mCherry cells over the course of 9.75 hours with 17 z-planes taken at a spacing of 1.6 μm per plane, at 4 frames per hour. (A) The same video is shown with the top 51.2% of the axial plane of the cell removed to show interior actin concentration and actin localization upon digestion. (B) Volume view of entire cell.

https://doi.org/10.1371/journal.pbio.3001858.s011

(AVI)

S3 Video. Corresponds to Fig 1C.

Confocal imaging and animation of volume view reconstruction of doxorubicin-induced senescent MCF-7-LifeAct-GFP cell contacting and overtopping an untreated MCF-7 mCherry cell.

https://doi.org/10.1371/journal.pbio.3001858.s012

(AVI)

S4 Video. Corresponds to Fig 1D.

Confocal imaging and animation of volume view reconstruction of doxorubicin-induced senescent MCF-7-LifeAct-GFP cell that has overtopped an untreated MCF-7 mCherry cell.

https://doi.org/10.1371/journal.pbio.3001858.s013

(AVI)

S5 Video. Corresponds to Fig 1E.

Confocal imaging and animation of volume view reconstruction of doxorubicin-induced senescent MCF-7-LifeAct-GFP cell rounding and beginning to internalize an untreated MCF-7 mCherry cell.

https://doi.org/10.1371/journal.pbio.3001858.s014

(AVI)

S6 Video. Corresponds to Fig 3A (lower).

Confocal imaging and animation of volume view reconstruction of doxorubicin-induced senescent MCF-7-PLCD1-GFP cell that has partially engulfed an untreated MCF-7 mCherry cell.

https://doi.org/10.1371/journal.pbio.3001858.s015

(AVI)

S7 Video. Corresponds to Fig 3B (upper).

Confocal imaging and animation of volume view reconstruction of doxorubicin-induced senescent MCF-7-P4M-SidMx2-GFP cell that has partially engulfed an untreated MCF-7 mCherry cell.

https://doi.org/10.1371/journal.pbio.3001858.s016

(AVI)

S8 Video. Corresponds to Fig 3C (upper).

(A) Confocal imaging and animation of volume view reconstruction of doxorubicin-induced senescent MCF-7-2xFYVE-GFP cell that has internalized an untreated MCF-7 mCherry cell. Corresponds to Fig 3C (lower). (B) Zoomed-in volume view reconstruction of the same 2xFYVE-GFP structure.

https://doi.org/10.1371/journal.pbio.3001858.s017

(AVI)

Acknowledgments

The authors wish to acknowledge support form Dorota Wyczechowska for flow cytometry, Osama El Demerdash for guide RNA design, Crystal Tonnessen-Murray for assistance in counting engulfment, and Rodrigo Jacamo for IncuCyte assistance.

References

  1. 1. Bertheau P, Plassa F, Espie M, Turpin E, de Roquancourt A, Marty M, et al. Effect of mutated TP53 on response of advanced breast cancers to high-dose chemotherapy. Lancet. 2002;360(9336):852–4. Epub 2002/09/24. [pii]. pmid:12243922.
  2. 2. Bertheau P, Turpin E, Rickman DS, Espie M, de Reynies A, Feugeas JP, et al. Exquisite sensitivity of TP53 mutant and basal breast cancers to a dose-dense epirubicin-cyclophosphamide regimen. PLoS Med. 2007;4(3):e90. Epub 2007/03/29. 06-PLME-RA-0470R3 [pii]. pmid:17388661; PubMed Central PMCID: PMC1831731.
  3. 3. Ungerleider NA, Rao SG, Shahbandi A, Yee D, Niu T, Frey WD, et al. Breast cancer survival predicted by TP53 mutation status differs markedly depending on treatment. Breast Cancer Res. 2018;20(1):115. pmid:30285883
  4. 4. Jackson JG, Pant V, Li Q, Chang LL, Quintas-Cardama A, Garza D, et al. p53-Mediated Senescence Impairs the Apoptotic Response to Chemotherapy and Clinical Outcome in Breast Cancer. Cancer Cell. 2012;21(6):793–806. Epub 2012/06/16. pmid:22698404; PubMed Central PMCID: PMC3376352.
  5. 5. Tonnessen-Murray C, Ungerleider NA, Rao SG, Wasylishen AR, Frey WD, Jackson JG. p53 Mediates Vast Gene Expression Changes That Contribute to Poor Chemotherapeutic Response in a Mouse Model of Breast Cancer. Transl Oncol. 2018;11(4):930–940. pmid:29852458
  6. 6. Shahbandi A, Rao SG, Anderson AY, Frey WD, Olayiwola JO, Ungerleider NA, et al. BH3 mimetics selectively eliminate chemotherapy-induced senescent cells and improve response in TP53 wild-type breast cancer. Cell Death Differ. 2020;27(11):3097–116. Epub 2020/05/26. PubMed Central PMCID: PMC7560696. pmid:32457483
  7. 7. Saleh T, Carpenter VJ, Tyutyunyk-Massey L, Murray G, Leverson JD, Souers AJ, et al. Clearance of therapy-induced senescent tumor cells by the senolytic ABT-263 via interference with BCL-XL -BAX interaction. Mol Oncol. 2020;14(10):2504–19. Epub 2020/07/12. pmid:32652830; PubMed Central PMCID: PMC7530780.
  8. 8. Maeda M, Scaglia N, Igal RA. Regulation of fatty acid synthesis and Delta9-desaturation in senescence of human fibroblasts. Life Sci. 2009;84(3–4):119–24. Epub 2008/12/09. pmid:19059270.
  9. 9. Kim JS, Kim EJ, Kim HJ, Yang JY, Hwang GS, Kim CW. Proteomic and metabolomic analysis of H2O2-induced premature senescent human mesenchymal stem cells. Exp Gerontol. 2011;46(6):500–10. Epub 2011/03/09. pmid:21382465.
  10. 10. Quijano C, Cao L, Fergusson MM, Romero H, Liu J, Gutkind S, et al. Oncogene-induced senescence results in marked metabolic and bioenergetic alterations. Cell Cycle. 2012;11(7):1383–92. Epub 2012/03/17. pmid:22421146; PubMed Central PMCID: PMC3350879.
  11. 11. Gey C, Seeger K. Metabolic changes during cellular senescence investigated by proton NMR-spectroscopy. Mech Ageing Dev. 2013. Epub 2013/02/19. pmid:23416267.
  12. 12. Dorr JR, Yu Y, Milanovic M, Beuster G, Zasada C, Dabritz JH, et al. Synthetic lethal metabolic targeting of cellular senescence in cancer therapy. Nature. 2013;501(7467):421–425. PubMed pmid:23945590.
  13. 13. Laberge RM, Sun Y, Orjalo AV, Patil CK, Freund A, Zhou L, et al. MTOR regulates the pro-tumorigenic senescence-associated secretory phenotype by promoting IL1A translation. Nat Cell Biol. 2015;17(8):1049–1061. pmid:26147250; PubMed Central PMCID: PMC4691706
  14. 14. Novakova Z, Hubackova S, Kosar M, Janderova-Rossmeislova L, Dobrovolna J, Vasicova P, et al. Cytokine expression and signaling in drug-induced cellular senescence. Oncogene. 2010;29(2):273–84. Epub 2009/10/06. onc2009318 [pii]. pmid:19802007.
  15. 15. Cahu J, Bustany S, Sola B. Senescence-associated secretory phenotype favors the emergence of cancer stem-like cells. Cell Death Dis. 2012;3:e446. Epub 2012/12/21. pmid:23254289; PubMed Central PMCID: PMC3542619.
  16. 16. Canino C, Mori F, Cambria A, Diamantini A, Germoni S, Alessandrini G, et al. SASP mediates chemoresistance and tumor-initiating-activity of mesothelioma cells. Oncogene. 2012;31(26):3148–63. Epub 2011/10/25. pmid:22020330.
  17. 17. Coppe JP, Patil CK, Rodier F, Sun Y, Munoz DP, Goldstein J, et al. Senescence-associated secretory phenotypes reveal cell-nonautonomous functions of oncogenic RAS and the p53 tumor suppressor. PLoS Biol. 2008;6(12):2853–68. Epub 2008/12/05. [pii]. pmid:19053174; PubMed Central PMCID: PMC2592359.
  18. 18. Rao SG, Jackson JG. SASP: Tumor Suppressor or Promoter? Yes! Trends Cancer. 2016;2(11):676–687. pmid:28741506
  19. 19. Velarde MC, Demaria M, Campisi J. Senescent cells and their secretory phenotype as targets for cancer therapy. Interdiscip Top Gerontol. 2013;38:17–27. pmid:23503512; PubMed Central PMCID: PMC4167737
  20. 20. Achuthan S, Santhoshkumar TR, Prabhakar J, Nair SA, Pillai MR. Drug-induced senescence generates chemoresistant stemlike cells with low reactive oxygen species. J Biol Chem. 2011;286(43):37813–29. Epub 2011/09/01. pmid:21878644; PubMed Central PMCID: PMC3199523.
  21. 21. Takamori H, Oades ZG, Hoch OC, Burger M, Schraufstatter IU. Autocrine growth effect of IL-8 and GROalpha on a human pancreatic cancer cell line, Capan-1. Pancreas. 2000;21(1):52–6. Epub 2000/07/06. pmid:10881932.
  22. 22. Dhawan P, Richmond A. Role of CXCL1 in tumorigenesis of melanoma. J Leukoc Biol. 2002;72(1):9–18. Epub 2002/07/09. pmid:12101257; PubMed Central PMCID: PMC2668262.
  23. 23. Yang X, Lu P, Fujii C, Nakamoto Y, Gao JL, Kaneko S, et al. Essential contribution of a chemokine, CCL3, and its receptor, CCR1, to hepatocellular carcinoma progression. Int J Cancer. 2006;118(8):1869–76. Epub 2005/11/15. pmid:16284949.
  24. 24. Begley LA, Kasina S, Mehra R, Adsule S, Admon AJ, Lonigro RJ, et al. CXCL5 promotes prostate cancer progression. Neoplasia. 2008;10(3):244–54. Epub 2008/03/06. pmid:18320069; PubMed Central PMCID: PMC2262133.
  25. 25. Karnoub AE, Weinberg RA. Chemokine networks and breast cancer metastasis. Breast Dis. 2006;26:75–85. Epub 2007/05/03. pmid:17473367.
  26. 26. Xu J, Patel NH, Gewirtz DA. Triangular Relationship between p53, Autophagy, and Chemotherapy Resistance. Int J Mol Sci. 2020;21(23). Epub 2020/12/02. pmid:33256191; PubMed Central PMCID: PMC7730978.
  27. 27. Tonnessen-Murray CA, Frey WD, Rao SG, Shahbandi A, Ungerleider NA, Olayiwola JO, et al. Chemotherapy-induced senescent cancer cells engulf other cells to enhance their survival. J Cell Biol. 2019;218(11):3827–3844. PubMed pmid:31530580.
  28. 28. Overholtzer M, Brugge JS. The cell biology of cell-in-cell structures. Nat Rev Mol Cell Biol. 2008;9(10):796–809. PubMed pmid:18784728.
  29. 29. Flannagan RS, Jaumouille V, Grinstein S. The cell biology of phagocytosis. Annu Rev Pathol. 2012;7:61–98. Epub 2011/09/14. pmid:21910624.
  30. 30. Riedl J, Crevenna AH, Kessenbrock K, Yu JH, Neukirchen D, Bista M, et al. Lifeact: a versatile marker to visualize F-actin. Nat Methods. 2008;5(7):605–7. Epub 2008/06/10. pmid:18536722; PubMed Central PMCID: PMC2814344.
  31. 31. Taunton J, Rowning BA, Coughlin ML, Wu M, Moon RT, Mitchison TJ, et al. Actin-dependent propulsion of endosomes and lysosomes by recruitment of N-WASP. J Cell Biol. 2000;148(3):519–30. Epub 2000/02/09. pmid:10662777; PubMed Central PMCID: PMC2174808.
  32. 32. Poirier MB, Fiorino C, Rajasekar TK, Harrison RE. F-actin flashes on phagosomes mechanically deform contents for efficient digestion in macrophages. J Cell Sci. 2020;133(12). Epub 2020/06/06. pmid:32501286.
  33. 33. Botelho RJ, Teruel M, Dierckman R, Anderson R, Wells A, York JD, et al. Localized biphasic changes in phosphatidylinositol-4,5-bisphosphate at sites of phagocytosis. J Cell Biol. 2000;151(7):1353–68. Epub 2001/01/03. pmid:11134066; PubMed Central PMCID: PMC2150667.
  34. 34. Marshall JG, Booth JW, Stambolic V, Mak T, Balla T, Schreiber AD, et al. Restricted accumulation of phosphatidylinositol 3-kinase products in a plasmalemmal subdomain during Fc gamma receptor-mediated phagocytosis. J Cell Biol. 2001;153(7):1369–80. Epub 2001/06/27. pmid:11425868; PubMed Central PMCID: PMC2150726.
  35. 35. Vieira OV, Botelho RJ, Rameh L, Brachmann SM, Matsuo T, Davidson HW, et al. Distinct roles of class I and class III phosphatidylinositol 3-kinases in phagosome formation and maturation. J Cell Biol. 2001;155(1):19–25. Epub 2001/10/03. pmid:11581283; PubMed Central PMCID: PMC2150784.
  36. 36. Levin R, Grinstein S, Schlam D. Phosphoinositides in phagocytosis and macropinocytosis. Biochim Biophys Acta. 2015;1851(6):805–23. Epub 2014/09/23. pmid:25238964.
  37. 37. Bohdanowicz M, Grinstein S. Role of phospholipids in endocytosis, phagocytosis, and macropinocytosis. Physiol Rev. 2013;93(1):69–106. Epub 2013/01/11. pmid:23303906.
  38. 38. Swanson JA. Phosphoinositides and engulfment. Cell Microbiol. 2014;16(10):1473–83. Epub 2014/07/31. pmid:25073505; PubMed Central PMCID: PMC4341895.
  39. 39. Varnai P, Balla T. Live cell imaging of phosphoinositide dynamics with fluorescent protein domains. Biochim Biophys Acta. 2006;1761(8):957–67. Epub 2006/05/17. pmid:16702024.
  40. 40. Hammond GR, Fischer MJ, Anderson KE, Holdich J, Koteci A, Balla T, et al. PI4P and PI(4,5)P2 are essential but independent lipid determinants of membrane identity. Science. 2012;337(6095):727–30. Epub 20120621. pmid:22722250; PubMed Central PMCID: PMC3646512.
  41. 41. Bohdanowicz M, Cosio G, Backer JM, Grinstein S. Class I and class III phosphoinositide 3-kinases are required for actin polymerization that propels phagosomes. J Cell Biol. 2010;191(5):999–1012. pmid:21115805; PubMed Central PMCID: PMC2995177
  42. 42. Vieira OV, Botelho RJ, Grinstein S. Phagosome maturation: aging gracefully. Biochem J. 2002;366(Pt 3):689–704. Epub 2002/06/14. pmid:12061891; PubMed Central PMCID: PMC1222826.
  43. 43. Ellson CD, Anderson KE, Morgan G, Chilvers ER, Lipp P, Stephens LR, et al. Phosphatidylinositol 3-phosphate is generated in phagosomal membranes. Curr Biol. 2001;11(20):1631–5. Epub 2001/10/26. pmid:11676926.
  44. 44. Yoshioka K. Class II phosphatidylinositol 3-kinase isoforms in vesicular trafficking. Biochem Soc Trans. 2021;49(2):893–901. Epub 2021/03/06. pmid:33666217; PubMed Central PMCID: PMC8106491.
  45. 45. Kim GH, Dayam RM, Prashar A, Terebiznik M, Botelho RJ. PIKfyve inhibition interferes with phagosome and endosome maturation in macrophages. Traffic. 2014;15(10):1143–63. Epub 2014/07/22. pmid:25041080.
  46. 46. Burke JE. Structural Basis for Regulation of Phosphoinositide Kinases and Their Involvement in Human Disease. Mol Cell. 2018;71(5):653–73. Epub 2018/09/08. pmid:30193094.
  47. 47. Katso RM, Pardo OE, Palamidessi A, Franz CM, Marinov M, De Laurentiis A, et al. Phosphoinositide 3-Kinase C2beta regulates cytoskeletal organization and cell migration via Rac-dependent mechanisms. Mol Biol Cell. 2006;17(9):3729–44. Epub 2006/06/16. pmid:16775008; PubMed Central PMCID: PMC1593155.
  48. 48. Domin J, Harper L, Aubyn D, Wheeler M, Florey O, Haskard D, et al. The class II phosphoinositide 3-kinase PI3K-C2beta regulates cell migration by a PtdIns3P dependent mechanism. J Cell Physiol. 2005;205(3):452–62. Epub 2005/08/23. pmid:16113997.
  49. 49. Aung KT, Yoshioka K, Aki S, Ishimaru K, Takuwa N, Takuwa Y. The class II phosphoinositide 3-kinases PI3K-C2alpha and PI3K-C2beta differentially regulate clathrin-dependent pinocytosis in human vascular endothelial cells. J Physiol Sci. 2019;69(2):263–80. Epub 2018/10/31. pmid:30374841.
  50. 50. Russo A, Okur MN, Bosland M, O’Bryan JP. Phosphatidylinositol 3-kinase, class 2 beta (PI3KC2beta) isoform contributes to neuroblastoma tumorigenesis. Cancer Lett. 2015;359(2):262–8. Epub 2015/01/28. pmid:25622909; PubMed Central PMCID: PMC4351744.
  51. 51. Maffucci T, Cooke FT, Foster FM, Traer CJ, Fry MJ, Falasca M. Class II phosphoinositide 3-kinase defines a novel signaling pathway in cell migration. J Cell Biol. 2005;169(5):789–99. Epub 2005/06/02. pmid:15928202; PubMed Central PMCID: PMC2171608.
  52. 52. Chikh A, Ferro R, Abbott JJ, Pineiro R, Buus R, Iezzi M, et al. Class II phosphoinositide 3-kinase C2beta regulates a novel signaling pathway involved in breast cancer progression. Oncotarget. 2016;7(14):18325–45. Epub 2016/03/05. pmid:26934321; PubMed Central PMCID: PMC4951291.
  53. 53. Liu Z, Sun C, Zhang Y, Ji Z, Yang G. Phosphatidylinositol 3-kinase-C2beta inhibits cisplatin-mediated apoptosis via the Akt pathway in oesophageal squamous cell carcinoma. J Int Med Res. 2011;39(4):1319–32. Epub 2011/10/12. pmid:21986133.
  54. 54. Cisse O, Quraishi M, Gulluni F, Guffanti F, Mavrommati I, Suthanthirakumaran M, et al. Downregulation of class II phosphoinositide 3-kinase PI3K-C2beta delays cell division and potentiates the effect of docetaxel on cancer cell growth. J Exp Clin Cancer Res. 2019;38(1):472. Epub 2019/11/23. pmid:31752944; PubMed Central PMCID: PMC6873561.
  55. 55. Schlam D, Bagshaw RD, Freeman SA, Collins RF, Pawson T, Fairn GD, et al. Phosphoinositide 3-kinase enables phagocytosis of large particles by terminating actin assembly through Rac/Cdc42 GTPase-activating proteins. Nat Commun. 2015;6:8623. Epub 2015/10/16. pmid:26465210; PubMed Central PMCID: PMC4634337.
  56. 56. Araki N, Johnson MT, Swanson JA. A role for phosphoinositide 3-kinase in the completion of macropinocytosis and phagocytosis by macrophages. J Cell Biol. 1996;135(5):1249–60. Epub 1996/12/01. pmid:8947549; PubMed Central PMCID: PMC2121091.
  57. 57. Serova M, de Gramont A, Tijeras-Raballand A, Dos Santos C, Riveiro ME, Slimane K, et al. Benchmarking effects of mTOR, PI3K, and dual PI3K/mTOR inhibitors in hepatocellular and renal cell carcinoma models developing resistance to sunitinib and sorafenib. Cancer Chemother Pharmacol. 2013;71(5):1297–307. Epub 2013/03/13. pmid:23479136.
  58. 58. Folkes AJ, Ahmadi K, Alderton WK, Alix S, Baker SJ, Box G, et al. The identification of 2-(1H-indazol-4-yl)-6-(4-methanesulfonyl-piperazin-1-ylmethyl)-4-morpholin-4-yl-t hieno[3,2-d]pyrimidine (GDC-0941) as a potent, selective, orally bioavailable inhibitor of class I PI3 kinase for the treatment of cancer. J Med Chem. 2008;51(18):5522–32. Epub 2008/08/30. pmid:18754654.
  59. 59. Wheeler M, Domin J. The N-terminus of phosphoinositide 3-kinase-C2beta regulates lipid kinase activity and binding to clathrin. J Cell Physiol. 2006;206(3):586–93. Epub 2005/10/14. pmid:16222711.
  60. 60. McMahon HT, Boucrot E. Molecular mechanism and physiological functions of clathrin-mediated endocytosis. Nat Rev Mol Cell Biol. 2011;12(8):517–33. Epub 2011/07/23. pmid:21779028.
  61. 61. Pizarro-Cerda J, Bonazzi M, Cossart P. Clathrin-mediated endocytosis: what works for small, also works for big. Bioessays. 2010;32(6):496–504. Epub 2010/05/21. pmid:20486136.
  62. 62. Mettlen M, Chen PH, Srinivasan S, Danuser G, Schmid SL. Regulation of Clathrin-Mediated Endocytosis. Annu Rev Biochem. 2018;87:871–96. Epub 2018/04/18. pmid:29661000; PubMed Central PMCID: PMC6383209.
  63. 63. Posor Y, Eichhorn-Gruenig M, Puchkov D, Schoneberg J, Ullrich A, Lampe A, et al. Spatiotemporal control of endocytosis by phosphatidylinositol-3,4-bisphosphate. Nature. 2013;499(7457):233–7. Epub 2013/07/05. pmid:23823722.
  64. 64. Dutta D, Donaldson JG. Search for inhibitors of endocytosis: Intended specificity and unintended consequences. Cell Logist. 2012;2(4):203–8. Epub 2013/03/30. pmid:23538558; PubMed Central PMCID: PMC3607622.
  65. 65. Macia E, Ehrlich M, Massol R, Boucrot E, Brunner C, Kirchhausen T. Dynasore, a cell-permeable inhibitor of dynamin. Dev Cell. 2006;10(6):839–50. Epub 2006/06/03. pmid:16740485.
  66. 66. Rothberg KG, Ying YS, Kamen BA, Anderson RG. Cholesterol controls the clustering of the glycophospholipid-anchored membrane receptor for 5-methyltetrahydrofolate. J Cell Biol. 1990;111(6 Pt 2):2931–8. Epub 1990/12/01. pmid:2148564; PubMed Central PMCID: PMC2116385.
  67. 67. Puri V, Watanabe R, Singh RD, Dominguez M, Brown JC, Wheatley CL, et al. Clathrin-dependent and -independent internalization of plasma membrane sphingolipids initiates two Golgi targeting pathways. J Cell Biol. 2001;154(3):535–47. Epub 2001/08/02. pmid:11481344; PubMed Central PMCID: PMC2196434.
  68. 68. Davidson AJ, Wood W. Macrophages Use Distinct Actin Regulators to Switch Engulfment Strategies and Ensure Phagocytic Plasticity In Vivo. Cell Rep. 2020;31(8):107692. Epub 2020/05/28. pmid:32460022; PubMed Central PMCID: PMC7262594.
  69. 69. Yam PT, Theriot JA. Repeated cycles of rapid actin assembly and disassembly on epithelial cell phagosomes. Mol Biol Cell. 2004;15(12):5647–58. Epub 2004/10/01. pmid:15456901; PubMed Central PMCID: PMC532043.
  70. 70. Lancaster CE, Fountain A, Dayam RM, Somerville E, Sheth J, Jacobelli V, et al. Phagosome resolution regenerates lysosomes and maintains the degradative capacity in phagocytes. J Cell Biol. 2021;220(9). Epub 2021/06/29. pmid:34180943; PubMed Central PMCID: PMC8241537.
  71. 71. Rizzo MA, Davidson MW, Piston DW. Fluorescent protein tracking and detection: fluorescent protein structure and color variants. Cold Spring Harb Protoc. 2009(12):pdb top63. Epub 2010/02/13. pmid:20150100.
  72. 72. Chisari M, Saini DK, Cho JH, Kalyanaraman V, Gautam N. G protein subunit dissociation and translocation regulate cellular response to receptor stimulation. PLoS ONE. 2009;4(11):e7797. Epub 2009/11/26. pmid:19936219; PubMed Central PMCID: PMC2777387.
  73. 73. Varnai P, Balla T. Visualization of phosphoinositides that bind pleckstrin homology domains: calcium- and agonist-induced dynamic changes and relationship to myo-[3H]inositol-labeled phosphoinositide pools. J Cell Biol. 1998;143(2):501–10. Epub 1998/10/24. pmid:9786958; PubMed Central PMCID: PMC2132833.
  74. 74. Padilla-Rodriguez M, Parker SS, Adams DG, Westerling T, Puleo JI, Watson AW, et al. The actin cytoskeletal architecture of estrogen receptor positive breast cancer cells suppresses invasion. Nat Commun. 2018;9(1):2980. Epub 2018/08/01. pmid:30061623; PubMed Central PMCID: PMC6065369.