“…PCA, UMAP, and Louvain clustering were performed with standard parameters in scanpy and velocyto using spliced counts (2,000 selected genes, 20 PCs, 50 nearest neighbors). Global dimensionality reduction of the entire dataset was performed using a profile of pluripotency (6 genes: LIN28A, NANOG, POU5F1, SALL4, SOX2, ZFP42), retinal progenitor (20 genes: BMP4, FGF3, FGF8, LHX2, NCAM1, OPTC, OTX1, OTX2, PAX2, PAX6, RAX, SIX3, SIX6, TBX2, TBX3, TBX5, VAX2, VSX2, ZIC1, ZIC2) and RPE (16 genes: BEST1, ELN, IGFBP5, MITF, PMEL, RGR, RLBP1, RPE65, SERPINF1, SLC7A8, TMEFF2, TRPM1, TRPM3, TTR, TYR, TYRP1) marker genes, and the first two principal components derived from normalized spliced counts were visualized (Bosze et al, 2020;Fuhrmann, 2010;Schmitt et al, 2009). Signature score for apoptosis and stress response were computed using the following 37 genes: CASP8, CASP10, E2F5, CASP2, PTEN, CASP7, CDKN2A, E2F3, E2F4, CASP6, TP53, E2F1, CDKN1B, CDKN1A, E2F2, CASP3, FHIT, CASP9, HSPH1 GPX3, SOCS3, FOS, HSPA1A, SOD2, SOD1, HSF1, SOD3, HSPA1B, HSPA4, CASP6, GPX4, CAT, TNF, IL6, HSPA1L, NR3C1, HSPB1.…”