Extrinsic apoptosis relies on TNF-family receptor activation by immune cells or receptor-activating drugs. Here, we monitored cell cycle progression at a resolution of minutes to relate apoptosis kinetics and cell-to-cell heterogeneities in death decisions to cell cycle phases. Interestingly, we found that cells in S phase delay TRAIL receptor-induced death in favour of mitosis, thereby passing on an apoptosis-primed state to their offspring. This translates into two distinct fates, apoptosis execution post mitosis or cell survival from inefficient apoptosis. Transmitotic resistance is linked to Mcl-1 upregulation and its increased accumulation at mitochondria from mid-S phase onwards, which allows cells to pass through mitosis with activated caspase-8, and with cells escaping apoptosis after mitosis sustaining sublethal DNA damage. Antagonizing Mcl-1 suppresses cell cycle-dependent delays in apoptosis, prevents apoptosis-resistant progression through mitosis and averts unwanted survival after apoptosis induction. Cell cycle progression therefore modulates signal transduction during extrinsic apoptosis, with Mcl-1 governing decision making between death, proliferation and survival. Cell cycle progression thus is a crucial process from which cell-to-cell heterogeneities in fates and treatment outcomes emerge in isogenic cell populations during extrinsic apoptosis.

This article has an associated First Person interview with the first author of the paper.

Healthy development and physiological organ functions rely on appropriately balanced cell proliferation and cell death, and consequently dysregulations in this balance manifest in both proliferative and degenerative diseases (Hanahan and Weinberg, 2011; Zhivotovsky and Orrenius, 2010). Proliferation signalling drives progression through the cell cycle and accumulation of cell mass, whereas apoptosis, the major form of programmed cell death, eliminates cells by triggering proteolytic cascades in which caspases act as cell death initiators and executors. Two main apoptosis pathways activate caspases. The extrinsic pathway is triggered by death ligands binding to their respective cell surface receptor(s), whereas the intrinsic, mitochondrial pathway is initiated in response to cytotoxic stress (Taylor et al., 2008). The extrinsic pathway can be triggered by immune cells to induce apoptosis in transformed and potentially harmful cells (Nagata and Tanaka, 2017), and the death ligand TRAIL (also known as TNFSF10) now serves as a template for the development of novel anti-cancer therapeutics and treatment strategies, with second-generation ligands recently having entered clinical trials (Phillips et al., 2021; Von Karstedt et al., 2017).

Following ligand binding, activated death receptors cluster and form the death-inducing signalling complex (DISC), the activation platform for initiator caspase-8. Active caspase-8 in select cell types sufficiently cleave and activate effector caspase-3, resulting in type I apoptosis (Scaffidi et al., 1998). Typically, however, the presence of the inactive caspase-8 homologue cFLIP (also known as CFLAR) as well as the caspase-3 inhibitor XIAP prevent this type of apoptosis. Instead, the majority of cells rely on the amplification of the apoptotic signal via the mitochondrial pathway (i.e. type II signalling; Danish et al., 2018; Jost et al., 2009; Wilson et al., 2009). Here, caspase-8 cleaves the BH3-only protein Bid, which then translocates to the mitochondria and stimulates Bax- and Bak (also known as BAX and BAK1, respectively)-dependent mitochondrial outer membrane permeabilisation (MOMP), a process that is inhibited by anti-apoptotic Bcl-2 family members such as Bcl-2, Bcl-xL (encoded by BCL2L1) and Mcl-1 (Kalkavan and Green, 2018). MOMP results in the release of cytochrome c and Smac (also known as DIABLO) from the mitochondria, driving efficient caspase-3 activation and apoptosis execution (Tait and Green, 2010).

Whereas prolonged cell cycle arrest, for example by cytotoxic drugs such as doxorubicin, CPT-11, bortezomib, simvastatin or the anti-microtubule agent taxol, evidently induces stress responses that sensitize cells to extrinsic apoptosis (Bullenkamp et al., 2014; Ehrhardt et al., 2013; Jin et al., 2002; Kim et al., 2008; Ray et al., 2007), it so far remains unsolved whether extrinsic apoptosis susceptibility is modulated during normal cell cycle progression. This may not be surprising, since combined and parallel monitoring of externally unperturbed cell cycle progression and apoptosis signalling is challenging. For example, bulk analyses of chemically synchronized cell populations provide little temporal resolution and suffer from interfering with normal cell cycle progression (Cooper, 2003). Single-cell time-lapse analyses, instead, so far remain restricted to studies conducted in the presence of apoptosis sensitizers, such as the translation inhibitor cycloheximide (Rehm et al., 2009; Spencer et al., 2009). Given that cell cycle progression heavily relies on concerted protein production and degradation, and since short-lived modulators of extrinsic apoptosis responsiveness, such as cFLIP, are rapidly lost upon translation inhibition (Poukkula et al., 2005), the question of whether cell cycle progression modulates extrinsic apoptosis signal transduction remains largely unaddressed so far.

Here, we therefore combined mathematical modelling and live-cell imaging to quantitatively study extrinsic apoptosis in TRAIL-susceptible cell lines during unperturbed cell cycle progression, with a specific focus on investigating whether and how cell cycle progression influences how and when cells commit to apoptotic cell death. Interestingly, we found that cells progressing towards mitosis delay apoptotic signal transduction, which was attributable to an intrinsic, cell cycle-associated upregulation of Mcl-1 and, to a lesser extent, Bcl-xL. These cells transition through mitosis with activated caspase-8 and frequently escape cell death, but show signs of sublethal DNA damage. Importantly, Mcl-1-targeting BH3-mimetics suffice to prevent unwanted transmitotic escape from extrinsic apoptosis.

TRAIL-treated cells delay apoptosis execution in favour of mitosis

To analyse whether cell cycle progression influences the timing of apoptosis execution in response to TRAIL treatment, we employed a combined approach of mathematical modelling and experimental observation of cell divisions and cell death events. To non-invasively monitor cell cycle phases, we generated NCI-H460 and HCT-116 cells stably expressing the fluorescently tagged degron of geminin (Kuritz et al., 2017; Sakaue-Sawano et al., 2008), which intrinsically labels cells throughout S, G2 and M phases, but not in G1 (Fig. S1A). Treatment with a translationally relevant TRAIL receptor (TRAIL-R) agonist (Fc-scTRAIL) (Hutt et al., 2017) dose-dependently reduced cell viability, accompanied by cleavage of caspase-8 and caspase-3, as well as PARP (Fig. S1B,C). Addition of the pan-caspase inhibitor QVD-OPh averted cell death, demonstrating that cell death proceeded entirely by apoptosis (Fig. S1D). Following TRAIL addition, we microscopically recorded whether cells executed apoptosis prior to cell division or instead underwent mitosis before executing apoptosis. These data were then compared to expected patterns of apoptosis execution obtained from a mathematical cell population model in which progression towards mitosis and the timing of apoptosis were implemented as independent processes (Fig. 1A) (Imig et al., 2020). Although modelling predicted that ∼20% of the NCI-H460 and HCT-116 cell populations would be expected to divide before executing apoptosis, significantly more cells underwent cell division prior to apoptosis in these experiments (Fig. 1B). Closer examination of the imaging data indeed revealed that geminin-positive cells in S/G2 phases frequently first underwent mitosis, with daughter cells dying soon after cytokinesis (Fig. 1C; Fig. S1E and Movie 1). Correspondingly, when both model and experimental analyses were constrained to cells in S/G2/M phases at the time of TRAIL addition, the mismatch in the timing between apoptosis execution and mitosis was even more pronounced (Fig. 1D). The comparison of mathematical predictions and experimental observations of cell death and cell division events therefore suggests that cell cycle progression might delay or decelerate the kinetics of apoptosis signal transduction, and that such delays would be expected to manifest predominantly during the S/G2/M phases of the cell cycle.

Fig. 1.

TRAIL-treated cells delay apoptosis execution in favour of mitosis. (A) Schematic representation of comparing mathematically modelled and experimentally observed cell cycle progression and apoptosis in response to TRAIL receptor stimulation. (B) Expected and observed cell death events in relation to cell division. Cells were treated with Fc-scTRAIL (0.06 nM) and monitored for 24 h. The percentage of apoptosis events before or after division was analysed for mathematically modelled or experimentally observed cell populations. Data are from at least 80–90 cells from three independent model runs or experiments and are shown as mean±s.d. (C) Representative time-lapse images from five experiments of NCI-H460/geminin cells treated with Fc-scTRAIL (0.06 nM). Scale bar: 50 µm. Cells dying prior to division are highlighted by black arrowheads. Red arrowhead tracks a cell dying after division. (D) Expected and observed cell death events of S/G2/M cells in relation to cell division. Analyses were performed as in B. *P<0.05, **P<0.01 (unpaired two-tailed t-tests).

Fig. 1.

TRAIL-treated cells delay apoptosis execution in favour of mitosis. (A) Schematic representation of comparing mathematically modelled and experimentally observed cell cycle progression and apoptosis in response to TRAIL receptor stimulation. (B) Expected and observed cell death events in relation to cell division. Cells were treated with Fc-scTRAIL (0.06 nM) and monitored for 24 h. The percentage of apoptosis events before or after division was analysed for mathematically modelled or experimentally observed cell populations. Data are from at least 80–90 cells from three independent model runs or experiments and are shown as mean±s.d. (C) Representative time-lapse images from five experiments of NCI-H460/geminin cells treated with Fc-scTRAIL (0.06 nM). Scale bar: 50 µm. Cells dying prior to division are highlighted by black arrowheads. Red arrowhead tracks a cell dying after division. (D) Expected and observed cell death events of S/G2/M cells in relation to cell division. Analyses were performed as in B. *P<0.05, **P<0.01 (unpaired two-tailed t-tests).

Commitment to mitosis and delay of apoptosis are decided at mid to end of the S phase

We next experimentally studied whether signal transduction kinetics towards apoptosis execution indeed differed between cells exposed to TRAIL at different times of the cell cycle and, if this is the case, at which point during cell cycle progression the decision is made to delay apoptosis in favour of mitosis. To this end, we monitored untreated cells for at least 20 h, so that cell divisions could be observed, and the age of each individual cell was known at a resolution of minutes at the time of TRAIL addition. In parallel, we tracked the durations of cell cycle phases based on geminin intensity and cell morphology, as well as the times from TRAIL addition to apoptosis execution (Fig. 2A). At the population level, NCI-H460 and HCT-116 cells underwent apoptosis as soon as 1 h following TRAIL addition, with a median tdeath of ∼3 h (2.8 h and 3.7 h, respectively) (Fig. 2B,C, orange). When separating cell populations based on the cell cycle phase in which they were exposed to TRAIL (Fig. 2A), we noted that cells exposed to TRAIL in S/G2/M phases required significantly longer to die than cells exposed to TRAIL in G1 phase (Fig. 2B,C). The delays observed in the S/G2/M populations could be ascribed to cells dying after mitosis (tdeath calculated from f0 and f1 cell generations) (Fig. 2B,C). Similar results were obtained in HeLa and HT1080 cells treated with TRAIL (Fig. 2D,E), in NCI-H460 cells treated with FasL (also known as FASLG) as an alternative death ligand (Fig. S2A), and in two primary cell isolates obtained from melanoma metastases (Fig. 2F).

Fig. 2.

Commitment to mitosis and delay of apoptosis are decided at mid to end of the S phase. (A) Schematic representation of kinetic analyses. Unsynchronized cell populations were monitored for prolonged times and separated into G1 and S/G2/M subpopulations at the time of Fc-scTRAIL addition, based on geminin intensities. The time required to die (tdeath) were recorded for each cell, taking into account that mitosis might precede cell death. Subsequent analyses aimed at characterizing death kinetics and identifying the decision point at which cells preferentially undergo mitosis prior to executing apoptosis. (B,C) tdeath of NCI-H460/geminin (B) or HCT-116/geminin (C) cells following treatment with Fc-scTRAIL (0.06 nM). Cells were monitored for 20 h prior to Fc-scTRAIL addition to ensure the age of each cell was known (time since previous division, mitotic history). Randomly chosen cells in different cell cycle stages were tracked and tdeath of individual cells was plotted according the population separation described in A. Plots show medians, and 2nd and 3rd quartiles from at least 180 cells from at least three independent experiments. ****P<0.0001, *P<0.05 (Mann–Whitney U-test). (D–F) Death times of HeLa (D) or HT1080 cells (E) or primary melanoma metastases cell isolates M51_2 or M54 (F) following treatment with Fc-scTRAIL (0.06 nM in D,F; 0.6 nM for HT1080). Cells were pre-imaged and exposed to Fc-scTRAIL. Randomly chosen cells were tracked and tdeath of individual cells dying before (f0) or after mitosis (f0+f1) was plotted. Plots show medians, 2nd and 3rd quartiles from at least 50 cells from two to three independent experiments (****P<0.0001 Mann–Whitney U-test). (G,H) tdeath times were plotted against cell cycle positions at the time of Fc-scTRAIL exposure. White areas represent G1, green areas mark S/G2/M phases. Dashed green line denotes cell division. Triangles and squares represent individual cells, with red error bars indicating uncertainty estimates (mean±s.d.). n=271 and n=183 cells were analysed, respectively. Cells that progress through mitosis frequently show a significant delay in apoptosis execution (red boxes).

