- Open Access
Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues
Molecular Cancer volume 5, Article number: 29 (2006)
MicroRNAs (miRNAs) are short non-coding RNA molecules playing regulatory roles by repressing translation or cleaving RNA transcripts. Although the number of verified human miRNA is still expanding, only few have been functionally described. However, emerging evidences suggest the potential involvement of altered regulation of miRNA in pathogenesis of cancers and these genes are thought to function as both tumours suppressor and oncogenes.
In our study, we examined by Real-Time PCR the expression of 156 mature miRNA in colorectal cancer. The analysis by several bioinformatics algorithms of colorectal tumours and adjacent non-neoplastic tissues from patients and colorectal cancer cell lines allowed identifying a group of 13 miRNA whose expression is significantly altered in this tumor. The most significantly deregulated miRNA being miR-31, miR-96, miR-133b, miR-135b, miR-145, and miR-183. In addition, the expression level of miR-31 was correlated with the stage of CRC tumor.
Our results suggest that miRNA expression profile could have relevance to the biological and clinical behavior of colorectal neoplasia.
MicroRNAs (miRNAs) are 19- to 25-nt non coding RNAs that are cleaved from 70- to 100-nt hairpin-shaped precursors [1, 2]. Initial estimates, relaying mostly on evolutionary conservation, suggested there were up to 255 humans miRNAs. More recent analysis have demonstrated there are numerous non conserved humans miRNAs and suggest this number may be significantly larger. Although the precise biological are not yet fully understood, miRNAs seems to be crucial factors of diverse regulation pathways, including development, cell differentiation, proliferation and apoptosis [3–6]. Moreover, miss-regulation of miRNA expression might contribute to human disease [7–10].
A more recent link between miRNA function and cancer pathogenesis is supported by studies examining the expression of miRNA in clinical samples. Calin et al reported the first evidence and showed a down-regulation of miRNA-15 and miRNA-16 in a majority of chronic lymphatic leukemia (CLL) . Then, altered miRNA expression has been reported, in lung cancer , breast cancer , glioblastoma , hepatocellular carcinoma , papillary thyroid carcinoma  and more recently colorectal cancer . These results could indicate that miRNA may be a new class of genes involved in human oncogenesis.
Up until very recently, the most common method for quantifying miRNA was Northern blotting. Over the past year, a number of different approaches to quantify miRNAs have been described, including cDNA arrays [17, 18], a modified Invader assay , a bead-based flow cytometric assay  and Real-time PCR . Arrays, invader assay and bead-base miRNA expression do not amplify miRNA and thus the sensitivity is often compromised. The main advantage of real-time PCR is that is more quantitative and more sensitive that other high-throughput assays. However, it could be an important disadvantage if the number of miRNA increase as expected. In this case, Real-Time PCR will be less practical than microarrays.
In our study, we analyze by Real-Time PCR the expression of 156 mature miRNA in colorectal cancer (CRC) cell lines, tumoral and normal-paired tissues from clinical samples. CRC is one of the major causes of cancer death worldwide. At a molecular level, much progress has been made in the last two decades in the identification and characterization of the genetic changes involved in the malignant colorectal transformation process . A number of molecular studies have shown that colon carcinogenesis results from an accumulation of epigenetic and genetic alterations, including activating mutations of the K-ras proto-oncogene and inactivating mutations of APC and TP53 tumor suppressor genes or of DNA repair genes. However, this stepwise model of colorectal tumorigenesis has been mainly validated conceptually, and there is mounting evidence that alternative genetic events may occur during colorectal carcinogenesis, sometimes preferentially, sometimes randomly, and sometimes with an overlap. miRNA expression regulation could help to identify mRNA targets associated with different colorectal carcinogenesis pathways and their role as potential therapeutic targets.
In the present study, we examined by Real-time PCR the expression of 156 mature miRNA in a panel of 16 CRC cell lines and 12 matched-pair of tumoral and non-tumoral tissues from patients. We identified a subset of 13 miRNAs differentially expressed in CRC cell lines and clinical samples.
Results and discussion
miRNA expression in CRC cell lines
In order to investigate miRNA differential expression in human colorectal cancer, we analyzed by real-time PCR using TaqMan MicroRNA Assay kit (Applied Biosystems), the expression of 156 mature miRNAs in total RNA extracted from 15 CRC cell lines. We compared their miRNA expression profile with those of CCD-18Co (human normal colon cell line).
It is generally accepted that gene-expression levels should be normalized by a carefully selectable stable internal control gene. However, to validate the presumed stable expression of a given control gene, prior knowledge of a reliable measure to normalize this gene in order to remove any non specific variation is required. To address this problem we assessed the normalization data using three different approaches: let-7a (a miRNA that manufacturer suggests may be useful as an endogenous microRNA control according their preliminary data across several human tissues and cell lines), 18s rRNA (the most stable housekeeping gene in our CRC samples) and global median-normalization, similar to microarray analysis.
After normalization, data were transformed as log10 of relative quantity (RQ) of target miRNA relative to control sample. As shown Additional file 1, the different normalization approach reveals similar results. Analysis of k-means clustering (k = 3) identify a group of 22 and 22 miRNA homogeneously up-regulated and down-regulated respectively in all colorectal cancer cell lines and commonly detected with the three different normalization approach used. Figure 1 shows patterns of expression of these 44 miRNA after normalization with median-global normalization. Remarkably this classification only include those miRNA whose expression are most prominently altered and in addition the expression of this group of miRNAs is highly reproducible in all cell lines analyzed. Interestingly, clustering analysis divided CRC cell lines in two different groups. Analysis of different common genetic alterations described in colorectal cancer including activation of oncogenes (KRAS, BRAF) and inactivation of tumor supressor genes (TP53) and microsatellite instability status (MSI) showed that these groups could be differentiate according the presence of mutation in KRAS and BRAF genes. One group included DLD1, SW1116, SW620, SW480, HCT116, Lovo, Colo320, LS174, LS513 and LS411 CRC cell lines. All of them, except for LS411 and Colo320, harbor mutation in KRAS gene. On the other hand, the other group includes mainly CRC cell lines with BRAF mutation (WiDR, SW1417, Caco2 and RKO). SAM analysis between both groups identified 6 miRNA differentially expressed. Colorectal cancer cell lines with KRAS mutations showed an over-expression of miR-9, miR-9*, miR-95, miR-148a, miR-190 and miR-372, in relation to the human normal colon cell line. This over-expression was lower in whose colorectal cancer cell lines with mutations in BRAF. The presences of both mutations were mutually excluding. It is interesting to note that the predicted miRNA for BRAF regulation (using miRANDA, TargetScan and PicTar algorithms) included miR-9. This miRNA was just over-expressed in CRC cell lines with BRAF wild-type. Moreover, miR-372 has been recently described as potential oncogene that collaborate with oncogenic RAS in cellular transformation .
In human colorectal cancers, KRAS mutation has been considered an early event in the development of adenomas . This genetic event is more common in large adenomas than small ones, suggesting that it may be required for the activation of adenoma progression. Recently, the activation of BRAF has been reported to occur by somatic mutation in many human cancers, particularly in human malignant melanoma (over 60%) , human colorectal cancers (5–15%)  and a small fraction of other cancers [27, 28]. The majority of the BRAF mutations each represent a single nucleotide change of T-A at nucleotide 1796, resulting in the change of valine to glutamic acid at codon 599 within the activation segment of BRAF. Although BRAF mutations were found in about 5–15% of colorectal carcinomas, colorectal carcinomas with BRAF mutations tended to be in lower clinical tumor stages. However, it has been suggested that alteration in the BRAF gene may cause the activation of the RAS/RAF/MEK/ERK pathway , consequently increasing cell proliferation but suppressing the inhibition of apoptosis. Differential miRNA expression between CRC samples could help to identify different mechanisms of CRC carcinogenesis associate with alterations of the RAS/RAF/MEK/ERK pathway.
As shown table 1, the fold-change observed in CRC cell lines in relation to CDC18Co differed between -4.5 to -1.5 log10 for down-regulation and 1.4 and 3.8 log10 for up-regulation. Some of the genes encoding miRNA that are modulated in CRC cell lines are located in determined chromosome segments, suggesting that their tumor-specific expression could be due to DNA abnormalities. In this context, we observed a preferential down-regulation in region 14q32.31 including miRNA miR-127, miR-370, miR-299, miR-154, miR-154*, miR-323, miR-134, miR-368 and miR-337. By using a computer-assisted approach, Seitz et al.  have identified 46 potential miRNA located in human 14q32 domain, 40 of which are organized as a large cluster. Although some of these clustered miRNA genes appear to be encoded by a single-copy DNA sequence, most of then are arranged in tandem arrays of closely related sequences.
However, 14q it is not a region usually deleted in CRC cancers although their loss have been associated with disease progression and worse prognosis . On the opposite, we can hypothesize that differential expression could be regulated by modulation of its transcription. We think that this hypothesis may be supported by the observation that different "isoforms" of some down-regulated and up-regulated miRNAin CRC cell lines are located in different chromosomes, and their coordinated expression might reflect the existence of a common target. The expression of mir-200a, mir-200b and mir-200c, located in two different chromosomes (1 and 12) and with a high sequence-homology, are up-regulated in all CRC cell lines. The analysis of their putative targets showed MLH1 and MSH2 as two candidate genes whose transcription could be down-regulated by miRNA.
Our findings indicate that miRNA expression patterns are closely related to characteristics of tumor derived cell lines. These patterns may either mark specific biologic characteristics or may mediate specific biologic activities important for the pathobiology of malignant tumors.
miRNA expression in colorectal tumours and adjacent non-tumor tissues
In order to investigate whether miRNAs are differentially expressed in CRC versus normal colon tissues, we analyzed miRNA expression in 12 matched-pairs of tumoral and non-tumoral tissues. After testing three different approaches to normalize the Ct raw data in CRC cell lines, median-normalization was selected as method for clinical samples since normal distribution was not required. Meanwhile in our study in CRC cell lines no differences were found in let-7a expression between tumoral and normal cell line, recent evidences identify let7-family as differentially expressed in CRC  and lung cancer . Moreover, global median normalization could provide results more easily comparable with those already published with microarray technology.
To identify miRNA with significantly differential expression among CRC samples, two multivariate permutation test provided in BRB-ArrayTools were performed: Class Comparison between Groups of Arrays and SAM (Significance Analysis of Microarrays). In both cases we selected paired t-test options and a FDR (False Discovery Rate) less to 10%. Fifty-nine miRNAs were significant when Class Comparison test was applied, 68 miRNA were significant using SAM test and 53 miRNA are common in both test. As expected, fold-change observed in clinical samples is less homogeneously distributed among samples that already obtained in CRC cell lines. It is not surprising regarding that patients samples are composed of mixed populations, whereas cell lines are clearly more uniform.
Interestingly, our results in CRC samples are in agreement with recent data published in CRC using direct miRNA cloning and SAGE (miRAGE). In this context, we detected an over-expression of miR-19a, miR-21, miR-29a, miR-92, miR-148a, miR-200b, and a down-regulation of miR-30c, miR-133a and miR-145 (figure 2). Moreover, change of expression of some of these miRNA has been previously reported in lung and breast cancer, B-cell lymphomas, and glioblastoma. miR-19a and miR-20 are including in the cluster miR-17-92 and it is located at intron 3 of C13orf25. The transfection of C13orf25 in lung cancer cell line enhancing cell growth and the introduction of miR-17-92 into hematopoietic stem cells in Eu-myc transgenic mice accelerates the formation of lymphoid malignancies. Furthermore, miR-21 has been described as an antiapoptotic factor in human glioblastoma cell lines. In contrast, other authors report that miR-21 suppression increase growth in HeLa cells without affecting their apoptosis. The different biologic effects of any particular miRNA in different cells could be dependent of the cell-specific repertoire in target genes. Some of miRNA differentially expressed in CRC samples have been associated with clinical parameters in other cancers. In particular, miR-145 is progressively down-regulated from normal breast to cancer with high proliferation index and miR-21 is progressively up-regulated with high grade tumor stage.
To identify the smallest set of predictive miRNAs differentiating normal versus cancer tissues, we have used support vector machines (SVMs) techniques. We attempted to use the class prediction tool (BRB-Array tools) that creates a multivariate predictor for determining to which of the two classes a given sample belongs. Several multivariate classification methods are available, including the Compound Covariate Predictor, Diagonal Linear Discriminate Analysis, Nearest Neighbor Predictor, Nearest Centroid Predictor, and Support Vector Machine Predictor. The classifier is composed for 18 miRNA, 10 down-regulated and 8 up-regulated, all of them significantly different by Class Comparison and SAM tests.
When we compared expression of these miRNA in CRC cell lines, 5 of 18 miRNA were revealed in the k-means analysis as those of highest fold-changes (in relation to CDC18Co). However, Class Comparison analysis between 15 CRC cell lines and the 12 non-tumoral colon tissues identifies 13 miRNA altered in both systems, CRC patients samples and CRC cell lines (Table 2). These results could indicate that miRNA profile in CRC cell lines can not be used to infer miRNA expression in clinical samples meanwhile cell lines can be used as model to validate and performed functional assays of data obtained in clinical samples. In this sense, hierarchical clustering of expression of these 13 miRNAs in CRC patients samples and CRC cell lines clearly separate samples in two groups: in one branch, the most different of non-tumoral samples were included all the CRC cell lines, and in the other branch tumoral samples of patients (Figure 3).
The expression of 5 of 13 miRNA is already described altered in CRC, lung and breast cancer, glioblastoma, B-cell lymphoma and CLL. Among the differentially expressed miRNAs, miR-31, miR-96, miR-133b, miR-135b, miR-145 and miR-183 as the most consistently deregulated in CRC. Two of them, miR-133b and miR-145 were down-regulated and the remaining four, miR-31, miR-96, miR-135b and miR-183, were up-regulated, suggesting that they may potentially act as tumor suppressor genes or oncogenes, respectively.
miR-145 was identified as a specific miRNA down-regulated in colorectal neoplasia and analysis of their pre-miRNA indicate that this reduction is due to posttranscriptional process . Recently, Cummins et al obtained similar results in CRC  and down-regulation of miR-145 have also reported in lung  and breast cancer . In our study, expression of miR-145 was not detected in any of 15 CRC cell lines tested and down-regulation was detected in all tumor samples. Other important down-regulated miRNA in our study was miR-133b. In our knowledge, this miRNA has not previously identified deregulated in cancer. For both down-regulated miRNAs (miR-145 and miR-133b), it may be expected that potential targets could include oncogenes or genes encoding proteins with potential oncogenic functions. Indeed, among putative targets for miR-145 with potential oncogenic functions, Iorio et al  described MYCN, FOS, YES, and FLI, cell cycle promoters such as cyclins D2 and L1; and MAPK transduction proteins such as MAP3K3 and MAPK4K4. Among putative targets of miR-133b, the most notable oncogenic target is KRAS. Interestingly, the proto-oncogen YES1 and the transduction protein MAP3K3 were potential targets of both miR-145 and miR-133b.
For the up-regulated miRNAs, miR-135b, miR-31, miR-96 and miR-183, it may be expected that gene targets belong to the class of tumor suppressor genes. miR-96, miR-182 and miR-183 are located in the same chromosomal region, 7q32.2. miR-182 was not detected as preferentially over-expressed with the most restricted analysis, but their up-regulation was clearly observed in CRC cell lines analysis (table 1). CHES1 protein was identified as a potential target of both miR-96 and miR-182. CHES 1 is a member of the forkhead family of transcription factors that repress genes involved in apoptosis. Other members of this family, including FOXF2, FOXK2, FOXO1A, FOXO3A and FOXQ 1, were also found as putative targets of miR-182, miR-183 and miR-96.
Finally, our analysis of a small number of CRC samples compared miRNA expression in tumors according to pathological stage (stage II versus stage IV). The up-regulation of miR-31 was significantly higher in stage IV than in CRC samples stage II (p = 0.028) (Figure 4). The expression levels of miR-31 were higher in the tumor samples and CRC cell lines in comparison to the non-tumoral samples and was related to pathological stage, suggesting that this miRNA could contribute to both, the tumorogenesis and the acquisition of a more aggressive phenotype in CRC. Other members of the forkhead family transcription factors, such as FOXC2 and FOXP3, were identified as putative targets of miR-31. Future studies will determine the correlation between of these miRNAs and their host genes in CRC.
In summary, our results by Real-time PCR identify alterations of miRNA expression in CRC that may deregulate cancer-related genes and would provide potential mechanisms that underly in the carcinogenesis and further acquisition of a more aggressive phenotype in colon cancer.
Materials and methods
Cell lines and tissues
The following 16 cancer cell lines were used: CDC18Co (human normal colon), HCT15, RKO, DLD1, Lovo, LS411, SW1417, Caco2, LS513, SW1116, HCT116, SW480, Colo320, SW620, WiDR and LS174. Colorectal cell lines were cultures in a humidified atmosphere of 95% air, 5% CO2 using recommended medium and 10% FBS. Colorectal tumours and paired-adjacent normal colorectal tissues were received from our Institutional Bank of Tumors. Collection and distribution of colorectal tissues were approved by the appropriate Institution Review Board.
RNA extraction, reverse transcription and Real-Time PCR quantification
Total RNA was extracted from cells with a cell density of 75% confluent using Trizol® total RNA isolation reagent (Gibco BRL, Life Technologies, Gaitherburg, MD, USA) as per the manufacturer's protocol. Total RNA was isolated from the frozen tissues disrupted using an Ultra Turrax T25 homogenizer and using Trizol. The concentration was quantified using NanoDrop Specthophotometer (NanoDrop Technologies, USA).
cDNA was synthesized from total RNA using gene-specific primers according to the TaqMan MicroRNA Assay protocol (PE Applied Biosystems, Foster City, CA). Reverse transcriptase reactions contained 10 ng of RNA samples, 50 nM stem-loop RT primer, 1 × RT buffer, 0.25 mM each of dNTPs, 3.33 U/μl MultiScribe reverse transcriptase and 0.25 U/μl RNase Inhibitor (all purchased from cDNA Archive kit of Applied Biosystems). The 7.5 μl reactions were incubated in an Applied Biosystems 9800 ThermaCycler in a 96-well plate for 30 min at 16°C, 30 min at 42°C, 5 min at 85°C and then held at 4°C.
Real-time PCR was performed using an Applied Biosystems 7300 Sequence Detection system. The 10 μl PCR included 0.67 μl RT product, 1× TaqMan Universal PCR master mix and 1 μl of primers and probe mix of the TaqMan MicroRNA Assay protocol (PE Applied Biosystems). The reactions were incubated in a 96-well optical plate at 95°C for 10 min, followed by 40 cycles of 95°C for 15s and 60° for 10 min. The Ct data was determinate using default threshold settings. The threshold cycle (Ct) is defined as the fractional cycle number at which the fluorescence passes the fixed threshold.
Normalization and data analysis
Careful normalization is essential for the accurate quantification of mRNA levels and commonly, normalization of the target gene with an endogenous standard, mainly housekeeping genes, is applied. However for miRNA, there is not data about the expression of miRNA as normalization control.
In our study, we have tried different approach for normalization data. First, miRNA expression data was normalized to let7-a miRNA (according to the manufacturer's suggestions). Relative quantification of miRNA expression was calculated with the 2-ΔΔCt method (Applied Biosystems User Bulletin N°2 (P/N 4303859)). The data were presented as log10 of relative quantity (RQ) of target miRNA, normalized respect to miR-let-7a and relative to a calibrator sample. As calibrator, we are used for colorectal cancer cell lines CDC18Co (human normal colon) and for CRC samples the paired-normal tissues.
The TaqMan Human Endogenous Control Plate (Applied Biosystems) precoated with lyophilized primers and probes for 11 differently commonly human control genes was also used to asses gene expression in two tumours and two normal colon tissues. PCR was set up according to the manufacturer's instructions. Among 11 genes analyzed, the variability in expression of 18s rRNA was shown to be the lowest (data not shown). Therefore, our second normalization approach showed the data as log10 of relative quantity (RQ) of target miRNA, normalized to 18s rRNA and relative to control sample.
Finally, similar to microarray data, raw data Ct was normalized and analyzed using BRB ArrayTools version 3.3.2. (Richard Simon and Amy Peng Lam, National Cancer Institute, Bethesda). After global median normalization, normalized data were presented as log10 of relative quantity (RQ) of target miRNA relative to a control sample. Class Comparison and Significant analysis of microarrays (SAM) was performed to identify differentially expressed miRNAs between tumors and normal samples.
Visualization of results was performed with the different normalized data using average linkage and Euclidean distance as a measurement of similarity using GENESIS Software (Alexander Sturn, Institute for Genomics and Bioinformatics, Graz University of Technology).
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5
He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004, 5 (7): 522-531. 10.1038/nrg1379
Miska EA: How microRNAs control cell division, differentiation and death. Curr Opin Genet Dev. 2005, 15 (5): 563-568. 10.1016/j.gde.2005.08.005
O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435 (7043): 839-843. 10.1038/nature03677
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 (5): 635-647. 10.1016/j.cell.2005.01.014
Ambros V: The functions of animal microRNAs. Nature. 2004, 431 (7006): 350-355. 10.1038/nature02871
Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6 (4): 259-269. 10.1038/nrc1840
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 (3): 189-198. 10.1016/j.ccr.2006.01.025
Cummins JM, He Y, Leary RJ, Pagliarini R, Diaz LAJ, Sjoblom T, Barad O, Bentwich Z, Szafranska AE, Labourier E, Raymond CK, Roberts BS, Juhl H, Kinzler KW, Vogelstein B, Velculescu VE: The colorectal microRNAome. Proc Natl Acad Sci U S A. 2006, 103 (10): 3687-3692. 10.1073/pnas.0511155103
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 (7): 2257-2261. 10.1073/pnas.0510565103
Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, Aldler H, Rattan S, Keating M, Rai K, Rassenti L, Kipps T, Negrini M, Bullrich F, Croce CM: Frequent deletions and down-regulation of micro- RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. 2002, 99 (24): 15524-15529. 10.1073/pnas.242606799
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 (21): 9628-9632. 10.1158/0008-5472.CAN-05-2352
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 (16): 7065-7070. 10.1158/0008-5472.CAN-05-1783
Chan JA, Krichevsky AM, Kosik KS: MicroRNA-21 is an antiapoptotic factor in human glioblastoma cells. Cancer Res. 2005, 65 (14): 6029-6033. 10.1158/0008-5472.CAN-05-0137
Murakami Y, Yasuda T, Saigo K, Urashima T, Toyoda H, Okanoue T, Shimotohno K: Comprehensive analysis of microRNA expression patterns in hepatocellular carcinoma and non-tumorous tissues. Oncogene. 2005.
He H, Jazdzewski K, Li W, Liyanarachchi S, Nagy R, Volinia S, Calin GA, Liu CG, Franssila K, Suster S, Kloos RT, Croce CM, de la Chapelle A: The role of microRNA genes in papillary thyroid carcinoma. Proc Natl Acad Sci U S A. 2005, 102 (52): 19075-19080. 10.1073/pnas.0509603102
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 (26): 9740-9744. 10.1073/pnas.0403293101
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, Bartel DP, Linsley PS, Johnson JM: Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005, 433 (7027): 769-773. 10.1038/nature03315
Allawi HT, Dahlberg JE, Olson S, Lund E, Olson M, Ma WP, Takova T, Neri BP, Lyamichev VI: Quantitation of microRNAs using a modified Invader assay. Rna. 2004, 10 (7): 1153-1161. 10.1261/rna.5250604
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 (7043): 834-838. 10.1038/nature03702
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33 (20): e179- 10.1093/nar/gni178
Fearon ER, Vogelstein B: A genetic model for colorectal tumorigenesis. Cell. 1990, 61 (5): 759-767. 10.1016/0092-8674(90)90186-I
Voorhoeve PM, le Sage C, Schrier M, Gillis AJ, Stoop H, Nagel R, Liu YP, van Duijse J, Drost J, Griekspoor A, Zlotorynski E, Yabuta N, De Vita G, Nojima H, Looijenga LH, Agami R: A Genetic Screen Implicates miRNA-372 and miRNA-373 As Oncogenes in Testicular Germ Cell Tumors. Cell. 2006, 124 (6): 1169-1181. 10.1016/j.cell.2006.02.037
Tuveson DA, Shaw AT, Willis NA, Silver DP, Jackson EL, Chang S, Mercer KL, Grochow R, Hock H, Crowley D, Hingorani SR, Zaks T, King C, Jacobetz MA, Wang L, Bronson RT, Orkin SH, DePinho RA, Jacks T: Endogenous oncogenic K-ras(G12D) stimulates proliferation and widespread neoplastic and developmental defects. Cancer Cell. 2004, 5 (4): 375-387. 10.1016/S1535-6108(04)00085-6
Davies H, Bignell GR, Cox C, Stephens P, Edkins S, Clegg S, Teague J, Woffendin H, Garnett MJ, Bottomley W, Davis N, Dicks E, Ewing R, Floyd Y, Gray K, Hall S, Hawes R, Hughes J, Kosmidou V, Menzies A, Mould C, Parker A, Stevens C, Watt S, Hooper S, Wilson R, Jayatilake H, Gusterson BA, Cooper C, Shipley J, Hargrave D, Pritchard-Jones K, Maitland N, Chenevix-Trench G, Riggins GJ, Bigner DD, Palmieri G, Cossu A, Flanagan A, Nicholson A, Ho JW, Leung SY, Yuen ST, Weber BL, Seigler HF, Darrow TL, Paterson H, Marais R, Marshall CJ, Wooster R, Stratton MR, Futreal PA: Mutations of the BRAF gene in human cancer. Nature. 2002, 417 (6892): 949-954. 10.1038/nature00766
Yuen ST, Davies H, Chan TL, Ho JW, Bignell GR, Cox C, Stephens P, Edkins S, Tsui WW, Chan AS, Futreal PA, Stratton MR, Wooster R, Leung SY: Similarity of the phenotypic patterns associated with BRAF and KRAS mutations in colorectal neoplasia. Cancer Res. 2002, 62 (22): 6451-6455.
Singer G, Oldt R, Cohen Y, Wang BG, Sidransky D, Kurman RJ, Shih Ie M: Mutations in BRAF and KRAS characterize the development of low-grade ovarian serous carcinoma. J Natl Cancer Inst. 2003, 95 (6): 484-486.
Kim IJ, Park JH, Kang HC, Shin Y, Park HW, Park HR, Ku JL, Lim SB, Park JG: Mutational analysis of BRAF and K-ras in gastric cancers: absence of BRAF mutations in gastric cancers. Hum Genet. 2003, 114 (1): 118-120. 10.1007/s00439-003-1027-0
Kolch W: Meaningful relationships: the regulation of the Ras/Raf/MEK/ERK pathway by protein interactions. Biochem J. 2000, 351 Pt 2: 289-305. 10.1042/0264-6021:3510289
Seitz H, Royo H, Bortolin ML, Lin SP, Ferguson-Smith AC, Cavaille J: A large imprinted microRNA gene cluster at the mouse Dlk1-Gtl2 domain. Genome Res. 2004, 14 (9): 1741-1748. 10.1101/gr.2743304
Tsafrir D, Bacolod M, Selvanayagam Z, Tsafrir I, Shia J, Zeng Z, Liu H, Krier C, Stengel RF, Barany F, Gerald WL, Paty PB, Domany E, Notterman DA: Relationship of gene expression and chromosomal abnormalities in colorectal cancer. Cancer Res. 2006, 66 (4): 2129-2137. 10.1158/0008-5472.CAN-05-2569
Michael MZ, SM OC, van Holst Pellekaan NG, Young GP, James RJ: Reduced accumulation of specific microRNAs in colorectal neoplasia. Mol Cancer Res. 2003, 1 (12): 882-891.
The author(s) declare that they have no competing interests.
All authors participated in the design of experiments. EB was responsible for qRT-PCR studies and drafted the manuscript. EC and XA assisted with analysis, and contributed to drafting the manuscript. RM and RZ contributed to methods development and qRT-PCR analysis. NR, AA, and AN assisted with methods development and data analysis. IM and MM provided the CRC tumor samples and assisted with critical examination of the manuscript. JGF designed and coordinated of the study, and drafted the manuscript All authors read and approved the final manuscript.