1. Developmental Biology
  2. Stem Cells and Regenerative Medicine

Fish primary embryonic pluripotent cells assemble into retinal tissue mirroring in vivo early eye development

  1. Lucie Zilova
  2. Venera Weinhardt
  3. Tinatini Tavhelidse
  4. Christina Schlagheck
  5. Thomas Thumberger
  6. Joachim Wittbrodt  Is a corresponding author
  1. Centre for Organismal Studies Heidelberg, Heidelberg University, Germany
  2. Heidelberg International Biosciences Graduate School HBIGS and HeiKa Graduate School on “Functional Materials”, Germany
https://doi.org/10.7554/eLife.66998
Share this article
Cite this article
  1. Lucie Zilova
  2. Venera Weinhardt
  3. Tinatini Tavhelidse
  4. Christina Schlagheck
  5. Thomas Thumberger
  6. Joachim Wittbrodt
(2021)
Fish primary embryonic pluripotent cells assemble into retinal tissue mirroring in vivo early eye development
eLife 10:e66998.
https://doi.org/10.7554/eLife.66998
  2. Download BibTeX
  3. Download .RIS

Abstract

Organoids derived from pluripotent stem cells promise the solution to current challenges in basic and biomedical research. Mammalian organoids are however limited by long developmental time, variable success, and lack of direct comparison to an in vivo reference. To overcome these limitations and address species-specific cellular organization, we derived organoids from rapidly developing teleosts. We demonstrate how primary embryonic pluripotent cells from medaka and zebrafish efficiently assemble into anterior neural structures, particularly retina. Within 4 days, blastula-stage cell aggregates reproducibly execute key steps of eye development: retinal specification, morphogenesis, and differentiation. The number of aggregated cells and genetic factors crucially impacted upon the concomitant morphological changes that were intriguingly reflecting the in vivo situation. High efficiency and rapid development of fish-derived organoids in combination with advanced genome editing techniques immediately allow addressing aspects of development and disease, and systematic probing of impact of the physical environment on morphogenesis and differentiation.

Introduction

Organ development is a complex process of orchestrated events of tissue specification, morphogenesis, and differentiation in specialized cell types. In vertebrates, retinal development is initiated during early neurulation, when the eye field is specified by the coordinated expression of eye field transcription factors Rx, Lhx2, Pax6, Six3, and Otx2 within the region of the anterior neural plate (Li et al., 1997; Loosli et al., 1999; Zuber et al., 2003). The first morphogenetic sign of retinal development is the formation of the optic vesicle (OV) that evaginates from the wall of the developing diencephalon (Figure 1a, 0–1 day post-fertilization [dpf]). Although the molecular mechanism of OV formation is not completely understood, studies performed in different vertebrate models indicate that the retinal-specific homeodomain transcription factor Rx is involved in this process (Loosli et al., 2003; Mathers et al., 1997; Medina-Martinez et al., 2009; Rembold et al., 2006; Stigloher et al., 2006). As one of the earliest genes indicative for the retinal lineage, Rx is expressed in the anterior neural plate and later in the neuroepithelium of the OV. In Rx3-null mutants of several vertebrate species, the OV fails to evaginate (Bailey et al., 2004; Loosli et al., 2003; Mathers et al., 1997; Voronina, 2003). Additionally, Rx3-deficient cells are excluded from OV domains in fish (Loosli et al., 2003) and embryonic mouse chimeras (Medina-Martinez et al., 2009), demonstrating the crucial role of Rx3 in early retina development and morphogenesis.

Figure 1 with 2 supplements see all
Download asset Open asset
Generation of medaka fish primary pluripotent cell-derived aggregates.

(a) Scheme representing stages and timing of the medaka fish retinal development. The retinal domain is indicated in green. Establishment of the eye field within the anterior neural plate is followed by the formation of optic vesicles at 1 day post-fertilization (dpf). Optic vesicle evagination is followed by the morphogenesis of a bi-layered optic cup formed by retina surrounded by retinal pigmented epithelium (RPE) and subsequent onset of retinal differentiation at 2 dpf. By 4 dpf, the major retinal cell types – retinal ganglion cells (yellow), amacrine cells (orange), bipolar cells (red), horizontal cells (cyan), and photoreceptor cells (blue) – are generated. (b) Schematic representation of aggregate generation, its timeline and culture conditions. At day 0, primary pluripotent cells were harvested from blastula-stage medaka embryos and re-aggregated in low binding U-shape 96-well plates. At day 1, the culture media was supplemented with Matrigel. At day 2, the aggregates were transferred to a low binding culture plate and maintained in 3D suspension culture conditions in DMEM/F12 media supplemented with 5% FBS, 5% FEE, and N2 supplement. The gross morphology of the aggregates was analyzed at days 1, 2, and 3. KSR, knockout serum replacement; FBS, fetal bovine serum; FEE, fish embryonic extract. (c) Dark-field images of a blastula-stage embryo, a blastula-derived cell suspension, and re-aggregated cells and the gross morphology of aggregates at days 1, 2, and 3 after re-aggregation. (d) Optical section showing aggregate organization at day 2 visualized by immunostaining against neuroepithelium-specific markers, N-cadherin (Ncad), and acetylated tubulin (AcTub), co-stained with DAPI nuclear stain. Scale bars: 100 and 50 μm (enlargement in (d)).

Soon after evagination, two major domains are being specified within the OV (Figure 1a, 2 dpf). The distal/ventral region expresses retina-specific genes, marking prospective retinal territory, while the dorsal/outer region expresses the transcription factor Otx2, marking prospective retinal pigmented epithelium (RPE) territory (Bovolenta et al., 1997; Fuhrmann, 2010; Hatakeyama et al., 2001; Hirashima et al., 2008). Following a precisely coordinated morphogenesis (Heermann et al., 2015), the OV further forms the two-layered optic cup: the outer layer giving rise to the RPE and the inner layer forming retina populated with mitotically active retinal progenitor cells (reviewed by Casey et al., 2021; Chow and Lang, 2001; Fuhrmann, 2010). Soon after optic cup formation, retinal progenitor cells start to differentiate into seven retinal cell types that together form the structure of the adult retina: retinal ganglion cells, amacrine cells, bipolar cells, Müller glia cells, horizontal cells, and (rod, cone) photoreceptors (Young, 1985; Figure 1a, 4 dpf).

Although many aspects of retinal development are known, its mechanistic analysis in vivo is challenging due to the complex environment of an embryo, including signals from the surrounding tissues. Taking advantage of cell assembly into various organ-like structures, in vitro systems provide a unique opportunity to study the basic principles of organ formation. Already in the 1940s, experiments in amphibians showed that early embryonic tissues can assemble to form neuronal tissues when cultured under specific conditions (Simian and Bissell, 2017). Animal cap ectoderm of Ambystoma maculatum salamander embryos was shown to differentiate into forebrain and occasionally retinal tissues in the absence of inductive signals, when cultured in a simple saline solution (Barth, 1941; Holtfreter, 1944; Hurtado and De Robertis, 2007). More recently, mouse and human embryonic stem (ES) cells have been shown to form retinal tissue when aggregated and cultured under 3D suspension culture conditions (Eiraku et al., 2011; Kuwahara et al., 2015; Nakano et al., 2012). However, so far, the organoid field has been restricted to mammalian species, leaving the ability of ES cells derived from different species to assemble in retinal tissue unresolved.

ES cells derived from teleosts (medaka and zebrafish) blastula-stage embryos are able to contribute to all germ layers of chimeric embryos, indicating that the ES cell pluripotency is conserved across vertebrate species (Ho et al., 2014; Hong et al., 1998; Hong et al., 1996; Peng et al., 2019; Robles et al., 2011; Yi et al., 2010). It has recently been shown that zebrafish blastula explants can form polarized structures that recapitulate aspects of embryonic patterning and morphogenesis when cultured in the absence of yolk (Fulton et al., 2020; Schauer et al., 2020). However, under no condition reported so far, fish-derived explants or re-aggregates have formed highly organized neural structures.

Here, we demonstrate that primary embryonic pluripotent cells derived from medaka and zebrafish form aggregates of neuroepithelial identity. These aggregates adopt retinal fate following the species-specific pace of development. Strikingly, the addition of extracellular matrix components (Matrigel) allows those to form retinal organoids with multiple cell layers. By adjusting the cell-seeding density, we directed organoids to form multiple OVs and non-retinal domains mimicking architecture of a ‘head with two eyes’. Based on 3D in vivo imaging, we show that cell migration driving OV morphogenesis in these organoids is intriguingly reminiscent of the in vivo situation. Additionally, fish-derived retinal organoids retain the genetic constraints and the developmental program displayed in the embryo. Altogether, within only 4 days of culture, fish-derived primary embryonic pluripotent cells efficiently proceed through retinal differentiation, OV morphogenesis, and the onset of retinal differentiation, offering new ways to address multiple aspects of development and to systematically probe the dependence of organogenesis on variable physical environments.

Results

Generation of fish primary embryonic pluripotent cell-derived aggregates

Here, we used primary embryonic pluripotent cells derived from blastula-stage embryos of medaka (Oryzias latipes) as a source of pluripotent cells and established the conditions to generate the anterior neural structures, particularly retinal tissue. Previous studies performed with mouse and human ES cell aggregates have shown that low serum concentration in combination with 3D suspension culture support retinal specification and differentiation. In particular, the addition of extracellular matrix components such as laminin-rich Matrigel promotes retinal formation (Eiraku et al., 2011; Nakano et al., 2012). We dissociated blastula-stage embryos and used U-shaped low adhesive wells to re-aggregate 1000–2000 cells in basic media supplemented with 5% KSR (knockout serum replacement) (Figure 1b,c). The aggregation and compaction of cells was highly efficient (100%; n = ~1500 from 54 independent experiments) and resulted in a smooth and compacted morphology by day 1 (Figure 1c; Video 1). The potential of blastomeres to form aggregates was restricted to the blastula stage (1000–2000 cells) and cells derived from earlier developmental stages such as early and late morula failed to form stable aggregates or aggregated only partially (Figure 1—figure supplement 1). Thus, the exact time point of cell dissociation plays a decisive role for the success and efficiency of aggregation.

Video 1
Download asset
Aggregation and compaction of blastula-derived cells.

Time-lapse imaging of medaka-derived blastula cells going through the process of aggregation and compaction. Imaging was performed with 30 min intervals on all aggregates (n = 30) in the batch. Scale bar: 100 μm.

After 1 day of culture, media were supplemented with 2% Matrigel (Figure 1b). The aggregates then formed an organized and layered epithelium surrounding the aggregate from day 2 onward (Figure 1c) and the cells of the peripheral layer expressed neuroepithelial markers, for example, N-cadherin and acetylated tubulin, indicating that cells adopted neuroepithelial identity (Figure 1d, Figure 1—figure supplement 2). 3D imaging of aggregates at day 2 revealed a complex internal morphology, consisting of lumens and cells of neuroepithelial identity (Figure 1—figure supplement 2). Based on the higher order epithelial organization, we call the aggregates from day 2 on organoids.

Fish-derived organoids form retinal neuroepithelium under the control of Rx3

