Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Highly Accessed Research article

Identifying MicroRNA-mRNA regulatory network in colorectal cancer by a combination of expression profile and bioinformatics analysis

Jihong Fu1, Wentao Tang2, Peng Du1, Guanghui Wang1, Wei Chen1, Jingming Li1, Yunxiang Zhu12, Jun Gao3* and Long Cui1*

Author Affiliations

1 Colorectal Surgery Department, Xinhua Hospital, School of Medicine, Shanghai Jiaotong University, Kongjiang Road 1665, Shanghai, 200092, China

2 Second Affiliated Hospital, Yangzhou University, Jiangsu, 225001, China

3 Department of Gastroenterology, Changhai Hospital, Second Military Medical University, Changhai Road 168, Shanghai, 200433, China

For all author emails, please log on.

BMC Systems Biology 2012, 6:68  doi:10.1186/1752-0509-6-68

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/6/68


Received:28 October 2011
Accepted:15 June 2012
Published:15 June 2012

© 2012 Fu et al.; licensee BioMed Central Ltd

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

MicroRNAs (miRNAs) are involved in carcinogenesis and tumor progression by regulating post-transcriptional gene expression. However, the miRNA-mRNA regulatory network is far from being fully understood. The objective of this study is to identify the colorectal cancer (CRC) specific miRNAs and their target mRNAs using a multi-step approach.

Results

A multi-step approach combining microarray miRNA and mRNA expression profile and bioinformatics analysis was adopted to identify the CRC specific miRNA-mRNA regulatory network. First, 32 differentially expressed miRNAs and 2916 mRNAs from CRC samples and their corresponding normal epithelial tissues were identified by miRNA and mRNA microarray, respectively. Secondly, 22 dysregulated miRNAs and their 58 target mRNAs (72 miRNA-mRNA pairs) were identified by a combination of Pearson’s correlation analysis and prediction by databases TargetScan and miRanda. Bioinformatics analysis revealed that these miRNA-mRNAs pairs were involved in Wnt signaling pathway. Additionally, 6 up-regulated miRNAs (mir-21, mir-223, mir-224, mir-29a, mir-29b, and mir-27a) and 4 down-regulated predicted target mRNAs (SFRP1, SFRP2, RNF138, and KLF4) were selected to validate the expression level and their anti-correlationship in an extended cohort of CRC patients by qRT-PCR. Except for mir-27a, the differential expression and their anti-correlationship were proven. Finally, a transfection assay was performed to validate a regulatory relationship between mir-29a and KLF4 at both RNA and protein levels.

Conclusions

Seventy-two miRNA-mRNA pairs combined by 22 dysregulated miRNAs and their 58 target mRNAs identified by the multi-step approach appear to be involved in CRC tumorigenesis. The results in our study were worthwhile to further investigation via a functional study to fully understand the underlying regulatory mechanisms of miRNA in CRC.

Keywords:
MicroRNAs; Target mRNAs; Colorectal cancer; Microarray

Background

MicroRNAs (miRNAs) are one family of small (~22 nucleotides), non-coding RNAs. They regulate gene expression post-transcriptionally through binding to the complementary sites of target mRNAs in the 3’UTR [1], and play an important role in regulating diverse biological processes, such as stem cell maintenance, development, metabolism, proliferation, differentiation and apoptosis [2]. Cumulative evidence indicates that some miRNAs act as either oncogenes or tumor suppressors in the tumorigenesis of CRC and possess the tumor marker potential for diagnostic, therapeutic, prognostic exploration [3-9]. To elucidate the mechanism of miRNA regulation, it is critical to discover the differentially expressed miRNA and identify the target mRNAs of miRNAs.

Many aberrant expressed miRNAs were found in CRC tissues by miRNA microarray assay [10-15]. The target relationship between particular miRNA and mRNA were identified by in-vitro experiments such as transfection and luciferase activity assays [16-19]. Although in-vitro experiments provide invaluable information on miRNA regulation, those experiments are not possible to be applied to clinical samples, and it is also impossible to obtain results under a genome-wide circumstance. Many databases were explored to predict the targets of miRNAs by either seed sequence complementary or thermodynamic algorithm [20-23], however these databases only provide the possibility of predicting direct target-relationship miRNA-mRNA pairs and can not predict disease-specific pairs in a specific disease.

MiRNAs regulate target gene expression in two ways, mRNA degradation or translation inhibition [24,25]. The target genes which were regulated by miRNAs through mRNA degradation are negatively correlated with the miRNA regulators. Therefore an approach for the identification of miRNA-mRNA regulatory modules was introduced by constructing the paired miRNA-mRNA expression profiles and predicting the putative target genes of miRNAs [26,27]. However, this novel multi-step approach has not been tested in CRC, especially in clinical matched tumor/normal mucosa samples.

In this study, we applied the multi-step approach by constructing the aberrant paired miRNA-mRNA expression profiles and the clusters of miRNA target genes to investigate the role of miRNA in CRC.

Results

Strategy

As shown in Figure 1, a multi-step approach was adopted to identify the mRNA targets of CRC-specific miRNAs in colorectal cancer. First, in the same cohort of eight cases of pathologically confirmed CRC and their corresponding adjacent normal tissues, the significantly aberrant expression profiles of both miRNA and mRNA were obtained by microarray test. Second, the significant anti-correlation between miRNA-mRNA pairs was determined based on the expression level of each aberrantly expressed miRNA and each aberrantly expressed mRNA. Third, only the target-relationship miRNA-mRNA pairs, which were also predicted by both TargetScanHuman V5.1 [20] and miRanda [23], were retained. Fourth, some of the predicted direct target-relationship miRNA-mRNA pairs were selected to be validated in an extended cohort of CRC tissues by qRT-PCR detection.

thumbnailFigure 1. Step 1, both miRNA and mRNA microarray tests for identifying significantly dysregulated miRNAs and mRNAs in colorectal cancer. Step 2, Pearson correlation analysis for anti-correlationship of dysregulated miRNAs and mRNAs. Step 3, TargetScan and Mirnada predicting the target-relationship miRNA-mRNA pairs.

Further in detailed explanation, the target genes of dysregulated miRNA in CRC tissues should meet the following requirements: 1)the microRNAs and mRNAs were differentially expressed between CRC tissues and corresponding normal samples 2) the expression level of miRNAs and mRNAs were significantly anti-correlated (P < 0.05) 3) the mRNA is a putative target of miRNA predicted by both TargetScan [20] and miRanda [23].

miRNA and mRNA aberrant expression profiles in CRC samples

The miRNA expression profiles in 8 CRC tissues and their corresponding adjacent normal tissues were determined using miRNA microarray analysis. Histologically, the 8 CRC cases were all moderately differentiated, non-mucinous adeno-carcinoma, micro-satellite stable. The clinical data of the 8 CRC cases are listed in Additional file 1: Table S1. After between-array normalization and filtering on flags (IsGeneDetected, WellAboveNeg) carried out by RMA algorithm [28], 125 miRNAs were selected for further analysis. Using Student’s t-test and Significant Analysis of Microarray analyses (SAM) [29] (FDR < 0.05 from student’s t-test and q-value < 0.05 from SAM), 32 miRNAs were found to be significantly differentially expressed between CRC tissues and their corresponding adjacent normal tissues, including 16 up-regulated and 16 down-regulated which were shown in Additional file 1: Table S2. In addition, the aberrantly different miRNA expression levels were illustrated by Heatmap in Figure 2.

Additional file 1. Supplemental tables and figures

Format: PDF Size: 190KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 2. FDR < 0.05 for students t-test and q-value < 0.05 for SAM were considered as significantly differentially expressed. The data were z-score transformed and were indicated by the color bar below the heatmap.

The mRNA expression profiles of the same 8 paired-samples in miRNA expression profiling were analyzed by Whole Human Genome Agilent 4 × 44 K Oligonucleotide Microarray (Agilent). After between-array normalization performed by Linear Models for Microarray Data (R-package Limma) [30], 2916 genes were identified as significantly differentially expressed in CRC tissues using the criteria similar to that in miRNA microarray assays (FDR < 0.05 from student’s t-test and q-value < 0.05 from SAM) and was included in Additional file 2: Table S1.

Additional file 2. Differentially expressed mRNA in CRC tissues

Format: XLS Size: 693KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

The microarray data was deposited to the NCBI Gene Expression Omnibus (accession GSE35982)

Identifying target genes by anti-correlation and target prediction

Pearson’s correlation analysis was applied to the 32 significantly differentially expressed miRNAs and significantly differentially expressed 2916 mRNAs. Since miRNAs act as negative regulators, up-regulated miRNAs resulted in down-regulated target mRNAs, and vice versa. A total of 17360 miRNA-mRNA pairs (Additional file 3: Table S1) were obtained by two criteria: Pearson’s correlation less than -0.5 was considered as anti-correlation and a p-value less than 0.05 as statistically significant. A multiple-hypothesis correction was carried out by Benjamini-Horchberg algorithm. The FDR value of all the 17630 miRNA-mRNA pairs were less than 0.3 and were selected for further analysis. To predict the target mRNAs for the 32 dysregulated miRNAs, TargetScanHuman (V5.1.) and miRanda database were used. To eliminate false positive rates of the target prediction databases, only the miRNA-mRNA pairs simultaneously predicted by both TargetScanHuman V5.1 and miRanda were considered as target relationship miRNA-mRNA pairs in our study. Seventy-two miRNA-mRNA pairs combined by 22 dysregulated miRNAs and their 58 target mRNAs were predicted by those two prediction databases and simultaneously presented an anti-correlationship between the expression level of miRNAs and mRNAs. These 58 predicted target mRNAs included many genes previously reported to be involved in tumorigenesis of CRC or other malignant tumors such as KLF4 [31], SFRP1 [32], SFRP2 [33], RNF138 [34]. These 72 predicted miRNA-mRNA pairs contained 14 up-regulated miRNAs targeting 43 down-regulated mRNA (Additional file 4) and 8 down-regulated miRNAs targeting 15 up-regulated mRNAs (Additional file 5). The expression levels of 58 target mRNAs were illustrated by Heatmap (Additional file 1: Figure S1).

Additional file 3. Significantly anti-correlated miRNA-mRNA pairs

Format: XLS Size: 3.6MB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 4. Target mRNAs of up-regulated miRNAs

Format: XLS Size: 16KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 5. Target mRNAs of down-regulated miRNAs

Format: XLS Size: 12KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Pathway analysis of the miRNA targets

To obtain a better understanding of biological function of miRNAs dysregulated in CRC, function analysis of the miRNA target genes was carried out by Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 [35]. As a result of finding which biological processes of GO or pathway from the KEGG database were significantly enriched, 38 GO categories were found enriched with a p-value < 0.05 (Additional file 6), however using a cutoff of FDR < 0.05, no GO category were found significantly enriched. Using a cutoff of FDR < 0.05, the Wnt signaling pathway (hsa04310) was the only significantly enriched KEGG pathway (Table 1).

Additional file 6. Significantly enriched GO biological processes

Format: XLS Size: 17KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Table 1. Wnt signaling pathway was enriched by the miRNA targets

Validation of the expression level of miRNA and mRNA in an extended cohort of CRC tissues by qRT-PCR

To validate the reliability of the microarray test,6 up-regulated miRNAs (mir-21, mir223, mir224, mir29a, mir29b, and mir27a) and their 4 down-regulated target mRNAs (SFRP1, SFRP2, RNF138, and KLF4) were selected for qRT-PCR examination in 40 matched CRC/normal mucosa tissues whose clinical data were listed in Additional file 1: Table S3. The rationale for selecting these up-regulated miRNAs and down-regulated mRNAs was the down-regulated expression of SFRP1 [32,36], SFRP2 [32,33,37,38], and RNF138 [34] were reported to be involved in the activity of the Wnt receptor signaling pathway. Further it was found that the expression of KLF4 [31] was down-regulated in CRC tissues compared with normal mucosa.

As demonstrated in Figure 3A, the five miRNAs (mir-21, mir-223, mir-224, mir-29A and mir-29B) were significantly up-regulated in CRC tissues, except for mir-27a due to non-efficient amplification (CT value in most of the samples was either higher than 35 or even cannot be defined). As shown in Figure 3B, the expression levels of the 4 mRNAs (SFRP1, SFRP2, RNF138, KLF4)were significantly down-regulated in CRC tissues.

thumbnailFigure 3. A. The expression levels of 5 miRNAs were significantly up-regulated. B. The expression levels of 4 mRNA targets in CRC samples and adjacent normal tissues. The expression levels were expressed as –ΔΔCT after normalized with beta-actin. P-value was calculated by paired t-test.

Additionally, the Pearson’s correlation coefficient between the expression levels of miRNA-mRNA pairs were calculated, and two identified miRNA-mRNA pairs, mir-224 and SFRP2, mir-29a and KLF4 were found significantly anti-correlated in qRT-PCR tests. However, 2 miRNA-mRNA pairs (mir-29a and SFRP1; mir-21 and SFRP1) not predicted by the two exist between miRNAs and SFRP1 (Table 2).

Table 2. Validation of the anti-correlationship between miRNAs and mRNAsexpression levels by qRT-PCR

In vitro analysis revealed expression of KLF4 is regulated by mir-29a

To further validate identified miRNA-mRNA interactions, the transfection assay was carried out to demonstrate the relationship between mir-29a and its identified target, KLF4. HCT-116 cells were transfected with mir-29a mimics and mir-29a inhibitor. The qRT-PCR and Western-blot revealed that the KLF4 mRNA and protein expression were significantly down-regulated by the mir-29a mimics (Figure 4a) while the KLF4 expression levels were significantly up-regulated by mir-29a inhibitor (Figure 4b), compared with the control group and mock group. These data indicated that the KLF4 mRNA and protein expression levels are controlled by mir-29a, suggesting that KLF4 probably was the target of mir-29a.

thumbnailFigure 4. A. The transfection of mir-29a mimics into HCT-116 cells down-regulates KLF4 mRNA expression (Left, QRT-PCR) and protein expression (Right, Western blot). B. The transfection of mir-29a inhibitor into HCT-116 cells up-regulates KLF4 mRNA expression (Left, QRT-PCR) and protein expression (Right, Western blot). The experiments were independently repeated three times and GADPH served as control. ***, P < 0.001.

Discussion

miRNAs, acting as negative regulators of gene expression, regulate various biological processes via inhibiting the expression of their target genes. Increasing evidence shows that some miRNAs play critical roles in tumorigenesis of CRC and have a potential clinical value in diagnosis, treatment and prognosis evaluation for CRC [4-9]. It is of great importance to identify the targets of miRNAs and elucidate their complex regulatory networks.

In this study, a multi-step approach to identify the mRNA targets of dysregulated miRNAs in CRC was adopted. Firstly, in the same cohort of samples containing CRC tissues and their corresponding adjacent tissues, the expression profiles of miRNAs and mRNAs were screened. Secondly, all possible pairs of miRNA-mRNA with an anti-correlationship were identified by Pearson’s correlation analysis, implying that these pairs of miRNA-mRNA may contain both the direct target-relationship and indirect anti-correlationship pairs. Thirdly, the reliability of the identified target-relationship pairs of miRNAs-mRNAs was assessed and confirmed by using both TargetScan and miRanda databases. Finally, the differential expression and their anti-correlationship of 6 miRNAs and 4 mRNAs were validated by qRT-PCR and the regulatory effect of mir-29a on the expression of target of KLF4 mRNA was validated by in-vitro analysis. Through the multi-step approach, 72 miRNA-mRNA pairs containing 22 dysregulated miRNAs and their 58 target mRNAs in CRC were identified, and the Wnt signaling pathway was the only affected pathway based on pathway analysis for the miRNA targets.

Compared with previous studies on the miRNA expression profile of CRC, including miRNA microarray [12-15], miRNA cloning and SAGE [39], most dysregualted miRNAs detected in this study were consistent with at least one of the previous studies. However, some discrepancy and inconsistency of the miRNA expression profile between studies were found. For instance, the elevated expression of mir-29b in CRC tissues in this study and the recent report [14] contradicted findings in Cummins’ report [39] and the elevated expression of mir-424 in CRC tissues in this study was in consistent with the other two studies based on Chinese patients’ samples[11,40]. However, no other studies reported the altered expression level of mir-424 in CRC. The discrepancy of miRNA expression profile between studies might be attributed to the sample size, genetic background, clinical characteristics or different technique platforms. Strikingly, it was found that mir-7, an up-regulated miRNA with high fold -change in this study (fold-change = 3.60) was not reported as a dysregulated miRNA in CRC in previous studies. In studies with other tumors, an elevated expression level of mir-7 was an indicator of poor disease-free survival in non-small cell lung carcinoma [41], and correlated with a larger size of tumor and poor tumor grade in breast cancer tissues [42] while some reports stated that mir-7 inhibits EGFR and the AKT pathway, and the expression level of mir-7 was down-regulated expression in gilobastoma [43,44]. It seems that the role of mir-7 differs in various types of cancer, suggesting that the expression and function of miRNAs is tissue cell specific. The expression level of mir-7 and its potential targets in CRC are worthwhile for further investigation

In this study, an important note should be added that a very strict criterion was adopted to identify the miRNA targets. Briefly, the target mRNA should be differentially expressed, significantly anti-correlated with its miRNA regulator, and predicted by two databases, TargetScan and miRanda. Based on this strict criterion,our final identified miRNA-mRNA pairs have a high probability of involvement in CRC tumorigenesis. The following analysis will support this high probability. On one hand, comparing the mRNA targets (n = 4338) of the dysregualted miRNAs (n = 32) predicted just by TargetScan, only a small portion of the target genes (244/4338, 5.6 %) were significantly differentially expressed in CRC, and were selected for further analysis, while a great number of mRNAs (4094/4338, 94.4 %) were not differentially expressed in CRC tissue, and were not selected for further analysis. The fact that these 4094 genes were predicted as miRNAs targets but not significantly differentially expressed, suggested that the expression of these 4094 mRNAs might not be regulated by the 32 dysregulated miRNAs, but by other mechanisms such as regulation by other miRNA, methylation, cis-regulatory elements, etc. which may act and counteract the regulation of miRNA. On the other hand, Pearson’s correlation analysis selected 17360 miRNA-mRNA pairs whose expression levels were significantly inversely correlated, however only 72 target relationships of miRNA were finally identified. The existence of the rest of 17288 miRNA-mRNA pairs may be caused by the indirect regulation.

Another interesting finding was that the expression values of mir-29a showed a narrower range of change than those of its predicted target, KLF4. To rationalize this finding, it is hypothesized that the KLF4 mRNA expression was regulated by many other factors including other miRNAs. Besides mir-29a, KLF4 was also identified as a target gene of mir-7, an up-regulated miRNA with high fold-change. The unparalleled change of miRNA and its target mRNA expression levels could be caused by the two miRNA (mir-29a and mir-7) targeting one gene (KLF4). Inspired by this finding, the relationship between fold-change of mRNAs expression and the number of corresponding miRNA regulators were further analyzed. In the identified 72 miRNA-mRNA pairs, three groups can be classified: 1) 2 mRNAs were regulated by 3 miRNAs, 2)10 mRNAs were regulated by 2 miRNAs, and 3) the rest (60 mRNAs) were regulated by 1 miRNA. However, the fold-changes of the gene mRNA expression between these three groups do not reach the statistical significance, indicating the gene expression regulation is very complex and dynamic with impact from multi-factors.

miRNAs regulate diverse biological processes via regulating the expression of target genes. In this study, only the Wnt signaling pathway, which had been proved to be aberrantly activated in CRC [45,46], was found to be affected through pathway analysis of the miRNA targets. Our study revealed that 4 down-regulated mRNA and an up-regulated miRNA regulated by 5 dysregulated miRNA in CRC respectively were involved in the Wnt signaling pathway, and shed light on the miRNA regulation of Wnt signaling pathway.

Nevertheless, the orchestration of these regulatory pairs of miRNA-mRNA differed in different diseases, presenting the disease-specific expression profiles and composing a complex regulatory network. Although many methods had been used to identify the regulatory pairs of miRNA-mRNA in diseases such as microarray and software-prediction, no method possessed the capacity to determine all miRNA-mRNA regulatory pairs. Moreover, since the anti-correlationship of the miRNAs-mRNA expression levels was used as the screening criteria, only those target genes regulated through mRNA degradation could be identified here, target genes regulated by miRNAs through translation inhibition would be lost.

In summary, these new findings of pairs of miRNA-mRNA provide some hints for the mechanisms of CRC tumorigenesis and were worthwhile for further functional analysis. The more complex regulatory networks call for more research in the future.

Conclusions

Based on this research, seventy-two miRNA-mRNA pairs combined by 22 dysregulated miRNAs and their 58 target mRNAs identified by the multi-step approach appear to be involved in CRC tumorigenesis. The identified miRNA-mRNA pairs are worthwhile for further functional studies.

Methods

Patients and clinical samples

Forty CRC tissues and matched adjacent normal colorectal mucosa tissues (at least 5 cm from the cancer) were obtained from patients who underwent surgical excision for CRC at Xinhua Hospital, School of Medicine, Shanghai Jiaotong University. The samples were collected according to the following histological prerequisites: the CRC tissues had more than 80 % tumor cells; the matched adjacent normal colorectal mucosa tissues had normal mucosal structure without dysplastic cells. The samples were snap frozen in liquid nitrogen in 5 minutes after surgical resection and were then stored at −80°C temperature until RNA extraction. Informed consents were obtained from all the patients and the study was approved by hospital ethics committee.

RNA extraction

Total RNA was extracted and isolated using mirVana TM RNA Isolation Kit (Applied Biosystem p/n AM1556) following the manufacturer’s instructions. The concentration of RNA was measured by NanoDrop spectrophotometer. Total RNA concentration ranged form 100 ng/ul to 1ug/ul. The RNA quality was evaluated by Agilent 2100 Bioanalyzer; RNAs with RIN ≥ 7.0 and 28 S/18 S > 0.7 were accepted for microarray analysis.

MicroRNA microarray procedure

8 colorectal cancer tissues and their corresponding adjacent normal tissues ere selected for microRNA microarray analysis from the 40 matched tissues pairs of CRC tissues and their adjacent normal colorectal mucosa tissues. These 8 pairs of samples shared similar clinical and pathological features that they were all the moderately differentiated, non-mucinous adeno-carcinoma and micro-satellite stable. The total RNA was analyzed by microarray (Agilent V12.0) according to the manufacture’s protocol. In brief, 100 ng of total RNA was labeled with Cy3, and then hybridized on a miRNA array (8*15 K,Agilent V12.0) using Agilent’s miRNA Complete Labeling and Hyb Kit (p/n 5190–0456). Microarray after hybridization was washed by GE Wash Buffer 1 and GE Wash Buffer 2 (Gene Expression Wash Buffer kit, Agilent p/n 5188–5327). The slides were then scanned in microarray scanner by Agilent scan control software Version A7.0. The data collection, background subtraction and with-in array normalization were carried out by Agilent Feature Extraction (FE) software version 9.5.3. After Between-array normalization and signal filtration (IsGeneDetected, WellAboveNeg) performed by RMA algorithm [28], the data was analyzed by T-test and Significant Analysis of Microarray analyses(SAM) [29] to identify differentially expressed miRNAs. miRNAs with false discovery rate (FDR) < 0.05 using student’s t-test and q-value < 0.05 for SAM were identified as differentially expressed genes.

mRNA microarray procedure

An amount of 2 ug of total RNA isolated from CRC tissues and corresponding normal tissues were reverse-transcribed into cDNA, and labeled with Cy3 or Cy5. The labeled cDNA was then hybridized on a Whole Human Genome Agilent 4X44K Oligonucleotide Microarray (Agilent) according to the manufacture’s protocol using miRNA Complete Labeling and Hyb Kit. The chips were washed using Agilent Gene Expression Wash Buffer Kit (Agilent 5188–5327). Scanning of the microarrays was carried out immediately using the Agilent G2565BA Microarray Scanner System. Images were auto gridded, analyzed and data extracted using Agilent Feature Extraction software (version 9.5.3). Raw expression data generated by the Feature Extraction software was imported into R using the LIMMA package in Bioconductor [30]. The intensity distributions within and between arrays were normalized using the quantile scaling algorithm in LIMMA, and the log-2 based data was then filtered on flags (IsFound, IsWellAboveBG, IsSaturated) to identify differentially expressed features between colorectal cancer sample and normal colon mucosa. T-test and Significant Analysis of Microarray analyses (SAM) [29] were carried out to identify dysregulated mRNAs.

Similar to the criteria of identifying dysregulated miRNAs, mRNAs with FDR <0.05 for student’s t-test and q-value < 0.05 for SAM were identified as differentially expressed genes.

Target predicting

TargetScanHuman v. 5.1 [20] (http://www.targetscan.org/ webcite; Whitehead institute for Biomedical Research) and Mirnada [23] were used to predict the target genes of miRNAs

Pathway enrichment analysis

The pathway enrichment analysis was performed using online database. The Database for Annotation, Visualization and Integrated Discovery (DAVID) v 6.7 [35].

qRT-PCR

Both mRNA and miRNA expression level were detected in 40 CRC tissues and their corresponding adjacent normal tissues by qRT-PCR method using the ABI Detection Kit. For mRNA expression detection, 2 ug of total RNA was reverse transcribed by MMLV transcription according to manufacture’s protocol. The cDNA was analyzed by TaqMan assays. qRT-PCR was performed by ONE-STEP PLUS system (ABI). The primers and probes of qRT-PCR were listed in Additional file 1: Table S3. For miRNAs expression detection, a 15ul of reverse transcription mixture containing 3ul of RT primer, 5 ng of total RNA and 50U of MultiScribe Reverse Transcriptase, 1.5ul Reverse Transcription Buffer (10×), 0.15ul 1 dNTPs (00 mM), 4U of RNase Inhibitor was incubated in PCR tube for 16°C 30 min, 42°C 30 min, 85°C 5 min. The cDNA product were then stored at −4°C. The real-time PCR was performed by ONE-STEP PLUS system (ABI). A total volume of 20 ul reaction mixture included 1.0ul TaqMan Small RNA Assay (20×) ,1.33 ul of Product from RT reaction, 10.0 ul TaqMan Universal PCR Master Mix II (2×), 7.67 ul Nuclease-free water. The reactions were incubated in 96-well optical plate for 95°C 10 min; 95°C 15 s, 60°C 60s with 40 cycles.

The expression level of each mRNA and miRNA was calculated by threshold cycle (CT). The relative expression level was calculated using –ΔΔCT. A median expression sample among all samples was chosen as a calibrator, and U6 as endogenous control for miRNAs expression, beta-actin for mRNAs expression.

Cell culture

Human colorectal cell line HCT 116 was purchased from the Type Culture Collection of the Chinese Academy of Sciences, Shanghai, China. The cell line was incubated in Dulbecco's modified Eagle's medium (DMEM, Invitrogen, Carlsbad, CA, USA) with 10 % FBS.

Transfection

The hsa-mir-29a mimics, mimics control, hsa-mir-29a inhibitor (anti–miR-95, chemically modified antisense oligonucleotides designed to target specifically against mature miR-29a), miRNA inhibitor control were designed and synthesized from Ribobio (Guangzhou, China). The 4 × 106 of HCT 116 cell were grown in 6-well plates for 24 hours before transfection. The transfection was performed by Roche X-tremeGENE DNA Transfection Reagents according to manufacturer’s protocol. Every experiment was performed in triplicates in 6-well cell culture plates with the appropriate controls. Cells were harvested 48 h after transfection.

Western blot

Total proteins from cultured cells were extracted by RIPA Lysis Buffer (Beyotine, Shanghai, China). 40 ug of proteins were electrophoresed through 10 % SDS polyacrylamide gels and were transferred to a PVDF membranes (Millipore) The membranes were incubated with anti-KLF4 Rabbit polyclonal antibody (ab26648,1:1000, abcam), anti-GADPH antibody was chosen as an endogenous control. Membranes were probed with the appropriate secondary antibodies linked to horseradish peroxidase. Antibody binding was visualized using ECL.

Statistical analysis

Statistical analysis was performed using SPSS15.0. Values were expressed as mean ± standard deviation for parametric data. Difference between groups was compared using student’s t-test. Correlation between miRNA and mRNA expression profiles was calculated using Pearson’s correlation. P value less than 0.05 was considered as statistically significant. The statistical analysis of microarray analysis was described elsewhere.

Authors’ contributions

JF, JG and LC designed the study and authored the manuscript. PD, JL and WC collected the samples. GW and WT performed the qRT-PCR analysis. WT carried out transfection assays. PD, JF and JG carried out the microarray test and the bioinformatics analysis. LC, PD, JG and JF analyzed data. All authors read and approved the final manuscript. JF and WT contributed equally to this work.

Acknowledgements

The project was supported by the Major Program for the Fundamental Research of the Shanghai Committee of Science and Technology (Grant No.10DJ1400504), the Natural Science Foundation of Shanghai (Grant No. 11ZR1424200) and the Key Program for the Fundamental Research of the Shanghai Committee of Science and Technology (Grant No. 09JC1410900). The support on transfection assay from Professor Jiuhong Kang of laboratory of stem cell and epigenetic research, School of Life Science and Technology, Tongji University, Shanghai is appreciated. Thanks for the technical support of microarray from Oebiotech Company and qRT-PCR test from Biotecan Company

References

  1. Gregory RI, Shiekhattar R: MicroRNA biogenesis and cancer.

    Cancer Res 2005, 65(9):3509-3512. PubMed Abstract | Publisher Full Text OpenURL

  2. Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer.

    Nat Rev Cancer 2006, 6(4):259-269. PubMed Abstract | Publisher Full Text OpenURL

  3. Schee K, Fodstad O, Flatmark K: MicroRNAs as biomarkers in colorectal cancer.

    Am J Pathol 2010, 177(4):1592-1599. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Xi Y, Formentini A, Chien M, Weir DB, Russo JJ, Ju J, Kornmann M: Prognostic Values of microRNAs in Colorectal Cancer.

    Biomark Insights 2006, 2:113-121. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Slaby O, Svoboda M, Michalek J, Vyzula R: MicroRNAs in colorectal cancer: translation of molecular biology into clinical application.

    Mol Cancer 2009, 8:102. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  6. Ahmed FE, Jeffries CD, Vos PW, Flake G, Nuovo GJ, Sinar DR, Naziri W, Marcuard SP: Diagnostic microRNA markers for screening sporadic human colon cancer and active ulcerative colitis in stool and tissue.

    Cancer Genomics Proteomics 2009, 6(5):281-295. PubMed Abstract | Publisher Full Text OpenURL

  7. Faber C, Kirchner T, Hlubek F: The impact of microRNAs on colorectal cancer.

    Virchows Arch 2009, 454(4):359-367. PubMed Abstract | Publisher Full Text OpenURL

  8. Huang Z, Huang D, Ni S, Peng Z, Sheng W, Du X: Plasma microRNAs are promising novel biomarkers for early detection of colorectal cancer.

    Int J Cancer 2010, 127(1):118-126. PubMed Abstract | Publisher Full Text OpenURL

  9. Li JM, Zhao RH, Li ST, Xie CX, Jiang HH, Ding WJ, Du P, Chen W, Yang M, Cui L: Down-regulation of fecal miR-143 and miR-145 as potential markers for colorectal cancer.

    Saudi Med J 2012, 33(1):24-29. PubMed Abstract OpenURL

  10. Motoyama K, Inoue H, Takatsuno Y, Tanaka F, Mimori K, Uetake H, Sugihara K, Mori M: Over- and under-expressed microRNAs in human colorectal cancer.

    Int J Oncol 2009, 34(4):1069-1075. PubMed Abstract | Publisher Full Text OpenURL

  11. Wang YX, Zhang XY, Zhang BF, Yang CQ, Chen XM, Gao HJ: Initial study of microRNA expression profiles of colonic cancer without lymph node metastasis.

    J Dig Dis 2010, 11(1):50-54. PubMed Abstract | Publisher Full Text OpenURL

  12. Mosakhani N, Sarhadi VK, Borze I, Karjalainen-Lindsberg ML, Sundstrom J, Ristamaki R, Osterlund P, Knuutila S: MicroRNA profiling differentiates colorectal cancer according to KRAS status.

    Genes Chromosomes Cancer 2012, 51(1):1-9. PubMed Abstract | Publisher Full Text OpenURL

  13. Chen X, Guo X, Zhang H, Xiang Y, Chen J, Yin Y, Cai X, Wang K, Wang G, Ba Y, et al.: Role of miR-143 targeting KRAS in colorectal tumorigenesis.

    Oncogene 2009, 28(10):1385-1392. PubMed Abstract | Publisher Full Text OpenURL

  14. Monzo M, Navarro A, Bandres E, Artells R, Moreno I, Gel B, Ibeas R, Moreno J, Martinez F, Diaz T, et al.: Overlapping expression of microRNAs in human embryonic colon and colorectal cancer.

    Cell Res 2008, 18(8):823-833. PubMed Abstract | Publisher Full Text OpenURL

  15. Ng EK, Chong WW, Jin H, Lam EK, Shin VY, Yu J, Poon TC, Ng SS, Sung JJ: Differential expression of microRNAs in plasma of patients with colorectal cancer: a potential marker for colorectal cancer screening.

    Gut 2009, 58(10):1375-1381. PubMed Abstract | Publisher Full Text OpenURL

  16. Huang Z, Huang S, Wang Q, Liang L, Ni S, Wang L, Sheng W, He X, Du X: MicroRNA-95 promotes cell proliferation and targets sorting Nexin 1 in human colorectal carcinoma.

    Cancer Res 2011, 71(7):2582-2589. PubMed Abstract | Publisher Full Text OpenURL

  17. Chang KH, Miller N, Kheirelseid EA, Ingoldsby H, Hennessy E, Curran CE, Curran S, Smith MJ, Regan M, McAnena OJ, et al.: MicroRNA-21 and PDCD4 expression in colorectal cancer.

    Eur J Surg Oncol 2011, 37(7):597-603. PubMed Abstract | Publisher Full Text OpenURL

  18. Sun K, Wang W, Zeng JJ, Wu CT, Lei ST, Li GX: MicroRNA-221 inhibits CDKN1C/p57 expression in human colorectal carcinoma.

    Acta Pharmacol Sin 2011, 32(3):375-384. PubMed Abstract | Publisher Full Text OpenURL

  19. Wang H, Wu J, Meng X, Ying X, Zuo Y, Liu R, Kang T, Huang W: MicroRNA-342 Inhibits Colorectal Cancer Cell Proliferation and Invasion by Directly Targeting DNA Methyltransferase 1.

    Carcinogenesis 2011, 7:1033-42. OpenURL

  20. Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets.

    Cell 2005, 120(1):15-20. PubMed Abstract | Publisher Full Text OpenURL

  21. Maziere P, Enright AJ: Prediction of microRNA targets.

    Drug Discov Today 2007, 12(11–12):452-458. PubMed Abstract | Publisher Full Text OpenURL

  22. John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets.

    PLoS Biol 2004, 2(11):e363. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org resource: targets and expression.

    Nucleic Acids 2008, 36:149-153. OpenURL

  24. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function.

    Cell 2004, 116(2):281-297. PubMed Abstract | Publisher Full Text OpenURL

  25. Cannell IG, Kong YW, Bushell M: How do microRNAs regulate gene expression?

    Biochem Soc Trans 2008, 36(Pt 6):1224-1231. PubMed Abstract | Publisher Full Text OpenURL

  26. Peng X, Li Y, Walters KA, Rosenzweig ER, Lederer SL, Aicher LD, Proll S, Katze MG: Computational identification of hepatitis C virus associated microRNA-mRNA regulatory modules in human livers.

    BMC Genomics 2009, 10:373. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  27. Liu H, Brannon AR, Reddy AR, Alexe G, Seiler MW, Arreola A, Oza JH, Yao M, Juan D, Liou LS, et al.: Identifying mRNA targets of microRNA dysregulated in cancer: with application to clear cell Renal Cell Carcinoma.

    BMC Syst Biol 2010, 4:51. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. Lopez-Romero P, Gonzalez MA, Callejas S, Dopazo A, Irizarry RA: Processing of Agilent microRNA array data.

    BMC Res Notes 2010, 3:18. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response.

    Proc Natl Acad Sci U S A 2001, 98(9):5116-5121. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Smyth GK, Speed TP: Normalization of cDNA microarray data.

    Methods 2003, 31:265-273. PubMed Abstract | Publisher Full Text OpenURL

  31. Ghaleb AM, Yang VW: The Pathobiology of Kruppel-like Factors in Colorectal Cancer.

    Curr Colorectal Cancer Rep 2008, 4(2):59-64. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Caldwell GM, Jones C, Gensberg K, Jan S, Hardy RG, Byrd P, Chughtai S, Wallis Y, Matthews GM, Morton DG: The Wnt antagonist sFRP1 in colorectal tumorigenesis.

    Cancer Res 2004, 64(3):883-888. PubMed Abstract | Publisher Full Text OpenURL

  33. Qi J, Zhu YQ, Luo J, Tao WH: Hypermethylation and expression regulation of secreted frizzled-related protein genes in colorectal tumor.

    World J Gastroenterol 2006, 12(44):7113-7117. PubMed Abstract | Publisher Full Text OpenURL

  34. Yamada M, Ohnishi J, Ohkawara B, Iemura S, Satoh K, Hyodo-Miura J, Kawachi K, Natsume T, Shibuya H: NARF, an nemo-like kinase (NLK)-associated ring finger protein regulates the ubiquitylation and degradation of T cell factor/lymphoid enhancer factor (TCF/LEF).

    J Biol Chem 2006, 281(30):20749-20760. PubMed Abstract | Publisher Full Text OpenURL

  35. da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.

    Nat Protoc 2009, 4(1):44-57. PubMed Abstract | Publisher Full Text OpenURL

  36. Caldwell GM, Jones CE, Taniere P, Warrack R, Soon Y, Matthews GM, Morton DG: The Wnt antagonist sFRP1 is downregulated in premalignant large bowel adenomas.

    Br J Cancer 2006, 94(6):922-927. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Huang D, Yu B, Deng Y, Sheng W, Peng Z, Qin W, Du X: SFRP4 was overexpressed in colorectal carcinoma.

    J Cancer Res Clin Oncol 2010, 136(3):395-401. PubMed Abstract | Publisher Full Text OpenURL

  38. MacDonald BT, Tamai K, He X: Wnt/beta-catenin signaling: components, mechanisms, and diseases.

    Dev Cell 2009, 17(1):9-26. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Cummins JM, He Y, Leary RJ, Pagliarini R, Diaz LA: Sjoblom T, Barad O, Bentwich Z, Szafranska AE, Labourier E et al: The colorectal microRNAome.

    Proc Natl Acad Sci U S A 2006, 103(10):3687-3692. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Wang X, Wang J, Ma H, Zhang J, Zhou X: Downregulation of miR-195 correlates with lymph node metastasis and poor prognosis in colorectal cancer.

    Med Oncol 2011, 29(2):919-27. PubMed Abstract | Publisher Full Text OpenURL

  41. Chou YT, Lin HH, Lien YC, Wang YH, Hong CF, Kao YR, Lin SC, Chang YC, Lin SY, Chen SJ, et al.: EGFR promotes lung tumorigenesis by activating miR-7 through a Ras/ERK/Myc pathway that targets the Ets2 transcriptional repressor ERF.

    Cancer Res 2010, 70(21):8822-8831. PubMed Abstract | Publisher Full Text OpenURL

  42. Foekens JA, Sieuwerts AM, Smid M, Look MP, de Weerd V, Boersma AW, Klijn JG, Wiemer EA, Martens JW: Four miRNAs associated with aggressiveness of lymph node-negative, estrogen receptor-positive human breast cancer.

    Proc Natl Acad Sci U S A 2008, 105(35):13021-13026. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Kefas B, Godlewski J, Comeau L, Li Y, Abounader R, Hawkinson M, Lee J, Fine H, Chiocca EA, Lawler S, et al.: microRNA-7 inhibits the epidermal growth factor receptor and the Akt pathway and is down-regulated in glioblastoma.

    Cancer Res 2008, 68(10):3566-3572. PubMed Abstract | Publisher Full Text OpenURL

  44. Webster RJ, Giles KM, Price KJ, Zhang PM, Mattick JS, Leedman PJ: Regulation of epidermal growth factor receptor signaling in human cancer cells by microRNA-7.

    J Biol Chem 2009, 284(9):5731-5741. PubMed Abstract | Publisher Full Text OpenURL

  45. Sugimura R, Li L: Noncanonical Wnt signaling in vertebrate development, stem cells, and diseases.

    Birth Defects Res C Embryo Today 2010, 90(4):243-256. PubMed Abstract | Publisher Full Text OpenURL

  46. Polakis P: Wnt signaling and cancer.

    Genes Dev 2000, 14(15):1837-1851. PubMed Abstract | Publisher Full Text OpenURL