Skip to main content

mRNA/microRNA gene expression profile in microsatellite unstable colorectal cancer



Colorectal cancer develops through two main genetic instability pathways characterized by distinct pathologic features and clinical outcome.


We investigated colon cancer samples (23 characterized by microsatellite stability, MSS, and 16 by high microsatellite instability, MSI-H) for genome-wide expression of microRNA (miRNA) and mRNA. Based on combined miRNA and mRNA gene expression, a molecular signature consisting of twenty seven differentially expressed genes, inclusive of 8 miRNAs, could correctly distinguish MSI-H versus MSS colon cancer samples. Among the differentially expressed miRNAs, various members of the oncogenic miR-17-92 family were significantly up-regulated in MSS cancers. The majority of protein coding genes were also up-regulated in MSS cancers. Their functional classification revealed that they were most frequently associated with cell cycle, DNA replication, recombination, repair, gastrointestinal disease and immune response.


This is the first report that indicates the existence of differences in miRNA expression between MSS versus MSI-H colorectal cancers. In addition, the work suggests that the combination of mRNA/miRNA expression signatures may represent a general approach for improving bio-molecular classification of human cancer.


Colorectal cancer develops through two main genetic pathways characterized by different forms of genomic instability [1]. Most tumors are generated by the chromosomal instability (CIN) pathway and display marked cytogenetic abnormalities, aneuploidy and allelic losses at multiple chromosomal arms. CIN is probably caused by various molecular mechanisms, but the underlying genetic alterations are still poorly defined. About 15% of colorectal carcinomas develop through the microsatellite instability (MSI) pathway. MSI tumors show stable karyotype, low frequencies of allelic losses and diploid nuclear DNA content. MSI results from defects in the DNA mismatch repair system (MMR) [2]. In HNPCC, MSI is produced by germline mutations of one of the MMR genes (MLH1, MSH2 and less frequently MSH6 and PMS2) with somatic inactivation of the wild-type allele [3, 4]. In sporadic tumors, MMR deficiency is near always determined by epigenetic inactivation of the MLH1 gene by biallelic promoter methylation [57]. MSI colorectal adenocarcinomas display distinctive pathologic features, such as proximal location, poor differentiation, frequent mucinous and medullary phenotype, and marked peritumoral and intratumoral lymphocytic infiltration [810]. MSI carcinomas have a more favorable clinical outcome than non-MSI tumors and the survival advantage conferred by the MSI phenotype is independent of tumor stage and other clinical and pathological variables [1113]. In addition, MMR-deficient cancer cells are thought to be less responsive to 5-fluorouracil and other chemotherapeutic agents in vitro and in vivo [1416].

Gene expression analysis based on genome-wide microarrays has been largely used to characterize human cancers. This approach allowed the identification of genes important in tumorigenesis. Furthermore, the discovery of gene expression signatures characteristic of distinctive clinico-pathological features suggested that expression profiles could be used for molecular classification of human cancer [1720]. Microarray tools have been recently enriched by the development of platforms for the analysis of microRNAs (miRNA) expression [21, 22]. miRNAs are a class of small non-coding RNAs involved in temporal and tissue-specific eukaryotic gene regulation [23]. Comparison between human cancers and their normal counterparts revealed that miRNAs exhibit differential expression profiles in normal versus cancer tissues [2429]. These studies revealed that some human miRNAs are consistently deregulated in human cancer, suggesting a role in tumorigenesis either as oncogenes or tumor suppressor genes [3032]. Unique miRNA expression signatures were found to be associated with bio-molecular and prognostic characteristics of human lung cancer and chronic lymphocytic leukemia [24, 33], indicating that miRNA signatures could be used to define biological or clinical features of human cancers. Known function of mammalian miRNAs is to post-transcriptionally regulate target mRNAs, implying that the combination of miRNAs and mRNAs expression may better represent the transcriptional program that dictates normal and tumor cell characteristics. Here, we identified differentially expressed miRNAs and mRNAs able to distinguish colon cancers with or without microsatellite instability.


We analyzed miRNA and mRNA expression using microarrays in a set of microsatellite stable (MSS) and unstable (MSI-H) colorectal cancers (23 MSS and 16 MSI-H) (Table 1) with the aim of recognizing the most significant differences in gene expression. In particular, this is the first study that investigates the differences in miRNAs expression between MSS and MSI tumors.

Table 1 Clinical and bio-pathological features of colorectal carcinomas employed in the study.

