- Open Access
The circular RNA circDLG1 promotes gastric cancer progression and anti-PD-1 resistance through the regulation of CXCL12 by sponging miR-141-3p
Molecular Cancer volume 20, Article number: 166 (2021)
Dysregulation of circular RNAs (circRNAs) plays an important role in the development of gastric cancer; thus, revealing the biological and molecular mechanisms of abnormally expressed circRNAs is critical for identifying novel therapeutic targets in gastric cancer.
A circRNA microarray was performed to identify differentially expressed circRNAs between primary and distant metastatic tissues and between gastric cancer tissues sensitive or resistant to anti-programmed cell death 1 (PD-1) therapy. The expression of circRNA discs large homolog 1 (DLG1) was determined in a larger cohort of primary and distant metastatic gastric cancer tissues. The role of circDLG1 in gastric cancer progression was evaluated both in vivo and in vitro, and the effect of circDLG1 on the antitumor activity of anti-PD-1 was evaluated in vivo. The interaction between circDLG1 and miR-141-3p was assessed by RNA immunoprecipitation and luciferase assays.
circDLG1 was significantly upregulated in distant metastatic lesions and gastric cancer tissues resistant to anti-PD-1 therapy and was associated with an aggressive tumor phenotype and adverse prognosis in gastric cancer patients treated with anti-PD-1 therapy. Ectopic circDLG1 expression promoted the proliferation, migration, invasion, and immune evasion of gastric cancer cells. Mechanistically, circDLG1 interacted with miR-141-3p and acted as a miRNA sponge to increase the expression of CXCL12, which promoted gastric cancer progression and resistance to anti-PD-1-based therapy.
Overall, our findings demonstrate how circDLG1 promotes gastric cancer cell proliferation, migration, invasion and immune evasion and provide a new perspective on the role of circRNAs during gastric cancer progression.
Gastric cancer is one of the most commonly diagnosed cancers and a cause of cancer-related death worldwide . Although great progress has been made in the treatment of gastric cancer in recent years, the prognosis is still poor, especially for advanced-stage patients. Tumor progression and distant metastasis are the main causes of death. For metastatic gastric cancer patients, chemotherapy is still the main therapeutic strategy. In recent years, increasing evidence has shown that immune evasion is essential for tumor survival and development [2, 3]. It has been reported that in the tumor microenvironment, tumor cells can recruit immunosuppressive cells, such as CD4+ T cells, to disrupt the cytotoxic functions of CD8+ T cells [4, 5]. In addition, programmed death ligand 1 (PD-L1), a B7 family ligand, can bind to its receptor programmed death-1 (PD-1) to influence tumor-specific T cells, induce apoptosis and inhibit the activity of CD8+ T cells, leading to immune evasion in tumors. Clinically, an increasing number of studies have shown that blocking the PD-1 checkpoint with an anti-PD-1 antibody is an efficient immunotherapy approach in different cancers [6, 7]. Our group confirmed the clinical benefit of PD-1 antibody therapy in gastric cancer patients . However, most gastric cancer patients are resistant to anti-PD-1 therapy, and only some patients benefit from this treatment . PD-L1 has been proven to be a predictive biomarker of anti-PD-1 efficacy in several tumor types. However, the predictive role of PD-L1 expression in gastric cancer is controversial . The KENOTE-061 and KENOTE-062 trials showed better survival in patients with PD-L1-positive tumors after pembrolizumab treatment [11, 12]. On the other hand, data from the Checkmate032, JAVELIN Gastric 300, and ATTRCTION-2 trials did not support the concept of PD-L1 positivity as a predictive biomarker of anti-PD-1 efficacy [13,14,15]. Therefore, a better understanding of the molecular mechanisms that contribute to gastric cancer progression and immune evasion is critical for developing effective therapeutic strategies for this disease.
Circular RNAs (circRNAs) are a class of recently identified noncoding RNAs or protein-coding RNAs that are characterized by covalently closed loops without 5′ to 3′ polarity or polyadenylation (polyA) tails . Previous reports have found that circRNAs are generally derived from the back-splicing of pre-mRNA transcripts and are mainly located in the cytoplasm [17, 18]. Studies have demonstrated that circRNAs are conserved, stable and abundant in different tumor cells and tissues . Increasing studies have indicated that circRNAs participate in multiple physiological and pathological processes by modulating gene expression, apoptosis, the cell cycle, cell migration and invasion . Our previous study identified a circRNA signature that can predict postoperative recurrence in stage II-III colon cancer . Mechanistically, circRNAs can exert their functions in different ways, such as sponging miRNAs, interacting with RNA-binding proteins, and translating proteins [22,23,24]. For instance, our previous study found that circUBXN7 expression was downregulated in bladder cancer, and forced circUBXN7 expression could suppress cell growth and invasion by sponging miR-1247 to enhance B4GALT3 expression in bladder cancer cells . Rong et al.  demonstrated that circPSMC3 could inhibit the proliferation and metastasis of gastric cancer. Han et al.  reported that circMTO1 modulated the progression of hepatocellular carcinoma through the regulation of p21 expression by sponging miR-9. Hsiao et al.  verified that the circRNA CCDC66 promoted the progression and metastasis of colon cancer. Although the molecular mechanisms of dysregulated circRNA-associated pathways have been extensively explored, the role of circRNAs in gastric cancer progression and immune evasion remains poorly understood.
In this study, to identify circRNAs involved in gastric cancer progression and immune evasion, we performed a circRNA array using primary and distant metastatic tumor tissues as well as tissues sensitive or resistant to anti-PD-1 therapy. We found that the circRNA circDLG1 (hsa_circ_0008583), which is derived from the DLG1 gene, was significantly upregulated in distant metastatic tissues and primary gastric cancer tissues resistant to anti-PD-1 therapy. More importantly, circDLG1 expression was significantly associated with an aggressive tumor phenotype and adverse prognosis in gastric cancer patients who received anti-PD-1 therapy. In addition, ectopic expression of circDLG1 promoted in vitro cell proliferation, invasion, immune evasion, and in vivo tumorigenesis and metastasis in immunocompetent mice. Mechanistically, circDLG1 could directly interact with miR-141-3p, acting as a miRNA sponge to increase the expression of the miR-141-3p target gene chemokine 12 (CXCL12), thereby promoting the progression of gastric cancer. Thus, circDLG1 might be a promising therapeutic target and biomarker for gastric cancer.
Human tissue samples and cell lines
Fresh-frozen or paraffin-embedded tissues were obtained from advanced-stage gastric cancer patients who underwent gastroscopy and biopsy, fine needle biopsy, or palliative surgery and received anti-PD-1 therapy at Sun Yat-sen University Cancer Center from August 2018 to October 2019. These patients were all included in a real-world study of PD-1 antibody therapy in gastric cancer (No: NCT04086888). A total of 126 gastric cancer patients were enrolled. Eighty-two of these patients had primary tumor tissues available. Seventy-three patients could be evaluated for the efficacy of PD-1 antibody therapy, among which 30 patients had tissues available (including primary tissues, adjacent normal tissues, and distant metastatic tissues). These patients were evaluated every 3-4 cycles of therapy and followed up regularly. The response rate was evaluated based on RECIST 1.1 guidelines. The tumor response included complete response (CR), partial response (PR), stable disease (SD), and progressive disease (PD). The objective response rate (ORR) was defined as the percentage of patients who achieved CR and PR. The clinical and pathological parameters, including age, sex, tumor size, tumor cell differentiation, peritoneal metastasis, and Lauren’s classification, were collected from patient records. This study was approved by the ethics committee of the Sun Yat-sen University Cancer Center, and informed consent was obtained from all patients. Overall survival (OS) was calculated from the start of PD-1 antibody treatment to the date of death or last contact, and progression-free survival (PFS) was calculated from the start of PD-1 antibody treatment to the date of progression or death.
Human gastric cancer cell lines (HGC27, BGC823, MKN45, MKN28, SGC7901, AGS), the normal gastric epithelial cell line GES-1, the murine gastric cancer cell line MFC and human embryonic kidney (HEK) 293 T cells were purchased from the Shanghai Institute of Cell Biology, Chinese Academy of Science (Shanghai, China). All cell lines were cultured according to the provider’s instructions.
CircRNA microarray assays
To identify the potential circRNAs involved in gastric cancer progression and anti-PD-1 resistance, 3 patients who responded to anti-PD-1 (experienced PR) therapy with PFS of more than 10 months and 6 patients who were nonresponsive to anti-PD-1 (2 SD and 4 PD) therapy with PFS of less than 5 months were selected for circRNA microarray analysis.
The circRNA microarray was performed by Kangchen Biotech (Shanghai, China). All sample preparation and microarray hybridization steps were conducted according to the instructions of Arraystar (Rockville, MD, USA). In brief, circRNAs were enriched by removing the linear RNAs with RNase R treatment, and the enriched circRNAs were then amplified and transcribed into fluorescent cRNA with a random priming method (Arraystar Super RNA Labeling Kit). Then, the labeled cRNAs were hybridized onto Human circRNA Array V2 (8 × 15 k, Arraystar). Finally, the arrays were scanned using the Agilent Scanner G2505C and analyzed with Agilent Feature Extraction software (version 126.96.36.199). The random variance model was used to identify the differentially expressed genes. The paired t-test was used to calculate the P value. The threshold for differentially expressed genes was fold change ≥2.0 and a P value ≤0.05.
Quantitative real-time polymerase chain reaction (qRT–PCR)
Total RNA was isolated from human gastric cancer tissues and cell lines using TRIzol reagent (Sigma–Aldrich, St. Louis, MO, USA) according the manufacturer’s instructions. Reverse transcription of mRNA and miRNA was performed using random primers and stem–loop primers, respectively. qRT–PCR was conducted using a TaqMan Universal Master Mix II kit on a Bio–Rad CFX96 qPCR system, and fold changes were determined by using the relative quantification 2-ΔΔCT method. The nuclear and cytoplasmic fractions of cells were separated using the PARIS Kit (Life Technologies) according to the manufacturer’s instructions. RNA was extracted from both fractions. Then, qRT–PCR was performed to determine the expression ratios of specific RNA molecules between the nuclear and cytoplasmic fractions. GAPDH and U6 served as cytoplasmic and nuclear markers, respectively. The primers used are presented in Additional file 1: Table S1.
Western blotting analysis
Western blotting analysis was performed according to a previously described method . Briefly, proteins were extracted from gastric cancer cell lines, and the protein concentration was calculated using the Pierce BCA protein assay kit (Thermo Scientific, Rockford, IL, USA). The following antibodies were used in this study: anti-GAPDH (CST, #2118 L) and anti-CXCL12 (CST, #3740S).
Immunohistochemistry (IHC) analysis
The IHC analysis was performed according to the method we described previously . Briefly, the paraffin-embedded tissue blocks were cut into 4-μm slides. Immunostaining images were captured using a microscope (Leica, Germany). The immunoreactivity in each tissue section was assessed by two pathologists, and the degree of positivity was evaluated according to the percentage of positive tumor cells. The following antibodies were used: anti-CXCL12 (CST, 97958S; 1:200 dilution), anti-PD-L1 (22C3 pharmDX; Dako, Carpinteria, CA, USA), and anti-CD33 (clone PWS44; 1:200 dilution; Leica Biosystems, Nussloch, Germany). PD-L1 positivity was defined as a membrane staining intensity ≥1% in the tumor cells or tumor infiltrating immune cells. The median value of PD-L1 expression was used as the cutoff value to stratify high and low PD-L1 expression.
The detailed flow cytometry procedure was performed as previously described . The following antibodies were used: anti-mouse IFNγ (BD Biosciences), anti-mouse CD45 (BD Biosciences), anti-mouse Gr-1 (TONBO Biosciences), anti-mouse Ly6G (TONBO Biosciences), anti-mouse ly6C (BD Biosciences), anti-mouse CD11b (TONBO Biosciences), anti-mouse F4/80 (BioLegend), anti-human CD11b (TONBO Biosciences), anti-human CD33 (BD Biosciences), and anti-human HLA-DR (BD Biosciences). Detailed information on the antibodies is provided in Additional file 2: Table S2. Briefly, the cells were digested and suspended as single cells, washed with PBS, and then resuspended in cell stabilizing buffer (BioLegend cat. No. 420201); the supernatant was discarded, and the sample was centrifuged at 350 g for 5 min. Then, 5 μl of Human TruStain FcX™ (Fc Receptor Blocking Solution, BioLegend Cat. No. 422301) and the antibodies were added and cultured for 15-20 min at room temperature. Finally, the cells were assessed by flow cytometry. For intracellular staining of IFNγ, fixation buffer (BioLegend Cat. No. 420801) was added and incubated for 20 min at room temperature. Then, the cells were washed with Intracellular Staining Perm Wash Buffer (BioLegend Cat. No. 421002), resuspended and centrifuged. IFNγ antibody was added, and the cells were cultured for 20 min and finally assessed by flow cytometry.
Fluorescence in situ hybridization (FISH) assays
The FISH assay was performed according to previously described methods . Briefly, cells were hybridized with Cy3-labeled circDLG1 probe and Cy5-labeled miR-141-3p probe (GenePharma, China) for 12 h at 37 °C. Then, the nuclei were counterstained using 4′,6-diamidino-2-phenylindole (DAPI) (Yeasen, Shanghai, China). Finally, the cells were viewed and captured under a ZEISS LSM800 confocal microscope (Carl Zeiss AG, Germany).
Cell proliferation, colony formation, migration and invasion assays
Cell counting kit-8 (CCK-8) and colony formation assays were performed to assess the cell proliferation ability. Transwell assays were conducted to evaluate cell migration and invasion abilities. The details of these assays were described in our previous study .
RNA immunoprecipitation (RIP) and luciferase activity assays
The RIP assay was conducted using the Magna RIP RNA-bing Protein Immunoprecipitation kit (Millipore, USA) according to the provider’s protocol. In brief, cell lysates were cultured with Dynabeads-coated IgG antibody (Millipore, USA) or AGO2 antibody (Cell Signaling Technology, USA) for 12 h at 4 °C. The purified RNA was subjected to qRT–PCR to detect the enriched circDLG1 and miRNA.
For the luciferase activity assay, potential binding sites were predicted using StarBase v3.0 and TargetScanHuman 7.2. Gastric cells were cotransfected with pGL-luc-circDLG1, pGL-luc-CXCL12 3′-UTR, and miR-141-3p mimics or negative control mimics for 48 h, and luciferase activity was detected using the dual-luciferase reporter assay system (Promega, USA) according to the manufacturer’s instructions.
circDLG1 knockdown or overexpression transfection experiment
Small hairpin RNAs (shRNAs) targeting the junction region of the circDLG1 sequence and the circDLG1-overexpressing lentiviral vector were synthesized by Geneseed Biotech Co., Ltd. (Guangzhou, China). Gastric cancer cell lines were transfected with circDLG1 shRNA or the circDLG1-overexpressing lentiviral vector following the manufacturer’s instructions according to a previously described method .
In vivo tumorigenesis and metastasis assays
The in vivo tumorigenesis and metastasis experiments in BALB/c nude mice and C57BL/6 mice were approved by the Animal Experiment Ethics Committee of Sun Yat-sen University Cancer Center. The mice were purchased from the Shanghai Institute of Material Medicine (Shanghai, China) and were fed according to the provider’s instructions. The experiments were conducted as described in our previous study . Briefly, to evaluate the in vivo tumorigenesis effect of circDLG1, MFC-sh-circDLG1 and MFC-sh-NC cells (1 × 106 cells/mouse) were inoculated subcutaneously into the flanks of two groups of BALB/c nude mice and C57BL/6 mice (ten for each cell group). Tumor size was measured every 4 days, and tumor volume was estimated. After 4 weeks, the mice were sacrificed, and the tumors were excised. To investigate the effect of circDLG1 on tumor metastasis, MFC-sh-circDLG1 and MFC-sh-NC cells (2 × 106 cells/mouse) were inoculated into the tail vein of two groups of BALB/c nude mice and C57BL/6 mice (ten for each cell group). Four weeks later, the mice were sacrificed, and the lungs were excised and paraffin embedded. Consecutive sections (4 μm) were made and stained with hematoxylin-eosin. The micrometastases in the lungs were examined and counted under a dissecting microscope.
Mouse xenograft anti-PD-1 therapy experiment
To evaluate the effect of CXCL12 on immune evasion and the efficacy of anti-PD-1 therapy in vivo, an in vivo xenograft experiment was conducted in C57BL/6 mice. In brief, four groups of mice were implanted subcutaneously in the left flanks with MFC-sh-NC, and two groups of mice were implanted subcutaneously in the left flanks with MFC-sh-CXCL12 cells (ten mice for each group, 1 × 106 cells/mouse). The MFC-sh-NC implanted mice were treated with 0.9% normal saline (NS), anti-PD-1 antibody, AMD3100 (a CXCR4 inhibitor), and the combination of AMD3100 and anti-PD-1 antibody. MFC-sh-CXCL12-implanted mice were treated with 0.9% normal saline (NS) and an anti-PD-1 antibody. A mouse anti-PD-1 antibody was purchased from Bio X Cell (West Lebanon, NH, USA). AMD3100 was purchased from Selleck (Shanghai, China) and dissolved in phosphate-buffered saline (PBS). The mice were intraperitoneally injected with 0.9% NS, AMD3100, or anti-PD-1 antibody according to a previously described method . The mice were sacrificed at 4 weeks or before becoming moribund. The survival time was defined as the date of first therapy to the date of sacrifice.
RNA in situ hybridization (ISH), sphere formation assay, and RNA sequencing (RNA-seq) analysis
The methods for ISH, sphere formation assay and RNA-seq are described in Additional file 3: Supplementary materials and methods.
All data are presented as the mean ± standard deviation unless otherwise noted. Statistical analyses were performed with SPSS 17.0 software (SPSS Inc., Chicago, IL, USA) or GraphPad Prism 7 (GraphPad Software, Inc., La Jolla, CA, USA). For comparisons between groups, Student’s t-test, chi-squared test, one-way ANOVA, and Pearson correlation analysis were used, as appropriate. PFS was analyzed using the Kaplan–Meier method with the log-rank test, and Cox’s proportional hazard regression model was performed to evaluate the independent prognostic indicators. A P-value of < 0.05 was considered significant.
Identification of gastric cancer progression- and immune evasion-associated circRNAs through microarray
To identify the circRNA expression profiles involved in gastric cancer progression and immune evasion, we performed a circRNA microarray analysis on 9 paired gastric cancer tissues (primary tissues and paired distant metastatic lesions) from patients who received anti-PD-1 treatment. The circus plot shows the distribution and expression profiles of the detected and differentially expressed circRNAs on human chromosomes, as well as potential miRNAs interacting with differentially expressed circRNAs (Fig. 1a). The heatmap presents the significantly expressed circRNAs in gastric cancer and paired distant metastatic tissues (Fig. 1b). To identify the potential circRNAs involved in anti-PD-1 resistance, we compared circRNA profiling in the primary tissues of 3 patients who were responsive (experienced PR) and 4 patients who were nonresponsive (experienced PD) to anti-PD-1 therapy (Fig. 1c). To select a circRNA that plays a critical role in gastric cancer progression and immune evasion, circRNAs (hsa_circ_0008583, has_circ_0002387) that were upregulated in both distant metastatic lesions and tissues resistant to anti-PD-1 therapy were identified. However, only hsa_circ_0008583 (derived from exons 13, 14, 15 and 16 of DLG1, thus designated circDLG1) was confirmed to be upregulated in distant metastatic lesions and primary gastric cancer tissues resistant to anti-PD-1 therapy in a larger cohort of patients (Fig. 1d and e). Moreover, circDLG1 expression was higher in primary tumor tissues than in adjacent normal tissues, although the difference was not significant (Additional file 4: Fig. S1). Kaplan–Meier analysis showed that circDLG1 expression was significantly associated with PFS in gastric cancer patients treated with anti-PD-1 therapy, and patients with low circDLG1 expression presented with significantly better prognosis than those with high circDLG1 expression (Fig. 1f). The patients were divided into two groups (high expression and low expression) based on the mean circDLG1 expression level from qRT–PCR. circDLG1 expression was significantly associated with tumor size and peritoneal metastasis. However, no association was observed between circDLG1 expression and age, sex, tumor cell differentiation, or Lauren’s classification (Additional file 5: Table S3). Several previous reports indicated that PD-L1 expression and tumor mutation burden (TMB) were associated with the anti-PD-1 response in gastric cancer [8, 9]. We questioned whether circDLG1 expression is associated with PD-L1 expression and TMB. ISH analysis was also used to confirm circDLG1 expression in gastric cancer tissues (blue staining indicates positive expression; red staining indicates negative), and IHC was performed to detect PD-L1 expression in gastric cancer tissues (Fig. 1g). However, no significant correlation was found between circDLG1 expression and PD-L1 expression or TMB (Fig. 1h and i).
Taken together, these results demonstrate that circDLG1 might play an important role in gastric cancer progression and immune evasion.
Characterization of circDLG1 in gastric cancer
According to the circBase database and UCSC, circDLG1 (circRNA ID: hsa_circ_0008583) is derived from exons 13, 14, 15 and 16 of DLG1, which is located on chromosome 3 and has a length of 496 nucleotides (Fig. 2a). CircDLG1 was amplified by using divergent primers and confirmed by Sanger sequencing in AGS and SGC7901 gastric cancer cells (Fig. 2b), and PCR analysis of reverse-transcribed RNA (cDNA) and genomic DNA (gDNA) was performed. The results showed that the divergent primers could amplify products from cDNA but not from gDNA (Fig. 2b). To further confirm the characteristics of circDLG1, we used random hexamer or oligo (dT)18 primers for reverse transcription of RNA from gastric cancer cells. The relative expression of circDLG1 was significantly decreased, while the expression of DLG1 mRNA was not decreased when oligo (dT)18 was used (Fig. 2c). This finding indicates that circDLG1 has no poly A tail. Moreover, circDLG1 was resistant to RNase R (a highly processive 3′ to 5′ exoribonuclease that degrades linear RNAs but not circRNAs), while DLG1 mRNA was sensitive to RNase R (Fig. 2d). In addition, SGC7901 cells were treated with actinomycin D (an inhibitor of transcription), and RNA was extracted at the indicated times for the detection of circDLG1. The results demonstrated that the half-life of circDLG1 was significantly longer than that of DLG1 mRNA (Fig. 2e). Furthermore, the qRT–PCR and FISH results indicated that circDLG1 is mainly located in the cytoplasm (Fig. 2f and g). These results demonstrate that circDLG1 is a stable and abundant circRNA transcript.
CircDLG1 promotes gastric cancer progression
To explore the biological function of circDLG1 in gastric cancer, we first measured circDLG1 expression in gastric cancer cell lines and the gastric normal epithelial cell line GES-1. The results showed that circDLG1 expression was upregulated in gastric cancer cell lines compared with GES-1 gastric normal epithelial cells (Fig. 3a). We designed three shRNA plasmids to target the unique back-splice junction. Back-splice junction-specific shRNAs (sh-1, sh-2, sh-3) effectively knocked down circDLG1 expression in HGC27 and SGC7901 cells (Fig. 3b). In vitro CCK-8 and colony formation assays indicated that knockdown of circDLG1 significantly inhibited the proliferation and colony formation ability of HGC27 and SGC7901 cells (Fig. 3c and d). In addition, transwell assays indicated that knockdown of circDLG1 significantly inhibited the invasion ability of HGC27 and SGC7901 cells (Fig. 3e). Moreover, considering that epithelial-mesenchymal transition (EMT) plays an important role in the process of cancer progression and invasion, we explored the effects of circDLG1 on the EMT and stemness of gastric cancer cells. Immunofluorescence analysis showed that knockdown of circDLG1 significantly increased the expression of E-cadherin while decreasing the expression of N-cadherin (Fig. 3f). qRT–PCR analysis indicated that knockdown of circDLG1 markedly increased the expression of epithelial markers (E-cadherin, α-catenin, β-catenin) while decreasing the expression of mesenchymal markers (N-cadherin, vimentin, Snail, Slug) (Fig. 3g). IHC showed that the expression of E-cadherin decreased, whereas vimentin expression increased in tissues with high expression of circDLG1 (Additional file 6: Fig. S2). In addition, sphere formation assays showed that knockdown of circDLG1 markedly inhibited sphere formation (Fig. 3h), and qRT–PCR analysis indicated that knockdown of circDLG1 significantly inhibited the expression of stemness-associated markers (CD44, CD133, Oct4, CD24, CD155, CD166, Nanog, Sox2) (Fig. 3i).
CircDLG1 promotes tumor progression via increased infiltration of myeloid-derived suppressor cells (MDSCs)
We generated MFC-sh-circDLG1 cell lines through the stable transduction of circDLG1-targeting shRNA vectors into the MFC murine gastric cancer cell line. qRT–PCR analysis confirmed that knockdown of circDLG1 significantly decreased circDLG1 expression in MFC cells (Additional file 7: Fig. S3). To elucidate the function of circDLG1 in vivo, MFC-sh-circDLG1 cells were injected into immunocompetent mice subcutaneously or via the tail vein. The results showed that knockdown of circDLG1 significantly inhibited tumor progression and metastasis compared with the control treatment (Fig. 4a). In contrast, no significant difference was observed in tumor progression or metastasis between the MFC-sh-circDLG1 and MFC-control groups of immunodeficient mice (Fig. 4b). These results demonstrated that circDLG1 promotes gastric cancer progression and metastasis and might be related to tumor immune evasion. IHC analysis showed that MFC-sh-circDLG1-derived tumors were associated with increased CD8+ tumor-infiltrating lymphocytes (TILs) and decreased infiltration of Gr-1- and CD11b-positive MDSCs (Gr-1 and CD11b are mouse MDSC markers) (Fig. 4c, Additional file 8: Fig. S4). Moreover, flow cytometric analysis of immune cells in subcutaneous tumors showed that knocking down circDLG1 significantly increased the number of CD8+ T cells and IFNγ+ cells but significantly decreased the number of MDSCs and granulocytic (G)-MDSCs (Fig. 4d). We further evaluated the relationship between circDLG1 and MDSCs in gastric cancer tissues. The expression of CD33 (a marker of human MDSCs) was detected by IHC, and the results showed a significant correlation between circDLG1 and CD33+ cells (Fig. 4e and f). These results demonstrated that circDLG1 promotes gastric cancer progression by inducing MDSC infiltration.
CircDLG1 upregulates CXCL12 to promote cell progression and immune evasion in gastric cancer
To further explore the underlying molecular mechanism by which circDLG1 promotes gastric cancer progression and immune evasion, gene expression profiling was performed in gastric cancer cells with or without knockdown of circDLG1. SGC7901 cells were transfected with circDLG1 shRNA, and the expression levels of mRNAs were detected by RNA sequencing in these cells (Additional file 9: Fig. S5a). Gene annotation enrichment analysis showed that differentially expressed genes involved in T cell migration, EMT and immune responses were significantly enriched based on Gene Ontology (GO) analysis, and TNF signaling and the PI3K-Akt signaling pathway were significantly enriched based on Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis (Additional file 9: Fig. S5b). The circRNA/mRNA microarray dataset was analyzed, and the correlation between circDLG1 and mRNAs was determined. The results showed that CXCL12 had the highest correlation coefficient. In addition, CXCL12 was one of the most significantly decreased genes by knockdown of circDLG1 in gastric cancer cells based on the RNA sequencing assay. Thus, we speculated that there might be an association between circDLG1 and CXCL12. qRT–PCR analysis showed that knockdown of circDLG1 significantly decreased the expression of CXCL12, whereas ectopic expression of circDLG1 significantly increased CXCL12 expression in gastric cancer cell lines (Additional file 9: Fig. S5c and d). Western blot analysis confirmed that knockdown of circDLG1 significantly decreased the protein level of CXCL12 in SGC7901 and HGC27 gastric cancer cells (Additional file 9: Fig. S5e). Moreover, the decreased level of CXCL12 caused by CXCL12 knockdown could be restored by circDLG1 ectopic expression (Additional file 9: Fig. S5f). Accordingly, the suppressed cell migration ability caused by CXCL12 knockdown could be restored by circDLG1 ectopic expression in gastric cancer cells (Additional file 9: Fig. S5g). Furthermore, circDLG1 expression was significantly associated with CXCL12 in gastric cancer tissues (Additional file 9: Fig. S5h). These data indicated that circDLG1 exerts its effect at least partly by regulating CXCL12 expression.
CircDLG1 functions as a sponge for miR-141-3p to regulate CXCL12 expression
Previous studies have reported that circRNAs can act as miRNA sponges to regulate the expression of protein-coding genes. In the present study, we found that circDLG1 was mainly located in the cytoplasm, which indicated that it might function as a miRNA sponge. Interestingly, by using the online bioinformatics tools miRanda and TargetScan, 46 miRNAs were predicted to be potential targets of circDLG1 (Additional file 10: Table S4). We first explored the occupancy of AGO2 in the region of circDLG1. The results of RNA immunoprecipitation for AGO2 in SGC7901 cells stably expressing Flag-AGO2 and Flag-GFP showed that endogenous circDLG1 pulled down by Flag-AGO2 was significantly enriched, suggesting that circDLG1 is incorporated into the RNA-induced silencing complex (RISC) (Fig. 5a). To identify the miRNAs that can interact with circDLG1, a circDLG1 fragment was amplified and inserted downstream of a luciferase reporter gene. Transfection of circDLG1 shRNA significantly decreased the luciferase activity of the circDLG1 reporter (Fig. 5b). To further identify the key miRNAs that can bind to circDLG1, a luciferase screen with a miRNA mimic library was performed. A total of 46 miRNA mimics were cotransfected with the circDLG1 luciferase reporter into SGC7901 cells (Fig. 5c). The luciferase activity of the circDLG1 reporter was significantly reduced by four miRNAs (miR-141-3p, miR-138-5p, miR-651-3p, miR-155-3p) (Fig. 5d). In addition, RNA precipitation (RIP) analysis with a biotin-labeled circDLG1 probe was conducted to evaluate the key miRNAs that can interact with circDLG1. Compared with the negative control, miR-141-3p, but not the other three miRNAs, was enriched by the biotin-labeled circDLG1 probe (Fig. 5e), and the enrichment of miR-141-3p by circDLG1 was reduced when circDLG1 was knocked down in gastric cancer cells (Fig. 5f). In addition, an inverse association was observed between the expression of circDLG1 and miR-141-3p in gastric cancer tissues (Fig. 5g). Furthermore, the RNA FISH assay demonstrated that miR-141-3p colocalized with both endogenous and exogenous circDLG1 in gastric cancer cells (Fig. 5h). These data demonstrated that circDLG1 could function as a miR-141-3p sponge in gastric cancer cells.
To evaluate whether circDLG1 could regulate CXCL12 expression by sponging miR-141-3p, we first confirmed the regulation of CXCL12 by miR-141-3p. The results showed that ectopic expression of miR-141-3p could significantly decrease the mRNA level of CXCL12 in gastric cancer cells (Fig. 5i). Moreover, Western blot analysis showed that ectopic expression of miR-141-3p or knockdown of circDLG1 markedly reduced the protein level of CXCL12, inhibition of miR-141-3p significantly increased the protein level of CXCL12, and the increased level of CXCL12 by miR-141-3p inhibition could be partly restored by circDLG1 knockdown (Fig. 5j). Accordingly, the luciferase activity assay showed that ectopic expression of miR-141-3p or knockdown of circDLG1 significantly reduced the luciferase activity of wild-type CXCL12-3′-UTR, inhibition of miR-141-3p significantly increased the luciferase activity of wild-type CXCL12-3′-UTR, and the increased luciferase activity of wild-type CXCL12-3′-UTR by miR-141-3p inhibition could be partly restored by circDLG1 knockdown; however, no significant change in luciferase activity was observed when the CXCL12-3′-UTR was mutated (Fig. 5k). In addition, miR-141-3p expression was inversely associated with CXCL12 expression in gastric cancer tissues (Fig. 5l). These results indicated that circDLG1 could function as a miR-141-3p sponge to regulate CXCL12 expression.
CXCL12 is associated with tumor progression and anti-PD-1 resistance in gastric cancer
It has been reported that CXCL12 plays an important role in cancer progression and immune evasion [35, 36], and we questioned whether CXCL12 is implicated in gastric cancer progression and anti-PD-1 therapy resistance. The potential role of CXCL12 was further explored in gastric cancer cells. In vitro experiments showed that knockdown of CXCL12 significantly inhibited the invasion ability of SGC7901 and HGC27 cells, whereas ectopic expression of CXCL12 significantly increased the invasion ability of AGS cells (Additional file 11: Fig. S6a). Moreover, knockdown of CXCL12 markedly decreased the expression of N-cadherin but increased the expression of E-cadherin (Additional file 11: Fig. S6b).
To evaluate the effect of CXCL12 on anti-PD-1 efficacy, MFC-sh-NC or MFC-sh-CXCL12 cells were injected into C57BL/6 mice subcutaneously or via the tail vein. Then, the mice were treated with 0.9% normal saline (NS), AMD3100 (CXCR4 inhibitor), or anti-PD-1. The results showed that compared with those in the control group, knockdown of CXCL12 significantly inhibited tumor growth and metastasis, and anti-PD-1 treatment further reduced tumor growth and metastasis in C57BL/6 mice (Fig. 6a and b). Likewise, blocking CXCL12/CXCR4 using AMD3100 significantly inhibited tumor growth and metastasis, and the combination of AMD3100 and anti-PD-1 further reduced tumor growth and metastasis (Fig. 6a and b). Moreover, Kaplan–Meier analysis indicated that knockdown of CXCL12 or AMD3100 treatment could significantly prolong the OS of mice, and anti-PD-1 therapy could further increase the OS in mice (Fig. 6c). In addition, flow cytometric analysis of immune cells in subcutaneous tumors showed that knockdown of CXCL12, AMD3100 treatment and anti-PD-1 therapy significantly increased the number of CD8+ T cells and IFNγ+ cells but decreased the number of MDSCs (Fig. 6d-f). We then evaluated the association of CXCL12 expression and the anti-PD-1 response in gastric cancer patients. The results showed that high expression of CXCL12 had a strong trend to be associated with resistance to anti-PD-1 therapy, and patients with high CXCL12 expression presented with adverse survival after anti-PD-1 therapy (Fig. 6g-i).
Increasing evidence has indicated that circRNAs are involved in different diseases [37, 38]. Moreover, many circRNAs are expressed in a cell type-specific or tissue-specific manner , indicating that they might play important biological functions. The dysregulation of circRNAs has been found in several pathological processes, including cardiac hypertrophy, neurological disorders and tumor development [38, 40]. By analyzing the circRNA expression profiles in primary and paired distant metastatic gastric cancer tissues as well as tissues that respond to or are resistant to anti-PD-1 therapy, we found that circDLG1 (circRNA ID: hsa_circ_0008583) was upregulated in distant metastatic lesions and tissues resistant to anti-PD-1 therapy, and high circDLG1 expression was associated with worse PFS in gastric cancer patients treated with anti-PD-1 therapy. This finding indicated that circDLG1 might play an important role in gastric cancer progression and anti-PD-1 resistance.
Further biological experiments demonstrated that circDLG1 could promote cell proliferation, migration, EMT, and stem cell formation. Moreover, an in vivo assay showed that circDLG1 induced gastric cancer progression and metastasis in immunocompetent mice but not in immunodeficient mice. This finding implies that circDLG1 might be involved in the regulation of immunity. Indeed, further analysis showed that knockdown of circDLG1 increased CD8+ T cells while decreasing MDSCs. Previous reports have indicated that EMT characteristics are associated with immune evasion and anti-PD-1 resistance in different tumors [41,42,43]. Our results showed that circDLG1 is critical for the two hallmark processes of tumors, EMT and immune evasion.
Previous studies have shown that several biomarkers, including MSI, PD-L1 expression, TMB, EBV status and Pold/Pole, might be associated with the response to anti-PD-1 therapy [9, 44]. Our previous study showed that the circulating tumor DNA landscape could predict the efficacy of immune checkpoint inhibitors in some gastric cancer patients . However, more biomarkers are needed to identify gastric cancer patients who can benefit from anti-PD-1 therapy. In this study, we found that circDLG1 upregulated the expression of CXCL12. Indeed, previous studies have indicated that CXCL12 is involved in the progression and metastasis of gastric cancer . Moreover, increasing evidence indicates that CXCL12 can induce immune evasion by recruiting infiltrating MDSCs into the tumor microenvironment . In this study, we showed that knockdown of CXCL12 could inhibit the invasion ability and improve the sensitivity of gastric cancer to anti-PD-1 therapy. CXCL12 expression is associated with the PFS of patients who receive anti-PD-1 therapy. This finding is in accordance with previous studies of other tumor types [34, 47]. Taken together, these data demonstrated that CXCL12 might be used as a biomarker to predict gastric cancer resistance to anti-PD-1 therapy. Furthermore, the CXCL12/CXCR4 axis has been found to play an important role in anti-PD-1 resistance in other cancer types . Blocking CXCL12/CXCR4 using AMD3100 (a CXCR4 inhibitor) could improve the efficacy of anti-PD-1 in melanoma . Interestingly, we found that AMD3100 treatment could significantly improve the efficacy of anti-PD-1 therapy in gastric cancer in vivo. This finding indicates that targeting CXCL12 might be a promising strategy to improve anti-PD-1 efficacy in gastric cancer.
Increasing evidence indicates that circRNAs can act as miRNA sponges to regulate the expression of protein-coding genes. For example, the circRNA ACVR2A suppresses bladder cancer cell proliferation and metastasis through the miR-626/EYA4 axis ; the circRNA circ-RanGAP1 regulates VEGFA expression by targeting miR-877-3p to facilitate gastric cancer invasion and metastasis , and the circRNA circCCDC9 acts as a sponge for miR-6792-3p to inhibit gastric cancer progression . In the present study, we found that circDLG1 was mainly located in the cytoplasm of gastric cancer cells, and RIP and RNA pulldown assays showed that circDLG1 could interact with miR-141-3p. Luciferase activity assays confirmed the direct relationship between miR-141-3p and circDLG1. In addition, RNA FISH assays indicated that miR-141-3p and circDLG1 could colocalize in the cytoplasm. These data demonstrate that circDLG1 could sponge miR-141-3p to regulate the expression of CXCL12. Our study has several limitations. First, the sample size was relatively small, and further validation is needed using a larger cohort of patients. Second, we did not detect circDLG1 in the plasma, which might be due to the degradation of RNA after a long storage time of the samples. Further validation is needed to confirm the existence of circDLG1 in plasma or exosomes. Third, the mouse model could not fully reflect the tumor immune microenvironment of the patients, and further study using a patient-derived xenograft model is needed to confirm our results.
In conclusion, we identified a novel circRNA, circDLG1, that is upregulated in gastric cancer distant metastatic lesions and tissues resistant to anti-PD-1 therapy. Further functional studies showed that circDLG1 could promote gastric cancer proliferation, invasion, immune evasion and anti-PD-1 therapy resistance. Clinically, circDLG1 expression was associated with adverse prognosis in gastric cancer patients treated with anti-PD-1 therapy. Additional experiments indicated that circDLG1 acts as a miR-141-3p sponge to regulate the expression of CXCL12, which could induce MDSC infiltration to impair the function of CD8+ T cells. Our study revealed the significance of the circDLG1/CXCL12 axis in the immune evasion and anti-PD-1 resistance of gastric cancer.
Availability of data and materials
All the data are included in this published article and additional files.
Discs large homolog 1
Fluorescence in situ hybridization
Myeloid-derived suppressor cells
Quantitative real-time polymerase chain reaction
Programmed death ligand 1
In situ hybridization
Cell Counting Kit-8
High-power visual field
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.
Lu C, Paschall AV, Shi H, Savage N, Waller JL, Sabbatini ME, et al. The MLL1-H3K4me3 axis-mediated PD-L1 expression and pancreatic cancer immune evasion. J Natl Cancer Inst. 2017;109(6):djw283.
Vinay DS, Ryan EP, Pawelec G, Talib WH, Stagg J, Elkord E, et al. Immune evasion in cancer: mechanistic basis and therapeutic strategies. Semin Cancer Biol. 2015;35(Suppl):S185–98.
Spranger S. Mechanisms of tumor escape in the context of the T-cell-inflamed and the non-T-cell-inflamed tumor microenvironment. Int Immunol. 2016;28(8):383–91.
Wu AA, Drake V, Huang HS, Chiu S, Zheng L. Reprogramming the tumor microenvironment: tumor-induced immunosuppressive factors paralyze T cells. Oncoimmunology. 2015;4(7):e1016700.
Rizvi NA, Mazieres J, Planchard D, Stinchcombe TE, Dy GK, Antonia SJ, et al. Activity and safety of nivolumab, an anti-PD-1 immune checkpoint inhibitor, for patients with advanced, refractory squamous non-small-cell lung cancer (CheckMate 063): a phase 2, single-arm trial. Lancet Oncol. 2015;16(3):257–65.
Ansell SM, Lesokhin AM, Borrello I, Halwani A, Scott EC, Gutierrez M, et al. PD-1 blockade with nivolumab in relapsed or refractory Hodgkin’s lymphoma. N Engl J Med. 2015;372(4):311–9.
Wang F, Wei XL, Wang FH, Xu N, Shen L, Dai GH, et al. Safety, efficacy and tumor mutational burden as a biomarker of overall survival benefit in chemo-refractory gastric cancer treated with toripalimab, a PD-1 antibody in phase Ib/II clinical trial NCT02915432. Ann Oncol. 2019;30(9):1479–86.
Kim ST, Cristescu R, Bass AJ, Kim KM, Odegaard JI, Kim K, et al. Comprehensive molecular characterization of clinical responses to PD-1 inhibition in metastatic gastric cancer. Nat Med. 2018;24(9):1449–58.
Davis AA, Patel VG. The role of PD-L1 expression as a predictive biomarker: an analysis of all US Food and Drug Administration (FDA) approvals of immune checkpoint inhibitors. J Immunother Cancer. 2019;7(1):278.
Shitara K, Ozguroglu M, Bang YJ, Di Bartolomeo M, Mandala M, Ryu MH, et al. Pembrolizumab versus paclitaxel for previously treated, advanced gastric or gastro-oesophageal junction cancer (KEYNOTE-061): a randomised, open-label, controlled, phase 3 trial. Lancet. 2018;392(10142):123–33.
Shitara K, Van Cutsem E, Bang YJ, Fuchs C, Wyrwicz L, Lee KW, et al. Efficacy and safety of pembrolizumab or pembrolizumab plus chemotherapy vs chemotherapy alone for patients with first-line, advanced gastric Cancer: the KEYNOTE-062 phase 3 randomized clinical trial. JAMA Oncol. 2020;6(10):1571–80.
Janjigian YY, Bendell J, Calvo E, Kim JW, Ascierto PA, Sharma P, et al. CheckMate-032 study: efficacy and safety of nivolumab and nivolumab plus ipilimumab in patients with metastatic esophagogastric cancer. J Clin Oncol. 2018;36(28):2836–44.
Bang YJ, Ruiz EY, Van Cutsem E, Lee KW, Wyrwicz L, Schenker M, et al. Phase III, randomised trial of avelumab versus physician’s choice of chemotherapy as third-line treatment of patients with advanced gastric or gastro-oesophageal junction cancer: primary analysis of JAVELIN Gastric 300. Ann Oncol. 2018;29(10):2052–60.
Kang YK, Boku N, Satoh T, Ryu MH, Chao Y, Kato K, et al. Nivolumab in patients with advanced gastric or gastro-oesophageal junction cancer refractory to, or intolerant of, at least two previous chemotherapy regimens (ONO-4538-12, ATTRACTION-2): a randomised, double-blind, placebo-controlled, phase 3 trial. Lancet. 2017;390(10111):2461–71.
Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, et al. Circular RNA: a new star of noncoding RNAs. Cancer Lett. 2015;365(2):141–8.
Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32(5):453–61.
Chen I, Chen CY, Chuang TJ. Biogenesis, identification, and function of exonic circular RNAs. Wiley Interdiscip Rev RNA. 2015;6(5):563–79.
Qu S, Zhong Y, Shang R, Zhang X, Song W, Kjems J, et al. The emerging landscape of circular RNA in life processes. RNA Biol. 2017;14(8):992–9.
Dong R, Ma XK, Chen LL, Yang L. Increased complexity of circRNA expression during species evolution. RNA Biol. 2017;14(8):1064–74.
Ju HQ, Zhao Q, Wang F, Lan P, Wang Z, Zuo ZX, et al. A circRNA signature predicts postoperative recurrence in stage II/III colon cancer. EMBO Mol Med. 2019;11(10):e10168.
Thomson DW, Dinger ME. Endogenous microRNA sponges: evidence and controversy. Nat Rev Genet. 2016;17(5):272–83.
Zhang Y, Zhang XO, Chen T, Xiang JF, Yin QF, Xing YH, et al. Circular intronic long noncoding RNAs. Mol Cell. 2013;51(6):792–806.
Li Z, Huang C, Bao C, Chen L, Lin M, Wang X, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22(3):256–64.
Liu H, Chen D, Bi J, Han J, Yang M, Dong W, et al. Circular RNA circUBXN7 represses cell growth and invasion by sponging miR-1247-3p to enhance B4GALT3 expression in bladder cancer. Aging (Albany NY). 2018;10(10):2606–23.
Rong D, Lu C, Zhang B, Fu K, Zhao S, Tang W, et al. CircPSMC3 suppresses the proliferation and metastasis of gastric cancer by acting as a competitive endogenous RNA through sponging miR-296-5p. Mol Cancer. 2019;18(1):25.
Han D, Li J, Wang H, Su X, Hou J, Gu Y, et al. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017;66(4):1151–64.
Hsiao KY, Lin YC, Gupta SK, Chang N, Yen L, Sun HS, et al. Noncoding effects of circular RNA CCDC66 promote colon cancer growth and metastasis. Cancer Res. 2017;77(9):2339–50.
Lu YX, Ju HQ, Wang F, Chen LZ, Wu QN, Sheng H, et al. Inhibition of the NF-kappaB pathway by nafamostat mesilate suppresses colorectal cancer growth and metastasis. Cancer Lett. 2016;380(1):87–97.
Chen DL, Wang DS, Wu WJ, Zeng ZL, Luo HY, Qiu MZ, et al. Overexpression of paxillin induced by miR-137 suppression promotes tumor progression and metastasis in colorectal cancer. Carcinogenesis. 2013;34(4):803–11.
Taki M, Abiko K, Baba T, Hamanishi J, Yamaguchi K, Murakami R, et al. Snail promotes ovarian cancer progression by recruiting myeloid-derived suppressor cells via CXCR2 ligand upregulation. Nat Commun. 2018;9(1):1685.
Chen DL, Wang ZQ, Zeng ZL, Wu WJ, Zhang DS, Luo HY, et al. Identification of microRNA-214 as a negative regulator of colorectal cancer liver metastasis by way of regulation of fibroblast growth factor receptor 1 expression. Hepatology. 2014;60(2):598–609.
Chen DL, Lu YX, Zhang JX, Wei XL, Wang F, Zeng ZL, et al. Long non-coding RNA UICLM promotes colorectal cancer liver metastasis by acting as a ceRNA for microRNA-215 to regulate ZEB2 expression. Theranostics. 2017;7(19):4836–49.
Zeng Y, Li B, Liang Y, Reeves PM, Qu X, Ran C, et al. Dual blockade of CXCL12-CXCR4 and PD-1-PD-L1 pathways prolongs survival of ovarian tumor-bearing mice by prevention of immunosuppression in the tumor microenvironment. FASEB J. 2019;33(5):6596–608.
Feng W, Huang W, Chen J, Qiao C, Liu D, Ji X, et al. CXCL12-mediated HOXB5 overexpression facilitates colorectal cancer metastasis through transactivating CXCR4 and ITGB3. Theranostics. 2021;11(6):2612–33.
Daniel SK, Seo YD, Pillarisetty VG. The CXCL12-CXCR4/CXCR7 axis as a mechanism of immune resistance in gastrointestinal malignancies. Semin Cancer Biol. 2020;65:176–88.
Zhou R, Wu Y, Wang W, Su W, Liu Y, Wang Y, et al. Circular RNAs (circRNAs) in cancer. Cancer Lett. 2018;425:134–42.
Shang Q, Yang Z, Jia R, Ge S. The novel roles of circRNAs in human cancer. Mol Cancer. 2019;18(1):6.
Chen LL. The biogenesis and emerging roles of circular RNAs. Nat Rev Mol Cell Biol. 2016;17(4):205–11.
Arnaiz E, Sole C, Manterola L, Iparraguirre L, Otaegui D, Lawrie CH. CircRNAs and cancer: biomarkers and master regulators. Semin Cancer Biol. 2019;58:90–9.
Kolijn K, Verhoef EI, Smid M, Bottcher R, Jenster GW, Debets R, et al. Epithelial-mesenchymal transition in human prostate cancer demonstrates enhanced immune evasion marked by IDO1 expression. Cancer Res. 2018;78(16):4671–9.
Hu B, Tian X, Li Y, Liu Y, Yang T, Han Z, et al. Epithelial-mesenchymal transition may be involved in the immune evasion of circulating gastric tumor cells via downregulation of ULBP1. Cancer Med. 2020;9(8):2686–97.
Hong W, Xue M, Jiang J, Zhang Y, Gao X. Circular RNA circ-CPA4/ let-7 miRNA/PD-L1 axis regulates cell growth, stemness, drug resistance and immune evasion in non-small cell lung cancer (NSCLC). J Exp Clin Cancer Res. 2020;39(1):149.
De Rosa S, Sahnane N, Tibiletti MG, Magnoli F, Vanoli A, Sessa F, et al. EBV+ and MSI gastric cancers harbor high PD-L1/PD-1 expression and high CD8+ Intratumoral lymphocytes. Cancers (Basel). 2018;10(4):102.
Jin Y, Chen DL, Wang F, Yang CP, Chen XX, You JQ, et al. The predicting role of circulating tumor DNA landscape in gastric cancer patients treated with immune checkpoint inhibitors. Mol Cancer. 2020;19(1):154.
Izumi D, Ishimoto T, Miyake K, Sugihara H, Eto K, Sawayama H, et al. CXCL12/CXCR4 activation by cancer-associated fibroblasts promotes integrin beta1 clustering and invasiveness in gastric cancer. Int J Cancer. 2016;138(5):1207–19.
Wei CY, Zhu MX, Lu NH, Liu JQ, Yang YW, Zhang Y, et al. Circular RNA circ_0020710 drives tumor progression and immune evasion by regulating the miR-370-3p/CXCL12 axis in melanoma. Mol Cancer. 2020;19(1):84.
Dong W, Bi J, Liu H, Yan D, He Q, Zhou Q, et al. Circular RNA ACVR2A suppresses bladder cancer cells proliferation and metastasis through miR-626/EYA4 axis. Mol Cancer. 2019;18(1):95.
Lu J, Wang YH, Yoon C, Huang XY, Xu Y, Xie JW, et al. Circular RNA circ-RanGAP1 regulates VEGFA expression by targeting miR-877-3p to facilitate gastric cancer invasion and metastasis. Cancer Lett. 2020;471:38–48.
Luo Z, Rong Z, Zhang J, Zhu Z, Yu Z, Li T, et al. Circular RNA circCCDC9 acts as a miR-6792-3p sponge to suppress the progression of gastric cancer through regulating CAV1 expression. Mol Cancer. 2020;19(1):86.
We thank all of the patients enrolled in this study.
This study was supported by the National Key R&D Program of China (2018YFC1313300); Science and Technology Program of Guangdong (2019B020227002); CAMS Innovation Fund for Medical Sciences (CIFMS) (2019-I2M-5-036); National Natural Science Foundation of China (No.81372570); Natural Science Foundation of Guangdong Province (No. 2014A030312015; No. 2016A030310195); and The Outstanding Youth Cultivation Fund of SYSUCC (No. PT04141002).
Ethics approval and consent to participate
This study was approved by the ethics committee of the Sun Yat-sen University Cancer Center, and informed consent was obtained from all patients.
Consent for publication
Informed consent was obtained from all patients before the use of these clinical materials for research purposes.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers used in the paper were listed.
The information for antibodies used in flow cytometry.
RT–qPCR analysis of circDLG1 expression in primary gastric cancer tissues and adjacent normal tissues (n = 30, P = 0.14). P = 0.14, Paired Student’s t-test. The data are representative of three technical replicates.
The association between clinicopathologic parameters and circDLG1 expression in 30 gastric cancer patients.
IHC analysis of E-cadherin and vimentin expression in primary gastric cancer tissues with high and low circDLG1 expression. Scale bar, 100 μm.
circDLG1 expression in MFC cells treated with shRNAs. *P < 0.05, one-way ANOVA test for the comparison of the three groups, followed by Student’s t-test for the comparison of sh-1 and sh-2 group with sh-NC group, respectively; multigroup comparisons were adjusted using the Bonferroni method. The experiment was repeated three times.
IHC analysis of CD4, CD8, Gr-1, CD11b expression in MFC-sh-circDLG1-derived tumors and MFC-sh-NC-derived tumors. Scale bar, 100 μm.
CircDLG1 upregulates CXCL12 to promote cell progression and immune evasion in gastric cancer. a, A gene expression profile was performed in gastric cancer cells with or without circDLG1 knockdown. SGC7901 cells were transfected with circDLG1 shRNA, and the expression level of mRNAs was detected by RNA sequencing in SGC7901 cells. b, Gene annotation enrichment analysis showed that differentially expressed genes involved in T cell migration, EMT and immune responses were significantly enriched based on the GO analysis, and TNF signaling and the PI3K-Akt signaling pathway were significantly enriched based on the KEGG analysis. c, qRT–PCR analysis of CXCL12 expression in SGC7901 gastric cancer cells after knockdown of circDLG1. **P < 0.01, Student’s t-test; the experiment was repeated three times. d, qRT–PCR analysis of CXCL12 expression in SGC7901 gastric cancer cells after ectopic expression of circDLG1. **P < 0.01, Student’s t-test; the experiment was repeated three times. e, Western blot analysis of CXCL12 in SGC7901 and HGC27 gastric cancer cells after knockdown of circDLG1. f, Western blot analysis of CXCL12 in SGC7901 cells treated with different vectors. g, The cell migration ability of SGC7901 gastric cancer cells treated with different vectors. Scale bar, 100 μm. h, The association between circDLG1 expression and CXCL12 expression as measured by qRT–PCR in gastric cancer tissues. n = 30, P = 0.002, Pearson correlation analysis.
The potential miRNAs that could bind to has_circ_0008583 (circDLG1).
CXCL12 is associated with an aggressive tumor phenotype in gastric cancer cells. a, The invasion ability of the gastric cancer cell lines SGC7901 and HGC27 after knockdown of CXCL12. *P < 0.05, one-way ANOVA test for the comparison of the three groups, followed by Student’s t-test for the comparison of sh-1 and sh-2 group with sh-NC group, respectively; multigroup comparisons were adjusted using the Bonferroni method. The invasion ability of AGS cells treated with the CXCL12 ectopic expression vector. *P < 0.05, Student’s t-test. Scale bar, 100 μm. b, FISH analysis of N-cadherin and E-cadherin expression after knockdown of CXCL12 in SGC7901 cells. Scale bar, 50 μm.
About this article
Cite this article
Chen, DL., Sheng, H., Zhang, DS. et al. The circular RNA circDLG1 promotes gastric cancer progression and anti-PD-1 resistance through the regulation of CXCL12 by sponging miR-141-3p. Mol Cancer 20, 166 (2021). https://doi.org/10.1186/s12943-021-01475-8
- Gastric cancer
- Immune evasion