Fig. 2.

Commitment to mitosis and delay of apoptosis are decided at mid to end of the S phase. (A) Schematic representation of kinetic analyses. Unsynchronized cell populations were monitored for prolonged times and separated into G1 and S/G2/M subpopulations at the time of Fc-scTRAIL addition, based on geminin intensities. The time required to die (tdeath) were recorded for each cell, taking into account that mitosis might precede cell death. Subsequent analyses aimed at characterizing death kinetics and identifying the decision point at which cells preferentially undergo mitosis prior to executing apoptosis. (B,C) tdeath of NCI-H460/geminin (B) or HCT-116/geminin (C) cells following treatment with Fc-scTRAIL (0.06 nM). Cells were monitored for 20 h prior to Fc-scTRAIL addition to ensure the age of each cell was known (time since previous division, mitotic history). Randomly chosen cells in different cell cycle stages were tracked and tdeath of individual cells was plotted according the population separation described in A. Plots show medians, and 2nd and 3rd quartiles from at least 180 cells from at least three independent experiments. ****P<0.0001, *P<0.05 (Mann–Whitney U-test). (D–F) Death times of HeLa (D) or HT1080 cells (E) or primary melanoma metastases cell isolates M51_2 or M54 (F) following treatment with Fc-scTRAIL (0.06 nM in D,F; 0.6 nM for HT1080). Cells were pre-imaged and exposed to Fc-scTRAIL. Randomly chosen cells were tracked and tdeath of individual cells dying before (f0) or after mitosis (f0+f1) was plotted. Plots show medians, 2nd and 3rd quartiles from at least 50 cells from two to three independent experiments (****P<0.0001 Mann–Whitney U-test). (G,H) tdeath times were plotted against cell cycle positions at the time of Fc-scTRAIL exposure. White areas represent G1, green areas mark S/G2/M phases. Dashed green line denotes cell division. Triangles and squares represent individual cells, with red error bars indicating uncertainty estimates (mean±s.d.). n=271 and n=183 cells were analysed, respectively. Cells that progress through mitosis frequently show a significant delay in apoptosis execution (red boxes).

We next sought to identify when during the S/G2/M phases cells decide to delay apoptosis in favour of mitosis. Following data normalisation (Imig et al., 2020), we plotted the tdeath of individual cells against their cell cycle position. This analysis revealed that cells notably delaying apoptosis execution following TRAIL exposure had completed approximately mid-to-late S phase (Fig. 2G,H).

Taken together, these results demonstrate that many cells that would have had sufficient time to die prior to mitosis instead undergo cell division and delay extrinsic apoptosis, a behaviour that manifests once cells have entered mid to late S phase.

Caspase-8 remains active throughout mitosis, but efficient caspase-3 activation is suppressed during G2/M phases

We next studied at which point during signal transduction the propagation of extrinsic apoptotic signalling is delayed. Surface amounts of TRAIL-R1 and -2 (TNFRSF10A and TNFRSF10B) continuously increased from G1 throughout S/G2 and M phases in both NCI-H460 and HCT-116 cells, with M phase amounts approximately twice as high as in G1 (Fig. S2B,C), as would be expected for proteins whose steady state amounts correlate with cell growth (Lin and Amir, 2018). TRAIL binding correlated with these changes (Fig. S2D,E). These results therefore exclude that delays in apoptosis signalling emanate from lack of receptors or impaired ligand binding in later phases of the cell cycle. Next, we studied caspase-8 processing and activity in cells exposed to TRAIL at different times of the cell cycle. Since the activation of downstream caspases interferes with measurements of initiator caspase processing and activities (McStay et al., 2008), we used HCT-116 cells deficient for Bax and Bak expression [hereafter denoted HCT-116 (Bax/Bak)−/−] (Fig. S3A), which are unable to undergo MOMP and apoptosis execution, in our initial analyses. HCT-116 (Bax/Bak)−/− cells sorted by geminin fluorescence for G1 or G2/M phases similarly processed pro-caspase-8 into p43/41 and p18 subunits (Fig. 3A). We next studied whether caspase-8 processing results in notable activity during S/G2/M phases. To this end, we introduced a FRET probe carrying the optimal recognition motif (IETD) for caspase-8-dependent proteolysis (Hellwig et al., 2008) into HCT-116 (Bax/Bak)−/− cells or NCI-H460/Bcl-2 cells (which overexpress Bcl-2; Fig. S3B). Time-lapse imaging revealed comparable caspase-8 activities in all stages of the cell cycle, and that cells can enter and pass through mitosis with activated caspase-8 (Fig. 3B; Fig. S3C,D). Corresponding to these observations, cFLIP, the first substrate directly cleaved by caspase-8 within the DISC, is equally cleaved in G1 or G2/M cells (Fig. 3C). Cleavage of Bid, the substrate of caspase-8 that links caspase-8 activity to the mitochondrial apoptosis pathway and which is essential for TRAIL-induced viability loss in NCI-H460 and HCT-116 cells (Fig. 3D), was indistinguishable between cells in G1 or G2/M phases (Fig. 3E).

Fig. 3.

Caspase-8 remains active throughout mitosis, but efficient caspase-3 activation is suppressed during G2/M phases. (A) Cells were treated with Fc-scTRAIL (0.6 nM) as indicated, and G1 and G2/M subpopulations were separated by FACS. Populations were analysed by immunoblotting. Blots are representative of three independent experiments. The level of cleaved caspase-8 fragments in G1 versus G2/M subpopulations was quantified by densitometry, na, not applicable. (B) HCT-116 (Bax/Bak)−/− cells expressing an IETD FRET probe were exposed to Fc-scTRAIL (0.6 nM) and monitored by time-lapse fluorescence microscopy. Traces show FRET disruption, indicative of caspase-8 activity, in individual cells that are transitioning through mitosis. Data show responses of cells from n=3 independent experiments. The relative mean caspase-8 activity in G1 (n=11) or S/G2/M (n=12) cells was calculated by linear regression analysis , P=0.89 (unpaired two-tailed t-test). The box represents the 25–75th percentiles, and the median is indicated. The whiskers show the range. (C) Cells were treated and sorted as in A and analysed by immunoblotting. Blots are representative of two independent experiments. (D) Cells were treated with Fc-scTRAIL [0.06 nM (NCI-H460) or 0.6 nM (HCT-116)] for 24 h, harvested and stained with Annexin V–GFP and PI. Cells were analysed by flow cytometry and viable fractions are shown as mean±s.d. of one representative experiment from three independent repeats. ****P<0.0001 (unpaired two-tailed t-tests). (E) Cells were treated and sorted as in A and analysed by immunoblotting. Blots are representative of three independent experiments. (F) Cells were treated with Fc-scTRAIL [0.06 nM (NCI-H460) or 0.6 nM (HCT-116)]. Cell lysates were analysed by immunoblotting. Blots are representative of two independent experiments. (G) Control and Bid-depleted (shBid) cells were exposed to Fc-scTRAIL [0.06 nM (NCI-H460) or 0.6 nM (HCT-116)]. Cell lysates were analysed by immunoblotting. Blots shown are representative of two independent experiments. (H) Upper panel shows a schematic representation of caspase-3 processing in type II cells in response to TRAIL receptor activation. Subpanels i–iii show caspase-3 cleavage in the indicated cell lines between 0–6 h after treatment with Fc-scTRAIL (0.6 nM for HCT-116; 0.06 nM for NCI-H460). Population distributions of one representative experiment out of three are shown. Peaks representing cells with pro-caspase-3, limited pre-MOMP cleavage of caspase-3 or full post-MOMP caspase-3 processing are highlighted by asterisks. (I) Quantification of pro-caspase-3 and cleaved caspase-3 signals in distinct cell cycle phases. A population mixture model with three subpopulations was fitted to data shown in H.

Fig. 3.

Caspase-8 remains active throughout mitosis, but efficient caspase-3 activation is suppressed during G2/M phases. (A) Cells were treated with Fc-scTRAIL (0.6 nM) as indicated, and G1 and G2/M subpopulations were separated by FACS. Populations were analysed by immunoblotting. Blots are representative of three independent experiments. The level of cleaved caspase-8 fragments in G1 versus G2/M subpopulations was quantified by densitometry, na, not applicable. (B) HCT-116 (Bax/Bak)−/− cells expressing an IETD FRET probe were exposed to Fc-scTRAIL (0.6 nM) and monitored by time-lapse fluorescence microscopy. Traces show FRET disruption, indicative of caspase-8 activity, in individual cells that are transitioning through mitosis. Data show responses of cells from n=3 independent experiments. The relative mean caspase-8 activity in G1 (n=11) or S/G2/M (n=12) cells was calculated by linear regression analysis , P=0.89 (unpaired two-tailed t-test). The box represents the 25–75th percentiles, and the median is indicated. The whiskers show the range. (C) Cells were treated and sorted as in A and analysed by immunoblotting. Blots are representative of two independent experiments. (D) Cells were treated with Fc-scTRAIL [0.06 nM (NCI-H460) or 0.6 nM (HCT-116)] for 24 h, harvested and stained with Annexin V–GFP and PI. Cells were analysed by flow cytometry and viable fractions are shown as mean±s.d. of one representative experiment from three independent repeats. ****P<0.0001 (unpaired two-tailed t-tests). (E) Cells were treated and sorted as in A and analysed by immunoblotting. Blots are representative of three independent experiments. (F) Cells were treated with Fc-scTRAIL [0.06 nM (NCI-H460) or 0.6 nM (HCT-116)]. Cell lysates were analysed by immunoblotting. Blots are representative of two independent experiments. (G) Control and Bid-depleted (shBid) cells were exposed to Fc-scTRAIL [0.06 nM (NCI-H460) or 0.6 nM (HCT-116)]. Cell lysates were analysed by immunoblotting. Blots shown are representative of two independent experiments. (H) Upper panel shows a schematic representation of caspase-3 processing in type II cells in response to TRAIL receptor activation. Subpanels i–iii show caspase-3 cleavage in the indicated cell lines between 0–6 h after treatment with Fc-scTRAIL (0.6 nM for HCT-116; 0.06 nM for NCI-H460). Population distributions of one representative experiment out of three are shown. Peaks representing cells with pro-caspase-3, limited pre-MOMP cleavage of caspase-3 or full post-MOMP caspase-3 processing are highlighted by asterisks. (I) Quantification of pro-caspase-3 and cleaved caspase-3 signals in distinct cell cycle phases. A population mixture model with three subpopulations was fitted to data shown in H.

Caspase-8 additionally cleaves and activates effector pro-caspase-3. In cells relying on the mitochondrial apoptosis pathway, X-linked inhibitor of apoptosis protein (XIAP) suppresses full autoproteolytic maturation of the large caspase-3 subunit (Deveraux et al., 1997; Jost et al., 2009). Consequently, in the absence of MOMP and amplification of apoptosis signalling, caspase-3 processing is limited and only the initial cleavage product of the large caspase-3 subunit can be detected. Indeed, this was seen when comparing caspase-3 processing between HCT-116 (Bax/Bak)−/− cells, parental HCT-116 cells and NCI-H460 cells (Fig. 3F). Similar patterns of caspase-3 processing were observed in Bid-depleted (shBid) cells (Fig. 3G). Unfortunately, parental cells induced to undergo apoptosis were too fragile to be sorted into G1 and G2/M subpopulations, but instead could be fixed and analysed by flow cytometry. Staining HCT-116 (Bax/Bak)−/− cells for DNA and processed caspase-3, revealed identical pro-caspase-3 processing in G1, S and G2/M phases (Fig. 3Hi). In MOMP-competent HCT-116 cells, instead, additional peaks for cells with fully cleaved caspase-3 pools could be detected. Interestingly, these were restricted to cells in G1 and S phases, indicating that cells in G2/M phases suppress full processing of caspase-3 (Fig. 3Hii). NCI-H460 cells similarly limited caspase-3 processing in G2/M phases, indicative of evading the post-MOMP phase of apoptosis execution (Fig. 3Hiii). When quantifying the kinetics of pro-caspase-3 cleavage events in the distinct subpopulations, we noted that initial pro-caspase-3 cleavage occurred independently of the cell cycle in both NCI-H460 and HCT-116 cells (Fig. 3I). Pre- and post-MOMP caspase-3 processing followed very similar kinetics in cells from G1 and S phases, whereas the post-MOMP subpopulation remained very low in cells analysed in G2/M phases (Fig. 3I).

