Targeting the spliceosome through RBM39 degradation results in exceptional responses in high-risk neuroblastoma models

Description


INTRODUCTION
High-risk neuroblastoma has a poor prognosis and leads to a disproportionate number of cancer-related deaths in children, despite the use of increasingly intensive, multimodal therapies (1). Survivors of high-risk disease have a significant risk of developing long-term complications from cytotoxic chemotherapy and radiotherapy (2,3). Safer, more effective therapies are needed for highrisk neuroblastoma patients. However, the low number of targetable recurrent mutations in neuroblastoma remains an obstacle to the development of targeted therapies with high therapeutic indices.
MYC is one of the nuclear factors most often associated with dysregulated transcription, as genetic abnormalities of the MYC family of oncogenes (C-MYC, MYCN, and MYCL1) are among the most common alterations in human cancer (4,5). Neuroblastoma, a malignancy of arrested development of the peripheral nervous system (1,6), is a MYC-driven cancer (7)(8)(9). Half of high-risk neuroblastomas harbor MYCN amplification, and the other half overexpress C-MYC (10). Without a strategy to target the MYC protein directly, investigation has focused on identifying the pathways that MYC requires to drive cancer growth. To that end, pre-mRNA splicing has emerged as a potential vulnerability in MYC-driven cancers (11)(12)(13)(14)(15)(16)(17)(18), including neuroblastoma (19). Splicing abnormalities are present in both human neuroblastomas and a mouse model of MYCN-driven neuroblastoma (19)(20)(21)(22). The dependence of MYC-driven cancers on components of the spliceosome and the presence of splicing abnormalities in MYC-driven neuroblastoma raise the possibility that splicing factors could be effective targets in neuroblastoma.
The aryl sulfonamide drug indisulam has been shown to act as a "molecular glue" that bridges the U2AF-related splicing factor RBM39 to the CRL4-DCAF15 E3 ubiquitin ligase for proteasomal degradation of RBM39 (23)(24)(25)(26)(27)(28). Recruitment of RBM39 to CRL4-DCAF15 causes RBM39 ubiquitination followed by proteasomal degradation, leading to splicing abnormalities within treated cells. The amount of RBM39 degradation is proportional to the expression of DCAF15, and therefore, the sensitivity of cells to indisulam is expected to depend on both RBM39 dependence and DCAF15 expression. Here, we demonstrate that neuroblastoma is both dependent on RBM39 and expresses relatively high levels of DCAF15, providing a mechanistic rationale for indisulam treatment. Accordingly, when compared to other lineages, cultured neuroblastoma cells are significantly more sensitive to indisulam. Moreover, indisulam-mediated RBM39 degradation leads to complete tumor regression in multiple, independent mouse models of high-risk neuroblastoma. These observations provide a scientific rationale for treating high-risk neuroblastoma with small-molecule RBM39 degraders like indisulam.

RBM39 is required for neuroblastoma cell survival and tumor growth
DEMETER2 is a model that scores gene dependencies in 712 cancer cell lines that have been evaluated in three different large-scale pooled RNA interference (RNAi) screens (30), providing an additional opportunity to assess the dependency of neuroblastoma on RBM39. Although RBM39 is essential to most cancer cell lines, neuroblastoma was the lineage most significantly sensitive to RBM39 RNAi compared to all other lineages (P = 0.0023) ( Fig. 1D and table  S3). Among the splicing factors directly interacting with RBM39 ( fig. S1A), genetic knockdown of many of them led to greater inhibition of neuroblastoma than other cancer lineages ( fig. S2), supporting the hypothesis that neuroblastoma is relatively more dependent on splicing factors related to RBM39. We further validated these findings by stable expression of two different short hairpin RNAs (shRNAs) designed to target RBM39 in four different neuroblastoma models, all of which are MYC-driven as a result of either MYCN amplification (BE2C, KELLY, and SIMA) or C-MYC overexpression (SK-N-AS) (Fig. 1E). RBM39 depletion led to significantly reduced neuroblastoma cell colony formation in vitro (Fig. 1F) and tumor growth in vivo (Fig. 1G). An analysis of the transcriptome from BE2C and SK-N-AS neuroblastoma cells with RBM39 depletion demonstrated alterations in pre-mRNA splicing (particularly exon skipping) in comparison to vector controls (Fig. 1H, fig. S3, A to H, and tables S4 and S5). For example, the EED gene that encodes a critical component of polycomb-repressive complex 2 that plays a critical step in neuroblastoma growth (31,32) showed exon 8 and 9 skipping after RBM39 depletion (Fig. 1I). By comparing with the cancer genes from COSMIC (Catalogue Of Somatic Mutations In Cancer) database, we found that a number of misspliced genes by RBM39 knockdown were overlapped with COSMIC cancer genes (table S6), which are involved in a variety of key cancer pathways such as cell cycle (CDK4 and CCND3), apoptosis (BAX and BIRC6), DNA repair (ATM, ATR, and MLH1), tyrosine kinase receptormediated signal transduction (NTRK1, FGFR1, FGFR4, RET, and MAP2K2), epigenic regulation (EED, EZH2, DNMT3A, ARID1B, and KDM6A), and mitochondrial metabolism (SDHC and SDHAF2), suggesting that RBM39 is a critical factor in the oncogenesis of neuroblastoma. Real-time PCR for specific exons spliced by RBM39 knockdown validated the RNA sequencing (RNA-seq) results (fig. S3, I to L). Western blot showed that the expression of these targets was reduced by RBM39 depletion (fig. S3, M and N). The splicing changes we observed for cells expressing each individual shRNA in both cell lines showed a large number of common splicing events and significant correlation with each other, as indicated by the Spearman correlation and Venn diagram analyses ( fig. S4), suggesting that the observed alterations were specific to RBM39 suppression. Gene set enrichment analysis (GSEA) for the differentially expressed genes showed that RBM39 depletion significantly disrupts pathways that regulate RNA and protein homeostasis (fig. S5, A and B, and table S7). Collectively, these data indicate that neuroblastoma has an inordinate dependency on RBM39-mediated RNA splicing, which may represent a vulnerability in MYC-driven high-risk neuroblastoma ( Table 1).

The RBM39 degrader indisulam is selective and highly active against neuroblastoma in vitro and in vivo
Indisulam and other related aryl sulfonamides represent a class of compounds known as "molecular glues" because of their ability to recruit the splicing factor RBM39 to the E3 ubiquitin ligase CRL4-DCAF15 for proteasomal degradation (23,28). We reasoned that indisulam might be an effective agent for targeting the splicing program in MYC-driven neuroblastoma by degrading RBM39. Profiling Relative Inhibition Simultaneously in Mixtures (PRISM) and the Cancer Therapeutic Response Portal (CTRP) are two independent large-scale datasets that tested the impact of many drugs, including indisulam, on hundreds of different cell lines, providing a resource to identify neuroblastoma-selective compounds. Primary PRISM evaluated the impact of 4686 different small molecules on the proliferation of up to 560 human cancer cell lines. A secondary analysis measured the sensitivity of 148 compounds using a dose titration and quantified sensitivity as assessed by determining the area under the dose-response curve (AUC). Ranking the compounds in order of selectivity for neuroblastoma revealed indisulam to be third out of 4686 and first out of 148 compounds in the primary and secondary analyses, respectively (Fig. 2, A and B; fig. S6, A and B; and tables S8 and S9). In the CTRP analysis of 709 cell lines, the indisulam AUC for neuroblastoma showed it to be significantly smaller (more sensitive) than all other cell lineages (Fig. 2C, fig.  S6C, and table S10). Consistent with these predictions, we found that indisulam induced dose-dependent reduction of RBM39 (Fig. 2D) and widespread splicing anomalies. The splicing anomalies were due to RBM39 degradation, because genome-wide pairwise comparisons between cells expressing RBM39 shRNA and those treated with indisulam showed a significant correlation (Fig. 2, E to G), in line with the  (F) Colony formation for neuroblastoma cells transduced with lentiviral shRNA control and two different RBM39 shRNAs followed by puromycin selection for 7 days. Cells stained with crystal violet for visualization. (G) Tumor volume comparison between shRNA control and RBM39 knockdown by two different shRNAs in four tumor models. P values calculated using unpaired Student's t test. (H) RNA splicing analyses for different events altered by RBM39 knockdown in BE2C and SK-N-AS cells. The y axis indicates numbers for each splicing event induced by RBM39 knockdown in comparison with shRNA control. A3SS, alternative 3′ splicing; A5SS, alternative 5′ splicing; MXE, mutually exclusive exon; RI, intron retention; SE, skipped exon. (I) Sashimi plot using IGV program shows the splicing changes of EED gene in BE2C and SK-N-AS cells after RBM39 knockdown, which indicates the missplicing among exons 7 to 10. The numbers indicate the read counts of exon-intron junctions of RNA-seq. common splicing events by Venn diagram analysis ( fig. S7). A PrestoBlue assay for viability showed that neuroblastoma cell lines were more sensitive than select normal human cell lines BJ and HS68 (Fig. 2H). In addition, indisulam reduced colony formation in 10 different neuroblastoma cell lines (Fig. 2I). Flow cytometry analysis showed that indisulam induces G 2 -M phase arrest and cell death, as indicated by an increase in sub-G 1 phase cells ( fig. S8). Pathway analysis of the genes with altered splicing by indisulam and RBM39 knockdown revealed that loss of RBM39 affected pathways involved in G 2 -M phase cell cycle, gene transcription, DNA repair, chromatin-modifying enzymes, metabolic pathways, and cancer genes ( fig. S9). Real-time PCR and Western blot validated that the splicing events induced by indisulam were consistent with those induced by RBM39 knockdown (fig. S10, A to N).
We then tested the in vivo efficacy and selectivity of indisulam in various xenograft models, including neuroblastoma (BE2C and SK-N-AS), rhabdomyosarcoma (Rh30), and colorectal cancer (HCT116). Mice harboring tumors of approximately 200 mm 3 were dosed intravenously with 25 mg/kg for two 1-week cycles (5 days on and 2 days off). Although all tumors in all four models underwent an initial complete response following treatment, the time to relapse, following discontinuation of therapy, in neuroblastoma was substantially longer in comparison to the Rh30 and HCT116 models (Fig. 2, J to M). Together, these data suggest that neuroblastoma is more sensitive to RBM39 degrader indisulam in comparison with other cancer lineages.

Indisulam efficacy against neuroblastoma is specifically through targeting RBM39
To test whether the antitumor activity of indisulam in neuroblastoma is an RBM39-dependent effect rather than due to off-target activity, we engineered each of three neuroblastoma cell lines (SIMA, BE2C, and SK-N-AS), using CRISPR-Cas9 knock-in, to create indisulamresistant RBM39 variants (RBM39 G268V) (Fig. 3A). A glycine to valine substitution in RBM39 (G268V) prevents its recruitment to DCAF15 and renders it resistant to degradation without affecting its essential role in pre-mRNA splicing (23). G268V-mutant cell lines were resistant to indisulam-mediated RBM39 degradation (Fig. 3B). Moreover, these cells were more than 100-fold resistant to the antiproliferative activity of indisulam (Fig. 3C). RBM39 G268V mutations also rendered neuroblastoma cells resistant to indisulam in colony formation assays (Fig. 3, D and E) and xenograft tumor regression (Fig. 3, F and G). The on-target action of indisulam in tumor cells was further demonstrated by resistance to indisulaminduced RBM39 degradation in vivo (Fig. 3H). On the basis of these results, we concluded that indisulam-mediated degradation of RBM39 is the mechanism of action that accounts for its antitumor effect in neuroblastoma both in culture and in vivo.

High expression of DCAF15 in neuroblastoma is correlated with indisulam activity
A prerequisite for indisulam-mediated RBM39 degradation is DCAF15 expression, and the degree of RBM39 degradation is proportional to the level of DCAF15 expression (23). After establishing the dependency of neuroblastoma on RBM39, we investigated what role DCAF15 expression levels play in determining the high activity of indisulam in neuroblastoma over other cancer lineages. Among cell lines representing 26 different cancer lineages, neuroblastoma cells expressed relatively high levels of DCAF15 mRNA, comparable to blood and lymphoma lineages ( Fig. 4A and table S11). In agreement, DCAF15 expression levels are higher in neuroblastoma tumors when compared to multiple adult solid tumors (fig. S11). We further compared DCAF15 expression between neuroblastomas and other types of pediatric cancers, among 2222 cases analyzed by

of 17
RNA-seq (St. Jude Cloud) (33), and found that neuroblastoma had the second highest levels of DCAF15, behind only retinoblastoma ( Fig. 4B and table S12). Relatively high DCAF15 expression may, in part, explain the exquisite sensitivity of neuroblastoma cancers to indisulam. Consistent with this hypothesis, among 709 different cancer cells in the CTRP, we found a significant correlation between indisulam sensitivity and DCAF15 expression (Fig. 4C). To directly test whether the antiproliferative effect of indisulam in neuroblastoma cells depends on DCAF15 expression, we engineered neuroblastoma cell lines to either overexpress DCAF15, silence DCAF15, or rescue DCAF15 in DCAF15 knockout (KO) cells (Fig. 4D). The genetic perturbations  , and their derivatives (DCAF15 3XF and DCAF15 KO) that were treated with indisulam (5 mg/kg; n = 5 per group) for 5 days on and 2 days off, two cycles via tail injection. The wild-type group data were generated from the same group of mice as used in Fig. 3 (F and G). Cyan arrow indicates when therapy was discontinued.
of DCAF15 did not affect expression of endogenous RBM39 (fig. S12). However, indisulam-mediated RBM39 degradation did correlate with DCAF15 levels. DCAF15 KO reduced degradation, whereas DCAF15 overexpression enhanced degradation (Fig. 4D). The effects of indisulam on the viability of these cells in vitro mirrored the changes observed with RBM39 degradation (Fig. 4, E and F). The efficacy of indisulam in reducing tumor growth in xenograft studies was higher in cells ectopically expressing DCAF15 and lower in DCAF15-silenced cells (Fig. 4, G and H). Together, these data demonstrate that DCAF15 expression levels determine the efficacy of indisulam with high levels of DCAF15 in neuroblastoma conferring more sensitivity to indisulam treatment.

Indisulam treatment results in significant response in high-risk neuroblastoma models
On the basis of above results, we concluded that indisulam-mediated degradation of RBM39 is a viable strategy to target pre-mRNA splicing in MYC-driven neuroblastoma. To further gauge the therapeutic efficacy of indisulam, we administered indisulam to immunodeficient mice harboring xenografts derived from patients with relapsed neuroblastoma, which have never been cultured in vitro, with three different doses of indisulam (5, 12.5, and 25 mg/kg). Indisulam led to a dose-dependent regression of tumor growth, although animals treated with the lowest dose, 5 mg/kg, had evidence of early tumor relapse (Fig. 5, A and B). The comparative studies demonstrate the conservation of alternative (61%) and constitutive (74%) splice junctions in mouse and human genomes, and splicing mechanism is largely conserved between human and mouse (34). We therefore hypothesized that mouse neuroblastomas driven by MYC will have a dependency on RBM39 analogous to human neuroblastoma. To test this hypothesis, we evaluated the efficacy of indisulam in an immune-competent, transgenic mouse neuroblastoma model driven by MYCN and ALK F1178L (TH-MYCN/ALK F1178L ) by using ultrasound imaging to quantify tumor response. Human ALK F1174L (analogous to mouse ALK F1178L ) is the most frequent somatic mutation in neuroblastoma and is associated with MYCN amplification, conferring a worse prognosis than MYCN amplification alone (35). Indisulam led to a complete response shortly after treatment of TH-MYCN/ALK F1178L mice even when tumor volume at initiation of treatment exceeded 2000 mm 3 (Fig. 5C). Tumors collected from treated mice before regression displayed a substantial reduction in RBM39 levels, consistent with an on-target effect (Fig. 5C, inlet). This response has been durable for over 50 weeks in four of six mice (Fig. 5, D and E, and fig. S13), despite only having received therapy for 2 weeks. We further tested the efficacy of indisulam in a syngeneic C-MYC neuroblastoma model, in which tumors derived from Dbh-iCre/CAG-C-MYC mice were engrafted in immune-competent mice of the same genetic background. Indisulam was administered when tumor sizes were between 2000 and 3500 mm 3 . Again, indisulam treatment resulted in durable complete response in all tumors (Fig. 5F). Together, these data reveal that indisulam induces a complete and durable antitumor response in different MYC-driven high-risk neuroblastoma models.

Indisulam shows no overt toxicity in mouse models
To monitor target engagement in patient-derived xenograft (PDX) mouse models, we used Western blot analysis to determine the levels of RBM39 in tumor tissues as well as normal organs including heart, lungs, liver, and kidney 120 hours following indisulam administration (24 hours after the fourth dosing). In tumors, we observed a reduction in RBM39 at all three doses administered in two different PDX models (Fig. 6, A and B). In contrast, there was no consistent reduction in RBM39 in heart, lungs, and liver, although there was some reduction in the kidney (Fig. 6, A and B), which may be a reflection of DCAF15 expression levels. Among all tested animals, we did not observe a significant change in mouse body weight during treatment (Fig. 6C). Although a phase 1 study of indisulam showed that nearly all severe adverse events were hematological toxicities (thrombocytopenia, neutropenia, and anemia) (36,37), we found no substantial changes in blood cell counts (neutrophils, platelets, red blood cells, and hemoglobin) during either the first or second cycle of treatment when using the dosing schedule of 25 mg/kg (Fig. 6, D and E). These data indicate that indisulam efficiently targets RBM39 in tumors but, at least in mice, not in most normal tissues examined.

DISCUSSION
Spliceosome-mediated pre-mRNA editing is an essential biological process that produces mature RNAs that serve as templates for protein synthesis. Alterations in RNA splicing pathways can lead to dysregulated gene splicing and a tumor-specific dependence (38,39). The MYC proto-oncogene is an amplifier of transcription of numerous genes including those involved in cell growth, cell cycle, and metabolism. During MYC-driven oncogenic transformation, cells must adapt to the resulting increased burden in RNA processing. Here, we show that pathways involving RNA metabolism and splicing are the most differentially expressed and required in high-risk neuroblastoma. We observed that one of these splicing factors, in particular, RBM39, is essential to RNA splicing in neuroblastoma cells. In genome-wide RNAi screens, suppression of RBM39 had a greater impact on neuroblastoma proliferation than other lineages. The dependency of neuroblastoma on RBM39 is consistent with the concept that RNA splicing may be an "Achilles' heel" of MYC-driven cancer (40). The mechanism for why RBM39 is especially essential in MYC-driven neuroblastoma is an important question to address in the future and has the potential to reveal new targets. We noticed that RBM39 depletion resulted in missplicing of a number of cancer genes that are involved in cell cycle, apoptosis, signal transduction, DNA repair, metabolism, and epigenetic modifications. Missplicing of these candidate genes may explain the greater degree of RBM39 dependence in neuroblastoma than other lineages. Alternatively, these genes could be synthetically lethal with MYC activation in neuroblastoma. The spliceosome involves nearly 300 proteins that form the megadalton complex, which undergoes dynamic assembly and disassembly of at least five subcomplexes in a stepwise manner during pre-mRNA splicing process (41,42). While RBM39 has been found to physically interact with U2AF65 (43,44), an important component of subcomplex E, A, and B of spliceosome, it also interacts with numerous splicing factors whose functions remain to be defined (Fig. 1C) (23). RBM39 degradation leads to selective effects on pre-mRNA splicing affecting only a subset of genes, and the precise biochemical complex in which RBM39 regulates these events is currently unknown. Defining the other components of this complex mechanism has the potential to unveil more targets in neuroblastoma. Therefore, which step and which subcomplex(es) of spliceosome will be affected by RBM39 depletion and indisulam treatment is a critical question and warrants further studies.

of 17
The dependence of neuroblastoma on RBM39 provides a rationale to target the altered splicing machinery in neuroblastoma using an RBM39 degrader as a therapeutic approach. Targeted protein degradation enables conditional and temporal control over the levels of proteins. This approach has emerged as a new therapeutic modality owing to its potential advantages over traditional occupancy-based inhibitors with respect to dosing, side effects, drug resistance, and modulating "undruggable" targets (45). Indisulam is one such molecule that recruits RBM39 as a selective neosubstrate to CRL4-DCAF15 E3 ubiquitin ligase, leading to proteasomal degradation (23). In this study, we show that indisulam is particularly effective against neuroblastoma in comparison with other cancer lineages. We have observed significant and durable responses to indisulam in multiple high-risk neuroblastoma models including PDXs derived from relapsed patient samples, transgenic MYCN/ALK F1178L mice, and tumors derived from C-MYC-overexpressing mice. We validated the on-target mechanism of indisulam in vitro and in vivo in several ways. First, we showed that cells expressing the RBM39 mutant G268V, which is incapable of binding indisulam but maintains its functional role in the spliceosome, were about 100-fold more resistant to indisulam in vitro than the parental cells and unresponsive to indisulam treatment in vivo. Second, genetic knockdown and pharmacologic degradation of RBM39 induced similar splicing changes, mainly leading to genome-wide skipped exons. Nevertheless, we noticed that the indisulam treatment induced more robust changes in splicing than the shRNA knockdown of RBM39 (figs. S3,  S7, and S10), which could be attributed to several reasons. First, indisulam might have other "off-target" effects. Second, indisulam may more potently destabilize the splicing complexes than shRNA knockdown due to insufficient RBM39 depletion. Third, the shRNA knockdown effect lagged behind indisulam-mediated degradation because it takes longer time for knockdown to effect in cells.
We further demonstrate that DCAF15 expression levels correlate with the efficacy of indisulam. The high expression levels of DCAF15 in neuroblastoma may explain, in combination with the aforementioned dependence on RBM39, why neuroblastoma is exceptionally sensitive to indisulam treatment in comparison with other cancer lineages. Thus, our data support the rationale for patient stratification based on the expression of DCAF15, as a biomarker in future clinical trials. The high levels of DCAF15 and RBM39 dependence in neuroblastoma might make indisulam especially efficacious against this MYC-driven cancer.
A phase 1 clinical study of indisulam showed that thrombocytopenia was the major dose-limiting toxicity followed by neutropenia (36). This was confirmed by a phase 2 clinical study showing that nearly all severe adverse events were hematological toxicities (thrombocytopenia, neutropenia, and anemia). The most common grade 3 and 4 nonhematologic toxicities were transient electrolyte abnormalities (34). In our pharmacodynamic studies, all PDXs showed substantial reduction of RBM39 in tumor tissues in response to three different indisulam doses. In contrast, there was no comparable reduction of RBM39 in heart, lungs, and liver, although there was some reduction in the kidney, which may be a reflection of higher DCAF15 expression levels in this tissue. In addition, we observed no whole-body weight loss of treated animals, which is a common indicator of adverse effects, during the treatment course. Whole-blood counts revealed that the three different dosing schedules had no significant adverse effect on the number of neutrophils, platelets, and hemoglobin or red blood cells, in contrast to human patients. The different toxicity profiles may be due to an off-target effect of indisulam in human patients or a mechanismrelated effect that is species specific. Although neuroblastoma expresses high levels of DCAF15 and is highly responsive to indisulam in comparison with other cancer lineages, it might be necessary in the future to stratify patients who are most likely to have the best response. It will be interesting to see whether any genetic surrogates could serve as biomarkers for indisulam treatment. It is worth noting that the neuroblastoma mouse model driven by MYCN/ALK F1174L and C-MYC had even better outcomes following indisulam treatment than the PDX and cell line-generated models, which are conducted in immunodeficient mice. This observation requires follow-up investigation because it raises the possibility that indisulam's efficacy might be enhanced by the immune system.
In clinical trials performed before understanding the mechanism of action of indisulam, only a small population of patients responded to indisulam. However, patients with neuroblastoma were not included in these trials. Here, we show that the RBM39 dependence and high levels of DCAF15 expression in neuroblastoma account for exceptional efficacy and a high therapeutic index of indisulam in a variety of high-risk neuroblastoma models. Thus, a clinical trial of indisulam for high-risk neuroblastoma is warranted, and preclinical studies are needed to evaluate indisulam in other solid tumors such as ovarian cancers with high MYCN levels (46).

Derivation of indisulam-resistant SK-N-AS cells from indisulam-resistant SK-N-AS tumors
The vehicle-treated SK-N-AS tumor (n = 3) and indisulam-resistant SK-N-AS tumors (n = 3) were excised and placed in a sterile tube containing phosphate-buffered saline (PBS) on wet ice during transport from the animal research facility to the research laboratory. Tumor samples were manually minced using a sterile scalpel and underwent an enzymatic digestion with collagenase IV (2 mg/ml; in 25 ml of RPMI medium) for 1 hour in a 37°C rotor (Robbins Scientific Corporation, model 2000). After digestion, cells were filtered using a 70-m sterile strainer and cultured in RPMI medium with 10% FBS and 1% penicillin and streptomycin.

C A A T G G C A A A C A A T T T A C A A A A G G G A A G T G C T G -GACCTATGAGGCTTTATGTGGGCTCATTACACTTCAA-CATAACTGAAGATATGCTTCGCGTGATCTTTGAG-C C T T T T G G A A G A G T A A G T C C A G G T T C T T C A A -
TGAATCTTCAGTAGGTTGTTGATCTGAGTATAACTA-CATAATTTCTGTATGTCTTGCT-3′. The oligonucleotide is engineered to introduce a synonymous mutation at R267, allowing us to confirm that the genomic sequence was edited with this template. Afterward, cells were exposed to 2 M indisulam for 2 weeks to select for cells with RBM39 G268V knock-in. Cells that survived indisulam treatment were recovered and confirmed to have the correct G268V genomic conversion via amplicon sequencing. Amplicon sequencing of PCR products using genomic DNA as template revealed BE2C (62%), SK-N-AS (75%), and SIMA (55%) reads with the repair template sequence, consistent with a heterozygous effect. HCT116 cells expressing RBM39 mutations were previously described by Han et al. (23).

Ectopic expression of 3xFLAG-tagged DCAF15
Full-length complementary DNA (cDNA) encoding DCAF15 with a 3xFLAG epitope at the C-terminal end was cloned into pLVX IRES Blast. The PAM (protospacer adjacent motif) site was mutated to prevent Cas9 from editing ectopic expression of 3xFLAG-tagged DCAF15 in SK-N-AS, BE2C, and SIMA cells generated by CRISPRmediated deactivation of DCAF15. Lentiviral packaging of the resulting plasmid was performed by cotransfecting with psPAX2 (Addgene, 12260) and pMD2.G (Addgene, 12259) into lenti-X 293T cells. The medium collected from transfected lenti-X 293T cells was used to infect SK-N-AS, BE2C, and SIMA cells. Cells stably expressing 3xFLAG-tagged DCAF15 were selected for 7 to 10 days with blasticidin (10 g/ml). HCT116 cells were described by Han et al. (23). Cells were split each into one six-well plate for an 18-hour incubation with indisulam at doses 0 [dimethyl sulfoxide (DMSO)], 0.05, 0.1, 0.5, 1, and 5 M to validate for 3xFLAGtagged DCAF15 and RBM39 degradation via Western blot.

Sanger DNA sequencing of RBM39
Total RNA was isolated from tissue samples as per the manufacturer's instructions of the RNeasy Mini Kit (Qiagen, catalog no. 74106). cDNA was prepared from 1 g of total RNA using a cDNA synthesis kit according to the manufacturer's protocol (Thermo Fisher Scientific, catalog no. 18091050). RBM39 gene was amplified by PCR using the cDNA as the template along with the following primers (forward: 5′-atttctagagccaccatggcagacgatattgatattgaagcaatgc-3′; reverse: 5′-attggatcctcatcgtctacttggaaccagtagc-3′) using the Phusion High Fidelity PCR Kit (catalog no. E0553S). The resultant PCR product was gel-purified using a Qiagen gel extraction kit (Qiagen, catalog no. 28706) according to the manufacturer's protocol. The DNA was sent for Sanger sequencing at Hartwell Center in St. Jude using the following primers: primer 288F, ACAGAAGTCCT-TACTCCGGACC; primer 388R, ACTTTTGCTTCGGGAACGTCG; primer 602F, GTCGATGTTAGCTCAGTGCCTC. Primer 602F is used to detect the RBM39 mutation in the RRM2 motif.

SDS-polyacrylamide gel electrophoresis and Western blot
Cells were washed twice with ice-cold PBS and directly lysed on ice with 2× sample loading buffer [0.1 M tris HCl (pH 6.8), 200 mM dithiothreitol, 0.01% bromophenol blue, 4% SDS, and 20% glycerol]. On ice, cell lysates were sonicated once with a 5-s bursts at 40% amplitude output (Sonics, VIBRA CELL) followed by 25-min heating at 95°C. After the cell lysates were centrifuged at 13,000g at room temperature for 2 min, 10 to 20 l of the cell lysates were separated on 4 to 15% Mini-PROTEAN TGX Stain-Free Protein Gels from Bio-Rad and transferred to methanol-soaked polyvinylidene difluoride membranes (Millipore). Lysates for RBM39 G268V mutant cell lines and DCAF15 genetically modified cells were generated as previously described (23). Membranes were blocked in PBS buffer supplemented with 0.1% Tween 20 and 5% skim milk (PBS-T) and incubated for 1 hour at room temperature under gentle horizontal shaking. Membranes were incubated overnight at 4°C with the primary antibodies. The next day, membranes were washed three times (for 5 min) with PBS-T at room temperature. Protected from light, membranes were then incubated with goat anti-mouse or goat anti-rabbit horseradish peroxidase-conjugated secondary antibodies (1:5000) for 1 hour at room temperature, followed by three 5-min washes with PBS-T at room temperature. Last, membranes were incubated for 1 min at room temperature with SuperSignal West Pico PLUS Chemiluminescent Substrate (Thermo Fisher Scientific, 34580), and the bound antigen-antibody complexes were visualized using the Odyssey Fc Imaging System (LI-COR Corp., Lincoln, NE).

Crystal violet staining
After removing medium, cells were washed with Dulbecco's PBS (DPBS) without calcium or magnesium (DPBS, Lonza) and treated with 4% formaldehyde in PBS [paraformaldehyde (PFA)] for 20 min. Once PFA was removed, cells were stained with 0.1% crystal violet stain for 1 hour.

Cell viability assay
Twelve-point dose responses were performed on 96-well assay plates with cell plating (1500 to 4000 cells per well) on day 1, compound addition on day 2, and survival measurement on day 5. Compounds were diluted in DMSO before adding to the cells. Final DMSO concentration was 0.5%. Cell survival assay was performed using a CellTiter-Glo Luminescent Cell Viability Assay kit (Promega) that measures cellular adenosine triphosphate (ATP) content. Luminescence was recorded with an EnVision multimode plate reader (PerkinElmer). Median inhibitory concentration (IC 50 ) was determined with GraphPad Prism using baseline correction (by normalizing to DMSO control), the asymmetric (five parameter) equation, and least-squares fit.

PrestoBlue assay
Cell viability was determined by PrestoBlue Cell Viability Reagent (Invitrogen, catalog no. A13262) according to the manufacturer's protocol. Briefly, cells were plated in white 96-well plates (PerkinElmer, catalog no. 6005680) at a density of 1000 to 3000 cells per well and incubated at 37°C with 5% CO 2 for 24 hours. Afterward, cells were treated with drugs at different final concentrations (10,000, 500, 250, 50, 25, 12.5, 5, 2.5, 0.5, 0.05, and 0 nM) for 5 days. Each concentration was tested in eight replicates. After completion of the treatment period, 10 l of PrestoBlue reagent was added to each well and plates were incubated at 37°C with 5% CO 2 for 30 min. Then, fluorescence (560-nm excitation/590-nm emission) was measured using a BioTek Synergy H1 microplate reader. GraphPad Prism was used to plot the curves and calculate the IC 50 values.

Cell cycle analysis
Cells (100, 000 per well) were plated in six-well culture plates and, next day, treated with indisulam for 48 hours. Cells were trypsinized and centrifuged at 1000 rpm for 5 min at 4°C and washed once with ice-cold PBS. The samples were then resuspended in the propidium iodide solution [propidium iodide (0.05 mg/ml), 0.1% sodium citrate, and 0.1% Triton X-100] and incubated at room temperature for 20 min. The cells were subjected to flow cytometric analysis by using an LSR II (BD Biosciences, CA, USA) flow cytometer, and the data were analyzed using BD FACSDiva (BD Biosciences, CA, USA).

Differential gene expression and GSEA for RNA-seq experiments
Total RNA from cells and tumor tissues was performed using the RNeasy Mini Kit (Qiagen) according to the manufacturer's instructions. Paired-end sequencing was performed using the High-Seq platform with 100-base pair (bp) read length. Total stranded RNA-seq data were processed by the internal AutoMapper pipeline. Briefly, the raw reads were first trimmed (Trim-Galore version 0.60) and mapped to human genome assembly (GRCh38) (STAR v2.7), and then the gene-level values were quantified (RSEM, v1.31) on the basis of GENCODE annotation (v31). Low-count genes were removed from analysis using a CPM (counts-per-million) cutoff corresponding to a count of 10 reads and only confidently annotated (level 1 and 2 gene annotation), and protein-coding genes are used for differential expression analysis. Normalization factors were generated using the TMM (Trimmed Mean of the M-values) method, counts were then transformed using Voom, and transformed counts were analyzed using the lmFit and eBayes functions (R limma package version 3.42.2). The significantly up-and down-regulated genes were defined by at least twofold changes and adjusted P value of <0.05. Then, GSEA was carried out using gene-level log 2 fold changes from differential expression results against gene sets in the Molecular Signatures Database (MSigDB 6.2) (gsea2 version 2.2.3) for the data from GSE62564.

RNA splicing analysis
After mapping RNA-seq data, rMATS v4.1.0 was used for RNA alternative splicing analysis by using the mapped BAM files as input. Specifically, five different kinds of alternative splicing events were identified, i.e., skipped exon, alternative 5′ splicing site, alternative 3′ splicing site, mutually exclusive exon, and intron retention. To keep consistent, the same GTF (General Feature Format) annotation reference file for mapping was used for rMATS. For stranded RNA-seq data, the argument "--libType fr-firststrand" was applied. To process reads with variable lengths, the argument "--variableread-length" was also used for rMATS. To select statistically significantly differential splicing events, the following thresholds were used: false discovery rate (FDR) < 0.05 and absolute value of IncLevel-Difference > 0.1. For visualization, the Integrated Genomic Viewer (IGV) Genome Browser was used to show the sashimi plots of splicing events. To investigate the genome-wide correlations of differential splicing between two genotypes (e.g., shRNA knockdown of RBM39 and indisulam-mediated degradation of RBM39 in cells), we extracted splice junctions for all samples of both genotypes of interest from the STAR (47) output files suffixed with "SJ.out.tab," which contain high-confidence collapsed splice junctions. Only those unique mapped reads crossing the junctions were considered. By extracting the union of the unique junction positions, we constructed a unified junction-read feature vector for each sample. Then, we normalized the junction-read vectors of each sample by their corresponding total read counts and multiplied them by 10 6 . Next, we averaged the junction-read vectors for samples of the same genotype. We further excluded those junctions for which the total normalized junction-read counts were smaller than 10. Then, we calculated the Pearson correlation coefficients for each pair of genotypes. The two-sided t tests were used to compute the significant levels of the correlations by the function chart.Correlation() in R.
To investigate the genome-wide correlations of differential splicing events between shRNA knockdown of RBM39 and indisulam treatment in different cell lines, we first used rMATS to find the differential alternative splicing events for each pair of comparisons in five splicing categories including exon skipping, mutually exclusive exons, alternative 3′ splicing site, alternative 5′ splicing site, and intron retention. Then, we selected the significant splicing events by using the criteria of FDR < 0.05 and abs(IncLevelDifference) > 0.1. Next, on the basis of these different splicing events, we used Venn diagrams to show their common and unique splicing events in the five splicing categories mentioned above.

RNA extraction and RT-PCR for RBM39
RNA was extracted using the RNeasy Plus Mini Kit (Qiagen, reference no. 74136) following the manufacturer's protocol. cDNA was prepared in 20-l reaction from 500 ng of total RNA using the Superscript IV First Strand Synthesis System (Invitrogen, reference no. 1809105). Real-time PCRs were run in triplicates (n = 3) in the 7500 Real-Time PCR System by Applied Biosystems (Thermo Fisher Scientific) using Power SYBR Green PCR Master Mix (Applied Biosystems, reference no. 4367660). C T methods were applied to analyze the results. The following primers were used to perform the quantitative real-time PCR: GAPDH (forward, AACGGGAAGCTTGTCAT-CAATGGAAA; reverse, GCATCAGCAGAGGGGGCAGAG) and RBM39 (forward, AGATGGACAACTGCCTCATTAC; reverse, GCCTCCCAGTGTTCACATATAC).

Assay for transposase-accessible chromatin using sequencing
Library preparations for ATAC-seq (assay for transposase-accessible chromatin using sequencing) were based on a published protocol with minor modifications (48,49). Briefly, freshly cultured SKNAS cells (100,000 per sample, naïve and two resistant ones) were harvested and washed with 150 l of ice-cold DPBS containing protease inhibitor (PI). Nuclei were collected by centrifugation at 500g for 10 min at 4°C after cell pellets were resuspended in lysis buffer

ATAC-seq alignment, peak calling, and annotation
The ATAC-seq raw reads were aligned to the human reference genome (hg19) using BWA (50) and then marked duplicated reads with Picard (version 1.65), with only high-quality reads kept by samtools (version 1.3.1, parameter "-q 1 -F 1024") (51). Reads mapping to mitochondrial DNA were excluded from the analysis. All mapped reads were offset by +4 bp for the + strand and −5 bp for the -strand (48). Peaks were called for each sample using MACS2 (52) with parameters "-q 0.01 -nomodel -extsize 200 -shift 100." Peaks were merged for the same cell types using BEDtools (53). Peak annotation was performed using HOMER (54). All sequencing tracks were viewed using IGV 2.3.82 (55).

Pathway network analysis
The 1585 essential fitness genes to neuroblastoma cell survival (n = 7, dependency score > 40%) identified through genome-wide CRISPR-Cas9 library screen by Sanger institute (24) were uploaded into STRING program (https://string-db.org) for network interaction analysis with confidence threshold 0.15. The resulting network was then uploaded into Cytoscape program (57), followed by gene ontology analysis using the BiNGO plugin program (58).

Analysis of PRISM data
We loaded the PRISM Repurposing 19Q4 Drug Sensitivity Primary Screen Replicate Collapsed log 2 fold change data (primary-screenreplicate-collapsed-logfold-change.csv, DepMap.org). The log 2 fold change is quantified with respect to DMSO and is a measure of cell viability in response to the chemical perturbations. Cell lines containing the string "FAILED_STR" (e.g., ACH-000010_FAILED_ STR) were removed, leaving 568 cell lines and 4686 compounds. Twelve of these cell lines are neuroblastoma lineage. Not all compounds have data available for the 568 cell lines. In the case of indisulam, data for 5 neuroblastoma cell lines and 536 nonneuroblastoma cell lines are available.

Differential sensitivity analysis of neuroblastoma cell lines
For each compound in the primary screen, we calculate the mean log 2 fold change of neuroblastoma cell lines minus the mean of log 2 fold change for non-neuroblastoma cell lines. To calculate P value of mean difference, we conducted a permutation test (i.e., Monte Carlo sampling) in which we randomly sampled five values from the indisulam log 2 fold change data and then took the difference in means of those 5 against the remaining 536 values. We randomly sampled in this manner 10 5 times. On the basis of the distribution of these 10 5 means, we concluded that indisulam mean difference of −1.30 has a P value less than 0.003 and was highly statistically significant. We performed the same calculation for the PRISM Repurposing 19Q4 Drug Sensitivity Secondary Screen AUC data (secondary-screen-dose-response-curve-parameters.csv). Again, we removed any cell lines in the data labeled with FAILED STR. For a given cell line and compound, there are AUC data from multiple screens. The documentation provided with the data (the README file) recommends using screen "MTS010" when available. The figures contain data from only this screen. Restricting to screen MTS010, we have 473 cell lines and 147 compounds, although any given compound may not have coverage across all 473 cell lines. In the case of indisulam, data were available for 4 neuroblastoma cell lines and 438 non-neuroblastoma cell lines. For each compound in the MTS010 screen, we calculated the mean AUC of neuroblastoma cell lines minus the mean AUC of non-neuroblastoma cell lines. We conducted a permutation test to calculate P value of mean difference AUC in which we randomly sampled four values from the indisulam AUC data and then took the difference in means of those 4 against the remaining 438 values. We randomly sampled in this manner 10 5 times. On the basis of the distribution of these 10 5 means, we conclude that indisulam mean difference AUC of −0.226 has a P value less than 0.001 and is highly statistically significant.
shRNA xenograft studies Neuroblastoma cells were transduced with shRNA lentiviral particles targeting RBM39. Forty-eight hours later, puromycin (1 g/ml) was added for selection for additional 72 hours (BE2C) or 48 hours (SKNAS, KELLY, and SIMA). Cancer cells (2.2 × 10 6 for BE2C, 5 × 10 6 for SIMA and SKNAS, and 2.5 × 10 6 for KELLY) were mixed with Matrigel (1:1 ratio in volume) and subcutaneously injected into the flank sites of NSG (NOD scid gamma) mice. Mice were sacrificed when they reached the humane endpoint for BE2C (day 12), SKNAS (day 14), SIMA (day 8), and KELLY (day 8) from the day measurement was started (day 1). Tumors were measured by using electronic calipers, and volumes were calculated as width /6 × d 3 , where d is the mean of two diameters taken at right angles.

Transgenic TH-MYCN/ALK F1178L mouse model
This model was provided by P. Murray, P. Thomas, and L.-A. Van de Velde, the details of which have been reported (59). The St. Jude Children's Research Hospital Institutional Animal Care and Use Committee approved all studies performed. Following genotyping, TH-MYCN/ALK F1178L mice were gender-segregated and assigned to imaging groups. Inclusion criteria were the presence of the TH-MYCN allele, heterozygosity for one ALK F1178L allele, and an age range from 1 to 3 weeks postweaning. Either gender was used and noted. Criteria for analysis were the presence or absence of an ultrasound-detectable tumor (quantified by staff not directly involved in this study) at any point within a 7-week period, where mice were imaged once per week.

Ultrasound imaging of TH-MYCN/ALK F1178L mouse
Fur was removed from the ventral side of each animal using Nair. Technicians in the St. Jude Center for In Vivo Imaging and Therapeutics performed ultrasound scanning on mice weekly using VEVO-2100 and determined tumor volumes using VevoLAB 3.0.0 software. All ultrasound data were acquired in a blinded fashion.

C-MYC syngeneic neuroblastoma model
Primary tumors derived from Dbh-iCre; CAG-MYC LSL-KI/LSL-KI were provided by R. Wang (this mouse model will be reported separately), which were subcutaneously engrafted in CB57BL6 mice. When tumor sizes reached about 2000 to 3500 mm 2 but below 20% of mouse body weight, indisulam treatment was started.

Indisulam efficacy studies in in vivo models
All murine experiments were done in accordance with a protocol approved by the Institutional Animal Care and Use Committee of St. Jude Children's Research Hospital. Subcutaneous xenografts were established in CB17 SCID mice (CB17 scid, Taconic) or NOD. Cg-Prkdcscid Il2rgtm1Wjl/SzJ (NSG; The Jackson Laboratory) mice by implanting 5 × 10 6 cells in Matrigel. We used both genders in NSG or transgenic TH-MYCN/ALK F1178L mice models to exclude sex bias. Tumor measurements were done weekly using electronic calipers, and volumes were calculated as width /6 × d 3 , where d is the mean of two diameters taken at right angles. Mice were randomly assigned to experimental groups at 4 to 6 weeks of age. A 2-week treatment commenced when the implanted tumor volume reached above 200 mm 3 . As a standard dosing schedule, indisulam (25 mg/kg) was administered via tail vein injection, once daily, for 5 days on and 2 days off. Tumor volume and body weight were measured twice weekly. Mice were euthanized when the tumor volume reached 20% of the body weight or the mice became moribund. Tumor response: For individual mice, progressive disease (PD) was defined as <50% regression from initial volume during the study period and >25% increase in initial volume at the end of study period. Stable disease (SD) was defined as <50% regression from initial volume during the study period and ≤25% increase in initial volume at the end of the study. Partial response (PR) was defined as a tumor volume regression of ≥50% for at least one time point but with measurable tumor (≥0.10 cm 3 ). Complete response (CR) was defined as a disappearance of measurable tumor mass (<0.10 cm 3 ) for at least one time point.

Pharmacodynamic studies
Tumor-bearing mice (SJNB11 and SJNB19 PDXs) were treated with indisulam (5, 12.5, and 25 mg/kg) for five consecutive days via tail injection. At the end of the experiment, mice were euthanized by CO 2 . Tumor and organ samples were excised and immediately snap-frozen in liquid nitrogen and preserved at −80°C until further use. Tissue lysates were made using 2× lysis buffer [1 M tris-HCl (pH 6.8), 10% SDS, glycerol, -mercaptoethanol, and bromophenol blue]. Briefly, a 5-mm 3 piece of tissue was added to 500 l of 2× lysis buffer, homogenized for 30 s (PRO200, Biogen, USA), heated at 95°C, and centrifuged at 12,000 rpm for 5 min at 4°C. Supernatant was further diluted 5×, sonicated for 30 s, and heated at 95°C. The sample was centrifuged at 13,000 rpm for 5 min and used as whole-tissue protein lysate.

Complete blood counting
Complete blood count was analyzed with the Forcyte analyzer hematology machine, according to the manufacturer's instructions. Briefly, blood was collected in Eppendorf tubes containing 10 l of 10% EDTA via retroorbital bleed using 200-l heparinized capillary tubes (catalog no. 22-362-566, Fisher brand). Blood samples were processed within 2 hours to avoid hemolysis. The number of leukocytes (white blood cell), erythrocytes (red blood cell), lymphocytes, neutrophils, monocytes, eosinophils, and platelets was counted. Proprietary lysing agent was added to liberate hemoglobin and ultimately convert it to cyanmethemoglobin to calculate the value.

Statistical analysis
Wilcoxon rank sum test (two-sided) was used to compare the volume between two groups at every time point. P values across multiple time points were adjusted for multiple comparison using the Benjamin-Hochberg method. Kaplan-Meier survival analyses were used with log-rank (Mantel-Cox) method in Prism program.

SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at https://science.org/doi/10.1126/ sciadv.abj5405 View/request a protocol for this paper from Bio-protocol.