Background:

Despite smoking being a well-established risk factor for pancreatic cancer, there is a need to further characterize pancreatic cancer risk according to lifespan smoking patterns and other smoking features, such as tobacco type. Our aim was to deeply investigate them within a large European case–control study.

Methods:

Tobacco smoking habits and other relevant information were obtained from 2,009 cases and 1,532 controls recruited in the PanGenEU study using standardized tools. Multivariate logistic regression analysis was performed to evaluate pancreatic cancer risk by smoking characteristics and interactions with other pancreatic cancer risk factors. Fractional polynomials and restricted cubic splines were used to test for nonlinearity of the dose–response relationships and to analyze their shape.

Results:

Relative to never-smokers, current smokers [OR = 1.72; 95% confidence interval (95% CI), 1.39–2.12], those inhaling into the throat (OR = 1.48; 95% CI, 1.11–1.99) or chest (OR = 1.33; 95% CI, 1.12–1.58), and those using nonfiltered cigarettes (OR = 1.69; 95% CI, 1.10–2.61), were all at an increased pancreatic cancer risk. Pancreatic cancer risk was highest in current black tobacco smokers (OR = 2.09; 95% CI, 1.31–3.41), followed by blond tobacco smokers (OR = 1.43; 95% CI, 1.01–2.04). Childhood exposure to tobacco smoke relative to parental smoking was also associated with increased pancreatic cancer risk (OR = 1.24; 95% CI, 1.03–1.49). Dose–response relationships for smoking duration, intensity, cumulative dose, and smoking cessation were nonlinear and showed different shapes by tobacco type. Effect modification by family history of pancreatic cancer and diabetes was likely.

Conclusions:

This study reveals differences in pancreatic cancer risk by tobacco type and other habit characteristics, as well as nonlinear risk associations.

Impact:

This characterization of smoking-related pancreatic cancer risk profiles may help in defining pancreatic cancer high-risk populations.

Pancreatic cancer is one of the deadliest cancer types worldwide (5-year relative survival rate in the range 5%–15%; ref. 1). Disastrously, estimates of pancreatic cancer incidence are increasing both in United States and Europe (2). Despite the etiology of pancreatic cancer is relatively unknown, it is estimated that 10% to 30% of all pancreatic cancer cases are caused by smoking (3). Prevention of smoking is therefore a strong measure to reduce the burden of pancreatic cancer in the population.

Although the association between smoking and pancreatic cancer is well-established, a detailed characterization of tobacco smoking habits in relation to pancreatic cancer risk is still lacking. A meta-analysis including 10,490 cases and 526,813 controls showed that being a current smoker, jointly with a longer smoking duration and a higher smoking intensity, were associated with an increase in pancreatic cancer risk (4). However, the authors assumed a linear trend for pancreatic cancer risk associated with increasing smoking exposure, a fact that was disputed by Zou and colleagues in an updated analysis combining 9,044 cases and 32,039 controls that showed a nonlinear dose–response relationship between several smoking characteristics and pancreatic cancer risk (5). In addition, the pooled analysis within the Pancreatic Cancer Case–Control Consortium (PanC4), including 6,507 cases and 12,890 controls, indicated that after a certain amount of smoking exposure pancreatic cancer risk levelled-off (6), shedding a different perspective on the dose–response relationship of smoking in relation to pancreatic cancer risk. However, in the aforementioned studies, an exploration of the shape of the association between smoking measures and pancreatic cancer risk was not further pursued. The shape of the dose–response relationship between cigarette smoking and pancreatic cancer risk was investigated in a recent meta-analysis of 38 case–control and 40 cohort studies (7). Risk patterns of pancreatic cancer in current versus smokers were compared in this study for smoking intensity and duration, ignoring the contribution to risk of former smokers. To understand multidimensional aspects of smoking in pancreatic cancer etiology, there is a need to provide consistent risk estimates for all smoking groups and to address the mutual influence of smoking intensity and duration.

Moreover, several aspects of tobacco smoking habits have not been considered until now. For instance, differences in pancreatic cancer risk by either black or blond tobacco use have not been explored despite the presumed differences in their chemical composition and damaging effects (8, 9). In fact, several studies have shown that black tobacco is associated with higher risk of bladder (8), colorectal (10), esophageal (11), and head-and-neck cancer (12, 13), than blond tobacco.

Therefore, we set out to investigate the association and dose–response relationship between tobacco smoking and pancreatic cancer risk in a large European population, considering every aspect of the smoking habit including use of black versus blond tobacco.

Study design and participants

The PanGenEU is an ongoing multicenter case–control study initiated in 2007, recruiting participants from six European countries (Germany, Ireland, Italy, Spain, Sweden, and United Kingdom) across 28 centers. Newly diagnosed patients with pancreatic cancer >18 years old and controls matched by age (±10 years), gender, and geographical area were included if they had lived in the study area for at least 6 months. A rapid ascertainment approach was applied: pancreatic cancer cases with a suspicion of the disease were recruited and remained in the study if the diagnosis was verified by the treating physician. Controls, sex-, age-, and centre-individually matched to cases, were mostly hospital-based and eligible if principal diagnosis at admission was unrelated to known risk factors of pancreatic cancer. Conditions of admission of controls are reported in Supplementary Materials and Methods. Population-based controls (Sweden and Ireland) were eligible if history of pancreatic cancer was absent. Participation rates were 86.3% for cases and 77.8% for controls. The study was approved by the IRB of all participant centers and all subjects gave written informed consent. More details of the study are provided elsewhere (14, 15).

