
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Occupational Health Program, Departments of Environmental Health [W. Z., L. L. X., D. P. M., L. S., D. C. C.] and Biostatistics [S. W. T.], Harvard School of Public Health, Boston, Massachusetts, 02115, and Department of Hematology-Oncology [G. L., T. J. L.] and Pulmonary and Critical Care Unit [D. C. C.], Department of Medicine, and Thoracic Surgery Unit, Department of Surgery [J. C. W.], Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, 02114
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
There are seven widely studied genetic polymorphisms associated with decreased enzyme activity and/or variable stability for the human NAT2 gene: 191G
A, 282C
T, 341T
C, 481C
T, 590G
A, 803A
G, and 857G
A (2)
. For mEH, there are two polymorphisms that affect enzyme activity: 113T
C (exon 3) and 139A
G (exon 4), where the variant genotypes confer decreased and increased activity, respectively (4)
. Because tobacco smoking is the major environmental risk factor for lung cancer, these genetic polymorphisms may modify host susceptibility to lung cancer through their roles in metabolizing tobacco procarcinogens. Most genotype or phenotype studies found no overall lung cancer risk related to NAT2 rapid or slow acetylators (5, 6, 7, 8, 9)
, whereas the slow acetylator genotype (10, 11, 12)
and the homozygous rapid acetylator genotype (13)
were associated with increased lung cancer risk in specific studies or populations. One study suggested that the rapid acetylator genotype was a protective factor in never-smokers but a risk factor in smokers, but the results were not conclusive, and the study population was too small for formal interaction analyses (n = 185, P > 0.05; Ref. 14
). In the case of mEH polymorphisms, inconsistent relationships have also been found with respect to lung cancer risk (15, 16, 17, 18, 19, 20)
.
We recently reported a gene-environment interaction between mEH polymorphisms and cumulative smoking exposure in lung cancer development, where the higher activity genotypes (compared with the very low activity genotype) were protective against lung cancer in nonsmokers but were a significant risk factor in heavy smokers (21) . Because NAT2 and mEH have similar putative dual functions and genetically polymorphic variants that affect enzyme function, a gene-smoking interaction may also modify the relationship between NAT2 polymorphisms and lung cancer risk.
The primary aim of this study was to determine whether cumulative smoking exposure altered the relationship between NAT2 polymorphisms and lung cancer risk in a gene-environment interaction analysis. If NAT2 genotypes were found to interact with cigarette smoking history, a secondary aim of this study was to determine whether the combined effects of NAT2 and mEH polymorphisms on lung cancer risk were additive, because their corresponding enzymes metabolized different tobacco-related substrates (1, 2, 3) .
| Materials and Methods |
|---|
|
|
|---|
Cumulative Smoking Exposure.
The questionnaire asked detailed smoking questions, including the age of onset of regular smoking, the average number of cigarettes smoked during different time periods, duration of smoking, and details on periods of smoking cessation. Pack-years were calculated as the product of the number of years an individual smoked and the average number of cigarettes smoked per day (converted into a standard pack of 20 cigarettes). Total periods of smoking cessation for periods >6 months were excluded from the total number of years of smoking in all calculations. Nonsmokers were defined as individuals who smoked fewer than 100 cigarettes in their lifetime.
NAT2 and mEH Genotyping.
Peripheral blood samples were obtained from each case or control, and DNA was extracted using the Puregene DNA Isolation Kit (Gentra Systems, Minneapolis, MN). NAT2 polymorphisms were detected using previously reported PCR-RFLP methods (22
, 23)
. In brief, the four most common NAT2 polymorphisms, 481C
T, 590G
A, 803A
G, and 857G
A, were identified by different restriction fragment analysis on digestion of the 547-bp PCR product, using restriction enzymes KpnI (481), TaqI (590), DdeI (803), and BamHI (857; New England BioLabs, Beverly, MA). DNA-digested PCR products were separated by size on a 2.5% agarose gel (Sigma Chemical Co., St. Louis, MO).
The genotyping methods for mEH polymorphisms have been described in detail (21) . In brief, two separate PCR assays were used to detect the polymorphisms in exons 3 and 4 of mEH, and EcoRV and RsaI enzyme digestion (New England BioLabs) were used respectively for restriction fragment analyses.
For quality control, a random 5% of the samples were repeated to assess the reproducibility of results, and a second investigator rechecked all data entry and analyses.
Statistical Analysis.
NAT2 and mEH genotypes were categorized according to their corresponding genetic polymorphisms. For NAT2, the wild type of each genetic polymorphism corresponds with rapid acetylation, and the variant genotypes correspond to slow acetylation. Thus, individuals carrying wild-type alleles (i.e., two "rapid" alleles, NAT2 *4/*4) were designated as carrying the homozygous rapid acetylator genotype; individuals with only one copy of a rapid allele (heterozygous variant genotype, NAT2 *4) were identified as heterozygous rapid acetylators, and those with two "slow" alleles (homozygous variant genotypes; NAT2 *5A, *5B, *5C, *6, or *7) were labeled slow acetylators. A similar system was used to categorize the mEH genotype into very low, low, intermediate, and high activity groups, according to different mEH genotypes in exons 3 and 4 (21)
.
Gene-environment interaction and gene-gene joint effects analyses involving multiple strata require enormous sample sizes. Combining different genotype strata on the basis of biological or epidemiological evidence has been a reasonable and practical approach. We combined NAT2 homozygous and heterozygous rapid acetylators into one group, comparing it with the slow acetylators, similar to approaches taken by others (9 , 11 , 12 , 14 , 23) . On the basis of our previous epidemiological research on mEH genotypes and lung cancer risk (21) , mEH was dichotomized into the genotype representing very low enzyme activity (homozygous for histidine in both exons and 4) versus all other genotypes combined.
Although individuals of all races were recruited for this study, we restricted our analyses to Caucasians to minimize the possibility of a different gene effect, i.e., effect modification by race or unmeasured confounders that vary by race. We included all Caucasians with complete information on age, gender, smoking status (non, ex-, and current smokers), pack-years of smoking, and years since smoking cessation (for ex-smokers). Logistic regression and GAM (24) were used to examine the relationship between the log odds of lung cancer and each covariate, after adjusting for possible confounding factors, such as age (years), gender, pack-years of smoking, smoking status, and years since smoking cessation (defined for ex-smokers; this variable was defined as zero for nonsmokers and current smokers). Some variables were transformed to approach normality. GAM allows the relationship between lung cancer risk and each covariate to be an unspecified smooth function. GAM plots showed a linear association between lung cancer risk and SR-PY but not for the untransformed variable. Subsequent analyses were based on logistic regression models, using SR-PY, age, and years since smoking cessation as continuous variables. In the NAT2 genotype-smoking primary interaction analysis, we included in our initial model all possible interactions among NAT2 genotype and two smoking variables (smoking status and SR-PY). To examine the combined effects of the NAT2 and mEH genotypes with smoking, we constructed a model that included all of the possible interactions among NAT2 and mEH genotypes and SR-PY, as well as any interactions between smoking status and either NAT2, mEH, or SR-PY that were included in the final models of NAT2 (or mEH) genotype-smoking interaction analysis. When the highest order interaction term was not significant, we considered models with only lower order interactions. Interaction terms were removed from the model when the likelihood ratio tests (comparing a candidate smaller model with all larger models in which it was nested) were not significant. Where appropriate, OR and 95% CI for the risk of lung cancer were calculated from these models. A lack of fit test was performed to summarize the goodness of fit of each logistic regression model (25) . Statistical analyses were all undertaken using the S-plus (MathSoft, Inc, Cambridge, MA) and SAS statistical packages (SAS Institute, Cary, NC).
| Results |
|---|
|
|
|---|
The mean age (and SD) of all cases was 65 ± 10.8 years (range 2691) and of controls, 58.5 ± 12.3 years (range 19100). Females represented 46.5% of cases and 54.2% of controls. Among cases, the mean pack-years of ex-smokers was 55.2 ± 35.9 and of current smokers, 64.3 ± 35.7; among controls, the numbers were 29.3 ± 27.7 and 38.1 ± 25.2, respectively. The years since smoking cessation for ex-smokers were 14.1 ± 11 and 19.3 ± 12 for cases and controls, respectively.
Of the seven widely studied genetic NAT2 polymorphisms, we identified the most frequently occurring four: 481C
T, 590G
A, 803A
G, and 857G
A (11
, 22
, 23)
. We also determined two other polymorphisms (341T
C and 282C
T) in 312 random DNA samples, using standard methods (27)
, and found that 99.7% of the 341T
C polymorphism was linked with the 481C
T or 803A
G polymorphisms, and 100% of the 282C
T polymorphism results were concordant with the 590G
A or 857G
A polymorphisms (2
, 7
, 13
, 14)
. The seventh polymorphism (191G
A), generally restricted to African populations, was not evaluated in our Caucasian-restricted samples (7
, 13)
.
Distribution of NAT2 and mEH Polymorphisms Among Cases and Controls.
All four polymorphic alleles of NAT2 were consistent with the Hardy-Weinberg equilibrium in our control population (P > 0.05,
2 goodness of fit). The frequencies of genotypes representing slow, heterozygous rapid, and homozygous rapid acetylators were 60.2, 35.8, and 4% in cases; these frequencies were very similar in controls (59.9, 35.8, and 4.2%, respectively; P > 0.05). As pack-years increased, the frequencies of slow acetylator genotypes decreased in cases (
2 trend test, P < 0.05) but not in controls (P > 0.05), e.g., the frequencies were 67% in cases and 56% in controls for nonsmokers and 61% in cases and 55% in controls with pack-years between 60 and 79, respectively.
A total of 896 cases and 1131 controls were genotyped for both NAT2 and mEH polymorphisms. The frequencies of the genotype combinations of NAT2 and mEH by pack-years and smoking status are shown in Table 1
. Baseline characteristics of mEH genotypes have been reported previously (21)
.
|
Interaction Between NAT2 Genotype and Cumulative Cigarette Smoking in Lung Cancer Risk.
NAT2 genotypes representing homozygous and heterozygous rapid acetylators were combined into one group for this analysis. In the gene-tobacco smoking interaction analysis, the three- and two-way interactions between smoking status and genotype were not significant, and these were removed from the model. The final model included interactions between NAT2 genotype and SR-PY (P = 0.04) and between smoking status and SR-PY (P = 0.03) and was adjusted for age, gender, NAT2 genotype, SR-PY, smoking status, and years since smoking cessation. In this model, the ORs of the rapid versus slow acetylator genotypes increased as SR-PY increased (Table 2)
. At 0 pack-years (nonsmokers), the fitted OR was 0.66 (95% CI, 0.440.99), indicating that the rapid acetylator genotype was protective against lung cancer. In contrast, at 80 pack-years, the fitted OR increased to 1.22 (95% CI, 0.891.67). No statistically significant results were found when cases of different histological types or disease stages were compared separately to all controls (P > 0.05 for all comparisons).
|
|
| Discussion |
|---|
|
|
|---|
Our study has several important strengths, including a large sample size that permitted formal analysis of interactions between NAT2 genetic polymorphisms and smoking. This large sample size permitted the use of lifelong nonsmokers as a stable baseline comparison group in gene-environment interaction analyses. In addition, complete smoking data were collected systematically, diagnoses of lung cancer were histologically confirmed, and females were well represented in both cases and controls. In an attempt to avoid overmatching, we did not individually match our controls to cases on sex, race, or smoking variables (28 , 29) . Our method resulted in a reasonable balance of these characteristics between cases and controls, because the distributions of smoking variables in our controls were similar to the general Massachusetts population overall and in age- and gender-specific strata (26) .4 ,5
The use of fitted ORs can occasionally lead to the identification of artifactual relationships among variables. Before we carried out the gene-smoking interaction analysis, we used a two-step test to identify possible misclassification of subjects with respect to cigarette smoking exposure (30) . An association between NAT2 genotype and SR-PY was found in the case-only analysis (P = 0.04) but not in the control-only analysis (P > 0.05). Because the association between NAT2 genotype and SR-PY was not statistically significant in controls, the true smoking exposure is independent of the genotype among controls; if there is misclassification of smoking exposure, it is likely that the possible misclassification is nondifferential with respect to the genotype. According to Garcia-Closas et al. (30) , the pattern of association observed in our study may actually underestimate the true gene-smoking interaction but would not lead to an artifactual statistical interaction. In addition to fitted ORs, we obtained ORs from separate analyses that stratified at different pack-years of cigarette smoking (heavy versus non or light smokers and at different quartiles of pack-years). The crude and adjusted ORs for different strata were similar in magnitudes and directions to the fitted ORs derived from the interaction models (data not shown). Additionally, we decomposed pack-years into its component parts of smoking intensity (mean number of cigarettes per day) and duration (in years) and found similar gene-smoking interaction associations when either cigarettes per day or smoking years (both as continuous variables) were substituted for pack-years in the regression models (data not shown).
Because the correlation between NAT2 genotype and phenotype has been reported to be as high as 9398% (31 , 32) , genotyping is generally accepted as an accurate and efficient means to determine acetylator status in large epidemiological studies (33) . In addition, we provide data on the concordance of nonevaluated NAT2 polymorphisms to suggest that the four common polymorphisms truly represent the enzyme activities of NAT2 in our Caucasian samples. Compared with other studies (8 , 12 , 14) , the frequencies of rapid and slow acetylators (as defined by genotype) in our study were very similar.
There are biological reasons to consider an interaction between NAT2 polymorphisms and smoking variables, because the purported mechanism for modifying lung cancer risk by these polymorphisms is dependent on their role in metabolizing tobacco-related substrates. This interaction is particularly important in light of evidence for a dual role of activation and detoxification of tobacco procarcinogens and carcinogens by the NAT2 enzyme. For nonsmokers, NAT2 appears to provide a competitive alternative pathway to the cytochrome P4501A2 and P4501A1-catalyzed N-oxidation pathway for reactive hydroxylamines without an alternative mechanism. These intermediates are then further metabolized into DNA-binding electrophiles (34) . NAT2 N-acetylates and detoxifies these aromatic amine compounds (35 , 36) and is a possible explanation for why the rapid acetylator genotypes showed a possible protective effect in nonsmokers.
By contrast in smokers, the amount of inhaled cigarette smoke appears to differentially induce cytochrome P4501A2 and 1A1-catalyzed N-oxidation, whereby heavy smoking significantly induces enzyme activity and increases the number of these reactive intermediates (37) . NAT2 can then O-acetylate these hydroxylamines to acetoxy esters, which can spontaneously disintegrate into highly reactive aryl nitrenium (38) . These reactive products have been associated with other genetic derangements, such as p53 mutations in cancer tissues (39 , 40) . Thus, higher NAT2 activity may result in an increased lung cancer risk in heavy smokers.
The varying relationship between NAT2 genotypes and lung cancer risk at different pack-years is supported by indirect evidence from several adduct studies. Adducts (4-aminobiphenyl-hemoglobin) were higher in slow acetylators at low or null nicotine-cotinine levels, whereas this difference was less evident at increasing nicotine-cotinine levels, and adduct levels were higher in rapid acetylators at the highest nicotine-cotinine levels (41 , 42) . In a Caucasian study of subjects >60 years old, rapid acetylator smokers had significantly higher frequencies of chromosome aberration than slow acetylator smokers (43) .
Specific evaluations of the smoking habits of the underlying lung cancer study populations in other studies provide some corroborative evidence. The possible protective effect of the rapid acetylator genotypes for lung cancer observed among nonsmokers in this population has also been reported in two previous studies with nonsmoker predominant populations (11
, 14)
. In contrast, the homozygous rapid acetylator genotype was reported to be a risk factor for lung cancer in another study of Caucasians; 77% of the 155 cases in this particular study were smokers with >20 pack-years of smoking (13)
. However, the lack of uniformity in the design, analysis, smoking data, and reporting may explain discrepancies across previous reports of NAT2 polymorphisms and lung cancer risk (Table 3)
.
|
The results of this study provide empiric evidence for two important interrelated paradigms. Firstly, genetic polymorphism-environmental relationships can influence cancer risk. Our results illustrate how negative overall results can hide important relationships that are only uncovered when an environmental factor integral to the proposed pathogenetic mechanism (such as smoking) is incorporated into the analysis. Secondly, combinations of different genetic polymorphisms can have additive effects on cancer risk. Cigarette smoke contains >4000 constituents, including
50 carcinogens (44)
; the carcinogenic effects of these compounds are mediated by multiple biotransformation enzymes in complex interrelated pathways. Our results highlight the potential for individual risk profiling through the evaluation of multiple genetic polymorphisms.
In conclusion, cumulative cigarette smoking exposure appears to play a critical role in altering the direction and magnitude of the association between NAT2 polymorphisms and the risk of lung cancer. Lung cancer risk by combined genotypes of both NAT2 and mEH is dependent on cumulative smoking exposures. The increase in risk conferred by carrying the higher risk mEH genotypes (low, intermediate, and high activity genotypes) and the higher risk rapid acetylator NAT2 genotypes appears to be independent and additive.
| Acknowledgments |
|---|
| Footnotes |
|---|
1 Supported by NIH Grants CA74386, ES/CA 06409, and ES00002. Dr. Zhou was supported by Fogarty Training Grant TW00828; Dr. Liu was supported by Noah Herndon Charles A. King Fellowships; and Dr. Xu was supported by NIH Training Grant T32ES07069. ![]()
2 To whom requests for reprints should be addressed, at Occupational Health Program, Harvard School of Public Health, 665 Huntington Avenue, Boston, MA 02115. Phone: (617) 432-3323; Fax: (617) 432-6981; E-mail: dchris{at}hohp.harvard.edu ![]()
3 The abbreviations used are: NAT2, N-acetyltrasferase-2; mEH, microsomal epoxide hydrolase; GAM, generalized additive model; SR-PY, square root of pack-years; OR, odds ratio; CI, confidence interval. ![]()
4 Internet address: http://www.state.ma.us/dph/mtcp/report/mats.htm. ![]()
5 W. Zhou et al., unpublished data. ![]()
Received 7/ 6/01; revised 10/26/01; accepted 11/ 5/01.
| References |
|---|
|
|
|---|
]pyrene metabolism, activation and carcinogenesis: role and regulation of mixed-function oxidases and related enzymes. Physiol. Rev., 60: 1107-1166, 1980.This article has been cited by other articles:
![]() |
A. G. Schwartz, G. M. Prysak, C. H. Bock, and M. L. Cote The molecular epidemiology of lung cancer Carcinogenesis, March 1, 2007; 28(3): 507 - 518. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. W. Thurston, G. Liu, D. P. Miller, and D. C. Christiani Modeling Lung Cancer Risk in Case-Control Studies Using a New Dose Metric of Smoking Cancer Epidemiol. Biomarkers Prev., October 1, 2005; 14(10): 2296 - 2302. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Suk, S. Gurubhagavatula, S. Park, W. Zhou, L. Su, T. J. Lynch, J. C. Wain, D. Neuberg, G. Liu, and D. C. Christiani Polymorphisms in ERCC1 and Grade 3 or 4 Toxicity in Non-Small Cell Lung Cancer Patients Clin. Cancer Res., February 15, 2005; 11(4): 1534 - 1538. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Zhou, G. Liu, S. Park, Z. Wang, J. C. Wain, T. J. Lynch, L. Su, and D. C. Christiani Gene-Smoking Interaction Associations for the ERCC1 Polymorphisms in the Risk of Lung Cancer Cancer Epidemiol. Biomarkers Prev., February 1, 2005; 14(2): 491 - 496. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Ryan, W. Huang, S. W Thurston, K. T Kelsey, J. K Wiencke, and D. C Christiani On the use of biomarkers for environmental health research Statistical Methods in Medical Research, June 1, 2004; 13(3): 207 - 225. [Abstract] [PDF] |
||||
![]() |
W. Zhou, G. Liu, D. P. Miller, S. W. Thurston, L. L. Xu, J. C. Wain, T. J. Lynch, L. Su, and D. C. Christiani Polymorphisms in the DNA Repair Genes XRCC1 and ERCC2, Smoking, and Lung Cancer Risk Cancer Epidemiol. Biomarkers Prev., April 1, 2003; 12(4): 359 - 365. [Abstract] [Full Text] [PDF] |
||||
![]() |
G E Goodman Lung cancer * 1: Prevention of lung cancer Thorax, November 1, 2002; 57(11): 994 - 999. [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 | Meeting Abstracts Online |