Abstract
Breast cancer is known to be a heterogeneous disease driven by a large repertoire of molecular abnormalities, which contribute to its diverse clinical behaviour. Despite the success of targeted therapy approaches for breast cancer patient management, there is still a lack of the molecular understanding of aggressive forms of the disease and clinical management of these patients remains difficult. The advent of high-throughput sequencing technologies has paved the way for a more complete understanding of the molecular make-up of the breast cancer genome. As such, it is becoming apparent that disruption of canonical splicing within breast cancer governs its clinical progression. In this review, we discuss the role of dysregulation of spliceosomal component genes and associated factors in the progression of breast cancer, their role in therapy resistance and the use of quantitative isoform expression as potential prognostic and predictive biomarkers with a particular focus on oestrogen receptor-positive breast cancer.
Introduction
Dysregulation of alternative splicing (AS) is widely considered a new hallmark of cancer and its products are being acknowledged as potentially useful biomarkers (Ladomery 2013). Canonical RNA splicing takes place in all mammalian cells and during this process, pre-mRNA becomes mature mRNA via the excision of introns and pasting together of exons (Fig. 1). AS affects about 90% of human genes resulting in a diverse selection of isoforms from one gene, each having different structural and functional properties that lead to a larger and more diverse cellular proteome. Indeed throughout evolution, AS has been used to propel species development evidenced by an increase in AS in higher eukaryotes compared to lower (Keren et al. 2010).
Splicing is performed by the spliceosome, which is a multi-protein complex called a ‘metalloribozyme’ that is made up of five small nuclear riboproteins (snRNPs) that contain snRNAs and a large number of accessory proteins to recognise the pre-mRNA being spliced. Assembly of this complex takes place during transcription suggesting that transcription and splicing machineries are space restricted as they happen closely in time (Herzel et al. 2017). The most commonly occurring spliceosome is the U2-dependent spliceosome that is assembled from the U1, U2, U5 and U4/U6 snRNPs and is responsible for the splicing of 99% of human introns as reviewed by Dvinge et al. (2016). Splicing is a two-step reaction involving transesterification occurring between two RNA nucleotides. The spliceosome recognises introns containing the consecutive nucleotides GU at the 5′ splice site (SS) by U1 snRNP binding and an AG sequence at the 3′ SS by U2AF1 binding. In order to properly position the splicing machinery a key adenine (also referred to as the branch point (BP)) must be recognised by the splicing factor 1 protein (SF1) as well as recognition of the polypyrimidine tract (poly-Y) by the U2 small nuclear RNA auxiliary factor 2 (U2AF2) (Pandya-Jones 2011). When the spliceosome complex is correctly bound to the mRNA, it can carry out intron excision (Fig. 1A). Exons are subsequently joined together and release a lariat intron, which is then degraded. The spliceosome components are then released and recycled for use in subsequent rounds of splicing. Splicing factors such as serine/arginine-rich (SR) proteins (SRSF) and the splicing factor 3b complex (SF3B) work in association with the splicing core complex to coordinate canonical and AS. The expression levels and binding affinities of the different splicing factors play a stoichiometric role in determining the final isoform of the protein that is to be expressed (da Luz et al. 2017).
Dysregulation of the normal splicing process governs many aspects of cancer cell biology such as managing cellular proliferation, angiogenesis, resisting apoptosis, adapting cell metabolism, enhancing the ability to invade and metastasise and plays a role in resistance to cancer therapy (David & Manley 2010, Lee & Abdel-Wahab 2016). The role of AS in disease can result from aberrant splicing of a gene due to incorrect 5′ or 3′ SS recognition leading to intron retention, exon skipping or exon inclusion (Fig. 1A). AS may then lead to a premature stop codon resulting from a frame-shift, whereby these transcripts are subsequently degraded by nonsense-mediated decay (NMD) (Fig. 1B). There are multiple ways in which a cancer cell can induce aberrant splicing including: (1) when there is a mutation in the exon or surrounding introns that compromises the canonical splicing signal thereby allowing an alternative signal to dominate and an aberrant mRNA to be made; (2) a mutation in one of the splicing regulators interrupts SS selection and results in a pattern of AS in multiple genes; (3) changes in histone acetylation of alternative exons (Khan et al. 2014) and (4) alterations in other RNA-binding proteins, splicing enhancers and suppressors or lncRNAs. Such splicing errors can lead to alterations in relative isoform expression of a particular mRNA or lead to an aberrant protein that has a change of function. A more detailed discussion on points 1 and 3 are detailed elsewhere (Martinez-Montiel et al. 2017). Aberrantly spliced apoptotic genes such as the RNA-binding protein RBM5 have been implicated in breast cancers as having an opposing role because the resulting isoform is more anti-apoptotic (Fushimi et al. 2008). Another example is the B-cell lymphoma gene, Bcl-x, which can be spliced into two different isoforms, long and short. Bcl-x(L) has anti-apoptotic properties whereas Bcl-x(s) has pro-apoptotic properties. High levels of Bcl-x(L) are seen in various types of cancer (Boise et al. 1993, Takehara et al. 2001, Fushimi et al. 2008). A similar situation is seen with the myeloid cell leukaemia-1 gene and its two isoforms MCL-1(S) and MCL-1(L). The long isoform is anti-apoptotic and seen frequently increased compared to the short isoform in breast and ovarian cancer cells and is linked to gene amplification of MCL-1 itself (Bae et al. 2000, Bingle et al. 2000, Gautrey & Tyson-Capper 2012). The choice between the long and short isoform is influenced by the splicing factors SRSF1 and SRSF5, which are also frequently upregulated in breast cancer (Gautrey & Tyson-Capper 2012).
Managing key cellular processes such as epithelial-to-mesenchymal differentiation (EMT) is a clear advantage of being able to manipulate the expression of different isoforms of a certain gene (Shapiro et al. 2011). As such, acquiring the ability to hijack these processes is critical in the evolution of a cancer cell in order to provide a fitness advantage. Given this, it is reasonable to postulate that the characterisation of the splicing programme of a cancer cell could predict its genomic and mutational status and potentially treatment outcome (Danan-Gotthold et al. 2015). Indeed, differential expression of AS transcripts in specific subtypes of breast cancer may add additional prognostic information in addition to canonical gene expression or protein expression biomarkers.
Evidence of splicing dysregulation in breast cancer
Since the seminal studies from Perou and colleagues describing the intrinsic subtypes of breast cancer (Perou et al. 2000), it is now widely accepted that the molecular make-up of breast cancer is heterogeneous and governed by differences in transcriptional make up. Inevitably, this also applies to the degree of isoform usage in cancer cells as well. For instance, well-known driver oncogenes and tumour suppressor genes such as ERBB2 and BRCA1 are known to be differentially spliced in different subtypes of breast cancer (as reviewed by Martinez-Montiel et al. 2017). BRCA1, which is involved in homologous recombination DNA repair, is alternatively spliced in breast cancer to exclude exon 11 that contains the nuclear localisation signal (Thakur et al. 1997). The ∆11q isoform produces a protein that is absent from the nucleus and is therefore unable to assist in DNA damage repair. Studies have shown that downregulation of the full-length nuclear BRCA1 isoform and overexpression of the cytoplasmic ∆11q isoform is evident in subsets of breast cancer and is potentially mediated through the presence of a non-functional TRA2β splicing factor (Raponi et al. 2014, Wiener et al. 2015). Another example is the ERBB2 tyrosine kinase signalling receptor, which is often found as alternatively spliced in breast cancer as the ∆16HER2 isoform. ∆16HER2 is constitutively active as a homodimer and promotes transformation in the mammary gland (Marchini et al. 2011). BRCA1 and ERBB2 splicing, as well as splicing of BCL-X and MCL-1 as described earlier, are examples of common driver oncogenes and tumour suppressor genes that can be aberrantly spliced in breast cancer. AS has also been shown to regulate protein diversity of the oestrogen receptor itself. In particular, previous studies have shown the ERαΔ5 splice variant has a positive effect on activation of transcription in the absence of oestrogen leading to constitutive transcriptional activation (Fuqua et al. 1991, Bollig & Miksicek 2000). ESR1 aberrant splicing events have also been identified in circulating tumour cells from metastatic breast cancer patients that have progressed on endocrine therapy, suggesting a role in mediating resistance (Beije et al. 2018). Current data sets describing AS events in the context of spliceosomal gene mutations, however, do not show changes in splicing of the oestrogen receptor itself (Darman et al. 2015, Maguire et al. 2015). Alternatively spliced isoforms of genes known to be transcriptionally regulated by the oestrogen receptor such as Cyclin D1 (cyclin D1b) and FGFR1 (FGFR1-beta) are also associated with poor prognosis in ER+ breast cancer (Wei et al. 2011, Wendt et al. 2014).
Alternatively-spliced transcripts as prognostic and predictive biomarkers in breast cancer
The recent advent of RNA-sequencing technologies has revolutionised our view of the molecular make up of breast cancer. These advances now allow accurate global quantification of the transcriptional isoform make-up in individual tumours rather than relative quantification that is based on microarray probe design. Indeed, a number of studies have shown that alternative isoform usage can be specific to different breast cancer molecular subtypes (Menon et al. 2014, Sebestyen et al. 2015, Zhao et al. 2016, Gracio et al. 2017, Stricker et al. 2017). For instance, Sebestyen et al. identified a specific 7 gene isoform signature that accurately identified basal-like breast cancers, including a number of known driver genes such as CTNND1 (Sebestyen et al. 2015). Analysis of the splicing balance (relative ratios of isoforms produced) in breast tumours revealed changes in isoform usage in oncogenic and tumour suppressive pathways that was not apparent when looking solely at gene expression data (Gracio et al. 2017). Importantly, it was found that the balance of different transcript isoforms was associated with patient prognosis. A subset of genes including the proto-oncogene MYB were identified to correlate with basal-like breast cancer patient survival based on varying isoform levels but not on whole gene expression analyses (Gracio et al. 2017). Additionally, splicing but not gene expression levels of immune-related genes CCR7 and FCRL3 were found to determine the immune control of the tumour. This has potential relevance given the role of lymphocytic infiltration in prognosis in breast cancer. Differential isoform usage can also stratify between different molecular subtypes of breast cancer. Indeed, global dysregulation of splicing specific to individual subtypes may drive the heterogeneous nature of breast cancer due to variation in the cellular proteome. Stricker et al. (2017) looked at the global isoform differences between ER+ and triple-negative breast cancer (TNBC) and identified a signature of subtype-specific alternatively spliced transcripts. Interestingly, around 63% of the genes that were found to be differentially expressed, between subtypes were also alternatively spliced. The particular type of splicing that occurred between the subtypes (exon skipping, intron retention, alternative acceptor or donor), however, was not significantly different indicating the unique splicing programmes of each intrinsic subtype is not necessarily due to the activity of one general splicing mechanism but more likely due to target gene selection (Stricker et al. 2017). Interestingly, this study also identified a significant difference in the total expression of some spliceosomal component genes themselves, such as YBX1 and MAGOH suggesting dysregulation of spliceosomal component proteins governs splicing dysregulation.
Although clear differences in transcript isoforms have been identified in different molecular subtypes of breast cancer, to date, no study has assessed the value of alternatively spliced transcripts as prognostic and predictive clinical biomarkers for patient stratification and of treatment response to both standard chemotherapy and targeted endocrine therapy. Assessment of differences in transcript isoform expression could add much needed biomarkers for patients who are most likely to relapse on standard-of-care therapy. Ideally, this would need to be tested in the context of randomised clinical trial cohorts, where good-quality RNA-sequencing data at sufficient depth are acquired.
Dysregulation of spliceosomal factors in breast cancer
Molecular alterations affecting spliceosomal component genes themselves are also known to be involved in breast cancer tumourigenesis. There is evidence that mutations, copy number alterations and differential expression of spliceosomal component genes and their interacting proteins are associated with specific molecular and histological subtypes of breast cancer as well as being associated with aggressive disease and resistance to therapy in multiple tumour types (Stark et al. 2009, Ng et al. 2012, Sotillo et al. 2015, Siegfried & Karni 2017). These alterations are thought to drive breast cancer progression through specific or novel isoform selectivity of key genes (Vanharanta et al. 2014, Anczukow et al. 2015, Gokmen-Polar et al. 2015, Maguire et al. 2015, Silipo et al. 2015, da Luz et al. 2017, Martinez-Montiel et al. 2017).
Mutations in spliceosomal component genes
Mutations affecting different components of the spliceosome have been identified in a range of solid and non-solid malignancies (Papaemmanuil et al. 2011, Quesada et al. 2011, Biankin et al. 2012, Furney et al. 2013, Yoshida & Ogawa 2014). Mutations in the splicing factor SF3B1 are the most common across multiple tumour types, and are found at particularly high frequencies in myelodysplastic syndrome (MDS), chronic myeloid leukaemia (CLL), uveal melanoma (UV), pancreatic cancer and breast cancer. Mutations generally cluster at hotspot amino acid residues K700, R625, K666 and H662 (Cerami et al. 2012, Gao et al. 2013). However, each cancer type harbours a different variation of hotspot mutations. For example, K700E mutations are invariably found in breast cancer, pancreatic cancer and CLL, whereas UV and endometrial cancers harbour the R625, R666 and R662 hotspots, suggesting some tissue specificity of the mutations. SF3B1 hotspot mutations in CLL are associated with a poor prognosis. However, in UV and MDS the prognosis is better with the presence of an SF3B1 mutation (Quesada et al. 2011, Furney et al. 2013). Interestingly, additional spliceosomal component genes are also recurrently mutated at high frequencies particularly in MDS, including U2AF1, which has a distinct S34F/Y hotspot mutation and mutations in SRSF2 that are associated with a poor outcome in MDS (Thol et al. 2012) and ZRSR2. Both SRSF2 and ZRSR2 harbour mutations spread throughout the gene, suggestive of a tumour suppressive function (Yoshida et al. 2011). Mutations in these genes including SF3B1 occur in a mutually exclusive manner in MDS, suggesting that cells may tolerate only a partial deviation from normal splicing activity. Indeed, these genes are all involved in the 3′-SS recognition during pre-mRNA processing, inducing abnormal RNA splicing and compromised haematopoiesis (Yoshida et al. 2011), implicating splicing dysregulation as a major driving force behind the development of MDS.
Our group has explored the mutational repertoire of spliceosomal component genes in breast cancer from a meta-analysis of whole genome and exome sequencing data (Maguire et al. 2015) (Fig. 2). This analysis identified that around 5.6% of unselected breast cancers have mutations in spliceosome component genes at low frequencies. The most common spliceosomal gene mutation is SF3B1, which is associated with ER+ breast cancer and seen in around 3% of ER+ tumours (Pereira et al. 2016), whereas mutations in SON and SAP130 appear to be associated with ER− disease (Maguire et al. 2015). Interestingly, we identified SF3B1 K700E mutations at higher frequencies in some rarer histological subtypes of breast cancer including 16% of papillary carcinomas and 8% of mucinous carcinomas of the breast, suggesting they may underpin their biology (Maguire et al. 2015). SF3B1 K700E mutations were also found to associate with losses of 16q11-q13 and gains of 16q12-q13 indicating a distinct mechanism of breast cancer progression independent of the canonical early event of 1q gain and 16q loss (Maguire et al. 2015).
The association of SF3B1 mutations and breast cancer clinical prognosis, however, is unclear, although mutations are being increasingly seen in metastatic disease (Lefebvre et al. 2016, Pereira et al. 2016). Further studies however are needed in order to truly assess the effect of SF3B1 hotspot mutations on outcome. Of note, SF3B1 mutations have been observed in adenoid cystic carcinomas of the breast (an ER-negative special histological subtype) that has an excellent clinical outcome and at increased frequency in ER+ mucinous and papillary carcinomas of the breast. These data perhaps suggest that SF3B1 mutations maybe associated with a good prognosis (Maguire et al. 2015, Martelotto et al. 2015).
SF3B is a complex that is part of the U2 spliceosome and controls 3′ SS recognition. Its core is required for alignment of the branch site proteins, which allows for correct branch site selection during the splicing process (Cretu et al. 2016). SF3B1 (SF3B155) is the largest component of the SF3B complex and contains the HEAT superhelix domain consisting of 20 tandem repeats of two alpha helices joined by a short loop. Mutations in the HEAT domains, which are responsible for interacting with pre-mRNA and other pre-mRNA-binding proteins, result in a change in the tertiary structure that causes the selection of an alternative branch site (Darman et al. 2015, Alsafadi et al. 2016, Kesarwani et al. 2017). It is not known, however, whether the SF3B1 mutant protein has a stronger affinity for the newly exposed BP sequence or if it is coping with a disruption in binding to the canonical BP sequence (Darman et al. 2015). Indeed, mutations in SF3B1 lead to alternative branchpoint usage and subsequent usage of a 3′ cryptic SS. This leads to aberrant transcript expression and subsequent NMD of around half the aberrant transcripts and hence leads to protein downregulation (Darman et al. 2015, Alsafadi et al. 2016, Kesarwani et al. 2017).
Although present as hotspot single amino acid changes, SF3B1 mutations are thought to lead to a change in function. This is because knockdown or overexpression of the mutant protein does not recapitulate the aberrant splice pattern seen in mutant vs WT patients (Alsafadi et al. 2016). Additional evidence suggests that these mutations may actually be loss of canonical function. For instance, using the Degron-knock-in approach to inactivate mutant or WT alleles specifically, Zhou et al. found that degradation of only the mutant SF3B1 allele in heterozygous SF3B1-mutant cells had no effect on growth, whereas degradation of only the WT allele resulted in a decrease in viability of the cells (Zhou et al. 2015). This suggests that SF3B1 is not likely to be haploinsufficient given the cells are solely relying on the WT copy of the gene to survive. This observation helps explain why SF3B1 mutations are uniformly heterozygous, as two copies of the mutant allele would likely be lethal.
The most common SF3B1 mutation in breast cancer is the K700E variant akin to CLL but K666Q and K666E are also observed, albeit at much lower frequencies (Maguire et al. 2015). Gene expression analysis in ER-positive disease shows that SF3B1 mutations affect regulators of the cell cycle, metabolism and motility as well as protein degradation and apoptosis, and splicing regulation itself (Maguire et al. 2015). Commonly differentially spliced mRNAs have been associated with SF3B1 mutations across tumour types including UV, chronic lymphocytic leukaemia, pancreatic cancer and breast cancer. Although a large number of transcripts have been identified to be aberrantly spliced and some are cancer specific (e.g. ABCB7 AS is only observed in MDS and gives rise to increased mitochondrial iron accumulation found in MDS patients with ring sideroblasts (Dolatshad et al. 2016), the overlap is rather strikingly consistent between tumour types, suggesting that there is a distinct signature of genes that are alternatively spliced and furthermore can be used as markers of the mutation status (Quesada et al. 2011, Biankin et al. 2012, Furney et al. 2013, Maguire et al. 2015, Dolatshad et al. 2016, Obeng et al. 2016, Wang et al. 2016). However, it has not yet been identified which of the many differentially spliced genes is/are responsible for the tumorigenic phenotype and if these are different between different cancer types. In our study, we used siRNA to silence different genes that had been identified as alternatively spliced in our data set as well as across multiple cancer types. Silencing eight different genes (ABCC5, ANKHD1, DYNLL1, F8, RPL31, TMEM14C, UQCC and CRNDE) did not show any changes in viability (Maguire et al. 2015). Given around half of all aberrantly expressed transcripts are subjected to NMD, they could be acting as tumour suppressors rather than in an oncogenic manner and will need to be explored in the future.
Spliceosomal component genes as oncoproteins in breast cancer
As well as mutations, alterations in components of the spliceosome, such as deletions or amplifications, are commonly seen across breast cancer (Fig. 2; Tables 1 and 2). In a similar vein to spliceosomal component mutations, they may lead to dysregulation of canonical splicing. SF3B3 (SF3B130) a component of the SF3B complex has been found to be significantly overexpressed in ER+ breast cancers and is associated with aggressive disease and resistance to tamoxifen therapy (Gokmen-Polar et al. 2015). SF3B3 is positioned closely to SF3B1 in the U2 complex and helps maintain the HEAT domain’s structural plasticity and has the ability to alter pre-mRNA splicing hence affecting gene expression in the cell (Garcia-Blanco et al. 2004). Overexpression of SF3B3 has thus been postulated to contribute to spicing aberrations in cancer cells. In clear cell renal cell carcinoma, SF3B3 overexpression was found to increase the expression of the pro-proliferative full-length isoform of EZH2 and not the commonly expressed EZH2∆14 that is found in normal tissue (Chen et al. 2017), thus promoting tumourigenicity in vivo. It could be that EZH2 AS plays a role in mediating the aggressive behaviour in endocrine resistant ER+ breast cancer as well; however, this is yet to be elucidated. SF3B3 has also been found to be amplified and highly expressed at the transcript level in basal-like breast cancers (Srihari et al. 2016). Overall, the level is actually higher in ER− than ER+ disease, perhaps highlighting the higher proliferative rate of these tumours.
Summary of spliceosome component genes and RNA-binding proteins found altered in breast cancer.
Splicing factor/RNA-binding protein | Gene name | Alteration | Occurrence in BrCa (%) | Functional impact |
---|---|---|---|---|
SF3B1 | Splicing factor 3B subunit 1 | Mutation and CNA | 3 | Change of function, oncogenic |
SF3B3 | Splicing factor 3B subunit 3 | CNA | 1.7 | Oncogenic |
SRSF1 | Serine/arginine rich splicing factor 1 | CNA | 8 | Oncogenic |
SRSF2 | Serine/arginine rich splicing factor 2 | CNA | 6 | Oncogenic |
SRSF3 | Serine/arginine rich splicing factor 3 | CNA | 1.1 | Oncogenic |
SRSF4 | Serine/arginine rich splicing factor 4 | CNA | 0.6 | Oncogenic |
RBFOX2 | RNA-binding protein fox-1 homolog 2 | CNA | 0.7 | EMT regulator |
ESRP1 | Epithelial splicing regulatory protein 1 | CNA | 18 | EMT regulator |
RBM47 | RNA-binding motif protein 47 | CNA | 1.4 | Downregulation |
LIN28A | Lin-28 Homolog A | CNA | 0.4 | Loss of function |
Sourced from all breast cancer studies available in cBioportal. n = 4587 sequenced cases.
CNA, copy number alteration.
The SRSF family of proteins are serine-arginine-rich splicing factors that are commonly found to be mutated or dysregulated in cancer (Das & Krainer 2014). These proteins contain RNA recognition motif (RPM) domains that contact the mRNA and also interact with other splicing machinery (Das & Krainer 2014). SRSF1 also referred to as SF2/ASF is the most common protein of this family to play a role in breast cancer and overexpression is associated with a poor prognosis in ER+ breast cancers (Anczukow et al. 2012). Overexpressing SRSF1 in 3D mammary organotypic assays is associated with larger acini structures indicating its oncogenic phenotype (Anczukow et al. 2012). This study also highlighted specific isoform dysregulation of the tumour suppressors BIM and BIN1, which resulted in loss of their pro-apoptotic functions (Karni et al. 2007, Anczukow et al. 2012). SRSF1 upregulation is thought to play a role in EMT through AS modulation of its transcriptional target genes (Valacca et al. 2010). Mechanistically, this is linked back to the splicing regulator Sam68, which modulates levels of SRSF1 (Valacca et al. 2010). It was found that SRSF1 is more likely to facilitate exon inclusion when it binds closer to the 5′ site of the splice junction and promotes exon skipping or inclusion when it binds to the 3′ end (Anczukow et al. 2015). SRSF1 was found to alternatively splice CASC4 by including exon 9, resulting in a longer protein. When tested alone, overexpression of this isoform of CASC4 phenocopied the tumorigenic abilities of SRSF1 overexpression by increasing proliferation and acinar size and decreasing apoptosis (Anczukow et al. 2015). These data highlight promising targets for therapeutic development in patients with SRSF1 overexpression.
Number and percentage of patients pertaining to each subtype with an alteration in the specified spliceosome component genes.
SF3B1 | SF3B3 | SRSF1 | SRSF2 | SRSF3 | SRSF4 | RBFOX2 | ESRP1 | RBM47 | LIN28A | |
---|---|---|---|---|---|---|---|---|---|---|
Basaln = 391 | 12 (3.07) | 10 (2.56) | 11 (2.81) | 22 (5.63) | 15 (3.84) | 6 (1.53) | 5 (1.28) | 84 (21.48) | 13 (3.32) | 1 (0.26) |
Her2n = 287 | 9 (3.14) | 3 (1.05) | 39 (13.59) | 25 (8.71) | 2 (0.70) | 1 (0.35) | 2 (0.70) | 99 (34.49) | 6 (2.09) | 1 (0.35) |
Luminal An = 909 | 40 (4.40) | 17 (1.87) | 31 (3.41) | 26 (2.86) | 5 (0.55) | 2 (0.22) | 2 (0.22) | 103 (11.33) | 5 (0.55) | 3 (0.33) |
Luminal Bn = 590 | 19 (3.22) | 5 (0.85) | 99 (16.78) | 50 (8.47) | 3 (0.51) | 2 (0.34) | 4 (0.68) | 159 (26.95) | 2 (0.34) | 1 (0.17) |
Normal liken = 179 | 4 (2.23) | 3 (1.68) | 10 (5.59) | 7 (3.91) | 0 (0.00) | 1 (0.56) | 0 (0.00) | 23 (12.85) | 1 (0.56) | 0 (0.00) |
Data were derived from METABRIC and TCGA samples with available PAM50 subtype scores (n = 2363). Percentages in brackets.
Other members of the SRSF family have also been implicated in breast cancer. For instance, SRSF2 gene amplification at 17q25 has been observed in 6% of breast cancers, although it is uncertain whether this plays an oncogenic role, given evidence that mutations are loss of function in this gene (Chung et al. 2017). Finally, SRSF4 overexpression has been identified in a small subset of breast cancer and its expression has been found responsible for cisplatin-induced AS that leads to apoptosis. Experiments where SRSF4 was silenced showed a decrease in apoptosis upon treatment with cisplatin and highlight the possibility of modulating splicing to regulate chemotherapy sensitivity (Gabriel et al. 2015).
Dysregulation of spliceosomal accessory proteins
Along with the major components of the spliceosome that were described earlier, there are also other regulators of splicing that have been found to be mutated or dysregulated in breast cancer. LIN28A has been identified specifically in HER2-positive breast cancer as being a regulator of AS through interactions with hnRNPA1 (Yang et al. 2015, Xiong et al. 2017). Loss of LIN28A in breast cancer results in isoform switching of the ENAH gene, which is overexpressed in some primary breast tumours (Yang et al. 2015, Xiong et al. 2017). It has also been identified as a feature of the malignant phenotype in a model of breast cancer progression and has been correlated with an unfavourable outcome in HER2-positive breast cancer (Du et al. 2012). Other examples are the epithelial splicing regulatory proteins (ESRP1 and ESRP2), which are splicing factors that have been found to regulate the AS that governs EMT and are amplified in breast cancers (Warzecha et al. 2009, Brown et al. 2011, Bebee et al. 2015) and regulate EMT in breast tumours by activating AKT signalling (Brown et al. 2011). The RNA-binding protein RBFOX2 is also involved in cellular transition, whose upregulation can perturb splicing events in breast cancer (Du et al. 2012, Lapuk et al. 2010). During EMT, RBFOX2-regulated splicing shifts from EMT-specific events, subsequently leading to a higher degree of tissue invasiveness (Braeutigam et al. 2014). Another RNA-binding protein, RBM47, has the ability to alter splicing by binding to introns and 3′ UTRs and loss of expression has been shown to prevent breast cancer progression and metastasis (Vanharanta et al. 2014). Taken together, these lines of evidence point to a fundamental role triggered by splicing dysregulation in breast cancer cells that can cause detrimental effects and lead to the progression of disease.
Evidence of oncogene-induced dependency on the spliceosome
Aside from alterations in spliceosomal component genes themselves, there is emerging evidence that oncogene activation imparts a functional dependency on SF3B1 and other components in breast cancer. A number of spliceosomal component proteins are known transcriptional targets of the oncoprotein MYC (including SF3B1 and SRSF1) and have been shown to both contribute to and cooperate with MYC in malignant transformation (Das et al. 2012, Koh et al. 2015). For instance, MYC addicted TNBCs cells have been shown to impart a specific dependency on the spliceosome via BUD31 and SF3B1 (Hsu et al. 2015) and impaired tumourigenesis was observed when SF3B1 was knocked down or pharmacologically inhibited in breast cancer cells MYC hyperactivation (Hsu et al. 2015). This could be explained due to the increased burden put on the spliceosome when the rate of transcription is increased due to MYC signalling. Recently, knockdown of SF3B1 was found to result in apoptosis in TNBC with MCL-1 inactivation being a likely mechanistic explanation, given MCL-1 is a SF3B1 splicing target (Gao & Koide 2013, Sridhar et al. 2017). Interestingly, MYC and MCL-1 have been shown to cooperate in chemoresistant TNBCs (Lee et al. 2017). This could be further support for the intricate co-operation of MYC with the spliceosome and the resulting changes in isoform dominance that allow the manipulation of cancer cells. In addition, SRSF1 is a known direct target of MYC. MYC induction leads to SRSF1-mediated AS of key protein isoforms involved in proliferation and anchorage-independent growth such as MKNK2 and TEAD1 (Anczukow et al. 2012, Das et al. 2012), which is in part through potentiating eIF4E activation (Anczukow et al. 2012, Das et al. 2012). Together, these studies suggest that multiple spliceosomal proteins are critical MYC targets that contribute to its oncogenic potential by enabling MYC to regulate the expression of specific protein isoforms via AS.
Therapeutic targeting of the spliceosome
There is emerging evidence that disruption of spliceosomal proteins induces selectivity to inhibitors that target the spliceosome. Indeed a number of these inhibitors have been developed including Spliceostatin A, Pladienolides (including E7107) and meyamycin analogues that are all specific SF3B inhibitors as reviewed in (Lee & Abdel-Wahab 2016) that inhibit canonical splicing (Kaida et al. 2007). We, and others, have shown that SF3B1-mutant cells are selectively sensitive to spliceosomal inhibitors (Maguire et al. 2015, Obeng et al. 2016). Moreover, SF3b inhibition in SF3B1 mutant cells resulted in a change in the reversal of the conserved splicing signature, suggesting that SF3B1 mutations are change of function rather than loss of function and that these alterations in aberrant isoforms could be used as biomarkers of therapeutic response (Maguire et al. 2015). There is additional evidence that other spliceosomal gene mutations can be therapeutically targeted with spliceosomal inhibitors. These include SRSF2 mutations, whereby genetically modified mice expressing the Srsf2(P95H) mutation, were sensitive to treatment with the spliceosome inhibitor E7107, which decreased leukaemic burden (Lee et al. 2016). Similar selective sensitivity in mutant U2AF1 cells to sudemycins has also been reported in in vitro and in vivo (Shirai et al. 2015). In addition, MYC-addicted TNBCs have been shown to be more sensitive to inhibition with the spliceosome inhibitor SD6 than MYC non-addicted cells are (Hsu et al. 2015), a mechanism that is likely due to the increased stress and dependency on SF3B1 (as discussed earlier). Further functional studies in the context of clear cell renal carcinoma show that knockdown of SF3B3 in SF3B3-overexpressing cells in vivo reduced tumour growth, highlighting the potential utility of SF3b inhibitors as a therapeutic agent for patients with SF3B3 amplification and/or overexpression (Chen et al. 2017). These lines of evidence raise the possible clinical utility of SF3b inhibitors in patients with additional spliceosomal gene mutations as well as other indirect reliance on the spliceosome. Further studies are warranted to ascertain if overexpression of spliceosomal genes also confers sensitivity to these compounds in breast cancer.
Phase one clinical trials have been performed for E7107 in patients with solid tumours and although the drug has been shown to be on target in patients (i.e. perturbs splicing), the US and European trials were suspended due to an unexpected toxicity involving bilateral optic neuritis (Eskens et al. 2013, Hong et al. 2014). Further studies to understand the causes of toxicity as well as new clinical trials will be necessary to take advantage of splicing’s therapeutic vulnerability in cancer. Currently, H3 biomedicine is testing the compound H3B-8800, which inhibits the SF3b complex and was successful in preclinical studies treating a range of spliceosomal mutant cancers (Buonamici et al. 2016). The compound is now in phase one studies (Nbib2841540) for MDS, acute myeloid leukaemia and chronic myelomonocytic leukaemia.
Conclusions
Mutations and changes in expression of splicing factors that lead to aberrant splicing is a hallmark of cancer that is also relevant to breast cancer. Development of prognostic and predictive aberrant splicing signatures specifically to predict patients that will respond to endocrine (or indeed CDK4/6 inhibitor) therapy could be particularly useful going forward. The increasing technical advances in sequencing methodologies, particularly those that aim to increase RNA read lengths, will undoubtedly enhance the ability to detect these events in the future and further increase our understanding of aberrant transcript expression on breast cancer tumourigenesis and therapy resistance. There is increasing evidence that spliceosomal component genes themselves are dysregulated in breast cancer, through mutations in SF3B1 that are also observed in metastatic disease and upregulation of SF3B3 and SRSF1 in particular, which are associated with resistance to endocrine therapy. Dissecting the function of the expression of the consequent alternatively spliced transcripts would give insight into the mechanism of these alterations and the role they play in therapy resistance. Indeed with the development of spliceosome inhibitors themselves, and exciting preclinical data in other tumour types highlight a potential novel treatment strategy in combination with endocrine therapy and CDK4/6 inhibitors for patients with metastatic disease with spliceosomal gene alterations.
Declaration of interest
The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of this review.
Funding
The authors thank Breast Cancer Now for funding this work as part of Programme Funding to the Breast Cancer Now Toby Robins Research Centre.
References
Alsafadi S, Houy A, Battistella A, Popova T, Wassef M, Henry E, Tirode F, Constantinou A, Piperno-Neumann S, Roman-Roman S, et al. 2016 Cancer-associated SF3B1 mutations affect alternative splicing by promoting alternative branchpoint usage. Nature Communications 7 10615. (https://doi.org/10.1038/ncomms10615)
Anczukow O, Rosenberg AZ, Akerman M, Das S, Zhan L, Karni R, Muthuswamy SK & Krainer AR 2012 The splicing factor SRSF1 regulates apoptosis and proliferation to promote mammary epithelial cell transformation. Nature Structural and Molecular Biology 19 220–228. (https://doi.org/10.1038/nsmb.2207)
Anczukow O, Akerman M, Clery A, Wu J, Shen C, Shirole NH, Raimer A, Sun S, Jensen MA, Hua Y, et al. 2015 SRSF1-regulated alternative splicing in breast cancer. Molecular Cell 60 105–117. (https://doi.org/10.1016/j.molcel.2015.09.005)
Bae J, Leo CP, Hsu SY & Hsueh AJ 2000 MCL-1S, a splicing variant of the antiapoptotic BCL-2 family member MCL-1, encodes a proapoptotic protein possessing only the BH3 domain. Journal of Biological Chemistry 275 25255–25261. (https://doi.org/10.1074/jbc.M909826199)
Bebee TW, Park JW, Sheridan KI, Warzecha CC, Cieply BW, Rohacek AM, Xing Y & Carstens RP 2015 The splicing regulators Esrp1 and Esrp2 direct an epithelial splicing program essential for mammalian development. eLife 15 4. (https://doi.org/10.7554/eLife.08954)
Beije N, Sieuwerts AM, Kraan J, Van NM, Onstenk W, Vitale SR, van der Vlugt-Daane M, Dirix LY, Brouwer A, Hamberg P, et al. 2018 Estrogen receptor mutations and splice variants determined in liquid biopsies from metastatic breast cancer patients. Molecular Oncology 12 48–57. (https://doi.org/10.1002/1878-0261.12147)
Biankin AV, Waddell N, Kassahn KS, Gingras MC, Muthuswamy LB, Johns AL, Miller DK, Wilson PJ, Patch AM, Wu J, et al. 2012 Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes. Nature 491 399–405. (https://doi.org/10.1038/nature11547)
Bingle CD, Craig RW, Swales BM, Singleton V, Zhou P & Whyte MK 2000 Exon skipping in Mcl-1 results in a bcl-2 homology domain 3 only gene product that promotes cell death. Journal of Biological Chemistry 275 22136–22146. (https://doi.org/10.1074/jbc.M909572199)
Boise LH, Gonzalez-Garcia M, Postema CE, Ding L, Lindsten T, Turka LA, Mao X, Nunez G & Thompson CB 1993 bcl-x, a bcl-2-related gene that functions as a dominant regulator of apoptotic cell death. Cell 74 597–608. (https://doi.org/10.1016/0092-8674(93)90508-N)
Bollig A & Miksicek RJ 2000 An estrogen receptor-alpha splicing variant mediates both positive and negative effects on gene transcription. Molecular Endocrinology 14 634–649.
Braeutigam C, Rago L, Rolke A, Waldmeier L, Christofori G & Winter J 2014 The RNA-binding protein Rbfox2: an essential regulator of EMT-driven alternative splicing and a mediator of cellular invasion. Oncogene 33 1082–1092. (https://doi.org/10.1038/onc.2013.50)
Brown RL, Reinke LM, Damerow MS, Perez D, Chodosh LA, Yang J & Cheng C 2011 CD44 splice isoform switching in human and mouse epithelium is essential for epithelial-mesenchymal transition and breast cancer progression. Journal of Clinical Investigation 121 1064–1074. (https://doi.org/10.1172/JCI44540)
Buonamici S, Thomas M, Seiler M, Chan B, Caleb B, Darman R, Fekkes P, Karr C, Keaney C, Klimek V, et al. 2016 H3B-8800, an orally bioavailable modulator of the sf3b complex, shows efficacy in spliceosome-mutan. Nature Medicine 24 497–504. (https://doi.org/10.1038/nm.4493)
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. 2012 The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discovery 2 401–404. (https://doi.org/10.1158/2159-8290.CD-12-0095)
Chen K, Xiao H, Zeng J, Yu G, Zhou H, Huang C, Yao W, Xiao W, Hu J, Guan W, et al. 2017 Alternative splicing of EZH2 pre-mRNA by SF3B3 contributes to the tumorigenic potential of renal cancer. Clinical Cancer Research 23 3428–3441. (https://doi.org/10.1158/1078-0432.CCR-16-2020)
Chung FF, Tan PF, Raja VJ, Tan BS, Lim KH, Kam TS, Hii LW, Tan SH, See SJ, Tan YF, et al. 2017 Jerantinine A induces tumor-specific cell death through modulation of splicing factor 3b subunit 1 (SF3B1). Scientific Reports 7 42504. (https://doi.org/10.1038/srep42504)
Cretu C, Schmitzova J, Ponce-Salvatierra A, Dybkov O, De Laurentiis EI, Sharma K, Will CL, Urlaub H, Luhrmann R & Pena V 2016 Molecular architecture of SF3b and structural consequences of its cancer-related mutations. Molecular Cell 64 307–319. (https://doi.org/10.1016/j.molcel.2016.08.036)
da Luz FAC, Brigido PC, Moraes AS, Araujo RA & Silva MJB 2017 In Splicing Factors in Breast Cancer: Drivers of the Breast Tumor Fate.London, UK: InTech Open. (https://doi.org/10.5772/66162)
Danan-Gotthold M, Golan-Gerstl R, Eisenberg E, Meir K, Karni R & Levanon EY 2015 Identification of recurrent regulated alternative splicing events across human solid tumors. Nucleic Acids Research 43 5130–5144. (https://doi.org/10.1093/nar/gkv210)
Darman RB, Seiler M, Agrawal AA, Lim KH, Peng S, Aird D, Bailey SL, Bhavsar EB, Chan B, Colla S, et al. 2015 Cancer-associated SF3B1 hotspot mutations induce cryptic 3′ splice site selection through use of a different branch point. Cell Reports 13 1033–1045. (https://doi.org/10.1016/j.celrep.2015.09.053)
Das S & Krainer AR 2014 Emerging functions of SRSF1, splicing factor and oncoprotein, in RNA metabolism and cancer. Molecular Cancer Research 12 1195–1204. (https://doi.org/10.1158/1541-7786.MCR-14-0131)
Das S, Anczukow O, Akerman M & Krainer AR 2012 Oncogenic splicing factor SRSF1 is a critical transcriptional target of MYC. Cell Reports 1 110–117. (https://doi.org/10.1016/j.celrep.2011.12.001)
David CJ & Manley JL 2010 Alternative pre-mRNA splicing regulation in cancer: pathways and programs unhinged. Genes and Development 24 2343–2364. (https://doi.org/10.1101/gad.1973010)
Dolatshad H, Pellagatti A, Liberante FG, Llorian M, Repapi E, Steeples V, Roy S, Scifo L, Armstrong RN, Shaw J, et al. 2016 Cryptic splicing events in the iron transporter ABCB7 and other key target genes in SF3B1-mutant myelodysplastic syndromes. Leukemia 30 2322–2331. (https://doi.org/10.1038/leu.2016.149)
Du JW, Xu KY, Fang LY & Qi XL 2012 Clinical significance of Mena and Her-2 expression in breast cancer. European Journal of Gynaecological Oncology 33 455–458.
Dvinge H, Kim E, Abdel-Wahab O & Bradley RK 2016 RNA splicing factors as oncoproteins and tumour suppressors. Nature Reviews Cancer 16 413–430. (https://doi.org/10.1038/nrc.2016.51)
Eskens FA, Ramos FJ, Burger H, O’Brien JP, Piera A, de Jonge MJ, Mizui Y, Wiemer EA, Carreras MJ, Baselga J, et al. 2013 Phase I pharmacokinetic and pharmacodynamic study of the first-in-class spliceosome inhibitor E7107 in patients with advanced solid tumors. Clinical Cancer Research 19 6296–6304. (https://doi.org/10.1158/1078-0432.CCR-13-0485)
Fuqua SA, Fitzgerald SD, Chamness GC, Tandon AK, McDonnell DP, Nawaz Z, O’Malley BW & McGuire WL 1991 Variant human breast tumor estrogen receptor with constitutive transcriptional activity. Cancer Research 51 105–109.
Furney SJ, Pedersen M, Gentien D, Dumont AG, Rapinat A, Desjardins L, Turajlic S, Piperno-Neumann S, de la Grange P, Roman-Roman S, et al. 2013 SF3B1 mutations are associated with alternative splicing in uveal melanoma. Cancer Discovery 3 1122–1129. (https://doi.org/10.1158/2159-8290.CD-13-0330)
Fushimi K, Ray P, Kar A, Wang L, Sutherland LC & Wu JY 2008 Up-regulation of the proapoptotic caspase 2 splicing isoform by a candidate tumor suppressor, RBM5. PNAS 105 15708–15713. (https://doi.org/10.1073/pnas.0805569105)
Gabriel M, Delforge Y, Deward A, Habraken Y, Hennuy B, Piette J, Klinck R, Chabot B, Colige A & Lambert C 2015 Role of the splicing factor SRSF4 in cisplatin-induced modifications of pre-mRNA splicing and apoptosis. BMC Cancer 15 227. (https://doi.org/10.1186/s12885-015-1259-0)
Gao Y & Koide K 2013 Chemical perturbation of Mcl-1 pre-mRNA splicing to induce apoptosis in cancer cells. ACS Chemical Biology 8 895–900. (https://doi.org/10.1021/cb300602j)
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. 2013 Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Science Signaling 6 pl1.
Garcia-Blanco MA, Baraniak AP & Lasda EL 2004 Alternative splicing in disease and therapy. Nature Biotechnology 22 535–546. (https://doi.org/10.1038/nbt964)
Gautrey HL & Tyson-Capper AJ 2012 Regulation of Mcl-1 by SRSF1 and SRSF5 in cancer cells. PLoS ONE 7 e51497. (https://doi.org/10.1371/journal.pone.0051497)
Gokmen-Polar Y, Neelamraju Y, Goswami CP, Gu X, Nallamothu G, Janga SC & Badve S 2015 Expression levels of SF3B3 correlate with prognosis and endocrine resistance in estrogen receptor-positive breast cancer. Modern Pathology 28 677–685. (https://doi.org/10.1038/modpathol.2014.146)
Gracio F, Burford B, Gazinska P, Mera A, Mohd Noor A, Marra P, Gillett C, Grigoriadis A, Pinder S, Tutt A, et al. 2017 Splicing imbalances in basal-like breast cancer underpin perturbation of cell surface and oncogenic pathways and are associated with patients’ survival. Scientific Reports 7 40177. (https://doi.org/10.1038/srep40177)
Herzel L, Ottoz DSM, Alpert T & Neugebauer KM 2017 Splicing and transcription touch base: co-transcriptional spliceosome assembly and function. Nature Reviews Molecular Cell Biology 18 637–650. (https://doi.org/10.1038/nrm.2017.63)
Hong DS, Kurzrock R, Naing A, Wheler JJ, Falchook GS, Schiffman JS, Faulkner N, Pilat MJ, O’Brien J & LoRusso P 2014 A phase I, open-label, single-arm, dose-escalation study of E7107, a precursor messenger ribonucleic acid (pre-mRNA) splicesome inhibitor administered intravenously on days 1 and 8 every 21 days to patients with solid tumors. Investigational New Drugs 32 436–444. (https://doi.org/10.1007/s10637-013-0046-5)
Hsu TY, Simon LM, Neill NJ, Marcotte R, Sayad A, Bland CS, Echeverria GV, Sun T, Kurley SJ, Tyagi S, et al. 2015 The spliceosome is a therapeutic vulnerability in MYC-driven cancer. Nature 525 384–388. (https://doi.org/10.1038/nature14985)
Kaida D, Motoyoshi H, Tashiro E, Nojima T, Hagiwara M, Ishigami K, Watanabe H, Kitahara T, Yoshida T, Nakajima H, et al. 2007 Spliceostatin A targets SF3b and inhibits both splicing and nuclear retention of pre-mRNA. Nature Chemical Biology 3 576–583. (https://doi.org/10.1038/nchembio.2007.18)
Karni R, de Stanchina E, Lowe SW, Sinha R, Mu D & Krainer AR 2007 The gene encoding the splicing factor SF2/ASF is a proto-oncogene. Nature Structural and Molecular Biology 14 185–193. (https://doi.org/10.1038/nsmb1209)
Keren H, Lev-Maor G & Ast G 2010 Alternative splicing and evolution: diversification, exon definition and function. Nature Reviews Genetics 11 345–355. (https://doi.org/10.1038/nrg2776)
Kesarwani AK, Ramirez O, Gupta AK, Yang X, Murthy T, Minella AC & Pillai MM 2017 Cancer-associated SF3B1 mutants recognize otherwise inaccessible cryptic 3′ splice sites within RNA secondary structures. Oncogene 36 1123–1133. (https://doi.org/10.1038/onc.2016.279)
Khan DH, Gonzalez C, Cooper C, Sun JM, Chen HY, Healy S, Xu W, Smith KT, Workman JL, Leygue E, et al. 2014 RNA-dependent dynamic histone acetylation regulates MCL1 alternative splicing. Nucleic Acids Research 42 1656–1670. (https://doi.org/10.1093/nar/gkt1134)
Koh CM, Bezzi M, Low DH, Ang WX, Teo SX, Gay FP, Al-Haddawi M, Tan SY, Osato M, Sabo A, et al. 2015 MYC regulates the core pre-mRNA splicing machinery as an essential step in lymphomagenesis. Nature 523 96–100. (https://doi.org/10.1038/nature14351)
Ladomery M 2013 Aberrant alternative splicing is another hallmark of cancer. International Journal of Cell Biology 2013 463786. (https://doi.org/10.1155/2013/421606)
Lapuk A, Marr H, Jakkula L, Pedro H, Bhattacharya S, Purdom E, Hu Z, Simpson K, Pachter L, Durinck S, et al. 2010 Exon-level microarray analyses identify alternative splicing programs in breast cancer. Molecular Cancer Research 8 961–974. (https://doi.org/10.1158/1541-7786.MCR-09-0528)
Lee SC & Abdel-Wahab O 2016 Therapeutic targeting of splicing in cancer. Nature Medicine 22 976–986. (https://doi.org/10.1038/nm.4165)
Lee SC, Dvinge H, Kim E, Cho H, Micol JB, Chung YR, Durham BH, Yoshimi A, Kim YJ, Thomas M, et al. 2016 Modulation of splicing catalysis for therapeutic targeting of leukemia with mutations in genes encoding spliceosomal proteins. Nature Medicine 22 672–678. (https://doi.org/10.1038/nm.4097)
Lee KM, Giltnane JM, Balko JM, Schwarz LJ, Guerrero-Zotano AL, Hutchinson KE, Nixon MJ, Estrada MV, Sanchez V, Sanders ME, et al. 2017 MYC and MCL1 cooperatively promote chemotherapy-resistant breast cancer stem cells via regulation of mitochondrial oxidative phosphorylation. Cell Metabolism 26 633.e637–647.e637.
Lefebvre C, Bachelot T, Filleron T, Pedrero M, Campone M, Soria JC, Massard C, Levy C, Arnedos M, Lacroix-Triki M, et al. 2016 Mutational profile of metastatic breast cancers: a retrospective analysis. PLoS Medicine 13 e1002201. (https://doi.org/10.1371/journal.pmed.1002201)
Maguire SL, Leonidou A, Wai P, Marchio C, Ng CK, Sapino A, Salomon AV, Reis-Filho JS, Weigelt B & Natrajan RC 2015 SF3B1 mutations constitute a novel therapeutic target in breast cancer. Journal of Pathology 235 571–580. (https://doi.org/10.1002/path.4483)
Marchini C, Gabrielli F, Iezzi M, Zenobi S, Montani M, Pietrella L, Kalogris C, Rossini A, Ciravolo V, Castagnoli L, et al. 2011 The human splice variant Delta16HER2 induces rapid tumor onset in a reporter transgenic mouse. PLoS ONE 6 e18727. (https://doi.org/10.1371/journal.pone.0018727)
Martelotto LG, De Filippo MR, Ng CK, Natrajan R, Fuhrmann L, Cyrta J, Piscuoglio S, Wen HC, Lim RS, Shen R, et al. 2015 Genomic landscape of adenoid cystic carcinoma of the breast. Journal of Pathology 237 179–189. (https://doi.org/10.1002/path.4573)
Martinez-Montiel N, Anaya-Ruiz M, Perez-Santos M & Martinez-Contreras RD 2017 Alternative splicing in breast cancer and the potential development of therapeutic tools. Genes 8 E217. (https://doi.org/10.3390/genes8100217)
Menon R, Im H, Zhang EY, Wu SL, Chen R, Snyder M, Hancock WS & Omenn GS 2014 Distinct splice variants and pathway enrichment in the cell-line models of aggressive human breast cancer subtypes. Journal of Proteome Research 13 212–227. (https://doi.org/10.1021/pr400773v)
Ng KP, Hillmer AM, Chuah CT, Juan WC, Ko TK, Teo AS, Ariyaratne PN, Takahashi N, Sawada K, Fei Y, et al. 2012 A common BIM deletion polymorphism mediates intrinsic resistance and inferior responses to tyrosine kinase inhibitors in cancer. Nature Medicine 18 521–528. (https://doi.org/10.1038/nm.2713)
Obeng EA, Chappell RJ, Seiler M, Chen MC, Campagna DR, Schmidt PJ, Schneider RK, Lord AM, Wang L, Gambe RG, et al. 2016 Physiologic expression of Sf3b1(K700E) causes impaired erythropoiesis, aberrant splicing, and sensitivity to therapeutic spliceosome modulation. Cancer Cell 30 404–417. (https://doi.org/10.1016/j.ccell.2016.08.006)
Pandya-Jones A 2011 Pre-mRNA splicing during transcription in the mammalian system. Wiley Interdisciplinary Reviews: RNA 2 700–717. (https://doi.org/10.1002/wrna.86)
Papaemmanuil E, Cazzola M, Boultwood J, Malcovati L, Vyas P, Bowen D, Pellagatti A, Wainscoat JS, Hellstrom-Lindberg E, Gambacorti-Passerini C, et al. 2011 Somatic SF3B1 mutation in myelodysplasia with ring sideroblasts. New England Journal of Medicine 365 1384–1395. (https://doi.org/10.1056/NEJMoa1103283)
Pereira B, Chin SF, Rueda OM, Vollan HK, Provenzano E, Bardwell HA, Pugh M, Jones L, Russell R, Sammut SJ, et al. 2016 The somatic mutation profiles of 2,433 breast cancers refines their genomic and transcriptomic landscapes. Nature Communications 7 11479. (https://doi.org/10.1038/ncomms11479)
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al. 2000 Molecular portraits of human breast tumours. Nature 406 747–752. (https://doi.org/10.1038/35021093)
Quesada V, Conde L, Villamor N, Ordonez GR, Jares P, Bassaganyas L, Ramsay AJ, Bea S, Pinyol M, Martinez-Trillos A, et al. 2011 Exome sequencing identifies recurrent mutations of the splicing factor SF3B1 gene in chronic lymphocytic leukemia. Nature Genetics 44 47–52. (https://doi.org/10.1038/ng.1032)
Raponi M, Smith LD, Silipo M, Stuani C, Buratti E & Baralle D 2014 BRCA1 exon 11 a model of long exon splicing regulation. RNA Biology 11 351–359. (https://doi.org/10.4161/rna.28458)
Sebestyen E, Zawisza M & Eyras E 2015 Detection of recurrent alternative splicing switches in tumor samples reveals novel signatures of cancer. Nucleic Acids Research 43 1345–1356. (https://doi.org/10.1093/nar/gku1392)
Shapiro IM, Cheng AW, Flytzanis NC, Balsamo M, Condeelis JS, Oktay MH, Burge CB & Gertler FB 2011 An EMT-driven alternative splicing program occurs in human breast cancer and modulates cellular phenotype. PLoS Genetics 7 e1002218. (https://doi.org/10.1371/journal.pgen.1002218)
Shirai C, Tripathi M, Ley J, Ndonwi M, White BS, Tapia R, Saez B, Bertino A, Shao J, Kim S, et al. 2015 Preclinical activity of splicing modulators in U2AF1 mutant MDS/AML. Blood 126 1653
Siegfried Z & Karni R 2017 The role of alternative splicing in cancer drug resistance. Current Opinion in Genetics and Development 48 16–21. (https://doi.org/10.1016/j.gde.2017.10.001)
Silipo M, Gautrey H & Tyson-Capper A 2015 Deregulation of splicing factors and breast cancer development. Journal of Molecular Cell Biology 7 388–401. (https://doi.org/10.1093/jmcb/mjv027)
Sotillo E, Barrett DM, Black KL, Bagashev A, Oldridge D, Wu G, Sussman R, Lanauze C, Ruella M, Gazzara MR, et al. 2015 Convergence of acquired mutations and alternative splicing of CD19 enables resistance to CART-19 immunotherapy. Cancer Discovery 5 1282–1295. (https://doi.org/10.1158/2159-8290.CD-15-1020)
Sridhar P, Chan S, Lock YJ & Petrocca F 2017 Preclinical evaluation of the SF3B1 inhibitor E7107 in triple negative breast cancer. Cancer Research 77 (13 Suppl) abstract 420. (https://doi.org/10.1158/1538-7445.AM2017-420)
Srihari S, Kalimutho M, Lal S, Singla J, Patel D, Simpson PT, Khanna KK & Ragan MA 2016 Understanding the functional impact of copy number alterations in breast cancer using a network modeling approach. Molecular BioSystems 12 963–972. (https://doi.org/10.1039/C5MB00655D)
Stark M, Wichman C, Avivi I & Assaraf YG 2009 Aberrant splicing of folylpolyglutamate synthetase as a novel mechanism of antifolate resistance in leukemia. Blood 113 4362–4369. (https://doi.org/10.1182/blood-2008-08-173799)
Stricker TP, Brown CD, Bandlamudi C, McNerney M, Kittler R, Montoya V, Peterson A, Grossman R & White KP 2017 Robust stratification of breast cancer subtypes using differential patterns of transcript isoform expression. PLoS Genetics 13 e1006589. (https://doi.org/10.1371/journal.pgen.1006589)
Takehara T, Liu X, Fujimoto J, Friedman SL & Takahashi H 2001 Expression and role of Bcl-xL in human hepatocellular carcinomas. Hepatology 34 55–61. (https://doi.org/10.1053/jhep.2001.25387)
Thakur S, Zhang HB, Peng Y, Le H, Carroll B, Ward T, Yao J, Farid LM, Couch FJ, Wilson RB, et al. 1997 Localization of BRCA1 and a splice variant identifies the nuclear localization signal. Molecular and Cellular Biology 17 444–452. (https://doi.org/10.1128/MCB.17.1.444)
Thol F, Kade S, Schlarmann C, Loffeld P, Morgan M, Krauter J, Wlodarski MW, Kolking B, Wichmann M, Gorlich K, et al. 2012 Frequency and prognostic impact of mutations in SRSF2, U2AF1, and ZRSR2 in patients with myelodysplastic syndromes. Blood 119 3578–3584. (https://doi.org/10.1182/blood-2011-12-399337)
Valacca C, Bonomi S, Buratti E, Pedrotti S, Baralle FE, Sette C, Ghigna C & Biamonti G 2010 Sam68 regulates EMT through alternative splicing-activated nonsense-mediated mRNA decay of the SF2/ASF proto-oncogene. Journal of Cell Biology 191 87–99. (https://doi.org/10.1083/jcb.201001073)
Vanharanta S, Marney CB, Shu W, Valiente M, Zou Y, Mele A, Darnell RB & Massague J 2014 Loss of the multifunctional RNA-binding protein RBM47 as a source of selectable metastatic traits in breast cancer. eLife 3. (https://doi.org/10.7554/eLife.02734)
Wang L, Brooks AN, Fan J, Wan Y, Gambe R, Li S, Hergert S, Yin S, Freeman SS, Levin JZ, et al. 2016 Transcriptomic characterization of SF3B1 mutation reveals its pleiotropic effects in chronic lymphocytic leukemia. Cancer Cell 30 750–763. (https://doi.org/10.1016/j.ccell.2016.10.005)
Warzecha CC, Shen S, Xing Y & Carstens RP 2009 The epithelial splicing factors ESRP1 and ESRP2 positively and negatively regulate diverse types of alternative splicing events. RNA Biology 6 546–562. (https://doi.org/10.4161/rna.6.5.9606)
Wei M, Zhu L, Li Y, Chen W, Han B, Wang Z, He J, Yao H, Yang Z, Zhang Q, et al. 2011 Knocking down cyclin D1b inhibits breast cancer cell growth and suppresses tumor development in a breast cancer model. Cancer Science 102 1537–1544. (https://doi.org/10.1111/j.1349-7006.2011.01969.x)
Wendt MK, Taylor MA, Schiemann BJ, Sossey-Alaoui K & Schiemann WP 2014 Fibroblast growth factor receptor splice variants are stable markers of oncogenic transforming growth factor beta1 signaling in metastatic breast cancers. Breast Cancer Research 16 R24. (https://doi.org/10.1186/bcr3717)
Wiener D, Gajardo-Meneses P, Ortega-Hernandez V, Herrera-Cares C, Diaz S, Fernandez W, Cornejo V, Gamboa J, Tapia T, Alvarez C, et al. 2015 BRCA1 and BARD1 colocalize mainly in the cytoplasm of breast cancer tumors, and their isoforms show differential expression. Breast Cancer Research and Treatment 153 669–678. (https://doi.org/10.1007/s10549-015-3575-0)
Xiong H, Zhao W, Wang J, Seifer BJ, Ye C, Chen Y, Jia Y, Chen C, Shen J, Wang L, et al. 2017 Oncogenic mechanisms of Lin28 in breast cancer: new functions and therapeutic opportunities. Oncotarget 8 25721–25735.
Yang J, Bennett BD, Luo S, Inoue K, Grimm SA, Schroth GP, Bushel PR, Kinyamu HK & Archer TK 2015 LIN28A modulates splicing and gene expression programs in breast cancer cells. Molecular and Cellular Biology 35 3225–3243.
Yoshida K & Ogawa S 2014 Splicing factor mutations and cancer. Wiley Interdisciplinary Reviews: RNA 5 445–459. (https://doi.org/10.1002/wrna.1222)
Yoshida K, Sanada M, Shiraishi Y, Nowak D, Nagata Y, Yamamoto R, Sato Y, Sato-Otsubo A, Kon A, Nagasaki M, et al. 2011 Frequent pathway mutations of splicing machinery in myelodysplasia. Nature 478 64–69. (https://doi.org/10.1038/nature10496)
Zhao W, Hoadley KA, Parker JS & Perou CM 2016 Identification of mRNA isoform switching in breast cancer. BMC Genomics 17 181. (https://doi.org/10.1186/s12864-016-2521-9)
Zhou Q, Derti A, Ruddy D, Rakiec D, Kao I, Lira M, Gibaja V, Chan H, Yang Y, Min J, et al. 2015 A chemical genetics approach for the functional assessment of novel cancer genes. Cancer Research 75 1949–1958. (https://doi.org/10.1158/0008-5472.CAN-14-2930)