Meat-related mutagen exposure, xenobiotic metabolizing gene polymorphisms and the risk of advanced colorectal adenoma and cancer

Meat mutagens, including heterocyclic amines (HCAs), polycyclic aromatic hydrocarbons (PAHs) and N -nitroso compounds (NOCs), may be involved in colorectal carcinogenesis depending on their activation or detoxification by phase I and II xenobiotic metabolizing enzymes (XME). Using unconditional logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI), we examined the intake of five meat mutagens and >300 single nucleotide polymorphisms (SNPs) in 18 XME genes in relation to advanced colorectal adenoma (1205 cases and 1387 controls) and colorectal cancer (370 cases and 401 controls) within the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Dietary intake of meat mutagens was assessed using a food frequency questionnaire with a detailed meat-cooking module. An interaction was observed between 2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline (MeIQx) intake and the NAT1 polymorphism rs6586714 in the ade- noma study ( P interaction = 0.001). Among individuals carrying a GG genotype, high MeIQx intake was associated with a 43% increased risk of adenoma (95% CI 1.11–1.85, P trend = 0.07), whereas the reverse was observed among carriers of the A variant (OR = 0.50, 95% CI 0.30–0.84, P trend = 0.01). In addition, we observed some sug- gestive ( P < 0.05) modifying effects for SNPs in other XME genes ( UGT1A , CYP2E1 , EPHX1 , AHR and GSTM3 ), but these were not significant after adjustment for multiple testing. This large and comprehensive study of XME genes, meat mutagens and the risk of colorectal tumours found that a NAT1 polymorphism modified the association between MeIQx intake and colorectal adenoma risk. NOC, N -nitroso compound; OR, odds ratio; PAH, polycyclic aromatic hydrocar- bon; PhIP, 2-amino-1-methyl-6-phenyl-imidazo[4,5-b]pyridine; PLCO, Prostate, Lung, Colorectal, and Ovarian; SNP, single nucleotide polymor- phisms; UGT, UDP-glucuronosyltransferases; XME, xenobiotic metabolizing enzymes.


Introduction
Experimental evidence suggests that the carcinogenic potential of several meat-specific mutagens may be one of the underlying causal factors for the well-established epidemiological association between red and processed meat consumption and colorectal cancer risk (1). Heterocyclic amines (HCAs) are formed when meat is cooked welldone at high temperature (2,3). Grilled meat may also contain polycyclic aromatic hydrocarbons (PAHs) from smoke that coats the meat. PAHs are also present in cigarette smoke, in the environment and in other food; although concentrations in foods not prepared by grilling or smoking are minor by comparison (2,3). N-nitroso compounds (NOCs) are formed from nitrate and nitrite that are added as preservatives to processed meat (4). Despite the strong carcinogenic potential of HCAs, PAHs and NOCs observed in animal studies (5)(6)(7)(8)(9), evidence in humans remains inconsistent (1).
HCAs, PAHs and NOCs undergo a series of chemical reactions in the human body during which they can be activated or detoxified by phase I and phase II xenobiotic metabolizing enzymes (XMEs) (10,11). Single nucleotide polymorphisms (SNPs) in the genes encoding these XMEs may modify the ability to activate or detoxify carcinogens. Previous studies examining interactions between XME polymorphisms, meat consumption and the risk of colorectal adenomas or carcinomas have reported mixed results (12)(13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)(26). In addition, most previous studies had a small number of cases or examined only a small set of SNPs from a limited number of candidate genes. Because the balance of activating and detoxifying enzymes is thought to influence carcinogen metabolism (27), comprehensive studies including numerous markers across multiple genes involved in xenobiotic metabolism are essential for studying this complex association. Furthermore, the inconsistencies in the data may partly result from the inability of most studies to estimate specific HCAs, PAHs and nitrite/nitrate due to lack of information on both cooking technique and doneness level. Although many studies ask about the consumption of well-done meat, one needs information on both cooking methods and doneness levels to get accurate intake estimates. Studies in different populations need specific databases relevant for their consumption patterns (e.g. Japan, Sweden and USA). The Computerized Heterocyclic Amines Resource for Research in Epidemiology of Disease software application was specifically developed for the US population, therefore, appropriate for our study (3).
We conducted two nested case-control studies within the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial to examine the effect of meat-related mutagens and their interaction with 513 tag SNPs from 18 selected genes involved in phase I [cytochrome P450s (CYP)] and phase II metabolism [sulfotransferases (SULT), N-acetyltransferases (NAT), UDP-glucuronosyltransferases (UGT) and glutathione S-transferases (GST)] in relation to both colorectal adenoma and colorectal cancer risk.

Study population
The PLCO Cancer Screening Trial is a randomized, multicenter clinical trial investigating the efficacy of screening for prostate, lung, colorectal and ovarian cancer (28,29). Participants aged 55-74 years were recruited from 10 centres in the USA from 1993 to 2001 and randomly assigned to the screened or the control arm of the trial. The study was approved by the Institutional Review Board at the National Cancer Institute and the 10 study centres, and all participants provided written informed consent. The present investigation is restricted to the individuals in the screening arm of the trial (n = 77 469), who completed a self-administered risk factor questionnaire at baseline, filled out a food frequency questionnaire (FFQ) with detailed meat-cooking questions, and provided a blood sample.
We conducted two separate nested case-control studies with two disease outcomes: (i) advanced colorectal adenoma and (ii) colorectal cancer. To be eligible for these analyses, participants had to have no history of cancer prior to completion of the FFQ, and no prior history of self-reported colon disease (Crohn's disease, ulcerative colitis, familial polyposis or Gardner's syndrome). Individuals with a prior history of colorectal polyps were also excluded from the analyses of adenoma.
Advanced adenoma sample. Adenoma cases were participants found to have at least one advanced colorectal adenoma (≥1 cm in size, containing villous/ tubulovillous characteristics, or had severe dysplasia) of the distal colon or rectum at baseline. Controls were subjects who underwent a successful sigmoidoscopy examination (defined as insertion to at least 50 cm with ≥90% of mucosa visible or suspicious lesion found) at baseline and had no evidence of a left-sided polyp. Controls were matched to cases on gender, ethnicity and for a subset, age. A total of 1243 advanced adenoma cases and 1419 controls with DNA were available for this study.
Colorectal cancer sample. Colorectal cancer cases were identified through self-report from the annual study update questionnaire, death certificates or physician reports, and confirmed by review of pathology reports and medical records. Cases were identified through 31 December 2006. Controls were subjects without a diagnosis of colorectal cancer at the time the case was diagnosed, matched on age, gender, ethnicity and year of randomization. A total of 371 cases and 405 controls with DNA were identified and genotyped in this study. Five cases from the adenoma sample later on developed colorectal cancer and were also included in the cancer sample.

Dietary assessment
Participants completed a 137-item FFQ with a detailed meat-cooking module that ascertained usual diet during the previous 12 months. Most (89%) study participants completed the FFQ prior to or on the same day as the baseline sigmoidoscopy. Using the Computerized Heterocyclic Amines Resource for Research in Epidemiology of Disease (www.charred.cancer.gov) software application (3), we generated intake estimates of three HCAs (ng/day): [4,5-f]quinoxaline (DiMeIQx) and 2-amino-1-methyl-6phenyl-imidazo [4,5-b]pyridine (PhIP), as well as benzo[a]pyrene (B[a]P), which is a marker of PAH exposure. We estimated nitrate and nitrite from processed meats using a nitrate/nitrite database based on laboratory measured values from 10 types of processed meat samples (bacon; sausage; hot dogs; roast beef; pork chops; ham; lunch meats, including bologna, salami and processed ham) that represented 90% of the processed meat consumed in the USA (3,17). We multiplied the frequency of consumption of each processed meat item by the portion size and by the nitrate or nitrite concentration of the respective item to estimate nitrate and nitrite intake (mg/day). Participants were excluded from these analyses if they had ≥8 missing or invalid responses of the FFQ or missing data on the meat mutagens (n = 71 and n = 18 for the adenoma and cancer samples, respectively), leaving 1386 advanced adenoma cases [171 rectal, 462 distal colon, 566 multiple adenoma (either side) and 6 with unknown location] and 1205 controls and 364 cancer cases (182 proximal colon, 102 distal colon, 79 rectal and 1 unknown location), and 394 controls for analyses.
Whole blood or buffy coat DNA was extracted with QIAamp DNA Blood Midi or Maxi Kits. For quality control purposes, replicate samples from 195 individuals (~7% of the population) were interspersed randomly within the plates. Genotyping was conducted at the National Cancer Institute Core Genotyping Facility, National Institutes of Health. We required all cases and controls to have a 90% call rate for analysis. The overall concordance rate was >99% for replicated samples. We excluded SNPs with a call rate <90%, minor allele frequency <1% in the adenoma set and 5% in the cancer set (due to lower power), or Hardy-Weinberg equilibrium P-value <1 × 10 −6 among Caucasian controls. Of the 513 SNPs we selected, 380 SNPs remained for the adenoma analyses and 325 SNPs remained for the cancer analyses. There were a few highly correlated SNPs among the Caucasians in our study as the tag SNPs were selected for multiple populations. SNPs with an r 2 > 0.95 among Caucasians were considered to be equivalent for interaction testing.

Statistical analysis
All analyses were conducted separately for the prevalent colorectal adenoma and colorectal cancer study populations. Differences in baseline characteristics between cases and controls were assessed using chi-square tests for categorical variables, such as sex, and t-tests for continuous variables, such as age. Due to a skewed intake distribution of the meat mutagens, we transformed these variables prior to performing t-tests, using the inverse rank transformation.

Colorectal Cancer
Step UGT1A UGT1A * Fig. 1. Schematic representation of the two-step test for G×E interaction of Murcray et al. (31). We first excluded SNPs with call rate <90%, minor allele frequency (MAF) <1% in the adenoma set and 5% in the cancer set, or Hardy-Weinberg Equilibrium (HWE) P-value <1 × 10 −6 among Caucasian controls. In Step 1 of the test G×E interaction, G-E association were first tested in the combined case and control sample. Only SNPs with P < 0.10 in Step 1 are then tested for G×E interaction using a standard case-control analysis in Step 2. We further explored the risk pattern for all interactions that met the unadjusted significance threshold of 0.05 (see Tables 3 and 4). *Significant after adjustment for multiple testing at P < 0.05.

Main effects: meat mutagens and SNPs
In order to confirm what has been shown previously in PLCO subgroup analyses, we examined the main effects of each mutagen and SNP on colorectal adenoma and cancer risk using unconditional multivariable logistic regression. The main effects of the meat mutagens were adjusted for all relevant covariates (listed as footnotes in the tables). Meat mutagens were categorized as low (0-39th percentile), medium (40-79th percentile) and high (≥80th percentile) intake, as the distribution of these mutagens is highly skewed with the majority of individuals consuming small amounts, and previous analyses, both within the PLCO study population and within other studies, have identified the top quintile as potentially the most important with regard to cancer risk (33,34) (e.g. PhIP intake in adenoma dataset: quintile 1: range, 0-19.5 mg/day, median, 10.1 mg/ day; quintile 5: range, 166-3069 mg/day, median, 268 mg/day). The main effects of the SNPs were examined using PLINK, a whole-genome association analysis toolset (35) assuming a log-additive model for the genotype, and were adjusted for age, sex and ethnicity. Results were adjusted for multiple testing using the false discovery rate (36). We also used a Bonferroni correction for the total number of tag SNPs for each individual gene or gene region (gene-based correction).

Gene × environment (G×E) interactions
Because the experimental literature on the role of XME genes in relation to HCA, PAH and NOC metabolism is far from complete, we did not restrict our analyses to those gene-mutagen interactions shown previously, but rather used an exploratory approach when testing for interactions between all the SNPs and mutagens under study. We used a two-step test for G×E interactions as described by Murcray et al. (37), to identify SNPs involved in a G×E interaction. In the first step, we examined the association of each SNP with each mutagen in a linear regression model among cases and controls combined, adjusting for age, sex and ethnicity. Due to a skewed intake distribution, we transformed the exposure variables using the inverse normal rank transformation. Linear regressions were calculated using PLINK. The subset of m SNPs that exceeded the significance threshold of P < 0.10 was taken forward to Step 2 for the actual case-control test for a G×E interaction. In Step 2, an unconditional logistic regression model, assuming a dominant genetic model, adjusting for age, sex and ethnicity, was calculated with a multiplicative interaction term of the SNP with the respective meat mutagen (categorical variable assigning median values), as well as the individual main effect terms. Further adjustment for all confounders selected in the mutagen main effects analyses did not change the results appreciably. Significance at this step was defined as a P-value < α/m, where α was the desired overall Type I error rate of 10%. When multiple SNPs from the same gene exceeded the P < 0.10 threshold level, a M eff was calculated to determine the effective number of independent comparisons (38).
All interactions meeting the significance threshold of P < 0.05 were further evaluated using logistic regression models to examine the association between the categorical meat mutagen intake and colorectal adenoma or cancer, stratified by genotype assuming a dominant model. We also conducted stratified analyses by ever and never smokers, rectal and non-rectal adenomas, and single versus multiple adenomas. The top interactions in the cancer set were also examined by proximal, distal and rectal cancer. Models for Step 2 and the models for the association of the meat mutagens and colorectal adenoma or cancer were calculated using STATA version 9.0.

Results
A total of 1205 cases and 1386 controls were available for the advanced colorectal adenoma analysis, and 364 cases and 394 controls were included in the colorectal cancer population (Table I). Over 90% of the study subjects were Caucasian and over 60% were males; in the adenoma study, cases were older than controls.

Main effects: meat mutagens
Consistent with a previous analysis of prevalent adenoma in the PLCO trial (34), there was an elevated risk of colorectal adenoma in the top quintile of MeIQx intake when compared with the bottom quintile (OR = 1.17, 95% CI 0.94-1.45); however, in this study population with ~2000 fewer adenoma cases, the risk did not reach statistical significance and was attenuated after adjustment for confounders (Supplementary Table I is available at Carcinogenesis Online). We observed a borderline significant positive association between nitrate and nitrite intake from meat and adenoma risk (P trend = 0.05), but this was attenuated with covariate adjustment. B[a]P, PhIP and DiMeIQx intake were not associated with adenoma risk in this population. Furthermore, none of the meat mutagens were associated with colorectal cancer. Due to subjects with missing values, the number of cases and controls in the multivariable models were lower (adenoma: 1183 cases and 1358 controls; cancer: 387 cases and 358 controls) than in the age-, ethnicity-and sex-adjusted model; however, this did not explain the observed attenuation. The ORs for MeIQx and colon cancer in particular changed considerably after adjustment for confounders; level of education and total daily energy intake contributed the most to this change.

Main effects: SNPs
We found that SNPs from the UGT1A, EPHX1 and NAT1 genes were nominally associated with both colorectal adenoma (n = 11 SNPs) and cancer risk (n = 20 SNPs) (P < 0.05) (Supplementary Table II is available at Carcinogenesis Online). In addition, one SNP from the GSTA1 locus was associated with adenoma risk, whereas one SNP at the NQO1 gene was associated with colorectal cancer. After adjusting for multiple testing for all the SNPs tested using the false discovery rate, none of these findings remained statistically significant in either the adenoma or cancer analysis. However, UGT1A rs7569014 remained associated with colorectal cancer risk at P ≤ 0.05 after the gene-based multiple testing correction.

Gene × environment interactions
A schematic representation of the results from the two-step test for G×E interaction for both the adenoma and cancer study populations is presented in Figure 1. The subset and number of SNPs that exceeded the significant threshold of P < 0.10 in Step 1, differed by outcome and mutagen under study. SNPs that passed this initial screening step were taken forward to the actual case-control test of G×E interaction analyses in Step 2.
Advanced prevalent adenoma. The G×E interaction analyses performed in Step 2 yielded 15 interactions with P ≤ 0.10 between B[a]P, HCAs and several SNPs, including those in the GSTM3, UGT1A, AHR, EPHX, NAT1 and CYP2C9 gene regions (Table II and Figure 1). The most significant interaction was between NAT1 (rs6586714) and MeIQx (P interaction = 0.001), which was the only finding that remained significant after adjustment for multiple testing. Various SNPs from the UGT1A locus appeared to modify the association between B[a]P (3 SNPs, r > 0.79), PhIP (4 SNPs, r −0.12 to 0.89), and MeIQx intake (1 SNP) and adenoma risk (P < 0.10). Two SNPs (r = 0.39) located within the EPHX1 region (rs2671272 and rs868966) appeared to modify the effect of PhIP and DiMeIQx on colorectal adenoma risk (P ≤ 0.05).
We further explored the risk pattern for all interactions that met the unadjusted significance threshold of 0.05. Age-, ethnicity-and sexadjusted associations between mutagen intake and colorectal adenoma stratified by genotype are presented in Table III. There was an increased risk of colorectal adenoma with increasing intake of MeIQx among individuals carrying a GG genotype for NAT1 rs6586714 (OR = 1.43, 95% CI 1.11-1.85 for high intake compared with low Meat mutagens classified as low (0-39th percentile), medium (40th-79th percentile) and high (≥80th percentile) intake. a P-value is for a test for interaction using the likelihood ratio test comparing models with and without the cross-product of SNP for the given gene and the median of the categories of meat mutagen intake. intake; P trend = 0.007), whereas the risk was decreased among GA and AA carriers (OR = 0.50 for high intake compared with low intake, 95% CI 0.30-0.84; P trend = 0.01; Table III). Increased intake of DiMeIQx was positively associated with adenoma risk in individuals carrying an AA genotype for EPHX1 rs868966 (OR= 1.45, 95% CI 0.98, 2.15; P trend = 0.06).
We additionally ran the G×E interaction analyses for total mutagenic activity, a measure incorporating mutagenicity of all meat-related mutagens. We found five interactions in the adenoma set, whereas one interaction (rs28969701-UGT1A locus) was observed in the cancer data. However, these SNPs came from genes similar to those interacting with the individual mutagens (CYP2C9, EPHX1, GSTM3, UGT1A locus) (data not shown). The correlation between total mutagenic activity and the individual HCAs and PAHs for adenomas was highest for MeIQx (r = 0.84), followed by PhIP (r = 0.78), DiMeIQx (r = 0.76) and B[a]P (r = 0.51).
Subgroup analyses indicated that the top interactions were similar for the subgroups of any rectal (n = 377) and any distal adenoma (n = 957), and single (n = 639) versus multiple adenomas (n = 566). However, the interaction between PhIP intake and rs2018985 was only observed with respect to multiple adenomas but not for single adenomas (data not shown). Moreover, our top interaction were similar in direction and in magnitude for ever (n = 743) and never (n = 462) smokers, except for the interaction between DiMeIQx and rs868966 that was only observed among ever smokers (data not shown).
Colorectal cancer. In the carcinoma dataset, seven interactions with P ≤ 0.10 were observed (Table II and Figure 1). We found a suggestive interaction between several SNPs in the UGT1A region and intake of B[a]P (rs6714486 and rs17868299), MeIQx (rs2011404), DiMeIQx (rs6717546) and combined nitrate and nitrite from meat (rs12466997). The association between MeIQx and colorectal cancer risk seemed to be modified by a SNP in the CYP2E1 gene (rs915908) (P interaction = 0.05). However, neither these nor other interactions for colorectal cancer were deemed significant after adjusting for multiple testing.
Age-, ethnicity-and sex-adjusted associations between mutagen intake and colorectal cancer stratified by genotype are presented in Table III (SNPs with P interaction ≤ 0.05). High consumption of B[a]P was associated with an increased risk of colorectal cancer among individuals carrying the CT genotype for UGT1A rs17868299 (OR = 6.04 for high intake compared with low intake, 95% CI 1.39-24.2; P trend = 0.06). A suggestive inverse association between MeIQx intake and colorectal cancer risk was observed among the wild-type variant of rs915908 (CYP2E1) (OR = 0.68 for high intake compared with low intake, 95% CI 0.40-1.08; P trend = 0.08).
Interactions between SNPs in the UGT1A locus and various meat mutagens were observed for both the adenoma and cancer endpoints. No strong correlation (r < 0.25) between SNPs found in the adenoma and cancer analyses were found (results not shown). However, the interaction observed between rs7571337 and MeIQx intake for adenoma could be replicated for the cancer endpoint (P interaction = 0.09), and the ORs for the SNP-stratified results were in the same direction, though attenuated. The pooled OR (95% CI) combining the adenoma and cancer results for high intake compared with low intake was 1.72 (1.03-2.39) for the TT genotype and 0.95 (0.70-1.13) for the TC+CC genotype. Nevertheless, all other interactions that were significant in the adenoma dataset were not found to be so in the cancer study (P > 0.10), and vice versa (Supplementary Table III is available at Carcinogenesis Online).
The OR and 95% CI for the main effects of the SNPs that were found to be involved in a G×E interaction (P < 0.10) are presented in Supplementary Table IV, available at Carcinogenesis Online . Two SNPs from the UGT1A locus (rs871514 and rs2018985) were marginally associated with colorectal adenoma (P = 0.06), but did not remain statistically significant after adjustment for multiple testing using either the false discovery rate or the gene-based approach.
We examined the top interactions (as described in Table III) by proximal, distal and rectal cancer subsite and observed that the interaction effect for rs6717546 remained present in all endpoints. However, the interaction between B[a]P and rs17868299 was only observed in the distal colon, and the magnitude of the SNP stratified associations for the variant allele was larger than that in the combined group of cancers (CC genotype: ORmedium = 4.36 95% CI 1.02-18.55 and ORhigh = 8.82 95% CI 1.60-48.72) (data not shown); however, since the numbers were small, this may be a chance finding.

Discussion
We observed that the effect of MeIQx intake on advanced colorectal adenoma risk differed by NAT1 genotype. Additional evidence of effect modification of the association between various meat-related HCAs, PAHs and nitrate/nitrite, and advanced colorectal adenoma and cancer by variants in UGT1A, CYP2E1, EPHX1, AHR and GSTM3 was observed; however, correction for multiple comparisons indicated that these observations may be due to chance, despite that this is the largest and most comprehensive study to date.
Our finding that NAT1 modified the effect of MeIQx on the risk of advanced colorectal adenoma is in agreement with results from animal and in vitro research. NAT1 is primarily expressed in hepatic and colonic tissue and has been extensively shown to be an important enzyme in the bioactivation of multiple HCAs via O-acetylation processes (39). NAT1 enzyme activity is characterized by strong between-population heterogeneity that probably results in interindividual differences in acetylation rates (40). Functionality of the variant observed in our analyses (rs6586714), which is located in the intron region of NAT1, is lacking and not strongly correlated with any known functional NAT1 alleles in the iSelect panel (r 2 < 0.20). One previous analysis studied the interaction between various NAT1 alleles (*3, *4, *10, *11; categorizing participants into slow, normal and fast acetylaters) and meat intake in relation to colorectal adenoma risk but reported null results (41). Nevertheless, evidence from the few epidemiological studies that examined interactions between NAT1 and meat (15,23,24) or HCA intake (42) and colorectal cancer risk, points towards an increased risk among rapid acetylaters. However, these findings may have bias towards the null due to the generally small study samples, lack of detailed HCA exposure information and the debatable classification of NAT1 alleles. From our prior research, there is evidence that MeIQx is estimated with greater accuracy than other HCAs (43). In addition, although MeIQx is consumed at lower levels in the diet, it is a more potent mutagen than PhIP is (44). Moreover, it could be that this interaction is observed due to linkage disequilibrium with another variant in the region that has an effect on a neighbouring gene.
Although we observed some suggestive modifying effects in SNPs in other XME genes, the findings were not significant after adjusting for multiple testing. For example, multiple UGT1A polymorphisms were suggestively associated with colorectal adenoma and cancer or appeared to interact with meat mutagens and tumour risk. UGTs are phase II conjugation enzymes that are primarily expressed in the luminal cells of the gastrointestinal tract, and have been implicated in the detoxification of HCAs and B[a]P (45,46). Previous epidemiological studies suggest that alterations in UGT genes may indeed influence the elimination of meat-specific mutagens; one case-control reported an interaction between B[a]P intake and variants in UGT1A1 promoter polymorphisms at position −53 and −3156 (47), whereas a suggestive interaction between DiMeIQx and UGT1A7 (48) acetylation status was reported in a different study. However, the multiple SNPs located across the UGT1A region that were involved in interactions with specific meat mutagens in this study were found in overlapping intron regions of UGT1A1 genes (UGT1A3-A10), and mostly highly correlated, making it impossible to differentiate the associations of each statistically. Future research is needed to unravel this linkage disequilibrium and focus on individual UGT1A isoenzymes to characterize their etiological relevance with meat mutagen-induced colorectal carcinogenesis.

Meat mutagen metabolism and risk of colorectal neoplasia
We found some other suggestive interactions between SNPs in the CYP2E1, AHR and EPHX1 gene, although results after multiple testing suggest that these findings may be due to chance. Although these genes are implicated in the metabolism of meat mutagens, there is little evidence in the literature for the functionality of the observed SNPs. However, rs2066853 (AHR) encodes for an amino acid change in the protein (Arg554Lys), whereas rs15864 (GSTM3) is correlated (r 2 = 0.81) with rs7483, a missense mutation in the same gene responsible for an amino acid substitution (Val224Ile) that significantly alters enzyme activity (49).
Despite the biological rationale for interactions between meatspecific mutagens and polymorphisms in phase I and II XME genes in relation to colorectal carcinogenesis, the epidemiological evidence remains inconclusive. This lack of consistency may be due to inadequate sample size, the high potential for chance findings as a result of multiple comparisons and the over-interpretation of subgroup findings in many gene-environment interaction studies. Although this study is the largest and most comprehensive to examine these associations to date, we were still underpowered to detect modest effects of single gene variants, especially for the colorectal cancers. To gain power, we applied a two-step approach (37), incorporating a preliminary screening step to identify the SNPs involved in possible gene-environment interactions, and hence protect against false positive results. Moreover, the iSelect genotyping panel was specifically designed to capture only tagging SNPs from a priori selected genes involved in carcinogen metabolism to reduce the number of tests performed.
Examination of the same set of SNPs and mutagens yielded different interactions in the adenoma and cancer analyses. Although it may be that some of our findings are false positives, the differences we observed may reflect variation in gene-environmental interactions between initiation and progression in the carcinogenesis process. There is tremendous heterogeneity among colorectal adenomas, with only a small subset of adenomas progressing to cancer. Therefore, interactions related to the formation of adenoma may not be observed for cancer if they do not also stimulate progression to cancer. In contrast, interactions associated with malignant transformation may only be observed for cancer. Nevertheless, we did observe multiple interactions between SNPs at the UGTA1 locus and both colorectal adenoma and cancer, suggesting that this locus may be important for progression and transformation. Many of the XME genes under study are part of larger multigene families. As a result there could be redundancy among these enzymes, such that if the activity of one is modified as a result of a specific genetic variation others may continue to carry out the metabolism of meat mutagens. We were not able to examine differences in etiology between proximal and distal adenomas and advanced and non-advanced adenomas; this study was restricted to advanced distal and rectal colorectal adenomas only.
We made use of a comprehensive assessment of meat-cooking methods that enabled detailed estimation of HCAs, PAHs and nitrate/ nitrite, but this is nonetheless susceptible to measurement error. We did not collect information concerning all aspects of meat-cooking methods that could have influenced carcinogen production, such as microwaving prior to cooking (50), marinating the meat (51) or flipping burgers more often (52). Moreover, during meat cooking, numerous other mutagenic compounds could be formed that may be correlated with the mutagens under study, and could thus potentially explain the observed interactions (e.g. between NAT1 and MeIQx). Our article focuses on meat carcinogenesis and, therefore, we did not consider other potential dietary sources of nitrite and nitrate. Although the majority of dietary nitrite comes from processed meats (53), other food items and drinking water may contribute to intake.
Given that this study was nested within the screening arm of the PLCO Cancer Screening Trial, we examined possible interactions using colorectal adenomas as an endpoint, largely asymptomatic precursors of colorectal cancer. Moreover, selection and surveillance bias is minimal since both cases and controls had an equal opportunity to have a colorectal adenoma or cancer detected. Our colorectal cancer analyses were based on a prospective study of incident cases, eliminating the potential for recall bias. However, inherent to the nature of the study population, participants were generally more educated, less likely to smoke, more physically active and more likely to be Caucasian than the general US population, thereby limiting the generalizability of the study results. All adenoma cases and controls and the majority of cancer cases and controls underwent at least one sigmoidoscopy screening as part of the trial. Since nearly all polyps found were subsequently removed in this heavily screened population, it is possible that the colorectal cancers observed were somehow different than what would be observed in the general population, perhaps more likely to arise de novo. It is also possible that some subjects changed their diet following screening. However, this is unlikely to have affected the results of this study since diet was assessed prior to screening and cancer cases were diagnosed a median of 3.2 years later, which is shorter than the predicted time span (10-20 years) of the adenoma-carcinoma sequence.
In conclusion, we found some evidence that common variants in XME genes may modify the association between meat mutagens and colorectal neoplasia. The strongest evidence for an interaction was observed between MeIQx intake and a NAT1 polymorphism in relation to colorectal adenoma risk. Future pooled initiatives that would have larger sample sizes should further evaluate this complex interplay to better understand these relationships.