Abstract
BACKGROUND:
Colorectal cancer (CRC) is a digestive tract malignancy microRNAs (miRNAs) have attracted much attention as biomarkers in tumor studies.
OBJECTIVE:
This work focused on the predictive potential and mechanism of miR-4310 in CRC.
METHODS:
The miRNA expression profile sets were obtained from the Gene Expression Omnibus (GEO) database, and the appropriate miRNA was screened by GEO2R. The CRC tissues and control tissues of 88 patients with CRC were collected, and the expression of miR-4310 was detected by quantitative real-time PCR, and the efficacy of miR-4310 in diagnosing CRC was evaluated by the receiver operating characteristic curve (ROC). The effects of miR-4310 on the proliferation, migration, and invasion of CRC cells were explored by Cell Counting Kit-8 (CCK-8) and Transwell experiments. Predicting the potential binding sites of miR-4310 and Runt-related transcription factor 1 (RUNX1) by four predictive websites. The relationship between miR-4310 and RUNX1 was confirmed by a double luciferase reporter gene experiment.
RESULTS:
The bioinformatics analysis found that miR-4310 was differentially expressed in CRC tissues and this finding was certified by the expression of miR-4310 in CRC tissues of collected patients and cultured CRC cell lines. The expression of miR-4310 had a predictive possibility for CRC patients. MiR-4310/RUNX1 pathway had effects on CRC viability, migration, and invasion.
CONCLUSION:
MiR-4310 had the potential to be a biomarker for early screening of CRC. MiR-4310 and RUNX1 participated in the regulation of CRC cells.
Introduction
Colorectal cancer (CRC) is located in the colon or rectum, including colon cancer and rectal cancer, which is a common malignant tumor in the digestive system [1]. Patients with CRC have no typical symptoms in the early stage, and the symptoms will become more and more obvious with the gradual development of the tumor, which will spread from local to the whole body [2]. The specific cause of CRC is not clear, which may be caused by the comprehensive effects of personal eating habits, the external environment, and genetic factors [4]. The main treatment methods for CRC are surgical resection, endoscopic submucosal dissection, drug chemotherapy, radiotherapy, and immune-targeted therapy [6]. Despite the continuous development of treatment technology, the incidence and mortality of CRC are still high, and the treatment effect is not satisfactory. Therefore, it is necessary to analyze the potential pathogenesis to help the clinical diagnosis of CRC.
The mechanism of occurrence and development of CRC is complicated. Multiple public papers document that miRNAs not only participate in the proliferation, migration, invasion, and apoptosis of CRC, and also function in tumor angiogenesis, immune system, and tumor-matrix interaction. MiR-21-5p, miR-27b-3p, and miR-15a-3p are involved in the pathogenesis of CRC [7, 8, 9]. MiR-452 is involved in the CRC progression via binding the ERK/MAPK pathway, reflecting its role in the development of CRC [10]. All these publications indicate that several miRNAs were involved in the development of CRC. MiR-4310 is a broadly researched miRNA in tumors, and its effects is revealed in several tumors MiR-4310 is a miRNA that took part in hepatocellular carcinoma, glioma, and lung cancer [11, 12, 13]. MiR-4310 is related to the regulation of lncRNA CASC15 in colon cancer [14]. These previous papers indicated that miR-4310 might be a marker in the development of some malignancies. However, information about miR-4310 on CRC is not elucidated.
In this work, we aimed to study the roles of miR-4310 in CRC. Therefore, the expression and diagnostic role of miR-4310 were evaluated. The impacts of miR-4310 on malignant phenotypes of CRC cells were analyzed.
Materials and methods
Screening of significantly different miRNA
In this study, the miRNA expression profile data sets including the CRC group and control group were screened from Gene Expression Omnibus (GEO) database. Finally, GSE156719 and GSE72281 were selected for the subsequent analysis. After providing and updating the GSE156719 dataset, it contains miRNA sequencing data of 3 tissue samples and 3 paired control tissues from patients. The GSE72281 dataset was submitted and updated, in which we chose miRNA sequencing results from 6 CRC tissues and 6 matched normal mucosa. Differential gene expression analysis was used to identify colorectal tumor samples by using the “GEO2R” interactive tool. The difference of miRNA expression levels between samples was selected by
Participants recruitment
Patients with CRC in People’s Hospital Affiliated to Fujian University of Traditional Chinese Medicine were included. All 88 patients were examined by colonoscopy, and polyethylene glycol electrolyte was used for intestinal cleaning. The inclusion criteria of the patient group were histopathological diagnosis of primary CRC. The exclusion criteria are recurrent CRC, other types of metastatic cancer, pregnant women and lactating women, accompanied by other gastrointestinal diseases, receiving radiotherapy and chemotherapy, and those with cognitive or communication obstacles. CRC tissue and normal tissues were obtained in the surgical resection of tumor tissues. All subjects signed written informed consent and agreed to use their samples. The research plan was reviewed and approved by the Ethics Committee of our hospital.
Cell foster and transfection
The human CRC cell lines (LOVO, HT29, SW620, Caco-2, and HCT116) and a normal control line (NCM460) are applied in this investigation. All these cell lines were purchased from A.T.C.C (Manassas, America). These cells were cultured in an incubator with 37 C, 95% minimum relative humidity, 5% carbon dioxide content and 95% oxygen content, and grew on the wall. The RPMI 1640 medium was changed every 2 to 3 days on average.
The centrifugal tubes filled with miR-4310-mimic or mimic control (con) dry powder, respectively were obtained from GenePharm (Shanghai, China). The expression vector pcDNA3.1 carrying RUNX1 or empty carriers was purchased from the same company. The Lipofectamine 3000 reagent and sequence to be transfected were mixed, respectively, and then the artificial sequences were transfected into cells with the help of Lipofectamine 3000 for 48 hours.
Verification of miRNA expression
The total RNA samples were gained from all tissue specimens using Trizol LS reagent (Thermo, San Jose, America). Hifair super mix (YEASEN, Shanghai, China) was prepared for 20
Verification of cell normal ability
Firstly, the CCK-8 assay was performed as following steps. A density of 2
Prediction and certification of target gene
The target prediction webs (miRDB, TargetScan, miRTarBase, and DIANA TOOLS) were used to download target genes. The Venn diagram of gene intersection was plotted by VENNY2.1. The sequence of the target gene was cloned in the pmirGLO carriers as the wide type of RUNX1. The mutant target site was cloned into pmirGLO luciferase reporter vectors. These obtained vectors and the artificial sequences of miR-4310 were transfected into cells. After 48 hours’ post-transfection, the intensity of luciferase in renilla and firefly was detected using the dual luciferase reporter assay system (Promega Corp., Madison, America).
Statistical analysis
The normal distribution was calculated by Shapiro-Wilk test. Data in accordance with the normal distribution were analyzed by paired Ttest. One-way and two-way ANOVA was used for comparing the discrepancy among three or more groups. The calculations were conducted by IBM SPSS, and the charts are plotted using GraphPad Prism. The clinical predictive possibility of miR-4310 was measured by plotting the receive operator characteristic (ROC) curve. The significant difference was represented by
Results
Differentially expressed miRNA in CRC
According to the screening criteria of
The differentially expressed miRNAs in CRC. (A) Differentially expressed miRNAs of GSE156719. (B) Differentially expressed miRNAs in GSE72281. (C) The Venn diagram was plotted using the data from GSE156719 and GSE72281.
The relative expression of miR-4310 and its potential to be a diagnostic biomarker. (A) The expression of miR-4310 was decreased in CRC tissues. (B) The high AUC value of the ROC curve showed that miR-4310 might act as a diagnostic indicator. ***P < 0.001, compared to control tissue.
The effect of miR-4310 on CRC cells. (A) The decreased miR-4310 was observed in CRC cell lines. (B) Transfection of mimics into LOVO and HT29 cells increased the expression of miR-4310. (C–D) Increased miR-4310 expression inhibited cell viability. (E) The migratory cells were suppressed by the overexpression of miR-4310. (F) Upregulation of miR-4310 inhibited the invasion of CRC cells. **P < 0.01, ***P < 0.001, compared to the mimic-con group.
Target gene of miR-4310. (A) The Venn diagram of predicted targets. (B) Bioinformatics analysis predicted the targeted sequence. (C) The reduction of luciferase activity proved the interaction of RUNX1 with miR-4310. (D) Increased miR-4310 restricted the expression of RUNX1. ***P < 0.001, compared to the mimic-con group.
RUNX1/miR-4310 participated in the bioactivity of CRC cells. (A) The expression of RUNX1 was reduced in miR-4310 mimic group and its expression was reversed by the transfection of p-RUNX1. (B–C) The inhibited cell viability was reversed by the upregulation of RUNX1. (D–E) RUNX1 affected the suppressed function of miR-4310 on cell migration and invasion. **P < 0.01, ***P < 0.001, compared to the mimic-con group; #P < 0.05, ###P < 0.001, compared to the miR-4310 mimic 
As expressed in Fig. 2A, we observed that circulating miR-4310 was remarkably lower in CRC tissues compared with control tissues.
Then, the value of miR-4310 in diagnosing CRC was evaluated by ROC analysis (Fig. 2B). MiR4310 had a high efficiency in screening CRC, with an AUC value of 0.908, a specificity of 0.830 and a sensitivity of 0.841. It can be seen from the above results that miR-4310 could effectively predict CRC, and could be used as new markers for CRC screening.
The impacts of miR-4310 on CRC cell lines
The decreased miR-4310 expression was observed in all CRC cells compared with NCM460 cells (Fig. 3A,
Target gene of miR-4310 in CRC
The intersection of target genes of miR-4310 in four websites was shown in Fig. 4A. The binding site between miR-4310 and RUNX1 was shown in Fig. 4B. The relative luciferase activity was suppressed in the group with miR-4310 mimic and wide-type RUNX1 in HT29 cells and LOVO cells (Fig. 4C,
RUNX1 participated in the antitumor effects of miR-4310
The transfection effect of this overexpression plasmid was shown in Fig. 5A, and in this picture, the mRNA expression of RUNX1 was inhibited by the increased expression of miR-4310 and improved by the cotransfection of miR-4310 mimics and p-RUNX1 (
Discussion
In the current study, we screened differently expressed miRNAs in patients with CRC. Then, the expression and diagnostic possibility of miR-4310 in patients with CRC were estimated. Finally, the potential effects of miR-4310 on CRC cells were researched. The expression of miR-4310 showed a reduction in CRC tissues and diagnostic value for patients with CRC. Overexpression of miR-4310 inhibited the viability, migration, and invasion of CRC cells by targeting RUNX1.
With the continuous innovation of clinical diagnosis and treatment technology, the mortality rate of CRC patients has effectively reduced, but there are still some defects in the existing CRC screening measures [15]. Scholars provide new insights for the screening of patients with CRC and the pathogenesis of CRC with the in-depth study of miRNA [16, 17]. In a clinical study of CRC, miR-99a, miR-1246, miR-1299, and miR-1539 are abnormally expressed in patients with CRC, and their possible predictive roles for distinguishing CRC patients are researched [18, 19, 20]. With the improvement of miRNA sequencing technology, scientists’ knowledge of malignant has been greatly promoted [21]. More and more public genome databases, such as the GEO, are commonly used to research potential tumor diagnostic markers [22]. In this study, firstly, we screened the abnormally expressed miRNAs in the tissues of CRC from the profiles GSE156719 and GSE72281. Then, we assessed the expression levels of circulating miR-4310 in tumor and control tissues from patients with CRC by qRT-PCR and analyzed their values as screening markers of CRC by ROC. Finally, the function and potential mechanism of miR-4310 in CRC cell lines were explored.
In this work, our team analyzed the data of GSE156719 and GSE72281 according to the data of CRC tissues and control tissues through NCBI online analysis software GEO2R and Venn. Finally, it was concluded that the expression of miR-4310, an obviously different miRNA, decreased in CRC tissues from patients. In different tumors, miR-4310 may play a role as a tumor suppressor or carcinogen. Previous studies showed that the levels of miR-4310 were declined in patients with hepatocellular carcinoma [12]. But in the tissues from glioma, Wu et al. report that the content of miR-4310 is overexpressed compared to the non-tumor brain tissues, underlining that miR-4310 might function as an activator in the carcinogenesis of glioma [13]. The main reasons causing this may be contributed to the different mechanisms of malignant progression in tissues. In this work, we observed that the expression level of miR-4310 in CRC tissue was obviously decreased, indicating miR-4310 acted as a tumor inhibitory role in CRC. Furthermore, miR-4310 could screen CRC patients with high sensitivity and specificity through the analysis of ROC.
Recent papers find that miRNAs are not only involved in tumorigenesis and development but also regulate the expression of oncogenes or suppressor indicators MiR-139-5p, miR-135b, and miR-93 are involved in the development and mechanism of CRC [23, 24, 25]. MiR-503-5p influences the phenotype of CRC cells by targeting SALL1, thus participating in the progression of CRC [26]. Yu et al. report that miR-17-5p has a suppressive effect on CRC cells’ normal abilities, and the rescue assays indicate that HAPB2 reversed the impacts of MiR-17-5p on CRC [27]. In this work, we verified that the expression of miR-4310 in CRC cell lines obviously declined, which was similar to the concentration of CRC tissues. To clarify the biological influence of miR-4310, we transfected miR-4310 mimics into LOVO cells and HT29 cells respectively and verified that the mimics could effectively increase the content of miR-4310 in the two cell lines, which indicated the feasibility of transfection experiment. CCK-8 experiment was used to explore the effect of miR-4310 on the activity of tumor cells. The results showed that miR-4310 mimics can effectively reduce the activity of tumor cells, reflecting that miR-4310 may influence the growth of CRC. Migration and invasion ability is an important biological phenotype of malignant tumors, so we used the Transwell experiment to explore the influence of miR-4310 on the migration and invasion ability of CRC cells. Overexpression of miR-4310 inhibited the migration and invasion ability of CRC cells. Hence, miR-4310 is a tumor suppressor of CRC, which can effectively inhibit the activity, migration and invasion of tumor cells, and may become a potential target of CRC drug therapy in the future.
Four predictive websites were performed to forecast the downstream targets of miR-4310, and through the intersection of all findings, we choose RUNX1 for the following study. According to the statement of an article of 2019, the expression of RUNX1 was elevated and it played roles in the exacerbation and metastasis process [28]. Lu et al. also underline that RUNX1 is upregulated in CRC tumor tissues compared to normal epithelial tissues and affects the migration of CRC cells, indicating that RUNX1 is associated with the mechanism of CRC [29]. The results of the current work suggested that overexpression of miR-4310 can significantly reduce the luciferase activity of wild-type RUNX1, but has no obvious effect on the luciferase activity of mutant RUNX1. In addition, we up-regulated the contents of miR-4310 and then detected the expression of RUNX1 gene. The results suggest that the overexpression of miR-4310 reduced the expression of RUNX1 in CRC cells. In in vitro function acquisition experiments, the miR-4310 could inhibit cell viability, migration, and invasion, while the increased RUNX1 expression reversed the function of miR-4310 on CRC cells, indicating that miR-4310/RUNX1 was correlated to the mechanism of CRC development. Limitations exist in this current study, Firstly, this work did not include patients with other intestines problems, for instance inflammatory bowel diseases or colorectal adenomas, and analyze the prognostic possibility of miR-4310. Secondly, the small sample size of included CRC patients might decrease the reliability of diagnostic analysis of miR-4310. Further study with a larger sample size needs perform to certified that whether miR-4310 expression can be used in normal clinical diagnosis.
Conclusion
The miR-4310 was screened through GEO database and we found that it was downexpressed in CRC tissues according to the GEO2R analysis. Based on our certification in patients with CRC, the expression of miR-4310 was decreased in the CRC tissues and it had the potential to be a biomarker for diagnosing CRC. MiR-4310 might function as a tumor suppressor in CRC by ameliorating the viability, migration, and invasion of CRC cells. RUNX1 is a downstream target participating in the regulation of miR-4310 in CRC.
Footnotes
Acknowledgments
The authors have no acknowledgments.
Conflict of interest
The authors declare that they have no conflict of interest.
Funding
The authors report no funding.
