Circular RNA hsa_circ_0000745 may serve as a diagnostic marker for gastric cancer

Abstract

AIM

To determine whether circular RNAs (circRNAs) are involved in pathological processes of gastric cancer (GC).

METHODS

Three circRNAs with differential expression in GC and colorectal cancer were randomly selected for validation by quantitative reverse transcription-polymerase chain reaction (qRT-PCR), using 20 pairs of gastric tissues and normal tissues. Based on the predicted circRNA-miRNA network, we then focused on hsa_circ_0000745, which was found to be down-regulated in 20 GC tissues compared with normal tissues. The hsa_circ_0000745 levels were further analyzed by qRT-PCR in 60 GC tissues and paired adjacent non-tumor tissues, as well as 60 plasma samples from GC patients and 60 plasma samples from healthy controls. The associations between the levels of hsa_circ_0000745 and the clinicopathological features of GC patients were statistically assessed. A receiver operating characteristic (ROC) curve was used to evaluate the diagnostic value of hsa_circ_0000745 in GC.

RESULTS

Hsa_circ_0000745 was down-regulated in GC tissues vs non-tumorous tissues (P < 0.001) and in plasma samples from patients with GC vs healthy controls (P < 0.001). The expression level of hsa_circ_0000745 in GC tissues correlated with tumor differentiation, while the expression level in plasma correlated with tumor-node-metastasis stage. The area under the ROC curve (AUC) of hsa_circ_0000745 in plasma was 0.683, suggesting good diagnostic value. Plasma hsa_circ_0000745 level combined with carcinoembryogenic antigen (CEA) level increased the AUC to 0.775.

CONCLUSION

Hsa_circ_0000745 plays an important role in GC and its expression level in plasma in combination with CEA level is a promising diagnostic marker for this malignancy.

Key Words: Circular RNA, circ_0000745, Gastric cancer, Diagnosis, Noncoding RNA

Core tip: This study demonstrates that expression of hsa_circ_0000745 is down-regulated in human gastric cancer (GC), with significantly lower levels detected in human GC tissues compared to paired adjacent non-tumor tissues, and in plasma of GC patients compared to plasma from normal healthy controls. Detection of this stable circular RNA in plasma may represent a clinically useful and convenient marker for GC diagnosis, with higher sensitivity and specificity when analyzed in combination with carcinoembryogenic antigen.

INTRODUCTION

Gastric cancer (GC) is one of the most common malignancies worldwide, and is the most frequently diagnosed cancer among East Asian populations[1-3]. Despite extensive research attention, the mortality rate of GC remains high due to the fact that it is usually diagnosed at late stages. To improve early diagnosis and reduce the high mortality rate of GC[4-6], it is critical to continue elucidating the molecular mechanisms of GC, which may help identify markers for earlier diagnosis and novel therapeutic targets.

Endogenous noncoding RNAs, including microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and the recently identified circular RNAs (circRNAs), play key regulatory roles in various cellular physiological processes and may contribute to pathological processes such as cancer formation and progression[7,8]. Unlike linear RNAs which terminate with 5' caps and 3' tails, circRNAs form covalently closed loop structures with neither 5' to 3' polarity nor a polyadenylated tail. They are not simply by-products of mis-spliced RNAs or splicing errors, but rather they are the products of regulated back-spliced RNAs with distinct sets of cis-elements and/or trans-factors[9]. Moreover, the formation of circRNAs involves the occurrence of back-splicing by the canonical spliceosome[10].

Recent research has demonstrated that circRNAs are widely involved in myriad physiological and pathological processes, including neuronal differentiation, synaptogenesis, neurological disorders, angiogenesis, prion disease, and cancer[11-16]. Likewise, circRNAs have been described as a class of biomarkers related to aging in Drosophila[17] and as putative disease biomarkers in human saliva[18]. Li et al[19] reported that hsa_circ_002059, a typical circRNA, was significantly down-regulated in GC tissues compared with paired adjacent non-tumorous tissues, suggesting its potential as a novel and stable biomarker for the diagnosis of GC.

The collective findings from the research on circRNAs indicate that they may play crucial roles in the molecular pathways of cancer. However, only a few publications have reported on the relationship between circRNAs and cancer. Whether and how circRNAs expression is dysregulated in GC remains unknown.

In this study, we randomly selected three differentially expressed circRNAs (hsa_circ_0000745, hsa_circ_0085616, and chr16:30740286-30740893|+) from two circRNA databases: CircBase (http://circbase.org/) and circ2Traits (http://gyanxet-beta.com/circdb/), which had been characterized using RNA-Seq analysis of ribosomal RNA-depleted total RNA from three GC or colorectal cancerous tissues and paired normal tissues[20,21] (Table 1). We then focused our attention on circRNA hsa_circ_0000745, which is encoded by the sperm antigen with calponin homology and coiled-coil domains 1 gene (SPECC1) that is located at chr17:20107654-20109225. We selected hsa_circ_0000745 for investigation based on preliminary broad sequencing analysis and prediction of its binding to various miRNAs with known relations to GC. By expanding the numbers of human GC and control samples for investigation, we were able to identify a statistically significant down-regulation of hsa_circ_0000745 expression in both GC tissues and plasma samples from patients with GC. Moreover, we were able to determine that the down-regulated expression correlated with tumor differentiation. Overall, the data identifies hsa_circ_0000745 as a novel potential biomarker for diagnosing GC.

Table 1 The three differently expressed circRNAs in human gastric cancer and colorectal cancer samples.

CircRNA	Chrom	TxStart	TxEnd	circBaseID	Best_transcript	Gene name	Catalog
Chr16: 30740287-30740893|+	16	30740286	30740893	266	NM_006662	SRCAP	Exonic
Hsa_circ_0000745	17	20107645	20109225	790	NM_152904	SPECC1	Exonic
Hsa_circ_0085616	8	131370262	131374017	127	NM_001247996	ASAP1	Exonic

Chrom: Chromosome; Tx: Transcription.

MATERIALS AND METHODS

Patient samples

This study was approved by the Institutional Ethics Review Board of Anhui Provincial Hospital Affiliated to Anhui Medical University (China) and was conducted according to the Ethical Guidelines for Human Genome/Gene Research issued by the Chinese Government. Human samples of GC tissues and paired adjacent non-tumor tissues were prospectively collected from 60 patients at the Anhui Provincial Hospital from January 2016 to January 2017. Upon removal, the respective tissues were immediately submerged in RNAlater RNA stabilization reagent (Qiagen, Hilden, Germany) and stored at -80 °C until analysis. None of the patients had undergone chemotherapy or radiotherapy prior to the operation (resection).

The GC diagnosis was confirmed by histopathology. Tumor stages were determined according to the American Joint Committee on Cancer tumor-node-metastasis (TNM) staging criteria.

Before any treatment was applied, each study participant’s peripheral blood (2-3 mL) was collected in an ethylenediaminetetraacetic acid tube, centrifuged at 3000 × g for 10 min to obtain plasma, and stored at -80 °C. Following the age- and sex-matching criteria of the study design, fresh normal plasma samples were collected from 60 healthy volunteers.

All research subjects provided written informed consent prior to study enrollment.

RNA extraction and quality control

Total RNA was isolated from three samples of gastric tumor and paired adjacent non-tumor (normal) tissues, as well as from plasma samples, using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s protocol. RNA integrity was assessed using standard denaturing agarose gel electrophoresis. Total RNA from each specimen was quantified and quality assurance was provided using a NanoDrop ND-1000 spectrophotometer (NanoDrop, Wilmington, DE, United States).

Transcriptome high-throughput sequencing and subsequent bioinformatics analysis

Transcriptome high-throughput sequencing and subsequent bioinformatics analysis were performed by Cloud-Seq Biotech (Shanghai, China). The paired-end reads were harvested on a HiSeq 4000 System (Illumina, San Diego, CA, United States). The high-quality reads were aligned to the reference genome/transcriptome with STAR software[22], and circRNAs were detected and identified with DCC software[23]. EdgeR software[24] was used to normalize the data and perform analysis to determine differential expression among the circRNAs.

Quantitative reverse transcription-polymerase chain reaction validation assay

Total RNA samples were reverse-transcribed into cDNA with a random primer using SuperScript™ III Reverse Transcriptase (Invitrogen) according to the manufacturer’s instructions.

The expression of circRNAs was measured using quantitative polymerase chain reaction (qPCR) SYBR Green Master Mix (Takara, Tokyo, Japan) in a ViiA 7 Real-time PCR System (Applied Biosystems Inc, Foster City, CA, United States). The sequences of the divergent primers for the detection of hsa_circ_0000745 by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) were 5’-GTTGAAAGTAGCCCGAGCAG-3’ and 5’-ACGTGGCACAGACCTCTCTC-3’. The primer sequences for glyceraldehyde-3-phosphate dehydrogenase (GAPDH; used as a control) were 5’-CTCGCTTCGGCAGCACA-3’ and 5’-AACGCTTCACGAATTTGCGT-3’. These primers were synthesized by Invitrogen (Shanghai, China). The reaction conditions were as follows: 95 °C for 10 min, and 40 cycles of 95 °C for 10 s and 60 °C for 60 s. The RNA levels were normalized to human GAPDH. The expression levels were analyzed by the 2^-ΔΔCt method.

