Cross-Talk between Alternatively Spliced UGT1A Isoforms and Colon Cancer Cell Metabolism

A Comprehensive Bioinformatic Analysis of RNA-seq Datasets Reveals a Differential and Variable Expression of Wildtype and Variant UGT1A Transcripts in Human Tissues and Their Deregulation in Cancers

The UGT1A locus generates over 60 different alternatively spliced transcripts and 30 circular RNAs. To date, v2 and v3 transcripts are the only variant UGT1A transcripts that have been functionally characterized. Both v2 and v3 transcripts encode the same inactive variant UGT1A proteins (i2s) that can negatively regulate glucuronidation activity and influence cancer cell metabolism. However, the abundance and interindividual variability in the expression of v2 and v3 transcripts in human tissues and their potential deregulation in cancers have not been comprehensively assessed. To address this knowledge gap, we quantified the expression levels of v1, v2, and v3 transcripts using RNA-seq datasets with large cohorts of normal tissues and paired normal and tumor tissues from patients with six different cancer types (liver, kidney, colon, stomach, esophagus, and bladder cancer). We found that v2 and v3 abundance varied significantly between different tissue types, and that interindividual variation was also high within the same tissue type. Moreover, the ratio of v2 to v3 variants varied between tissues, implying their differential regulation. Our results showed higher v2 abundance in gastrointestinal tissues than liver and kidney tissues, suggesting a more significant negative regulation of glucuronidation by i2 proteins in gastrointestinal tissues than in liver and kidney tissues. We further showed differential deregulation of wildtype (v1) and variant transcripts (v2, v3) in cancers that generally increased the v2/v1 and/or v3/v1 expression ratios in tumors compared to normal tissues, indicating a more significant role of the variants in tumors. Finally, we report ten novel UGT1A transcripts with novel 3' terminal exons, most of which encode variant proteins with a similar structure to UGT1A_i2 proteins. These findings further emphasize the diversity of the UGT1A transcriptome and proteome.

Targeted splicing therapy: new strategies for colorectal cancer

RNA splicing is the process of forming mature mRNA, which is an essential phase necessary for gene expression and controls many aspects of cell proliferation, survival, and differentiation. Abnormal gene-splicing events are closely related to the development of tumors, and the generation of oncogenic isoform in splicing can promote tumor progression. As a main process of tumor-specific splicing variants, alternative splicing (AS) can promote tumor progression by increasing the production of oncogenic splicing isoforms and/or reducing the production of normal splicing isoforms. This is the focus of current research on the regulation of aberrant tumor splicing. So far, AS has been found to be associated with various aspects of tumor biology, including cell proliferation and invasion, resistance to apoptosis, and sensitivity to different chemotherapeutic drugs. This article will review the abnormal splicing events in colorectal cancer (CRC), especially the tumor-associated splicing variants arising from AS, aiming to offer an insight into CRC-targeted splicing therapy.

A Non-Canonical Role for the Glycosyltransferase Enzyme UGT2B17 as a Novel Constituent of the B Cell Receptor Signalosome

In chronic lymphocytic leukemia (CLL), an elevated glycosyltransferase UGT2B17 expression (UGT2B17HI) identifies a subgroup of patients with shorter survival and poor drug response. We uncovered a mechanism, possibly independent of its enzymatic function, characterized by an enhanced expression and signaling of the proximal effectors of the pro-survival B cell receptor (BCR) pathway and elevated Bruton tyrosine kinase (BTK) phosphorylation in B-CLL cells from UGT2B17HI patients. A prominent feature of B-CLL cells is the strong correlation of UGT2B17 expression with the adverse marker ZAP70 encoding a tyrosine kinase that promotes B-CLL cell survival. Their combined high expression levels in the treatment of naïve patients further defined a prognostic group with the highest risk of poor survival. In leukemic cells, UGT2B17 knockout and repression of ZAP70 reduced proliferation, suggesting that the function of UGT2B17 might involve ZAP70. Mechanistically, UGT2B17 interacted with several kinases of the BCR pathway, including ZAP70, SYK, and BTK, revealing a potential therapeutic vulnerability. The dual SYK and JAK/STAT6 inhibitor cerdulatinib most effectively compromised the proliferative advantage conferred by UGT2B17 compared to the selective BTK inhibitor ibrutinib. Findings point to an oncogenic role for UGT2B17 as a novel constituent of BCR signalosome also connected with microenvironmental signaling.

The Uridine diphosphate (UDP)-glycosyltransferases (UGTs) superfamily: the role in tumor cell metabolism

