| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
1 Biostatistics Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute and 2 Genetic Epidemiology Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Rockville, Maryland
Requests for reprints: Jay H. Lubin, Biostatistics Branch, Division of Cancer Epidemiology and Genetics, National Cancer Institute, 6120 Executive Boulevard, Rockville, MD 20852. E-mail: lubinj{at}mail.nih.gov
| Abstract |
|---|
|
|
|---|
20 cigarettes per day and decreases with intensity for subjects who smoke >20 cigarettes per day. The intensity patterns are homogeneous by histologic type of lung cancer, suggesting that observed differences in risks by histologic type are related to total smoking exposure or smoking duration and not smoking intensity. At lower smoking intensities, there is an "exposure enhancement" effect such that for equal total exposure, the excess OR per pack-year increases with intensity. At higher smoking intensities, there is a "reduced potency" or "wasted exposure" effect such that for equal total exposure, the excess OR per pack-year decreases with intensity (i.e., smoking at a lower intensity for longer duration is more deleterious than smoking at a higher intensity for shorter duration). (Cancer Epidemiol Biomarkers Prev 2006;15(3):51723) | Introduction |
|---|
|
|
|---|
Investigators typically assess the association between cancer occurrence and cigarette smoking using duration of smoking and smoking intensity, as measured by the number of cigarettes smoked per day (1, 2). Faced with the complexity of the association between duration and intensity and disease risk, investigators may apply nonparametric or semiparametric models, such as splines (3) or generalized additive models (4), or create a single comprehensive smoking index (5, 6). However, problems of interpreting multiple characteristics of an exposure remain (7, 8). Comparisons of intensity based on odds ratios (OR) for categories of cigarettes per day either adjusted for, or jointly with, duration of smoking are influenced by differences in total exposure to cigarette smoke. For example, to assess the effect of smoking intensity, one would typically compare ORs of smoking 30 years and smoking 40 years among individuals who smoke 20 cigarettes per day to the ORs of smoking 30 years and smoking 40 years among individuals who smoke 30 cigarettes per day. Differences in the patterns of these ORs for duration would then be ascribed to the effects of intensity. However, differences in total exposure influence this comparison; i.e., 30 and 40 pack-years of total exposure for the 20 cigarettes per day individuals, respectively, compared with 45 and 60 pack-years for the 30 cigarettes per day individuals.
Studies of smoking and cancers of the lung and the bladder show a leveling of risk for smokers consuming more than 20 to 40 cigarettes per day (9, 10). The precise implication of this pattern is, however, uncertain. Heavy smokers inhaling less deeply, exposure-dependent biases, or behavioral factors are possible reasons for the leveling of risk at higher intensities. Another plausible explanation is a nonlinear relationship of risk and intensity of exposure. Using data from a large case-control study of lung cancer, we apply a model for total pack-years of exposure, which enables a more direct assessment of smoking intensity. Our primary focus is the delivery of total exposure; i.e., the risk associated with total exposure delivered at low intensity for a long period of time compared with the risk with the same total exposure delivered at high intensity for a short period of time. Our interest involves how intensity of smoking influences the association between pack-years and disease or, correspondingly, to what extent the disease-exposure relationship is modified by smoking intensity.
| Materials and Methods |
|---|
|
|
|---|
The primary measure of total exposure is pack-years, which is computed as the product of duration of smoking and mean number of packs (20 cigarettes) smoked per day.
This study was done under approval of all applicable institutional review boards.
Models
We first apply a simple linear model for the OR of lung cancer and pack-years of exposure, d,
![]() | (A) |
![]() | (B) |
For nonsmokers, d = 0 and OR(0) = 1. The
parameter measures downward concavity (
< 0) or upward convexity (
> 0) with d, and thus the degree of departure from linearity. A test of the null hypothesis
= 0 is a test of no departure from the linear model in pack-years. We considered other variants for modeling nonlinearity, but Eq. (B) proved satisfactory.
In initial analyses, we find that model (A) fits the data poorly, unless intensity of smoking is included. We use two approaches to model intensity of smoking. The first approach defines I intensity categories and indicator variables, ni, i = 1, ..., I, where ni = 1 if a subject's intensity level occurs within the ith category; 0 otherwise. The model is
![]() | (C) |
1 set to zero for identifiability,
2,...,
I represent category-specific effects relative to the i = 1 level. However, our primary interest is ßexp(
i), which represents the EOR/pack-year for the ith intensity category. For convenience in the fitting, exp(ß*) replaces ß in Eq. (C) to avoid range restrictions on the ß parameter. With this reparametrization, SEs for the parameter estimates translate into multiplicative SEs on the original scale. The model is easily extended to other categorical factors. We assess variation in the EOR/pack-year parameter for continuous intensity (n) with
![]() | (D) |
1n +
2n2}; g(n) = exp{
1ln(n) +
2n}; and g(n) = exp{
1ln(n) +
2ln(n)2}. We select these forms primarily for their flexibility and convenience in fitting. All models stratify on center (seven levels), sex (two levels), and age (five levels: 50-54, 55-59, 60-64, 65-70, and 70-74). We use likelihood ratio tests for comparisons of nested models. For all modeling, we use the Epicure software package (13).
| Results |
|---|
|
|
|---|
40 cigarettes per day (OR = 18.4) is complicated by differences in total exposure to tobacco smoke.
|
|
|
i = 0 for model (B) applied simultaneously to all intensity categories, whereas only one of 12 individual tests of
i = 0 is rejected at the 0.05 level, P = 0.04). Estimates of EOR/pack-year increase for lower intensities, then decline at higher intensities (Fig. 3). We fit model (D) to the case-control data to characterize the effect of intensity on the relationship between pack-years and lung cancer risk (Table 1). Changes in deviance indicate that model M1 with g(n) = exp{
1log(n) +
2log(n)2} provides the best fit although its improvement in fit over model M2 is slight. The locus of points for model M1 closely follows the category-specific EOR/pack-year estimates (Fig. 3, solid line). The maximum EOR/pack-year based on model M1 occurs at exp(
1/2
2) or 17.6 cigarettes smoked per day. Below the maximum, there is a "direct exposure-rate effect"; i.e., lung cancer risk per pack-year of exposure increases with increasing intensity or, more specifically, a given total exposure imparted at a higher intensity is more harmful than the same exposure imparted at a lower intensity. Above the maximum, there is an "inverse exposure-rate effect"; i.e., for equal total exposure, the exposure-response is inversely related to intensity or, more specifically, total exposure imparted at a lower intensity for a longer period of time is more deleterious than the equivalent exposure imparted at a higher intensity for a shorter period of time. For categories 1-19, 20-29, 30-39, and 40+ cigarettes smoked per day, mean intensities are 13.0, 22.0, 32.6, and 47.3, respectively. The predicted EOR/pack-year based on model M1 closely corresponds to the observed patterns for the four intensity categories (Fig. 2, dashed lines).
|
Histologic Type of Lung Cancer and Exposure
We fit model M1 with cases restricted to a major histologic type, including squamous cell carcinoma, small-cell carcinoma, large-cell carcinoma, and adenocarcinoma using the same control group for each analysis (Table 2; Fig. 4, solid lines). The estimate of ß varies with histologic type; however, tests of homogeneity (Table 2, column 7) indicate that intensity parameters (
1 and
2) for each histologic type do not vary significantly and are consistent with estimates from all data combined (Fig. 4, dashed lines). Adjusted for intensity, EOR/pack-years are highest for cases with squamous cell carcinoma histology, followed by small-cell carcinoma, large-cell carcinoma, and adenocarcinoma histologies. For squamous cell carcinoma, small-cell carcinoma, large-cell carcinoma, and adenocarcinoma, fitted maxima occur at 16.4, 19.1, 13.9, and 19.1 cigarettes per day, respectively, which are similar to the 17.6 cigarettes per day for all data, resulting in an estimated maximum EOR/pack-year estimates of 0.58, 0.38, 0.31, and 0.09/pack-year for the histologic types.
|
|
1 and
2; P < 0.001). However, use of filter and nonfilter cigarettes varies considerably across study areas. After controlling for type of cigarette smoked and center, we find a significant interaction of the ß parameter and center (P < 0.001) but do not reject homogeneity of the intensity parameters by center [P = 0.15 for the test of no interaction of g(n) and center]. | Discussion |
|---|
|
|
|---|
20 cigarettes per day and decreases with intensity for subjects who smoke >20 cigarettes per day. At lower smoking intensities, the data support an "exposure enhancement" effect, such that for equal total exposure, the EOR/pack-year of smoking increases with intensity. At higher smoking intensities, data support a "reduced potency" or "wasted exposure" effect, such that for equal total exposure, smoking at a lower intensity for longer duration is more deleterious than smoking at a higher intensity for shorter duration. It is important to note that these patterns reflect effects of intensity and not lung cancer risk. For example, the inverse intensity effect implies that an increase in smoking intensity decreases risk per pack-year and does not imply a decrease in the overall risk of lung cancer, which depends on both total pack-years of exposure and smoking intensity. Results from epidemiologic studies indicate that risks of lung cancer, as well as bladder cancer, tend to level off with increased smoking intensity (10). The EOR/pack-year patterns observed in the European Smoking and Health Study data are consistent with these previous findings, at least at higher intensities, but precise comparisons are difficult because previous results do not control for total pack-years of exposure.
It remains to be determined whether and to what extent the patterns of variation of EOR/pack-year with intensity reflect nicotine dependency and its effect on smoking inhalation and smoking intensity or underlying biological processes, such as changes in activation and detoxification capacities for carcinogenic compounds in cigarette smoke or in DNA repair capacity, or biases in exposure assessment. Patterns of variation of the EOR/pack-year by intensity of smoking may reflect modulation of inhalation practices such that lower-intensity and higher-intensity smokers ingest relatively fewer carcinogens per cigarette smoked compared with moderate-intensity smokers. This would result in reduced risks at lower intensities and leveling or declining risks at higher intensities (10). Thus, variation in frequency or depth of inhalation is a possible explanation for observed intensity patterns. A recent study of 190 smokers found increased plasma cotinine and nicotine levels with increased cigarettes smoked per day, but a marginally significant (P = 0.08) decline in "nicotine boost"; i.e., the increase in blood plasma nicotine after smoking one cigarette (14). We can directly evaluate the association between inhalation practices and intensity in the European Smoking and Health Study data. Table 3 shows percentages of smoking controls who inhale moderately or deeply or who inhale most or all of the time by categories of cigarettes per day and total exposure. Subjects with greater total exposure are more likely to inhale their cigarettes deeply or frequently; however, within categories of pack-years, there is little indication that depth or frequency of inhalation is related to smoking intensity.
|
Polycyclic aromatic hydrocarbons and DNA and protein adducts have limitations in analyses of tobacco effects because they can arise from sources other than smoking. The compound 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) is a tobacco-specific carcinogen and levels of its metabolites directly represent markers of tobacco effects (20, 21). If the patterns of intensity effects in the European Smoking and Health Study data derive from biological processes, then we would expect a nonlinear relationship between NNK metabolites and urinary cotinine levels, a marker of tobacco exposure, with relatively less NNK metabolites produced at high, as opposed to low, cotinine levels. Whereas this analysis has not formally been carried out, this pattern can be observed in Fig. 5 of Carmella et al. (20).
|
The increasing EOR/pack-year at lower intensities and the decreasing EOR/pack-year at higher intensities could be due, at least in part, to misclassification of cigarettes smoked per day. Observed variations of the EOR/pack-year with intensity would require a very specific pattern of misclassification; i.e., a progressively increasing amount of misclassification with increasing intensity (above 15 to 20 cigarettes smoked per day), resulting in an increasing bias towards the null and a decreasing exposure-response relationship, and a decreasing amount of misclassification of intensity up to 15 to 20 cigarettes smoked per day, resulting in decreasing bias towards the null and an increasing exposure-response relationship. The latter pattern is plausible if an increasing intensity results in a smoker consuming a proportionally smaller fraction of each cigarette. Nonetheless, if low and high intensities were differentially misclassified compared with intermediate intensities, then pack-years of smoking would reflect the differential misclassification because we expect that duration of smoking is relatively accurate. We would expect to observe a progressively more concave relationship between the ORs for lung cancer and pack-years with increasing intensity and with decreasing intensity. However, exposure-response relationships for ORs and pack-years are consistent with linearity over the full range of intensities, suggesting that any differential misclassification by intensity is minimal and not likely to induce the observed patterns.
The best fitting model in Table 1 (model M1) involves an intensity function of the form g(n) = exp{
1ln(n) +
2ln(n)2}. Because pack-years, d, is yxn and OR = 1 + ßdxg(n), model M1 can be rewritten in terms of duration and intensity as OR = ßyxexp{
ln(n) +
2ln(n)2}, where
=
1 + 1. Thus, model M1 with the above g(n) is equivalent to a model which is linear in duration for fixed intensity.
It is well recognized that the association between smoking and lung cancer varies with histologic type (1). Our analyses suggest that whereas levels of lung cancer risk vary by histologic classification, effects of smoking intensity (i.e., the curvatures) do not vary by histologic type. This suggests that stochastic factors which influence pathways that predispose towards a specific histology are not associated with smoking intensity but more closely linked with total exposure or smoking duration. The precise implication of this observation is unclear but may suggest that determinants of histologic type are more likely linked to early transformational processes.
The Armitage-Doll multistage model for carcinogenesis is based on the observation that the logarithm of the lung cancer rate increases approximately linearly with the logarithm of age and is consistent with a multistage model for carcinogenesis, whereby a normal epithelial cell undergoes multiple, sequential, heritable, and nonreversible transformations to become malignant. Doll and Peto (28) used this model in their analysis of the British Doctors' Study and found that lung cancer rates in continuing smokers increased with the square of intensity plus 6 and the 4.5 power of duration minus 3.5. In an analysis of lung cancer mortality in the American Cancer Society's Cancer Prevention Study I, a cohort study of 1,078,894 subjects, Knoke et al. (29) found significant improvement in model fit, in comparison with the Doll and Peto model, by including a factor for either age at smoking initiation or attained age. The model for lung cancer mortality rate in smokers was rate =
(n + 6)ß x (y 3.5)
x f
, where n is number cigarettes smoked per day, y is duration of smoking, and f is a factor representing either age at start of smoking or attained age (models C or D, respectively, in ref. 29), and where
, ß,
and
are parameters. There was a separate model for lung cancer mortality rate in nonsmokers with the form: rate =
(age 3.5)ß. As in the Doll and Peto analysis, Knoke et al. added 6 cigarettes per day to account for the background risk in the absence of smoking and added 3.5 years to reflect a time lag from the appearance of a malignant cell to lung cancer mortality. To compare our model M1 to Knoke's models, we assume our OR model for lung cancer incidence approximates the relative risk for lung cancer mortality, and multiply model M1 by the Cancer Prevention Study I lung cancer mortality rate model for nonsmokers to obtain lung cancer (mortality) rates. For the comparison, we assume individuals start smoking at age 19 years. Figure 5 shows that predicted lung cancer mortality for males who smoke 10, 20, or 30 cigarettes per day based on our exposure and exposure-rate model are very similar to predictions based on models C and D in Knoke et al.
Finally, patterns of intensity effects are subject to substantial uncertainty, particularly at lower intensities. This uncertainty is due to the limited range of pack-years of exposure and an increased variability in estimating EOR/pack-years. For example, the median and interquartile range for total exposure are 6.7 pack-years and 4.9 to 8.4 pack-years, respectively, for subjects smoking under 5 cigarettes per day, and 14.4 pack-years and 11.4 to 18.2 pack-years for subjects smoking 5 to 10 cigarettes per day, compared with 40 pack-years and 28 to 53.8 pack-years for all smokers. Thus, any conclusion about whether the estimated EOR/pack-year approaches zero remains constant or, indeed, increases for intensities below 5 to 8 cigarettes per day must be made cautiously. Results are subject to the additional uncertainty from use of the overall mean number of cigarettes smoked per day, rather than accounting for changes in smoking intensity throughout life.
In summary, our analysis of a large case-control study of lung cancer using a novel exposure and exposure-rate model reveals a direct intensity effect at low smoking intensities; i.e., an intensity enhancement effect, resulting in an increasing EOR/pack-year, and an inverse intensity-rate effect (i.e., reduced potency or wasted exposure effect), resulting in a decreasing EOR/pack-year at higher intensities. Our modeling approach is applicable in other epidemiologic settings where data on exposure and exposure rate are available.
| Acknowledgments |
|---|
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 11/ 8/05; revised 12/28/05; accepted 1/18/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N. Pandeya, G. M. Williams, S. Sadhegi, A. C. Green, P. M. Webb, and D. C. Whiteman Associations of Duration, Intensity, and Quantity of Smoking with Adenocarcinoma and Squamous Cell Carcinoma of the Esophagus Am. J. Epidemiol., May 15, 2008; (2008) kwn091v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. H. Lubin, J. Virtamo, S. J. Weinstein, and D. Albanes Cigarette Smoking and Cancer: Intensity Patterns in the Alpha-Tocopherol, Beta-Carotene Cancer Prevention Study in Finnish Men Am. J. Epidemiol., April 15, 2008; 167(8): 970 - 975. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. H. Lubin, N. Caporaso, D. K. Hatsukami, A. M. Joseph, and S. S. Hecht The Association of a Tobacco-Specific Biomarker and Cigarette Consumption and Its Dependence on Host Characteristics Cancer Epidemiol. Biomarkers Prev., September 1, 2007; 16(9): 1852 - 1857. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. H. Lubin, M. C. R. Alavanja, N. Caporaso, L. M. Brown, R. C. Brownson, R. W. Field, M. Garcia-Closas, P. Hartge, M. Hauptmann, R. B. Hayes, et al. Cigarette Smoking and Cancer Risk: Modeling Total Exposure and Intensity Am. J. Epidemiol., August 15, 2007; 166(4): 479 - 489. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. A. Melikian, M. V. Djordjevic, S. Chen, J. Richie Jr., and S. D. Stellman Effect of Delivered Dosage of Cigarette Smoke Toxins on the Levels of Urinary Biomarkers of Exposure Cancer Epidemiol. Biomarkers Prev., July 1, 2007; 16(7): 1408 - 1415. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. H Lubin, M. Kogevinas, D. Silverman, N. Malats, M. Garcia-Closas, A. Tardon, D. W Hein, R. Garcia-Closas, C. Serra, M. Dosemeci, et al. Evidence for an intensity-dependent interaction of NAT2 acetylation genotype and cigarette smoking in the Spanish Bladder Cancer Study Int. J. Epidemiol., February 1, 2007; 36(1): 236 - 241. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Cell Growth & Differentiation |