Bioinformatics analysis-based screening of circRNA gene with mainstream expression trend in colorectal cancer and construction of a coexpression regulatory network

PLoS One. 2023 Dec 8;18(12):e0295126. doi: 10.1371/journal.pone.0295126. eCollection 2023.

Abstract

Objective: Since circRNA can be utilized as a potential diagnostic marker for cancer, to explore the regulatory mechanism of colorectal cancer (CRC) using bioinformatics, the public database of circRNA was mined.

Methods: CRC differentially expressed miRNAs were screened in the Cancer Genome Atlas (TCGA) database, CRC differentially expressed circRNAs were searched in the Gene Expression Omnibus (GEO) database, the two databases were combined to identify CRC differentially expressed mRNAs, and a circRNA-miRNA‒mRNA regulatory network was constructed by combining a plurality of target prediction databases to identify key genes. The upstream circRNA and regulatory axis of the key genes were identified for gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis to explore the biological functions of circRNA in CRC using the regulatory axis.

Results: After the screening of the GSE21815 dataset, a total of 22 differentially expressed circRNAs were obtained, with 12 upregulated and 10 downregulated genes. Similarly, the GSE126094 dataset yielded 104 differentially expressed circRNAs, comprising 56 upregulated and 48 downregulated genes. Among the differentially expressed circRNAs, five were identified, with VDAC3 and SETD2 showing downregulated expression, while RAD23B, RPPH1, and MYBL2 exhibited upregulated expression. Following the selection process, five DEcircRNAs, eight target miRNAs, and 105 target DEmRNAs were identified. The protein-protein interaction (PPI) network revealed close relationships among the mRNAs, with E2F2, E2F3, CCND1, TNRC6A, and KAT2B identified as key genes. Notably, CCND1 emerged as a critical gene in the PPI network. Through the upregulation of has-circ-0087862, which binds to miR-892b, the translation inhibition of CCND1 by miR-892b was attenuated, leading to enhanced CCND1 expression. Functional enrichment analysis indicated that CCND1 was involved in protein binding and positive regulation of cellular processes, among other functions.

Conclusion: The differentially expressed genes (DEGs) in CRC markedly affected the survival time of patients. CircRNAs could be utilized as diagnostic markers of CRC, and the key genes in CRC could be screened out by bioinformatics, which would be helpful to understand the drug targets for the treatment of human immunodeficiency virus (HIV)-related CRC patients.

MeSH terms

  • Colorectal Neoplasms* / genetics
  • Colorectal Neoplasms* / metabolism
  • Computational Biology
  • Early Detection of Cancer
  • Gene Regulatory Networks
  • Humans
  • MicroRNAs* / genetics
  • RNA, Circular / genetics
  • RNA, Messenger / genetics

Substances

  • RNA, Circular
  • MicroRNAs
  • RNA, Messenger

Grants and funding

This work was supported by the National Natural Science Foundation of China (NO.81470881 and NO.8217112503), the Jiangsu Commission of Health (LGY2017031 and BRA2015473). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. There was no additional external funding received for this study.