Next Article in Journal
Asymmetric Dimethylarginine versus Proton Pump Inhibitors Usage in Patients with Stable Coronary Artery Disease: A Cross-Sectional Study
Next Article in Special Issue
Identification of More Feasible MicroRNA–mRNA Interactions within Multiple Cancers Using Principal Component Analysis Based Unsupervised Feature Extraction
Previous Article in Journal
Hybrid Macro-Porous Titanium Ornamented by Degradable 3D Gel/nHA Micro-Scaffolds for Bone Tissue Regeneration
Previous Article in Special Issue
miRTargetLink—miRNAs, Genes and Interaction Networks
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 17
Issue 4
10.3390/ijms17040568
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prognostic Value of MicroRNAs in Preoperative Treated Rectal Cancer

by
Azadeh Azizian
1,
Ingo Epping
1,
Frank Kramer
2,
Peter Jo
1,
Markus Bernhardt
1,
Julia Kitz
3,
Gabriela Salinas
4,
Hendrik A. Wolff
5,
Marian Grade
1,
Tim Beißbarth
2,
B. Michael Ghadimi
1 and
Jochen Gaedcke
1,*
1
Department of General, Visceral, and Pediatric Surgery, University Medical Center Göttingen, Robert-Koch-Str. 40, Göttingen 37075, Germany
2
Department of Medical Statistics, University Medical Center Göttingen, Robert-Koch-Str. 40, Göttingen 37075, Germany
3
Department of Pathology, University Medical Center Göttingen, Robert-Koch-Str. 40, Göttingen 37075, Germany
4
Department of Developmental Biochemistry, University of Göttingen, Göttingen 37075, Germany
5
Medical Practice Radiotherapy München, Burgstraße 7, München 80331, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(4), 568; https://doi.org/10.3390/ijms17040568
Submission received: 9 February 2016 / Revised: 1 April 2016 / Accepted: 6 April 2016 / Published: 15 April 2016
(This article belongs to the Special Issue MicroRNA Regulation)
Browse Figures
Versions Notes

Abstract

:
Background: Patients with locally advanced rectal cancer are treated with preoperative chemoradiotherapy followed by surgical resection. Despite similar clinical parameters (uT2-3, uN+) and standard therapy, patients’ prognoses differ widely. A possible prediction of prognosis through microRNAs as biomarkers out of treatment-naïve biopsies would allow individualized therapy options. Methods: Microarray analysis of 45 microdissected preoperative biopsies from patients with rectal cancer was performed to identify potential microRNAs to predict overall survival, disease-free survival, cancer-specific survival, distant-metastasis-free survival, tumor regression grade, or nodal stage. Quantitative real-time polymerase chain reaction (qPCR) was performed on an independent set of 147 rectal cancer patients to validate relevant miRNAs. Results: In the microarray screen, 14 microRNAs were significantly correlated to overall survival. Five microRNAs were included from previous work. Finally, 19 miRNAs were evaluated by qPCR. miR-515-5p, miR-573, miR-579 and miR-802 demonstrated significant correlation with overall survival and cancer-specific survival (p < 0.05). miR-573 was also significantly correlated with the tumor regression grade after preoperative chemoradiotherapy. miR-133b showed a significant correlation with distant-metastasis-free survival. miR-146b expression levels showed a significant correlation with nodal stage. Conclusion: Specific microRNAs can be used as biomarkers to predict prognosis of patients with rectal cancer and possibly stratify patients’ therapy if validated in a prospective study.

1. Introduction

Colorectal cancer is the third most common cancer worldwide accounting for 1.36 million newly diagnosed cases in 2012 [1]. About 30% of all colorectal cancers are rectal cancers. The treatment of rectal cancer patients depends on the stage of the disease. While patients with Union for International Cancer Control (UICC) stage I and II receive a primary surgery, others with locally advanced rectal cancer are treated first with preoperative chemoradiotherapy (CRT) followed by surgery (total mesorectal excision (TME)). However, within the respective UICC stages, the prognosis of patients may differ. This is partly due to the fact that the UICC staging relies mostly on imaging techniques that are not always able to discriminate appropriately. In addition, patients´ response to the preoperative treatment is heterogeneous. An adequate prediction of response and prognosis in addition to the UICC staging could serve to adapt the therapy modalities individually within the UICC stages and also to individualize the follow-up care to the predicted prognosis. Patients with poor prognosis could receive a more aggressive CRT (e.g., additional oxaliplatin) possibly accepting more severe side effects in order to achieve a better outcome. Furthermore, they could benefit from more frequent follow-up visits to detect early possible disease recurrence as defined by the occurrence of local recurrence or distant metastasis. To predict the prognosis of patients, reliable biomarkers are needed. These can be, for instance, specific DNA mutations, RNA sequences, or proteins. While proteins are subject to many post-transcriptional modifications and DNA sequences do not necessarily reflect the gene expression, some RNA molecules are promising candidates: in contrast to mRNA and tRNA, mature microRNAs (miRNAs) are highly conserved between species [2], very robust towards temperature and pH-values [3]. Due to their apparent proximity to chromosomal breakpoints [4] and their dysregulated expression levels in many malignancies [5], miRNAs were linked to tumorigenesis early after their discovery.
miRNA are small single-stranded non-coding RNAs, about 20–22 nucleotides in length, which play an important role in the post-transcriptional regulation of mRNA [6]. They bind to the targeting mRNA sequence within a complex of specific proteins, called RNA-induced silencing complex (RISC), and induce a cleavage or repression of the target gene [7]. Therefore miRNAs play an important role in the gene expression, and influence many physiological and pathophysiological processes. These qualities make them promising candidates for biomarkers in tumor specimen and liquid biopsies. In 2003, Michael et al., found, for the first time, miRNAs being associated with colorectal cancer [8]. Since then, several studies used miRNAs trying to predict different clinical parameters [9,10,11]. In 2012, we identified a rectal cancer specific miRNA-panel characterizing differences between tumors and normal adjacent mucosa [12]. Eleven of the identified microRNAs were also found by Li et al. [13] in the same year. While many following studies investigated a specific miRNA-panel for colorectal cancer and others tried to find miRNAs specific for each UICC stage, a large part of research was investigating miRNAs for predicting response to CRT [14,15,16,17].
Prognosis of patients depends on several known factors: UICC stage, tumor regression grade (TRG), nodal stage, and surgical margins. A large part of these factors, namely TRG, nodal stage, and the quality of the surgical margins, is not known until after CRT and surgery. A that point, the large part of the treatment is already performed and possible side-effects of the preoperative CRT and surgery can not be undone. Therefore, these factors can only have little impact on individualizing the therapy. Our aim is to predict patients’ prognoses in advance, before any treatment. This way, we would be able to stratify the therapy in a way patients would benefit from the most. In the present work, we aimed to explore the impact of miRNAs as biomarkers to predict the patients’ prognosis and response to CRT analyzed in biopsies, which were taken prior to any treatment. All patients were enrolled or treated according to the CAO/ARO/AIO-94 [18,19,20] and CAO/ARO/AIO-04 [21,22] trial of the German Rectal Cancer Study Group.

2. Results

2.1. Microarray Analysis Identified 14 miRNAs Significantly Associated with Overall Survival

First, we performed microarray analyses of 45 microdissected pretherapeutic biopsies from patients with rectal cancer to identify potential microRNAs with a prognostic value. The clinical data of the patients, including gender, age, UICC stage, cancer-specific survival (CSS), local recurrence (LR), distant-metastasis-free survival (DMS), and disease-free survival (DFS) are summarized in Supplementary Table S1A.
Time-to-event analyses were performed on a gene-by-gene basis, associating the survival times of patients with the expression level of each feature on the microarray chips using Cox Proportional Hazard Ratio. For 14 miRNAs, we could find a significant association to at least two of four survival parameter (overall survival (OS), DFS, CSS or DMS) with p < 0.05 (illustrated in Supplementary Figure S1). Those were chosen for further validation as they were considered as promising candidates with probable prognostic or predictive value in rectal cancer patients due to their expression level.
Further, this list of miRNAs was supplemented by five miRNAs (miR-198, miR-223, miR-320a, miR-34b, and miR-497), which showed a possible predictive value in rectal cancer in previous unpublished microarray analysis results of our group and literature research [12,23,24]. In order to further validate our findings, we collected 147 samples and tested the 19 miRNAs using qPCR.