Taken together, we conclude that caspase-8 is activated independently of the cell cycle phase, that caspase-8 remains active throughout mitosis and that activity can be transmitted to daughter cells. The natural caspase-8 substrates cFLIP, Bid and pro-caspase-3 can be processed in all phases of the cell cycle. However, cells analysed in G2/M phases can prevent full caspase-3 processing, suggesting that MOMP is not triggered in these subpopulations.

Cells evading TRAIL-induced apoptosis after transitioning through mitosis suffer sublethal damage

Given that cells are capable of evading MOMP and that daughter cells inherit caspase-8 activity following TRAIL treatment, we studied whether such caspase-8 activities inevitably result in apoptosis execution after mitosis. Even though we observed cells dying after mitosis (Fig. 2), we also noted that substantial numbers of cells remained viable after long treatment times (Fig. S1B). As expected, cells dying after mitosis lost their mitochondrial membrane potential, a surrogate marker for MOMP, just prior to displaying apoptotic phenotypes (Fig. 4A,B). Moreover, we occasionally observed that sibling cells may experience distinct cell fates after mitosis, with mitochondrial depolarization and apoptosis observed in one cell but not in the other (Fig. 4B). It is therefore conceivable that cells might escape apoptosis induction subsequent to caspase-8 activation. We consequently studied whether caspase activities could be observed in cells not showing overt apoptotic phenotypes following TRAIL treatment and mitosis. As a reference, TRAIL treatment caused dramatic caspase-dependent FRET disruption in non-dividing, apoptosis-competent NCI-H460 and HCT116 cells expressing the IETD FRET probe (Fig. 4C). Analysing cells passing through mitosis and subsequently either dying or failing to show an apoptotic phenotype revealed notable differences. Cells dying after mitosis displayed FRET loss that was similar to that in non-dividing cells, whereas cells that escaped apoptosis cleaved significantly less FRET probe (Fig. 4D). Importantly, FRET efficiency in apoptosis-averting cells was significantly different from that in untreated control cells, indicating that caspases were activated in sublethal amounts (Fig. 4D). Apparently, many of these cells suffer DNA damage, as detected by phosphorylation of H2A.X (pH2A.X; also known as γH2AX), a well-established marker for DNA double strand breaks (Fig. 4E,F). Addition of QVD-OPh or siRNA-mediated silencing of caspase-activated DNase (CAD; also known as DFFB), the major executioner of apoptotic DNA fragmentation, suppressed pH2A.X signals, demonstrating that DNA damage was caspase dependent (Fig. 4F,G).

Fig. 4.

Cells evading TRAIL-induced apoptosis after transitioning through mitosis suffer sublethal damage. (A) Geminin-expressing NCI-H460 cells were loaded with TMRM (30 nM) and treated with Fc-scTRAIL (0.06 nM). The white arrowhead indicates a cell progressing through mitosis, followed by mitochondrial depolarization and apoptosis of daughter cells. The yellow arrowhead indicates a cell progressing through mitosis without daughter cells committing apoptotic cell death. Scale bar: 20 µm. (B) Quantification of TMRM intensities in Fc-scTRAIL-treated cells passing through mitosis. Solid and dashed lines indicate two mother cells and their corresponding daughter cells. Asterisks indicate times of apoptotic phenotypes. Blue line indicates a surviving cell. Traces are representative for n=10 mother cells analysed from two independent experiments; au, arbitrary units. (C) NCI-H460 or HCT-116 cells expressing an IETD FRET probe were exposed to Fc-scTRAIL (0.06 nM or 0.6 nM, respectively) alone or in the presence of QVD-OPh (25 µM) and monitored by confocal microscopy. FRET efficiencies were calculated after acceptor photobleaching. Data show mean±s.d. of five cells for each condition. ****P<0.0001; ns, not significant, compared with control (one-way ANOVA followed by Dunnett's multiple comparison). Micrographs show FRET probe expressing HCT-116 cells treated with Fc-scTRAIL (0.6 nM), with the upper panel showing pseudo-coloured CFP/FRET emission ratios. Scale bars: 50 µm. (D) Cells were treated and monitored as in C. FRET efficiencies were analysed in individual cells that progressed through mitosis and for which daughter cells either showed or did not show apoptotic or non-apoptotic phenotypes. Data show mean±s.d. values of control cells (n=3–5), apoptotic cells (n=6) and non-apoptotic cells (n=10–12). ****P<0.0001, **P<0.01, *P<0.05 (one-way ANOVA with Dunnett's multiple comparison test). (E) NCI-H460/geminin cells were treated with Fc-scTRAIL (0.06 nM). After 24 h, adherent surviving cells were stained for pH2A.X (Ser139) and DNA and analysed by confocal microscopy. Images representative of three experiments. Scale bars: 10 µm. (F) Cells were treated as indicated. After 24 h, surviving cells were analysed by western blotting. Tubulin served as loading control. Blots are representative of three independent experiments. (G) Control and CAD-depleted cells were treated with Fc-scTRAIL (0.3 nM) for 24 h. Cell lysates were analysed as indicated. Tubulin was used as loading control. Blots are representative of three independent experiments.

Fig. 4.

Cells evading TRAIL-induced apoptosis after transitioning through mitosis suffer sublethal damage. (A) Geminin-expressing NCI-H460 cells were loaded with TMRM (30 nM) and treated with Fc-scTRAIL (0.06 nM). The white arrowhead indicates a cell progressing through mitosis, followed by mitochondrial depolarization and apoptosis of daughter cells. The yellow arrowhead indicates a cell progressing through mitosis without daughter cells committing apoptotic cell death. Scale bar: 20 µm. (B) Quantification of TMRM intensities in Fc-scTRAIL-treated cells passing through mitosis. Solid and dashed lines indicate two mother cells and their corresponding daughter cells. Asterisks indicate times of apoptotic phenotypes. Blue line indicates a surviving cell. Traces are representative for n=10 mother cells analysed from two independent experiments; au, arbitrary units. (C) NCI-H460 or HCT-116 cells expressing an IETD FRET probe were exposed to Fc-scTRAIL (0.06 nM or 0.6 nM, respectively) alone or in the presence of QVD-OPh (25 µM) and monitored by confocal microscopy. FRET efficiencies were calculated after acceptor photobleaching. Data show mean±s.d. of five cells for each condition. ****P<0.0001; ns, not significant, compared with control (one-way ANOVA followed by Dunnett's multiple comparison). Micrographs show FRET probe expressing HCT-116 cells treated with Fc-scTRAIL (0.6 nM), with the upper panel showing pseudo-coloured CFP/FRET emission ratios. Scale bars: 50 µm. (D) Cells were treated and monitored as in C. FRET efficiencies were analysed in individual cells that progressed through mitosis and for which daughter cells either showed or did not show apoptotic or non-apoptotic phenotypes. Data show mean±s.d. values of control cells (n=3–5), apoptotic cells (n=6) and non-apoptotic cells (n=10–12). ****P<0.0001, **P<0.01, *P<0.05 (one-way ANOVA with Dunnett's multiple comparison test). (E) NCI-H460/geminin cells were treated with Fc-scTRAIL (0.06 nM). After 24 h, adherent surviving cells were stained for pH2A.X (Ser139) and DNA and analysed by confocal microscopy. Images representative of three experiments. Scale bars: 10 µm. (F) Cells were treated as indicated. After 24 h, surviving cells were analysed by western blotting. Tubulin served as loading control. Blots are representative of three independent experiments. (G) Control and CAD-depleted cells were treated with Fc-scTRAIL (0.3 nM) for 24 h. Cell lysates were analysed as indicated. Tubulin was used as loading control. Blots are representative of three independent experiments.

Overall, these experiments therefore indicate that caspase-8 activity, which can persist through mitosis, does not inevitably induce cell death in daughter cells. Cells surviving TRAIL treatment, however, display overt signs of caspase-dependent DNA damage.

Mcl-1 accumulates in cells preparing to enter mitosis

Cells delaying or escaping from TRAIL-induced apoptosis do so by averting MOMP in late stages of the cell cycle. We therefore studied which molecular mechanisms might confer this increased MOMP resistance. The anti-apoptotic Bcl-2 family members Mcl-1 and Bcl-xL, key regulators of MOMP sensitivity, have previously been reported to prevent or suppress apoptosis during extended mitotic arrest induced by anti-mitotic drugs such as paclitaxel and taxol (Bah et al., 2014; Harley et al., 2010; Haschka et al., 2015; Shi et al., 2011; Sloss et al., 2016). First, we therefore determined absolute expression quantities for both proteins in NCI-H460 and HCT-116 cells, in comparison to HeLa cells for which nanomolar expression amounts were reported previously (Lindner et al., 2013). Both cell lines expressed higher levels of Bcl-xL than Mcl-1, and NCI-H460 cells expressed ∼9-fold higher amounts of Mcl-1 than HCT-116 cells (Fig. 5A). The combined expression of Bcl-xL and Mcl-1 was, however, comparable between both cell lines. Studying expression of both proteins across the cell cycle indicated that both cell lines substantially elevated Mcl-1 in M phase (i.e. above the expected doubling throughout the cell cycle; Lin and Amir, 2018; Soltani et al., 2016), as determined by flow cytometry (Fig. 5B). This finding is in line with the Mcl-1 accumulation observed across the cell cycle in synchronised U2OS cell populations (Harley et al., 2010). Higher than expected amounts of Bcl-xL were observed only in HCT-116 cells (Fig. 5B). Immunoblotting independently confirmed the upregulation of Mcl-1 accompanied by increased Noxa (also known as PMAIP1) in later stages of the cell cycle (Fig. 5C). Importantly, TRAIL treatment did not change amounts of Mcl-1 or Bcl-xL (Fig. S4A). A more detailed kinetic analysis of Bcl-xL and Mcl-1 expression across the cell cycle likewise confirmed the increase in Mcl-1, beginning in mid-to-end S phase and accumulating substantially in late G2/M phases (Fig. 5D,E).

Fig. 5.

Mcl-1 accumulates in cells preparing to enter mitosis. (A) Cell lysates were immunoblotted as indicated. Tubulin served as loading control. Blot is representative of three independent experiments. Concentrations of Bcl-xL and Mcl-1 were quantified from signal intensities in comparison to HeLa cells, for which protein concentrations were reported (Lindner et al., 2013). (B) Cells were co-stained with antibodies against Bcl-xL or Mcl-1 and VioBlue-pH3-S10, and analysed by flow cytometry. Median fluorescence intensities were normalized against signals obtained from cells in G1 phase. Data show fold expression (mean±s.e.m.) from at least four independent experiments. Dashed line indicates expected doubling throughout the cell cycle. (C) Cells were FACS-separated into G1 and G2/M populations and analysed by immunoblotting. COX IV or actin served as loading control. Blots are representative of two to three experiments. (D) Cells were stained for Bcl-xL or Mcl-1 and DNA, followed by flow cytometric analysis. Representative heat maps of cell cycle-dependent expression of Bcl-xL (upper panels) and Mcl-1 (lower panels) are shown as a function of DNA content and geminin expression. Three independent experiments were performed. (E) HeLa cells were incubated with EdU (10 µM, 60 min), stained for Bcl-xL or Mcl-1 and DNA and analysed by flow cytometry. Representative heat maps of cell cycle-dependent expression of Bcl-xL (left panel) or Mcl-1 (right panel) are shown as a function of DNA content and EdU intensity. Two independent experiments were performed. (F) Cells were FACS-separated into G1 and G2/M subpopulations, followed mRNA extraction. Transcript levels of BCL2L1 (Bcl-xL), MCL1 (Mcl-1), and CCNB1 (cyclin B1) were analysed by qPCR. Data show means±s.d. from three independent experiments. (G) NCI-H460/geminin cells were treated with DMSO, bortezomib (BTZ, 25 µg/ml) or cycloheximide (CHX, 5 µg/ml) for 1.5 h, co-stained for Mcl-1 and DNA, and analysed by flow cytometry. Mcl-1 amounts are shown as means±s.d. from three independent experiments. *P<0.05; ns, not significant (one-way ANOVA plus Sidak's multiple comparisons test). (H) Cells were analysed for mitochondrial and non-mitochondrial Mcl-1 or Bcl-xL intensities. Shown are median±s.d. calculated from G1 and S/G2 cells from at least two independent experiments ****P<0.0001; **P<0.01; ns, not significant (Welch's unequal variances t-test). au, arbitrary units.

Fig. 5.

Mcl-1 accumulates in cells preparing to enter mitosis. (A) Cell lysates were immunoblotted as indicated. Tubulin served as loading control. Blot is representative of three independent experiments. Concentrations of Bcl-xL and Mcl-1 were quantified from signal intensities in comparison to HeLa cells, for which protein concentrations were reported (Lindner et al., 2013). (B) Cells were co-stained with antibodies against Bcl-xL or Mcl-1 and VioBlue-pH3-S10, and analysed by flow cytometry. Median fluorescence intensities were normalized against signals obtained from cells in G1 phase. Data show fold expression (mean±s.e.m.) from at least four independent experiments. Dashed line indicates expected doubling throughout the cell cycle. (C) Cells were FACS-separated into G1 and G2/M populations and analysed by immunoblotting. COX IV or actin served as loading control. Blots are representative of two to three experiments. (D) Cells were stained for Bcl-xL or Mcl-1 and DNA, followed by flow cytometric analysis. Representative heat maps of cell cycle-dependent expression of Bcl-xL (upper panels) and Mcl-1 (lower panels) are shown as a function of DNA content and geminin expression. Three independent experiments were performed. (E) HeLa cells were incubated with EdU (10 µM, 60 min), stained for Bcl-xL or Mcl-1 and DNA and analysed by flow cytometry. Representative heat maps of cell cycle-dependent expression of Bcl-xL (left panel) or Mcl-1 (right panel) are shown as a function of DNA content and EdU intensity. Two independent experiments were performed. (F) Cells were FACS-separated into G1 and G2/M subpopulations, followed mRNA extraction. Transcript levels of BCL2L1 (Bcl-xL), MCL1 (Mcl-1), and CCNB1 (cyclin B1) were analysed by qPCR. Data show means±s.d. from three independent experiments. (G) NCI-H460/geminin cells were treated with DMSO, bortezomib (BTZ, 25 µg/ml) or cycloheximide (CHX, 5 µg/ml) for 1.5 h, co-stained for Mcl-1 and DNA, and analysed by flow cytometry. Mcl-1 amounts are shown as means±s.d. from three independent experiments. *P<0.05; ns, not significant (one-way ANOVA plus Sidak's multiple comparisons test). (H) Cells were analysed for mitochondrial and non-mitochondrial Mcl-1 or Bcl-xL intensities. Shown are median±s.d. calculated from G1 and S/G2 cells from at least two independent experiments ****P<0.0001; **P<0.01; ns, not significant (Welch's unequal variances t-test). au, arbitrary units.

We next studied which processes might contribute to the accumulation of Mcl-1 in late cell cycle stages. In contrast to cyclin B1 (CCNB1) mRNA, Mcl-1 mRNA did not accumulate in G2/M phases, so that transcriptional upregulation could be excluded (Fig. 5F). After short periods of proteasome inhibition induced by bortezomib (BTZ) treatment, Mcl-1 accumulated in cells from both G1 and S/G2/M phases (Fig. 5G). In reverse, Mcl-1 amounts decreased in both groups of cells when inhibiting protein translation through cycloheximide (CHX) treatment (Fig. 5G). Mcl-1 is therefore translated and degraded in early and in later stages of the cell cycle, so that a combination of alterations in both translation and/or degradation rates would contribute to a net increase in Mcl-1 amounts prior to mitosis. Multi-domain Bcl-2 family members continuously shuttle between mitochondrial membranes and the cytosol (Edlich et al., 2011). Given that proteasomal degradation is limited to the soluble protein fraction, increased Mcl-1 stability could also be achieved by a higher membrane association during later cell cycle stages. Studying subcellular distributions indeed indicated a significant redistribution of Mcl-1 towards the mitochondrial fraction in later cell cycle stages, which was not observed for Bcl-xL (Fig. 5H). Overall, these data suggest that the pronounced accumulation of Mcl-1 towards mitosis and its redistribution towards mitochondrial membranes could contribute to delaying or preventing apoptosis in cells transitioning through mitosis.

Antagonizing Mcl-1 restores apoptosis sensitivity during cell cycle progression

To assess whether anti-apoptotic Bcl-2 family members interfere with apoptosis signal transduction during later stages of the cell cycle, we next studied cells in which Bcl-2, Bcl-xL or Mcl-1 were pharmacologically inhibited. In NCI-H460 cells, Bcl-2 inhibition with ABT-199 or Bcl-xL inhibition by WEHI-539 moderately reduced tdeath times in cells passing through mitosis, whereas Mcl-1 inhibition by S63845 was significantly more effective (Fig. 6A,B compare to Fig. 2G). A similar shortening of tdeath times was observed in HCT-116 (Fig. 6A), HeLa and HT1080 (Fig. S5A,B) cells upon Mcl-1 inhibition. Additionally, treatment with the Bcl-xL inhibitor WEHI-539 also prevented prolonged tdeath times in HCT-116 (Fig. 6A), HeLa or HT1080 (Fig. S5A,B) cells passing through mitosis, possibly due to their high dependency on Bcl-xL expression (Fig. 5A). As expected, Bcl-2 family antagonists also slightly affected Fc-scTRAIL-induced death kinetics in cells that did not progress through mitosis (Fig. S5C,D). Since Mcl-1 inhibition averts delays in apoptosis signalling in cells preparing to enter mitosis, it would be expected that fewer cells overall are able to progress through mitosis under these treatment conditions. Indeed, the frequency of mitotic events observed for TRAIL-treated NCI-H460 and HCT-116 cells was significantly reduced when co-treated with S63845 (Fig. 6C).

Fig. 6.

Antagonizing Mcl-1 restores apoptosis sensitivity during cell cycle progression. (A) Cells were monitored by time-lapse imaging and stimulated with Fc-scTRAIL in presence or absence of Bcl-2 family antagonists (10 µM ABT-199 (Bcl-2), 10 µM WEHI-539 (Bcl-xL), 1 or 10 µM S63845 (Mcl-1). tdeath times were determined as in Fig. 2. Box and whisker plots show tdeath for cells that die after (f0+f1) mitosis. Death times of Fc-scTRAIL-treated cells dying before mitosis (f0, DMSO) are plotted for comparison. Shown are medians with interquartile ranges, plus the minimum to maximum range from n=60–115 cells observed in two to three independent experiments. ****P<0.0001, ***P<0.01, *P<0.05 (Kruskal–Wallis test followed by Dunn's multiple comparison test). (B) Death times were plotted against cell cycle positions at the time of Fc-scTRAIL/S63845 (0.06 nM/1 µM) exposure. White areas represent G1, green areas mark S/G2/M phases. Dashed green line denotes cell division. Triangles and squares represent individual cells, with red error bars indicating uncertainty estimates (mean±s.d.). Mcl-1 inhibition prevents delayed apoptosis execution for cells that progress through mitosis. (C) The percentage of S/G2/M cells that progress through mitosis upon Fc-scTRAIL exposure was determined from cells analysed in A. Data show means±s.d. ***P<0.001; **P<0.01; ns, not significant (one-way ANOVA with Tukey's multiple comparison test). (D) NCI-H460/geminin cells were transfected for 24 h and the amount of Bcl-xL or Mcl-1 was determined by immunoblotting. Tubulin served as loading control. Blot is representative of three independent experiments. (E) Death times after depletion of Mcl-1. NCI-H460/geminin cells were treated with Fc-scTRAIL (0.06 nM). Box-and-whisker plots show time to death (tdeath) for cells that die after mitosis (f0+f1). Shown are medians with interquartile ranges, plus minimum to maximum ranges calculated from at least 75 cells from three independent experiments. ****P<0.0001 (one-way ANOVA with Tukey's multiple comparison test). (F) The percentage of S/G2/M cells undergoing mitosis was calculated from experiments analysed in (D). Data are means±s.d. from n=3 experiments. *P<0.05 (unpaired two-tailed t-test). (G) As in Fig. 3H, panels show caspase-3 cleavage between 0–6 h after treatment with Fc-scTRAIL (0.6 nM for HCT-116 or 0.06 nM for NCI-H460) in presence of Mcl-1 inhibitor A-1201477 (10 µM). Peaks representing cells with pro-caspase-3, limited pre-MOMP cleavage of caspase-3 or full post-MOMP caspase-3 processing are highlighted by asterisks. (H) Quantification of pro-caspase-3 and cleaved caspase-3 signals in distinct cell cycle phases. A population mixture model with three subpopulations was fitted to time course data shown in G. Note the absence of kinetic differences between the subpopulations (compare also to Fig. 3I). (I,J) Cells were treated with Fc-scTRAIL (0.06 nM) and QVD-OPh (25 μM) or the indicated BH3 mimetic (10 μM) for 24 h. Surviving cells were cultivated in fresh medium for another 7 days and stained with Crystal Violet. One of two independent experiments is shown.

Fig. 6.

Antagonizing Mcl-1 restores apoptosis sensitivity during cell cycle progression. (A) Cells were monitored by time-lapse imaging and stimulated with Fc-scTRAIL in presence or absence of Bcl-2 family antagonists (10 µM ABT-199 (Bcl-2), 10 µM WEHI-539 (Bcl-xL), 1 or 10 µM S63845 (Mcl-1). tdeath times were determined as in Fig. 2. Box and whisker plots show tdeath for cells that die after (f0+f1) mitosis. Death times of Fc-scTRAIL-treated cells dying before mitosis (f0, DMSO) are plotted for comparison. Shown are medians with interquartile ranges, plus the minimum to maximum range from n=60–115 cells observed in two to three independent experiments. ****P<0.0001, ***P<0.01, *P<0.05 (Kruskal–Wallis test followed by Dunn's multiple comparison test). (B) Death times were plotted against cell cycle positions at the time of Fc-scTRAIL/S63845 (0.06 nM/1 µM) exposure. White areas represent G1, green areas mark S/G2/M phases. Dashed green line denotes cell division. Triangles and squares represent individual cells, with red error bars indicating uncertainty estimates (mean±s.d.). Mcl-1 inhibition prevents delayed apoptosis execution for cells that progress through mitosis. (C) The percentage of S/G2/M cells that progress through mitosis upon Fc-scTRAIL exposure was determined from cells analysed in A. Data show means±s.d. ***P<0.001; **P<0.01; ns, not significant (one-way ANOVA with Tukey's multiple comparison test). (D) NCI-H460/geminin cells were transfected for 24 h and the amount of Bcl-xL or Mcl-1 was determined by immunoblotting. Tubulin served as loading control. Blot is representative of three independent experiments. (E) Death times after depletion of Mcl-1. NCI-H460/geminin cells were treated with Fc-scTRAIL (0.06 nM). Box-and-whisker plots show time to death (tdeath) for cells that die after mitosis (f0+f1). Shown are medians with interquartile ranges, plus minimum to maximum ranges calculated from at least 75 cells from three independent experiments. ****P<0.0001 (one-way ANOVA with Tukey's multiple comparison test). (F) The percentage of S/G2/M cells undergoing mitosis was calculated from experiments analysed in (D). Data are means±s.d. from n=3 experiments. *P<0.05 (unpaired two-tailed t-test). (G) As in Fig. 3H, panels show caspase-3 cleavage between 0–6 h after treatment with Fc-scTRAIL (0.6 nM for HCT-116 or 0.06 nM for NCI-H460) in presence of Mcl-1 inhibitor A-1201477 (10 µM). Peaks representing cells with pro-caspase-3, limited pre-MOMP cleavage of caspase-3 or full post-MOMP caspase-3 processing are highlighted by asterisks. (H) Quantification of pro-caspase-3 and cleaved caspase-3 signals in distinct cell cycle phases. A population mixture model with three subpopulations was fitted to time course data shown in G. Note the absence of kinetic differences between the subpopulations (compare also to Fig. 3I). (I,J) Cells were treated with Fc-scTRAIL (0.06 nM) and QVD-OPh (25 μM) or the indicated BH3 mimetic (10 μM) for 24 h. Surviving cells were cultivated in fresh medium for another 7 days and stained with Crystal Violet. One of two independent experiments is shown.

Importantly, none of the inhibitors reduced cell viability or induced notable apoptosis in absence of TRAIL treatment (Fig. S5E,F). Using A-1210477 as an alternative Mcl-1 inhibitor gave similar findings (Fig. S5G,H). Similar to what was seen upon pharmacological Mcl-1 inhibition, depleting Mcl-1 expression by siRNA also abolished prolonged tdeath times and reduced the frequency of cell divisions following TRAIL treatment (Fig. 6D–F). In line with these findings, flow cytometric studies revealed that TRAIL-treated cells no longer suppress caspase-3 processing in G2/M phases (Fig. 6G,H, compare to Fig. 3H,I). Instead, the kinetics of caspase-3 processing were largely identical in cells from G1, S and G2/M phases.

Together, these findings demonstrate that Mcl-1 confers transient resistance to TRAIL-induced extrinsic apoptosis in cells entering and passing through mitosis, and that this resistance can be antagonized by inhibiting or depleting Mcl-1. Consequently, the combination of TRAIL and Mcl-1 inhibition almost completely abrogated long-term survival of NCI-H460, HCT-116 (Fig. 6I,J), HeLa or HT1080 (Fig. S5I,J) cells, with HCT-116 cell survival also markedly reduced by co-treatment with Bcl-xL inhibitor (Fig. 6I).

Here, we showed that cells progressing towards and passing through mitosis delay TRAIL-induced extrinsic apoptosis. Importantly, cells transition through mitosis with active caspase-8, and are protected from mitochondrial apoptosis, but nevertheless accumulate caspase-dependent sublethal DNA damage. Upregulation of Mcl-1 and, to a lesser extent, Bcl-xL confers transmitotic apoptosis resistance, with BH3-mimetics targeting these proteins prevent unwanted escape from extrinsic apoptosis.

To obtain insight into the relationship between extrinsic apoptosis signalling and cell cycle progression, we employed fluorescence-based time-lapse imaging of non-synchronized cell populations. This circumvented the limited temporal resolution and accuracy of traditional chemical synchronization (Cooper, 2003). Furthermore, cell cycle arrest for synchronization purposes not only would have interfered with normal cell cycle progression but could have further confounded the analyses by inducing stress responses that modify apoptosis sensitivities in distinct cell cycle phases.

We found that cells delay or prevent apoptosis execution when proceeding towards cell division, with the decision on whether to favour mitosis over apoptosis being made approximately half way through S phase. Decision making therefore substantially precedes entry into mitosis and gives rise to prolonged periods of increased apoptosis resistance, which also explains why cell death events are more frequently observed during G1 phase. The cell-to-cell signalling heterogeneity during the pre-MOMP phase so far was solely attributed to protein expression noise across cell populations or to variation in survival pathway activities (Baskar et al., 2019; Spencer et al., 2009), whereas here we show that cell cycle progression adds an additional and significant non-random layer of cell death regulation. Since we observed similar delays in extrinsic apoptosis signalling when stimulating NCI-H460 cells with FasL, transmitotic apoptosis resistance seems independent of the choice of death ligand.

Mcl-1 is the primary factor conferring transmitotic resistance to extrinsic apoptosis. A key role for Mcl-1 in delaying or suppressing cell death has been reported previously in conditions of mitotic arrest (Allan et al., 2018; Harley et al., 2010; Haschka et al., 2015; Sloss et al., 2016; Wertz et al., 2011), and various E3 ligases have been reported that could contribute to eliminating Mcl-1 (Senichkin et al., 2020). In our setting, where mitosis proceeds unperturbed, cells rapidly deplete Mcl-1 upon entry into G1 phase, with many daughter cells regaining the competency to die due to their inherited caspase-8 activity. Whether one or several of the known Mcl-1 E3 ligases is implicated in rapid Mcl-1 depletion upon regular exit from mitosis remains to be studied.

We found that Mcl-1 steadily accumulates from mid S phase on, matching trends reported for synchronized cell populations (Harley et al., 2010), but also to additionally sharply increase at the very end of the cell cycle. Apoptosis resistance must be very high during this period, since we never observed apoptosis execution during mitosis itself, despite analysing hundreds of cells. Since transcription is largely blocked upon entry into mitosis, Mcl-1 accumulation could arise from continued translation during M phase, possibly combined with suppressed degradation, continuing trends we observed when analysing the pooled fraction of cells in G2/M phases. Indeed, Mcl-1 mRNA continues to be translated during mitosis, as shown in other studies (Sloss et al., 2016; Zhou et al., 2013). Besides translation and proteasomal degradation contributing to defining Mcl-1 amounts throughout the cell cycle, we observed a significant redistribution of Mcl-1 towards mitochondrial membranes in later cell cycle phases. This offers a possibly underappreciated mechanism by which Mcl-1 could be depleted from the degradable soluble fractions and which thereby could contribute to protein stabilisation. Bcl-xL likewise accumulates across the cell cycle but lacks the sharp increase in its amounts during mitosis. While Bcl-xL can reduce apoptosis sensitivity during mitotic arrest by paclitaxel in breast cancer cells (Bah et al., 2014), in our cell line models Mcl-1 inhibition more consistently avoided transmitotic apoptosis resistance. Mcl-1 inhibition likewise outperforms Bcl-xL inhibition in apoptosis induction by mitotic driver drugs, such as aurora kinase inhibitors (Bennett et al., 2016). The third major anti-apoptotic Bcl-2 family member, Bcl-2 itself, is not regulated by cell cycle progression (Harley et al., 2010) and instead is apparently inactivated by phosphorylation prior to/upon entry into mitosis (Terrano et al., 2010). Similar activity-suppressing phosphorylation has also been reported for Bcl-xL (Bah et al., 2014; Basu and Haldar, 2003; Terrano et al., 2010). Overall, this indicates that in particular during mitosis, apoptosis resistance would be expected to strongly rely on Mcl-1. Earlier studies have reported the suppression of caspase-8 activation by Ser-387 phosphorylation in cells arrested in mitosis (Matthess et al., 2010, 2014), whereas we observed that caspase-8 can be activated and remains active in all phases of the cell cycle. The phosphorylation-dependent inhibition of caspase-8 activation therefore might be restricted to cells arrested in mitosis rather than manifesting prominently during regular cell cycle progression or might represent a layer of regulation that does not prominently display in the cellular model systems studied here. Likewise, additional regulatory steps downstream of MOMP, such as the inactivation of caspase-9 by phosphorylation during mitosis (Allan and Clarke, 2007), could further contribute to establishing transiently higher apoptosis resistance.

A more general question arises regarding why cells would need to raise their apoptosis resistance in preparation for and during mitosis. During S phase, cells expend substantial amounts of energy on synthesis and growth processes, which in themselves might impose increased basal stress against which cells protect themselves. Alternatively or in addition, proliferation might be associated with mild or moderate intrinsic pro-apoptotic activities that are not supposed to lead to cell death. Indeed, various apoptotic proteins have been ascribed roles that are independent of cell death. For example, cellular remodelling during development and cell fate determination of stem cells, as well as the regulation of cytoskeletal dynamics partially depend on non-apoptotic caspase signalling (Nakajima and Kuranaga, 2017). Interestingly, RIPK1-containing complexes form during mitosis to activate caspase-8, which contributes to orderly chromosome segregation, without these cells showing signs of apoptosis (Liccardi et al., 2019). Furthermore, caspase-8 activity seems to be required to overcome the G2/M checkpoint by p53 destabilization (Müller et al., 2020). This would imply that apoptosis thresholds must be elevated during mitosis to withstand unwanted apoptosis. In line with this, cells in which Mcl-1 and Bcl-xL are pharmacologically inhibited spontaneously and preferentially die during mitosis (Bennett et al., 2016). Evasion from extrinsic apoptosis thus would arise as an unwanted side effect, which however could have substantial implications for the efficacy by which the immune system can clear transformed cells and, by extension, the risk to relapse from immune stimulatory or death ligand-based therapy (Brown et al., 2018; Seidel et al., 2018; Von Karstedt et al., 2017). Inhibitors of Mcl-1 and Bcl-xL might therefore have a space as co-treatments in such scenarios, especially since Mcl-1 is highly amplified in various cancers (Beroukhim et al., 2010). Results from a clinically relevant TRAIL variant, used in combination with ABT-199 as a Bcl-2 inhibitor, and active clinical studies (Phillips et al., 2021; NCT03082209), could also serve as frameworks in which similar Mcl-1 and Bcl-xL antagonists could be studied. Clinically relevant scenarios might also include treatments with MAPK inhibitors, where Mcl-1 amounts appear to co-determine response variabilities in isogenic cell populations (Inde et al., 2020). Moreover, the Mcl-1 dependency of malignant neoplasms generally might be elevated, since recent CRISPR-Cas9 screening identified Mcl-1 and, with a lower score, Bcl-xL as attractive priority targets in colorectal cancers (Behan et al., 2019).

We noted caspase-dependent DNA damage in cells escaping TRAIL-induced apoptosis by passing through mitosis. DNA damage in TRAIL-surviving cells has been previously described (Lovric and Hawkins, 2010; Miles and Hawkins, 2017), and such damage could arise from moderate activation of caspase-3 either through direct cleavage by caspase-8 and/or the inefficient activation of MOMP, resulting in minority MOMP, a process of inefficient MOMP induction that can contribute to further cell transformation in response to Bcl-2 family antagonism and to inflammatory cytokine production in infection scenarios (Brokatzky et al., 2019; Ichim et al., 2015). Given that caspase-8 can also cleave DEVD sequences, such as found in PARP (Boldin et al., 1996; McStay et al., 2008), it cannot be excluded that caspase-8 itself likewise contributes to DNA damage.

Taken together, our study shows that regular cell cycle progression modulates signal transduction during extrinsic apoptosis, with Mcl-1 modulating cellular decision making between death, proliferation and survival from inefficient apoptosis. These findings not only extend our understanding of the temporal regulation of extrinsic apoptosis susceptibility but also are relevant for better understanding cellular mechanisms contributing to escape from apoptosis induction.

Maintenance of cell lines

Cells were cultured in RPMI 1640 medium supplemented with 10% fetal bovine serum (FBS) and kept at 37°C in 5% CO2 humidified atmosphere. HCT-116 cells were obtained from Interlab Cell Line Collection (Italy). NCI-H460, HeLa and HT1080 cells were obtained from ATCC (LGC Standards GmbH, Germany). Cell isolates from melanoma metastases were kindly provided by Dagmar Kulms (Department of Dermatology, University of Dresden, Germany) and obtained as part of routine resections at University Hospital Dresden, under the auspices of the local Ethics Committee (ethical approval number EK335082018). Platinum-E (Plat-E) cells were kindly provided by Monilola Olayioye, University of Stuttgart, Germany. Cell lines were routinely tested for mycoplasma infection and verified by multiplex human cell line authentication test (Multiplexion GmbH, Germany). NCI-H460/geminin and NCI-H460/Bcl-2 cells were described previously (Danish et al., 2018; Kuritz et al., 2017), HCT-116 (Bax/Bak)−/− cells were kindly provided by Richard J. Youle (National Institutes of Health, USA).

Antibodies and reagents

Primary antibodies were as follows: mouse anti-β-Actin (8H10D10, Cell Signaling #3700, 1:2000), rabbit anti-Bak (Cell Signaling #3814, 1:1000), rabbit anti-Bax (Cell Signaling #2772, 1:1000), rabbit anti-Bcl-xL [Cell Signaling #2764, 1:1000 western blotting (WB), 1:200 flow cytometry (FC) and immunofluoresence (IF)], rabbit anti-Bid (Cell Signaling #2002, 1:1000), rabbit anti-Bim (Cell Signaling #2933, 1:1000), mouse anti-caspase-8 (1C12, Cell Signaling #9746, 1:1000), rabbit anti-caspase-3 (Cell Signaling #9662, 1:1000), rabbit anti-COX IV (3E11, Cell Signaling #4850, 1:1000), rabbit anti-cyclin B1 (D5C10, Cell Signaling #12231, 1:1000 WB, 1:200 FC), rabbit anti-FLIP (D5J1E, Cell Signaling #56343, 1:1000), mouse anti-GAPDH (D4C6R, Cell Signaling #97166, 1:2000), rabbit anti-Mcl-1 (D2W9E, Cell Signaling #94296, 1:1000 WB, 1:50 FC, 1:400 IF), mouse anti-Noxa (F-3, Santa Cruz sc-515840, 1:500), mouse anti-PARP (4C10-5, BD Biosciences 556494, 1:1000), rabbit anti-phospho-histone H2A.X (Ser139) (20E3, Cell Signaling #9718, 1:1000 WB, 1:400 IF), rabbit anti-phospho-histone H3 (Ser10) Pacific Blue™ Conjugate (D2C8, Cell Signaling #8552, 1:50 FC), mouse anti-α-tubulin (DM1A, Cell Signaling #3873, 1:2000), rabbit anti-cleaved caspase-3 AF647 Conjugate (Asp175) (D3E9, Cell Signaling #9602, 1:50 FC), mouse anti-TRAIL-R1 (R&D Systems MAB347, 5 µg/ml FC), mouse anti-TRAIL-R2 (R&D Systems MAB6311, 5 µg/ml FC), mouse IgG1 κ isotype control (BD Biosciences 554121, 5 µg/ml FC), mouse IgG2a κ isotype control (BD Biosciences 554126, 5 µg/ml FC), mouse IgG2b κ isotype control (BD Biosciences 555740, 5 µg/ml FC), mouse anti-FLAG-PE conjugate (BioLegend #637310, 1 µg/ml FC), rabbit (DA1E) mAb IgG XP® isotype control (Cell Signaling #3900, 0.5 µg/ml FC), rabbit IgG isotype control Alexa Fluor 647 conjugate (Cell Signaling #3452, 0.5 µg/ml FC). Secondary antibodies were goat anti-rabbit IgG (H+L) secondary antibody, Alexa Fluor 647 (Thermo Fisher Scientific, #A-21244, 1:500 IF), Alexa Fluor 647-conjugated goat anti-mouse IgG (Dianova #1036-31, 1:500 FC), IRDye®680RD goat anti-rabbit/mouse IgG (LI-COR Biosciences, 1:10,000), IRDye®800CW goat-anti-rabbit/mouse IgG (LI-COR Biosciences, 1:10,000) and peroxidase-conjugated goat anti-rabbit/mouse IgG (Dianova #111-035-144/115-035-062, 1:10,000).

Fc-scTRAIL was produced as described previously (Hutt et al., 2017), FasL-Fc was kindly provided by Harald Wajant (University Hospital Würzburg, Germany) and Annexin V–GFP was produced in-house. DMSO was purchased from Carl Roth, ABT-199 was from Active Biochem, A-1210477 and Q-VD-Oph from Selleckchem, S63845 and WEHI-539 from APExBIO and bortezomib (PS-341) from UBPBio. TMRM, MitoTracker™ Red CMXRos and DAPI were obtained from Invitrogen, propidium iodide (PI) and cycloheximide (CHX) were purchased from Sigma. Puromycin and G-418 were obtained from AppliChem GmbH.

Cell death measurements

Cells were treated as indicated and viability was determined by Crystal Violet staining after 24 h of treatment. Absorbance was measured in a micro plate reader at λ=550 nm and normalized to the value of control cells. Alternatively, cells were harvested, stained with Annexin V-GFP and PI and viability was analysed by flow cytometry (MACSQuant, Miltenyi Biotec, Germany). To analyse the long-term survival, cells were treated as indicated for 24 h. Cells were washed to remove dead cells and surviving cells were cultivated in fresh medium for 7 days. Viability was then determined by Crystal Violet staining.

Plasmids and transfections

Plasmid pEFpuro-hGem(1/110) was obtained by subcloning cDNA encoding mAG-hGem(1/110) (Sakaue-Sawano et al., 2008) into pEF-pGKpuro (kindly provided by Andreas Strasser, WEHI, Australia). The plasmid pSCAT8 encoding the caspase-8 FRET probe CFP-IETD-Venus was described previously (Hellwig et al., 2008), and a SCAT8 variant in a pIRESpuro3 backbone was obtained by subcloning the CFP-IETD-Venus sequence. Cells were transfected using Lipofectamine®LTX (Thermo Fisher Scientific) according to the manufacturer's protocol. At 24 h post transfection, cells were cultured in medium supplemented with appropriate selection antibiotics and fluorescent clones were isolated or enriched after 2 weeks.

The pQCXIPN/ecoR plasmid encoding the receptor for ecotropic murine retroviruses was a kind gift of Monilola Olayioye, University Stuttgart, Germany (Herr et al., 2018). Cells were transfected using Lipofectamine®LTX. At 24 h post transfection, cells were incubated with supernatant from Plat-E cells transfected with pSuper RV or pSUPER RV Bid (#2) (kindly provided by Stephen Tait, University of Glasgow, UK) (Tait et al., 2007). Cells were then selected in medium with puromycin and analyzed for Bid expression.

Silencer®Select Negative Control siRNA #2 (4390846, Thermo Fisher Scientific), Silencer®Select siRNA targeting Mcl-1 (s8583) or Silencer®Select siRNAs targeting CAD (s4060, s223407) were transfected using Lipofectamine®RNAiMAX (Thermo Fisher Scientific) according to the manufacturer's protocol.

Western blotting

Cells were lysed in solubilization buffer [50 mM Tris-HCl pH 7.5, 150 mM NaCl, 1 mM EDTA, 1% (v/v) TritonX-100 plus Complete Protease Inhibitors (Roche) and 1 mM DTT], incubated on ice for 10 min and centrifuged at 16,100 g and 4°C. Protein concentrations were determined by Bradford assay (Bio-Rad). Proteins were separated by SDS-PAGE (4-12% NuPAGE® Novex Bis-Tris gels, Thermo Fisher Scientific) and transferred to nitrocellulose membranes (iBlot® Gel Transfer Stacks, Thermo Fisher Scientific). Membranes were then blocked with 1% blocking reagent (Roche) or 5% BSA in TBS containing 0.05% Tween-20 (TBST) and incubated with primary antibodies followed by POD-conjugated secondary antibodies, and detection was performed using Pierce ECL Western Blotting Substrate (Thermo Fisher Scientific) using an Amersham Imager 600 (GE Healthcare). Alternatively, membranes were incubated with IRDye® secondary antibodies and detection was performed using an Odyssey Imaging system (Li-COR Biosciences).

Quantitative RT-PCR analysis

Total RNA was extracted from cells using the NucleoSpin RNA kit (Machery & Nagel). Synthesis of cDNA and amplification of BCL2L1 (Bcl-xL), MCL1 (Mcl-1), CCNB1 (Cyclin B1) or RPLP0 (60S acidic ribosomal protein P0) was performed using enzyme mix (Power SYBR™ Green RNA-to-CT™ 1-Step kit) and the following primers: Bcl-xL-for, 5′-AGGCGGATTTGAATCTCTTTC-3′; Bcl-xL-rev, 5′-CCCGGTTGCTCTGAGACATT-3′; Mcl-1-for, 5′-AACTGGGGCAGGATTGTGAC-3′; Mcl-1-rev, 5′-CAAACCCATCCCAGCCTCTT-3′; CCNB1-for, 5′-ACCTGTGTCAGGCTTTCTCTG-3′; CCNB1-rev, 5′-CTGACTGCTTGCTCTTCCTCA-3′; RPLP0-for, 5′-CTCTGCATTCTCGCTTCCTGGAG-3′; RPLP0-rev, 5′-CAGATGGATCAGCCAAGAAGG-3′. Relative expression was calculated using the ΔΔCt method and was plotted as fold of G1 mRNA amounts.

Coupled analysis of cell cycle and caspase-3 processing

Cells were harvested, fixed, permeabilised and incubated with Alexa Fluor 647-labelled caspase-3 antibody or isotype control antibody for 60 min, followed by DNA staining (1 µg/ml DAPI) for 10 min. Cells were then washed and analysed by flow cytometry (MACSQuant). Cleaved caspase-3 subpopulations in distinct cell cycle phases were quantified by multi-experiment mixture model analysis based on three lognormal distributions, using the MATLAB Toolbox MEMO (Geissen et al., 2016). Mixture model parameters, comprising subpopulation mean and variance, and subpopulation weights for all time points and cell cycle stages were derived by minimizing the negative logarithm of the likelihood to observe the experimental data.

Live-cell imaging, immunostaining and image analysis

For live-cell imaging, cells were plated on 35 mm glass-bottom dishes (CellView Cell Culture Dish, Greiner Bio One) in Phenol Red-free RPMI 1640 containing 10% FBS. Images were acquired for 20 h in 15 min intervals at 37°C and 5% CO2 on a Zeiss Cell Observer microscope equipped with an Axiocam MRm CCD camera, a Plan-Apochromat 20×/0.8 objective. Fluorescence of mAG-hGeminin(1/110) was acquired with a 470 nm LED module combined with a 62 HE filter set (Carl Zeiss). Then, medium containing Fc-scTRAIL alone or in combination with DMSO or the indicated inhibitors was added and cells were imaged for another 10 h. Randomly chosen cells were tracked and the time of division, onset of geminin expression, subsequent division and tdeath of cells were recorded. Cells were defined as apoptotic upon onset of membrane blebbing.

FRET-based imaging was performed on a Zeiss Cell Observer SD Spinning Disk microscope, equipped with an Axiocam 503 mono CCD camera. For ratiometric FRET imaging, CFP and FRET channels were acquired with a 445 nm diode excitation laser combined with 485/30 nm (CFP) and 562/45 nm (FRET) emission filters. Venus was excited with a 514 nm laser combined with a 562/45 nm emission filter. Images were acquired in 2×2 binning mode at 37°C and 5% CO2 every 15 min for 20 h. Cells were then treated with Fc-scTRAIL and imaged for another 8–10 h. Subsequently, Venus was bleached in the final frame with a 515 nm bleaching laser using a UGA42-Firefly illumination system (Rapp Optoelectronics, Germany). The CFP:YFP emission ratio upon background subtraction over time, the FRET efficiency of the final frame and the FRET probe cleavage were calculated using Zen blue 2.3 software (Zeiss). To record mitochondrial membrane potentials, cells were stained with 30 nM TMRM for 30 min at 37°C and imaged (561 nm excitation laser combined with a 600/50 nm emission filter) in 3 min intervals for 6 h post treatment. Image analysis was performed with Zen Blue 2.3 software (Zeiss). Cell divisions and the time from Fc-scTRAIL addition to apoptotic cell death were analysed in parallel. For pH2A.X staining, cells grown on coverslips were treated for 24 h, dead cells were removed by washing, and surviving cells were fixed, permeabilised and immunostained for pH2A.X followed by DAPI-based nuclear staining. Images were acquired on a Zeiss Axio Observer SD Spinning Disk microscope using 405, 488, and 638 nm excitation with a PlanApochromat 63×/1.4 NA oil objective and an Axiocam 503 Mono CCD camera. Images were processed with the ZEN software (Zeiss).

To analyse Mcl-1 or Bcl-xL subcellular localization, cells grown on coverslips were incubated with MitoTracker Red CMXRos (100 nM) for 30 min before fixation, permeabilisation and immunostaining for Mcl-1 or Bcl-xL. Nuclei were counterstained using DAPI. Images were acquired on a Zeiss Axio Observer SD Spinning Disk microscope equipped with a PlanApochromat 40×/1.4 NA oil objective and an Axiocam 503 Mono CCD camera. Geminin was excited with a 488 nm diode laser using a 525/50 nm emission filter, the Alexa Fluor 647 dye was excited with a 638 nm diode laser using the 690/50 nm emission filter and DAPI was excited with a 405 nm diode laser using the 450/50 nm emission filter. MitoTracker Red was excited with a 561 nm diode laser using a 575/50 nm emission filter. Maximum intensity projections of the acquired z-stacks were calculated and used for analysis with Cell Profiler 3.1.9 (McQuin et al., 2018). In brief, masks for nuclei, mitochondria and cytoplasm (cell mask excluding mitochondria) were segmented, mean fluorescence intensities were measured and the ratio for non-mitochondria to mitochondria was plotted for geminin-negative (G1)/DAPI-low or geminin-positive (late S/G2)/DAPI-high cells.

Mathematical modelling

The mathematical model serving as a reference for independent progression of apoptosis and cell cycle signalling was implemented in MATLAB (The MathWorks, UK). Detailed information on model development and parameterization are reported elsewhere (Imig et al., 2020).

Flow cytometry

For intracellular protein staining, cells were fixed with 4% paraformaldehyde, washed with PBS supplemented with 3% (w/v) FBS, and permeabilised with FACSTM permeabilizing solution 2 (BD Biosciences, Germany). Cells were incubated with primary antibodies diluted in PBS with 3% FBS for 60 min on ice followed by two washing steps. Cells were then incubated with Alexa Fluor 647-conjugated secondary antibody diluted in PBS with 3% FBS. After washing, cells were incubated with phospho-Histone H3 (Ser10) Pacific Blue™ for 30 min and subsequently analysed by flow cytometry.

For DNA labelling using EdU incorporation, cells were incubated in medium containing 10 µM EdU for 60 min and subsequently fixed, permeabilized and stained according the recommended protocol (Click-iT™ EdU AF488 Flow Cytometry assay kit, Thermo Fisher Scientific). Flow cytometry data was processed using custom scripts in MATLAB (The MathWorks, UK).

To analyse the effect of bortezomib (BTZ) or cycloheximide (CHX) on Mcl-1 levels, cells were treated with DMSO, BTZ (25 µg/ml) or CHX (5 µg/ml) for 1.5 h before fixation. For staining of surface TRAIL-R1 and TRAIL-R2, cells were resuspended in PBA (PBS with 0.25% BSA and 0.02% NaN3) containing primary antibody or corresponding isotype control and incubated for 45 min on ice. After washing, cells were incubated with Alexa Fluor 647-conjugated secondary antibody in PBA for 45 min on ice. Cells were fixed, permeabilised and incubated with phospho-histone H3 (Ser10) Pacific Blue™ for 30 min, washed again and analysed by flow cytometry (MACSQuant). To analyse the binding of Fc-scTRAIL, cells were incubated with 0.06 nM Fc-scTRAIL for 15 min, washed and incubated with anti-FLAG–PE conjugate on ice for 30 min. Cells were washed and analysed by flow cytometry.

Preparative cell sorting by FACS

To analyse protein amounts in G1 (geminin-negative) or G2/M (representing the highest 20% of geminin-positive cells) phases, cells were harvested, pelleted and washed with PBS supplemented with 10% FBS. Cells (2×106/ml) were then separated according their mAG-geminin intensity (excitation 494 nm/emission 519 nm) using a FACSAria™III cell sorter (BD Biosciences) at 4°C. Sorted subpopulations were centrifuged and analysed by Western blotting or used for RNA isolation and quantitative RT-PCR analysis.

Statistical analysis

Statistical analysis was performed using PRISM 7.05 (GraphPad Software).

The authors acknowledge Drs. D. Kulms (University of Dresden), M. Olayioye (University of Stuttgart), S. Tait (University of Glasgow), H. Wajant (University of Würzburg) and R. J. Youle (NIH Bethesda) for providing crucial materials. The authors thank Beate Budai and Nathalie Peters for expert technical assistance.

Author contributions

Conceptualization: N.P., P.S., M.R.; Methodology: N.P., S.E., D.I., K.K.; Formal analysis: N.P., D.I., K.K.; Investigation: N.P., A.L., J.S.F., L.D., I.H., J.S., S.E., D.S.; Data curation: N.P.; Writing - original draft: N.P., M.R.; Writing - review & editing: N.P., M.R.; Visualization: N.P., D.I., K.K.; Supervision: P.S., F.A., M.R.; Project administration: N.P., M.R.; Funding acquisition: F.A., P.S., M.R.

Funding Statement

M.R. receives funding from the Deutsche Forschungsgemeinschaft (DFG) under Germany's Excellence Strategy (EXC 2075 – 390740016) and through DFG grants MO 3226/1-1 and MO 3226/4-1. P.S. received funding from the DFG within the project SCHE349/10-1 and the Cluster of Excellence in Simulation Technology (EXC 310/2). N.P. thanks H. Friedrich for financial support.

The peer review history is available online at https://journals.biologists.com/jcs/article-lookup/doi/10.1242/jcs.258966.

Allan
,
L. A.
and
Clarke
,
P. R.
(
2007
).
Phosphorylation of Caspase-9 by CDK1/Cyclin B1 protects mitotic cells against apoptosis
.
Mol. Cell
26
,
301
-
310
.
Allan
,
L. A.
,
Skowyra
,
A.
,
Rogers
,
K. I.
,
Zeller
,
D.
and
Clarke
,
P. R.
(
2018
).
Atypical APC/C–dependent degradation of Mcl–1 provides an apoptotic timer during mitotic arrest
.
EMBO J.
37
,
e96831
.
Bah
,
N.
,
Maillet
,
L.
,
Ryan
,
J.
,
Dubreil
,
S.
,
Gautier
,
F.
,
Letai
,
A.
,
Juin
,
P.
and
Barillé-Nion
,
S.
(
2014
).
Bcl-xL controls a switch between cell death modes during mitotic arrest
.
Cell Death Dis.
5
,
e1291
.
Baskar
,
R.
,
Fienberg
,
H. G.
,
Khair
,
Z.
,
Favaro
,
P.
,
Kimmey
,
S.
,
Green
,
D. R.
,
Nolan
,
G. P.
,
Plevritis
,
S.
and
Bendall
,
S. C.
(
2019
).
TRAIL-induced variation of cell signaling states provides nonheritable resistance to apoptosis
.
Life Sci. Alliance
2
,
e201900554
.
Basu
,
A.
and
Haldar
,
S.
(
2003
).
Identification of a novel Bcl-xL phosphorylation site regulating the sensitivity of taxol- or 2-methoxyestradiol-induced apoptosis
.
FEBS Lett.
538
,
41
-
47
.
Behan
,
F. M.
,
Iorio
,
F.
,
Picco
,
G.
,
Gonçalves
,
E.
,
Beaver
,
C. M.
,
Migliardi
,
G.
,
Santos
,
R.
,
Rao
,
Y.
,
Sassi
,
F.
,
Pinnelli
,
M.
et al. 
(
2019
).
Prioritization of cancer therapeutic targets using CRISPR–Cas9 screens
.
Nature
568
,
511
-
516
.
Bennett
,
A.
,
Sloss
,
O.
,
Topham
,
C.
,
Nelson
,
L.
,
Tighe
,
A.
and
Taylor
,
S. S.
(
2016
).
Inhibition of Bcl-xL sensitizes cells to mitotic blockers, but not mitotic drivers
.
Open Biol.
6
,
160134
.
Beroukhim
,
R.
,
Mermel
,
C. H.
,
Porter
,
D.
,
Wei
,
G.
,
Raychaudhuri
,
S.
,
Donovan
,
J.
,
Barretina
,
J.
,
Boehm
,
J. S.
,
Dobson
,
J.
,
Urashima
,
M.
et al. 
(
2010
).
The landscape of somatic copy-number alteration across human cancers
.
Nature
463
,
899
-
905
.
Boldin
,
M. P.
,
Goncharov
,
T. M.
,
Goltsev
,
Y. V.
and
Wallach
,
D.
(
1996
).
Involvement of MACH, a novel MORT1/FADD-interacting protease, in Fas/APO-1-and TNF receptor-induced cell death
.
Cell
85
,
803
-
815
.
Brokatzky
,
D.
,
Dörflinger
,
B.
,
Haimovici
,
A.
,
Weber
,
A.
,
Kirschnek
,
S.
,
Vier
,
J.
,
Metz
,
A.
,
Henschel
,
J.
,
Steinfeldt
,
T.
,
Gentle
,
I. E.
et al. 
(
2019
).
A non-death function of the mitochondrial apoptosis apparatus in immunity
.
EMBO J.
38
,
e100907
.
Brown
,
J. S.
,
Sundar
,
R.
and
Lopez
,
J.
(
2018
).
Combining DNA damaging therapeutics with immunotherapy: More haste, less speed
.
Br. J. Cancer
118
,
312
-
324
.
Bullenkamp
,
J.
,
Raulf
,
N.
,
Ayaz
,
B.
,
Walczak
,
H.
,
Kulms
,
D.
,
Odell
,
E.
,
Thavaraj
,
S.
and
Tavassoli
,
M.
(
2014
).
Bortezomib sensitises TRAIL-resistant HPV-positive head and neck cancer cells to TRAIL through a caspase-dependent, E6-independent mechanism
.
Cell Death Dis.
5
,
e1489
.
Cooper
,
S.
(
2003
).
Rethinking synchronization of mammalian cells for cell cycle analysis
.
Cell. Mol. Life Sci.
60
,
1099
-
1106
.
Danish
,
L.
,
Imig
,
D.
,
Allgöwer
,
F.
,
Scheurich
,
P.
and
Pollak
,
N.
(
2018
).
Bcl-2-mediated control of TRAIL-induced apoptotic response in the non-small lung cancer cell line NCI-H460 is effective at late caspase processing steps
.
PLoS ONE
13
,
e0198203
.
Deveraux
,
Q. L.
,
Takahashi
,
R.
,
Salvesen
,
G. S.
and
Reed
,
J. C.
(
1997
).
X-linked IAP is a direct inhibitor of cell-death proteases
.
Nature
388
,
300
-
304
.
Edlich
,
F.
,
Banerjee
,
S.
,
Suzuki
,
M.
,
Cleland
,
M. M.
,
Arnoult
,
D.
,
Wang
,
C.
,
Neutzner
,
A.
,
Tjandra
,
N.
and
Youle
,
R. J.
(
2011
).
Bcl-xL retrotranslocates Bax from the mitochondria into the cytosol
.
Cell
145
,
104
-
116
.
Ehrhardt
,
H.
,
Wachter
,
F.
,
Grunert
,
M.
and
Jeremias
,
I.
(
2013
).
Cell cycle-arrested tumor cells exhibit increased sensitivity towards TRAIL-induced apoptosis
.
Cell Death Dis.
4
,
e661
.
Geissen
,
E.-M.
,
Hasenauer
,
J.
,
Heinrich
,
S.
,
Hauf
,
S.
,
Theis
,
F. J.
and
Radde
,
N. E.
(
2016
).
MEMO: Multi-experiment mixture model analysis of censored data
.
Bioinformatics
32
,
2464
-
2472
.
Hanahan
,
D.
and
Weinberg
,
R. A.
(
2011
).
Hallmarks of cancer: the next generation
.
Cell
144
,
646
-
674
.
Harley
,
M. E.
,
Allan
,
L. A.
,
Sanderson
,
H. S.
and
Clarke
,
P. R.
(
2010
).
Phosphorylation of Mcl-1 by CDK1-cyclin B1 initiates its Cdc20-dependent destruction during mitotic arrest
.
EMBO J.
29
,
2407
-
2420
.
Haschka
,
M. D.
,
Soratroi
,
C.
,
Kirschnek
,
S.
,
Häcker
,
G.
,
Hilbe
,
R.
,
Geley
,
S.
,
Villunger
,
A.
and
Fava
,
L. L.
(
2015
).
The NOXA-MCL1-BIM axis defines lifespan on extended mitotic arrest
.
Nat. Commun.
6
,
6891
.
Hellwig
,
C. T.
,
Kohler
,
B. F.
,
Lehtivarjo
,
A.-K.
,
Dussmann
,
H.
,
Courtney
,
M. J.
,
Prehn
,
J. H. M.
and
Rehm
,
M.
(
2008
).
Real time analysis of tumor necrosis factor-related apoptosis-inducing ligand/cycloheximide-induced caspase activities during apoptosis initiation
.
J. Biol. Chem.
283
,
21676
-
21685
.
Herr
,
R.
,
Halbach
,
S.
,
Heizmann
,
M.
,
Busch
,
H.
,
Boerries
,
M.
and
Brummer
,
T.
(
2018
).
BRAF inhibition upregulates a variety of receptor tyrosine kinases and their downstream effector Gab2 in colorectal cancer cell lines
.
Oncogene
37
,
1576
-
1593
.
Hutt
,
M.
,
Marquardt
,
L.
,
Seifert
,
O.
,
Siegemund
,
M.
,
Müller
,
I.
,
Kulms
,
D.
,
Pfizenmaier
,
K.
and
Kontermann
,
R. E.
(
2017
).
Superior properties of Fc-comprising scTRAIL fusion proteins
.
Mol. Cancer Ther.
16
,
2792
-
2802
.
Ichim
,
G.
,
Lopez
,
J.
,
Ahmed
,
S. U.
,
Muthalagu
,
N.
,
Giampazolias
,
E.
,
Delgado
,
M. E.
,
Haller
,
M.
,
Riley
,
J. S.
,
Mason
,
S. M.
,
Athineos
,
D.
et al. 
(
2015
).
Limited mitochondrial permeabilization causes DNA damage and genomic instability in the absence of cell death
.
Mol. Cell
57
,
860
-
872
.
Imig
,
D.
,
Pollak
,
N.
,
Allgöwer
,
F.
and
Rehm
,
M.
(
2020
).
Sample-based modeling reveals bidirectional interplay between cell cycle progression and extrinsic apoptosis
.
PLoS Comput. Biol.
16
,
e1007812
.
Inde
,
Z.
,
Forcina
,
G. C.
,
Denton
,
K.
and
Dixon
,
S. J.
(
2020
).
Kinetic heterogeneity of cancer cell fractional killing
.
Cell Rep.
32
,
107845
.
Jin
,
Z.
,
Dicker
,
D. T.
and
El-Deiry
,
W. S.
(
2002
).
Enhanced sensitivity of G1 arrested human cancer cells suggests a novel therapeutic strategy using a combination of simvastatin and TRAIL
.
Cell Cycle
1
,
82
-
89
.
Jost
,
P. J.
,
Grabow
,
S.
,
Gray
,
D.
,
McKenzie
,
M. D.
,
Nachbur
,
U.
,
Huang
,
D. C. S.
,
Bouillet
,
P.
,
Thomas
,
H. E.
,
Borner
,
C.
,
Silke
,
J.
et al. 
(
2009
).
XIAP discriminates between type i and type II FAS-induced apoptosis
.
Nature
460
,
1035
-
1039
.
Kalkavan
,
H.
and
Green
,
D. R.
(
2018
).
MOMP, cell suicide as a BCL-2 family business
.
Cell Death Differ.
25
,
46
-
55
.
Kim
,
M.
,
Liao
,
J.
,
Dowling
,
M. L.
,
Voong
,
K. R.
,
Parker
,
S. E.
,
Wang
,
S.
,
El-Deiry
,
W. S.
and
Kao
,
G. D.
(
2008
).
TRAIL inactivates the mitotic checkpoint and potentiates death induced by microtubule-targeting agents in human cancer cells
.
Cancer Res.
68
,
3440
-
3449
.
Kuritz
,
K.
,
Stöhr
,
D.
,
Pollak
,
N.
and
Allgöwer
,
F.
(
2017
).
On the relationship between cell cycle analysis with ergodic principles and age-structured cell population models
.
J. Theor. Biol.
414
,
91
-
102
.
Liccardi
,
G.
,
Ramos Garcia
,
L.
,
Tenev
,
T.
,
Annibaldi
,
A.
,
Legrand
,
A. J.
,
Robertson
,
D.
,
Feltham
,
R.
,
Anderton
,
H.
,
Darding
,
M.
,
Peltzer
,
N.
et al. 
(
2019
).
RIPK1 and Caspase-8 ensure chromosome stability independently of their role in cell death and inflammation
.
Mol. Cell
73
,
413
-
428
.
Lin
,
J.
and
Amir
,
A.
(
2018
).
Homeostasis of protein and mRNA concentrations in growing cells
.
Nat. Commun.
9
,
4496
.
Lindner
,
A. U.
,
Concannon
,
C. G.
,
Boukes
,
G. J.
,
Cannon
,
M. D.
,
Llambi
,
F.
,
Ryan
,
D.
,
Boland
,
K.
,
Kehoe
,
J.
,
McNamara
,
D. A.
,
Murray
,
F.
et al. 
(
2013
).
Systems analysis of BCL2 protein family interactions establishes a model to predict responses to chemotherapy
.
Cancer Res.
73
,
519
-
528
.
Lovric
,
M. M.
and
Hawkins
,
C. J.
(
2010
).
TRAIL treatment provokes mutations in surviving cells
.
Oncogene
29
,
5048
-
5060
.
Matthess
,
Y.
,
Raab
,
M.
,
Sanhaji
,
M.
,
Lavrik
,
I. N.
and
Strebhardt
,
K.
(
2010
).
Cdk1/Cyclin B1 controls Fas-mediated apoptosis by regulating caspase-8 activity
.
Mol. Cell. Biol.
30
,
5726
-
5740
.
Matthess
,
Y.
,
Raab
,
M.
,
Knecht
,
R.
,
Becker
,
S.
and
Strebhardt
,
K.
(
2014
).
Sequential Cdk1 and Plk1 phosphorylation of caspase-8 triggers apoptotic cell death during mitosis
.
Mol. Oncol.
8
,
596
-
608
.
McQuin
,
C.
,
Goodman
,
A.
,
Chernyshev
,
V.
,
Kamentsky
,
L.
,
Cimini
,
B. A.
,
Karhohs
,
K. W.
,
Doan
,
M.
,
Ding
,
L.
,
Rafelski
,
S. M.
,
Thirstrup
,
D.
et al. 
(
2018
).
CellProfiler 3.0: next-generation image processing for biology
.
PLoS Biol.
16
,
e2005970
.
McStay
,
G. P.
,
Salvesen
,
G. S.
and
Green
,
D. R.
(
2008
).
Overlapping cleavage motif selectivity of caspases: Implications for analysis of apoptotic pathways
.
Cell Death Differ.
15
,
322
-
331
.
Miles
,
M. A.
and
Hawkins
,
C. J.
(
2017
).
Executioner caspases and CAD are essential for mutagenesis induced by TRAIL or vincristine
.
Cell Death Dis.
8
,
e3062
.
Müller
,
I.
,
Strozyk
,
E.
,
Schindler
,
S.
,
Beissert
,
S.
,
Oo
,
H. Z.
,
Sauter
,
T.
,
Lucarelli
,
P.
,
Raeth
,
S.
,
Hausser
,
A.
,
Al Nakouzi
,
N.
et al. 
(
2020
).
Cancer cells employ nuclear caspase-8 to overcome the p53-dependent G2/M checkpoint through cleavage of USP28
.
Mol. Cell
77
,
970
-
984
.
Nagata
,
S.
and
Tanaka
,
M.
(
2017
).
Programmed cell death and the immune system
.
Nat. Rev. Immunol.
17
,
333
-
340
.
Nakajima
,
Y.
and
Kuranaga
,
E.
(
2017
).
Caspase-dependent non-apoptotic processes in development
.
Cell Death Differ.
24
,
1422
-
1430
.
Phillips
,
D. C.
,
Buchanan
,
F. G.
,
Cheng
,
D.
,
Solomon
,
L. R.
,
Xiao
,
Y.
,
Xue
,
J.
,
Tahir
,
S. K.
,
Smith
,
M. L.
,
Zhang
,
H.
,
Widomski
,
D.
et al. 
(
2021
).
Hexavalent TRAIL fusion protein eftozanermin alfa optimally clusters apoptosis-inducing TRAIL receptors to induce on-target antitumor activity in solid tumors
.
Cancer Res.
81
,
3402
-
3414
.
Poukkula
,
M.
,
Kaunisto
,
A.
,
Hietakangas
,
V.
,
Denessiouk
,
K.
,
Katajamäki
,
T.
,
Johnson
,
M. S.
,
Sistonen
,
L.
and
Eriksson
,
J. E.
(
2005
).
Rapid turnover of c-FLIPshort is determined by its unique C-terminal tail
.
J. Biol. Chem.
280
,
27345
-
27355
.
Ray
,
S.
,
Shyam
,
S.
,
Fraizer
,
G. C.
and
Almasan
,
A.
(
2007
).
S-phase checkpoints regulate Apo2 ligand/TRAIL and CPT-11-induced apoptosis of prostate cancer cells
.
Mol. Cancer Ther.
6
,
1368
-
1378
.
Rehm
,
M.
,
Huber
,
H. J.
,
Hellwig
,
C. T.
,
Anguissola
,
S.
,
Dussmann
,
H.
and
Prehn
,
J. H. M.
(
2009
).
Dynamics of outer mitochondrial membrane permeabilization during apoptosis
.
Cell Death Differ.
16
,
613
-
623
.
Sakaue-Sawano
,
A.
,
Kurokawa
,
H.
,
Morimura
,
T.
,
Hanyu
,
A.
,
Hama
,
H.
,
Osawa
,
H.
,
Kashiwagi
,
S.
,
Fukami
,
K.
,
Miyata
,
T.
,
Miyoshi
,
H.
et al. 
(
2008
).
Visualizing spatiotemporal dynamics of multicellular cell-cycle progression
.
Cell
132
,
487
-
498
.
Scaffidi
,
C.
,
Fulda
,
S.
,
Srinivasan
,
A.
,
Friesen
,
C.
,
Li
,
F.
,
Tomaselli
,
K. J.
,
Debatin
,
K. M.
,
Krammer
,
P. H.
and
Peter
,
M. E.
(
1998
).
Two CD95 (APO-1/Fas) signaling pathways
.
EMBO J.
17
,
1675
-
1687
.
Seidel
,
J. A.
,
Otsuka
,
A.
and
Kabashima
,
K.
(
2018
).
Anti-PD-1 and anti-CTLA-4 therapies in cancer: mechanisms of action, efficacy, and limitations
.
Front. Oncol.
8
,
86
.
Senichkin
,
V. V.
,
Streletskaia
,
A. Y.
,
Gorbunova
,
A. S.
,
Zhivotovsky
,
B.
and
Kopeina
,
G. S.
(
2020
).
Saga of Mcl-1: regulation from transcription to degradation
.
Cell Death Differ.
27
,
405
-
419
.
Shi
,
J.
,
Zhou
,
Y.
,
Huang
,
H. C.
and
Mitchison
,
T. J.
(
2011
).
Navitoclax (ABT-263) accelerates apoptosis during drug-induced mitotic arrest by antagonizing Bcl-xL
.
Cancer Res.
71
,
4518
-
4526
.
Sloss
,
O.
,
Topham
,
C.
and
Taylor
,
S.
(
2016
).
Mcl-1 dynamics influence mitotic slippage and death in mitosis
.
Eur. J. Cancer
7
,
5176
-
5192
.
Soltani
,
M.
,
Vargas-Garcia
,
C. A.
,
Antunes
,
D.
and
Singh
,
A.
(
2016
).
Intercellular variability in protein levels from stochastic expression and noisy Cell Cycle Processes
.
PLoS Comput. Biol.
12
,
e1004972
.
Spencer
,
S. L.
,
Gaudet
,
S.
,
Albeck
,
J. G.
,
Burke
,
J. M.
and
Sorger
,
P. K.
(
2009
).
Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis
.
Nature
459
,
428
-
432
.
Tait
,
S. W. G.
and
Green
,
D. R.
(
2010
).
Mitochondria and cell death: Outer membrane permeabilization and beyond
.
Nat. Rev. Mol. Cell Biol.
11
,
621
-
632
.
Tait
,
S. W. G.
,
De Vries
,
E.
,
Maas
,
C.
,
Keller
,
A. M.
,
D'Santos
,
C. S.
and
Borst
,
J.
(
2007
).
Apoptosis induction by Bid requires unconventional ubiquitination and degradation of its N-terminal fragment
.
J. Cell Biol.
179
,
1453
-
1466
.
Taylor
,
R. C.
,
Cullen
,
S. P.
and
Martin
,
S. J.
(
2008
).
Apoptosis: controlled demolition at the cellular level
.
Nat. Rev. Mol. Cell Biol.
9
,
231
-
241
.
Terrano
,
D. T.
,
Upreti
,
M.
and
Chambers
,
T. C.
(
2010
).
Cyclin-dependent kinase 1-mediated Bcl-xL/Bcl-2 phosphorylation acts as a functional link coupling mitotic arrest and apoptosis
.
Mol. Cell. Biol.
30
,
640
-
656
.
Von Karstedt
,
S.
,
Montinaro
,
A.
and
Walczak
,
H.
(
2017
).
Exploring the TRAILs less travelled: TRAIL in cancer biology and therapy
.
Nat. Rev. Cancer
17
,
352
-
366
.
Wertz
,
I. E.
,
Kusam
,
S.
,
Lam
,
C.
,
Okamoto
,
T.
,
Sandoval
,
W.
,
Anderson
,
D. J.
,
Helgason
,
E.
,
Ernst
,
J. A.
,
Eby
,
M.
,
Liu
,
J.
et al. 
(
2011
).
Sensitivity to antitubulin chemotherapeutics is regulated by MCL1 and FBW7
.
Nature
471
,
110
-
114
.
Wilson
,
T. R.
,
McEwan
,
M.
,
McLaughlin
,
K.
,
Le Clorennec
,
C.
,
Allen
,
W. L.
,
Fennell
,
D. A.
,
Johnston
,
P. G.
and
Longley
,
D. B.
(
2009
).
Combined inhibition of FLIP and XIAP induces Bax-independent apoptosis in type II colorectal cancer cells
.
Oncogene
28
,
63
-
72
.
Zhivotovsky
,
B.
and
Orrenius
,
S.
(
2010
).
Cell cycle and cell death in disease: past, present and future
.
J. Intern. Med.
268
,
395
-
409
.
Zhou
,
T.
,
Li
,
G.
,
Cao
,
B.
,
Liu
,
L.
,
Cheng
,
Q.
,
Kong
,
H.
,
Shan
,
C.
,
Huang
,
X.
,
Chen
,
J.
and
Gao
,
N.
(
2013
).
Downregulation of Mcl-1 through inhibition of translation contributes to benzyl isothiocyanateinduced cell cycle arrest and apoptosis in human leukemia cells
.
Cell Death Dis.
4
,
e515
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information