Variables

Personal interviews to the study subjects were conducted by trained monitors using standardized protocols and questionnaires to obtain detailed information on lifetime smoking habits, among other pancreatic cancer risk factors. The smoking status of the participants was categorized into never-smokers if they smoked <100 cigarettes during their lifetime; occasional smokers if they smoked ≥1 cigarette/day for ≥6 months; former smokers if they quitted smoking for >1 year; and current smokers otherwise (>100 cigarettes during lifetime without permanent smoking cessation). Information on smoking habits by tobacco type (only black, blond, or both) was only collected in the Spanish centers. Smoking exposure was further assessed by the age at smoking initiation (years), age when last smoked (years), cigarettes/cigar/pipe-use (yes, no), the amount of cigarettes/cigars/pipes smoked in units of time (days, months, years), depth of inhalation (mouth, throat, chest), filter-use (filtered cigarettes, nonfiltered, both), and smoking status of the parents (never- or ever-smoker). From these characteristics, data on smoking duration (years), smoking intensity for cigarettes (per day) and cigars/pipes (per week), and time since cessation (years) was derived. Number of pack-years, representing cumulative dose, was calculated as (cigarettes per day/20) × smoking duration in years. Smoking variables by use of tobacco type were generated likewise. Environmental tobacco smoke (ETS) exposure during childhood was categorized according to the smoking status of the parents (none, one, or both).

Statistical analysis

Imputation of missing values, assumed to be at random, was performed using the Random Forest algorithm (R package missForest). Predictor variables such as center, country, and case–control status were kept in the imputation set. The performance of the imputation (Supplementary Table S1) was assessed by calculating the out of bag mean square error (OOB), representing the mean of squared differences between each observed value and its prediction, based on random forest trees (n = 100 was applied). The average OOB for all smoking variables was 5.27, with categorical variables presenting a markedly lower estimate (OOB = 0.04), indicating a better imputation performance of the latter. Use of unimputed data of all continuous variables, for which the proportion of missing values was relatively low (6.7%), was therefore deemed more appropriate for dose–response analyses. The performance of the imputation was also assessed with concordance rates between the observed and imputed data, considering a test dataset consisting of only subjects with complete data and missing values introduced by following the missingness rates of the original data. The concordance of all categorical variables was 94.4%.

Differences between cases and controls regarding smoking characteristics were evaluated by χ2 and Student t test (or Kruskal–Wallis test, where appropriate). Unconditional logistic regression analysis was performed to estimate OR and 95% confidence intervals (95% CI). Never-smokers were chosen as the reference category, except for the variables “age when last smoked” and “time since smoking cessation,” where current smokers were taken as the reference. Tertiles were created for the continuous variables based on the distribution of controls. A P value for trend was calculated by assuming ordinal variables in linear regression models. Age (≤54, 55–64, 65–74, ≥75 years), gender, and country-adjusted models (aOR) were considered (Model 1). For the tobacco type-specific analyses within the Spanish PanGenEU study population, the same model was applied, but replacing country by region (East, Central, and Northern Spain). The attributable risk (AR) of smoking in relation to pancreatic cancer (population exposed: 59%) was calculated from the fitted multivariate adjusted logistic regression models (R packages attribrisk and epiR). Because heterogeneity by country (P < 0.05) was evident for all smoking variables (e.g., P value for interaction by smoking status = 0.007; Supplementary Fig. S1), random effects for country were applied in mixed effects models. Because of absence of heterogeneity by region in the Spanish study population, logistic regression models without random effects were considered.

The influence of confounding factors or effect modification on the association was assessed for several variables: gender (female, male), age (<65 years, ≥65 years), obesity (body mass index >30: yes, no), diabetes (no, yes less than 2 years, yes more than 2 years), asthma (yes, no), chronic pancreatitis (yes, no), alcohol status/consumption (never, former, current), presence of periodontitis (yes, no) and recession (yes, no), educational level as a proxy for socioeconomic status (low, medium, high), and family history of pancreatic cancer (yes, no). Variables changing estimates by more than 10% or having a significant influence in the model (diabetes and family history of pancreatic cancer in some smoking-related variables) were considered as potential confounders. The le Cessie-van Houwelingen-Copas-Hosmer unweighted sum of squares test indicated a high goodness of fit of the models (16). Effect modification was assessed in interaction and stratified analyses. Additive interaction by time-related variables such as smoking duration was also evaluated by the relative excess risk due to interaction (RERI) and Delta-method CIs (17, 18).

To test for interaction, a likelihood ratio (LR) test was performed comparing models with and without an interaction term between the smoking variables and the covariates (e.g., age, gender, BMI and obesity, diabetes, asthma, alcohol, periodontitis, recession, educational level, and family history of pancreatic cancer). Effect modification was tested further via stratified analyses. To assess interaction by time-related variables, we explored the combined effect of smoking duration and other smoking characteristics such as tobacco type on pancreatic cancer risk. Smoking duration was categorized into <20, 20 to 30, and ≥30 years of smoking and stratified further by tobacco type considering never-smokers as the reference category.

To assess the dose–response relationships, pancreatic cancer risk estimates were calculated per 1-unit of increase in continuous smoking exposure variables considering linear and nonlinear models if so indicated by fractional polynomials (R package mfp; ref. 19). In addition, restricted cubic splines were used to confirm nonlinear associations and for modelling the shape of the dose–response relationships (R package splines; ref. 20). Nonlinearity of the models was tested via the likelihood-ratio test comparing the model with and without restricted cubic splines. Knots were set at the 10%, 50%, and 90% percentile of the exposure distribution, as comparable results were obtained with five knots (21).

Sensitivity analyses were performed comparing the risk estimates in magnitude and trend regarding the unimputed and imputed data, and the PanGenEU study population with and without Italy (since Italy provided cases only). As information bias could be induced by neglecting the quantity of smoking exposure, adjustment for cumulative dose (pack-years) was considered, thereby accounting for both smoking duration and smoking intensity. Additional adjustments were made also for smoking intensity and duration separately, to assess both the individual and joined effects of smoking characteristics independent of smoking duration and/or intensity. These adjustment variables were considered on the continuous scale and modeled as fractional polynomials to account for nonlinear effects. To further assess the performance of the restricted cubic splines, additional smoothing was applied by varying the degrees of freedom, allowing more flexibility into the model (22).

The threshold for statistical significance in two-sided tests was set at P < 0.05. Data were analyzed with R-project (version 3.4.1; ref. 23).

Table 1 shows the characteristics of the 2,009 cases and 1,532 controls included in this analysis. The Spanish centers contributed the most to both cases (N = 876) and controls (N = 762). Pancreatic cancer cases presented more frequently with a family history of pancreatic cancer and a diagnosis of diabetes or chronic pancreatitis.

Table 2 shows pancreatic cancer risk associated with smoking characteristics. The prevalence of smoking was higher in cases (27.4%) than in controls (17.6%), with a corresponding aOR of 1.72 (95% CI, 1.39–2.12) for current smokers compared with never-smokers. Furthermore, a statistically significant increased trend (P < 0.001) in pancreatic cancer risk was observed for longer smoking duration, higher smoking intensity and higher cumulative dose. The use of nonfiltered cigarettes increased risk of pancreatic cancer more prominently (aOR = 1.69; 95% CI, 1.10–2.61), although use of filtered cigarettes was also associated with an increased pancreatic cancer risk (aOR = 1.25; 95% CI, 1.06–1.48). Marked increases in pancreatic cancer risk were also observed for inhalation into the throat (aOR = 1.48; 95% CI, 1.11–1.99) and chest (aOR = 1.33; 95% CI, 1.12–1.58). Childhood exposure to ETS by smoking parents (vs. nonparental exposure) was also associated with a 24% (95% CI, 1.03–1.49) increased pancreatic cancer risk. Risk for former smokers decreased progressively with longer time since smoking cessation when compared with current smokers (aOR for 14–28 years after cessation = 0.67; 95% CI, 0.51–0.88). A negative trend of the risk was also observed if compared with never-smokers (pancreatic cancer risk was diminished from 14 years of cessation), and when considering smoking cessation time at 5-year intervals (Supplementary Table S2). No significant associations between pancreatic cancer risk and pipe/cigar-use or other smoking variables were observed (Supplementary Table S3). Additional adjustment for diabetes and family history of pancreatic cancer led to minimal differences in risk estimates (Supplementary Table S3). Effect modification was apparent only for family history of pancreatic cancer and diabetes status (Supplementary Table S4), pointing towards a higher pancreatic cancer risk among current smokers with family history of the disease (aOR = 2.24; 95% CI, 0.66–7.61) and former smokers with diabetes (aOR = 1.44; 95% CI, 0.91–2.28; P value for interaction <0.001).

Table 3 shows pancreatic cancer risk estimates by tobacco type in PanGenEU-Spain. Compared with never-smokers, pancreatic cancer risk was significantly increased for smokers of only black tobacco (aOR = 1.55; 95% CI, 1.13–2.12) and of both tobacco types (aOR = 1.58; 95% CI, 1.14–2.17). Considering smokers of only blond tobacco, pancreatic cancer risk tended to be increased (aOR = 1.23; 95% CI, 0.94–1.62), though without reaching statistical significance. When further stratifying by smoking status, a significant increase in risk was observed for current smokers of only black tobacco (aOR = 2.09; 95% CI, 1.31–3.41) and blond tobacco (aOR = 1.43; 95% CI, 1.01–2.04). Former smokers of only black tobacco were at increased, though milder, pancreatic cancer risk (aOR = 1.40; 95% CI, 0.98–1.99).

Table 4 shows the combined effect of smoking duration and type of tobacco on pancreatic cancer risk. Compared with never-smokers, smoking for ≥30 years of both tobacco types was associated with a higher pancreatic cancer risk than smoking only black or blond tobacco (aOR = 2.05; 95% CI, 1.25–3.36 and RERI = 0.206; 95% CI, –0.49–0.91).

