Retinoic acid signaling regulates spatiotemporal specification of human green and red cones.

Trichromacy is unique to primates among placental mammals, enabled by blue (short/S), green (medium/M), and red (long/L) cones. In humans, great apes, and Old World monkeys, cones make a poorly understood choice between M and L cone subtype fates. To determine mechanisms specifying M and L cones, we developed an approach to visualize expression of the highly similar M- and L-opsin mRNAs. M-opsin was observed before L-opsin expression during early human eye development, suggesting that M cones are generated before L cones. In adult human tissue, the early-developing central retina contained a mix of M and L cones compared to the late-developing peripheral region, which contained a high proportion of L cones. Retinoic acid (RA)-synthesizing enzymes are highly expressed early in retinal development. High RA signaling early was sufficient to promote M cone fate and suppress L cone fate in retinal organoids. Across a human population sample, natural variation in the ratios of M and L cone subtypes was associated with a noncoding polymorphism in the NR2F2 gene, a mediator of RA signaling. Our data suggest that RA promotes M cone fate early in development to generate the pattern of M and L cones across the human retina.

Synchronized long-read genome, methylome, epigenome, and transcriptome for resolving a Mendelian condition.

Resolving the molecular basis of a Mendelian condition (MC) remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome, and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion, and structural variant calling and diploid de novo genome assembly, and permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility, and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network (UDN) participant with a chromosome X;13 balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1, and RB1) previously associated with single-gene MCs. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included nonsense-mediated decay, fusion transcript formation, enhancer adoption, transcriptional readthrough silencing, and inappropriate X chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.

Developmental adaptation of rod photoreceptor number via photoreception in melanopsin (OPN4) retinal ganglion cells.

Photoreception, a form of sensory experience, is essential for normal development of the mammalian visual system. Detecting photons during development is a prerequisite for visual system function - from vision's first synapse at the cone pedicle and maturation of retinal vascular networks, to transcriptional establishment and maturation of cell types within the visual cortex. Consistent with this theme, we find that the lighting environment regulates developmental rod photoreceptor apoptosis via OPN4-expressing intrinsically photosensitive retinal ganglion cells (ipRGCs). Using a combination of genetics, sensory environment manipulations, and computational approaches, we establish a molecular pathway in which light-dependent glutamate release from ipRGCs is detected via a transiently expressed kainate receptor (GRIK3) in immature rods localized to the inner retina. Communication between ipRGCs and nascent inner retinal rods appears to be mediated by unusual hybrid neurites projecting from ipRGCs that sense light before eye-opening. These structures, previously referred to as outer retinal dendrites (ORDs), span the ipRGC-immature rod distance over the first postnatal week and contain the machinery for sensory detection (melanopsin, OPN4) and axonal/anterograde neurotransmitter release (Synaptophysin, and VGLUT2). Histological and computational assessment of human mid-gestation development reveal conservation of several hallmarks of an ipRGC-to-immature rod pathway, including displaced immature rods, transient GRIK3 expression in the rod lineage, and the presence of ipRGCs with putative neurites projecting deep into the developing retina. Thus, this analysis defines a retinal retrograde signaling pathway that links the sensory environment to immature rods via ipRGC photoreceptors, allowing the visual system to adapt to distinct lighting environments priory to eye-opening.

Single-cell sequencing of individual retinal organoids reveals determinants of cell-fate heterogeneity.

With a critical need for more complete in vitro models of human development and disease, organoids hold immense potential. Their complex cellular composition makes single-cell sequencing of great utility; however, the limitation of current technologies to a handful of treatment conditions restricts their use in screens or studies of organoid heterogeneity. Here, we apply sci-Plex, a single-cell combinatorial indexing (sci)-based RNA sequencing (RNA-seq) multiplexing method to retinal organoids. We demonstrate that sci-Plex and 10× methods produce highly concordant cell-class compositions and then expand sci-Plex to analyze the cell-class composition of 410 organoids upon modulation of critical developmental pathways. Leveraging individual organoid data, we develop a method to measure organoid heterogeneity, and we identify that activation of Wnt signaling early in retinal organoid cultures increases retinal cell classes up to 6 weeks later. Our data show sci-Plex's potential to dramatically scale up the analysis of treatment conditions on relevant human models.

Single-cell sequencing of individual retinal organoids reveals determinants of cell fate heterogeneity.

With a critical need for more complete in vitro models of human development and disease, organoids hold immense potential. Their complex cellular composition makes single-cell sequencing of great utility; however, the limitation of current technologies to a handful of treatment conditions restricts their use in screens or studies of organoid heterogeneity. Here, we apply sci-Plex, a single-cell combinatorial indexing (sci)-based RNA-seq multiplexing method to retinal organoids. We demonstrate that sci-Plex and 10x methods produce highly concordant cell class compositions and then expand sci-Plex to analyze the cell class composition of 410 organoids upon modulation of critical developmental pathways. Leveraging individual organoid data, we develop a method to measure organoid heterogeneity, and we identify that activation of Wnt signaling early in retinal organoid cultures increases retinal cell classes up to six weeks later. Our data show sci-Plex's potential to dramatically scale-up the analysis of treatment conditions on relevant human models.

