Identification and functional characterization of a novel UDP-glucuronosyltransferase 2A1 splice variant: potential importance in tobacco-related cancer susceptibility

Activation of Cryptic Donor Splice Sites within the UDP-Glucuronosyltransferase (UGT)1A First-Exon Region Generates Variant Transcripts That Encode UGT1A Proteins with Truncated Aglycone-Binding Domains

The human UDP-glucuronosyltransferases (UGTs) have crucial roles in metabolizing and clearing numerous small lipophilic compounds. The UGT1A locus generates nine UGT1A mRNAs, 65 spliced transcripts, and 34 circular RNAs. In this study, our analysis of published UGT-RNA capture sequencing (CaptureSeq) datasets identified novel splice junctions that predict 24 variant UGT1A transcripts derived from ligation of exon 2 to unique sequences within the UGT1A first-exon region using cryptic donor splice sites. Of these variants, seven (1A1_n1, 1A3_n3, 1A4_n4, 1A5_n1, 1A8_n2, 1A9_n2, 1A10_n7) are predicted to encode UGT1A proteins with truncated aglycone-binding domains. We assessed their expression profiles and deregulation in cancer using four RNA sequencing (RNA-Seq) datasets of paired normal and cancerous drug-metabolizing tissues from large patient cohorts. Variants were generally coexpressed with their canonical counterparts with a higher relative abundance in tumor than in normal tissues. Variants showed tissue-specific expression with high interindividual variability but overall low abundance. However, 1A8_n2 showed high abundance in normal and cancerous colorectal tissues, with levels that approached or surpassed canonical 1A8 mRNA levels in many samples. We cloned 1A8_n2 and showed expression of the predicted protein (1A8_i3) in human embryonic kidney (HEK)293T cells. Glucuronidation assays with 4-methylumbelliferone (4MU) showed that 1A8_i3 had no activity and was unable to inhibit the activity of 1A8_i1 protein. In summary, the activation of cryptic donor splice sites within the UGT1A first-exon region expands the UGT1A transcriptome and proteome. The 1A8_n2 cryptic donor splice site is highly active in colorectal tissues, representing an important cis-regulatory element that negatively regulates the function of the UGT1A8 gene through pre-mRNA splicing. SIGNIFICANT STATEMENT: The UGT1A locus generates nine canonical mRNAs, 65 alternately spliced transcripts, and 34 different circular RNAs. The present study reports a series of novel UDP-glucuronosyltransferase (UGT)1A variants resulting from use of cryptic donor splice sites in both normal and cancerous tissues, several of which are predicted to encode variant UGT1A proteins with truncated aglycone-binding domains. Of these, 1A8_n2 shows exceptionally high abundance in colorectal tissues, highlighting its potential role in the first-pass metabolism in gut through the glucuronidation pathway.

Integrated clinical genomic analysis reveals xenobiotic metabolic genes are downregulated in meningiomas of current smokers

INTRODUCTION

Meningiomas are the most common primary intracranial tumor. Recently, various genetic classification systems for meningioma have been described. We sought to identify clinical drivers of different molecular changes in meningioma. As such, clinical and genomic consequences of smoking in patients with meningiomas remain unexplored.

METHODS

88 tumor samples were analyzed in this study. Whole exome sequencing (WES) was used to assess somatic mutation burden. RNA sequencing data was used to identify differentially expressed genes (DEG) and genes sets (GSEA).

RESULTS

Fifty-seven patients had no history of smoking, twenty-two were past smokers, and nine were current smokers. The clinical data showed no major differences in natural history across smoking status. WES revealed absence of AKT1 mutation rate in current or past smokers compared to non-smokers (p = 0.046). Current smokers had increased mutation rate in NOTCH2 compared to past and never smokers (p < 0.05). Mutational signature from current and past smokers showed disrupted DNA mismatch repair (cosine-similarity = 0.759 and 0.783). DEG analysis revealed the xenobiotic metabolic genes UGT2A1 and UGT2A2 were both significantly downregulated in current smokers compared to past (Log2FC = - 3.97, padj = 0.0347 and Log2FC = - 4.18, padj = 0.0304) and never smokers (Log2FC = - 3.86, padj = 0.0235 and Log2FC = - 4.20, padj = 0.0149). GSEA analysis of current smokers showed downregulation of xenobiotic metabolism and enrichment for G2M checkpoint, E2F targets, and mitotic spindle compared to past and never smokers (FDR < 25% each).

CONCLUSION

