922
Views
3
CrossRef citations to date
0
Altmetric
Research Paper

Long non-coding RNA long intergenic non-coding 00641 mediates cell progression with stimulating cisplatin-resistance in osteosarcoma cells via microRNA-320d/myeloid cell leukemia-1 axis

, , , , , & show all
Pages 7238-7252 | Received 26 Oct 2021, Accepted 10 Feb 2022, Published online: 10 Mar 2022

ABSTRACT

As a staple chemotherapy medicine, cisplatin (DDP) is extensively applied in cancer patients, but its drug resistance is limited. Numerous studies have elucidated that long non-coding RNA (lncRNA) performs as a pivotal agent in osteosarcoma (OS). Nevertheless, lncRNA long intergenic non-coding 00641 (LINC00641)’s functions in DDP resistance for OS remain obscure. The purpose of this study was to investigate the effect and mechanism of LINC00641 on drug resistance of OS. The tissues of both clinical cancer patients and the normal control were gathered. Detection of LINC00641, microRNA-320d (miR-320d) and myeloid cell leukemia-1 (MCL1) was conducted. After the selection of OS cell lines, the detection of cell advancement was applied. Series of experiments were conducted to verify the interaction of LINC00641, miR-320d and MCL1. Xenografted tumor model in vivo was utilized to determine the function of LINC00641. The data displayed, LINC00641 was prominently elevated in OS tissues and cells, especially in DDP-resistant tumors and cell lines. Knock-down LINC00641 was able to attenuate progression of DDP-resistant OS cells thus dampening their drug resistance toward DDP. Moreover, knock-downing LINC00641 gene was also able to manifest antagonism toward DDP-resistance in vivo. On the grounds of bioinformatics prediction, a direct binding of LINC00641 with miR-320d existed, whose target was MCL1. Meanwhile, LINC00641 modulated MCL1 via targeting miR-320d. Additionally, repressive LINC00641 blocked MCL1 via emulative interaction with miR-320d, thus expediting DDP-sensitivity of OS cells. All in all, it is found that LINC00641 is available to escalate drug resistance of DDP-resistant OS cells via mediation of miR-320d/MCL1 axis.

Graphical abstract

1 Introduction

Osteosarcoma (OS) is the most dominating malignant bone cancer among children and adolescents [Citation1]. With the progress of radiotherapy, adjuvant chemotherapy and surgery (whose scope of tumor resection is more far-ranging), the global survival rate of most OS patients has remained ascending notably [Citation2]. The universal treatment for OS patients is surgery combined with diverse medicine chemotherapy, like methotrexate, doxorubicin, cisplatin (DDP) and ifosfamide [Citation3]. Even though the combination chemotherapy has achieved elevated survival rate, unfortunately, OS is a relatively drug-resistant disease, and the incidence rate of local recurrence and distant metastasis of OS patients after operation still stays augmented [Citation3]. For the sake of exploring specific biomarkers and therapeutic targets, it is fairly indispensable to learn about the molecular mechanism of OS.

As a kind of non-coding RNA, long non-coding RNAs (lncRNAs) participate in miscellaneous physiological and pathological functions like cell progression, metabolism, differentiation [Citation4–7]. LncRNAs perform as a crucial agent in the process of OS. For instance, LncRNA MRPL23-AS1 activates Wnt/β-catenin signaling pathway by repressing microRNA (miR)-30b and upregulating myosin heavy chain 9 (MYH9) to promote the progress of OS and cancerization [Citation8]; CRISPR/cas9-mediated overexpression of lncRNA SRY-box transcription factor 21 antisense divergent transcript 1 modulates the proliferation of OS by increasing the expression of mechanistic target of rapamycin kinase and Kruppel like factor 4 [Citation9]. lncRNA OIP5-AS1 induces LPAATβ/PI3K/AKT/mTOR signaling pathway via sponging miR-340-5p, which contributes to DDP-resistance in OS [Citation10]. LncRNA OIP5-AS1 modulates adriamycin-resistance via the mediation of miR-137-3p/PTN axis in OS [Citation11]. LncRNA LINC00641 is located on chromosome 14q11.2, and it frequently descends under normal conditions [Citation12]. Latest studies have clarified that LINC00641 is implicated in the tumorigenicity of bladder cancer, non-small cell lung cancer (NSCLC) and acute myeloid leukemia [Citation13–15]. Nevertheless, the biological function of LINC00641 in OS and its resistance toward DDP have still stayed unclear.

MicroRNAs (miRNAs), a type of endogenous non-coding small RNA, demonstrate mediation of gene in the post-transcription level. A great many cellular pathways and functions, involving cell metabolism, differentiation and apoptosis, are under the modulation of miRNAs, hence their disorder will generate numerous human diseases [Citation16]. An ascending number of studies have uncovered miRNAs take a pivotal effect in chemotherapy-resistance of cancer, and have verified that miRNAs have underlying effect in determining drug sensitivity or drug resistance [Citation17]. With its abnormal manifestation in diverse kinds of tumors, it has been testified that miR-320d is momentous in NSCLC [Citation18], cardiac adenocarcinoma [Citation19]. Yet its function in DDP resistance of OS cells has remained still undefined.

Myeloid cell leukemia −1 (MCL1) belongs to Bcl-2 family [Citation20]. A great many literatures and studies have expounded that MCL1 is rather vital in the survival and death of miscellaneous cells [Citation21,Citation22]. It has been indicated that miR-26a reverses multidrug resistance of OS via targeting MCL1 [Citation23]. Yet it has remained unknown whether LINC00641 signifies the mediation of MCL1 in OS.

This study was to figure out how LINC00641 manifests in DDP-resistant tissues and cell lines in OS, and to explore the influence of LINC00641 on DDP-resistance in OS patients in vivo and in vitro. Apart from that, the author also discussed the latent molecular mechanism of LINC00641 in DDP-resistance. It was hypothesized that LINC00641 regulated cell proliferation, invasion and apoptosis and promoted DDP resistance of OS cells through miR-320d/MCL1 axis. This study was inclined to be conducive to providing an underlying therapeutic method for OS.

2 Materials and methods

2.1 Patients and specimens

The study gained authorization from the Ethics Committee of Huai’an Second People’s Hospital, and collection of written informed consent was from both Huai’an Second People’s Hospital and all patients. After the resection operation of patients in Huai’an Second People’s Hospital, both OS tissues and para-cancerous normal tissues were gathered from 58 OS patients. All cancer tissue specimens were diagnosed as OS via pathological examination. The quick refrigeration and preservation in liquid nitrogen of all fresh samples were exerted until further experiments. Then, 58 patients with OS, on the grounds of their sensitivity to chemotherapy medicine, were assigned into two groups. One was the chemotherapy sensitive (the tumor was relieved after six cycles of chemotherapy, n = 24), the other was the chemotherapy resistant (the tumor emerged stability or progression after six cycles of chemotherapy, n = 34). Additionally, 58 patients were assigned into two groups in line with LINC00641 expression, and calculation of the total survival period of all participants was in different stages (0, 20, 40, 60 months) after DDP therapy.

2.2 Cell culture and transfection

The culture of OS cell lines Saos2 and MG-63 (all American Type Culture Collection, Manassas, VA, USA) was exerted in elevated concentration of DDP (Sigma, St.Louis, MO, USA) for more than 6 months, and then DDP-resistant cell lines Saos2/DDP and MG-63/DDP were constructed (Attached ). Normal cervical epithelial cell line HaCaT (Shanghai Institute of Biochemistry and Cell Biology, Shanghai, China) was applied. Then all the cells were cultivated in Dulbecco’s Modified Eagle Medium (DMEM, Gibco, Carlsad, CA, CA) comprising 10% fetal bovine serum (FBS) (Gibco), 100 U/mL penicillin and 100 U/mL streptomycin (Sigma Aldrich, Shanghai, China).

Figure 1. LINC00641 is augmented in DDP-resistant OS tissues and cell lines, and its elevation heralds unpleasing prognosis. (a) RT-qPCR to detect LINC00641 in normal tissues and DDP-sensitive or drug-resistant cancer tissues; (b) Kaplan-Meier method to analyze the function of LINC00641 in OS prognosis; (c) RT-qPCR to examine LINC00641 in hFOB1.19 cells, OS cells and OS DDP-resistant cell lines. *P < 0.05.

Figure 1. LINC00641 is augmented in DDP-resistant OS tissues and cell lines, and its elevation heralds unpleasing prognosis. (a) RT-qPCR to detect LINC00641 in normal tissues and DDP-sensitive or drug-resistant cancer tissues; (b) Kaplan-Meier method to analyze the function of LINC00641 in OS prognosis; (c) RT-qPCR to examine LINC00641 in hFOB1.19 cells, OS cells and OS DDP-resistant cell lines. *P < 0.05.

Synthesis of shRNA short hairpin RNA targeting LINC00641 (sh-LINC00641) and shRNA interference control (sh-NC) contraposing LINC00641, pcDNA and pcDNA-MCL1 elevated vector (MCL1) was via Genepharma (Shanghai, China). MiR-320d mimic (miR-320d), mimic negative control (miR-NC), miR-320d inhibitor (anti-miR-320d) and inhibitor NC (anti-NC) (all RiboBio Co., Ltd., Guangzhou, China) were put into practice. When the cell concentration of Saos-2/DDP and MG-63/DDP attained 50 ~ 60%, TE transfection of miRNA mimic (10 nM) or vector was adopted via LipofectamineSaos2/DDP 2000 reagent (Carlsbad, California, USA). Ten cells were gathered for 48 h for subsequent analysis [Citation24].

2.3 Quantitative real time polymerase chain reaction (RT-qPCR)

Trizol reagent (Invitgen) was utilized to lyse treated cells and extract the total RNA. Reverse transcription into complementary DNA (cDNA) was via reverse transcription kit (Invitgen). RT-qPCR was implemented on the applied biological system 7300 (Invitgen) via SYBR Green RT-PCR kit (Invitgen). Detection of RNA expression was conducted. The calculation of relative expression amount of indicator genes was via 2−ΔΔCT method, and glyceraldehyde-3-phosphate dehydrogenase (GAPDH)/U6 was adopted as standardized internal references [Citation25].

2.4 Cell viability assay

The 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay was utilized to examine cell viability. Next, the seeding of transfected DDP-resistant cells was into 96-well plates at a density of 5 × 103 cells/well, and incubation together with DDP was conducted with different doses. In different time points (24, 48, 72 h), addition of 20 μL thiazole blue solution was to each well, and then addition of formaldehyde generated via dimethyl sulfoxide dissolving was manifested. At last, the measurement of absorbance was exerted at 490 nm via a microplate reader (Bio-Rad, Hercules, CA, USA).

2.5 Cell migration and invasion assay

After the detection of migration and invasion abilities in vitro of DDP-resistant cells, the seeding of transfected DDP-resistant cells was utilized into a 24-well plate upper chamber (Becton-Dickinson, Franklin Lake, New Jersey, USA) with or without matrix, and then the culture in serum-free DMEM was carried out. Then the lower chamber in each well was added with complete culture solution comprising 10% FBS, which was regarded as a chemical inducer. After incubation, inactive cells on the upper surface were removed from with a dry cotton swab. Then fixation of the cells was conducted via methanol with ventilating appropriately the chamber. Next, after staining with 0.1% crystal violet, non-transferring cells on the upper layer were lightly wiped off with a cotton swab, and were immediately observed in five visual fields under a phosphate buffer saline (PBS) 400-fold microscope for counting [Citation26].

2.6 Cell apoptosis analysis

After transfection, analysis of the apoptosis rate of DDP-resistant cells was via Annexin V- fluoresceinisothiocyanat/ propidium iodide apoptosis detection kit (Solarbio, Beijing). All specimens were in triplicate [Citation10].

2.7 Total RNA extraction and quantitative reverse transcriptase polymerase chain (RT-qPCR)

With extraction of total RNA, TRIzol reagent was obtained from Invitrogen Company, and applied following the manufacturer’s instructions. After cDNA synthesis, Power SYBR Green (TaKaRa, Shiga, Japan) was applied on ABI Step-One real-time system (Thermo Fisher, Waltham, MA, USA) for performing qPCR. The relative quantitative method (2−ΔΔCT) was used to calculate gene expression amount. U6 and GAPDH were applied as normalized genes. The primer sequences were listed in .

Table 1. Primer sequences

2.8 Western blot assay

