Association studies among admixed populations pose many challenges including confounding of genetic effects due to population substructure and heterogeneity due to different patterns of linkage disequilibrium (LD). We use simulations to investigate controlling for confounding by indicators of global ancestry and the impact of including a covariate for local ancestry. In addition, we investigate the use of an interaction term between a single-nucleotide polymorphism (SNP) and local ancestry to capture heterogeneity in SNP effects. Although adjustment for global ancestry can control for confounding, additional adjustment for local ancestry may increase power when the induced admixture LD is in the opposite direction as the LD in the ancestral population. However, if the induced LD is in the same direction, there is the potential for reduced power because of overadjustment. Furthermore, the inclusion of a SNP by local ancestry interaction term can increase power when there is substantial differential LD between ancestry populations. We examine these approaches in genome-wide data using the University of Southern California's Children's Health Study investigating asthma risk. The analysis highlights rs10519951 (P = 8.5 × 10−7), a SNP lacking any evidence of association from a conventional analysis (P = 0.5).
confounding; genetic association studies; genome-wide association studies; heterogeneity; linkage disequilibrium; population stratification
CYP2B6 variation predicts pharmacokinetic characteristics of its substrates. Consideration for underlying genetic structure is critical to protect against spurious associations with the highly polymorphic CYP2B6 gene.
The effect of CYP2B6 variation on response to its substrates, nonnucleoside reverse transcriptase inhibitors (NNRTIs), was explored in the Women's Interagency HIV Study.
Five putative functional polymorphisms were tested for associations with virologic suppression within one year after NNRTI initiation in women naïve to antiretroviral agents (n=91). Principal components (PCs) were generated to control for population substructure. Logistic regression was used to test the joint effect of rs3745274 and rs28399499, which together indicate slow, intermediate, and extensive metabolizers.
Rs3745274 was significantly associated with virologic suppression (OR=3.61, 95% CI 1.16-11.22, p trend=0.03); the remaining polymorphisms tested were not significantly associated with response. Women classified as intermediate and slow metabolizers were 2.90 (95% CI 0.79-12.28) and 13.44 (95% CI 1.66-infinity) times as likely to achieve virologic suppression compared to extensive metabolizers after adjustment for PCs (p trend=0.005). Failure to control for genetic ancestry resulted in substantial confounding of the relationship between the metabolizer phenotype and treatment response.
The CYP2B6 metabolizer phenotype was significantly associated with virologic response to NNRTIs; this relationship would have been masked by simple adjustment for self-reported ethnicity. Given the appreciable genetic heterogeneity that exists within self-reported ethnicity, these results exemplify the importance of characterizing underlying genetic structure in pharmacogenetic studies. Further follow-up of the CYP2B6 metabolizer phenotype is warranted given the potential clinical importance of this finding.
CYP2B6; population substructure; women; NNRTIs; confounding
We are interested in investigating the involvement of multiple rare variants within a given region by conducting analyses of individual regions with two goals: (1) to determine if regional rare variation in aggregate is associated with risk; and (2) conditional upon the region being associated, to identify specific genetic variants within the region that are driving the association. In particular, we seek a formal integrated analysis that achieves both of our goals. For rare variants with low minor allele frequencies, there is very little power to statistically test the null hypothesis of equal allele or genotype counts for each variant. Thus, genetic association studies are often limited to detecting association within a subset of the common genetic markers. However, it is very likely that associations exist for the rare variants that may not be captured by the set of common markers. Our framework aims at constructing a risk index based on multiple rare variants within a region. Our analytical strategy is novel in that we use a Bayesian approach to incorporate model uncertainty in the selection of variants to include in the index as well as the direction of the associated effects. Additionally, the approach allows for inference at both the group and variant-specific levels. Using a set of simulations, we show that our methodology has added power over other popular rare variant methods to detect global associations. In addition, we apply the approach to sequence data from the WECARE Study of second primary breast cancers.
genetic association studies; Bayesian model uncertainty; Bayes factors; multiplicity correction; sequence analysis; WECARE
Evaluate nicotinic acetycholine receptor (nAChR) single nucleotide polymorphism (SNP) association with seven day point prevalence abstinence (abstinence) in randomized clinical trials of smoking cessation therapies (RCTs) in individuals grouped by pharmacotherapy randomization to inform the development of personalized smoking cessation therapy.
We quantified association of four SNPs at three nAChRs with abstinence in eight RCTs. Participants were 2,633 outpatient treatment-seeking, self-identified European ancestry individuals smoking ≥10 cigarettes per day, recruited via advertisement, prescribed pharmacotherapy, and provided with behavioral therapy. Interventions included nicotine replacement therapy (NRT), bupropion, varenicline, placebo or combined NRT and bupropion, and five modes of group and individual behavioral therapy. Outcome measures tested in multivariate logistic regression were end of treatment (EOT) and six month (6MO) abstinence, with demographic, behavioral and genetic covariates.
“Risk” alleles previously associated with smoking heaviness were significantly (P<0.05) associated with reduced abstinence in the placebo pharmacotherapy group (PG) at 6MO [for rs588765 OR (95%CI) 0.41 (0.17–0.99)], and at EOT and at 6MO [for rs1051730, 0.42 (0.19–0.93) and 0.31 (0.12–0.80)], and with increased abstinence in the NRT PG at 6MO [for rs588765 2.07 (1.11–3.87) and for rs1051730 2.54 (1.29–4.99)]. We observed significant heterogeneity in rs1051730 effects (F=2.48, P=0.021) between PGs.
chr15q25.1 nAChR SNP risk alleles for smoking heaviness significantly increase relapse with placebo treatment and significantly increase abstinence with NRT. These SNP-PG associations require replication in independent samples for validation, and testing in larger sample sizes to evaluate whether similar effects occur in other PGs.
logistic regression; mediation analysis; nAChR variation; nicotine dependence; pharmacotherapy; randomized clinical trials
In a variety of taxa, males deploy alternative reproductive tactics to secure fertilizations. In many species, small “sneaker” males attempt to steal fertilizations while avoiding encounters with larger, more aggressive, dominant males. Sneaker males usually face a number of disadvantages, including reduced access to females and the higher likelihood that upon ejaculation, their sperm face competition from other males. Nevertheless, sneaker males represent an evolutionarily stable strategy under a wide range of conditions. Game theory suggests that sneaker males compensate for these disadvantages by investing disproportionately in spermatogenesis, by producing more sperm per unit body mass (the “fair raffle”) and/or by producing higher quality sperm (the “loaded raffle”). Here, we test these models by competing sperm from sneaker “jack” males against sperm from dominant “hooknose” males in Chinook salmon. Using two complementary approaches, we reject the fair raffle in favor of the loaded raffle and estimate that jack males were ∼1.35 times as likely as hooknose males to fertilize eggs under controlled competitive conditions. Interestingly, the direction and magnitude of this skew in paternity shifted according to individual female egg donors, suggesting cryptic female choice could moderate the outcomes of sperm competition in this externally fertilizing species.
Hooknose; jack; salmon; sexual selection; sneaker male; sperm competition
Although previous investigations have indicated a role for genetic factors in smoking initiation, the underlying genetic mechanisms are still unknown. In 2,339 adolescents from a Chinese Han population in the Wuhan Smoking Prevention Trial (Wuhan, China, 1998–1999), the authors explored the association of 57 genes in the dopamine pathway with smoking initiation. Using a conservative approach for declaring significance, positive findings were further examined in an independent sample of 603 Caucasian adolescents followed for up to 10 years as part of the Children's Health Study (Southern California, 1993–2009). The authors identified 1 single nucleotide polymorphism (rs2298122) in the calcyon neuron-specific vesicular protein gene (CALY) that was positively associated with smoking initiation in females (odds ratio = 2.21, 95% confidence interval: 1.49, 3.27; P = 8.4 × 10−5) in the Wuhan Smoking Prevention Trial cohort, and they replicated the association in females from the Children's Health Study cohort (hazard rate ratio = 2.05, 95% confidence interval: 1.27, 3.31; P = 0.003). These results suggest that the CALY gene may influence smoking initiation in adolescents, although the potential roles of underlying psychological characteristics that may be components of the smoking-initiation phenotype, such as impulsivity or novelty-seeking, remain to be explored.
adolescent; dopamine; genetic association studies; smoking
Nicotine metabolism and genetic variation have an impact on nicotine addiction and smoking abstinence, but further research is required. The nicotine metabolite ratio (NMR) is a robust biomarker of nicotine metabolism used to categorize slow and normal nicotine metabolizers (lower 25th quartile cutoff). In two randomized clinical trials of smoking abstinence treatments, we conducted NMR-stratified analyses on smoking abstinence across 13 regions coding for nicotinic acetylcholine receptors and proteins involved in the dopamine reward system. Gene × NMR interaction P-values were adjusted for multiple correlated tests, and we used a Bonferroni-corrected α-level of 0.004 to determine system-wide significance. Three SNPs in DRD1 (rs11746641, rs2168631, rs11749035) had significant interactions (0.001 ≤ adjusted P-values ≤ 0.004), with increased odds of abstinence within slow metabolizers (ORs=3.1–3.5, 95% CI 1.7–6.7). Our findings support the role of DRD1 in nicotine dependence, and identify genetic and nicotine metabolism profiles that may interact to impact nicotine dependence.
Genetic association studies; heterogeneity; smoking abstinence; nicotine metabolism; nicotine metabolite ratio; DRD1
We conducted gender-stratified analyses on a systems-based candidate gene study of 53 regions involved in nicotinic response and the brain-reward pathway in two randomized clinical trials of smoking cessation treatments (placebo, bupropion, transdermal and nasal spray nicotine replacement therapy). We adjusted P-values for multiple correlated tests, and used a Bonferroni corrected α-level of 5 × 10−4 to determine system-wide significance. Four SNPs (rs12021667, rs12027267, rs6702335, rs12039988; r2>0.98) in erythrocyte membrane protein band 4.1 (EPB41) had a significant male-specific marginal association with smoking abstinence (OR=0.5; 95% CI 0.3–0.6) at end of treatment (adjusted P<6 × 10−5). rs806365 in cannabinoid receptor 1 (CNR1) had a significant male-specific gene-treatment interaction at 6-month follow-up (adjusted P=3.9 × 10−5); within males using nasal spray, rs806365 was associated with a decrease in odds of abstinence (OR=0.04; 95% CI 0.01–0.2). While the role of CNR1 in substance abuse has been well studied, we report EPB41 for the first time in the nicotine literature.
Genetic association studies; heterogeneity; smoking cessation
There is little information regarding associations between suspected bladder cancer risk factors and tumor subtypes at diagnosis. Some, but not all, studies have found that bladder cancer among smokers is often more invasive than it is among nonsmokers. This population-based case-control study was conducted in Los Angeles, California, involving 1,586 bladder cancer patients and their individually matched controls. Logistic regression was used to conduct separate analyses according to tumor subtypes defined by stage and grade. Cigarette smoking increased risk of both superficial and invasive bladder cancer, but the more advanced the stage, the stronger the effect. The odds ratios associated with regular smokers were 2.2 (95% confidence intervals, 1.8-2.8), 2.7 (2.1-3.6) and 3.7 (2.5-5.5) for low-grade superficial, high-grade superficial and invasive tumors respectively. This pattern was consistently observed regardless of the smoking exposure index under examination. Women had higher risk of invasive bladder cancer than men even they smoked comparable amount of cigarettes as men. There was no gender difference in the association between smoking and risk of low-grade superficial bladder cancer. The heterogeneous effect of cigarette smoking was attenuated among heavy users of NSAIDs. Our results indicate that cigarette smoking was more strongly associated with increased risk of invasive bladder cancer than with low-grade superficial bladder cancer.
cigarette smoking; bladder cancer; tumor subtypes; non-steroidal anti-inflammatory drugs; Los Angeles
Colorectal cancer is the second leading cause of cancer death in developed countries. Genome-wide association studies (GWAS) have successfully identified novel susceptibility loci for colorectal cancer. To follow-up on these findings, and try to identify novel colorectal cancer susceptibility loci, we present results for genome-wide association studies (GWAS) of colorectal cancer (2,906 cases, 3,416 controls) that have not previously published main associations. Specifically, we calculated odds ratios (ORs) and 95% confidence intervals (CIs) using log-additive models for each study. In order to improve our power to detect novel colorectal cancer susceptibility loci, we performed a meta-analysis combining the results across studies. We selected the most statistically significant single nucleotide polymorphisms (SNPs) for replication using 10 independent studies (8,161 cases and 9,101 controls). We again used a meta-analysis to summarize results for the replication studies alone, and for a combined analysis of GWAS and replication studies. We measured 10 SNPs previously identified in colorectal cancer susceptibility loci and found eight to be associated with colorectal cancer (p-value range: 0.02 to 1.8 × 10−8). When we excluded studies that have previously published on these SNPs, five SNPs remained significant at p<0.05 in the combined analysis. No novel susceptibility loci were significant in the replication study after adjustment for multiple testing, and none reached genome-wide significance from a combined analysis of GWAS and replication. We observed marginally significant evidence for a second independent SNP in the BMP2 region at chromosomal location 20p12 (rs4813802; replication p-value 0.03; combined p-value 7.3 × 10−5). In a region on 5p33.15, which includes the coding regions of the TERT-CLPTM1L genes and has been identified in GWAS to be associated with susceptibility to at least seven other cancers, we observed a marginally significant association with rs2853668 (replication p-value 0.03; combined p-value 1.9 × 10−4). Our study suggests a complex nature of the contribution of common genetic variants to risk for colorectal cancer.
We utilized a cohort of 828 treatment seeking self-identified white cigarette smokers (50% female) to rank candidate gene single nucleotide polymorphisms (SNPs) associated with the Fagerström Test for Nicotine Dependence (FTND), a measure of nicotine dependence which assesses quantity of cigarettes smoked and time- and place-dependent characteristics of the respondent’s smoking behavior. 1123 SNPs at 55 autosomal candidate genes, nicotinic acetylcholine receptors and genes involved in dopaminergic function, were tested for association to baseline FTND scores adjusted for age, depression, education, sex and study site. SNP P values were adjusted for the number of transmission models, the number of SNPs tested per candidate gene, and their intragenic correlation. DRD2, SLC6A3 and NR4A2 SNPs with adjusted P values < 0.10 were considered sufficiently noteworthy to justify further genetic, bioinformatic and literature analyses. Each independent signal among the top-ranked SNPs accounted for ~1% of the FTND variance in this sample. The DRD2 SNP appears to represent a novel association with nicotine dependence. The SLC6A3 SNPs have previously been shown to be associated with SLC6A3 transcription or dopamine transporter density in vitro, in vivo and ex vivo. Analysis of SLC6A3 and NR4A2 SNPs identified a statistically significant gene-gene interaction (P=0.001), consistent with in vitro evidence that the NR4A2 protein product (NURR1) regulates SLC6A3 transcription. A community cohort of N=175 multiplex ever smoking pedigrees (N=423 ever smokers) provided nominal evidence for association with the FTND at these top ranked SNPs, uncorrected for multiple comparisons.
dopamine transporter; Fagerström Test for Nicotine Dependence; single nucleotide polymorphism; candidate gene association scan; gene-gene interaction
The conventional method of detecting gene-environment interactions, the case-control analysis, suffers from low statistical power. In contrast, the case-only analysis/design can be powerful in certain scenarios, although violation of the assumption of independence between the genetic and environmental factors can greatly bias the results. As an alternative, Bayes model averaging may be used to combine the case-control and case-only analyses. This approach first frames the case-control and case-only analyses as variations of a log-linear model. The weighting between these 2 models is then a function of the data and prior beliefs on the independence of the 2 potentially interacting factors. In this paper, the authors demonstrate via simulations that when there is no prior information on the independence of the genetic and environmental factors, this approach tends to be more powerful than the case-control analysis. Additionally, when the genetic and environmental factors are not independent in the population, bias is substantially reduced, with a corresponding reduction in type I error in comparison with the case-only analysis. Increased power or increased robustness to violations of the independence assumption may be obtained with more appropriate prior specification. The authors use an example data analysis to demonstrate the advantages of this approach.
Bayesian estimation; Bayesian model; case-control studies; epidemiologic methods; interaction
Asthma is a common disease with a complex risk architecture including both genetic and environmental factors. We performed a meta-analysis of North American genome-wide association studies (GWAS) of asthma in 5,416 asthma cases representing European Americans, African Americans/African Caribbeans, and Latinos, and replicated five regions among the most significant signals in 12,649 individuals from the same ethnic groups. Four were at previously reported loci on 17q21, and near the IL1RL1, TSLP, and IL33, genes, but we report for the first time that these loci are associated with asthma risk in three ethnic groups. In addition, we identified a novel association with asthma in the PYHIN1, gene that was specific to individuals of African descent (p=3.9×10−9). These results suggest that some asthma susceptibility loci are robust to differences in ancestry when sufficiently large samples sizes are investigated, and that ancestry-specific associations also contribute to the complex genetic architecture of asthma.
To assess the impact of change in visual field (VF) on change in health related quality of life (HRQoL) at the population level.
Prospective cohort study
3,175 Los Angles Latino Eye Study (LALES) participants
Objective measures of VF and visual acuity and self-reported HRQoL were collected at baseline and 4-year follow-up. Analysis of covariance was used to evaluate mean differences in change of HRQoL across severity levels of change in VF and to test for effect modification by covariates.
Main outcome measures
General and vision-specific HRQoL.
Of 3,175 participants, 1430 (46%) showed a change in VF (≥1 decibel [dB]) and 1651, 1715 (54%) reported a clinically important change (≥5 points) in vision-specific HRQoL. Progressive worsening and improvement in the VF were associated with increasing losses and gains in vision-specific HRQoL for the composite score and 10 of its 11 subscales (all Ptrends<0.05). Losses in VF > 5 dB and gains > 3 dB were associated with clinically meaningful losses and gains in vision-specific HRQoL, respectively. Areas of vision-specific HRQoL most affected by greater losses in VF were driving, dependency, role-functioning, and mental health. The effect of change in VF (loss or gain) on mean change in vision-specific HRQoL varied by level of baseline vision loss (in visual field and/or visual acuity) and by change in visual acuity (all P-interactions<0.05). Those with moderate/severe VF loss at baseline and with a > 5 dB loss in visual field during the study period had a mean loss of vision-specific HRQoL of 11.3 points, while those with no VF loss at baseline had a mean loss of 0.97 points Similarly, with a > 5 dB loss in VF and baseline visual acuity impairment (mild/severe) there was a loss in vision-specific HRQoL of 10.5 points, whereas with no visual acuity impairment at baseline there was a loss of vision-specific HRQoL of 3.7 points.
Both losses and gains in VF produce clinically meaningful changes in vision-specific HRQoL. In the presence of pre-existing vision loss (VF and visual acuity), similar levels of visual field change produce greater losses in quality of life.
The insulin-like growth factor (IGF) signaling axis plays an important role in cancer biology. We hypothesized that genetic variation in this pathway may influence risk of ovarian cancer. A three-center study of non-Hispanic whites including 1880 control women, 1135 women with invasive epithelial ovarian cancer and 321 women with borderline epithelial ovarian tumors was carried out to test the association between tag single-nucleotide polymorphisms (tSNPs) (n=58) in this pathway and risk of ovarian cancer. We found no association between variation in IGF1, IGFBP1 or IGFBP3 and risk of invasive disease, whereas five tSNPs in IGF2 were associated with risk of invasive epithelial ovarian cancer at P<0.05 and followed-up one of the associated SNPs. We conducted genotyping in 3216 additional non-Hispanic white cases and 5382 additional controls and were able to independently replicate our initial findings. In the combined set of studies, rs4320932 was associated with a 13% decreased risk of ovarian cancer per copy of the minor allele carried (95% confidence interval 0.81–0.93, P-trend=7.4 × 10−5). No heterogeneity of effect across study centers was observed (phet=0.25). IGF2 is emerging as an important gene for ovarian cancer; additional genotyping is warranted to further confirm these associations with IGF2 and to narrow down the region harboring the causal SNP.
Rare deleterious mutations in BRCA1 and BRCA2 are associated with an elevated risk of breast and ovarian cancer. Whether or not common variants in these genes are independently associated with risk of breast cancer remains unclear. In this study, we included 632 Caucasian women with asynchronous contralateral breast cancer (CBC, cases) and 1,221 women with unilateral breast cancer (controls) from the WECARE (Women’s Environment, Cancer And Radiation Epidemiology) Study. BRCA1 and BRCA2 deleterious mutation status was measured using denaturing high-performance liquid chromatography followed by direct sequencing, yielding including 88 BRCA1 and 60 BRCA2 deleterious mutation carriers. We also genotyped samples on the Illumina Omni1-Quad platform. We assessed the association between CBC risk and common (minor allele frequency>0.05) single nucleotide polymorphisms (SNPs) in BRCA1 (n SNPs=22) and BRCA2 (n SNPs=30) and haplotypes using conditional logistic regression accounting for BRCA1/BRCA2 mutation status. We found no significant associations between any single-SNPs or haplotypes of BRCA1 or BRCA2 and risk of CBC among all women. When we stratified by BRCA1 and BRCA2 mutation carrier status, we found suggestive evidence that risk estimates for selected SNPs in BRCA1 (rs8176318, rs1060915 and rs16940) and BRCA2 (rs11571686, rs206115 and rs206117) may differ in non-carriers and carriers of deleterious mutations in BRCA1 and BRCA2. One common haplotype on BRCA1 was inversely significantly associated with risk only among non-BRCA1 and BRCA2 carriers. The association between common variants in BRCA1 and BRCA2 and risk of CBC may differ depending on BRCA1 and BRCA2 mutation carrier status.
BRCA1; BRCA2; haplotype; polymorphism; breast cancer; contralateral
Genome-wide association studies (GWAS) have led to the identification of a number of common susceptibility loci for colorectal cancer (CRC); however, none of these GWAS have considered gene-environment (GxE) interactions. Therefore, it is unclear whether current hits are modified by environmental exposures or whether there are additional hits whose effects are dependent on environmental exposures.
We conducted a systematic search for GxE interactions using genome wide data from the Colon Cancer Family Registry that included 1,191 cases of microsatellite stable (MSS) or microsatellite instability (MSI)-low CRC and 999 controls genotyped using either the Illumina Human1M or Human1M-Duo BeadChip. We tested for interactions between genotypes and 14 environmental factors using three methods: a traditional case-control test, a case-only test, and the recently proposed two-step method by Murcray et al. All potentially significant findings were replicated in the ARCTIC Study.
No GxE interactions were identified that reached genome-wide significance by any of the three methods. When analyzing previously reported susceptibility loci, seven significant GxE interactions were found at a 5% significance level. We investigated these seven interactions in an independent sample and none of the interactions were replicated.
Identifying GxE interactions will present challenges in a GWAS setting. Our power calculations illustrate the need for larger sample sizes; however, since CRC is a heterogeneous disease, a tradeoff between increasing sample size and heterogeneity needs to be considered.
The results from this first genome-wide analysis of GxE in CRC identify several challenges, which may be addressed by large consortium efforts.
GWAS; colorectal cancer; gene-environment interactions
BACKGROUND & AIMS
Heritable factors contribute to the development of colorectal cancer. Identifying the genetic loci associated with colorectal tumor formation could elucidate the mechanisms of pathogenesis.
We conducted a genome-wide association study that included 14 studies, 12,696 cases of colorectal tumors (11,870 cancer, 826 adenoma), and 15,113 controls of European descent. The 10 most statistically significant, previously unreported findings were followed up in 6 studies; these included 3056 colorectal tumor cases (2098 cancer, 958 adenoma) and 6658 controls of European and Asian descent.
Based on the combined analysis, we identified a locus that reached the conventional genome-wide significance level at less than 5.0 × 10−8: an intergenic region on chromosome 2q32.3, close to nucleic acid binding protein 1 (most significant single nucleotide polymorphism: rs11903757; odds ratio [OR], 1.15 per risk allele; P = 3.7 × 10−8). We also found evidence for 3 additional loci with P values less than 5.0 × 10−7: a locus within the laminin gamma 1 gene on chromosome 1q25.3 (rs10911251; OR, 1.10 per risk allele; P = 9.5 × 10−8), a locus within the cyclin D2 gene on chromosome 12p13.32 (rs3217810 per risk allele; OR, 0.84; P = 5.9 × 10−8), and a locus in the T-box 3 gene on chromosome 12q24.21 (rs59336; OR, 0.91 per risk allele; P = 3.7 × 10−7).
In a large genome-wide association study, we associated polymorphisms close to nucleic acid binding protein 1 (which encodes a DNA-binding protein involved in DNA repair) with colorectal tumor risk. We also provided evidence for an association between colorectal tumor risk and polymorphisms in laminin gamma 1 (this is the second gene in the laminin family to be associated with colorectal cancers), cyclin D2 (which encodes for cyclin D2), and T-box 3 (which encodes a T-box transcription factor and is a target of Wnt signaling to β-catenin). The roles of these genes and their products in cancer pathogenesis warrant further investigation.
Colon Cancer; Genetics; Risk Factors; SNP
Many complex diseases are likely to be a result of the interplay of genes and environmental exposures. The standard analysis in a genome-wide association study (GWAS) scans for main effects and ignores the potentially useful information in the available exposure data. Two recently proposed methods that exploit environmental exposure information involve a two-step analysis aimed at prioritizing the large number of SNPs tested to highlight those most likely to be involved in a G×E interaction. For example, Murcray et al (2009) proposed screening on a test that models the G-E association induced by an interaction in the combined case-control sample. Alternatively, Kooperberg et al (2008) suggested screening on genetic marginal effects. In both methods, SNPs that pass the respective screening step at a pre-specified significance threshold are followed up with a formal test of interaction in the second step. We propose a hybrid method that combines these two screening approaches by allocating a proportion of the overall genomewide significance level to each test. We show that the Murcray et al. approach is often the most efficient method, but that the hybrid approach is a powerful and robust method for nearly any underlying model. As an example, for a GWAS of 1 million markers including a single true disease SNP with minor allele frequency of 0.15, and a binary exposure with prevalence 0.3, the Murcray, Kooperberg and hybrid methods are 1.90, 1.27, and 1.87 times as efficient, respectively, as the traditional case-control analysis to detect an interaction effect size of 2.0.
G×E interaction; case-control; genome-wide association study; efficiency
Genome-wide association studies have associated common variations at chromosomes 5p15 and 15q25 with lung cancer risk. The 5p15 locus has also been associated with increased bladder cancer risk in a recent report. The 15q25 locus has been associated with nicotine dependence and self-reported number of cigarettes smoked per day in some studies and it was proposed that its association with lung cancer may be mediated through differences in smoking behavior. Here, we investigated the roles of variations at 5p15 (rs401681, rs402710, rs2736098 and rs2736100) and 15q25 (rs1051730 and rs8034191) in bladder cancer etiology in two case–control studies conducted separately in Los Angeles County, CA, USA (498 cases and 588 controls) and in Shanghai, China (506 cases and 530 controls). We replicated the association between the 5p15 locus and bladder cancer among non-Hispanic whites (NHW) in Los Angeles [for rs2736100, per C allele odds ratio (OR) = 1.23; 95% confidence interval (CI), 1.02–1.48; P = 0.029] and among Chinese in Shanghai (OR = 1.22; 95% CI, 1.02–1.47; P = 0.033). Both rs1051730 and rs8034191 at 15q25 were rare among Chinese. Among NHW, a significant association was found between rs8034191 and bladder cancer which persisted after adjustment for cigarette smoking status, number of cigarettes smoked per day and number of years of smoking (per C allele OR = 1.26; 95% CI, 1.04–1.54; P = 0.017). Our results support 5p15 and 15q25 as susceptibility regions for bladder cancer risk.
Despite effective therapies for smoking cessation, most smokers find quitting difficult and most successful quitters relapse. Considerable evidence supports a genetic risk for nicotine dependence; however, less is known about the pharmacogenetics of smoking cessation. In the first pharmacogenetic investigation of the efficacy of varenicline and bupropion, we examined whether genes important in the pharmacodynamics and pharmacokinetics of these drugs and nicotine predict medication efficacy and adverse events. Subjects participated in randomized, double-blind, placebo-controlled smoking cessation clinical trials, comparing varenicline, a nicotinic acetylcholine receptor (nAChR) partial agonist, with bupropion, a norepinephrine/dopamine reuptake inhibitor, and placebo. Primary analysis included 1175 smokers of European ancestry, and 785 single nucleotide polymorphisms from 24 genes, representing 254 linkage disequilibrium (LD) bins (genes included nAChR subunits, additional varenicline-specific genes, and genes involved in nicotine or bupropion metabolism). For varenicline, continuous abstinence (weeks 9–12) was associated with multiple nAChR subunit genes (including CHRNB2, CHRNA5, and CHRNA4) (OR=1.76; 95% CI: 1.23–2.52) (p<0.005); for bupropion, abstinence was associated with CYP2B6 (OR=1.78; 95% CI: 1.27–2.50) (p<0.001). Incidence of nausea was associated with several nAChR subunit genes (OR=0.50; 95% CI: 0.36–0.70) (p<0.0001) and time to relapse after quitting was associated with HTR3B (HR=1.97; 95% CI: 1.45–2.68) (p<0.0001). These data provide evidence for multiple genetic loci contributing to smoking cessation and therapeutic response. Different loci are associated with varenicline vs bupropion response, suggesting that additional research may identify clinically useful markers to guide treatment decisions.
varenicline; bupropion; pharmacogenetics; nicotine; nicotinic receptor; CYP2B6; pharmacogenetics/pharmacogenomics; addiction and substance abuse; clinical pharmacology/clinical trials; neuropharmacology; varenicline; bupropion; nicotine; smoking cessation; nicotinic receptors
Survivors of pediatric Hodgkin lymphoma (HL) are at significant risk for radiation therapy (RT)-induced second malignant neoplasms (SMNs). We identified two variants at chromosome 6q21 associated with SMNs in HL survivors treated with RT as children but not as adults. The variants comprise a risk locus associated with decreased basal PRDM1 expression and impaired induction of PRDM1 by radiation exposure. These data suggest a novel gene-exposure interaction that may implicate PRDM1 in the etiology of RT-induced SMNs.
The aim of this study is to investigate the relationships between hypertension, hypertension medication and bladder cancer risk in a population-based case–control study conducted in Los Angeles. Non-Asians between the ages of 25 and 64 years with histologically confirmed bladder cancers diagnosed between 1987 and 1996 were identified through the Los Angeles County Cancer Surveillance Program. A total of 1585 cases and their age-, gender- and race-matched neighborhood controls were included in the analyses. Conditional logistic regression models were used to examine the relationship between history of hypertension, medication use and bladder cancer risk. A history of hypertension was not related to bladder cancer; however, among hypertensive individuals, there was a significant difference in bladder cancer risk related to the use of diuretics or antihypertensive drugs (P for heterogeneity = 0.004). Compared with individuals without hypertension, hypertensive individuals who regularly used diuretics/antihypertensives had a similar risk [odds ratio (OR) 1.06; 95% confidence interval (CI) 0.86–1.30], whereas untreated hypertensive subjects had a 35% reduction in risk (OR: 0.65; 95% CI: 0.48–0.88). A greater reduction in bladder cancer risk was observed among current-smokers (OR: 0.43; 95% CI: 0.27–0.71) and carriers of GSTM1-null (homozygous absence) genotypes (OR: 0.43; 95% CI: 0.22–0.85). Similarly, among smokers with GSTM1-null genotype, levels of 4-aminobiphenyl-hemoglobin adducts were significantly lower among untreated hypertensive individuals (45.7 pg/g Hb) compared with individuals without hypertension (79.8 pg/g Hb) (P = 0.009). In conclusion, untreated hypertension was associated with a reduced risk of bladder cancer.
Unraveling complex interactions has been a challenge in epidemiologic research. We introduce a pathway modeling framework that discovers plausible pathways from observational data, and allows estimation of both the net effect of the pathway and the types of interactions occurring among genetic or environmental risk factors. Each discovered pathway structure links combinations of observed variables through intermediate latent nodes to a final node, the outcome. Biologic knowledge can be readily applied in this framework as a prior on pathway structure to give preference to more biologically plausible models, thereby providing more precise estimation of Bayes factors for pathways of greatest interest by Markov Chain Monte Carlo (MCMC) methods.
Data was simulated for binary inputs of which only a subset was involved in different pathway topologies. Our algorithm was then used to recover the pathway from the simulated data. The posterior distributions of inputs, pair-wise and higher order interactions, and topologies were obtained by MCMC methods. The evidence in favor of a particular pathway or interaction was summarized using Bayes factors. Our method can correctly identify the risk factors and interactions involved in the simulated pathway. We apply our framework to an asthma case-control dataset with polymorphisms in 12 genes.
Genetic variants in the calcium/vitamin D metabolic pathway may be related to risk for colorectal cancer. While several investigations of vitamin D receptor (VDR) polymorphisms and colorectal cancer have been conducted, no studies to date have evaluated the association of genetic variation in the heterodimer partner for VDR, the retinoid X receptor (RXR). Another important gene in this pathway is the calcium-sensing receptor (CASR). Employing a discordant-sibship case–control design, we examined the association between single nucleotide polymorphisms (SNPs) in RXRA and CASR and risk for colorectal cancer overall and by colorectal subsite and microsatellite instability (MSI) status using data from the Colon Cancer Family Registry. No gene-level relationships between RXRA or CASR and colorectal cancer overall were observed. However, for RXRA SNP rs7861779, a high-interest SNP selected for study a priori, there was a statistically significantly increased risk for proximal colorectal cancer among those with at least one A allele [odds ratio (OR) = 1.42; 95% confidence interval (CI) = 1.03–1.97]. Another selected RXRA SNP, rs12004589, was significantly associated with risk of MSI-high cancers (OR = 2.27; 95% CI = 1.13–4.56). Additionally, CASR SNP rs1801726 was significantly associated with a reduced risk for rectal cancer (OR = 0.53; 95% CI = 0.29–0.96). These results provide support that RXRA SNPs rs7861779 and rs12004589 and CASR SNP rs1801726 may be important markers for colorectal neoplasia. Further work is needed to elucidate their role in the carcinogenic pathway.