Transcriptional regulation of human UDP-glucuronosyltransferase genes

CYP P450 and non-CYP P450 Drug Metabolizing Enzyme Families Exhibit Differential Sensitivities towards Proinflammatory Cytokine Modulation

Compromised hepatic drug metabolism in response to proinflammatory cytokine release is primarily attributed to downregulation of cytochrome P450 (CYP) enzymes. However, whether inflammation also affects other phase I and phase II drug metabolizing enzymes (DMEs), such as the flavin monooxygenases (FMOs), carboxylesterases (CESs), and UDP glucuronosyltransferases (UGTs), remains unclear. This study aimed to decipher the impact of physiologically relevant concentrations of proinflammatory cytokines on expression and activity of phase I and phase II enzymes, to establish a hierarchy of their sensitivity as compared with the CYPs. Hereto, HepaRG cells were exposed to interleukin-6 and interleukin-1β to measure alterations in DME gene expression (24 h) and activity (72 h). Sensitivity of DMEs toward proinflammatory cytokines was evaluated by determining IC50 (potency) and Imax (maximal inhibition) values from the concentration-response curves. Proinflammatory cytokine treatment led to nearly complete downregulation of CYP3A4 (∼98%) but was generally less efficacious at reducing gene expression of the non-CYP DME families. Importantly, FMO, CES, and UGT family members were less sensitive toward interleukin-6 induced inhibition in terms of potency, with IC50 values that were 4.3- to 7.4-fold higher than CYP3A4. Similarly, 18- to 31-fold more interleukin-1β was required to achieve 50% of the maximal downregulation of FMO3, FMO4, CES1, UGT2B4, and UGT2B7 expression. The differential sensitivity persisted at enzyme activity level, highlighting that alterations in DME gene expression during inflammation are predictive for subsequent alterations in enzyme activity. In conclusion, this study has shown that FMOs, CESs, and UGTs enzymes are less impacted by IL-6 and IL-1β treatment as compared with CYP enzymes. SIGNIFICANCE STATEMENT: While the impact of proinflammatory cytokines on CYP expression is well established, their effects on non-CYP phase I and phase II drug metabolism remains underexplored, particularly regarding alterations in drug metabolizing enzyme (DME) activity. This study provides a quantitative understanding of the sensitivity differences to inflammation between DME family members, suggesting that non-CYP DMEs may become more important for the metabolism of drugs during inflammatory conditions due to their lower sensitivity as compared with the CYPs.

Cytochrome P450 and UDP-Glucuronosyltransferase Expressions, Activities, and Induction Abilities in 3D-Cultured Human Renal Proximal Tubule Epithelial Cells

The role of the kidney as an excretory organ for exogenous and endogenous compounds is well recognized, but there is a wealth of data demonstrating that the kidney has significant metabolizing capacity for a variety of exogenous and endogenous compounds that in some cases surpass the liver. The induction of drug-metabolizing enzymes by some chemicals can cause drug-drug interactions and intraindividual variability in drug clearance. In this study, we evaluated the expression and induction of cytochrome P450 (P450) and UDP-glucuronosyltransferase (UGT) isoforms in 3D-cultured primary human renal proximal tubule epithelial cells (RPTEC) to elucidate their utility as models of renal drug metabolism. CYP2B6, CYP2E1, CYP3A4, CYP3A5, and all detected UGTs (UGT1A1, UGT1A4, UGT1A6, UGT1A9, and UGT2B7) mRNA levels in 3D-RPTEC were significantly higher than those in 2D-RPTEC and HK-2 cells and were close to the levels in the human kidney cortex. CYP1B1 and CYP2J2 mRNA levels in 3D-RPTEC were comparable to those in 2D-RPTEC, HK-2 cells, and the human kidney cortex. Midazolam 1'-hydroxylation, trifluoperazine N-glucuronidation, serotonin O-glucuronidation, propofol O-glucuronidation, and morphine 3-glucuronidation in the 3D-RPTEC were significantly higher than the 2D-RPTEC and comparable to those in the HepaRG cells, although bupropion, ebastine, and calcitriol hydroxylations were not different between the 2D- and 3D-RPTEC. Treatment with ligands of the aryl hydrocarbon receptor and farnesoid X receptor induced CYP1A1 and UGT2B4 expression, respectively, in 3D-RPTEC compared with 2D-RPTEC. We provided information on the expression, activity, and induction abilities of P450s and UGTs in 3D-RPTEC as an in vitro human renal metabolism model. SIGNIFICANCE STATEMENT: This study demonstrated that the expression of cytochrome P450s (P450s) and UDP-glucuronosyltransferases (UGTs) in 3D-cultured primary human renal proximal tubule epithelial cells (3D-RPTEC) was higher than those in 2D-cultured primary human renal proximal tubule epithelial cells and HK-2 cells. The results were comparable to that in the human kidney cortex. 3D-RPTEC are useful for evaluating the induction of kidney P450s, UDP-glucuronosyltransferases, and human renal drug metabolism in cellulo.

Glucuronidation dynamics of curcumin and tetrahydrocurcumin for differential structures and chemical reactivities in human liver microsome and uridine diphosphate glucuronosyltransferase 2B7

THC is the main metabolite of curcumin with better bioactivity. This study aimed to explore the factors that cause differences in the bioactivity of curcumin and THC. We analyzed the metabolic activities of curcumin and THC and the factors responsible for the differences in their activities by glucuronidation activity assay, LC-MS, HPLC, homologous sequence comparisons, and molecular docking. Curcumin has higher metabolic activity than THC in HLM and UGT2B7, while the keto-enol isomers of curcumin and THC were distinctly different under different pH, and their structural transformations were hypothesized. Furthermore, UGT1A and UGT2B are differential sequences of curcumin and THC in UGTs. The binding sites and patterns of curcumin and THC in UGT2B7 are markedly different. In summary, the difference in keto-enolic interconversion isomerism between curcumin and THC is the main factor causing the difference in their activities, which provides a scientific basis for the development of curcumin.

Role of blood metabolites in mediating the effect of gut microbiome on the mutated-RAS/BRAF metastatic colorectal cancer-specific survival

BACKGROUND

Recent studies have linked alterations in the gut microbiome and metabolic disruptions to the invasive behavior and metastasis of colorectal cancer (CRC), thus affecting patient prognosis. However, the specific relationship among gut microbiome, metabolite profiles, and mutated-RAS/BRAF metastatic colorectal cancer (M-mCRC) remains unclear. Furthermore, the potential mechanisms and prognostic implications of metabolic changes induced by gut microbiome alterations in patients with M-mCRC still need to be better understood.

METHODS

We conducted Mendelian randomization (MR) to evaluate the causal relationship of genetically predicted 196 gut microbiome features and 1400 plasma metabolites/metabolite ratios on M-mCRC-specific survival. Additionally, we identified significant gut microbiome-metabolites/metabolite ratio associations based on M-mCRC. Metabolite information was annotated, and functional annotation and pathway enrichment analyses were performed on shared proteins corresponding to significant metabolite ratios, aiming to reveal potential mechanisms by which gut microbiome influences M-mCRC prognosis via modulation of human metabolism.

RESULTS

We identified 11 gut microbiome features and 49 known metabolites/metabolite ratios correlated with M-mCRC-specific survival. Furthermore, we identified 17 gut microbiome-metabolite/metabolite ratio associations specific to M-mCRC, involving eight lipid metabolites and three bilirubin degradation products. The shared proteins corresponding to significant metabolite ratios were predominantly localized within the integral component of the membrane and exhibited enzymatic activities such as glucuronosyltransferase and UDP-glucuronosyltransferase, crucial in processes such as glucuronidation, bile secretion, and lipid metabolism. Moreover, these proteins were significantly enriched in pathways related to ascorbate and aldarate metabolism, pentose and glucuronate interconversions, steroid hormone biosynthesis, and bile secretion.

CONCLUSION

Our study offers novel insights into the potential mechanisms underlying the impact of the gut microbiome on the prognosis of M-mCRC. These findings serve as a meaningful reference for exploring potential therapeutic targets and strategies in the future.

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.

Non-canonical transcriptional regulation of the poor prognostic factor UGT2B17 in chronic lymphocytic leukemic and normal B cells

BACKGROUND

