Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Modelling the dynamics of senescence spread

Modelling the dynamics of senescence spread AbbreviationsABCapproximate bayesian computationEGFepidermal growth factorSASPsenescence‐associated secretory phenotypeINTRODUCTIONSenescence is a cellular stress response that can be activated by telomere shortening, DNA damage and the up‐regulation of oncogenes (d'Adda di Fagagna et al., 2003; Hayflick & Moorhead, 1961; Serrano et al., 1997). Senescent cells display cell cycle arrest and several changes to the cell, including large‐scale chromatin rearrangement and the production of secretory proteins, which can induce senescence in nearby cells (the SASP, senescence associate secretory phenotype) (Acosta et al., 2008; Basisty et al., 2020; Copp'e et al., 2008; Kuilman et al., 2008; Narita et al., 2003). Senescence is necessary for many processes in tissues including development, wound healing and defence against tumour formation (Demaria et al., 2014; Kirschner et al., 2020; Muñoz‐Espín & Serrano, 2014; Passos et al., 2009). When senescent cells accumulate, the molecules they release can create a pro‐tumourigenic and inflammatory microenvironment, which has been linked to both the ageing process and the recurrence of cancerous tumours following treatment (Wang et al., 2020).SASP molecules can spread the senescence phenotype through paracrine signalling, with molecules diffusing from a secreting senescent cell and binding to non‐senescent cells. Several experiments have demonstrated paracrine transmission of senescence, for example, senescent cells contained within a transwell induce senescence in a population of growing cells sharing the same media (Acosta et al., 2013). SASP signalling also recruits immune cells which can clear senescent cells (Prata et al., 2018).The terms “primary” and “secondary” senescence are used to distinguish between cells in which senescence has been induced cell‐intrinsically (e.g., through DNA damage) and extrinsically (through signals from primary senescent cells), respectively (Figure 1b). Recent evidence shows that these subtypes of senescence are transcriptomically distinct (Rattanavirotkul et al., 2020; Teo et al., 2019).1FIGUREModeling the spread of senescence. (a) A controlled senescence spread is beneficial in vivo, but it is unclear how a limited spread of senescence is achieved. (b) Illustration of the cellular phenotypes considered in our model. (c) A cartoon of the simulated diffusion of ligands through a 3D layer of media over a 2D layer of cells. (d) The cumulative probability of a ligand binding at a distance, r, from the cell that emitted it. (e) An illustration of the minimal model of senescence spread developed in this paper, where we consider the probability of inducing senescence in the orange annulus. (f) An illustration of selected model parameters. (g) The ways senescent cells can induce senescence in the cells around them.In addition to paracrine signalling (SASP), senescence can also spread via juxtacrine communication between cells that are in direct contact (NOTCH‐mediated) (Hoare et al., 2016; Teo et al., 2019). Secondary senescent cells created through juxtacrine signalling may differ from secondary senescent cells created through paracrine signalling, with juxtacrine secondary cells producing fewer SASP proteins, alongside other transcriptional differences (Hoare et al., 2016; Rattanavirotkul et al., 2020). Whether these two types of secondary senescent cells play different roles in the spread of senescence is as yet unclear (Kirschner et al., 2020).Here, we define the controlled spread of senescence as a local (finite) senescence spread over an entire human lifespan. In vivo, a limited and controlled senescence spread is thought to play an important role in at least two scenarios. First, any secondary senescent cells created around a primary senescent lesion may help to amplify the signal to the immune system. Second, as a failsafe mechanism, damaged cells that failed to undergo primary senescence will become senescent in response to the primary cells (Figure 1a). However, multiple lines of evidence suggest that an accumulation of senescent cells and uncontrolled spread will ultimately lead to tissue dysfunction (Baker et al., 2011, 2016; Burd et al., 2013; Katzir et al., 2021).Currently, there are no hypothesised mechanisms to contain the spread of senescence in the absence of the immune system. However, it has been shown that the spread of senescence is contained when immune cells are not present (Acosta et al., 2013). In vitro experiments have shown that from a seeded circle of senescent cells, the spread of senescence is finite (senescence was found at a distance of approximately 1 mm from seeded cells) (Acosta et al., 2013). Although the medium taken from primary senescent fibroblasts is capable of inducing secondary senescence in growing fibroblasts, it is uncertain if tertiary senescence is induced when the process is repeated (Acosta et al., 2013). Controlled senescence spread is also observed in vivo, where some senescent lesions persist over time (melanocytic naevi) and are not cleared by the immune system, but the senescence phenotype does not continue to spread (Bartkova et al., 2006). This is not due to SASP production halting, as cells continue to produce a SASP even months after the onset of senescence and enter a deep senescence state, which can lead to chronic senescence in vivo (van Deursen, 2014).In this paper, we use both a minimal mathematical model and stochastic simulation to investigate the mechanisms that determine the dynamics of senescence spread. We identify that juxtacrine‐induced senescent cells can limit senescence spread by acting as a fire break in the wildfire of paracrine signalling and that this effect is more pronounced when SASP‐signalling is delayed compared with juxtacrine signalling. We propose an in vitro experiment and demonstrate that the parameters of senescence spread can be inferred using Bayesian statistics. With the emergence of large‐scale descriptive efforts such as the SenNet consortium (US) (Lee et al., 2022) and the advancement of spatial transcriptomics (Kiss et al., 2022), our model can serve as a framework for further exploration as new data become available.RESULTSSenescence spread modelWe considered senescence spread in a 2D layer of cells covered by a 3D layer of media so that we could compare our results with senescence spread in vitro, as senescence is best characterised in adherent human fibroblasts (Cristofalo et al., 2003; Hayflick & Moorhead, 1961). Despite the 2D nature of our model, we can use it to infer parameters that are unchanged between 2D and 3D models of senescence spread, such as the number of ligands needed to induce senescence in a cell.To simplify the system, we described the cells as circular and non‐motile, evenly and randomly distributed on the 2D plane with a specified density. Following work by Batsilas et al. (2003), we considered the cells to be covered by a layer of media of height h (as shown in Figure 1c). The density of the cells in the 2D plane is given by σ, and the radius of each cell is given by rcell. The number of receptors per cell is Rtot.Paracrine signalling is mediated by ligands, and in our model, the ligands are SASP molecules. Ligands are emitted by senescent cells and diffuse through the surrounding media. DL is the diffusion coefficient for the ligand, and a larger diffusion coefficient results in ligands travelling further in a given time frame. Three rates are used to describe a ligand binding to a receptor; κon gives the rate at which ligands bind to a receptor, κe the rate of endocytosis of bound ligand‐receptor complexes and κoff the rate at which ligands dissociate from receptors.Batsilas et al. describe the cumulative probability that a ligand binds at a distance, r, from the cell that emitted it, as1Pr≈rr+1.1DL/κeff.Here, κeff is the effective binding rate of ligands to cells (defined in Table 2). As κeff depends on several parameters, including the number of receptors per cell, the cumulative binding probability changes with Rtot and other system parameters (Figure 1d).Building on this work, we assume that cells have a fixed number of receptors for detecting ligands, senescent cells secrete ligands at a fixed rate (NE) and that a cell must bind a certain number of ligands (ND) to become senescent. We calculate the rate that ligands bind to a cell at distance r from the emitting cell and the rate that senescence is induced in cells at distance r from a senescent cell.Throughout this paper, we assume that juxtacrine secondary senescent cells do not produce a SASP and cannot, in turn, induce juxtacrine senescence in their neighbours as a blunted SASP response in juxtacrine secondary cells has been reported (Hoare et al., 2016; Teo et al., 2019). We explore different scenarios of senescence spread, such as secondary or tertiary induction of juxtacrine senescence from paracrine senescent cells, as summarised in Figure 1g. We assume that juxtacrine‐mediated senescence induction occurs when cells are in contact, that is, when they are in close spatial proximity.In a minimal model, we calculate the probability of senescence spread from a single senescent cell, by considering the probability of inducing senescence in an annulus of non‐senescent cells around the primary senescent cell, as shown in Figure 1e. By considering the probability of senescence spread to a series of annuli, each with the same width, but an increasing radius, we determine the likelihood that senescence spreads from a single cell.To model more realistic experimental scenarios we built a stochastic simulation of senescence spreading from a collection of cells over time. We use this to investigate the role of secondary senescent cells, by varying parameters such as their SASP secretion, time delays involved in senescence induction and lesion size.Detailed descriptions of the minimal model and stochastic simulations can be found in Methods (Section 4).Model parametersOur model of the spread of senescence contains several parameters, such as the constants relating to ligand–receptor interactions, where exact values have not been determined experimentally. Where this was the case we used values from closely related biological systems, for example, the parameters for EGF and its receptor suggested in Batsilas et al. (2003). These parameters are shown in Table 1.1TABLETable showing the estimated parameters for the spread of EGF, as described in Batsilas et al. (2003), originally from Wiley et al. (2003).ParameterMeaningEstimated valueNENumber of ligands emitted per senescent cell per unit time (h−1)VariableNDNumber of ligands required to induce senescence per unit time (h−1)UnknownDLDiffusion coefficient for ligands (cm2/s)10−6σDensity of cells (m2/m2)0.1–0.4rcellRadius of cells (m)10×10−6RtotTotal number of receptors per cell104–106κeEndocytosis rate (min−1)0.1–0.3κoffDissociation rate (min−1)0.1–0.3κonForward binding rate constant (M−1min−1)108We assume that a non‐senescent cell needs to bind a threshold number of ligands (ND) in a given time frame for senescence to be induced (Figure 1f) and that a single senescent cell produces NE ligands per unit time. We calculated the number of SASP proteins produced by a senescent cell per hour (Figure S7) based on a recent paper by Schafer et al. (2020) that measured the concentrations of SASP proteins produced by radiation‐induced senescent cells. For fibroblasts, the data suggest ~275,000 protein molecules per hour across all SASP proteins investigated. However, if we limit the proteins of interest to Interleukin 6 (IL‐6), Interleukin 8 (IL‐8) and Activin A (secreted proteins thought to be responsible for the induction of senescence; Acosta et al., 2008; Ito et al., 2017; Jochems et al., 2021; Prata et al., 2018; Putavet & de Keizer, 2021), the number produced per hour is much smaller, ~2887.For simplicity, we did not distinguish between the types of ligands produced. Although realistically, a specific ligand will bind to its complementary receptor, and multiple signalling molecules may be required to induce senescence, we did not include this complexity in the model.Cell communication distances are determined by the number of receptors per cell and the effective binding rateTo understand the impact of each model parameter on senescence spread we used our minimal model (described in Section 4.1) to scan through a range of parameters, noting the probability of senescence spread from a single primary senescent cell. We chose to fix the number of ligands emitted per hour by a cell (justification in Section 2.1.1), and vary ND, the number of ligands required to induce senescence. For small ND we observe senescence spread (blue region in Figure 2a) whereas high ND invariably prevents spread (red region in Figure 2a). Parameter regimes for a controlled (finite) spread are likely to be found at the transition between these regions. We, therefore, investigated how changing other parameters shifts the values of ND for which senescence spreads. Parameter values which increase the maximum of ND under which senescence can be induced in a nearby cell correspond to senescence being more likely to spread through the tissue (Figure 2a). Using this method we were able to draw several conclusions about the dynamics of paracrine signalling in a 2D layer of cells.2FIGUREA minimal model for senescence spread. (a) An illustration of the phase space of Rtot and ND, including a region where senescence spreads and one where senescence does not. (b) The division in phase space between senescence spread and no senescence spread for two values of κe and κoff, with and without juxtacrine senescence. The shaded region shows the parameter regime where senescence would spread without juxtacrine senescence, but not with juxtacrine senescence. (c) The range of parameters (shaded) for which senescence only spreads a fraction of the time. (d) The division in phase space between senescence spread and no senescence spread for a range of cell densities. (e) The division in phase space between senescence spread and no senescence spread for two values of Rtot, with and without juxtacrine senescence. The shaded region shows the parameter regime where senescence would spread without juxtacrine senescence, but not with juxtacrine senescence. (f) An illustration of the impact of receptor number on senescence spread. In this figure, unless specified, parameters were, DL=10−6cm2/s, rcell=10μm, κon=108, σ=0.35.We found that there is both an optimum number of receptors (Rtot) and an optimum forward binding rate (κon) for the spread of a signal (Figure 2a–e) in the regime where both κe and κoff (the endocytosis and dissociation rates) are negligible, such that ligands which bind to a receptor remain bound. Intuitively, it may seem that more receptors would increase the signal spread, but this is not the case. When ligands do not dissociate from a cell easily, once they are bound they are removed from the system. Therefore, in regimes with large numbers of receptors, ligands are more likely to bind to the cell that emitted them (known as autocrine signalling) or to bind in spatial proximity, so the signal does not spread. Conversely, in regimes with a small number of receptors, the signal diffuses to a much greater distance from the senescent cell, but it is unlikely that enough ligands will bind to a cell to signal effectively (Shvartsman et al., 2001). This is illustrated in Figure 2f.The same logic applies if κe>>κoff such that the probability or internalisation, ν, tends to 1, meaning that the rate of dissociation is relatively low. If the binding rate is large, many ligands will be sequestered in autocrine trajectories, spreading the signal ineffectively. However, if the binding rate is small it is unlikely that enough ligands will bind to a cell to induce a change. Figure 2e illustrates the effect of increasing the binding rate.When ligands can dissociate effectively and quickly from cells, the optimal number of receptors (for the signal to spread) is increased (Figure 2b; black vs. blue curves). In this regime, a single ligand can bind to multiple cells over the course of its lifetime. Thus, we assume only the last binding event is relevant for senescence induction. In this regime, senescence spreads further with a larger number of receptors and a greater binding rate.A similar effect is observed when increasing the diffusion coefficient (Figure 2c). With increased diffusivity a ligand can more easily move away from the cell that emitted it, there are fewer autocrine binding events, and so a larger number of receptors per cell can be accommodated without limiting the signal spread.Compared with the binding rate and receptor number, the cell density has only a small impact on the system's dynamics, with a greater cell density resulting in a greater chance of senescence spread (Figure 2d).Juxtacrine senescence cells limit the spread of senescenceTo investigate the effect of juxtacrine senescent cells on the spread of senescence, we first considered a ring of juxtacrine cells around the primary senescent cell in our minimal model. We assume juxtacrine secondary cells do not emit a SASP (Figure 1g). As expected, the layer of juxtacrine cells reduces the set of parameters for which senescence spreads (Figure 2b; solid vs. dotted curves), as ligands emitted by the primary senescent cell bind to the juxtacrine senescent cells, which are already senescent, but in this approximation do not emit ligands of their own. Thus, there is a parameter regime where senescence spread would occur without juxtacrine cells but does not when they are present (Figure 2b, shaded region).Using our minimal model we can find parameter regimes that lead to senescence spread from a single cell. However, we can not determine from this model whether the spread of senescence is controlled (finite) or uncontrolled. To further investigate the effect of juxtacrine senescent cells on the controlled spread of senescence over time, we created a stochastic simulation (described in Section 4.2) that also allowed us to investigate the spread from multiple primary senescent cells. Compared with the spread of senescence from a single primary senescent cell, senescence can spread at higher values of ND, and the secondary senescent cells are more densely distributed around the lesion (Figure S1). We found that juxtacrine senescent cells also slow the spread of senescence in this more complex setting (Figure 3), from hours to several days, for a given value of ND (Figure 3c, solid vs. dotted curves). However, once senescence starts to spread, it spreads rapidly and in an uncontrolled manner. As juxtacrine senescent cells in our model do not produce a SASP, their presence is analogous to removing trees to prevent the spread of a forest fire.3FIGUREStochastic simulation of the spread of senescence. (a) The simulation grid randomly populated with cells. (b) An illustration of the time delay in senescence induction. (c) The spread of senescence from a lesion of primary senescent cells with Rtot=105, σ=0.35, ke=koff=0.2, ND=62 and no delay in the induction of senescence. The left‐hand side plots show the spatial distribution of senescent cells and the right‐hand side shows the creation of paracrine secondary senescent cells over time, both with juxtacrine senescence (solid) and without (dotted). Vertical grey dashed lines mark days 1–7. (d) The same simulation as in (c), but with a 6‐day time delay before a cell can induce senescence in its neighbours. (e) The same simulation as in (d), but with immediate induction of juxtacrine senescence and delayed induction of paracrine senescence.A time delay in the senescence induction process limits spread and alters the dynamicsAfter a growing cell binds the required number of ligands to induce senescence, it takes a number of days before the cell becomes fully senescent and it is capable of inducing senescence in its neighbours. To investigate this effect, we introduced a 6‐day time delay between senescence induction and a cell becoming fully senescent in the stochastic simulation, as Young et al. (2009) suggest 6 days are required to establish full Ras‐induced senescence (Figure 3b).We found that introducing a time delay drastically increases the time taken for senescence to spread, as expected (Figure 3d and Figure S1b). An initial region of paracrine secondary senescent cells is created by the primary lesion, limited in diameter by the reach of the SASP from the primary senescent cells (Figure 3d). Once these cells become fully senescent and produce their own SASP, the signal spreads further, inducing a new wave of paracrine secondary senescent cells.The time delay also alters the density of juxtacrine secondary senescent cells, reducing their number at larger distances from the primary lesion (Figure S5c,d). This is due to the assumption that a cell must be fully senescent before it can induce either juxtacrine or paracrine senescence in its neighbours. The delay in juxtacrine induction increases the likelihood that neighbouring cells will become paracrine secondary senescent in the intervening time, increasing the influence of the SASP over senescence spread. The reduction in juxtacrine cells further from the senescent lesion means that as senescence spreads further, the rate of senescence spread increases due to the higher ratio of paracrine secondary senescent cells, which is inconsistent with experimental observations in vitro and in vivo (Acosta et al., 2013; Michaloglou et al., 2005).When ND, the number of ligands that a cell needs to bind to induce senescence, is decreased, the number of paracrine secondary senescent cells plateaus at a larger number for several days, before increasing sharply again once the cells become fully senescent. In an attempt to find a stable system, where the spread of senescence is limited and controlled, we increased ND. Instead of creating a small, stable number of paracrine senescent cells, this led to no senescence spread on a time scale of several years from this small primary lesion, suggesting that varying ND alone cannot lead to controlled senescence spread.A dynamic senescence induction may be required to create a stable senescent lesion. Several experiments have determined that senescence induction in a cell is a dynamic process, with evidence that NOTCH signalling (required for juxtacrine induction) is first up‐regulated before the cell starts to produce cytokines (for paracrine spread) (Hoare et al., 2016). Preliminary evidence suggests that primary cells induce juxtacrine senescence during this intermediate phase (days 3–5) of senescence activation (Teo et al., 2019). This time‐dependent senescence induction may ensure a larger number of juxtacrine secondary senescent cells, limiting the spread of senescence.A dynamic SASP is required to avoid exponential increases in cell numbers at late timesTo further investigate the effects of a dynamic SASP we simulated the limiting case, in which senescent cells can induce juxtacrine senescence in their neighbours as soon as they begin the journey towards becoming fully senescent. In reality, evidence suggests that juxtacrine induction may actually be possible after a 2‐day delay (JAG1, the NOTCH ligand responsible for juxtacrine senescence induction, is upregulated even after 1 or 2 days [Hoare et al., 2016]), but simulating the limiting case of immediate juxtacrine senescence induction clearly shows the effect of this change on the dynamics of senescence spread.The dynamic induction of senescence both slows the spread of senescence and stops the explosion in the number of senescent cells as the lesion grows, which was observed without it (Figure 3d vs. Figure 3e).Reduced SASP from secondary senescent cells can lead to controlled senescence spreadWe wanted to explore parameter regimes in which the spread of senescence is controlled (local over a human lifetime), but this was difficult with such a large number of unknown system parameters.We have previously outlined how senescence spread can be slowed with the presence of juxtacrine senescence, a time delay, and juxtacrine signalling, which precedes paracrine signalling. If the spread of senescence can be slowed drastically enough, this will lead to controlled senescence spread, as even over a human lifespan the senescence spread will be local and contained.Once senescence spread from a lesion has begun it can only stop if the balance of system parameters is right, for example, if the flux of SASP from the growing lesion is no longer enough to induce senescence. One way to achieve this is to ensure that the number of juxtacrine secondary senescent cells is large enough to dilute the SASP. This is also achieved if paracrine secondary senescent cells produce fewer SASP molecules than primary senescent cells, as suggested by experimental evidence (Acosta et al., 2013). In our simulations, we find that in this case, the primary cells in the initial lesion drive most of the senescence spread (Figure S2a). Some secondary senescent cells are generated over short time scales, but the spread eventually slows down to the point where it remains controlled, in some cases even after 500,000 hours (on the order of a human lifespan). When the lesion size is increased the size of the region of secondary cells will also increase (Figure S2b).Parameters of paracrine senescence spread can be determined through inferenceTo estimate the parameters of the model, we propose a similar experiment to those described in Acosta et al. (2013), seeding a circle of primary senescent cells and watching the phenotype spread (Section 4.3). By focusing on the first wave of senescence induction, we can forgo computationally intensive simulations. Due to the time delay before a cell becomes fully senescent, there is a period of time when a cell is unable to pass the senescent phenotype on to neighbouring cells through paracrine signalling. If it takes cells 6 days to become fully senescent, we can stop the experiment before 6 days. Using a marker to detect cells that are becoming senescent, we can determine the region to which senescence has spread through paracrine signalling (Figure 4a).4FIGURE(a) A cartoon showing the experimental setup proposed to infer the system parameters. (b) Posterior distribution of parameters from the Bayesian inference on simulated data. The plots below the diagonal show the kernel density estimation, and the plots above the diagonal show the final points used in the sampling. Red lines show the true parameter values.Focusing on short timescales allows us to build on the minimal model described in Section 4.1 to predict the fraction of senescent cells at a given distance from the primary lesion. By comparing this with the experimentally observed fraction of senescent cells we can determine NE, ND, Rtot and kon using Bayesian inference, fixing other model parameters to experimentally determined values (Section 4.3). As a proof of principle for this experiment, we used our model to generate synthetic data and attempted to infer the parameters used to generate the data. As a prior distribution of the parameters, we chose a uniform distribution over the parameter regime given in Table 1 (Figure 4b).We show that it is possible to infer some of the parameters accurately from such an experiment (Figure 4b). The posterior distribution for Rtot is narrow and centred on the true value, and although the posterior for ND is not as narrow it is well centred on the true value. This reflects the intuition we gained from the minimal model, that senescence spread is sensitive to ND and Rtot. The posterior distributions for NE and κon are much broader and less centred on the true values.DISCUSSIONCurrently, it is unclear how senescence can spread in a controlled way from a primary lesion. We systematically investigated senescence spread, using mathematical modelling informed by experimental evidence. We found that the number of receptors per cell and the binding rate of ligands were critical in determining the distance over which senescence spreads. We observed that an intermediate number of receptors led to the most effective signal spread and uncontrolled senescence spread. Several scenarios (such as juxtacrine‐induced cells, in turn, being capable of juxtacrine senescence induction) were not explored as they would lead to the uncontrolled spread of senescence in any parameter regime, highlighting that the current way of thinking about senescence spread is naive.Our results provide new insights into the role of juxtacrine senescent cells in controlling the spread of senescence. We demonstrate the crucial importance of the dynamic nature of senescence induction in generating juxtacrine secondary senescent cells before paracrine secondary senescent cells and limiting the size of senescent lesions over extended periods. As juxtacrine cells do not emit a SASP (or their SASP is blunted), the presence of juxtacrine senescent cells both reduces the concentration of SASP and limits the parameter regime in which senescence can spread, as their presence limits which cells can become paracrine senescent. Additionally, we propose an experiment aimed at determining four critical parameters involved in senescence spread, and we show that Bayesian inference is a viable tool for estimating these parameters from experimental data.Although determined from a 2D experiment, several of the inferred parameters (ND, NE, Rtot, κon) also describe paracrine senescence spread in 3D. We believe our model of senescence spread is well‐justified to study the in vitro scenario and of an appropriate complexity to identify which assumptions are important or affect senescence spread. Once parameterized with future in vitro experiments, the model can straightforwardly be expanded to 3D, to describe senescence spread in tissues, in vivo.While we can quantify parameters associated with the paracrine spread of senescence from a lesion, determining parameters associated with juxtacrine secondary senescent cells is more challenging.Single‐cell transcriptomic data may help to determine the difference between the senescent cell types. However, it is likely that differences between senescent subtypes will vary depending on the cell type and method of senescence induction, suggesting that senescence spread might differ depending on the location in the body and the cause of senescence.By considering a minimal model of senescence spread we have extracted an understanding of the dynamics of this system. Understanding the way senescence spreads through a population of cells can help us understand several biological processes, from the senescent lesion created in the radiation treatment of cancers to the accumulation of senescent cells as we age. Some questions remain unresolved, such as, whether paracrine senescence cells produce the same amount of SASP as primary senescent cells. These questions may be resolved by specially designed in vitro and in vivo experiments in combination with our modelling and inference framework.METHODSMinimal modelAutocrine and paracrine signalling in cell cultureTo understand the limiting factors in the dynamics of senescence spread, we first considered a minimal mathematical model.To mathematically describe signalling in a 2D layer of cells covered by a 3D layer of media (Figure 1c), we required knowledge of how a ligand produced by a cell will diffuse through a layer of media. Equation (1), from Batsilas et al. (2003), gives the cumulative probability of a ligand produced by a cell binding at distance r from the emitting cell.To arrive at Equation (1), Batsilas et al. assumed that ligands reflect from the upper surface of the media on contact and that ligands can bind to receptors on the surface of cells, but they reflect from the space between cells (the dish). Thus, they approximated the cell‐covered dish as a semi‐reflecting boundary layer (which is less accurate between 1<r/rcell<2, close to the emitting cell) and simulated the diffusion and binding of a large number of emitted ligands.We recreated the simulations described by Batsilas et al. and verified that Equation (1) predicts the binding location of ligands with a relative error <5%. For paracrine trajectories, we compared the time at which a simulated ligand first binds to a new cell to the square of the mean distance at which it binds (Figure S4b,c), showing that the limiting timescale in this process is diffusion (Müller & Schier, 2011). The shape of the cumulative binding probability (Pr) is shown in Figure 1d, for both small and large numbers of receptors per cell, Rtot.Equation (1) describes paracrine trajectories, not including ligands that bind to the cell that emitted them (autocrine trajectories). The probability of an autocrine trajectory is given as,2Pau=DaDa+4/π,where Da is defined in Table 2 (Batsilas et al., 2003).2TABLETable showing compound parameters from Batsilas et al. (2003).ParameterMeaningFormulaνProbability of internalisationκe/κe+κoffκRate constantκonRtotπrcell2NADaDamköhler numberrcellκDLκeffEffective trapping rate for a homogeneous surfaceκσ1+πDa/4Equations (1) and (2) do not account for the dissociation of ligands from a cell. If we assume ligands can be internalised by cells via endocytosis with rate κe and dissociate with rate κoff, the equations can be modified to describe the final distance at which ligands are bound,3Pinr≈rr+1.1DL/νκeff,andPauin=νDaνDa+4/π,where ν=κe/κe+κoff. The first binding event and internalisation are the same if κoff=0 and approach each other when the rate of endocytosis is large and the rate of dissociation is small.Between the first binding event and internalisation a ligand could bind to many cells. However, this binding may be only for a short period of time if the rate of dissociation is high. In this paper, we consider only the location of the final binding of the ligand to be important for senescence induction. This assumption breaks down if even fleeting binding of ligands to a cell produces a signal that persists intracellulary, and when combined with the signal from many other fleeting ligand binding events also induces senescence. A comparison of the binding and internalisation distances for an example set of parameters is given in Figure S6.The diffusion of ligands is clearly affected by the height of the media layer, h. However, to sustain cells in vitro, we would require a layer of media with h>2 mm. Using dimensional analysis and biologically reasonable system parameters, Batsilas et al. (2003) argue that the characteristic trajectory length, DL/κ, is much smaller than 2 mm, and so h does not need to be accounted for.Senescence spread from a single cell to an annulusExtending the work of Batsilas et al., we find that for each ligand produced, the probability of a ligand binding in an annulus of width 2rcell (illustrated in Figure 1e) at distance r from the emitting cell is,41−Pau×r+rcellr+rcell+a−r−rcellr−rcell+a=1−Pau2arcellr2−rcell2+2arcell+a2,where a=1.1DL/νKeff. The area of this annulus is given by 4πrrcell.To find the probability that a cell in this annulus will bind the ligand, we multiply this probability by the area of one cell relative to the annulus, to get,51−Pauarcell22rr2−rcell2+2arcell+a2.If one cell emits NE ligands per unit of time, the expected number of ligands bound to a cell at a distance r, per unit of time, is,6λr=1−PauNEarcell22rr2−rcell2+2arcell+a2,assuming that diffusion has reached a steady state, and NE is constant so that the rate of ligands reaching a cell at distance r is also constant. This is justified if the diffusion timescale is much faster than the other timescales in the system.To find the probability per unit time that a specific cell in an annulus will become senescent in any given second, λcellr, we assume that a cell must bind to a threshold number (ND) of ligands per unit of time to become senescent. Here, we use the binomial distribution to describe this process, where NE becomes the number of trials, and Equation (5) gives the probability of success.However, there is more than one cell in an annulus, on average there will be m=4rσrcell. Therefore, λspreadr=mλcellr is the rate of senescence induction in each annulus.Following this analysis through, we find the probability of creating at least one senescent cell at distance r from the primary lesion in time t is,7PrT≤t=1−e−λspreadrt,which is 1, minus the probability (given by a Poisson process) that no cells become senescent.Total senescence spread from a single cellWe now have the probability that one cell at distance r from the primary lesion will become senescent in a time period. We use the Poisson‐Binomial distribution to calculate the probability that a senescent cell is created in at least one annulus, as the probability of senescence induction in each annulus is different. Therefore, the probability of senescence induction over all annuli is given by a set of independent Bernoulli trials that are not identically distributed, where each annulus is a Bernoulli trial with p=PrT≤t, and we calculate the probability that one or more of these trials is successful.We use this simplified model of the system to determine whether the spread of senescence is uncontrollable, considering the limiting case of a single primary senescent cell. If a single senescent cell spreads senescence to one other cell before the body can remove it (we approximate that this happens in 2 days), with a probability of ~1, we define the spread of senescence as uncontrollable. This uncontrollable spread of senescence is not observed in vitro; therefore, this is an unrealistic parameter regime.Senescence spread with juxtacrine cellsWe considered two different regimes, with and without juxtacrine secondary senescence. Where juxtacrine signalling induces secondary senescence, we expect that cells surrounding the primary senescent cell will become juxtacrine senescent through contact. Assuming that juxtacrine secondary senescent cells do not emit a SASP, for the senescence to spread, a cell beyond this ring of juxtacrine cells must become senescent within 2 days. If juxtacrine cells instead produce a blunted SASP, this will lead to an underestimate of the spread of senescence (Rattanavirotkul et al., 2020).The distance of the annulus from the primary senescent cell which is considered in each case is given by,8r1=πrcell2σ,or,r2=2πrcell2σ.The first value, r1, is simply the average distance between cells, and r2=2r1, accounts for a ring of juxtacrine senescent cells.LimitationsThis minimal mathematical model provides a tractable way to consider senescence spread from a single cell but is difficult to extend to a system containing many cells. Furthermore, this model cannot be used to investigate the effect of a time delay in senescence induction, and it makes a number of assumptions which may limit the accuracy of the predictions.Stochastic simulationsTo deepen our understanding of the system and more closely match the reality of senescence spread in vivo, we created a stochastic model for the spread of senescence from multiple initial cells and over many subsequent senescence induction events. The stochastic model lets us introduce the added complexity of paracrine senescent cells which can in turn emit a SASP. We used the simulation to explore the change in dynamics resulting from changing the strength of the SASP produced by these secondary senescent cells. Finally, we introduce a time delay between senescence induction and SASP, more accurately modelling the cellular biology. Cells do not immediately start emitting a SASP once they have bound enough ligands to induce senescence, and instead, take several days to become fully senescent.Simulation implementationWe consider three possible cell states: senescent, in the process of becoming senescent, and non‐senescent. A cell starts to become paracrine secondary senescent when it binds to more than ND ligands per unit of time (see Section 2.1.1). The simulation is initialised on a square plane, and a 100 × 100 cell radii area is commonly used as larger areas increase the computational time taken for each simulation. Cells are then created so that they randomly populate this area with a specified density (Figure 3a). Primary senescent cells are seeded at time 0. The simulation proceeds by tracking the state of each cell.In the stochastic simulation, we approximate the number of ligands bound from each emitting senescent cell using the Poisson distribution. This approximation holds when the number of ligands emitted, NE, is large and the probability of binding (equation (5)) is small. Large binding probabilities are not of interest as they would lead to the fast and uncontrolled spread of senescence. The number of ligands bound from all senescent cells is given by a Poisson distribution whose mean is the sum of the Poisson distribution of ligands from each emitting cell (equation 6),9λtotr=∑i=1Sλir,where S is the number of senescent cells emitting a SASP.This approximation dictates the appropriate time frame for the simulation. Considering paracrine signalling using only IL‐6, IL‐8 and Activin A molecules, we can estimate NE≈2887h−1.We can now calculate the probability of a cell becoming senescent in one unit of time from its position relative to the SASP‐emitting cells. This probability can then be used to give the rate of senescence induction for that cell. Thus, each cell on the grid which is not yet senescent has a rate of senescence induction. By summing these rates we calculate the time at which the next cell becomes paracrine secondary senescent. We use the Gillespie algorithm, meaning that we first determine the time an event occurs based on creation rates, for example, a cell becomes paracrine secondary senescent at t = 1.5 h. We then determined which cell becomes senescent, again based on the rate of senescence induction associated with each cell (Gillespie, 1977).The situation is complicated further by the time delay: It may be that cells have finished becoming fully senescent before the next paracrine induction event, in which case the SASP that they produce must also be accounted for. This is achieved simply by tracking the number of cells that are in the process of becoming senescent. If a cell finishes becoming senescent before the next reaction (i.e., new senescence induction event) occurs, then the cell is added to the list of fully senescent cells, and the simulation run again from that time point, with the probabilities of each cell becoming senescent re‐calculated including any SASP from the newest senescent cell. This works as the system is Markovian.There are a number of algorithms that can be used to speed up the delayed reaction calculation (Cai, 2007). However, at any one time, we do not have a large number of cells waiting to become senescent and the bottleneck in the process is the repeated calculation of the probability of inducing senescence in each cell.The simulation code can be found on GitHub (https://github.com/lkmartin90/Senescence_Spread), and a flowchart of the stochastic simulation implementation is shown in Figure S10.Comparison of stochastic simulations to minimal modelWe compared our stochastic simulation to our minimal model to verify that the two agree. To do this, we considered the spread of senescence from a single cell, in the case where no juxtacrine secondary senescent cells are created.We used the stochastic simulation to simulate this scenario 100 times (fixing all parameters but initialising a new population of cells each time) and recorded the proportion of the simulations in which there was no senescence spread. We varied ND and Rtot, to ensure that the simulations agree over a range of parameters. We compared the results with our minimal model, which gives us the probability of senescence spread.From this, we see a good agreement between the two methods at both Rtot=105 and Rtot=106, both when cells in the stochastic simulation are created so that they lie on a lattice (Figure S3a,b), and when cells are randomly distributed (Figure S3c,d). There are some differences in the regime without juxtacrine senescence induction with Rtot=106, due to the large number of ligands which bind to nearby cells and the assumption we make when constructing the minimal model that cells are uniformly distributed (Figure S3b, dotted).In the stochastic simulation, with cells distributed randomly (Figure S3c), the range of parameters that result in a probability of senescence spread >0 and <1is greater. If a nearby cell happens to be closer than the average distance between cells then there is a higher probability that senescence will spread. This implies that the parameter regime of large Rtot is highly dependent on the initialisation of the simulation, in particular, the precise location of the cells. Although the average density in each simulation is the same, the signal spread is possible in some cases and not others. If there are no cells close enough to the primary senescent cell, senescence will not spread.We investigated whether the optimum number of receptors that we found for senescence spread in the minimal model, as shown in Figure 2, also appears in the stochastic simulation, or whether this is a consequence of the simplifying assumptions made in the minimal model. We see that the random nature of the cell placement in the stochastic simulation means that there is no longer a clear optimal number of receptors to ensure the spread of senescence, however, the general dynamics remain the same, with juxtacrine cells reducing the spread of senescence and the number of receptors altering the distance over which cells communicate (Figure S3d).LimitationsThe stochastic simulation is a useful tool to understand senescence spread from a small lesion, but computationally intensive to scale to larger spatial regions. These limitations, as well as the relationships between the models discussed in this paper, are shown in Figure S9.Parameter inferenceWe use Bayesian statistics to infer system parameters from a realistic in vitro experiment, where a circle of primary senescent cells is seeded in the centre of a dish, as described in Section 2.7. We calculate the distance senescence spreads from the circular lesion.The number of SASP molecules emitted per unit area by a primary senescent lesion is, NEσπrcell2, and the probability that a single test cell (Figure 4a) at distance x from a senescent cell binds a ligand emitted by a senescent cell is101−Pauarcell22xx2−rcell2+2arcell+a2.Integrating for the emission from the area of seeded cells to get the total expected number of ligands binding to the test cell gives111−PauNEσπrcell2arcell22∫02π∫0rseedr′xx2−rcell2+2arcell+a2dθdr′.Writing x in terms of r, r′ and θ, according to the geometry in Figure 4a, gives121−PauNEσπrcell2∫02π∫0rseed1r′2+r2−2rr′cosθ0.5×r′r′2+r2−2rr′cosθ−rcell2+2arcell+a2dθdr′.Here, rseed is the radius of the circle of primary senescent cells. We solve this integral numerically.Using the method previously described in Section 4.1, we calculate the probability that the cell that binds these ligands becomes senescent, using a Poisson distribution to determine the probability that more than ND ligands are bound per hour. We then calculate the probability that a cell becomes senescent within the duration of the experiment. Finally, we calculate the fraction of cells, which are senescent at a given distance from the seeded senescent region of cells. As this process is probabilistic, the fraction of senescent cells varies on each run of the model.Before comparing the model output with an experiment, we ought to consider whether the assumptions made in the model's construction are reasonable. In particular, in the model, we have not accounted for the time taken for ligands to diffuse away from the cell that emitted them. From Figure S4, we can see that the overall process of signal spread is diffusion‐limited. Previous experiments show senescence spread over a distance of the order of 1 mm, indicating that this is a meaningful distance to consider. The time taken for ligands to travel 1 mm is between 50 and 100 seconds for the parameter regimes shown, depending on the diffusion coefficient of the ligand. As long as the time taken for diffusion is small compared with the timescales associated with senescence spread (hours), we can assume that the diffusion occurs instantaneously.To determine the parameters that lead to the closest agreement between the model and the experimental data, we compare two summary statistics. First, the distance where the fraction of senescent cells crosses 0.5 (S). Second, the fraction of senescent cells at different distances from the lesion (F→). The distances need not be equally spaced, and we use more data points around the region of r=S, where the fraction of senescent cells is <1 and >0. We verified that these summary statistics are sensitive to changes in the parameters (Figure S8a).To quantify the agreement between predicted and observed summary statistics, use the distance function:13d=d12+d22,where d1=Smodel−Sdata and d2=∥F→model−F→data∥.We use pyABC (Klinger et al., 2018), a likelihood‐free inference software using a sequential Monte Carlo scheme, to estimate posterior distributions of model parameters. In this paper, we use synthetic data, created from the same model we use for the inference.Our inference ran on 60 cores for 13 populations, with 1000 particles, taking 7.5 h in total. While acceptance thresholds for d (epsilon values) have not converged at this stage, indicating that further iterations may improve the inference, the fraction of accepted samples was ≈10−3, making it computationally wasteful to run for more populations (Figure S8b).AUTHOR CONTRIBUTIONSL.S. and T.C. supervised the study. L.K.M., L.S. and T.C. wrote the manuscript. L.K.M. and L.S. developed the modelling framework. L.K.M. conducted the analysis.ACKNOWLEDGMENTSL.K.M is a cross‐disciplinary research fellow supported by funding from the CRUK Brain Tumour Centre of Excellence Award (C157/A27589). T.C. is supported through a Chancellor's Fellow at the University of Edinburgh and the MRC Human Genetics Unit. We thank M. Nicholson for the critical reading of the manuscript. We thank all members of the Chandra lab and Schumacher lab for their input. We thank Acosta et al. (2013) and Schafer et al. (2020) for sharing additional data related to their publications on request.CONFLICT OF INTEREST STATEMENTThe authors declare no conflict of interest.DATA AVAILABILITY STATEMENTCode used in this manuscript is available at https://github.com/lkmartin90/Senescence_Spread.REFERENCESAcosta, J. C., Banito, A., Wuestefeld, T., Georgilis, A., Janich, P., Morton, J. P., Athineos, D., Kang, T. W., Lasitschka, F., Andrulis, M., Pascual, G., Morris, K. J., Khan, S., Jin, H., Dharmalingam, G., Snijders, A. P., Carroll, T., Capper, D., Pritchard, C., … Gil, J. (2013). A complex secretory program orchestrated by the inflammasome controls paracrine senescence. Nature Cell Biology, 15(8), 978–990. https://doi.org/10.1038/ncb2784Acosta, J. C., O'Loghlen, A., Banito, A., Guijarro, M. V., Augert, A., Raguz, S., Fumagalli, M., Da Costa, M., Brown, C., Popov, N., Takatsu, Y., Melamed, J., d'Adda di Fagagna, F., Bernard, D., Hernando, E., & Gil, J. (2008). Chemokine signaling via the CXCR2 receptor reinforces senescence. Cell, 133(6), 1006–1018. https://doi.org/10.1016/j.cell.2008.03.038Baker, D. J., Childs, B. G., Durik, M., Wijers, M. E., Sieben, C. J., Zhong, J., Saltness, R., Jeganathan, K. B., Versoza, G. C., Pezeshki, A.‐M., Khazaie, K., Miller, J. D., & van Deursen, J. M. (2016). Naturally occurring p16ink4a‐positive cells shorten healthy lifespan. Nature, 530(7589), 184–189. https://doi.org/10.1038/nature16932Baker, D. J., Wijshake, T., Tchkonia, T., LeBrasseur, N. K., Childs, B. G., van de Sluis, B., Kirkland, J. L., & van Deursen, J. M. (2011). Clearance of p16ink4a‐positive senescent cells delays ageing‐associated disorders. Nature, 479(7372), 232–236. https://doi.org/10.1038/nature10600Bartkova, J., Rezaei, N., Liontos, M., Karakaidos, P., Kletsas, D., Issaeva, N., Vassiliou, L. V. F., Kolettas, E., Niforou, K., Zoumpourlis, V. C., Takaoka, M., Nakagawa, H., Tort, F., Fugger, K., Johansson, F., Sehested, M., Andersen, C. L., Dyrskjot, L., Ørntoft, T., … Gorgoulis, V. G. (2006). Oncogene‐induced senescence is part of the tumorigenesis barrier imposed by DNA damage checkpoints. Nature, 444(7119), 633–637. https://doi.org/10.1038/nature05268Basisty, N., Kale, A., Jeon, O. H., Kuehnemann, C., Payne, T., Rao, C., Holtz, A., Shah, S., Sharma, V., Ferrucci, L., Campisi, J., & Schilling, B. (2020). A proteomic atlas of senescence‐associated secretomes for aging biomarker development. PLoS Biology, 18(1), e3000599. https://doi.org/10.1371/journal.pbio.3000599Batsilas, L., Berezhkovskii, A. M., & Shvartsman, S. Y. (2003). Stochastic model of autocrine and paracrine signals in cell culture assays. Biophysical Journal, 85(6), 3659–3665. https://doi.org/10.1016/S0006‐3495(03)74783‐3Burd, C. E., Sorrentino, J. A., Clark, K. S., Darr, D. B., Krishnamurthy, J., Deal, A. M., Bardeesy, N., Castrillon, D. H., Beach, D. H., & Sharpless, N. E. (2013). Monitoring tumorigenesis and senescence in vivo with a p16ink4a‐luciferase model. Cell, 152, 340–351. https://doi.org/10.1016/j.cell.2012.12.010Cai, X. (2007). Exact stochastic simulation of coupled chemical reactions with delays. The Journal of Chemical Physics, 126(12), 124108. https://doi.org/10.1063/1.2710253Copp'e, J.‐P., Patil, C. K., Rodier, F., Sun, Y., Muñoz, D. P., Goldstein, J., Nelson, P. S., Desprez, P.‐Y., & Campisi, J. (2008). Senescence‐associated secretory phenotypes reveal cell‐nonautonomous functions of oncogenic RAS and the p53 tumor suppressor. PLoS Biology, 6(12), e301. https://doi.org/10.1371/journal.pbio.0060301Cristofalo, V. J., Volker, C., & Allen, R. G. (2003). Use of the fibroblast model in the study of cellular senescence. In Aging methods and protocols (pp. 23–52). Humana Press.d'Adda di Fagagna, F., Reaper, P. M., Clay‐Farrace, L., Fiegler, H., Carr, P., Von Zglinicki, T., Saretzki, G., Carter, N. P., & Jackson, S. P. (2003). A DNA damage checkpoint response in telomere‐initiated senescence. Nature, 426(6963), 194–198. https://doi.org/10.1038/nature02118Demaria, M., Ohtani, N., Youssef, S. A., Rodier, F., Toussaint, W., Mitchell, J. R., Laberge, R.‐M., Vijg, J., Van Steeg, H., Doll'e, M. E. T., Hoeijmakers, J. H. J., de Bruin, A., Hara, E., & Campisi, J. (2014). An essential role for senescent cells in optimal wound healing through secretion of PDGF‐AA. Developmental Cell, 31(6), 722–733. https://doi.org/10.1016/j.devcel.2014.11.012Gillespie, D. T. (1977). Exact stochastic simulation of coupled chemical reactions. Journal of Physical Chemistry, 81(25), 2340–2361. https://doi.org/10.1021/j100540a008Hayflick, L., & Moorhead, P. S. (1961). The serial cultivation of human diploid cell strains. Experimental Cell Research, 25(3), 585–621. https://doi.org/10.1016/0014‐4827(61)90192‐6Hoare, M., Ito, Y., Kang, T.‐W., Weekes, M. P., Matheson, N. J., Patten, D. A., Shetty, S., Parry, A. J., Menon, S., Salama, R., Antrobus, R., Tomimatsu, K., Howat, W., Lehner, P. J., Zender, L., & Narita, M. (2016). NOTCH1 mediates a switch between two distinct secretomes during senescence. Nature Cell Biology, 18(9), 979–992. https://doi.org/10.1038/ncb3397Ito, Y., Hoare, M., & Narita, M. (2017). Spatial and temporal control of senescence. Trends in Cell Biology, 27, 820–832. https://doi.org/10.1016/j.tcb.2017.07.004Jochems, F., Thijssen, B., De Conti, G., Jansen, R., Pogacar, Z., Groot, K., Wang, L., Schepers, A., Wang, C., Jin, H., Beijersbergen, R. L., Leite de Oliveira, R., Wessels, L. F., & Bernards, R. (2021). The cancer SENESCopedia: A delineation of cancer cell senescence. Cell Reports, 36(4). https://doi.org/10.1016/J.CELREP.2021.109441Katzir, I., Adler, M., Karin, O., Mendelsohn‐Cohen, N., Mayo, A., & Alon, U. (2021). Senescent cells and the incidence of age‐related diseases. Aging Cell, 20, e13314. https://doi.org/10.1111/acel.13314Kirschner, K., Rattanavirotkul, N., Quince, M. F., & Chandra, T. (2020). Functional heterogeneity in senescence. Biochemical Society Transactions, 48(3), 765–773. https://doi.org/10.1042/BST20190109Kiss, T., Nyúl‐Tóth, A., DelFavero, J., Balasubramanian, P., Tarantini, S., Faakye, J., Gulej, R., Ahire, C., Ungvari, A., Yabluchanskiy, A., Wiley, G., Garman, L., Ungvari, Z., & Csiszar, A. (2022). Spatial transcriptomic analysis reveals inflammatory foci defined by senescent cells in the white matter, hippocampi and cortical grey matter in the aged mouse brain. GeroScience, 2022, 1–21. https://doi.org/10.1007/S11357‐022‐00521‐7Klinger, E., Rickert, D., & Hasenauer, J. (2018). PyABC: distributed, likelihoodfree inference. Bioinformatics, 34(20), 3591–3593. https://doi.org/10.1093/bioinformatics/bty361Kuilman, T., Michaloglou, C., Vredeveld, L. C. W., Douma, S., van Doorn, R., Desmet, C. J., Aarden, L. A., Mooi, W. J., & Peeper, D. S. (2008). Oncogene‐induced senescence relayed by an interleukin‐dependent inflammatory network. Cell, 133(6), 1019–1031. https://doi.org/10.1016/j.cell.2008.03.039Lee, P. J., Benz, C. C., Blood, P., Börner, K., Campisi, J., Chen, F., Daldrup‐Link, H., De Jager, P., Ding, L., Duncan, F. E., Eickelberg, O., Fan, R., Finkel, T., Furman, D., Garovic, V., Gehlenborg, N., Glass, C., Heckenbach, I., Joseph, Z.‐B., … Yale, T. M. C. (2022). NIH SenNet consortium to map senescent cells throughout the human lifespan to understand physiological health. Nature Aging, 2(12), 1090–1100. https://doi.org/10.1038/s43587‐022‐00326‐5Michaloglou, C., Vredeveld, L. C. W., Soengas, M. S., Denoyelle, C., Kuilman, T., van der Horst, C. M. A. M., Majoor, D. M., Shay, J. W., Mooi, W. J., & Peeper, D. S. (2005). BRAFE600‐associated senescence‐like cell cycle arrest of human naevi. Nature, 436(7051), 720–724. https://doi.org/10.1038/nature03890Müller, P., & Schier, A. F. (2011). Extracellular movement of signaling. Molecules, 21, 145–158. https://doi.org/10.1016/j.devcel.2011.06.001Muñoz‐Espín, D., & Serrano, M. (2014). Cellular senescence: From physiology to pathology. Nature Reviews Molecular Cell Biology, 15(7), 482–496. https://doi.org/10.1038/nrm3823Narita, M., Nuñez, S., Heard, E., Narita, M., Lin, A. W., Hearn, S. A., Spector, D. L., Hannon, G. J., & Lowe, S. W. (2003). Rb‐mediated heterochromatin formation and silencing of e2f target genes during cellular senescence. Cell, 113(6), 703–716. https://doi.org/10.1016/s00928674(03)00401‐xPassos, J. F., Simillion, C., Hallinan, J., Wipat, A., & von Zglinicki, T. (2009). Cellular senescence: unravelling complexity. Age, 31(4), 353–363. https://doi.org/10.1007/s11357‐009‐9108‐1Prata, L. G. P. L., Ovsyannikova, I. G., Tchkonia, T., & Kirkland, J. L. (2018). Senescent cell clearance by the immune system: Emerging therapeutic opportunities. Seminars in Immunology, 40, 101275. https://doi.org/10.1016/j.smim.2019.04.003Putavet, D. A., & de Keizer, P. L. J. (2021). Residual disease in glioma recurrence: A dangerous liaison with senescence. Cancers, 13(7). https://doi.org/10.3390/CANCERS13071560Rattanavirotkul, N., Kirschner, K., & Chandra, T. (2020). Induction and transmission of oncogene‐induced senescence. Cellular and Molecular Life Sciences, 78, 843–852. https://doi.org/10.1007/s00018‐020‐03638‐0Schafer, M. J., Zhang, X., Kumar, A., Atkinson, E. J., Zhu, Y., Jachim, S., Mazula, D. L., Brown, A. K., Berning, M., Aversa, Z., Kotajarvi, B., Bruce, C. J., Greason, K. L., Suri, R. M., Tracy, R. P., Cummings, S. R., White, T. A., & LeBrasseur, N. K. (2020). The senescence‐associated secretome as an indicator of age and medical risk. JCI Insight, 5(12). https://doi.org/10.1172/jci.insight.133668Serrano, M., Lin, A. W., McCurrach, M. E., Beach, D., & Lowe, S. W. (1997). Oncogenic ras provokes premature cell senescence associated with accumulation of p53 and p16INK4a. Cell, 88(5), 593–602. https://doi.org/10.1016/S0092‐8674(00)81902‐9Shvartsman, S. Y., Steven Wiley, H., Deen, W. M., & Lauffenburger, D. A. (2001). Spatial range of autocrine signaling: Modeling and computational analysis. Biophysical Journal, 81(4), 1854–1867. https://doi.org/10.1016/S0006‐3495(01)75837‐7Teo, Y. V., Rattanavirotkul, N., Olova, N., Kirschner, K., Neretti, N., & Correspondence, C. (2019). Notch signaling mediates secondary senescence. Cell Reports, 27, 997–1007. https://doi.org/10.1016/j.celrep.2019.03.104van Deursen, J. M. (2014). The role of senescent cells in ageing. Nature, 509(7501), 439–446. https://doi.org/10.1038/nature13193Wang, B., Kohli, J., & Demaria, M. (2020). Senescent cells in cancer therapy: friends or foes? Trends in Cancer, 6, 838–857. https://doi.org/10.1016/j.trecan.2020.05.004Wiley, H. S., Shvartsman, S. Y., & Lauffenburger, D. A. (2003). Computational modeling of the EGF‐receptor system: A paradigm for systems biology. Trends in Cell Biology, 13(1), 43–50. https://doi.org/10.1016/S0962‐8924(02)00009‐0Young, A. R. J., Narita, M., Ferreira, M., Kirschner, K., Sadaie, M., Darot, J. F. J., Tavaré, S., Arakawa, S., Shimizu, S., Watt, F. M., & Narita, M. (2009). Autophagy mediates the mitotic senescence transition. Genes & Development, 23(7), 798–803. https://doi.org/10.1101/gad.519709 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Aging Cell Wiley