In this study, we conducted a comparative analysis of meningioma patients based on their smoking history, examining both their clinical trajectories and molecular changes. Meningiomas from current smokers were more likely to harbor NOTCH2 mutations, and AKT1 mutations were absent in current or past smokers. Moreover, both current and past smokers exhibited a mutational signature associated with DNA mismatch repair. Meningiomas from current smokers demonstrate downregulation of xenobiotic metabolic enzymes UGT2A1 and UGT2A2, which are downregulated in other smoking related cancers. Furthermore, current smokers exhibited downregulation xenobiotic metabolic gene sets, as well as enrichment in gene sets related to mitotic spindle, E2F targets, and G2M checkpoint, which are hallmark pathways involved in cell division and DNA replication control. In aggregate, our results demonstrate novel alterations in meningioma molecular biology in response to systemic carcinogens.

Genetic associations with lymphomas in Polish patients: A pooled-DNA genome-wide association analysis

Several single nucleotide polymorphisms (SNPs) associated with susceptibility to Hodgkin lymphoma (HL) and diffuse large B-cell lymphoma (DLBCL) have been identified. The aim of this study was to identify susceptibility loci for HL and DLBCL in Polish patients. Altogether, DLBCL (n = 218 and HL patients (n = 224) and healthy individuals (n = 1181) were recruited. Lymphoma diagnosis was based on standard criteria. Genome-wide association study (GWAS) was performed using pooled-DNA samples on llumina Infinium Omni2.5 Exome-8 v1.3, and selected loci were replicated by TaqMan SNP genotyping of individuals. GWAS detected thirteen and seven SNPs associated with DLBCL and HL, respectively. In the replication study, six and seven SNPs reached significance after correction for multiple testing in the DLBCL and HL cohorts, respectively. One and four SNPs associated with DLBCL and HL, respectively, were localized within, and two SNPs-near the major histocompatibility complex (MHC) region. In conclusion, the majority of loci associated with HL and DLBCL aetiology in previous studies have potential roles in immune function. Our pooled-DNA GWAS enabled the identification of several susceptibility loci for DLBCL and HL in the Polish population; some of them were mapped within or adjacent to the MHC, and other associated SNPs were located outside the MHC.

Regulation of human UDP-glycosyltransferase (UGT) genes by miRNAs

The human UGT gene superfamily is divided into four subfamilies (UGT1, UGT2, UGT3 and UGT8) that encodes 22 functional enzymes. UGTs are critical for the metabolism and clearance of numerous endogenous and exogenous compounds, including steroid hormones, bile acids, bilirubin, fatty acids, carcinogens, and therapeutic drugs. Therefore, the expression and activities of UGTs are tightly regulated by multiple processes at the transcriptional, post-transcriptional and post-translational levels. During recent years, nearly twenty studies have investigated the post-transcriptional regulation of UGT genes by miRNAs using human cancer cell lines (predominantly liver cancer). Overall, 14 of the 22 UGT mRNAs (1A1, 1A3, 1A4, 1A6, 1A8, 1A9, 1A10, 2A1, 2B4, 2B7, 2B10, 2B15, 2B17, UGT8) have been shown to be regulated by various miRNAs through binding to their respective 3' untranslated regions (3'UTRs). Three 3'UTRs (UGT1A, UGT2B7 and UGT2B15) contain the largest number of functional miRNA target sites; in particular, the UGT1A 3'UTR contains binding sites for 12 miRNAs (548d-5p, 183-5p, 214-5p, 486-3p, 200a-3p, 491-3p, 141-3p, 298, 103b, 376b-3p, 21-3p, 1286). Although all nine UGT1A family members have the same 3'UTR, these miRNA target sites appear to be functional in an isoform-specific and cellular context-dependent manner. Collectively, these observations demonstrate that miRNAs represent important post-transcriptional regulators of the UGT gene superfamily. In this article, we present a comprehensive review of reported UGT/miRNA regulation studies, describe polymorphisms within functional miRNA target sites that may affect their functionalities, and discuss potential cooperative and competitive regulation of UGT mRNAs by miRNAs through adjacently located miRNA target sites.

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.

Mixed-Up Sugars: Glycosyltransferase Cross-Reactivity in Cancerous Tissues and Their Therapeutic Targeting

The main categories of glycan changes in cancer are: (1) decreased expression of histo-blood group A and/or B antigens and increased Lewis-related antigens, (2) appearance of cryptic antigens, such as Tn and T, (3) emergence of genetically incompatible glycans, such as A antigen expressed in tumors of individuals of group B or O and heterophilic expression of Forssman antigen (FORS1), and (4) appearance of neoglycans. This review focuses on the expression of genetically incompatible A/B/FORS1 antigens in cancer. Several possible molecular mechanisms are exemplified, including missense mutations that alter the sugar specificity of A and B glycosyltransferases (AT and BT, respectively), restoration of the correct codon reading frame of O alleles, and modification of acceptor specificity of AT to synthesize the FORS1 antigen by missense mutations and/or altered splicing. Taking advantage of pre-existing natural immunity, the potential uses of these glycans for immunotherapeutic targeting will also be discussed.

Circular RNAs of UDP-Glycosyltransferase (UGT) Genes Expand the Complexity and Diversity of the UGT Transcriptome

The human UDP-glycosyltransferase (UGT) gene superfamily generates 22 canonical transcripts coding for functional enzymes and also produces nearly 150 variant UGT transcripts through alternative splicing and intergenic splicing. In the present study, our analysis of circRNA databases identified backsplicing events that predicted 85 circRNAs from UGT genes, with 33, 11, and 19 circRNAs from UGT1A, UGT2B4, UGT8, respectively. Most of these UGT circRNAs were reported by one database and had low abundance in cell- or tissue-specific contexts. Using reverse-transcriptase polymerase chain reaction with divergent primers and cDNA samples from human tissues and cell lines, we found 13 circRNAs from four UGT genes: UGT1A (three), UGT2B7 (one), UGT2B10 (one), and UGT8 (eight). Notably, all eight UGT8 circRNAs contain open reading frames that include the canonical start AUG codon and encode variant proteins that all have the common 274-amino acidN-terminal region of wild-type UGT8 protein. We further showed that one UGT8 circRNA (circ_UGT8-1) was broadly expressed in human tissues and cell lines, resistant to RNase R digestion, and predominately present in the cytoplasm. We cloned five UGT8 circRNAs into the Zinc finger with KRAB and SCAN domains 1 vector and transfected them into HEK293T cells. All these vectors produced both circRNAsand linear transcripts with varying circular/linear ratios (0.17-1.14).Western blotting and mass spectrometry assays revealed that only linear transcripts and not circRNAs were translated. In conclusion, our findings of nearly 100 circRNAs greatly expand the complexity and diversity of the UGT transcriptome; however, UGT circRNAs are expressed at a very low level in specific cellular contexts, and their biologic functions remain to be determined. SIGNIFICANCE STATEMENT: The human UGT gene transcriptome comprises 22 canonical transcripts coding for functional enzymes and approximately 150 alternatively spliced and chimeric variant transcripts. The present study identified nearly 100 circRNAs from UGT genes, thus greatly expanding the complexity and diversity of the UGT transcriptome. UGT circRNAs were expressed broadly in human tissues and cell lines; however, most showed very low abundance in tissue- and cell-specific contexts, and therefore their biological functions remain to be investigated.

The odorant metabolizing enzyme UGT2A1: Immunolocalization and impact of the modulation of its activity on the olfactory response

Odorant metabolizing enzymes (OMEs) are expressed in the olfactory epithelium (OE) where they play a significant role in the peripheral olfactory process by catalyzing the fast biotransformation of odorants leading either to their elimination or to the synthesis of new odorant stimuli. The large family of OMEs gathers different classes which interact with a myriad of odorants alike and complementary to olfactory receptors. Thus, it is necessary to increase our knowledge on OMEs to better understand their function in the physiological process of olfaction. This study focused on a major olfactory UDP-glucuronosyltransferase (UGT): UGT2A1. Immunohistochemistry and immunogold electronic microscopy allowed to localize its expression in the apical part of the sustentacular cells and originally at the plasma membrane of the olfactory cilia of the olfactory sensory neurons, both locations in close vicinity with olfactory receptors. Moreover, using electroolfactogram, we showed that a treatment of the OE with beta-glucuronidase, an enzyme which counterbalance the UGTs activity, increased the response to eugenol which is a strong odorant UGT substrate. Altogether, the results supported the function of the olfactory UGTs in the vertebrate olfactory perireceptor process.

Alternative splicing: An important regulatory mechanism in colorectal carcinoma

Alternative splicing (AS) is a process that produces various mRNA splicing isoforms via different splicing patterns of mRNA precursors (pre-mRNAs). AS is the primary mechanism for increasing the types and quantities of proteins to improve biodiversity and influence multiple biological processes, including chromatin modification, signal transduction, and protein expression. It has been reported that AS is involved in the tumorigenesis and development of colorectal carcinoma (CRC). In this review, we delineate the concept, types, regulatory processes, and technical advances of AS and focus on the role of AS in CRC initiation, progression, treatment, and prognosis. This summary of the current knowledge about AS will contribute to our understanding of CRC initiation and development. This study will help in the discovery of novel biomarkers and therapeutic targets for CRC prognosis and treatment.

The UGTome: The expanding diversity of UDP glycosyltransferases and its impact on small molecule metabolism

The UDP glycosyltransferase (UGT) superfamily of enzymes is responsible for the metabolism and clearance of thousands of lipophilic chemicals including drugs, toxins and endogenous signaling molecules. They provide a protective interface between the organism and its chemical-rich environment, as well as controlling critical signaling pathways to maintain healthy tissue function. UGTs are associated with drug responses and interactions, as well as a wide range of diseases including cancer. The human genome contains 22 UGT genes; however as befitting their exceptionally diverse substrate ranges and biological activities, the output of these UGT genes is functionally diversified by multiple processes including alternative splicing, post-translational modification, homo- and hetero-oligomerization, and interactions with other proteins. All UGT genes are subject to extensive alternative splicing generating variant/truncated UGT proteins with altered functions including the capacity to dominantly modulate/inhibit cognate full-length forms. Heterotypic oligomerization of different UGTs can alter kinetic properties relative to monotypic complexes, and potentially produce novel substrate specificities. Moreover, the recently profiled interactions of UGTs with non-UGT proteins may facilitate coordination between different metabolic processes, as well as providing opportunities for UGTs to engage in novel 'moonlighting' functions. Herein we provide a detailed and comprehensive review of all known modes of UGT functional diversification and propose a UGTome model to describe the resulting expansion of metabolic capacity and its potential to modulate drug/xenobiotic responses and cell behaviours in normal and disease contexts.

Regulation of UGT2A1 by miR-196a-5p and miR-196b-5p

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.

Evaluation of vitamin D biosynthesis and pathway target genes reveals UGT2A1/2 and EGFR polymorphisms associated with epithelial ovarian cancer in African American Women

An association between genetic variants in the vitamin D receptor (VDR) gene and epithelial ovarian cancer (EOC) was previously reported in women of African ancestry (AA). We sought to examine associations between genetic variants in VDR and additional genes from vitamin D biosynthesis and pathway targets (EGFR, UGT1A, UGT2A1/2, UGT2B, CYP3A4/5, CYP2R1, CYP27B1, CYP24A1, CYP11A1, and GC). Genotyping was performed using the custom-designed 533,631 SNP Illumina OncoArray with imputation to the 1,000 Genomes Phase 3 v5 reference set in 755 EOC cases, including 537 high-grade serous (HGSOC), and 1,235 controls. All subjects are of African ancestry (AA). Logistic regression was performed to estimate odds ratios (OR) and 95% confidence intervals (CI). We further evaluated statistical significance of selected SNPs using the Bayesian False Discovery Probability (BFDP). A significant association with EOC was identified in the UGT2A1/2 region for the SNP rs10017134 (per allele OR = 1.4, 95% CI = 1.2-1.7, P = 1.2 × 10-6 , BFDP = 0.02); and an association with HGSOC was identified in the EGFR region for the SNP rs114972508 (per allele OR = 2.3, 95% CI = 1.6-3.4, P = 1.6 × 10-5 , BFDP = 0.29) and in the UGT2A1/2 region again for rs1017134 (per allele OR = 1.4, 95% CI = 1.2-1.7, P = 2.3 × 10-5 , BFDP = 0.23). Genetic variants in the EGFR and UGT2A1/2 may increase susceptibility of EOC in AA women. Future studies to validate these findings are warranted. Alterations in EGFR and UGT2A1/2 could perturb enzyme efficacy, proliferation in ovaries, impact and mark susceptibility to EOC.

Functional and molecular characterization of UDP-glucuronosyltransferase 2 family in cynomolgus macaques

