Aspirin intake and head and neck cancer: A pooled analysis within the INHANCE consortium

BACKGROUND

Aspirin intake might be inversely associated with head and neck cancer (HNC). Thus, we investigated this relationship within the International Head and Neck Cancer Epidemiology (INHANCE) consortium.

METHODS

Four case-control studies within the INHANCE consortium were included (2024 cases, 4196 controls). Study-specific odds ratios (ORs) and 95% confidence intervals (CIs) were estimated using logistic regression and subsequently pooled with DerSimonian-Laird random-effects model. Nonlinearity of the relationship between duration of intake and HNC was modeled with fractional polynomials.

RESULTS

Aspirin was inversely associated with HNC overall (OR = 0.48; 95% CI: 0.26, 0.91). Results for laryngeal cancer were similar (OR = 0.54; 95% CI: 0.30, 0.96). Analysis on duration of intake confirmed findings for HNC overall, showing also inverse associations for oropharyngeal and laryngeal cancer.

CONCLUSIONS

This study suggests that aspirin intake may reduce the risk of HNC, driven mainly by decreases in risk for laryngeal and oropharyngeal cancer.

Smoking, Lung Cancer Stage, and Prognostic Factors-Findings from the National Lung Screening Trial

BACKGROUND

Low-dose computed tomography (LDCT) increases the early detection of lung cancer. Identifying modifiable behaviors that may affect tumor progression in LDCT-detected patients increases the likelihood of long-term survival and a good quality of life.

METHODS

We examined cigarette smoking behaviors on lung cancer stage, progression, and survival in 299 ever-smoking patients with low-dose CT-detected tumors from the National Lung Screening Trial. Univariate and multivariate Cox models were used to estimate the hazard ratio (HR) for smoking variables on survival time.

RESULTS

Current vs. former smokers and early morning smokers (≤5 min after waking, i.e., time to first cigarette (TTFC) ≤ 5 min) had more advanced-stage lung cancer. The adjusted HR for current vs. former smokers was 1.3 (95% confidence interval [CI] 0.911-1.98, p = 0.136) for overall survival (OS) and 1.3 (0.893-1.87, p = 0.1736) for progression-free survival (PFS). The univariate hazard ratios for TTFC ≤ 5 min vs. >5 min were 1.56 (1.1-2.2, p = 0.013) for OS and 1.53 (1.1-2.12, p = 0.01) for PFS. Among current smokers, the corresponding HRs for early TTFC were 1.78 (1.16-2.74, p = 0.0088) and 1.95 (1.29-2.95, p = 0.0016) for OS and PFS, respectively. In causal mediation analysis, the TTFC effect on survival time was mediated entirely through lung cancer stage.

CONCLUSION

The current findings indicate smoking behaviors at diagnosis may affect lung cancer stage and prognosis.

Cannabinoid-Induced Inhibition of Morphine Glucuronidation and the Potential for In Vivo Drug-Drug Interactions

Opioids are commonly prescribed for the treatment of chronic pain. Approximately 50% of adults who are prescribed opioids for pain co-use cannabis with their opioid treatment. Morphine is primarily metabolized by UDP-glucuronosyltransferase (UGT) 2B7 to an inactive metabolite, morphine-3-glucuronide (M3G), and an active metabolite, morphine-6-glucuronide (M6G). Previous studies have shown that major cannabis constituents including Δ9-tetrahydrocannabinol (THC) and cannabidiol (CBD) inhibit major UGT enzymes. To examine whether cannabinoids or their major metabolites inhibit morphine glucuronidation by UGT2B7, in vitro assays and mechanistic static modeling were performed with these cannabinoids and their major metabolites including 11-hydroxy-Δ9-tetrahydrocannabinol (11-OH-THC), 11-nor-9-carboxy-Δ9-tetrahydrocannabinol (11-COOH-THC), 7-hydroxy-cannabidiol (7-OH-CBD), and 7-carboxy-cannabidiol (7-COOH-CBD). In vitro assays with rUGT-overexpressing microsomes and human liver microsomes showed that THC and CBD and their metabolites inhibited UGT2B7-mediated morphine metabolism, with CBD and THC exhibiting the most potent Ki,u values (0.16 µM and 0.37 µM, respectively). Only 7-COOH-CBD exhibited no inhibitory activity against UGT2B7-mediated morphine metabolism. Static mechanistic modeling predicted an in vivo drug-drug interaction between morphine and THC after inhaled cannabis, and between THC, CBD, and 7-OH-CBD after oral consumption of cannabis. These data suggest that the co-use of these agents may lead to adverse drug events in humans.

Identification of genetically predicted DNA methylation markers associated with non-small cell lung cancer risk among 34,964 cases and 448,579 controls

BACKGROUND

Although the associations between genetic variations and lung cancer risk have been explored, the epigenetic consequences of DNA methylation in lung cancer development are largely unknown. Here, the genetically predicted DNA methylation markers associated with non-small cell lung cancer (NSCLC) risk by a two-stage case-control design were investigated.

METHODS

The genetic prediction models for methylation levels based on genetic and methylation data of 1595 subjects from the Framingham Heart Study were established. The prediction models were applied to a fixed-effect meta-analysis of screening data sets with 27,120 NSCLC cases and 27,355 controls to identify the methylation markers, which were then replicated in independent data sets with 7844 lung cancer cases and 421,224 controls. Also performed was a multi-omics functional annotation for the identified CpGs by integrating genomics, epigenomics, and transcriptomics and investigation of the potential regulation pathways.

RESULTS

Of the 29,894 CpG sites passing the quality control, 39 CpGs associated with NSCLC risk (Bonferroni-corrected p ≤ 1.67 × 10-6 ) were originally identified. Of these, 16 CpGs remained significant in the validation stage (Bonferroni-corrected p ≤ 1.28 × 10-3 ), including four novel CpGs. Multi-omics functional annotation showed nine of 16 CpGs were potentially functional biomarkers for NSCLC risk. Thirty-five genes within a 1-Mb window of 12 CpGs that might be involved in regulatory pathways of NSCLC risk were identified.

CONCLUSIONS

Sixteen promising DNA methylation markers associated with NSCLC were identified. Changes of the methylation level at these CpGs might influence the development of NSCLC by regulating the expression of genes nearby.

PLAIN LANGUAGE SUMMARY

The epigenetic consequences of DNA methylation in lung cancer development are still largely unknown. This study used summary data of large-scale genome-wide association studies to investigate the associations between genetically predicted levels of methylation biomarkers and non-small cell lung cancer risk at the first time. This study looked at how well larotrectinib worked in adult patients with sarcomas caused by TRK fusion proteins. These findings will provide a unique insight into the epigenetic susceptibility mechanisms of lung cancer.

Lung Cancer in Ever- and Never-Smokers: Findings from Multi-Population GWAS Studies

BACKGROUND

Clinical, molecular, and genetic epidemiology studies displayed remarkable differences between ever- and never-smoking lung cancer.

METHODS

We conducted a stratified multi-population (European, East Asian, and African descent) association study on 44,823 ever-smokers and 20,074 never-smokers to identify novel variants that were missed in the non-stratified analysis. Functional analysis including expression quantitative trait loci (eQTL) colocalization and DNA damage assays, and annotation studies were conducted to evaluate the functional roles of the variants. We further evaluated the impact of smoking quantity on lung cancer risk for the variants associated with ever-smoking lung cancer.

RESULTS

Five novel independent loci, GABRA4, intergenic region 12q24.33, LRRC4C, LINC01088, and LCNL1 were identified with the association at two or three populations (P < 5 × 10-8). Further functional analysis provided multiple lines of evidence suggesting the variants affect lung cancer risk through excessive DNA damage (GABRA4) or cis-regulation of gene expression (LCNL1). The risk of variants from 12 independent regions, including the well-known CHRNA5, associated with ever-smoking lung cancer was evaluated for never-smokers, light-smokers (packyear ≤ 20), and moderate-to-heavy-smokers (packyear > 20). Different risk patterns were observed for the variants among the different groups by smoking behavior.

CONCLUSIONS

We identified novel variants associated with lung cancer in only ever- or never-smoking groups that were missed by prior main-effect association studies.

IMPACT

Our study highlights the genetic heterogeneity between ever- and never-smoking lung cancer and provides etiologic insights into the complicated genetic architecture of this deadly cancer.

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.

