Transcriptome-wide analysis uncovers the targets of the RNA-binding protein MSI2 and effects of MSI2's RNA-binding activity on IL-6 signaling

The RNA-binding protein Musashi 2 (MSI2) has emerged as an important regulator in cancer initiation, progression, and drug resistance. Translocations and deregulation of the MSI2 gene are diagnostic of certain cancers, including chronic myeloid leukemia (CML) with translocation t(7;17), acute myeloid leukemia (AML) with translocation t(10;17), and some cases of B-precursor acute lymphoblastic leukemia (pB-ALL). To better understand the function of MSI2 in leukemia, the mRNA targets that are bound and regulated by MSI2 and their MSI2-binding motifs need to be identified. To this end, using photoactivatable ribonucleoside cross-linking and immunoprecipitation (PAR-CLIP) and the multiple EM for motif elicitation (MEME) analysis tool, here we identified MSI2's mRNA targets and the consensus RNA-recognition element (RRE) motif recognized by MSI2 (UUAG). Of note, MSI2 knockdown altered the expression of several genes with roles in eukaryotic initiation factor 2 (eIF2), hepatocyte growth factor (HGF), and epidermal growth factor (EGF) signaling pathways. We also show that MSI2 regulates classic interleukin-6 (IL-6) signaling by promoting the degradation of the mRNA of IL-6 signal transducer (IL6ST or GP130), which, in turn, affected the phosphorylation statuses of signal transducer and activator of transcription 3 (STAT3) and the mitogen-activated protein kinase ERK. In summary, we have identified multiple MSI2-regulated mRNAs and provided evidence that MSI2 controls IL6ST activity that control oncogenic signaling networks. Our findings may help inform strategies for unraveling the role of MSI2 in leukemia to pave the way for the development of targeted therapies.

expression of HOXA9 compared with healthy individuals (4,5). Although HOXA9 family of genes has been well-studied, less is known about the fusion partner MSI2. Musashi RNA-binding protein 2 (MSI2) is one of the two Musashi family of RNAbinding proteins (RBPs) containing two RNA recognition motifs (Fig. S1A) (6). MSI2 was found be highly expressed in multiple hematological malignancies including AML, CML, and B-ALL (2,7,8). MSI2 plays a major role as cell cycle regulator and maintainer of stem cell phenotype by balancing selfrenewal and differentiation in HSCs, thus balancing HSC homeostasis (9 -11). Despite some overall progress, the RNA-binding protein properties and the actually regulated mRNAs of MSI2 have not been studied in detail. Few publications have previously searched for the binding partners using techniques such as HITS-CLIP and RIP-CHIP, and although these techniques provide good insights into the interactome of the investigated RBP, they have several drawbacks related to the cross-linking efficiency and the use of UV-254 (12)(13)(14)(15). Therefore, we used PAR-CLIP (16) to gain insight into the transcriptome-wide mRNA targets and the RNA-binding motif of MSI2. PAR-CLIP utilizes a photoactivatable nucleoside, 4-S-uridine, to achieve higher cross-linking efficiency

MSI2 regulates IL-6 signaling
and determine binding sites with higher resolution compared with other techniques.
We here report targets regulated by MSI2 post-transcriptionally, which play roles in functions including DNA replication, DNA repair, and cell death and survival. Pathways regulated by MSI2 included EGF, HGF signaling, and eIF2 signaling pathways. The EGF and HGF pathways in turn regulate PI3K/ AKT, JAK/STAT, and ERK/MAPK pathways. We investigated one of the most highly regulated targets from our pathway analysis, IL6ST, a transmembrane signal transducer belonging to the IL-6 signaling cascade entity. IL-6 signaling includes multiple cytokines like IL-6, IL-11, LIF (leukemia inhibitory factor), OSM (oncostatin M) which together are called IL-6 type cytokines. These cytokines activate the receptor IL6ST which in turn phosphorylates downstream targets, thereby activating JAK/STAT, ERK, and PI3K/AKT signaling pathways (17,18). The list of targets and the pathways they regulate increase our understanding of the role of MSI2 in leukemia, and thus may pave the way for the development of targeted therapies.

Identification of transcriptome-wide direct mRNA targets of MSI2 by PAR-CLIP
We performed PAR-CLIP to identify the transcriptome-wide targets of MSI2 in human embryonic kidney cells (HEK293, Flp-In TREx) (Fig. 1B). The analysis of the MSI2 PAR-CLIP datasets revealed a total of 20,531 RBP-bound genomic regions after filtering for confident clusters (CCs). CCs contained two or more T-C conversion sites (cross-linked sites) and Ն10 overlapping reads (for cluster lengths and cluster site distances between CCs, see Fig. S1, B and C). Using gene annotation information for GRCh38 we identified 7378 distinct genes targeted by MSI2. We then examined the distribution of CCs across intronic and exonic regions of mRNA primary transcripts. MSI2 bound to intronic regions (39%) and 3ЈUTR regions (40%) equally, whereas only 7% bound to exonic regions (Fig. 1C).

Identification of the RNA-recognition element of the MSI2 protein
Next, we used the data from PAR-CLIP to identify the actual RNA sequence motif that was bound by the RBPs. Many RBPs bind RNA in a sequence-specific manner, with the recognized binding region usually referred to as the RNA-recognition element. This motif prediction was performed using MEME (multiple EM for motif elicitation). To elucidate the RRE of MSI2, we restricted our analysis to CCs with a T-C mutation frequency of Ն20% and two or more T-C positional changes. The genomic reference sequences of the 100 most abundantly sequenced clusters located in introns, exons, or 3ЈUTRs were used. MSI2 datasets revealed the motif as TTAG (E-value (MSI2) ϭ 9.3E-44) (Fig. 1D). We performed EMSAs to test a selection of synthetic RNAs representing the MSI2-binding sites using recombinant MSI2 protein. We tested a tetranucleotide repeat of the motif (UUAG) (length of 24 nucleotides) as well as a control with GC repeated 12 times to match the length of the motif probe. We observed a shift in the protein-RNA complex containing the predicted MSI2 motif, whereas a shift is completely absent using the control probe, indicating no binding between the RBP and the control probe (Fig. 1E). The exact motif was present in 76% of the top 100 MSI2 CCs, whereas the remaining clusters contained slight variations of the given motif.

Data integration of PAR-CLIP and RNA-Seq reveals direct mRNA target spectrum of MSI2
Next, we investigated which transcripts were preferentially bound and directly regulated by MSI2. Therefore, we performed siRNA-mediated knockdown (KD) of MSI2 in HEK293T cells ( Fig. 2A) and investigated the gene expression profiles with high-throughput RNA-sequencing (RNA-Seq) (Fig. 2B). We identified differentially expressed genes with stringent thresholds having a logFC (logarithm of fold change) Ͼ1.5 or ϽϪ1.5 and false discovery rate (FDR) Յ0.05 (Fig. 2C). Next, we integrated the PAR-CLIP data for MSI2 to identify physically bound and differentially regulated genes, which we term as the direct targets of MSI2. In total, we found 230 genes to be directly under regulation of MSI2 ( Fig. 2D; Table S1). Among those, 82% of the genes showed a down-regulation versus 18% to be up-regulated upon MSI2 knockdown. Previous RIP-CHIP analysis for MSI2 also revealed a similar pattern of more down-regulated genes versus up-regulated ones upon MSI2 knockdown (13).

Several growth signaling pathways are regulated by MSI2
Using our results from direct MSI2 targets, we evaluated the biological processes that are regulated by MSI2 by applying Ingenuity Pathway Analysis (IPA). The entire output of differentially regulated genes for MSI2 KD (FDR Յ0.05 and logFC Ն Ϯ 1) as well as only the genes, which were differentially expressed and showed at least one confident PAR-CLIP cluster (FDR Յ0.05, logFC Ն Ϯ 1.5, T-C conversion Ն0.25), were selected to be analyzed by IPA.
The entirety of differentially regulated genes revealed six canonical pathways to be significantly affected (FDR Յ0.05) by MSI2 (Table S2). These included but were not limited to growth factor pathways such as HGF and EGF signaling, which were activated upon MSI2 KD (z-scores of 2.4 and 1.7, respectively). The eIF2 signaling pathway was the most significantly affected pathway upon MSI2 KD and showed a total overlap of 43.6% genes of the pathway with the differentially expressed genes of MSI2 KD. It was predicted to be highly inhibited upon MSI2 knockdown with z-score of Ϫ5.217 for MSI2. The other pathways, including eIF4 (eukaryotic translation initiation factor 4) and p70S6 kinase pathway, actin cytoskeleton pathway and EGF, HGF growth factor pathways were all highly activated upon MSI2KD with positive z-scores. To investigate the downstream effects of MSI2 further, we searched for upstream targets in the direct targets list of MSI2 that are associated with these signaling pathways, and found IL6ST to be bound and transcriptionally regulated by MSI2 in a negative fashion as proposed by our genome-wide pathway analysis.

MSI2 regulates the expression of IL6ST
We eventually focused our attention onto interleukin 6 signal transducer, or IL6ST, which was also one of the most highly regulated targets. IL6ST expression was significantly increased (logFC ϭ 1.98; FDR ϭ 6.4E-17) upon MSI2 knockdown and thus negatively correlated with MSI2 expression (Fig. 3A). To

MSI2 regulates IL-6 signaling
confirm the result, we performed quantitative PCR (qPCR) on MSI2 knockdown cells and observed increased expression of IL6ST (Fig. 3B). In addition, IL6ST contained an MSI2-binding site in its 3ЈUTR, which we validated by EMSA (Fig. S1D). We also created stable knockdown of MSI2 in HEK293T cells using shRNA construct and observed an increase in the levels of IL6ST protein levels (Fig. 3C). To further confirm that MSI2 is a negative regulator of IL6ST, a control vector and a vector expressing MSI2 stably were transfected into HEK293T cells. The stable expression of MSI2 significantly decreased IL6ST mRNA levels (Fig. 3D). Similarly, we tested this regulatory link in cells with MSI2 expression under the control of an inducible tetracycline promoter. As controls, we used cells in which RBMS2 and EWSR1 expression were under the control of tetracycline promoter. Upon induction of these three cell lines with doxycycline, levels of IL6ST decreased in MSI2 cells, whereas no effect was observed in cells expressing RBMS2 and EWSR1 (Fig. 3E). To avoid cell-specific impact of MSI2 on IL6ST expression, we chose to examine Lenti GFPϩ Hep G2 cells, which also have high native IL6ST surface expression. We knocked down MSI2 in these cells and checked for surface expression of IL6ST using FACS and observed an increased surface expression of IL6ST (Fig. 3F). We also analyzed RNA-Seq data of 72 AML patient datasets and found that this reciprocal relationship between IL6ST and MSI2 holds true in AML context as well (see supporting information). Together, these data suggest that overexpression of MSI2 decreases IL6ST transcript and protein levels and vice versa. These effects on mRNA abundance of IL6ST by differential expression of MSI2 suggest that MSI2 likely regulates the stability of the IL6ST transcript.

MSI2 regulates IL6ST mRNA stability
To gain understanding of how MSI2 binding to the IL6ST 3ЈUTR affects IL6ST mRNA translation, we cloned the 3ЈUTR region of IL6ST that contains multiple MSI2 binding motifs into a luciferase reporter construct (psiCHECK2) bearing the Renilla luciferase coding region along with firefly luciferase, which served as an internal control for transfection (Fig. S2A). psiCHECK2-IL6ST (3ЈUTR) and the parent vector psiCHECK2 were transfected into HEK293Ts that expressed either normal or reduced levels of MSI2. As shown in Fig. 3G, MSI2 silencing clearly increased the luciferase expression of psiCHECK2-IL6ST, whereas parent vector levels did not change, indicating that the binding of MSI2 to IL6ST 3ЈUTR was responsible for promoting the degradation of the IL6ST.
To further test if the binding of MSI2 to IL6ST mRNA affects its stability, we performed mRNA stability assays by treating the cells with actinomycin D to inhibit de novo RNA synthesis. We found that the mRNA levels of IL6ST were higher in cells treated with MSI2 siRNA when compared with those in control cells (Fig. 3H). Together, these findings indicate that binding of MSI2 to the 3ЈUTR of IL6ST mRNA enhances the degradation of the IL6ST by destabilizing it.

MSI2 expression affects JAK/STAT and MAPK pathways
Having established the regulatory role of MSI2 on IL6ST stability and expression, we aimed to understand the molecular mechanisms modulated by this regulation. We therefore eval-uated the pathways, STAT3 phosphorylation, and ERK phosphorylation, which are induced by signal transduction through IL6ST as well as affected by EGF and HGF signaling. We checked for the phosphorylation status of STAT3 and ERK upon IL6ST overexpression in the stable MSI2 knockdown cells. For this, we stimulated the stable knockdown cells with different concentrations of IL-6 cytokine (0 -900 ng/l). IL-6 activates IL6ST which in turn leads to initiation of JAK/STAT signaling. The protein levels of IL6ST, total STAT3, and phosphorylated STAT3 were examined after IL-6 stimulation. We noticed a constant hyperphosphorylation of STAT3 upon increased IL-6 concentration whereas total protein content remained unchanged (Fig. 4A). Similarly, we also checked the protein levels of IL6ST, total ERK, and phosphorylated ERK after IL-6 stimulation and noticed similar hyperphosphorylation in ERK (Fig. 4B). These results indicated that post-transcriptional gene regulation of IL6ST by MSI2 affects the phosphorylation of STAT3 and ERK proteins, which in turn affects JAK/STAT and MAPK signaling pathways.

Discussion
In this study we investigated the RNA-binding partners of the RNA-binding protein MSI2 and confirmed the RRE by both PAR-CLIP and EMSA. Furthermore, we identified that IL6ST, a signal transducer for multiple cytokines, is bound and regulated by MSI2. We also showed that stimulating MSI2 knockdown cells with IL-6 led to hyperphosphorylation of STAT3 and ERK1/2. Thus, we identified targets bound and regulated by MSI2 and put forth a novel mechanism by which IL6ST is post-transcriptionally regulated by an RNA-binding protein.
Because IL6ST is part of IL-6 signaling mechanism, we described a novel pathway, by which this signaling mechanism can be regulated.
Although MSI2 has been implicated in the development of several aggressive cancers as well as exerting important cellular functions, the RNA-binding partners of this RNA-binding protein have not been thoroughly investigated. Applying PAR-CLIP and RNA-Seq to MSI2, we found several relevant RNAbinding partners and their respective RNA-binding motif. The high degree of conservation among MSI family members across broad ranges of animal phyla indicates that their motifs might be similar. Okano and co-workers identified a similar motif in Drosophila MSI and mouse MSI1 as (G/A)U1-3AGU (19). Also, the motif in mouse MSI1 was found to be 87% identical to human MSI2 (20). Although the motifs are identical, there is a differential binding specificity among Drosophila, mouse, and human MSI RRM family domains, suggesting a functional divergence. Despite this divergence, the core element driving the binding among MSI family has been evolutionarily conserved and was established by Zearfoss et al. as UAG using quantitative fluorescence anisotropy assay (20). Consistent with this study, our data also revealed a similar motif in human MSI2 to be U followed by UAG (UUAG) which was present in 76% of the top 100 binding sites in our PAR-CLIP datasets.
We have identified several novel targets of MSI2 which were previously shown to play an important role in hematopoiesis. Among them is KMT2A (alias ALL1, MLL), which is a transcriptional co-activator that plays an important role in hema-

MSI2 regulates IL-6 signaling
topoiesis and is also known to be frequently translocated in many leukemias (21). Another target, EGR1, was shown to be often deregulated in hematological malignancies such as AML, CML, multiple myelomas and B cell lymphomas (22). ASXL2, another target of MSI2 is a tumor suppressor which maintains normal HSC function (23,24). The functional link between MSI2 and these targets gives us a better idea about the role of MSI2 in hematopoiesis and leukemia. Similarly, targets such as SMG1 and RUNX1, which regulate cell growth in AML (25,26), might be affected upon overexpression of MSI2 in AML. This regulation therefore could explain the increased cell expansion seen in AML patients (27). Consistent with the previous finding that MSI2 modulates the TGF-␤ pathway by direct binding to genes such as TGFBR1 and SMAD3, our data add additional targets of MSI2 involved in this pathway such as RAF1, p300, CBP, SMAD5, and SMAD9 (13).
Our genome-wide pathway analysis revealed targets of MSI2 that regulate important cellular pathways such as eIF2, eIF4/ p70S6K, EGF, and HGF. Both the eukaryotic translation initiation factor pathways (eIF2, eIF4 and p70S6K) attribute a translational role to MSI2 protein. The genes included in these pathways, that are bound and regulated by MSI2, are ribosomal proteins including RPL11, RPS7, RPL19, and RPS20, which suggests that MSI2 might be an essential factor in recruiting these proteins for protein translation. Additionally, the two growth factor pathways HGF and EGF in turn regulate central signaling cascades including PI3K-Akt, JAK/STAT, or MAPK (28 -30), which have increased activity after MSI2 KD, thus being negatively regulated. The EGF and HGF signaling pathways activate the PI3K-Akt and MAPK/ERK pathways by phosphorylation/ dephosphorylation of downstream targets, indicating that the phosphorylation in general seems to be likely affected by MSI2 (31,32). Additionally, to support this idea, we have identified many targets in the pathway analysis, including JAK1, PIK3CA, ITPR2, and MAP3K1. Among several targets of MSI2, we here described a regulation of IL6ST, which is one of the upstream targets of these signaling pathways. We showed that regulation of IL6ST by MSI2 affected the phosphorylation of downstream targets such as ERK1/2 and STAT3 upon IL-6 stimulation. Because MSI2 was found to be up-regulated whereas IL6ST was down-regulated in AML patients, we hypothesize that STAT3 is hypophosphorylated upon IL-6 stimulation thus rendering the cells less susceptible to respond to immune modulatory attacks. We showed that ERK1/2 was hyperphosphorylated in stable MSI2 knockdown cells upon IL-6 stimulation, which is not in line with the previous findings that ERK phosphorylation decreased, and leads to inactivation of MAPK pathway (33). This ambiguity might be because of the cellular context in Antibodies and concentration of IL-6 are indicated. On the right is the densitometric quantification of phosphorylated STAT3 protein band, which is normalized to total STAT3 protein from three independent blots. Data were analyzed with the Student's t test and presented as the mean with S.D. (*, p Ͻ 0.05; n ϭ 3 independent blots). B, Western blotting showing the protein levels of phosphorylated ERK and total ERK in cells stably expressing scrambled (control) and MSI2 shRNA upon IL-6 treatment. Antibodies and concentration of IL-6 are indicated. On the right is the densitometric quantification of phosphorylated ERK protein band normalized to total ERK protein from three independent blots. Data were analyzed with the Student's t test and presented as the mean with S.D. (*, p Ͻ 0.05; n ϭ 3 independent blots).

MSI2 regulates IL-6 signaling
which the study was undertaken or that IL-6 -stimulated cells might react differently from nonstimulated cells.
Taken together, our study identified multiple novel target mRNAs regulated by human Musashi protein MSI2. We describe a novel mode of regulation of IL6ST (GP130) by MSI2. Our findings show that MSI2 broadly controls the post-transcriptional regulation of multiple mRNAs that control many oncogenic signaling networks, but a precise understanding of the underlying mechanisms will be vital to be able to translate this knowledge to clinical therapy.

Plasmids
Plasmid pENTR4_MSI2 was generated by PCR amplification of the respective coding sequence (CDS) followed by restriction digestion with SalI and NotI and ligation into pENTR4 (Invitrogen). Sequences were as follows: MSI2_Not_Rev, ATAGTTTAGCGGCCGCTCAATGGTATCCATTTGTAAA-GGC and MSI2_Sal_For, ACGCGTCGACATGGAGGCAAA-TGGGAGCCAAGGC. PENTR4_MSI2 was recombined into pFRT_DESTFLAGHA modified destination vector to allow for expression of FLAG/HA-tagged protein in stably transfected Flp-In TREx HEK293 cells (Invitrogen) (34).

PAR-CLIP
PAR-CLIP from Flp-In TREx HEK293 cells (Invitrogen) stably overexpressing FLAG/HA-tagged MSI2 protein was performed as described previously (16,35). Briefly, cells were grown in medium supplemented with 100 M 4SU (4-S-uridine) (Sigma-Aldrich) for 12 h prior to cross-linking. The medium was decanted and cells were washed in PBS and irradiated with 0.15 J per cm 2 total energy of 365 nm UV light (Spectrolinker XL-1500). Cells were then harvested and lysed and the cleared lysate was treated with RNase T1 (Thermo Fisher). FLAG/HA-tagged MSI proteins were immunoprecipitated with an anti-FLAG antibody (Sigma-Aldrich, F3165) conjugated to Dynabeads (Thermo Fisher). After second RNase T1 digestion, beads were washed in immunoprecipitation wash buffer and resuspended in dephosphorylation buffer and incubated with calf intestinal alkaline phosphatase beads. The cross-linked RNA was then radiolabeled with T4 polynucleotide kinase (Thermo Fisher) and washed, and the protein-RNA complex was released from the beads by incubating at 90°C and separated on SDS gels (NuPAGE Bis-Tris Midi Gels and XCell4 SureLock Midi-Cell chamber from Thermo Fisher Scientific). The protein-RNA complex was excised from the gel followed by proteinase K digestion. The leftover RNA was extracted using phenol-chloroform (Sigma-Aldrich) extraction and ethanol precipitation. The 5Ј-32 P-phosphorylated RNA14C was then subjected to a standard cDNA library preparation protocol (36). The RNA was ligated with adaptors at 3Ј and 5Ј ends followed by reverse transcription. The resulting cDNA was amplified by adaptor-specific primers and sequenced with an Illumina HiSeq 2000 (50 cycles).

siRNA knockdown, library preparation, and RNA-Seq
siRNA transfections (50 nM final concentration; duration of knockdown 72 h) of HEK293 cells were performed in 24-well format using Lipofectamine RNAiMAX (Invitrogen) as described by the manufacturer. The sequences of three different MSI2 siRNAs are given in the supporting information. For other knockdown experiments (apart from RNA-Seq), we used a mixture of these three different siRNAs. Total RNA was extracted using TRIzol (Invitrogen) and further purified using the RNeasy purification kit (Qiagen). Libraries were prepared according to manufacturer's instructions using TruSeq Total RNA sample preparation kit (Illumina) and further sequenced using Illumina HiSeq 2500 platform (single read, 100 cycles).

qRT-PCR
Total RNA was extracted using the RNeasy Mini kit (Qiagen). cDNA was synthesized using SuperScript III Reverse Transcriptase kit (Invitrogen) according to the manufacturer's instructions (oligo dT). qRT-PCR was carried out using Power SYBR Green Kit (Applied Biosystems). All reactions were run on an ABI 7500 Real-Time PCR System (Applied Biosystems) in triplicate. Data were acquired using the ABI SDS 2.0.1 software package. RNA isolated from the samples was tested for the expression levels of the chosen targets and their ⌬Ct values were subtracted from the respective beta-actin expression levels.

Recombinant protein expression and purification
A synthetized gene encoding the Musashi 2 protein containing optimized codon usages for Escherichia coli was purchased from GeneArt. Following digestion with PmlI and AvrII restriction enzymes (New England Biolabs), the gene was ligated into pET302/NT-His vector and transformed into E. coli BL21 (DE3) cells. Protein expression was induced by addition of isopropyl ␤-D-1-thiogalactopyranoside (IPTG) (to a final concentration of 1 mM) to the cells grown at 37°C, with A 600 of 0.7. After 4 h, cells were harvested by centrifugation at 5000 ϫ g. Cell pellets were resuspended in lysis buffer (500 mM NaCl, 20 sodium phosphate, 20 mM imidazol, pH 7.5, containing 5 mM DTT and EDTA-free protease inhibitor tablets (Roche Applied Sciences)) and lysed using a cell disrupter (Constant Systems). After a mild sonication step, insoluble materials were removed by centrifugation at 28,000 ϫ g. The soluble fraction was loaded onto a HisTrap FF ml column (GE Healthcare Life Sciences). The His-tagged Musashi protein was eluted with 200 mM imidazol in 20 mM sodium phosphate, 500 mM NaCl, pH 7.5 buffer containing 5 mM DTT. Eluted fractions were pooled and dialyzed against 20 mM Tris HCl (pH 7.5) containing 5 mM DTT and 1 M NaCl (molecular weight cut-off, 3 kDa) to remove RNA impurities. After an additional dialysis step against the same buffer without NaCl, the protein sample was loaded onto a cation exchange column (HiTrap SP FF 5 ml, GE Healthcare). Next, the purified protein was further purified by size-exclusion chromatography using a HiLoad 16/60 Superdex 75 (GE Healthcare), equilibrated in 20 mM Tris HCl, 150 mM NaCl, pH 8.0, containing 5 mM DTT. All purification steps were carried out with an Äkta purifier system and monitored by measuring the absorbance at 215 nm, 254 nm, and 280 nm. Protein sample analysis by 15% PAGE (SDS-PAGE) did not detect any impurities. Protein concentration was determined using UV spectroscopy.

Electrophoretic mobility shift assays (EMSA)
1 nM radioactively labeled (␥-ATP) RNA was incubated with protein concentrations varying from 0 -1 M for 1 h in 1.5 ml passivated (50 l 1 mg/ml acetylated BSA, 27°C, 1 h) prelubricated reaction tubes. The reaction was separated by native PAGE for 2 h at 300 V. The radioactive signal was visualized by phosphorimaging and the signal was quantified using Image Gauge software (Fujifilm).

Luciferase assays
200 ng of control plasmid psiCHECK2 (Promega) or psi-CHECK2 bearing the 3ЈUTR of IL6ST mRNA containing multiple binding motifs were transfected into HEK293T cells together with siRNAs for MSI2 and control using Lipofectamine 2000 (Invitrogen). 48 h later, Renilla luciferase and firefly luciferase were measured using Dual-Glo Luciferase Assay System (Promega) using Tecan 10 M Spark Luminometer.

Half-life assays
The mRNA half-life was determined by treating the cells with 3 g/l actinomycin D (Sigma-Aldrich). Cells were collected at different time points (0 -6 h) and total RNA was isolated followed by cDNA synthesis. Real-time PCR, to measure the percentage of RNA remaining for MSI2, IL6ST along with GAPDH as an internal control, was performed.

Statistical analysis
All experiments were performed at least in triplicates (biological triplicates). Numerical data were expressed as mean Ϯ S.D. Group comparisons were analyzed by two-way ANOVA or unpaired t test. p values were calculated and a p value of Ͻ0.05 was considered significant.

Bioinformatic analysis
Processing RNA-Seq reads-We obtained a total of 138,315, 336 sequencing reads for the regular RNA-Seq before and after MSI2 knockdown (each sequenced in biological triplicates). First, adapter sequences and low-quality ends of the raw sequencing reads were trimmed using seqtk (https://github. com/lh3/seqtk) 8 and cutadapt (37). All sequencing reads shorter than 25 bases after trimming were discarded. The trimmed and filtered sequencing reads were then aligned against the human genome sequence GRCh38 (38) with TopHat2 version 2.0.13 (39). Next, read counts on a gene basis were estimated by applying HTSeq (40) using gene annotations downloaded from Ensembl Genes version 77 (41). A total of 94,422,196 sequencing reads were mapped to annotated genes. The read counts represented an estimate of the respective transcript abundance and were normalized according to the sequencing library size (resulting in the so-called counts per million or cpm) using the Bioconductor package edgeR (42). Genes with low (no discernible) expression were discarded, as these were not accessible for statistical testing. Subsequently, statistical testing was performed with edgeR to identify differentially expressed genes between the controls and the MSI2 knockdown samples. The rate of differential expression is represented by the logarithm of the -fold change (logFC). To adjust for multiple testing, we used the FDR and chose genes as significantly differentially expressed if the FDR was below 5%. All sequencing data were deposited at European Nucleotide Archive (ENA), accession number PRJEB25085.
Processing of PAR-CLIP reads-The PAR-CLIP sequencing reads were preprocessed similarly to the RNA-Seq reads. However, because of some intrinsic properties compared with RNA-Seq reads (43), we kept sequencing reads Ն18 bases after trimming with cutadapt for further analyses. Next, the sequencing reads were aligned against the reference sequence of GRCh38 using Burrows-Wheeler Aligner (BWA) allowing for up to two mismatches between a single read and the reference sequence (44). All aligned reads were piled up into clusters by the PARAsuite (45) applying a hierarchical clustering algorithm. CCs were identified by filtering the obtained clusters for containing at least 10 reads, being shorter than 75 bases, showing at least 25% T-C conversions; 100% T-C conversion sites or T-C conversions reported as SNPs in the dbSNP version 142 (46) were excluded. To identify the genomic RBP-binding sites of the 8 Please note that the JBC is not responsible for the long-term archiving and maintenance of this site or any other third party hosted site.

MSI2 regulates IL-6 signaling
respective RBP, the CCs were overlapped with gene annotations downloaded from Ensembl Genes V77 using the Homer annotation tool (47).
Motif discovery-To predict the binding motifs from the PAR-CLIP data, we considered the 100 most abundant CCs in each dataset defined by the highest number of reads. The detection of the sequence motif of the RBP-binding sites was performed using MEME (48). We assumed a single occurrence of the sequence motif per cluster sequence, so we chose the one occurrence per sequence (OOPS) motif model with a minimum motif length of four nucleotides.