Table 5 shows risk estimates for continuous smoking variables associated with pancreatic cancer. Nonlinear associations were evident for smoking duration and intensity, cumulative dose, time since cessation, and age at smoking initiation. Adjusted fractional polynomials models suggested a statistically significantly higher pancreatic cancer risk per 1-unit increase in smoking duration, smoking intensity and cumulative dose, and decreasing pancreatic cancer risks for age at smoking initiation and time since smoking cessation. Linear associations were observed for other variables such as intensity of smoking cigars/pipes (data not shown). The restricted cubic splines approximating the shape of the dose–response relationships confirmed these nonlinear associations. Compared with never-smokers, smoking for >25 years (Fig. 1A and B) and smoking >20 cigarettes/day (Fig. 1D and E) was associated with a statistically significant increase of pancreatic cancer risk. Similarly, a cumulative dose of >14 pack-years was associated with increased pancreatic cancer risk (Fig. 1C). Visual inspection for smoking intensity and cumulative dose was suggestive of plateauing of pancreatic cancer risk, at approximately 27 cigarettes/day or pack-years. Concerning time since smoking cessation (Fig. 1FI), and relative to current smokers, risk appeared to decrease between 8 and 11 years of cessation and after around 18 years of cessation regardless of cumulative dose. In between these periods, the significant effect disappeared. By tobacco type, corresponding periods of significant decrease in pancreatic cancer risk were observed for black tobacco (after about 14 years since cessation) and for blond tobacco (between 2 and 8 years and after >20 years of cessation).

No relevant differences in the trend or magnitude of the estimates were found in sensitivity analyses (Supplementary Tables S3 and S5–S7), including further smoothing of the splines fit (Supplementary Fig. S2). In analyses adjusting for smoking intensity, risk estimates decreased in magnitude but showed a similar trend. By tobacco type, this adjustment did not affect either the associations nor the shapes of the relationships despite black tobacco smokers smoked heavier and for a longer time (Supplementary Table S8). Importantly, adjustment for smoking duration led to statistically nonsignificant risk estimates and change in the shape of the dose–response relationships (Supplementary Table S9). Joint effect analyses of smoking intensity and duration showed that long-lasting smoking together with intense smoking increase pancreatic cancer risk, whereas for less intense smoking the association weakened (Supplementary Table S10).

This study confirms that, in comparison to never-smokers, being a current smoker increases the risk of pancreatic cancer by 72%. In terms of AR, this study also endorses that around 16% (95% CI, 9.24–22.47) of all pancreatic cancer diagnoses could be avoided through tobacco preventive measures. A more detailed examination of smoking characteristics showed that the use of nonfiltered cigarettes, deep inhalation into the throat or chest, and exposure to tobacco smoke in the parental household were all associated with increased pancreatic cancer risk. Pancreatic cancer risk in black tobacco smokers was significantly higher compared with never-smokers, with blond tobacco smokers showing a less prominent risk pattern. Analysis of dose–response relationships corroborated that a higher smoking intensity, longer smoking duration, and increased levels of cumulative dose were associated further with an increased pancreatic cancer risk, whereas smoking cessation led to a gradual decline in pancreatic cancer risk, all in a nonlinear manner.

Our results are concordant with earlier studies on the same topic. Regarding the magnitude of pancreatic cancer risk associated with current versus never tobacco smoking, a meta-analysis and pooled analyses from the PanC4 and the Pancreatic Cancer Cohort Consortium showed similar estimates (RR = 1.74; 95% CI, 1.61–1.87; OR = 2.20; 95% CI, 1.71–2.83; and OR = 1.77; 95% CI, 1.38–2.26, respectively; refs. 4, 6, 24). Similarly, our study confirmed the trends and timing of tobacco smoking (4, 6), the excess risk conferred by tobacco smoking (4, 25, 26), the nonlinear tobacco– pancreatic cancer associations (5, 7), and risk due to childhood ETS (27). Compared with studies restricting ETS exposure to never-smokers, we also did not observe significant risk estimates (aOR = 1.24; 95% CI, 0.95–1.63; refs. 28, 29), suggesting that smokers, possibly more likely being exposed to childhood ETS, were driving this association in the overall analyses (aOR for current smokers exposed to parental smoking versus never smoking exposure = 2.01; 95% CI, 1.50–2.69). In contrast to the positive association between current cigar/pipe smokers and pancreatic cancer risk reported before (4, 30), we did not observe a significant associations in our study, probably due to low statistical power.

Effect-modification factors

The higher pancreatic cancer risk among smokers with family history of pancreatic cancer was previously described in our study population (14). Although statistical significance was not reached, former smoking diabetes patients tended to have a higher pancreatic cancer risk too. Were this true, lifestyle changes among diabetic patients including smoking cessation, which in turn may lead to weight gain and insulin resistance (31, 32), might explain this finding. Previous studies suggested differences in smoking effects on pancreatic cancer risk by gender (5, 6), although they failed to demonstrate effect modification by this variable. Similarly, nonsignificant differences by gender were found in our study, which included a large female sample with a relatively high smoking prevalence.

Dose–response relationships

Nonlinear relationships of the association between smoking variables and pancreatic cancer risk were supported by both fractional polynomials and restricted cubic splines approaches. Because fractional polynomials in regression models become imprecise with small sample sizes (22), we based the dose–response curves on results derived from restricted cubic splines, which allow a more flexible modeling (33). Concordant with the observation of nonlinear associations for smoking duration, intensity, and cumulative dose, a plateuing in the dose–response relationship was apparent. This observed pattern was previously reported (5, 7), and could be attributed to the saturation of the detoxification processes of tobacco carcinogens in the body (34), or to a presumably weaker inhalation of tobacco smoke but stronger DNA repair efficiency among heavy smokers (35, 36), among other factors. Nonlinear associations for smoking cessation, with decreased pancreatic cancer risk after 20 years of smoking cessation, were also suggested (5) and confirmed by other studies (7). However, in these earlier studies, consideration was not given to the influence of smoking intensity and duration on these associations. Patterns in pancreatic cancer risk in our study changed after adjusting for smoking duration mainly, whereby the magnitude of the risk estimates was affected (Supplementary Table S9).

Black versus blond tobacco use

Compared with never-smokers, black tobacco smokers showed a significantly higher pancreatic cancer risk, this tobacco type appearing to be more harmful than blond tobacco. This result is consistent with the few studies that examined the association between smoking by tobacco type in bladder (8, 37, 38) and other cancers (10–13). Smoking both black and blond tobacco for a long time (≥30 years) tended to be related to higher pancreatic cancer risk, this also being shown in previous studies on tobacco smoking and bladder cancer (8).

The difference between the two tobacco types could be explained by their smoke composition: black tobacco mostly contains early-stage carcinogens, such as N-nitrosamines and aromatic amines including 4-amino-biphenyl and 2-naphthylamine (39), whereas blond tobacco may mostly consist of late-stage carcinogens (37). It is conceivable that the two tobacco types contribute to pancreas carcinogenesis through different mechanisms: black tobacco may predominantly cause DNA mutations whereas blond tobacco may preferentially act through epigenetic changes, as has been shown for LINE-1 (9). As a consequence, an immediate and significantly higher increase in pancreatic cancer risk could be expected in black tobacco smokers, whereas blond tobacco might need a longer time to trigger pancreatic cancer development. This may also imply that following smoking cessation of blond tobacco pancreatic cancer risk can keep increasing for some time, slowing down after recovery of certain DNA methylation changes. In fact, methylation changes due to smoking seem to persist up to 22 years after smoking cessation (40). For black tobacco, instead, the pancreatic cancer risk reduction effects might not take place or might require longer since smoking cessation. Our results support these hypotheses to some extent. Compared with never-smokers, not only did black tobacco smoking have a more detrimental effect on pancreatic cancer risk, but also the risk tended to increase soon after smoking initiation, whereas downward risks were observed after smoking cessation for >10 years. A similar decreasing risk with long-term smoking cessation of black tobacco has been observed in bladder cancer in some (37, 38), but not all (8), studies. Among blond tobacco smokers, the trend towards a reduction in pancreatic cancer risk became evident shortly after smoking cessation (Supplementary Table S9). The shape of dose–response curves supported the aforementioned trends, specifically regarding smoking cessation. Thus, our study suggests that black tobacco consumption may play a role in several steps of the carcinogenic process with possibly both early and late-stage carcinogens being involved. For blond tobacco, our results point to a two-tier mechanism after smoking cessation driven by late-stage carcinogens, the first consisting of a sudden change in risk estimates with risk levels more akin to never-smokers likely due to desaturation of detoxification routes of tobacco-carcinogens (5, 41), the second showing risks leveling-off after approximately 20 years of smoking cessation, once alteration of DNA methylation levels of key genes regain the state of normalcy.

Among the limitations of the study, stratifying by tobacco type might have underpowered the analyses to detect any differences. As in any other study, subgroup analyses and multiple statistical tests are prone to chance findings due to increased type I error. Also, we could not consider potential differences in the content of carcinogens because we lacked information on tobacco brands, likely to contain varying amounts of heavy metals (42) and other carcinogens (39). Residual confounding can be therefore expected, also due to lack of, or imprecise, information on other relevant data such as ETS in adulthood. Extensive efforts have been made to adjust for as much confounding as possible, thereby alleviating the bias to the highest extent possible. Moreover, differential misclassification of the exposure due to recall bias of smoking habits among either the cases or controls is possible, or because use of only black or blond tobacco smoking might not have been reliably reported. Therefore, mixed effects due to alternate use of both tobacco types cannot be ruled out. We considered only smokers of black or blond tobacco to keep the effects by tobacco type separate, and considered switching from one type to the other in the group of users of both tobacco types.

Major strengths of the study are the large number of pancreatic cancer cases representing a European-wide pancreatic cancer population and the degree of detail in the information collected about smoking habits. This allowed us to undertake exhaustive and solid analyses considering many aspects of the habit in relation to pancreatic cancer risk. In fact, this is the first study assessing pancreatic cancer risk by black and blond tobacco. Also, as a novelty, the shapes of dose–response relationships have been fully characterized using different modeling strategies to account for nonlinear effects of smoking on pancreatic cancer risk.

