Introduction {#Sec1}
============

Circular RNAs comprise a new class of long noncoding RNAs characterized by their 5′ and 3′ ends covalently joined. They were misinterpreted as splicing errors for more than 20 years until their rediscovery in 2012 as diverse, highly abundant, conserved and naturally occurring RNAs in eukaryotes^[@CR1]--[@CR5]^.

About 90,000 different circular RNAs were described in human, which most are derived mainly from annotated exons (\~85%) and a smaller fraction from untranslated regions (UTRs), introns and unannotated regions of the genome. They are most commonly formed from two or three exons, comprising between a hundred and four thousand nucleotides in length^[@CR1],[@CR2],[@CR5]--[@CR7]^.

These RNA molecules are likely generated by a process known as back-splicing. This noncanonical splicing can produce three types of circular RNAs, in which they are classified: exonic circular RNAs (circRNAs), circular intronic RNAs (ciRNAs) and exon-intron circular RNAs (EIciRNAs)^[@CR1],[@CR8]--[@CR11]^. CircRNAs are predominantly cytoplasmic and were reported acting as microRNAs (miRNAs) and RNA-binding proteins (RBPs) sponges. CiRNAs and EIciRNAs are enriched in the nucleus and are RNA polymerase II-associated, suggesting that they promote the transcription of their parent genes^[@CR2],[@CR5],[@CR8],[@CR12],[@CR13]^.

Circular RNAs molecules are easily accessed and measured in body fluids and have distinct characteristics such as tissue-specificity and stability in both intra and extracellular environments. This suggest their potential as clinical markers that may provide new insights into the prevention and treatment of several diseases^[@CR14]^.

Although neither their biogenesis nor roles have been entirely understood, circular RNA expression has already been described as altered in human diseases such diabetes, atherosclerosis, Alzheimer's disease and cancer^[@CR15]--[@CR17]^. On cancer, they were associated with cellular proliferation, and some clinical features such as tumor size and presence of distal metastases^[@CR14],[@CR18]--[@CR21]^.

Among different cancers, gastric cancer remains the third leading cause of cancer-related death worldwide. Due the lack of specific symptoms, most gastric cancer patients are diagnosed in advanced-stage disease with a poor prognosis^[@CR22]^. Some reports have shown that recurrence of gastric cancer may be due the field cancerization (or field effect) in gastric mucosa. According to this theory, the tissue surrounding tumors, despite being histologically normal, shares molecular abnormalities that are present in fully developed tumors^[@CR23]--[@CR25]^. Multiple genetic and epigenetic alterations, mostly DNA methylation and miRNA abnormal expression, have been described as field effect biomarkers in gastric cancer, reinforcing the occurrence of a field effect in this type of cancer^[@CR26],[@CR27]^.

MiRNAs are a class of small nonconding RNAs involved in many biological processes by blocking target mRNAs translation^[@CR28]^. The epigenetic network in which the miRNAs participate is complex and dynamic since its involves not only target mRNAs, but also other types of noncoding RNAs such as the circular RNAs^[@CR14]^. Given that some circular RNAs act as miRNAs sponges, they may also have a potential epigenetic regulation role in gastric cancer.

The aim of this study was to identify, characterize and compare the entirety of all expressed circular RNAs in samples of patients without gastric cancer, gastric cancer samples and matched tumor-adjacent gastric tissue. Additionally, we correlated circular RNAs' expression data with miRNA expression.

Results {#Sec2}
=======

We performed RNA-Seq on ribosomal-depleted total RNA isolated from gastric tissue samples. Head-to-tail back-spliced junctions were detected by using two combined prediction algorithms (Supplementary Fig. [1](#MOESM1){ref-type="media"}).

In total, we detected 736 unique annotated circular RNAs in all three groups of gastric tissues. As shown in Fig. [1a](#Fig1){ref-type="fig"}, we identified 66 annotated circular RNAs in gastric tissue without gastric cancer, 620 in matched tumor-adjacent gastric tissue and 220 in gastric cancer samples.Figure 1Total of annotated circular RNAs detected in gastric tissue. (**a**) Number of expressed circular RNAs in each type of gastric tissue according to their origin. (**b**) Venn diagram of all expressed circular RNAs between the three types of gastric tissue. CDS: coding DNA sequence.

A previous study showed that most of human circular RNAs contain two or three exons^[@CR29]^. To further evaluate this data, we analyzed the number of exons per circular RNA in gastric tissue and found similar results (Supplementary Fig. [2](#MOESM1){ref-type="media"}). As shown in Table [1](#Tab1){ref-type="table"}, the number of exons is not necessarily related to the circular RNAs spliced lengths. A notable example is that hsa_circ_0004176, which harbors 26,767 nt in length, spans only two exons, while hsa_circ_0020397, which harbors 2,738 nt in length, spans 26 exons.Table 1Transcript features of the expressed circular RNAs in gastric tissues.Attribute/Type of gastric tissueWithout cancerTumor-adjacentGastric cancerLongest circular RNA\*hsa_circ_0000246 (*MCU*)\*hsa_circ_0000230 (*ZEB1*)hsa_circ_0004176 (*IFT43*)7620 nt-3 exons88220 nt-4 exons26767 nt-2 exonsShortest circular RNAhsa_circ_0055734 (*ANKRD36*)hsa_circ_0000439 (*ATXNK2*)hsa_circ_0000439 (*ATXNK2*)98 nt spliced97 nt97 ntCircular RNA with the highest number of exonshsa_circ_0020397 (*DOCK1*)hsa_circ_0001613 (*SENP6*)hsa_circ_0023923 (*PICALM*)26 exons-2738 nt12 exons-1722 nt11 exons-1128 ntGene with the highest number of expressed circular RNA isoforms*XPO1UBAP2UBAP2*hsa_circ_0001017 (307 nt-3 exons)hsa_circ_0001849 (119 nt -2 exons)hsa_circ_0001849 (119 nt -2 exons)hsa_circ_0001016 (2 exons)hsa_circ_0001851 (159 nt-2 exons)hsa_circ_0001851 (159 nt-2 exons)*UBAP2*hsa_circ_0001847 (377 nt-4 exons)hsa_circ_0001847 (377 nt-4 exons)hsa_circ_0001849 (119 nt -2 exons)hsa_circ_0005993 (187 nt-3 exons)hsa_circ_0001850 (278 nt-4 exons)hsa_circ_0001851 (159 nt 2 exons)\*Exon-intron circular RNA.

Interestingly, *UBAP2* gene presented five different circular RNA isoforms expressed in gastric tissue, suggesting that circular alternative splicing is also occurring in the stomach (Table [1](#Tab1){ref-type="table"}).

To examine the genomic localization of gastric circular RNAs, we analyzed the number of circular RNAs per chromosome, and found that most of them is derived from chromosome 1 of the human genome (Supplementary Fig. [3](#MOESM1){ref-type="media"}).

Although most gastric circular RNAs had less than 10 back-spliced junction reads of coverage, some highly expressed circular RNAs in matched adjacent gastric tissue had a read count of more than 35. Table [2](#Tab2){ref-type="table"} shows the most expressed circular RNAs in gastric tissue without gastric cancer, matched tumor-adjacent gastric tissue and gastric cancer samples.Table 2List of the most expressed circular RNAs in gastric tissue.Gastric tissueCircular RNANumber of back spliced junction readsGene symbolWithout cancerhsa_circ_00013406*TMCC1*hsa_circ_00063545*VAMP3*hsa_circ_00013803*UBXN7*hsa_circ_00004193*RAB3IP*hsa_circ_00024963*APPBP2*hsa_circ_00043683*REPS1*hsa_circ_00011123*DGKD*Tumor-adjacenthsa_circ_0001092\*41*CFLAR*hsa_circ_0000284\*37*HIPK3*hsa_circ_0000437\*37*CORO1C*hsa_circ_0001136\*31*ASXL1*hsa_circ_000144528*SMARCA5*hsa_circ_0000211\*26*SFMBT2*hsa_circ_000149821*WDR41*Gastric cancerhsa_circ_0000437\*17*CORO1C*hsa_circ_0001136\*15*ASXL1*hsa_circ_0000211\*14*SFMBT2*hsa_circ_0000284\*13*HIPK3*hsa_circ_0001092\*11*CFLAR*hsa_circ_000189711*POMT1*hsa_circ_00017278*ZKSCAN1*\*Circular RNAs in common between tumor-adjacent and gastric cancer tissues.

To further explore the potential function of the expressed circular RNAs in gastric tissue, we selected the gastric circular RNAs-derived genes to perform GO enrichment analysis (Fig. [2](#Fig2){ref-type="fig"}). The gastric tissue without gastric cancer and matched tumor-adjacent gastric tissue circular RNAs-derived genes were enriched in the process of bacterial invasion of epithelial cells, such as *Salmonella* sp., *Listeria* sp. and *Shigella* sp^[@CR30]^ (Supplementary Fig. [4](#MOESM1){ref-type="media"}). Tumor-adjacent gastric tissue circular RNAs-derived genes also were enriched in cancer-related processes, as well as gastric cancer's.Figure 2GO enrichment of the gastric circular RNAs-derived genes, evidencing the KEGG pathways and its scores. (**A**) Gastric tissue without gastric cancer. (**B**) Matched tumor-adjacent gastric tissue. (**C**) Gastric cancer.

Circular RNAs can regulate miRNAs by sequestering them by binding to their seed sequences^[@CR2],[@CR5]^. Given that, we identified candidate target miRNAs of the most expressed circular RNAs in gastric tissues. We realized that the seed sequence is the key that may link circRNAs, miRNAs, miRNAs target genes and circular RNAs-derived genes. Therefore, we searched for the candidate target miRNAs by identifying the miRNAs that regulates such circular RNA-derived gene and by confirming that the complementary seed sequence is present in the circRNA sequence.

After this analysis, to consolidate the candidate target miRNAs, we compared them with the differentially expressed miRNAs identified in the same samples of this study, which were obtained previously by RNA-Seq by our group \[data not published\]. We found five candidate miRNAs potentially regulated by five circRNAs. All of them were previously described in gastric cancer (Table [3](#Tab3){ref-type="table"}). In Fig. [3](#Fig3){ref-type="fig"}, we illustrated the interaction between *CORO1C*, hsa_circ_0000437 and *hsa-miR-1*.Table 3Candidate target microRNAs of some of the high expressed circular RNAs in gastric tissue.Circular RNAsTarget microRNAsNameGeneNameNumber of seed matchesRefs in gastric cancerhsa_circ_0001340*TMCC1hsa-miR-452--5p*1\[30,31\]hsa_circ_0000419*RAB3IPhsa-miR-145--5p*1\[32--34\]hsa_circ_0001112*DGKDhsa-miR-375*3\[35--37\]hsa_circ_0000284*HIPK3hsa-miR-224--5p*1\[31,38\]hsa_circ_0000437*CORO1Chsa-miR-1*1\[39\] Figure 3Simulation of the relation between *CORO1C*, hsa_circ_0000437 and *hsa-miR-1*. Pol II: RNA polymerase II.

Unlike CDR1as, some studies have demonstrated that most circRNAs would have only 1--2 miRNA binding sites^[@CR13],[@CR31]^. Our data corroborate to these studies given that most of the circRNA identified have only one miRNA-binding site, except for hsa_circ_0001112 that have three binding sites (Table [3](#Tab3){ref-type="table"}).

We also analyzed the distribution of the expressed circular RNAs in gastric tissue without cancer, matched tumor-adjacent gastric tissue and gastric cancer samples. The Fig. [1b](#Fig1){ref-type="fig"} shows that there are exclusive circular RNAs of each group, but also there are common circular RNAs between them. Differential expression analysis showed that of the 27 circular RNAs in common between the three groups, five are significantly different (Table [4](#Tab4){ref-type="table"}).Table 4List describing the five differentially expressed circular RNAs in gastric tissue. The differential expression was evaluated with negative binomial regression adjusting for common and tagwise variation, and p-values were adjusted for multiple testing using a FDR procedure.Circular RNAGeneGene official nameP-valueDriver gene\*hsa_circ_0001136*ASXL1*Additional sex combs like 18,3E-04Yeshsa_circ_0000284*HIPK3*Homeodomain interacting protein kinase 39,0E-04Nohsa_circ_0000211*SFMBT2*Scm-like with four mbt domains 29,8E-04Nohsa_circ_0004771*NRIP1*Nuclear receptor interacting protein 14,6E-05Nohsa_circ_0000524*RBM23*RNA binding motif protein 231,9E-04No\*According to Vogelstein *et al*.^[@CR41]^.

The differential expression analysis was performed by comparing the samples without cancer with both tumor-adjacent and gastric cancer samples combined. All five differentially expressed circular RNAs are exonic, and were found down regulated in samples without cancer (Fig. [4](#Fig4){ref-type="fig"}).Figure 4Expression of the five differentially expressed circular RNAs in gastric tissue. This analysis was performed by comparing the samples without cancer with both tumor-adjacent and gastric cancer samples combined.

Discussion {#Sec3}
==========

Circular RNAs are a novel class of regulatory noncoding RNAs with yet unknown impact on the cellular machinery. Our study is the first to investigate and describe all circular RNAs expressed in adult human gastric tissue, comprising patients without gastric cancer, matched tumor-adjacent gastric tissue and gastric cancer samples.

We found that the matched tumor-adjacent gastric samples were the group with the highest number of circular RNAs identified, followed by gastric cancer and samples of patients without gastric cancer (Fig. [1a](#Fig1){ref-type="fig"}). Most of the previous studies about circular RNAs global expression in human cancers used only the matched tumor-adjacent samples as normal control. In all these studies, the expression of circular RNAs in cancer is down-regulated in comparison to the matched tumor-adjacent tissue^[@CR13],[@CR18],[@CR19],[@CR32]^. These data suggest that the abundant expression of circular RNAs in tumor-adjacent tissue samples is a general pattern in several types of cancer, including gastric cancer.

Circular RNAs expressions were analyzed in gastric cancer in some previous studies. However, these studies used matched tumor-adjacent as control^[@CR13],[@CR33]--[@CR38]^. The use of adjacent tissue for comparison purposes can lead to biases since the evidences have demonstrated the field cancerization in gastric tissue surrounding the tumors^[@CR26],[@CR27]^. Thus, we chose to investigate the circular RNAs expression in patients without gastric cancer, matched tumor-adjacent gastric and gastric cancer samples.

Our data suggests that circular RNAs abundance in tumor-adjacent tissue may be somehow related to gastric carcinogenesis, given its similarity to gastric cancer tissue. Most of the highest expressed circRNA genes in gastric cancer samples are also present in tumor-adjacent tissue (*CFLAR*, *CORO1C*, *HIPK3*, *ASXL1* and *SFMBT2*) (Table [2](#Tab2){ref-type="table"}).

It is possible that the circular RNAs are not essential molecules in fully developed tumors, explaining their high expression in tumor-adjacent tissues. Bachmayr-Heyda *et al*.^[@CR18]^ showed that the expression of circular RNAs in colorectal cancer cell lines is even smaller than those in colorectal cancer tissue. The cancer cell lines have a higher proliferation rate and are pure cancer cells, indicating that cancerous cells do not require a high level of circular RNAs to maintain their malignant features.

Although most circular RNAs does not have its function completely understood, it is possible to estimate their cellular role by performing a functional enrichment analysis of their derived genes. GO enrichment indicated that the gastric tissue without gastric cancer circular RNAs-derived genes were enriched for the process of bacterial invasion of epithelial cells, which is a natural process in stomach (Fig. [2](#Fig2){ref-type="fig"}). This KEGG pathway was also enriched in tumor-adjacent samples, but not in gastric cancer samples, indicating the cellular loss of function typically found in cancer.

Previous studies have discussed the potential function of circRNAs as miRNA sponges. Memczak *et al*.^[@CR5]^ reported that the circRNA CDR1as (or ciRS-7) harbors about 70 binding sites for miR-7 seed. However, a deeper analysis showed that most circRNAs have less than 10 miRNA binding sites, indicating that miRNA sponging by circRNAs may not require a large number of binding sites^[@CR31]^. To further investigate this information, we identified the potential circRNAs target and found five candidate miRNAs, and most of them present only one target site (Table [3](#Tab3){ref-type="table"}).

All five candidate miRNAs were found differentially expressed between patients without gastric cancer, matched tumor-adjacent gastric and gastric cancer samples (data not shown), and previously described in association with gastric cancer in the literature. Their expressions were correlated with several features of gastric cancer, such as drug resistance, proliferation, invasion, migration and cell growth in gastric cancer^[@CR39]--[@CR46]^.

The Fig. [3](#Fig3){ref-type="fig"} illustrates how complex and dynamic is the interaction between circRNA, mRNA and circRNA-derived gene. *CORO1C* gene produces circRNA and mRNA by noncanonical and canonical splicing, respectively, and both types of RNA may interact with the same miRNA. Circ-CORO1C blocks *hsa-miR-1*, while CORO1C mRNA is blocked by *hsa-miR-1*. It suggests that the circRNA production may be a gene mechanism to ensure its own mRNA translation.

Given that, regarding the type (circRNA, ciRNA or EIciRNA), circular RNAs seem to be a positive self-mechanism of gene regulation by sponging miRNAs or by interacting with RNA polymerase II.

To identify circular RNAs with potential to become gastric cancer biomarkers, we performed differential expression analysis. Among the five differentially expressed circRNAs, hsa_circ_0001136 is derived from *ASXL1*, which is a driver gene involved in chromatin modelling^[@CR47]^ (Table [4](#Tab4){ref-type="table"}).

Hsa_circ_0000284 (*HIPK3* gene) was found differentially expressed in gastric tissues (Table [4](#Tab4){ref-type="table"}). Given that this circRNA is overexpressed in tumor-adjacent and gastric cancer samples, and also may regulate *hsa-miR-224--5p* (Table [3](#Tab3){ref-type="table"}), the interaction between hsa_circ_0000284 and *hsa-miR-224* is possibly involved in gastric carcinogenesis. In fact, this circRNA was found overexpressed in seven types of cancer, including gastric, and related to cell proliferation^[@CR13]^. *Hsa-miR-224* was also described in association to gastric cancer^[@CR45]^.

Circular RNAs have some particularities that made them potential biomarkers of both physiological and pathological processes. Besides being abundant, stable and resistant, their little invasiveness remarkably increases its potential, since their expression can be accessed by body fluids^[@CR14]^. Shao *et al*.^[@CR34]^ demonstrated that the expression of circular RNA can be accessed by gastric juice, suggesting their potential as biomarker for disease screening.

Overall, our results revealed that the circular RNAs is overexpressed in tumor-adjacent and in gastric cancer samples in comparison to samples without cancer. We showed the presence of field cancerization in gastric cancer, indicating that the tumor-adjacent tissue cannot be considered as normal tissue. We also found five differentially expressed circRNAs that may become novel biomarkers of gastric cancer and need to be further validated. Nevertheless, our results support the hypothesis of circular RNAs representing a novel factor in the dynamic epigenetic network of gene regulation, which involves the miRNAs and its mRNAs targets and the circular RNAs-derived genes. Further studies are needed to elucidate the roles and the functional relevance of the circular RNAs in human diseases.

Methods {#Sec4}
=======

Clinical samples {#Sec5}
----------------

We included tissue samples of patients without gastric cancer (n = 8), gastric cancer (n = 8) and matched tumor-adjacent (n = 8), from the Universitary Hospital of João de Barros Barreto of the Federal University of Pará. All samples were collected, stored in *RNAlater* (Thermo Fisher Scientific) and frozen in liquid nitrogen until RNA total isolation. The study including all experimental protocols was approved by the *Ethics Committee of the Center of Oncology Research of the Federal University of Pará* (No. 1.432.512). All study participants or their legal guardian provided informed written consent in accordance with the Helsinki Declaration. The methods were performed in accordance to the approved guidelines.

RNA isolation {#Sec6}
-------------

Total RNA was isolated from tissue samples by using TRIzol Reagent (Thermo Fisher Scientific) following the manufacture's protocol. Total RNA integrity and amount were evaluated by Qubit 2.0 Fluorometer (Thermo Fisher Scientific), NanoDrop ND-1000 (Thermo Fisher Scientific) and 2200 Tape Station System (Agilent). The integrity criteria were values between 1.8 and 2.2 (A ~260/280~), \>1.8 (A~260/230~), and RIN ≥ 5.

Circle-Seq sample treatment, library synthesis, sequencing and analysis {#Sec7}
-----------------------------------------------------------------------

First, a step of circular RNA enrichment was made by treating the total RNA with 3U of RNase R (Epicentre), followed by 15 minutes at 37 °C. After this, the treated RNA was re-quantified, and 1 μg of treated RNA per sample was used as input to prepare the libraries. We synthesized 24 libraries by using TruSeq Stranded Total RNA Library Prep with Ribo-Zero Gold (Illumina), which already has a step of rRNA depletion included. The libraries quality was controlled with 2200 TapeStation (Agilent), normalized to 10 nM and sequenced on a MiSeq Sequencing System (Illumina) by using the MiSeq Reagent Kit v3 (Illumina).

FASTQ was trimmed, cropped and adapters contaminant were removed (Trimmomatic v.0.36). The resulting reads were aligned to human genome (hg19) using both BWA (v.0.7) and STAR (v.2.5), which were processed by CIRI (v.2.0)^[@CR48]^ and CIRCexplorer2 (v.2.2)^[@CR49]^, respectively, to detect head-to-tail back-spliced junctions. We considered only the junctions detected by both tools to improve prediction accuracy^[@CR50]^.

The detected circRNA list was used to made a Venn diagram (Venny 2.1 - <http://bioinfogp.cnb.csic.es/tools/venny/index.html>) representing the distribution of the expressed circular RNAs among gastric tissue without gastric cancer, gastric cancer samples and matched tumor-adjacent samples. All other graphics and statistical analyses were performed by using R (v.3.3). The read count was normalized and compared between groups using edgeR (v.3.18) package (REF).

Circular RNAs functional analysis {#Sec8}
---------------------------------

Gastric circular RNAs-derived genes were selected to perform for the functional enrichment analysis. This analysis was performed by DAVID Bioinformatics Resources v6.8 (<https://david.ncifcrf.gov>). All enriched KEGG pathways were plotted. P-values were adjusted by using Bonferroni's correction.

Selection of the candidate target microRNAs {#Sec9}
-------------------------------------------

The candidate target miRNAs were predicted by searching which miRNA has the circular RNA-derived gene as a target. This search was performed by using the miRTarBase, an experimentally validated microRNA-target interactions database (<http://mirtarbase.mbc.nctu.edu.tw>). After that, we searched for a complementary region to miRNA seed sequence in circular RNA, and confirmed that the predicted miRNA was found differentially expressed in gastric cancer \[data not published\].

Electronic supplementary material
=================================

 {#Sec10}

Supplementary material

Amanda Ferreira Vidal and André M. Ribeiro-dos-Santos contributed equally to this work.

**Electronic supplementary material**

**Supplementary information** accompanies this paper at 10.1038/s41598-017-15061-w.

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

We thank Fabiano Moreira and Igor Schneider for the helpful suggestions. This work was financially supported by the Rede de Pesquisa em Genômica Populacional Humana (Biocomputacional-Protocol No. 3381/2013/CAPES), PROPESP/UFPA-FADESP; and CNPq (MCTI/CNPq/Universal 14/2014 - No. 457672/2014-5). Ândrea Ribeiro-dos-Santos is supported by CNPq/Produtividade (304413/2015-1). Sidney Santos is supported by CNPq/Produtividade (CNPq 305258/2013-3).

P.P.A., S.E.B.S. and Â. Ribeiro-dos-Santos designed and coordinated the study. S.D. and A.K.M.A collected and prepared the gastric samples. A.F.V., T.V., L.M. and P.P. performed the experimental analyses. A.M. Ribeiro-dos-Santos performed the bioinformatical and statistical analyses. A.F.V. drafted the manuscript, and designed the study. All authors read and approved the final manuscript.

Competing Interests {#FPar1}
===================

The authors declare that they have no competing interests.
