Molecular Vision 2019; 25:663-678
Received 26 August 2019 | Accepted 08 November 2019 | Published 11 November 2019
Koray D. Kaya,1 Holly Y. Chen,1 Matthew J. Brooks,1 Ryan A. Kelley,1 Hiroko Shimada,1 Kunio Nagashima,2 Natalia de Val,2 Charles T. Drinnan,1 Linn Gieser,1 Kamil Kruczek,1 Slaven Erceg,3 Tiansen Li,1 Dunja Lukovic,4 Yogita K. Adlakha,1,5 Emily Welby,1 Anand Swaroop1
The first three authors contributed equally to this work.
1Neurobiology, Neurodegeneration & Repair Laboratory, National Eye Institute, National Institutes of Health, Bethesda, MD; 2Electron Microscopy Laboratory, National Cancer Institute, Center for Cancer Research, Leidos Biomedical Research, Frederick National Laboratory, Frederick, MD; 3Stem Cell Therapies for Neurodegenerative Diseases Lab and National Stem Cell Bank – Valencia Node, Research Center Principe Felipe, Valencia, Spain; 4Retinal Degeneration Lab and National Stem Cell Bank – Valencia Node, Research Center Principe Felipe, Valencia, Spain; 5Department of Molecular and Cellular Neuroscience, National Brain Research Centre, Manesar, Haryana, India
Correspondence to: Anand Swaroop, National Eye Institute, NIH, Bldg 6/338, MSC0610, 6 Center Drive, Bethesda, MD 20892-0610; Phone: (301) 435 5754; FAX: (301) 480 9917 www.nei.nih.gov/intramural/nnrl, email@example.com.
Dr. Shimada is now at Department of Physiology, Keio University School of Medicine, Tokyo, 160-8582, Japan
Purpose: Retinal organoids generated from human pluripotent stem cells exhibit considerable variability during differentiation. Our goals are to assess developmental maturity of the neural retina in vitro and design improved protocols based on objective criteria.
Methods: We performed transcriptome analyses of developing retinal organoids from human embryonic and induced pluripotent stem cell lines and utilized multiple bioinformatic tools for comparative analysis. Immunohistochemistry, immunoblotting and electron microscopy were employed for validation.
Results: We show that the developmental variability in organoids was reflected in gene expression profiles and could be evaluated by molecular staging with the human fetal and adult retinal transcriptome data. We also demonstrate that the addition of 9-cis retinal, instead of the widely used all-trans retinoic acid, accelerated rod photoreceptor differentiation in organoid cultures, with higher rhodopsin expression and more mature mitochondrial morphology evident by day 120.
Conclusion: Our studies provide an objective transcriptome-based modality for determining the differentiation state of retinal organoids and for comparisons across different stem cell lines and platforms, which should facilitate disease modeling and evaluation of therapies in vitro.
Human development requires stringent and coordinated control of gene expression, signaling pathways, and cellular interactions that result in the generation of distinct cell types and tissues with complex morphological and functional phenotypes [1,2]. However, most of our current knowledge of the fundamental molecular events underlying cell-type specification and tissue differentiation has been derived from model organisms. Despite studies using human preimplantation embryos [3,4] and fetal tissue [5-7], the complexities of human organogenesis are poorly understood [8,9]. Pioneering advances in the generation of human embryonic stem cells (hESCs)  and induced pluripotent stem cells (iPSCs) , together with the development of three-dimensional (3D) organoid cultures [12-14], have revolutionized the studies of human development, facilitated individualized disease modeling, and rejuvenated the field of regenerative medicine [15-18].
Retinogenesis begins with specification of the forebrain neuroectoderm, and patterning of the early eye field is governed by finely tuned regulatory networks of signaling pathways and transcription factors . Distinct morphological changes in the rostral neuroectoderm involve lateral expansion of bilateral eye fields to produce the optic vesicles, which invaginate and become the optic cups . Landmark studies using fetal and neonatal tissues have provided unique insights, distinct from model organisms, into human retinal development [21-24]. By providing appropriate systemic and exogenous cues, human pluripotent stem cells can be directed to self-organize into 3D optic vesicle or optic cup structures [13,14]. Retinal organoids mimic early eye field development, with VSX2+ (also called CHX10) multipotent retinal progenitor cells differentiating into a polarized and laminated architecture harboring all types of retinal neurons and the Müller glia . Rod and cone photoreceptors in organoid culture express opsin and other phototransduction genes [26-28], and develop rudimentary outer segment-like structures at late stages of differentiation [29,30]. However, current methods for characterizing retinal organoids have largely relied upon expression of select cell-type specific markers and histology, which provide limited information about the precise differentiation status. Although live imaging modalities have been employed recently for characterization and developmental staging [31,32], we still lack molecular insights into retinal organoid differentiation and maturation on a global scale, and how different experimental conditions (e.g., cell lines and protocols) could impact organoid cultures.
In this study, we performed transcriptome profiling of developing retinal organoids generated from hESCs and hiPSCs, using modifications of a widely used protocol . Comparative transcriptome analyses with gene profiles of human fetal and adult retina revealed the molecular stages of retinal organoids and demonstrated their differentiation status and cellular composition more accurately. We also identified a specific role of 9-cis retinal (9CRAL) in expediting rod photoreceptor differentiation compared to the currently used all-trans retinoic acid (ATRA). Thus, these studies established a transcriptome-based molecular staging system using human fetal and adult data, enabling direct comparison of organoids under different experimental conditions for disease modeling and evaluation of therapies.
CRX-GFP H9 is a subclone of the H9 human ESC line, carrying a green fluorescent protein (GFP) gene under the control of the cone‐rod homeobox (CRX; Gene ID: 1406, associated OMIM numbers: 120970, 613829, 268000) promoter as previously reported . The human iPSC lines PEN8E and NEI377 were reprogrammed from skin biopsies using integration-free Sendai virus carrying the four Yamanaka factors, as described , and their genome integrity and pluripotency have been evaluated . The ESP1 and ESP2 lines were reprogrammed by Oct4, Klf4, Sox2, c-Myc, and Lin-28 mRNAs  and the Sendai virus (Ctrl1 FiPS4F1, Spanish National Stem Cell Bank, Valencia, Spain), respectively. The H9, PEN8E, and NEI377 cell lines were maintained in Essential 8 medium (E8; ThermoFisher Scientific, Waltham, MA) under hypoxia (5% O2), and the ESP1 and ESP2 lines in mTeSR1® (Stem Cell Technologies, Vancouver, Canada) under normoxia (20% O2). All lines were sustained on BD MatrigelTM human embryonic stem cell-qualified Basement Membrane Matrix (Corning, New York, NY)-coated plates. The PSCs were passaged every 3–4 days at 70–80% confluency using the EDTA-based protocol .
To start differentiation, all lines were detached and dissociated into small clumps with the EDTA dissociation protocol , and cultured in polyHEMA-coated or Ultra Low Attachment Culture Dishes (Corning) in E8 with 10 μM Y-27632 (Tocris) to form embryoid bodies (EBs). Neural-induction medium (NIM) consisting of Dulbecco's Modified Eagle's medium (DMEM/Ham's F-12 Nutrient Mix (F-12) (1:1), 1% N2 supplement (ThermoFisher Scientific), 1X Minimum essential medium (MEM) non-essential amino acids (NEAA; MilliporeSigma), and 2 µg/ml heparin (MilliporeSigma) was added at differentiation day (D)1 and D2 to reach a final ratio of 1:3 and 1:1 NIM:E8, respectively. Starting at D3, EBs were cultured in 100% NIM for an additional 4 days. EBs were collected at D7 and plated onto BD MatrigelTM Growth Factor Reduced (GFR) Basement Membrane Matrix (Corning)-coated plates. On D16, the medium was changed to photoreceptor induction medium (PIM) consisting of DMEM/F-12 (3:1) with 2% B27 without vitamin A, 1% antibiotic-antimycotic solution (ThermoFisher Scientific), and 1X NEAA. Medium was changed daily.
Upon the appearance of optic vesicles (OVs; typically, D21–28) with neuroepithelium morphology, the regions were excised with tungsten needles under magnification and transferred to suspension culture in polyHEMA-coated or Ultra Low Attachment 10 cm2 Culture Dishes (Corning) in PIM. After D42, PIM was supplemented with 10% fetal bovine serum (FBS; ThermoFisher Scientific), 100 mM taurine (MilliporeSigma), and 2 mM GlutaMAX (ThermoFisher Scientific). Starting at D63, the PIM was supplemented with 1 µM retinoid (MilliporeSigma) three times a week and switched to 0.5 µM at D92. The 1% N2 supplement was included starting at D92. The medium for the H9 and PEN8E-derived organoids was supplemented with 20 ng/ml insulin-like growth factor 1 (IGF1; ThermoFisher Scientific) and 55 nM beta-mercaptoethanol (2-ME; ThermoFisher Scientific) from dissection until the end of differentiation. 9-cis retinal (9CRAL; MilliporeSigma) was added at the described concentrations and corresponding time periods during the media change. Organoids for the data set PEN8E_2 were cultured with 10 ng/ml IGF1 from D42 to D200, and ATRA (MilliporeSigma) was supplemented to the PIM at 1 µM 5 times per week and at 0.5 µM twice a week for D63–92 and D92–180, respectively. The ESP-derived organoid cultures were maintained in the same manner as the PEN8E_2 organoids, except no IGF1 was used. For all experiments directly comparing 9CRAL and ATRA, media changes were performed in the dark under dim red light to reduce isomerization of the retinoids. Retinal organoids were collected in the dark and processed in the light as described for all other experiments.
hPSC-derived retinal organoids were collected throughout the differentiation process, fixed, and immunostained as previously described . Briefly, the organoids were fixed in 2% or 4% paraformaldehyde in 1X PBS (D-PBS, 1X, 138 mM NaCl, 2.6 mM KCl, 8 mM Na2HPO4-7H2O, 1.5 mM KH2PO4, pH 7.4) for 1 h at room temperature. Fixation was followed by cryoprotection in a sucrose gradient, 10–30%. Retinal tissues were then embedded in Shandon M1 embedding matrix and frozen on dry ice. Ten-micron sections were obtained using a Leica cryostat (Leica Biosystems, Buffalo Grove, IL) at −14 °C. The sections were placed on Superfrost Plus slides (Fisher Scientific, Hampton, NH) and stored at −20 °C until use. Specific antibodies and concentrations are summarized in Appendix 1.
The protocol for immunoblot analysis was as previously described . Individual organoids were lysed in 50 μl of 1% Triton X-100 (MilliporeSigma) in 1X PBS (ThermoFisher Scientific) supplemented with 1X protease inhibitors (Roche, Basel, Switzerland). Gentle trituration was then used to dissociate the organoid until a homogenous solution was obtained. The samples were then incubated on a nutator at 4 °C for 1 h and centrifuged at 4 °C at 4,000 ×g for 5 min. The protein concentration was measured via the Pierce bicinchoninic acid (BCA) protein assay (ThermoFisher). Approximately 15 μg supernatant protein was diluted 4:1 in reducing 4X Laemmeli buffer for 1 h at room temperature and separated at 100 V for 1 h on 10% sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) gels, which were then transferred to polyvinylidene fluoride (PVDF) membranes at 100 V for 1 h. After blocking in 5% milk in 1X Tris-Buffered Saline, 0.1% Tween® 20 Detergent (TBST; Santa Cruz Biotech, Dallas, TX) for 1 h at room temperature, the blots were incubated in an antibody cocktail (1:1,000 1D4, 3A6, and 4D2, a generous gift from Dr. Robert Molday, Vancouver, Canada) overnight in 1% milk in 1X TBST. The following morning, the membranes were washed in 1X TBST and shaken at room temperature four times for 10 min each. The membranes were then incubated in appropriate secondary species immunoglobulin G (IgG) conjugated to horseradish peroxidase (1:10,000) in 1X TBST shaking at room temperature for 1.5 h. The membranes were then washed again four times in 1X TBST for 10 min each. Before imaging, the membranes were exposed to SuperSignal® West Pico enhanced chemiluminescence (ECL) solution (ThermoFisher) for 5 min, and chemiluminescence was captured using a Bio-Rad ChemiDoc™ touch (Bio-Rad, Hercules, CA). Membranes were next dried overnight at room temperature to remove antibody binding to the protein. Dry membranes were reactivated in 100% methanol and incubated with gamma-tubulin (1:1,000, ab11317, Abcam, Cambridge, United Kingdom) overnight as described above. Images were analyzed in Image Lab (Bio-Rad) and exported to Adobe Photoshop (San Jose, CA) for figure generation.
Retinal organoids were processed for transmission electron microscopy (TEM) analysis as previously described . Briefly, the organoids were initially fixed in 4% formaldehyde and 2% glutaraldehyde in 0.1 M cacodylate buffer (pH 7.4; Tousimis, Rockville, MD) for 2 h, washed in cacodylate buffer three times, and then fixed in osmium tetroxide (1% v/v in 0.1 M cacodylate buffer; Electron Microscope Science, Hatfield, PA) for 1 h in a room temperature. The organoids were washed in the same buffer three times, followed by acetate buffer (0.1 M, pH 4.2), and en-bloc staining in uranyl acetate (0.5% w/v; Electron Microscope Science) in acetate buffer for 1 h. The samples were dehydrated in ethanol solution (35%, 50%, 75%, 95%, and 100%), followed by propylene oxide. Subsequently, the samples were infiltrated in a mixture of propylene oxide and epoxy resin (1:1) overnight, embedded in a pure epoxy resin in a flat mold, and cured in a 55 °C oven for 48 h. Thin sections (70 to 80 nm) were made with an ultramicrotome (UC 7) and diamond knife (Diatome, Hatfield, PA), attached on 200-mesh copper grid, and counter-stained in aqueous solution of uranyl acetate (0.5% w/v) followed by lead citrate solutions. The thin sections were stabilized with carbon evaporation under a vacuum evaporator before the EM examination. The digital images were taken in the electron microscope (H7650)-equipped digital camera (AMT, Rosemont, IL).
High-quality total RNA (100 ng, RIN >7) from at least two independent replicates for hPSCs at various differentiation time points were subjected to mRNA directional library construction as described previously [26,38]. Paired-end sequencing was performed to a length of 125 bases using HiSeq2500 (Illumina, San Diego, CA). Genome reference sequence GRCm38.p7 and Ensembl v82 annotation was used for alignment and quantitation. Quality control, sequence alignment, transcript, and gene-level quantification of primary RNA sequence (RNA-seq) data were accomplished using an established bioinformatics pipeline  Appendix 2. The transcriptome of adult human retina samples was obtained from NEI Commons (Brooks, MJ and Swaroop, A).
Data generated for the five different cell lines, transcriptome comparison and 9CRAL/ATRA analysis were analyzed separately. For each data set, genes were kept for further analysis only if there were ≥5 count per million (CPM) in all replicates of at least one group of the data sets. The data were subjected to trimmed mean of log expression ratios values (TMM) normalization with edgeR v3.20.9 [40,41], and principle component analysis (PCA) and Pearson correlation were performed with normalized CPM (log2) values. Differential expression analysis was performed using limma v3.34.9 , and genes that had ≥2-fold change at least between two time points in either transcriptome of comparison and a false discovery rate (FDR) ≤0.01 were considered statistically significantly differentially expressed.
Gene expression clustering was performed on gene-wise Z-scores using AP (apcluster v1.4.7 package) [43,44] with the k-means option, where k=12 by choosing the corSimMat function as the similarity measure. Determination of the cluster number k was accomplished by empirical observation of the heatmap produced from the PCA projections of the genes-wise z-score.
GO enrichment analysis was performed using clusterProfiler v3.6.0 . To reduce redundancy associated with GO analysis, we performed Wang semantic similarity comparisons of the enriched GO terms using GOSemSim v2.4.1  and semantic similarities clustered with AP. The term with the highest level in each cluster was taken as the enriched GO term. Heatmaps of GO term gene lists were generated using the log2 CPM values. The GO enrichment ratio was calculated with reference to the number of pathway genes found overlapping with the analysis.
Open-ended dynamic time warping (OE-DTW) analysis was used to compare the maturation state of the different retina organoid time-course stage expression data to the human fetal retina (GSE104827) and the human adult retina, as previously described . Retina organoid time-course stages were defined as Stage 1: 25–36 days, Stage 2: 50–75 days, Stage 3: 80–90 days, Stage 4: 105–125 days, Stage 5: 145–172 days, and Stage 6: 186–205. Mean gene expression values at each stage were used for OE-DTW analysis. The gene set used for the OE-DTW analysis consisted of 118 well-known retina-centric transcription factors and cell-type markers (see Appendix 3 and Appendix 4). Co-inertia analysis was performed using MADE4 v1.52.0 package in R.
Two human pluripotent stem cell lines, H9 (hESC) and PEN8E (hiPSC), were differentiated into retinal organoids using a widely used robust protocol , with minor modifications based on our mouse organoid differentiation method that included 9CRAL . Both lines could generate phase-bright neural retina with cone and rod photoreceptors of comparable morphology (Figure 1A–C). S opsin (OPN1SW) was polarized to the apical side of the organoids (outer surface of organoid, exposed to media) as early as D90, with a further increase in expression from D120 to D200. L/M cone photoreceptors (labeled with L/M opsin, OPN1L/MW) were barely evident at D120 but increased at D200 (Figure 1C, left). Rhodopsin (RHO) immunostaining was observed at the apical side as early as D120 in the H9- and PEN8E-derived organoids, and elongated as the photoreceptors matured (Figure 1C, right). Retinal ganglion cells (RGCs) were evident at early stages, but could not be maintained through the end of differentiation, while all other neural retina cell types showed similar morphology and development in the H9- and PEN8E-organoids (Appendix 5).
We then performed RNA-seq analysis of the developing retinal organoids to decipher the major differentiation stages in vitro. The H9- and PEN8E-derived retinal organoids showed a substantial overlap of a total of 3,851 differentially expressed (DE) genes that exhibited statistically significant changes in expression during differentiation (Figure 2A). The 1,041 DE genes unique to H9 and the 1,749 genes unique to PEN8E showed the same general trend in both data sets (Appendix 6). We noted that many of these genes were statistically significantly differentially expressed in both data sets under a less stringent cutoff of 1.5-fold and 5% FDR (data not shown). AP clustering analysis of the common DE genes between the two datasets yielded eight clusters (Cs) with either monotonically increasing or decreasing expression (Figure 2B). GO Biologic Process enrichment analysis uncovered key biologic pathways and genes that could be associated with specific stages of retinal development, and identified statistically significantly enriched terms for all clusters except C1 and C7 (Figure 2C). The C2 through C4 genes displayed a gradual decrease throughout differentiation, and were associated with negative regulation of neuronal differentiation, mitotic cell division, and neurodevelopmental processes and signaling (Figure 2D, Appendix 7). C4 also included axon guidance genes, which likely reflected the loss of ganglion cells in the organoid cultures. In contrast, the C5, C6, and C8 genes exhibited progressively increasing expression during organoid development, and included genes associated with neuronal or retinal differentiation and functions, including visual perception, synaptogenesis, and phototransduction.
To evaluate the impact of iPSC lines or modifications in the differentiation protocol, we compared the expression of the 3,851 DE genes with organoids produced from different hiPSC lines (ESP1, ESP2, and NEI377) and protocols (PEN8E_2). PCA revealed a similar temporal progression in all organoid transcriptomes with increments on the differentiation days, as evident from PC1 (Appendix 4). As the harvesting time points varied among individuals and laboratories, we organized the organoid samples into six groups, for ease of comparison, based on the collection day (chronological age) and state of differentiation (Figure 3A) and performed PCA (Figure 3B) using a selection of established, retina-centric genes (Appendix 3). Although all organoids clustered together at an early stage of differentiation (group 1), the H9- and PEN8E-derived organoids (which included 9CRAL supplementation) showed accelerated differentiation from group 3 and onward, compared to the ESP, NEI377, and PEN8E_2 organoids (which used ATRA as described by ). In concordance, the ESP, NEI377, and PEN8E_2 organoids displayed lower expression of genes associated with phototransduction and outer segments compared to the H9- and PEN8E-derived organoids (Appendix 4). The NEI377 and PEN8E_2 organoids were generated by an identical protocol, and they displayed a similar developmental pattern, except group 6 of PEN8E_2, which is likely due to iPSC line differences.
To determine the maturity of organoid cultures using objective parameters, we performed OE-DTW analysis to align the time-series expression data between the organoids and the human fetal and adult retina gene profiles (Figure 3C) using a set of retina-centric genes (Appendix 3). Across all stem cell lines, group 1 organoids (D25–37) matched with the earliest human fetal time points (D52–67). Subsequently, the distinct PSC organoids showed differences in their highest correlation to corresponding human retinal stages. While organoids derived from ESP1, ESP2, NEI377, and PEN8E_2 lines revealed higher correlation with fetal retina samples; groups 4–6 H9 and PEN8E organoids were more concordant with the late fetal and adult retina (Figure 3C). For example, the group 4 H9 and PEN8E-derived organoids (D145–172) were highly correlated with adult samples (17, 24, and 28 years of age), whereas the ESP1, ESP2, NEI377, and PEN8E_2 organoids corresponded with late fetal samples (D105–136) at a similar stage (Figure 3C).
To confirm the result of the DTW analysis, we compared the expression of mature retinal genes associated with phototransduction and synaptic function between organoids and in vivo retinas. The late stage (group 5 and 6) H9 and PEN8E organoids showed similar expression of phototransduction genes as that of the adult retina (Figure 3D). However, most genes from “Transmission across chemical synapses” pathways were downregulated in the organoid cultures, suggesting the retinal organoids were not yet fully mature (Appendix 4).
Much like its physiologically relevant isomer 11-cis retinal, 9CRAL can also bind opsin to form functioning rhodopsin, and retinoic acid (RA) can be generated as an oxidation product. Therefore, we hypothesized that replacing RA with 9CRAL in the modified protocol would result in improved rod differentiation. To evaluate this hypothesis, we performed a direct comparison of the D90 and D120 H9 organoids cultured with ATRA or 9CRAL from D64 onward. We then applied coinertia analysis to the transcriptomes of the ATRA- and 9CRAL-supplemented organoids and the developing human fetal retina (D53 to D136). The D90 and D120 of the 9CRAL organoid transcriptomes corresponded with the D80–94 and D125–136 human fetal retina gene profiles, respectively. However, the D90 and D120 ATRA transcriptomes matched with earlier fetal retina time points (D57–67 and D105–115, respectively; Figure 4A). Thus, treatment with 9CRAL expedited retinal development in organoid cultures and showed closer temporal transcriptome dynamics to the human fetal retina. The DE genes between D90 and D120 in the ATRA (1284 genes) and 9CRAL (852 genes) data presented an overlap of 594 genes (Figure 4B), which could, in turn, be grouped into five clusters with hierarchical clustering analysis (Figure 4C). Redundancy-reduced GO analysis of the common DE genes revealed a monotonic decrease in the C1 and C2 genes involved in the cell cycle, suggesting neural progenitor cells exit the cell cycle and become mature in retinal organoids (Figure 4D). The C3 through C5 genes were involved in energy metabolism and visual perception, which further characterized their progressive increase from D90 to D120. Although the ATRA and 9CRAL organoids showed comparable developmental patterns, changes in the DE genes during development varied between the two groups (Figure 4C). Day-matched comparison of DE genes revealed lower expression of genes involved in regulation of the RA signaling pathway, RA metabolic processing, and energy metabolism, whereas higher expression of visual perception and function genes was evident in the D90 and D120 9CRAL-treated organoids (Figure 4E). The heatmap of early eye field transcription factors and photoreceptor genes consistently showed expedited photoreceptor differentiation in the 9CRAL organoids, as revealed by the higher expression of these genes, compared to the ATRA organoids (Figure 4F).
To validate the findings that 9CRAL accelerated photoreceptor differentiation, we performed immunohistochemistry on the D90 and D120 organoids. The D90 9CRAL organoids showed lower expression of retinal progenitor cell marker VSX2 (also called CHX10) and higher expression of a pan-photoreceptor marker (recoverin, RCVRN) compared to the ATRA organoids (Figure 5A, upper panel). At D120, the 9CRAL-treated organoids showed distinctively more RHO+ cells, which were barely detected in the ATRA organoids (Figure 5A, bottom panel), consistent with the higher rhodopsin expression in 9CRAL (14.0±9.00 CPM) compared to the ATRA (3.8±1.7 CPM) RNA-seq data. No marked difference was observed in the expression or morphology of the L/M-opsin and S-opsin+ cells (see Figure 5A, bottom panel), as well as other neural retina cell types (Appendix 8). In concordance, immunoblot analysis demonstrated increased rhodopsin expression in the 9CRAL organoids compared to the ATRA organoids at D120 (Figure 5B). Ultrastructural analysis of D130 retinal organoids using transmission electron microscopy (TEM) showed an advanced stage of maturation of the 9CRAL photoreceptors compared to the ATRA photoreceptors (Figure 5C). The ellipsoid (apical) side of the photoreceptor inner segments in the 9CRAL organoids contained a higher number of mitochondria with typical morphology, which could explain the divergent expression of genes involved in energy metabolism between the 9CRAL and ATRA organoids. No differences were evident in the morphology of the photoreceptor connecting cilia between the two groups. Taken together, these data demonstrate that 9CRAL supplementation expedited differentiation and maturation of rod photoreceptors in retinal organoids.
Three-dimensional organoids derived from human pluripotent stem cells have been extensively applied for investigating organogenesis, modeling diseases, and developing therapies. However, substantial variations within organoids pose challenges for comparison among different cultures and studies. We generated transcriptomes of multiple distinct retinal organoids and compared them to human fetal and adult retina gene profiles for molecular staging of differentiation state of the cultures. The analysis revealed the advantage of using 9-cis retinal, instead of the widely used all-trans retinoic acid, in facilitating rod photoreceptor differentiation. Thus, a transcriptome-based comparison can provide an objective method for uncovering the maturity of organoid cultures across different lines and in various study platforms.
The generation of in vitro 3D organoids from pluripotent stem cells has permitted rapid advances in our understanding of human organogenesis, disease mechanisms, and therapeutic interventions [16,18]. Relatively easy access and promising applications of cell replacement therapy to alleviate blinding diseases have prompted an explosion of studies on human retinal organoids, which exhibit appropriate stratified architecture and differentiation of all relevant cell types [26,29-31,48-51]. However, a lack of appropriate photoreceptor outer segments and synaptic connectivity, and the loss of RGCs after prolonged cultures, has hampered progress. The variability in temporal differentiation, depending on protocols and iPSC lines, also demands objective criteria for staging of human retinal organoids. Live imaging and reporter quantification assays have been used to characterize organoid development and staging [31,52]. More recently, retinal organoids from 16 hPSC lines were examined using multiple structural criteria that were then employed for developing a staging system . In comparison, the present molecular staging provides a more objective measure of organoid maturity based on molecular staging with human retinal development. The present molecular staging profiles are concordant with a recent report describing similar temporal gene expression between in vivo human retina and in vitro cone-rich retinal organoid transcriptome data .
Human iPSCs harbor epigenetic memory of their somatic tissue of origin, which appears to favor their subsequent differentiation toward the lineage related to the donor cells and restrict other cell fates . The presence of epigenetic signatures of donor tissues depends on the reprogramming methods and passages of iPSCs ; therefore, the differentiation capacity and therapeutic potential of hESCs and hiPSCs are not easily comparable. Interneurons in mouse retinal organoids differentiated from iPSCs cannot be well differentiated due to epigenetic markers of fibroblasts , but this phenotype could be alleviated in cultures with additional nutrients [36,39], suggesting that culture conditions and cell lines impact the epigenome in iPSCs. In this report, we differentiated genetically unmatched hESC and hiPSC lines reprogrammed from fibroblasts into retinal organoids, and the data showed a similar differentiation capacity of hESCs and hiPSCs in retinal lineages. The present results are in agreement with a recent study showing equivalent gene expression and neuronal differentiation potential between hESCs and iPSCs .
Human organoid cultures manifest substantial variability that may result from multiple factors, including intrinsic genetic variations and the epigenome state of the iPSC lines, reprogramming method, and differentiation protocol. A recent report used a TaqMan array-based analysis of key marker genes to demonstrate differences in retinal organoids . The present comprehensive transcriptome-based molecular staging method utilizes global gene expression profiles and can robustly capture organoid variability and developmental status by comparing them to in vivo retinal transcriptome data. The gene profiles of ESP1, ESP2, NEI377, and PEN8E_2 showed broadly similar gene profiles and developmental trajectories, suggesting a minimal impact of cell-line intrinsic and protocol variations. However, the H9 and PEN8E lines clearly exhibited a more mature developmental trajectory and transcriptomes, resulting from the use of 9CRAL in the organoid cultures.
The studies reported here demonstrated accelerated rod photoreceptor differentiation and detection of rhodopsin protein with immunoblot analysis in 9CRAL-supplemented organoids as early as D120. We suggest that 9CRAL is driving more cells toward a photoreceptor cell fate. In the experimental protocol, all media changes were performed in the dark under dim red light to avoid isomerization of 9CRAL to all-trans retinal; therefore, at least some of the 9CRAL is likely oxidized into 9CRA inside the cells. 9-cis retinoic acid (9CRA) is a potent agonist for retinoid X receptors (RXRs) and retinoic acid receptors (RARs), while ATRA binds only to RARs [58-60]. Retinoic acid has been shown to promote photoreceptor development [61-63] and induce the expression of rod differentiation factor NRL . We should note that retinoid-related orphan receptor beta (RORβ) regulates rod development by activating NRL . In general, the interplay of retinoic acid receptors with other nuclear receptors has a substantial impact on transcriptional regulation of genes involved in photoreceptor development . Thus, expedited rod differentiation by 9CRAL could result from more potent activation of retinoic acid receptors that may induce rod genes through a NRL-regulated gene network . Nevertheless, further investigations are needed to elucidate underlying mechanisms of expediated rod photoreceptor differentiation via 9CRAL. We note that transcriptome profiling of retinal organoids derived from mouse pluripotent stem cells was recently used to identify missing regulatory signals and produce a more mature neural retina in vitro using docosahexaenoic acid (DHA) and fibroblast growth factor 1 (FGF1) .
Appendix 1. Antibody information for immunostaining.
Appendix 2. Sample sequencing depth and transcriptome alignment statistics.
Appendix 3. Retina-centric genes used for DTW analysis.
Appendix 4. Comparative transcriptome analysis of retinal organoids and human retinal samples (relates to Figure 3).
Appendix 5. Immunohistochemistry images showing differentiation of retinal cell types in H9 and PEN8E retinal organoids (relates to Figure 1).
Appendix 6. Genes uniquely expressed in H9 or PEN8E organoids (relates to Figure 2).
Appendix 7. Heatmaps of CPM (log2) gene expression values from selected enriched GO Biological Process terms within clusters C2, C3, and C5 from Figure 2.
Appendix 8. Immunohistochemistry images of H9 retinal organoids treated with all-trans retinoic acid (ARTA, left) or 9-cis retinal (9CRAL) (relates to Figure 5).
We are grateful to Drs. Samuel G. Jacobson and Brian Brooks for skin biopsy samples that were used for generating fibroblasts and subsequently iPSC lines at the Stem Cell Core facility of National Heart, Lung and Blood Institute. We thank Comparative Cytogenetics Core Facility of National Cancer Institute for karyotyping assay. We acknowledge Anupam Mondal, Ben Fadl, Lina Zelinger and Samantha Papal for insightful discussions and constructive comments, and Jacob Nellissery and John Wilson for technical maintenance. This research was supported by Intramural Research Program of the NEI (ZIAEY000450 and ZIAEY000456). Y.K.A. was partially supported by Department of Science and Technology (DST) INSPIRE Faculty grant (IFA13-LSBM-90). Stem cell work in Spain was supported by Spanish Ministry of Science, Innovation and Universities, Instituto de Salud Carlos III (ISCIII)-European Regional Developmental Fund (FEDER) PI16/00409, CP18/00033 and ISCIII-FEDER Platform for Proteomics, Genotyping and Cell Lines, PRB3, PT17/0019/0024. The EM work was funded by FNLCR Contract HHSN261200800001E. The bioinformatic analyses used the high-performance computational capabilities of the Biowulf Linux cluster at NIH. Data availability statement: All raw and processed data are available through Gene Expression Omnibus with accession GSE129104 and at https://neicommons.nei.nih.gov.