CRISPR Generated SIX6 and POU4F2 Reporters Allow Identification of Brain and Optic Transcriptional Differences in Human PSC-Derived Organoids.

Human pluripotent stem cells (PSCs) represent a powerful tool to investigate human eye development and disease. When grown in 3D, they can self-assemble into laminar organized retinas; however, variation in the size, shape and composition of individual organoids exists. Neither the microenvironment nor the timing of critical growth factors driving retinogenesis are fully understood. To explore early retinal development, we developed a SIX6-GFP reporter that enabled the systematic optimization of conditions that promote optic vesicle formation. We demonstrated that early hypoxic growth conditions enhanced SIX6 expression and promoted eye formation. SIX6 expression was further enhanced by sequential inhibition of Wnt and activation of sonic hedgehog signaling. SIX6 + optic vesicles showed RNA expression profiles that were consistent with a retinal identity; however, ventral diencephalic markers were also present. To demonstrate that optic vesicles lead to bona fide "retina-like" structures we generated a SIX6-GFP/POU4F2-tdTomato dual reporter line that labeled the entire developing retina and retinal ganglion cells, respectively. Additional brain regions, including the hypothalamus and midbrain-hindbrain (MBHB) territories were identified by harvesting SIX6 + /POU4F2- and SIX6- organoids, respectively. Using RNAseq to study transcriptional profiles we demonstrated that SIX6-GFP and POU4F2-tdTomato reporters provided a reliable readout for developing human retina, hypothalamus, and midbrain/hindbrain organoids.

Human retinal model systems: Strengths, weaknesses, and future directions.

The retina is a complex neuronal structure that converts light energy into visual perception. Many specialized aspects of the primate retina, including a cone rich macula for high acuity vision, ocular size, and cell type diversity are not found in other animal models. In addition, the unique morphologies and distinct laminar positions of cell types found in the retina make this model system ideal for the study of neuronal cell fate specification. Many key early events of human retinal development are inaccessible to investigation as they occur during gestation. For these reasons, it has been necessary to develop retinal model systems to gain insight into human-specific retinal development and disease. Recent advances in culturing retinal tissue have generated new systems for retinal research and have moved us closer to generating effective regenerative therapies for vision loss. Here, we describe the strengths, weaknesses, and future directions for different human retinal model systems including dissociated primary tissue, explanted primary tissue, retinospheres, and stem cell-derived retinal organoids.

Modeling binary and graded cone cell fate patterning in the mouse retina.

Nervous systems are incredibly diverse, with myriad neuronal subtypes defined by gene expression. How binary and graded fate characteristics are patterned across tissues is poorly understood. Expression of opsin photopigments in the cone photoreceptors of the mouse retina provides an excellent model to address this question. Individual cones express S-opsin only, M-opsin only, or both S-opsin and M-opsin. These cell populations are patterned along the dorsal-ventral axis, with greater M-opsin expression in the dorsal region and greater S-opsin expression in the ventral region. Thyroid hormone signaling plays a critical role in activating M-opsin and repressing S-opsin. Here, we developed an image analysis approach to identify individual cone cells and evaluate their opsin expression from immunofluorescence imaging tiles spanning roughly 6 mm along the D-V axis of the mouse retina. From analyzing the opsin expression of ~250,000 cells, we found that cones make a binary decision between S-opsin only and co-expression competent fates. Co-expression competent cells express graded levels of S- and M-opsins, depending nonlinearly on their position in the dorsal-ventral axis. M- and S-opsin expression display differential, inverse patterns. Using these single-cell data, we developed a quantitative, probabilistic model of cone cell decisions in the retinal tissue based on thyroid hormone signaling activity. The model recovers the probability distribution for cone fate patterning in the mouse retina and describes a minimal set of interactions that are necessary to reproduce the observed cell fates. Our study provides a paradigm describing how differential responses to regulatory inputs generate complex patterns of binary and graded cell fates.

Thyroid hormone signaling specifies cone subtypes in human retinal organoids.

The mechanisms underlying specification of neuronal subtypes within the human nervous system are largely unknown. The blue (S), green (M), and red (L) cones of the retina enable high-acuity daytime and color vision. To determine the mechanism that controls S versus L/M fates, we studied the differentiation of human retinal organoids. Organoids and retinas have similar distributions, expression profiles, and morphologies of cone subtypes. S cones are specified first, followed by L/M cones, and thyroid hormone signaling controls this temporal switch. Dynamic expression of thyroid hormone-degrading and -activating proteins within the retina ensures low signaling early to specify S cones and high signaling late to produce L/M cones. This work establishes organoids as a model for determining mechanisms of human development with promising utility for therapeutics and vision repair.

Left Habenular Activity Attenuates Fear Responses in Larval Zebrafish.

Fear responses are defensive states that ensure survival of an organism in the presence of a threat. Perception of an aversive cue causes changes in behavior and physiology, such as freezing and elevated cortisol, followed by a return to the baseline state when the threat is evaded [1]. Neural systems that elicit fear behaviors include the amygdala, hippocampus, and medial prefrontal cortex. However, aside from a few examples, little is known about brain regions that promote recovery from an aversive event [2]. Previous studies had implicated the dorsal habenular nuclei in regulating fear responses and boldness in zebrafish [3-7]. We now show, through perturbation of its inherent left-right (L-R) asymmetry at larval stages, that the dorsal habenulo-interpeduncular (dHb-IPN) pathway expedites the return of locomotor activity following an unexpected negative stimulus, electric shock. Severing habenular efferents to the IPN, or only those from the left dHb, prolongs the freezing behavior that follows shock. Individuals with a symmetric, right-isomerized dHb also exhibit increased freezing. In contrast, larvae that have a symmetric, left-isomerized dHb, or in which just the left dHb-IPN projection is optogenetically activated, rapidly resume swimming post shock. In vivo calcium imaging reveals a neuronal subset, predominantly in the left dHb, whose activation is correlated with resumption of swimming. The results demonstrate functional specialization of the left dHb-IPN pathway in attenuating the response to fear.

Factors regulating excystment of Alexandrium in Puget Sound, WA, USA.

Factors regulating excystment of a toxic dinoflagellate in the genus Alexandrium were investigated in cysts from Puget Sound, Washington State, USA. Experiments were carried out in the laboratory using cysts collected from benthic seedbeds to determine if excystment is controlled by internal or environmental factors. The results suggest that the timing of germination is not tightly controlled by an endogenous clock, though there is a suggestion of a cyclical pattern. This was explored using cysts that had been stored under cold (4 °C), anoxic conditions in the dark and then incubated for 6 weeks at constant favorable environmental conditions. Excystment occurred during all months of the year, with variable excystment success ranging from 31-90%. When cysts were isolated directly from freshly collected sediments every month and incubated at the in situ bottom water temperature, a seasonal pattern in excystment was observed that was independent of temperature. This pattern may be consistent with secondary dormancy, an externally modulated pattern that prevents excystment during periods that are not favorable for sustained vegetative growth. However, observation over more annual cycles is required and the duration of the mandatory dormancy period of these cysts must be determined before the seasonality of germination can be fully characterized in Alexandrium from Puget Sound. Both temperature and light were found to be important environmental factors regulating excystment, with the highest rates of excystment observed for the warmest temperature treatment (20 °C) and in the light.

Amphetamine-induced sensitization has little effect on multiple learning paradigms and fails to rescue mice with a striatal learning defect.

Behavioral sensitization to psychostimulants such as amphetamine (AMPH) is associated with synaptic modifications that are thought to underlie learning and memory. Because AMPH enhances extracellular dopamine in the striatum where dopamine and glutamate signaling are essential for learning, one might expect that the molecular and morphological changes that occur in the striatum in response to AMPH, including changes in synaptic plasticity, would affect learning. To ascertain whether AMPH sensitization affects learning, we tested wild-type mice and mice lacking NMDA receptor signaling in striatal medium spiny neurons in several different learning tests (motor learning, Pavlovian association, U-maze escape test with strategy shifting) with or without prior sensitization to AMPH. Prior sensitization had minimal effect on learning in any of these paradigms in wild-type mice and failed to restore learning in mutant mice, despite the fact that the mutant mice became sensitized by the AMPH treatment. We conclude that the changes in synaptic plasticity and many other signaling events that occur in response to AMPH sensitization are dissociable from those involved in learning the tasks used in our experiments.

Severely impaired learning and altered neuronal morphology in mice lacking NMDA receptors in medium spiny neurons.

The striatum is composed predominantly of medium spiny neurons (MSNs) that integrate excitatory, glutamatergic inputs from the cortex and thalamus, and modulatory dopaminergic inputs from the ventral midbrain to influence behavior. Glutamatergic activation of AMPA, NMDA, and metabotropic receptors on MSNs is important for striatal development and function, but the roles of each of these receptor classes remain incompletely understood. Signaling through NMDA-type glutamate receptors (NMDARs) in the striatum has been implicated in various motor and appetitive learning paradigms. In addition, signaling through NMDARs influences neuronal morphology, which could underlie their role in mediating learned behaviors. To study the role of NMDARs on MSNs in learning and in morphological development, we generated mice lacking the essential NR1 subunit, encoded by the Grin1 gene, selectively in MSNs. Although these knockout mice appear normal and display normal 24-hour locomotion, they have severe deficits in motor learning, operant conditioning and active avoidance. In addition, the MSNs from these knockout mice have smaller cell bodies and decreased dendritic length compared to littermate controls. We conclude that NMDAR signaling in MSNs is critical for normal MSN morphology and many forms of learning.