OSCC is one of the most common malignant tumors in the head and neck cancer. Its incidence is increasing year by year in the world, and the population of OSCC tends to be younger, which poses a serious threat to human health and life. Therefore, early detection and treatment of OSCC should not be delayed. So we build a dynamic animal model of OSCC in different stages, which can make up for the difficulty of obtaining samples in different clinical stages. In order to better prognosis treatment, it is important to elucidate the pathogenesis of OSCC. Non-coding RNA as a research hotspot in oncology, plays an important role in OSCC.
At present, the non-coding RNAs studied mainly include miRNA, lncRNA and circRNA [22]. Among them, miRNA and lncRNA are more mature, and circRNA is a hotspot in recent years. Compared with linear miRNA and mRNA, circRNA as a disease diagnostic marker has the unique advantages of high content and stable structure. Many studies on miRNAs and lncRNAs have been carried out and have demonstrated that these RNAs are related to oral squamous cell carcinoma proliferation, migration and invasion. Studies have shown that the miR-504 inhibiting cell proliferation, migration and invasion by regulating the expression of CDK6 in the process of OSCC [23]. Chen et al found that miRNA-10a promotes cancer cell proliferation in oral squamous cell carcinoma by regulating GLUT1 and promoting glucose metabolism [24]. The literature reports that the lncRNA LUCAT1 might promote OSCC migration and invasion by regulating PCNA [25]. Moreover, some factors also regulate the function of ncRNA, for instance, Chen et al has pointed out NKILA has been shown to bind to the nuclear factor-kB (NF-kB)/IkB complex, which masks the phosphorylation sites of IkB, thus preventing its degradation [26]. It is well known that circRNA is stably present in vivo and vitro due to its covalently closed-loop structure. Consequently, using circRNA as a therapy target would better than miRNA and lncRNA.
Based on the above considerations, we established and validated oral squamous cell carcinoma animal model using Chinese hamster. Then we used RNA-seq technology to analyze the expression profiles of circRNA in tissues of OSCC at three different pathological stages and normal control tissues, and screen significantly different circRNAs in OSCC tissues. Our results pointed out that the expression of circRNA is abundant in twelve samples, including 3485 existing circRNA molecules and 229 novel circRNA molecules, 38 circRNAs were up-regulated while 51 were down-regulated in model group compared with control group. We selected four up-regulated and four down-regulated circRNAs for qRT-PCR validated that accordance with the high-throughput sequencing results. To the best of our knowledge, most of circRNAs were differentially expressed, indicating that these differentially expresses circRNAs might be implicated in the incidence and progression of OSCC. GO function and KEGG pathway analysis of the host genes of the significantly different circRNAs were used to predict related functions through bioinformatics. Based on the correlation between circRNA and its host gene and its known function, GO annotation analysis results show that GO term is mainly enriched in epithelial cell proliferation, actin-based cell projection and protein domain specific binding; KEGG further analysis shows that KEGG term is mainly enriched in Wnt, PI3K-Akt, MAPK, TGF-β, p53 and mTOR signaling pathways, revealing that the host genes of circRNAs may participate in the regulation of OSCC development through related molecular functions or signaling pathways. Besides important function and pathway, more and more studies have revealed that the mechanism of circRNAs sponge miRNAs [27, 28]. Finally, 35 abnormally expression circRNAs were selected to predict the circRNA-miRNA regulatory network from the perspective of ceRNA, which the potential biological functions were analyzed from the perspective of possible molecular mechanisms and revealed the important regulatory functions of circRNA in OSCC pathology. The above bioinformatics prediction analysis shows that these significantly different circRNAs have important regulatory roles in the pathology of OSCC, and are involved in promoting or inhibiting the occurrence and development of OSCC.
Due to the conservation and tissue specificity of circRNA, we screened up-regulated expressed cgr_circ_MAN2A1 in tissues by high-throughput sequencing. Therefore, we combined with bioinformatics databases (circBase, circBank and UCSC databases) and bioinformatics analysis to screen the gene of MAN2A1 related circRNA, it was found that hsa_circ_0127523 is obtained by cleaving the pre-mRNA of MAN2A1 and the full-length sequence of hsa_circ_0127523 was gotten. By detecting the expression level of the gene hsa_circ_0127523, it was found that the expression level of hsa_circ_0127523 in oral squamous cell carcinoma cell lines was significantly higher than that of normal cells. It is necessary to further study whether hsa_circ_0127523 plays an important regulatory role in OSCC. According to Ouyang et al research reports that circRNA_0109291 is highly expressed in oral squamous cell carcinoma and the authors tested its functions and mechanisms by design to interfere with the circRNA_0109291 in OSCC [29]. Then, the behavioral function of hsa_circ_0127523 in CAL27 and Tca8113 cell lines shows that silencing hsa_circ_0127523 inhibited the proliferation, migration and invasion of OSCC cells in vitro.
It is well known that circRNAs play crucial roles in cancer by functioning as miRNA sponges, and affect the stability of target mRNA and regulate mRNA expression level [30, 31]. The latest literatures report the mechanism of circRNA to miRNA is a competitive endogenous mechanism, also called sponge mechanism [32], and the expression of circRNA and mRNA is positively correlated, the expression of circRNA and mRNA is negatively correlated with the expression of miRNA, which formed “circRNAs-miRNAs-mRNAs” target regulation network [33–35]. Studies have shown that hsa_circRNA_100533 act as a competitive endogenous RNA (ceRNA) and regulates the expression of GNAS by sponging has-miR-933, thereby modulating the proliferation, migration and apoptosis of OSCC [36]. By using a circRNA microarray, Gao et al demonstrated that circ-PKD2 as a ceRNA regulates APC2 expression by sponge miR-204-3p and participates in OSCC cell proliferation, migration and invasion, induced apoptosis and cell cycle arrest [37].
In our next study, the interaction network of hsa_circ_0127523 showed that it was most likely to regulate the expression of miR-515-5p according to the bioinformatics analysis. Further interaction analysis showed that the target gene of miR-515-5p most likely to affect the TRIP13 gene regulation. Therefore, we detected the expression of miR-515-5p and TRIP13 at the gene level after hsa_circ_0127523 was silenced, and found that the expression level of miR-515-5p decreased and TRIP13 increased after silencing hsa_circ_0127523. Moreover, Zhang et al found that TRIP13 as a direct target of miR-515-5p by binding to it’s 3’-UTR, which abolish the role of miR-515-5p in suppressing cell proliferation, migration and invasion that the miR-515-5p-TRIP13 axis may be a brightening prognostic marker and therapeutic target for prostate cancer [21]. TRIP13 (Thyriod hormone receptor interacting protein 13) is a nuclear factor, which is a protein-coding gene from the perspective of gene type. It is mainly involved in meiosis, cell proliferation and metastasis and is expressed in normal and cancerous tissues. Recent studies have found that it’s related to tumorigenesis, development and prognosis [38, 39]. Kurita et al [40] found through research that TRIP13 is highly expressed in cancer tissues of colorectal cancer patients and the deletion of the TRIP13 gene can inhibit cancer cell proliferation, migration and invasion, which may become a potential target for colorectal cancer treatment. In addition, TRIP13 gene is abnormally highly expressed in head and neck malignancies, multiple myeloma, chronic lymphocytic leukemia and lung adenocarcinoma [41–43], which may be related to the occurrence and development of malignant tumors. There are no reports about the study of TRIP13 in OSCC. Therefore, this study hopes to provide valuable reference for clinical treatment from the perspective of TRIP13. Our study forecasted the sponge mechanisms of hsa_circ_0127523/miR-515-5p/TRIP13 axis in OSCC, which verified the differential expression level that are coincidence with competitive endogenous mechanism.
In addition, the expression level of Notch1 and p21 were abnormal when inhibited the expression level of hsa_circ_0127523 in OSCC cells, current studies have suggested that hsa_circ_0127523 may regulate the occurrence and development of OSCC through Notch pathway. Role of TRIP13 has been reported in the literature, and related research validation has also suggested that TRIP13 gene affects the expression levels of Notch1 and p21 in the Notch pathway in ovarian cancer [21, 44]. Analysis by Yoshida et al reported that Notch1 plays a vital role in the development and progress of OSCC. Based on deep sequencing of the cancer genome, research has also suggested the anti-cancer effect of Notch1. Notch plays an indispensable role in OSCC because Notch is considered to play an oncogene role in a subset of OSCC and also has a suppressive effect [45]. Lv et al has found that miR-495 restricts the proliferation and invasion of OSCC cells by directly targeting Notch1. Notch1 knockdown shows an inhibitory effect similar to that induced by miR-495 overexpression in OSCC cells [46]. Cyclin-dependent kinase (CDK) and cyclin-dependent kinase inhibitor (CDKI) can directly control the cell cycle. p21 is a negative regulator of cyclin-dependent kinases and is a negative checkpoint regulator of the cell cycle in this function [47]. Previous studies have shown that p21 is over expressed in lung cancer and gastrointestinal cancer cells [48, 49].The results of Zhang's study show that the 5-year survival rate of OSCC patients with high p21 expression is significantly reduced. Kaplan-Meier analysis and Cox proportional hazard model showed that lymph node metastasis and p21 expression were independent prognostic factors of OSCC [50], which are consistent with our findings that silenced the hsa_circ_0127523, the TRIP13 expression level was decreases and the Notch1 decreases in the Notch pathway.
Increasing studies suggest that circRNA has a vital part in a variety of tumors as a tumor regulator. As reported by Su that hsa_circ_0055538 could regulate the malignant biological behavior of OSCC through the p53/Bcl-2/caspase signaling pathway [51]. And it has been reported that circDOCK1 suppresses cell apoptosis via inhibition of miR‑196a‑5p by targeting BIRC3 in OSCC [19]. Chen et al found that circRNA_100290 was bound to miR-29 family members directly and co-expressed with CDK6 influenced proliferation in oral cancer cell [52]. In non-small cell lung cancer, Huang et al demonstrated that hsa_circ_0001946 may serve as a novel biomarker for the diagnosis and prediction of platinum-based chemosensitivity in patients with NSCLC [53]. Sang et al has identificated that hsa_circ_0025202 served an anti-oncogenic role in HR-positive breast cancer, could suppress tumor growth and enhance tamoxifen efficacy, and it could be exploited as a novel marker for breast cancer [54]. Based on these, we propose a view that circRNA may participate in development and prognosis of OSCC and our study showed that hsa_circ_0127523/miR-515-5p/TRIP13 by modulated Notch pathway may play an important role in OSCC cancer development. However, hsa_circ_0127523 has not been reported in the any disease, so it is concluded that hsa_circ_0127523 may be a potential biomarker of OSCC. In the next study, we will further elucidate the regulatory mechanism of hsa_circ_0127523/miR-515-5p/TRIP13 in OSCC.
In conclusion, our study provided evidence that oral circRNAs are differentially expressed in response to OSCC. In the expression signatures, we uncovered the circRNAs expression profile in three dynamic stages by comparing it with normal tissue. Additionally, eight dysregulated circRNAs expression level in tissue were verified by qRT-PCR, which is consistent with the sequencing result. Additionally, we confirmed hsa_circ_0127523 was interfered significantly affected the proliferation, migration and invasion abilities of OSCC cells. More important, we predicted and confirmed hsa_circ_0127523 inhibit miR-515-5p expression and increase TRIP13 expression in the OSCC cells which provide potential axis for the next step diagnosis of OSCC and novel genetic insights, and the specific mechanism still needs further research. Hsa_circ_0127523 may be a potential marker and therapeutic target for the diagnosis of OSCC.