Objective. To develop a conceptual model for lung cancer rates to describe and quantify observed differences between Sweden and USA contra Japan. Method. A two-parameter lognormal distribution was used to describe the lung cancer rates over time after a 1-year period of smoking. Based on that risk function in combination with smoking prevalence, the calculated age-standardized rates were adjusted to fit reported data from Japan, Sweden, and the USA by parameter variation. Results. The risk of lung cancer is less in Japan than in Sweden and in the USA at the same smoking prevalence and intensity. Calculated age-specific rates did also fit well to reported rates without further parameter adjustments. Conclusions. This new type of cancer model appears to have high degree of predictive value. It is recommended that data from more countries are studied to identify important life-style factors related to lung cancer.
1. Introduction
Prior to 1955, the relationship between smoking tobacco and lung cancer was not recognised as it is today. A 1955 Swedish encyclopaedia describes nervous problems, stomach problems, and throat problems, but not lung cancer, as possible results of heavy smoking [1]. However, after 1955, lung cancer mortality began to suddenly increase in Sweden, and older people in particular were subjected to a fast increasing risk of lung cancer.
The aim of this study is to develop a conceptual model for predicting lung cancer rates. Such means may be useful in promoting successful preventive work against lung cancer. In cancer epidemiology, statistical models use linear and logistic regression to evaluate relationships between risk factors and cancer incidence. Less frequently used are biomathematical models translating a hypothesis about the biological process into mathematical terms [2].
There are differences between countries worldwide regarding lung cancer even if they have similar cigarette consumption per inhabitant and year [3]. This indicates that there are other causative factors, as well as different genetics, still not considered or yet identified in epidemiological studies of lung cancer.
The model developed for this study was based on a similar model, successfully used for estimating melanoma rates over time, taking reduced repair efficiency of sun-induced skin damages into consideration [4]. In the current study, the model developed for lung cancer had to account for varying smoking prevalence over time.
The main finding from this study is that the model with the parameters used was capable of predicting age-specific incidence for both men and women in Sweden and for men in the USA and Japan. The model study indicates that the lung cancer risk after one year of smoking is larger in Sweden and USA as compared to Japan. This conclusion was strengthened by the fact that western countries in general have twice the lung cancer rate as nonwestern countries at equal cigarette consumption per person and year. A study of associations between lung cancer rates and life-style factors other than smoking is thus recommended.
2. Materials and Methods
Data on smoking prevalence were collected for men and women from official sources [5]. Age-specific smoking prevalence data have been reported in Sweden since 1980 and the total smoking prevalence since 1916. Smoking prevalence data from the USA and Japan were retrieved from [6, 7]. The expected life duration for smokers and non-smokers in western countries were obtained from survival plots given by Doll et al. [8].
This study sought to define a one-year period of average smoking for men in Sweden as the basic damage-generating unit causing a lung cancer risk that would increase by time according to a statistical function. The smoking intensity per lung volume among Swedish male smokers was defined as the norm in this model. The smoking intensities for men in Japan and the USA and for women in Sweden were adjusted to best fit data and then compared with available reported data. The probability function was used to calculate cancer cases from each smoking period of one year over the life time of a birth cohort consisting of 100,000 men or women. The cancer cases from all these smoking periods were then summed over time (calendar years) in a life matrix (see Table 1) to give the primary calculated incidence. This primary incidence versus age was finally adjusted for Sweden and the USA to account for the reduced fraction of smokers alive at old age due to the longer life expectancy of nonsmokers in western countries [8].
A life matrix based on annual smoking damage (N) and cancer risk (F) to get the accumulated risk over life time (vertical sum).
Year
Year
Dmg
Y1
Y2
Y3
Y4
Y5
Y6
Y7
⋯
Y1
N1
N1F1
N1F2
N1F3
N1F4
N1F5
N1F6
N1F7
⋯
Y2
N2
N2F1
N2F2
N2F3
N2F4
N2F5
N2F6
⋯
Y3
N3
N3F1
N3F2
N3F3
N3F4
N3F5
⋯
Y4
N4
N4F1
N4F2
N4F3
N4F4
⋯
Y5
N5
N5F1
N5F2
N5F3
⋯
Y6
N6
N6F1
N6F2
⋯
Y7
N7
N7F1
⋯
⋮
⋮
∑Y1
∑Y2
∑Y3
∑Y4
∑Y5
∑Y6
∑Y7
⋯
This analysis was repeated for all birth cohorts to calculate the age-standardized incidence over time. The statistical functions used were characterized by two statistical parameters: the dispersion (in time, decades) and time to reach 0.1% (used to define the median value) of a lognormal distribution. The parameters were first varied to identify the best fit of calculated age-standardized rates to reported rates in Japan. To verify the usefulness of the model, calculated age-specific incidences were compared with reported data without further parameter variation. The dispersion giving best fit to reported rates in Japan was also used for Sweden and the USA, while the time to 0.1% was further fine-tuned to fit the Swedish and US data. We considered the dispersion mainly to be a characteristic of the disease (lung cancer) so that it should be used for all countries. The time to 0.1% risk of cancer is more likely to be a characteristic of the environment and should thus be used as an adjustable parameter to fit data.
3. Results
Figure 1 shows the basic lung cancer risks in Sweden and the USA and in Japan from one year of smoking. Figure 2 shows age-standardized rates for men in Japan after varying the two parameters (dispersion and time to 0.1%) to best fit calculated and reported data. Figure 3(a) shows that the calculated incidence for Swedish women would have been about twice as low as actually reported. In order to better fit reported data, the time to 0.1% had to be adjusted somewhat which is shown in Figure 3(b) and in Table 2. The same dispersion as for Japan was used also for men in Sweden and in the USA and for women in Sweden.
Log-normal parameters for lung cancer risk functions.
Population
Dispersion (decades)
Time to 0.1% (Years)
Japan male
0.30
45.13
Se male, Se female USA male
0.30; 0.30 0.30
41.22; 41.00 41.43
Accumulated lung cancer risk in Sweden and the USA and in Japan from one year of smoking.
Calculated and reported age-standardized incidence for men in Japan.
(a) Calculated and reported age-standardized incidence for women in Sweden using the same two parameters as were used for Japan in Figure 2. (b) Calculated and reported age-standardized incidence for women in Sweden after optimization of the parameter value for time to 0.1%.
We found the parameter values extracted from the statistics on women very similar to parameters to fit data on men in Sweden and in the USA (Figures 4 and 5). Table 2 gives the parameters for best fit to reported age-standardised rates, used to calculate all incidences in Figures 2–9. Figure 6 shows reported and calculated age-specific lung cancer incidence for men in Japan. Similar results were obtained for men in Sweden and for men in the USA (Figures 8 and 9). It is of interest to note that the rates for men and women older than 80 years since around 1995 have been increasing (Figures 7 and 8). The sharp increase after 1997 for corresponding older men in Sweden cannot be explained by the model, however.
Calculated and reported age-standardized incidence for men in Sweden after optimization of the parameter value for time to 0.1%.
Calculated and reported age-standardized incidence for men in the USA after optimization of the parameter value for time to 0.1%.
Calculated and reported age-specific lung cancer incidence for men over time in Japan. Data are given for age groups from 37 to 82 years, where R stands for reported data.
Calculated and reported age-specific lung cancer incidence for women over time in Sweden.
Calculated and reported age-specific lung cancer incidence for men over time in Sweden.
Calculated and reported age-specific lung cancer incidence for men over time in the USA.
A study of data from international cancer registries [9] in combination with reported cigarette consumption per capita [10] showed that western countries have twice the amount of lung cancer compared to nonwestern countries at similar smoking levels (See Figure 10).
Lung cancer incidence versus cigarette consumption in western and nonwestern countries.
4. Discussion
We have here been using a mechanistic model requiring just two parameters instead of the empirical methods normally used, for example, in age-period cohort (APC) models [11, 12]. The rate (R) at any age (a) and period (p) can be calculated asRap=exp(Aa+D*p+Pp+Cc),
Here Aa is a constant that varies by age, D is a constant that is multiplied by the period p, Pp is a specific constant added for each period p, and C is a specific constant for each cohort, for example, birth cohort.
In order to generate the age-specific incidence plots and age-standardised plots in Figures 2–9 using the APC model, quite a large number of parameters is necessary. Such a model of lung cancer mortality in relation to age, duration of smoking, and daily cigarette consumption was studied by Flanders et al. [13].
The connection between tobacco smoking and lung cancer has long been known (review in [14]). A comprehensive report on tobacco use and cancer was given in [6]. In order to take different smoking intensity patterns, defined as number of cigarettes per lung volume, between Sweden and USA into consideration, we allowed an intensity factor to be varied. Best fit was found when the smoking intensity was 3.3 times larger in the USA as compared to Sweden. The average annual weight of smoked tobacco per inhabitant during the period from 1922 to 2002 was 1.047 kg in Sweden and 3.318 kg in the USA or 3.2 times the Swedish weight. The computer also found a best fit for Swedish women using a smoking intensity of 1.19. Initially this was intriguing since Swedish smoking women tend to smoke 12 cigarettes per day compared with 15 for male smokers. However, taking their different lung volumes into consideration (4.2 vs. 6 litres), it turned out that 12/4.2*6/15=1.14, a value which was then used for the calculations (Table 3).
Smoking intensity relative to that of Swedish men.
Country, gender
Smoking intensity relative to Swedish men
Sweden, men
1
Sweden, women
1, 14
Japan, men
1, 24
USA, men
3, 30
We used information on life expectancy of smokers and nonsmokers in western countries to correct the model for decreasing ratios of smokers versus age among older people in Sweden and in the USA. The generally observed decline in cancer incidence for the oldest people was reviewed in [15]. Due to natural death, there will be an accumulation of survivors who have lived in good and healthy conditions and this will manifest itself in lower cancer rates. This also makes an argument to search for more specific information about what these good conditions really are.
The question why the overall cancer risk is higher in more developed countries was raised in [16]. One possible explanation was that the vulnerability has increased because the “immunity is affected by factors associated with economic progress and a western life-style.” Such factors may include diet, chemicals, or continuous exposure to immune-disturbing electromagnetic fields. Figure 10 shows lung cancer rates versus smoking habits in western and nonwestern countries.
Several studies report that there is a significant reduction in DNA repair capacity (DRC) found in lung cancer patients compared with healthy controls, [17, 18]. There are rapid methods for determination of DRC by blood tests to find candidates for intensive follow-up and possibly screening, [19]. By using such methods, it is also possible not only to compare lung cancer patients with healthy controls, but also to compare DRC among healthy smokers from different populations, for example, from Japan and the USA. One step in that direction was recently presented by Gandara et al. [20].
Our findings from this model study indicate a significant difference between the capability of coping with heavy smoking without getting lung cancer in Japan versus the USA, and Sweden. If DRC generally is better in Japan than in for example, the USA, then it is our firm belief that we immediately should study what environmental factors that might have the strongest influence on DRC and then use this information for better cancer prevention strategies.
Three major conclusions can be drawn from this study.
The aim of this study was to develop a conceptual model for predicting lung cancer rates. Such a model was indeed achieved and it will be an important task to follow it over time as well as applying it to statistical material from other countries.
A study of associations between lung cancer rates, DRC, and environmental and life-style factors other than smoking is strongly recommended.
Lung cancer rates in Sweden and the USA might become significantly reduced if we could identify and adapt to such relevant factors in Japan.
Conflict of Interests
No conflict of interests is known to the authors.
Acknowledgements
O. Johansson was supported by a grant from Mr. Einar Rasmussen, Kristiansand S, Norway, and by generous support from the Cancer and Allergy Foundation, Stockholm, Sweden and the Karolinska Institute. Mr. Brian Stein, Melton Mowbray, Leicestershire, UK, and the Irish Doctors Environmental Association (IDEA; Cumann Comhshaoil Dhoctúirí na hÉireann) are gratefully acknowledged for their general support to O. Johansson.
LidmanS.Kunskapens bok, term “Tobak”; Bokförlaget Natur och Kultur: Stockholm, Sweden, 1955KaldorJ.DayN.SchottenfeldD.FraumeniJ.Mathematical models in cancer epidemiology1996New York, NJ, USAOxford University Press127137WakaiK.AndoM.OzasaK.ItoY.SuzukiK.NishinoY.KuriyamaS. I.SekiN.KondoT.WatanabeY.OhnoY.TamakoshiA.MoriM.MotohashiY.TsujiI.NakamuraY.IsoH.MikamiH.InabaY.HoshiyamaY.SuzukiH.ShimizuH.ToyoshimaH.TokudomeS.HashimotoS.KikuchiS.KoizumiA.KawamuraT.MikiT.DateC.SakataK.NoseT.HayakawaN.YoshimuraT.ShibataA.OkamotoN.ShioH.KitagawaT.KurokiT.TajimaK.Updated information on risk factors for lung cancer: findings from the JACC study200515supplement 2S134S1392-s2.0-2624445481710.2188/jea.15.S134HallbergÖ.A reduced repair efficiency can explain increasing melanoma rates20081721471522-s2.0-3974919252010.1097/CEJ.0b013e3282b6fcebPerssonJ.SjöbergI.JohanssonS.-E.2004Örebro, SwedenStatistics SwedenThunM. J.HenleyS. J.CalleE. E.Tobacco use and cancer: an epidemiologic perspective for geneticists20022148730773252-s2.0-003715264910.1038/sj.onc.1205807Viewing the news from abroad2005281277577710.1093/jjco/28.12.775-aDollR.PetoR.BorehamJ.SutherlandI.Mortality in relation to smoking: 50 years' observations on male British doctors20043287455151915282-s2.0-3042700097International Agency for Research on Cancer, http://www-dep.iarc.frWHOTobacco Free Initiative, Economics and Statisticshttp://www.who.int/tobacco/economics/enOsmondC.Using age, period and cohort models to estimate future mortality rates19851411241292-s2.0-0021905690MøllerB.FekjærH.HakulinenT.TryggvadóttirL.StormH. H.TalbäckM.HaldorsenT.Prediction of cancer incidence in the Nordic countries up to the year 20202002111S1S962-s2.0-0036313711FlandersW. D.LallyC. A.ZhuB. P.HenleyS. J.ThunM. J.Lung cancer mortality in relation to age, duration of smoking, and daily cigarette consumption: results from cancer prevention study II20036319655665622-s2.0-0141884400SascoA. J.SecretanM. B.StraifK.Tobacco smoking and cancer: a brief review of recent epidemiological evidence200445supplement 2S3S92-s2.0-844423476210.1016/j.lungcan.2004.07.998ArbeevK. G.UkraintsevaS. V.ArbeevaL. S.YashinA. I.Mathematical models for human cancer incidence rates2005122372722-s2.0-26044469066UkraintsevaS. V.YashinA. I.Economic progress as cancer risk factor. II. Why is overall cancer risk higher in more developed countries?2005WP 2005-022Rostock, GermanyMax Planck Institute for Demographics ResearchWeiQ.ChengL.HongW. K.SpitzM. R.Reduced DNA repair capacity in lung cancer patients19965618410341072-s2.0-0029832418ChengL.SpitzM. R.HongW. K.WeiQ.Reduced expression levels of nucleotide excision repair genes in lung cancer: a case-control analysis2000218152715302-s2.0-0033854786El-ZeinR. A.MonroyC. M.CortesA.SpitzM. R.GreisingerA.EtzelC. J.Rapid method for determination of DNA repair capacity in human peripheral blood lymphocytes amongst smokers2010104392-s2.0-7795564865610.1186/1471-2407-10-439GandaraD. R.KawaguchiT.CrowleyJ.MoonJ.FuruseK.KawaharaM.TeramukaiS.OheY.KubotaK.WilliamsonS. K.GautschiO.LenzH. J.McLeodH. L.LaraP. N.ColtmanC. A.FukuokaM.SaijoN.FukushimaM.MackP. C.Japanese-US common-arm analysis of paclitaxel plus carboplatin in advanced non-small-cell lung cancer: a model for assessing population-related pharmacogenomics20092721354035462-s2.0-7024913764410.1200/JCO.2008.20.8793