UDP-glucuronosyltransferases (UGTs) are essential enzymes metabolizing endogenous and exogenous chemicals. However, characteristics of UGTs have not been fully investigated in molecular levels of cynomolgus macaques, one of non-human primates widely used in preclinical drug metabolism studies. In this study, three UGT2A cDNAs (UGT2A1, 2A2, and 2A3) were isolated and characterized along with seven UGT2Bs previously identified in cynomolgus macaques. Several transcript variants were found in cynomolgus UGT2A1 and UGT2A2, like human orthologs. Cynomolgus UGT2A and UGT2B amino acid sequences were highly identical (87-96%) to their human counterparts. By phylogenetic analysis, all these cynomolgus UGT2s were more closely clustered with their human homologs than with dog, rat, or mouse UGT2s. Especially, UGT2As showed orthologous relationships between humans and cynomolgus macaques. All the cynomolgus UGT2 mRNAs were expressed in livers, jejunum, and/or kidneys abundantly, except that UGT2A1 and UGT2A2 mRNAs were predominantly expressed in nasal mucosa, like human UGT2s. UGT2A and UGT2B genes together form a gene cluster in the cynomolgus and human genome. Among the seven cynomolgus UGT2Bs heterologously expressed in yeast, UGT2B9 and UGT2B30 showed activities in estradiol 17-O-glucuronidation and morphine 3-O-glucuronidation but did not show activities in estradiol 3-O-glucuronidation, similar to human UGT2Bs. In liver microsomes, cynomolgus macaques showed higher estradiol 17-O-glucuronidase and morphine 3-O-glucuronidase activities than humans, suggesting functional activities of the responsible UGT2B enzymes in cynomolgus macaques. Therefore, cynomolgus UGT2s had overall molecular similarities to human UGT2s, but also showed some differences in UGT2B enzyme properties.

Intergenic Splicing between Four Adjacent UGT Genes (2B15, 2B29P2, 2B17, 2B29P1) Gives Rise to Variant UGT Proteins That Inhibit Glucuronidation via Protein-Protein Interactions

Recent studies have investigated alternative splicing profiles of UDP-glucuronosyltransferase (UGT) genes and identified over 130 different alternatively spliced UGT transcripts. Although UGT genes are highly clustered, the formation of chimeric transcripts by intergenic splicing between two or more UGT genes has not yet been reported. This study identified 12 chimeric transcripts (chimeras A-L) containing exons from two or three genes of the four neighboring UGT genes (UGT2B15, UGT2B29P2, UGT2B17, and UGT2B29P1) in human liver and prostate cancer cells. These chimeras typically contain the first five exons of UGT2B15 or UGT2B17 (exons 1-5) spliced to a terminal exon (exon 6) from a downstream UGT gene. Hence they encode truncated UGTs with novel C-terminal peptides. Functional assays of representative chimeric UGT proteins (termed chimeric UGT2B15 and chimeric UGT2B17) showed that they are inactive and can repress the activity of wild-type UGTs. Coimmunoprecipitation assays demonstrated heterotypic interactions between chimeric UGT2B15 (or chimeric UGT2B17) and the UGT2B7 protein. Thus oligomerization of the chimeric UGTs with wild-type UGTs may explain their inhibitory activity. Studies in breast and prostate cancer cells showed that both wild-type and chimeric UGT2B15 and UGT2B17 transcripts are regulated in a similar way at the transcriptional level by sex hormones through their canonical promoters but are differentially regulated at the post-transcriptional level by micro-RNA 376c via their unique 3'-untranslated regions. In conclusion, the formation of chimeric transcripts by intergenic splicing among UGT genes represents a novel mechanism contributing to the diversity of the human UGT transcriptome and proteome. The differential post-transcriptional regulation of wild-type and variant transcripts by micro-RNAs may contribute to their deregulated expression in cancer.

Post-transcriptional Regulation of UGT2B10 Hepatic Expression and Activity by Alternative Splicing

The detoxification enzyme UDP-glucuronosyltransferase UGT2B10 is specialized in the N-linked glucuronidation of many drugs and xenobiotics. Preferred substrates possess tertiary aliphatic amines and heterocyclic amines, such as tobacco carcinogens and several antidepressants and antipsychotics. We hypothesized that alternative splicing (AS) constitutes a means to regulate steady-state levels of UGT2B10 and enzyme activity. We established the transcriptome of UGT2B10 in normal and tumoral tissues of multiple individuals. The highest expression was in the liver, where 10 AS transcripts represented 50% of the UGT2B10 transcriptome in 50 normal livers and 44 hepatocellular carcinomas. One abundant class of transcripts involves a novel exonic sequence and leads to two alternative (alt.) variants with novel in-frame C termini of 10 or 65 amino acids. Their hepatic expression was highly variable among individuals, correlated with canonical transcript levels, and was 3.5-fold higher in tumors. Evidence for their translation in liver tissues was acquired by mass spectrometry. In cell models, they colocalized with the enzyme and influenced the conjugation of amitriptyline and levomedetomidine by repressing or activating the enzyme (40%-70%; P < 0.01) in a cell context-specific manner. A high turnover rate for the alt. proteins, regulated by the proteasome, was observed in contrast to the more stable UGT2B10 enzyme. Moreover, a drug-induced remodeling of UGT2B10 splicing was demonstrated in the HepaRG hepatic cell model, which favored alt. variants expression over the canonical transcript. Our findings support a significant contribution of AS in the regulation of UGT2B10 expression in the liver with an impact on enzyme activity.

Quantitative profiling of the UGT transcriptome in human drug-metabolizing tissues

Alternative splicing as a mean to control gene expression and diversify function is suspected to considerably influence drug response and clearance. We report the quantitative expression profiles of the human UGT genes including alternatively spliced variants not previously annotated established by deep RNA-sequencing in tissues of pharmacological importance. We reveal a comprehensive quantification of the alternative UGT transcriptome that differ across tissues and among individuals. Alternative transcripts that comprise novel in-frame sequences associated or not with truncations of the 5'- and/or 3'- termini, significantly contribute to the total expression levels of each UGT1 and UGT2 gene averaging 21% in normal tissues, with expression of UGT2 variants surpassing those of UGT1. Quantitative data expose preferential tissue expression patterns and remodeling in favor of alternative variants upon tumorigenesis. These complex alternative splicing programs have the strong potential to contribute to interindividual variability in drug metabolism in addition to diversify the UGT proteome.

Endogenous Protein Interactome of Human UDP-Glucuronosyltransferases Exposed by Untargeted Proteomics

The conjugative metabolism mediated by UDP-glucuronosyltransferase enzymes (UGTs) significantly influences the bioavailability and biological responses of endogenous molecule substrates and xenobiotics including drugs. UGTs participate in the regulation of cellular homeostasis by limiting stress induced by toxic molecules, and by controlling hormonal signaling networks. Glucuronidation is highly regulated at genomic, transcriptional, post-transcriptional and post-translational levels. However, the UGT protein interaction network, which is likely to influence glucuronidation, has received little attention. We investigated the endogenous protein interactome of human UGT1A enzymes in main drug metabolizing non-malignant tissues where UGT expression is most prevalent, using an unbiased proteomics approach. Mass spectrometry analysis of affinity-purified UGT1A enzymes and associated protein complexes in liver, kidney and intestine tissues revealed an intricate interactome linking UGT1A enzymes to multiple metabolic pathways. Several proteins of pharmacological importance such as transferases (including UGT2 enzymes), transporters and dehydrogenases were identified, upholding a potential coordinated cellular response to small lipophilic molecules and drugs. Furthermore, a significant cluster of functionally related enzymes involved in fatty acid β-oxidation, as well as in the glycolysis and glycogenolysis pathways were enriched in UGT1A enzymes complexes. Several partnerships were confirmed by co-immunoprecipitations and co-localization by confocal microscopy. An enhanced accumulation of lipid droplets in a kidney cell model overexpressing the UGT1A9 enzyme supported the presence of a functional interplay. Our work provides unprecedented evidence for a functional interaction between glucuronidation and bioenergetic metabolism.

Transcriptomic analysis identified up-regulation of a solute carrier transporter and UDP glucuronosyltransferases in dogs with aggressive cutaneous mast cell tumours

Gene expression analyses have been recently used in cancer research to identify genes associated with tumorigenesis and potential prognostic markers or therapeutic targets. In the present study, the transcriptome of dogs that had died because of mast cell tumours (MCTs) was characterised to identify a fingerprint having significant influence on prognosis determination and treatment selection. A dataset (GSE50433) obtained using a commercial canine DNA microarray platform was used. The transcriptome of seven biopsies obtained from dogs with histologically confirmed, surgically removed MCTs, treated with chemotherapy, and dead for MCT-related causes, was compared with the transcriptional portrait of 40 samples obtained from dogs with histologically confirmed, surgically removed MCTs and that were still alive at the end of the follow-up period. Among the differentially expressed genes (DEGs), eight transcripts were validated by quantitative real time PCR and their mRNA levels were measured in a cohort of 22 additional MCTs. Statistical analysis identified 375 DEGs (fold change 2, false discovery rate 5%). The functional annotation analysis indicated that the DEGs were associated with drug metabolism and cell cycle pathways. Particularly, members of solute carrier transporter (SLC) and UDP glucuronosyltransferase (UGT) gene families were identified as dysregulated. Principal component analysis (PCA) of the 22 additional MCTs identified the separate cluster dogs dead for MCT-related causes. SLCs and UGTs have been recently recognised in human cancer as important key factors in tumour progression and chemo-resistance. An in-depth analysis of their roles in aggressive canine MCT is warranted in future studies.

Time-Dependent Metabolism of Luteolin by Human UDP-Glucuronosyltransferases and Its Intestinal First-Pass Glucuronidation in Mice

Luteolin is a well-known flavonoid with various pharmacological properties but has low bioavailability due to glucuronidation. This study investigated the time-course of luteolin glucuronidation by 12 human UDP-glucuronosyltransferases (UGTs) and its intestinal first-pass metabolism in mice. Six metabolites, including two novel abundant diglucuronides [3',7-O-diglucuronide (diG) and 4',7-diG] and four known ones, were identified. UGT1A6 and UGT1A9 generated almost only monoglucuronides (G's). The production of 3',7-diG followed a sequential time-dependent process along with decrease of 3'-G mainly by UGT1A1, indicating that 3',7-diG was produced from 3'-G. Metabolism in mice intestine differed from that in humans. Probenecid, a nonspecific UGT inhibitor, did not affect absorption but significantly inhibited production of 7-, 4'-, and 3'-G, and enhanced the formation of another novel metabolite, 5-G, in mice. In conclusion, diglucuronide formation is time-dependent and isoform-specific. UGT1A1 preferentially generates diG, whereas UGT1A6 and UGT1A9 share a preference for G production.

MicroRNA regulation of the major drug-metabolizing enzymes and related transcription factors

Identifying novel mechanisms contributing to patient variability of drug response is a major goal of personalized medicine. Epigenetic regulation of gene expression by microRNA (miRNA) impacts a broad range of cellular processes, but knowledge of its regulation of drug-metabolizing enzymes (DMEs) is more limited. This review provides an introduction to miRNA and their functionality and summarizes known miRNA regulation of DME families, including the cytochrome P450s, UDP-glucuronoslytransferases, glutathione-S-transferases, sulfotransferases and aldo-keto reductases, and the transcription factors known to be involved in DME regulation.

Unravelling the transcriptomic landscape of the major phase II UDP-glucuronosyltransferase drug metabolizing pathway using targeted RNA sequencing

A comprehensive view of the human UDP-glucuronosyltransferase (UGT) transcriptome is a prerequisite to the establishment of an individual's UGT metabolic glucuronidation signature. Here, we uncover the transcriptome landscape of the 10 human UGT gene loci in normal and tumoral metabolic tissues by targeted RNA next-generation sequencing. Alignment on the human hg19 reference genome identifies 234 novel exon-exon junctions. We recover all previously known UGT1 and UGT2 enzyme-coding transcripts and identify over 130 structurally and functionally diverse novel UGT variants. We further expose a revised genomic structure of UGT loci and provide a comprehensive repertoire of transcripts for each UGT gene. Data also uncover a remodelling of the UGT transcriptome occurring in a tissue- and tumor-specific manner. The complex alternative splicing program regulating UGT expression and protein functions is likely critical in determining detoxification capacity of an organ and stress-related responses, with significant impact on drug responses and diseases.

Regulation of UDP-glucuronosyltransferase 1A1 expression and activity by microRNA 491-3p

The UDP-glucuronosyltransferase (UGT) 1A enzymes are involved in the phase II metabolism of many important endogenous and exogenous compounds. The nine UGT1A isoforms exhibit high interindividual differences in expression, but their epigenetic regulation is not well understood. The purpose of the present study was to examine microRNA (miRNA) regulation of hepatic UGT1A enzymes and determine whether or not that regulation impacts enzymatic activity. In silico analysis identified miRNA 491-3p (miR-491-3p) as a potential regulator of the UGT1A gene family via binding to the shared UGT1A 3'-untranslated region common to all UGT1A enzymes. Transfection of miR-491-3p mimic into HuH-7 cells significantly repressed UGT1A1 (P < 0.001), UGT1A3 (P < 0.05), and UGT1A6 (P < 0.05) mRNA levels. For UGT1A1, this repression correlated with significantly reduced metabolism of raloxifene into raloxifene-6-glucuronide (ral-6-gluc; P < 0.01) and raloxifene-4'-glucuronide (ral-4'-gluc; P < 0.01). In HuH-7 cells with repressed miR-491-3p expression, there was a significant increase (~80%; P < 0.01) in UGT1A1 mRNA and a corresponding increase in glucuronidation of raloxifene into ral-6-gluc (50%; P < 0.05) and ral-4'-gluc (22%; P < 0.01). Knockdown of endogenous miR-491-3p in HepG2 cells did not significantly alter UGT1A1 mRNA levels but did increase the formation of ral-6-gluc (50%; P < 0.05) and ral-4'-gluc (34%; P < 0.001). A significant inverse correlation between miR-491-3p expression and both UGT1A3 (P < 0.05) and UGT1A6 (P < 0.01) mRNA levels was observed in a panel of normal human liver specimens, with a significant (P < 0.05) increase in UGT1A3 and UGT1A6 mRNA levels observed in miR-491-3p nonexpressing versus expressing liver specimens. These results suggest that miR-491-3p is an important factor in regulating the expression of UGT1A enzymes in vivo.

Modulation of the UGT2B7 enzyme activity by C-terminally truncated proteins derived from alternative splicing

The enzyme UGT2B7 is one of the most active UDP-glucuronosyltransferases (UGTs) involved in drug metabolism and in maintaining homeostasis of endogenous compounds. We recently reported the existence of 22 UGT2B7 mRNAs, two with a classic 5' region but alternative 3' ends namely UGT2B7_v5 (containing a novel terminal exon 6b) and _v7 (exon 5 excluded) that encode enzymatically inactive isoforms 2 and 4 (i2 and i4), respectively. The v1 mRNA encoding the UGT2B7 enzyme (renamed isoform 1 or i1) is coexpressed with the splice variants v5 and v7 in human liver, kidney, and small intestine and the hepatic cell lines HepG2 and C3A. The presence of alternate v5 and v7 transcripts in isolated polysomes from these hepatic cells further supports endogenous protein translation. Cellular fractionation of clonal HEK293 cell lines overexpressing UGT2B7 isoforms demonstrates that i1, i2, and i4 proteins colocalize in the microsomal/Golgi fraction, whereas i2 and i4 can also be found in the cytosol; a finding sustained by immunofluorescence experiments using tagged proteins. By modifying splice variant abundance in overexpression in HEK293 and HepG2 cells as well as RNA interference experiments in HepG2 and C3A cells, we observe drug glucuronidation phenotypes compatible with variant-mediated repression of UGT2B7 activity without consequent alteration of the apparent enzyme affinity (K(m)). Finally, coimmunoprecipitation experiments support a direct protein-protein interaction of i2 and i4 proteins with the functional UGT2B7 enzyme as a potential causative mechanism. These findings point toward a novel autoregulatory mechanism of the UGT2B7 glucuronidation pathway by naturally occurring alternative i2 and i4 proteins.

Importance of UDP-glucuronosyltransferases 2A2 and 2A3 in tobacco carcinogen metabolism

UDP-glucuronosyltransferase A1 (UGT2A1) is expressed in the lung and exhibits activity against polycyclic aromatic hydrocarbons (PAHs), suggesting UGT2A1 involvement in the local metabolism of PAH tobacco carcinogens. The goal of the present study was to investigate the importance of two additional UGT2A enzymes, UGT2A2 and UGT2A3, in tobacco carcinogen metabolism. Real-time polymerase chain reaction suggested that wild-type UGT2A2 had the highest expression in the breast, followed by trachea > larynx > kidney. A novel splice variant of UGT2A2 lacking exon 3 (termed UGT2A2Δexon3) was investigated, with UGT2A2Δexon3 expression determined to be 25-50% that of wild-type UGT2A2 in all tissues examined. UGT2A3 was determined to be well expressed in the liver and colon, followed by pancreas > kidney > lung > tonsil > trachea > larynx. Cell homogenates prepared from human embryonic kidney (HEK)293 cells overexpressing wild-type UGT2A2 (termed UGT2A2_i1) exhibited glucuronidation activity, as observed by reverse-phase ultra-pressure liquid chromatography, against 1-hydroxy-(OH)-pyrene, 1-naphthol, and hydroxylated benzo(a)pyrene metabolites, whereas homogenates prepared from HEK293 cells overexpressing UGT2A3 only showed activity against simple PAHs like 1-OH-pyrene and 1-naphthol. Activity assays showed the UGT2A2Δexon3 protein (termed UGT2A2_i2) exhibited no detectable glucuronidation activity against all substrates examined; however, coexpression studies suggested that UGT2A2_i2 negatively modulates UGT2A2_i1 activity. Both UGT2A2 and UGT2A3 exhibited no detectable activity against complex PAH proximate carcinogens, tobacco-specific nitrosamines, or heterocyclic amines. These data suggest that, although UGT2A1 is the only UGT2A enzyme active against PAH proximate carcinogens (including PAH diols), both UGTs 2A1 and 2A2 play an important role in the local detoxification of procarcinogenic monohydroxylated PAH metabolites.