Functional and Evolutionary Diversification of Otx2 and Crx in Vertebrate Retinal Photoreceptor and Bipolar Cell Development

Id2 transcriptional repressor promotes retinal Müller glial cell fate choice by increasing Notch activity

In the vertebrate retina, Müller glial cells and neurons originate from common progenitor cells. To identify the genes that regulate the cell fate choice of Müller glial cells from retinal progenitor cells, we first generated Rlbp1-EGFP transgenic mice in which differentiating Müller glia were specifically labeled with EGFP. By performing fluorescence-activated cell sorting and RNA sequencing analysis, we found that Id1, Id2, and Id3, were highly enriched in developing and differentiated Müller glial cells. While Id1/2/3 are inhibitory helix-loop-helix (HLH) transcription factors, the functional role and underlying mechanism of each Id protein in Müller glial cell development are not fully understood. In situ hybridization revealed Id1/2/3 expression in the layer where retinal progenitor cells and Müller glial cells reside in developing and adult retinas, respectively. Cell lineage analysis using retroviruses expressing Id1, Id2, or Id3 showed that Id2 possesses a stronger inducing activity in Müller glial cell fate choice than Id1 and Id3. Co-immunoprecipitation and luciferase reporter analyses suggested that Id2 interacts with the Notch intracellular region and promotes Notch activity. Taken together, these results suggest that Id2 promotes Müller glial cell fate choice by increasing Notch activity in the developing mouse retina.

The advancements in precision medicine for Leber congenital amaurosis: Breakthroughs from genetic diagnosis to therapy

Leber congenital amaurosis (LCA) is a hereditary retinal disease, typically manifesting as severe vision impairment in infancy. With the advancement of precision medicine, genetic diagnosis and targeted therapies offer new hope for LCA patients, significantly improving both diagnostic accuracy and therapeutic efficacy. We summarize the epidemiological characteristics, clinical manifestations, and molecular genetics underlying LCA. It also highlights recent developments in precision treatment strategies, including gene replacement therapy, CRISPR/Cas9-mediated gene editing, and antisense oligonucleotide therapies. In addition, we discuss the applications of induced pluripotent stem cells and retinal organoids in LCA treatment research. Furthermore, we explore preventive strategies and future treatment directions for LCA, including the development of novel gene therapy vectors, the optimization of combinatorial treatment strategies, and the formulation of personalized treatment approaches. These advancements hold significant potential to offer improved treatment options and enhance the quality of life for LCA patients.

Generation of a Double Reporter mES Cell Line to Simultaneously Trace the Generation of Retinal Progenitors and Photoreceptors

Three-dimensional retinal culture systems help to understand eye development and the pathology of disorders. There is a need for reporter stem cell lines to allow in vitro studies on retinal progenitors and photoreceptors and their developmental dynamics or properties and to test therapeutic approaches. The isolation of pure progenitor populations or photoreceptor precursors may serve for drug, gene, and cell therapy development. Here, we generated a dual-reporter mouse embryonic stem cell line Crx-GFP;Rax-mCherry enabling the visualization or isolation of photoreceptors and retinal progenitors from retinal organoid settings. From day 4 organoids, we isolated mCherry-positive cells to assess their early retinal progenitor identity with proliferation tests as well as transcriptomic and proteomic profiling. The timing of eye field transcription factor expression at the transcriptomic and protein levels is in accordance with mouse retinogenesis. This new line will be helpful for rapidly investigating biological questions or testing therapeutics before using human induced pluripotent stem cells (iPSCs), which require a much longer time for retinal organoid formation.

Comprehensive transcription factor perturbations recapitulate fibroblast transcriptional states

Cell atlas projects have nominated recurrent transcriptional states as drivers of biological processes and disease, but their origins, regulation, and properties remain unclear. To enable complementary functional studies, we developed a scalable approach for recapitulating cell states in vitro using CRISPR activation (CRISPRa) Perturb-seq. Aided by a novel multiplexing method, we activated 1,836 transcription factors in two cell types. Measuring 21,958 perturbations showed that CRISPRa activated targets within physiological ranges, that epigenetic features predicted activatable genes, and that the protospacer seed region drove an off-target effect. Perturbations recapitulated in vivo fibroblast states, including universal and inflammatory states, and identified KLF4 and KLF5 as key regulators of the universal state. Inducing the universal state suppressed disease-associated states, highlighting its therapeutic potential. Our findings cement CRISPRa as a tool for perturbing differentiated cells and indicate that in vivo states can be elicited via perturbation, enabling studies of clinically relevant states ex vivo.

Abnormal eyes and spine development in zebrafish (Danio rerio) embryos and larvae induced by triphenyltin

Triphenyltin (TPT) is widely used in crop pest control and ship antifouling coatings, which leads to its entry into aquatic environment and poses a threat to aquatic organisms. However, the effects of TPT on the early life stages of wild fish in natural water environments remains unclear. The aim of this study was to assess the toxic effects of TPT on the early life stages of fish under two different environments: field investigation and laboratory experiment. The occurrence of deformities in wild fish embryos and larvae in the Three Gorges Reservoir (TGR) and the developmental toxicity of TPT at different concentrations (0, 0.15, 1.5 and 15 μg Sn/L) to zebrafish embryos and larvae were observed. The results showed that TPT content was higher in wild larvae, reaching 27.21 ng Sn/g w, and the malformation of wild fish larvae mainly occurred in the eyes and spine under natural water environment. Controlled experiment exposure of zebrafish larvae to TPT also resulted in eye and spinal deformities. Gene expression analysis showed that compared with the control group, the expression levels of genes related to eye development (sox2, otx2, stra6 and rx1) and spine development (sox9a and bmp2b) were significantly up-regulated in the 15 μg Sn/L exposure group, which may be the main cause of eye and spine deformity in the early development stage of fish. In addition, the molecular docking results further elucidate that the strong hydrophobic and electrostatic interactions between TPT and protein residues are the main mechanism of TPT induced abnormal gene expression. Based on these results, it can be inferred that TPT is one of the teratogenic factors of abnormal eye and spine development in the early life stage of fish in the TGR. These findings have important implications for understanding the toxicity of TPT on fish.

Cross-species single-cell landscapes identify the pathogenic gene characteristics of inherited retinal diseases

Introduction

Inherited retinal diseases (IRDs) affect ∼4.5 million people worldwide. Elusive pathogenic variants in over 280 genes are associated with one or more clinical forms of IRDs. It is necessary to understand the complex interaction among retinal cell types and pathogenic genes by constructing a regulatory network. In this study, we attempt to establish a panoramic expression view of the cooperative work in retinal cells to understand the clinical manifestations and pathogenic bases underlying IRDs.

Methods

Single-cell RNA sequencing (scRNA-seq) data on the retinas from 35 retina samples of 3 species (human, mouse, and zebrafish) including 259,087 cells were adopted to perform a comparative analysis across species. Bioinformatic tools were used to conduct weighted gene co-expression network analysis (WGCNA), single-cell regulatory network analysis, cell-cell communication analysis, and trajectory inference analysis.

Results

The cross-species comparison revealed shared or species-specific gene expression patterns at single-cell resolution, such as the stathmin family genes, which were highly expressed specifically in zebrafish Müller glias (MGs). Thirteen gene modules were identified, of which nine were associated with retinal cell types, and Gene Ontology (GO) enrichment of module genes was consistent with cell-specific highly expressed genes. Many IRD genes were identified as hub genes and cell-specific regulons. Most IRDs, especially the retinitis pigmentosa (RP) genes, were enriched in rod-specific regulons. Integrated expression and transcription regulatory network genes, such as congenital stationary night blindness (CSNB) genes GRK1, PDE6B, and TRPM1, showed cell-specific expression and transcription characteristics in either rods or bipolar cells (BCs). IRD genes showed evolutionary conservation (GNAT2, PDE6G, and SAG) and divergence (GNAT2, MT-ND4, and PDE6A) along the trajectory of photoreceptors (PRs) among species. In particular, the Leber congenital amaurosis (LCA) gene OTX2 showed high expression at the beginning of the trajectory of both PRs and BCs.

Conclusion

We identified molecular pathways and cell types closely connected with IRDs, bridging the gap between gene expression, genetics, and pathogenesis. The IRD genes enriched in cell-specific modules and regulons suggest that these diseases share common etiological bases. Overall, mining of interspecies transcriptome data reveals conserved transcriptomic features of retinas across species and promising applications in both normal retina anatomy and retina pathology.

Disease modeling and pharmacological rescue of autosomal dominant retinitis pigmentosa associated with RHO copy number variation

Retinitis pigmentosa (RP), a heterogenous group of inherited retinal disorder, causes slow progressive vision loss with no effective treatments available. Mutations in the rhodopsin gene (RHO) account for ~25% cases of autosomal dominant RP (adRP). In this study, we describe the disease characteristics of the first-ever reported mono-allelic copy number variation (CNV) in RHO as a novel cause of adRP. We (a) show advanced retinal degeneration in a male patient (68 years of age) harboring four transcriptionally active intact copies of rhodopsin, (b) recapitulated the clinical phenotypes using retinal organoids, and (c) assessed the utilization of a small molecule, Photoregulin3 (PR3), as a clinically viable strategy to target and modify disease progression in RP patients associated with RHO-CNV. Patient retinal organoids showed photoreceptors dysgenesis, with rod photoreceptors displaying stunted outer segments with occasional elongated cilia-like projections (microscopy); increased RHO mRNA expression (quantitative real-time PCR [qRT-PCR] and bulk RNA sequencing); and elevated levels and mislocalization of rhodopsin protein (RHO) within the cell body of rod photoreceptors (western blotting and immunohistochemistry) over the extended (300 days) culture time period when compared against control organoids. Lastly, we utilized PR3 to target NR2E3, an upstream regulator of RHO, to alter RHO expression and observed a partial rescue of RHO protein localization from the cell body to the inner/outer segments of rod photoreceptors in patient organoids. These results provide a proof-of-principle for personalized medicine and suggest that RHO expression requires precise control. Taken together, this study supports the clinical data indicating that RHO-CNV associated adRPdevelops as a result of protein overexpression, thereby overloading the photoreceptor post-translational modification machinery.

Homeobox gene-encoded transcription factors in development and mature circadian function of the rodent pineal gland

Homeobox genes encode transcription factors that are widely known to control developmental processes. This is also the case in the pineal gland, a neuroendocrine brain structure devoted to nighttime synthesis of the hormone melatonin. Thus, in accordance with high prenatal gene expression, knockout studies have identified a specific set of homeobox genes that are essential for development of the pineal gland. However, as a special feature of the pineal gland, homeobox gene expression persists into adulthood, and gene product abundance exhibits 24 h circadian rhythms. Recent lines of evidence show that some homeobox genes even control expression of enzymes catalyzing melatonin synthesis. We here review current knowledge of homeobox genes in the rodent pineal gland and suggest a model for dual functions of homeobox gene-encoded transcription factors in developmental and circadian mature neuroendocrine function.

Pathogenic variants in CRX have distinct cis-regulatory effects on enhancers and silencers in photoreceptors

Dozens of variants in the gene for the homeodomain transcription factor (TF) cone-rod homeobox (CRX) are linked with human blinding diseases that vary in their severity and age of onset. How different variants in this single TF alter its function in ways that lead to a range of phenotypes is unclear. We characterized the effects of human disease-causing variants on CRX cis-regulatory function by deploying massively parallel reporter assays (MPRAs) in mouse retina explants carrying knock-ins of two variants, one in the DNA-binding domain (p.R90W) and the other in the transcriptional effector domain (p.E168d2). The degree of reporter gene dysregulation in these mutant Crx retinas corresponds with their phenotypic severity. The two variants affect similar sets of enhancers, and p.E168d2 has distinct effects on silencers. Cis-regulatory elements (CREs) near cone photoreceptor genes are enriched for silencers that are derepressed in the presence of p.E168d2. Chromatin environments of CRX-bound loci are partially predictive of episomal MPRA activity, and distal elements whose accessibility increases later in retinal development are enriched for CREs with silencer activity. We identified a set of potentially pleiotropic regulatory elements that convert from silencers to enhancers in retinas that lack a functional CRX effector domain. Our findings show that phenotypically distinct variants in different domains of CRX have partially overlapping effects on its cis-regulatory function, leading to misregulation of similar sets of enhancers while having a qualitatively different impact on silencers.

Transplanted human photoreceptors transfer cytoplasmic material but not to the recipient mouse retina

BACKGROUND

The discovery of material transfer between transplanted and host mouse photoreceptors has expanded the possibilities for utilizing transplanted photoreceptors as potential vehicles for delivering therapeutic cargo. However, previous research has not directly explored the capacity for human photoreceptors to engage in material transfer, as human photoreceptor transplantation has primarily been investigated in rodent models of late-stage retinal disease, which lack host photoreceptors.

METHODS

In this study, we transplanted human stem-cell derived photoreceptors purified from human retinal organoids at different ontological ages (weeks 10, 14, or 20) into mouse models with intact photoreceptors and assessed transfer of human proteins and organelles to mouse photoreceptors.

RESULTS

Unexpectedly, regardless of donor age or mouse recipient background, human photoreceptors did not transfer material in the mouse retina, though a rare subset of donor cells (< 5%) integrated into the mouse photoreceptor cell layer. To investigate the possibility that a species barrier impeded transfer, we used a flow cytometric assay to examine material transfer in vitro. Interestingly, dissociated human photoreceptors transferred fluorescent protein with each other in vitro, yet no transfer was detected in co-cultures of human and mouse photoreceptors, suggesting that material transfer is species specific.

CONCLUSIONS

While xenograft models are not a tractable system to study material transfer of human photoreceptors, these findings demonstrate that human retinal organoid-derived photoreceptors are competent donors for material transfer and thus may be useful to treat retinal degenerative disease.

Transcriptional precision in photoreceptor development and diseases - Lessons from 25 years of CRX research

The vertebrate retina is made up of six specialized neuronal cell types and one glia that are generated from a common retinal progenitor. The development of these distinct cell types is programmed by transcription factors that regulate the expression of specific genes essential for cell fate specification and differentiation. Because of the complex nature of transcriptional regulation, understanding transcription factor functions in development and disease is challenging. Research on the Cone-rod homeobox transcription factor CRX provides an excellent model to address these challenges. In this review, we reflect on 25 years of mammalian CRX research and discuss recent progress in elucidating the distinct pathogenic mechanisms of four CRX coding variant classes. We highlight how in vitro biochemical studies of CRX protein functions facilitate understanding CRX regulatory principles in animal models. We conclude with a brief discussion of the emerging systems biology approaches that could accelerate precision medicine for CRX-linked diseases and beyond.

Transcription factor interactions explain the context-dependent activity of CRX binding sites

The effects of transcription factor binding sites (TFBSs) on the activity of a cis-regulatory element (CRE) depend on the local sequence context. In rod photoreceptors, binding sites for the transcription factor (TF) Cone-rod homeobox (CRX) occur in both enhancers and silencers, but the sequence context that determines whether CRX binding sites contribute to activation or repression of transcription is not understood. To investigate the context-dependent activity of CRX sites, we fit neural network-based models to the activities of synthetic CREs composed of photoreceptor TFBSs. The models revealed that CRX binding sites consistently make positive, independent contributions to CRE activity, while negative homotypic interactions between sites cause CREs composed of multiple CRX sites to function as silencers. The effects of negative homotypic interactions can be overcome by the presence of other TFBSs that either interact cooperatively with CRX sites or make independent positive contributions to activity. The context-dependent activity of CRX sites is thus determined by the balance between positive heterotypic interactions, independent contributions of TFBSs, and negative homotypic interactions. Our findings explain observed patterns of activity among genomic CRX-bound enhancers and silencers, and suggest that enhancers may require diverse TFBSs to overcome negative homotypic interactions between TFBSs.