The extraction of cell protein was conducted from 6-well plate via Radioimmunoprecipitation Assay lysis buffer solution consisting of protease inhibitor. After the quantification via bicinchoninic acid kit (Thermo Fisher Science), the cell protein was diluted to identical concentration in loading buffer and denatured at 95°C, then the electroblot of 20 μg protein was onto Polyvinylidene fluoride membrane via sodium dodecyl sulfate polyacrylamide gel electrophoresis. After seal with skim milk, the protein was co-cultured separately with anti- MCL1 (1: 1000, 5453), caspase-3 (1: 1000; 9661), PARP (1: 1000; 9542) and GAPDH (1: 1000, 2118), and all antibodies were from Cell Signaling Technology. Then the detection of the protein was utilized via chemiluminescence detection system (GE Healthcare, Chicago, IL, USA) [Citation27].

2.9 The luciferase activity assay

The cloning of LINC00641 mRNA consisting of miR-320d wild-type (WT) and mutant (MUT) binding sequences was separately exerted into psiCHECK™-2 luciferase plasmid (Promega, Shanghai, China). Saos2/DDP and MG-63/DDP cells were seeded in a 24-well plate, and then respectively co-transfected with miR-320d mimic/miR-NC and LINC00641-WT/MUT. Finally, analysis of the luciferase activity was carried out via dual luciferase assay kit (Promega) [Citation11].

2.10 RNA pull-down assay

The association between LINC00641 and miR-320d was evaluated by biotinylated RNA pull down assay. Biotin-labeled miR-361-3p probes and NC probes were provided by GenePharma. In short, 1 × 105 HOS cells were obtained, then lysed and treated with ultrasound. Probes were incubated at 25°C for 3 h using C-1 magnetic beads (Life Technologies). The cell lysates with miR-361-3p or oligo probes were incubated at 4°C for 24 h and then washed with PBS buffer. Binding RNA was purified with Trizol reagent and RT-qPCR was performed [Citation28].

2.11 RNA immunoprecipitation (RIP) assay

RIP assay was conducted using human Argonaute2 (Ago2) antibody (anti-Ago2) with the help of the Magna RIP Kit (EMD Millipore, Billerica, MA, USA). Magnetic beads were incubated with indicated antibodies. The cells were lysed in RIP lysis buffer. Cell lysates were obtained and magnetic beads were combined to incubate with indicated antibodies. RNA was purified from the precipitate and analyzed by RT-qPCR [Citation29].

2.12 Tumor xenograft test

Lentivirus vector (GenePharma, Shanghai, China), that is lentit-short hairpin LINC00641, named sh-LINC00641, was adopted to facilitate stabilization of LINC00641, and the control of messy code of shRNA (sh-NC) was put into use. Meanwhile, 6-week-old male Balb/c mice (five in each group) were collected from National Laboratory Animal Center (Beijing, China). This study had obtained ratification from Animal Research Committee, Huai’an Second People’s Hospital. All animal studies were executed on the grounds of Animal Research: Reporting of In Vivo Experiments Guidelines and the Basel Declaration. All animals were taken humane care in line with the guidelines of the National Institutes of Health. After steady transfection with sh-LINC00641 or sh-NC, subcutaneous injection of the cells (5 × 106) was utilized into the left abdomen of nude mice. After 7-d injection, the mice were given PBS (Invitrogen) or 3 mg/kg DDP (DDP) (Sigma) and their tumor volume was measured weekly. All mice were implemented euthanasia after 28 d, and their resected tumors were weighed (Attached ).

2.13 Terminal deoxynucleotidyl transferase-mediated dUTP-biotin nick end labeling (TUNEL)

Xenografted tumor sections were dewaxed with xylene and rehydrated with gradient alcohol. The 4% paraformaldehyde was applied for fixing the sections. Then the sections were treated with proteinase K solution (20 μg/mL) and TUNEL reaction buffer. Then the sections were utilized staining with diaminobenzidine and hematoxylin. Next the sections were dehydrated with gradient ethanol, treated with xylene, and sealed with neutral resin. The observation of TUNEL positive cells (brown granules) was conducted under the microscope.

2.14 Statistical analysis

All statistical analysis was adopted via GraphPad Prism 7(GraphPad Inc, San Diego, California, USA). Manifestation of the data was as mean ± standard deviation. The pertinence among LINC00641, miR-320d and MCL1 was analyzed via Pearson test. Kaplan-Meier survival curve was drawn and logarithmic rank test were put into use to analyze the survival distinctions between two groups. Student’s t-test or one-way analysis of variance (ANOVA) was exerted to evaluate the prominent distinctions between groups. P < 0.05 was considered the apparent difference.

3 Results

Here, it was aimed to investigate the role of LINC00641 and its downstream molecular mechanisms in OS cells. A series of in vitro experiments was conducted and it was found that LINC00641 regulated cell proliferation, invasion and apoptosis, and promoted DDP resistance in OS cells through miR-320d/MCL1 axis. Therefore, the function and mechanism of LINC00641 in OS were first investigated in the data, providing new insights into the pathogenesis of OS.

3.1 LINC00641 manifests elevation in DDP-resistant OS tissues and cancer cell lines, which portends unpleasing prognosis

For the sake of exploring the function of LINC00641 in DDP-resistant OS, detection of LINC00641 was conducted in DDP-sensitive and -resistant OS tissues and cells. The results clarified that LINC00641 was saliently augmented in OS tissues and cells, especially in DDP-resistant tumors and cell lines (Saos2/DDP and MG-63/DDP) (). Clinical correlation analysis manifested higher LINC00641 level was associated with higher Enneking stage, larger tumor size, and positive distant metastasis in patients with OS (P < 0.05, ). Additionally, 58 patients who had been under DDP therapy were executed 60-month follow-up. And it came out that the overall survival rate of patients with ascending LINC00641 was inferior to that of patients with depressive one (), which insinuated that the augmentation of LINC00641 was a predictor of unpleasing prognosis of OS patients.

Table 2. The relationship between LINC00641 expression level and clinical characteristics in patients with OS

3.2 Inhibition of LINC00641 restrains the proliferation, migration and invasion, but promotes apoptosis of OS cells

