- Open Access
MicroRNA-181a promotes tumor growth and liver metastasis in colorectal cancer by targeting the tumor suppressor WIF-1
Molecular Cancer volume 13, Article number: 86 (2014)
Given the emerging role of microRNA in tumor disease progression, we investigated the association between microRNA expression, liver metastasis and prognosis of colorectal cancer.
Colorectal cancer tissues from patients with or without liver metastases were profiled to identify differentially expressed microRNA. Expression profile was further assessed using quantitative reverse transcription PCR and in situ hybridization. Correlation between miR-181a expression, the most differentially expressed microRNA, between patients with and without liver metastasis, and its downstream target genes were investigated using qRT-PCR. Luciferase reporter assay was conducted to establish functional association between miR-181a and its target genes. Manipulation of miR-181a expression and its consequences in tumor growth and metastasis were demonstrated in various in vitro and in vivo models.
miR-181a was revealed being the most elevated in CRC with liver metastases. miR-181a expression correlated with advanced stage, distant metastasis, and served as an independent prognostic factor of poor overall survival. Stable transfection of CRC cell lines with miR-181a promoted cell motility and invasion, as well as tumor growth and liver metastasis,while silencing its expression resulted in reduced migration and invasion. Additionally, we identified WIF-1 as direct and functional targets of miR-181a. Ectopic expression of miR-181a suppressed the epithelial markers E-cadherin and β-catenin, while enhanced the mesenchymal markers vimentin.
Our data demonstrate that miR-181a expression is associated with CRC liver metastasis and survival. miR-181a has strong tumor-promoting effects through inhibiting the expression of WIF-1, and its potential role in promoting epithelial-mesenchymal transition.
Colorectal cancer (CRC) is the third most prevalent cancer worldwide . Approximately 50% of CRC patients die from distant metastases, in particular liver metastases. Most CRC patients with liver metastasis are not candidates for surgical treatment, with 5-year survival rate below 10% . New diagnostic means to detect early metastasis and to predict metastatic risk in CRC are in urgent demand.
microRNAs (miRNAs) are small, non-coding RNA that silence specific target genes in mammalian cells by repressing translation . miRNAs play important roles in tumor invasion and metastasis . Depending on its target genes, an miRNA can function either as an oncogene or a tumor suppressive gene . Dysfunction of miRNA is associated with CRC tumorigenesis and progression. Furthermore, aberrant expression of specific miRNA may be used as potential prognostic and predictive markers in CRC [6–8]. However, more extensive research is required to elucidate molecular mechanisms of miRNAs that mediate CRC liver metastasis and to identify those miRNAs that may be employed as novel prognostic predictors.
In the present study, microarray-based strategy was applied to identify differentially expressed miRNA in CRC by comparing miRNA profiles between colorectal cancer tissues from patients with and without liver metastases. miR-181a was most significantly up-regulated in cancerous tissues with liver metastasis compared to without liver metastasis. Its tumor promoting function and mechanism were further characterized.
miRNAs are differentially expressed in CRC tumor tissues from patients with or without liver metastasis
To investigate the role of miRNA in CRC liver metastasis, miRNA expression was profiled in CRC tumor tissue from patients with (n = 5) or without (n = 5) liver metastases. Among the 743 miRNA probes on Mammalian miRNA arrayV2.0, 214 miRNAs were detected in CRC tissues. hsa-miR-181a was detected to be most elevated with > 2-fold between CRC with and without liver metastasis. (Additional file 1: Table S1).
The microarray results were validated with quantitative reverse transcription PCR (qRT-PCR) by confirming the expression of hsa-miR-181a in the same 10 human CRC tissues. All further analyses in this study focused on miR-181a.
Deregulated expression of miR-181a is associated with liver metastasis in CRC patients
To verify the microarray results of differentiated miRNA expression, we performed qRT-PCR analyses in a cohort of 137 paired CRC and normal tissues (Validation Cohort 1).
qRT-PCR analysis revealed a significant increase in miR-181a expression in tumor tissue (n = 137) compared with matched adjacent normal tissues (n = 137, fold change: 2.46, p = 0.0002). We then analyzed the relationship between miR-181a expression and liver metastasis. The relative expression ratio of miR-181a to U6 in patients with synchronous liver metastasis (n = 63) is significantly higher than those without metastasis (n = 60, fold change: 2.56, p = 0.0003). In addition, its expression in patients with metachronous liver metastasis (n = 14) is remarkably higher than those without metastasis (fold change: 1.93, p = 0.0114). These data demonstrate that miR-181a expression is elevated in CRC tumor tissues compared to surrounding normal tissue, and this elevated expression level of miR-181a correlates with liver metastasis in CRC patients (Figure 1).
Given that miR-181a is part of a family that includes miR-181a-d, we also checked the expression of three other miR-181 family members, miR-181b, miR-181c and miR-181d, to see whether they were also associated with metastases. Our results show that miR-181b-d expressions are not associated with CRC liver metastasis (Figure 1G-I). These results suggest that the association of miR-181a with CRC metastases is specific. The miRBase accession number and the mature sequence of the microRNAs studied are available in Additional file 1: Table S2.
Elevated miR-181a expression correlates with advanced clinical stage, distant metastasis, and serves as an independent prognosis factor of poor survival
The clinical significance of miR-181a up regulation in CRC was further investigated in Validation Cohort 1. The expression of miR-181a increased remarkably with the advanced CRC TNM stages. Elevated miR-181a expression was significantly correlated with local tumor invasion (T stage), lymph node metastases or liver metastases, but not with other clinical-pathological parameters (Table 1, Figure 1C-E). Kaplan-Meier analysis showed that high miR-181a expression (above the median value in this 137 patient cohort) was associated with poor overall survival in CRC (Figure 1F). These results suggested that the up regulation of miR-181a may play an important role in the development and progression of CRC.
Furthermore, multivariate analysis was performed to determine the prognostic value of miR-181a expression using the Cox proportional hazard model. The risk variables examined included miR-181a expression level, as well as factors known to significantly affect the outcome of CRC such as age and gender of patients, differentiation, lymph node metastasis, surgical-pathological staging, histological type and venous invasion. In the univariate analysis, high miR-181a expression level in tumors (HR 2.22; 95% CI 1.307 to 3.77; p = 0.003), TNM staging (HR 11.057; 95% CI 3.995 to 30.605; p < 0.0001), venous invasion (HR 2.462; 95% CI 1.476 to 4.108; p = 0.001), and lymph node metastasis (HR 4.27; 95% CI 2.156 to 8.454; p < 0.0001) were significantly associated with survival; while age, gender, differentiation and histological type were not (Table 2). In the final multivariate Cox regression model, high miR-181a expression in tumors was associated with a poor survival prognosis (HR 1.871; 95% CI 1.078 to 3.246; p = 0.026) independent of other clinical covariates (Table 2).
miR-181a localizes in colorectal epithelial cells and serves as an independent prognostic marker of CRC overall survival
To further elucidate the functional role of miR-181a in CRC, we utilized in situ hybridization (ISH) to study the expression pattern of miR-181a in tumor and corresponding non-tumor tissue. A high level expression of miR-181a was detected in epithelial cells with elevated expression in tumor tissue compared with normal tissue. In cancerous tissue, miR-181a expression was detected only in the cytoplasm of tumor cells but not in tumor stroma, consistent with its role in tumor cells during colorectal carcinogenesis. Low levels of miR-181a expression were also seen in normal mucosa (Figure 2). To assess whether miR-181a expression is correlated with miR-181a ISH score, we assayed miR-181a expression by qPCR in 30 fresh CRC samples frozen in liquid nitrogen that were matched to the samples used for ISH. Our results showed that miR-181a expression was correlated with miR-181a ISH score (Figure 2H).
Expression of miR-181a was characterized in Validation Cohort 2, micro tissue arrays consisting of tumor samples from 294 CRC patients with different liver metastasis status, e.g., none, synchronous or metachronous liver metastasis to further validate the association between miR-181a expression and liver metastasis. Consistent with qRT-PCR results derived from Validation Cohort 1, miR-181a was significant less frequently expressed in PNLM (16/97, 16.5%) and normal tissue (14/72, 19.4%) than in PSLM (63/160, 39.4%) and PMLM (22/37, 59.5%) (p < 0.001)
Moreover, elevated miR-181a expression in tumor tissue detected by ISH was significantly associated with poor survival (p < 0.001) (Figure 2G). In univariate analysis of Validation Cohort 2, high expression of miR-181a in tumor tissue (HR 1.514; 95% CI 1.223 to 1.874; p < 0.001), TNM staging (HR 6.509; 95% CI 4.02 to 10.54; p < 0.001), venous invasion (HR 2.375; 95% CI 1.758 to 3.21; p < 0.001), differentiation (HR 0.438; 95% CI 0.224 to 0.858; p = 0.016) and lymphnode metastasis (HR 2.916; 95% CI 2.062 to 4.125; p < 0.001) were significantly associated with survival; while age, gender, and histological type were not. In the final multivariate Cox regression model, high miR-181a expression in tumor tissue was associated with a poor survival prognosis (HR 1.382; 95% CI 1.11 to 1.72; p = 0.004) independent of other clinical covariates consistent with results derived from Validation Cohort 1 (Table 2).
miR-181a promotes CRC cell growth, invasion and liver metastasis
To explore the potential mechanism of miR-181a in promoting CRC progression, HT29 CRC cells were infected with pri-miR-181a lentivirus to determine its effects on cell growth in vitro and in vivo. Stable expression of miR-181a was confirmed using quantitative PCR (Figure 3A). Elevated miR-181a expression results in stimulated cell growth in culture (Figure 3B). Consistent with in-vitro results, the tumor growth curves showed that the growth of lentiviral miR-181a HT29 tumors was markedly increased compared with tumors derived from parental cell line HT-29, and lentiviral empty vector-infected cells (HT29-NC). Elevated miR-181a expression in HT29 CRC cells stimulated tumor growth in NOD-SCID mice resulting in increased tumor weight 1.66 ± 0.08 g (HT29-181a), 0.72 ± 0.13 (HT29-NC), 0.81 ± 0.10 g (HT29) (p < 0.05) (Figure 3C).
The effect of miR-181a on tumor invasion and metastasis was assessed. Cell wound healing assay showed that elevated miR-181a expression significantly enhanced tumor cell mobility in HT29-181a cells compared with HT29-NC or the parental HT 29 cells (p = 0.004) (Figure 3D). Furthermore, the ability of cell migration and cell invasion through matrigel were significantly enhanced by overexpressing miR-181a in HT29 cells (HT29-181a) compared with controls (HT29-NC and HT29) as demonstrated in Boyden Chamber assays (migration 2.74-fold, p = 0.03, invasion 34.7-fold, p = 0.049 respectively) (Figure 3E). The migration and invasion through matrigel of lenti–miR-181a-infected SW480 cells were dramatically enhance by 3.75-fold (p = 0.0008) and 4.67-fold (p = 0.0003) respectively, compared with controls (Figure 3G). A wound-healing assay revealed that the spreading of lenti–miR-181a-infected SW480 cells into the wound was much faster than the control cells (p = 0.0083) (Figure 3F). These results indicate that miR-181a stimulates colorectal cancer cell spreading, migration and invasion, suggesting its potential role in cancer metastasis.
The promoting effect of miR-181a on liver metastasis was demonstrated applying an experimental metastasis model by injecting HT29 overexpressing miR-181a cells into mouse spleens. HT29 cells with elevated miR-181a expression resulted in increased number of liver metastatic nodules as compared to that induced by HT29 cells transfected with control vector or the parental HT29 cells with average liver surface metastatic nodules of 6 ± 1.6; 2.4 ± 1.3; 1.2 ± 1.3 for HT29-181a, HT29-NC, and HT29 respectively (p < 0.05) (Figure 3H).
These combined in vitro and in vivo results illustrate the role of miR-181a in promoting tumor metastasis consistent with its clinical association with liver metastases in CRC patients.
miR-181a targets the 3′-UTR of tumor suppressor gene WIF-1
To elucidate the biological mechanisms underlying the role of miR-181a in promoting tumor cell growth and metastasis, we investigated the potential targets of miR-181a. Target prediction programs, miRanda and TargetScan, were applied to identify WIF-1 as putative miR-181a target. The 3′-UTR of WIF-1 mRNA contains a complementary site for the seed region of miR-181a (Figure 4A). To verify this finding, WIF-1 3′-UTRs and its mutant containing the putative miR-181a binding sites were cloned downstream of the luciferase open reading frame. These reporter constructs were co-transfected into HEK293T cells with either miR-SC or miR-181a mimics. Increased expression of miR-181a upon infection of miR-181a mimics, significantly suppressed luciferase expression derived from reporter constructs containing wild type WIF-1 3′-UTRs with inhibition rates 40% (p < 0.05) comparing to cells co-transfected with miR-SC. This suppressive effect was abolished when mutated 3′-UTR of WIF-1 mRNAs, in which the binding sites for miR-181a were inactivated by site-directed mutagenesis, were co-infected with miR-181a (Figure 4B). These results support WIF-1 as putative target of miR-181a.
Functional regulation of WIF-1 expression by miR-181a was further analyzed by modulating miR-181a levels via overexpression or knockdown in three CRC cell lines, HT29, SW480 and SW620. WIF-1 mRNA levels were substantially suppressed in HT29 overexpressing miR-181a and SW480 overexpressing miR-181a cells as compared with that in control cells (Figure 4C, E). Meanwhile, the protein levels of WIF-1 were also suppressed after ectopic overexpression of miR-181a in HT29 and SW480 cell lines (Figure 4F, G). On the other hand, knock down of miR-181a via RNA interference in HT29 andSW620 cells resulted in increased mRNA and protein levels of WIF-1 (Figure 4C-G).
Collectively, these data support the bioinformatics prediction of WIF-1 as direct targets of miR-181a and established a functional association.
Knockdown of miR-181a suppresses tumour growth and metastasis
To confirm further the tumour-promoting function of miR-181a, we infected the highly metastatic cell line SW620 with miR-181a-RNAi lentivirus and measured its effects on cell spreading, migration and invasion in vitro. Consistent with the results of over-expressing miR-181a in HT-29 cells, the inhibition of miR-181a in SW620 cells remarkably decreased the cell motility (p = 0.02) and cell invasion through matrigel (p = 0.03), as demonstrated by Boyden chamber assays (Figure 5A). A wound-healing assay revealed that the spreading of lenti–miR-181a-RNAi infected SW620 cells into the wound was much slower than the control cells (Figure 5B).
Knockdown of WIF-1 enhances tumour growth and metastasis
To explore the biological functions of WIF-1 in CRC cells, we determined whether inhibition of WIF-1 promoted CRC cell migration and invasion. In HT-29 cells, in vitro knockdown of WIF-1 promoted cell migration (p < 0.0001) and cell invasion (p = 0.0018) (Figure 5C). To further test this, we transfected SW620-181a-RNAi cells with siRNA for WIF-1 mRNA and found that the effect of miR-181a RNAi was partially attenuated by siRNA for WIF-1 mRNA (Figure 5D-E). These data confirmed that the prometastasis effect of miR-181a was mediated by inhibiting target genes WIF-1.
miR-181a levels are inversely correlated with mRNA and protein expression of WIF-1 in CRC tissues
The regulation of WIF-1 mRNA expression by miR-181a was further characterized in Validation Cohort 1 of 137 primary CRC tissues using qRT-PCR. A significant inverse correlation between the levels of miR-181a and mRNA expression of WIF-1 was confirmed with Pearson correlation = −0.5291; and p < 0.0001 (Figure 6A).
The protein expression of WIF-1 was also examined in paired CRC and normal tissues by Western Blot analysis. A significant decrease in WIF-1 expression was seen in tumor tissues compared with the matched adjacent normal tissue (Figure 6B), while their miR-181a expressions were shown an inverse trend (Figure 6C). The same inverse correlation was further demonstrated when protein expression of WIF-1 was examined in Validation Cohort 2 of TMA. High levels of miR-181a in tumor cells correlated with low levels of WIF-1 protein expression. Conversely, low levels of miR-181a in adjacent normal tissue correlated with high levels of WIF-1 protein expression (Figure 6D). We also investigated the prognostic significance of WIF-1 in CRC Cohort 2, Kaplan-Meier analysis was performed illustrating that low WIF-1 expression was associated with poorer overall survival in the CRC cohort 2(p = 0.004, Figure 6E).
Overexpression of miR-181a induces epithelial-mesenchymal transition (EMT) in CRC
The cell morphology of the stable SW480 hsa-miR-181a–overexpressing cell line was significantly altered; the cells had long and thin processes, cell-to-cell contact was reduced. The cell morphology of the miR-181a-knockdown SW620 was also altered, and the cell-to-cell contaction was increased (Figure 7A).
When miR-181a was overexpressed in SW480 cells, two epithelial markers, E-cadherin and β-catenin, were suppressed; while the mesenchymal marker, vimentin, was induced. On the other hand, knockdown of miR-181a in SW620 cells induced E-cadherin and β-catenin expression and suppressed vimentin expression (Figure 7B). Consistent results were also obtained when these CRC cells expressing different levels of miR-181a were examined with immunofluorescence staining (Figure 7C). We also assessed the epithelial and mesenchymal markers in HT29 overexpressing miR-181a cells. The results are consistent with results derived from SW480 cells (Figure 7B). We examined the effect of WIF1 on EMT in HT29 cell lines. The expression levels of E-cadherin and β-catenin were decreased; while vimentin was increased in HT29-WIF-1 siRNA cells (Figure 7B).
Taken together, these results suggest that elevated miR-181a levels in CRC may be involved in inducing EMT and therefore further promoting invasion and metastasis.
Recent studies have demonstrated that miRNA play important roles in tumor invasion and metastasis. In this study, miRNA profiles between five colorectal cancer tissues with synchronous liver metastases and five colorectal cancer tissues without liver metastases were compared. miR-181a demonstrated overexpression in colorectal cancer with liver metastasis.
Following microarray, we validated the deregulated miRNAs in two cohorts. We observed that up regulation of miR-181a is a frequent event in CRC tissues. High-level expression of miR-181a was significantly associated with colorectal cancer liver metastasis and especially with metachronous liver metastasis. Due to the largest absolute fold change, miR-181a was chosen for further study. The up regulation miR-181a was significantly associated with advanced clinical stage, distant metastasis and poorer overall survival in CRC.
Our results are similar with Nishimura’s report . It is reported that high miR-181a expression group had a significantly poorer prognosis. High miR-181a expression was an independent significant prognostic factor for CRC. However, in Nishimura’s report, they did not find miR-181a to be differentially expressed between the normal colon and colorectal cancer tissue, nor did they see any associations between miR-181a and any clinicopathological parameters investigated. In the present study, qRT-PCR and ISH analysis revealed a significant increase in miR-181a expression in tumor tissue compared with matched adjacent normal tissues. We think the different result is associated with the different normal tissues origin and number in two studies. We think the difference of sample cohort’scomposition, especially the clinical stage distribution and the proportion of liver metastasis samples results in different analysis in the correlations between miR-181a expression and clinicopathological parameters.
The role of miR-181a in cancer development is complicated. For instance, some studies demonstrated that it is up-regulated in pancreatic cancer and breast cancer [10, 11]. However, different groups showed that miR-181a was down-regulated in gliomas  and aggressive CLL.  Several studies have shown that miR-181 family members were up-regulated in hepatocellular cancer stem cells [14, 15].
In this study, the effects of miR-181a on CRC development and progression were also investigated by both in-vitro and in-vivo assays. Overexpression of miR-181a could effectively promote CRC cell growth rate, migration and invasion in vitro and tumor growth and liver metastasis in vivo, whereas knockdown of miR-181a expression reduced CRC cell migration and invasion.
These findings suggest that miR-181a plays a critical role in the invasive and/or metastatic potential of CRC.
Each miRNA can potentially downregulate many target genes through binding their 3′-UTR. One important task is to identify the downstream target genes regulated by the dysregulated miRNA. In the present study, quantitative PCR and western blot analysis demonstrated that miR-181a could decrease the expression of WIF-1 in both mRNA and protein levels. In addition, inhibiting miR-181a expression could increase WIF-1 expression. Luciferase assay showed that miR-181a could interact with the 3′-UTR of WIF-1. We also found that knockdown of WIF-1 gene expression by siRNA reversed miR-181a-RNAi-mediated suppression of tumor cell migration. Moreover, the expression of WIF-1 was negatively correlated with miR-181a in clinical CRC specimens. These observations indicate that miR-181a mechanistically acts via the regulation of WIF-1. Our results also showed that low-level expression of WIF-1 is positively correlated with tumor metastasis and/or poor prognosis. This is consistent with our finding that up regulation of miR-181a is associated with CRC liver metastasis and poorer overall survival.
WIF-1, Wnt inhibitory factor-1, is a secreted antagonist that can bind to Wnt proteins directly and inhibit Wnt signaling pathway. It has been reported that WIF-1 expression is down regulated in several human cancer including prostate, bladder, breast, nonsmall-cell lung carcinomas and colorectal cancer [16–18]. Several groups have reported that restoring WIF-1 expression in cancer cells could inhibit cancer cell growth [19, 20]. Our results, along with the roles of WIF-1 in tumor progression, indicate the involvement of the miR-181a-WIF-1 regulation chain in CRC progression.
It is reported that miR-181 directly targets nemo-like kinase (NLK), an inhibitor of wnt/β-catenin signaling in Hepatocellular cancer . We speculate that miR-181a could affect colorectal cancer cells EMT. EMT appears to be a key event in tumor invasion and metastasis [21, 22]. Thus, epithelial cells lose their epithelial adherence and tight junction proteins, lose their polarity and cell-cell contact, and undergo remarkable remodeling of the cytoskeleton, all of which facilitates cell motility and invasion. In present study, we found that the overexpression of miR-181a induces EMT, as shown by the decreased expression of the epithelial markers E-cadherin and β-catenin, and enhanced expression of the mesenchymal markers vimentin.
These findings suggest that miR-181a plays a critical role in regulating epithelial-emesenchymal cell transition and, ultimately, promotes the invasive and/or metastatic potential of CRC, probably by its direct target on WIF-1 gene. We further assessed whether the miR-181a induced EMT is due to the direct regulation of WIF1. We found that inhibition of WIF1 could lead to EMT in HT29 cell line, which suggesting that the miR-181a induced EMT is due to the direct regulation of WIF1.
In summary, we demonstrated in this study that the higher expression of miR-181a in CRC specimens is associated with liver metastasis and poor survival, and in particular, the higher expression of miR-181a is associated with metachronous liver metastasis. In this study, we also investigated the potential role of miR-181a in CRC liver metastasis and its underlying mechanisms. Our data suggest that up regulation of miR-181a plays an important role in CRC cell metastasis, and that the tumor-promoting function of miR-181a is through repressing its downstream target gene WIF-1. miR-181a could be employed as a new prognostic marker and/or as an effective therapeutic target for CRC.
Materials and methods
Cell line and experimental animals
Human CRC cell lines HT29, SW480 and SW620 were purchased from American Type Culture Collection and were maintained in Dulbecco’s Modified Eagle medium (Gibco, Carlsbad, CA) with 10% fetal bovine serum, 100 U/ml penicillin sodium and 100 mg/ml streptomycin sulfate in humidified 5% CO2 at 37°C. Male NOD-SCID mice (4–6 weeks old) were purchased from Beijing HFK Bio-technology Co. Ltd. (Beijing, China). Mice were maintained in a pathogen-free facility and used in accordance with the institutional guidelines for animal care.
Patients and samples
CRC samples were collected from surgical resection at Peking University Cancer Hospital & Institute (Beijing, China). All samples were immediately frozen in liquid nitrogen and stored at −80°C or fixed in 10% formalin for paraffin embedding. Samples collection and usage in the present study were approved by the Ethics Review Committees of Peking University Cancer Hospital & Institute.
The Testing Cohort consisted of CRC tumor specimens from patients with synchronous liver metastases (n = 5) or without (n = 5) liver metastases (follow up > 6 years). Samples were analyzed using microarray to generate differential miRNA profiling between CRC tumor tissues with and without liver metastases. Validation Cohort 1 consisting of fresh CRC and surrounding non-tumor tissue were obtained from 137 patients who underwent surgical resections from 2001 to 2006. These samples were analyzed by quantitative reverse transcription PCR (qRT-PCR). Validation Cohort 2, independent from Validation Cohort 1, consisting of CRC and surrounding non-tumor tissues were obtained from 294 patients who underwent surgical resections from 1999 to 2006. These samples were paraffin-embedded and micro-dissected to generate tissue microarrays (TMAs) for in situ hybridization (ISH) and immunohistochemstry. All samples were collected from patients prior to receiving any preoperative chemotherapy or radiotherapy, and were verified to contain at least 80% tumor cells. Cases with familial adenomatous polyposis CRC were excluded from the study. Final cut off date of follow-up was January 30, 2012. All samples have follow-up data. A summary of the clinical characteristics of these patients is shown in Additional file 1: Table S3 and S4.
RNA isolation and miRNA profiling
Total RNA was extracted from 5 CRC tissues with and 5 without liver metastases using TRIZOL reagent (Invitrogen) according to the manufacturer’s instructions. miRNA microarray profiling was performed at CapitalBio Corporation (Beijing, China) following standard procedures detailed on http://www.capitalbio.com. Fluorescein-labeled miRNA were used for hybridization on Mammalian miRNA arrayV2.0 (CapitalBio, Beijing, China) containing 743 probes in triplicate, corresponding to 576 human miRNAs (including 122 predicted miRNAs). Hybridization signals were detected and arrays were scanned with a LuxScanTM 10 K/A laser confocal scanner (CapitalBio, Beijing, China). The raw pixel intensities were extracted using the LuxScan3.0 software (CapitalBio, Beijing, China). The data were normalized between slides using the quantile normalization method proposed by Bolstad et al. . The differentially expressed genes, classified as those with Fold changes above 2, were selected using the SAM software, version 2.1.
Quantitative real-time RT-PCR (qRT-PCR)
Total RNA was extracted using miRNeasy Mini Kit (Qiagen, Hilden, Germany). For miRNA quantification, the miRNAs levels were assayed using TaqMan MicroRNA assays, according to the manufacturer’s protocol (Applied Biosystems, Carlsbad, CA). The u6 small nuclear B noncoding RNA (RNU6B, Applied Biosystems, Carlsbad, CA) level was used as an internal normalization control.
For mRNA detection, total RNA (2 μg) was used in 20 μl of reverse transcriptase reaction to synthesize cDNA by using Moloney murine leukemia virus reverse transcriptase (M-MLV RT) for RT-PCR (Invitrogen, Carlsbad, CA), according to the manufacturer’s protocol. Real-time PCR was performed using the above RT products with the SYBR Green PCR Master Mix (Toyobo Co. Ltd., Osaka, Japan) on an ABI7500 PCR machine. PCR cycling conditions were: 95°C for 1 minute, 40 cycles of 95°C for 15 seconds, 60°C for 15 seconds and 72°C for 45 seconds. Data are presented as relative quantification (RQ) to U6 or GAPDH based on calculations of 2-ΔCt where ΔCt = Ct(Target)-Ct(Reference), or as ΔCt. Fold changes were calculated by the 2-ΔΔCt. All primers are listed in Additional file 1: Table S5.
Paraffin-embedded samples were microdissected and tissue microarrays (TMAs) were constructed in quadruplicates using a previously described method . Seven TMAs were constructed to include 160 Primary CRC tumors with Synchronous Liver Metastasis (PSLM), 97 Primary CRC tumors with Non-Liver Metastasis (PNLM), 37 Primary CRC tumors with metachronous liver metastasis (PMLM), and 72 paired adjacent normal colorectal mucosa and 32 matched liver metastasis tissues.
In situ hybridization analysis
Digoxigenin (DIG)-labeled mercury locked nucleic acid probes for human miR-181a, U6 (positive control) and scramble DNA (negative control) were purchased from TaKaRa (Dalian, China). The in situ hybridization (ISH) was performed with a modified version of the protocol for formalin-fixed paraffin-embedded tissue by Wigard Kloosterman (http://www.exiqon.com/ls/Documents/Scientific/FFPE%20in%20situ%20hybridization.pdf) on human colorectal cancer tissues. The sequence and hybridization temperature of probes were showed in Additional file 1: Table S5.
The slides were then scored by two pathologists independently as negative (−), weak or focally positive (1+), or strongly positive (2+) . Both pathologists were blinded to the patient’s clinical outcome.
Oligonucleotide, lentiviral vector construction and cell infection
miRNA mimics, the miRNA inhibitor and negative control miRNA oligonucleotides of has-miR-181a were from RiboBio Co. Ltd (Guangzhou, China). The small interfering RNAs (siRNAs) targeting WIF-1 were synthesized by RiboBio Co. Ltd. An unrelated sequence was used as a negative control (provided by RiboBio). The sequence was WIF1 siRNA, 5′ GAGUACUCAUAGGAUUUGA dTdT -3′ (sense). Stable transfectants overexpressing or knockdown miR-181a were generated by lentiviral transduction using a pGCsil-GFP vector (GeneChem Co., Ltd, Shanghai,China). A lentiviral vector expressing green fluorescent protein alone (LV-GFP) was used as a control. Transfection of oligonucleotides or lentivirus construction was performed using the Lipofectamine 2000 reagent (Invitrogen) according to the manufacturer’s instructions. HT-29, SW480 and SW620 cells were infected with recombinant lentivirus-transducing units plus 8 μg/ml Polybrene (Sigma, St Louis, Missouri, USA).
In vitro cell growth, spreading, motility and invasion assays
Cell growth was assessed using CCK-8 (Dojindo, Tokyo, Japan) according to the manufacturer’s instructions.
Cell spreading was evaluated by a wound healing assay. The wound healing assay was performed by the CellPlayer™ 96-Well Kinetic Cell Migration assay (The Essen BioScience IncuCyteTM Live-Cell Imaging System). Cells were seeded into 96-well plate and cultured until attached. The woundMaker™ and wounding procedure were used to create precise and reproducible wounds in all wells of the 96-well plate. The plate was placed into the IncuCyte™ ZOOM to scan every hour for migration assays. The data was analyzed by wound confluence, Wound confluence (%) represents the fractional area of the wound that is occupied by cells.
Cell motility or invasion was measured using a Boyden chamber assay with or without Matrigel. Photographs of three randomly selected fields of the fixed cells were taken and cells were counted. Experiments were repeated three times independently. These assays were carried out as previously described .
In vivo tumor growth and liver metastasis assay
For the tumorigenecity study, HT29 cells, has-miR-181a-expressing cells (HT29-181a+) and empty vector cells (HT29-NC) were injected subcutaneously into the back of NOD-SCID mice (8 mice/group). Over a 4-week period, tumor formation in mice was observed by measuring the tumor volume calculated by the formula V = 0.5xLxW2. Tumors were then excised and weighed.
HT29 cells, has-miR-181a-expressing cells (HT29-181a+) and empty vector cells (HT29-NC) were assayed for metastasis in 5- to 6-week-old NOD-SCID mice using the splenic injection model as described previously . Cells (1 × 106) were injected into the spleen subcapsular through an abdominal incision under sterile conditions. After injection, the spleen was returned gently back to abdominal cavity and homoeostasis was assured. The area was thoroughly irrigated with warm sterile water and the abdominal cavity was closed in appropriate layers. Ten weeks after surgery, mice were euthanized and the livers were obtained to determine the liver weight and hepatic metastatic foci.
All animal experiments were reviewed and approved by the Ethics Review Committee at the Peking University School of Oncology.
Luciferase reporter assay
The 3′-UTR of WIF-1 was amplified via polymerase chain reaction (PCR) from human genomic DNA constructed and cloned into pmiR-RB-REPORT™Luciferase, downstream of the firefly luciferase gene to form pMIR-3′-UTR WIF-1 and pMIR-3′-UTR WIF-1-Mut constructs (primers are listed in Additional file 1: Table S3).
For the luciferase assay, HEK 293 T cells were seeded in 96-well plates, pMIR-3′-UTR constructs were co-transfected with miR-181a mimics or miR scramble (miR-SC) (50 nM, synthesized by Ribobio Co. Guangzhou, China) using the Lipofectamine2000 method. The activities of firefly luciferase and Renilla luciferase were measured at 24 h post-transfection using the dual-luciferase reporter assay kit (Promega) following the manufacturer’s protocol. Firefly luciferase activity was normalized to that of the Renilla luciferase for each sample. Each experimental group consisted of three wells and the experiment was repeated three times.
Western blot analysis
Total protein lysates extracted from samples were quantitated using BCA protein assay kit (Pierce, Rockford, IL). 20 μg total protein lysates extracted from samples were separated with 10% sodium dodecyl sulfate-polyacryl-amide gels and transferred to a polyvinylidene fluoridemembrane. The membrane was blocked with 3% BSA, followed by incubation with mouse monoclonal anti-human WIF-1 (Abcam, Cambridge, MA), rabbit polyclonal anti-human E-cadherin, β-catenin and vimentin (Cell Signaling Technology Inc., Danvers, MA), anti-human β-Actin antibody (Epitomics, Burlingame, CA). The membrane was then incubated with secondary horseradish peroxidase-conjugated goat anti-rabbit or anti-mouse antibody (Jackson Immune Research Laboratories Inc., West Grove, PA), and visualized using Immobilon™ Western Chemiluminescent HRP substrate (Millipore). β-Actin was used as the loading control.
Immunohistochemical analysis on Tissue Microarray (TMA)
Immunohistochemical studies of WIF-1 were done on tissue microarray using the Polymer Detection System for immune-histological staining (PV9000 kit, ZSJQ-Bio, China). All images were examined by two experienced pathologists independently. Immunostaining was evaluated as described  previous.
SW480, SW480-NC, SW480-181a, SW620, SW620-NC, SW620-181a RNAi cells were plated in a 96-well microplate and were cultured overnight until cells were adherent. Cells were fixed and stained with anti-E-cadherin, −β-catenin and-vimentin primary antibody and DyLight549 Goat Anti-rabbit Secondary Antibody. Cells were also labeled with Hoechst 33342 Dye for nuclear staining. Expression of E-cadherin, β-catenin and vimentin was imaged using a laser scanning confocal microscopy (TCS-SP2, Leica Microsystems).
All data are presented as the mean ± standard error (SE) and were analyzed using the statistical package SPSS13.0. Data are expressed as median values (interquartile range) in supplementary tables or as means ± SE in the figures. Comparisons were performed using the nonparametric Mann–Whitney U test or Kruskal-Walis test for continuous variables. Spearman’s rho correlation test (two tailed) was used to estimate the correlation between the expression in CRC tissue of each marker and clinical-pathological variables. Survival curves were plotted using the Kaplan-Meier analysis and the log-rank test was used to test for significant differences. Cox univariate and multivariate proportional hazard models were used to estimate the hazard ratio for each marker. P ≤ 0.05 was considered statistically significant.
Jemal A, Siegel R, Xu J, Ward E: Cancer statistics, 2010. CA. 2010, 60 (5): 277-300.
Manfredi S, Lepage C, Hatem C, Coatmeur O, Faivre J, Bouvier AM: Epidemiology and management of liver metastases from colorectal cancer. Ann Surg. 2006, 244 (2): 254-259. 10.1097/01.sla.0000217629.94941.cf
Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nat Rev Genet. 2008, 9 (2): 102-114.
Ma L, Weinberg RA: Micromanagers of malignancy: role of microRNAs in regulating metastasis. Trends Genet. 2008, 24 (9): 448-456. 10.1016/j.tig.2008.06.004
Zhang W, Dahlberg JE, Tam W: MicroRNAs in tumorigenesis: a primer. Am J Pathol. 2007, 171 (3): 728-738. 10.2353/ajpath.2007.070070
Arndt GM, Dossey L, Cullen LM, Lai A, Druker R, Eisbacher M, Zhang C, Tran N, Fan H, Retzlaff K, Bittner A, Raponi M: Characterization of global microRNA expression reveals oncogenic potential of miR-145 in metastatic colorectal cancer. BMC Cancer. 2009, 9: 374- 10.1186/1471-2407-9-374
Ma Y, Zhang P, Wang F, Zhang H, Yang J, Peng J, Liu W, Qin H: miR-150 as a potential biomarker associated with prognosis and therapeutic outcome in colorectal cancer. Gut. 2012, 61 (10): 1447-1453. 10.1136/gutjnl-2011-301122
Schee K, Fodstad O, Flatmark K: MicroRNAs as biomarkers in colorectal cancer. Am J Pathol. 2010, 177 (4): 1592-1599. 10.2353/ajpath.2010.100024
Nishimura J, Handa R, Yamamoto H, Tanaka F, Shibata K, Mimori K, Takemasa I, Mizushima T, Ikeda M, Sekimoto M, Ishii H, Doki Y, Mori M: microRNA-181a is associated with poor prognosis of colorectal cancer. Oncol Rep. 2012, 28 (6): 2221-2226.
Lee EJ, Gusev Y, Jiang J, Nuovo GJ, Lerner MR, Frankel WL, Morgan DL, Postier RG, Brackett DJ, Schmittgen TD: Expression profiling identifies microRNA signature in pancreatic cancer. Int J Canc. 2007, 120 (5): 1046-1054.
Wang Y, Yu Y, Tsuyada A, Ren X, Wu X, Stubblefield K, Rankin-Gee EK, Wang SE: Transforming growth factor-beta regulates the sphere-initiating stem cell-like feature in breast cancer through miRNA-181 and ATM. Oncogene. 2011, 30 (12): 1470-1480. 10.1038/onc.2010.531
Shi L, Cheng Z, Zhang J, Li R, Zhao P, Fu Z, You Y: hsa-mir-181a and hsa-mir-181b function as tumor suppressors in human glioma cells. Brain Res. 2008, 1236: 185-193.
Zhu DX, Miao KR, Fang C, Fan L, Zhu W, Zhu HY, Zhuang Y, Hong M, Liu P, Xu W, Li JY: Aberrant microRNA expression in Chinese patients with chronic lymphocytic leukemia. Leuk Res. 2011, 35 (6): 730-734. 10.1016/j.leukres.2010.11.005
Ji J, Yamashita T, Budhu A, Forgues M, Jia HL, Li C, Deng C, Wauthier E, Reid LM, Ye QH, Qin LX, Yang W, Wang HY, Tang ZY, Croce CM, Wang XW: Identification of microRNA-181 by genome-wide screening as a critical player in EpCAM-positive hepatic cancer stem cells. Hepatology. 2009, 50 (2): 472-480. 10.1002/hep.22989
Meng F, Glaser SS, Francis H, DeMorrow S, Han Y, Passarini JD, Stokes A, Cleary JP, Liu X, Venter J, Kumar P, Priester S, Hubble L, Staloch D, Sharma J, Liu CG, Alpini G: Functional analysis of microRNAs in human hepatocellular cancer stem cells. J Cell Mol Med. 2012, 16 (1): 160-173. 10.1111/j.1582-4934.2011.01282.x
Lee BB, Lee EJ, Jung EH, Chun HK, Chang DK, Song SY, Park J, Kim DH: Aberrant methylation of APC, MGMT, RASSF2A, and Wif-1 genes in plasma as a biomarker for early detection of colorectal cancer. Clin Canc Res. 2009, 15 (19): 6185-6191. 10.1158/1078-0432.CCR-09-0111.
Wissmann C, Wild PJ, Kaiser S, Roepcke S, Stoehr R, Woenckhaus M, Kristiansen G, Hsieh JC, Hofstaedter F, Hartmann A, Knuechel R, Rosenthal A, Pilarsky C: WIF1, a component of the Wnt pathway, is down-regulated in prostate, breast, lung, and bladder cancer. J Pathol. 2003, 201 (2): 204-212. 10.1002/path.1449
Yee DS, Tang Y, Li X, Liu Z, Guo Y, Ghaffar S, McQueen P, Atreya D, Xie J, Simoneau AR, Hoang BH, Zi X: The Wnt inhibitory factor 1 restoration in prostate cancer cells was associated with reduced tumor growth, decreased capacity of cell migration and invasion and a reversal of epithelial to mesenchymal transition. Mol Cancer. 2010, 9: 162- 10.1186/1476-4598-9-162
Hu J, Dong A, Fernandez-Ruiz V, Shan J, Kawa M, Martinez-Anso E, Prieto J, Qian C: Blockade of Wnt signaling inhibits angiogenesis and tumor growth in hepatocellular carcinoma. Cancer Res. 2009, 69 (17): 6951-6959. 10.1158/0008-5472.CAN-09-0541
Kawakami K, Hirata H, Yamamura S, Kikuno N, Saini S, Majid S, Tanaka Y, Kawamoto K, Enokida H, Nakagawa M, Dahiya R: Functional significance of Wnt inhibitory factor-1 gene in kidney cancer. Cancer Res. 2009, 69 (22): 8603-8610.
Thiery JP: Epithelial-mesenchymal transitions in tumour progression. Nat Rev Canc. 2002, 2 (6): 442-454. 10.1038/nrc822.
Yang J, Mani SA, Donaher JL, Ramaswamy S, Itzykson RA, Come C, Savagner P, Gitelman I, Richardson A, Weinberg RA: Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell. 2004, 117 (7): 927-939. 10.1016/j.cell.2004.06.006
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185
Packeisen J, Korsching E, Herbst H, Boecker W, Buerger H: Demystified…tissue microarray technology. Mol Pathol. 2003, 56 (4): 198-204. 10.1136/mp.56.4.198
Dillhoff M, Liu J, Frankel W, Croce C, Bloomston M: MicroRNA-21 is overexpressed in pancreatic cancer and a potential predictor of survival. J Gastrointest Surg. 2008, 12 (12): 2171-2176. 10.1007/s11605-008-0584-x
Kong KL, Kwong DL, Chan TH, Law SY, Chen L, Li Y, Qin YR, Guan XY: MicroRNA-375 inhibits tumour growth and metastasis in oesophageal squamous cell carcinoma through repressing insulin-like growth factor 1 receptor. Gut. 2012, 61 (1): 33-42. 10.1136/gutjnl-2011-300178
Buchanan FG, Gorden DL, Matta P, Shi Q, Matrisian LM, DuBois RN: Role of beta-arrestin 1 in the metastatic progression of colorectal cancer. Proc Nat Acad Sci USA. 2006, 103 (5): 1492-1497. 10.1073/pnas.0510562103
Ma Y, Peng J, Liu W, Zhang P, Huang L, Gao B, Shen T, Zhou Y, Chen H, Chu Z, Zhang M, Qin H: Proteomics identification of desmin as a potential oncofetal diagnostic and prognostic biomarker in colorectal cancer. Mol Cell Proteomics. 2009, 8 (8): 1878-1890. 10.1074/mcp.M800541-MCP200
This work was supported by the National Natural Science Foundation (30872467, 81030040, 81201965), the Program for New Century Excellent Talents in Universities (NCET-BMU20100010), Beijing Natural Science Foundation (7132052), Research Fund for the Doctoral Program of Higher Education of China (RFDP) 20100001120127 the National High Technology Research and Development Program of China (863 Program) (No.2012AA02A506 ), and Beijing Municipal Administration of Hospitals Clinical Medicine Development of Special Funding Support (code ZY201410). The authors would like to thank Dr Bin Dong, the Department of Pathology, Peking University Cancer Hospital & Institute for her technical assistance.
The authors declare that they have no competing interests.
D-BJ contributed to study design, data collection, analysis, and interpretation; Z-GC performed tissue microarrays and in situ hybridization analysis; ML carried out immunoassays; T-CZ contributed to validation of target gene; Y-FY contributed to statistical analysis; J-ZX, LY and Z-QZ contributed to the discussion of the results; JG outlined the experimental design of the study, contributed to the discussion of the results and revised the manuscript. All authors contributed to writing, critical review, and final approval to submit the paper for publication.
Dengbo Ji, Zhiguo Chen, Ming Li contributed equally to this work.
Electronic supplementary material
Additional file 1: Table S1: All detected miRNAs expression in CRC patients without liver metastasis compared to patients with liver metastasis by microarray profiling screen. . The miRBase accession number and the mature sequence of the microRNA studied. Table S3. Clinical characteristics of colorectal cancer patients of Testing Cohort. Table S4. Clinical characteristics of colorectal cancer patients of Validation Cohort 1 and 2. Table S5. Primer sequences.” We highlight the paragraph in the resubmitted version. (DOC 364 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ji, D., Chen, Z., Li, M. et al. MicroRNA-181a promotes tumor growth and liver metastasis in colorectal cancer by targeting the tumor suppressor WIF-1. Mol Cancer 13, 86 (2014). https://doi.org/10.1186/1476-4598-13-86
- Colorectal cancer
- Liver metastasis