Skip to main content
  • Research Article
  • Open access
  • Published:

LncRNA EBLN3P attributes methotrexate resistance in osteosarcoma cells through miR-200a-3p/O-GlcNAc transferase pathway



Osteosarcoma is highly malignant. The migration, invasion, and chemoresistance contribute to poor prognosis of osteosarcoma. Research reported that endogenous bornavirus-like nucleoprotein 3 pseudogene (EBLN3P) promotes the progression of osteosarcoma.


In this study, the expression of EBLN3P in osteosarcoma tissue with different methotrexate (MTX) treatment responses was measured. Osteosarcoma cell lines with MTX resistance were constructed, and bioinformatic analysis was performed to explore the potential involved targets and pathways.


Higher EBLN3P was associated with MTX resistance. Downregulation of LncEBLN3P decreased the MTX resistance of osteosarcoma cells by sponging miR-200a-3p, an important microRNA that affects epithelial-mesenchymal transition (EMT). The decreased miR-200a-3p resulted in the upregulation of its target gene O-GlcNAc transferase (OGT), which in turn promoted the EMT process of osteosarcoma cells. Further analysis confirmed that the loss of OGT and over-expression of miR-200a-3p could partly abolish the MTX resistance induced by LncEBLN3P.


LncEBLN3P is upregulated in osteosarcoma and increases the MTX resistance in osteosarcoma cells through downregulating miR-200a-3p, which in turn promoted the EMT process of osteosarcoma cells by increasing the OGT.


Osteosarcoma is one of the malignant orthopedic tumors that occur in young under the age of 30. The advance in neoadjuvant chemotherapy increases the 5-year overall survival rate by over 60% [1]. Currently, the commonly used treatment is surgery plus chemotherapy. However, the resistance to chemotherapy affects the curative effect and prognosis of osteosarcoma [2]. The treatment of new adjuvant chemotherapy drugs and immunotherapy has not achieved great breakthroughs yet [3], and it is particularly important to find the molecular targets involved in drug resistance to osteosarcoma. The commonly used reagents for chemotherapy include cisplatin, methotrexate (MTX), ifosfamide, and doxorubicin [4]. Increasing the dose intensity of MTX during chemotherapy can improve patients' 5-year tumor-free survival rate, but high-dose MTX has potential toxic side effects and even induces drug resistance [5]. Metastasis and chemotherapy resistance become the main factors restricting the overall survival rate in patients with osteosarcoma.

LncRNA is a group of noncoding RNA with a length of more than 200 bp. The transcription of lncRNA is generated from independent promoters by RNA polymerase II. They participate in the progression of various cancers and drug resistance. As an important epigenetic regulator, LncRNAs mainly regulate the expression of related functional proteins via miRNA [6]. And the miRNA-based therapies have great potential compared with traditional therapies for diseases treatment, including musculoskeletal diseases [7,8,9]. LncEBLN3P is a pseudogene of endogenous bornavirus-like nucleoprotein 3 (EBLN3) located on chromosome 9. Research has demonstrated that EBLN3P is over-expressed in tumors [10,11,12], including osteosarcoma [13].

Increased O-GlcNAc transferase (OGT) expression has been reported in many human cancers [14]. A domestic study in China reported the expression of OGT in osteosarcoma tissues increased and downregulating its expression increased the sensitivity of osteosarcoma cells to cisplatin [15]. Deng et al. [16] reported that OGT degradation contributed to tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) resistance. The epithelial-to-mesenchymal transition (EMT) increases the motility of cancer cells. EMT also plays important role in maintaining stemness properties, which contributes to the chemoresistance of cancer cells [17]. Loss of miR-200 expression is correlated with the EMT process of tumors [18, 19]. An in vitro analysis indicated that overexpression miR-200a could inhibit osteosarcoma cell proliferation, promotes apoptosis, and increases cellular radiosensitivity [20]. In this study, we used clinical specimens, further with the help of in vitro molecular methods, upregulation or knocking-down LcnEBLN3P expression in osteosarcoma cells, to investigate the biofunction of LncEBLN3P in invasion, migration, and MTX resistance of osteosarcoma cells, and finally explored whether it play its function via the miR-200a-3p/O-GlcNAc transferase signaling axis.


Osteosarcoma tissue samples and cell lines

The osteosarcoma tissues and the adjacent non-tumor tissues were obtained from the surgical resection sample of our hospital. The freshly resected tissues were preserved at − 80 °C refrigerators or fixed for pathological examination. All patients were diagnosed with osteosarcoma by pathological examination and did not receive any anti-tumor treatment before the operation. The informed consent was obtained from each patient and the study was approved by the ethics committee of our hospital.

In this study, patients who received cisplatin + ifosfamide + high dose methotrexate regimen after the operation and had preserved samples were screened. Patients who received at least 2 cycles of treatment were selected. A total of 53 pairs of samples and the medical data of these patients were collected. The sensitivity of methotrexate (MTX) was classified by the treatment response: MTX sensitive group (complete response or partial response, and the maintenance time ≥ 4 weeks), and MTX-resistant group (stable disease or progressive disease) according to the Response Evaluation Criteria in Solid Tumors (RECIST) [21].


Fresh tissue samples were fixed with 4% neutral paraformaldehyde buffer, routinely dehydrated, waxed, embedded, and then cut into a 4-µm slide. Paraffin sections were baked, dewaxed, and antigen repaired with sodium citrate for 20 min. The slide was sealed with bovine serum albumin (BSA), and OGT antibody (#PA5-83154, Invitrogen, Foster, USA) diluted at 1:50 was dropped and incubated overnight in a 4 °C refrigerator. Rabbit secondary antibody (#PA5-59497) was added and incubated at room temperature for 1 h. Diaminobenzidine (DAB) was used to develop color, hematoxylin was used for counterstain. The slide was sealed with neutral gum and observed and photographed under an optical microscope (× 400 magnification, Olympus, Japan).

Induction of MTX-resistant osteosarcoma cells

Human osteosarcoma cell lines MNNG/HOS, Saos-2, U2OS, and MG63 were purchased from cell bank of the Chinese Academy of Sciences, Shanghai. MTX-resistant cells were inducted through gradually increasing MTX concentration. Briefly, MTX with an initial concentration of 50 ng/ml was added to the logarithmic proliferating osteosarcoma cells. After treated for 24 h, cells were washed and cultured with McCoy's 5A (#16600082, Thermo Fisher Scientific, USA) or DMEM medium (#10566016, Thermo Fisher Scientific) without drugs to restore the normal proliferation state and repeat the same concentration until the cells can adapt to this concentration and escalated to the next concentration. The increased concentrations of MTX were 50 ng/ml, 100 ng/ml, 200 ng/ml, 500 ng/ml, 1000 ng/ml, and 2000 ng/ml. The stably growing MTX (2000 ng/ml)-resistant cell lines were termed MTXr, such as MNNG/HOS-MTXr.

Measurement of MTX sensitivity in cells

The logarithmic growth cells were digested with 0.25% trypsin, resuspended, and prepared into single-cell suspension with culture medium. A total of 200 µL cells with a concentration of 5 × 104/ml were inoculated into 96 well plates and cultured for 24 h. MTX with gradient concentrations (including one control well) were prepared according to the pre-experiment results. Three multiple wells were set for each concentration. After cultured for 24 h, 20 µL MTT (#C0009M, Beyotime, China) were added to each well and cultured for another 4 h. The optical density (OD) of each well at 490 nm was detected by an enzyme labeling instrument to calculate the cell viability. Cell viability (%) = OD of experimental group / OD of control group × 100%. Detailed MTX concentration in osteosarcoma cells was: MNNG/HOS: 0, 1, 3, 4, 5 µM; MNNG/HOS-MTXr: 0, 2, 5, 10, 20 µM; MG63: 0, 1, 2, 3, 4 µM; MG63-MTXr: 0, 2, 5, 10, 20 µM.

Real-time quantitative PCR (RT-qPCR)

RT-qPCR was used to detect the expression of EBLN3P, miR-200a-3p, and OGT. Briefly, total RNA was extracted from preserved tissues or cultured cells by a commercial kit (#R1200, Solarbio Life Science, Beijing, China). The extracted RNA was reverse transcribed into cDNA. PCR was carried out according to the instructions of the quantitative kit (#RR037A, TaKaRa Bio Inc., Japan) and performed on 7500 Real-Time PCR Systems (Applied Biosystems, Foster, CA, USA). GADPH and U6 were used as internal references, respectively. The relative expression was calculated by the 2(− δ C(T)) method prescribed previously [10]. Used primers are listed in Table 1.

Table 1 Primers for RT-qPCR

RNA Oligoribonucleotides and cell transfection

The transfection was performed using Lipofectamine™ 2000, according to the manufacturer’s instructions. Cells were plated in a 24-well plate with a density of 1 × 104 cells per well. The transfected vector includes lentivirus harboring an shRNA of LncRNA EBLN3P (EBLN3P-si), miR-200a-3p mimic, OGT overexpression (OGT OE, #HG17892-NM, Sino Biological Inc., Beijing, China), OGT interfering sequence (OGT-si) and related negative control sequence, respectively. The EBLN3P-si and overexpression of EBLN3P (pcDNA-EBLN3P, EBLN3P OE) were used to knock down and upregulation of LncRNA EBLN3P expression, respectively. The other related sequences used in this study were synthesized by RIBOBIO Co., Ltd. The transfection efficiency was determined by qRT-PCR methods [22].

Wound healing assay measurement of migration and invasion

Cells with the concentration of 5 × 104 were inoculated into a 6-well plate. After cells were adherent to the culture medium, a sterilized two-milliliter pipette tip was used to generate wounding across the cell monolayer. The debris of cells was washed with PBS twice, and the cells were observed under a microscope at 0, 24, and 48 h after treatment. The scratch width was calculated by ImageJ software (NIH, USA).

Matrigel invasion assay

Two hundred microliter cells with a concentration of 5 × 104 cells/ml were inoculated into the Transwell chamber (Corning, CA, USA) and then moved into a 24-well plate. Six milliliter mediums containing 10% fetal bovine serum were added to the lower compartment of the chamber and cells were cultured for 48 h. The cells on the compartment of the chamber were fixed with formaldehyde for 20 min and stained with 0.5% crystal violet 20 min. The number of invaded cells in 5 separate fields was counted under the microscope, and the average value was calculated to evaluate the number of invaded cells.

Microarray data collection

The microarray data used in this study were retrieved from the GEO database ( The following strategy was used: (Homo sapiens) and (Expression profiling by array) and (Osteosarcoma). We selected data according to the following criteria: Each data set included osteosarcoma sample and normal sample and contained more than 20 samples. GSE65071 data set (plasma from 20 osteosarcoma samples, including 10 localized, 10 metastatic, and 15 healthy controls) and GSE79181 (all came from pre-treatment osteosarcoma tissues, including 9 patients without relapse and 14 patients with relapse) were analyzed. No ethical approval or informed consent was required in this study due to the public availability of data in the GEO databases.

Dual-luciferase reporter assay

Bioinformatic analysis was performed in online tools TargetScan ( and ENCORI ( To verify the binding results in cells, regions containing the binding site of the 3′-UTR targeted were generated by PCR amplification and subcloned into the NheI/XbaI sites of the pmirGLO dual-luciferase reporter plasmid (Cat.#E1330, Promega). The mutations were generated according to the manufacturer’s instructions. The miR-200a-3p mimic or negative control sequence, pmirGLO-EBLN3P-3′UTR plasmid, or mutation plasmid was co-transfected into cells using Lipofectamine 2000 Reagent (Invitrogen). Luciferase activities were measured 48 h after transfection. Firefly luciferase (luc2) was used as the primary reporter to monitor miRNA regulation and Renilla luciferase (hRluc-neo) acted as a control reporter for normalization and selection. Reduced firefly luciferase expression indicates the binding of introduced miRNAs to the cloned mRNA target sequence.

Western blot analysis

The total protein of tissue homogenate or cells was extracted from lysed cells by a protein extraction kit (#R0010, Solarbio Life Science). A total of 40 µg protein was loaded into each well after protein concentration was determined by BCA test kit (#PC0020, Solarbio Life Science). The proteins were separated by SDS-PAGE and transferred to a nitrocellulose membrane. After being blocked with 5% skimmed milk at room temperature for 120 min, the primary antibody for OGT (#ab250489), Vimentin (#ab52942), N-cadherin (#ab207608), E-cadherin (#ab239883), and GAPDH (#ab9484, all 1:1000 dilution) was added and incubated overnight at 4 °C, then incubated with secondary antibody, and cultured at room temperature for 2 h. ECL developer was loaded, and the expression of the proteins was analyzed by Image J software (NIH, USA).

Statistical analysis

All data were analyzed by GraphPad 8.0 software. The measurement data were expressed as mean ± standard deviation (SD). The t-test or nonparametric test was used to compare values between groups. One-way ANOVA was used for multi-group comparison, and LSD-t-test was used for pairwise comparison. P < 0.05 was considered statistically significant.


Highly expressed LncRNA EBLN3P is related to methotrexate (MTX) resistance in osteosarcoma

The qRT-PCR analysis showed expression of EBLN3P in osteosarcoma was significantly higher than in adjacent tissues (Fig. 1A). EBLN3P expression in MTX sensitive osteosarcoma tumor tissues was significantly lower than resistant osteosarcoma tumor tissues. (Fig. 1B). Four osteosarcoma cell lines were used in the current study, and their IC50 of MTX were screened by MTT assay, which is shown in Fig. 1C. Meanwhile, the expression of LncRNA EBLN3P in each cell lines was also examined by qRT-PCR, which is shown in Fig. 1D. As shown in Fig. 1E, the MTX IC50 for each cell lines were positively correlated with RNA expression of EBLN3P, also suggest that higher expression of EBLN3P might be associated with MTX resistance. We selected 2 osteosarcoma cell lines, MNNG/HOS and MG63 with relatively lower EBLN3P expressions used for further MTX-resistant cell establishment, and finally, the MTX-resistant cell line was successfully established, and IC50 increased from 0.2 to 1 µM for MNNG/HOS, and 0.35 to 1.3 µM for MG63 (Fig. 1F, G).

Fig. 1
figure 1

Increased EBLN3P is related to methotrexate (MTX) resistance in osteosarcoma. A Expression of EBLN3P in osteosarcoma was significantly higher than adjacent tissues. B Patients with MTX resistance have higher EBLN3P expression than MTX-sensitive patients. C The IC50 of MTX in human osteosarcoma cell line MNNG/HOS, Saos-2, U2OS, and MG63. D The expression of LncRNA EBLN3P in each cell lines. E The positive correlation of MTX IC50 and RNA expression of EBLN3P in each cell lines. F and G The cell viability of the 2-cell line increased with the induction of MTX, which demonstrated by increased IC50 after induction. **p < 0.01

Knockdown EBLN3P increases MTX sensitivity of osteosarcoma cells

Three complementary sequences were designed and synthesized to knock down EBLN3P (EBLN3P-si) in osteosarcoma cells, and the sequence with the highest interfering efficacy was used. The interfering efficacy was measured (Fig. 2A). The results showed that knockdown of EBLN3P decreased the migration ability of MNNG/HOS-MTXr and MG63-MTXr cells significantly (Fig. 2B, C). Furthermore, the cell viability of the 2 MTX-resistant cell line decreased after the knocking down of EBLN3P (Fig. 2D, E).

Fig. 2
figure 2

Knockdown EBLN3P increases MTX sensitivity of osteosarcoma cells. A qRT-PCR comparison of EBLN3P expression after knocking down EBLN3P (EBLN3P-si). **p < 0.01. B and C The migration ability of MNNG/HOS-MTXr and MG63-MTXr cells decreased significantly after knockdown of EBLN3P. **p < 0.01. D and E Comparison of cell viability of the 2 MTX-resistant cell line after knocking down of EBLN3P. * < 0.05, **p < 0.01, compared with EBLN3P NC at the same dose

EBLN3P sponge miR-200a-3p in MTX resistance osteosarcoma

Bioinformatics analysis demonstrated that miR-200a-3p is one of the potential binding genes with EBLN3P (Fig. 3A). qPCR results demonstrated that the expression of miR-200a-3p was significantly decreased in the osteosarcoma tissue compared to the adjacent tissues (Fig. 3B). To verify our results, the expression of miR-200a-3p in osteosarcoma patients from public data, GSE65071 data set and GSE79181 were analyzed. The results demonstrated expression of miR-200a-3p decreased significantly in osteosarcomas patients and patients with relapse (Fig. 3C). Moreover, knockdown of EBLN3P increased miR-200a-3p expression in MTX-resistant osteosarcoma cells (Fig. 3D). In vitro dual-luciferase activity analysis also demonstrated the luciferase activity of the WT- EBLN3P (wide type EBLN3P) sequence was suppressed by the mimic of miR-200a-3p (miR-200a-3p mimic) while the activity did not affect the Mut- EBLN3P (mutant EBLN3P sequence) sequence compared with negative control sequence of miR-200a-3p (miR-200a-3p NC) (Fig. 3E).

Fig. 3
figure 3

EBLN3P sponge miR-200a-3p in MTX resistance osteosarcoma. A The potential binding site of miR-200a-3p with EBLN3P. B The expression of miR-200a-3p was significantly decreased in the osteosarcoma tissue compared to adjacent tissues. C The expression of miR-200a-3p in osteosarcoma patients of GSE65071 and GSE79181 demonstrated significantly downregulated expression in osteosarcoma patients or osteosarcoma patients with relapse. D Knockdown of EBLN3P (EBLN3P-si) increased miR-200a-3p expression in MTX-resistant osteosarcoma cell line MNNG/HOS-MTXr and MG63-MTXr. E The luciferase activity of the WT- EBLN3P sequence was suppressed by the mimic of miR-200a-3p while the activity did not affect the Mut- EBLN3P (mutant EBLN3P) sequence compared with negative control sequence of miR-200a-3p (miR-200a-3p NC). **p < 0.01

miR-200a-3p increases MTX sensitivity via O-GlcNAc transferase (OGT)

To further identify the targeted regulating gene of miR-200a-3p, the TargetScan was used to search for the target gene. The O-GlcNAc transferase (OGT), a proven elevated protein in many tumors [23], was found to be one of the binding genes of miR-200a-3p (Fig. 4A). The immunohistochemistry and western analysis showed that increased expression of OGT is obvious when compared to adjacent tissues (Fig. 4B, C). We observed a negative relationship between miR-200a-3p and OGT expression levels in osteosarcoma tissues (Fig. 4D). Furthermore, upregulation of OGT almost counteracts the sensitization effect of miR-200a-3p to MTX in osteosarcoma cells (Fig. 4E, F).

Fig. 4
figure 4

miR-200a-3p increases MTX sensitivity via OGT. A The binding site between miR-200a-3p and OGT. B and C The immunohistochemistry and western analysis of OGT expression in osteosarcoma and adjacent tissues. D The correlation analysis showed that the expression of OGT is negatively correlated with the miR-200a-3p. E and F Upregulation of miR-200a-3p increases the sensitization of osteosarcoma cells to MTX, while upregulation of OGT almost abolished the effect of miR-200a-3p. *p < 0.05, **p < 0.01, compared with miR-200a-3p mimic; #p < 0.05, ##p < 0.01, compared with miR-200a-3p mimic + OGT OE

EBLN3P increased MTX-resistance of osteosarcoma cells by enhancing miR-200a-3p/OGT axis

The effect of EBLN3P on osteosarcoma cells was measured through transfecting pcDNA- EBLN3P to osteosarcoma MNNG/HOS and MG63 cells. We observed that upregulation of EBLN3P increased the invasion and MTX resistance of both cells (Fig. 5A–D). Blocking OGT by OGT-si crippled the resistant effect of EBLN3P on MTX in osteosarcoma cells (Fig. 5C, D). Moreover, miR-200a-3p could not increase the sensitivity of MTX in osteosarcoma cells after OGT was overexpressed (Fig. 5C, D). Research has demonstrated that OGT promotes EMT and invasion of cancer cells mainly via the junction with skeleton proteins [24,25,26]. We further detected the expression of invasion markers of EMT under the influence of loss-of-function of OGT. Western blot analysis also demonstrated that upregulation of EBLN3P increased protein expression of OGT, Vimentin, and N-cadherin while the expression of E-cadherin decreased (Fig. 5E).

Fig. 5
figure 5

EBLN3P increased MTX-resistance of osteosarcoma cells by enhancing miR-200a-3p/OGT axis. A Upregulation of EBLN3P (EBLN3P OE) increased the invasion ability of MNNG/HOS and MG63 cells. B Column comparison of invaded cells in MNNG/HOS and MG63 cells. **p < 0.01. C and D Effect of blocking OGT on EBLN3P upregulated cells was evaluated. Blocking OGT by OGT-si decreased the resistant effect of EBLN3P on MTX in MNNG/HOS and MG63 cells. *p < 0.05, **p < 0.01, compared with vector; #p < 0.05, ##p < 0.01, compared with EBLN3P OE + OGT-si. E Western blot analyzing the effect of EBLN3P, miR-200a-3p, and OGT on MNNG/HOS cells. Upregulation of EBLN3P leads to increased expression of OGT, Vimentin, and N-cadherin while decreasing the expression of E-cadherin. Further upregulation of miR-200a-3p or knockdown OGT (OGT-si) decreased the expression of OGT, Vimentin, and N-cadherin. *p < 0.05, **p < 0.01, compared with vector; #p < 0.05, ##p < 0.01, compared with EBLN3P OE


Osteosarcoma is highly malignant and tends to metastasize at the early stage. For patients with osteosarcoma, most current treatment programs usually involve intensive neoadjuvant chemotherapy [27]. However, chemotherapy resistance is an important factor affecting the disease progression and survival of patients with osteosarcoma [28]. Our study shows that EBLN3P can affect the MTX resistance of osteosarcoma cells by regulating miR-200a-3p, an important microRNA that affects EMT. Reversing chemotherapy resistance is of great significance to cancer management strategies. The mechanisms of drug resistance in tumors are complex and diverse, which are mainly related to ABC transporter, Topo II, oncogenes Bcl-2, p53, NF-κB, and other factors [29, 30]. MDR1 (P-gp) is one of the main causes of drug resistance in a variety of tumors, including osteosarcoma [31, 32]. It can expel chemotherapeutic drugs out of cells against the concentration gradient, thereby increasing the tolerance of tumor cells to chemotherapeutic drugs. The ability of tumors to reduce the uptake of drug molecules is also considered to be a mechanism of drug resistance. This mechanism is similar to the way of increasing drug efflux (or drug transport to intracellular concentration decreased), which will reduce the concentration of drug molecules in the cellular environment, thereby limiting its efficacy on target tumors. The chemotherapy drug molecules are most susceptible to this resistance mechanism, such as MTX, 5-fluorouracil, 8-azaguanine, and cisplatin. They have been proved to use transporters such as solute carriers (SLCs) to enter cells [33, 34].

MTX is an anti-folate cellular immunosuppressant. It mainly inhibits purine synthesis by blocking dihydrofolate reductase, thus inhibiting thymine synthesis, reducing the chemotaxis of neutrophils, and inhibiting the release of inflammatory cytokines. In osteosarcoma cells, MTX binds to the dihydrofolate reductase (DHFR) enzyme, leading to inhibition of DNA synthesis and replication and ultimately causing apoptosis [33], which induced apoptosis or programmed cell death due to induced DNA damage [35]. As a response, some cell lines can form clonal expansion colonies, express DNA repair enzymes, and produce strong resistance to MTX, cisplatin, and other drugs. The drug resistance of MTX in osteosarcoma is related to the above factors in different ways, we speculate there might be genes involved in MTX transport and metabolism might be regulated by genes of the miR-200a-3p/OGT axis.

EBLN3P is a novel lncRNA, reports show that expression of EBLN3P is upregulated in a variety of tumors. Existing research evidence shows that it plays as oncogene in tumor cells through different miRNA/target gene axes [10, 13]. Our study shows that EBLN3P promotes the metastasis and drug resistance of osteosarcoma through miR-200a-3p/OGT axis. The miRNAs also are closely related to the occurrence of bone-related diseases [8, 9, 36]. Acted as effective regulatory genes, miRNAs have been demonstrated binding their target gene such as ERK 2, TGF- β1, and Smad3 to regulate tendon homeostasis [36]. miR-200a-3p is a member of the miR-200 family, and this family controls EMT process both in development and tumorigenesis [37]. miR-200 also mediates the chemoresistance of tumor cells. Pan et al. [38] reported that miR-200b, a member of the miR-200 family regulates autophagy associated with chemoresistance in human lung adenocarcinoma. This is consistent with Jin’s report [39] that downregulated miR-200b is involved in the chemoresistance process. Same as miR-200b, miR-200a also act as a tumor-suppressive gene in hepatocellular carcinoma [40]. Our study showed that miR-200a-3p was downregulated in osteosarcoma, resulting in the downregulation of downstream target gene OGT, which in turn promoted the EMT process of osteosarcoma cells, leading to the metastasis and MTX resistance of osteosarcoma cell.

EMT is a key mechanism of tissue remodeling during multicellular morphogenesis. The role of EMT in drug resistance has been confirmed [41]. Loss of the epithelial phenotype is the most significant feature of EMT; meanwhile, the cells obtained the interstitial phenotype, and the expressions of vimentin, N-cadherin increased. O-GlcNAc glycosylation modification is widely involved in many cells’ life processes, such as cell cycle, gene transcription, protein translation and processing, and signal transduction. Abnormal modification of O-GlcNAc glycosylation occurred in many diseases [42]. Studies have confirmed that OGT is highly expressed in liver tumors [43], colon tumors, bladder cancer, and breast cancers [44].

Inconsistent with the reported results of Sombutthaweesri et al. [45], our study showed that the expression of OGT was increased in surgical resection samples of osteosarcoma, and the upregulated expression was associated with MTX resistance. We speculate that inconsistency of OGT expression might be related to different clinical characteristics of patients (stage, age, course of disease, and background treatment). Moreover, our study also showed that decreasing the expression of OGT in osteosarcoma cells could significantly enhance the sensitivity of osteosarcoma cells to MTX. In tumor cells, O-GlcNAc proteins could promote vimentin expression and cell migration [46], while OGT catalyzes O-GlcNAc proteins to exert its effect [47]. Chen et al. demonstrated that OGT could restrain the expansion of DNA damage [48]. Moreover, Barkovskaya et al. [49] reported that inhibition of OGT could activate the expression of the tumor-suppressor gene in tamoxifen-resistant breast cancer cells. These results indicate that OGT contributes the drug resistance in some tumors. By repairing damaged DNA and increasing the EMT of tumor cells, tumor cells strengthen the discharge of these chemotherapy drugs and increase the resistance to these agents.

Our study has several limitations. First, the number of tumor samples involved in the current study is small, and tumor tissues from more patients are further needed to confirm the relationship between LncEBLN3P with clinical parameters. Second, the major hypothesis in the present research has not been verified by in vivo research, so the current conclusion is preliminary. In vivo investigation and deeper clinical research are necessary to explore how does LncEBLN3P/miR-200a-3p/OGT axis regulate the MTX sensitivity in osteosarcoma.

In conclusion, our research demonstrated that LncEBLN3P is upregulated in osteosarcoma and affects the MTX resistance in osteosarcoma cells through downregulating miR-200a-3p, which in turn promoted the EMT process of osteosarcoma cells by increased OGT, and this needs to be further confirmed in the follow-up study.


  1. Zhao X, Wu Q, Gong X, Liu J, Ma Y. Osteosarcoma: a review of current and future therapeutic approaches. Biomed Eng Online. 2021;20(1):24.

    Article  Google Scholar 

  2. Cheng C, Ding Q, Zhang Z, Wang S, Zhong B, Huang X, Shao Z. PTBP1 modulates osteosarcoma chemoresistance to cisplatin by regulating the expression of the copper transporter SLC31A1. J Cell Mol Med. 2020;24(9):5274–89.

    Article  CAS  Google Scholar 

  3. Harris MA, Hawkins CJ. Recent and ongoing research into metastatic osteosarcoma treatments. Int J Mol Sci. 2022;23(7):3817.

    Article  CAS  Google Scholar 

  4. Marchandet L, Lallier M, Charrier C, Baud’huin M, Ory B, Lamoureux F. Mechanisms of resistance to conventional therapies for osteosarcoma. Cancers (Basel). 2021;13(4):683.

    Article  CAS  Google Scholar 

  5. Sun L, Li Y, Li H, Zhang J, Li B, Ye Z. Analysis of chemotherapy dosage and dosage intensity and survival outcomes of high-grade osteosarcoma patients younger than 40 years. Clin Ther. 2014;36(4):567–78.

    Article  CAS  Google Scholar 

  6. Zhang HD, Jiang LH, Zhong SL, Li J, Sun DW, Hou JC, Wang DD, Zhou SY, Tang JH. The role of long non-coding RNAs in drug resistance of cancer. Clin Genet. 2021;99(1):84–92.

    Article  CAS  Google Scholar 

  7. Giordano L, Porta GD, Peretti GM, Maffulli N. Therapeutic potential of microRNA in tendon injuries. Br Med Bull. 2020;133(1):79–94.

    Article  CAS  Google Scholar 

  8. Oliviero A, Della Porta G, Peretti GM, Maffulli N. MicroRNA in osteoarthritis: physiopathology, diagnosis and therapeutic challenge. Br Med Bull. 2019;130(1):137–47.

    Article  CAS  Google Scholar 

  9. Gargano G, Oliva F, Oliviero A, Maffulli N. Small interfering RNAs in the management of human rheumatoid arthritis. Br Med Bull. 2022;142(1):34–43.

    Article  Google Scholar 

  10. Xu XH, Song W, Li JH, Huang ZQ, Liu YF, Bao Q, Shen ZW. Long non-coding RNA EBLN3P regulates UHMK1 expression by sponging miR-323a-3p and promotes colorectal cancer progression. Front Med (Lausanne). 2021;8:651600.

    Article  Google Scholar 

  11. Mathias C, Muzzi JCD, Antunes BB, Gradia DF, Castro MAA, de Oliveira JC. Unraveling immune-related lncRNAs in breast cancer molecular subtypes. Front Oncol. 2021;11:692170.

    Article  Google Scholar 

  12. Li H, Wang M, Zhou H, Lu S, Zhang B. Long noncoding RNA EBLN3P promotes the progression of liver cancer via alteration of microRNA-144-3p/DOCK4 signal. Cancer Manag Res. 2020;12:9339–49.

    Article  CAS  Google Scholar 

  13. Dai S, Li N, Zhou M, Yuan Y, Yue D, Li T, Zhang X. LncRNA EBLN3P promotes the progression of osteosarcoma through modifying the miR-224-5p/Rab10 signaling axis. Sci Rep. 2021;11(1):1992.

    Article  CAS  Google Scholar 

  14. Itkonen HM, Loda M, Mills IG. O-GlcNAc transferase: an auxiliary factor or a full-blown oncogene? Mol Cancer Res. 2021;19(4):555–64.

    Article  CAS  Google Scholar 

  15. Zhu K, Wang T, Chen H, Chen T, Zhong Z. Expression of O-GlcNAc transferase in osteosarcoma and its effect on proliferation and sensitivity to cisplatin. China Oncol. 2020;30(9):682–8.

    Article  CAS  Google Scholar 

  16. Deng X, Yi X, Huang D, Liu P, Chen L, Du Y, Hao L. ROCK2 mediates osteosarcoma progression and TRAIL resistance by modulating O-GlcNAc transferase degradation. Am J Cancer Res. 2020;10(3):781–98.

    CAS  Google Scholar 

  17. Dudas J, Ladanyi A, Ingruber J, Steinbichler TB, Riechelmann H. Epithelial to mesenchymal transition: a mechanism that fuels cancer radio/chemoresistance. Cells. 2020;9(2):428.

    Article  CAS  Google Scholar 

  18. Lu Y, Lu J, Li X, Zhu H, Fan X, Zhu S, Wang Y, Guo Q, Wang L, Huang Y, et al. MiR-200a inhibits epithelial-mesenchymal transition of pancreatic cancer stem cell. BMC Cancer. 2014;14:85.

    Article  CAS  Google Scholar 

  19. Pichler M, Ress AL, Winter E, Stiegelbauer V, Karbiener M, Schwarzenbacher D, Scheideler M, Ivan C, Jahn SW, Kiesslich T, et al. MiR-200a regulates epithelial to mesenchymal transition-related gene expression and determines prognosis in colorectal cancer patients. Br J Cancer. 2014;110(6):1614–21.

    Article  CAS  Google Scholar 

  20. Tao X, Cheng J, Wang X. Effect of miRNA-200a on radiosensitivity of osteosarcoma cells by targeting Bone morphogenetic protein receptor 2. Bioengineered. 2021;12(2):12625–35.

    Article  CAS  Google Scholar 

  21. Guenther LM, Rowe RG, Acharya PT, Swenson DW, Meyer SC, Clinton CM, Guo D, Sridharan M, London WB, Grier HE, et al. Response evaluation criteria in solid tumors (RECIST) following neoadjuvant chemotherapy in osteosarcoma. Pediatr Blood Cancer. 2018;65(4):26896.

    Article  CAS  Google Scholar 

  22. Mauro C, Pacifico F, Lavorgna A, Mellone S, Iannetti A, Acquaviva R, Formisano S, Vito P, Leonardi A. ABIN-1 binds to NEMO/IKKgamma and co-operates with A20 in inhibiting NF-kappaB. J Biol Chem. 2006;281(27):18482–8.

    Article  CAS  Google Scholar 

  23. Akella NM, Le Minh G, Ciraku L, Mukherjee A, Bacigalupa ZA, Mukhopadhyay D, Sodi VL, Reginato MJ. O-GlcNAc transferase regulates cancer stem-like potential of breast cancer cells. Mol Cancer Res. 2020;18(4):585–98.

    Article  CAS  Google Scholar 

  24. Ge X, Peng X, Li M, Ji F, Chen J, Zhang D. OGT regulated O-GlcNacylation promotes migration and invasion by activating IL-6/STAT3 signaling in NSCLC cells. Pathol Res Pract. 2021;225:153580.

    Article  CAS  Google Scholar 

  25. Gao J, Yang Y, Qiu R, Zhang K, Teng X, Liu R, Wang Y. Proteomic analysis of the OGT interactome: novel links to epithelial-mesenchymal transition and metastasis of cervical cancer. Carcinogenesis. 2018;39(10):1222–34.

    Article  CAS  Google Scholar 

  26. Li X, Wu Z, He J, Jin Y, Chu C, Cao Y, Gu F, Wang H, Hou C, Liu X, et al. OGT regulated O-GlcNAcylation promotes papillary thyroid cancer malignancy via activating YAP. Oncogene. 2021;40(30):4859–71.

    Article  CAS  Google Scholar 

  27. Benjamin RS. Adjuvant and neoadjuvant chemotherapy for osteosarcoma: a historical perspective. Adv Exp Med Biol. 2020;1257:1–10.

    Article  CAS  Google Scholar 

  28. Alemi F, Malakoti F, Vaghari-Tabari M, Soleimanpour J, Shabestani N, Sadigh AR, Khelghati N, Asemi Z, Ahmadi Y, Yousefi B. DNA damage response signaling pathways as important targets for combination therapy and chemotherapy sensitization in osteosarcoma. J Cell Physiol. 2022.

    Article  Google Scholar 

  29. Zheng HC. The molecular mechanisms of chemoresistance in cancers. Oncotarget. 2017;8(35):59950–64.

    Article  Google Scholar 

  30. Ozes AR, Miller DF, Ozes ON, Fang F, Liu Y, Matei D, Huang T, Nephew KP. NF-kappaB-HOTAIR axis links DNA damage response, chemoresistance and cellular senescence in ovarian cancer. Oncogene. 2016;35(41):5350–61.

    Article  CAS  Google Scholar 

  31. Lee YH, Yang HW, Yang LC, Lu MY, Tsai LL, Yang SF, Huang YF, Chou MY, Yu CC, Hu FW. DHFR and MDR1 upregulation is associated with chemoresistance in osteosarcoma stem-like cells. Oncol Lett. 2017;14(1):171–9.

    Article  CAS  Google Scholar 

  32. Qin T, Zhu W, Kan X, Li L, Wu D. Luteolin attenuates the chemoresistance of osteosarcoma through inhibiting the PTN/beta-catenin/MDR1 signaling axis by upregulating miR-384. J Bone Oncol. 2022;34:100429.

    Article  Google Scholar 

  33. Hattinger CM, Patrizio MP, Fantoni L, Casotti C, Riganti C, Serra M. Drug Resistance in osteosarcoma: emerging biomarkers, therapeutic targets and treatment strategies. Cancers (Basel). 2021;13(12):2878.

    Article  CAS  Google Scholar 

  34. Li Q, Shu Y. Role of solute carriers in response to anticancer drugs. Mol Cell Ther. 2014;2:15.

    Article  Google Scholar 

  35. Holohan C, Van Schaeybroeck S, Longley DB, Johnston PG. Cancer drug resistance: an evolving paradigm. Nat Rev Cancer. 2013;13(10):714–26.

    Article  CAS  Google Scholar 

  36. Gargano G, Oliviero A, Oliva F, Maffulli N. Small interfering RNAs in tendon homeostasis. Br Med Bull. 2021;138(1):58–67.

    Article  CAS  Google Scholar 

  37. Hill L, Browne G, Tulchinsky E. ZEB/miR-200 feedback loop: at the crossroads of signal transduction in cancer. Int J Cancer. 2013;132(4):745–54.

    Article  CAS  Google Scholar 

  38. Pan B, Feng B, Chen Y, Huang G, Wang R, Chen L, Song H. MiR-200b regulates autophagy associated with chemoresistance in human lung adenocarcinoma. Oncotarget. 2015;6(32):32805–20.

    Article  Google Scholar 

  39. Jin HF, Wang JF, Song TT, Zhang J, Wang L. MiR-200b inhibits tumor growth and chemoresistance via targeting p70S6K1 in lung cancer. Front Oncol. 2020;10:643.

    Article  Google Scholar 

  40. Feng J, Wang J, Chen M, Chen G, Wu Z, Ying L, Zhuo Q, Zhang J, Wang W. miR-200a suppresses cell growth and migration by targeting MACC1 and predicts prognosis in hepatocellular carcinoma. Oncol Rep. 2015;33(2):713–20.

    Article  CAS  Google Scholar 

  41. Du B, Shim JS. Targeting epithelial-mesenchymal transition (EMT) to overcome drug resistance in cancer. Molecules. 2016;21(7):965.

    Article  CAS  Google Scholar 

  42. Burt RA, Alghusen IM, John Ephrame S, Villar MT, Artigues A, Slawson C. Mapping the O-GlcNAc modified proteome: applications for health and disease. Front Mol Biosci. 2022;9:524.

    Article  CAS  Google Scholar 

  43. Xu W, Zhang X, Wu JL, Fu L, Liu K, Liu D, Chen GG, Lai PB, Wong N, Yu J. O-GlcNAc transferase promotes fatty liver-associated liver cancer through inducing palmitic acid and activating endoplasmic reticulum stress. J Hepatol. 2017;67(2):310–20.

    Article  CAS  Google Scholar 

  44. Fardini Y, Dehennaut V, Lefebvre T, Issad T. O-GlcNAcylation: A new cancer hallmark? Front Endocrinol (Lausanne). 2013;4:99.

    Article  Google Scholar 

  45. Sombutthaweesri T, Wu S, Chamusri N, Settakorn J, Pruksakorn D, Chaiyawat P, Sastraruji T, Krisanaprakornkit S, Supanchart C. Relationship between O-GlcNAcase expression and prognosis of patients with osteosarcoma. Appl Immunohistochem Mol Morphol. 2022;30(1):e1–10.

    Article  CAS  Google Scholar 

  46. Komura K, Ise H, Akaike T. Dynamic behaviors of vimentin induced by interaction with GlcNAc molecules. Glycobiology. 2012;22(12):1741–59.

    Article  CAS  Google Scholar 

  47. Yang X, Qian K. Protein O-GlcNAcylation: emerging mechanisms and functions. Nat Rev Mol Cell Biol. 2017;18(7):452–65.

    Article  CAS  Google Scholar 

  48. Chen Q, Yu X. OGT restrains the expansion of DNA damage signaling. Nucleic Acids Res. 2016;44(19):9266–78.

    Article  CAS  Google Scholar 

  49. Barkovskaya A, Seip K, Prasmickaite L, Mills IG, Moestue SA, Itkonen HM. Inhibition of O-GlcNAc transferase activates tumor-suppressor gene expression in tamoxifen-resistant breast cancer cells. Sci Rep. 2020;10(1):16992.

    Article  CAS  Google Scholar 

Download references



Author information

Authors and Affiliations



BB was responsible for the conception and design of the study. MXS and HYA were responsible for the acquisition and analysis of data. YBS and YBS conducted bioinformatics and statistical analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bing Bai.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sun, MX., An, HY., Sun, YB. et al. LncRNA EBLN3P attributes methotrexate resistance in osteosarcoma cells through miR-200a-3p/O-GlcNAc transferase pathway. J Orthop Surg Res 17, 557 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: