- Open Access
circLMTK2 acts as a sponge of miR-150-5p and promotes proliferation and metastasis in gastric cancer
Molecular Cancer volume 18, Article number: 162 (2019)
As a novel class of non-coding RNAs, circular RNAs (circRNAs) are key regulators of the development and progression of different cancers. However, little is known about the function and biological mechanism of circLMTK2, also named hsa_circ_0001725, in gastric cancer (GC) tumourigenesis.
circLMTK2 was identified in ten paired cancer specimens and adjacent normal tissues by RNA sequencing and genome-wide bioinformatic analysis and verified by quantitative real-time PCR (qRT-PCR). Knockdown or exogenous expression of circLMTK2 combined with in vitro and in vivo assays were performed to prove the functional significance of circLMTK2. The molecular mechanism of circLMTK2 was demonstrated by searching the CircNet database and confirmed by RNA in vivo precipitation assays, western blotting, luciferase assays and rescue experiments.
circLMTK2 was frequently upregulated in GC tissues, and high circLMTK2 expression was associated with poor prognosis, lymph node metastasis and poor TNM stage in GC patients. Functionally, circLMTK2 overexpression promoted GC cell proliferation and tumourigenicity in vitro and in vivo. Furthermore, ectopic circLMTK2 expression enhanced GC cell migration and invasion in vitro and tumour metastasis in vivo. In addition, we demonstrated that circLMTK2 could sponge miR-150-5p, thus indirectly regulating the c-Myc expression and contributing to GC tumourigenesis.
Our findings demonstrate that circLMTK2 functions as a tumour promoter in GC through the miR-150-5p/c-Myc axis and could thus be a prognostic predictor and therapeutic target for GC.
Gastric cancer (GC) remains one of the most common and lethal malignancies worldwide, with a particularly high disease incidence rate in East Asia . Although the incidences and mortality trends for GC have declined in recent years , the outcomes of this disease are still among the poorest of all solid-organ tumours, predominantly due to the frequent presence of advanced stage disease with lymphatic or distant metastasis . Because there are limited therapeutic approaches for treating advanced GC, it is urgent that we search for novel biomarkers and prognostic indicators to reflect the disease status and develop more therapeutic targets for this lethal disease.
A poorly characterized component of the GC transcriptome is circular transcripts (circRNAs), which have been implicated in other diseases [4,5,6]. circRNAs are a type of RNA formed by back-splicing , and compared to their linear counterparts, they are highly stable due to their covalently closed loop . Advances in high-throughput sequencing technology and novel bioinformatics algorithms have facilitated the systematic detection of circRNAs [9, 10]. circRNAs have been posited to function as sponges of microRNAs (miRNAs) and decoys of RNA-binding proteins (RBPs), which affect biological processes such as proliferation, metastasis, and apoptosis in tumour cells . Therefore, circRNAs may be potential biomarkers or therapeutic targets .
In this study, we generated de-ribosome RNA sequencing data from GC tissues from ten patients, and identified approximately 35,350 circRNA candidates (at least two unique back-splicing reads). We characterized one circRNA derived from the LMTK2 gene locus, termed circLMTK2 (hsa_circ_0001725), that was upregulated in GC patients. In vitro and in vivo experiments showed that circLMTK2 can promote cell growth and metastasis. Moreover, a clinical analysis showed a significant negative correlation between circLMTK2 and the patient’s prognosis. Our study provides a new insight into the pathogenesis of GC.
Materials and methods
We retrospectively collected ten paired cancer specimens and adjacent normal tissues from patients with gastric cancer who had surgically proven primary GC and received a D2 radical gastrectomy (R0 resection) at the Department of General Surgery of Northern Jiangsu People’s Hospital between November 2008 and December 2011. None of these patients received preoperative chemotherapy or radiotherapy. Clinicopathological features, which included age, sex, tumour site, tumour size, differentiation grade, TNM stage (American Joint Committee on Cancer classification, AJCC), lymphatic invasion and neural invasion, are shown in Table 1. The median follow-up time was 25.0 months (range: 1–85 months). The follow-up interval began on the date of surgery and ended on the date of disease progression, death or the last clinical investigation. This study was approved by the Medical Ethics Committee of Northern Jiangsu People’s Hospital. Written informed consent was obtained from all participants.
Total RNA was isolated using TRIzol reagent (Life Technologies, Carlsbad, CA, USA). Approximately 3 μg of total RNA from each sample was subjected to the RiboMinus Eukaryote Kit (Qiagen, Valencia, CA) to remove ribosomal RNA prior to RNA-seq library construction. Strand-specific RNA-seq libraries were prepared using a NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB, Beverly, MA, USA). Briefly, approximately 50 ng of ribosome-depleted RNA samples was fragmented and then used for first- and second-strand cDNA synthesis with random hexamer primers. A dUTP mix was used for second-strand cDNA synthesis. An End-It DNA End Repair Kit was used to repair the ends of the double-stranded cDNA fragments, which were then modified by the Klenow fragment so that an A was added to the 3′ end of the DNA fragments and were finally ligated to adapters. The ligated products were purified and treated with uracil DNA glycosylase (UDG) to remove the second-strand cDNA. Purified first-strand cDNA was subjected to 13–15 cycles of PCR amplification, followed by library analysis with a Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA); the cDNA was then sequenced using a HiSeq 2000 system (Illumina, San Diego, CA, USA) and a 100-bp paired-end run.
Identification and quantification of circRNAs
The RNA-seq fastq reads were first mapped to the human reference genome (GRCh38/hg38) obtained from the UCSC genome database (http://genome.ucsc.edu/) using TopHat2. The unmapped reads were then used to identify circRNAs as previously described . Briefly, the unmapped reads were processed to 20-nucleotide anchors from both ends of the sequencing read. Anchors that aligned in the reverse orientation (head-to-tail) represent a back-spliced junction. Anchor alignment was extended such that the complete read alignment and the breakpoint were flanked by a GT/AG splice site. The total number of reads that spanned back-spliced junctions was used as an absolute measure of circRNA abundance. The genomic regions that mapped to inferred circRNAs were annotated according to RefSeq. Gene coordinates were downloaded from the RefGene tables in the UCSC Genome Browser. The host genes of circRNAs were determined using a custom script. For each circRNA, we searched for the longest transcript fragment whose boundaries (5′ end or 3′ end) exactly matched both ends of the circRNA in the same strand and then defined the corresponding gene of the transcript fragment as the host gene of the circRNA.
Cell culture and treatments
BGC-823 and AGS cells were cultured in RPMI 1640 supplemented with 10% foetal bovine serum (FBS) and 1% penicillin-streptomycin at 37 °C and 5% CO2. HEK 293 T and MGC-803 cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% FBS and antibiotics. Transcription was blocked by the addition of 2 μg/ml actinomycin D or DMSO (Sigma-Aldrich, St. Louis, MO, USA), which served as a control for the cell culture medium.
RNA fluorescence in situ hybridization (FISH)
In situ hybridization was performed using specific probes to the circLMTK2 sequence. PCR fragments with the T7 promoter were amplified with specific primers for the back-splice region of circLMTK2. Primers are listed in Additional file 1: Table S3. Digoxin or biotin-labelled RNA probes were transcribed from PCR fragments using the DIG or biotin RNA labelling mix and T7 RNA polymerase (Roche) according to the manufacturers’ instructions. AGS cells were grown to the exponential phase and were 80–95% confluent at the time of fixation. After pre-hybridization (1 × PBS/0.5% Triton X-100), the cells were hybridized in hybridization buffer (40% formamide, 10% dextran sulfate, 1 × Denhardt’s solution, 4× SSC, 10 mM DDT, 1 mg ml− 1 yeast transfer RNA, 1 mg ml−1sheared salmon sperm DNA) with DIG-labelled probes specific to circLMTK2 at 60 °C overnight. Signals were detected using a tyramide-conjugated Alexa 488 fluorochrome TSA kit (Life Technologies). Nuclei were counterstained with 4,6-diamidino-2-phenylindole. Images were acquired on a Leica SP5 confocal microscope (Leica Micosystems, Mannheim, Germany).
circRNA in vivo precipitation (circRIP)
A biotin-labelled circLMTK2 probe (5′-CTACCTGTTTGACCAGGGTCTCTGGGTGT-3′-biotin) was designed and synthesized by RiboBio (Guangzhou, China), and a circRIP assay was performed as described. circLMTK2-overexpressing BGC-823 cells were seeded in a 10-cm dish. After reaching sufficient confluency, the cells were transfected with the specific biotin-tagged probe or control probe at a final concentration of 200 nmol/L. Then, the cells were fixed with 1% formaldehyde for 10 min, lysed and sonicated. After centrifugation, 50 μl of the supernatant was retained as input, and the remaining cell lysis solution was incubated with a circLMTK2-specific probe-streptavidin-dynabead (M-280, Invitrogen) mixture overnight at 30 °C. The next day, the M-280 dynabead-probe-circRNA mixture was washed and incubated with 200 μl lysis buffer and proteinase K to reverse the formaldehyde crosslinking. Finally, total RNA was extracted from the mixture using an miRNeasy Mini Kit according to the manufacturer’s instructions (Qiagen).
Nucleic acid preparation and quantitative real-time polymerase chain reaction (qRT-PCR)
Genomic DNA was isolated with a QIAamp DNA Mini Kit (Qiagen, Valencia, CA, USA), and total RNA was isolated using TRIzol reagent (Life Technologies, Carlsbad, CA, USA).The nuclear and cytoplasmic fractions were isolated using NE-PER Nuclear and Cytoplasmic Extraction Reagents (Thermo Scientific). Total RNA from the nuclear and cytoplasmic fractions was isolated with TRIzol. Complementary DNA was synthesized using a PrimeScript RT reagent kit (Takara Bio Inc., Dalian, China), and RT-PCR was performed using SYBR Premix Ex Taq (Takara Bio Inc.). For miRNA measurements, mature miR-150-5p was reverse-transcribed and quantified with TaqMan® RT primers and probes, and the data were normalized to U6 small nuclear RNA expression using predesigned TaqMan assays (Applied Biosystems, Foster City, USA). The primers are listed in Additional file 1: Table S3.
The circLMTK2 genomic region and its wild-type flanking introns were amplified from the cDNA using PrimerSTAR Max DNA Polymerase Mix (Takara) and were subcloned into a pcDNA3.0 vector. In the luciferase reporter assay, circLMTK2 was amplified from the cDNA and was inserted into the region directly downstream of a cytomegalovirus (CMV) promoter-driven firefly luciferase cassette in the pCDNA3.0 vector. The 3′-UTR sequence of c-Myc was amplified and inserted into the psiCHECK-2 vector (Promega, Madison, USA). Mutations in the miRNA binding sites in the c-Myc 3′-UTR sequence were generated using a Mut Express II Fast Mutagenesis Kit (Vazyme, NanJing, China). The constructs were verified by sequencing. The primers are listed in Additional file 1: Table S3.
SiRNA and miRNA mimics and inhibitors were synthesized by RiboBio (Guangzhou, China). The sequences that were used are shown in Additional file 1: Table S3. The cells were transfected using Lipofectamine RNAiMax (Life Technologies).
Luciferase reporter assay
HEK 293 T cells (5 × 103) were seeded into 96-well plates and were cotransfected with a mixture of 50 ng of psiCHECK-2 vector, and miRNA mimics. After 48 h of incubation, the firefly and Renilla luciferase activities were quantified with a dual-luciferase reporter assay (Promega, Madison, WI, USA).
A Matrigel Cell Migration Assay and Invasion System was used to measure cell migration and invasion in vitro as described previously. For cell migration assays, AGS, MGC-803 and BGC-823 cells (5 × 104 cells/well) were suspended in 50 μL serum-free RPMI 1640 and DMEM and placed in the upper collagen-coated chambers (8-mm pore size; Millipore, Temecula, CA, USA) of each transwell insert (Matrigel: serum-free medium 1:5) respectively. Next, 800 μL of media with 10% FBS was placed in the lower chamber and incubated for 24 and 48 h. After incubation, cells adhering to the upper surface of the membrane were removed with a cotton swab. The migrated cells that adhered to the lower surface were fixed with 4% formaldehyde in PBS, followed by staining with a 1% crystal violet solution. All samples were then examined and photographed under light microscopy at × 200. The total cells were counted, and the total percentage of inhibition based on the cells in each image was measured as described previously. For the invasion assays, which were performed almost the same as the cell migration assays, Matrigel was used instead of collagen on the filter membrane as described previously.
5-Ethynyl-2′-deoxyuridine (EdU) incorporation assay
EdU assays were performed with a Cell-Light EdU DNA Cell Proliferation Kit (RiboBio, Shanghai, China). Cells (1 × 104) were seeded in each well of 96-well plates. After incubation with 50 μM EdU for 2 h, the cells were fixed in 4% paraformaldehyde and stained with Apollo Dye Solution. Hoechst 33342 was used to stain the nucleic acids within the cells. Images were obtained with an Olympus FSX100 microscope (Olympus, Tokyo, Japan), and the number of EdU-positive cells was counted.
Cell proliferation was assessed by Cell Counting Kit-8 assays (Dojindo Laboratories, Kumamoto, Japan). Cells (1 × 103) were seeded into 96-well plates. Then, 10 μl of CCK-8 solution was added to each well on days 1, 2, 3, 4 and 5. After 2 h of incubation at 37 °C, the absorbance at 450 nM was measured using an automatic microplate reader (Synergy4; BioTek, Winooski, VT, USA). The experiment was repeated three times.
Colony formation assay
For the colony formation assays, cells were trypsinized, and 1 × 103 cells were plated in 6-well plates and incubated at 37 °C for 14 days. Colonies were dyed with a 0.1% crystal violet and 20% methanol solution. Cell colonies were then counted and analysed.
BGC-823 cells that stably expressed circLMTK2 and control cells were harvested and suspended in DMEM without FBS. Sixteen mice (male BALB/c-nu/nu, 6 weeks old) were divided randomly into two groups, and each mouse was injected subcutaneously in the lower back with 2 × 106 cells in 200 μl of DMEM without FBS. The mice were monitored weekly for tumour weight and tumour volume (volume = length ×width2/2). At approximately 4 weeks after injection, the mice were sacrificed, and the tumours were dissected and weighed. For the in vivo tumour metastasis studies, eighteen mice (male BALB/c-nu/nu, 6 weeks old) were divided randomly into two groups. BGC-823 cells stably transfected with circLMTK2 lentiviruses or control vector were injected into the lateral tail veins of the nude mice (2 × 106 cells per mouse). Fifty days later, the mice were sacrificed and examined for the numbers of lung metastatic colonies. Paraffin sections were imaged with a Leica Microsystems Microscope (Leica Biosystems, Wetzlar, Germany). The mouse experiments were conducted according to the Guide for the Care and Use of Laboratory Animals of Yangzhou University. The protocol was approved by the Committee on the Ethics and Welfare of Laboratory Animal Science of Yangzhou University.
Statistical analyses were performed using SPSS 20.0 (IBM, SPSS, Chicago, IL, USA) and GraphPad Prism for Windows, version 6.00 (GraphPad Software, La Jolla, USA). Unless otherwise stated, Student’s t-test and one-way ANOVA were used to determine the statistical significance for comparisons of 2 or more groups. The Pearson correlation coefficient was used to analyse the correlations. Overall survival (OS) was measured from the date of surgery. For OS, patients known to be alive at the time of the last follow-up were censored on the last date of contact. OS curves were calculated with the Kaplan-Meier method and were analysed with the log-rank test. Univariate analysis and multivariate models were constructed using a Cox proportional hazards regression model. P values < 0.05 were considered statistically significant.
Identification of circular RNAs by RNA-seq in GC
We first characterized circular RNA transcripts using an RNA-seq analysis of de-ribosome RNA-seq from ten paired GC/adjacent tissues (Fig. 1a). Each sample was sequenced on an Illumina HiSeq and yielded ~ 60 million reads, which were mapped to the human reference genome (GRCh38/hg38) by TopHat2 . A computational pipeline based on the anchor alignment of unmapped reads was used to identify circRNAs without reliance on gene annotations . Collectively, 35,350 distinct circRNA candidates, which contained at least two unique back-splicing reads, were found in these tissues (Fig. 1b, Additional file 2: Table S1). Among these circRNAs, there are 3450 intergenic circRNAs, and 31,900 overlapped with known genes (Fig. 1c). The expression analysis of these circRNA transcripts revealed that a series of circRNAs was differentially expressed in cancerous tissues compared with those in matched normal tissues. Among the 142 differentially expressed circRNAs (filtered by |FC(fold change)| ≥ 2 and P<0.05), 105 were upregulated, and 37 were downregulated in GC compared with those in normal tissues (Fig. 1d, Additional file 3: Table S2). These circRNAs and their host genes were located in diverse genomic regions (Fig. 1e). Then, we further confirmed the RNA-seq results of four circRNAs in 25 paired normal and cancerous gastric tissues by qRT-PCR analysis (Additional file 1: Figure S1).
Characterization of circLMTK2 in GC
We noted that one of the most differentially expressed circRNAs (chr7:98190727–98,194,572, GRCh38/hg38, Fig. 2a) was derived from a protein-coding gene locus, LMTK2, located on chromosome 7q21.3. Thus, we termed this circRNA as “circLMTK2” (hsa_circ_0001725). The genomic structure shows that circLMTK2 is looped by the tenth and eleventh exons of the LMTK2 gene (NM_014916). We designed outward-facing primers (Additional file 1: Table S3) and confirmed the junction site of circLMTK2 by Sanger sequencing (Fig. 2a). Resistance to digestion by RNase R exonuclease further confirmed that this RNA species exists in a circular form (Fig. 2b). Total RNA was harvested at the indicated time points after treatment with actinomycin D, a transcription inhibitor. Analysis of circLMTK2 and LMTK2 mRNA revealed that the circRNA isoform circLMTK2 is highly stable, as its transcript half-life exceeded 24 h, whereas the linear transcript of LMTK2 mRNA exhibited a half-life of < 4 h in AGS cells (Fig. 2c). Next, qRT-PCR analysis of nuclear and cytoplasmic RNA and FISH against circLMTK2 demonstrated that the circular form of LMTK2 preferentially localized within the cytoplasm in AGS cells (Fig. 2d and e). These results suggest that circLMTK2 is a stable and cytoplasmic circRNA derived from the LMTK2 mRNA.
circLMTK2 promotes GC cell proliferation and tumourigenicity in vitro and in vivo
To investigate the potential role of circLMTK2 in GC cells, we initially explored the effect of blocking circLMTK2 on cell growth. We designed two siRNAs to target the back-splice sequence. As expected, siRNA directed against the back-splice sequence inhibited only the circular transcript of circLMTK2 and did not affect the expression of the LMTK2 linear species (Additional file 1: Figure S3C). Silencing circLMTK2 expression significantly suppressed cell proliferation rates (Fig. 3a and b) and colony formation abilities in both GC cell lines (Fig. 3g and h) and impaired nucleotide synthesis (Fig. 3d and e). In contrast, the ectopic expression of circLMTK2 in BGC-823 cells (designated circLMTK2-OE), which was induced by lentiviruses with circular frames and circLMTK2 sequences (Additional file 1: Figure S3A and B), dramatically promoted cell growth (Fig. 3c and i) and increased the EdU incorporation rate (Fig. 3f). Furthermore, circLMTK2-OE and empty vector cells were inoculated subcutaneously into the flanks of nude mice, and these mice were monitored closely for tumour growth for 4 weeks. Our results illustrated that tumours derived from circLMTK2-OE cells were significantly larger than those derived from empty vector cells, both in terms of tumour volumes and weights (Fig. 3j–l). These results suggest that circLMTK2 significantly promotes GC cell proliferation and tumourigenicity in vitro and in vivo.
circLMTK2 enhances GC cell migration and invasion in vitro and tumour metastasis in vivo
Next, we studied the effects of circLMTK2 expression on GC cell migration and invasion. Our migration and invasion assay results showed that circLMTK2 knockdown significantly eliminated the migratory and invasive capacities of AGS and MGC-803 cells when compared with the indicated controls conditions (Fig. 4a-d). In contrast, the migration and invasion were greater in BGC-823 cells overexpressing circLMTK2 than in empty vector cells (Fig. 4e). These results indicate that circLMTK2 increases cell invasion and metastasis, which we further validated in vivo. The in vivo metastasis assay was performed by injecting circLMTK2-OE and empty vector cells into nude mice through the lateral tail vein to examine the cells’ lung metastasis ability. Metastatic nodules in the lungs were confirmed histologically. The number of metastatic nodules was significantly higher in mice injected with circLMTK2-OE BGC-823 cells than in mice injected with empty vector BGC-823 cells (Fig. 4f). Consistently, these data demonstrate that circLMTK2 significantly promotes GC cell migration and invasion in vitro and tumour metastasis in vivo.
circLMTK2 abundantly sponges miR-150-5p in GC cells
Given that circRNA has been shown to act as a miRNA sponge and that circLMTK2 is abundant and stable in the cytoplasm, we next investigated the ability of circLMTK2 to bind to miRNAs. CircNet (http://syslab5.nchu.edu.tw/CircNet/) and TargetScan (http://www.targetscan.org) were used to predict the potential target miRNAs that could bind with the circLMTK2 sequence (Fig. 5c), and miR-150-5p was selected as the best potential target of circLMTK2 (Fig. 5a). Next, we designed a specific biotin-labelled circLMTK2 probe to perform RNA in vivo precipitation (RIP) to confirm whether miR-150-5p can interact with circLMTK2 in GC cells, which has been reported in several studies [14,15,16]. We purified circLMTK2-associated RNA and performed qPCR to measure circLMTK2 and miR-150-5p expression. The results showed a significant enrichment of circLMTK2 and miR-150-5p compared to the controls (Fig. 5b), indicating that miR-150-5p was sponged by circLMTK2 in GC cells. We further constructed a circLMTK2 fragment and inserted it immediately downstream of the luciferase reporter gene (LUC+ circLMTK2). Then, miR-150-5p or some potential miRNAs mimics were transfected with reporter gene into HEK 293 T cells. A significant reduction in luciferase reporter activity was observed relative to co-transfection with control mimic or other miRNAs (Fig. 5d). In addition, there was a significant inverse correlation between miR-150-5p and circLMTK2 expression levels in 120 GC samples (Fig. 5e). Together, these results suggest that circLMTK2 may serve as a binding platform for miR-150-5p.
circLMTK2 promotes GC cell growth and metastasis by sponging miR-150-5p and downregulating c-myc
Previous studies have described the tumour suppressor role of miR-150-5p in various cancers [17,18,19]; however, the potential mechanisms of miR-150-5p regulation of GC progression remain unclear. We detected that miR-150-5p overexpression in MGC-803 cells blocked the proliferation rate (Fig. 6a), DNA synthesis (Fig. 6b), and migration and invasion abilities (Fig. 6c), while miR-150-5p blockade had the opposite effects on GC cells. To further address whether circLMTK2 executes its function by interacting with miR-150-5p, we cotransfected miR-150-5p mimics and circLMTK2 expression constructs into GC cells. The effects on GC cell growth and motility suppression induced by miR-150-5p were reversed when the cells overexpressed circLMTK2 (Fig. 6a-c and Additional file 1: Figure S5A-D). These data reveal that circLMTK2 promotes GC cell growth and metastasis by sponging miR-150-5p.
We previously studied that c-Myc-associated circular RNAs, and we noted that circLMTK2 overexpression could increase c-Myc expression at both the mRNA and protein levels, while circLMTK2 knockdown and miR-150-5p overexpression exerted the opposite effects on GC cells (Fig. 7a-b). As expected, we validated that c-Myc could be a potential target gene of miR-150-5p using the TargetScan database. 3′-UTR luciferase reporter assays further confirmed that miR-150-5p could bind directly to a site in the 3′-UTR of c-Myc (Fig. 7c). Although a correlation analysis showed a weak positive correlation between circLMTK2 expression and c-Myc mRNA levels in cancerous tissues (Fig. 7d), it also indicated that circLMTK2 could indirectly regulate c-Myc expression to some extent. These results suggested that c-Myc might be a putative target of miR-150-5p.
circLMTK2 expression is upregulated in GC
To assess the clinical impact of the circLMTK2 expression on GC, we collected clinical data from the aforementioned patients. As presented in Fig. 7e, circLMTK2 levels were significantly higher in tumour tissues than in non-cancerous tissues. More importantly, increased circLMTK2 expression in GC tissues was significantly correlated with poor prognosis in GC patients, as shown by the Kaplan-Meier survival curve using the median value of circLMTK2 expression as the cut-off value (median survival of 16 months vs 45 months, P = 0.0114, log-rank test; Fig. 7f). Significantly higher numbers of late T stage (T4) or poor TNM stage (III-IV) tumours and cases of positive lymph node status were found in the high-circLMTK2 group (P < 0.05, Table 1). No significant differences were observed in the other clinical and pathological characteristics between the high and low circLMTK2 groups. We also performed univariate and multivariate Cox proportional hazards analyses, which included several known prognostic markers (sex, age, tumour size, differentiation grade, TNM stage, Lauren classification, lymphatic invasion and nerve invasion). The results showed that circLMTK2 level and TNM stage were independent prognostic factors for OS in patients with GC (Table 2).
Recently, with the advent of next-generation sequencing, numerous circular RNAs have been identified from various animal genomes. Many of these highly stable circRNAs are abundantly expressed and play a role in many diseases, especially in tumours, via acting as miRNA sponges, decoying proteins, and affecting translation . RNA-seq results provide useful information for revealing general circRNA expression trends and helping select candidate circRNAs for further research. In this study, we successfully identified thousands of circRNAs in human GC tissues and normal gastric tissues, and hundreds of them were differentially expressed. Nevertheless, the expression of each candidate circRNA still needs to be verified in a cohort of clinical samples and cell lines because most of the circRNAs have low abundance (Additional file 1: Figure S2).
circRNAs have been characterized as vital factors in cancer biology [21, 22]. It was previously reported that circRNAs could act as a miRNA “sponge” to decrease the abundance of miRNAs . Shen et al. conducted a study showing that circRNA_001569 could significantly increase cell viability and inhibit cell apoptosis in GC via the miR-145/NR4A2 axis . Another research group found that the novel circRNA-100,269 could target miR-630, leading to the inhibition of tumour cell growth . circRNA_0023642 was shown to induce apoptosis and suppress cell proliferation, migration and invasion in GC via the EMT signalling pathway . Additionally, circLARP4 was found to sponge miR-424, thus inhibiting the biological behaviours of GC . Several studies revealed that circRNAs might also be valuable factors for the diagnosis of GC. Li et al. have assessed the plasma circRNA_002059 levels and determined a qRT-PCR method for detecting circRNA levels in GC patients . Lower levels of circRNA_0000190 were also found in 104 gastric patient plasma samples and served as a novel non-invasive diagnostic biomarker for GC . One recent study showed that GC cell-derived exosomes could promote the transformation of preadipocytes into brown-like cells by delivering ciRS-133 to suppress miR-133 and activate PRDM16 . These studies show that the primary function of circRNAs is to deregulate miRNAs in pathogenic processes, which is important for demonstrating their potential role in cancer .
Here, we demonstrated that circLMTK2, which is an exonic circRNA originating from exons 10 and 11 of LMTK2 mRNA, was upregulated in GC tissues and predominantly localized in the cytoplasm. circLMTK2 is derived from the human lemur tyrosine kinase 2 gene, which is known as LMTK2. Previous studies have focused mainly on the biological function of LMTK2 in prostate cancer. The rs6465657 variant of LMTK2 was evidently related to the development of prostate cancer . Additionally, LMTK2 has been identified as a driver mutation in lung adenocarcinoma by a large-scale RNA-seq analysis .. Recently, He et al. reported that circLMTK2 could act as a novel tumour suppressor in GC; however, they did not explain the source of circLMTK2 or submit any RNA-seq data in their study . In addition, we validated the expression of circLMTK2 in our 25-patient cohort by qRT-PCR using He’s primer (Additional file 1: Table S3). We further confirmed that circLMTK2 was indeed upregulated in tumours compared with that in non-cancerous tissues (Additional file 1: Figure S4). Further analyses were conducted to determine the biological function and molecular mechanism of circLMTK2 in GC cells. We found that circLMTK2 promoted GC cell proliferation and tumourigenicity in vitro and in vivo. Furthermore, ectopic circLMTK2 expression effectively enhanced GC cell migration and invasion in vitro and tumour metastasis in vivo.
circRNAs have been suggested to function as miRNA sponges. We also demonstrated that circLMTK2 could bind to miR-150-5p and that five miR-150-5p-binding elements in circLMTK2 are essential for their interaction. MiR-150 has been implicated as either an oncogene or a tumour suppressor in various types of solid tumours [35,36,37,38]. MiR-150-5p was first considered as the main miRNA in immune and haematopoietic cells . It has recently been shown that Myc upregulates LIN28 expression in myeloid cells, which inhibits miR-150 maturation from its precursor pri-miR-150 . However, the function of miR-150-5p in the pathogenesis of GC is unknown. We also demonstrated that miR-150-5p inhibited GC cell proliferation and motility in vitro. More importantly, we identified that circLMTK2 exerts oncogenic functions by sponging miR-150-5p.
Intriguingly, we found that c-Myc is a potential functional target of miR-150-5p. c-Myc is a key basic helix-loop-helix leucine zipper transcription factor that is frequently upregulated in various human cancers, including GC . Its overexpression contributes to malignant transformation by regulating the expression of a number of genes participating in multiple aspects of tumourigenesis, such as cell cycle progression, cell invasion, migration, metastasis, and angiogenesis [42,43,44,45,46]. c-Myc has been implicated in controlling miRNA expression, and c-Myc-regulated miRNAs affect virtually all aspects of the hallmarks of Myc-driven diseases [47, 48]. Increasing evidence has indicated that there is significant crosstalk between c-Myc and miRNAs. In fact, c-Myc regulates the expression of a number of miRNAs, resulting in widespread miRNA repression, and c-Myc is regulated by miRNAs, leading to sustained Myc activity . Our results showed that circLMTK2 could regulate c-Myc expression by sponging miR-150-5p, indicating that c-Myc may act as an essential component of the regulatory circuit and providing further mechanistic evidence to support the notion that c-Myc is a promising therapeutic target in the treatment of GC.
In conclusion, our study reveals that circLMTK2 is upregulated in GC tissues and that high circLMTK2 expression is associated with poor prognosis, lymph node metastasis and poor TNM stage in GC patients. Moreover, enhanced circLMTK2 expression promoted GC cell growth and motility in vitro and in vivo through sponging miR-150-5p to upregulate c-Myc. Consequently, circLMTK2 may have considerable potential as a prognostic predictor and therapeutic target for GC.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its Additional files.
Cell Counting Kit-8
CircRNA in vivo precipitation
MYC proto-oncogene, bHLH transcription factor
Fluorescence in situ hybridization
Lemur tyrosine kinase 2
Quantitative real-time polymerase chain reaction
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69(1):7–34.
Siegel R, Ma J, Zou Z, Jemal A. Cancer statistics, 2014. CA Cancer J Clin. 2014;64(1):9–29.
Chen D, Wang H, Chen J, Li Z, Li S, Hu Z, Huang S, Zhao Y, He X. MicroRNA-129-5p regulates glycolysis and cell proliferation by targeting the glucose transporter SLC2A3 in gastric Cancer cells. Front Pharmacol. 2018;9:502.
Guarnerio J, Bezzi M, Jeong JC, Paffenholz SV, Berry K, Naldini MM, Lo-Coco F, Tay Y, Beck AH, Pandolfi PP. Oncogenic role of fusion-circRNAs derived from Cancer-associated chromosomal translocations. Cell. 2016;165(2):289–302.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.
Weng W, Wei Q, Toden S, Yoshida K, Nagasaka T, Fujiwara T, Cai S, Qin H, Ma Y, Goel A. Circular RNA ciRS-7-a promising prognostic biomarker and a potential therapeutic target in colorectal Cancer. Clinical cancer research : an official journal of the American Association for Cancer Research. 2017;23(14):3918–28.
Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, Sun W, Dou K, Li H. Circular RNA: a new star of noncoding RNAs. Cancer Lett. 2015;365(2):141–8.
Memczak S, Papavasileiou P, Peters O, Rajewsky N. Identification and characterization of circular RNAs as a new class of putative biomarkers in human blood. PLoS One. 2015;10(10):e0141214.
Shang Q, Yang Z, Jia R, S G. The novel roles of circRNAs in human cancer. Mol Cancer. 2019;18(1):6.
Salzman J, Chen RE, Olsen MN, Wang PL, Brown PO. Cell-type specific features of circular RNA expression. PLoS Genet. 2013;9(9):e1003777.
Salzman J. Circular RNA expression: its potential regulation and function. Trends in genetics : TIG. 2016;32(5):309–16.
Vo JN, Cieslik M, Zhang Y, Shukla S, Xiao L, Zhang Y, Wu YM, Dhanasekaran SM, Engelke CG, Cao X, et al. The landscape of circular RNA in Cancer. Cell. 2019;176(4):869–81 e813.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Han D, Li J, Wang H, Su X, Hou J, Gu Y, Qian C, Lin Y, Liu X, Huang M, et al. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017;66(4):1151–64.
Su X, Wang H, Ge W, Yang M, Hou J, Chen T, Li N, Cao X. An in vivo method to identify microRNA targets not predicted by computation algorithms: p21 targeting by miR-92a in Cancer. Cancer Res. 2015;75(14):2875–85.
Liu H, Liu Y, Bian Z, Zhang J, Zhang R, Chen X, Huang Y, Wang Y, Zhu J. Circular RNA YAP1 inhibits the proliferation and invasion of gastric cancer cells by regulating the miR-367-5p/p27 (Kip1) axis. Mol Cancer. 2018;17(1):151.
Chen X, Xu X, Pan B, Zeng K, Xu M, Liu X, He B, Pan Y, Sun H, Wang S. miR-150-5p suppresses tumor progression by targeting VEGFA in colorectal cancer. Aging. 2018;10(11):3421–37.
Shu L, Wang Z, Wang Q, Wang Y, Zhang X. Signature miRNAs in peripheral blood monocytes of patients with gastric or breast cancers. Open Biol. 2018:8(10).
Lu W, Zhang H, Niu Y, Wu Y, Sun W, Li H, Kong J, Ding K, Shen HM, Wu H, et al. Long non-coding RNA linc00673 regulated non-small cell lung cancer proliferation, migration, invasion and epithelial mesenchymal transition by sponging miR-150-5p. Mol Cancer. 2017;16(1):118.
Li Y, Zheng F, Xiao X, Xie F, Tao D, Huang C, Liu D, Wang M, Wang L, Zeng F, et al. CircHIPK3 sponges miR-558 to suppress heparanase expression in bladder cancer cells. EMBO Rep. 2017;18(9):1646–59.
LH J, DW S, Hou JC, Ji ZL. CircRNA: a novel type of biomarker for cancer. Breast Cancer. 2018;25(1):1–7.
Meng S, Zhou H, Feng Z, Xu Z, Tang Y, Li P, Wu M. CircRNA: functions and properties of a novel potential biomarker for cancer. Mol Cancer. 2017;16(1):94.
Chen L-L, Yang L. Regulation of circRNA biogenesis. RNA Biol. 2015;12(4):381–8.
Shen F, Liu P, Xu Z, Li N, Yi Z, Tie X, Zhang Y, Gao L. CircRNA_001569 promotes cell proliferation through absorbing miR-145 in gastric cancer. J Biochem. 2019;165(1):27–36.
Zhang Y, Liu H, Li W, Yu J, Li J, Shen Z, Ye G, Qi X, Li G. CircRNA_100269 is downregulated in gastric cancer and suppresses tumor cell growth by targeting miR-630. Aging. 2017;9(6):1585–94.
Zhou LH, Yang YC, Zhang RY, Wang P, Pang MH, Liang LQ. CircRNA_0023642 promotes migration and invasion of gastric cancer cells by regulating EMT. Eur Rev Med Pharmacol Sci. 2018;22(8):2297–303.
Zhang J, Liu H, Hou L, Wang G, Zhang R, Huang Y, Chen X, Zhu J. 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.
Li P, Chen S, Chen H, Mo X, Li T, Shao Y, Xiao B, Guo J. Using circular RNA as a novel type of biomarker in the screening of gastric cancer. Clinica chimica acta; international journal of clinical chemistry. 2015;444:132–6.
Chen S, Li T, Zhao Q, Xiao B, Guo J. Using circular RNA hsa_circ_0000190 as a new biomarker in the diagnosis of gastric cancer. Clinica chimica acta; international journal of clinical chemistry. 2017;466:167–71.
Wang KW, Dong M. Role of circular RNAs in gastric cancer: recent advances and prospects. World journal of gastrointestinal oncology. 2019;11(6):459–69.
Dragomir M, Calin GA. Circular RNAs in Cancer - lessons learned from microRNAs. Front Oncol. 2018;8:179.
Harries LW, Perry JR, McCullagh P, Crundwell M. Alterations in LMTK2, MSMB and HNF1B gene expression are associated with the development of prostate cancer. BMC Cancer. 2010;10:315.
Seo JS, Ju YS, Lee WC, Shin JY, Lee JK, Bleazard T, Lee J, Jung YJ, Kim JO, Shin JY, et al. The transcriptional landscape and mutational profile of lung adenocarcinoma. Genome Res. 2012;22(11):2109–19.
He J, Chen J, Ma B, Jiang L, Zhao G. CircLMTK2 acts as a novel tumor suppressor in gastric cancer. Biosci Rep. 2019:39(5).
Li YJ, Zhang YX, Wang PY, Chi YL, Zhang C, Ma Y, Lv CJ, Xie SY. Regression of A549 lung cancer tumors by anti-miR-150 vector. Oncol Rep. 2012;27(1):129–34.
Srivastava SK, Bhardwaj A, Singh S, Arora S, Wang B, Grizzle WE, Singh AP. MicroRNA-150 directly targets MUC4 and suppresses growth and malignant behavior of pancreatic cancer cells. Carcinogenesis. 2011;32(12):1832–9.
Wu Q, Jin H, Yang Z, Luo G, Lu Y, Li K, Ren G, Su T, Pan Y, Feng B, et al. MiR-150 promotes gastric cancer proliferation by negatively regulating the pro-apoptotic gene EGR2. Biochem Biophys Res Commun. 2010;392(3):340–5.
Zhang J, Luo N, Luo Y, Peng Z, Zhang T, Li S. microRNA-150 inhibits human CD133-positive liver cancer stem cells through negative regulation of the transcription factor c-Myb. Int J Oncol. 2012;40(3):747–56.
He Y, Jiang X, Chen J. The role of miR-150 in normal and malignant hematopoiesis. Oncogene. 2014;33(30):3887–93.
Jiang X, Huang H, Li Z, Li Y, Wang X, Gurbuxani S, Chen P, He C, You D, Zhang S, et al. Blockade of miR-150 maturation by MLL-fusion/MYC/LIN-28 is required for MLL-associated leukemia. Cancer Cell. 2012;22(4):524–35.
Calcagno DQ, Leal MF, Assumpcao PP, Smith MA, Burbano RR. MYC and gastric adenocarcinoma carcinogenesis. World J Gastroenterol. 2008;14(39):5962–8.
Gabay M, Li Y, Felsher DW. MYC activation is a hallmark of cancer initiation and maintenance. Cold Spring Harbor perspectives in medicine. 2014:4(6).
Qian K, Mao B, Zhang W, Chen H. MicroRNA-561 inhibits gastric cancercell proliferation and invasion by downregulating c-Myc expression. Am J Transl Res. 2016;8(9):3802–11.
Dang CV. MYC on the path to cancer. Cell. 2012;149(1):22–35.
Ben-Yosef T, Yanuka O, Halle D, Benvenisty N. Involvement of Myc targets in c-myc and N-myc induced human tumors. Oncogene. 1998;17(2):165–71.
Kalkat M, De Melo J, Hickman KA, Lourenco C, Redel C, Resetca D, Tamachi A, Tu WB, Penn LZ. MYC Deregulation in Primary Human Cancers. Genes. 2017:8(6).
Slack GW, Gascoyne RD. MYC and aggressive B-cell lymphomas. Adv Anat Pathol. 2011;18(3):219–28.
Ott G, Rosenwald A, Campo E. Understanding MYC-driven aggressive B-cell lymphomas: pathogenesis and classification. Blood. 2013;122(24):3884–91.
Tao J, Zhao X, Tao J. C-MYC-miRNA circuitry: a central regulator of aggressive B-cell malignancies. Cell Cycle. 2014;13(2):191–8.
We thank Prof. Lu Chun (Department of Microbiology and Immunology, Nanjing Medical University, China) for kindly providing the lentiviral packaging system consisting of psiCHECK-2, psPAX2, and pMD2.G.
This work was supported by grants from 333 Project for High-level Talents Cultivation of Jiangsu Province (BRA2017153).
Ethics approval and consent to participate
The research was approved by the medical ethics committee of Northern Jiangsu People’s Hospital.
Consent for publication
All authors have agreed to publish this manuscript.
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers and RNA sequences used in this study. Figure S1. Validation of 4 differentially expressed circRNA candidates in GC tissues using qRT-PCR in 25 paired GC tissues and matched normal gastric tissues. Figure S2. Expression levels of circLMTK2 in 8 GC cell lines. Figure S3. Schematic illustration of the circLMTK2 overexpression vector. Figure S4. Detection of circLMTK2 in GC tissues by qRT-PCR using He’s primer in 25 paired GC tissues and matched normal gastric tissues. Figure S5. circLMTK2 promotes GC cell growth and metastasis by sponging miR-150-5p.
35,350 distinct circRNA candidates.
142 differentially expressed circRNAs in GC compared with those in normal tissues.
About this article
Cite this article
Wang, S., Tang, D., Wang, W. et al. circLMTK2 acts as a sponge of miR-150-5p and promotes proliferation and metastasis in gastric cancer. Mol Cancer 18, 162 (2019). https://doi.org/10.1186/s12943-019-1081-4
- Gastric cancer