Factors associated with genetic markers for rotator cuff disease in patients with atraumatic rotator cuff tears

For atraumatic rotator cuff tears, genetics contributes to symptomatic tear risk and may influence rotator cuff healing after surgical repair. But little is known about how genetic factors influence rotator cuff tear patient characteristics at presentation. We collected saliva samples for genotyping from atraumatic rotator cuff tear patients. We examined nine single nucleotide polymorphisms (SNPs) associated with cuff tears in prior literature. We estimated associations of SNP dosage with (1) age at tear diagnosis, (2) bilateral atraumatic tear prevalence, and (3) tear size. Linear regression was used to estimate associations with diagnosis age adjusted for sex and principal components. Logistic regression and ordinal logistic regression were used to estimate associations with bilateral tear prevalence and tear size category, respectively, adjusting for age, sex, and principal components. Of 344 eligible patients, 336 provided sufficient samples for genotyping. Median age at tear diagnosis was 61, 22% (N = 74) had bilateral atraumatic tears, and 9% (N = 29) had massive tears. SNP rs13107325 in the SLC39A8 gene and rs11850957 in the STXBP6 gene were associated with younger diagnosis age even after accounting for multiple comparisons (rs13107325: -4 years, 95% CI = -6.5, -1.4; rs11850957: -2.7 years, 95% CI = -4.3, -1.1). No other significant associations were observed with diagnosis age, tear size, or bilateral tear prevalence. SLC39A8 encodes a Mn transporter. STXBP6 may play a role in inflammatory responses by altering phagocytosis and antigen presentation of monocytes and macrophages. Further research is needed to determine if genetic markers can be used alongside patient characteristics to aid in identifying optimal surgical repair candidates.

Environmental and genetic risk factors associated with total knee arthroplasty following cruciate ligament surgery

Aims

The purpose of this study is to determine an individual's age-specific prevalence of total knee arthroplasty (TKA) after cruciate ligament surgery, and to identify clinical and genetic risk factors associated with undergoing TKA.

Methods

This study was a retrospective case-control study using the UK Biobank to identify individuals reporting a history of cruciate ligament surgery. Data from verbal history and procedural codes recorded through the NHS were used to identify instances of TKA. Patient clinical and genetic data were used to identify risk factors for progression from cruciate ligament surgery to TKA. Individuals without a history of cruciate ligament reconstruction were used for comparison.

Results

A total of 2,576 individuals with a history of cruciate ligament surgery were identified, with 290 (11.25%) undergoing TKA. In patients with prior cruciate ligament surgery, prevalence of TKA was 0.75% at age 45 years, 9.10% at age 65 years, and 20.43% at age 80 years. Patients with prior cruciate ligament surgery were 4.6 times more likely to have undergone TKA by age 55 years than individuals without prior cruciate ligament surgery. In the cruciate ligament surgery cohort, BMI > 30 kg/m2 (odds ratio (OR) 4.01 (95% confidence interval (CI) 2.74 to 5.87)), a job that always involved heavy manual or physical labour (OR 2.72 (95% CI 1.57 to 4.71)), or a job that always involved walking and standing (OR 2.58 (95% CI 1.58 to 4.20)) were associated with greater TKA odds. No single-nucleotide polymorphism (SNP) was associated with risk of TKA following cruciate ligament surgery.

Conclusion

Patients with a history of prior cruciate ligament surgery have substantially higher risk of TKA and undergo arthroplasty at a relatively younger age than individuals without a history of prior cruciate ligament surgery. Physically demanding work and obesity were associated with higher odds of TKA after cruciate ligament surgery, but no SNP was associated with risk of TKA.

Evaluation of Genetic and Nongenetic Risk Factors for Degenerative Cervical Myelopathy

STUDY DESIGN

Cohort study.

OBJECTIVE

We aimed to evaluate the associations of genetic and nongenetic factors with degenerative cervical myelopathy (DCM).

SUMMARY OF BACKGROUND DATA

There is mounting evidence for an inherited predisposition for DCM, but uncertainty remains regarding specific genetic markers involved. Similarly, nongenetic factors are thought to play a role.

MATERIALS AND METHODS

Using diagnosis codes from hospital records linked to the UK Biobank cohort, patients with cervical spondylosis were identified followed by the identification of a subset with DCM. Nongenetic variables evaluated included age, sex, race, Townsend deprivation index, body mass index, occupational demands, osteoporosis, and smoking. Genome-wide association analyses were conducted using logistic regression adjusted for age, sex, population principal components, and follow-up.

RESULTS

A total of 851 DCM cases out of 2787 cervical spondylosis patients were identified. Several nongenetic factors were independently associated with DCM including age [odds ratio (OR)=1.11, 95% CI=1.01-1.21, P =0.024], male sex (OR=1.63, 95% CI=1.37-1.93, P <0.001), and relative socioeconomic deprivation (OR=1.03, 95% CI=1.00-1.06, P =0.030). Asian race was associated with lower DCM risk (OR=0.44, 95% CI=0.22-0.85, P =0.014). We did not identify genome-wide significant (≤5×10 -8 ) single-nucleotide polymorphisms (SNPs) associated with DCM. The strongest genome-wide signals were at SNP rs67256809 in the intergenic region of the genes LINC02582 and FBXO15 on chromosome 18 ( P =1.12×10 -7 ) and rs577081672 in the GTPBP1 gene on chromosome 22 ( P =2.9×10 -7 ). No SNPs reported in prior DCM studies were significant after adjusting for replication attempts.

CONCLUSIONS

Increasing age, male sex, and relative socioeconomic deprivation were identified as independent risk factors for DCM, whereas Asian race was inversely associated. SNPs of potential interest were identified in GTPBP1 and an intergenic region on chromosome 18, but these associations did not reach genome-wide significance. Identification of genetic and nongenetic DCM susceptibility markers may guide understanding of DCM disease processes, inform risk, guide prevention and potentially inform surgical outcomes.

LEVEL OF EVIDENCE

Prognostic level III.

Occupational demands associated with rotator cuff disease surgery in the UK Biobank

OBJECTIVES

Physically-demanding occupations may increase rotator cuff disease (RCD) risk and need for surgery. We linked a job-exposure matrix (JEM) to the UK Biobank cohort study to measure physical occupational exposures and estimate associations with RCD surgery.

METHODS

Jobs and UK Standard Occupational Classification (SOC) codes were recorded during the UK Biobank verbal interview. Lifetime job histories were captured through a web-based survey. UK SOC codes were linked to a JEM based on the US O*NET database. O*NET-based scores [static strength, dynamic strength, general physical activities, handling/moving objects (range=1-7), time spent using hands, whole body vibration, and cramped/awkward positions (range=1-5)] were assigned to jobs. RCD surgeries were identified through linked national hospital inpatient records. Multivariable Cox regression was used to calculate hazard ratios (HR) as estimates of associations with RCD surgery. Among those with lifetime job histories, associations were estimated for duration of time with greatest exposure (top quartile of exposure).

RESULTS

Of 277 808 people reporting jobs, 1997 (0.7%) had an inpatient RCD surgery. After adjusting for age, sex, race, education, area deprivation, and body mass index, all O*NET variables considered were associated with RCD surgery (HR per point increase range=1.10-1.45, all P<0.005). A total of 100 929 people reported lifetime job histories, in which greater exposures were significantly associated with RCD surgery after >10 years of work (eg, HR for 11-20 versus 0 years with static strength score ≥4 = 2.06, 95% confidence interval 1.39-3.04).

CONCLUSIONS

Workplace physical demands are an important risk factor for RCD surgery, particularly for workers with more than a decade of exposure.

Genetic diversity fuels gene discovery for tobacco and alcohol use

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Functional characterization of enhancer activity during a long terminal repeat's evolution

Many transposable elements (TEs) contain transcription factor binding sites and are implicated as potential regulatory elements. However, TEs are rarely functionally tested for regulatory activity, which in turn limits our understanding of how TE regulatory activity has evolved. We systematically tested the human LTR18A subfamily for regulatory activity using massively parallel reporter assay (MPRA) and found AP-1- and CEBP-related binding motifs as drivers of enhancer activity. Functional analysis of evolutionarily reconstructed ancestral sequences revealed that LTR18A elements have generally lost regulatory activity over time through sequence changes, with the largest effects occurring owing to mutations in the AP-1 and CEBP motifs. We observed that the two motifs are conserved at higher rates than expected based on neutral evolution. Finally, we identified LTR18A elements as potential enhancers in the human genome, primarily in epithelial cells. Together, our results provide a model for the origin, evolution, and co-option of TE-derived regulatory elements.

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

BACKGROUND

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

METHODS

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

RESULTS

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

CONCLUSIONS

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

Studying the Utility of Using Genetics to Predict Smoking-Related Outcomes in a Population-Based Study and a Selected Cohort

INTRODUCTION

The purpose of this study is to examine the predictive utility of polygenic risk scores (PRSs) for smoking behaviors.

AIMS AND METHODS

Using summary statistics from the Sequencing Consortium of Alcohol and Nicotine use consortium, we generated PRSs of ever smoking, age of smoking initiation, cigarettes smoked per day, and smoking cessation for participants in the population-based Atherosclerosis Risk in Communities (ARIC) study (N = 8638), and the Collaborative Genetic Study of Nicotine Dependence (COGEND) (N = 1935). The outcomes were ever smoking, age of smoking initiation, heaviness of smoking, and smoking cessation.

RESULTS

In the European ancestry cohorts, each PRS was significantly associated with the corresponding smoking behavior outcome. In the ARIC cohort, the PRS z-score for ever smoking predicted smoking (odds ratio [OR]: 1.37; 95% confidence interval [CI]: 1.31, 1.43); the PRS z-score for age of smoking initiation was associated with age of smoking initiation (OR: 0.87; 95% CI: 0.82, 0.92); the PRS z-score for cigarettes per day was associated with heavier smoking (OR: 1.17; 95% CI: 1.11, 1.25); and the PRS z-score for smoking cessation predicted successful cessation (OR: 1.24; 95% CI: 1.17, 1.32). In the African ancestry cohort, the PRSs did not predict smoking behaviors.

CONCLUSIONS

Smoking-related PRSs were associated with smoking-related behaviors in European ancestry populations. This improvement in prediction is greatest in the lowest and highest genetic risk categories. The lack of prediction in African ancestry populations highlights the urgent need to increase diversity in research so that scientific advances can be applied to populations other than those of European ancestry.

IMPLICATIONS

This study shows that including both genetic ancestry and PRSs in a single model increases the ability to predict smoking behaviors compared with the model including only demographic characteristics. This finding is observed for every smoking-related outcome. Even though adding genetics is more predictive, the demographics alone confer substantial and meaningful predictive power. However, with increasing work in PRSs, the predictive ability will continue to improve.

Identification of a Novel Genetic Marker for Risk of Degenerative Rotator Cuff Disease Surgery in the UK Biobank

BACKGROUND

While evidence indicates that familial predisposition influences the risk of developing degenerative rotator cuff disease (RCD), knowledge of specific genetic markers is limited. We conducted a genome-wide association study of RCD surgery using the UK Biobank, a prospective cohort of 500,000 people (40 to 69 years of age at enrollment) with genotype data.

METHODS

Cases with surgery for degenerative RCD were identified using linked hospital records. The cases were defined as an International Classification of Diseases, Tenth Revision (ICD-10) code of M75.1 determined by a trauma/orthopaedic specialist and surgery consistent with RCD treatment. Cases were excluded if a diagnosis of traumatic injury had been made during the same hospital visit. For each case, up to 5 controls matched by age, sex, and follow-up time were chosen from the UK Biobank. Analyses were limited to European-ancestry individuals who were not third-degree or closer relations. We used logistic regression to test for genetic association of 674,405 typed and >10 million imputed markers, after adjusting for age, sex, population principal components, and follow-up.

RESULTS

We identified 2,917 RCD surgery cases and 14,158 matched controls. We observed 1 genome-wide significant signal (p < 5 × 10-8) for a novel locus tagged by rs2237352 in the CREB5 gene on chromosome 7 (odds ratio [OR] = 1.17, 95% confidence interval [CI] = 1.11 to 1.24). The single-nucleotide polymorphism (SNP) rs2237352 was imputed with a high degree of confidence (info score = 0.9847) and is common, with a minor allele frequency of 47%. After expanding the control sample to include additional unmatched non-cases, rs2237352 and another SNP in the CREB5 gene, rs12700903, were genome-wide significant. We did not detect genome-wide significant signals at loci associated with RCD in previous studies.

CONCLUSIONS

We identified a novel association between a variant in the CREB5 gene and RCD surgery. Validation of this finding in studies with imaging data to confirm diagnoses will be an important next step.

CLINICAL RELEVANCE

Identification of genetic RCD susceptibility markers can guide understanding of biological processes in rotator cuff degeneration and help inform disease risk in the clinical setting.

LEVEL OF EVIDENCE

Prognostic Level III. See Instructions for Authors for a complete description of levels of evidence.

Current Smoking Raises Risk of Incident Hypertension: Hispanic Community Health Study-Study of Latinos

BACKGROUND

Hypertension has been implicated as a smoking-related risk factor for cardiovascular disease but the dose-response relationship is incompletely described. Hispanics, who often have relatively light smoking exposures, have been understudied in this regard.

METHODS

We used data from a 6-year follow-up study of US Hispanic adults aged 18-76 to address the dose-response linking cigarette use with incident hypertension, which was defined by measured blood pressure above 140/90 mm Hg or initiation of antihypertensive medications. Adjustment was performed for potential confounders and mediators, including urinary albumin-to-creatinine ratio which worsened over time among smokers.

RESULTS

Current smoking was associated with incident hypertension, with a threshold effect above 5 cumulative pack-years of smoking (vs. never smokers, hazard ratio for hypertension [95% confidence interval] of 0.95 [0.67, 1.35] for 0-5 pack-years, 1.47 [1.05, 2.06] for 5-10 pack-years, 1.40 [1.00, 1.96] for 10-20 pack-years, and 1.34 [1.09, 1.66] for ≥20 pack-years, P = 0.037). In contrast to current smokers, former smokers did not appear to have increased risk of hypertension, even at the highest cumulative pack-years of past exposure.

CONCLUSIONS

The results confirm that smoking constitutes a hypertension risk factor in Hispanic adults. A relatively modest cumulative dose of smoking, above 5 pack-years of exposure, raises risk of hypertension by over 30%. The increased hypertension risk was confined to current smokers, and did not increase further with higher pack-year levels. The lack of a smoking-hypertension association in former smokers underscores the value of smoking cessation.

A large-scale genome-wide association study meta-analysis of cannabis use disorder

BACKGROUND

Variation in liability to cannabis use disorder has a strong genetic component (estimated twin and family heritability about 50-70%) and is associated with negative outcomes, including increased risk of psychopathology. The aim of the study was to conduct a large genome-wide association study (GWAS) to identify novel genetic variants associated with cannabis use disorder.

METHODS

To conduct this GWAS meta-analysis of cannabis use disorder and identify associations with genetic loci, we used samples from the Psychiatric Genomics Consortium Substance Use Disorders working group, iPSYCH, and deCODE (20 916 case samples, 363 116 control samples in total), contrasting cannabis use disorder cases with controls. To examine the genetic overlap between cannabis use disorder and 22 traits of interest (chosen because of previously published phenotypic correlations [eg, psychiatric disorders] or hypothesised associations [eg, chronotype] with cannabis use disorder), we used linkage disequilibrium score regression to calculate genetic correlations.

FINDINGS

We identified two genome-wide significant loci: a novel chromosome 7 locus (FOXP2, lead single-nucleotide polymorphism [SNP] rs7783012; odds ratio [OR] 1·11, 95% CI 1·07-1·15, p=1·84 × 10-9) and the previously identified chromosome 8 locus (near CHRNA2 and EPHX2, lead SNP rs4732724; OR 0·89, 95% CI 0·86-0·93, p=6·46 × 10-9). Cannabis use disorder and cannabis use were genetically correlated (rg 0·50, p=1·50 × 10-21), but they showed significantly different genetic correlations with 12 of the 22 traits we tested, suggesting at least partially different genetic underpinnings of cannabis use and cannabis use disorder. Cannabis use disorder was positively genetically correlated with other psychopathology, including ADHD, major depression, and schizophrenia.

INTERPRETATION

These findings support the theory that cannabis use disorder has shared genetic liability with other psychopathology, and there is a distinction between genetic liability to cannabis use and cannabis use disorder.

FUNDING

National Institute of Mental Health; National Institute on Alcohol Abuse and Alcoholism; National Institute on Drug Abuse; Center for Genomics and Personalized Medicine and the Centre for Integrative Sequencing; The European Commission, Horizon 2020; National Institute of Child Health and Human Development; Health Research Council of New Zealand; National Institute on Aging; Wellcome Trust Case Control Consortium; UK Research and Innovation Medical Research Council (UKRI MRC); The Brain & Behavior Research Foundation; National Institute on Deafness and Other Communication Disorders; Substance Abuse and Mental Health Services Administration (SAMHSA); National Institute of Biomedical Imaging and Bioengineering; National Health and Medical Research Council (NHMRC) Australia; Tobacco-Related Disease Research Program of the University of California; Families for Borderline Personality Disorder Research (Beth and Rob Elliott) 2018 NARSAD Young Investigator Grant; The National Child Health Research Foundation (Cure Kids); The Canterbury Medical Research Foundation; The New Zealand Lottery Grants Board; The University of Otago; The Carney Centre for Pharmacogenomics; The James Hume Bequest Fund; National Institutes of Health: Genes, Environment and Health Initiative; National Institutes of Health; National Cancer Institute; The William T Grant Foundation; Australian Research Council; The Virginia Tobacco Settlement Foundation; The VISN 1 and VISN 4 Mental Illness Research, Education, and Clinical Centers of the US Department of Veterans Affairs; The 5th Framework Programme (FP-5) GenomEUtwin Project; The Lundbeck Foundation; NIH-funded Shared Instrumentation Grant S10RR025141; Clinical Translational Sciences Award grants; National Institute of Neurological Disorders and Stroke; National Heart, Lung, and Blood Institute; National Institute of General Medical Sciences.

Review and Consensus on Pharmacogenomic Testing in Psychiatry

The implementation of pharmacogenomic (PGx) testing in psychiatry remains modest, in part due to divergent perceptions of the quality and completeness of the evidence base and diverse perspectives on the clinical utility of PGx testing among psychiatrists and other healthcare providers. Recognizing the current lack of consensus within the field, the International Society of Psychiatric Genetics assembled a group of experts to conduct a narrative synthesis of the PGx literature, prescribing guidelines, and product labels related to psychotropic medications as well as the key considerations and limitations related to the use of PGx testing in psychiatry. The group concluded that to inform medication selection and dosing of several commonly-used antidepressant and antipsychotic medications, current published evidence, prescribing guidelines, and product labels support the use of PGx testing for 2 cytochrome P450 genes (CYP2D6, CYP2C19). In addition, the evidence supports testing for human leukocyte antigen genes when using the mood stabilizers carbamazepine (HLA-A and HLA-B), oxcarbazepine (HLA-B), and phenytoin (CYP2C9, HLA-B). For valproate, screening for variants in certain genes (POLG, OTC, CSP1) is recommended when a mitochondrial disorder or a urea cycle disorder is suspected. Although barriers to implementing PGx testing remain to be fully resolved, the current trajectory of discovery and innovation in the field suggests these barriers will be overcome and testing will become an important tool in psychiatry.

Expanding the genetic architecture of nicotine dependence and its shared genetics with multiple traits

Cigarette smoking is the leading cause of preventable morbidity and mortality. Genetic variation contributes to initiation, regular smoking, nicotine dependence, and cessation. We present a Fagerström Test for Nicotine Dependence (FTND)-based genome-wide association study in 58,000 European or African ancestry smokers. We observe five genome-wide significant loci, including previously unreported loci MAGI2/GNAI1 (rs2714700) and TENM2 (rs1862416), and extend loci reported for other smoking traits to nicotine dependence. Using the heaviness of smoking index from UK Biobank (N = 33,791), rs2714700 is consistently associated; rs1862416 is not associated, likely reflecting nicotine dependence features not captured by the heaviness of smoking index. Both variants influence nearby gene expression (rs2714700/MAGI2-AS3 in hippocampus; rs1862416/TENM2 in lung), and expression of genes spanning nicotine dependence-associated variants is enriched in cerebellum. Nicotine dependence (SNP-based heritability = 8.6%) is genetically correlated with 18 other smoking traits (rg = 0.40-1.09) and co-morbidities. Our results highlight nicotine dependence-specific loci, emphasizing the FTND as a composite phenotype that expands genetic knowledge of smoking.

Dissecting the genetic overlap of smoking behaviors, lung cancer, and chronic obstructive pulmonary disease: A focus on nicotinic receptors and nicotine metabolizing enzyme

Smoking is a major contributor to lung cancer and chronic obstructive pulmonary disease (COPD). Two of the strongest genetic associations of smoking-related phenotypes are the chromosomal regions 15q25.1, encompassing the nicotinic acetylcholine receptor subunit genes CHRNA5-CHRNA3-CHRNB4, and 19q13.2, encompassing the nicotine metabolizing gene CYP2A6. In this study, we examined genetic relations between cigarettes smoked per day, smoking cessation, lung cancer, and COPD. Data consisted of genome-wide association study summary results. Genetic correlations were estimated using linkage disequilibrium score regression software. For each pair of outcomes, z-score-z-score (ZZ) plots were generated. Overall, heavier smoking and decreased smoking cessation showed positive genetic associations with increased lung cancer and COPD risk. The chromosomal region 19q13.2, however, showed a different correlational pattern. For example, the effect allele-C of the sentinel SNP (rs56113850) within CYP2A6 was associated with an increased risk of heavier smoking (z-score = 19.2; p = 1.10 × 10-81 ), lung cancer (z-score = 8.91; p = 5.02 × 10-19 ), and COPD (z-score = 4.04; p = 5.40 × 10-5 ). Surprisingly, this allele-C (rs56113850) was associated with increased smoking cessation (z-score = -8.17; p = 2.52 × 10-26 ). This inverse relationship highlights the need for additional investigation to determine how CYP2A6 variation could increase smoking cessation while also increasing the risk of lung cancer and COPD likely through increased cigarettes smoked per day.

Variants in the CHRNA5-CHRNA3-CHRNB4 Region of Chromosome 15 Predict Gastrointestinal Adverse Events in the Transdisciplinary Tobacco Use Research Center Smoking Cessation Trial

INTRODUCTION

Reducing adverse events from pharmacologic treatment is an important goal of precision medicine and identifying genetic predictors of adverse events is a step toward this goal. In 2012, King et al. reported associations between genetic variants and adverse events in a placebo-controlled smoking cessation trial of varenicline and bupropion. Strong associations were found between gastrointestinal adverse events and 11 variants in the CHRNA5-CHRNA3-CHRNB4 region of chromosome 15, a region repeatedly associated with smoking-related phenotypes. Our goal was to replicate, in an independent sample, the impact of variants in the CHRNA5-CHRNA3-CHRNB4 region on gastrointestinal adverse events and to extend the analyses to adherence and smoking cessation.

METHODS

The University of Wisconsin Transdisciplinary Tobacco Use Research Center (TTURC) conducted a multiarmed, placebo-controlled smoking cessation trial of bupropion and nicotine replacement therapy that included 985 genotyped European-ancestry participants. We evaluated relationships between our key variables using logistic regression.

RESULTS

Gastrointestinal adverse events were experienced by 31.6% TTURC participants. Each of the CHRNA5-CHRNA3-CHRNB4 associations from the King et al. study was found in TTURC, with the same direction of effect. Neither these variants nor the gastrointestinal adverse events themselves were associated with adherence to medication or successful smoking cessation.

CONCLUSIONS

Variants in the CHRNA5-CHRNA3-CHRNB4 region of chromosome 15 are associated with gastrointestinal adverse events in smoking cessation. Additional independent variants in this region strengthen the association. The consistency between the results of these two independent studies supports the conclusion that these findings reflect biological response to the use of smoking cessation medication.

IMPLICATIONS

The fact that our findings from the TTURC smoking cessation trial support the independent findings of King et al. suggest that associations of variants in the CHRNA5-CHRNA3-CHRNB4 region of chromosome 15 with gastrointestinal adverse events while taking medications for smoking cessation reflect biology. However, although adherence to medication was a strong predictor of successful smoking cessation in TTURC, neither adverse events nor the genetic variants associated with them predicted either adherence or successful cessation in this study. Thus, although we should strive to minimize adverse events during treatment, we should not expect that to increase successful smoking cessation substantially.

Genetic Variant in CHRNA5 and Response to Varenicline and Combination Nicotine Replacement in a Randomized Placebo-Controlled Trial

