- Open Access
Circular RNA hsa_circ_0004872 inhibits gastric cancer progression via the miR-224/Smad4/ADAR1 successive regulatory circuit
Molecular Cancer volume 19, Article number: 157 (2020)
Emerging evidence has shown that circular RNAs (circRNAs) play a crucial regulatory role in the occurrence and development of cancer. Exploring the roles and mechanisms of circRNAs in tumorigenesis and progression may help to identify new diagnostic markers and therapeutic targets. In the present study, we investigated the role and regulatory mechanism of hsa_circ_0004872 in gastric cancer (GC).
qRT-PCR was used to determine the expression of hsa_circ_0004872 in GC tissues and cells. EdU, CCK-8, transwell and scratch wound healing assays were used to assess the role of hsa_circ_0004872 in GC cell proliferation, invasion and migration, respectively. Subcutaneous and tail vein tumor injections in nude mice were used to assess the role of hsa_circ_0004872 in vivo. RIP assay, biotin-coupled probe pull-down assay, FISH and luciferase reporter assay were performed to confirm the relationship between hsa_circ_0004872 and the identified miRNA. ChIP assay, luciferase reporter assay and western blot were used to determine the direct binding of Smad4 to the promoter of the ADAR1 gene.
In this study, we found that hsa_circ_0004872 was dramatically downregulated in GC tissues compared with adjacent noncancerous tissues. The expression level of hsa_circ_0004872 was associated with tumor size and local lymph node metastasis. Enforced expression of hsa_circ_0004872 inhibited the proliferation, invasion and migration of GC cells, whereas knockdown of hsa_circ_0004872 had the opposite effects. Nude mice experiments showed that ectopic expression of hsa_circ_0004872 dramatically inhibited tumor growth and metastasis in vivo. Moreover, we demonstrated that hsa_circ_0004872 acted as a “molecular sponge” for miR-224 to upregulate the expression of the miR-224 downstream targets p21 and Smad4. Importantly, we found that the RNA-editing enzyme ADAR1 inhibited hsa_circ_0004872 expression and further led to the upregulation of miR-224. Smad4, the downstream target of miR-224, could further affect hsa_circ_0004872 levels by directly binding to the promoter region of ADAR1 to inhibit ADAR1 expression.
Our findings showed that hsa_circ_0004872 acted as a tumor suppressor in GC by forming a negative regulatory loop consisting of hsa_circ_0004872/miR-224/Smad4/ADAR1. Thus, hsa_circ_0004872 may serve as a potential biomarker and therapeutic target for GC.
Gastric cancer (GC) is one of the most common malignancies in the world and ranks third in terms of cancer-related deaths . Despite improvements in clinical diagnosis and therapeutic strategies, the prognosis of GC patients remains poor due to a high recurrence rate and distant metastasis [2, 3]. Therefore, understanding the molecular mechanism of GC progression and identifying effective molecular diagnostic markers and therapeutic targets are extremely important to improve the survival rate of GC patients.
Circular RNA (circRNA) is a class of single-stranded RNAs with development/ tissue-specific expression patterns in eukaryotic cells . CircRNAs are generated from the ‘back-splicing’ of precursor mRNA (pre-mRNA) transcripts to form a covalently closed circular structure without 5′ to 3′ polarity or 3′ polyA tail . The unique circular structure of circRNAs confers inherent resistance to degradation by exonuclease and makes them more stable than their linear parental genes [4,5,6]. In addition, the expression of circRNA is evolutionarily conserved, and the expression levels of some circRNAs are much higher than those of their parental genes . Recently, evidence has shown that circRNAs are aberrantly expressed in tumor tissues and cells and play indispensable roles during the occurrence and development of various tumors [8, 9], such as hepatocellular carcinoma (HCC) [10, 11], lung cancer , GC , breast cancer , bladder cancer , and glioblastoma . Therefore, circRNAs have great potential as valuable diagnostic biomarkers or therapeutic targets in tumors . The function of circRNAs is usually associated with their source and location in the cells. Most circRNAs are derived from the exon(s) of genes, called exonic circRNAs (ecRNAs), which are located in the cytoplasm  and usually exert their biological functions at the posttranscriptional level by acting as miRNA “sponges”  or binding to RNA-binding proteins [20, 21]. The other two kinds of circRNAs, intron circRNAs (ciRNAs) and exon-intron circular RNAs (EIciRNAs), are mostly located in the nucleus and play a regulatory role in the transcription of their parental genes [22, 23] or act as templates for protein translation [24,25,26].
In this study, we used next-generation sequencing of circRNAs to screen the differentially expressed circRNAs in GC and adjacent nontumor tissues and identified the downregulated hsa_circ_0004872 (also known as circMAPK1 [27, 28]), which is derived from the exons of the MAPK1 gene. Then, we investigated the role and regulatory mechanism of hsa_circ_0004872 in GC. We found that hsa_circ_0004872 inhibited GC growth, invasion and metastasis both in vitro and in vivo by acting as a “sponge” of miR-224 to upregulate p21 and Smad4 expression. Importantly, we verified that Smad4 could further regulate hsa_circ_0004872 expression by directly binding to the promoter region of the RNA-editing enzyme ADAR1 and inhibiting ADAR1 expression in GC. Therefore, hsa_circ_0004872 exerted its functions via the miR-224/Smad4/ADAR1/hsa_circ_0004872 negative regulatory circuit. This is the first report on the expression, role and regulatory mechanism of hsa_circ_0004872 in GC. Our findings indicate that hsa_circ_0004872 may act as a promising diagnostic and prognostic marker for GC patients.
Materials and methods
Patient samples and ethical approval
We collected 76 pairs of GC tissues and corresponding adjacent noncancerous tissues from Shandong Cancer Hospital in 2012–2013 and Taian City Central Hospital in 2016–2017. None of the patients enrolled in the study received chemotherapy or radiotherapy before surgery, and there was no evidence of any other malignancies. For each tumor, age, sex, tumor size, invasiveness and regional and distant metastases were evaluated and described in supplementary Table S1. The samples were excised from patients, immediately frozen in liquid nitrogen and stored until use. All samples were collected for research use only. The research was approved by the ethics committee of the School of Basic Medical Sciences, Shandong University.
RNA sequencing analysis
The total RNA was isolated from the GC and adjacent nontumor tissue using TRIzol reagent (Invitrogen, USA) according to the manufacturer’s instructions. The RNA purity and concentration were determined with NanoDrop ND-1000 (NanoDrop Thermo). The RNA integrity of the samples was assesed by denaturing agarose gel electrophoresis.The Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, USA) was used to remove the rRNA. The high-throughput whole transcriptome sequencing and subsequent bioinformatics analysis were performed by Cloud-Seq Biotech (Shanghai, China) as previously reported . Paired-end reads were acquired from the Illumina HiSeq 4000 sequencer. Edger software (v3.16.5) was used to identify the differentially expressed circRNAs.
Cell culture and transfection
The human GC cell lines AGS, BGC-823, SGC-7901 and gastric epithelial immortalized GES-1 cells were purchased from the Cell Resource Center, Institute of Biochemistry and Cell Biology at the Chinese Academy of Science (Shanghai, China). All cells were cultured in RPMI-1640 medium (GIBCO Thermo Fisher Scientific, USA) with 10% fetal bovine serum (Biological Industries, USA) and 100 U/mL penicillin-streptomycin (Life Technologies, USA). All cells were cultured at 37 °C with 5% CO2.
Plasmids were transfected into GC cells using X-tremeGENE HP Transfection Reagent (Roche Applied Science, Germany). SiRNAs, miRNA mimics or inhibitors were transfected into cells using Lipofectamine 2000 (Invitrogen, USA). The experiment was performed according to the manufacturer’s instructions. To construct the hsa_circ_0004872 stably overexpressing cell lines, we transfected the control vector pLCDH-circ or the hsa_circ_0004872 overexpression vector pLCDH-hsa_circ_0004872 into BGC-823 cells and then selected them with puromycin (Sigma, USA) for 2–3 weeks until hsa_circ_0004872 was stably overexpressed in the cells.
miRNAs, siRNAs and plasmid construction
The hsa_circ_0004872 overexpression vector pLCDH-hsa_circ_0004872 was constructed by BioSune (Shanghai, China) by inserting the sequence of hsa_circ_0004872 into the pLCDH-circ expression vector. The mutated vector was constructed by changing the binding site of miR-224 in hsa_circ_0004872“GTGACTT” into “ACAGTCC” and inserting the mutated sequence into the pLCDH-circ expression vector. SiRNAs targeting hsa_circ_0004872 and the control siRNA were synthesized by Genepharma (Shanghai, China). All the siRNA sequences are listed in Table S2. Mimics and inhibitor of miR-224 were synthesized by RiboBio (Guangzhou, China). The pMIRGLO-hsa_circ_0004872 luciferase reporter plasmid was constructed by inserting the sequence of hsa_circ_0004872 into pMIRGLO (Promega, USA) between the Xho I and Sal I sites. The pMIR-p21 and pMIR-Smad4 luciferase reporter plasmids were constructed by inserting the 3’UTR fragment of p21 or Smad4 into the pMIR reporter vector (Promega, USA) between the Spe I and Hind III sites. The miR-224 complementary sequence “GTGACTT” in hsa_circ_0004872 and the 3′UTRs of p21 and Smad4 were mutated to remove the complementarity. The pGL3-ADAR1 luciferase reporter plasmid was constructed by inserting the promoter sequence of ADAR1 (+ 12 to − 1999) into the pGL3-basic vector (Promega, USA) between the Mlu I and Hind III sites. All of the primer sequences are listed in Table S3. The p-3 × flag-CMV-Smad4 plasmid and control plasmid (p-3 × flag-CMV) were kindly provided by Assistant researcher Lidong Zu (Center of Pathology, School of Medicine, Shanghai Jiao Tong University, China). pmGFP-ADAR1-p150 (Plasmid #117927) and pmGFP-ADAR1-p110 (Plasmid #117928) were purchased from Addgene (Cambridge, USA). All of the constructs were verified by sequencing.
RNA extraction, qRT-PCR, RNase R treatment and actinomycin D treatment
The total RNA was isolated from tissues and cells using TRIzol reagent (Invitrogen, USA) according to the manufacturer’s protocol. cDNA was synthesized with random primers or miRNA-specific stem-loop primers using a Revert Aid First Strand cDNA Synthesis kit (Thermo Scientific, USA). qRT-PCR was performed on the Bio-Rad CFX96TM Real-Time PCR System (Bio-Rad). Primers for mRNAs were synthesized by BioSune (Shanghai, China), and miRNA primers were synthesized by RiboBio (Guangzhou, China). The relative RNA expression levels were analyzed using the 2-△△Ct method. β2-M and U6 were used as the internal control genes for mRNA and miRNA, respectively. For the RNase R treatment, 2 μg total RNA was digested at 37 °C for 20 min and 70 °C for 5 min with 3 U/μg RNase R (Epicentre Technologies, USA). Then, the expression levels of linear mRNA and circular RNA were determined by qRT-PCR, RT-PCR and Northern blot. For actinomycin D (ActD) treatment, 2 μg/mLActD (Sigma-Aldrich, USA) for BGC-823 cells and 1 μg/mL ActD for SGC-7901 cells were used. Cells were collected in a series of time intervals, and the expression of linear mRNA and circRNA was detected by qRT-PCR.
Total proteins from GC cells were extracted with RIPA lysis buffer containing proteinase inhibitor. The protein concentrations were measured by the BCA reagent kit (Beyotime, China). The proteins were separated by SDS-PAGE and transferred to polyvinylidene difluoride membranes (Merck Millipore, Germany), which were blocked in 5% nonfat milk and then incubated with primary antibodies against p21 (1:1000, Proteintech, 10355-1-AP, China), Smad4 (1:2500, Abcam, ab40759, USA), HOXD10 (1:5000, Abcam, ab138508, USA), PTEN (1:1000, Cell Signaling Technology, #9559, USA), ADAR1 (1:4000, Abcam, ab126745, USA), myc (1:2000, Proteintech, 16286-1-AP, China), β-actin (1:4000, Cell Signaling Technology, #12262, USA) and GAPDH (1:3000, Proteintech, 10494-1-AP, China) at 4 °C overnight. Then, the membranes were washed in Tris-buffered saline with Tween, incubated with anti-mouse or anti-rabbit horseradish peroxidase-conjugated secondary antibody for 1 h at room temperature, and developed with the enhanced chemiluminescence method (Millipore, USA). β-Actin or GAPDH served as a loading control.
IHC was performed by using the kit PV-9000 (ZSGB-BIO, China) according to the manufacturer’s protocol. Paraffin sections were incubated with primary antibodies against p21 (1:300, Proteintech, 10355-1-AP, China) or Smad4 (1:100, Abcam, ab40759, USA) at 4 °C overnight. After washing in PBS, the sections were incubated with anti-mouse or rabbit horseradish peroxidase-conjugated secondary antibodies for 1 h at room temperature. Then paraffin sections were stained with DAB and hematoxylin. Finally, the paraffin sections were covered with coverslips for microscopic (Nikon, Japan) observation.
RNA fluorescence in situ hybridization (FISH)
The Cy3-labeled hsa_circ_0004872 probes and FITC-labeled hsa-miR-224 probes were designed and synthesized by Geneseed Biotechnology Co., Ltd. (Guangzhou, China).The sequences of the probes are listed in Table S4. BGC-823 cells were grown on round coverslips, fixed, permeabilized in PBS with 0.5% Triton X-100 and dehydrated in ethanol. FISH probes were diluted (1:50), denatured, balanced and added to cells at 37 °C overnight. After hybridization, cells were stained with DAPI-Antifade for 10 min at room temperature, and the slides were sealed with rubber cement and placed in the dark for more than 20 min. Finally, the results were observed with a TCS SP2 AOBS laser confocal microscope (Leica Microsystems, Germany).
RNA-binding protein immunoprecipitation (RIP) assay
The RIP assay was performed using the EZ-Magna RIP RNA-binding protein immunoprecipitation kit 17-701 (Merck Millipore, Germany) according to the manufacturer’s instructions. Briefly, BGC-823 cells were collected and lysed in RIP lysis buffer with protease and RNase inhibitors. The cell lysates were incubated with magnetic beads conjugated with Ago2 antibodies or IgG (Millipore) at 4 °C overnight. Then, the beads were washed and incubated with proteinase K to remove proteins. Finally, RNA was extracted and subjected to RT-PCR and agarose gel electrophoresis analysis.
Biotin-coupled probe pull-down assay
The biotin-labeled hsa_circ_0004872 probe and the negative control probe were synthesized by RiboBio Biotech (Guangzhou, China), and the sequences are listed in Table S4. The circRIP assay was performed as previously described  with minor modifications. Briefly, hsa_circ_0004872-overexpressing BGC-823 cells were fixed with 1% formaldehyde for 10 min at 37 °C, and then the cells were lysed and sonicated. After centrifugation, 50 μL of the supernatant was used as input, and the remaining part was incubated with the biotin-labeled probe against hsa_circ_0004872 or negative control for 4 h at 25 °C. Then, hsa_circ_0004872 miRNAs complexes were captured with streptavidin-coupled Dynabeads (Invitrogen) overnight. On the next day, the hsa_circ_0004872/miRNAs/bead complexes were incubated with RIP wash buffer containing proteinase K. Finally, TRIzol was added to the mixture for RNA extraction and detection.
The probe which spanned the back-splice junction of hsa_circ_0004872 was labeled with Digoxin. The probe sequences were listed in Table S4. Twenty micrograms total RNA was denatured in formaldehyde and then electrophoresed in a 1% agarose-formaldehyde gel. Then the RNAs were transferred onto a Hybond-N+ nylon membrane (Millipore).The membranes were crosslinked, pre-hybridized and hybridized with DIG-labeled DNA probes overnight. After washing, the membranes was incubated with alkaline phosphatase (AP)-conjugated anti-DIG antibodies (Roche). Immunoreactive bands were visualized after adding the chemiluminescent substrate CSPD (Roche) followed by exposure in the exposure apparatus (Tanon-520, Tianneng, China).
Luciferase reporter assay
BGC-823 and SGC-7901 cells were seeded in 24-well plates (3 × 104 cells/well). Negative control mimics or miR-224 mimics were cotransfected with the reporter plasmid into GC cells using Lipofectamine 2000. After 48 h, the cells were collected and lysed with passive lysis buffer. The luciferase activities were assessed using the Dual Luciferase Assay Kit (Promega) according to the manufacturer’s protocol. Firefly luciferase activity was normalized to Renilla luciferase activity.
Chromatin immunoprecipitation (ChIP) assay
Treated GC cells were cross-linked in 1% (vol/vol) formaldehyde-containing medium for 10 min at 37 °C and then lysed and sonicated to disrupt chromatin DNA into fragments between 200 and 1000 bp. ChIP was performed using the ChIP Assay 17–295 kit (Merck Millipore, Germany) according to the manufacturer’s protocol. An antibody against Flag (Sigma-Aldrich, USA) was used to immunoprecipitate the DNA fragments. The protein-DNA complexes were collected with protein A Sepharose beads, eluted and reverse cross-linked. The samples were extracted with Dr.GenTLE™ Precipitation Carrier (TaKaRa Biotechnology Co., Ltd., Dalian, China). The recovered DNA was resuspended in DDW and used as the template to amplify the ADAR1 promoter. The primer sequences are listed in Table S3.
Cell proliferation assay
Cell proliferation was detected with EdU and CCK-8 assays. The EdU assay was performed according to the protocol of the Cell-Light™ EdU Apollo®567 In Vitro Imaging Kit (RiboBio, Guangzhou, China). Briefly, the treated cells were seeded in 96-well plates and incubated with 50 μM EdU for 2 h at 37 °C. After being fixed with 4% paraformaldehyde, the cells were exposed to 100 μL of 1 × Apollo® reaction cocktail and then incubated with 5 μg/mL Hoechst 33342 to stain cell nuclei. Images were captured using a fluorescence microscope (Olympus, Tokyo, Japan). The percentage of EdU-positive cells was defined as the proliferation rate. For the CCK-8 assay, the treated cells were seeded in 96-well plates and incubated with 100 μL 10% CCK-8 solution for 4 h at 37 °C. The absorbance was measured at 450 nm with an Infinite M200 spectrophotometer (Tecan). All of the experiments were repeated three times in triplicate.
Scratch wound healing assay
GC cells with different transfections were cultured in six-well plates at 37 °C. Scratch wounds were created by using the fine end of 10-μL pipette tips. Images of migrated cells were captured under phase-contrast microscopy at different times.
Invasion and migration assays
The invasion and migration assays were performed with a Transwell chamber coated with Matrigel (BD Biosciences, USA) (for the invasion assay) or without Matrigel (for the migration assay). The transfected GC cells were collected and resuspended in serum-free RPMI-1640 medium. A total of 5 × 104 cells for the migration assay or 1 × 105 cells for the invasion assay were seeded into the upper compartment of a 24-well chamber with an 8.0 μm pore (Corning, USA). RPMI-1640 medium containing 20% FBS was added to the lower chambers as a chemoattractant. The cells were incubated for 48 h for the invasion assay or 24 h for the migration assay. Then, the cells on the upper surface of the polycarbonate membrane were removed with cotton swabs, while the cells on the lower side were fixed with 100% methanol, stained with 0.05% crystal violet and imaged under a microscope. The number of migrating cells was counted from three independent experiments.
Xenograft nude mouse model
Five-week-old female BALB/c nude mice were purchased from the Nanjing Biomedical Research Institute of Nanjing University (Nanjing, China). For the tumor formation experiment, BGC-823 cells (4 × 105) that were stably transfected with the hsa_circ_0004872 overexpression vector (pLCDH-hsa_circ_0004872) or control vector (pLCDH_circ) were subcutaneously injected into either side of the back of each mouse. Tumor size was monitored by measuring the length (L) and width (W) of the tumor every 3 days with a caliper, and the tumor volume (V) was calculated with the formula V = 1/2 × L × W2. Twenty days after injection, the mice were euthanized, and the tumors were weighed.
To investigate the effect of hsa_circ_0004872 on the metastatic ability of GC cells in nude mice, GC cells (6 × 105) stably transfected with pLCDH-hsa_circ_0004872 or pLCDH_circ were injected into the tail veins of nude mice. After six weeks, the mice were sacrificed. The lungs of the nude mice were collected and stained with hematoxylin-eosin staining. All mice experiments were approved by the ethics committee of the School of Basic Medical Sciences, Shandong University and were performed according to the guidance of animal experiments in the Laboratory Animal Center of Shandong University.
hsa_circ_0004872 sequence data were obtained from circBase (http://www.circbase.org/). The target miRNAs of hsa_circ_0004872 were predicted with circular RNA interactome (https://circinteractome.nia.nih.gov), BIOINF (http://www.bioinf.com.cn/), and starBase (http://starbase.sysu.edu.cn/index.php). The binding sites of Smad4 in the ADAR1 promoter region were predicted with Jaspar (http://jaspar.genereg.net/).
Comparisons between different groups were analyzed with GraphPad Prism v7.0 software (GraphPad Software, La Jolla, CA, USA) using a paired or unpaired t-test. Linear regression analysis was used to analyze the correlation between hsa_circ_0004872 and miR-224 expression and between Smad4 and ADAR1 expression in GC samples. Overall survival comparison between the high ADAR1 group and the low Smad4 group was conducted by log-rank (Mantel–Cox) test in the Kaplan–Meier plots. p < 0.05 was considered statistically significant.
hsa_circ_0004872 is downregulated in GC tissues
To explore circRNAs expression profile in GC, we performed circRNA-seq analysis of GC tissue and corresponding nontumor tissue from Shanghai Cloud-seq Biotech Co.Ltd. (Shanghai, China).The differentially expressed circRNAs were listed in Table S5. Then, we chose the most differentially expressed circRNAs, hsa_circ_0004872, hsa_circ_0002483, hsa_circ_0000847 and hsa_circ_0001566, to be further detected with qRT-PCR in 42 paired GC and adjacent nontumor tissues. We found that hsa_circ_0004872 showed the most significant downregulation (36/42 × 100% = 85.7%, p < 0.0001) in GC tissues (Fig. S1). Therefore, we chose hsa_circ_0004872 as our candidate for further experiments.
hsa_circ_0004872 contains 490 nucleotides and is derived from the exon 2, 3 and 4 of human MAPK1 in chr22 (22153300-22162135) (Fig. 1a). It is also known as circMAPK1 [27, 28]. To verify the qRT-PCR results of hsa_circ_0004872 in GC tissues, we performed Sanger sequencing on the PCR products of hsa_circ_0004872. The result confirmed the head-to-tail splicing in the PCR products with the expected size and predicted splicing site (Fig. 1b). Given that head-to-tail splicing is not only the result of reverse splicing of cDNA but also the result of genomic rearrangements, we designed convergent primers and divergent primers for hsa_circ_0004872 to perform PCR using cDNA or genomic DNA from BGC-823 cells as the template. The results indicated that hsa_circ_0004872 was amplified by the divergent primers only from cDNA but not from gDNA (Fig. 1c). We used the divergent primers to further detect the expression of hsa_circ_0004872 in another 34 paired GC and nontumor tissues. The statistical analysis from the 34 pairs of samples and the first 42 pairs of samples showed that hsa_circ_0004872 was downregulated in 86.8% of the GC tissues (66/76 × 100% = 86.8%), and the expression of hsa_circ_0004872 was significantly decreased in GC tissues (p < 0.0001) (Fig. 1d). The clinical and pathological characteristics of these patients are shown in Table S1. Statistical analysis revealed that low hsa_circ_0004872 expression levels were correlated with advanced T stage and N stage and big tumor size in GC patients (Table 1). It has been reported that stability is one of the crucial characteristics of circRNAs [5,6,7]. To confirm the stability of hsa_circ_0004872, we used RNase R to treat the RNA extracted from BGC-823 and SGC-7901 cells. qRT-PCR and RT-PCR results showed that the levels of linear MAPK1 decreased sharply under RNase R treatment, but RNase R failed to degrade hsa_circ_0004872 (Fig. 1e and f). Northern blot result also verified the stability of hsa_circ_0004872 (Fig. 1g). To further confirm the circular characteristics of hsa_circ_0004872, we used the transcription inhibitor ActD to treat BGC-823 and SGC-7901 cells for different times. As shown in Fig. 1h, the mRNA level of linear MAPK1 decreased gradually with time, while the level of hsa_circ_0004872 was more stable and resistant to ActD treatment.
hsa_circ_0004872 inhibits the proliferation, invasion and migration of GC cells
To investigate the biological role of hsa_circ_0004872 in GC cells, we constructed the hsa_circ_0004872 overexpression vector pLCDH-hsa_circ_0004872 and transfected the overexpression vector or the control vector (pLCDH-circ) into the GC cells BGC-823 and SGC-7901. The transfection efficiency was validated by Northern blot and qRT-PCR (Fig. 2a and Fig. S2A). Then, we used EdU and CCK-8 assays to detect the proliferation ability and used the scratch wound healing and transwell assays to detect the invasion and migration abilities of the GC cells. The results showed that enforced expression of hsa_circ_0004872 significantly suppressed both the proliferation ability (Fig. 2b-d) and invasion and migration ability of GC cells (Fig. 2e-h).
To further determine whether hsa_circ_0004872 knockdown has the opposite effects, we designed two siRNAs to specifically target the junction site of hsa_circ_0004872 (Fig. S3 A). We determined the expression level of hsa_circ_0004872 in different GC cells and found that the expression of hsa_circ_0004872 was relatively low in SGC-7901 cells (Fig. S4). Therefore, we only used BGC-823 cells to perform the knockdown experiment. We transfected the siRNAs into BGC-823 cells. The qRT-PCR results showed that the expression of hsa_circ_0004872 was significantly decreased in hsa_circ_0004872 siRNA-transfected cells, while the expression level of the parental gene MAPK1 showed no obvious changes (Fig. S3B). Functionally, the EdU assay revealed that the proliferation of BGC-823 cells increased in the hsa_circ_0004872 siRNA group compared with the control group (p < 0.01) (Fig. S3C and D). Moreover, the scratch wound healing assay and transwell assay showed that hsa_circ_0004872 siRNAs significantly promoted the invasion and migration ability of GC cells (Fig. S3 E-H).
hsa_circ_0004872 acts as a “molecular sponge” for miR-224
Since the functions of circRNAs are usually related to their different localization in cells , we detected the subcellular localization of hsa_circ_0004872 in GC cells via the FISH assay. The results indicated that hsa_circ_0004872 mainly existed in the cytoplasm of GC cells (Fig. 3a). It has been reported that circRNAs in the cytoplasm usually act as miRNA “molecular sponges” . Therefore, we performed a RIP assay with an AGO2 antibody to explore whether hsa_circ_0004872 could bind to miRNAs in GC cells. As shown in Fig. 3b, hsa_circ_0004872 could be amplified from the immunoprecipitate pulled down by the AGO2 antibody, suggesting that hsa_circ_0004872 binds with miRNAs via AGO2.
Then, we used different databases (circular RNA interactome, BIOINF and starBase) to predict the potential miRNAs that bind with hsa_circ_0004872 (Table S6) and found 4 miRNAs (hsa-miR-224, hsa-miR-654-3p, hsa-miR-513a-5p and hsa-miR-599) from the overlap of the three databases (Fig. 3c). Among these miRNAs, we found that only hsa-miR-224 (miR-224) exerts oncogenic functions according to previous reports. Therefore, we performed a pull-down assay with a biotin-coupled hsa_circ_0004872 probe to determine whether hsa_circ_0004872 could interact with miR-224. The results showed that the biotin-coupled hsa_circ_0004872 probe group exhibited an obvious enrichment of miR-224 compared with the negative control group (Fig. 3d and e), suggesting that hsa_circ_0004872 could directly interact with miR-224. FISH assays showed that hsa_circ_0004872 colocalized with miR-224 in the cytoplasm of GC cells (Fig. 3f). qRT-PCR results showed that hsa_circ_0004872 knockdown increased miR-224 levels (Fig.3g) and hsa_circ_0004872 overexpression suppressed miR-224 expression, while miR-224 binding site mutated hsa-circ_0004872 overexpression has no obvious effect on the expression of miR-224 (Fig. S2B and Fig. 3h). To further determine whether miR-224 could directly target hsa_circ_0004872 in GC cells, we performed a dual luciferase assay. We constructed an hsa_circ_0004872 wild-type luciferase reporter vector (pMIRGLO-hsa_circ_ 0004872-WT) containing the binding sequences of miR-224 and a mutant luciferase reporter vector (pMIRGLO- hsa_circ_0004872-Mut) in which the binding sequences were mutated (Fig. S5A). BGC-823 and SGC-7901 cells were cotransfected with the WT or Mut luciferase reporter and miR-224 mimics. The results showed that hsa_circ_0004872 wild-type luciferase reporter activity was significantly suppressed by miR-224 mimics, while the activity of the Mut luciferase reporter was not affected by miR-224 mimics (Fig. 3i). These results suggested that miR-224 was a direct target of hsa_circ_0004872 in GC cells.
Additionally, we investigated the expression of miR-224 in paired GC and corresponding nontumor tissues and found that the expression level of miR-224 was significantly increased in GC tissues (p = 0.0098) (Fig. 3j and Fig. S6A). Statistical analysis showed that miR-224 expression was negatively associated with hsa_circ_0004872 in these GC samples (p = 0.0144) (Fig. 3k). Taken together, these results suggest that hsa_circ_0004872 functions as a “molecular sponge” for miR-224 in GC cells.
miR-224 promotes the proliferation, invasion and migration of GC cells by targeting the 3′-UTR of p21 and Smad4
Having confirmed the interaction of hsa_circ_0004872 and miR-224, we next sought to explore the biological role of miR-224 in gastric carcinogenesis and progression. We transfected the miR-224 mimics or inhibitor into GC cells. The transfection efficiency was verified by qRT-PCR (Fig. 4a; Fig. S7A). Then, EdU and CCK-8 assays revealed that miR-224 mimics promoted, whereas its inhibitor suppressed, GC cell proliferation (Fig. 4b-e; Fig. S7B-D). Moreover, scratch wound healing assays and transwell assays showed that miR-224 mimics promoted, whereas its inhibitor suppressed, cell invasion and migration (Fig. 4f-i; Fig. S7E-H).
Smad4 [30, 31], p21 , HOXD10  and PTEN  are validated targets of miR-224 in other tumors. To determine whether miR-224 can regulate these targets in GC cells, we transfected the miR-224 mimics or inhibitor into GC cells and used western blot and qRT-PCR to determine the expression levels of these genes. As shown in Fig. 4j, miR-224 mimics significantly decreased, whereas miR-224 inhibitor increased, the protein levels of p21 and Smad4, although miR-224 mimics or inhibitor had no obvious effect on the expression levels of HOXD10 and PTEN. qRT-PCR results further showed that miR-224 mimics significantly decreased, whereas miR-224 inhibitor increased, the mRNA levels of p21 and Smad4 (Fig. 4k, l).
To determine whether miR-224 could directly target the 3’UTRs of Smad4 and p21 in GC cells, we cloned the wild-type 3′-UTR sequences and the mutant sequences of p21 and Smad4 to construct wild-type luciferase reporter vectors (pMIR-Smad4(p21)-3’UTR-WT) and mutant vectors (Fig. S5B). GC cells were cotransfected with WT or Mut luciferase reporter plasmids and miR-224 mimics. Dual luciferase reporter assay results showed that miR-224 mimics strongly reduced the wild-type luciferase reporter activity of Smad4 and p21, while the mimics had no obvious effect on the Mut luciferase reporter activity of Smad4 and p21 (Fig. 4m, n). These results suggested that Smad4 and p21 were the direct targets of miR-224 in GC cells.
We further determined the effects of hsa_circ_0004872 on Smad4 and p21 expression and found that hsa_circ_0004872 knockdown reduced the expression of Smad4 and p21 in GC cells (Fig. 4o) and the enforced expression of hsa_circ_0004872 significantly increased the p21 and Smad4 protein levels, while mutated hsa_circ_0004872 has no obvious effect on the expression of p21 and Smad4. (Fig. 4p).
hsa_circ_0004872 exerts its antitumor effect via miR-224
To determine whether hsa_circ_0004872 functions as a tumor suppressor gene via miR-224, we performed rescue experiments. We cotransfected the hsa_circ_0004872 overexpression vector (pLCDH-circ_0004872) and miR-224 mimics into GC cells to determine whether the tumor-suppressing effect of hsa_circ_0004872 could be blocked by miR-224 mimics. The results showed that miR-224 mimics could partly attenuate the inhibition of cell proliferation (Fig. 5a-c), invasion and migration (Fig. 5d-g) mediated by hsa_circ_0004872 overexpression.
We next investigated whether hsa_circ_0004872-mediated upregulation of p21 and Smad4 can be abolished by miR-224 mimics. As shown in Fig. 5h, the upregulation of Smad4 and p21 mediated by hsa_circ_0004872 was partly rescued by miR-224 mimics (Fig. 5h). Therefore, these data demonstrated that hsa_circ_0004872 inhibited GC proliferation, invasion and migration, at least partly through miR-224.
hsa_circ_0004872 impaired tumorigenesis and metastasis in vivo
To investigate whether hsa_circ_0004872 suppressed the tumorigenesis and metastasis of GC cells in vivo, we performed subcutaneous injection and tail vein injection in nude mice. The subcutaneous injection results showed that the mean tumor volume in the hsa_circ_0004872 overexpression group was much smaller than that in the control group (Fig. 6a, b), and the average tumor weight in the hsa_circ_0004872 overexpression group was much lighter than that in the control group (Fig. 6c). HE staining showed that all the tumors were solid tumors (Fig. 6d). qRT-PCR analysis demonstrated that the expression of hsa_circ_0004872 was elevated (Fig. 6e), whereas miR-224 expression was decreased (Fig. 6f), in the tumors from the hsa_circ_0004872 overexpression group, which was consistent with what we observed in vitro. IHC results showed that the expression of p21 and Smad4 was increased in the tumors from the hsa_circ_0004872 overexpression group (Fig. 6g,h). In the tail vein injection assay, we observed that both the mice injected with the control group and the pLCDH-hsa_circ_0004872 group formed metastatic nodules at the surface of the lungs, but the lungs from mice injected with pLCDH-hsa_circ_0004872 cells were significantly smaller and lighter than those from mice injected with pLCDH-circ cells (Fig. 6i). Additionally, mice injected with pLCDH-hsa_circ_0004872 cells had much fewer metastatic nodules in the lungs than mice injected with pLCDH-circ cells (Fig. 6j, k). Taken together, these results demonstrate that hsa_circ_0004872 efficiently inhibited GC cells growth and metastasis in vivo.
hsa_circ_0004872 is regulated by the RNA-editing enzyme ADAR1
We next sought to explore which factors lead to the downregulation of hsa_circ_0004872 in GC. According to previous reports, the formation of circRNA is regulated by RNA binding proteins , such as MBl , QKI [36, 37], and ADAR1 [38, 39]. We analyzed the public NCBI GEO databases GSE27342 and GSE66229 and found that only ADAR1 expression was significantly higher in GC tissues than in nontumor tissues in both GSE27342 and GSE66229 (p < 0.0001) (Fig. S8A-F). Our qRT-PCR results also confirmed the higher expression level of ADAR1 in GC tissues compared with the adjacent nontumor tissues (Fig. 7a, Fig. S6B), which was consistent with the results from Dou et al. . Kaplan-Meier Plotter Database (201786_s_at) analysis of overall survival (Fig. 7b) and survival with lymph node metastasis (Fig. 7c) indicated that higher expression of ADAR1 in GC tissues was associated with shorter survival period. ADAR1 is an RNA editing enzyme that catalyzes the deamination of adenosine into inosine (A-to-I) and prevents the link between exon ends and circRNA formation [38, 41]. Thus, we wondered whether ADAR1 contributed to the downregulation of hsa_circ_0004872 in GC. Therefore, we transfected ADAR1 siRNA or the ADAR1 overexpressing vector (pmGFP- ADAR1-p110 or pmGFP-ADAR1-p150) into GC cells. The transfection efficiency was validated by western blot and qRT-PCR (Fig. 7d, e, g, h). Then, qRT-PCR was used to examine the expression of hsa_circ_0004872. The results showed that ADAR1 knockdown increased the expression of hsa_circ_0004872 (Fig. 7f) and that overexpression of ADAR1 led to the downregulation of hsa_circ_0004872 in GC cells (Fig. 7i). We further investigated whether ADAR1 can regulate miR-224 expression via hsa_circ_0004872. As shown in Fig. 7j, ADAR1 knockdown led to the downregulation of miR-224, and the downregulation of miR-224 mediated by ADAR1 siRNA was partly reversed by hsa_circ_0004872 siRNA. Therefore, these data demonstrated that ADAR1 regulated miR-224 expression via hsa_circ_0004872.
Smad4 regulates hsa_circ_0004872 expression via ADAR1
Given that hsa_circ_0004872 acts as a “molecular sponge” of miR-224 to upregulate the miR-224 target Smad4, we asked whether Smad4, as a transcription factor, can regulate ADAR1 to further lead to the dysregulation of hsa_circ_0004872. We used software to analyze the promoter of ADAR1 and found 5 potential binding sites of Smad4 at the region from nucleotide − 1993 to − 1981 (Site A), − 1959 to − 1947 (Site B), − 1590 to − 1578 (Site C), − 1280 to − 1268 (Site D) and − 760 to − 748 (Site E). Then, we performed a ChIP assay to determine the binding of Smad4 on the promoter of ADAR1. We transfected BGC-823 cells with the overexpression vector of Smad4 (p-3 × flag-CMV-Smad4) or the control vector (p-3 × flag-CMV) and used the flag antibody to precipitate the promoter fragment spanning these five sites. The results illustrated in Fig. 8a show a clear signal in Site D of the ADAR1 promoter in the cells transfected with p-3 × flag-CMV-Smad4.
To further determine whether Smad4 regulates ADAR1 and hsa_circ_0004872 expression, we transfected p-3 × flag-CMV-Smad4 or p-3 × flag-CMV-14 into GC cells. Western blot and qRT-PCR results showed that Smad4 overexpression led to a significant downregulation of the protein and mRNA levels of ADAR1 (Fig. 8b, c), while hsa_circ_0004872 was upregulated after transfection (Fig. 8d). We further constructed a luciferase reporter vector containing the fragment from − 1999 to + 12 in the ADAR1 promoter (pGL3-ADAR1-promoter-WT), which contained the sequence of Site D, and a mutant luciferase reporter vector (pGL3-ADAR1-promoter-Mut) in which the Site D sequence was mutated (Fig. S5C). BGC-823 cells were cotransfected with the WT or Mut luciferase reporter and the Smad4 overexpression vector. The luciferase reporter assay showed that Smad4 overexpression decreased the WT promoter activity of ADAR1, whereas it had no effect on the Site D Mut promoter activity of ADAR1 (Fig. 8e).
We next analyzed the expression of Smad4 in the public database NCBI GEO (GSE66229) and found that the expression of Smad4 was lower in GC tissues than in adjacent nontumor tissues (p < 0.0001) (Fig. 8f). Statistical analysis showed that there was a highly significant negative correlation between ADAR1 and Smad4 in these samples (p < 0.0001) (Fig. 8g). Kaplan-Meier Plotter Database (202527_s_at) analysis of overall survival (Fig. 8h) and survival with lymph node metastasis (Fig. 8i) indicated that lower expression of Smad4 in GC tissues resulted in a shorter survival period.
Circular RNA (circRNA) is a kind of single-stranded RNA that is widely found in eukaryotic cells . In contrast to linear RNAs, circRNAs form a covalently closed loop structure, which makes them more stable than their linear host genes. Recently, evidence has accumulated that circRNAs play important roles in different biological processes, especially in the tumorigenesis, invasion and metastasis [9,10,11,12,13,14,15,16]. However, the function and regulatory mechanism of circRNAs in GC remain largely unknown. In this study, we focused on hsa_circ_0004872, which originates from the back-splice of MAPK1 (mitogen-activated protein kinase 1) exons 2, 3 and 4. MAPK1 is a signaling molecule of the MAPK signaling pathway that plays important roles in regulating cell proliferation, differentiation, invasion and migration . However, the role of hsa_circ_0004872 in disease is unknown. In this study, we found that the expression of hsa_circ_0004872 was dramatically decreased in GC tissues compared with adjacent nontumor tissues. Consistent with previous reports, hsa_circ_0004872 was more stable than the linear RNA and not easily degraded by RNase R. The biological function assays in vitro showed that hsa_circ_0004872 inhibited GC cells proliferation, invasion and migration. Moreover, in vivo nude mice experiments also verified the tumor-suppressor activity of hsa_circ_0004872. In addition, Zhu et al.  also found that hsa_circ_0004872 was downregulated in glioblastoma, revealing the tumor suppressor role of hsa_circ_0004872 in different tumors.
It is well known that exonic circRNAs are usually located in the cytoplasm. Our FISH experiment also showed that hsa_circ_0004872 mainly existed in the cytoplasm. The most frequently reported function of exonic circRNAs is their role as an “miRNA sponge” to bind with miRNAs and inhibit the function of the miRNAs, thus protecting the target genes from miRNA-mediated degradation . For example, a report from Han et al. showed that exonic circMTO1 suppressed hepatocellular carcinoma progression by acting as a sponge of oncogenic miR-9 and increasing p21 expression . A study from Wang et al. showed that circNT5E acted as a miR-422a sponge to regulate the proliferation, apoptosis, invasion and migration of glioblastoma cells . In this study, we confirmed the “miRNA sponge” role of hsa_circ_0004872 by an RNA binding protein immunoprecipitation experiment (RIP), in which we used an antibody against Ago2, an important component of the miRNA-induced silencing complex , to immunoprecipitate Ago2-miRNA-circRNA complexes and then amplified hsa_circ_0004872 from the complexes. Then, we used different databases to predict the potential miRNAs and selected potential oncogenic miRNAs to investigate their interaction with hsa_circ_0004872 by biotin-coupled probe pull-down assay. Our results showed that hsa_circ_0004872 was able to bind with miR-224. Furthermore, FISH and dual luciferase reporter assays confirmed the colocalization and interaction between hsa_circ_0004872 and miR-224 in GC cells. Further statistical analysis of the expression of hsa_circ_0004872 and miR-224 in GC tissues showed a negative correlation between them.
We next sought to elucidate the biological role and mechanism of miR-224 in GC cells and found that miR-224 exerted critical functions in GC progression and metastasis by targeting p21 and Smad4. p21 is an inhibitor of the cell cycle and affects the progression of the cell cycle, leading to the inhibition of cell proliferation . Smad4 is a member of the TGF-β signaling pathway and plays an important role in cell growth, apoptosis, invasion and migration . To further investigate whether the biological function of hsa_circ_0004872 in GC cells was associated with miR-224, we performed a rescue experiment. Our findings suggested that hsa_circ_0004872 influences the biological functions of GC by regulating p21 and Smad4 via miR-224.
It has been reported that the formation of circRNAs is regulated by the regulatory proteins MBl , QKI [36, 37], adenosine deaminase 1 (ADAR1) [38, 39] and other RNA binding proteins. To determine which factors may affect the level of hsa_circ_0004872 in GC, we analyzed the public database NCBI GEO (GSE27342 and GSE66229) and found only ADAR1 was elevated in GC tissues compared with adjacent nontumor tissues in the both GEO databases. ADAR1, a member of the ADAR family, has two subtypes p110 and p150 [45, 46]. ADAR1 catalyzes the deamination of adenosine and converts adenosine (A) to inosine (I). ADAR1 was reported to strongly inhibit the production of circRNA, and its A-to-I editing process usually occurs near the reverse complementary matches (RCMs) region, a structure believed to be critical for circRNA formation . Shi et al. reported that ADAR1 in hepatocellular carcinoma suppresses the expression of CircARSP91 . Therefore, we investigated the regulatory effect of ADAR1 on hsa_circ_0004872 expression and confirmed that ADAR1 could decrease the level of hsa_circ_0004872 in GC cells. Taken together, our findings demonstrated that the overexpression of ADAR1 contributed to the downregulation of hsa_circ_0004872 in GC and that the low expression of hsa_circ_0004872 promoted the tumorigenesis and progression of GC by regulating p21 and Smad4 via miR-224.
We further investigated whether Smad4, as a transcription factor, was able to regulate ADAR1 expression, thereby forming a regulatory feedback loop of Smad4/ADAR1/hsa_circ_0004872/miR-224/Smad4. We used the Jaspar database to predict the binding sites of Smad4 in the promoter region of ADAR1 and found five potential sites. As expected, our results showed that Smad4 overexpression decreased the expression level of ADAR1 and increased the level of hsa_circ_0004872 in GC cells. ChIP and dual luciferase reporter experiments verified that Smad4 could bind to the promoter region of ADAR1 to regulate its expression in GC cells. Furthermore, according to the public database NCBI GEO (GSE66229), the expression level of Smad4 was significantly low in GC tissues. Statistical analysis showed that there was a remarkable negative correlation between ADAR1 and Smad4 in GC tissues. The Kaplan-Meier Plotter Database indicated that higher expression of ADAR1 or lower expression of Smad4 in GC tissues resulted in a shorter survival period.
In summary, in this study, we investigated the expression, role and regulatory mechanism of the tumor-suppressor circular RNA hsa_circ_0004872 in GC. We found that hsa_circ_0004872 was dramatically downregulated in GC tissues, which was due at least partially to the overexpression of ADAR1. hsa_circ_0004872 inhibited the proliferation, invasion and migration of GC cells by acting as a “miRNA sponge” to bind with miR-224 and increase the expression of the endogenous miR-224 targets p21 and Smad4. Moreover, Smad4, as a transcription factor, could also regulate hsa_circ_0004872 expression by directly binding to the promoter region of ADAR1 and decreasing its expression. Therefore, we revealed a novel regulatory feedback loop formed by hsa_circ_0004872/miR-224/Smad4/ADAR1 in GC. Our findings may provide new ideas and targets for the diagnosis and treatment of GC.
Availability of data and materials
All data generated or analyzed during this study are included either in this article or in the supplementary files.
3′- untranslated region
Cell counting kit- 8
Fluorescence in situ hybridization
Hematoxylin and eosine
Messenger ribonucleic acid
Real-time quantification PCR
RNA-binding protein immunoprecipitation
Reverse transcriptase PCR
Small interfering RNA
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.
Sasako M, Sakuramoto S, Katai H, Kinoshita T, Furukawa H, Yamaguchi T, Nashimoto A, Fujii M, Nakajima T, Ohashi Y. Five-year outcomes of a randomized phase III trial comparing adjuvant chemotherapy with S-1 versus surgery alone in stage II or III gastric cancer. J Clin Oncol. 2011;29(33):4387–93.
Correa P. Is gastric cancer preventable? Gut. 2004;53(9):1217–9.
Li J, Sun D, Pu W, Wang J, Peng Y. Circular RNAs in cancer: biogenesis, function, and clinical significance. Trends Cancer. 2020;6(4):319–36.
Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32:453–61.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M, Loewer A, Ziebold U, Landthaler M, Kocks C, le Noble F, Rajewsky N. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333–8.
Barrett SP, Salzman J. Circular RNAs: analysis, expression and potential functions. Development. 2016;143(11):1838–47.
Qu S, Liu Z, Yang X, Zhou J, Yu H, Zhang R, Li H. The emerging functions and roles of circular RNAs in cancer. Cancer Lett. 2018;414:301–9.
Wang Y, Mo Y, Gong Z, Yang X, Yang M, Zhang S, Xiong F, Xiang B, Zhou M, Liao Q, Zhang W, Li X, Li X, Li Y, Li G, Zeng Z, Xiong W. Circular RNAs in human cancer. Mol Cancer. 2017;16(1):25.
Han D, Li J, Wang H, Su X, Hou J, Gu Y, Qian C, Lin Y, Liu X, Huang M, Li N, Zhou W, Yu Y, Cao X. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017;66(4):1151–64.
Wang L, Long H, Zheng Q, Bo X, Xiao X, Li B. Circular RNA circRHOT1 promotes hepatocellular carcinoma progression by initiation of NR2F6 expression. Mol Cancer. 2019;18(1):119.
Cheng Z, Yu C, Cui S, Wang H, Jin H, Wang C, Li B, Qin M, Yang C, He J, Zuo Q, Wang S, Liu J, Ye W, Lv Y, Zhao F, Yao M, Jiang L. circTP63 functions as a ceRNA to promote lung squamous cell carcinoma progression by upregulating FOXM1. Nat Commun. 2019;10(1):3200.
Wang S, Zhang X, Li Z, Wang W, Li B, Huang X, Sun G, Xu J, Li Q, Xu Z, Xia Y, Wang L, Zhang Q, Li Q, Zhang L, Chen J, Wu Y, Cao J, Xu P, Zhang D, Xu H, Xu Z. Circular RNA profile identifies circOSBPL10 as an oncogenic factor and prognostic marker in gastric cancer. Oncogene. 2019;38(44):6985–7001.
Xu JZ, Shao CC, Wang XJ, Zhao X, Chen JQ, Ouyang YX, Feng J, Zhang F, Huang WH, Ying Q, Chen CF, Wei XL, Dong HY, Zhang GJ, Chen M. circTADA2As suppress breast cancer progression and metastasis via targeting miR-203a-3p/SOCS3 axis. Cell Death Dis. 2019;10(3):175.
Lu Q, Liu T, Feng H, Yang R, Zhao X, Chen W, Jiang B, Qin H, Guo X, Liu M, Li L, Guo H. Circular RNA circSLC8A1 acts as a sponge of miR-130b/miR-494 in suppressing bladder cancer progression via regulating PTEN. Mol Cancer. 2019;18(1):111.
Wang R, Zhang S, Chen X, Li N, Li J, Jia R, Pan Y, Liang H. CircNT5E acts as a sponge of miR-422a to promote glioblastoma tumorigenesis. Cancer Res. 2018;78(17):4812–25.
Vo JN, Cieslik M, Zhang Y, Shukla S, Xiao L, Zhang Y, Wu YM, Dhanasekaran SM, Engelke CG, Cao X, Robinson DR, Nesvizhskii AI, Chinnaiyan AM. The landscape of circular RNA in cancer. Cell. 2019;176(4):869–881. e13.
Zhang XO, Wang HB, Zhang Y, Lu X, Chen LL, Yang L. Complementary sequence-mediated exon circularization. Cell. 2014;159(1):134–47.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.
Petkovic S, Muller S. RNA circularization strategies in vivo and in vitro. Nucleic Acids Res. 2015;43(4):2454–65.
Du WW, Yang W, Liu E, Yang Z, Dhaliwal P, Yang BB. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016;44(6):2846–58.
Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, Evantal N, Memczak S, Rajewsky N, Kadener S. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. 2014;56(1):55–66.
Li Z, Huang C, Bao C, Chen L, Lin M, Wang X, Zhong G, Yu B, Hu W, Dai L, Zhu P, Chang Z, Wu Q, Zhao Y, Jia Y, Xu P, Liu H, Shan G. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22(3):256–64.
Zhang M, Zhao K, Xu X, Yang Y, Yan S, Wei P, Liu H, Xu J, Xiao F, Zhou H, Yang X, Huang N, Liu J, He K, Xie K, Zhang G, Huang S, Zhang N. A peptide encoded by circular form of LINC-PINT suppresses oncogenic transcriptional elongation in glioblastoma. Nat Commun. 2018;9(1):4475.
Pamudurti NR, Bartok O, Jens M, Ashwal-Fluss R, Stottmeister C, Ruhe L, Hanan M, Wyler E, Perez-Hernandez D, Ramberger E, Shenzis S, Samson M, Dittmar G, Landthaler M, Chekulaeva M, Rajewsky N, Kadener S. Translation of CircRNAs. Mol Cell. 2017;66(1):9–21.e7.
Legnini I, Di Timoteo G, Rossi F, Morlando M, Briganti F, Sthandier O, Fatica A, Santini T, Andronache A, Wade M, Laneve P, Rajewsky N, Bozzoni I. Circ-ZNF609 is a circular RNA that can be translated and functions in myogenesis. Mol Cell. 2017;66(1):22–37.
Wang S, Zhang Y, Cai Q, Ma M, Jin L, Weng M, Zhou D, Tang Z, Wang J, Quan Z. Circular RNA FOXP1 promotes tumor progression and Warburg effect in gallbladder cancer by regulating PKLR expression. Mol Cancer. 2019;18(1):145.
Xu K, Dong C, Wang Z, Ma J, Zhou J, Chen N, Lv L, Zheng Y, Hu X, Zhang Y, Li J. Annotation and functional clustering of circRNA expression in rhesus macaque brain during aging. Cell Discov. 2018;4:48.
Xu H, Wang C, Song H, Xu Y, Guang J. RNA-Seq profiling of circular RNAs in human colorectal cancer liver metastasis and the potential biomarkers. Mol Cancer. 2019;18:8.
Cui R, Meng W, Sun HL, Kim T, Ye Z, Fassan M, Jeon YJ, Li B, Vicentini C, Peng Y, Lee TJ, Luo Z, Liu L, Xu D, Tili E, Jin V, Middleton J, Chakravarti A, Lautenschlaeger T, Croce CM. MicroRNA-224 promotes tumor progression in nonsmall cell lung cancer. Proc Natl Acad Sci U S A. 2015;112(31):E4288–97.
Zhou J, Hu M, Wang F, Song M, Huang Q, Ge B. miR-224 controls human colorectal cancer cell line HCT116 proliferation by targeting Smad4. Int J Med Sci. 2017;14(10):937–42.
Wang H, Zhu LJ, Yang YC, Wang ZX, Wang R. MiR-224 promotes the chemoresistance of human lung adenocarcinoma cells to cisplatin via regulating G1/S transition and apoptosis by targeting p21(WAF1/CIP1). Br J Cancer. 2014;111(2):339–54.
Li Q, Ding C, Chen C, Zhang Z, Xiao H, Xie F, Lei L, Chen Y, Mao B, Jiang M, Li J, Wang D, Wang G. miR-224 promotion of cell migration and invasion by targeting Homeobox D 10 gene in human hepatocellular carcinoma. J Gastroenterol Hepatol. 2014;29(4):835–42.
Yang C, Yuan W, Yang X, Li P, Wang J, Han J, Tao J, Li P, Yang H, Lv Q, Zhang W. Circular RNA circ-ITCH inhibits bladder cancer progression by sponging miR-17/miR-224 and regulating p21, PTEN expression. Mol Cancer. 2018;17(1):19.
Chen LL. The biogenesis and emerging roles of circular RNAs. Nat Rev Mol Cell Biol. 2016;17(4):205–11.
Conn SJ, Pillman KA, Toubia J, Conn VM, Salmanidis M, Phillips CA, Roslan S, Schreiber AW, Gregory PA, Goodall GJ. The RNA binding protein quaking regulates formation of circRNAs. Cell. 2015;160(6):1125–34.
Han J, Meng J, Chen S, Wang X, Yin S, Zhang Q, Liu H, Qin R, Li Z, Zhong W, Zhang C, Zhang H, Tang Y, Lin T, Gao W, Zhang X, Yang L, Liu Y, Zhou HG, Sun T, Yang C. YY1 complex promotes quaking expression via super-enhancer binding during EMT of hepatocellular carcinoma. Cancer Res. 2019;79(7):1451–64.
Ivanov A, Memczak S, Wyler E, Torti F, Porath HT, Orejuela MR, Piechotta M, Levanon EY, Landthaler M, Dieterich C, Rajewsky N. Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals. Cell Rep. 2015;10(2):170–7.
Hong H, An O, Chan THM, Ng VHE, Kwok HS, Lin JS, Qi L, Han J, Tay DJT, Tang SJ, Yang H, Song Y, Bellido Molias F, Tenen DG, Chen L. Bidirectional regulation of adenosine-to-inosine (A-to-I) RNA editing by DEAH box helicase 9 (DHX9) in cancer. Nucleic Acids Res. 2018;46(15):7953–69.
Dou N, Yu S, Ye X, Yang D, Li Y, Gao Y. Aberrant overexpression of ADAR1 promotes gastric cancer progression by activating mTOR/p70S6K signaling. Oncotarget. 2016;7(52):86161–73.
Shi L, Yan P, Liang Y, Sun Y, Shen J, Zhou S, Lin H, Liang X, Cai X. Circular RNA expression is suppressed by androgen receptor (AR)-regulated adenosine deaminase that acts on RNA (ADAR1) in human hepatocellular carcinoma. Cell Death Dis. 2017;8(11):e3171. https://doi.org/10.1038/cddis.2017.556.
Kamiyama M, Naguro I, Ichijo H. In vivo gene manipulation reveals the impact of stress-responsive MAPK pathways on tumor progression. Cancer Sci. 2015;106(7):785–96.
Zhu J, Ye J, Zhang L, Xia L, Hu H, Jiang H, Wan Z, Sheng F, Ma Y, Li W, Qian J, Luo C. Differential expression of circular RNAs in glioblastoma multiforme and its correlation with prognosis. Transl Oncol. 2017;10(2):271–9.
Winter J, Jung S, Keller S, Gregory RI, Diederichs S. Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol. 2009;11:228–34.
George CX, Samuel CE. Human RNA-specific adenosine deaminase ADAR1 transcripts possess alternative exon 1 structures that initiate from different promoters, one constitutively active and the other interferon inducible. Proc Natl Acad Sci U S A. 1999;96:4621–6.
George CX, Wagner MV, Samuel CE. Expression of interferon-inducible RNA adenosine deaminase ADAR1 during pathogen infection and mouse embryo development involves tissue-selective promoter utilization and alternative splicing. J Biol Chem. 2005;280:15020–8.
We thank assistant researcher Lidong Zu in the Center of Pathology, School of Medicine, Shanghai Jiao Tong University, China for providing the plasmids p-3 × flag-CMV-Smad4 and p-3 × flag-CMV-14.
This study was supported by the National Natural Science Foundation of China [81871621,81671981, 81772143].The Department of Science and Technology of Shandong Province (Nos. 2018CXGC1208).
Ethics approval and consent to participate
The present study was approved the ethics committee of School of Basic Medical Sciences, Shandong University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
qRT-PCR analysis of the expression of hsa_circ_0004872, hsa_circ_0002483,hsa_circ_0000847, hsa_circ_0001566 in 42 paired GC tissues and corresponding nontumor tissues.
qRT-PCR analysis of the expression of hsa_circ_0004872 in BGC-823 and SGC-7901 cells transfected with hsa_circ_0004872 overexpression vector (pLCDH-circ_4872) (A) or the miR-224 binding site mutated hsa_circ_0004872 overexpression vector ((pLCDH-circ_4872-Mut) (B).
hsa_circ_0004872 siRNAs promote the proliferation, invasion and migration of GC cells. (A) Schematic representation of the siRNA sequences specifically target the junction site of hsa_circ_0004872. (B) qRT-PCR analysis of hsa_circ_0004872 and MAPK1 mRNA level in BGC-823 cells transfected with hsa_circ_0004872 siRNAs (si-circ_4872) or the control siRNA. (C) EdU analysis of the cell proliferation ability in BGC-823 cells transfected with the si-circ_4872 or the control siRNA.Representative images are shown. Scale bar: 20 μm. (D) Statistical analysis of the EdU-positive cell ratio in the cells transfected with si-circ_4872 or the control siRNA. (E) The scratch wound healing assays in BGC-823 cells transfected with the si-circ_4872 or the control siRNA. Scale bar: 500 μm. (F) Statistical analysis of the cell migration in the scratch wound healing assays.The data are expressed as the means±SD from three experiments. (G) Transwell invasion and migration assay in BGC-823 cells transfected with the si-hsa_circ_0004872 or the control siRNA. Scale bar: 100 μm. (H) Statistical analysis of the cell numbers passing through the transwell chamber in the transfected BGC-823 cells. The data are expressed as the means±SD from three experiments.
qRT-PCR analysis of the expression of hsa_circ_0004872 in different GC cells.
Schematic diagam of dual luciferase vector. (A) Schematic diagam of dual luciferase vector pMIRGLO-circ_4872-WT/Mut. Upper: diagram of the luciferase reporter construct containing the sequences of hsa_circ_0004872. The mutations were generated at the predicted miR-224 binding sites in the hsa_circ_0004872 sequences. Lower: the predicted complementary sequences of miR-224 in the sequences of hsa_circ_0004872. (B) Schematic diagam of dual luciferase vector pMIR-Smad4(p21)-WT/Mut. Upper: diagram of the luciferase reporter construct containing 3’UTR sequences of Smad4 (p21). The mutations were generated at the predicted miR-224 binding sites located in the 3’UTR of Smad4(p21). Lower: the predicted complementary sequences of miR-224 in the 3’UTR of Smad4 (p21). (C) Schematic diagram of dual luciferase vector pGL3-ADAR1-WT/Mut. Upper: diagram of the luciferase reporter construct containing promoter sequence of ADAR1. The mutations were generated at the predicted Smad4 binding sites located in promoter sequence of ADAR1. Lower: the predicted complementary sequences of Smad4 in promoter sequence of ADAR1.
qRT-PCR analysis of the expression of miR-224 (A) and ADAR1 (B) in 39 paired GC tissues and corresponding nontumor tissues.
miR-224 inhibitor inhibited the proliferation, invasion and migration in GC cells (A) The expression level of miR-224 was analyzed with qRT-PCR in BGC-823 and SGC-7901 cells transfected with miR-224 inhibitor or the control inhibitor. (B) EdU analysis of the cell proliferation ability in BGC-823 and SGC-7901 cells transfected with miR-224 inhibitor or the control inhibitor. Scale bar: 20 μm. (C) Statistical analysis of the EdU-positive cell ratio in the transfected cells. (D) CCK-8 analysis of the cell proliferation ability in BGC-823 and SGC-7901 cells transfected with miR-224 inhibitor or the control inhibitor. (E) The scratch wound healing assays of the migration ability in transfected BGC-823 and SGC-7901 cells. Scale bar: 500 μm. (F) Statistical analysis of the scratch wound healing assays. (G) Transwell assay of the migration (without matrigel) and invasion ability (with matrigel) in BGC-823 and SGC-7901 cells transfected with miR-224 inhibitor or the control inhibitor. Scale bar: 100 μm. (H) Statistical analysis of the cell numbers passing through the transwell chamber in the transfected BGC-823 and SGC-7901 cells. All datas were the means ± SD.
The expression of ADAR1, MBl and QKI were analyzed in NCBI GEO database GSE27342 and GSE66229. (A) The expression level of ADAR1 was analyzed with paired t-tests (n = 80, p < 0.0001) in GEO database GSE27342. (B) The expression level of ADAR1 was analyzed with unpaired t-tests (n = 400, p < 0.0001) in GSE66229. (C) The expression level of MBL was analyzed with paired t-tests (n = 80, p = 0.7912) in GSE27342. (D) The expression level of MBL was analyzed with unpaired t-tests (n = 400, p = 0.7566) in GSE66229. (E) The expression level of QKI was analyzed with paired t-tests (n = 80, p = 0.3101) in GSE27342. (F) The expression level of QKI was analyzed with unpaired t-tests (n = 400, p < 0.0001) in GSE66229.
Patients, tumor characteristics and hsa_circ_0004872 expression in GC samples.
siRNA sequences in the research.
Primer sequences in this study.
Probe sequences of hsa_circ_0004872 and miR-224 in this study.
CircRNA-seq analysis of the differentially expressed circRNAs in GC tissue and corresponding nontumor tissue.
. Predicted miRNAs with potential binding ability with hsa_circ_0004872 in different databases.
About this article
Cite this article
Ma, C., Wang, X., Yang, F. et al. Circular RNA hsa_circ_0004872 inhibits gastric cancer progression via the miR-224/Smad4/ADAR1 successive regulatory circuit. Mol Cancer 19, 157 (2020). https://doi.org/10.1186/s12943-020-01268-5
- Gastric cancer (GC)