A single-cell guide to retinal development: Cell fate decisions of multipotent retinal progenitors in scRNA-seq

Unraveling the developmental heterogeneity within the human retina to reconstruct the continuity of retinal ganglion cell maturation and stage-specific intrinsic and extrinsic factors

Tissue development is a complex spatiotemporal process with multiple interdependent components. Anatomical, histological, sequencing, and evolutional strategies can be used to profile and explain tissue development from different perspectives. The introduction of scRNAseq methods and the computational tools allows to deconvolute developmental heterogeneity and draw a decomposed uniform map. In this manuscript, we decomposed the development of a human retina with a focus on the retinal ganglion cells (RGC). To increase the temporal resolution of retinal cell classes maturation state we assumed the working hypothesis that that maturation of retinal ganglion cells is a continuous, non-discrete process. We have assembled the scRNAseq atlas of human fetal retina from fetal week 8 to week 27 and applied the computational methods to unravel maturation heterogeneity into a uniform maturation track. We align RGC transcriptomes in pseudotime to map RGC developmental fate trajectories against the broader timeline of retinal development. Through this analysis, we identified the continuous maturation track of RGC and described the cell-intrinsic (DEGs, maturation gene profiles, regulons, transcriptional motifs) and -extrinsic profiles (neurotrophic receptors across maturation, cell-cell interactions) of different RGC maturation states. We described the genes involved in the retina and RGC maturation, including de novo RGC maturation drivers. We demonstrate the application of the human fetal retina atlas as a reference tool, allowing automated annotation and universal embedding of scRNAseq data. Altogether, our findings deepen the current knowledge of the retina and RGC maturation by bringing in the maturation dimension for the cell class vs. state analysis. We show how the pseudotime application contributes to developmental-oriented analyses, allowing to order the cells by their maturation state. This approach not only improves the downstream computational analysis but also provides a true maturation track transcriptomics profile.

Identification and characterization of early human photoreceptor states and cell-state-specific retinoblastoma-related features

Human cone photoreceptors differ from rods and serve as the retinoblastoma cell-of-origin, yet the developmental basis for their distinct behaviors is poorly understood. Here, we used deep full-length single-cell RNA-sequencing to distinguish post-mitotic cone and rod developmental states and identify cone-specific features that contribute to retinoblastomagenesis. The analyses revealed early post-mitotic cone- and rod-directed populations characterized by higher THRB or NRL regulon activities, an immature photoreceptor precursor population with concurrent cone and rod gene and regulon expression, and distinct early and late cone and rod maturation states distinguished by maturation-associated declines in RAX regulon activity. Unexpectedly, both L/M cone and rod precursors co-expressed NRL and THRB RNAs, yet they differentially expressed functionally antagonistic NRL and THRB isoforms and prematurely terminated THRB transcripts. Early L/M cone precursors exhibited successive expression of several lncRNAs along with MYCN, which composed the seventh most L/M-cone-specific regulon, and SYK, which contributed to the early cone precursors' proliferative response to RB1 loss. These findings reveal previously unrecognized photoreceptor precursor states and a role for early cone-precursor-intrinsic SYK expression in retinoblastoma initiation.

Single cell dual-omic atlas of the human developing retina

The development of the retina is under tight temporal and spatial control. To gain insights into the molecular basis of this process, we generate a single-nuclei dual-omic atlas of the human developing retina with approximately 220,000 nuclei from 14 human embryos and fetuses aged between 8 and 23-weeks post-conception with matched macular and peripheral tissues. This atlas captures all major cell classes in the retina, along with a large proportion of progenitors and cell-type-specific precursors. Cell trajectory analysis reveals a transition from continuous progression in early progenitors to a hierarchical development during the later stages of cell type specification. Both known and unrecorded candidate transcription factors, along with gene regulatory networks that drive the transitions of various cell fates, are identified. Comparisons between the macular and peripheral retinae indicate a largely consistent yet distinct developmental pattern. This atlas offers unparalleled resolution into the transcriptional and chromatin accessibility landscapes during development, providing an invaluable resource for deeper insights into retinal development and associated diseases.

Systematic mapping of TF-mediated cell fate changes by a pooled induction coupled with scRNA-seq and multi-omics approaches

Transcriptional regulation controls cellular functions through interactions between transcription factors (TFs) and their chromosomal targets. However, understanding the fate conversion potential of multiple TFs in an inducible manner remains limited. Here, we introduce iTF-seq as a method for identifying individual TFs that can alter cell fate toward specific lineages at a single-cell level. iTF-seq enables time course monitoring of transcriptome changes, and with biotinylated individual TFs, it provides a multi-omics approach to understanding the mechanisms behind TF-mediated cell fate changes. Our iTF-seq study in mouse embryonic stem cells identified multiple TFs that trigger rapid transcriptome changes indicative of differentiation within a day of induction. Moreover, cells expressing these potent TFs often show a slower cell cycle and increased cell death. Further analysis using bioChIP-seq revealed that GCM1 and OTX2 act as pioneer factors and activators by increasing gene accessibility and activating the expression of lineage specification genes during cell fate conversion. iTF-seq has utility in both mapping cell fate conversion and understanding cell fate conversion mechanisms.

A framework to identify functional interactors that contribute to disrupted early retinal development in Vsx2 ocular retardation J mice

BACKGROUND

A goal of developmental genetics is to identify functional interactions that underlie phenotypes caused by mutations. We sought to identify functional interactors of Vsx2, which when mutated, disrupts early retinal development. We utilized the Vsx2 loss-of-function mouse, ocular retardation J (orJ), to assess interactions based on principles of positive and negative epistasis as applied to bulk transcriptome data. This was first tested in vivo with Mitf, a target of Vsx2 repression, and then to cultures of orJ retina treated with inhibitors of Retinoid-X Receptors (RXR) to target Rxrg, an up-regulated gene in the orJ retina, and gamma-Secretase, an enzyme required for Notch signaling, a key mediator of retinal proliferation and neurogenesis.

RESULTS

Whereas Mitf exhibited robust positive epistasis with Vsx2, it only partially accounts for the orJ phenotype, suggesting other functional interactors. RXR inhibition yielded minimal evidence for epistasis between Vsx2 and Rxrg. In contrast, gamma-Secretase inhibition caused hundreds of Vsx2-dependent genes associated with proliferation to deviate further from wild-type, providing evidence for convergent negative epistasis with Vsx2 in regulating tissue growth.

CONCLUSIONS

Combining in vivo and ex vivo testing with transcriptome analysis revealed quantitative and qualitative characteristics of functional interaction between Vsx2, Mitf, RXR, and gamma-Secretase activities.

Notch pathway mutants do not equivalently perturb mouse embryonic retinal development

In the vertebrate eye, Notch ligands, receptors, and ternary complex components determine the destiny of retinal progenitor cells in part by regulating Hes effector gene activity. There are multiple paralogues for nearly every node in this pathway, which results in numerous instances of redundancy and compensation during development. To dissect such complexity at the earliest stages of eye development, we used seven germline or conditional mutant mice and two spatiotemporally distinct Cre drivers. We perturbed the Notch ternary complex and multiple Hes genes to understand if Notch regulates optic stalk/nerve head development; and to test intracellular pathway components for their Notch-dependent versus -independent roles during retinal ganglion cell and cone photoreceptor competence and fate acquisition. We confirmed that disrupting Notch signaling universally blocks progenitor cell growth, but delineated specific pathway components that can act independently, such as sustained Hes1 expression in the optic stalk/nerve head. In retinal progenitor cells, we found that among the genes tested, they do not uniformly suppress retinal ganglion cell or cone differentiation; which is not due differences in developmental timing. We discovered that shifts in the earliest cell fates correlate with expression changes for the early photoreceptor factor Otx2, but not with Atoh7, a factor required for retinal ganglion cell formation. During photoreceptor genesis we also better defined multiple and simultaneous activities for Rbpj and Hes1 and identify redundant activities that occur downstream of Notch. Given its unique roles at the retina-optic stalk boundary and cone photoreceptor genesis, our data suggest Hes1 as a hub where Notch-dependent and -independent inputs converge.

Cell-type specification in the retina: Recent discoveries from transcriptomic approaches

Understanding the formation of the complex nervous system hinges on decoding the mechanism that specifies a vast array of neuronal types, each endowed with a unique morphology, physiology, and connectivity. As a pivotal step towards addressing this problem, seminal work has been devoted to characterizing distinct neuronal types. In recent years, high-throughput, single-cell transcriptomic methods have enabled a rapid inventory of cell types in various regions of the nervous system, with the retina exhibiting complete molecular characterization across many vertebrate species. This invaluable resource has furnished a fresh perspective for investigating the molecular principles of cell-type specification, thereby advancing our understanding of retinal development. Accordingly, this review focuses on the most recent transcriptomic characterizations of retinal cells, with a particular focus on amacrine cells and retinal ganglion cells. These investigations have unearthed new insights into their cell-type specification.

Single-Cell RNA Sequencing: Opportunities and Challenges for Studies on Corneal Biology in Health and Disease

The structure and major cell types of the multi-layer human cornea have been extensively studied. However, various cell states in specific cell types and key genes that define the cell states are not fully understood, hindering our comprehension of corneal homeostasis, related diseases, and therapeutic discovery. Single-cell RNA sequencing is a revolutionary and powerful tool for identifying cell states within tissues such as the cornea. This review provides an overview of current single-cell RNA sequencing studies on the human cornea, highlighting similarities and differences between them, and summarizing the key genes that define corneal cell states reported in these studies. In addition, this review discusses the opportunities and challenges of using single-cell RNA sequencing to study corneal biology in health and disease.

Time to see: How temporal identity factors specify the developing mammalian retina

Understanding how retinal progenitor cells (RPCs) give rise to the variety of neural cell types of the retina has been a question of major interest over the last few decades. While environmental cues and transcription factor networks have been shown to control specific cell fate decisions, how RPCs alter fate output over time to control proper histogenesis remains poorly understood. In recent years, the identification of "temporal identity factors (TIFs)", which control RPC competence states to ensure that the right cell types are produced at the right time, has contributed to increasing our understanding of temporal patterning in the retina. Here, we review the different TIFs identified to date in the mammalian retina and discuss the underlying mechanisms by which they are thought to operate. We conclude by speculating on how identification of temporal patterning mechanisms might support the development of new therapeutic approaches against visual impairments.

Transcriptomics of CD29+/CD44+ cells isolated from hPSC retinal organoids reveals a single cell population with retinal progenitor and Müller glia characteristics

Müller glia play very important and diverse roles in retinal homeostasis and disease. Although much is known of the physiological and morphological properties of mammalian Müller glia, there is still the need to further understand the profile of these cells during human retinal development. Using human embryonic stem cell-derived retinal organoids, we investigated the transcriptomic profiles of CD29+/CD44+ cells isolated from early and late stages of organoid development. Data showed that these cells express classic markers of retinal progenitors and Müller glia, including NFIX, RAX, PAX6, VSX2, HES1, WNT2B, SOX, NR2F1/2, ASCL1 and VIM, as early as days 10-20 after initiation of retinal differentiation. Expression of genes upregulated in CD29+/CD44+ cells isolated at later stages of organoid development (days 50-90), including NEUROG1, VSX2 and ASCL1 were gradually increased as retinal organoid maturation progressed. Based on the current observations that CD24+/CD44+ cells share the characteristics of early and late-stage retinal progenitors as well as of mature Müller glia, we propose that these cells constitute a single cell population that upon exposure to developmental cues regulates its gene expression to adapt to functions exerted by Müller glia in the postnatal and mature retina.

Identity, lineage and fates of a temporally distinct progenitor population in the embryonic olfactory epithelium

We defined a temporally and transcriptionally divergent precursor cohort in the medial olfactory epithelium (OE) shortly after it differentiates as a distinct tissue at mid-gestation in the mouse. This temporally distinct population of Ascl1+ cells in the dorsomedial OE is segregated from Meis1+/Pax7+ progenitors in the lateral OE, and does not appear to be generated by Pax7+ lateral OE precursors. The medial Ascl1+ precursors do not yield a substantial number of early-generated ORNs. Instead, they first generate additional proliferative precursors as well as a distinct population of frontonasal mesenchymal cells associated with the migratory mass that surrounds the nascent olfactory nerve. Parallel to these in vivo distinctions, isolated medial versus lateral OE precursors in vitro retain distinct proliferative capacities and modes of division that reflect their in vivo identities. At later fetal stages, these early dorsomedial Ascl1+ precursors cells generate spatially restricted subsets of ORNs as well as other non-neuronal cell classes. Accordingly, the initial compliment of ORNs and other OE cell types is derived from at least two distinct early precursor populations: lateral Meis1/Pax7+ precursors that generate primarily early ORNs, and a temporally, spatially, and transcriptionally distinct subset of medial Ascl1+ precursors that initially generate additional OE progenitors and apparent migratory mass cells before yielding a subset of ORNs and likely supporting cell classes.