Pathogenic variants in Crx have distinct cis-regulatory effects on enhancers and silencers in photoreceptors

Dozens of variants in the photoreceptor-specific transcription factor (TF) CRX are linked with human blinding diseases that vary in their severity and age of onset. It is unclear how different variants in this single TF alter its function in ways that lead to a range of phenotypes. We examined the effects of human disease-causing variants on CRX cis-regulatory function by deploying massively parallel reporter assays (MPRAs) in live mouse retinas carrying knock-ins of two variants, one in the DNA binding domain (p.R90W) and the other in the transcriptional effector domain (p.E168d2). The degree of reporter gene dysregulation caused by the variants corresponds with their phenotypic severity. The two variants affect similar sets of enhancers, while p.E168d2 has stronger effects on silencers. Cis-regulatory elements (CREs) near cone photoreceptor genes are enriched for silencers that are de-repressed in the presence of p.E168d2. Chromatin environments of CRX-bound loci were partially predictive of episomal MPRA activity, and silencers were notably enriched among distal elements whose accessibility increases later in retinal development. We identified a set of potentially pleiotropic regulatory elements that convert from silencers to enhancers in retinas that lack a functional CRX effector domain. Our findings show that phenotypically distinct variants in different domains of CRX have partially overlapping effects on its cis-regulatory function, leading to misregulation of similar sets of enhancers, while having a qualitatively different impact on silencers.

The Rax homeoprotein in Müller glial cells is required for homeostasis maintenance of the postnatal mouse retina

Müller glial cells, which are the most predominant glial subtype in the retina, play multiple important roles, including the maintenance of structural integrity, homeostasis, and physiological functions of the retina. We have previously found that the Rax homeoprotein is expressed in postnatal and mature Müller glial cells in the mouse retina. However, the function of Rax in postnatal and mature Müller glial cells remains to be elucidated. In the current study, we first investigated Rax function in retinal development using retroviral lineage analysis and found that Rax controls the specification of late-born retinal cell types, including Müller glial cells in the postnatal retina. We next generated Rax tamoxifen-induced conditional KO (Rax iCKO) mice, where Rax can be depleted in mTFP-labeled Müller glial cells upon tamoxifen treatment, by crossing Raxflox/flox mice with Rlbp1-CreERT2 mice, which we have produced. Immunohistochemical analysis showed a characteristic of reactive gliosis and enhanced gliosis of Müller glial cells in Rax iCKO retinas under normal and stress conditions, respectively. We performed RNA-seq analysis on mTFP-positive cells purified from the Rax iCKO retina and found significantly reduced expression of suppressor of cytokinesignaling-3 (Socs3). Reporter gene assays showed that Rax directly transactivates the Socs3 promoter. We observed decreased expression of Socs3 in Müller glial cells of Rax iCKO retinas by immunostaining. Taken together, the present results suggest that Rax suppresses inflammation in Müller glial cells by transactivating Socs3. This study sheds light on the transcriptional regulatory mechanisms underlying retinal Müller glial cell homeostasis.

Single-cell multiomics of the human retina reveals hierarchical transcription factor collaboration in mediating cell type-specific effects of genetic variants on gene regulation

BACKGROUND

Systematic characterization of how  genetic variation modulates gene regulation in a cell type-specific context is essential for understanding complex traits. To address this question, we profile gene expression and chromatin accessibility in cells from healthy retinae of 20 human donors through single-cell multiomics and genomic sequencing.

RESULTS

We map eQTL, caQTL, allelic-specific expression, and allelic-specific chromatin accessibility in major retinal cell types. By integrating these results, we identify and characterize regulatory elements and genetic variants effective on gene regulation in individual cell types. The majority of identified sc-eQTLs and sc-caQTLs display cell type-specific effects, while the cis-elements containing genetic variants with cell type-specific effects are often accessible in multiple cell types. Furthermore, the transcription factors whose binding sites are perturbed by genetic variants tend to have higher expression levels in the cell types where the variants exert their effects, compared to the cell types where the variants have no impact. We further validate our findings with high-throughput reporter assays. Lastly, we identify the enriched cell types, candidate causal variants and genes, and cell type-specific regulatory mechanism underlying GWAS loci.

CONCLUSIONS

Overall, genetic effects on gene regulation are highly context dependent. Our results suggest that cell type-dependent genetic effect is driven by precise modulation of both trans-factor expression and chromatin accessibility of cis-elements. Our findings indicate hierarchical collaboration among transcription factors plays a crucial role in mediating cell type-specific effects of genetic variants on gene regulation.

Disease modeling and pharmacological rescue of autosomal dominant Retinitis Pigmentosa associated with RHO copy number variation

Retinitis pigmentosa (RP), a heterogenous group of inherited retinal disorder causes slow progressive vision loss with no effective treatments available. Mutations in the rhodopsin gene (RHO), account for ~25% cases of autosomal dominant RP (adRP). In this study, we describe the disease characteristics of the first ever reported mono-allelic copy number variation (CNV) in RHO as a novel cause of adRP. We (1) show advanced retinal degeneration in a male patient (60-70 year old) harboring four transcriptionally active intact copies of rhodopsin, (2) recapitulated the clinical phenotypes using retinal organoids, and (3) assessed the utilization of a small molecule, Photoregulin3 (PR3), as a clinically viable strategy to target and modify disease progression in RP patients associated with RHO-CNV. Patient retinal organoids showed photoreceptors dysgenesis, with rod photoreceptors displaying stunted outer segments with occasional elongated cilia-like projections (microscopy); increased RHO mRNA expression (qRT-PCR and bulk RNA-sequencing); and elevated levels and mislocalization of rhodopsin protein (RHO) within the cell body of rod photoreceptors (western blotting and immunohistochemistry) over the extended (300-days) culture time period when compared against control organoids. Lastly, we utilized PR3 to target NR2E3, an upstream regulator of RHO, to alter RHO expression and observed a partial rescue of RHO protein localization from the cell body to the inner/outer segments of rod photoreceptors in patient organoids. These results provide a proof-of-principle for personalized medicine and suggest that RHO expression requires precise control. Taken together, this study supports the clinical data indicating that adRP due to RHO-CNV develops due protein overexpression overloading the photoreceptor post-translational modification machinery.

Reliability of human retina organoid generation from hiPSC-derived neuroepithelial cysts

The possible applications for human retinal organoids (HROs) derived from human induced pluripotent stem cells (hiPSC) rely on the robustness and transferability of the methodology for their generation. Standardized strategies and parameters to effectively assess, compare, and optimize organoid protocols are starting to be established, but are not yet complete. To advance this, we explored the efficiency and reliability of a differentiation method, called CYST protocol, that facilitates retina generation by forming neuroepithelial cysts from hiPSC clusters. Here, we tested seven different hiPSC lines which reproducibly generated HROs. Histological and ultrastructural analyses indicate that HRO differentiation and maturation are regulated. The different hiPSC lines appeared to be a larger source of variance than experimental rounds. Although previous reports have shown that HROs in several other protocols contain a rather low number of cones, HROs from the CYST protocol are consistently richer in cones and with a comparable ratio of cones, rods, and Müller glia. To provide further insight into HRO cell composition, we studied single cell RNA sequencing data and applied CaSTLe, a transfer learning approach. Additionally, we devised a potential strategy to systematically evaluate different organoid protocols side-by-side through parallel differentiation from the same hiPSC batches: In an explorative study, the CYST protocol was compared to a conceptually different protocol based on the formation of cell aggregates from single hiPSCs. Comparing four hiPSC lines showed that both protocols reproduced key characteristics of retinal epithelial structure and cell composition, but the CYST protocol provided a higher HRO yield. So far, our data suggest that CYST-derived HROs remained stable up to at least day 200, while single hiPSC-derived HROs showed spontaneous pathologic changes by day 200. Overall, our data provide insights into the efficiency, reproducibility, and stability of the CYST protocol for generating HROs, which will be useful for further optimizing organoid systems, as well as for basic and translational research applications.

Molecular mechanisms controlling vertebrate retinal patterning, neurogenesis, and cell fate specification

This review covers recent advances in understanding the molecular mechanisms controlling neurogenesis and specification of the developing retina, with a focus on insights obtained from comparative single cell multiomic analysis. We discuss recent advances in understanding the mechanisms by which extrinsic factors trigger transcriptional changes that spatially pattern the optic cup (OC) and control the initiation and progression of retinal neurogenesis. We also discuss progress in unraveling the core evolutionarily conserved gene regulatory networks (GRNs) that specify early- and late-state retinal progenitor cells (RPCs) and neurogenic progenitors and that control the final steps in determining cell identity. Finally, we discuss findings that provide insight into regulation of species-specific aspects of retinal patterning and neurogenesis, including consideration of key outstanding questions in the field.

The importance of thyroid hormone signaling during early development: Lessons from the zebrafish model

The zebrafish is an optimal experimental model to study thyroid hormone (TH) involvement in vertebrate development. The use of state-of-the-art zebrafish genetic tools available for the study of the effect of gene silencing, cell fate decisions and cell lineage differentiation have contributed to a more insightful comprehension of molecular, cellular, and tissue-specific TH actions. In contrast to intrauterine development, extrauterine embryogenesis observed in zebrafish has facilitated a more detailed study of the development of the hypothalamic-pituitary-thyroid axis. This model has also enabled a more insightful analysis of TH molecular actions upon the organization and function of the brain, the retina, the heart, and the immune system. Consequently, zebrafish has become a trendy model to address paradigms of TH-related functional and biomedical importance. We here compilate the available knowledge regarding zebrafish developmental events for which specific components of TH signaling are essential.

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

Cellular and Molecular Mechanisms of Pathogenesis Underlying Inherited Retinal Dystrophies

Inherited retinal dystrophies (IRDs) are congenital retinal degenerative diseases that have various inheritance patterns, including dominant, recessive, X-linked, and mitochondrial. These diseases are most often the result of defects in rod and/or cone photoreceptor and retinal pigment epithelium function, development, or both. The genes associated with these diseases, when mutated, produce altered protein products that have downstream effects in pathways critical to vision, including phototransduction, the visual cycle, photoreceptor development, cellular respiration, and retinal homeostasis. The aim of this manuscript is to provide a comprehensive review of the underlying molecular mechanisms of pathogenesis of IRDs by delving into many of the genes associated with IRD development, their protein products, and the pathways interrupted by genetic mutation.

Disease-causing mutations in genes encoding transcription factors critical for photoreceptor development

Photoreceptor development of the vertebrate visual system is controlled by a complex transcription regulatory network. OTX2 is expressed in the mitotic retinal progenitor cells (RPCs) and controls photoreceptor genesis. CRX that is activated by OTX2 is expressed in photoreceptor precursors after cell cycle exit. NEUROD1 is also present in photoreceptor precursors that are ready to specify into rod and cone photoreceptor subtypes. NRL is required for the rod fate and regulates downstream rod-specific genes including the orphan nuclear receptor NR2E3 which further activates rod-specific genes and simultaneously represses cone-specific genes. Cone subtype specification is also regulated by the interplay of several transcription factors such as THRB and RXRG. Mutations in these key transcription factors are responsible for ocular defects at birth such as microphthalmia and inherited photoreceptor diseases such as Leber congenital amaurosis (LCA), retinitis pigmentosa (RP) and allied dystrophies. In particular, many mutations are inherited in an autosomal dominant fashion, including the majority of missense mutations in CRX and NRL. In this review, we describe the spectrum of photoreceptor defects that are associated with mutations in the above-mentioned transcription factors, and summarize the current knowledge of molecular mechanisms underlying the pathogenic mutations. At last, we deliberate the outstanding gaps in our understanding of the genotype-phenotype correlations and outline avenues for future research of the treatment strategies.

Single-cell transcriptomics of the goldfish retina reveals genetic divergence in the asymmetrically evolved subgenomes after allotetraploidization

The recent whole-genome duplication (WGD) in goldfish (Carassius auratus) approximately 14 million years ago makes it a valuable model for studying gene evolution during the early stages after WGD. We analyzed the transcriptome of the goldfish retina at the level of single-cell (scRNA-seq) and open chromatin regions (scATAC-seq). We identified a group of genes that have undergone dosage selection, accounting for 5% of the total 11,444 ohnolog pairs. We also identified 306 putative sub/neo-functionalized ohnolog pairs that are likely to be under cell-type-specific genetic variation at single-cell resolution. Diversification in the expression patterns of several ohnolog pairs was observed in the retinal cell subpopulations. The single-cell level transcriptome analysis in this study uncovered the early stages of evolution in retinal cell of goldfish after WGD. Our results provide clues for understanding the relationship between the early stages of gene evolution after WGD and the evolution of diverse vertebrate retinal functions.

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

Mutations in BCOR, a co-repressor of CRX/OTX2, are associated with early-onset retinal degeneration

Many transcription factors regulating the production, survival, and function of photoreceptor cells have been identified, but little is known about transcriptional co-regulators in retinal health and disease. Here, we show that BCL6 co-repressor (BCOR), a Polycomb repressive complex 1 factor mutated in various cancers, is involved in photoreceptor degenerative diseases. Using proteomics and transcription assays, we report that BCOR interacts with the transcription factors CRX and OTX2 and reduces their ability to activate the promoters of photoreceptor-specific genes. CUT&RUN sequencing further shows that BCOR shares genome-wide binding profiles with CRX/OTX2, consistent with a general co-repression activity. We also identify missense mutations in human BCOR in five families that have no evidence of cancer but present severe early-onset X-linked retinal degeneration. Last, we show that the human BCOR mutants cause degeneration when expressed in the mouse retina and have enhanced repressive activity on OTX2. These results uncover a role for BCOR in photoreceptors in both health and disease.

Timed Notch Inhibition Drives Photoreceptor Fate Specification in Human Retinal Organoids

Purpose

Transplanting photoreceptors from human pluripotent stem cell-derived retinal organoids have the potential to reverse vision loss in affected individuals. However, transplantable photoreceptors are only a subset of all cells in the organoids. Hence, the goal of our current study was to accelerate and synchronize photoreceptor differentiation in retinal organoids by inhibiting the Notch signaling pathway at different developmental time-points using a small molecule, PF-03084014 (PF).

Methods

Human induced pluripotent stem cell- and human embryonic stem cells-derived retinal organoids were treated with 10 µM PF for 3 days starting at day 45 (D45), D60, D90, and D120 of differentiation. Organoids were collected at post-treatment days 14, 28, and 42 and analyzed for progenitor and photoreceptor markers and Notch pathway inhibition by immunohistochemistry (IHC), quantitative PCR, and bulk RNA sequencing (n = 3-5 organoids from three independent experiments).

Results

Retinal organoids collected after treatment showed a decrease in progenitor markers (KI67, VSX2, PAX6, and LHX2) and an increase in differentiated pan-photoreceptor markers (OTX2, CRX, and RCVRN) at all organoid stages except D120. PF-treated organoids at D45 and D60 exhibited an increase in cone photoreceptor markers (RXRG and ARR3). PF treatment at D90 revealed an increase in cone and rod photoreceptors markers (ARR3, NRL, and NR2E3). Bulk RNA sequencing analysis mirrored the immunohistochemistry data and quantitative PCR confirmed Notch effector inhibition.

Conclusions

Timing the Notch pathway inhibition in human retinal organoids to align with progenitor competency stages can yield an enriched population of early cone or rod photoreceptors.

RB1 loss triggers dependence on ESRRG in retinoblastoma

Retinoblastoma (Rb) is a deadly childhood eye cancer that is classically initiated by inactivation of the RB1 tumor suppressor. Clinical management continues to rely on nonspecific chemotherapeutic agents that are associated with treatment resistance and toxicity. Here, we analyzed 103 whole exomes, 20 whole transcriptomes, 5 single-cell transcriptomes, and 4 whole genomes from primary Rb tumors to identify previously unknown Rb dependencies. Several recurrent genomic aberrations implicate estrogen-related receptor gamma (ESRRG) in Rb pathogenesis. RB1 directly interacts with and inhibits ESRRG, and RB1 loss uncouples ESRRG from negative regulation. ESRRG regulates genes involved in retinogenesis and oxygen metabolism in Rb cells. ESRRG is preferentially expressed in hypoxic Rb cells in vivo. Depletion or inhibition of ESRRG causes marked Rb cell death, which is exacerbated in hypoxia. These findings reveal a previously unidentified dependency of Rb cells on ESRRG, and they implicate ESRRG as a potential therapeutic vulnerability in Rb.

An enhancer located in a Pde6c intron drives transient expression in the cone photoreceptors of developing mouse and human retinas

How cone photoreceptors are formed during retinal development is only partially known. This is in part because we do not fully understand the gene regulatory network responsible for cone genesis. We reasoned that cis-regulatory elements (enhancers) active in nascent cones would be regulated by the same upstream network that controls cone formation. To dissect this network, we searched for enhancers active in developing cones. By electroporating enhancer-driven fluorescent reporter plasmids, we observed that a sequence within an intron of the cone-specific Pde6c gene acted as an enhancer in developing mouse cones. Similar fluorescent reporter plasmids were used to generate stable transgenic human induced pluripotent stem cells that were then grown into three-dimensional human retinal organoids. These organoids contained fluorescently labeled cones, demonstrating that the Pde6c enhancer was also active in human cones. We observed that enhancer activity was transient and labeled a minor population of developing rod photoreceptors in both mouse and human systems. This cone-enriched pattern argues that the Pde6c enhancer is activated in cells poised between rod and cone fates. Additionally, it suggests that the Pde6c enhancer is activated by the same regulatory network that selects or stabilizes cone fate choice. To further understand this regulatory network, we identified essential enhancer sequence regions through a series of mutagenesis experiments. This suggested that the Pde6c enhancer was regulated by transcription factor binding at five or more locations. Binding site predictions implicated transcription factor families known to control photoreceptor formation and families not previously associated with cone development. These results provide a framework for deciphering the gene regulatory network that controls cone genesis in both human and mouse systems. Our new transgenic human stem cell lines provide a tool for determining which cone developmental mechanisms are shared and distinct between mice and humans.

Transcriptomic evidence for visual adaptation during the aquatic to terrestrial metamorphosis in leopard frogs

BACKGROUND

Differences in morphology, ecology, and behavior through ontogeny can result in opposing selective pressures at different life stages. Most animals, however, transition through two or more distinct phenotypic phases, which is hypothesized to allow each life stage to adapt more freely to its ecological niche. How this applies to sensory systems, and in particular how sensory systems adapt across life stages at the molecular level, is not well understood. Here, we used whole-eye transcriptomes to investigate differences in gene expression between tadpole and juvenile southern leopard frogs (Lithobates sphenocephalus), which rely on vision in aquatic and terrestrial light environments, respectively. Because visual physiology changes with light levels, we also tested the effect of light and dark exposure.

RESULTS

We found 42% of genes were differentially expressed in the eyes of tadpoles versus juveniles and 5% for light/dark exposure. Analyses targeting a curated subset of visual genes revealed significant differential expression of genes that control aspects of visual function and development, including spectral sensitivity and lens composition. Finally, microspectrophotometry of photoreceptors confirmed shifts in spectral sensitivity predicted by the expression results, consistent with adaptation to distinct light environments.

CONCLUSIONS

Overall, we identified extensive expression-level differences in the eyes of tadpoles and juveniles related to observed morphological and physiological changes through metamorphosis and corresponding adaptive shifts to improve vision in the distinct aquatic and terrestrial light environments these frogs inhabit during their life cycle. More broadly, these results suggest that decoupling of gene expression can mediate the opposing selection pressures experienced by organisms with complex life cycles that inhabit different environmental conditions throughout ontogeny.

New Insights on the Regulatory Gene Network Disturbed in Central Areolar Choroidal Dystrophy-Beyond Classical Gene Candidates

Central areolar choroidal dystrophy (CACD) is a rare hereditary disease that mainly affects the macula, resulting in progressive and usually profound visual loss. Being part of congenital retinal dystrophies, it may have an autosomal dominant or recessive inheritance and, until now, has no effective treatment. Given the shortage of genotypic information about the disease, this work systematically reviews the literature for CACD-causing genes. Three independent researchers selected 33 articles after carefully searching and filtering the Scielo, Pubmed, Lilacs, Web of Science, Scopus, and Embase databases. Mutations of six genes (PRPH2, GUCA1A, GUCY2D, CDHR1, ABCA4, and TTLL5) are implicated in the monogenic dominant inheritance of CACD. They are functionally related to photoreceptors (either in the phototransduction process, as in the case of GUCY2D, or the recovery of retinal photodegradation in photoreceptors for GUCA1A, or the formation and maintenance of specific structures within photoreceptors for PRPH2). The identified genetic variants do not explain all observed clinical features, calling for further whole-genome and functional studies for this disease. A network analysis with the CACD-related genes identified in the systematic review resulted in the identification of another 20 genes that may influence CACD onset and symptoms. Furthermore, an enrichment analysis allowed the identification of 13 transcription factors and 4 long noncoding RNAs interacting with the products of the previously mentioned genes. If mutated or dysregulated, they may be directly involved in CACD development and related disorders. More than half of the genes identified by bioinformatic tools do not appear in commercial gene panels, calling for more studies about their role in the maintenance of the retina and phototransduction process, as well as for a timely update of these gene panels.

Generation of an RCVRN-eGFP Reporter hiPSC Line by CRISPR/Cas9 to Monitor Photoreceptor Cell Development and Facilitate the Cell Enrichment for Transplantation

Stem cell-based cell therapies are considered to be promising treatments for retinal disorders with dysfunction or death of photoreceptors. However, the enrichment of human photoreceptors suitable for transplantation has been highly challenging so far. This study aimed to generate a photoreceptor-specific reporter human induced pluripotent stem cell (hiPSC) line using CRISPR/Cas9 genome editing, which harbored an enhanced green fluorescent protein (eGFP) sequence at the endogenous locus of the pan photoreceptor marker recoverin (RCVRN). After confirmation of successful targeting and gene stability, three-dimensional retinal organoids were induced from this reporter line. The RCVRN-eGFP reporter faithfully replicated endogenous protein expression of recoverin and revealed the developmental characteristics of photoreceptors during retinal differentiation. The RCVRN-eGFP specifically and steadily labeled photoreceptor cells from photoreceptor precursors to mature rods and cones. Additionally, abundant eGFP-positive photoreceptors were enriched by fluorescence-activated cell sorting, and their transcriptome signatures were revealed by RNA sequencing and data analysis. Moreover, potential clusters of differentiation (CD) biomarkers were extracted for the enrichment of photoreceptors for clinical applications, such as CD133 for the positive selection of photoreceptors. Altogether, the RCVRN-eGFP reporter hiPSC line was successfully established and the first global expression database of recoverin-positive photoreceptors was constructed. These achievements will provide a powerful tool for dynamically monitoring photoreceptor cell development and purification of human photoreceptors, thus facilitating photoreceptor cell therapy for advanced retinal disorders.