Cannabinoid-Induced Stereoselective Inhibition of R-S-Oxazepam Glucuronidation: Cannabinoid-Oxazepam Drug Interactions

Benzodiazepines (BZDs) such as oxazepam are commonly prescribed depressant drugs known for their anxiolytic, hypnotic, muscle relaxant, and anticonvulsant effects and are frequently used in conjunction with other illicit drugs including cannabis. Oxazepam is metabolized in an enantiomeric-specific manner by glucuronidation, with S-oxazepam metabolized primarily by UGT2B15 and R-oxazepam glucuronidation mediated by both UGT 1A9 and 2B7. The goal of the present study was to evaluate the potential inhibitory effects of major cannabinoids, Δ9-tetrahydrocannabinol (THC) and cannabidiol (CBD), and major THC metabolites, 11-hydroxy-Δ9-tetrahydrocannabinol (11-OH-THC) and 11-nor-9-carboxy-Δ9-tetrahydrocannabinol (11-COOH-THC), on the UGT-mediated metabolism of R- and S-oxazepam. The cannabinoids and metabolites were screened as inhibitors of R- and S-oxazepam glucuronidation in microsomes isolated from HEK293 cells overexpressing individual UGT enzymes (rUGTs). The IC50 values were determined in human liver microsomes (HLM), human kidney microsomes (HKM), and rUGTs and utilized to estimate the nonspecific, binding-corrected Ki (Ki,u) values and predict the area under the concentration-time curve ratio (AUCR). The estimated Ki,u values observed in HLM for S- and R-oxazepam glucuronidation by CBD, 11-OH-THC, and THC were in the micromolar range (0.82 to 3.7 µM), with the Ki,u values observed for R-oxazepam glucuronidation approximately 2- to 5-fold lower as compared to those observed for S-oxazepam glucuronidation. The mechanistic static modeling predicted a potential clinically significant interaction between oral THC and CBD with oxazepam, with the AUCR values ranging from 1.25 to 3.45. These data suggest a pharmacokinetic drug-drug interaction when major cannabinoids like CBD or THC and oxazepam are concurrently administered.

Impact of individual level uncertainty of lung cancer polygenic risk score (PRS) on risk stratification

BACKGROUND

Although polygenic risk score (PRS) has emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS), the individual-level accuracy of lung cancer PRS and the extent to which its impact on subsequent clinical applications remains largely unexplored.

METHODS

Lung cancer PRSs and confidence/credible interval (CI) were constructed using two statistical approaches for each individual: (1) the weighted sum of 16 GWAS-derived significant SNP loci and the CI through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes), among 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium. Individuals were classified into different genetic risk subgroups based on the relationship between their own PRS mean/PRS CI and the population level threshold.

RESULTS

Considerable variances in PRS point estimates at the individual level were observed for both methods, with an average standard deviation (s.d.) of 0.12 for PRS-16-CV and a much larger s.d. of 0.88 for PRS-Bayes. Using PRS-16-CV, only 25.0% of individuals with PRS point estimates in the lowest decile of PRS and 16.8% in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individuals. Only 19% of the individuals were concordantly identified as having high genetic risk (> 90th percentile) using the two PRS estimators. An increased relative risk of lung cancer comparing the highest PRS percentile to the lowest was observed when taking the CI into account (OR = 2.73, 95% CI: 2.12-3.50, P-value = 4.13 × 10-15) compared to using PRS-16-CV mean (OR = 2.23, 95% CI: 1.99-2.49, P-value = 5.70 × 10-46). Improved risk prediction performance with higher AUC was consistently observed in individuals identified by PRS-16-CV CI, and the best performance was achieved by incorporating age, gender, and detailed smoking pack-years (AUC: 0.73, 95% CI = 0.72-0.74).

CONCLUSIONS

Lung cancer PRS estimates using different methods have modest correlations at the individual level, highlighting the importance of considering individual-level uncertainty when evaluating the practical utility of PRS.

Hydrocodone, Oxycodone, and Morphine Metabolism and Drug-Drug Interactions

Awareness of drug interactions involving opioids is critical for patient treatment as they are common therapeutics used in numerous care settings, including both chronic and disease-related pain. Not only do opioids have narrow therapeutic indexes and are extensively used, but they have the potential to cause severe toxicity. Opioids are the classical pain treatment for patients who suffer from moderate to severe pain. More importantly, opioids are often prescribed in combination with multiple other drugs, especially in patient populations who typically are prescribed a large drug regimen. This review focuses on the current knowledge of common opioid drug-drug interactions (DDIs), focusing specifically on hydrocodone, oxycodone, and morphine DDIs. The DDIs covered in this review include pharmacokinetic DDI arising from enzyme inhibition or induction, primarily due to inhibition of cytochrome p450 enzymes (CYPs). However, opioids such as morphine are metabolized by uridine-5'-diphosphoglucuronosyltransferases (UGTs), principally UGT2B7, and glucuronidation is another important pathway for opioid-drug interactions. This review also covers several pharmacodynamic DDI studies as well as the basics of CYP and UGT metabolism, including detailed opioid metabolism and the potential involvement of metabolizing enzyme gene variation in DDI. Based upon the current literature, further studies are needed to fully investigate and describe the DDI potential with opioids in pain and related disease settings to improve clinical outcomes for patients. SIGNIFICANCE STATEMENT: A review of the literature focusing on drug-drug interactions involving opioids is important because they can be toxic and potentially lethal, occurring through pharmacodynamic interactions as well as pharmacokinetic interactions occurring through inhibition or induction of drug metabolism.

Biochemical and genetic biomarkers associated with nicotine dependence in Mexican smokers

Cigarette smoking remains an important health concern and is still a leading cause of preventable mortality. Nicotine is the substance responsible for sustained tobacco use and dependence. Identification of biomarkers underlying nicotine dependence behavior is important to identify people at risk for this dependence. In the present study, we identified biochemical and genetic biomarkers of nicotine dependence detected by the Fagerström Test for Nicotine Dependence (FTDN) in Mexican smokers. The nicotine metabolites nicotine-N'-oxide, trans-3'-hydroxycotinine-glucuronide (3HC-O-Gluc), and nicotine-N-Gluc (Gluc) were useful to differentiate nicotine-dependent from non-dependent subjects (p < .0001) with an area under the curve (AUC) of 0.7818. Genetic variants in CYP2A6, FMO3, and UGT2B7 (rs2431413, rs28363545, and rs7439326, respectively) were associated with nicotine dependence (p = .03, p = .01, p = .01, respectively). Variations in the enzymatic activity of CYP2A6 were associated with altered nicotine-N'-oxide and 3HC-O-Gluc levels. Decreased urinary levels of 3HC-O-Gluc and increased nicotine-N'-oxide were associated with a decrease in the functional activity of CYP2A6. A strong positive correlation was observed between the ratio of urinary 3HC/cotinine, a measure of CYP2A6 activity, and the levels of 3HC-O-Gluc (p < .0001, r = .6835), while a strong negative correlation was observed with nicotine-N'-oxide (p < .0001, r = .6522) in nicotine-dependent subjects. No correlations were observed in non-nicotine-dependent subjects. These data suggest that particular urinary nicotine metabolites and genetic variants involved in nicotine metabolism are useful to identify subjects with nicotine dependence in the Mexican population.

Candidate pathway analysis of surfactant proteins identifies CTSH and SFTA2 that influences lung cancer risk

Pulmonary surfactant is a lipoprotein synthesized and secreted by alveolar type II cells in lung. We evaluated the associations between 200,139 single nucleotide polymorphisms (SNPs) of 40 surfactant-related genes and lung cancer risk using genotyped data from two independent lung cancer genome-wide association studies. Discovery data included 18,082 cases and 13,780 controls of European ancestry. Replication data included 1,914 cases and 3,065 controls of European descent. Using multivariate logistic regression, we found novel SNPs in surfactant-related genes CTSH [rs34577742 C > T, odds ratio (OR) = 0.90, 95% confidence interval (CI) = 0.89-0.93, P = 7.64 × 10-9] and SFTA2 (rs3095153 G > A, OR = 1.16, 95% CI = 1.10-1.21, P = 1.27 × 10-9) associated with overall lung cancer in the discovery data and validated in an independent replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.80-0.96, P = 5.76 × 10-3) and SFTA2 (rs3095153 G > A, OR = 1.14, 95% CI = 1.01-1.28, P = 3.25 × 10-2). Among ever smokers, we found SNPs in CTSH (rs34577742 C > T, OR = 0.89, 95% CI = 0.85-0.92, P = 1.94 × 10-7) and SFTA2 (rs3095152 G > A, OR = 1.20, 95% CI = 1.14-1.27, P = 4.25 × 10-11) associated with overall lung cancer in the discovery data and validated in the replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.79-0.97, P = 1.64 × 10-2) and SFTA2 (rs3095152 G > A, OR = 1.15, 95% CI = 1.01-1.30, P = 3.81 × 10-2). Subsequent transcriptome-wide association study using expression weights from a lung expression quantitative trait loci study revealed genes most strongly associated with lung cancer are CTSH (PTWAS = 2.44 × 10-4) and SFTA2 (PTWAS = 2.32 × 10-6).

Mosaic Chromosomal Alterations Are Associated With Increased Lung Cancer Risk: Insight From the INTEGRAL-ILCCO Cohort Analysis

INTRODUCTION

Mosaic chromosomal alterations (mCAs) detected in white blood cells represent a type of clonal hematopoiesis (CH) that is understudied compared with CH-related somatic mutations. A few recent studies indicated their potential link with nonhematological cancers, especially lung cancer.

METHODS

In this study, we investigated the association between mCAs and lung cancer using the high-density genotyping data from the OncoArray study of INTEGRAL-ILCCO, the largest single genetic study of lung cancer with 18,221 lung cancer cases and 14,825 cancer-free controls.

RESULTS

We identified a comprehensive list of autosomal mCAs, ChrX mCAs, and mosaic ChrY (mChrY) losses from these samples. Autosomal mCAs were detected in 4.3% of subjects, in addition to ChrX mCAs in 3.6% of females and mChrY losses in 9.6% of males. Multivariable logistic regression analysis indicated that the presence of autosomal mCAs in white blood cells was associated with an increased lung cancer risk after adjusting for key confounding factors, including age, sex, smoking status, and race. This association was mainly driven by a specific type of mCAs: copy-neutral loss of heterozygosity on autosomal chromosomes. The association between autosome copy-neutral loss of heterozygosity and increased risk of lung cancer was further confirmed in two major histologic subtypes, lung adenocarcinoma and squamous cell carcinoma. In addition, we observed a significant increase of ChrX mCAs and mChrY losses in smokers compared with nonsmokers and racial differences in certain types of mCA events.

CONCLUSIONS

Our study established a link between mCAs in white blood cells and increased risk of lung cancer.

Genome-wide association study of lung adenocarcinoma in East Asia and comparison with a European population

Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (Pinteraction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

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.

Impact of Genetic Variants in the Nicotine Metabolism Pathway on Nicotine Metabolite Levels in Smokers

BACKGROUND

Nicotine metabolism is a major factor in nicotine dependence, with approximately 70% to 80% of nicotine metabolized to cotinine in Caucasians. Cotinine formation is catalyzed primarily by CYP2A6, which also converts cotinine to trans-3'-hydroxycotinine (3HC). The goal of the present study was to examine the effects of CYP2A6 deficiency on nicotine metabolism profiles in vivo and the importance of genetic variants in nicotine-metabolizing enzyme genes on urinary nicotine metabolites levels.

METHODS

Urine samples from 722 smokers who participated in the Singapore Chinese Health Study were analyzed using UPLC-MS/MS to detect nicotine and eight of its urinary metabolites, and a total of 58 variants in 12 genes involved in nicotine metabolism were investigated in 475 of these subjects with informative genotyping data.

RESULTS

Urine samples stratified by the ratio of 3HC/cotinine exhibited a 7-fold increase in nicotine-N'-oxide, a 6-fold increase in nicotine-Glucuronide (Gluc), and a 5-fold decrease in 3HC-Gluc when comparing the lower versus upper 3HC/cotinine ventiles. Significant (P < 0.0001) associations were observed between functional metabolizing enzyme genotypes and levels of various urinary nicotine metabolites, including CYP2A6 genotype and levels of nicotine, nicotine-Gluc, nicotine-N'-oxide and 3HC, UGT2B10 genotype and levels of cotinine, nicotine-Gluc and cotinine-Gluc, UGT2B17 genotype and levels of 3HC-Gluc, FMO3 genotype and levels of nicotine-N'-oxide, and CYP2B6 genotype and levels of nicotine-N'-oxide and 4-hydroxy-4-(3-pyridyl)-butanoic acid.

CONCLUSIONS

These data suggest that several pathways are important in nicotine metabolism.

IMPACT

Genotype differences in several nicotine-metabolizing enzyme pathways may potentially lead to differences in nicotine dependence and smoking behavior and cessation.

CYP2C19 Plays a Major Role in the Hepatic N-Oxidation of Cotinine

The primary mode of metabolism of nicotine is via the formation of cotinine by the enzyme CYP2A6. Cotinine undergoes further CYP2A6-mediated metabolism by hydroxylation to 3-hydroxycotinine and norcotinine, but can also form cotinine-N-glucuronide and cotinine-N-oxide (COX). The goal of this study was to investigate the enzymes that catalyze COX formation and determine whether genetic variation in these enzymes may affect this pathway. Specific inhibitors of major hepatic cytochrome P450 (P450) enzymes were used in cotinine-N-oxidation reactions using pooled human liver microsomes (HLMs). COX formation was monitored by ultrahigh-pressure liquid chromatography-tandem mass spectrometry and enzyme kinetic analysis was performed using microsomes from P450-overexpressing human embryonic kidney 293 (HEK293) cell lines. Genotype-phenotype analysis was performed in a panel of 113 human liver specimens. Inhibition of COX formation was only observed in HLMs when using inhibitors of CYP2A6, CYP2B6, CYP2C19, CYP2E1, and CYP3A4. Microsomes from cells overexpressing CYP2A6 or CYP2C19 exhibited similar N-oxidation activity against cotinine, with maximum reaction rate over Michaelis constant values (intrinsic clearance) of 4.4 and 4.2 nL/min/mg, respectively. CYP2B6-, CYP2E1-, and CYP3A4-overexpressing microsomes were also active in COX formation. Significant associations (P < 0.05) were observed between COX formation and genetic variants in CYP2C19 (*2 and *17 alleles) in HLMs. These results demonstrate that genetic variants in CYP2C19 are associated with decreased COX formation, potentially affecting the relative levels of cotinine in the plasma or urine of smokers and ultimately affecting recommended smoking cessation therapies. SIGNIFICANCE STATEMENT: This study is the first to elucidate the enzymes responsible for cotinine-N-oxide formation and genetic variants that affect this biological pathway. Genetic variants in CYP2C19 have the potential to modify nicotine metabolic ratio in smokers and could affect pharmacotherapeutic decisions for smoking cessation treatments.

Quantification of Plasma 8-Isoprostane by High-Performance Liquid Chromatography with Tandem Mass Spectrometry in a Case-Control Study of Lung Cancer

AIM

8-iso-prostaglandin F2α is a biomarker of lipid peroxidation, and one of the most commonly used measures of oxidative stress. It is an established biomarker of lung cancer risk. It is commonly measured by enzyme-linked immunosorbent assay. Given its importance, we developed a stable isotope dilution UPLC-tandem mass spectrometric method for the rapid determination of 8-isoprostane in blood.

METHODS

We tested the discriminatory capability of the method in 49 lung cancer patients, 55 benign lung nodule patients detected by chest X-ray, and 41 patients with chronic obstructive pulmonary disease (COPD) or asthma.

RESULTS

Significant differences were found in mean 8-isoprostane levels between the three groups (p = 0.027), and post-hoc tests found higher levels in the lung cancer patients than in patients with benign nodules (p = 0.032) and COPD/asthma (p = 0.014). The receiving operating characteristic area under the curve (AUC) was 0.69 for differentiating the lung cancer group from the benign nodule group, and 0.7 for differentiating from the COPD/asthma group.

CONCLUSIONS

The UPLC-MS/MS-based method is an efficient analytical tool for measuring 8-isoprostane plasma concentrations. The results suggest exploring its utility as a marker for early lung cancer screening.

Genome-wide interaction analysis identified low-frequency variants with sex disparity in lung cancer risk

Differences by sex in lung cancer incidence and mortality have been reported which cannot be fully explained by sex differences in smoking behavior, implying existence of genetic and molecular basis for sex disparity in lung cancer development. However, the information about sex dimorphism in lung cancer risk is quite limited despite the great success in lung cancer association studies. By adopting a stringent two-stage analysis strategy, we performed a genome-wide gene-sex interaction analysis using genotypes from a lung cancer cohort including ~ 47 000 individuals with European ancestry. Three low-frequency variants (minor allele frequency < 0.05), rs17662871 [odds ratio (OR) = 0.71, P = 4.29×10-8); rs79942605 (OR = 2.17, P = 2.81×10-8) and rs208908 (OR = 0.70, P = 4.54×10-8) were identified with different risk effect of lung cancer between men and women. Further expression quantitative trait loci and functional annotation analysis suggested rs208908 affects lung cancer risk through differential regulation of Coxsackie virus and adenovirus receptor gene expression in lung tissues between men and women. Our study is one of the first studies to provide novel insights about the genetic and molecular basis for sex disparity in lung cancer development.

Genetic variants in CYP2B6 and HSD17B12 associated with risk of squamous cell carcinoma of the head and neck

Polycyclic aromatic hydrocarbons (PAH) and tobacco-specific nitrosamines (TSNA) metabolism-related genes play an important role in the development of cancers. We assessed the associations of genetic variants in genes involved in the metabolism of PAHs and TSNA with risk of squamous cell carcinoma of the head and neck (SCCHN) in European populations using two published genome-wide association study datasets. In the single-locus analysis, we identified two SNPs (rs145533669 and rs35246205) in CYP2B6 to be associated with risk of SCCHN (P = 1.57 × 10-4 and .004, respectively), two SNPs (EPHX1 rs117522494 and CYP2B6 rs145533669) to be associated with risk of oropharyngeal cancer (P = .001 and .004, respectively), and one SNP (rs4359199 in HSD17B12) to be associated with risk of oral cancer (P = .006). A significant interaction effect was found between rs4359199 and drinking status on risks of SCCHN and oropharyngeal cancer (P < .05). eQTL and sQTL analyzes revealed that two SNPs (CYP2B6 rs35246205 and HSD17B12 rs4359199) were correlated with alternative splicing or mRNA expression levels of the corresponding genes in liver cells (P < .05 for both). In silico functional annotation suggested that these two SNPs may regulate mRNA expression by affecting the binding of transcription factors. Results from phenome-wide association studies presented significant associations between these genes and risks of other cancers, smoking behavior and alcohol dependence (P < .05). Thus, our study provided some insight into the underlying genetic mechanism of head and neck cancer, which warrants future functional validation.

Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage.

Epigenetics in drug disposition & drug therapy: symposium report of the 24th North American meeting of the International Society for the Study of Xenobiotics (ISSX)

The 24th North American International Society for the Study of Xenobiotics (ISSX) meeting, held virtually from September 13 to 17, 2021, embraced the theme of "Broadening Our Horizons." This reinforces a key mission of ISSX: striving to share innovative science related to drug discovery and development. Session speakers and the ISSX New Investigators Group, which supports the scientific and professional development of student and early career ISSX members, elected to highlight the scientific content presented during the captivating session titled, "Epigenetics in Drug Disposition & Drug Therapy." The impact genetic variation has on drug response is well established; however, this session underscored the importance of investigating the role of epigenetics in drug disposition and drug discovery. Session speakers, Drs. Ning, McClay, and Lazarus, detailed mechanisms by which epigenetic players including long non-coding RNA (lncRNAs), microRNA (miRNAs), DNA methylation, and histone acetylation can alter the expression of genes involved in pharmacokinetics, pharmacodynamics, and toxicity. Dr. Ning detailed current knowledge about miRNAs and lncRNAs and the mechanisms by which they can affect the expression of drug metabolizing enzymes (DMEs) and nuclear receptors. Dr. Lazarus discussed the potential role of miRNAs on UDP-glucuronosyltransferase (UGT) expression and activity. Dr. McClay provided evidence that aging alters methylation and acetylation of DMEs in the liver, affecting gene expression and activity. These topics, compiled by the symposium organizers, presenters, and the ISSX New Investigators Group, are herein discussed, along with exciting future perspectives for epigenetics in drug disposition and drug discovery research.

Drug-Drug Interaction Between Orally Administered Hydrocodone-Acetaminophen and Inhalation of Cannabis Smoke: A Case Report

Objective

To determine if a 2-day protocol measuring pharmacokinetic and pharmacodynamic characteristics can demonstrate drug-drug interactions when smoked cannabis is added to orally administered hydrocodone/acetaminophen combination products.

Case Summary

A 51-year-old non-Hispanic white male with chronic pain diagnoses participated in a 2-day pilot protocol. The participant attended two 7-hour in-lab days where he received 10 blood draws each day and completed self-administered pain and anxiety surveys. For both days, the participant took his prescribed dose of hydrocodone/acetaminophen (1/2 tablet of 7.5 mg/325 mg combination product) with the addition of 1 smoked pre-rolled marijuana cigarette (labeled as 0.5 g; 22.17% Δ9-tetrahydrocannabinol; 0.12% cannabidiol) on Day 2. Blood specimens were analyzed using mass spectrometry to quantify the difference of plasma hydrocodone levels between Day 1 and Day 2.

Results

Compared to Day 1, lower levels of pain and anxiety were reported during Day 2 with the addition of cannabis to oral hydrocodone/acetaminophen. Day 2 pharmacokinetic analysis also revealed more rapid absorption and overall lower levels of hydrocodone in plasma.

Discussion

Lower hydrocodone plasma levels in Day 2 may indicate cannabis's effect on metabolism and reduce the risk of opioid toxicity. The quicker absorption rate of hydrocodone could explain lower pain and anxiety scores reported on the second day.

Conclusion and Relevance

A 2-day protocol was able to capture differences across time in pharmacokinetic and pharmacodynamic measurements. Larger studies can be designed to better characterize the potential drug-drug interaction of cannabis and opioids.

A Large-Scale Genome-Wide Gene-Gene Interaction Study of Lung Cancer Susceptibility in Europeans With a Trans-Ethnic Validation in Asians

INTRODUCTION

Although genome-wide association studies have been conducted to investigate genetic variation of lung tumorigenesis, little is known about gene-gene (G × G) interactions that may influence the risk of non-small cell lung cancer (NSCLC).

METHODS

Leveraging a total of 445,221 European-descent participants from the International Lung Cancer Consortium OncoArray project, Transdisciplinary Research in Cancer of the Lung and UK Biobank, we performed a large-scale genome-wide G × G interaction study on European NSCLC risk by a series of analyses. First, we used BiForce to evaluate and rank more than 58 billion G × G interactions from 340,958 single-nucleotide polymorphisms (SNPs). Then, the top interactions were further tested by demographically adjusted logistic regression models. Finally, we used the selected interactions to build lung cancer screening models of NSCLC, separately, for never and ever smokers.

RESULTS

With the Bonferroni correction, we identified eight statistically significant pairs of SNPs, which predominantly appeared in the 6p21.32 and 5p15.33 regions (e.g., rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.17, p = 6.57 × 10-13; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.17, p = 2.43 × 10-13; rs2858859HLA-DQA1 and rs9275572HLA-DQA2, ORinteraction = 1.15, p = 2.84 × 10-13; rs2853668TERT and rs62329694CLPTM1L, ORinteraction = 0.73, p = 2.70 × 10-13). Notably, even with much genetic heterogeneity across ethnicities, three pairs of SNPs in the 6p21.32 region identified from the European-ancestry population remained significant among an Asian population from the Nanjing Medical University Global Screening Array project (rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.13, p = 0.008; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.11, p = 5.23 × 10-4; rs3135369BTNL2 and rs9271300HLA-DQA1, ORinteraction = 0.89, p = 0.006). The interaction-empowered polygenetic risk score that integrated classical polygenetic risk score and G × G information score was remarkable in lung cancer risk stratification.

CONCLUSIONS

Important G × G interactions were identified and enriched in the 5p15.33 and 6p21.32 regions, which may enhance lung cancer screening models.

Gene-gene interaction of AhRwith and within the Wntcascade affects susceptibility to lung cancer

BACKGROUND

Aberrant Wnt signalling, regulating cell development and stemness, influences the development of many cancer types. The Aryl hydrocarbon receptor (AhR) mediates tumorigenesis of environmental pollutants. Complex interaction patterns of genes assigned to AhR/Wnt-signalling were recently associated with lung cancer susceptibility.

AIM

To assess the association and predictive ability of AhR/Wnt-genes with lung cancer in cases and controls of European descent.

METHODS

Odds ratios (OR) were estimated for genomic variants assigned to the Wnt agonist and the antagonistic genes DKK2, DKK3, DKK4, FRZB, SFRP4 and Axin2. Logistic regression models with variable selection were trained, validated and tested to predict lung cancer, at which other previously identified SNPs that have been robustly associated with lung cancer risk could also enter the model. Furthermore, decision trees were created to investigate variant × variant interaction. All analyses were performed for overall lung cancer and for subgroups.

RESULTS

No genome-wide significant association of AhR/Wnt-genes with overall lung cancer was observed, but within the subgroups of ever smokers (e.g., maker rs2722278 SFRP4; OR  = 1.20; 95% CI 1.13-1.27; p  = 5.6 × 10-10) and never smokers (e.g., maker rs1133683 Axin2; OR  = 1.27; 95% CI 1.19-1.35; p  = 1.0 × 10-12). Although predictability is poor, AhR/Wnt-variants are unexpectedly overrepresented in optimized prediction scores for overall lung cancer and for small cell lung cancer. Remarkably, the score for never-smokers contained solely two AhR/Wnt-variants. The optimal decision tree for never smokers consists of 7 AhR/Wnt-variants and only two lung cancer variants.

CONCLUSIONS

The role of variants belonging to Wnt/AhR-pathways in lung cancer susceptibility may be underrated in main-effects association analysis. Complex interaction patterns in individuals of European descent have moderate predictive capacity for lung cancer or subgroups thereof, especially in never smokers.

Iam hiQ-a novel pair of accuracy indices for imputed genotypes

BACKGROUND

Imputation of untyped markers is a standard tool in genome-wide association studies to close the gap between directly genotyped and other known DNA variants. However, high accuracy with which genotypes are imputed is fundamental. Several accuracy measures have been proposed and some are implemented in imputation software, unfortunately diversely across platforms. In the present paper, we introduce Iam hiQ, an independent pair of accuracy measures that can be applied to dosage files, the output of all imputation software. Iam (imputation accuracy measure) quantifies the average amount of individual-specific versus population-specific genotype information in a linear manner. hiQ (heterogeneity in quantities of dosages) addresses the inter-individual heterogeneity between dosages of a marker across the sample at hand.

RESULTS

Applying both measures to a large case-control sample of the International Lung Cancer Consortium (ILCCO), comprising 27,065 individuals, we found meaningful thresholds for Iam and hiQ suitable to classify markers of poor accuracy. We demonstrate how Manhattan-like plots and moving averages of Iam and hiQ can be useful to identify regions enriched with less accurate imputed markers, whereas these regions would by missed when applying the accuracy measure info (implemented in IMPUTE2).

CONCLUSION

We recommend using Iam hiQ additional to other accuracy scores for variant filtering before stepping into the analysis of imputed GWAS data.

Inhibition of UDP-Glucuronosyltransferase Enzymes by Major Cannabinoids and Their Metabolites

The UDP-glucuronosyltransferase (UGT) family of enzymes play a central role in the metabolism and detoxification of a wide range of endogenous and exogenous compounds. UGTs exhibit a high degree of structural similarity and display overlapping substrate specificity, often making estimations of potential drug-drug interactions difficult to fully elucidate. One such interaction yet to be examined may be occurring between UGTs and cannabinoids, as the legalization of recreational and medicinal cannabis and subsequent co-usage of cannabis and therapeutic drugs increases in the United States and internationally. In the present study, the inhibition potential of the major cannabinoids Δ9-tetrahydrocannabinol (THC), cannabidiol (CBD), and cannabinol (CBN), as well as their major metabolites, was determined in microsomes isolated from HEK293 cells overexpressing individual recombinant UGTs and in microsomes from human liver and kidney specimens. The highest inhibition was seen by CBD against the glucuronidation activity of UGTs 1A9, 2B4, 1A6, and 2B7, with binding-corrected IC50 values of 0.12 ± 0.020 µM, 0.22 ± 0.045 µM, 0.40 ± 0.10 µM, and 0.82 ± 0.15 µM, respectively. Strong inhibition of UGT1A9 was also demonstrated by THC and CBN, with binding-corrected IC50 values of 0.45 ± 0.12 μM and 0.51 ± 0.063 μM, respectively. Strong inhibition of UGT2B7 was also observed for THC and CBN; no or weak inhibition was observed with cannabinoid metabolites. This inhibition of UGT activity suggests that in addition to playing an important role in drug-drug interactions, cannabinoid exposure may have important implications in patients with impaired hepatic or kidney function. SIGNIFICANCE STATEMENT: Major cannabinoids found in the plasma of cannabis users inhibit several UDP-glucuronosyltransferase (UGT) enzymes, including UGT1A6, UGT1A9, UGT2B4, and UGT2B7. This study is the first to show the potential of cannabinoids and their metabolites to inhibit all the major kidney UGTs as well as the two most abundant UGTs present in liver. This study suggests that as all three major kidney UGTs are inhibited by cannabinoids, greater drug-drug interaction effects might be observed from co-use of cannabinods and therapeutics that are cleared renally.

Characterization of Cytosolic Glutathione S-Transferases Involved in the Metabolism of the Aromatase Inhibitor, Exemestane

Exemestane (EXE) is a hormonal therapy used to treat estrogen receptor-positive breast cancer by inhibiting the final step of estrogen biosynthesis catalyzed by the enzyme aromatase. Cysteine conjugates of EXE and its active metabolite 17β-dihydro-EXE (DHE) are the major metabolites found in both the urine and plasma of patients taking EXE. The initial step in cysteine conjugate formation is glutathione conjugation catalyzed by the glutathione S-transferase (GST) family of enzymes. The goal of the present study was to identify cytosolic hepatic GSTs active in the GST-mediated metabolism of EXE and 17β-DHE. Twelve recombinant cytosolic hepatic GSTs were screened for their activity against EXE and 17β-DHE, and glutathionylated EXE and 17β-DHE conjugates were detected by ultra-performance liquid chromatography tandem mass spectrometry. GST α (GSTA) isoform 1, GST μ (GSTM) isoform 3 and isoform 1 were active against EXE, whereas only GSTA1 exhibited activity against 17β-DHE. GSTM1 exhibited the highest affinity against EXE with a Michaelis-Menten constant (KM) value that was 3.8- and 7.1-fold lower than that observed for GSTA1 and GSTM3, respectively. Of the three GSTs, GSTM3 exhibited the highest intrinsic clearance against EXE (intrinsic clearance = 0.14 nl·min-1·mg-1). The KM values observed for human liver cytosol against EXE (46 μM) and 17β-DHE (77 μM) were similar to those observed for recombinant GSTA1 (53 and 30 μM, respectively). Western blot analysis revealed that GSTA1 and GSTM1 composed 4.3% and 0.57%, respectively, of total protein in human liver cytosol; GSTM3 was not detected. These data suggest that GSTA1 is the major hepatic cytosolic enzyme involved in the clearance of EXE and its major active metabolite, 17β-DHE. SIGNIFICANCE STATEMENT: Most previous studies related to the metabolism of the aromatase inhibitor exemestane (EXE) have focused mainly on phase I metabolic pathways and the glucuronidation phase II metabolic pathway. However, recent studies have indicated that glutathionylation is the major metabolic pathway for EXE. The present study is the first to characterize hepatic glutathione S-transferase (GST) activity against EXE and 17β-dihydro-EXE and to identify GST α 1 and GST μ 1 as the major cytosolic GSTs involved in the hepatic metabolism of EXE.

Cannabinoid Metabolites as Inhibitors of Major Hepatic CYP450 Enzymes, with Implications for Cannabis-Drug Interactions

The legalization of cannabis in many parts of the United States and other countries has led to a need for a more comprehensive understanding of cannabis constituents and their potential for drug-drug interactions. Although (-)-trans-Δ9-tetrahydrocannabinol (THC), cannabidiol (CBD), and cannabinol (CBN) are the most abundant cannabinoids present in cannabis, THC metabolites are found in plasma at higher concentrations and for a longer duration than that of the parent cannabinoids. To understand the potential for drug-drug interactions, the inhibition potential of major cannabinoids and their metabolites on major hepatic cytochrome P450 (P450) enzymes was examined. In vitro assays with P450-overexpressing cell microsomes demonstrated that the major THC metabolites 11-hydroxy-Δ9-tetra-hydrocannabinol and 11-nor-9-carboxy-Δ9-THC-glucuronide competitively inhibited several major P450 enzymes, including CYP2B6, CYP2C9, and CYP2D6 (apparent Ki,u values = 0.086 ± 0.066 µM and 0.90 ± 0.54 µM, 0.057 ± 0.044 µM and 2.1 ± 0.81 µM, 0.15 ± 0.067 µM and 2.3 ± 0.54 µM, respectively). 11-Nor-9-carboxy-Δ9- tetrahydrocannabinol exhibited no inhibitory activity against any CYP450 tested. THC competitively inhibited CYP1A2, CYP2B6, CYP2C9, and CYP2D6; CBD competitively inhibited CYP3A4, CYP2B6, CYP2C9, CYP2D6, and CYP2E1; and CBN competitively inhibited CYP2B6, CYP2C9, and CYP2E1. THC and CBD showed mixed-type inhibition for CYP2C19 and CYP1A2, respectively. These data suggest that cannabinoids and major THC metabolites are able to inhibit the activities of multiple P450 enzymes, and basic static modeling of these data suggest the possibility of pharmacokinetic interactions between these cannabinoids and xenobiotics extensively metabolized by CYP2B6, CYP2C9, and CYP2D6. SIGNIFICANCE STATEMENT: Major cannabinoids and their metabolites found in the plasma of cannabis users inhibit several P450 enzymes, including CYP2B6, CYP2C9, and CYP2D6. This study is the first to show the inhibition potential of the most abundant plasma cannabinoid metabolite, THC-COO-Gluc, and suggests that circulating metabolites of cannabinoids play an essential role in CYP450 enzyme inhibition as well as drug-drug interactions.

Polymorphic tandem DNA repeats activate the human telomerase reverse transcriptase gene

Multiple independent sequence variants of the hTERT locus have been associated with telomere length and cancer risks in genome-wide association studies. Here, we identified an intronic variable number tandem repeat, VNTR2-1, as an enhancer-like element, which activated hTERT transcription in a cell in a chromatin-dependent manner. VNTR2-1, consisting of 42-bp repeats with an array of enhancer boxes, cooperated with the proximal promoter in the regulation of hTERT transcription by basic helix-loop-helix transcription factors and maintained hTERT expression during embryonic stem-cell differentiation. Genomic deletion of VNTR2-1 in MelJuSo melanoma cells markedly reduced hTERT transcription, leading to telomere shortening, cellular senescence, and impairment of xenograft tumor growth. Interestingly, VNTR2-1 lengths varied widely in human populations; hTERT alleles with shorter VNTR2-1 were underrepresented in African American centenarians, indicating its role in human aging. Therefore, this polymorphic element is likely a missing link in the telomerase regulatory network and a molecular basis for genetic diversities of telomere homeostasis and age-related disease susceptibilities.

Comprehensive functional annotation of susceptibility variants identifies genetic heterogeneity between lung adenocarcinoma and squamous cell carcinoma

Although genome-wide association studies have identified more than eighty genetic variants associated with non-small cell lung cancer (NSCLC) risk, biological mechanisms of these variants remain largely unknown. By integrating a large-scale genotype data of 15 581 lung adenocarcinoma (AD) cases, 8350 squamous cell carcinoma (SqCC) cases, and 27 355 controls, as well as multiple transcriptome and epigenomic databases, we conducted histology-specific meta-analyses and functional annotations of both reported and novel susceptibility variants. We identified 3064 credible risk variants for NSCLC, which were overrepresented in enhancer-like and promoter-like histone modification peaks as well as DNase I hypersensitive sites. Transcription factor enrichment analysis revealed that USF1 was AD-specific while CREB1 was SqCC-specific. Functional annotation and gene-based analysis implicated 894 target genes, including 274 specifics for AD and 123 for SqCC, which were overrepresented in somatic driver genes (ER = 1.95, P = 0.005). Pathway enrichment analysis and Gene-Set Enrichment Analysis revealed that AD genes were primarily involved in immune-related pathways, while SqCC genes were homologous recombination deficiency related. Our results illustrate the molecular basis of both well-studied and new susceptibility loci of NSCLC, providing not only novel insights into the genetic heterogeneity between AD and SqCC but also a set of plausible gene targets for post-GWAS functional experiments.

Assessing Lung Cancer Absolute Risk Trajectory Based on a Polygenic Risk Model

Lung cancer is the leading cause of cancer-related death globally. An improved risk stratification strategy can increase efficiency of low-dose CT (LDCT) screening. Here we assessed whether individual's genetic background has clinical utility for risk stratification in the context of LDCT screening. On the basis of 13,119 patients with lung cancer and 10,008 controls with European ancestry in the International Lung Cancer Consortium, we constructed a polygenic risk score (PRS) via 10-fold cross-validation with regularized penalized regression. The performance of risk model integrating PRS, including calibration and ability to discriminate, was assessed using UK Biobank data (N = 335,931). Absolute risk was estimated on the basis of age-specific lung cancer incidence and all-cause mortality as competing risk. To evaluate its potential clinical utility, the PRS distribution was simulated in the National Lung Screening Trial (N = 50,772 participants). The lung cancer ORs for individuals at the top decile of the PRS distribution versus those at bottom 10% was 2.39 [95% confidence interval (CI) = 1.92-3.00; P = 1.80 × 10-14] in the validation set (P trend = 5.26 × 10-20). The OR per SD of PRS increase was 1.26 (95% CI = 1.20-1.32; P = 9.69 × 10-23) for overall lung cancer risk in the validation set. When considering absolute risks, individuals at different PRS deciles showed differential trajectories of 5-year and cumulative absolute risk. The age reaching the LDCT screening recommendation threshold can vary by 4 to 8 years, depending on the individual's genetic background, smoking status, and family history. Collectively, these results suggest that individual's genetic background may inform the optimal lung cancer LDCT screening strategy. SIGNIFICANCE: Three large-scale datasets reveal that, after accounting for risk factors, an individual's genetics can affect their lung cancer risk trajectory, thus may inform the optimal timing for LDCT screening.

Causal relationships between body mass index, smoking and lung cancer: Univariable and multivariable Mendelian randomization

At the time of cancer diagnosis, body mass index (BMI) is inversely correlated with lung cancer risk, which may reflect reverse causality and confounding due to smoking behavior. We used two-sample univariable and multivariable Mendelian randomization (MR) to estimate causal relationships of BMI and smoking behaviors on lung cancer and histological subtypes based on an aggregated genome-wide association studies (GWASs) analysis of lung cancer in 29 266 cases and 56 450 controls. We observed a positive causal effect for high BMI on occurrence of small-cell lung cancer (odds ratio (OR) = 1.60, 95% confidence interval (CI) = 1.24-2.06, P = 2.70 × 10-4 ). After adjustment of smoking behaviors using multivariable Mendelian randomization (MVMR), a direct causal effect on small cell lung cancer (ORMVMR = 1.28, 95% CI = 1.06-1.55, PMVMR = .011), and an inverse effect on lung adenocarcinoma (ORMVMR = 0.86, 95% CI = 0.77-0.96, PMVMR = .008) were observed. A weak increased risk of lung squamous cell carcinoma was observed for higher BMI in univariable Mendelian randomization (UVMR) analysis (ORUVMR = 1.19, 95% CI = 1.01-1.40, PUVMR = .036), but this effect disappeared after adjustment of smoking (ORMVMR = 1.02, 95% CI = 0.90-1.16, PMVMR = .746). These results highlight the histology-specific impact of BMI on lung carcinogenesis and imply mediator role of smoking behaviors in the association between BMI and lung cancer.

Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers

Squamous cell carcinomas (SqCC) of the aerodigestive tract have similar etiological risk factors. Although genetic risk variants for individual cancers have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. To identify novel and pleotropic SqCC risk variants, we performed a meta-analysis of GWAS data on lung SqCC (LuSqCC), oro/pharyngeal SqCC (OSqCC), laryngeal SqCC (LaSqCC) and esophageal SqCC (ESqCC) cancers, totaling 13,887 cases and 61,961 controls of European ancestry. We identified one novel genome-wide significant (Pmeta<5x10-8) aerodigestive SqCC susceptibility loci in the 2q33.1 region (rs56321285, TMEM273). Additionally, three previously unknown loci reached suggestive significance (Pmeta<5x10-7): 1q32.1 (rs12133735, near MDM4), 5q31.2 (rs13181561, TMEM173) and 19p13.11 (rs61494113, ABHD8). Multiple previously identified loci for aerodigestive SqCC also showed evidence of pleiotropy in at least another SqCC site, these include: 4q23 (ADH1B), 6p21.33 (STK19), 6p21.32 (HLA-DQB1), 9p21.33 (CDKN2B-AS1) and 13q13.1(BRCA2). Gene-based association and gene set enrichment identified a set of 48 SqCC-related genes rel to DNA damage and epigenetic regulation pathways. Our study highlights the importance of cross-cancer analyses to identify pleiotropic risk loci of histology-related cancers arising at distinct anatomical sites.

Assessment of Regular and NPH Insulin Concentration via Two Methods of Quantification: The Washington State Insulin Concentration Study (WICS)

BACKGROUND

Recent reports have suggested that insulin vials purchased in community pharmacies do not meet the minimum required intact insulin concentration (≥95 U/mL) as defined by the United States Pharmacopeia. We sought to independently obtain multidose human insulin vials from a variety of community pharmacies across the state of Washington and quantitatively measure intact insulin.

METHODS

Sixty 10-mL vials of insulin (n = 30 regular human insulin and n = 30 neutral protamine Hagedorn insulin) were purchased and assayed. To ensure random selection of lots and supply chain sources, insulin samples were purchased on a variety of calendar dates from various pharmacy locations across Washington State, inclusive of both chain and independent pharmacies. All samples were assessed for intact insulin concentration via both Ultra Performance Liquid Chromatography coupled with UV detection (UPLC-UV) and Ultra Performance Liquid Chromatography-Mass Spectrometry (UPLC-MS).

RESULTS

When considering all samples (N = 60), the mean concentration was 101.8 ± 4.4and 91.5 ± 1.9 U/mL as determined by UPLC-UV and UPLC-MS, respectively. Measured concentrations ranged from 90.0 to 108.4 U/mL when assayed by UV UPLC and 86.1 to 95.4 U/mL for UPLC-MS.

CONCLUSION

To our knowledge, this is the first study following the report by Carter et al that assessed human insulin concentrations by both UPLC-UV and UPLC-MS. These findings are important because they demonstrate that the results obtained from these two methods differ and that the method used must be considered when interpreting findings.

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.

Integration of multiomic annotation data to prioritize and characterize inflammation and immune-related risk variants in squamous cell lung cancer

Clinical trial results have recently demonstrated that inhibiting inflammation by targeting the interleukin-1β pathway can offer a significant reduction in lung cancer incidence and mortality, highlighting a pressing and unmet need to understand the benefits of inflammation-focused lung cancer therapies at the genetic level. While numerous genome-wide association studies (GWAS) have explored the genetic etiology of lung cancer, there remains a large gap between the type of information that may be gleaned from an association study and the depth of understanding necessary to explain and drive translational findings. Thus, in this study we jointly model and integrate extensive multiomics data sources, utilizing a total of 40 genome-wide functional annotations that augment previously published results from the International Lung Cancer Consortium (ILCCO) GWAS, to prioritize and characterize single nucleotide polymorphisms (SNPs) that increase risk of squamous cell lung cancer through the inflammatory and immune responses. Our work bridges the gap between correlative analysis and translational follow-up research, refining GWAS association measures in an interpretable and systematic manner. In particular, reanalysis of the ILCCO data highlights the impact of highly associated SNPs from nuclear factor-κB signaling pathway genes as well as major histocompatibility complex mediated variation in immune responses. One consequence of prioritizing likely functional SNPs is the pruning of variants that might be selected for follow-up work by over an order of magnitude, from potentially tens of thousands to hundreds. The strategies we introduce provide informative and interpretable approaches for incorporating extensive genome-wide annotation data in analysis of genetic association studies.

Altered Metabolism of Polycyclic Aromatic Hydrocarbons by UDP-Glycosyltransferase 3A2 Missense Variants

The UDP-glycosyltransferase (UGT) family of enzymes are important in the metabolism of a variety of exogenous substances including polycyclic aromatic hydrocarbons (PAHs), a potent class of environmental carcinogens. As compared to the majority of UGT enzymes, which utilize UDP-glucuronic acid as a cosubstrate, UGT3A2 utilizes alternative cosubstrates (UDP-glucose and UDP-xylose). UGT3A2 is expressed in aerodigestive tract tissues and was highly active against multiple PAHs with both cosubstrates. The goal of the present study was to assess the functional effects of UGT3A2 missense variants (MAF ≥ 0.005) on PAH metabolism and the utilization of cosubstrates. The glycosylation activity (V_max/K_m) of all variants against simple PAHs using both cosubstrates was significantly (P < 0.05) decreased by 42-100% when compared to wild-type UGT3A2. When utilizing UDP-glucose, the variant isoforms exhibited up to a 362-fold decrease in V_max/K_m when compared to wild-type UGT3A2, with a 3.1- to 14-fold decrease for D140N, A344T, and S435Y, a 24- and 43-fold decrease for A436T and R445C, respectively, and a 147- and 362-fold decrease for Y474C and Y74N, respectively. When utilizing UDP-xylose, the variants exhibited up to a 4.0-fold decrease in V_max/K_m when compared to wild-type UGT3A2; Y74N did not exhibit activity, and Y474C did not reach saturation (K_m > 4000 μM). Additionally, both wild-type and variant UGT3A2 exhibited a significant (P < 0.05) difference in their utilization of UDP-glucose vs UDP-xylose as cosubstrates using 1-OH-pyrene as substrate. These data suggest that UGT3A2 missense variants decrease the detoxification of PAHs, potentially resulting in altered individual risk for PAH-related cancers.

Alcohol drinking and head and neck cancer risk: the joint effect of intensity and duration

BACKGROUND

Alcohol is a well-established risk factor for head and neck cancer (HNC). This study aims to explore the effect of alcohol intensity and duration, as joint continuous exposures, on HNC risk.

METHODS

Data from 26 case-control studies in the INHANCE Consortium were used, including never and current drinkers who drunk ≤10 drinks/day for ≤54 years (24234 controls, 4085 oral cavity, 3359 oropharyngeal, 983 hypopharyngeal and 3340 laryngeal cancers). The dose-response relationship between the risk and the joint exposure to drinking intensity and duration was investigated through bivariate regression spline models, adjusting for potential confounders, including tobacco smoking.

RESULTS

For all subsites, cancer risk steeply increased with increasing drinks/day, with no appreciable threshold effect at lower intensities. For each intensity level, the risk of oral cavity, hypopharyngeal and laryngeal cancers did not vary according to years of drinking, suggesting no effect of duration. For oropharyngeal cancer, the risk increased with durations up to 28 years, flattening thereafter. The risk peaked at the higher levels of intensity and duration for all subsites (odds ratio = 7.95 for oral cavity, 12.86 for oropharynx, 24.96 for hypopharynx and 6.60 for larynx).

CONCLUSIONS

Present results further encourage the reduction of alcohol intensity to mitigate HNC risk.

Variation in the UGT2B17 genotype, exemestane metabolism and menopause-related toxicities in the CCTG MAP.3 trial

PURPOSE

To examine associations between the UGT2B17 gene deletion and exemestane metabolites, and commonly reported side effects (fatigue, hot flashes, and joint pain) among postmenopausal women participating in the MAP.3 chemoprevention trial.

METHODS

The analytical samples for the UGT2B17 analysis comprised 1752 women on exemestane and 1721 women on placebo; the exemestane metabolite analysis included 1360 women on exemestane with one-year serum samples. Both the UGT2B17 gene deletion and metabolites were measured in blood. The metabolites were conceptualized as a ratio (17-DHE-Gluc:17-DHE). Symptoms were assessed using the CTCAE v4.0 at approximately 1-year intervals. Log-binomial regression was used to examine the associations between UGT2B17 deletion, exemestane metabolites and each side effect at 1 and up to 5-year follow-up, adjusting for potential confounders.

RESULTS

Among individuals on exemestane with the UGT2B17 gene deletion (i.e., lower detoxification), a higher risk of severe fatigue (RR = 2.59 95% CI: 1.14-5.89) was observed at up to 5-year follow-up. Among individuals on placebo, those with the UGT2B17 gene deletion had a higher risk of any fatigue (RR = 1.39, 95% CI: 1.02-1.89) at year 1. A lower metabolite ratio (poor detoxification) was associated with a higher risk of any fatigue, hot flashes and joint pain at year 1 (fatigue: RR = 1.89, 95% CI: 1.16-3.09; hot flashes: RR = 1.77, 95% CI: 1.40-2.24; joint pain: RR = 2.05, 95% CI: 1.35-3.12); similar associations were observed at 5-year follow-up.

CONCLUSION

Variation in the metabolism of exemestane through the UGT2B17-mediated pathway is associated with subsequent risk of commonly reported symptoms in MAP.3.

Association Analysis of Driver Gene-Related Genetic Variants Identified Novel Lung Cancer Susceptibility Loci with 20,871 Lung Cancer Cases and 15,971 Controls

BACKGROUND

A substantial proportion of cancer driver genes (CDG) are also cancer predisposition genes. However, the associations between genetic variants in lung CDGs and the susceptibility to lung cancer have rarely been investigated.

METHODS

We selected expression-related single-nucleotide polymorphisms (eSNP) and nonsynonymous variants of lung CDGs, and tested their associations with lung cancer risk in two large-scale genome-wide association studies (20,871 cases and 15,971 controls of European descent). Conditional and joint association analysis was performed to identify independent risk variants. The associations of independent risk variants with somatic alterations in lung CDGs or recurrently altered pathways were investigated using data from The Cancer Genome Atlas (TCGA) project.

RESULTS

We identified seven independent SNPs in five lung CDGs that were consistently associated with lung cancer risk in discovery (P < 0.001) and validation (P < 0.05) stages. Among these loci, rs78062588 in TPM3 (1q21.3) was a new lung cancer susceptibility locus (OR = 0.86, P = 1.65 × 10-6). Subgroup analysis by histologic types further identified nine lung CDGs. Analysis of somatic alterations found that in lung adenocarcinomas, rs78062588[C] allele (TPM3 in 1q21.3) was associated with elevated somatic copy number of TPM3 (OR = 1.16, P = 0.02). In lung adenocarcinomas, rs1611182 (HLA-A in 6p22.1) was associated with truncation mutations of the transcriptional misregulation in cancer pathway (OR = 0.66, P = 1.76 × 10-3).

CONCLUSIONS

Genetic variants can regulate functions of lung CDGs and influence lung cancer susceptibility.

IMPACT

Our findings might help unravel biological mechanisms underlying lung cancer susceptibility.

Comparison of tobacco-specific nitrosamine levels in smokeless tobacco products: High levels in products from Bangladesh

Bangladesh exhibits the second highest rate of smokeless tobacco (SLT) product usage in the world, and this has been associated with the high upper aerodigestive tract cancer incidence in this country. The goal of the present study was to examine the levels of the highly carcinogenic tobacco-specific nitrosamines (TSNAs) in Bangladeshi SLT products and compare these levels to that observed in SLT brands from southeast Asia and the USA. The levels of TSNAs and nicotine were determined by LC-MS/MS in twenty-eight SLT brands and several tobacco additives from Bangladesh, as well as several SLT brands from India, Pakistan and the USA. The levels of N-nitrosonornicotine (NNN), 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), N-nitrosoanatabine (NAT) and N-nitrosoanabasine (NAB) in Bangladeshi SLT brands ranged from 1.1-59, 0.15-34, 0.79-45, and 0.037-13 μg/g SLT powder, respectively. The mean levels of the highly carcinogenic TSNAs (NNN+NNK) were 7.4-, 2.4-, and 63-fold higher in Bangladeshi SLT products as compared to SLT brands from the USA, India and Pakistan, respectively; these trends were also observed for NAT and NAB. Similar mean levels of nicotine were observed in the Bangladeshi brands (31 mg/g powder) versus brands from the USA (25 mg/g powder) and India (20 mg/g powder); they were 3-fold higher than brands from Pakistan (10 mg/g powder). Gul SLT brands exhibited the highest pH and the highest levels of unprotonated nicotine. The high levels of TSNAs in Bangladeshi SLT brands may be an important factor contributing to the high rates of upper aerodigestive tract cancer in Bangladesh.

Genome-wide association study of INDELs identified four novel susceptibility loci associated with lung cancer risk

Genome-wide association studies (GWAS) have identified 45 susceptibility loci associated with lung cancer. Only less than SNPs, small insertions and deletions (INDELs) are the second most abundant genetic polymorphisms in the human genome. INDELs are highly associated with multiple human diseases, including lung cancer. However, limited studies with large-scale samples have been available to systematically evaluate the effects of INDELs on lung cancer risk. Here, we performed a large-scale meta-analysis to evaluate INDELs and their risk for lung cancer in 23,202 cases and 19,048 controls. Functional annotations were performed to further explore the potential function of lung cancer risk INDELs. Conditional analysis was used to clarify the relationship between INDELs and SNPs. Four new risk loci were identified in genome-wide INDEL analysis (1p13.2: rs5777156, Insertion, OR = 0.92, p = 9.10 × 10-8 ; 4q28.2: rs58404727, Deletion, OR = 1.19, p = 5.25 × 10-7 ; 12p13.31: rs71450133, Deletion, OR = 1.09, p = 8.83 × 10-7 ; and 14q22.3: rs34057993, Deletion, OR = 0.90, p = 7.64 × 10-8 ). The eQTL analysis and functional annotation suggested that INDELs might affect lung cancer susceptibility by regulating the expression of target genes. After conducting conditional analysis on potential causal SNPs, the INDELs in the new loci were still nominally significant. Our findings indicate that INDELs could be potentially functional genetic variants for lung cancer risk. Further functional experiments are needed to better understand INDEL mechanisms in carcinogenesis.

