Schetter AJ, Leung SY, Sohn JJ, Zanetti KA, Bowman ED, Yanaihara N, Yuen ST, Chan TL, Kwong DLW, Au GKH, Liu C, Calin GA, Croce CM, Harris CC. MicroRNA Expression Profiles Associated With Prognosis and Therapeutic Outcome in Colon Adenocarcinoma. JAMA. 2008;299(4):425–436. doi:10.1001/jama.299.4.425
Author Affiliations: Laboratory of Human Carcinogenesis, Center for Cancer Research (Drs Schetter, Sohn,
Zanetti, Yanaihara, and Harris and Ms Bowman) and Cancer Prevention Fellowship Program, Office of Preventive Oncology (Drs Schetter and Zanetti), National Cancer Institute, National Institutes of Health,
Bethesda, Maryland; Departments of Pathology (Drs Leung, Yuen, and Chan) and Clinical Oncology (Drs Kwong and Au), University of Hong Kong, Queen Mary Hospital, Pokfulam, Hong Kong; Department of Molecular Virology, Immunology, and Medical Genetics, Ohio State University Comprehensive Cancer Center, Columbus, Ohio (Drs Liu and Croce); and Department of Experimental Therapeutics, M. D. Anderson Cancer Center,
Houston, Texas (Dr Calin).
Context MicroRNAs have potential as diagnostic biomarkers and therapeutic targets in cancer. No study has evaluated the association between microRNA expression patterns and colon cancer prognosis or therapeutic outcome.
Objective To identify microRNA expression patterns associated with colon adenocarcinomas, prognosis, or therapeutic outcome.
Design, Setting, and Patients
MicroRNA microarray expression profiling of tumors and paired nontumorous tissues was performed on a US test cohort of 84 patients with incident colon adenocarcinoma, recruited between 1993 and 2002.
We evaluated associations with tumor status, TNM staging, survival prognosis, and response to adjuvant chemotherapy. Associations were validated in a second, independent Chinese cohort of 113 patients recruited between 1991 and 2000, using quantitative reverse transcription polymerase chain reaction assays. The final date of follow-up was December 31, 2004, for the Maryland cohort and August 16, 2004, for the Hong Kong cohort.
Main Outcome Measures MicroRNAs that were differentially expressed in tumors and microRNA expression patterns associated with survival using cancer-specific death as the end point.
Thirty-seven microRNAs were differentially expressed in tumors from the test cohort. Selected for validation were miR-20a, miR-21, miR-106a, miR-181b, and miR-203,
and all 5 were enriched in tumors from the validation cohort (P < .001). Higher miR-21 expression was present in adenomas (P = .006)
and in tumors with more advanced TNM staging (P < .001).
In situ hybridization demonstrated miR-21 to be expressed at high levels in colonic carcinoma cells. The 5-year cancer-specific survival rate was 57.5% for the Maryland cohort and was 49.5% for the Hong Kong cohort. High miR-21 expression was associated with poor survival in both the training (hazard ratio,
2.5; 95% confidence interval, 1.2-5.2) and validation cohorts (hazard ratio, 2.4; 95% confidence interval, 1.4-3.9), independent of clinical covariates, including TNM staging, and was associated with a poor therapeutic outcome.
Expression patterns of microRNAs are systematically altered in colon adenocarcinomas. High miR-21 expression is associated with poor survival and poor therapeutic outcome.
Colon adenocarcinoma is a major cause of cancer mortality worldwide.1
Colorectal cancer is the third most common and second leading cause of cancer death in the United States.2
Even though 5-year mortality rates have modestly declined over the last 3 decades,3
there is still a need to identify new prognostic biomarkers and therapeutic targets for this disease. Currently, chemotherapy has significant therapeutic value, but surgery is the only curative form of treatment.4
Ideal therapeutic targets should be causally associated with disease and amenable to designing therapeutic interventions, whereas ideal biomarkers should be easy to measure and have strong associations with clinical outcomes. MicroRNAs could match both criteria.5- 8
MicroRNAs are 18- to 25-nucleotide, noncoding RNA molecules that regulate the translation of many genes.9
Since their discovery,10,11
microRNAs have been found to regulate a variety of cellular processes including apoptosis,12- 14
and cell proliferation.16
MicroRNAs may also have a causal role in carcinogenesis.5,17,18
MicroRNA expression levels are altered in most tumor types,19,20
including colon tumors.20- 23
Experimental manipulation of specific microRNAs modulates tumor development in mouse-model systems.16,24- 26
The prognostic potential of microRNAs has also been demonstrated for chronic lymphocytic leukemia,6
If aberrant microRNA expression is causal to carcinogenesis,
inhibiting specific microRNAs may have therapeutic implications. Modified antisense oligonucleotides can easily be designed to specifically inhibit microRNA function.29
Antagomirs are one type of antisense oligonucleotide that has proven effective at inhibiting microRNA function in vivo in mice.30 The ease of designing specific inhibitors of microRNA function makes them candidates for therapeutic targets.
Given the therapeutic and prognostic potential for microRNAs in cancer, we evaluated microRNA profiles of colon tumors and paired nontumorous tissue to study their potential role in tumor formation,
diagnosis, and response to chemotherapy in colon carcinoma.
Pairs of primary colon tumor and adjacent nontumorous tissues came from 84 patients recruited from the University of Maryland Medical Center or Baltimore Veterans Affairs Medical Center between 1993 and 2002, and from 113 patients recruited from Queen Mary Hospital in Hong Kong between 1991 and 2000. Cases with familial adenomatous polyposis or human nonpolyposis colorectal cancer were excluded from this study.
Tissues were flash frozen after surgery. Detailed backgrounds for each tissue donor, including age, sex, clinical staging, tumor location, survival times from diagnosis, and receipt of adjuvant chemotherapy have been collected. The final date of follow-up was December 31, 2004 for the Maryland cohort and August 16, 2004, for the Hong Kong cohort. Tumor histopathology was classified according to the World Health Organization Classification of Tumor system.1 Adenoma tissue was obtained from the Cooperative Human Tissue Network. This study was approved by the Institutional Review Board of the National Institutes of Health, the Institutional Review Board of the University of Hong Kong/Hospital Authority Hong Kong West Cluster, and the Institutional Review Board for Human Subject Research at the University of Maryland. Race was self-reported as white or black.
RNA from frozen tissue samples was extracted using standard TRIZOL (Invitrogen, Carlsbad, California) methods. MicroRNA microarray profiling was performed as previously described.31 Briefly, 5 μg of total RNA was labeled and hybridized to each microRNA microarray (Ohio State microRNA microarray version 2.0, Columbus) containing quadruplicates of 389 human microRNA probes. Tumor/nontumorous pairs of tissues were profiled at the same time. Slides were scanned using a PerkinElmer ScanArray LX5K scanner (Perkin Elmer, Waltham, Massachusetts).
The data discussed in this publication have been deposited in National Center for Biotechnology Information's (NCBI’s) Gene Expression Omnibus (NCBIGEO GSE 7828).
The data were preprocessed by the statistical software R 2.5.0 (R Foundation for Statistical Computing, Vienna, Austria) to remove probes with higher background intensities than foreground and probes with inconsistent measurements across the quadruplicates. The data were normalized by locally weighted scatter plot smoothing LOESS and imported into Biometric Research Branch (BRB) array tools 3.5.0 (http://linus.nci.nih.gov/BRB-ArrayTools.html) for subsequent microarray. Probes with values missing from more than 20% of the arrays were removed from the analysis leaving 230 probes. This filtering method was decided a priori to eliminate probes whose microRNAs expression levels were thought to be unreliable.
Class comparison analysis using paired t tests identified microRNAs that were differentially expressed in tumors (P < .001). Class prediction algorithms in BRB array tools were used to determine whether microRNA microarray expression patterns could accurately differentiate between tumor and paired nontumor tissue. For these analyses, 3 nearest neighbors and nearest centroid algorithms were arbitrarily chosen and percent accuracy reports the percentage of tissues that were correctly identified.
These algorithms were also used for quantitative reverse transcription polymerase chain reaction (RT-PCR) data in the Hong Kong validation cohort.
To initially search for microRNAs associated with poor survival,
tumor:nontumor (T:N) microRNA expression ratios were analyzed in the Maryland cohort using microarray data. Tumor:nontumor expression ratios for microRNAs were created by subtracting the log2 nontumor from the log2 tumor expression values. MicroRNAs missing more than 25% of T:N ratios were filtered out leaving 208. Expression data were dichotomized into clearly defined high and low groups to examine associations with microRNA expression and survival. Tumor:nontumor expression ratios were dichotomized with the highest tertile classified as high and the lower 2 tertiles classified as low. This cutoff was set based on associations within the test cohort prior to analyzing the validation cohort. Once set, this high-low cutoff was used universally throughout this study. To analyze associations with tumor expression and nontumor expression with survival using microarray data, the array data had to be normalized based on the day of microarray profiling to remove systematic bias introduced from the day-to-day variability observed in the microarray data acquisition. To do this, for each given day, the highest one-third expressing values were labeled high and the lowest two-thirds were labeled low, consistent with the predetermined cutoff that was used for this study.
Quantitative RT-PCR of microRNAs was performed using Taqman MicroRNA assays (Applied Biosystems, Foster City, California) according to the manufacturer's instructions with the 7500 real-time RT-PCR system (Applied Biosystems, Foster City) using expression levels of the small nuclear RNA, U6B, as the normalization control. All assays were performed in duplicate (miR-20a, miR-203) or triplicate (miR-21, miR-106a, miR-181b). Quantitative RT-PCR for miR-21, miR-106a,
and miR-181b was performed by one of the investigators (A.J.S.) who was blinded to the survival outcomes and clinical data for members of the validation cohort.
In situ hybridization was performed with probes for human miR-21, scramble, and U6 (Exiqon, Woburn, Massachusetts)
with a modified version of the manufacturer's protocol for formalin-fixed paraffin-embedded tissue written by W. Kloosterman (http://www.exiqon.com/uploads/LNA_52-_FFPE_miRNA_in_situ_protocol.pdf) on human colon tissue. Modifications included the use of polyclonal rabbit anti-DIG/HRP-conjugated antibody and DakoCytomation GenPoint Tyramide Signal Amplification System (DakoCytomation, Carpinteria,
California), and VECTOR NovaRed substrate (Vector Laboratories, Burlingame,
California). Images were taken on an Olympus BX40 microscope using the Olympus DP70 digital camera and DP controller software (Olympus,
Expression graphs and Wilcoxon matched-pairs tests were used to analyze differences in microRNA expression between tumors and paired nontumorous tissue as well as differences between adenoma and paired nonadenoma tissue for all quantitative RT-PCR data using Graphpad Prism 4.0 (Graphpad Software Inc, San Diego, California). All trend tests reported are Cuzick nonparametric test for trend across ordered groups32
and were performed using Stata 9.2 (StataCorp LP, College Station, Texas). Associations with prognosis in the validation cohort were considered statistically significant only if the P value were less than .01 to adjust for multiple comparisons testing (5 tests using a Bonferroni correction).
KaplanMeier analysis was performed with WINSTAT 2001 (R Fitch Software, Bad Krozingen, Germany). Multivariate Cox regression analysis was performed using StataCorp 9.2. For these models, we dichotomized age as 50 years or older vs younger than 50 years because the recommended screening age for colon cancer is at age 50 years; TNM staging was dichotomized based on metastasic vs nonmetastasic disease. One patient in the Maryland cohort died on the day of surgery and was included in Kaplan-Meier analysis but removed for Cox regression analysis.
Analyses involving response to adjuvant therapy included only TNM stage II and III cases because treatment in stage IV is palliative care and TNM stage I cases have excellent survival prognosis regardless of therapy. Univariate Cox regression was performed on each clinical covariate to examine influence of each on patient survival. Final multivariate models were based on stepwise addition and removal of clinical covariates found to be associated with poor survival in univariate models (P < .10). A Wald statistic of P < .05 was used as the criterion for inclusion in final multivariate models. All stepwise addition models gave the same final models as stepwise removal models. All P values reported are 2-sided. All univariate and multivariate Cox regression models were tested for proportional hazards assumptions based on Schoenfeld residuals, and no model violated these assumptions.
The characteristics of the patients with incident colon adenocarcinoma in the test cohort (from Baltimore) and the validation cohort (from Hong Kong) are shown in Table 1.
The median follow-up time was 68.0 months for the Baltimore cohort and 84.6 months for the Hong Kong cohort. The 2 cohorts were similar in TNM staging, tumor histology, and cancer-specific mortality rates.
The 5-year survival rate was 57.5% for the US cohort and 49.5% for the Hong Kong cohort and were not significantly different from one another (P = .49, Kaplan-Meier test). In addition to the racial, geographic, and cultural differences between these 2 cohorts, the Baltimore cohort was considerably older (average 64.6 years vs 55.8 years) had a higher percentage of men (79% vs 50%).
We compared microRNA profiles of 84 pairs of colon tumor and adjacent nontumorous tissues in the Baltimore cohort using microRNA microarrays.31
Tumor microRNA profiles were distinctly different from nontumor profiles. Using class comparison analysis in BRB array tools, 37 independent microRNAs were found to be differentially expressed in tumors (P < .001
with a false-discovery rate <0.5%; Table 2). Twenty-six microRNAs were expressed at higher levels in tumors with miR-21 enriched the most at 1.8-fold. Global microRNA profiles distinguish between tumor and paired nontumorous tissue with 89% accuracy using either the 3 nearest neighbors or nearest centroid class prediction algorithms within BRB array tools (10-fold cross validation repeated 100 times), suggesting a systematic change in microRNA expression patterns during tumor formation.
We next performed a preliminary analysis to identify whether any of the 37 differentially expressed (P < .001)
microRNAs were associated with cancer survival. We analyzed individual microRNA T:N expression ratios for associations with poor prognosis in the Maryland test cohort. Tumor:nontumor microRNA expression ratios were classified as high based on highest tertile. We searched for any microRNA for which high T:N ratios were associated with cancer survival (P < .05) using Cox regression. Five microRNAs satisfied these criteria. High expression of miR-20a (hazard ratio [HR], 2.2; 95% confidence interval [CI], 1.1-4.6; P = .03), miR-21 (HR, 2.5; 95% CI, 1.2-5.0; P = .01), miR-106a (HR, 2.3; 95% CI, 1.1-4.5; P = .02), miR-181b (HR,
2.0; 95% CI, 1.0-3.9; P = .04),
and miR-203 (HR, 3.1; 95% CI, 1.5-6.4; P = .003) were each associated with poor survival and were selected for further analysis.
To validate overall differences in microRNA expression between tumor and nontumorous tissue, we measured the expression levels of these 5 microRNAs with quantitative RT-PCR in tumor and paired nontumorous tissue in the independent validation cohort of 113 patients with incident colon cancer recruited from Hong Kong, China (Table 1). MiR-20a (2.3-fold), miR-21 (2.8-fold), miR106a (2.4-fold), miR-181b (1.4-fold), and miR-203 (1.8-fold)
were all expressed at higher levels in tumors (P < .001,
Wilcoxon matched-pairs test; Table 3).
Most tumors (89% for miR-20a, 87% for miR-21, 90% for miR-106a,
71% for miR-181b, and 74% for miR-203) had higher expression of these microRNAs than paired nontumorous tissue. Expression patterns for these 5 microRNAs distinguish tumor vs paired nontumor status with 96% or 98% accuracy based on 3 nearest neighbors or nearest centroid algorithms, respectively (10-fold cross validation, repeated 100 times).
Colon adenocarcinomas from 89% to 93% of the patients in this study were of a typical histology. A minority of tumors were of mucinous adenocarcinoma (8 of 84 patients [10%] in the US cohort; 7 of 113
patients [6%] in the Hong Kong cohort), adenosquamous carcinoma (1
of 84 patients [1%] in the US cohort), or signet ring cell carcinoma (1 of 113 patients [1%] in the Hong Kong cohort) histologies (Table 1). Different subtypes of adenocarcinomas can be associated with different clinical outcomes, including survival prognosis.33 To remove potential confounding associated with histology, we excluded all patients (9 in the US cohort;
8 in the Hong Kong cohort) with mucinous adenocarcinomas, adenosquamous carcinomas, and signet ring cell carcinomas from the initial analysis.
We found high T:N expression ratios for miR-20a, miR-21, miR-106a, miR-181b, and miR-203 to be associated with poor survival in the Maryland test cohort. These associations could be due to microRNA expression levels in the tumor tissue, the surrounding nontumorous tissue, or a combination of both.
To distinguish these possibilities, we analyzed the association of microRNA expression in tumors and paired nontumorous tissues separately.
High expression levels in tumors (based on highest tertile) for miR-20a (HR, 2.7; 95% CI, 1.3-5.8; P = .01), miR-21 (HR, 2.5;
95% CI, 1.2-5.2; P = .01), miR-106a (HR, 2.4; 95% CI, 1.2-5.1; P = .02), miR-181b (HR,
3.2; 95% CI, 1.6-6.7; P = .002),
and miR-203 (HR, 3.3; 95% CI, 1.5-7.1; P = .001) were each associated with a poor survival in the Maryland test cohort. No significant association with microRNA expression in nontumorous tissue was observed for any of the 5 microRNAs.
The associations with microRNA expression and survival in the test cohort were made in the context of a microarray experiment in which we were evaluating the dichotomized expression of 37 microRNAs.
To validate these findings, we used quantitative RT-PCR to measure tumor- and nontumor-expression levels for these 5 microRNAs in the Hong Kong validation cohort and analyzed associations with prognosis.
We dichotomized high and low expression for each microRNA based on highest tertiles, consistent with our methods in the test cohort.
High miR-21 tumor expression was associated with poor prognosis in the Hong Kong validation cohort (P = .001, Kaplan-Meier log-rank test) while expression in nontumorous tissue was not (Figure 1), consistent with associations in the Maryland test cohort.
We did not find statistically significant associations with prognosis and expression of miR-20a, miR-106a, miR-181b, or miR-203 in this cohort.
Multivariate Cox proportional hazards analysis was used to further evaluate the association of miR-21 expression in tumors with prognosis in both the Maryland test cohort and the Hong Kong validation cohort (Table 4)
to evaluate the potential for miR-21 expression as a prognostic biomarker. The dichotomized miR-21 expression values were not associated with age, sex, race, or tumor location (Fisher exact test). In univariate analysis for the Maryland test cohort, high expression of miR-21 in tumors (HR, 2.5; 95% CI, 1.2-5.2; P = .01)
and TNM staging (HR, 3.5; 95% CI, 1.6-7.9; P = .002)
was associated with prognosis while age, sex, race, and tumor location were not. In the final multivariate model, which included miR-21 expression and TNM staging, high miR-21 expression in tumors was associated with a poor survival prognosis independent of tumor staging (HR, 2.7; 95%
CI, 1.3-5.5; P = .008; Table 4).
In the Hong Kong validation cohort, the dichotomized values for miR-21 expression in tumors were not significantly associated with age, sex, tumor histology, or tumor location (Fisher exact test). High miR-21 expression in tumors (HR, 2.4; 95% CI, 1.4-3.9; P = .002)
and TNM staging (HR, 4.7; 95% CI, 2.4-9.5; P < .001)
were significantly associated with survival in univariate models while age, sex, and tumor location were not (Table 4). In the final multivariate Cox regression model, including miR-21 expression and TNM staging, high miR-21 expression in tumors was associated with poor survival prognosis (HR, 2.4; 95% CI, 1.4-4.1; P = .002)
independent of other clinical covariates, consistent with findings in the Maryland test cohort.
Adenomas represent a precursor stage for colon adenocarcinomas.34
We tested miR-20a, miR-21, miR-106a, miR-181b, and miR-203 expression levels by quantitative RT-PCR in 18 pairs of adenoma and adjacent nonadenoma tissue that were acquired from the Cooperative Human Tissue Network. Using only 18 pairs of tissues reduced the power to detect differences in these analyses and may result in false-negatives. However, miR-21 was significantly enriched at 1.6-fold higher (P = .006, Wilcoxon matched-pairs test; Table 3). Adenoma tissue expressed higher levels of in 15 of 18 matched pairs.
If miR-21 expression is causal to the progression of colon cancer, expression of miR-21 may be associated with more advanced stages of the disease. Patients were stratified based on the diagnosis of adenoma and TNM staging, in which adenoma was considered the least advanced and TNM stage IV was most advanced. Adenomas expressed lower levels of miR-21 than tumors from the validation cohort (P < .001,
Mann-Whitney test). More advanced tumors expressed higher levels of miR-21 using either microarray data from the Maryland test cohort (P = .04, test for trend) or the quantitative RT-PCR data from the Hong Kong validation cohort (test for trend, P < .001; Figure 2).
Although we found that high expression of miR-21 in tumors was associated with a worse survival outcome, these experiments did not identify the cells within a tumor that expressed miR-21. To identify these cells, we used in situ hybridization to visualize miR-21 expression in tumor and adjacent nontumor tissue (Figure 3). The miR-21 is expressed at higher levels in colonic epithelial cells in human tumor tissue, consistent for a role for miR-21 overexpression within tumor cells during colon carcinogenesis.
We analyzed associations with miR-21 expression and therapeutic outcomes in stage II and III cancer patients treated with adjuvant chemotherapy. Information on the administration of adjuvant chemotherapy was available for 47 of 65 stage II or III patients in the Maryland test cohort and for all patients in the Hong Kong validation cohort. In both cohorts chemotherapy regimens were primarily fluorouracil based (in forms of either intravenous 5-fluorouracil or oral drugs including tegafur with uracil [UFT]) with or without generics. Only patients with typical adenocarcinoma histology were used for this analysis, leaving 20 of 42 stage II or III patients who had received chemotherapy in the Maryland cohort. For those who received chemotherapy,
high miR-21 expression in tumors predicted worse overall survival (HR, 4.3; 95% CI, 1.1-16.4; P = .03)
giving preliminary support that high miR-21 is associated with poor therapeutic outcome.
For the Hong Kong validation cohort, 77 individuals with stage II or III cancer with typical adenocarcinoma histology were used for this analysis. Among the 36 patients who received adjuvant chemotherapy,
high miR-21 expression in tumors was associated with a poor response to therapy (HR, 3.5; 95% CI, 1.1-11.6; P = .04), consistent with observations in the Maryland cohort. Additionally, among the 25 patients with stage III cancer who received adjuvant chemotherapy, high miR-21 expression was associated with poor survival (HR, 3.9; 95%
CI, 1.2-12.9; P = .03). Analyses using cancer relapse as an end point resulted in similar associations with high miR-21 expression in tumors predicting a more rapid disease recurrence in patients with TNM stage III cancer who received adjuvant chemotherapy (HR, 3.5; 95% CI, 1.0-11.5; P = .04).
Both cohorts showed similar associations of high miR-21 expression in tumors with therapeutic outcomes, although neither cohort was large enough to perform a fully stratified analysis.
Although these are 2 different populations, the similarity of the associations with miR-21 expression, prognosis,
and therapeutic outcome allowed for pooled analysis using both cohorts.
Kaplan-Meier analysis of the pooled cohorts demonstrated that high miR-21 expression was associated with a poor prognosis in either stage II (P = .02) or stage III (P = .004) patients (Figure 4) further indicating its potential as a prognostic biomarker. Receipt of adjuvant chemotherapy was beneficial for patients with either stage II or stage III cancer, although this was only significant for stage III patients. For individuals who received adjuvant therapy, high miR-21 expression was associated with a poor therapeutic outcome in patients with stage II or III cancer (P = .003, Kaplan-Meier log rank) or in patients with stage III cancer alone (P = .007, Kaplan-Meier log rank; Figure 4). Multivariate Cox regression demonstrated that high miR-21 expression predicted poor prognosis (HR, 3.0; 95% CI, 1.7-5.4; P < .001)
and treatment with adjuvant chemotherapy was associated with better survival (HR, 0.4; 95% CI, 0.2-0.8; P = .004)
independent of other clinical covariates (Table 5).
We performed the largest study to date analyzing microRNA profiles in colon cancer tissues and the first, to our knowledge, using 2 independent cohorts. Thirty-seven microRNAs were differentially expressed in tumor tissues by microRNA microarray analysis in the Maryland test cohort.
Expression patterns of all 5 tested microRNAs were validated in the Hong Kong cohort. The discriminatory power of 5 microRNAs to differentiate between tumor and nontumorous tissue suggests that predictable and systematic changes of microRNA expression patterns may occur during tumorigenesis and may be representative of sporadic colon adenocarcinomas.
miR-20a, miR-21, miR-106a, miR-181b, and miR-203 were all found to be expressed at higher levels in colon tumors, although it is uncertain whether these changes in microRNA expression patterns are merely associated with colon cancer or causal to the histological progression to cancer. Our data are consistent with published studies that provide evidence for changes in microRNA expression promoting tumor formation, especially for miR-20a and miR-21. miR-20a is part of the miR-17-92 polycistronic microRNA cluster.35
Overexpression of this cluster enhances cell proliferation in vitro36
and accelerates tumor formation in animal models.16
Enforced expression of the miR-17-92 cluster causes increased tumor size and tumor vascularization in mice by negatively regulating the anti-angiogenic thrombospondin 1 (Tsp1) protein.24
Experimental evidence also suggests that increased miR-21 expression promotes tumor development.
miR-21 is expressed at high levels in most solid tumors.20,37
Overexpression of miR-21 acts as an antiapoptotic factor in human glioblastoma cells.13
Inhibition of miR-21 inhibits cell growth in vitro and inhibits tumor growth in xenograft mouse models through an indirect down-regulation of the antiapoptotic factor, B-cell lymphoma 2 (Bcl-2).38
Studies in human cell lines have shown miR-21 can also target the tumor suppressor genes,
phosphatase and tensin homolog (PTEN)39
and tropomyosin 1 (TPM1).40 These data, taken together, support a causal role for altered microRNA expression during tumorigenesis.
Adenomas represent a precursor stage of adenocarcinoma. Adenomas express high levels of miR-21. If increased miR-21 expression promotes colon tumor progression,
increased expression in adenomas may be an early cellular event in the progression to cancer. Inhibiting miR-21 activity may help prevent tumor promotion in populations at high risk of colon cancer, such as individuals with familial adenomatous polyposis.41
In this study, we also demonstrated an association with microRNA expression patterns with colon cancer prognosis and therapeutic outcome.
A robust association of high miR-21 expression in tumors with poor survival was observed in the Maryland test cohort and the Hong Kong validation cohort, separately. In each cohort, these associations were independent of other clinical covariates indicating that miR-21 expression may be a useful prognostic indicator, in addition to TNM staging and other clinical parameters,
to help identify patients at a higher risk of terminal cancer. These observations were made in 2 independent cohorts with different racial and geographical compositions. Therefore, our observations may be broadly applicable to other populations.
High miR-21 expression in tumors was associated with a poor therapeutic outcome in both cohorts. This association may help predict the benefits of therapy in individuals whose miR-21 expression status is known and identify patients who are candidates for more aggressive initial therapies.
But, if high miR-21 expression is causal to the poor therapeutic outcome, antagomirs30,42
or other antisense therapeutics that target miR-21 may prove to have therapeutic benefits in patients with tumors with high expression of miR-21.
Additional studies are required to demonstrate a causal link with miR-21 and the progression of colon cancer to determine the potential of miR-21 as either a biomarker or therapeutic target. Although a causal role is still uncertain, many of Hill's criteria for causation43
have already been met. There is a temporal relationship such that a high expression of miR-21 in tumors precedes the progression, therapeutic response, and subsequent death due to cancer. The strength and consistency of these associations has been found in 2 independent cohorts. There is a dose-response relationship such that more advanced tumors express higher levels of miR-21 in both cohorts. Published research has demonstrated biological plausibility of miR-21 causing the progression of tumors using in vitro and in vivo models. All of these are consistent with a role for miR-21 in colon carcinogenesis.
In conclusion, we found systematic differences in microRNA expression patterns between colon tumors and paired nontumorous tissue. Tumors with high expression of miR-21 was associated with poor survival outcome and poor response to adjuvant chemotherapy in 2 independent cohorts, independent of staging and other clinical covariates suggesting that miR-21 may be a useful diagnostic biomarker for colon adenocarcinomas and survival prognosis including response to therapy.
Corresponding Author: Curtis C.
Harris, MD, Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, Bldg 37, Room 3068A, Bethesda,
MD 20892 (Curtis_Harris@nih.gov).
Author Contributions: Drs Schetter and Harris had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Schetter,
Leung, Yanaihara, Yuen, Harris.
Acquisition of data: Schetter, Leung,
Sohn, Bowman, Yuen, Chan, Kwong, Au, Liu, Calin, Croce, Harris.
Analysis and interpretation of data:
Schetter, Leung, Sohn, Zanetti, Liu, Croce, Harris.
Drafting of the manuscript: Schetter,
Sohn, Croce, Harris.
Critical revision of the manuscript for important intellectual content: Schetter, Leung, Zanetti, Bowman, Yanaihara,
Yuen, Chan, Kwong, Au, Liu, Calin, Harris.
Statistical analysis: Schetter, Zanetti,
Obtained funding: Leung, Yuen, Chan.
Administrative, technical, or material support:
Leung, Sohn, Bowman, Yuen, Chan, Kwong, Au, Liu.
Study supervision: Harris.
Financial Disclosures: None reported.
Funding/Support: This research was supported by the Intramural Research Program of the National Cancer Institute, Center for Cancer Research, National Institutes of Health.
Drs Leung, Chan, and Yuen were supported by the Research Grants Council of the Hong Kong Special Administrative Region (HKU 7498/06M). Drs Schetter and Zanetti were supported by the Cancer Prevention Fellowship Program, Office of Preventive Oncology, National Cancer Institute.
Role of the Sponsors: The funding agencies had no role in the design and conduct of the study; collection, management,
analysis, and interpretation of the data; and preparation, review,
or approval of the manuscript.
Additional Contributions: We thank Xin Wang, PhD, Stefan Ambs, PhD, Kensuke Kumamoto, MD, PhD, Ewy Mathe, PhD, and Rosemary Braun, PhD, for helpful discussion; Megan Hames,
BS, and Jin Jen, MD, PhD, for advice on in situ hybridization techniques;
Judith Welsh, BS, for technical assistance; and Dorothea Dudek-Creaven,
MSC, for editorial assistance, all of the Center for Cancer Research,
National Cancer Institute. We thank Kathleen Meyer, BA, Center for Information Technology, National Institutes of Health, for uploading microarray data files to Gene Omnibus Profiles. We also than Raymond Jones, PhD, Audrey Salabes, BS, and John Cottrell, MS, of the University of Maryland Medical Center, the clinicians of the departments of surgery and pathology in Queen Mary Hospital, at University of Maryland Hospital, and at the Baltimore Veterans Administration Medical Center for contributions to this study.
This article was corrected online for typographical errors on 10/17/2008.