Rod genesis driven by mafba in an nrl knockout zebrafish model with altered photoreceptor composition and progressive retinal degeneration

Neural retina leucine zipper (NRL) is an essential gene for the fate determination and differentiation of the precursor cells into rod photoreceptors in mammals. Mutations in NRL are associated with the autosomal recessive enhanced S-cone syndrome and autosomal dominant retinitis pigmentosa. However, the exact role of Nrl in regulating the development and maintenance of photoreceptors in the zebrafish (Danio rerio), a popular animal model used for retinal degeneration and regeneration studies, has not been fully determined. In this study, we generated an nrl knockout zebrafish model via the CRISPR-Cas9 technology and observed a surprising phenotype characterized by a reduced number, but not the total loss, of rods and over-growth of green cones. We discovered two waves of rod genesis, nrl-dependent and -independent at the embryonic and post-embryonic stages, respectively, in zebrafish by monitoring the rod development. Through bulk and single-cell RNA sequencing, we characterized the gene expression profiles of the whole retina and each retinal cell type from the wild type and nrl knockout zebrafish. The over-growth of green cones and mis-expression of green-cone-specific genes in rods in nrl mutants suggested that there are rod/green-cone bipotent precursors, whose fate choice between rod versus green-cone is controlled by nrl. Besides, we identified the mafba gene as a novel regulator of the nrl-independent rod development, based on the cell-type-specific expression patterns and the retinal phenotype of nrl/mafba double-knockout zebrafish. Gene collinearity analysis revealed the evolutionary origin of mafba and suggested that the function of mafba in rod development is specific to modern fishes. Furthermore, the altered photoreceptor composition and abnormal gene expression in nrl mutants caused progressive retinal degeneration and subsequent regeneration. Accordingly, this study revealed a novel function of the mafba gene in rod development and established a working model for the developmental and regulatory mechanisms regarding the rod and green-cone photoreceptors in zebrafish.

Vision is mediated by two types of light-sensing cells named rod and cone photoreceptors in animal eyes. Abnormal generation, dysfunction or death of photoreceptor cells all cause irreversible vision problems. NRL is an essential gene for the generation and function of rod cells in mice and humans. Surprisingly, we found that in the zebrafish, a popular animal model for human diseases and therapeutic testing, there are two types of rod cells, and eliminating the function of nrl gene affects the rod cell formation at the embryonic stage but not at the juvenile and adult stages. The rod cell formation at the post-embryonic is driven by the mafba gene, which has not been reported to play a role in rod cells. In addition to the reduced number of rod cells, deletion of nrl also results in the emergence of rod/green-cone hybrid cells and an increased number of green cones. The ensuing cellular and molecular alterations collectively lead to retinal degeneration. These findings expand our understanding of photoreceptor development and maintenance and highlight the underlying conserved and species-specific regulatory mechanisms.

The RNA-binding protein and stress granule component ATAXIN-2 is expressed in mouse and human tissues associated with glaucoma pathogenesis

Polyglutamine repeat expansions in the Ataxin-2 (ATXN2) gene were first implicated in Spinocerebellar Ataxia Type 2, a disease associated with degeneration of motor neurons and Purkinje cells. Recent studies linked single nucleotide polymorphisms in the gene to elevated intraocular pressure in primary open angle glaucoma (POAG); yet, the localization of ATXN2 across glaucoma-relevant tissues of the vertebrate eye has not been thoroughly examined. This study characterizes ATXN2 expression in the mouse and human retina, and anterior eye, using an antibody validated in ATXN2^-/- retinas. ATXN2-ir was localized to cytosolic sub compartments in retinal ganglion cell (RGC) somata and proximal dendrites in addition to GABAergic, glycinergic, and cholinergic amacrine cells in the inner plexiform layer (IPL) and displaced amacrine cells. Human, but not mouse retinas showed modest immunolabeling of bipolar cells. ATXN2 immunofluorescence was prominent in the trabecular meshwork and pigmented and nonpigmented cells of the ciliary body, with analyses of primary human trabecular meshwork cells confirming the finding. The expression of ATXN2 in key POAG-relevant ocular tissues supports the potential role in autophagy and stress granule formation in response to ocular hypertension.

Hagfish to Illuminate the Developmental and Evolutionary Origins of the Vertebrate Retina

The vertebrate eye is a vital sensory organ that has long fascinated scientists, but the details of how this organ evolved are still unclear. The vertebrate eye is distinct from the simple photoreceptive organs of other non-vertebrate chordates and there are no clear transitional forms of the eye in the fossil record. To investigate the evolution of the eye we can examine the eyes of the most ancient extant vertebrates, the hagfish and lamprey. These jawless vertebrates are in an ideal phylogenetic position to study the origin of the vertebrate eye but data on eye/retina development in these organisms is limited. New genomic and gene expression data from hagfish and lamprey suggest they have many of the same genes for eye development and retinal neurogenesis as jawed vertebrates, but functional work to determine if these genes operate in retinogenesis similarly to other vertebrates is missing. In addition, hagfish express a marker of proliferative retinal cells (Pax6) near the margin of the retina, and adult retinal growth is apparent in some species. This finding of eye growth late into hagfish ontogeny is unexpected given the degenerate eye phenotype. Further studies dissecting retinal neurogenesis in jawless vertebrates would allow for comparison of the mechanisms of retinal development between cyclostome and gnathostome eyes and provide insight into the evolutionary origins of the vertebrate eye.

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.

Photoreceptor cKO of OTX2 Enhances OTX2 Intercellular Transfer in the Retina and Causes Photophobia

In the mature mouse retina, Otx2 is expressed in both retinal pigmented epithelium (RPE) and photoreceptor (PR) cells, and Otx2 knock-out (KO) in the RPE alone results in PR degeneration. To study the cell-autonomous function of OTX2 in PRs, we performed PR-specific Otx2 KO (cKO) in adults. As expected, the protein disappears completely from PR nuclei but is still observed in PR inner and outer segments while its level concomitantly decreases in the RPE, suggesting a transfer of OTX2 from RPE to PRs in response to Otx2 ablation in PRs. The ability of OTX2 to transfer from RPE to PRs was verified by viral expression of tagged-OTX2 in the RPE. Transferred OTX2 distributed across the PR cytoplasm, suggesting functions distinct from nuclear transcription regulation. PR-specific Otx2 cKO did not alter the structure of the retina but impaired the translocation of PR arrestin-1 on illumination changes, making mice photophobic. RNA-seq analyses following Otx2 KO revealed downregulation of genes involved in the cytoskeleton that might account for the arrestin-1 translocation defect, and of genes involved in extracellular matrix (ECM) and signaling factors that may participate in the enhanced transfer of OTX2. Interestingly, several RPE-specific OTX2 target genes involved in melanogenesis were downregulated, lending weight to a decrease of OTX2 levels in the RPE following PR-specific Otx2 cKO. Our study reveals a new role of endogenous OTX2 in PR light adaptation and demonstrates the existence of OTX2 transfer from RPE to PR cells, which is increased on PR-specific Otx2 ablation and might participate in PR neuroprotection.

Downregulation of CRX, a Group 3-specific oncogenic transcription factor, inhibits TGF-β/activin signaling in medulloblastoma cells

Medulloblastoma, the most common malignant brain tumor in children, consists of four molecular subgroups WNT, SHH, Group 3, and Group 4. Group 3 has the worst survival rate among the four subgroups and is characterized by the expression of retina-specific genes. CRX, the master regulator of the photoreceptor differentiation, is aberrantly expressed in Group 3 medulloblastomas. CRX expression increased the proliferation, anchorage-independent growth, invasion potential, and tumorigenicity of medulloblastoma cells indicating the oncogenic role of CRX in medulloblastoma pathogenesis. CRX knockdown resulted in the downregulation of expression of several retina-specific genes like IMPG2, PDC, RCVRN. and Group 3 specific genes like GABRA5, MYC, PROM1. Thus, CRX plays a major role not only in the expression of retina-specific genes but also in defining Group 3 identity. Increased expression of several pro-apoptotic genes upon CRX knockdown suggests that CRX could protect Group 3 medulloblastoma cells from cell death. Several negative regulators of the TGF-β signaling pathway like SMAD7, PMEPA1, KLF2 were upregulated upon the CRX knockdown. Western blot analysis showed a decrease in the levels of (Phospho)-SMAD2, total levels of SMAD2, SMAD4, and an increase in the levels of SMAD7 indicating inhibition of the TGF-β signaling pathway upon CRX knockdown. Copy number variations in several genes involved in the TGF-β signaling pathway occur in a subset of Group 3 tumors. Autocrine TGF-β/activin signaling has recently been reported to be active in a subset of Group 3 medulloblastomas. CRX knockdown resulting in the inhibition of the TGF-β/activin signaling pathway demonstrates an interaction between the two Group 3 specific oncogenic pathways and suggests simultaneous targeting of both CRX and TGF-β signaling as a possible therapeutic strategy.

The role of homeobox gene-encoded transcription factors in regulation of phototransduction: Implementing the primary pinealocyte culture as a photoreceptor model

Homeobox genes encode transcription factors controlling development; however, a number of homeobox genes are expressed postnatally specifically in melatonin-producing pinealocytes of the pineal gland and photoreceptors of the retina along with transcripts devoted to melatonin synthesis and phototransduction. Homeobox genes regulate melatonin synthesis in pinealocytes, but some homeobox genes also seem to be involved in regulation of retinal phototransduction. Due to the lack of photoreceptor models, we here introduce the rat pinealocyte culture as an in vitro model for studying retinal phototransduction. Systematic qPCR analyses were performed on the rat retina and pineal gland in 24 hour in vivo series and on primary cultures of rat pinealocytes: All homeobox genes and melatonin synthesis components, as well as nine out of ten phototransduction genes, were readily detectable in all three experimental settings, confirming molecular similarity between cultured pinealocytes and in vivo retinal tissue. 24 hours circadian expression was mostly confined to transcripts in the pineal gland, including a novel rhythm in arrestin (Sag). Individual knockdown of the homeobox genes orthodenticle homeobox 2 (Otx2), cone-rod homeobox (Crx) and LIM homeobox 4 (Lhx4) in pinealocyte culture using siRNA resulted in specific downregulation of transcripts representing all levels of phototransduction; thus, all phototransduction genes studied in culture were affected by one or several siRNA treatments. Histological colocalization of homeobox and phototransduction transcripts in the rat retinal photoreceptor was confirmed by RNAscope in situ hybridization, thus suggesting that homeobox gene-encoded transcription factors control postnatal expression of phototransduction genes in the retinal photoreceptor.

Features of Retinal Neurogenesis as a Key Factor of Age-Related Neurodegeneration: Myth or Reality?

Age-related macular degeneration (AMD) is a complex multifactorial neurodegenerative disease that constitutes the most common cause of irreversible blindness in the elderly in the developed countries. Incomplete knowledge about its pathogenesis prevents the search for effective methods of prevention and treatment of AMD, primarily of its "dry" type which is by far the most common (90% of all AMD cases). In the recent years, AMD has become "younger": late stages of the disease are now detected in relatively young people. It is known that AMD pathogenesis-according to the age-related structural and functional changes in the retina-is linked with inflammation, hypoxia, oxidative stress, mitochondrial dysfunction, and an impairment of neurotrophic support, but the mechanisms that trigger the conversion of normal age-related changes to the pathological process as well as the reason for early AMD development remain unclear. In the adult mammalian retina, de novo neurogenesis is very limited. Therefore, the structural and functional features that arise during its maturation and formation can exert long-term effects on further ontogenesis of this tissue. The aim of this review was to discuss possible contributions of the changes/disturbances in retinal neurogenesis to the early development of AMD.

Initiation of Otx2 expression in the developing mouse retina requires a unique enhancer and either Ascl1 or Neurog2 activity

During retinal development, a large subset of progenitors upregulates the transcription factor Otx2, which is required for photoreceptor and bipolar cell formation. How these retinal progenitor cells initially activate Otx2 expression is unclear. To address this, we investigated the cis-regulatory network that controls Otx2 expression in mice. We identified a minimal enhancer element, DHS-4D, that drove expression in newly formed OTX2+ cells. CRISPR/Cas9-mediated deletion of DHS-4D reduced OTX2 expression, but this effect was diminished in postnatal development. Systematic mutagenesis of the enhancer revealed that three basic helix-loop-helix (bHLH) transcription factor-binding sites were required for its activity. Single cell RNA-sequencing of nascent Otx2+ cells identified the bHLH factors Ascl1 and Neurog2 as candidate regulators. CRISPR/Cas9 targeting of these factors showed that only the simultaneous loss of Ascl1 and Neurog2 prevented OTX2 expression. Our findings suggest that Ascl1 and Neurog2 act either redundantly or in a compensatory fashion to activate the DHS-4D enhancer and Otx2 expression. We observed redundancy or compensation at both the transcriptional and enhancer utilization levels, suggesting that the mechanisms governing Otx2 regulation in the retina are flexible and robust.

Single Cell Transcriptomic Analyses Reveal the Impact of bHLH Factors on Human Retinal Organoid Development

The developing retina expresses multiple bHLH transcription factors. Their precise functions and interactions in uncommitted retinal progenitors remain to be fully elucidated. Here, we investigate the roles of bHLH factors ATOH7 and Neurog2 in human ES cell-derived retinal organoids. Single cell transcriptome analyses identify three states of proliferating retinal progenitors: pre-neurogenic, neurogenic, and cell cycle-exiting progenitors. Each shows different expression profile of bHLH factors. The cell cycle-exiting progenitors feed into a postmitotic heterozygous neuroblast pool that gives rise to early born neuronal lineages. Elevating ATOH7 or Neurog2 expression accelerates the transition from the pre-neurogenic to the neurogenic state, and expands the exiting progenitor and neuroblast populations. In addition, ATOH7 and Neurog2 significantly, yet differentially, enhance retinal ganglion cell and cone photoreceptor production. Moreover, single cell transcriptome analyses reveal that ATOH7 and Neurog2 each assert positive autoregulation, and both suppress key bHLH factors associated with the pre-neurogenic and states and elevate bHLH factors expressed by exiting progenitors and differentiating neuroblasts. This study thus provides novel insight regarding how ATOH7 and Neurog2 impact human retinal progenitor behaviors and neuroblast fate choices.

Functional analysis of Samd11, a retinal photoreceptor PRC1 component, in establishing rod photoreceptor identity

Establishing correct neuronal cell identity is essential to build intricate neural tissue architecture and acquire precise neural function during vertebrate development. While it is known that transcription factors play important roles in retinal cell differentiation, the contribution of epigenetic factors to establishing cell identity during retinal development remains unclear. We previously reported that Samd7, a rod photoreceptor cell-specific sterile alpha motif (SAM) domain protein, functions as a Polycomb repressive complex 1 component (PRC1) that is essential for establishing rod identity. In the current study, we analyzed a functional role of Samd11, another photoreceptor-enriched SAM-domain protein, in photoreceptor differentiation and maturation. We observed that Samd11 interacts with Phc2 and Samd7, suggesting that Samd11 is a component of PRC1 in photoreceptor cells. We generated Samd11-null allele and established Samd7/11 double knock-out (DKO) mouse. The Samd7/11 DKO retina exhibits shortened photoreceptor outer segments by electron microscopy analysis. Microarray analysis revealed that Samd7/11 DKO up-regulated more retinal genes than Samd7^-/- alone, partial functional redundancy of Samd7 and Samd11. Taken together, the current results suggest that Samd7 and Samd11 are PRC1 components and that Samd7 is the major regulator while Samd11 is an accessory factor used for the establishment of precise rod photoreceptor identity.

Influence of Aging on the Retina and Visual Motion Processing for Optokinetic Responses in Mice

The decline in visual function due to normal aging impacts various aspects of our daily lives. Previous reports suggest that the aging retina exhibits mislocalization of photoreceptor terminals and reduced amplitudes of scotopic and photopic electroretinogram (ERG) responses in mice. These abnormalities are thought to contribute to age-related visual impairment; however, the extent to which visual function is impaired by aging at the organismal level is unclear. In the present study, we focus on the age-related changes of the optokinetic responses (OKRs) in visual processing. Moreover, we investigated the initial and late phases of the OKRs in young adult (2-3 months old) and aging mice (21-24 months old). The initial phase was evaluated by measuring the open-loop eye velocity of OKRs using sinusoidal grating patterns of various spatial frequencies (SFs) and moving at various temporal frequencies (TFs) for 0.5 s. The aging mice exhibited initial OKRs with a spatiotemporal frequency tuning that was slightly different from those in young adult mice. The late-phase OKRs were investigated by measuring the slow-phase velocity of the optokinetic nystagmus evoked by sinusoidal gratings of various spatiotemporal frequencies moving for 30 s. We found that optimal SF and TF in the normal aging mice are both reduced compared with those in young adult mice. In addition, we measured the OKRs of 4.1G-null (4.1G ^-/-) mice, in which mislocalization of photoreceptor terminals is observed even at the young adult stage. We found that the late phase OKR was significantly impaired in 4.1G ^{- / -} mice, which exhibit significantly reduced SF and TF compared with control mice. These OKR abnormalities observed in 4.1G ^{- / -} mice resemble the abnormalities found in normal aging mice. This finding suggests that these mice can be useful mouse models for studying the aging of the retinal tissue and declining visual function. Taken together, the current study demonstrates that normal aging deteriorates to visual motion processing for both the initial and late phases of OKRs. Moreover, it implies that the abnormalities of the visual function in the normal aging mice are at least partly due to mislocalization of photoreceptor synapses.

The potential role of Arhgef33 RhoGEF in foveal development in the zebra finch retina

The fovea is a pit formed in the center of the retina that enables high-acuity vision in certain vertebrate species. While formation of the fovea fascinates many researchers, the molecular mechanisms underlying foveal development are poorly understood. In the current study, we histologically investigated foveal development in zebra finch (Taeniopygia guttata) and found that foveal pit formation begins just before post-hatch day 14 (P14). We next performed RNA-seq analysis to compare gene expression profiles between the central (foveal and parafoveal) and peripheral retina in zebra finch at P14. We found that the Arhgef33 expression is enriched in the middle layer of the inner nuclear layer at the parafovea, suggesting that Arhgef33 is dominantly expressed in Müller glial cells in the developing parafovea. We then performed a pull-down assay using Rhotekin-RBD and observed GEF activity of Arhgef33 against RhoA. We found that overexpression of Arhgef33 in HEK293 cells induces cell contraction and that Arhgef33 expression inhibits neurite extension in Neuro 2A cells, which is partially recovered by a Rho-kinase (ROCK) inhibitor. Taken together, we used zebra finch as a model animal to investigate foveal development and identified Arhgef33 as a candidate protein possibly involved in foveal development through modulating RhoA activity.

Mechanistic insights into the evolution of the differential expression of tandemly arrayed cone opsin genes in zebrafish

The genome of many organisms contains several loci consisting of duplicated genes that are arrayed in tandem. The daughter genes produced by duplication typically exhibit differential expression patterns with each other or otherwise experience pseudogenization. Remarkably, opsin genes in fish are preserved after many duplications in different lineages. This fact indicates that fish opsin genes are characterized by a regulatory mechanism that could intrinsically facilitate the differentiation of the expression patterns. However, little is known about the mechanisms that underlie the differential expression patterns or how they were established during evolution. The loci of green (RH2)- and red (LWS)-sensitive cone opsin genes in zebrafish have been used as model systems to study the differential regulation of tandemly arrayed opsin genes. Over a decade of studies have uncovered several mechanistic features that might have assisted the differentiation and preservation of duplicated genes. Furthermore, recent progress in the understanding of the transcriptional process in general has added essential insights. In this article, the current understanding of the transcriptional regulation of differentially expressed tandemly arrayed cone opsin genes in zebrafish is summarized and a possible evolutionary scenario that could achieve this differentiation is discussed.

Simultaneous deletion of Prdm1 and Vsx2 enhancers in the retina alters photoreceptor and bipolar cell fate specification, yet differs from deleting both genes

The transcription factor OTX2 is required for photoreceptor and bipolar cell formation in the retina. It directly activates the transcription factors Prdm1 and Vsx2 through cell type-specific enhancers. PRDM1 and VSX2 work in opposition, such that PRDM1 promotes photoreceptor fate and VSX2 bipolar cell fate. To determine how OTX2+ cell fates are regulated in mice, we deleted Prdm1 and Vsx2 or their cell type-specific enhancers simultaneously using a CRISPR/Cas9 in vivo retina electroporation strategy. Double gene or enhancer targeting effectively removed PRDM1 and VSX2 protein expression. However, double enhancer targeting favored bipolar fate outcomes, whereas double gene targeting favored photoreceptor fate. Both conditions generated excess amacrine cells. Combined, these fate changes suggest that photoreceptors are a default fate outcome in OTX2+ cells and that VSX2 must be present in a narrow temporal window to drive bipolar cell formation. Prdm1 and Vsx2 also appear to redundantly restrict the competence of OTX2+ cells, preventing amacrine cell formation. By taking a combinatorial deletion approach of both coding sequences and enhancers, our work provides new insights into the complex regulatory mechanisms that control cell fate choice.

Prdm1 overexpression causes a photoreceptor fate-shift in nascent, but not mature, bipolar cells

The transcription factors Prdm1 (Blimp1) and Vsx2 (Chx10) work downstream of Otx2 to regulate photoreceptor and bipolar cell fates in the developing retina. Mice that lack Vsx2 fail to form bipolar cells while Prdm1 mutants form excess bipolars at the direct expense of photoreceptors. Excess bipolars in Prdm1 mutants appear to derive from rods, suggesting that photoreceptor fate remains mutable for some time after cells become specified. Here we tested whether bipolar cell fate is also plastic during development. To do this, we created a system to conditionally misexpress Prdm1 at different stages of bipolar cell development. We found that Prdm1 blocks bipolar cell formation if expressed before the fate choice decision occurred. When we misexpressed Prdm1 just after the decision to become a bipolar cell was made, some cells were reprogrammed into photoreceptors. In contrast, Prdm1 misexpression in mature bipolar cells did not affect cell fate. We also provide evidence that sustained misexpression of Prdm1 was selectively toxic to photoreceptors. Our data show that bipolar fate is malleable, but only for a short temporal window following fate specification. Prdm1 and Vsx2 act by stabilizing photoreceptor and bipolar fates in developing OTX2+ cells of the retina.