Protein-altering germline mutations implicate novel genes related to lung cancer development

Few germline mutations are known to affect lung cancer risk. We performed analyses of rare variants from 39,146 individuals of European ancestry and investigated gene expression levels in 7,773 samples. We find a large-effect association with an ATM L2307F (rs56009889) mutation in adenocarcinoma for discovery (adjusted Odds Ratio = 8.82, P = 1.18 × 10-15) and replication (adjusted OR = 2.93, P = 2.22 × 10-3) that is more pronounced in females (adjusted OR = 6.81 and 3.19 and for discovery and replication). We observe an excess loss of heterozygosity in lung tumors among ATM L2307F allele carriers. L2307F is more frequent (4%) among Ashkenazi Jewish populations. We also observe an association in discovery (adjusted OR = 2.61, P = 7.98 × 10-22) and replication datasets (adjusted OR = 1.55, P = 0.06) with a loss-of-function mutation, Q4X (rs150665432) of an uncharacterized gene, KIAA0930. Our findings implicate germline genetic variants in ATM with lung cancer susceptibility and suggest KIAA0930 as a novel candidate gene for lung cancer risk.

Risk Prediction Models for Head and Neck Cancer in the US Population From the INHANCE Consortium

Head and neck cancer (HNC) risk prediction models based on risk factor profiles have not yet been developed. We took advantage of the large database of the International Head and Neck Cancer Epidemiology (INHANCE) Consortium, including 14 US studies from 1981-2010, to develop HNC risk prediction models. Seventy percent of the data were used to develop the risk prediction models; the remaining 30% were used to validate the models. We used competing-risk models to calculate absolute risks. The predictors included age, sex, education, race/ethnicity, alcohol drinking intensity, cigarette smoking duration and intensity, and/or family history of HNC. The 20-year absolute risk of HNC was 7.61% for a 60-year-old woman who smoked more than 20 cigarettes per day for over 20 years, consumed 3 or more alcoholic drinks per day, was a high school graduate, had a family history of HNC, and was non-Hispanic white. The 20-year risk for men with a similar profile was 6.85%. The absolute risks of oropharyngeal and hypopharyngeal cancers were generally lower than those of oral cavity and laryngeal cancers. Statistics for the area under the receiver operating characteristic curve (AUC) were 0.70 or higher, except for oropharyngeal cancer in men. This HNC risk prediction model may be useful in promoting healthier behaviors such as smoking cessation or in aiding persons with a family history of HNC to evaluate their risks.