Journal of Cancer Science and Research

Journal of Cancer Science and Research
Open Access

ISSN: 2576-1447

+44 1478 350008

Research Article - (2018) Volume 0, Issue 0

Differential Expression of Mirnas in Colorectal Cancer.

Nihat Dilsiz1*, Ahmet Balik2, Filiz Mutaf3, Cemile Yesil4 and Ersin Borazan2
1Department of Molecular Biology and Genetics, University of Harran, Sanliurfa, Turkey
2School of Medicine, Gaziantep University, Gaziantep, Turkey
3Ay-Ka Ltd. Şti, Ankara, Turkey
4Synevo Medical Lab., Pathology Department, Turkey
*Corresponding Author: Nihat Dilsiz, Faculty of Arts and Sciences, Department of Molecular Biology and Genetics, University of Harran, Sanlıurfa, Turkey Email:

Keywords: Colon cancer; miRNAs; Microarray; Transcriptomics

Introduction

Colorectal cancer (CRC) is one of the deadliest diseases in the world and it is the second most common cause of cancer-related deaths in women, and the third in men. Although CRC mortality rates have declined in Western population, incidence rate of CRC has been increasing in Asia. Identified risk factors for CRC include: aging, gender, family history, lifestyle, pre-existing conditions such as, lynch syndrome, inflammatory bowel disease. CRC is most of the time asymptomatic until the late stages and can be only diagnosed when the cancer is already metastasized [1].

CRC was used to be considered as a disease caused by genetic mutations. However, it is currently viewed as a complex malignancy including epigenetic abnormalities in addition to genetic mutations. Three different pathogenic pathways have been identified which play role in CRC development: chromosomal instability, microsatellite instability and CpG island methylator phenotype. Chromosomal instability (CI) means deregulation of oncogenes such as, KRAS (functions as a transmitter of key extracellular signals) and TP53 (p53 functions as a key transcriptional regulator in cell cycle regulation) and also tumor suppressor genes. Microsatellite instability (MI) is caused by the inactivation of mismatch repair genes such as, hMLH1, hMSH2, hMSH6, and hPMS2. MI can result in Lynch syndrome and sporadic tumors. The CpG island methylator phenotype is caused by hypermethylation of CpG islands at tumor suppressors gene promotors [2,3].

The most common screening methods for CRC are colonoscopy, flexible sigmoidoscopy, guaiac-based fecal occult blood tests, and fecal immunochemical tests. Having the fact that these methods are still limited by unsatisfactory sensitivity and specificity; researchers have investigated the potential use of miRNAs as diagnostic and prognostic biomarkers for CRC screening [1].

During the past decade, it has been discovered that miRNAs play important roles in cancer initiation, progression and metastasis. Functionally, miRNAs can mediate cell proliferation, cell cycle progression, differentiation, metabolism, apoptosis, invasion to trigger tumor formation and angiogenesis [4].

MiRNA is a ≈ 20 nucleotides long non-protein-coding singlestranded RNA molecule. They are involved in almost all biological pathways, including gene expression, cell cycle regulation, cellular development, proliferation, differentiation and apoptosis [5-9].

MiRNAs are transcribed by RNA polymerase II in nucleus into long primary transcripts of miRNA as hairpins (pri-miRNAs, 1-3 kilobases). Nuclear RNase III enzyme (also called Drosha) cleaves primiRNA into a precursor (pre-miRNA, around ≈ 70 nucleotides in length) hairpin-shaped dsRNA. Pre-miRNAs are then exported from nucleus into cytoplasm by a nuclear export factor, Exportin-5 (XPO5). Then, protein complex of RNase III, Dicer, trims pre-miRNA in the cytoplasm to produce mature double-stranded (miRNA/miRNA complex) miRNA duplexes (≈ 20 nucleotides in length). Next, these two strands of the duplexes are unwounded by helicase, and then combined with RNA-induced silencing complex (RISC) containing an argonaut protein. One strand of the duplex is removed and degraded while the other one (single-stranded miRNA, ssmiRNA) is actively interacts with target mRNA to regulate the expression [6,8-13]. The seed region (6-8 nucleotides) at 5’ ends of miRNAs mainly bind to the 3’-untranslated region (3’-UTR) of their target mRNAs with imperfect sequence complementarity. However, recent studies have reported that miRNAs can also bind to 5’UTR, or open reading frame (ORF) of the target mRNA [8]. By binding to their target mRNA, miRNAs can inhibit expression of proteins from mRNA (suppression of translation) or causing mRNA degradation [12,14-18].

Since miRNAs can function as oncogenic (oncomir) as well as tumor suppressor depending on the function of their target genes, miRNAs are known to play an important role in tumorigenesis. A particular miRNA may be found up regulated in some cancer types as oncogenic function, but down regulated in other cancers, indicative of tumor suppressor function [19]. It is possible that the function of a single miRNA is cell or tissue specific. Recent studies indicate that miRNAs may have more oncogenic than tumor suppressive function in CRC [11,20,21]. 1.900 precursor and 2.600 mature miRNA have been identified up to now and these miRNAs are able to regulate the expression of almost 60% of all protein-coding genes in human genome [6,8,9,14,22-24]. Due to its small size, a single miRNA may act on expression of thousands of mRNAs. However, a single gene can be regulated only by multiple different miRNAs [11,13,23,25]. Therefore, miRNA expression profiles have been shown as potential signatures which are highly tissue specific for the classification, diagnosis and progression of cancer. Early detection of CRC by using miRNAs as diagnostic biomarkers provides the best chance for predictive diagnosis and successful treatment of diseases [1,21,26].