UDP-glycosyltransferases (UGTs), important enzymes in biotransformation, control the levels and distribution of numerous endogenous signaling molecules and the metabolism of a wide range of endogenous and exogenous chemicals. The UGT superfamily in mammals consists of the UGT1, UGT2, UGT3, and UGT8 families. UGTs are rate-limiting enzymes in the glucuronate pathway, and in tumors, they are either overexpressed or underexpressed. Alterations in their metabolism can affect gluconeogenesis and lipid metabolism pathways, leading to alterations in tumor cell metabolism, which affect cancer development and prognosis. Glucuronidation is the most common mammalian conjugation pathway. Most of its reactions are mainly catalyzed by UGT1A, UGT2A and UGT2B. The body excretes UGT-bound small lipophilic molecules through the bile, urine, or feces. UGTs conjugate a variety of tiny lipophilic molecules to sugars, such as galactose, xylose, acetylglucosamine, glucuronic acid, and glucose, thereby inactivating and making water-soluble substrates, such as carcinogens, medicines, steroids, lipids, fatty acids, and bile acids. This review summarizes the roles of members of the four UGT enzyme families in tumor function, metabolism, and multiple regulatory mechanisms, and its Inhibitors and inducers. The function of UGTs in lipid metabolism, drug metabolism, and hormone metabolism in tumor cells are among the most important topics covered.

The Glycosyltransferase Pathway: An Integrated Analysis of the Cell Metabolome

Nucleotide sugar-dependent glycosyltransferases (UGTs) are critical to the homeostasis of endogenous metabolites and the detoxification of xenobiotics. Their impact on the cell metabolome remains unknown. Cellular metabolic changes resulting from human UGT expression were profiled by untargeted metabolomics. The abundant UGT1A1 and UGT2B7 were studied as UGT prototypes along with their alternative (alt.) splicing-derived isoforms displaying structural differences. Nineteen biochemical routes were modified, beyond known UGT substrates. Significant variations in glycolysis and pyrimidine pathways, and precursors of the co-substrate UDP-glucuronic acid were observed. Bioactive lipids such as arachidonic acid and endocannabinoids were highly enriched by up to 13.3-fold (p < 0.01) in cells expressing the canonical enzymes. Alt. UGT2B7 induced drastic and unique metabolic perturbations, including higher glucose (18-fold) levels and tricarboxylic acid cycle (TCA) cycle metabolites and abrogated the effects of the UGT2B7 canonical enzyme when co-expressed. UGT1A1 proteins promoted the accumulation of branched-chain amino acids (BCAA) and TCA metabolites upstream of the mitochondrial oxoglutarate dehydrogenase complex (OGDC). Alt. UGT1A1 exacerbated these changes, likely through its interaction with the OGDC component oxoglutarate dehydrogenase-like (OGDHL). This study expands the breadth of biochemical pathways associated with UGT expression and establishes extensive connectivity between UGT enzymes, alt. proteins and other metabolic processes.

Identification of the prognostic signature based on genomic instability-related alternative splicing in colorectal cancer and its regulatory network

Background: Colorectal cancer (CRC) is a heterogeneous disease with many somatic mutations defining its genomic instability. Alternative Splicing (AS) events, are essential for maintaining genomic instability. However, the role of genomic instability-related AS events in CRC has not been investigated.

Methods: From The Cancer Genome Atlas (TCGA) program, we obtained the splicing profiles, the single nucleotide polymorphism, transcriptomics, and clinical information of CRC. Combining somatic mutation and AS events data, a genomic instability-related AS signature was constructed for CRC. Mutations analyses, clinical stratification analyses, and multivariate Cox regression analyses evaluated this signature in training set. Subsequently, we validated the sensitivity and specificity of this prognostic signature using a test set and the entire TCGA dataset. We constructed a nomogram for the prognosis prediction of CRC patients. Differentially infiltrating immune cells were screened by using CIBERSORT. Inmmunophenoscore (IPS) analysis was used to evaluate the response of immunotherapy. The AS events-related splicing factors (SF) were analyzed by Pearson’s correlation. The effects of SF regulating the prognostic AS events in proliferation and migration were validated in Caco2 cells.

Results: A prognostic signature consisting of seven AS events (PDHA1-88633-ES, KIAA1522-1632-AP, TATDN1-85088-ES, PRMT1-51042-ES, VEZT-23786-ES, AIG1-77972-AT, and PHF11-25891-AP) was constructed. Patients in the high-risk score group showed a higher somatic mutation. The genomic instability risk score was an independent variable associated with overall survival (OS), with a hazard ratio of a risk score of 1.537. The area under the curve of receiver operator characteristic curve of the genomic instability risk score in predicting the OS of CRC patients was 0.733. Furthermore, a nomogram was established and could be used clinically to stratify patients to predict prognosis. Patients defined as high-risk by this signature showed a lower proportion of eosinophils than the low-risk group. Patients with low risk were more sensitive to anti-CTLA4 immunotherapy. Additionally, HSPA1A and FAM50B were two SF regulating the OS-related AS. Downregulation of HSPA1A and FAM50B inhibited the proliferation and migration of Caco2 cells.

Conclusion: We constructed an ideal prognostic signature reflecting the genomic instability and OS of CRC patients. HSPA1A and FAM50B were verified as two important SF regulating the OS-related AS.

Consensus molecular subtype differences linking colon adenocarcinoma and obesity revealed by a cohort transcriptomic analysis

Colorectal cancer (CRC) is the third-leading cause of cancer-related deaths in the United States and worldwide. Obesity—a worldwide public health concern—is a known risk factor for cancer including CRC. However, the mechanisms underlying the link between CRC and obesity have yet to be fully elucidated in part because of the molecular heterogeneity of CRC. We hypothesized that obesity modulates CRC in a consensus molecular subtype (CMS)-dependent manner. RNA-seq data and associated tumor and patient characteristics including body weight and height data for 232 patients were obtained from The Cancer Genomic Atlas–Colon Adenocarcinoma (TCGA-COAD) database. Tumor samples were classified into the four CMSs with the CMScaller R package; body mass index (BMI) was calculated and categorized as normal, overweight, and obese. We observed a significant difference in CMS categorization between BMI categories. Differentially expressed genes (DEGs) between obese and overweight samples and normal samples differed across the CMSs, and associated prognostic analyses indicated that the DEGs had differing associations on survival. Using Gene Set Enrichment Analysis, we found differences in Hallmark gene set enrichment between obese and overweight samples and normal samples across the CMSs. We constructed Protein-Protein Interaction networks and observed differences in obesity-regulated hub genes for each CMS. Finally, we analyzed and found differences in predicted drug sensitivity between obese and overweight samples and normal samples across the CMSs. Our findings support that obesity impacts the CRC tumor transcriptome in a CMS-specific manner. The possible associations reported here are preliminary and will require validation using in vitro and animal models to examine the CMS-dependence of the genes and pathways. Once validated the obesity-linked genes and pathways may represent new therapeutic targets to treat colon cancer in a CMS-dependent manner.

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.

A FACS-Free Purification Method to Study Estrogen Signaling, Organoid Formation, and Metabolic Reprogramming in Mammary Epithelial Cells

Few in vitro models are used to study mammary epithelial cells (MECs), and most of these do not express the estrogen receptor α (ERα). Primary MECs can be used to overcome this issue, but methods to purify these cells generally require flow cytometry and fluorescence-activated cell sorting (FACS), which require specialized instruments and expertise. Herein, we present in detail a FACS-free protocol for purification and primary culture of mouse MECs. These MECs remain differentiated for up to six days with >85% luminal epithelial cells in two-dimensional culture. When seeded in Matrigel, they form organoids that recapitulate the mammary gland's morphology in vivo by developing lumens, contractile cells, and lobular structures. MECs express a functional ERα signaling pathway in both two- and three-dimensional cell culture, as shown at the mRNA and protein levels and by the phenotypic characterization. Extracellular metabolic flux analysis showed that estrogens induce a metabolic switch favoring aerobic glycolysis over mitochondrial respiration in MECs grown in two-dimensions, a phenomenon known as the Warburg effect. We also performed mass spectrometry (MS)-based metabolomics in organoids. Estrogens altered the levels of metabolites from various pathways, including aerobic glycolysis, citric acid cycle, urea cycle, and amino acid metabolism, demonstrating that ERα reprograms cell metabolism in mammary organoids. Overall, we have optimized mouse MEC isolation and purification for two- and three-dimensional cultures. This model represents a valuable tool to study how estrogens modulate mammary gland biology, and particularly how these hormones reprogram metabolism during lactation and breast carcinogenesis.

Emerging roles for UDP-glucuronosyltransferases in drug resistance and cancer progression

The best-known role of UDP-glucuronosyltransferase enzymes (UGTs) in cancer is the metabolic inactivation of drug therapies. By conjugating glucuronic acid to lipophilic drugs, UGTs impair the biological activity and enhance the water solubility of these agents, driving their elimination. Multiple clinical observations support an expanding role for UGTs as modulators of the drug response and in mediating drug resistance in numerous cancer types. However, accumulating evidence also suggests an influence of the UGT pathway on cancer progression. Dysregulation of the expression and activity of UGTs has been associated with the progression of several cancers, arguing for UGTs as possible mediators of oncogenic pathways and/or disease accelerators in a drug-naive context. The consequences of altered UGT activity on tumour biology are incompletely understood. They might be associated with perturbed levels of bioactive endogenous metabolites such as steroids and bioactive lipids that are inactivated by UGTs or through non-enzymatic mechanisms, thereby eliciting oncogenic signalling cascades. This review highlights the evidence supporting dual roles for the UGT pathway, affecting cancer progression and drug resistance. Pharmacogenomic testing of UGT profiles in patients and the development of therapeutic options that impair UGT actions could provide useful prognostic and predictive biomarkers and enhance the efficacy of anti-cancer drugs.

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.

Advances in the Understanding of Protein-Protein Interactions in Drug Metabolizing Enzymes through the Use of Biophysical Techniques

In recent years, a growing appreciation has developed for the importance of protein-protein interactions to modulate the function of drug metabolizing enzymes. Accompanied with this appreciation, new methods and technologies have been designed for analyzing protein-protein interactions both in vitro and in vivo. These technologies have been applied to several classes of drug metabolizing enzymes, including: cytochrome P450's (CYPs), monoamine oxidases (MAOs), UDP-glucuronosyltransferases (UGTs), glutathione S-transferases (GSTs), and sulfotransferases (SULTs). In this review, we offer a brief description and assessment of the impact of many of these technologies to the study of protein-protein interactions in drug disposition. The still expanding list of these techniques and assays has the potential to revolutionize our understanding of how these enzymes carry out their important functions in vivo.

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.