2.2. Expression Levels of 19 miRNAs Were Analyzed in 147 Samples: miR-515-5p, miR-573, miR-579, and miR-802 Were Significantly Correlated to Overall Survival and Cancer-Specific Survival

The expression levels of the 19 miRNAs (listed in Table 1) in biopsies of rectal cancer tumor tissue (n = 147 patients, one sample per patient) were analyzed via qPCR and compared with the clinical parameters OS, CSS, DFS, DMS, and postoperative nodal stage (ypN). Four miRNAs, namely miR-515-5p, miR-573, miR-579 and miR-802, were associated significantly with OS and CSS (p < 0.05). Of these four miRNAs, only miR-573 was also associated significantly with the TRG (p = 0.0416), which has a known association to the survival of patients [25]. miR-515-5p, miR-573, miR-579, and miR-802 were not able to discriminate between good and poor DFS, DMS, or ypN. All p-values of the analyzed miRNAs are summarized in Table 1. Figure 1 shows the survival curves miR-515-5p, miR-573, miR-579, and miR-802 separately. Figure 2 shows the discrimination between TRG 1-4 for miR-573.

2.3. miR-133b Is Significantly Associated with Distant-Metastasis-Free Survival

For miR-133b expression in treatment-naïve rectal cancer tumor biopsies, a significant association with DMS has been shown (p = 0.032). Patients with a low expression level of miR-133b develop more frequently distant-metastasis. miR-133b is the only miRNA (among the investigated miRNAs in this study) being significantly associated to distant-metastasis-free survival, while it does not show any significant association with other clinical parameters according to our analysis. Figure 3 shows distant-metastasis-free survival dependent on a high and low expression of miR-133b in the tumor biopsies, respectively.

2.4. miR-146b Expression Levels Show a Significant Association with a Negative Post-Therapeutic Nodal Stage (ypN0)

While miR-146b was not associated with OS, DFS, CSS, DMS, or TRG, it showed a significant association with a negative nodal stage after preoperative CRT and surgery, assessed by the pathologist. However, since we do not know the nodal stage before CRT exactly—this can only be estimated by imaging techniques and radiologists—there is no comparison possible. Therefore, the association between the miR-146b expression and ypN does not make any statement concerning the response to CRT.

3. Discussion

Since the treatment of locally advanced rectal cancer includes a preoperative CRT, to which patients respond differently, the aim to find valid biomarkers to predict response in advance has driven many studies. Not only the response towards CRT is heterogeneous but also the prognoses of patients within the UICC stages. Possible biomarkers in blood or cancer biopsies taken prior to any therapy may help to stratify the therapy or adapt the postoperative follow-up examinations in terms of frequency and invasively. Here, miRNAs are one of several possible approaches: Different possible biomarkers were associated with therapy response: fibroblast growth factor receptor 2 (FGFR2) [26], β-catenin, vascular endothelial growth factor, and apoptotic protease activating factor 1 [27], expression of DNAJC12 [28], altered DNA methylation [29], and neutrophil-lymphocyte ratio [30]. Concerning survival prognosis, most studies rely on clinical parameters [31,32,33], but also some molecular parameters were evaluated [34,35,36,37,38].
In the present work, we showed that miR-515-5p, miR-573, miR-579, and miR-802 were correlated significantly with overall survival and cancer-specific survival (p < 0.05) of patients with locally advanced rectal cancer treated with a preoperative CRT followed by surgery with quality controlled nerve-sparing TME (according to the study protocol).
These four miRNAs are expressed independently from each other, meaning they do not belong to a common cluster. Although miR-515-5p, located on chromosome 19, belongs to a cluster with 42 members, none of the other three miRNAs (miR-573, miR-579 or miR-802) belongs to that cluster. miR-515 is yet investigated only in few studies and has not been reported in rectal cancer: in 2010, miR-515 was one of 13 miRNAs found overexpressed in oral carcinoma [39]. Follicular tumors showed a distinct overexpression of several members of miR-515 family [40]. These studies together with our work suggest a tumor promoting role of miR-515. But a recent study showed that miR-515-5p inhibits cancer cell migration and metastasis by targeting microtubule affinity regulating kinase 4 (MARK4) in breast cancer [41]. A different breast cancer study revealed IGF-1R gene as another target for 515-5p [42].
Recently, miR-573 is found to down-regulate the oncogene Tetraspanin 1 (TSPAN1) in gastric cancer [43]. Also, in a study concerning breast cancer, miR-573 was found to be down-regulated in BRCA 1/2-related breast cancer [44]. Wang et al. found that miR-573 targets the melanoma cell adhesion molecule (MCAM) and regulates this way melanoma progression [45]. They showed in vivo using a murine xenograft model that overexpression groups had lower rates of tumor growth compared with the control group. Taking these studies together, miR-573 seems to perform a tumor suppressive function in diverse cancers. However the results of our work, where we investigated miR-573 for the first time in rectal cancer, show a higher expression level of miR-573 associated with poor overall survival and cancer-specific survival (see Figure 1), suggesting a tumor promoting role. This might show tissue specific effects of miR-573, which can be oppositional in different tumor entities.
The role of miR-579 in cancer cells is poorly investigated yet. Quinn et al. concluded in their review about the role of miRNA in regulation of endotoxin tolerance that miR-579 plays a role in regulating the TLR4 signaling pathway during the development of endotoxin tolerance at receptor, signaling pathway, gene transcription and translational levels [46]: miR-579 blocks the translation of TNF-α [47]. According to our results, a high expression of miR-579 in rectal cancer tissue is significantly associated with a poor overall survival and cancer-specific survival (see Figure 1). Further investigations concerning target genes or function in cancer cells needs to be done for miR-579.
Being located on chromosome 21, miR-802 is consequently of great interest in down syndrome associated research [48,49], where it is shown to target the methyl-CpG-binding protein (MeCP2), which is underexpressed in down syndrome brains [50]. Also, an association between miR-802 and diabetes was found by Kornfeld et al., who showed that obesity-induced overexpression of miR-802 compromises glucose metabolism through silencing of Hnf1b [51]. Nevertheless, it also seems to be of importance in several cancer types: It is supposed to suppress breast cancer proliferation by targeting and suppressing FoxM1 [52] and showed a moderate tumor suppressive activity in liver cancer cells [53]. On the other hand, Cao et al., claim that miR-802 promotes osteosarcoma cell proliferation by targeting p27 [54]. Another study shows that miRNA‑802 targets the tumor suppressor menin and promotes lung carcinoma proliferation [55]. According to our work, a high expression of miR-802 in rectal cancer tissue is associated significantly with a poor overall survival and cancer-specific survival (Figure 1).
We showed a significant association of miR-133b expression in rectal cancer tumor tissue with a distant-metastasis-free survival of the patients. miR-133b is located on chromosome 6, with the seed sequence TTGGTCC, which is unique and not shared (as far as miRNAs are sequenced till today) between other miRNAs in a so called “seed-family”. In addition, miR-133b does not belong to any bigger cluster of microRNAs (as far as present studies demonstrate), indicating a unique role of this miRNA. miR-133b is declared to have tumor suppressor functions in several studies [56,57,58]. The association of a high expression of miR-133b in rectal cancer with good distant-metastasis-free survival, as shown in our data, is supported by different studies: Qiu et al. showed that miR-133b (and to other miRNAs) suppressed the proliferation, migration, invasion and cell cycle progression in gastric cancer cells through decreasing expression of the transcription factor specificity protein 1 (Sp1) and its downstream proteins using human gastric cancer cell lines SGC7901, MKN45 and BCG823 [59]. Also, Zhao et al. showed for gastric cancer that miR-133b is frequently down-regulated in these tumor biopsies and that its overexpression in gastric cancer cell lines reduces the metastatic potential [60]. They identified the transcriptional factor Gli1 as a direct target of miR-133b: level of Gli1 protein but not mRNA was decreased by miR-133b. In colon cancer cell lines (HT-29 and SW-620) miR-133b was shown be essential for the inhibitory effects of TAp63 on RhoA, E-cadherin and vimentin [61], the authors concluded that miR-133b is able to suppress the metastasis of colon cancer. Also in non-cancer cells, miR-133b seems to correlate inversely to proliferation markers (e.g., Ki67 and cleaved-CK18) suggesting an anti-proliferative role for miR-133b [62]. This data along with our results reveal that miR-133b has an anti-proliferative effect in general, which can be shown in cancer but also in cardiac diseases [63], muscle regeneration [64,65] and embryonic development.
For miR-146b we showed a significant association with the post-CRT nodal stage (ypN) diagnosed by experienced pathologists. A high miR-146b expression is associated with a negative post-CRT nodal stage. According to several studies miR-146b seems to play an important role papillary thyroid carcinoma: Geraldo et al. [66] showed a significantly poorer overall survival for patients with higher miR-146b expression levels compared to those with lower miR-146b levels. Using BCPAP, a human papillary thyroid cancer cell line model, they showed that an overexpression of miR-146b significantly increases cell migration and invasiveness, along with increased resistance to chemotherapy-induced apoptosis [66]. Hardin et al. [67] analyzed plasma-derived miRNAs in patients with human papillary carcinoma and showed a higher level of miR-146b in patients compared to those with benign lesions. Furthermore, they also showed that the levels of miR-146b increased in proportion to tumor size [67]. Other studies confirmed its role in thyroid cancer [68,69,70,71]. In other cancer types, miR-146b seems to be involved: Li et al. reported that miR-146b promotes cell proliferation and inhibits cell apoptosis in esophageal cancer cell lines [72], whereas in glioma miR-146b also seems to perform a tumor suppressor function [73,74].
Other studies have also analyzed the value of miRNAs to predict prognosis or therapy response. miR-21-5p was shown to be overexpressed in complete responders [17]. However, the cohort was rather small (n = 43). Another study investigated miR-21 for the same purpose but could not confirm the predictive value of miR-21 [75]. High levels of miR-125b and miR-137 in tumor tissue are shown to be associated with a worse response to CRT [15], although in a small cohort (n = 31). Also miR-16, miR-590-5p and miR-153 were associated with therapy response by Kheirelseid et al. [76]. Other studies identified miRNA profiles that seem to correlate with response to neoadjuvant chemoradiotherapy [77]. A possible reason for the different identified miRNAs may be the difference between the tissues preservation methods (formalin-fixed paraffin-embedded (FFPE), RNAlater, fresh frozen biopsies) and analyzing methods (qPCR, microarray). Several studies focused on miRNAs and response to targeted anti epidermal growth factor receptor (anti-EGFR) therapy, which is a well-established therapy in KRAS-wildtype metastatic CRC. However, even in KRAS-wildtype patients only approximately 40% respond to this therapy. miR-31 and miR-181 were associated with response to anti-EGFR therapy [78,79]; a let-7 microRNA-binding site polymorphism in KRAS, the let-7 complementary site 6 (LCS6)-variant, also seems to have an impact on response to anti-EGFR therapy, though in this regard the study results remain contradictory [80,81].
Taken our results together, we identified six miRNAs with a predictive value concerning patients’ prognosis. High expression levels of miR-515-5p, miR-573, miR-579 and miR-802 showed all a significant correlation to a poor OS and CSS. These four miRNAs could serve together as a miRNA profile to identify patients with a poor outcome, if validated prospectively.
We used two different types of preservation (FFPE for the first set of 45 patients and RNAlater for the validation set with 147 patients), to overcome biases arising from the one or the other type of preservation. Also, we used microarray analysis first and than qPCR to confirm our results. Thus, the identified miRNAs showed their prognostic value in both FFPE and RNAlater, and also in both microarray analysis and qPCR, which strengthens their impact.
However, the identified miRNAs seems to have an impact in other cancer types too, therefore none of them is a rectal cancer specific miRNA. Besides, in some cases, the investigated miRNAs play more of an oppositional role in rectal cancer than in other cancer types. Furthermore, for some of these miRNAs the target genes and functions still need to be investigated.

4. Experimental Section

4.1. Patients and Samples

Pretherapeutic tumor samples were collected from patients with locally advanced rectal cancer, who were treated at the Department of General, Visceral and Pediatric Surgery, University Medical Center, Göttingen, Germany. For the genome-wide screening, 45 patients were enrolled. These patients were identified in the database of the CAO/ARO/AIO-94 [18,19,20] and CAO/ARO/AIO-04 [21,22] trial of the German Rectal Cancer Study Group, and selected based on the criteria: locally advanced tumor, preoperative CRT, R0 resection, follow-up of at least 150 days, and availability of a pretherapeutic Formalin-Fixed Paraffin-Embedded biopsy (FFPE) showing invasive cancer. These patients were then split into two comparable groups based on the occurrence of distant metastases. For validation, 147 patients from the prospective collection were identified. Inclusion criteria were comparable to the screening set, but here the biomaterial used was preserved in RNAlater (in contrast to the biomaterial used for microarray analysis in the first set, which was FFPE material). Patients that were already analyzed in the screening set were excluded from validation analysis. Figure 4 illustrates the study design.
According to the CAO/ARO/AIO-94 [18,19,20] and CAO/ARO/AIO-04 [21,22] trial of the German Rectal Cancer Study Group, they received a preoperative CRT with a total radiation dose of 50.4 Gy (single dose of 1.8 Gy) accompanied by either an intravenous (iv) application of 5-FU or a combination of an iv-infusion of oxaliplatin and 5-FU according to the study protocol. Four to six weeks after preoperative CRT, surgery with quality controlled nerve-sparing TME was performed. Tumor regression grade (TRG) and posttherapeutic nodal stage (ypN) were assessed by experienced pathologists at the University Medical Center Göttingen, Germany. Written informed consent was obtained from all patients according to the guidelines approved by the local ethic committee respecting the principles of the Declaration of Helsinki. The clinical data are summarized in Supplementary Table S1. The follow-up period accounts for minimum 158 days, maximum 4111 days and mean 1620 days.

4.2. Tumor Biopsies, RNA Isolation and Microarray Analysis

Tumor biopsies, taken during the index rectoscopy, were immediately buffered formalin for the first set (n = 45); those for the validation set (n = 147) were transferred into RNAlater (Qiagen, Hilden, Germany). Tissue in RNAlater was stored over night at 4 °C to allow saturation of the entire biopsy. The following day, they were stored at −20 °C for long-term storage. Biopsies for FFPE tissue were stored overnight in 4% formalin and were processed the next day. RNA isolation from FFPE samples was carried out using the Qiagen AllPrep® DNA/RNA FFPE kit (Qiagen, Hilden, Germany) according to the manufactures protocol after microdissection of the tumor tissue as previously described. RNA extraction from RNAlater biopsies was performed using TRIZOL (cat. no. 15596-018, Life Technologies, Rockeville, MD, USA) as previously described [12]. For these samples nucleic acid quantity, quality, and purity was determined using a spectrophotometer (NanoDrop ND1000, Thermo Fisher Scientific Inc., Waltham, MA, USA) and a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Expression profiling for miRNAs was carried out on the 60K Human microRNA Microarray (Agilent Technologies) as previously performed [82].

4.3. Semi-Quantitative Real-Time PCR

In an independent validation set of 147 patients, tumor samples were collected prospectively and processed. RNA was isolated as described above. Using the miScript Reverse Transcription Kit (Qiagen, Hilden, Germany) cDNA was generated from total RNA. Forward primers for miR-133b, miR-146b, miR-198, miR-223, miR-224, miR-23b, miR-3133, miR-320a, miR-34b, miR-3941, miR-4263, miR-450b-3p, miR-497, miR-515-5p, miR-518f, miR-573, miR-579, miR-612, and miR-802 were obtained from Qiagen’s database. For reverse primer, we used the miScript Universal Primer (Qiagen). Semi-quantitative real-time PCR was performed using the CFX 384™-Real-Time System (BIO-RAD, Hercules, CA, USA). The Quantification was carried out using QuantiTect SBYR Green PCR Master Mix (Qiagen, Hilden, Germany). Probes were analyzed in triplicates and assays were performed according to the manufacturer’s instruction.

4.4. Statistical Analysis

The miRNA array expression levels were analyzed using log2 transformation and quantile normalization [83]. Except for control spots, all 807 features were used without any a priori filtering. The data were analyzed with regard to clinical parameters overall survival (OS), disease-free survival (DFS), distant-metastasis-free survival (DMS), cancer-specific survival (CSS), TRG, and post-therapeutic nodal stage (ypN). For 14 miRNAs, we could find a significant association to at least two of four survival parameter (OS, DFS, CCS or DMS) with p < 0.05.
The miRNA PCR values of 19 miRNAs of 147 patients were normalized using housekeepers miR-202, miR-874 and small nuclear RNA U44 (SNORD44). These miRNAs were selected from the genome-wide screening and showed the least level of variability over all samples and represented miRNAs of high, medium, and low expression. Survival analysis of this data was performed as well as analyses of variance in order to assess two linear models of miRNA expression levels and clinical parameters [84]. These linear models included the expression level as dependent variable and the independent variables nodal state and tumor regression grade respectively. A significant outcome with p < 0.05 suggests an association of miRNA expression and the respective clinical data.
Survival data were visualized using Kaplan–Meier plots and the effect of the individual miRNAs on survival was assessed using Cox proportional hazards regression [85]. A Cox regression model was calculated miRNA-wise, correlating miRNA expression levels and time to event data. The overall survival was calculated as time from surgery until death. Cancer-specific survival time was measured from the date of resection to the date of death due to rectal cancer. Disease-free survival was calculated as the time from surgery to local recurrence or distant metastasis. Finally, the distant metastasis–free survival is the period until metastasis is detected. For visualization on Kaplan-Meier plots, patients were grouped depending whether they had an expression level above or below median expression level for a particular miRNA.
All analyses were performed using the free statistical software R (version 3.1, R Core Team, Wien, Austria) [86]. Linear models were computed using the “limma” package. Survival analysis was conducted using the R package “survival”. p-values < 0.05 were considered significant. In order to not exceed a false discovery rate of 5%, p-values were adjusted for multiple testing using the Benjamini–Hochberg method [87].

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/17/4/568/s1.

Acknowledgments

This work was supported by the Deutsche Forschungsgemeinschaft KFO 179 and the Else Kröner Forschungskolleg. We would like to gratefully and sincerely thank Chang-Rong Lai and Jessica Eggert for their technical support. We did not receive any funds for covering the costs to publish in open access. We acknowledge support by the German Research Foundation and the Open Access Publication Funds of the Göttingen University.

Author Contributions

Azadeh Azizian, Jochen Gaedcke, Marian Grade and Peter Jo conceived and designed the experiments; Ingo Epping, Peter Jo and Markus Bernhardt performed the experiments; Gabriela Salinas, Azadeh Azizian, Frank Kramer and Tim Beißbarth analyzed the data; Julia Kitz analyzed the pathological results; and Hendrik A. Wolff and B. Michael Ghadimi contributed samples, reagents, materials and analysis tools. All authors wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CRTchemoradiotherapy
TMEtotal mesorectal excision
TRGtumor regression grade
miRNAmicroRNA
OSoverall-survival
CSScancer-specific-survival
LRlocal recurrence
DMdistant metastasis
DMSdistant-metastasis-free survival
DF disease-free-survival
ypNpost-chemoradiotherapy nodal status
ffemale
mmale