In 2003, it was first reported by Michael that differentially expressed upregulated and downregulated miRNAs were associated with CRC development and progression [27]. Since then more than 100 miRNAs have been shown to be associated with CRC since then. It was also shown that downregulation of let-7a-1 in CRC samples. MiR-34a has been shown to block cell proliferation and it is downregulated in colon cancer cells, whereas miR-31, miR-96, miR-31, miR-135b, and miR-183 were up regulated in CRC cells in the same study [28].

MiRNA associations with specific tumor markers also have been identified. Loss of miR-34 expression has been shown to deteriorate TP53-mediated cell death while overexpression of miR-34 results in apoptosis. MiRNAs also have been exhibited to be associated with regulation of oncogenes and with tumor suppressor genes. MiR-143 and let-7 expression levels have been associated with KRAS2 mutations [21]. Furthermore, miRNA-21 and miRNA-183 were found to be upregulated in CRC and miRNA-143, miRNA-145 and miRNA-497 are down-regulated when compared to healthy subjects as stated by different studies [3,6,7,14,17,21,29-32]. Compared with those in normal tissue, miRNAs expression levels are down-regulated in malignant tissue, only several miRNAs are up regulated; most of them play oncogenic roles [8,15,33]. On the other hand, tumor suppressor’s miRNA-497 and miRNA-378c are mainly downregulated in colorectal cancer [29,30].

The first miRNA mimic entered the clinic for cancer therapy is synthetic miRNA-34 (tumor suppressor) loaded in liposomal nanoparticles [34]. In general, there are two strategies to developing miRNA-based therapeutics either by direct inhibition of the interaction between oncogenic miRNA and mRNA or replacement of target tumor suppressor miRNA genes to restore a loss of function [15]. Oncogenic miRNAs could be therapeutically targeted by repression. A simple method to block oncogenic miRNAs is the use of anti-miRNA oligonucleotides complementary to the sequence of the targeted mature oncogenic miRNAs (antagomirs). These oligonucleotides disrupt miRISC complex and prevent degradation of an mRNA which can then be translated. Another approach to more specifically inhibit the miRNA function is to use of miRNA masks which are complementary to binding sites in 3’UTR of target mRNA. MiRNA replacement therapy aims to substitution of tumor suppressor miRNAs absent or expressed at lower levels by using oligonucleotide mimics containing the same sequence as the mature miRNA. Restoration of silenced tumor suppressor miRNAs might produce beneficial effects on cell migration, invasion and increasing sensitivity to therapeutic agents [6,15,21,35].

Materials and Methods

Ethics statement

This study was approved by the Medical Ethical Committee of the Harran University, Turkey and has been performed in accordance with the ethical principles of the 2008 revised Declaration of Helsinki. The samples were used with written informed consents from patients and the approval of the Turkish Academy of Medical Sciences. Tissue samples (~0.5 cm3) from CRC patients were collected at surgery after the pathologist had confirmed the histopathology at the Hospital of Gaziantep Medical University. Normal tissue adjacent to the cancer tissue was used as control for this study. Subjects were signed by participants in a written informed consent form for this study. None of this CRC patients received chemotherapy or radiotherapy before the surgical resection. Each sample was placed in a cryovial and covered with the RNALater (Applied Biosystems) that prevents fragmentation of the fragile mRNA. Due to small sizes of miRNA (20-22 nt.), miRNA levels are remarkably stable in tissue samples when compared to mRNA. Tissue samples were homogenized in liquid nitrogen by using precellys homogenizer (Bertin) and stored at -86°C for subsequent analysis.

Total RNA isolation and reverse transcription of mature miRNAs

Total RNA containing miRNA was purified from homogenate by using the microRNAsy mini isolation kit (Qiagen) according to manufacturer’s instructions. The purity and concentration of total RNA was assessed by using NanoDrop spectrophotometer p360 (Implen). Total RNA samples were aliquoted and stored at -86°C until used.

Microarray profiling of miRNAs

MiRNA expression was profiled by using GeneChip miRNA 3.0 Array (Affymetrix) for five controls and seven CRC with stage I-IV samples. About 1,000 ng RNA from each sample were labeled by biotin using the FlashTag Biotin HSR RNA labeling kit (Affymetrix). Then, the labeled samples were hybridized with the human probes on the microarray chips for 16 hours at 48°C with a rotation speed of 60 rpm. After hybridization, the probe arrays were washed and stained by using the Fluidics Station 450 according to AGCC Fluidics Control Software. Then the fluorescence on the array was scanned using the Affymetrix® GeneChip® Scanner 3000 with a high resolution 6 g patch. The probe cell intensity files (*CEL files) generated by Affymetrix GeneChip® Command Console® software were transferred into probe level summarization files (*CHP files) using a robust multi-array (RMA) detection algorithm workflow. The *CHP files were further analyzed using Transcriptome Analysis Console (TAC) software, version 3.0, to identify and visualize the differentially expressed genes.

We identified the level of expressed miRNA with statistical significance from a volcano plot filtering between the colorectal cancerous and normal miRNAs from the experiment. Expression data of all probe sets detected by microarray analysis have been deposited in ArrayExpress, E-MTAB-4573.

Results and Discussion

