Abstract
ATP-binding cassette (ABC) transporters represent a large group of efflux pumps that are strongly involved in the pharmacokinetics of various drugs and nutrient distribution. It was recently shown that micro-RNAs (miRNAs) may significantly alter their expression as proven, e.g., for miR-379 and ABCC2. However, alternative mRNA polyadenylation may result in expression of 3′-untranslated regions (3′-UTRs) with varying lengths. Thus, length variants may result in presence or absence of miRNA binding sites for regulatory miRNAs with consequences on posttranscriptional control. In the present study, we report on 3′-UTR variants of ABCC1, ABCC2, and ABCC3 mRNA. Applying in vitro luciferase reporter gene assays, we show that expression of short ABCC2 3′-UTR variants leads to a significant loss of miR-379/ABCC2 interaction and subsequent upregulation of ABCC2 expression. Furthermore, we show that expression of ABCC2 3′-UTR lengths varies significantly between human healthy tissues but is not directly correlated to the respective protein level in vivo. In conclusion, the presence of altered 3′-UTR lengths in ABC transporters could lead to functional consequences regarding posttranscriptional gene expression, potentially regulated by alternative polyadenylation. Hence, 3′-UTR length variability may be considered as a further mechanism contributing to variability of ABCC transporter expression and subsequent drug variation in drug response.
SIGNIFICANCE STATEMENT micro-RNA (miRNA) binding to 3′-untranslated region (3′-UTR) plays an important role in the control of ATP-binding cassette (ABC)-transporter mRNA degradation and translation into proteins. We disclosed various 3′-UTR length variants of ABCC1, C2, and C3 mRNA, with loss of mRNA seed regions partly leading to varying and tissue-dependent interaction with miRNAs, as proven by reporter gene assays. Alternative 3′-UTR lengths may contribute to variable ABCC transporter expression and potentially explains inconsistent findings in miRNA studies.
Introduction
ATP-binding cassette (ABC) transporters mediate active transport of diverse endogenous compounds, xenobiotics, and drugs across barriers, thereby contributing to endogenous distribution of nutrients, detoxification, and drug elimination (Rees et al., 2009). Having 13 members, the C-subfamily is the largest among the seven ABC-transporter subfamilies (Moitra and Dean, 2011). Within the ABCC subfamily, mainly ABCC1 [multidrug resistance-associated protein (MRP) 1], ABCC2 (MRP2), and ABCC3 (MRP3) are involved in the bioavailability of various drugs, influencing their absorption, distribution, and elimination or restricting the permeability of blood-tissue barriers (Kunická and Souček, 2014; van der Schoor et al., 2015; Chen et al., 2016). Next to their contribution to drug response, ABCC transporters are involved in the etiology of several human pathologies. ABCC1, transporting conjugated and unconjugated organic anions, plays a role in immunologic and cardiovascular diseases as well as neurologic disorders and tumor progression (Cole, 2014). ABCC2 also mediates the transport of various conjugated organic anions and is mainly responsible for the excretion of bilirubin from the bile. A genetic defect in ABCC2 is associated with Dubin-Johnson syndrome, a recessively inherited disorder characterized by conjugated hyperbilirubinemia (van der Schoor et al., 2015). Similar to ABCC2, ABCC3 transports bile acids, organic anions, and numerous xenobiotics, including anticancer drugs, and acts as an alternative transporter for the export of bile acids and glucuronides from cholestatic hepatocytes. Because of the importance of these transporters in the pharmacokinetics of drugs and transport of endogenous organic anions, increased knowledge of regulatory mechanisms is of high interest.
On the transcriptional level, the genes of ABCC transporters are regulated by nuclear receptor signaling as a consequence of xenobiotic sensing and/or hormonal regulation (Kast et al., 2002; Miller, 2015). Moreover, hereditary genetic variants significantly affect transcriptional and posttranscriptional regulation and function, as shown for ABCC2 haplotypes (Laechelt et al., 2011). In addition, ABCC transporters underlie extensive posttranscriptional regulation through interaction with micro-RNAs (miRNAs) (Haenisch et al., 2014).
miRNAs are small noncoding RNAs, which act through RNA interference by forming imperfect hybrids with the 3′-UTR of their target mRNAs, leading to mRNA degradation or translation inhibition (Ambros, 2001). Various miRNAs were reported to bind to the 3′-UTRs of ABCC mRNAs, currently nine for ABCC1 (Liang et al., 2010; Pan et al., 2013; Kunická and Souček, 2014; Liu et al., 2015; Ma et al., 2015; Pei et al., 2016; Sun et al., 2017; Hu et al., 2018; Li et al., 2018a,b), four for ABCC2 (Haenisch et al., 2011; Zhan et al., 2013; He et al., 2017; Tian et al., 2017), and two for ABCC3 (Bruckmueller et al., 2017; An et al., 2018). Functional confirmation of miRNA binding was predominantly executed by reporter gene assays, i.e., binding of miRNA-379 to the ABCC2 3′-UTR (Haenisch et al., 2011). However, alternative polyadenylation was not considered in these studies.
Alternative polyadenylation is an RNA-processing mechanism generating mRNAs with distinct 3′-UTR lengths. Being widespread across all eukaryotic species, it is considered as a major mechanism of gene regulation (Tian and Manley, 2017). About 50%–75% of all human genes underlie alternative polyadenylation, resulting in multiple possible poly-A sites of one respective mRNA (Tian et al., 2005; Shi, 2012). Moreover, more than 50% of conserved miRNA target sites in the 3′-UTRs are affected by alternative polyadenylation (Sandberg et al., 2008). Furthermore, alternative polyadenylation is tissue specific and an important determinant for protein synthesis. Thereby, it plays a role in cell proliferation and differentiation as well as cancer cell progression (Mayr and Bartel, 2009; Mayr, 2016). Several studies have observed that ABC transporters are substantially affected by alternative polyadenylation. For ABCB1 (P-glycoprotein) and ABCG2 (breast cancer related protein), it was shown that alternative polyadenylation influences their posttranscriptional regulation and protein synthesis, subsequently contributing to pharmacoresistances of cancer cells (To et al., 2009; Bruhn et al., 2016). Despite their importance in drug disposition, nutrient transport, and health, the extent of alternative polyadenylation in ABCC transporters remains unknown. The present study aimed to determine alternative polyadenylation of the ABC transporters ABCC1, ABCC2, and ABCC3 and to investigate its potential regulatory implications or if 3′-UTR shortening leads to loss of miRNA binding sites and respective miRNA binding. Furthermore, the expression of different ABCC2 3′-UTR variants in human liver, colon, and gall bladder tissues were analyzed and compared with ABCC2 protein levels.
Materials and Methods
Cell Cultivation.
HepG2, SK-Hep-1, and Caco-2 cells were purchased from the German Collection of Microorganisms and Cell Cultures (DSMZ, Braunschweig, Germany). HepG2 cells were grown in RPMI1640 medium (PAA, Pasching, Germany) supplemented with 10% v/v heat-inactivated FBS (PAA) and 1% v/v penicillin (10,000 U/ml)/streptomycin (10 mg/ml) (P/S; PAA), SK-Hep-1 cells were grown in RPMI1640 medium supplemented with 20% v/v FBS and 1% v/v P/S, and Caco-2 cells were grown in DMEM medium (with L-glutamine and 4.5 g/l glucose) (PAA) supplemented with 10% v/v FBS enriched with 1% v/v nonessential amino acids (PAA) and 1% v/v P/S. The medium was replaced and cells were passaged every 2 or 3 days, depending on the cell density. For RNA extraction, 5 × 106 cells were solubilized in 600 µl RLT buffer (Qiagen, Hilden, Germany) supplemented with 1% v/v β-mercaptoethanol (Roth, Karlsruhe, Germany) and homogenized using a rotor-stator homogenizer (Polytron PT 3000; Kinematica AG, Littau, Switzerland). The RNA was extracted using the RNeasy Mini Kit (Qiagen) according to the manufacturer’s specifications. Prior to reverse transcription, residual DNA was digested using the Turbo DNA-free Kit (Ambion, Austin) according to the manufacturer’s specifications.
Human Tissue Samples.
Human tissue samples (liver, colon, and gall bladder) were obtained from patients undergoing tumor surgery at the University Hospital Schleswig-Holstein, Campus Kiel (Table 1). Patients gave their written informed consent, and the study was approved by the ethics committee of the Medical Faculty of the University of Kiel. Specimens of nonmalignant peritumoral tissue were obtained from the resectate after verification by a pathologist and were stored overnight at 4°C in RNAlater reagent (Qiagen), prior to long-term storage at −80°C. Tissue RNA was extracted using the Precellys 24 Tissue Homogenizer, applying Precellys ceramic beads (1.4 mm) (PEQLAB, Erlangen, Germany), and the RNeasy Mini Kit (Qiagen) according to the manufacturer’s specifications.
Sample and patient data
Reverse Transcription.
Five micrograms total RNA was reversely transcribed using the SuperScript II reverse transcriptase (Invitrogen, Carlsbad) and the poly-dT-Primer QT (Table 2) according to the protocol of Scotto-Lavino et al. (2006). Residual RNA fragments were digested using RNAse H (Invitrogen), and the generated cDNA was diluted with 75 µl Tris-EDTA buffer (Roth).
Primer and primer specifications
Primer Generation, 3′-RACE Reaction, and Purification of PCR Products.
All gene-specific forward primers were generated using the Primer-BLAST online tool with standard settings (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) based on the NCBI ABCC1/2/3 sequences (NM_004996.3, NM_000392.4, and NM_003786.3). A schematic representation of primers used in the present study is shown in Fig. 1, and primer specifications and sequences are shown in Table 2. 3′-RACE reactions were performed in a two-step procedure according to the protocol of Scotto-Lavino, using various primer combinations as shown in Fig. 1 (Scotto-Lavino et al., 2006). In a first amplification round, gene-specific forward primers with binding site upstream the stop codon in combination with the QO reverse primer, complementary to a rear artificial tail located at the 3′-end of the cDNA generated by the QT poly-dT primer, were used. In a second amplification round, diluted amplification products (1:20 in Tris-EDTA buffer) of the first round were used as a template in a nested PCR approach. For this, gene-specific forward primers binding downstream the forward primers used for the first amplification round were combined with R-QI primers, complementary to the anterior artificial end of the cDNA 3′-end. A schematic representation of all 3′-RACE reactions is shown in Fig. 1. All experiments were performed on a GeneAmp PCR system 9700 (Applied Biosystems, Darmstadt, Germany) using the HotStarTaqMaster Mix Kit (Qiagen). All products used for sequencing were purified by gel electrophoresis (1.5% w/v agarose, supplemented with 1 mM guanosine) and subsequent gel extraction using the QIAquick Gel Extraction Kit (Qiagen) according to the manufacturer’s recommendations.
Schematic representation of all forward primers used in the present study. The numbers in the primer names specify the first nucleotide of their binding site within ABCC1 (at the top), ABCC2 (in the middle), and ABCC3 (at the bottom). The 3′-UTR is shown in dark gray; the coding sequence upstream the 3′-UTR and the genetic sequence downstream the 3′-UTR is shown in light gray. Additional sequence position information is shown at the top of ABCC1, ABCC2, and ABCC3 (end of CDS, start and end of the respective 3′-UTR according to NCBI, and the accession numbers of mRNA and gene sequences). Primers named with “gen” bind in genomic regions behind the given 3′-UTR (NCBI database) to observe if longer unknown 3′-UTRs exists in the respective mRNA.
Sequencing and Sequence Analysis.
PCR products were verified by Sanger sequencing at the Institute of Clinical Molecular Biology of the University Hospital Kiel, Campus Kiel, Germany. Obtained sequences were analyzed and aligned using the BioEdit Sequence Alignment Editor version 7.2.6 (http://www.mbio.ncsu.edu/BioEdit/).
Quantitative Real-Time PCR.
To quantify different ABCC2 3′-UTR lengths present in the human tissue samples, quantitative real-time PCR experiments were performed, applying three different primer combinations (Table 2) with binding sites in distinct ABCC2 3′-UTR isoforms by using the SYBR Green SYBR SelectTM Master Mix (Applied Biosystems) on an ABI Prism 7900 HT Sequence Detection System (Applied Biosystems). The relative 3′-UTR amount was calculated by normalization with β-actin.
Luciferase Reporter-Gene Assays.
Primers to amplify the ABCC2 3′-UTR variants with 178/429/944 bp were generated using the TaKaRa In-Fusion Cloning primer designing tool (TaKaRa Bio Group; https://www.takarabio.com/learning-centers/cloning/in-fusion-cloning-tools), and PCR products were amplified with the forward primer Fhd (5′- GCTCGCTAGCCTCGAGCTGGCATTGAGAATGTGAA-3′) combined with the reverse primers Rhd178 (5′- CGACTCTAGACTCGACAATCGAGGGGTTTCTC-3′), Rhd429 (5′-CGACTCTAGACTCGATGCACCTATTTGCATCACCA-3′), and Rhd944 (5′- CGACTCTAGACTCGAAAAAATTCACAAGACATACAAGGAA-3′) using the HotStarTaq Master Mix Kit (Qiagen). The pmirGLO dual luciferase vector (Promega, Mannheim, Germany) was digested overnight with the restriction enzyme XhoI (NEB, Ipswich, MA). Amplicons were cloned into the pmirGLO vector using the In-Fusion HD Cloning Kit (Takara Bio) according to the manufacturer specifications. Correct insertion and orientation of the amplicons into the vector was confirmed by Sanger sequencing. A total of 105 HepG2 cells were cotransfected with 100 ng of the target vector or empty control vector and 10/25/50 nM of hsa-miR-379-5p Pre-miR miRNA Precursor (PM10316; Thermo Fisher Scientific, Germany) using the siPORT NeoFX transfection agent (Ambion). Pre-miR micro-RNA precursor negative control number 1 (AM17110; Thermo Fisher) served as negative control. Forty-eight hours after transfection, HepG2 cells were lysed, and firefly and renilla luciferase activity was measured using the Dual-Luciferase Reporter Assay System (Promega) and a Veritas microplate luminometer (Turner Biosystems, CA). Firefly luciferase activity was normalized to renilla luciferase activity (internal standard), and normalized firefly luciferase activity obtained with pre-miR-379 was compared with the activity obtained with the pre-miR negative control.
Immunoblotting.
Human tissue samples were cut into 15-µm slices using the CryoStar NX50 (Thermo Fisher), collected into Precellys Tubes with 2.8-mm ceramic beads, and homogenized at 6000 rpm for 20 seconds in the Precellys homogenizer (VWR) in 500 µl radioimmunoprecipitation assay lysis buffer buffer supplemented with 0.2% w/v BSA. Subsequent lysis was performed by sonication for 13 seconds. Enrichment of membrane proteins was accomplished using the Plasma Membrane Protein Extraction Kit (Abcam, Berlin, Germany) with 500 µl Homogenize Buffer Mix according to the manufacturer’s protocol. Protein quantification and Western blotting were executed according to standard protocols as previously described (Waetzig et al., 2019). Blots were probed with the following antibodies: ABCC2 (c-518048, dilution 1:100; Santa Cruz Biotechnology), GAPDH (sc-47724, dilution 1:5000; Santa Cruz Biotechnology), and anti-mouse (IRDye 800CW 926-32210, dilution 1:10,000; LiCOR, Bad Homburg, Germany). Primary antibodies were diluted in Odyssey Blocking Solution and secondary antibodies in TBS with 0.2% v/v Tween. Blots were visualized using the Odyssey CLx imager (LiCOR). Densitometry was performed using Empiria Studio Software 1.1 (LiCOR).
Results
3′-RACE and Agarose Gel Electrophoresis
Using HepG2, SK-Hep-1, or Caco-2 mRNA as templates, numerous amplicons of ABCC1, ABCC2, and ABCC3 mRNAs were generated using various forward primers combined with the reverse primers QO and/or QI (Fig. 1; Table 2). Some of the primer combinations generated side products with either no sequence homologies to ABCC1, ABCC2, or ABCC3 or to products that have been unevaluable because of sequence overlays (black arrows, Fig. 2), or they had the expected 3′-UTR sequence but no detectable poly-A tail. However, the majority of primer combinations, especially for ABCC2 and ABCC3, led to PCR products with evaluable 3′-UTR sequences, which are required for exact determination of 3′-UTR variants (white arrows, Fig. 2). Most of the products were generated by the second amplification round (Fig. 2A, right gel; Fig. 2B; Fig. 2C, middle and right gel), and only a few products were directly evaluable after sequencing of the first amplification round (Fig. 2A, left gel; Fig. 2C, left gel). All evaluable sequences with detectable poly-A tails at the 3′-UTR ends were analyzed. The obtained consensus sequences are shown in Supplemental Fig. 1. No detectable amplicons were generated with the forward primers ABCC1_Fgen6629, ABCC1_Fgen6674, ABCC1_Fgen6727, ABCC3_Fgen6105, ABCC3_Fgen6107, and ABCC3_Fgen6193 (Fig. 1; Table 2), as sufficiently long 3′-UTR sequences were missing. The ABCC1 forward primers ABCC1_F4884, ABCC1_F4943, and ABCC1_F5074 (Fig. 1; Table 2) generated products homologous to the ABCC1 3′-UTR but with no detectable poly-A tail, indicating no further 3′-UTR variants between the two observed ABCC1 3′-UTRs. Five different 3′-UTR variants for ABCC2 and four 3′-UTR variants for ABCC3 were identified.
3′-RACE and agarose gel electrophoresis. Agarose gels are shown that resulted or partly resulted in evaluable ABCC1/ABCC2/ABCC3 3′-UTR sequences with a visible poly-A tail. Arrows indicate gel extracted and sequenced amplicons. A white arrow indicates that sequencing of the respective amplicon resulted in evaluable ABCC1/ABCC2/ABCC3 3′-UTR sequences with a detectable poly-A tail. Black arrows indicate artifacts (side products) unevaluable because of sequence overlays, or they represent the expected sequence without a detectable poly-A tail. If amplicons were generated in the first amplification round, this is noted at the bottom end of the gel. All other gel runs represent results generated with the second amplification round (nested PCR). In these cases, the respective primer combination used in the first amplification round is shown at the top of the gels. Primer combinations used to generate the amplicons are shown below the gels. Agarose gels of (A) ABCC1, (B) ABCC2, and (C) ABCC3. Total RNA of the following cell cultures was used as template: ABCC1, Lane 1, 2, 5, 6, 9, 10, 13, 14, 17, 18, SK-Hep-1; Lane 3, 4, 7, 8, 11, 12, 15, 16: Caco-2; ABCC2 (Lane 1–16) and ABCC3 (Lane 1–20): HepG2. M, marker (100 bp DNA ladder; Roth and 1 kb DNA ladder; Roth); ampl., amplification round.
Determination of mRNA 3′-UTR Lengths in ABCC1, ABCC2, and ABCC3
ABCC1.
According to the NCBI reference sequence NM_004996.3, the ABCC1 gene consists of 6564 bp with a coding sequence of 4595 bp. The ABCC1 3′-UTR starts at position 4772 after the TGA stop codon and has a length of 1793 bp ending at position 6564 (poly-A site). However, our 3′-RACE experiments on cDNA from hepatoblastoma and colon carcinoma cell lines demonstrated that ABCC1 is expressed with two different mRNA 3′-UTRs, a short one of 105 bp (counted from the first base downstream the stop codon up to the last base upstream the poly-A tail) and one markedly prolonged 3′-UTR of 1771 bp (Fig. 3). A 3′-UTR of 1793 bp, as noted in the NCBI database, could not be confirmed in our study.
Representation of the observed 3′-UTR variants of ABCC1, ABCC2, and ABCC3. 3′-UTRs observed in the present study are shown in dark gray, and a small part of the coding sequence upstream the 3′-UTR is shown in light gray. The start of the 3′-UTR is marked by an interrupted vertical line, and the poly-A tail is indicated by four A’s. An asterisk at the indicated poly-A tails suggests the most common 3′-UTR variant. The exact length of the observed 3′-UTRs is given within or to the right of the bars representing the respective mRNAs.
ABCC2.
According to the latest version in the NCBI database, the ABCC2 gene spans 5806 bp with a coding sequence of 4637 bp. No polyadenylation signals or polyadenylation sites are noted in the database. The ABCC2 consensus sequence from our 3′-RACE experiments was identical to the sequence of the NCBI database entry. Nevertheless, the last 22 nucleotides of the NCBI sequence, denoted as part of exon 31, were not present in the 3′-UTRs identified in this study. Overall, a total of five different 3′-UTRs with lengths of 178, 363, 429, 840, and 944 bp, respectively, were discovered (Fig. 3).
ABCC3.
Twelve different putative mRNA 3′-UTRs for ABCC3 are listed in the NCBI AceView database, ranging from 46 to 2077 bp (spliced mRNA variants). It should be noted that two different gene variants exist for ABCC3. According to the NCBI database, variant 1 encodes the longer isoform (NM_003786.3), whereas variant 2 lacks multiple 3′ exons and has an alternative 3′ sequence, and the resulting isoform is much shorter with a different C-terminus (NM_001144070.1). However, as the fully functional ABCC3 transporter is based on variant 1 (isoform NM_003786.3), the present study refers only to this variant. In the NCBI nucleotide database, the ABCC3 entry (NM_003786.3) has a 3′-UTR length of 502 bp, with the most common AATAAA polyadenylation signal 24 bp upstream of the polyadenylation site. This 3′-UTR was confirmed in the present study. In addition, three alternative 3′-UTR variants were identified for ABCC3 with lengths of 306, 1053, and 1427 bp, respectively (Fig. 3).
Functional Verification
To prove the concept that miRNA-mediated posttranscriptional regulation is dependent on alternative polyadenylation and can be impaired by shortening of the respective 3′-UTR, we cloned three ABCC2 3′-UTR variants (178-, 429-, and 944-bp–long 3′-UTRs) in a luciferase vector and cotransfected the respective vector constructs along with miR-379 in HepG2 cells. As expected, the shortest variant (178 bp) was not regulated by miR-379, whereas the longer variants (429 and 944 bp) both showed a significant reduction of reporter-gene expression in the presence of miR-379 (Fig. 4). These results confirm the influence of alternative polyadenylation on the posttranscriptional regulation of ABCC2 through miR-379.
Proof of principle using luciferase reporter-gene assays applying three ABCC2 3′-UTRs and miR-379. The binding site of miR-379 within the ABCC2 3′-UTR is shown at the top. The upper diagram shows the results of the reporter-gene assay determined with miR-379 and the 178 bp ABCC2 3′-UTR variant (dark gray; no reduction of reporter-gene activity) compared with the 429 bp variant (light gray; concentration-dependent reduced reporter-gene activity) normalized to the negative control (black). The lower diagram shows the results determined with miR-379 and the 178 bp ABCC2 3′-UTR variant (dark gray; no reduction of reporter-gene activity) compared with the 944 bp variant (light gray; concentration-dependent reduced reporter-gene activity) normalized to the negative control (black). The pre-miR miRNA Precursor Negative Control 1 (Thermo Fisher Scientific) served as negative control. miR-379 was applied in three different concentrations (10, 25, and 50 nM) as indicated below the diagrams. Data were analyzed using ANOVA and Tukey Honestly Significant Difference post-hoc tests to detect statistically significant differences of 1) the negative control compared with different miR-379 concentrations and 2) matching miR concentrations of the long and short fragments. Comparisons are indicated with horizontal lines on top; asterisks indicate significance. Error bars show ± one S.D. *P < 0.05; **P < 0.01; ***P < 0.001.
Distribution of ABCC2 3′-UTR Variants in Human Tissue
The distribution of the three ABCC2 mRNA 3′-UTR variants (178, 363, 840 bp) were investigated in human peritumoral nonmalignant liver, gall bladder, and colon tissues by quantitative real-time PCR using SYBR Green and primer combinations with binding sites within the respective 3′-UTR, generating amplicons with approximately equal lengths (Fig. 5). Taking into account that primer pair 1 binds to all possible 3′-UTRs, primer combination 2 binds to the 363/429/840/944 bp variants, and primer combination 3 binds to the two longest 3′-UTRs (840 and 944 bp), the results show that the long 3′-UTR variants (840 and 944 bp ABCC2 3′-UTR isoforms) are predominantly expressed in the liver sample compared with colon (P = 0.03) or gall bladder (P = 0.04).
Distribution of ABCC2 3′-UTR variants in human tissue. The binding regions of used primers pairs for quantitative real-time PCR experiments are shown at the top. Three primer pairs were used with binding sites within the 178, 363, and 429 bp 3′-UTR variants of ABCC2. The relative expression to β-actin, normalized to the expression of the 178 bp 3′-UTR variant, is shown for human liver (white), colon (light gray), and gall bladder samples (dark gray) as box whisker plots. Error bars depict minimum and maximum values. Data were analyzed using ANOVA with subsequent Tukey post-hoc test; *P < 0.05.
ABCC2 Protein Quantification
Western blotting revealed that only the liver samples exhibit a consistent ABCC2 protein synthesis (Fig. 6). Four of the five analyzed gall bladder samples show a minor but detectable ABCC2 abundance, whereas there was lack of ABCC2 protein within the colon samples. Calculation of the mean intensity optical densitometry (IOD) of the bands confirmed significant ABCC2 synthesis in the liver samples, absence of ABCC2 protein in the colon samples, and marginal ABCC2 synthesis in the gall bladder .
ABCC2 protein quantification by immunoblotting and IOD calculation. Western blots of ABCC2 using gall bladder, liver, and colon whole cell lysates and membrane enrichment of one respective liver sample are shown at the top. GAPDH was used as positive control. The calculated IOD (mean IOD) of ABCC2 protein in the cell lysates normalized to GAPDH is shown at the bottom. In accordance with the antibodies datasheet (https://datasheets.scbt.com/sc-518048.pdf) and the literature, ABCC2 shows a size variance in immunoblots, probably because of posttranslational modifications. C, cytosolic fraction; M, membrane fraction.
Discussion
Alternative polyadenylation occurs in the majority of eukaryotic cells and is an important regulatory mechanism of gene expression (Tian and Manley, 2013). The 3′-UTR of mRNAs serves as a major regulatory region for trans-acting factors controlling stability, cellular localization, and translation efficiency of the target mRNAs (Elkon et al., 2013). The latter is especially true for miRNAs. The importance of alternative polyadenylation was shown for various biologic processes, including cell development, differentiation, and proliferation (Elkon et al., 2013; Carpenter et al., 2014; Jia et al., 2017). Widespread 3′-UTR shortening by alternative polyadenylation was reported in cancer cell development, leading to oncogene activation (Mayr and Bartel, 2009). Alternative polyadenylation also occurs in ABC transporters, as reported for ABCB1 (Bruhn et al., 2016) and ABCG2 (To et al., 2009), both highly relevant for drug bioavailability contributing to anticancer drug resistance (Fohner et al., 2017; Yakusheva and Titov, 2018). Besides ABCG2 and ABCB1, ABC transporters of the C-family, especially ABCC1, ABCC2, and ABCC3, are relevant for pharmacokinetics as well as internal nutrient distribution. Because of their biological importance, it is essential to understand the regulation of these transporters and the extent of alternative polyadenylation they underlie. Here, we identified two different 3′-UTR length variants in ABCC1 mRNA, five in ABCC2, and four in ABCC3.
A comparison of the resulting ABCC1 mRNA sequence with the NCBI database entry revealed no sequence varieties, except for the last 22 bp of the 3′-UTR, which are absent in the consensus sequence. The ABCC1 3′-UTR starts with a GCC-triplet downstream the TGA stop codon. A TATATC motif at the end of the 1771 bp variant and an ACCAAA motif at the end of the short 105 bp 3′-UTR was observed. Potential polyadenylation signals for the short ABCC1 3′-UTR might be GCCTCC or CCTCCC located 22 or 21 nucleotides upstream the poly-A tail (Gruber et al., 2016). However, these are rare polyadenylation signals in eukaryotes, which are not experimentally confirmed. For the long 3′-UTR, potential polyadenylation sites are TAAAAA, AAAAAT, or AAAATA (21, 17, and 16 nucleotides upstream the poly-A tail). Three different ABCC1 3′-UTR fragments were previously described in the experimental NCBI AceView Database (aAug10-gAug10) with lengths of 1793, 244, and 227 bp. Compared with the present analysis, these entries could not be confirmed despite the usage of 12 different forward primers for the 3′-RACE experiments.
The shortest ABCC2 3′-UTR is 178 bp long and has a common AATAAA polyadenylation signal 23 nucleotides upstream of the polyadenylation site (Tian et al., 2005; Gruber et al., 2016). The 363 bp 3′-UTR has putative polyadenylation sites 36 and 38 bp upstream (TTTTTT and CCTTTT) the polyadenylation site, and the 429 bp 3′-UTR has three potential polyadenylation signals (TTCTTT, TCTTTT, and TTTTGT) 47, 46, and 44 nucleotides upstream the polyadenylation site that are rare in eukaryotes (Gruber et al., 2016). The putative polyadenylation signals of the 840 and 944 bp ABCC2 3′-UTRs are more common. These are TTTTTT or TTTAAA for the 840 bp 3′-UTR (30 and 27 nucleotides upstream the poly-A site) and TTTTTA or TTTATT for the 944 bp fragment (28 and 26 nucleotides upstream the poly-A site). The AceView database reveals six entries for ABCC2 3′-UTR variants (13, 56, 292, 750, and 945 bp). Only the 945 bp fragment (944 bp in the present study) was confirmed. Neither the short ABCC2 3′-UTRs variants with a length of 13 or 56 bp nor the 292- or 750-bp–long variants were found.
Except for the 502 bp 3′-UTR, ABCC3 entries in AceView were not confirmed. The 306 bp ABCC3 3′-UTR variant found here has a common polyadenylation signal (AAATAA) 22 nucleotides upstream the poly-A tail. In contrast, the 1053 and 1427 bp variants have rare polyadenylation signals, potentially ATTAAA for the 1053 bp 3′-UTR (28 nucleotides upstream the polyadenylation site) and ATAAAG for the 1427 bp fragment (24 nucleotides upstream the polyadenylation site). The ABCC3 3′-UTR starts with an AAT triplet following the TAA stop codon.
Various miRNAs were identified to bind to the 3′-UTR of the three ABCC transporters (Table 3). It should be noted that for three of them (miR-9 for ABCC1, miR-205-5p for ABCC2, and miR-143 for ABCC3), the reported miRNA binding sites are not present in the current NCBI mRNA sequence entries. However, the majority of miRNAs have binding sites located in prolonged 3′-UTRs of ABCC1, ABCC2, and ABCC3 but not in the short 3′-UTR variantsFig 7. Thus, these miRNAs might be affected by alternative polyadenylation, and their binding sites disappear when shortening of the respective 3′-UTR occurs. As a consequence, the posttranscriptional downregulation of expression may be reduced, resulting in higher protein synthesis. Hence, likewise for ABCG2 (To et al., 2009) or ABCB1 (Bruhn et al., 2016), the variation of 3′-UTR lengths in members of the ABCC family may also contribute to the disparity of tissue expression, as exemplified here for ABCC2 by performing reporter-gene assays, revealing an impaired miRNA interaction and subsequent impact on gene expression dependent on the 3′-UTR length. Whereas the expression of the short ABCC2 178 bp 3′-UTR was not affected in the reporter gene assay, the longer 429 and 944 bp 3′-UTRs were downregulated when cotransfected with miR-379, an miRNA binding to the ABCC2 3′-UTR, as proven before by our group (Haenisch et al., 2011). In this earlier study, we verified the binding of miRNA-379 to the ABCC2 3′-UTR and an miRNA-379–dependent regulation of ABCC2 mRNA and protein. Here, we observed that miRNA-379–dependent regulation of ABCC2 3′UTR is lost by the presence of shortened 3′UTR lengths. This shows that the ABCC2 mRNA may lose posttranscriptional control through alternative polyadenylation by 3′-UTR shortening. Hence, a possibly higher ABCC2 protein content may contribute to altered drug response.
miRNAs that have been reported and confirmed by luciferase reporter-gene assays to bind to the 3′-UTRs of ABCC1, ABCC2, and ABCC3
Detailed sequences of ABCC1, ABCC2 and ABCC3 showing regulatory site tags and miRNA bindings sites.. The stop codons are written in bold letters and are underlined at the beginning of each 3′-UTR. The binding sites of reported miRNAs for ABCC1 (at the top), ABCC2 (in the middle), and ABCC3 (at the bottom) are highlighted in light gray, and the name of the respective miRNA is given below. Putative poly-A signals and corresponding poly-A sites are highlighted in dark gray. The length of the observed 3′-UTRs with ABCC1, ABCC2, and ABCC3 is given in brackets at every indicated poly-A site.
The importance of ABCC2 in pharmacokinetics underlines the need for a better understanding of its regulation because of its numerous substrates, i.e., cytostatics and antibiotics (Nies and Keppler, 2007; Bruhn and Cascorbi, 2014). According to the extent of possible ABCC2 3′-UTR variants shown in the present study, alternative polyadenylation must be taken into account as one of the major regulation mechanisms of ABCC2. Here, only one miRNA interference was analyzed, but it confirmed the hypothesis that the shorter the 3′-UTR is, the lower is the miRNA-dependent suppression of mRNA expression. Further mRNA/miRNA interactions of ABC transporters should be performed to analyze the extent of posttranscriptional regulation of these important transport proteins.
By performing quantitative real-time PCR experiments, we observed a higher amount of prolonged ABCC2 3′-UTRs in the liver compared with gall bladder and colon. These findings suggest a high likelihood of miRNA-mediated posttranscriptional ABCC2 control in the liver and, in comparison, a potentially reduced regulation in gall bladder and colon. It must be taken into account that all tissue samples were obtained from cancer patients undergoing tumor surgery with a possibly changed homeostatic state of the respective site. It was observed that ABCC2 and ABCG2 expression was altered in moderate dysplasia in colorectal carcinogenesis, suggesting an involvement of ABC transporters in early carcinogenesis (Andersen et al., 2015). Alternative polyadenylation may contribute to ABCC2 expression changes in colorectal carcinogenesis (Mayr and Bartel, 2009).
Immunoblotting revealed the highest ABCC2 protein abundance in the liver, minor ABCC2 protein in the gall bladder, and almost no detectable ABCC2 protein in the colon. The high amount in the liver is in line with the role of hepatic ABCC2 as an efflux pump to extrude xenobiotics or metabolites, i.e., bilirubin conjugates (König et al., 1999; Cascorbi, 2006). Nevertheless, ABCC2 is known to be low expressed in the colon compared with the small intestine and underlies a large interindividual variation (Berggren et al., 2007). In gall bladder, ABCC2 is detectable but upregulated in gall bladder carcinoma (Kim et al., 2013). Here, quantitative real-time PCR analysis of ABCC2 3′-UTR length polymorphisms in tissue samples revealed long ABCC2 3′-UTR variants in the liver. However, the observation of long 3′-UTRs in liver tissue does not imply a downregulation in vivo but potentially mirrors an important mechanism for miRNA-mediated fine-tuning of posttranscriptional gene expression. This suggests that the tissue abundance of ABCC transporters does not necessarily reflect 3′-UTR lengths and/or miRNA interactions but is mainly defined by other regulatory mechanisms, such as transcription factor expression.
Concluding Remarks
The aim of our study was to identify 3′-UTR lengths of three members of the ABCC transporter family and its functional implications. First, we identified 3′-UTR length polymorphisms for ABCC1, ABCC2, and ABCC3. In addition, the functional significance of 3′-UTR length variants on miRNA binding were shown for ABCC2.
Regarding cell type, cell line, or cellular status–dependent alterations in alternative polyadenylation, this might explain conflicting results in miRNA studies. Therefore, alternative polyadenylation must be taken into account when analyzing miRNAs or miRNA/target gene interactions. This is also true for widely used ready-to-use plasmids applied in reporter-gene assays, representing only limited possibilities in observing posttranscriptional regulation.
The multitude of different 3′-UTR variants found for ABCB1, ABCG2, and members of the ABCC family suggest that besides transcriptional control, alternative polyadenylation is a regulatory mechanism of protein synthesis. Alternative polyadenylation should also be taken into account in pharmacogenetic studies of ABC transporters, particularly in malignant tissue exhibiting wide variation in miRNA expression profiles.
Acknowledgments
We thank Britta Schwarten for her excellent technical assistance.
Authorship Contributions
Participated in research design: Bruhn, Lindsay, Cascorbi.
Conducted experiments: Bruhn, Lindsay, Wiebel, Kaehler, Nagel.
Contributed new reagents or analytic tools: Bruhn, Lindsay, Kaehler, Nagel, Röder.
Performed data analysis: Bruhn, Lindsay, Kaehler, Nagel, Böhm.
Wrote or contributed to the writing of the manuscript: Bruhn, Lindsay, Kaehler, Nagel, Cascorbi.
Footnotes
- Received March 2, 2019.
- Accepted November 12, 2019.
↵1 O.B and M.L. contributed equally to this work.
This work was supported by the Werner und Klara Kreitz-Stiftung Kiel.
↵
This article has supplemental material available at molpharm.aspetjournals.org.
Abbreviations
- ABC
- ATP-binding cassette
- HepG2
- human liver carcinoma cell line
- IOD
- intensity optical densitometry
- miRNA
- micro-RNA
- MRP
- multidrug resistance-associated protein
- P/S
- penicillin/streptomycin
- RACE
- rapid amplification of CDNA-ends
- UTR
- untranslated regions
- Copyright © 2020 by The American Society for Pharmacology and Experimental Therapeutics