- Open Access
LncRNA CBR3-AS1 regulates of breast cancer drug sensitivity as a competing endogenous RNA through the JNK1/MEK4-mediated MAPK signal pathway
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 41 (2021)
Adriamycin (ADR) resistance is one of the main obstacles to improving the clinical prognosis of breast cancer patients. Long noncoding RNAs (lncRNAs) can regulate cell behavior, but the role of these RNAs in the anti-ADR activity of breast cancer remains unclear. Here, we aim to investigate the imbalance of a particular long noncoding RNA, lncRNA CBR3 antisense RNA 1 (CBR3-AS1), and its role in ADR resistance.
Microarray analysis of ADR-resistant breast cancer cells was performed to identify CBR3-AS1. CCK-8 and colony formation assays were used to detect the sensitivity of breast cancer cells to ADR. Dual-luciferase reporter, RNA pulldown, IHC and western blot analyses were used to verify the relationship between the expression of CBR3-AS1, miRNA and target genes. For in vivo experiments, the effect of CBR3-AS1 on breast cancer resistance was observed in a xenograft tumor model. The role of CBR3-AS1 in influencing ADR sensitivity was verified by clinical breast cancer specimens from the TCGA, CCLE, and GDSC databases.
We found that CBR3-AS1 expression was significantly increased in breast cancer tissues and was closely correlated with poor prognosis. CBR3-AS1 overexpression promoted ADR resistance in breast cancer cells in vitro and in vivo. Mechanistically, we identified that CBR3-AS1 functioned as a competitive endogenous RNA by sponging miR-25-3p. MEK4 and JNK1 of the MAPK pathway were determined to be direct downstream proteins of the CBR3-AS1/miR-25-3p axis in breast cancer cells.
In summary, our findings demonstrate that CBR3-AS1 plays a critical role in the chemotherapy resistance of breast cancer by mediating the miR-25-3p and MEK4/JNK1 regulatory axes. The potential of CBR3-AS1 as a targetable oncogene and therapeutic biomarker of breast cancer was identified.
Worldwide, breast cancer is the most common malignant tumor in women [1,2,3]. Every year, an estimated 1 million new breast cancer cases are diagnosed worldwide. Adriamycin (ADR) is considered to be one of the most effective drugs for breast cancer, especially after tamoxifen treatment fails. However, its efficacy as a chemotherapeutic drug is greatly reduced due to intrinsic and acquired drug resistance during treatment [4, 5]. The underlying mechanisms resulting in ADR resistance in breast cancer remain poorly understood.
LncRNAs are a group of noncoding RNA molecules more than 200 nucleotides in length . According to reports, lncRNAs are involved in many biological processes related to carcinogenesis and drug resistance [7,8,9,10]. For example, TINCR promotes epithelial-mesenchymal transition and trastuzumab resistance by targeting miR-125b . HOTAIR mediates tamoxifen resistance by promoting estrogen receptor transcriptional activity . DSCAM-AS1 mediates tumor progression and tamoxifen resistance by regulating the estrogen receptor (ER) . However, the role of lncRNAs in the process of ADR resistance in breast cancer is still not fully explored.
LncRNA CBR3-AS1 is located in the antisense region of carbonyl reductase 3 (CBR3). It was initially discovered that CBR3-AS1 was abnormally expressed in prostate cancer cell lines and tissues [14,15,16]. Further studies showed that up-regulation of CBR3-AS1 expression promotes the proliferation of prostate cancer cells and inhibit apoptosis. Subsequent studies found that overexpression of CBR3-AS1 was observed in gastric carcinoma, esophageal squamous cell carcinoma, osteosarcoma, colorectal cancer and retinoblastoma [17,18,19,20,21]. However, the biological function and clinical significance of CBR3-AS1 in breast cancer drug resistance remain unclear.
In our study, we screened differentially expressed lncRNAs between ADR-resistant breast cancer cells and their parent cells. We found that the expression of CBR3-AS1 was up-regulated in ADR-resistant breast cancer cells, and its high expression was associated with poor prognosis of breast cancer. The overexpression of CBR3-AS1 promoted ADR resistance whereas knockdown of CBR3-AS1 expression significantly repressed it both in vivo and in vitro. Further studies showed that CBR3-AS1 functions as a ceRNA (competing endogenous RNA) by binding miR-25-3p. JNK1 and MEK4 are downstream proteins of the CBR3-AS1/miR-25-3p axis in breast cancer cells. Our results confirm for the first time that CBR3 promotes ADR resistance in breast cancer cells via regulation of the CBR3-AS1/miR-25-3p/MEK4/JNK1 axis, indicating that CBR3-AS1 has potential as a targetable oncogene and biomarker for breast cancer treatment.
Access and analysis of public data
GSE20685 was downloaded from the GEO dataset (http://www.ncbi.nlm.nih.gov/geo/). The whole genome expression profiles and clinicopathological information of human cancer were downloaded from The Cancer Genome Atlas (TCGA) (https://tcga-data.nci.nih.gov/). The expression of lncRNAs in the cell lines was obtained from the Cancer Cell Line Encyclopedia (CCLE) database (http://www.broadinstitute.org/ccle). The IC50 values of drug in multiple cell lines were obtained from the Genomics of Drug Sensitivity in Cancer (GDSC) database. All transcripts were normalized by log2 transformation. The Diana lncbase (http://carolina.imis.athena-innovation.gr/diana_tools) was used to predict the binding of miRNAs and lncRNAs. The bioinformatics analysis tool miRdb (http://www.miRdb.org/) was used to predict the binding of miRNA and mRNA 3′UTR. The RNAhybrid website (https://bibiserv.cebitec.uni-bielefeld.de/rnahybrid/) was used to predict the binding sites of miRNA and lncRNA. The bioinformatics website lncLocator (http://www.csbio.sjtu.edu.cn/cgi-bin/lncLocator.py) was used to predict the subcellular localization of lncRNAs.
The correlation between genes was evaluated through the Pearson correlation test. The differences between breast cancer tumor and normal samples were evaluated by unpaired t-test. The log-rank test was used to test the difference in survival rate among different groups of patients.
Breast cancer specimens (n = 96) were obtained from hospitalized patients from November 2008 to June 2009. All patients were diagnosed with breast cancer at the First Affiliated Hospital of China Medical University. Before surgery, the patient did not receive chemotherapy or radiotherapy. This study was approved by the Ethics Committee of China Medical University (Approval number: AF-SOP- 07-1.1-01).
MCF-7, T47D, MDA-MB-231 and HEK-293 T cell lines were purchased from the Institute of Biochemistry and Cell Biology at the Chinese Academy of Sciences (Shanghai, China), and MCF-7/ADR cells were purchased from the Zhen Shanghai and Shanghai Industrial Co., Ltd. (Shanghai, China). The cells were cultured in L-15 medium, DMEM or RPMI 1640 medium (HyClone, USA) and cultured in a humidified incubator at 37 °C. Except for cells cultured in L-15 medium, all cells were incubated in an environment containing 5% carbon dioxide.
Microarray analysis was performed by Shanghai OE Biotech (China) with the microarray Agilent Human lncRNA 4*180 K (Design ID: 076500). The raw microarray data were submitted to the GEO database (accession number: GSE155478). The lncRNAs were identified by the Gencode database (https://www.gencodegenes.org). The p-value and fold change value of the t-test were used to screen differential lncRNAs. The screening criteria were up-regulated or down-regulated log2FC ≥ 2.0 and FDR ≤ 0.05.
Gene set enrichment analysis
Data on CBR3-AS1 expression were obtained from breast cancer patients with TCGA. Taking the median CBR3-AS1 expression as the cut-off point, breast cancer cases were divided into two groups. The index of gene sequencing was set as “Pearson”, the curve of the top set was set to 150, and all other parameters were set at the default values.
The cells were seeded on a 6-well plate and incubated for 24 h. The cells reached 70–80% confluence. Lipofectamine 3000 (Thermo Fisher Scientific, MA, USA) transfection reagent (5 μl) was dissolved in 125 μl of serum-free medium. Two micrograms of DNA was diluted with 125 μl of serum-free medium to prepare the DNA (Hanheng, China) premix, and then 10 μl of p3000 reagent was added and mixed well. After incubation for 5 min, the DNA solution was added to the cells. After 4 h, the medium was replaced with fresh medium containing 10% serum. The final transfection concentration of miRNA mimics (RiboBio, Guangzhou, China) was 50 μM, and the final concentration of siRNA (RiboBio, China) transfection was 50 μM. In order not to decompose RNA, p3000 was not used in the transfection process, and the other steps were equivalent to DNA transfection.
In the animal experiments, the plasmid was stably transfected into cells as described above. Then, 1 mg/L puromycin was added to the culture for 1 month, and stably transfected cell lines were selected for the experiment.
Cells were lysed on ice in RIPA buffer (Beyotime, China). The lysate was centrifuged at 15000×g at 4 °C for 15 min. SDS polyacrylamide gel electrophoresis was used to isolate proteins before they were transferred onto PVDF membranes. The membrane was blocked in 5% skim milk and incubated with the primary antibody in TBST at 4 °C overnight. Then, the cells were incubated with a horseradish peroxidase-labelled secondary antibody at a concentration of 1:10000 for 1 h. Finally, chemiluminescence was carried out using a Microchemi 4.2 system (DNR Bioimaging System, Israel). The details of all antibodies used are shown in Table S1.
HEK-293 T cells were transfected with MEK4–3′UTR-WT, MEK4–3′UTR-MU, JNK1–3′UTR-WT and JNK1–3′UTR-MU plasmids. At the same time, the four groups of cells were transfected with either NC mimic or miR-25-3p mimic to form eight experimental groups. Finally, each group of cells was transfected with the Renilla luciferase plasmid. After 36 h, the cells were harvested for luciferase analysis using a dual luciferase reporting and detection system (Promega, USA). The results were normalized against Renilla luciferase activity.
TRIzol (Invitrogen, USA) was used to extract total RNA from cultured cells. For mRNA and lncRNAs, a ReverTra Ace qPCR RT Kit (Toyobo, Japan) was used, and 200 ng of total RNA was used to synthesize cDNA in a 10 μl reaction volume. For miRNAs, the Bulge-Loop™ miRNA qRT-PCR Starter Kit (RiboBio, China) was used, and 1 μg of total RNA was used to synthesize cDNA in a 10 μl reaction volume. qRT-PCR of mRNA and lncRNAs was performed using SYBR Green Real-time PCR Master Mix (Toyobo, Japan) in a 12.5 μl reaction volume. The qRT-PCR of miRNAs used Bulge-LoopTM miRNA qRT-PCR Starter Kit (Ribobio, China) in a 20 μl reaction volume. The primers used are listed in Table S2.
Colony formation assay
MCF-7/ADR and MCF-7 cells (1 × 103/well) in the logarithmic growth phase were resuspended and seeded into a 60 mm cell culture dish. For the drug treatments, 100 nM ADR was added to MCF-7 cells and 30 μM ADR was added to MCF-7/ADR cells for 48 h. Afterward, the medium in the culture was replaced with normal medium. Two weeks later, the cells were washed with PBS twice, fixed with methanol, and stained with 0.1% crystal violet (Beyotime, China) for 30 min. The concentration of JNK1 inhibitor (SP600125, MedChemExpress, USA) used in the cell experiment was 40 nM.
MCF-7/ADR and MCF-7 cells that were transiently transfected were seeded into 96-well plates at 3 × 103 cells per well. After 24 h, the cells were treated with ADR (Sigma Chemical Co, St. Louis, MO) for 48 h. Next, cells were incubated with 10 μl of CCK-8 reagent (Dojindo, Japan) for 1 h before the absorbance was measured at 450 nm on a spectrophotometer.
ADR accumulation assay
According to the method of our previous articles , a cumulative ADR measurement assay was performed. Briefly, the cells were exposed to 5 μM ADR for 2 h and then washed with PBS. Next, flow cytometry was used to measure ADR fluorescence to determine the intracellular ADR level.
Flow cytometric analysis of apoptosis
MCF-7/ADR and MCF-7 cell apoptosis after ADR treatment was evaluated using the Annexin-V APC/7AAD Iodide Detection Kit (BD Biosciences, USA) according to the manufacturer’s instructions. Cells were then analysed by MACSQuant (Miltenyi Biotec, Germany).
The negative control plasmid (GFP alone) or OE-CBR3-AS1 plasmid were transfected into MCF-7 cells, which were then treated with puromycin to obtain stably expressing cell lines. All BALB/c mice (4 weeks old) were purchased from Shanghai Laboratory Animal Centre (Shanghai, China). MCF-7/control and MCF-7/CBR3-AS1 cells (1 × 107 cells/100 μl/nude mice) were subcutaneously injected into the axillary region of nude mice. When the subcutaneous tumors grew to 200–250 mm3 (16 days later), nude mice in the treatment group were intraperitoneally injected with ADR (2 mg/kg) and/or JNK1 inhibitor (15 mg/kg; MedChemExpress, USA) once every 3 days. All nude mice were treated 6 times. The diameter of the transplanted tumor was measured with digital calipers every week. The tumor volume was calculated as follows: volume = 0.5 × width2 × length. All nude mice were humanly sacrificed at the end of the experiment, and the tumor tissues were collected for further study by qRT-PCR, IHC and western blotting.
ISH (in situ hybridization) and IHC (immunohistochemistry)
The assay was carried out according to the method mentioned in a previous study . We collected breast cancer tissues from patients from the First Affiliated Hospital of China Medical University. The ISH probes were ordered from BOSTER Biological Technology Co., Ltd. (USA). For the IHC assays, we collected breast cancer tissues from different groups of subcutaneous. The details of all antibodies used are shown in Table S1.
MCF-7 cell lysates were incubated with CBR3-AS1 or control probe (Sangon, China) conjugated to streptavidin agarose resin beads (Thermo Fisher Scientific, USA) to generate probe-bound Dynabeads. After the samples were washed with buffer, DNase I (20 U) was added, and purified RNA was collected. The purified RNA was analysed by qRT-PCR. The whole experimental process and the buffer formulation were described by Hsu et al. .
Quantitative data are expressed as the means ± SD of at least three independent experiments. All experimental values were evaluated using GraphPad Prism 8.3.0 (GraphPad software). The unpaired t-test was used for statistical analysis between the two experimental groups. The relationship between CBR3-AS1 expression and clinicopathological factors was tested by the chi square test. In all cases, P < 0.05 was considered statistically significant.
CBR3-AS1 is up-regulated in drug-resistant breast cancer cells and is related to drug sensitivity and poor prognosis of breast cancer
To identify whether the differential expression of lncRNAs can be induced during the course of breast cancer drug resistance, we investigated the expression of lncRNAs in MCF-7 cells and drug-resistant MCF-7/ADR breast cancer cells. LncRNAs with log2FC > 1 and an FDR (false discovery rate) < 0.05 were selected from the microarray data. The results showed that, compared with MCF-7 cells, MCF-7/ADR cells had 949 up-regulated lncRNAs and 915 downregulated lncRNAs (Fig. 1a-b). The top 20 lncRNAs with the most obvious up-regulation and downregulation in MCF-7/ADR cells are listed in Fig. 1c. Moreover, in the three groups of MCF-7/ADR cells and MCF-7 cells, differential expression of lncRNAs was used for principal component analysis (PCA). The results showed that the differential expression of lncRNAs could effectively distinguish MCF-7/ADR and MCF-7 cells (Fig. S1 A). We further selected five lncRNAs (DSCR8, TP53TG1, CBR3-AS1, LINC01006, HAGLROS) that were not previously reported to be related to the drug sensitivity of breast cancer and performed qRT-PCR to verify their expression levels (Fig. 1d). Among them, CBR3-AS1 showed the highest expression in MCF-7/ADR cells. To confirm the relationship between the expression of the five lncRNAs and the sensitivity of breast cancer cells to ADR, we observed the expression patterns of lncRNAs in various breast cancer cells from the CCLE database and then obtained the IC50 values of ADR in these breast cancer cells from the GDSC database. The results showed that the expression of CBR3-AS1 was positively correlated with the resistance of breast cancer cells to ADR (p < 0.05, r > 0.3) (Fig. 1e, Fig. S1B). To verify this result, we performed CCK-8 experiments on four breast cancer cell lines. The IC50 value of ADR in each cell line was calculated (Fig. 1f). The qRT-PCR results showed that CBR3-AS1 was highly expressed in breast cancer cell lines with high IC50 values but showed reduced expression in cell lines with a low IC50 value (Fig. 1g). In addition, we obtained the mRNA expression profile and clinical information of breast cancer patients from the GSE20685 dataset in the GEO dataset and used the Kaplan-Meier plotter website to calculate whether the survival of patients with breast cancer with different expression patterns of the five lncRNAs was different. The results showed that in breast cancer patients receiving chemotherapy, the prognosis of patients with high expression of CBR3-AS1 was worse (Fig. S1C). We also found that in the TCGA database, CBR3-AS1 was up-regulated in breast cancer tissues compared with adjacent normal tissues (Fig. 1h). Gene set enrichment analysis (GSEA) showed that CBR3-AS1 was enriched in the ABC transporter signal transduction pathway (Fig. 1i). In addition, the expression of CBR3-AS1 was significantly correlated with the mRNA expression of p-glycoprotein (P-gp, ABCB1) and breast cancer resistance protein (BCRP, ABCG2) (Fig. 1j, k). The log-rank test of the OS curve of breast cancer patients in the TCGA also indicated that compared with low CBR3-AS1 expression, high CBR3-AS1 expression was significantly correlated with worse OS (Fig. 1l). We evaluated the expression of CBR3-AS1 via ISH in 96 breast cancer samples from the First Affiliated Hospital of China Medical University. High and low expression is shown in Fig. 1m. There was a significant correlation between CBR3-AS1 level and both tumor size and pathological staging (Table 1). The log-rank test of the overall survival curve of these breast cancer patients showed that elevated expression of CBR3-AS1 was significantly related to the worse OS of breast cancer (Fig. 1n).
Effect of CBR3-AS1 on the drug sensitivity of breast cancer cells
To identify the effect of CBR3-AS1 on the drug sensitivity of breast cancer cells, we overexpressed CBR3-AS1 in MCF-7 cells with low expression of CBR3-AS1 and silenced CBR3-AS1 in MCF-7/ADR cells with high expression of CBR3-AS1. The overexpression or silencing efficiency of CBR3-AS1 was evaluated by qRT-PCR. After transfection of the CBR3-AS1 overexpression plasmid, CBR3-AS1 levels were significantly increased compared with those in the control group (Fig. S2A). Upon completion of si-CBR3-AS1 siRNA transfection, CBR3-AS1 expression was significantly downregulated compared with that in the si-control group (Fig. S2B). The results showed that in cells overexpressing CBR3-AS1, the protein expression levels of P-gp increased significantly (Fig. 2a). After CBR3-AS1 was silenced, the protein expression levels of P-gp decreased significantly (Fig. 2b). Then, we used flow cytometry to detect drug accumulation in cells. The overexpression of CBR3-AS1 significantly reduced the accumulation of ADR in cells (Fig. 2c), while silencing of CBR3-AS1 significantly increased the accumulation of ADR in cells (Fig. 2d). Next, we used the CCK-8 cytotoxicity test to detect changes in the drug sensitivity of cells. After overexpressing CBR3-AS1, cells exhibited a decrease in drug sensitivity (Fig. 2e, Table 2). However, after CBR3-AS1 was silenced, the cells’ drug sensitivity was significantly increased (Fig. 2f). We performed CCK-8 cytotoxicity experiments with other anticancer drugs in cells overexpressing CBR3-AS1. The results showed that after overexpressing CBR3-AS1, MCF-7 cells showed decreased sensitivity to paclitaxel decreased, while their sensitivity to cisplatin did not change (Fig. S2C, D). The above results showed that the effects of CBR3-AS1 on drug sensitivity are not specific to ADR. Colony formation experiments showed that colony-forming ability of cells was significantly enhanced by overexpression of CBR3-AS1 after adding ADR (Fig. 2g). However, compared with the control cells, si-CBR3-AS1 cells showed reduced colony-forming ability (Fig. 2h). At the same time, the IC50(48 h) concentration of ADR was used to induce apoptosis, and flow cytometry showed that CBR3-AS1 overexpression significantly reduced apoptosis (Fig. 2i).
CBR3-AS1 regulates MEK4 / JNK1 expression through miR-25-3p
The bioinformatics analysis tool Diana lncbase revealed that lncRNA CBR3-AS1 targeted miRNAs (Fig. 3a, Fig. S3A). Research has shown that activation of miR-25-3p can overcome chemotherapy resistance . qRT-PCR analysis showed that miR-25-3p expression was downregulated in MCF-7/ADR cells (Fig. 3b). The bioinformatics analysis tool miRdb was used to discover the 3′UTR binding site for miR-25-3p with JNK1/MEK4 in the MAPK pathway (Table S3). CBR3-AS1 showed a significant negative correlation with miR-25-3p and positive correlation with JNK1/MEK4 mRNA expression in the TCGA breast cancer data (Fig. 3c). To determine whether JNK1 and MEK4 are possible targets for miR-25-3p, a luciferase reporter gene assay was conducted. We constructed luciferase reporter vectors Luc-JNK1 and Luc-MEK4 containing the 3′UTRs of JNK1 and MEK4 (Fig. 3d). Overexpression of miR-25-3p significantly reduced the luc-jnk1 luciferase activity in MCF-7 cells (Fig. 3e), indicating that JNK1 and MEK4 are direct targets of miR-25-3p. The RNAhybrid website predicted that CBR3-AS1 has a binding site with miR-25-3p (Fig. 3f). RNA pulldown assays showed that CBR3-AS1 could directly bind to miR-25-3p in breast cancer cells (Fig. 3g). Since the function of lncRNAs depends on their subcellular distribution, we examined the subcellular localization of CBR3-AS1. The bioinformatics website lncLocator predicted that CBR3-AS1 was expressed in both the nucleus and cytoplasm (Fig. S3B). The results were verified by qRT-PCR (Fig. 3h). qRT-PCR analysis and western blot analysis showed that the overexpression of CBR3-AS1 increased the mRNA and protein expression of JNK1 and MEK4 in MCF-7 cells after ADR was added (Fig. 3i, j), while silencing CBR3-AS1 reduced the mRNA and protein expression of JNK1 and MEK4 in MCF-7/ADR cells (Fig. 3k, l). Moreover, overexpression of miR-25-3p decreased JNK1 and MEK4 expression in MCF-7/ADR cells (Fig. S3C, D), while silencing miR-25-3p expression may increase JNK1 and mek4 expression in MCF-7/ADR cells (Fig. S3E, F). We next used ISH and IHC to examine 96 breast cancer tissues, sliced into serial sections (Fig. 3m). Consistent with the results in breast cancer cells, a negative correlation was confirmed to exist between CBR3-AS1 and miR-25-3p within this cohort (Fig. 3n), the expression of CBR3-AS1 was positively correlated with JNK1/MEK4 (Fig. S3G) and miR-25-3p was negatively correlated with JNK1/MEK4 (Fig. S3H). In conclusion, CBR3-AS1 regulates the expression of MEK4/JNK1 through miR-25-3p.
Silencing JNK1 or overexpressing miR-25-3p restored ADR drug sensitivity in breast cancer cells after overexpressing CBR3-AS1
To determine whether miR-25-3p and JNK1 are involved in CBR3-AS1-mediated chemotherapy resistance in breast cancer cells, a CCK-8 assay was conducted with MCF-7 cells, and the results showed that overexpression of miR-25-3p or inhibition of JNK1 reversed the decrease in ADR sensitivity caused by overexpression of CBR3-AS1 (Fig. 4a). Colony-forming experiments confirmed that overexpression of miR-25-3p or inhibition of JNK1 reversed the increased cell proliferation by overexpression of CBR3-AS1 in cells treated with ADR (Fig. 4b). Apoptosis experiments confirmed that overexpression of CBR3-AS1 led to a decrease in the number of apoptotic cells, which was reversed by overexpressing miR-25-3p or inhibiting JNK1 (Fig. 4c). qRT-PCR analysis and western blot analysis showed that overexpression of miR-25-3p reversed the activation of JNK1 and MEK4 caused by overexpression of CBR3-AS1 (Fig. 4d-e).
In vivo experiments showed that overexpression of CBR3-AS1 promoted the drug resistance of breast cancer cells to ADR
To determine whether the overexpression of CBR3-AS1 can reduce the sensitivity of tumors to ADR in vivo and whether this effect is achieved by regulating JNK1, MCF-7 cells stably transfected with CBR3-AS1-cDNA or NC-cDNA were subcutaneously injected into nude mice. A schematic of the experiment is shown in Fig. 5a. After 34 days, the expression level of CBR3-AS1 in subcutaneous tumor tissue was detected. The expression level of CBR3-AS1 in the transfected group was still significantly higher than that in the control group, confirming the construction of the model (Fig. 5b). Next, the expression level of miR-25-3p in the tumor was detected, and it was found that the expression level of miR-25-3p in the OE-CBR3-AS1 group was significantly reduced (Fig. 5c). In vivo imaging in mice showed that overexpression of CBR3-AS1 resisted the inhibitory effect of ADR on tumor cells, which was reversed by inhibition of JNK1 (Fig. 5d). Furthermore, the effect of CBR3-AS1 on increasing tumor size and weight could be reversed by a JNK1 inhibitor (Fig. 5e-g). Western blot analysis showed that overexpression of CBR3-AS1 increased the expression of JNK1 and MEK4 (Fig. 5h). JNK1 and MEK4 expression was also detected by IHC (immunohistochemistry) (Fig. 5i). These results showed that CBR3-AS1 promotes ADR resistance in breast cancer through JNK1/MEK4 in vivo.
In the past few decades, breast cancer has been treated with standard chemotherapy after surgical resection, which greatly improves the overall survival of patients. However, cancer-related deaths are still on the rise due to breast cancer metastasis and drug resistance. The MAPK pathway plays a key role in the invasion, migration and metastasis of various malignant tumors [26,27,28,29]. Shen et al. found that overexpression of lncRNA PCAT-1 can enhance cell division and inhibit apoptosis in myeloma cells through the p38 and JNK MAPK pathways . Xu et al. found that ROR2 facilitates epithelial-mesenchymal transition of breast cancer by regulating the MAPK/p38 signaling pathway . Another study reported that circ-MAPK4 plays a key role in the survival and apoptosis of glioma cells by regulating the miR-125a-3p and p38/MAPK signaling pathways . Another study reported that circ-MAPK4 plays a key role in the survival and apoptosis of glioma cells by regulating the miR-125a-3p and p38/MAPK signaling pathways .
JNK1 is one of the important genes in the MAPK pathway. MAPK pathway genes act as integration points of various biochemical signals and participate in many biological processes, such as differentiation, proliferation, apoptosis and transcriptional regulation [34, 35]. .MEK4 can phosphorylate the tyrosine residues of cytoplasmic JNK, resulting in its activation. Activated JNK quickly and significantly accumulates in the nucleus and causes the phosphorylation of the downstream transcription factors c-Jun, ATF-22, and Elk21; this in turn promotes the expression of related genes .
Some studies have shown that CBR3-AS1 is involved in and regulates the occurrence and development of cancers. It has been reported that CBR3-AS1 promotes osteosarcoma tumorigenesis and predicts poor prognosis . CBR3-AS1 promotes the development of colorectal cancer by regulating the PI3K/Akt signaling pathway . CBR3-AS1 plays an important role in the proliferation of esophageal squamous cells . CBR3-AS1 is highly expressed in retinoblastoma and mediates the proliferation of retinoblastoma cells by regulating CBR3 . However, it is not clear how CBR3-AS1 affects the drug resistance of breast cancer cells.
The ceRNA (competitive endogenous RNA) hypothesis suggests that a large-scale transcriptional regulatory network may be formed among RNAs, which can affect biological processes . Our study is based on the ceRNA hypothesis, that is, cbr3-as1 directly binds to miR-25-3p, thus targeting JNK1/MEK4, activating the MAPK pathway, and ultimately affecting the development of chemoresistance in breast cancer (Fig. 6).
In this study, we determined that CBR3-AS1 plays a novel role in breast cancer cell resistance. Overexpression of CBR3-AS1 confers chemotherapy resistance to breast cancer cells in vitro and in vivo. It is worth noting that CBR3-AS1 expression is up-regulated in breast cancer, while patients with high levels of CBR3-AS1 show poor overall survival. These findings indicate that CBR3-AS1 can be used as a predictor of drug response in breast cancer patients and may be a useful biomarker for precision medicine in the clinic.
Our study has some limitations. We did not detect the effect of CBR-AS1 on MCF-7/ADR apoptosis. This is because to induce apoptosis in these cells apoptotic, a very large dose of ADR is needed. Because ADR itself has strong fluorescence, large amounts of ADR will interfere with commonly used fluorescent dyes such as FITC, 7-AAD, PE and APC, and it is difficult to interpret the results. Similarly, we did not use MCF-7/ADR in animal experiments. Due to the cardiotoxicity of ADR, it is difficult for mice to bear such a large dose. MCF-7/ADR cells also have difficulty forming tumors in mice.
CBR3-AS1 might be an important biomarker for evaluating the prognosis of breast cancer patients. CBR3-AS1 increased the resistance of breast cancer to ADR by competitively binding miR-25-3p with JNK1/MEK4 and enhancing the MAPK signaling pathway. The CBR3-AS1/miR-25-3p/JNK1/MEK4 axis provides insight into breast cancer drug resistance mechanisms and theoretical support in the search for new breast cancer diagnostic markers and specific therapeutic targets.
Availability of data and materials
The datasets generated and analysed during the current study are available in the Gene Expression Omnibus data sets (https://www.ncbi.nlm.nih.gov/gds) (GSE155487).
Cancer Genome Atlas
Cell Counting Kit-8
Gene Expression Omnibus
Mitogen-activated protein kinase
c-Jun N-terminal kinase
Mitogen-activated extracellular signal-regulated kinase
In situ hybridization
Half maximal inhibitory concentration
Carbonyl reductase 3 antisense RNA 1
Global Burden of Disease Cancer C, Fitzmaurice C, Abate D, Abbasi N, Abbastabar H, Abd-Allah F, et al. Global, regional, and National Cancer Incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the global burden of disease study. JAMA Oncol. 2019;5:1749.
Liang Y, Zhang H, Song X, Yang Q. Metastatic heterogeneity of breast cancer: molecular mechanism and potential therapeutic targets. Semin Cancer Biol. 2020;60:14–27.
Yeo SK, Guan JL. Breast cancer: multiple subtypes within a tumor? Trends Cancer. 2017;3(11):753–60.
Batist G, Harris L, Azarnia N, Lee LW, Daza-Ramirez P. Improved anti-tumor response rate with decreased cardiotoxicity of non-pegylated liposomal doxorubicin compared with conventional doxorubicin in first-line treatment of metastatic breast cancer in patients who had received prior adjuvant doxorubicin: results of a retrospective analysis. Anti-Cancer Drugs. 2006;17(5):587–95.
Howard GR, Johnson KE, Rodriguez Ayala A, Yankeelov TE, Brock A. A multi-state model of chemoresistance to characterize phenotypic dynamics in breast cancer. Sci Rep. 2018;8(1):12058.
Lee H, Zhang Z, Krause HM. Long noncoding RNAs and repetitive elements: junk or intimate evolutionary partners? Trends Genet. 2019;35(12):892–902.
Lu C, Wei Y, Wang X, Zhang Z, Yin J, Li W, et al. DNA-methylation-mediated activating of lncRNA SNHG12 promotes temozolomide resistance in glioblastoma. Mol Cancer. 2020;19(1):28.
Peng WX, Koirala P, Mo YY. LncRNA-mediated regulation of cell signaling in cancer. Oncogene. 2017;36(41):5661–7.
Bhan A, Soleimani M, Mandal SS. Long noncoding RNA and cancer: a new paradigm. Cancer Res. 2017;77(15):3965–81.
Zhang S, Zheng F, Zhang L, Huang Z, Huang X, Pan Z, et al. LncRNA HOTAIR-mediated MTHFR methylation inhibits 5-fluorouracil sensitivity in esophageal cancer cells. J Exp Clin Cancer Res. 2020;39(1):131.
Dong H, Hu J, Zou K, Ye M, Chen Y, Wu C, et al. Activation of LncRNA TINCR by H3K27 acetylation promotes Trastuzumab resistance and epithelial-mesenchymal transition by targeting MicroRNA-125b in breast cancer. Mol Cancer. 2019;18(1):3.
Xue X, Yang YA, Zhang A, Fong KW, Kim J, Song B, et al. LncRNA HOTAIR enhances ER signaling and confers tamoxifen resistance in breast cancer. Oncogene. 2016;35(21):2746–55.
Niknafs YS, Han S, Ma T, Speers C, Zhang C, Wilder-Romans K, et al. The lncRNA landscape of breast cancer reveals a role for DSCAM-AS1 in breast cancer progression. Nat Commun. 2016;7:12791.
Fang Z, Xu C, Li Y, Cai X, Ren S, Liu H, et al. A feed-forward regulatory loop between androgen receptor and PlncRNA-1 promotes prostate cancer progression. Cancer Lett. 2016;374(1):62–74.
Zhang A, Zhang J, Kaipainen A, Lucas JM, Yang H. Long non-coding RNA: a newly deciphered “code” in prostate cancer. Cancer Lett. 2016;375(2):323–30.
Cui Z, Ren S, Lu J, Wang F, Xu W, Sun Y, et al. The prostate cancer-up-regulated long noncoding RNA PlncRNA-1 modulates apoptosis and proliferation through reciprocal regulation of androgen receptor. Urol Oncol. 2013;31(7):1117–23.
Zhang Y, Meng W, Cui H. LncRNA CBR3-AS1 predicts unfavorable prognosis and promotes tumorigenesis in osteosarcoma. Biomed Pharmacother. 2018;102:169–74.
Song W, Mei JZ, Zhang M. Long noncoding RNA PlncRNA-1 promotes colorectal cancer cell progression by regulating the PI3K/Akt signaling pathway. Oncol Res. 2018;26(2):261–8.
Wang CM, Wu QQ, Li SQ, Chen FJ, Tuo L, Xie HW, et al. Upregulation of the long non-coding RNA PlncRNA-1 promotes esophageal squamous carcinoma cell proliferation and correlates with advanced clinical stage. Dig Dis Sci. 2014;59(3):591–7.
Baratieh Z, Khalaj Z, Honardoost MA, Emadi-Baygi M, Khanahmad H, Salehi M, et al. Aberrant expression of PlncRNA-1 and TUG1: potential biomarkers for gastric cancer diagnosis and clinically monitoring cancer progression. Biomark Med. 2017;11(12):1077–90.
Wang S, Liu J, Yang Y, Hao F, Zhang L. PlncRNA-1 is overexpressed in retinoblastoma and regulates retinoblastoma cell proliferation and motility through modulating CBR3. IUBMB Life. 2018;70(10):969–75.
Zhao L, Wang Y, Jiang L, He M, Bai X, Yu L, et al. MiR-302a/b/c/d cooperatively sensitizes breast cancer cells to adriamycin via suppressing P-glycoprotein(P-gp) by targeting MAP/ERK kinase kinase 1 (MEKK1). J Exp Clin Cancer Res. 2016;35:25.
Song X, Zhang X, Wang X, Chen L, Jiang L, Zheng A, et al. LncRNA SPRY4-IT1 regulates breast cancer cell stemness through competitively binding miR-6882-3p with TCF7L2. J Cell Mol Med. 2020;24(1):772–84.
Hsu RJ, Yang HJ, Tsai HJ. Labeled microRNA pull-down assay system: an experimental approach for high-throughput identification of microRNA-target mRNAs. Nucleic Acids Res. 2009;37(10):e77.
Niu M, Feng Y, Zhang N, Shao T, Zhang H, Wang R, et al. High expression of miR-25 predicts favorable chemotherapy outcome in patients with acute myeloid leukemia. Cancer Cell Int. 2019;19:122.
Al-Khayal K, Vaali-Mohammed MA, Elwatidy M, Bin Traiki T, Al-Obeed O, Azam M, et al. A novel coordination complex of platinum (PT) induces cell death in colorectal cancer by altering redox balance and modulating MAPK pathway. BMC Cancer. 2020;20(1):685.
Alimirah F, Pulido T, Valdovinos A, Alptekin S, Chang E, Jones E, et al. Cellular senescence promotes skin carcinogenesis through p38MAPK and p44/p42MAPK signaling. Cancer Res. 2020;80:3606.
Franklin DA, Sharick JT, Gonzalez-Ericsson PI, Sanchez V, Dean P, Opalenik SR, et al. MEK activation drives glycolysis and supports suppressive myeloid cells in TNBC. JCI Insight. 2020;5(15):e134290.
Chakraborty R, Abdel-Wahab O, Durham BH. MAP-kinase-driven hematopoietic neoplasms: a decade of Progress in the molecular age. Cold Spring Harb Perspect Med. 2020;1:a034892.
Shen X, Shen P, Yang Q, Yin Q, Wang F, Cong H, et al. Knockdown of long non-coding RNA PCAT-1 inhibits myeloma cell growth and drug resistance via p38 and JNK MAPK pathways. J Cancer. 2019;10(26):6502–10.
Xu J, Shi J, Tang W, Jiang P, Guo M, Zhang B, et al. ROR2 promotes the epithelial-mesenchymal transition by regulating MAPK/p38 signaling pathway in breast cancer. J Cell Biochem. 2020;121(10):4142.
He J, Huang Z, He M, Liao J, Zhang Q, Wang S, et al. Circular RNA MAPK4 (circ-MAPK4) inhibits cell apoptosis via MAPK signaling pathway by sponging miR-125a-3p in gliomas. Mol Cancer. 2020;19(1):17.
Sheng W, Shi X, Lin Y, Tang J, Jia C, Cao R, et al. Musashi2 promotes EGF-induced EMT in pancreatic cancer via ZEB1-ERK/MAPK signaling. J Exp Clin Cancer Res. 2020;39(1):16.
Barr RK, Bogoyevitch MA. The c-Jun N-terminal protein kinase family of mitogen-activated protein kinases (JNK MAPKs). Int J Biochem Cell Biol. 2001;33(11):1047–63.
Bogoyevitch MA, Boehm I, Oakley A, Ketterman AJ, Barr RK. Targeting the JNK MAPK cascade for inhibition: basic science and therapeutic potential. Biochim Biophys Acta. 2004;1697(1–2):89–101.
Cano E, Mahadevan LC. Parallel signal processing among mammalian MAPKs. Trends Biochem Sci. 1995;20(3):117–22.
We thank Hai Shanzhao and Weifan Yao for their great help in our experiment.
This work was supported by grants from National Natural Science Foundation of China (No.81703560, 81573462, 81673475, U1608281, 81903658), Liaoning Revitalization Talents Program (No.XLYC1807201), Major Special S&T Projects in Liaoning Province [2019JH1/10300005], Liaoning Province Scientific Research Foundation (ZF2019038) and Shenyang S&T Projects (No. 19–109–4-09).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the first affiliated hospital of China Medical University (Approval number: AF-SOP-07-1.1-01). All animal experiments were approved by the Institutional Animal Care and Use Committee of China Medical University (Approval number: CMU2020250).
Consent for publication
All contributing authors agree to the publication of this article.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
CBR3-AS1 is related to breast cancer drug resistance and poor prognosis. A PCA analysis of the microarray of MCF-7/ADR cells and MCF-7/ADR cells. B The correlation between the expression of lncRNAs and the drug resistance of breast cancer cells to ADR. C Kaplan–Meier analysis were performed on patients of the expression of 5 lncRNAs in GSE20685.
Construction of the cell models that interferes with the expression of CBR3-AS1. A After transfection of CBR3-AS1 plasmid, the expression of CBR3-AS1 in MCF-7 cells. B After transfection of CBR3-AS1 siRNAs, the expression of CBR3-AS1 in MCF-7/ADR cells. C The relative viability of the MCF-7 cells in cisplatin was detected by CCK-8 assays 48 h after transfection with Control and OE-CBR3-AS1. D The relative viability of the MCF-7 cells in paclitaxel was detected by CCK-8 assays 48 h after transfection with Control and OE-CBR3-AS1. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
The relationship between CBR3-AS1, miR-25-3p and JNK1/MEK4. A Diana lncbase predicted the miRNAs combined with CBR3-AS1. B Subcellular localization of CBR3-AS1 predicted by the lnclocator website. C-D JNK1/MEK4 expression was measured by qRT-PCR after inhibited miR-25-3p in MCF-7 cells. E-F JNK1/MEK4 expression was measured by western blot and qRT-PCR after overexpressed miR-25-3p in MCF-7/ADR cells. G Linear correlation pattern showing a positive relationship between the expression of CBR3-AS1 and JNK1/MEK4. H Linear correlation pattern showing a negative relationship between the expression of miR-25-3p and JNK1/MEK4. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
The sequences for primers used in the study.
Antibodies used for IHC and WB.
The mRNAs banding to miR-25-3p.
About this article
Cite this article
Zhang, M., Wang, Y., Jiang, L. et al. LncRNA CBR3-AS1 regulates of breast cancer drug sensitivity as a competing endogenous RNA through the JNK1/MEK4-mediated MAPK signal pathway. J Exp Clin Cancer Res 40, 41 (2021). https://doi.org/10.1186/s13046-021-01844-7
- Breast cancer
- MAPK signaling pathway