We carried out initial studies on miRNA and mRNA-chip data separately. The analysis of miRNA expression profiles of MSI-H versus MSS tumors identified 14 differentially expressed miRNAs (p < 0.05). Additional file 1 provides the list of deregulated miRNAs and their mean expression values and standard errors in the two groups. The analysis of mRNA expression profiles of MSI-H versus MSS identified a large number of differentially expressed genes (451 genes at p < 0.05; Additional file 2). By using the more stringent Bonferroni multiple testing correction at p < 0.05, a sub-list of 72 differentially expressed mRNAs was identified (Additional file 3).

To visually assess the ability of the differentially expressed miRNAs and mRNAs to distinguish MSS versus MSI-H tumors, we applied the combination of the 86 differentially expressed miRNAs and mRNAs to perform a supervised cluster analysis. Using this combined set of genes, a perfect classification of the tumor samples resulted (Figure 1). The combination of the two sets of genes resulted also in an improvement in the separation of the MSS versus the MSI tumor samples in comparison to the individual sets (see Additional files 4 and 5).

Figure 1
figure 1

Supervised hierarchical clustering of 39 colon cancer samples (23 MSS and 16 MSI-H) using expression data from the list of 86 differentially expressed mRNA/miRNA genes described in Supplemental Tables 1 and 3. Columns display the clustering of tumor samples; rows the clustering of genes. MSI status is indicated by colored squares: red for MSS and yellow for MSI-H tumor samples. All sample were properly assigned to the correct class. The expression intensity of each gene in each sample varies from red to green: red color means an expression value over the average across samples, green color the opposite. For quantification of color codes, see color-bar on the right.

A smaller selected signature of best predictors of microsatellite status was generated by applying Support Vector Machine (SVM) (GeneSpring software) and Prediction Analysis of Microarrays (PAM) algorithms to the 86 differentially expressed miRNAs and mRNAs: 27 genes (inclusive of mRNAs and miRNAs) were identified as predictors by at least one of the two algorithms (Table 2). Cluster analysis with these 27 miRNA/mRNAs still produced a perfect separation between the two tumor classes, and a PAM plot displayed a correct predictions for all the samples (Figure 2). The prediction of MSI status obtained with 72 mRNAs or 14 miRNAs signatures did not perform equally well when tested independently. Of the 39 samples, the miRNA-only signature could correctly predict 36 samples using SVM or 33 using PAM algorithms, and the mRNA-only signature could correctly predict 38 samples using SVM or 39 using PAM.

Table 2 Predictors of microsatellite status identified by Support Vector Machine and PAM algorithms.
Figure 2
figure 2

Prediction of microsatellite instability. (A) Prediction of MSI status for 39 samples using the PAM cross-validation procedure. The graphical representation shows the probabilities (0.0 to 1.0) of each tumor for being MSS (pink dot) or MSI-H (blue dot). All tumors were correctly assigned. (B) Supervised hierarchical clustering of colon cancer samples using expression data from the 27 predictor mRNA/miRNA genes described in Table 2. Using the predictors for the cluster analysis, a good separation among colorectal cancers in MSS (red) and MSI-H (yellow) groups was achieved. Expression intensity is shown as described in Figure 1.

Confirmation of the differential expression for some of the signature genes was achieved by Northern blot to validate miRNA expression and Real time RT-PCR to validate expression of mRNAs. Figure 3A shows the expression of miR-25 and miR-92 as detected by Northern blot. This analysis confirmed that the expression of these miRNAs is increased in MSS samples as detected by microarray analysis. Similarly, the evaluation of expression of 5 protein-coding genes showed an overlap with the expression values obtained by microarray analysis (Figure 3B).

Figure 3
figure 3

Confirmation of microarray data. (A) Northern blot analysis of miRNAs miR-25 and miR-92. As expected from microarray data, Northern blot analysis confirmed that these miRNAs were up-regulated in MSS samples. miRNA/U6 ratio of densitometric intensities is shown below each miRNA lane. (B) Quantitative Real time RT-PCR validation of microarray results. Microarray expression values of a selected panel of 5 genes differentially expressed between MSS and MSI-H groups of colorectal cancers were validated with quantitative real-time PCR. Values for microarray data were normalized on chip 50th percentile while for real-time PCR values are equal to ΔCt*10^6.


We report the analysis of the combined miRNA/mRNA expression for the discrimination of MSI-H versus MSS human colon cancer. Because colorectal tumors characterized by MSI are distinct from MSS tumors in many molecular aspects, such the association with the methylator phenotype, which is responsible for MLH1 methylation, the higher frequency of BRAF mutations and the lower frequencies of KRAS, APC and TP53 mutations, MSI and MSS colon cancers represent tumors with a different molecular background. Thus, it is reasonable that their overall gene expression pattern (including both mRNAs and microRNAs) might be affected by any of the above mentioned mechanisms. Indeed, our analyses identified 14 miRNAs and 451 mRNAs differentially expressed between the two genetically distinct colon cancer classes. These results not only indicate the existence of a mRNA/miRNA gene expression profile able to distinguish MSS versus MSI colon cancers, but also indicate that the combination of miRNA/mRNA gene expression data can distinguish microsatellite instability status of human colorectal cancers better than each RNA class independently, suggesting that the combination of mRNA and miRNA expression may potentially represent a general approach for improving characterization and classification of bio-molecular and possibly clinical traits associated with human cancer.

We also examined the prediction power of published lists of genes differentially expressed between MSI-H and MSS colorectal cancers [3436] against our data set. Two published lists describe differentially expressed genes (100 genes in the report from Di Pietro et al. 2005 [36] and 542 genes from Banerjea et al. 2004 [35]) and one is a short list of predictors (9 genes in Kruhoffer et al. 2005 [34]). Although these studies were all performed on Affymetrix platforms, they generated only partially overlapping results. The list of genes from Di Pietro et al. [36] was the one that better correlated with our results. Twenty two of their 100 genes were present in our list of differentially expressed genes and, for all the remaining genes, expression tendency completely matched our data. A cluster analysis of our data set using the 100 genes from the Di Pietro's paper could generate a good, although not perfect, separation between MSS versus MSI-H tumors and performed well in the prediction of our samples (36 of 39 samples were correctly classified using the SVM prediction) (data not shown). These results indicate that a similar set of genes emerged as differentially expressed between MSS and MSI-H colon cancers in the Di Pietro's and our studies, with discrepancies possibly due to small statistical differences. Discrepancies between results obtained with different microarray platforms are not uncommon and are mainly determined by differences in microarray probe sequences used to detect mRNA transcripts [37] and in algorithms used for predictor genes identification. Hence, the parallel between Di Pietro's and our results is significant given that data were generated using a different set of samples and a different microarray platform, and suggests that the commonly identified set of genes may represent the most significant differences between MSI-H and MSS colon cancers (Additional file 6). Unfortunately, data sets from published reports were not available for cross confirmation of our list of predictive genes on published microarray data.

In this report, quantitative real time PCRs for differentially expressed mRNAs and Northern blot for differentially expressed miRNAs were performed to validate microarray expression data. These methods confirmed the differential expression detected by microarray-based expression methods. Additional evidence for the robustness of data was the finding that the MLH1 gene is among the most significantly (P < 0.001) down-regulated genes in MSI-H tumors versus the MSS cancer set (Additional file 2). Indeed, this gene represents a sort of internal control, since it is well-known that the loss of MLH1 function, which confers the microsatellite instability phenotype to tumor samples, is caused by transcriptional silencing due to promoter methylation in MSI tumors [5, 6].

In addition to their role as discriminating markers of MSS versus MSI-H tumors, some of the gene-associated functions may possibly be involved in the different phenotypes that characterize the two types of colon cancers. In fact, analysis of the functions associated with 212 (for which annotation was present) of the 451 differentially expressed genes revealed that the most frequently associated classes were cell cycle, DNA replication, recombination, repair, gastrointestinal disease and immune response (Additional file 7), suggesting that these molecular differences may be responsible for traits that distinguish MSS and MSI-H tumors. Additional studies on the molecular and biological functions of these differentially expressed genes will be required to substantiate this hypothesis.

Among the differentially expressed genes, it is interesting to note the detection of the up-regulation of several members of the mir-17-92 family in the MSS colon cancers. This family includes fourteen homologous miRNAs organized in three gene clusters [38]. Our study revealed that, among these, miR-17-5p, miR-20, miR-25, miR-92-1, miR-92-2, miR-93-1 and miR-106a were significantly up-regulated in MSS versus MSI-H colon cancer. The chromosomes 13 and X gene clusters were previously found up-regulated in B-cell lymphoma [39]. It was also shown that c-MYC promotes their transcription [40] and, interestingly, enforced expression of the mir-17-92 cluster acted with c-myc to accelerate tumour development in a mouse B-cell lymphoma model [39]. In human solid tumors, the chromosome 13 mir-17-92 cluster was found up-regulated in small-cell lung cancer [41] and its ectopic over-expression enhanced lung cancer cell growth [41]. Since these data indicate that members of the mir-17-92 family can act as oncogenes to promote cell growth and inhibit apoptosis, our data suggests that up-regulation of these miRNAs may have a role in the more aggressive clinical behavior of MSS versus MSI-H tumors.


This report provides the first study on microRNA expression in MSI versus MSS colon cancer. We identified microRNAs that are differentially expressed between these two classes of tumors; moreover, the addition of microRNAs in the molecular classifier improves the separation between MSI and MSS cancer samples, suggesting that the mRNA/miRNA combination could provide an improved stratification of tumor-associated characters. Interestingly, the most prominent class of differentially expressed miRNAs includes various members of the oncogenic miR-17-92 family, suggesting that these microRNAs have a role in bio-pathologic characteristics that distinguish MSS versus MSI colon cancers.


Colorectal cancer samples

Samples of colorectal cancer tissue and matched normal colonic mucosa were obtained from fresh surgical specimens, frozen in liquid nitrogen, and stored at -80°C. Thirty nine carcinomas, 23 microsatellite stable (MSS) and 16 with high-frequency MSI (MSI-H) were analyzed in the course of the study. The clinico-pathological features of the tumors are detailed in Table 1. MSI status was determined with a fluorescence-based PCR method, using the five markers of the panel of Bethesda (BAT25, BAT26, D2S123, D5S346, and D17S250) [42]. PCR products were run in an ABI PRISM 377 DNA sequencer (Perkin-Elmer Applied Biosystems Division, Foster City, CA) and analyzed by the GeneScan 3.1 version software (Perkin-Elmer) [43]. According to the guidelines of the Workshop of Bethesda [42], tumors showing instability at two or more loci were classified as MSI-H and tumors without instability at any locus as MSS. None of the tumors included in this study exhibited instability at a single locus (low-frequency MSI or MSI-L). All MSI-H carcinomas displayed instability at mononucleotide DNA sequences (BAT25 and BAT26 markers). Tumors were also examined for expression of the DNA mismatch repair proteins MLH1 and MSH2, using the immunohistochemical analytic procedure previously described [43]. Carcinomas showing complete loss of nuclear MLH1 or MSH2 expression were classified as MLH1 or MSH2 negative, whereas cancers with normal expression of MLH1 and MSH2 gene products were classified as MLH1 and MSH2 positive. As reported in Table 1, all MSS carcinomas demonstrated normal nuclear expression of both MLH1 and MSH2 proteins (MLH1/MSH2 positive). By contrast, of the 16 MSI-H tumors 15 were MLH1 negative and one was MSH2 negative. As expected, MSI-H carcinomas were more frequently located in the proximal colon (P < 0.001) and poorly differentiated (P < 0.001) and more often contained a mucinous component (P < 0.001) with respect to MSS tumors (Table 1). Total RNA was isolated using Trizol (Invitrogen) following manufacturer's indications. RNA integrity was assessed on a 2100 Agilent Bioanalyzer. Only samples with intact RNA were used for the gene expression analysis.

RNA labeling and microarray hybridization

RNA labeling and hybridization on microRNA microarray chips was performed as previously described. [21, 22] Briefly, 5 μg of total RNA from each sample was biotin-labeled by reverse transcription using 5' biotin end labeled random examers oligo primer. Hybridization of biotin-labeled cDNA was carried out on our miRNA microarray chip (ArrayExpress accession number: A-MEXP-86), which contains 230 human microRNA probes, in triplicate. Hybridization signals were detected by biotin binding of a Streptavidin – Alexa 647 conjugate using Axon Scanner 4000B (Axon Instrument Inc. CA).

For mRNA-chip hybridization, 5 ug of total RNA were used for the production of biotinylated cRNA. Labeled cRNA was synthesized following the Ambion cRNA Synthesis Protocol and purified using an RNeasy® Kit (Qiagen). cRNA yield was quantified by measuring the UV absorbance at 260 nm. Unfragmented biotinylated cRNA (about 50 ug) was used for hybridization on custom made high density oligonucleotide arrays. The array was the human 18.5K Expression Bioarray (Compugen Human Oligo Set 1.0,, which contains 18,861 probes corresponding to approximately 17,260 unique human gene clusters and 18 bacterial control probes. All probes on these microarrays are 65-mer oligonucleotides spotted by contacting technologies and covalently attached to a polymeric matrix. Microarrays were hybridized in 6 × SSPE/40% formamide using a Tecan HS4800 hybridization station at 37°C for 20 hours, post-hybridization washed in 46°C pre-warmed 0.75 × TNT (1 × TNT:0.1 M Tris-HCL, pH7.5/0.15 M NaCL/0.05% Tween-20) at 46°C for one hour, and processed using a direct detection method of the biotin-containing transcripts by a Streptavidin-Alexa647 conjugate (1:500 diluted) in TNB (0.1 M Tris-HCL, pH7.5/0.15 M NaCL/0.5% Blocking Reagent-PerkinElmer) at RT for 30 min. Post-staining washing in 1 × TNT for one hour. Processed slides were scanned using Axon 4000B scanner (Molecular Device, CA). Image were quantified by GenePix Pro 6.0 software.

Microarray data analysis

MiRNA-chip and mRNA chip raw data were normalized separately using the GeneSpring software version 7.2 (Agilent). Both on chip and on gene median methods were used to normalize gene expression data. Microarray data were then joined into one GeneSpring genome and samples were assigned to one of two groups. The comparative analysis between MSI-H and MSS samples was carried out using the Welch's t-test and the Benjamini & Hochberg or Bonferroni (for a more stringent analysis) False Discovery Rate correction. Cluster analysis was performed using the Pearson correlation as a measure of similarity. Predictions were made using both Support Vector Machine algorithm and PAM software [44]. The Gene Ontology analysis of the gene lists of interest was generated by using the web-delivered tools of Ingenuity Pathway Analysis. Data have been submitted to Array Express (Accession number E-MEXP-326).

Northern blot of miRNAs

RNA samples (10 μg each) were electrophoresed on 15% acrylamide, 7 M urea Criterion precasted gels (Bio-Rad, Hercules, CA) and transferred onto Hybond N+ membrane (Amersham Bioscience, Piscataway, NJ). Membranes were hybridized with oligonucleotide probes, corresponding to the complementary sequences of the following mature miRNAs: miR-25 (TCA GAC CGA GAC AAG TGC AAT G) and miR-92 (CAG GCC GGG ACA AGT GCA ATA). Probes were 5'-end labeled using the polynucleotide kinase in the presence of P32-gamma-ATP. Hybridization was performed at 37°C in 7% SDS/0.2 mol/L Na2PO4 (pH 7.0) for 16 hours. Membranes were washed at 42°C, twice with 2 × standard saline phosphate (0.18 mol/L NaCl/10 mmol/L phosphate pH 7.4), 1 mmol/L EDTA (saline-sodium phosphate-EDTA, SSPE), and 0.1% SDS and twice with 0.5 × SSPE/0.1% SDS. Northern blots were re-hybridized after stripping the oligonucleotides used as probes in boiling 0.1% SDS for 10 minutes. As a control for normalization of RNA expression levels, we hybridized blots with an oligonucleotide probe complementary to the U6 RNA (5'-GCA GGG GCC ATG CTA ATC TTC TCT GTA TCG-3').

Quantitative Real time PCR for mRNAs

The RT reaction was perfomed using 200 ng of total RNA for each sample according to the manufacturer instructions (SuperScript First Strand Invitrogen). The Real Time reactions were performed using Taqman Gene Expression Assay on a ABI PRISM 7900HT. For each gene one set of primers and a probe were chosen from the Applied Biosystems list of TaqMan® Gene Expression Assays. Hs02379661_g1 for MTA2, Hs00228732_m1 for CTCL tumor antigen, Hs00228336_m1 for C13orf18, Hs00196125_m1 for VAV3 oncogene, Hs00914163_m1 for GGH, Hs00214886_m1 for FLJ20315. Expression analysis was performed in triplicate for each sample. Expression of 18S rRNA, which displayed the most constant expression among tested housekeeping genes between different samples, was used as endogenous reference control. The fold difference for each sample was obtained using the following equation 2-dCt. Ct is the Threshold Cycle, the cycle number at which the fluorescence generated within a reaction crosses the threshold; dCt = Ct average sample gene - Ct average 18S.


  1. Kinzler KW, Vogelstein B: Lessons from hereditary colorectal cancer. Cell. 1996, 87: 159-170. 10.1016/S0092-8674(00)81333-1

    Article  CAS  PubMed  Google Scholar 

  2. Peltomaki P: Role of DNA mismatch repair defects in the pathogenesis of human cancer. J Clin Oncol. 2003, 21: 1174-1179. 10.1200/JCO.2003.04.060

    Article  CAS  PubMed  Google Scholar 

  3. Lynch HT, de la Chapelle A: Hereditary colorectal cancer. N Engl J Med. 2003, 348: 919-932. 10.1056/NEJMra012242

    Article  CAS  PubMed  Google Scholar 

  4. de la Chapelle A: Genetic predisposition to colorectal cancer. Nat Rev Cancer. 2004, 4: 769-780. 10.1038/nrc1453

    Article  CAS  PubMed  Google Scholar 

  5. Kane MF, Loda M, Gaida GM, Lipman J, Mishra R, Goldman H, Jessup JM, Kolodner R: Methylation of the hMLH1 promoter correlates with lack of expression of hMLH1 in sporadic colon tumors and mismatch repair-defective human tumor cell lines. Cancer Res. 1997, 57: 808-811.

    CAS  PubMed  Google Scholar 

  6. Esteller M, Levine R, Baylin SB, Ellenson LH, Herman JG: MLH1 promoter hypermethylation is associated with the microsatellite instability phenotype in sporadic endometrial carcinomas. Oncogene. 1998, 17: 2413-2417. 10.1038/sj.onc.1202178

    Article  CAS  PubMed  Google Scholar 

  7. Herman JG, Umar A, Polyak K, Graff JR, Ahuja N, Issa JP, Markowitz S, Willson JK, Hamilton SR, Kinzler KW, Kane MF, Kolodner RD, Vogelstein B, Kunkel TA, Baylin SB: Incidence and functional consequences of hMLH1 promoter hypermethylation in colorectal carcinoma. Proc Natl Acad Sci U S A. 1998, 95: 6870-6875. 10.1073/pnas.95.12.6870

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Kim H, Jen J, Vogelstein B, Hamilton SR: Clinical and pathological characteristics of sporadic colorectal carcinomas with DNA replication errors in microsatellite sequences. Am J Pathol. 1994, 145: 148-156.

    PubMed Central  CAS  PubMed  Google Scholar 

  9. Jass JR, Do KA, Simms LA, Iino H, Wynter C, Pillay SP, Searle J, Radford-Smith G, Young J, Leggett B: Morphology of sporadic colorectal cancer with DNA replication errors. Gut. 1998, 42: 673-679.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Gafa R, Maestri I, Matteuzzi M, Santini A, Ferretti S, Cavazzini L, Lanza G: Sporadic colorectal adenocarcinomas with high-frequency microsatellite instability. Cancer. 2000, 89: 2025-2037. 10.1002/1097-0142(20001115)89:10<2025::AID-CNCR1>3.0.CO;2-S

    Article  CAS  PubMed  Google Scholar 

  11. Gryfe R, Kim H, Hsieh ET, Aronson MD, Holowaty EJ, Bull SB, Redston M, Gallinger S: Tumor microsatellite instability and clinical outcome in young patients with colorectal cancer. N Engl J Med. 2000, 342: 69-77. 10.1056/NEJM200001133420201

    Article  CAS  PubMed  Google Scholar 

  12. Popat S, Hubner R, Houlston RS: Systematic review of microsatellite instability and colorectal cancer prognosis. J Clin Oncol. 2005, 23: 609-618. 10.1200/JCO.2005.01.086

    Article  CAS  PubMed  Google Scholar 

  13. Lanza G, Gafa R, Santini A, Maestri I, Guerzoni L, Cavazzini L: Immunohistochemical test for MLH1 and MSH2 expression predicts clinical outcome in stage II and III colorectal cancer patients. J Clin Oncol. 2006, 24: 2359-2367. 10.1200/JCO.2005.03.2433

    Article  CAS  PubMed  Google Scholar 

  14. Arnold CN, Goel A, Boland CR: Role of hMLH1 promoter hypermethylation in drug resistance to 5-fluorouracil in colorectal cancer cell lines. Int J Cancer. 2003, 106: 66-73. 10.1002/ijc.11176

    Article  CAS  PubMed  Google Scholar 

  15. Carethers JM, Chauhan DP, Fink D, Nebel S, Bresalier RS, Howell SB, Boland CR: Mismatch repair proficiency and in vitro response to 5-fluorouracil. Gastroenterology. 1999, 117: 123-131. 10.1016/S0016-5085(99)70558-5

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Ribic CM, Sargent DJ, Moore MJ, Thibodeau SN, French AJ, Goldberg RM, Hamilton SR, Laurent-Puig P, Gryfe R, Shepherd LE, Tu D, Redston M, Gallinger S: Tumor microsatellite-instability status as a predictor of benefit from fluorouracil-based adjuvant chemotherapy for colon cancer. N Engl J Med. 2003, 349: 247-257. 10.1056/NEJMoa022289

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Brennan DJ, O'Brien SL, Fagan A, Culhane AC, Higgins DG, Duffy MJ, Gallagher WM: Application of DNA microarray technology in determining breast cancer prognosis and therapeutic response. Expert Opin Biol Ther. 2005, 5: 1069-1083. 10.1517/14712598.5.8.1069

    Article  CAS  PubMed  Google Scholar 

  18. Mischel PS, Cloughesy TF, Nelson SF: DNA-microarray analysis of brain cancer: molecular classification for therapy. Nat Rev Neurosci. 2004, 5: 782-792. 10.1038/nrn1518

    Article  CAS  PubMed  Google Scholar 

  19. Calvo A, Gonzalez-Moreno O, Yoon CY, Huh JI, Desai K, Nguyen QT, Green JE: Prostate cancer and the genomic revolution: Advances using microarray analyses. Mutat Res. 2005, 576: 66-79.

    Article  CAS  PubMed  Google Scholar 

  20. Dyrskjot L: Classification of bladder cancer by microarray expression profiling: towards a general clinical use of microarrays in cancer diagnostics. Expert Rev Mol Diagn. 2003, 3: 635-647. 10.1586/14737159.3.5.635

    Article  CAS  PubMed  Google Scholar 

  21. Liu CG, Calin GA, Meloon B, Gamliel N, Sevignani C, Ferracin M, Dumitru CD, Shimizu M, Zupo S, Dono M, Alder H, Bullrich F, Negrini M, Croce CM: An oligonucleotide microchip for genome-wide microRNA profiling in human and mouse tissues. Proc Natl Acad Sci U S A. 2004, 101: 9740-9744. 10.1073/pnas.0403293101

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Calin GA, Liu CG, Sevignani C, Ferracin M, Felli N, Dumitru CD, Shimizu M, Cimmino A, Zupo S, Dono M, Dell'Aquila ML, Alder H, Rassenti L, Kipps TJ, Bullrich F, Negrini M, Croce CM: MicroRNA profiling reveals distinct signatures in B cell chronic lymphocytic leukemias. Proc Natl Acad Sci U S A. 2004, 101: 11755-11760. 10.1073/pnas.0404432101

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel genes coding for small expressed RNAs. Science. 2001, 294: 853-858. 10.1126/science.1064921

    Article  CAS  PubMed  Google Scholar 

  24. Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T, Calin GA, Liu CG, Croce CM, Harris CC: Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006, 9: 189-198. 10.1016/j.ccr.2006.01.025

    Article  CAS  PubMed  Google Scholar 

  25. Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, Prueitt RL, Yanaihara N, Lanza G, Scarpa A, Vecchione A, Negrini M, Harris CC, Croce CM: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006, 103: 2257-2261. 10.1073/pnas.0510565103

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702

    Article  CAS  PubMed  Google Scholar 

  27. Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6: 259-269. 10.1038/nrc1840

    Article  CAS  PubMed  Google Scholar 

  28. Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, Menard S, Palazzo JP, Rosenberg A, Musiani P, Volinia S, Nenci I, Calin GA, Querzoli P, Negrini M, Croce CM: MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005, 65: 7065-7070. 10.1158/0008-5472.CAN-05-1783

    Article  CAS  PubMed  Google Scholar 

  29. Croce CM, Calin GA: miRNAs, cancer, and stem cell division. Cell. 2005, 122: 6-7. 10.1016/j.cell.2005.06.036

    Article  CAS  PubMed  Google Scholar 

  30. Chen CZ: MicroRNAs as oncogenes and tumor suppressors. N Engl J Med. 2005, 353: 1768-1771. 10.1056/NEJMp058190

    Article  CAS  PubMed  Google Scholar 

  31. Cimmino A, Calin GA, Fabbri M, Iorio MV, Ferracin M, Shimizu M, Wojcik SE, Aqeilan RI, Zupo S, Dono M, Rassenti L, Alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: miR-15 and miR-16 induce apoptosis by targeting BCL2. Proc Natl Acad Sci U S A. 2005, 102: 13944-13949. 10.1073/pnas.0506654102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647. 10.1016/j.cell.2005.01.014

    Article  CAS  PubMed  Google Scholar 

  33. Calin GA, Ferracin M, Cimmino A, Di Leva G, Shimizu M, Wojcik SE, Iorio MV, Visone R, Sever NI, Fabbri M, Iuliano R, Palumbo T, Pichiorri F, Roldo C, Garzon R, Sevignani C, Rassenti L, Alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. N Engl J Med. 2005, 353: 1793-1801. 10.1056/NEJMoa050995

    Article  CAS  PubMed  Google Scholar 

  34. Kruhoffer M, Jensen JL, Laiho P, Dyrskjot L, Salovaara R, Arango D, Birkenkamp-Demtroder K, Sorensen FB, Christensen LL, Buhl L, Mecklin JP, Jarvinen H, Thykjaer T, Wikman FP, Bech-Knudsen F, Juhola M, Nupponen NN, Laurberg S, Andersen CL, Aaltonen LA, Orntoft TF: Gene expression signatures for colorectal cancer microsatellite status and HNPCC. Br J Cancer. 2005, 92: 2240-2248. 10.1038/sj.bjc.6602621

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Banerjea A, Ahmed S, Hands RE, Huang F, Han X, Shaw PM, Feakins R, Bustin SA, Dorudi S: Colorectal cancers with microsatellite instability display mRNA expression signatures characteristic of increased immunogenicity. Mol Cancer. 2004, 3: 21. 10.1186/1476-4598-3-21

    Article  PubMed Central  PubMed  Google Scholar 

  36. di Pietro M, Bellver JS, Menigatti M, Bannwart F, Schnider A, Russell A, Truninger K, Jiricny J, Marra G: Defective DNA mismatch repair determines a characteristic transcriptional profile in proximal colon cancers. Gastroenterology. 2005, 129: 1047-1059. 10.1053/j.gastro.2005.06.028

    Article  CAS  PubMed  Google Scholar 

  37. Carter SL, Eklund AC, Mecham BH, Kohane IS, Szallasi Z: Redefinition of Affymetrix probe sets by sequence overlap with cDNA microarray probes reduces cross-platform inconsistencies in cancer-associated gene expression measurements. BMC Bioinformatics. 2005, 6: 107. 10.1186/1471-2105-6-107

    Article  PubMed Central  PubMed  Google Scholar 

  38. Tanzer A, Stadler PF: Molecular evolution of a microRNA cluster. J Mol Biol. 2004, 339: 327-335. 10.1016/j.jmb.2004.03.065

    Article  CAS  PubMed  Google Scholar 

  39. He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, Hammond SM: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435: 828-833. 10.1038/nature03552

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435: 839-843. 10.1038/nature03677

    Article  PubMed  Google Scholar 

  41. Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T: A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res. 2005, 65: 9628-9632. 10.1158/0008-5472.CAN-05-2352

    Article  CAS  PubMed  Google Scholar 

  42. Boland CR, Thibodeau SN, Hamilton SR, Sidransky D, Eshleman JR, Burt RW, Meltzer SJ, Rodriguez-Bigas MA, Fodde R, Ranzani GN, Srivastava S: A National Cancer Institute Workshop on Microsatellite Instability for cancer detection and familial predisposition: development of international criteria for the determination of microsatellite instability in colorectal cancer. Cancer Res. 1998, 58: 5248-5257.

    CAS  PubMed  Google Scholar 

  43. Lanza G, Gafa R, Maestri I, Santini A, Matteuzzi M, Cavazzini L: Immunohistochemical pattern of MLH1/MSH2 expression is related to clinical and pathological features in colorectal adenocarcinomas with microsatellite instability. Mod Pathol. 2002, 15: 741-749. 10.1097/01.MP.0000018979.68686.B2

    Article  PubMed  Google Scholar 

  44. Tibshirani R, Hastie T, Narasimhan B, Chu G: Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci U S A. 2002, 99: 6567-6572. 10.1073/pnas.082099299

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


We wish to thank Augusto Bevilaqua, Pietro Zucchini, Annalisa Peverati and Iva Pivanti for the excellent technical support. This work was supported by funding from Associazione Italiana per la Ricerca sul Cancro (Regional grant), from Ministero della Salute (Program Ricerca Oncologica 2006) and from Ministero dell'Università e Ricerca to MN and GL, and by Program Project Grants from the National Cancer Institute to CMC and by a Kimmel Foundation Scholar award to GAC. MF is a recipient of a fellowship from Fondazione Italiana per la Ricerca sul Cancro.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Massimo Negrini.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

GL, RG, FP samples histo-pathological characterization; MF, AV, C-gL, GAC microarray preparation, hybridization and data analysis; RS, quantitative RT-PCR; CMC, MN, work planning and manuscript preparation. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: MicroRNAs differentially expressed between MSS and MSI-H colorectal cancers at P < 0.05. List of differentially expressed microRNAs (PDF 52 KB)


Additional file 2: Protein-coding genes differentially expressed between MSS and MSI-H colorectal cancers at P < 0.05. List of differentially expressed protein-coding genes (PDF 116 KB)


Additional file 3: Protein-coding genes differentially expressed between MSS and MSI-H colorectal cancers with Bonferroni's correction. List of differentially expressed protein-coding genes, selected from Additional file 2 using the more stringent Bonferroni's statistical correction (PDF 67 KB)


Additional file 4: Classification of tumors according to expression of microRNAs. Cluster analysis based on differentially expressed microRNAs (PDF 180 KB)


Additional file 5: Classification of tumors according to expression of 72 differentially expressed mRNAs. Cluster analysis based on differentially expressed mRNAs (PDF 278 KB)


Additional file 6: Differentially expressed protein-coding genes shared with Di Pietro's study. List of protein-coding genes shared by this study and Di Pietro's study (PDF 68 KB)


Additional file 7: Gene Ontology classification based on Ingenuity Pathway Analysis. List of functional classes present in the differentially expressed protein-coding genes (PDF 70 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lanza, G., Ferracin, M., Gafà, R. et al. mRNA/microRNA gene expression profile in microsatellite unstable colorectal cancer. Mol Cancer 6, 54 (2007).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: