Deregulation of the miR-16-KRAS axis promotes colorectal cancer

KRAS plays a significant role in the etiology and progression of colorectal cancer (CRC), but the mechanism underlying this process has not been fully elucidated. In this study, we found that the KRAS protein levels were higher in CRC tissues than in the normal adjacent tissues, whereas its mRNA levels varied irregularly, suggesting that a post-transcriptional mechanism is involved in the regulation of KRAS. Then, we performed bioinformatic analyses to search for miRNAs that potentially target KRAS. We predicted and experimentally validated that miR-16 directly recognizes the 3′-UTR of the KRAS transcript and regulates KRAS expression. Furthermore, the in vitro results showed that the repression of KRAS by miR-16 suppressed the proliferation and invasion and induced the apoptosis of CRC cells, and the in vivo results revealed that miR-16 exerted a tumor-suppressive effect by negatively regulating KRAS in xenograft mice. Taken together, our findings provide evidence supporting the role of miR-16 as a tumor suppressor in CRC by targeting KRAS.

RNA isolation and quantitative RT-PCR. Total RNA was extracted from the cultured cells and human tissue specimens using Trizol reagent (Invitrogen, CA, USA) according to the manufacturer's instructions. Assays to quantify miRNAs were processed using TaqMan miRNA probes (Applied Biosystems, CA, USA) according to the manufacturer's instructions. Briefly, 1 μ g of total RNA was reverse-transcribed to cDNA using a stem-loop RT primer (Applied Biosystems) and AMV reverse transcriptase (TaKaRa, Dalian, China). The reaction conditions were as follows: 16 °C for 30 min, 42 °C for 30 min and 85 °C for 5 min. Quantitative real-time PCR was processed using a TaqMan PCR kit on an Applied Biosystems 7300 Sequence Detection System (Applied Biosystems). The reaction conditions were 95 °C for 5 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. All the reactions were processed in triplicate. After all reactions were completed, the cycle threshold (C T ) data were gathered using fixed threshold settings, and the mean C T values were determined from triplicate PCRs. U6 snRNA was used as an internal control, and the amount of miRNA normalized to the U6 level was calculated using the formula 2 −ΔΔCT , in which Δ Δ C T = (C TmiRNA − C TU6 ) tumor − (C TmiRNA − C TU6 ) control .
Overexpression or knockdown of miR-16. Overexpression of miR-16 was achieved by transfecting cells with pre-miR-16 (a synthetic RNA oligonucleotide mimicking the miR-16 precursor). Knockdown was achieved by transfecting a miRNA inhibitor (a chemically modified single-stranded antisense oligonucleotide designed to specifically sequester the mature miRNA). Synthetic pre-miR-16, anti-miR-16 and scrambled negative control RNA (pre-miR-control and anti-miR-control) were purchased from GenePharma (Shanghai, China). SW480, HT-29 and Caco2 cells were seeded in 6-well plates. On the following day, the cells were transfected using Lipofectamine 2000 (Invitrogen) when the cells were approximately 70% confluent. In each well, 100 pmol of pre-miR-16, pre-miR-control, anti-miR-control and anti-miR-16 were used. After 4-6 h, the media were changed to RPMI-1640 or RPMI-DMEM medium supplemented with 2% FBS. The cells were harvested at 24 or 48 h after transfection for the isolation of total RNA or protein, respectively.

Plasmid construction and siRNA interference assay. A mammalian expression plasmid designed to
specifically express the open reading frame (ORF) of human KRAS without the miR-16-responsive 3′ -UTR was purchased from GeneCopoeia (Germantown, MD, USA). An empty plasmid served as a negative control. Four siRNA sequences targeting different sites of human KRAS cDNA (si-KRAS) were synthesized by GenePharm. siRNA sequences were designed according to previous reports 29,30 . A scrambled siRNA (GenePharma) was used as a negative control. The siRNA sequences were as follows: si-KRAS#1: 5′ -GGUGACUUAGGUUCUAGAUTT-3′ ; si-KRAS#2: 5′ -GGAAGCAAGUAGUAAUUGATT-3′ ; si-KRAS#3: 5′ -CGAAUAUGAUCCAACAAUATT-3′ ; and si-KRAS#4: 5′ -UAAGGACUCUGAAGAUGUATT-3′ . The KRAS overexpression plasmid or siRNA was transfected into SW480 cells using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions. Total RNA or protein was isolated 24 or 48 h post-transfection. The KRAS mRNA and protein expression levels were assessed by quantitative RT-PCR and Western blot, respectively. Luciferase reporter assay. The entire 3′ -UTR of human KRAS was amplified via PCR with human genomic DNA as a template. Then, the PCR products were inserted into the pMIR-REPORT plasmid (Ambion, USA), and successful insertion was confirmed by DNA sequencing. To evaluate the binding specificity, the sequences that interacted with the seed region of miR-16 were mutated (from TGCTGCT to ACGACGA), and the mutant KRAS 3′ -UTR was inserted into an equivalent luciferase reporter. For luciferase reporter assays, SW480 or Caco2 cells were seeded in 24-well plates and co-transfected with 1 μ g of the firefly luciferase reporter plasmid, 1 μ g of the β -galactosidase (β -gal) expression plasmid (Ambion), and equal amounts (100 pmol) of pre-miR-control, pre-miR-16 or anti-miR-control, anti-miR-16 using Lipofectamine 2000 (Invitrogen). A β -gal plasmid was used as a transfection efficiency control. The cells were harvested 24 h post-transfection and were assayed with a luciferase assay kit (Promega, Madison, WI, USA).
Protein isolation and Western blot. Cells or tissues were lysed in RIPA lysis buffer (Beyotime, China) containing a protease inhibitor cocktail (Thermo Scientific 78440) for 30 min on ice and were then centrifuged at 12,000 × g at 4 °C for 10 min. The supernatant was collected, and the protein concentration was calculated using a BCA protein assay kit (Thermo Scientific, Rockford, USA). The protein levels were analyzed via Western blotting using the corresponding antibodies and normalized by probing the same blots using a GAPDH antibody. The antibodies against KRAS and GAPDH were purchased from Santa Cruz Biotechnology (sc-30  Cell proliferation assay. To assess cell proliferation, SW480 and Caco2 cells were seeded in triplicate at a density of 1 × 10 4 cells per well in 96-well plates and then incubated overnight in 100 μ l RPMI-1640 or RPMI-DMEM medium supplemented with 10% FBS. The cell proliferation index was measured using the Cell Counting Kit-8 (CK04-500, Dojindo, Japan) at 12, 24, 36, 48 and 60 h after transfection according to the manufacturer's instructions. All experiments were performed in triplicate.

Cell invasion assay. The invasion ability of SW480 and Caco2 cells was tested using Matrigel Invasion
Chambers (BD Biosciences, MA, USA) with inserts containing an 8-μ m pore-size membrane. The membranes on the bottom of the upper compartment of the Transwells were coated with a thin layer of Matrigel. The cells were harvested 24 h after transfection, suspended in FBS-free RPMI-1640 or RPMI-DMEM culture medium and added to the upper chamber (4 × 10 4 cells per well). At the same time, 0.6 ml of RPMI-1640 or RPMI-DMEM with 20% FBS was added to the lower compartment. The cells were allowed to invade for 24 h in a 5% CO 2 atmosphere that was saturated with H 2 O at 37 °C. After incubation, cells that entered the lower surface of the filter membrane were fixed with 4% paraformaldehyde for 20 min at room temperature, washed 2 times with 1 × PBS and stained with 0.5% crystal violet solution for 15 min at room temperature. The cells remaining in the upper surface of the filter membrane (non-migrant) were scraped out gently with a cotton swab. The lower surfaces (with cells that had migrated) were imaged using a photomicroscope (10 × fields per chamber) (BX51, Olympus, Japan). All experiments were performed in triplicate.
Apoptosis assays. The apoptosis of SW480 and Caco2 cells was tested using an Annexin V-FITC/propidium iodide (PI) staining assay. The SW480 cells were cultured in 12-well plates and transfected with pre-miR-16, KRAS siRNA or the KRAS overexpression plasmid and the Caco2 cells transfected with anti-miR-16 to induce apoptosis. Pre-miR-control, anti-miR-control, control siRNA and a control plasmid served as negative controls. The cells were cultured for 48 h in FBS-free RPMI-1640 medium or RPMI-DMEM, and then the attached and floating cells were harvested. Flow cytometric analysis of the apoptotic cells was performed using an Annexin

Establishment of tumor xenografts in mice.
Five-week-old male C57BL/6J nude mice were purchased from the Model Animal Research Center of Nanjing University (Nanjing, China) and housed under specific pathogen-free conditions at Nanjing University. SW480 cells were infected with the miR-16 overexpression lentivirus alone, transfected with the KRAS overexpression plasmid alone, or co-transfected with the miR-16 overexpression lentivirus and the KRAS overexpression plasmid. After infection and transfection, SW480 cells were subcutaneously injected into C57BL/6J nude mice (2 × 10 6 cells per mouse, 6 mice per group). Mice injected with untreated SW480 cells served as controls. The mice were sacrificed one month later. The tumors were separated from the animals, and the tumor weights were measured. Parts of the tissues were used for protein and total RNA extraction, and the remaining tissues were fixed in 4% paraformaldehyde for 24 h and then processed for Hematoxylin and eosin (H&E) staining and immunohistochemical staining for KRAS and Ki-67. All experiments were approved by the Institutional Review Board of Nanjing University (Nanjing, China) and performed in accordance with the guidelines of the National Institutes.
Statistical analysis. All the Western blot images are representative of at least three independent experiments. The quantitative RT-PCR assays, luciferase reporter assays, proliferation assays, invasion assays and apoptosis assays were performed in triplicate, and each experiment was repeated several times. The results are shown as the means ± SD of at least three independent experiments. The numerical data were statistically analyzed by two-tailed Student's t-test. Bivariate correlation between two independent variables was calculated by Spearman's rank correlation co-efficient. All the differences were considered statistically significant at p < 0.05.

Results
Upregulation of the KRAS protein, but not mRNA, in CRC tissues. We first determined the expression patterns of the KRAS protein in human CRC tissues. We found that the KRAS protein expression levels were significantly higher in the CRC tissues than in the normal adjacent tissue (NAT) samples ( Fig. 1A and B). The mean expression levels of KRAS protein in the CRC and NAT samples are shown in Supplementary Fig. 1A. However, although the KRAS protein was consistently upregulated in the CRC samples, KRAS mRNA levels appeared to be irregular between the tumor and paired non-tumor samples (Fig. 1C). The mean expression levels of KRAS mRNA in the CRC and NAT samples are shown in Supplementary Fig. 1B. This disparity between the KRAS protein and mRNA expression in the CRC tissues suggests that a post-transcriptional mechanism is involved in the regulation of KRAS. However, it is noted that a very large difference was observed in expression of the KRAS in different samples (e.g., NAT#2 vs. NAT#4). Such inconsistency may be due to different age, gender and pathological stages of the CRC patients.
Prediction of conserved miR-16 binding sites within the 3′-UTR of KRAS. One important mode of post-transcriptional regulation is the repression of mRNA transcripts through miRNAs. Three algorithms (TargetScan 31 , miRanda 32 and PicTar 33 ) were used in combination to identify potential miRNAs that target KRAS. All three algorithms predicted miR-16 as a candidate miRNA that targets KRAS. The predicted interaction between miR-16 and KRAS mRNA is illustrated in Fig. 2A. One potential binding site in the 3′ -UTR of KRAS mRNA was observed. There was perfect base-pairing between the seed region (the core sequence that encompasses the first 2-8 bases of the mature miRNA) and the cognate target. The minimum free energy value of the hybridization between miR-16 and KRAS mRNA was − 24.2 kcal/mol, which is well within the range of genuine miRNA target pairs 34,35 . Furthermore, the miR-16 binding sequence in the KRAS 3′ -UTR was highly conserved across species. We then investigated whether the expression levels of miR-16 were inversely correlated with the levels of the KRAS protein in CRC tissues. We measured the miR-16 levels in the same paired CRC and NAT samples and found that miR-16 was significantly decreased in CRC samples (Fig. 2B), which is consistent with the notion that the levels of miRNAs are opposite to that of their targets. The mean expression levels of miR-16 in the CRC and NAT samples are shown in Supplementary Fig. 1C. The inverse correction between miR-16 and KRAS protein levels (Fig. 2C) and the disparity between the miR-16 and KRAS mRNA levels (   Fig. 3A. As anticipated, the overexpression of miR-16 dramatically reduced the KRAS protein levels in SW480 and HT29 cells, whereas the knockdown of miR-16 significantly increased the KRAS protein levels in Caco2 cells ( Fig. 3B and C). We also determined the level at which miR-16 regulates KRAS expression. By measuring the expression of KRAS mRNA after transfecting CRC cells with pre-miR-16 or anti-miR-16, we found that miR-16 did not significantly affect the mRNA levels of KRAS (Fig. 3D). These results demonstrate that miR-16 specifically regulates KRAS expression at the post-transcriptional level.
To determine whether the negative regulatory effect of miR-16 on KRAS expression was mediated by the binding of miR-16 to the predicted target sites in the 3′ -UTR of the KRAS mRNA, the KRAS 3′ -UTR containing the predicted miR-16 binding site was inserted the downstream of the firefly luciferase gene in a reporter plasmid. The resulting plasmid was transfected into SW480 cells along with either pre-miR-16 or pre-miR-control, or into Caco2 cells along with either anti-miR-16 or anti-miR-control. As expected, the luciferase activity was remarkably reduced in SW480 cells transfected with pre-miR-16 compared to cells transfected with pre-miR-control, whereas the luciferase activity was significantly increased in Caco2 cells transfected with anti-miR-16 compared to cells transfected with anti-miR-control (Fig. 3E). Furthermore, we introduced point mutations into the corresponding complementary sites in the 3′ -UTR of KRAS to eliminate the predicted miR-16 binding site. The mutated luciferase reporter was unaffected by the overexpression or knockdown of miR-16 (Fig. 3E). This result suggests that the putative miR-16 binding site of KRAS strongly contributes to the miRNA-mRNA interaction.  and 60 h after the transfection of SW480 cells with pre-miR-control plus a control plasmid, pre-miR-control plus a KRAS overexpression plasmid, pre-miR-16 plus a control plasmid, or pre-miR-16 plus a KRAS overexpression plasmid. (D and E) Transwell analysis of SW480 cells transfected with pre-miR-16 or pre-miR-control, or with pre-miR-control plus a control plasmid, pre-miR-control plus a KRAS overexpression plasmid, pre-miR-16 plus a control plasmid, or pre-miR-16 plus a KRAS overexpression plasmid. At the same time, Caco2 cells were transfected with anti-miR-16 or anti-miR-control and then subjected to Transwell analysis. D: representative image; E: quantitative analysis. (F and G) An apoptosis assay was performed 48 h after the transfection of SW480 cells with pre-miR-16 or pre-miR-control, or with pre-miR-control plus a control plasmid, pre-miR-control plus a KRAS overexpression plasmid, pre-miR-16 plus a control plasmid, or pre-miR-16 plus a KRAS overexpression plasmid. At the same time, Caco2 cells were transfected with anti-miR-16 or anti-miRcontrol and then subjected to apoptosis analysis. Cell apoptosis profiles were analyzed by flow cytometry. F: representative image; G: quantitative analysis. (mean ± S.D.; *p < 0.05; **p < 0.01; ***p < 0.001).
opposite effect on cell proliferation in Caco2 cells (Fig. 4A and B). Then, we examined whether the overexpression or knockdown of KRAS would impact the proliferation of SW480 cells. To knockdown KRAS, four siRNA sequences targeting different sites of KRAS cDNA were synthesized, and the sequence with the best interfering effect (si-KRAS#2) was selected. To overexpress KRAS, an expression plasmid designed to specifically express the full-length KRAS ORF without the miR-16-responsive 3′ -UTR was constructed. As expected, the cells transfected with KRAS siRNA showed a reduction of KRAS mRNA and protein levels, whereas the cells transfected with a KRAS overexpression plasmid showed increased KRAS mRNA and protein levels (Supplementary Fig. 2A-C). Consequently, cells transfected with KRAS siRNA proliferated at a significantly lower rate than the control siRNA-transfected cells (Supplementary Fig. 2D), whereas the cell proliferation rate was increased in the cells transfected with the KRAS overexpression plasmid compared to the cells transfected with an empty control plasmid ( Supplementary Fig. 2E). Thus, miR-16 and KRAS have opposite effects on cell proliferation. Subsequently, we investigated whether the overexpression of miR-16-resistant KRAS (KRAS ORF) was sufficient to rescue the suppression of KRAS by miR-16 and attenuated the anti-proliferative effect of miR-16 in CRC cells. Cells transfected with pre-miR-16 and the KRAS overexpression plasmid showed significantly higher proliferation rates than cells transfected with pre-miR-16 alone (Fig. 4C), suggesting that the overexpression of KRAS rescued the miR-16-mediated downregulation of the proliferation rates of SW480 cells. Taken together, the results indicate that miR-16 can inhibit cell proliferation by silencing KRAS.
Furthermore, we assessed the effect of miR-16 and KRAS on the invasion ability of SW480 and Caco2 cells. The invasion assays indicated that the invasion rate of SW480 cells transfected with pre-miR-16 was significantly decreased compared to control cells, whereas anti-miR-16 increased the invasion ability of Caco2 cells (Fig. 4D and E). Transfection of KRAS siRNA dramatically reduced the number of SW480 cells that pass through the Transwell membrane, whereas transfection of the KRAS overexpression plasmid increased the invasion rate ( Supplementary Fig. 2F and G). Thus, the inhibition of cell invasion by KRAS knockdown was similar to that elicited by miR-16 overexpression, further indicating that miR-16 and KRAS have opposite effects on cell invasion. Moreover, the overexpression of miR-16-resistant KRAS remarkably attenuated the anti-invasion effect of miR-16 when SW480 cells were co-transfected with pre-miR-16 and the KRAS overexpression plasmid (Fig. 4D and E). In summary, these results suggest that miR-16 may suppress cell invasion by silencing KRAS.
We finally investigated apoptosis in SW480 and Caco2 cells via flow cytometric analysis. In SW480 cells with enhanced miR-16 expression, the percentage of apoptotic cells was significantly higher than in the control cells, whereas anti-miR-16 had the opposite effect on cell apoptosis in Caco2 cells (Fig. 4F and G). In addition, KRAS interference remarkably increased the percentage of apoptotic cells, whereas KRAS overexpression decreased cell apoptosis ( Supplementary Fig. 2H and I). Moreover, when SW480 cells were co-transfected with pre-miR-16 and the KRAS overexpression plasmid, KRAS dramatically attenuated the miR-16-induced apoptosis effect ( Fig. 4E and F), suggesting that KRAS might reverse the pro-apoptosis effect of miR-16. The results indicate that miR-16 may modulate cell apoptosis by downregulating KRAS.

The effect of miR-16 and KRAS on the growth of CRC cells in vivo.
We next evaluated the effects of miR-16 and KRAS on the growth of CRC xenografts in mice. The use of viral constructs provides the possibility to rapidly produce high intracellular levels of mature miRNA via the endogenous miRNA processing pathway. To generate a viral expression construct, a 300-bp fragment containing the genomic sequences of miR-16 was obtained by PCR amplification of human genomic DNA and cloned into a lentiviral expressing vector. SW480 cells were infected with the miR-16 overexpression lentivirus. The expression of mature miR-16 was observed to be 8-10-fold higher after infection ( Supplementary Fig. 3A), and KRAS protein expression was consequently inhibited (Supplementary Fig. 3B and C). Subsequently, SW480 cells were infected with the miR-16-overexpressing lentivirus, or transfected with the KRAS overexpression plasmid, or co-transfected with the miR-16-overexpressing lentivirus and the KRAS overexpression plasmid, and then the untreated or transfected cells were subcutaneously injected into C57BL/6J nude mice. After one month of xenograft growth in vivo, the mice were sacrificed and the weights of the tumors were measured. A marked reduction in the size and weight of the tumors was observed in the miR-16-overexpressing group compared to the control group, whereas the size and weight of the tumors in the KRAS-overexpressing group was dramatically increased. Additionally, KRAS overexpression attenuated the suppressive effect of miR-16 on tumor growth ( Fig. 5A and B), suggesting that miR-16 may inhibit tumor growth by silencing KRAS. Next, total RNA and protein were isolated from the xenograft tumors and analyzed. Tumors from the miR-16-overexpressing group showed a significant increase in the expression of mature miR-16 compared to tumors from the control group (Fig. 5C). Likewise, KRAS mRNA levels were unchanged in tumors from the miR-16-overexpressing group but increased in tumors from the KRAS-overexpressing group (Fig. 5D). Moreover, tumors from the miR-16-overexpressing group displayed reduced KRAS protein levels compared to tumors from the control group, whereas tumors from the KRAS-overexpressing group showed elevated KRAS protein levels ( Fig. 5E and F). More importantly, tumors with both miR-16 and KRAS overexpression exhibited significantly higher levels of KRAS than tumors with miR-16 overexpression (Fig. 5E and F), suggesting that KRAS overexpression could rescue the KRAS suppression caused by miR-16. Additionally, H&E staining of xenograft tissues showed less cell mitosis in the group implanted with the miR-16 lentivirus, whereas more cell mitosis was observed in the KRAS overexpression group (Fig. 5G). Xenografts with both miR-16 and KRAS overexpression exhibited more cell mitosis than xenografts with miR-16 overexpression (Fig. 5G). Immunohistochemical staining also revealed the presence of lower levels of KRAS in the tumors from mice implanted with miR-16-overexpressing cells, whereas the tumors from the KRAS-overexpressing mice showed increased KRAS protein levels ( Fig. 5G and H). Finally, the proliferative activity of the tumor cells was assessed via immunocytochemistry with the mouse monoclonal antibody Ki-67. The cell proliferation rate, which was measured by the percentage of Ki-67-positive tumor cells, was decreased in the group implanted with the miR-16 lentivirus and increased in the group treated with the KRAS plasmid ( Fig. 5G and I). Likewise, KRAS overexpression attenuated the anti-proliferative effect caused by miR-16 Scientific RepoRts | 6:37459 | DOI: 10.1038/srep37459 overexpression ( Fig. 5G and I). These results were consistent with the findings of the in vitro assays, which firmly validated the tumor-suppressive role of miR-16 in CRC tumorigenesis through the targeting of KRAS.

Discussion
CRC is one of the most common malignant human cancers worldwide. The development and progression of CRC is a complicated process that involves the dysregulation of a variety of genes that are essential for cellular processes 36 . Three members of the GTPases of the RAS family, KRAS, HRAS and NRAS, are well known for their ability to cause neoplasia. Among the RAS family, KRAS is one of the most prominent oncogenes due to its ability to transform human cells into malignant tumor cells, particularly when harboring an activating mutation in codon 12 or 13 37 . KRAS mutations frequently occur in many types of human cancers, for example in 70-90% of pancreatic cancer, 30-60% of colon cancer and 15-50% of lung cancer patients 12 . Moreover, activating oncogenic KRAS mutations are often associated with resistance to chemotherapy and targeted therapies in CRC 38,39 . Due to the poor prognosis for cancer patients with mutated KRAS, much effort has been spent on developing specific therapies for targeting oncogenic KRAS. However, technical limitations make it difficult to specifically inhibit KRAS in vivo. Apart from specific RNA interference methods, currently, there are no small molecules available that can specifically target KRAS. In this study, we showed that silencing KRAS expression could inhibit cell proliferation and invasion and induce apoptosis in CRC cells, while overexpressing KRAS had the opposite effects on CRC cells, validating the role of KRAS as a crucial oncogene during CRC tumorigenesis. More importantly, this study identified miR-16 as a novel link between the KRAS regulatory pathway and CRC and pointed the important role of miR-16 as a tumor suppressor in CRC through the inhibition of KRAS translation. Considering that miR-16 is an upstream regulator of KRAS, it is possible to restore miR-16 expression to inhibit KRAS expression in vivo. Currently, the overexpression of miRNAs can be silenced using antagomirs, and the re-expression of miRNAs that are lost in cancers can be achieved by the overexpression of miRNA mimics. Greater research emphasis is needed to characterize the feasibility of targeting miR-16 in CRC therapy and develop simplified and cost-effective manipulation methods. miR-16 was first found to be associated with chronic lymphocytic leukemia 40,41 . Subsequent studies demonstrated that miR-16 is downregulated and serves as a tumor suppressor in a variety of cancer types. For example, Rivas et al. showed that downregulation of miR-16 via progestin-mediated oncogenic signaling contributed to breast cancer development 42 . Ma et al. reported that miR-16 could inhibit proliferation and induce the apoptosis of CRC cells by regulating the P53/survivin signaling pathway 43 . Ke et al. found that miR-16 was downregulated in non-small cell lung cancer cells and had an inverse correlation with the protein expression of hepatoma-derived growth factor 44 . Liang et al. reported that miR-16 showed an inverse correlation with the protein expression of FEAT in breast cancer, lung cancer and hepatocellular cancer tissues and that the overexpression of miR-16 promoted the apoptosis of cancer cells by targeting FEAT 45 . In this study, we identified a novel regulatory network that employs miR-16 and KRAS to regulate cell proliferation, invasion and apoptosis in CRC cells. We also provided evidence that the restoration of KRAS expression could partly rescue miR-16-suppressed cell proliferation and invasion and miR-16-induced apoptosis, suggesting that the targeting of KRAS is an important mechanism by which miR-16 exerts its tumor-suppressive function. Because the detailed roles that miR-16 plays in the initiation and progression of CRC remain not fully understood, the results of this study may explain, at least in part, why the downregulation of miR-16 during CRC carcinogenesis can promote cancer progression.

Conclusions
In this study, we found that the expression levels of KRAS were significantly higher in CRC clinical tissues than in the non-tumor adjacent tissues. In addition, we showed for the first time that KRAS is a direct target of miR-16. We also provided evidence that miR-16 could inhibit the proliferation and invasion and induce the apoptosis of CRC cells by silencing KRAS. Furthermore, experiments using tumor xenografts in mice validated the tumor-suppressive role of miR-16 in CRC tumorigenesis through the targeting of KRAS. Taken together, this study highlights an important role for miR-16 in the regulation of KRAS in CRC cells and may open new avenues for future CRC therapy.