Analysis of gene expression in wild-type and Notch1 mutant retinal cells by single cell profiling

Application of Human Stem Cell Derived Retinal Organoids in the Exploration of the Mechanisms of Early Retinal Development

The intricate neural circuit of retina extracts salient features of the natural world and forms bioelectric impulse as the origin of vision. The early development of retina is a highly complex and coordinated process in morphogenesis and neurogenesis. Increasing evidence indicates that stem cells derived human retinal organoids (hROs) in vitro faithfully recapitulates the embryonic developmental process of human retina no matter in the transcriptome, cellular biology and histomorphology. The emergence of hROs greatly deepens on the understanding of early development of human retina. Here, we reviewed the events of early retinal development both in animal embryos and hROs studies, which mainly comprises the formation of optic vesicle and optic cup shape, differentiation of retinal ganglion cells (RGCs), photoreceptor cells (PRs) and its supportive retinal pigment epithelium cells (RPE). We also discussed the classic and frontier molecular pathways up to date to decipher the underlying mechanisms of early development of human retina and hROs. Finally, we summarized the application prospect, challenges and cutting-edge techniques of hROs for uncovering the principles and mechanisms of retinal development and related developmental disorder. hROs is a priori selection for studying human retinal development and function and may be a fundamental tool for unlocking the unknown insight into retinal development and disease.

Cellular and Molecular Determinants of Retinal Cell Fate

The vertebrate retina is regarded as a simple part of the central nervous system (CNS) and thus amenable to investigations of the determinants of cell fate. Its five neuronal cell classes and one glial cell class all derive from a common pool of progenitors. Here we review how each cell class is generated. Retinal progenitors progress through different competence states, in each of which they generate only a small repertoire of cell classes. The intrinsic state of the progenitor is determined by the complement of transcription factors it expresses. Thus, although progenitors are multipotent, there is a bias in the types of fates they generate during any particular time window. Overlying these competence states are stochastic mechanisms that influence fate decisions. These mechanisms are determined by a weighted set of probabilities based on the abundance of a cell class in the retina. Deterministic mechanisms also operate, especially late in development, when preprogrammed progenitors solely generate specific fates.

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

Recent advances in high throughput single-cell RNA sequencing (scRNA-seq) technology have enabled the simultaneous transcriptomic profiling of thousands of individual cells in a single experiment. To investigate the intrinsic process of retinal development, researchers have leveraged this technology to quantify gene expression in retinal cells across development, in multiple species, and from numerous important models of human disease. In this review, we summarize recent applications of scRNA-seq and discuss how these datasets have complemented and advanced our understanding of retinal progenitor cell competence, cell fate specification, and differentiation. Finally, we also highlight the outstanding questions in the field that advances in single-cell data generation and analysis will soon be able to answer.

Development and diversification of bipolar interneurons in the mammalian retina

The bipolar interneurons of the mammalian retina have evolved as a diverse set of cells with distinct subtype characteristics, which reflect specialized contributions to visual circuitry. Fifteen subtypes of bipolar interneurons have been identified in the mouse retina, each with characteristic gene expression, morphology, and light responses. This review provides an overview of the developmental events that underlie the generation of the diverse bipolar cell class, summarizing the current knowledge of genetic programs that establish and maintain bipolar subtype fates, as well as the events that shape the final distribution of bipolar subtypes. With much left to be discovered, bipolar interneurons present an ideal model system for studying the interplay between cell-autonomous and non-cell-autonomous mechanisms that influence neuronal subtype development within the central nervous system.

Regulatory Mechanisms of Retinal Photoreceptors Development at Single Cell Resolution

Photoreceptors are critical components of the retina and play a role in the first step of the conversion of light to electric signals. With the discovery of the intrinsically photosensitive retinal ganglion cells, which regulate non-image-forming visual processes, our knowledge of the photosensitive cell family in the retina has deepened. Photoreceptor development is regulated by specific genes and proteins and involves a series of molecular processes including DNA transcription, post-transcriptional modification, protein translation, and post-translational modification. Single-cell sequencing is a promising technology for the study of photoreceptor development. This review presents an overview of the types of human photoreceptors, summarizes recent discoveries in the regulatory mechanisms underlying their development at single-cell resolution, and outlines the prospects in this field.

Notch signaling represses cone photoreceptor formation through the regulation of retinal progenitor cell states

