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.

Genotyping, characterization, and imputation of known and novel CYP2A6 structural variants using SNP array data.

CYP2A6 metabolically inactivates nicotine. Faster CYP2A6 activity is associated with heavier smoking and higher lung cancer risk. The CYP2A6 gene is polymorphic, including functional structural variants (SV) such as gene deletions (CYP2A6*4), duplications (CYP2A6*1 × 2), and hybrids with the CYP2A7 pseudogene (CYP2A6*12, CYP2A6*34). SVs are challenging to genotype due to their complex genetic architecture. Our aims were to develop a reliable protocol for SV genotyping, functionally phenotype known and novel SVs, and investigate the feasibility of CYP2A6 SV imputation from SNP array data in two ancestry populations. European- (EUR; n = 935) and African- (AFR; n = 964) ancestry individuals from smoking cessation trials were genotyped for SNPs using an Illumina array and for CYP2A6 SVs using Taqman copy number (CN) assays. SV-specific PCR amplification and Sanger sequencing was used to characterize a novel SV. Individuals with SVs were phenotyped using the nicotine metabolite ratio, a biomarker of CYP2A6 activity. SV diplotype and SNP array data were integrated and phased to generate ancestry-specific SV reference panels. Leave-one-out cross-validation was used to investigate the feasibility of CYP2A6 SV imputation. A minimal protocol requiring three Taqman CN assays for CYP2A6 SV genotyping was developed and known SV associations with activity were replicated. The first domain swap CYP2A6-CYP2A7 hybrid SV, CYP2A6*53, was identified, sequenced, and associated with lower CYP2A6 activity. In both EURs and AFRs, most SV alleles were identified using imputation (>70% and >60%, respectively); importantly, false positive rates were <1%. These results confirm that CYP2A6 SV imputation can identify most SV alleles, including a novel SV.

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.

Functional characterization of novel rare CYP2A6 variants and potential implications for clinical outcomes.

CYP2A6 activity, phenotyped by the nicotine metabolite ratio (NMR), is a predictor of several smoking behaviors, including cessation and smoking-related disease risk. The heritability of the NMR is 60-80%, yet weighted genetic risk scores (wGRSs) based on common variants explain only 30-35%. Rare variants (minor allele frequency <1%) are hypothesized to explain some of this missing heritability. We present two targeted sequencing studies where rare protein-coding variants are functionally characterized in vivo, in silico, and in vitro to examine this hypothesis. In a smoking cessation trial, 1687 individuals were sequenced; characterization measures included the in vivo NMR, in vitro protein expression, and metabolic activity measured from recombinant proteins. In a human liver bank, 312 human liver samples were sequenced; measures included RNA expression, protein expression, and metabolic activity from extracted liver tissue. In total, 38 of 47 rare coding variants identified were novel; characterizations ranged from gain-of-function to loss-of-function. On a population level, the portion of NMR variation explained by the rare coding variants was small (~1%). However, upon incorporation, the accuracy of the wGRS was improved for individuals with rare protein-coding variants (i.e., the residuals were reduced), and approximately one-third of these individuals (12/39) were re-assigned from normal to slow metabolizer status. Rare coding variants can alter an individual's CYP2A6 activity; their integration into wGRSs through precise functional characterization is necessary to accurately assess clinical outcomes and achieve precision medicine for all. Investigation into noncoding variants is warranted to further explain the missing heritability in the NMR.

The Role of Pharmacogenetics in Smoking.

Smoking continues to be the leading preventable contributor to death worldwide. Twin studies have suggested a significant genetic contribution underlying most smoking behaviors (40-70% heritability estimates). Candidate gene studies of smoking phenotypes have identified several pharmacogenes implicated in nicotine's pharmacokinetics (CYP2A6, CYP2B6, CYP2A13, FMOs, UGTs, and OCT2), and nicotine's pharmacodynamic response in the central nervous system (nicotinic acetylcholine receptors, as well as through the dopaminergic and serotonergic systems). Subsequent genome-wide association studies (GWAS) have confirmed the role of certain pharmacogenes through hypothesis-free approaches. Furthermore, pharmacogenes that alter the efficacy of smoking cessation pharmacotherapies, including nicotine replacement therapies, bupropion, and varenicline, may also impact quitting success. In this brief review we highlight the role of pharmacogenes in smoking behaviors, such as smoking status, consumption, nicotine dependence, spontaneous quitting, and altered abstinence to pharmacotherapies; We provide examples from initial candidate gene associations and subsequent GWAS. The genes CYP2A6 and the CHRNA5-A3-B4 confer the most replicated sources of genetic variation in smoking behaviors, likely due to their importance in nicotine's pharmacology. We will also provide examples of genetic scoring approaches, and the role of rare variants in explaining a portion of the missing heritability in smoking behaviors.

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.

PharmVar GeneFocus: CYP2B6.

The Pharmacogene Variation Consortium (PharmVar) catalogs star (*) allele nomenclature for the polymorphic human CYP2B6 gene. Genetic variation within the CYP2B6 gene locus impacts the metabolism or bioactivation of clinically important drugs. Of particular importance are efficacy and safety concerns regarding: efavirenz, which is used for the treatment of HIV type-1 infection; methadone, a mainstay in the treatment of opioid use disorder and as an analgesic; ketamine, used as an antidepressant and analgesic; and bupropion, which is prescribed to treat depression and for smoking cessation. This GeneFocus provides a comprehensive overview and summary of CYP2B6 and describes how haplotype information catalogued by PharmVar is utilized by the Pharmacogenomics Knowledgebase (PharmGKB) and the Clinical Pharmacogenetics Implementation Consortium (CPIC).

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.

Evaluation of a weighted genetic risk score for the prediction of biomarkers of CYP2A6 activity.

The nicotine metabolite ratio (NMR; 3-hydroxycotinine/cotinine) is an index of CYP2A6 activity. CYP2A6 is responsible for nicotine's metabolic inactivation and variation in the NMR/CYP2A6 is associated with several smoking behaviors. Our aim was to integrate established alleles and novel genome-wide association studies (GWAS) signals to create a weighted genetic risk score (wGRS) for the CYP2A6 gene for European-ancestry populations. The wGRS was compared with a previous CYP2A6 gene scoring approach designed for an alternative phenotype (C2/N2; cotinine-d2/(nicotine-d2 + cotinine-d2)). CYP2A6 genotypes and the NMR were assessed in European-ancestry participants. The wGRS training set included N = 933 smokers recruited to the Pharmacogenetics of Nicotine Addiction and Treatment clinical trial [NCT01314001]. The replication cohort included N = 196 smokers recruited to the Quit 2 Live clinical trial [NCT01836276]. Comparisons between the two CYP2A6 phenotypes and with fractional clearance were made in a laboratory-based pharmacokinetic study (N = 92 participants). In both the training and replication sets, the wGRS, which included seven CYP2A6 variants, explained 33.8% (P < 0.001) of the variance in NMR, providing improved predictive power to the NMR phenotype when compared with other CYP2A6 gene scoring approaches. NMR and C2/N2 were strongly correlated to nicotine clearance (ρ = 0.70 and ρ = 0.79, respectively; P < 0.001), and to one another (ρ = 0.82; P < 0.001); however reduced function genotypes occurred in slow NMR but throughout C2/N2. The wGRS was able to predict smoking quantity and nicotine intake, to discriminate between NMR slow and normal metabolizers (AUC = 0.79; P < 0.001), and to replicate previous NMR-stratified cessation outcomes showing unique treatment outcomes between metabolizer groups.

Pharmacogenomics of Nicotine Metabolism: Novel CYP2A6 and CYP2B6 Genetic Variation Patterns in Alaska Native and American Indian Populations.

INTRODUCTION: Alaska Native and American Indian (AN/AI) populations have higher tobacco use prevalence than other ethnic/racial groups. Pharmacogenetic testing to tailor tobacco cessation treatment may improve cessation rates. This study characterized polymorphic variations among AN/AI people in genes associated with metabolism of nicotine and drugs used for tobacco cessation. METHODS: Recruitment of AN/AI individuals represented six subgroups, five geographic subgroups throughout Alaska and a subgroup comprised of AIs from the lower 48 states living in Alaska. We sequenced the CYP2A6 and CYP2B6 genes to identify known and novel gain, reduced, and loss-of-function alleles, including structural variation (eg, gene deletions, duplications, and hybridizations). RESULTS: Variant allele frequencies differed substantially between AN/AI subgroups. The gene deletion CYP2A6*4 and reduced function CYP2A6*9 alleles were found at high frequency in Northern/Western subgroups and in Lower 48/Interior subgroups, respectively. The reduced function CYP2B6*6 allele was observed in all subgroups and a novel, predicted reduced function CYP2B6 variant was found at relatively high frequency in the Southeastern subgroup. CONCLUSIONS: Diverse CYP2A6 and CYP2B6 variation among the subgroups highlight the need for comprehensive pharmacogenetic testing to guide tobacco cessation therapy for AN/AI populations. IMPLICATIONS: Nicotine metabolism is largely determined by CYP2A6 genotype, and variation in CYP2A6 activity has altered the treatment success in other populations. These findings suggest pharmacogenetic-guided smoking cessation drug treatment could provide benefit to this unique population seeking tobacco cessation therapy.