Modelling the dynamics of senescence spread

Aging Cell , Volume 22 (8) – Aug 1, 2023

Loading next page...
 
/lp/wiley/modelling-the-dynamics-of-senescence-spread-b4M2SYFvV0

References (44)

Publisher
Wiley
Copyright
Copyright © 2023 The Anatomical Society and John Wiley & Sons Ltd
eISSN
1474-9726
DOI
10.1111/acel.13892
Publisher site
See Article on Publisher Site

Abstract

AbbreviationsABCapproximate bayesian computationEGFepidermal growth factorSASPsenescence‐associated secretory phenotypeINTRODUCTIONSenescence is a cellular stress response that can be activated by telomere shortening, DNA damage and the up‐regulation of oncogenes (d'Adda di Fagagna et al., 2003; Hayflick & Moorhead, 1961; Serrano et al., 1997). Senescent cells display cell cycle arrest and several changes to the cell, including large‐scale chromatin rearrangement and the production of secretory proteins, which can induce senescence in nearby cells (the SASP, senescence associate secretory phenotype) (Acosta et al., 2008; Basisty et al., 2020; Copp'e et al., 2008; Kuilman et al., 2008; Narita et al., 2003). Senescence is necessary for many processes in tissues including development, wound healing and defence against tumour formation (Demaria et al., 2014; Kirschner et al., 2020; Muñoz‐Espín & Serrano, 2014; Passos et al., 2009). When senescent cells accumulate, the molecules they release can create a pro‐tumourigenic and inflammatory microenvironment, which has been linked to both the ageing process and the recurrence of cancerous tumours following treatment (Wang et al., 2020).SASP molecules can spread the senescence phenotype through paracrine signalling, with molecules diffusing from a secreting senescent cell and binding to non‐senescent cells. Several experiments have demonstrated paracrine transmission of senescence, for example, senescent cells contained within a transwell induce senescence in a population of growing cells sharing the same media (Acosta et al., 2013). SASP signalling also recruits immune cells which can clear senescent cells (Prata et al., 2018).The terms “primary” and “secondary” senescence are used to distinguish between cells in which senescence has been induced cell‐intrinsically (e.g., through DNA damage) and extrinsically (through signals from primary senescent cells), respectively (Figure 1b). Recent evidence shows that these subtypes of senescence are transcriptomically distinct (Rattanavirotkul et al., 2020; Teo et al., 2019).1FIGUREModeling the spread of senescence. (a) A controlled senescence spread is beneficial in vivo, but it is unclear how a limited spread of senescence is achieved. (b) Illustration of the cellular phenotypes considered in our model. (c) A cartoon of the simulated diffusion of ligands through a 3D layer of media over a 2D layer of cells. (d) The cumulative probability of a ligand binding at a distance, r, from the cell that emitted it. (e) An illustration of the minimal model of senescence spread developed in this paper, where we consider the probability of inducing senescence in the orange annulus. (f) An illustration of selected model parameters. (g) The ways senescent cells can induce senescence in the cells around them.In addition to paracrine signalling (SASP), senescence can also spread via juxtacrine communication between cells that are in direct contact (NOTCH‐mediated) (Hoare et al., 2016; Teo et al., 2019). Secondary senescent cells created through juxtacrine signalling may differ from secondary senescent cells created through paracrine signalling, with juxtacrine secondary cells producing fewer SASP proteins, alongside other transcriptional differences (Hoare et al., 2016; Rattanavirotkul et al., 2020). Whether these two types of secondary senescent cells play different roles in the spread of senescence is as yet unclear (Kirschner et al., 2020).Here, we define the controlled spread of senescence as a local (finite) senescence spread over an entire human lifespan. In vivo, a limited and controlled senescence spread is thought to play an important role in at least two scenarios. First, any secondary senescent cells created around a primary senescent lesion may help to amplify the signal to the immune system. Second, as a failsafe mechanism, damaged cells that failed to undergo primary senescence will become senescent in response to the primary cells (Figure 1a). However, multiple lines of evidence suggest that an accumulation of senescent cells and uncontrolled spread will ultimately lead to tissue dysfunction (Baker et al., 2011, 2016; Burd et al., 2013; Katzir et al., 2021).Currently, there are no hypothesised mechanisms to contain the spread of senescence in the absence of the immune system. However, it has been shown that the spread of senescence is contained when immune cells are not present (Acosta et al., 2013). In vitro experiments have shown that from a seeded circle of senescent cells, the spread of senescence is finite (senescence was found at a distance of approximately 1 mm from seeded cells) (Acosta et al., 2013). Although the medium taken from primary senescent fibroblasts is capable of inducing secondary senescence in growing fibroblasts, it is uncertain if tertiary senescence is induced when the process is repeated (Acosta et al., 2013). Controlled senescence spread is also observed in vivo, where some senescent lesions persist over time (melanocytic naevi) and are not cleared by the immune system, but the senescence phenotype does not continue to spread (Bartkova et al., 2006). This is not due to SASP production halting, as cells continue to produce a SASP even months after the onset of senescence and enter a deep senescence state, which can lead to chronic senescence in vivo (van Deursen, 2014).In this paper, we use both a minimal mathematical model and stochastic simulation to investigate the mechanisms that determine the dynamics of senescence spread. We identify that juxtacrine‐induced senescent cells can limit senescence spread by acting as a fire break in the wildfire of paracrine signalling and that this effect is more pronounced when SASP‐signalling is delayed compared with juxtacrine signalling. We propose an in vitro experiment and demonstrate that the parameters of senescence spread can be inferred using Bayesian statistics. With the emergence of large‐scale descriptive efforts such as the SenNet consortium (US) (Lee et al., 2022) and the advancement of spatial transcriptomics (Kiss et al., 2022), our model can serve as a framework for further exploration as new data become available.RESULTSSenescence spread modelWe considered senescence spread in a 2D layer of cells covered by a 3D layer of media so that we could compare our results with senescence spread in vitro, as senescence is best characterised in adherent human fibroblasts (Cristofalo et al., 2003; Hayflick & Moorhead, 1961). Despite the 2D nature of our model, we can use it to infer parameters that are unchanged between 2D and 3D models of senescence spread, such as the number of ligands needed to induce senescence in a cell.To simplify the system, we described the cells as circular and non‐motile, evenly and randomly distributed on the 2D plane with a specified density. Following work by Batsilas et al. (2003), we considered the cells to be covered by a layer of media of height h (as shown in Figure 1c). The density of the cells in the 2D plane is given by σ, and the radius of each cell is given by rcell. The number of receptors per cell is Rtot.Paracrine signalling is mediated by ligands, and in our model, the ligands are SASP molecules. Ligands are emitted by senescent cells and diffuse through the surrounding media. DL is the diffusion coefficient for the ligand, and a larger diffusion coefficient results in ligands travelling further in a given time frame. Three rates are used to describe a ligand binding to a receptor; κon gives the rate at which ligands bind to a receptor, κe the rate of endocytosis of bound ligand‐receptor complexes and κoff the rate at which ligands dissociate from receptors.Batsilas et al. describe the cumulative probability that a ligand binds at a distance, r, from the cell that emitted it, as1Pr≈rr+1.1DL/κeff.Here, κeff is the effective binding rate of ligands to cells (defined in Table 2). As κeff depends on several parameters, including the number of receptors per cell, the cumulative binding probability changes with Rtot and other system parameters (Figure 1d).Building on this work, we assume that cells have a fixed number of receptors for detecting ligands, senescent cells secrete ligands at a fixed rate (NE) and that a cell must bind a certain number of ligands (ND) to become senescent. We calculate the rate that ligands bind to a cell at distance r from the emitting cell and the rate that senescence is induced in cells at distance r from a senescent cell.Throughout this paper, we assume that juxtacrine secondary senescent cells do not produce a SASP and cannot, in turn, induce juxtacrine senescence in their neighbours as a blunted SASP response in juxtacrine secondary cells has been reported (Hoare et al., 2016; Teo et al., 2019). We explore different scenarios of senescence spread, such as secondary or tertiary induction of juxtacrine senescence from paracrine senescent cells, as summarised in Figure 1g. We assume that juxtacrine‐mediated senescence induction occurs when cells are in contact, that is, when they are in close spatial proximity.In a minimal model, we calculate the probability of senescence spread from a single senescent cell, by considering the probability of inducing senescence in an annulus of non‐senescent cells around the primary senescent cell, as shown in Figure 1e. By considering the probability of senescence spread to a series of annuli, each with the same width, but an increasing radius, we determine the likelihood that senescence spreads from a single cell.To model more realistic experimental scenarios we built a stochastic simulation of senescence spreading from a collection of cells over time. We use this to investigate the role of secondary senescent cells, by varying parameters such as their SASP secretion, time delays involved in senescence induction and lesion size.Detailed descriptions of the minimal model and stochastic simulations can be found in Methods (Section 4).Model parametersOur model of the spread of senescence contains several parameters, such as the constants relating to ligand–receptor interactions, where exact values have not been determined experimentally. Where this was the case we used values from closely related biological systems, for example, the parameters for EGF and its receptor suggested in Batsilas et al. (2003). These parameters are shown in Table 1.1TABLETable showing the estimated parameters for the spread of EGF, as described in Batsilas et al. (2003), originally from Wiley et al. (2003).ParameterMeaningEstimated valueNENumber of ligands emitted per senescent cell per unit time (h−1)VariableNDNumber of ligands required to induce senescence per unit time (h−1)UnknownDLDiffusion coefficient for ligands (cm2/s)10−6σDensity of cells (m2/m2)0.1–0.4rcellRadius of cells (m)10×10−6RtotTotal number of receptors per cell104–106κeEndocytosis rate (min−1)0.1–0.3κoffDissociation rate (min−1)0.1–0.3κonForward binding rate constant (M−1min−1)108We assume that a non‐senescent cell needs to bind a threshold number of ligands (ND) in a given time frame for senescence to be induced (Figure 1f) and that a single senescent cell produces NE ligands per unit time. We calculated the number of SASP proteins produced by a senescent cell per hour (Figure S7) based on a recent paper by Schafer et al. (2020) that measured the concentrations of SASP proteins produced by radiation‐induced senescent cells. For fibroblasts, the data suggest ~275,000 protein molecules per hour across all SASP proteins investigated. However, if we limit the proteins of interest to Interleukin 6 (IL‐6), Interleukin 8 (IL‐8) and Activin A (secreted proteins thought to be responsible for the induction of senescence; Acosta et al., 2008; Ito et al., 2017; Jochems et al., 2021; Prata et al., 2018; Putavet & de Keizer, 2021), the number produced per hour is much smaller, ~2887.For simplicity, we did not distinguish between the types of ligands produced. Although realistically, a specific ligand will bind to its complementary receptor, and multiple signalling molecules may be required to induce senescence, we did not include this complexity in the model.Cell communication distances are determined by the number of receptors per cell and the effective binding rateTo understand the impact of each model parameter on senescence spread we used our minimal model (described in Section 4.1) to scan through a range of parameters, noting the probability of senescence spread from a single primary senescent cell. We chose to fix the number of ligands emitted per hour by a cell (justification in Section 2.1.1), and vary ND, the number of ligands required to induce senescence. For small ND we observe senescence spread (blue region in Figure 2a) whereas high ND invariably prevents spread (red region in Figure 2a). Parameter regimes for a controlled (finite) spread are likely to be found at the transition between these regions. We, therefore, investigated how changing other parameters shifts the values of ND for which senescence spreads. Parameter values which increase the maximum of ND under which senescence can be induced in a nearby cell correspond to senescence being more likely to spread through the tissue (Figure 2a). Using this method we were able to draw several conclusions about the dynamics of paracrine signalling in a 2D layer of cells.2FIGUREA minimal model for senescence spread. (a) An illustration of the phase space of Rtot and ND, including a region where senescence spreads and one where senescence does not. (b) The division in phase space between senescence spread and no senescence spread for two values of κe and κoff, with and without juxtacrine senescence. The shaded region shows the parameter regime where senescence would spread without juxtacrine senescence, but not with juxtacrine senescence. (c) The range of parameters (shaded) for which senescence only spreads a fraction of the time. (d) The division in phase space between senescence spread and no senescence spread for a range of cell densities. (e) The division in phase space between senescence spread and no senescence spread for two values of Rtot, with and without juxtacrine senescence. The shaded region shows the parameter regime where senescence would spread without juxtacrine senescence, but not with juxtacrine senescence. (f) An illustration of the impact of receptor number on senescence spread. In this figure, unless specified, parameters were, DL=10−6cm2/s, rcell=10μm, κon=108, σ=0.35.We found that there is both an optimum number of receptors (Rtot) and an optimum forward binding rate (κon) for the spread of a signal (Figure 2a–e) in the regime where both κe and κoff (the endocytosis and dissociation rates) are negligible, such that ligands which bind to a receptor remain bound. Intuitively, it may seem that more receptors would increase the signal spread, but this is not the case. When ligands do not dissociate from a cell easily, once they are bound they are removed from the system. Therefore, in regimes with large numbers of receptors, ligands are more likely to bind to the cell that emitted them (known as autocrine signalling) or to bind in spatial proximity, so the signal does not spread. Conversely, in regimes with a small number of receptors, the signal diffuses to a much greater distance from the senescent cell, but it is unlikely that enough ligands will bind to a cell to signal effectively (Shvartsman et al., 2001). This is illustrated in Figure 2f.The same logic applies if κe>>κoff such that the probability or internalisation, ν, tends to 1, meaning that the rate of dissociation is relatively low. If the binding rate is large, many ligands will be sequestered in autocrine trajectories, spreading the signal ineffectively. However, if the binding rate is small it is unlikely that enough ligands will bind to a cell to induce a change. Figure 2e illustrates the effect of increasing the binding rate.When ligands can dissociate effectively and quickly from cells, the optimal number of receptors (for the signal to spread) is increased (Figure 2b; black vs. blue curves). In this regime, a single ligand can bind to multiple cells over the course of its lifetime. Thus, we assume only the last binding event is relevant for senescence induction. In this regime, senescence spreads further with a larger number of receptors and a greater binding rate.A similar effect is observed when increasing the diffusion coefficient (Figure 2c). With increased diffusivity a ligand can more easily move away from the cell that emitted it, there are fewer autocrine binding events, and so a larger number of receptors per cell can be accommodated without limiting the signal spread.Compared with the binding rate and receptor number, the cell density has only a small impact on the system's dynamics, with a greater cell density resulting in a greater chance of senescence spread (Figure 2d).Juxtacrine senescence cells limit the spread of senescenceTo investigate the effect of juxtacrine senescent cells on the spread of senescence, we first considered a ring of juxtacrine cells around the primary senescent cell in our minimal model. We assume juxtacrine secondary cells do not emit a SASP (Figure 1g). As expected, the layer of juxtacrine cells reduces the set of parameters for which senescence spreads (Figure 2b; solid vs. dotted curves), as ligands emitted by the primary senescent cell bind to the juxtacrine senescent cells, which are already senescent, but in this approximation do not emit ligands of their own. Thus, there is a parameter regime where senescence spread would occur without juxtacrine cells but does not when they are present (Figure 2b, shaded region).Using our minimal model we can find parameter regimes that lead to senescence spread from a single cell. However, we can not determine from this model whether the spread of senescence is controlled (finite) or uncontrolled. To further investigate the effect of juxtacrine senescent cells on the controlled spread of senescence over time, we created a stochastic simulation (described in Section 4.2) that also allowed us to investigate the spread from multiple primary senescent cells. Compared with the spread of senescence from a single primary senescent cell, senescence can spread at higher values of ND, and the secondary senescent cells are more densely distributed around the lesion (Figure S1). We found that juxtacrine senescent cells also slow the spread of senescence in this more complex setting (Figure 3), from hours to several days, for a given value of ND (Figure 3c, solid vs. dotted curves). However, once senescence starts to spread, it spreads rapidly and in an uncontrolled manner. As juxtacrine senescent cells in our model do not produce a SASP, their presence is analogous to removing trees to prevent the spread of a forest fire.3FIGUREStochastic simulation of the spread of senescence. (a) The simulation grid randomly populated with cells. (b) An illustration of the time delay in senescence induction. (c) The spread of senescence from a lesion of primary senescent cells with Rtot=105, σ=0.35, ke=koff=0.2, ND=62 and no delay in the induction of senescence. The left‐hand side plots show the spatial distribution of senescent cells and the right‐hand side shows the creation of paracrine secondary senescent cells over time, both with juxtacrine senescence (solid) and without (dotted). Vertical grey dashed lines mark days 1–7. (d) The same simulation as in (c), but with a 6‐day time delay before a cell can induce senescence in its neighbours. (e) The same simulation as in (d), but with immediate induction of juxtacrine senescence and delayed induction of paracrine senescence.A time delay in the senescence induction process limits spread and alters the dynamicsAfter a growing cell binds the required number of ligands to induce senescence, it takes a number of days before the cell becomes fully senescent and it is capable of inducing senescence in its neighbours. To investigate this effect, we introduced a 6‐day time delay between senescence induction and a cell becoming fully senescent in the stochastic simulation, as Young et al. (2009) suggest 6 days are required to establish full Ras‐induced senescence (Figure 3b).We found that introducing a time delay drastically increases the time taken for senescence to spread, as expected (Figure 3d and Figure S1b). An initial region of paracrine secondary senescent cells is created by the primary lesion, limited in diameter by the reach of the SASP from the primary senescent cells (Figure 3d). Once these cells become fully senescent and produce their own SASP, the signal spreads further, inducing a new wave of paracrine secondary senescent cells.The time delay also alters the density of juxtacrine secondary senescent cells, reducing their number at larger distances from the primary lesion (Figure S5c,d). This is due to the assumption that a cell must be fully senescent before it can induce either juxtacrine or paracrine senescence in its neighbours. The delay in juxtacrine induction increases the likelihood that neighbouring cells will become paracrine secondary senescent in the intervening time, increasing the influence of the SASP over senescence spread. The reduction in juxtacrine cells further from the senescent lesion means that as senescence spreads further, the rate of senescence spread increases due to the higher ratio of paracrine secondary senescent cells, which is inconsistent with experimental observations in vitro and in vivo (Acosta et al., 2013; Michaloglou et al., 2005).When ND, the number of ligands that a cell needs to bind to induce senescence, is decreased, the number of paracrine secondary senescent cells plateaus at a larger number for several days, before increasing sharply again once the cells become fully senescent. In an attempt to find a stable system, where the spread of senescence is limited and controlled, we increased ND. Instead of creating a small, stable number of paracrine senescent cells, this led to no senescence spread on a time scale of several years from this small primary lesion, suggesting that varying ND alone cannot lead to controlled senescence spread.A dynamic senescence induction may be required to create a stable senescent lesion. Several experiments have determined that senescence induction in a cell is a dynamic process, with evidence that NOTCH signalling (required for juxtacrine induction) is first up‐regulated before the cell starts to produce cytokines (for paracrine spread) (Hoare et al., 2016). Preliminary evidence suggests that primary cells induce juxtacrine senescence during this intermediate phase (days 3–5) of senescence activation (Teo et al., 2019). This time‐dependent senescence induction may ensure a larger number of juxtacrine secondary senescent cells, limiting the spread of senescence.A dynamic SASP is required to avoid exponential increases in cell numbers at late timesTo further investigate the effects of a dynamic SASP we simulated the limiting case, in which senescent cells can induce juxtacrine senescence in their neighbours as soon as they begin the journey towards becoming fully senescent. In reality, evidence suggests that juxtacrine induction may actually be possible after a 2‐day delay (JAG1, the NOTCH ligand responsible for juxtacrine senescence induction, is upregulated even after 1 or 2 days [Hoare et al., 2016]), but simulating the limiting case of immediate juxtacrine senescence induction clearly shows the effect of this change on the dynamics of senescence spread.The dynamic induction of senescence both slows the spread of senescence and stops the explosion in the number of senescent cells as the lesion grows, which was observed without it (Figure 3d vs. Figure 3e).Reduced SASP from secondary senescent cells can lead to controlled senescence spreadWe wanted to explore parameter regimes in which the spread of senescence is controlled (local over a human lifetime), but this was difficult with such a large number of unknown system parameters.We have previously outlined how senescence spread can be slowed with the presence of juxtacrine senescence, a time delay, and juxtacrine signalling, which precedes paracrine signalling. If the spread of senescence can be slowed drastically enough, this will lead to controlled senescence spread, as even over a human lifespan the senescence spread will be local and contained.Once senescence spread from a lesion has begun it can only stop if the balance of system parameters is right, for example, if the flux of SASP from the growing lesion is no longer enough to induce senescence. One way to achieve this is to ensure that the number of juxtacrine secondary senescent cells is large enough to dilute the SASP. This is also achieved if paracrine secondary senescent cells produce fewer SASP molecules than primary senescent cells, as suggested by experimental evidence (Acosta et al., 2013). In our simulations, we find that in this case, the primary cells in the initial lesion drive most of the senescence spread (Figure S2a). Some secondary senescent cells are generated over short time scales, but the spread eventually slows down to the point where it remains controlled, in some cases even after 500,000 hours (on the order of a human lifespan). When the lesion size is increased the size of the region of secondary cells will also increase (Figure S2b).Parameters of paracrine senescence spread can be determined through inferenceTo estimate the parameters of the model, we propose a similar experiment to those described in Acosta et al. (2013), seeding a circle of primary senescent cells and watching the phenotype spread (Section 4.3). By focusing on the first wave of senescence induction, we can forgo computationally intensive simulations. Due to the time delay before a cell becomes fully senescent, there is a period of time when a cell is unable to pass the senescent phenotype on to neighbouring cells through paracrine signalling. If it takes cells 6 days to become fully senescent, we can stop the experiment before 6 days. Using a marker to detect cells that are becoming senescent, we can determine the region to which senescence has spread through paracrine signalling (Figure 4a).4FIGURE(a) A cartoon showing the experimental setup proposed to infer the system parameters. (b) Posterior distribution of parameters from the Bayesian inference on simulated data. The plots below the diagonal show the kernel density estimation, and the plots above the diagonal show the final points used in the sampling. Red lines show the true parameter values.Focusing on short timescales allows us to build on the minimal model described in Section 4.1 to predict the fraction of senescent cells at a given distance from the primary lesion. By comparing this with the experimentally observed fraction of senescent cells we can determine NE, ND, Rtot and kon using Bayesian inference, fixing other model parameters to experimentally determined values (Section 4.3). As a proof of principle for this experiment, we used our model to generate synthetic data and attempted to infer the parameters used to generate the data. As a prior distribution of the parameters, we chose a uniform distribution over the parameter regime given in Table 1 (Figure 4b).We show that it is possible to infer some of the parameters accurately from such an experiment (Figure 4b). The posterior distribution for Rtot is narrow and centred on the true value, and although the posterior for ND is not as narrow it is well centred on the true value. This reflects the intuition we gained from the minimal model, that senescence spread is sensitive to ND and Rtot. The posterior distributions for NE and κon are much broader and less centred on the true values.DISCUSSIONCurrently, it is unclear how senescence can spread in a controlled way from a primary lesion. We systematically investigated senescence spread, using mathematical modelling informed by experimental evidence. We found that the number of receptors per cell and the binding rate of ligands were critical in determining the distance over which senescence spreads. We observed that an intermediate number of receptors led to the most effective signal spread and uncontrolled senescence spread. Several scenarios (such as juxtacrine‐induced cells, in turn, being capable of juxtacrine senescence induction) were not explored as they would lead to the uncontrolled spread of senescence in any parameter regime, highlighting that the current way of thinking about senescence spread is naive.Our results provide new insights into the role of juxtacrine senescent cells in controlling the spread of senescence. We demonstrate the crucial importance of the dynamic nature of senescence induction in generating juxtacrine secondary senescent cells before paracrine secondary senescent cells and limiting the size of senescent lesions over extended periods. As juxtacrine cells do not emit a SASP (or their SASP is blunted), the presence of juxtacrine senescent cells both reduces the concentration of SASP and limits the parameter regime in which senescence can spread, as their presence limits which cells can become paracrine senescent. Additionally, we propose an experiment aimed at determining four critical parameters involved in senescence spread, and we show that Bayesian inference is a viable tool for estimating these parameters from experimental data.Although determined from a 2D experiment, several of the inferred parameters (ND, NE, Rtot, κon) also describe paracrine senescence spread in 3D. We believe our model of senescence spread is well‐justified to study the in vitro scenario and of an appropriate complexity to identify which assumptions are important or affect senescence spread. Once parameterized with future in vitro experiments, the model can straightforwardly be expanded to 3D, to describe senescence spread in tissues, in vivo.While we can quantify parameters associated with the paracrine spread of senescence from a lesion, determining parameters associated with juxtacrine secondary senescent cells is more challenging.Single‐cell transcriptomic data may help to determine the difference between the senescent cell types. However, it is likely that differences between senescent subtypes will vary depending on the cell type and method of senescence induction, suggesting that senescence spread might differ depending on the location in the body and the cause of senescence.By considering a minimal model of senescence spread we have extracted an understanding of the dynamics of this system. Understanding the way senescence spreads through a population of cells can help us understand several biological processes, from the senescent lesion created in the radiation treatment of cancers to the accumulation of senescent cells as we age. Some questions remain unresolved, such as, whether paracrine senescence cells produce the same amount of SASP as primary senescent cells. These questions may be resolved by specially designed in vitro and in vivo experiments in combination with our modelling and inference framework.METHODSMinimal modelAutocrine and paracrine signalling in cell cultureTo understand the limiting factors in the dynamics of senescence spread, we first considered a minimal mathematical model.To mathematically describe signalling in a 2D layer of cells covered by a 3D layer of media (Figure 1c), we required knowledge of how a ligand produced by a cell will diffuse through a layer of media. Equation (1), from Batsilas et al. (2003), gives the cumulative probability of a ligand produced by a cell binding at distance r from the emitting cell.To arrive at Equation (1), Batsilas et al. assumed that ligands reflect from the upper surface of the media on contact and that ligands can bind to receptors on the surface of cells, but they reflect from the space between cells (the dish). Thus, they approximated the cell‐covered dish as a semi‐reflecting boundary layer (which is less accurate between 1<r/rcell<2, close to the emitting cell) and simulated the diffusion and binding of a large number of emitted ligands.We recreated the simulations described by Batsilas et al. and verified that Equation (1) predicts the binding location of ligands with a relative error <5%. For paracrine trajectories, we compared the time at which a simulated ligand first binds to a new cell to the square of the mean distance at which it binds (Figure S4b,c), showing that the limiting timescale in this process is diffusion (Müller & Schier, 2011). The shape of the cumulative binding probability (Pr) is shown in Figure 1d, for both small and large numbers of receptors per cell, Rtot.Equation (1) describes paracrine trajectories, not including ligands that bind to the cell that emitted them (autocrine trajectories). The probability of an autocrine trajectory is given as,2Pau=DaDa+4/π,where Da is defined in Table 2 (Batsilas et al., 2003).2TABLETable showing compound parameters from Batsilas et al. (2003).ParameterMeaningFormulaνProbability of internalisationκe/κe+κoffκRate constantκonRtotπrcell2NADaDamköhler numberrcellκDLκeffEffective trapping rate for a homogeneous surfaceκσ1+πDa/4Equations (1) and (2) do not account for the dissociation of ligands from a cell. If we assume ligands can be internalised by cells via endocytosis with rate κe and dissociate with rate κoff, the equations can be modified to describe the final distance at which ligands are bound,3Pinr≈rr+1.1DL/νκeff,andPauin=νDaνDa+4/π,where ν=κe/κe+κoff. The first binding event and internalisation are the same if κoff=0 and approach each other when the rate of endocytosis is large and the rate of dissociation is small.Between the first binding event and internalisation a ligand could bind to many cells. However, this binding may be only for a short period of time if the rate of dissociation is high. In this paper, we consider only the location of the final binding of the ligand to be important for senescence induction. This assumption breaks down if even fleeting binding of ligands to a cell produces a signal that persists intracellulary, and when combined with the signal from many other fleeting ligand binding events also induces senescence. A comparison of the binding and internalisation distances for an example set of parameters is given in Figure S6.The diffusion of ligands is clearly affected by the height of the media layer, h. However, to sustain cells in vitro, we would require a layer of media with h>2 mm. Using dimensional analysis and biologically reasonable system parameters, Batsilas et al. (2003) argue that the characteristic trajectory length, DL/κ, is much smaller than 2 mm, and so h does not need to be accounted for.Senescence spread from a single cell to an annulusExtending the work of Batsilas et al., we find that for each ligand produced, the probability of a ligand binding in an annulus of width 2rcell (illustrated in Figure 1e) at distance r from the emitting cell is,41−Pau×r+rcellr+rcell+a−r−rcellr−rcell+a=1−Pau2arcellr2−rcell2+2arcell+a2,where a=1.1DL/νKeff. The area of this annulus is given by 4πrrcell.To find the probability that a cell in this annulus will bind the ligand, we multiply this probability by the area of one cell relative to the annulus, to get,51−Pauarcell22rr2−rcell2+2arcell+a2.If one cell emits NE ligands per unit of time, the expected number of ligands bound to a cell at a distance r, per unit of time, is,6λr=1−PauNEarcell22rr2−rcell2+2arcell+a2,assuming that diffusion has reached a steady state, and NE is constant so that the rate of ligands reaching a cell at distance r is also constant. This is justified if the diffusion timescale is much faster than the other timescales in the system.To find the probability per unit time that a specific cell in an annulus will become senescent in any given second, λcellr, we assume that a cell must bind to a threshold number (ND) of ligands per unit of time to become senescent. Here, we use the binomial distribution to describe this process, where NE becomes the number of trials, and Equation (5) gives the probability of success.However, there is more than one cell in an annulus, on average there will be m=4rσrcell. Therefore, λspreadr=mλcellr is the rate of senescence induction in each annulus.Following this analysis through, we find the probability of creating at least one senescent cell at distance r from the primary lesion in time t is,7PrT≤t=1−e−λspreadrt,which is 1, minus the probability (given by a Poisson process) that no cells become senescent.Total senescence spread from a single cellWe now have the probability that one cell at distance r from the primary lesion will become senescent in a time period. We use the Poisson‐Binomial distribution to calculate the probability that a senescent cell is created in at least one annulus, as the probability of senescence induction in each annulus is different. Therefore, the probability of senescence induction over all annuli is given by a set of independent Bernoulli trials that are not identically distributed, where each annulus is a Bernoulli trial with p=PrT≤t, and we calculate the probability that one or more of these trials is successful.We use this simplified model of the system to determine whether the spread of senescence is uncontrollable, considering the limiting case of a single primary senescent cell. If a single senescent cell spreads senescence to one other cell before the body can remove it (we approximate that this happens in 2 days), with a probability of ~1, we define the spread of senescence as uncontrollable. This uncontrollable spread of senescence is not observed in vitro; therefore, this is an unrealistic parameter regime.Senescence spread with juxtacrine cellsWe considered two different regimes, with and without juxtacrine secondary senescence. Where juxtacrine signalling induces secondary senescence, we expect that cells surrounding the primary senescent cell will become juxtacrine senescent through contact. Assuming that juxtacrine secondary senescent cells do not emit a SASP, for the senescence to spread, a cell beyond this ring of juxtacrine cells must become senescent within 2 days. If juxtacrine cells instead produce a blunted SASP, this will lead to an underestimate of the spread of senescence (Rattanavirotkul et al., 2020).The distance of the annulus from the primary senescent cell which is considered in each case is given by,8r1=πrcell2σ,or,r2=2πrcell2σ.The first value, r1, is simply the average distance between cells, and r2=2r1, accounts for a ring of juxtacrine senescent cells.LimitationsThis minimal mathematical model provides a tractable way to consider senescence spread from a single cell but is difficult to extend to a system containing many cells. Furthermore, this model cannot be used to investigate the effect of a time delay in senescence induction, and it makes a number of assumptions which may limit the accuracy of the predictions.Stochastic simulationsTo deepen our understanding of the system and more closely match the reality of senescence spread in vivo, we created a stochastic model for the spread of senescence from multiple initial cells and over many subsequent senescence induction events. The stochastic model lets us introduce the added complexity of paracrine senescent cells which can in turn emit a SASP. We used the simulation to explore the change in dynamics resulting from changing the strength of the SASP produced by these secondary senescent cells. Finally, we introduce a time delay between senescence induction and SASP, more accurately modelling the cellular biology. Cells do not immediately start emitting a SASP once they have bound enough ligands to induce senescence, and instead, take several days to become fully senescent.Simulation implementationWe consider three possible cell states: senescent, in the process of becoming senescent, and non‐senescent. A cell starts to become paracrine secondary senescent when it binds to more than ND ligands per unit of time (see Section 2.1.1). The simulation is initialised on a square plane, and a 100 × 100 cell radii area is commonly used as larger areas increase the computational time taken for each simulation. Cells are then created so that they randomly populate this area with a specified density (Figure 3a). Primary senescent cells are seeded at time 0. The simulation proceeds by tracking the state of each cell.In the stochastic simulation, we approximate the number of ligands bound from each emitting senescent cell using the Poisson distribution. This approximation holds when the number of ligands emitted, NE, is large and the probability of binding (equation (5)) is small. Large binding probabilities are not of interest as they would lead to the fast and uncontrolled spread of senescence. The number of ligands bound from all senescent cells is given by a Poisson distribution whose mean is the sum of the Poisson distribution of ligands from each emitting cell (equation 6),9λtotr=∑i=1Sλir,where S is the number of senescent cells emitting a SASP.This approximation dictates the appropriate time frame for the simulation. Considering paracrine signalling using only IL‐6, IL‐8 and Activin A molecules, we can estimate NE≈2887h−1.We can now calculate the probability of a cell becoming senescent in one unit of time from its position relative to the SASP‐emitting cells. This probability can then be used to give the rate of senescence induction for that cell. Thus, each cell on the grid which is not yet senescent has a rate of senescence induction. By summing these rates we calculate the time at which the next cell becomes paracrine secondary senescent. We use the Gillespie algorithm, meaning that we first determine the time an event occurs based on creation rates, for example, a cell becomes paracrine secondary senescent at t = 1.5 h. We then determined which cell becomes senescent, again based on the rate of senescence induction associated with each cell (Gillespie, 1977).The situation is complicated further by the time delay: It may be that cells have finished becoming fully senescent before the next paracrine induction event, in which case the SASP that they produce must also be accounted for. This is achieved simply by tracking the number of cells that are in the process of becoming senescent. If a cell finishes becoming senescent before the next reaction (i.e., new senescence induction event) occurs, then the cell is added to the list of fully senescent cells, and the simulation run again from that time point, with the probabilities of each cell becoming senescent re‐calculated including any SASP from the newest senescent cell. This works as the system is Markovian.There are a number of algorithms that can be used to speed up the delayed reaction calculation (Cai, 2007). However, at any one time, we do not have a large number of cells waiting to become senescent and the bottleneck in the process is the repeated calculation of the probability of inducing senescence in each cell.The simulation code can be found on GitHub (https://github.com/lkmartin90/Senescence_Spread), and a flowchart of the stochastic simulation implementation is shown in Figure S10.Comparison of stochastic simulations to minimal modelWe compared our stochastic simulation to our minimal model to verify that the two agree. To do this, we considered the spread of senescence from a single cell, in the case where no juxtacrine secondary senescent cells are created.We used the stochastic simulation to simulate this scenario 100 times (fixing all parameters but initialising a new population of cells each time) and recorded the proportion of the simulations in which there was no senescence spread. We varied ND and Rtot, to ensure that the simulations agree over a range of parameters. We compared the results with our minimal model, which gives us the probability of senescence spread.From this, we see a good agreement between the two methods at both Rtot=105 and Rtot=106, both when cells in the stochastic simulation are created so that they lie on a lattice (Figure S3a,b), and when cells are randomly distributed (Figure S3c,d). There are some differences in the regime without juxtacrine senescence induction with Rtot=106, due to the large number of ligands which bind to nearby cells and the assumption we make when constructing the minimal model that cells are uniformly distributed (Figure S3b, dotted).In the stochastic simulation, with cells distributed randomly (Figure S3c), the range of parameters that result in a probability of senescence spread >0 and <1is greater. If a nearby cell happens to be closer than the average distance between cells then there is a higher probability that senescence will spread. This implies that the parameter regime of large Rtot is highly dependent on the initialisation of the simulation, in particular, the precise location of the cells. Although the average density in each simulation is the same, the signal spread is possible in some cases and not others. If there are no cells close enough to the primary senescent cell, senescence will not spread.We investigated whether the optimum number of receptors that we found for senescence spread in the minimal model, as shown in Figure 2, also appears in the stochastic simulation, or whether this is a consequence of the simplifying assumptions made in the minimal model. We see that the random nature of the cell placement in the stochastic simulation means that there is no longer a clear optimal number of receptors to ensure the spread of senescence, however, the general dynamics remain the same, with juxtacrine cells reducing the spread of senescence and the number of receptors altering the distance over which cells communicate (Figure S3d).LimitationsThe stochastic simulation is a useful tool to understand senescence spread from a small lesion, but computationally intensive to scale to larger spatial regions. These limitations, as well as the relationships between the models discussed in this paper, are shown in Figure S9.Parameter inferenceWe use Bayesian statistics to infer system parameters from a realistic in vitro experiment, where a circle of primary senescent cells is seeded in the centre of a dish, as described in Section 2.7. We calculate the distance senescence spreads from the circular lesion.The number of SASP molecules emitted per unit area by a primary senescent lesion is, NEσπrcell2, and the probability that a single test cell (Figure 4a) at distance x from a senescent cell binds a ligand emitted by a senescent cell is101−Pauarcell22xx2−rcell2+2arcell+a2.Integrating for the emission from the area of seeded cells to get the total expected number of ligands binding to the test cell gives111−PauNEσπrcell2arcell22∫02π∫0rseedr′xx2−rcell2+2arcell+a2dθdr′.Writing x in terms of r, r′ and θ, according to the geometry in Figure 4a, gives121−PauNEσπrcell2∫02π∫0rseed1r′2+r2−2rr′cosθ0.5×r′r′2+r2−2rr′cosθ−rcell2+2arcell+a2dθdr′.Here, rseed is the radius of the circle of primary senescent cells. We solve this integral numerically.Using the method previously described in Section 4.1, we calculate the probability that the cell that binds these ligands becomes senescent, using a Poisson distribution to determine the probability that more than ND ligands are bound per hour. We then calculate the probability that a cell becomes senescent within the duration of the experiment. Finally, we calculate the fraction of cells, which are senescent at a given distance from the seeded senescent region of cells. As this process is probabilistic, the fraction of senescent cells varies on each run of the model.Before comparing the model output with an experiment, we ought to consider whether the assumptions made in the model's construction are reasonable. In particular, in the model, we have not accounted for the time taken for ligands to diffuse away from the cell that emitted them. From Figure S4, we can see that the overall process of signal spread is diffusion‐limited. Previous experiments show senescence spread over a distance of the order of 1 mm, indicating that this is a meaningful distance to consider. The time taken for ligands to travel 1 mm is between 50 and 100 seconds for the parameter regimes shown, depending on the diffusion coefficient of the ligand. As long as the time taken for diffusion is small compared with the timescales associated with senescence spread (hours), we can assume that the diffusion occurs instantaneously.To determine the parameters that lead to the closest agreement between the model and the experimental data, we compare two summary statistics. First, the distance where the fraction of senescent cells crosses 0.5 (S). Second, the fraction of senescent cells at different distances from the lesion (F→). The distances need not be equally spaced, and we use more data points around the region of r=S, where the fraction of senescent cells is <1 and >0. We verified that these summary statistics are sensitive to changes in the parameters (Figure S8a).To quantify the agreement between predicted and observed summary statistics, use the distance function:13d=d12+d22,where d1=Smodel−Sdata and d2=∥F→model−F→data∥.We use pyABC (Klinger et al., 2018), a likelihood‐free inference software using a sequential Monte Carlo scheme, to estimate posterior distributions of model parameters. In this paper, we use synthetic data, created from the same model we use for the inference.Our inference ran on 60 cores for 13 populations, with 1000 particles, taking 7.5 h in total. While acceptance thresholds for d (epsilon values) have not converged at this stage, indicating that further iterations may improve the inference, the fraction of accepted samples was ≈10−3, making it computationally wasteful to run for more populations (Figure S8b).AUTHOR CONTRIBUTIONSL.S. and T.C. supervised the study. L.K.M., L.S. and T.C. wrote the manuscript. L.K.M. and L.S. developed the modelling framework. L.K.M. conducted the analysis.ACKNOWLEDGMENTSL.K.M is a cross‐disciplinary research fellow supported by funding from the CRUK Brain Tumour Centre of Excellence Award (C157/A27589). T.C. is supported through a Chancellor's Fellow at the University of Edinburgh and the MRC Human Genetics Unit. We thank M. Nicholson for the critical reading of the manuscript. We thank all members of the Chandra lab and Schumacher lab for their input. We thank Acosta et al. (2013) and Schafer et al. (2020) for sharing additional data related to their publications on request.CONFLICT OF INTEREST STATEMENTThe authors declare no conflict of interest.DATA AVAILABILITY STATEMENTCode used in this manuscript is available at https://github.com/lkmartin90/Senescence_Spread.REFERENCESAcosta, J. C., Banito, A., Wuestefeld, T., Georgilis, A., Janich, P., Morton, J. P., Athineos, D., Kang, T. W., Lasitschka, F., Andrulis, M., Pascual, G., Morris, K. J., Khan, S., Jin, H., Dharmalingam, G., Snijders, A. P., Carroll, T., Capper, D., Pritchard, C., … Gil, J. (2013). A complex secretory program orchestrated by the inflammasome controls paracrine senescence. Nature Cell Biology, 15(8), 978–990. https://doi.org/10.1038/ncb2784Acosta, J. C., O'Loghlen, A., Banito, A., Guijarro, M. V., Augert, A., Raguz, S., Fumagalli, M., Da Costa, M., Brown, C., Popov, N., Takatsu, Y., Melamed, J., d'Adda di Fagagna, F., Bernard, D., Hernando, E., & Gil, J. (2008). Chemokine signaling via the CXCR2 receptor reinforces senescence. Cell, 133(6), 1006–1018. https://doi.org/10.1016/j.cell.2008.03.038Baker, D. J., Childs, B. G., Durik, M., Wijers, M. E., Sieben, C. J., Zhong, J., Saltness, R., Jeganathan, K. B., Versoza, G. C., Pezeshki, A.‐M., Khazaie, K., Miller, J. D., & van Deursen, J. M. (2016). Naturally occurring p16ink4a‐positive cells shorten healthy lifespan. Nature, 530(7589), 184–189. https://doi.org/10.1038/nature16932Baker, D. J., Wijshake, T., Tchkonia, T., LeBrasseur, N. K., Childs, B. G., van de Sluis, B., Kirkland, J. L., & van Deursen, J. M. (2011). Clearance of p16ink4a‐positive senescent cells delays ageing‐associated disorders. Nature, 479(7372), 232–236. https://doi.org/10.1038/nature10600Bartkova, J., Rezaei, N., Liontos, M., Karakaidos, P., Kletsas, D., Issaeva, N., Vassiliou, L. V. F., Kolettas, E., Niforou, K., Zoumpourlis, V. C., Takaoka, M., Nakagawa, H., Tort, F., Fugger, K., Johansson, F., Sehested, M., Andersen, C. L., Dyrskjot, L., Ørntoft, T., … Gorgoulis, V. G. (2006). Oncogene‐induced senescence is part of the tumorigenesis barrier imposed by DNA damage checkpoints. Nature, 444(7119), 633–637. https://doi.org/10.1038/nature05268Basisty, N., Kale, A., Jeon, O. H., Kuehnemann, C., Payne, T., Rao, C., Holtz, A., Shah, S., Sharma, V., Ferrucci, L., Campisi, J., & Schilling, B. (2020). A proteomic atlas of senescence‐associated secretomes for aging biomarker development. PLoS Biology, 18(1), e3000599. https://doi.org/10.1371/journal.pbio.3000599Batsilas, L., Berezhkovskii, A. M., & Shvartsman, S. Y. (2003). Stochastic model of autocrine and paracrine signals in cell culture assays. Biophysical Journal, 85(6), 3659–3665. https://doi.org/10.1016/S0006‐3495(03)74783‐3Burd, C. E., Sorrentino, J. A., Clark, K. S., Darr, D. B., Krishnamurthy, J., Deal, A. M., Bardeesy, N., Castrillon, D. H., Beach, D. H., & Sharpless, N. E. (2013). Monitoring tumorigenesis and senescence in vivo with a p16ink4a‐luciferase model. Cell, 152, 340–351. https://doi.org/10.1016/j.cell.2012.12.010Cai, X. (2007). Exact stochastic simulation of coupled chemical reactions with delays. The Journal of Chemical Physics, 126(12), 124108. https://doi.org/10.1063/1.2710253Copp'e, J.‐P., Patil, C. K., Rodier, F., Sun, Y., Muñoz, D. P., Goldstein, J., Nelson, P. S., Desprez, P.‐Y., & Campisi, J. (2008). Senescence‐associated secretory phenotypes reveal cell‐nonautonomous functions of oncogenic RAS and the p53 tumor suppressor. PLoS Biology, 6(12), e301. https://doi.org/10.1371/journal.pbio.0060301Cristofalo, V. J., Volker, C., & Allen, R. G. (2003). Use of the fibroblast model in the study of cellular senescence. In Aging methods and protocols (pp. 23–52). Humana Press.d'Adda di Fagagna, F., Reaper, P. M., Clay‐Farrace, L., Fiegler, H., Carr, P., Von Zglinicki, T., Saretzki, G., Carter, N. P., & Jackson, S. P. (2003). A DNA damage checkpoint response in telomere‐initiated senescence. Nature, 426(6963), 194–198. https://doi.org/10.1038/nature02118Demaria, M., Ohtani, N., Youssef, S. A., Rodier, F., Toussaint, W., Mitchell, J. R., Laberge, R.‐M., Vijg, J., Van Steeg, H., Doll'e, M. E. T., Hoeijmakers, J. H. J., de Bruin, A., Hara, E., & Campisi, J. (2014). An essential role for senescent cells in optimal wound healing through secretion of PDGF‐AA. Developmental Cell, 31(6), 722–733. https://doi.org/10.1016/j.devcel.2014.11.012Gillespie, D. T. (1977). Exact stochastic simulation of coupled chemical reactions. Journal of Physical Chemistry, 81(25), 2340–2361. https://doi.org/10.1021/j100540a008Hayflick, L., & Moorhead, P. S. (1961). The serial cultivation of human diploid cell strains. Experimental Cell Research, 25(3), 585–621. https://doi.org/10.1016/0014‐4827(61)90192‐6Hoare, M., Ito, Y., Kang, T.‐W., Weekes, M. P., Matheson, N. J., Patten, D. A., Shetty, S., Parry, A. J., Menon, S., Salama, R., Antrobus, R., Tomimatsu, K., Howat, W., Lehner, P. J., Zender, L., & Narita, M. (2016). NOTCH1 mediates a switch between two distinct secretomes during senescence. Nature Cell Biology, 18(9), 979–992. https://doi.org/10.1038/ncb3397Ito, Y., Hoare, M., & Narita, M. (2017). Spatial and temporal control of senescence. Trends in Cell Biology, 27, 820–832. https://doi.org/10.1016/j.tcb.2017.07.004Jochems, F., Thijssen, B., De Conti, G., Jansen, R., Pogacar, Z., Groot, K., Wang, L., Schepers, A., Wang, C., Jin, H., Beijersbergen, R. L., Leite de Oliveira, R., Wessels, L. F., & Bernards, R. (2021). The cancer SENESCopedia: A delineation of cancer cell senescence. Cell Reports, 36(4). https://doi.org/10.1016/J.CELREP.2021.109441Katzir, I., Adler, M., Karin, O., Mendelsohn‐Cohen, N., Mayo, A., & Alon, U. (2021). Senescent cells and the incidence of age‐related diseases. Aging Cell, 20, e13314. https://doi.org/10.1111/acel.13314Kirschner, K., Rattanavirotkul, N., Quince, M. F., & Chandra, T. (2020). Functional heterogeneity in senescence. Biochemical Society Transactions, 48(3), 765–773. https://doi.org/10.1042/BST20190109Kiss, T., Nyúl‐Tóth, A., DelFavero, J., Balasubramanian, P., Tarantini, S., Faakye, J., Gulej, R., Ahire, C., Ungvari, A., Yabluchanskiy, A., Wiley, G., Garman, L., Ungvari, Z., & Csiszar, A. (2022). Spatial transcriptomic analysis reveals inflammatory foci defined by senescent cells in the white matter, hippocampi and cortical grey matter in the aged mouse brain. GeroScience, 2022, 1–21. https://doi.org/10.1007/S11357‐022‐00521‐7Klinger, E., Rickert, D., & Hasenauer, J. (2018). PyABC: distributed, likelihoodfree inference. Bioinformatics, 34(20), 3591–3593. https://doi.org/10.1093/bioinformatics/bty361Kuilman, T., Michaloglou, C., Vredeveld, L. C. W., Douma, S., van Doorn, R., Desmet, C. J., Aarden, L. A., Mooi, W. J., & Peeper, D. S. (2008). Oncogene‐induced senescence relayed by an interleukin‐dependent inflammatory network. Cell, 133(6), 1019–1031. https://doi.org/10.1016/j.cell.2008.03.039Lee, P. J., Benz, C. C., Blood, P., Börner, K., Campisi, J., Chen, F., Daldrup‐Link, H., De Jager, P., Ding, L., Duncan, F. E., Eickelberg, O., Fan, R., Finkel, T., Furman, D., Garovic, V., Gehlenborg, N., Glass, C., Heckenbach, I., Joseph, Z.‐B., … Yale, T. M. C. (2022). NIH SenNet consortium to map senescent cells throughout the human lifespan to understand physiological health. Nature Aging, 2(12), 1090–1100. https://doi.org/10.1038/s43587‐022‐00326‐5Michaloglou, C., Vredeveld, L. C. W., Soengas, M. S., Denoyelle, C., Kuilman, T., van der Horst, C. M. A. M., Majoor, D. M., Shay, J. W., Mooi, W. J., & Peeper, D. S. (2005). BRAFE600‐associated senescence‐like cell cycle arrest of human naevi. Nature, 436(7051), 720–724. https://doi.org/10.1038/nature03890Müller, P., & Schier, A. F. (2011). Extracellular movement of signaling. Molecules, 21, 145–158. https://doi.org/10.1016/j.devcel.2011.06.001Muñoz‐Espín, D., & Serrano, M. (2014). Cellular senescence: From physiology to pathology. Nature Reviews Molecular Cell Biology, 15(7), 482–496. https://doi.org/10.1038/nrm3823Narita, M., Nuñez, S., Heard, E., Narita, M., Lin, A. W., Hearn, S. A., Spector, D. L., Hannon, G. J., & Lowe, S. W. (2003). Rb‐mediated heterochromatin formation and silencing of e2f target genes during cellular senescence. Cell, 113(6), 703–716. https://doi.org/10.1016/s00928674(03)00401‐xPassos, J. F., Simillion, C., Hallinan, J., Wipat, A., & von Zglinicki, T. (2009). Cellular senescence: unravelling complexity. Age, 31(4), 353–363. https://doi.org/10.1007/s11357‐009‐9108‐1Prata, L. G. P. L., Ovsyannikova, I. G., Tchkonia, T., & Kirkland, J. L. (2018). Senescent cell clearance by the immune system: Emerging therapeutic opportunities. Seminars in Immunology, 40, 101275. https://doi.org/10.1016/j.smim.2019.04.003Putavet, D. A., & de Keizer, P. L. J. (2021). Residual disease in glioma recurrence: A dangerous liaison with senescence. Cancers, 13(7). https://doi.org/10.3390/CANCERS13071560Rattanavirotkul, N., Kirschner, K., & Chandra, T. (2020). Induction and transmission of oncogene‐induced senescence. Cellular and Molecular Life Sciences, 78, 843–852. https://doi.org/10.1007/s00018‐020‐03638‐0Schafer, M. J., Zhang, X., Kumar, A., Atkinson, E. J., Zhu, Y., Jachim, S., Mazula, D. L., Brown, A. K., Berning, M., Aversa, Z., Kotajarvi, B., Bruce, C. J., Greason, K. L., Suri, R. M., Tracy, R. P., Cummings, S. R., White, T. A., & LeBrasseur, N. K. (2020). The senescence‐associated secretome as an indicator of age and medical risk. JCI Insight, 5(12). https://doi.org/10.1172/jci.insight.133668Serrano, M., Lin, A. W., McCurrach, M. E., Beach, D., & Lowe, S. W. (1997). Oncogenic ras provokes premature cell senescence associated with accumulation of p53 and p16INK4a. Cell, 88(5), 593–602. https://doi.org/10.1016/S0092‐8674(00)81902‐9Shvartsman, S. Y., Steven Wiley, H., Deen, W. M., & Lauffenburger, D. A. (2001). Spatial range of autocrine signaling: Modeling and computational analysis. Biophysical Journal, 81(4), 1854–1867. https://doi.org/10.1016/S0006‐3495(01)75837‐7Teo, Y. V., Rattanavirotkul, N., Olova, N., Kirschner, K., Neretti, N., & Correspondence, C. (2019). Notch signaling mediates secondary senescence. Cell Reports, 27, 997–1007. https://doi.org/10.1016/j.celrep.2019.03.104van Deursen, J. M. (2014). The role of senescent cells in ageing. Nature, 509(7501), 439–446. https://doi.org/10.1038/nature13193Wang, B., Kohli, J., & Demaria, M. (2020). Senescent cells in cancer therapy: friends or foes? Trends in Cancer, 6, 838–857. https://doi.org/10.1016/j.trecan.2020.05.004Wiley, H. S., Shvartsman, S. Y., & Lauffenburger, D. A. (2003). Computational modeling of the EGF‐receptor system: A paradigm for systems biology. Trends in Cell Biology, 13(1), 43–50. https://doi.org/10.1016/S0962‐8924(02)00009‐0Young, A. R. J., Narita, M., Ferreira, M., Kirschner, K., Sadaie, M., Darot, J. F. J., Tavaré, S., Arakawa, S., Shimizu, S., Watt, F. M., & Narita, M. (2009). Autophagy mediates the mitotic senescence transition. Genes & Development, 23(7), 798–803. https://doi.org/10.1101/gad.519709

Journal

Aging CellWiley

Published: Aug 1, 2023

Keywords: cell signalling; mathematical model; SASP; senescence

There are no references for this article.