Association Between Genetic Variants on Chromosome 15q25 Locus and Objective Measures of Tobacco Exposure
JNCI-JOURNAL OF THE NATIONAL CANCER INSTITUTE
Authors: Munafo, Marcus R.; Timofeeva, Maria N.; Morris, Richard W.; Prieto-Merino, David; Sattar, Naveed; Brennan, Paul; Johnstone, Elaine C.; Relton, Caroline; Johnson, Paul C. D.; Walther, Donna; Whincup, Peter H.; Casas, Juan P.; Uhl, George R.; Vineis, Paolo; Padmanabhan, Sandosh; Jefferis, Barbara J.; Amuzu, Antoinette; Riboli, Elio; Upton, Mark N.; Aveyard, Paul; Ebrahim, Shah; Hingorani, Aroon D.; Watt, Graham; Palmer, Tom M.; Timpson, Nicholas J.; Smith, George Davey
Abstract
Background Two single-nucleotide polymorphisms, rs1051730 and rs16969968, located within the nicotinic acetylcholine receptor gene cluster on chromosome 15q25 locus, are associated with heaviness of smoking, risk for lung cancer, and other smoking-related health outcomes. Previous studies have typically relied on self-reported smoking behavior, which may not fully capture interindividual variation in tobacco exposure. Methods We investigated the association of rs1051730 and rs16969968 genotype (referred to as rs1051730-rs16969968, because these are in perfect linkage disequilibrium and interchangeable) with both self-reported daily cigarette consumption and biochemically measured plasma or serum cotinine levels among cigarette smokers. Summary estimates and descriptive statistical data for 12 364 subjects were obtained from six independent studies, and 2932 smokers were included in the analyses. Linear regression was used to calculate the per-allele association of rs1051730-rs16969968 genotype with cigarette consumption and cotinine levels in current smokers for each study. Meta-analysis of per-allele associations was conducted using a random effects method. The likely resulting association between genotype and lung cancer risk was assessed using published data on the association between cotinine levels and lung cancer risk. All statistical tests were two-sided. Results Pooled per-allele associations showed that current smokers with one or two copies of the rs1051730-rs16969968 risk allele had increased self-reported cigarette consumption (mean increase in unadjusted number of cigarettes per day per allele = 1.0 cigarette, 95% confidence interval [CI] = 0.57 to 1.43 cigarettes, P = 5.22 x 10(-6)) and cotinine levels (mean increase in unadjusted cotinine levels per allele = 138.72 nmol/L, 95% CI = 97.91 to 179.53 nmol/L, P = 2.71 x 10(-11)). The increase in cotinine levels indicated an increased risk of lung cancer with each additional copy of the rs1051730-rs16969968 risk allele (per-allele odds ratio = 1.31, 95% CI = 1.21 to 1.42). Conclusions Our data show a stronger association of rs1051730-rs16969968 genotype with objective measures of tobacco exposure compared with self-reported cigarette consumption. The association of these variants with lung cancer risk is likely to be mediated largely, if not wholly, via tobacco exposure.
Comparison of Pathway Analysis Approaches Using Lung Cancer GWAS Data Sets
PLOS ONE
Authors: Fehringer, Gordon; Liu, Geoffrey; Briollais, Laurent; Brennan, Paul; Amos, Christopher I.; Spitz, Margaret R.; Bickeboeller, Heike; Wichmann, H. Erich; Risch, Angela; Hung, Rayjean J.
Abstract
Pathway analysis has been proposed as a complement to single SNP analyses in GWAS. This study compared pathway analysis methods using two lung cancer GWAS data sets based on four studies: one a combined data set from Central Europe and Toronto (CETO); the other a combined data set from Germany and MD Anderson (GRMD). We searched the literature for pathway analysis methods that were widely used, representative of other methods, and had available software for performing analysis. We selected the programs EASE, which uses a modified Fishers Exact calculation to test for pathway associations, GenGen (a version of Gene Set Enrichment Analysis (GSEA)), which uses a Kolmogorov-Smirnov-like running sum statistic as the test statistic, and SLAT, which uses a p-value combination approach. We also included a modified version of the SUMSTAT method (mSUMSTAT), which tests for association by averaging chi(2) statistics from genotype association tests. There were nearly 18000 genes available for analysis, following mapping of more than 300,000 SNPs from each data set. These were mapped to 421 GO level 4 gene sets for pathway analysis. Among the methods designed to be robust to biases related to gene size and pathway SNP correlation (GenGen, mSUMSTAT and SLAT), the mSUMSTAT approach identified the most significant pathways (8 in CETO and 1 in GRMD). This included a highly plausible association for the acetylcholine receptor activity pathway in both CETO (FDR <= 0.001) and GRMD (FDR = 0.009), although two strong association signals at a single gene cluster (CHRNA3-CHRNA5-CHRNB4) drive this result, complicating its interpretation. Few other replicated associations were found using any of these methods. Difficulty in replicating associations hindered our comparison, but results suggest mSUMSTAT has advantages over the other approaches, and may be a useful pathway analysis tool to use alongside other methods such as the commonly used GSEA (GenGen) approach.