It is unclear if genetic variants affect smoking cessation treatment response. This study tested whether variants in the cholinergic receptor nicotinic alpha 5 subunit (CHRNA5) predict response to smoking cessation medication by directly comparing the two most effective smoking cessation pharmacotherapies. In this genotype-stratified randomized, double-blind, placebo-controlled clinical trial (May 2015-August 2019 in St Louis, Missouri), smokers were randomized by genotype in blocks of six (1:1:1 ratio) to three conditions: 12 weeks of placebo (n = 273), combination nicotine patch and lozenge (combination nicotine replacement therapy, cNRT, n = 275), or varenicline (n = 274). All participants received counseling and were followed for 12 months. The primary end point was biochemically verified 7-day point prevalence abstinence at the end of treatment (EOT, week 12). Trial registration and eligibility criteria are on clinicaltrials.gov (https://clinicaltrials.gov/) (NCT02351167). We conducted the genetic analyses separately for 516 European ancestry (EA) smokers and 306 non-EA smokers (including 270 African American smokers). In African American smokers, there was a genotype-by-treatment interaction for EOT abstinence (χ²  = 10.7, degrees of freedom = 2. P = 0.0049): specifically, cNRT was more effective in smokers with rs16969968 GG genotype than was placebo, while varenicline was more effective in smokers of GA/AA genotypes. In EA ancestry smokers, there was no significant genotype-by-treatment interaction. In the whole sample, although both were effective at EOT, only varenicline, and not cNRT, was significantly effective relative to placebo at 6-month follow-up. Importantly, this study suggests that genetic information can further enhance smoking cessation treatment effectiveness.

Risk factors for surgery due to rotator cuff disease in a population-based cohort

AIMS

Few risk factors for rotator cuff disease (RCD) and corresponding treatment have been firmly established. The aim of this study was to evaluate the relationship between numerous risk factors and the incidence of surgery for RCD in a large cohort.

METHODS

A population-based cohort of people aged between 40 and 69 years in the UK (the UK Biobank) was studied. People who underwent surgery for RCD were identified through a link with NHS inpatient records covering a mean of eight years after enrolment. Multivariate Cox proportional hazards regression was used to calculate hazard ratios (HRs) as estimates of associations with surgery for RCD accounting for confounders. The risk factors which were considered included age, sex, race, education, Townsend deprivation index, body mass index (BMI), occupational demands, and exposure to smoking.

RESULTS

Of the 421,894 people who were included, 47% were male. The mean age at the time of enrolment was 56 years (40 to 69). A total of 2,156 people were identified who underwent surgery for RCD. Each decade increase in age was associated with a 55% increase in the incidence of RCD surgery (95% confidence interval (CI) 46% to 64%). Male sex, non-white race, lower deprivation score, and higher BMI were significantly associated with a higher risk of surgery for RCD (all p < 0.050). Greater occupational physical demands were significantly associated with higher rates of RCD surgery (HR = 2.1, 1.8, and 1.4 for 'always', 'usually', and 'sometimes' doing heavy manual labour vs 'never', all p < 0.001). Former smokers had significantly higher rates of RCD surgery than those who had never smoked (HR 1.23 (95% CI 1.12 to 1.35), p < 0.001), while current smokers had similar rates to those who had never smoked (HR 0.94 (95% CI 0.80 to 1.11)). Among those who had never smoked, the risk of surgery was higher among those with more than one household member who smoked (HR 1.78 (95% CI 1.08 to 2.92)). The risk of RCD surgery was not significantly related to other measurements of secondhand smoking.

CONCLUSION

Many factors were independently associated with surgery for RCD, including older age, male sex, higher BMI, lower deprivation score, and higher occupational physical demands. Several of the risk factors which were identified are modifiable, suggesting that the healthcare burden of RCD might be reduced through the pursuit of public health goals, such as reducing obesity and modifying occupational demands. Cite this article: Bone Joint J 2020;102-B(3):352-359.

International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci

The risk of posttraumatic stress disorder (PTSD) following trauma is heritable, but robust common variants have yet to be identified. In a multi-ethnic cohort including over 30,000 PTSD cases and 170,000 controls we conduct a genome-wide association study of PTSD. We demonstrate SNP-based heritability estimates of 5-20%, varying by sex. Three genome-wide significant loci are identified, 2 in European and 1 in African-ancestry analyses. Analyses stratified by sex implicate 3 additional loci in men. Along with other novel genes and non-coding RNAs, a Parkinson's disease gene involved in dopamine regulation, PARK2, is associated with PTSD. Finally, we demonstrate that polygenic risk for PTSD is significantly predictive of re-experiencing symptoms in the Million Veteran Program dataset, although specific loci did not replicate. These results demonstrate the role of genetic variation in the biology of risk for PTSD and highlight the necessity of conducting sex-stratified analyses and expanding GWAS beyond European ancestry populations.

E-cigarette Usage Is Associated With Increased Past-12-Month Quit Attempts and Successful Smoking Cessation in Two US Population-Based Surveys

INTRODUCTION

We examined past-12-month quit attempts and smoking cessation from 2006 to 2016 while accounting for demographic shifts in the US population. In addition, we sought to understand whether the current use of electronic cigarettes was associated with a change in past-12-month quit attempts and successful smoking cessation at the population level.

METHODS

We analyzed data from 25- to 44-year-olds from the National Health Interview Survey (NHIS) from 2006 to 2016 (N = 26,354) and the Tobacco Use Supplement to the Current Population Survey (TUS-CPS) in 2006-2007, 2010-2011, and 2014-2015 (N = 33,627). Data on e-cigarette use were available in the 2014-2016 NHIS and 2014-2015 TUS-CPS surveys.

RESULTS

Past-12-month quit attempts and smoking cessation increased in recent years compared with 2006. Current e-cigarette use was associated with higher quit attempts (adjusted odds ratio [aOR] = 2.29, 95% confidence interval [CI] = 1.87 to 2.81, p < .001) and greater smoking cessation (aOR = 1.64, 95% CI = 1.21 to 2.21, p = .001) in the NHIS. Multivariable logistic regression of the TUS-CPS data showed that current e-cigarette use was similarly significantly associated with increased past-12-month quit attempts and smoking cessation. Significant interactions were found for smoking frequency (everyday and some-day smoking) and current e-cigarette use for both outcomes (p < .0001) with the strongest positive effects seen in everyday smokers.

CONCLUSIONS

Compared with 2006, past-12-month quit attempts and smoking cessation increased among adults aged 25-44 in recent years. Current e-cigarette use was associated with increased past-12-month quit attempts and successful smoking cessation among established smokers. These findings are relevant to future tobacco policy decisions.

IMPLICATIONS

E-cigarettes were introduced into the US market over the past decade. During this period, past-12-month quit attempts and smoking cessation have increased among US adults aged 25-44. These trends are inconsistent with the hypothesis that e-cigarette use is delaying quit attempts and leading to decreased smoking cessation. In contrast, current e-cigarette use was associated with significantly higher past-12-month quit attempts and past-12-month cessation. These findings suggest that e-cigarette use contributes to a reduction in combustible cigarette use among established smokers.

Genome-Wide Association Study of Heavy Smoking and Daily/Nondaily Smoking in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL)

Introduction

Genetic variants associated with nicotine dependence have previously been identified, primarily in European-ancestry populations. No genome-wide association studies (GWAS) have been reported for smoking behaviors in Hispanics/Latinos in the United States and Latin America, who are of mixed ancestry with European, African, and American Indigenous components.

Methods

We examined genetic associations with smoking behaviors in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) (N = 12 741 with smoking data, 5119 ever-smokers), using ~2.3 million genotyped variants imputed to the 1000 Genomes Project phase 3. Mixed logistic regression models accounted for population structure, sampling, relatedness, sex, and age.

Results

The known region of CHRNA5, which encodes the α5 cholinergic nicotinic receptor subunit, was associated with heavy smoking at genome-wide significance (p ≤ 5 × 10-8) in a comparison of 1929 ever-smokers reporting cigarettes per day (CPD) > 10 versus 3156 reporting CPD ≤ 10. The functional variant rs16969968 in CHRNA5 had a p value of 2.20 × 10-7 and odds ratio (OR) of 1.32 for the minor allele (A); its minor allele frequency was 0.22 overall and similar across Hispanic/Latino background groups (Central American = 0.17; South American = 0.19; Mexican = 0.18; Puerto Rican = 0.22; Cuban = 0.29; Dominican = 0.19). CHRNA4 on chromosome 20 attained p < 10-4, supporting prior findings in non-Hispanics. For nondaily smoking, which is prevalent in Hispanic/Latino smokers, compared to daily smoking, loci on chromosomes 2 and 4 achieved genome-wide significance; replication attempts were limited by small Hispanic/Latino sample sizes.

Conclusions

Associations of nicotinic receptor gene variants with smoking, first reported in non-Hispanic European-ancestry populations, generalized to Hispanics/Latinos despite different patterns of smoking behavior.

Implications

We conducted the first large-scale genome-wide association study (GWAS) of smoking behavior in a US Hispanic/Latino cohort, and the first GWAS of daily/nondaily smoking in any population. Results show that the region of the nicotinic receptor subunit gene CHRNA5, which in non-Hispanic European-ancestry smokers has been associated with heavy smoking as well as cessation and treatment efficacy, is also significantly associated with heavy smoking in this Hispanic/Latino cohort. The results are an important addition to understanding the impact of genetic variants in understudied Hispanic/Latino smokers.

Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use

Tobacco and alcohol use are leading causes of mortality that influence risk for many complex diseases and disorders1. They are heritable2,3 and etiologically related4,5 behaviors that have been resistant to gene discovery efforts6-11. In sample sizes up to 1.2 million individuals, we discovered 566 genetic variants in 406 loci associated with multiple stages of tobacco use (initiation, cessation, and heaviness) as well as alcohol use, with 150 loci evidencing pleiotropic association. Smoking phenotypes were positively genetically correlated with many health conditions, whereas alcohol use was negatively correlated with these conditions, such that increased genetic risk for alcohol use is associated with lower disease risk. We report evidence for the involvement of many systems in tobacco and alcohol use, including genes involved in nicotinic, dopaminergic, and glutamatergic neurotransmission. The results provide a solid starting point to evaluate the effects of these loci in model organisms and more precise substance use measures.

Transancestral GWAS of alcohol dependence reveals common genetic underpinnings with psychiatric disorders

Liability to alcohol dependence (AD) is heritable, but little is known about its complex polygenic architecture or its genetic relationship with other disorders. To discover loci associated with AD and characterize the relationship between AD and other psychiatric and behavioral outcomes, we carried out the largest genome-wide association study to date of DSM-IV-diagnosed AD. Genome-wide data on 14,904 individuals with AD and 37,944 controls from 28 case-control and family-based studies were meta-analyzed, stratified by genetic ancestry (European, n = 46,568; African, n = 6,280). Independent, genome-wide significant effects of different ADH1B variants were identified in European (rs1229984; P = 9.8 × 10-13) and African ancestries (rs2066702; P = 2.2 × 10-9). Significant genetic correlations were observed with 17 phenotypes, including schizophrenia, attention deficit-hyperactivity disorder, depression, and use of cigarettes and cannabis. The genetic underpinnings of AD only partially overlap with those for alcohol consumption, underscoring the genetic distinction between pathological and nonpathological drinking behaviors.

Use of polygenic risk scores of nicotine metabolism in predicting smoking behaviors

AIM

This study tests whether polygenic risk scores (PRSs) for nicotine metabolism predict smoking behaviors in independent data.

MATERIALS & METHODS

Linear regression, logistic regression and survival analyses were used to analyze nicotine metabolism PRSs and nicotine metabolism, smoking quantity and smoking cessation.

RESULTS

Nicotine metabolism PRSs based on two genome wide association studies (GWAS) meta-analyses significantly predicted nicotine metabolism biomarkers (R² range: 9.2-16%; minimum p = 7.6 × 10^-8). The GWAS top hit variant rs56113850 significantly predicted nicotine metabolism biomarkers (R² range: 14-17%; minimum p = 4.4 × 10^-8). There was insufficient evidence for these PRSs predicting smoking quantity and smoking cessation.

CONCLUSION

Results suggest that nicotine metabolism PRSs based on GWAS meta-analyses predict an individual's nicotine metabolism, so does use of the top hit variant. We anticipate that PRSs will enter clinical medicine, but additional research is needed to develop a more comprehensive genetic score to predict smoking behaviors.

Genome-wide association study across European and African American ancestries identifies a SNP in DNMT3B contributing to nicotine dependence

Cigarette smoking is a leading cause of preventable mortality worldwide. Nicotine dependence, which reduces the likelihood of quitting smoking, is a heritable trait with firmly established associations with sequence variants in nicotine acetylcholine receptor genes and at other loci. To search for additional loci, we conducted a genome-wide association study (GWAS) meta-analysis of nicotine dependence, totaling 38,602 smokers (28,677 Europeans/European Americans and 9925 African Americans) across 15 studies. In this largest-ever GWAS meta-analysis for nicotine dependence and the largest-ever cross-ancestry GWAS meta-analysis for any smoking phenotype, we reconfirmed the well-known CHRNA5-CHRNA3-CHRNB4 genes and further yielded a novel association in the DNA methyltransferase gene DNMT3B. The intronic DNMT3B rs910083-C allele (frequency=44-77%) was associated with increased risk of nicotine dependence at P=3.7 × 10-8 (odds ratio (OR)=1.06 and 95% confidence interval (CI)=1.04-1.07 for severe vs mild dependence). The association was independently confirmed in the UK Biobank (N=48,931) using heavy vs never smoking as a proxy phenotype (P=3.6 × 10-4, OR=1.05, and 95% CI=1.02-1.08). Rs910083-C is also associated with increased risk of squamous cell lung carcinoma in the International Lung Cancer Consortium (N=60,586, meta-analysis P=0.0095, OR=1.05, and 95% CI=1.01-1.09). Moreover, rs910083-C was implicated as a cis-methylation quantitative trait locus (QTL) variant associated with higher DNMT3B methylation in fetal brain (N=166, P=2.3 × 10-26) and a cis-expression QTL variant associated with higher DNMT3B expression in adult cerebellum from the Genotype-Tissue Expression project (N=103, P=3.0 × 10-6) and the independent Brain eQTL Almanac (N=134, P=0.028). This novel DNMT3B cis-acting QTL variant highlights the importance of genetically influenced regulation in brain on the risks of nicotine dependence, heavy smoking and consequent lung cancer.

Is the Fagerström test for nicotine dependence invariant across secular trends in smoking? A question for cross-birth cohort analysis of nicotine dependence

BACKGROUND

The Fagerström Test for Nicotine Dependence (FTND), a derivation of the Fagerström Tolerance Questionnaire, was first published in 1991. The FTND remains one of the most widely used measures of nicotine dependence for studying genetic and epidemiological risk factors and the likelihood of smoking cessation. However, it is unclear whether secular trends in patterns of smoking alter the psychometric properties of the FTND and its interpretation.

METHODS

We examined measurement invariance in the lifetime and current FTND scores across birth cohorts using participants drawn from six study samples (N = 13,775).

RESULTS

We found significant (p < 0.05) measurement non-invariance in means and factor loadings of most FTND items by birth cohort, but effect sizes, ranging from r2 = 0.0001 to r2 = 0.0035, indicated that less than 0.5% of the model variance was explained by the measurement non-invariance for each factor loading. To assess its impact, we regressed the lifetime FTND latent variable on well-established factors associated with nicotine dependence (quitting smoking and the nicotinic acetylcholine receptor gene [CHRNA5] variant rs16969968, separately), and we observed that the regression coefficients were unchanged between models with and without adjustment for measurement non-invariance.

CONCLUSIONS

These findings suggest that possible FTND non-invariance that occurs across study samples of various birth years has a negligible impact on study results.

The Value of Biosamples in Smoking Cessation Trials: A Review of Genetic, Metabolomic, and Epigenetic Findings

Introduction

Human genetic research has succeeded in definitively identifying multiple genetic variants associated with risk for nicotine dependence and heavy smoking. To build on these advances, and to aid in reducing the prevalence of smoking and its consequent health harms, the next frontier is to identify genetic predictors of successful smoking cessation and also of the efficacy of smoking cessation treatments ("pharmacogenomics"). More broadly, additional biomarkers that can be quantified from biosamples also promise to aid "Precision Medicine" and the personalization of treatment, both pharmacological and behavioral.

Aims and Methods

To motivate ongoing and future efforts, here we review several compelling genetic and biomarker findings related to smoking cessation and treatment.

Results

These Key results involve genetic variants in the nicotinic receptor subunit gene CHRNA5, variants in the nicotine metabolism gene CYP2A6, and the nicotine metabolite ratio. We also summarize reports of epigenetic changes related to smoking behavior.

Conclusions

The results to date demonstrate the value and utility of data generated from biosamples in clinical treatment trial settings. This article cross-references a companion paper in this issue that provides practical guidance on how to incorporate biosample collection into a planned clinical trial and discusses avenues for harmonizing data and fostering consortium-based, collaborative research on the pharmacogenomics of smoking cessation.

Implications

Evidence is emerging that certain genotypes and biomarkers are associated with smoking cessation success and efficacy of smoking cessation treatments. We review key findings that open potential avenues for personalizing smoking cessation treatment according to an individual's genetic or metabolic profile. These results provide important incentive for smoking cessation researchers to collect biosamples and perform genotyping in research studies and clinical trials.

Leveraging Genomic Data in Smoking Cessation Trials in the Era of Precision Medicine: Why and How

Implications

This article outlines a framework for the consistent integration of biological data/samples into smoking cessation pharmacotherapy trials, aligned with the objectives of the recently unveiled Precision Medicine Initiative. Our goal is to encourage and provide support for treatment researchers to consider biosample collection and genotyping their existing samples as well as integrating genetic analyses into their study design in order to realize precision medicine in treatment of nicotine dependence.

Largest GWAS of PTSD (N=20 070) yields genetic overlap with schizophrenia and sex differences in heritability

The Psychiatric Genomics Consortium-Posttraumatic Stress Disorder group (PGC-PTSD) combined genome-wide case-control molecular genetic data across 11 multiethnic studies to quantify PTSD heritability, to examine potential shared genetic risk with schizophrenia, bipolar disorder, and major depressive disorder and to identify risk loci for PTSD. Examining 20 730 individuals, we report a molecular genetics-based heritability estimate (h2SNP) for European-American females of 29% that is similar to h2SNP for schizophrenia and is substantially higher than h2SNP in European-American males (estimate not distinguishable from zero). We found strong evidence of overlapping genetic risk between PTSD and schizophrenia along with more modest evidence of overlap with bipolar and major depressive disorder. No single-nucleotide polymorphisms (SNPs) exceeded genome-wide significance in the transethnic (overall) meta-analysis and we do not replicate previously reported associations. Still, SNP-level summary statistics made available here afford the best-available molecular genetic index of PTSD-for both European- and African-American individuals-and can be used in polygenic risk prediction and genetic correlation studies of diverse phenotypes. Publication of summary statistics for ∼10 000 African Americans contributes to the broader goal of increased ancestral diversity in genomic data resources. In sum, the results demonstrate genetic influences on the development of PTSD, identify shared genetic risk between PTSD and other psychiatric disorders and highlight the importance of multiethnic/racial samples. As has been the case with schizophrenia and other complex genetic disorders, larger sample sizes are needed to identify specific risk loci.

Toward the implementation of genomic applications for smoking cessation and smoking-related diseases

The incorporation of genomic information into routine care settings is a burgeoning area for investigation in behavioral medicine. The past decade has witnessed rapid advancements in knowledge of genetic biomarkers associated with smoking behaviors and tobacco-related morbidity and mortality, providing the basis for promising genomic applications in clinical and community settings. We assessed the current state of readiness for implementing genomic applications involving variation in the α5 nicotinic cholinergic receptor subunit gene CHRNA5 and smoking outcomes (behaviors and related diseases) using a process that could be translatable to a wide range of genomic applications in behavioral medicine. We reviewed the scientific literature involving CHRNA5 genetic variation and smoking cessation, and then summarized and synthesized a chain of evidence according to analytic validity, clinical validity, clinical utility, and ethical, legal, and social implications (ACCE), a well-established set of criteria used to evaluate genomic applications. Our review identified at least three specific genomic applications for which implementation may be considered, including the use of CHRNA5 genetic test results for informing disease risk, optimizing smoking cessation treatment, and motivating smoking behavior change. For these genomic applications, we rated analytic validity as convincing, clinical validity as adequate, and clinical utility and ethical, legal, and social implications as inadequate. For clinical genomic applications involving CHRNA5 variation and smoking outcomes, research efforts now need to focus on establishing clinical utility. This approach is compatible with pre-implementation research, which is also needed to accelerate translation, improve innovation design, and understand and refine system processes involved in implementation. This study informs the readiness to incorporate smoking-related genomic applications in real-world settings and facilitates cross-disciplinary collaboration to accelerate the integration of evidence-based genomics in behavioral medicine.

CYP2A6 metabolism in the development of smoking behaviors in young adults

Cytochrome P450 2A6 (CYP2A6) encodes the enzyme responsible for the majority of nicotine metabolism. Previous studies support that slow metabolizers smoke fewer cigarettes once nicotine dependent but provide conflicting results on the role of CYP2A6 in the development of dependence. By focusing on the critical period of young adulthood, this study examines the relationship of CYP2A6 variation and smoking milestones. A total of 1209 European American young adults enrolled in the Collaborative Study on the Genetics of Alcoholism were genotyped for CYP2A6 variants to calculate a previously well-validated metric that estimates nicotine metabolism. This metric was not associated with the transition from never smoking to smoking initiation nor with the transition from initiation to daily smoking (P > 0.4). But among young adults who had become daily smokers (n = 506), decreased metabolism was associated with increased risk of nicotine dependence (P = 0.03) (defined as Fagerström Test for Nicotine Dependence score ≥4). This finding was replicated in the Collaborative Genetic Study of Nicotine Dependence with 335 young adult daily smokers (P = 0.02). Secondary meta-analysis indicated that slow metabolizers had a 53 percent increased odds (OR = 1.53, 95 percent CI 1.11-2.11, P = 0.009) of developing nicotine dependence compared with normal metabolizers. Furthermore, secondary analyses examining four-level response of time to first cigarette after waking (>60, 31-60, 6-30, ≤5 minutes) demonstrated a robust effect of the metabolism metric in Collaborative Study on the Genetics of Alcoholism (P = 0.03) and Collaborative Genetic Study of Nicotine Dependence (P = 0.004), illustrating the important role of this measure of dependence. These findings highlight the complex role of CYP2A6 variation across different developmental stages of smoking behaviors.

A putative causal relationship between genetically determined female body shape and posttraumatic stress disorder

Background

The nature and underlying mechanisms of the observed increased vulnerability to posttraumatic stress disorder (PTSD) in women are unclear.

Methods

We investigated the genetic overlap of PTSD with anthropometric traits and reproductive behaviors and functions in women. The analysis was conducted using female-specific summary statistics from large genome-wide association studies (GWAS) and a cohort of 3577 European American women (966 PTSD cases and 2611 trauma-exposed controls). We applied a high-resolution polygenic score approach and Mendelian randomization analysis to investigate genetic correlations and causal relationships.

Results

We observed an inverse association of PTSD with genetically determined anthropometric traits related to body shape, independent of body mass index (BMI). The top association was related to BMI-adjusted waist circumference (WC_adj; R = –0.079, P < 0.001, Q = 0.011). We estimated a relative decrease of 64.6% (95% confidence interval = 27.5–82.7) in the risk of PTSD per 1-SD increase in WC_adj. MR-Egger regression intercept analysis showed no evidence of pleiotropic effects in this association (P_pleiotropy = 0.979). We also observed associations of genetically determined WC_adj with age at first sexual intercourse and number of sexual partners (P = 0.013 and P < 0.001, respectively).

Conclusions

There is a putative causal relationship between genetically determined female body shape and PTSD, which could be mediated by evolutionary mechanisms involved in human sexual behaviors.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0491-4) contains supplementary material, which is available to authorized users.

Targeted sequencing identifies genetic polymorphisms of flavin-containing monooxygenase genes contributing to susceptibility of nicotine dependence in European American and African American

BACKGROUND

Smoking is a leading cause of preventable death. Early studies based on samples of twins have linked the lifetime smoking practices to genetic predisposition. The flavin-containing monooxygenase (FMO) protein family consists of a group of enzymes that metabolize drugs and xenobiotics. Both FMO1 and FMO3 were potentially susceptible genes for nicotine metabolism process.

METHODS

In this study, we investigated the potential of FMO genes to confer risk of nicotine dependence via deep targeted sequencing in 2,820 study subjects comprising 1,583 nicotine dependents and 1,237 controls from European American and African American. Specifically, we focused on the two genomic segments including FMO1,FMO3, and pseudo gene FMO6P, and aimed to investigate the potential association between FMO genes and nicotine dependence. Both common and low-frequency/rare variants were analyzed using different algorithms. The potential functional significance of SNPs with association signal was investigated with relevant bioinformatics tools.

RESULTS

We identified different clusters of significant common variants in European (with most significant SNP rs6674596, p = .0004, OR = 0.67, MAF_EA = 0.14, FMO1) and African Americans (with the most significant SNP rs6608453, p = .001, OR = 0.64, MAF_AA = 0.1, FMO6P). No significant signals were identified through haplotype-based analyses. Gene network investigation indicated that both FMO1 and FMO3 have a strong relation with a variety of genes belonging to CYP gene families (with combined score greater than 0.9). Most of the significant variants identified were SNPs located within intron regions or with unknown functional significance, indicating a need for future work to understand the underlying functional significance of these signals.

CONCLUSIONS

Our findings indicated significant association between FMO genes and nicotine dependence. Replications of our findings in other ethnic groups were needed in the future. Most of the significant variants identified were SNPs located within intronic regions or with unknown functional significance, indicating a need for future work to understand the underlying functional significance of these signals.

KAT2B polymorphism identified for drug abuse in African Americans with regulatory links to drug abuse pathways in human prefrontal cortex

Drug abuse is a common and heritable set of disorders, but the underlying genetic factors are largely unknown. We conducted genome-wide association studies of drug abuse using 7 million imputed single nucleotide polymorphisms (SNPs) and insertions/deletions in African Americans (AAs; n = 3742) and European Americans (EAs; n = 6845). Cases were drawn from the Urban Health Study of street-recruited people, who injected drugs and reported abusing opioids, cocaine, marijuana, stimulants and/or other drugs 10 or more times in the past 30 days, and were compared with population controls. Independent replication testing was conducted in 755 AAs and 1131 EAs from the Genetic Association Information Network. An intronic SNP (rs9829896) in the K(lysine) acetyltransferase 2B (KAT2B) gene was significantly associated with drug abuse in AAs (P = 4.63 × 10^-8 ) and independently replicated in AAs (P = 0.0019). The rs9829896-C allele (frequency = 12%) had odds ratios of 0.68 and 0.53 across the AA cohorts: meta-analysis P = 3.93 × 10^-10 . Rs9829896-C was not associated with drug abuse across the EA cohorts: frequency = 36% and meta-analysis P = 0.12. Using dorsolateral prefrontal cortex data from the BrainCloud cohort, we found that rs9829896-C was associated with reduced KAT2B expression in AAs (n = 113, P = 0.050) but not EAs (n = 110, P = 0.39). KAT2B encodes a transcriptional regulator in the cyclic adenosine monophosphate and dopamine signaling pathways, and rs9829896-C was associated with expression of genes in these pathways: reduced CREBBP expression (P = 0.011) and increased OPRM1 expression (P = 0.016), both in AAs only. Our study identified the KAT2B SNP rs9829896 as having novel and biologically plausible associations with drug abuse and gene expression in AAs but not EAs, suggesting ancestry-specific effects.

Rare, low frequency and common coding variants in CHRNA5 and their contribution to nicotine dependence in European and African Americans

The common nonsynonymous variant rs16969968 in the α5 nicotinic receptor subunit gene (CHRNA5) is the strongest genetic risk factor for nicotine dependence in European Americans and contributes to risk in African Americans. To comprehensively examine whether other CHRNA5 coding variation influences nicotine dependence risk, we performed targeted sequencing on 1582 nicotine-dependent cases (Fagerström Test for Nicotine Dependence score⩾4) and 1238 non-dependent controls, with independent replication of common and low frequency variants using 12 studies with exome chip data. Nicotine dependence was examined using logistic regression with individual common variants (minor allele frequency (MAF)⩾0.05), aggregate low frequency variants (0.05>MAF⩾0.005) and aggregate rare variants (MAF<0.005). Meta-analysis of primary results was performed with replication studies containing 12 174 heavy and 11 290 light smokers. Next-generation sequencing with 180 × coverage identified 24 nonsynonymous variants and 2 frameshift deletions in CHRNA5, including 9 novel variants in the 2820 subjects. Meta-analysis confirmed the risk effect of the only common variant (rs16969968, European ancestry: odds ratio (OR)=1.3, P=3.5 × 10(-11); African ancestry: OR=1.3, P=0.01) and demonstrated that three low frequency variants contributed an independent risk (aggregate term, European ancestry: OR=1.3, P=0.005; African ancestry: OR=1.4, P=0.0006). The remaining 22 rare coding variants were associated with increased risk of nicotine dependence in the European American primary sample (OR=12.9, P=0.01) and in the same risk direction in African Americans (OR=1.5, P=0.37). Our results indicate that common, low frequency and rare CHRNA5 coding variants are independently associated with nicotine dependence risk. These newly identified variants likely influence the risk for smoking-related diseases such as lung cancer.

Replication of ZNF804A gene variant associations with risk of heroin addiction

Heroin addiction is heritable, but few specific genetic variants have been reproducibly associated with this disease. The zinc finger protein 804A (ZNF804A) gene is a biologically plausible susceptibility gene for heroin addiction, given its function as a transcription factor in human brain. Novel associations of two common ZNF804A single nucleotide polymorphisms (SNPs), rs7597593 and rs1344706, with heroin addiction have been reported in Han Chinese. Both SNPs have also been implicated for regulating ZNF804A expression in human brain, including the addiction-relevant dorsolateral prefrontal cortex. In this independent replication study, we tested the rs7597593 and rs1344706 SNP genotypes and their corresponding haplotypes for association with heroin addiction using cases drawn from the Urban Health Study and population controls: total N = 10 757 [7095 European Americans (EAs) and 3662 African Americans (AAs)]. We independently replicated both ZNF804A SNP associations in EAs: the rs7597593-T (P = 0.016) and rs1344706-A (P = 0.029) alleles both being associated with increased risk of heroin addiction, consistent with the prior report. Neither SNP was associated in AAs alone, but meta-analysis across both ancestry groups resulted in significant associations for rs1344706-A [P = 0.016, odds ratio (95% confidence interval) = 1.13 (1.02-1.25)] and its haplotype with rs7597593-T [P = 0.0067, odds ratio (95% confidence interval) = 1.16 (1.04-1.29)]. By showing consistent associations across independent studies and diverse ancestry groups, our study provides evidence that these two ZNF804A SNPs and their risk haplotype are among the few replicable genetic associations with heroin addiction.

Genome-wide meta-analysis reveals common splice site acceptor variant in CHRNA4 associated with nicotine dependence

We conducted a 1000 Genomes-imputed genome-wide association study (GWAS) meta-analysis for nicotine dependence, defined by the Fagerström Test for Nicotine Dependence in 17 074 ever smokers from five European-ancestry samples. We followed up novel variants in 7469 ever smokers from five independent European-ancestry samples. We identified genome-wide significant association in the alpha-4 nicotinic receptor subunit (CHRNA4) gene on chromosome 20q13: lowest P=8.0 × 10(-9) across all the samples for rs2273500-C (frequency=0.15; odds ratio=1.12 and 95% confidence interval=1.08-1.17 for severe vs mild dependence). rs2273500-C, a splice site acceptor variant resulting in an alternate CHRNA4 transcript predicted to be targeted for nonsense-mediated decay, was associated with decreased CHRNA4 expression in physiologically normal human brains (lowest P=7.3 × 10(-4)). Importantly, rs2273500-C was associated with increased lung cancer risk (N=28 998, odds ratio=1.06 and 95% confidence interval=1.00-1.12), likely through its effect on smoking, as rs2273500-C was no longer associated with lung cancer after adjustment for smoking. Using criteria for smoking behavior that encompass more than the single 'cigarettes per day' item, we identified a common CHRNA4 variant with important regulatory properties that contributes to nicotine dependence and smoking-related consequences.

Cis-Expression Quantitative Trait Loci Mapping Reveals Replicable Associations with Heroin Addiction in OPRM1

BACKGROUND

No opioid receptor, mu 1 (OPRM1) gene polymorphisms, including the functional single nucleotide polymorphism (SNP) rs1799971, have been conclusively associated with heroin/other opioid addiction, despite their biological plausibility. We used evidence of polymorphisms altering OPRM1 expression in normal human brain tissue to nominate and then test associations with heroin addiction.

METHODS

We tested 103 OPRM1 SNPs for association with OPRM1 messenger RNA expression in prefrontal cortex from 224 European Americans and African Americans of the BrainCloud cohort. We then tested the 16 putative cis-expression quantitative trait loci (cis-eQTL) SNPs for association with heroin addiction in the Urban Health Study and two replication cohorts, totaling 16,729 European Americans, African Americans, and Australians of European ancestry.

RESULTS

Four putative cis-eQTL SNPs were significantly associated with heroin addiction in the Urban Health Study (smallest p = 8.9 × 10(-5)): rs9478495, rs3778150, rs9384169, and rs562859. Rs3778150, located in OPRM1 intron 1, was significantly replicated (p = 6.3 × 10(-5)). Meta-analysis across all case-control cohorts resulted in p = 4.3 × 10(-8): the rs3778150-C allele (frequency = 16%-19%) being associated with increased heroin addiction risk. Importantly, the functional SNP allele rs1799971-A was associated with heroin addiction only in the presence of rs3778150-C (p = 1.48 × 10(-6) for rs1799971-A/rs3778150-C and p = .79 for rs1799971-A/rs3778150-T haplotypes). Lastly, replication was observed for six other intron 1 SNPs that had prior suggestive associations with heroin addiction (smallest p = 2.7 × 10(-8) for rs3823010).

CONCLUSIONS

Our findings show that common OPRM1 intron 1 SNPs have replicable associations with heroin addiction. The haplotype structure of rs3778150 and nearby SNPs may underlie the inconsistent associations between rs1799971 and heroin addiction.

A multiancestry study identifies novel genetic associations with CHRNA5 methylation in human brain and risk of nicotine dependence

Nicotine dependence is influenced by chromosome 15q25.1 single nucleotide polymorphisms (SNPs), including the missense SNP rs16969968 that alters function of the α5 nicotinic acetylcholine receptor (CHRNA5) and noncoding SNPs that regulate CHRNA5 mRNA expression. We tested for cis-methylation quantitative trait loci (cis-meQTLs) using SNP genotypes and DNA methylation levels measured across the IREB2-HYKK-PSMA4-CHRNA5-CHRNA3-CHRNB4 genes on chromosome 15q25.1 in the BrainCloud and Brain QTL cohorts [total N = 175 European-Americans and 65 African-Americans (AAs)]. We identified eight SNPs that were significantly associated with CHRNA5 methylation in prefrontal cortex: P ranging from 6.0 × 10(-10) to 5.6 × 10(-5). These SNP-methylation associations were also significant in frontal cortex, temporal cortex and pons: P ranging from 4.8 × 10(-12) to 3.4 × 10(-3). Of the eight cis-meQTL SNPs, only the intronic CHRNB4 SNP rs11636753 was associated with CHRNA5 methylation independently of the known SNP effects in prefrontal cortex, and it was the most significantly associated SNP with nicotine dependence across five independent cohorts (total N = 7858 European ancestry and 3238 AA participants): P = 6.7 × 10(-4), odds ratio (OR) [95% confidence interval (CI)] = 1.11 (1.05-1.18). The rs11636753 major allele (G) was associated with lower CHRNA5 DNA methylation, lower CHRNA5 mRNA expression and increased nicotine dependence risk. Haplotype analyses showed that rs11636753-G and the functional rs16969968-A alleles together increased risk of nicotine dependence more than each variant alone: P = 3.1 × 10(-12), OR (95% CI) = 1.32 (1.22-1.43). Our findings identify a novel regulatory SNP association with nicotine dependence and connect, for the first time, previously observed differences in CHRNA5 mRNA expression and nicotine dependence risk to underlying DNA methylation differences.

Novel genetic locus implicated for HIV-1 acquisition with putative regulatory links to HIV replication and infectivity: a genome-wide association study

Fifty percent of variability in HIV-1 susceptibility is attributable to host genetics. Thus identifying genetic associations is essential to understanding pathogenesis of HIV-1 and important for targeting drug development. To date, however, CCR5 remains the only gene conclusively associated with HIV acquisition. To identify novel host genetic determinants of HIV-1 acquisition, we conducted a genome-wide association study among a high-risk sample of 3,136 injection drug users (IDUs) from the Urban Health Study (UHS). In addition to being IDUs, HIV- controls were frequency-matched to cases on environmental exposures to enhance detection of genetic effects. We tested independent replication in the Women’s Interagency HIV Study (N=2,533). We also examined publicly available gene expression data to link SNPs associated with HIV acquisition to known mechanisms affecting HIV replication/infectivity. Analysis of the UHS nominated eight genetic regions for replication testing. SNP rs4878712 in FRMPD1 met multiple testing correction for independent replication (P=1.38x10^-4), although the UHS-WIHS meta-analysis p-value did not reach genome-wide significance (P=4.47x10^-7 vs. P<5.0x10^-8) Gene expression analyses provided promising biological support for the protective G allele at rs4878712 lowering risk of HIV: (1) the G allele was associated with reduced expression of FBXO10 (r=-0.49, P=6.9x10^-5); (2) FBXO10 is a component of the Skp1-Cul1-F-box protein E3 ubiquitin ligase complex that targets Bcl-2 protein for degradation; (3) lower FBXO10 expression was associated with higher BCL2 expression (r=-0.49, P=8x10^-5); (4) higher basal levels of Bcl-2 are known to reduce HIV replication and infectivity in human and animal in vitro studies. These results suggest new potential biological pathways by which host genetics affect susceptibility to HIV upon exposure for follow-up in subsequent studies.

Multiple distinct CHRNB3-CHRNA6 variants are genetic risk factors for nicotine dependence in African Americans and European Americans

AIMS

Studies have shown association between common variants in the α6-β3 nicotinic receptor subunit gene cluster and nicotine dependence in European ancestry populations. We investigate whether this generalizes to African Americans, whether the association is specific to nicotine dependence and whether this region contains additional genetic contributors to nicotine dependence.

DESIGN

We examined consistency of association across studies and race between the α6β3 nicotinic receptor subunit locus and nicotine, alcohol, marijuana and cocaine dependence in three independent studies.

SETTING

United States of America.

PARTICIPANTS

European Americans and African Americans from three case-control studies of substance dependence.

MEASUREMENTS

Subjects were evaluated using the Semi-Structured Assessment for the Genetics of Alcoholism. Nicotine dependence was determined using the Fagerström Test for Nicotine Dependence.

FINDINGS

The single nucleotide polymorphism rs13273442 was associated significantly with nicotine dependence across all three studies in both ancestry groups [odds ratio (OR) = 0.75, P = 5.8 × 10(-4) European Americans; OR = 0.80, P = 0.05 African Americans]. No other substance dependence was associated consistently with this variant in either group. Another SNP in the region, rs4952, remains modestly associated with nicotine dependence in the combined data after conditioning on rs13273442.

CONCLUSIONS

The common variant rs13273442 in the CHRNB3-CHNRA6 region is associated significantly with nicotine dependence in European Americans and African Americans across studies recruited for nicotine, alcohol and cocaine dependence. Although these data are modestly powered for other substances, our results provide no evidence that correlates of rs13273442 represent a general substance dependence liability. Additional variants probably account for some of the association of this region to nicotine dependence.

Distinct loci in the CHRNA5/CHRNA3/CHRNB4 gene cluster are associated with onset of regular smoking

