Carcinogen metabolism, cigarette smoking, and breast cancer risk: a Bayes model averaging approach
© Stephenson et al. 2010
Received: 26 May 2010
Accepted: 16 November 2010
Published: 16 November 2010
Standard logistic regression with or without stepwise selection has the disadvantage of not incorporating model uncertainty and the dependency of estimates on the underlying model into the final inference. We explore the use of a Bayes Model Averaging approach as an alternative to analyze the influence of genetic variants, environmental effects and their interactions on disease.
Logistic regression with and without stepwise selection and Bayes Model Averaging were applied to a population-based case-control study exploring the association of genetic variants in tobacco smoke-related carcinogen pathways with breast cancer.
Both regression and Bayes Model Averaging highlighted a significant effect of NAT1*10 on breast cancer, while regression analysis also suggested a significant effect for packyears and for the interaction of packyears and NAT2.
Bayes Model Averaging allows incorporation of model uncertainty, helps reduce dimensionality and avoids the problem of multiple comparisons. It can be used to incorporate biological information, such as pathway data, into the analysis. As with all Bayesian analysis methods, careful consideration must be given to prior specification.
Logistic regression and regression with stepwise selection are standard approaches to assess individual and joint effects of genetic and environmental factors on disease risk. However, one drawback is that the resulting estimates depend on the choice of the underlying causal model, and that hence a different set of covariates may lead to different effect estimates and potentially a different pattern of significance. Moreover, standard regression approaches do not incorporate the uncertainty about our choice of the assumed causal model into the final inference.
An alternative approach to analyze such data in combination is Bayes Model Averaging (BMA) , which explicitly accounts for uncertainty with respect to the causal model. BMA specifies prior distributions for model parameters and uses Markov Chain Monte Carlo (MCMC) methods to infer posterior estimates from the priors and from the data. Its inherent model selection feature evaluates different submodels and inference is obtained by averaging over all models considered. By selecting and evaluating a range of submodels, BMA provides a means to reduce dimensionality in the presence of many predictors, when including all variables and their pairwise or higher-order interactions into a logistic model might lead to unstable estimates and bias due to sparse data and correlation . Model selection methods like stepwise regression achieve a similar goal, but do so in a mechanical way, often leading to globally suboptimal and unstable estimates.
Data were derived from a population-based matched case-control study on breast cancer conducted in the study regions "Rhein-Neckar-Odenwald" and "Freiburg" of Southern Germany between 1992 and 1995, as described previously [3, 4]. Cases were diagnosed by age 50 with invasive or in-situ breast cancer, and two controls were matched to cases by age at diagnosis and study region. Participants completed a self-administered questionnaire assessing demographic factors, anthropometric measures and other known or putative risk factors, including smoking history. All study participants gave written informed consent, and the study was reviewed and approved by the ethics committee of the University of Heidelberg, Heidelberg, Germany.
Smoking behavior was assessed over a lifetime, accounting up to eight different phases of active smoking habits. Cumulative cigarette smoking was quantified in packyears, defined as the number of packs of cigarettes smoked per day multiplied by the number of years the individual has smoked.
For the present study we analyzed polymorphisms in the genes NAT1, NAT2, CYP1A1, CYP1B1, as well as the GSTM1 and GSTT1 deletion polymorphisms. Specifically, NAT1 and CYP1B1 genotypes were coded as the number of NAT1*10 and CYP1B1*3 alleles, respectively. NAT2 was coded rapid acetylating conditional on the presence of at least one NAT2*4 allele, which is characterized by the absence of four point mutations (as previously described in ), and slow acetylating otherwise. CYP1A1 was either homozygote for the wild-type allele CYP1A1*1, defined by the absence of three point mutations (rs1056827, rs1056836, rs1800440), or otherwise. GSTM1 and GSTT1 were characterized by the absence of their gene product.
Analyses were adjusted for age, family history of breast cancer in terms of number of affected first-degree relatives, and menopausal status classified as either pre- or postmenopausal, or unknown for women with previous hysterectomy not accompanied by bilateral oophorectomy. Menopausal status was assigned according to the reported state a year before the reference date. Study region showed no effect in an earlier analysis  and was hence not considered in the model.
Study characteristics of the breast cancer case-control study in Germany and variable definition.
42.5 +/- 5.7
42.6 +/- 5.7
first-degree relatives with breast cancer
at least 1
packyears over a lifetime
8.27 +/- 12.20
6.96 +/- 10.88
number *10 alleles
number *3 alleles
presence of at least one *4 allele
homozygote for *1 allele
absence of gene product
absence of gene product
Interactions were only considered if all constituting main effects were in the model. We further restricted the domain of possible gene-gene interactions to polymorphisms in the same pathway. All variables were treated as continuous and were centered on 0.
Logistic regression and backward selection
All logistic regression models contained terms for age, family history, and menopausal status. We tested (i) main effects of smoking and the six polymorphisms separately, (ii) interaction of smoking with each polymorphism, and (iii) gene-gene interactions for polymorphisms in the same pathway. Stepwise regression with backwards selection was used to identify subsets of variables that best explained the data according to the Akaike Information Criterion . We applied stepwise regression to (i) the model containing smoking and all six polymorphisms, and (ii) the model containing all main effects as well as interactions of smoking with all polymorphisms and all gene-gene interactions within pathways.
Bayes Model Averaging (BMA)
The following prior distributions were specified for the model parameters. The probability p ν = Pr(I ν = 1) was beta-distributed with parameters if X ν was a main effect. This prior corresponded to a marginal inclusion probability of 0.25 for main effects and was chosen to reflect our prior emphasis on models with fewer main effects. Moreover, since the inclusion of interactions was limited by the hierarchical dependency on the presence of the main effects, we encouraged inclusion of interaction terms by specifying a greater marginal probability of including term X ν = X rs , provided that I r = I s = 1. Specifically, we set and therefore, the conditional prior probability of including term X rs was E(I rs = 1|I r = I s = 1) = 0.5, and the marginal prior probability was E(I rs = 1) = 0.03 = 0.25 × 0.25 × 0.5.
Thus, if I ν = 1, the variance of β ν was the product of a fixed component σ 2 and of a component depending on whether X ν was a main effect or an interaction. Specifically, we set , such that σ 2corresponded to the variance in main effect coefficients, and modeled the change in that variance for interactions terms. To allow for updating via the Gibbs sampler when no interaction term was present in the current model, a slightly informative prior was chosen for via log(ψ int ) ~ N(2.3,1).
An uninformative prior distribution was placed on σ 2. Since an inverse gamma prior, often chosen for the sake of its conditional conjugacy, leads to improper posterior distributions and sensitivity of inference to hyperparameter choice , we specified a half-Cauchy prior distribution with scale 100, as proposed by Gelman .
Note that via the above prior choices we a priori distinguished main and interaction effects via their model inclusion probabilities and the assumed variance of their coefficients. However, within these two groups variables were treated as exchangeable due to lack of prior evidence suggesting otherwise.
The approach was implemented using the software WinBUGS , running 20,000 iterations and discarding the first 2,000 as burn-in to ensure independence of the results from the initial values. Another 2,000 iterations were discarded through the default settings of WinBUGS to allow the Markov Chain to converge; hence inference was based on 16,000 iterations. We visually inspected plots of the sampled values to ensure convergence of the chain (data not shown). WinBUGS code for our analysis is provided [see Additional file 1].
Selected results from logistic regression and BMA.
Model main a
Model all b
OR (95% CI)
OR (95% CI)
OR (95% CI)
OR (95% CI) &
OR| I = 1 (95% CI) $
BF( I = 1)
1.18 (0.98-1.42) #
1.21 (1.00-1.46) #
NAT2 slow vs fast1
packyears × NAT2
packyears × GSTM1
For BMA, posterior odds ratios and 95%-confidence intervals (CI) were calculated based on the coefficients β ν across all models (marginal odds ratios) and across all models with I ν = 1 (conditional odds ratios). We also report the posterior probability that a variable is included in the model, expecting that variables harboring an association with disease will be included more frequently. Significance of findings was assessed via Bayes factors (BF), the ratio of posterior to prior odds that a variable was included in the model . Thus evaluation of support of a non-zero coefficient took into account the specified prior distribution. Two model Bayes factors were computed in a similar fashion to evaluate support of a selected model versus (i) competing models, and (ii) the null model.
In our prior specifications we emphasized sparse models, resulting in coefficient estimates of 0 for many terms. Hence, expected values of marginal posterior odds ratios showed considerable shrinkage towards 1 with tight confidence intervals due to the hierarchical model (Table 2). The largest effect was observed for NAT1 with OR NAT1 = 1.05 (95%-CI = 1.05-1.05). To assess magnitude of effects conditional on model inclusion we also tabulate expected odds ratios and confidence intervals conditional on model inclusion. The resulting values were similar to the three logistic regression scenarios, but again with much tighter confidence intervals.
The latent indicator variable I was used to compute the posterior probability of model inclusion for each variable. The most frequently selected predictors were NAT1 (Prposterior = 0.26), NAT2 (Prposterior = 0.13) and packyears (Prposterior = 0.13). The posterior probability for the interaction term of packyears and NAT2 was decreased at 0.01, due to the additional restriction that both main effects had to be present in the model.
We used Bayes factors (BFs), the ratio of posterior and prior odds that a variable was selected into the model, to assess the significance of a result in relation to the prior that had been assumed before the analysis. The following calibration has been proposed by Kass and Raftery  to interpret Bayes factors: between 1 and 3 suggests very mild evidence, between 3 and 20 positive evidence, between 20 and 150 strong, and above 150 very strong evidence for an association. Based on these guidelines, very mild evidence was found for NAT1 (BF = 1.05), while all other terms exhibited Bayes factors below 1.
Model results for selected models.
In the comparison of one model to all remaining ones, we found positive evidence for the null model (BFall = 3.8) and for the model containing only NAT1 (BFall = 3.4). Very mild evidence was suggested for the single-effect models of packyears (BFall = 1.4), NAT2 (BFall = 1.3) and GSTT1 (BFall = 1.1). Very mild evidence was also found for the combination of NAT1 with each GSTT1, GSTM1, and CYP1B1 (BFall = 1.4, 1.1, and 1.1, respectively). When taking into account interactions, none of the models exhibited a Bayes factor greater than 1. The same was true for Bayes factors versus the null model.
Hyperparameter scenarios for the sensitivity analysis.
Variation of p ν :
E(p main )
E(p int )
Variation of ψ int :
Variation of σ 2 :
scale( σ 2 )
Estimates of posterior odds ratios showed little variation for different expected prior values of p ν . The posterior probability of model inclusion changed according to the changes in the prior parameters, i.e. doubling the prior probability of including an effect typically led to twice the posterior probability of actually including it.
Varying the mean μ int in the log-normal distribution of ψ int showed no effect on the results, neither did choosing a different scale of the prior for σ 2.
Both logistic regression and BMA highlighted a significant effect of NAT1. Furthermore, logistic regression showed significant effects of packyears and of the interaction of packyears with NAT2 on breast cancer risk. The role of NAT1 as strongest effect is supported by the Bayesian analysis of selected models. Stepwise regression analysis indicated the additional involvement of CYP1B1 and of the interaction of packyears and GSTM1 in breast carcinogenesis.
On a biological level, NAT1 was initially implicated in breast cancer susceptibility through a report of a positive association of the NAT1*11 allele with breast cancer risk as well as combined effects with cigarette smoking and meat consumption , which was, however, not confirmed in a subsequent study . The inconsistent results could be attributed to sample size requirements necessary for assessing effects of NAT1*11, which occurs in approximately only 3% of the general population . We studied the NAT1*10 allele, which occurs with much greater frequency in the Caucasian population than the NAT1*11 allele, and may be rapid acetylating. NAT1*10 has been reported to be associated with higher NAT1 activity in both bladder and colon tissue [13–15]. However, the association between the NAT1*10 allele and increased NAT1 activity in vivo has not been confirmed in other studies [16–18]. For breast cancer, no significant effect of NAT1*10 has been found in several studies [10, 11, 19].
Detection of a gene effect with odds ratio in the order of magnitude that we have found for NAT1 with 80% power at a significance level of 0.05 (assuming allele frequency 0.17, population risk 10%, log-additive disease model and unmatched 1:2 case-control design) requires 1,088 cases and twice the number of controls . Thus the previous studies, as well as our own, would not have enough power to consistently detect such an effect.
Our results from logistic regression analysis regarding the association of NAT2 with breast cancer risk, as previously reported , are in line with findings from other studies. In a meta- and pooled analysis including 13 studies, NAT2 was not independently associated with breast cancer risk but smoking was found to be associated with increased risk in NAT2 slow acetylators but not in rapid acetylators .
The GSTM1 null genotype has not been found to confer susceptibility to breast cancer . However, smokers carrying the GSTM1 null genotype were at significantly elevated risk for breast cancer overall in a meta-analysis of seven studies . An earlier pooled analysis of another seven smaller almost non-overlapping studies, however, did not show clear effect modification in the association between GSTM1 and smoking . Our results from stepwise regression showed a non-significant effect modification by GSTM1, with higher risk of breast cancer associated with smoking among those with the GSTM1 null genotype.
Results from regression and Bayesian analyses differed in that univariate BMA analysis identified only NAT1 as significant and did not yield significant findings for packyears and the interaction of packyears and NAT2. One possible explanation is that inference from BMA is based on posterior and prior probabilities instead of p-values. Thereby it avoids the problem of multiple comparisons inherent in pointwise testing of coefficients in a logistic model. In fact, none of the findings from logistic regression remain significant when Bonferroni-corrected for multiple testing. However, there is no simple one-to-one correspondence between frequentist and Bayesian analyses, since the latter explicitly depend on the specified priors. Our results were stable for different hyperparameter choices. However, adequate prior specification always needs to be kept in mind before starting any Bayesian analysis. In our case, mostly uninformative prior distributions were specified to reflect the lack of sufficient external information justifying an a priori distinction of variables. We therefore allowed the data more weight versus prior information in estimation and model selection. However, if desired, BMA provides a framework for the explicit inclusion of biological prior information, like pathway characteristics, into the analysis through prior specification. If one is confident about biological prior information, stronger prior assumptions may be helpful to guide the analysis. However, bias will be introduced at the same time, so that this trade-off must be carefully considered.
The strength of BMA is its explicit statement of the prior assumptions given by the prior distributions for model parameters, and its consideration of model uncertainty by obtaining results averaged over a multitude of possible models. It evaluates single variables and a range of models at the same time, yielding stabilized estimates based on a set of potential data-generating models. Moreover, it provides a means to reduce dimensionality and avoids the problem of multiple comparisons.
In our study both BMA and regression analyses yielded a significant effect of NAT1*10, while BMA attenuated other significant findings from logistic regression. Since all Bayesian inference depends on the specified prior information, prior choice must be carefully considered when conducting a Bayesian analysis.
NS is funded by the DKFZ International PhD Program and LB by DFG grant BE 3906/2-2. The authors thank two anonymous reviewers for helpful suggestions on the manuscript.
- Conti DV, Cortessis V, Molitor J, Thomas DC: Bayesian modeling of complex metabolic pathways. Hum Hered 2003, 56:83–93.View ArticlePubMedGoogle Scholar
- Greenland S: Bayesian perspectives for epidemiological research. II. Regression analysis. Int J Epidemiol 2007, 36:195–202.View ArticlePubMedGoogle Scholar
- Chang-Claude J, Kropp S, Jager B, Bartsch H, Risch A: Differential effect of NAT2 on the association between active and passive smoke exposure and breast cancer risk. Cancer Epidemiol Biomarkers Prev 2002, 11:698–704.PubMedGoogle Scholar
- Kropp S: The Association between Tobacco Exposure and Breast Cancer Risk considering potentially modifying Effects of specific Genotypes. In PhD Thesis. University of Heidelberg; 2001.Google Scholar
- Kropp S, Chang-Claude J: Active and passive smoking and risk of breast cancer by age 50 years among German women. Am J Epidemiol 2002, 156:616–626.View ArticlePubMedGoogle Scholar
- Hastie TJ, Pregibon D: Generalized linear models. 1992.Google Scholar
- Gelman A: Prior distributions for variance parameters in hierarchical models. Bayesian Analysis 2006, 1:515–533.View ArticleGoogle Scholar
- Spiegelhalter DJ: WinBUGS version 2.10. 2005.Google Scholar
- Raftery AE: Bayes Factors. J Am Statist Assoc 1995, 90:773–795.View ArticleGoogle Scholar
- Zheng W, Deitz AC, Campbell DR, Wen WQ, Cerhan JR, Sellers TA, Folsom AR, Hein DW: N-acetyltransferase 1 genetic polymorphism, cigarette smoking, well-done meat intake, and breast cancer risk. Cancer Epidemiol Biomarkers Prev 1999, 8:233–239.PubMedGoogle Scholar
- Millikan RC: NAT1*10 and NAT1*11 polymorphisms and breast cancer risk. Cancer Epidemiol Biomarkers Prev 2000, 9:217–219.PubMedGoogle Scholar
- Loktionov A, Moore W, Spencer SP, Vorster H, Nell T, O'Neill IK, Bingham SA, Cummings JH: Differences in N-acetylation genotypes between Caucasians and Black South Africans: implications for cancer prevention. Cancer Detect Prev 2002, 26:15–22.View ArticlePubMedGoogle Scholar
- Hughes NC, Janezic SA, McQueen KL, Jewett MA, Castranio T, Bell DA, Grant DM: Identification and characterization of variant alleles of human acetyltransferase NAT1 with defective function using p-aminosalicylate as an in-vivo and in-vitro probe. Pharmacogenetics 1998, 8:55–66.View ArticlePubMedGoogle Scholar
- Payton MA, Sim E: Genotyping human arylamine N-acetyltransferase type 1 (NAT1): the identification of two novel allelic variants. Biochem Pharmacol 1998, 55:361–366.View ArticlePubMedGoogle Scholar
- Bruhn C, Brockmoller J, Cascorbi I, Roots I, Borchert HH: Correlation between genotype and phenotype of the human arylamine N-acetyltransferase type 1 (NAT1). Biochem Pharmacol 1999, 58:1759–1764.View ArticlePubMedGoogle Scholar
- Yang M, Katoh T, Delongchamp R, Ozawa S, Kohshi K, Kawamoto T: Relationship between NAT1 genotype and phenotype in a Japanese population. Pharmacogenetics 2000, 10:225–232.View ArticlePubMedGoogle Scholar
- Badawi AF, Hirvonen A, Bell DA, Lang NP, Kadlubar FF: Role of aromatic amine acetyltransferases, NAT1 and NAT2, in carcinogen-DNA adduct formation in the human urinary bladder. Cancer Res 1995, 55:5230–5237.PubMedGoogle Scholar
- Bell DA, Stephens EA, Castranio T, Umbach DM, Watson M, Deakin M, Elder J, Hendrickse C, Duncan H, Strange RC: Polyadenylation polymorphism in the acetyltransferase 1 gene (NAT1) increases risk of colorectal cancer. Cancer Res 1995, 55:3537–3542.PubMedGoogle Scholar
- van der Hel OL, Bueno-de-Mesquita HB, van Gils CH, Roest M, Slothouber B, Grobbee DE, Peeters PH: Cumulative genetic defects in carcinogen metabolism may increase breast cancer risk (The Netherlands). Cancer Causes Control 2005, 16:675–681.View ArticlePubMedGoogle Scholar
- Gauderman J, Morrison J, QUANTO: A computer program for power and sample size calculations for genetic-epidemiology studies. (1.2.4). 2006.Google Scholar
- Ambrosone CB, Kropp S, Yang J, Yao S, Shields PG, Chang-Claude J: Cigarette smoking, N-acetyltransferase 2 genotypes, and breast cancer risk: pooled analysis and meta-analysis. Cancer Epidemiol Biomarkers Prev 2008, 17:15–26.View ArticlePubMedGoogle Scholar
- Vogl FD, Taioli E, Maugard C, Zheng W, Pinto LF, Ambrosone C, Parl FF, Nedelcheva-Kristensen V, Rebbeck TR, Brennan P, et al.: Glutathione S-transferases M1, T1, and P1 and breast cancer: a pooled analysis. Cancer Epidemiol Biomarkers Prev 2004, 13:1473–1479.PubMedGoogle Scholar
- Terry PD, Goodman M: Is the association between cigarette smoking and breast cancer modified by genotype? A review of epidemiologic studies and meta-analysis. Cancer Epidemiol Biomarkers Prev 2006, 15:602–611.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.