In microarray analysis, we identified 279 aberrantly expressed miRNA genes that either have increased (250 miRNAs), or reduced (29 miRNAs) expression levels as shown in Table 1. We found that MiR-1201, miR-181-a, miR-7, miR-188, miR-552 and miR-183 were highly expressed in colon carcinoma tissues compared with the normal tissues (Table 2a). Among the overexpressed miRNAs, miR-1201 is the highest expressed one in CRC tissues. MiR-181, miR-183, 34a, miR-92, miR-21, miR-431 and miR-487a were identified as on co mirRNA in various tumor tissues [36-39]. MiR-21 has been shown to be overexpressed in many cancer types [18,24,37,39,40]. On the other hand, five miRNAs (miR-451, miR-486, miR-550, miR-635 and miR-10b) had decreased expression significantly in CRC tissue samples compared with the normal samples (Table 2b). These down regulated miRNAs were found to be tumor suppressors in CRC tissues in previous studies [29,30,39,41,42]. MiR-486 was found to be the lowest expressed one CRC tissues compared with the normal tissues. We also found some new aberrant expression miRNAs in colorectal cancer tissues (Tables 1-2b). Finally, we could expect that miRNAs have important role as tumor suppressors or oncogenic factors in the network of carcinogenesis. The genes with similar expression patterns are grouped together and connected by a series of branches (clustering tree or dendrogram). The threshold we used to screen up-regulated or down-regulated miRNAs is fold change ≥2.0 and p value˂0.05. The results of hierarchical clustering show aberrantly expressed miRNAs among CRC tissues and normal CR tissues (Figure 1). Principal component analysis was also used to find out the similarities and differences between CRC samples and normal tissues (Figure 2). Total RNA was obtained from tissues of control and CRC and hybridized to Affymetrix GeneChip® miRNA 3.0 Arrays. After normalization, differential miRNA expression data was analyzed by unsupervised hierarchical clustering.

Number Systematic name Tumor Bi-weight Avg
Signal (log2)
Control Bi-weight Avg
Signal (log2)
FoldChange (linear) (Tumor vs. Control) ANOVA p-value
(Tumor vs. Control)
1 hsa-mir-1201 6.65 1.77 29.64 0.000641
2 hsa-miR-181a-star 6.03 1.56 22.21 0.000461
3 hsa-miR-7 6.17 1.77 21.16 0.020387
4 hsa-mir-188 7.05 2.75 19.63 0.001616
5 hsa-mir-552 7.67 3.49 18.11 0.014192
6 hsa-mir-183 8.17 4.02 17.67 0.021005
7 hsa-miR-941 6.36 2.31 16.64 0.000003
8 U71d 6.29 2.33 15.57 0.021340
9 ACA3-2 7.68 3.81 14.67 0.000048
10 hsa-mir-34c 4.74 1.15 12.06 0.000309
11 hsa-miR-7-1-star 5.91 2.34 11.9 0.000627
12 hsa-mir-148a 5.38 1.81 11.9 0.001089
13 hsa-mir-409 6.08 2.54 11.61 0.000224
14 hsa-mir-769 5.19 1.67 11.49 0.000261
15 hsa-mir-331 6.2 2.68 11.43 0.007081
16 hsa-mir-10a 5.75 2.26 11.27 0.017370
17 ACA9 5.55 2.09 11.03 0.000002
18 hsa-mir-181d 8.24 4.81 10.75 0.004900
19 U71d 6.43 3.02 10.58 0.015358
20 hsa-mir-30e 5.55 2.15 10.55 0.047163
21 hsa-mir-455 5.84 2.47 10.36 0.000016
22 hsa-mir-495 5.29 2.01 9.74 0.000285
23 hsa-mir-941-3 5.37 2.11 9.52 0.007112
24 hsa-mir-146b 5.13 1.89 9.44 0.003417
25 hsa-mir-181c 5.51 2.3 9.3 0.003289
26 hsa-mir-299 5.04 1.87 9.03 0.036420
27 ACA9 6.14 2.96 9.02 0.000125
28 hsa-mir-18a 7.04 3.87 9 0.004278
29 hsa-mir-493 5.64 2.52 8.73 0.001531
30 hsa-miR-92a-1-star 5.82 2.73 8.49 0.005380
31 hsa-mir-301a 5.25 2.19 8.34 0.000267
32 ACA6 6.26 3.2 8.32 0.001425
33 HBII-99 8.66 5.61 8.28 0.020411
34 hsa-mir-431 5.44 2.43 8.05 0.000265
35 U70 7.25 4.25 7.99 0.013214
36 hsa-mir-487a 5.34 2.35 7.91 0.000055
37 hsa-mir-505 6.69 3.72 7.85 0.003072
38 ACA34 7.03 4.1 7.65 0.004997
39 hsa-mir-34a 6.1 3.18 7.54 0.000132
40 U65 6.29 3.38 7.5 0.013446
41 hsa-mir-224 7.95 5.06 7.44 0.006594
42 hsa-mir-154 6.18 3.29 7.41 0.014806
43 hsa-mir-181c 7.5 4.63 7.29 0.001522
44 hsa-mir-629 5.61 2.75 7.28 0.033213
45 hsa-mir-542 5.09 2.26 7.1 0.005752
46 hsa-mir-183 6.08 3.29 6.95 0.000749
47 hsa-mir-454 5.22 2.42 6.93 0.000161
48 U71c 4.83 2.07 6.79 0.001399
49 ACA55 5.51 2.74 6.79 0.005687
50 ENSG00000206903 10.38 7.62 6.76 0.005521
51 hsa-miR-181a-2-star 7.28 4.53 6.73 0.019146
52 hsa-mir-654 5.61 2.86 6.72 0.003405
53 hsa-mir-337 6.98 4.25 6.66 0.007153
54 hsa-miR-29b 8.03 5.29 6.66 0.012884
55 hsa-mir-622 4.84 2.11 6.65 0.046174
56 hsa-mir-135b 4.53 1.8 6.64 0.027261
57 U71a 6.63 3.91 6.56 0.011939
58 ACA24 5.95 3.3 6.29 0.000034
59 hsa-mir-429 7.3 4.65 6.26 0.001175
60 hsa-mir-411 5.79 3.15 6.21 0.000330
61 hsa-mir-192 9.41 6.83 5.97 0.013712
62 ACA52 8.2 5.64 5.87 0.010732
63 ACA5 5.05 2.51 5.84 0.000047
64 ENSG00000201199 5.57 3.04 5.8 0.001431
65 U53 7.44 4.9 5.78 0.008812
66 U67 5.21 2.74 5.55 0.000761
67 U71c 4.61 2.14 5.55 0.018617
68 hsa-mir-374b 5.65 3.2 5.46 0.049707
69 ENSG00000201592 5.64 3.21 5.4 0.000223
70 U46 8.64 6.21 5.39 0.020504
71 hsa-let-7i 5.38 2.95 5.38 0.000545
72 hsa-mir-376c 6.83 4.41 5.37 0.002284
73 hsa-mir-485 4.04 1.66 5.21 0.011868
74 ENSG00000207130 10.59 8.22 5.15 0.008658
75 hsa-mir-936 3.92 1.56 5.12 0.000137
76 hsa-mir-377 4.04 1.68 5.12 0.012927
77 hsa-mir-489 5.71 3.35 5.11 0.021337
78 hsa-mir-502 4.25 1.92 5.01 0.003544
79 hsa-mir-941-4 5.76 3.45 4.94 0.034572
80 hsa-mir-34c 4.73 2.44 4.9 0.000645
81 hsa-miR-550-star 6.61 4.34 4.82 0.032494
82 hsa-let-7g 3.66 1.39 4.8 0.000264
83 hsa-mir-361 5.31 3.05 4.77 0.002347
84 hsa-mir-362 3.85 1.6 4.75 0.000054
85 SNORA11B 3.94 1.71 4.69 0.009034
86 ACA66 3.49 1.29 4.58 0.002080
87 hsa-mir-217 3.54 1.35 4.56 0.046813
88 U71b 4.59 2.4 4.55 0.007406
89 hsa-mir-127 4.88 2.7 4.53 0.022632
90 HBII-316 7.19 5.03 4.48 0.023347
91 hsa-mir-370 6.01 3.85 4.47 0.011705
92 U107 5.22 3.07 4.46 0.019918
93 hsa-mir-433 5.7 3.56 4.43 0.004833
94 U103 4.83 2.69 4.4 0.014816
95 ACA10 7.17 5.03 4.39 0.002936
96 hsa-mir-543 4.86 2.73 4.39 0.006302
97 hsa-mir-140 6.87 4.74 4.37 0.024446
98 U46 6.65 4.54 4.34 0.002078
99 hsa-mir-501 7.43 5.31 4.34 0.009283
100 hsa-mir-326 4.34 2.23 4.32 0.000050
101 HBI-100 4.78 2.68 4.28 0.000311
102 hsa-mir-1254 4.86 2.76 4.28 0.001672
103 ACA62 5.54 3.45 4.26 0.001151
104 hsa-mir-196b 6.9 4.86 4.13 0.020518
105 ACA33 6.89 4.84 4.12 0.001894
106 hsa-mir-941-2 5.68 3.64 4.12 0.007328
107 hsa-mir-17 8.91 6.87 4.09 0.000788
108 hsa-mir-147b 3.4 1.4 4.02 0.049155
109 U67 4.25 2.25 3.99 0.017309
110 ACA3-2 9.62 7.63 3.96 0.000053
111 U31 9.14 7.16 3.96 0.002425
112 ACA3 7.61 5.62 3.95 0.005823
113 ACA24 10.34 8.36 3.95 0.011863
114 ACA5 5.73 3.75 3.93 0.000789
115 hsa-mir-299 4.47 2.52 3.87 0.001478
116 HBII-99B 5.39 3.44 3.87 0.008777
117 hsa-mir-96 3.45 1.5 3.87 0.012393
118 hsa-mir-1291 3.91 1.97 3.85 0.002163
119 hsa-mir-425 8.46 6.52 3.82 0.008893
120 ACA26 4.8 2.87 3.8 0.002390
121 hsa-mir-3200 5.02 3.1 3.79 0.011422
122 hsa-miR-3130-5p 3.48 1.58 3.73 0.004953
123 hsa-mir-200a 10.98 9.09 3.7 0.004609
124 ACA31 4.53 2.64 3.7 0.009703
125 ACA17 3.57 1.68 3.69 0.000557
126 ACA16 7.48 5.62 3.64 0.001931
127 hsa-mir-199b 5.25 3.39 3.62 0.004455
128 hsa-mir-221 5.35 3.49 3.62 0.033553
129 hsa-mir-3138 3.82 1.97 3.61 0.035768
130 U58C 6.04 4.2 3.57 0.003545
131 ACA43 5.15 3.32 3.56 0.000473
132 U68 8.18 6.35 3.56 0.007408
133 hsa-mir-21 8.47 6.63 3.56 0.036590
134 U19 6.04 4.21 3.54 0.015532
135 hsa-mir-1201 4.95 3.13 3.53 0.004735
136 hsa-let-7d 3.5 1.71 3.45 0.049116
137 ENSG00000206913 9.34 7.58 3.4 0.000693
138 ACA23 3.48 1.73 3.36 0.000177
139 hsa-miR-103-2-star 3.69 1.94 3.35 0.000714
140 ACA46 3.92 2.19 3.3 0.004512
141 U15B 7.04 5.32 3.3 0.005781
142 ENSG00000206603 3.64 1.93 3.27 0.031852
143 ACA66 3.39 1.69 3.26 0.000431
144 ACA16 4.97 3.27 3.25 0.003536
145 hsa-miR-1285 5.76 4.06 3.25 0.008728
146 hsa-mir-660 8.76 7.07 3.22 0.000667
147 U18A 4.04 2.36 3.22 0.011729
148 hsa-mir-1292 4.7 3.02 3.22 0.020186
149 hsa-miR-24-1-star 3.14 1.46 3.2 0.002881
150 snR38B 5.43 3.77 3.15 0.005555
151 hsa-mir-1183 4.44 2.78 3.15 0.013904
152 HBII-180C 8.22 6.56 3.14 0.010323
153 HBI-6 6.33 4.7 3.1 0.000781
154 hsa-miR-103-as 2.93 1.29 3.1 0.003721
155 hsa-mir-424 3.39 1.76 3.1 0.013516
156 ACA41 9.86 8.24 3.09 0.043374
157 hsa-mir-4317 4.24 2.66 3 0.007898
158 U14B 5.12 3.55 2.98 0.000241
159 hsa-mir-18b 8.01 6.44 2.97 0.048861
160 hsa-mir-1259 4 2.43 2.96 0.005410
161 E2 5.41 3.84 2.95 0.011803
162 U23 3.44 1.91 2.9 0.042505
163 U71b 3.72 2.19 2.89 0.000953
164 14qI-4 3.61 2.09 2.88 0.003175
165 hsa-mir-2276 3.69 2.17 2.87 0.000172
166 U99 6.76 5.24 2.87 0.019247
167 hsa-mir-4257 2.88 1.37 2.86 0.008870
168 ACA34 7.2 5.7 2.83 0.038452
169 v49_ENSG00000206633 3.43 1.93 2.83 0.039065
170 hsa-miR-181b 11.76 10.27 2.8 0.019370
171 ACA15 4.6 3.11 2.8 0.030689
172 ENSG00000238956 7.95 6.47 2.79 0.006921
173 U68 10.87 9.39 2.79 0.009064
174 ENSG00000222489 5.95 4.48 2.77 0.005702
175 hsa-mir-500 9.15 7.69 2.75 0.013823
176 ENSG00000252213 4.59 3.13 2.75 0.013963
177 U48 9.58 8.12 2.75 0.026235
178 U64 4.88 3.43 2.74 0.000422
179 U96b 4.67 3.22 2.74 0.005474
180 ENSG00000200879 8.8 7.35 2.73 0.004711
181 U54 9.23 7.78 2.73 0.006385
182 U27 8.39 6.95 2.71 0.009984
183 U44 11.36 9.92 2.7 0.005082
184 U23 8.9 7.47 2.7 0.039786
185 hsa-mir-532 9.26 7.83 2.69 0.001247
186 ACA67B 3.23 1.8 2.69 0.015256
187 ACA7B 9.29 7.88 2.66 0.011805
188 U60 5.58 4.18 2.65 0.000133
189 U18C 3.27 1.86 2.65 0.001589
190 hsa-mir-148a 9.72 8.31 2.65 0.002850
191 U42B 3.7 2.3 2.64 0.000125
192 ENSG00000207187 7.21 5.81 2.64 0.002058
193 ACA32 5.95 4.55 2.64 0.007992
194 U103B 5.17 3.79 2.61 0.001744
195 U17a 5.01 3.63 2.6 0.001582
196 U106 4.59 3.22 2.6 0.012716
197 snR39B 8.08 6.7 2.6 0.015181
198 hsa-mir-93 3.41 2.04 2.58 0.001656
199 SNORA11B 3.63 2.27 2.56 0.006011
200 14qII-26 4.28 2.93 2.56 0.023833
201 U35A 9.92 8.56 2.56 0.033458
202 ACA2b 3.83 2.48 2.54 0.005727
203 mgU6-53 3.85 2.51 2.53 0.000996
204 snR38A 4.86 3.52 2.53 0.001270
205 ACA18 10.27 8.94 2.52 0.017484
206 SNORA84 4.26 2.93 2.52 0.021507
207 hsa-mir-22 5.96 4.62 2.52 0.034913
208 HBII-180A 9.58 8.25 2.5 0.013371
209 hsa-mir-185 4.1 2.78 2.5 0.023328
210 U46 7.63 6.32 2.49 0.011107
211 hsa-mir-141 9.85 8.55 2.47 0.033946
212 ACA41 4.48 3.18 2.47 0.043068
213 hsa-mir-132 2.56 1.26 2.46 0.013196
214 hsa-mir-493 2.47 1.17 2.46 0.030161
215 hsa-mir-654 3.16 1.87 2.45 0.002769
216 U92 4.26 2.97 2.45 0.003824
217 ACA47 3.02 1.73 2.45 0.013111
218 ACA32 5.71 4.42 2.44 0.016438
219 mgh28S-2411 11.22 9.95 2.42 0.002499
220 hsa-mir-744 2.99 1.72 2.41 0.017194
221 HBI-61 4.28 3.01 2.41 0.025907
222 U13 6.73 5.48 2.38 0.035739
223 ACA50 4.89 3.65 2.37 0.000397
224 hsa-miR-376a 3.39 2.15 2.37 0.029751
225 hsa-mir-612 2.9 1.66 2.36 0.000370
226 hsa-mir-671 3.19 1.95 2.36 0.004397
227 U72 2.43 1.21 2.34 0.007780
228 U18A 4.46 3.27 2.29 0.016004
229 ACA48 8.5 7.31 2.27 0.040439
230 U42B 5.91 4.74 2.24 0.006167
231 U30 11.43 10.27 2.22 0.009547
232 U75 11.24 10.1 2.2 0.003475
233 hsa-mir-421 8.25 7.11 2.2 0.034311
234 U56 10.08 8.95 2.18 0.011403
235 hsa-mir-492 2.79 1.68 2.16 0.019203
236 hsa-mir-182 2.6 1.5 2.15 0.034690
237 HBII-95 2.71 1.62 2.12 0.005868
238 U109 4.11 3.03 2.12 0.009534
239 ACA11 5.71 4.62 2.12 0.023544
240 U36A 6.28 5.2 2.11 0.001139
241 hsa-mir-106b 12.64 11.58 2.1 0.002558
242 hsa-mir-711 2.97 1.9 2.09 0.009272
243 U106 5.36 4.3 2.09 0.041728
244 U47 3.68 2.62 2.08 0.028603
245 ACA53 3.31 2.26 2.07 0.002782
246 hsa-mir-1291 3.42 2.37 2.07 0.008497
247 U58C 8.26 7.23 2.05 0.025618
248 hsa-mir-222 2.56 1.53 2.05 0.030812
249 hsa-mir-212 7.03 6 2.05 0.041135
250 hsa-mir-222 13.12 12.09 2.04 0.014307
251 ENSG00000238807 1.23 2.23 -2.01 0.005883
252 ENSG00000238804 1.08 2.1 -2.02 0.018261
253 ENSG00000238430 1.35 2.38 -2.04 0.013001
254 14qI-1 1.46 2.5 -2.06 0.010180
255 v49_ENSG00000201733 1.36 2.51 -2.22 0.022257
256 ENSG00000239140 1.22 2.38 -2.24 0.023096
257 ENSG00000201619 4.34 5.52 -2.26 0.021703
258 hsa-mir-744 8.21 9.39 -2.27 0.015659
259 hsa-mir-497 9.66 10.86 -2.3 0.021457
260 hsa-mir-4314 1.74 2.95 -2.31 0.012926
261 hsa-mir-3141 8.28 9.5 -2.33 0.029632
262 hsa-mir-631 1.24 2.49 -2.38 0.006627
263 ENSG00000252921 2.81 4.06 -2.38 0.010667
264 ENSG00000252921 1.69 2.98 -2.45 0.033509
265 ENSG00000221345 1.41 2.71 -2.46 0.036777
266 hsa-mir-548x 1 2.3 -2.47 0.002773
267 hsa-mir-1281 7.03 8.42 -2.6 0.025311
268 ENSG00000239080 1.72 3.1 -2.6 0.047346
269 hsa-mir-606 1.29 2.73 -2.72 0.041951
270 ENSG00000238549 1.7 3.15 -2.74 0.003969
271 hsa-miR-548a-3p 2.29 3.8 -2.84 0.037466
272 ENSG00000238430 1.26 2.78 -2.85 0.024357
273 hsa-mir-550-2 2.09 4.04 -3.87 0.000059
274 hsa-mir-10b 9.04 11.05 -4.04 0.002456
275 hsa-mir-635 1.28 3.48 -4.58 0.019899
276 hsa-mir-550-1 4.85 7.12 -4.82 0.001534
277 HBII-85-6 6.54 9.08 -5.83 0.045268
278 hsa-mir-486 8.5 11.05 -5.89 0.024529
279 hsa-mir-451 8.32 11.04 -6.62 0.036823