Retinal specification in vertebrates is governed by the action of retina-specific transcription factors such as Rx, Pax6, Six3, Sox2, Six6, and Lhx2 within the anterior neural plate (Li et al., 1997; Loosli et al., 1999; Zuber et al., 2003) and the formation of the OV neuroepithelium (Chow and Lang, 2001; Fuhrmann, 2010; Martinez-Morales et al., 2017). Rx genes, in teleosts represented by three paralogous genes (Rx1, Rx2, and Rx3), are the earliest genes expressed by the retinal lineage (Chuang et al., 1999; Deschet et al., 1999; Loosli et al., 2003; Loosli et al., 2001). To analyze the efficiency of retinal fate acquisition, we monitored the earliest expressed Rx paralogue, Rx3, and employed a transgenic reporter line (Rx3::H2B-GFP) in which nuclear GFP is expressed under the control of the Rx3 regulatory elements (Figure 2a,b; Rembold et al., 2006). Rx3::H2B-GFP drives the expression of H2B-GFP already in the anterior neural plate and is subsequently found in the neuroepithelium of the OV and prospective forebrain of fish embryos 1 dpf (Figure 2a). We addressed the onset of GFP expression in Rx3::H2B-GFP aggregates by time-lapse imaging over a time span of 18 hr post-aggregation (hpa) and compared it to corresponding Rx3::H2B-GFP embryos for reference (Figure 2c,d; Video 2). All aggregates started expressing GFP at 15.5 ± 0.36 hpa (n = 17) with GFP expression in control Rx3::H2B-GFP embryos starting at 13.5 ± 1.03 hr (n = 6) after the blastula stage (Figure 2c), corresponding to the onset of Rx3 expression at around 20 hrs post-fertilization (hpf) (Loosli et al., 2001). In comparison to the embryo, the organoids showed a delay of GFP expression of approximately 2–3 hr which corresponds to the time of re-aggregation (Video 1). This data shows that following our protocol, the re-aggregation of primary embryonic pluripotent cells results in an efficient acquisition of retinal fate.

Figure 2 with 4 supplements see all
Download asset Open asset
Medaka-derived organoids form retinal neuroepithelium under the control of Rx3.

(a) Schematic representation of the Rx3::H2B-GFP transgenic construct and the corresponding expression domain of GFP in the optic vesicles of a developing medaka embryo at 1 dpf. (b) Scheme of organoid generation from Rx3::H2B-GFP transgenic fish. (c) Bright-field and fluorescence images of aggregates derived from Rx3::H2B-GFP transgenic fish at 0.5, 6, and 16 hpa. (d) Analysis of the onset of GFP expression in Rx3::H2B-GFP-derived fish embryos (n = 6) and organoids (n = 17). ***p<0.001. (e) Optical sections of day 1 (before the addition of Matrigel) and day 2 organoids derived from Rx3::H2B-GFP transgenic fish stained with antibodies against N-cadherin (Ncad) and GFP, co-stained with DAPI nuclear stain. (f) Generation of Rx3KO (Rx3saGFP) line – schematic representation of the Rx3 locus with integrated saGFP-OPT cassette. An open reading frame-adjusted gene trap cassette comprising a splice acceptor and a GFP sequence (saGFP) followed by a polyA and a strong terminator sequence derived from the ocean pout (OPT; STOP) were inserted into the Rx3 locus. (g) Scheme of organoid generation from Rx3-deficient single blastulae. (h) Bright-field and fluorescence images of phenotypes of Rx3saGFP/+(Rx3 +/- heterozygote) and Rx3saGFP/saGFP (Rx3 -/- homozygote) mutants at 1 dpf and corresponding organoids at day 2. dpf, days post-fertilization; hpa, hours post-aggregation; hpf, hours post-fertilization; wt, wild type. Scale bar: 100 μm.

Video 2
Download asset
Acquisition of retinal fate within medaka blastula-derived aggregates.

Time-lapse imaging of Rx3::H2B-GFP medaka-derived blastula cells going through the process of aggregation, compaction, and acquisition of retinal fate (GFP expression). Imaging was performed with 30 min intervals on all aggregates (n = 54). Scale bar: 100 μm.

To address whether cell aggregation is a prerequisite for retinal specification or, alternatively, cells can acquire retinal fate autonomously, we followed the onset of Rx3 expression in individual blastula-stage cells. Rx3::H2B-GFP-derived primary embryonic pluripotent cells were cultured under suspension culture condition, that is, gentle rocking to prevent cell aggregation and analyzed for GFP expression 24 hr later (Figure 2—figure supplement 1a,b). Under those conditions, cells formed GFP-expressing clusters, indicating that cells acquired retinal fate even without aggregation. To further refine the analysis and address whether isolated single primary embryonic pluripotent cells give rise to GFP-expressing clusters, cells were cultured individually in the 96-well plate and imaged over time (Figure 2—figure supplement 1a,c; Video 3).

Video 3
Download asset
Acquisition of retinal fate within individual medaka – derived primary embryonic pluripotent cells.

Time-lapse imaging of medaka-derived blastula cells derived from Rx3::H2B-GFP transgenic fish. Imaging was performed with 30 min intervals. Scale bar: 100 μm.

Also under those conditions, individual cells gave rise to GFP-expressing clones, showing that cells acquire retinal fate autonomously. Similar to GFP expression in aggregated primary embryonic pluripotent cells (Figure 2c,d; Video 2), single cell-derived clones started to express GFP at about 16 hr (Video 3), indicating that the onset of retinal fate is genetically timed and independent of cellular environment.

Interestingly, Rx3 expression was initiated also in the absence of Matrigel already at day 0, indicating that the acquisition of retinal fate occurs spontaneously and independently of extracellular matrix components. Further analysis showed that GFP-expressing cells are localized in disperse manner before the addition of Matrigel, that is, on day 1 (Figure 2e). Matrigel addition triggered the formation of a neuroepithelium by day 2 and confined GFP-expressing cells to the peripheral region of the organoid. The formation of a neuroepithelium was clearly dependent on the action of Matrigel (Figure 2—figure supplement 2).

To address whether GFP-positive aggregates followed retinal specification, we analyzed the expression of retinal markers expressed subsequently during retinal specification in vivo such as Rx2, Sox2, and Lhx2. The expression of these transcription factors in the epithelium of day 2 organoids resembled their expression in the neuroepithelium of OV in the 1 dpf embryo (Figure 2—figure supplement 3), indicating that the aggregates adopted the complex retinal fate. All organoids (n = ~1500) from 52 independent experiments formed retinal neuroepithelia, indicating that retinal specification under the established conditions is robust and highly efficient.

The formation of the evaginated OV strictly depends on expression of the transcription factor Rx such that in Rx3-null mutants, OVs fail to form leading to an eyeless phenotype (Loosli et al., 2003; Loosli et al., 2001; Mathers et al., 1997; Rembold et al., 2006; Stigloher et al., 2006). Organoids derived from Rx3-null mutants can be used to investigate the dependence of the formation of the retinal neuroepithelium on Rx3. Unfortunately, null alleles resulting from mutagenesis screens do not contain traceable markers. We therefore followed a Crispr/Cas targeted gene-trapping approach into an intron of Rx3. We inserted an open reading frame-adjusted gene trap cassette comprising a splice acceptor and a GFP sequence (saGFP) for visual readout followed by a polyA and strong terminator sequence derived from the ocean pout (OPT; Clark et al., 2011, Figure 2—figure supplement 4). Using the CRISPR/Cas9 system for targeted integration, we generated traceable Rx3 mutants (Rx3saGFP; Figure 2—figure supplement 4). While heterozygotes Rx3saGFP/+ develop normally and give rise to GFP labeled OVs, homozygous-mutant embryos (Rx3saGFP/saGFP) fail to evaginate the OVs from the lateral wall of the diencephalon (Figure 2g, Figure 2—figure supplement 4). In contrast to the eyeless mutant described (Loosli et al., 2001), the mutant phenotype, that is the inability to evaginate the OV, was not temperature sensitive. To address whether fish-derived organoids retain the genetic constraints and the developmental program displayed in the embryo, we derived organoids from individual Rx3saGFP/saGFP and Rx3saGFP/+ blastulae (Figure 2g). While a retinal neuroepithelium formed in Rx3saGFP/+ (n = 9) organoids by day 2, aggregates derived from individual Rx3-mutant homozygote blastulae (Rx3saGFP/saGFP) (n = 4) did not form a neuroepithelium and aggregates were lacking regular morphology by day 2 (Figure 2g) and thus resembled organoids cultured in the absence of Matrigel. These data indicate that the formation of retinal neuroepithelium is hard-wired also in the organoid and follows the same developmental program as in vivo, also upon loss-of-function (Loosli et al., 2001; Rembold et al., 2006; Stigloher et al., 2006; Winkler et al., 2000).

We next asked whether the spontaneous acquisition of retinal fate is generally applicable to other fish species, for example, the evolutionarily distant zebrafish (Danio rerio). We used zebrafish blastulae as the source for primary embryonic pluripotent cells to establish zebrafish retinal organoids. To monitor the acquisition of retinal fate, we used the early retinal marker gene Rx3 as detailed above. We used a transient assay and injected a Rx3::H2B-GFP reporter construct into one cell stage embryos and generated cell aggregates when embryos reached blastula stage (Figure 3a).

Figure 3
Download asset Open asset
Zebrafish blastula-derived primary embryonic pluripotent cells forming retinal organoids.

(a) Scheme of aggregate generation from Rx3::H2B-GFP DNA-injected blastulae. (b) Bright-field images of zebrafish-derived aggregates at 1 hpa and 1 dpa. (c) Analysis of the onset of GFP expression in Rx3::H2B-GFP-derived zebrafish embryos (n = 6) and organoids (n = 15). **p<0.01. (d) Fluorescent images of mCherry and GFP expression domains in Rx3::H2B-GFP-injected embryos at 1 dpf and corresponding aggregates at day 1. Dashed lines indicate the embryonic retina and organoid retinal neuroepithelium, respectively. wt, wild type; dpf, days post-fertilization; dpa, days post-aggregation; hpa, hours post-aggregation. Scale bar: 100 μm.

Zebrafish cells aggregated as efficiently as medaka cells and displayed an epithelial organization already by day 1 (Figure 3b). Considering the faster early development of zebrafish compared to medaka (Iwamatsu, 2004; Kimmel et al., 1995), we asked whether the difference in embryonic development impacts on the re-aggregation and subsequent differentiation of aggregates. We performed time-lapse imaging in zebrafish aggregates over 18 hpa and addressed the onset of retinal differentiation by Rx3-driven GFP expression. Aggregates generated from zebrafish embryos formed organoids that started expressing GFP at 9.75 ± 0.21 hpa (n = 15). As in medaka, the onset of Rx3 expression in the forming organoids was delayed by 2–3 hr in comparison to the corresponding Rx3::H2B-GFP reporter embryos (Figure 3c; Video 4) that started to express GFP at 7.88 ± 0.77 hr (n = 6) after the blastula stage. The cells followed their endogenous program consistently and were only delayed by the dissociation/aggregation procedure. Thus, the formation of a neuroepithelium in aggregates reflects the species-specific difference in the relative pace of development.

Video 4
Download asset
Acquisition of retinal fate within zebrafish blastula-derived aggregates.

Time-lapse imaging of Rx3::H2B-GFP zebrafish-derived blastula cells going through the process of aggregation, compaction, and acquisition of retinal fate (GFP expression). Imaging was performed with 30 min intervals. Scale bar: 100 μm.

Rx3 reporters for retina formation (Rx3::H2B-GFP) displayed H2B-GFP expression in the developing retina at 1 dpf and aggregates formed organoids showing retinal identity as indicated by Rx3 expression in the outer cellular layer (Figure 3d) by day 1. This data shows that under the established conditions, fish primary embryonic pluripotent cells are guided to differentiate into a retinal neuroepithelium – a feature found over a wide evolutionary distant medaka and zebrafish, indicating a deep conservation of early retinal development irrespective of the embryonic environment.

Fish organoids form OV-like structures

Since primary embryonic pluripotent cells from both, medaka and zebrafish, efficiently form retinal organoids, we subsequently focused our experiments on one organism – medaka.

Development of a functional retina is accompanied by several morphological transitions that ultimately result in the formation of a functional retina. One of them is the formation of the OV. When organoids were generated by the aggregation of 1000–2000 cells (referred to as >1000 cells), the approximate number of cells in a single blastula, OV formation was not favored (Figure 4a) and instead these organoids formed a continuous neuroepithelium on their entire surface. This neuroepithelium expressed retina-specific markers (Rx2 and Rx3::H2B-GFP) all around the surface of organoid (Figure 4a–b) indicating the formation of a single retinal anlage. It has been previously reported that in some in vitro self-organizing systems, the number of interacting cells can generate a bias toward particular morphological processes (van den Brink et al., 2014; Fulton et al., 2020; Völkner et al., 2016). Thus, we reduced the cell-seeding density by 50%. Interestingly, those organoids (‘small organoids’, generated by the aggregation of <1000 cells) reproducibly displayed expression of Rx2 and Rx3 only in restricted regions (Figure 4a,b) reminiscent of the forming OV in vivo. While organoids generated by aggregation of >1000 cells established a single, uniform, domain of retinal neuroepithelium (n = 35/35), small organoids displayed a more complex and variable morphology, forming one to four isolated OV-like retinal domains (n = 92/123 one domain, n = 27/123 two domains, n = 3/123 three domains, and n = 1/123 four domains) (Figure 4c, Figure 4—figure supplements 1 and 2). The size of these isolated retinal domains in small organoids (154.5 ± 28.6 μm; n = 56) was reminiscent of the OVs in embryos at 1 dpf (162 ± 10.6 μm; n = 16) (Figure 4d).

Figure 4 with 2 supplements see all
Download asset Open asset
Medaka primary embryonic pluripotent cells form optic vesicle-like structures.

(a) Fluorescence and bright-field images of day 2 organoids produced by aggregation of >1000 cells and <1000 cells stained with anti-Rx2 antibody. (b) Analysis of the area of Rx2 (wild-type organoids stained with anti-Rx2 antibody) and Rx3 (Rx3::H2B-GFP-derived organoids stained with anti-GFP antibody) expression area (% of total organoid area) in day 2 organoids. ****p<0.0001. (c) Number of organoids forming (1–4) individual retinal regions produced by aggregation of >1000 (n = 26 for Rx2, n = 9 for Rx3) and <1000 (n = 57 for Rx2, n = 66 for Rx3) cells from nine independent experiments. (d) Size, measured as largest circumference, of the optic vesicle of 1 dpf embryos and optic vesicle-like structures formed by day 2 organoids (n = 16 embryos, n = 56 aggregates from six independent experiments). ns, non-significant. (e) Bright-field and fluorescent images of day 2 Rx3::H2B-GFP organoids stained with anti-GFP antibody. Optical sections of an organoid (day 2) (n = 9/10) and an embryo (1 dpf) stained with anti-Rx2 and anti-Otx2 antibodies. (f) Maximal projection of day 2 organoids and 1 dpf embryo generated from Rx3::H2B-GFP transgenic or wild-type blastulae and stained with neural tissue-specific anti-Sox2 (n = 12/12) and anti-Otx2 (n = 10/10) antibodies, co-stained with anti-Rx2 and DAPI nuclear stain. hpa, hours post-aggregation; dpf, days post-fertilization. Scale bars: 100 μm.

Interestingly, retinal domains formed in small organoids displayed OV-like morphology with specifically localized expression of Rx2 and Otx2 transcription factors, marking prospective retinal and RPE territories respectively (Figure 4e) and indicating further compartmentalization of the organoid-derived OV.

Besides the retinal tissue, small aggregates also established an Rx2-negative, non-retinal domain, expressing the general neuronal marker Sox2 (Pevny and Placzek, 2005; Wegner and Stolt, 2005) and the forebrain marker Otx2 (Acampora et al., 1995; Simeone et al., 1993; Tian et al., 2002), strongly indicating that the neighboring tissue adopted anterior neuronal identity (Figure 4f). In a fraction (22%) of organoids, intriguingly, retinal and non-retinal domains showed an arrangement, resembling a ‘head with two eyes’ (n = 27/123) (Figure 4e–f, Figure 4—figure supplements 1 and 2).

To address the mode of OV formation in fish organoids, we employed live imaging of Rx3::H2B-GFP-derived organoids from day 1 to day 2 of development. Figure 5a displays the development of the organoid forming two opposing OV structures. At the initial time point (24 hpa), GFP-positive cells were distributed throughout the whole volume of the organoid indicating that cells acquired retinal identity individually. Single cell tracking and analysis of cell displacement showed that initially cells behaved as freely diffusing particles (Figure 5b). This behavior changed at about 34 hpa, when GFP-positive cells showed directed migration – non-linear growth of mean squared displacement – with an estimated speed of 0.413–0.479 μm/min resulting in the morphogenesis of the OV-like structures by 40 hpa (Figure 5a,b; Video 5). A similar transition from linear to nonlinear diffusion, due to increasing density of nuclear packing, was recently demonstrated in vivo in the growing retina of 24 hpf zebrafish (Azizi et al., 2020). Directionality analysis of single cell tracks (Figure 5c) showed that 51% of the Rx3-expressing cells were moving from the center to the periphery (outward) of the organoid and contributed to the formation of the OV (see dark blue, red, orange, and yellow tracks in Figure 5a). About 49% of the initially laterally located cells moved in the opposite direction – toward the interior of the organoid (inward) (light green track in Figure 5a). These outward- and inward-directed cell movements were highly reminiscent of the migratory behavior of Rx3-expressing cells during OV formation in medaka fish (Rembold et al., 2006). Notably, a similar behavior was observed in the organoids forming single and multiple OVs (Video 6). Altogether these results show that the cell behavior in fish organoids closely recapitulates the corresponding behavior during development in vivo.

Figure 5
Download asset Open asset
Dynamics of the optic vesicle formation in medaka organoids.

(a) In vivo time-lapse images acquired with MuVi SPIM of optic vesicle-like structure evagination in Rx3::H2B-GFP retinal organoids. From top to bottom: maximal projection of quasi bright field, GFP, and tracks of exemplary cells. (b) Normalized velocity autocorrelation for all tracks (n = 4600; from two optic vesicles of one organoid). (c) Directional histogram of tracks (n = 4600; from two optic vesicles of one organoid) separated for inward and outward movements. hpa, hours post-aggregation; dpf, days post-fertilization. Scale bar: 100 μm.

Video 5
Download asset
Optic vesicle formation within blastula-derived retinal organoid.

Time-lapse imaging of optic vesicle-like structure evagination in the medaka Rx3::H2B-GFP retinal organoid. Maximal projection of quasi bright field, GFP, and tracks of exemplary cells. Imaging was performed with 15 min intervals. Scale bar: 100 μm.

Video 6
Download asset
Example of organoids forming single and multiple optic vesicle-like structures.

Time-lapse imaging of optic vesicle-like structure evagination in medaka Rx3::H2B-GFP retinal organoids. Maximal projection of GFP and tracks of exemplary cells. Imaging was performed with 15 min intervals. Scale bar: 100 μm.

It is worth mentioning that the organization of retina-committed cells into the OV neuroepithelium is strictly dependent on the presence of extracellular matrix proteins as Rx3-expressing cells attempt but fail to form OV structures in the absence of Matrigel (Video 7). Consistently, laminin-1 (the main component of the Matrigel) is highly enriched specifically at the basal surface of the forming OV and has been reported to play an essential role in cell polarization and maintenance of neuroepithelial cell morphology during OV evagination (Ivanovitch et al., 2013).

Video 7
Download asset
Extracellular matrix controls formation of optic vesicle neuroepithelium.

Time-lapse imaging of optic vesicle-like structure evagination in medaka Rx3::H2B-GFP retinal organoids in presence or absence of Matrigel (extracellular matrix component) shown as maximal projection of GFP and bright field. Imaging was performed with 30 min intervals. Scale bar: 100 μm.

Fish organoids show onset of retinal differentiation

Retinal specification is followed by the process of differentiation which leads to the generation of seven retinal cell types (retinal ganglion cells, amacrine cells, horizontal cells, bipolar cells, rod and cone photoreceptors, and Müller glia cells) organized in three nuclear layers. One of the first hallmarks of retinal differentiation is the expression of the transcription factor Atoh7 (Brown et al., 1998; Del Bene et al., 2007; Kay et al., 2001). Atoh7-positive progenitors have been found to give rise to retinal ganglion cells, amacrine cells, horizontal and photoreceptor cells during fish retinal development (Poggi et al., 2005). We monitored Atoh7 expression in medaka-derived retinal organoids, using an Atoh7::EGFP transgenic line (Del Bene et al., 2007; Figure 6a–c). In embryos, EGFP expression was located specifically in the differentiating retina at 2 dpf (Figure 6c). Day 3 organoids generated from Atoh7::GFP blastulae (generated by the aggregation of <1000 cells) showed specific morphology with EGFP-negative non-retinal and EGFP-expressing retinal domains (Figure 6b; Video 8), indicating that the retinal differentiation program was successfully initiated. Furthermore, non-retinal regions proceeded through the process of neuronal differentiation as indicated by expression of HuC/D (Figure 6c; Video 8), a marker of early post-mitotic neurons (Good, 1995; Kim et al., 1996; Park et al., 2000).

Figure 6 with 2 supplements see all
Download asset Open asset
Medaka-derived organoids show onset of retinal differentiation.

(a) Scheme of organoid generation from Atoh7::EGFP transgenic fish. (b) Fluorescent images of day 3 Atoh7::EGFP organoids (generated by aggregation of <1000 cells). (c) Optical sections and maximal projections showing EGFP expression in the eye of the developing embryo at 2 dpf and the retinal organoid at day 3 co-stained with antibody against HuC/D and DAPI. (d) Optical sections showing expression of HuC/D (amacrine and ganglion cells), Otx2 (bipolar and photoreceptor cells), and Prox1 (horizontal cells) in day 4 organoid. (e) Sketch of the arrangement of cellular layers in the organoid and the embryonic retina. dpf, days post-fertilization. Scale bar: 100 μm.

Video 8
Download asset
Organoids show complex morphology with the onset of retinal differentiation.

3D rendering of medaka Atoh7::EGFP day 3 organoid stained with anti-GFP (green), anti-β-catenin (cyan), and anti HuC/D (magenta) antibodies.

Furthermore, retinal domains of day 4 organoids consisted of differentiating retinal neurons organized in multiple layers of cells differentiating toward amacrine and ganglion (expressing HuC/D) (Kay et al., 2001; Link et al., 2000; Park et al., 2000), photoreceptor and bipolar (expressing Otx2) (Fossat et al., 2007; Glubrecht et al., 2009; Koike et al., 2007; Nishida et al., 2003), as well as horizontal (expressing Prox1) (Dyer et al., 2003) cell lineages (Figure 6d). Interestingly, the arrangement of the respective cellular layers was inverted when compared to normal retinal organization (Figure 6e). This can be attributed to the morphological differences, particularly the OV to optic cup transition, between retinal morphogenesis in vivo and organoid. In organoids, retinal differentiation is initiated without the optic cup formation, most probably leading to the inverted retinal organization.

Organoids generated by the aggregation of >1000 cells, which form retinal neuroepithelium around entire organoid surface, proceed through the process of retinal differentiation as well (Figure 6—figure supplement 1), suggesting that the OV formation is not a prerequisite for the retinal differentiation.

Although the organoid-derived OVs compartmentalize to prospective retinal and RPE domains by day 2 (Figure 4e), none of the organoids analyzed between day 3 and day 4 (n = 88, data not shown) showed significant level of pigmentation, the hallmark of RPE differentiation. This indicates that, similarly to human and mouse retinal organoid cultures (Eiraku et al., 2011; Kuwahara et al., 2015; Nakano et al., 2012), RPE fate was not favored under the established culture conditions. Previous studies showed that RPE fate determination in vivo, as well as in retinal organoids, is promoted by the activation of the Wnt/β-catenin signaling pathway (Eiraku et al., 2011; Fujimura et al., 2009; Kuwahara et al., 2015; Nakano et al., 2012; Westenskow et al., 2009). We thus promoted the Wnt/β-catenin pathway by treatment of the organoids from day 1 on with the GSK3β inhibitor CHIR-99021. Analysis at day 2 and day 4 indicated the presence of RPE (Figure 6—figure supplement 2). Compared to control (DMSO-treated) organoids, CHIR-99021-treated day 2 organoids showed a RPE-specific Otx2 expression on the expense of the retina-specific Rx2 expression (Figure 6—figure supplement 2b). By day 4, retinal domains in treated organoids (13/13) became pigmented with most cells expressing Otx2 (Figure 6—figure supplement 2c), demonstrating that those cells acquired RPE fate.

Taken together, our findings show that primary embryonic pluripotent cells derived from blastula-stage fish embryos form aggregates with the retinal differentiation potential. By adjusting numbers of aggregating cells and the addition of Matrigel, these aggregates are efficiently directed to the retinal fate, the OV morphogenesis, and the retinal differentiation.

Discussion

During embryogenesis, individual cells interact with each other and the environment to establish 3D structures of high complexity – tissues, organs, and ultimately the entire embryo. In the absence of immediate access to mammalian embryos, development of 3D in vitro cultures provides a unique platform to study development and the roots of pathologies under tightly controlled conditions (Hofer and Lutolf, 2021; Kretzschmar and Clevers, 2016; Takebe and Wells, 2019). In less than a decade, retinal organoids from human and mouse pluripotent stem cells have been shown to recapitulate in vivo retinogenesis, model several retinal diseases, and serve as a potential source for transplantation therapy (Artero Castro et al., 2019; Brooks et al., 2019; Fligor et al., 2018; Gao et al., 2020; Kruczek and Swaroop, 2020; Kuwahara et al., 2015; Völkner et al., 2016). While human retinal organoids are the most desired, they also take the longest time to develop and their systematic analysis is limited by technical challenges in genetic manipulations and, last not least, by the lack of direct comparison to the in vivo biological processes. In developmental biology, these challenges are tackled by use of other, more accessible genetic models. Similarly, the derivation of retinal organoids from non-mammalian models, for example medaka and zebrafish, may provide an easier and more immediate way to address different aspects of eye development. Although fish blastula-derived cells have been used to generate ES cultures (Ho et al., 2014; Hong et al., 1998; Hong et al., 1996; Peng et al., 2019; Robles et al., 2011; Yi et al., 2010), our current understanding of their ability to differentiate and assemble into tissues is very limited. Zebrafish blastula explants have been found to form polarized aggregates capable of specification of all three germ layers (Fulton et al., 2020; Schauer et al., 2020), resembling mammalian ES cell-derived gastruloids (Beccari et al., 2018; van den Brink et al., 2014; Moris et al., 2020; Turner et al., 2017). Besides that, fish primary ES cells have been directed to differentiate to cardiomyocytes forming beating cell sheets with contractile kinetics and electrophysiological features when cultured under defined chemical and mechanical conditions (Xiao et al., 2016). However, under no condition reported so far, fish-derived 3D in vitro cultures have formed highly organized neural structures.

Here, we demonstrate that similar to mammalian ES cells, primary embryonic pluripotent cells derived from fish can organize organizeinto anterior neural structures, particularly the retina, consistent with the proposed default differentiation program followed by pluripotent stem cells (Eiraku et al., 2011; Hemmati-Brivanlou and Melton, 1997; Kuwahara et al., 2015; Levine and Brivanlou, 2007; Nakano et al., 2012; Tropepe et al., 2001; Turner et al., 2014). In addition, inter-species blastula cell transplantation experiments between zebrafish and medaka result in ectopic retina formation composed exclusively of donor cells (Fuhrmann et al., 2020), further reflecting the intrinsic tendency of primary embryonic pluripotent cells to form retinal structures. Our results show that retinal fate specification in fish-derived organoids is executed spontaneously.

We found that the starting size (cell number) of the initial aggregate crucially impacts on the morphogenesis of the resulting organoid. While aggregates generated from <1000 cells differentiated into anterior neuroectoderm that subsequently underwent morphogenesis, bulging out retinal primordia that closely resembled OVs; aggregates generated from >1000 cells entirely differentiated into a single giant retina. This appeared counterintuitive, since the ‘classical’ reaction-diffusion model of pattern formation (Meinhardt, 2012; Turing, 1952) predicts that the number of repeating patterns – OVs – will increase proportionally with tissue size. For fish retinal organoids, the effect is however inverted, offering fish retinal organoids as a model to address alternative pattering scenarios.

Remarkably, re-aggregation of dissociated fish primary pluripotent cells was fast (about 2–3 hr) and highly efficient (100%) when derived from blastula-stage embryos. With OV formation by day 2 and the onset of differentiation by day 3, fish organoids proceed through development almost 10 times faster than their murine counterparts (Eiraku et al., 2011; Kuwahara et al., 2015; Nakano et al., 2012). The formation of fish retinal organoids followed the in vivo developmental timing, reflecting the difference in relative speed of development of the corresponding species.

This rapid development in combination with the ease of highly efficient genome editing (Gutierrez-Triana et al., 2018; Prykhozhij and Berman, 2018; Stemmer et al., 2015), small size, and full transparency granting direct comparison to in vivo processes, pose fish-derived organoids as a rapid and efficient prototyping platform for systematic drug tests in disease models, which then can be translated to more demanding mammalian systems.

Although eye development is highly stereotypic across vertebrates and the players involved are highly conserved, there are certain species-specific features that are also manifested in organoids. One such example is the cellular mechanism of OV evagination. In mammalian eye development, paralleled in mouse and human retinal organoid cultures, the OV is formed by a set of coordinated collective epithelial cell movements driving vesicle evagination from the initially formed neuroepithelium (Chauhan et al., 2015; Eiraku et al., 2011; Martinez-Morales et al., 2017). In teleosts, in contrast, OV formation is driven by the individual migration of Rx3-expressing retinal progenitor cells (Rembold et al., 2006). Fish-derived retinal organoids exhibit the same mechanism of OV formation, as recapitulated by the migratory trajectories and individual movements of Rx3-expressing cells (Figure 4g,j). Our analysis indicates that this fundamental property does not depend on the environment (embryo vs. culture) and rather reflects a ‘hard-wired’ feature of the OV formation.

While the concept of evolutionary conservation allows to deduce fundamental mechanisms by comparison of embryos and the crucially contributing genes, it is not immediately transferable to the comparison of organoids from different species. Selection acts only on development and consequently relies on the entire context of the developing embryo. Morphological and molecular resemblances of organoids to the corresponding organs may indicate a self-organizing capacity that is irrespective of the environment or organismal constraints. The comparison of organoids derived from a wide range of evolutionarily diverse species provides the opportunity to systematically address and distinguish intrinsic and extrinsic constraints contributing to self-organization and organ morphogenesis, and to ultimately tackle the question of the ‘conservation’ of cellular self-organization, a feature that has apparently not been selected for in evolution.

Materials and methods

Key resources table
Reagent type
(species) or resource		DesignationSource or referenceIdentifiersAdditional information
Gene (Oryzia latipes)rx3EnsemblENSORLG00000027320–
Strain, strain background (Oryzia latipes)CabLoosli et al., 2000––
Strain, strain background (Danio rerio)ABZIRCZFIN: ZBD-GENO-960809–7–
Genetic reagent (Oryzia latipes)Atoh7::EGFPDel Bene et al., 2007––
Genetic reagent (Oryzia latipes)Rx3::H2B-GFPRembold et al., 2006––
Genetic reagent (Oryzia latipes)Rx3saGFPThis study–Insertion of splice acceptor followed by GFP (saGFP) andocean pout polyA terminator (OPT) in the first intron of Rx3 gene
AntibodyAnti-Otx2 (polyclonal goat IgG)R&D SystemsCat#:AF1979 RRID:AB_2157172IHC (1:200)
AntibodyAnti-Rx2 (polyclonal rabbit IgG)In-house;
Reinhardt et al., 2015		–IHC (1:200)
AntibodyAnti-Lhx2 (polyclonal rabbit IgG)GeneTexCat#:GTX129241 RRID:AB_2783558IHC (1:500)
AntibodyAnti-β-catenin (polyclonal rabbit IgG)AbcamCat#:Ab6302 RRID:AB_305407IHC (1:500)
AntibodyAnti-acetylated tubulin (monoclonal mouse IgG2b)MerckCat#:T7451 RRID:AB_609894IHC (1:200)
AntibodyAnti-Sox2 (polyclonal rabbit IgG)GeneTexCat#:GTX124477 RRID:AB_11178063IHC (1:200)
AntibodyAnti-GFP (polyclonal chicken IgY)Thermo Fisher ScientificCat#:A10262 RRID:AB_2534023IHC (1:500)
AntibodyAnti-N-cadherin (monoclonal rabbit IgG)AbcamCat#:ab76011 RRID:AB_1310479IHC (1:200)
AntibodyAnti-Prox1 (polyclonal rabbit IgG)MerckCat#:AB5475 RRID:AB_177485IHC (1:500)
AntibodyAnti-HuC/HuD (monoclonal mouse IgG2b)Thermo Fisher ScientificCat#:A21271 RRID:AB_221448IHC (1:200)
Recombinant DNA reagentpGBT-RP2
(plasmid)		Clark et al., 2011–Ocean pout polyA terminator (OPT) sequence-carrying plasmid
Recombinant DNA reagentpGGEV_5_linker (plasmid)AddgeneRRID:Addgene_49285–
Recombinant DNA reagentpGGDestSC-ATG (plasmid)AddgeneRRID:Addgene_49322–
Recombinant DNA reagentpGGEV_4_linker (plasmid)AddgeneRRID:Addgene_49284–
Recombinant DNA reagentpGGEV_7’_linker (plasmid)AddgeneRRID:Addgene_49293–
Recombinant DNA reagentRx3::H2B-GFP
(plasmid)		Rembold et al., 2006––
Recombinant DNA reagentpGGD(saGFP-OPT-MCS+2) (plasmid)This study–Splice acceptor and GFP (saGFP); ocean pout polyA terminator (OPT) sequence-carrying plasmid
Sequence-based reagentRx3_FThis studyPCR primerTCCTTTTTAGACAAATGTGGCTCC
Sequence-based reagentGFP_RThis studyPCR primerGCTCGACCAGGATGGGCA
Sequence-based reagentpDest_FThis studyPCR primerATTACCGCCTTTGAGTGAGC
Sequence-based reagentRx3_RThis studyPCR primerGACAGGTATCCGGTAAGCGG
Sequence-based reagentrx3_T1This studysgRNAAGCAGAGCGCGCAAAGAACC[AGG]
Sequence-based reagentrx3_T2This studysgRNAAGCGCGCAAAGAACCAGGCA[GGG]
Peptide, recombinant proteinQ5 High-Fidelity DNA PolymeraseNEBCat#:M0491L–
Peptide, recombinant proteinI-SceI meganucleaseNEBCat#:R0694L–
Commercial assay or kitInnuPREP DOUBLEpure KitAnalyticJenaCat#:845-KS-5050250–
Chemical compound, drugCHIR-99021MerckCat#:SML1046–
Software, algorithmGeneious R8.1Biomatters––
Software, algorithmCCTopStemmer et al., 2015RRID:SCR_016890–
Software, algorithmFiji distribution of ImageJSchindelin et al., 2012RRID:SCR_002285–
Software, algorithmHyper stack generator (Fiji plugin)10.5281/zenodo.3368134–Generate hyperstack for images acquired with the Acquifer machine
Software, algorithmElastixWrapper (Fiji plugin)Tischer, 2019; Klein et al., 2010––
Software, algorithmMATLABMathWorksRRID:SCR_001622–
Software, algorithmMSD analysisTinevez and Herbert, 2020–MATLAB script for MSD analysis
Software, algorithmDirectionality analysisThis studyhttps://github.com/VeneraW/DirectionalityANalysisOrganoids (copy archived atswh:1:rev:31a89aead3e83ac774e13c0161e44deafce58f05; Zilova, 2021)–MATLAB script for organoid directionality analysis
Software, algorithmTrackMateTinevez et al., 2017––
Software, algorithmStatannot package, Jupyter notebookhttps://github.com/webermarcolivier/statannot, Weber, 2020––
OtherMatrigelCorningCat#:356238–
OtherDRAQ5Thermo Fisher ScientificCat#:65-0880-921:1000
OtherCellMaskThermo Fisher ScientificCat#:C100451:1000
OtherTissue Freezing MediaLeicaCat#:14020108926–

