Sunitinib efficacy with minimal toxicity in patient-derived retinoblastoma organoids
Journal of Experimental & Clinical Cancer Research volume 42, Article number: 39 (2023)
Recurrence of retinoblastoma (RB) following chemoreduction is common and is often managed with local (intra-arterial/intravitreal) chemotherapy. However, some tumors are resistant to even local administration of maximum feasible drug dosages, or effective tumor control and globe preservation may be achieved at the cost of vision loss due to drug-induced retinal toxicity. The aim of this study was to identify drugs with improved antitumor activity and more favorable retinal toxicity profiles via screening of potentially repurposable FDA-approved drugs in patient-derived tumor organoids.
Genomic profiling of five RB organoids and the corresponding parental tissues was performed. RB organoids were screened with 133 FDA-approved drugs, and candidate drugs were selected based on cytotoxicity and potency. RNA sequencing was conducted to generate a drug signature from RB organoids, and the effects of drugs on cell cycle progression and proliferative tumor cone restriction were examined. Drug toxicity was assessed with human embryonic stem cell-derived normal retinal organoids. The efficacy/toxicity profiles of candidate drugs were compared with those of drugs in clinical use.
RB organoids maintained the genomic features of the parental tumors. Sunitinib was identified as highly cytotoxic against both classical RB1-deficient and novel MYCN-amplified RB organoids and inhibited proliferation while inducing differentiation in RB. Sunitinib was a more effective suppressor of proliferative tumor cones in RB organoids and had lower toxicity in normal retinal organoids than either melphalan or topotecan.
The efficacy and retinal toxicity profiles of sunitinib suggest that it could potentially be repurposed for local chemotherapy of RB.
Retinoblastoma (RB), a tumor of the developing retina, is the most frequent primary intraocular cancer in children and has an estimated worldwide incidence of 8000 cases annually . The advent of chemotherapy, first systemically and then locally via intra-arterial or intravitreal administration, has led to improved treatment outcomes and great success in globe retention [2,3,4,5,6,7]. However, some tumors are resistant to even local administration of maximum feasible drug dosages, and surgical removal of the eye becomes necessary to prevent death associated with extraocular dissemination [2, 3, 6]. In other cases, tumor control is effective and the globe can be preserved, but at the cost of drug-induced retinal toxicity and loss of vision [8,9,10,11,12]. Thus, there is an urgent need to identify drugs with improved antitumor activity and more favorable retinal toxicity profiles for the treatment of RB.
Drug repurposing (also known as repositioning or reprofiling) is receiving increased interest as an alternative and more cost-effective strategy compared with de novo drug development. However, preclinical testing of drugs is still necessary to avoid unexpected lack of efficacy and/or unnecessary side effects during treatment for a new indication. Tumor organoids derived from patient tissues have been reported to faithfully model many aspects of the parent tumor and to serve as effective in vitro models for testing of anticancer agents [13,14,15,16,17]. In many studies, organoid cultures have been shown to recapitulate the genomic and transcriptomic profiles of donor tumors and be predictive of patient-drug responses [13, 16]. Additionally, differences in drug responses associated with intratumoral heterogeneity has been demonstrated in patient-derived organoids [14, 15, 17].
Tumor organoids representing two known subtypes of RB have been reported and shown to exhibit tumor-relevant genotypes and phenotypes [18, 19]. Classical RB is initiated by biallelic inactivation of the tumor suppressor gene RB1 and identified in most RB cases, whereas novel subtype of RB is initiated by high focal MYCN amplification and accounts for < 2% of all cases [20, 21]. RB1-deficient and MYCN-amplified RBs appear to differ with respect to genetic and epigenetic characteristics that lead to differences in expression levels of cone markers and in stemness [20, 22,23,24]. Given these differences, it is possible that the RB subtypes may display differing responses to the same chemotherapy regimen, and whether a single drug could be developed to treat both RB subtypes is currently unknown.
In this study we aimed to discover whether one or more FDA-approved anticancer drugs could potentially be repurposed for the treatment of RB. We performed high-throughput screening of 133 FDA-approved drugs on fully characterized patient-derived RB organoids and conducted detailed mechanistic investigations of selected candidates. We compared the efficacy and toxicity profiles of candidate drugs in RB1-deficient and MYCN-amplified RB organoids and in retinal organoids derived from human embryonic stem cells (hESCs) with those of clinically used anticancer drugs.
RB was located via transillumination of the enucleated eyes with a fiber-optic light source; a puncture was created in the sclera at the tumor shadow for tissue excision. A tumor biopsy with the total volume of 12–18 mm3 was obtained directly from the enucleated globes of patients. Fresh surgical specimens were collected and processed as described previously . Blood was drawn from patients for germline RB1 mutation screening. All experimental protocols were approved by the Institutional Review Board at the Faculty of Medicine Ramathibodi Hospital, Mahidol University (protocol number ID11-58-53 and ID12-61-32) and performed in accordance with the relevant guidelines and regulations. Informed consent was obtained from a parent of each patient before the samples were collected.
Four lines of organoids were previously established (RB668, RB654, RB170, and RB187) [18, 19]; the other line (RB394) was generated in this study from RB tissues. Derivation of RB organoids was performed as previously described . In brief, RB tissues were dissociated, and tumor cells were resuspended in organoid medium and embedded in growth factor-reduced Matrigel solution (Corning) at a ratio of 1:1.8 by volume. Aliquots of 20 µL/well of the mixed cell-gel solution were added to 6-well plates via 5–7 drops/well and allowed to solidify at 37 °C. Organoid medium was added to cover the gel drops, and cultures were maintained in a humidified incubator 5% CO2 incubator at 37 °C. Every 3–4 weeks, organoids were manually dissociated and passaged at a 1:3 or 1:4 ratio by embedding in fresh Matrigel.
The Approved Oncology Drugs Set VIII, consisting of 133 FDA-approved chemotherapeutic drugs was kindly provided by the National Cancer Institute Division of Cancer Treatment and Diagnosis/ Developmental Therapeutics Program (http://dtp.cancer.gov) of the National Institutes of Health (Bethesda, MD, USA). Drugs were provided at concentrations of 10 mM in 100% dimethyl sulfoxide (DMSO). Stock plates were diluted in 96- or 384-well polypropylene V-bottom microplates using an automated JANUS liquid handling platform (Perkin Elmer).
Primary drug screening and imaging-based drug response analysis
Tumor organoids were manually dissociated, incubated with TrypLE Express for 15 min at 37 °C, and resuspended in organoid medium with 5% Matrigel. Droplets of the cell/gel mixture (8 × 104 cells/10 µL) were placed in the center of wells of a 96-well plate. After a thin gel had formed, medium was added and the cultures were maintained for 7 days. RB organoids were then exposed to drugs (10 µM) or vehicle (0.1% DMSO) for 72 h. To evaluate the drug response, viable and dead cells were distinguished by incubating the organoids with a mixture of the fluorescent DNA-binding dyes Hoechst 33342 (Invitrogen) and ethidium homodimer-1 (Abcam, 1:1000) for 15 min. Organoids were imaged before and after drug treatment using an Operetta plate imaging system (Perkin Elmer). The area of live organoid cells (Hoechst 33342-positive, ethidium homodimer-1−negative) was determined using a Columbus image data storage and analysis system. Normalized growth rate inhibition (GR) was calculated as previously described  (see Supplementary Methods for details). GR values range from − 1 to + 1, where negative values indicated cytotoxicity, a value of 0 reflected complete cytostasis, and positive values indicated partial growth inhibition.
Secondary drug screening and dose-response curve analysis
Organoid cells were seeded in 384-well culture plates at 7.5 × 103 cells/15 µL/well in organoid medium containing 5% Matrigel. A 6-point, 10-fold dilution series ranging from 0 to 100 µM or 200 µM was generated and added to the cells in duplicate using a JANUS liquid handler workstation. After 72 h, cell viability was quantified using a CellTiter-Glo (ATP) luminescent assay according to the manufacturer’s instructions (Promega). Dose-response curves were fitted using nonlinear least squares regression with variable slope (four parameters); half-maximal inhibitory concentration (IC50), the fraction of viable cells at the highest drug concentration (Emax), and curve steepness (Hill Slope, HS) were computed using GraphPad Prism 8.
Differentiation of hESCs into retinal organoids
Retinal organoids were generated according to a previously published protocol . Briefly, the hESC cell line KhES-1 (RIKEN, Japan) was dissociated and rapidly aggregated in mTeSR1 medium (STEMCELL Technologies) containing 10 µM Y-27632 (Sigma-Aldrich) at 3 ×103 cells/well in 96-well ultra-low adhesion V bottom Lipidure-coated plates (S-BIO). Differentiation media were added at specific time points as previously described, without the addition of γ-secretase inhibitor (DAPT) . Retinal organoids were maintained for 65–70 days, examined for the expression of photoreceptor-associated genes, and used to evaluate retinal toxicity of the drugs.
Sunitinib, melphalan, topotecan, or vehicle (0.1% DMSO) was added at the indicated concentrations to organoid cultures and incubated for 24 h. Drug efficacy and toxicity were evaluated by immunostaining.
Cell cycle analysis
Cell cycle was analyzed using 647 EdU click proliferation kit (BD Biosciences) in accordance with the manufacturer’s instructions. Briefly, organoid cells were incubated with 10 µM EdU for 90 min. Cells were fixed, permeabilized, and stained with EdU detection cocktail for 30 min. Cell suspensions were washed and incubated with propidium iodide for 15 min. Stained cells were immediately analyzed by flow cytometry (BD Biosciences). The data were analyzed using FlowJo software (BD Biosciences).
Retinal and tumor organoids were fixed with 3.7% and 2% paraformaldehyde solution, respectively for 10 min, washed with phosphate-buffered saline (PBS) and incubated in 30% (w/v) sucrose overnight, then embedded in OCT compound and snap frozen. Cryosections (10 µm) of samples were mounted on SuperFrost Plus slides for immunostaining. Cryosections were washed with PBS and stained with primary (overnight) and secondary antibodies. The following antibodies were used for staining: RXR gamma (1:100, mouse, sc-365252, Santa Cruz Biotechnology), Ki67 (1:100, rabbit, RB1510P0, Thermo Scientific™ Lab Vision™), c-Jun (1:400, rabbit, 9165, Cell Signaling Technology), and cleaved caspase 3 (1:400, rabbit, 9661, Cell Signaling Technology). Secondary antibodies (1:500) included Alexa Fluor 488 goat anti-mouse IgG and Alexa Fluor 568 goat anti-rabbit IgG (Invitrogen). Nuclei were counterstained by 4’,6-diamidino-2-phenylindole (DAPI). Fluorescent images were acquired by confocal laser scanning microscopy, and Z-stacking was performed with NIS-Element AR (Nikon).
Tumor organoids were lysed using a lysis buffer containing 1x Halt protease and phosphatase inhibitor cocktail (Thermo Fisher Scientific). Briefly, 20 µg total protein was separated using 10% SDS-PAGE gels and transferred to PVDF membranes. The membranes were blocked with 5% non-fat dry milk solution and incubated overnight at 4°C with anti-c-Jun rabbit antibodies (1:1000, 9165, Cell Signaling Technology), followed by incubation with anti-rabbit secondary antibodies conjugated to HRP. The target proteins were detected using a chemiluminescent HRP detection reagent (Millipore). GAPDH on the same membranes was detected using anti-GAPDH HRP-conjugated antibodies. The signal intensity was analyzed using ImageLab software (Bio-Rad). The relative densitometry of c-Jun bands was expressed as a function of GAPDH band densitometry and reported.
Differences in response to drugs (IC50) of RB organoids were evaluated using the extra-sum of squares F-test. Cell cycle and death analyses were evaluated with an unpaired t-test. Multiple groups were compared using one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparisons test or Welch’s ANOVA (for unequal variance) followed by Tamhane’s T2 multiple comparisons test. Significance was considered at P < 0.05. All statistical data analysis was performed using Prism 8. Hierarchical clustering and visualization were performed using heatmap.2.
Full details of the genomic and RNA sequencing methods and analyses, reverse transcription-quantitative polymerase chain reaction (RT-qPCR), and cell death analysis are provided in Supplementary Methods.
Tumor organoids maintain the genomic landscapes of the parental RB tissues
Microarray data analysis of copy number variation showed that genomic alterations in the RB organoids were similar to those of the corresponding parental tissue from which they were derived (Fig. 1 A and Supplementary Fig. S1). The most common RB-associated alterations  identified in both tumor tissues and organoids were gains at chromosome 1q and/or 2p (RB668, RB170, RB187, and RB394), 16q loss (RB170, RB187, and RB394), and 6p gain (RB668). Unlike the other RB tissues/organoids, RB170, which was previously reported to have high-level focal MYCN amplification (90 copies) as the disease initiator , had chromosomal gain in 7q, 13q, and 14q. Gain or loss of smaller fragments was detected in each chromosome (Fig. 1 A and Supplementary Fig. S1). In contrast to the other four organoid/tissue pairs, RB654 exhibited minimal genomic alterations in tissue and organoid.
Of note, many subclonal mutations were further enriched in the organoids and were evident in chromosome 10q for RB688, chromosomes 2 and 16 for RB187, and chromosomes 1, 11q, 15, and 18 for RB394 (Supplementary Fig. S1). Furthermore, some lesions that were not detected in tumor tissues were identified in organoids (2p and 16q for RB668 and 8 and 9q for RB170), and conversely, some lesions previously detected in tumors were virtually absent in organoids, as shown for 7p and 11p of RB394 tissue and organoids (Fig. 1A and Supplementary Fig. S1). These findings imply that clonal evolution that results in tumor heterogeneity may occur in organoids, as previously described . Analysis of whole-exome sequencing data indicated that more than 70% of pathogenic variants and mutant genes in the parental tumors were also present in the corresponding organoids for four of the pairs (Fig. 1B, C), but only 50% were maintained in RB668 organoids, possibly due to deletion of 13q and subsequent induction of significant genomic changes (Fig. 1A–C).
The five organoids were derived from patients with both heritable and non-heritable RB, as indicated by germline RB1 mutation analysis (Fig. 1D). Four organoids (RB668, RB654, RB187 and RB394) represented RB1-deficient RB, and the other (RB170) was MYCN-amplified RB . Many genes commonly altered in RB were identified, including those associated with 1q, 2p, and 6p gains and 16q loss (Fig. 1D) . Additionally, mutations were identified in genes associated with cone signal circuitry that is oncogenically activated in RB ; namely, MYCN (RB668, RB170, RB187 and RB394), RXRG (RB170, RB187 and RB394) and MDM2 (RB170). Taken together, these data demonstrate that RB organoids recapitulated the genomic landscapes of the parental tissues. We selected three organoid lines representing RB initiated by MYCN amplification (RB170) and biallelic RB1 inactivation (RB654 with minimal lesions and RB668 with 13q deletion) (Fig. 1A–D) for the initial drug screening.
High-throughput screening discriminates between cytotoxic and cytostatic effects of anticancer drugs on RB organoids
Organoids were incubated for 72 h with vehicle (0.1% DMSO) or 133 FDA-approved compounds at 10 µM (Fig. 2A). Areas of live and dead organoid cells were quantified before and after drug exposure and used to calculate the normalized GR on a scale of − 1 to + 1, where negative and positive numbers represent cytotoxicity and partial growth inhibition, respectively, and 0 reflects complete cytostasis  (Fig. 2B). RB organoids exhibited heterogenous drug responses; ~78% of drugs were cytotoxic in RB170 organoids compared with 49% and 65% of drugs in RB668 and RB654 organoids, respectively (Fig. 2B). Among the drugs commonly used to treat RB, at least one organoid line was resistant to carboplatin, etoposide, and vincristine, whereas all three organoid lines were inhibited by melphalan and topotecan (Fig. 2C).
We then selected the drugs that were highly cytotoxic based on GR ≤ the mean GR (-0.84) in all organoid lines. We identified 11 drugs that were consistently highly cytotoxic in all three organoid lines; these included five antitumor antibiotics, two antifolates, two tyrosine kinase receptor inhibitors, an alkylating agent, and a poly (ADP-ribose) polymerase inhibitor (Fig. 2D). However, melphalan and topotecan appeared less cytotoxic in RB668 and RB654, respectively, as indicated by high GR compared with the mean GR and thus were excluded (Fig. 2C).
Sunitinib is a potent inhibitor of RB organoid growth
To evaluate the potency and efficacy of the 11 selected drugs, full dose-response matrices were generated for the three organoid lines (Supplementary Fig. S2A–C). Three drugs failed to achieve 50% cytotoxicity over the range of concentrations tested and were excluded based on low potency (Supplementary Fig. S2A–C). The dose-response curves of the remaining 8 drugs (Fig. 3A–H) were used to compute IC50, Emax, and HS, which are measures of the potency, efficacy, and the degree of homogeneity of the drug response, respectively  (Fig. 3I–K). The matrices for the 8 drugs showed positive correlations between the IC50 and Emax and negative correlations between HS and both IC50 and Emax, indicating that the drugs with high potency and efficacy induced a homogeneous response in the organoid cultures (Supplementary Fig. S2D).
We found that sunitinib, a multiple tyrosine kinase inhibitor, was the only drug that had an IC50 of < 10 µM in all three organoid lines [RB170 (MYCN amplification), RB654 (RB1-/-) and RB668 (RB1-/-)]; the drug responses were similar across the organoids (Fig. 3A, I). The responses to afatinib were also similar across the organoids, but the IC50 were > 10 µM in RB668 and RB654 (Fig. 3B, I). Additionally, sunitinib and afatinib had HS values > 1 in all three organoid lines, indicative of homogeneous responses (Fig. 3J). Sunitinib, afatinib, and epirubicin had low Emax values (< 10%) in the three organoid lines (Fig. 3K). Collectively, these data showed that sunitinib outperformed the other drugs tested in terms of potency (IC50), efficacy (Emax), and response (HS) in all three organoids. We also compared sunitinib with the clinically used drugs carboplatin, melphalan, and topotecan in five organoid lines (RB187, RB394, RB170, RB654, and RB668). Sunitinib had consistently superior potency and efficacy compared with carboplatin and melphalan in four lines (RB187, RB394, RB170, and RB654) (Fig. 3L–P). When compared with topotecan sunitinib showed higher potency in three lines (RB187, RB394, and RB668) and higher efficacy in all lines (Fig. 3L–P). Additionally, sunitinib treatment provided for the most homogeneity of the drug response as indicated by curve steepness resulting in high HS values (Fig. 3L–P).
Sunitinib signatures indicate the drug effects on proliferation and differentiation in RB
To examine the action of sunitinib in RB tumors RB654 (RB1-/-) and RB170 (MYCN amplification) organoids reflecting on the case study of two RB subtypes initiated with different cancer driver genes were selected for downstream analyses. Sunitinib signatures were generated by evaluating the differentially expressed genes (DEGs) in sunitinib-treated vs. vehicle-treated RB organoids. The concentration used was about IC10, which retained 90% cell viability, to avoid excessive cell death (Supplementary Fig.S3A) and thus reflect the changes in gene expression of sunitinib-treated tumor cells. The sunitinib signatures of RB170 and RB654 comprised 422 and 194 DEGs, respectively. Gene Ontology (GO) analysis revealed that the most significantly enriched terms for both organoids were strongly associated with mitogen-activated protein kinase (MAPK) signaling (Fig. 4A, B). This finding was supported by additional Kyoto Encyclopedia of Genes and Genomes pathway analysis (Supplementary Fig. S3B, C).
Next, we analyzed gene-biological concept networks of the top five GO terms for RB170 and showed interconnection between “regulation of the MAPK cascade” and other terms related to negative regulation of proliferation and positive regulation of organ morphogenesis (Fig. 4 C). JUN, a component of the AP-1 transcription factor that regulates cell proliferation and survival , and CD44, a target of AP-1 , were dramatically downregulated by sunitinib treatment (Fig. 4C), suggesting that sunitinib exerted antiproliferative on RB organoids via regulation of MAPK signaling. Additionally, NOTCH1 was connected to all five GO terms, and its downregulation is known to be associated with impairment of RB proliferation [30, 31].
The same analysis was performed for RB654 and showed that the term “inactivation of MAPK activity” was connected to three GO terms associated with “sensory organ/camera-type eye development” and “peptidyl-threonine dephosphorylation” (Fig. 4D). JUN was downregulated in RB654 by sunitinib treatment, whereas NRL, a transcription factor of rod photoreceptor fate  that is usually depleted in RB , was upregulated (Fig. 4D). Both JUN and NRL were included in the GO terms related to development, prompting us to speculate that sunitinib may promote tumor differentiation via inactivation of MAPK activity in RB654 organoids. Notably, JUN was the most significantly downregulated gene from tumor organoids of both RB subtypes (Fig. 4E, F). Downregulation of the JUN expression was confirmed by mRNA and protein analyses (Supplementary Fig. S3D and Fig. 4G, H). We found that c-Jun was co-expressed with RXRG cone marker; the levels of c-Jun were reduced in sunitinib-treated RB170 and RB654 cells (Fig. 4I, J). Examining retinal genes in a set of DEGs revealed that cone (RXRG and PDE6C) and amacrine (SLC6A9 also expressed in cones) genes in RB170 and rod gene (NRL) in RB654 were upregulated after sunitinib treatment (Fig. 4E, F). In contrast, bipolar (SYT2, CABP5) and rod (GNGT1)-enriched genes were downregulated in organoid cells after sunitinib treatment (Fig. 4E, F). Taken together, sunitinib signatures suggest that the drug affected cell proliferation and differentiation in RB.
Sunitinib suppresses RB proliferative tumor cones
To support the findings from sunitinib signatures the effect of drug treatment on cell cycle progression was determined. The cell cycle analysis indicated that treatment of RB170 and RB654 organoids for 24 or 72 h with sunitinib induced significant G1 arrest and a reduction in cell numbers in S and M phases (Fig. 5 A, B), which is in line with the GO term related to antiproliferation (Fig. 4A, C).
Most tumor cells in RB have characteristics of cone photoreceptors, including the subpopulations that have ability to self-renew and initiate new RB in vivo [33,34,35]. To determine whether drugs suppressed proliferative tumor cones, thus inhibiting tumor growth, we examined co-localization of the cone-specific transcription factor RXRG , which is expressed in tumor cells of origin in both RB1–/–- and MYCN-driven RB subtypes [33, 37,38,39], and the proliferative marker Ki67 in RB170 and RB654 organoids treated with vehicle or selected drugs at their IC50 for 24 h (Fig. 5 C). We compared the effects of sunitinib on proliferative tumor cones with melphalan and topotecan.
Immunostaining of organoids revealed that a significant reduction in the proportion of proliferative tumor cones (RXRG+ Ki67+: RXRG+) was observed after treatment of RB170 organoids with sunitinib (49.29% ± 2.02%, mean ± SEM) compared with vehicle (62.71% ±1.57%), topotecan (57.43% ±1.67%) and melphalan (65.46% ± 1.23%; Fig. 5D–H). In contrast, topotecan and melphalan had no significant effects on proliferative tumor cones compared with vehicle (Fig. 5D, F–H). There was a slight increase in the number of non-cycling cells (RXRG+) after treatment of RB170 with the tested drugs, but only sunitinib significantly decreased the number of cycling cells (RXRG+ Ki67+), resulting in a reduced proportion of proliferative tumor cones (Fig. 5H).
The same analysis in RB654 organoids showed significant reductions in the proportion of proliferative tumor cones (RXRG+ Ki67+: RXRG+) following treatment with sunitinib (42.74% ± 2.72%) and topotecan (51.92% ± 1.28%) but no significant change with melphalan treatment (57.07% ± 1.71%) compared with vehicle (59.03% ± 1.79%) (Fig. 5I–M). Interestingly, the non-cycling cells (RXRG+) increased, whereas the cycling cells (RXRG+ Ki67+) decreased after only sunitinib treatment (Fig. 5M), making sunitinib more effective than either topotecan or melphalan as also observed with RB170. The findings from immunostaining analysis agreed with the results from cell cycle and drug signature analyses and indicated that sunitinib inhibited proliferation as shown by the number of cycling and non-cycling cells in RB organoids. Altogether, sunitinib was superior to the other tested drugs in terms of suppression of proliferative tumor cones.
Sunitinib exhibits low toxicity in normal hESC-derived retinal organoids and has high efficacy against tumor organoids
To examine the potential retinal toxicity of sunitinib, we generated retinal organoids derived from hESCs, which represent a useful model of otherwise inaccessible human retina . Retinal organoids staining positively with the RXRG cone marker (day 65–70; Fig. 6A–H) were incubated with vehicle, sunitinib, topotecan, or melphalan [at the IC50 for tumor organoids for 24 h (Fig. 5C)], and cell death was examined by immunostaining of cleaved caspase 3 (CC3), a marker of apoptosis. Under the identical conditions used for evaluation of drug efficacy, the percentage CC3+ cells in retina organoids were 0.66% ± 0.11% (mean ± SEM) for vehicle; 6.57% ± 0.44% and 9.15% ± 0.79% for 4 µM and 6 µM sunitinib, respectively; 12.82% ± 0.78% and 21.91% ± 0.83% for 0.15 µM and 1 µM topotecan, respectively; and 13.71% ± 0.68% and 22.63% ± 0.70% for 7 µM and 17 µM melphalan, respectively (Fig. 6B–I). Topotecan and melphalan were clearly highly toxic to the retinal organoids compared with both vehicle and sunitinib despite treatment with lower drug dose (< IC50 for 7 µM melphalan) (Fig. 6B–I).
Immunostaining of CC3 was conducted to examine the tumor killing effect of the same set of drugs under the identical conditions used for evaluating drug toxicity. The percentage CC3+ cells in RB170 organoids were 4.59 ± 0.31% (mean ± SEM) for vehicle; 8.16% ± 0.28% for sunitinib; 6.99% ± 0.56% for topotecan; and 4.47% ± 0.31% for melphalan (Fig. 6J–N). In RB654, the percentage CC3+ cells were 7.28 ± 0.29% for vehicle; 17.27% ± 0.59% for sunitinib; 15.11% ± 0.50% for topotecan; and 12.74% ± 0.68% for melphalan (Fig. 6O–S). Sunitinib was more effective than other drugs in killing RB tumor cells (Fig. 6N, S). Death cells from sunitinib treatment were more in tumor organoids but fewer in retinal organoids than those from topotecan and melphalan treatment (Fig. 6I, N, S). Thus, sunitinib displayed high efficacy against RB with minimal retinal toxicity compared with topotecan and melphalan.
A drug with improved antitumor activity and a favorable retinal safety profile is urgently needed to advance globe-sparing therapies in RB. In the present study, we sought to identify FDA-approved drugs that could potentially be repurposed for RB treatment. We performed high-throughput screening of a panel of 133 drugs and successfully identified sunitinib as a highly cytotoxic agent for RB. Sunitinib effectively suppressed proliferative tumor cones in RB organoids but had minimal toxicity towards normal retinal organoids, and outperformed melphalan and topotecan, which are used for local chemotherapy of RB.
RB organoids derived from both RB1-deficient and MYCN-amplified subtypes were used as a model for drug screening. MYCN-amplified organoids were generally more susceptible than RB1-deficient organoids to inhibition by the panel of anticancer drugs. High MYCN expression confers a rapid growth and proliferation phenotype [41, 42], which may explain this sensitivity to drugs that typically kill rapidly dividing cells. In line with this possibility, poorly differentiated cells in MYCN-amplified RB have been shown to be more responsive to chemotherapy than differentiated tumor cells in RB1-deficient RB with featured rosettes [4, 22]. Nevertheless, we showed that sunitinib exhibited anticancer activity against RB organoids of both subtypes in the present study.
Sunitinib has shown broad and potent antitumor and/or antiangiogenic activity in preclinical studies of multiple cancer types [43,44,45]. We demonstrated the antitumor activity of sunitinib in RB, consistent with the previous report on direct effect of sunitinib on inhibition of RB invasion and metastasis , although the effect on angiogenesis is inconclusive in orthotopic model of RB in zebrafish . Sunitinib has been shown to induce G1-phase arrest in acute myeloid leukemia (AML) cells , similar to our observations here with RB, but this contrasts with the G2/M-phase arrest seen in sunitinib-treated glioblastoma cultures .
We speculate that sunitinib affects tumor proliferation and differentiation via downregulation of JUN for suppression of the RB growth. c-Jun is overexpressed in a large number of cancers , and its silencing has been shown to inhibit proliferation, migration, and invasion of cancer cells [48, 49]. Downregulation of JUN with a concomitant G1-phase arrest is consistent with the function of c-Jun in promoting G0 to G1 transition [49, 50] and indicates the antiproliferative effect of sunitinib in tumor organoids of both RB subtypes. Furthermore, high expression levels of c-Jun are linked to poorly differentiated cells in a high grade of malignancy in glial tumors , suggesting the role of c-Jun in inhibition of cancer differentiation. Downregulation of JUN with concomitant changes in expression of photoreceptor genes indicates sunitinib-induced RB differentiation in MYCN-amplified and RB1-deficient organoids. As indicated by the upregulation of cone-specific transcription factor RXRG , sunitinib induced cone differentiation in MYCN-amplified RB170 which has been classified as aggressive RB exhibiting poor differentiation and dedifferentiated cone states (subtype 2) . RB1-deficient RB654, classified as differentiated tumors with high expression of cone genes (subtype 1) , underwent further differentiation by sunitinib treatment as shown by the upregulation of a critical transcription factor NRL which is essential for rod differentiation . An additional supporting evidence for sunitinib-induced RB differentiation is an increase in non-cycling cells and decrease in cycling cells that result in a reduced proportion of proliferative tumor cones in RB organoids. RB that is turned into more differentiated states might be less aggressive by sunitinib treatment, as previously described in human AML blast cells undergoing sunitinib-elicited monocytic differentiation , providing a potential strategy for the cancer treatment.
Intravitreal or intra-arterial chemotherapy has increased the rate of globe retention in globe-preservation treatment and significantly reduces drug-induced systemic toxicity [5, 10]. However, drug-induced retinal toxicity remains a major drawback of local chemotherapy. Melphalan and topotecan have been reported to cause retinal toxicity, as measured by reductions in electroretinography [8,9,10,11,12], particularly when administered intravitreally alone or in combination with intra-arterial chemotherapy for the management of vitreous disease in RB. This means that the choice of chemotherapeutic agents is crucial to advance the treatment. Our study suggests that sunitinib has potential utility for local chemotherapy for RB given its high efficacy in terms of suppressing proliferative tumor cones and minimal retinal toxicity compared with other drugs currently in use.
Patient-derived tumor organoids and hESC-derived retinal organoids as a 3-dimensional model are of use to drug discovery. Additionally, previous work has shown that RB-like retinal organoids can be generated from RB1-/- hESC and iPSC cells and allow for the possibility of testing therapeutic agents [38, 39]. However, in vivo models may be needed to verify antineoplastic effects and retinal safety of drugs. Our study provides the groundwork for and encourages further investigation of the therapeutic potential of sunitinib for local chemotherapy of RB. In a preclinical model, sunitinib has been shown to have a good retinal safety profile when administered by intravitreal injection , which supports the potential utility of locally delivered sunitinib for treating RB.
In the present study, we utilized RB organoids derived from different RB subtypes to screen for drugs that could be repurposed for RB. We identified sunitinib as a possible drug for RB. Sunitinib exhibits high efficacy with minimal retinal toxicity, suggesting that it may be suitable for local therapy. Sunitinib is FDA-approved for the treatment of renal carcinoma and imatinib-refractory gastrointestinal stromal tumors, and our study identifies RB as a potential new indication.
Availability of data and materials
All data generated or analyzed during this study are included in this published.
article (and its supplementary information files). The datasets generated during and/or analyzed during the current study are available in the Gene Expression Omnibus (GEO) repository, [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE201740]
- IC50 :
Half-maximal inhibitory concentration
- Emax :
The fraction of viable cells at the highest drug concentration
Growth rate inhibition
Human embryonic stem cells
Differentially expressed genes
- RB1 :
RB transcriptional corepressor 1
- MYCN :
MYCN proto-oncogene, BHLH transcription factor
- RXRG :
Retinoid X receptor gamma
- NR L :
Neural retina leucine zipper
- MDM2 :
- JUN :
Jun proto-oncogene, AP-1 transcription factor subunit
Marker of proliferation Ki-67
Chantada GL, Qaddoumi I, Canturk S, Khetan V, Ma Z, Kimani K, et al. Strategies to manage retinoblastoma in developing countries. Pediatr Blood Cancer. 2011;56(3):341–8.
Abramson DH, Daniels AB, Marr BP, Francis JH, Brodie SE, Dunkel IJ, et al. Intra-arterial chemotherapy (ophthalmic artery chemosurgery) for Group D retinoblastoma. PLoS One. 2016;11(1):e0146582.
Francis JH, Roosipu N, Levin AM, Brodie SE, Dunkel IJ, Gobin YP, et al. Current treatment of bilateral retinoblastoma: the impact of intraarterial and intravitreous chemotherapy. Neoplasia. 2018;20(8):757–63.
Kaewkhaw R, Rojanaporn D, Retinoblastoma. Etiology, modeling, and treatment. Cancers (Basel). 2020;12(8):2304.
Munier FL, Beck-Popovic M, Chantada GL, Cobrinik D, Kivela TT, Lohmann D, et al. Conservative management of retinoblastoma: challenging orthodoxy without compromising the state of metastatic grace. “Alive, with good vision and no comorbidity”. Prog Retin Eye Res. 2019;73:100764.
Rojanaporn D, Chanthanaphak E, Boonyaopas R, Sujirakul T, Hongeng S, Ayudhaya SSN. Intra-arterial chemotherapy for retinoblastoma: 8-year experience from a tertiary referral institute in Thailand. Asia Pac J Ophthalmol (Phila). 2019;8(3):211–7.
Shields CL, Mashayekhi A, Au AK, Czyz C, Leahey A, Meadows AT, et al. The International classification of Retinoblastoma predicts chemoreduction success. Ophthalmology. 2006;113(12):2276–80.
Francis JH, Brodie SE, Marr B, Zabor EC, Mondesire-Crump I, Abramson DH. Efficacy and toxicity of intravitreous chemotherapy for retinoblastoma: four-year experience. Ophthalmology. 2017;124(4):488–95.
Francis JH, Iyer S, Gobin YP, Brodie SE, Abramson DH. Retinoblastoma vitreous seed clouds (class 3): a comparison of treatment with ophthalmic artery chemosurgery with or without intravitreous and periocular chemotherapy. Ophthalmology. 2017;124(10):1548–55.
Francis JH, Schaiquevich P, Buitrago E, Del Sole MJ, Zapata G, Croxatto JO, et al. Local and systemic toxicity of intravitreal melphalan for vitreous seeding in retinoblastoma: a preclinical and clinical study. Ophthalmology. 2014;121(9):1810–7.
Liao A, Hsieh T, Francis JH, Lavery JA, Mauguen A, Brodie SE, et al. Toxicity and efficacy of intravitreal melphalan for retinoblastoma: 25 µg versus 30 µg. Retina. 2021;41(1):208–12.
Yousef YA, Al Jboor M, Mohammad M, Mehyar M, Toro MD, Nazzal R, et al. Safety and efficacy of intravitreal chemotherapy (Melphalan) to treat vitreous seeds in retinoblastoma. Front Pharmacol. 2021;12:696787.
Arena S, Corti G, Durinikova E, Montone M, Reilly NM, Russo M, et al. A subset of colorectal cancers with cross-sensitivity to olaparib and oxaliplatin. Clin Cancer Res. 2020;26(6):1372–84.
de Witte CJ, Espejo Valle-Inclan J, Hami N, Lohmussaar K, Kopper O, Vreuls CPH, et al. Patient-derived ovarian cancer organoids mimic clinical response and exhibit heterogeneous inter- and intrapatient drug responses. Cell Rep. 2020;31(11):107762.
Lee SH, Hu W, Matulay JT, Silva MV, Owczarek TB, Kim K, et al. Tumor evolution and drug response in patient-derived organoid models of bladder cancer. Cell. 2018;173(2):515–28e17.
Ooft SN, Weeber F, Dijkstra KK, McLean CM, Kaing S, van Werkhoven E, et al. Patient-derived organoids can predict response to chemotherapy in metastatic colorectal cancer patients. Sci Transl Med. 2019;11(513):eaay2574.
Schumacher D, Andrieux G, Boehnke K, Keil M, Silvestri A, Silvestrov M, et al. Heterogeneous pathway activation and drug response modelled in colorectal-tumor-derived 3D cultures. PLoS Genet. 2019;15(3):e1008076.
Saengwimol D, Chittavanich P, Laosillapacharoen N, Srimongkol A, Chaitankar V, Rojanaporn D, et al. Silencing of the long noncoding RNA MYCNOS1 suppresses activity of MYCN-amplified retinoblastoma without RB1 mutation. Invest Ophthalmol Vis Sci. 2020;61(14):8.
Saengwimol D, Rojanaporn D, Chaitankar V, Chittavanich P, Aroonroch R, Boontawon T, et al. A three-dimensional organoid model recapitulates tumorigenic aspects and drug responses of advanced human retinoblastoma. Sci Rep. 2018;8(1):15664.
Rushlow DE, Mol BM, Kennett JY, Yee S, Pajovic S, Theriault BL, et al. Characterisation of retinoblastomas without RB1 mutations: genomic, gene expression, and clinical studies. Lancet Oncol. 2013;14(4):327–34.
Zugbi S, Ganiewich D, Bhattacharyya A, Aschero R, Ottaviani D, Sampor C, et al. Clinical, genomic, and pharmacological study of MYCN-amplified RB1 wild-type metastatic retinoblastoma. Cancers (Basel). 2020;12(9):2714.
Kooi IE, Mol BM, Moll AC, van der Valk P, de Jong MC, de Graaf P, et al. Loss of photoreceptorness and gain of genomic alterations in retinoblastoma reveal tumor progression. EBioMedicine. 2015;2(7):660–70.
Liu J, Ottaviani D, Sefta M, Desbrousses C, Chapeaublanc E, Aschero R, et al. A high-risk retinoblastoma subtype with stemness features, dedifferentiated cone states and neuronal/ganglion cell gene expression. Nat Commun. 2021;12(1):5578.
Roohollahi K, de Jong Y, van Mil SE, Fabius AWM, Moll AC, Dorsman JC. High-level MYCN-amplified RB1-proficient retinoblastoma tumors retain distinct molecular signatures. Ophthalmol Sci. 2022;2(3):100188.
Hafner M, Niepel M, Chung M, Sorger PK. Growth rate inhibition metrics correct for confounders in measuring sensitivity to cancer drugs. Nat Methods. 2016;13(6):521–7.
Eldred KC, Hadyniak SE, Hussey KA, Brenerman B, Zhang PW, Chamling X, et al. Thyroid hormone signaling specifies cone subtypes in human retinal organoids. Science. 2018;362(6411):eaau6348.
Fallahi-Sichani M, Honarnejad S, Heiser LM, Gray JW, Sorger PK. Metrics other than potency reveal systematic variation in responses to cancer drugs. Nat Chem Biol. 2013;9(11):708–14.
Shaulian E, Karin M. AP-1 in cell proliferation and survival. Oncogene. 2001;20(19):2390–400.
Lamb RF, Hennigan RF, Turnbull K, Katsanakis KD, MacKenzie ED, Birnie GD, et al. AP-1-mediated invasion requires increased expression of the hyaluronan receptor CD44. Mol Cell Biol. 1997;17(2):963–76.
Asnaghi L, Tripathy A, Yang Q, Kaur H, Hanaford A, Yu W, et al. Targeting notch signaling as a novel therapy for retinoblastoma. Oncotarget. 2016;7(43):70028–44.
Xiao W, Chen X, He M. Inhibition of the Jagged/Notch pathway inhibits retinoblastoma cell proliferation via suppressing the PI3K/Akt, src, p38MAPK and Wnt/β–catenin signaling pathways. Mol Med Rep. 2014;10(1):453–8.
Mears AJ, Kondo M, Swain PK, Takada Y, Bush RA, Saunders TL, et al. Nrl is required for rod photoreceptor development. Nat Genet. 2001;29(4):447–52.
Singh HP, Wang S, Stachelek K, Lee S, Reid MW, Thornton ME, et al. Developmental stage-specific proliferation and retinoblastoma genesis in RB-deficient human but not mouse cone precursors. Proc Natl Acad Sci U S A. 2018;115(40):E9391.
Xu XL, Fang Y, Lee TC, Forrest D, Gregory-Evans C, Almeida D, et al. Retinoblastoma has properties of a cone precursor tumor and depends upon cone-specific MDM2 signaling. Cell. 2009;137(6):1018–31.
Xu XL, Singh HP, Wang L, Qi DL, Poulos BK, Abramson DH, et al. Rb suppresses human cone-precursor-derived retinoblastoma tumours. Nature. 2014;514(7522):385–8.
Roberts MR, Hendrickson A, McGuire CR, Reh TA. Retinoid X receptor (gamma) is necessary to establish the S-opsin gradient in cone photoreceptors of the developing mouse retina. Invest Ophthalmol Vis Sci. 2005;46(8):2897–904.
Singh HP, Shayler DWH, Fernandez GE, Thornton ME, Craft CM, Grubbs BH, et al. An immature, dedifferentiated, and lineage-deconstrained cone precursor origin of N-Myc-initiated retinoblastoma. Proc Natl Acad Sci U S A. 2022;119(28):e2200721119.
Kanber D, Woestefeld J, Döpper H, Bozet M, Brenzel A, Altmüller J et al. RB1-negative retinal organoids display proliferation of cone photoreceptors and loss of retinal differentiation.Cancers (Basel). 2022;14(9).
Rozanska A, Cerna-Chavez R, Queen R, Collin J, Zerti D, Dorgau B, et al. pRB-depleted pluripotent stem cell retinal organoids recapitulate cell state transitions of retinoblastoma development and suggest an important role for pRB in retinal cell differentiation. Stem Cells Transl Med. 2022;11(4):415–33.
Kaewkhaw R, Kaya KD, Brooks M, Homma K, Zou J, Chaitankar V, et al. Transcriptome dynamics of developing photoreceptors in three-dimensional retina cultures recapitulates temporal sequence of human cone and rod differentiation revealing cell surface markers and gene networks. Stem Cells. 2015;33(12):3504–18.
Boon K, Caron HN, van Asperen R, Valentijn L, Hermus MC, van Sluis P, et al. N-myc enhances the expression of a large set of genes functioning in ribosome biogenesis and protein synthesis. EMBO J. 2001;20(6):1383–93.
Shostak A, Ruppert B, Ha N, Bruns P, Toprak UH, Eils R, et al. MYC/MIZ1-dependent gene repression inversely coordinates the circadian clock with cell cycle and proliferation. Nat Commun. 2016;7:11807.
Chow LQ, Eckhardt SG. Sunitinib: from rational design to clinical efficacy. J Clin Oncol. 2007;25(7):884–96.
de Bouard S, Herlin P, Christensen JG, Lemoisson E, Gauduchon P, Raymond E, et al. Antiangiogenic and anti-invasive effects of sunitinib on experimental human glioblastoma. Neuro Oncol. 2007;9(4):412–23.
Teng CL, Yu CT, Hwang WL, Tsai JR, Liu HC, Hwang GY, et al. Effector mechanisms of sunitinib-induced G1 cell cycle arrest, differentiation, and apoptosis in human acute myeloid leukaemia HL60 and KG-1 cells. Ann Hematol. 2013;92(3):301–13.
Chen X, Wang J, Cao Z, Hosaka K, Jensen L, Yang H, et al. Invasiveness and metastasis of retinoblastoma in an orthotopic zebrafish tumor model. Sci Rep. 2015;5:10351.
Brennan A, Leech JT, Kad NM, Mason JM. Selective antagonism of c-Jun for cancer therapy. J Exp Clin Cancer Res. 2020;39(1):184.
Lin G, Yu B, Liang Z, Li L, Qu S, Chen K, et al. Silencing of c-Jun decreases cell migration, invasion, and EMT in radioresistant human nasopharyngeal carcinoma cell line CNE-2R. Onco Targets Ther. 2018;11:3805–15.
Smith MJ, Prochownik EV. Inhibition of c-Jun causes reversible proliferative arrest and withdrawal from the cell cycle. Blood. 1992;79(8):2107–15.
Mechta-Grigoriou F, Gerald D, Yaniv M. The mammalian Jun proteins: redundancy and specificity. Oncogene. 2001;20(19):2378–89.
Blau L, Knirsh R, Ben-Dror I, Oren S, Kuphal S, Hau P, et al. Aberrant expression of c-Jun in glioblastoma by internal ribosome entry site (IRES)-mediated translational activation. Proc Natl Acad Sci U S A. 2012;109(42):E2875–84.
Tsujinaka H, Fu J, Shen J, Yu Y, Hafiz Z, Kays J, et al. Sustained treatment of retinal vascular diseases with self-aggregating sunitinib microparticles. Nat Commun. 2020;11(1):694.
We gratefully acknowledge the help and assistance of the staffs of the Excellent Center for Drug Discovery.
This work was supported by Mahidol University (136/2562 to RK), Health Systems Research Institute (HSRI 63 − 072 to RK), and Thailand Research Fund (TRF) and Office of the Higher Education Commission (OHEC) (MRG6280083 to RK)
Ethics approval and consent to participate
All experimental protocols were approved by the Institutional Review Board at the Faculty of Medicine Ramathibodi Hospital, Mahidol University (protocol number ID11-58-53 and ID12-61-32) and performed in accordance with the relevant guidelines and regulations. Informed consent was obtained from a parent of each patient before the samples were collected.
Consent for publication
All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Methods. for genomic and RNA sequencing analyses, cell death analysis, and RT-qPCR. Supplementary Figure S1. Copy number alterations and loss of heterozygosity by using the comparative genomic hybridization/single nucleotide polymorphism array. B-allelic frequency (BAF) and log R-ratio (LRR) plots indicate genomic alterations in tumor tissue (T) and the corresponding organoid (O). Supplementary Figure S2. Secondary drug screening assay: A–C Dose-response curves of 11 drugs tested in RB organoids. D Correlation coefficients of dose-response matrixes of 8 drugs tested in three organoid lines were computed by Pearson’s correlation at a significant level of 0.05. *P = 0.0337 for IC50 and Emax; * P = 0.0352 for Emax and HS; *** P = 4.885 × 10− 15 for Emax and area under the curve (AUC). Supplementary Figure S3. Cell death and pathway analyses: A Cell death analysis by annexin V-FITC and PI staining of RB170 and RB654 organoid cells. B, C KEGG pathway analysis of gene signature of sunitinib treatment in RB170 (B) and RB654 (C) organoids. D JUN mRNA expression in organoids treated with vehicle or sunitinib at IC10 and IC50 for 72 h.
About this article
Cite this article
Srimongkol, A., Laosillapacharoen, N., Saengwimol, D. et al. Sunitinib efficacy with minimal toxicity in patient-derived retinoblastoma organoids. J Exp Clin Cancer Res 42, 39 (2023). https://doi.org/10.1186/s13046-023-02608-1
- Drug screening
- Retinal toxicity