To investigate the role of LINC00641 in OS, LINC00641 was first down-regulated in two parent strains Saos2 and MG-63 to explore its role. Saos2 and MG-63 were transfected with sh-LINC00641 and its NC, and via qPCR was detected that the expression of LINC00641 was clearly decreased in cells (). Via MTT assay was detected cell proliferation, and it was found that the proliferation ability of Saos2 and MG-63 was apparently reduced after LINC00641 was inhibited (). Apoptosis level was detected by flow cytometry, and it was found that repression of LINC00641 promoted cell apoptosis (). In addition, apoptosis-related proteins Bax, Bcl-2, Cleaved caspase3 and PARP expression levels were detected by qPCR or Western blot, and the results manifested Bax, cleaved-PARP, cleaved-caspase3 expression was elevated, while Bcl-2 expression was decreased clearly after LINC00641 down-regulation (). Transwell was applied to detect the migration and invasion abilities of cells and it was found that low expression of LINC00641 effectively restrained the migration and invasion of Saos2 and MG-63 cells (). In summary, it was found that inhibition of LINC00641 constrained the proliferation, migration and invasion, but promoted apoptosis of OS cells.

Figure 2. Inhibition of LINC00641 represses the proliferation, migration and invasion, but promotes apoptosis of OS cells. (a) qPCR to detect the expression of LINC00641 in Saos2 and MG-63 cells; (B/C) MTT to detect the proliferation ability of Saos2 and MG-63 cells; (d) Flow cytometry to detect apoptosis of Saos2 and MG-63 cells; (E/F) qPCR to detect the mRNA expression of Bax and Bcl-2 in Saos2 and MG-63 cells. (G/H) Western blot to detect the protein expression of Caspase-3 and PARP in Saos2 and MG-63 cells; (I/J) Transwell to detect the migration and invasion of Saos2 and MG-63 cells. * vs. the sh-NC, P < 0.05.

Figure 2. Inhibition of LINC00641 represses the proliferation, migration and invasion, but promotes apoptosis of OS cells. (a) qPCR to detect the expression of LINC00641 in Saos2 and MG-63 cells; (B/C) MTT to detect the proliferation ability of Saos2 and MG-63 cells; (d) Flow cytometry to detect apoptosis of Saos2 and MG-63 cells; (E/F) qPCR to detect the mRNA expression of Bax and Bcl-2 in Saos2 and MG-63 cells. (G/H) Western blot to detect the protein expression of Caspase-3 and PARP in Saos2 and MG-63 cells; (I/J) Transwell to detect the migration and invasion of Saos2 and MG-63 cells. * vs. the sh-NC, P < 0.05.

3.3 Knock-down LINC00641 attenuates progression of DDP-resistant OS cells

In order to further investigate the biological effect of LINC00641 on DDP-resistant OS cells, sh-NC/LINC00641 was transfected into Saos2/DDP and MG-63/DDP cells, with validating the transfection efficiency. The results manifested that LINC00641 was declined in cells transfected with sh-LINC00641 (). After 48 h of DDP treatment with different concentrations, it was found that LINC00641 loss reduced the IC50 value of Saos2/DDP and MG-63/DDP cells, and repressed the proliferation ability of DDP resistant cells (). In addition, flow cytometry results showed that LINC00641 knockdown clearly induced apoptosis of drug-resistant cells (). Meanwhile, qPCR results manifested LINC00641 deletion promoted the expression of Bax in Saos2/DDP and MG-63/DDP cells, while restrained the expression of Bcl-2, further suggesting that LINC00641 deletion could promote apoptosis of drug-resistant cells (). Moreover, cleaved-PARP, cleaved-caspase3 results of Western blot detection showed that LINC00641 inhibition clearly increased cleaved-PARP, cleaved-caspase3 protein (). Transwell assay was then performed, and the results showed that LINC00641 knockdown could repress the migration and invasion abilities of Saos2/DDP and MG-63/DDP cells (). In conclusion, knockdown of LINC00641 reduced the resistance of DDP resistant cells to DDP by restraining the proliferation and inducing apoptosis of DDP resistant OS cells.

Figure 3. Knock-down of LINC00641 gene is able to restrain the survival and advancement of DDP-resistant OS cells. (a) After transfection of sh-LINC00641, detection of LINC00641; (b) IC50 values; (c-d) MTT assay to detect the survival rate and cell proliferation; (e-f) Flow cytometry to analyze the apoptosis rate; (g-h) RT-qPCR to examine Bcl-2 and Bax mRNA; (i-j) Western blotto detect caspase-3 and PARP protein expression; (k-l) Transwell method to analyze the migration and invasion abilities. In Saos2/DDP and MG-63/DDP cells. * vs. the sh-NC, P < 0.05.

Figure 3. Knock-down of LINC00641 gene is able to restrain the survival and advancement of DDP-resistant OS cells. (a) After transfection of sh-LINC00641, detection of LINC00641; (b) IC50 values; (c-d) MTT assay to detect the survival rate and cell proliferation; (e-f) Flow cytometry to analyze the apoptosis rate; (g-h) RT-qPCR to examine Bcl-2 and Bax mRNA; (i-j) Western blotto detect caspase-3 and PARP protein expression; (k-l) Transwell method to analyze the migration and invasion abilities. In Saos2/DDP and MG-63/DDP cells. * vs. the sh-NC, P < 0.05.

3.4 LINC00641 emerges direct binding and negative mediation to miR-320d

On the grounds of the detection of miR-320d, it was observed that miR-320d was silenced in OS tissues and cell lines in comparison with normal one, especially in DDP-resistant tumors and cell lines (). Then, miR-320d and LINC00641 were negatively absolutely correlated in OS patients (). The above results suggested that miR-320d was likely to be relevant to the resistance of LINC00641 toward DDP. In order to elaborate this hypothesis, bioinformatics was conducted to analyze and prefigure that miR-320d was the target of LINC00641 and they had hypothetic binding site (). Meanwhile, it was clear that transfection of miR-320d mimic reduced the luciferase activity of WT-LINC00641 reporter vector, but could not reduce that of empty vector or MUT-LINC00641 reporter vector in Saos2/DDP and MG-63/DDP cells (). Meanwhile, the further verification was gained that a direct reciprocity of miR-320d with LINC00641 existed, since noteworthy enrichment of LINC00641 was examined in both Saos2/DDP and MG-63/DDP cells (). RIP test found that LINC00641 and miR-320d were significantly enriched in the anti-ago2 group compared with the anti-igg group (). Furthermore, it was discovered that deletion of LINC00641 accelerated miR-320d in Saos2/DDP and MG-63/DDP cells (). Therefore, LINC00641 was able to negatively modulate miR-320d.

Figure 4. LINC00641 directly binds to miR-320d to negatively control it. (a-b) RT-qPCR to detect miR-320d in DDP-sensitive or drug-resistant cancer tissues and cell lines; (c) The analysis of the pertinence of LINC00641 with miR-320d; (d) The binding sequence of LINC00641 and miR-320d was predicted; (e-f) LINC00641-WT/MUT or miR-320d mimic/miRNC co-transfected into Saos2/DDP and MG-63/DDP cells for the luciferase activity detection; (g) RT-qPCR to examine LINC00641 after biotinylated miR-320d/NC; (H/I) RIP to detect the enrichment of LINC00641 or miR-320d in anti-ago2 combination precipitation. (j) RT-qPCR to detect miR-320d in Saos2/DDP and MG-63/DDP cells transfected with sh-NC/LINC00641. *P < 0.05.

Figure 4. LINC00641 directly binds to miR-320d to negatively control it. (a-b) RT-qPCR to detect miR-320d in DDP-sensitive or drug-resistant cancer tissues and cell lines; (c) The analysis of the pertinence of LINC00641 with miR-320d; (d) The binding sequence of LINC00641 and miR-320d was predicted; (e-f) LINC00641-WT/MUT or miR-320d mimic/miRNC co-transfected into Saos2/DDP and MG-63/DDP cells for the luciferase activity detection; (g) RT-qPCR to examine LINC00641 after biotinylated miR-320d/NC; (H/I) RIP to detect the enrichment of LINC00641 or miR-320d in anti-ago2 combination precipitation. (j) RT-qPCR to detect miR-320d in Saos2/DDP and MG-63/DDP cells transfected with sh-NC/LINC00641. *P < 0.05.

3.5 Knock-down of LINC00641 gene represses DDP-resistance of DDP-resistant OS cells via modulating miR-320d

Saos2/DDP and MG-63/DDP cells were transfected with miR-NC/320d, sh-NC/LINC00641 or sh-LINC00641 + anti-miR-NC/miR-320d to explore whether miR-320d participated in drug resistance of LINC00641 (). After transfection, the progression of DDP-resistant cells was observed. The results showed that overexpression of miR-320d inhibited the IC50 value (), proliferation (), migration () and invasion (), but promoted apoptosis of DDP drug-resistant cells (). In addition, low expression of miR-320d reversed the therapeutic effect of repressing LINC00641 on DDP resistant cells, significantly increasing IC50 value (), proliferation (), migration (), invasion () and clearly decreasing apoptosis (). As has been noted, it became clear that knock-down of LINC00641 was capable of restraining DDP-resistance of DDP-resistant OS cells via targeting miR-320d.

Figure 5. Knock-down LINC00641 is capable of reducing DDP-resistance of DDP-resistant OS cells via mediating miR-320d. (a) qPCR to detect miR-320d in cells; (b) IC50 value of cells in each group; (c-d) MTT assay to detect the survival rate and proliferation; (e) Flow cytometry to analyze apoptosis rates; (f-g) qPCR to examine Bcl-2 and Bax mRNA; (h-i) Western blot to detect caspase-3 and PARP protein expression; (j-k) Transwell method to analyze the migration and invasion abilities. In C-I, in Saos2/DDP and MG-63/DDP cells. * vs. the miR-NC, P < 0.05; + vs. the sh-LINC00641 + anti-NC, P < 0.05.

Figure 5. Knock-down LINC00641 is capable of reducing DDP-resistance of DDP-resistant OS cells via mediating miR-320d. (a) qPCR to detect miR-320d in cells; (b) IC50 value of cells in each group; (c-d) MTT assay to detect the survival rate and proliferation; (e) Flow cytometry to analyze apoptosis rates; (f-g) qPCR to examine Bcl-2 and Bax mRNA; (h-i) Western blot to detect caspase-3 and PARP protein expression; (j-k) Transwell method to analyze the migration and invasion abilities. In C-I, in Saos2/DDP and MG-63/DDP cells. * vs. the miR-NC, P < 0.05; + vs. the sh-LINC00641 + anti-NC, P < 0.05.

3.6 MCL1 is the alternative target of miR-320d

Detection of MCL1 was in DDP-sensitive and -resistant OS tissues and cells. The results clarified that MCL1 was saliently elevated in OS tissues and cells, especially in DDP-resistant tumors and cell lines (Saos2/DDP and MG-63/DDP) ( and b). Moreover, it was found that clinical data manifested miR-320d was negatively associated with MCL1, while LINC00641 was positively linked with MCL1 (). Therefore, it was presumed that a certain pertinence of MCL1 with miR-320d existed. Then starBase v3.0 was exerted to predict underlying targets of miR-320d. The results illustrated that MCL1 3ʹuntranslated region (UTR) and miR-320d had complementary binding sites (). In MG63/DDP and Saos-2/DDP cells transfected with miR-320d, the luciferase activity of MCL1 3ʹUTR-WT reporter gene emerged acute declination, while in any group that of MCL1 3ʹUTR-MUT reporter gene did not present conspicuous fluctuation (). RIP test found that miR-320d and MCL1 were significantly enriched in the anti-ago2 group compared with the anti-igg group ().Additionally, MCL1 were dramatically silenced in MG63/DDP and Saos-2/DDP cells transfected with sh-LINC00641 or miR-320d ( and k). The above results concluded that miR-320d is capable of interacting with MCL1, to attenuate MCL1 in MG63/DDP and Saos-2/DDP cells.

Figure 6. MCL1 is the alternative target of miR-320d. (a) RT-qPCR to detect MCL1 in normal tissues and DDP-sensitive or drug-resistant cancer tissues; (b) RT-qPCR to examine MCL1 in hFOB1.19 cells, OS cells and OS DDP-resistant cell lines; (c-d) Pearson test to analyze the pertinence of MCL1, LINC00641 with miR-320d; (e) The binding site of MCL1 with miR-320d predicted via bioinformatics website; (f-g) The luciferase activity assay to testify the targeting of MCL1 with miR-320d; (h-i) RIP to detect the enrichment of LINC00641 or miR-320d in anti-ago2 combination precipitation; J-K. qPCR and Western blot to detect MCL1. * vs. the sh-NC, P < 0.05; + vs. the miR-NC, P < 0.05.

Figure 6. MCL1 is the alternative target of miR-320d. (a) RT-qPCR to detect MCL1 in normal tissues and DDP-sensitive or drug-resistant cancer tissues; (b) RT-qPCR to examine MCL1 in hFOB1.19 cells, OS cells and OS DDP-resistant cell lines; (c-d) Pearson test to analyze the pertinence of MCL1, LINC00641 with miR-320d; (e) The binding site of MCL1 with miR-320d predicted via bioinformatics website; (f-g) The luciferase activity assay to testify the targeting of MCL1 with miR-320d; (h-i) RIP to detect the enrichment of LINC00641 or miR-320d in anti-ago2 combination precipitation; J-K. qPCR and Western blot to detect MCL1. * vs. the sh-NC, P < 0.05; + vs. the miR-NC, P < 0.05.

3.7 Refrained LINC00641 gene antagonizes DDP-resistance of DDP-resistant OS cells via modulating MCL1

For the sake of verifying whether MCL1 was implicated in DDP-resistance triggered via LINC00641, sh-NC/LINC00641 and sh-LINC00641 + MCL1/pcDNA were transfected into Saos2/DDP and MG-63/DDP cells. After validating transfection efficiency, the results uncovered that MCL1 was ascending in cells transfected with MCL1 (). Examination of the advancement of DDP-resistant cells was conducted. As a matter of reality, the results indicated that transfection of MCL1 was able to weaken repressive effect in the progression of Saos2/DDP and MG-63/DDP cells triggered via deletion of LINC00641, repressing IC50 value (), proliferation (), migration (), invasion (), but promoting apoptosis ().It was validated that knockdown of LINC00641 gene was inclined to emerge antagonism toward drug resistance of DDP-resistant OS cells via the mediation of MCL1. Besides, on the basis of the above results, knock-down of LINC00641 gene tended to block MCL1 via emulative interaction with miR-320d, thus expediting the sensitivity of OS toward DDP.

Figure 7. Restrained LINC00641 gene emerges antagonism toward DDP resistance of DDP-resistant OS cells via controlling MCL1. (a) RT-qPCR to detect MCL1; (b) IC50 value of cells in each group; (c-d) MTT assay to examine the proliferation ability; (e) Flow cytometry to detect apoptosis rates; (f-g) RT-qPCR to examine Bcl-2 and Bax mRNA; (h-i) Western blot to detect caspase-3 and PARP protein expression; (j-k) Transwell method to analyze the migration and invasion abilities. Apart from B, in Saos2/DDP and MG-63/DDP cells. * vs. the pcDNA, P < 0.05; # vs. the sh-LINC00641 + pcDNA, P < 0.05.

Figure 7. Restrained LINC00641 gene emerges antagonism toward DDP resistance of DDP-resistant OS cells via controlling MCL1. (a) RT-qPCR to detect MCL1; (b) IC50 value of cells in each group; (c-d) MTT assay to examine the proliferation ability; (e) Flow cytometry to detect apoptosis rates; (f-g) RT-qPCR to examine Bcl-2 and Bax mRNA; (h-i) Western blot to detect caspase-3 and PARP protein expression; (j-k) Transwell method to analyze the migration and invasion abilities. Apart from B, in Saos2/DDP and MG-63/DDP cells. * vs. the pcDNA, P < 0.05; # vs. the sh-LINC00641 + pcDNA, P < 0.05.

3.8 Depressive LINC00641 constrains OS tumor and enhances DDP sensitive growth

Next, OS mice xenograft tumor model was established in vivo to further figure out the impact of LINC00641 deletion on tumor growth. As shown in (), DDP therapy or declination of LINC00641 gave rise to a decrease in tumor volume and weight, which indicated that DDP therapy or silenced LINC00641 repressed the growth of OS tumor in vivo. In the meantime, repressive effect on tumor growth was manifested more observably via the combination of sh-LINC00641 with DDP therapy. Apart from that, the apoptosis of tumor tissues was detected and it displayed that DDP therapy or silenced LINC00641 boosted the apoptosis of cells. The effect of boosting cell apoptosis was exhibited more signally via the combination of sh-LINC00641 with DDP (), which suggested that deletion of LINC00641 was able to facilitate anti-tumor effect of DDP in vivo.

Figure 8. Depressive LINC00641 suppresses OS tumor and escalates development of DDP sensitivity. (a) Tumor volume of mice in each group; (b) tumor images and weight variations of mice in each group; (c) TUNEL staining to examine apoptotic cells. * vs. the sh-NC + PBS, P < 0.05; # vs. the sh-LINC00641 + PBS, P < 0.05.

Figure 8. Depressive LINC00641 suppresses OS tumor and escalates development of DDP sensitivity. (a) Tumor volume of mice in each group; (b) tumor images and weight variations of mice in each group; (c) TUNEL staining to examine apoptotic cells. * vs. the sh-NC + PBS, P < 0.05; # vs. the sh-LINC00641 + PBS, P < 0.05.

4 Discussion

DDP has been identified as the pivotal medicine for cancer chemotherapy since its discovery [Citation30], has been widely used in the treatment of various solid tumors. It is a metal-containing drug, which does not make it cross-resistant to other drugs in some key pathways [Citation31]. Therefore, only elucidating the mechanism of OS resistance to DDP has important therapeutic value. The purpose of this study was to investigate the action mechanism of LINC00641 in DDP resistance in OS. The results showed that LINC00641 knockdown reduced DDP resistance of Saos2/DDP and MG-63/DDP cells, repressed cell proliferation, migration and invasion, but promoted apoptosis through miR-320d/MCL1 axis.

In recent years, a great deal of studies hasve verified the pertinence of maladjustment of lncRNAs with DDP-resistance in tumors involving cancer, and investigated numerous lncRNAs linked with DDP-resistance. For instance, ROR is implicated in OS resistance toward DDP via miR-153-3p/ABCB1 axis and unlocks novel mediative axis of OS cells toward DDP, which is likely to provide therapeutic targets for OS patients [Citation32]; and the study of LINC00641 in OS and its influences on OS prognosis and DDP resistance have remained still obscure. In this study, it was found that LINC00641 was elevated in OS tissues and cells. At the same time, LINC00641 was dramatically augmented in DDP-resistant tumors and cell lines, insinuating that LINC00641 took on a latent regulatory effect in DDP-resistance of OS. On this basis, advancement of DDP-resistant OS cells was detected, and it was discovered that knock-down of LINC00641 gene was able to restrain advancement and metastasis of DDP-resistant cells thus weakening their DDP resistance. In the meantime, steadily transfected MG-63/DDP cells were utilized to construct xenografted tumor model in OS mice. The results manifested that DDP therapy or silenced LINC00641 contributed to a declination in tumor volume and weight, which signified that DDP therapy or depressive LINC00641 attenuated the development of OS tumor in vivo; Meanwhile, the combination of sh-LINC00641 and DDP had a more remarkable repressive effect on tumor growth, which was consistent with the experimental results in vitro. Taking the above biological behavior of LINC00641 into consideration, the pertinence of LINC00641 with overall survival period was evaluated, and it was determined that ascending LINC00641 was capable of pointing out the unpleasing prognostic of OS patients.

It is reported that lncRNA is utilized as competitive endogenous RNA (ceRNA) of miRNAs, and to modulate gene indirectly in miscellaneous biological processes involving cancer [Citation33]. Therefore, it was guessed that whether LINC00641 performed as a ceRNA to mediate miR-320d/MCL1 axis and taking on effect in stimulating cancer in OS. In line with bioinformatics prediction program, miR-320d manifested as the target of LINC00641, and an immediate binding of LINC00641 to miR-320d was verified. Then, the analysis of miR-320d was conducted and it was observed that miR-320d was silenced in OS tissues and cell lines. In OS patients, miR-320d and LINC00641 were immediately validated to be completely negatively correlated.

Studies have elucidated members of miR-320 family have been reported to be implicated in multiple malignant tumors, including glioma [Citation34], colon cancer [Citation35], large B-cell lymphoma [Citation36]. Nevertheless, the research on the effect and mechanism of miR-320d in the development of OS and its DDP-resistance has remained still scarce. For the sake of verifying whether miR-320d was relevant to DDP-resistance incurred via LINC00641, augmentation of miR-320d and rescue experiments were carried out. The results uncovered that miR-320d dampened the activity and advancement of DDP-resistant OS cells, while transfection of miR-320d inhibitor could restore the repressive effect of DDP-resistance of DDP-resistant OS cells triggered via deletion of LINC00641. Prior to this, LINC00641 ameliorated DDP-resistance of DDP-resistant OS cells via targeting miR-320d.

Many evidences have expounded that MCL1 also participates in tumor development and is under the mediation of miRNA. For instance, miR-153-3p modulates development of ovarian cancer in vivo and in vitro via targeting MCL1 gene [Citation37]. In this study, it was found that MCL1 was signally augmented in OS tissues and cells and was targeted via miR-320d. Moreover, it was observed via transfection of elevated MLC1 plasmid into cells with suppressive LINC00641 that enhancive MLC1 could reverse the impact of LINC00641 deletion on DDP-resistant OS cells, indicating that knock-down of LINC00641 gene was inclined to antagonize resistance of DDP-resistant OS cells via controlling MCL1.

Nevertheless, there are still some limitations in this study. In vivo experiments were worth conducting to further explore the function of miR-320d/MCL1, the downstream gene of LINC00641, in OS DDP-resistance. Besides, limited patient specimens might not completely validate accuracy of the results.

5 Conclusion

In a word, the study illustrates that LINC00641 is elevated in OS tissues and cells, and its augmentation is connected with unpleasing prognosis of OS patients. Moreover, declining LINC00641 could restrain the activity of OS cells, and further mechanism study uncovers that LINC00641 principally took on effect in modulating miR-320d/MCL1 axis. The discoveries suggest that LINC00641 is likely to be a novel molecular marker and therapeutic target of OS. To sum up, this study testifies that LINC00641 facilitates the advancement and metastasis of OS cells via negatively modulating miR-320d/MCL1, thus expediting the occurrence and development of OS.

Supplemental material

Supplemental Material

Download Zip (156.6 KB)

Disclosure statement

No potential conflict of interest was reported by the author(s).

Supplementary material

Supplemental data for this article can be accessed here.

Additional information

Funding

The author(s) reported there is no funding associated with the work featured in this article.

References

  • Zhu KP, Ma XL, Zhang CL. Overexpressed circPVT1, a potential new circular RNA biomarker, contributes to doxorubicin and cisplatin resistance of osteosarcoma cells by regulating ABCB1. Int J Biol Sci. 2018;14(3):321–330.
  • Liu K, Hou Y, Liu Y, et al. LncRNA SNHG15 contributes to proliferation, invasion and autophagy in osteosarcoma cells by sponging miR-141. J Biomed Sci. 2017;24(1):46.
  • Serra M, Hattinger CM. The pharmacogenomics of osteosarcoma. Pharmacogenomics J. 2017;17(1):11–20.
  • Wang H, Qin R, Guan A, et al. HOTAIR enhanced paclitaxel and doxorubicin resistance in gastric cancer cells partly through inhibiting miR-217 expression. J Cell Biochem. 2018;119(9):7226–7234.
  • Li Z, Liang Z, Wang Q. Overexpression of long non-coding RNA HOTTIP increases chemoresistance of osteosarcoma cell by activating the Wnt/b-catenin pathway. Am J Transl Res. 2016;8(5):2385–2393.
  • Liu C, Lin J. Long noncoding RNA -AZS1 acts as an oncogene in osteosarcoma by epigenetically activating KCNC4. Am J Transl Res. 2016;8(10):4095–4105.
  • Sun L, Sun P, Zhou Q, et al. Long noncoding RNA MALAT1 promotes uveal melanoma cell growth and invasion by silencing of miR-140. Am J Transl Res. 2016;8(10):3939–3946.
  • Hanwen Z, Shuya L, Lian T, et al. Long non-coding RNA (LncRNA) MRPL23-AS1 promotes tumor progression and carcinogenesis in osteosarcoma by activating Wnt/β-catenin signaling via inhibiting microRNA miR-30b and upregulating myosin heavy chain 9 (MYH9). Bioengineered. 2021;12:162–171.
  • Weiying Z, Qiang W, Du H, et al. CRISPR/cas9-mediated overexpression of long non-coding RNA SRY-box transcription factor 21 antisense divergent transcript 1 regulates the proliferation of osteosarcoma by increasing the expression of mechanistic target of rapamycin kinase and Kruppel like factor 4. Bioengineered. 2021. DOI:10.1080/21655979.2021.1995106
  • Lei S, Zhigang Z, Yibo G, et al. Long noncoding RNA OIP5-AS1 causes cisplatin resistance in osteosarcoma through inducing the LPAATβ/PI3K/AKT/mTOR signaling pathway by sponging the miR-340-5p. J Cell Biochem. 2019;120:9656–9666.
  • Xingxing S, Cong T, Hui Z, et al. Long noncoding RNA OIP5-AS1 mediates resistance to doxorubicin by regulating miR-137-3p/PTN axis in osteosarcoma. Biomed Pharmacother. 2020;128:110201.
  • Jianping Z, Shengming J, Wenjun X, et al. Long noncoding RNA LINC00641 promotes renal cell carcinoma progression via sponging microRNA-340-5p. Cancer Cell Int. 2021;21:210.
  • Zhijia L, Shizhe H, Zhihu L. LncRNA LINC00641 predicts prognosis and inhibits bladder cancer progression through miR-197-3p/KLF10/PTEN/PI3K/AKT cascade. Biochem Biophys Res Commun. 2018;503:1825–1829.
  • Yi L, Lily Z, Pei Z, et al. Long non-coding RNA LINC00641 suppresses non-small-cell lung cancer by sponging miR-424-5p to upregulate PLSCR4. Cancer Biomark. 2019;26:79–91.
  • Wang J, Liu Z, YuL-J H. Long non-coding RNA LINC00641 promotes cell growth and migration through modulating miR-378a/ZBTB20 axis in acute myeloid leukemia. Eur Rev Med Pharmacol Sci. 2019;23:7498–7509.
  • Bao L, Li X. MicroRNA-32 targeting PTEN enhances M2 macrophage polarization in the glioma microenvironment and further promotes the progression of glioma. Mol Cell Biochem. 2019;460(1–2):67–79.
  • Donzelli S, Mori F, Biagioni F, et al. MicroRNAs: short non-coding players in cancer chemoresistance. Mol Cell Ther. 2014;2:16–27.
  • Xin L, Yingjie L, Zhanye L, et al. Exosomal long non-coding RNA LINC00662 promotes non-small cell lung cancer progression by miR-320d/E2F1 axis. Aging. 2021;13:6010–6024.
  • Xiaojie C, Shegan G, Zhiwei Z, et al. MicroRNA-320d regulates tumor growth and invasion by promoting FoxM1 and predicts poor outcome in gastric cardiac adenocarcinomax. Aging. 2020;10:80.
  • Ivan D, Alexis D, You-Wen H. The anti-apoptotic Bcl-2 family member Mcl-1 promotes T lymphocyte survival at multiple stages. J Immunol. 2008;181:521–528.
  • Inoue S, Walewska R, Dyer MJS, et al. Downregulation of Mcl-1 potentiates HDACi-mediated apoptosis in leukemic cells. Leukemia. 2008;22:819–825.
  • Céline G, White E. BH3-only proteins in control: specificity regulates MCL-1 and BAK-mediated apoptosis. Genes Dev. 2005;19:1263–1268.
  • Ming L, Ma W. miR-26a reverses multidrug resistance in osteosarcoma by targeting MCL1. Front Cell Dev Biol. 2021;9:645381.
  • Jiankuo G, Dongmei D, Tianlun Z, et al. HOTAIR promotes cisplatin resistance of osteosarcoma cells by regulating cell proliferation, invasion, and apoptosis via miR-106a-5p/STAT3 axis. Cell Transplant. 2020;29:963689720948447.
  • Yali Y, Leiming W, Zijun L, et al. Long noncoding RNA CRNDE functions as a diagnostic and prognostic biomarker in osteosarcoma, as well as promotes its progression via inhibition of miR-335-3p. J Biochem Mol Toxicol. 2021;35:e22734.
  • Xinxin C, Chengyong Z, Xiao W. Long noncoding RNA DLEU1 aggravates osteosarcoma carcinogenesis via regulating the miR-671-5p/DDX5 axis. Artif Cells Nanomed Biotechnol. 2019;47:3322–3328.
  • Jiang L, Quan-Min W, Xiao-Qing W, et al. Long noncoding RNA miR210HG sponges miR-503 to facilitate osteosarcoma cell invasion and metastasis. DNA Cell Biol. 2017;36:1117–1125.
  • Yang G, Chengtao L, Xiaoling Z, et al. hsa_circ_0000006 induces tumorigenesis through miR-361-3p targeting immunoglobulin-like domains protein 1 (LRIG1) in osteosarcoma. Ann Transl Med. 2021;9:1242.
  • Wenpeng Z, Jiang G, Honglu L, et al. FAM83H-AS1/miR-485-5p/MEF2D axis facilitates proliferation, migration and invasion of hepatocellular carcinoma cells. BMC Cancer. 2021;21:1310.
  • Shen DW, Pouliot LM, Hall MD, et al. Cisplatin resistance: a cellular self-defense mechanism resulting from multiple epigenetic and genetic changes. Pharmacol Rev. 2012;64(3):706–721.
  • Shaloam D, Bernard Tchounwou P. Tchounwou Paul Bernard, Cisplatin in cancer therapy: molecular mechanisms of action. Eur J Pharmacol. 2014;740:364–378.
  • Cheng F-H, Zhao Z-S, Liu W-D. Long non-coding RNA ROR regulated ABCB1 to induce cisplatin resistance in osteosarcoma by sponging miR-153-3p. Eur Rev Med Pharmacol Sci. 2019;23(17):7256–7265.
  • Jin-Yan W, Yang Y, Yajun M, et al. Potential regulatory role of lncRNA-miRNA-mRNA axis in osteosarcoma. Biomed Pharmacother. 2020;121:109627.
  • Chong-Zhen Q, Qiao-Li L, Yan-Tao Y, et al. Downregulation of MicroRNA-320d predicts poor overall survival and promotes the growth and invasive abilities in glioma. Chem Biol Drug Des. 2017;89:806–814.
  • Toshihiro T, Yoichi K, Shin H, et al. MicroRNA-320 family is downregulated in colorectal adenoma and affects tumor proliferation by targeting CDK6. World J Gastrointest Oncol. 2016;8:532–542.
  • Peng Yan W, Xu Dong Z, Jiang Z, et al. Low expression of microRNA-146b-5p and microRNA-320d predicts poor outcome of large B-cell lymphoma treated with cyclophosphamide, doxorubicin, vincristine, and prednisone. Hum Pathol. 2014;45:1664–1673.
  • Cui L, Yinmin Z, Zhao W, et al. miR-153-3p regulates progression of ovarian carcinoma in vitro and in vivo by targeting MCL1 gene. J Cell Biochem. 2019;120:19147–19158.