Deletion of histone demethylase Lsd1 (Kdm1a) during retinal development leads to defects in retinal function and structure

Purpose

The purpose of this study was to investigate the role of Lysine specific demethylase 1 (Lsd1) in murine retinal development. LSD1 is a histone demethylase that can demethylate mono- and di-methyl groups on H3K4 and H3K9. Using Chx10-Cre and Rho-iCre75 driver lines, we generated novel transgenic mouse lines to delete Lsd1 in most retinal progenitor cells or specifically in rod photoreceptors. We hypothesize that Lsd1 deletion will cause global morphological and functional defects due to its importance in neuronal development.

Methods

We tested the retinal function of young adult mice by electroretinogram (ERG) and assessed retinal morphology by in vivo imaging by fundus photography and SD-OCT. Afterward, eyes were enucleated, fixed, and sectioned for subsequent hematoxylin and eosin (H&E) or immunofluorescence staining. Other eyes were plastic fixed and sectioned for electron microscopy.

Results

In adult Chx10-Cre Lsd1fl/fl mice, we observed a marked reduction in a-, b-, and c-wave amplitudes in scotopic conditions compared to age-matched control mice. Photopic and flicker ERG waveforms were even more sharply reduced. Modest reductions in total retinal thickness and outer nuclear layer (ONL) thickness were observed in SD-OCT and H&E images. Lastly, electron microscopy revealed significantly shorter inner and outer segments and immunofluorescence showed modest reductions in specific cell type populations. We did not observe any obvious functional or morphological defects in the adult Rho-iCre75 Lsd1fl/fl animals.

Conclusion

Lsd1 is necessary for neuronal development in the retina. Adult Chx10-Cre Lsd1fl/fl mice show impaired retinal function and morphology. These effects were fully manifested in young adults (P30), suggesting that Lsd1 affects early retinal development in mice.

Not all Notch pathway mutations are equal in the embryonic mouse retina

In the vertebrate retina, combinations of Notch ligands, receptors, and ternary complex components determine the destiny of retinal progenitor cells by regulating Hes effector gene activity. Owing to reiterated Notch signaling in numerous tissues throughout development, there are multiple vertebrate paralogues for nearly every node in this pathway. These Notch signaling components can act redundantly or in a compensatory fashion during development. To dissect the complexity of this pathway during retinal development, we used seven germline or conditional mutant mice and two spatiotemporally distinct Cre drivers. We perturbed the Notch ternary complex and multiple Hes genes with two overt goals in mind. First, we wished to determine if Notch signaling is required in the optic stalk/nerve head for Hes1 sustained expression and activity. Second, we aimed to test if Hes1, 3 and 5 genes are functionally redundant during early retinal histogenesis. With our allelic series, we found that disrupting Notch signaling consistently blocked mitotic growth and overproduced ganglion cells, but we also identified two significant branchpoints for this pathway. In the optic stalk/nerve head, sustained Hes1 is regulated independent of Notch signaling, whereas during photoreceptor genesis both Notch-dependent and -independent roles for Rbpj and Hes1 impact photoreceptor genesis in opposing manners.

Comprehensive characterization of fetal and mature retinal cell identity to assess the fidelity of retinal organoids

Characterizing cell identity in complex tissues such as the human retina is essential for studying its development and disease. While retinal organoids derived from pluripotent stem cells have been widely used to model development and disease of the human retina, there is a lack of studies that have systematically evaluated the molecular and cellular fidelity of the organoids derived from various culture protocols in recapitulating their in vivo counterpart. To this end, we performed an extensive meta-atlas characterization of cellular identities of the human eye, covering a wide range of developmental stages. The resulting map uncovered previously unknown biomarkers of major retinal cell types and those associated with cell-type-specific maturation. Using our retinal-cell-identity map from the fetal and adult tissues, we systematically assessed the fidelity of the retinal organoids in mimicking the human eye, enabling us to comprehensively benchmark the current protocols for retinal organoid generation.

Emergence of neuron types

Neuron types are the building blocks of the nervous system, and therefore, of functional circuits. Understanding the origin of neuronal diversity has always been an essential question in neuroscience and developmental biology. While knowledge on the molecular control of their diversification has largely increased during the last decades, it is now possible to reveal the dynamic mechanisms and the actual stepwise molecular changes occurring at single-cell level with the advent of single-cell omics technologies and analysis with high temporal resolution. Here, we focus on recent advances in the field and in technical and analytical tools that enable detailed insights into the emergence of neuron types in the central and peripheral nervous systems.

Cell Fate of Retinal Progenitor Cells: In Ovo UbC-StarTrack Analysis

Clonal cell analysis outlines the ontogenic potential of single progenitor cells, allowing the elucidation of the neural heterogeneity among different cell types and their lineages. In this work, we analyze the potency of retinal stem/progenitor cells through development using the chick embryo as a model. We implemented in ovo the clonal genetic tracing strategy UbC-StarTrack for tracking retinal cell lineages derived from individual progenitors of the ciliary margin at E3.5 (HH21-22). The clonal assignment of the derived-cell progeny was performed in the neural retina at E11.5-12 (HH38) through the identification of sibling cells as cells expressing the same combination of fluorophores. Moreover, cell types were assessed based on their cellular morphology and laminar location. Ciliary margin derived-cell progenies are organized in columnar associations distributed along the peripheral retina with a limited tangential dispersion. The analysis revealed that, at the early stages of development, this region harbors multipotent and committed progenitor cells.

A complete temporal transcription factor series in the fly visual system

The brain consists of thousands of neuronal types that are generated by stem cells producing different neuronal types as they age. In Drosophila, this temporal patterning is driven by the successive expression of temporal transcription factors (tTFs)1-6. Here we used single-cell mRNA sequencing to identify the complete series of tTFs that specify most Drosophila optic lobe neurons. We verify that tTFs regulate the progression of the series by activating the next tTF(s) and repressing the previous one(s), and also identify more complex mechanisms of regulation. Moreover, we establish the temporal window of origin and birth order of each neuronal type in the medulla and provide evidence that these tTFs are sufficient to explain the generation of all of the neuronal diversity in this brain region. Finally, we describe the first steps of neuronal differentiation and show that these steps are conserved in humans. We find that terminal differentiation genes, such as neurotransmitter-related genes, are present as transcripts, but not as proteins, in immature larval neurons. This comprehensive analysis of a temporal series of tTFs in the optic lobe offers mechanistic insights into how tTF series are regulated, and how they can lead to the generation of a complete set of neurons.

Single-cell RNA-sequencing analysis of the ciliary epithelium and contiguous tissues in the mouse eye

The ciliary epithelium plays a central role in ocular homeostasis but cells of the pigmented and non-pigmented layers are difficult to isolate physically and study. Here we used single-cell RNA-sequencing (scRNA-seq) to analyze the transcriptional signatures of cells harvested from the ciliary body and contiguous tissues. Microdissected tissue was dissociated by collagenase digestion and the transcriptomes of individual cells were obtained using a droplet-based scRNA-seq approach. In situ hybridization was used to verify the expression patterns of selected differentially-expressed genes. High quality transcriptomes were obtained from 10,024 cells and unsupervised clustering distinguished 22 cell types. Although efforts were made to specifically isolate the ciliary body, approximately half of the sequenced cells were derived from the adjacent retina. Cluster identities were assigned using expression of canonical markers or cluster-specific genes. The transcriptional signature of cells in the PCE and NPCE were distinct from each other and from cells in contiguous tissues. PCE cell transcriptomes were characterized by genes involved in melanin synthesis and transport proteins such as Slc4a4. Among the most differentially expressed genes in NPCE cells were those encoding members of the Zic family of transcription factors (Zic1, 2, 4), collagen XVIII (Col18a1), and corticotrophin-releasing hormone-binding protein (Crhbp). The ocular melanocyte population was distinguished by expression of the gap junction genes Gjb2 and Gjb6. Two fibroblast signatures were detected in the ciliary body preparation and shown by in situ hybridization to correspond to uveal and scleral populations. This cell atlas for the ciliary body and contiguous layers represents a useful resource that may facilitate studies into the development of the ciliary epithelium, the production of the aqueous and vitreous humors, and the synthesis of the ciliary zonule.