Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA. E-mail: moolgavkar{at}earthlink.net
Fifty years after its publication, the multistage model proposed by Armitage and Doll1 continues to influence biological and epidemiological thinking on the processes underlying carcinogenesis. I am pleased that the International Journal of Epidemiology has chosen to republish Armitage and Doll's landmark paper, and I am honoured to write a commentary on it.
With the rapid advances in molecular biology over the last 50 years, we know a lot more about the process of carcinogenesis than we did when Armitage and Doll wrote their paper. Yet, the fundamental predicate on which their model is based, that malignant transformation of a normal cell results from the accumulation of a few critical mutations, remains unchallenged. To be sure, the Armitage-Doll model needs to be embellished in various ways to accommodate our current thinking on carcinogenesis. Perhaps the most important of these embellishments is the allowance for clonal expansion of partially altered cells on the pathway to cancer. Even this development was anticipated by Armitage and Doll in their follow-up paper in 1957.2
Building on ideas of Fisher and Holloman3 and Nordling,4 Armitage and Doll proposed that the remarkable regularity observed in the age-specific mortality rates for many adult carcinomas could be explained by a multistage model for carcinogenesis. But they went a lot further. They recognized that some carcinomas, particularly the hormone-dependent carcinomas, showed departures from this behaviour and considered reasons for such departures. In particular, they investigated the implications of carcinogens acting on various stages of the process on the temporal evolution of cancer risk. This work foreshadowed and set the stage for later work by Whittemore,5 Day and Brown,6 and Brown and Chu,7 among others.
Following Fisher and Holloman and Nordling, Armitage and Doll restricted their attention to cancer mortality between the ages of 25 and 74 because the data were considered to be unreliable for older age groups. In this range of ages many of the carcinomas they examined exhibited log-log behaviour, i.e. the logarithm of age-specific mortality increased linearly with the logarithm of age with a slope of about 6, which accorded well with a model requiring seven stages for malignant conversion of a normal cell. In later analyses, however, it was noted8 that at higher ages the age-specific rates showed substantial departures below the rates predicted by log-log behaviour. Doll and Peto9 reported that lung cancer risk among men aged 8084 appeared to be half that among men aged 7579. They offered four possible explanations for this finding. These were under-diagnosis, selective survival, unreported cohort differences of smoking patterns in early life, and the possibility that the biology of extreme old age reduces the risk of carcinoma. A more general explanation is that human populations are heterogeneous with respect to their susceptibility to cancer, either for genetic or environmental reasons. As populations age, the more susceptible individuals get cancer so that at extreme old ages one observes cancer in only the least susceptible individuals. In a later paper, Stevens and Moolgavkar10 reported that after adjustment for cigarette smoking and other temporal trends, lung cancer mortality rates among British males 8084 were not lower than the rates among those 7579. Nevertheless age-specific rates showed significant departures below the predictions of log-log behaviour. In fact, any general multistage model makes precisely this prediction as I discuss below.
In their original paper describing the multistage model, Armitage and Doll used a simple approximation to the exact hazard (incidence) function. While this approximation was perfectly adequate for their purposes, the exact hazard functions of fairly general multistage models have important properties that the approximations do not possess. First, the (exact) incidence functions predicted by the models approach a finite asymptote, i.e. they cannot increase without bound.11,12 For the Armitage-Doll model it is easy to see that this asymptote is Nmin, where N is the total number of susceptible cells and
min is the minimum of the transition (mutation) rates between stages.12 For the log-log cancers this fact implies that with increasing age, the age-specific incidence will increasingly deviate below the line defined by the log-log relationship.
The second consequence following from consideration of the exact hazard functions of multistage models has to do with the behaviour of the age-specific incidence curve after exposure to an environmental carcinogen ceases. Suppose that exposure to an environmental agent increases the rate of one or more of the steps in the carcinogenic pathway, e.g. by increasing the rates of specific mutations or of clonal expansion of initiated cells. Suppose also that these rates revert to background rates after exposure stops. Then the exact (but not the approximate) hazard functions asymptotically approach the hazard function in those individuals who were never exposed.11 This phenomenon has been observed among ex-smokers, for example. In many epidemiological studies incidence rates of lung cancer in ex- smokers appear to approach the incidence rates in non-smokers about 15 years after quitting smoking. One possible explanation for this phenomenon is that repair of damaged tissue has occurred after exposure stopped. This could well be true, but this explanation need not be invoked. The reversion of the incidence function to background rates is a mathematical consequence of multistage carcinogenesis.
The third consequence has to do with the incidence of second malignant tumours among individuals in whom one has already occurred. A computation using the exact solution shows that the age-specific incidence of second malignant tumours is higher than the age-specific incidence of the first malignant tumour (at the same age). While the incidence of second malignant tumours is difficult to study in human populations because of the treatment intervention after the occurrence of the first tumour, animal experiments appear to show that this is indeed true.13 The explanation has been advanced that physiological and immunological changes after the first malignancy renders the animal susceptible to a second tumour. This may be true but the higher incidence of second tumours is a logical consequence of multistage carcinogenesis.14
It is a rare paper that continues to be widely cited 50 years after publication. The paper describing the Armitage-Doll modeloften referred to as THE multistage modelis one such rare paper. I congratulate Professors Armitage and Doll on the occasion of the republication of this important and influential paper.
![]() |
Acknowledgments |
---|
![]() |
References |
---|
![]() ![]() |
---|
2 Armitage P, Doll R. A two-stage theory of carcinogenesis in relation to the age distribution of human cancer. Br J Cancer 1957;11:16169.[ISI][Medline]
3 Fisher JC, Holloman JH. A hypothesis for the origin of cancer foci. Cancer 1953;4:91618.[ISI]
4 Nordling CO. A new theory of the cancer inducing mechanism. Br J Cancer 1953;7:6872.[ISI][Medline]
5 Whittemore AS. The age distribution of human cancer for carcinogenic exposure of varying intensity. Am J Epidemiol 1977;106:41832.[ISI][Medline]
6 Day NE, Brown CC. Multistage models and the primary prevention of cancer. J Natl Cancer Inst 1980;64:97789.[ISI][Medline]
7 Brown CC, Chu KC. Implications of the multistage theory of carcinogenesis applied to occupational arsenic exposure. J Natl Cancer Inst 1983;70:45563.[ISI][Medline]
8 Cook PJ, Doll R, Fellingham SA. A mathematical model for the age distribution of cancer in man. Int J Cancer 1969;4:93112.[ISI][Medline]
9 Doll R, Peto R. Cigarette smoking and bronchial carcinoma: dose and time relationships among regular smokers and life-long non-smokers. J Epidemiol Community Health 1978;32:30313.[Abstract]
10 Stevens RG, Moolgavkar SH. A cohort analysis of lung cancer and smoking in British males. Am J Epidemiol 1984;119:62441.[Abstract]
11 Heidenreich W, Luebeck EG, Moolgavkar SH. Some properties of the hazard function of the two-mutation clonal expansion model. Risk Anal 1997;17:39199.[ISI][Medline]
12 Moolgavkar SH, Krewski D, Schwarz M. Mechanisms of carcinogenesis and biologically-based models for quantitative estimation and prediction of cancer risk. In: Moolgavkar SH, Krewski D, Zeise L, Cardis E, Moller H (eds). Quantitative Estimation and Prediction of Cancer Risk. Lyon: IARC Scientific Publications, 1999, pp. 179238.
13 Prehn RT. Nongenetic variability in susceptibility to oncogenesis. Science 1975;90:109596.
14 Dewanji A, Moolgavkar SH, Luebeck EG. Two-mutation model for carcinogenesis: joint analysis of premalignant and malignant lesions. Math Biosci 1991;104:97109.[CrossRef][ISI][Medline]