Research Paper Volume 12, Issue 6 pp 5259—5279
Whole-transcriptome analysis reveals a potential hsa_circ_0001955/hsa_circ_0000977-mediated miRNA-mRNA regulatory sub-network in colorectal cancer
- 1 Department of Breast Surgery, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang Province, China
- 2 Program of Innovative Cancer Therapeutics, Division of Hepatobiliary and Pancreatic Surgery, Department of Surgery, First Affiliated Hospital, College of Medicine, Zhejiang University, Key Laboratory of Organ Transplantation, Hangzhou, Zhejiang Province, China
- 3 Key Laboratory of Organ Transplantation, Hangzhou, Zhejiang Province, China
- 4 Key Laboratory of Combined Multi-organ Transplantation, Ministry of Public Health, Hangzhou, Zhejiang Province, China
received: December 24, 2019 ; accepted: March 9, 2020 ; published: March 28, 2020 ;https://doi.org/10.18632/aging.102945
How to Cite
Copyright © 2020 Ding et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Background: Circular RNAs (circRNAs), a novel class of non-coding RNAs, have been found to act as microRNA (miRNA) sponges and thus play key roles in biological processes and pathogenesis. However, studies regarding circRNAs in colorectal cancer (CRC) remain inadequate.
Results: By differential expression analysis, 10 candidate circRNAs (6 upregulated and 4 downregulated circRNAs) were chosen. 9 of 10 circRNAs were available on CSCD and their structure showed the binding potential of miRNA. Intersection analysis revealed that miR-145-5p, miR-3127-5p, miR-761, miR-4766-3p, miR-135a-5p, miR-135b-5p, miR-374a-3p and miR-330-3p were 8 miRNAs with the most potential in binding circRNAs. Further expression validation and correlation analysis demonstrated hsa_circ_0001955/miR-145-5p and hsa_circ_0000977/miR-135b-5p axes as key pathways in CRC. Subsequently, target gene prediction, differential expression analysis, intersection analysis and correlation analysis showed that CDK6, MMP12 and RAB3IP were the three potential downstream targets of hsa_circ_0001955/miR-145-5p axis and FOXO1, MBNL1, MEF2C, RECK, PPM1E, TTLL7 and PCP4L1 were the seven potential downstream targets of hsa_circ_0000977/miR-135b-5p axis in CRC. Finally, we also confirmed that expression of hsa_circ_0001955 or hsa_circ_0000977 was significantly positively correlated with their individual targets in CRC.
Conclusions: In the present work, we constructed a potential hsa_circ_0001955/hsa_circ_0000977-mediated circRNA-miRNA-mRNA regulatory network in CRC by a series of in silico analysis and experimental validation.
Methods: Whole-transcriptome microarrays from CRC and matched normal samples were obtained from GEO. The structure of circRNA was identified by CSCD. starBase and miRNet were successively used to predict miRNA of circRNA and target gene of miRNA. Expression correlation between RNA-RNA interactions was assessed using GEO and TCGA data. Finally, a potential circRNA-miRNA-mRNA network was established based on competing endogenous RNA (ceRNA) hypothesis.
Colorectal cancer (CRC) ranks as the third most common malignancy and is the third leading cause of cancer-associated deaths in the United States . A lot of factors, including diet, polyps and chronic inflammation, have been documented to link to CRC pathogenesis, consisting of multiple processes. However, the alterations and roles of genomics and epigenetics in CRC initiation and progression are still rarely known. Moreover, despite received standard therapies, like surgery combined with chemotherapy and radiotherapy, many CRC patients will eventually recur or metastasize . Therefore, it is urgent to further investigate molecular mechanism of CRC carcinogenesis, and seek and develop effective therapeutic targets for patients with CRC.
Circular RNAs (circRNAs) are a class of novel endogenous non-coding RNAs (ncRNAs) with a covalently closed continuous loop, and they are generated by the selective cleavage of premature RNAs . The loop structure produced by linking 5’ and 3’ end renders circRNAs resistant to exonucleases and more stable than their corresponding linear transcripts . In 2011, Salmena L et al. proposed a competing endogenous RNA (ceRNA) hypothesis, in which ncRNAs could positively regulate mRNA expression through competitively binding to shared microRNAs (miRNAs) . Recent studies have well-demonstrated that circRNAs increase oncogene or tumor suppressor expression by ceRNA mechanism, and thereby participate in cancer development. For example, Xue D et al. found that circ-AKT3 inhibited clear cell renal cell carcinoma metastasis by modulating miR-296-3p/E-cadherin axis ; Han D et al. suggested that circMTO1 acted as the sponge of miR-9 to suppress progression of hepatocellular carcinoma ; Wang C et al. showed that circFOXO3 induced cell apoptosis in urothelial carcinoma by interaction with miR-191-5p . circRNAs have also been confirmed to play critical roles in various biological behaviors of CRC, including proliferation, migration and invasion [9, 10]. Moreover, circRNAs can serve as diagnostic or prognostic biomarkers for CRC patients [11, 12].
Establishment of circRNA-miRNA-mRNA regulatory networks may provide key clues for us to explore the detailed molecular mechanism of human disorders including malignancies. For example, Liu K et al. uncovered the possible mechanism of Hsp90 inhibitor-induced cell death of CRC by identification of a circRNA-miRNA-mRNA regulatory network . Yuan W et al. identified circRNAs as ceRNAs for miRNA-mRNA in CRC ; Song W et al. also constructed a circRNA-associated ceRNA network in CRC . The two studies constructed circRNA-miRNA-mRNA network using RNA data from different CRC patients.
In this study, we obtained and analyzed the circRNA-miRNA-mRNA ceRNA regulatory network using 10 CRC patients’ circRNA, miRNA and mRNA expression profile data. Furthermore, TCGA data were employed to validate the analytic results. A hsa_circ_0001955/hsa_circ_0000977-mediated miRNA-mRNA regulatory sub-network in colorectal cancer was successfully established. Finally, experimental validation was given. Those findings from the present work may provide novel insight into the molecular mechanism of CRC carcinogenesis.
Selection of 9 potential circRNAs in CRC
To find the potential circRNAs in CRC, circRNA dataset of 10 CRC patients (GSE126094) from GEO database was selected. GEO2R tool provided by GEO database was used to perform differential expression analysis for GSE126094. As listed in Figure 1A and Supplementary Table 1, a total of 1850 significant DECs, containing 1831 upregulated DECs and 19 downregulated DECs, were screened out. Next, we further selected some more potential DECs from the 1850 significant DECs based on the selection criterion of |log2FC| > 4, and 10 circRNAs, including 6 upregulated and 4 downregulated circRNAs, were finally identified (Figure 1B and Table 1). The circBase names and parental genes of the 10 circRNAs were presented in Table 2. Subsequently, the structure of 9 circRNAs, hsa_circ_0072088 (Figure 2A), hsa_circ_0000512 (Figure 2B), hsa_circ_0000511 (Figure 2C), hsa_circ_0001955 (Figure 2D), hsa_circ_0008274 (Figure 2E), hsa_circ_0001666 (Figure 2F), hsa_circ_0006220 (Figure 2G), hsa_circ_0000977 (Figure 2H) and hsa_circ_0043278 (Figure 2I) (hsa_circ_0000981 not available) was described in Figure 2 using the data from CSCD, indicating that all the 9 circRNAs had microRNA response elements (MREs). Taken together, the 9 circRNAs may be the key circRNAs in CRC by acting as miRNA sponges.
Figure 1. Identification of potential circRNAs in colorectal cancer. (A) The volcano plot of differentially expressed circRNAs (DECs) in colorectal cancer from GSE126094 dataset. The red dots and green dots represent upregulated DECs and downregulated DECs with significance (adjusted P-value < 0.05 and |log2FC| > 1), respectively. The black dots are those DECs without significance. (B) The heatmap of 10 potential DECs (|log2FC| > 4).
Table 1. The potential differentially expressed circRNAs (DECs) between colorectal cancer tissues and adjacent normal tissues.
Table 2. The parental genes of 10 potential circRNAs.
|circRNA ID||circBase ID||Parental gene|
Figure 2. Structural patterns of the 9 circRNAs in CSCD database. (A) The structural pattern of hsa_circ_0072088. (B) The structural pattern of hsa_circ_0000512. (C) The structural pattern of hsa_circ_0000511. (D) The structural pattern of hsa_circ_0001955. (E) The structural pattern of hsa_circ_0008274. (F) The structural pattern of hsa_circ_0001666. (G) The structural pattern of hsa_circ_0006220. (H) The structural pattern of hsa_circ_0000977. (I) The structural pattern of hsa_circ_0043278. (J) The representation of MRE, RBP and ORF.
Prediction and analysis of binding miRNAs of potential circRNAs in CRC
Next, miRNAs that might bind to the 9 circRNAs were predicted through starBase database. As listed in Supplementary Table 2, 14, 9, 9, 25, 6, 52, 8, 10 and 11 miRNAs were found to potentially bind to hsa_circ_0072088, hsa_circ_0000512, hsa_circ_0000511, hsa_circ_0001955, hsa_circ_0008274, hsa_circ_0001666, hsa_circ_0006220, hsa_circ_0000977 and hsa_circ_0043278, respectively. For better visualization, the circRNA-miRNA connecting network was established by Cytoscape software as presented in Figure 3. According to ceRNA hypothesis, circRNA expression should be theoretically inversely linked to miRNA expression. Thus, we firstly screened out DEmiRNAs between CRC tissues and adjacent normal tissues using miRNA dataset GSE126093 from 10 CRC patients. Identically, GEO2R tool was employed to perform differential expression analysis for GSE126093, and 178 significant DEmiRNAs, containing 148 upregulated DEmiRNAs and 30 downregulated DEmiRNAs, were identified. The 178 significant DEmiRNAs were listed in Supplementary Table 3 and vividly presented in Figure 4A. Subsequently, by intersection of predicted targeting miRNAs of circRNAs and significant DEmiRNAs, we discovered that miR-145-5p were commonly appeared in target miRNAs of upregulated DECs set and downregulated DEmiRNAs set (Figure 4B), and 7 miRNAs (miR-3127-5p, miR-761, miR-4766-3p, miR-135a-5p, miR-135b-5p, miR-374a-3p and miR-330-3p) were collectively appeared in target miRNAs of downregulated DECs set and upregulated DEmiRNAs set (Figure 4C). Followingly, expression levels of the 8 miRNAs were also validated using TCGA data as shown in Figure 4D–4K. The results demonstrated that only miR-145-5p (Figure 4D), miR-4766-3p (Figure 4G), miR-135b-5p (Figure 4I) and miR-374a-3p (Figure 4J) expression were in accordance with the previous analytic findings. Thus, miR-145-5p, miR-4766-3p, miR-135b-5p and miR-374a-3p may be the most potential downstream miRNAs of the key circRNAs in CRC. By search of the predicted circRNA-miRNA binding pairs in Figure 3 or Supplementary Table 2, miR-145-5p potentially bound to upregulated hsa_circ_0001955, and miR-4766-3p, miR-135b-5p and miR-374a-3p were three binding miRNAs of one common downregulated circRNA, hsa_circ_0000977 (Figure 5A). The expression correlation of the four circRNA-miRNA pairs was also determined using data from GSE126094 and GSE126093 (Figure 5B–5E). All the 4 miRNAs were negatively correlated with their corresponding paired circRNA, which fitted the ceRNA hypothesis. However, among the four pairs, only hsa_circ_0001955/miR-145-5p and hsa_circ_0000977/miR-135b-5p possessed statistically significant. Altogether, hsa_circ_0001955/miR-145-5p and hsa_circ_0000977/miR-135b-5p might be two important axes in carcinogenesis of CRC.
Figure 4. Identification of potential miRNAs that bind to circRNAs in colorectal cancer. (A) The volcano plot of differentially expressed miRNAs (DEmiRNAs) in colorectal cancer from GSE126093 dataset. The red dots and green dots represent upregulated DEmiRNAs and downregulated DEmiRNAs with significance (adjusted P-value < 0.05 and |log2FC| > 1), respectively. The black dots are those DEmiRNAs without significance. (B) The intersection analysis of target miRNAs of upregulated circRNAs and downregulated DEmiRNA. (C) The intersection analysis of target miRNAs of downregulated circRNAs and upregulated DEmiRNA. (D) The expression level of hsa-miR-145-5p in COAD. (E) The expression level of hsa-miR-3127-5p in COAD. (F) The expression level of hsa-miR-761 in COAD. (G) The expression level of hsa-miR-4766-3p in COAD. (H) The expression level of hsa-miR-135a-5p in COAD. (I) The expression level of hsa-miR-135b-5p in COAD. (J) The expression level of hsa-miR-374a-3p in COAD. (K) The expression level of hsa-miR-330-3p in COAD. “p<0.05” represents significant difference and “ns” represents no significance.
Figure 5. Correlation analysis for potential circRNA-miRNA pairs in colorectal cancer. (A) The circRNA-miRNA sub-network of interest. (B) The correlation of hsa_circ_0001955 with hsa-miR-145-5p in colorectal cancer. (C) The correlation of hsa_circ_0000977 with hsa-miR-4766-3p in colorectal cancer. (D) The correlation of hsa_circ_0000977 with hsa-miR-135b-5p in colorectal cancer. (E) The correlation of hsa_circ_0000977 with hsa-miR-374a-3p in colorectal cancer. P-value < 0.05 was considered as statistically significant.
Identification of downstream target genes of miRNAs in CRC
It has been widely acknowledged that miRNAs exert their biological functions by negatively regulated expression of downstream target genes. Therefore, the next step was to ascertain the downstream target genes of hsa_circ_0001955/miR-145-5p and hsa_circ_0000977/miR-135b-5p axes. Target genes of miR-145-5p and miR-135b-5p were predicted through miRNet database, and 238 and 83 target genes were forecasted to potentially bind to miR-145-5p and miR-135b-5p, respectively (Supplementary Table 4). DEGs between CRC and normal controls were obtained by GEO2R analysis using GSE126092 from the same 10 CRC patients with GSE126093 and GSE126094. As presented in Figure 6A and Supplementary Table 5, a total of 2787 DEGs, including 1194 upregulated DEGs and 1593 downregulated DEGs, were finally found. According to action mechanism of miRNA, miRNA expression should be negatively associated with expression of target gene. Thus, we performed intersection analysis for target genes set and DEGs set. As shown in Figure 6B, 19 genes (CBFB, CD44, CDK4, CDK6, CFTR, MMP1, MMP12, MYC, PIGF, POU5F1, ROBO2, SNTB1, TGFBI, VEGFA, PSAT1, LMNB2, RAB3IP, C11orf65 and CRNDE) were commonly appeared in miR-145-5p’s target genes set and upregulated DEGs set. And 10 genes (APC, FOXO1, MBNL1, MEF2C, RECK, KLF4, PPM1E, ARC, TTLL7 and PCP4L1) were commonly appeared in miR-135b-5p’s target genes set and downregulated DEGs set (Figure 6C). Next, expression levels of the 19 upregulated genes and 10 downregulated genes were further validated using TCGA data by starBase database. The results demonstrated that 17 of 19 upregulated genes (except CFTR and POU5F1) and 10 of 10 downregulated genes fitted the previous analytic findings as shown in Figure 6D and 6E, respectively. The 27 genes were considered as the downstream target genes of miR-145-5p and miR-135b-5p in CRC and selected for subsequent analysis.
Figure 6. Identification of potential targets of hsa-miR-145-5p and hsa-miR-135b-5p in colorectal cancer. (A) The volcano plot of differentially expressed genes (DEGs) in colorectal cancer from GSE126092 dataset. The red dots and green dots represent upregulated DEGs and downregulated DEGs with significance (adjust P-value < 0.05 and |log2FC| > 1), respectively. The black dots are those DEGs without significance. (B) The intersection analysis of target genes of downregulated hsa-miR-145-5p and upregulated DEGs. (C) The intersection analysis of target genes of upregulated hsa-miR-135b-5p and downregulated DEGs. (D) Validation of expression of 17 potential upregulated DEGs using starBase. (E) Validation of expression of 10 potential downregulated DEGs using starBase. “p<0.05” represents significant difference.
Establishment of a potential circRNA-miRNA-mRNA regulatory sub-network in CRC
Expression correlation of between miR-145-5p or miR-135b-5p and their corresponding target genes in CRC was assessed using the data from GSE126093 and GSE126092. The miR-145-5p-target gene network, containing 17 miRNA-target gene pairs, was constructed in Figure 7A. As presented in Figure 7B–7R, miR-145-5p expression was inversely linked to all the 17 target gene expression whereas only miR-145-5p/CBFB, miR-145-5p/CDK4, miR-145-5p/MMP1, miR-145-5p/MMP12, miR-145-5p/ROBO2, miR-145-5p/SNTB1, miR-145-5p/TGFBI, miR-145-5p/LMNB2, miR-145-5p/RAB3IP and miR-145-5p/CRNDE pairs had statistically significant. The miR-135b-5p-target gene network, involving 10 miRNA-target gene pairs, was shown in Figure 8A. Next, the expression correlation of miR-135b-5p with 10 target genes was also detected (Figure 8B–8K). The results suggested that 7 of 10 pairs (miR-135b-5p/FOXO1, miR-135b-5p/MBNL1, miR-135b-5p/MEF2C, miR-135b-5p/RECK, miR-135b-5p/PPM1E, miR-135b-5p/TTLL7 and miR-135b-5p/PCP4L1) possessed statistically significant. We further determine the expression relationship between miR-145-5p or miR-135b-5p and their target genes using TCGA data and the results were listed in Tables 3, 4, respectively. By intersection of analytic findings from GEO and TCGA database, we found that three genes (CDK6, MMP12 and RAB3IP) and seven genes (FOXO1, MBNL1, MEF2C, RECK, PPM1E, TTLL7 and PCP4L1) were the potential downstream targets of miR-145-5p and miR-135b-5p in CRC, respectively. According to ceRNA mechanism, circRNA expression should be positively associated with target gene expression. Therefore, we finally evaluated the expression correlation of hsa_circ_0001955 with CDK6/MMP12/RAB3IP or hsa_circ_0000977 with FOXO1/MBNL1/MEF2C/RECK/PPM1E/TTLL7/PCP4L1 in CRC (Figure 9A–9J). Intriguingly, all the 10 genes were significantly positively correlated with their corresponding circRNAs. Taken all these findings together, a potential circRNA-miRNA-mRNA regulatory sub-network in CRC was successively established, containing two pathways hsa_circ_0001955-miR-145-5p-CDK4/MMP12/RAB3IP and hsa_circ_0000977-miR-135b-5p-FOXO1/MBNL1/MEF2C/RECK/PPM1E/TTLL7/PCP4 L1 (Figure 10). To further affirm the analytic results, we collected 20 CRC cancer samples and adjacent normal samples from the First Affiliated Hospital of Zhejiang University. The expression levels of all RNAs in the constructed network in CRC were determined by qRT-PCR and expression relationships among them were also assessed. As presented in Figure 11A–11E and Figure 12A–12I, hsa_circ_0001955, hsa-miR-135b-5p, CDK4, MMP12 and RAB3IP were significantly upregulated whereas hsa_circ_0000977, hsa-miR-145-5p, FOXO1, MBNL1, MEF2C, RECK, PPM1E, TTLL7 and PCP4L1 were markedly downregulated in CRC when compared with normal controls, which were in accordance with GEO and TCGA data. Moreover, Figure 11F–11L and Figure 12J–12X showed a similar pattern of results compared to previous in silico analytic results, partially supporting the accuracy of in silico analysis in this study.
Figure 7. Correlation analysis for potential hsa-miR-145-5p-target gene pairs in colorectal cancer. (A) The visual hsa-miR-145-5p-target gene network. (B) The expression correlation of hsa-miR-145-5p with CBFB in colorectal cancer. (C) The expression correlation of hsa-miR-145-5p with CD44 in colorectal cancer. (D) The expression correlation of hsa-miR-145-5p with CDK4 in colorectal cancer. (E) The expression correlation of hsa-miR-145-5p with CDK6 in colorectal cancer. (F) The expression correlation of hsa-miR-145-5p with MMP1 in colorectal cancer. (G) The expression correlation of hsa-miR-145-5p with MMP12 in colorectal cancer. (H) The expression correlation of hsa-miR-145-5p with MYC in colorectal cancer. (I) The expression correlation of hsa-miR-145-5p with PIGF in colorectal cancer. (J) The expression correlation of hsa-miR-145-5p with ROBO2 in colorectal cancer. (K) The expression correlation of hsa-miR-145-5p with SNTB1 in colorectal cancer. (L) The expression correlation of hsa-miR-145-5p with TGFBI in colorectal cancer. (M) The expression correlation of hsa-miR-145-5p with VEGFA in colorectal cancer. (N) The expression correlation of hsa-miR-145-5p with PSAT1 in colorectal cancer. (O) The expression correlation of hsa-miR-145-5p with LMNB2 in colorectal cancer. (P) The expression correlation of hsa-miR-145-5p with RAB3IP in colorectal cancer. (Q) The expression correlation of hsa-miR-145-5p with C11orf65 in colorectal cancer. (R) The expression correlation of hsa-miR-145-5p with CRNDE in colorectal cancer. P-value < 0.05 was considered as statistically significant.
Figure 8. Correlation analysis for potential hsa-miR-135b-5p-target gene pairs in colorectal cancer. (A) The visual hsa-miR-145-5p-target gene network. (B) The expression correlation of hsa-miR-135b-5p with APC in colorectal cancer. (C) The expression correlation of hsa-miR-135b-5p with FOXO1 in colorectal cancer. (D) The expression correlation of hsa-miR-135b-5p with MBNL1 in colorectal cancer. (E) The expression correlation of hsa-miR-135b-5p with MEF2C in colorectal cancer. (F) The expression correlation of hsa-miR-135b-5p with RECK in colorectal cancer. (G) The expression correlation of hsa-miR-135b-5p with KLF4 in colorectal cancer. (H) The expression correlation of hsa-miR-135b-5p with PPM1E in colorectal cancer. (I) The expression correlation of hsa-miR-135b-5p with ARC in colorectal cancer. (J) The expression correlation of hsa-miR-135b-5p with TTLL7 in colorectal cancer. (K) The expression correlation of hsa-miR-135b-5p with PCP4L1 in colorectal cancer. P-value < 0.05 was considered as statistically significant.
Table 3. The correlation analysis of miR-145-5p and potential target gene in colorectal cancer using starBase.
|miRNA ID||gene ID||R||P-value|
Table 4. The correlation analysis of miR-135b-5p and potential target gene in colorectal cancer using starBase.
|miRNA ID||gene ID||R||P-value|
Figure 9. Correlation analysis for potential circRNA-target gene pairs in colorectal cancer. (A) The expression correlation of hsa_circ_0001955 with CDK4 in colorectal cancer. (B) The expression correlation of hsa_circ_0001955 with MMP12 in colorectal cancer. (C) The expression correlation of hsa_circ_0001955 with RAB3IP in colorectal cancer. (D) The expression correlation of hsa_circ_0000977 with FOXO1 in colorectal cancer. (E) The expression correlation of hsa_circ_0000977 with MBNL1 in colorectal cancer. (F) The expression correlation of hsa_circ_0000977 with MEF2C in colorectal cancer. (G) The expression correlation of hsa_circ_0000977 with RECK in colorectal cancer. (H) The expression correlation of hsa_circ_0000977 with PPM1E in colorectal cancer. (I) The expression correlation of hsa_circ_0000977 with TTLL7 in colorectal cancer. (J) The expression correlation of hsa_circ_0000977 with PCP4L1 in colorectal cancer. P-value < 0.05 was considered as statistically significant.
Figure 10. Construction of a potential circRNA-mediated ceRNA regulatory network in colorectal cancer.
Figure 11. Validation of hsa_circ_0001955-mediated miRNA-mRNA sub-network in CRC using clinical tissue samples. Expression of hsa_circ_0001955 (A), hsa-miR-145-5p (B), CDK4 (C), MMP12 (D) and RAB3IP (E) in CRC cancer samples compared with normal samples. Expression correlation between hsa_circ_0001955/hsa-miR-145-5p (F), hsa-miR-145-5p/CDK4 (G), hsa-miR-145-5p/MMP12 (H), hsa-miR-145-5p/RAB3IP (I), hsa_circ_0001955/CDK4 (J), hsa_circ_0001955/MMP12 (K) and hsa_circ_0001955/RAB3IP (L) pairs in CRC. *p<0.05 represents significant difference.
Figure 12. Validation of hsa_circ_0000977-mediated miRNA-mRNA sub-network in CRC using clinical tissue samples. Expression of hsa_circ_0000977 (A), hsa-miR-135b-5p (B), FOXO1 (C), MBNL1 (D), MEF2C (E), RECK (F), PPM1E (G), TTLL7 (H) and PCP4L1 (I) in CRC cancer samples compared with normal samples. Expression correlation between hsa_circ_0000977/ hsa-miR-135b-5p (J), hsa-miR-135b-5p/FOXO1 (K), hsa-miR-135b-5p/MBNL1 (L), hsa-miR-135b-5p/MEF2C (M), hsa-miR-135b-5p/RECK (N), hsa-miR-135b-5p/PPM1E (O), hsa-miR-135b-5p/TTLL7 (P), hsa-miR-135b-5p/PCP4L1 (Q), hsa_circ_0000977/FOXO1 (R), hsa_circ_0000977/MBNL1 (S), hsa_circ_0000977/MEF2C (T), hsa_circ_0000977/RECK (U), hsa_circ_0000977/PPM1E (V), hsa_circ_0000977/TTLL7 (W) and hsa_circ_0000977/PCP4L1 (X) pairs in CRC. *p<0.05 represents significant difference.
Increasing evidences have demonstrated that circRNAs serve as crucial players in initiation and progression of human cancers, including CRC. In spite of some studies have reported their constructed circRNA-miRNA-mRNA regulatory network in CRC [14, 15]. However, the current knowledge of circRNA-associated ceRNA network in CRC remain inadequate and need to be further investigated. In the present study, by a series of in silico analysis, we identified a potential circRNA-miRNA-mRNA ceRNA regulatory network in CRC pathogenesis.
After performing differential expression analysis, 10 circRNAs were selected for further research. Most of the 10 circRNAs have been found to be closely linked to cancer development and progression. For example, hsa_circRNA_103809 regulates CRC proliferation and migration via miR-532-3p/FOXO4 axis ; hsa_circRNA_101555 is involved in carcinogenesis of breast cancer  hsa_circRNA_102619 played a suppressive role in progression of pancreatic ductal adenocarcinoma . By using circBase and CSCD databases, 9 of 10 circRNAs were finally identified.
Multiple studies have confirmed that circRNAs positively regulate expression of downstream genes by acting as miRNA sponges [19, 20]. Therefore, potential miRNAs of the 9 circRNAs were predicted. Followed by a series of analyses including miRNA expression analysis and correlation analysis between circRNA and miRNA in CRC, two pathways (hsa_circ_0001955/miR-145-5p and hsa_circ_0000977/miR-135b-5p) were regarded as the key axes in carcinogenesis of CRC. The functions of miR-145-5p and miR-135b-5p have been investigated, which were in accordance with ceRNA hypothesis. For example, miR-145-5p inhibited growth, migration and invasion of CRC [21, 22]. Silencing miR-135b-5p sensitized CRC cells to oxaliplatin-induced apoptosis by upregulating FOXO1 expression .
Next, we further explored the detailed downstream molecular mechanisms of hsa_circ_0001955/miR-145-5p and hsa_circ_0000977/miR-135b-5p axes in CRC. As is known to all, miRNAs exert their functions by negatively modulating gene expression [24–27]. Firstly, target genes of miR-145-5p and miR-135b-5p were predicted by a comprehensive database, namely miRNet. Then, the DEGs between CRC tissues and normal tissues were screened. 29 targets were selected by intersection of target gene set and DEG set, after which expression of the 29 targets were validated using TCGA data by starBase. Subsequently, correlation analysis of target gene with miRNA was performed through both GEO and TCGA data. Finally, we found that three targets of miR-145-5p (CDK4, MMP12 and RAB3IP) and seven targets of miR-135b-5p (FOXO1, MBNL1, MEF2C, RECK, PPM1E, TTLL7 and PCP4L1) possessed the most potential. These targets have been well-documented to closely link to cancer progression, including CRC. For instance, miR-143-3p-suppressed CDK4 enhanced CRC growth .
According to ceRNA mechanism, there should be positive expression relationship between circRNA and target gene. Intriguingly, hsa_circ_0001955 expression was positively correlated with expression of CDK4, MMP12 and RAB3IP, and hsa_circ_0000977 was also positively associated with expression of FOXO1, MBNL1, MEF2C, RECK, PPM1E, TTLL7 and PCP4L1 in CRC.
Taken all these data together, we successfully established a circRNA-miRNA-mRNA triple ceRNA sub-network mediated by two circRNAs (hsa_circ_0001955 and hsa_circ_0000977) in CRC. Moreover, the expression levels of RNAs in the network and expression correlation were followingly experimentally validated using clinical CRC samples, partially supporting the analytic accuracy of in silico analysis. The current work provides novel insight into the detailed molecular mechanism of CRC carcinogenesis, and the research route may also be applied to explore the ceRNA molecular interaction mechanism of circRNA-miRNA-mRNA network in initiation and progression of other types of human cancer. However, more efforts should be given to elucidate the function of the established network in CRC with in vitro and in vivo experiments. Meanwhile, the diagnostic and prognostic values of each RNA in the established network should be evaluated using a large number of clinical CRC data in the future, which will help to develop promising diagnostic and prognostic biomarkers for patients with CRC.
In summary, this study reveals a potential hsa_circ_0001955/hsa_circ_0000977-mediated circRNA-miRNA-mRNA ceRNA sub-network in CRC by whole-transcriptome analysis, including differential expression analysis, intersection analysis and correlation analysis. Despite need more experimental and clinical validations, targeting the molecules in this sub-network may represent a promising treatment for patients with CRC.
Materials and Methods
Inclusion of datasets
In this study, we intended to construct a key circRNA-miRNA-mRNA regulatory network in CRC. Firstly, we selected the datasets for analysis using NCBI GEO database (http://www.ncbi.nlm.nih.gov/geo/). Selection criteria: only datasets that focused on human CRC tissue samples were included; datasets regarding CRC cell lines or animal CRC tissue samples were excluded; only datasets with sample count more than ten were included; only datasets containing circRNA expression, miRNA expression and mRNA expression were included. Finally, only dataset GSE126095 met all these criteria mentioned above. GSE126095 had three sub-datasets (GSE126094, GSE126093 and GSE126092) from 10 CRC tissues and 10 matched normal tissues. GSE126094, based on the platform of GPL19978 Agilent-069978 Arraystar Human CircRNA microarray V1, investigated circRNA expression profile; GSE126093, based on the platform of GPL18058 Exiqon miRCURY LNA microRNA array G7, investigated miRNA expression profile; and GSE126092, based on the platform of GPL21047 Agilent-074348 Human LncRNA v6 4X180K.
Differential analysis for circRNA, miRNA and mRNA
Differential analysis for GSE126094, GSE126093 and GSE126092 was successively conducted to obtain differentially expressed circRNAs (DECs), differentially expressed miRNAs (DEmiRNAs) and differentially expressed genes (DEGs) between CRC tissues and matched normal tissues by using an online tool GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r/) which is provided by NCBI GEO database . Adjust P-vlaue (adj.P.Val) < 0.05 and |log2FC (Fold change)| > 1 were set as the thresholds for identifying DECs, DEmiRNAs and DEGs.
circBase and CSCD (Cancer-Specific CircRNA Database) analysis
circBase, a database for circRNA-associated studies, provides scripts to acquire known and novel circRNAs in sequencing data . circBase was employed to identify parental genes of potential circRNAs. CSCD (http://gb.whu.edu.cn/CSCD/), which is an online tool to investigate cancer-specific circRNAs, was utilized to obtain structure of potential circRNAs.
miRNAs that might bind to circRNAs were predicted by starBase database (http://starbase.sysu.edu.cn/). starBase database is an open-source platform that provides the ncRNA interactions from CLIP-seq, degradome-seq and RNA-RNA interactome data [31, 32]. The circBase ID of circRNA was entered into starBase database, after which potential miRNAs of circRNA were directly presented on the webpage.
Validation of miRNA expression
The miRNAs that commonly appeared in target miRNAs set and DEmiRNAs set were regarded as potential miRNAs of circRNAs in CRC. Expression of these miRNAs was further validated using TCGA data by starBase database [31, 32]. P < 0.05 was considered as statistically significant.
Target gene prediction
miRNet (https://www.mirnet.ca/), an integrated platform linking miRNAs, targets and functions, was introduced to perform miRNA’s target gene prediction . The miRNA-gene interactions were directly downloaded from the webpage.
Validation of gene expression
The genes that commonly appeared in target genes set and DEGs set were considered as potential target genes of miRNAs in CRC. Expression of these target genes were further confirmed using TCGA data using starBase database [31, 32]. P < 0.05 was considered as statistically significant.
Correlation analysis for RNA-RNA interactions
The expression correlation of hsa_circ_0001955 or hsa_circ_0000977 with miRNA in CRC was evaluated using the data from GSE126094 and GSE126093; the expression correlation of miR-145-5p or miR-135b-5p with target gene in CRC was assessed using the data from GSE126093 and GSE126092; the expression correlation of hsa_circ_0001955 or hsa_circ_0000977 with target gene in CRC was determined using the data from GSE126094 and GSE126092. Then, starbase database was also employed to validate the expression correlation of miR-145-5p or miR-135b-5p with target gene. |R| > 0.1 and P-value < 0.05 were set as the thresholds for identifying significant RNA-RNA interactions.
Clinical tissues and qRT-PCR
Fresh frozen CRC cancer tissues and adjacent normal tissues were obtained from 20 CRC patients who received surgical treatment at the First Affiliated Hospital of Zhejiang University (Hangzhou, China), and were frozen and stored in liquid nitrogen. This study was approved by the Ethical Committee of First Affiliated Hospital of Zhejiang University. Each informed consent was acquired from every patient. The expression of potential circRNAs, miRNAs and genes in these clinical samples was detected by the method of qRT-PCR as previously described . The Primers used in this study were designed and purchased from RiboBio Co. Ltd (Guangzhou, China).
Most of statistical analyses have been calculated by the online databases or tools. Expression correlation for RNA-RNA interactions was perform by Person correlation coefficient using Graphpad Prism software (version 7). Paired Student’s t-test was used to evaluate expression differences of circRNA, miRNA or gene between normal group and cancer group. Results with P-value < 0.05 were considered as statistically significant.
This study was approved by the Ethical Committee of First Affiliated Hospital of Zhejiang University. Each informed consent was acquired from every patient.
All authors contributed to the conception and design of this manuscript, collection and assembly of data, manuscript writing, and final approval.
Thanks for the anonymous reviewers for their valuable comments and suggestions that helped improve the quality of our manuscript.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work was supported by National Natural Science Foundation of China (81572987, 81874225).
- 1. Siegel R, Desantis C, Jemal A. Colorectal cancer statistics, 2014. CA Cancer J Clin. 2014; 64:104–17. https://doi.org/10.3322/caac.21220 [PubMed]
- 2. Shen T, Cheng X, Liu X, Xia C, Zhang H, Pan D, Zhang X, Li Y. Circ_0026344 restrains metastasis of human colorectal cancer cells via miR-183. Artif Cells Nanomed Biotechnol. 2019; 47:4038–45. https://doi.org/10.1080/21691401.2019.1669620 [PubMed]
- 3. Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, Sun W, Dou K, Li H. Circular RNA: A new star of noncoding RNAs. Cancer Lett. 2015; 365:141–48. https://doi.org/10.1016/j.canlet.2015.06.003 [PubMed]
- 4. Hu X, Wu D, He X, Zhao H, He Z, Lin J, Wang K, Wang W, Pan Z, Lin H, Wang M. circGSK3β promotes metastasis in esophageal squamous cell carcinoma by augmenting β-catenin signaling. Mol Cancer. 2019; 18:160. https://doi.org/10.1186/s12943-019-1095-y [PubMed]
- 5. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011; 146:353–58. https://doi.org/10.1016/j.cell.2011.07.014 [PubMed]
- 6. Xue D, Wang H, Chen Y, Shen D, Lu J, Wang M, Zebibula A, Xu L, Wu H, Li G, Xia L. Circ-AKT3 inhibits clear cell renal cell carcinoma metastasis via altering miR-296-3p/E-cadherin signals. Mol Cancer. 2019; 18:151. https://doi.org/10.1186/s12943-019-1072-5 [PubMed]
- 7. Han D, Li J, Wang H, Su X, Hou J, Gu Y, Qian C, Lin Y, Liu X, Huang M, Li N, Zhou W, Yu Y, Cao X. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017; 66:1151–64. https://doi.org/10.1002/hep.29270 [PubMed]
- 8. Wang C, Tao W, Ni S, Chen Q. Circular RNA circ-Foxo3 induced cell apoptosis in urothelial carcinoma via interaction with miR-191-5p. Onco Targets Ther. 2019; 12:8085–94. https://doi.org/10.2147/OTT.S215823 [PubMed]
- 9. Xie H, Ren X, Xin S, Lan X, Lu G, Lin Y, Yang S, Zeng Z, Liao W, Ding YQ, Liang L. Emerging roles of circRNA_001569 targeting miR-145 in the proliferation and invasion of colorectal cancer. Oncotarget. 2016; 7:26680–91. https://doi.org/10.18632/oncotarget.8589 [PubMed]
- 10. Hsiao KY, Lin YC, Gupta SK, Chang N, Yen L, Sun HS, Tsai SJ. Noncoding Effects of Circular RNA CCDC66 Promote Colon Cancer Growth and Metastasis. Cancer Res. 2017; 77:2339–50. https://doi.org/10.1158/0008-5472.CAN-16-1883 [PubMed]
- 11. Lin J, Cai D, Li W, Yu T, Mao H, Jiang S, Xiao B. Plasma circular RNA panel acts as a novel diagnostic biomarker for colorectal cancer. Clin Biochem. 2019; 74:60–68. https://doi.org/10.1016/j.clinbiochem.2019.10.012 [PubMed]
- 12. Ye DX, Wang SS, Huang Y, Chi P. A 3-circular RNA signature as a noninvasive biomarker for diagnosis of colorectal cancer. Cancer Cell Int. 2019; 19:276. https://doi.org/10.1186/s12935-019-0995-7 [PubMed]
- 13. Liu K, Guo Y, Zheng K, Zou C, Wu H, Wang S, Ou L, Wang Y, Huang B, Wang X. Identification of the circRNA-miRNA-mRNA regulatory network of Hsp90 inhibitor-induced cell death in colorectal cancer by integrated analysis. Gene. 2020; 727:144232. https://doi.org/10.1016/j.gene.2019.144232 [PubMed]
- 14. Yuan W, Peng S, Wang J, Wei C, Ye Z, Wang Y, Wang M, Xu H, Jiang S, Sun D, Dai C, Jiang L, Li X. Identification and characterization of circRNAs as competing endogenous RNAs for miRNA-mRNA in colorectal cancer. PeerJ. 2019; 7:e7602. https://doi.org/10.7717/peerj.7602 [PubMed]
- 15. Song W, Fu T. Circular RNA-Associated Competing Endogenous RNA Network and Prognostic Nomogram for Patients With Colorectal Cancer. Front Oncol. 2019; 9:1181. https://doi.org/10.3389/fonc.2019.01181 [PubMed]
- 16. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, Yefanov A, Lee H, Zhang N, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013; 41:D991–95. https://doi.org/10.1093/nar/gks1193 [PubMed]
- 17. Glažar P, Papavasileiou P, Rajewsky N. circBase: a database for circular RNAs. RNA. 2014; 20:1666–70. https://doi.org/10.1261/rna.043687.113 [PubMed]
- 18. Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014; 42:D92–97. https://doi.org/10.1093/nar/gkt1248 [PubMed]
- 19. Yang JH, Li JH, Shao P, Zhou H, Chen YQ, Qu LH. starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 2011; 39:D202–09. https://doi.org/10.1093/nar/gkq1056 [PubMed]
- 20. Fan Y, Siklenka K, Arora SK, Ribeiro P, Kimmins S, Xia J. miRNet - dissecting miRNA-target interactions and functional associations through network-based visual analysis. Nucleic Acids Res. 2016; 44:W135–41. https://doi.org/10.1093/nar/gkw288 [PubMed]
- 21. Lou W, Liu J, Ding B, Chen D, Xu L, Ding J, Jiang D, Zhou L, Zheng S, Fan W. Identification of potential miRNA-mRNA regulatory network contributing to pathogenesis of HBV-related HCC. J Transl Med. 2019; 17:7. https://doi.org/10.1186/s12967-018-1761-7 [PubMed]
- 22. Bian L, Zhi X, Ma L, Zhang J, Chen P, Sun S, Li J, Sun Y, Qin J. Hsa_circRNA_103809 regulated the cell proliferation and migration in colorectal cancer via miR-532-3p / FOXO4 axis. Biochem Biophys Res Commun. 2018; 505:346–52. https://doi.org/10.1016/j.bbrc.2018.09.073 [PubMed]
- 23. Afzali F, Salimi M. Unearthing Regulatory Axes of Breast Cancer circRNAs Networks to Find Novel Targets and Fathom Pivotal Mechanisms. Interdiscip Sci. 2019; 11:711–22. https://doi.org/10.1007/s12539-019-00339-6 [PubMed]
- 24. Huang WJ, Wang Y, Liu S, Yang J, Guo SX, Wang L, Wang H, Fan YF. Silencing circular RNA hsa_circ_0000977 suppresses pancreatic ductal adenocarcinoma progression by stimulating miR-874-3p and inhibiting PLK1 expression. Cancer Lett. 2018; 422:70–80. https://doi.org/10.1016/j.canlet.2018.02.014 [PubMed]
- 25. Liu Z, Yu Y, Huang Z, Kong Y, Hu X, Xiao W, Quan J, Fan X. CircRNA-5692 inhibits the progression of hepatocellular carcinoma by sponging miR-328-5p to enhance DAB2IP expression. Cell Death Dis. 2019; 10:900. https://doi.org/10.1038/s41419-019-2089-9 [PubMed]
- 26. Liu Z, Pan HM, Xin L, Zhang Y, Zhang WM, Cao P, Xu HW. Circ-ZNF609 promotes carcinogenesis of gastric cancer cells by inhibiting miRNA-145-5p expression. Eur Rev Med Pharmacol Sci. 2019; 23:9411–17. https://doi.org/10.26355/eurrev_201911_19433 [PubMed]
- 27. Shen X, Jiang H, Chen Z, Lu B, Zhu Y, Mao J, Chai K, Chen W. MicroRNA-145 Inhibits Cell Migration and Invasion in Colorectal Cancer by Targeting TWIST. Onco Targets Ther. 2019; 12:10799–809. https://doi.org/10.2147/OTT.S216147 [PubMed]
- 28. Niu Y, Zhang J, Tong Y, Li J, Liu B. miR-145-5p restrained cell growth, invasion, migration and tumorigenesis via modulating RHBDD1 in colorectal cancer via the EGFR-associated signaling pathway. Int J Biochem Cell Biol. 2019; 117:105641. https://doi.org/10.1016/j.biocel.2019.105641 [PubMed]
- 29. Qin Y, Li L, Wang F, Zhou X, Liu Y, Yin Y, Qi X. Knockdown of Mir-135b Sensitizes Colorectal Cancer Cells to Oxaliplatin-Induced Apoptosis Through Increase of FOXO1. Cell Physiol Biochem. 2018; 48:1628–37. https://doi.org/10.1159/000492284 [PubMed]
- 30. Lou W, Chen J, Ding B, Chen D, Zheng H, Jiang D, Xu L, Bao C, Cao G, Fan W. Identification of invasion-metastasis-associated microRNAs in hepatocellular carcinoma based on bioinformatic analysis and experimental validation. J Transl Med. 2018; 16:266. https://doi.org/10.1186/s12967-018-1639-8 [PubMed]
- 31. Lou W, Ding B, Xu L, Fan W. Construction of Potential Glioblastoma Multiforme-Related miRNA-mRNA Regulatory Network. Front Mol Neurosci. 2019; 12:66. https://doi.org/10.3389/fnmol.2019.00066 [PubMed]
- 32. Lou W, Liu J, Ding B, Xu L, Fan W. Identification of chemoresistance-associated miRNAs in breast cancer. Cancer Manag Res. 2018; 10:4747–57. https://doi.org/10.2147/CMAR.S172722 [PubMed]
- 33. Zhu X, Ma SP, Yang D, Liu Y, Wang YP, Lin T, Li YX, Yang SH, Zhang WC, Wang XL. miR-142-3p Suppresses Cell Growth by Targeting CDK4 in Colorectal Cancer. Cell Physiol Biochem. 2018; 51:1969–81. https://doi.org/10.1159/000495721 [PubMed]