High expression of the glycosyltransferase UGT2B17 represents an independent adverse prognostic marker in chronic lymphocytic leukemia (CLL). It also constitutes a predictive marker for therapeutic response and a drug resistance mechanism. The key determinants driving expression of the UGT2B17 gene in normal and leukemic B-cells remain undefined. The UGT2B17 transcriptome is complex and is comprised of at least 10 alternative transcripts, identified by previous RNA-sequencing of liver and intestine. We hypothesized that the transcriptional program regulating UGT2B17 in B-lymphocytes is distinct from the canonical expression previously characterized in the liver.

RESULTS

RNA-sequencing and genomics data revealed a specific genomic landscape at the UGT2B17 locus in normal and leukemic B-cells. RNA-sequencing and quantitative PCR data indicated that the UGT2B17 enzyme is solely encoded by alternative transcripts expressed in CLL patient cells and not by the canonical transcript widely expressed in the liver and intestine. Chromatin accessible regions (ATAC-Seq) in CLL cells mapped with alternative promoters and non-coding exons, which may be derived from endogenous retrotransposon elements. By luciferase reporter assays, we identified key cis-regulatory STAT3, RELA and interferon regulatory factor (IRF) binding sequences driving the expression of UGT2B17 in lymphoblastoid and leukemic B-cells. Electrophoretic mobility shift assays and pharmacological inhibition demonstrated key roles for the CLL prosurvival transcription factors STAT3 and NF-κB in the leukemic expression of UGT2B17.

CONCLUSIONS

UGT2B17 expression in B-CLL is driven by key regulators of CLL progression. Our data suggest that a NF-κB/STAT3/IRF/UGT2B17 axis may represent a novel B-cell pathway promoting disease progression and drug resistance.

Role of microsomal metabolism in bromfenac-induced cytotoxicity

This study delves into the intricate mechanisms underlying drug-induced liver injury (DILI) with a specific focus on bromfenac, the withdrawn nonsteroidal anti-inflammatory drug. DILI is a pervasive concern in drug development, prompting market withdrawals and posing significant challenges to healthcare. Despite the withdrawal of bromfenac due to DILI, the exact role of its microsomal metabolism in inducing hepatotoxicity remains unclear. Herein, employing HepG2 cells with human liver microsomes and UDP-glucuronic acid (UDPGA), our investigation revealed a substantial increase in bromfenac-induced cytotoxicity in the presence of UDPGA, pointing to the significance of UDP-glucuronosyltransferase (UGT)-dependent metabolism in augmenting toxicity. Notably, among the recombinant UGTs examined, UGT2B7 emerged as a pivotal enzyme in the metabolic activation of bromfenac. Metabolite identification studies disclosed the formation of reactive intermediates, with bromfenac indolinone (lactam) identified as a potential mediator of hepatotoxic effects. Moreover, in cytotoxicity experiments, the toxicity of bromfenac lactam exhibited a 34-fold increase, relative to bromfenac. The toxicity of bromfenac lactam was mitigated by nicotinamide adenine dinucleotide phosphate-dependent metabolism. This finding underscores the role of UGT-dependent metabolism in generating reactive metabolites that contribute to the observed hepatotoxicity associated with bromfenac. Understanding these metabolic pathways and the involvement of specific enzymes, such as UGT2B7, provides crucial insights into the mechanisms of bromfenac-induced liver injury. In conclusion, this research sheds light on the metabolic intricacies leading to cytotoxicity induced by bromfenac, especially emphasizing the role of UGT-dependent metabolism and the formation of reactive intermediates like bromfenac lactam. These findings offer insight into the mechanistic basis of DILI and emphasize the importance of understanding metabolism-mediated toxicity.

Identification of an ADME-related gene for forecasting the prognosis and responding to immunotherapy in sarcomas

There are more than 170 subtypes of sarcomas (SARC), which pose a challenge for diagnosis and patient management. Relatively simple or complex karyotypes play an indispensable role in the early diagnosis and effective treatment of SARC. The genes related to absorption, distribution, metabolism, and excretion (ADME) of a drug can serve as prognostic biomarkers of cancer and potential drug targets. In this study, a risk score signature was created. The SARC cohort was downloaded from The Cancer Genome Atlas (TCGA) database, and divided into high-risk group and low-risk group according to the median value of risk score. Compared with high-risk group, low-risk group has a longer survival time, which is also verified in osteosarcoma cohort from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. In addition, the relationship between the signature and immunophenotypes, including status of immune cell infiltration and immune checkpoint expression, was explored. Then, we found that high-risk group is in immunosuppressive status. Finally, we verified that PPARD played a role as a carcinogen in osteosarcoma, which provided a direction for targeted treatment of osteosarcoma in the future. Generally speaking, the signature can not only help clinicians predict the prognosis of patients with SARC, but also provide a theoretical basis for developing more effective targeted drugs in the future.

Androgen Receptor/AP-1 Activates UGT2B15 Transcription to Promote Esophageal Squamous Cell Carcinoma Invasion

Esophageal squamous cell carcinoma (ESCC) is an aggressive epithelial malignancy with poor prognosis. Interestingly, ESCC is strongly characterized by a male-predominant propensity. Our previous study showed that androgen receptor (AR) orchestrated a transcriptional repression program to promote ESCC growth, but it remains unclear whether AR can also activate oncogenic signaling during ESCC progression. In this study, by analyzing our previous AR cistromes and androgen-regulated transcriptomes, we identified uridine diphosphate glucuronosyltransferase family 2 member B15 (UGT2B15) as a bona fide target gene of AR. Mechanistically, AP-1 cofactors played important and collaborative roles in AR-mediated UGT2B15 upregulation. Functional studies have revealed that UGT2B15 promoted invasiveness in vitro and lymph node metastasis in vivo. UGT2B15 was partially responsible for the AR-induced invasive phenotype in ESCC cells. Importantly, simultaneous blocking of AP-1 and AR resulted in stronger inhibition of cell invasiveness compared to inhibiting AP-1 or AR alone. In conclusion, our study reveals the molecular mechanisms underlying the AR-driven ESCC invasion and suggests that the AR/AP1/UGT2B15 transcriptional axis can be potentially targeted in suppressing metastasis in male ESCC patients.

Quantitative analysis of the UDP-glucuronosyltransferase transcriptome in human tissues

UDP-glucuronosyltransferases (UGTs) are phase II drug metabolizing enzymes that play important roles in the detoxification of endogenous and exogenous substrates. The 22 human UGTs belong to four families (UGT1, UGT2, UGT3, and UGT8) and differ in their expression, substrate specificity, UDP-sugar preference, and physiological functions. Differential expression/activity of the UGTs contributes to interperson variability in drug responses and toxicity, hormone homeostasis, and disease/cancer risks. However, in normal tissues, the tissue-specific expression profiles and transcriptional regulation of the UGTs are still not fully understood. In this study, we comprehensively analyzed the transcriptome of 22 UGTs in 54 human tissues/regions using RNAseq data from GTEx. We then validated the findings in the liver and small intestine samples using real-time PCR. Our results showed large interindividual variability across tissues in the expression of each UGT and the overall composition of UGT pools, consisting of different UGTs and their splice isoforms. Our results also revealed coexpression of the UGTs, Cytochrome P450s, and many transcription factors in the liver, suggesting potential coregulation or functional coordination. Our results provide the groundwork for future studies to detail further the regulation of the expression and activity of the UGTs.

Drug-drug interactions that alter the exposure of glucuronidated drugs: Scope, UDP-glucuronosyltransferase (UGT) enzyme selectivity, mechanisms (inhibition and induction), and clinical significance

Drug-drug interactions (DDIs) arising from the perturbation of drug metabolising enzyme activities represent both a clinical problem and a potential economic loss for the pharmaceutical industry. DDIs involving glucuronidated drugs have historically attracted little attention and there is a perception that interactions are of minor clinical relevance. This review critically examines the scope and aetiology of DDIs that result in altered exposure of glucuronidated drugs. Interaction mechanisms, namely inhibition and induction of UDP-glucuronosyltransferase (UGT) enzymes and the potential interplay with drug transporters, are reviewed in detail, as is the clinical significance of known DDIs. Altered victim drug exposure arising from modulation of UGT enzyme activities is relatively common and, notably, the incidence and importance of UGT induction as a DDI mechanism is greater than generally believed. Numerous DDIs are clinically relevant, resulting in either loss of efficacy or an increased risk of adverse effects, necessitating dose individualisation. Several generalisations relating to the likelihood of DDIs can be drawn from the known substrate and inhibitor selectivities of UGT enzymes, highlighting the importance of comprehensive reaction phenotyping studies at an early stage of drug development. Further, rigorous assessment of the DDI liability of new chemical entities that undergo glucuronidation to a significant extent has been recommended recently by regulatory guidance. Although evidence-based approaches exist for the in vitro characterisation of UGT enzyme inhibition and induction, the availability of drugs considered appropriate for use as 'probe' substrates in clinical DDI studies is limited and this should be research priority.

Quantitative Proteomic Analysis of Tibetan Pig Livers at Different Altitudes

In this study, the differences in protein profiles between the livers of Shannan Tibetan pigs (SNT), Linzhi Tibetan pigs (LZT) and Jiuzhaigou Tibetan pigs (JZT) were comparatively analyzed by tandem mass spectrometry-labeling quantitative proteomics. A total of 6804 proteins were identified: 6471 were quantified and 1095 were screened as differentially expressed proteins (DEPs). Bioinformatics analysis results show that, compared with JZT livers, up-regulated DEPs in SNT and LZT livers mainly promoted hepatic detoxification through steroid hormone biosynthesis and participated in lipid metabolism to maintain body energy homeostasis, immune response and immune regulation, while down-regulated DEPs were mainly involved in lipid metabolism and immune regulation. Three proteases closely related to hepatic fatty acid oxidation were down-regulated in enzymatic activity, indicating higher levels of lipid oxidation in SNT and LZT livers than in JZT livers. Down-regulation of the expression of ten immunoglobulins suggests that JZT are more susceptible to autoimmune diseases. It is highly likely that these differences in lipid metabolism and immune-related proteins are in response to the ecological environment at different altitudes, and the findings contribute to the understanding of the potential molecular link between Tibetan pig livers and the environment.

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 Somatic Mutation Landscape of UDP-Glycosyltransferase (UGT) Genes in Human Cancers

The human UDP-glycosyltransferase (UGTs) superfamily has a critical role in the metabolism of anticancer drugs and numerous pro/anti-cancer molecules (e.g., steroids, lipids, fatty acids, bile acids and carcinogens). Recent studies have shown wide and abundant expression of UGT genes in human cancers. However, the extent to which UGT genes acquire somatic mutations within tumors remains to be systematically investigated. In the present study, our comprehensive analysis of the somatic mutation profiles of 10,069 tumors from 33 different TCGA cancer types identified 3427 somatic mutations in UGT genes. Overall, nearly 18% (1802/10,069) of the assessed tumors had mutations in UGT genes with huge variations in mutation frequency across different cancer types, ranging from over 25% in five cancers (COAD, LUAD, LUSC, SKCM and UCSC) to less than 5% in eight cancers (LAML, MESO, PCPG, PAAD, PRAD, TGCT, THYM and UVM). All 22 UGT genes showed somatic mutations in tumors, with UGT2B4, UGT3A1 and UGT3A2 showing the largest number of mutations (289, 307 and 255 mutations, respectively). Nearly 65% (2260/3427) of the mutations were missense, frame-shift and nonsense mutations that have been predicted to code for variant UGT proteins. Furthermore, about 10% (362/3427) of the mutations occurred in non-coding regions (5' UTR, 3' UTR and splice sites) that may be able to alter the efficiency of translation initiation, miRNA regulation or the splicing of UGT transcripts. In conclusion, our data show widespread somatic mutations of UGT genes in human cancers that may affect the capacity of cancer cells to metabolize anticancer drugs and endobiotics that control pro/anti-cancer signaling pathways. This highlights their potential utility as biomarkers for predicting therapeutic efficacy and clinical outcomes.

A transcriptional regulatory network of HNF4α and HNF1α involved in human diseases and drug metabolism

HNF4α and HNF1α are core transcription factors involved in the development and progression of a variety of human diseases and drug metabolism. They play critical roles in maintaining the normal growth and function of multiple organs, mainly the liver, and in the metabolism of endogenous and exogenous substances. The twelve isoforms of HNF4α may exhibit different physiological functions, and HNF4α and HNF1α show varying or even opposing effects in different types of diseases, particularly cancer. Additionally, the regulation of CYP450, phase II drug-metabolizing enzymes, and drug transporters is affected by several factors. This article aims to review the role of HNF4α and HNF1α in human diseases and drug metabolism, including their structures and physiological functions, affected diseases, regulated drug metabolism genes, influencing factors, and related mechanisms. We also propose a transcriptional regulatory network of HNF4α and HNF1α that regulates the expression of target genes related to disease and drug metabolism.

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 Hub Genes in Colorectal Adenocarcinoma by Integrated Bioinformatics

An improved understanding of the molecular mechanism of colorectal adenocarcinoma is necessary to predict the prognosis and develop new target gene therapy strategies. This study aims to identify hub genes associated with colorectal adenocarcinoma and further analyze their prognostic significance. In this study, The Cancer Genome Atlas (TCGA) COAD-READ database and the gene expression profiles of GSE25070 from the Gene Expression Omnibus were collected to explore the differentially expressed genes between colorectal adenocarcinoma and normal tissues. The weighted gene co-expression network analysis (WGCNA) and differential expression analysis identified 82 differentially co-expressed genes in the collected datasets. Enrichment analysis was applied to explore the regulated signaling pathway in colorectal adenocarcinoma. In addition, 10 hub genes were identified in the protein–protein interaction (PPI) network by using the cytoHubba plug-in of Cytoscape, where five genes were further proven to be significantly related to the survival rate. Compared with normal tissues, the expressions of the five genes were both downregulated in the GSE110224 dataset. Subsequently, the expression of the five hub genes was confirmed by the Human Protein Atlas database. Finally, we used Cox regression analysis to identify genes associated with prognosis, and a 3-gene signature (CLCA1–CLCA4–GUCA2A) was constructed to predict the prognosis of patients with colorectal cancer. In conclusion, our study revealed that the five hub genes and CLCA1–CLCA4–GUCA2A signature are highly correlated with the development of colorectal adenocarcinoma and can serve as promising prognosis factors to predict the overall survival rate of patients.

Gill and Liver Transcript Expression Changes Associated With Gill Damage in Atlantic Salmon (Salmo salar)

Gill damage represents a significant challenge in the teleost fish aquaculture industry globally, due to the gill's involvement in several vital functions and direct contact with the surrounding environment. To examine the local and systemic effects accompanying gill damage (which is likely to negatively affect gill function) of Atlantic salmon, we performed a field sampling to collect gill and liver tissue after several environmental insults (e.g., harmful algal blooms). Before sampling, gills were visually inspected and gill damage was scored; gill scores were assigned from pristine [gill score 0 (GS0)] to severely damaged gills (GS3). Using a 44K salmonid microarray platform, we aimed to compare the transcriptomes of pristine and moderately damaged (i.e., GS2) gill tissue. Rank Products analysis (5% percentage of false-positives) identified 254 and 34 upregulated and downregulated probes, respectively, in GS2 compared with GS0. Differentially expressed probes represented genes associated with functions including gill remodeling, wound healing, and stress and immune responses. We performed gill and liver qPCR for all four gill damage scores using microarray-identified and other damage-associated biomarker genes. Transcripts related to wound healing (e.g., neb and klhl41b) were significantly upregulated in GS2 compared with GS0 in the gills. Also, transcripts associated with immune and stress-relevant pathways were dysregulated (e.g., downregulation of snaclec 1-like and upregulation of igkv3) in GS2 compared with GS0 gills. The livers of salmon with moderate gill damage (i.e., GS2) showed significant upregulation of transcripts related to wound healing (i.e., chtop), apoptosis (e.g., bnip3l), blood coagulation (e.g., f2 and serpind1b), transcription regulation (i.e., pparg), and stress-responses (e.g., cyp3a27) compared with livers of GS0 fish. We performed principal component analysis (PCA) using transcript levels for gill and liver separately. The gill PCA showed that PC1 significantly separated GS2 from all other gill scores. The genes contributing most to this separation were pgam2, des, neb, tnnt2, and myom1. The liver PCA showed that PC1 significantly separated GS2 from GS0; levels of hsp70, cyp3a27, pparg, chtop, and serpind1b were the highest contributors to this separation. Also, hepatic acute phase biomarkers (e.g., serpind1b and f2) were positively correlated to each other and to gill damage. Gill damage-responsive biomarker genes and associated qPCR assays arising from this study will be valuable in future research aimed at developing therapeutic diets to improve farmed salmon welfare.

