3.1. Generation of CiPSCs by defined small-molecule compounds
To effectively trace the generation of CiPSCs, Oct4 promoter-driven green fluorescent protein (GFP) reporter MEFs (OG-MEFs) were employed. According to the previously developed method, it takes about 40 days to induce MEFs into CiPSCs by small-molecule compounds[13] (Fig. 1A). As early as day 5 (D5, same abbreviation in below), there were a number of epithelial colonies formed. These colonies rapidly expanded and some of them expressed GFP (Fig. 1B). Subsequently, GFP+ colonies were picked up and seeded on feeder cells to establish CiPSC lines for further characterization and differentiation. These CiPSCs maintained GFP expression and ESC-like morphology throughout passages (Fig. 1C, Figure S1A).
First, CiPSCs were identified to have normal karyotypes (Fig. 1D). Next, the pluripotenct cell markers of CiPSCs were determined at mRNA and protein levels. The pluripotency-associated proteins Oct4, Nanog, Sox2, and SSEA4 could be identified by immunostaining of CiPSCs colonies (Fig. 1E). Meanwhile, the CiPSCs expressed pluripotency-associated transcripts such as Oct4, Nanog, Sox2, Lin28a, and Sall4 and their expression levels were similar to those of mESCs (Fig. 1F). Last, the differentiation potentials of CiPSCs were further examined. In vitro differentiation test showed that CiPSCs could differentiate into different cell types of three germ layers (Figure S1B). In addition, CiPSCs were able to form teratoma after injected into immune-deficient SCID mice (Fig. 1G). These findings demonstrated that CiPSCs were fully pluripotent as other typical iPSC lines.
3.2. Generation of polarized retinal organoids from CiPSCs
We optimized previously reported methods to differentiate 3D ROs (Fig. 2A)[4, 19]. During the formation of embryoid bodies (EBs), the majority of cell aggregates developed neuroepithelium at D5, and a few optic vesicle-like structures appeared at D7 (Fig. 2B). For ROs differentiation, the D7 EBs were transferred into retinal maturation medium, and more optic cup-like structures formed at D10 (Fig. 2B). The optic cups contained retinal progenitors (Rax+, Pax6+, Vsx2+) and retinal pigment epithelium (RPE)progenitors (Otx2+, Rax−) (Fig. 2C). Furthermore, the neural retinae were polarized with the inner side to be basal (Laminin+) and the outer side apical (ZO-1+) (Fig. 2D). When tracing the differentiation of EBs from Oct4-GFP+ CiPSCs, we observed that the GFP+ signal gradually decayed and completely disappeared at D8 (Figure S1C).
To further examine the cell proliferation during differentiation, we performed the immunostaining of organoids for PH3 and Ki67 expression (Fig. 2E). Both PH3+ and Ki67+ cells were gradually decreased from D14 to D21. The proportion of positive cells at D14, 18 and 21 is 9.00 ± 0.73%, 6.56 ± 0.41% and 4.40 ± 0.40% for PH3, and 67.67 ± 3.00%, 53.67 ± 3.01% and 47.13 ± 2.18% for Ki67 (means ± SEM. N > 8), respectively (Fig. 2F-G). Taken together, these results demonstrated that the cells in organoids gradually exited the cell cycles and differentiated into various retinal cell types.
3.3. Morphological and molecular characterization of CiPSC-derived retinal organoids
At the late stage of neuroretinal differentiation, ROs formed big epithelial structures with clear stratification (Fig. 3A). Additionally, outer segment (OS)-like structures were observed at the outer surface of organoids (Fig. 3A), which corresponded to the photoreceptor outer segments in vivo[20].
To investigate the morphology of photoreceptors, we immunostained ROs for photoreceptor markers Recoverin and Rhodopsin. It showed that Recoverin+ and Rhodopsin+ cells were detectable at D21, and both Recoverin+ and Rhodopsin+ cells increased as differentiation proceeded (Fig. 3B). Interestingly, Recoverin and Rhodopsin were located at the apical side of the neural retinal (Fig. 3B), consistent with photoreceptor morphogenesis in vivo. The proportion of positive cells at D18, 21, 25 and 28 was 0.0 ± 0.0%, 14.70 ± 0.73%, 33.56 ± 2.33% and 45.38 ± 1.98% for Recoverin, and 0.0 ± 0.0%, 5.90 ± 0.64%, 8.56 ± 0.73% and 14.22 ± 1.45% (means ± SEM. N > 8) for Rhodopsin, respectively (Fig. 3C), indicating the gradual maturation of the photoreceptors during differentiation. At the transcriptomic level, qRT-PCR analysis of photoreceptor cell-related genes Rcvrn, Rho and Nrl showed a continuous increase of their expression during differentiation (Fig. 3D, Figure S2A-B), which was also in concordance with the developing mouse retina in vivo (Figure S2A-B).
Overall, these findings suggested that CiPSCs could differentiate into retinal organoids with polarized photoreceptor layers, which highly resembled in vivo photoreceptor genesis.
3.4. Development and maturation of retinal organoids
Retinal neurons form stratified structures called inner and outer plexiform layers which are essential for visual signal processing and transferring. On D18, Pou4f1, a marker for RGCs, was detected at the basal side of organoids (Fig. 4A). However, RGCs could not survive throughout the later stages of differentiation (Fig. 4A), probably due to the lack of axon formation[21] and/or the insufficiency of oxygen and nutrients[22] during differentiation, which also happened in other 3D ROs differentiation systems[23].
Prox1 (a marker for horizontal and amacrine cells) and Calretinin (a marker for amacrine and ganglion cells) were detected on D25. Both signals were polarized at the basal side of organoids, and the Calretinin+ cells were located at the inner side of Prox1+ cells, which was the same as the distribution of horizontal cells, amacrine cells, and RGCs in the mouse retina (Fig. 4B). At the same time, the Prox1+;Calretinin+ double-positive cells indicated the generation of amacrine cells (Fig. 4B). As the bipolar cells are late-born cell type during mouse retinal development[20], Prkca, a marker of rod bipolar cells, can only be detected at the late stage of differentiation (Fig. 4C). Moreover, most Prkca+ signals were located at the basal side of the neural retina and the inner side of Rhodopsin+ cells, reflecting the correct distribution of bipolar cells and rod photoreceptors in ROs (Fig. 4C).
Müller glia cells also have essential tasks in maintaining retinal homeostasis and photoreceptor function. On D21, Rlbp1, a marker for Müller glia cells, was detected at the basal side of organoids. The positive signals spanned the whole neural retina and formed the outer limiting membrane (OLM) at the apical side by D28 (Fig. 4D), which was similar to in vivo retinogenesis. Rhodopsin+ rod photoreceptors grew beyond the OLM (Fig. 4D), indicating the formation of the cilia and outer segments.
To examine the maturation of photoreceptors, we conducted immunostaining of organoids for cilium and synapse-related proteins. The connecting cilia marker Arl13b was detected at the apical side of the neural retina, co-labeled with Rhodopsin+ rod photoreceptors (Fig. 4E). Meanwhile, some Rhodopsin+ signals grew beyond the Arl13b+ cilia (Fig. 4E), indicating the formation of the outer segments. Synaptophysin is a marker for synaptic vesicles, normally labeling the outer plexiform layer (OPL) and inner plexiform layer (IPL) in the neural retina. The immunostaining result showed a stratified distribution of synaptophysin staining, representing the OPL and IPL (Fig. 4F) and reflecting the maturation of photoreceptors. In addition, Ctbp2 is a presynaptic protein that typically stains at the rod terminal in the ribbon synapse[24]. Ctbp2 staining was detected at the basal side of the neural retina on D32 and co-localization of Ctbp2 with Recoverin indicated the formation of synaptic connections between photoreceptors and bipolar cell types (Fig. 4G).
The ultrastructures of photoreceptors were further characterized by transmission electron microscopy (TEM). Numerous components of mature photoreceptors such as OLM, connecting cilia (CC), basal cilia (BC), mitochondria-rich inner segment (IS), and outer segment (OS)-like structure could be identified in the TEM images (Fig. 4H).
At the mRNA level, the qRT-PCR analysis of RGCs (Pou4f2, Atoh7), interneurons (Neurod1), and Müller glia (Rlbp1) related genes showed an increase of gene expression during differentiation, which conformed to the developmental mouse retina in vivo (Figure S2C-2E).
These results suggested that CiPSC-derived 3D ROs contained all types of examined retina cell types, including photoreceptors, RGCs, horizontal, amacrine, bipolar and Müller glial cells, and the genesis and distribution of different cell types recapitulated the mouse retinogenesis in vivo. More importantly, the CiPSC-derived photoreceptors should have the ability to receive and transfer visual information by the successful formation of intact structures.
3.5. Dynamic transcriptomic profiling of retinal organoids
To characterize the transcriptome during retinal development, RNA-seq of CiPSC-derived ROs (D0, 14, 21, 28, and 32) and in vivo mouse retinae (embryonic day 18 (E18), postnatal day 1(P1), P7, P14, and P30) was performed (Fig. 5A-5C and Figure S3). The principal component analysis (PCA) of these samples indicated that developmental timing is the most important factor, which, as the first principal component (PC1), can explain 72% and 67.5% of the variance in the two datasets (Fig. 5A, 5B).
The expression of selected cell-type-specific genes in ROs and in vivo mouse retinae was further analyzed (Fig. 5C, Figure S3). Based on IHC data, gene expression patterns during organoid differentiation showed a remarkable resemblance to those of the developing mouse retina in vivo (Fig. 5C, Figure S3). Rod and cone photoreceptor genes showed a progressive increase with differentiation progress, which was consistent with in vivo retinogenesis. Changes in horizontal, amacrine, bipolar, and Müller glial cell-specific gene expression in organoids revealed a striking coincidence with in vivo retinae, suggesting that the CiPSC-derived ROs closely recapitulated in vivo development. Moreover, expression of RGC-specific genes was downregulated with differentiation and was hardly detectable in organoids after D21 (Fig. 5C, Figure S3), consistent with their disappearance in organoids.
3.6. Generation of Crx-tdTomato labeled photoreceptors from CiPSCs
The Crx plays an important role in specifying the fates of photoreceptor cells. Crx expression is first detected in photoreceptor precursors and then continues to express in all photoreceptors as the differentiation proceeds[25]. In order to track the development of photoreceptor cells, we constructed a Crx-reporter CiPSC line by inserting the coding sequence of the tdTomato fluorescent protein prior to the stop codon of Crx gene with CRISPR/Cas9 technology (Fig. 6A, Figure S4A-4C).
From D14 onwards, both the number of tdTomato+ photoreceptors (Fig. 6B, Figure S4D) and the tdTomato intensity (Figure S4E) kept increased, especially at the apical side. The tdTomato+ signals colocalized with Crx+ cells by immunostaining (Fig. 6B), indicating that insertion of tdTomato coding sequence did not alter the expression pattern of Crx gene. Meanwhile, the percentage of tdTomato+ cells constituted around 85% of all cells in the ROs at D32 (Figure S4E), which was consistent with the mature mouse retina in vivo[20]. Additionally, qRT-PCR analysis of Crx also showed a continuous increase of Crx expression during differentiation, which was also in concordance with the developmental mouse retina in vivo (Figure S4F). Furthermore, almost all organoids strongly expressed Crx-tdTomato on D21, demonstrating the consistency in generating tdTomato+ photoreceptors (Fig. 6C). Therefore, the photoreceptors in retinal organoids could be reflected and recorded by tdTomato reporter during differentiation.
3.7. Transplantation of CiPSC-derived photoreceptors into rd mouse model
To determine whether CiPSC-derived photoreceptors (CiPSCs-PR) can restore visual function, single Crx-tdTomato+ photoreceptor cell suspension was transplanted into the subretinal space of P24 rd10 mice, a mouse model of retinal degeneration (Fig. 6A). For comparison, the P24 rd10 mice received Crx-tdTomato+ donor cells in one eye and sham injections (DPBS) in the other (N = 9). Four weeks after transplantation, the pupillary light reflex (PLR) test was conducted to measure photoreceptor function because the light-induced pupil constriction is a behavioral response that depends on photoreceptors’ functional connections with central brainstem targets[18, 26](Fig. 6D). The results showed that the eyes received Crx-tdTomato+ cells had a more sensitive pupil reflex compared with the sham-injected eyes, especially at low light intensities (Fig. 6E).
Three months after CiPSCs-PR transplantation, the eyes were examined by immunofluorescence for evidence of functional integration into the host retina. The tdTomato+ or tdTomato+ Recoverin+ cells were distributed in the ONL (Fig. 6F, 6G). The Prkca+ rod bipolar cells were in close proximity to tdTomato+ cells, indicating the potential formation of synaptic connections between graft photoreceptor cells and host rod bipolar cells (Fig. 6H). Moreover, CiPSCs-PR were found to form synaptic terminals that express the synaptic vesicle protein synaptophysin, and synaptophysin+ synapses were consistently located in the area of the host-graft interface as well as the IPL within the host retina (Fig. 6I), which is essential for transmitting the light signals into the inner retina.
These results demonstrated that some of the transplanted CiPSCs-PR integrated into the host retina, formed functional synaptic connections with downstream retinal neurons, and contributed to visual function.