Genome-wide association meta-analysis of nicotine metabolism and cigarette consumption measures in smokers of European descent

PharmVar GeneFocus: CYP2A6

The Pharmacogene Variation Consortium (PharmVar) provides nomenclature for the human CYP2A gene locus containing the highly polymorphic CYP2A6 gene. CYP2A6 plays a role in the metabolism of nicotine and various drugs. Thus, genetic variation can substantially contribute to the function of this enzyme and associated efficacy and safety. This GeneFocus provides an overview of the clinical significance of CYP2A6, including its genetic variation and function. We also highlight and discuss caveats in the identification and characterization of allelic variation of this complex pharmacogene, a prerequisite for accurate genotype determination and prediction of phenotype status.

Proteomic Networks and Related Genetic Variants Associated with Smoking and Chronic Obstructive Pulmonary Disease

Background

Studies have identified individual blood biomarkers associated with chronic obstructive pulmonary disease (COPD) and related phenotypes. However, complex diseases such as COPD typically involve changes in multiple molecules with interconnections that may not be captured when considering single molecular features.

Methods

Leveraging proteomic data from 3,173 COPDGene Non-Hispanic White (NHW) and African American (AA) participants, we applied sparse multiple canonical correlation network analysis (SmCCNet) to 4,776 proteins assayed on the SomaScan v4.0 platform to derive sparse networks of proteins associated with current vs. former smoking status, airflow obstruction, and emphysema quantitated from high-resolution computed tomography scans. We then used NetSHy, a dimension reduction technique leveraging network topology, to produce summary scores of each proteomic network, referred to as NetSHy scores. We next performed genome-wide association study (GWAS) to identify variants associated with the NetSHy scores, or network quantitative trait loci (nQTLs). Finally, we evaluated the replicability of the networks in an independent cohort, SPIROMICS.

Results

We identified networks of 13 to 104 proteins for each phenotype and exposure in NHW and AA, and the derived NetSHy scores significantly associated with the variable of interests. Networks included known (sRAGE, ALPP, MIP1) and novel molecules (CA10, CPB1, HIS3, PXDN) and interactions involved in COPD pathogenesis. We observed 7 nQTL loci associated with NetSHy scores, 4 of which remained after conditional analysis. Networks for smoking status and emphysema, but not airflow obstruction, demonstrated a high degree of replicability across race groups and cohorts.

Conclusions

In this work, we apply state-of-the-art molecular network generation and summarization approaches to proteomic data from COPDGene participants to uncover protein networks associated with COPD phenotypes. We further identify genetic associations with networks. This work discovers protein networks containing known and novel proteins and protein interactions associated with clinically relevant COPD phenotypes across race groups and cohorts.

CYP2A6 Activity and Cigarette Consumption Interact in Smoking-Related Lung Cancer Susceptibility

Cigarette smoke, containing both nicotine and carcinogens, causes lung cancer. However, not all smokers develop lung cancer, highlighting the importance of the interaction between host susceptibility and environmental exposure in tumorigenesis. Here, we aimed to delineate the interaction between metabolizing ability of tobacco carcinogens and smoking intensity in mediating genetic susceptibility to smoking-related lung tumorigenesis. Single-variant and gene-based associations of 43 tobacco carcinogen-metabolizing genes with lung cancer were analyzed using summary statistics and individual-level genetic data, followed by causal inference of Mendelian randomization, mediation analysis, and structural equation modeling. Cigarette smoke-exposed cell models were used to detect gene expression patterns in relation to specific alleles. Data from the International Lung Cancer Consortium (29,266 cases and 56,450 controls) and UK Biobank (2,155 cases and 376,329 controls) indicated that the genetic variant rs56113850 C>T located in intron 4 of CYP2A6 was significantly associated with decreased lung cancer risk among smokers (OR = 0.88, 95% confidence interval = 0.85-0.91, P = 2.18 × 10-16), which might interact (Pinteraction = 0.028) with and partially be mediated (ORindirect = 0.987) by smoking status. Smoking intensity accounted for 82.3% of the effect of CYP2A6 activity on lung cancer risk but entirely mediated the genetic effect of rs56113850. Mechanistically, the rs56113850 T allele rescued the downregulation of CYP2A6 caused by cigarette smoke exposure, potentially through preferential recruitment of transcription factor helicase-like transcription factor. Together, this study provides additional insights into the interplay between host susceptibility and carcinogen exposure in smoking-related lung tumorigenesis.

SIGNIFICANCE

The causal pathway connecting CYP2A6 genetic variability and activity, cigarette consumption, and lung cancer susceptibility in smokers highlights the need for behavior modification interventions based on host susceptibility for cancer prevention.

CYP2A6 associates with respiratory disease risk and younger age of diagnosis: a phenome-wide association Mendelian Randomization study

CYP2A6, a genetically variable enzyme, inactivates nicotine, activates carcinogens, and metabolizes many pharmaceuticals. Variation in CYP2A6 influences smoking behaviors and tobacco-related disease risk. This phenome-wide association study examined associations between a reconstructed version of our weighted genetic risk score (wGRS) for CYP2A6 activity with diseases in the UK Biobank (N = 395 887). Causal effects of phenotypic CYP2A6 activity (measured as the nicotine metabolite ratio: 3'-hydroxycotinine/cotinine) on the phenome-wide significant (PWS) signals were then estimated in two-sample Mendelian Randomization using the wGRS as the instrument. Time-to-diagnosis age was compared between faster versus slower CYP2A6 metabolizers for the PWS signals in survival analyses. In the total sample, six PWS signals were identified: two lung cancers and four obstructive respiratory diseases PheCodes, where faster CYP2A6 activity was associated with greater disease risk (Ps < 1 × 10-6). A significant CYP2A6-by-smoking status interaction was found (Psinteraction < 0.05); in current smokers, the same six PWS signals were found as identified in the total group, whereas no PWS signals were found in former or never smokers. In the total sample and current smokers, CYP2A6 activity causal estimates on the six PWS signals were significant in Mendelian Randomization (Ps < 5 × 10-5). Additionally, faster CYP2A6 metabolizer status was associated with younger age of disease diagnosis for the six PWS signals (Ps < 5 × 10-4, in current smokers). These findings support a role for faster CYP2A6 activity as a causal risk factor for lung cancers and obstructive respiratory diseases among current smokers, and a younger onset of these diseases. This research utilized the UK Biobank Resource.

An in-depth association analysis of genetic variants within nicotine-related loci: Meeting in middle of GWAS and genetic fine-mapping

In the last two decades of Genome-wide association studies (GWAS), nicotine-dependence-related genetic loci (e.g., nicotinic acetylcholine receptor - nAChR subunit genes) are among the most replicable genetic findings. Although GWAS results have reported tens of thousands of SNPs within these loci, further analysis (e.g., fine-mapping) is required to identify the causal variants. However, it is computationally challenging for existing fine-mapping methods to reliably identify causal variants from thousands of candidate SNPs based on the posterior inclusion probability. To address this challenge, we propose a new method to select SNPs by jointly modeling the SNP-wise inference results and the underlying structured network patterns of the linkage disequilibrium (LD) matrix. We use adaptive dense subgraph extraction method to recognize the latent network patterns of the LD matrix and then apply group LASSO to select causal variant candidates. We applied this new method to the UK biobank data to identify the causal variant candidates for nicotine addiction. Eighty-one nicotine addiction-related SNPs (i.e.,-log(p) > 50) of nAChR were selected, which are highly correlated (average r2>0.8) although they are physically distant (e.g., >200 kilobase away) and from various genes. These findings revealed that distant SNPs from different genes can show higher LD r2 than their neighboring SNPs, and jointly contribute to a complex trait like nicotine addiction.

Dimensional and transdiagnostic phenotypes in psychiatric genome-wide association studies

Genome-wide association studies (GWAS) provide biological insights into disease onset and progression and have potential to produce clinically useful biomarkers. A growing body of GWAS focuses on quantitative and transdiagnostic phenotypic targets, such as symptom severity or biological markers, to enhance gene discovery and the translational utility of genetic findings. The current review discusses such phenotypic approaches in GWAS across major psychiatric disorders. We identify themes and recommendations that emerge from the literature to date, including issues of sample size, reliability, convergent validity, sources of phenotypic information, phenotypes based on biological and behavioral markers such as neuroimaging and chronotype, and longitudinal phenotypes. We also discuss insights from multi-trait methods such as genomic structural equation modelling. These provide insight into how hierarchical 'splitting' and 'lumping' approaches can be applied to both diagnostic and dimensional phenotypes to model clinical heterogeneity and comorbidity. Overall, dimensional and transdiagnostic phenotypes have enhanced gene discovery in many psychiatric conditions and promises to yield fruitful GWAS targets in the years to come.

Associations between smoke exposure and kidney stones: results from the NHANES (2007-2018) and Mendelian randomization analysis

Purpose

It is currently controversial whether smoke exposure is associated with the risk of kidney stones. Herein, publicly available databases were combined to explore relationships with the risk of nephrolithiasis in terms of smoking status and serum cotinine concentrations.

Materials and methods

First, we conducted an observational study using data from 2007 to 2018, based on the National Health and Nutrition Examination Survey (NHANES) database. Univariate analysis, multivariate logistic regression, trend testing, restricted cubic spline (RCS), and multiple imputation (MI) were the main analytical methods of our study. Then, A Mendelian randomization (MR) analysis was performed to explore the causal relationship between serum cotinine and nephrolithiasis. Genetic instruments for serum cotinine and pooled data for kidney stones were derived from publicly available large-scale genome-wide association studies (GWAS). Inverse-variance weighting (IVW) was the primary method for our MR analysis.

Results

A total of 34,657 and 31,352 participants were included in the observational study based on smoking status and serum cotinine concentrations, respectively. Under full adjustment of covariates, current smokers had an increased risk of kidney stones compared to non-smokers [OR = 1.17 (1.04-1.31), P = 0.009, P for trend = 0.010]. Compared with serum cotinine of <0.05 ng/ml, serum cotinine levels of 0.05-2.99 ng/ml [OR = 1.15 (1.03-1.29), P = 0.013] and ≥3.00 ng/ml [OR = 1.22 (1.10-1.37), P < 0.001] were observed to have a higher risk of nephrolithiasis (P for trend < 0.001). In addition, a non-linear relationship between log2-transformed serum cotinine and the risk of nephrolithiasis was found (P for non-linearity = 0.028). Similar results were found when serum cotinine (log2 transformation) was used as a continuous variable [OR = 1.02 (1.01-1.03), P < 0.001] or complete data was used to analyze after MI. In the MR analysis, genetically predicted high serum cotinine was causally related to the high risk of nephrolithiasis [IVW: OR = 1.09 (1.00-1.19), P = 0.044].

Conclusion

Current smoking and high serum cotinine concentrations may be associated with an increased risk of kidney stones. Further research is needed to validate this relationship and explore its underlying mechanisms.

Inhibition of Nicotine Metabolism by Cannabidiol (CBD) and 7-Hydroxycannabidiol (7-OH-CBD)

Cannabis-based products have experienced notable increases in co-usage alongside tobacco products. Several cannabinoids exhibit inhibition of a number of cytochrome P450 (CYP) and UDP glucuronosyltransferase (UGT) enzymes, but few studies have examined their inhibition of enzymes involved in nicotine metabolism. The goal of the present study was to examine potential drug-drug interactions occurring in the nicotine metabolism pathway perpetrated by cannabidiol (CBD) and its active metabolite, 7-hydroxy-CBD (7-OH-CBD). The inhibitory effects of CBD and 7-OH-CBD were tested in microsomes from HEK293 cells overexpressing individual metabolizing enzymes and from human liver tissue. Assays with overexpressing microsomes demonstrated that CBD and 7-OH-CBD inhibited CYP-mediated nicotine metabolism. Binding-corrected IC50,u values for CBD inhibition of nicotine metabolism to cotinine and nornicotine, and cotinine metabolism to trans-3'-hydroxycotinine (3HC), were 0.27 ± 0.060, 0.23 ± 0.14, and 0.21 ± 0.14 μM, respectively, for CYP2A6; and 0.26 ± 0.17 and 0.029 ± 0.0050 μM for cotinine and nornicotine formation, respectively, for CYP2B6. 7-OH-CBD IC50,u values were 0.45 ± 0.18, 0.16 ± 0.08, and 0.78 ± 0.23 μM for cotinine, nornicotine, and 3HC formation, respectively, for CYP2A6, and 1.2 ± 0.44 and 0.11 ± 0.030 μM for cotinine and nornicotine formation, respectively, for CYP2B6. Similar IC50,u values were observed in HLM. Inhibition (IC50,u = 0.37 ± 0.06 μM) of 3HC to 3HC-glucuronide formation by UGT1A9 was demonstrated by CBD. Significant inhibition of nicotine metabolism pathways by CBD and 7-OH-CBD suggests that cannabinoids may inhibit nicotine metabolism, potentially impacting tobacco addiction and cessation.

The Use of Biomarkers to Guide Precision Treatment for Tobacco Use

This review summarizes the evidence to date on the development of biomarkers for personalizing the pharmacological treatment of combustible tobacco use. First, the latest evidence on FDA-approved medications is considered, demonstrating that, while these medications offer real benefits, they do not contribute to smoking cessation in approximately two-thirds of cases. Second, the case for using biomarkers to guide tobacco treatment is made based on the potential to increase medication effectiveness and uptake and reduce side effects. Next, the FDA framework of biomarker development is presented along with the state of science on biomarkers for tobacco treatment, including a review of the nicotine metabolite ratio, electroencephalographic event-related potentials, and other biomarkers utilized for risk feedback. We conclude with a discussion of the challenges and opportunities for the translation of biomarkers to guide tobacco treatment and propose priorities for future research.

Genetically predicted tobacco consumption and risk of intracranial aneurysm: a Mendelian randomization study

Several observational studies have suggested that tobacco consumption is a risk factor for intracranial aneurysms (IAs). We here genetically predict the causal association between specific smoking features and biomarkers for smokers and IA risk. The Mendelian randomization (MR) analysis considered summary statistics from the largest current genome-wide association studies of smoking and IA. The inverse-variance weighted (IVW) method, weighted median method, MR-RAPS, and multiple variants Mendelian randomization (MVMR) were performed to estimate the effect of different smoking features and drinking in IA. We observed significant causal effects of smoking on the risk of both aneurysmal subarachnoid hemorrhage (aSAH) and unruptured IA (uIA). The ORs of IAs based on the IVW method were 1.890 (95% CI 1.486-2.405) of ever smoking regularly. MVMR analysis afforded odds ratios of 1.685 (95% CI 1.136-2.501). In the further subgroup analysis, a similar causal relationship was observed in aSAH. Moreover, our analyses suggested that higher blood cotinine level and cadmium increases aSAH risk, and ORs were 1.235 (95%CI 1.009-1.186) and 1.235 (95%CI 1.046-1.458), respectively. Our study suggests that ever smoking regularly is associated with the IA risk, which includes both uIA and aSAH. Besides, higher blood cadmium and cotinine level may increases IA and aSAH risk. Thus, tobacco control should be promoted as primordial prevention for IAs, and screening for patients with a smoking history is emphasized.

Predicting nicotine metabolism across ancestries using genotypes

BACKGROUND

There is a need to match characteristics of tobacco users with cessation treatments and risks of tobacco attributable diseases such as lung cancer. The rate in which the body metabolizes nicotine has proven an important predictor of these outcomes. Nicotine metabolism is primarily catalyzed by the enzyme cytochrone P450 (CYP2A6) and CYP2A6 activity can be measured as the ratio of two nicotine metabolites: trans-3'-hydroxycotinine to cotinine (NMR). Measurements of these metabolites are only possible in current tobacco users and vary by biofluid source, timing of collection, and protocols; unfortunately, this has limited their use in clinical practice. The NMR depends highly on genetic variation near CYP2A6 on chromosome 19 as well as ancestry, environmental, and other genetic factors. Thus, we aimed to develop prediction models of nicotine metabolism using genotypes and basic individual characteristics (age, gender, height, and weight).

RESULTS

We identified four multiethnic studies with nicotine metabolites and DNA samples. We constructed a 263 marker panel from filtering genome-wide association scans of the NMR in each study. We then applied seven machine learning techniques to train models of nicotine metabolism on the largest and most ancestrally diverse dataset (N=2239). The models were then validated using the other three studies (total N=1415). Using cross-validation, we found the correlations between the observed and predicted NMR ranged from 0.69 to 0.97 depending on the model. When predictions were averaged in an ensemble model, the correlation was 0.81. The ensemble model generalizes well in the validation studies across ancestries, despite differences in the measurements of NMR between studies, with correlations of: 0.52 for African ancestry, 0.61 for Asian ancestry, and 0.46 for European ancestry. The most influential predictors of NMR identified in more than two models were rs56113850, rs11878604, and 21 other genetic variants near CYP2A6 as well as age and ancestry.

CONCLUSIONS

We have developed an ensemble of seven models for predicting the NMR across ancestries from genotypes and age, gender and BMI. These models were validated using three datasets and associate with nicotine dosages. The knowledge of how an individual metabolizes nicotine could be used to help select the optimal path to reducing or quitting tobacco use, as well as, evaluating risks of tobacco use.

Accuracy and applications of sequencing and genotyping approaches for CYP2A6 and homologous genes

OBJECTIVES

We evaluated multiple genotyping/sequencing approaches in a homologous region of chromosome 19, and investigated associations of two common 3'-UTR CYP2A6 variants with activity in vivo.

METHODS

Individuals (n = 1704) of European and African ancestry were phenotyped for the nicotine metabolite ratio (NMR), an index of CYP2A6 activity, and genotyped/sequenced using deep amplicon exon sequencing, SNP array, genotype imputation and targeted capture sequencing. Amplicon exon sequencing was the gold standard to which other methods were compared within-individual for CYP2A6, CYP2A7, CYP2A13, and CYP2B6 exons to identify highly discordant positions. Linear regression models evaluated the association of CYP2A6*1B and rs8192733 genotypes (coded additively) with logNMR.

RESULTS

All approaches were ≤2.6% discordant with the gold standard; discordant calls were concentrated at few positions. Fifteen positions were discordant in >10% of individuals, with 12 appearing in regions of high identity between homologous genes (e.g. CYP2A6 and CYP2A7). For six, allele frequencies in our study and online databases were discrepant, suggesting errors in online sources. In the European-ancestry group (n = 935), CYP2A6*1B and rs8192733 were associated with logNMR (P < 0.001). A combined model found main effects of both variants on increasing logNMR. Similar trends were found in those of African ancestry (n = 506).

CONCLUSION

Multiple genotyping/sequencing approaches used in this chromosome 19 region contain genotyping/sequencing errors, as do online databases. Gene-specific primers and SNP array probes must consider gene homology; short-read sequencing of related genes in a single reaction should be avoided. Using improved sequencing approaches, we characterized two gain-of-function 3'-UTR variants, including the relatively understudied rs8192733.

Examining the role of mitochondrial genetic variation in nicotine dependence

Nicotine dependence (ND) has a heritability rate of ∼50%, suggesting genetic factors contribute to underlying mechanisms. Here, we aimed to examine variants within both mtDNA and the nuclear genome to determine if mitochondrial genes are associated with ND. A total of 129 mtDNA SNPs and 1136 nuclear-encoded mitochondrial genes in a sample of N = 374 Caucasians were selected for analysis. Age of onset of first, occasional, and daily smoking and Fagerström Test for Nicotine Dependence were used as outcomes for the analysis. Linear regression was used to test common variants. Gene analyses were performed using MAGMA. One nuclear mitochondrial SNP, rs78417112 found in the HSD17B4 gene, was significantly associated with the age of onset of occasional smoking. Additionally, one nuclear mitochondrial gene, PRKACA, was significantly associated with age of onset of both first and occasional smoking. Replication testing of the mtDNA m.1700T>C SNP, nominally associated with age of onset of daily smoking, was available in the PNAT2 clinical trial (N = 930 Caucasians). A meta-analysis showed this SNP was associated with age of onset of daily smoking (p-value = 0.004). Overall, the findings suggest mitochondrial genetic variation may contribute to variability in smoking phenotypes, although replication in larger samples is required.

The Promise of Polygenic Risk Prediction in Smoking Cessation: Evidence From Two Treatment Trials

INTRODUCTION

Tobacco use disorder is a complex behavior with a strong genetic component. Genome-wide association studies (GWAS) on smoking behaviors allow for the creation of polygenic risk scores (PRSs) to approximate genetic vulnerability. However, the utility of smoking-related PRSs in predicting smoking cessation in clinical trials remains unknown.

AIMS AND METHODS

We evaluated the association between polygenic risk scores and bioverified smoking abstinence in a meta-analysis of two randomized, placebo-controlled smoking cessation trials. PRSs of smoking behaviors were created using the GWAS and Sequencing Consortium of Alcohol and Nicotine use (GSCAN) consortium summary statistics. We evaluated the utility of using individual PRS of specific smoking behavior versus a combined genetic risk that combines PRS of all four smoking behaviors. Study participants came from the Transdisciplinary Tobacco Use Research Centers (TTURCs) Study (1091 smokers of European descent), and the Genetically Informed Smoking Cessation Trial (GISC) Study (501 smokers of European descent).

RESULTS

PRS of later age of smoking initiation (OR [95% CI]: 1.20, [1.04-1.37], p = .0097) was significantly associated with bioverified smoking abstinence at end of treatment. In addition, the combined PRS of smoking behaviors also significantly predicted bioverified smoking abstinence (OR [95% CI] 0.71 [0.51-0.99], p = .045).

CONCLUSIONS

PRS of later age at smoking initiation may be useful in predicting smoking cessation at the end of treatment. A combined PRS may be a useful predictor for smoking abstinence by capturing the genetic propensity for multiple smoking behaviors.

IMPLICATIONS

There is a potential for polygenic risk scores to inform future clinical medicine, and a great need for evidence on whether these scores predict clinically meaningful outcomes. Our meta-analysis provides early evidence for potential utility of using polygenic risk scores to predict smoking cessation amongst smokers undergoing quit attempts, informing further work to optimize the use of polygenic risk scores in clinical care.

The Vulnerability to Methamphetamine Dependence and Genetics: A Case-Control Study Focusing on Genetic Polymorphisms at Chromosomal Region 5q31.3

OBJECTIVES

Methamphetamine (METH) is a central nervous psychostimulant and one of the most frequently used illicit drugs. Numerous genetic loci that influence complex traits, including alcohol abuse, have been discovered; however, genetic analyses for METH dependence remain limited. An increased histone deacetylase 3 (HDAC3) expression has been detected in Fos-positive neurons in the dorsomedial striatum following withdrawal after METH self-administration. Herein, we aimed to systematically investigate the contribution of HDAC3 to the vulnerability to METH dependence in a Han Chinese population.

METHODS

In total, we recruited 1,221 patients with METH dependence and 2,328 age- and gender-matched controls. For genotyping, we selected 14 single nucleotide polymorphisms (SNPs) located within ± 3 kb regions of HDAC3. The associations between genotyped genetic polymorphisms and the vulnerability to METH dependence were examined by single marker- and haplotype-based methods using PLINK. The effects of expression quantitative trait loci (eQTLs) on targeted gene expressions were investigated using the Genotype-Tissue Expression (GTEx) database.

RESULTS

The SNP rs14251 was identified as a significant association signal (χ² = 9.84, P = 0.0017). An increased risk of METH dependence was associated with the A allele (minor allele) of rs14251 [odds ratio (95% CI) = 1.25 (1.09-1.43)]. The results of in silico analyses suggested that SNP rs14251 could be a potential eQTL signal for FCHSD1, PCDHGB6, and RELL2, but not for HDAC3, in various human tissues.

CONCLUSION

We demonstrated that genetic polymorphism rs14251 located at 5q31.3 was significantly associated with the vulnerability to METH dependence in Han Chinese population.

Meta-Analysis of Transcriptome-Wide Association Studies across 13 Brain Tissues Identified Novel Clusters of Genes Associated with Nicotine Addiction

Genome-wide association studies (GWAS) have identified and reproduced thousands of diseases associated loci, but many of them are not directly interpretable due to the strong linkage disequilibrium among variants. Transcriptome-wide association studies (TWAS) incorporated expression quantitative trait loci (eQTL) cohorts as a reference panel to detect associations with the phenotype at the gene level and have been gaining popularity in recent years. For nicotine addiction, several important susceptible genetic variants were identified by GWAS, but TWAS that detected genes associated with nicotine addiction and unveiled the underlying molecular mechanism were still lacking. In this study, we used eQTL data from the Genotype-Tissue Expression (GTEx) consortium as a reference panel to conduct tissue-specific TWAS on cigarettes per day (CPD) over thirteen brain tissues in two large cohorts: UK Biobank (UKBB; number of participants (N) = 142,202) and the GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN; N = 143,210), then meta-analyzing the results across tissues while considering the heterogeneity across tissues. We identified three major clusters of genes with different meta-patterns across tissues consistent in both cohorts, including homogenous genes associated with CPD in all brain tissues; partially homogeneous genes associated with CPD in cortex, cerebellum, and hippocampus tissues; and, lastly, the tissue-specific genes associated with CPD in only a few specific brain tissues. Downstream enrichment analyses on each gene cluster identified unique biological pathways associated with CPD and provided important biological insights into the regulatory mechanism of nicotine dependence in the brain.

Multiethnic Prediction of Nicotine Biomarkers and Association With Nicotine Dependence

INTRODUCTION

The nicotine metabolite ratio and nicotine equivalents are measures of metabolism rate and intake. Genome-wide prediction of these nicotine biomarkers in multiethnic samples will enable tobacco-related biomarker, behavioral, and exposure research in studies without measured biomarkers.

AIMS AND METHODS

We screened genetic variants genome-wide using marginal scans and applied statistical learning algorithms on top-ranked genetic variants, age, ethnicity and sex, and, in additional modeling, cigarettes per day (CPD), (in additional modeling) to build prediction models for the urinary nicotine metabolite ratio (uNMR) and creatinine-standardized total nicotine equivalents (TNE) in 2239 current cigarette smokers in five ethnic groups. We predicted these nicotine biomarkers using model ensembles and evaluated external validity using dependence measures in 1864 treatment-seeking smokers in two ethnic groups.

RESULTS

The genomic regions with the most selected and included variants for measured biomarkers were chr19q13.2 (uNMR, without and with CPD) and chr15q25.1 and chr10q25.3 (TNE, without and with CPD). We observed ensemble correlations between measured and predicted biomarker values for the uNMR and TNE without (with CPD) of 0.67 (0.68) and 0.65 (0.72) in the training sample. We observed inconsistency in penalized regression models of TNE (with CPD) with fewer variants at chr15q25.1 selected and included. In treatment-seeking smokers, predicted uNMR (without CPD) was significantly associated with CPD and predicted TNE (without CPD) with CPD, time-to-first-cigarette, and Fagerström total score.

CONCLUSIONS

Nicotine metabolites, genome-wide data, and statistical learning approaches developed novel robust predictive models for urinary nicotine biomarkers in multiple ethnic groups. Predicted biomarker associations helped define genetically influenced components of nicotine dependence.

IMPLICATIONS

We demonstrate development of robust models and multiethnic prediction of the uNMR and TNE using statistical and machine learning approaches. Variants included in trained models for nicotine biomarkers include top-ranked variants in multiethnic genome-wide studies of smoking behavior, nicotine metabolites, and related disease. Association of the two predicted nicotine biomarkers with Fagerström Test for Nicotine Dependence items supports models of nicotine biomarkers as predictors of physical dependence and nicotine exposure. Predicted nicotine biomarkers may facilitate tobacco-related disease and treatment research in samples with genomic data and limited nicotine metabolite or tobacco exposure data.

Genetic contributions to alcohol use disorder treatment outcomes: a genome-wide pharmacogenomics study

Naltrexone can aid in reducing alcohol consumption, while acamprosate supports abstinence; however, not all patients with alcohol use disorder (AUD) benefit from these treatments. Here we present the first genome-wide association study of AUD treatment outcomes based on data from the COMBINE and PREDICT studies of acamprosate and naltrexone, and the Mayo Clinic CITA study of acamprosate. Primary analyses focused on treatment outcomes regardless of pharmacological intervention and were followed by drug-stratified analyses to identify treatment-specific pharmacogenomic predictors of acamprosate and naltrexone response. Treatment outcomes were defined as: (1) time until relapse to any drinking (TR) and (2) time until relapse to heavy drinking (THR; ≥ 5 drinks for men, ≥4 drinks for women in a day), during the first 3 months of treatment. Analyses were performed within each dataset, followed by meta-analysis across the studies (N = 1083 European ancestry participants). Single nucleotide polymorphisms (SNPs) in the BRE gene were associated with THR (min p = 1.6E-8) in the entire sample, while two intergenic SNPs were associated with medication-specific outcomes (naltrexone THR: rs12749274, p = 3.9E-8; acamprosate TR: rs77583603, p = 3.1E-9). The top association signal for TR (p = 7.7E-8) and second strongest signal in the THR (p = 6.1E-8) analysis of naltrexone-treated patients maps to PTPRD, a gene previously implicated in addiction phenotypes in human and animal studies. Leave-one-out polygenic risk score analyses showed significant associations with TR (p = 3.7E-4) and THR (p = 2.6E-4). This study provides the first evidence of a polygenic effect on AUD treatment response, and identifies genetic variants associated with potentially medication-specific effects on AUD treatment response.

Lymphoma-Associated Biomarkers Are Increased in Current Smokers in Twin Pairs Discordant for Smoking

Simple Summary

Smoking is associated with a moderate increased risk of Hodgkin and follicular lymphoma. To help understand why, we examined lymphoma-related biomarker levels among 134 smoking and non-smoking twins (67 pairs) ascertained from the Finnish Twin Cohort. We validated self-reported smoking history by measuring serum cotinine, a metabolite of nicotine, from previously collected frozen serum samples. In total, 27 immune biomarkers were assayed using the Luminex Multiplex platform (R & D Systems). We found that four immune response biomarkers were higher and one was lower among smoking compared to non-smoking twins. The strongest association was observed for CCL17/TARC, a biomarker elevated in Hodgkin lymphoma patients. Immune biomarker levels were similar in former smokers and non-smokers. Current smoking may increase levels of immune proteins that could partially explain the association between smoking and risk of certain lymphomas.

Abstract

Smoking is associated with a moderate increased risk of Hodgkin and follicular lymphoma. To understand why, we examined lymphoma-related biomarker levels among 134 smoking and non-smoking twins (67 pairs) ascertained from the Finnish Twin Cohort. Previously collected frozen serum samples were tested for cotinine to validate self-reported smoking history. In total, 27 immune biomarkers were assayed using the Luminex Multiplex platform (R & D Systems). Current and non-current smokers were defined by a serum cotinine concentration of >3.08 ng/mL and ≤3.08 ng/mL, respectively. Associations between biomarkers and smoking were assessed using linear mixed models to estimate beta coefficients and standard errors, adjusting for age, sex and twin pair as a random effect. There were 55 never smokers, 43 current smokers and 36 former smokers. CCL17/TARC, sgp130, haptoglobin, B-cell activating factor (BAFF) and monocyte chemoattractant protein-1 (MCP1) were significantly (p < 0.05) associated with current smoking and correlated with increasing cotinine concentrations (P_trend < 0.05). The strongest association was observed for CCL17/TARC (P_trend = 0.0001). Immune biomarker levels were similar in former and never smokers. Current smoking is associated with increased levels of lymphoma-associated biomarkers, suggesting a possible mechanism for the link between smoking and risk of these two B-cell lymphomas.

Analyses of nicotine metabolism biomarker genetics stratified by sex in African and European Americans

Nicotine is inactivated by the polymorphic CYP2A6 enzyme to cotinine and then to 3′hydroxycotinine. The Nicotine Metabolite Ratio (NMR; 3′hydroxycotinine/cotinine) is a heritable nicotine metabolism biomarker, varies with sex and ancestry, and influences smoking cessation and disease risk. We conducted sex-stratified genome-wide association studies of the NMR in European American (EA) and African American (AA) smokers (NCT01314001, NCT00666978). In EA females (n = 389) and males (n = 541), one significant (P < 5e−8) chromosome 19 locus was found (top variant: rs56113850, CYP2A6 (intronic), for C vs. T: females: beta = 0.67, P = 7.5e−22, 21.8% variation explained; males: beta = 0.75, P = 1.2e−37, 26.1% variation explained). In AA females (n = 503) and males (n = 352), the top variant was found on chromosome 19 but differed by sex (females: rs11878604, CYP2A6 (~ 16 kb 3′), for C vs. T: beta = − 0.71, P = 6.6e−26, 16.2% variation explained; males: rs3865454, CYP2A6 (~ 7 kb 3′), for G vs. T: beta = 0.64, P = 1.9e−19, 18.9% variation explained). In AA females, a significant region was found on chromosome 12 (top variant: rs12425845: P = 5.0e−9, TMEM132C (~ 1 Mb 5′), 6.1% variation explained) which was not significant in AA males. In AA males, significant regions were found on chromosomes 6 (top variant: rs9379805: P = 4.8e−9, SLC17A2 (~ 8 kb 5′), 8.0% variation explained) and 16 (top variant: rs77368288: P = 3.5e−8, ZNF469 (~ 92 kb 5′), 7.1% variation explained) which were not significant in AA females. Further investigation of these associations outside of chromosome 19 is required, as they did not replicate. Understanding how sex and ancestry influence nicotine metabolism genetics may improve personalized approaches for smoking cessation and risk prediction for tobacco-related diseases.

Functional validation of a finding from a mouse genome-wide association study shows that Azi2 influences the acute locomotor stimulant response to methamphetamine

In a previous genome-wide association study (GWAS) using outbred Carworth Farms White (CFW) mice, we identified a locus that influenced the stimulant response to methamphetamine and colocalized with an eQTL for Azi2. Based on those findings, we hypothesized that heritable differences in Azi2 expression were causally related to the differential response to methamphetamine. To test that hypothesis, we created a mutant Azi2 allele on an inbred C57BL/6J background. The mutant allele enhanced the locomotor response to methamphetamine. However, the GWAS had suggested that lower Azi2 would decrease the locomotor response to methamphetamine. We also sought to explore the mechanism by which Azi2 influenced methamphetamine sensitivity. A recent publication reported that the 3'UTR of Azi2 mRNA downregulates the expression of Slc6a3, which encodes the dopamine transporter, which is a key target of methamphetamine. We evaluated the relationship between Azi2, Azi2 3'UTR and Slc6a3 expression in the ventral tegmental area of wildtype, mutant Azi2 heterozygotes and mutant Azi2 homozygotes and in a new cohort of outbred CFW mice where both allele mapped in our prior GWAS were segregating. We did not observe any correlation between Azi2 and Slc6a3 in either cohort. However, RNA sequencing confirmed that the Azi2 mutation altered Azi2 expression and also revealed a number of potentially important genes and pathways that were regulated by Azi2, including the metabotropic glutamate receptor group III pathway and nicotinic acetylcholine receptor signaling pathway. Our results support a role for Azi2 in methamphetamine sensitivity; however, the exact mechanism does not appear to involve regulation of Slc6a3.

Genetic Susceptibility to Pneumonia: A GWAS Meta-Analysis Between the UK Biobank and FinnGen

Pneumonia is a respiratory condition with complex etiology. Host genetic variation is thought to contribute to individual differences in susceptibility and symptom manifestation. Here, we analyze pneumonia data from the UK Biobank (14,780 cases and 439,096 controls) and FinnGen (9980 cases and 86,519 controls) and perform a genomewide association study meta-analysis. We use gene-based tests, colocalization, genetic correlation, latent causal variable (LCV) and polygenic prediction in an independent Australian sample (N = 5595) to draw insights into the etiology of pneumonia risk. We identify two independent loci on chromosome 15 (lead single-nucleotide polymorphisms rs2009746 and rs76474922) to be associated with pneumonia (p < 5e-8). Gene-based tests revealed 18 genes in chromosomes 15, 16 and 9, including IL127, PBX3, ApoB receptor (APOBR) and smoking related genes CHRNA3/5, statistically associated with pneumonia. We observed genetic correlations between pneumonia and cardiorespiratory, psychiatric and inflammatory related traits. LCV analysis suggests a strong genetic causal relationship with cardiovascular health phenotypes. Polygenic risk scores for pneumonia significantly predicted self-reported pneumonia in an independent sample, albeit with a small effect size (OR = 1.11 95% CI [1.04, 1.19], p < .05). Sensitivity analyses suggested the associations in chromosome 15 are mediated by smoking history, but the associations in chromosomes 16 and 9, and polygenic prediction were robust to adjustment for smoking. Altogether, our results highlight common genetic variants, genes and potential pathways that contribute to individual differences in susceptibility to pneumonia, and advance our understanding of the genetic factors underlying heterogeneity in respiratory medical outcomes.

Effect of race and glucuronidation rates on the relationship between nicotine metabolite ratio and nicotine clearance

OBJECTIVES

To investigate if the nicotine metabolite ratio (NMR, the ratio of nicotine metabolites 3'-hydroxycotinine/cotinine) is a reliable phenotypic biomarker for nicotine clearance across races, and as a function of differences in the rate of nicotine, cotinine and 3'-hydroxycotinine glucuronidation and UGT genotypes.

METHODS

Participants [Caucasians (Whites), African Americans (Blacks) and Asian-Americans (Asians)] received an oral solution of deuterium-labeled nicotine and its metabolite cotinine. Plasma and saliva concentrations of nicotine and cotinine were used to determine oral clearances. Rates of glucuronidation were assessed from urine glucuronide/parent ratios, and UGT2B10 and UGT2B17 genotypes from DNA.

RESULTS

Among the 227 participants, 96 (42%) were White, 67 (30%) Asian and 64 (28%) Black. Compared to the other two races, Whites had higher nicotine and cotinine total oral clearance, Blacks had lower nicotine and cotinine glucuronidation rates and Asians had lower 3'-hydroxycotinine glucuronidation rates. A strong positive correlation (correlations coefficients 0.77-0.84; P < 0.001) between NMR and nicotine oral clearance was found for all three races, and NMR remained a strong predictor for the nicotine oral clearance while adjusting for race, sex and age. Neither the metabolite glucuronidation ratios nor the UGT genotypes had significant effects on the ability of NMR to predict nicotine oral clearance.

CONCLUSION

NMR appears to be a reliable phenotypic biomarker for nicotine clearance across races, glucuronidation phenotypes and genotypes. Racial differences in the relationships between NMR, smoking behaviors and addiction are unlikely to be related to an inadequate estimation of nicotine clearance on the basis of NMR.

Biochemistry of nicotine metabolism and its relevance to lung cancer

Nicotine is the key addictive constituent of tobacco. It is not a carcinogen, but it drives smoking and the continued exposure to the many carcinogens present in tobacco. The investigation into nicotine biotransformation has been ongoing for more than 60 years. The dominant pathway of nicotine metabolism in humans is the formation of cotinine, which occurs in two steps. The first step is cytochrome P450 (P450, CYP) 2A6–catalyzed 5′-oxidation to an iminium ion, and the second step is oxidation of the iminium ion to cotinine. The half-life of nicotine is longer in individuals with low P450 2A6 activity, and smokers with low activity often decrease either the intensity of their smoking or the number of cigarettes they use compared with those with “normal” activity. The effect of P450 2A6 activity on smoking may influence one's tobacco-related disease risk. This review provides an overview of nicotine metabolism and a summary of the use of nicotine metabolite biomarkers to define smoking dose. Some more recent findings, for example, the identification of uridine 5′-diphosphoglucuronosyltransferase 2B10 as the catalyst of nicotine N-glucuronidation, are discussed. We also describe epidemiology studies that establish the contribution of nicotine metabolism and CYP2A6 genotype to lung cancer risk, particularly with respect to specific racial/ethnic groups, such as those with Japanese, African, or European ancestry. We conclude that a model of nicotine metabolism and smoking dose could be combined with other lung cancer risk variables to more accurately identify former smokers at the highest risk of lung cancer and to intervene accordingly.

Integration of evidence across human and model organism studies: A meeting report

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

Clarifying the Genetic Influences on Nicotine Dependence and Quantity of Use in Cigarette Smokers

Nicotine dependence and smoking quantity are both robustly associated with the CHRNA5-A3-B4 gene cluster in the 15q25 region, and SNP rs16969968 in particular. The purpose of this paper is to use structural equation modeling techniques (SEM) to disentangle the complex pattern of relationships between rs16969968, nicotine quantity (as measured by the number of cigarettes an individual smokes per day; CPD) and nicotine dependence (as measured by the Fagerström Test for Nicotine Dependence; FTND). CPD is an indicator, but also a potential cause, of FTND, complicating the interpretation of associations between these constructs and requires a more detailed investigation than standard GWAS or general linear regression models can provide. FTND items and genotypes were collected in four samples, with a combined sample size of 5,373 respondents. A mega-analysis was conducted using a multiple group SEM approach to test competing hypotheses regarding the relationships between the SNP rs16969968, FTND and CPD. In the best fitting model, the FTND items loaded onto two correlated factors. The first, labeled "maintenance," assesses the motivation to maintain constant levels of nicotine through out the day. The second was labeled "urgency" as its items concern the urgency to restore nicotine levels after abstinence. We focus our attention on the "maintenance" factor, of which CPD was an indicator. The best fitting model included a negative feedback loop between the Maintenance factor and CPD. Accordingly, the motivation to maintain higher levels of nicotine increased the quantity of nicotine consumed, which subsequently decreases the maintenance motivation. The fact that the Maintenance-CPD feedback model fits the data best implies that there are at least two biological pathways that lead from rs16969968 to smoking behaviors. The model is consistent with a supply and demand system, which allows individuals to achieve a homeostatic equilibrium for their nicotine concentration.

Impact of CYP2A6 Activity on Nicotine Reinforcement and Cue-Reactivity in Daily Smokers

INTRODUCTION

Variation in CYP2A6, the primary enzyme responsible for nicotine metabolism, is associated with nicotine dependence, cigarette consumption, and abstinence outcomes in smokers. The impact of CYP2A6 activity on nicotine reinforcement and tobacco cue-reactivity, mechanisms that may contribute to these previous associations, has not been fully evaluated.

AIMS AND METHODS

CYP2A6 activity was indexed using 3 genetic approaches in 104 daily smokers completing forced-choice and cue-induced craving tasks assessing nicotine reinforcement and tobacco cue-reactivity, respectively. First, smokers were stratified by the presence or absence of reduced/loss-of-function CYP2A6 gene variants (normal vs. reduced metabolizers). As nicotine metabolite ratio (NMR) is a reliable biomarker of CYP2A6 activity, our second and third approaches used additional genetic variants identified in genome-wide association studies of NMR to create a weighted genetic risk score (wGRS) to stratify smokers (fast vs. slow metabolizers) and calculate a wGRS-derived NMR.

RESULTS

Controlling for race and sex, normal metabolizers (vs. reduced) selected a greater proportion of puffs from nicotine-containing cigarettes (vs. denicotinized) on the forced-choice task (p = .031). In confirmatory analyses, wGRS-based stratification (fast vs. slow metabolizers) produced similar findings. Additionally, wGRS-derived NMR, which correlated with actual NMR assessed in a subset of participants (n = 55), was positively associated with the proportion of puffs from nicotine-containing cigarettes controlling for race and sex (p = .015). None of the CYP2A6 indices were associated with tobacco cue-reactivity in minimally deprived smokers.

CONCLUSIONS

Findings suggest increased nicotine reinforcement is exhibited by smokers with high CYP2A6 activity, which may contribute to heavier smoking and poorer cessation outcomes previously reported in faster metabolizers.

IMPLICATIONS

CYP2A6 activity is a key determinant of smoking behavior and outcomes. Therefore, these findings support the targeting of CYP2A6 activity, either therapeutically or as a clinically relevant biomarker in a precision medicine approach, for tobacco use disorder treatment.

Genomics and epigenomics of addiction

Recent progress in the genomics and epigenomics of addiction has contributed to improving our understanding of this complex mental disorder's etiology, filling the gap between genes, environment, and behavior. We review the behavioral genetic studies reporting gene and environment interactions that explain the polygenetic contribution to the resilience and vulnerability to develop addiction. We discuss the evidence of polymorphic candidate genes that confer susceptibility to develop addiction as well as the studies of specific epigenetic marks that contribute to vulnerability and resilience to addictive-like behavior. A particular emphasis has been devoted to the miRNA changes that are considered potential biomarkers. The increasing knowledge about the technology required to alter miRNA expression may provide promising novel therapeutic tools. Finally, we give future directions for the field's progress in disentangling the connection between genes, environment, and behavior.

Pharmacogenetics factors influencing smoking cessation success; the importance of nicotine metabolism

Introduction: Smoking remains a worldwide epidemic, and despite an increase in public acceptance of the harms of tobacco use, it remains the leading cause of preventable death. It is estimated that up to 70% of all smokers express a desire to quit, but only 3-5% of them are successful.Areas covered: The goal of this review was to evaluate the current status of smoking cessation treatments and the feasibility of implementing personalized-medicine approaches to these pharmacotherapies. We evaluated the genetics associated with higher levels of nicotine addiction and follow with an analysis of the genetic variants that affect the nicotine metabolic ratio (NMR) and the FDA approved treatments for smoking cessation. We also highlighted the gaps in the process of translating current laboratory understanding into clinical practice, and the benefits of personalized treatment approaches for a successful smoking cessation strategy.Expert opinion: Evidence supports the use of tailored therapies to ensure that the most efficient treatments are utilized in an individual's smoking cessation efforts. An understanding of the genetic effects on the efficacy of individualized smoking cessation pharmacotherapies is key to smoking cessation, ideally utilizing a polygenetic risk score that considers all genetic variation.

Transferability of Ancestry-Specific and Cross-Ancestry CYP2A6 Activity Genetic Risk Scores in African and European Populations

The Nicotine Metabolite Ratio (NMR; 3-hydroxycotinine/cotinine), a highly heritable index of nicotine metabolic inactivation by the CYP2A6 enzyme, is associated with numerous smoking behaviors and diseases, as well as unique cessation outcomes. However, the NMR cannot be measured in nonsmokers, former smokers, or intermittent smokers, for example, in evaluating tobacco-related disease risk. Traditional pharmacogenetic groupings based on CYP2A6 * alleles capture a modest portion of NMR variation. We previously created a CYP2A6 weighted genetic risk score (wGRS) for European (EUR)-ancestry populations by incorporating independent signals from genome-wide association studies to capture a larger proportion of NMR variation. However, CYP2A6 genetic architecture is unique to ancestral populations. In this study, we developed and replicated an African-ancestry (AFR) wGRS, which captured 30-35% of the variation in NMR. We demonstrated model robustness against known environmental sources of NMR variation. Furthermore, despite the vast diversity within AFR populations, we showed that the AFR wGRS was consistent between different US geographical regions and unaltered by fine AFR population substructure. The AFR and EUR wGRSs can distinguish slow from normal metabolizers in their respective populations, and were able to reflect unique smoking cessation pharmacotherapy outcomes previously observed for the NMR. Additionally, we evaluated the utility of a cross-ancestry wGRS, and the capacity of EUR, AFR, and cross-ancestry wGRSs to predict the NMR within stratified or admixed AFR-EUR populations. Overall, our findings establish the clinical benefit of applying ancestry-specific wGRSs, demonstrating superiority of the AFR wGRS in AFRs.

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.