- Primary research
- Open Access
Silencing of SNHG6 induced cell autophagy by targeting miR-26a-5p/ULK1 signaling pathway in human osteosarcoma
Cancer Cell Internationalvolume 19, Article number: 82 (2019)
lncRNAs have been proved to play crucial parts in various human cytopathology and cell physiology, including tumorigenesis. Down-regulated lncRNAs SNHG6 have shown great cell proliferation inhibitory effects in cancer development. Here we investigated how SNHG6 effected human osteosarcoma (OS) development and progression. Methods: Reverse transcription-quantitative PCR was performed to detect SNHG6 mRNA level in both OS tissues and cell lines. MTT and colony formation assays were used to determine the growth impact of SNHG6. Wound healing and trans-well assay were performed to measure the invasion effect of SNHG6. Western blotting were utilized to dissect molecular mechanisms.
We identified SNHG6 as a lncRNAs that significantly up-regulated in OS tissues and cells, patients with high SNHG6 expression suffered more malignant metastasis and shorter survival times. Furthermore, silencing of SNHG6 in OS significantly inhibited OS cell growth, weakened cell invasion capacity, arrested cell cycle at G0/G1 phase, and induced cell apoptosis. Additionally, mechanism assays suggested that SNHG6 could competitively sponging miR-26a-5p thereby regulating ULK1, and induced cell apoptosis and autophagy by targeting caspase3 and ATF3. Conclusions: Our findings demonstrated that SNHG6 acted as an oncogene in osteosarcoma cells through regulating miR-26a-5p/ULK1 at a post-transcriptional level. SNHG6 might serve as a candidate prognostic biomarker and a target for novel therapies of osteosarcoma patients.
Osteosarcoma (OS), which remains the most frequent human primary solid malignancy bone cancer, has been suggested to be a differentiation disease caused by malignant mesenchymal cells differentiation interrupt [1, 2]. OS mainly presents in children and young adult at 15–19 years of age, accounting for almost 60% of malignant bone tumors, moreover, the incidence of OS was about 8–11/million/year throughout the world [3, 4]. OS has highly aggressive and early metastatic potential . Approximately 10%–25% of patients occur lung metastasis and die of pulmonary metastasis within 1 year . The use of effective adjuvant chemotherapies plays an important role in controlling OS. Combination of chemotherapeutics and surgery resection have been used for human OS treatment, which significantly increased OS patients’ survival ratefrom 20% at 1979s to 65%–75% at 2012, however less than 30% patients with metastasis disease can survive [7,8,9,10]. Thus, it was necessary to develop more effective adjuvant treatments.
Long non-coding RNAs (lncRNAs) were numerous transcripts that consist of over 200 nucleotides. Though 90% human genome DNA were widely transcribed, 98% of them encode none protein, which were identified as non-coding RNAs (ncRNAs) [11, 12]. LncRNAs make up more than 70% of ncRNAs,  and they have received more and more attention for the new function in tumorigenesis, metastasis and chemoresistance . Therefore, more researchers were seeking for lncRNAs which might contribute to the pathogenesis or chemoresistance of OS.
During our research, we found that nucleolar RNA host gene 6 (SNHG6) was over-expressed in primary tumor samples and OS cell lines (SOSP-9607 and MG63). Furthermore, SNHG6 was identified as a tumor promoter targeting miR-26a-5p, as reduction of SNHG6 expression in OS cell lines was able to reduce cell viability and migration in vitro, arrest cell cycle at G0/G1 phase, induce cell apoptosis. Additionally, unc-51 like autophagy activating kinase (ULK1), a key component of autophagy pathway, was identified to be a novel direct target of miR-26a-5p. And mechanism experiments uncovered that SNHG6 might functioned through sponging miR-26a-5p and thus regulating cell apoptosis and cell autophagy by targeting ULK1. Collectedly, we might find a novel therapy strategy based on lncRNA-directed against human OS.
Patients and tumor samples
Osteosarcoma tumor tissues and adjacent normal tissues were collected from 30 primary osteosarcoma patients among 2006 and 2009 at Shangqiu First People’s Hospital. Following surgery, tissues samples were stored at liquid nitrogen until subsequent experiments.
Informed consents were obtained from all individual patients involved in this research. All procedures performed in studies involving human participants were in accordance with the ethical standards of the Ethics Review Committee of Shangqiu First People’s Hospital (Shangqiu, China). The characteristics of these patients were shown in Table 1.
Human normal osteoblastic cell line hFOB, human osteosarcoma cell lines SOSP-9607 and MG63 cells were obtained from American Type Culture Collection (ATCC). SOSP-9607 cells was maintained in RPMI 1640 medium (HyClone), 2.0 mM l-glutamine, 100 U/ml penicillin, and 100 ug/ml streptomycin, while other cells were maintained in the same conditions, except that Dulbecco’s modified Eagle’s medium (DMEM, Sigma) was used. All the media were supplemented with 10% dialyzed fetal bovine serum (FBS; HyClone) at 37 ̊C in a humidified incubator containing 5% CO2.
hFOB, SOSP-9607 and MG63 cells were seeded into each well of 6-well plate of 5 × 103 cells, respectively. Then, transfected with human SNHG6 small interfering RNA (siRNAs) (Cat. No.: iV022927) or negative control (NC) siRNA, which were purchased from Applied Biological Materials Inc. (abm) and the corresponding sequences were as follows: siRNA-SNHG6, CGCGAAGAGCCGTTAGTCATGCCGGTGTG; siRNA-NC, UUCUC CGAACGUGUCACGUTT. All constructs were transfected into cells using Lipofectamine 2000 (Invitrogen), according to the manufacturer’s instructions.miR-26a-5p mimic and miR-26a-5p inhibitor were obtained from Beyotime (Beijing, China). Following incubation for 48 h, the transfected cells were har- vested and used for further experiments.
RNA extraction and qRT-PCR
After 48 h post transfection, total RNA was extracted by using Trizol reagent (Invitrogen), respectively. RNA was resuspended in DEPC-treated H2O, and the concentration and purity were confirmed at 260 nm. For reverse transcription PCR, 10 mg of total RNA was converted to cDNA using cDNA Synthesis Kit (Bio-Rad), then PCR were performed. PCR was run at 25–30 cycles at 95 °C for 30 s, 56 °C for 30 s, and 72 °C for 1 min. U6 was used as an internal control. The relative expression levels were calculated by the 2−ΔΔCt method. The primer sequences were designed by Beyotime (Beijing) and listed as follows: SNHG6, Forward primer: 5′-CCTACTGACAACATCGACGTTGAAG-3′; Reverse primer: 5′-GGAGAAAACGCTTAGCCATACAG-3′. miR- 26a-5p forward: 5′-GCGCGCGTAACAGTCTACAGC-3′ and reverse: 5′-GTCGTATCCAGTGCAGGGTCC-3′ U6: forward primer: 5′-CTCGCTTCGGCA GCA CA-3′ and reverse primer: 5′-AACGCTTCACGAATTTGCGT-3′. All reactions were performed at least three time.
Western blot analysis
After 48 h post transfection, transfected cells were washed by ice-cold PBS and subjected to lysis in a RIPA buffer with 0.5% sodium dodecyl sulfate (SDS), 3% proteinase inhibitor cocktail (Sigma). Protein concentration was detected by BCA (bicinchoninic acid) protein assay. Then 15% sodium dodecyl sulfate polyacryl-amide gel electrophoresis gels (SDS-PAGE) and PVDF membranes (Millipore) were prepared to separate different protein. After blocking in 5% (w/v) non-fat dry milk for 1 h, PVDF membranes were incubated with primary antibodies against ULK1 (Cat. No.: ab203207), or Caspase 3 (Cat. No.: ab13585), or Activating transcription factor 3 (ATF3) (Cat. No.: ab207434) or GAPDH (Cat. No.: ab8245; abcam) overnight at 4 °C. Then membranes were incubated with horseradish peroxidase conjugated anti-mouse or rabbit secondary antibody (abcam) for 2 h. Finally, membranes were detected and analyzed using ECL Western Blotting Substrate (Thermo Fisher Scientific).
Cell viability assay
3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT; Sigma) assay were used to detect cell viability. Cells transfected with SNHG6 siRNA or NC siRNA were seeded into 96 well plates and cultured overnight. Cells transfected nothing (Blank group) were used as control. Every 24 h interval, 20 ml of 5 mg/ml MTT was added into each well and cells were incubated for 2 h at 37 °C. Then 200 ml of DMSO (dimethyl sulfoxide; Sigma) were added to replace the medium. Optical density (OD) was evaluated at wavelength of 490 nm by measuring the absorbance. Cells treated with DMSO alone were performed to be a control. % viability survival = optical density of each group/optical density of DMSO group.
Each test was performed for five consecutive days and repeated in five wells. The experiments were repeated three times independently.
Colony formation assay
Cells (500 cells per well) transfected with indicated vector were plated in 6-well plates and incubated at 37 °C. After 2 weeks, cells were fixed and stained with 0.1% crystal violet. Numbers of visible colonies were then counted manually after washing.
Migration scratch assay
Transfected OS cells were placed on 12-well plates, then a cell‐scratch was made. Cell debris were washed by phosphate buffered saline (PBS) and cells were incubated at 37 °C for another 24 h. The rate of cell migration was determined by comparison of the images at 0 and 24 h after scratching. Each experiment was repeated three times.
Trans-well migration assays
Boyden chamber with 8.0 μm pore (Corning) were used for this assay. 2 × 105 transfected OS cells were seeded into the upper compartment in serum-free media. Media containing 10% FBS was placed in the lower chamber. Following 48 h incubation at 37 °C in 5% CO2, cells remaining in upper surface of the filter were wiped off with a cotton swab, while cells that migrated were fixed with methanol and stained with Giemsa (Sigma). Invasive cells were counted in five randomly light-microscopy fields. Each experiment was repeated in triple.
Luciferase reporter assay
Wild-type (WT) fragments from SNHG6 (SNHG6-Wt) containing relative predicted miR-26a-5p binding sites or mutated-type (MUT) SNHG6 (SNHG6-Mut) sequence were synthesized and cloned downstream of the Dual-luciferase gene in the pGL3-basic vector (Promega, Madison, WI, USA). The cells in the logarithmic growth were inoculated into the 96-well plate and transfected using Lipofectamine 2000 (Invitrogen, NY, USA) when the cell density was about 70%. Dual-luciferase reporter plasmid containing SNHG6-Wt or SNHG6-Mut was co-transfected into HEK293T cells along miR-26a-5p mimics in parallel with an empty plasmid vector (miR-NC). After 24 h of transfection, relative luciferase activity was determined using a Dual luciferase assay kit (BioVision, San Francisco, CA, USA). The firefly luciferase activity was detected and normalized to Renilla activity. The experiment was repeated three times.
Similarly, WT or MUT ULK1 sequence containing putative miR-26a-5p recognition sites was cloned downstream of the luciferase gene in the Dual-luciferase Vector to form the reporter vector ULK1-wild-type (ULK1-Wt) or ULK1-mutated-type (ULK1-Mut).
Flow cytometry assay
For cell-cycle analysis, OS cells were seeded in 6 cm dishes at 5 × 104/dish after transfection. After 48 h of cultivation, cells were washed and suspended in PBS containing 50 μg/ml RNase A (Sigma-Aldrich) and 50 μg/ml propidium iodide (PI; BD Biosciences, Erembodegem, Belgium) for 1 h in the dark. Fluorescence associated cellsorting (FACS) Verse™ flow cytometer (Becton–Dickinson, CA, USA) was used to count cell populations in each cell cycle phases and the results were analyzed with the ModFit (Verity Software House, Maine, USA) software.
For apoptosis analysis, cells (1 × 105/well) were collected after 48 h of transfection and stained with both Annexin V–fluorescein isothiocyanate (FITC) and PI according to the manufacturer’s instructions (BD). Apoptosis was detected by FACSVerse™ flow cytometer.
Data were analyzed by GraphPad Prism 6.0 (GraphPad) and presented as mean ± standard deviation (SD). All P values were calculated using unpaired Student’s t test or ANOVA test with Tukey’s post hoc test, unless otherwise stated. Paired Student’s t-test (for parametric data) or Wilcoxon signed-rank test (for non-parametric data) were applied to analyze the differences of SNHG6 expression levels between OS tissues and adjacent normal tissues. The Pearson’s Chi square were performed to determine the difference between SNHG6 expression levels and clinic pathological factors. The overall survival probability was analyzed using Kaplan–Meier methods and evaluated using the log-rank test. P < 0.05 was considered as statistically significant. 3. Results
SNHG6 was over-expressed in human osteosarcoma tissues
To investigate whether SNHG6 was correlated with the progression of OS patients, we explored SNHG6 expression in The Cancer Genome Atlas (TCGA) Data Portal from starBASE v2.0 . We found that SNHG6 level was much higher in tumor tissues than in normal tissues among Pan-Cancer including 14 cancer types (breast cancer, clear cell kidney carcinoma, lung cancer, and so on) (Fig. 1a). Then we analyzed 30 pairs of human OS tissues samples (30 human OS tissues and 30 corresponding adjacent normal tissues) by using qRT-PCR. Comparing with matched adjacent normal tissues, SNHG6 was significantly over-expressed in human OS cancer tissues by fourfold (Fig. 1b, P < 0.001). Besides, over expression of SNHG6 level was also observed in two human OS cell lines (SOSP-9607 and MG63 cells) comparing with human hFOB cells (Fig. 1c, P < 0.001).
High level SNHG6 was correlated with poor overall survival in OS patients
To further understand the association of SNHG6 in OS, we divided all OS tissue samples into 2 groups using median level of SNHG6 as a cutoff value (high SNHG6 expression group: SNHG6 expression ratio > median level; low SNHG6 expression group: SNHG6 expression ratio ≤ median level). This finding demonstrated that higher SNHG6 level was highly correlated with tumor invasion depth (Table 1, P = 0.001), lymph node metastasis (Table 1, P = 0.011) and TNM stage (Table 1, P = 0.001), while there was no significant correlation with gender, age and histological grade (P > 0.05). Additionally, Kaplan–Meier method analysis and log-rank test showed that patients with high level of SNHG6 processed a shorter overall survival (Fig. 1d, P < 0.01). All these data implicated that SNHG6 might implicated a poor prognosis.
SNHG6 knockdown suppressed OS cell proliferation
To clarify the function of SNHG6 in OS, we reduced SNHG6 level in SOSP-9607 and MG63 cells by transfecting siRNA specifically targeting SNHG6 (siRNA-SNHG6). qRT-PCR assay was performed to confirmed the transfection effect (Fig. 2a, p < 0.01). Knockdown of SNHG6 decreased the SNHG6 expression levels by fourfold. MTT and colony formation assay were performed after transfection, and which suggested that knockdown of SNHG6 markedly suppressed cell proliferation (Fig. 2b, P < 0.05) and colony formation (Fig. 2c, P < 0.05) in both OS cell lines.
SNHG6 knockdown suppressed OS cell migration and invasion
Migration and invasion experiments were performed on SOSP-9607 and MG63 cells transfected with siRNA-SNHG6 or siRNA-NC. SNHG6 silencing significantly inhibited OS cells migration (Fig. 3a, P < 0.05) and invasion (Fig. 3b, P < 0.05), respectively.
SNHG6 knockdown arrested cell cycle and cell apoptosis
Next, flow cytometric analysis were carried out to characterize the functions of SNHG6 on OS cell cell-cycle or apoptosis, flow cytometric analysis was performed. Based on these results, a significant arrested at the G1/G0 phase were observed on OS cells transfected with siRNA-SNHG6, comparing with cells transfected with siRNA-NC (Fig. 4, P < 0.01). In the meantime, apoptosis cells increased significantly in OS cells transfected with siRNA-SNHG6 (Fig. 5, P < 0.01). Thus, we confirmed that SNHG6 silencing suppressed OS cell proliferation by inhibited cell cycle progression and induced cell apoptosis.
miR-26a-5p was an inhibitory target of SNHG6
Recently, various researches proved that lncRNAs mainly functioned as miRNA sponges by binding functional miRNAs and inhibited the transcription affection [16, 17]. To confirm whether SNHG6 played the similar role in OS progression, we predicted the potential miRNAs associated with SNHG6 by the online software starBase v2.0, 5 miRNAs, including miR-101-3p, miR-26a-5p, miR-1297, miR-26b-5p and miR-4465 were the potential targets of SNHG6, while none of them have been implicated in OS progression. To further confirmed this prediction, a Dual-luciferase reporter assay was used in 293T cells. The pGL3 Basic Luciferase Reporter Vectors was shown in Fig. 6a. The predicted target sites of miR-26a-5p binding to the 5′-UTR of SNHG6 mRNA was showed in Fig. 6b. Co-transfection of miR-26a-5p and SNHG6-Wt but not the SNHG6-Mut induced a significantly inhibition on luciferase activity (Fig. 6a, P < 0.05), indicating miR-26a-5p might be a direct target of SNHG6 in OS progression.
To further support this conclusion, qRT-PCR assay were performed. SNHG6 silencing induced a significantly increased expression of miR-26a-5p (Fig. 6c, P < 0.01), while over-expressed of miR-26a-5p did not act on the expression of SNHG6 (Fig. 6d, P > 0.05). In a summary, our research suggested that miR-26a-5p was a direct inhibitory target for SNHG6 in OS progression.
ULK1 was a target gene of miR-26a-5p and was regulated by SNHG6
ULK1, a serine/threonine protein kinase as the key component of autophagy pathway, came out as one of the target genes for miR-26a-5p . To confirm the regulation of ULK1 by miR-26a-5p, Dual-luciferase reporter assay were performed in 293T cells. The predicted target sites of miR-26a-5p binding to the 5′-UTR of ULK1 mRNA was showed in Fig. 7a. The luciferase activity was significantly inhibited while co-transfected with miR-26a-5p mimic and ULK1–3′UTR-wt plasmids (Fig. 7a, P < 0.05). To further confirm the regulation of ULK1 by miR-26a-5p, we examined the protein expression of ULK1 in OS cells transfected with miR-26a-5p, miR-26a-5p inhibitor, siRNA-SNHG6 or both miR-26a-5p inhibitor and siRNA-SNHG6. Our results in Fig. 7b showed that increased expression of miR-26a-5p markedly increased the protein levels of ULK1 and this could been reversed by miR-26a-5p inhibitor. Knockdown of SNHG6 reduced the protein expression of ULK1, while knockdown of SNHG6 was no longer able to increase ULK1 protein level with miR-26a-5p inhibition. Collectively, these data might imply that ULK1 mRNA was the target of miR-26a-5p.
SNHG6 regulated autophagy by targeting miR-26a-5p/ULK1
Autophagy played opposed survival-supporting and death-promoting roles in cancer, which might be a tumor suppressor for preventing tumorigenesis, or act as an oncogene to promote cancer progression . ULK1 and its complex has been reported as an essential role in autophagy induction [20, 21]. Activating transcription factor 3 (ATF3) was an autophagy-related protein, and caspase3 was known as the apoptotic executor, both were modulated by ULK1 in the downstream signals . We then explored the relationship between SNHG6 and the ULK1 signaling pathway in OS cells by western blot analysis. As Fig. 8a shown, cells transfected with siRNA-SNHG6 had significantly decreased Caspase 3 exposure levels. Moreover, the expression of ATF3 in OS cells was significantly suppressed after SNHG6 siRNA transfection. All these results demonstrated that caspase3 and ATF3 were all involved in human OS cells death, which may be regulated by SNHG6/miR-26a-5p/ULK1 signaling pathway (Fig. 8b).
In the past decades, a substantial body of research have indicated the critical role of lncRNAs in carcinogenesis, included in osteosarcoma [23,24,25]. We focused on SNHG6 function because previous studies have demonstrated the potential oncogene function of SNHG6 to tumorigenesis. Recently, Chang et al. demonstrated that SNHG6 served as an oncogene by competitively binding miR-101-3p in human hepatocellular carcinoma . Soon afterwards, Kai et al. demonstrated that SNHG6 promoted gastric cancer cells proliferation via miR-101-3p/ZEB1 regulation which performed at a post-transcriptional level . However, more role of SNHG6 in human OS urgently need to be discovered.
Among this research, function and mechanism experiments were performed. Firstly, we found that SNHG6 expression was significantly increased in OS tissues comparing with corresponding normal tissues. The increased SNHG6 expression was determined to be significantly correlated with invasion depth and TNM stage, and positively correlated with poor prognosis of OS patients. Then loss-of-function assay was performed, and we found that silenced SNHG6 could decreased OS cell proliferation, migration and invasion. This was the first report of SNHG6 function in human OS, which may present a novel therapeutic strategy.
lncRNAs have been reported to regulate tumor development and progression by modulating different targets . The most common regulation patten of lncRNAs was to be a direct targets of miRNAs as miRNA molecular sponge during the last few year, moreover, miRNA altered the expression of target mRNA [29, 30]. The lncRNAs-miRNA-mRNA axis seemed to be an extensive regulatory network in gene expression, and the extraordinary gene expression always indicated a disease progression, including cancer development. Likewise, in this investigation, SNHG6 silencing was attributed to its influence depending on its regulation on ULK1 by sponging suppress miR-26a-5p. Collectively, this is the first time that SNHG6 was identified to be a member of miR-26a-5p/ULK1 regulation axis in human OS development and progression, which probably offered an optional therapeutics against OS.
Autophagy was a process of cell self-digestion, which help to degrade and re-use those unwanted or dysfunctional components in cells, and turn out to help to kept cells in a state of dynamic equilibrium . Autophagy acted as a double-edged sword in cancer. On the one hand, autophagy seemed to play a role as a suppressor by clearing mutant or misfolded proteins and alleviating cellular stress at the beginning of tumor development,  on the other hand, autophagy enabled tumor cells to survive under nutrient deficiency, hypoxia and other stresses, or to form a barrier to treatment with chemotherapeutic drugs. Thus, we speculated autophagy block might be beneficial to cancer therapy. However, most studies have focused on using common autophagy inhibitors, and most inhibiting targets were far downstream of the autophagy process . ULK1, a serine/threonine protein kinase, which was one of the crucial factors for autophagy initiation and played a cardinal role in the regulation of autophagy . Accumulating evidence suggestedULK1 act as a key initiator of autophagy. Zheng et al. firstly found that MiR-26a-5p regulate autophagic pathway in cardiac broblasts by targeting ULK1 . Jin et al. demonstrated for the first time that miR-26a/b can promote apoptosis and sensitize hepatocellular carcinoma to chemotherapy via suppressing the expression of autophagy initiator ULK1 .
As the most common form of childhood and adolescent cancer, osteosarcoma suffered a rare survival rate for drug resistance. Inhibition of autophagy has been proved to increase the drug sensitivity of osteosarcoma cells in vivo and in vitro . So we speculated whether SNHG6 promoted OS cells proliferation and reduced cell apoptosis through cell autophagy. Caspase-3 has been considered as crucial modulators of apoptosis. In our research, the expression level of Caspase-3 increased significantly following SNHG6 siRNA treatment, indicating that SNHG6 suppression promoted cell death by activating caspase-3-dependent apoptosis. Furthermore, SNHG6 siRNA significantly up-regulated the expressions of ATF3, implicating caspase3 and ATF3 were involved in SNHG6 siRNA induced cell death, which may be simultaneously regulated by miR-26a-5p/ULK1.
In summary, we have newly identified the SNHG6/miR-26a-5p/ULK1 regulatory mechanism of autophagy in human osteosarcoma, and provided an optional choose for human OS treatment in the future.
nucleolar RNA host gene 6
long non-coding RNAs
activating transcription factor 3
Sobin L, Wittekind C. TNM classification of malignant tumours, (UICC). New York: Springer-Verlag; 2002.
Fletcher CDM, Unni KK, Mertens F. World Health Organization classification of tumours. Pathology and genetics of tumours of soft tissue and bone. Lyon: IARC Press; 2004. p. 360–2
Stiller CA, Stevens MC, Magnani C, Corazziari I. Survival of children with soft-tissue sarcoma in Europe since 1978: results from the EUROCARE study. Eur J Cancer. 2001;37(6):767–74.
Stiller CA, Bielack SS, Jundt G, Steliarova-Foucher E. Bone tumours in European children and adolescents, 1978–1997. Report from the Automated Childhood Cancer Information System project. Eur J Cancer. 2006;42(13):2124–35.
Tan ML, Choong PF, Dass CR: Osteosarcoma: Conventional treatment vs. gene therapy. Cancer Biology & Therapy 2009, 8(2):106-117.
Miller BJ, Cram P, Lynch CF, Buckwalter JA. Risk factors for metastatic disease at presentation with osteosarcoma: an analysis of the SEER database. J Bone Joint Surg Am. 2013;95(13):e89.
Daw NC, Billups CA, Rodriguezgalindo C, Mccarville MB, Rao BN, Cain AM, Jenkins JJ, Neel MD, Meyer WH. Metastatic osteosarcoma. Cancer. 2006;106(2):403.
Iwamoto Y, Tanaka K, Isu K, Kawai A, Tatezaki S, Ishii T, Kushida K, Beppu Y, Usui M, Tateishi A. Multiinstitutional phase II study of neoadjuvant chemotherapy for osteosarcoma (NECO study) in Japan: NECO-93J and NECO-95J. J Orthopaed Sci. 2009;14(4):397–404.
Jaffe N. Osteosarcoma: review of the past, impact on the future the american experience. Cancer Treatment Res. 2009;152(152):239.
Allison DC, Carney SC, Ahlmann ER, Hendifar A, Chawla S, Fedenko A, Angeles C, Menendez LR. A Meta-analysis of osteosarcoma outcomes in the modern medical era. Sarcoma. 2012;2012(4):704872.
Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, Tanzer A, Lagarde J, Lin W, Schlesinger F. Landscape of transcription in human cells. Nature. 2012;489(7414):101–8.
Martens-Uzunova ES, Böttcher R, Croce CM, Jenster G, Visakorpi T, Calin GA. Long noncoding RNA in prostate, bladder, and kidney cancer. Eur Urol. 2014;65(6):1140.
Yang Q, Zhang S, Liu H, Wu J, Xu E, Peng B, Jiang Y. Oncogenic role of long noncoding RNA AF118081 in anti-benzo[a]pyrene-trans-7,8-dihydrodiol-9,10-epoxide-transformed 16HBE cells. Toxicol Lett. 2014;229(3):430–9.
Chen R, Wang G, Zheng Y, Hua Y, Cai Z. Long non-coding RNAs in osteosarcoma. Oncotarget. 2017;8(12):20462.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v20: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92.
Paraskevopoulou MD, Georgakilas G, Kostoulas N, Reczko M, Maragkakis M, Dalamagas TM, Hatzigeorgiou AG. DIANA-LncBase: experimentally verified and computationally predicted microRNA targets on long non-coding RNAs. Nucleic Acids Res. 2013;41(1):D239.
Han D, Li J, Wang H, Su X, Jin H, Yan G, Cheng Q, Lin Y, Liu X, Huang M. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017;66(4):1151.
Zheng L, Lin S, Lv C. MiR-26a-5p regulates cardiac fibroblasts collagen expression by targeting ULK1. Sci Rep. 2018;8:2104.
Mathew R, Karantzawadsworth V, White E. Role of autophagy in cancer. Med Recapitulate. 2010;7(12):961.
Kim J, Kundu M, Viollet B, Guan KL. AMPK and mTOR regulate autophagy through direct phosphorylation of Ulk1. Nat Cell Biol. 2011;13(2):132–41.
Yuan J, Zhao X, Hu Y, Sun H, Gong G, Huang X, Chen X, Xia M, Sun C, Huang Q. Autophagy regulates the degeneration of the auditory cortex through the AMPK-mTOR-ULK1 signaling pathway. Int J Mol Med. 2018;41(4):2086–98.
Tattoli I, Sorbara MT, Philpott DJ, Girardin SE. Bacterial autophagy: the trigger, the target and the timing. Autophagy. 2012;8(12):1848.
Sun XH, Yang LB, Geng XL, Wang R, Zhang ZC. Increased expression of lncRNA HULC indicates a poor prognosis and promotes cell metastasis in osteosarcoma. Int J Clin Exp Pathol. 2015;8(3):2994–3000.
Tao Z, Gao Y. Increased expression of LncRNA BANCR and its prognostic significance in human hepatocellular carcinoma. World J Surg Oncol. 2015;14(1):1–8.
Wei L, Peng X, Ruan WH. Overexpression of lncRNA UCA1 promotes osteosarcoma progression and correlates with poor prognosis. J Bone Oncol. 2016;5(2):80–5.
Chang L, Yuan Y, Li C, Guo T, Qi H, Xiao Y, Dong X, Liu Z, Liu Q. Upregulation of SNHG6 regulates ZEB1 expression by competitively binding miR-101-3p and interacting with UPF1 in hepatocellular carcinoma. Cancer Lett. 2016;383(2):183–94.
Kai Y, Jie T, Shi W, Hao X, Zhu Y. LncRNA SNHG6 is Associated with Poor Prognosis of Gastric Cancer and Promotes Cell Proliferation and EMT through Epigenetically Silencing p27 and Sponging miR-101-3p. Cell Physiol Biochem. 2017;42(3):999.
Du WW, Yang W, Liu E, Yang Z, Dhaliwal P, Yang BB. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016;44(6):2846.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.
Mizushima N, Komatsu M. Autophagy: renovation of cells and tissues. Cell. 2011;147(4):728–41.
Yang ZJ, Cheng EC, Huang S, Sinicrope FA. The role of autophagy in cancer: therapeutic implications. Mol Cancer Ther. 2011;10(9):1533–41.
Kabeya Y, Mizushima N, Ueno T, Yamamoto A, Kirisako T, Noda T, Kominami E, Ohsumi Y, Yoshimori T. LC3, a mammalian homologue of yeast Apg8p, is localized in autophagosome membranes after processing. EMBO J. 2000;19(21):5720–8.
Tannous P, Zhu H, Nemchenko A, Berry JM, Johnstone JL, Shelton JM, Miller FJ, Rothermel BA, Hill JA. Intracellular protein aggregation is a proximal trigger of cardiomyocyte autophagy. Circulation. 2008;117(24):3070–8.
Jin F, Wang Y, Li M, Zhu Y, Liang H, Chen W, Feng W, Zhang CY, Ke Z, Li L. MiR-26 enhances chemosensitivity and promotes apoptosis of hepatocellular carcinoma cells through inhibiting autophagy. Cell Death Dis. 2017;8(1):e2540.
Huang J, Liu K, Yu Y, Xie M, Kang R, Vernon P, Cao L, Tang D, Ni J. Targeting HMGB1-mediated autophagy as a novel therapeutic strategy for osteosarcoma. Autophagy. 2012;8(2):275–7.
All authors contributed equally in this study. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the Ethics Review Committee of Shangqiu First People’s Hospital (Shangqiu, China).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.