Mutagenesis of the RB1 Gene in hESCs.

To date, more than 3,000 RB1 mutations have been reported in the online RB1-LOVD database (23). The nonsense mutation p.R320X (c.958C > T) is the mutation most frequently associated with predisposition to Rb (SI Appendix, Fig. S1A). Thus, we generated p.R320X mutation knock-in hESCs using CRISPR/Cas9 genome-editing technology (Fig. 1 A and B and SI Appendix, Fig. S1 B–D). Three clones with homozygous mutation (RB1^Mut/Mut) in the RB1 gene were successfully established (SI Appendix, Fig. S1 E and F).

• Download figure
• Open in new tab
• Download powerpoint

Fig. 1.

Induction of 3D hRBOs from RB1^Mut/Mut and RB1^−/− hESCs generated by the CRISPR-Cas9 system. (A) Scheme of the stepwise development of hRBOs. (B) Schematic view of the targeting strategy for generation of RB1^Mut/Mut hESCs. (C) Targeting strategy for generation of RB1^−/− hESCs. (D) Representative images of developing hRBOs in different organoid stages. (Insets) Images of the whole organoids. Wild-type hESC-derived hROs served as controls. Images in the top row show a schematic overview of developing hRBOs. Asterisks indicate the abnormal proliferating “primary foci” of the tumor. Dashed lines denote the tumor-like structures. (Scale bar: 400 µm.) (E) Representative images show the rapid expansion of tumor-like structures (dashed lines) in developing hRBOs. (Scale bar: 2,000 µm.) (F) Quantification of RB1^Mut/Mut hESC-derived organoids with tumor-like structures in different organoid stages. Data are mean ± SD. n = 5. (G) Imaging of superior 3D structures of hRBOs and hROs at day 90. (i) Microscope view of screening sites in organoids with extension line. (Scale bar: 500 µm.) (ii) Imaging of screening sites in 3D rotation. (iii) Imaging of screening sites at high x-axis. Dashed lines denote Flexner–Wintersteiner rosettes. See also Movies S2 and S3. (H) TEM image of cells from day 105 hRBOs. Dashed lines denote fleurettes (i) and Flexner–Wintersteiner rosettes (ii); the pentagram marks the rosette cavity. Mitotic figures (arrowheads) represent proliferating cells, which are also represented by invagination in the nucleus (arrow). (iii) Asterisks indicate abundant mitochondria with heterogeneous morphology. (Scale bars: i, 10 µm; ii, 5 µm; iii, 2 µm.)

To further investigate the consequence of RB1 mutagenesis, we created two clones with a biallelic knockout (RB1^−/−) by targeting the first exon of the RB1 gene (Fig. 1C and SI Appendix, Fig. S1F). Western blot assays showed significantly reduced pRb expression in all lines (SI Appendix, Fig. S1G), which was validated by immunocytochemistry (SI Appendix, Fig. S1H).

The RB protein family (RB1, p107, and p130) plays a critical role in maintaining hESC self-renewal (24). To investigate the effect of RB1 inactivation on hESCs, we examined the pluripotency and cell cycle properties of hESCs using flow cytometry. In RB1-mutated clones, neither the proportion of OCT3/4⁺ cells (SI Appendix, Fig. S2A) nor the cell cycle (SI Appendix, Fig. S2 B and C) were altered. Furthermore, RNA sequencing (RNA-Seq) analysis showed that the expression of pluripotency, self-renewal, and cell cycle-associated genes remained unaffected by the mutagenesis of RB1 in hESCs (SI Appendix, Fig. S2D). Karyotyping of these genetically engineered hESCs showed a normal 46,XX karyotype (SI Appendix, Fig. S2E). Taken together, these results demonstrate that the RB1^Mut/Mut and RB1^−/− hESC lines sustained the primordial state without changing pluripotency and the cell cycle of the hESC lines.

Tumorigenesis in Retinal Organoids Derived from Mutant hESC Lines.

Inspired by the speculation that Rb originates from developing retinal neurons (2, 12, 25), we attempted to investigate Rb tumorigenesis through differentiation of retinal organoids from mutant hESC lines. We generated human retinal organoids using stepwise induction (Fig. 1A and SI Appendix, Fig. S3A) as described previously (16, 17, 19). As expected, the molecular and morphogenetic properties of the organoids recapitulated the developing human retina (SI Appendix, Fig. S3 B–D) (26, 27). There was no difference between hRBOs and hROs before differentiation day 45 (Fig. 1D and SI Appendix, Fig. S3D). Remarkably, at day 75, tumor-like “primary foci” were clearly visible in hRBOs derived from both RB1^−/− and RB1^Mut/Mut hESCs (Fig. 1D; SI Appendix, Fig. S3F; and Movie S1). In hRBOs, tumor-like structures exhibited obvious uneven interior density, ill-defined edges, and larger size (Fig. 1D and SI Appendix, Fig. S3F). Of note, the tumor-like foci expanded rapidly from the masses thereafter (Fig. 1 E and F and SI Appendix, Fig. S3E).

To further characterize hRBOs, we examined surface topography and internal structure in live organoids with deep optical coherence tomography. As shown in Fig. 1G, ill-defined edges, uneven reflectance, and notable inner cavities resembling Flexner–Wintersteiner rosettes were observed (Movies S2 and S3). Furthermore, transmission electron microscopy (TEM) demonstrated the same ultrastructural features as seen in primary Rb tumors, including a high nucleocytoplasmic ratio, abundant mitochondria with heterogeneous morphology, and tightly packed clusters of cells and rosettes (including Flexner–Wintersteiner rosettes characteristic of Rb and fleurettes characteristic of retinoma), as well as the presence of condensed chromosomes (mitotic figures) (Fig. 1H and SI Appendix, Fig. S3G). Taken together, these results demonstrate that tumorigenesis occurred in hRBOs, indicating that organoid Rbs differentiated from each mutated hESC line could be used for disease modeling.

Rb-Related Molecular Signatures in Induced hRBOs.

To further elucidate the molecular characteristics of hRBOs, we analyzed the gene expression and DNA methylation profiles of hRBOs by RNA-Seq and whole-genome bisulfite sequencing (WGBS). The transcriptomic profiles of hRBOs at different stages clearly demonstrated that proto-oncogenes and Rb signature genes (MKI67, spleen tyrosine kinase [SYK], proliferating cell nuclear antigen [PCNA], DEK, and TOP2A) (28⇓–30) were markedly up-regulated (Fig. 2A). In contrast, the tumor-suppressor genes or low-expressed genes (S100B, VIM, and LATS2) in primary Rb tumors were significantly down-regulated compared with normal hROs (Fig. 2A). Notably, transcriptomic dysregulation was more significant after day 75 (Fig. 2A).

• Download figure
• Open in new tab
• Download powerpoint

Fig. 2.

Molecular characterization of hRBOs. (A) Heatmaps showing the average log2-fold difference in the indicated genes selected from previously published transcriptome data between the hRBOs and normal hROs at different organoid stages. Representative signature genes (including the known oncogenes or tumor-suppressor genes) of human primary Rbs are shown on the left. Expression data of the indicated genes between the Y79 cell line and normal hROs older than 75 d are also shown in heatmaps. Rb, published transcriptomes of human primary Rbs (28). The top row shows correlation heatmaps of hRBOs at different organoid stages and human primary Rbs based on the expression of indicated genes. The larger the colored fan area, the greater the correlation. (B) Heatmaps of the selected top 10 KEGG pathway enrichments for differentially expressed genes (adjusted log2 fold change >1) between hRBOs and hROs at different organoid stages. Each row represents the enriched pathway terms. (C) Volcano plot visualization of significant differential methylated genes of hRBOs compared to hROs at day 120. Each point represents a distinct gene, and interesting genes are indicated. The lines represent average methylation differences of 0.25 and adjusted P value < 0.05, respectively. (D) Genomic snapshots of RNA-Seq and WGBS for SYK in hRBOs and parental hROs. A representative methylation level of SYK from human primary Rb sample (9) is also shown. An enlarged view of the specific site is shown below. (E) Immunostaining of Rb signature genes, Ki67, SYK, DEK, NSE, and p16^INK4a, in hRBOs and hROs at day 105. (Insets) Magnified views of the boxed region. Arrows indicate Ki67⁺ tumor cells with weak cytoplasmic staining for p16^INK4a. (Scale bar: 50 µm.)

When we compared the transcriptomes of hRBOs, primary Rb tumors (28), and the Y79 cell line, we found that 1,014 differentially expressed genes (574 up-regulated and 440 down-regulated) in hRBOs (after differentiation for >75 d) were similar to those reported in primary Rb tumors (Fig. 2A). In the gene expression profile, the hRBOs were more closely correlated with human primary Rb tumors than the Y79 cell line, and this correlation markedly increased across the developmental stages (Fig. 2A), consistent with the results of principal component analysis (SI Appendix, Fig. S4 A and B). A global comparative transcriptomic analysis of all genes between hROs and hRBOs revealed that the number of genes changed by more than twofold and gradually increased with the differentiation stage (SI Appendix, Fig. S4 C and D). Interestingly, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed dysregulation of key cancer-related signaling pathways involved in human Rb tumorigenesis, including the cell cycle and PI3K-Akt pathways (Fig. 2B and SI Appendix, Fig. S4E).

To further understand transcriptional regulation in Rb tumorigenesis, we focused on alterations in DNA methylation, an important epigenetic modification that contributes to transcription dysregulation in tumors. A total of 8,968 differentially methylated regions (DMRs) (4,597 hypermethylation and 4,371 hypomethylation) were identified (Fig. 2C). As expected, the DNA methylation states of key Rb-related genes, such as SYK, RXRG, MKI67, CCNE2, VIM, and CDH11, were significantly changed in hRBOs compared with control hROs (average methylation difference >0.25; adjusted P < 0.05) (Fig. 2C). Of note, SYK in the PI3K-Akt pathway, a pivotal marker gene, was identified as one of the most up-regulated genes during Rb tumorigenesis in hRBOs (Fig. 2A). The methylation level was markedly decreased at the CpG island shore located in the SYK promoter region in hRBOs (average methylation difference = 0.504’ adjusted P = 1.13e-07) and human primary Rb tumors (9) compared with control hROs (Fig. 2D), while the methylation level of VIM, a non–Rb-related gene, was significantly increased (SI Appendix, Fig. S5A), in line with the expression change. Moreover, we found that the DMRs were distributed along the whole genome and enriched in the noncoding regions compared with the coding regions (SI Appendix, Fig. S5 B and C). We also inspected the colocalization between DMR and CGI elements (CGI, CGI shore, and CGI shelf). Most (77%) of the DMRs were distal (>4 kb) to CGI elements, and more than one-half of CGI-proximal DMRs were significantly enriched in the CGI shore regions (χ² test, P = 2.9e-222) (SI Appendix, Fig. S5D). KEGG pathway analysis identified DMRs closely related to Axon guidance, Pathways in cancer, and the PI3K-Akt signaling pathway (SI Appendix, Fig. S5E).

Immunostaining of organoid sections further confirmed that pRb was readily detected in control hROs, but not in derived hRBOs (SI Appendix, Fig. S5F). In line with the findings from transcriptomic analysis, >70% of cells in hRBOs were highly proliferative (Ki67⁺), even after week 15 (Fig. 2E and SI Appendix, Fig. S5 F–H), while few cells in normal organoids continued to divide. p16^INK4a (CDKN2A), which has been used as a potential marker to distinguish retinoma cells (Rb precursors, Ki67⁻) from Rb and retina (29), stained strongly positive in hRBOs (Fig. 2E and SI Appendix, Fig. S5H). We also confirmed markedly increased expression of other oncogene markers (SYK, DEK, and NSE) in hRBOs, strongly indicating that Rbs and retinoma-like cells were transformed in the cancerous organoids (Fig. 2E and SI Appendix, Fig. S5H). We did not observe any significant change in the expression of p53 (SI Appendix, Fig. S5F), an important tumor repressor. These results support the previous assumption that loss of RB1 function may initiate retinoma and subsequently lead to malignant transformation and progression (29, 31). Taken together, these results demonstrate that organoid Rb recapitulates human Rb in gene expression, DNA methylation, and protein markers.

Tumorigenicity of hRBOs In Vivo.

Since hRBOs retained viability and expanded rapidly after 10 wk in vitro, we attempted to assess the tumorigenicity of cancerous cells in vivo. As shown in Fig. 3A, excessively proliferating Rb-like cells in 120-d-old hRBOs could migrate to the medium in suspension conditions (Fig. 3 A and B and SI Appendix, Fig. S6A). We engrafted these cells into the subretinal or vitreous space of immunodeficient NSG mice. Fundography and optical coherence tomography (OCT) imaging were used to monitor the intraocular grafts (Fig. 3C). At 2 mo postinjection, 9 out of 12 of subretinally injected mice and 9 out of 11 vitreous cavity-injected mice with hRBO-derived Rb-like cells showed visible tumors in the intraocular grafts, but an intraocular tumor was not observed in any eyes injected with hRO-derived retinal cells (Fig. 3D and SI Appendix, Fig. S6 B and D), as further confirmed by hematoxylin and eosin (H&E) staining (Fig. 3E and SI Appendix, Fig. S6C). No significant differences in tumor formation were seen between the subretinal and vitreous cavity grafts. Moreover, immunocytochemistry revealed extensive tumor foci positive for Ki67 and oncogene markers (NSE, SYK, and DEK), indicating the tumor characteristics of Rb (Fig. 3F). TEM also demonstrated a high nucleocytoplasmic ratio, abundant mitochondria with heterogeneous morphology, tightly packed clusters of cells and rosettes, and numerous mitotic figures (Fig. 3G and SI Appendix, Fig. S6E), in line with the findings observed in the in-dish hRBOs. These results demonstrate that hRBOs could retain viability and expand after becoming engrafted into immunodeficient mice.

• Download figure
• Open in new tab
• Download powerpoint

Fig. 3.

Orthotopic engraftment of hRBOs in NSG mice. (A) Representative bright-field image of 120-d-old hRBOs and hROs. Arrows indicate excessive outgrowth of Rb cells, subsequently leading to detachment from the organoids. (Scale bar: 400 µm.) (B) Rb cells derived from hRBOs (hRBO-Rb cells) can be maintained in long-term culture in suspension. (Scale bar: 200 µm.) (C) Schematic view of the procedure for generating tumor engraftments in NSG mice. (D) Ophthalmic examination of mouse eyes after subretinal engraftment. Representative photograph of engrafted eyes at 2 mo after subretinal injection (Left). hRBO-Rb cell-engrafted eye showing leukocoria-like features (black arrow). Representative fundus photography and OCT images of an hRBO-Rb cell-engrafted eye showing an obvious Rb tumor in the retina. Retinal cells derived from hROs (hRO-retinal cells) served as normal controls (Bottom). (E) Representative H&E staining of hRBO-Rb cell-engrafted eyes at 2 mo after subretinal injection. The boxed region is magnified on the right. (Scale bar: 200 µm.) (F) Immunostaining of Rb signature genes Ki67, SYK, DEK, and p16^INK4a in hRBO-Rb cell-engrafted eyes at 3 mo after subretinal injection. (Scale bar: 50 µm.) (G) Representative TEM image of Rb engraftment. The dashed line denotes the typical rosette structure. (Scale bar: 10 µm.)

Cellular Identity of hRBOs by Single-Cell Transcriptomes Analysis.

To define heterogeneity, cell types, and tumor-initiating cells (32, 33), single-cell RNA-Seq analysis was performed in both hRBOs and control hROs at day 90. On average, 4,274 hRBO and 3,080 hRO UMIs and 1,784 hRBO and 1,297 hRO genes were detected per cell. As a result, 9,665 and 12,218 single-cell transcriptomes were obtained from hRBOs and control hROs, respectively. As expected, based on cluster-specific marker genes, there were nine distinct retinal cell types, including retinal progenitor, precursor, cone, rod, bipolar, horizontal/amacrine, ganglion and Müller glial cells in control hROs (Fig. 4 A–C). Notably, we found four additional cell clusters—Rb cells, retinoma-like cells, unfolded protein response (UPR)-related cells, and excessive cone precursors—in hRBOs (Fig. 4 A, B, and D). Among these clusters, the Rb cell cluster highly expressed MKI67, CDK1, TOP2A, KIF14, and CDC25C but not VIM, GLUL, CLU, and FRZB (Fig. 4C and SI Appendix, Fig. S7A). The retinoma-like cell cluster strongly coexpressed Rb-enriched genes, including PCNA and DNA replication licensing factors (MCM2-7) and cone precursor marker genes (Fig. 4C and SI Appendix, Fig. S7A). Different from Rbs, these retinoma-like cells have greatly increased expression of CDKN2A and CDCA7 and reduced expression of MKI67 and KLF14 (SI Appendix, Fig. S7B). Thus, we speculated that retinoma-like cells could be an intermediate cell stage between premalignant cone precursors and tumor cells. Interestingly, the number of normal retinal cells was dramatically reduced while the amount of cone precursors was increased in hRBOs (Fig. 4 A, B, and D). In addition, we identified a cluster of abundantly expressed UPR-related genes (Fig. 4 A–C), indicative of endoplasmic reticulum (ER) stress and essential for tumorigenesis and progression in many cancer types (34⇓–36). Notably, we did not detect up-regulation of proapoptotic genes in this cluster, suggesting that these cells were viable and could represent transformed malignant cells. RNA-Seq and immunostaining analysis also confirmed the presence of these UPR-related cells in hRBOs (SI Appendix, Fig. S7 C and D). In addition, an excessive cluster of cone precursors present only in hRBOs expressed abundant cone precursor markers, including ARR3 and RXRG, implicating the cell of origin of Rbs (Fig. 4C and SI Appendix, Fig. S7A).

• Download figure
• Open in new tab
• Download powerpoint

Fig. 4.

Single-cell RNA-Seq analysis of hRBOs and hROs. (A) t-distributed stochastic neighbor embedding (tSNE) visualization of 21,883 cells in 90-d-old hRBOs and control hROs profiled, with cells color-coded according to cell type annotation. (B) tSNE visualization of 21,883 cells profiled, with cells color-coded according to sample type. (C) Expression of marker genes for the main cell clusters on the tSNE map. (D) Fraction of cells in each cell cluster. (E) Monocle pseudospace trajectory revealing the Rb progression, color-coded according to cell clusters.

To check the effects of cell cycle heterogeneity in the single-cell RNA-Seq data, we reanalyzed the single-cell RNA-Seq data by regressing out cell cycle scores during data scaling. The cell types did not show any obvious changes (SI Appendix, Fig. S7E), suggesting that the cell cycle was not driving cell cluster generation. The violin plots of UMI, mitochondrial, and ribosome counts shown in SI Appendix, Fig. S8 A–C indicate that these counts did not strongly affect clusters. Together, these single-cell transcriptomes reveal distinct cellular diversity and compositions of hRBOs, indicating that the cancerous origin of Rbs may be cone precursors.

Maturing Cone Precursor Is the Cell of Origin of Rbs in hRBOs.

