Tamoxifen (TAM) is an antiestrogen widely used in the treatment and prevention of breast cancer in women. One of the major mechanisms of metabolism of TAM and one of its major active metabolites, 4-hydroxytamoxifen (4-OH-TAM), is via glucuronidation. In the present study, the glucuronidating activities of three common variant isoforms encoded by the human UDP-glucuronosyltransferase (UGT) 1A4 gene were examined against TAM, trans-4-OH-TAM and cis-4-OH-TAM.

Genome-wide association studies (GWAS) identified the chromosome 15q25.1 locus as a leading susceptibility region for lung cancer. However, the pathogenic pathways, through which susceptibility SNPs within chromosome 15q25.1 affects lung cancer risk, have not been explored. We analyzed three cohorts with GWAS data consisting 42,901 individuals and lung expression quantitative trait loci (eQTL) data on 409 individuals to identify and validate the underlying pathways and to investigate the combined effect of genes from the identified susceptibility pathways. The KEGG neuroactive ligand receptor interaction pathway, two Reactome pathways, and 22 Gene Ontology terms were identified and replicated to be significantly associated with lung cancer risk, with P values less than 0.05 and FDR less than 0.1. Functional annotation of eQTL analysis results showed that the neuroactive ligand receptor interaction pathway and gated channel activity were involved in lung cancer risk. These pathways provide important insights for the etiology of lung cancer.

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) is the most abundant and carcinogenic tobacco-specific nitrosamine in tobacco and tobacco smoke. The major metabolic pathway for NNK is carbonyl reduction to form the (R) and (S) enantiomers of 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanol (NNAL) which, like NNK, is a potent lung carcinogen. The goal of this study was to characterize NNAL enantiomer formation in human lung and identify the enzymes responsible for this activity. While (S)-NNAL was the major enantiomer of NNAL formed in incubations with NNK in lung cytosolic fractions, (R)-NNAL comprised ~60 and ~95% of the total NNAL formed in lung whole cell lysates and microsomes, respectively. In studies examining the role of individual recombinant cytosolic reductase enzymes in lung NNAL enantiomer formation, AKR1C1, AKR1C2, AKR1C3, AKR1C4 and CBR1 all exhibited (S)-NNAL-formation activity. To identify the microsomal enzymes responsible for (R)-NNAL formation, 28 microsomal reductase enzymes were screened for expression by real-time PCR in normal human lung. HSD17β6, HSD17β12, KDSR, NSDHL, RDH10, RDH11 and SDR16C5 were all expressed at levels ≥HSD11β1, the only previously reported microsomal reductase enzyme with NNK-reducing activity, with HSD17β12 the most highly expressed. Of these lung-expressing enzymes, only HSD17β12 exhibited activity against NNK, forming primarily (>95%) (R)-NNAL, a pattern consistent with that observed in lung microsomes. siRNA knock-down of HSD17β12 resulted in significant decreases in (R)-NNAL-formation activity in HEK293 cells. These data suggest that both cytosolic and microsomal enzymes are active against NNK and that HSD17β12 is the major active microsomal reductase that contributes to (R)-NNAL formation in human lung.

The UDP-glucuronosyltransferase (UGT) family of enzymes is important in the metabolic elimination of a variety of endogenous compounds such as bile acids, steroids, and fat-soluble vitamins, as well as exogenous compounds including many pharmaceuticals. The UGT2B subfamily is a major family of UGT enzymes expressed in human liver. The identification of novel mechanisms including post-transcriptional regulation by microRNA (miRNA) contributes to interindividual variability in UGT2B expression and is a crucial component in predicting patient drug response. In the present study, a high-resolution liquid chromatography-tandem mass spectrometry method was employed to measure UGT2B protein levels in a panel of human liver microsomal samples (n = 62). Concurrent in silico analysis identified eight candidate miRNAs as potential regulators of UGT2B enzymes. Comparison of UGT2B protein expression and candidate miRNA levels from human liver samples demonstrated a significant inverse correlation between UGT2B10 and UGT2B15 and one of these candidate miRNAs, miR-485-5p. A near-significant correlation was also observed between UGT2B7 and miR-485-5p expression. In vitro analysis using luciferase-containing vectors suggested an interaction of miR-485-5p within the UGT2B10 3'-untranslated region (UTR), and significant reduction in luciferase activity was also observed for a luciferase vector containing the UGT2B7 3'-UTR; however, none was observed for the UBT2B15 3'-UTR. UGT2B10 and UGT2B7 activities were probed using nicotine and 3'-azido-3'-deoxythymidine, respectively, and significant decreases in glucuronidation activity were observed for both substrates in HuH-7 and Hep3B cells upon overexpression of miR-485-5p mimic. This is the first study demonstrating a regulatory role of miR-485-5p for multiple UGT2B enzymes. SIGNIFICANCE STATEMENT: The purpose of this study was to identify novel epigenetic miRNA regulators of the UGT2B drug-metabolizing enzymes in healthy human liver samples. Our results indicate that miRNA 485-5p is a novel regulator of UGT2B7 and UGT2B10, which play an important role in the metabolism of many commonly prescribed medications, carcinogens, and endogenous compounds. This study identified potential miRNA-UGT2B mRNA interactions using a novel proteomic approach, with in vitro experiments undertaken to validate these interactions.

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

Cellular stress response mechanisms normally function to enhance survival and allow for cells to return to homeostasis following an adverse event. Cancer cells often co-opt these same mechanisms as a means to evade apoptosis and mitigate a state of constant cellular stress. Activating transcription factor 5 (ATF5) is upregulated under diverse stress conditions and is overexpressed in a variety of cancers. It was demonstrated ATF5 is a survival factor in transformed, but not normal cells. However, the regulation of ATF5 is not fully understood. The purpose of the present study was to investigate miRNA regulation at the 3' untranslated region (UTR) of ATF5, with the goal of demonstrating a reversal of the upregulation of ATF5 induced under diverse cellular stress in cancer cells. A multifactorial approach using in silico analysis was employed to identify miRNAs 433-3p, 520b-3p, and 129-5p as potential regulators of ATF5, based on their predicted binding sites over the span of the ATF5 3' UTR. Luciferase reporter assay data validated all three miRNA candidates by demonstrating direct binding to the target ATF5 3'. However, functional studies revealed miR-520b-3p as the sole candidate able to reverse the upregulation of ATF5 protein under diverse cellular stress. Additionally, miR-520b-3p levels were inversely related to ATF5 mRNA under endoplasmic reticulum stress and amino acid deprivation. This is the first evidence that regulation at the 3' UTR is involved in modulating ATF5 levels under cellular stress and suggests an important role for miRNA-520b-3p in the regulation of ATF5.

Exemestane (EXE) is used to treat postmenopausal women diagnosed with estrogen receptor positive (ER+) breast cancer. A major mode of metabolism of EXE and its active metabolite, 17β-dihydroexemestane, is via glutathionylation by glutathione-S-transferase (GST) enzymes. The goal of the present study was to investigate the effects of genetic variation in EXE-metabolizing GST enzymes on overall EXE metabolism. Ex vivo assays examining human liver cytosols from 75 subjects revealed the GSTA1 *B*B genotype was associated with significant decreases in S-(androsta-1,4-diene-3,17-dion-6α-ylmethyl)-L-glutathione (P = 0.034) and S-(androsta-1,4-diene-17β-ol-3-on-6α-ylmethyl)-L-gutathione (P = 0.014) formation. In the plasma of 68 ER+ breast cancer patients treated with EXE, the GSTA1 *B*B genotype was associated with significant decreases in both EXE-cysteine (cys) (29%, P = 0.0056) and 17β-DHE-cys (34%, P = 0.032) as compared with patients with the GSTA1*A*A genotype, with significant decreases in EXE-cys (Ptrend = 0.0067) and 17β-DHE-cys (Ptrend = 0.028) observed in patients with increasing numbers of the GSTA1*B allele. A near-significant (Ptrend = 0.060) trend was also observed for urinary EXE-cys levels from the same patients. In contrast, plasma and urinary 17β-DHE-Gluc levels were significantly increased (36%, P = 0.00097 and 52%, P = 0.0089; respectively) in patients with the GSTA1 *B*B genotype. No significant correlations were observed between the GSTM1 null genotype and EXE metabolite levels. These data suggest that the GSTA1*B allele is associated with interindividual differences in EXE metabolism and may play a role in interindividual variability in overall response to EXE. SIGNIFICANCE STATEMENT: The present study is the first comprehensive pharmacogenomic investigation examining the role of genetic variability in GST enzymes on exemestane metabolism. The GSTA1 *B*B genotype was found to contribute to interindividual differences in the metabolism of EXE both ex vivo and in clinical samples from patients taking EXE for the treatment of ER+ breast cancer. Since GSTA1 is a major hepatic phase II metabolizing enzyme in EXE metabolism, the GSTA1*B allele may be an important biomarker for treatment outcomes and toxicities.

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).

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.

Electronic cigarette (E-cigarettes) emissions present a potentially new hazard to neonates through inhalation, dermal and oral contact. Exposure to nicotine containing E-cigarettes may cause significant systemic absorption in neonates due to the potential for multi-route exposure. Systemic absorption of nicotine and constituents of E-cigarette emissions may adversely impact weight and lung development in the neonate. To address these questions we exposed neonatal mice to E-cigarette emissions and measured systemic cotinine levels and alveolar lung growth.

There is a paucity of data describing the impact of type of beverage (coffee versus energy drink), different rates of consumption and different temperature of beverages on the pharmacokinetic disposition of caffeine. Additionally, there is concern that inordinately high levels of caffeine may result from the rapid consumption of cold energy drinks.

UDP-glucuronosyltransferase (UGT) 2A1 is an important enzyme in the detoxification of polycyclic aromatic hydrocarbons found in cigarette smoke. This enzyme is expressed in aerodigestive tract tissues including lung as both its wild-type and exon 4-deleted splice variant isoforms, with the latter acting as a negative regulator of wild-type UGT2A1 activity. UGT2A1 regulation may also be mediated by microRNA (miRNA). To identify miRNA important in the regulation of UGT2A1, expression analysis in tandem with in silico analysis suggested miR-196a-5p and miR-196b-5p as potential top candidates. Significant reductions in firefly luciferase activity were observed in human embryonic kidney cell line 293 cells cotransfected with the wild-type UGT2A1 3'-untranslated region (UTR)-containing luciferase plasmid and either miR-196a-5p (62%, P = 0.00080) or miR-196b-5p (60%, P = 0.00030) mimics. In pull-down assays, there was a 3.4- and 5.2-fold increase in miR-196a-5p (P = 0.054) and miR-196b-5p (P = 0.035), respectively, using the UGT2A1 3'-UTR biotinylated mRNA probe as compared with the β-actin coding region control mRNA probe. UGT2A1 mRNA was reduced by 25% (P = 0.058) and 35% (P = 0.023) in H146 and H1944 cells, respectively, after overexpression of the miR196a-5p mimic. A similar 32% (P = 0.030) and 41% (P = 0.016) reduction was observed after over-expression of the miR-196b-5p mimic. In H146 cells transfected with miRNA mimic together with a small interfering RNA (siRNA) specific for the UGT2A1 splice variant, a significant reduction in 3-hydroxy-benzo[a]pyrene-glucuronide formation was observed. The miR-196a-5p- and miR-196b-55p-treated cells exhibited reductions of 35% (P = 0.047) and 44% (P = 0.0063), respectively. These data suggest that miR-196a-5p and miR-196b-5p play an important role in UGT2A1 regulation within the lung and potentially other aerodigestive tract tissues.

Nicotine is the major pharmacologically active substance in tobacco. Several studies have examined genotypes related to nicotine metabolism, but few studies have been performed in the Mexican population. The objective was to identify associations between gene variants in metabolizing enzymes and the urinary levels of nicotine metabolites among Mexican smokers. The levels of nicotine and its metabolites were determined in the urine of 88 young smokers from Mexico, and 167 variants in 24 genes associated with nicotine metabolism were genotyped by next-generation sequencing (NGS). Trans-3'-hydroxy-cotinine (3HC) and 4-hydroxy-4-(3-pyridyl)-butanoic acid were the most abundant metabolites (35 and 17%, respectively). CYP2A6*12 was associated with 3HC (p = 0.014). The rs145014075 was associated with creatinine-adjusted levels of nicotine (p = 0.035), while the rs12471326 (UGT1A9) was associated to cotinine-N-glucuronide (p = 0.030). CYP2A6 and UGT1A9 variants are associated to nicotine metabolism. 4HPBA metabolite was an abundant urinary metabolite in young Mexican smokers.

Polycyclic aromatic hydrocarbons (PAHs) are potent carcinogens and are a primary risk factor for the development of lung and other aerodigestive tract cancers in smokers. The detoxification of PAHs by glucuronidation is well-characterized for the UDP-glycosyltransferase (UGT) 1A, 2A, and 2B subfamilies; however, the role of the UGT3A subfamily in PAH metabolism remains poorly understood. UGT3A enzymes are functionally distinct from other UGT subfamilies (which use UDP-glucuronic acid as a cosubstrate) due to their utilization of alternative cosubstrates (UDP-N-acetylglucosamine for UGT3A1, and UDP-glucose and UDP-xylose for UGT3A2). The goal of the present study was to characterize UGT3A glycosylation activity against PAHs and examine their expression in human aerodigestive tract tissues. In vitro metabolism assays using UGT3A2-overexpressing cell microsomes indicated that UGT3A2 exhibits glycosylation activity against all of the simple and complex PAHs tested. The V max/K m ratios for UGT3A2 activity with UDP-xylose versus UDP-glucose as the cosubstrate ranged from 0.65 to 4.4 for all PAHs tested, demonstrating that PAH glycosylation may be occurring at rates up to 4.4-fold higher with UDP-xylose than with UDP-glucose. Limited glycosylation activity was observed against PAHs with UGT3A1-overexpressing cell microsomes. While UGT3A2 exhibited low levels of hepatic expression, it was shown by western blot analysis to be widely expressed in aerodigestive tract tissues. Conversely, UGT3A1 exhibited the highest expression in liver with lower expression in aerodigestive tract tissues. These data suggest that UGT3A2 plays an important role in the detoxification of PAHs in aerodigestive tract tissues, and that there may be cosubstrate-dependent differences in the detoxification of PAHs by UGT3A2. SIGNIFICANCE STATEMENT: UGT3A2 is highly active against PAHs with either UDP-glucose or UDP-xylose as a cosubstrate. UGT3A1 exhibited low levels of activity against PAHs. UGT3A1 is highly expressed in liver while UGT3A2 is well expressed in extrahepatic tissues. UGT3A2 may be an important detoxifier of PAHs in humans.

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.

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.

Exemestane (EXE) is an aromatase inhibitor used to treat hormone-dependent breast cancer. EXE is extensively metabolized, with unchanged EXE and its active metabolite 17-dihydroexemestane (17-DHE) accounting for 17 and 12%, respectively, of total plasma EXE in vivo The major circulating EXE metabolites are the cysteine conjugates of EXE and 17-DHE, and the 17-DHE glucuronide, which together account for 70% of total plasma EXE in vivo The goal of the present study was to examine the inhibition potential of major metabolites of EXE through inhibition assays using aromatase-overexpressing cells and pooled ovarian tissues. Estrone formation was used as a measure of aromatase activity and was detected and quantified using UPLC-MS. EXE-cys, 17β-DHE, and 17β-DHE-cys all exhibited inhibition of estrone formation at both 1 µM and 10 µM concentrations, with 17β-DHE and EXE-cys showing significant inhibition of estrone formation (63% each) at 10 µM. In contrast, 17β-DHE-Gluc displayed minimal inhibition (5-8%) at both concentrations. In ovarian tissue, EXE-cys and 17β-DHE showed similar patterns of inhibition, with 49% and 47% inhibition, respectively, at 10 µM. The IC50 value for EXE-cys (16 {plus minus} 10 µM) was similar to 17β-DHE (9.2 {plus minus} 2.7 µM) and higher than EXE (1.3 {plus minus} 0.28 µM), and all three compounds showed time-dependent inhibition with IC50 shifts of 13 {plus minus} 10, 5.0 {plus minus} 2.5 and 36 {plus minus} 12-fold, respectively. Given its high circulating levels in patients taking EXE, these results suggest that EXE-cys may contribute to the pharmacologic effect of EXE in vivo Significance Statement The current study is the first to examine the major phase II metabolites of EXE (EXE-cys, 17β-DHE-cys, and 17β-DHE-Gluc) for inhibition potential against the target enzyme, aromatase (CYP19A1). EXE-cys was found to significantly inhibit aromatase in a time dependent manner. Given its high circulating levels in patients taking EXE, this phase II metabolite may play an important role in reducing circulating estrogen levels in vivo.

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.

Lung cancer has several genetic associations identified within the major histocompatibility complex (MHC); although the basis for these associations remains elusive. Here, we analyze MHC genetic variation among 26,044 lung cancer patients and 20,836 controls densely genotyped across the MHC, using the Illumina Illumina OncoArray or Illumina 660W SNP microarray. We impute sequence variation in classical HLA genes, fine-map MHC associations for lung cancer risk with major histologies and compare results between ethnicities. Independent and novel associations within HLA genes are identified in Europeans including amino acids in the HLA-B*0801 peptide binding groove and an independent HLA-DQB1*06 loci group. In Asians, associations are driven by two independent HLA allele sets that both increase risk in HLA-DQB1*0401 and HLA-DRB1*0701; the latter better represented by the amino acid Ala-104. These results implicate several HLA-tumor peptide interactions as the major MHC factor modulating lung cancer susceptibility.

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