Introduction

Short- and long-term biologically based models

Biologically motivated mathematical modeling of background and ionizing radiation-induced carcinogenesis has a history spanning more than 50 years (Nordling 1953). Many of the models can be characterized as short-term, in that they focus on processes occurring during and shortly after irradiation (Hahnfeldt and Hlatky 1998; Radivoyevitch et al. 2001; Mebust et al. 2002; Schollnberger et al. 2002; Sachs and Brenner 2005; Hofmann et al. 2006; Shuryak et al. 2006; Little 2007; Sachs et al. 2007; Schneider and Walsh 2008). The main advantage of such models is that they provide a detailed initial dose–response relation for short-term endpoints, which are used as surrogates for carcinogenesis. The main disadvantage is that the possibly substantial modulations of the magnitude and shape of this initial dose–response during the lengthy period (multiple years-decades) between irradiation and manifestation of typical solid tumors are not considered mechanistically.

In contrast, another class of biologically motivated models can be characterized as long-term, in the sense that they track carcinogenesis rates throughout the entire life span, e.g., the Armitage-Doll model (Armitage and Doll 1954; Armitage 1985), the Moolgavkar-Venzon-Knudson two-stage clonal expansion (TSCE) model (Moolgavkar 1978, 1980; Moolgavkar and Knudson 1981), the two-stage logistic model (Sachs et al. 2005), and many others (Yakovlev and Polig 1996; Pierce and Mendelsohn 1999; Wheldon et al. 2000; Little and Wright 2003; Pierce and Vaeth 2003; Ritter et al. 2003; Little and Li 2007). The main advantage of long-term models is the more detailed treatment of slow carcinogenesis processes, including the modulation of the radiation dose–response during the long latency period. The main disadvantage is that the initial dose–response is typically treated in a non-mechanistic, phenomenological manner.

The need for a new approach: integration of long- and short-term models

The lack of detailed treatment of radiation-specific effects typically limits risk predictions from long-term models to exposure conditions where a simple dose–response relationship holds. Exposures where this relationship is more complex, such as high fractionated radiotherapeutic doses that can lead to treatment-induced cancers in nearby organs (Travis et al. 1996, 1997, 2002, 2003; van Leeuwen et al. 2003; Travis et al. 2005) are difficult to describe with current long-term models. Conversely, the more detailed dose–responses produced by sophisticated short-term models can be scaled to cancer risks only by considering the effects of factors such as background risks, age at exposure, and time since exposure, which are not explicitly taken into account by the short-term formalisms. A new, unified approach of integrating short- and long-term techniques is needed, where a detailed initial dose–response for pre-malignant cell numbers is produced over a wide range of doses, and changes to the shape of this dose–response during the latency period before the development of cancer are also analyzed in detail. A schematic representation of model unification is provided in Fig. 1.

Fig. 1
figure 1

A general scheme of short- and long-term processes governing the total number of pre-malignant cells. Details are discussed throughout the main text. As the individual ages, the number of viable pre-malignant cells grows, but the curve may turn over and decrease at very old age (blue line in the main graph). This pattern parallels background cancer incidence, since cancer risk is assumed to be proportional to the number of pre-malignant cells (discussed in the text). Radiation exposure (e.g., radiotherapy for an existing cancer) initially causes the number of pre-malignant cells to decrease due to cell killing (red line). After exposure stops, the irradiated tissues recover, allowing pre-malignant cells to repopulate and reach a number higher than was present before irradiation, i.e., a net excess radiogenic cancer risk is produced. Fluctuations in the number of pre-malignant cells throughout the irradiation and recovery periods (i.e., the short-term processes) are shown in the inset graph

Motivations for the specific unification used here

Here, we provide a specific example of integrating short- and long-term radiation carcinogenesis models. In that a variety of short- and long-term models have been proposed, our approach should be considered as illustrative, intended to investigate the practicality of integrating the two model types. The goal was to produce a novel formalism, which can describe typical patterns of background and radiogenic carcinogenesis with the smallest possible number of adjustable parameters. This goal required multiple simplifying assumptions about the complex multi-step carcinogenesis process. We suggest that for the purposes of generating a preliminary integrated model, the consequent reduction in biological realism is compensated by the reduction in model complexity and number of parameters. More complicated examples of unifying long- and short-term models are certainly possible, e.g., models analyzing multiple pre-malignant cell stages, analyzing genomic instability, and/or analyzing stochastic effects at both time scales, rather than just at short-time scales as we shall do.

Short-term model

The short-term part of our model (Fig. 1a) is the more mathematically intensive part of the integrated formalism. It is based on existing initiation, inactivation and repopulation (iir) models originally designed for analyzing second cancers induced by radiotherapy (Lindsay et al. 2001; Sachs and Brenner 2005; Shuryak et al. 2006; Little 2007; Sachs et al. 2007). This short-term part analyzes normal and pre-malignant stem cells during complex radiation exposure regimens and during the following weeks of tissue recovery. It tracks individual pre-malignant cell clones rather than just the total number of pre-malignant cells. The probability that a pre-malignant clone becomes extinct during radiotherapy can be substantial, so a fully stochastic formalism is used for the population dynamics of pre-malignant cells.

Long-term model

The long-term part of our model (Fig. 1b) approximates carcinogenesis by a process where normal target (e.g., stem) cells can be initiated by spontaneous or radiation-induced mutations to become pre-malignant cells, which can clonally expand and perhaps produce by mutation a fully malignant cell, which then gives rise to cancer after some lag period. The basic assumptions are somewhat similar to those of several long-term models cited above, including the TSCE model. Specifically, our long-term formalism uses the following main assumptions:

  • Normal and pre-malignant stem cells are located in specific tissue niches (Slack 2000; Borthwick et al. 2001; Potten and Booth 2002; Bennett et al. 2003; Fuchs et al. 2004; Ghazizadeh and Taichman 2005; Li and Xie 2005). We here have in mind a general concept of a niche, including, e.g., a clone whose further growth is prevented or slowed by any kind of micro-environmental constraints. The number of niches per organ, and the number of normal or pre-malignant stem cells per niche, is homeostatically regulated (Fuchs et al. 2004; Li and Xie 2005). We allow for the possibility that pre-malignant stem cells are regulated somewhat less stringently than their normal counterparts, so that a niche filled with pre-malignant cells can contain more such cells on average than a niche filled with normal cells.

  • An initiated stem cell either dies out, or grows into a pre-malignant clone, which quickly (e.g., within a year, Campbell et al. 1996) fills the entire stem cell niche in which it originated, i.e., “initiates” the whole niche. To expand beyond the first niche, the pre-malignant clone needs to invade an adjacent niche; alternatively, the niche containing the clone can divide into two daughter niches, e.g., colon crypt fission (Greaves et al. 2006; Johnston et al. 2007; Edwards and Chapman 2007). These processes require years or decades, and may involve acquisition of new mutations (Spencer et al. 2006).

  • Pre-malignant cells in all niches are assumed to gradually lose their carcinogenic potential with age, so that at old age they have a progressively smaller probability of being transformed to malignant cells. This assumption is suggested by a decrease in incidence of most cancers at very old age, e.g., >80 years in humans [Surveillance Epidemiology and End Results (SEER) database, http://seer.cancer.gov] and >800 days in mice (Pompei et al. 2001). The likely mechanism is senescence of stem cells and/or deterioration of stem cell function, or niche function, with age (Brunet and Rando 2007; Carlson and Conboy 2007; Sharpless and DePinho 2007).

  • To reduce the number of model parameters, clonal expansion of pre-malignant stem cell niches is treated deterministically, using a net proliferation rate. In other words, it is estimated by a deterministic “exponential model”, instead of the stochastic TSCE model. Given realistic parameter values, the deterministic and stochastic versions of the two-stage model are numerically very similar from birth until old age (Heidenreich and Hoogenveen 2001). In our analyses, by the time irradiation is over and cell repopulation in the exposed organ has occurred, surviving pre-malignant clones have grown to substantial sizes, making extinction very unlikely. This reduces the need for stochastic treatment of clone dynamics during the decades following radiation exposure. At old age, the deterministic approximation predicts an increasing hazard, whereas the stochastic TSCE model predicts a plateau. Considering the evidence cited above about a decrease in cancer incidence at old age, neither behavior is fully realistic. We here instead model the downturn in incidence at old age explicitly, according to the assumption above.

Unification of short- and long-term models

The integration of long- and short-term models will have two features typical of multi-timescale modeling (Engquist and Runborg 2005): (1) information is passed in both directions between the two components; and (2) a formally infinite time interval in the short-term model represents a short time interval, here typically several months, in the long-term model (Fig. 1).

The deterministic long-term equations provide the mean number of niches filled with pre-malignant stem cells and the mean number of pre-malignant cells per niche just before irradiation. The stochastic short-term equations then provide the number of these niches that are eradicated by radiation as well as the number of pre-malignant clones that are induced by, and survive, the radiation exposure. Each of these clones is assumed to be independent (i.e., located in a different part of the organ), and to be capable of quickly filling a niche with pre-malignant cells. Therefore, the total number of pre-malignant niches soon after irradiation can be calculated by the stochastic short-term model. The mean of this number is the initial condition for the deterministic long-term equations, which are applied from this point onwards until old age (Fig. 1). Cancer incidence is assumed to be proportional to the total number of pre-malignant cells in all niches, shifted by a lag time. A mathematical description of the model is provided below.

Materials and methods

The mathematical techniques for implementing the assumptions for our specific example of integrating long- and short-term modeling are discussed next. A schematic representation of the concepts is provided in Fig. 2. The notation and interpretations for model parameters are listed in Table 1. Data on both background and radiation-induced cancers are used to estimate the parameters.

Fig. 2
figure 2

Schematic representation of proposed pre-malignant cell kinetics. Each square represents a tissue niche. Cyan color normal cells; other colors pre-malignant cells—each color indicates a different clone. Black arrows invasion of an adjacent niche. The order of panels denotes a hypothetical time sequence: a only niches filled with normal cells are present at a young age. b Some niches filled with pre-malignant cells are spontaneously initiated (red and blue). ce The pre-malignant clones expand over time by taking over adjacent niches. f Radiation exposure (lightning symbols) perturbs the system: some normal and pre-malignant cells are inactivated (niches become smaller), and some niches are completely inactivated (shown by gaps). Some new pre-malignant clones are initiated (green and brown). g After exposure, all niches filled with pre-malignant cells (red, blue, green, brown) expand in cell number (undergo promotion, shown by an increase in the size of the squares). Possible promotion of normal niches is neglected. Surviving normal and pre-malignant niches replace the niches inactivated by radiation (shown by disappearance of the gaps). h Over many years/decades after irradiation, roughly normal tissue architecture is restored. The promoting effect may disappear, as all niches return to the default size. Meanwhile, replication of existing pre-malignant niches continues (shown by arrows), and some new clones appear by spontaneous initiation (shown by the black square). i Subsequent mutations in one of the pre-malignant clones produce the first fully malignant cell, which grows into a malignant tumor (gray mass)

Table 1 Summary of model parameters

Long-term model for background cancers

The long-term model was intended mainly to place the results of the stochastic short-term calculations in an appropriate context, enabling estimation of the effects of age at exposure and time since exposure on predicted cancer risk. Simplicity and parsimony were emphasized. The long-term model in the absence of radiation consists just of a few rather simple deterministic equations, as follows.

Each pre-malignant stem cell in any niche has a certain small probability per unit time (q units = cells−1 × time−1) of transforming into a fully malignant cell, and, after a fixed lag time (L), into clinical cancer. These are common assumptions made in many long-term carcinogenesis models, e.g., in the TSCE model. The average number of new fully malignant cells per unit time (A units = time−1) is the product of q and three other variables: the number of stem cell niches filled with pre-malignant cells (M units = niches), the average number of pre-malignant stem cells per niche (ρ units = cells × niche−1), and the probability that the pre-malignant cells remain non-senescent and capable of producing cancer (P). The number of pre-malignant stem cells per niche (ρ) is assumed to be homeostatically regulated (in a qualitatively similar manner to the regulation of normal stem cells), so that it always tends towards a constant number—a carrying capacity Z (units = cells × niche−1). For convenience, the number of pre-malignant stem cells per niche can be redefined as a dimensionless normalized fraction C = ρ/Z. The proportionality constant qZ can be removed from the expression for A by defining N = qZM (units = time−1), where N is proportional to the number of pre-malignant niches. Consequently, A can be written in the following simplified form:

$$ A = NCP $$
(1)

Equation (1) is an approximation for the incidence hazard function: it gives the mean expected number of new malignant cells (i.e., eventual cancers) per individual per year, whereas the hazard (H) given in (2) is the yearly probability that a malignant cell occurs in a previously healthy individual at age t:

$$ H = A/\left( {1 - \int_{0}^{t} {A{\text{d}}u} } \right) $$
(2)

The functions A and H are numerically very similar for realistic parameter values; they diverge substantially only if cancer incidence is high. We use the exact hazard H for data fitting in the companion paper, but use the simpler expression for A in the equations below, keeping in mind its interpretation and limitations.

The function N is described by the following differential equation (3), where t is patient age, the constant a (units = time−2) is proportional to the spontaneous stem cell initiation rate, and the constant b (units = time−1) is the pre-malignant niche replication rate:

$$ {\text{d}}N/{\text{d}}t = a + bN $$
(3)

Equation (3) can be adjusted to accommodate an assumption that more than one mutation is needed to initiate a stem cell by replacing the constant a with the composite term at κ−1, where κ is the necessary number of mutations. However, doing so does not substantially improve the fit of the model to the data sets analyzed, and so was not used as the default because it introduces an extra adjustable parameter κ. The results for fitting multi-stage extensions of our model (i.e., where κ > 1) to SEER data for female breast and male lung cancers are shown in Fig. 3.

Fig. 3
figure 3

Fits of the long-term part of our model to spontaneous incidence data for female breast and male lung cancers from SEER (the data for young ages were excluded because cancer incidence in these age groups is likely to be dominated by genetic predisposition). Error bars represent 95% confidence intervals. The different curves represent our default two-stage formalism and its multi-stage extensions described in the main text. These extended versions of the formalism do not alter the fit substantially

The function C is regulated by the following logistic differential equation, where δ (units = time−1) is a rate constant representing the strength of homeostatic control of the number of pre-malignant cells per niche:

$$ {\text{d}}C/{\text{d}}t = \delta \, C\left( {1 - C} \right) $$
(4)

For background carcinogenesis the solution C = 1 is used.

The probability that the pre-malignant cells remain non-senescent and capable of producing cancer (P) is described by a Gaussian function with an adjustable parameter c (units = time−2), which cannot become negative even as age approaches infinity, contrary to the expression used by Pompei et al. (2001) and Pompei and Wilson (2002):

$$ P = \exp \left[ { - ct^{2} } \right] $$
(5)

For convenience later, when radiation-induced excess relative risk (ERR) will be calculated, it is useful to express patient age (t) as the sum of age at exposure (T x ) and time since exposure (T y ). The background cancer risk function without radiation (A bac(T x , T y )) follows from (1) and (35), assuming that no pre-malignant cells are present at birth (i.e., that N(0) = 0):

$$ A_{\text{bac}} \left( {T_{x} , \, T_{y} } \right) = \left( {a/b} \right)\left( {\exp \left[ {b\left( {T_{x} + T_{y} } \right)} \right] - \, 1} \right)\exp \left[ { - c\left( {T_{x} + T_{y} } \right)^{2} } \right] $$
(6)

Here, the duration of radiation exposure is ignored, because it is considered to be short compared with the multi-year time scale of the long-term carcinogenesis processes.

Effects of radiation

Short-term processes

The exposure scenario analyzed here is radiotherapy for an existing malignancy, where there are K daily dose-fractions, all of equal size d in some nearby organ, with treatment gaps during the weekends. Straightforward generalizations to variable doses per fraction, to other temporal patterns, and/or to cases where one must consider dose-volume histograms are omitted for brevity. A single acute dose exposure is a simple special case, where the number of fractions is K = 1.

The short-term part of the model considers initiation, inactivation, and repopulation (iir) of normal and pre-malignant stem cells during the radiation regimen and a recovery period of about a month (Fig. 1). Niche takeover by pre-malignant cells that survive the irradiation and modulation of niche size by radiogenic promotion (discussed below) are here also considered short-term processes. They are assumed to be essentially completed by a few months following exposure, before the long-term processes of spontaneous initiation and niche replication, which operate on the time scale of multiple years, start to have an appreciable effect (Fig. 1).

Normal stem cells are treated deterministically because their number per organ is assumed to be large. In contrast, the number of pre-malignant stem cells per clone is much smaller, and extinction of some pre-malignant clones is a real possibility. Consequently, a stochastic approach described below is used to estimate the average number of live pre-malignant niches when post-irradiation long-term processes begin.

For simplicity, both normal and pre-malignant stem cells are assumed to be equally radiosensitive. The cell surviving fraction (S) is described by the standard linear-quadratic (LQ) function with parameters α and β for dose-fractions of size d:

$$ S = \exp \left[ { - \alpha d - \beta d^{2} } \right] $$
(7)

The normal stem cell number (n) in the entire organ just before the kth dose-fraction is denoted by n (k), and the number just after the kth dose-fraction by n +(k). The reduction in n due to initiation of a few normal stem cells is neglected. Inactivation is calculated as follows:

$$ n^{ + } \left( k \right) = Sn^{ - } \left( k \right) $$
(8)

Repopulation of normal stem cells during radiotherapy is assumed to be regulated by homeostatic mechanisms—if some stem cells are killed by radiation, surviving stem cells are induced to proliferate (Dörr and Kummermehr 1990; Pabst et al. 2004) with the goal of restoring the total number of normal stem cells in the organ (n) to the number present before irradiation (ν). It should be noted that ν is conceptually and numerically distinct from the carrying capacity Z introduced above for pre-malignant stem cells per pre-malignant niche. The homeostatic regulation of cell proliferation during the time gaps between dose-fractions and after the last fraction is described by the Logistic differential equation, where λ (with units = time−1) is the maximum net proliferation rate:

$$ {\text{d}}n/{\text{d}}t = \lambda n\left[ {1 - n/\nu } \right] $$
(9)

Equations (79) allow the calculation of the normal stem cell number n(t) in the entire organ at all relevant times t throughout the radiation exposure and recovery periods. Explicit results for n(t) are provided in the Appendix.

We assume that the long-term growth advantage of pre-malignant cells manifests itself only on the scale of years and decades, and is negligible on the much shorter time scale of a few weeks of radiotherapy. Consequently, the maximum net proliferation rates for normal and pre-malignant stem cells are assumed to be equal, described by the parameter λ.

As discussed above, radiation initiates some normal stem cells, making them pre-malignant. The number of cells initiated in the kth dose-fraction is assumed to be Poisson distributed, with average aXI(k), where X is an adjustable parameter (units = time × dose−1) and the function I(k) (units = dose) is given by:

$$ I\left( k \right) = dSn^{ - } \left( k \right)/\nu = dn^{ + } \left( k \right)/\nu $$
(10)

Our assumptions allow the birth, death, and initiation rates for normal stem cells calculated deterministically by (710) to be used as parameters for a stochastic formalism for pre-malignant cell clones. Each such clone is initiated as a single cell in some random niche by the kth dose-fraction, and can fluctuate in cell number during subsequent radiotherapy due to the opposing effects of inactivation and repopulation. We count the number of clones that contain at least one viable cell when radiotherapy ends, because only these surviving clones are capable of eventually taking over their niches. Their number is determined by the probability F(k) that a live stem cell initiated by the kth dose-fraction produces a clonal lineage, which survives all subsequent dose-fractions. Using analytic results on stochastic birth–death processes with variable rates (Tan 2002; Hanin 2004) the Appendix derives an equation, (22), for F(k), which is repeated here:

$$ F(k) = {\frac{1}{\Updelta (k)}} ,\quad {\text{where }}\Updelta (k) = n^{ + } (k)\left[ {1 + \sum\limits_{j = k}^{K} {{\frac{{n^{ - } (j + 1) - n^{ + } (j)}}{{n^{ - } (j + 1)n^{ + } (j)}}}} } \right]. $$

Here, F(0) is the probability that a pre-malignant cell that was present before irradiation began produces a lineage which survives all dose-fractions.

Using the facts that a random thinning of a Poisson distribution is Poisson and the sum of independent Poisson distributions is Poisson, one can show that the total number of surviving pre-malignant clones at the end of the short-term period is Poisson distributed. The above arguments give the mean value as N init = aXISf(D), where D is the total radiation dose (i.e., the sum of all doses per fraction, dK) and ISf(D) (units = dose) represents a net outcome of initiation, inactivation, and cell repopulation during exposure, \( ISf\left( D \right) = \sum\nolimits_{k = 1}^{K} {I(k)} F(k). \) The probability that a pre-malignant niche that was present before exposure is not inactivated by radiation (i.e., that at least one pre-malignant stem cell in the niche survives) is given by:

$$ Sf\left( {Z,D} \right) = 1 - \left[ {1 - F\left( 0 \right)} \right]^{Z} $$
(11)

Unifying short- and long-term processes

The above short-term processes are regarded as effectively instantaneous relative to the long-term ones. The stochastic results for the short-term exposure period, in the form of the functions ISf(D) and Sf(Z, D), are inserted into the deterministic equations for long-term carcinogenesis processes. On the long-term time scale, niches already pre-malignant at the end of irradiation and recovery then increase in number by replication; we will refer to all the resulting niches as “old” niches. Other, “new” niches are formed by spontaneous initiation after the end of the exposure period. The contributions of niches in these two categories to the cancer risk are called N rad E(T x , T y ) and N rad N(T x , T y ), respectively. At a given time (T y ) after irradiation they are given by the following solutions for (3):

$$ \begin{aligned} N_{\text{rad}} E\left( {T_{x} ,T_{y} } \right) = & \left[ {\left( {a/b} \right)\left( {\exp \left[ {bT_{x} } \right] - 1} \right)Sf\left( {Z,D} \right) + aXISf\left( D \right)} \right]\exp \left[ {bT_{y} } \right] \\ N_{\text{rad}} N\left( {T_{x} ,T_{y} } \right) = & \left( {a/b} \right)\left( {\exp \left[ {bT_{y} } \right] - 1} \right) \\ \end{aligned} $$
(12)

Radiation is commonly assumed to promote hyper-proliferation of pre-malignant cells (Heidenreich et al. 2007). Here, we interpret this effect as an increase in the number of pre-malignant stem cells per surviving niche above the niche carrying capacity Z, i.e., C becomes >1. For simplicity, the initial excess of C is assumed to be linearly dependent on radiation dose with the coefficient Y (units = dose−1). Promotion may be eventually reversed, because the number of pre-malignant cells per pre-malignant niche may gradually return to pre-irradiation carrying capacity Z. This process may occur concurrently with extinction of some radiation-induced niches (Zhang et al. 2001; Ullrich 1986). Since only the product of the number of niches and the number of cells per niche is relevant for cancer risk (1), extinction of some niches and/or shrinkage of niche size do not need to be modeled separately. The net effect—i.e., a gradual reversal of promotion—can be modeled using only one adjustable parameter (δ), as done here. According to these assumptions, at any given time (T y ) after exposure to radiation the average normalized number of pre-malignant cells per surviving niche (C rad(T x , T y )) can be calculated by solving (4):

$$ C_{\text{rad}} \left( {T_{x} ,T_{y} } \right) = \left( {1 + YD} \right)/\left[ {1 + YD \, \left( {1 - \exp \left[ { - \delta T_{y} } \right]} \right)} \right] $$
(13)

Calculation of absolute and relative cancer risks

The approximation for the absolute cancer risk after radiation (A rad(T x , T y )) can now be calculated at any age at exposure and time since exposure by using the equations above:

$$ \begin{gathered} \, A_{\text{rad}} (T_{x} , \, T_{y} ) = \left[ {N_{\text{rad}} E(T_{x} ,T_{y} )C_{\text{rad}} (T_{x} ,T_{y} ) + N_{\text{rad}} N(T_{x} ,T_{y} )} \right]P(T_{x} ,T_{y} ) = \\ {\frac{a}{b}}\left[ {{\frac{{(\exp [b \, T_{x} ] - 1)Sf(Z,D) + bX \, ISf(D)\exp [bT_{y} ](1 + YD)}}{{1 + YD(1 - \exp [ - \delta T_{y} ]) \, }}} + \exp [bT_{y} ] - 1} \right]\exp \left[ { - c(T_{x} + T_{y} )^{2} } \right] \\ \end{gathered} $$
(14)

The excess relative cancer risk (ERR) follows from (6) and (14): ERR = [A rad(T x , T y )/A bac(T x , T y )] − 1. By substitution and simplification, a more explicit expression for the ERR can be obtained:

$$ \begin{aligned} {\text{ERR}} = & \left[ {\left( {{\text{Q}}_{1} {\text{Q}}_{2} + {\text{Q}}_{3} } \right)/{\text{Q}}_{4} } \right] - 1,\quad {\text{where}} \\ {\text{Q}}_{1} = & \left( {1 + YD} \right)/\left[ {1 + YD \, \left( {1 - \exp \left[ { - \delta T_{y} } \right]} \right)} \right]; \\ {\text{Q}}_{2} = & \left[ {\left( {\exp \left[ {bT_{x} } \right] - 1} \right)Sf\left( {Z,D} \right) + bXISf\left( D \right)} \right]\exp \left[ {bT_{y} } \right]; \\ Q_{3} = & \exp \left[ {bT_{y} } \right] - 1; \\ Q_{4} = & \exp \left[ {b \, \left( {T_{x} + T_{y} } \right)} \right] - \, 1 \\ \end{aligned} $$
(15)

Note that the senescence parameter (c) cancels out of this ERR expression; due to the way we have defined our radiation initiation parameter (X), the spontaneous initiation parameter (a) also cancels out. The term Q 1 can be interpreted as the normalized size of old pre-malignant stem cell niches. Q 2 is proportional to the number of such niches. Q 3 is proportional to the number of new pre-malignant niches, and Q 4 is proportional to the total number of pre-malignant niches under background conditions.

Results

To display the properties of the model, we produced Figs. 47 using generic parameter values guided by the best-fit values for specific cancer sites found in the accompanying paper (Shuryak et al. 2009).

Fig. 4
figure 4

The typical shape for age dependence of background incidence for a specific type of adult-onset solid cancer is reasonably reproduced by our model. The curve was generated using the following parameter values: a = 1.0 × 10−8 y−2, b = 0.25 y−1, c = 1.75 × 10−3 y−2, L = 10 y

Figure 4 shows an example of a model-generated age-dependent background incidence curve for an adult-onset solid cancer, using the four relevant parameters a, b, c, and L. The shape of this curve agrees well with the data for many cancers (e.g., SEER database, http://seer.cancer.gov), including the downturn in incidence at old age, which is not described as well by standard models.

Figure 5 shows a typical radiation dose–response, which is determined by the balance of cell initiation, inactivation (killing), and repopulation. Fractionation of the dose increases cancer ERR because repopulation of both normal and pre-malignant cells during the inter-fraction intervals compensates for much of the cell killing.

Fig. 5
figure 5

The effect of dose-fractionation on predicted excess relative cancer risk (ERR): as the same total radiation dose is split into fractions (one fraction per day, with gaps on weekends), thereby protracting it over a longer time, cancer risk predicted by the model increases. This occurs because cell repopulation during prolonged exposure partially compensates for cell killing by radiation. The doses refer to a given organ, such as the lungs or female breast, and not to whole body exposures. The following parameter values were used: b = 0.25 y−1, X = 10.0 y × Gy−1, Y = 0.5 Gy−1, δ = 0.01 y−1, Z = 10.0 cells/niche, α = 0.3 Gy−1, β = 0.0 Gy−2, λ = 0.35 day−1, T x  = 30 y, T y  = 0.0 y, and L = 10 y

Figures 6 and 7 show the effects of age at exposure and time since exposure. The components of cancer ERR produced by radiation-induced initiation and radiation-induced promotion exhibit very different dependences on age at exposure. Radiation is assumed to initiate the same number of cells per unit dose independent of age, whereas the background number of pre-malignant cells, which is essentially the denominator of ERR, grows with age. Consequently, the initiation-driven component of ERR decreases with age at exposure. This process dominates the ERR for ages <20 years in Fig. 6. In contrast, promotion is assumed to be a multiplicative amplification of the background number of pre-malignant cells per niche. Consequently, the promotion-driven component of ERR is approximately constant over most ages at exposure. This process dominates the ERR for older ages in Fig. 6.

Fig. 6
figure 6

Combined effects of age at exposure and time after exposure on radiation-induced excess relative cancer risk (ERR). The ERR is calculated at age 70, which is intended to approximate lifetime risk. As discussed in the text, only age at exposure matters if parameter δ = 0, but if δ > 0, time after exposure affects the ERR as well. The parameter values were the same as in Fig. 5, except for T x and T y , and a reference radiation dose of 1 Gy

Fig. 7
figure 7

Modulation of promotion-driven excess relative cancer risk (ERR) over time after exposure: a low doses (either acute or fractionated produce similar results, since cell inactivation is negligible), b high fractionated doses (fractions are given once daily, with gaps on weekends). It is seen that, contrary to a proportional hazards assumption, not only the amplitude but also the shape of the dose–response curve is altered by long-term processes. The parameters were the same as in previous figures

Promotion-driven ERR can also be modulated by time after exposure. This occurs due to the assumption that the number of pre-malignant cells per niche is homeostatically regulated (by parameter δ), so that radiation-induced hyper-proliferation of cells within their niches can be reversed, as pre-irradiation cell birth/death rates, and hence niche sizes, are restored. If δ > 0, ERR due to promotion will decrease over time following exposure. This effect is seen in Figs. 6, 7: If the risk is measured at some constant age (e.g., 70 years), which is a sum of age at exposure and time since exposure, a decrease in ERR with time since exposure due to a δ > 0 will appear as an increase in the ERR with age at exposure (Fig. 6). A decrease in promotion-driven ERR with time since exposure also can have a conceptually important effect on the dose–response—at longer times after exposure, not only the magnitude, but also the shape, of the dose–response can change (Fig. 7).

Discussion

The formalism presented here is the first comprehensive attempt to unify short- and long-term modeling approaches. The short-term part of the model belongs to the previously discussed iir category. It tracks the numbers of pre-malignant cells throughout irradiation stochastically. The long-term part of the model builds on the concepts developed in previous two-stage formalisms by adding an analysis of some aspects of tissue architecture (i.e., stem cell niches/compartments) and aging of pre-malignant stem cells. The particular short and long-term models that we have chosen to use are not crucial—the real issue is the integration. Certainly other long-term models could perfectly well be integrated into this short-/long-term framework—and we hope they will be.

The unified formalism has a number of advantages. The short-term part can generate reasonable predictions even at high doses, such as those in cancer radiotherapy. The long-term part analyzes the entire lifetime of the individual, putting the short-term predictions in an appropriate context by estimating the effects of age at exposure and time since exposure. The combined approach therefore allows the dose–response for the number of pre-malignant cells to be examined at any time point, from the start of irradiation until development of cancer years to decades later, which is not possible using either short- or long-term models alone.

Our model can be used for estimating risks of second malignancies induced by radiotherapy. This issue is growing in importance (Brenner et al. 2000; Ron 2006) as patients are treated earlier in life and the number of cancer survivors increases (Editorial 2004); the lifetime risk of radiation-induced second cancers is not negligible (Brenner et al. 2007). Direct measurement of second cancer incidence requires decades of follow-up because the latency period for radiogenic solid tumors is long (Tokunaga et al. 1979; Brenner et al. 2000; Ivanov et al. 2004, 2009). Meanwhile, radiotherapy protocols are rapidly changing. Our model, calibrated using data from older protocols, can produce risk predictions for any modern or prospective radiotherapeutic protocol. This application of the model is discussed in the accompanying paper (Shuryak et al. 2009).