Notch signaling is required to repress the formation of vertebrate cone photoreceptors and to maintain the proliferative potential of multipotent retinal progenitor cells. However, the mechanism by which Notch signaling controls these processes is unknown. Recently, restricted retinal progenitor cells with limited proliferation capacity and that preferentially generate cone photoreceptors have been identified. Thus, there are several potential steps during cone genesis that Notch signaling could act. Here we use cell type specific cis-regulatory elements to localize the primary role of Notch signaling in cone genesis to the formation of restricted retinal progenitor cells from multipotent retinal progenitor cells. Localized inhibition of Notch signaling in restricted progenitor cells does not alter the number of cones derived from these cells. Cell cycle promotion is not a primary effect of Notch signaling but an indirect effect on progenitor cell state transitions that leads to depletion of the multipotent progenitor cell population. Taken together, this suggests that the role of Notch signaling in cone photoreceptor formation and proliferation are both mediated by a localized function of Notch in multipotent retinal progenitor cells to repress the formation of restricted progenitor cells.

The Role of Small Molecules and Their Effect on the Molecular Mechanisms of Early Retinal Organoid Development

Early in vivo embryonic retinal development is a well-documented and evolutionary conserved process. The specification towards eye development is temporally controlled by consecutive activation or inhibition of multiple key signaling pathways, such as the Wnt and hedgehog signaling pathways. Recently, with the use of retinal organoids, researchers aim to manipulate these pathways to achieve better human representative models for retinal development and disease. To achieve this, a plethora of different small molecules and signaling factors have been used at various time points and concentrations in retinal organoid differentiations, with varying success. Additions differ from protocol to protocol, but their usefulness or efficiency has not yet been systematically reviewed. Interestingly, many of these small molecules affect the same and/or multiple pathways, leading to reduced reproducibility and high variability between studies. In this review, we make an inventory of the key signaling pathways involved in early retinogenesis and their effect on the development of the early retina in vitro. Further, we provide a comprehensive overview of the small molecules and signaling factors that are added to retinal organoid differentiation protocols, documenting the molecular and functional effects of these additions. Lastly, we comparatively evaluate several of these factors using our established retinal organoid methodology.

Transcriptional profiling of murine retinas undergoing semi-synchronous cone photoreceptor differentiation

Uncovering the gene regulatory networks that control cone photoreceptor formation has been hindered because cones only make up a few percent of the retina and form asynchronously during development. To overcome these limitations, we used a γ-secretase inhibitor, DAPT, to disrupt Notch signaling and force proliferating retinal progenitor cells to rapidly adopt neuronal identity. We treated mouse retinal explants at the peak of cone genesis with DAPT and examined tissues at several time-points by histology and bulk RNA-sequencing. We found that this treatment caused supernumerary cone formation in an overwhelmingly synchronized fashion. This analysis revealed several categorical patterns of gene expression changes over time relative to DMSO treated control explants. These were placed in the temporal context of the activation of Otx2, a transcription factor that is expressed at the onset of photoreceptor development and that is required for both rod and cone formation. One group of interest had genes, such as Mybl1, Ascl1, Neurog2, and Olig2, that became upregulated by DAPT treatment before Otx2. Two other groups showed upregulated gene expression shortly after Otx2, either transiently or permanently. This included genes such as Mybl1, Meis2, and Podxl. Our data provide a developmental timeline of the gene expression events that underlie the initial steps of cone genesis and maturation. Applying this strategy to human retinal organoid cultures was also sufficient to induce a massive increase in cone genesis. Taken together, our results provide a temporal framework that can be used to elucidate the gene regulatory logic controlling cone photoreceptor development.

A Novel Approach to Single Cell RNA-Sequence Analysis Facilitates In Silico Gene Reporting of Human Pluripotent Stem Cell-Derived Retinal Cell Types

Cell type-specific investigations commonly use gene reporters or single-cell analytical techniques. However, reporter line development is arduous and generally limited to a single gene of interest, while single-cell RNA (scRNA)-sequencing (seq) frequently yields equivocal results that preclude definitive cell identification. To examine gene expression profiles of multiple retinal cell types derived from human pluripotent stem cells (hPSCs), we performed scRNA-seq on optic vesicle (OV)-like structures cultured under cGMP-compatible conditions. However, efforts to apply traditional scRNA-seq analytical methods based on unbiased algorithms were unrevealing. Therefore, we developed a simple, versatile, and universally applicable approach that generates gene expression data akin to those obtained from reporter lines. This method ranks single cells by expression level of a bait gene and searches the transcriptome for genes whose cell-to-cell rank order expression most closely matches that of the bait. Moreover, multiple bait genes can be combined to refine datasets. Using this approach, we provide further evidence for the authenticity of hPSC-derived retinal cell types. Stem Cells 2018;36:313-324.

Photoreceptor Fate Determination in the Vertebrate Retina

Photoreceptors are highly specialized primary sensory neurons that sense light and initiate vision. This critical role is well demonstrated by the fact that visual impairment accompanies photoreceptor loss or dysfunction in many human diseases. With the remarkable advances in stem cell research, one therapeutic approach is to use stem cells to generate photoreceptors and then engraft them into diseased eyes. Knowledge of the molecular mechanisms that control photoreceptor genesis during normal development can greatly aid in the production of photoreceptor cells for this approach. This article will discuss advances in our understanding of the molecular mechanisms that regulate photoreceptor fate determination during development. Recent lineage studies have shown that there are distinct retinal progenitor cells (RPCs) that produce specific combinations of daughter cell types, including photoreceptors and other types of retinal cells. Gene regulatory networks, in which transcription factors interact via cis-regulatory DNA elements, have been discovered that operate within distinct RPCs, and/or newly postmitotic cells, to direct the choice of photoreceptor fate.

Brg1 coordinates multiple processes during retinogenesis and is a tumor suppressor in retinoblastoma

Retinal development requires precise temporal and spatial coordination of cell cycle exit, cell fate specification, cell migration and differentiation. When this process is disrupted, retinoblastoma, a developmental tumor of the retina, can form. Epigenetic modulators are central to precisely coordinating developmental events, and many epigenetic processes have been implicated in cancer. Studying epigenetic mechanisms in development is challenging because they often regulate multiple cellular processes; therefore, elucidating the primary molecular mechanisms involved can be difficult. Here we explore the role of Brg1 (Smarca4) in retinal development and retinoblastoma in mice using molecular and cellular approaches. Brg1 was found to regulate retinal size by controlling cell cycle length, cell cycle exit and cell survival during development. Brg1 was not required for cell fate specification but was required for photoreceptor differentiation and cell adhesion/polarity programs that contribute to proper retinal lamination during development. The combination of defective cell differentiation and lamination led to retinal degeneration in Brg1-deficient retinae. Despite the hypocellularity, premature cell cycle exit, increased cell death and extended cell cycle length, retinal progenitor cells persisted in Brg1-deficient retinae, making them more susceptible to retinoblastoma. ChIP-Seq analysis suggests that Brg1 might regulate gene expression through multiple mechanisms.

Summary: The SWI/SNF protein Brg1 controls cell cycle length, cell cycle exit and cell survival, and is required for cell differentiation and retinal lamination, in the developing mouse retina.

The Determination of Rod and Cone Photoreceptor Fate

Photoreceptors have been the most intensively studied retinal cell type. Early lineage studies showed that photoreceptors are produced by retinal progenitor cells (RPCs) that produce only photoreceptor cells and by RPCs that produce both photoreceptor cells and other retinal cell types. More recent lineage studies have shown that there are intrinsic, molecular differences among these RPCs and that these molecular differences operate in gene regulatory networks (GRNs) that lead to the choice of the rod versus the cone fate. In addition, there are GRNs that lead to the choice of a photoreceptor fate and that of another retinal cell type. An example of such a GRN is one that drives the binary fate choice between a rod photoreceptor and bipolar cell. This GRN has many elements, including both feedforward and feedback regulatory loops, highlighting the complexity of such networks. This and other examples of retinal cell fate determination are reviewed here, focusing on the events that direct the choice of rod and cone photoreceptor fate.

Gene expression changes during retinal development and rod specification

PURPOSE

Retinitis pigmentosa (RP) typically results from individual mutations in any one of >70 genes that cause rod photoreceptor cells to degenerate prematurely, eventually resulting in blindness. Gene therapies targeting individual RP genes have shown efficacy at clinical trial; however, these therapies require the surviving photoreceptor cells to be viable and functional, and may be economically feasible for only the more commonly mutated genes. An alternative potential treatment strategy, particularly for late stage disease, may involve stem cell transplants into the photoreceptor layer of the retina. Rod progenitors from postnatal mouse retinas can be transplanted and can form photoreceptors in recipient adult retinas; optimal numbers of transplantable cells are obtained from postnatal day 3-5 (P3-5) retinas. These cells can also be expanded in culture; however, this results in the loss of photoreceptor potential. Gene expression differences between postnatal retinas, cultured retinal progenitor cells (RPCs), and rod photoreceptor precursors were investigated to identify gene expression patterns involved in the specification of rod photoreceptors.

METHODS

Microarrays were used to investigate differences in gene expression between cultured RPCs that have lost photoreceptor potential, P1 retinas, and fresh P5 retinas that contain significant numbers of transplantable photoreceptors. Additionally, fluorescence-activated cell sorting (FACS) sorted Rho-eGFP-expressing rod photoreceptor precursors were compared with Rho-eGFP-negative cells from the same P5 retinas. Differential expression was confirmed with quantitative polymerase chain reaction (q-PCR).

RESULTS

Analysis of the microarray data sets, including the use of t-distributed stochastic neighbor embedding (t-SNE) to identify expression pattern neighbors of key photoreceptor specific genes, resulted in the identification of 636 genes differentially regulated during rod specification. Forty-four of these genes when mutated have previously been found to cause retinal disease. Although gene function in other tissues may be known, the retinal function of approximately 61% of the gene list is as yet undetermined. Many of these genes' promoters contain binding sites for the key photoreceptor transcription factors Crx and Nr2e3; moreover, the genomic clustering of differentially regulated genes appears to be non-random.

CONCLUSIONS

This study aids in understanding gene expression differences between rod photoreceptor progenitors versus cultured RPCs that have lost photoreceptor potential. The results provide insights into rod photoreceptor development and should expedite the development of cell-based treatments for RP. Furthermore, the data set includes a large number of retinopathy genes; less-well-characterized genes within this data set are a resource for those seeking to identify novel retinopathy genes in patients with RP (GEO accession: GSE59201).

Transcriptional profiling at whole population and single cell levels reveals somatosensory neuron molecular diversity

The somatosensory nervous system is critical for the organism's ability to respond to mechanical, thermal, and nociceptive stimuli. Somatosensory neurons are functionally and anatomically diverse but their molecular profiles are not well-defined. Here, we used transcriptional profiling to analyze the detailed molecular signatures of dorsal root ganglion (DRG) sensory neurons. We used two mouse reporter lines and surface IB4 labeling to purify three major non-overlapping classes of neurons: 1) IB4⁺SNS-Cre/TdTomato⁺, 2) IB4⁻SNS-Cre/TdTomato⁺, and 3) Parv-Cre/TdTomato⁺ cells, encompassing the majority of nociceptive, pruriceptive, and proprioceptive neurons. These neurons displayed distinct expression patterns of ion channels, transcription factors, and GPCRs. Highly parallel qRT-PCR analysis of 334 single neurons selected by membership of the three populations demonstrated further diversity, with unbiased clustering analysis identifying six distinct subgroups. These data significantly increase our knowledge of the molecular identities of known DRG populations and uncover potentially novel subsets, revealing the complexity and diversity of those neurons underlying somatosensation.

DOI:http://dx.doi.org/10.7554/eLife.04660.001

In the nervous system, a network of specialized neurons—known as the somatosensory system—carries information about sensations including touch, muscle position, temperature and pain. Distinct sets of somatosensory neurons are thought to carry information about the different types of sensations. In young animals, the precise switching on, or ‘expression’, of genes controls the formation of the network of neurons. However, it is not known exactly which genes are expressed in what types of neurons, where, or when.

Here, Chiu et al. used a technique called flow cytometry using different fluorescent markers to isolate a group of cells called Dorsal Root Ganglion (DRG) neurons in mice. These neurons have long thread-like fibers that extend from the spinal cord to the skin, muscles and joints all over the body. These fibers carry sensory information to the spinal cord, where it can be relayed to the brain and processed. The experiments compared three distinct types of DRG neuron and found that they differed in their ability to send information to other cells.

Chiu et al. analyzed the expression of all the genes in the three types of DRG neurons. Each type of neuron had distinct groups of genes that were being expressed. Also, several genes that are known to be important for sensation were expressed at different levels in the different types of cells. Next, large numbers of single cells were analyzed to find out the finer details about the three types of neuron. These findings made it possible to further divide the DRG neurons into six distinct subsets that matched previously known groups of somatosensory neurons, and also identified new ones.

Chiu et al.'s findings reveal the complexity and diversity of the neurons involved in carrying information about sensations towards the brain. This is an important step in classifying the nervous system, and uncovers many genes previously not linked to sensation. The next challenges lie in understanding how the expression of these genes in each type of neuron relates to their unique roles.

DOI:http://dx.doi.org/10.7554/eLife.04660.002

A gene regulatory network controls the binary fate decision of rod and bipolar cells in the vertebrate retina

Gene regulatory networks (GRNs) regulate critical events during development. In complex tissues, such as the mammalian central nervous system (CNS), networks likely provide the complex regulatory interactions needed to direct the specification of the many CNS cell types. Here, we dissect a GRN that regulates a binary fate decision between two siblings in the murine retina, the rod photoreceptor and bipolar interneuron. The GRN centers on Blimp1, one of the transcription factors (TFs) that regulates the rod versus bipolar cell fate decision. We identified a cis-regulatory module (CRM), B108, that mimics Blimp1 expression. Deletion of genomic B108 by CRISPR/Cas9 in vivo using electroporation abolished the function of Blimp1. Otx2 and RORβ were found to regulate Blimp1 expression via B108, and Blimp1 and Otx2 were shown to form a negative feedback loop that regulates the level of Otx2, which regulates the production of the correct ratio of rods and bipolar cells.

Gene expression is dynamically regulated in retinal progenitor cells prior to and during overt cellular differentiation

The retina is comprised of one glial and six neuronal populations that are generated from a multipotent pool of retinal progenitor cells (RPCs) during development. To give rise to these different cell types, RPCs undergo temporal identity transitions, displaying distinct gene expression profiles at different stages of differentiation. Little, however, is known about temporal differences in RPC identities prior to the onset of overt cellular differentiation, during the period when a retinal identity is gradually acquired. Here we examined the sequential onset of expression of regional markers (i.e., homeodomain transcription factors) and cell fate determinants (i.e., basic-helix-loop-helix transcription factors and neurogenic genes) in RPCs from the earliest appearance of a morphologically-distinct retina. By performing a comparative analysis of the expression of a panel of 27 homeodomain, basic-helix-loop-helix and Notch pathway genes between embryonic day (E) 8.75 and postnatal day (P) 9, we identified six distinct RPC molecular profiles. At E8.75, the earliest stage assayed, murine RPCs expressed five homeodomain genes and a single neurogenic gene (Pax6, Six3, Six6, Rx, Otx2, Hes1). This early gene expression profile was remarkably similar to that of 'early' RPCs in the amphibian ciliary marginal zone (CMZ), where RPCs are compartmentalised according to developmental stage, and homologs of Pax6, Six3 and Rx are expressed in the 'early' stem cell zone. As development proceeds, expression of additional homeodomain, bHLH and neurogenic genes was gradually initiated in murine RPCs, allowing distinct genetic profiles to also be defined at E9.5, E10.5, E12.5, E15.5 and P0. In addition, RPCs in the postnatal ciliary margin, where retinal stem cells are retained throughout life, displayed a unique molecular signature, expressing all of the early-onset genes as well as several late-onset markers, indicative of a 'mixed' temporal identity. Taken together, the identification of temporal differences in gene expression in mammalian RPCs during pre-neurogenic developmental stages leads to new insights into how regional identities are progressively acquired during development, while comparisons at later stages highlight the dynamic nature of gene expression in temporally distinct RPC pools.

Notch1 is required in newly postmitotic cells to inhibit the rod photoreceptor fate

Several models of cell fate determination can be invoked to explain how single retinal progenitor cells (RPCs) produce different cell types in a terminal division. To gain insight into this process, the effects of the removal of a cell fate regulator, Notch1, were studied in newly postmitotic cells using a conditional allele of Notch1 (N1-CKO) in mice. Almost all newly postmitotic N1-CKO cells became rod photoreceptors, whereas wild-type (WT) cells achieved a variety of fates. Single cell profiling of wild-type and N1-CKO retinal cells transitioning from progenitor to differentiated states revealed differential expression of inhibitor of DNA binding factors Id1 and Id3, as well as Notch-regulated ankyrin repeat protein (Nrarp). Misexpression of Id1 and Id3 was found to be sufficient to drive production of Müller glial cells and/or RPCs. Moreover, Id1 and Id3 were shown to partially rescue the production of bipolar and Müller glial cells in the absence of Notch1 in mitotic and newly postmitotic cells. Misexpression of Nrarp, a downstream target gene and inhibitor of the Notch signaling pathway, resulted in the overproduction of rod photoreceptors at the expense of Müller glial cells. These data demonstrate that cell fate decisions can be made in newly postmitotic retinal cells, and reveal some of the regulators downstream of Notch1 that influence the choice of rod and non-rod fates. Taken together, our results begin to address how different signals downstream from a common pathway lead to different fate outcomes.