Table 1: Microarray Analysis: Differentially expressed miRNAs in CRC samples compared with samples of controls in the training and validation sets. Statistical analysis was performed using TAC software one way analysis of variance (ANOVA) was used. The statistical significance level was set at p-value<0.05. Greater than 2-fold changes were analyzed for up/down regulated miRNAs. A total of 279 miRNAs were found to be differentially expressed in the above table.

Number Systematic name Tumor Bi-weight Avg Signal (log2) Control Bi-weight Avg Signal (log2) Fold Change (linear)
1 hsa-mir-1201 6.65 1.77 29.64
2 hsa-miR-181a-star 6.03 1.56 22.21
3 hsa-miR-7 6.17 1.77 21.16
4 hsa-mir-188 7.05 2.75 19.63
5 hsa-mir-552 7.67 3.49 18.11
6 hsa-mir-183 8.17 4.02 17.67
7 hsa-miR-941 6.36 2.31 16.64
8 U71d 6.29 2.33 15.57
9 ACA3-2 7.68 3.81 14.67
10 hsa-mir-34c 4.74 1.15 12.06

Table 2a: The most important significant increase and decrease rates of expression in miRNA transcripts, some upregulated miRNA transcripts.

Number Systematic name Tumor Bi-weight Avg
Signal (log2)
Control Bi-weight Avg Signal (log2) Fold Change (linear)
1 hsa-mir-606 1.29 2.73 -2.72
2 ENSG00000238549 1.7 3.15 -2.74
3 hsa-miR-548a-3p 2.29 3.8 -2.84
4 ENSG00000238430 1.26 2.78 -2.85
5 hsa-mir-550-2 2.09 4.04 -3.87
6 hsa-mir-10b 9.04 11.05 -4.04
7 hsa-mir-635 1.28 3.48 -4.58
8 hsa-mir-550-1 4.85 7.12 -4.82
9 HBII-85-6 6.54 9.08 -5.83
10 hsa-mir-486 8.5 11.05 -5.89

Table 2b: The most important significant increase and decrease rates of expression in miRNA transcripts, some down-regulated miRNA transcripts.

cancer-science-and-research-hierarchical-clustering

Figure 1: Hierarchical Clustering of Gene Expression.

cancer-science-and-research-component-analyses

Figure 2: Principle Component Analyses (PCA). Each red dot represents a tumor sample and each blue dot represents control sample. According to the graph tumor and control samples were clustered at different localities.

Each row represents individual miRNA. In hierarchical clustering, genes with similar expression patterns are grouped together and connected by a series of branches (clustering tree or dendrogram). The red colors on the dendrogram indicate expression levels higher than the median. While green colors show expression levels lower than the median.

Conclusion

Our results show that most of miRNAs differentially expressed are up regulated miRNAs in CRC tissues compared to control tissues. Evidences have been accumulating that some miRNAs are overexpressed in CRC cells and might function as inhibitors of different tumor suppressor genes. There are three strategies for miRNA-based therapies are to inhibit oncogenic miRNAs or restoring tumor suppressor miRNAs. The first strategy is to direct inhibition of oncogenic microRNAs and this can be achieved by using singlestranded antisense oligonucleotides (approximately 20-22 nucleotides in long) that act through complementar base-pairing with target miRNAs. The second strategy is tumor suppressor miRNA replacement. This involves reintroducing synthetic miRNA or expression vectors that will produce the miRNA of interest to restore a loss of function. The third strategy is to develop drugs that decrease the levels of oncogenic miRNAs or increase the levels of tumor-suppressor miRNA for cancer prevention [21,43].

The most important challenge is the identification of definitive miRNA signatures for CRC by large and comprehensive profiling studies. This will allow the identification of certain diagnostic and prognostic biomarkers that can help physicians with patient evaluation and also prospective markers for the early therapy [8,12,15,44]. In summary, miRNAs are extremely important regulators of oncogenes and tumor suppressor genes that are responsible for pathologic processes associated with malignant progression. Although the roles of many different miRNA have been identified in various tumors, more study is still needed to fully understand the role of each miRNA and assessing their roles in personalized miRNA targeted cancer therapy [24]. In future cancer therapy, miRNAs may be important in deciding which drugs are selected for a patient and in determination of whether the patient has responded to the drug [23]. In conclusion, further efforts to get more sensitive, fast and effective methods are needed to address the role of miRNAs biomarkers in clinical diagnostic, prognostic and therapy. Finally, miRNAs will emerge as a powerful resource to advance the diagnosis and management of cancer including colorectal cancer.

Acknowledgment

We thank C. Gover for data management.

References

  1. Ren A, Dong Y, Tsoi H, Yu J (2015) Detection of miRNA as non-invasive biomarkers of colorectal cancer.Int J MolSci 16: 2810-2823.
  2. Pancione M, Remo A, Colantuoni V (2012) Genetic and epigenetic events generate multiple pathways in colorectal cancer progression.Patholog Res Intpp:
  3. Corté H, Manceau G, Blons H, Laurent-Puig P (2012) MicroRNA and colorectal cancer.Dig Liver Dis 44: 195-200.
  4. Lin S, Gregory RI (2015) MicroRNA biogenesis pathways in cancer. Nat Rev Cancer 15: 321-333.
  5. Liu M, Chen H (2010) The role of microRNAs in colorectal cancer.J Genet Genomics 37: 347-358.
  6. Bonfrate L, Altomare DF, Di Lena M, Travaglio E, Rotelli MT, et al. (2013) MicroRNA in colorectal cancer: new perspectives for diagnosis, prognosis and treatment.J Gastrointestin Liver Dis 22: 311-320.
  7. Shen J, Stass SA, Jiang F (2013) MicroRNAs as potential biomarkers in human solid tumors.Cancer Lett 329: 125-136.
  8. Di Leva G, Garofalo M, Croce CM (2014) MicroRNAs in cancer.Annu Rev Pathol 9: 287-314.
  9. Acunzo M, Romano G, Wernicke D, Croce CM (2015) MicroRNA and cancer-a brief overview.AdvBiolRegul 57: 1-9.
  10. Abbott AL (2011) Introduction to Gene Regulation by Small RNAs. MicroRNAs in Development and Cancer 1: 1.
  11. Yi B, Piazza GA, Su X, Xi Y (2013) MicroRNA and cancer chemoprevention. Cancer Prev Res (Phila) 6: 401-409.
  12. Xing Z, Li D, Yang L, Xi Y, Su X (2014) MicroRNAs and anticancer drugs. ActaBiochimBiophys Sin (Shanghai) 46: 233-239.
  13. Mazeh H, Mizrahi I, Ilyayev N, Halle D, Brücher B, et al. (2013) The Diagnostic and Prognostic Role of microRNA in Colorectal Cancer-a Comprehensive review. J Cancer 4: 281-295.
  14. Bouyssou JM, Manier S, Huynh D, Issa S, Roccaro AM, et al. (2014) Regulation of microRNAs in cancer metastasis.BiochimBiophysActa 1845: 255-265.
  15. Zhu W, Xu B (2014) MicroRNA-21 identified as predictor of cancer outcome: a meta-analysis. PLoS One 9: e103373.
  16. Chen J, Wang BC, Tang JH (2012) Clinical significance of microRNA-155 expression in human breast cancer.J SurgOncol 106: 260-266.
  17. Esquela-Kerscher A, Slack FJ (2006) Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer 6: 259-269.
  18. Schetter AJ, Okayama H, Harris CC (2012) The role of microRNAs in colorectal cancer. Cancer J 18: 244-252.
  19. Pérez-Rivas LG, Jerez JM, Carmona R, de Luque V, Vicioso L, et al. (2014) A microRNA signature associated with early recurrence in breast cancer.PLoS One 9: e91884.
  20. Hayes J, Peruzzi PP, Lawler S (2014) MicroRNAs in cancer: biomarkers, functions and therapy. Trends Mol Med 20: 460-469.
  21. Sethi S, Ali S, Sethi S, Sarkar FH (2014) MicroRNAs in personalized cancer therapy.Clin Genet 86: 68-73.
  22. Huat TJ, Khan AA, Abdullah JM, Idris FM, Jaafar H (2015) MicroRNA Expression Profile of Neural Progenitor-Like Cells Derived from Rat Bone Marrow Mesenchymal Stem Cells under the Influence of IGF-1, bFGF and EGF. Int J MolSci16: 9693-9718.
  23. dePlanell-Saguer M, Rodicio MC(2013) Detection methods for microRNAs in clinic practice. J clinbiochem46: 869-878.
  24. Michael MZ, Susan M O'Connor, Nicholas G, van Holst Pellekaan(2003) Reduced Accumulation of Specific MicroRNAs in Colorectal Neoplasia. Mol Cancer Res 1: 882-891.
  25. Yang L, Belaguli N, Berger DH (2009) MicroRNA and colorectal cancer.World J Surg 33: 638-646.
  26. Hamfjord J, Stangeland AM, Hughes T, Skrede ML, Tveit KM, et al. (2012) Differential expression of miRNAs in colorectal cancer: comparison of paired tumor tissue and adjacent normal mucosa using high-throughput sequencing. PloS one7: e34150.
  27. Guo ST, Jiang CC, Wang GP, Li YP, Wang CY, et al. (2013) MicroRNA-497 targets insulin-like growth factor 1 receptor and has a tumour suppressive role in human colorectal cancer. Oncogene 32: 1910-1920.
  28. Ye JJ, Cao J (2014) MicroRNAs in colorectal cancer as markers and targets: Recent advances. World J Gastroenterol 20: 4288-4299.
  29. Zhou T, Zhang GJ, Zhou H, Xiao HX, Li Y (2014) Overexpression of microRNA-183 in human colorectal cancer and its clinical significance. Eur J GastroenterolHepatol 26: 229-233.
  30. Ogata-Kawata H, Izumiya M, Kurioka D, Honma Y, Yamada Y, et al. (2014) Circulating exosomal microRNAs as biomarkers of colon cancer. PLoS One 9: e92921.
  31. Bouchie A (2013) First microRNA mimic enters clinic. Nat Biotechnol 31: 577.
  32. Garzon R, Marcucci G, Croce CM(2010) Targeting microRNAs in cancer: rationale, strategies and challenges. Nature reviews Drug discovery 9: 775-789.
  33. Tanaka T, Sugaya S, Kita K, Arai M, Kanda T, et al. (2012) Inhibition of cell viability by human IFN-β is mediated by microRNA-431.Int J Oncol 40: 1470-1476.
  34. Brunet Vega A, Pericay C, Moya I, Ferrer A, Dotor E, et al. (2013) MicroRNA expression profile in stage III colorectal cancer: circulating miR-18a and miR-29a as promising biomarkers.Oncol Rep 30: 320-326.
  35. Ke TW, Hsu HL, Wu YH, Chen WT, Cheng YW, et al. (2014) MicroRNA-224 suppresses colorectal cancer cell migration by targeting Cdc42. Dis Markers.
  36. Slattery M, Herrick JS, Pellatt DF, Stevens JR, Mullany LE, et al. (2016) MicroRNA profiles in colorectal carcinomas, adenomas, and normal colonic mucosa: variations in miRNA expression and disease progression. Carcinogenesis37: 245-261.
  37. Shibuya H, Iinuma H, Shimada R, Horiuchi A, Watanabe T (2010) Clinicopathological and prognostic value of microRNA-21 and microRNA-155 in colorectal cancer.Oncology 79: 313-320.
  38. Faltejskova P, Svoboda M, Srutova K, Mlcochova J, Besse A, et al. (2012) Identification and functional screening of microRNAs highly deregulated in colorectal cancer. J Cell Mol Med16: 2655-2666.
  39. Menigatti M, Staiano T, Manser CN, Bauerfeind P, Komljenovic A,et al. (2013) Epigenetic silencing of monoallelically methylated miRNA loci in precancerous colorectal lesions. Oncogenesis2: e56.
  40. Chi Y, Zhou D (2016) MicroRNAs in colorectal carcinoma--from pathogenesis to therapy.J ExpClin Cancer Res 35: 43.
  41. Andreoli SC, Gasparini NJ, de Carvalho GP, Garicochea B, Pogue RE, et al. (2014) Use of microRNAs in directing therapy and evaluating treatment response in colorectal cancer.Einstein (Sao Paulo) 12: 256-258.
Citation: Dilsiz N, Balik A, Mutaf F, Yesil C, Borazan E (2016) Differential Expression of Mirnas in Colorectal Cancer. J Can Sci Res 3: 01.

Copyright: © 2016 Dilsiz N, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Top