Mutant U2AF1-expressing cells are sensitive to pharmacological modulation of the spliceosome

Abstract

Somatic mutations in spliceosome genes are detectable in ∼50% of patients with myelodysplastic syndromes (MDS). We hypothesize that cells harbouring spliceosome gene mutations have increased sensitivity to pharmacological perturbation of the spliceosome. We focus on mutant U2AF1 and utilize sudemycin compounds that modulate pre-mRNA splicing. We find that haematopoietic cells expressing mutant U2AF1(S34F), including primary patient cells, have an increased sensitivity to in vitro sudemycin treatment relative to controls. In vivo sudemycin treatment of U2AF1(S34F) transgenic mice alters splicing and reverts haematopoietic progenitor cell expansion induced by mutant U2AF1 expression. The splicing effects of sudemycin and U2AF1(S34F) can be cumulative in cells exposed to both perturbations--drug and mutation--compared with cells exposed to either alone. These cumulative effects may result in downstream phenotypic consequences in sudemycin-treated mutant cells. Taken together, these data suggest a potential for treating haematological cancers harbouring U2AF1 mutations with pre-mRNA splicing modulators like sudemycins.

ARTICLE

Received 25 Mar 2016 | Accepted 24 Nov 2016 | Published 9 Jan 2017

DOI: 10.1038/ncomms14060 OPEN

Mutant U2AF1-expressing cells are sensitive to pharmacological modulation of the spliceosome

Cara Lunn Shirai1,*, Brian S. White1,*, Manorama Tripathi1,*, Roberto Tapia1, James N. Ley1, Matthew Ndonwi1, Sanghyun Kim1, Jin Shao1, Alexa Carver1, Borja Saez2, Robert S. Fulton3, Catrina Fronick3, Michelle OLaughlin3, Chandraiah Lagisetti4, Thomas R. Webb4, Timothy A. Graubert2 & Matthew J. Walter1

Somatic mutations in spliceosome genes are detectable in B50% of patients with myelodysplastic syndromes (MDS). We hypothesize that cells harbouring spliceosome gene mutations have increased sensitivity to pharmacological perturbation of the spliceosome. We focus on mutant U2AF1 and utilize sudemycin compounds that modulate pre-mRNA splicing. We nd that haematopoietic cells expressing mutant U2AF1(S34F), including primary patient cells, have an increased sensitivity to in vitro sudemycin treatment relative to controls. In vivo sudemycin treatment of U2AF1(S34F) transgenic mice alters splicing and reverts haematopoietic progenitor cell expansion induced by mutant U2AF1 expression. The splicing effects of sudemycin and U2AF1(S34F) can be cumulative in cells exposed to both perturbationsdrug and mutationcompared with cells exposed to either alone. These cumulative effects may result in downstream phenotypic consequences in sudemycin-treated mutant cells. Taken together, these data suggest a potential for treating haematological cancers harbouring U2AF1 mutations with pre-mRNA splicing modulators like sudemycins.

1 Division of Oncology, Washington University School of Medicine, St Louis, Missouri 63110, USA. 2 Massachusetts General Hospital Cancer Center, Boston, Massachusetts 02114, USA. 3 McDonnell Genome Institute, Washington University, St Louis, Missouri 63108, USA. 4 SRI International, Bioscience Division, Menlo Park, California 94025, USA. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to M.J.W. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060

Myelodysplastic syndromes (MDS) are the most common adult myeloid malignancy with up to 40,000 new cases diagnosed each year in the United States1,2. MDS are a

heterogeneous group of clonal haematopoietic stem cell disorders characterized by peripheral blood cytopaenias and progenitor expansion; approximately one-third of patients will transform to a secondary acute myeloid leukaemia (AML) that has a poor prognosis3. The only curative therapy is bone marrow transplantation, which is often not an option because of patient comorbidities3. New treatment approaches are greatly needed. At least half of all MDS patient bone marrow samples harbour a mutation in one of several spliceosome genes410, highlighting a potential genetic vulnerability. In addition, spliceosome gene mutations often occur in the founding clones of MDS tumours, providing an attractive target for elimination of all tumour cells10,11. Spliceosome gene mutations are mutually exclusive of each other in patients4,1012, implying either a redundancy in pathogenic function or that a cell cannot tolerate two spliceosome perturbations at once. With this in mind, we hypothesized that cells harbouring a spliceosome gene mutation would have increased sensitivity to further perturbation of the spliceosome by splicing modulator drugs. To examine this, we utilized sudemycin compounds that bind the SF3B1 spliceosome protein and modulate pre-mRNA splicing1315. We used sudemycin D1 and D6, which are synthetic compounds that have been optimized by several rounds of medicinal chemistry for their potent in vivo antitumour activity13. We examined the sensitivity of spliceosome mutant cells to sudemycin treatment, focusing on mutations in the spliceosome gene U2AF1, which have been identied in 11% of MDS patients, utilizing the S34F missense mutation most commonly found in our studies4,5. Mutant U2AF1(S34F) expression has been shown by our group and others to cause altered pre-mRNA splicing in a variety of cell types, as well as altered haematopoiesis and pre-mRNA splicing in mice4,5,1619.

In this manuscript, we provide evidence that U2AF1(S34F)-expressing cells are sensitive to the splicing modulator drug sudemycin. Haematopoietic cells expressing mutant U2AF1 show reduced survival and altered cell cycle in response to sudemycin D6 in vitro. In vivo treatment of U2AF1(S34F) transgenic mice with sudemycin results in an attenuation of mutant U2AF1-induced haematopoietic progenitor cell expansion that is associated with increased cell death. In addition, unsupervised analysis of whole-transcriptome sequencing (RNA-seq) nds that sudemycin D6 perturbs RNA splicing in both mutant U2AF1(S34F)- and U2AF1(WT)-expressing bone marrow cells; however, sudemycin D6 treatment further modulates mutant U2AF1(S34F)-induced splicing changes to create cumulative effects on cells in vivo. The cumulative RNA-splicing effects of sudemycin and mutant U2AF1 may contribute to the downstream phenotypic consequences we observe in vivo.

Results

Sudemycin alters RNA splicing in primary human CD34 cells.

We rst examined the pre-mRNA splicing alterations induced by sudemycin D6 in primary human haematopoietic cells. We treated CD34 haematopoietic progenitor cells isolated

from human umbilical cord blood with 1,000 nM of sudemycin D6 or dimethylsulphoxide (DMSO) vehicle control for 6 h in vitro and performed whole-transcriptome (RNA-seq) analysis (n 6 each, Supplementary Fig. 1). We identied robustly altered

gene expression and pre-mRNA splicing patterns induced by sudemycin, as shown by unsupervised clustering of samples using expressed genes (Fig. 1a) and pre-mRNA splice junctions

(Fig. 1b), respectively. Our analysis identied 1,030 differentially expressed genes (FDRo5%, |log2FC|41) and 18,833 dysregulated splicing events (FDRo5%, |delta per cent spliced in or PSI (DC)|410%, Supplementary Data 1 and 2, respectively) that discriminated between sudemycin D6-treated samples and controls. Sudemycin D6 treatment induced altered pre-mRNA splicing with a bias towards increased exon skipping and intron retention (Fig. 1b). However, there was no apparent bias in the sequence motif surrounding splice acceptor sites of cassette exons that were alternatively spliced (Fig. 1c), in contrast to previously observed biases in sequences surrounding alternatively spliced junctions induced by expression of mutant spliceosome proteins U2AF1, SF3B1 and SRSF2 (refs 1625).

To determine whether particular pathways are enriched for splicing perturbations, we applied GOseq to 6,278 genes with junctions signicantly altered by sudemycin D6 treatment (FDRo5%, |log2FC|42). While pathway enrichment was minimal (enrichment scores o2), GOseq analysis indicated that pathways involved in pre-mRNA splicing, RNA processing and transport, cell cycle, as well as ATPase and helicase activity were enriched in splice junctions altered by sudemycin D6 treatment (FDRo10%; Supplementary Data 3). Genes with sudemycin-altered expression were enriched in pathways involved in receptor and signal transduction activities (FDRo10%; Supplementary

Data 4).

Mutant U2AF1 cells have increased sensitivity to sudemycin. To examine the effects of sudemycin D6 on haematopoietic cells expressing mutant U2AF1, we generated K562 human erythroleukaemia and OCI-AML3 AML cell lines that have stably integrated doxycycline-inducible, FLAG-tagged U2AF1(S34F) or FLAG-tagged U2AF1(WT) to control for U2AF1 overexpression (Supplementary Fig. 2a,b for K562; Fig. 2c,d for OCI-AML3). Mutant U2AF1(S34F)-expressing K562 cells showed reduced survival and lower IC50 (Po0.0001, extra sum-of-squares F-test)

relative to uninduced mutant U2AF1(S34F) and U2AF1(WT)-expressing control cells (Fig. 2a). These effects were also observed in human OCI-AML3 cell lines expressing mutant U2AF1(S34F) compared with U2AF1(WT)-expressing cells and other control cells (Po0.003, extra sum-of-squares F-test; Fig. 2b). Reduced survival of K562 cells in the presence of sudemycin D6 is associated with an altered cell cycle prole: U2AF1(S34F)-expressing K562 cells had a decrease of cells in the S-phase and an increase of cells in the sub-G0/G1 and G2/M phases (Fig. 2c). Furthermore, MDS or AML cells with U2AF1(S34F) mutations treated in vitro with sudemycin D1, a sudemycin compound very similar to D6, showed an increased sensitivity to sudemycin (reduced S-phase) relative to control MDS/AML cells without spliceosome gene mutations (Fig. 2d). In contrast, treatment of MDS/AML patient cells with the chemotherapeutic drug daunorubicin (not predicted to disrupt splicing) showed no specicity for mutant U2AF1(S34F) samples compared with controls (Supplementary Fig. 2e). In addition, human CD34

cells expressing U2AF1(S34F) showed increased sensitivity to another splicing modulator drug (E7107) similar to sudemycin (Supplementary Fig. 2f).

Sudemycin reduces mutant U2AF1 progenitor expansion in vivo. We next examined the effect of sudemycin treatment in vivo on mutant U2AF1(S34F)-induced phenotypes using our previously described U2AF1(S34F) transgenic mouse model19. We induced U2AF1(S34F) or U2AF1(WT) transgenes for 7 days in the bone marrow cells of transplanted mice (to study haematopoietic cell-intrinsic effects) and treated mice concurrently with sudemycin D6 (50 mg kg 1 per day) or vehicle for 5 of

2 NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060 ARTICLE

a b

2

3 Event type

RI

SE

1

2

0

1

1

0

2

1

3

2

Vehicle

n = 9,898 n = 134 n = 90

Sudemycin

Vehicle Sudemycin

c

2

1.5

1

2

1.5

1

2

1.5

1

Information content

0.5

0

Information content

0.5

0

Information content

0.5

0

50

48

46

44

42

40

38

36

34

32

30

28

26

24

22

20

18

16

14

12

10

8

6

4

2

1 3

50

48

46

44

42

40

38

36

34

32

30

28

26

24

22

20

18

16

14

12

10

8

6

4

2

1 3

50

48

46

44

42

40

38

36

34

32

30

28

26

24

22

20

18

16

14

12

10

8

6

4

2

1 3

Position

Position

Position

Figure 1 | Sudemycin D6 alters gene expression and pre-mRNA splicing in primary human CD34 haematopoietic cells. Whole-transcriptome

(that is, RNA-seq) analysis was performed on CD34 cells isolated from human umbilical cord blood following treatment of samples with 1,000 nM

Sudemycin D6 or DMSO vehicle for 6 h (n 6). Unsupervised hierarchical clustering of (a) expressed genes and (b) splice junctions. Skipped exons (SE,

green) and retained introns (RI, purple) event types are visualized. Values are z-scores computed from regularized logarithm values for genes and from per cent spliced in (PSI or C) values for splicing events. (c) Intronic sequence contexts of cassette exon 30 splice sites skipped more often in sudemycin- relative to vehicle-treated cells (FDRo5%, |DC|410%, left panel) or skipped more often in vehicle-relative to sudemycin-treated cells (FDRo5%, |DC|410%, middle panel), along with a context of unperturbed control exons (FDR450%, |DC|o0.1%, right panel). Position is relative to the rst base in the exon.

those days; see schema (Fig. 3a). Sudemycin D6 treatment of transplanted mice showed an attenuation of the previously described19 mutant U2AF1(S34F)-induced haematopoietic progenitor cell expansion by colony-forming unit (CFU-C) assay (Fig. 3b) and by ow cytometry for lineage-, c-Kit ,

Sca1 (KLS) cells (Fig. 3c) when compared with control U2AF1

mutant mice treated with vehicle and mice transplanted with U2AF1(WT)-expressing bone marrow. The attenuation of mutant U2AF1-induced progenitor expansion by sudemycin-treated mice is associated with increased Annexin V staining of

KLS cells (Fig. 3d).

Sudemycin and U2AF1 (S34F) splicing effects can be cumulative. To investigate the potential genotype-specic effects of sudemycin treatment on splice isoform expression, we performed whole-transcriptome sequencing (RNA-seq) on U2AF1(S34F)- and U2AF1(WT)-recipient mouse bulk bone marrow cells following in vivo U2AF1 transgene induction and treatment with sudemycin

D6 (50 mg kg 1 per day for 5 days) or vehicle (Supplementary Fig. 3a,b). RNA was harvested 18 h after the last drug treatment (similar to described above; schema shown in Fig. 3a). Sudemycin D6 treatment at this dose and schedule does not markedly skew the mature lineage distribution within bulk bone marrow of mutant or wild-type (WT) U2AF1 transgenic mice (Supplementary Fig. 3c). Using an unsupervised approach, we observed that sudemycin D6 perturbs splicing in both mutant U2AF1(S34F) and U2AF1(WT)-expressing bone marrow cells (Supplementary Data 59); this is visualized by the segregation of samples according to genotype and treatment within a principal component analysis (PCA) of cassette exon (Fig. 4a) and retained intron (Supplementary Fig. 4a) splicing events. Furthermore, the splicing bias observed in human cells treated with sudemycin (described above) is recapitulated in U2AF1(WT) mouse cells: sudemycin D6 induces exon skipping more often than exon inclusion relative to vehicle (388 of 657 signicant (FDRo10%, |DC|41%) events;

Po2 10 6, one-sided binomial test), as well as intron retention

more often than removal (98 of 145 signicant (FDRo10%,

NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060

a b

WT no dox (IC50=246.5 nM)

MUT no dox (IC50=244.1 nM) MUT+dox (IC50=97.7 nM)

WT+dox (IC50=223.1 nM)

WT no dox (IC50=1214 nM)

MUT no dox (IC50=904.3 nM) MUT+dox (IC50=373.5 nM)

WT+dox (IC50=929.8 nM)

WT no dox

MUT no dox MUT+dox

WT+dox

Surviving fraction

(relative to DMSO)

1.0

0.5

Surviving fraction

(relative to DMSO)

0.0

1.0

0.5

IC50 IC50

0.0

1 2 3 4

1 2 3 4

Log [sudemycin D6 (nM)]

Log [sudemycin D6 (nM)]

c

Sub-G0/G1 phase

S phase

G0/G1 phase

G2/M phase

Cells in G2/M (%)Cells in G0/G1 (%)

Sub-G0/G1 cells (%)

Cells in S (%)

EdU incorporation (%)

4 2 0

8 25

20

15

10

5

0

6

***

** *

***

DMSO 1 100 250

[Sudemycin D6 (nM)]

DMSO 1 100 250

[Sudemycin D6 (nM)]

80 60 40 20

0.0

60

40

20

100 250

[Sudemycin D6 (nM)]

0 DMSO 1 100 250

[Sudemycin D6 (nM)]

0 DMSO 1

d

U2AF1 (S34F)

U2AF1 (WT)

UCB CD34+

1.5

1.0

0.5

0 1 2 3 4 5

Log [sudemycin D1, nM]

Figure 2 | Mutant U2AF1(S34F)-expressing cells display increased sensitivity to sudemycin D in vitro. (a) K562 cells (n 6 for control groups, n 9 for

U2AF1(S34F) treated with doxycycline) or (b) OCI-AML3 cells (n 3 for all groups) with stably integrated, doxycycline-inducible U2AF1(WT) or mutant

U2AF1(S34F) were cultured with increasing concentrations of sudemycin D6 concurrently with doxycycline (250 ng ml 1, where indicated) for 5 days following 2 days of initial induction of mutant or WT U2AF1; total cell numbers were measured. The surviving fraction of cells is shown.

IC50, inhibitory concentration at 50% of maximum cell survival. (c) K562 cell cycle phases were determined using BrdU/7AAD (n 3, representative of

two experiments; *Po0.05, **Po0.01, ***Po0.001 statistics calculated with two-tailed t-tests of MUT dox samples compared with each control group

at a given concentration of sudemycin D6, and the least signicant value is given for each group; mean values with s.d. shown). (d) Primary human MDS or AML cells (both mutant U2AF1(S34F) samples (n 3) and those wild type for U2AF1 (n 6)) or normal umbilical cord blood CD34 cells (n 1) were

cultured on irradiated HS27 stroma, and proliferation (EdU incorporation) was measured after 3 days of exposure to increasing concentrations of sudemycin D1.

|DC|41%) events; Po1.4 10 5, one-sided binomial test).

As in human CD34 cells, the sudemycin-induced

changes were not associated with an apparent sequence motif (Supplementary Fig. 4b); however, we did observe the previously reported increase in a T in the 3 position of the intronic 30 splice

acceptor site of exons more commonly skipped in mutant U2AF1(S34F) cells1619 (Supplementary Fig. 4c). In addition, we dened high-condence sets of U2AF1(S34F) and sudemycin targets, and subsets of those had a high validation rate in orthogonal experimental (NanoString26) and statistical (edgeR27) platforms (Supplementary Information and Supplementary Data 10 and 11).

Next, we examined potential interactions between the drug and mutation within cells, focusing on exon skipping events. Along these lines, we observed that the exon skipping effects

induced by sudemycin D6 (relative to vehicle) within U2AF1(WT)-expressing cells are highly correlated with the drug effects in U2AF1(S34F)-expressing cells (R2 0.8,

Po2.2 10 16, F-test, events signicant in both comparisons

(FDRo10%), Fig. 4b). The vast majority of these events are concordant (in the same direction of induced change with similar magnitude) across genotypes (slope of the regression line 0.75), suggesting that sudemycin treatment results

in similar splicing alterations in these targets in both mutant U2AF1 and WT cells (Fig. 4b). We further assessed drug/genotype interaction using a statistical linear model: of 32,529 dysregulated splicing events (across all event types), only 136 showed statistically signicant evidence of interaction (that is, synergy or antagonism; DEXSeq, FDRo10%). However, when sudemycin D6 and mutant U2AF1(S34F) dysregulate a

4 NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060 ARTICLE

a

IV catheter placement (IJ), begin transgene induction

Begin drug treatment

Euthanized and analysis

BM txp

100

Number of colonies/

10,000 cells

KLS cells annexin V+ (%)

6 wks Day 0 Day 1 Day 2 Day 3 Day 4 Day 5 Day 6 Day 7 Day 8

Doxycycline chow

Sudemycin D6 (or vehicle)4 hr. IV infusion daily x 5 days

b c

p<0.01 p<0.001

p=0.07

p<0.02

p=0.07

0.8

80

60

40

20

0.6

Percent of donor

bone marrow cells

0.4

0.2

0.0

0 Vehicle

Sudemycin D6

Vehicle

Sudemycin D6

U2AF1(WT)/rtTA

U2AF1(S34F)/rtTA

15

10

5

d

p=0.02

U2AF1(WT)/rtTA

U2AF1(S34F)/rtTA

0 Vehicle

Sudemycin D6

Vehicle

Sudemycin D6

Figure 3 | Sudemycin D6 treatment attenuates mutant U2AF1-induced progenitor cell expansion in U2AF1(S34F) transgenic mice. (a) Schema of sudemycin treatment of transgenic mice in vivo. Doxycycline-inducible U2AF1(S34F) or U2AF1(WT) transgenic mouse bone marrow was transplanted into recipient mice. Intrajugular (IJ) catheters were placed for IV drug infusions, which were performed over 4 h for 5 days. (b) Haematopoietic progenitor CFU-C colony-forming assay and (c) ow cytometry for haematopoietic stem and progenitor (HSPC) cell surface markers (c-Kit , lineage , Sca1 , KLS,

right panel) on U2AF1(WT)- or (S34F) mutant-recipient mouse bone marrow following treatment is shown (n 6 for both U2AF1(WT) conditions, n 7

for vehicle-treated U2AF1(S34F), n 11 for sudemycin D6-treated U2AF1(S34F)). (d) Annexin V KLS cells were quantied in the bone marrow of mice

following in vivo sudemycin treatment as described above (n 79). Data pooled from two independent experiments (bd); statistics calculated using

two-tailed t-tests for each comparison shown; mean values with s.d. shown.

junction in the same direction (for example, both increasing exon skipping), this results in a cumulative effect in a sudemycin-treated mutant cell that is greater than the effect induced by sudemycin treatment of WT cells (indicated by red and blue colour in Fig. 4bd). As an example, U2AF1(S34F) expression induces increased exon skipping in a 4932438A13Rik splice junction compared with U2AF1(WT) in vehicle-treated cells (DC 0.205, middle two columns in Fig. 4c). Sudemycin

D6 also increases exon skipping of the same junction in U2AF1(WT) (DC 0.188) and U2AF1(S34F) (DC 0.204) cells

(that is, sudemycin-induced exon skipping is independent of

genotype; Fig. 4c). Therefore, as expected, the cumulative effect induced by both mutation and drug treatment (DC 0.409) exceeds the individual effects of mutant U2AF1

expression or sudemycin treatment alone, ultimately resulting in different levels of exon skipping in sudemycin-treated cells expressing mutant U2AF1 versus WT cells (Fig. 4c). The cumulative effect of mutation and drug can also be observed in splicing events that result in increased exon inclusion (Fig. 4d). Together, these data signify that the effects of sudemycin and U2AF1(S34F) on splicing in the same cell can be cumulative that is, greater than the separate effects of drug or mutant.

NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060

a b

0.6

0.0

20

PC2 (13.77% variance)

Condition

MT Sud

WT Sud

MT Veh

WT Veh

(MT Sud) - (MT Veh)

0.3

(MT Sud) -(WT Sud)

10

10

0.100.05

0.05 0.10

0.00

0

0.3

0.6

0.6

0.3

0.0

0.3

0.6

20

(WT Sud) - (WT Veh)

30

20

10 10 20

0

PC1 (18.39% variance)

c d

0.0

relative to WT Veh

0.4

0.3

0.2

0.1

0.0

relative to WT Veh

0.1

0.2

0.3

WT Sud

WT Veh

MTVeh WT

Sud

MT Sud

WT Veh

MT Veh

MT Sud

Condition

Condition

e f

log2 fold change (MT SUD vs WT SUD)

40

1

PC3 (10.91% variance)

log2 fold change (MT Sud vs MT Veh)

Condition

0

0

MT Sud

WT Sud

MT Veh

WT Veh

1

0

1

40

log2 fold change (WT Sud vs WT Veh)

1.00.5

0.5 1.0

0.0

100 50 50

1

0

PC1 (30.64% of variance)

Figure 4 | Sudemycin D6 treatment alters splicing and gene expression in mutant U2AF1-haematopoietic cells. RNA sequencing and transcriptome analysis was performed on RNA harvested from mouse bone marrow cells expressing mutant U2AF1(S34F) or U2AF1(WT) following treatment of mice with sudemycin D6 (50 mg kg 1) or vehicle control for 5 days (n 5 per genotype and treatment). (a) PCA of normalized expression of skipped

cassette exon events. (b) Correlation between sudemycin-induced changes relative to vehicle treatment in cassette exons of U2AF1(WT) cells DCWTSud CWT;Sud CWT;Veh

(FDRo10%) and of mutant U2AF1(S34F) cells DCMTSud CMT;Sud CMT;Veh

(FDRo10%; R2 0.8; Po2.2 10 16, F-test).

Dashed line indicates similar effects in U2AF1(WT) and U2AF1(S34F) cells DCWTSudDCMTSud

. Colour scale indicates cumulative splicing changes of U2AF1(S34F) expression with sudemycin treatment (that is, DCSudMT CMT;Sud CWT;Sud), with red being a positive change, blue being a negative change, and white being no difference between sudemycin-treated U2AF1(S34F) and U2AF1(WT) cells. (c,d) Delta PSI (DC) for each condition on the horizontal axis relative to vehicle-treated U2AF1(WT) cells ((CWT;Veh) 0) for events that are signicantly dysregulated across all ve comparisons (see

Supplementary Methods), have concordant sudemycin-induced dysregulation across genotype (that is, DCWTSud and DCMTSud have the same direction), and in which the sudemycin effect is exacerbated by mutation. The cumulative effect on cassette exons induced by both mutant U2AF1 and sudemycin D6 treatment relative to vehicle-treated WT cells, DCcum, may result in increased exon skipping (positive values, c) or increased exon inclusion (negative values, d). (e) PCA of normalized expression of expressed genes. (f) Correlation between sudemycin-induced changes relative to vehicle treatment (expressed as log2 fold changes, FDRo10%) in U2AF1(WT) and U2AF1(S34F) cells (R2 0.8; Po2.2 10 16, F-test). Dashed line indicates log2 fold

changes induced by sudemycin are the same in U2AF1(WT) and U2AF1(S34F) cells. Colour scale indicates cumulative contribution of mutant U2AF1 expression to sudemycin-induced gene expression changes, that is, log2 fold change of gene expression altered by sudemycin in mutant U2AF1(S34F) cells relative to U2AF1(WT) cells, with red being a positive change, blue being a negative change and white being no difference. U2AF1 mutant (MT), sudemycin D6 (Sud), U2AF1 WT, vehicle (Veh), principal component 1 (PC1), principal component 2 (PC2), principal component 3 (PC3), per cent spliced in (PSI, C).

6 NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060 ARTICLE

Sudemycin induces gene expression alterations. Sudemycin treatment in vivo also results in altered gene expression in mutant U2AF1 mouse bone marrow cells compared with U2AF1(WT)-expressing control cells (Supplementary Data 12). As with splicing alterations, mouse bone marrow samples segregate by genotype and by treatment in an unsupervised analysis of gene expression (Fig. 4e). As seen at the junction level, there is a high-degree of correlation between sudemycin-induced, gene-level effects (direction and magnitude relative to vehicle) across mutant and WT genotypes (R2 0.8,

Po2.2 10 6, F-test, events signicant (FDRo10%) in both

comparisons, Fig. 4f). We also observed cumulative effects of mutation and drug on gene expression (colour in Fig. 4f), which may result in downstream cellular pathway changes. Using GOseq, we identied pathways that were most enriched for differentially expressed genes in mutant cells treated with sudemycin D6 relative to the other genotype and treatment groups (FDRo10%). We found mutant cells treated with sudemycin D6 were enriched in biologic pathways related to immune and inammatory responses, antigen processing and presentation, cytokine production, leukocyte differentiation, cell death and apoptotic processes, when compared with the other genotype and treatment groups (Supplementary Data 13).

DiscussionWe provide evidence suggesting that U2AF1(S34F)-expressing cells are sensitive to the splicing modulator drug sudemycin. Haematopoietic cells expressing mutant U2AF1 have reduced survival and cell cycle changes following sudemycin D6 treatment in vitro. In vivo, sudemycin D6 is capable of attenuating mutant U2AF1-associated expansion of haematopoietic progenitors in transgenic mice that is associated with increased cell death. In addition, while the effects of sudemycin D6 treatment on splicing and gene expression can be independent of the effects of mutant U2AF1 expression, sudemycin D6 treatment can also modulate splicing changes induced by mutant U2AF1(S34F) to create cumulative effects on cells in vivo. Together, these data indicate that mutant U2AF1(S34F)-expressing cells may have a therapeutic vulnerability to splicing modulator drugs such as sudemycin.

Previous studies using splicing modulator drugs, as well as non-pharmacological methods to target splicing factors, have indicated that the spliceosome is a promising therapeutic target for many cancer cell types13,14,2831. Recent studies have suggested that cancers with spliceosome gene mutations may have an increased sensitivity to splicing modulator drugs like sudemycin;3234 this study provides further evidence to support this hypothesis. Specically, we show that mutant U2AF1 cells are sensitive to in vivo treatment with sudemycin. Using unbiased RNA-seq, we show that sudemycin D6 has effects on mutant U2AF1-associated splicing changes, resulting in different levels of transcript isoforms in cells expressing mutant U2AF1 compared with WT following sudemycin treatment. Ultimately, cumulative changes in isoform expression could be the cause of the mutant U2AF1-specic responses to sudemycin treatment that we observe in vitro and in vivo. Identication of critical downstream targets of sudemycin in vivo will likely require examination of various cellular populations within the bone marrow (including haematopoietic stem and progenitor cells) and harvesting cells at more immediate time points following drug treatmentboth limitations of the current study. Future studies will focus on generating mutant U2AF1-expressing leukaemias in mice to test the efcacy of sudemycin on fully transformed haematopoietic tumours, as has

recently been reported using E7107 to treat Srsf2(P95H)-expressing leukaemias33.

Our data and others highlight several possible mechanisms for the sensitivity of spliceosome mutant samples to splicing modulator therapy. Pathway analysis of differentially expressed genes induced by sudemycin revealed enrichment in inamma-tory signalling pathways. This is consistent with the enrichment in biological pathways related to cytokine and immune signalling observed following in vivo treatment of Srsf2(P95H) mutant mice with the splicing modulator E7107 (ref. 33), raising the possibility that spliceosome mutant-specic phenotypes observed in mice following splicing modulator drug treatment may be driven by an altered inammatory response in mutant cells. Whether the altered inammatory response in mutant cells treated with splicing modulator drugs is a direct result of mutant-altered splicing or gene expression could be explored in future studies. Alternatively, it is also possible that the cumulative effect of sudemycin and mutant U2AF1 expression on pre-mRNA splicing may simply create a state of spliceosome sickness in cells by exceeding a tolerable threshold of splicing perturbations and their downstream consequences. Along these lines, other splicing modulator drugs have been shown to cause increased intron retention, R loop formation, DNA damage and cell death28,29,31,35. Spliceosome gene mutations and splicing modulator drugs may both induce these consequences in a cell, and the cumulative effect of the drug and mutation may create a toxic intracellular milieu.

Exploring the clinical utility of splicing modulator therapies in MDS patients with spliceosome mutations who have failed current therapies is warranted and currently being pursued (NCT02841540), as these patients have few treatment options. Whether excessive toxicity will occur in WT cells treated with newer splicing modulators is a major question. Initial phase I clinical trial studies using a splicing modulator, E7107, showed toxicities in some patients with solid tumours, including ocular toxicity (NCT00459823, NCT00499499). The mechanism for the toxicity is not known. Moving forward, mutant U2AF1 and other preclinical models of spliceosome mutations will be valuable in further testing the in vivo efcacy and toxicity of drugs that modulate splicing through various mechanisms. Collectively, our data suggest that a mutant splicing factor together with a splicing modulator drug like sudemycin may create a unique cellular toxicity that could be exploited for therapeutic purposes.

Methods

Isolation and drug treatment of CD34 haematopoietic cells. Mononuclear

cells from human umbilical cord blood were separated by Ficoll gradient centrifugation (each n value is a pooled set of four individual umbilical cord blood samples). CD34 cells were enriched using autoMACS-positive

selection (CD34 MicroBead Kit, Miltenyi Biotec) according to the manufacturers instructions to achieve 490% purity of CD34 cells.

Following isolation, CD34 cells were cultured in X-Vivo media

(Lonza) with cytokines (SCF, FLT3L, IL3 and TPO) overnight. For whole-transcriptome sequencing with sudemycin, either sudemycin D6 (1,000 nM) or the vehicle DMSO was added to cells for 6 h, and cells were then harvested for RNA. Tissue acquisition was performed per protocol approved by the Washington University School of Medicine Institutional Review Board.

RNA-seq of CD34 cells treated with sudemycin. RNA was isolated

from human CD34 cells using the miRNEasy Kit (Qiagen), and removal

of genomic DNA was performed via a Turbo DNA-Free Kit (Ambion). Ribosomal RNA was removed using Ribozero (Epicenter), followed by cDNA preparation and generation of stranded libraries using the TruSeq Stranded Total RNA Sample Prep Kit (Illumina). Sequencing was performed on the HiSeq2500 platform (Illumina) to generate 2 125 bp paired-end

reads. RNA-seq data were deposited in NCBI dbGAP (phs000159.v9).

NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060

RNA-seq analysis of human CD34 cells treated with sudemycin. Analysis of

(stranded) RNA-seq data generated from primary human CD34 haematopoietic

cells treated with sudemycin D6 was performed using the genome modelling system36. Reads were aligned to the human genome (hg19/NCBI build 37)using TopHat37 version 2.0.8, with annotations provided by Ensembl38version 67. All downstream bioinformatic and statistical analyses, including calculation of P values using Fishers exact test and simulation, were performed in R39 and python. The heatmap of gene expression was created using normalized expression values output by DESeq2 package (version 1.6.3 (ref. 40)) that was subsequently z-scored. DESeq2 was used to detect differentially expressedgenes, and gene-set enrichment analysis was then performed by applying GOseq41 (version 1.18.0) to differentially expressed genes (FDRo5%, |log2FC|41,

DESeq2). Additional details of RNA-seq analysis of gene expression are further described in the Supplementary Methods.

Alternative splicing events were computed via rMATS42 (version 3.0.8).

We generated the heatmap of splice junctions using per cent spliced in(PSI or C) values for exon skipping and intron retention events. For downstream analyses, namely reported dysregulated events and visualization of sequence contexts below, P values for events passing the lter were re-adjusted for multiple hypothesis testing using the method of Benjamini and Hochberg43.

A positive DPSI value indicates an event that was spliced in more often in the vehicle-treated relative to the sudemycin-treated samples. Sequence contexts of splice sites in exon skipping and intron retention events were visualized using seqLogo version 1.32.1. Dysregulated events visualized were those with |DPSI|410% and having a post-ltering re-adjusted FDRo5% from rMATS.

Dysregulated splicing junctions were determined by DEXSeq4 version 1.12.2. Only expressed junctions were analysed, and gene-set enrichment analysis was performed by applying GOseq to dysregulated junctions (FDRo5%, |log2FC|41,

DEXSeq). All of these are described in more detail in the Supplementary Methods.

Creation of inducible U2AF1(S34F) and U2AF1(WT) cell lines. Doxycyclineinducible FLAG-tagged WT U2AF1 or FLAG-tagged mutant U2AF1(S34F) lentiviral expression plasmids were previously described4, and lentivirus was generated in 293T cells with the packaging plasmids pMD-G, pMD-Lg and REV. Concentrated virus (multiplicity of infection (MOI) of 3 for K562, MOI of 5 for OCI-AML3) was used to transduce K562 cells (ATCC, CCL243) or OCI-AML3 (DSMZ, ACC 582). Transduced cells (marked by green uorescent protein) were isolated by ow cytometry cell sorting. Expression of mutant or WT U2AF1 was induced using the indicated concentrations of doxycycline hyclate (Sigma, St Louis, MO) in water.

Cell line counting and BrdU incorporation assays. K562 (or OCI-AML3) cells in culture were seeded at B200,000 cells per ml with different concentrations of Sudemycin D6 (or with DMSO control) following 48 h of initial doxycycline treatment to induce mutant U2AF1(S34F) or U2AF1(WT) expression. Cells were counted on Days 0, 3 and 5 of drug treatment using ow cytometry counting particles (Spherotech) as per the manufacturers recommendations, along with propidium iodide (Millipore) to exclude dead cells from counts. Day 5 data were graphed in Graphpad Prism (Graphpad Software Inc.) using a nonlinear regression best t line of the log (drug) versus response for each genotype and treatment. Statistical differences between curves were examined by the extra sum-of-squares F-test function in Graphpad Prism. For K562 cells, 5-bromodeoxyuridine (BrdU) incorporation was performed using the BrdU Flow Kit (APC, #559619, BD Biosciences) as per the manufacturers instructions using a 45 min pulse of BrdU of cells in culture on Day 5 of drug treatment.

Culture of MDS/AML cells and EdU incorporation assay. Primary MDS and AML cells were cultured on irradiated (4,000 cGy) HS27 stroma as described previously45. Briey, MDS or AML cells were cultured on stroma for 23 days and then treated with increasing doses of sudemycin D1 (or daunorubicin in Supplementary Information). MDS and AML cells were cultured with drugfor 3 days and analysed for cell proliferation by EdU incorporation assay (Invitrogen) via ow cytometry. All experiments were performed on 96-well plates. Studies were performed per protocol approved by the Washington University School of Medicine Institutional Review Board and with patients consent for sample use.

Murine bone marrow transplant and drug infusion. To generate mice for each experiment, 1 106 transgenic mouse donor bone marrow cells from two

to three mice pooled (CD45.2) were transplanted into at least ve lethally irradiated (1,100 rads) congenic WT recipient mice (C57BL/6 129S4Sv/Jae)F1

(CD45.1/CD45.2) per genotype, as previously described19. Donor mice were between 8 and 12 weeks of age, and recipient mice ranged from 6 to 12 weeksof age; donor and recipient mice were sex-matched (both sexes were used in experiments). Donor chimerism was conrmed Z6 weeks post transplantation to ensure engraftment of transgenic bone marrow. Post-engraftment, intrajugular catheters were surgically placed for intravenous drug infusions of mice over4 h daily for 5 days with either sudemycin D6 (50 mg kg 1) or vehicle control (HP-b-CD (2-hydroxypropyl)-b-cyclodextrin in phosphate buffer pH 7.4). This

dose was determined by prior studies13. Mice received doxycycline chow to induce U2AF1(S34F) or U2AF1(WT) transgene on the day of surgery, and drug infusion began 2 days later. Mice of each group were randomly selected for treatment with sudemycin or vehicle control. Investigators were not blinded to the group allocation during the experiment or analysis. Animals were excluded from analysis if their catheters did not remain patent during the entire treatment period. Mice were euthanized for analysis the day following the last drug infusion. Sample sizes for experiments were chosen to allow for statistical comparison between groups. All mouse procedures were performed according to the protocols approved by the Washington University Animal Studies Committee.

Mouse haematopoietic progenitor cell assay. Methylcellulose progenitor CFU-C assays were performed using Methocult GF M3434 (Stem Cell Technologies). Bulk bone marrow cells were obtained from experimental mice, and red blood cells were lysed before plating of 10,000 bone marrow cells per 1.3 ml media; each sample was evaluated in duplicate. Progenitor colonies (dened as

Z40 cells per colony) were counted following 7 days culture at 37 C with 5% CO2.

Mouse haematopoietic cell ow cytometry. For ow cytometry, all antibodies are from eBioscience (unless indicated), and catalogue number provided (if available). For haematopoietic progenitor/stem cells, we used the following antibodies (volume of antibody used in 200 ml of uorescence-activated cell sorting (FACS) buffer for staining is also indicated): CD45.1-APC (#17-0453, 3 ml), CD45.2-PE (#12-0454, 3 ml), Biotin-conjugated lineage (Gr-1 (#13-5931,0.25 ml), Cd3e (#13-0032, 0.5 ml), B220 (#13-0452, 0.5 ml), Ter119 (#13-5921, 0.5 ml) and CD41(#13-0411, 1 ml)), streptavidin secondary-eFluor605NC, c-Kit-APCeFluor780 (#47-1172, 1.5 ml), Sca1-PerCP-Cy5.5 (#45-5981, 0.5 ml). Flow cytometry for Annexin V staining of KLS cells was performed

following initial staining of KLS cells as described above using Annexin V-APC (2.5 ml, BD Biosciences, #550474) incubated in 1 Annexin V binding

buffer (BD Biosciences). All other incubations occurred in FACS buffer. All ow cytometry analyses were performed using FACScan or Gallios cytometers (BD Biosciences) and analysed using the FlowJo software (FlowJo, LLC, Ashland, OR, USA).

RNA-seq of mouse bone marrow cells in vivo. RNA was isolated and prepared as described above for human cells, using the miRNEasy Kit (Qiagen) followed by removal of genomic DNA via a Turbo DNA-Free Kit (Ambion) and ribosomal RNA using Ribozero (Epicenter). The cDNA preparation and generation of stranded libraries was performed using the TruSeq Stranded Total RNASample Prep Kit (Illumina). Sequencing was performed on the HiSeq2500 platform (Illumina) to generate 2 126 bp paired-end reads. The RNA-seq

data have been deposited in NCBIs Gene Expression Omnibus46 and are accessible through the GEO Series accession number GSE89834.

RNA-seq analysis of mouse bone marrow cells in vivo. Analysis of (stranded) RNA-seq data generated from U2AF1(S34F) and U2AF1(WT) murine cells treated with sudemycin D6 or vehicle was performed similarly to the human CD34 cell analysis described above, with some differences as follows.

Alignment again utilized genome model system, and reads were aligned to the mouse genome (mm9/NCBI build 37). Per cent spliced in (PSI or C) values were calculated for all four conditions ({U2AF1(S34F), U2AF1(WT)}

{sudemycin D6, vehicle}) using rMATS. PCA was performed independently on events annotated by rMATS as skipping cassette exons or retaining introns using the z-scored C values of these events.

To quantitate the simultaneous effect of treatment and genotype for each splicing event, we calculated the change in per cent spliced in values for thefour pairwise comparisons in which the genotype (alternately, treatment) was the same in the pair, but in which the treatment (alternately, genotype) differed.

We refer to the unchanged condition as the context and to the two conditions that differ as A and B and denote the corresponding change in the per cent spliced in values as DCcontextA B C

context;A

context;B. In addition, we dened the cumulative effect in a mutant, drug-treated cell relative to a WT, vehicle-treated cell as DCcum CSud;MT CVeh;WT. Each of the ve comparisons described

were evaluated using rMATS, and P values for these events were then adjusted for multiple hypothesis testing as described above. Scatterplots were then plotted of cassette exon-skipping events dysregulated by sudemycin in both U2AF1(WT) and U2AF1(S34F) contexts ( DCWTSudVeh

41%, DCMTSudVeh

41%, FDRo10%

in both comparisons). To highlight cumulative cassette exon-skipping effects, trajectories comparing C values to the baseline CVeh;WT were plotted. Events

that were signicant in all ve comparisons (that is, DCWTSudVeh

41%,

DCMTSud Veh

41%, DCVehMT WT

41%, DCSudMT WT

C

41%, DCCum

j j41%,

with FDRo10% in all comparisons) were plotted; the latter condition ensures that visually discernible differences are statically signicant and not attributable to statistical noise. Additional details of this can be found in the Supplementary Methods.

8 NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060 ARTICLE

Non-additive (that is, super-additive synergistic or sub-additive antagonistic) interactions between drug and mutation were assessed in DEXSeq by comparing a generalized linear model that included main effects for drug and mutation with a second model that additionally included an interaction term representing drug/mutation synergy or antagonism.

As above, splice site sequence contexts of exon-skipping events dysregulated by sudemycin D6 ( DCWTSud Veh

41%, FDRo10%) and corresponding

unperturbed control events ( DCWTSud Veh

40.1%, FDR450%) were plotted

using seqLogo. Similarly, sequence contexts were displayed for events dysregulated ( DCVehMT WT

41%, FDRo10%) or unaffected ( DCVehMT WT

41%,

FDR450%) by U2AF1(S34F).

At the gene-level, PCA was performed as described above for human expressed genes. Genes differentially expressed within the above-described ve conditions were determined via DESeq2. Scatterplots were then made comparing log2 fold changes of genes dysregulated by sudemycin in both U2AF1(WT) and

U2AF1(S34F) contexts (FDRo10% in both comparisons). Pathwaysenriched (FDRo10%) for genes differentially expressed between U2AF1(S34F), sudemycin-treated cells and U2AF1(WT) and/or vehicle-treated cells were independently determined in a pairwise manner using GOseq. Again, expanded details of this approach can be found in the Supplementary Methods.

Data availability. All relevant data generated in this study are available at data-deposition sites. For human CD34 cells treated with sudemycin D6 in vitro,

data are available at NCBI dbGAP (phs000159.v9). For transgenic mice expressing mutant U2AF1(S34F) or U2AF1(WT) and treated with sudemycin D6 in vivo, data are available at Gene Expression Omnibus (GSE89834).

References

1. Ma, X. Epidemiology of myelodysplastic syndromes. Am. J. Med. 125, S2S5 (2012).

2. Cogle, C. R., Craig, B. M., Rollison, D. E. & List, A. F. Incidence of the myelodysplastic syndromes using a novel claims-based algorithm: high number of uncaptured cases by cancer registries. Blood 117, 71217125 (2011).

3. Troy, J. D., Atallah, E., Geyer, J. T. & Saber, W. Myelodysplastic syndromes in the United States: an update for clinicians. Ann. Med. 46, 283289 (2014).

4. Yoshida, K. et al. Frequent pathway mutations of splicing machinery in myelodysplasia. Nature 478, 6469 (2011).

5. Graubert, T. A. et al. Recurrent mutations in the U2AF1 splicing factor in myelodysplastic syndromes. Nat. Genet. 44, 5357 (2011).

6. Papaemmanuil, E. et al. Somatic SF3B1 mutation in myelodysplasia with ring sideroblasts. N Engl. J. Med. 365, 13841395 (2011).

7. Visconte, V. et al. SF3B1, a splicing factor is frequently mutated in refractory anemia with ring sideroblasts. Leukemia 26, 542545 (2011).

8. Damm, F. et al. Mutations affecting mRNA splicing dene distinct clinical phenotypes and correlate with patient outcome in myelodysplastic syndromes. Blood 119, 32113218 (2012).

9. Thol, F. et al. Frequency and prognostic impact of mutations in SRSF2, U2AF1, and ZRSR2 in patients with myelodysplastic syndromes. Blood 119, 35783584 (2012).

10. Walter, M. J. et al. Clonal diversity of recurrently mutated genes in myelodysplastic syndromes. Leukemia 27, 12751282 (2013).

11. Papaemmanuil, E. et al. Clinical and biological implications of driver mutations in myelodysplastic syndromes. Blood 122, 36163627 quiz 3699 (2013).

12. Haferlach, T. et al. Landscape of genetic lesions in 944 patients with myelodysplastic syndromes. Leukemia 28, 241247 (2014).

13. Lagisetti, C. et al. Optimization of antitumor modulators of pre-mRNA splicing. J. Med. Chem. 56, 1003310044 (2013).

14. Fan, L., Lagisetti, C., Edwards, C. C., Webb, T. R. & Potter, P. M. Sudemycins, novel small molecule analogues of FR901464, induce alternative gene splicing. ACS Chem. Biol. 6, 582589 (2011).

15. Convertini, P. et al. Sudemycin E inuences alternative splicing and changes chromatin modications. Nucleic Acids Res. 42, 49474961 (2014).16. Okeyo-Owuor, T. et al. U2AF1 mutations alter sequence specicity of pre-mRNA binding and splicing. Leukemia 29, 909917 (2014).

17. Brooks, A. N. et al. A pan-cancer analysis of transcriptome changes associated with somatic mutations in U2AF1 reveals commonly altered splicing events. PLoS ONE 9, e87361 (2014).

18. Ilagan, J. O. et al. U2AF1 mutations alter splice site recognition in hematological malignancies. Genome Res. 25, 1426 (2014).

19. Shirai, C. L. et al. Mutant U2AF1 expression alters hematopoiesis and pre-mRNA splicing in vivo. Cancer Cell 27, 631643 (2015).

20. Przychodzen, B. et al. Patterns of missplicing due to somatic U2AF1 mutations in myeloid neoplasms. Blood 122, 9991006.

21. DeBoever, C. et al. Transcriptome sequencing reveals potential mechanism of cryptic 30 splice site selection in SF3B1-mutated cancers. PLoS Comput. Biol.

11, e1004105 (2015).

22. Darman, R. B. et al. Cancer-associated SF3B1 hotspot mutations induce cryptic 30 splice site selection through use of a different branch point. Cell Rep. 13, 10331045 (2015).

23. Alsafadi, S. et al. Cancer-associated SF3B1 mutations affect alternative splicing by promoting alternative branchpoint usage. Nat. Commun. 7, 10615 (2016).

24. Kim, E. et al. SRSF2 mutations contribute to myelodysplasia by mutant-specic effects on exon recognition. Cancer Cell 27, 617630 (2015).

25. Zhang, J. et al. Disease-associated mutation in SRSF2 misregulates splicing by altering RNA-binding afnities. Proc. Natl Acad. Sci. USA 112, E4726E4734 (2015).

26. Geiss, G. K. et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat. Biotechnol. 26, 317325 (2008).

27. Robinson, M. D. & Smyth, G. K. Moderated statistical tests for assessing differences in tag abundance. Bioinformatics 23, 28812887 (2007).

28. Kaida, D. et al. Spliceostatin A targets SF3b and inhibits bothsplicing and nuclear retention of pre-mRNA. Nat. Chem. Biol. 3, 576583 (2007).

29. Kotake, Y. et al. Splicing factor SF3b as a target of the antitumor natural product pladienolide. Nat. Chem. Biol. 3, 570575 (2007).

30. Bonnal, S., Vigevani, L. & Valcarcel, J. The spliceosome as a target of novel antitumour drugs. Nat. Rev. Drug Discov. 11, 847859 (2012).

31. Hsu, T. Y. et al. The spliceosome is a therapeutic vulnerability in MYC-driven cancer. Nature 525, 384388 (2015).

32. Xargay-Torrent, S. et al. The splicing modulator sudemycin induces a specic antitumor response and cooperates with ibrutinib in chronic lymphocytic leukemia. Oncotarget 6, 2273422749 (2015).

33. Lee, S. C. et al. Modulation of splicing catalysis for therapeutic targeting of leukemia with mutations in genes encoding spliceosomal proteins. Nat. Med. 22, 672678 (2016).

34. Obeng, E. A. et al. Physiologic expression of Sf3b1(K700E) causes impaired erythropoiesis, aberrant splicing, and sensitivity to therapeutic spliceosome modulation. Cancer Cell 30, 404417 (2016).

35. Wan, Y. et al. Splicing function of mitotic regulators links R-loop-mediated DNA damage to tumor cell killing. J. Cell Biol. 209, 235246 (2015).

36. Grifth, M. et al. Genome modeling system: a knowledge management platform for genomics. PLoS Comput. Biol. 11, e1004274 (2015).37. Kim, D. et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36 (2013).

38. Flicek, P. et al. Ensembl 2013. Nucleic Acids Res. 41, D48D55 (2013).39. Gentleman, R. C. et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 5, R80 (2004).

40. Anders, S. & Huber, W. Differential expression analysis for sequence count data. Genome Biol. 11, R106 (2010).

41. Young, M. D., Wakeeld, M. J., Smyth, G. K. & Oshlack, A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11, R14 (2010).

42. Shen, S. et al. rMATS: robust and exible detection of differential alternative splicing from replicate RNA-Seq data. Proc. Natl Acad. Sci. USA 111, E5593E5601 (2014).

43. Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate-a practical and powerful approach to multiple testing. J. R. Stat. Soc. B Methodol. 57, 289300 (1995).

44. Anders, S., Reyes, A. & Huber, W. Detecting differential usage of exons from RNA-seq data. Genome Res. 22, 20082017 (2012).

45. Klco, J. M. et al. Functional heterogeneity of genetically dened subclones in acute myeloid leukemia. Cancer Cell 25, 379392 (2014).

46. Edgar, R., Domrachev, M. & Lash, A. E. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 30, 207210 (2002).

Acknowledgements

Support was provided by NIH/NHLBI (T32HL007088 to C.L.S.), Barnes-Jewish Hospital Foundation (B.S.W., T.A.G. and M.J.W.), an NIH/NCI SPORE in Leukemia (P50CA171963 to C.L.S., B.S.W., T.A.G. and M.J.W.), an NIH/NCI grant (K12CA167540 to C.L.S. and B.S.W.), a Clinical and Translational Award from the NIH National Center for Advancing Translational Sciences (UL1 TR000448 to B.S.W.), the Edward P Evans Foundation (T.A.G. and M.J.W.), the Lottie Caroline Hardy Trust (T.A.G. and M.J.W.), a Leukemia and Lymphoma Society Scholar Award (M.J.W.) and Translational Research Award (T.A.G.), by NIH grant CA140474 (T.R.W.) and Department of Defense (BM120018; M.J.W.). Support for procurement of human samples was provided by an NIH/NCI grant (P01 CA101937). Technical assistance was provided by the Alvin J Siteman Cancer Center High Speed Cell Sorting Core, the Tissue Procurement Core supported by an NCI Cancer Center Support Grant (P30CA91842), Carla Weinheimer and Mouse Cardiovascular Phenotyping Core in the Center for Cardiovascular Research at Washington University School of Medicine, and the McDonnell Genome Institute (Director, Richard Wilson and Co-Director, Elaine Mardis) and Chris Markovic for sequencing and NanoString experiments, respectively.

NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14060

E7107 and technical assistance were kindly provided by Silvia Buonamici and Peter Smith (H3 Biomedicine). We are grateful to Drs Tim Ley, Dan Link, and John DiPersio for helpful scientic discussions.

Author contributions

The study was designed by: C.L.S., B.S.W., M.T., T.A.G. and M.J.W. Primary cell/cell line experiments performed by: C.L.S., M.T., R.T., S.K., J.S., A.C. and B.S. Mouse model experiments by: C.L.S., R.T., J.N.L., M.N. and S.K. RNA sequencing by: R.S.F., C.F. and M.O. Bioinformatics analysis by: B.S.W. Sudemycin drug development and synthesis by: C.L. and T.R.W. The manuscript was written and edited by: C.L.S., B.S.W., T.A.G. and M.J.W. All co-authors reviewed and approved the submission.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Shirai, C. L. et al. Mutant U2AF1-expressing cells are sensitive to pharmacological modulation of the spliceosome. Nat. Commun. 8, 14060 doi: 10.1038/ncomms14060 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

10 NATURE COMMUNICATIONS | 8:14060 | DOI: 10.1038/ncomms14060 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications