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
There is increasing evidence that response to pharmacological treatment for nicotine dependence may be moderated by genetic polymorphisms. However, the feasibility, acceptability, and impact of genetically tailoring treatment in real-world clinical settings are unknown.
We conducted a multiphased, mixed-methods feasibility study with current smokers to develop and evaluate a patient-centered, theoretically grounded personalized medicine treatment protocol. The initial research phase included formative work to develop intervention materials. The second phase included a randomized pilot trial to evaluate the intervention. Trial participants (n = 36) were genotyped for ANKK1 rs1800497 and were randomized to receive genetic feedback (GF) plus standard behavioral counseling (BC) for smoking cessation or BC without GF. All participants received genetically tailored pharmacotherapy (nicotine patch or bupropion).
The intervention was feasible to implement and was acceptable to participants based on satisfaction ratings and objective measures of participation. There was no evidence that the GF resulted in adverse psychological outcomes (e.g., depression, fatalism, reduced perceived control over quitting, differential motivation for quitting) based on quantitative or qualitative outcomes.
Study results suggest that it is feasible to offer treatment within a health care setting that includes genetically tailored pharmacotherapy and doing so had no apparent adverse psychological impacts. Further evaluation of pharmacogenetically tailored smoking cessation interventions appears warranted.
We investigated whole saliva as a source of biomarkers to distinguish individuals who have, and who have not, been chronically exposed to severe and threatening life difficulties. We evaluated RNA and DNA metrics, expression of 37 candidate genes, and cortisol release in response to the Trier Social Stress Test, as well as clinical characteristics, from 48 individuals stratified on chronic exposure to psychosocial stressors within the last year as measured by the Life Events and Difficulties Schedule. Candidate genes were selected based on their differential gene expression ratio in circulating monocytes from a published genome-wide analysis of adults experiencing different levels of exposure to a chronic stressor.
In univariate analyses, we observed significantly decreased RNA integrity (RIN) score (P = 0.04), and reduced expression of glucocorticoid receptor-regulated genes (Ps < 0.05) in whole saliva RNA from individuals exposed to chronic stressors, as compared to those with no exposure. In those exposed, we observed significantly decreased BMI (P < 0.001), increased ever-smoking and increased lifetime alcohol abuse or dependence (P ≤ 0.03), and a reduction of cortisol release. In post hoc multivariate analyses including clinical and biospecimen-derived variables, we consistently observed significantly decreased expression of IL8 (Ps < 0.05) in individuals exposed, with no significant association to RIN score. Alcohol use disorders, tobacco use, a reduced acute stress response and decreased salivary IL8 gene expression characterize emerging adults chronically exposed to severe and threatening psychosocial stressors.
Human; Saliva; Gene expression; IL8; qPCR
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
To date, most genetic association studies of tobacco use have been conducted in European American subjects using the phenotype of smoking quantity (cigarettes per day). However, smoking quantity is a very imprecise measure of exposure to tobacco smoke constituents. Analyses of alternate phenotypes and populations may improve our understanding of tobacco addiction genetics. Cotinine is the major metabolite of nicotine, and measuring serum cotinine levels in smokers provides a more objective measure of nicotine dose than smoking quantity. Previous genetic association studies of serum cotinine have focused on individual genes. We conducted a genetic association study of the biomarker in African American (N=365) and European American (N=315) subjects from the Coronary Artery Risk Development in Young Adults study using a chip containing densely-spaced tag SNPs in ∼2100 genes. We found that rs11187065, located in the non-coding region (intron 1) of insulin-degrading enzyme (IDE), was the most strongly associated SNP (p=8.91 × 10−6) in the African American cohort, whereas rs11763963, located on chromosome 7 outside of a gene transcript, was the most strongly associated SNP in European Americans (p=1.53 × 10−6). We then evaluated how the top variant association in each population performed in the other group. We found that the association of rs11187065 in IDE was also associated with the phenotype in European Americans (p=0.044). Our top SNP association in European Americans, rs11763963 was non-polymorphic in our African American sample. It has been previously shown that psychostimulant self-administration is reduced in animals with lower insulin because of interference with dopamine transmission in the brain reward centers. Our finding provides a platform for further investigation of this, or additional mechanisms, involving the relationship between insulin and self-administered nicotine dose.
cotinine; nicotine; IDE; CARDIA; MORF4L1; IMAT-Broad-CARe; behavioral science; neurogenetics; addiction & substance abuse; pharmacogenetics/pharmacogenomics; cotinine; nicotine; IDE; CARDIA; MORF4L1
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
Nicotine withdrawal symptoms are related to smoking cessation. A Rasch model has been used to develop a unidimensional sensitivity score representing multiple correlated measures of nicotine withdrawal. A previous autosome-wide screen identified a nonparametric linkage (NPL) log-likelihood ratio (LOD) score of 2.7 on chromosome 6q26 for the sum of nine withdrawal symptoms.
The objectives of these analyses are: a) to assess the influence of nicotine withdrawal sensitivity on relapse, b) conduct autosome-wide NPL analysis of nicotine withdrawal sensitivity among 158 pedigrees with 432 individuals with microsatellite genotypes and nicotine withdrawal scores, and c) explore family-based association of single nucleotide polymorphism (SNPs) at the mu opioid receptor (MOR) candidate gene (OPRM1) to nicotine withdrawal sensitivity in 172 nuclear pedigrees with 419 individuals with both SNP genotypes and nicotine withdrawal scores.
An increased risk for relapse was associated with nicotine withdrawal sensitivity score (odds ratio, OR=1.25, 95% confidence interval, 95%CI=1.10,1.42). A maximal NPL LOD score of 3.15, suggestive of significant linkage, was identified at chr6q26 for nicotine withdrawal sensitivity. Evaluation of 18 OPRM1 SNPs via the family based association test (FBAT) with the nicotine withdrawal sensitivity score identified eight tagging SNPs with global P-values<0.05 and false discovery rate Q-values<0.06.
An increased risk of relapse, suggestive linkage at chr6q26, and nominally significant association with multiple OPRM1 SNPs was found with Rasch modeled nicotine withdrawal sensitivity score in a multiplex smoking pedigree sample. Future studies should attempt to replicate these findings and investigate the relationship between nicotine withdrawal symptoms and variation at OPRM1.
Follow-up studies of eating disorders (EDs) suggest outcomes ranging from recovery to chronic illness or death, but predictors of outcome have not been consistently identified. We tested 5151 single-nucleotide polymorphisms (SNPs) in approximately 350 candidate genes for association with recovery from ED in 1878 women. Initial analyses focused on a strictly defined discovery cohort of women who were over age 25 years, carried a lifetime diagnosis of an ED, and for whom data were available regarding the presence (n=361 ongoing symptoms in the past year, ie, ‘ill') or absence (n=115 no symptoms in the past year, ie, ‘recovered') of ED symptoms. An intronic SNP (rs17536211) in GABRG1 showed the strongest statistical evidence of association (p=4.63 × 10−6, false discovery rate (FDR)=0.021, odds ratio (OR)=0.46). We replicated these findings in a more liberally defined cohort of women age 25 years or younger (n=464 ill, n=107 recovered; p=0.0336, OR=0.68; combined sample p=4.57 × 10−6, FDR=0.0049, OR=0.55). Enrichment analyses revealed that GABA (γ-aminobutyric acid) SNPs were over-represented among SNPs associated at p<0.05 in both the discovery (Z=3.64, p=0.0003) and combined cohorts (Z=2.07, p=0.0388). In follow-up phenomic association analyses with a third independent cohort (n=154 ED cases, n=677 controls), rs17536211 was associated with trait anxiety (p=0.049), suggesting a possible mechanism through which this variant may influence ED outcome. These findings could provide new insights into the development of more effective interventions for the most treatment-resistant patients.
GABA; anorexia nervosa; recovery from eating disorders; genetic association; single nucleotide polymorphisms; eating/metabolic disorders; GABA; eating/metabolic disorders; neurogenetics; biological psychiatry; genetic association; anorexia nervosa; recovery from eating disorders; single-nucleotide polymorphisms; phenomic association
Affordable early screening in subjects with high risk of lung cancer has great potential to improve survival from this deadly disease. We measured gene expression from lung tissue and peripheral whole blood (PWB) from adenocarcinoma cases and controls to identify dysregulated lung cancer genes that could be tested in blood to improve identification of at-risk patients in the future. Genome-wide mRNA expression analysis was conducted in 153 subjects (73 adenocarcinoma cases, 80 controls) from the Environment And Genetics in Lung cancer Etiology (EAGLE) study using PWB and paired snap-frozen tumor and non-involved lung tissue samples. Analyses were conducted using unpaired t-tests, linear mixed effects and ANOVA models. The area under the receiver operating characteristic curve (AUC) was computed to assess the predictive accuracy of the identified biomarkers. We identified 50 dysregulated genes in stage I adenocarcinoma versus control PWB samples (False Discovery Rate ≤0.1, fold change ≥1.5 or ≤0.66). Among them, eight (TGFBR3, RUNX3, TRGC2, TRGV9, TARP, ACP1, VCAN, and TSTA3) differentiated paired tumor versus non-involved lung tissue samples in stage I cases, suggesting a similar pattern of lung cancer-related changes in PWB and lung tissue. These results were confirmed in two independent gene expression analyses in a blood-based case-control study (n=212) and a tumor-non tumor paired tissue study (n=54). The eight genes discriminated patients with lung cancer from healthy controls with high accuracy (AUC=0.81, 95% CI=0.74–0.87). Our finding suggests the use of gene expression from PWB for the identification of early detection markers of lung cancer in the future.
microarray gene expression; peripheral blood; lung cancer; stage I
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
This analysis is a follow-up to an earlier investigation of 182 genes selected as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN). As those initial case-control results revealed no statistically significant differences in single nucleotide polymorphisms, herein we investigate alternative phenotypes associated with AN. In 1762 females using regression analyses we examined: (1) lowest illness-related attained body mass index; (2) age at menarche; (3) drive for thinness; (4) body dissatisfaction; (5) trait anxiety; (6) concern over mistakes; and (7) the anticipatory worry and pessimism vs. uninhibited optimism subscale of the harm avoidance scale. After controlling for multiple comparisons, no statistically significant results emerged. Although results must be viewed in the context of limitations of statistical power, the approach illustrates a means of potentially identifying genetic variants conferring susceptibility to AN because less complex phenotypes associated with AN are more proximal to the genotype and may be influenced by fewer genes.
covariates; eating disorders; association studies; personality; genetic
The molecular drivers that determine histology in lung cancer are largely unknown. We investigated whether microRNA (miR) expression profiles can differentiate histological subtypes and predict survival for non-small cell lung cancer.
We analyzed miR expression in 165 adenocarcinoma (AD) and 125 squamous cell carcinoma (SQ) tissue samples from the Environmental And Genetics in Lung cancer Etiology (EAGLE) study using a custom oligo array with 440 human mature antisense miRs. We compared miR expression profiles using t-tests and F-tests and accounted for multiple testing using global permutation tests. We assessed the association of miR expression with tobacco smoking using Spearman correlation coefficients and linear regression models, and with clinical outcome using log-rank tests, Cox proportional hazards and survival risk prediction models, accounting for demographic and tumor characteristics.
MiR expression profiles strongly differed between AD and SQ (global p<0.0001), particularly in the early stages, and included miRs located on chromosome loci most often altered in lung cancer (e.g., 3p21-22). Most miRs, including all members of the let-7 family, were down-regulated in SQ. Major findings were confirmed by QRT-PCR in EAGLE samples and in an independent set of lung cancer cases. In SQ, low expression of miRs down-regulated in the histology comparison was associated with 1.2 to 3.6-fold increased mortality risk. A 5-miR signature significantly predicted survival for SQ.
We identified a miR expression profile that strongly differentiated AD from SQ and had prognostic implications. These findings may lead to histology-based therapeutic approaches.
The authors investigated the relation between alcohol consumption and lung cancer risk in the Environment and Genetics in Lung Cancer Etiology (EAGLE) Study, a population-based case-control study. Between 2002 and 2005, 2,100 patients with primary lung cancer were recruited from 13 hospitals within the Lombardy region of Italy and were frequency-matched on sex, area of residence, and age to 2,120 randomly selected controls. Alcohol consumption during adulthood was assessed in 1,855 cases and 2,065 controls. Data on lifetime tobacco smoking, diet, education, and anthropometric measures were collected. Adjusted odds ratios and 95% confidence intervals for categories of mean daily ethanol intake were calculated using unconditional logistic regression. Overall, both nondrinkers (odds ratio = 1.42, 95% confidence interval: 1.03, 2.01) and very heavy drinkers (≥60 g/day; odds ratio = 1.44, 95% confidence interval: 1.01, 2.07) were at significantly greater risk than very light drinkers (0.1–4.9 g/day). The alcohol effect was modified by smoking behavior, with no excess risk being observed in never smokers. In summary, heavy alcohol consumption was a risk factor for lung cancer among smokers in this study. Although residual confounding by tobacco smoking cannot be ruled out, this finding may reflect interplay between alcohol and smoking, emphasizing the need for preventive measures.
alcohol drinking; case-control studies; ethanol; lung neoplasms; risk factors; smoking
We performed association studies with 5,151 SNPs that were judged as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN) based on location under reported linkage peaks, previous results in the literature (182 candidate genes), brain expression, biological plausibility, and estrogen responsivity. We employed a case–control design that tested each SNP individually as well as haplotypes derived from these SNPs in 1,085 case individuals with AN diagnoses and 677 control individuals. We also performed separate association analyses using three increasingly restrictive case definitions for AN: all individuals with any subtype of AN (All AN: n = 1,085); individuals with AN with no binge eating behavior (AN with No Binge Eating: n = 687); and individuals with the restricting subtype of AN (Restricting AN: n = 421). After accounting for multiple comparisons, there were no statistically significant associations for any individual SNP or haplotype block with any definition of illness. These results underscore the importance of large samples to yield appropriate power to detect genotypic differences in individuals with AN and also motivate complementary approaches involving Genome-Wide Association (GWA) studies, Copy Number Variation (CNV) analyses, sequencing-based rare variant discovery assays, and pathway-based analysis in order to make up for deficiencies in traditional candidate gene approaches to AN.
single nucleotide polymorphisms; probands; anorexia nervosa; bulimia nervosa
This study addressed the hypothesis that variation in genes associated with dopamine function (SLC6A3, DRD2, DRD4), serotonin function (SLC6A4), and regulation of monoamine levels (MAOA) may be predictive of BMI categories (obese and overweight + obese) in young adulthood and of changes in BMI as adolescents transition into young adulthood. Interactions with gender and race/ethnicity were also examined.
Research Methods and Procedures
Participants were a subsample of individuals from The National Longitudinal Study of Adolescent Health (Add Health), a nationally representative sample of adolescents followed from 1995 to 2002. The sample analyzed included a subset of 1584 unrelated individuals with genotype data. Multiple logistic regressions were conducted to evaluate associations between genotypes and obesity (BMI > 29.9) or overweight + obese combined (BMI > 25) with normal weight (BMI = 18.5–24.9) as a referent. Linear regression models were used examine change in BMI from adolescence to young adulthood.
Significant associations were found between SLC6A4 5HTTLPR and categories of BMI, and between MAOA promoter VNTR among males and categories of BMI. Stratified analyses revealed that the association between these two genes and excess BMI was significant for males overall, and for White and Hispanic males specifically. Linear regression models indicated a significant effect of SLC6A4 5HTTLPR on change in BMI from adolescence to young adulthood.
Our findings lend further support to the involvement of genes implicated in dopamine and serotonin regulation on energy balance.
Adolescents; Genetic Epidemiology; Serotonin; Neuro Transmitter
Developing, targeting, and evaluating genomic strategies for population-based disease prevention require population-based data. In response to this urgent need, genotyping has been conducted within the Third National Health and Nutrition Examination (NHANES III), the nationally-representative household-interview health survey in the U.S. However, before these genetic analyses can occur, family relationships within households must be accurately ascertained. Unfortunately, reported family relationships within NHANES III households based on questionnaire data are incomplete and inconclusive with regards to actual biological relatedness of family members. We inferred family relationships within households using DNA fingerprints (Identifiler®) that contain the DNA loci used by law enforcement agencies for forensic identification of individuals. However, performance of these loci for relationship inference is not well understood. We evaluated two competing statistical methods for relationship inference on pairs of household members: an exact likelihood ratio relying on allele frequencies to an Identical By State (IBS) likelihood ratio that only requires matching alleles. We modified these methods to account for genotyping errors and population substructure. The two methods usually agree on the rankings of the most likely relationships. However, the IBS method underestimates the likelihood ratio by not accounting for the informativeness of matching rare alleles. The likelihood ratio is sensitive to estimates of population substructure, and parent-child relationships are sensitive to the specified genotyping error rate. These loci were unable to distinguish second-degree relationships and cousins from being unrelated. The genetic data is also useful for verifying reported relationships and identifying data quality issues. An important by-product is the first explicitly nationally-representative estimates of allele frequencies at these ubiquitous forensic loci.
Forensics; allele sharing; population structure; CODIS; IBS; IBD
Characterizing cotinine pharmacokinetics is a useful way to study nicotine metabolism because the same liver enzyme is primarily responsible for the metabolism of both, and the clearances of nicotine and cotinine are highly correlated. We conducted a whole-genome linkage analysis to search for candidate regions influencing quantitative variation in cotinine pharmacokinetics in a large-scale pharmacokinetic study with 61 families containing 224 healthy adult participants. The strongest linkage signal was identified at 135 cM of chromosome 9 with LOD=2.81 and P=0.0002; two other suggestive linkage peaks appear at 31.4 and 73.5 cM of chromosome 11 with LOD=1.96 (P=0.0013) and 1.94 (P=0.0014). The confidence level of the linkage between the three genome regions and cotinine pharmacokinetics is statistically significant with a genome-wide empirical probability of P=0.029.
pharmacokinetics; nicotine; dependence; linkage analysis
In the genetic study of complex traits, especially behavior related ones, such as smoking and alcoholism, usually several phenotypic measurements are obtained for the description of the complex trait, but no single measurement can quantify fully the complicated characteristics of the symptom because of our lack of understanding of the underlying etiology. If those phenotypes share a common genetic mechanism, rather than studying each individual phenotype separately, it is more advantageous to analyze them jointly as a multivariate trait in order to enhance the power to identify associated genes. We propose a multilocus association test for the study of multivariate traits. The test is derived from a partially linear tree-based regression model for multiple outcomes. This novel tree-based model provides a formal statistical testing framework for the evaluation of the association between a multivariate outcome and a set of candidate predictors, such as markers within a gene or pathway, while accommodating adjustment for other covariates. Through simulation studies we show that the proposed method has an acceptable type I error rate and improved power over the univariate outcome analysis, which studies each component of the complex trait separately with multiple-comparison adjustment. A candidate gene association study of multiple smoking-related phenotypes is used to demonstrate the application and advantages of this new method. The proposed method is general enough to be used for the assessment of the joint effect of a set of multiple risk factors on a multivariate outcome in other biomedical research settings.
Generalized estimating equation; Genetic association study; Model selection; Multiple-comparison adjustment; Tree-based model
The ratio of trans-3’hydroxycotinine/cotinine (3HC/COT) is a marker of CYP2A6 activity, an important determinant of nicotine metabolism. This analysis sought to conduct a combined genetic epidemiologic and pharmacogenetic investigation of the 3HC/COT ratio in plasma and urine.
One hundred thirty nine twin pairs (110 monozygotic [MZ] and 29 dizygotic [DZ]) underwent a 30-minute infusion of stable isotope-labeled nicotine and its major metabolite, cotinine, followed by an 8-hour in-hospital stay. Blood and urine samples were taken at regular intervals for analysis of nicotine, cotinine, and metabolites. DNA was genotyped to confirm zygosity and for variation in the gene for the primary nicotine metabolic enzyme, CYP2A6 (variants genotyped: *1B, *1×2, *2, *4, *9, *12). Univariate biometric analyses quantified genetic and environmental influences on each measure in the presence and absence of covariates, including measured CYP2A6 genotype.
There was a substantial amount of variation in the free 3HC/COT ratio in plasma (6 hours post-infusion) attributable to additive genetic influences (67.4%, 95% CI = 55.9–76.2%). The heritability estimate was reduced to 61.0% and 49.4%, respectively, after taking into account the effect of covariates and CYP2A6 genotype. In urine (collected over 8 hours), the estimated amount of variation in the 3HC/COT ratio attributable to additive genetic influences was smaller (47.2%, 95% CI = 0–67.2%) and decreased to 44.6% and 42.0% after accounting for covariates and genotype.
Additive genetic factors are prominent in determining variation in plasma 3HC/COT variation but less so in determining variation in urine 3HC/COT.
pharmacogenetics; nicotine; cotinine; metabolism; CYP2A6; twins; genetics; heritability
It is increasingly recognized that pathway analyses—a joint test of association between the outcome and a group of single nucleotide polymorphisms (SNPs) within a biological pathway—could potentially complement single-SNP analysis and provide additional insights for the genetic architecture of complex diseases. Building upon existing P-value combining methods, we propose a class of highly flexible pathway analysis approaches based on an adaptive rank truncated product (ARTP) statistic that can effectively combine evidence of associations over different SNPs and genes within a pathway. The statistical significance of the pathway-level test-statistics is evaluated using a highly efficient permutation algorithm that remains computationally feasible irrespective of the size of the pathway and complexity of the underlying test-statistics for summarizing SNP- and gene-level associations. We demonstrate through simulation studies that a gene-based analysis, that treats the underlying genes, as opposed to the underlying SNPs, as the basic units for hypothesis testing, is a very robust and powerful approach to pathway-based association testing. We also illustrate the advantage of the proposed methods using a study of the association between the nicotinic receptor pathway and cigarette smoking behaviors.
Pathway analysis; genetic association study; permutation procedure
Clinical trial and epidemiological studies need high quality biospecimens from a representative sample of participants to investigate genetic influences on treatment response and disease. Obtaining blood biospecimens presents logistical and financial challenges. As a result, saliva biospecimen collection is becoming more frequent because of the ease of collection and lower cost. This article describes an assessment of saliva biospecimen samples collected through the mail, trial participant demographic and behavioral characteristics, and their association with saliva and DNA quantity and quality.
Saliva biospecimens were collected using the Oragene® DNA Self-Collection Kits from participants in a National Cancer Institute funded smoking cessation trial. Saliva biospecimens from 565 individuals were visually inspected for clarity prior to and after DNA extraction. DNA samples were then quantified by UV absorbance, PicoGreen®, and qPCR. Genotyping was performed on 11 SNPs using TaqMan® SNP assays and two VNTR assays. Univariate, correlation, and analysis of variance analyses were conducted to observe the relationship between saliva sample and participant characteristics.
The biospecimen kit return rate was 58.5% among those invited to participate (n = 967) and 47.1% among all possible COMPASS participants (n = 1202). Significant gender differences were observed with males providing larger saliva volume (4.7 vs. 4.5 ml, p = 0.019), samples that were more likely to be judged as cloudy (39.5% vs. 24.9%, p < 0.001), and samples with greater DNA yield as measured by UV (190.0 vs. 138.5, p = 0.002), but reduced % human DNA content (73.2 vs. 77.6 p = 0.005) than females. Other participant characteristics (age, self-identified ethnicity, baseline cigarettes per day) were associated with saliva clarity. Saliva volume and saliva and DNA clarity were positively correlated with total DNA yield by all three quantification measurements (all r > 0.21, P < 0.001), but negatively correlated with % human DNA content (saliva volume r = -0.148 and all P < 0.010). Genotyping completion rate was not influenced by saliva or DNA clarity.
Findings from this study show that demographic and behavioral characteristics of smoking cessation trial participants have significant associations with saliva and DNA metrics, but not with the performance of TaqMan® SNP or VNTR genotyping assays.
COMPASS; registered as NCT00301145 at clinicaltrials.gov.
Genetic variation may influence initial sensitivity to nicotine (i.e. during early tobacco exposure), perhaps helping to explain differential vulnerability to nicotine dependence. This study explored associations of functional candidate gene polymorphisms with initial sensitivity to nicotine in 101 young adult nonsmokers of European ancestry. Nicotine (0, 5, 10 μg/kg) was administered via nasal spray followed by mood, nicotine reward (e.g. “liking”) and perception (e.g. “feel effects”) measures, physiological responses, sensory processing (pre-pulse inhibition of startle), and performance tasks. Nicotine reinforcement was assessed in a separate session using a nicotine vs. placebo spray choice procedure. For the dopamine D4 receptor (DRD4 VNTR), presence of the 7 repeat allele was associated with greater aversive responses to nicotine (decreases in “vigor”, positive affect, and rapid information processing; increased cortisol) and reduced nicotine choice. Individuals with at least one DRD4 7-repeat allele also reported increased “feel effects” and greater startle response, but in men only. Also observed in men but not women were other genetic associations, such as greater “feel effects” and anger, and reduced fatigue, in the dopamine D2 receptor (DRD2 C957T SNP) TT versus CT or CC genotypes. Very few or no significant associations were seen for the DRD2/ANKK1 TaqIA polymorphism, the serotonin transporter promoter VNTR or 5HTTLPR (SLC6A4), the dopamine transporter 3’ VNTR (SLC6A3), and the mu opioid receptor A118G SNP (OPRM1). Although these results are preliminary, this study is the first to suggest that genetic polymorphisms related to function in the dopamine D4, and perhaps D2, receptor may modulate initial sensitivity to nicotine prior to the onset of dependence and may do so differentially between men and women.
nicotine; sensitivity; genetics; dopamine; reward; reinforcement
Negative mood increases smoking reinforcement and risk of relapse. We explored associations of gene variants in the dopamine, opioid, and serotonin pathways with smoking reward (“liking”) and reinforcement (latency to first puff, total puffs) as a function of negative mood and expected vs. actual nicotine content of the cigarette. Smokers of European ancestry (n=72) were randomized to one of four groups in a 2 × 2 balanced-placebo design, corresponding to manipulation of actual (0.6 mg vs. 0.05 mg) and expected (told nicotine, told denicotinized) nicotine “dose” in cigarettes during each of two sessions (negative vs. positive mood induction). Following mood induction and expectancy instructions, they sampled and rated the assigned cigarette, and then smoked additional cigarettes ad lib during continued mood induction. The increase in smoking amount due to negative mood was associated with: DRD2 C957T (CC>TT or CT), SLC6A3 (presence of 9 repeat > absence of 9), and among those given a nicotine cigarette, DRD4 (presence of 7 repeat > absence of 7) and DRD2/ANKK1 TaqIA (TT or CT > CC). SLC6A3 and DRD2/ANKK1 TaqIA were also associated with smoking reward and smoking latency. OPRM1 (AA > AG or GG) was associated with smoking reward, but SLC6A4 VNTR was unrelated to any of these measures. These results warrant replication but provide the first evidence for genetic associations with the acute increase in smoking reward and reinforcement due to negative mood.
smoking reward; reinforcement; mood; genetics; dopamine
To investigate the contribution of the dopamine transporter to dopaminergic reward-related behaviors and anthropometry, we evaluated associations between polymorphisms at the dopamine transporter gene(SLC6A3) and body mass index (BMI), among participants in the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial.
Four polymorphisms (rs6350, rs6413429, rs6347 and the 3' variable number of tandem repeat (3' VNTR) polymorphism) at the SLC6A3 gene were genotyped in 2,364 participants selected from the screening arm of PLCO randomly within strata of sex, age and smoking history. Height and weight at ages 20 and 50 years and baseline were assessed by questionnaire. BMI was calculated and categorized as underweight, normal, overweight and obese (<18.5, 18.5–24.9, 25.0–29.9, or ≥ 30 kg/m2, respectively). Odds ratios (ORs) and 95% confidence intervals (CIs) of SLC6A3 genotypes and haplotypes were computed using conditional logistic regression.
Compared with individuals having a normal BMI, obese individuals at the time of the baseline study questionnaire were less likely to possess the 3' VNTR variant allele with 9 copies of the repeated sequence in a dose-dependent model (** is referent; OR*9 = 0.80, OR99 = 0.47, ptrend = 0.005). Compared with individuals having a normal BMI at age 50, overweight individuals (A-C-G-* is referent; ORA-C-G-9 = 0.80, 95% CI 0.65–0.99, p = 0.04) and obese individuals (A-C-G-* is referent; ORA-C-G-9 = 0.70, 95% CI 0.49–0.99, p = 0.04) were less likely to possess the haplotype with the 3'variant allele (A-C-G-9).
Our results support a role of genetic variation at the dopamine transporter gene, SLC6A3, as a modifier of BMI.