We next sought to confirm the cell of origin of Rbs in the hRBOs using single-cell pseudo-time trajectory analysis and immunocytochemistry. As expected in the retinal developmental trajectory of both the hROs and hRBOs, the population of retinal progenitors was present at the start of the trajectory (Fig. 4E and SI Appendix, Fig. S8 E and F). Notably, across the developmental trajectory specific to hRBOs, ARR3⁺ maturing cone precursors were present at the branch point, followed sequentially by retinoma-like and Rb cells (Fig. 4E and SI Appendix, Fig. S8E). which was consistent with the results of RNA velocity analysis (SI Appendix, Fig. S8D), a recently developed cell fate tracing and cell lineage reconstructing approach (37). These results suggest that the maturing cone precursor cells might develop into premalignant cone precursors, become converted to retinoma-like cells, and eventually terminate at Rb cells (Fig. 4E and SI Appendix, Fig. S8 D, E, and G).

To further elucidate the cell of origin, we compared the retinal cell type-specific gene expression at different organoid stages based on the RNA-Seq data. As shown in Fig. 5 A and B, genes enriched in retinal progenitor, rod, Müller glial, amacrine, horizontal, ganglion, and bipolar cells were significantly down-regulated in hRBOs after day 75. In contrast, the derived hRBOs exhibited dramatic and consistent up-regulation of cone precursor-enriched genes (Fig. 5 A and B). Consistent with these findings, immunostaining revealed an increasing proportion of cone precursor cells (RXRG⁺ or ARR3⁺) and a decreasing number of other retinal neurons in hRBOs with tumorigenesis (Fig. 5 C and D and SI Appendix, Fig. S9A). RNA-Seq data showed no difference between hRBOs and hROs in the expression levels of genes specific for proliferative cells (MKI67), early photoreceptor precursors (THRB and RXRG for cones, NR2E3 for rods, CRX and RORB for cones or rods), or retinal progenitors (OTX2) before day 60 (SI Appendix, Fig. S9B). Remarkably, ARR3 expression was noted at day 60 and showed an abnormal increasing trend thereafter, accompanied by significant up-regulation of MKI67 (SI Appendix, Fig. S9B), suggesting that the abnormal expression of MKI67 may begin with the initial appearance of ARR3⁺ maturing cone precursors in hRBO. Notably, the developmental time window of postmitotic maturing cone precursors (38) exactly matched the tumorigenesis occurring in derived hRBOs from day 60 to day 90 (SI Appendix, Fig. S3B).

• Download figure
• Open in new tab
• Download powerpoint

Fig. 5.

Rbs in hRBOs originating from cone precursors. (A) Gene–gene correlation heatmaps of the indicated cell type signature genes from Rbs, cone precursors/cones, and other retinal cells in hRBOs. The signature genes of Rbs in hRBOs had the greatest correlation with most genes specific for cone precursors. Node size and color indicate the degree of correlation. (B) The fold change of cell type-specific marker gene expression in hRBOs and hROs at different organoid stages on days 30, 45, 60, 75, 90, 105, and 120. Most of the genes specific for cone precursors showed elevated expression in hRBOs compared with hROs from day 60. (C) Immunostaining of maturing cone precursor marker ARR3 and pan-cone precursor marker RXRG in day 105 hRBOs and hROs. (Scale bar: 50 µm.) (D) Quantification of ARR3⁺ or RXRG⁺ cone precursors, NRL⁺ rods, SOX2⁺/SOX9⁺ Müller cells, SOX2⁺/PAX6⁺ retinal progenitors, SOX2⁻/PAX6⁺⁺ ganglion/horizontal/amacrine cells, PROX1⁺ or CALB2⁺ horizontal/amacrine cells, and CHX10⁺⁺ bipolar cells in hRBOs and parental hROs. Data are mean ± SD; n = 5. ***P < 0.001; **P < 0.01, unpaired Student’s t test. Markers measured in cells shown in C and SI Appendix, Fig. S9A. (E) Coimmunostaining of ARR3 and RXRG with the proliferative marker Ki67 in day 105 hRBOs. Magnified views of the boxed regions are shown in the corresponding panel on the right. White arrows indicate Ki67⁺ tumor cells coexpressing ARR3 or RXRG. (Scale bar: 50 µm.) (F) Coimmunostaining of rod, retinal progenitor, ganglion/horizontal/amacrine cell, Müller glial cell, horizontal and amacrine cell, or bipolar cell markers with Ki67 in day 105 hRBOs. Magnified views of the dashed-line boxed regions are shown in the corresponding panels on the right, and the solid-line boxed region is shown in the left panels. Arrows indicate markers specific for different retinal cell types with (white arrows) or without (yellow arrows) Ki67 coexpression. (Scale bar: 50 µm.) (G) Quantification of the percentage of Ki67⁺ cells expressing the indicated markers specific for cone precursors and other retinal cell types in hRBOs as analyzed by immunostaining. Data are mean ± SD; n = 5. ***P < 0.001, unpaired Student’s t test. (H) Schematic diagram summarizing the cell of origin and development of Rbs in hRBOs.

As postmitotic cells, maturing cone precursors (ARR3⁺) are nonproliferative and do not express Ki67 (26). Consistently, the few Ki67⁺ cells detected in the normal hROs after day 75 represented a small number of Müller glial and retinal progenitor cells rather than cone precursor cells (Fig. 5 F and G). In contrast, most ARR3⁺ maturing cone precursors within hRBOs simultaneously expressed Ki67, indicating the mitotic property of the maturing cone precursors in cancerous organoids (Fig. 5 E and G). We did not detect Ki67 in rod, bipolar, ganglion, amacrine, or horizontal cells in hRBOs (Fig. 5 F and G). In addition, most RXRγ⁺/ARR3⁺ cone precursors were positive for Ki67, while RXRγ⁺ cone precursors lacking ARR3 expression in the hRBOs were negative for Ki67 (SI Appendix, Fig. S9C). Notably, CRX⁺ early photoreceptor precursors were postmitotic and negative for Ki67 (17). We found that CRX⁺/ARR3⁻ early photoreceptor precursors were detected in the hRBOs at day 60, and none was positive for Ki67 (SI Appendix, Fig. S9D). In contrast, all ARR3⁺ cone precursors in the hRBOs at day 60 were strongly Ki67⁺, indicative of their cell cycle entry. These Ki67⁺/ARR3⁺ cone precursors were significantly increased after day 60 (SI Appendix, Fig. S9D). As expected, most ARR3⁺ maturing cone precursors were also positive for mitosis-associated phosphohistone H3, PCNA, and proto-oncogenic SYK (SI Appendix, Fig. S9 E–G). Only few UPR-related cells (CHOP⁺) expressed Ki67 in hRBOs (SI Appendix, Fig. S9H). Moreover, in hRBO-derived intraocular tumors, while genes enriched in retinal progenitor, rod, amacrine, horizontal, ganglion, and bipolar cells were barely detectable, ARR3 was detected in the majority of tumor cells and coexpressed with Ki67 (SI Appendix, Fig. S9I). Overall, our results suggest that the ARR3⁺ maturing cone precursor may serve as the originating cell of tumorigenesis in Rb organoids (Fig. 5H). Thus, we propose that RB1 inactivation induces unusual activation of the cell cycle in maturing cone precursors and drives these postmitotic cells to reenter the cell cycle to induce the Rb precursor (retinoma-like cells, represented by up-regulation of PCNA and CDKN2A), after which the cells undergo undue proliferation (indicated by high expression of Ki67) to eventually transform into Rbs (Fig. 5H).

hRBO-Based Drug Testing Identifies Potential Therapeutic Agents.

