Identification of functional circRNA/miRNA/mRNA regulatory network for exploring prospective therapy strategy of colorectal cancer

Circular RNA (circRNA) has been reported to have great scientific significance and clinical value in multiple cancers including colorectal cancer (CRC). However, the biological function of most circRNAs in CRC is still in its infancy. Herein, we discovered the differential expressed circRNAs (DECs) between CRC tissues and matched adjacent using deep RNA sequencing and further confirmed the DECs expression by combining with another Gene Expression Omnibus dataset. Furthermore, we validated the expression of the top four upregulated circRNAs (hsa_circ_0030632, hsa_circ_0004887, hsa_circ_0001550, and hsa_circ_0001681) in both of paired CRC tissues and CRC cell lines. Then, a circRNA/microRNA/messenger RNA regulatory network was established and the Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed these four circRNAs participated in various biological processed including apoptotic process and multiple metabolic processes. Moreover, based on the regulatory network, three bioactive compounds (pergolide, pivampicillin, and methylergometrine) for the treatment of CRC were also found. In conclusion, this study improved our understanding of circRNAs and may also facilitate the finding of promising targets and biomarkers in CRC.     

Circular RNA hsa_circ_0006401 promotes proliferation and metastasis in colorectal carcinoma

Dysregulation of circular RNA (circRNA) expression is involved in the progression of cancer. Here, we aimed to study the potential function of hsa_circ_0006401 in colorectal cancer (CRC). CircRNA hsa_circ_0006401 expression levels in CRC and adjacent nontumor tissues were analyzed by real-time quantitative PCR (qRT-PCR) and circRNA in situ hybridization (RNA-ISH). Then, CRC cell proliferation was assessed by cell counting. Wound-healing and transwell assays were utilized to detect the effect of hsa_circ_0006401 on CRC migration. A circRNA-ORF construct was created, and a specific antibody against the splice junction of hsa_circ_0006401 was prepared. Finally, the proteins directly binding to hsa_circ_0006401 peptides were identified by immunoprecipitation combined with mass spectrometry. In our study, we found hsa_circ_0006401 was closely related to CRC metastasis and exhibited upregulated expression in metastatic CRC tissue samples. Proliferation and migration were inhibited in vitro when hsa_circ_0006401 expression was silenced. Downregulation of hsa_circ_0006401 expression decreased CRC proliferation and liver metastasis in vivo. A 198-aa peptide was encoded by sequences of the splice junction absent from col6a3. Hsa_circ_0006401 promoted CRC proliferation and migration by encoding the hsa_circ_0006401 peptide. Hsa_circ_0006401 peptides decreased the mRNA and protein level of the host gene col6a3 by promoting col6a3 mRNA stabilation. In conclusion, our study revealed that circRNAs generated from col6a3 that contain an open-reading frame (ORF) encode a novel 198-aa functional peptide and hsa_circ_0006401 peptides promote stability of the host gene col6a3 mRNA to promote CRC proliferation and metastasis.Subject terms: Cancer, Cell biology     

Systematic identification of lncRNAs and circRNAs-associated ceRNA networks in human lumbar disc degeneration

ABSTRACT

Lumbar disc degeneration (LDD) is a common cause of low back and neck pain. The molecular mechanisms underlying LDD, however, are unclear. Noncoding RNAs have been reported to participate in human diseases. We investigated a series of public datasets (GSE67566, GSE56081 and GSE63492) and identified 568 mRNAs, 55 microRNAs (miRNAs), 765 long noncoding RNAs (lncRNAs), and 586 circular RNAs (circRNAs) that were expressed differently in LDD than in normal discs. We constructed lncRNAs and circRNAs regulated competing endogenous RNAs (ceRNA) networks in LDD. Four lncRNAs, DANCR, CASK-AS1, SCARNA2, and LINC00638), and three circRNAs, hsa_circ_0005139, hsa_circ_0037858, and hsa_circ_0087890, were identified as key regulators of LDD progression. We found that hsa-miR-486-5p regulated the crosstalk among circRNA hsa_circ_0000189, lncRNA DANCR and 6 mRNAs, PYCR2, TOB1, ARHGAP5, RBPJ, CD247, SLC34A1. Gene ontology (GO) analysis demonstrated that these differently expressed lncRNAs and circRNAs were involved in cellular component organization or biogenesis, gene expression and negative regulation of metabolic processes. Our findings provide useful information for exploring new mechanisms for LDD and candidates for therapeutic targets.     

Comprehensive characterization of circular RNAs in osteosarcoma cell lines

Circular RNA (circRNA) is a looped noncoding RNA with a stable structure and tissue-specific expression and widely reported to regulate cancer initiation and progression. However, the circRNA expression patterns and their roles in osteosarcoma initiation and progression are still poorly understood. In this study, we characterized the landscape of circRNAs in osteosarcoma (OS) cell lines, and calculated the epithelial-mesenchymal transition (EMT) scores for OS cell lines. The differential expression analysis revealed that the EMT-related genes were significantly upregulated in the OS cell lines with higher metastatic potentials, and some inflammation-related pathways and pathways involved in cell-cell communications were enriched by these upregulated genes. Furthermore, we constructed a circRNA-based competing endogenous RNA (ceRNA) network, which consisted of 5 circRNAs, 17 miRNAs, and 73 mRNAs. Particularly, hsa_circ_0085360, which had the highest correlation with TRPS1, were characterized by some cancer-related pathways, and TRPS1 and its target gene FGFR3 were closely associated with both event-free survival and overall survival of OS, indicating that hsa_circ_0085360 might have the potential to predict the OS prognosis. In summary, we profiled the circRNA expression patterns in OS, predicted their functionality, and explored the underlying mechanism and prognostic values, which might provide some evidences for OS-related circRNA researches.     

Profiling and integrated analysis of differentially expressed circRNAs in cervical cancer

Circular RNAs (circRNAs) are a new type of regulatory RNAs, which have been identified to play critical role in various tumors. However, the profiles and roles of circRNAs in cervical cancer (CCa) have not been fully understood and need to be further explored. In the present study, we performed circRNA array and mRNA-sequencing (mRNA-Seq) to profile the differentially expressed circRNAs and mRNAs in CCa tissues. A total of 397 differentially expressed circRNAs and 2138 differentially expressed mRNAs were detected, respectively. Subsequently, a circRNA-miRNA-mRNA regulatory network was constructed and indicated that hsa_circ_0026377 was downregulated in CCa. Overexpression of hsa_circ_0026377 inhibited HeLa and SiHa cells proliferation, migration and invasion. Collectively, this study provided new insights into the circRNA profiles in CCa and suggested that hsa_circ_0026377 might play important roles in CCa development.     

Identification of downregulated circRNAs from tissue and plasma of patients with gastric cancer and construction of a circRNA-miRNA-mRNA network

The connection between circular RNAs (circRNAs) and gastric cancer has been reported widely in recent years. However, previous studies have focused mainly on circRNAs from gastric cancer tissue. The objectives of the present study were to detect dysregulated circRNAs from both tissue and plasma of patients with gastric cancer and to explore their potential roles in the pathogenesis of gastric cancer. Expression profiles of circRNAs were obtained from the Gene Expression Omnibus (GEO) and analyzed using the GEO2R tool to identify differential expressed circRNAs. The significance threshold was set as |log2 (fold change)| > 2 and adjusted P < .05. The microRNA (miRNA) binding sites of the differentially expressed circRNAs were predicted using the Circular RNA Interactome web tool. TargetScan and the miRNet database were used to predict the miRNA target genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed using Database for Annotation Visualization and Integrated Discovery. Hub genes were identified and a network was constructed with Cytoscape. The overall survival rates for the selected miRNAs and messenger RNAs were evaluated by Kaplan-Meier Plotter. A total of three downregulated circRNAs (hsa_circ_0001190, hsa_circ_0036287, and hsa_circ_0048607) were identified in this study. Six miRNAs and eight hub genes met the significance criteria and were selected for further analysis. A circRNA-miRNA-hub gene network was constructed based on three circRNAs, six miRNAs, and eight hub genes. Evaluation of overall survival rates for the hub genes showed that low expression levels of GADD45A, PPP1CB, PJA2, and KLF2 were associated with poor overall survival. This study identified potential novel plasma circRNA biomarkers and provides insights into the underlying mechanisms of gastric cancer pathogenesis.     

Analysis of co-expression networks for circular RNAs and mRNAs reveals that circular RNAs hsa_circ_0047905, hsa_circ_0138960 and has-circRNA7690-15 are candidate oncogenes in gastric cancer

Accumulating evidence has suggested that circular RNAs (circRNAs) play important roles in oncogenesis and tumor progression. However, our knowledge of circRNAs in gastric cancer (GC) remains limited. To investigate circRNAs involved in GC oncogenesis, we examined differentially-expressed circRNAs and mRNAs in GC tissues and paired noncancerous mucosa tissues using circRNA and mRNA microarrays. Next, we built gene co-expression networks according to the degree of correlation to predict the critical circRNAs in GC. Through bioinformatics analysis, we observed three newly identified circRNAs that are substantially upregulated in GC: hsa_circ_0047905, hsa_circ_0138960 and has-circRNA7690-15. Additionally, hsa_circ_0047905 and hsa_circ_0138960 positively correlated with their parental gene mRNA. Knockdown of hsa_circ_0047905, hsa_circ_0138960 and has-circRNA7690-15 in GC cells, resulted in downregulation of parental gene expression. Functional assays suggested that inhibition of these three circular RNAs suppresses GC cell proliferation and invasion in vitro. Those findings suggest that hsa_circ_0047905, hsa_circ_0138960 and has-circRNA7690-15 might act as tumor promoters in the pathogenesis of gastric cancer.     

Circular RNA expression profiles and features in human tissues: a study using RNA-seq data

Background

Circular RNA (circRNA) is one type of noncoding RNA that forms a covalently closed continuous loop. Similar to long noncoding RNA (lncRNA), circRNA can act as microRNA (miRNA) ‘sponges’ to regulate gene expression, and its abnormal expression is related to diseases such as atherosclerosis, nervous system disorders and cancer. So far, there have been no systematic studies on circRNA abundance and expression profiles in human adult and fetal tissues.

Results

We explored circRNA expression profiles using RNA-seq data for six adult and fetal normal tissues (colon, heart, kidney, liver, lung, and stomach) and four gland normal tissues (adrenal gland, mammary gland, pancreas, and thyroid gland). A total of 8120, 25,933 and 14,433 circRNAs were detected by at least two supporting junction reads in adult, fetal and gland tissues, respectively. Among them, 3092, 14,241 and 6879 circRNAs were novel when compared to the published results. In each adult tissue type, we found at least 1000 circRNAs, among which 36.97–50.04% were tissue-specific. We reported 33 circRNAs that were ubiquitously expressed in all the adult tissues we examined. To further explore the potential “housekeeping” function of these circRNAs, we constructed a circRNA-miRNA-mRNA regulatory network containing 17 circRNAs, 22 miRNAs and 90 mRNAs. Furthermore, we found that both the abundance and the relative expression level of circRNAs were higher in fetal tissue than adult tissue. The number of circRNAs in gland tissues, especially in mammary gland (9665 circRNA candidates), was higher than that of other adult tissues (1160–3777).

Conclusions

We systematically investigated circRNA expression in a variety of human adult and fetal tissues. Our observation of different expression level of circRNAs in adult and fetal tissues suggested that circRNAs might play their role in a tissue-specific and development-specific fashion. Analysis of circRNA-miRNA-mRNA network provided potential targets of circRNAs. High expression level of circRNAs in mammary gland might be attributed to the rich innervation.

     

Genome-wide circular RNA (circRNA) and mRNA profiling identify a circMET-miR-410-3p regulatory motif for cell growth in colorectal cancer

Circular RNA (circRNA) is a non-coding RNA molecule that lacks polyadenylated tails and is highly stable, abundant, and conserved in human cells. CircRNAs can serve as a competing endogenous RNA (ceRNA) to sponge microRNAs (miRNA) and block their effects on target mRNA expression. CircRNAs also have possible relevance to cancer and therefore may be considered as ideal biomarkers for monitoring cancer progression. Of the about 300,000 predicted human circRNAs, only a few have validated biological functions related to cancer. To better understand the ceRNA role of circRNAs in colorectal cancer (CRC), we performed genome-wide circRNA-based RNA-sequencing (RNA-Seq) on nine CRC tumor samples and their paired histologically normal adjacent tissue samples. By profiling the mRNA expression in the same patients, we further explored the expression correlation between circRNAs and mRNAs generated from the same parental gene. Focusing on the concordant differential expression between circRNAs and mRNAs, we substantially reduced the regulatory noise. In total, we identified 694 circRNA-mRNA pairs that were consistently up or downregulated between tumor and normal tissues. These 694 circRNA-mRNA pairs are from 182 protein-coding genes associated with hormone responses and chemotaxis. Of these 182 genes, 43 are downstream targets of three highly conserved miRNAs (miR-410-3p, miR-135a, and miR-30a). Interestingly, these 43 genes are highly mutated in another cohort from eight independent CRC studies, which have significant effects on patient survival time. Focusing on miR-410-3p and its target oncogene MET, we experimentally validated the ceRNA regulatory motif of circMET. Notably, circMET is substantially upregulated in CRC cell lines and could promote cell proliferation and growth. By confirming the regulatory relationship between miR-410-3p and circMET, we propose a new mechanism for the observed sustained activation of MET in CRC. In conclusion, our work identifies a novel regulatory role of circMET and provides a potential diagnostic biomarker for CRC.     

Hsa_circ_0136666 promotes the proliferation and invasion of colorectal cancer through miR-136/SH2B1 axis

Colorectal cancer (CRC) is one of the leading causes of cancer-related death worldwide. Currently, an increasing evidence showed that circular RNAs (circRNAs) play important roles in tumor progression. However, the effects and underlying mechanisms of circRNAs in CRC progression remain unclear. In the present study, through circRNA high-throughput sequencing and quantitative real-time polymerase chain reaction, we identified that hsa_circ_0136666 was significantly overexpressed in CRC tissues and cell lines. High hsa_circ_0136666 expression was associated with poor overall survival of patients with CRC. In vitro function assays showed that hsa_circ_0136666 inhibition suppressed CRC cell proliferation, migration, invasion, and arrested CRC cells in the G0/G1 phase. Furthermore, we showed that hsa_circ_0136666 inhibition reduced CRC cell growth in vivo. Mechanistically, we revealed that hsa_circ_0136666 could increase SH2B1 expression via competitively binding miR-136 in CRC cells. In addition, SH2B1 overexpression could reverse the effects of hsa_circ_0136666 inhibition on CRC cell progression. In conclusion, our data suggested that hsa_circ_0136666 could promote CRC cell progression via the miR-136/SH2B1 axis, elucidating a novel approach to improve the effectiveness of CRC treatment.     

Circular RNA expression profile of knee condyle in osteoarthritis by illumina HiSeq platform

Osteoarthritis (OA) is the utmost commonly arising joint disease. Knee condyles play an essential role during OA progression. Circular RNA (or circRNA) is a novel kind of RNA, which, unlike the well-known linear RNA, plays an important regulatory role in OA on the basis of a previous research. In our study, expression of circRNAs in OA knee condyle was measured by illumine sequencing platform. A total of 197 differentially expressed circRNAs, such as hg38_circ_0007474 and hg38_circ_0000118 were identified, and 21 target miRNAs, 2466 source genes and 166 394 circRNA-miRNA-mRNA pairs were predicted. Further analysis was applied on three OA-related circRNAs (hsa_circ_0045714, hsa_circ_0002485, and hsa_circ_0005567). The results were partly verified by previous studies. Further biological research is needed to unfold the possible pathway and therapeutic target of OA.     

Genome‐Wide Analysis of circular RNAs and validation of hsa_circ_0006719 as a potential novel diagnostic biomarker in congenital scoliosis patients

Congenital scoliosis (CS) is a form of spinal curvature resulting from anomalous development of vertebrae. Recent studies demonstrated that circRNAs could serve as potential biomarkers of disease diagnosis. Genome‐wide circRNAs expression in seven CS patients and three healthy controls was initially detected. Bioinformatics analysis was conducted to explore the potential pathological pathway of CS. Quantitative PCR (qPCR) was performed to validate the selected circRNAs in the replication cohort with 32 CS patients and 30 healthy controls. Logistic regression controlling for gender was conducted to compare the expression difference. Receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic value. Twenty‐two differentially expressed circRNAs were filtered from genome‐wide circRNA sequencing. Seven circRNAs were validated by qPCR. Only hsa_circ_0006719 was confirmed to have a higher expression level in the CS group than the healthy control group (P = 0.036). Receiver operating characteristic curve also suggested that hsa_circ_0006719 had significant diagnostic value for CS (AUC = 0.739, P = 0.001). We described the first study of circRNAs in CS and validated hsa_circ_0006719 as a potential novel diagnostic biomarker of CS.     

Profiling and integrated analysis of differentially expressed circRNAs as novel biomarkers for breast cancer

Breast cancer (BC) is a globally common cancer with the highest and increasing morbidity and mortality among females. Novel biomarkers are warranted to be discovered for the early detection, treatment, and prognosis of BC. In this study, we investigated the profiles of differentially expressed (DE) circular RNAs (circRNAs) by competing endogenous RNAs (ceRNA) microarray to construct a genome-wide circRNA profile. Then, we performed Gene Ontology (GO) analysis and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis of the host genes (HGs) of circRNAs. A total of 4,370 DE circRNAs were detected and GO and KEGG analysis showed that they were significantly associated with cell cycle, DNA replication, BC, and familial BC. We validated the differential circRNAs and relevant HGs through quantitative real-time polymerase chain reaction and constructed a putative circRNA–microRNA–messenger RNA regulatory network. Eight circRNAs, including hsa_circ_0069094, hsa_circ_0062558, hsa_circ_0074026, hsa_circ_0079876, hsa_circ_0017536, hsa_circ_0023302, hsa_circ_0017650, and hsa_circ_0017545, were validated significantly DE in BC tissue and associated with TNM staging, lymph node infiltration, and Ki67. Hsa_circ_0069094, hsa_circ_0079876, hsa_circ_0017650, and hsa_circ_0017526 were upregulated in plasma. This study revealed the general expression characteristics of specific DE circRNAs in BC and hsa_circ_0069094, hsa_circ_0079876, hsa_circ_0017650, and hsa_circ_0017526 might be promising candidate targets.     

环状RNA：胃癌诊治的新星

环状RNA (circular RNA,circRNA)是一类闭合环状结构的RNA分子，广泛分布于各种组织中，它比线性RNA更稳定。circRNA分为外显子circRNA、外显子-内含子circRNA和内含子circRNA等3类。circRNA的主要功能为充当微RNA海绵、与RNA结合蛋白结合、翻译成蛋白质和调节转录等。近年来，大量研究表明，circRNA的异常表达在胃癌发生发展过程中起着至关重要的作用。circPTPN22、hsa＿circ＿0001772、circCYFIP2、hsa＿circ＿0017639和circPIP5K1A等的上调以及hsa＿circ＿002059、hsa＿circ＿0000190和circMTO1等的下调与胃癌的增殖和转移密切相关；而hsa＿circ＿0001313等影响胃癌细胞的顺铂耐药性。组织、血浆及外泌体中circPTPN22、hsa＿circ＿102958、hsa＿circ＿0141633、hsa＿circ＿0065149和hsa＿circ＿0026344等是胃癌新型诊断标志物；而hsa＿circ＿0005529、circ-RanGAP1、cir...     

Identification and characterization of tumorigenic circular RNAs in cervical cancer

Circular RNAs (circRNAs) are a distinctive family of ncRNAs, and they function as key regulators in the initiation, development and progression of various diseases. However, the regulatory roles of circRNAs in the tumorigenesis of cervical cancer (CC) have not been fully understood. In this study, we identified a set of circRNAs in CC and paired normal tissues, using RNA sequencing data, and found that the cancer and normal tissues could be told apart by those differentially expressed (DE) circRNAs, indicating that circRNA expression profiles in CC were significantly different from those in the normal tissues. Meanwhile, the upregulated genes in CC were enriched in inflammation-related pathways, and the correlation analysis between the DE circRNAs and genes revealed that the abundance of DE circRNAs was positively related to the expression of their host genes. However, the expression of TGFBR2 and KDM4C were found to exhibit a negative correlation with their corresponding circRNAs. Furthermore, we also predicted the interactions between circRNAs and proteins, and constructed a competing endogenous RNA (ceRNA) network. Specifically, hsa_circ_0001495 was predicted to interact with ESRP2, and acted as a sponge by competing for miRNAs with TBL1XR1. Functionally, hsa_circ_0001495 was predicted to regulate epithelial cell proliferation and NOTCH signaling via ESRP2 and TBL1XR1, respectively. We also evaluated the prognostic values of downstream target genes of selected circRNAs, using clinical records of CC patients. In summary, the present study provided some regulatory circRNAs involved in CC tumorigenesis based on bioinformatics approaches, which brought strong evidences for researchers to further explore their biological and clinical values.     

An evaluation of gastric adenocarcinoma-associated CircRNAs based on microarray meta-analysis and ceRNA networks

Gastric cancer is the fourth leading cause of cancer-related mortality and one of the most commonly diagnosed malignancies worldwide. Gastric adenocarcinoma (GAC) accounts for the majority of gastric cancer cases. Circular RNAs (circRNAs) have been shown to be associated with carcinogenesis and cancer progression. This research aims to investigate GAC-associated circRNAs and the underlying mechanisms of circRNA-miRNA-mRNA networks in the development and progression of GAC. Differentially expressed miRNAs and mRNAs (DEMs and DEGs) were identified in Gene Expression Omnibus (GEO) microarray datasets using the R package Limma. A microarray meta-analysis was performed to identify potential GAC-associated circRNAs with high statistical power, resulting in 13 up-regulated and 19 down-regulated circRNAs. CircRNA-miRNA-mRNA networks were constructed by combining predicted and experimentally validated databases and considering differentially expressed miRNAs and mRNAs. The constructed ceRNA networks revealed the potential regulatory effect of hsa_circ_0002019 and hsa_circ_0074736 on key survival-related genes. The expression levels of these two circRNAs were measured in plasma samples from GAC patients and healthy controls using SYBR Green-based real-time PCR. Axon guidance, cellular senescence, AGE-RAGE signaling pathway in diabetic complications, and AMPK signaling pathway were among the major significant (P-value <0.05) enriched pathways of "main mRNAs" in the constructed ceRNA networks. In conclusion, we identified strongly correlated circRNAs and their likely mechanisms of action in GAC, which may improve the knowledge of regulatory networks underlying GAC formation and contribute to developing better strategies for early diagnosis, prognosis, and treatment.