Induction by Phenobarbital of Phase I and II Xenobiotic-Metabolizing Enzymes in Bovine Liver: An Overall Catalytic and Immunochemical Characterization

In cattle, phenobarbital (PB) upregulates target drug-metabolizing enzyme (DME) mRNA levels. However, few data about PB's post-transcriptional effects are actually available. This work provides the first, and an almost complete, characterization of PB-dependent changes in DME catalytic activities in bovine liver using common probe substrates and confirmatory immunoblotting investigations. As expected, PB increased the total cytochrome P450 (CYP) content and the extent of metyrapone binding; moreover, an augmentation of protein amounts and related enzyme activities was observed for known PB targets such as CYP2B, 2C, and 3A, but also CYP2E1. However, contradictory results were obtained for CYP1A, while a decreased catalytic activity was observed for flavin-containing monooxygenases 1 and 3. The barbiturate had no effect on the chosen hydrolytic and conjugative DMEs. For the first time, we also measured the 26S proteasome activity, and the increase observed in PB-treated cattle would suggest this post-translational event might contribute to cattle DME regulation. Overall, this study increased the knowledge of cattle hepatic drug metabolism, and further confirmed the presence of species differences in DME expression and activity between cattle, humans, and rodents. This reinforced the need for an extensive characterization and understanding of comparative molecular mechanisms involved in expression, regulation, and function of DMEs.

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.

Establishment of non-small-cell lung cancer risk prediction model based on prognosis-associated ADME genes

PURPOSE

ADME genes are those involved in the absorption, distribution, metabolism, and excretion (ADME) of drugs. In the present study, a non-small-cell lung cancer (NSCLC) risk prediction model was established using prognosis-associated ADME genes, and the predictive performance of this model was evaluated and verified. In addition, multifaceted difference analysis was performed on groups with high and low risk scores.

METHODS

An NSCLC sample transcriptome and clinical data were obtained from public databases. The prognosis-associated ADME genes were obtained by univariate Cox and lasso regression analyses to build a risk model. Tumor samples were divided into high-risk and low-risk score groups according to the risk score. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of the differentially expressed genes and the differences in the immune infiltration, mutation, and medication reactions in the two groups were studied in detail.

RESULTS

A risk prediction model was established with seven prognosis-associated ADME genes. Its good predictive ability was confirmed by studies of the model's effectiveness. Univariate and multivariate Cox regression analyses showed that the model's risk score was an independent prognostic factor for patients with NSCLC. The study also showed that the risk score closely correlated with immune infiltration, mutations, and medication reactions.

CONCLUSION

The risk prediction model established with seven ADME genes in the present study can predict the prognosis of patients with NSCLC. In addition, significant differences in immune infiltration, mutations, and therapeutic efficacy exist between the high- and low-risk score groups.

Circulating microRNA-30a-5p, microRNA-101-3p, microRNA-140-3p and microRNA-141-3p as potential biomarkers for dexmedetomidine response in pediatric patients

PURPOSE

The aim of this study was to investigate the expression levels of plasma miR-30a-5p, miR-101-3p, miR-140-3p and miR-141-3p and their relationship to dexmedetomidine efficacy and adverse effects in pediatric patients.

METHODS

The expression levels of miR-30a-5p, miR-101-3p, miR-140-3p and miR-141-3p were measured by qRT-PCR in plasma of 133 pediatric patients receiving dexmedetomidine for preoperative sedation. We analyzed the relationship between miRNA abundance and dexmedetomidine response, including sedative effect and adverse effects, and assessed the predictive power of miRNAs for drug response.

RESULTS

Among 133 pediatric patients, 111 patients were dexmedetomidine responders (UMSS ≥ 2) and 22 patients were non-responders (UMSS < 2). We observed higher expression levels of miR-101-3p and miR-140-3p in dexmedetomidine responders compared with non-responders (P < 0.05, P < 0.0001). In contrast, there was no significant difference in the expression levels of miR-30a-5p and miR-141-3p between responders and non-responders (P > 0.05). The plasma levels of miR-101-3p and miR-30a-5p were markedly downregulated in patients who experienced hypotension and bradycardia, respectively (P < 0.05). MiR-101-3p and miR-140-3p demonstrated a potential discriminatory ability between dexmedetomidine responders and non-responders, with AUC of 0.64 (P < 0.05) and 0.77 (P < 0.0001), respectively. The AUC of miR-101-3p in distinguishing patients without hypotension was 0.63 (P < 0.05). The AUC of miR-30a-5p in distinguishing patients without bradycardia was 0.74 (P < 0.05).

CONCLUSION

Our study demonstrated that circulating miR-101-3p, miR-140-3p and miR-30a-5p might be used as a blood-based marker for dexmedetomidine efficacy and safety in pediatric patients.

Identification and validation of ADME genes as prognosis and therapy markers for hepatocellular carcinoma patients

Purpose: ADME genes are genes involved in drug absorption, distribution, metabolism, and excretion (ADME). Previous studies report that expression levels of ADME-related genes correlate with prognosis of hepatocellular carcinoma (HCC) patients. However, the role of ADME gene expression on HCC prognosis has not been fully explored. The present study sought to construct a prediction model using ADME-related genes for prognosis of HCC.

Methods: Transcriptome and clinical data were retrieved from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC), which were used as training and validation cohorts, respectively. A prediction model was constructed using univariate Cox regression and Least Absolute Shrinkage and Selection Operator (LASSO) analysis. Patients were divided into high- and low-risk groups based on the median risk score. The predictive ability of the risk signature was estimated through bioinformatics analyses.

Results: Six ADME-related genes (CYP2C9, ABCB6, ABCC5, ADH4, DHRS13, and SLCO2A1) were used to construct the prediction model with a good predictive ability. Univariate and multivariate Cox regression analyses showed the risk signature was an independent predictor of overall survival (OS). A single-sample gene set enrichment analysis (ssGSEA) strategy showed a significant relationship between risk signature and immune status. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed differentially expressed genes (DEGs) in the high- and low-risk groups were enriched in biological process (BP) associated with metabolic and cell cycle pathways.

Conclusion: A prediction model was constructed using six ADME-related genes for prediction of HCC prognosis. This signature can be used to improve HCC diagnosis, treatment, and prognosis in clinical use.

A validation study of the UGT1A4 rs2011404 variant and the risk of anti-tuberculosis drug-induced hepatotoxicity in an Eastern Chinese Han population

WHAT IS KNOWN AND OBJECTIVE

Anti-tuberculosis (anti-TB) drug-induced hepatotoxicity (ATDH) is a serious adverse drug reaction. A recent study found that the rs2011404 variant of uridine 5'-diphospho-glucuronosyl-transferase 1A4 (UGT1A4) is a marker of susceptibility to ATDH. The present study aimed to validate this relationship in an Eastern Chinese Han anti-TB treatment population.

METHODS

A 1:4 matched case-control study was conducted among anti-TB treatment patients in four regions of Jiangsu. ATDH was diagnosed based on the criteria from the Chinese Society of Hepatology and the updated Roussel Uclaf Causality Assessment Method. A conditional logistic regression model was used to estimate the association between rs2011404 genotypes and the risk of ATDH using odds ratios (ORs) with 95% confidence intervals (95% CIs) and smoking, drinking, hepatoprotectant use and liver diseases as covariates.

RESULTS AND DISCUSSION

A total of 202 ATDH cases and 808 controls were matched according to age, sex and treatment history. After correcting for potential confounding factors, conditional logistic regression analysis indicated no significant differences in genotypes between the two groups (CC vs. TC: OR = 0.933, 95% CI: 0.457-1.905, p = 0.849). Subgroup analysis suggested that patients carrying the CC genotype at rs2011404 in UGT1A4 were at a reduced risk of moderate or severe liver injury (OR = 0.293, 95% CI: 0.093-0.921, p = 0.036).

WHAT IS NEW AND CONCLUSION

Based on a 1:4 individual matched case-control study, possessing the CC genotype at rs2011404 of the UGT1A4 gene reduces the risk of moderate or severe liver injury in Eastern Chinese Han patients receiving anti-TB treatment. Further research is warranted to explain the role of the UGT1A4 gene and its contribution to individual differences in susceptibility to ATDH.

Epigenetics and microRNAs in UGT1As

UDP-glucuronosyltransferases (UGTs) are the main phase II drug-metabolizing enzymes mediating the most extensive glucuronidation-binding reaction in the human body. The UGT1A family is involved in more than half of glucuronidation reactions. However, significant differences exist in the distribution of UGT1As in vivo and the expression of UGT1As among individuals, and these differences are related to the occurrence of disease and differences in metabolism. In addition to genetic polymorphisms, there is now interest in the contribution of epigenetics and noncoding RNAs (especially miRNAs) to this differential change. Epigenetics regulates UGT1As pretranscriptionally through DNA methylation and histone modification, and miRNAs are considered the key mechanism of posttranscriptional regulation of UGT1As. Both epigenetic inheritance and miRNAs are involved in the differences in sex expression and in vivo distribution of UGT1As. Moreover, epigenetic changes early in life have been shown to affect gene expression throughout life. Here, we review and summarize the current regulatory role of epigenetics in the UGT1A family and discuss the relationship among epigenetics and UGT1A-related diseases and treatment, with references for future research.

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.

Regulation of uridine diphosphate-glucuronosyltransferase 1A expression by miRNA-214-5p and miRNA-486-3p

Aim: This study aimed to identify novel miRNAs (miRs) as regulators of UGT1A gene expression and to evaluate them as potential risk factors for the development of liver fibrosis/cirrhosis. Materials & methods: miRNA target sites in UDP-glucuronosyltransferase 1A (UGT1A) 3'-UTR were predicted and confirmed by luciferase assays, quantitative real-time PCR and western blot using HEK293, HepG2 and Huh7 cells. UGT1A and miRNA expression were analyzed in cirrhotic patients and a mouse model of alcoholic liver fibrosis. Results: miR-214-5p and miR-486-3p overexpression reduced UGT1A mRNA, protein levels and enzyme activity in HepG2 and Huh7 cells. miR-486-3p was upregulated in cirrhotic patients and fibrotic mice livers, whereas UGT1A mRNA levels were reduced. Conclusion: In conclusion, we identified two novel miRNAs capable to repress UGT1A expression in vitro and in vivo. Furthermore, miR-486-3p may represent a potential risk factor for the development or progression of liver fibrosis/cirrhosis by means of a reduced UGT1A-mediated detoxification activity.

TBBPA, TBBPS, and TCBPA disrupt hESC hepatic differentiation and promote the proliferation of differentiated cells partly via up-regulation of the FGF10 signaling pathway

Halogenated flame retardants (HFRs), including Tetrabromobisphenol A (TBBPA), Tetrabromobisphenol S (TBBPS), and Tetrachlorobisphenol A (TCBPA), are widely applied in the manufacturing industry to improve fire safety and can be detected in pregnant women's serum at nanomolar levels. Thus, it is necessary to pay attention to the three HFR potential development toxicity, which has not been conclusively addressed yet. The liver is the main organ that detoxifies our body; TBBPA exposure may lead to increased liver weight in rodents. Therefore, in this study, we assessed the developmental hepatic toxicity of the three HFRs with a human embryonic stem cell hepatic differentiation-based system and transcriptomics analyses. We mostly evaluated lineage fate alterations and demonstrated the three HFRs may have common disruptive effects on hepatic differentiation, with TCBPA being significantly more potent. More specifically, the three HFRs up-regulated genes related to cell cycle and FGF10 signaling, at late stages of the hepatic differentiation. This indicates the three chemicals promoted hepatoblast proliferation likely via up-regulating the FGF10 cascade. At the same time, we also presented a powerful way to combine in vitro differentiation and in silico transcriptomic analyses, to efficiently evaluate hazardous materials' adverse effects on lineage fate decisions during early development.

The N6-methyladenosine modification posttranscriptionally regulates hepatic UGT2B7 expression

UDP-glucuronosyltransferases (UGTs) are enzymes catalyzing the glucuronidation of various endogenous and exogenous compounds. In this study, we examined the possibility that N⁶-methyladenosine (m⁶A) modification affects hepatic UGT expression. Treatment of HepaRG cells with 3-deazaadenosine, an inhibitor of RNA methylation, significantly increased UGT1A1, UGT1A3, UGT1A4, UGT1A9, UGT2B7, UGT2B10, and UGT2B15 mRNA levels (1.3- to 2.6-fold). Among them, we focused on UGT2B7 because it most highly contributes to glucuronidation of clinically used drugs. Methylated RNA immunoprecipitation assays revealed that UGT2B7 mRNA in HepaRG cells and human livers is subjected to m⁶A modification mainly at the 5' untranslated region (UTR) and secondarily at the 3'UTR. UGT2B7 mRNA and protein levels in Huh-7 cells were significantly increased by double knockdown of methyltransferase-like 3 (METTL3) and METTL14, whereas those were decreased by knockdown of fat mass and obesity-associated protein (FTO) or alkB homolog 5, RNA demethylase (ALKBH5), suggesting that m⁶A modification downregulates UGT2B7 expression. By experiments using actinomycin D, an inhibitor of transcription, it was demonstrated that ALKBH5-mediated demethylation would attenuate UGT2B7 mRNA degradation, whereas METTL3/METTL14 or FTO-mediated m⁶A modification would alter the transactivity of UGT2B7. Luciferase assays revealed that the promoter region at -118 to -106 has a key role in the decrease in transactivity of UGT2B7 by FTO knockdown. We found that hepatocyte nuclear factor 4α (HNF4α) expression was significantly decreased by knockdown of FTO, indicating that this would be the underlying mechanism of the decreased transactivity of UGT2B7 by knockdown of FTO. Interestingly, treatment with entacapone, which is used for the treatment of Parkinson's disease and is an inhibitor of FTO, decreased HNF4α and UGT2B7 expression. In conclusion, this study clarified that RNA methylation posttranscriptionally controls hepatic UGT2B7 expression.

Aryl hydrocarbon receptor (AHR), integrating energy metabolism and microbial or obesity-mediated inflammation

Aryl hydrocarbon receptor (AHR) has been characterized as multifunctional sensor, integrator and ligand-activated transcription factor of the bHLH/PAS family. Regulation of inflammatory diseases and energy metabolism are among the putative functions of AHR. Challenges in AHR research include marked species differences, and cell, tissue and context dependence of AHR functions. The commentary is focused on AHR's role in the integration between energy expenditure and microbial and non-infectious inflammation, the latter exemplified by obesity-mediated nonalcoholic fatty liver disease. One of the mechanisms controlling energy-consuming inflammation is represented by a signalsome that is involved in retinoic acid-triggered neutrophil differentiation and regulation of the NADPH oxidase complex (NOX). Established signalsome components are AHR, CD38, multiple protein kinases and adaptors. To prevent chronic inflammatory diseases, the complex interplay between a range of inflammatory responses and energy expenditure must be precisely regulated. Surviving an infection requires both pathogen clearance and tissue protection from inflammatory damage. Defenses are energy-consuming anabolic programs. Therefore, anti-inflammatory, catabolic tolerance programs by metabolic reprogramming of macrophages have evolved. Therapeutic options of AHR agonists to reduce chronic inflammatory diseases are discussed.

The Expression Profiles of ADME Genes in Human Cancers and Their Associations with Clinical Outcomes

ADME genes are a group of genes that are involved in drug absorption, distribution, metabolism, and excretion (ADME). The expression profiles of ADME genes within tumours is proposed to impact on cancer patient survival; however, this has not been systematically examined. In this study, our comprehensive analyses of pan-cancer datasets from the Cancer Genome Atlas (TCGA) revealed differential intratumoral expression profiles for ADME genes in 21 different cancer types. Most genes also showed high interindividual variability within cancer-specific patient cohorts. Using Kaplan-Meier plots and logrank tests, we showed that intratumoral expression levels of twenty of the thirty-two core ADME genes were associated with overall survival (OS) in these cancers. Of these genes, five showed significant association with unfavourable OS in three cancers, including SKCM (ABCC2, GSTP1), KIRC (CYP2D6, CYP2E1), PAAD (UGT2B7); sixteen showed significant associations with favourable OS in twelve cancers, including BLCA (UGT2B15), BRCA (CYP2D6), COAD (NAT1), HNSC (ABCB1), KIRC (ABCG2, CYP3A4, SLC22A2, SLC22A6), KIRP (SLC22A2), LIHC (CYP2C19, CYP2C8, CYP2C9, CYP3A5, SLC22A1), LUAD (SLC15A2), LUSC (UGT1A1), PAAD (ABCB1), SARC (ABCB1), and SKCM (ABCB1, DYPD). Overall, these data provide compelling evidence supporting ADME genes as prognostic biomarkers and potential therapeutic targets. We propose that intratumoral expression of ADME genes may impact cancer patient survival by multiple mechanisms that can include metabolizing/transporting anticancer drugs, activating anticancer drugs, and metabolizing/transporting a variety of endogenous molecules involved in metabolically fuelling cancer cells and/or controlling pro-growth signalling pathways.

Expression of UDP Glucuronosyltransferases 2B15 and 2B17 is associated with methylation status in prostate cancer cells

Studies have suggested that abrogated expression of detoxification enzymes, UGT2B15 and UGT2B17, are associated with prostate tumour risk and progression. We investigated the role of EGF on the expression of these enzymes since it interacts with signalling pathways to also affect prostate tumour progression and is additionally associated with decreased DNA methylation. The expression of UGT2B15, UGT2B17, de novo methyltransferases, DNMT3A and DNMT3B was assessed in prostate cancer cells (LNCaP) treated with EGF, an EGFR inhibitor PD16893, and the methyltransferase inhibitor, 5-azacytidine, respectively. The results showed that EGF treatment decreased levels of expression of all four genes and that their expression was reversed by PD16893. Treatment with 5-azacytidine, markedly decreased expression of UGT2B15 and UGT2B17 over 85% as well as significantly decreased expression of DNMT3B, but not the expression of DNMT3A. DNMT3B siRNA treated LNCaP cells had decreased expression of UGT2B15 and UGT2B17, while DNMT3A siRNA treated cells had only moderately decreased UGT2B15 expression. Treatment with DNMT methyltransferase inhibitor, RG108, significantly decreased UGT2B17 expression. Additionally, methylation differences between prostate cancer samples and benign prostate samples from an Illumina 450K Methylation Array study were assessed. The results taken together suggest that hypomethylation of the UGT2B15 and UGT2B17 genes contributes to increased risk of prostate cancer and may provide a putative biomarker or epigenetic target for chemotherapeutics. Mechanistic studies are warranted to determine the role of the methylation marks in prostate cancer.

Influence of breast cancer risk factors and intramammary biotransformation on estrogen homeostasis in the human breast

Understanding intramammary estrogen homeostasis constitutes the basis of understanding the role of lifestyle factors in breast cancer etiology. Thus, the aim of the present study was to identify variables influencing levels of the estrogens present in normal breast glandular and adipose tissues (GLT and ADT, i.e., 17β-estradiol, estrone, estrone-3-sulfate, and 2-methoxy-estrone) by multiple linear regression models. Explanatory variables (exVARs) considered were (a) levels of metabolic precursors as well as levels of transcripts encoding proteins involved in estrogen (biotrans)formation, (b) data on breast cancer risk factors (i.e., body mass index, BMI, intake of estrogen-active drugs, and smoking) collected by questionnaire, and (c) tissue characteristics (i.e., mass percentage of oil, oil%, and lobule type of the GLT). Levels of estrogens in GLT and ADT were influenced by both extramammary production (menopausal status, intake of estrogen-active drugs, and BMI) thus showing that variables known to affect levels of circulating estrogens influence estrogen levels in breast tissues as well for the first time. Moreover, intratissue (biotrans)formation (by aromatase, hydroxysteroid-17beta-dehydrogenase 2, and beta-glucuronidase) influenced intratissue estrogen levels, as well. Distinct differences were observed between the exVARs exhibiting significant influence on (a) levels of specific estrogens and (b) the same dependent variables in GLT and ADT. Since oil% and lobule type of GLT influenced levels of some estrogens, these variables may be included in tissue characterization to prevent sample bias. In conclusion, evidence for the intracrine activity of the human breast supports biotransformation-based strategies for breast cancer prevention. The susceptibility of estrogen homeostasis to systemic and tissue-specific modulation renders both beneficial and adverse effects of further variables associated with lifestyle and the environment possible.

Regional distributions of UDP-glucuronosyltransferase activities toward estradiol and serotonin in the liver and small intestine of cynomolgus macaques

The cynomolgus macaque is a nonhuman primate species that is often used in drug metabolism studies during drug development. However, the localization of UDP-glucuronosyltransferases (UGTs), essential drug-metabolizing enzymes, has not been fully investigated in the liver and small intestine of cynomolgus macaques. In this study, UGT activities were analyzed in liver (five lobes) and small intestine (the duodenum and six sections from the proximal jejunum to the distal ileum) using typical probe substrates of human UGTs: 7-hydroxycoumarin, estradiol, serotonin, propofol, and zidovudine. In liver, UGT activities with respect to all substrates were detected, and the activity levels were similar in all liver lobes of the cynomolgus macaques tested. In contrast, in the small intestine, UGT activities toward all substrates were detected, but their levels generally decreased from jejunum to ileum in cynomolgus macaques. The localization of estradiol 3-O-glucuronosyltransferases and serotonin O-glucuronosyltransferases (which are mainly UGT1A enzymes) appear to be different in liver and small intestine. These results collectively suggest that, in cynomolgus macaques, UGT1As are differentially localized in the small intestine but are relatively homogeneously distributed in the liver.

Pharmacologic Treatment of Transplant Recipients Infected With SARS-CoV-2: Considerations Regarding Therapeutic Drug Monitoring and Drug-Drug Interactions

COVID-19 is a novel infectious disease caused by the severe acute respiratory distress (SARS)-coronavirus-2 (SARS-CoV-2). Several therapeutic options are currently emerging but none with universal consensus or proven efficacy. Solid organ transplant recipients are perceived to be at increased risk of severe COVID-19 because of their immunosuppressed conditions due to chronic use of immunosuppressive drugs (ISDs). It is therefore likely that solid organ transplant recipients will be treated with these experimental antivirals.

A Gilbert syndrome-associated haplotype protects against fatty liver disease in humanized transgenic mice

UDP-glucuronosyltransferases 1 A (UGT1A) enzymes are capable of detoxifying a broad range of endo- and xenobiotic compounds, which contributes to antioxidative effects, modulation of inflammation and cytoprotection. In the presence of low-function genetic UGT1A variants fibrosis development is increased in various diseases. This study aimed to examine the role of common UGT1A polymorphisms in NASH. Therefore, htgUGT1A-WT mice and htgUGT1A-SNP mice (carrying a common human haplotype present in 10% of the white population) were fed a high-fat Paigen diet for 24 weeks. Serum aminotransferase activities, hepatic triglycerides, fibrosis development and UGT1A expression were assessed. Microscopic examination revealed higher hepatic fat deposition and a significant induction of UGT1A gene expression in htgUGT1A-WT mice. In agreement with these observations, lower serum aminotransferase activities and lower expression levels of fibrosis-related genes were measured in htgUGT1A-SNP mice. This was accompanied by reduced PPARα protein levels in htgUGT1A-WT but not in SNP mice. Our data demonstrate a protective effect of a UGT1A SNP haplotype, leading to milder hepatic steatosis and NASH. Higher PPARα protein levels in animals with impaired UGT1A activity are the likely result of reduced glucuronidation of ligands involved in PPARα-mediated fatty acid oxidation and may lead to the observed protection in htgUGT1A-SNP mice.

Network Pharmacology-Based Prediction and Verification of Qingluo Tongbi Formula to Reduce Liver Toxicity of Tripterygium wilfordii via UGT2B7 in Endoplasmic Reticulum

BACKGROUND The hepatotoxicity of Tripterygium wilfordii Hook. f. (TWHF) limits its clinic utilization. Qingluo Tongbi formula (QTF) was formulated based on a basic Chinese medicine theory. Previous studies have confirmed the safety and efficacy of QTF in treating rheumatoid arthritis. Therefore, we considered that TWHF could be detoxified based on its reasonable compatibility with QTF. We investigated the detoxicity mechanism of QTF in reducing the liver toxicity of TWHF. MATERIAL AND METHODS We used network pharmacology to determine the relevant metabolism targets of TWHF, focusing on the phase II metabolic enzymes uridine diphosphate-glucuronosyltransferase 1A1 (UGT1A1), UGT1A6, and UGT2B7. Based on the molecular mechanisms of these predictions and the results of the network analysis, we designed experiments to verify our hypothesis in vivo. We used western blotting, real-time quantitative polymerase chain reaction (RT-qPCR), double immunofluorescence, and laser confocal microscopy to detect the expression of UGTs. Finally, we used transmission electron microscopy to observe the endoplasmic reticulum structure. RESULTS The results confirmed that QTF reversed the TWHF-induced reduction of UGT content in liver microsomes, upregulated UGT1A1 and UGT1A6 but not UGT2B7 in the liver tissue. UGT2B7 expression in the liver and liver microsomes was inconsistent. QTF upregulated the expression of UGT2B7 in the endoplasmic reticulum, and QTF upregulated UGT2B7 expression levels in the endoplasmic reticulum compared with TWHF, which reduced liver toxicity. Structural changes were observed in the endoplasmic reticulum. CONCLUSIONS The Chinese traditional medicine compound QTF can achieve the effect of detoxification by upregulating the expression of UGT2B7 in the endoplasmic reticulum.

Differential Role of Liver X Receptor (LXR) α and LXRβ in the Regulation of UDP-Glucuronosyltransferase 1A1 in Humanized UGT1 Mice

Liver X receptors (LXRs), LXRα and LXRβ, are nuclear receptors that regulate the metabolism of cholesterol and bile acids and are activated by oxysterols. Humanized UGT1 (hUGT1) mice express the 9-human UGT1A genes associated with the UGT1 locus in a Ugt1-null background. The expression of UGT1A1 is developmentally delayed in the liver and intestines, resulting in the accumulation of serum bilirubin during the neonatal period. Induction of UGT1A1 in newborn hUGT1 mice leads to rapid reduction in total serum bilirubin (TSB) levels, a phenotype measurement that allows for an accurate prediction on UGT1A1 expression. When neonatal hUGT1 mice were treated by oral gavage with the LXR agonist T0901317, TSB levels were dramatically reduced. To determine the LXR contribution to the induction of UGT1A1 and the lowering of TSB levels, experiments were conducted in neonatal hUGT1/Lxrα ^-/- , hUGT1/Lxrβ ^-/- , and hUGT1/Lxrαβ ^-/- mice treated with T0901317. Induction of liver UGT1A1 was dependent upon LXRα, with the induction pattern paralleling induction of LXRα-specific stearoyl CoA desaturase 1. However, the actions of T0901317 were also shown to display a lack of specificity for LXR, with the induction of liver UGT1A1 in hUGT1/Lxrαβ ^-/- mice, a result associated with activation of both pregnane X receptor and constitutive androstane receptor. However, the LXR agonist GW3965 was highly selective toward LXRα, showing no impact on lowering TSB values or inducing UGT1A1 in hUGT1/Lxrα ^-/- mice. An LXR-specific enhancer site on the UGT1A1 gene was identified, along with convincing evidence that LXRα is crucial in maintaining constitutive expression of UGT1A1 in adult hUGT1 mice. SIGNIFICANCE STATEMENT: It has been established that activation of LXRα, and not LXRβ, is responsible for the induction of liver UGT1A1 and metabolism of serum bilirubin in neonatal hUGT1 mice. Although induction of the human UGT1A1 gene is initiated at a newly characterized LXR enhancer site, allelic deletion of the Lxrα gene drastically reduces the constitutive expression of liver UGT1A1 in adult hUGT1 mice. Combined, these findings indicate that LXRα is critical for the developmental expression of UGT1A1.

Molecular characterization of functional UDP-glucuronosyltransferases 1A and 2B in common marmosets

UDP-glucuronosyltransferases (UGTs) are essential drug-conjugation enzymes that metabolize a variety of endobiotic and xenobiotic substrates. The molecular characteristics of UGTs have been extensively investigated in humans, but remain to be investigated in common marmosets, a nonhuman primate species widely used in drug metabolism studies. In this study, 11 UGT cDNAs (UGT1A1, 1A3, 1A4, 1A6, 1A7, and 1A9; and UGT2B49, 2B50, 2B51, 2B52, and 2B53) were isolated and characterized in marmosets. Marmoset UGT1As had high sequence identities (89-93%) with human UGT1As, but the sequence identities of marmoset UGT2Bs were lower (82-86%). Marmoset UGTs were found to be phylogenetically close to human UGTs. Just as human UGT1As do, marmoset UGT1A genes shared exons 2-5 and contained a variable exon 1 unique to each gene; in contrast, marmoset UGT2B genes contained six unique exons. Moreover, marmoset and human UGT1A and UGT2B gene clusters were located in corresponding regions in their respective genomes. Among the five tissue types tested, marmoset UGT mRNAs were most abundantly expressed in liver, jejunum, and/or kidney, i.e., in tissues important for drug metabolism, just as human UGTs are. Among the 11 marmoset UGT mRNAs investigated, marmoset UGT1A9, 1A4, and 1A6 mRNAs were the most abundantly expressed in liver, small intestine, and kidney, respectively. Marmoset liver microsomes and recombinant UGT1A proteins catalyzed the glucuronidation of the same substrates that human UGT1As catalyze, including estradiol, trifluoperazine, 4-methylumbelliferone, serotonin, 4-nitrophenol, and propofol. Trifluoperazine was glucuronidated by marmoset liver microsomes, but not by any of the UGT1A isoforms examined under the present conditions. These results collectively suggest that functional marmoset UGTs have generally similar molecular characteristics to human UGTs.

Polymorphic Expression of UGT1A9 is Associated with Variable Acetaminophen Glucuronidation in Neonates: A Population Pharmacokinetic and Pharmacogenetic Study

INTRODUCTION

Acetaminophen (paracetamol, APAP) is widely used as an analgesic and antipyretic drug in children and neonates. A number of enzymes contribute to the metabolism of acetaminophen, and genetic factors might be important to explain variability in acetaminophen metabolism among individuals.

METHODS

The current investigation utilized a previously published parent-metabolite population pharmacokinetic model describing acetaminophen glucuronidation, sulfation, and oxidation to examine the potential role of genetic variability on the relevant metabolic pathways. Neonates were administered 30-min intravenous infusions of acetaminophen 15 mg/kg every 12 h (< 28 weeks' gestational age [GA]) or every 8 h (≥ 28 weeks GA) for 48 h. A total of 18 sequence variations (SVs) in UDP-glucuronosyltransferase (UGT), sulfotransferase (SULT), and cytochrome P450 (CYP) genes from 33 neonates (aged 1-26 days) were examined in a stepwise manner for an effect on the metabolic formation clearance of acetaminophen by glucuronidation (UGT), sulfation (SULT), and oxidation (CYP). The stepwise covariate modeling procedure was performed using NONMEM^® version 7.3.

RESULTS

Incorporation of genotype as a covariate for one SV located in the UGT1A9 gene promoter region (rs3832043, - 118 > insT, T₉ > T₁₀) significantly improved model fit (likelihood ratio test, p < 0.001) and reduced between-subject variability in glucuronide formation clearance. Individuals with the UGT1A9 T₁₀ polymorphism, indicating insertion of an additional thymidine nucleotide, had a 42% reduction in clearance to APAP-glucuronide as compared to their wild-type counterparts.

CONCLUSION

This study shows a pharmacogenetic effect of an SV in the UGT1A9 promoter region on the metabolism of acetaminophen in neonates.

Current trends in drug metabolism and pharmacokinetics

Pharmacokinetics (PK) is the study of the absorption, distribution, metabolism, and excretion (ADME) processes of a drug. Understanding PK properties is essential for drug development and precision medication. In this review we provided an overview of recent research on PK with focus on the following aspects: (1) an update on drug-metabolizing enzymes and transporters in the determination of PK, as well as advances in xenobiotic receptors and noncoding RNAs (ncRNAs) in the modulation of PK, providing new understanding of the transcriptional and posttranscriptional regulatory mechanisms that result in inter-individual variations in pharmacotherapy; (2) current status and trends in assessing drug-drug interactions, especially interactions between drugs and herbs, between drugs and therapeutic biologics, and microbiota-mediated interactions; (3) advances in understanding the effects of diseases on PK, particularly changes in metabolizing enzymes and transporters with disease progression; (4) trends in mathematical modeling including physiologically-based PK modeling and novel animal models such as CRISPR/Cas9-based animal models for DMPK studies; (5) emerging non-classical xenobiotic metabolic pathways and the involvement of novel metabolic enzymes, especially non-P450s. Existing challenges and perspectives on future directions are discussed, and may stimulate the development of new research models, technologies, and strategies towards the development of better drugs and improved clinical practice.

In vitro glucuronidation of 7-hydroxycoumarin derivatives in intestine and liver microsomes of Beagle dogs

Beagle dog is a standard animal model for evaluating nonclinical pharmacokinetics of new drug candidates. Glucuronidation in intestine and liver is an important first-pass drug metabolic pathway, especially for phenolic compounds. This study evaluated the glucuronidation characteristics of several 7-hydroxycoumarin derivatives in beagle dog's intestine and liver in vitro. To this end, glucuronidation rates of 7-hydroxycoumarin (compound 1), 7-hydroxy-4-trifluoromethylcoumarin (2), 6-methoxy-7-hydroxycoumarin (3), 7-hydroxy-3-(4-tolyl)coumarin (4), 3-(4-fluorophenyl)coumarin (5), 7-hydroxy-3-(4-hydroxyphenyl)coumarin (6), 7-hydroxy-3-(4-methoxyphenyl)coumarin (7), and 7-hydroxy-3-(1H-1,2,4-tirazole)coumarin (8) were determined in dog's intestine and liver microsomes, as well as recombinant dog UGT1A enzymes. The glucuronidation rates of 1, 2 and 3 were 3-10 times higher in liver than in small intestine microsomes, whereas glucuronidation rates of 5, 6, 7 and 8 were similar in microsomes from both tissues. In the colon, glucuronidation of 1 and 2 was 3-5 times faster than in small intestine. dUGT1A11 glucuronidated efficiently all the substrates and was more efficient catalyst for 8 than any other dUGT1A. Other active enzymes were dUGT1A2 that glucuronidated efficiently 2, 3, 4, 5, 6 and 7, while dUGT1A10 glucuronidated efficiently 1, 2, 3, 4, 5 and 7. Kinetic analyses revealed that the compounds' K_m values varied between 1.1 (dUGT1A10 and 2) and 250 µM (dUGT1A7 and 4). The results further strengthen the concept that dog intestine has high capacity for glucuronidation, and that different dUGT1As mediate glucuronidation with distinct substrates selectivity in dog and human.

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.

Qualitative and quantitative differences in estrogen biotransformation in human breast glandular and adipose tissues: implications for studies using mammary biospecimens

Because of its assumed role in breast cancer etiology, estrogen biotransformation (and interaction of compounds therewith) has been investigated in human biospecimens for decades. However, little attention has been paid to the well-known fact that large inter-individual variations exist in the proportion of breast glandular (GLT) and adipose (ADT) tissues and less to adequate tissue characterization. To assess the relevance of this, the present study compares estrogen biotransformation in GLT and ADT. GLT and ADT were isolated from 47 reduction mammoplasty specimens derived from women without breast cancer and were characterized histologically and by their percentages of oil. Levels of 12 unconjugated and five conjugated estrogens were analyzed by GC- and UHPLC-MS/MS, respectively, and levels of 27 transcripts encoding proteins involved in estrogen biotransformation by Taqman^® probe-based PCR. Unexpectedly, one-third of specimens provided neat GLT only after cryosection. Whereas 17β-estradiol, estrone, and estrone-3-sulfate were detected in both tissues, estrone-3-glucuronide and 2-methoxy-estrone were detected predominately in GLT and ADT, respectively. Estrogen levels as well as ratios 17β-estradiol/estrone and estrone-3-sulfate/estrone differed significantly between GLT and ADT, yet less than between individuals. Furthermore, estrogen levels in GLT and ADT correlated significantly with each other. In contrast, levels of most transcripts encoding enzymes involved in biotransformation differed more than between individuals and did not correlate between ADT and GLT. Thus, mixed breast tissues (and plasma) will not provide meaningful information on local estrogen biotransformation (and interaction of compounds therewith) whereas relative changes in 17β-estradiol levels may be investigated in the more abundant ADT.

The UDP-Glycosyltransferase (UGT) Superfamily: New Members, New Functions, and Novel Paradigms

UDP-glycosyltransferases (UGTs) catalyze the covalent addition of sugars to a broad range of lipophilic molecules. This biotransformation plays a critical role in elimination of a broad range of exogenous chemicals and by-products of endogenous metabolism, and also controls the levels and distribution of many endogenous signaling molecules. In mammals, the superfamily comprises four families: UGT1, UGT2, UGT3, and UGT8. UGT1 and UGT2 enzymes have important roles in pharmacology and toxicology including contributing to interindividual differences in drug disposition as well as to cancer risk. These UGTs are highly expressed in organs of detoxification (e.g., liver, kidney, intestine) and can be induced by pathways that sense demand for detoxification and for modulation of endobiotic signaling molecules. The functions of the UGT3 and UGT8 family enzymes have only been characterized relatively recently; these enzymes show different UDP-sugar preferences to that of UGT1 and UGT2 enzymes, and to date, their contributions to drug metabolism appear to be relatively minor. This review summarizes and provides critical analysis of the current state of research into all four families of UGT enzymes. Key areas discussed include the roles of UGTs in drug metabolism, cancer risk, and regulation of signaling, as well as the transcriptional and posttranscriptional control of UGT expression and function. The latter part of this review provides an in-depth analysis of the known and predicted functions of UGT3 and UGT8 enzymes, focused on their likely roles in modulation of levels of endogenous signaling pathways.

Deregulation of the Genes that Are Involved in Drug Absorption, Distribution, Metabolism, and Excretion in Hepatocellular Carcinoma

Genes involved in drug absorption, distribution, metabolism, and excretion (ADME) are called ADME genes. Currently, 298 genes that encode phase I and II drug metabolizing enzymes, transporters, and modifiers are designated as ADME genes by the PharmaADME Consortium. ADME genes are highly expressed in the liver and their levels can be influenced by liver diseases such as hepatocellular carcinoma (HCC). In this study, we obtained RNA-sequencing and microRNA (miRNA)-sequencing data from 371 HCC patients via The Cancer Genome Atlas liver hepatocellular carcinoma project and performed ADME gene-targeted differential gene expression analysis and expression correlation analysis. Two hundred thirty-three of the 298 ADME genes (78%) were expressed in HCC. Of these genes, almost one-quarter (58 genes) were significantly downregulated, while only 6% (15) were upregulated in HCC relative to healthy liver. Moreover, one-half (14/28) of the core ADME genes (CYP1A2, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2E1, CYP3A4, NAT1, NAT2, UGT2B7, SLC22A1, SLCO1B1, and SLCO1B3) were downregulated. In addition, about one-half of the core ADME genes were positively correlated with each other and were also positively (AHR, ARNT, HNF4A, PXR, CAR, PPARA, and RXRA) or negatively (PPARD and PPARG) correlated with transcription factors known as ADME modifiers. Finally, we show that most miRNAs known to regulate core ADME genes are upregulated in HCC. Collectively, these data reveal 1) an extensive transcription factor-mediated ADME coexpression network in the liver that efficiently coordinates the metabolism and elimination of endogenous and exogenous compounds; and 2) a widespread deregulation of this network in HCC, most likely due to deregulation of both transcriptional and post-transcriptional (miRNA) pathways.