The high-fidelity Rb model enabled us to perform drug testing. We tested four chemotherapeutic drugs currently in clinical use—vincristine (5 nM), etoposide (0.5 µM), carboplatin (10 µM) and topotecan (10 nM)—as a proof of feasibility. As mentioned above, the activator of the PI3K-Akt signaling pathway, SYK, was aberrantly up-regulated in Rb organoids. Thus, we chose three candidate drugs: SYK inhibitors R406 (5 µM) and BAY 61-3606 (5 µM) (39) and the mTOR (downstream of PI3K-Akt signaling) inhibitor rapamycin (10 µM). We investigated the response to these drugs in hRBOs at the early stage (at day 60, referring to hRBOs at the onset of tumorigenesis) and late stage (at day 90, referring to hRBOs with significant tumorigenesis) (Fig. 6A). After 7 d of treatment and continuous culture to day 120, drug responses were determined by flow cytometry and/or immunocytochemistry of Ki67 and cleaved caspase 3-positive cells (Fig. 6A). As expected, flow cytometry showed that all tested drugs except for etoposide caused a significant reduction in the number of Ki67⁺ cells in hRBOs (Fig. 6 B and C). These results were further confirmed by immunocytochemistry (Fig. 6D and SI Appendix, Fig. S10A). In addition, immunocytochemistry showed a significant increase in cleaved caspase 3 in hRBOs treated with the candidate drugs (Fig. 6D and SI Appendix, Fig. S10A), suggesting that a 7-d administration of these candidate drugs caused remarkable cell apoptosis in cancerous organoids.

• Download figure
• Open in new tab
• Download powerpoint

Fig. 6.

Chemotherapeutic agents for Rbs tested in hRBOs. (A) Schematic workflow of the chemotherapy treatment and assessment using hRBOs. Timelines of the culture and treatment of hRBOs are shown. (B) Representative fluorescence-activated cell sorting (FACS) profile of Ki67⁺ cells in day 120 hRBOs after the indicated early-stage treatment with DMSO; four currently used clinical drugs, including vincristine (VCR; final concentration: 5 nM), etoposide (ETO; 0.5 µM), carboplatin (CAB; 10 µM), and topotecan (TOP; 10 nM); and three candidate drugs, including rapamycin (RAP; 10 µM), R406 (5 µM) and BAY 61-3606 (BAY; 5 µM) (Right). Normal day 120 hROs served as a blank control. The dashed line demarcates the Ki67⁺ gate. (Left) Statistical analysis of the percentage of Ki67⁺ cells in hRBOs after the indicated treatments as analyzed by FACS. Data are mean ± SD; n = 5. *P < 0.05; **P < 0.01; ***P < 0.001, one-way analysis of variance (ANOVA) with Dunnett’s test. (C) Representative FACS profile of Ki67⁺ cells in day 120 hRBOs after the indicated late-stage treatment with DMSO, VCR, ETO, CAB, TOP, RAP, R406, and BAY. Normal day 120 hROs served as a blank control (Right). The dashed line demarcates the Ki67⁺ gate. (Left) Statistical analysis of the percentage of Ki67⁺ cells in hRBOs after the indicated treatments as analyzed by FACS. Data are mean ± SD; n = 5. **P < 0.01; ***P < 0.001, one-way ANOVA with Dunnett’s test. (D) Immunostaining of Ki67 (green) and apoptosis marker cleaved caspase 3 (red) in day 120 hRBOs after treatment with the indicated drugs. (Scale bar: 50 µm.)

To provide evidence for the specificity of tumor cell killing by these candidate drugs, we also performed parallel toxicity testing with normal hROs. The results suggested that carboplatin, rapamycin, R406, and BAY 61-3606 restored normal cell proliferation, but vincristine and topotecan affected cell viability to some extent in normal hROs (SI Appendix, Fig. S10B). Together, these results suggest that our model can be used to test the effect of potential therapeutic drugs on Rb tumors. The inhibitors of SYK, a key therapeutic target of Rb development, could be considered potential therapeutic drugs for Rb.