Fish handling and maintenance

Request a detailed protocol

Medaka (O. latipes) and Zebrafish (D. rerio) stocks were maintained according to the local animal welfare standards (Tierschutzgesetz §11, Abs. 1, Nr. 1, husbandry permit AZ35-9185.64/BH, line generation permit number 35–9185.81/G-145/15 Wittbrodt). The following medaka lines were used in this study: Cab strain as a wild type (Loosli et al., 2000), Rx3::H2B-GFP (Rembold et al., 2006), Atoh7::EGFP (Del Bene et al., 2007), Rx3saGFP (this study). The following zebrafish lines were used in this study: AB zebrafish as a wild type (ZIRC, ZFIN: ZBD-GENO-960809–7).

Cloning of pGGD(saGFP-OPT-MCS) and generation of Rx3saGFP knock-in line

Request a detailed protocol

The ocean pout polyA terminator (OPT) sequence was released from plasmid pGBT-RP2 (Clark et al., 2011) via restriction digest with BfaI-FD (Thermo Fisher Scientific), ligated into pGGEV_5_linker (Addgene #49285), and fused into pGGDestSC-ATG (#49322 Addgene) using the Golden GATEway cloning system (Kirchmaier et al., 2013) and the following sequences: target site sequence of sgRNA GFP_T1 (Stemmer et al., 2015) inserted in pGGEV_1, a multiple cloning site in pGGEV_2, a strong AD splice acceptor (Centanin et al., 2011) fused with a GFP variant not targeted by sgRNA GFP_T1 (Stemmer et al., 2015) in pGGEV_3, pGGEV_4_linker (#49284 Addgene), a multiple cloning cassette in pGGEV_6 and pGGEV_7’_linker (#49293 Addgene). Adjustment of open reading frame following the splice acceptor was accomplished via Q5 (NEB, Cat#:M0491L) mutagenesis by inserting a single or two nucleotides (+1 or +2, respectively) to yield gene trap vectors for all three forward frames.

sgRNA rx3_T1 (AGCAGAGCGCGCAAAGAACC[AGG], PAM in brackets) and rx3_T2 (AGCGCGCAAAGAACCAGGCA[GGG], PAM in brackets) were designed using CCTop and cloned and transcribed as described previously (Stemmer et al., 2015). The saGFP-OPT-MCS+2 cassette was inserted into the first intron of Rx3 by non-homologous end joining via microinjection into the cytoplasm of one cell stage medaka zygotes. Injection mix contained 15 ng/µl of each sgRNA rx3_T1, rx3_T2 and GFP_T1, 150 ng/µl Cas9 mRNA and 10 ng/µl pGGD(saGFP-OPT-MCS+2) in nuclease-free water. Embryos were raised and maintained at 28°C in 1× Embryo Rearing Medium (ERM, 17 mM NaCl, 40 mM KCl, 0.27 mM CaCl2, 0.66 mM MgSO4, 17 mM HEPES) and screened for ocular GFP expression 1 dpf on a Nikon SMZ18. Genotyping-PCR was performed with Q5 High-Fidelity DNA Polymerase (NEB) with 98°C initial denaturation for 2 min, followed by 30 cycles of 98°C denaturation for 20 s, 66°C annealing for 30 s, 72°C extension for 25 s. Primers used to amplify 5’ integration: Rx3_F 5’-TCCTTTTTAGACAAATGTGGCTCC, GFP_R 5’-GCTCGACCAGGATGGGCA; 3’ integration: pDest_F 5’-ATTACCGCCTTTGAGTGAGC, Rx3_R 5’-GACAGGTATCCGGTAAGCGG. Following gel electrophoresis, amplicons were gel-purified (InnuPrep, AnalyticJena), sequenced (Eurofins Genomics), and analyzed using Geneious R8.1 (Biomatters).

Injection in zebrafish embryos

Request a detailed protocol

Rx3::H2B-GFP DNA (10 ng/μl) was co-injected with Meganuclease (I-SceI) (NEB, Cat#:R0694L) and mCherry mRNA (10 ng/μl) into the cytoplasm of one cell stage zebrafish embryos as previously described (Thermes et al., 2002). Only morphologically intact and brightest mCherry expressing embryos were used for aggregate formation.

Generation of organoids

Request a detailed protocol

For medaka, blastula-stage (6 hpf) embryos (Iwamatsu, 2004) were collected, dechorionated using hatching enzyme, and washed in ERM (17 mM NaCl, 40 mM KCl, 0.27 mM CaCl2, 0.66 mM MgSO4, 17 mM HEPES). Cell mass was separated from yolk, washed three times with sterile PBS (Thermo Fisher Scientific, Cat#:10010023), and dissociated by gentle pipetting with 200 μl pipet tip. Cell suspension was pelleted (180 × g for 3 min) and re-suspended in differentiation media: GMEM (Glasgow’s Minimal Essential Medium, Gibco Cat#:11710035), 5% KSR (Gibco Cat#:10828028), 0.1 mM non-essential amino acids, sodium pyruvate, 0.1 mM β-mercaptoethanol, 50 U/ml penicillin-streptomycin to desired cell density, that is, for aggregation of >1000 cells, 10–20 cells/µl and for aggregation of <1000 cells, 5–8 cells/µl. Cell suspension (100 μl per well – per single aggregate) was transferred to low binding 96-well plate (Nunclon Sphera U-Shaped Bottom Microplate, Thermo Fisher Scientific) and centrifugated (for 3 min at 180 × g) to speed up cell aggregation. The aggregates were incubated over night at 26°C in the incubator without CO2 control. The following day (day 1) aggregates were washed with differentiation media, transferred to fresh wells, and Matrigel (Corning, Cat#:356238) was added to the media to a final concentration of 2%. From day 1 onward, aggregates were incubated at 26°C and 5% CO2. Alternatively, when CO2 control was not possible, media was supplemented with 20 mM HEPES, pH=7.4. From day 2 onward, organoids were kept in DMEM/F12 supplemented with 5% FBS (Sigma Aldrich, Cat#:12103C), 5% FEE (fish embryonic extract) (https://zfin.org/zf_info/zfbook/chapt6.html), 20 mM HEPES pH=7.4, N2 supplement (Gibco, Cat#:17502048) and 50 U/ml penicillin-streptomycin.

For zebrafish cell aggregation, blastula-stage embryos (at high stage) were processed according to the same protocol as for medaka but re-suspended in Leibowitz’s L-15 media (Gibco, Cat#:11415064) supplemented with 5% KSR and 50 U/ml penicillin-streptomycin. Matrigel (Corning, Cat#:356238) was added at 2 hpa.

Treatment of medaka organoids with GSK3 inhibitor CHIR-99021

Request a detailed protocol

After 6 hr of incubation with Matrigel at 26°C, day 1 aggregates were washed and supplemented with differentiation media containing 5 μM CHIR-99021 (Merck, Cat#:SML1046) in DMSO or DMSO only and incubated till the point of analysis (day 2 or day 4).

Fluorescent labeling

Request a detailed protocol

To stain the plasma membranes, organoids were incubated in CellMask Orange Plasma Membrane Stain (Thermo Fisher Scientific, Cat#:C10045; 1:1000) for 30 min at 26°C, washed and fixed in 4% PFA.

Immunohistochemistry was performed as previously described (Inoue and Wittbrodt, 2011) with slight modifications. Embryos and organoids were fixed in 4% PFA overnight at 4°C, and washed in PTW (PBS with 0.05% Tween20). For sectioning, fixed samples were cryopreserved in 30% (w/w) sucrose over night at 4°C. Samples were equilibrated in the 1:1 mixture of sucrose and Tissue Freezing Media (Leica; Cat#:14020108926) overnight at 4°C, frozen in Tissue Freezing Media and sectioned to 10–12 μm. Sections were re-fixed for 10 min in 4% PFA, washed in PTW, blocked for 2 hr in 10% BSA, incubated with primary antibody (anti-Otx2, R&D Systems, Cat#:AF1979) overnight and washed 5 times 10 min in PTW. After 2 hr incubation with secondary antibody, samples were washed in PTW, co-stained with DAPI nuclear stain (1:500), mounted in 60% glycerol, and imaged with Leica Sp8 confocal microscope.

For whole-mount staining, fixed samples were heated in 50 mM Tris-HCl at 70°C for 15 min, permeabilized 15 min in acetone at –20°C, and blocked in 10% BSA in PTW for 1 hr. Samples were incubated with primary antibody (1:200) overnight (for Otx2 antibody incubation was prolonged to 3 days) at 4°C. The following antibodies were used: rabbit anti-Rx2 (Reinhardt et al., 2015), rabbit anti-Lhx2 (GeneTex, Cat#:GTX129241), rabbit anti-β-catenin (Abcam, Cat#:Ab6302), mouse anti-acetylated tubulin (Merck, Cat#:T7451), rabbit anti-Sox2 (GeneTex, Cat#:GTX124477), chicken anti-GFP (Thermo Fisher Scientific, Cat#:A10262), goat anti-Otx2 (R&D Systems, Cat#:AF1979), rabbit anti-N-cadherin (Abcam, Cat#:ab76011), rabbit anti-Prox1 (Merck, Cat#:AB5475), and mouse anti-HuC/D (Thermo Fisher Scientific, Cat#:A21271). Samples were washed six times 10 min in PTW, incubated with secondary antibody (1:750) (Invitrogen) with DAPI (1:500) or DRAQ5 (Thermo Fisher Scientific, Cat#:65-0880-92; 1:1000) nuclear stain overnight at 4°C and washed five times 10 min in PTW. All samples were mounted in 1% low melting agarose in PTW and imaged with Leica Sp8 confocal, Acquifer, SPIM microscopes.

Imaging

Request a detailed protocol

Fixed organoid and embryonic samples were imaged with Leica Sp8 confocal microscope. Gross morphology of embryos and organoids was assessed by Nikon SMZ18 and Leica DMi8 microscope.

Time-lapse imaging of aggregation and Rx3 expression profiles in organoids (Videos 1–4) was performed on the ACQUIFER Imaging Machine (ACQUIFER Imaging GmbH, Heidelberg, Germany) (Pandey et al., 2019). Aggregates, single cells, or dechorionated embryos were loaded into 96-well plates and placed in the plate holder of the ACQUIFER machine at 26°C for medaka and 28°C for zebrafish. For each well plate, a set of 10 z-slices (75 µm step size) was acquired in the bright field (50% LED intensity, 50 ms exposure time) and 470 nm fluorescence (50% LED excitation source, FITC channel, 200 ms exposure time) channels with a 4× NA 0.13 objective (Nikon, Düsseldorf, Germany). Imaging was performed over 20 hr of organoid development with 30 min intervals.

3D imaging of fixed Atoh7::EGFP-derived organoids (Video 8) and N-cadherin and DRAQ5-stained day 2 organoids (Figure 1—figure supplement 2) was performed on multiview selective-plane illumination MuVi SPIM Multiview light-sheet microscope (Luxendo Light-sheet, Bruker Corporation) (Krzic et al., 2012). The organoids were mounted in 1% low melting agarose (StarPure Low Melt Agarose, Cat#N3103-0100, StarLab GmbH) inside an FEP tube (Karl Schupp AG) fixed on glass capillaries. Four volumes (two cameras and two rotation angles) were acquired with the 25× detection setup (Caroti et al., 2018), 1 µm z step size for three channels: GFP (Atoh7) 488 nm excitation laser at 40% intensity, 525/50 nm emission filter, 600 ms exposure time; RFP (β-catenin) 561 nm excitation laser at 40% intensity, 607/70 nm emission filter, 600 ms exposure time; far Red (HuC/D) 642 nm excitation laser at 50% intensity, 579/40 nm emission filter, 1000 ms exposure time. The volumes were fused with Luxendo Image fusion software. Fused volumes were visualized by 3D rendering with Amira 6.2 Main software.

Live imaging of Rx3::H2B-GFP-derived organoids (Video 5) was performed on the 16× detection MuVi SPIM Multiview light-sheet microscope (Luxendo Light-sheet, Bruker Corporation). To assure normal development of organoids, the FEP tube was partially filled with 2% low melting agarose. After solidification of agarose, the tube was further filled with differentiation media and an organoid was scooped inside the tube. The tube was positioned vertically such that the organoid fell onto the solid support of agarose. Two volumes with 1.6 µm z step size for two channels, that is, GFP 488 nm excitation laser at 10% intensity, 525/50 nm emission filter, 50 ms exposure time and quasi bright field 561 nm excitation laser at 5% intensity, 568LP nm emission filter, 100 ms exposure time, were acquired every 15 min.

Live imaging of Rx3::H2B-GFP-derived organoids (Videos 6 and 7) was performed with Sp8 confocal microscope (Leica). Organoids were imaged from day 1 to day 2 at room temperature directly in low binding 96-well plate (Nunclon Sphera U-Shaped Bottom Microplate, Thermo Fisher Scientific). The volume of the organoids was acquired with 2.24 µm z step size using 488 nm excitation laser at 10% intensity in 30 min intervals.

Quantitative analysis

Request a detailed protocol

To analyze the dynamics of retinal cell fate acquisition, for each time point, a maximum projection for fluorescence and a single focused bright-field image were calculated using Fiji distribution of ImageJ (Schindelin et al., 2012) and hyper stack plugin for ACQUIFER machine (https://doi.org/10.5281/zenodo.3368134). The onset of GFP expression was determined by detection of maximum in first derivative of GFP expression. For better representation, the sum of fluorescence intensity for each time point was normalized from 0 to 1.

Analysis of gene expression in <1000 cell and >1000 cell organoids was performed on images acquired with Leica DMi8 microscope. Automatic thresholding (minimum) on bright channel was used to determine area of an organoid and on fluorescence image to identify areas corresponding to Rx3 or Rx2 expression.

To analyze cell migration behavior in organoids of day 2, time-lapse volumes acquired with SPIM were registered with ElastixWrapper for Fiji (Tischer, 2019; Klein et al., 2010) using Euler transformation, 1000 iterations, full data points. The cells within organoids were tracked with TrackMate (Tinevez et al., 2017), using simple linear tracker. The analysis of tracks’ directionality was performed with custom-written MATLAB script (https://github.com/VeneraW/DirectionalityAnalysisOrganoids; Zilova, 2021). The MSD analysis including velocity autocorrelation was performed with MATLAB tool for analysis of particle trajectories (Tinevez and Herbert, 2020). For demonstration purposes, we chose a few representative tracks with high total displacement (total distance traveled) and duration (time required for travel).

Statistical analysis and plots were prepared with Jupyter notebook (Perez and Granger, 2007), using statannot package (https://github.com/webermarcolivier/statannot) to compute statistical tests (Wilcoxon-Mann-Whitney) and add statistical annotations. For all figures, ns 0.05<p<1, *0.01<p<0.05, **0.001<p<0.01, ***0.0001<p<0.001, ****p<0.0001. Figures were assembled with Adobe Illustrator CS6.

Data availability

Raw datasets from time-lapse experiments (Figure 5, Video 5) are deposited in publicly available repository HeiData (https://heidata.uni-heidelberg.de/).

The following data sets were generated
    1. Weinhardt V
    2. Zilova L
    3. Tavhelidse T
    4. Schlagheck C
    5. Thumberger T
    6. Wittbrodt J
    (2021) Heidelberg University (HeiData)
    Fish primary embryonic pluripotent cells assemble into retinal tissue mirroring in vivo early eye development - in vivo imaging of OV formation.
    https://doi.org/10.11588/data/AOSUS8

References

  13. Book
    1. Chauhan B
    2. Plageman T
    3. Lou M
    4. Lang R
    (2015) Chapter Eleven - Epithelial Morphogenesis: The Mouse Eye as a Model System
    In: Trainor P. A, editors. Current Topics in Developmental Biology. Academic Press. pp. 375–399.
    https://doi.org/10.1016/bs.ctdb.2014.11.011
    1. Loosli F
    2. Winkler S
    3. Burgtorf C
    4. Wurmbach E
    5. Ansorge W
    6. Henrich T
    7. Grabher C
    8. Arendt D
    9. Carl M
    10. Krone A
    11. Grzebisz E
    12. Wittbrodt J
    (2001)
    Medaka eyeless is the key factor linking retinal determination and eye growth
    Development 128:4035–4044.
  88. Book
    1. Tinevez J-Y
    2. Herbert S
    (2020) The NEMO Dots Assembly: Single-Particle Tracking and Analysis
    In: Miura K, Sladoje N, editors. Bioimage Data Analysis Workflows, Learning Materials in Biosciences. Springer. pp. 67–96.
    https://doi.org/10.1007/978-3-030-22386-1_4
  89. Software
    1. Tischer C
    (2019) ElastixWrapper: Fiji Plugin for 3D Image Registration with Elastix, Zenodo
    ElastixWrapper: Fiji Plugin for 3D Image Registration with Elastix, Zenodo.
    https://zenodo.org/record/2602549
    1. Turing AM
    (1952) The chemical basis of morphogenesis
    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 237:37–72.
    https://doi.org/10.1098/rstb.1952.0012

Decision letter

  1. Alfonso Martínez Arias
    Reviewing Editor; Universitat Pompeu Fabra, Spain
  2. Didier YR Stainier
    Senior Editor; Max Planck Institute for Heart and Lung Research, Germany
  3. Alfonso Martínez Arias
    Reviewer; Universitat Pompeu Fabra, Spain

Our editorial process produces two outputs: i) public reviews designed to be posted alongside the preprint for the benefit of readers; ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Acceptance summary:

In this manuscript, Zilova et al. show that primary embryonic cells derived from blastula-stage Medaka and Zebrafish embryos can self-organize into retinal organoids. This is a novel and original piece of work that reveals the capacity of fish primary embryonic pluripotent cells to behave like mammalian embryonic stem cells and organize optic cup organoids.

Decision letter after peer review:

Thank you for submitting your article "Fish primary embryonic stem cells self-assemble into retinal tissue mirroring in vivo early eye development" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Alfonso Martínez Arias as the Reviewing Editor and Reviewer #1. and the evaluation has been overseen by Didier Stainier as the Senior Editor.

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) The introduction is very sparse and should be expanded. At the moment it reads as it had been cut off. A more thorough introduction of the issues will help highlighting the novelty of their observations. In it they should refer to PMID 17540356 and, in particular Figure 2 and references therein about old newt experiments hinting at what they observe here. They should also avoid referring to their cells as 'stem cell', they are (as they say elsewhere) primary embryonic pluripotent cells.

2) It would be helpful to have a timeline of retinal development at the beginning of the paper (figure 1). In this way, it will be easier for the reader to follow the figures in the paper that go from the neuroepithelium to the optic-vesicle like structures in the organoids.

3) In Figure 1d, they should show markers of epithelial cells e.g Cadherins or ZO-1.

4) In Figure 2, It is important to show Rx3 expression before Matrigel embedding. It seems to be polarized but if this is the case, it is important to show the frequency of this event. It would also be good to have a video of the development of the organoids in Matrigel which, at the moment, is missing (or we have not found it). It is important to show the transition from 'polarised' Rx3 expression in the organoids go to the Rx3 expression in the neuroepithelium?

5) In Figure 4, the authors should specify that this is a Medaka experiment. Figure 4a is very small. It would be important to have a large view of the organoids and the existent variability in these.

6) In the smaller organoids (figure 4a, e, g), is there epithelium formation similar to the regular organoids? Because it does not seem to. So, is epithelium necessary to progress to retinal differentiation?

7) In Figure 4, the authors show that the optic vesicle organoids are organized as in vivo with cells expressing RPE markers. These cells are no longer present in Figure 5. What happens to them? There is no mention of this problem in the text. This should be addressed or a least discussed. The RPE absence may be a reason why the retina differentiates with an inverted organization.

8) The authors variably interpret their observations as the result of self-assembly or self-organization. At the moment, the data does not allow distinguishing whether the observed phenomena result from cells following largely cell-autonomous differentiation paths and come together through cell sorting, or whether dissociation and aggregation generates a condition that leads to (spatially restricted) retinal differentiation in cells that would not normally adopt this fate. I would say that the first scenario is consistent with self-assembly, while the second one is more self-organized in the sense that the new cell-cell interactions resulting from the aggregation result in emergent cellular behaviours. A first step to distinguish between these possibilities would be to quantitatively demonstrate that aggregation biases cell differentiation towards neural and retinal fates at the expense of other cell types, compared to the intact embryo. The examples shown in Figure 2 and 3d seem to indicate an overrepresentation of neural cells, but it would be good to see a quantitative comparison to the embryo.

9) The authors claim that their system is highly reproducible. Unfortunately, they do not give an indication of the success rate of aggregate formation in figure 1. Figure 4 shows the most complex patterns, but I realize that there is quite a bit of variability in between the aggregates – they are just as likely to have one or two Rx2-expressing areas (panel b). I also could not find information how many aggregates show the patterns in panels e and f, and from how many aggregates the data in panels g – i has been collected.

Reviewer #1 (Recommendations for the authors):

In this manuscript, Zilova et al. show that primary embryonic cells derived from blastula-stage Medaka and Zebrafish embryos can self-organize into retinal organoids. When aggregates of 1000-2000 primary embryonic cell are embedded in Matrigel addition, they form a neuroepithelium under the control of Rx3 which develops into a retinal organoid. The process mirrors some aspects of embryo development. Moreover, another interesting finding is that Rx3 expression is initiated in the absence of Matrigel at day 0, which indicates that the retinal fate occurs by default and is not dependent on extracellular matrix components. The authors compare the ability of cells from Mesaka and zebra fish and show that both are competent to form organoids, though each does it with the time scale of the embryo of origin. The authors show that by reducing the number of Medaka cells to aggregate (500-800 cells), Rx2 and Rx3 are expressed only in restricted regions of the small aggregates, presumably where they organize into discrete circular Rx2 and Rx3 positive neuroepithelial units that develop into structure resembling retinal epitjhelia with some diversity of retinal cell types including amacrine, ganglion, photoreceptor, bipolar and horizontal cells.

This is a novel and original piece of work that reveals the capacity of fish primary embryonic pluripotent cells to behave like mammalian embryonic stem cells and organize optic cup organoids. It is worth considering it for publication but first the authors should address a number of issues that could improve the presentation of their findings and also expand their observations beyond the description of the phenomenon.

Details are indicated below but there are two pieces of additional information that would unwrap the potential of their observations. The first one is the, these days almost obligatory, single cell analysis of their retinal organoids compared to the embryo, particularly of the retinal epithelium. Of course, were they able to deepen the analysis in Figure 5 this would not be necessary, but a more complete account of the cell types present in the organoids is necessary as is whether (or how much brain tissue there is). In addition, the observations of the different timings of the ex vivo development of Medaka and zebrafish begs the experiment of doing the experiment of generating chimeric organoids and see how the combinations of cells from different species affects their development. A more detailed characterization of the development of the organoids in comparison with the embryo and, importantly, a detailed reporting of the frequency of the events, will be appreciated.

The present manuscript also requires clarification of several points

The introduction is very sparse and should be expanded. At the moment it reads as it had been cut off. A more thorough introduction of the issues will help highlighting the novelty of their observations. In it they should refer to PMID 17540356 and, in particular Figure 2 and references therein about old newt experiments hinting at what they observe here. They should also avoid referring to their cells as 'stem cell', they are (as they say elsewhere) primary embryonic pluripotent cells.

It would be helpful to have a timeline of retinal development at the beginning of the paper (figure 1). In this way, it will be easier for the reader to follow the figures in the paper that go from the neuroepithelium to the optic-vesicle like structures in the organoids.

Other issues they should address:

At the beginning, they should not refer to the structures as 'organoids' at this stage as they are only epithelial cysts at this point.

In Figure 1, they should discuss a comparison of figures 1c, d with figure 4a, c. Are these hollow in the middle or do they have cells? Is the middle GFP expression autofluorescence?

In Figure 1d, they should show markers of epithelial cells e.g Cadherins or ZO-1.

In Figure 2, It is important to show Rx3 expression before Matrigel embedding. It seems to be polarized but if this is the case, it is important to show the frequency of this event. It would also be good to have a video of the development of the organoids in Matrigel which, at the moment, is missing (or we have not found it). It is important to show the transition from 'polarised' Rx3 expression in the organoids go to the Rx3 expression in the neuroepithelium?

Figure 2g – It would be easier to read and understand by depicting Rx3-/- for the mutants rather than Rx3::saGFP+/+.

On the zebrafish experiments, Figure 3, need to state in-text the exact stage of the zebrafish from which the organoids where then generated. In methodology is stated: ' For zebrafish cell aggregation, blastula-stage embryos were processed according to the same protocol and re-suspended in Leibowitz's L-15 media supplemented with 5% KSR and penicillin-streptomycin. Matrigel (Corning, 356238) was added 2 hpa'. But the authors should point the exact number of cells here. How many where they added? Have the authors tested aggregating different number of cells as with Medaka and examined Rx3 expression and morphogenesis? Important to have more comparisons with the medaka.

It would also be helpful to depict the protocol of the organoids generation from Zebrafish in a schematic as done with the Medaka. And indicate until which time point and developmental stage can they progress. Do we see similar retinal differentiation as with the medaka (figure 5)?

Compare Figure 3b with 3d: Looking at the organoids, one can see major differences here. The organoid shown at figure 3b is circular with an outer neuroepithelium (assuming also that Rx3 expression is circular as with medaka) whereas organoids in figure 3d are variable with Rx3 expression being polarised. Also, again is this here neuroepithelium? It would be important to stain for Cadherin.

In Figure 4, the authors should specify that this is a Medaka experiment. Figure 4a is very small. It would be important to have a large view of the organoids and the existent variability in these.

In Figure 4c it is stated that the 100% of the regions of regular organoids express Rx2. But in figure 1d we can see that the center might be hollow. So, again, is there autofluorescence in figure 4a (also see figure 2g) or not? Need to clarify.

Figure 4e: how many organoids show this pattern of gene expression? Is there variability when we have 2 poles of Rx3 expression? The issue of frequencies is very important in this field and the authors should report it.

In the smaller organoids (figure 4a, e, g), is there epithelium formation similar to the regular organoids? Because it does not seem to. So, is epithelium necessary to progress to retinal differentiation?

In Figure 5, compare figure 5b with 5c: One can see organoids with one pole of Atoh7. So, in figure 7c how often is this differentiation phenomenon?

Figure 5d: need to show an embryo to compare with the organoids' results.

Reviewer #2 (Recommendations for the authors):

In Figure 4, the authors show that the optic vesicle organoids are organized as in vivo with cells expressing RPE markers. These cells are no longer present in Figure 5. What happens to them? There is no mention of this problem in the text. This should be addressed or a least discussed. The RPE absence may be a reason why the retina differentiates with an inverted organization.

The discussion is generally informative but somehow fails to provide real advantages of using teleost organoids vs the fish per se or vs for example human organoids. Indeed, obtaining a fish organoid is faster that a human one, but more expensive and time consuming than using fish embryos. The author should perhaps provide more information in the introduction to explain what has push them to undertake this effort, besides the technical challenge. What can we learn from fish vs mammalian organoids?

Reviewer #3 (Recommendations for the authors):

1) It was not immediately obvious to me what motivated the authors to study retinal differentiation in the aggregates. (line 95 ff). This could be discussed more clearly.

2) It is unclear to me what the quantity "level of organization" in figure 2h means. Is there a more specific way to express this?

3) I was a bit confused with the positioning of the zebrafish experiments in the manuscript. It took me a while to realize that the experiments in Figure 4 were performed in the medaka system again. Perhaps the zebrafish data could be discussed elsewhere in the manuscript, or the authors could indicate more clearly in the text which system they used for the respective experiments.

4) The discussion dwells extensively on potential uses of the fish aggregate system as an alternative to organoids from mammalian cells. I am not convinced that this is a strong point, as potential advantages due to faster development are offset by evolutionary differences. Rather, I think that these systems are of interest in their own to investigate developmental mechanisms of cell differentiation and morphogenesis.

5) The reference to Fuhrmann et al., 2020 in line 389 seems to be wrong.

6) I find the last section of the discussion hard to grasp. If the authors are referring to previously described concepts, it might be helpful to add some references here.

https://doi.org/10.7554/eLife.66998.sa1

Author response

Essential revisions:

1) The introduction is very sparse and should be expanded. At the moment it reads as it had been cut off. A more thorough introduction of the issues will help highlighting the novelty of their observations. In it they should refer to PMID 17540356 and, in particular Figure 2 and references therein about old newt experiments hinting at what they observe here. They should also avoid referring to their cells as 'stem cell', they are (as they say elsewhere) primary embryonic pluripotent cells.

Indeed, the introduction was very concise. We happily responded to the suggestion of the referees and have extended the introduction and have taken care to incorporate and quote all the suggested work. Throughout the revised version of the manuscript we have taken care to consistently refer to the “starting material” as “primary embryonic pluripotent cells”.

2) It would be helpful to have a timeline of retinal development at the beginning of the paper (figure 1). In this way, it will be easier for the reader to follow the figures in the paper that go from the neuroepithelium to the optic-vesicle like structures in the organoids.

We followed the referee’s suggestion and have included a scheme of retina development in Figure 1 (panel a) for reference and introduce retina development in the extended introduction of the revised version of the manuscript.

3) In Figure 1d, they should show markers of epithelial cells e.g Cadherins or ZO-1.

We repeated the respective analyses and have now included N-cadherin staining on day 2 organoids presented in the new Figure 1 panel d, as well as in the linked supplements (Figure 1—figure supplement 2) and Figure 2 panel e.

4) In Figure 2, It is important to show Rx3 expression before Matrigel embedding. It seems to be polarized but if this is the case, it is important to show the frequency of this event. It would also be good to have a video of the development of the organoids in Matrigel which, at the moment, is missing (or we have not found it). It is important to show the transition from 'polarised' Rx3 expression in the organoids go to the Rx3 expression in the neuroepithelium?

To address the question concerning distribution of Rx3-expressing cells in relation to the addition of Matrigel, we have performed extended time-lapse imaging on all organoids (n=54) of a single batch, see video 2 (new panel c, Figure 2). Our analysis shows that there is no polarisation of Rx3 expression prior to the addition of Matrigel. Rx3 expression emerges as salt and pepper pattern and only subsequently continuous Rx3 expressing domains emerge as highlighted in the revised Figure 2. We have also included panels addressing Rx3 expression before and after Matrigel addition in panel e, Figure 2 showing the distribution of Rx3-expressing cells in the context of forming neuroepithelium.

Additionally, we included the data showing the impact of Matrigel on the formation of the neuroepithelium (Figure 2—figure supplement 2). Furthermore, Video 7 addresses the behaviour of Rx3-expressing cells in the presence or absence of Matrigel in a time frame of more than 28 hours. Videos 5 and 6 show behaviour of Rx3-expressing cells after Matrigel addition up to optic vesicle formation.

The referees had referred to an apparent polarisation of one organoid presented in figure 2 panel c of the initially submitted manuscript. This panel was misleading and as stated above we have in fact not observed any polarized Rx3 expression. On the contrary, Rx3 expression emerged in individual, isolated cells in a salt and pepper fashion. We have taken care to avoid the impression of polarized Rx3 expression (due to maximum projection and high gain) in the revised version, both in the figure as well as in the video.

With the addition of Matrigel and following epithelialisation, a continuous Rx3 expression domain is forming.

We have refined our statements following the feedback of the referees in the revised version of the manuscript.

5) In Figure 4, the authors should specify that this is a Medaka experiment. Figure 4a is very small. It would be important to have a large view of the organoids and the existent variability in these.

We have revised Figure 4 to include large overview images for panel a and Figure 4—figure supplement 1.

Both figures are excellent examples of the morphological variability of organoids generated by the aggregation of less than 1,000 cells.

For a better representation of the variable number of optic vesicles forming, apparent as evaginated Rx3-positive retinal regions, we have revised Figure 4 and present an improved panel c presenting the absolute number of organoids with 1, 2, 3 or 4 optic vesicles/retinal regions.

This allows to instantly appreciate the impact of the starting cell number on the number of forming optic vesicles/retinal regions.

Details on the morphological variability of 72 organoids derived from less than 1,000 cells are now presented in Figure 4—figure supplement 1.

6) In the smaller organoids (figure 4a, e, g), is there epithelium formation similar to the regular organoids? Because it does not seem to. So, is epithelium necessary to progress to retinal differentiation?

The formation of an epithelium does not depend on the number of starting cells. We have addressed this question, in Figure 6—figure supplement 1. Here we show that the onset of retinal differentiation in organoids does not depend on the starting size of the initial aggregate. There is no difference in the onset of retinal differentiation (analysed by the expression of Atoh7::EGFP) and layering of retinal cells. Although organoids generated by aggregation of > 1,000 and < 1,000 cells show different morphologies and distribution of retinal domains/optic vesicles, the forming retinal neuroepithelium follows the process or retinal differentiation irrespective of these morphological differences. Epithelium formation (in both, > 1,000 and < 1,000 organoids) is a prerequisite for the survival of organoids as is now presented in Figure 2—figure supplement 2.

7) In Figure 4, the authors show that the optic vesicle organoids are organized as in vivo with cells expressing RPE markers. These cells are no longer present in Figure 5. What happens to them? There is no mention of this problem in the text. This should be addressed or a least discussed. The RPE absence may be a reason why the retina differentiates with an inverted organization.

The referees relate to an important point: As in human and mouse retinal organoids, also in fish retinal organoids the differentiation of the neuroretina is favoured over the formation of RPE. In established systems, this can be induced by supplementing the culture media with Wnt/β-catenin pathway-inducing agents (Eiraku et al., 2011; Nakano et al., 2012; Kuwahara et al., 2015). We demonstrate that this approach also induces RPE in fish-derived organoid. In the revised version of the manuscript we have included our results on RPE differentiation by treatment with the canonical Wnt/β-catenin pathway agonist CHIR99021, see Figure 6—figure supplement 2 and discussed in the manuscript on page 21.

In the absence of a shielding RPE, ECM components provided by the addition of Matrigel polarizes the forming epithelium in an “inverted” fashion.

8) The authors variably interpret their observations as the result of self-assembly or self-organization. At the moment, the data does not allow distinguishing whether the observed phenomena result from cells following largely cell-autonomous differentiation paths and come together through cell sorting, or whether dissociation and aggregation generates a condition that leads to (spatially restricted) retinal differentiation in cells that would not normally adopt this fate. I would say that the first scenario is consistent with self-assembly, while the second one is more self-organized in the sense that the new cell-cell interactions resulting from the aggregation result in emergent cellular behaviours. A first step to distinguish between these possibilities would be to quantitatively demonstrate that aggregation biases cell differentiation towards neural and retinal fates at the expense of other cell types, compared to the intact embryo. The examples shown in Figure 2 and 3d seem to indicate an overrepresentation of neural cells, but it would be good to see a quantitative comparison to the embryo.

Here the referees open a debate that is controversial in the field and we are happy to contribute. We followed the referees’ suggestions on how to address the difference between self-assembly and self-organization and are providing two additional experiments tackling this point.

To address whether the cells would autonomously follow the retinal differentiation path we dissociated blastulae and cultured individual dissociated blastula-stage cells (Rx3::H2B-GFP reporter line) in suspension or as a single cells.

Our results added as Figure 2—figure supplement 1 show that, consistent with a default neuronal state, individual blastula-derived cells stochastically acquire retinal fate autonomously, in the absence of cell-cell contacts established by cell aggregation. The clones derived from individual cells eventually show a distribution of retinal (Rx3+) and non-retinal cells.

In the absence of a fully conclusive answer for either self-organization or self-assembly (or likely a combination of both) we have taken care not to confuse the reader by an imprecision in terminology in the crucial point. In light of the historical debate about the differences between self-assembly and self-organisation “Self-assembly, Self-Organization: A philosophical perspective on a major challenge of nanotechnology” by Bensaude-Vincent https://halshs.archives-ouvertes.fr/halshs-00350831/document, we have revised our manuscript accordingly.

9) The authors claim that their system is highly reproducible. Unfortunately, they do not give an indication of the success rate of aggregate formation in figure 1. Figure 4 shows the most complex patterns, but I realize that there is quite a bit of variability in between the aggregates – they are just as likely to have one or two Rx2-expressing areas (panel b). I also could not find information how many aggregates show the patterns in panels e and f, and from how many aggregates the data in panels g – i has been collected.

The reviewer touches a crucial point and actually one to the strengths of our system. Overall the derivation of retinal organoids from primary embryonic pluripotent cells is highly reproducible (in 100% of all experiments performed) with routinely almost 100 % efficiency per experiment (video 1 and video 2). We clearly state and present those key aspects in the revised version of the manuscript.

With our protocol on derivation of organoids, all combined primary embryonic pluripotent cells proceed through aggregation (84 out of 84 aggregates in video 1 and 2). In an independent experiment we show the highly synchronous onset of Rx3 expression (54 out of 54 aggregates in video 2). There is a variability in size (panel d Figure 4) and number of optic vesicles (panel c Figure 4 and Figure 4—figure supplement 1), all of which will eventually initiate retinal differentiation (panel b Figure 6—figure supplement 1). To allow the reader to get an immediate impression of the efficiency of the procedure, we have now added the number of organoids used in every experiment in figure captions.

https://doi.org/10.7554/eLife.66998.sa2

Article and author information

Author details

  1. Lucie Zilova

    Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
    Contribution
    Conceptualization, Data curation, Formal analysis, Supervision, Validation, Investigation, Visualization, Methodology, Writing - original draft
    Contributed equally with
    Venera Weinhardt
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6404-9119

  2. Venera Weinhardt

    Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Supervision, Validation, Investigation, Visualization, Methodology, Writing - original draft
    Contributed equally with
    Lucie Zilova
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9774-3833

  3. Tinatini Tavhelidse

    Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
    Contribution
    Resources, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6103-9019

  4. Christina Schlagheck

    1. Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
    2. Heidelberg International Biosciences Graduate School HBIGS and HeiKa Graduate School on “Functional Materials”, Heidelberg, Germany
    Contribution
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1311-5945

  5. Thomas Thumberger

    Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
    Contribution
    Resources, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8485-457X

  6. Joachim Wittbrodt

    Centre for Organismal Studies Heidelberg, Heidelberg University, Heidelberg, Germany
    Contribution
    Conceptualization, Resources, Supervision, Funding acquisition, Writing - original draft, Project administration
    For correspondence
    jochen.wittbrodt@cos.uni-heidelberg.de
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8550-7377

Funding

H2020 European Research Council (810172)

  • Joachim Wittbrodt

Deutsche Forschungsgemeinschaft (3DMM2O)

  • Joachim Wittbrodt

Deutsche Forschungsgemeinschaft (WE 6221/2-1)

  • Venera Weinhardt

Carl Zeiss Foundation (HeiKa Graduate School)

  • Christina Schlagheck

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The authors thank L Centanin, E Tsingos, N Sokolova, and Z Kozmik for valuable discussions and comments on the manuscript. We are grateful to I Thomas, A Sarvari, M Pandya, and N Priya for the help at the initial stage of the project. We thank Darius Balciunas for providing the GBT-RP2 plasmid containing the OPT sequence. VW was supported by the German Research Foundation research fellowship WE 6221/2–1. This work was supported by grants of the Excellence Cluster ‘3D Matter Made to Order’ (3DMM2O) funded through the German Excellence Strategy via Deutsche Forschungsgemeinschaft (DFG), by the Carl Zeiss Foundation and by the ERC Synergy Grant IndiGene (Number 810172) to JW.

Senior Editor

  1. Didier YR Stainier, Max Planck Institute for Heart and Lung Research, Germany

Reviewing Editor

  1. Alfonso Martínez Arias, Universitat Pompeu Fabra, Spain

Reviewer

  1. Alfonso Martínez Arias, Universitat Pompeu Fabra, Spain

Publication history

  1. Received: January 28, 2021
  2. Accepted: June 24, 2021
  3. Version of Record published: July 12, 2021 (version 1)

Copyright

© 2021, Zilova et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 3,298
    Page views
  • 328
    Downloads
  • 7
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Lucie Zilova
  2. Venera Weinhardt
  3. Tinatini Tavhelidse
  4. Christina Schlagheck
  5. Thomas Thumberger
  6. Joachim Wittbrodt
(2021)
Fish primary embryonic pluripotent cells assemble into retinal tissue mirroring in vivo early eye development
eLife 10:e66998.
https://doi.org/10.7554/eLife.66998

Categories and tags

Further reading

    1. Stem Cells and Regenerative Medicine
    2. Developmental Biology

    Integrated transcriptome and proteome analysis reveals posttranscriptional regulation of ribosomal genes in human brain organoids

    Jaydeep Sidhaye, Philipp Trepte ... Jürgen A Knoblich
    Research Article

    During development of the human cerebral cortex, multipotent neural progenitors generate excitatory neurons and glial cells. Investigations of the transcriptome and epigenome have revealed important gene regulatory networks underlying this crucial developmental event. However, the posttranscriptional control of gene expression and protein abundance during human corticogenesis remains poorly understood. We addressed this issue by using human telencephalic brain organoids grown using a dual reporter cell line to isolate neural progenitors and neurons and performed cell class and developmental stage-specific transcriptome and proteome analysis. Integrating the two datasets revealed modules of gene expression during human corticogenesis. Investigation of one such module uncovered mTOR-mediated regulation of translation of the 5’TOP element-enriched translation machinery in early progenitor cells. We show that in early progenitors partial inhibition of the translation of ribosomal genes prevents precocious translation of differentiation markers. Overall, our multiomics approach proposes novel posttranscriptional regulatory mechanisms crucial for the fidelity of cortical development.

    1. Developmental Biology
    2. Neuroscience

    Origin of wiring specificity in an olfactory map revealed by neuron type-specific, time-lapse imaging of dendrite targeting

    Kenneth Kin Lam Wong, Tongchao Li ... Liqun Luo
    Research Article

    How does wiring specificity of neural maps emerge during development? Formation of the adult Drosophila olfactory glomerular map begins with patterning of projection neuron (PN) dendrites at the early pupal stage. To better understand the origin of wiring specificity of this map, we created genetic tools to systematically characterize dendrite patterning across development at PN type-specific resolution. We find that PNs use lineage and birth order combinatorially to build the initial dendritic map. Specifically, birth order directs dendrite targeting in rotating and binary manners for PNs of the anterodorsal and lateral lineages, respectively. Two-photon- and adaptive optical lattice light-sheet microscope-based time-lapse imaging reveals that PN dendrites initiate active targeting with direction-dependent branch stabilization on the timescale of seconds. Moreover, PNs that are used in both the larval and adult olfactory circuits prune their larval-specific dendrites and re-extend new dendrites simultaneously to facilitate timely olfactory map organization. Our work highlights the power and necessity of type-specific neuronal access and time-lapse imaging in identifying wiring mechanisms that underlie complex patterns of functional neural maps.

    1. Developmental Biology

    Histone deacetylase 1 maintains lineage integrity through histone acetylome refinement during early embryogenesis

    Jeff Jiajing Zhou, Jin Sun Cho ... Ken WY Cho
    Research Article

    Histone acetylation is a pivotal epigenetic modification that controls chromatin structure and regulates gene expression. It plays an essential role in modulating zygotic transcription and cell lineage specification of developing embryos. While the outcomes of many inductive signals have been described to require enzymatic activities of histone acetyltransferases and deacetylases (HDACs), the mechanisms by which HDACs confine the utilization of the zygotic genome remain to be elucidated. Here, we show that histone deacetylase 1 (Hdac1) progressively binds to the zygotic genome from mid blastula and onward. The recruitment of Hdac1 to the genome at blastula is instructed maternally. Cis-regulatory modules (CRMs) bound by Hdac1 possess epigenetic signatures underlying distinct functions. We highlight a dual function model of Hdac1 where Hdac1 not only represses gene expression by sustaining a histone hypoacetylation state on inactive chromatin, but also maintains gene expression through participating in dynamic histone acetylation-deacetylation cycles on active chromatin. As a result, Hdac1 maintains differential histone acetylation states of bound CRMs between different germ layers and reinforces the transcriptional program underlying cell lineage identities, both in time and space. Taken together, our study reveals a comprehensive role for Hdac1 during early vertebrate embryogenesis.