In conclusion, findings of this study support and add to the previous evidence that smoking increases pancreatic cancer risk and demonstrates, for the first time, that both blond and black tobacco smoke are key in pancreatic cancer etiology, although probably acting through different genetic mechanisms. Considering these smoking-related pancreatic cancer risk profiles may help to refine the definition of high-risk pancreatic cancer population towards screening interventions implementation. Future studies should confirm our findings on type of tobacco and shed light on the mechanisms underlying their differential association with pancreatic cancer risk.

M. Hidalgo is a consultant for Agenus, Pharmacy, InxMed, Tolero, and Takeda; reports receiving a commercial research grant from Bioline; and has ownership interest (including patents) in Agenus, Pharmacy, and BioOncotech. No potential conflicts of interest were disclosed by the other authors.

Conception and design: E. Molina-Montes, L. Sharp, M. Márquez, N. Malats

Development of methodology: E. Molina-Montes, L. Sharp, N. Malats

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): M. Löhr, L. Sharp, X. Molero, M. Márquez, C.W. Michalski, A. Farré, J. Perea, M. O'Rorke, W. Greenhalf, L. Ilzarbe, A. Tardon, T.M. Gress, V.M. Barberà, T. Crnogorac-Jurcevic, L. Muñoz-Bellvis, E. Domínguez-Muñoz, J. Balsells, E. Costello, M. Iglesias, J. Kleeff, B. Kong, J. Mora, D. O'Driscoll, I. Poves, A. Scarpa, J. Yu, W. Ye, M. Hidalgo, A. Carrato, R. Lawlor, N. Malats

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): E. Molina-Montes, L. Van Hoogstraten, P. Gomez-Rubio, X. Molero, M. Hidalgo, F.X. Real, N. Malats

Writing, review, and/or revision of the manuscript: E. Molina-Montes, L. Van Hoogstraten, P. Gomez-Rubio, M. Löhr, L. Sharp, C.W. Michalski, J. Perea, M. O'Rorke, W. Greenhalf, L. Ilzarbe, A. Tardon, T.M. Gress, V.M. Barberà, T. Crnogorac-Jurcevic, L. Muñoz-Bellvis, E. Domínguez-Muñoz, J. Balsells, E. Costello, J. Kleeff, I. Poves, A. Scarpa, W. Ye, M. Hidalgo, A. Carrato, F.X. Real, N. Malats

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): E. Molina-Montes, X. Molero, M. Márquez, D. O'Driscoll

Study supervision: E. Molina-Montes, X. Molero, M. Márquez, D. O'Driscoll, N. Malats

Other (acquisition of funding): L. Sharp

Other (study design): F.X. Real

The full list of PanGenEU centers and investigators is available in the supplementary materials for this article. The work was partially supported by Fondo de Investigaciones Sanitarias (FIS), Instituto de Salud Carlos III, Spain (#PI11/01542, #PI0902102, #PI12/01635, #PI12/00815, #PI15/01573); Red Temática de Investigación Cooperativa en Cáncer, Spain (#RD12/0036/0034, #RD12/0036/0050, RD12/0036/0073); European Cooperation in Science and Technology – COST Action #BM1204: EUPancreas EU-6FP Integrated Project (#018771-MOLDIAG-PACA), EU-FP7-HEALTH (#259737-CANCERALIA, #256974-EPC-TM-Net); Associazione Italiana Ricerca sul Cancro (12182); Cancer Focus Northern Ireland and Department for Employment and Learning; and ALF (#SLL20130022), Sweden. The work of E. Molina-Montes was supported by a grant from WCR (15-0391). The authors thank the coordinators, field and administrative workers, technicians, and study participants of the European Study into Digestive Illnesses and Genetics (PanGenEU).

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.

1.
Allemani
C
,
Matsuda
T
,
Di Carlo
V
,
Harewood
R
,
Matz
M
,
Niksic
M
, et al
Global surveillance of trends in cancer survival 2000-14 (CONCORD-3): analysis of individual records for 37 513 025 patients diagnosed with one of 18 cancers from 322 population-based registries in 71 countries
.
Lancet
2018
;
391
:
1023
75
.
2.
Rahib
L
,
Smith
BD
,
Aizenberg
R
,
Rosenzweig
AB
,
Fleshman
JM
,
Matrisian
LM
. 
Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the united states
.
Cancer Res
2014
;
74
:
2913
21
.
3.
Maisonneuve
P
,
Lowenfels
AB
. 
Risk factors for pancreatic cancer: a summary review of meta-analytical studies
.
Int J Epidemiol
2015
;
44
:
186
98
.
4.
Iodice
S
,
Gandini
S
,
Maisonneuve
P
,
Lowenfels
AB
. 
Tobacco and the risk of pancreatic cancer: a review and meta-analysis
.
Langenbeck's Arch Surg
2008
;
393
:
535
45
.
5.
Zou
L
,
Zhong
R
,
Shen
N
,
Chen
W
,
Zhu
B
,
Ke
J
, et al
Non-linear dose-response relationship between cigarette smoking and pancreatic cancer risk: evidence from a meta-analysis of 42 observational studies
.
Eur J Cancer
2014
;
50
:
193
203
.
6.
Bosetti
C
,
Lucenteforte
E
,
Silverman
DT
,
Petersen
G
,
Bracci
PM
,
Ji
BT
, et al
Cigarette smoking and pancreatic cancer: an analysis from the international pancreatic cancer case-control consortium (Panc4)
.
Ann Oncol Off J Eur Soc Med Oncol
2012
;
23
:
1880
8
.
7.
Lugo
A
,
Peveri
G
,
Bosetti
C
,
Bagnardi
V
,
Crippa
A
,
Orsini
N
, et al
Strong excess risk of pancreatic cancer for low frequency and duration of cigarette smoking: a comprehensive review and meta-analysis
.
Eur J Cancer
2018
;
104
:
117
26
.
8.
Samanic
C
,
Kogevinas
M
,
Dosemeci
M
,
Malats
N
,
Real
FX
,
Garcia-Closas
M
, et al
Smoking and bladder cancer in Spain: effects of tobacco type, timing, environmental tobacco smoke, and gender
.
Cancer Epidemiol Biomarkers Prev
2006
;
15
:
1348
54
.
9.
Tajuddin
SM
,
Amaral
AFS
,
Fernandez
AF
,
Rodriguez-Rodero
S
,
Rodriguez
RM
,
Moore
LE
, et al
Genetic and non-genetic predictors of LINE-1 methylation in leukocyte DNA
.
Environ Health Perspect
2013
;
121
:
650
6
.
10.
Peppone
LJ
,
Hyland
A
,
Moysich
KB
,
Reid
ME
,
Piazza
KM
,
Purnell
JQ
, et al
Examining the association between cigarette smoking and colorectal cancer using historical case-control data
.
Cancer Epidemiol
2009
;
33
:
182
8
.
11.
Launoy
G
,
Milan
C
,
Faivre
J
,
Pienkowski
P
,
Gignoux
M
. 
Tobacco type and risk of squamous cell cancer of the oesophagus in males: a French multicentre case-control study
.
Int J Epidemiol
2000
;
29
:
36
42
.
12.
Sancho-Garnier
H
,
Theobald
S
. 
Black (air-cured) and blond (flue-cured) tobacco and cancer risk II: pharynx and larynx cancer
.
Eur J Cancer
1993
;
29A
:
273
6
.
13.
De Stefani
E
,
Boffetta
P
,
Oreggia
F
,
Mendilaharsu
M
,
Deneo-Pellegrini
H
. 
Smoking patterns and cancer of the oral cavity and pharynx: a case-control study in Uruguay
.
Oral Oncol
1998
;
34
:
340
6
.
14.
Molina-Montes
E
,
Gomez-Rubio
P
,
Marquez
M
,
Rava
M
,
Lohr
M
,
Michalski
CW
, et al
Risk of pancreatic cancer associated with family history of cancer and other medical conditions by accounting for smoking among relatives
.
Int J Epidemiol
2018
;
47
:
473
83
.
15.
Gomez-Rubio
P
,
Zock
J-P
,
Rava
M
,
Marquez
M
,
Sharp
L
,
Hidalgo
M
, et al
Reduced risk of pancreatic cancer associated with asthma and nasal allergies
.
Gut
2017
;
66
:
314
22
.
16.
Hosmer
DW
,
Hosmer
T
,
Cessie
SLE
,
Lemeshow
S
. 
A comparison of goodness of-fit tests for the logistic regression model
.
Stat Med
1998
;
16
:
965
80
.
17.
Rothman
K
,
Greenland
S
. 
Modern epidemiology. 2nd ed
.
Philadelphia: Lippincott-Raven
;
1998
.
18.
Hosmer
DW
,
Lemeshow
S
. 
Confidence interval estimation of interaction
.
Epidemiology
1992
;
3
:
452
6
.
19.
Sauerbrei
W
,
Meier-Hirmer
C
,
Benner
A
,
Royston
P
. 
Multivariable regression model building by using fractional polynomials: description of SAS
,
STATA and R programs. Comput Stat Data Anal
2006
;15:3464–85.
20.
Durrleman
S
,
Simon
R
. 
Flexible regression models with cubic splines
.
Stat Med
1989
;
8
:
551
61
.
21.
Harrell
FE
.
Regression modeling strategies: with applications to linear models, logistic regression, and survival analysis
.
New York
:
Springer
; 
2001
.
22.
Strasak
AM
,
Umlauf
N
,
Pfeiffer
RM
,
Lang
S
. 
Comparing penalized splines and fractional polynomials for flexible modelling of the effects of continuous predictor variables
.
Comput Stat Data Anal
2011
;
55
:
1540
51
.
23.
Core Team
.
R: a language and environment for statistical computing
.
Vienna
:
R Foundation for Statistical Computing
; 
2013
[cited 2018 Jan 30]. Available from
: http://www.r-project.org/.
24.
Lynch
SM
,
Vrieling
A
,
Lubin
JH
,
Kraft
P
,
Mendelsohn
JB
,
Hartge
P
, et al
Cigarette smoking and pancreatic cancer: a pooled analysis from the pancreatic cancer cohort consortium
.
Am J Epidemiol
2009
;
170
:
403
13
.
25.
Wienecke
A
,
Barnes
B
,
Lampert
T
,
Kraywinkel
K
. 
Changes in cancer incidence attributable to tobacco smoking in Germany, 1999–2008
.
Int J Cancer
2014
;
134
:
682
91
.
26.
Rosato
V
,
Polesel
J
,
Bosetti
C
,
Serraino
D
,
Negri
E
,
La Vecchia
C
. 
Population attributable risk for pancreatic cancer in Northern Italy
.
Pancreas
2015
;
44
:
216
20
.
27.
Chuang
S-C
,
Gallo
V
,
Michaud
D
,
Overvad
K
,
Tjonneland
A
,
Clavel-Chapelon
F
, et al
Exposure to environmental tobacco smoke in childhood and incidence of cancer in adulthood in never smokers in the European prospective investigation into cancer and nutrition
.
Cancer Causes Control
2011
;
22
:
487
94
.
28.
Vrieling
A
,
Bueno-de-Mesquita
HB
,
Boshuizen
HC
,
Michaud
DS
,
Severinsen
MT
,
Overvad
K
, et al
Cigarette smoking, environmental tobacco smoke exposure and pancreatic cancer risk in the european prospective investigation into cancer and nutrition
.
Int J Cancer
2010
;
126
:
2394
403
.
29.
Zhou
J
,
Wellenius
GA
,
Michaud
DS
. 
Environmental tobacco smoke and the risk of pancreatic cancer among non-smokers: a meta-analysis
.
Occup Environ Med
2012
;
69
:
853
7
.
30.
Bertuccio
P
,
La Vecchia
C
,
Silverman
DT
,
Petersen
GM
,
Bracci
PM
,
Negri
E
, et al
Cigar and pipe smoking, smokeless tobacco use and pancreatic cancer: an analysis from the international pancreatic cancer case-control consortium (PanC4)
.
Ann Oncol
2011
;
22
:
1420
6
.
31.
Ohkuma
T
,
Iwase
M
,
Fujii
H
,
Kaizu
S
,
Ide
H
. 
Dose- and time-dependent association of smoking and its cessation with glycemic control and insulin resistance in male patients with type 2 diabetes mellitus: the fukuoka diabetes registry
.
PLoS Med
2015
;
10
:
e0122023
.
32.
Hu
Y
,
Zong
G
,
Liu
G
,
Wang
M
,
Rosner
B
,
Pan
A
, et al
Smoking cessation, weight change, type 2 diabetes, and mortality
.
N Engl J Med
2018
;
379
:
623
32
.
33.
Harald
B
,
Willi
S
,
Patrick
R
. 
Comparison between splines and fractional polynomials for multivariable model building with continuous covariates: a simulation study with continuous response
.
Stat Med
2012
;
32
:
2262
77
.
34.
Vineis
P
,
Kogevinas
M
,
Simonato
L
,
Brennan
P
,
Boffetta
P
. 
Levelling-off of the risk of lung and bladder cancer in heavy smokers: an analysis based on multicentric case-control studies and a metabolic interpretation
.
Mutat Res
2000
;
463
:
103
10
.
35.
Ruano-Ravina
A
,
Figueiras
A
,
Montes-Martinez
A
,
Barros-Dios
JM
. 
Dose-response relationship between tobacco and lung cancer: new findings
.
Eur J Cancer Prev
2003
;
12
:
257
63
.
36.
Wei
Q
,
Cheng
L
,
Amos
CI
,
Wang
LE
,
Guo
Z
,
Hong
WK
, et al
Repair of tobacco carcinogen-induced DNA adducts and lung cancer risk: a molecular epidemiologic study
.
J Natl Cancer Inst
2000
;
92
:
1764
72
.
37.
Vineis
P
,
Esteve
J
,
Hartge
P
,
Hoover
R
,
Silverman
DT
,
Terracini
B
. 
Effects of timing and type of tobacco in cigarette-induced bladder cancer
.
Cancer Res
1988
;
48
:
3849
52
.
38.
Clavel
J
,
Cordier
S
,
Boccon-Gibod
L
,
Hemon
D
. 
Tobacco and bladder cancer in males: increased risk for inhalers and smokers of black tobacco
.
Int J Cancer
1989
;
44
:
605
10
.
39.
Patrianakos
C
,
Hoffmann
D
. 
Chemical studies of tobacco smoke. LXIV. on the analysis of aromatic amines in cigarette smoke
.
Chem
1979
;
3
:
150
4
.
40.
Ambatipudi
A
,
Cuenin
C
,
Hernandez-Vargas
H
,
Ghantous
A
,
Le Calvez-Kelm
F
,
Kaaks
R
, et al
Tobacco smoking-associated genome-wide DNA methylation changes in the EPIC study
.
Epigenomics
2016
;
8
:
599
618
.
41.
Slikker
WJ
,
Andersen
ME
,
Bogdanffy
MS
,
Bus
JS
,
Cohen
SD
,
Conolly
RB
, et al
Dose-dependent transitions in mechanisms of toxicity
.
Toxicol Appl Pharmacol
2004
;
201
:
203
25
.
42.
Rubio Armendariz
C
,
Garcia
T
,
Soler
A
,
Gutierrez Fernandez
AJ
,
Glez-Weller
D
,
Luis Gonzalez
G
, et al
Heavy metals in cigarettes for sale in Spain
.
Environ Res
2015
;
143
:
162
9
.