ORIGINAL RESEARCH article

Front. Oncol., 23 July 2020

Sec. Cancer Molecular Targets and Therapeutics

Volume 10 - 2020 | https://doi.org/10.3389/fonc.2020.01104

Circular RNA DGKB Promotes the Progression of Neuroblastoma by Targeting miR-873/GLI1 Axis

  • 1. Department of General Surgery, Second Affiliated Hospital of Nanchang University, Nanchang, China

  • 2. Department of General Surgery, Jiangxi Provincial Children's Hospital, Nanchang, China

  • 3. Department of Ultrasound, Jiangxi Provincial Maternal and Child Health Hospital, Nanchang, China

Abstract

Accumulated evidences suggested that circular RNAs (circRNA) played critical roles in tumorigenesis and progression. To our knowledge, no study reported the function of circular RNA DGKB (circDGKB, circRNA ID: hsa_circ_0133622) on progression of neuroblastoma (NB). Here, we showed that circDGKB was upregulated in NB tissues compared to the normal dorsal root ganglia. Moreover, the expression level of circDGKB was negatively correlated with the survival rate of NB patients. Mechanically, overexpression of circDGKB promoted the proliferation, migration, invasion, and tumorigenesis of NB cells and reduced cell apoptosis, and vice versa. In addition, qRT-PCR and/or Western blot results showed that circDGKB overexpression inhibited the expression level of miR-873 and enhanced GLI1 expression. Moreover, miR-873 functioned an opposite role to circDGKB and significantly weakened circDGKB role in promoting NB progression. Furthermore, GLI1 upregulation also rescued the miR-873 role in inhibiting NB progression. In conclusion, our work proved that circDGKB promoted NB progression via targeting miR-873/GLI1 axis in vitro and in vivo. Our study provided a new target for NB treatment and indicated that circDGKB could act as a novel diagnostic marker for NB.

Introduction

Neuroblastoma (NB) is a malignant tumor mainly occurring in children (1), which commonly arises in adrenal glands, and also occurs in the abdomen, neck, spine, or chest (2). The most malignant tumors have amplification of the MYCN oncogene encoding a transcription factor, which is associated with poor survival of NB, even in localized diseases (3). NB is divided into two classifications in clinic, low-risk disease and high-risk disease. Patients with low-risk disease respond positively to therapy and even can recover without any treatment. Unfortunately, 40% of NB patients suffer from high-risk disease which is easy to relapse and lethal (1, 4). Various kinds of therapeutic methods have been adopted to treat high-risk NB, including surgical treatment, chemotherapy, radiotherapy, immunotherapy, cis-retinoic acid, and proton therapy (5), but the prognosis of NB is still poor due to the complex mechanisms underlying NB occurrence and development. Therefore, it is urgent to further reveal the molecular mechanisms of NB in order to find novel and potent targets for the treatment of NB.

Circular RNAs (circRNAs) are a class of non-coding RNA molecules generated from exons or introns by forming covalently-closed loops (6). Previous studies have suggested that circRNAs can serve as the sponges of microRNAs (miRNAs) to inhibit the function of miRNAs. For example, circATP82A2 repressed the expression of miR-443 via binding to miR-443 directly (7). circRNAs are abundantly expressed in eukaryotic organisms and closely implicated in a variety of disease, such as diabetes, cardiovascular diseases (CVDs), and cancers (8, 9). Numerous studies demonstrated that circRNA/miRNA axis played critical roles in the progression of cancers. For instance, circ0103552 significantly promoted cell proliferation and invasion via sponging miR-1236 in breast cancer (10). Mao et al. (11) demonstrated that circ0068871 could accelerate the progression of bladder cancer by targeting miR-181a-5p and then regulate the expression of FGFR3 and activate STAT3. Huang et al. (12) identified that circ100338 was highly expressed in hepatocellular carcinoma (HCC) and functioned as an inhibitor of miR-141-3p to promote the proliferation of HCC cells. However, the functions of lots of circRNAs are waiting for us to explore.

MiRNAs are short non-coding transcripts with 22 nucleotides in length (13). Numbers of studies have demonstrated that miRNAs played significant roles in RNA translation and degradation via binding to the 3′UTR of mRNA (14). MiRNAs were proved to participate in progress of many kinds of cancers by acting as the targets of circRNAs and then regulating the expression of oncogenes or tumor suppressive genes (15). In the present study, we found that circDGKB (circRNA ID: hsa_circ_0133622) was highly expressed in NB. In addition, circDGKB could act as an oncogene in NB by targeting miR-873/GLI1 (Glioma-associated oncogene 1) axis. Our study provided a novel potential target for the diagnosis and treatment of NB in clinical.

Materials and Methods

Tissue Collection

Thirty NB tissues were obtained from 30 NB patients at Second Affiliated Hospital of Nanchang University during 2015 to 2017. Meanwhile, 10 human normal dorsal root ganglia were collected from interrupted pregnancies, which were used as a normal control group. Patients participated in this study never received any other therapeutic method but surgery and the informed consents were signed by their parents. This study was approved by the Research Ethics Committee of Second Affiliated Hospital of Nanchang University.

circRNA Microarray

After being obtained from surgical specimens, 3 NB tissues and 3 normal dorsal root ganglia were immediately frozen using liquid nitrogen. Sample preparation and microarray hybridization were carried out in accordance with the protocols of Arraystar (Rockville, MD, USA). circRNAs were enriched through removing linear RNAs with Rnase R (Epicentre, Madison, WI, USA), followed by amplification and labeling with Arraystar Super RNA Labeling Kit (Arraystar). Then, the Arraystar Human circRNA Array (8 × 15 K) was used for hybridization and scanned by the Agilent Scanner G2505C (Jamul, CA, USA). circRNAs demonstrating fold-changes of ≥2 and P-values of < 0.05 were considered as significantly differentially expressed.

Cell Culture

Two NB cell lines SK-N-SH and SH-SY5Y were bought from American Type Culture Collection (ATCC, VA, USA) and cultured in Dulbecco's Modified Eagle Medium (DMEM, Gibco, CA, USA), with 10% fetal serum (FBS, Gibco, CA, USA) and 1% penicillin/streptomycin (Invitrogen, CA, USA).

Cell Infection

The lentiviral vectors used to upregulate circDGKB (Lentiv-circDGKB), miR-873 (mimic-miR-873), and GLI1 (Lentiv-GLI1) in SK-N-SH and SH-SY5Y, and the lentiviral vectors used to downregulate circDGKB (shRNA-circDGKB), miR-873 (inhibitor-miR-873), and GLI1 (shRNA-GLI1), in addition to the negative control vector (Lentiv-NC, inhibitor-NC, mimic-NC), were all designed and synthesized by GenePharma (Shanghai, China). These vectors were introduced into SK-N-SH and SH-SY5Y cells by cell infection with the help of polybrene (Hanbio, Inc, Shanghai, China), followed by incubation with G418 or puromycin for 14 days to establish the stable transfection cell lines used in in vivo experiments.

Cell Growth Assay

Cell Counting Kit-8 (CCK-8) (Beyotime, Beijing, China) was adopted to assess the proliferation of NB cells. The cells were seeded into 96-well plates at the density of 2,000 cell/well and incubated at 37°C. Following 24 h, the cells were infected with lentiviral vectors instantaneously. Forty-eight hours later, 10 μl CCK-8 solution was added into each well and incubation at 37°C for another 2 h. The absorbance was detected by a plate reader (model 680; Bio-Rad, Hertfordshire, UK) at 450 nm.

Cell Migration Assay

NB cells were seeded into six-well plates at the density of 3,000 cells/group. The wounds were scraped by pipettes following cell infection. The width of the wound was captured at 0 and 24 h following scratch by DM2500 bright field microscope (LEICA, Wetzlar, Germany) and the imageJ software was used to quantify the migration distance.

Transwell Invasion Assay

The upper chamber of the transwell chambers (8 μm, BD Bioscience, USA) were first coated with matrigel. Then, the infected SK-N-SH and SH-SY5Y cells were suspended with FBS-free medium and seeded into the upper chamber at the density of 3,000 cells/group. At the same time, in the bottom chamber were added 600 μl complete medium containing 10% FBS (Gibco). Following 48 h of incubation at 37°C, the inversed cells were fixed with 4% paraformaldehyde and stained with crystal violet. The number of inversed cells from 6 randomly selected fields were counted.

Cell Apoptosis Assay

Cell apoptosis detection was performed by using flow cytometry with the help of Annexin V-FITC/PI Apoptosis Detection Kit (Vazyme, Biotech, Jiangsu, China) according to the instructions of the manufacturer.

Luciferase Reporter Assay

The wild type (circDGKB-WT) and mutated type (circDGKB-MUT) vectors between the binding sites in circDGKB were synthesized by GenePharma. Then, SK-N-SH cells were co-transfected with mimic-NC/mimic-miR-873 and circDGKB-WT/circDGKB-MUT using Lipofectamine (Invitrogen) or polybrene. Following 48 h, the relative luciferase activity was detected by Dual Luciferase Assay System (Promega, Madison, WI, USA) according to the manufacturer's instructions.

Western Blot

The cells were lysed with RIPA lysis buffer (Thermo Fisher Scientific, NYC, USA) supplemented with protease inhibitor cocktail (APExBIO, AL, and USA) on ice. Then, the supernatant was transferred into new microphage tubes (Eppendorf, Hamburg, Germany), followed by centrifugation at 12,000 rpm for 30 min at 4°C. The concentrations of protein samples were measured with BCA Protein Quantification Kit (YEASEN, Shanghai, China). Next, 30 μg protein of each sample was used for electrophoresis. After transferring the protein to PDVF membrane, the membranes were blocked with 5% milk in TBST for 1 h. Primary antibodies included anti-ZBE1 antibody (No. ab203829, Abcam, Cambridge, MA, USA), anti-GLI1 (No. ab49314, Abcam), and anti-GAPDH antibody (No. ab011-040, Multi sciences, Zhejiang, China). The secondary antibodies was used as follows: anti-Mouse IgG-HRP (No. RA230188, Thermo Fisher Scientific).

Reverse Real-Time Quantitative Polymerase Chain Reaction (qRT-PCR)

Total RNA extraction from cells and tissues was carried out with Trizol reagent (Invitrogen). Then, the RNAs were reversed into cDNAs by using the high capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA). TaqMan Universal Master MixII (Takara, Dalian, China) was adopted to perform real-time quantitative PCR and the internal reference was GAPDH/U6.

Primers are shown as follows:

  • circDGKB-forward: 5′-AGACTCTGCCCACTTCAGGA-3′,

  • circDGKB-reverse: 5′-AGGCACTGGGTCTCCTTTCT-3′;

  • miR-873-forward: 5′-TGTGCATTTGCAGGAACTTGT-3′,

  • miR-873-reverse: 5′-GGGAACTCATCAGTCTCCTGTTC-3′;

  • GAPDH-forward: 5′-CCACCCCCAATGTCTCTGTT-3′,

  • GAPDH-reverse: 5′-ATGGATGAACGGCAATCCCC-3′.

In vivo Tumorigenicity Assay

Animal study was performed in accordance with the Research Ethics Committee of Second Affiliated Hospital of Nanchang University. In brief, 4-week nude mice were purchased from Shanghai Slac Laboratory Animal Company (Shanghai, China) and the flanks of mice were injected with 1 × 107 stably infected SK-N-SH cells, including shRNA-circDGKB, Lentiv-circDGKB, shRNA-NC, Lentiv-NC, mimic-NC, mimic-miR-873, inhibitor-miR-873, inhibitor-NC, Lentiv-GLI1, shRNA-GLI1. After 4 weeks, mice were euthanized and tumor weights were measured.

Statistical Analysis

Date from three independent experiments are shown as mean ± standard deviation (SD). Comparison among 2 groups or > 2 groups was executed by using the Student's t-test or one-way analysis of variance (ANOVA) using SPSS 17.0 software. P < 0.05 was considered as statistically significant. ROC (receiver operating characteristic) curve was established to evaluate clinical diagnostic value of circDGKB in NB, with larger area under the curve meaning the higher prognostic. The relationship between circDGKB expression levels and the overall survival rates of NB patients was assessed by using Kaplan-Meier analysis with log rank test.

Result

circDGKB Was Highly Expressed in NB Tissue Specimens

To find the specific circRNAs which are differentially expressed in NB specimens, 3 NB tissues obtained from NB patients and 3 normal dorsal root ganglia collected from interrupted pregnancies were subjected to circRNA microarray assay. A total of 235 circRNAs were found to be differentially expressed in NB specimens compared with the dorsal root ganglia, among which 138 circRNAs were upregulated and 97 circRNAs were downregulated (Figures 1A,B). And circDGKB expression was increased, highest among the 138 upregulated circRNAs. The qRT-PCR assay was carried out to confirm the expression pattern of circDGKB in 30 paired NB tissues and 10 normal dorsal root ganglia. The results showed that the expression level of circDGKB in NB tissues was significantly increased (Figure 1C). Figure 1D demonstrated the relative expression levels of circDGKB in 30 NB tissues to normal dorsal root ganglia. Taken together, these results demonstrated that circDGKB was overexpressed in NB tissues.

Figure 1

circDGKB Accelerated the Progression of NB

To investigate the role of circDGKB in the progression of NB, gain of function and loss of function experiments were conducted in NB cell lines SK-N-SH and SH-SY5Y. QRT-PCR assay was taken to test the infection efficiencies of shRNA-circDGKB and Lentiv-circDGKB. Compared with the control group, circDGKB expression was obviously reduced after the cells were infected with shRNA-circDGKB in both SK-N-SH and SH-SY5Y cell lines (Figures 2A,B). CCK-8 assay demonstrated that downregulation of circDGKB prominently attenuated cell proliferation in SK-N-SH and SH-SY5Y cells (Figures 2C,D). Flow cytometry showed that circDGKB downregulation increased cell apoptosis (Figures 2E,F). In addition, we found that circDGKB downregulation could lead to S phase reduction (Figures 2G,H). Moreover, the wound healing assay and transwell invasion assay were performed to examine the migration and invasion of SK-N-SH and SH-SY5Y cells. As shown in Figures 2I–L, downregulation of circDGKB significantly impaired cell invasion and migration abilities in SK-N-SH and SH-SY5Y cells. Furthermore, the in vivo assay was carried out to explore the function of circDGKB in tumor formation. The result showed that downregulation of circDGKB reduced tumor formation in SK-N-SH cells (Figure 2M).

Figure 2

To the contrary of shRNA-circDGKB infection, circDGKB expression was markedly increased when the cells were infected with Lentiv-circDGKB in both SK-N-SH and SH-SY5Y cell lines (Figures 3A,B). CCK-8 assay demonstrated that overexpression of circDGKB obviously enhanced the proliferative capacity of SK-N-SH and SH-SY5Y cells (Figures 3C,D). Flow cytometry showed that circDGKB overexpression significantly inhibited cell apoptosis (Figures 3E,F). In addition, we found that overexpression of circDGKB could induce cell S phase arrest (Figures 3G,H). Moreover, the wound healing assay and transwell invasion assay showed that overexpression of circDGKB significantly promoted cell invasion and migration in SK-N-SH and SH-SY5Y cells (Figures 3I–L). Furthermore, the in vivo assay showed that overexpression of circDGKB obviously promoted the tumorigenesis of SK-N-SH cells (Figure 3M). These results showed that circDGKB could accelerate the progression of NB.

Figure 3

miR-873/GL1-3 Axis Was a Target of circDGKB in NB Cells

Previous studies proved that circRNA could act as a sponge of miRNA. To explore the underlying targets of circDGKB, we performed miRNA prediction with miRanda. As shown in Figure 4A, results of the luciferase reporter assay showed that miR-873 overexpression dramatically reduced the luciferase activity of circDGKB-WT, whereas the mutation of binding sites abolished this effect. The further bioinformatics analysis predicted that GLI1 was a target for miR-873. To verify the foreclosed results, we conducted “gain” and “lost” experiment in SK-N-SH cells. The results of western blot indicated that circDGKB increased the expression of GLI1 and ZEB1 under the regulation of GLI1 (16) (Figure 4B), while circDGKB downregulation reduced GLI1 and ZEB1 expression (Figure 4C). In addition, miR-873 overexpression downregulated GLI1 and ZEB1 expression (Figure 4D), while knockdown of miR-837 upregulated the expression of ZEB1 and GLI1 (Figure 4E). These results suggested that circDGKB increased the expression of ZEB1 and GLI1 via targeting miR-873.

Figure 4

miR-873 Played a Vital Role in the Progression of NB

To explore the role of miR-873 in the progression of NB, the gain of function and loss of function experiments were carried out in SK-N-SH cells. MiR-873 overexpression significantly reduced the proliferation of SK-N-SH cells, while downregulation led to the opposite results (Figure 5A). Cell apoptosis and S phase arrest were induced by miR-873 upregulation in SK-N-SH cells, and knockdown of miR-873 caused the opposite effects (Figures 5B,C). The results of transwell experiments and wound healing assay showed that overexpression of miR-873 promoted the invasion and migration of SK-N-SH cells, and vice versa (Figures 5D,E). In addition, the in vivo assay demonstrated that overexpression of miR-873 inhibited the tumor formation ability of SK-N-SH cells, and vice versa (Figure 5F). These results suggested that miR-873 inhibited the progression of NB.

Figure 5

circDGKB Promoted NB Progression via Targeting miR-873

Then, we explored whether miR-873 was involved in circDGKB-mediated NB progression. The results showed that miR-873 overexpression significantly weakened circDGKB roles in promoting cell proliferation (Figure 6A), S phase arrest (Figure 6C), invasion (Figure 6D), and migration (Figure 6E), as well as inhibiting cell apoptosis (Figure 6B). These results demonstrated that circDGKB promoted NB progression via targeting miR-873.

Figure 6

miR-873 Promoted NB Progression via Targeting GLI1

Next, we explored GLI1 role in the progression of NB in SK-N-SH cells. Western blot assay was used to verify the efficiency of lentiv-GLI1 and shRNA-GLI1 (Figures 7A,B). As expected, overexpression of GLI1 obviously promoted the proliferation of SK-N-SH cells, while downregulation of GLI1 caused the opposite effect (Figure 7C). GLI1 upregulation significantly inhibited apoptosis and induced S phase arrest in SK-N-SH cells, and vice versa (Figures 7D,E), as well as promoted cell invasion, migration, and tumorigenesis in vivo, and vice versa (Figures 7F–H). These results indicated that GLI1 accelerated the progression of NB.

Figure 7

To explore whether GLI1 was implicated in miR-873-mediated NB progression inhibition, the rescued assays were carried out. Compared with the mimic-miR-873+lentive-NC group, cell proliferation (Figure 8A), invasion (Figure 8D), and migration (Figure 8E) abilities were significantly enhanced, while cell apoptosis were reduced (Figure 8B) and cell numbers in S phase were increased (Figure 8C) in mimic-miR-873+lentiv-GLI1 group. These findings suggested that miR-873 promoted NB progression via targeting GLI1.

Figure 8

circDGKB Acted as a Potential Marker for NB Diagnosis

To investigate whether circDGKB could act as a diagnostic marker for NB, we detected the expression of circDGKB in the blood samples of NB patients and healthy subjects. The results of qRT-PCR showed that circDGKB expression in blood was prominently upregulated in NB patients compared with normal subjects (Figure 9A). The results of ROC analyses demonstrated that the level of circDGKB had an important clinical value in the diagnosis of NB with an AUC (area under the curve) of 0.7778 (Figure 9B). The results of Kaplan-Meier analysis showed that patients with high level of circDGKB had lower survival rates (Figure 9C).

Figure 9

Discussion

circRNA produced from exons or introns is stable existing and was first found in RNA viruses in the 1970s (17, 18), at which time the circRNA was thought to have no biological function (1923). To date, a great deal of studies proved that circRNAs were ubiquitous in mammalian cells and played critical roles in regulating gene expression by targeting miRNA or other molecules including linear slicing and proteins (2426). Growing evidences have proved that circRNAs play a key role in many diseases, including systemic lupus erythematous (27), atherosclerosis (28), Neuropsychiatric Disorders (29), Alzheimer's disease (30), cardiovascular diseases (31), and cancers (32). Qiu et al. (33) identified that the expression level of circ000911 could promote the proliferation and invasion of non-small-cell lung cancer (NSCLC) cells by interacting with miR-22-3p directly and then enhanced the expression level of galectin-1 (Gal-1), p-AKT, and p-ERK1/2. Fang et al. (34) indicated that circ100290 was highly expressed in colorectal cancer (CRC) tissues and cell lines; silence of circ100290 inhibited the proliferation, migration, and invasion of CRC cells, but promoted cell apoptosis. In addition, Fang et al. (34) also proved that circ100290 expression level showed positive correlation with CRC metastasis and negative correlation with patients' prognosis. Gao et al. (35) illustrated that circ0006528 promoted the proliferation, migration, and invasion of breast cancer cells by targeting miR-7-5p/raf1/MAPK/ERK pathways. In the present study, we found that circDGKB expression was prominently unregulated in NB tissues; overexpression of circDGKB promoted the proliferation, migration, and invasion of SK-N-SH cells, and inhibited cell apoptosis and induced S phase arrest, while circDGKB downregulation showed the reverse effects.

MiRNAs, a class of non-coding RNA with ~22 nucleotides, have been proved to play a vital role in tumorigenesis by regulating the expression of oncogenes or tumor suppressive genes. For instance, miR-776 was downregulated in colon cancer, and low expression level of miR-776 was associated with low survival rate of patients. In addition, the authors proved that overexpression of miR-776 obviously inhibited cell proliferation and initiated cell apoptosis in colon cancer cells by targeting MDM4/p53 pathway (36). miR-20a-5p was identified to be highly expressed in triple-negative breast cancer (TNBC) tissues, which promoted TNBC progression by targeting Runt-related transcription factor 3 (RUNX3) (37). In our work, we proved that miR-873 was a direct target of circDGKB, and inhibition of the expression of miR-873 enhanced the proliferation, migration, invasion, and S phase arrest of SK-N-SH cells and reduced cell apoptosis. This results were consistent with previous studies, which indicated the suppressive role of miR-873 in cancers including CRC (38, 39), breast cancer (40), NSCLC (41), and gastric cancer (42). Furthermore, we found that miR-873 overexpression significantly rescued circDGKB role in promoting NB progression, illustrating that circDGKB promoted NB progression by targeting miR-873.

It's well-known that Hedgehog (HH) signaling pathway plays critical roles in amount of cancers, including NB (4346). GLI1, a terminal effectors of the HH factors, was reported to serve as an oncogene in wide ranges of cancers, such as breast cancer (47), CRC (48), prostate cancer (49), glioma (50), pancreatic cancer (51), gastric cancer (52), and cervical cancer (53). Noticeably, GLI1 was reported to be a direct target of miR-873-5p in gastric cancer (42). Similarly, in our work, GLI1 was also demonstrated as a direct target of miR-873 in NB cells, and GLI1 overexpression could markedly weaken miR-873 roles in inhibiting cell proliferation, migration, and invasion of SK-N-SH cells and inducing cell apoptosis and S phase arrest, suggesting that miR-873 inhibited NB progression via targeting GLI1.

Collectively, our work demonstrated that circDGKB acted as an oncogene in NB progression by targeting miR-873/GLI1 axis. In addition, we found the high expression level of circDGKB was associated with the poor prognosis of NB patients, suggesting that circDGKB could be used as a diagnostic marker for NB. Our study suggested that circDGKB was a potential therapeutic target and diagnostic marker for NB.

Statements

Data availability statement

The datasets generated for this study can be found in the National Center for Biotechnology Information. Circular RNA DGKB promotes the progression of neuroblastoma by targeting miR-873-GLI1 axis:PRJNA554935; Human sample from Homo sapiens:SAMN12287376, SAMN12287375, SAMN12287374, SAMN12287373, SAMN12287372, and SAMN12287371; circRNA seq of neuroblastma:SRR9694951, SRR9694952, SRR9694948, SRR9694947, SRR9694949, and SRR9694950.

Ethics statement

The studies involving human participants were reviewed and approved by the Research Ethics Committee of Second Affiliated Hospital of Nanchang University. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin. The animal study was reviewed and approved by the Research Ethics Committee of Second Affiliated Hospital of Nanchang University.

Author contributions

SZ and XY came up with the conception and design of study. JYang and LY performed the experiment and acquired the data. JYan, YX, and WL analyzed the data. JX, JL, and DX performed some operations such as specimen acquisition and disposal. JYang drafted the manuscript. SZ and XY revised the manuscript critically for important intellectual content. All authors agreed to publish the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 81960101 and 81360338).

Acknowledgments

All authors are very grateful to those persons who have guided us.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.01104/full#supplementary-material

References

  • 1.

    SzemesMGreenhoughAMalikK. Wnt Signaling is a major determinant of neuroblastoma cell lineages. Front Mol Neurosci. (2019) 12:90. 10.3389/fnmol.2019.00090

  • 2.

    NagaichU. Neuroblastoma cancer in children. J Adv Pharm Technol Res. (2019) 10:45. 10.4103/japtr.JAPTR_37_19

  • 3.

    MatthayKKMarisJMSchleiermacherGNakagawaraAMackallCLDillerLet al. Neuroblastoma. Nat Rev Dis Primers. (2016) 2:16078. 10.1038/nrdp.2016.78

  • 4.

    Vivancos StalinLGualandiMSchulteJHRenellaRShakhovaOMuhlethaler-MottetA. Expression of the neuroblastoma-associated ALK-F1174L activating mutation during embryogenesis impairs the differentiation of neural crest progenitors in sympathetic ganglia. Front Oncol. (2019) 9:275. 10.3389/fonc.2019.00275

  • 5.

    BagleyAFGrosshansDRPhilipNVFosterJMcAleerMFMcGovernSLet al. Efficacy of proton therapy in children with high-risk and locally recurrent neuroblastoma. Pediatr Blood Cancer. (2019) 66:e27786. 10.1002/pbc.27786

  • 6.

    LuWY. Roles of the circular RNA circ-Foxo3 in breast cancer progression. Cell Cycle. (2017) 16:58990. 10.1080/15384101.2017.1278935

  • 7.

    DingLZhangH. Circ-ATP8A2 promotes cell proliferation and invasion as a ceRNA to target EGFR by sponging miR-433 in cervical cancer. Gene. (2019) 705:1038. 10.1016/j.gene.2019.04.068

  • 8.

    BurdCEJeckWRLiuYSanoffHKWangZSharplessNE. Expression of linear and novel circular forms of an INK4/ARF-associated non-coding RNA correlates with atherosclerosis risk. PLoS Genet. (2010) 6:e1001233. 10.1371/journal.pgen.1001233

  • 9.

    WangYMoYGongZYangXYangMZhangSet al. Circular RNAs in human cancer. Mol Cancer. (2017) 16:25. 10.1186/s12943-017-0598-7

  • 10.

    YangLSongCChenYJingGSunJ. Circular RNA circ_0103552 forecasts dismal prognosis and promotes breast cancer cell proliferation and invasion by sponging miR-1236. J Cell Biochem. (2019) 120:1555360. 10.1002/jcb.28822

  • 11.

    MaoWHuangXWangLZhangZLiuMLiYet al. Circular RNA hsa_circ_0068871 regulates FGFR3 expression and activates STAT3 by targeting miR-181a-5p to promote bladder cancer progression. J Exp Clin Cancer Res. (2019) 38:169. 10.1186/s13046-019-1136-9

  • 12.

    HuangXYHuangZLXuYHZhengQChenZSongWet al. Comprehensive circular RNA profiling reveals the regulatory role of the circRNA-100338/miR-141–3p pathway in hepatitis B-related hepatocellular carcinoma. Sci Rep. (2017) 7:5428. 10.1038/s41598-017-05432-8

  • 13.

    de FoucherTSbeihMUzanJBendifallahSLefevreMChabbert-BuffetNet al. Identification of micro-RNA expression profile related to recurrence in women with ESMO low-risk endometrial cancer. J Transl Med. (2018) 16:131. 10.1186/s12967-018-1515-6

  • 14.

    GuoHIngoliaNTWeissmanJSBartelDP. Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature. (2010) 466:83540. 10.1038/nature09267

  • 15.

    LiuHLeiCHeQPanZXiaoDTaoY. Nuclear functions of mammalian MicroRNAs in gene regulation, immunity and cancer. Mol Cancer. (2018) 17:64. 10.1186/s12943-018-0765-5

  • 16.

    LiJHeYCaoYYuYChenXGaoXet al. Upregulation of twist is involved in Gli1 induced migration and invasion of hepatocarcinoma cells. Biol Chem. (2018) 399:9119. 10.1515/hsz-2018-0131

  • 17.

    WiluszJESharpPA. Molecular biology. A circuitous route to noncoding RNA. Science. (2013) 340:4401. 10.1126/science.1238522

  • 18.

    WuQWangYCaoMPantaleoVBurgyanJLiWXet al. Homology-independent discovery of replicating pathogenic circular RNAs by deep sequencing and a new computational algorithm. Proc Natl Acad Sci USA. (2012) 109:393843. 10.1073/pnas.1117815109

  • 19.

    CocquerelleCMascrezBHetuinDBailleulB. Mis-splicing yields circular RNA molecules. FASEB J. (1993) 7:15560. 10.1096/fasebj.7.1.7678559

  • 20.

    LiPFChenSCChenHLMoXYLiTWShaoYFet al. Using circular RNA as a novel type of biomarker in the screening of gastric cancer. Clin Chim Acta. (2015) 444:1326. 10.1016/j.cca.2015.02.018

  • 21.

    HentzeMWPreissT. Circular RNAs: splicing's enigma variations. EMBO J. (2013) 32:9235. 10.1038/emboj.2013.53

  • 22.

    JeckWRSorrentinoJAWangKSlevinMKBurdCELiuJet al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. (2013) 19:14157. 10.1261/rna.035667.112

  • 23.

    SuzukiHZuoYHWangJHZhangMQMalhotraAMayedaA. Characterization of RNase R-digested cellular RNA source that consists of lariat and circular RNAs from pre-mRNA splicing. Nucleic Acids Res. (2006) 34:e63. 10.1093/nar/gkl151

  • 24.

    HansenTBJensenTIClausenBHBramsenJBFinsenBDamgaardCKet al. Natural RNA circles function as efficient microRNA sponges. Nature. (2013) 495:3848. 10.1038/nature11993

  • 25.

    GuoJUAgarwalVGuoHBartelDP. Expanded identification and characterization of mammalian circular RNAs. Genome Biol. (2014) 15:409. 10.1186/s13059-014-0409-z

  • 26.

    Cortes-LopezMMiuraP. Emerging functions of circular RNAs. Yale J Biol Med. (2016) 89:52737.

  • 27.

    CortesRFornerMJ. Circular RNAS: novel biomarkers of disease activity in systemic lupus erythematosus?Clin Sci. (2019) 133:104952. 10.1042/CS20180826

  • 28.

    LiASunYDrummerCtLuYYuDZhouYet al. Increasing upstream chromatin long-range interactions may favor induction of circular RNAs in LysoPC-activated human aortic endothelial cells. Front Physiol. (2019) 10:433. 10.3389/fphys.2019.00433

  • 29.

    NieJHLiTXZhangXQLiuJ. Roles of non-coding RNAs in normal human brain development, brain tumor, and neuropsychiatric disorders. Noncoding RNA. (2019) 5:36. 10.3390/ncrna5020036

  • 30.

    LukiwWJ. Circular RNA (circRNA) in Alzheimer's disease (AD). Front Genet. (2013) 4:307. 10.3389/fgene.2013.00307

  • 31.

    LiMDuanLLiYLiuB. Long noncoding RNA/circular noncoding RNA-miRNA-mRNA axes in cardiovascular diseases. Life Sci. (2019) 233:116440. 10.1016/j.lfs.2019.04.066

  • 32.

    ZhangHDJiangLHSunDWHouJCJiZL. CircRNA: a novel type of biomarker for cancer. Breast Cancer. (2018) 25:17. 10.1007/s12282-017-0793-9

  • 33.

    QiuBQZhangPFXiongDXuJJLongXZhuSQet al. CircRNA fibroblast growth factor receptor 3 promotes tumor progression in non-small cell lung cancer by regulating Galectin-1-AKT/ERK1/2 signaling. J Cell Physiol. (2019) 234:1125664. 10.1002/jcp.27783

  • 34.

    FangGYeBLHuBRRuanXJShiYX. CircRNA_100290 promotes colorectal cancer progression through miR-516b-induced downregulation of FZD4 expression and Wnt/beta-catenin signaling. Biochem Biophys Res Commun. (2018) 504:1849. 10.1016/j.bbrc.2018.08.152

  • 35.

    GaoDQiXZhangXFangKGuoZLiL. hsa_circRNA_0006528 as a competing endogenous RNA promotes human breast cancer progression by sponging miR-7–5p and activating the MAPK/ERK signaling pathway. Mol Carcinog. (2019) 58:55464. 10.1002/mc.22950

  • 36.

    ZhouJChenQ. Poor expression of microRNA-135b results in the inhibition of cisplatin resistance and proliferation and induces the apoptosis of gastric cancer cells through MST1-mediated MAPK signaling pathway. FASEB J. (2019) 33:342036. 10.1096/fj.201800618RRR

  • 37.

    BaiXHanGLiuYJiangHHeQ. MiRNA-20a-5p promotes the growth of triple-negative breast cancer cells through targeting RUNX3. Biomed Pharmacother. (2018) 103:14829. 10.1016/j.biopha.2018.04.165

  • 38.

    CuiJYangYLiHLengYQianKHuangQet al. MiR-873 regulates ERalpha transcriptional activity and tamoxifen resistance via targeting CDK3 in breast cancer cells. Oncogene. (2015) 34:3895907. 10.1038/onc.2014.430

  • 39.

    ZhangYYangLLingCHengW. HuR facilitates cancer stemness of lung cancer cells via regulating miR-873/CDK3 and miR-125a-3p/CDK3 axis. Biotechnol Lett. (2018) 40:62331. 10.1007/s10529-018-2512-9

  • 40.

    GaoLGuoQLiXYangXNiHWangTet al. MiR-873/PD-L1 axis regulates the stemness of breast cancer cells. EBioMedicine. (2019) 41:395407. 10.1016/j.ebiom.2019.02.034

  • 41.

    JinSHeJLiJGuoRShuYLiuP. MiR-873 inhibition enhances gefitinib resistance in non-small cell lung cancer cells by targeting glioma-associated oncogene homolog 1. Thorac Cancer. (2018) 9:126270. 10.1111/1759-7714.12830

  • 42.

    CaoDYuTOuX. MiR-873–5P controls gastric cancer progression by targeting hedgehog-GLI signaling. Pharmazie. (2016) 71:6036. 10.1691/ph.2016.6618

  • 43.

    DiaoYRahmanMFVyatkinYAzatyanASt LaurentGKapranovPet al. Identification of novel GLI1 target genes and regulatory circuits in human cancer cells. Mol Oncol. (2018) 12:171834. 10.1002/1878-0261.12366

  • 44.

    WickstromMDybergCShimokawaTMilosevicJBaryawnoNFuskevagOMet al. Targeting the hedgehog signal transduction pathway at the level of GLI inhibits neuroblastoma cell growth in vitro and in vivo. Int J Cancer. (2013) 132:151624. 10.1002/ijc.27820

  • 45.

    XuLWangXWanJLiTGongXZhangKet al. Sonic hedgehog pathway is essential for neuroblastoma cell proliferation and tumor growth. Mol Cell Biochem. (2012) 364:23541. 10.1007/s11010-011-1222-6

  • 46.

    SchiapparelliPShahiMHEnguita-GermanMJohnsenJIKognerPLazcozPet al. Inhibition of the sonic hedgehog pathway by cyplopamine reduces the CD133+/CD15+ cell compartment and the in vitro tumorigenic capability of neuroblastoma cells. Cancer Lett. (2011) 310:22231. 10.1016/j.canlet.2011.07.005

  • 47.

    XuLKwonYJFrolovaNStegADYuanKJohnsonMRet al. Gli1 promotes cell survival and is predictive of a poor outcome in ERalpha-negative breast cancer. Breast Cancer Res Treat. (2010) 123:5971. 10.1007/s10549-009-0617-5

  • 48.

    ChengJDengRWuCZhangPWuKShiLet al. Inhibition of SALL4 suppresses carcinogenesis of colorectal cancer via regulating Gli1 expression. Int J Clin Exp Pathol. (2015) 8:10092101.

  • 49.

    ChenGGotoYSakamotoRTanakaKMatsubaraENakamuraMet al. GLI1, a crucial mediator of sonic hedgehog signaling in prostate cancer, functions as a negative modulator for androgen receptor. Biochem Biophys Res Commun. (2011) 404:80915. 10.1016/j.bbrc.2010.12.065

  • 50.

    ClementVSanchezPde TriboletNRadovanovicIRuiziAltabaA. HEDGEHOG-GLI1 signaling regulates human glioma growth, cancer stem cell self-renewal, and tumorigenicity. Curr Biol. (2007) 17:16572. 10.1016/j.cub.2006.11.033

  • 51.

    GanHLiuHZhangHLiYXuXXuXet al. SHh-Gli1 signaling pathway promotes cell survival by mediating baculoviral IAP repeat-containing 3 (BIRC3) gene in pancreatic cancer cells. Tumour Biol. (2016) 37:994350. 10.1007/s13277-016-4898-0

  • 52.

    SunLYaoYPanGZhanSShiWLuTet al. Small interfering RNA-mediated knockdown of fatty acid synthase attenuates the proliferation and metastasis of human gastric cancer cells via the mTOR/Gli1 signaling pathway. Oncol Lett. (2018) 16:594602. 10.3892/ol.2018.8648

  • 53.

    NayakASatapathySRDasDSiddharthSTripathiNBharatamPVet al. Nanoquinacrine induced apoptosis in cervical cancer stem cells through the inhibition of hedgehog-GLI1 cascade: role of GLI-1. Sci Rep. (2016) 6:20600. 10.1038/srep20600

Summary

Keywords

neuroblastoma, circDGKB, miR-873, GLI1, diagnostic marker

Citation

Yang J, Yu L, Yan J, Xiao Y, Li W, Xiao J, Lei J, Xiang D, Zhang S and Yu X (2020) Circular RNA DGKB Promotes the Progression of Neuroblastoma by Targeting miR-873/GLI1 Axis. Front. Oncol. 10:1104. doi: 10.3389/fonc.2020.01104

Received

09 December 2019

Accepted

02 June 2020

Published

23 July 2020

Volume

10 - 2020

Edited by

Jian-Guo Zhou, University of Erlangen Nuremberg, Germany

Reviewed by

Miguel F. Segura, Vall d'Hebron Research Institute (VHIR), Spain; Malin Wickstrom, Karolinska Institutet (KI), Sweden

Updates

Copyright

*Correspondence: Shouhua Zhang Xin Yu

This article was submitted to Cancer Molecular Targets and Therapeutics, a section of the journal Frontiers in Oncology

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics