==== Front Int J OncolInt. J. OncolIJOInternational Journal of Oncology1019-64391791-2423D.A. Spandidos 3006684710.3892/ijo.2018.4499ijo-53-04-1529ArticlesCompeting endogenous RNA analysis reveals the regulatory potency of circRNA_036186 in HNSCC Wang Wen-Long 12Yang Zhi 12Zhang Yi-Juan 3Lu Ping 12Ni You-Kang 12Sun Chang-Fu 12Liu Fa-Yu 12 1 Departments of Oromaxillofacial-Head and Neck Surgery 2 Departments of Oral and Maxillofacial Surgery, School of Stomatology, China Medical University, Shenyang, Liaoning 110002 3 Department of Anesthesiology, The First Affiliated Hospital, China Medical University, Shenyang, Liaoning 110002, P.R. ChinaCorrespondence to: Dr Fa-Yu Liu or Dr Chang-Fu Sun, Department of Oromaxillofacial-Head and Neck Surgery, School of Stomatology, China Medical University, 117 Nanjing Bei Jie, Heping, Shenyang, Liaoning 110002, P.R. China, E-mail: lfyhjk@126.com, E-mail: changfusun@hotmail.com10 2018 24 7 2018 24 7 2018 53 4 1529 1543 26 1 2018 26 6 2018 Copyright: © Wang et al.2018This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.This study aimed to characterize circular RNA (circRNA) expression profiles and biological functions in head and neck squamous cell carcinoma (HNSCC). Differentially expressed circRNAs were screened using an Arraystar Human CircRNA Array and verified by reverse transcription-quantitative polymerase chain reaction. Multiple bioinformatics methods and a hypergeometric test were employed to predict the interactions between RNAs and the functional circRNA-microRNA (miRNA)-mRNA axes in HNSCC. As a result, 287 circRNAs and 1,053 mRNAs were determined to be differentially expressed in HNSCC compared with the adjacent tissue. In addition, the expression levels of circRNA_036186 and tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, ζ polypeptide (14-3-3ζ) were identified to be significantly different. A competing endogenous RNA (ceRNA) network was constructed, consisting of 5 circRNAs, 385 miRNAs and 96 mRNAs. Furthermore, we predicted that miR-193b-3p exerts a significant effect on 14-3-3ζ, and was significantly associated with the Hippo signaling pathway in HNSCC. On the whole, these findings suggest that circRNA_036186 likely regulates 14-3-3ζ expression by functioning as a ceRNA in the development and progression of HNSCC. circular RNAmicroRNAcompeting endogenous RNAmicroarrayhead and neck squamous cell carcinoma ==== Body Introduction Head and neck squamous cell carcinoma (HNSCC) is the sixth most common type of cancer worldwide, with ~600,000 new cases diagnosed each year. Treatment protocols include surgery, chemotherapy and radiation, all of which may result in major physical and psychological effects on the patient (1). Despite continuing research and advances in treatment, the clinical outcomes and overall survival rates of patients with HNSCC have not improved significantly over the past several decades, with an overall 5-year survival rate of 60% (2,3). A high rate of metastasis is a characteristic of HNSCC, which highlights the urgent need for the identification of novel biomarkers to more accurately predict recurrence and metastasis in patients with HNSCC. MicroRNAs (miRNAs or miRs) are endogenous short non-coding RNAs (ncRNAs) that bind to the complementary seed sequence on the 3′-untranslated region (UTR) of specific target protein-coding genes, known as miRNA response elements (MREs), thereby resulting in the inhibition of target gene expression (4). Increasing evidence suggests that miRNA dysregulation is associated with the tumorigenesis and progression of various cancer types (5-8). It has been demonstrated that certain protein-coding genes and their pseudogene sequences contain the same conserved MREs in their 3′UTRs, and their respective expression levels may be regulated through competing for the same miRNAs through these MREs (9). A hypothesis regarding these RNA transcripts, termed competing endogenous RNAs (ceRNAs), has been proposed (10). According to this hypothesis, MREs can be seen as an 'RNA language' through which transcripts can communicate with each other to regulate their respective expression levels. Any RNA transcript with MREs may function as a ceRNA and compete with other RNA transcripts with similar MREs for the same miRNAs. Circular RNAs (circRNAs) are widespread and diverse endogenous ncRNAs with covalently closed circular structures (11,12). There are three main types of circRNAs based on their biogenesis patterns, including exonic circRNAs (13), intronic circRNAs (14,15) and exon-intron circRNAs (15-17), which are transcribed from pre-mRNA sequences by RNA polymerase (Pol) II (18). circRNAs are predominantly generated by reverse splicing between the specific and conservative sequences in upstream and downstream regions, respectively, a process known as back-splicing (19,20). Of the three types, exonic circRNAs account for the largest proportion, and are abundant in the cytoplasm (21-23). There are multiple conserved MREs on circRNAs, known as miRNA sponges, that sequester miRNAs from binding to targets on mRNAs, thereby enhancing the expression of certain genes (24,25). However, studies have demonstrated that only a few circRNAs have multiple MREs for specific miRNAs, most of which contain only one or two MREs (26,27). Intronic circRNAs are mainly present in the nucleus and have almost no enrichment in MREs. They largely accumulate around the transcription initiation sites of their parental genes through interaction with elongation RNA Pol II machinery, serving a cis-regulatory role on their parental coding genes (14). Exon-intron circRNAs are abundant in the nucleus, and interact with small nuclear ribonucleoprotein (U1 snRNP) and Pol II at the promoter regions of genes to promote the transcription of their parental genes (16). At present, studies on circRNAs are mainly focused on exonic circRNAs. Their interactions with disease-related miRNAs indicate that circular RNAs are important factors in disease regulation (25,28). For instance, the circRNA ciRS-7 contains multiple and conservative miRNA-7 binding sites, thus acting as an efficient miRNA sponge to adsorb miRNA-7, and thus inhibiting its function (29). Increasing evidence has demonstrated that ciRS-7 plays a crucial role in the tumorigenesis and progression of cancers through its effect on the level of miRNA-7, thus mediating gene expression via ceRNA mechanisms (30). Given the widespread involvement of miR-7 in various cancer-related pathways and its potential regulatory functions in Parkinson's and Alzheimer's disease, ciRS-7 may play a key role in other neurological disorders and tumorigenesis (24). Exploring the mechanisms underlying the interactions between these ceRNAs may add a new dimension to our understanding of the mechanisms of cancer, and may potentially provide a novel approach to cancer treatment. The molecular mechanisms of circRNAs in the carcinogenesis and progression of HNSCC are largely uncharacterized thus far. Herein, we hypothesized that differentially expressed circRNAs could participate in the tumorigenesis and metastasis of HNSCC. To explore the roles of circRNAs, we examined the different expression patterns of circRNAs and mRNAs in HNSCC using microarray analyses. Subsequently, based on our microarray data and the circRNA-miRNA-mRNA interactions predicted with Arraystar's miRNA target prediction software based on TargetScan (31) and miRanda (32), a network was constructed using Cytoscape (33,34). OncomiR (35) and UALCAN (36) are online resources from The Cancer Genome Atlas (TCGA) specifically for miRNA and protein-coding gene analysis across many cancer types, respectively. DIANA-miRPath v3.0 (37) is an efficient tool to analyze the combinatorial effect of miRNAs on target pathways by P-value. According to the results obtained using these three tools, we predicted that circRNA_036186 may promote the expression of tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, ζ polypeptide (14-3-3ζ) by functioning as a sponge for the inhibitory miR-193b-3p. Our data indicate that circRNA_036186 may be suitable as a novel diagnostic marker and therapeutic target in HNSCC, as well as providing a basis for further studies regarding the role of circRNAs in HNSCC. Materials and methods Patients and samples This study was approved and supervised by The Ethics Committee of China Medical University (Shenyang, China). Five pairs of HNSCC and para-carcinoma tissues were collected from patients that had not received pre-operative chemotherapy or radiotherapy who underwent resective surgery at the Head and Neck Tumor Center (School of Stomatology, China Medical University). Written informed consent was obtained from all patients. All samples were obtained during surgery and immediately stored in liquid nitrogen. The diagnosis of HNSCC was confirmed by a pathological examination (Table I). RNA extraction Total RNA was extracted from the HNSCC and para-carcinoma tissue samples using TRIzol reagent (Invitrogen/Thermo Fisher Scientific, Inc., Waltham, MA, USA), according to the manufacturer's instructions. Subsequently, the RNA quantity and quality was measured using an ND-1000 system (NanoDrop; Thermo Fisher Scientific, Inc., Wilmington, DE, USA). RNA integrity was assessed via standard denaturing agarose gel electrophoresis methods. Microarray assay Five pairs of carcinoma and para-carcinoma tissues were used for a microarray assay, with the intention of identifying differentially expressed circRNAs and mRNAs by comparing the HNSCC samples with the adjacent control samples. Sample preparation and microarray hybridization were performed based on Arraystar's standard protocols. Briefly, total RNAs were digested with RNase R (Epicentre; Illumina, Inc., San Diego, CA, USA) to remove linear RNAs and enrich for circRNAs. Subsequently, the enriched circRNAs were amplified and transcribed into fluorescent cRNA, utilizing a random priming method (Arraystar Super RNA Labeling kit; Arraystar Inc., Rockville, MD, USA). The labeled cRNAs were hybridized onto the Arraystar Human circRNA Array V2 (8×15K; Arraystar Inc.). After washing the slides, the arrays were scanned using an Agilent G2505C scanner (Agilent Technologies, Inc., Santa Clara, CA, USA). Agilent Feature Extraction software (version 11.0.1.1) was used to analyze the acquired array images. The raw data were quantile normalized, and subsequent data processing was performed using the limma package in R. Significantly differentially expressed circRNAs and mRNAs between the two groups were identified through volcano plot filtering with the thresholds of P<0.05 and the Benjamini-Hochberg false discovery rate (FDR) of <0.5. Differentially expressed circRNAs and mRNAs between the HNSCC and adjacent control groups were identified through fold change filtering (circRNAs, |log fold change| >1.5; mRNAs, |log fold change| >2). Hierarchical clustering was performed to show the characteristics of the expression profiles based on the values of all expressed transcripts, and differentially expressed transcripts. Gene Ontology (GO) and pathway analyses The GO project provides a controlled vocabulary to describe gene and gene product attributes in any organism (http://www.geneontology.org/). The GO categories cover three different domains: Biological process (BP), cellular component (CC) and molecular function (MF). The P-value denotes the significance of GO terms enrichment in the differentially expressed genes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was used for harvesting pathway clusters covering the molecular interaction and reaction networks in differential gene expression profile. The P-value denotes the significance of the pathway association with the specified genes. Competing endogenous RNA network analysis The potential miRNA response elements were searched on the sequences of circRNAs and mRNAs, besides a measure with the number of common miRNAs, a hypergeometric test was executed for each ceRNA pair separately (38). The overlap of the same miRNA seed sequence binding site both on the circRNA and mRNA was predicted with the following parameters: miRNA coverage ≥0.1; context +< -4.999999977648258e-2; species = human; common Num ≥1; ceRNA type = protein coding; structure >140; P<0.05; context <-4.999999977648258e-2; energy <-10. miRNA coverage describes the threshold value for both SeqMiRNA coverage and CeMiRNA coverage. Context+ is the threshold value of the sum of the context+ scores used in TargetScan from version 6.0. Common Num is the threshold value of the number of common miRNAs between the ceRNA gene and the gene of interest. Structure is the threshold value of the sum of the structure scores used in miRanda. P-value is the threshold for the P-values calculated using hypergeometric tests. Context is the threshold value of the sum of the context scores used in TargetScan before version 5.x. Energy is the threshold value of the sum of the free energy predicted by miRanda. Exploring miRNA dysregulation in tumor development and progression miRNA dysregulation in tumor development and progression was predicted using OncomiR (www.oncomir.org (35) based on TCGA across 30 cancer types. Significance in tumor development was determined through a paired Student's t-test between normal and tumor tissues; and for survival, an unpaired Student's t-test and univariate Cox proportional hazards analysis between living and deceased patients. Prediction of potentially functional circRNA-miRNA-mRNA axes The miRNA pathway investigation was carried out based on DIANA-miRPath v3.0 (http://www.microrna.gr/miRPathv3) (37). Firstly, we searched miRNA candidates in DIANA-miRPath v3.0-DIANA-TarBase v7.0 (39) and microT-CDS v5.0 (40,41) to predict the potential target genes and pathways with thresholds of P<0.05 and MicroT score >0.8. Subsequently, elaborated signaling pathway networks based on KEGG database were constructed gathering target genes and interacting proteins. Finally, the expression differences for target genes in malignant tumors were analyzed using UALCAN (http://ualcan.path.uab.edu) (36). Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) validation The microarray results were confirmed by RT-qPCR with the same RNA samples as those used for microarray analysis. Following RNA isolation, SuperScript™ III Reverse Transcriptase (Invitrogen/Thermo Fisher Scientific, Inc., Waltham, MA, USA) was used to synthesize the cDNA according to the manufacturer's instructions. Subsequently, qPCR was performed using a ViiA 7 Real-time PCR System (Applied Biosystems/Thermo Fisher Scientific, Inc., Foster City, CA, USA) with a total reaction volume of 10 μl, including 0.5 μl PCR forward primer (10 μM), 0.5 μl PCR reverse primer (10 μM), 2 μl cDNA, 5 μl 2X Master Mix and 2 μl double distilled water. The thermocycling conditions were as follows: 95°C for 10 min, then 10 sec at 95°C and 60 sec at 60°C for a total of 40 cycles. β-actin was used as an internal reference. The experiments were run using three independent wells per condition. For quantitative analyses, the relative expression levels of circRNAs and mRNAs were calculated using the 2−∆∆Cq (42) method. The primers used were as follows: circRNA_036186 forward, 5′-CTGAAGCACCGCCCAGCT-3′ and reverse, 5′-GACGAGCCACATTCATTCCAG-3′; 14-3-3ζ forward, 5′-TGTTGTAGGAGCCCGTAG-3′ and reverse, 5′-GCAACCTCAGCCAAGTAA-3′; β-actin (H) forward, 5′-GTGGCCGAGGACTTTGATTG-3′ and reverse, 5′-CCTGTAACAACGCATCTCATATT-3′. Statistical analysis All data are expressed as the means ± SD. The statistical significance of the differentially expressed genes as determined by RT-qPCR was estimated using the Student's t-test. P<0.05 was considered to indicate a statistically significant difference. Results Identification of differentially expressed circRNA and mRNA profiles by microarray assay High-throughput microarray analysis was performed to detect circRNA and mRNA expression in the HNSCC tissues. In the present study, a total of 12,366 circRNAs and 35,252 mRNAs were detected by microarray probes in five paired HNSCC and normal tissues. As a result, 287 dysregulated circRNAs were detected. Among these, 146 and 141 circRNAs were upregulated and downregulated, respectively. Additionally, 1,053 mRNAs were found to be differentially regulated, among which 377 mRNAs were upregulated, while 676 mRNAs were downregulated. Scatter and Volcano Plots, tools for visualizing differential expression between two compared groups, were produced; the former was constructed using fold change values (Fig. 1A and B), while the latter considered both fold change and statistical significance (Fig. 1C and D). Hierarchical clustering analysis indicated that the circRNA and mRNA expression patterns were distinguishable between the groups (Fig. 2). The distribution of the dysregulated circRNAs on human chromosomes was summarized, and the results revealed that these circRNAs were distributed among all chromosomes, apart from the Y sex chromosome (Fig. 3B). According to the associations with the protein-coding genes, the dysregulated circRNAs were classified into five categories: 214 were exonic, 35 were intronic, 7 were antisense, 3 were intergenic and 28 were sense overlapping (Fig. 3A). Thus, the expression profiles of the circRNAs and mRNAs in the HNSCC tissues were determined to differ from those in the matched tumor-adjacent tissues. GO term and KEGG and pathway enrichment analyses To date, the functions of the majority circRNAs have not been annotated. The functional prediction of circRNAs is largely based on the annotations of their interacting protein-coding genes. Therefore, the GO term and KEGG pathway enrichment analyses of dysregulated mRNAs can reveal the roles of significantly differentially expressed circRNAs, to a certain extent. According to the results of a GO term enrichment analysis, the principally enriched significant BP terms included immune system process, immune response and signaling transduction, such as 'antigen processing and presentation of exogenous peptide antigen' (GO:0002478), 'lymphocyte mediated immunity' (GO:0002449) and G-protein coupled receptor signaling pathway (GO:0007186). The most enriched CC terms were mostly regarding membranes, such as 'integral component of membrane' (GO:0016021), 'plasma membrane' (GO:0005886) and 'extracellular membrane-bounded organelle' (GO:0065010). Furthermore, the most enriched MF terms were associated with receptor activity and receptor binding, including 'signaling receptor activity' (GO:0038023), 'signal transducer activity' (GO:0004871), 'activating transcription factor binding' (GO:0033613) and 'calcium-dependent protein binding' (GO:0048306). Fig. 4A shows the top 10 most significantly enriched GO terms in the BP, CC and MF categories for the up- and downregulated mRNAs. Pathway analysis was performed with the KEGG database. A total of 31 pathways associated with the upregulated mRNAs, and 6 related to the downregulated mRNAs were identified, including 'antigen processing and presentation' (hsa04612), 'natural killer cell mediated cytotoxicity' (hsa04650) and 'calcium signaling pathway' (hsa04020), among others. The most enriched pathways for the up- and downregulated protein-coding genes are listed in Fig. 4B. These biological processes, molecular functions and pathways likely contribute to the occurrence and development of HNSCC. Construction of circRNA-miRNA-mRNA network According to the ceRNA hypothesis, the members of each ceRNA group compete for the same MREs to regulate each other. We selected the top five upregulated and downregulated circRNAs as ranked by fold change from the microarray data (Table II), and then identified the potential ceRNA interactions via various bioinformatics methods. A summary of ceRNAs, including 5 circRNAs, 385 miRNAs and 5,148 mRNAs, was produced by merging the common targeting miRNAs. To further investigate the interactions between these ceRNAs in HNSCC tissues, we used only the differentially expressed mRNAs from the mRNA microarray data, and a specific ceRNA network, consisting of five circRNAs, 385 miRNAs and 96 mRNAs, for HNSCC was established (Fig. 5). Exploring miRNAs associated with tumor development and patient survival miRNAs have been shown to play a crucial role in the pathogenesis and development of cancer (5-8). OncomiR online software (http://www.oncomir.org/) was used to systematically identify the miRNAs associated with development, staging and overall survival in 30 cancer types, and these data were used to identify the miRNAs associated with tumor formation and patient survival in HNSCC. The results revealed that there were 376 miRNAs associated with tumor development and 172 miRNAs related to tumor prognosis in HNSCC. Among these, 80 overlapping miRNAs were found in both the tumor development and prognosis groups, suggesting that these miRNAs were likely to be involved in the regulation of both HNSCC occurrence and metastasis. Finally, we retained the miRNAs associated with both tumor development and patient survival from the ceRNA network (Fig. 6A and Table III), and thus constructed a network consisting of 4 circRNAs, 8 miRNAs and 51 mRNAs (Fig. 6B). Prediction of responsible and functional miRNA-mRNA axes To explore the molecular mechanisms of action of circRNAs in HNSCC, we predicted the circRNA-miRNA axes associated with cancer-related pathways. DIANA-miRPath v3.0 allows the analysis of miRNA-associated pathways as well as the prediction of miRNA targets using DIANA-microT-CDS v5.0 and/or DIANA-TarBase v7.0. We queried with the 8 miRNAs from our network, including miR-377-5p, -432-5p, -125b-5p, -629-3p, -532-3p, -181c-5p, -193b-3p and -135b-5p, in DIANA-miRPath v3.0-DIANA-TarBase v7.0 and microT-CDS v5.0 to predict the potential target genes and pathways, with thresholds of P<0.05 and MicroT score >0.8. The overlapping results from DIANA-TarBase v7.0 and DIANA-microT-CDS v5.0, and the aforementioned analyses on ceRNAs revealed that miR-193b-3p was matched with 14-3-3ζ according to the results of all three methods, indicating that miR-193b-3p had a significant potential to interact with 14-3-3ζ (Fig. 7A). Therefore, 14-3-3ζ and miR-193b-3p were analyzed, and the results from pathway intersection analysis indicated that miR-193b-3p was likely to participate in the Hippo signaling pathway (hsa04390) (P=5.442667e-06) by binding to 14-3-3ζ in HNSCC (Fig. 7B). Prediction of the potential function of circRNA_036186 regarding 14-3-3ζ expression According to GO, 14-3-3ζ may act through 'protein targeting' (GO:0006605), 'anti-apoptosis' (GO:0006916) and 'signal transduction' (GO:0007165). All these GO biological processes are associated with the occurrence and development of cancer. We analyzed the expression differences for miR-193b-3p and 14-3-3ζ in malignant tumors using OncomiR and UALCAN, respectively. The results revealed that the expression levels of miR-193b-3p and 14-3-3ζ were both upregulated (P<0.01) in the six malignant tumor samples compared with in the adjacent tissues (Fig. 8A and B). However, evidence has indicated that the mRNA expression level is negatively associated with its corresponding miRNA (4), and we hypothesized that this effect may be attributed to the participation of circRNA expression. As a miRNA sponge, circRNAs can counteract, and even reverse, the functional differences caused by an alteration in miRNA expression. In this study, according to the results of ceRNA analysis, there was an evolutionarily conserved regulatory region (7mer-m8) in circRNA_036186 for binding to miR-193b-3p (Fig. 8C), and the expression of circRNA_036186 was significantly higher in the HNSCC tumor tissues compared with the para-carcinoma tissues in the microarray data. Therefore, we proposed that the circRNA_036186-miR-193b-3p-14-3-3ζ axis was likely crucial in the occurrence and prognosis of a variety of cancer types, including HNSCC. Furthermore, our results from microarray analysis revealed that the level of 14-3-3ζ was upregulated in HNSCC, compared with the control, which was consistent with the data from TCGA (36). This indicated that circRNA_036186 probably serves as a crucial and positive regulator of 14-3-3ζ and is significantly engaged in the Hippo signaling pathway in HNSCC. Validation of the differential expression levels of circRNA_036186 and 14-3-3ζ The results of microarray analysis of circRNA_036186 and 14-3-3ζ were verified by RT-qPCR with the same RNA samples used for the micro-array analysis. The results of RT-qPCR assay revealed that the expression levels of circRNA_036186 and 14-3-3ζ were each upregulated in HNSCC compared with the controls. The results are consistent with those of the microarray assay (Table IV); therefore, the RT-qPCR data verified the accuracy of the microarray results. These results provide evidence for our hypothesis that circRNA_036186 participates in the pathogenesis and development of HNSCC. Discussion As a group of functional molecules with abundant biological information, RNAs have particular value in the prevention and treatment of human diseases. In eukaryotic cells, protein coding RNA accounts for only 2% of the genome, and the vast majority of transcripts are ncRNAs (43). In recent years, the biogenesis and functions of ncRNAs have attracted increasing attention (44). The malignant behaviors of cancer cells are usually regulated by multiple gene products, and, theoretically, these genes can be grouped into 'networks' based on their interactions (10,33,45,46). These interactions may provide new insight into the potential mechanisms of neoplasia and lay the groundwork for a new approach to cancer treatment. miRNAs are the most well-studied ncRNAs, and are a class of small ncRNAs ~22 nucleotides (nt) in length that, by binding to MREs, can block protein translation or modulate mRNA stability on a post-transcriptional level (4,47,48). Unlike miRNAs, the functions and regulatory mechanisms of circRNAs remain largely unknown. Fortunately, the application of circRNA microarray techniques has allowed for the identification of differentially expressed circRNAs and the exploration of their functions (49). An increasing number of circRNAs have been found to be cell type- or tissue-specific, rendering them more predictive of the potential biological functions (50,51). Evidence has indicated that the functions of circRNAs stem from three main aspects: Absorbing miRNA sponge (24,29,30,52), binding RNA-binding proteins (RBPs) (16) and translating peptides (53,54). To date, the miRNA sponge theory remains the most popular hypothesis. HNSCC is an aggressive tumor that often occurs with recurrence and metastasis, and the survival rate has changed little in recent decades. Over the past few years, an increasing number of experimental and clinical studies have suggested that ncRNAs play an important role in HNSCC pathophysiology (55 and refs. therein). For instance, miR-150-5p and miR-150-3p, antitumor miRNAs, have been shown to significantly inhibit cell aggressiveness by directly targeting SPARC (osteonectin), cwcv and kazal like domains proteoglycan 1 (SPOCK1) in HNSCC (56). miR-206 inhibited cell proliferation by arresting the cell cycle in S-phase and targeting histone deacetylase 6 (HDAC6) via the phosphatase and tensin homolog (PTEN)/AKT/mammalian target of rapamycin (mTOR) pathway (57). Additionally, by evaluating the associations between miRNA polymorphisms and HNSCC risk according to the cancer site, miR-605 rs2043556 and miR-196a2 rs11614913 have been shown as likely to affect genetic susceptibility to oral squamous cell carcinoma (58). Evidence has suggested that circRNA_100290 may serve as a sponge of the miR-29 family to regulate cyclin-dependent kinase 6 (CDK6) expression in oral squamous cell carcinoma (59). In the present study, microarray screening and RT-qPCR verification were used to identify circRNA and mRNA expression profiles, and ceRNA analysis was performed to evaluate the circRNA-miRNA-mRNA interactions. The results of microarray analysis and RT-qPCR revealed that circRNA_036186 and 14-3-3ζ were highly expressed in HNSCC compared with the controls. Using bioinformatics analysis, we found that circRNA_036186 and 14-3-3ζ had the same seed sequence binding site to for interaction with miR-193b-3p. It has been previously reported that miR-193b-3p and 14-3-3ζ are closely related to the occurrence and prognosis of cancer (60-64). Existing evidence for miR-193b-3p has been largely focused on the negative regulation of the expression of oncogenes in cancer development (65-67), and the effect on carcinogenesis and prognosis is the most well-characterized aspect of 14-3-3ζ (68-70). Furthermore, 14-3-3ζ mRNA and protein expression have been reported to be upregulated in HNSCC tissue samples, and the overexpression of 14-3-3ζ promotes cell growth, as well as morphological changes. By contrast, reduced 14-3-3ζ levels increase the relative proportion of cells in the G1/G0-phase (71). The results from DIANA-miRPath v3.0 revealed that miR-193b-3p significantly engaged in the Hippo signaling pathway via 14-3-3ζ with the function of cell contact inhibition and organ size control by combining with yes-associated protein 1 (YAP)/WW domain-containing transcription regulator protein 1 (TAZ). A growing body of evidence has suggested that the Hippo pathway may act as a critical role in organ size control by promoting cell growth, cell proliferation and preventing apop-tosis (72-75). A recent study reported that the interaction between 14-3-3ζ and YAP, as a negative regulatory feedback loop, may play an important role in cell proliferation and survival in gastric cancer (76). Moreover, using TCGA data, we identified a group of miRNAs that were relevant to both cancer development and patient survival in HNSCC. The result revealed that the expression of miR-193b-3p in HNSCC tissues was significantly higher than that in the adjacent tissues. However, to date, at least to the best of our knowledge, there have been no reports on the functions and mechanisms of the circRNA_036186-miR-193b-3p-14-3-3ζ pathway in patients with HNSCC. Considering the data from previous studies in addition to the results of the current study, the emerging endogenous circRNA molecule circRNA_036186 probably serves as a crucial regulator of 14-3-3ζ by acting as a miRNA sponge to inhibit miR-193b-3p function, and hence participates in the Hippo signaling pathway to exert an effect on the course of HNSCC. In this study, due to the limited number of available tissues samples, we only analyzed circRNA and mRNA expression profiles in five paired HNSCC tissues and their normal controls. More samples should be obtained and analyzed in the future to verify our results. Additionally, our results were mainly based on bioinformatics predictions; whether they can be extended to wider applications remains undefined. The main limitation of the present study was that at present, there is no convincing functional study to verify that circRNA_036186 inhibited miR-193b-3p to result in regulation of the Hippo signaling pathway; further functional studies to verify the roles of circRNA_036186 are required. In conclusion, the findings of this study revealed the circRNA and mRNA expression characteristics in HNSCC and highlighted the possibility that circRNA-036186 may serve as a diagnostic and therapeutic biomarker in HNSCC. Acknowledgments Not applicable. Abbreviations circRNAcircular RNA miRNA or miRmicroRNA ncRNAnon-coding RNA RT-qPCRreverse transcription-quantitative polymerase chain reaction GOGene Ontology KEGGKyoto Encyclopedia of Genes and Genomes MREmiRNA response element ceRNAcompeting endogenous RNA UTRuntranslated region RBPRNA-binding protein TCGAThe Cancer Genome Atlas BLCAbladder urothelial carcinoma BRCAbreast invasive carcinoma HNSCChead and neck squamous cell carcinoma LUADlung adenocarcinoma LUSClung squamous cell carcinoma UCECuterine corpus endometrial carcinoma YAPyes-associated protein 1 TAZWW domain-containing transcription regulator protein 1 14-3-3ζtyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, ζ polypeptide HDAC6histone deacetylase 6 SPOCK1sparc/osteonectin, cwcv and kazal-like domains proteoglycan 1 CDK6cyclin-dependent kinase 6 Funding This study was supported by grants from the National Natural Science Foundation of China (no. 81372877), the National Young Scholars Science Foundation of China (no. 81102058), the Foundation of Education Bureau of Liaoning Province (no. L2014317), the Natural Science Foundation of Liaoning Province (no. 2014021096), the Excellent Talent Fund Project of Higher Education of Liaoning Province (no. LJQ2014087), the Doctoral Scientific Research Foundation of Liaoning Province (no. 201501002), and the Young and Middle-aged Science and Technology Innovation Talent Support Project of Shenyang (no. RC170489). Availability of data and materials The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Authors' contributions WLW, CFS and FYL conceived and designed the experiments. ZY and YJZ carried out the experiments. WLW analyzed the data, prepared the figures and wrote the manuscript. PL and YKN were mainly involved in the collection and pathological examinations of patient specimens. CFS supervised the whole project. All authors have read and approved the final manuscript. Ethics approval and consent to participate This study was approved and supervised by The Ethics Committee of China Medical University (Shenyang, China). Written informed consent was obtained from all patients. Patient consent for publication Not applicable. Competing interests The authors declare that they have no competing interests. Figure 1 Differences in circRNA and mRNA expression profiles between HNSCC tissues and tumor-adjacent tissues. (A) Scatter plot of circRNAs: The values on the x- and y-axes are the averaged normalized signal values of groups of samples (log2-scaled). The green lines are fold change lines. The circRNAs above the top line and below the bottom line indicate >1.5-fold change between the two compared samples. (B) Scatter plot of mRNAs: The fold change cut-off was 2.0. (C) Volcano plot of circRNAs: The vertical lines correspond to 1.5-fold up and down, respectively, and the horizontal line represents the P-value cut-off of 0.05. The red points in the plot represent the differentially expressed circRNAs with statistical significance. (D) Volcano plot of mRNAs: The red and green points represent the differentially expressed mRNAs based on the thresholds of fold change >2.0 and P<0.05. circRNA, circular RNA; HNSCC, head and neck squamous cell carcinoma. Figure 2 Heatmap showing the circRNA and mRNA expression profiles. Maps on the left are based on expression values of all circRNAs and mRNAs detected by microarray probes. Maps on the right correspond to normalized expression values of significantly changed circRNAs and mRNAs. The expression values are depicted in line with the color scale. The intensity increases from green to red. Each column represents one sample, and each row indicates a transcript. C, C3, C6, C8 and C9 represent cancer tissues, and N, N3, N6, N8 and N9 represent normal control tissues. (A) Hierarchical cluster analysis of all expressed circRNAs and dysregulated circRNAs with fold change >1.5 and P<0.05. (B) Hierarchical cluster analysis of all expressed mRNAs and significantly changed mRNAs with fold change >2.0 and P<0.05. circRNA, circular RNA. Figure 3 Characterizations of circRNAs in HNSCC. (A) Bar graph showing the types and counts of differentially regulated circRNAs. The circRNAs were classified into five types according to the relationship with their associated coding genes. (B) Distribution of differentially expressed circRNAs in human chromosomes. circRNA, circular RNA; HNSCC, head and neck squamous cell carcinoma. Figure 4 GO and pathway analyses. (A) GO annotation of up- and downregulated mRNAs with the top 10 enrichment scores, covering the domains of biological process, cellular component and molecular function, respectively. (B) KEGG pathway enrichment analysis of dysregulated mRNAs. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes. Figure 5 circRNA-miRNA-mRNA network. In this ceRNA network, each gene corresponds to a node, and 2 interacting genes based on seed sequence pairings are connected by a solid line. The 'V' shapes represent circRNAs; the diamond shapes represent miRNAs; the ellipses represent mRNAs. The blue and gray lines represent circRNA-miRNA and miRNA-mRNA interactions, respectively. circRNA, circular RNA; miRNA, microRNA; ceRNA, competing endogenous RNA. Figure 6 miRNAs associated with tumor development and patient survival. (A) Venn diagram indicating that the predicted miRNAs derived from the results of OncomiR and ceRNA analyses. miR-377-5p, miR-432-5p, miR-125b-5p, miR-629-3p, miR-532-3p, miR-181c-5p, miR-193b-3p and miR-135b-5p were intersecting. (B) The circRNA-miRNA-mRNA network consists of 4 circRNAs, 8 microRNAs and 51 mRNAs. circRNA, circular RNA; miRNA, microRNA; ceRNA, competing endogenous RNA. Figure 7 miRNA target prediction and miRNA-mediated signaling pathways. (A) Venn diagrams showing the target genes of miRNAs predicted by DIANA-microT-CDS v5.0, DIANA-TarBase v7.0 and miRNA target prediction software based on TargetScan and miRanda. The results indicated that NM_003406 (red) was targeted by miR-193b-3p is intersected. (B) Mapping of the Hippo signaling pathway indicated that 14-3-3ζ (red) mediated by miR-193b-3p serves as a crucial factor significantly engaged in the function of cell-contact inhibition and organ size control by combining with yes-associated protein 1/WW domain-containing transcription regulator protein 1 (YAP/TAZ). miRNA, microRNA. Figure 8 Expression changes of RNA transcripts and circRNA-miRNA interaction. (A) Bar graph showing the log2 transform of the mean miR-193b-3p expression in both tumor (red) and non-tumor tissues (blue). **P<0.01. (B) Boxplot showing relative expression of 14-3-3ζ in normal and cancer samples. **P<0.01. (C) The interaction of circRNA_036186-miR-193b-3p predicted based on data from TargetScan and miRanda. BRCA, breast invasive carcinoma; HNSCC, head and neck squamous cell carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamouscell carcinoma; UCEC, uterine corpus endometrial carcinoma. Table I Clinicopathological characteristics of the five patients. Patient no. Sex Age Position Pathological diagnosis Histological grade TNM stage 1 Male 53 Abdomen of the tongue and bottom of the oral cavity HNSCC G2 T2N0M0 3 Male 33 Abdomen of the tongue HNSCC G2 T2N0M0 6 Female 63 Lateral border of tongue HNSCC G2 T2N1M0 8 Male 59 Lateral border of tongue HNSCC G2 T2N1M0 9 Male 43 Lateral border of tongue HNSCC G1-G2 T2N1M0 HNSCC, head and neck squamous cell carcinoma. Table II Top 5 up- and downregulated circRNAs ranked by fold change >3.0, P<0.01 and FDR <0.05 in the microarray data. circRNA Regulation Absolute fold change P-value FDR Source Alias Chrom circRNA type hsa_circRNA_014280 Up 4.005569 4.01E-05 0.017333277 circBase hsa_circ_0014280 chr1 Exonic hsa_circRNA_402089 Up 3.2867219 4.64E-05 0.018507674 25242744 _ chr19 Exonic hsa_circRNA_036186 Up 3.2503065 7.05E-05 0.022943867 circBase hsa_circ_0036186 chr15 Exonic hsa_circRNA_404474 Up 3.1449976 8.22E-05 0.024789243 25070500 _ chr1 Exonic hsa_circRNA_102485 Down 4.3172347 1.71E-05 0.009591013 circBase hsa_circ_0050102 chr19 Exonic 'Regulation' describes that the experimental group has greater or lower intensity values than the control group. 'Absolute fold change' indicates the absolute ratio (no log scale) of normalized intensities between two conditions. 'P-value' indicates the P-value calculated from the paired t-test. 'FDR' indicates that FDR was calculated from the Benjamini-Hochberg false discovery rate. 'Alias' indicates the circRNA ID in circBase (http://circbase.mdc-berlin.de). 'circRNA type' indicates the circRNAs classified into 5 types (namely exonic, intronic, antisense, intergenic and sense overlapping. circRNA, circular RNA. Table III miRNAs associated with both tumor development and overall survival in the ceRNA network. miRNA Cancer Tumor development Tumor progression t-test P-value t-test FDR Upregulated in: Log-rank P-value Log-rank FDR Z-score hsa-miR-432-5p HNSCC 0.00177a 0.00458 Normal 0.000457a 0.143 3.496 hsa-miR-181c-5p HNSCC 0.00723a 0.0166 Normal 0.00159a 0.106 3.18 hsa-miR-377-5p HNSCC 0.000000119a 0.000000861 Normal 0.00334a 0.154 3.017 hsa-miR-125b-5p HNSCC 7.03E-11a 1.23E-09 Normal 0.00775a 0.416 2.672 hsa-miR-193b-3p HNSCC 3.41E-13a 1.47E-11 Tumor 0.0285a 0.686 2.163 hsa-miR-532-3p HNSCC 0.000438a 0.00127 Tumor 0.0314a 0.174 2.156 hsa-miR-629-3p HNSCC 2.38E-10a 3.32E-09 Tumor 0.0389a 0.183 2.059 The Student's paired t-test was used to compare expression between tumor and normal tissues. Log-rank P-values and Z-scores obtained from univariate Cox proportional hazards analysis were used to show the relative effect of miRNA expression levels in regards to patient survival time. a P≤0.05. Table IV Microarray and RT-qPCR analyses of selected circRNAs and mRNAs. Detection method hsa_circRNA_036186 14-3-3ζ Absolute fold change Regulation P-value Absolute fold change Regulation P-value Microarray 3.250307 Up 0.000071a 2.318862 Up 0.021603a RT-qPCR 2.836309 Up 0.001582a 1.7610711 Up 0.011234a The expression levels of circRNA_036186 and 14-3-3ζ were determined by RT-qPCR. Each RT-qPCR assay was performed at least three times. The RT-qPCR and microarray data were in accord. a P<0.05. circRNA, circular RNA. ==== Refs References 1 Marur S Forastiere AA Head and neck squamous cell carcinoma: Update on epidemiology, diagnosis, and treatment Mayo Clin Proc 91 386 396 2016 10.1016/j.mayocp.2015.12.017 26944243 2 Siegel RL Miller KD Jemal A Cancer statistics, 2016 CA Cancer J Clin 66 7 30 2016 10.3322/caac.21332 26742998 3 Chen W Zheng R Baade PD Zhang S Zeng H Bray F Jemal A Yu XQ He J Cancer statistics in China, 2015 CA Cancer J Clin 66 115 132 2016 10.3322/caac.21338 26808342 4 Fabian MR Sonenberg N Filipowicz W Regulation of mRNA translation and stability by microRNAs Annu Rev Biochem 79 351 379 2010 10.1146/annurev-biochem-060308-103103 20533884 5 Li H Xie S Liu M Chen Z Liu X Wang L Li D Zhou Y The clinical significance of downregulation of mir-124-3p, mir-146a-5p, mir-155-5p and mir-335-5p in gastric cancer tumorigenesis Int J Oncol 45 197 208 2014 10.3892/ijo.2014.2415 24805774 6 Cao J Cai J Huang D Han Q Yang Q Li T Ding H Wang Z miR-335 represents an invasion suppressor gene in ovarian cancer by targeting Bcl-w Oncol Rep 30 701 706 2013 10.3892/or.2013.2482 23708561 7 Croce CM Causes and consequences of microRNA dysregulation in cancer Nat Rev Genet 10 704 714 2009 10.1038/nrg2634 19763153 8 Esquela-Kerscher A Slack FJ Oncomirs - microRNAs with a role in cancer Nat Rev Cancer 6 259 269 2006 10.1038/nrc1840 16557279 9 Poliseno L Salmena L Zhang J Carver B Haveman WJ Pandolfi PP A coding-independent function of gene and pseudogene mRNAs regulates tumour biology Nature 465 1033 1038 2010 10.1038/nature09144 20577206 10 Salmena L Poliseno L Tay Y Kats L Pandolfi PP A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell 146 353 358 2011 10.1016/j.cell.2011.07.014 21802130 11 Rybak-Wolf A Stottmeister C Glažar P Jens M Pino N Giusti S Hanan M Behm M Bartok O Ashwal-Fluss R Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed Mol Cell 58 870 885 2015 10.1016/j.molcel.2015.03.027 25921068 12 Qu S Yang X Li X Wang J Gao Y Shang R Sun W Dou K Li H Circular RNA: A new star of noncoding RNAs Cancer Lett 365 141 148 2015 10.1016/j.canlet.2015.06.003 26052092 13 Jeck WR Sorrentino JA Wang K Slevin MK Burd CE Liu J Marzluff WF Sharpless NE Circular RNAs are abundant, conserved, and associated with ALU repeats RNA 19 141 157 2013 10.1261/rna.035667.112 23249747 14 Zhang Y Zhang XO Chen T Xiang JF Yin QF Xing YH Zhu S Yang L Chen LL Circular intronic long noncoding RNAs Mol Cell 51 792 806 2013 10.1016/j.molcel.2013.08.017 24035497 15 Shen T Han M Wei G Ni T An intriguing RNA species - perspectives of circularized RNA Protein Cell 6 871 880 2015 10.1007/s13238-015-0202-0 26349458 16 Li Z Huang C Bao C Chen L Lin M Wang X Zhong G Yu B Hu W Dai L Exon-intron circular RNAs regulate transcription in the nucleus Nat Struct Mol Biol 22 256 264 2015 10.1038/nsmb.2959 25664725 17 Chen LL The biogenesis and emerging roles of circular RNAs Nat Rev Mol Cell Biol 17 205 211 2016 10.1038/nrm.2015.32 26908011 18 Ashwal-Fluss R Meyer M Pamudurti NR Ivanov A Bartok O Hanan M Evantal N Memczak S Rajewsky N Kadener S circRNA biogenesis competes with pre-mRNA splicing Mol Cell 56 55 66 2014 10.1016/j.molcel.2014.08.019 25242144 19 Wilusz JE Sharp PA Molecular biology. A circuitous route to noncoding RNA Science 340 440 441 2013 10.1126/science.1238522 23620042 20 Vicens Q Westhof E Biogenesis of circular RNAs Cell 159 13 14 2014 10.1016/j.cell.2014.09.005 25259915 21 Zhang XO Wang HB Zhang Y Lu X Chen LL Yang L Complementary sequence-mediated exon circularization Cell 159 134 147 2014 10.1016/j.cell.2014.09.001 25242744 22 Westholm JO Miura P Olson S Shenker S Joseph B Sanfilippo P Celniker SE Graveley BR Lai EC Genome-wide analysis of drosophila circular RNAs reveals their structural and sequence properties and age-dependent neural accumulation Cell Reports 9 1966 1980 2014 10.1016/j.celrep.2014.10.062 25544350 23 Suzuki H Aoki Y Kameyama T Saito T Masuda S Tanihata J Nagata T Mayeda A Takeda S Tsukahara T Endogenous multiple exon skipping and back-splicing at the DMD mutation hotspot Int J Mol Sci 17 17 2016 10.3390/ijms17101722 24 Hansen TB Jensen TI Clausen BH Bramsen JB Finsen B Damgaard CK Kjems J Natural RNA circles function as efficient microRNA sponges Nature 495 384 388 2013 10.1038/nature11993 23446346 25 Zheng Q Bao C Guo W Li S Chen J Chen B Luo Y Lyu D Li Y Shi G Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs Nat Commun 7 11215 2016 10.1038/ncomms11215 27050392 26 Guo JU Agarwal V Guo H Bartel DP Expanded identification and characterization of mammalian circular RNAs Genome Biol 15 409 2014 10.1186/s13059-014-0409-z 25070500 27 Tay Y Rinn J Pandolfi PP The multilayered complexity of ceRNA crosstalk and competition Nature 505 344 352 2014 10.1038/nature12986 24429633 28 Ghosal S Das S Sen R Basak P Chakrabarti J Circ2Traits: A comprehensive database for circular RNA potentially associated with disease and traits Front Genet 4 283 2013 10.3389/fgene.2013.00283 24339831 29 Lukiw WJ Circular RNA (circRNA) in Alzheimer's disease (AD) Front Genet 4 307 2013 10.3389/fgene.2013.00307 24427167 30 Hansen TB Kjems J Damgaard CK Circular RNA and miR-7 in cancer Cancer Res 73 5609 5612 2013 10.1158/0008-5472.CAN-13-1568 24014594 31 Enright AJ John B Gaul U Tuschl T Sander C Marks DS MicroRNA targets in Drosophila Genome Biol 5 R1 2003 10.1186/gb-2003-5-1-r1 14709173 32 Pasquinelli AE MicroRNAs and their targets: Recognition, regulation and an emerging reciprocal relationship Nat Rev Genet 13 271 282 2012 10.1038/nrg3162 22411466 33 Assenov Y Ramírez F Schelhorn SE Lengauer T Albrecht M Computing topological parameters of biological networks Bioinformatics 24 282 284 2008 10.1093/bioinformatics/btm554 18006545 34 Shannon P Markiel A Ozier O Baliga NS Wang JT Ramage D Amin N Schwikowski B Ideker T Cytoscape: A software environment for integrated models of biomolecular interaction networks Genome Res 13 2498 2504 2003 10.1101/gr.1239303 14597658 35 Wong NW Chen Y Chen S Wang X OncomiR: An online resource for exploring pan-cancer microRNA dysregulation Bioinformatics 34 713 715 2018 10.1093/bioinformatics/btx627 29028907 36 Chandrashekar DS Bashel B Balasubramanya SAH Creighton CJ Ponce-Rodriguez I Chakravarthi BVSK Varambally S UALCAN: A portal for facilitating tumor subgroup gene expression and survival analyses Neoplasia 19 649 658 2017 10.1016/j.neo.2017.05.002 28732212 37 Vlachos IS Zagganas K Paraskevopoulou MD Georgakilas G Karagkouni D Vergoulis T Dalamagas T Hatzigeorgiou AG DIANA-miRPath v3.0: Deciphering microRNA function with experimental support Nucleic Acids Res 43W1 W460 W466 2015 10.1093/nar/gkv403 38 Li JH Liu S Zhou H Qu LH Yang JH starBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data Nucleic Acids Res 42D1 D92 D97 2014 10.1093/nar/gkt1248 39 Vlachos IS Paraskevopoulou MD Karagkouni D Georgakilas G Vergoulis T Kanellos I Anastasopoulos IL Maniou S Karathanou K Kalfakakou D DIANA-TarBase v7.0: indexing more than half a million experimentally supported miRNA:mRNA interactions Nucleic Acids Res 43D D153 D159 2015 10.1093/nar/gku1215 40 Paraskevopoulou MD Georgakilas G Kostoulas N Vlachos IS Vergoulis T Reczko M Filippidis C Dalamagas T Hatzigeorgiou AG DIANA-microT web server v5.0: Service integration into miRNA functional analysis workflows Nucleic Acids Res 41W W169 W173 2013 10.1093/nar/gkt393 41 Reczko M Maragkakis M Alexiou P Grosse I Hatzigeorgiou AG Functional microRNA targets in protein coding sequences Bioinformatics 28 771 776 2012 10.1093/bioinformatics/bts043 22285563 42 Livak KJ Schmittgen TD Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method Methods 25 402 408 2001 10.1006/meth.2001.1262 11846609 43 Alexander RP Fang G Rozowsky J Snyder M Gerstein MB Annotating non-coding regions of the genome Nat Rev Genet 11 559 571 2010 10.1038/nrg2814 20628352 44 Mercer TR Dinger ME Mattick JS Long non-coding RNAs: Insights into functions Nat Rev Genet 10 155 159 2009 10.1038/nrg2521 19188922 45 Mei Y Yang JP Qian CN For robust big data analyses: A collection of 150 important pro-metastatic genes Chin J Cancer 36 16 2017 10.1186/s40880-016-0178-z 28109319 46 Qian CN Mei Y Zhang J Cancer metastasis: Issues and challenges Chin J Cancer 36 38 2017 10.1186/s40880-017-0206-7 28372569 47 Carthew RW Sontheimer EJ Origins and mechanisms of miRNAs and siRNAs Cell 136 642 655 2009 10.1016/j.cell.2009.01.035 19239886 48 He L Hannon GJ MicroRNAs: Small RNAs with a big role in gene regulation Nat Rev Genet 5 522 531 2004 10.1038/nrg1379 15211354 49 Caiment F Gaj S Claessen S Kleinjans J High-throughput data integration of RNA-miRNA-circRNA reveals novel insights into mechanisms of benzo[a]pyrene-induced carcinogenicity Nucleic Acids Res 43 2525 2534 2015 10.1093/nar/gkv115 25690898 50 Xia S Feng J Lei L Hu J Xia L Wang J Xiang Y Liu L Zhong S Han L Comprehensive characterization of tissue-specific circular RNAs in the human and mouse genomes Brief Bioinform 18 984 992 2017 27543790 51 Salzman J Chen RE Olsen MN Wang PL Brown PO Cell-type specific features of circular RNA expression PLoS Genet 9 e1003777 2013 10.1371/journal.pgen.1003777 24039610 52 Memczak S Jens M Elefsinioti A Torti F Krueger J Rybak A Maier L Mackowiak SD Gregersen LH Munschauer M Circular RNAs are a large class of animal RNAs with regulatory potency Nature 495 333 338 2013 10.1038/nature11928 23446348 53 Granados-Riveron JT Aquino-Jarquin G The complexity of the translation ability of circRNAs Biochim Biophys Acta 1859 1245 1251 2016 10.1016/j.bbagrm.2016.07.009 27449861 54 Wang Y Wang Z Efficient backsplicing produces translatable circular mRNAs RNA 21 172 179 2015 10.1261/rna.048272.114 25449546 55 Sannigrahi MK Sharma R Panda NK Khullar M Role of non-coding RNAs in head and neck squamous cell carcinoma: A narrative review Oral Dis 9 21 2017 Epub ahead of print 10.1111/odi.12782 28941018 56 Koshizuka K Hanazawa T Kikkawa N Katada K Okato A Arai T Idichi T Osako Y Okamoto Y Seki N Antitumor miR-150-5p and miR-150-3p inhibit cancer cell aggressiveness by targeting SPOCK1 in head and neck squamous cell carcinoma Auris Nasus Larynx 45 854 865 2018 10.1016/j.anl.2017.11.019 29233721 57 Liu F Zhao X Qian Y Zhang J Zhang Y Yin R miR-206 inhibits head and neck squamous cell carcinoma cell progression by targeting HDAC6 via PTEN/AKT/mTOR pathway Biomed Pharmacother 96 229 237 2017 10.1016/j.biopha.2017.08.145 28987947 58 Miao L Wang L Zhu L Du J Zhu X Niu Y Wang R Hu Z Chen N Shen H Association of microRNA polymorphisms with the risk of head and neck squamous cell carcinoma in a Chinese population: A case-control study Chin J Cancer 35 77 2016 10.1186/s40880-016-0136-9 27515039 59 Chen L Zhang S Wu J Cui J Zhong L Zeng L Ge S circRNA_100290 plays a role in oral cancer by functioning as a sponge of the miR-29 family Oncogene 36 4551 4561 2017 10.1038/onc.2017.89 28368401 60 Rauhala HE Jalava SE Isotalo J Bracken H Lehmusvaara S Tammela TL Oja H Visakorpi T miR-193b is an epigenetically regulated putative tumor suppressor in prostate cancer Int J Cancer 127 1363 1372 2010 10.1002/ijc.25162 20073067 61 Leivonen SK Rokka A Ostling P Kohonen P Corthals GL Kallioniemi O Perälä M Identification of miR-193b targets in breast cancer cells and systems biological analysis of their functional impact Mol Cell Proteomics 10 M110.005322 2011 10.1074/mcp.M110.005322 21512034 62 Mitra AK Chiang CY Tiwari P Tomar S Watters KM Peter ME Lengyel E Microenvironment-induced downregulation of miR-193b drives ovarian cancer metastasis Oncogene 34 5923 5932 2015 10.1038/onc.2015.43 25798837 63 Zhao GY Ding JY Lu CL Lin ZW Guo J The overexpression of 143-3ζ and Hsp27 promotes non-small cell lung cancer progression Cancer 120 652 663 2014 10.1002/cncr.28452 24804299 64 Rehman SK Li SH Wyszomierski SL Wang Q Li P Sahin O Xiao Y Zhang S Xiong Y Yang J 14-3-3ζ orchestrates mammary tumor onset and progression via miR-221-mediated cell proliferation Cancer Res 74 363 373 2014 10.1158/0008-5472.CAN-13-2016 24197133 65 Li J Kong F Wu K Song K He J Sun W miR-193b directly targets STMN1 and uPA genes and suppresses tumor growth and metastasis in pancreatic cancer Mol Med Rep 10 2613 2620 2014 10.3892/mmr.2014.2558 25215905 66 Mets E Van der Meulen J Van Peer G Boice M Mestdagh P Van de Walle I Lammens T Goossens S De Moerloose B Benoit Y MicroRNA-193b-3p acts as a tumor suppressor by targeting the MYB oncogene in T-cell acute lymphoblastic leukemia Leukemia 29 798 806 2015 10.1038/leu.2014.276 25231743 67 Zhang J Qin J Su Y miR-193b-3p possesses anti-tumor activity in ovarian carcinoma cells by targeting p21-activated kinase 3 Biomed Pharmacother 96 1275 1282 2017 10.1016/j.biopha.2017.11.086 29169729 68 Matta A Bahadur S Duggal R Gupta SD Ralhan R Over-expression of 14-3-3zeta is an early event in oral cancer BMC Cancer 7 169 2007 10.1186/1471-2407-7-169 17764575 69 Nishimura Y Komatsu S Ichikawa D Nagata H Hirajima S Takeshita H Kawaguchi T Arita T Konishi H Kashimoto K Overexpression of YWHAZ relates to tumor cell proliferation and malignant outcome of gastric carcinoma Br J Cancer 108 1324 1331 2013 10.1038/bjc.2013.65 23422756 70 Murata T Takayama K Urano T Fujimura T Ashikari D Obinata D Horie-Inoue K Takahashi S Ouchi Y Homma Y 14-3-3zeta, a novel androgen-responsive gene, is upregulated in prostate cancer and promotes prostate cancer cell proliferation and survival Clin Cancer Res 18 5617 5627 2012 10.1158/1078-0432.CCR-12-0281 22904106 71 Lin M Morrison CD Jones S Mohamed N Bacher J Plass C Copy number gain and oncogenic activity of YWHAZ/14-3-3zeta in head and neck squamous cell carcinoma Int J Cancer 125 603 611 2009 10.1002/ijc.24346 19405126 72 Yao CB Zhou X Chen CS Lei QY The regulatory mechanisms and functional roles of the Hippo signaling pathway in breast cancer Yi Chuan 39 617 629 2017 28757476 73 Patel SH Camargo FD Yimlamai D Hippo signaling in the liver regulates organ size, cell fate, and carcinogenesis Gastroenterology 152 533 545 2017 10.1053/j.gastro.2016.10.047 28003097 74 Ji XY Zhong G Zhao B Molecular mechanisms of the mammalian Hippo signaling pathway Yi Chuan 39 546 567 2017 28757470 75 Pan D The hippo signaling pathway in development and cancer Dev Cell 19 491 505 2010 10.1016/j.devcel.2010.09.011 20951342 76 Zhang B Gong A Shi H Bie Q Liang Z Wu P Mao F Qian H Xu W Identification of a novel YAP-14-3-3ζ negative feedback loop in gastric cancer Oncotarget 8 71894 71910 2017 29069755