- 1Department of Clinical Laboratory, Qilu Hospital of Shandong University, Jinan, China
- 2Department of Laboratory Medicine, Weifang People’s Hospital, Weifang, China
Objective: At present, there are still no effective diagnosis methods for gastric cancer (GC). Increasing evidences indicate that Extracellular Vesicle circular RNAs (EV circRNAs) play a crucial role in several diseases. However, their correlations with GC are not clarified. This study aims to investigate the expression profile of serum EV circRNAs in GC and evaluate its potential clinical value.
Methods: High-throughput RNA sequencing (RNA-seq) was used to assess circRNA expression profiles between 4 patients with GC and 4 healthy controls. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were employed to determine the biological functions of differentially expressed (DE) circRNAs. A circRNA-miRNA-mRNA network was constructed using bioinformatics tools. Reverse transcription-quantitative polymerase chain reaction (RT-q)PCR was used to validate the dysregulated circRNAs. Receiver operating characteristic (ROC) curves were used to evaluate the diagnostic value of circRNAs for GC.
Results: A total of 4692 circRNAs were detected in the serum EVs of healthy controls and patients with GC, most of which were novel (98%) and intergenic (52%). 7 circRNAs were upregulated and 4 circRNAs were downregulated (|log2Fold Change| > 2, P < 0.05). GO and KEGG pathway enrichment analyses revealed that DE circRNAs were primarily involved in glutathione metabolism, protein folding, and drug metabolism-cytochrome P450. Of these, 3 circRNAs (Chr10q11, Chr1p11, and Chr7q11) were identified to be significantly overexpressed in patients with GC compared with healthy controls using RT-qPCR. The combination of 3 EV circRNAs and carcinoembryonic antigen (CEA) produced an area under the curve (AUC) of 0.866 (95%CI: 0.803-0.915) with a sensitivity and specificity of 80.4% and 81.8%, respectively. Additionally, the expression levels of 3 EV circRNAs were significantly correlated with tumor size, lymph node metastasis, and TNM stage. The circRNA-miRNA-mRNA network showed that the 3 identified circRNAs were predicted to interact with 13 miRNAs and 91 mRNAs.
Conclusion: Our results illustrate that the panel of EV circRNAs in serum are aberrantly expressed and may act as the suitable biomarkers for gastric cancer.
Introduction
Gastric cancer (GC) is the fifth most common type of cancer and the third most common cause of cancer-related death globally. Global cancer statistics show more than one million new cases are diagnosed every year (1–3). Due to the lack of effective markers for GC, the 5-year survival rate is more than 90% in the early stages. Whereas it is only 5-20% in patients with advanced stage, and the median overall survival (OS) is 10 months (1, 4, 5). Therefore, early detection of GC is particularly important for improving the survival rates (6, 7). At present, the diagnostic methods of GC include endoscopic ultrasonography, tissue biopsy, imaging examination, and screening of tumor markers and so on (8). But these methods still have defects to some extent. For example, tissue biopsy is invasive, complicated, and expensive. Improper preparation of samples can lead to misdiagnosis or missed diagnosis (9). X-ray and CT scans have difficulty in detecting lesions less than 1cm (10). The traditional tumor markers for GC detection include carcinoembryonic antigen (CEA), carbohydrate antigen 19−9 (CA19−9), and carbohydrate antigen 72−4 (CA72−4), but they are insufficient (11). For example, our previous study has showed that the expression of CA19-9 and CA72-4 do not differ between groups (12). Therefore, identifying novel markers for GC remains a challenge.
Circular RNAs (circRNAs) are a unique class of non-coding RNAs (ncRNAs) that form a covalently closed loop-structure without 5’ and 3’ ends, making it resistant to RNase R digestion (13, 14).. CircRNAs can be divided into exon, intron and exon-intron according to its origin. It is abundant, stable, highly conserved and tissue-specific (15). As circRNAs do not typically encode protein and its function is not clear, they are initially considered to be the by-products of post-transcriptional modification of mRNA (16). Recently, circRNAs can be detected in body fluid, serum, plasma, and tissue samples, and play an important regulatory role in the development of most diseases (17). Therefore, they have the potential to become promising biomarkers. Chi Hin Wong et al. found that circRTN4 was significantly upregulated in primary tumors from pancreatic ductal adenocarcinoma patients (18). The levels of plasma circFARSA in lung cancer patients was higher than that in healthy controls, suggesting that circFARSA may be used as a potential biomarker for non-invasive diagnosis (19). Three circRNAs (hsa_circ_0000370, hsa_circ_0082812 and hsa_circ_0035445) were abnormally expressed in the plasma of colorectal cancer patients, with an area under the curve (AUC) of 0.815, 0.737 and 0.703, respectively (20). These results suggest that circRNAs can be served as a non-invasive biomarker for cancer diagnosis. Zhu Jin et al. found that the absence of circ-BIRC6 abrogated the progression of non-small cell lung cancer by targeting miR-4491 (21). A study from Zong et al. (22) showed that the expression of circRNA_102231 was significantly upregulated in lung adenocarcinoma and was related to lymph node metastasis. Chen et al. (23) found that patients with lower expression levels of circRNA cRAPGEF5 had a poorer prognosis and inhibit the proliferation and migration of renal cell carcinoma through miR-27a-3p/TXNIP pathway. Interesting, circRNAs also play an indispensable role in GC. A study showed that the upregulated expression of circPVT1 in GC tissues competitively bound to miR-125, promoting the proliferation of cancer cells (24). Zang et al. found that the downregulated expression of circEIF4G3 was observed in GC and was associated with poor clinical outcomes (25). Another study showed that circMAN1A2 could regulate MTA2 through sponging miR-1236-3p, promoting the progression of H. pylori-induced GC (26). These above studies give us a hint that the value of circRNAs in GC remains to be further clarified.
Accumulating studies have shown that extracellular vesicles (EVs) are important carriers of circRNAs. The protective effects of EV membranes and the unique structure of circRNAs extend the half-life of circRNAs and increase RNase R resistance compared to linear RNAs (27, 28). A study showed that there were significant differences between hsa-circRNA-0005795 and hsa-circRNA-0088088 in serum exosomes (29). Chen et al. (30) found that serum exocrine circPRMT5 was increased in bladder cancer, which was closely related to stage and prognosis, and a MiR-30c/SNAIL1/E-cadherin pathway was involved in the process of tumor EMT. Some studies have explored its value for biomarkers. For example, Lu et al. (31) observed that the high expression of circ-RanGAP1 was closely related to stage, metastasis, and prognosis. EV circRNAs have provided a novel perspective for GC diagnosis. However, the research on EV circRNAs as tumor marker in GC are considerably limited.
There are difficulties in detecting and identifying novel EVs circRNAs. In the past decade, RNA-seq has become an indispensable method for transcriptome analysis (32, 33). And now it is most commonly method in analyzing differentially expressed (DE) genes. Zhu et al. (34) performed high-throughput RNA-seq and identified 84 circular RNAs, 41 miRNAs and 398 mRNAs in CRC. Long et al. (35) found that the expression of hsa_circ_0007694 in thyroid papillary carcinoma was downregulated through high-throughput sequencing. Further study found that PI3K/AKT/mTOR and Wnt signaling pathway were involved. In fact, RNA-seq technology can be used to identify novel biomarkers for improving tumor diagnosis. However, there are few studies on EV circRNAs in GC.
In the current study, RNA-seq technology was used to determine the expression profiles of circRNAs in serum EV between patients with GC and healthy controls. GO and KEGG pathway analyses indicated that circRNAs may serve as potential therapeutic targets in GC. CeRNA networks showed a potential role of the circRNAs in cancer progression. The expression and diagnostic value of 3 candidate circRNAs were further evaluated. In addition, the association of these 3 circRNAs with clinicopathological variables was assessed.
Materials and methods
Patients and control subjects
A total of 260 serum samples were collected from Qilu Hospital of Shandong University from December 2016 to December 2017. Of these, serum samples of 4 patients with GC and 4 healthy controls were used for high-throughput sequencing. Subsequently, the remaining 252 samples were randomly assigned to training set (10 GC and 10 healthy controls) and validation set, including 44 healthy controls (age ranged 29-66years), chronic gastritis (age ranged 32-88 years), typical hyperplasia (age ranged 33-76 years), and GC (age ranged 25-82 years). Participants’ demographics and clinicopathological parameters were also recorded. The inclusion criteria were the same as our previous study (36). This study has been approved by the Medical Ethics Committee of Qilu Hospital of Shandong University (KYLL-2015-097).
Extraction of extracellular vesicles
The extraction of extracellular vesicles were performed as described in our previous study (12, 36).
High-throughput sequencing
Total RNAs were extracted from serum EVs, and then Ribo-Zero rRNA Removal (Illumina, USA) was used to deplete ribosomal RNA. RNA libraries were constructed using TruSeq Stranded Total RNA Library Prep Kit (Illumina). BioAnalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA) was used to quantify and quality-control libraries. Then, 10 pM libraries were denatured to single-stranded DNA molecules, captured on Illumina flow cells, amplified in situ as clusters, and sequenced for 150 cycles using an Illumina HiSeq platform.
CircRNAs sequencing data analysis
Briefly, paired-end reads were harvested from Illumina HiSeq 4000 platform and controlled quality by Q30. After 3’ adaptor-trimming and low-quality read were removed by cutadapt software (v1.9.3), the high-quality reads were aligned to the reference genome/transcriptome with STAR software. CircRNAs were detected and identified with DCC software. EdgeR software was used to normalize the data and perform DE circRNA analysis.
Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis
To predict the function of circRNAs, 20 miRNAs with the strongest binding to DE circRNAs, and its target gene 50 mRNAs were selected for GO and KEGG pathway analysis. P-value< 0.05 was considered as the threshold of significant enrichment. The differential genes were grouped into hierarchical categories of molecular function, biological process, and cellular component based on the Gene Ontology database (http://www.geneontology.org/). Pathway analysis was performed to assess the potential pathways of the differentially expressed genes according to KEGG (http://www.genome.jp/kegg/).
Prediction of circRNA-miRNA-mRNA interactions
We predicted the miRNA-binding sites of circRNA using two different algorithms: TargetScan and miRANDA based on database from miRBase release 21.0. MiRNA-mRNA regulatory relationships was identified using miRTarBase database. Cytoscape was applied to draw the graph of a circRNA-miRNA-mRNA interaction network.
RNA extraction and quantitative real-time RT-PCR
Serum EVs RNA was extracted using exoRNeasy Serum/Plasma Midi Kit (Cat:77044, Qiagen, Germany) according to the manufacturer’s instructions. ND-1000 Nanodrop (Thermo Fisher, USA) was applied to test the purity of total RNA. All-in-one first-strand cDNA synthesis kit (Cat: QP006, GeneCopoeia company, Rockville, Maryland) was used to reversely transcribe total RNA to cDNA according to instructions. cDNA was diluted at 1:5. PCR reactions were performed using CFX96 (BIO-RAD, USA) with All-in-one™ qPCR Mix (Cat: QP001, GeneCopoeia company). Amplification cycle conditions were as follows: an initial incubation at 95°C 10min, then 95°C 15s, 62°C 20s, 72°C 10s, and the cycle was 40. The sequences of primers are listed in Table 1. Data was normalized using 18S, and relative expression levels were evaluated using the 2−ΔΔCT method.
Fecal occult blood tests
Bleeding of digestive tracts was measured by the Colloidal gold-based fecal occult blood diagnostic Kit (Chemtron Biotech Co. Shanghai, China).
Statistical analysis
Kolmogorov Smirnov test was used to determine the distribution of data in each group. Data was presented as median and interquartile ranges. Serum EV circRNAs levels were evaluated using Mann Whitney U or Kruskal Wallis test in different groups, including clinicopathological variables. ROC curves were used to evaluate the diagnostic value of the potential biomarkers in GC. The cut-off value was calculated according to Youden index (sensitivity + specificity-1). P<0.05 was considered statistically significant. SPSS 25.0 (Chicago, Illinois, USA) and Medcalc 8.0 (Korea) were used in this study.
Results
Characteristics of the total circRNAs expression
The A260/280 ratio range of EVs RNA of 4 healthy controls and 4 GC were 1.8-2.1. This indicated that the purity of the extracted RNA was high and could be used for high-throughput sequencing (Table 2). Results showed that a total of 4692 circRNAs (reads ≥ 2 in a single sample) were identified, of which 4575 (98%) were novel, and 117 (2%) have been reported in circBase database (Figure 1A). According to the location of their parent genes, circRNAs were divided into five categories. Intergenic, intronic, antisense, sense overlapping, and exonic accounted for 52%, 17%, 17%, 11%, 3%, respectively (Figure 1B). In addition, circRNAs were observed to be widely distributed to all chromosomes according to the location of its transcriptional gene (Figure 1C). Further analysis demonstrated that the length of circRNAs was predominantly < 5,000 bp (Figure 1D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/962831/fonc-12-962831-HTML/image_m/fonc-12-962831-g001.jpg)
Figure 1 The characteristics of circRNAs. (A) The classification of circRNAs identified by sequencing. (B) The categories of circRNAs. (C) The location distribution of circRNAs on chromosomes. (D) The length distribution of circRNAs.
Differential circRNA profiles in serum extracellular vesicles in patients with GC and healthy controls
11 differential circRNAs were initially selected according to the criteria |log2Fold Change| > 2, P < 0.05). 7 were significantly upregulated and 4 were downregulated (Table 3). A cluster heat map (Figure 2A), scatter map (Figure 2B), and volcano map (Figure 2C) were used to visualize the differences of circRNAs in GC and healthy controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/962831/fonc-12-962831-HTML/image_m/fonc-12-962831-g002.jpg)
Figure 2 Differential circRNA expression profiles in serum extracellular vesicles. (A) Cluster heat map. (B) Scatter map. (C) Volcano map.
GO enrichment analysis of DE circular RNA
The enrichment results of molecular function (MF) showed that the upregulated circRNAs were mainly related to glutathione transferase activity and transferase activity (Figure 3A). The major downregulated circRNAs were associated with heat shock protein (HSP) binding and Hsp90 protein binding (Figure 3B). In the biological process (BP), the upregulated circRNAs were involved in glutathione derivative metabolic process and glutathione derivative biosynthetic process (Figure 3C). The downregulated circRNAs were primarily enriched in chaperone-mediated protein folding and protein folding (Figure 3D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/962831/fonc-12-962831-HTML/image_m/fonc-12-962831-g003.jpg)
Figure 3 GO and KEGG analysis of the differentially expressed circRNAs. GO analysis of the upregulated (A) and the downregulated (B) DE circRNAs in terms of molecular function. GO analysis of the upregulated (C) and the downregulated (D) DE circRNAs in biological processes. KEGG enrichment analysis of the differentially expressed circRNAs (E, F). X-axis represented GeneRatio, and it equals the ratio of the number of DE genes associated with the gene ontology term to the total number of DE genes.
KEGG analysis of DE circRNAs
The results of KEGG were shown in a column chart (Figure 3E) and bubble chart (Figure 3F). DE circRNAs were enriched in glutathione metabolism, drug metabolism-cytochrome P450, metabolism of xenobiotics by cytochrome P450, and chemical carcinogenesis.
Preliminary screening of DE circRNAs in training set
First of all, 11 DE circRNAs were verified in training set including 10 healthy controls and 10 GC using RT-qPCR. Results showed that the expression levels of circRNA Chr10q11, Chr1p11, Chr7q11 in GC were significantly higher than those in healthy controls (Figures 4A–C) (all P < 0.05). Whereas circRNA chr6:52650927-52691423+ was lower in GC (P = 0.022) (Figure 4D). Interesting, circRNA chr21:9827153-9827349+ and chr17:33478113-33478328- were not significantly different between two groups (Figures 4E, F). The remaining 5 circRNAs were not amplified or nonspecific amplification. Remarkably, the expression of circRNA chr6:52650927-52691423+ in training set was inconsistent with the results of RNA-seq, and it was not considered in the following experiments.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/962831/fonc-12-962831-HTML/image_m/fonc-12-962831-g004.jpg)
Figure 4 Verification of differentially expressed circRNAs in training set. The expression of 6 EVs circRNAs Chr10q11 (A), Chr1p11 (B), Chr7q11 (C) chr6:52650927-52691423+ (D), chr21:9827153-9827349+ (E) and chr17:33478113-33478328- (F) between healthy controls and gastric cancer.
CircRNAs feature verification
The 3 DE circRNAs, Chr10q11, Chr1p11, and Chr7q11, were amplified by RT-qPCR, and PCR products were identified using Sanger sequencing. Results showed that 3 circRNA products had specific sequences of cyclization sites, and the bases were TT, AC, and TG, respectively (Supplementary Figures 1A–C). To further observe the stability of 3 circRNAs, RT-qPCR was used to detect total RNA samples treated with or without RNase R, respectively. Results showed that after RNase R digestion, the expression levels of circRNAs did not differ significantly, but linear RNA GAPDH were significantly downregulated (Supplementary Figures 2, P=0.01).
DE circRNAs in validation set
In validation set, the expression of circRNA Chr10q11 in GC [3.395 (1.455-4.743)] was higher than that in healthy controls [1.000 (0.378-2.395)], chronic gastritis [1.060 (0.495-2.025)] and atypical hyperplasia [1.950 (1.180-3.900)], respectively. Atypical hyperplasia was higher than that of healthy controls and chronic gastritis (Figure 5A). Similarly, the expression of circRNA Chr1p11 in GC [3.600 (2.325-5.625)] was higher than that in healthy controls [0.995 (0.495-1.965)], chronic gastritis [1.100 (0.325-4.105)] and atypical hyperplasia [2.430 (1.200-3.660)], respectively. Chr1p11 was higher in atypical hyperplasia than in healthy controls and chronic gastritis (Figure 5B). Figure 5C showed that circRNA Chr7q11 in GC [2.403 (1.225-4.639)] was higher than that in healthy control [1.001 (0.279-2.061)], chronic gastritis [1.183 (0.419-3.130)] and atypical hyperplasia [1.721 (0.520-5.247)]. Additionally, we investigated the expression of circRNAs in the different TNM stages. The levels of serum EVs circRNA Chr10q11 (Figure 5D), Chr1p11 (Figure 5E) and Chr7q11 (Figure 5F) in stage III + IV were all significantly higher than that in stage I + II (all P < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/962831/fonc-12-962831-HTML/image_m/fonc-12-962831-g005.jpg)
Figure 5 Expression of differentially expressed circRNAs in validation set. The expressions of circRNA Chr10q11 (A), Chr1p11 (B), and Chr7q11 (C) in different groups including healthy control, chronic gastritis, atypical hyperplasia and gastric cancer. Meanwhile, the expression of circRNA Chr10q11 (D), Chr1p11 (E), Chr7q11 (F) in stage I + II and stage III + IV were also compared.
Relationship between DE circRNAs and clinicopathological parameters
Results from Table 4 revealed that tumor size, lymph node metastasis, distant metastasis, and TNM stage were all correlated with 3 circRNAs including Chr10q11, Chr1p11, Chr7q11 (all P < 0.05). Sex, age, blood pressure, location, or fecal occult blood were all not related to the above 3 circRANs (all P > 0.05). In addition, circRNA Chr10q11 and Chr1p11 were both associated with invasion depth (all P < 0.05). But circRNA Chr10q11 was observed to be not significantly correlated with Bormann classification (P = 0.035).
Diagnostic value of DE circRNAs
The AUC of circRNA Chr10q11, Chr1p11, Chr7q11 and CEA were 0.726, 0.822, 0.749 and 0.691, respectively. The sensitivity and specificity of these above molecular were 76.8% and 65.9%, 82.1% and 77.3%, 79.5% and 59.1%, 70.5% and 63.6%, respectively (Figures 6A–D, Table 5). The AUC of the combined 3 circRNAs was 0.839 (95%CI: 0.772-0.893) with sensitivity 73.2% and specificity 84.1% (Figure 6E). When combining 3 circRANs and CEA, the AUC was 0.866 (95%CI: 0.803-0.915). The sensitivity and specificity were 80.4% and 81.8% (Figure 6F; Table 5), achieving a high classification power for GC and healthy controls.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/962831/fonc-12-962831-HTML/image_m/fonc-12-962831-g006.jpg)
Figure 6 Diagnostic value of differentially expressed circRNAs. ROC curves were used to evaluate the diagnostic value of circRNA Chr10q11 (A), Chr1p11 (B), Chr7q11 (C), CEA (D), 3 circRNAs combination (E), 3 circRNAs + CEA (F).
Predicted circRNA-miRNA-mRNA network based on the DE circRNAs
The 3 selected DE circRNAs, Chr10q11, Chr1p11, and Chr7q11, were used to predict the ceRNA network of circRNA-miRNA-mRNA. First, we analyzed mRNA expression levels in sequencing data (Figure 7A). A total of 20,308 mRNAs (13,991 upregulated mRNAs and 6,317 downregulated mRNAs) were obtained. 704 DE mRNAs were screened according to FC > 1.5, P < 0.05. Next, TargetScan and miRANDA were used to predict the binding miRNAs, and the top five miRNAs with the strongest binding to each circRNA were selected. Based on this, the top 200 mRNAs with the strongest binding to each miRNA were also predicted. Subsequently, these target genes were intersected with the DE mRNAs in sequencing data. Finally, 37 upregulated mRNAs and 54 downregulated mRNAs were obtained (Figure 7B), and then the circRNA-miRNA-mRNA network was drawn. Figure 7C showed that DE circRNAs were interacted with 13 miRNAs (hsa-miR-3198, hsa-miR-4677-3p, hsa-miR-450b-3p, hsa-miR-944, hsa-miR-513b-5p, hsa-miR-6888-3p, hsa-miR-580-3p, hsa-miR-219a-1-3p, hsa-miR-4251, hsa-miR-4457, hsa-miR-136-3p, hsa-miR-3174, and hsa-miR-616-3p) and 91 mRNAs. By visualizing the ceRNA network, we can directly observe the center of ceRNA network and the relationship between these components. This provides a reference for future research on the identified circRNAs. In addition, the function of mRNAs were analyzed by GO and KEGG. Results showed that GO (Supplementary Figures 3A–F) and KEGG (Supplementary Figures 4A–D) were mainly enriched in protein binding, biological regulation, endocytosis, and PI3K-Akt signaling pathway.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/962831/fonc-12-962831-HTML/image_m/fonc-12-962831-g007.jpg)
Figure 7 Construction of circRNA-miRNA-mRNA network. (A) showed the heatmap of EV mRNAs in sequencing results. (B) The intersection of mRNAs from sequencing results and mRNAs from prediction by DE circRNAs. (C) Construction of circRNA-miRNA-mRNA network. Circles represent circRNAs, arrows represent miRNAs, and hexagons represent mRNAs.
Discussion
With the improvement of living standards and irregular diets, patients with GC is increasing annually. GC has become one of the most common high-incidence malignancies, and over 40% of new cases occur in China (37, 38). Recently, EVs have attracted increasing attention as a novel pathway of intracellular communication. Tumor-derived EVs regulate gastric carcinogenesis and drug resistance by delivering some molecules such as lncRNAs, miRNAs, circRNAs to recipient cells (39). CircRNAs has been proved to be able to participate in the occurrence and development of cancer through competing for miRNAs or via other mechanisms (40). Most patients with GC are diagnosed in the first instance with advanced stage, thus the survival rate is low. In this study, 4692 circRNAs were identified, of which 11 were differentially expressed. GO and KEGG analyses of DE circRNAs showed that they were involved in several molecular biological pathways such as glutathione metabolism. Of these, the expression of 3 GC-associated DE circRNA, Chr10q11, Chr1p11, and Chr7q11 were upregulated in GC, and has a good diagnostic value. The construction of a circRNA-miRNA-mRNA network showed 13 miRNA and 91 mRNA interactions. Taken together, our findings provide a new direction for improving the early diagnosis of GC.
EVs are a signaling vector for the exchange of ‘information’ from one cell to another. The cell-derived EVs release into circulation and play a complex role in tumor development, including angiogenesis and inflammation, immunosuppression, and so on (41, 42). Zhang et al. (43) demonstrated that EVs from GC cells induced neutrophil polarization to N2, promoting GC cell migration. EVs are able to function as an intercellular communication medium by various contents (proteins, lipids, DNAs, mRNAs, ncRNAs, etc.). Ye et al. found that hsa_circ_0000069 was significantly upregulated and inhibited the proliferation and migration of human pancreatic ductal epithelial cells by inhibiting STIL expression (44). Studies have shown that EVs are the important carriers of circRNAs. CircRNAs are more stable in serum due to the protection provided by EVs (45). Our results showed that circRNAs were more stable in serum. This provides a basis for EVs circRNAs to be used as biomarkers of GC.
CircRNAs have been found to play roles in several biological processes such as miRNA binding, protein binding, and transcriptional and post-transcriptional regulation. To explore the potential function of circRNAs in this study, we performed GO and KEGG pathway analysis. Results of GO analysis indicated that the most significantly enriched terms were heat shock protein binding and glutathione derivative metabolic process, indicating an association with the basic pathophysiological processes of cancer. It has been shown that the high expression of HSPs is significantly associated with the reduced OS in cancer patients (46). KEGG pathway analysis also revealed that the DE circRNAs were mainly involved in pathways such as glutathione metabolism, drug metabolism-cytochrome P450, and chemical carcinogenesis. The cytochrome P450, a monooxygenase, is the oxidative metabolism of exogenous substances. It has also been found that cytochrome P450 is overexpressed in cancer tissues and can predict the response to chemotherapy (47). This was also consistent with the findings of Luo et al. (48) that higher cytochrome P450 2U1 (CYP2U1) levels were predictive of a poor 5-year OS and were associated with the histopathological grade of BCa. These results suggest that dysregulation of circRNAs in GC may be involved in tumorigenesis.
In clinical practice, the diagnosis of certain tumor-related diseases often requires histopathological confirmation, leading to unnecessary injury. There is an urgent need for simpler and noninvasive diagnostic methods for GC. With the development of liquid biopsy technology, EV circRNAs have become a research hotspot in the field of tumor markers, and has a potential in the diagnosis of various diseases (49, 50). In the current study, serum EV circRNA Chr10q11, Chr1p11, and Chr7q11 expression levels were significantly higher in GC, and related to tumor size, lymph node metastasis, distant metastasis, and TNM stage. Several studies also suggest that EV circRNAs in fluids is associated with pathological features of tumor vascular invasion and TNM stage. Wang et al. found that circ-ITCH was downregulated in serum exosomes, and acted as a sponge of miR-199a-5p, inhibiting proliferation and epithelial-mesenchymal transformation of GC (51). The development of GC is a complex process. The roles of EVs circRNAs in GC need to be further explored.
Patients with GC lack typical symptoms in the early stages, and the sensitivity and specificity of the traditional biomarkers is insufficient. Serum EV circRNAs are easy to be obtained and are widely present in various body fluids, providing a novel approaches for GC diagnosis. Our results showed the AUC of 3 DE circRNAs combined with CEA was 0.866. The sensitivity and specificity were 80.4%, and 81.8%, indicating a better value for GC diagnosis. Previous studies have reported that hsa_circ_0065149, circ-KIAA1244, and hsa_circ_0000419 can be used as tumor markers for GC screening. The sensitivity and specificity are 74.8% and 84% (52–54). These results show that EV circRNAs may help to distinguish patients with GC from healthy controls.
Notably, circRNAs contain miRNA binding sites. It can serve as miRNA sponges to alleviate the inhibitory effects. However, the underlying mechanisms in GC remains unclear. To further explore the potential role of circRNAs in GC, we constructed circRNA-miRNA-mRNA network. Results revealed that the DE circRNAs interacted with 13 miRNAs and 91 mRNAs, which were ultimately involved in the regulation of gene expression. MiRNAs have been confirmed to function by binding to mRNAs. For example, in colon cancer, miR-7 binding to the 3 ′ UTR of TFF3 via the PI3K-AKT signaling pathway to inhibit proliferation and migration (55). The abnormally expressed circRNAs may participate in the pathogenesis of GC by regulating miRNAs and their target mRNAs. Hu et al. found that CircPIP5K1A activates KRT80 and PI3K/AKT pathway to promote GC development through sponging miR-671-5p (56). In addition, a study showed that circRAB31, a sponge of mir-885-5p, inhibits GC procession via PTEN/PI3K/AKT pathway (57). These results suggest that circRNAs can play a significant role in the occurrence and development of GC through circRNA-miRNA-mRNA networks. However, the present study has still some limitations, including the relatively small sample size and lacking of mechanism research. Further study is guaranteed to determine its clinical predictive value and to elucidate the underlying mechanism in gastric cancer.
Conclusions
Our findings showed that EVs circRNAs in serum were significantly upregulated in gastric cancer, and related to TNM stage and lymphatic metastasis. It has a good diagnostic value for gastric cancer, which will provide novel insights into early diagnosis for gastric cancer.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/, GSE165394.
Ethics statement
The studies involving human participants were reviewed and approved by Medical Ethics Committee of Qilu Hospital of Shandong University (KYLL-2015-097). The patients/participants provided their written informed consent to participate in this study.
Author contributions
KX, ZD and YZ designed the experiments; KX and JD performed the experiments; DW collected the samples; ZW and XC analyzed the data; KX and ZD created the figure and wrote the manuscript; SL and YZ provided comments and suggestions. All the authors reviewed and approved the final version of the manuscript.
Funding
This research was funded by the National Natural Science Foundation of China (82072358, 81972005, 81572070). Latitudinal research project of Shandong University (6010119015). The Major Scientific and Technological Innovation Project of Shandong Province (2021CXGC010603).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.962831/full#supplementary-material
Supplementary Figure 1 | Sanger sequencing of circRNAs. Sanger sequencing of CircRNAs Chr10q11 (A), Chr1p11 (B), and Chr7q11 (C). The direction of arrows represented splice junction of circRNAs.
Supplementary Figure 2 | The effect of RNase R on EV circRNAs. The expression of circRNAs Chr10q11, Chr1p11, Chr7q11 and GAPDH with or without RNase R treatment.
Supplementary Figure 3 | GO analysis for mRNA from intersection. The CC, MF, and BP of the upregulated mRNAs (A–C) and the downregulated mRNAs (D–F).
Supplementary Figure 4 | KEGG analysis for mRNA from intersection. KEGG enrichment analysis of the upregulated mRNAs (A, B) and the downregulated mRNAs (C, D).
References
1. Smyth EC, Nilsson M, Grabsch HI, van Grieken NC, Lordick F. Gastric cancer. Lancet (2020) 396(10251):635–48. doi: 10.1016/S0140-6736(20)31288-5
2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: Globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2018) 68(6):394–424. doi: 10.3322/caac.21492
3. Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin (2016) 66(2):115–32. doi: 10.3322/caac.21338
4. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin (2015) 65(2):87–108. doi: 10.3322/caac.21262
5. Tan Z. Recent advances in the surgical treatment of advanced gastric cancer: A review. Med Sci Monit (2019) 25:3537–41. doi: 10.12659/MSM.916475
6. Kim GH, Liang PS, Bang SJ, Hwang JH. Screening and surveillance for gastric cancer in the united states: Is it needed? Gastrointest Endosc (2016) 84(1):18–28. doi: 10.1016/j.gie.2016.02.028
7. Fu M, Gu J, Jiang P, Qian H, Xu W, Zhang X. Exosomes in gastric cancer: Roles, mechanisms, and applications. Mol Cancer (2019) 18(1):41. doi: 10.1186/s12943-019-1001-7
8. Cummings D, Wong J, Palm R, Hoffe S, Almhanna K, Vignesh S. Epidemiology, diagnosis, staging and multimodal therapy of esophageal and gastric tumors. Cancers (Basel) (2021) 13(3):582. doi: 10.3390/cancers13030582
9. Liu T, Jiao Y. [Pathological diagnosis of gastric intraepithelial neoplasia and early gastric cancer ingastric biopsy specimens]. Zhonghua Bing Li Xue Za Zhi (2014) 43(9):644–6. doi: 10.3760/cma.j.issn.0529-5807.2014.09.020
10. Ito M, Tanaka S, Chayama K. Characteristics and early diagnosis of gastric cancer discovered after helicobacter pylori eradication. Gut Liver (2021) 15(3):338–45. doi: 10.5009/gnl19418
11. Thanh Huong P, Gurshaney S, Thanh Binh N, Gia Pham A, Hoang Nguyen H, Thanh Nguyen X, et al. Emerging role of circulating tumor cells in gastric cancer. Cancers (Basel) (2020) 12(3):695. doi: 10.3390/cancers12030695
12. Dong Z, Sun X, Xu J, Han X, Xing Z, Wang D, et al. Serum membrane type 1-matrix metalloproteinase (Mt1-mmp) mrna protected by exosomes as a potential biomarker for gastric cancer. Med Sci Monit (2019) 25:7770–83. doi: 10.12659/MSM.918486
13. Goodall GJ, Wickramasinghe VO. Rna in cancer. Nat Rev Cancer (2021) 21(1):22–36. doi: 10.1038/s41568-020-00306-0
14. Xie S, Chang Y, Jin H, Yang F, Xu Y, Yan X, et al. Non-coding rnas in gastric cancer. Cancer Lett (2020) 493:55–70. doi: 10.1016/j.canlet.2020.06.022
15. Zhou X, Zhan L, Huang K, Wang X. The functions and clinical significance of circrnas in hematological malignancies. J Hematol Oncol (2020) 13(1):138. doi: 10.1186/s13045-020-00976-1
16. Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, et al. Circular rna: A new star of noncoding rnas. Cancer Lett (2015) 365(2):141–8. doi: 10.1016/j.canlet.2015.06.003
17. Ren S, Lin P, Wang J, Yu H, Lv T, Sun L, et al. Circular rnas: Promising molecular biomarkers of human aging-related diseases Via functioning as an mirna sponge. Mol Ther Methods Clin Dev (2020) 18:215–29. doi: 10.1016/j.omtm.2020.05.027
18. Wong CH, Lou UK, Fung FK, Tong JHM, Zhang CH, To KF, et al. Circrtn4 promotes pancreatic cancer progression through a novel circrna-Mirna-Lncrna pathway and stabilizing epithelial-mesenchymal transition protein. Mol Cancer (2022) 21(1):10. doi: 10.1186/s12943-021-01481-w
19. Hang D, Zhou J, Qin N, Zhou W, Ma H, Jin G, et al. A novel plasma circular rna circfarsa is a potential biomarker for non-small cell lung cancer. Cancer Med (2018) 7(6):2783–91. doi: 10.1002/cam4.1514
20. Ye DX, Wang SS, Huang Y, Chi P. A 3-circular rna signature as a noninvasive biomarker for diagnosis of colorectal cancer. Cancer Cell Int (2019) 19:276. doi: 10.1186/s12935-019-0995-7
21. Jin Z, Gao B, Gong Y, Guan L. Depletion of circ-Birc6, a circular rna, suppresses non-small cell lung cancer progression by targeting mir-4491. Biosci Trends (2021) 14(6):399–407. doi: 10.5582/bst.2020.03310
22. Zong L, Sun Q, Zhang H, Chen Z, Deng Y, Li D, et al. Increased expression of Circrna_102231 in lung cancer and its clinical significance. BioMed Pharmacother (2018) 102:639–44. doi: 10.1016/j.biopha.2018.03.084
23. Chen Q, Liu T, Bao Y, Zhao T, Wang J, Wang H, et al. Circrna Crapgef5 inhibits the growth and metastasis of renal cell carcinoma Via the mir-27a-3p/Txnip pathway. Cancer Lett (2020) 469:68–77. doi: 10.1016/j.canlet.2019.10.017
24. Chen J, Li Y, Zheng Q, Bao C, He J, Chen B, et al. Circular rna profile identifies Circpvt1 as a proliferative factor and prognostic marker in gastric cancer. Cancer Lett (2017) 388:208–19. doi: 10.1016/j.canlet.2016.12.006
25. Zang X, Jiang J, Gu J, Chen Y, Wang M, Zhang Y, et al. Circular rna Eif4g3 suppresses gastric cancer progression through inhibition of beta-catenin by promoting delta-catenin ubiquitin degradation and upregulating Sik1. Mol Cancer (2022) 21(1):141. doi: 10.1186/s12943-022-01606-9
26. Guo R, Cui X, Li X, Zang W, Chang M, Sun Z, et al. Circman1a2 is upregulated by helicobacter pylori and promotes development of gastric cancer. Cell Death Dis (2022) 13(4):409. doi: 10.1038/s41419-022-04811-y
27. Saaoud F, Drummer IVC, Shao Y, Sun Y, Lu Y, Xu K, et al. Circular rnas are a novel type of non-coding rnas in ros regulation, cardiovascular metabolic inflammations and cancers. Pharmacol Ther (2021) 220:107715. doi: 10.1016/j.pharmthera.2020.107715
28. Long F, Lin Z, Li L, Ma M, Lu Z, Jing L, et al. Comprehensive landscape and future perspectives of circular rnas in colorectal cancer. Mol Cancer (2021) 20(1):26. doi: 10.1186/s12943-021-01318-6
29. Yang SJ, Wang DD, Zhou SY, Zhang Q, Wang JY, Zhong SL, et al. Identification of circrna-mirna networks for exploring an underlying prognosis strategy for breast cancer. Epigenomics (2020) 12(2):101–25. doi: 10.2217/epi-2019-0058
30. Chen X, Chen RX, Wei WS, Li YH, Feng ZH, Tan L, et al. Prmt5 circular rna promotes metastasis of urothelial carcinoma of the bladder through sponging mir-30c to induce epithelial-mesenchymal transition. Clin Cancer Res (2018) 24(24):6319–30. doi: 10.1158/1078-0432.CCR-18-1270
31. Lu J, Wang YH, Yoon C, Huang XY, Xu Y, Xie JW, et al. Circular rna circ-Rangap1 regulates vegfa expression by targeting mir-877-3p to facilitate gastric cancer invasion and metastasis. Cancer Lett (2020) 471:38–48. doi: 10.1016/j.canlet.2019.11.038
32. Lister R, O'Malley RC, Tonti-Filippini J, Gregory BD, Berry CC, Millar AH, et al. Highly integrated single-base resolution maps of the epigenome in arabidopsis. Cell (2008) 133(3):523–36. doi: 10.1016/j.cell.2008.03.029
33. Wang Y, Mashock M, Tong Z, Mu X, Chen H, Zhou X, et al. Changing technologies of rna sequencing and their applications in clinical oncology. Front Oncol (2020) 10:447. doi: 10.3389/fonc.2020.00447
34. Zhu M, Dang Y, Yang Z, Liu Y, Zhang L, Xu Y, et al. Comprehensive rna sequencing in adenoma-cancer transition identified predictive biomarkers and therapeutic targets of human crc. Mol Ther Nucleic Acids (2020) 20:25–33. doi: 10.1016/j.omtn.2020.01.031
35. Long MY, Chen JW, Zhu Y, Luo DY, Lin SJ, Peng XZ, et al. Comprehensive circular rna profiling reveals the regulatory role of Circrna_0007694 in papillary thyroid carcinoma. Am J Transl Res (2020) 12(4):1362–78.
36. Xiao K, Dong Z, Wang D, Liu M, Ding J, Chen W, et al. Clinical value of lncrna Ccat1 in serum extracellular vesicles as a potential biomarker for gastric cancer. Oncol Lett (2021) 21(6):447. doi: 10.3892/ol.2021.12708
37. Zhang J, Liu H, Hou L, Wang G, Zhang R, Huang Y, et al. Circular Rna_Larp4 inhibits cell proliferation and invasion of gastric cancer by sponging mir-424-5p and regulating Lats1 expression. Mol Cancer (2017) 16(1):151. doi: 10.1186/s12943-017-0719-3
38. Jin G, Lv J, Yang M, Wang M, Zhu M, Wang T, et al. Genetic risk, incident gastric cancer, and healthy lifestyle: A meta-analysis of genome-wide association studies and prospective cohort study. Lancet Oncol (2020) 21(10):1378–86. doi: 10.1016/S1470-2045(20)30460-5
39. Liu X, Lu Y, Xu Y, Hou S, Huang J, Wang B, et al. Exosomal transfer of mir-501 confers doxorubicin resistance and tumorigenesis Via targeting of blid in gastric cancer. Cancer Lett (2019) 459:122–34. doi: 10.1016/j.canlet.2019.05.035
40. Huang XX, Zhang Q, Hu H, Jin Y, Zeng AL, Xia YB, et al. A novel circular rna Circfn1 enhances cisplatin resistance in gastric cancer via sponging mir-182-5p. J Cell Biochem (2020) 122(9):1009–20. doi: 10.1002/jcb.29641
41. van Niel G, D'Angelo G, Raposo G. Shedding light on the cell biology of extracellular vesicles. Nat Rev Mol Cell Biol (2018) 19(4):213–28. doi: 10.1038/nrm.2017.125
42. Hoshino A, Costa-Silva B, Shen TL, Rodrigues G, Hashimoto A, Tesic Mark M, et al. Tumour exosome integrins determine organotropic metastasis. Nature (2015) 527(7578):329–35. doi: 10.1038/nature15756
43. Zhang X, Shi H, Yuan X, Jiang P, Qian H, Xu W. Tumor-derived exosomes induce N2 polarization of neutrophils to promote gastric cancer cell migration. Mol Cancer (2018) 17(1):146. doi: 10.1186/s12943-018-0898-6
44. Ye Z, Zhu Z, Xie J, Feng Z, Li Y, Xu X, et al. Hsa_Circ_0000069 knockdown inhibits tumorigenesis and exosomes with downregulated Hsa_Circ_0000069 suppress malignant transformation Via inhibition of stil in pancreatic cancer. Int J Nanomedicine (2020) 15:9859–73. doi: 10.2147/IJN.S279258
45. Li Y, Zheng Q, Bao C, Li S, Guo W, Zhao J, et al. Circular rna is enriched and stable in exosomes: A promising biomarker for cancer diagnosis. Cell Res (2015) 25(8):981–4. doi: 10.1038/cr.2015.82
46. Rodina A, Wang T, Yan P, Gomes ED, Dunphy MP, Pillarsetty N, et al. The epichaperome is an integrated chaperome network that facilitates tumour survival. Nature (2016) 538(7625):397–401. doi: 10.1038/nature19807
47. Elfaki I, Mir R, Almutairi FM, Duhier FMA. Cytochrome P450: Polymorphisms and roles in cancer, diabetes and atherosclerosis. Asian Pac J Cancer Prev (2018) 19(8):2057–70. doi: 10.22034/APJCP.2018.19.8.2057
48. Luo B, Chen C, Wu X, Yan D, Chen F, Yu X, et al. Cytochrome P450 2u1 is a novel independent prognostic biomarker in breast cancer patients. Front Oncol (2020) 10:1379. doi: 10.3389/fonc.2020.01379
49. Wang S, Dong Y, Gong A, Kong H, Gao J, Hao X, et al. Exosomal circrnas as novel cancer biomarkers: Challenges and opportunities. Int J Biol Sci (2021) 17(2):562–73. doi: 10.7150/ijbs.48782
50. Cheng J, Meng J, Zhu L, Peng Y. Exosomal noncoding rnas in glioma: Biological functions and potential clinical applications. Mol Cancer (2020) 19(1):66. doi: 10.1186/s12943-020-01189-3
51. Wang Y, Wang H, Zheng R, Wu P, Sun Z, Chen J, et al. Circular rna itch suppresses metastasis of gastric cancer via regulating mir-199a-5p/Klotho axis. Cell Cycle (2021) 20(5-6):522–36. doi: 10.1080/15384101.2021.1878327
52. Tang W, Fu K, Sun H, Rong D, Wang H, Cao H. Circrna microarray profiling identifies a novel circulating biomarker for detection of gastric cancer. Mol Cancer (2018) 17(1):137. doi: 10.1186/s12943-018-0888-8
53. Tao X, Shao Y, Lu R, Ye Q, Xiao B, Ye G, et al. Clinical significance of Hsa_Circ_0000419 in gastric cancer screening and prognosis estimation. Pathol Res Pract (2020) 216(1):152763. doi: 10.1016/j.prp.2019.152763
54. Shao Y, Tao X, Lu R, Zhang H, Ge J, Xiao B, et al. Hsa_Circ_0065149 is an indicator for early gastric cancer screening and prognosis prediction. Pathol Oncol Res (2020) 26(3):1475–82. doi: 10.1007/s12253-019-00716-y
55. Guo J, Xu L, Teng X, Sun M. Microrna-7-5p regulates the proliferation and migration of intestinal epithelial cells by targeting trefoil factor 3 Via inhibiting the phosphoinositide 3-Kinase/Akt signalling pathway. Int J Mol Med (2017) 40(5):1435–43. doi: 10.3892/ijmm.2017.3120
56. Song H, Xu Y, Xu T, Fan R, Jiang T, Cao M, et al. Circpip5k1a activates Krt80 and Pi3k/Akt pathway to promote gastric cancer development through sponging mir-671-5p. BioMed Pharmacother (2020) 126:109941. doi: 10.1016/j.biopha.2020.109941
Keywords: gastric cancer, circRNA, extracellular vesicle, diagnosis, biomarker
Citation: Xiao K, Li S, Ding J, Wang Z, Wang D, Cao X, Zhang Y and Dong Z (2022) Expression and clinical value of circRNAs in serum extracellular vesicles for gastric cancer. Front. Oncol. 12:962831. doi: 10.3389/fonc.2022.962831
Received: 06 June 2022; Accepted: 28 July 2022;
Published: 17 August 2022.
Edited by:
Xianda Zhao, University of Minnesota, United StatesReviewed by:
Lianjun Ma, Jilin University, ChinaGuojun Zhang, Beijing Tiantan Hospital, Capital Medical University, China
Copyright © 2022 Xiao, Li, Ding, Wang, Wang, Cao, Zhang and Dong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhaogang Dong, zhaogang.dong@email.sdu.edu.cn; Yi Zhang, yizhang@sdu.edu.cn