Circular RNA expression in pancreatic ductal adenocarcinoma
- Authors:
- Published online on: July 16, 2019 https://doi.org/10.3892/ol.2019.10624
- Pages: 2923-2930
-
Copyright: © Zhang et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
Pancreatic cancer is one of the top five causes of cancer-associated mortality with a fatality rate of 95%; the majority of patients with pancreatic cancer have a delayed diagnosis, as early diagnosis is difficult (1). Recent data showed that pancreatic cancer has the lowest 5-year relative survival rate (9%) among all cancer types (2). Pancreatic ductal adenocarcinoma (PDAC) accounts for a majority of pancreatic cancer cases and is an aggressive and difficult malignancy to treat (3). Despite advancing knowledge of the mechanisms of carcinogenicity and improvement in diagnosis and management, the prognosis remains poor for patients with PDAC (3). Therefore, new mechanisms and biomarkers are needed to improve the prognosis and detection of PDAC in patients.
MicroRNAs (miRNAs) are a class of small non-coding RNAs (ncRNAs), 18–25 nucleotides long, that negatively control gene expression at the mRNA and protein level (4). Accumulating evidence has suggested that miRNAs are dysregulated in various types of cancers, including PDAC (5). As the absence of symptoms in early disease results in a late diagnosis of PDAC, numerous miRNAs, such as miR-21, miR-155, miR-196 and miR-210, have been identified as biomarkers for PDAC (5). Understanding the mechanisms contributing to miRNA dysregulation may help to improve PDAC diagnosis and treatment.
Another class of ncRNAs are circular RNAs (circRNAs), which function as miRNA sponges, thus serving as competitive inhibitors and suppressing the ability of the miRNA to bind its mRNA targets (6,7). Although circRNAs were discovered in RNA viruses in the 1970s, these molecules were considered as molecular flukes or products of aberrant RNA splicing (7). Recent evidence indicates that circRNAs may serve important roles in cancers, including bladder cancer (8), esophageal squamous cell carcinoma (9), breast cancer (10), basal cell carcinoma (11), cutaneous squamous cell carcinoma (12), colorectal cancer (13) and PDAC (14). In the present study, microarray data downloaded from Gene Expression Omnibus (GEO) were used to explore the roles of circRNAs and miRNAs in PDAC and to provide novel insights into PDAC biology.
Materials and methods
Microarray data
Raw gene expression profile datasets GSE69362, GSE79634 and GSE43796 were downloaded from the GEO database (http://www.ncbi.nlm.nih.gov/geo). The datasets GSE69362 and GSE79634 were analyzed on the GPL19978 Agilent-069978 Arraystar Human CircRNA microarray V1 platform, which analyzes 5,396 human circRNAs at the same time. The dataset GSE43796 was analyzed on the GPL15159 Agilent-031181 Unrestricted Human miRNA V16.0 Microarray platform, which can detect 1,205 mature human miRNAs.
The GSE69362 dataset accessed the expression of circRNAs in six pairs of PDAC and adjacent normal tissues. The GSE79634 dataset accessed the expression of circRNAs in 20 pairs of samples. GSE43796 detected the expression of miRNAs in solid-pseudopapillary neoplasm of pancreas (n=14), ductal adenocarcinoma (n=6), neuroendocrine tumor (n=6) and non-neoplastic pancreas (n=5). The six PDAC tissues and the five non-neoplastic pancreases from GSE43796 were included in the present study.
Data analysis
Differentially expressed genes were identified using the limma package in the R software (version 3.3.2; http://www.r-project.org) (15). Raw data was read and background was corrected and normalized. Probes that were ≥10% brighter compared with the negative controls on at least one-half of arrays were averaged for each gene. The array weights were estimated and used in the linear model analysis. Paired or unpaired Student's t-test was used to identify differentially expressed genes between tumor tissues and normal tissues. The Benjamini-Hochberg method was used to correct for multiple comparisons (16). Differentially expressed genes between tumor tissues and normal tissues were identified, using fold-change >2 and adjusted P<0.05 as selection criteria.
miRNA prediction
As certain sequences of circRNAs in circBase may differ from the mature sequence of circRNAs (17), sequences of the circRNAs were extracted using circPrimer (version 1.2; http://www.bioinf.com.cn) (18), which was developed by our laboratory. As exon-derived circRNAs have been demonstrated to serve as miRNA sponges, the 256 differentially expressed circRNAs were annotated using circPrimer and only 214 exon-derived circRNAs were used for further analysis.
The sequences of the 214 circRNAs and 20 differentially expressed miRNAs were saved as FASTA format. circRNA/miRNA interactions were predicted using circMir (version 1.0; http://www.bioinf.com.cn), which used the miRanda (19) and RNAhybrid (20) algorithms to predict miRNAs that bind to spliced junctions of circRNAs. As the prediction algorithm often suffers from high false positive rates, only the circRNA/miRNA interactions predicted by the two independent algorithms were used. Cytoscape software version 3.1.1 (http://cytoscape.org) was utilized to construct a possible functional network of the differentially expressed circRNAs and miRNAs.
Pathway analysis of circRNA gene symbols
circRNA host genes were submitted to the web-based tool DAVID (http://david.abcc.ncifcrf.gov) to perform a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The following default thresholds were used: Count=2 and EASE (a modified Fisher's Exact test P-value)=0.1. The enriched pathways and gene symbols were imported into Cytoscape software version 3.1.1 to construct a network.
Results
Overview of the datasets
The distribution of normalized signal intensities of the datasets indicated that microarray quality for circRNA and miRNA was good (Fig. 1A and B, respectively). The expression levels of the majority of genes were correlated between PDAC and normal tissues, although there were a number of differentially expressed circRNAs and miRNAs (Fig. 1C and D, respectively). Volcano plots were generated to visualize the relationship of fold-change and statistical significance for circRNAs and miRNAs (Fig. 1E and F, respectively).
Identification of differentially expressed genes
Compared with adjacent normal tissues, there were 256 differentially expressed circRNAs in PDAC tissues, including 115 upregulated circRNAs and 141 downregulated circRNAs. The top 40 differentially expressed circRNAs in PDAC are presented in Table I, and the full list is presented in Table SI. Eight upregulated miRNAs and 12 downregulated miRNAs were identified in PDAC (Table II).
Prediction of circRNA/miRNA interaction
circRNAs may upregulate the expression of target genes of miRNAs by sequestering the miRNAs. In the present study, a total of 51 circRNA/miRNA interactions were identified and a possible functional network of the differentially expressed circRNAs and miRNAs was generated (Fig. 2). The detailed results including miRNA binding sites are presented in Table SII.
KEGG pathway analysis
Previous studies have shown that circRNAs modulate the expression of their host genes (6,7). Thus, KEGG pathway enrichment analysis was performed using host genes of differentially expressed circRNAs. A total of 41 host genes were enriched in 17 pathways (Fig. 3; Table SIII). B-Raf proto-oncogene, serine/threonine kinase (BRAF) and mitogen-activated protein kinase kinase 2 (MAP2K2) interacted with the highest number of pathways, which suggested that these two genes may serve important roles in PDAC.
Since apoptosis and proliferation are associated with cancer development and progression, the subnetworks of these pathways were extracted from the 17 enriched KEGG pathways. Subsequently, a new network was constructed with these extracted subnetworks, which included 16 host genes that potentially affect cell survival directly or indirectly (Fig. 4). Thus, target circRNAs were obtained for future studies.
Discussion
Two former studies from one research group have profiled the expression of circRNAs in PDAC (14,21). However, only six pairs of PDAC and adjacent normal tissues were profiled and the sample size was relatively small. In the present study, comprehensive analysis with 26 samples from GEO was performed to explore the putative roles of circRNAs in PDAC and to provide novel insights into PDAC biology. The regulatory roles of circRNAs in cancer are attracting increasing attention. Li et al (14) have identified a number of dysregulated circRNAs in PDAC. Another study further investigated the expression level of one dysregulated circRNA termed circ-LDLRAD3, which was upregulated in pancreatic cancer cell lines as well as tumors and plasma from patients with pancreatic cancer, and circ-LDLRAD3 was demonstrated to be associated with invasion and metastasis of pancreatic cancer (22). circRNA exerts its regulatory roles through sequestering cancer-related miRNAs (13,23–26). Therefore, differentially expressed circRNAs and miRNAs in PDAC were used to predict circRNA/miRNA interactions. circRNAs may upregulate the expression of target genes of miRNAs by sequestering the miRNAs. In addition, if circRNAs and their host genes share the same miRNA binding sites, the circRNAs may upregulate the expression of their host genes. In the present study, a total of 51 circRNA/miRNA interactions were identified using 256 differentially expressed circRNAs and 20 differentially expressed miRNAs and a network of the circRNA/miRNA interactions was produced to illustrate the relationships of these genes. The network contained 41 circRNA nodes, 14 miRNA nodes and 51 edges. The number of connections indicates gene involvement in various pathways in PDAC; therefore, the nodes connecting multiple genes may be worth further investigation using experimental methods.
circRNAs are encoded in exons and/or introns of their parental genes (27). Previous studies have shown that circRNAs modulate the expression of their parental genes (6,7,28,29). To explore the potential roles of circRNAs, the predominant pathways of the gene symbols of the differentially expressed circRNAs were identified by a pathway-mapping tool. BRAF and MAP2K2 interacted with the highest number of pathways. BRAF and MAP2K2 are two important members of MAPK signaling pathway, which is associated with PDAC (30,31). The results of the present study indicated that the majority of the inspected genes achieved their effects through the PI3K/AKT signaling pathway. Pancreatic acinar-to-ductal metaplasia is an initiating event that is induced through the PI3K/AKT signaling pathway and can progress to PDAC (32). Another pathway involving cell proliferation is the Wnt/β-catenin signaling pathway. Sox15 exerts its tumor-suppressive effects in pancreatic cancer by suppressing the Wnt/β-catenin signaling pathway (33). Recent data has suggested that the WNT/β-catenin signaling pathway upregulated the expression of cellular communication network factor 1 and consequently enhanced pancreatic cancer development and malignancy (34). In the present study, the aim was to identify circRNA-associated pathways that are involved in the development of PDAC.
As well as acting as miRNA sponges, circRNAs may exert their function through protein binding (35), protein coding (36,37), modulating transcriptional activity of RNA polymerase II (38) and competing with linear splicing (39). Cellular localization has been linked with physiological function and functional mechanism of circRNAs (40,41). Most circRNAs comprise of exonic sequences and are located in the cytoplasm, whereas a small number of intron-derived circRNAs are located in the nucleus (42). Therefore, it is necessary to confirm the cellular localization of a circRNA before studying its roles.
The potential limitations of the present study need to be considered when interpreting the results. Firstly, the results were not verified in clinical samples; differentially expressed circRNAs and miRNAs may be verified using reverse transcription-quantitative PCR in larger clinical samples. Secondly, the circRNA/miRNA interactions were predicted using a bioinformatics method and should be confirmed by the luciferase reporter assay. Thirdly, host gene symbols of circRNAs were used to carry out the KEGG pathway enrichment analysis; whether these circRNAs exert an effect in the enriched pathways should be confirmed using experimental methods.
In conclusion, the potential roles of circRNAs in PDAC were explored using bioinformatics. In total, 256 differentially expressed circRNAs and 20 differentially expressed miRNAs were identified in PDAC tissues compared with normal tissues. The roles of circRNAs in PDAC were explored by prediction of circRNA/miRNA interactions and KEGG pathway enrichment analysis. The circRNA-associated pathways presented in present study could help uncover the effect of circRNAs on PDAC development. Further study confirming the effects of these circRNAs on PDAC development are required in order to improve PDAC diagnosis and therapy.
Supplementary Material
Supporting Data
Acknowledgements
Not applicable
Funding
This study was funded by The National Natural Science Foundation of China (grant no. 81602551) and the young talents program of Jiangsu Cancer Hospital (grant no. 2017YQL-10).
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Authors' contributions
QZ and JYW interpreted the data and wrote the paper. SLZ, SYZ and SJY collected the data and conducted the analysis. SLZ designed this study and revised the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Not applicable.
Patient consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
References
Wachsmann MB, Pop LM and Vitetta ES: Pancreatic ductal adenocarcinoma: A review of immunologic aspects. J Investig Med. 60:643–663. 2012. View Article : Google Scholar : PubMed/NCBI | |
Siegel RL, Miller KD and Jemal A: Cancer statistics, 2019. CA Cancer J Clin. 69:7–34. 2019. View Article : Google Scholar : PubMed/NCBI | |
Stark A and Eibl G: Pancreatic ductal adenocarcinoma. Pancreapedia: Exocrine Pancreas Knowledge Base, Version 1.0. 2015.https://www.pancreapedia.org/reviews/pancreatic-ductal-adenocarcinomaMay 23–2015 | |
He L and Hannon GJ: MicroRNAs: Small RNAs with a big role in gene regulation. Nat Rev Genet. 5:522–531. 2004. View Article : Google Scholar : PubMed/NCBI | |
Hernandez YG and Lucas AL: MicroRNA in pancreatic ductal adenocarcinoma and its precursor lesions. World J Gastrointest Oncol. 8:18–29. 2016. View Article : Google Scholar : PubMed/NCBI | |
Li J, Yang J, Zhou P, Le Y, Zhou C, Wang S, Xu D, Lin HK and Gong Z: Circular RNAs in cancer: Novel insights into origins, properties, functions and implications. Am J Cancer Res. 5:472–480. 2015.PubMed/NCBI | |
Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, Sun W, Dou K and Li H: Circular RNA: A new star of noncoding RNAs. Cancer Lett. 365:141–148. 2015. View Article : Google Scholar : PubMed/NCBI | |
Zhong Z, Lv M and Chen J: Screening differential circular RNA expression profiles reveals the regulatory role of circTCF25-miR-103a-3p/miR-107-CDK6 pathway in bladder carcinoma. Sci Rep. 6:309192016. View Article : Google Scholar : PubMed/NCBI | |
Xia W, Qiu M, Chen R, Wang S, Leng X, Wang J, Xu Y, Hu J, Dong G, Xu PL and Yin R: Circular RNA has_circ_0067934 is upregulated in esophageal squamous cell carcinoma and promoted proliferation. Sci Rep. 6:355762016. View Article : Google Scholar : PubMed/NCBI | |
Nair AA, Niu N, Tang X, Thompson KJ, Wang L, Kocher JP, Subramanian S and Kalari KR: Circular RNAs and their associations with breast cancer subtypes. Oncotarget. 7:80967–80979. 2016. View Article : Google Scholar : PubMed/NCBI | |
Sand M, Bechara FG, Sand D, Gambichler T, Hahn SA, Bromba M, Stockfleth E and Hessam S: Circular RNA expression in basal cell carcinoma. Epigenomics. 8:619–632. 2016. View Article : Google Scholar : PubMed/NCBI | |
Sand M, Bechara FG, Gambichler T, Sand D, Bromba M, Hahn SA, Stockfleth E and Hessam S: Circular RNA expression in cutaneous squamous cell carcinoma. J Dermatol Sci. 83:210–218. 2016. View Article : Google Scholar : PubMed/NCBI | |
Xie H, Ren X, Xin S, Lan X, Lu G, Lin Y, Yang S, Zeng Z, Liao W, Ding YQ and Liang L: Emerging roles of circRNA_001569 targeting miR-145 in the proliferation and invasion of colorectal cancer. Oncotarget. 7:26680–26691. 2016.PubMed/NCBI | |
Li H, Hao X, Wang H, Liu Z, He Y, Pu M, Zhang H, Yu H, Duan J and Qu S: Circular RNA expression profile of pancreatic ductal adenocarcinoma revealed by microarray. Cell Physiol Biochem. 40:1334–1344. 2016. View Article : Google Scholar : PubMed/NCBI | |
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W and Smyth GK: Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43:e472015. View Article : Google Scholar : PubMed/NCBI | |
Benjamini Y, Drai D, Elmer G, Kafkafi N and Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 125:279–284. 2001. View Article : Google Scholar : PubMed/NCBI | |
Yang Y, Gao X, Zhang M, Yan S, Sun C, Xiao F, Huang N, Yang X, Zhao K, Zhou H, et al: Novel role of FBXW7 circular RNA in repressing glioma tumorigenesis. J Natl Cancer Inst. 110:2018. View Article : Google Scholar | |
Zhong S, Wang J, Zhang Q, Xu H and Feng J: CircPrimer: A software for annotating circRNAs and determining the specificity of circRNA primers. BMC Bioinformatics. 19:2922018. View Article : Google Scholar : PubMed/NCBI | |
Enright AJ, John B, Gaul U, Tuschl T, Sander C and Marks DS: MicroRNA targets in Drosophila. Genome Biol. 5:R12003. View Article : Google Scholar : PubMed/NCBI | |
Krüger J and Rehmsmeier M: RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res 34 (Web Server Issue). W451–W454. 2006. View Article : Google Scholar | |
Qu S, Song W, Yang X, Wang J, Zhang R, Zhang Z, Zhang H and Li H: Microarray expression profile of circular RNAs in human pancreatic ductal adenocarcinoma. Genom Data. 5:385–387. 2015. View Article : Google Scholar : PubMed/NCBI | |
Yang F, Liu DY, Guo JT, Ge N, Zhu P, Liu X, Wang S, Wang GX and Sun SY: Circular RNA circ-LDLRAD3 as a biomarker in diagnosis of pancreatic cancer. World J Gastroenterol. 23:8345–8354. 2017. View Article : Google Scholar : PubMed/NCBI | |
Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, Luo Y, Lyu D, Li Y, Shi G, et al: Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 7:112152016. View Article : Google Scholar : PubMed/NCBI | |
Xu H, Guo S, Li W and Yu P: The circular RNA Cdr1as, via miR-7 and its targets, regulates insulin transcription and secretion in islet cells. Sci Rep. 5:124532015. View Article : Google Scholar : PubMed/NCBI | |
Nan A, Chen L, Zhang N, Liu Z, Yang T, Wang Z, Yang C and Jiang Y: A novel regulatory network among LncRpa, CircRar1, MiR-671 and apoptotic genes promotes lead-induced neuronal cell apoptosis. Arch Toxicol. 91:1671–1684. 2017. View Article : Google Scholar : PubMed/NCBI | |
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK and Kjems J: Natural RNA circles function as efficient microRNA sponges. Nature. 495:384–388. 2013. View Article : Google Scholar : PubMed/NCBI | |
Chen LL and Yang L: Regulation of circRNA biogenesis. RNA Biol. 12:381–388. 2015. View Article : Google Scholar : PubMed/NCBI | |
He J, Xie Q, Xu H, Li J and Li Y: Circular RNAs and cancer. Cancer Lett. 396:138–144. 2017. View Article : Google Scholar : PubMed/NCBI | |
Li Z, Huang C, Bao C, Chen L, Lin M, Wang X, Zhong G, Yu B, Hu W, Dai L, et al: Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 22:256–264. 2015. View Article : Google Scholar : PubMed/NCBI | |
Papademetrio DL, Lompardía SL, Simunovich T, Costantino S, Mihalez CY, Cavaliere V and Álvarez É: Inhibition of survival pathways MAPK and NF-κB triggers apoptosis in pancreatic ductal adenocarcinoma cells via suppression of autophagy. Target Oncol. 11:183–195. 2016. View Article : Google Scholar : PubMed/NCBI | |
Herner A, Sauliunaite D, Michalski CW, Erkan M, De Oliveira T, Abiatari I, Kong B, Esposito I, Friess H and Kleeff J: Glutamate increases pancreatic cancer cell invasion and migration via AMPA receptor activation and Kras-MAPK signaling. Int J Cancer. 129:2349–2359. 2011. View Article : Google Scholar : PubMed/NCBI | |
Chen Y, Bai X, Zhang Q, Wen L, Su W, Fu Q, Sun X, Lou Y, Yang J, Zhang J, et al: The hepatitis B virus X protein promotes pancreatic cancer through modulation of the PI3K/AKT signaling pathway. Cancer Lett. 380:98–105. 2016. View Article : Google Scholar : PubMed/NCBI | |
Thu KL, Radulovich N, Becker-Santos DD, Pikor LA, Pusic A, Lockwood WW, Lam WL and Tsao MS: SOX15 is a candidate tumor suppressor in pancreatic cancer with a potential role in Wnt/β-catenin signaling. Oncogene. 33:279–288. 2014. View Article : Google Scholar : PubMed/NCBI | |
Sano M, Driscoll DR, DeJesus-Monge WE, Quattrochi B, Appleman VA, Ou J, Zhu LJ, Yoshida N, Yamazaki S, Takayama T, et al: Activation of WNT/β-catenin signaling enhances pancreatic cancer development and the malignant potential Via Up-regulation of Cyr61. Neoplasia. 785–794. 2016. View Article : Google Scholar : PubMed/NCBI | |
Du WW, Yang W, Liu E, Yang Z, Dhaliwal P and Yang BB: Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 44:2846–2858. 2016. View Article : Google Scholar : PubMed/NCBI | |
Legnini I, Di Timoteo G, Rossi F, Morlando M, Briganti F, Sthandier O, Fatica A, Santini T, Andronache A, Wade M, et al: Circ-ZNF609 Is a circular RNA that can be translated and functions in myogenesis. Mol Cell. 66:22–37.e9. 2017. View Article : Google Scholar : PubMed/NCBI | |
Yang Y, Fan X, Mao M, Song X, Wu P, Zhang Y, Jin Y, Yang Y, Chen LL, Wang Y, et al: Extensive translation of circular RNAs driven by N6-methyladenosine. Cell Res. 27:626–641. 2017. View Article : Google Scholar : PubMed/NCBI | |
Zhang Y, Zhang XO, Chen T, Xiang JF, Yin QF, Xing YH, Zhu S, Yang L and Chen LL: Circular intronic long noncoding RNAs. Mol Cell. 51:792–806. 2013. View Article : Google Scholar : PubMed/NCBI | |
Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, Evantal N, Memczak S, Rajewsky N and Kadener S: circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. 56:55–66. 2014. View Article : Google Scholar : PubMed/NCBI | |
Chen L and Shan G: Circular RNAs remain peculiarly unclear in biogenesis and function. Sci China Life Sci. 58:616–618. 2015. View Article : Google Scholar : PubMed/NCBI | |
Chen L, Huang C, Wang X and Shan G: Circular RNAs in eukaryotic cells. Curr Genomics. 16:312–318. 2015. View Article : Google Scholar : PubMed/NCBI | |
Chen S and Zhao Y: Circular RNAs: Characteristics, function, and role in human cancer. Histol Histopathol. 33:887–893. 2018.PubMed/NCBI |