Identification of miR-193b Targets in Breast Cancer Cells and Systems Biological Analysis of Their Functional Impact*

Identification of protein targets for microRNAs (miRNAs) is a significant challenge due to the complexity of miRNA-mediated regulation. We have previously demonstrated that miR-193b targets estrogen receptor-α (ERα) and inhibits estrogen-induced growth of breast cancer cells. Here, we applied a high-throughput strategy using quantitative iTRAQ (isobaric tag for relative and absolute quantitation) reagents to identify other target proteins regulated by miR-193b in breast cancer cells. iTRAQ analysis of pre-miR-193b transfected MCF-7 cells resulted in identification of 743 unique proteins, of which 39 were down-regulated and 44 up-regulated as compared with negative control transfected cells. Computationally predicted targets of miR-193b were highly enriched (sevenfold) among the proteins whose level of expression decreased after miR-193b transfection. Only a minority of these (13%) showed similar effect at the mRNA level illustrating the importance of post-transcriptional regulation. The most significantly repressed proteins were selected for validation experiments. These data confirmed 14–3-3ζ (YWHAZ), serine hydroxyl transferase (SHMT2), and aldo-keto reductase family 1, member C2 (AKR1C2) as direct, previously uncharacterized, targets of miR-193b. Functional RNAi assays demonstrated that specific combinations of knockdowns of these target genes by siRNAs inhibited growth of MCF-7 cells, mimicking the effects of the miR-193b overexpression. Interestingly, the data imply that besides targeting ERα, the miR-193b effects include suppression of the local production of estrogens and other steroid hormones mediated by the AKR1C2 gene, thus provoking two separate molecular mechanisms inhibiting steroid-dependent growth of breast cancer cells. In conclusion, we present here a proteomic screen to identify targets of miR-193b, and a systems biological approach to mimic its effects at the level of cellular phenotypes. This led to the identification of multiple genes whose combinatorial knock-down likely mediates the strong anti-cancer effects observed for miR-193b in breast cancer cells.


Identification of protein targets for microRNAs (miRNAs)
is a significant challenge due to the complexity of miRNAmediated regulation. We have previously demonstrated that miR-193b targets estrogen receptor-␣ (ER␣) and inhibits estrogen-induced growth of breast cancer cells. Here, we applied a high-throughput strategy using quantitative iTRAQ (isobaric tag for relative and absolute quantitation) reagents to identify other target proteins regulated by miR-193b in breast cancer cells.
iTRAQ analysis of pre-miR-193b transfected MCF-7 cells resulted in identification of 743 unique proteins, of which 39 were down-regulated and 44 up-regulated as compared with negative control transfected cells. Computationally predicted targets of miR-193b were highly enriched (sevenfold) among the proteins whose level of expression decreased after miR-193b transfection. Only a minority of these (13%) showed similar effect at the mRNA level illustrating the importance of post-transcriptional regulation. The most significantly repressed proteins were selected for validation experiments. These data confirmed 14 -3-3 (YWHAZ), serine hydroxyl transferase (SHMT2), and aldo-keto reductase family 1, member C2 (AKR1C2) as direct, previously uncharacterized, targets of miR-193b. Functional RNAi assays demonstrated that specific combinations of knockdowns of these target genes by siRNAs inhibited growth of MCF-7 cells, mimicking the effects of the miR-193b overexpression. Interestingly, the data imply that besides targeting ER␣, the miR-193b effects include suppression of the local production of estrogens and other steroid hormones mediated by the AKR1C2 gene, thus provoking two separate molecular mechanisms inhibiting steroid-dependent growth of breast cancer cells.
In conclusion, we present here a proteomic screen to identify targets of miR-193b, and a systems biological approach to mimic its effects at the level of cellular phe-notypes. This led to the identification of multiple genes whose combinatorial knock-down likely mediates the strong anti-cancer effects observed for miR-193b in breast cancer cells. MicroRNAs (miRNAs) 1 regulate gene expression post-transcriptionally by binding primarily to the 3Јuntranslated region (3ЈUTR) of their target mRNAs, resulting in mRNA destabilization or translational repression (1,2). Genes encoding 1048 human miRNAs have so far been identified (miRBase v. 16.0) (3), and miRNAs are predicted to regulate the expression of up to 60% of all human protein-encoding genes (4). A large number of bioinformatic methods have been developed for miRNA target prediction (1), but these are still highly unspecific and inaccurate because miRNAs typically contain an imperfect match to their target sequences. In addition, a single miRNA can target hundreds of proteins and a single protein can be influenced by multiple miRNAs (1,5,6). Hence, comprehensive understanding of the phenotypic effects of miRNAs at the level of the entire cell is currently difficult. mRNA profiling by microarrays has been widely used for miRNA target identification, but microarrays only detect the effects of miRNAs at the transcriptional level, and will miss targets repressed solely at the translational level. However, compared with mRNA profiling techniques changes in protein levels are difficult to assess in a high-throughput fashion. Recently, we introduced a protein lysate microarray technique for high-throughput identification of miRNAs targeting estrogen receptor-␣ (ER␣) protein (7). In the present study, we have carried out a systematic analysis of target protein and mRNA levels in living cells after transfection with a pre-miRNA con- 1 The abbreviations used are: miRNA, microRNA; AKR, aldo-keto reductase; ER, estrogen receptor; LC-ESI-MS/MS, liquid chromatography electrospray ionization tandem mass spectrometry; HSP90AB1, heat-schock protein 90 ␤; iTRAQ, isobaric tag for relative and absolute quantitation; MCM7, DNA replication licensing factor MCM7; MSigDB, molecular signatures database; MHTFD1, C-1-tetrahydrofolate synthase; PTPLB, Protein-tyrosine phosphatase-like member B; SHMT2, serine hydroxyl transferase; UTR, untranslated region; YWHAZ, 14 -3-3 protein .
struct leading to forced up-regulation of the miRNA molecule. We compared quantitative proteomics data obtained by using iTRAQ (isobaric tag for relative quantitation) reagents with mRNA profiling data obtained from microarrays to characterize miRNA targets in breast cancer cells. We focused on dissecting the effects of miR-193b, which is known to directly target ER␣, and to inhibit growth of breast cancer cells by inducing a cell cycle arrest in G 1 phase (7). Using iTRAQ based quantitative mass spectrometry we found that the expression of 39 out of 743 investigated proteins was repressed by miR-193b, and that the predicted miR-193b targets were highly enriched among the down-regulated proteins. The impact of miR-193b on the protein levels of key targets were validated, and individual siRNA mediated knock-downs were performed to demonstrate that specific combinations of these could mimick the phenotypic effects observed after miR-193b transfection.

EXPERIMENTAL PROCEDURES
Cell Culture and Transfections-MCF-7 cells were obtained from Interlab Cell Line Collection (ICLC, Italy), and cultured in Dulbecco's modified Eagles medium (DMEM) (1 g/l glucose) supplemented with 10% fetal bovine serum (FBS), 2 mM L-glutamine and 1% penicillin/ streptomycin. Cells were regularly screened for Mycoplasma with PCR using Mycoplasma Plus™ PCR Primer Set from Agilent Technologies Inc. (Santa Clara, CA).
Human Pre-miR™ miRNA Precursors, Anti-miR™ inhibitors, and an siRNA for human estrogen receptor-␣ (GGCCAAAUUCA-GAUAAUCGTT) were purchased from Ambion (Austin, TX). miRCURY LNA™ inhibitor for miR-193b was obtained from Exiqon (Vedbaek iTRAQ Labeling-MCF-7 cells were transfected with 20 nM pre-miR negative control and pre-miR-193b (Ambion) in six-well plates. Two biological replicate experiments were performed. After 48 h incubation, the cells were harvested, and lysed in 150 l of 50 mM Tris-HCl pH 8.0, 100 mM NaCl and 0.5% Triton X-100. Complete mini EDTAfree protease inhibitor mixture (Promega Corp, Madison, WI) was added in the lysis buffer before use. Unsolubilized material was removed by centrifugation at ϩ4°C at 15,000 ϫ g for 20 min. Protein concentration was measured by DC Protein Assay kit (Bio-Rad Laboratories, Hercules, CA).
For iTRAQ labeling, 100 g of proteins were precipitated with six volumes of cold acetone at Ϫ20°C for 2 h. Protein digestion and labeling were done according to iTRAQ 4-plex Reagent Kit (Applied Biosystems, Foster City, CA). Protein pellets were suspended in dissolution buffer (0.5 M triethyammonium bicarbonate), and denatured and disulfide bonds were cleaved with 50 mM Tris-(2-carboxyethyl)phosphine. Each sample was then digested by 10 g of sequence grade modified trypsin (Promega) at ϩ37°C overnight, and labeled with the iTRAQ reagents.
Before pI-based peptide prefractionation, the samples were desalted by C18 Empore disks (3 M) and dried peptides were dissolved in IPG strip rehydration solution (4 M urea, 2% (v/v) IPG buffer pH 3-10). Peptides were loaded on a 13-cm IPG strip with a pH gradient of 3-10. Electrofocusing of peptides was performed with IPGphor (Amersham Biosciences, GE Healthcare, Piscataway, NJ) until 20 kVh. After focusing, the strip was cut into 12 fractions, peptides were extracted as described previously (8) and desalted again to remove urea.
Protein identification and quantitation were carried out using the ProteinPilot software (version 2.0.1, Applied Biosystems, MDS-Sciex) with the Paragon TM search algorithm (9). The software compares relative intensity of proteins present in samples based on the intensity of reporter ions released from each labeled peptide and automatically calculates protein ratios and p values for each protein. The data from 12 liquid chromatography tandem MS (LC-MS/MS) analyses were merged and searched against combined human Swiss-Prot (version 54.6) -TrEMBL (version 37.6) protein sequence database containing 76,122 protein entries. The following search parameters were used: iTRAQ 4-plex as the sample type, digestion with trypsin, and cystein alkylation with methyl methane thiosulfate. The precursor and fragment ion mass accuracy, number of missed cleavages permitted, fixed and variable modifications as well the peak list generating parameters are built-in functions of ProteinPilot. For protein identification, 95% confidence was used. Only those proteins with p Ͻ0.05 were accepted for quantitation and further validation experiments.
Real-time Quantitative PCR Analysis-Total cellular RNAs were isolated with miRVana™ Total RNA isolation kit (Ambion Inc.). For cDNA synthesis, 200 ng of total RNA was reverse transcribed with High Capacity cDNA Reverse Transcription kit (Applied Biosystems, Foster City, CA). Thereafter, the cDNAs were diluted 1/10 and the Taqman quantitative real-time PCR (qRT-PCR) analysis performed with Applied Biosystems 7900HT instrument using specific primers and probes for YWHAZ, SHMT2, AKR1C2, PTPLB, MCM7, HSP90AB1, and glyseraldehyde 3-phosphate-dehydrogenase (GAPDH) designed by the Universal Probe Library Assay Design Center (Roche Applied Biosciences, Basel, Switzerland) (supplementary Table S1). The results were analyzed with SDS 2.3 and RQ manager softwares (Applied Biosystems), and the expression of the mRNAs determined by relative quantitation method using GAPDH as an endogenous control. The data were collected from three separate biological experiments, which were each run with triplicates.
UTR Reporter Constructs and Luciferase Assays-Fragments covering the wild type and mutated miR-193b target sites in the YHWAZ, SHMT2, PTPLB, and MTHFD1 3Јuntranslated region (UTR) regions, and AKR1C2 5ЈUTR region (supplementary Table S2) were cloned into the Spe/MluI sites of a pMIR-REPORT™ Luciferase vector (Am-bion Inc.) downstream of a luciferase gene. For luciferase assays, MCF-7 cells (15,000 per well) were plated 24 h before transfection onto white, clear-bottom 96-well plates in normal culture medium without antibiotics. The cells were cotransfected with 50 ng of the 3ЈUTR reported plasmids, 50 ng Renilla luciferase plasmid, and with 50 nM pre-miR construct with Lipofectamine 2000 (Invitrogen Inc.) as previously described (7). The luciferase activity was assayed 24 h after transfection with Dual-Glo Luciferase Assay System (Promega Corp, Madison, WI) and measured with Envision Plate-reader (Perkin Elmer Inc, Wellesley, MA).
Proliferation Assays-For assaying the cell growth, MCF-7 cells (8000/well) were reverse transfected with 5 nM, 10 nM, and 20 nM miRNAs or siRNAs in 96-well plates using SiLentFect (Bio-Rad). The growth curves were generated during 88 h period of time by using the INCUCYTE TM Live-Cell Imaging System (Essen BioSciences, Ann Arbor, MI), which provides a time-lapse method for quantifying cell growth inside the cell culture incubator. Cell viability was measured with CellTiter-Glo Cell Viability Assay (Promega), as previously described (7).
microRNA Target Predictions-Predicted targets for miR-193b were obtained from TargetScan (10), PicTar (11), and miRBase (3). The lists of predicted targets from each prediction program were compared with the lists of miR-193b regulated genes separately. Enrichment of the predictions was calculated by comparing the proportion of predicted targets among the miR-193b down-regulated proteins to predicted targets among all proteins detected in the control sample, or by comparing the proportion of predicted targets among the miR-193b down-regulated transcripts to predicted targets among all genes on the array.

Identification of miR-193b Targets by iTRAQ Proteomics-
We have previously shown, that miR-193b directly targets ER␣, induces cell cycle arrest, and inhibits growth when overexpressed in breast cancer cells (7) (supplementary Fig. S1). To study the effects of miR-193b in MCF-7 cells in more detail and to identify additional protein targets for miR-193b, we used a proteomics approach. For that purpose, MCF-7 cells were transfected with a pre-miR construct for miR-193b and with a negative control scrambled pre-miR, and incubated for 48 h. Thereafter, the cells were harvested, lysed, and subjected to iTRAQ labeling. Relative protein expression differences between negative control pre-miR and pre-miR-193b-transfected MCF-7 cells were then measured by mass spectrometry (Fig. 1A). Protein identification and quantitation were done using the ProteinPilot software (version 2.0.1) with the Paragon TM search algorithm (9). The software compares relative intensity of proteins present in samples based on the intensity of reporter ions released from each labeled peptide and automatically calculates protein ratios and p values for each protein. For protein identification, 95% confidence was used. Only those proteins with p Ͻ0.05 were accepted for quantitation and further validation experiments.
With iTRAQ, a total of ϳ1000 proteins were identified, and 743 of those were unique, common proteins in two biological replicate analyses. Of those, 39 were down-regulated and 44 up-regulated 48 h post miR-193b transfection (Fig. 1B and supplementary Tables S3 and S4). We next compared the protein data with gene expression microarray data obtained from our previous study, where miR-193b targets were analyzed at the mRNA level 24 h post-transfection (7). For the microarray experiment, 24 h time point was selected in order to detect the more rapidly occurring changes at the mRNA level and to avoid indirect changes occurring at later time points. Most of the identified proteins (96%) were detected also at the mRNA level (Fig. 1B). However, only five (13%) of the candidate protein targets for miR-193b (SHMT2, YWHAZ, HN1L, PIP4K2C, and PTPLB) were significantly (Ͼtwofold) inhibited at the mRNA level according to the microarray profiling data (supplementary Table S3) (7). Overall, the correlation of the protein levels with the mRNA levels was moderate (r ϭ 0.30) (supplementary Fig. S2).
We showed previously, that ER␣-related gene signatures are highly enriched among the transcripts down-regulated by miR-193b (7). Here, we did similar gene set enrichment anal- yses for the proteomics data using the Molecular Signatures Database (MSigDB), and the results indicated that, in addition to the ER␣ signaling, cell cycle related pathways and genes are significantly enriched among miR-193b target proteins (supplementary Table S5).
Next, we tested how miRNA target predictions correlate with our data. miR-193b target predictions were obtained with three different algorithms: TargetScan (10), PicTar (11), and miRBase (3). As shown in Fig. 2A, the predicted targets for miR-193b were highly enriched (sevenfold) among the down-regulated proteins, as compared with all the proteins identified in a control MCF-7 sample. Similarly, in the microarray data, the predicted targets were highly enriched (5.5-fold) among the transcripts down-regulated by miR-193b, when compared with all genes on the array (Fig. 2B).
Validation of Putative miR-193b Targets by Western Blot and qRT-PCR Analyses-For validation experiments, we selected proteins that were identified in two replicate iTRAQ assays with a 95% confidence for protein identification (p Ͻ0.05 in both replicates), as calculated by the ProteinPilot software (Table I). Four of the most significantly down-regulated proteins (YWHAZ, SHMT2, PTPLB, and MTHFD1) were predicted to be targeted by miR-193b by TargetScan or miRBase, and three of those (YWHAZ, SHMT2, and PTPLB) were also among the most prominently down-regulated genes at the mRNA level (Table I). Western blot analyses of the selected proteins corroborated the potent down-regulation of AKR1C2, YWHAZ, PTPLB, SHMT2, MCM7, and HSP90AB1 (HSP90␤) proteins by miR-193b 48 h after transfection (Fig.  3A). qRT-PCR analyses showed that miR-193b down-regulated AKR1C2, YWHAZ, PTPLB, SHMT2, MCM7, and HSP90AB1 also at the mRNA level (Fig. 3B). These data were mainly in line with the microarray data, except for HSP90AB1, whose transcript was not down-regulated in the microarray analysis (Table I).
YWHAZ, SHMT2, and AKR1C2 are Direct Targets of miR-193b-YWHAZ, SHMT2, PTPLB, and MTHFD1 contain predicted miR-193b target sites at their 3ЈUTR region, whereas AKR1C2 contains a putative miR-193b site in its 5ЈUTR region (Fig. 4A). To analyze, whether miR-193b directly interacts with these predicted sites, we subcloned fragments covering the miR-193b binding sites into a luciferase reporter vector. In addition, we designed constructs with mutated miR-193b target sites. The resulting 3ЈUTR-luciferase reporter plasmids were transfected in MCF-7 cells together with pre-miR-193b and with a Renilla luciferase control plasmid. Transfection of miR-193b significantly inhibited the reporter activity of SHMT2 and YWHAZ 3ЈUTR and AKR1C2 5ЈUTR wild-type constructs, as compared with pre-miR negative controls (Fig.  4B). In contrast, miR-193b had no effect on the luciferase activity of the mutated constructs, indicating that miR-193b directly targets YWHAZ and SHMT2 by interacting with their 3ЈUTR, and AKR1C2 by interacting its 5ЈUTR. The luciferase activities of the PTPLB or MTHFD1 3ЈUTR constructs were not affected by miR-193b, indicating that miR-193b does not target its predicted sites in the PTPLB or MTHFD1 3ЈUTR.
To further study the importance of miR-193b in regulating SHMT2, AKR1C2, and YWHAZ, we inhibited miR-193b with specific anti-miRNAs and analyzed the effect on the protein level 96 h after transfection. As shown in Fig. 4C, miR-193b inhibition by 5 nM anti-miR-193b or LNA-anti-miR-193b resulted in up-regulation of SHMT2 and AKR1C2 protein levels (by 1.6-to 1.8-fold). With YWHAZ, the change was modest (1.4-fold), but consistently observed. Conceivably, there are also other factors regulating YWHAZ and inhibition of miR-193b is not enough to up-regulate YWHAZ.
siRNA-mediated Knockdown of miR-193b Regulated Genes Inhibits Breast Cancer Cell Growth-We have previously shown that miR-193b potently suppresses the growth and viability of MCF-7 cells in vitro ((7) and supplementary Fig. S1). To study   FIG. 2. Enrichment of predicted targets among the miR-193b down-regulated proteins and transcripts. (A and B) Target predictions were obtained with TargetScan, PicTar, and miRBase. Enrichment of the predictions was calculated by comparing the proportion of predicted targets among the miR-193b down-regulated proteins to predicted targets among all proteins detected in the control sample (A), or by comparing the proportion of predicted targets among the miR-193b down-regulated transcrips to predicted targets among all genes on the array (B). the functional significance of the newly-identified miR-193b targets, we tested the effect of siRNA-mediated knockdown of these genes on breast cancer cell growth. For that purpose, MCF-7 cells were transfected with siRNAs for YWHAZ, SHMT2, PTPLB, MCM7, AKR1C2, and HSP90AB1. In addition, we tested an siRNA for AKR1C1, which belongs to the same family as AKR1C2, and which was significantly down-regulated by miR-193b in one of the two replicate iTRAQ assays (data not shown). The siRNAs were transfected in MCF-7 cells alone or in various combinations, and the cell growth was monitored over 88 h period of time in Incucyte Live-Cell Imaging System, which provides a time-lapse method for quantifying cell growth inside the cell culture incubator. As expected, miR-193b significantly inhibited the growth of MCF-7 cells (Fig. 5A). In addition, knocking down of AKR1C1 or AKR1C2 potently inhibited cell growth, and combination of these two siRNAs further enhanced this effect, giving an additive effect ( Fig. 5A and  supplementary Fig. S3A). siRNA-mediated silencing of PT-PLB, YWHAZ, SHMT2, HSP90AB1, or MCM7 alone had no significant effect on MCF-7 cell growth (data not shown). However, simultaneous silencing of two of the miR-193b targets (PTPLBϩSHMT2, YWHAZϩMCM7, YWHAZϩSHMT2, or SHMT2ϩMCM7) suppressed cell growth, suggesting that these genes synergistically mediate miR-193b effects on MCF-7 cell growth ( Fig. 5A and supplementary Fig. S3A). In addition, most of the siRNAs were effective already at 5 nM concentration, as determined by siRNA dose response assays, which further highlights the importance of miR-193b and its targets on MCF-7 cell growth (supplementary Fig. S3A). The knock-downs with the siRNAs were verified by Western blotting (Fig. 5B) and by qRT-PCR (supplementary Fig. S3B).

DISCUSSION
Identification of miRNA targets is central to understanding the biologic roles of miRNAs. Computational and experimental approaches to identify miRNA targets have proven remarkably difficult because of the complexity of the miRNA-mediated regulation. Proteomic profiling of human cells with SILAC following miRNA transfections revealed how individual miRNAs modulated the levels of hundreds of proteins (5,6). Therefore, high-throughput methods for miRNA target identification are required.
In the present study, we have performed quantitative iTRAQ experiments in combination with gene expression microarray profiling, bioinformatic predictions, and functional assays to identify protein targets for miR-193b in breast cancer cells. miR-193b directly targets ER␣, inhibits breast cancer cell proliferation and induces cell cycle arrest in G 1 (7). With the iTRAQ analysis, we identified 39 proteins to be down-regulated by miR-193b at 48 h time point. Of these, only 13% were significantly repressed at the mRNA level according to microarray gene expression data obtained 24 h after miR-193b transfection, emphasizing the importance of protein (2010) demonstrated that miRNAs predominantly act by decreasing mRNA levels as a consequence of mRNA destabilization. The discrepancies here could partly be due to different time-points used in the microarray and proteomics analyses, 24 h and 48 h, respectively. Furthermore, the steady-state-levels of mRNA and protein rarely correlate directly because of complicated post-transcriptional control mechanisms (12,13). Although proteomic strategies are currently limited to the identification of only a fraction of the cell's proteome, they complement the microarray experiments as they offer a read-out of the functional component of the transcript activity. The exact mechanism of miRNA mediated repression is still unclear, and a combination of both mRNA and corresponding protein data are expected to reveal deeper insights into this mode of biological regulation. Furthermore, all methods for identification of miRNA target genes are susceptible to detection of indirect targets that are repressed through secondary effects, and experimental validation of the targets is required. The protein concentration changes measured by the iTRAQ method were not as high as in the Western blot validation data. This is interesting in itself, and it is not entirely clear as to why this occurs. Information on the accuracy of iTRAQ based methods has recently emerged, which enables clarification to some extent. Two reports have shown that the iTRAQ ratios tend to be underestimated and typically the iTRAQ experiments report fold changes of less than two or-ders of magnitude (14,15). However, although underestimation of protein ratios have been demonstrated, the direction of change (differential up-, down-, or unregulated) is itself unaffected. Therefore, the iTRAQ technology could be considered suitable as a primary screening method for miRNA target identification. However, we cannot exclude the possibility that there would be additional targets for miR-193b. The number of detected proteins in two experiments was 743, which does not cover the whole proteome. The MS-based proteomics is more biased toward more abundant proteins and may thus miss low-copy-number proteins. For example, Cyclin D1, a previously published miR-193b target (16), was not among the identified proteins in our study.
We verified YWHAZ, SHMT2, and AKR1C2 as novel, direct targets of miR-193b. miRNA regulation primarily occurs via interaction of miRNAs with the 3ЈUTR of their target genes (1). Interestingly, we demonstrate here that miR-193b regulates AKR1C2 by binding to its 5ЈUTR region. This is in accordance with recent reports suggesting that miRNAs can also target other regions of the transcripts besides the 3ЈUTR (17)(18)(19)(20). The prediction algorithms mainly use 3ЈUTR databases as the source of predictions. Therefore, they will miss putative miRNA targeting beyond the 3ЈUTR regions. For instance, AKR1C2 was not predicted to be targeted by miR-193b by TargetScan, PicTar, or miRBase.
Li et al. (21) have shown that during breast cancer cell metastasis, the expression of miR-193b is down-regulated, which in turn up-regulate urokinase-type plasminogen activa- Therafter, cell lysates were analyzed with Western blotting. Equal loading was confirmed by probing the same filters with ␤-actin. B, MCF-7 breast cancer cells were transfected with miR-193b or negative control miRNA and incubated for 24 h or 48 h. Thereafter, the total cellular RNAs were isolated and analyzed with Taqman qRT-PCR for the expression of AKR1C2, MCM7, PTPLB, SHMT2, YWHAZ, and HSP90AB1 mRNAs. The results normalized for GAPDH are shown as relative expression using negative control miRNA as a reference. The data were collected from three independent experiments, which were run with triplicates (Mean Ϯ S.D.). FIG. 4. YWHAZ, SHMT2, and AKR1C2 are direct targets for miR-193b. A, The complementarity of miR-193b sequence to the five down-regulated genes. B, Activity of luciferase reporters containing wild type (wt) or mutated (mut) putative miR-193b target sites in the YWHAZ, SHMT2, PTPLB, and MTHFD1 3ЈUTR, or AKR1C2 5ЈUTR. The reporters (50 ng) were co-transfected with Renilla luciferace control vector (50 ng) and with pre-miR-193b (50 nM), and the luciferase activities measured after 24 h incubation. Firefly luciferase activity was normalized to that of Renilla luciferase. The data shown are mean Ϯ S.E. of at least four independent experiments, each performed with five replicates. ** p Ͻ 0.01. C, The effect of miR-193b inhibition on the expression of SHMT2, AKR1C2 and YWHAZ. MCF-7 cells were transfected with negative control anti-miRNA (anti-Scr) or with miR-193b inhibitors (anti-193b or LNA-193b) at 5 nM concentrations. After 96 h incubation, cells were harvested and the cell lysates were analyzed with Western blotting. The numbers below each panel show quantitations relative to anti-Scr, which was set to 1.0. tor expression and contribute to the development of breast cancer. In addition, miR-193b expression was related to tumor stage and tumor grade, and correlated with the metastasis-free survival of breast cancer patients. Here, the identified miR-193b targets included other known breast cancerassociated genes, such as AKR1C2 and YWHAZ. YWHAZ belongs to the 14 -3-3 family of proteins, and it has been shown to be overexpressed in Ͼ40% of advanced stage breast cancers, confer cancer cell apoptosis resistance, and correlate with poor survival in breast cancer patients (22,23). AKR1C2 belongs to the aldo-keto reductase (AKR) superfamily of enzymes which participate in the local production of steroid hormones, such as estrogens, progesterones, and androgens (24). AKRs have been suggested as potential therapeutic targets in breast cancer due to their positive effects on proliferative signaling, which stimulates the growth of hormone-dependent and -independent breast cancer. Silencing AKR1C2 and another member of the same family, AKR1C1, very efficiently inhibited MCF-7 cell growth, suggesting that these genes are essential for MCF-7 cell survival. Recently, we demonstrated that ER␣ signaling is directly regulated by miR-193b in breast cancer cells (7). Together, these findings suggest that miR-193b regulates estrogen signaling by multiple mechanisms: by directly inhibiting the estrogen receptor as well as regulating the production of estrogen. This results in suppression of cell growth.
Our functional assays demonstrated that knocking down miR-193b targets by siRNAs phenocopied the anti-proliferative effects of miR-193b. However, with siRNA-mediated silencing potent growth inhibitory effects were mostly evident only after synergistic inhibition of multiple miR-193b target genes. This highlights the fact that miRNAs regulate the expression of several genes and their cell phenotypic effects are often not mediated by a single target.
To conclude, we present here that proteomic expression profiling is an efficient strategy to functionally identify protein FIG. 5. Inhibition of miR-193b targets suppresses MCF-7 cell growth. A, MCF-7 cells were transfected with miR-193b or with siRNAs for YWHAZ, SHMT2, PTPLB, MCM7, AKR1C2, AKR1C1, and HSP90AB1 at 20 nM concentration alone or in different combinations, as indicated. Cell growth was monitored over 88 h period of time in Incucyte Live-Cell Imaging System, which provides a time-lapse method for quantifying cell growth inside the cell culture incubator. As a negative control, Allstars negative control siRNA (Qiagen) was used, and as a positive control, an siRNA for KIF11 was used. The data represent average of two independent assays, which were both done with three technical replicates (mean Ϯ S.D.). B, The levels of SHMT2, YWHAZ, PTPLB, AKR1C2, and MCM7 proteins after siRNA transfections (72 h) were analyzed by Western blotting. targets and mechanisms of action for miRNAs that have anticancer effects when transfected into breast cancer cells. In addition, we provide functional data on novel targets for a breast cancer associated miRNA, miR-193b, whose effects likely involve several distinct steps of ER␣ signaling, such as steroid synthesis and down-regulation of the ER␣ receptor.