References

  1. Ferlay, J.; Steliarova-Foucher, E.; Lortet-Tieulent, J.; Rosso, S.; Coebergh, J.W.W.; Comber, H.; Forman, D.; Bray, F. Cancer incidence and mortality patterns in Europe: Estimates for 40 countries in 2012. Eur. J. Cancer 2013, 49, 1374–1403. [Google Scholar] [CrossRef] [PubMed]
  2. Ambros, V. MicroRNA pathways in flies and worms: Growth, death, fat, stress, and timing. Cell 2003, 113, 673–676. [Google Scholar] [CrossRef]
  3. Li, Y.; Jiang, Z.; Xu, L.; Yao, H.; Guo, J.; Ding, X. Stability analysis of liver cancer-related microRNAs. Acta Biochim. Biophys. Sin. 2011, 43, 69–78. [Google Scholar] [CrossRef] [PubMed]
  4. Calin, G.A.; Sevignani, C.; Dumitru, C.D.; Hyslop, T.; Noch, E.; Yendamuri, S.; Shimizu, M.; Rattan, S.; Bullrich, F.; Negrini, M.; et al. Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc. Natl. Acad. Sci. USA 2004, 101, 2999–3004. [Google Scholar] [CrossRef] [PubMed]
  5. Lu, J.; Getz, G.; Miska, E.A.; Alvarez-Saavedra, E.; Lamb, J.; Peck, D.; Sweet-Cordero, A.; Ebert, B.L.; Mak, R.H.; Ferrando, A.A.; et al. MicroRNA expression profiles classify human cancers. Nature 2005, 435, 834–838. [Google Scholar] [CrossRef] [PubMed]
  6. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef]
  7. Lai, E.C. Micro RNAs are complementary to 3′ UTR sequence motifs that mediate negative post-transcriptional regulation. Nat. Genet. 2002, 30, 363–364. [Google Scholar] [CrossRef] [PubMed]
  8. Michael, M.Z.; O’ Connor, S.M.; van Holst Pellekaan, N.G.; Young, G.P.; James, R.J. Reduced accumulation of specific microRNAs in colorectal neoplasia. Mol. Cancer Res. 2003, 1, 882–891. [Google Scholar] [PubMed]
  9. Xu, X.-H.; Wu, X.-B.; Wu, S.-B.; Liu, H.-B.; Chen, R.; Li, Y. Identification of miRNAs differentially expressed in clinical stages of human colorectal carcinoma—An investigation in Guangzhou, China. PLoS ONE 2014, 9, e94060. [Google Scholar] [CrossRef] [PubMed]
  10. Slattery, M.L.; Herrick, J.S.; Mullany, L.E.; Valeri, N.; Stevens, J.; Caan, B.J.; Samowitz, W.; Wolff, R.K. An evaluation and replication of miRNAs with disease stage and colorectal cancer-specific mortality. Int. J. Cancer 2015, 137, 428–438. [Google Scholar] [CrossRef] [PubMed]
  11. Nielsen, B.S.; Jørgensen, S.; Fog, J.U.; Søkilde, R.; Christensen, I.J.; Hansen, U.; Brünner, N.; Baker, A.; Møller, S.; Nielsen, H.J. High levels of microRNA-21 in the stroma of colorectal cancers predict short disease-free survival in stage II colon cancer patients. Clin. Exp. Metastasis 2011, 28, 27–38. [Google Scholar] [CrossRef] [PubMed]
  12. Gaedcke, J.; Grade, M.; Camps, J.; Søkilde, R.; Kaczkowski, B.; Schetter, A.J.; Difilippantonio, M.J.; Harris, C.C.; Ghadimi, B.M.; Møller, S.; et al. The rectal cancer microRNAome--microRNA expression in rectal cancer and matched normal mucosa. Clin. Cancer Res. 2012, 18, 4919–4930. [Google Scholar] [CrossRef] [PubMed]
  13. Li, X.; Zhang, G.; Luo, F.; Ruan, J.; Huang, D.; Feng, D.; Xiao, D.; Zeng, Z.; Chen, X.; Wu, W. Identification of aberrantly expressed miRNAs in rectal cancer. Oncol. Rep. 2012, 28, 77–84. [Google Scholar] [PubMed]
  14. Azizian, A.; Kramer, F.; Jo, P.; Wolff, H.A.; Beissbarth, T.; Skarupke, R.; Bernhardt, M.; Grade, M.; Ghadimi, B.M.; Gaedcke, J. Preoperative prediction of lymph node status by circulating miR-18b and MiR-20a during chemoradiotherapy in patients with rectal cancer. World J. Surg. 2015, 39, 2329–2335. [Google Scholar] [CrossRef] [PubMed]
  15. Svoboda, M.; Izakovicova Holla, L.; Sefr, R.; Vrtkova, I.; Kocakova, I.; Tichy, B.; Dvorak, J. Micro-RNAs miR125b and miR137 are frequently upregulated in response to capecitabine chemoradiotherapy of rectal cancer. Int. J. Oncol. 2008, 33, 541–547. [Google Scholar] [CrossRef] [PubMed]
  16. Svoboda, M.; Sana, J.; Fabian, P.; Kocakova, I.; Gombosova, J.; Nekvindova, J.; Radova, L.; Vyzula, R.; Slaby, O. MicroRNA expression profile associated with response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer patients. Radiat. Oncol. 2012, 7. [Google Scholar] [CrossRef] [PubMed]
  17. Lopes-Ramos, C.M.; Habr-Gama, A.; Quevedo, B.S.; Felício, N.M.; Bettoni, F.; Koyama, F.C.; Asprino, P.F.; Galante, P.A.; Gama-Rodrigues, J.; Camargo, A.A.; et al. Overexpression of miR-21–5p as a predictive marker for complete tumor regression to neoadjuvant chemoradiotherapy in rectal cancer patients. BMC Med. Genom. 2014, 7. [Google Scholar] [CrossRef] [PubMed]
  18. Sauer, R.; Becker, H.; Hohenberger, W.; Rödel, C.; Wittekind, C.; Fietkau, R.; Martus, P.; Tschmelitsch, J.; Hager, E.; Hess, C.F.; et al. German Rectal Cancer Study Group Preoperative versus postoperative chemoradiotherapy for rectal cancer. N. Engl. J. Med. 2004, 351, 1731–1740. [Google Scholar] [CrossRef] [PubMed]
  19. Sauer, R.; Liersch, T.; Merkel, S.; Fietkau, R.; Hohenberger, W.; Hess, C.; Becker, H.; Raab, H.-R.; Villanueva, M.-T.; Witzigmann, H.; et al. Preoperative versus postoperative chemoradiotherapy for locally advanced rectal cancer: Results of the German CAO/ARO/AIO-94 randomized phase III trial after a median follow-up of 11 years. J. Clin. Oncol. 2012, 30, 1926–1933. [Google Scholar] [CrossRef] [PubMed]
  20. Fokas, E.; Liersch, T.; Fietkau, R.; Hohenberger, W.; Beissbarth, T.; Hess, C.; Becker, H.; Ghadimi, M.; Mrak, K.; Merkel, S.; et al. Tumor regression grading after preoperative chemoradiotherapy for locally advanced rectal carcinoma revisited: Updated results of the CAO/ARO/AIO-94 trial. J. Clin. Oncol. 2014, 32, 1554–1562. [Google Scholar] [CrossRef] [PubMed]
  21. Rödel, C.; Liersch, T.; Becker, H.; Fietkau, R.; Hohenberger, W.; Hothorn, T.; Graeven, U.; Arnold, D.; Lang-Welzenbach, M.; Raab, H.-R.; et al. German Rectal Cancer Study Group Preoperative chemoradiotherapy and postoperative chemotherapy with fluorouracil and oxaliplatin versus fluorouracil alone in locally advanced rectal cancer: Initial results of the German CAO/ARO/AIO-04 randomised phase 3 trial. Lancet Oncol. 2012, 13, 679–687. [Google Scholar] [PubMed]
  22. Rödel, C.; Graeven, U.; Fietkau, R.; Hohenberger, W.; Hothorn, T.; Arnold, D.; Hofheinz, R.-D.; Ghadimi, M.; Wolff, H.A.; Lang-Welzenbach, M.; et al. German Rectal Cancer Study Group Oxaliplatin added to fluorouracil-based preoperative chemoradiotherapy and postoperative chemotherapy of locally advanced rectal cancer (the German CAO/ARO/AIO-04 study): Final results of the multicentre, open-label, randomised, phase 3 trial. Lancet Oncol. 2015, 16, 979–989. [Google Scholar] [PubMed]
  23. Hotchi, M.; Shimada, M.; Kurita, N.; Iwata, T.; Sato, H.; Morimoto, S.; Yoshikawa, K.; Higashijima, J.; Miyatani, T. microRNA expression is able to predict response to chemoradiotherapy in rectal cancer. Mol. Clin. Oncol. 2013, 1, 137–142. [Google Scholar] [CrossRef]
  24. Toyota, M.; Suzuki, H.; Sasaki, Y.; Maruyama, R.; Imai, K.; Shinomura, Y.; Tokino, T. Epigenetic silencing of microRNA-34b/c and B-cell translocation gene 4 is associated with CpG island methylation in colorectal cancer. Cancer Res. 2008, 68, 4123–4132. [Google Scholar] [CrossRef] [PubMed]
  25. Quah, H.-M.; Chou, J.F.; Gonen, M.; Shia, J.; Schrag, D.; Saltz, L.B.; Goodman, K.A.; Minsky, B.D.; Wong, W.D.; Weiser, M.R. Pathologic stage is most prognostic of disease-free survival in locally advanced rectal cancer patients after preoperative chemoradiation. Cancer 2008, 113, 57–64. [Google Scholar] [CrossRef] [PubMed]
  26. Li, C.-F.; He, H.-L.; Wang, J.-Y.; Huang, H.-Y.; Wu, T.-F.; Hsing, C.-H.; Lee, S.-W.; Lee, H.-H.; Fang, J.-L.; Huang, W.-T.; et al. Fibroblast growth factor receptor 2 overexpression is predictive of poor prognosis in rectal cancer patients receiving neoadjuvant chemoradiotherapy. J. Clin. Pathol. 2014, 67, 1056–1061. [Google Scholar] [CrossRef] [PubMed]
  27. García-Flórez, L.J.; Gómez-Álvarez, G.; Frunza, A.M.; Barneo-Serra, L.; Martínez-Alonso, C.; Fresno-Forcelledo, M.F. Predictive markers of response to neoadjuvant therapy in rectal cancer. J. Surg. Res. 2015, 194, 120–126. [Google Scholar] [CrossRef] [PubMed]
  28. He, H.-L.; Lee, Y.-E.; Chen, H.-P.; Hsing, C.-H.; Chang, I.-W.; Shiue, Y.-L.; Lee, S.-W.; Hsu, C.-T.; Lin, L.-C.; Wu, T.-F.; et al. Overexpression of DNAJC12 predicts poor response to neoadjuvant concurrent chemoradiotherapy in patients with rectal cancer. Exp. Mol. Pathol. 2015, 98, 338–345. [Google Scholar] [CrossRef] [PubMed]
  29. Tsang, J.S.; Vencken, S.; Sharaf, O.; Leen, E.; Kay, E.W.; McNamara, D.A.; Deasy, J.; Mulligan, E.D. Global DNA methylation is altered by neoadjuvant chemoradiotherapy in rectal cancer and may predict response to treatment—A pilot study. Eur. J. Surg. Oncol. 2014, 40, 1459–1466. [Google Scholar] [CrossRef] [PubMed]
  30. Kim, I.Y.; You, S.H.; Kim, Y.W. Neutrophil-lymphocyte ratio predicts pathologic tumor response and survival after preoperative chemoradiation for rectal cancer. BMC Surg. 2014, 14. [Google Scholar] [CrossRef] [PubMed]
  31. Memon, S.; Lynch, A.C.; Akhurst, T.; Ngan, S.Y.; Warrier, S.K.; Michael, M.; Heriot, A.G. Systematic review of FDG-PET prediction of complete pathological response and survival in rectal cancer. Ann. Surg. Oncol. 2014, 21, 3598–3607. [Google Scholar] [CrossRef] [PubMed]
  32. Yiu, R.; Wong, S.K.; Cromwell, J.; Madoff, R.D.; Rothenberger, D.A.; Garcia-Aguilar, J. Pelvic wall involvement denotes a poor prognosis in T4 rectal cancer. Dis. Colon Rectum 2001, 44, 1676–1681. [Google Scholar] [CrossRef] [PubMed]
  33. Hwang, K.; Park, I.J.; Yu, C.S.; Lim, S.-B.; Lee, J.L.; Yoon, Y.S.; Kim, C.W.; Kim, J.C. Impression of prognosis regarding pathologic stage after preoperative chemoradiotherapy in rectal cancer. World J. Gastroenterol. 2015, 21, 563–570. [Google Scholar] [CrossRef] [PubMed]
  34. Slattery, M.L.; Herrick, J.S.; Bondurant, K.L.; Wolff, R.K. Toll-like receptor genes and their association with colon and rectal cancer development and prognosis. Int. J. Cancer 2012, 130, 2974–2980. [Google Scholar] [CrossRef] [PubMed]
  35. Huh, J.W.; Lee, W.Y.; Kim, S.H.; Park, Y.A.; Cho, Y.B.; Yun, S.H.; Kim, H.C.; Park, H.C.; Choi, D.H.; Park, J.O.; et al. Immunohistochemical detection of p53 expression in patients with preoperative chemoradiation for rectal cancer: Association with prognosis. Yonsei Med. J. 2015, 56, 82–88. [Google Scholar] [CrossRef] [PubMed]
  36. Gaedcke, J.; Leha, A.; Claus, R.; Weichenhan, D.; Jung, K.; Kitz, J.; Grade, M.; Wolff, H.A.; Jo, P.; Doyen, J.; et al. Identification of a DNA methylation signature to predict disease-free survival in locally advanced rectal cancer. Oncotarget 2014, 5, 8123–8135. [Google Scholar] [CrossRef] [PubMed]
  37. Jo, P.; Jung, K.; Grade, M.; Conradi, L.-C.; Wolff, H.A.; Kitz, J.; Becker, H.; Rüschoff, J.; Hartmann, A.; Beissbarth, T.; et al. CpG island methylator phenotype infers a poor disease-free survival in locally advanced rectal cancer. Surgery 2012, 151, 564–570. [Google Scholar] [CrossRef] [PubMed]
  38. Ghadimi, B.M.; Grade, M.; Difilippantonio, M.J.; Varma, S.; Simon, R.; Montagna, C.; Füzesi, L.; Langer, C.; Becker, H.; Liersch, T.; et al. Effectiveness of gene expression profiling for response prediction of rectal adenocarcinomas to preoperative chemoradiotherapy. J. Clin. Oncol. 2005, 23, 1826–1838. [Google Scholar] [CrossRef] [PubMed]
  39. Scapoli, L.; Palmieri, A.; Lo Muzio, L.; Pezzetti, F.; Rubini, C.; Girardi, A.; Farinella, F.; Mazzotta, M.; Carinci, F. MicroRNA expression profiling of oral carcinoma identifies new markers of tumor progression. Int. J. Immunopathol. Pharmacol. 2010, 23, 1229–1234. [Google Scholar] [PubMed]
  40. Mancikova, V.; Castelblanco, E.; Pineiro-Yanez, E.; Perales-Paton, J.; de Cubas, A.A.; Inglada-Perez, L.; Matias-Guiu, X.; Capel, I.; Bella, M.; Lerma, E.; et al. MicroRNA deep-sequencing reveals master regulators of follicular and papillary thyroid tumors. Mod. Pathol. 2015, 28, 748–757. [Google Scholar] [CrossRef] [PubMed]
  41. Pardo, O.E.; Castellano, L.; Munro, C.E.; Hu, Y.; Mauri, F.; Krell, J.; Lara, R.; Pinho, F.G.; Choudhury, T.; Frampton, A.E.; et al. miR-515–5p controls cancer cell migration through MARK4 regulation. EMBO Rep. 2016, 17, 570–584. [Google Scholar] [CrossRef] [PubMed]
  42. Gilam, A.; Edry, L.; Mamluk-Morag, E.; Bar-Ilan, D.; Avivi, C.; Golan, D.; Laitman, Y.; Barshack, I.; Friedman, E.; Shomron, N. Involvement of IGF-1R regulation by miR-515–5p modifies breast cancer risk among BRCA1 carriers. Breast Cancer Res. Treat. 2013, 138, 753–760. [Google Scholar] [CrossRef] [PubMed]
  43. Lu, Z.; Luo, T.; Nie, M.; Pang, T.; Zhang, X.; Shen, X.; Ma, L.; Bi, J.; Wei, G.; Fang, G.; et al. TSPAN1 functions as an oncogene in gastric cancer and is downregulated by miR-573. FEBS Lett. 2015, 589, 1988–1994. [Google Scholar] [CrossRef] [PubMed]
  44. Danza, K.; de Summa, S.; Pinto, R.; Pilato, B.; Palumbo, O.; Merla, G.; Simone, G.; Tommasi, S. MiR-578 and miR-573 as potential players in BRCA-related breast cancer angiogenesis. Oncotarget 2015, 6, 471–483. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, H.-F.; Chen, H.; Ma, M.-W.; Wang, J.-A.; Tang, T.-T.; Ni, L.-S.; Yu, J.-L.; Li, Y.-Z.; Bai, B.-X. miR-573 regulates melanoma progression by targeting the melanoma cell adhesion molecule. Oncol. Rep. 2013, 30, 520–526. [Google Scholar] [PubMed]
  46. Quinn, E.M.; Wang, J.; Redmond, H.P. The emerging role of microRNA in regulation of endotoxin tolerance. J. Leukoc. Biol. 2012, 91, 721–727. [Google Scholar] [CrossRef] [PubMed]
  47. El Gazzar, M.; McCall, C.E. MicroRNAs distinguish translational from transcriptional silencing during endotoxin tolerance. J. Biol. Chem. 2010, 285, 20940–20951. [Google Scholar] [CrossRef] [PubMed]
  48. Kuhn, D.E.; Nuovo, G.J.; Martin, M.M.; Malana, G.E.; Pleister, A.P.; Jiang, J.; Schmittgen, T.D.; Terry, A.V.; Gardiner, K.; Head, E.; et al. Human chromosome 21-derived miRNAs are overexpressed in down syndrome brains and hearts. Biochem. Biophys. Res. Commun. 2008, 370, 473–477. [Google Scholar] [CrossRef] [PubMed]
  49. Kuhn, D.E.; Nuovo, G.J.; Terry, A.V.; Martin, M.M.; Malana, G.E.; Sansom, S.E.; Pleister, A.P.; Beck, W.D.; Head, E.; Feldman, D.S.; et al. Chromosome 21-derived microRNAs provide an etiological basis for aberrant protein expression in human Down syndrome brains. J. Biol. Chem. 2010, 285, 1529–1543. [Google Scholar] [CrossRef] [PubMed]
  50. Bofill-De Ros, X.; Santos, M.; Vila-Casadesús, M.; Villanueva, E.; Andreu, N.; Dierssen, M.; Fillat, C. Genome-wide miR-155 and miR-802 target gene identification in the hippocampus of Ts65Dn Down syndrome mouse model by miRNA sponges. BMC Genom. 2015, 16. [Google Scholar] [CrossRef]
  51. Kornfeld, J.W.; Baitzel, C.; Könner, A.C.; Nicholls, H.T.; Vogt, M.C.; Herrmanns, K.; Scheja, L.; Haumaitre, C.; Wolf, A.M.; Knippschild, U.; et al. Obesity-induced overexpression of miR-802 impairs glucose metabolism through silencing of Hnf1b. Nature 2013, 494, 111–115. [Google Scholar] [CrossRef] [PubMed]
  52. Yuan, F.; Wang, W. MicroRNA-802 suppresses breast cancer proliferation through downregulation of FoxM1. Mol. Med. Rep. 2015, 12, 4647–4651. [Google Scholar] [CrossRef] [PubMed]
  53. Tao, J.; Ji, J.; Li, X.; Ding, N.; Wu, H.; Liu, Y.; Wang, X.W.; Calvisi, D.F.; Song, G.; Chen, X. Distinct anti-oncogenic effect of various microRNAs in different mouse models of liver cancer. Oncotarget 2015, 6, 6977–6988. [Google Scholar] [CrossRef] [PubMed]
  54. Cao, Z.-Q.; Shen, Z.; Huang, W.-Y. MicroRNA-802 promotes osteosarcoma cell proliferation by targeting p27. Asian Pac. J. Cancer Prev. 2013, 14, 7081–7084. [Google Scholar] [CrossRef] [PubMed]
  55. Wang, L.Q.; Chen, G.; Liu, X.Y.; Liu, F.Y.; Jiang, S.Y.; Wang, Z. MicroRNA‑802 promotes lung carcinoma proliferation by targeting the tumor suppressor menin. Mol. Med. Rep. 2014, 10, 1537–1542. [Google Scholar] [CrossRef] [PubMed]
  56. Xiang, K.-M.; Li, X.-R. MiR-133b acts as a tumor suppressor and negatively regulates TBPL1 in colorectal cancer cells. Asian Pac. J. Cancer Prev. 2014, 15, 3767–3772. [Google Scholar] [CrossRef] [PubMed]
  57. Kano, M.; Seki, N.; Kikkawa, N.; Fujimura, L.; Hoshino, I.; Akutsu, Y.; Chiyomaru, T.; Enokida, H.; Nakagawa, M.; Matsubara, H. miR-145, miR-133a and miR-133b: Tumor-suppressive miRNAs target FSCN1 in esophageal squamous cell carcinoma. Int. J. Cancer 2010, 127, 2804–2814. [Google Scholar] [CrossRef] [PubMed]
  58. Bhattacharjya, S.; Roy, K.S.; Ganguly, A.; Sarkar, S.; Panda, C.K.; Bhattacharyya, D.; Bhattacharyya, N.P.; Roychoudhury, S. Inhibition of nucleoporin member Nup214 expression by miR-133b perturbs mitotic timing and leads to cell death. Mol. Cancer 2015, 14. [Google Scholar] [CrossRef] [PubMed]
  59. Qiu, T.; Zhou, X.; Wang, J.; Du, Y.; Xu, J.; Huang, Z.; Zhu, W.; Shu, Y.; Liu, P. miR-145, miR-133a and miR-133b inhibit proliferation, migration, invasion and cell cycle progression via targeting transcription factor Sp1 in gastric cancer. FEBS Lett. 2014, 588, 1168–1177. [Google Scholar] [CrossRef] [PubMed]
  60. Zhao, Y.; Huang, J.; Zhang, L.; Qu, Y.; Li, J.; Yu, B.; Yan, M.; Yu, Y.; Liu, B.; Zhu, Z. MiR-133b is frequently decreased in gastric cancer and its overexpression reduces the metastatic potential of gastric cancer cells. BMC Cancer 2014, 14. [Google Scholar] [CrossRef] [PubMed]
  61. Lin, C.W.; Li, X.R.; Zhang, Y.; Hu, G.; Guo, Y.H.; Zhou, J.Y.; Du, J.; Lv, L.; Gao, K.; Deng, H. TAp63 suppress metastasis via miR-133b in colon cancer cells. Br. J. Cancer 2014, 110, 2310–2320. [Google Scholar] [CrossRef] [PubMed]
  62. Boštjančič, E.; Jerše, M.; Glavač, D.; Zidar, N. miR-1, miR-133a/b, and miR-208a in human fetal hearts correlate to the apoptotic and proliferation markers. Exp. Biol. Med. 2015, 240, 211–219. [Google Scholar]
  63. Ferreira, L.R.P.; Frade, A.F.; Santos, R.H.B.; Teixeira, P.C.; Baron, M.A.; Navarro, I.C.; Benvenuti, L.A.; Fiorelli, A.I.; Bocchi, E.A.; Stolf, N.A.; et al. MicroRNAs miR-1, miR-133a, miR-133b, miR-208a and miR-208b are dysregulated in Chronic Chagas disease Cardiomyopathy. Int. J. Cardiol. 2014, 175, 409–417. [Google Scholar] [CrossRef] [PubMed]
  64. Koutsoulidou, A.; Mastroyiannopoulos, N.P.; Furling, D.; Uney, J.B.; Phylactou, L.A. Expression of miR-1, miR-133a, miR-133b and miR-206 increases during development of human skeletal muscle. BMC Dev. Biol. 2011, 11. [Google Scholar] [CrossRef] [PubMed]
  65. Boettger, T.; Wüst, S.; Nolte, H.; Braun, T. The miR-206/133b cluster is dispensable for development, survival and regeneration of skeletal muscle. Skelet. Muscle 2014, 4. [Google Scholar] [CrossRef] [PubMed]
  66. Chou, C.-K.; Yang, K.D.; Chou, F.-F.; Huang, C.-C.; Lan, Y.-W.; Lee, Y.-F.; Kang, H.-Y.; Liu, R.-T. Prognostic implications of miR-146b expression and its functional role in papillary thyroid carcinoma. J. Clin. Endocrinol. Metab. 2013, 98, E196–E205. [Google Scholar] [CrossRef] [PubMed]
  67. Lee, Y.S.; Lim, Y.S.; Lee, J.-C.; Wang, S.-G.; Park, H.-Y.; Kim, S.Y.; Lee, B.-J. Differential expression levels of plasma-derived miR-146b and miR-155 in papillary thyroid cancer. Oral Oncol. 2015, 51, 77–83. [Google Scholar] [CrossRef] [PubMed]
  68. Sun, M.; Fang, S.; Li, W.; Li, C.; Wang, L.; Wang, F.; Wang, Y. Associations of miR-146a and miR-146b expression and clinical characteristics in papillary thyroid carcinoma. Cancer Biomark. 2015, 15, 33–40. [Google Scholar] [PubMed]
  69. Deng, X.; Wu, B.; Xiao, K.; Kang, J.; Xie, J.; Zhang, X.; Fan, Y. MiR-146b-5p promotes metastasis and induces epithelial-mesenchymal transition in thyroid cancer by targeting ZNRF3. Cell. Physiol. Biochem. 2015, 35, 71–82. [Google Scholar] [CrossRef] [PubMed]
  70. Geraldo, M.V.; Yamashita, A.S.; Kimura, E.T. MicroRNA miR-146b-5p regulates signal transduction of TGF-β by repressing SMAD4 in thyroid cancer. Oncogene 2012, 31, 1910–1922. [Google Scholar] [CrossRef] [PubMed]
  71. Hardin, H.; Guo, Z.; Shan, W.; Montemayor-Garcia, C.; Asioli, S.; Yu, X.-M.; Harrison, A.D.; Chen, H.; Lloyd, R.V. The roles of the epithelial-mesenchymal transition marker PRRX1 and miR-146b-5p in papillary thyroid carcinoma progression. Am. J. Pathol. 2014, 184, 2342–2354. [Google Scholar] [CrossRef] [PubMed]
  72. Li, J.; Shan, F.; Xiong, G.; Wang, J.-M.; Wang, W.-L.; Xu, X.; Bai, Y. Transcriptional regulation of miR-146b by C/EBPβ LAP2 in esophageal cancer cells. Biochem. Biophys. Res. Commun. 2014, 446, 267–271. [Google Scholar] [CrossRef] [PubMed]
  73. Li, Y.; Wang, Y.; Yu, L.; Sun, C.; Cheng, D.; Yu, S.; Wang, Q.; Yan, Y.; Kang, C.; Jin, S.; et al. miR-146b-5p inhibits glioma migration and invasion by targeting MMP16. Cancer Lett. 2013, 339, 260–269. [Google Scholar] [CrossRef] [PubMed]
  74. Katakowski, M.; Buller, B.; Zheng, X.; Lu, Y.; Rogers, T.; Osobamiro, O.; Shu, W.; Jiang, F.; Chopp, M. Exosomes from marrow stromal cells expressing miR-146b inhibit glioma growth. Cancer Lett. 2013, 335, 201–204. [Google Scholar] [CrossRef] [PubMed]
  75. Drebber, U.; Lay, M.; Wedemeyer, I.; Vallböhmer, D.; Bollschweiler, E.; Brabender, J.; Mönig, S.P.; Hölscher, A.H.; Dienes, H.P.; Odenthal, M. Altered levels of the onco-microRNA 21 and the tumor-supressor microRNAs 143 and 145 in advanced rectal cancer indicate successful neoadjuvant chemoradiotherapy. Int. J. Oncol. 2011, 39, 409–415. [Google Scholar] [CrossRef] [PubMed]
  76. Kheirelseid, E.A.H.; Miller, N.; Chang, K.H.; Curran, C.; Hennessey, E.; Sheehan, M.; Newell, J.; Lemetre, C.; Balls, G.; Kerin, M.J. miRNA expressions in rectal cancer as predictors of response to neoadjuvant chemoradiation therapy. Int. J. Colorectal Dis. 2013, 28, 247–260. [Google Scholar] [CrossRef] [PubMed]
  77. Della Vittoria Scarpati, G.; Falcetta, F.; Carlomagno, C.; Ubezio, P.; Marchini, S.; de Stefano, A.; Singh, V.K.; D'Incalci, M.; de Placido, S.; Pepe, S. A specific miRNA signature correlates with complete pathological response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer. Int. J. Radiat. Oncol. Biol. Phys. 2012, 83, 1113–1119. [Google Scholar] [CrossRef] [PubMed]
  78. Mlcochova, J.; Faltejskova-Vychytilova, P.; Ferracin, M.; Zagatti, B.; Radova, L.; Svoboda, M.; Nemecek, R.; John, S.; Kiss, I.; Vyzula, R.; et al. MicroRNA expression profiling identifies miR-31–5p/3p as associated with time to progression in wild-type RAS metastatic colorectal cancer treated with cetuximab. Oncotarget 2015, 6, 38695–38704. [Google Scholar] [PubMed]
  79. Pichler, M.; Winter, E.; Ress, A.L.; Bauernhofer, T.; Gerger, A.; Kiesslich, T.; Lax, S.; Samonigg, H.; Hoefler, G. miR-181a is associated with poor clinical outcome in patients with colorectal cancer treated with EGFR inhibitor. J. Clin. Pathol. 2014, 67, 198–203. [Google Scholar] [CrossRef] [PubMed]
  80. Saridaki, Z.; Weidhaas, J.B.; Lenz, H.J.; Laurent-Puig, P.; Jacobs, B.; de Schutter, J.; de Roock, W.; Salzman, D.W.; Zhang, W.; Yang, D.; et al. A let-7 microRNA-binding site polymorphism in KRAS predicts improved outcome in patients with metastatic colorectal cancer treated with salvage cetuximab/panitumumab monotherapy. Clin. Cancer Res. 2014, 20, 4499–4510. [Google Scholar] [CrossRef] [PubMed]
  81. Sebio, A.; Paré, L.; Páez, D.; Salazar, J.; González, A.; Sala, N.; del Río, E.; Martín-Richard, M.; Tobeña, M.; Barnadas, A.; et al. The LCS6 polymorphism in the binding site of let-7 microRNA to the KRAS 3′-untranslated region: Its role in the efficacy of anti-EGFR-based therapy in metastatic colorectal cancer patients. Pharmacogenet. Genom. 2013, 23, 142–147. [Google Scholar] [CrossRef] [PubMed]
  82. Salendo, J.; Spitzner, M.; Kramer, F.; Zhang, X.; Jo, P.; Wolff, H.A.; Kitz, J.; Kaulfuß, S.; Beissbarth, T.; Dobbelstein, M.; et al. Identification of a microRNA expression signature for chemoradiosensitivity of colorectal cancer cells, involving miRNAs-320a, -224, -132 and let7g. Radiother. Oncol. 2013, 108, 451–457. [Google Scholar] [CrossRef] [PubMed]
  83. Bolstad, B.M.; Irizarry, R.A.; Astrand, M.; Speed, T.P. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003, 19, 185–193. [Google Scholar] [CrossRef] [PubMed]
  84. Smyth, G.K. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 2004, 3. [Google Scholar] [CrossRef] [PubMed]
  85. Cox, D.R.; Oakes, D. Analysis of Survival Data; CRC Press: London, UK, 1984. [Google Scholar]
  86. The R Project for Statistical Computing. Available online: http://www.r-project.org (accessed on 13 April 2016).
  87. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar]
Ijms 17 00568 g001 1024
Figure 1. Correlation of miR-515-5p, miR-573, miR-579, and miR-802 to overall survival an cancer-specific survival. Overall survival (left) and cancer-specific survival (right) of a set of 147 patients are illustrated as Kaplan-Meyer-curves for miR-515-5p, -573, -579, and -802. Red curve represents patients with a high expression of the mentioned miRNAs, blue curve those with a low expression. In all shown cases a high expression of the corresponding miRNA is significantly associated with a worse overall survival and cancer-specific survival.
Figure 1. Correlation of miR-515-5p, miR-573, miR-579, and miR-802 to overall survival an cancer-specific survival. Overall survival (left) and cancer-specific survival (right) of a set of 147 patients are illustrated as Kaplan-Meyer-curves for miR-515-5p, -573, -579, and -802. Red curve represents patients with a high expression of the mentioned miRNAs, blue curve those with a low expression. In all shown cases a high expression of the corresponding miRNA is significantly associated with a worse overall survival and cancer-specific survival.
Ijms 17 00568 g001
Ijms 17 00568 g002 1024
Figure 2. Correlation of miR-573 expression to tumor regression grade (TRG). A high tumor regression grade (TRG 4), which shows a good response to chemoradiotherapy (CRT), is significantly associated to a low expression of miR-573 (p < 0.05). Patients with a high expression of miR-573 show a low response to CRT (TRG 1).
Figure 2. Correlation of miR-573 expression to tumor regression grade (TRG). A high tumor regression grade (TRG 4), which shows a good response to chemoradiotherapy (CRT), is significantly associated to a low expression of miR-573 (p < 0.05). Patients with a high expression of miR-573 show a low response to CRT (TRG 1).
Ijms 17 00568 g002
Ijms 17 00568 g003 1024
Figure 3. miR-133b distant-metastasis-free survival plot.
Figure 3. miR-133b distant-metastasis-free survival plot.
Ijms 17 00568 g003
Ijms 17 00568 g004 1024
Figure 4. Study design: On the left, the initial set of 45 patients’ formalin-fixed paraffin-embedded (FFPE) samples analyzed via microarray analysis is shown. The results of that analysis showed a significant correlation (p < 0.05) of 14 miRNAs to overall survival (OS). These miRNAs together with 5 further miRNAs from previous work were analyzed in 147 patients’ samples (preserved in RNAlater), which were analyzed via quantitative real-time polymerase chain reaction (qPCR). The results are explained in the result paragraph.
Figure 4. Study design: On the left, the initial set of 45 patients’ formalin-fixed paraffin-embedded (FFPE) samples analyzed via microarray analysis is shown. The results of that analysis showed a significant correlation (p < 0.05) of 14 miRNAs to overall survival (OS). These miRNAs together with 5 further miRNAs from previous work were analyzed in 147 patients’ samples (preserved in RNAlater), which were analyzed via quantitative real-time polymerase chain reaction (qPCR). The results are explained in the result paragraph.
Ijms 17 00568 g004
Table
Table 1. p-values of the 19 microRNAs concerning correlation to the clinical parameters.
Table 1. p-values of the 19 microRNAs concerning correlation to the clinical parameters.
MicroRNAOverall SurvivalDisease-Free SurvivalCancer-Specific SurvivalDistant-Metastasis-Free SurvivalypN 1TRG 2
miR-133b0.19810.17970.10540.032 * 0.88690.5895
miR-146b0.59120.40650.79990.42820.0465 * 0.9552
miR-1980.51480.83780.46010.77710.30380.8873
miR-2230.84910.88870.28340.98010.13270.8693
miR-224 *0.10970.82560.23890.6640.49270.736
miR-23c0.64340.84180.4450.89750.23710.8907
miR-31330.75810.60150.59510.45160.1460.2699
miR-320a0.20160.49670.30110.29850.20240.8778
miR-34b0.96690.79040.90920.55060.68190.4109
miR-39410.06370.56560.05180.90550.72890.2858
miR-42630.37520.85410.30790.37680.91540.5357
miR-450b-3p0.3210.93340.37940.82340.47220.0878
miR-4970.28830.87040.48520.69420.88240.4389
miR-515-5p0.0021 * 0.66660.0132 * 0.97540.42350.867
miR-518f *0.34820.83580.52930.51790.11160.0849
miR-5730.0364 * 0.86720.0223 * 0.46980.09370.0416 *
miR-5790.0092 * 0.6950.0186 * 0.8960.09840.1266
miR-6120.120.69070.07340.9790.59310.4797
miR-8020.0231 * 0.17450.0168 * 0.4110.1520.1402
1 ypN = pathological nodal stage post preoperative chemoradiotherapy; 2 TRG = tumor regression grade according to Dworak; * = p < 0.05. Italic: miRNAs with significant correlation to patients’ survival in microarray analysis, which were analyzed here in an independent set of 147 samples via quantitative real-time polymerase chain reaction (qPCR) for validation.

