The integrative multi-omics approach identifies the novel competing endogenous RNA (ceRNA) network in colorectal cancer

Abstract

Circular RNAs (circRNA) are known to function as competing endogenous RNA (ceRNA) in various cancers by regulating microRNAs (miRNA). However, in colorectal cancer (CRC), the precise pathological role of circ000240/miRNA/mRNA remains indeterminate. The expression level of hsa_circ_000240 was evaluated using qRT-PCR in matching pairs of CRC tumor and adjacent normal tissue samples in our laboratory. Then, to determine whether hsa_circ_000240 acted as a ceRNA in CRC, the linked miRNAs and gene targets were retrieved. Topological analysis of candidate genes using a network approach identified the most critical hub genes and subnetworks related to CRC disease. Microarray and bulk RNA sequencing analyses were utilized to comprehensively evaluate the expression levels of both miRNA and mRNA in CRC. Single-cell RNA-seq analysis was also used to evaluate the significant overall survival (OS) genes at the cellular level. ATAC-seq data provided insights into candidate genes' accessible chromatin regions. The research uncovered a considerable upregulation of hsa_circ_000240 in CRC tissues. Three miRNAs interacted with the target circRNA. One thousand six hundred eighty intersected genes regulated by three miRNAs were further identified, and the relevant functionality of identified neighbor genes highlighted their relevance to cancer. The topological analysis of the constructed network has identified 33 hub genes with notably high expression in CRC. Among these genes, eight, including CHEK1, CDC6, FANCI, GINS2, MAD2L1, ORC1, RACGAP1, and SMC4, have demonstrated a significant impact on overall survival. The utilization of single-cell RNA sequencing unequivocally corroborated the augmented expression levels of CDC6 and ORC1 in individuals with CRC, alongside their noteworthy connection with the infiltration of immune cells. ATAC-seq analyses revealed altered accessibility regions in Chr2, 4, and 12 for CDC6 and ORC1 high-expression. Correlation analysis of CDC6 and ORC1 further highlighted the association of candidate gene expression with exhaustion markers such as CTLA4, CD247, TIGIT, and CD244. The candidate genes exhibit a positive correlation with chromatin remodeling and histone acetylation. These epigenetic modifications play a significant role in influencing the cancer progression following expression of CDC6 and ORC1 in CRC. Additionally, results showed that the methylation rate of the promoter region of CDC6 was elevated in CRC disease, confirming the functional importance of CDC6 and their interaction with hsa_circ_000240 and associated ceRNA in CRC. In conclusion, this study highlights hsa_circ_000240's role as a ceRNA in CRC. It opens new avenues for further dissection of CDC6, ORC1, and underlying novel epigenetics and immunotherapy targets for CRC therapy.