Next Article in Journal
Molecular Insights of MAP4K4 Signaling in Inflammatory and Malignant Diseases
Previous Article in Journal
Improvement of Lung NET Management through Standardized Care—A Swiss Nationwide Observational Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SRSF5 Regulates the Expression of BQ323636.1 to Modulate Tamoxifen Resistance in ER-Positive Breast Cancer

1
Department of Pathology, School of Clinical Medicine, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
2
Department of Clinical Oncology, School of Clinical Medicine, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
3
Department of Clinical Oncology, Queen Mary Hospital, Hong Kong SAR, China
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(8), 2271; https://doi.org/10.3390/cancers15082271
Submission received: 31 January 2023 / Revised: 11 April 2023 / Accepted: 11 April 2023 / Published: 13 April 2023

Abstract

:

Simple Summary

Overexpression of BQ can confer tamoxifen resistance in ER +ve breast cancer. BQ is a variant of NCOR2 spliced through the exclusion of exon 11. SRSF5 is found to be involved in regulating the splicing of NCOR2. A low expression of SRSF5 enhances the chance of excluding exon 11 from mRNA, leading to the production of mRNA translating for BQ. Therefore, a low expression of SRSF5 can enhance BQ expression and, thus, tamoxifen resistance. SRPK1 can alter the activity of SRSF5 for BQ expression. Inhibition of SRPK1 can reduce TAM resistance by suppressing BQ expression in ER +ve breast cancer. Therefore, targeting the upstream pathway of BQ will be a possible strategy for reversing TAM resistance.

Abstract

About 70% of breast cancer patients are oestrogen receptor-positive (ER +ve). Adjuvant endocrine therapy using tamoxifen (TAM) is an effective approach for preventing local recurrence and metastasis. However, around half of the patients will eventually develop resistance. Overexpression of BQ323636.1 (BQ) is one of the mechanisms that confer TAM resistance. BQ is an alternative splice variant of NCOR2. The inclusion of exon 11 generates mRNA for NCOR2, while the exclusion of exon 11 produces mRNA for BQ. The expression of SRSF5 is low in TAM-resistant breast cancer cells. Modulation of SRSF5 can affect the alternative splicing of NCOR2 to produce BQ. In vitro and in vivo studies confirmed that the knockdown of SRSF5 enhanced BQ expression, and conferred TAM resistance; in contrast, SRSF5 overexpression reduced BQ expression and, thus, reversed TAM resistance. Clinical investigation using a tissue microarray confirmed the inverse correlation of SRSF5 and BQ. Low SRSF5 expression was associated with TAM resistance, local recurrence and metastasis. Survival analyses showed that low SRSF5 expression was associated with poorer prognosis. We showed that SRPK1 can interact with SRSF5 to phosphorylate it. Inhibition of SRPK1 by a small inhibitor, SRPKIN-1, suppressed the phosphorylation of SRSF5. This enhanced the proportion of SRSF5 interacting with exon 11 of NCOR2, reducing the production of BQ mRNA. As expected, SRPKIN-1 reduced TAM resistance. Our study confirms that SRSF5 is essential for BQ expression. Modulating the activity of SRSF5 in ER +ve breast cancer will be a potential approach to combating TAM resistance.

1. Introduction

Breast cancer is a heterogeneous disease with distinctive molecular subtypes. Based on the information from WHO, in 2020, there were 2.3 million women diagnosed with breast cancer, with 685,000 deaths globally, making it the world’s most prevalent cancer (https://www.who.int/news-room/fact-sheets/detail/breast-cancer; accessed on 27 June 2022). Many therapeutic options are now available, which are dependent upon the biological subtyping of breast cancers. Cancers that express oestrogen receptor α (ER +ve) and/or the progesterone receptor (PR +ve) are likely to respond to endocrine therapies such as tamoxifen (TAM). About 70% of breast cancers belong to this class. TAM is usually taken orally for 5–10 years, and can reduce the chance of the cancer recurrence by nearly half [1]. ER-driven signalling pathways are the dominant mechanism for driving the development and progression of ER +ve breast cancer. Therefore, suppressing ER-mediated signalling can effectively suppress breast cancer. TAM is an antagonist of ER. Once TAM binds to the ligand-binding domain of ER, it can prevent ER activation and therefore block ER activity.
Although TAM can effectively prevent breast cancer from recurrence and metastasis in ER +ve breast cancer, half of the patients will eventually develop resistance [2]. Endocrine-resistant patients represent up to 25% of all breast cancer patients, presenting a substantial clinical challenge. Possible mechanisms that confer TAM resistance have been intensely studied. These include the hypermethylation of the ER promoter region, leading to the downregulation of ER expression [3], hyperactivation of signalling cascades such as the PI3K/AKT/mTOR signalling pathway [4], the involvement of SOX9 and SOX2 transcription factors [5,6], and MYC [7]. Inhibiting the PI3K/AKT pathway with drugs targeting PI3K, mTOR, or AKT may improve the effectiveness of chemotherapy and endocrine therapy for breast cancer patients [8,9]. However, given the complexity of the pathway, pathway inhibition at any level may elicit compensatory mechanisms, limiting its effectiveness [10].
Alternative splicing (AS) can generate distinct mRNA isoforms from a single gene, and is shown to play a vital role in cancer progression and therapy [11]. We previously identified BQ323636.1 (BQ), a novel splice variant of NCOR2 that induces TAM resistance [12] and generates a monoclonal antibody that recognizes the unique epitope for the BQ protein (US patent no. US 10,823,735) [13]. BQ results in the exclusion of exon 11 of NCOR2, with exons 10 and 12 joined together to produce mRNA for BQ synthesis. We confirmed that BQ could induce TAM resistance through multiple signalling mechanisms such as HIF-1α signalling [14], the IL-6/STAT3 axis [14], and the IL-8/CXCR1 axis [15]. TAM resistance was reduced by reducing the expression of BQ [16]. Therefore, identifying upstream modulators of BQ would be essential to uncover novel mechanisms for TAM resistance in breast cancer.
Most precursor mRNAs are subject to alternative splicing through which gene expression can be regulated and the diversity of the proteome can be increased [17,18]. Splicing factors (SFs) that bind to specific sequences in an exon can regulate the choice of splice sites [19]. These SFs can act as splice enhancers to enhance the inclusion of an exon or as splice silencers which mediate the exclusion of an exon [20]. The serine/arginine-rich splicing factors (SRSF1-12) are a family of SFs that contribute significantly to regulating pre-mRNA alternative splicing [21]. SRSFs have been demonstrated to be associated with cancer [22,23]. SRSF1 overexpression has been shown to cause an oncogenic effect in various cancers [24]. SRSF2 knockdown has been shown to enhance the expression of GCH1-L and of STK39-L, which are the isoforms involved in cell cycle regulation and DNA repair, to cause an oncogenic effect [25]. SRSF3 overexpression has been shown to enhance the expression of the interleukin enhancer-binding factor 3 (ILF3) isoform 1 and isoform 2, which could promote cell proliferation in lung cancer [26]. Overexpression of SRSF5 has been shown to up-regulate the expression of oncogenic MCM2 and MCM4 [27]. The activities of SRSFs are regulated by phosphorylation mediated by several kinases, such as the serine-arginine protein kinases 1, 2, and 3 (SRPK1-3) [28], SRPK1 being the most well-documented. Aberrant SRPK1 expression has been documented in cancers, breast cancer included [29,30]. In addition, it has been reported that SRPK1 could modulate cisplatin resistance [31] and metastasis [32] in breast cancer. Therefore, the SRPK–SRSF regulatory axis should be considered important in cancer development and drug resistance.
TAM resistance remains a significant obstacle that hinders breast cancer treatment. Identifying possible strategies to combat TAM resistance is urgently needed. In this study, we illustrate a novel strategy to achieve this. We confirm that the down-regulation of SRSF5 results in the skipping of exon 11 of NCOR2, which thus favours BQ production, resulting in TAM resistance. Overexpression of SRSF5, on the other hand, reduces the expression of BQ and compromises TAM resistance. We also confirm that the function of SRSF5 is regulated by SRPK1. Inhibiting SRPK1, while decreasing SRSF5 phosphorylation, increases SRSF5 protein interaction with exon 11 of NCOR2. SRSF5 proteins favour the incorporation of exon 11, producing wild-type NCOR2 rather than BQ. A low expression of BQ compromises TAM resistance. Therefore, inhibiting SRPK1 in ER +ve breast cancer could reduce TAM resistance. Our investigation highlights the feasibility of targeting SRPK1 to reverse TAM resistance.

2. Materials and Methods

2.1. Mammalian Cell Culture, Transfection and Stable Cell Line Establishment

All cell lines were incubated in a tissue culture incubator with 5% CO2 at 37 °C. MCF-10A (a non-cancerous breast cell line), MCF-7 (tamoxifen-sensitive) and ZR-75 (tamoxifen-sensitive) were obtained from American Type Culture Collection (ATCC, Manassas, VA, USA). LCC2 is a tamoxifen-resistant cell line derived from MCF-7, while AK47 is a tamoxifen-resistant cell line derived from ZR-75. LCC2 and AK-47 were kindly provided by Professor Robert Clarke from Georgetown University Medical School, Washington, DC, USA. All cell lines were re-authenticated by short tandem repeat profiling [13].These cell lines were used in our previous studies [13,14,15,16]. MCF-10A was maintained in DMEM/F12 (11330032; Thermo Fisher Scientific, Waltham, MA, USA) supplemented with horse serum (5%; 16050122; Thermo Fisher Scientific, Waltham, MA, USA), EGF (20 ng/mL; PHG0313; Thermo Fisher Scientific, Waltham, MA, USA), hydrocortisone (0.5 mg/mL; H-0888; Sigma-Aldrich, St. Louis, MO, USA), cholera toxin (100 ng/mL; C-8052; Sigma-Aldrich, St. Louis, MO, USA), insulin (10 μg/mL; I-1882; Sigma-Aldrich, St. Louis, MO, USA), and P/S (1%; 10378016; Thermo Fisher Scientific, Waltham, MA, USA). MCF-7, ZR-75, LCC2, and AK-47 cells were cultured in DMEM (12100046; Thermo Fisher Scientific, Waltham, MA, USA) supplemented with FBS (10%; 26140079; Thermo Fisher Scientific, Waltham, MA, USA), and P/S (1%; 10378016; Thermo Fisher Scientific, Waltham, MA, USA). Mycoplasma examination was performed to confirm that all cell cultures were mycoplasma-free. The assay was conducted by the Core Facility provided by the Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China.
Lipofectamine 2000 (11668019; Thermo Fisher Scientific, Waltham, MA, USA) and Oligofectamine (12252011; Thermo Fisher Scientific, Waltham, MA, USA) were employed for plasmid and siRNA transfection, respectively, in accordance to the manufacturer’s instructions. An amount of 1 μg/mL of puromycin (A1113803; Thermo Fisher Scientific, Waltham, MA, USA) or 500 μg/mL of G418 (10131027; Thermo Fisher Scientific, Waltham, MA, USA) was employed for stable cell line selection. The selection process took 6 weeks.

2.2. Chemicals, siRNA, shRNA and Plasmids

The substances 4-hydroxytamoxifen (TAM; S7827; Selleck Chemicals LLC, Houston, TX, USA), SRPKIN-1 (HY-116856; MedChemExpress LLC; Monmouth Junction, NJ, USA) and TG003 (HY-15338; MedChemExpress LLC; Monmouth Junction, NJ, USA) were purchased and dissolved in DMSO. SRSF3-targeting siRNA (siSRSF3; sc-38338) and non-targeting siRNA (siCtrl; sc-37007) were obtained from Santa Cruz Biotechnology, Dallas, TX, USA. SRSF5 Human shRNA Plasmid Kit (TL309487; OriGene Technologies, Inc., Rockville, MD, USA) was used. SRSF5 Human Tagged ORF Clone (RC207271) was purchased from OriGene Technologies, Inc., Rockville, MD, USA. pCMV6-Entry Mammalian Expression Vector (PS100001; OriGene Technologies, Inc., Rockville, MD, USA) was used for control overexpression.

2.3. Real-Time qPCR and PCR-Based Splicing Assay

TRIzol was employed to extract the total RNA from mammalian cells (15596026; ThermoFisher, Waltham, MA, USA). cDNA was synthesised by PrimeScriptTM RT Master Mix (RR036; Takara Biomedical Technology Co., Ltd., Changping District, Beijing, China). An amount of 0.5 µg of RNA was used for cDNA synthesis. cDNA was analysed by PCR using the ProFlex PCR system (Thermo Fisher Scientific, Waltham, MA, USA). qPCR was analysed using SYBR Green Master Mix (A25742; ThermoFisher, Waltham, MA, USA). qPCR was performed using the StepOne Real-Time PCR system (Thermo Fisher Scientific, Waltham, MA, USA). The relative gene expression level was determined by the ΔΔCT method with actin as an internal control. GoTaq® G2 DNA Polymerase (M784A; Promega, Madison, WI, USA) was used to detect the expressions of BQ (352 bp) and NCOR2 (443 bp) amplicons with primers forward-splicing NCOR2-BQ and reverse-splicing NCOR2-BQ (refer to Figure S1 for the details). The following PCR conditions (35 cycles) were used: 95 °C for 30 s, 54 °C for 15 s and 72 °C for 15 s. Actin was used as the loading control. The DNA sequences of the primers used are as shown in Table 1.

2.4. Cell Viability Assay

All experiments were performed in triplicate. A total of 5000 cells were seeded onto 96-well plates, and cell viability was determined by Cell Counting Kit-8 (C0038; Beyotime, Shanghai, China) and an MTT assay (M6494; Thermo Fisher Scientific, Waltham, MA, USA). The microplate reader Infinite F200 (Tecan, Seestrasse, Switzerland) was employed to determine the absorbance at 570 nm. The reference wavelength was 650 nm. For the clonogenic assay, 2000 cells were seeded onto 2-well plates, with 0.01% crystal violet in PBS (C0775; Sigma, St. Louis, MO, USA) being used to stain the cell colonies. A colony was defined to contain an aggregate of more than 50 cells. For the parental cell line MCF-7, since the colonies were easily distinguishable, the clonogenicity (%) of MCF-7 shCtrl, MCF-7 shSRSF5.1 and MCF-7 shSRSF5.2 was determined. On the other hand, for the parental ZR-75 cell line, colonies were not easily distinguished, hence the solvent extraction method was applied to determine the cell viability of ZR-75 shCtrl, ZR-75 shSRSF5.1 and ZR-75 shSRSF5.2. An amount of 500 μL of isopropanol was added to extract the crystal violet to form a purple solution. An amount of 100 μL of the solution was transferred to a 96-well plate. Triplicate samples were used. The absorbance values at 570 nm and 650 nm were determined using the microplate reader Infinite F200 (Tecan, Seestrasse, Switzerland).

2.5. Western Blot and Co-Immunoprecipiation

Proteins from the cells were extracted using a RIPA buffer (50 mM Tris-Cl pH7.4, 150 mM NaCl, 1% NP-40, 0.1% SDS, 1% sodium deoxycholate) supplemented with phenylmethanesulfonyl fluoride (PMSF; P7626; Sigma, St. Louis, MO, USA) and PhosSTOP EASYPack tablets (4906837001; Sigma, St. Louis, MO, USA). A protein assay (5000112; Bio-Rad, Hercules, CA, USA) was employed to determine protein concentrations. An amount of 20 μg of proteins was analysed using SDS-PAGE and transferred to PVDF (1620177; Bio-Rad, Hercules, CA, USA). The following primary antibodies were employed: anti-BQ (1:1000; D-12; Versitech Limited, Hong Kong SAR, China), anti-SRSF5 (1:1000; ab67175; Abcam; Cambridge, UK); anti-SRSF3 (1:1000; #51039; Cell Signaling Technology, Inc.; Danvers, MA, USA); anti-NCOR2 (1:1000; #62370; Cell Signaling Technology, Inc.; Danvers, MA, USA); anti-phosphoepitope SR proteins (MABE50; 1:1000; Merck, Rahway, NJ, USA); anti-GAPDH (1:10,000; sc-32233; Santa Cruz Biotechnology, Dallas, TX, USA); and anti-HSP90 (1:4000; #4874; Cell Signaling Technology, Inc.; Danvers, MA, USA). The secondary antibodies used were anti-rabbit HRP (1:5000; P0260; Agilent Dako, Santa Clara, CA, USA) and anti-mouse HRP (1:5000; P0447; Agilent Dako, Santa Clara, CA, USA). The experiments were repeated at least three times, and representative images were shown.
Cells were lysed in 200 μL of a Co-IP buffer (20 mM Tris-Cl, pH 7.4, 150 mM NaCl, 0.5% NP-40, and 10% glycerol) supplemented with phenylmethanesulfonyl fluoride (PMSF; P7626; Sigma, St. Louis, MO, USA) and PhosSTOP EASYPack tablets (4906837001; Sigma, St. Louis, MO, USA). An amount of 10 μL of the lysate was saved (input fraction). An amount of 90 μL of the lysate was incubated with anti-SRPK (1:200; sc-100443; Santa Cruz Biotechnology, Dallas, TX, USA) or anti-mouse IgG (1:200; X0931; Agilent Dako, Santa Clara, CA, USA) at 4 °C overnight. In total, 50 µL of Dynabead Protein A (10002D; Thermo Fisher Scientific, Waltham, MA, USA) was used. The beads were incubated with the immunoprecipitant at 4 °C for 2 h. The loaded beads were incubated in 1 mL of the Co-IP buffer 3 times at room temperature, followed by a 10 min wash each time. A total of 50 µL of the sample buffer (2 × SDS) was used to incubate the loaded beads at 99 °C for 5 min. Candidate proteins were detected by Western blot. HRP-conjugated Protein A (1:4000; 18–160; Sigma-Aldrich, St. Louis, MO, USA) was used as the secondary antibody.

2.6. RNA–Protein Interaction

1 × 107 were resuspended in a 200 μL lysis buffer (with 20 mM HEPES, a pH of 7.4, 150 mM NaCl, 5 mM MgCl2, 0.5% Nonidet P-40, and a 40-unit RNaseOUT™ Recombinant Ribonuclease Inhibitor; 10777019; ThermoFisher, Waltham, MA, USA), and 10 µg/mL of yeast tRNA (AM7119; ThermoFisher, Waltham, MA, USA). The suspended cells were lysed by sonication (70% power; 30 s on/off cycle for 6 cycles) on ice using ultrasonic cell disruptor Branson Digital Sonifier (SFX 150; Branson Ultrasonics Corp; Brookfield CT, USA) with a microtip (109-122-1066; Branson Ultrasonics Corp; Brookfield, CT, USA). The RNA-protein complex was immunoprecipitated by anti-Myc (1:400; #2276; Cell Signaling Technology, Inc.; Danvers, MA, USA) or anti-rabbit IgG (1:400; X0903; Agilent Dako, Santa Clara, CA, USA) at 4 °C overnight. No antibody control was included. The immunoprecipitant was incubated with 50 μL of Protein A agarose beads (20333; ThermoFisher, Waltham, MA, USA) at 4 °C for 2 h. The beads were incubated with 1 mL of the lysis buffer three times with rotation at 4 °C, followed by 10 min washing each time. An amount of 1 mL of the Trizol reagent was employed to extract associated RNA. RT-qPCR was performed to detect the target mRNA. NCOR2-exon 11-F (5′-AGA AGA AGG TGG AGC GCA T-3′) and NCOR2-exon 11-R (5′-GCT GCT TGC GGA TCT CAG-3′) were used. The sequences were as follows. HSP70-3′-UTR-F: AGG CCA GGT TCT TAG CAC AT; HSP70-3′-UTR-R: GGC CAG TTG TAA AGG GTT GG; β-globin-3′-UTR-F: CTC GCT TTC TTG CTG TCC AA; β-globin-3′-UTR-R: CAA GGC CCT TCA TAA TAT CCC C. ∆CT = CTelute − CTinput and ∆∆CT = ∆CTanti-SRSF5 or IgG − ∆CTno Ab. Expression = 2−(∆∆CT).

2.7. Xenograft

The 1 × 107 cells of ZR-75-shCtrl, ZR-75-shSRSF5, LCC2-OE-Ctrl, LCC2-OE-SRSF5 were implanted into the mammary fatpad of nude mice (BALB/cAnN-nu) at 6 weeks of age. When tumours became palpable, the mice were randomly assigned into treatment and control groups. TAM dissolved in peanut oil or saline was used. Each mouse received TAM (0.5 mg) through subcutaneous injection twice per week. The treatment window was 6 weeks. The mice were administered TAM or saline 12 times. An electronic caliper was employed to measure tumour dimensions (the longest diameter (D) and the shortest parameter (d)). The volume (mm3) of tumour was calculated according to the formula D × d2/2. The procedures received official approval from the University of Hong Kong Committee on the Use of Live Animals in Teaching and Research (5103-19).

2.8. Tissue Microarray and Immunohistochemistry

Approval for the clinical study was obtained from the Institutional Review Board of HKU and HA Hong Kong West Cluster (identification number: HKU/HA HKW IRB, no. UW 08-147). Clinical information about the patients was retrieved from the database of Queen Mary Hospital, Hong Kong SAR. Representative histological sections were reviewed by pathologists. Regions from each of the tumours were selected and used for the construction of the tissue microarray (TMA), totalling 137 cases (Table 2). Each tumour was constructed in triplicate, the average score being taken for each case. TMA slides were deparaffinised with xylene and rehydrated in ethanol. Antigen retrieval involved heating them in the Tris-EDTA buffer at pH 8.0 using a pressure cooker for 10 min. The slides were then immersed in 3% H2O2 solution, and then rinsed twice with TBS. Incubation in primary monoclonal antibodies, anti-SRSF5 (SRP40) (1:50; RN082PW; Medical and Biological Laboratories Co. Ltd., Tokyo, Japan) or monoclonal BQ323636.1 (1:50; D-12, Versitech Ltd., Hong Kong SAR, China) was carried out at 4 °C for 20 h. The slides were rinsed in TBST and then incubated with secondary anti-rabbit (K4003; Agilent Dako, Santa Clara, CA, USA) or anti-mouse (K4001; Agilent Dako, Santa Clara, CA, USA) antibodies for 30 min at room temperature. A further wash with TBST was performed before chromogen DAB/substrate reagent incubation for 10 min. The images of TMA slides were captured using the Aperio ScanScope system (Leica Biosystems, Wetzlar, Germany) for further evaluation. Scoring was performed by two independent individuals. The H-score was used for SRSF5 and BQ323636.1 nuclear expression, calculated as follows: 1 × % of nuclei stained at low intensity + 2 × % of nuclei stained at moderate intensity + 3 × % of nuclei stained at high intensity. The median of the H-score was the threshold, which was 110 for the BQ nuclear score and 110 for the SRSF5 nuclear score.

2.9. Statistical Analysis

All numeric data were recorded, maintained and processed in Excel (Microsoft). Statistical analyses were performed using the internal functions of Prism5 (GraphPad) or SPSS25 (IBM). All results were shown as mean value ± SD obtained from at least three independent experiments in a bar chart format. For in vitro studies and a Students’ t-test were used to determine the significance between the two groups. A one-way ANOVA with Bonferroni’s post hoc test was employed for comparisons with multiple groups. For animal experiments, a two-way ANOVA with Bonferroni’s post hoc test was employed to determine the statistical significance between groups. For in vivo studies, a Mann–Whitney U test was employed to determine the significance between the clinical samples. TAM resistance was defined as ER +ve breast cancer patients who had been treated with TAM in the adjuvant setting but who eventually developed local cancer recurrence or metastases. All tests were two-sided. A chi-square test was used for hypothesis testing. Overall and disease-specific survival analyses were performed using Kaplan–Meier estimates followed by a log-rank test. The association between clinical–pathological parameters with nuclear BQ expression and with nuclear SRSF5 expression was determined by univariate and multivariate Cox-regression analyses, with relative risk (RR) and the 95% confidence intervals (CI) given. The Omnibus test in SPSS was employed to test the proportional hazards assumption. *, **, and *** represent p < 0.05, p < 0.01 and p < 0.001, respectively.

3. Results

3.1. SRSF5 Modulated the Expression of BQ

We previously showed that the overexpression of BQ323636.1 (BQ) is involved in the development of tamoxifen (TAM) resistance in ER +ve breast cancer [13]. Understanding the molecular mechanism through which BQ is overexpressed would help us gain insight into developing a novel strategy to combat TAM resistance. From the PCR-based splicing assay, we confirmed that BQ mRNA was enhanced in the two TAM-resistant cell lines LCC2 and AK-47 compared with their parental cell lines, MCF-7 and ZR-75 (Figure S1). Serine and arginine-rich splicing factors (SRSFs) are known to bind to the exonic splicing enhancers (ESS) located on exons, giving rise to the inclusion of that particular exon in the mature mRNA [33]. Hence, when the ESS is not bound with any SRSF, an exon will be excluded from the mature mRNA. We hypothesized that the down-regulation of particular SRSFs would increase the chance of excluding exon 11 of NCOR2, and thus generate mRNA for BQ synthesis. The TAM-resistant cell lines LCC2 and AK-47 have high endogenous BQ expression. We compared the expression of SRSF1-12 between TAM-resistant LCC2 and AK-47 cell lines of high BQ expression, with TAM-sensitive parental cell lines MCF-7 and ZR-75, which have low BQ expression. Through qPCR, we found that SRSF3 and SRSF5 were consistently down-regulated in the TAM-resistant cell lines (Figure 1A and Figure S2). We confirmed by Western blot that the protein expression of both SRSF3 and SRSF5 was reduced in LCC2 and AK-47 (Figure 1B and Figure S3). We employed shRNA to generate two independent clones (shSRSF5.1 and shSRSF5.2) with a SRSF5 stable knockdown using MCF-7 and ZR-75 cell lines. The knockdown efficiency was confirmed by RT-qPCR (Figure S4). Subsequently, we confirmed that only the downregulation of SRSF5 could modulate the expression of BQ in MCF-7 and ZR-75 (Figure 1C and Figure S5) but not SRSF3 (Figures S6 and S7). In addition, we found that the overexpression of SRSF5 could reduce the expression of BQ in LCC2 and AK-47 cell lines (Figure 1D and Figure S8). Our results, therefore, indicate that SRSF5 can modulate BQ expression.

3.2. SRSF5 Bound to Exon 11 of NCOR2 to Favour the Synthesis of Wild-Type NCOR2 mRNA

To produce BQ, skipping exon 11 of NCOR2 is necessary, as verified from clinical samples [12]. We proposed that SRSF5 binds to exon 11, mediating the splicing process to include this exon in the resultant mRNA, producing wild-type NCOR2. A low expression of SRSF5 would reduce the chance of exon 11 being incorporated into NCOR2 mRNA; thus, BQ would be produced instead of NCOR2 (Figure 2A). First, we confirmed that SRSF5 could bind to exon 11 of NCOR2 in MCF-7 and ZR-75 by studying RNA–protein complexes with NCOR2 exon 11-F and NCOR2 exon 11-R. The cells were transfected with a pCMV-myc-SRSF5 mammalian expression plasmid, and immunoprecipitation of the SRSF5/RNA complex with an anti-myc antibody was carried out. The presence of exon 11 was determined by qPCR (Figure 2B). Furthermore, we found that the ectopic SRSF5 could not interact with 3′-UTR of HSP70 and β-globin (Figure S9). The results, therefore, indicate that the interaction between SRSF5 and exon 11 of NCOR2 should be specific. Next, we found that the knockdown of SRSF5 could enhance the BQ mRNA level while reducing NCOR2 mRNA in MCF-7 and ZR-75, as revealed by the qPCR (Figure 2C) with the use of BQ targeting primers, BQ-F and BQ-R, and NCOR2 targeting primers, NCOR2-F and NCOR2-R (Table 1). In addition, the splicing assay was performed using the primers NCOR2-BQ-splicing-F and NCOR2-BQ-splicing-R, which could amplify both NCOR2 and BQ (Figure S10). The opposite results were found when SRSF5 was overexpressed in LCC2 and AK-47 (Figure 2D and Figure S11). These findings suggest that SRSF5 expression can modulate the expression of NCOR2 and BQ in breast cancer.

3.3. SRSF5 Could Modulate TAM Response

Since BQ overexpression is responsible for TAM resistance [14], we hypothesised that the down-regulation of SRSF5 would induce TAM resistance in TAM-sensitive MCF-7 and ZR-75 cell lines. Indeed, the knockdown of SRSF5 in MCF-7 and ZR-75 made the cells more tolerant to TAM (Figure 3A). Similar results were obtained from a clonogenic assay (Figure 3B). Next, we employed the xenograft model. We found that the knockdown of SRSF5 could make the cells TAM-resistant, as revealed by the tumour volume (Figure 3C). Finally, we found that the expression of BQ in tumour samples from the SRSF5 knockdown group was significantly enhanced (Figure 3D and Figure S12).
Conversely, we also showed that overexpression of SRSF5 in LCC2 could reduce TAM resistance. As expected, the results from the cell viability (Figure 4A) and clonogenic assays (Figure 4B) suggested that SRSF5 overexpression could recover TAM sensitivity. Results from the nude mice model also showed similar results (Figure 4C). The expression of BQ was significantly reduced in tumours with SRSF5 overexpression (Figure 4D and Figure S13). Our findings thus also indicate that SRSF5 expression can modulate the TAM response in ER +ve breast cancer.

3.4. Clinical Significance of SRSF5 in ER +ve Breast Cancer

Immunohistochemistry was employed to determine SRSF5 and BQ expression in 137 breast cancer samples. The samples were assigned into low and high nuclear expression groups based on the median value obtained (Figure 5A). The expressions of SRSF5 and BQ were found to be negatively correlated (Spearman’s R value = −0.3387; p < 0.001), which is in keeping with our in vitro findings. Previously, we had already confirmed that high BQ expression was associated with TAM resistance [13]. As expected, in ER +ve cases, the low expression of SRSF5 was significantly associated with TAM resistance (p = 0.015; Figure 5C) and local cancer recurrence (p = 0.024; Figure 5D) and metastasis (p = 0.002; Figure 5E) as indicated by a chi-square test. A low expression of SRSF5 was also a poor prognostic factor for both overall (p < 0.001; Figure 5F) and disease-specific (p = 0.001; Figure 5G) survival, as revealed by a Kaplan–Meier survival analysis with the log-rank test. The univariate cox regression analysis also showed that cases with high SRSF5 expression were significantly associated with both reduced overall (p < 0.001; RR = 0.182; 95% CI 0.071, 0.466; Table 3) and disease-specific survival (p = 0.003; RR = 0.142; 95% CI 0.039, 0.519; Table 3), which remained significant in the multivariate cox-regression analysis; cases with high SRSF5 expression showed both reduced overall (p = 0.041; RR = 0.087; 95% CI 0.008, 0.901; Table 4) and disease-specific survival (p = 0.039; RR = 0.178; 95% CI 0.035, 0.920; Table 4). These findings suggest that low SRSF5 expression is a poor prognostic factor in ER +ve breast cancer.

3.5. Inhibition of SRPK1 Could Reduce BQ Expression and TAM Resistance by Altering SRSF5 Activity

Given the above findings, we proposed that TAM resistance could be reduced by altering the activity of SRSF5. Cdc2-like kinase 1 (CLK1) has been reported to regulate the activity of SRSF5 [34]. TG003 is an inhibitor of CLK1 [35]. We identified 50 nM of TG003 as the maximum non-lethal dosage in the non-cancerous breast cell line MCF-10A (Figure S14A). However, when we treated LCC2 and AK-47 with 50 nM of TG003, it could not reduce the mRNA expression of BQ and NCOR2 based on the results from qPCR (Figure S14B). As expected, this treatment could not reduce TAM resistance, as revealed by the cell viability assay (Figure S14C). Using the STRING database [36], we found that serine/arginine-rich splicing factor (SRSF) protein kinase-1 (SRPK1) would interact with SRSF5 (Figure S15). By co-immunoprecipitation, we confirmed that SRPK1 could interact with SRSF5 in LCC2 (Figure 6A). Next, we proposed that the inhibition of SRPK1 might alter the splicing activity of SRSF5 and, thus, BQ expression. We employed SRPKIN-1 to target SRPK1 [37]. We confirmed that 100 nM of SRPKIN-1 was the maximum non-lethal dosage on MCF-10A (Figure S16). Next, we found that 100 nM of SRPKIN-1 could reduce BQ expression while enhancing NCOR2 at both mRNA (Figure 6B) and protein (Figure 6C and Figure S17) levels, with SRSF5 protein expression remained the same. Next, we evaluated the effect of SRPKIN-1 on SRSF5 activity at the phosphorylation level. Since no available antibody specifically recognises phosphorylated SRSF5, we needed to perform an indirect method using co-immunoprecipitation and an anti-phosphoepitope SR protein antibody. This antibody can detect phosphorylated SRSF proteins but is not specific to SRSF5. Hence, immunoprecipitation was needed to isolate the Myc-tagged SRSF5 and thus detect the degree of SRSF5 phosphorylation. The results showed that SRPKIN-1 could reduce the degree of phosphorylation on SRSF5 (Figure 6D and Figure S18).
Reduced SRSF5 phosphorylation resulted in the proportion of SRSF5 interacting with exon 11 of NCOR2 to be increased (Figure 6E), in keeping with the decreased BQ and increased NCOR2, mRNA and protein expression demonstrated on SRPK1 treatment (Figure 6B,C). As our earlier findings had shown, increased binding of SRSF5 to exon 11 would result in preferential NCOR2 expression over BQ expression, leading to decreased TAM resistance. In keeping with this, we also confirmed that treatment with SRPKIN-1 could reverse TAM resistance (Figure 6F). Our findings, therefore, indicate that modulation of SRSF5 activity would be a potential approach to reduce the expression of BQ and, thus, reverse TAM resistance.

4. Discussion

Tamoxifen (TAM) resistance is an obstacle to adjuvant therapy for breast cancer treatment. Although ER +ve breast cancer patients initially show a positive response to TAM, up to 50% of patients eventually develop resistance [10]. Loss of the ER during the treatment window is one of the reasons for not responding to TAM. It has been demonstrated that the CpG island of the ER promoter region is hypermethylated in TAM-resistant breast cancer cells [38]. In addition to DNA methylation, histone modification and miRNA have been implicated in TAM resistance [39,40]. Therefore, epigenetics contribute significantly to the development of TAM resistance. However, it would be less feasible to develop a therapy to combat TAM resistance through interfering with DNA methylation, histone modification or miRNA as this would systematically impact healthy cells. Cytotoxicity cannot be neglected.
Receptor-tyrosine kinases are a class of enzyme-linked receptors, which include the epidermal growth factor receptor, insulin and insulin-like growth factor-1 (IGF1) receptor, vascular endothelial growth factor receptor, etc., which have been found to confer tamoxifen resistance [4]. HER2, an EGF receptor, can activate the PI3K/AKT signalling pathway, which is one of the critical mechanisms of TAM resistance [41]. The IGF1 receptor activating PAK2 and VEGF overexpression, which leads to the activation of MAPK/ERK signalling, modulates the development of tamoxifen resistance [42,43,44,45]. These pathways can also crosstalk with each other to make the mechanisms of TAM resistance more complex. Therefore, targeting a single pathway might not effectively combat TAM resistance. However, the cytotoxicity would be too strong to manage when targeting these pathways altogether.
Transcriptome changes have also been shown to contribute to TAM resistance [46]. The transcriptome is governed by transcription factors [47], alternative splicing [48], non-coding RNAs [49] and the epigenome [50]. Alterations in the transcriptome would impact the proteomics [51], with altered protein expression affecting signalling mechanisms. Splicing increases the complexity and diversity of the genome, as a single pre-mRNA can generate various alternative spliced isoforms and thus translate to different protein isoforms [52]. Dysregulation of alternative splicing can lead to aberrant protein isoforms, which may contribute to cancer development [53]. Changing the splicing pattern of the spleen tyrosine kinase gene (SYK) has been shown to impair cell cycle progression [54]. Studies have demonstrated that alternative splicing isoforms could be novel therapeutic candidates for human diseases such as breast cancer [55,56,57,58].
Our previous studies identified a novel splice isoform of NCOR2, called BQ323636.1 (BQ), which is associated with TAM resistance in ER +ve breast cancer [12]. BQ has been shown to modulate various molecular mechanisms, such as HIF-1α-, NF-ĸB-, and AR-driven pathways, to confer TAM resistance [14,15,16]. Knockdown of BQ can compromise TAM resistance, highlighting the importance of BQ in breast cancer, especially in TAM resistance [15]. In this study, we identified that SRSF5 could modulate the alternative splicing of NCOR2 (Figure 1). SRSF5 was found to interact with exon 11 of NCOR2 (Figure 2). The presence of SRSF5 favours the inclusion of the exon 11, resulting in mRNA translating for wild-type NCOR2. Reduced expression of SRSF5 favours exon 11 exclusion, resulting in mRNA translating for BQ. As a functional consequence, TAM resistance would result from low SRSF5 expression. In our xenograft model, we confirmed that the expression of SRSF5 could modulate the TAM response (Figure 3 and Figure 4), indicating the functional significance of SRSF5 in TAM resistance. The negative correlation of SRSF5 expression with BQ overexpression was further examined using the primary cancer samples in the tissue microarray (TMA) to validate this notion. The protein expression of SRSF5 and BQ was determined by immunohistochemistry (IHC). The results confirmed that BQ and SRSF5 were negatively correlated (Figure 5B). As expected, a low expression of SRSF5 was associated with TAM resistance (Figure 5C), local cancer recurrence (Figure 5D) and metastasis (Figure 5E) since ineffective TAM treatment can result in cancer recurrence and metastasis [59]. Our results validated our hypothesis and suggested that low SRSF5 expression should be of clinical significance. Survival analyses confirmed that patients with low SRSF5 expression had poorer outcomes (Figure 5F,G), as also shown by univariate and multivariate cox regression analyses. Conversely, ER +ve breast cancer patients with high SRSF5 expression had a reduced risk ratio (Table 3 and Table 4). These findings confirm that SRSF5 should be involved in modulating NCOR2 production. Reducing SRSF5 expression would favour the production of BQ and induce the development of TAM resistance.
Previously, through a PCR array, we found that BQ overexpression can activate various cancer-related pathways [14]. BQ is an N-terminal part of NCOR2. NCOR2 binds to itself in an anti-parallel fashion, forming a homodimer and giving rise to a functional gene repressor complex [60]. Through mechanistic study, we showed that BQ can interact with NCOR2 to disrupt the repressor activity of NCOR2. In ER +ve breast cancer, NCOR2 represses the activity of ER. However, this repressive activity is compromised in the presence of BQ overexpression. BQ can modulate the interaction of NCOR2 with other transcription factors [14,15,16,61]. Therefore, it would not be surprising that BQ can modulate different mechanisms to confer TAM resistance. Thus, by reducing the expression of BQ, it might be possible to suppress various pathways to overcome TAM resistance. Checkpoint kinase 2 (CHK2) has been shown to regulate the protein stability of BQ via post-translational modification; suppressing CHK2 could reduce the expression of BQ and thus TAM resistance [62]. This study confirms that SRSF5 is an upstream mediator of BQ.
SRSF5 has been shown to modulate apoptosis in lung cancer by controlling the alternative splicing of CCAR1 [63]. Chen et al. [63] demonstrated that acetylation mediated by TIP60 and ubiquitination mediated by SMURF1 could modulate the expression and activity of SRSF5. Cdc2-like kinase 1 (CLK1) has been shown to enhance the phosphorylation of SRSF5 on serine 250 [34]. SRSF5 phosphorylation could modulate the splicing of METTL14 by inhibiting exon skipping while it promoted Cyclin L2 exon skipping [34]. Thus it is also possible that the inhibition of SRSF5 modulators such as CLK1 might alter BQ expression. However, we found that the inhibition of CLK1 was unable to affect BQ expression (Figure S12B). Fortunately, SRPK1 was found to interact with SRSF5 (Figure 6A and Figure S13). Inhibition of SRPK1 by SRPKIN-1 could reduce BQ expression but enhance NCOR2 expression (Figure 6B,C). Our results suggest that SRPKIN-1 could reduce the degree of phosphorylation on SRSF5 (Figure 6D) and enhance the proportion of SRSF5 interacting with exon 11 of NCOR2 (Figure 6E). This would favour the production of mRNA for NCOR2 expression. Therefore, as illustrated by Western blot, SRPKIN-1 could reduce the protein level of BQ (Figure 6C). BQ overexpression leads to TAM resistance [12,13]. Therefore, SRPKIN-1 could reverse TAM resistance by inhibiting the exon splicing event, leading to a reduced production of BQ mRNA.
Targeting epigenetic mechanisms or various kinase signalling pathways is not a practical approach to combat TAM resistance in breast cancer. Our work indicates that BQ was differentially expressed to induce TAM resistance through multiple mechanisms [13,14]. Reducing BQ expression could compromise TAM resistance [15,62], highlighting the importance of BQ in breast cancer. Therefore, any approach that can reduce the expression of BQ may be a potential method to reduce TAM resistance. In the current study, we confirmed that SRSF5 is an upstream factor that governs BQ expression. Altering either SRSF5 expression or its phosphorylation could modulate BQ expression and, thus, TAM resistance. These are novel approaches for developing strategies for reducing TAM resistance. As an exploratory study, we have confirmed that SRPKIN-1 could do so. It is possible that other types of post-translational modification of SRSF5 by different modulators can give rise to similar effects. Identifying modulators that would suppress alternative splicing for BQ production would help identify novel targets for combating TAM resistance.

5. Conclusions

Our results indicated that BQ expression was governed by SRSF5 expression. SRSF5 could bind to exon 11 of NCOR2 to ensure the inclusion of this exon in the mature mRNA. When SRSF5 was reduced, the chance of excluding exon 11 would be enhanced, favouring the product of mRNA for BQ synthesis. TAM response would be altered by SRSF5 expression. An in vivo study indicated that a low expression of SRSF5 was associated with high BQ expression, TAM resistance, local recurrence and metastasis. Furthermore, ER +ve breast cancer patients with low SRSF5 would have poorer overall and disease-specific survival outcomes.

6. Patents

Anti-BQ antibody (D12; US patent no. US 10,823,735).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15082271/s1, Figure S1: PCR-based splicing assay of NCOR2 and BQ. The DNA sequence of exons 10-12 of NCOR2 is shown. Using NCOR2-BQ splicing forward and reverse primers would yield two different PCR products with different product sizes, 443 bp and 352 bp for NCOR2 and BQ amplicons, respectively. Actin was the loading control. Figure S2: Expression levels of SRSF1-12. qPCR was performed. ∆∆CT method was used and MCF-10A was used as the reference. Actin was used as the internal control. Results were shown as mean ± sd from three independent experiments. Figure S3: Quantification of western blot results in Figure 1B. Results were shown as mean ± SD from four independent experiments. Students’ t-test was employed for statistical analysis. *** represents p < 0.001. Figure S4: Knockdown efficiency of shRNA targeting SRSF5. qPCR was employed to determine the expression of SRSF5. Actin was employed as the internal control. ∆∆CT method was used to evaluate the relative gene expression level. Results were shown as mean ± SD from three independent experiments. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance between shCtrl and shSRSF5.1 or shSRSF5.2 groups. *** represents p < 0.001. Figure S5: Quantification of western blot results in Figure 1C. Results were shown as mean ± SD from four independent experiments. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance between shCtrl and shSRSF5.1 or shSRSF5.2 groups. *** represents p < 0.001. Figure S6: The effect of SRSF3 siRNA. SRSF3 expression in (A) MCF-7 and (B) ZR-75. BQ expression in (C) MCF7 and (D) ZR-75. NCOR2 expression in (E) MCF-7 and (F) ZR-75. 25 nM of siSRSF3 and siCtrl was employed. qPCR was performed 48 h post-transfection. Actin was employed as the internal control. ∆∆CT method was used to evaluate the relative gene expression level. Results were shown as mean ± SD from four independent experiments. Students’ t-test was employed for statistical analysis. *** represents p < 0.001. Figure S7: SRSF3 knockdown could not affect BQ expression. 25 nM of siSRSF3 and siCtrl was employed. Western blot was performed 48 h post-transfection. HSP90 was the loading control. The bar charts are the quantification of the western blots. Results were shown as mean ± SD from four independent experiments. Students’ t-test was employed for statistical analysis. *** represents p < 0.001. Figure S8: Quantification of western blot results in Figure 1D. Results were shown as mean ± SD from four independent experiments. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance between untransfected and SRSF5 OE groups. *, and *** represent p < 0.05 and p < 0.001, respectively. Figure S9: SRSF5 could not bind to 3’-UTR of HSP70 and β-globin. MCF-7 and ZR-75 were transfected and overexpressed with myc-tagged SRSF5. Immunoprecipitation was performed using anti-SRSF5 or anti-IgG. RNA associated with SRSF5 was extracted by Trizol reagent. RT-qPCR was performed to detect the expression of 3’-UTR of HSP70 and β-globin. The expression of these candidates in the input fraction was used as reference. 3’-UTR of HSP70 and β-globin were absent in the immunoprecipitants. Results were shown as mean ± SD from three independent experiments. Figure S10: Knockdown of SRSF5 could enhance BQ expression but reduce NCOR2 expression. PCR-based splicing assay was performed on the control (shCtrl) and SRSF5 knockdown (shSRSF5.1 or shSRSF5.2). Actin was the loading control. Figure S11: Overexpression of SRSF5 could reduce BQ expression but enhance NCOR2 expression. PCRbased splicing assay was performed on LCC2 and AK-47 stably transfected with pCMV-Myc-SRSF5 (SRSF5 OE) or pCMV (Ctrl OE). Actin was the loading control. Figure S12: Quantification of western blot results in Figure 3D. Results were shown as mean ± SD from three independent tumours. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance compared with the control (ZR-75 shCtrl + Saline). *** represents p < 0.001. Figure S13: Quantification of western blot results in Figure 4D. Results were shown as mean ± SD from three independent tumours. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance compared with the control (LCC2 Ctrl OE + Saline). *** represents p < 0.001. Figure S14: TG003 could not affect BQ expression and TAM resistance. (A) the dose-dependent effect of CLK1 inhibitor TG003. MCF-10A was treated with different concentrations of TG003 for 96 h. Cell viability assay was performed. Results were shown as mean ± SD from four independent experiments. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance between DMSO and treatment groups. (B) the effect of TG003 on BQ and NCOR mRNA expressions in LCC2 and AK47. The cells were treated with 50 nM of TG003 for 72 h. qPCR was performed. Results were shown as mean ± SD from four independent experiments. Students’ t-test was employed for statistical analysis. (C) TG003 could not reverse TAM resistance. LCC2 and AK-47 were treated with 50 nM of TG003 and 5 μM of TAM for 96 h. Cell viability assay was performed. Results were shown as mean ± SD from four independent experiments. Two-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance. *, **, and *** represents p < 0.05, p < 0.01 and p < 0.001, respectively. Figure S15: Protein-protein interaction network of SRPK1. The information was obtained from STRING database. Figure S16: The dose-dependent effect of SRPK1 inhibitor SRPKIN-1. MCF-10A was treated with different concentrations of SRPKIN-1 for 96 h. Cell viability assay was performed. Results were shown as mean ± SD from four independent experiments. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance between DMSO and treatment groups. ** and *** represents p < 0.01 and p < 0.001, respectively. Figure S17: Quantification of western blot results in Figure 6C. Results were shown as mean ± SD from three independent experiments. Students’ t-test was employed for statistical analysis. *** represents p < 0.001. Figure S18: Quantification of western blot results in Figure 6D. SRSF5-Myc and pSR from elute fraction were normalised using SRSF5-Myc from input fraction. The relative expression of normalised SRSF5-Myc and pSR was determined. Results were shown as mean ± SD from three independent experiments. Oneway ANOVA with Bonferroni post hoc test was used to determine the statistical significance between untreated and SRPKIN-1 groups. * represents p < 0.05.

Author Contributions

Conceptualisation, H.T. and U.-S.K.; writing—original draft preparation, H.T.; writing—review and editing, H.T., C.-P.Y., M.-H.L., E.P.S.M. and U.-S.K.; in vitro and animal studies: H.T., C.-P.Y., M.-H.L., N.N.-C.F., E.P.S.M. and U.-S.K.; visualisation, H.T., C.-P.Y., M.-H.L. and U.-S.K.; histological examination: U.-S.K.; clinical investigation (samples preparation; IHC; data entry; data analysis), W.-L.C., S.-Y.C., E.P.S.M. and U.-S.K.; supervision, H.T. and U.-S.K.; project administration, E.P.S.M. and U.-S.K.; funding acquisition, H.T. and U.-S.K. All authors have read and agreed to the published version of the manuscript.

Funding

The General Research Fund (GRF) (grant number: 17113020) and the Health and Medical Research Fund (grants number: 07182026; 07182046) funded and supported the study.

Institutional Review Board Statement

The HKU Committee on the Use of Live Animals in Teaching and Research (5103-19) approved all the procedures for the animal study. The Institutional Review Board of The University of Hong Kong/Hospital Authority Hong Kong West Cluster (HKU/HA HKW IRB No. UW 08-147) approved the clinical study.

Informed Consent Statement

Not applicable.

Data Availability Statement

The materials and raw data will be made available to authorised researchers from the corresponding author upon reasonable request.

Acknowledgments

U.S. Khoo is an Ada Chan-Endowed Professor in Oncological Pathology. The Endowed Professorship is supported by the University of Hong Kong.

Conflicts of Interest

U.S Khoo is the inventor of the anti-BQ antibody (US patent no. US 10,823,735). Other authors declare no conflicts of interest.

References

  1. Yussof, I.; Mohd Tahir, N.A.; Hatah, E.; Mohamed Shah, N. Factors influencing five-year adherence to adjuvant endocrine therapy in breast cancer patients: A systematic review. Breast 2022, 62, 22–35. [Google Scholar] [CrossRef] [PubMed]
  2. Davies, C.; Godwin, J.; Gray, R.; Clarke, M.; Darby, S.; McGale, P.; Wang, Y.C.; Peto, R.; Pan, H.C.; Cutter, D.; et al. Relevance of breast cancer hormone receptors and other factors to the efficacy of adjuvant tamoxifen: Patient-level meta-analysis of randomised trials. Lancet 2011, 378, 771–784. [Google Scholar] [CrossRef] [Green Version]
  3. Ring, A.; Dowsett, M. Mechanisms of tamoxifen resistance. Endocr. Relat. Cancer 2004, 11, 643–658. [Google Scholar] [CrossRef] [PubMed]
  4. Mansouri, S.; Farahmand, L.; Teymourzadeh, A.; Majidzadeh-A, K. Clinical Evidence on the Magnitude of Change in Growth Pathway Activity in Relation to Tamoxifen Resistance is Required. Curr. Cancer Drug Targets 2018, 18, 668–676. [Google Scholar] [CrossRef]
  5. Jeselsohn, R.; Cornwell, M.; Pun, M.; Buchwalter, G.; Nguyen, M.; Bango, C.; Huang, Y.; Kuang, Y.N.; Paweletz, C.; Fu, X.Y.; et al. Embryonic transcription factor SOX9 drives breast cancer endocrine resistance. Proc. Natl. Acad. Sci. USA 2017, 114, E4482–E4491. [Google Scholar] [CrossRef] [Green Version]
  6. Piva, M.; Domenici, G.; Iriondo, O.; Rabano, M.; Simoes, B.M.; Comaills, V.; Barredo, I.; Lopez-Ruiz, J.A.; Zabalza, I.; Kypta, R.; et al. Sox2 promotes tamoxifen resistance in breast cancer cells. EMBO Mol. Med. 2014, 6, 66–79. [Google Scholar] [CrossRef]
  7. Yu, P.; Li, A.X.; Chen, X.S.; Tian, M.; Wang, H.Y.; Wang, X.L.; Zhang, Y.; Wang, K.S.; Cheng, Y. PKM2-c-Myc-Survivin Cascade Regulates the Cell Proliferation, Migration, and Tamoxifen Resistance in Breast Cancer. Front. Pharmacol. 2020, 11, 550469. [Google Scholar] [CrossRef] [PubMed]
  8. Lui, A.; New, J.; Ogony, J.; Thomas, S.; Lewis-Wambi, J. Everolimus downregulates estrogen receptor and induces autophagy in aromatase inhibitor-resistant breast cancer cells. BMC Cancer 2016, 16, 487. [Google Scholar] [CrossRef] [Green Version]
  9. Zhu, Y.H.; Liu, Y.J.; Zhang, C.; Chu, J.J.; Wu, Y.Q.; Li, Y.D.; Liu, J.Q.; Li, Q.; Li, S.Y.; Shi, Q.F.; et al. Tamoxifen-resistant breast cancer cells are resistant to DNA-damaging chemotherapy because of upregulated BARD1 and BRCA1. Nat. Commun. 2018, 9, 1595. [Google Scholar] [CrossRef] [Green Version]
  10. Yao, J.; Deng, K.; Huang, J.; Zeng, R.; Zuo, J. Progress in the Understanding of the Mechanism of Tamoxifen Resistance in Breast Cancer. Front. Pharmacol. 2020, 11, 592912. [Google Scholar] [CrossRef]
  11. Sette, C.; Ladomery, M.; Ghigna, C. Alternative splicing: Role in cancer development and progression. Int. J. Cell Biol. 2013, 2013, 421606. [Google Scholar] [CrossRef] [Green Version]
  12. Zhang, L.D.; Gong, C.; Lau, S.L.Y.; Yang, N.; Wong, O.G.W.; Cheung, A.N.Y.; Tsang, J.W.H.; Chan, K.Y.K.; Khoo, U.S. SpliceArray Profiling of Breast Cancer Reveals a Novel Variant of NCOR2/SMRT That Is Associated with Tamoxifen Resistance and Control of ER alpha Transcriptional Activity. Cancer Res. 2013, 73, 246–255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gong, C.; Man, E.P.S.; Tsoi, H.; Lee, T.K.W.; Lee, P.; Ma, S.T.; Wong, L.S.; Luk, M.Y.; Rakha, E.A.; Green, A.R.; et al. BQ323636.1, a Novel Splice Variant to NCOR2, as a Predictor for Tamoxifen-Resistant Breast Cancer. Clin. Cancer Res. 2018, 24, 3681–3691. [Google Scholar] [CrossRef] [Green Version]
  14. Tsoi, H.; Man, E.P.; Leung, M.H.; Mok, K.C.; Chau, K.M.; Wong, L.S.; Chan, W.L.; Chan, S.Y.; Luk, M.Y.; Cheng, C.N.; et al. KPNA1 regulates nuclear import of NCOR2 splice variant BQ323636.1 to confer tamoxifen resistance in breast cancer. Clin. Transl. Med. 2021, 11, e554. [Google Scholar] [CrossRef]
  15. Tsoi, H.; Shi, L.; Leung, M.-H.; Man, E.P.S.; So, Z.-Q.; Chan, W.-L.; Khoo, U.-S. Overexpression of BQ323636.1 Modulated AR/IL-8/CXCR1 Axis to Confer Tamoxifen Resistance in ER-Positive Breast Cancer. Life 2022, 12, 93. [Google Scholar] [CrossRef] [PubMed]
  16. Tsoi, H.; Man, E.P.S.; Chau, K.M.; Khoo, U.S. Targeting the IL-6/STAT3 Signalling Cascade to Reverse Tamoxifen Resistance in Estrogen Receptor Positive Breast Cancer. Cancers 2021, 13, 1511. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, E.T.; Sandberg, R.; Luo, S.; Khrebtukova, I.; Zhang, L.; Mayr, C.; Kingsmore, S.F.; Schroth, G.P.; Burge, C.B. Alternative isoform regulation in human tissue transcriptomes. Nature 2008, 456, 470–476. [Google Scholar] [CrossRef] [Green Version]
  18. Pan, Q.; Shai, O.; Lee, L.J.; Frey, B.J.; Blencowe, B.J. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat. Genet. 2008, 40, 1413–1415. [Google Scholar] [CrossRef]
  19. Dvinge, H. Regulation of alternative mRNA splicing: Old players and new perspectives. FEBS Lett. 2018, 592, 2987–3006. [Google Scholar] [CrossRef] [Green Version]
  20. Chao, Y.; Jiang, Y.; Zhong, M.; Wei, K.; Hu, C.; Qin, Y.; Zuo, Y.; Yang, L.; Shen, Z.; Zou, C. Regulatory roles and mechanisms of alternative RNA splicing in adipogenesis and human metabolic health. Cell Biosci. 2021, 11, 66. [Google Scholar] [CrossRef]
  21. Shepard, P.J.; Hertel, K.J. The SR protein family. Genome Biol. 2009, 10, 242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Cerasuolo, A.; Buonaguro, L.; Buonaguro, F.M.; Tornesello, M.L. The Role of RNA Splicing Factors in Cancer: Regulation of Viral and Human Gene Expression in Human Papillomavirus-Related Cervical Cancer. Front. Cell Dev. Biol. 2020, 8, 474. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, Y.; Qian, J.; Gu, C.; Yang, Y. Alternative splicing and cancer: A systematic review. Signal Transduct. Target. Ther. 2021, 6, 78. [Google Scholar] [CrossRef] [PubMed]
  24. Karni, R.; de Stanchina, E.; Lowe, S.W.; Sinha, R.; Mu, D.; Krainer, A.R. The gene encoding the splicing factor SF2/ASF is a proto-oncogene. Nat. Struct. Mol. Biol. 2007, 14, 185–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Luo, C.L.; Cheng, Y.M.; Liu, Y.G.; Chen, L.L.; Liu, L.N.; Wei, N.; Xie, Z.Q.; Wu, W.W.; Feng, Y. SRSF2 Regulates Alternative Splicing to Drive Hepatocellular Carcinoma Development. Cancer Res. 2017, 77, 1168–1178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Jia, R.; Ajiro, M.; Yu, L.L.; McCoy, P.; Zheng, Z.M. Oncogenic splicing factor SRSF3 regulates ILF3 alternative splicing to promote cancer cell proliferation and transformation. RNA 2019, 25, 630–644. [Google Scholar] [CrossRef] [PubMed]
  27. Yang, S.S.; Jia, R.; Bian, Z. SRSF5 functions as a novel oncogenic splicing factor and is upregulated by oncogene SRSF3 in oral squamous cell carcinoma. Biochim. Biophys. Acta Mol. Cell Res. 2018, 1865, 1161–1172. [Google Scholar] [CrossRef] [PubMed]
  28. Aubol, B.E.; Wu, G.W.; Keshwani, M.M.; Movassat, M.; Fattet, L.; Hertel, K.J.; Fu, X.D.; Adams, J.A. Release of SR Proteins from CLK1 by SRPK1: A Symbiotic Kinase System for Phosphorylation Control of Pre-mRNA Splicing. Mol. Cell 2016, 63, 218–228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. van Roosmalen, W.; Le Devedec, S.E.; Golani, O.; Smid, M.; Pulyakhina, I.; Timmermans, A.M.; Look, M.P.; Zi, D.; Pont, C.; de Graauw, M.; et al. Tumor cell migration screen identifies SRPK1 as breast cancer metastasis determinant. J. Clin. Investig. 2015, 125, 1648–1664. [Google Scholar] [CrossRef] [Green Version]
  30. Nikas, I.P.; Themistocleous, S.C.; Paschou, S.A.; Tsamis, K.I.; Ryu, H.S. Serine-Arginine Protein Kinase 1 (SRPK1) as a Prognostic Factor and Potential Therapeutic Target in Cancer: Current Evidence and Future Perspectives. Cells 2020, 9, 19. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, C.; Zhou, Z.H.; Subhramanyam, C.S.; Cao, Q.; Heng, Z.S.L.; Liu, W.; Fu, X.D.; Hu, Q.D. SRPK1 acetylation modulates alternative splicing to regulate cisplatin resistance in breast cancer cells. Commun. Biol. 2020, 3, 268. [Google Scholar] [CrossRef] [PubMed]
  32. Malvi, P.; Janostiak, R.; Chava, S.; Manrai, P.; Yoon, E.; Singh, K.; Harigopal, M.; Gupta, R.; Wajapeyee, N. LIMK2 promotes the metastatic progression of triple-negative breast cancer by activating SRPK1. Oncogenesis 2020, 9, 77. [Google Scholar] [CrossRef] [PubMed]
  33. Gehring, N.H.; Roignant, J.Y. Anything but Ordinary—Emerging Splicing Mechanisms in Eukaryotic Gene Regulation. Trends Genet. 2021, 37, 355–372. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, S.; Yang, C.; Wang, Z.W.; Hu, J.F.; Pan, J.J.; Liao, C.Y.; Zhang, J.Q.; Chen, J.Z.; Huang, Y.; Huang, L.; et al. CLK1/SRSF5 pathway induces aberrant exon skipping of METTL14 and Cyclin L2 and promotes growth and metastasis of pancreatic cancer. J. Hematol. Oncol. 2021, 14, 60. [Google Scholar] [CrossRef]
  35. Muraki, M.; Ohkawara, B.; Hosoya, T.; Onogi, H.; Koizumi, J.; Koizumi, T.; Sumi, K.; Yomoda, J.; Murray, M.V.; Kimura, H.; et al. Manipulation of alternative splicing by a newly developed inhibitor of Clks. J. Biol. Chem. 2004, 279, 24246–24254. [Google Scholar] [CrossRef] [Green Version]
  36. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef]
  37. Hatcher, J.M.; Wu, G.W.; Zeng, C.Y.; Zhu, J.; Meng, F.; Patel, S.; Wang, W.Q.; Ficarro, S.B.; Leggett, A.L.; Powell, C.E.; et al. SRPKIN-1: A Covalent SRPK1/2 Inhibitor that Potently Converts VEGF from Pro-angiogenic to Anti-angiogenic Isoform. Cell Chem. Biol. 2018, 25, 460–470.e6. [Google Scholar] [CrossRef] [Green Version]
  38. Williams, K.E.; Anderton, D.L.; Lee, M.P.; Pentecost, B.T.; Arcaro, K.F. High-density array analysis of DNA methylation in Tamoxifen-resistant breast cancer cell lines. Epigenetics 2014, 9, 297–307. [Google Scholar] [CrossRef] [Green Version]
  39. Wu, Y.M.; Zhang, Z.; Cenciarini, M.E.; Proietti, C.J.; Amasino, M.; Hong, T.; Yang, M.; Liao, Y.J.; Chiang, H.C.; Kaklamani, V.G.; et al. Tamoxifen Resistance in Breast Cancer Is Regulated by the EZH2-ER alpha-GREB1 Transcriptional Axis. Cancer Res. 2018, 78, 671–684. [Google Scholar] [CrossRef] [Green Version]
  40. Barazetti, J.F.; Jucoski, T.S.; Carvalho, T.M.; Veiga, R.N.; Kohler, A.F.; Baig, J.; Al Bizri, H.; Gradia, D.F.; Mader, S.; de Oliveira, J.C. From Micro to Long: Non-Coding RNAs in Tamoxifen Resistance of Breast Cancer Cells. Cancers 2021, 13, 3688. [Google Scholar] [CrossRef]
  41. Dong, C.; Wu, J.; Chen, Y.; Nie, J.N.; Chen, C.S. Activation of PI3K/AKT/mTOR Pathway Causes Drug Resistance in Breast Cancer. Front. Pharmacol. 2021, 12, 628690. [Google Scholar] [CrossRef] [PubMed]
  42. Apte, R.S.; Chen, D.S.; Ferrara, N. VEGF in Signaling and Disease: Beyond Discovery and Development. Cell 2019, 176, 1248–1264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Heckler, M.M.; Thakor, H.; Schafer, C.C.; Riggins, R.B. ERK/MAPK regulates ERR gamma expression, transcriptional activity and receptor-mediated tamoxifen resistance in ER plus breast cancer. FEBS J. 2014, 281, 2431–2442. [Google Scholar] [CrossRef] [Green Version]
  44. Peng, W.X.; Huang, J.G.; Yang, L.; Gong, A.H.; Mo, Y.Y. Linc-RoR promotes MAPK/ERK signaling and confers estrogen-independent growth of breast cancer. Mol. Cancer 2017, 16, 161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Zhang, Y.H.; Wester, L.; He, J.C.; Geiger, T.; Moerkens, M.; Siddappa, R.; Helmijr, J.A.; Timmermans, M.M.; Look, M.P.; van Deurzen, C.H.M.; et al. IGF1R signaling drives antiestrogen resistance through PAK2/PIX activation in luminal breast cancer. Oncogene 2018, 37, 1869–1884. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Rahem, S.M.; Epsi, N.J.; Coffman, F.D.; Mitrofanova, A. Genome-wide analysis of therapeutic response uncovers molecular pathways governing tamoxifen resistance in ER+ breast cancer. Ebiomedicine 2020, 61, 103047. [Google Scholar] [CrossRef] [PubMed]
  47. Lago, S.; Nadai, M.; Cernilogar, F.M.; Kazerani, M.; Moreno, H.D.; Schotta, G.; Richter, S.N. Promoter G-quadruplexes and transcription factors cooperate to shape the cell type-specific transcriptome. Nat. Commun. 2021, 12, 3885. [Google Scholar] [CrossRef]
  48. Vaquero-Garcia, J.; Barrera, A.; Gazzara, M.R.; Gonzalez-Vallinas, J.; Lahens, N.F.; Hogenesch, J.B.; Lynch, K.W.; Barash, Y. A new view of transcriptome complexity and regulation through the lens of local splicing variations. Elife 2016, 5, e11752. [Google Scholar] [CrossRef]
  49. Nowakowski, T.J.; Rani, N.; Golkaram, M.; Zhou, H.J.R.; Alvarado, B.; Huch, K.; West, J.A.; Leyrat, A.; Pollen, A.A.; Kriegstein, A.R.; et al. Regulation of cell-type-specific transcriptomes by microRNA networks during human brain development. Nat. Neurosci. 2018, 21, 1784–1792. [Google Scholar] [CrossRef]
  50. Moroney, J.B.; Vasudev, A.; Pertsemlidis, A.; Zan, H.; Casali, P. Integrative transcriptome and chromatin landscape analysis reveals distinct epigenetic regulations in human memory B cells. Nat. Commun. 2020, 11, 5435. [Google Scholar] [CrossRef]
  51. Duan, H.M.; Ge, W.H.; Zhang, A.F.; Xi, Y.; Chen, Z.H.; Luo, D.D.; Cheng, Y.; Fan, K.S.; Horvath, S.; Sofroniew, M.V.; et al. Transcriptome analyses reveal molecular mechanisms underlying functional recovery after spinal cord injury. Proc. Natl. Acad. Sci. USA 2015, 112, 13360–13365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Yang, X.P.; Coulombe-Huntington, J.; Kang, S.L.; Sheynkman, G.M.; Hao, T.; Richardson, A.; Sun, S.; Yang, F.; Shen, Y.A.; Murray, R.R.; et al. Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing. Cell 2016, 164, 805–817. [Google Scholar] [CrossRef] [Green Version]
  53. Zhang, J.; Manley, J.L. Misregulation of Pre-mRNA Alternative Splicing in Cancer. Cancer Discov. 2013, 3, 1228–1237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Prinos, P.; Garneau, D.; Lucier, J.F.; Gendron, D.; Couture, S.; Boivin, M.; Brosseau, J.P.; Lapointe, E.; Thibault, P.; Durand, M.; et al. Alternative splicing of SYK regulates mitosis and cell survival. Nat. Struct. Mol. Biol. 2011, 18, 673–680. [Google Scholar] [CrossRef] [PubMed]
  55. Lin, J.C. Therapeutic Applications of Targeted Alternative Splicing to Cancer Treatment. Int. J. Mol. Sci. 2018, 19, 75. [Google Scholar] [CrossRef] [Green Version]
  56. Miura, K.; Fujibuchi, W.; Unno, M. Splice isoforms as therapeutic targets for colorectal cancer. Carcinogenesis 2012, 33, 2311–2319. [Google Scholar] [CrossRef] [Green Version]
  57. Gao, D.D.; Morini, E.; Salani, M.; Krauson, A.J.; Chekuri, A.; Sharma, N.; Ragavendran, A.; Erdin, S.; Logan, E.M.; Li, W.C.; et al. A deep learning approach to identify gene targets of a therapeutic for human splicing disorders. Nat. Commun. 2021, 12, 3332. [Google Scholar] [CrossRef]
  58. Koedoot, E.; Wolters, L.; van de Water, B.; Devedec, S.E.L. Splicing regulatory factors in breast cancer hallmarks and disease progression. Oncotarget 2019, 10, 6021–6037. [Google Scholar] [CrossRef] [Green Version]
  59. Damen, M.P.F.; van Rheenen, J.; Scheele, C.L.G.J. Targeting dormant tumor cells to prevent cancer recurrence. FEBS J. 2021, 288, 6286–6303. [Google Scholar] [CrossRef]
  60. Varlakhanova, N.; Hahm, J.B.; Privalsky, M.L. Regulation of SMRT corepressor dimerization and composition by MAP kinase phosphorylation. Mol. Cell. Endocrinol. 2011, 332, 180–188. [Google Scholar] [CrossRef] [Green Version]
  61. Leung, M.H.; Tsoi, H.; Gong, C.; Man, E.P.S.; Zona, S.; Yao, S.; Lam, E.W.F.; Khoo, U.S. A Splice Variant of NCOR2, BQ323636.1, Confers Chemoresistance in Breast Cancer by Altering the Activity of NRF2. Cancers 2020, 12, 533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Tsoi, H.; Tsang, W.C.; Man, E.P.S.; Leung, M.H.; You, C.P.; Chan, S.Y.; Chan, W.L.; Khoo, U.S. Checkpoint Kinase 2 Inhibition Can Reverse Tamoxifen Resistance in ER-Positive Breast Cancer. Int. J. Mol. Sci. 2022, 23, 12290. [Google Scholar] [CrossRef] [PubMed]
  63. Chen, Y.H.; Huang, Q.Y.; Liu, W.; Zhu, Q.; Cui, C.P.; Xu, L.; Guo, X.; Wang, P.; Liu, J.W.; Dong, G.L.; et al. Mutually exclusive acetylation and ubiquitylation of the splicing factor SRSF5 control tumor growth. Nat. Commun. 2018, 9, 2464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Expression of BQ modulated by SRSF5 in breast cancer. (A) SRSF5 and SRSF3 down-regulated in both LCC2 and AK47 cell lines. LCC2 and AK-47 are high BQ-expressing cell lines with tamoxifen (TAM) resistance, while MCF-7 and ZR-75 are TAM-sensitive cell lines with low BQ expression. Abbreviations for TAM-resistant and TAM-sensitive are TamR and TamS, respectively. TamR/TamS indicates the comparison between TamR cell lines (LCC2 and AK-47) and TamS cell lines (MCF-7 and ZR-75). qPCR was used to determine the expression of SRSF1-12. Actin was employed as the internal control. ∆∆CT method was used to evaluate the relative gene expression level. The expression of each SRSF was compared between TAM-resistant and TAM-sensitive cell lines. LCC2 and AK-47 were derived from MCF-7 and ZR-75, respectively. The comparison was made between LCC2 vs. MCF-7 and ZK-47 vs. ZR-75. Results are shown as mean ± SD from three independent experiments. (B) Protein expressions of BQ, SRSF3 and SRSF5 determined by Western blot in different cell lines. HSP90 was the loading control. (C) Knockdown of SRSF5 enhanced protein expression of BQ in MCF-7 and ZR-75. GAPDH was the loading control. (D) Overexpression of SRSF5 reduced BQ expression. LCC2 and AK-47 were stably transfected with pCMV-Myc-SRSF5 (SRSF5 OE) or pCMV (Ctrl OE). The expression of the protein candidates was determined by Western blot with HSP90 or GAPDH as the loading control. Results were obtained from four independent experiments.
Figure 1. Expression of BQ modulated by SRSF5 in breast cancer. (A) SRSF5 and SRSF3 down-regulated in both LCC2 and AK47 cell lines. LCC2 and AK-47 are high BQ-expressing cell lines with tamoxifen (TAM) resistance, while MCF-7 and ZR-75 are TAM-sensitive cell lines with low BQ expression. Abbreviations for TAM-resistant and TAM-sensitive are TamR and TamS, respectively. TamR/TamS indicates the comparison between TamR cell lines (LCC2 and AK-47) and TamS cell lines (MCF-7 and ZR-75). qPCR was used to determine the expression of SRSF1-12. Actin was employed as the internal control. ∆∆CT method was used to evaluate the relative gene expression level. The expression of each SRSF was compared between TAM-resistant and TAM-sensitive cell lines. LCC2 and AK-47 were derived from MCF-7 and ZR-75, respectively. The comparison was made between LCC2 vs. MCF-7 and ZK-47 vs. ZR-75. Results are shown as mean ± SD from three independent experiments. (B) Protein expressions of BQ, SRSF3 and SRSF5 determined by Western blot in different cell lines. HSP90 was the loading control. (C) Knockdown of SRSF5 enhanced protein expression of BQ in MCF-7 and ZR-75. GAPDH was the loading control. (D) Overexpression of SRSF5 reduced BQ expression. LCC2 and AK-47 were stably transfected with pCMV-Myc-SRSF5 (SRSF5 OE) or pCMV (Ctrl OE). The expression of the protein candidates was determined by Western blot with HSP90 or GAPDH as the loading control. Results were obtained from four independent experiments.
Cancers 15 02271 g001
Figure 2. SRSF5 bound to exon 11 of NCOR2 to facilitate the splicing process to generate NCOR2 mRNA. (A) Schematic diagram of splicing event that produces NCOR2 and BQ. The arrows indicate the positions of the specific primers used for detecting NCOR2-BQ splicing (NCOR2-BQ splicing -F and -R), NCOR2 and BQ. The diagram is not in scale. (B) SRSF5 could bind to exon 11 of NCOR2. MCF-7 and ZR-75 cells were transfected with pCMV-myc-SRSF5 mammalian expression plasmid. The cells were harvested 48 h post-transfection. Anti-myc antibody was employed to immunoprecipitate the SRSF5/RNA complex. The presence of exon 11 was determined by qPCR. The expression level was relative to the input. Results were shown as mean ± SD from three independent experiments. Students’ t-test was employed to determine the statistical significance. (C) Knockdown of SRSF5 could enhance the expression of BQ while reducing that of NCOR2 in MCF7 and ZR-75 cells. qPCR was employed to determine the mRNA expression of BQ and NCOR2. Actin was the internal control. The bar charts showed results as mean ± SD from three independent experiments. A one-way ANOVA with a Bonferroni post hoc test determined the statistical significance between shCtrl and shSRSF5.1 or shSRSF5.2 groups. (D) SRSF5 overexpression could reduce the expression of BQ while enhancing that of NCOR2 in LCC2 and AK-47 cells. BQ and NCOR2 expression levels were determined by qPCR. The results showed the mean expression of BQ and NCOR2 ± SD from three independent experiments. The statistical significance between groups was determined by students’ t-test. *, **, and *** indicate p < 0.05, p < 0.01 and p < 0.001, respectively.
Figure 2. SRSF5 bound to exon 11 of NCOR2 to facilitate the splicing process to generate NCOR2 mRNA. (A) Schematic diagram of splicing event that produces NCOR2 and BQ. The arrows indicate the positions of the specific primers used for detecting NCOR2-BQ splicing (NCOR2-BQ splicing -F and -R), NCOR2 and BQ. The diagram is not in scale. (B) SRSF5 could bind to exon 11 of NCOR2. MCF-7 and ZR-75 cells were transfected with pCMV-myc-SRSF5 mammalian expression plasmid. The cells were harvested 48 h post-transfection. Anti-myc antibody was employed to immunoprecipitate the SRSF5/RNA complex. The presence of exon 11 was determined by qPCR. The expression level was relative to the input. Results were shown as mean ± SD from three independent experiments. Students’ t-test was employed to determine the statistical significance. (C) Knockdown of SRSF5 could enhance the expression of BQ while reducing that of NCOR2 in MCF7 and ZR-75 cells. qPCR was employed to determine the mRNA expression of BQ and NCOR2. Actin was the internal control. The bar charts showed results as mean ± SD from three independent experiments. A one-way ANOVA with a Bonferroni post hoc test determined the statistical significance between shCtrl and shSRSF5.1 or shSRSF5.2 groups. (D) SRSF5 overexpression could reduce the expression of BQ while enhancing that of NCOR2 in LCC2 and AK-47 cells. BQ and NCOR2 expression levels were determined by qPCR. The results showed the mean expression of BQ and NCOR2 ± SD from three independent experiments. The statistical significance between groups was determined by students’ t-test. *, **, and *** indicate p < 0.05, p < 0.01 and p < 0.001, respectively.
Cancers 15 02271 g002
Figure 3. Knockdown of SRSR5 could induce tamoxifen (TAM) resistance. (A) Knockdown of SRSF5 could make MCF-7 and ZR-75 become TAM-resistant. The cells were treated with 5 μM of TAM, and a cell viability assay was conducted after 96 h of the treatment. (B) Results of the clonogenic assay. The cells were treated with 5 μM of TAM for 21 days. Results are shown as mean ± SD from five independent experiments. A one-way ANOVA with a Bonferroni post hoc test was used to determine the statistical significance between shCtrl and shSRSF5.1 or shSRSF5.2 groups. Clonogenicity (%) was determined in the MCF-7 panel, while cell viability (%) was determined in the ZR-75 panel due to the lack of a countable colony. (C) Knockdown of SRSF5 could confer TAM resistance in vivo. ZR-75 shCtrl and ZR-75 shSRSF5.1 cells were engrafted onto the mammary fatpad of nude female mice. TAM or saline was given to the mice via subcutaneous injection. A one-way ANOVA with a Bonferroni post hoc test was used to determine the statistical significance between ZR-75 shCtrl and ZR-75 shSRSF5.1 groups treated with TAM. (D) Knockdown of SRSF5 could enhance BQ expression in the tumour tissues. Proteins were extracted from the tumour tissues, and Western blot was employed to determine BQ expression. Protein lysates from 3 independent tumours from each of the groups were analysed. GAPDH was the loading control. *** represents p < 0.001.
Figure 3. Knockdown of SRSR5 could induce tamoxifen (TAM) resistance. (A) Knockdown of SRSF5 could make MCF-7 and ZR-75 become TAM-resistant. The cells were treated with 5 μM of TAM, and a cell viability assay was conducted after 96 h of the treatment. (B) Results of the clonogenic assay. The cells were treated with 5 μM of TAM for 21 days. Results are shown as mean ± SD from five independent experiments. A one-way ANOVA with a Bonferroni post hoc test was used to determine the statistical significance between shCtrl and shSRSF5.1 or shSRSF5.2 groups. Clonogenicity (%) was determined in the MCF-7 panel, while cell viability (%) was determined in the ZR-75 panel due to the lack of a countable colony. (C) Knockdown of SRSF5 could confer TAM resistance in vivo. ZR-75 shCtrl and ZR-75 shSRSF5.1 cells were engrafted onto the mammary fatpad of nude female mice. TAM or saline was given to the mice via subcutaneous injection. A one-way ANOVA with a Bonferroni post hoc test was used to determine the statistical significance between ZR-75 shCtrl and ZR-75 shSRSF5.1 groups treated with TAM. (D) Knockdown of SRSF5 could enhance BQ expression in the tumour tissues. Proteins were extracted from the tumour tissues, and Western blot was employed to determine BQ expression. Protein lysates from 3 independent tumours from each of the groups were analysed. GAPDH was the loading control. *** represents p < 0.001.
Cancers 15 02271 g003
Figure 4. Overexpression of SRSF5 could compromise tamoxifen (TAM) resistance. (A) Overexpression of SRSF5 could make LCC2 sensitive to TAM. An amount of 5 μM of TAM was employed to treat the cells. After 96 h of treatment, cell viability was assayed. Results are shown in the bar chart as mean value of cell viability ± SD from five independent experiments. Statistical significance compared with the untreated group was determined using a one-way ANOVA with a Bonferroni post hoc test. (B) Clonogenic assay indicated the effect of SRSF5 overexpression on TAM resistance in LCC2. An amount of 5 μM of TAM was employed to treat LCC2 cells for 3 weeks. The bar chart shows the results of mean value of cell viability ± SD from five independent experiments. Statistical significance between Ctrl OE and SRSF5 OE was determined by students’ t-test. (C) Overexpression of SRSF5 could recover the sensitivity to TAM in vivo. LCC2 Ctrl OE and LCC2 SRSF5 OE cells were engrafted onto the mammary fatpad of female nude mice. TAM or saline was given to the mice via subcutaneous injection. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance between indicated groups. (D) Overexpression of SRSF5 could reduce BQ expression in the tumour tissues. Proteins were extracted from the tumour tissues, and Western blot was employed to determine BQ expression. Protein lysates from 3 independent tumours from each of the groups were analysed. GAPDH was the loading control. *** represents p < 0.001.
Figure 4. Overexpression of SRSF5 could compromise tamoxifen (TAM) resistance. (A) Overexpression of SRSF5 could make LCC2 sensitive to TAM. An amount of 5 μM of TAM was employed to treat the cells. After 96 h of treatment, cell viability was assayed. Results are shown in the bar chart as mean value of cell viability ± SD from five independent experiments. Statistical significance compared with the untreated group was determined using a one-way ANOVA with a Bonferroni post hoc test. (B) Clonogenic assay indicated the effect of SRSF5 overexpression on TAM resistance in LCC2. An amount of 5 μM of TAM was employed to treat LCC2 cells for 3 weeks. The bar chart shows the results of mean value of cell viability ± SD from five independent experiments. Statistical significance between Ctrl OE and SRSF5 OE was determined by students’ t-test. (C) Overexpression of SRSF5 could recover the sensitivity to TAM in vivo. LCC2 Ctrl OE and LCC2 SRSF5 OE cells were engrafted onto the mammary fatpad of female nude mice. TAM or saline was given to the mice via subcutaneous injection. One-way ANOVA with Bonferroni post hoc test was used to determine the statistical significance between indicated groups. (D) Overexpression of SRSF5 could reduce BQ expression in the tumour tissues. Proteins were extracted from the tumour tissues, and Western blot was employed to determine BQ expression. Protein lysates from 3 independent tumours from each of the groups were analysed. GAPDH was the loading control. *** represents p < 0.001.
Cancers 15 02271 g004
Figure 5. The clinical significance of SRSF5 in ER +ve breast cancer. (A) Expressions of SRSF5 and BQ in tumour tissues examined by immunohistochemistry. Representative images are shown. The scale bar represents 20 μm. (B) The expressions of SRSF5 and BQ showed a negative correlation. Linear regression was performed to determine the relationship, and the Spearman coefficient was determined. Low expression of SRSF5 was associated with (C) tamoxifen (TAM) resistance, (D) local recurrence and (E) metastasis. Chi-square test was employed. Kaplan–Meier survival analysis showed that the cases of ER +ve breast cancer with low expression of SRSF5 had poorer outcomes of (F) overall and (G) disease-specific survival.
Figure 5. The clinical significance of SRSF5 in ER +ve breast cancer. (A) Expressions of SRSF5 and BQ in tumour tissues examined by immunohistochemistry. Representative images are shown. The scale bar represents 20 μm. (B) The expressions of SRSF5 and BQ showed a negative correlation. Linear regression was performed to determine the relationship, and the Spearman coefficient was determined. Low expression of SRSF5 was associated with (C) tamoxifen (TAM) resistance, (D) local recurrence and (E) metastasis. Chi-square test was employed. Kaplan–Meier survival analysis showed that the cases of ER +ve breast cancer with low expression of SRSF5 had poorer outcomes of (F) overall and (G) disease-specific survival.
Cancers 15 02271 g005
Figure 6. SRPKIN-1 could reverse TAM resistance by modulating the activity of SRSF5 towards exon 11 of NCOR2. (A) SRPK1 could interact with SRSF5 in LCC2. SRPK1 was immunoprecipitated. Anti-SRPK1 and anti-SRSF5 were employed to analyse the immunoprecipitant by Western blot. Results were obtained from three independent experiments. SRPKIN-1 could reduce BQ expression but enhance NCOR2 expression at (B) mRNA and (C) protein levels. qPCR was employed to examine the mRNAs. Actin was the internal control. ∆∆CT method was used to determine the relative expression of the candidates. Results are shown as mean ± SD from three independent experiments. (D) SRPKIN-1 could reduce the degree of phosphorylation on SRSF5. LCC2 was transfected with pCMV_Myc_SRSF5. Immunoprecipitation with anti-Myc was performed 48 h post-transfection. pSR represents phosphorylated SRSF proteins. Results were obtained from three independent experiments. *** represents p < 0.001. (E) SRPKIN-1 could enhance the proportion of SRSF5 interacting with exon 11 of NCOR2. Anti-myc antibody was employed to immunoprecipitate the SRSF5/RNA complex. The presence of exon 11 was determined by qPCR. The expression level was relative to the input. Results are shown as mean ± SD from three independent experiments. Students’ t-test was employed to determine the statistical significance between DMSO and SRPKIN-1 groups. (F) SRPKIN-1 could reverse TAM resistance. LCC2 was treated with 5 μM of TAM and 100 nM of SRPKIN-1 for 96 h. Cell viability was determined. Results are shown as mean ± SD from four independent experiments. Students’ t-test was employed to determine the statistical significance between DMSO and SRPKIN-1 groups. *** represents p < 0.001.
Figure 6. SRPKIN-1 could reverse TAM resistance by modulating the activity of SRSF5 towards exon 11 of NCOR2. (A) SRPK1 could interact with SRSF5 in LCC2. SRPK1 was immunoprecipitated. Anti-SRPK1 and anti-SRSF5 were employed to analyse the immunoprecipitant by Western blot. Results were obtained from three independent experiments. SRPKIN-1 could reduce BQ expression but enhance NCOR2 expression at (B) mRNA and (C) protein levels. qPCR was employed to examine the mRNAs. Actin was the internal control. ∆∆CT method was used to determine the relative expression of the candidates. Results are shown as mean ± SD from three independent experiments. (D) SRPKIN-1 could reduce the degree of phosphorylation on SRSF5. LCC2 was transfected with pCMV_Myc_SRSF5. Immunoprecipitation with anti-Myc was performed 48 h post-transfection. pSR represents phosphorylated SRSF proteins. Results were obtained from three independent experiments. *** represents p < 0.001. (E) SRPKIN-1 could enhance the proportion of SRSF5 interacting with exon 11 of NCOR2. Anti-myc antibody was employed to immunoprecipitate the SRSF5/RNA complex. The presence of exon 11 was determined by qPCR. The expression level was relative to the input. Results are shown as mean ± SD from three independent experiments. Students’ t-test was employed to determine the statistical significance between DMSO and SRPKIN-1 groups. (F) SRPKIN-1 could reverse TAM resistance. LCC2 was treated with 5 μM of TAM and 100 nM of SRPKIN-1 for 96 h. Cell viability was determined. Results are shown as mean ± SD from four independent experiments. Students’ t-test was employed to determine the statistical significance between DMSO and SRPKIN-1 groups. *** represents p < 0.001.
Cancers 15 02271 g006
Table 1. DNA sequences of primers.
Table 1. DNA sequences of primers.
Forward (5′ → 3′)Reverse (5′ → 3′)
SRSF1GCG GTC TGA AAA CAG AGT GGACA AAC TCC ACG ACA CCA GT
SRSF2TCC AAG TCC AAG TCC TCG TCGCG ACC TGG ATT TGG ATT CC
SRSF3CGA AGT GTG TGG GTT GCT AGAGT TCC ACT CTT ACA CGG CA
SRSF4GGA AGG TCG AGG AGA GAG TGGCT GAC TTT GAT CTG GAG CG
SRSF5GGT GGA AGA GGT AGA GGA CGTTT GAG ATC CTG CCA GCT GA
SRSF6CGT GCT TTG GAC AAA CTG GATCG AGA CCT GGA TCT GCT TC
SRSF7CGA AGA AGA AGC AGG TCA CGGGA GAT GCT GAC CTT GAC CT
SRSF8GTC TCA CTC GAA GTC TGG GTTAG ATG AAG ACC TGG ACC GC
SRSF9GAG AAG CTG GGG ATG TCT GTTCC AGT TTA CGC AGG GCA TA
SRSF10ATT TCT ACA CTC GCC GTC CAGCC GTC CAC AAA TCC ACT TT
SRSF11GGG GCT CCT ACT CTT GAT CCTGA AAC GAG ACC AGC AGC TA
SRSF12AAT AGG AGG CGG TCA GAC AGGTT CTT GAC TGC CTT GCT GA
ActinATC GTG CGT GAC ATT AAG GAG AAGAGG AAG GAA GGC TGG AAG AGT G
BQGGA GCG CAT GCA GAG AAC CCTG GCG GTC TTT GTA CAC CT
NCOR2GGT GGA GCG CAT CGA GAA CCCC GCT GGC CCA CCC TCT GCA TG
NCOR2-BQ splicingGAG CAGA AGT TCT GCC AGC GCTTC TCC CGG AAG GTC TCC TTC
Table 2. Clinical information about breast cancer patients.
Table 2. Clinical information about breast cancer patients.
Number of CasesPercentage (%)
ER +ve breast cancer 137100
Age<566547.4
≥5652.652.6
T stageI, II4129.9
III, IV107.3
Lymph node statusPositive6748.9
Negative5640.9
Tumour grade1, 24935.8
37454.0
Tumour size<2 cm1410.2
≥2 cm7957.7
ER statusPositive7353.3
Negative2518.2
PR statusPositive4935.8
Negative3727.0
HER2 receptor statusPositive3525.5
Negative3727.0
Triple negative statusYes139.5
No7151.8
Table 3. Univariate cox regression analysis of overall and disease-specific survival in ER +ve breast cancer.
Table 3. Univariate cox regression analysis of overall and disease-specific survival in ER +ve breast cancer.
Overall SurvivalDisease-Specific Survival
Clinical–Pathological ParametersNo. of CasesRR (95% CI)p-ValueRR (95% CI)p-Value
Age711.736 (0.811, 3.718)0.1551.421 (0.564, 3.584)0.456
T-stage325.897 (1.311, 26.531)0.0213.931 (0.714, 21.647)0.116
Lymph node involvement650.948 (0.423, 2.214)0.8961.323 (0.478, 3.662)0.589
Tumour grade701.336 (0.613, 2.913)0.4663.319 (1.092, 10.087)0.034
Histological type711.342 (0.317, 5.689)0.6900.791 (0.182, 3.447)0.755
PR status590.237 (0.096, 0.582)0.0020.111 (0.038, 0.330)<0.001
HER2 status491.110 (0.427, 2.888)0.8301.715 (0.499, 5.893)0.392
Tumour size541.691 (0.390, 7.327)0.48226.528 (0.035, 20173.302)0.333
Cases with high SRSF5580.182 (0.071, 0.466)<0.0010.142 (0.039, 0.519)0.003
Cases with high BQ647.450 (2.559, 21.692)<0.00121.188 (2.802, 160.238)0.003
Either high SRSF5/low BQ or low SRSF5/high BQ400.122 (0.035, 0.421)0.0010.060 (0.008, 0.466)0.007
Table 4. Multivariate cox regression analysis of overall and disease-specific survival in ER +ve breast cancer.
Table 4. Multivariate cox regression analysis of overall and disease-specific survival in ER +ve breast cancer.
Overall SurvivalDisease-Specific Survival
Clinical–Pathological ParametersCasesRR (95% CI)p-ValueRR (95% CI)p-Value
T-stage2112.931 (0.708, 236.235)0.0842.583 (0.537, 12.436)0.237
PR status210.361 (0.053, 2.437)0.2960.288 (0.084, 0.986)0.047
Cases with high SRSF5210.087 (0.008, 0.901)0.0410.178 (0.035, 0.920)0.039
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tsoi, H.; Fung, N.N.-C.; Man, E.P.S.; Leung, M.-H.; You, C.-P.; Chan, W.-L.; Chan, S.-Y.; Khoo, U.-S. SRSF5 Regulates the Expression of BQ323636.1 to Modulate Tamoxifen Resistance in ER-Positive Breast Cancer. Cancers 2023, 15, 2271. https://doi.org/10.3390/cancers15082271

AMA Style

Tsoi H, Fung NN-C, Man EPS, Leung M-H, You C-P, Chan W-L, Chan S-Y, Khoo U-S. SRSF5 Regulates the Expression of BQ323636.1 to Modulate Tamoxifen Resistance in ER-Positive Breast Cancer. Cancers. 2023; 15(8):2271. https://doi.org/10.3390/cancers15082271

Chicago/Turabian Style

Tsoi, Ho, Nicholas Nok-Ching Fung, Ellen P. S. Man, Man-Hong Leung, Chan-Ping You, Wing-Lok Chan, Sum-Yin Chan, and Ui-Soon Khoo. 2023. "SRSF5 Regulates the Expression of BQ323636.1 to Modulate Tamoxifen Resistance in ER-Positive Breast Cancer" Cancers 15, no. 8: 2271. https://doi.org/10.3390/cancers15082271

APA Style

Tsoi, H., Fung, N. N. -C., Man, E. P. S., Leung, M. -H., You, C. -P., Chan, W. -L., Chan, S. -Y., & Khoo, U. -S. (2023). SRSF5 Regulates the Expression of BQ323636.1 to Modulate Tamoxifen Resistance in ER-Positive Breast Cancer. Cancers, 15(8), 2271. https://doi.org/10.3390/cancers15082271

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop