Skip to main content

Expression profile of circular RNA and construction of circular RNA-Micro RNA network in salivary adenoid cystic carcinoma

Abstract

Background

Circular RNAs (circRNAs) is a newly discovered type of non-coding RNA, the abnormal expression of which has been demonstrated in many types of human tumors. So they have been considered as promising candidates as diagnostic and therapeutic targets in cancer. This research aimed to screen the profile of circRNA expression in salivary adenoid cystic carcinoma (SACC).

Methods

Using the threshold of FDR < 0.05 and fold change > 2 or < 0.5, 5 up-regulated and 26 down-regulated circRNAs were identified. The reliability of sequencing was verified by the expression detection of randomly selected circRNAs via qRT-PCR.

Results

Moreover, the circRNA-miRNA system was established by bioinformatics approaches and successfully identified an interaction between circRNA ABCA13 and a cancer-related miRNA (miR-138-5p), which was also verified by qRT-PCR. Moreover, the predicted molecular interaction proved that circRNA ABCA13 may promote SACC through inhibition of miR-138-5p.

Conclusions

Collectively, this study has offered the first report about the circRNA expression profile and circRNA-miRNA network in SACC. All of the above could benefit the exploration of novel therapeutic target in SACC treatment.

Background

Salivary adenoid cystic carcinoma (SACC) is regarded as one of the most frequently diagnosed malignant salivary gland tumors, holding 1% of all head and neck malignant tumors, and accounting for 30.42% of all malignant epithelial tumors of the salivary gland [1, 2]. Clinicopathological studies have suggested that SACC is extremely invasive, inducing distal hematogenous metastasis through blood vessel invasion, inducing perineural invasion through the infiltration of nerves. Statistics showed that the organ metastasis rate of SACC was as high as 40.9% [3, 4]. Despite some remarkable achievements have been made recently, the etiology and pathogenesis of SACC remain largely unclear [5]. Over the past years, many studies concentrated in seeking biomarkers to improve the clinical prognosis of SACC, but the results showed that current chemotherapy and molecular targeting drugs can only achieve temporary partial response or stable disease in advanced SACC with local recurrence or distant metastasis [6, 7]. Accordingly, currently the mortality of SACC still remains high and the prognosis still remains poor [3]. Therefore, a deeper understanding of SACC-related biomarkers is of great importance in revealing the pathogenesis of SACC and could provide basis for the progression of novel targeted therapies with higher efficacy.

The emergence of new generation of sequencing technology, high-throughput sequencing, is a revolutionary improvement in sequencing technology, through which more and more important genes in the tumorigenesis and development of human cancers have been screened out [8,9,10,11]. Moreover, high-throughput transcriptome sequencing provides a promising method to study the pathogenesis of SACC. It is of great significance to elucidate the molecular mechanism of SACC and find specific molecular markers for early diagnosis as well as intervention of tumor.

Up to now, the abnormal expression of various non-coding RNAs including lncRNAs and miRNAs has been demonstrated to be closely correlated to the pathogenesis of human cancers [12,13,14]. Circular RNA (circRNA) is one of the most recently discovered non-coding RNAs forming a closed-loop circRNA via binding between 3′ and 5′ ends [15, 16]. Researches have suggested that circRNAs could evade the degradation of nucleic acid exonucleases, thus having a high degree of stability and resistance to RNA degradation pathways, which endows them the promising potential to act as cancer biomarkers [17,18,19,20]. Previous studies have confirmed that circRNAs were differentially expressed in tumors, thereby playing the role of oncogene or tumor suppressor gene, and are closely related to biological characteristics of tumors. Moreover, circRNAs are able to inhibit of many cancer-related miRNAs expression, and the circRNA-miRNA-mRNA interaction cascade acts a regulatory role in several cancer-related pathways, which can either inhibit or promote the development and progress of cancer [21]. As far as we konw, circRNA-related studies have rarely been reported in the field of SACC, and therefore no detailed description of the relationship between circRNA-miRNA-mRNA cascade and SACC has been made. This study aimed to identify novel circRNA biomarkers by recognizing SACC-related abnormal expression of circRNA by next generation RNA sequencing analysis. The further analysis of the relationship between these circRNAs and their potential association with miRNAs could further promote our understanding of the SACC pathogenesis and propose potential therapy targets for SACC treatment.

Methods

Patient enrollment and tissues collection

Samples came from SACC patients who had finished the operation in the Department of Maxillofacial Head and Neck Oncology, Ninth People’s Hospital Affiliated to Medical College of Shanghai Jiaotong University from March 2016 to January 2017. People receiving any cancer treatment before admission were kept out from this study. After tumors resection, tumors and adjacent tissues were frozen in liquid nitrogen, and stored at − 80 °C. Histopathological examination as well as diagnosis of tumors were performed independently by two pathologists. Six matched SACC tissue samples and corresponding para-carcinoma tissues were examined by secondary sequencing, which were also investigated by real-time qPCR. The research was approved by the Ethics Committee of Shanghai Ninth People's Hospital and all participants in this research has signed the informed consent form.

Sample preparation

Total RNA from cell was extracted by Trizol reagent (Invitrogen, CA, USA). The total RNA was quantified using Bioanalyzer 2200 (Agilent Technologies, CA, USA), and then stored at −  80 °C.

cDNA library construction

The cDNA libraries were established by VAHTSTM Total RNA-seq (H/M/R). Briefly, divalent cations was used for depletion of rRNA and fragmenting into 150–200 bp at 94 °C for 8 min, and RNA fragments were then reverse-transcribed into first-strand cDNA, and second-strand cDNA was synthesized, fragments were finally end repaired VAHTSTM DNA Clean Beads were used for harvesting target bands. The products were purified and enriched by PCR which were used for creating the final cDNA libraries then quantified by Agilent2200, and then pooling the labeled cDNA libraries in equal proportion, sequenced in a single lane of the Illumina HiSeqTM 2500 with 51 plus 7 cycles by NovelBio Corp. Laboratory, Shanghai, China to generate 150 bp paired-end.

RNA sequencing mapping

Clean reads were obtained from the original reads via expurgating the adaptor sequences before mapping, over 20% of bases with qualities of < 20 were contained by low-quality reads as well as reads with > 5% vague bases (noted as N), and the clean reads were then aligned to Human genome (version: GRCh38) through the HISAT2 program [22].

CircRNA identification and quantification

The pipeline ‘acfs’ which was able to be used widely at https://code.google.com/p/acfs/ was used for identifying each samples’ circRNA comprising of the reported steps [23]. BOWTIE2 version 2.2.5 was used for mapping the respective reference genome [GRCH37.p13 NCBI]. Unmapped Reads were used for identifying the circRNA through BWA mem (bwa mem -t 1-k 16-T 20): local alignments of fragments in a single read that mapped to (i) regions of the same chromosome which are less than 1 Mb apart from each other (ii) on the same chain (iii) however, in the opposite order, kept them as candidates which supported junction between head to tail. Intensity of potential splice sites was then estimated, which were supported by these junction reads via MaxEntScan33, and determined the accurate juncture by choosing the donor and receptor sites with the highest splicing intensity score. Finally we reported candidate circRNAs if at least 2 reads held out the head-to-tail junction and the splicing score ≥ 10.

Expression analysis

To evaluate circRNA expression, the BWA-mem was used for realigning all the unmapped reads to circRNA candidates under the following parameters (bwa mem-t 1-k 16-T 20). We concatenated sequences at the 5′ to 3′ end to form circular junctures, and counted the readings of each candidate which were mapped to the juncture (with at least 6 nt overhang).

Dif-gene-finder

DESeq algorithm was utilized for filtering the differentially expressed genes, and the following are the threshold represented by FDR and P value analysis:

  1. i.

    Fold Change > 2 or < 0.5;

  2. ii.

    FDR < 0.05.

GO analysis

The GO annotations from UniProt (http://www.uniprot.org/), the Gene Ontology (http://www.geneontology.org/) and NCBI (http://www.ncbi.nlm.nih.gov/) were used. In addition, we applied Fisher’s exact test for identifying the significant GO categories and used FDR for correcting the P values.

Pathway analysis

The pathway analysis was used for looking for the differential genes’ significant pathway in accordance with KEGG database. Our group selected the significant pathway through passing the Fisher’s exact test, and made use of P value and FDR for defining the threshold of significance.

GO-tree

We built GO-tree on the basis of the Gene Ontology Directed Acyclic Graph for providing user with friendly visualization and data navigation. The significant GO-Term (P Value < 0.01) was selected in GO Analysis in accordance with the significantly up- and down-regulated genes expression for constructing the GO-Tree in order to sum up the functions which were influenced in this assay.

Path-act-network

KEGG was comprised of membrane transport, signal transduction, cell cycle pathways and metabolism. Genes were picked in abound biological pathway and Cytoscape was used to graphically represent the pathways, and then we made use of the KEGG database for building a genetic network based on the proteins, compounds in the database and the relationship among the genes.

Target analysis

We utilize the Miranda as the tools for predicting differentially expressed miRNA Target on circRNA, lncRNA and mRNA.

Real time-PCR

Total RNA was extracted via Trizol (Invitrogen, CA, USA) and Direct-zol™ RNA MiniPrep (Zymo Research, Orange, CA, USA) according to the instructions of the manufacturer, and quantified the RNA using Nanodrop 2000 (Thermo Fisher Scientific, MA, USA). We obtained first-strand cDNA by M-MLV reverse transcriptase (RT) (Promega), arbitrary primers and 1 µg of total RNA and the cDNA templates were stored at − 20 °C. The PCR products were separated by 20 g/L agarose gel, stained by ethidium bromide, observed under ultraviolet light. qPCR was performed using AceQ qPCR SYBR Green Master Mix (Vazyme). 2−∆∆Ct method was used for the quantification of expressions. In this experiment, β-actin was chosen for the internal control.

Statistical analysis

Our group expressed the data as mean ± SD (n = 6) and analyzed it by GraphPad Prism 6 software (GraphPad Software Inc., San Diego, CA, USA). In addition, t-test was used for comparing the difference. It was considered to be statistical significance that P value < 0.05.

Results

Identification of differentially expressed circRNAs

Our group performed the high-throughput sequencing for evaluating circRNA expression in the SACC samples (T group) in comparison with the matched normal samples (G group). The results showed that 49475 circRNAs were detected through high-throughput sequencing. The normalized log2 scales were shown by the scatter plot to illustrate the differences between the two groups of samples (Fig. 1a) and we also displayed the identified circRNAs with deregulated expression in the volcano plot (Fig. 1b). In Table 1, the differentially circRNAs expression and the related information were listed. As demonstrated, 31 differentially expressed circRNAs were detected through using the fold-change filtering (absolute log2 (fold change) > 1, and FDR < 0.05). Compared to the matched normal tissues, 26 circRNAs were markedly down-regulated and 5 circRNAs were up-regulated in SACC tissues. Moreover, the sequencing results showed that, among the 31 differentially expressed circRNAs, 21 of them (17 down-regulated and 4 up-regulated) were covered in the exon of genomic, 2 of them (1 down-regulated and 1 up-regulated) were in the intronic genomic, while the rest were unknown (Fig. 1c). Furthermore, the chromosome distribution of the detected circRNA and differentially expressed circRNA was shown in Fig. 1d. Finally, our group utilized hierarchical clustering to illustrate the expression patterns of both all the circRNAs and the differentially circRNAs expression in the samples (Fig. 2).

Fig. 1
figure1

Characterizations in circRNA expression patterns of SACC and matched normal tissues. a The difference in the expression of circRNAs in SACC and matched normal tissues was shown by scattered plot. b Volcano plot shows the differential expression of circRNA between the two groups. The red dots indicated the up-regulated circRNAs, the blue dots indicated the down-regulated circRNAs. c The type of differentially expressed circRNAs. d The distribution of all circRNAs and the differentially expressed circRNAs in human chromosomes

Table 1 The differentially expressed circRNAs identified by sequencing
Fig. 2
figure2

Hierarchical clusters of all circRNAs and differentially expressed circRNAs in two groups. a All circRNAs. b Differentially expressed circRNAs

Validation of the reliability of the sequencing data

For the aim to validate the reliability of the RNA sequencing, 3 up-regulated circRNA ABCA13, circRNA CASC15 and circRNA VCAN, and 1 down-regulated circRNA LIFR were selected randomly as the target of detection whose expression levels were discovered by qRT-PCR in the SACC tissues and matched normal tissues. As illustrated by Fig. 3, the qRT-PCR results were concordant with the sequencing data, indicating that the circRNA profile had high reliability.

Fig. 3
figure3

Verification of the reliability of RNA sequencing. qRT-PCR was performed for detecting the expression of circRNA circRNA ABCA13, circRNA CASC15, circRNA VCAN and circRNA LIFR in SACC and matched normal tissues. Data are shown as mean ± SD. * P < 0.05, ** P < 0.01, *** P < 0.001

GO and KEGG analysis of differentially expressed circRNAs

In order for further exploring the functionalities of the differentially expressed circRNAs, our group performed GO enrichment analysis and KEGG pathway analysis on the ceRNA-related mRNAs with the dysregulated circRNAs based on negative correlation between circRNA and miRNA, the negative correlation between the mRNA and the miRNA, and thus the positive correlation between circRNA and mRNA. As shown in Fig. 4a, GO analysis demonstrated that canonical Wnt signaling pathway, which plays important role in human cancers, and multicellular organismal development were two of the biological processes that were most likely related to the differentially expressed circRNAs. Moreover, the analysis indicated that the dysregulated circRNAs could also have impact in the molecular functions such as proline-rich region binding and phospholipid scramblase activity, and cellular components such as cell junction and tight junction (Fig. 4b, c). Subsequently, the KEGG analysis for the circRNAs showed the strong enrichment in proteoglycan in cancer and pathways in cancer such as Wnt signaling pathway and Hippo signaling pathway (Fig. 4d). Moreover, the bulb map of KEGG analysis for the up-regulated circRNAs also showed distinct enrichment in pathways in cancer and that for the down-regulated circRNAs suggested a significant enrichment in salivary secretion (Fig. 5). Taken together, the results of GO and KEGG analysis provided the potential mechanism of circRNAs in SACC.

Fig. 4
figure4

GO and KEGG pathway analysis. a Biological processes of GO database. b Cellular component of GO database. c Molecular functions of GO database. d KEGG pathway analysis

Fig. 5
figure5

Bulb map of KEGG analysis. a Bulb map of KEGG analysis for up-regulated circRNAs. b Bulb map of KEGG analysis for down-regulated circRNAs

Establishment of the circRNA-miRNA network

On the basis of the results of RNA sequencing and bioinformatics analysis, miRANDA and TargetScan databases were used for predicting miRNA targets of each conserved seed matching sequence. The results showed that there is a fairly complex regulation network between circRNA and miRNA, and Fig. 6a showed the interaction network constructed by cytoscape based on the selected pairs of circRNA-miRNA with one-to-one correspondence. Then, the interaction pairs involving exonic circRNAs were further screened and shown in Fig. 6b. Among them, we found that circRNA ABCA13 and circRNA LIFR, whose expression levels have both been verified by qRT-PCR, were predicted to interact with miR-138-5p and miR-520a-5p, respectively. In the following verification of the miRNAs by qRT-PCR, the down-regulation of miR-138-5p in SACC tissues was consistent with the predicted interaction (Fig. 6c), while the simultaneously up-regulated expression of circRNA LIFR and miR-520a-5p was not in accordance with the prediction (data not shown). The predicted molecular interaction between circRNA ABCA13 and miR-138-5p was displayed in Fig. 7a. Moreover, we also predict the regulation network of miR-138-5p and showed the molecular interaction between miR-138-5p and Notch1, which has been reported to have important role in SACC (Fig. 7b, c).

Fig. 6
figure6

The circRNA-miRNA network. a The general view of circRNA-miRNA network for the differentially expressed circRNAs. b The circRNA-miRNA network for the exonic differentially expressed circRNAs. c The expression of miR-138-5p was detected by qRT-PCR in SACC and matched normal tissues. Data are shown as mean ± SD. * P < 0.05, ** P < 0.01, *** P < 0.001

Fig. 7
figure7

The circRNA-miRNA interaction. a The predicted interaction between circRNA ABCA13 and miR-138-5p. b The predicted ceRNA network of miR-138-5p. c The predicted interaction between miR-138-5p and Notch1

Discussion

CircRNA, a novel kind of non-coding RNA with regulatory ability, is different from traditional mRNA and was first discovered in RNA viruses in the 1970s [15]. As the high-throughput RNA sequence technology and bioinformatics further developed, researchers have found that circRNA is ubiquitous in various species and has important physiological functions [15, 24]. Moreover, previous results have manifested that some circRNAs may play crucial parts in the development and progression of human cancers [25]. Sequencing analysis of the expression profiles of circRNA in tumor tissues and the matched adjacent tissues showed that abnormal circRNAs expression was common in tumors [25]. For example, Lu et al. identified 1155 differentially expressed circRNAs through the microarray detection of human breast cancer tissues and corresponding para-carcinoma tissues. They further verified that hsa_circ_100219 was differentially expressed in breast cancer tissues and normal tissues, and of important value in the diagnosis of breast cancer [26]. Another research showed that hsa_circ_0001649 expression was inhibited in hepatocellular carcinoma and positively related to the growth of tumor. Subsequent study indicated that it could regulate various genes expression, thus participating in the incidence and metastasis of hepatocellular carcinoma [27]. Moreover, researchers found a great amount of complete and stable circRNAs in human serum exosomes, and identified the significantly different expression profile of circRNAs in serum of cancer patents and normal people [18]. Collectively, all the results indicated that aberrant circRNAs expression may be taken part in the incidence and development of various human cancer.

Current studies have confirmed that miRNAs are related to almost all cellular functions and play a crucial role in the development of diseases [28]. CircRNA has abundant microRNA (miRNA) response element (MRE), which plays the role of miRNA sponge in cells, and can eliminate the regulation of miRNA on its target genes [29]. Thus, circRNAs may involve in the physiological regulation and human diseases through the interaction with miRNA. Because of the post-transcriptional regulatory effect of circRNA, circRNA has become an ideal therapeutic target, which can be used in the development of new drugs through the simulation of its molecular characteristics, providing new means for treating diseases [30]. However, the circRNA research in oncology is still in its initial stage and the relationship between circRNA and SACC is rarely reported.

In this research, we present the expression profile of circRNA in SACC for the first time by using secondary sequencing, attempting to further understand the molecular mechanism of SACC. Based on the results of RNA sequencing, it revealed that 31 circRNAs were dysregulated in SACC tissues compared to the matched normal tissues based on the filtration by the absolute of log2 fold change > 1 and FDR < 0.05. Among them, 26 circRNAs were significantly down-regulated and 5 circRNAs were significantly up-regulated. 4 differentially expressed circRNAs were randomly chosen for verifying the dependability of the sequencing data through detecting their expression levels in the tissue samples by qRT-PCR. The results of qRT-PCR were in perfectly consistent with the RNA sequencing data, clearly suggesting the validation of RNA sequencing.

Considering that accumulating evidence has proved that circRNA was able to regulate their parent genes expression [4], we performed GO and KEGG pathway analysis for evaluating the function of circRNAs in SACC. GO analysis revealed that axon guidance, proline-rich region binding and cell junction were the functions that are most likely regulated by the differentially expressed circRNAs among biological processes (BP), molecular functions (MF) and cellular components (CC), respectively. Notably, the Wnt signaling pathway, which has been proved to take part in various human cancers, was also potentially regulated by the differentially expressed circRNAs. On the other hand, the KEGG analysis indicated the relationship between the dysregulated circRNAs and several tumor-associated signaling pathways such as Wnt pathway, Hippo pathway and Ras pathway. Moreover, the bulb map of the up-regulated circRNAs showed the most potential involvement of pathways in cancer and that for the down-regulated circRNAs was salivary secretion, indicating their potential regulation mechanism of SACC.

In this study, we also constructed the circRNA-miRNA interaction network in SACC for the first time to further investigate the functional roles of circRNAs. Through screening the circRNA-miRNA interaction pair with one-to-one correspondence, we found that miR-138-5p and miR-520a-5p were the targets of circRNA ABCA13 and circRNA LIFR, respectively. The verification by qRT-PCR also proved the down-regulation of miR-138-5p expression in SACC tissues, which was in accordance with the negative regulation effect between circRNA and miRNA. Notably, miR-138-5p was recently reported to be able to act as a tumor inhibitor in colorectal cancer through targeting PD-L1 [31]. Therefore, circRNA ABCA13 may promote the development and progression of SACC through the inhibition of miR-138-5p expression. Interestingly, the molecular interaction analysis demonstrated not only the interaction between circRNA ABCA13 and miR-138-5p, but also the interaction between miR-138-5p and Notch1, which was recently suggested to play a critical role in the cell growth, anti-apoptosis and metastasis of SACC [32]. Despite of all the results, the detailed regulation effect and underlying mechanism of circRNA ABCA13/miR-138-5p axis still need further investigation and would be our focus in the future work.

Conclusion

In conclusion, this study illustrated the first circRNA signature of SACC. Through the high throughput transcriptome sequencing and bioinformatics analysis, we recognized the differentially expressed circRNAs, analyzed the potential functional roles of them and constructed the circRNA-miRNA interaction network. Moreover, we demonstrated that circRNA ABCA13 may promote the development and progression of SACC through the inhibition of miR-138-5p expression. Our findings benefit for the exploration of novel therapeutic target and provide potential diagnostic markers for SACC treatment.

Availability of data and materials

The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.

References

  1. 1.

    Bell D, Bell AH, Bondaruk J, Hanna EY, Weber RS. In-depth characterization of the salivary adenoid cystic carcinoma transcriptome with emphasis on dominant cell type. Cancer-Am Cancer Soc. 2016;122:1513–22.

    CAS  Google Scholar 

  2. 2.

    Tian Z, Li L, Wang L, Hu Y, Li J. Salivary gland neoplasms in oral and maxillofacial regions: a 23-year retrospective study of 6982 cases in an eastern Chinese population. Int J Oral Max Surg. 2010;39:235–42.

    CAS  Article  Google Scholar 

  3. 3.

    Dillon PM, Chakraborty S, Moskaluk CA, Joshi PJ, Thomas CY. Adenoid cystic carcinoma: a review of recent advances, molecular targets, and clinical trials. Head Neck. 2016;38:620–7.

    PubMed  Article  Google Scholar 

  4. 4.

    Zhang Y, Zhang X, Chen T, Xiang J, Yin Q, Xing Y, Zhu S, Yang L, Chen L. Circular Intronic Long Noncoding RNAs. Mol Cell. 2013;51:792–806.

    CAS  Article  Google Scholar 

  5. 5.

    Wysocki PT, Izumchenko E, Meir J, Ha PK, Sidransky D, Brait M. Adenoid cystic carcinoma: emerging role of translocations and gene fusions. Oncotarget. 2016;7:66239–54.

    PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Andrs CP, Rodrigo JP, Bradley PJ, Vincent VP, Asterios T, Hunt JL, Primo S, Alessandra R, Missak H, Takes RP. Adenoid cystic carcinoma of the head and neck–An update. Oral Oncol. 2015;51:652–61.

    Article  Google Scholar 

  7. 7.

    Dodd RL, Slevin NJ. Salivary gland adenoid cystic carcinoma: a review of chemotherapy and molecular therapies. Oral Oncol. 2006;42:759–69.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Giordano TJ, Rork K, Tobias E, Gauger PG, Michelle V, Juliane B, Donita S, Thomas DG, Gerard D, Gary H. Molecular classification and prognostication of adrenocortical tumors by transcriptome profiling. Clin Cancer Res. 2009;15:668–76.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Peng G, Lin CJ, Mo W, Dai H, Park YY, Kim SM, Peng Y, Mo Q, Siwko S, Hu R. Genome-wide transcriptome profiling of homologous recombination DNA repair. Nat Commun. 2014;5:3361.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Reuter J, Spacek DV, Snyder M. High-Throughput Sequencing Technologies. Mol Cell. 2015;58:586–97.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Underwood JG, Uzilov AV, Katzman S, Onodera CS, Mainzer JE, Mathews DH, Lowe TM, Salama SR, Haussler D. FragSeq: transcriptome-wide RNA structure probing using high-throughput sequencing. Nat Methods. 2010;7:995–1001.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Alessandro F, Irene B. Long non-coding RNAs: new players in cell differentiation and development. Nat Rev Genet. 2014;15:7–21.

    Google Scholar 

  13. 13.

    Li W, Notani D, Rosenfeld MG. Enhancers as non-coding RNA transcription units: recent insights and future perspectives. Nat Rev Genet. 2016;17:207–23.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Matsui M, Corey DR. Non-coding RNAs as drug targets. Nat Rev Drug Discov. 2016;16:167.

    PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Sebastian M, Marvin J, Antigoni E, Francesca T, Janna K, Agnieszka R, Luisa M, Sebastian DM, Lea HG, Mathias M. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333–8.

    Article  Google Scholar 

  16. 16.

    Hsiao KY, Sun HS, Tsai SJ. Circular RNA-new member of noncoding RNA with novel functions. Exp Biol Med. 2017;242:1136.

    CAS  Article  Google Scholar 

  17. 17.

    Lei Y, Gong X, Lei S, Zhou Q, Lu B, Zhu L. The circular RNA Cdr1as act as an oncogene in hepatocellular carcinoma through targeting miR-7 expression. PLoS ONE. 2016;11:e158347.

    Google Scholar 

  18. 18.

    Li Y, Zheng Q, Bao C, Li S, Guo W, Zhao J, Chen D, Gu J, He X, Huang S. Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis. Cell Res. 2015;25:981–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Yang W, Du WW, Li X, Yee AJ, Yang BB. Foxo3 activity promoted by non-coding effects of circular RNA and Foxo3 pseudogene in the inhibition of tumor growth and angiogenesis. Oncogene. 2015;35:3919.

    PubMed  Article  Google Scholar 

  20. 20.

    Zheng Q, Bao C, Guo W, Li S, Jie C, Bing C, Luo Y, Lyu D, Yan L, Shi G. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 2016;7:11215.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    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. 2015;365:141–8.

    CAS  Article  Google Scholar 

  22. 22.

    Qu S, Song W, Yang X, Wang J, Zhang R, Zhang Z, Zhang H, Li H. Microarray expression profile of circular RNAs in human pancreatic ductal adenocarcinoma. Genomics Data. 2015;5:385–7.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Xintian Y, Irena V, Ana B, Tristan W, Irina E, Georgi T, Güney A, Mantian W, Caspar G, Claudia Q. Neural circular RNAs are derived from synaptic genes and regulated by development and plasticity. Nat Neurosci. 2015;18:603–10.

    Article  Google Scholar 

  24. 24.

    Miri D, Schraga S, Sarit E, Rotem S. Transcriptome-wide discovery of circular RNAs in Archaea. Nucleic Acids Res. 2012;40:3131–42.

    Article  Google Scholar 

  25. 25.

    Jingqiu L, Jie Y, Ping Z, Yanping L, Chengwei Z, Shaomin W, Dazhi X, Hui-Kuan L, Zhaohui G. Circular RNAs in cancer: novel insights into origins, properties, functions and implications. Am J Cancer Res. 2015;5:472–80.

    Google Scholar 

  26. 26.

    Lü L, Sun J, Shi P, Kong W, Xu K, He B, Zhang S, Wang J. Identification of circular RNAs as a promising new class of diagnostic biomarkers for human breast cancer. Oncotarget. 2017;8:44096–107.

    PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Qin M, Liu G, Huo X, Tao X, Sun X, Ge Z, Yang J, Fan J, Liu L, Qin W. Hsa_circ_0001649: A circular RNA and potential novel biomarker for hepatocellular carcinoma. Cancer Biomark. 2016;16:161–9.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Li M, Marinmuller C, Bharadwaj U, Chow KH, Yao Q, Chen C. MicroRNAs: Control and Loss of Control in Human Physiology and Disease. World J Surg. 2009;33:667–84.

    PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495:384–8.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Wei L, Feng T, Xing F, You Y, Yang J, Zhao ZA, Liu J, Shen Z, Deng W, Hu S. Signature of circular RNAs in human induced pluripotent stem cells and derived cardiomyocytes. Stem Cell Res Ther. 2018;9:56.

    Article  Google Scholar 

  31. 31.

    Zhao L, Yu H, Yi S, Peng X, Su P, Xiao Z, Liu R, Tang A, Li X, Liu F. The tumor suppressor miR-138-5p targets PD-L1 in colorectal cancer. Oncotarget. 2016;7:45370–84.

    PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Su BH, Qu J, Song M, Huang XY, Hu XM, Xie J, Zhao Y, Ding LC, She L, Chen J. NOTCH1 signaling contributes to cell growth, anti-apoptosis and metastasis in salivary adenoid cystic carcinoma. Oncotarget. 2014;5:6885–95.

    PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

This work was financially supported by National Natural Science Foundation of China (Grant number 81901049 and 81772870).

Funding

National Natural Science Foundation of China.

Author information

Affiliations

Authors

Contributions

MR and JL made substantial contributions to the concept and design of the present study. JH and NH conducted the experiments. JH, CZ and ZX conducted data analysis. JH and JL produced the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jiannan Liu or Min Ruan.

Ethics declarations

Ethics approval and consent to participate

The research was approved by the Ethics Committee of Shanghai Ninth People's Hospital and all participants in this research has signed the informed consent form.

Consent for publication

Not applicable.

Competing interests

The authors have no conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Han, J., Han, N., Xu, Z. et al. Expression profile of circular RNA and construction of circular RNA-Micro RNA network in salivary adenoid cystic carcinoma. Cancer Cell Int 21, 28 (2021). https://doi.org/10.1186/s12935-020-01681-2

Download citation

Keywords

  • Circular RNA
  • Salivary adenoid cystic carcinoma
  • microRNA
  • RNA sequencing
  • circRNA ABCA13