Objective: This study aimed to discover the ceRNAs network in the pathophysiological development of human colorectal cancer (CRC) and to screen biomarkers for target therapy and prognosis by using integrated bioinformatics analysis.
Methods: Data on gene expressions of mRNAs, miRNAs, and circRNAs and clinical information were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus databases, respectively. Differentially expressed mRNAs (DEmRNAs) were identified by using the DESeq2 package of R software. Functional enrichment analysis was conducted using the ClusterProfiler package of R software. The protein-protein interaction (PPI) network was shown by the STRING website. Survival analysis of hub genes was performed using the survival package in R software. Interactions among hub genes, differentially expressed miRNAs (DEmiRNAs), and differentially expressed circRNAs (DEcircRNAs) were used to construct the ceRNAs network.
Results: A total of 412 DEmRNAs including 82 upregulated and 330 downregulated genes were screened out between 473 CRC and 41 normal samples. Two hundred and sixty DEcircRNAs including 253 upregulated and 7 downregulated genes were altered between 23 CRC and 23 normal samples. One hundred and ninety DEmiRNAs including 82 upregulated and 108 downregulated genes were obtained between 450 CRC and 8 normal samples. A ceRNAs and PPI network were successfully constructed, and TIMP1 associated with prognosis was employed.
Conclusion: The present study identified a novel circRNAs-miRNAs-mRNA ceRNAs network, which implied that TIMP1 and related miRNAs, circRNAs were potential biomarkers underlying the development of CRC, providing new insights for survival predictions and therapeutic targets.
Keywords: TIMP1; bioinformatic analysis; ceRNAs; circRNAs; colorectal cancer.
© 2021 Qin et al.