Advertisement
Research ArticleMetabolism, Transport, and Pharmacogenomics

Regulation of UDP-Glucuronosyltransferase 1A1 Expression and Activity by MicroRNA 491-3p

Douglas F. Dluzen, Dongxiao Sun, Anna C. Salzberg, Nate Jones, Ryan T. Bushey, Gavin P. Robertson and Philip Lazarus
Journal of Pharmacology and Experimental Therapeutics March 2014, 348 (3) 465-477; DOI: https://doi.org/10.1124/jpet.113.210658

Abstract

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

Introduction

MicroRNA (miRNA) are dynamic epigenetic regulators of gene expression in both normal tissue development and tumor biology. MicroRNA post-transcriptionally repress protein expression primarily by binding to the 3′-untranslated region (UTR) of target mRNA and by the inhibition of protein synthesis (Carthew and Sontheimer, 2009; Guo et al., 2010). Acting as a genetic switch and/or fine-tuner of protein expression, miRNA regulation can account for not only mild (e.g., 2-fold) protein repression, but can also cause mRNA destabilization to dramatically reduce protein output (Baek et al., 2008; Selbach et al., 2008; Mukherji et al., 2011).

MicroRNA are short, ∼22-nucleotide long, single-stranded RNA species that are gene-encoded and transcribed as a primary transcript (pri-miRNA) by RNA polymerase II. Pri-miRNA form a hairpin loop structure that is recognized and processed in the nucleus and then exported to the cytoplasm as a precursor miRNA (pre-miRNA) consisting of a double-stranded duplex and stem-loop. The precursor is further processed to cleave off the stem-loop structure, and the RNA duplex is then recognized by the RNA-induced silencing complex (RISC) consisting of Argonaute-2 and associated proteins, which degrades the passenger strand of the RNA duplex and leaves the mature miRNA (Carthew and Sontheimer, 2009). The mature miRNA guides RISC to target mRNA 3′-UTRs (less commonly, the mRNA open-reading frame or 5′-UTR), with target recognition primarily driven by the hybridization of the miRNA “seed sequence” to the mRNA target. This seed sequence of a target gene is defined as nucleotides 2–8 from the 5′ end of the mature miRNA strand and forms perfect complementation to its mRNA target (Bartel, 2009).

The role of miRNA in tissue development is continually being investigated, but research has shown miRNA regulation to be a major contributor to tissue and cell homeostasis (Khraiwesh et al., 2010; Liu et al., 2012), signaling (Valadi et al., 2007), and oncogenesis (Hatziapostolou et al., 2011). The impact of miRNAs on drug response has not been studied as extensively. MicroRNA regulate expression of several human cytochrome P450 phase I metabolic enzymes, including the major drug and hormone metabolizing CYPs 3A4, 2E1, 1B1, and 24 (Tsuchiya et al., 2006; Komagata et al., 2009; Pan et al., 2009; Mohri et al., 2010), but there are few studies on other phase I or II enzymes.

The UDP-glucuronosyltransferase (UGT) phase II metabolic enzyme family primarily consists of two large subfamilies, the UGT1As and 2Bs. The UGT1A family, located on chromosome 2q37, codes for nine functional protein isoforms, sharing alternative first exons spliced to common exons 2–5. The common fifth exon contains a shared 3′-UTR (Fig. 1A). The UGT2B family members, located on chromosome 4q13, are expressed as individual genes encoded by six exons, each with a unique 3′-UTR (Mackenzie et al., 2005; Nagar and Remmel, 2006). The UGTs are generally highly expressed in the liver, the primary organ of metabolism, but are also expressed extrahepatically in tissues such as colon, prostate, lung, pancreas, kidney, and tissues of the head and neck and other aerodigestive tract tissues (Nakamura et al., 2008; Bushey and Lazarus, 2012; Jones et al., 2012). Each UGT1A gene is regulated by an individual promoter, driving tissue-specific expression. UGTs 1A1, 1A3, 1A4, 1A5, 1A6, and 1A9 are the primary hepatic isoforms; UGTs 1A7, 1A8, and 1A10 are exclusively expressed extrahepatically (Nagar and Remmel, 2006; Nakamura et al., 2008; Izukawa et al., 2009). The UGTs’ prime numerous endogenous compounds include bilirubin (Bosma et al., 1994) and steroid hormones (Belanger et al., 2003), as well as exogenous compounds, including drugs, chemotherapeutic agents, and carcinogens (Court et al., 2001; Kemp et al., 2002; Nagar and Blanchard, 2006; Nagar and Remmel, 2006; Chen et al., 2007, 2008, 2010, 2012; Balliet et al., 2009; Lazarus et al., 2009; Mizuma, 2009; Rowland et al., 2013; Sun et al., 2013) for excretion from the body by catalyzing the conjugation of a sugar moiety (e.g., glucuronide) to its substrate. The exon 1-encoded amino-termini of individual UGT enzymes determines their unique substrate-specificities; however, there is a high degree of substrate overlap between UGT family members (Nagar and Remmel, 2006; Mackenzie et al., 2010).

Fig. 1.
Fig. 1.
Fig. 1.

The UGT1A gene locus and potential miR-491-3p MRE within the common UGT1A 3′-UTR. (A) The nine alternative UGT1A first exons code for unique enzymes. UGTs 1A1, 1A6, and 1A9 (dark boxes) are further highlighted in B. (B) Alignment of three representative UGT1A mRNAs (1A1, 1A6, and 1A9) to highlight the common 3′-UTR region. The miR-491-3p MRE begins +92 nucleotides from the stop codon for all UGT enzymes, and the region spans 24 nucleotides in length. The miR-148a (spanning 22 nucleotides) and miR-136 (spanning 24 nucleotides) MREs begin at +512 and +610 nucleotides, respectively, from the stop codon. The predicted hybridization structure between the UGT1A 3′-UTR and the miR-491-3p, miR-148a, and miR-136 MREs are enlarged at the bottom. Highlighted above the black bars are the MRE seed sequences. (C) Schematic of the pGL3-Promoter luciferase vector with the SV40 promoter and the SV40 poly-A tail, including the XbaI 3′-UTR restriction digest cloning site. The pGL3-UGT1A 3′-UTR luciferase vector contains the first 683 nucleotides of the UGT1A family 3′-UTR, encompassing the wild-type miR-491-3p MRE. The pGL3-UGT1A 3′-UTR/Seed Deletion vector contains all sequences of the pGL3-UGT1A 3′-UTR vector except for a four nucleotide deletion within the miR-491-3p MRE seed sequence. X indicates deleted base positions. (D) Luciferase activity of UGT1A 3′-UTR luciferase reporter vectors cotransfected with miR-491-3p mimic or scrambled miRNA control in HEK293 cells. Columns represent mean ± S.E. of three independent experiments. *P < 0.05.

The UGT1A family members exhibit extensive interindividual variability of expression that contributes to variability in patient response and toxicity (Court et al., 2001). Although polymorphisms in the UGT genes and altered transcriptional regulation can contribute to altered expression and/or activity of UGT1A enzymes (Wiener et al., 2004; Hoskins et al., 2007; Chen et al., 2010; Gallagher et al., 2010), these mechanisms do not explain the observations that UGT1A mRNA and protein levels often do not correlate (Izukawa et al., 2009; Oda et al., 2012; Ohtsuki et al., 2012) and that there is a high degree of variability in UGT mRNA and protein expression within different tissues and between different individuals (Izukawa et al., 2009; Ohno and Nakajin, 2009). Such observations suggest an epigenetic mechanism of posttranscriptional regulation may be affecting UGT protein expression. We hypothesized that miRNA regulate expression of UGT1A isoforms and that this may contribute to interindividual variability in UGT1A expression.

Our present study is the first to examine miRNA interactions within the UGT1A family of enzymes. In silico approaches were used to screen for possible miRNA interactions with the UGT1A family of enzymes, and functional studies were then performed to examine the potential role of miRNA on UGT1A gene expression. Evidence is presented demonstrating that microRNA 491-3p (miR-491-3p) may be an important regulator of hepatic UGT1A expression and activity.

Materials and Methods

Chemicals and Reagents.

The pGL3-Promoter and pRL-TK Renilla plasmids were obtained from Promega (Madison, WI). All synthesized DNA oligos used for 3′-UTR amplification, site-directed mutagenesis, and polymerase chain reaction (PCR) analysis were from Integrated DNA Technologies (Coralville, IA). Lipofectamine 2000 transfection reagent was from Life Technologies (Carlsbad, CA). The miRVana miRNA mimic miR-491-3p (#4464066), negative control #1 mimic (#4464058), miRVana miRNA inhibitor miR-491-3p (#4464084), and negative control #1 inhibitor (#4464076) were purchased from Ambion (Austin, TX). Uridine diphosphate glucuronic acid, raloxifene, alamethicin, and bovine serum albumin were obtained from Sigma-Aldrich (St. Louis, MO). Raloxifene-6-glucuronide (ral-6-gluc), raloxifene-4′-glucuronde (ral-4′-gluc), raloxifene-d4, ral-6-gluc-d4, ral-4′gluc-d4, and epirubicin hydrochloride were purchased from Toronto Research Chemical (Toronto, ON, Canada). Rabbit anti-UGT2B7 and anti-calnexin antibodies were from BD Biosciences (San Jose, CA) and Cell Signaling Technology (Danvers, MA), respectively. Goat anti-rabbit secondary antibody conjugated to hydrogen peroxidase was from Thermo Scientific (Waltham, MA). All other chemicals used were purchased from Fisher Scientific (Pittsburgh, PA) unless otherwise specified.

Cell Lines and Culture Conditions.

Human embryonic kidney (HEK) cell line 293, human hepatocellular carcinoma cell lines HepG2 and Hep3B, human liver adenocarcinoma cell line SK-HEP-1, human prostate carcinoma cell line LNCaP, human colon adenocarcinoma cell line Caco-2, human breast adenocarcinoma cell line MCF-7, and human lung carcinoma cell line A-549 were purchased from the American Type Culture Collection (ATCC, Manassas, VA). The human hepatocellular carcinoma cell line HuH-7 was a kind gift from Dr. Jianming Hu (Penn State Hershey College of Medicine, Hershey, PA). Hep3B, HEK293, A-549, and HuH-7 cells were cultured in Dulbecco’s modified Eagle’s medium (GIBCO, Carlsbad, CA) supplemented with 10% fetal bovine serum (FBS) (Atlanta Biologicals, Lawrenceville, GA) and 1% penicillin/streptomycin (GIBCO). HepG2 and Caco-2 cells were cultured in Dulbecco’s modified Eagle’s medium supplemented with 10% FBS, 1% penicillin/streptomycin, and 1% nonessential amino acids (Lonza, Basel, Switzerland). SK-HEP-1 and LNCaP cells were cultured in RPMI 1640 (GIBCO) supplemented with 10% FBS, 1% penicillin/streptomycin, and 1% nonessential amino acid. MCF-7 cells were cultured in RPMI 1640 supplemented with 10% FBS and 1% penicillin/streptomycin. All cells were grown and maintained at 37°C with 5% CO2.

Tissues and miRNA Isolation.

Normal colon and endometrium specimens (n = 5 each) were obtained from the tissue bank at Penn State University College of Medicine, and normal liver specimens and their matching total RNA were obtained from the H. Lee Moffitt Cancer Center Tissue Procurement Facility (n = 39). All protocols involving the collection and analysis of tissue specimens from these tissue banks were approved by the institutional review boards at their respective institutions and were in accordance with assurances filed with and approved by the U.S. Department of Health and Human Services. Normal jejunum tissue specimens were purchased from the Sun Health Research Institute (Sun City, AZ). All tissue samples were isolated and frozen at −70°C within 2 hours after surgery. Colon, liver, jejunum, endometrium, and cell line total RNA were extracted using the RNeasy Mini Kit (Qiagen, Hilden, Germany). All total RNA samples from cell lines were subject to on-column DNAse digestion during RNA purification (Qiagen). Pooled normal breast total RNA was purchased from the Biochain Institute (Hayward, CA) while normal lung, pancreas, larynx, trachea, and kidney RNA was purchased from Clontech (Mountain View, CA) or Agilent Technologies (Santa Clara, CA). Small RNA (<200 nt) containing the miRNA fraction was isolated and purified from total RNA using the mirVana miRNA Isolation Kit (Ambion). All RNA concentrations were ascertained using a Nanodrop ND-1000 spectrophotometer (Thermo Scientific) and was eluted and stored in RNase-, DNase-free water in a −80°C freezer.

MicroRNA Binding Site Predictions.

The 3′-UTR of the UGT1A enzyme locus was obtained from the University of California–Santa Cruz Genome Browser (hg18 assembly). MicroRNA binding site predictions were obtained using 1) TargetScan (Lewis et al., 2005), scored with the Total Context+ score as described in Garcia et al. (2011), and 2) miRanda, algorithm v3.0 (Betel et al., 2008) with the following parameters: Gap Open Penalty, −8.00; Gap Extend, −2.00; Score Threshold, 50.00; Energy Threshold, −20.00 kcal/mol; Scaling Parameter, 4.00.

Quantitative Real-Time PCR.

We synthesized cDNA from mRNA using total RNA and the SuperScript First-Strand Synthesis Kit (Invitrogen, Carlsbad, CA). MicroRNA cDNAs were synthesized using the Taqman MicroRNA Reverse Transcription Kit (Ambion). Taqman gene expression assays (Applied Biosystems, Carlsbad, CA) were used to amplify UGTs 1A1 (Hs02511055_s1), 1A3 (Hs04194492_g1), 1A4 (Hs016555285_s1), 1A5 (Hs01374521_s1), 1A6 (Hs01592477_m1), 1A9 (Hs02516855_gH), RPLPO (Hs99999902_m1), and 2B7 (Hs02556232_s1) in HuH-7 and HepG2 cells. Taqman miRNA assays (Ambion) were used to amplify the expression of miR-491-3p (cat. no. 4427974, ID 002360) and RNU6B (cat. no. 4427975, ID 001093) in all cell lines and tissue samples. PCR reactions were performed in 10-µl reactions in 384-well plates using an ABI 7900HT Sequence Detection System with incubations performed at 50°C for 2 minutes; 95°C for 10 minutes; and 40 cycles of 95°C for 15 seconds, 60°C for 1 minute. Reactions included 2× Universal PCR Master Mix (Applied Biosystems), Taqman gene expression primers or Taqman miRNA primers, and corresponding cDNA according to the manufacturer’s protocol. Each plate was run with a negative control (no DNA template), and all assays were performed in quadruplicate. Gene expression was compared with an endogenous internal control [ribosomal protein, large, P0 (RFLPO) for mRNA or U6 small nuclear RNA 2 (RNU6B) for miRNA] using the 2−ΔΔCt method (Livak and Schmittgen, 2001). Cycle threshold (Ct) values were determined using the SDS 2.4 software (Applied Biosystems), and amplification Ct values higher than 36 cycles were designated as below the limit of detection (BLD). Samples lacking any amplification curves were also designated BLD.

Construction of Reporter Plasmids.

The luciferase report plasmids used in this study were constructed by inserting the common UGT1A 3′-UTR into the XbaI restriction site located downstream of the luciferase reporter gene in the pGL3-Promoter vector. Briefly, primers (sense, 5′-GCTATCTAGAGAAGTGGGTGGGAAATAAGGTAA-3′; and antisense, 5′-GCTATCTAGAGAAACTTGCCCAGCACTTCATAGCT-3′) modified with the XbaI digestion site at both ends were used to amplify the UGT1A subfamily 3′-UTR (corresponding to nucleotides +1618 to +2301 of the human UGT1A1 mRNA, and 683 nucleotides in length) using genomic DNA isolated from the HEK293 cell line. The PCR-amplified region was cloned into the XbaI restriction site of the pGL3-promoter vector using standard protocols. The pGL3-seed deletion plasmid was created by performing site-directed mutagenesis using the QuikChange II Site Directed Mutagenesis Kit (Agilent Technologies) and the sense and antisense primers ′5-TCATTTTATTCTTATTAAGGAAATACTTTAAATTAATCAGCCCCAGAGTGCTT-3′ and 5′-AAGCACTCTGGGGCTGATTAATTTAAAGTATTTCCTTAATAAGAATAAAATGA-3′, respectively. Nucleotide sequences of all plasmids used in this study were confirmed by DNA sequencing analysis performed at the Penn State University Nucleic Acid Facility (State College, PA).

Luciferase Assays.

The pGL3-Promoter vector cloned with the UGT1A 3′-UTR (termed pGL3-UGT1A 3′-UTR) was cotransfected with the pRL-TK Renilla control vector into HEK293 cells. The day before transfection, HEK293 cells were seeded onto 24-well plates at a density of approximately 50,000 cells/well; after 24 hours, 380 ng of pGL3-UGT1A 3′-UTR plasmid and 20 ng pRL-TK plasmid were cotransfected together with various concentrations of either scrambled miRNA control or miRNA mimic using Lipofectamine 2000 transfection reagent. HEK293 cells were harvested 24 hours after transfection using passive lysis buffer, and luciferase activity was measured with a luminometer (Bio-tek Synergy HT; BioTek Instruments, Winooski, VT) using the Dual-Luciferase Reporter Assay System (Promega). Renilla luciferase served as the internal control.

Glucuronidation Assays.

Preparation of cell homogenates for glucuronidation activity assays was performed as previously described (Dellinger et al., 2007) after cells were transfected with 50 μl of Lipofectamine 2000 in 10 cm dishes for 48 hours according to manufacturer’s protocol. HuH-7 or HepG2 cells were transfected with either 50 nM scrambled miRNA controls, 50 nM miR-491-3p mimic, or 50 nM miR-491-3p inhibitor. Collected cell pellets were freeze-thawed three times in liquid nitrogen and subjugated to 3 × 10-second pulses using a handheld Bio-Vortexer (Biospec, Bartlesville, OK) before storing of 50-µl aliquots at −80°C. Protein concentrations within the cell homogenates were quantified using the BCA Protein Assay Kit (Pierce, Rockford, IL) and measured using an Appliskan Luminometor and SkanIT Software v2.3 (Thermo Scientific).

The glucuronidation assays using homogenates from HuH-7 and HepG2 cells were performed essentially as described elsewhere (Sun et al., 2007, 2013). HuH-7 and HepG2 cell homogenate (100–300 µg) was incubated with 50 µM raloxifene for 90 minutes, or 500 µM epirubicin hydrochloride for 60 minutes. Raloxifene glucuronidation reactions were terminated by the addition of 25 µl of cold 100% acetonitrile and centrifuged before the drying of the supernatant in a speed vac and reconstitution in 40 µl of 1:1 water/acetonitrile; epirubicin reactions were terminated in the presence of 25 µl of cold 100% acetonitrile, centrifuged for 20 minutes at 13,000g, and the supernatant was collected for analysis on ultra-pressure liquid chromatography–tandem mass spetrometry.

Stock solutions of raloxifene, ral-6-gluc, and ral-4′-gluc and their deuterated internal standards were prepared in dimethylsulfoxide. Raloxifene, ral-6-gluc, and ral-4′-gluc were combined into a standard stock solution with final concentrations of 25, 100, and 100 µg/ml, respectively. Deuterated internal standards were combined in dimethylsulfoxide to final working concentrations of 50 µg/ml for raloxifene, ral-6-gluc, and ral-4′-gluc. The combined standard solution was then serially diluted in 1:1 water/acetonitrile to make standard working solutions from 50 ng/ml to 51.2 µg/ml for raloxifene, and 1.25 ng/ml to 1.28 µg/ml for ral-6-gluc and ral-4′-gluc. Standard calibration samples were prepared daily by spiking 1.0 µl of the combined internal standard into each of the serial-diluted working solutions (500 ng/ml final concentrations for each deuterated internal standard). All solutions were kept at −20°C before use.

Calibration standards, as well as formed raloxifene and epirubicin glucuronides, were analyzed using a Waters Acquity ultra-pressure liquid chromatography-UV detector system (Waters Corporation, Milford, MA) with a 1.7 µ Acquity UPLIC BEH C18 analytical column (2.1 mm × 50 mm; Waters Corporation, Dublin, Ireland) in series with a 0.2 µm Waters assay frit filter (2.1 mm). Raloxifene gradient elution conditions were performed as previously described elsewhere (Sun et al., 2013). Epirubicin gradient elution conditions were performed using a flow rate of 0.5 ml/min, starting with 95% buffer A (0.1% formic acid in water) and 5% buffer B (0.1% formic acid in acetonitrile) for 30 seconds, a subsequent linear gradient to 50% buffer B over 2.5 minutes, and then 100% buffer B maintained over the next 2 minutes. Raloxifene and epirubicin glucuronides were confirmed by their sensitivity to the treatment of β-glucuronidase.

Quantification of raloxifene, ral-6-gluc, and ral-4′-gluc was performed as previously described elsewhere (Sun et al., 2013). Quantification of epirubicin and epirubicin-glucuronide was determined based on the ratio of epirubicin-glucuronide versus unconjugated epirubicin after calculating the area under the curve for the epirubicin-glucuronide and epirubicin peaks using a known amount of epirubicin in each reaction as a reference. Data were quantified using the MassLynxTM NT 4.1 software with QuanLynx program (Waters Corporation). All standard curves were constructed by plotting the ratio of analyte peak area to corresponding internal standard peak area. All experiments were performed in triplicate in independent assays.

Western Blot Analysis.

UGT2B7 protein levels were determined via immunoblotting. HuH-7 and HepG2 protein homogenate was adjusted to equal volumes of loading buffer and heated at 90°C for 10 minutes. Samples were run at 90 V on a 10% acrylamide gel and transferred to a polyvinylidene difluoride membrane for 2 hours at 33 V. The polyvinylidene difluoride membranes were blocked in 5% milk in Tris-buffered saline with Tween-20 for 1 hour, probed with UGT2B7 or calnexin primary antibody (1:1000 dilution for each) overnight at 4°C, and washed three times in Tris-buffered saline with Tween-20, followed by goat anti-rabbit secondary antibody (1:5000 dilution). Protein bands were visualized using the Amersham ECL Prime Western Blotting Detection Reagent (GE Healthcare, Pittsburgh, PA) and Hyoblot CL autoradiography film (Denville Scientific, Metuchen, NJ). Calnexin expression levels served as a loading control.

Statistical Analysis.

Statistical analysis was performed using Graphpad Prism 5 (GraphPad Software, La Jolla, CA). For studies involving miR-491-3p miRNA mimic or miR-491-3p inhibitors, the Student’s t test (two-tailed) was used to compare experimental groups with scrambled miRNA controls. For statistical analysis of expression levels of genes labeled BLD, a Ct value of 40 was assigned to generate the necessary 2−ΔΔCt values needed for statistical comparison. For the analysis of UGT1A1, 1A3, 1A6, and 1A9 versus miR-491-3p miRNA expression in human liver specimens, a one-tailed Spearman correlation was used. In this analysis, for those specimens with no miR-491-3p expression, they were assigned a Ct value half that of the lowest miR-491-3p expressing liver specimen. For comparison of UGT1A1, 1A3, 1A6, and 1A9 mRNA levels in miR-491-3p–expressing versus nonexpressing human liver specimens, the one-tailed Mann–Whitney t test was used. P < 0.05 was considered statistically significant.

Results

In Silico Analysis of miRNA Binding to UGT1A mRNA.

The human UGT1A gene locus consists of nine alternative first exons spliced to common exons 2–5, with all nine UGT1A enzymes sharing the same 3′-UTR (Fig. 1A). The mRNA corresponding to the UGT1A 3′-UTR was analyzed by the miRanda (v3.0) and TargetScan algorithms to identify miRNAs that could potentially target UGT1A mRNAs (Enright et al., 2003; John et al., 2004; Lewis et al., 2005). miRanda is an excellent model algorithm due to its emphasis on stringent base-pairing of the seed region supplemented by downstream binding near the 3′ end of the miRNA gene (Betel et al., 2008). TargetScan is a similarly effective model algorithm due to its emphasis on contributions from binding position, 3′ pairing, and miRNA target site abundance (Garcia et al., 2011). MicroR-491-3p was the highest scoring prediction by miRanda to bind the UGT1A common 3′-UTR, followed by miR-148a and miR-136. TargetScan also identified these three miRNAs as strong candidates to bind the UGT1A 3′-UTR, with miR-491-3p having the highest overall Total Context+ score of the three miRNAs. MicroR-491-3p targets the common UGT1A 3′-UTR at a position 92–114 base pairs (bp) downstream of the UGT1A stop codon [termed the miR-491-3p miRNA recognition element (MRE)] (Fig. 1B). This predicted MRE contains the canonical 8-nucleotide miRNA seed sequence of perfect complementation between miR-491-3p and the UGT1A 3′-UTR (Fig. 1B, bottom, underlined). The MREs for miRs 148a and miR-136 target regions corresponding to bp 512–533 and 610–633, respectively, downstream of the UGT1A stop codon. Like that observed for miR-491-3p, both contain a 8-nucleotide seed sequence of perfect complementation with the UGT1A 3′-UTR (Fig. 1B, bottom, underlined).

Effect of miR-491-3p on UGT1A 3′-UTR–Mediated Luciferase Activity.

To determine whether any of the predicted miRNAs described (miR-491-3p, miR-148a, and miR-136) could potentially affect UGT1A expression in vitro, the UGT1A 3′-UTR was cloned into the luciferase pGL3-promoter vector immediately 3′ of the luciferase open-reading frame (Fig. 1C). The pGL3-UGT1A 3′-UTR plasmid was transiently cotransfected into HEK293 cells with scrambled negative control or miR-491-3p, miR-148a, or miR-136 miRNA mimics together with the Renilla luciferase plasmid (serving as a transfection control). HEK293 cells are an ideal cell line for studying heterologous expression of UGT1A enzymes because they do not exhibit endogenous UGT1A expression (Nakamura et al., 2008). Compared with the scrambled miRNA-transfected control, luciferase activity was significantly repressed at both 1 nM (P < 0.05) and 2 nM (P < 0.05) miR-491-3p mimic concentrations, corresponding to a 43% and 48% loss, respectively, of luciferase activity (Fig. 1D). To confirm that this inhibition was due to miR-491-3p binding to the UGT1A 3′-UTR miR-491-3p MRE, a 4-bp deletion mutation was created within the MRE seed sequence within the pGL3-UGT1A 3′-UTR (termed pGL3-UGT1A3′-UTR/SeedDel; Fig. 1C). Although no significant difference in luciferase activity was observed between the seed deletion mutant and the negative scrambled miRNA control, a significant (P < 0.05) 1.5-fold increase in luciferase expression was observed between cells overexpressing the wild-type pGL3-UGT1A 3′-UTR versus cells overexpressing the mutant pGL3-UGT1A 3′-UTR/SeedDel plasmids using 2 nM miR-491-3p mimic (Fig. 1D). No significant alteration in luciferase activity was observed in cotransfections with either the miR-148a and miR-136 mimics (data not shown).

MicroR-491-3p Expression in UGT-Expressing Tissues and Cell Lines.

To determine the levels of miR-491-3p expression in cells with known UGT expression, total RNA was screened from a variety of normal human tissues and cell lines by quantitative real-time reverse transcription (qRT)-PCR (Fig. 2). MicroR-491-3p expression was highest in colon (25.6 ± 5.9) followed by trachea (12.0 ± 4.0) > lung (6.5 ± 0.81) > breast (3.0 ± 1.1) > liver (1.0 ± 0.11; set as reference; Fig. 2A). Of the 39 livers analyzed, 12 exhibited no expression of miR-491-3p; of the remaining samples, up to 4.0-fold differences in miR-491-3p expression were observed for 26 of the samples (a 7.8-fold difference was observed for the remaining sample; Fig. 2B). High miR-491-3p expression was also observed in Hep3B, HuH-7, HepG2, and SK-HEP-1 cells, exhibiting between 6.2- to 15.0-fold higher levels of miR-491-3p expression as compared with the lowest expressing individual liver specimen (Fig. 2B). Low levels of miR-491-3p expression were observed in the Caco-2, MCF-7, HEK293, and LNCaP cell lines; no expression was detected in jejunum, pancreas, kidney, larynx, and endometrium, or in the A-549 lung carcinoma cell line (data not shown).

Fig. 2.
Fig. 2.

Tissue and cell line expression of miR-491-3p and UGT1A family members. (A) Relative expression levels of miR-491-3p in several UGT-expressing tissues. Expression levels were quantified using qRT-PCR and set relative to the lowest-expressing tissue (i.e., liver). (B) Expression levels of mature miR-491-3p in individual normal livers and in liver cancer cell lines. Expression levels were quantified by qRT-PCR. MicroR-491-3p expression was adjusted to the RNU6B endogenous control gene in all tissues and cell lines. Expression is shown relative to the lowest overall miR-491-3p–expressing sample (i.e., liver #13; set to 1.0). (C) Comparison of UGT1A mRNA levels between HepG2 and HuH-7 cell lines. mRNA expression was adjusted to the RPLPO endogenous control gene and shown relative to the highest expressing UGT enzyme in both cell lines (UGT1A1 in HuH-7 cells; set to 1.0). The y axis is broken into two segments, with a gap between 0.08 and 0.2 to better adjust for the high expression of UGT1A1 in HuH-7 cells. Columns represent the mean ± S.E. of three independent replicates. ***P < 0.001.

UGT1A Expression in miR-491-3p–Expressing Cell Lines.

To identify a hepatic cell culture model that could be used to monitor the effects of miR-491-3p on hepatic UGT1A expression, UGT1A expression was quantified by qRT-PCR in both the HuH-7 and HepG2 cell lines. UGT1A1 exhibited the highest level of expression of any hepatically expressed UGT in both cell lines, with the levels in HuH-7 cells 40-fold higher than that observed in HepG2 cells (P < 0.001; Fig. 2C). UGT1A1 exhibited > 400- and 15-fold higher levels of expression than the next-most well-expressed UGT enzyme, UGT1A6, in HuH-7 and HepG2 cells, respectively. HuH-7 and HepG2 cells also expressed the hepatic UGTs 1A3, 1A4, and 1A9. Neither cell line expressed UGT1A5 (Fig. 2C).

Effects of miR-491-3p on UGT1A Expression and Glucuronidation Activity.

To identify which hepatic UGT1A enzyme is regulated by miR-491-3p, HuH-7 cells were transiently transfected with miR-491-3p mimic, and UGT1A mRNA levels were assayed by qRT-PCR. MicroR-491-3p–transfected HuH-7 cells exhibited a significant (P < 0.001) 48% decrease in UGT1A1 mRNA levels as compared with scrambled miRNA-transfected control cells (Fig. 3A). In addition, UGTs 1A3 (P < 0.05) and 1A6 (P < 0.05) mRNA were also decreased in the presence of miR-491-3p, with UGT transcript levels reduced by 35 and 27%, respectively. Serving as a negative control, the levels of UGT2B7 mRNA, with its own unique 3′-UTR different from the common UGT1A 3′-UTR and not predicted to bind miR-491-3p, did not significantly change after miR-491-3p transfection (Fig. 3A). The expression of UGTs 1A4 and 1A9 were not significantly repressed after transfection of miR-491-3p (data not shown).

Fig. 3.
Fig. 3.

Quantitative analysis of UGT1A mRNA and enzymatic activity in HuH-7 cells. (A) mRNA levels of UGT1A1, 1A3, 1A6, and 2B7 in miR-491-3p transfected HuH-7 cells. 50 nM of scrambled miRNA or miR-491-3p mimic were transiently transfected into HuH-7 cells for 48 hours. The mRNA was isolated, and levels were compared with scrambled control (set at 1.0) for each gene. Genes were normalized to the RPLPO endogenous control gene using the 2−ΔΔCt method. UGT2B7 mRNA levels served as a negative control. (Inset) A representative semiquantitative electrophoresis gel of UGT1A1 mRNA. (B) UGT1A1 glucuronidation activity against raloxifene in HuH-7 cells. Ultra-pressure liquid chromatography–tandem mass spetrometry was used to quantify the formation of ral-6-gluc (top) and ral-4′-gluc (bottom) from raloxifene substrate in miRNA-transfected HuH-7 homogenates. (C) UGT2B7 protein levels and glucuronidation activity in the presence of 50 nM miR-491-3p miRNA mimic or scrambled control. Columns represent the mean ± S.E. of three independent experiments. *P < 0.05; **P < 0.01; ***P < 0.001.

Because UGT1A1 was the only UGT1A family member that was relatively well expressed in the miR-491-3p–expressing HuH-7 and HepG2 cell lines and because UGT1A1 exhibited the largest variation in UGT1A expression after miR-491-3p overexpression, the role of miR-491-3p in regulating UGT1A1 expression and activity was examined specifically. Raloxifene, a chemotherapeutic agent used in the prevention of breast cancer in postmenopausal women (Jordan et al., 2001), is exclusively metabolized by the UGT1A family of enzymes, with the major raloxifene glucuronide ral-6-gluc formed primarily by UGT1A1 in human liver, whereas both UGTs 1A1 and 1A9 contribute to hepatic ral-4′-gluc formation (Kemp et al., 2002; Mizuma, 2009; Trdan Lusin et al., 2011; Sun et al., 2013). Because UGT1A1 mRNA expression levels are ∼469- and ∼15-fold higher than UGT1A9 mRNA levels in HuH-7 and HepG2 cells, respectively (Fig. 2C), levels of raloxifene glucuronidation activity can serve as a specific marker for UGT1A1 activity and overall protein expression.

As a marker of UGT1A1 activity, the levels of ral-6-gluc and ral-4′-gluc formation were monitored in homogenates from HuH-7 cells transiently transfected with miR-491-3p mimic versus scrambled miRNA-transfected controls. In HuH-7 cell homogenates, both ral-6-gluc and ral-4′-gluc formation were significantly reduced by 40% (P < 0.01) and 38%, (P < 0.01), respectively, after transfection of miR-491-3p mimic as compared with scrambled miRNA-transfected controls (Fig. 3B). No change in UGT2B7 protein levels were observed as determined by Western blot analysis, and no difference in UGT2B7 glucuronidation activity was observed against epirubicin, a specific substrate for UGT2B7 (Innocenti et al., 2001; Zaya et al., 2006), in the presence of miR-491-3p mimic as compared with scrambled controls (Fig. 3C). Interestingly, no effect on UGT1A1 expression or activity was observed in HepG2 cells after transfection of miR-491-3p mimic as compared with scrambled miRNA-transfected controls (results not shown).

To further examine the endogenous effect of miR-491-3p on the expression and downstream activity of UGT1A1, HuH-7 cells were transiently transfected with 100 nM miR-491-3p inhibitor for 48 hours. Endogenous miR-491-3p levels were significantly (P < 0.01) reduced to below the level of detection in the presence of inhibitor as compared with scrambled miRNA inhibitor control (Fig. 4A). UGT1A1 mRNA levels were correspondingly significantly (P < 0.01) elevated (∼80%) in the presence of inhibitor as compared with scrambled control (Fig. 4B). This increase in UGT1A1 mRNA levels corresponded with significant increases in ral-6-gluc (50%; P < 0.05) and ral-4′-gluc (22%; P < 0.01) formation for HuH-7 homogenates with repressed miR-491-3p levels as compared with scrambled control inhibitor-treated cell homogenates (Fig. 4C). No differences in 1) UGT2B7 mRNA levels (Fig. 4B), 2) UGT2B7 protein levels (Fig. 4D), or 3) glucuronidation activity against epirubicin (Fig. 4D) were observed in HuH-7 homogenates with repressed miR-491-3p levels.

Fig. 4.
Fig. 4.

The effect of miR-491-3p inhibitor on UGT1A1 expression and activity in HuH-7 cells. (A) 100 nM of scrambled inhibitor control or miR-491-3p inhibitor were transiently transfected into HuH-7 cells for 48 hours. Quantitative analysis of endogenous mature miR-491-3p expression levels were measured following miRNA isolation from transfected HuH-7 cells using qRT-PCR with RNU6B as an endogenous control. (B) Quantification of UGT1A1 and UGT2B7 mRNA levels by qRT-PCR in HuH-7 cells in the presence of 100 nM miR-491-3p inhibitor or scrambled control. Quantification was performed with RPLPO serving as the endogenous expression control gene using the 2−ΔΔCt method. (Inset) A representative electrophoresis gel of the RT-PCR of UGT1A1 mRNA levels. (C) UGT glucuronidation activity against raloxifene in HuH-7 cells with inhibited miR-491-3p. Ultra-pressure liquid chromatography–tandem mass spetrometry was used to quantify the formation of ral-6-gluc (top) and ral-4′-gluc (bottom) from raloxifene substrate in miR-491-3p–inhibited HuH-7 cell homogenates. (D) UGT2B7 protein levels and glucuronidation activity in the presence of 100 nM miR-491-3p miRNA mimic or scrambled control. Columns represent the mean ± S.E. of three independent experiments. *P < 0.05; **P < 0.01.

A similar pattern was observed for the same miRNA inhibitor experiments performed in HepG2 cells. Endogenous miR-491-3p levels in HepG2 cells were significantly (P < 0.001) repressed in the presence of inhibitor when compared with scrambled control (Fig. 5A). Although UGT1A1 mRNA levels were not significantly altered in the presence of repressed miR-491-3p (Fig. 5B), miR-491-3p inhibition in HepG2 cells significantly increased ral-6-gluc and ral-4′-gluc formation by nearly 50% (P < 0.05) and 34% (P < 0.001), respectively (Fig. 5C). No difference in 1) UGT2B7 mRNA levels (Fig. 5B), 2) UGT2B7 protein levels (Fig. 5D), or 3) glucuronidation activity against epirubicin (Fig. 5D) were observed in HepG2 homogenates with repressed miR-491-3p levels.

Fig. 5.
Fig. 5.

The effect of miR-491-3p inhibitor on UGT1A1 expression and activity in HepG2 cells. (A) One hundred nM scrambled inhibitor control or miR-491-3p inhibitor were transiently transfected into HepG2 cells for 48 hours. Quantitative analysis of endogenous mature miR-491-3p expression levels were measured after miRNA isolation from transfected HepG2 cells using qRT-PCR and the RNU6B endogenous control. (B) Quantification of UGT1A1 and UGT2B7 mRNA levels by qRT-PCR in HepG2 cells in the presence of 100 nM miR-491-3p inhibitor or scrambled inhibitor control. Quantification was performed with RPLPO serving as the endogenous expression control gene using the 2−ΔΔCt method. (C) UGT glucuronidation activity against raloxifene in HepG2 cells with inhibited miR-491-3p. Ultra-pressure liquid chromatography–tandem mass spetrometry was used to quantify the formation of ral-6-gluc (top) and ral-4′-gluc (bottom) from raloxifene substrate in miR-491-3p- inhibited HepG2 cell homogenates. (D) UGT2B7 protein levels and glucuronidation activity in the presence of 100 nM miR-491-3p miRNA inhibitor or scrambled inhibitor control. Columns represent the mean ± S.E. of three independent experiments. *P < 0.05; ***P < 0.001.

Hepatic UGT versus miR-491-3p Expression in Human Liver Specimens.

To assess whether there was any functional effect of miR-491-3p expression on UGT1A expression, the expression of UGTs 1A1, 1A3, 1A6, and 1A9 was examined by qRT-PCR in a panel of normal human liver specimens. A significant inverse correlation was observed between UGT1A3 mRNA and miR-491-3p expression levels in 38 liver specimens that expressed UGT1A3 (Spearman r = −0.296; P < 0.05; Fig. 6A). A statistically significantly (P < 0.05) higher level of UGT1A3 expression was observed in miR-491-3p nonexpressers as compared with miR-491-3p expressers in these specimens (Fig. 6B). Similar to that observed for UGT1A3, a significant inverse correlation was also observed between UGT1A6 mRNA and miR-491-3p expression in the same liver specimens expressing UGT1A6 mRNA (n = 37; Spearman r = −0.487; P < 0.01; Fig. 6C). Also similar to that observed for UGT1A3, a statistically significantly higher (P < 0.05) level of UGT1A6 expression was observed in the miR-491-3p nonexpressers as compared with expressers (Fig. 6D). No correlation existed between UGT1A1 or UGT1A9 mRNA levels and miR-491-3p expression was observed in the same panel of liver specimens, and no difference in UGT1A1 or UGT1A9 mRNA levels was observed between miR-491-3p–expressing versus miR-491-3p–nonexpressing liver specimens (results not shown).

Fig. 6.
Fig. 6.

Expression of UGT1A3 and UGT1A6 mRNA versus miR-491-3p in human liver specimens. Levels of UGT1A3 (n = 38; A and B) and UGT1A6 (n = 37; C and D) mRNA versus miR-491-3p were quantified in normal liver samples via qRT-PCR and normalized to RPLPO and RNU6B, respectively. (A and C) Expression levels are shown relative to the lowest expressing sample (set to 1.0) for each gene, with those specimens with no miR-491-3p expression assigned a Ct value half that of the lowest miR-491-3p–expressing liver specimen. (A) UGT1A3 mRNA versus miR-491-3p expression was examined using the Spearman ranking method. (B) UGT1A3 mRNA expression levels were examined in liver specimens stratified by expression versus no expression of miR-491-3p. (C) UGT1A6 mRNA versus miR-491-3p expression was examined using the Spearman ranking method. (D) UGT1A6 mRNA expression levels were examined in liver specimens stratified by expression versus no expression of miR-491-3p. Dots represent the mean ± S.E. of three independent replicates.

Discussion

This is the first study to examine the regulation of the UGT family of enzymes by miRNA. In silico prediction modeling using both miRanda and TargetScan indicated that miR-491-3p was the top candidate to bind to the UGT1A 3′-UTR. Although miR-491-3p and two other top candidate miRNAs (miR-148a and miR-136) were tested using the luciferase reporter assay, only miR-491-3p repressed luciferase activity in vitro. These data were further confirmed in the UGT1A1-expressing HuH-7 cell line, as miR-491-3p overexpression was shown to significantly inhibit UGT1A1 expression and glucuronidation activity whereas inhibition of miR-491-3p was shown to significantly induce UGT1A1 expression and glucuronidation activity. In addition, deletion of the seed sequence within the 3′-UTR miR-491-3p MRE of the UGT1A1 gene resulted in an elimination of these effects, suggesting specificity of action.

The mRNA expression of UGTs 1A3 and 1A6 were significantly and inversely correlated with miR-491-3p expression levels in a panel of normal human liver specimens examined in this study. Furthermore, a significant difference in UGT1A3 and UGT1A6 mRNA expression was observed in liver specimens stratified by miR-491-3p expression versus miR-491-3p nonexpression. Together, these data strongly suggest that miR-491-3p plays a role in the regulation of UGT1A expression and activity and may contribute to the observed interindividual variability of several UGT1A genes in humans.

Although a significant increase in UGT1A1 expression and enzyme activity was observed in HepG2 cells after inhibition of miR-491-3p, no effect on UGT1A1 mRNA levels or activity was observed in HepG2 cells after miR-491-3p overexpression. This contrasts with the inhibitory effect observed for both UGT1A1 mRNA expression and enzyme activity after miR-491-3p overexpression in HuH-7 cells. This difference is likely due to the higher miR-491-3p:UGT1A1 expression ratio observed in HepG2 versus HuH-7 cells. This ratio difference is due mainly to the difference in expression of UGT1A1 in the two cell lines—UGT1A1 is expressed in HepG2 cells at levels well below (∼40-fold less) that observed in HuH-7 cells; there is a ∼5-fold difference in miR-491-3p expression levels between the two cell lines. In HepG2 cells where the miR-491-3p:UGT1A1 expression ratio is high, it is likely that further expression of miR-491-3p has little effect on UGT1A1 expression and subsequent activity—it may already be at saturating levels, a situation that is less likely to be observed for HuH-7 cells where the miR-491-3p:UGT1A1 expression ratio is significantly lower. When miR-491-3p is inhibited, a significant effect on UGT1A1 expression and subsequent activity would be expected in both cell lines.

Significant inverse correlations were observed between both UGT1A3 and UGT1A6 mRNA levels and miR-491-3p expression in normal human liver specimens, but no such association was observed for UGT1A1. Similarly, although a significant difference in UGT1A3 and UGT1A6 mRNA expression was observed in miR-491-3p expressing versus miR-491-3p–nonexpressing liver specimens in this study, no difference was observed for UGT1A1 expression. Additionally, no association was observed between UGT1A9 mRNA levels and miR-491-3p expression, and this is consistent with the fact that UGT1A9 was not shown to be affected by miR-491-3p in vitro. Previous studies have shown that UGTs 1A3 and 1A6 are expressed at levels up to 10- and 3-fold lower, respectively, than UGT1A1 in normal human liver (Izukawa et al., 2009; Ohno and Nakajin, 2009). Similar to that described above for cell lines, the miR-491-3p:UGT1A expression ratio is therefore higher for UGT1A3 and UGT1A6 than that for UGT1A1 in human liver, and changes in hepatic miR-491-3p expression may only significantly affect the lower mRNA expression of hepatic UGTs such as UGTs 1A3 and 1A6. These results do not rule out the possibility that miR-491-3p regulates UGT1A1 protein expression in human liver. Future studies are required to examine the impact of miR-491-3p expression on the enzymatic activity of UGT1A1 in liver specimens overexpressing miR-491-3p.

The levels of both UGT1A1 mRNA and enzyme activity were significantly altered in HuH-7 cells whereas in HepG2 cells only UGT1A1 enzyme activity was affected after overexpression of miR-491-3p. This may reflect different cell-specific mechanisms of miR-491-3p regulation of UGT1A1 protein expression. It is likely that the alterations in UGT1A1 activity observed in HuH-7 cells are via a reduction in mRNA levels, presumably through mRNA degradation. In HepG2 cells, with no observable reduction in UGT1A1 mRNA levels after miR-491-3p overexpression, the primary mechanism of protein repression could be through translational inhibition or relocation and sequestering of UGT1A1 mRNA into P-bodies, where mRNA can be stored with RISC for later usage or until degraded (Macfarlane and Murphy, 2010). Cell-specific miRNA regulation of protein expression is observed for other genes. For example, hepatocyte nuclear factor 4 α (HNF4α) is regulated in HepG2 cells by miR-24 and miR-34a. Mechanistically, miR-24 reduces HNF4α protein expression primarily via mRNA degradation whereas miR-34a primarily regulates HNF4α expression via translational inhibition (Takagi et al., 2010). In HuH-7 cells, peroxisome proliferator–activated receptor α is primarily regulated by miR-21 and miR-27b via translation inhibition (Kida et al., 2011), while in HepG2 cells peroxisome proliferator–activated receptor α is regulated by miR-141 via reduced mRNA expression levels (Hu et al., 2012). In both cases, however, the differential effect observed between cell lines was manifested by two different miRNA species. For UGT1A1, this differential cell line effect was linked to the same miRNA (miR-491-3p).

It is possible that the UGT1A 3′-UTR is not the sole region where miRNA may play a role in UGT1A regulation, as there are other genes that are regulated by miRNA predicted to bind within the 5′-untranslated or coding regions. For example, miR-24 regulates HNF4α protein expression primarily by binding within the mRNA coding region (Takagi et al., 2010). In the present study, no effect on UGT1A glucuronidation activity was observed when levels of miR-125a-3p were altered in either HepG2 or Huh7 cells (results not shown). As miR-125a-3p was the miRNA predicted by miRanda to bind within non–3′-UTR UGT1A sequences with the most efficacy, this suggests that miRNA regulation of UGT1A expression is primarily via binding to the UGT1A 3′-UTR. It is also possible that the prediction programs used in this analysis (miRanda and TargetScan) failed to adequately predict the binding of other miRNA within these UGT1A sequences. MicroRNA-pulldown techniques using the UGT1A 3′-UTR as a probe could be used to better assess this possibility (Roff et al., 2013).

The expression of the hepatic UGT1A4 and UGT1A9 mRNA levels were not significantly regulated by miR-491-3p overexpression in the present studies. This may be due to UGT1A mRNA secondary structure potentially interfering with miR-491-3p and/or protein binding. We analyzed the mRNA secondary structures for UGTs 1A1, 1A3, 1A4, 1A6, and 1A9 using the RNA Fold software (Gruber et al., 2008) to better assess this possibility (Supplemental Fig. 1). Although UGTs 1A3 and 1A4 had the most similar structure, each UGT1A isoform exhibited a unique secondary structure, and this may impact miR-491-3p binding to each isoform (Supplemental Fig. 1). The mRNA secondary structure can play an important role in miRNA regulation and targeting, as previous reports have identified polymorphisms within mRNA 3′-UTRs that alter miRNA regulation due to changes in mRNA secondary structure (Haas et al., 2012). Further studies will be necessary to better assess this possibility.

Although miR-491-3p was shown to be expressed hepatically, it was also shown to be expressed in several other normal human tissues including the colon, lung, trachea, and breast. In fact, the expression of miR-491-3p was higher in these tissues than was observed in the liver. Therefore, in addition to playing a potential role in whole-body metabolism via hepatic mechanisms, miR-491-3p could be playing an important role in regulating local metabolism within these extrahepatic tissues. Further studies are required examining the effects of miR-491-3p on UGT1A expression for nonhepatic human tissues.

In summary, we present the first evidence demonstrating that miRNA regulate phase II UGT metabolic enzymes. Results from both in vitro studies in cell lines and human tissue phenotype data demonstrate that UGTs including 1A3, 1A6, and possibly other UGT1A enzymes including UGT1A1, are regulated by miR-491-3p. Alterations in the expression of UGT1A enzymes by miR-491-3p may be an important mechanism controlling phase II metabolism and may play an important role in modulating drug response in humans.

Acknowledgments

The authors thank the Pennsylvania State University Nucleic Acids Research Core Facility for their support and help with sequencing and PCR expression analysis, and the authors acknowledge the assistance and help of the Penn State Hershey Cancer Institute Biorepository.

Authorship Contributions

Participated in research design: Dluzen, Sun, Salzberg, Jones, Bushey, Robertson, Lazarus.

Conducted experiments: Dluzen, Salzberg.

Contributed new reagents or analytic tools: Dluzen, Sun, Salzberg, Jones, Bushey, Robertson, Lazarus.

Performed data analysis: Dluzen, Sun, Salzberg, Jones, Bushey, Lazarus.

Wrote or contributed to the writing of the manuscript: Dluzen, Sun, Salzberg, Jones, Bushey, Robertson, Lazarus.

Footnotes

    • Received October 18, 2013.
    • Accepted January 3, 2014.

  • This work was funded in part by a student fellowship provided by AstraZeneca (to D.F.D.); grants from the National Institutes of Health National Institute of Dental and Craniofacial Research [Grant R01-DE13158] (to P.L.); the Pennsylvania Department of Health, Health Research Formula Funding Program [Grant 4100038714] (to P.L.); and the National Institutes of Health National Center for Research Resources and National Center for Advancing Translational Sciences [Grant UL1-TR000127].

  • dx.doi.org/10.1124/jpet.113.210658.

  • ↵Embedded ImageThis article has supplemental material available at jpet.aspetjournals.org.

Abbreviations

BLD
below the limit of detection
bp
base pair
Ct
cycle threshold
FBS
fetal bovine serum
HEK
human embryonic kidney
HNF4α
hepatocyte nuclear factor 4 α
miRNA
microRNA
PCR
polymerase chain reaction
MRE
miRNA recognition elements
qRT-PCR
quantitative real-time reverse-transcription polymerase chain reaction
ral-4′-gluc
raloxifene-4′-glucuronide
ral-6-gluc
raloxifene-6-glucuronide
RFLPO
ribosomal protein, large, P0
RNU6B
U6 small nuclear RNA 2
UGT
UDP-glucuronosyltransferase
RISC
RNA-induced silencing complex
UTR
untranslated region

References

  1. ↵
    1. Baek D,
    2. Villén J,
    3. Shin C,
    4. Camargo FD,
    5. Gygi SP, and
    6. Bartel DP
    (2008) The impact of microRNAs on protein output. Nature 455:64–71.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Balliet RM,
    2. Chen G,
    3. Gallagher CJ,
    4. Dellinger RW,
    5. Sun D, and
    6. Lazarus P
    (2009) Characterization of UGTs active against SAHA and association between SAHA glucuronidation activity phenotype with UGT genotype. Cancer Res 69:2981–2989.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Bartel DP
    (2009) MicroRNAs: target recognition and regulatory functions. Cell 136:215–233.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bélanger A,
    2. Pelletier G,
    3. Labrie F,
    4. Barbier O, and
    5. Chouinard S
    (2003) Inactivation of androgens by UDP-glucuronosyltransferase enzymes in humans. Trends Endocrinol Metab 14:473–479.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Betel D,
    2. Wilson M,
    3. Gabow A,
    4. Marks DS, and
    5. Sander C
    (2008) The microRNA.org resource: targets and expression. Nucleic Acids Res 36 (Database issue):D149–D153.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Bosma PJ,
    2. Seppen J,
    3. Goldhoorn B,
    4. Bakker C,
    5. Oude Elferink RP,
    6. Chowdhury JR,
    7. Chowdhury NR, and
    8. Jansen PL
    (1994) Bilirubin UDP-glucuronosyltransferase 1 is the only relevant bilirubin glucuronidating isoform in man. J Biol Chem 269:17960–17964.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Bushey RT and
    2. Lazarus P
    (2012) Identification and functional characterization of a novel UDP-glucuronosyltransferase 2A1 splice variant: potential importance in tobacco-related cancer susceptibility. J Pharmacol Exp Ther 343:712–724.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Carthew RW and
    2. Sontheimer EJ
    (2009) Origins and mechanisms of miRNAs and siRNAs. Cell 136:642–655.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Chen G,
    2. Blevins-Primeau AS,
    3. Dellinger RW,
    4. Muscat JE, and
    5. Lazarus P
    (2007) Glucuronidation of nicotine and cotinine by UGT2B10: loss of function by the UGT2B10 Codon 67 (Asp>Tyr) polymorphism. Cancer Res 67:9024–9029.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Chen G,
    2. Dellinger RW,
    3. Sun D,
    4. Spratt TE, and
    5. Lazarus P
    (2008) Glucuronidation of tobacco-specific nitrosamines by UGT2B10. Drug Metab Dispos 36:824–830.
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Chen G,
    2. Giambrone NE Jr.,
    3. Dluzen DF,
    4. Muscat JE,
    5. Berg A,
    6. Gallagher CJ, and
    7. Lazarus P
    (2010) Glucuronidation genotypes and nicotine metabolic phenotypes: importance of functional UGT2B10 and UGT2B17 polymorphisms. Cancer Res 70:7543–7552.
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Chen G,
    2. Giambrone NE, and
    3. Lazarus P
    (2012) Glucuronidation of trans-3′-hydroxycotinine by UGT2B17 and UGT2B10. Pharmacogenet Genomics 22:183–190.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Court MH,
    2. Duan SX,
    3. von Moltke LL,
    4. Greenblatt DJ,
    5. Patten CJ,
    6. Miners JO, and
    7. Mackenzie PI
    (2001) Interindividual variability in acetaminophen glucuronidation by human liver microsomes: identification of relevant acetaminophen UDP-glucuronosyltransferase isoforms. J Pharmacol Exp Ther 299:998–1006.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Dellinger RW,
    2. Chen G,
    3. Blevins-Primeau AS,
    4. Krzeminski J,
    5. Amin S, and
    6. Lazarus P
    (2007) Glucuronidation of PhIP and N-OH-PhIP by UDP-glucuronosyltransferase 1A10. Carcinogenesis 28:2412–2418.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Enright AJ,
    2. John B,
    3. Gaul U,
    4. Tuschl T,
    5. Sander C, and
    6. Marks DS
    (2003) MicroRNA targets in Drosophila. Genome Biol 5:R1.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Gallagher CJ,
    2. Balliet RM,
    3. Sun D,
    4. Chen G, and
    5. Lazarus P
    (2010) Sex differences in UDP-glucuronosyltransferase 2B17 expression and activity. Drug Metab Dispos 38:2204–2209.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Garcia DM,
    2. Baek D,
    3. Shin C,
    4. Bell GW,
    5. Grimson A, and
    6. Bartel DP
    (2011) Weak seed-pairing stability and high target-site abundance decrease the proficiency of lsy-6 and other microRNAs. Nat Struct Mol Biol 18:1139–1146.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Gruber AR,
    2. Lorenz R,
    3. Bernhart SH,
    4. Neuböck R, and
    5. Hofacker IL
    (2008) The Vienna RNA websuite. Nucleic Acids Res 36(Web Server issue):W70–74.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Guo H,
    2. Ingolia NT,
    3. Weissman JS, and
    4. Bartel DP
    (2010) Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature 466:835–840.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Haas U,
    2. Sczakiel G, and
    3. Laufer SD
    (2012) MicroRNA-mediated regulation of gene expression is affected by disease-associated SNPs within the 3′-UTR via altered RNA structure. RNA Biol 9:924–937.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Hatziapostolou M,
    2. Polytarchou C,
    3. Aggelidou E,
    4. Drakaki A,
    5. Poultsides GA,
    6. Jaeger SA,
    7. Ogata H,
    8. Karin M,
    9. Struhl K,
    10. Hadzopoulou-Cladaras M,
    11. et al.
    (2011) An HNF4α-miRNA inflammatory feedback circuit regulates hepatocellular oncogenesis. Cell 147:1233–1247.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Hoskins JM,
    2. Goldberg RM,
    3. Qu P,
    4. Ibrahim JG, and
    5. McLeod HL
    (2007) UGT1A1*28 genotype and irinotecan-induced neutropenia: dose matters. J Natl Cancer Inst 99:1290–1295.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Hu W,
    2. Wang X,
    3. Ding X,
    4. Li Y,
    5. Zhang X,
    6. Xie P,
    7. Yang J, and
    8. Wang S
    (2012) MicroRNA-141 represses HBV replication by targeting PPARA. PLoS ONE 7:e34165.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Innocenti F,
    2. Iyer L,
    3. Ramírez J,
    4. Green MD, and
    5. Ratain MJ
    (2001) Epirubicin glucuronidation is catalyzed by human UDP-glucuronosyltransferase 2B7. Drug Metab Dispos 29:686–692.
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Izukawa T,
    2. Nakajima M,
    3. Fujiwara R,
    4. Yamanaka H,
    5. Fukami T,
    6. Takamiya M,
    7. Aoki Y,
    8. Ikushiro S,
    9. Sakaki T, and
    10. Yokoi T
    (2009) Quantitative analysis of UDP-glucuronosyltransferase (UGT) 1A and UGT2B expression levels in human livers. Drug Metab Dispos 37:1759–1768.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. John B,
    2. Enright AJ,
    3. Aravin A,
    4. Tuschl T,
    5. Sander C, and
    6. Marks DS
    (2004) Human MicroRNA targets. PLoS Biol 2:e363.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Jones NR,
    2. Sun D,
    3. Freeman WM, and
    4. Lazarus P
    (2012) Quantification of Hepatic UDP glucuronosyltransferase 1A splice variant expression and correlation of UDP glucuronosyltransferase 1A1 variant expression with glucuronidation activity. J Pharmacol Exp Ther 342:720–729.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Jordan VC,
    2. Gapstur S, and
    3. Morrow M
    (2001) Selective estrogen receptor modulation and reduction in risk of breast cancer, osteoporosis, and coronary heart disease. J Natl Cancer Inst 93:1449–1457.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Kemp DC,
    2. Fan PW, and
    3. Stevens JC
    (2002) Characterization of raloxifene glucuronidation in vitro: contribution of intestinal metabolism to presystemic clearance. Drug Metab Dispos 30:694–700.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Khraiwesh B,
    2. Arif MA,
    3. Seumel GI,
    4. Ossowski S,
    5. Weigel D,
    6. Reski R, and
    7. Frank W
    (2010) Transcriptional control of gene expression by microRNAs. Cell 140:111–122.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Kida K,
    2. Nakajima M,
    3. Mohri T,
    4. Oda Y,
    5. Takagi S,
    6. Fukami T, and
    7. Yokoi T
    (2011) PPARα is regulated by miR-21 and miR-27b in human liver. Pharm Res 28:2467–2476.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Komagata S,
    2. Nakajima M,
    3. Takagi S,
    4. Mohri T,
    5. Taniya T, and
    6. Yokoi T
    (2009) Human CYP24 catalyzing the inactivation of calcitriol is post-transcriptionally regulated by miR-125b. Mol Pharmacol 76:702–709.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Lazarus P,
    2. Blevins-Primeau AS,
    3. Zheng Y, and
    4. Sun D
    (2009) Potential role of UGT pharmacogenetics in cancer treatment and prevention: focus on tamoxifen. Ann N Y Acad Sci 1155:99–111.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Lewis BP,
    2. Burge CB, and
    3. Bartel DP
    (2005) Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 120:15–20.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Liu DZ,
    2. Ander BP,
    3. Tian Y,
    4. Stamova B,
    5. Jickling GC,
    6. Davis RR, and
    7. Sharp FR
    (2012) Integrated analysis of mRNA and microRNA expression in mature neurons, neural progenitor cells and neuroblastoma cells. Gene 495:120–127.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Livak KJ and
    2. Schmittgen TD
    (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25:402–408.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Macfarlane LA and
    2. Murphy PR
    (2010) MicroRNA: Biogenesis, Function and Role in Cancer. Curr Genomics 11:537–561.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Mackenzie PI,
    2. Bock KW,
    3. Burchell B,
    4. Guillemette C,
    5. Ikushiro S,
    6. Iyanagi T,
    7. Miners JO,
    8. Owens IS, and
    9. Nebert DW
    (2005) Nomenclature update for the mammalian UDP glycosyltransferase (UGT) gene superfamily. Pharmacogenet Genomics 15:677–685.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Mackenzie PI,
    2. Hu DG, and
    3. Gardner-Stephen DA
    (2010) The regulation of UDP-glucuronosyltransferase genes by tissue-specific and ligand-activated transcription factors. Drug Metab Rev 42:99–109.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Mizuma T
    (2009) Intestinal glucuronidation metabolism may have a greater impact on oral bioavailability than hepatic glucuronidation metabolism in humans: a study with raloxifene, substrate for UGT1A1, 1A8, 1A9, and 1A10. Int J Pharm 378:140–141.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Mohri T,
    2. Nakajima M,
    3. Fukami T,
    4. Takamiya M,
    5. Aoki Y, and
    6. Yokoi T
    (2010) Human CYP2E1 is regulated by miR-378. Biochem Pharmacol 79:1045–1052.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Mukherji S,
    2. Ebert MS,
    3. Zheng GX,
    4. Tsang JS,
    5. Sharp PA, and
    6. van Oudenaarden A
    (2011) MicroRNAs can generate thresholds in target gene expression. Nat Genet 43:854–859.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Nagar S and
    2. Blanchard RL
    (2006) Pharmacogenetics of uridine diphosphoglucuronosyltransferase (UGT) 1A family members and its role in patient response to irinotecan. Drug Metab Rev 38:393–409.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Nagar S and
    2. Remmel RP
    (2006) Uridine diphosphoglucuronosyltransferase pharmacogenetics and cancer. Oncogene 25:1659–1672.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Nakamura A,
    2. Nakajima M,
    3. Yamanaka H,
    4. Fujiwara R, and
    5. Yokoi T
    (2008) Expression of UGT1A and UGT2B mRNA in human normal tissues and various cell lines. Drug Metab Dispos 36:1461–1464.
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Oda S,
    2. Nakajima M,
    3. Hatakeyama M,
    4. Fukami T, and
    5. Yokoi T
    (2012) Preparation of a specific monoclonal antibody against human UDP-glucuronosyltransferase (UGT) 1A9 and evaluation of UGT1A9 protein levels in human tissues. Drug Metab Dispos 40:1620–1627.
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Ohno S and
    2. Nakajin S
    (2009) Determination of mRNA expression of human UDP-glucuronosyltransferases and application for localization in various human tissues by real-time reverse transcriptase-polymerase chain reaction. Drug Metab Dispos 37:32–40.
    OpenUrlAbstract/FREE Full Text
  48. ↵
    1. Ohtsuki S,
    2. Schaefer O,
    3. Kawakami H,
    4. Inoue T,
    5. Liehner S,
    6. Saito A,
    7. Ishiguro N,
    8. Kishimoto W,
    9. Ludwig-Schwellinger E,
    10. Ebner T,
    11. et al.
    (2012) Simultaneous absolute protein quantification of transporters, cytochromes P450, and UDP-glucuronosyltransferases as a novel approach for the characterization of individual human liver: comparison with mRNA levels and activities. Drug Metab Dispos 40:83–92.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Pan YZ,
    2. Gao W, and
    3. Yu AM
    (2009) MicroRNAs regulate CYP3A4 expression via direct and indirect targeting. Drug Metab Dispos 37:2112–2117.
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Roff AN,
    2. Panganiban RP,
    3. Bond JS, and
    4. Ishmael FT
    (2013) Post-transcriptional regulation of meprin α by the RNA-binding proteins Hu antigen R (HuR) and tristetraprolin (TTP). J Biol Chem 288:4733–4743.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Rowland A,
    2. Miners JO, and
    3. Mackenzie PI
    (2013) The UDP-glucuronosyltransferases: their role in drug metabolism and detoxification. Int J Biochem Cell Biol 45:1121–1132.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Selbach M,
    2. Schwanhäusser B,
    3. Thierfelder N,
    4. Fang Z,
    5. Khanin R, and
    6. Rajewsky N
    (2008) Widespread changes in protein synthesis induced by microRNAs. Nature 455:58–63.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Sun D,
    2. Jones NR,
    3. Manni A, and
    4. Lazarus P
    (2013) Characterization of raloxifene glucuronidation: potential role of UGT1A8 genotype on raloxifene metabolism in vivo. Cancer Prev Res (Phila) 6:719–730.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Sun D,
    2. Sharma AK,
    3. Dellinger RW,
    4. Blevins-Primeau AS,
    5. Balliet RM,
    6. Chen G,
    7. Boyiri T,
    8. Amin S, and
    9. Lazarus P
    (2007) Glucuronidation of active tamoxifen metabolites by the human UDP glucuronosyltransferases. Drug Metab Dispos 35:2006–2014.
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Takagi S,
    2. Nakajima M,
    3. Kida K,
    4. Yamaura Y,
    5. Fukami T, and
    6. Yokoi T
    (2010) MicroRNAs regulate human hepatocyte nuclear factor 4α, modulating the expression of metabolic enzymes and cell cycle. J Biol Chem 285:4415–4422.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Trdan Lusin T,
    2. Trontelj J, and
    3. Mrhar A
    (2011) Raloxifene glucuronidation in human intestine, kidney, and liver microsomes and in human liver microsomes genotyped for the UGT1A1*28 polymorphism. Drug Metab Dispos 39:2347–2354.
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Tsuchiya Y,
    2. Nakajima M,
    3. Takagi S,
    4. Taniya T, and
    5. Yokoi T
    (2006) MicroRNA regulates the expression of human cytochrome P450 1B1. Cancer Res 66:9090–9098.
    OpenUrlAbstract/FREE Full Text
  58. ↵
    1. Valadi H,
    2. Ekström K,
    3. Bossios A,
    4. Sjöstrand M,
    5. Lee JJ, and
    6. Lötvall JO
    (2007) Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol 9:654–659.
    OpenUrlCrossRefPubMed
  59. ↵
    1. Wiener D,
    2. Fang JL,
    3. Dossett N, and
    4. Lazarus P
    (2004) Correlation between UDP-glucuronosyltransferase genotypes and 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone glucuronidation phenotype in human liver microsomes. Cancer Res 64:1190–1196.
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Zaya MJ,
    2. Hines RN, and
    3. Stevens JC
    (2006) Epirubicin glucuronidation and UGT2B7 developmental expression. Drug Metab Dispos 34:2097–2101.
    OpenUrlAbstract/FREE Full Text
Back to top

In this issue

Download PDF
Article Alerts
Email Article
Citation Tools

Research ArticleMetabolism, Transport, and Pharmacogenomics

Regulation of UGT1A1 Expression and Activity by miR-491-3p

Douglas F. Dluzen, Dongxiao Sun, Anna C. Salzberg, Nate Jones, Ryan T. Bushey, Gavin P. Robertson and Philip Lazarus
Journal of Pharmacology and Experimental Therapeutics March 1, 2014, 348 (3) 465-477; DOI: https://doi.org/10.1124/jpet.113.210658
Reddit logo Twitter logo Facebook logo Mendeley logo

Jump to section

Advertisement