All of the quantitative PCR reactions were conducted in triplicate. The appearance of a single-peak in the melt-curve suggested the specificity of the PCR products.

Statistical analysis

All experimental data were analyzed using SPSS software (version 22.0; IBM, Armonk, NY, United States) and GraghPad Prism 5.0 (GraphPad Software, La Jolla, CA, United States). The expression level of each circRNA was represented as fold-change using the 2^-ΔΔCt method. Differences of hsa_circ_0000745 levels between GC tissues and adjacent non-tumor tissues, and between plasma samples from patients with GC and healthy controls were assessed using t-test. The correlation between hsa_circ_0000745 levels and clinicopathological factors was further analyzed by one-way analysis of variance. The receiver operating characteristic (ROC) curve was established to evaluate the diagnostic value. The cut-off value of hsa_circ_0000745 was analyzed with SPSS software. A P value < 0.05 was considered statistically significant.

RESULTS

Validation of differentially expressed circRNAs by qRT-PCR

Among the three differentially expressed circRNAs examined (hsa_circ_0000745, hsa_circ_0085616, and chr16:30740286-30740893|+), hsa_circ_0085616 was reported to be up-regulated in colorectal tissues compared to that in normal tissues[20], while hsa_circ_0000745 and chr16:30740286-30740893|+ were found to be down-regulated in GC tissues compared to paired adjacent non-tumor tissues[21]. The observed changes in expression levels were validated by qRT-PCR analysis using 20 sets of gastric tissues and paired adjacent non-tumor tissues.

Divergent primers, rather than the more commonly used convergent primers, were designed for quantitative amplification of the circRNAs (Figure 1 and Table 2). The qRT-PCR results showed that circ_0085616 was up-regulated, while circ_0000745 and chr16:30740286-30740893|+ were down-regulated in the GC tissues compared to the paired adjacent non-tumor tissues (Figure 2), suggesting that the general expression patterns of the three circRNAs in GC were consistent with those in colorectal cancer.

Figure 1 Schematic diagram of the polymerase chain reaction primers used to specifically detect circular transcripts. Divergent polymerase chain reaction primers, rather than the more commonly used convergent primers, were designed for the circRNAs.

Table 2 Primers used for quantitative reverse transcription-polymerase chain reaction analysis of circular RNA and mRNA levels.

Target ID	Primer sequence, 5’-3’	Product size in bp
Circ_0000745	F: GTTGAAAGTAGCCCGAGCAG	204
	R: ACGTGGCACAGACCTCTCTC
Circ_0085616	F: GACTACAACTCGCCCACCAC	200
	R: TCCATTTCTGGGCCATAATC
Chr16:30740286-30740893|+	F: CCTTTGCACCGTATTGTGTG	199
	R: GGCAGGGTACAGAAATCCAG
GAPDH	F: CTCGCTTCGGCAGCACA	122
	R: AACGCTTCACGAATTTGCGT

bp: Base pairs; F: Forward; R: Reverse.

Figure 2 Validation of the expression levels of the three circRNAs using quantitative reverse transcription-polymerase chain reaction analysis. The three differently expressed circRNAs were validated by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) in 20 gastric tumor tissues and paired adjacent non-tumor tissues. The qRT-PCR results are presented as mean and SE of the mean.

CircRNA binding miRNA prediction by sequence analysis

Given that miRNAs play important roles in the progression of GC, some circRNAs could be involved in GC, likely via interacting with miRNAs. We analyzed the potential binding miRNAs for the three circRNAs using sequence analysis. The association of miRNAs with GC indicated that circRNAs may have a regulatory role in GC. A tree diagram of circRNAs and their potential binding miRNAs is shown in Figure 3. Based on these potential circRNA/miRNA interactions, hsa_circ_0000745 was predicted to be able to bind a spectrum of miRNAs, with the known functions suggesting its potential role in the regulation of GC growth.

Figure 3 Network of circular RNAs and the predicted binding miRNAs.

Hsa_circ_0000745 expression is down-regulated in GC tissues and plasma samples from patients with GC

We tested hsa_circ_0000745 levels in GC tissues and plasma samples from patients with GC. As shown in Figure 4A, hsa_circ_0000745 expression was found to be down-regulated in the GC tissues (P < 0.001). Furthermore, hsa_circ_0000745 levels in the plasma samples of patients with GC were lower than those in plasma samples from healthy controls (P < 0.001; Figure 4B).

Figure 4 Hsa_circ_0000745 expression levels. A: Hsa_circ_0000745 expression levels in gastric cancer tissues and paired adjacent non-tumor tissues; B: Hsa_circ_0000745 expression levels in plasma samples from patients with gastric cancer and healthy volunteers. P < 0.001, n = 60. The expression is shown as the 2^{- ΔΔCt} value.

Potential diagnostic values of hsa_circ_0000745 in GC

Since we found that hsa_circ_0000745 expression levels were lower in GC tissues and patients’ plasma samples, we further analyzed their association with clinicopathological features of patients with GC. As shown in Table 3, hsa_circ_0000745 levels in the GC tissues were significantly related to tumor differentiation (P = 0.012). However, we found no association between its levels with other clinicopathological features, such as age, sex, tumor diameters, lymphatic metastasis, TNM stage and carcinoembryogenic antigen (CEA) levels.

Table 3 Relationship of hsa_circ_0000190 expression levels (2^-ΔΔCt) in gastric cancer tissues with clinicopathological factors of patients with gastric cancer.

Clinicopathological factor		Tissue hsa_circ_0000745		Plasma hsa_circ_0000745
	n	Mean ± SD	P value	mean ± SD	P value
Age (yr)			0.757		0.195
< 60	33	0.77 ± 0.83		0.90 ± 0.96
≥ 60	27	0.58 ± 0.60		0.69 ± 0.53
Sex			0.203		0.252
Male	43	0.59 ± 0.67		0.68 ± 0.60
Female	17	0.92 ± 0.85		1.13 ± 1.12
Diameter in cm			0.168		0.127
≥ 5	26	0.76 ± 0.77		0.87 ± 0.76
< 5	34	0.63 ± 0.71		0.76 ± 0.84
Differentiation			0.012a		0.087
Well	2	0.96 ± 1.10		0.89 ± 0.44
Moderate	25	0.65 ± 0.71		0.70 ± 0.66
Poor	33	0.70 ± 0.77		0.88 ± 0.92
Lymphatic metastasis			0.064		0.087
Present	44	0.71 ± 0.71		0.85 ± 0.78
Absent	16	0.61 ± 0.81		0.69 ± 0.88
TNM stage			0.056		0.046a
I-II	14	0.63 ± 0.81		0.60 ± 0.53
III-IV	46	0.70 ± 0.72		0.87 ± 0.87
CEA in ng/mL			0.077		0.130
≥ 5	10	0.79 ± 0.72		0.65 ± 0.59
< 5	50	0.66 ± 0.74		0.85 ± 0.78

^aindicated statistical significance. CEA: Carcinoembryogenic antigen.

Meanwhile, hsa_circ_0000745 levels in plasma from GC patients were significantly related to TNM stage (P = 0.046). Nevertheless, we did not find any association between hsa_circ_0000745 and the other clinicopathological features.

To evaluate the potential diagnostic value, a ROC curve was generated for hsa_circ_0000745 levels in plasma. We found that the area under the ROC curve (AUC) was 0.683 (Figure 5). The sensitivity and specificity were 0.855 and 0.45, respectively. When the expression level of plasma hsa_circ_0000745 was combined with CEA level, the AUC increased to 0.775, while the sensitivity and specificity changed to 0.800 and 0.633, respectively.

Figure 5 Receiver operating characteristic curve of plasma hsa_circ_0000745 and serum carcinoembryogenic antigen levels. CEA: Carcinoembryogenic antigen.

DISCUSSION

Gastric carcinogenesis is a complicated biological process that involves a wide array of molecular abnormalities. In recent decades, the molecular mechanisms of GC have been extensively investigated. Expression arrays for genes and noncoding RNAs, such as miRNAs and lncRNAs, have been applied and recognized as feasible and useful approaches to profiling the molecular signatures of GC[25,26].

Today, circRNAs are the newest type of noncoding RNAs, and their discovery expanded our knowledge of the complexity of noncoding RNAs. The circRNAs were originally thought to be by-products of splicing errors or the mRNA process[27]. They represent a stable, diverse and conserved class of RNA molecules[11,28]. Ongoing research is continuing to provide evidence that circRNAs are involved in the development and progression of diseases, especially cancer[29].

Hsa_circ_002059 was first found to be significantly down-regulated in GC and proposed as a potential novel and stable biomarker for the diagnosis of this particular cancer type[19]. According to other reports, expression of the circular RNA-ITCH is typically down-regulated in esophageal squamous cell carcinoma (ESCC) and colorectal carcinoma, and may exert an inhibitory effect on ESCC and colorectal cancer by suppressing the Wnt/β-catenin pathway[30,31]. Collectively, these findings have suggested that circRNAs play a crucial role in cancers, and that they may serve as the basis for developing new strategies for diagnosis and therapy. However, the tumorigenesis mechanisms of GC are far from being completely elucidated and this idea is still in its infancy.

In this study, we randomly selected three differently expressed circRNAs, which we characterized using RNA-Seq analysis of ribosomal RNA-depleted total RNA from three GC or colorectal cancer tumor tissues and paired normal tissues. Circ_0085616 was up-regulated in colorectal cancer tissues compared with non-cancerous tissues, while circ_0000745 and chr16:30740286-30740893|+ were down-regulated in GC tissues compared with normal tissues. We further confirmed the RNA-Seq results of three circRNAs in 20 paired normal and cancerous gastric tissues by qRT-PCR analysis; circ_0085616 up-regulation and circ_0000745 and chr16:30740286-30740893|+ down-regulation in the GC tissues suggested that the circRNA sequencing results of GC are consistent with the expression patterns observed in colorectal cancer.

Recent studies have demonstrated that circRNAs can function as miRNA sponges or potent competitive endogenous RNA molecules[27,32-35]. Given that miRNAs play important roles in the progression of GC, it is possible and even likely that some circRNAs are involved in GC via interacting with miRNAs. We analyzed the potential binding miRNAs of the three circRNAs by sequence analysis. The association of miRNAs with GC indicated that circRNAs may have a regulatory role in GC. Based on the potential circRNA/miRNA interactions identified, hsa_circ_0000745 is potentially able to bind various miRNAs that may play critical roles in regulating GC growth and migration. However, additional studies are required to clarify whether these circRNAs are truly involved in GC, and whether circRNAs acting as miRNA sponges are a general phenomenon.

In this study, we found that hsa_circ_0000745 was down-regulated in GC tissues as well as in the plasma samples from patients with GC. It was suggested that hsa_circ_0000745 may be involved in other diseases besides GC, including various cancers (i.e., breast cancer, colorectal cancer, hepatocellular carcinoma, lung cancer, ovarian cancer, and rhabdomyosarcoma) and cerebellar neurodegeneration, Waldenstrom macroglobulinemia, cardiac hypertrophy, cardiomyopathy, coronary artery disease, heart failure, myeloproliferative disorder, myocardial infarction, retinitis pigmentosa, and schizophrenia, according to the previous analysis (http://gyanxet-beta.com/circdb/). That is the main reason why we choose hsa_circ_0000745 for further study.

We found that expression levels of hsa_circ_0000745 in GC tissues were associated with tumor differentiation. However, its expression did not show any statistically significant relationship with tumor size, lymphatic metastasis or TNM stage (Table 3). Moreover, we determined its potential diagnostic value, reported herein for the first time. When compared with serum CEA levels, hsa_circ_0000745 had a higher sensitivity and specificity in the screening of GC (Figure 5). The combination of plasma hsa_circ_0000745 and serum CEA levels, however, increased the diagnostic value. In summary, our study revealed that the circular RNA hsa_circ_0000745 may serve as a diagnostic marker for GC.

ACKNOWLEDGMENTS

We are grateful to the members of our laboratory for helpful discussions regarding this study. We thank Cloud-Seq Biotech Ltd. Co. (Shanghai, China) for the transcriptome sequencing service and the subsequent bioinformatics analysis.

COMMENTS

Background

Circular RNA (circRNA) is a recently identified type of noncoding RNA that plays an important role in the biological processes of tumor development and progression. The stability of circRNA makes it an attractive candidate as a serological diagnostic marker for cancer. However, the circRNA expression profiles in gastric cancer (GC) have not been elucidated.

Research frontiers

Recent studies have shown that different tumors have distinctive differential expression profiles of various circRNAs, such as hsa_circ_002059 in GC and hsa_circ_0004277 in acute myeloid leukemia.

Innovations and breakthroughs

This is the first study to demonstrate that hsa_circ_0000745 is differentially expressed in plasma samples of GC patients and normal healthy controls, as well as in gastric tumor tissues and paired adjacent non-tumor tissues. Moreover, the data suggest the possibility that hsa_circ_0000745 may serve as a novel diagnostic marker for GC.

Applications

Expression of hsa_circ_0000745 was found to correlate with tumor differentiation and tumor-node-metastasis stage. When plasma hsa_circ_0000745 expression was combined with serum level of carcinoembryogenic antigen, the diagnostic ability (especially specificity) for GC was improved.

Terminology

Circular RNA is a type of RNA which, unlike the better known linear RNA, forms a covalently closed continuous loop. GC is also known as stomach cancer, developing from the lining of the stomach. Diagnostic biomarkers are found in blood, urine or tissues and are abnormally expressed, in a trackable manner, in the presence of diseases.

Peer-review

This is the first study to show the promising diagnostic value of measuring circulating hsa_circ_0000745 levels in GC.