The DNA methylation-regulated miR-193a-3p dictates the multi-chemoresistance of bladder cancer via repression of SRSF2/PLAU/HIC2 expression

Chemoresistance hinders the curative cancer chemotherapy. To define the role of the DNA methylation-regulated microRNA (miR) genes in the chemoresistance of bladder cancer, we performed both DNA methylomic and miRomic analyses of a multi-chemosensitive (5637) versus a multi-chemoresistant (H-bc) cell line and found that miR-193a-3p is hypermethylated/silenced in 5637 and hypomethylated/expressed in H-bc cells. A forced reversal of its level turned around the chemoresistance in the cultured cells and the tumor xenografts in nude mice. Three of its targets: SRSF2, PLAU and HIC2, work in concert to relay the miR-193a-3p's impact on the bladder cancer chemoresistance by modulating the activities of the following five signaling pathways: DNA damage, Notch, NF-κB, Myc/Max, and Oxidative Stress. In addition to the mechanistic insights in how the newly identified miR-193a-3p/SRSF2,PLAU,HIC2/five signaling pathway axis regulates the chemoresistance of bladder cancer cells, our study provides a new set of diagnostic targets for the guided personalized chemotherapy of bladder cancer.

Bladder cancer (BCa) is the second most common and the most life-threatening urogenital cancer in men. 1,2 In addition to the high recurring rate, BCa is highly refractory to drug therapy and regarded as one of the most costly and difficult type of cancer to be contained. 3 Although an eradication of primary lesions by chemotherapy is achievable in about half of the BCa patients, the diseases often recur and become refractory to the second round of chemotherapy. 4 Once the chemoresistance is established, the cancer cells are resistant to multiple chemo-therapeutics, irrespective of the nature of the previous drug exposure. 5 As other hallmarks of cancer, the chemoresistance varies drastically among the patients suffered from any given type of cancer, making both the chemoresistance prediction and the rational personalized chemotherapy very difficult.
In addition to the extensive genetic abnormality, cancer cells suffer from a wide range of the epigenetic defects that are more tightly linked to the gene expression pattern and therefore the cancerous phenotypes. The following entities constitute the epigenetic interface of the human genome: DNA methylation, histone modification, noncoding RNA (including microRNAs (miRs)), the nucleosome organization, and high-order chromatin structure. The miRs that aberrantly expressed in cancer make a significant contribution to cancerous behaviors, including chemoresistance. 6 Taking the miRs in the chemoresistance state of BCa as examples, the following miRs have been implicated: miR-34a, 7 miR-21, 8 miR-30d, miR-181, and miR-199a-5p. 9 DNA methylation is the best-characterized epigenetic mechanism, underlying the faithful cross-cell-division transmission of the gene transcription memory. 10 The hypermethylated state of the promoter and enhancer regions tightly correlates with the transcriptionally silenced state of both protein-coding and non-coding RNA genes. Therefore determining of the DNA methylation state of the promoter regions, rather than the level of the corresponding RNAs or proteins, in patient samples promises a better way for both early detection and rationale personalized therapy of cancer. Approximately, 10% of the miR genes are co-localized with the CpG-rich region (CpG islands) in genome, the expression of which is negatively regulated by DNA methylation. Hence we confined our attention to this subset of miRs for their roles in BCa's chemoresistance. In this report, we began with a bi-omic (miRomic and methylomic) analysis of a multi-chemosensitive (5637) versus a multi-chemoresistant BCa cell line (H-bc) for the differentially expressed/methylated miR genes. The miR-193a-3p, one of the most differentially expressed/methylated CpG island containing miRs identified, was subjected to a systematic investigation for its role and the underlying mechanisms in the BCa's multi-chemoresistance.

Result
Expression of the DNA methylation-regulated miR-193a gene positively correlates with the multi-chemoresistance of BCa cells. The drug dose for 50% cells killed by the following drugs: Pirarubicin (Pi), Paclitaxel (Pa), Adriamycin (Ad), and Epirubicin Hydrochloride (EH) after a treatment of 72 h, was determined in the following five BCa cell lines: 5637, T24, EJ, H-bc, and Biu87 (Figures 1a and b). Judged by the fold difference over the lowest IC 50 , 5637 was the most multi-chemosensitive cell line, with the lowest IC 50 to four of five drugs, while H-bc was the most resistant cell line with its relative IC 50 by 8.32 to 36.96 folds higher than 5637 cells (Figure 1c). From a RNA-seq based miR-omic analysis of 5637 and H-bc cell lines (Figure 2b and data not shown), 83 miRs were found differentially expressed by no less than two folds: 37 higher and 45 lower in 5637 than H-bc cells, respectively. Among the 20 miR genes that are colocalized with the CpG island, the miR-193a-3p differentially expressed between 5637 and H-bc cells at the top range in both RNA-seq omic ( Figure 2b) and a qRT-PCR analysis (Figures 2c and d). The hypermethylated state in 5637 and hypomethylated state in H-bc cells of the miR-193a gene suggested by the methyl-capture seq analysis 11 (Figure 2e, omic data not shown) were confirmed by a bisulfite conversion sequencing (BSP) analysis (Figures 2f and g): It is hypermethylated in 5637 (the average CpG methylation: 78.6%) and barely methylated in H-bc (the average CpG methylation: 4.3%). In conclusion, the differential state of the miR-193a gene at both DNA methylation and expression levels tightly correlates with the multi-chemoresistance of BCa cells.
SRSF2, PLAU and HIC2 are direct targets of miR-193a-3p in BCa cells. The number of the genes potentially regulated by one miR ranges from several hundreds to a couple of thousands. Among the predicted target genes of miR-193a-3p that were picked up by no less than three out of four commonly used predicting methods: miRDB (166 entries identified), miR base (380 entries), targets can (443 entries), target miner (243 entries), HIC2 (hypermethylated in cancer 2) gene, along with two known miR-193a-3p's targets: SRSF2 (serine/arginine-rich splicing factor 2) 12 and PLAU (plasminogen activator, urokinase), [13][14][15] were found to be expressed in an opposite manner of the miR-193a-3p at both RNA (RNA-seq based omic: Figure 3b and qRT-PCR analysis, Figure 3c) and protein levels ( Figure 3d). Furthermore, a miR-193a-3p mimic transfection brought down the levels of all three genes in 5637 cells, and their levels were raised in the antagomiR-transfected H-bc cells (Figures 3e  and f).
For the proof that HIC2 and PLAU genes are direct targets of miR-193a-3p, both 3 0 -UTR regions were put at the downstream of the firefly luciferase gene in pGL3 (Promega, Madison, WI, USA) to create pGL3-HIC2 UTR and pGL3-PLAU UTR constructs ( Figure 3g). Both constructs and pGL3 were transfected into 5637 and H-bc cells, respectively, for the functional state of miR-193a-3p in cells. Both UTR-contained constructs gave a significantly higher luciferase activity in 5637 than H-bc cells in an opposite pattern of miR-193a-3p's expression, in contrast to the no difference state of the pGL-3resulted luciferase activities in both cell lines (Figures 3h and  i). Moreover, the miR-193a-3p mimic transfection repressed the luciferase activities of both 3 0 -UTR containing but not pGL-3 constructs in 5637 cells, and the reverse was found in the antagomiR-transfected H-bc cells (Figures 3h and i). In conclusion, similar to SRSF2, 12 both PLAU and HIC2 are the true direct targets of miR-193a-3p in BCa cells. The miR-193a-3p's effect on the BCa chemoresistance is likely to be realized by its repression of these three genes.  Table   miR-193a-3p regulates bladder cancer drug resistance L Lv et al miR-193a-3p rather than -5p dictates the multi-chemoresistance of BCa cells. We have previously reported the positive involvement of miR-193a-3p via repression of SRSF2 expression in the 5-FU chemoresistance of hepatocellular carcinoma (HCC). 12 The miR-193a-5p's regulatory role in both cell survival and chemosensitivity was also reported, probably via repression of TP73 expression. 16 Both SRSF2 and TP73 proteins were indeed significantly lower in H-bc than in 5637 cells (not shown) and were upregulated by the mimic transfection in 5637 and downregulated by the antagomiR transfection in H-bc cells (Figures 4c and e), well matching the level (functional state too) of miR-193a-3p and -5p ( Figure 1). A miR-193a-3p but not -5p mimic transfection are capable to reduce the drug-triggered cell death in 5637 cells by 0.2-1.3-fold in a drug-specific manner (Figure 4b).
The same conclusion was drawn from the analysis in the antagomiR-transfected H-bc cells ( Figure 4d). Therefore, miR-193a-3p but not -5p contributes to the BCa's multi-chemoresistance, although both express in a same pattern. We then individually repressed each miR-193a-3p's target genes by a siRNA transfection (Figure 4g) in 5637 cells, followed by the measurement of the Pa-or Ad-triggered cell death. With the miR-193a-3p mimic transfection as a reference, a compatible impact was observed. Repression of either HIC2 or PLAU level gave a stronger effect than the repression of SRSR2 gene (Figure 4f). The siRNA-mediated co-repression of both SRSF2 and PLAU or all three genes caused a further reduction of Pa-triggered but not Ad-triggered cell death ( Figure 4l). Intriguingly, only the co-repression of both SRSF2 and HIC2 led to a further reduction of the Ad-triggered (but not Pa-triggered) cell death  (Figures 4h and k). Therefore, the synergistic effect of miR-193a-3p's target genes on the BCa's chemoresistance was in both gene-specific and drug-type-specific fashion. Conversely, a forced upregulation of both SRSF2 and PLAU proteins by transfection of the constructs where the GFP-tagged genes are under the potent mammalian promoter sensitized H-bc cells to the Ad-triggered cell death, but neither mock nor the GFP vector transfection did. Compatible to the 3PA's effect, the forced upregulation of both genes had no effect on the Pa-triggered cell death (Figure 4i). An expected expression of the exogenously introduced SRSF2-GFP and HIC2-GFP proteins in H-bc cells was confirmed by examination under the fluorescent microscope (Supplementary Figure S1B). In conclusion, the miR-193a-3p's effect on the BCa chemoresistance is indeed significantly accomplished by repression of these three genes, the subset of genes involved contributing to the drug-specific nature of the BCa chemoresistance.   Compatible to the effect on the drug-triggered cell death, the influence of the miR-193a-3p's target genes on the pathways are in a gene-specific manner.
For the further mechanistic understanding, we searched for the interactions among the three miR-193a-3p's targets (PLAU, SRSF2 and HIC2 genes) and the master transcription factor genes for these five signaling pathways (Figure 5d): NRF1 for Oxidative stress, MYC for Myc/Max, RBPJ for Notch, TP53 for DNA Damage and NFKB1 for NF-kB pathways from the STRING database of the known proteinprotein interactions (PPI) interactions. 17 Using one-step direct path algorithm, a direct link between PLAU and TP53 (DNAdamage response pathway) or NFKB1 (NF-kB pathway) was  (Figure 5e). Via both literature mining and the expression analysis of these eight genes in the RNA-seq data sets of 5637 and H-bc cells (not shown), we built a more sophisticated connection network. Three hub genes connecting both SRSF2 and PLAU to NRF1, RBPJ and TP53 genes, respectively, were all involved in the regulation of the acetylation state of both histones and non-histone proteins: PCAF (histone acetyltransferase b) and two histone deacetylases (HDAC6 and HDAC9). The hub gene connecting HIC2 protein to MYC (Myc/Max) pathway is CCNT1, encoding a CDK9-associated C-Type protein with a positive regulatory role in transcript elongation. 18 The recent finding that HIC2 forms a complex with SRSF2 via both CCNT1 and Cdk9 19-21 may explain why both HIC2 and SRSF2 genes share the repressing role on both Notch and Oxidative pathways.  Figures  6c and d). Again, neither miR-193a-5p agomiR nor antagomiR altered the growth of the Pa-treated 5637 or H-bc tumors. Further confirmation of the miR-193a's role in Pa resistance of BCa came from the immuno-histological analysis of SRSF2 (a miR-193a-3p's target), TP73 (a -5p's target) and Ki67 (an indicator for cell proliferation) in the tumor sessions of the Pa-treated versus PBS-treated mice (Figure 6e and Supplementary Figure S4). The intratumor injection of either miR-193a-3p or -5p's agomiR (into 5637)/antagomiR (into H-bc tumor) indeed led the expected changes of both SRSF2 or TP73 levels in tumor slides (Figure 6e), which consolidate the conclusion that miR-193a-3p but not miR-193a-5p has a profound positive effect on both the growth and chemoresistance of the BCa cell-derived tumor xenografts in nude mice.

Discussion
Aberrant miR-193a expression has been reported in all the types of cancer examined, including colorectal cancer, 22 non-small cell lung cancer (NSCLC), 23 myeloid leukemia 24 and Wilms' tumor blastema. 25 The transcription factors implicated in the regulation of miR-193a transcription include XB130, 26 p63, 16 Max and RXR. 13 As a typical CpG island-containing gene, its transcription is also under a regulation of the DNA methylation state of the promoter region and activated by a 5' aza-2'-deoxycytocine-mediated demethylation. 12 The tumorsuppressor role of miR-193a-3p has been reported in NSCLC 27 and epithelial ovarian cancer cells. 28 On the other hand, miR-193a-3p can also promote both in vivo growth and chemoresistance of HCC 12 and of BCa in the present study. The functional disparity in cancer biology of miR-193a-3p in literature is likely attributed to the system difference of studies concerning the type of cancer, the experimental designs and whether the in vivo studies involved. It is also very likely that miR-193a-3p can act as both a tumor suppressor and oncogene in a biological and pathological content-dependent manner.
A miR executes its biological function via repression in a sequence-specific manner of up to a couple of thousands protein-coding genes at both stability and translation levels of mRNAs. In this study, we have defined the roles of PLAU, HIC2 and SRSF2 genes, three direct targets of miR-193a-3p, that work in concert to relay the miR-193a-3p's impact on the BCa multi-chemoresistance ( Figure 4). PLAU encodes the urokinase-type plasminogen-activator protein and actively involves in cell migration, invasion and angiogenesis processes. Its level is much higher in an Ad-resistant than its sensitive counterpart in breast cancer. 29 There was a report that the siRNAmediated repression of PLAU sensitizes HEK-293 cells to the UV-, cisplatin-and detachment-triggered apoptosis through DNA-damage pathway. 30,31 However, the opposite was found in our study, that is, a siRNA-mediated repression of PLAU enhanced rather than repressed the multi-chemoresistance of 5637 cells (Figure 4). HIC2 encodes a putative transcriptional factor, the genetic aberration of which is associated with both digeorge syndrome, 32 and systemic lupus erythematosus 33 in humans. However, its role in chemoresistance of cancer is, for miR-193a-3p regulates bladder cancer drug resistance L Lv et al the first time, described here. SRSF2 protein is a member of the serine/arginine (SR)-rich family of pre-mRNA splicing factors, a critical component of the spliceosome executing both splicing and nuclear-cytoplasmic transportation of mature mRNA. 34 Its role in the miR-193a-3p's effect on the 5-FU chemoresistance in HCC 12 is essentially confirmed in the BCa multi-chemoresistance in the present study. We also demonstrated a genespecific as well as a drug-type specific nature of the additive (or synergistic) effect of these three genes on the BCa chemoresistance (Figure 4h).
Although the signaling pathways affected by the target therapeutics are well established, the pathways challenged by the conventional chemotherapeutics remain undefined. Pa, a cytoskeleton drug is expected to disturb different pathways from Ad, a genotoxic drug. The assumption that a key set of pathways involved in the BCa chemoresistance have significantly different activities between the multi-chemosensitive 5637 and the multi-chemoresistant H-bc cells was confirmed; 7 of the 16 signaling pathways tested displayed different activities by more than two-fold in this pair of BCa cells (Figure 5b), 5 of which were regulated by miR-193a-3p ( Figure 5c) and its target genes (Figure 5d) in an expected fashion. By a network analysis using three miR-193a-3p's target genes and five master transcription factors as cues, a number of significant interactions among these eight genes have been found (Figure 5e). HIC2's capability to form a complex with SRSF2 via a physical interaction with both CCNT1 and Cdk9 19-21 may explain why these two miR-193a-3p's target genes repress both Notch and Oxidative pathways in a similar manner. In addition, HDAC6 deacetylates, stabilizes SRSF2 protein and therefore enhances its function. 35 Three hub genes connecting SRSF2 to the master transcription factors RBPJ (Notch) and NRF1 (Oxidative stress) are two histone deacetylase (HDAC6 and 9) and one histone acetyl-transferase (PCAF). These three proteins work in concert to delicately control the acetylation state of both core histones and non-histone proteins. 36 PCAF can also interact with RBPJ 37 or NRF1 38 for an activation of both pathways. HIC2, but not SRSF2, acts as a positive regulator of MYC/Max pathway, probably via HIC2's interaction with CCNT1 and in turn activates the positive transcription elongation factor b (P-TEFb) 21 . 39 The transcription of the miR-193a gene is regulated by Max and RXR. 13 The regulatory role of both HIC2 and PLAU on MYC/max pathway in the context of the multi-chemoresistance of BCa cells (Figures 4 and 5) suggests a negative regulating loop between the miR-193a-3p and MYC/Max pathway.
As summarized in Figure 7, we have identified a miR-193a-3p centered axis that dictates the BCa's multi-chemoresistance, expression of which was under the negative control of the DNA methylation. Via its repressive effect on three target genes (and others) and in turn the five signaling pathways, the miR-193a-3p promotes the multi-chemoresistance (Pa, Pi, EH and Ad studied in this report) and the in vivo tumor growth of BCa cells. For further mechanistic insights, we have informatically identified the key links that connect miR-193a-3p via three of its target genes (SRSF2, PLAU and HIC2) to the five signaling pathways and therefore the chemoresistant phenotype of BCa cells (Figures 5e and 7b). Our study provides a new set of diagnostic targets for the guided personalized chemotherapy of BCa. The mimic/antagomiR/siRNA/overexpression plasmids transfection: all the mimic, antagomiR, siRNA, and the scramble sequence control (NC) as well as riboFECT CP transfection kit were supplied by Guangzhou Ribobio (Guangzhou, China). The mammalian expression constructs for PLAU (EX-F0073-M98-5) and SRSF2 (EX-F0622-M98-5) with GFP tag were supplied by Guangzhou Fulengen (Guangzhou, China). Transfection of both ribonucleic acid reagents or plasmids mentioned above and the reporter plasmids in a Cignal Finder Pathway Reporter package (Qiagen, Hilden, Germany) was performed according to the manufacturer's instruction.
The luciferase reporter assay: A full length of the human PLAU 3 0 -untranslated region (932 bp) and partial length of HIC2 3 0 -UTR (1131 bp, 1862-2992 from a 4719-bp full length) with the miR-193a-3p targeting sequence were cloned at the downstream of the firefly luciferase gene in pGL3 (Invitrogen) to construct pGL3-luc-PLAU and pGL3-luc-HIC2, respectively. All the constructs were confirmed by restriction digestion.
Cells were seeded into 96-well plates at around 1 Â 10 4 cells per well and transfected with a mixture of 50 ng pGL3-luc-PLAU or pGL3-luc-HIC2, 5 ng Renilla plus 5 pmol mimic or scramble control (NC) nucleotides, with the riboFECT CP transfection reagents according to the manufacturer's instruction. Both firefly and Renilla luciferase activities were measured 18 h after transfection by the Dual-Luciferase Reporter Assay System (Promega) using a Promega GloMax 20/20 luminometer. The relative firefly luciferase activities were normalized with the Renilla luciferase activities as a for transfection efficiency.
The pathway luciferase reporter constructs (Supplementary Figure S2A): (1) the negative control construct: the firefly luciferase gene is under the control of the minimal promoter. (2) the pathway reporter construct: a tandem repeat of the cognate consensus motif that is recognized by each master transcription factor for the corresponding pathway was placed at the upstream of the minimal promoter in the construct 1. (3) the positive control construct: the firefly luciferase gene is under the control of CMV promoter and (4) the internal control construct. The firefly luciferase gene in construct 3 was replaced with the Renilla luciferase gene. The analysis was carried out according to the manufacturer's instruction (Qiagen). 42 Briefly, the cells were transfected in triplet with each firefly luciferase reporter construct in combination with the Renilla luciferase construct using ribo FECT CP transfection reagent, and both luciferase activities in cell extracts at 24 h after transfection were measured by a Promega Dual-Luciferase Reporter assay (Promega) using a Promega GloMax 20/20 luminometer. 43 Firefly luciferase activities from each set were normalized to the activity of Renilla luciferase to control the inter-transfection bias. The relative luciferase activities (luciferase unit) of the pathway reporter over the negative control in the transfected cells were calculated as a measurement of the pathway activity.
Chemotherapeutics: All the chemotherapeutic drugs used are of the clinic grade 44,45  Chemoresistance profiling (IC 50 measurements): Cells in the logarithmic phase of growth were seeded in triplicate in 96-well plates at the density of 0.5 Â 10 4 /well and treated with four-fold serially diluted drugs for 72 h. Cell survival was then measured by a thiazolyl blue tetrazolium bromide (MTT, 490 nm reading)-based cell proliferation assay. 46 Both the linear regression parameters and the IC 50 (the concentration of drug required for 50% of cells to be killed) with the no-drug control as the reference were calculated. The relative chemoresistance was presented as the fold for each of the cell line over the lowest IC 50 .
RNA analysis. Total RNA was isolated using the TRIzol reagent (Tiangen Biotech Co., Ltd., Beijing, China). For mRNA analysis, the cDNA was made from total RNA by oligo-dT priming with a primeScript RT reagent kit (Tiangen Biotech Co., Ltd.) and the mRNA level of the following three genes (SRSF2, HIC2 and PLAU) was measured by qRT-PCR with gene-specific fluorescent Taqman probe together with the b-actin using a different fluorescence-labeled probe (provided by ShingGene, Shanghai, China) in the FTC-3000P (FUNGLYN BIOTECH INC, Toronto, ON, Canada). For miR analysis, the cDNA was synthesized with the specific stem-loop primer and quantified by SYBR Green-based real-time PCR assay. Using the 2 À DDCt method, the normalization with the U6 reads for miR or with b-actin for mRNA was performed before each's relative level between 5637 and H-bc was calculated. 47  Western blotting analysis. Cell lysates with 1 Â SDS loading buffer (60 mM Tris-HCl, pH6.8, 2% SDS, 20% glycerol, 0.25% bromphenol blue, 1.25%2mercaptoethanol) were incubated at 100 1C for 10 min to facilitate the sample loading for the conventional western blotting analysis. The anti-SRSF2 (AP2800a), anti-HIC2 (AP18558c), anti-PLAU (AP8161b) and anti-GAPDH (AM1020a) were provided by Wuxi phama, Shanghai, China. The target proteins were then probed with anti-rabbit IgG peroxidase-conjugated antibody (LP1001b), or HRP goat antimouse IgG antibody (LP1002a) (All antibodies are from Abgent, San Diego, CA, USA) followed by an enhanced chemiluminescence reaction (Thermo Fisher Scientific, Waltham, MA, USA). The relative levels of proteins were quantified using densitometry with the Gel-Pro Analyzer (Media Cybernetics, Rockville, MD, USA). The target bands over the GAPDH band were densitometrically quantified and indicated under each band.
BSP analysis: Genomic DNA was isolated by a standard phenol/chloroform purification method, qualified by electrophoresis on an agarose gel, treated by ammonium bisulfate-based bisulfite conversion method 48,49 and the PCR fragments from the converted DNA was sequenced and analyzed. Raw sequence data files was processed and the area ratio (%) of C over C þ T of the primary CpG dinucleotide was calculated as the percenatge of methylation and plotted. 50 The in vivo studies. Animal experiments were undertaken in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals. BALB/c male nude mice aged of 8-12 weeks were used for this study. 5637 or H-bc cells were embedded in BD Matrigel Matrix (Becton, Dickinson, Franklin Lakes, NJ, USA) 51 and subcutaneously injected into the four sites at the back of mice as following: 1.7 Â 10 7 cells/site for 5637, 0.7 Â 10 7 cells/site for H-bc, 4 sites/mouse, 6 mice for 5637, 6 mice for H-bc, respectively. From the fourth day after cell injection, all 5637-generated tumors were intratumorally injected with 2nM miR-193a-3p agomiR or miR-193a-5p agomiR, while H-bcgenerated tumors were injected with 4 nM miR-193a-3p/-5p/Mock antagomiR/PBS in a similar manner. From the sixth day after cell injection, three mice from 5637 and three from H-bc intraperitoneally received Pa (45 ug/ mouse) once in 2 days. The remaining six mice (three from 5637 and three from H-bc) received PBS as a miR-193a-3p regulates bladder cancer drug resistance L Lv et al mock treatment control. Mice were humanely killed on day 25, and the tumors were weighed and photographed. The tumor weight was described as the mean ± S.D. Expression levels of SRSF2, TP73 and Ki67 proteins were measured using immunochemical analysis on 5-mm slices of formalin-fixed paraffin-embedded tumor xenografts in nude mice. To avoid inter-treatment bias, the tissue slides from all the six groups were made on a single slide and subjected to the same immunostaining simultaneously. Antigens were retrieved by pretreating dewaxed sections in a microwave oven at 750 Watts for 5 min in a citrate buffer (pH 6) processed with the Super Sensitive Link-Labeled Detection System (Biogenex, Menarini, Florence, Italy). The enzymatic activities were developed using 3-amino-9-ethylcarbazole (Dako, Milan, Italy) as a chromogenic substrate. Following counter staining with Mayer hematoxylin (Invitrogen), slides were mounted in an aqueous mounting medium (glycergel, Dako).Pictures were taken using the LEICA DM 4000B microscope (Wetzlar, Germany), while the relative level of each protein was calculated using the LEICA software (Wetzlar, Germany), and the percentage of the mock over the chemotherapeutic treated tumors was calculated and plotted.
Bioinformatics analysis. R package pheatmap version 0.77 (http://cran. r-project.org/web/packages/pheatmap/index.html) were used to create the diagrams. CpG island information for the human genome (hg19) was obtained from UCSC database. 52 Differentially expressed miRs were calculated with in-house developed perl scripts. The networking analysis: the following eight genes: PLAU, HIC2, and SRSF2, together with the master transcription factor genes in each pathways (Oxidative stress: RBPJ, MYC/Max: MYC, Notch: NRF1, DNA damage: TP53 and NF À kB pathways: NFKB1) were used as seeds to extract from the known PPI from the STRING database version 9.10. 17 The build-in one-step expand algorithm was used until a fully connected network of all seed genes were obtained. To simplify the model, the shortest paths between seeds were extracted from the networks. RNA-seq data (not shown) were used to identify the expression pattern of the model. To obtain the true functional role of each interaction in the connected network of eight seeds and four hub genes, the literature mining was performed by the STRING via full text reading, with both standard gene name and alternative name considered. Edges that are not consistent with our experimental data were removed from the network.
Statistical analysis. Data are presented as means, and error bars indicate the S.D. or S.E. All statistical analyses were performed with Excel (Microsoft, Redmond, WA, USA) or Prism (GraphPad Software Inc., La Jolla, CA, USA). Twotailed Student's t-test, a one-way analysis of variance or Mann-Whitney U test was used to calculate statistical significance. A P-value of o0.05 was considered to be significant.

Conflict of Interest
Dr. Xiao Liu is the founder of MHBI (Shanghai) Biotechnology Incorporation. The other authors declare no conflict of interest.