Share and Cite

MDPI and ACS Style

Azizian, A.; Epping, I.; Kramer, F.; Jo, P.; Bernhardt, M.; Kitz, J.; Salinas, G.; Wolff, H.A.; Grade, M.; Beißbarth, T.; et al. Prognostic Value of MicroRNAs in Preoperative Treated Rectal Cancer. Int. J. Mol. Sci. 2016, 17, 568. https://doi.org/10.3390/ijms17040568

AMA Style

Azizian A, Epping I, Kramer F, Jo P, Bernhardt M, Kitz J, Salinas G, Wolff HA, Grade M, Beißbarth T, et al. Prognostic Value of MicroRNAs in Preoperative Treated Rectal Cancer. International Journal of Molecular Sciences. 2016; 17(4):568. https://doi.org/10.3390/ijms17040568

Chicago/Turabian Style

Azizian, Azadeh, Ingo Epping, Frank Kramer, Peter Jo, Markus Bernhardt, Julia Kitz, Gabriela Salinas, Hendrik A. Wolff, Marian Grade, Tim Beißbarth, and et al. 2016. "Prognostic Value of MicroRNAs in Preoperative Treated Rectal Cancer" International Journal of Molecular Sciences 17, no. 4: 568. https://doi.org/10.3390/ijms17040568

APA Style

Azizian, A., Epping, I., Kramer, F., Jo, P., Bernhardt, M., Kitz, J., Salinas, G., Wolff, H. A., Grade, M., Beißbarth, T., Ghadimi, B. M., & Gaedcke, J. (2016). Prognostic Value of MicroRNAs in Preoperative Treated Rectal Cancer. International Journal of Molecular Sciences, 17(4), 568. https://doi.org/10.3390/ijms17040568

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop