| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Clinical Research |
1 Molecular Diagnostic Laboratory, Department of Clinical Biochemistry, Aarhus University Hospital, Aarhus N, Denmark; 2 Wilhelm Johannsen Centre for Functional Genome Research, Department of Cellular and Molecular Medicine, University of Copenhagen, Copenhagen N, Denmark; 3 Bioinformatics Research Center, University of Aarhus and 4 Department of Surgery P, Aarhus University Hospital, THG, Aarhus C, Denmark; and 5 Santaris Pharma, Hørsholm, Denmark
Requests for reprints: Claus L. Andersen, Molecular Diagnostic Laboratory, Department of Clinical Biochemistry, Aarhus University Hospital, Skejby, DK-8200 Aarhus N, Denmark. Phone: 45-89495186; Fax: 45-89496018; E-mail: cla{at}ki.au.dk.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
13% of all cancers and is the second most common cause of cancer death in the Western world (1–3). CRC is staged according to the extent of primary organ involvement and metastatic spread to lymph nodes or distant organs (4–7). Stage II (no lymph node or distant metastases) cancer patients have a 5-year survival of 75% to 80% when surgically treated. Despite surgical resection being highly effective for localized disease, a significant proportion of stage II patients (20–25%) develops recurrence and dies from the disease. At present, it is not possible to accurately differentiate between good and poor prognosis stage II patients. Consequently, there is an acute need for markers capable of distinguishing these two groups (8) as the poor prognosis group should receive chemotherapy. Molecularly, CRC can be discerned into two major subgroups: microsatellite stable (MSS) and microsatellite unstable cancers (MSI), respectively. Histologically, the two molecular subgroups are similar, but, clinically, they behave differently. The MSI cancers have a better prognosis: In general, they follow a more benign disease course; however, they respond less well to chemotherapy (9–11). Interestingly, recent studies have documented a link between the aberrant expression of a class of small noncoding RNAs (17–28 nucleotides), termed microRNAs (miRNA), and the pathogenesis of several cancer types including CRC (12–14). MiRNAs negatively regulate gene expression at the posttranscriptional level by cleavage and/or translational repression of their mRNA targets and their dysregulation is believed to promote the malignant behavior of tumors (12). Although exciting, the association between miRNAs and CRC warrants further investigation as most studies have been limited to a modest number of cancer tissue samples (15–20) and uncertainty exist as to the precise nature of the disturbed miRNA expression pattern observed in cancer tissue, which has been reported to be severely skewed toward an overall up- or down-regulation relative to normal tissue (14, 20). Differences between methodologies (14, 18–20) and clinicopathologic characteristics of the tissue used may explain this discrepancy. This latter issue, however, is difficult to account for because only a minority of the articles provides information such as tumor stage and location (16, 17). Accordingly, the association between miRNA expression and clinicopathologic variables, including metastatic relapse, remains to be investigated. In addition, no reports have described how the expression of miRNAs in MSS and MSI tumor subgroups relates to the expression of miRNAs in normal mucosa. Finally, in contrast to crude differences between cancerous and healthy colon tissue, virtually nothing is known about which specific cells contribute to the miRNA expression patterns observed as the location of only one miRNA, miR-143, has been described in colon tissue sections (21). All of these issues are important for providing a deeper understanding of the roles of miRNAs in CRC pathogenesis and as potential prognostic and diagnostic markers in CRC.
In this study, we have begun to address the issues raised above. We profiled the expression of 315 unique human miRNAs in 49 colon cancers and 10 normal mucosa samples using spotted locked nucleic acid (LNA)–based oligonucleotide microarrays and identified several miRNAs that may serve as potential diagnostic and prognostic disease markers. We confirmed the expression of selected miRNAs using real-time reverse transcription-PCR (RT-PCR) and in situ hybridization (ISH), and provide data suggesting that altered expression levels of miRNAs have a functional consequence in colon cancer.
| Materials and Methods |
|---|
|
|
|---|
MiRNA microarray expression profiling and data analysis. Microarrays were produced using a LNA-based oligonucleotide probe library (miRCURY LNA array ready to spot v.7.1; Exiqon A/S). Oligonucleotides were spotted in duplicates in a 10 µmol/L phosphate buffer on CodeLink microarray glass slides (GE Health Care) using a VersArray Chipwriter Pro system (Bio-Rad) as previously described (22). Two micrograms of sample RNA were directly labeled with Hy3 using the miRCURY LNA array labeling kit (Exiqon). As reference, we used a pool of RNA extracted from bladder, prostate, and colon tumors. For each experiment, we labeled 2 µg of the reference RNA with Hy5 using the LNA array labeling kit (Exiqon). Hybridization and washing of the microarray slides were performed as recommended by Exiqon. Scanning was performed using a ScanArray 4000 scanner (GSI Lumonics). After scanning of the microarrays, we used TIGR Spotfinder 2.23 software to generate raw intensity data, which were LOWESS (global) normalized using TIGR MIDAS 2.19 software (23). Average log 2 ratios were calculated from the normalized data based on the two measurements of each miRNA.
Hierarchical cluster analysis. For hierarchical cluster analysis, human miRNA probes on the array were median centered and normalized. Samples were then clustered hierarchically using Cluster 2.0. To evaluate the robustness of the formed clusters, the clustering analysis was performed using three different similarity metrics (correlation centered, correlation uncentered, and Spearman's rank) and two different clustering algorithms (average linkage clustering and complete linkage clustering). Treeview 2.0 was used for visualization. The cluster dendrograms shown in this article were produced using the correlation-centered metric and average linkage clustering algorithm.
Statistical analysis. For statistical analysis of microarray data, we used the significance analysis of microarrays (SAM) implemented in TIGR MEV 3.1 software (23). For all analyses, we performed 1,000 permutations of the data to select differentially expressed miRNAs. The significance cutoff,
value, was adjusted so as to set the median number of falsely discovered miRNAs to <0.00001, resulting in a median false discovery rate of <0.0001%.
For recurrence-free survival analysis, patients were dichotomized according to the expression level of a given miRNA. High and low expressions were defined according to the median expression of the miRNA in the tumor samples. Censored recurrence-free survival curves were plotted according to Kaplan-Meier. Univariate analysis was performed using the log-rank test. Multivariate analysis was performed using the Cox proportional hazards regression model. The software used for survival analysis was STATA 9.2 (StataCorp).
Construction of predictive biomarkers based on miRNA expression. Molecular predictors were built using the Prophet tool (24) with a forward filtering approach for gene selection based on the F ratio (between-group to within-group sums of squares). Three prediction algorithms were chosen for the analyses: prediction analysis of microarrays (PAM), diagonal linear discriminant analysis (DLDA), and k-nearest neighbor (k-NN). For the k-NN algorithm, the value of k (number of neighbors) was varied to achieve more robustness against noise; specifically, k was varied across the following values: 1, 2, 3, 4, and 5. The average of sensitivity, TP / (TP + FN), and specificity, TN / (TN + FP), was used to evaluate the performance of the predictors in leave-one-out cross-validation. TP, TN, FN, and FP denote true positives, true negatives, false negatives, and false positives, respectively. Accuracy was defined as (TP + TN) / (TP + FP + FN + TN).
Real-time RT-PCR. cDNA was synthesized from total RNA using gene-specific primers according to the TaqMan MicroRNA Assay protocol (Applied Biosystems) in a 15 µL reaction volume with 10 ng of RNA template. Reverse transcription was performed using the following program: 30 min at 16°C, 30 min at 42°C, 5 min at 85°C and then held at 4°C. Reverse transcription products were diluted 20-fold, and 2 µL were used in a total reaction volume of 20 µL for relative quantification by real-time PCR using an Applied Biosystems 7500 Sequence Detection system. Thermal cycling program used for quantification was as follows: 50°C for 2 min and 95°C for 10 min, followed by 50 cycles of 95°C for 15 s and 60° for 1 min. For normalization, random nonamer-primed cDNA synthesis was made in parallel, and the relative expression of the Ubiquitin gene (UBC) was quantified as previously described (25). Each measurement was performed in triplicate and no-template controls were included for each assay. Relative expression values were obtained using 5-point 10-fold dilution curves. The dilution curve was created using a cDNA pool of each of the test cDNAs.
In situ detection of miRNAs. The in situ detection was performed on 5-µm paraffin sections from matched normal mucosa and adenocarcinoma material. Sections were deparaffinized and subsequently treated with 1 mol/L HCl for 15 min at room temperature to inhibit endogenous alkaline phosphatases. The sections were then digested with proteinase K [50 µg/mL in 50 mmol/L Tris-HCl (pH 7.5)] for 30 min at room temperature. After proteinase K digestion, the sections were fixed in 4% paraformaldehyde at 4°C for 10 min and prehybridized in prehybridization buffer [50% formamide, 50 µg/mL heparin, 5 mmol/L EDTA, 5x SSC, 50 µg/mL yeast tRNA, 0.1% Tween 20, and 2% blocking reagent (Roche)] for 3 h and 15 min at 57°C. Subsequently, the slides were hybridized with 10 pmol probe (LNA-modified and DIG-labeled oligonucleotide; Exiqon) complementary to miR-145 and miR-320 in hybridization buffer [50% formamide, 0.3 mol/L NaCl, 20 mmol/L Tris-HCl (pH 8.0), 5 mmol/L EDTA, 10 mmol/L NaPO4 (pH 8.0), 10% dextran sulfate, 0.5 mg/mL yeast tRNA, and 1x Denhardt's solution] overnight at 57°C. After incubation with anti–DIG-AP Fab fragments conjugated to alkaline phosphatase diluted 1:250 in blocking solution (Roche), the hybridized probes were detected by applying nitroblue tetrazolium/5-bromo-4-chloro-3-indolyl phosphate color substrate (Roche) to the slides. Slides were mounted in AquaTex (Merck) and analyzed using a Zeiss Axiovert 40 CFL microscope. No-probe controls were included for each hybridization procedure, and specificity was verified using probes that were mismatched at two (miR-145) or three (miR-320) nucleotide positions.
Cell culture and transfection. LS174T and DLD1 cell lines were maintained in RPMI 1640 (+ 25 mmol/L HEPES; Invitrogen/Life Technologies), whereas McCoy's 5A medium was used for HCT116 cells. The medium was supplemented with 10% fetal bovine serum (FBS; Invitrogen/Life Technologies) and 1% penicillin-streptomycin (Invitrogen/Life Technologies) for all cell lines but DLD1 cells, which grew in 5% FBS.
Antisense LNA oligonucleotides (Exiqon) and precursor miRNAs (pre-miRs; Ambion/Applied Biosystems) were transfected using Lipofectamine 2000 reagent (Invitrogen), whereas cells were in suspension (reverse transfection) according to the manufacturer's instructions. Transfected cells were seeded in 96-well plates, and the medium was replaced 12 h posttransfection. Final concentrations of antisense LNA oligonucleotides were 10, 20, or 50 nmol/L, whereas final concentrations of 2, 5, and 20 nmol/L were used for pre-miRs. The following LNA oligonucleotides and pre-miRs were used for the experiments: LNA-miR-92, LNA-miR-20a, LNA-scrambled, pre-miR-145, and scrambled pre-miR.
Cell viability and cell death assays. Cell viability was assessed using the Cell Proliferation kit I [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT)] from Sigma-Aldrich according to the manufacturer's instructions. Cell death was quantified using the Cytotoxicity Detection kit [lactate dehydrogenase (LDH)] from Roche essentially as described by the manufacturer; however, cells were lysed in a buffer consisting of 0.5% Triton X-100, 25 mmol/L HEPES (pH 7.0), 5 mmol/L MgCl2, 1 mol/L EGTA supplemented with serine protease inhibitor, Pefabloc (Roche) 0.5 mmol/L. Transfected cells were seeded in quadruplicate wells for each condition in 96-well plates and assayed 72 h posttransfection. Spectrophotometry for the MTT and LDH assays was performed on a microplate reader (Labsystems Multiscan MCC/340) at 540 and 492 nm, respectively. Absorbance at 690 nm was used as reference.
| Results |
|---|
|
|
|---|
To investigate whether miRNA profiles were informative with regard to tissue type and tumor subgroup, we performed an unsupervised hierarchical cluster analysis. Samples were essentially partitioned into two groups—the first containing the normal mucosa biopsies and the vast majority of the MSI tumors and the other containing most of the MSS tumors (Fig. 1A ). When clustering tumor samples separately, the MSS and MSI tumor subtypes segregated into two clusters with statistical significance (Fig. 1B; Fisher's exact test, P < 0.005). Collectively, these observations indicated that the clinically relevant MSS and MSI tumor subtypes differed with regard to miRNA expression patterns, and that MSI tumors seemed to resemble normal mucosa samples more so than MSS tumors. Yet, MSI tumors were readily distinguishable from normal tissue as indicated by a perfect split between the two groups when restricting the cluster analysis to normal mucosa and MSI samples only (Fig. 1C).
|
Differentially expressed miRNAs between normal mucosa and colon cancer subtypes, MSS and MSI. We next used SAM to identify which specific miRNAs were differentially expressed between normal mucosa and each tumor subtype, MSI and MSS. When comparing normal mucosa to MSI colon cancers, we found 25 miRNAs to be differentially expressed (7 down; 18 up), and 54 miRNAs were differentially expressed when comparison was made to MSS colon cancers (29 down; 25 up). The most significantly dysregulated miRNAs are listed in Table 1 (for a complete list, see Supplementary Table S2). In total, 60 miRNAs were dysregulated in cancer with 19 miRNAs being dysregulated in both MSS and MSI subtypes. Notably, the majority of dysregulated miRNAs, 63% (38 of 60), had not previously been reported to be dysregulated in colon cancer, thus providing a novel association of several miRNAs with this disease. Of the miRNAs previously reported to be dysregulated in colon cancer, 59% (13 of 22) exhibited concordant patterns of dysregulation (Supplementary Table S3). Therefore, in addition to identifying numerous novel miRNAs associated with colon cancer, we successfully validated several miRNAs previously reported to be dysregulated in colon cancer. When comparing MSS and MSI subtypes directly using SAM, 11 miRNAs were differentially expressed with the average expression being highest in MSI tumors (Supplementary Table S2; Table 1). In agreement with MSI cancers exhibiting a miRNA expression pattern closer to normal mucosa than MSS cancers, most of these miRNAs 91% (10 of 11) were not significantly differentially expressed between MSI and normal mucosa, whereas the majority 64% (7 of 11) of the same miRNAs were significantly down-regulated in MSS cancers relative to normal mucosa.
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
Intriguingly, tumor subtypes distinguished by microsatellite status, MSI or MSS, exhibited distinct miRNA expression patterns, and, accordingly, the expression of specific miRNAs could be used to correctly classify the majority of colon cancers as either MSI or MSS. These observations are consistent with a recent study by Lanza and colleagues (31) who also found that miRNA profiles could accurately predict microsatellite status in a set of 39 colon cancers. Together, these findings support the notion that the type of genetic instability in cancer tissue is reflected at the miRNA level, which is also consistent with results from our unsupervised cluster analyses showing that miRNA expression profiles from MSI tumors, known to generally be near diploid, have a tendency to cluster more closely with normal tissue than do MSS tumors, known to exhibit aneuploidy with varying chromosome numbers (32). Thus, chromosomal abnormalities may directly or indirectly contribute to the altered expression of several miRNAs in colon cancer. It is noteworthy, however, that we did not identify the same specific miRNAs that distinguished MSI from MSS cancers as Lanza and colleagues (Supplementary Table S3; ref. 31). Proper comparison of the results is obscured by differences in analytic platforms. Specifically, the probes on the array platform used by Lanza and colleagues (31) were designed against precursor miRNAs, and this platform has previously been reported to yield data that do not correlate well with real-time RT-PCR data measuring the expression of mature miRNAs (33). By contrast, our array platform with probes designed against mature miRNAs produces data that for the majority of investigated miRNAs correlate well with real-time RT-PCR data.
We found several miRNAs to be associated with metastatic recurrence of stage II MSS colon cancers. By predicting recurrence status of the primary tumor samples based on miRNA expression levels, we found 17 miRNAs to give the best performance of 81% accuracy, 77% sensitivity, and 83% specificity. These results are highly similar to the predictive power obtained by other investigators when classification was based on mRNA microarray expression profiles (34). It may prove beneficial to combine miRNA/mRNA gene expression data, thereby possibly enhancing the accuracy of the molecular predictors, as shown by Lanza and colleagues (31), when predicting microsatellite status. Our analysis of the probability of recurrence-free survival as a function of the expression levels of either of the two top-ranked predictive miRNAs from our predictive recurrence signature, miR-320 and miR-498, underpins the prognostic potential of miRNAs in relation to stage II colon cancer. Both miRNAs independently of known risk factors such as age, pT stadium, differentiation grade, and histologic subtype significantly predict metastatic recurrence. Ultimately, it may be possible to use miRNA expression profiles, either alone or in combination with other biomarkers including mRNA expression data, to identify the stage II colon cancer patients who are at high risk of metastatic recurrence and thus likely to benefit from adjuvant chemotherapy. Although our study is the first to investigate whether miRNA expression levels can predict recurrence of stage II colon cancer, the prognostic value of miRNAs have also been suggested from other studies investigating disease course in lung cancer, neuroblastoma, and leukemia (26–28). It is interesting to note that, in general, a few hundred miRNAs seem to be at least as informative with respect to diagnostic and prognostic characteristics as the abundant class of mRNAs containing thousands of variants when analyzing microarray expression profiles. However, as has been suggested, an effort to achieve universal standards for miRNA arrays would be required to allow optimal comparison of miRNA array data from multiple laboratories (35). Furthermore, an independent validation study is needed before firm conclusions can be made on clinical utilization.
Although the literature describing miRNA expression and function is exponentially growing, little is still known about the in vivo localization of miRNAs. This issue is important, not only for assessing whether differences in miRNA expression levels between crude tissue biopsies arise from variations in cellular composition, or true differences between normal and transformed cells, but also for understanding the biological roles of miRNAs linked to particular cell types. By performing ISH with probes against miR-320 and miR-145, we found that both miRNAs are expressed in the epithelial cells of normal colon tissue. Interestingly, miR-320 expression was highest in the proliferative compartment of the colonic crypts, and in the small panel of tumors investigated, miR-320 was highly expressed in cancer cells. These observations are consistent with our microarray data showing that miR-320 is more abundantly expressed in cancer relative to normal tissue. It may therefore be speculated that miR-320 expression is involved in maintaining the undifferentiated and proliferative tumor phenotype. Conversely, microarray expression data for miR-145 showed that expression levels varied across cancer samples, but, overall, miR-145 expression was drastically reduced relative to normal tissue. ISH analysis showed that in four of five tumors, miR-145 expression seemed to be reduced compared with healthy colon tissue.
Strikingly, by restoring miR-145 levels in a panel of three different colon carcinoma cell lines, we found cell growth to be markedly retarded and the reduction in cell viability was subsequently found to be caused by an increased rate of cell death. In contrast, in vitro inhibition of two miRNAs highly expressed in tumor tissue, miR-20a and miR-92, did not affect cell viability. This was unexpected, considering that these miRNAs are part of the miR-17-92 cluster reported to be oncogenic in other cancer types (36, 37). Although our data indicated that the viability of colon cancer cell lines were unaffected by inhibition of miR-20a and miR-92, other cancer hall marks such as motility, invasion, and angiogenesis might have been influenced. This, however, was beyond the scope of this study to investigate. Importantly, the growth inhibition we observed by miR-145 are in agreement with a recent study reporting miR-143 and miR-145 to negatively influence cell growth in two colon carcinoma cell lines, of which one overlapped with the ones we tested (17). This miRNA has also been reported to be down-regulated in colon carcinomas by other investigators (15, 18), and recent work by Shi and colleagues (38) showed that one target of miR-145 is the insulin receptor substrate-1 conveying mitogenic and antiapoptotic signaling. Thus, increasing the levels of miR-145 in transformed cells may prove to be an efficient way of inhibiting tumor growth, and our results support the suggestion that miR-145 is an interesting candidate for the development of future RNA medicine against cancer as has been put forward for hematopoietic malignancies (39) and colon cancer (38).
In conclusion, we showed significant correlations between the expression of miRNAs and histologic as well as clinical variables in colon cancer. Our findings provide the basis for an in-depth study of the role of certain miRNAs as clinical prospective markers and as new pharmaceutical entities.
| Disclosure of Potential Conflicts of Interest |
|---|
|
|
|---|
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank Inge Lis Thorsen, Gitte Høj, and Lisbet Kjeldsen for their excellent technical assistance, and the staff at the Department of Pathology at Aarhus University Hospital.
| Footnotes |
|---|
Received 11/ 4/07. Revised 3/25/08. Accepted 5/22/08.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T Paranjape, F J Slack, and J B Weidhaas MicroRNAs: tools for cancer diagnostics Gut, November 1, 2009; 58(11): 1546 - 1554. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. VALERI, C. M. CROCE, and M. FABBRI Pathogenetic and Clinical Relevance of MicroRNAs in Colorectal Cancer Cancer Genomics Proteomics, July 1, 2009; 6(4): 195 - 204. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Dyrskjot, M. S. Ostenfeld, J. B. Bramsen, A. N. Silahtaroglu, P. Lamy, R. Ramanathan, N. Fristrup, J. L. Jensen, C. L. Andersen, K. Zieger, et al. Genomic Profiling of MicroRNAs in Bladder Cancer: miR-129 Is Associated with Poor Outcome and Promotes Cell Death In vitro Cancer Res., June 1, 2009; 69(11): 4851 - 4860. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Friedman, G. Liang, C.-C. Liu, E. M. Wolff, Y. C. Tsai, W. Ye, X. Zhou, and P. A. Jones The Putative Tumor Suppressor microRNA-101 Modulates the Cancer Epigenome by Repressing the Polycomb Group Protein EZH2 Cancer Res., March 15, 2009; 69(6): 2623 - 2629. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Sachdeva, S. Zhu, F. Wu, H. Wu, V. Walia, S. Kumar, R. Elble, K. Watabe, and Y.-Y. Mo p53 represses c-Myc through induction of the tumor suppressor miR-145 PNAS, March 3, 2009; 106(9): 3207 - 3212. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. E. AHMED, P. W. VOS, C. JEFFRIES, J. E. WILEY, D. A. WEIDNER, H. MOTA, C. BONNERUP, C. SIBATA, and R. R. ALLISON Differences in mRNA and microRNA Microarray Expression Profiles in Human Colon Adenocarcinoma HT-29 Cells Treated with either Intensity-modulated Radiation Therapy (IMRT), or Conventional Radiation Therapy (RT) Cancer Genomics Proteomics, March 1, 2009; 6(2): 109 - 127. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |