Estimation of the cancer risk induced by therapies targeting stem cell replication and treatment recommendations

Rejuvenation of stem cell activity might increase life expectancy by prolonging functionality of organs. Higher stem cell replication rates also bear the risk of cancer. The extent of this risk is not known. While it is difficult to evaluate this cancer risk in experiments, it can be estimated using a mathematical model for tissue homeostasis by stem cell replication and associated cancer risk. The model recapitulates the observation that treatments targeting stem cell replication can induce a substantial delay of organ failure. The model predicts that the cancer risk is minor under particular conditions. It depends on the assumed implications for cell damage repair during treatment. The benefit of rejuvenation therapy and its impact on cancer risk depend on the biological age at the time of treatment and on the overall cell turnover rate of the organs. Different organs have to be considered separately in the planning of systemic treatments. In recent years, the transfer of blood from young to old individuals was shown to bear the potential of rejuvenation of stem cell activity. In this context, the model predicts that the treatment schedule is critical for success and that schedules successful in animal experiments are not transferable to humans. Guidelines for successful protocols are proposed. The model presented here may be used as a guidance for the development of stem cell rejuvenation treatment protocols and the identification of critical parameters for cancer risk.

predictive power of the mathematical model it has to be constructed on the basis of rather generic assumptions. The model complexity is chosen such that the driving question, whether cancer risk is under control during rejuvenation therapy or not, can be addressed, avoiding speculations on molecular mechanisms. The model (see Methods: Model equations and analytical solution) distinguishes stem (S) and tissue cells (T). It further distinguishes two reasons for death: organ failure and cancer. More specifically, the model is based on the following set of assumptions: • Tissue cells (T) have a natural life time (death rate δ) and are replenished by stem cell (S) division (proliferation rate p). • The proliferation rate p is either constant or age-dependent according to Eq. (15).
• During treatment, the division rate is changed and reverted back to its normal value at the end of treatment sessions. • The total population of stem cells (σ 0 ) remains constant, ignoring both, stem cell loss and homeostatic symmetric division. • The ratio of stem to tissue cells is 1:1000.
• The life time (Θ) is 80 years in humans and 2 years in mice.
• Death is induced by organ failure, which happens when the total number of tissue cells reaches a fraction (f) of 50%. • Both, tissue and stem cells accumulate damage with rate γ. The kind of damage is not distinguished and only the number of accumulated damage events is considered. • Damaged stem cells pass their damage state to tissue progeny in the course of homeostatic division.
• Cancer is induced at a critical number of damage events per cell (c).
• Cancer risk at the age of 80 is ξ = 30%.
An important source of cell damage is mutations upon cell division. The rate of point mutations in human adult stem cells was measured in the range of 40 mutations per year 11 , where an approximately linear relationship of the number of mutations with age was found. This supports a constant damage rate (γ). While the number of stem cell mutations correlates with cancer incidence 12 , this absolute frequency of mutations is still difficult to relate to cancer risk because not all of these mutations happen in oncogenic regions of the genome and induce cell damage. Furthermore, the immune system is removing pathogenic cells, which makes the absolute number of mutations associated with cancer unclear. In order to keep the model simple, a linear accumulation of damage with age was used initially. However, non-linear age-and division-dependent damage rates were also considered. The cancer-inducing number of cell damage events was fixed to c = 10 per cell. This choice can be considered as a worst case scenario for the cancer risk induced by rejuvenation therapy because, on the one hand, cancer is expected to be induced by more than 10 accumulated damage events, and, on the other hand, with > c 10 the cancer risk induced by the therapy is getting smaller in the model, as was revealed by a robustness analysis (see Methods: Robustness of the statements).
As the division (p) and the cell damage rate (γ) are determined by the life expectancy (Θ; see Methods: Determination of the division rate) and the cancer risk (ξ; see Methods: Determination of the damage rate) at this age, respectively, this model has only a single free parameter: The cell death rate δ, which controls the overall turnover rate of tissue cells. The turnover rate is a tissue specific property and different values are considered in the following. Parameter sets for low, intermediate, and high turnover are defined in Table 1 and used as reference parameter sets.
Starting from a young individual without any damaged cell, the time course of cell damage over age is depicted in Fig. 1. In this setting, the age of death is determined by the time, at which the total number of cells constituting the virtual organ reaches the fraction f = 50%. This time will be denoted as age of organ failure in the following. The age of cancer onset is determined by the time when the sum of all cells with c = 10 or more damage events exceeds the limit of 1 cell (purple lines in Fig. 1). Assuming ergodicity, it can be stated that a number of cells with 10 or more damage events below 1 reflects the cancer risk at a particular age. In the following, the age of cancer risk ξ will be discussed, which denotes the age at which the number of cells with c = 10 or more damage reaches Organ failure is more sensitive to constant stem cell replication rates than cancer risk. For life long fixed stem cell replication and damage rates, the model can be solved analytically (Eq. (2)). In order to determine the sensitivity of organ failure and cancer risk to changes in the division rate, the constant stem cell replication rate p is varied. Figure 2 shows how the ages of organ failure Eq. (14) and cancer risk ξ Eq. (12) depend on the assumed division rate. The normal division rate determined in Eq. (6) (vertical dashed line), by construction, induces organ failure and a cancer risk of ξ = 0.3 at the age of Θ = 80 years. The division rate for perfect tissue homeostasis determined in Eq. (11) prohibits death by organ failure. The result shows that the age of organ failure is highly sensitive to the division rate (red line in Fig. 2) while the age of cancer risk ξ is hardly changing (blue line in Fig. 2). For a stem cell replication rate inducing organ failure at the age of 400 (this is achieved by doubling the stem cell replication rate), the age at which the cancer risk reaches 30% is shifted from 80 to 78. Note that this result relies on life long fixed stem cell replication and damage rates.
Organ failure is delayed with limited cancer risk. Rejuvenation  a duration of 10 years starting at different ages (see Methods: Treatment at particular ages), a similar picture in terms of sensitivity of the ages of organ failure and cancer risk emerges (full lines in Fig. 3). The overall effect on the age of organ failures is smaller than for life long treatments (Fig. 3a,b). With an induced division rate of perfect tissue homeostasis Eq. (11), organ failure happens at ages between 85 and 100 years, and the cancer risk of ξ = 30% is reached at age 79 ± 1 years. When the same treatment is repeated three times every 10 years with pausing intervals of 10 years (dotted lines, Fig. 3a,b), organ failure happens at an age between 110 and 150 years and the cancer risk of 30% is reached at the age of 78 ± 1 years. While the absolute numbers have to be interpreted carefully, the difference in the sensitivity of both read-outs is a stable finding.
Organ failure is delayed the more the later the treatment is started. Cancer risk is increased correspondingly. Surprisingly, when the first of the three treatments is started at an age of 42 or later, cancer risk is getting lower when the treatment is started later (see Fig. 3a,b; blue and shorter wavelength curves). In this regime, rejuvenation treatment has a strong effect on the delay of organ failure with a comparably small increased risk of cancer. However, overall organ functionality is lower when treatment is started later. The treated individual faces a  physiological state that is kept above but close to the organ failure threshold, which implies a lower life quality compared to treatments started at younger age.
High-turnover organs are more sensitive to rejuvenation treatment but less stable. The investigation was built on organs with a low turnover rate (see Table 1, low). For high turnover organs (Fig. 3c,d), the sensitivity of the age of organ failure to a single rejuvenation treatment is even higher. With an induced division rate of perfect tissue homeostasis Eq. (11), organ failure happens at an age between 90 and 160 years, while the age of cancer risk ξ is kept between 76 and 80 years. With three treatments of 10 years, the age of organ failure can be delayed to ages of 130 to 200 years, while the cancer risk is not further increased. Importantly, for high turnover organs, a single rejuvenation treatment can delay the age of organ failure to 120 years, while inducing no major change in the cancer risk (Fig. 3c,d; full green and longer wavelength lines).
Failure of high turnover organs cannot be further delayed by higher than perfect homeostasis p juvenile stem cell replication rates. However, cancer risk enters a non-linear regime at higher replication rates (Fig. 3d). Organs with higher turnover rates react with higher or lower cancer risks than for low turnover organs depending on how and when the treatment is applied. There is hardly any increased cancer risk when the treatment is done at young age, irrespective of the target replication rate of the treatment. There is a transition to a strong impact on cancer risk at the age of 35 to 42. Cancer risk is highest for treatments started at ages around 60 and decreases again at higher ages ( Fig. 4a,b).
As the treatment is applied systemically, it will always also affect high turnover organs. The presented results illustrate that the timing and the dosage have to be planned quantitatively in order to avoid adverse effects.

Limited cancer risk also in a model with age-dependent loss of stem-cell replication. A general
tendency of rejuvenation treatment is that the later the treatment is started, the more organ failure is delayed and the more cancer risk is increased (Fig. 3). However, this result might rely on the assumption that stem cell replication is constant over the life time but for the periods of treatment. While this assumption is a good approximation for skin tissue, in which stem cell replication is widely conserved at all ages 13 , other organs exhibit a decrease of the stem cell replication rate with age 1 . For example, hematopoietic stem cells lose power in the course of ageing 14,15 . The degree of this reduction is again organ specific. Assuming a Hill-function for the dependence of p on age (see Methods: Age-dependent stem cell division rate), the set of differential equations Eq. (1) cannot be solved analytically anymore and is solved numerically.
Induction of higher stem cell replication rates in low turnover organs exhibits a qualitatively similar result as with constant p (Fig. 5a,b). An increase of the replication rate to an absolute level has a stronger effect the later in life the treatment is done. The sensitivity of the age of organ failure to treatment is still substantially higher than the sensitivity of the age of cancer risk ξ.
The mechanism of rejuvenation therapy impacts on the hierarchy of best treatment age.
Treatments increasing stem cell activity at a particular age can either relatively improve the remaining stem cell replication activity or induce an absolute replication level. In the particular case of young blood transfer, the mechanism of improvements in stem cell replication is not known. The aforementioned results (Fig. 5a,b), which  14)) is highly sensitive to changes in the cell renewal rate (p) while the age of cancer risk ξ (blue line, Eq. (12)) remains comparably stable. Unmodified stem cell renewal p normal (see Eq. (6)) and perfect cell renewal p juvenile (see Eq. (11)) are marked by the vertical dashed and full lines, respectively. Between the red and the vertical black line, Eq. (14) has no solution, which means that organ failure never happens. Parameter set for low turnover in Table 1 were based on absolute renewal rates, are compared to results from a model with relative improvements of the renewal rate.
When a rejuvenation treatment induces an improved replication relative to the level of stem cell replication at the time of treatment, the situation changes in organs with low turnover rates (Fig. 5c,d). First, the overall impact on the ages of organ failure and cancer risk ξ is smaller. Secondly, the advantageous age of treatment is different: stem cell replication increased by a factor of ten at young age induces a substantial retardation of organ failure. However, the later the treatment is done, the weaker the effect, thus, inverting the result found with either constant stem cell replication (Fig. 3) or age-dependent stem cell replication treated to adopt absolute replication levels (Fig. 5a,b). 10 times no replication is still a weak replication. At advanced age, a treatment with a relative improvement of stem cell replication has weak effects on both, the age of organ failure (Fig. 5c) and the shift of cancer risk (Fig. 4c).
For organs with intermediate or high turnover (Supplementary Figure 3), the original sequence as found in (Figs 3 and 5a,b) is restored, i.e. the later the treatment, the more organ failure is retarded. Cancer risk is highest at ages around 60 years as it was for age-independent stem-cell replication (compare Fig. 4b,d). Whether treatment of young or old individuals is more advantageous in terms of retardation of organ failure depends on the organ turnover level with a switch-like behaviour. Hence, a turnover rate may exist, at which the treatment effect becomes widely independent of the age of treatment. This is, indeed, the case (Supplementary Figure 4) and holds true in all subsequently described model variants (Supplementary Figures 15, 27 and 39). This result suggests (2) together with 10 year treatments (Methods: Treatment at particular ages) is analysed for ages of organ failure (a,c) and cancer risk ξ (b,d) in dependence on the age of first treatment (colors) and with different induced renewal rate p (horizontal axes). A single treatment (full lines, circles) is compared to three treatments with interruptions of 10 years (dashed lines, crosses). Vertical lines denote the renewal rate of humans with a life expectancy of 80 years (dashed black line) and with perfect tissue homeostasis (full black line). Parameter set of low (a,b) and high (c,d) turnover in Table 1 that, provided a treatment induces a relative increase of stem cell activity, the best age for treatment is different for different organs.
Increasing cell damage rate with age reduces cancer risk. Similar to the division rate, the damage rate is considered age-dependent. Implicitly assuming weaker repair mechanisms at higher age, the damage rate is increased with age (see Eq. (16)). The cancer risk induced by treatments increasing stem cell replication is predicted even smaller than with age-dependent division rates alone (see Supplementary Figures 9 to 20). The impact on organ failure is not altered (Supplementary Figures 11a,c, 13a,c and 14a,c).
Division-dependent cell damage may induce high cancer risk during treatment. Next, the cell damage rate was envisaged dependent on the division frequency of stem cells according to Eq. (17). This corresponds to the view, that cell mutations associated with cell division would be the dominant cue of cancer pathogenesis. In this setting, the amount of damage per cell might be identified with the number of oncogenic mutations.
In a model with age-and division-dependent damage rates, treatments increasing stem cell replication lead to a delay of organ failure similar to constant damage (Supplementary Figures 23a,c, 25a,c and 26a,c). The impact on cancer risk highly depends on how the increased division rate impacts on the cell damage rate during treatment. If it is linearly translated into a higher damage rate, the cancer risk dramatically increases (Fig. 6a). Repeated sessions of strong treatment (factors of 4 and beyond) did not further increase the cancer risk (Fig. 6a) (dotted lines) despite a substantial delay in organ failure (Supplementary Figure 23a,c). At young ages between 10 and 20, the Optimal treatment age to minimise cancer risk. Cancer risk of single 10 year treatments depending on the age of treatment (horizontal axes). Constant (a,b) and age-dependent division rate (c,d) are compared for low (a,c) and high (b,d) turnover organs (see Table 1). The colors correspond to the strength of treatment as the factor of the division rate p. Same data as in Fig. 3 a,b ( cancer risk is increased most and gets milder with increasing age of treatment (Fig. 6c). This cancer risk analysis also holds true for higher turnover organs (Supplementary Figures 22b, 25b,d, and 26b,d).
If the cell damage rate during treatment is derived from the division rate associated with the age of the treated in silico individual, the delay of organ failure remains as effective as with constant damage rates ( Supplementary  Figures 35a,c, 37a,c and 38a,c). The cancer risk is hardly increased (Fig. 6b) relative to a division-independent damage rate. The highest cancer risk is found for individuals at young age (Fig. 6d). In contrast, for high turnover organs, the worst age of treatment is around the age of 60 years (Supplementary Figure 34b). This basically restores the situation as found for age-dependent division and constant damage rate (Fig. 4c,d).
Taking these results together, most of the treatment induced cancer risk is related to the treatment phase, because the simulations were identical during the non-treatment phases. The impact of treatments increasing stem cell division and tissue regeneration onto cell damage and repair is critical to ascertain the harmlessness of a rejuvenation treatment. Furthermore, it can flip the most advantageous age of treatment from young to old individuals.
Not all rejuvenation treatment reported in mice can be transfered to men. In C57BL/6 mice with a mutation associated with the development of Alzheimer's disease, parabiosis of young and old mice for 5 weeks is analysed for ages of organ failure (a,c) and cancer risk ξ (b,d) in dependence on the age of first treatment (colors) and on the strength of treatment (horizontal axes). Absolute (a,b) or relative improvements (c,d) of the renewal rate are distinguished. A single treatment (full lines, circles) is compared to three treatments with interruptions of 10 years (dashed lines, crosses). Vertical lines denote the renewal rate p juvenile inducing perfect tissue homeostasis (full black line). Parameter of low turnover in Table 1 with age-dependent stem-cell replication rate p in Eq. (15). For higher turnover rates compare Supplementary Figure 2  per year) such that organ failure can happen at the age of two years (Table 1). Judging from the wound healing dynamics, the turnover rate of mouse organs might really be higher than in humans. However, the imposed higher turnover rate has to be interpreted as a phenomenological parameter recollecting all factors associated with faster ageing of mice the model does not contain. In the model, the only way of making mice age faster is a higher turnover rate.
In the model with age-dependent loss of stem cell division, a 5 weeks rejuvenation treatment induces a substantial effect on organ failure (Fig. 7a). In contrast, a 5 weeks treatment of humans had only weak effects even if applied to high turnover organs (Fig. 7b). This remains true when applied to a human organ with even higher and unphysiological turnover rates. This implies that the positive results found in preclinical trials in mice 3 cannot be expected in humans when the same treatment protocol is applied.
The next question was under which conditions an effect of the treatment is expected to be detectable in humans. Assuming a lifelong constant stem cell replication in high turnover organs, a 5 weeks rejuvenation treatment (inducing p juvenile in Eq. (11)) leads to a substantial improvement (Fig. 7c). This implies that organs with (1) together with 10 year treatments (Methods: Treatment at particular ages) is analysed for cancer risk ξ in dependence on the strength of treatment by relative increase of the division rate (a,b; colors reflect the age of first treatment), and on the age of first treatment (c,d; colors reflect the strength of treatment). (a,c) The damage rate adapts to the increased division rate during treatment. (b,d) The damage rate is calculated based on the age-specific division and damage rates, also during treatment. Parameter of low turnover in Table 1 with age-dependent stem-cell replication rate p in Eq. (15). For absolute renewal rates compare Supplementary Figures 23 and 35, panels a,b in each. For higher turnover rates: Supplementary Figures 25 and 37  high turnover and with slow ageing like the skin would benefit from short therapies. In contrast, low turnover organs with ageing stem cells like the brain will only exhibit an effect comparable to the one seen in mice when the treatment is prolonged to the scale of years (Fig. 7d). Both statement remain true, if a relative instead of an absolute improvement of stem cell replication is assumed (Supplementary Figure 5). Consistent with the results in Fig. 5, the most advantageous age of treatment in terms of cancer risk is altered in this case. Both statements also remain true, when age-dependent or age-and division-dependent replication rates are assumed ( Supplementary  Figures 12 and 24 for absolute, and Supplementary Figures 16 and 28 for relative improvements of stem cell replication).

Discussion
The mathematical analysis of the impact of rejuvenation treatments improving stem cell replication and tissue regeneration predicts that organ failure is delayed at the cost of an increased risk of developing cancer. The degree of increased cancer risk depends on how the treatment modulates cell damage rates. The details of improved life time depend on the age of treatment, the turnover rate of the organ, and the mechanism of action of rejuvenation treatment.
The models with age-dependent and age-independent stem cell replication rate have many properties in common. The larger the turnover rate of the considered organ, the more efficiently organ failure is retarded by small improvements in stem cell replication. However, organ failure cannot be further delayed by stronger treatment, i.e. there is a natural limit of retardation of organ failure. In organs with higher turnover rates, the cancer risk is stable in response to moderate treatments performed at early age. This opens an opportunity to really delay organ failure without major risk to induce cancer. However, a conservative dosage of the treatment and a subtle analysis of the biological age of the individual are required. At some biological ages, the treatment can induce an unwanted increase of cancer risk.
Limited cancer risk is not found when the damage rate relied directly on the division rate. In this case, a substantially higher cancer risk is found irrespective of the turnover rate of the organs and the highest risk of cancer is found for treatments at young instead of old ages. A tolerable increase of cancer risk from 80 to 70 years of age is found for a doubling of the renewal rate, limiting possible treatment intensities. This situation corresponds to cancer pathogenesis being exclusively driven by the accumulation of mutations during stem cell division. As mutations by division are an important cause of cancer development, this is a relevant limitation. However, the increased cancer risk also relies on the assumption that the rejuvenation treatment does not alter the efficiency of DNA repair mechanisms. If the treatment not only increases tissue regeneration by stem cell replication but also rejuvenates the efficiency of repair mechanisms, cancer risk is limited as for division-independent damage rates. This implies, that for each therapy targeting stem cell activity, understanding its impact on repair mechanisms is pivotal for the estimation of the associated cancer risk.
Further, the treatment strategy has to consider the different turnover rates of the different organs. If the treatment is provided systemically, as in the case of young blood transfer, all organs are equally affected. The optimum of treatment dosage and treatment age, however, turned out to depend on the organ turnover rate. As a consequence, the treatment has to be chosen in a rather conservative way such that none of the organs bears an increased cancer risk, but as many of the organs as possible benefit from delayed organ failure.
It is important to note that the model predictions rely on averaged human data and a rather simple model of cancer pathogenesis. It only distinguishes cancer types developed from environmental damage or from division-associated mutations. The reality of cancer pathogenesis is rather variable. Different stem cell populations might be considered in future model extensions. In the current model, all results apply to the mean of these potential paths to cancer. An individual at particular risk for a specific type of cancer might be at higher or lower risk for developing cancer in response to rejuvenation treatment than estimated by the model. Therefore, the model can be used to estimate the impact of treatments on cancer risk on a population level but not as a tool for individualized therapy recommendations.
The model suggests that the best possible treatment strategy also depends on the mechanistic details of rejuvenation treatment. If the treatment induces an absolute replication rate of stem cells, the age-dependence is similar for all organ turnover rates. However, if the treatment induces an improvement of the replication rate relative to the replication potential of the stem cells at the age of treatment, treatment of older individuals can be without major effect on organ failure in low turnover organs. Furthermore, it is not clear whether treatments improve stem cell activity only during the treatment or keep the increased activity after treatment. The estimated effects would be enhanced if a treatment memory existed. Therefore, there is an urgent need to clarify the mechanistic details of the effects observed upon rejuvenation therapy, which, in particular, applies to young blood treatments 16 .
In the context of young blood rejuvenation therapy, improved tissue homeostasis was found in a mouse parabiosis model in response to 5-weeks treatments 4,5 . If the simple transfer of blood had an impact on tissue homeostasis, which is controversial 6,9 , this would open the possibility to apply young blood therapy to humans. According to the model results, the treatment schedule identified in mice with positive effect on tissue homeostasis cannot be transferred to the human case and would not induce the same degree of rejuvenation, in particular, in low turnover organs like the brain. Longer treatments in the range of years are required, which shall be considered in the interpretation of the outcome of ongoing clinical trials. A negative outcome would not necessarily imply that rejuvenation of human individuals is not possible by blood derived systemic factors. The code developed here may be used to plan clinical trials and to adapt the treatment schedule to the individual state of organs.
The strength of the presented analysis relies on the simplicity of the model based exclusively on generic assumptions and parameters determined by plausible side conditions. The analysis suggests that treatments targeting stem cell activity have a realistic chance to increase quality of life in the future. However, the treatment exhibits an age-and organ turnover rate-dependent effect and requires a decent planning. In particular, better knowledge on the induced rejuvenation mechanisms and the effect on DNA repair mechanisms are critical to avoid adverse effects. with initial values σ = S t ( )

Methods
. p is the stem cell replication rate, γ is the cell damage rate, and δ is the death rate of tissue cells, which corresponds to the turnover rate of the considered organ.
The analytical solution of the model Eq. (1)  Determination of the division rate. The division rate p can be fixed by the condition that organ failure happens at the age of Θ. Organ failure is assumed to happen when the initial number of tissue cells drops below a fraction f = 0.5. It is assumed independent of the number of damage events like mutations, thus, Eq. (1) can be solved in the limit of no damage, i.e. γ = 0, which reduces all equations to those indexed with 0. The differential equations for stem and tissue cells read: An implicit equation for γ of this type is solved by a Lambert W-function as This solution is just an upper bound for γ because all contributions of cells with more than c damage events were ignored. The condition Eq. (7)  Comparison of mouse and human. The life expectancy of mice is in the range of two years 20 . With the organ turnover rates δ for humans (see Table 1), organ failure at the age of two years cannot be induced in the model (see Eq. (13)). Substantially higher turnover rates with δ > 0.4 per year are required to get a solution of the model with positive division rates p. At such high turnover rates, organ failure at the age of human life expectancy Θ = 80 years is not achieved (see Eqs (6) and (13)).
Robustness of the statements. The sensitivity of the results to changes in the parameter values used in Table 1 and in Eq. (15) is analysed. The number of damage events per cell associated with induction of cancer c was chosen arbitrarily as c = 10. For very low c the system gets unstable. For larger c the behaviour and all described results remain qualitatively unchanged. There is strictly no impact of c on the age of organ failure. The impact of treatment on the age of cancer risk ξ is further reduced for larger values of c in all tested models (Supplementary Figures 6, 17, 29 and 41). The results based on c = 10 can be considered as an upper limit for the increase in cancer risk, as the number of damage events inducing cancer is likely not smaller than 10 but for extreme cases.
Larger steepness of the age-dependent division rate n p in Eq. (15) reduces while smaller n p increases the overall effect on both, the ages of organ failure and cancer risk ξ, in particular, for late treatments (see Supplementary  Figure 7). This dependence on n is similar in all settings (Supplementary Figures 18, 30 and 42). This is consistent with the observation that in the extreme case of constant p the effect was larger as well (see Fig. 3). A determination of the stem cell activity in human organs over life would allow for a more concrete attribution of therapy effects in silico to human organs.
A corresponding analysis of the steepness of damage rate dependence on age n γ with values of 1.5 and 3 shows a mild impact on the quantitative results (Supplementary Figures 19, 31 and 43). The same holds true for a variation of the age of half-maximum damage rate K γ = 50 years, which was shifted to 30 and 70 years ( Supplementary  Figures 20, 32 and 44). Qualitative results are not affected.
A change in the relative proportion of stem cells σ 0 and tissue cells τ 0 has strictly no influence on the results (see Supplementary Figure 8). The stem cell division rate, by construction, compensates for a different fraction of stem cells. However, the association of the virtual turnover rates with real tissues turnover rates changes. Proof of the analytical solution. The ordinary differential equation system Eq. (1) can be solved analytically (see Eq. (2)). The solution can be proven by complete induction, which is sketched here. For the stem cell compartment S i with ≥ i 0 the base case is explicitly calculated. The solutions for S 0,1 are trivial. Let us assume that the general solution Eq. (2)  A differential equation of the form is solved by