Study on the function and prognosis of circular RNA in colorectal cancer tissues based on high-throughput sequencing

doi:10.3969/j.issn.1674-8115.2021.02.010

Abstract

Abstract: Objective

·To reveal the mechanism of colorectal cancer and the association between differential expression of circular RNAs (circRNAs) and patients' prognosis.

Methods

·Colorectal cancer tissues and normal tumor-adjacent tissues were harvested from 12 patients undergoing radical resection of colorectal cancer in the Department of General Surgery, Shanghai General Hospital, Shanghai Jiao Tong University, and the expressions of circRNAs were detected by high-throughput sequencing. Function predictions of differentially expressed circRNAs were carried out by gene ontology (GO) annotation analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Two circRNAs were selected from the significantly up-regulated and down-regulated circRNAs respectively, and the accuracy of sequencing results was verified in 44 sets of colorectal cancer tissues by real-time fluorescent quantitative PCR (RT-qPCR). The circRNA-miRNA targeting binding network map was constructed by searching the public database for the four circRNAs. Finally, the relationships between differential expression of circRNAs and prognosis were determined based on the patients' clinical data.

Results

·A total of 20 924 differentially expressed circRNAs were detected after high-throughput sequencing. Compared with normal colorectal tissues, there were 373 circRNAs (︱log₂FC︱>1, P<0.05) abnormally and significantly expressed in colorectal cancer tissues, of which 243 were significantly up-regulated and 130 were significantly down-regulated. In 44 pairs of colorectal cancer samples, the results verified by RT-qPCR were basically consistent with the sequencing data. After the circRNA-miRNA-mRNA interactions being plotted, it was found that miRNA and target genes associated with circRNA may play an important role in the pathogenesis of colorectal cancer. Combined with the prognosis information, the patients with high expression of hsa_circ_0023984 or hsa_circ_0008192 had poorer prognosis while patients with high expression of hsa_circ_0020093 or hsa_circ_0069922 had better prognosis.

Conclusion

·Abnormally expressed circRNAs in colorectal cancer tissues lay a foundation for further studies on the occurrence and development of colorectal cancer, and provide valuable information for the development of new clinical diagnostic markers and treatments in the future.

Key words: colorectal cancer, circular RNA (circRNA), prognosis, regulatory network

CLC Number:

R735.3

Zhong-mao FU, Zai LUO, Ze-yin RONG, Jian-ming ZHANG, Teng-fei LI, Zhi-long YU, Chen HUANG. Study on the function and prognosis of circular RNA in colorectal cancer tissues based on high-throughput sequencing[J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2021, 41(2): 187-195.

Figures/Tables 6

TrendMD

References 29

1	Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2018: globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries[J]. CA Cancer J Clin, 2018, 68(6): 394-424.
2	Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020[J]. CA Cancer J Clin, 2020, 70(1): 7-30.
3	Chen W, Zheng R, Baade PD, et al. Cancer statistics in China, 2015[J]. CA Cancer J Clin, 2016, 66(2): 115-132.
4	Zhang H, Wang Z, Ma R, et al. MicroRNAs as biomarkers for the progression and prognosis of colon carcinoma[J]. Int J Mol Med, 2018, 42(4): 2080-2088.
5	Sun K, Han RJ, Han Y, et al. Accuracy of combined computed tomography colonography and dual energy iiodine map imaging for detecting colorectal masses using high-pitch dual-source CT[J]. Sci Rep, 2018, 8(1): 3790.
6	Thomsen M, Skovlund E, Sorbye H, et al. Prognostic role of carcinoembryonic antigen and carbohydrate antigen 19-9 in metastatic colorectal cancer: a BRAF-mutant subset with high CA 19-9 level and poor outcome[J]. Br J Cancer, 2018, 118(12): 1609-1616.
7	Gao Y, Wang J, Zhou Y, et al. Evaluation of serum CEA, CA19-9, CA72-4, CA125 and ferritin as diagnostic markers and factors of clinical parameters for colorectal cancer[J]. Sci Rep, 2018, 8(1): 2732.
8	Chen LL, Yang L. Regulation of circRNA biogenesis[J]. RNA Biol, 2015, 12(4): 381-388.
9	Memczak S, Jens M, Elefsinioti A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency[J]. Nature, 2013, 495(7441): 333-338.
10	Pamudurti NR, Bartok O, Jens M, et al. Translation of circRNAs[J]. Mol Cell, 2017, 66(1): 9-21.e7.
11	Barrett SP, Salzman J. Circular RNAs: analysis, expression and potential functions[J]. Dev Camb Engl, 2016, 143(11): 1838-1847.
12	Zhu Z, Rong Z, Luo Z, et al. Circular RNA circNHSL1 promotes gastric cancer progression through the miR-1306-3p/SIX1/vimentin axis[J]. Mol Cancer, 2019, 18(1): 126.
13	Siegel RL, Miller KD, Goding Sauer A, et al. Colorectal cancer statistics, 2020[J]. CA Cancer J Clin, 2020, 70(3): 145-164.
14	Wang Y, Wang Z. Efficient backsplicing produces translatable circular mRNAs[J]. RNA, 2015, 21(2): 172-179.
15	Jeck WR, Sorrentino JA, Wang K, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats[J]. RNA, 2013, 19(2): 141-157.
16	Schwanhäusser B, Busse D, Li N, et al. Global quantification of mammalian gene expression control[J]. Nature, 2011, 473(7347): 337-342.
17	Vo JN, Cieslik M, Zhang YJ, et al. The landscape of circular RNA in cancer[J]. Cell, 2019, 176(4): 869-881.e13.
18	Zhu Z, Yu Z, Rong Z, et al. The novel GINS4 axis promotes gastric cancer growth and progression by activating Rac1 and CDC42[J]. Theranostics, 2019, 9(26): 8294-8311.
19	Li Y, Zheng QP, Bao CY, et al. Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis[J]. Cell Res, 2015, 25(8): 981-984.
20	Zhu P, Zhu X, Wu J, et al. IL-13 secreted by ILC2s promotes the self-renewal of intestinal stem cells through circular RNA circPan3[J]. Nat Immunol, 2019, 20(2): 183-194.
21	Li WL, Xu YQ, Wang XD, et al. circCCT3 modulates vascular endothelial growth factor A and Wnt signaling to enhance colorectal cancer metastasis through sponging miR-613[J]. DNA Cell Biol, 2020, 39(1): 118-125.
22	Zhong Z, Huang M, Lv M, et al. Circular RNA MYLK as a competing endogenous RNA promotes bladder cancer progression through modulating VEGFA/VEGFR2 signaling pathway[J]. Cancer Lett, 2017, 403: 305-317.
23	Liang WC, Wong CW, Liang PP, et al. Translation of the circular RNA circβ-catenin promotes liver cancer cell growth through activation of the Wnt pathway[J]. Genome Biol, 2019, 20(1): 84.
24	Zhi X, Zhang J, Cheng Z, et al. circLgr4 drives colorectal tumorigenesis and invasion through Lgr4-targeting peptide[J]. Int J Cancer, 2019. DOI:10.1002/ijc.32549.
25	Rong DW, Sun HD, Li ZX, et al. An emerging function of circRNA-miRNAs-mRNA axis in human diseases[J]. Oncotarget, 2017, 8(42): 73271-73281.
26	Gallí M, van Gool F, Rongvaux A, et al. The nicotinamide phosphoribosyltransferase: a molecular link between metabolism, inflammation, and cancer[J]. Cancer Res, 2010, 70(1): 8-11.
27	Garten A, Schuster S, Penke M, et al. Physiological and pathophysiological roles of NAMPT and NAD metabolism[J]. Nat Rev Endocrinol, 2015, 11(9): 535-546.
28	Sampath D, Zabka TS, Misner DL, et al. Inhibition of nicotinamide phosphoribosyltransferase (NAMPT) as a therapeutic strategy in cancer[J]. Pharmacol Ther, 2015, 151: 16-31.
29	Lucena-Cacace A, Otero-Albiol D, Jiménez-García MP, et al. NAMPT is a potent oncogene in colon cancer progression that modulates cancer stem cell properties and resistance to therapy through Sirt1 and PARP[J]. Clin Cancer Res, 2018, 24(5): 1202-1215.

Gene	Primer sequence （5'→3'）
hsa_circ_0023984
Forward	5'-GTCAGTATATTACTCTACATTGTCCCAG-3'
Reverse	5'-TCTAAATTAGTTTGATACTTCAGCAGC-3'
hsa_circ_0008192
Forward	5'-TGTCAGTGCCGTCCAATC-3'
Reverse	5'-GCGAACAGGGAGGTCTAT-3'
hsa_circ_0020093
Forward	5'-CGGTTACTGTGACCTGACTGGAG -3'
Reverse	5'-TGGGCCATCTGTTAAATATCCAG-3'
hsa_circ_0069922
Forward	5'-GCTGATGAATCCAGCAGCAGC-3'
Reverse	5'-TCCAGGACAGCTTCATCCTCC-3'
GAPDH
Forward	5'-GTCTCCTCTGACTTCAACAGCG-3'
Reverse	5'-ACCACCCTGTTGCTGTAGCCAA-3'

Gene	Primer sequence （5'→3'）
hsa_circ_0023984
Forward	5'-GTCAGTATATTACTCTACATTGTCCCAG-3'
Reverse	5'-TCTAAATTAGTTTGATACTTCAGCAGC-3'
hsa_circ_0008192
Forward	5'-TGTCAGTGCCGTCCAATC-3'
Reverse	5'-GCGAACAGGGAGGTCTAT-3'
hsa_circ_0020093
Forward	5'-CGGTTACTGTGACCTGACTGGAG -3'
Reverse	5'-TGGGCCATCTGTTAAATATCCAG-3'
hsa_circ_0069922
Forward	5'-GCTGATGAATCCAGCAGCAGC-3'
Reverse	5'-TCCAGGACAGCTTCATCCTCC-3'
GAPDH
Forward	5'-GTCTCCTCTGACTTCAACAGCG-3'
Reverse	5'-ACCACCCTGTTGCTGTAGCCAA-3'

circRNA ID	FC	P value	Regulation	circRNA chromosome	Gene
hsa_circ_0079534	14.24	<0.001	Up	NC_000007.14:20154193_20161870_-	MACC1
hsa_circ_0023984	10.43	<0.001	Up	NC_000011.10:89400017_89402542_-	NOX4
hsa_circ_0008192	10.27	<0.001	Up	NC_000009.12:112262435_112297916_-	PTBP3
hsa_circ_0008694	8.55	<0.001	Up	NC_000020.11:51712966_51729978_-	ATP9A
hsa_circ_0004705	7.40	<0.001	Up	NC_000003.12:14443624_14447816_+	SLC6A6
hsa_circ_0020093	-8.27	<0.001	Down	NC_000010.11:115120185_115129535_+	ATRNL1
hsa_circ_0069922	-7.22	<0.001	Down	NC_000004.12:71236576_71255399_+	SLC4A4
hsa_circ_0073244	-6.54	<0.001	Down	NC_000005.10:84060300_84066606_-	EDIL3
hsa_circ_0005720	-5.31	0.001	Down	NC_000001.11:95143891_95173889_+	TMEM56
hsa_circ_0008676	-5.26	0.002	Down	NC_000004.12:71339370_71357187_+	SLC4A4

circRNA ID	FC	P value	Regulation	circRNA chromosome	Gene
hsa_circ_0079534	14.24	<0.001	Up	NC_000007.14:20154193_20161870_-	MACC1
hsa_circ_0023984	10.43	<0.001	Up	NC_000011.10:89400017_89402542_-	NOX4
hsa_circ_0008192	10.27	<0.001	Up	NC_000009.12:112262435_112297916_-	PTBP3
hsa_circ_0008694	8.55	<0.001	Up	NC_000020.11:51712966_51729978_-	ATP9A
hsa_circ_0004705	7.40	<0.001	Up	NC_000003.12:14443624_14447816_+	SLC6A6
hsa_circ_0020093	-8.27	<0.001	Down	NC_000010.11:115120185_115129535_+	ATRNL1
hsa_circ_0069922	-7.22	<0.001	Down	NC_000004.12:71236576_71255399_+	SLC4A4
hsa_circ_0073244	-6.54	<0.001	Down	NC_000005.10:84060300_84066606_-	EDIL3
hsa_circ_0005720	-5.31	0.001	Down	NC_000001.11:95143891_95173889_+	TMEM56
hsa_circ_0008676	-5.26	0.002	Down	NC_000004.12:71339370_71357187_+	SLC4A4

[1]	Sheng CHENG, Yi ZHAO, Yong-chen WANG, Ping HUANG. Meta-analysis of the effect of BRAF gene mutation on the prognosis of colorectal cancer patients with liver metastases after hepatectomy [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2021, 41(6): 786-792.
[2]	Xu-guang CHEN, Sheng-yi SHI, Lan HU, Yu CHEN, Yi-ming LU, Jing YE. Evaluative value of plasma fibrin degradation product in early prognosis of patients with hemorrhagic stroke [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2021, 41(5): 612-616.
[3]	Ru-juan BAO, Hui-fang CHEN, Yu DONG, You-qiong YE, Bing SU. Construction of prognostic risk score model of colorectal cancer gene signature based on transcriptome dysregulation [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2021, 41(3): 285-296.
[4]	Yan-ru MA, Lin-hua JI, Tian-ying TONG, Yu-qing YAN, Chao-qin SHEN, Xin-yu ZHANG, Ying-ying CAO, Jie HONG, Hao-yan CHEN. Establishment and validation of prognostic prediction model of colorectal cancer based on single-cell RNA sequencing [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2021, 41(2): 159-165.
[5]	Pei-kun HU, Jie HE, Lian-ming WU, Heng GE, Jian-rong XU, Jun PU. Effect of microvascular obstruction on left ventricle function and prognosis in patients with ST-segment elevation myocardial infarction [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2021, 41(2): 173-179.
[6]	TANG Dong-juan, XUE Xiao-mei, HE Bin. Early diagnosis and prognostic evaluation value of miR-133a in patients with acute myocardial infarction [J]. , 2020, 40(3): 339-.
[7]	YE Xin, ZHOU Xiao-yun, YANG Li, WANG Jie, HE Qi. Analysis of clinicopathological features and prognosis of young patients with breast cancer in different age groups [J]. , 2020, 40(3): 351-.
[8]	LIU Jin-yan, ZHAO Jun-tao, WANG Heng-ru, HUANG Jia-yi, WANG Yu-ting, XIANG Ming-jie. Effect of cytoskeleton-associated protein 4 on tumor diagnosis and prognosis [J]. , 2020, 40(2): 255-.
[9]	HE Chun-ming, YIN Hang, TANG Jian, DING Yi-zong, FU Yu-jie, ZHAO Xiao-jing. Establishment and internal validation of a prognostic model for elderly patients after lung cancer surgery based on SEER [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2020, 40(11): 1554-1561.
[10]	JIANG Ying-ying1, 2, QIN Xing1, ZHANG Jian-jun1, CHEN Wan-tao1. Expression of long noncoding RNA COL11A1-208 in oral squamous cell carcinoma and its clinical significance [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2020, 40(10): 1334-1339.
[11]	ZHANG Feng-chun, ZHANG Shuo-yuan, CHEN Tian-en, XU Ying-chun. Clinical application of PD-1/PD-L1 to prognosis prediction and treatment of triple-negative breast cancer [J]. , 2020, 40(1): 128-.
[12]	YAN Yu-qing, SHEN Chao-qin, CHEN Hao-yan, HONG Jie#, WANG Zhen-hua#. Expression of lnc-MTBP-5 in colorectal cancer and its effect on cell invasion [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2020, 40(09): 1193-1201.
[13]	LI Chao, MI Jian-qing, WANG Jin. Advances in Philadelphia chromosome-like acute lymphoblastic leukemia [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2020, 40(09): 1294-1301.
[14]	WANG Yang-yang1, ZHENG Xiao2, ZHU Chun-chao1, LIU Qiang3, ZHANG Zi-zhen1, ZHAO En-hao1. Analysis of Tim-3/Gal-9 expression and T cell infiltration in adenocarcinoma of the esophagogastric junction [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2020, 40(06): 847-855.
[15]	WANG Xiao-xia1, XU Tian-xiang2, LÜ Tang-feng3, XU Wu-jian3, LIU Jian-bo1, SONG Yong3. Clinical significance of geranylgeranyl diphosphate synthase expression levels in patients with acute respiratory distress syndrome [J]. JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY (MEDICAL SCIENCE), 2020, 40(05): 639-645.