microRNA expression is able to predict response to chemoradiotherapy in rectal cancer
- Authors:
- Published online on: August 9, 2012 https://doi.org/10.3892/mco.2012.9
- Pages: 137-142
Abstract
Introduction
Loco-regional recurrence subsequent to rectal cancer resection is difficult to treat and is associated with severe debilitating symptoms. The prognosis after a local recurrence is poor, with a median survival period of 12–18 months (1). Pre-operative chemoradiotherapy (CRT) has been widely used as the major treatment modality to improve local control of the disease, as well as to preserve anal sphincter (2). However, response to CRT differs in tumors.
microRNAs (miRNA) are a family of small non-coding RNA molecules that downregulate the expression of their protein-coding gene targets (3). Despite the relatively small number of miRNAs, computational and experimental studies have shown that a number of human protein-coding genes are collectively regulated by miRNAs (4–6). Thus, miRNAs are considered to be master regulators of several important biological processes, such as cell growth, apoptosis, viral infection and cancer development (3,7–10). miRNA expression profiles have been shown to be promising biomarkers for the classification or outcome prediction of a wide array of human cancers (11,12). While global miRNA expression patterns of several embryologic, physiologic and oncogenic processes have been thoroughly studied, no reports are available on the role of miRNAs in pre-operative CRT in rectal cancer. This is the first study using miRNA microarray for predicting response to CRT in rectal cancer. To predict response to CRT, the miRNA expression patterns of rectal cancer using miRNA microarray prior to pre-operative CRT were examined. In addition to miRNA microarray analysis, a candidate miRNA was examined using real-time reverse transcription-polymerase chain reaction (RT-PCR).
The present study aimed to define the expression patterns of miRNA for the prediction of response to CRT and to establish tailor-made therapies in rectal cancer.
Materials and methods
Patients and tissue samples
For miRNA expression profiling, rectal cancer samples were obtained from 43 patients, who were willing to receive pre-operative CRT between September 2007 and December 2009, at the Tokushima University Hospital. The 43 independent rectal tumor samples included 22 for training and 21 for testing the outcome prediction model, respectively. The patient characteristics and response to CRT are summarized in Table I. The study was approved by the ethics committee of the Tokushima University Hospital, while every patient provided written informed consent for their samples to be used. Biopsy specimens from rectal cancer were prospectively collected during colonoscopic examinations prior to pre-operative CRT. Parallel tumor specimens were formalin-fixed and paraffin-embedded for histologic examination, while other specimens were used for RNA extraction. Samples were used for RNA extraction when paralleled specimens contained at least 70% tumor cells. Samples were snap-frozen immediately in liquid nitrogen and stored at −80°C until RNA extraction.
The patients received CRT at a total dose of 4,000 cGy of pelvic irradiation, five times a week, with a daily fraction of 200 cGy utilizing a four-field technique. Radiation was delivered concomitantly with S-1, a novel oral fluoropyrimidine inhibitory for dihydropyrimidine dehydrogenase with a potent radiosensitizing property. The S-1 was administered on radiation days. Surgical treatment was performed 6–8 weeks after the completion of pre-operative CRT.
miRNA microarrays
Total RNA was isolated using the miRNeasy mini kit (Qiagen, Hilden, Germany). One hundred nanograms of each RNA sample were hybridized to Agilent Human miRNA Microarray v2.0, containing 821 miRNAs (G4470B; Agilent Technologies, Inc., Santa Clara, CA, USA). miRNA labeling, hybridization and washing were carried out according to the manufacturer's instructions. Images of hybridized microarrays were obtained with a DNA micro-array scanner (G2565BA; Agilent Technologies), while features were extracted using the AFE image analysis tool version A.9.5.3.1.
Real-time RT-PCR
Complementary DNA (cDNA) was synthesized from total RNA to quantify the expression levels of mature miR-223 using gene-specific primers, according to the TaqMan miRNA Assay protocol (Applied Biosystems, Carlsbad, CA, USA) in a 15 μl reaction volume with 10 ng of RNA template. RT 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, while 2 μl were used in a total reaction volume of 20 μl for relative quantification by RT-PCR using an Applied Biosystems 7500 Sequence Detection system. The 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 sec and 60°C for 1 min. Normalization was performed using the small nuclear RNA U6 (RNU6B; Applied Biosystems).
Data analysis
To identify genes that were differentially expressed in the two groups, the data sets were assigned to either responders or non-responders. Response to CRT was evaluated by three parameters [histopathological examination, Response Evaluation Criteria in Solid Tumors (RECIST) and downstaging], while miRNA microarray was analyzed in the three parameters. The histopathological examination of surgically resected specimens was based on a semi-quantitative classification system, as described in detail previously (13). Tumors were classified as responder when assigned to regression grade 2 or 3, and non-responder when assigned to grade 0 or 1. Tumors in RECIST were classified as responder when assigned to complete response (CR) or partial response (PR), and non-responder when stable disease (SD) or progressive disease (PD). Tumors in downstaging were classified as responder when assigned to ‘Yes’, and non-responder when ‘No’.
The samples of 43 patients were divided into a training set (22 samples) and a testing set (21 samples). In the miRNA microarray analysis to evaluate miRNA expression only training samples were used. The expression patterns were compared and the fold change value was determined to identify gene markers that best discriminated between responders and non-responders. Two-dimensional hierarchical clustering was applied to the log-transformed data with average-linkage clustering, with a standard correlation serving as the similarity metric for the discriminating genes identified as differentially expressed in responders and non-responders. Using real-time RT-PCR in a training set, a candidate miRNA detected by miRNA microarray analysis with regard to the histopathological examination of surgical specimens was evaluated. At the end. A candidate miRNA was then evaluated in a testing set using real-time RT-PCR with regard to the histopathological examination of surgical specimens.
Statistical analysis
Quantitative data were provided as the median (range). The statistical analysis was performed using statistical software (JMP 8.0.1., SAS Campus Drive, Cary, NC, USA). Comparisons of clinicopathological and PCR data were performed using the Fisher's exact, Chi-square and Mann-Whitney U tests, as appropriate. The expression patterns in the miRNA microarray were compared using unpaired t-tests (with Welch's correction for unequal variances). Receiver-operating characteristics (ROC) curves were established to evaluate the diagnostic value of candidate miRNA for differentiating between responders and non-responders. The statistical tests performed were two-sided and declared at the 5% significance level.
Results
miRNA expression patterns using microarrays in responders and non-responders
miRNA expression profiling was established using miRNA array in the training samples. No statistically significant difference was detected between the training and the testing sets with regards to clinicopathological factors, such as gender, age, histopathologicalal classification, pre-operative tumor stage and response to CRT. The patient characteristics and response to CRT are summarized in Table I. Of the 22 training samples, 15 were classified as responders and 7 as non-responders, according to the histopathologicalal examination of the surgical specimens. Initially, with regard to the histopathologicalal examination of surgically resected specimens, 2 genes differentially expressed at significant levels (P<0.05) in responders and non-responders (Table II) were identified. The 2 genes (miR-223 and miR-142-3p) showed a higher expression in responders compared to non-responders. Of the 22 training samples, 14 were classified as responders and 8 as non-responders, based on the RECIST. Nine genes were differentially expressed at significant levels in responders and non-responders, with regard to RECIST (Table III). One gene (miR-223) showed a higher, while 8 genes a lower expression in responders compared to non-responders. The list of under-expressed genes included miR-20b, miR-92a, let-7a*, miR-20a, miR-17*, miR-106a, miR-17 and miR-20a*. Additionally, of the 22 training samples, 13 were classified as responders and 9 as non-responders, based on downstaging. Three genes were differentially expressed at significant levels in responders and non-responders, with regard to downstaging (Table IV). The 3 genes (miR-223, miR-630 and miR-126*) showed a higher expression in responders compared to non-responders. Results of a hierarchical cluster analysis of the 12 genes detected by three parameters are shown in Fig. 1. Responders and non-responders were clustered into two distinct groups, except for four responder cases.
Table IIFrequent differentially expressed miRNAs in responders and non-responders with regard to histopathologicalal examination of surgically-resected specimens. |
Table IIIFrequent differentially expressed miRNAs in responders and non-responders with regard to RECIST. |
Table IVFrequent differentially expressed miRNAs in responders and non-responders with regard to downstaging. |
Evaluation of miRNA microarray data by real-time RT-PCR
miR-223, a candidate gene, showing a higher expression in responders compared to non-responders in the three parameters (histopathological examination, RECIST and downstaging) was chosen for the evaluation of miRNA microarray data using real-time RT-PCR. Initially, the miR-223 level was evaluated using real-time RT-PCR in a training set. The miR-223 level was significantly higher in responders compared to non-responders, with regard to the histopathological examination of surgical specimens (P<0.01) (Fig. 2). The miR-223 level in a testing set (21 independent rectal cancer samples) was evaluated. Of the 21 testing samples, 9 were classified as responders and 12 as non-responders, based on the histopathological examination of surgical specimens. The miR-223 level was higher in responders compared to non-responders (P<0.01) (Fig. 3). Additionally, ROC curve analyses showed that miR-223 might differentiate between responders and non-responders in 22 training samples with an area under the curve (AUC) of 0.768 [95% confidence internal (CI), 0.661–0.865]. At the cut-off value of 0.4 for miR-223, the sensitivity and the specificity in the 21 testing samples were 100 and 78.0%, respectively (Fig. 3).
Discussion
Although miRNA expression patterns have been applied when predicting the outcome of multiple cancers, few studies have been reported as yet on the application of miRNAs when predicting the response to CRT using pre-operative biopsy tissue samples in rectal cancer. The expression patterns of miRNAs for the prediction of response to CRT by miRNA microarray were defined, while a candidate miRNA (miR-223) was evaluated by real-time RT-PCR.
In order to determine whether miR-223 showed a high or low expression when predicting the response to CRT in rectal cancer, in their study Aslam et al(14) summarized the differential expression of miRNAs in colorectal cancer compared to healthy flanking tissue identified by various studies. The results of those studies suggested that the variability of individual miRNAs with tumor type and stage probably rendered the combination of miRNAs a reliable method for detecting cancer status. In this study, the miR-223 detected as a candidate gene for the prediction of response to CRT is reported as one of the overexpressed miRNAs in colorectal cancer compared to healthy tissue (15,16).
The potential value of miRNAs as prognostic and predictive biomarkers in colon cancer is demonstrated by Schetter et al(17), who compared miRNA expression patterns in stage II colonic adenocarcinoma and in adjacent healthy tissue. A high tumor:healthy expression ratio of miR-20a, miR-21, miR-106a, miR-181b and miR-203 was associated with poor survival. Underexpressed miR-20a in a tumor may be associated with good survival and responders to CRT.
Regarding the correlation between miRNAs and radiation therapy, higher rectal tumor expression levels of miR-125b and miR-137 are reported to be associated with a poorer response to CRT (18). In lung cancer and healthy lung cell lines, levels of 81 out of 440 miRNAs demonstrated statistically significant differences subsequent to irradiation with 23 miRNAs downregulated after treatment. The members of the let-7 family (except let-7g) decreased markedly by 2–8 h subsequent to irradiation in both cancer and healthy lung epithelium. Furthermore, let-60/RAS and genes in the DNA damage response pathway were identified as mechanisms potentially affected by let-7 (19). In addition to these results, radiotherapy of two different prostate cancer cell lines resulted in a differentiated expression of a high number of miRNAs. The expression of 48 miRNAs altered significantly subsequent to the application of radiotherapy, with 22 of them presenting a >3-fold change in levels. Of those 22 miRNAs hsa-miR-521, hsamiR-196a and hsa-miR-133b were demonstrated to be decreased, while hsa-miR-34c was shown to be increased in the two cell lines subsequent to radiation (20). Although inconsistent with the results of the present study, the data on the let-7 family are likely to be affected by the differences in cell culture experiments and clinical research.
Different miRNAs have been found to predict sensitivity to anticancer treatment. miR-30c, miR-130a and miR-335 have been reported to be downregulated in various chemoresistant cell lines (21), while the restoration of miR-34 in p53-deficient human gastric cancer cells induced chemosensitisation (22). Irrelevant to the results of this study, regarding the S-1 used for CRT in the present study, Nakajima et al(23) reported that let-7g and miR-181b were strongly associated with response to 5-fluorouracil-based antimetabolite S-1.
miR-223 has been reported to function as an important modulator of cellular differentiation (24,25). In bone cells, the enhanced expression of pre-miR-223 completely abrogated osteoclast differentiation (25). However, few reports are available on the dysregulation of miR-223 in epithelial cancers. Although miR-223 transcription in granulopoiesis has been reported to be regulated by an evolutionarily conserved system driven by the myeloid transcription factors, PU.1 and CAAT enhancer-binding protein (26), In their study, Wong et al(27) identified stathmin 1 (STMN1) as a putative target of miR-223 in hepatocellular carcinoma (HCC). STMN1 is a key microtubule-regulatory protein that controls the microtubule dynamics, cell proliferation and S-phase of the cell cycle (28,29). High-levels of STMN1 expression have been reported in leukemia, breast, ovarian and prostate cancers, in which its increased expression has been associated with increased histologic grading, shorter patient survival period and increased resistance to medication (30–33). By these mechanism, STMN1 as a putative target of miR-223, may be associated with responders to CRT.
Additionally, miR-223 has recently been found to be myeloid-specific, negatively regulating the progenitor proliferation, granulocyte differentiation and activation in mice (34). Guo et al(35) also reported that the CD133+ ovarian cancer stem cells in the OVACAR3 cell line were correlated with miR-223. miR-223, as a candidate miRNA to regulate stemness gene, may be associated with responders to CRT.
The present study defined the expression patterns of miRNAs for the prediction of response to CRT by miRNA microarray and evaluated miR-223 by real-time RT-PCR using pre-operative biopsy tissue samples in rectal cancer. Albeit the small samples, the lack of analysis of target mRNA, of miR-223 and functional analysis of miR-223, the expression patterns of miRNAs for the prediction of response to CRT by miRNA microarray were evaluated using several miRNAs (821 gene). Furthermore, miR-223 was evaluated by real-time RT-PCR as well as in independent samples. These results suggest that a candidate common miR-223 shown by three parameters (histopathological examination, RECIST and existence of downstaging) may be a new biomarker for the prediction of response to CRT in rectal cancer.
In conclusion, the candidate miR-223 is a potential new biomarker for the prediction of response to CRT and is likely to be useful for establishing tailor-made therapies for rectal cancer.
Acknowledgements
The present study was partly financed by the Research Support Foundation of the University of Tokushima and Taiho Pharmaceutical, Co., Ltd. (research funding provided to Professor Mitsuo Shimada), and the Grants-in-Aid for Scientific Research of the Japan Society for the Promotion of Science. A part of this study was presented at the Gastrointestinal Cancer Symposium, San Francisco, CA, USA (January 15–17, 2009).