Abstract
Transcripts of the UGT1A gene, encoding half of human UDP-glucuronosyltransferase (UGT) enzymes, undergo alternative splicing, resulting in active enzymes named isoforms 1 (i1s) and novel truncated isoforms 2 (i2s). Here, we investigated the effects of depleting endogenous i2 on drug response and attempted to unveil any additional biologic role(s) for the truncated novel UGT proteins. We used an integrated systems biology approach that combines RNA interference with unbiased global genomic and proteomic screens, and used HT115 colorectal cancer cells as a model. Consistent with previous evidence suggesting that i2s negatively regulate i1s through protein–protein interactions, i2-depleted cells were less sensitive to drug-induced cell death (IC50 of 0.45 ± 0.05 µM versus 0.22 ± 0.03 µM; P = 0.006), demonstrating that modulation of i2 levels meaningfully impacts drug bioavailability and cellular response. We also observed reduced production of reactive oxygen species by 30% (P < 0.05), and an enhanced expression (>1.2-fold; P < 0.05) of several proteins, such as hemoglobin α genes and superoxide dismutase 1, that have network functions associated with antioxidant properties. Interaction proteomics analysis of endogenous proteins from the cellular model, mainly in human intestine but also in kidney tissues, further uncovered interactions between i2s (but not i1s) and the antioxidant enzymes catalase and peroxiredoxin 1, which may influence antioxidant potential through sequestration of these novel partners. Our findings demonstrate for the first time dual roles for i2s in the cellular defense system as endogenous regulators of drug response as well as in oxidative stress.
Introduction
Conjugative cellular metabolism mediated by UDP-glucuronosyltransferases (UGTs) has a clear effect on xenobiotic and endobiotic bioavailability as well as an influence on drug response, treatment outcomes, and disease susceptibility and progression (Guillemette et al., 2010). Alternative splicing is one of the most significant components of the functional complexity of human UGTs. UGT1A is an ideal, yet complex, example of this protein repertoire–expanding mechanism because it gives rise to at least 27 unique mRNAs and 18 functional proteins (Bellemare et al., 2010b). These numerous alternative UGT1A transcripts are created by a series of splicing events involving 13 alternative first exons (1A1 through 1A13) and alternative terminal fifth exons (5a and 5b) (Girard et al., 2007). The alternative usage of first exons and the terminal exon 5a lead to classic variant mRNAs (v1) encoding nine enzymatically active UGT1A enzymes (isoforms 1, or i1s) with a distinct substrate specificity and tissue expression profile. The inclusion of the terminal exon 5b yields nine distinct UGT1A_v2 mRNAs (containing exon 5b but not 5a) and nine distinct UGT1A_v3 mRNAs (containing exons 5a and 5b), each having an identical open reading frame encoding nine shorter polypeptides named UGT1A isoforms 2 (i2s).
The i2 and i1 proteins have similar expression patterns in human tissues, are highly abundant in the gastrointestinal tract, and colocalize in the endoplasmic reticulum (ER); however, i2s are devoid of transferase activity (Lévesque et al., 2007; Bellemare et al., 2011). Further functional assays with human cells coexpressing both i1s and i2s clearly demonstrated a significant decrease in cellular glucuronidation activity for many different UGT1A substrates (Girard et al., 2007; Bellemare et al., 2010b). It was previously shown that transient repression of i2 by small interfering RNAs leads to an increased glucuronidation activity in human colon and liver cells, revealing a novel regulatory mechanism by i2s (Bellemare et al., 2010c; Jones et al., 2012). UGT oligomerization is likely the underlying mechanism for the negative modulation of transferase activity by splice variant i2. Indeed, several groups demonstrated that UGTs act as dimers/tetramers (Fujiwara et al., 2007; Operaña and Tukey, 2007), and experimental evidence demonstrated a direct interaction between i1–i2 proteins (Bellemare et al., 2010a,b).
We proposed that the novel i2 splice products may be biologically relevant with physiologic and pharmacological influences based on the following: 1) their wide expression profile (Girard et al., 2007; Bellemare et al., 2011); 2) their substantial expression levels particularly in the gastrointestinal tract and kidneys (Girard et al., 2007; Lévesque et al., 2007; Bellemare et al., 2011); 3) their often concurrent expression with the UGT1A_i1 enzymes in the same healthy tissue structures and cell types (Bellemare et al., 2011); 4) an evident lack of correlation between mRNA expression levels, UGT1A protein expression, and transferase activity (Ohtsuki et al., 2012), which might be explained by i2 alternative splicing events recently reported (Girard et al., 2007; Lévesque et al., 2007); 5) the modulatory function of i2 on cellular UGT activity in vitro that is comparable with the effect observed for genetic polymorphisms having an influence on drug response in vivo (Iyer et al., 2002); and 6) the observation that i2s do not reside exclusively in the ER (where i1s localize) and in fact can be detected in the cytoplasm (Lévesque et al., 2007), raising the possibility that i2s may have additional cellular functions perhaps independent of their role in drug metabolism.
To test this hypothesis, we developed a stable knockdown (KD) model of i2s using short hairpin RNA (shRNA) in HT115 colon cancer cells, which express both classes of i1 and i2 with an almost equivalent level of v2/v3 transcripts encoding i2 isoforms relative to v1 levels encoding UGT1A enzymes (0.8:1) (Bellemare et al., 2010c). We used, as a pharmacological model compound, the active metabolite of irinotecan (CPT-11), namely SN-38 (7-ethyl-10-hydroxycamptothecin), which is almost exclusively metabolized by the UGT1A pathway and is one of the most common cytotoxic agents used to treat metastatic colorectal cancer. The interindividual variability of the UGT1A pathway induced by genetic variants such as the UGT1A1*28 allele has well defined clinical relevance by considerably affecting the glucuronidation of SN-38 and risk for drug-induced adverse effects in cancer patients (Iyer et al., 2002). Using novel and unbiased approaches, we studied the impact of i2 depletion on cell metabolism and demonstrated that KD of i2 increases drug metabolism and decreases colon cancer cell susceptibility to the pharmacological compound, in addition to significantly affecting cell exposition to intracellular oxidative stress, demonstrating novel biologic function(s) for these spliced proteins.
Materials and Methods
Depletion of i2s by Exon 5b shRNA in Human Cells and Cellular Metabolism.
Primers for shRNAmir specific to exon 5b of UGT1A (Supplemental Table 1), were designed using RNAi Central (http://cancan.cshl.edu/RNAi_central/). Additional details are provided in the Supplemental Experimental Procedures. Human embryonic kidney 293T cells were transfected with pGIPZ (shRNAmir nontarget control or targeting exon 5b) and Trans-Lentiviral Packaging Mix using Arrest-In from the GIPZ system (Open Biosystem, Ottawa, ON, Canada). A multiplicity of infection of 1 was utilized to infect HT115 cells (European Collection of Cell Cultures, Salisbury, UK). Isoform-specific quantitative polymerase chain reaction strategies were used to confirm v2/v3 mRNA depletion (Bellemare et al., 2010c). Enzymatic metabolism assays for SN-38 glucuronidation were performed using 50 µg cell homogenates incubated for 60 minutes at 37°C with 200 µM SN-38 and rates of SN-38 glucuronide formation were measured by mass spectrometry (MS) (Gagné et al., 2002).
Cell Viability Assay and Intracellular Reactive Oxygen Species Detection.
HT115 cells were treated with increasing concentrations of SN-38 and AlamarBlue (Life Technologies, Burlington, ON, Canada) was added to assess cell viability. Fluorescence was measured (565/585 nm). IC50 values were calculated using Prism 5 (GraphPad Software, Inc., La Jolla, CA). To assess reactive oxygen species (ROS) levels, HT115 cells were treated with SN-38 for 48 hours and incubated with 5 µM CellROX Deep Red Reagent (Life Technologies) (30 minutes at 37°C). The fluorescence of labeled cells was measured by flow cytometry (BD FACSCanto II; BD Biosciences, San Jose, CA).
Global Gene Expression Analysis.
Whole-genome gene expression was assessed with total RNA from HT115 cells (KD versus reference). Total RNA was amplified, labeled, and hybridized to HumanHT-12 v4 Expression BeadChips (Illumina, San Diego, CA) at the Genome Quebec Innovation Centre at McGill University (Montreal, QC, Canada). Data were analyzed using FlexArray 1.6.1 software, and the Lumi algorithm was used to preprocess microarray data by variance-stabilized transformation and quantile normalization and used the Significance Analysis of Microarrays algorithm (Du et al., 2008). Differentially expressed genes were then subjected to functional clustering analysis using the Ingenuity Pathway Analysis system (Ingenuity, Redwood City, CA). Microarray data were submitted to Gene Expression Omnibus (accession number GSE48118).
Global Proteomic Analysis Using Stable Isotope Labeling of Amino Acids in Cell Culture–MS.
Detailed procedures of stable isotope labeling of amino acids in cell culture (SILAC) are provided in the Supplemental Experimental Procedures. HT115 were double SILAC labeled using stable isotopes of lysine and arginine (13C6 Lys and 13C6 15N4 Arg; Cambridge Isotope Laboratories, St. Leonard, QC, Canada). Equal amounts of light (nontarget) and heavy (KD) cells were mixed, washed, and lysed [150 mM NaCl, 50 mM Tris-HCl pH 7.4, 0.3% sodium deoxycholate, 1% (w/v) Igepal, 1 mM EDTA, complete protease inhibitor]. Trypsin digested peptides were analyzed using isoelectric focusing as described (Gagné et al., 2012) and identified using a TripleTOF 5600 (ABSciex, Concord, ON, Canada). Proteins were identified/quantified using Protein Pilot v4 software (Paragon and Progroup algorithms; Shilov/ABSciex) and searched using MASCOT (Matrix Science, Boston, MA) in the human protein database (Uniref May 2012). Data were submitted to the Proteomics Identifications database (accession number PXD000295).
Global Analysis of Protein–Protein Interactions.
Immunoprecipitation of UGT1A isoforms and protein complexes was carried out on S9 fractions of human kidney or intestine (Xenotech, Lenexa, KS) and HT115 cell lysates. Immunoprecipitation was performed on 1 mg total protein lysate (Supplemental Experimental Procedures) with 4 µg antibody (anti-UGT1A_i1 #9348 or _i2 #4863) and with G protein–coated magnetic beads (Dynabeads; Life Technologies) for 15 hours at 4°C. Protein identification was established using a TripleTOF 5600 after on-beads trypsin digestion. The observed endogenous interactions were confirmed by coimmunoprecipitation experiments with HT115 lysates. Protein interaction complexes were fixed using 0.125% formaldehyde. Cells were lysed [150 mM NaCl, 1% (w/v) Igepal, 1 mM dithiothreitol, 175 mM Tris-HCl pH 7.4, complete protease inhibitor] and immunoprecipitation was performed on supernatants using 200 µl beads (mix with 8 µg anti-UGT1A_i1 or i2). Protein complexes were eluted and subjected to SDS-PAGE, and the presence of interacting partners was revealed on a nitrocellulose membrane with a specific antibody.
Statistical Analysis.
All analyses were performed in triplicate, and data were derived from at least three independent experiments unless otherwise specified. Results are expressed as the mean ± S.D. Statistical analysis was performed with two-tailed t tests, with P < 0.05 considered statistically significant. Statistical analysis for microarrays was performed with the Significance Analysis of Microarrays algorithm, and P < 0.01 was considered statistically significant.
Results
Depletion of Endogenous i2 Splice Forms Modifies Cell Metabolism and Drug Response.
We developed a stable KD cellular model based on an i2-specific shRNAmir cloned into a lentiviral vector designed to deplete endogenous v2/v3 transcripts encoding i2 proteins in the HT115 colon cancer cell model. Eight shRNAmirs were tested (Fig. 1A; Supplemental Table 1) and the clone that provided the best KD was used. The shRNA p27 (start position at nucleotide 27 of exon 5b) led to a 49% decrease (P = 0.005) in v2/v3 mRNA levels without significantly affecting v1 mRNA levels (Fig. 1B). Depletion of i2 was thus specific and stable (from 20 to 200 days postinfection). Compared with nontarget shRNA-infected HT115 cells (corresponding to the reference cell line), drug glucuronidation activity increased 126% (P = 0.008) for i2-depleted cells (HT115 KD) (Fig. 1B). This increase in drug inactivation correlated with enhanced viability of cells, that is, the IC50 value of SN-38 for HT115 KD cells (0.45 ± 0.05 µM) was higher than that for the reference (0.22 ± 0.03 µM; P = 0.006) (Fig. 1C). This increased viability induced by i2 KD was particularly significant between 0.08 and 0.8 µM of SN-38, which correspond to doses required to achieve a 50% inhibition in cell growth.
Global Gene Expression Profiling of Functional KD of UGT1A_i2 and Altered Production of ROS.
Compared with the reference cells, microarray gene expression profiling of HT115 KD cells revealed significant modulation of 207 and 586 genes (fold change ≥1.2; P < 0.01), under standard culture conditions and upon treatment with SN-38 (0.5 µM, 48 hours), respectively. Figure 2 shows heatmap data for the 50 most differentially-expressed genes in cells treated or not with SN-38. Among these genes, hemoglobin α genes (HBA1/HBA2) were significantly upregulated (P < 0.01) (Supplemental Table 2). HBA1/HBA2 expression was confirmed in peripheral blood cells from few healthy donors and we also demonstrated their expression in tissues from the gastrointestinal tract, including sections of the intestine, healthy colon tissue, and colon tumors, as well as in HT115 cells (Fig. 3A). Compared with the reference cell line, HT115 KD cells further displayed upregulated expression of HBA2 under basal conditions (1.3-fold; P = 0.05) and of HBA1 after drug exposure (5.5-fold; P = 0.01), which was subsequently confirmed at the protein level (Fig. 3). Because the function of HBA proteins has been linked to cellular antioxidant potential (Widmer et al., 2010), we evaluated whether depletion of i2s influences cellular oxidative stress response by quantification of intracellular ROS. HT115 KD cells had greater antioxidant potential compared with the reference cell line, as ROS formation was reduced by >30% (P < 0.05) after SN-38 treatment (Fig. 3C), consistent with previous reports indicating that SN-38 induces cellular oxidative stress (Conklin, 2004; Timur et al., 2005).
Global Quantitative Protein Expression Profiling Associated with Functional KD of i2.
We used the unbiased SILAC-MS approach to quantify altered protein expression upon i2 depletion (Fig. 4A). We first used lysates of cells from SILAC to demonstrate essentially full incorporation (98.3%) of targeted isotopic labels. We next conducted a global MS study to quantitatively compare the proteomes of SILAC-labeled HT115 KD cells and unlabeled HT115 cells grown without drug treatment. We detected a significant increased expression for 11 proteins (>1.2-fold; P < 0.05) (Fig. 4B). We then applied a reliable bioinformatics approach to establish potential functional network based on the regulated proteins (Supplemental Fig. 1) and reveal that almost all proteins we identified are related to the intrinsic antioxidant defense network, including superoxide dismutase 1 (SOD1) and others (Table 1). Supplemental Figure 2 presents a representative MS spectrum of SOD1 peptides. Western blotting confirmed that SOD1 levels increased in HT115 KD cells (Fig. 4C).
Interaction Proteomics Experiments Reveal Antioxidant Enzymes Catalase and Peroxiredoxin 1 as Specific Partners of i2 Isoforms in Human Tissues and HT115 Cells.
Because i1 enzymes and i2 proteins function as oligomers, we tested the hypothesis that i2s interact with proteins other than i1s and specific partners not interacting with i1 enzymes, some of which may participate in oxidative stress pathways. We used immunoprecipitation coupled to MS using specific antibodies for each class of UGT1A isoforms to depict i1- and i2-associated interactomes (Fig. 5A), and observations were based on three independent immunoprecipitation–tandem MS experiments for each antibody. As specific partners of i2, notably in the intestine, we identified for the first time the enzyme catalase, which has the highest turnover number of all antioxidant enzymes (Valko et al., 2006), as well as the antioxidant enzyme peroxiredoxin 1 (Prdx1), which reduces hydrogen peroxide and alkyl hydroperoxides (Imlay, 2008). In i2 immunoprecipitates from the intestine, we detected several unique peptides specific to catalase and Prdx1 that were not identified in anti-i1 immunoprecipitates (Supplemental Table 3). Furthermore, coimmunoprecipitation using specific antibodies directed against catalase and Prdx1 confirmed their association with i2s (but not i1s) (Fig. 5B).
Discussion
UGTs are major mediators in conjugative metabolism, and variations of this pathway have been clearly shown to influence drug response and treatment outcomes (Lazarus et al., 2009; Guillemette et al., 2010). In keeping with this view, we sought to assess whether modulation of endogenous levels of regulatory i2 proteins, through RNA interference-mediated repression of i2 spliced isoforms, would impact cellular metabolism and pharmacological response to an anticancer agent (SN-38) primarily metabolized by UGT1A enzymes. A partial depletion of i2 isoforms was achieved and sufficient to significantly enhance glucuronidation of the drug. Furthermore, i2-depleted cells were also more viable than the reference cell line upon drug treatment and displayed increased resistance to this anticancer agent that cannot be attributed to a direct increase of active i1s. Because we previously showed that i2s act to repress glucuronidation activity through protein–protein interactions with i1s (Bellemare et al., 2010b), this partial depletion in i2 splice forms was expected to lead to increased availability of the active i1 enzymes, with a consequent overall increase in cellular UGT1A-mediated glucuronidation and a meaningful impact on the pharmacological response to a drug such as an improvement in cell viability. Our data are clearly supportive of a negative regulatory role for i2s upon drug glucuronidation and have the potential to affect the pharmacological response to drug and likely contribute to the observed interindividual variability in drug metabolism and response.
Our work also provides the first evidence that i2s participate in the regulation of oxygen and ROS metabolic processes through direct protein–protein interactions with antioxidant enzymes and by affecting their expression at the mRNA and protein levels. Microarrays revealed upregulated expression of HBA1 and HBA2 upon i2 depletion. HBA products are well known for their oxygen transporter function, but they also have an antioxidant function as they bind cellular ROS (Reeder, 2010; Widmer et al., 2010). These genes are primarily expressed in peripheral mononuclear blood cells (Haas et al., 2011) and, as demonstrated in our present work, are also expressed throughout the gastrointestinal tract with relatively higher expression in duodenum and jejunum, consistent with their role as antioxidant enzymes. Our data also demonstrate their expression in colon tumors and in the related HT115 cellular model. The increased expression of HBA in i2-depleted cells was further associated with a significant reduction in intracellular ROS formation, thus supporting a potential link between i2s and the oxidative stress pathway. Accumulation of ROS is a vital signal that initiates diverse cellular processes, and it is also plausible that altered ROS homeostasis might provide a critical signal for inducing antioxidant pathways, such as HBA. Likewise, a change in intracellular ROS levels similar to the one observed upon i2 depletion was previously shown to influence cellular autophagy in U2OS osteosarcoma cells (Bensaad et al., 2009). In addition, data from SILAC-MS also pointed to an enhanced expression of several other proteins that have network functions associated with antioxidant properties. This is in line with an interrelation among factors of the cellular defense system that is not only composed of phase I and II detoxification enzymes, such as UGTs, but of many of the proteins shown to be affected by i2 depletion, including SOD1 and others (Table 1). Defining the exact mechanisms by which these changes occur requires further experiments.
By applying an unbiased MS-based immunoprecipitation proteomics analysis to HT115 cells and human tissues, our experiments reveal for the first time a direct protein–protein interaction between i2s and antioxidant enzymes, such as catalase and Prdx1, further confirmed by coimmunoprecipitation. Whereas most protein complexes are usually isolated through overexpression of a tagged exogenous protein within cells, the strength of our approach relies on targeting endogenous proteins of human intestine and kidney where i2 expression was previously demonstrated (Lévesque et al., 2007; Bellemare et al., 2011). In addition, we show an i2-specific interaction with catalase and Prdx1, especially in the intestine. Both the catalase and Prdx1 pathways are involved in the elimination and modulation of local H2O2 levels that play a pivotal signaling role in several cellular processes. This previously unknown capacity of i2 proteins to interact with antioxidant proteins evokes their potential to regulate the activities of these enzymes. Both catalase and Prdx1 must form a tetramer to achieve activation, and thus it is possible that binding of i2s to these enzymes may serve as a way of inhibiting their activities by perturbing their oligomerization. Clearly such a property of i2s is biologically relevant, as we have shown that a portion of i2s is cytoplasmic in contrast to i1s that are resident of the ER (Lévesque et al., 2007; Bellemare et al., 2011). The cytoplasmic distribution of Prdx1, catalase, and i2 isoforms is therefore consistent with i2s interfering with their oligomerization and the biologic functions of these ROS scavengers (Leon et al., 2006; Rho et al., 2008).
We conclude that, in addition to regulating drug metabolism and response (Fig. 6), the endogenous levels of i2 splice isoforms may contribute to controlling oxidative stress and consequently numerous signaling pathways, thereby affecting several cellular phenotypes.
Acknowledgments
The authors thank Patrick Caron for MS analysis of drug metabolism and the proteomics platform of the Quebec Genomic Center for protein-related MS analyses.
Authorship Contributions
Participated in research design: Rouleau, Bellemare, Guillemette.
Conducted experiments: Rouleau, Roberge.
Performed data analysis: Rouleau, Roberge, Guillemette.
Wrote or contributed to the writing of the manuscript: Rouleau, Roberge, Guillemette.
Footnotes
- Received August 20, 2013.
- Accepted October 18, 2013.
This research was supported by the Canadian Institutes of Health Research (CIHR) [Grant MOP-84223/111086]; the National Sciences and Engineering Research Council of Canada [Grant CG086976]; and the Canada Research Chair Program [(to C.G.)]. M.R., J.R., and J.B. received a CIHR Frederick Banting and Charles Best Canada Graduate Scholarship studentship award, and J.R. received a studentship award from the Fonds de la recherche en santé du Québec. C.G. holds a Tier II Canada Research Chair in Pharmacogenomics.
↵This article has supplemental material available at molpharm.aspetjournals.org.
Abbreviations
- ER
- endoplasmic reticulum
- HBA
- hemoglobin α
- i1
- isoform 1
- i2
- isoform 2
- KD
- knockdown
- MS
- mass spectrometry
- Prdx1
- Peroxiredoxin 1
- ROS
- reactive oxygen species
- shRNA
- short hairpin RNA
- SILAC
- stable isotope labeling of amino acids in cell culture
- SN-38
- 7-ethyl-10-hydroxycamptothecin
- SOD1
- superoxide dismutase 1
- UGT
- UDP-glucuronosyltransferase
- Copyright © 2013 by The American Society for Pharmacology and Experimental Therapeutics