Abstract
The initiation and progression of breast cancer (BRCA) is associated with inflammation and immune-overactivation. The E3 ubiquitin ligase is known to subtly balance immune-overactivation and pro-tumorigenesis. Here, by global transcriptional profiling of BRCA patient tissues, we identified a signature expression profile of F-box factors, of which FBXL8 emerged as a novel key component of E3 ligase. Our ex vivo studies (n=1349) indicate that FBXL8 promotes tumorigenesis and its level escalates with BRCA advancement. Knockdown of FBXL8 caused: (i) accumulation of two tumor-suppressors,CCND2 and IRF5, suggesting their collaborative regulation of BRCA status, (ii) intrinsic apoptosis in BRCA, (iii) inhibition of cell migration and invasion and (iv) downregulation of cancer-promoting cytokines/ chemokines; all of these effects curtailed the tumor microenvironment and suppressed cancer progression. Our findings highlight the translational impact of exploiting FBXL8 and its interaction partners (CCND2 and IRF5) for developing anti-cancer strategies and potential therapeutics to limit BRCA progression.
Introduction
Breast cancer (BRCA) is the most common malignancy in women, being responsible for more than 500,000 deaths each year worldwide. In Western countries, BRCA constitutes 25% of female malignancies (1, 2). Increasing evidence suggest that the ubiquitin machinery plays a pivotal role in cancer development and progression (3, 4). However, the key factors and underlying mechanisms involved in BRCA are not fully explored. MLN4924 is a newly discovered investigational small molecule inhibitor of NAE (NEDD8 Activating Enzyme), currently under Phase I clinical trial. MLN4924 was reported to inactivate SCF E3 ligase by inhibiting cullin neddylation (5). Cumulative studies have demonstrated that MLN4924 inhibited cancer cell growth, including BRCA (6, 7). Both in vivo and in vitro studies showed that upon stimulation, the SCF E3 ubiquitin machinery can induce dynamic changes to a series of cytokine expression, which favors a pro-tumorigenic microenvironment (8, 9). Therefore, there may be vital members of the SCF E3 ubiquitin ligase family, which could drive pro-tumorigenesis. Being the largest family of E3 ubiquitin ligases, the SCF is composed of four components: SKP1, a cullin, an F-box protein, and RBX1 or RBX2. This multimeric complex is responsible for ubiquitination-mediated degradation of target proteins involved in various biological processes (10). A well-regulated SCF E3 ubiquitin ligase maintains homeostasis in cell proliferation and genome stability, hence preventing carcinogenesis (10).
Previously, we demonstrated that SCF (Skp1-cullin-F-box proteins) E3 ubiquitin ligase critically balances immune-overactivation and pro-tumorigenesis (11, 12), prompting us to hypothesize a vicious network of chronic inflammation-driven activation of SCF E3 ligase, which probably facilitates ubiquitination-mediated degradation of cancer-modulating factors and also perturbs the pro-tumorigenic cytokines. To test this hypothesis, here we aimed to first elucidate the specific SCF E3 ligases and their potential pathobiological partners, and then explore their mechanisms in driving BRCA progression, with a view to understanding how the SCF E3 ubiquitin ligases might act as significant risk factors in cancer advancement (13). Since the F-box is responsible for the functional specificity of SCF E3 tetrameric complex (14), we focused our attention on the expression profiles of F-box factors. Notably, there are 71 genes in the human “F-box gene family”, including: (1) 21 F-boxes grouped as leucine-rich repeat protein, (2) 11 F-boxes with WD repeat domain and (3) 39 F-boxes grouped as “others” (15-17).
By ex vivo characterization of the “signature profile” of F-box factors in patients’ BRCA tissues, we identified for the first time, a novel F-Box factor, FBXL8, which was significantly upregulated with disease advancement. FBXL8 was found to interact with two tumor-suppressors, CCND2 and IRF5. Knockdown of FBXL8 in BRCA cells resulted in the accumulation of CCND2 and IRF5, concomitant with cancer suppression. CCND2, a member of the D-type cyclins has been implicated in cell cycle regulation, differentiation and malignant transformation, and is seemingly inactivated in human cancers (18). IRF5 is a transcription factor, purported to be a cancer-suppressor which regulates apoptosis and immune activation (19, 20). Our findings will open new avenues for further translational research on E3 ubiquitin ligase-related regulation and intervention of cancer. We propose that manipulation of FBXL8 and its binding partners could beneficially dampen the pro-tumorigenic microenvironment, and curb metastasis.
Results
Clinicopathological parameters of BRCA patient tissues used for the global RNA sequencing
Clinical information of five representative BRCA patients with confirmed cancer staging, was obtained retrospectively from the medical records of BioBank, Taipei Medical University. The clinical stages of patients were selected based on tumor-node-metastasis (TNM): stage I (n = 1), stage II (n = 2) and stage III (n = 2). The age of the BRCA patients ranged from 38 to 69 years, with a mean of 59 years. The histology type was diagnosed as invasive lobular carcinoma (ILC) or invasive ductal carcinoma (IDC). BRCA resection pathology reports included primary tumor, categories and pathologic lymph node status. Figure 1A summarizes the clinicopathological parameters of the BRCA patients.
FBXL8 is significantly upregulated in human primary BRCA tissues
To explore the clinical significance of the key components of SCF multimeric complex which might respond to BRCA progression, and to identify the corresponding F-box factors involved, we performed an NGS-based RNA-Seq to profile the gene expression of the F-Box members in human primary breast carcinoma tissues. The five breast carcinoma tissues which have been clinically staged and characterized alongside corresponding normal counterpart tissues, were subjected to global mRNA profiling. We found that several SCF components were overexpressed in BRCA, including SKP1 (S-phase kinase-associated protein 1), SKP2, CUL1 (cullin 1), FBXW7 (F-box/WD repeat-containing protein 7), SAG (Sensitive to apoptosis gene), FBXO4 (F-box only protein 4) and RBX1 (ring-box 1) (Supplementary Figure S1). Red arrow (↑) highlights the upregulation of the identified F-box gene, FBXL8 (Figure 1A). Compared to normal control tissues, the mRNA level of FBXL8 in carcinoma tissues was significantly upregulated by up to 7.5-fold (*, False Discovery Rate, FDR<0.05). Scatter plot (Figure 1B) illustrates quantitative differences in the expression levels of FBXL8 mRNA in individual clinical samples, with 3-fold increase in breast carcinoma (expressed in log2). We found that amongst all 71 F-box family members, FBXL8, a novel F-box factor, was prominently upregulated in breast carcinoma tissues, compared to the corresponding paired normal tissues (FDR<0.05). Consistently, RNA-Sequencing data in TCGA (The Cancer Genome Atlas) showed that the mRNA expression of FBXL8 was significantly upregulated in BRCA tissues (Figure 1C, ** p<0.01). Since information on FBXL8 is lacking in literature, we performed sequence alignment of FBXL8 with several mammalian species using Clustal Omega and found that FBXL8 contains highly conserved amino acid sequence amongst mammals (Figure 1D), suggesting its fundamental biological roles.
FBXL8 is upregulated with the advancement of BRCA stages
To affirm the RNA-Seq data (Fig. 1) and to further understand the correlation between FBXL8 expression and the clinical status of BRCA patients, we performed qRT-PCR and IHC staining to retrospectively analyze more clinical BRCA tissues (Figure 2). The FBXL8 mRNA profiles of 32 paired samples (n=64) were determined. Supplementary Figure S2 shows the clinicopathological parameters of the BRCA patient breast tissue samples used for qRT-PCR analysis. Consistently, the mRNA levels of FBXL8 were significantly elevated (p<0.01) in primary BRCA tissues (Figure 2A). The expression levels of FBXL8 correlated significantly (p<0.01) with clinical staging. In addition, SAG (sensitive to apoptosis gene), another important component of SCF E3 ubiquitin ligase, which confers anti-apoptosis and promotes liver cancer (12, 21, 22), was also significantly upregulated (p<0.01) in all stages of BRCA (Figure 2B). This further corroborates the critical role of the members of SCF E3 ubiquitin ligase in BRCA progression. IHC of FBXL8 in 30 paired BRCA tissues (n=60) also showed increase in FBXL8 with disease progression from stages IA to IIA/B to IIIA/B/C (Figure 2C). The clinicopathological parameters of the corresponding tissue sections are shown in Supplementary Figures S3. The expression of both the mRNA and protein of FBXL8 in BRCA tissues were consistently increased by up to 4-fold compared to normal tissues, and the levels in BRCA tissues were found to elevate with the advancement of cancer (p<0.001) (Figure 2D). Clearly, the FBXL8, a newly identified F-box factor is associated with the progression of BRCA. Furthermore, a rise in the ubiquitin protein in BRCA tissues (Figure 2E) also corroborates the significance of the ubiquitin proteasome system (UPS) in BRCA development. Altogether, we have demonstrated simultaneous increases in the levels of FBXL8, SAG and ubiquitin, all of which are functionally related to the UPS protein degradation system, in BRCA formation.
FBXL8 is a novel anti-apoptosis factor in BRCA
To elucidate the mechanisms underlying the pathophysiological significance of FBXL8, we performed gene specific-siRNA knockdown of FBXL8 in breast carcinoma cells to assess the potential loss/gain-of-function. We first examined the endogenous levels of FBXL8 in two human BRCA cell lines, including a highly invasive MDA-MB231, a less invasive MCF7, and a control cell line, MCF10A which expresses very low basal level of FBXL8. Prior to knockdown, we observed 23- and 15-folds higher expression of FBXL8 in MCF7 and MDA-MB231, respectively, compared to that of the control MCF10A cells (p<0.01) (Figure 3A). With up to 95% knockdown efficiency of FBXL8 in MCF7 and MDA-MB231 cells (Figure 3B), we noted a reduction in cell viability by 52.5% within 48 h (an average taken from MTT and CTB assays) (Figure 3C, D), which was significantly different (p<0.001) from that of control siRNA-treated cells.
FBXL8-specific siRNA knockdown effected 2-3-fold reduction in cell proliferation of both MCF7 and MDA-MB231 cells (Figure 3E, F, p<0.001). This reduction in cell proliferation is supported by significant increases in early apoptosis of 20% and 24% in FBXL8-knocked down MCF7 and MDA-MB231 cells, respectively (Figure 3G). Representative histograms of apoptosis assay are shown in Supplementary Figure S4 (p<0.001). Further investigation of Caspase activity showed that both Caspases-9 and −3 were activated in FBXL8 knocked down cells (Figure 3H), indicating that the presence/ increase in FBXL8 prevents intrinsic apoptosis and confers survival advantage to BRCA, hence favoring cancer progression.
Knockdown of FBXL8 inhibited BRCA cell migration and invasion
To further investigate the impact of FBXL8 as a driver of BRCA metastasis, we examined the effects of FBXL8-knockdown on cell migration and invasion. When monitored over 30 h, the highly invasive MDA-MB231 cells showed 40% slowdown in cell migration compared to 14% with the less invasive MCF7 (both compared to control siRNA-treated cells, p<0.001) (Figure 4A, B, C, D). Consistently, cell invasion was suppressed by 4-8 fold in MDA-MB231 and MCF7 cells (p<0.001) (Figure 4E, F). Importantly, direct analysis of cell growth showed that the doubling time for MCF7 and MDA-MB231 cells was 26 h and 23 h, respectively (Fig. 3D). The 4-8 fold reduction in invasion observed at 24 h therefore cannot be solely attributable to the inhibition of cell growth or a difference in cell doubling time. These results again support that FBXL8 is a significant contributor to BRCA cell growth and metastasis.
Knockdown of FBXL8 suppressed pro-inflammatory cytokines/chemokines in BRCA
There is increasing evidence that the SCF E3 ubiquitin machinery favors cancer formation by controlling immunomodulators such as cytokines which are abundant in the tumor microenvironment (8, 9, 23). It was therefore pertinent for us to examine whether and how FBXL8 might regulate the cytokines. A cytokine array analysis of FBXL8-knocked down BRCA cells showed that the loss of FBXL8 caused a downregulation of the pro-tumorigenic factors, including MCP-1, I-TAC, TECK, CTACK, MIF, GM-CSF, NT-3, FGF-6, Angiogenin, ICAM-1, DtK and EGFR (Fig. 5A). Interestingly, we observed that Angiogenin, GM-CSF, ICAM-1 and EGFR were dramatically over-expressed (for example, up to 20-fold increase for EGFR) in the highly invasive MDA-MB231 cells compared to the less invasive MCF7 cells. Previous findings have also shown that these factors support the survival of breast cancer cells and promote angiogenesis, further triggering cancer invasion and metastatic spread (24-27). Additionally, NT-3, which was concomitantly decreased by knockdown of FBXL8, has been reported to modulate the tumor microenvironment to promote breast-to-brain metastasis (28). Suffice to say, knockdown of the pro-tumorigenic FBXL8 caused the loss of another pro-tumorigenic factor, NT-3, which plausibly further restrains BRCA metastasis. It would be interesting, in future, to examine the potential functional partnership between FBXL8 and NT-3, in their co-regulation of BRCA progression.
Quantification by ELISA further affirmed that depletion of FBXL8 significantly reduced the production of several inflammatory cytokines and chemokines e.g. GM-CSF, MCP-1, CTACK, EGFR, ICAM-1, I-TAC, MIF and TECK (Figure 5B, blue bars). Such cytokines/chemokines are reported to promote cancer, and are associated with cancer survival, growth and/or metastasis in BRCA (25, 29-36). The association of these cytokines to FBXL8 activity suggests that FBXL8 promotes BRCA survival by modulating the cytokine/chemokine-rich pro-tumorigenic microenvironment.
On the other hand, Osteoprotegerin (OPG, overexpressed in BRCA metastasis) is up-regulated in MCF7 and down-regulated in MDA-MB231 during FBXL8 knockdown. Although we found opposite trends of OPG in the two BRCA cell types, higher levels of OPG have been associated with an increased risk of ER-negative BRCA (37). Similarly, knockdown of FBXL8 increased the expression of pro-tumorigenic TNFR1. Inhibition of TNFR1 expression was reported to increase apoptosis, and inhibit cell proliferation and invasion in triple-negative and ER positive BRCA cells (38). IL-6R was negatively regulated by FBXL8 expression, although IL-6 appeared unaffected by FBXL8 knockdown. Such data highlight the complexity of FBXL8-associated immune network. It is plausible that the BRCA tissues contained heterogeneous tumor cell types. Therefore, our findings suggest the propensity for developing variant forms of FBXL8-based combinatorial treatments possibly with FBXL8-mediated cytokines or FBXL8-binding partners, for specific and targeted chemotherapy against different heterogeneous forms of BRCA. This prompted us to further search for and characterize FBXL8-binding partners.
Cancer-suppressors, CCND2 and IRF5, are FBXL8-binding partners in BRCA
FBXL8 is a leucine-rich repeat protein member of the E3 SCF multimeric complex. Our in silico modeling prediction revealed FBXL8-binding partners including CUL1, SKP1, CUL2 and CUL7 (Supplementary Figure S5). These results support the involvement of E3 ligase in promoting BRCA. In addition, several non-SCF E3 ligase-associated components were proposed, including cell-cycle associated proteins, cancer-promoting transcription factors and DNA repair/replication factors. Notably, two cancer suppressors, CCND2 (Cyclin D2) and IRF5 (Interferon Regulatory Factor 5), emerged prominently as potential interaction partners of FBXL8. Knockdown of CCND2 has been reported to increase the growth rate and migration ability of cancer cells (39). On the other hand, IRF5 overexpression is capable of reverting MDA-MB231 to normal acini-like structures (20). Notably, our IHC staining of breast carcinoma tissues demonstrated a reciprocal profile of FBXL8 versus CCND2 and IRF5 with cancer advancement over stages I - III. Clearly, the FBXL8 level increased while CCND2 and IRF5 levels decreased (Figure 6A, B, C), thus: (i) corroborating the significant role of CCND2 and IRF5 as BRCA-suppressors and (ii) highlighting the potential role of FBXL8 in specific degradation of CCND2 and IRF5 proteins as cancer progresses. The mRNA levels of CCND2 and IRF5 were further examined by both NGS-based RNA-Seq and TCGA database analyses. Figure 6D, E, F, G showed that the IRF5 mRNA was significantly upregulated in breast tumor tissues (Fig. 6G), although its protein level was reduced in BRCA tissues. This suggests potential post-translational modifications which led to the loss of IRF5 proteins in BRCA.
Gilt by association, we next investigated whether anti-cancer CCND2 and IRF5 are binding partners of FBXL8, with a view to providing explanations on how their interactions might impact BRCA progression. Co-immunoprecipitation (Co-IP) showed that FBXL8 pulled down both CCND2 and IRF5 in MCF-7 and MDA-MB231 cells (Figure 6H), indicating that CCND2 and IRF5 are conserved specific substrates of FBXL8 E3 ligase. IP with control IgG2b demonstrated the specificity of the assay. Importantly, we showed that knockdown of FBXL8 caused the accumulation of CCND2 and IRF5 (but not tubulin) proteins in both BRCA cell lines, suggesting that FBXL8 selectively regulates the protein degradation of CCND2 and IRF5 (Figure 6I), thus explaining why FBXL8 rises while CCND2 and IRF5 declines as BRCA progresses. Given the critical role of CCND2 and IRF5 in cell cycle regulation as growth-inhibitors (18, 20), our data here account the functional significance of FBXL8 as a tumor-promoter, whose knockdown spared the tumor suppressors, CCND2 and IRF5 from degradation, thus allowing them to accumulate and suppress BRCA progression.
Discussion
Due to the lack of early symptoms and a dearth of markers for dynamic progression of BRCA, breast cancer malignancy is often detected too late. BRCA appears to be associated with inflammation and immune-overactivation conditions, which perpetuates cancer recurrence and lowers patient survival rate (13). Concordant with our earlier findings that showed the involvement of SCF E3 in inflammation and tumorigenesis (11), treatment with non-steroidal anti-inflammatory drugs has been shown to reduce BRCA metastasis and recurrence of the disease (40, 41). Other clinical efforts to arrest BRCA have used Bortezomib (Velcade or PS-341), the first proteasome inhibitor to be approved by FDA. Unfortunately, Bortezomib was found to be cytotoxic due to its inhibition of protein degradation on a global scale (42). Therefore, we reasoned that the ideal approach is to target a specific protein degradation enzyme like an E3 ligase, where different F-box factor family members are selectively activated in a disease-specific manner (43). This approach will improve the specificity of anti-BRCA targeting. Here, by retrospective studies with BRCA patient tissues (n=134 from Taiwan BioBank and n=1215 from TCGA database), we report that FBXL8, a novel F-box E3 ligase family member, is significantly correlated with BRCA disease progression. For the first time, we revealed the pathobiological functions of FBXL8 and elucidated its functional significance ex vivo. We found that FBXL8 blocks apoptosis of BRCA cells and modulates cytokine profiles to promote a pro-inflammatory condition. Furthermore, we showed that knockdown of FBXL8 induced intrinsic apoptosis of the BRCA cells, suggesting that FBXL8 promotes BRCA advancement, as supported by the involvement of other F-box proteins in regulating the EMT process in cancers (14). Consistently, we demonstrated that the loss of FBXL8 decreased cell migration and invasion, affirming that FBXL8 is a key promoter of BRCA progression. This is strongly corroborated by the increase in FBXL8 level as the disease advanced in BRCA patients (Figure 2; **, p<0.01; ***, p<0.001). Congruent to these findings, knockdown of FBXL8 significantly down-regulated cancer-promoting cytokines and chemokines. Therefore, repression of FBXL8 is a potential anticancer strategy to activate apoptosis, to inhibit metastasis, and to modulate pro-tumorigenic cytokines in BRCA (Figure 7A).
Further insights on the translational impact of FBXL8 was accrued from in silico modeling prediction, which revealed various potential interaction partners of FBXL8, for example: Cyclin family members like cancer suppressor CCND2 (18), IRF5 (Interferon Regulatory Factor 5), and DNA repair /replication associated proteins (Nibrin and SMAD4) (Figure S5). Thus, FBXL8 appears rather promiscuous, associating with multiple binding partners, to plausibly promote cancer formation and progression through the turnover of: (i) cell-cycle associated proteins, (ii) cancer-suppressors and (iii) proteins responsible for DNA repair/replication. We further validated two of the in silico predicted binding partners, CCND2 and IRF5 (tumor-suppressors), and empirically demonstrated their functional roles which might countermeasure FBXL8, in cancer cells. Notably, the specificity of ubiquitination is mainly controlled at the level of E3 ligase (44), and it is possible that CCND2 and IRF5 are the specific substrates targeted for degradation by FBXL8, an E3 ligase. The molecular axes of interactions of FBXL8-CCND2, FBXL8-IRF5 and/or FBXL8-CCND2-IRF5 prompted us to speculate that their partnership(s) could fine-tune and maintain cellular homeostasis (Figure 7B). Should the balance of this partnership tip over in favor of FBXL8, the tumorigenic potential will rise. Since IRF-5 has been reported to stimulate inflammatory genes (20, 45), the involvement of IRF5 indicates how FBXL8 might indirectly influence the cytokine profiles, hence promoting a pro-inflammatory cytokine-rich cancer microenvironment, as our results alluded to (Figure 5). Additionally, since IRF5 itself can also acts as a pro-tumorigenic factor (46), it is plausible that besides E3, different E2s in the UPS may also respond to the pathophysiological conditions to modulate the FBXL8-IRF5 axis, which should enhance the tumorigenic and metastatic potential.
Since FBXL8-knockdown raised the levels of CCND2 and IRF5, which resulted in cancer suppression, it is conceivable that as a tumor-promoter, FBXL8 acts as a dominant member of a tripartite liaison, in which FBXL8 controls the two tumor suppressors, CCND2 and IRF5 (Figure 6). Of significance is that FBXL8 is a specific protein degradation enzyme which selectively targets cancer suppressors like CCND2 and IRF5. Unlike Bortezomib which is a general protein degradation enzyme, we propose that FBXL8-blockers can be developed to specifically unleash CCDN2 and IRF5, in BRCA treatment. Furthermore, antibodies against the pro-inflammatory cytokines, and/or binding partners identified in this study may be employed clinically in combination with FBXL8-blocker to treat heterogeneous variants of BRCA.
Data Availability
All data generated or analysed during this study are included in this published article (and its supplementary information files).
Author Contributions
S.C.C. designed, conducted the experiments and wrote the paper; W.H. and C.Y.S. verified the analytical methods; C.S.H. provided the human primary samples and designed the experiments; J.L.D. designed the experiments and wrote the paper.
Competing interests
The authors declare no conflict of interest.
Materials and Methods
Tissue samples
67 pairs of tissue samples (tumor and normal tissues; n=134) from BRCA patients were acquired from the BioBank, Taipei Medical University, Taiwan, after clinical diagnosis was confirmed by biopsy and histological evaluation. All patient samples were collected with informed consent. The study was approved by the TMU-JIRB (Taipei Medical University-Joint Institutional Review Board; IRB: N201803107). Adjacent non-cancerous breast tissues were obtained at least 2 cm away from the tumor node, to serve as paired normal tissue control. Ex vivo experiments included (i) NGS-based RNA-seq analysis, where five-pairs of the primary tissues (n=10) were obtained from: one stage I, three stage II and one stage III patient samples (the corresponding clinicopathological information are shown in Figure 1A); (ii) qRT-PCR analysis of 32 paired tissues (n=64), including fifteen stage I, nine stage II and eight stage III patient samples (the corresponding clinicopathological information are shown in Supplementary Figure S1); and (iii) immunofluorescence staining, where 30 pairs of the tissues (n=60) were used, including: nine stage I, eleven stage II and ten stage III patient samples (the corresponding clinicopathological information are shown in Supplementary Figure S2). H&E sections were obtained from the BioBank, Taipei Medical University, Taiwan. For immunofluorescence staining, paraffin-embedded tissues were sectioned at 5 μm thickness, dehydrated and blocked in Ultravision Protein Block (Thermo Scientific TA060PBQ) for 10 min at room temperature. Then, the sections were incubated with primary and secondary antibodies overnight and 1 h, respectively. The stained tissue sections were observed under a microscope (Olympus, Tokyo, Japan) and images were acquired using software (EOS Utility, Canon).
Cell lines and reagents
Human BRCA cell lines (MCF7, MDA-MB231) and an FBXL8-negative control breast cell line (MCF10A) were obtained from American Type Culture Collection (ATCC). MCF7 and MDA-MB231 were cultured in complete DMEM medium (Gibco). MCF10A was cultured in complete DMEM/F12 medium (Gibco). All the complete media were supplemented with 10 % FBS (Thermo Scientific) and 100 U/ml penicillin and 100 μg/ml streptomycin (Invitrogen). Antibodies used in IHC, immune blotting and immunoprecipitation analysis were anti-FBXL8 antibody (sc-390582, Santa Cruz Biotechnology, Dallas, TX, USA), anti-ubiquitin antibody (ab7780, Abcam, Cambridge, UK), anti-CCND2 antibody (ab226972, Abcam, Cambridge, UK), anti-IRF5 antibody (ab2932, Abcam, Cambridge, UK), normal mouse IgG2b (sc-3879, Santa Cruz Biotechnology, Dallas, TX, USA) and anti-Tubulin antibody (ab176560, Abcam, Cambridge, UK).
Next-generation sequencing (NGS)-based RNA-seq analysis
RNA isolation for NGS-based RNA-seq. Total RNA was extracted by Trizol® Reagent (Invitrogen, USA) according to the instruction manual. Purified RNA was quantified at OD260nm by using Qubit (Life Technologies, Qubit® 2.0 Fluorometer, USA) and analysed using a bioanalyzer (Agilent, Bioanalyzer 2100 system, USA) with RNA 6000 labchip kit (Agilent Technologies, USA).
Library preparation & sequencing for NGS-based RNA-Seq application. All procedures were carried out according to the manufacturer’s protocol from Illumina. Library construction of all samples was performed using Agilent’s SureSelect Strand Specific RNA Library Preparation Kit (100SE, Single-End) bp sequencing on Solexa platform. The sequences were determined using sequencing-by-synthesis technology via the TruSeq SBS Kit. Raw sequences were obtained from the Illumina Pipeline software bcl2fastq v2.0, which was expected to generate 5M (million reads) per sample.
RNA sequence analysis. Quality control of raw sequence data was performed with FastQC. Adapter and quality trimming was performed using cutadapt (phred cutoff: 20; minimum length: 25) (47). Reads were mapped to the human reference genome using STAR aligner indexed with GRCh37 assembly and GENCODE Release 26 v26lift37 annotation. Gene expression was quantified using RSEM. The expected counts were used as input for the differential gene expression analysis between tumor and adjacent-normal BRCA samples using edgeR in the R environment. A false discovery rate (FDR) <0.05 was used as a threshold to define genes that showed statistically significant differential expression. A log2 fold change (logFC) less than 0 indicates under-expression in the tumor samples, whereas logFC greater than 0 represents over-expression in the tumor samples.
TCGA database analysis
TCGA breast cancer transcript dataset was downloaded from R package, TCGAbiolinks (48). In total, 1215 primary tissues were analyzed, including 1102 carcinoma tissues and 113 normal tissues. All of the raw counts were used as input for the differential gene expression analysis between tumor and normal breast samples, using DESeq2. The P adjust value < 0.05 was used as threshold to define genes that showed statistically significant differential expression.
Quantitative real-time PCR
24 h after treatment of cell lines, total RNA was extracted with TRIzol reagent (Invitrogen) and the RNA was reverse transcribed using random hexamers and HiSenScript(tm) RH(-) RT-PCR PreMix Kit (iNtRON Biotechnology, Seongnam, South Korea), according to the manufacturer’s instructions. The synthesized cDNA was used for real-time PCR analysis using SYBR Green (Life Technologies) on a QuantStudio 3 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). The primers were human FBXL8 (83 bp product): sense, 5’-AATCAGTTGCGAATGTGAGC-3’; antisense, 5’-CCAGCCGTAGGTTGTGAATG-3’; human SAG (109 bp product): sense, 5’-CAGGCTCCAAGTCGGGAGGCG-3’, antisense, 5’-TGGACCCTGCAGATGGCACAGG-3’; human GAPDH (131 bp product): sense, 5’-GTCTCCTCTGACTTCAACAGCG-3’, antisense, 5’-ACCACCCTGTTGCTGTAGCCAA-3’. All expression values were normalized based on GAPDH as an endogenous control.
FBXL8 knockdown in BRCA cells
FBXL8 siRNA (siGENOME Human FBXL8-SMARTpool, M-017504-00-0005) was purchased from Dharmacon (Lafayette, CO, USA). Control (scrambled) siRNA were purchased from Invitrogen. siGENOME Human FBXL8-SMARTpool contains a mixture of 4 siRNAs provided as a single reagent; providing advantages in both potency and specificity, which was anticipated to achieve silencing effects > 75%. siRNA transfection into MCF7 or MDA-MB231 cells (at 5×105 cells per well of a 6-well plate) was conducted using 7.5 μl TransIT-X2 Transfection Reagent (Mirus) with 25 nM of siRNAs per well.
Cell viability assay
Both MTT {3-(4,5-dimethylathiazol-2-yl)-2,5-diphenyl tetrazolium bromide} and CTB (CellTiter Blue) obtained from Life Technologies (Carlsbad, CA, USA) and Promega (Madison, WA, USA), respectively, were used to determine the cell viability, according to the manufacturer’s instructions. To assess the cell viability the functional mitochondrial activity was measured by incubating the cells with either 10 μl MTT reagent in 100 μl culture medium at 37°C for 4 h, or 10 μl CTB reagent at 37°C for 4 h. The metabolic activity of viable cells was measured at 570 nm using a Synergy H4 microplate reader (BioTek, Winooski, VT, USA) at time points indicated. Samples from each time point were normalized with corresponding PBS controls.
Cell proliferation assay
After cells are transfected with FBXL8-specific siRNA or control scrambled siRNA for 16 h, both Alamar blue and Trypan Blue dye exclusion were used to measure the cell proliferation. The cell growth was measured by Alamar blue. The metabolic activity of the cells was determined at 570 nm using a Synergy H4 microplate reader (BioTek, Winooski, VT, USA) at time points indicated. For Trypan Blue dye exclusion, to assess total cell number, cells were scraped and resuspended in equal volumes of culture medium and trypan blue dye (0.4% solution; Gibco, Waltham, MA, USA) and counted using an improved Neubauer hemocytometer.
Cellular apoptosis assay
MCF7 and MDA-MB231 cells were transfected with FBXL8-specific siRNA or control scrambled siRNA for 24 h, followed by determination of apoptosis. Early apoptosis was measured using annexin V (BioLegend, San Diego, CA, USA) in conjunction with 7-AAD (BioLegend), according to manufacturer’s instructions. Annexin V identifies surface-exposed phosphatidylserine, and 7-AAD is retained in late apoptotic cells. Early apoptosis (Annexin V+/7-AAD-) was then analyzed using FACScan flow cytometer (FACSVerse, BD Biosciences, Franklin Lakes, NJ, USA) for a minimum of 10,000 events.
Caspases-9 and −3 assays
24 h after the transfection of FBXL8-specific siRNA or control scrambled siRNA, apoptosis of MCF7 and MDA-MB231 cells were examined by determining the caspase-specific cleavages of activated caspases-9 and −3 (Abcam). The p-nitroanilide (p-NA) light emission was measured at 405 nm using a Synergy H4 microplate reader (BioTek).
Cell migration and invasion assays
24 h after treatment of the MCF7 and MDA-MB231 cells with FBXL8-siRNA or controls, cell migration assay was performed. The cells were rinsed with PBS, harvested using 0.05 % Trypsin-EDTA (Invitrogen) and plated into a 2-well Culture-Insert (ibidi, Martinsried, Germany) according to the manufacturer’s instructions. The Culture-Inserts were removed after 16 h, and cell migration was monitored over indicate time-scale. For cell invasion assay, biocoat matrigel invasion chambers with 8-μm pores in 24-well plates (Corning, Discovery Labware, Inc., Bedford, MA, USA) were used. 24 h after siRNA transfection, the cells were plated onto the chambers. The cells were detached with 0.05 % Trypsin-EDTA, resuspended in conditioned medium (10 % FBS) and added to the upper compartment of the chambers, according to the manufacturer’s instructions. After 24 h of incubation at 37 °C, the cells on the upper surface were completely removed by wiping with a cotton swab, and then the filter was fixed with 100 % methanol and stained with crystal violet solution (0.5 % (w/v) crystal violet in 25 % (v/v, methanol). Cells that had migrated from the upper to the lower side of the filter were imaged and counted with a light microscope (5 fields/filter).
Cytokine arrays
The secretion of cytokines from BRCA cells effected by FBXL8-knockdown was analyzed using human cytokine antibody array (human cytokine antibody array C6 and C7, AAH-CYT-1000) purchased from RayBiotech (Peachtree Corners, GA, USA). Cells were cultured in complete medium with transfection of FBXL8-specific siRNA or control scrambled siRNA for 48 h. Culture supernatants were collected for cytokine antibody array analysis. Cytokine signal intensities were quantified by Scion Image software (Frederick, MD, USA). The fold-change of cytokines secreted from cells was calculated by normalization of FBXL8-specific siRNA-treated samples to control scrambled siRNA-treated samples.
ELISA quantification of cytokines
To determine the effects of FBXL8 on cytokine release, both MCF7 and MDA-MB231 were transiently transfected with FBXL8-specific siRNA or control siRNA for 48 h. The levels of MCP-1, EGFR, I-TAC, TECK, CTACK, MIF, GM-CSF and ICAM-1 (based on the cytokine array results), secreted by the cells, were quantified in supernatants derived from 48-h cultures, by using the respective ELISA (R&D systems, Minneapolis, MN, USA), according to manufacturer’s instructions.
Co-IP and immunodetection
For Co-IP study with MCF7 or MDA-MB231 total cell lysates, the cells were lysed in CHAPS lysis buffer containing 20 mM Tris-HCl, pH 7.5, 5 mM MgCl2, 137 mM KCl, 1 mM EDTA, 1 mM EGTA, 1% CHAPS and 1 protease inhibitors (complete EDTA-free cocktail, Roche). Then, the cell lysate was precleared by incubating with protein G Sepharose (GE Healthcare, Little Chalfont, UK) at 4 °C for 2 h. The supernatant was incubated overnight at 4 °C with FBXL8 or control IgG2b antibody, followed by a 3-h incubation with protein G Sepharose. The washed immunoprecipitates, resuspended in Laemmli buffer, was boiled at 95 °C for 5 min before immunodetection of FBXL8, CCND2 or IRF5.
Statistical analysis
Data were expressed as means ± S.D. from three independent experiments, with three replicates per sample/condition tested. Differences between averages were analyzed by two-tailed Student’s t-test. Significance was set at p-value of <0.01 (**, p<0.01; ***, p<0.001). The acquired data from FACS were analyzed with BD FACSuite™ software (BD Biosciences, CA, USA). The relative migration rate, indicated as % gap closure, was calculated using Image J analysis software. All target signals from IHC were quantified by HistoQuest software (TissueGnostics, Vienna, Austria).
Acknowledgements
The authors thank Dr. I-Hsuan Lin and Dr. Tsung-Han Hsieh, TMU Cancer Translational Core Facility for the NGS-based RNA-seq technology support. This work was jointly supported by the Ministry of Science and Technology, Taiwan (MOST105-2320-B-038-007-MY2, MOST107-2314-B-038-066); Taipei Medical University (TMU104-AE1-B10, 105TMU-TMUH-06); the Ministry of Education (R-154-000-A76-114) and the National University of Singapore Graduate School of Integrative Sciences and Engineering (NGS, C-154-000-017-091), Singapore.
Footnotes
↵Σ Co-senior authors: sc.chang{at}tmu.edu.tw, dbsdjl{at}nus.edu.sg