Neuronal nicotinic acetylcholine receptor (nAChR) genes (CHRNA5/CHRNA3/CHRNB4) have been reproducibly associated with nicotine dependence, smoking behaviors, and lung cancer risk. Of the few reports that have focused on early smoking behaviors, association results have been mixed. This meta-analysis examines early smoking phenotypes and SNPs in the gene cluster to determine: (1) whether the most robust association signal in this region (rs16969968) for other smoking behaviors is also associated with early behaviors, and/or (2) if additional statistically independent signals are important in early smoking. We focused on two phenotypes: age of tobacco initiation (AOI) and age of first regular tobacco use (AOS). This study included 56,034 subjects (41 groups) spanning nine countries and evaluated five SNPs including rs1948, rs16969968, rs578776, rs588765, and rs684513. Each dataset was analyzed using a centrally generated script. Meta-analyses were conducted from summary statistics. AOS yielded significant associations with SNPs rs578776 (beta = 0.02, P = 0.004), rs1948 (beta = 0.023, P = 0.018), and rs684513 (beta = 0.032, P = 0.017), indicating protective effects. There were no significant associations for the AOI phenotype. Importantly, rs16969968, the most replicated signal in this region for nicotine dependence, cigarettes per day, and cotinine levels, was not associated with AOI (P = 0.59) or AOS (P = 0.92). These results provide important insight into the complexity of smoking behavior phenotypes, and suggest that association signals in the CHRNA5/A3/B4 gene cluster affecting early smoking behaviors may be different from those affecting the mature nicotine dependence phenotype.

Protocol for a collaborative meta-analysis of 5-HTTLPR, stress, and depression

BACKGROUND

Debate is ongoing about what role, if any, variation in the serotonin transporter linked polymorphic region (5-HTTLPR) plays in depression. Some studies report an interaction between 5-HTTLPR variation and stressful life events affecting the risk for depression, others report a main effect of 5-HTTLPR variation on depression, while others find no evidence for either a main or interaction effect. Meta-analyses of multiple studies have also reached differing conclusions.

METHODS/DESIGN

To improve understanding of the combined roles of 5-HTTLPR variation and stress in the development of depression, we are conducting a meta-analysis of multiple independent datasets. This coordinated approach utilizes new analyses performed with centrally-developed, standardized scripts. This publication documents the protocol for this collaborative, consortium-based meta-analysis of 5-HTTLPR variation, stress, and depression.

STUDY ELIGIBILITY CRITERIA

Our goal is to invite all datasets, published or unpublished, with 5-HTTLPR genotype and assessments of stress and depression for at least 300 subjects. This inclusive approach is to minimize potential impact from publication bias.

DATA SOURCES

This project currently includes investigators from 35 independent groups, providing data on at least N = 33,761 participants.The analytic plan was determined prior to starting data analysis. Analyses of individual study datasets will be performed by the investigators who collected the data using centrally-developed standardized analysis scripts to ensure a consistent analytical approach across sites. The consortium as a group will review and interpret the meta-analysis results.

DISCUSSION

Variation in 5-HTTLPR is hypothesized to moderate the response to stress on depression. To test specific hypotheses about the role of 5-HTTLPR variation on depression, we will perform coordinated meta-analyses of de novo results obtained from all available data, using variables and analyses determined a priori. Primary analyses, based on the original 2003 report by Caspi and colleagues of a GxE interaction will be supplemented by secondary analyses to help interpret and clarify issues ranging from the mechanism of effect to heterogeneity among the contributing studies. Publication of this protocol serves to protect this project from biased reporting and to improve the ability of readers to interpret the results of this specific meta-analysis upon its completion.

A glimpse into the future - Personalized medicine for smoking cessation

The devastating consequences of tobacco smoking for individuals and societies motivate studies to identify and understand the biological pathways that drive smoking behaviors, so that more effective preventions and treatments can be developed. Cigarette smokers respond to nicotine in different ways, with a small number of smokers remaining lifelong low-level smokers who never exhibit any symptoms of dependence, and a larger group becoming nicotine dependent. Whether or not a smoker transitions to nicotine dependence has clear genetic contributions, and variants in the genes encoding the α5-α3-β4 nicotinic receptor subunits most strongly contribute to differences in the risk for developing nicotine dependence among smokers. More recent work reveals a differential response to pharmacologic treatment for smoking cessation based on these same genetic variants in the α5-α3-β4 nicotinic receptor gene cluster. We anticipate a continuing acceleration of the translation of genetic discoveries into more successful treatment for smoking cessation. Given that over 400,000 people in the United States and over 5 million people world-wide die each year from smoking related illnesses, an improved understanding of the mechanisms underlying smoking behavior and smoking cessation must be a high public health priority so we can best intervene at both the public health level and the individual level. This article is part of a Special Issue entitled 'NIDA 40th Anniversary Issue'.

ANKK1, TTC12, and NCAM1 polymorphisms and heroin dependence: importance of considering drug exposure

CONTEXT

The genetic contribution to liability for opioid dependence is well established; identification of the responsible genes has proved challenging.

OBJECTIVE

To examine association of 1430 candidate gene single-nucleotide polymorphisms (SNPs) with heroin dependence, reporting here only the 71 SNPs in the chromosome 11 gene cluster (NCAM1, TTC12, ANKK1, DRD2) that include the strongest observed associations.

DESIGN

Case-control genetic association study that included 2 control groups (lacking an established optimal control group).

SETTING

Semistructured psychiatric interviews.

PARTICIPANTS

A total of 1459 Australian cases ascertained from opioid replacement therapy clinics, 531 neighborhood controls ascertained from economically disadvantaged areas near opioid replacement therapy clinics, and 1495 unrelated Australian Twin Registry controls not dependent on alcohol or illicit drugs selected from a twin and family sample.

MAIN OUTCOME MEASURE

Lifetime heroin dependence.

RESULTS

Comparison of cases with Australian Twin Registry controls found minimal evidence of association for all chromosome 11 cluster SNPs (P ≥ .01); a similar comparison with neighborhood controls revealed greater differences (P ≥ 1.8 × 10(-4)). Comparing cases (n = 1459) with the subgroup of neighborhood controls not dependent on illicit drugs (n = 340), 3 SNPs were significantly associated (correcting for multiple testing): ANKK1 SNP rs877138 (most strongly associated; odds ratio = 1.59; 95% CI, 1.32-1.92; P = 9.7 × 10(-7)), ANKK1 SNP rs4938013, and TTC12 SNP rs7130431. A similar pattern of association was observed when comparing illicit drug-dependent (n = 191) and nondependent (n = 340) neighborhood controls, suggesting that liability likely extends to nonopioid illicit drug dependence. Aggregate heroin dependence risk associated with 2 SNPs, rs877138 and rs4492854 (located in NCAM1), varied more than 4-fold (P = 2.7 × 10(-9) for the risk-associated linear trend).

CONCLUSIONS

Our results provide further evidence of association for chromosome 11 gene cluster SNPs with substance dependence, including extension of liability to illicit drug dependence. Our findings highlight the necessity of considering drug exposure history when selecting control groups for genetic investigations of illicit drug dependence.

Interpreting joint SNP analysis results: when are two distinct signals really two distinct signals?

In genetic association studies, much effort has focused on moving beyond the initial single-nucleotide polymorphism (SNP)-by-SNP analysis. One approach is to reanalyze a chromosomal region where an association has been detected, jointly analyzing the SNP thought to best represent that association with each additional SNP in the region. Such joint analyses may help identify additional, statistically independent association signals. However, it is possible for a single genetic effect to produce joint SNP results that would typically be interpreted as two distinct effects (e.g., both SNPs are significant in the joint model). We present a general approach that can (1) identify conditions under which a single variant could produce a given joint SNP result, and (2) use these conditions to identify variants from a list of known SNPs (e.g., 1000 Genomes) as candidates that could produce the observed signal. We apply this method to our previously reported joint result for smoking involving rs16969968 and rs588765 in CHRNA5. We demonstrate that it is theoretically possible for a joint SNP result suggestive of two independent signals to be produced by a single causal variant. Furthermore, this variant need not be highly correlated with the two tested SNPs or have a large odds ratio. Our method aids in interpretation of joint SNP results by identifying new candidate variants for biological causation that would be missed by traditional approaches. Also, it can connect association findings that may seem disparate due to lack of high correlations among the associated SNPs.

Imputation across genotyping arrays for genome-wide association studies: assessment of bias and a correction strategy

A great promise of publicly sharing genome-wide association data is the potential to create composite sets of controls. However, studies often use different genotyping arrays, and imputation to a common set of SNPs has shown substantial bias: a problem which has no broadly applicable solution. Based on the idea that using differing genotyped SNP sets as inputs creates differential imputation errors and thus bias in the composite set of controls, we examined the degree to which each of the following occurs: (1) imputation based on the union of genotyped SNPs (i.e., SNPs available on one or more arrays) results in bias, as evidenced by spurious associations (type 1 error) between imputed genotypes and arbitrarily assigned case/control status; (2) imputation based on the intersection of genotyped SNPs (i.e., SNPs available on all arrays) does not evidence such bias; and (3) imputation quality varies by the size of the intersection of genotyped SNP sets. Imputations were conducted in European Americans and African Americans with reference to HapMap phase II and III data. Imputation based on the union of genotyped SNPs across the Illumina 1M and 550v3 arrays showed spurious associations for 0.2 % of SNPs: ~2,000 false positives per million SNPs imputed. Biases remained problematic for very similar arrays (550v1 vs. 550v3) and were substantial for dissimilar arrays (Illumina 1M vs. Affymetrix 6.0). In all instances, imputing based on the intersection of genotyped SNPs (as few as 30 % of the total SNPs genotyped) eliminated such bias while still achieving good imputation quality.