Dissecting the transcriptome landscape of the human fetal neural retina and retinal pigment epithelium by single-cell RNA-seq analysis

Single-cell RNA sequencing analysis of the retina under acute high intraocular pressure

JOURNAL/nrgr/04.03/01300535-202419110-00032/figure1/v/2024-03-08T184507Z/r/image-tiff High intraocular pressure causes retinal ganglion cell injury in primary and secondary glaucoma diseases, yet the molecular landscape characteristics of retinal cells under high intraocular pressure remain unknown. Rat models of acute hypertension ocular pressure were established by injection of cross-linked hyaluronic acid hydrogel (Healaflow®). Single-cell RNA sequencing was then used to describe the cellular composition and molecular profile of the retina following high intraocular pressure. Our results identified a total of 12 cell types, namely retinal pigment epithelial cells, rod-photoreceptor cells, bipolar cells, Müller cells, microglia, cone-photoreceptor cells, retinal ganglion cells, endothelial cells, retinal progenitor cells, oligodendrocytes, pericytes, and fibroblasts. The single-cell RNA sequencing analysis of the retina under acute high intraocular pressure revealed obvious changes in the proportions of various retinal cells, with ganglion cells decreased by 23%. Hematoxylin and eosin staining and TUNEL staining confirmed the damage to retinal ganglion cells under high intraocular pressure. We extracted data from retinal ganglion cells and analyzed the retinal ganglion cell cluster with the most distinct expression. We found upregulation of the B3gat2 gene, which is associated with neuronal migration and adhesion, and downregulation of the Tsc22d gene, which participates in inhibition of inflammation. This study is the first to reveal molecular changes and intercellular interactions in the retina under high intraocular pressure. These data contribute to understanding of the molecular mechanism of retinal injury induced by high intraocular pressure and will benefit the development of novel therapies.

Temporal Transcriptomic Profiling of the Developing Xenopus laevis Eye

Retinal progenitor cells (RPCs) are a multipotent and highly proliferative population that give rise to all retinal cell types during organogenesis. Defining their molecular signature is a key step towards identifying suitable approaches to treat visual impairments. Here, we performed RNA sequencing of whole eyes from Xenopus at three embryonic stages and used differential expression analysis to define the transcriptomic profiles of optic tissues containing proliferating and differentiating RPCs during retinogenesis. Gene Ontology and KEGG pathway analyses showed that genes associated with developmental pathways (including Wnt and Hedgehog signaling) were upregulated during the period of active RPC proliferation in early retinal development (Nieuwkoop Faber st. 24 and 27). Developing eyes had dynamic expression profiles and shifted to enrichment for metabolic processes and phototransduction during RPC progeny specification and differentiation (st. 35). Furthermore, conserved adult eye regeneration genes were also expressed during early retinal development, including sox2, pax6, nrl, and Notch signaling components. The eye transcriptomic profiles presented here span RPC proliferation to retinogenesis and include regrowth-competent stages. Thus, our dataset provides a rich resource to uncover molecular regulators of RPC activity and will allow future studies to address regulators of RPC proliferation during eye repair and regrowth.

Temporal Transcriptomic Profiling of the Developing Xenopus laevis Eye

Retinal progenitor cells (RPCs) are a multipotent and highly proliferative population that give rise to all retinal cell types during organogenesis. Defining their molecular signature is a key step towards identifying suitable approaches to treat visual impairments. Here, we performed RNA-sequencing of whole eyes from Xenopus at three embryonic stages and used differential expression analysis to define the transcriptomic profiles of optic tissues containing proliferating and differentiating RPCs during retinogenesis. Gene Ontology and KEGG pathway analyses showed that genes associated with developmental pathways (including Wnt and Hedgehog signaling) were upregulated during the period of active RPC proliferation in early retinal development (Nieuwkoop Faber st. 24 and 27). Developing eyes had dynamic expression profiles and shifted to enrichment for metabolic processes and phototransduction during RPC progeny specification and differentiation (st. 35). Furthermore, conserved adult eye regeneration genes were also expressed during early retinal development including sox2, pax6, nrl, and Notch signaling components. The eye transcriptomic profiles presented here span RPC proliferation to retinogenesis and included regrowth-competent stages. Thus, our dataset provides a rich resource to uncover molecular regulators of RPC activity and will allow future studies to address regulators of RPC proliferation during eye repair and regrowth.

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.

ID factors regulate the ability of Müller glia to become proliferating neurogenic progenitor-like cells

The purpose of this study was to investigate how ID factors regulate the ability of Müller glia (MG) to reprogram into proliferating MG-derived progenitor cells (MGPCs) in the chick retina. We found that ID1 is transiently expressed by maturing MG (mMG), whereas ID4 is maintained in mMG in embryonic retinas. In mature retinas, ID4 was prominently expressed by resting MG, but following retinal damage ID4 was rapidly upregulated and then downregulated in MGPCs. By contrast, ID1, ID2, and ID3 were low in resting MG and then upregulated in MGPCs. Inhibition of ID factors following retinal damage decreased numbers of proliferating MGPCs. Inhibition of IDs, after MGPC proliferation, significantly increased numbers of progeny that differentiated as neurons. In damaged or undamaged retinas inhibition of IDs increased levels of p21Cip1 in MG. In response to damage or insulin+FGF2 levels of CDKN1A message and p21Cip1 protein were decreased, absent in proliferating MGPCs, and elevated in MG returning to a resting phenotype. Inhibition of notch- or gp130/Jak/Stat-signaling in damaged retinas increased levels of ID4 but not p21Cip1 in MG. Although ID4 is the predominant isoform expressed by MG in the chick retina, id1 and id2a are predominantly expressed by resting MG and downregulated in activated MG and MGPCs in zebrafish retinas. We conclude that ID factors have a significant impact on regulating the responses of MG to retinal damage, controlling the ability of MG to proliferate by regulating levels of p21Cip1, and suppressing the neurogenic potential of MGPCs.

Human CRB1 and CRB2 form homo- and heteromeric protein complexes in the retina

Crumbs homolog 1 (CRB1) is one of the key genes linked to retinitis pigmentosa and Leber congenital amaurosis, which are characterized by a high clinical heterogeneity. The Crumbs family member CRB2 has a similar protein structure to CRB1, and in zebrafish, Crb2 has been shown to interact through the extracellular domain. Here, we show that CRB1 and CRB2 co-localize in the human retina and human iPSC-derived retinal organoids. In retina-specific pull-downs, CRB1 was enriched in CRB2 samples, supporting a CRB1-CRB2 interaction. Furthermore, novel interactors of the crumbs complex were identified, representing a retina-derived protein interaction network. Using co-immunoprecipitation, we further demonstrate that human canonical CRB1 interacts with CRB1 and CRB2, but not with CRB3, which lacks an extracellular domain. Next, we explored how missense mutations in the extracellular domain affect CRB1-CRB2 interactions. We observed no or a mild loss of CRB1-CRB2 interaction, when interrogating various CRB1 or CRB2 missense mutants in vitro. Taken together, our results show a stable interaction of human canonical CRB2 and CRB1 in the retina.

Agarose disk electroporation method for ex vivo retinal tissue cultured at the air-liquid interface reveals electrical stimulus-induced cell cycle reentry in retinal cells

It is advantageous to culture the ex vivo murine retina along with many other tissue types at the air-liquid interface. However, gene delivery to these cultures can be challenging. Electroporation is a fast and robust method of gene delivery, but typically requires submergence in a liquid buffer to allow electric current flow. We have developed a submergence-free electroporation technique using an agarose disk that allows for efficient gene delivery to the ex vivo murine retina. This method advances our ability to use ex vivo retinal tissue for genetic studies and can easily be adapted for any tissue cultured at an air-liquid interface. We found an increased ability to transfected Muller glia at 14 days ex vivo and an increase in BrdU incorporation in Muller glia following electrical stimulation. Use of this method has revealed valuable insights on the state of ex vivo retinal tissues and the effects of electrical stimulation on retinal cells.

Single-cell analyses reveal transient retinal progenitor cells in the ciliary margin of developing human retina

The emergence of retinal progenitor cells and differentiation to various retinal cell types represent fundamental processes during retinal development. Herein, we provide a comprehensive single cell characterisation of transcriptional and chromatin accessibility changes that underline retinal progenitor cell specification and differentiation over the course of human retinal development up to midgestation. Our lineage trajectory data demonstrate the presence of early retinal progenitors, which transit to late, and further to transient neurogenic progenitors, that give rise to all the retinal neurons. Combining single cell RNA-Seq with spatial transcriptomics of early eye samples, we demonstrate the transient presence of early retinal progenitors in the ciliary margin zone with decreasing occurrence from 8 post-conception week of human development. In retinal progenitor cells, we identified a significant enrichment for transcriptional enhanced associate domain transcription factor binding motifs, which when inhibited led to loss of cycling progenitors and retinal identity in pluripotent stem cell derived organoids.

Clinically compliant cryopreservation of differentiated retinal pigment epithelial cells

BACKGROUND AIMS

Age-related macular degeneration (AMD) is the most common cause of blindness in elderly patients within developed countries, affecting more than 190 million worldwide. In AMD, the retinal pigment epithelial (RPE) cell layer progressively degenerates, resulting in subsequent loss of photoreceptors and ultimately vision. There is currently no cure for AMD, but therapeutic strategies targeting the complement system are being developed to slow the progression of the disease.

METHODS

Replacement therapy with pluripotent stem cell-derived (hPSC) RPEs is an alternative treatment strategy. A cell therapy product must be produced in accordance with Good Manufacturing Practices at a sufficient scale to facilitate extensive pre-clinical and clinical testing. Cryopreservation of the final cell product is therefore highly beneficial, as the manufacturing, pre-clinical and clinical testing can be separated in time and location.

RESULTS

We found that mature hPSC-RPE cells do not survive conventional cryopreservation techniques. However, replating the cells 2-5 days before cryopreservation facilitates freezing. The replated and cryopreserved hPSC-RPE cells maintained their identity, purity and functionality as characteristic RPEs, shown by cobblestone morphology, pigmentation, transcriptional profile, RPE markers, transepithelial resistance and pigment epithelium-derived factor secretion. Finally, we showed that the optimal replating time window can be tracked noninvasively by following the change in cobblestone morphology.

CONCLUSIONS

The possibility of cryopreserving the hPSC-RPE product has been instrumental in our efforts in manufacturing and performing pre-clinical testing with the aim for clinical translation.

scRNA-Seq: First Atlas and Cellular Landscape of Lacrimal Sac: Implications in Primary Acquired Nasolacrimal Duct Obstruction Pathogenesis

Purpose

The aim of this study was to describe the transcriptional changes of individual cellular components in the lacrimal sac in patients with primary acquired nasolacrimal duct obstruction (PANDO) and attempt to construct the first lacrimal sac cellular atlas to elucidate the potential mechanisms that may drive the disease pathogenesis.

Methods

Lacrimal sac samples were obtained intra-operatively during the endoscopic dacryocystorhinostomy (EnDCR) procedure from five patients. Single-cell RNA sequencing was performed to analyze each individual cell population including epithelial and immune cells during the early inflammatory and late inflammatory phases of the disease.

Results

Eleven cell types were identified among 25,791 cells. T cells and B cells were the cell populations with the greatest variation in cell numbers between the two phases and were involved in immune response and epithelium migration-related pathways. The present study showed that epithelial cells highly expressed the genes of senescence-associated secretory phenotype (SASP) and were involved in influencing the inflammation, neutrophil chemotaxis, and migration during the late inflammatory stage. Enhanced activity of CXCLs-CXCRs between the epithelial cells and neutrophils was noted by the cell-cell communication analysis and is suspected to play a role in inflammation by recruiting more neutrophils.

Conclusions

The study presents a comprehensive single-cell landscape of the lacrimal sac cells in different phases of PANDO. The contribution of T cells, B cells, and epithelial cells to the inflammatory response, and construction of the intercellular signaling networks between the cells within the lacrimal sac has further enhanced the present understanding of the PANDO pathogenesis.

Recent Achievements in the Heterogeneity of Mammalian and Human Retinal Pigment Epithelium: In Search of a Stem Cell

Retinal pigment epithelium (RPE) cells are important fundamentally for the development and function of the retina. In this regard, the study of the morphological and molecular properties of RPE cells, as well as their regenerative capabilities, is of particular importance for biomedicine. However, these studies are complicated by the fact that, despite the external morphological similarity of RPE cells, the RPE is a population of heterogeneous cells, the molecular genetic properties of which have begun to be revealed by sequencing methods only in recent years. This review carries out an analysis of the data from morphological and molecular genetic studies of the heterogeneity of RPE cells in mammals and humans, which reveals the individual differences in the subpopulations of RPE cells and the possible specificity of their functions. Particular attention is paid to discussing the properties of "stemness," proliferation, and plasticity in the RPE, which may be useful for uncovering the mechanisms of retinal diseases associated with pathologies of the RPE and finding new ways of treating them.

Integrating mechanical cues with engineered platforms to explore cardiopulmonary development and disease

Mechanical forces provide critical biological signals to cells during healthy and aberrant organ development as well as during disease processes in adults. Within the cardiopulmonary system, mechanical forces, such as shear, compressive, and tensile forces, act across various length scales, and dysregulated forces are often a leading cause of disease initiation and progression such as in bronchopulmonary dysplasia and cardiomyopathies. Engineered in vitro models have supported studies of mechanical forces in a number of tissue and disease-specific contexts, thus enabling new mechanistic insights into cardiopulmonary development and disease. This review first provides fundamental examples where mechanical forces operate at multiple length scales to ensure precise lung and heart function. Next, we survey recent engineering platforms and tools that have provided new means to probe and modulate mechanical forces across in vitro and in vivo settings. Finally, the potential for interdisciplinary collaborations to inform novel therapeutic approaches for a number of cardiopulmonary diseases are discussed.

ASCL1 induces neurogenesis in human Müller glia

In mammals, loss of retinal cells due to disease or trauma is an irreversible process that can lead to blindness. Interestingly, regeneration of retinal neurons is a well established process in some non-mammalian vertebrates and is driven by the Müller glia (MG), which are able to re-enter the cell cycle and reprogram into neurogenic progenitors upon retinal injury or disease. Progress has been made to restore this mechanism in mammals to promote retinal regeneration: MG can be stimulated to generate new neurons in vivo in the adult mouse retina after the over-expression of the pro-neural transcription factor Ascl1. In this study, we applied the same strategy to reprogram human MG derived from fetal retina and retinal organoids into neurons. Combining single cell RNA sequencing, single cell ATAC sequencing, immunofluorescence, and electrophysiology we demonstrate that human MG can be reprogrammed into neurogenic cells in vitro.

Rare intercellular material transfer as a confound to interpreting inner retinal neuronal transplantation following internal limiting membrane disruption

Intercellular cytoplasmic material transfer (MT) occurs between transplanted and developing photoreceptors and ambiguates cell origin identification in developmental, transdifferentiation, and transplantation experiments. Whether MT is a photoreceptor-specific phenomenon is unclear. Retinal ganglion cell (RGC) replacement, through transdifferentiation or transplantation, holds potential for restoring vision in optic neuropathies. During careful assessment for MT following human stem cell-derived RGC transplantation into mice, we identified RGC xenografts occasionally giving rise to labeling of donor-derived cytoplasmic, nuclear, and mitochondrial proteins within recipient Müller glia. Critically, nuclear organization is distinct between human and murine retinal neurons, which enables unequivocal discrimination of donor from host cells. MT was greatly facilitated by internal limiting membrane disruption, which also augments retinal engraftment following transplantation. Our findings demonstrate that retinal MT is not unique to photoreceptors and challenge the isolated use of species-specific immunofluorescent markers for xenotransplant identification. Assessment for MT is critical when analyzing neuronal replacement interventions.

Dynamic human retinal pigment epithelium (RPE) and choroid architecture based on single-cell transcriptomic landscape analysis

The retinal pigment epithelium (RPE) and choroid are located behind the human retina and have multiple functions in the human visual system. Knowledge of the RPE and choroid cells and their gene expression profiles are fundamental for understanding retinal disease mechanisms and therapeutic strategies. Here, we sequenced the RNA of about 0.3 million single cells from human RPE and choroids across two regions and seven ages, revealing regional and age differences within the human RPE and choroid. Cell-cell interactions highlight the broad connectivity networks between the RPE and different choroid cell types. Moreover, the transcription factors and their target genes change during aging. The coding of somatic variations increases during aging in the human RPE and choroid at the single-cell level. Moreover, we identified ELN as a candidate for improving RPE degeneration and choroidal structure during aging. The mapping of the molecular architecture of the human RPE and choroid improves our understanding of the human vision support system and offers potential insights into the intervention targets for retinal diseases.

Identification of the immune-associated characteristics and predictive biomarkers of keratoconus based on single-cell RNA-sequencing and bulk RNA-sequencing

Background

Whether keratoconus (KC) is an inflammatory disease is currently debated. Hence, we aimed to investigate the immune-related features of KC based on single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing (bulk RNA-seq) data.

Methods

scRNA-seq data were obtained from the Genome Sequence Archive (GSA), bulk RNA-seq data were obtained from the Gene Expression Omnibus (GEO), and immune-associated genes(IAGs) were obtained from the ImmPort database. Cell clusters of KC were annotated, and different cell clusters were then selected. The IAG score of each cell was calculated using the AUCell package. Three bulk RNA-seq datasets were merged and used to identify the differentially expressed genes (DEGs), biological functions, and immune characteristics. Weighted gene coexpression network analysis (WGCNA) was used to select the IAG score-related hub genes. Based on scRNA-seq and bulk RNA-seq analyses, three machine learning algorithms, including random forest (RF), support vector machine (SVM), and least absolute shrinkage and selection operator (LASSO) regression analysis, were used to identify potential prognostic markers for KC. A predictive nomogram was developed based on prognostic markers.

Results

Six cell clusters were identified in KC, and decreased corneal stromal cell-5 (CSC-5) and increased CSC-6 were found in KC. CSC and immune cell clusters had the highest IAG scores. The bulk RNA-seq analysis identified 1362 DEGs (553 upregulated and 809 downregulated) in KC. We found different immune cell populations and differentially expressed cytokines in KC. More than three key IAG score-related modules and 367 genes were identified. By integrating the scRNA-seq and bulk RNA-seq analyses, 250 IAGs were selected and then incorporated into three machine learning models, and 10 IAGs (CEP112, FYN, IFITM1, IGFBP5, LPIN2, MAP1B, RNASE1, RUNX3, SMIM10, and SRGN) were identified as potential prognostic genes that were significantly associated with cytokine and matrix metalloproteinase(MMP)1-14 expression. Finally, a predictive nomogram was constructed and validated.

Conclusion

Taken together, our results identified CSCs and immune cell clusters that may play a key role during KC progression by regulating immunological features and maintaining cell stability.

Transcriptomic analysis of the ocular posterior segment completes a cell atlas of the human eye

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here, we used high-throughput single-nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in six extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases-for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ~151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the data presented here complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as the neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

The efficient induction of human retinal ganglion-like cells provides a platform for studying optic neuropathies

Retinal ganglion cells (RGCs) are essential for vision perception. In glaucoma and other optic neuropathies, RGCs and their optic axons undergo degenerative change and cell death; this can result in irreversible vision loss. Here we developed a rapid protocol for directly inducing RGC differentiation from human induced pluripotent stem cells (hiPSCs) by the overexpression of ATOH7, BRN3B, and SOX4. The hiPSC-derived RGC-like cells (iRGCs) show robust expression of various RGC-specific markers by whole transcriptome profiling. A functional assessment was also carried out and this demonstrated that these iRGCs display stimulus-induced neuronal activity, as well as spontaneous neuronal activity. Ethambutol (EMB), an effective first-line anti-tuberculosis agent, is known to cause serious visual impairment and irreversible vision loss due to the RGC degeneration in a significant number of treated patients. Using our iRGCs, EMB was found to induce significant dose-dependent and time-dependent increases in cell death and neurite degeneration. Western blot analysis revealed that the expression levels of p62 and LC3-II were upregulated, and further investigations revealed that EMB caused a blockade of lysosome-autophagosome fusion; this indicates that impairment of autophagic flux is one of the adverse effects of that EMB has on iRGCs. In addition, EMB was found to elevate intracellular reactive oxygen species (ROS) levels increasing apoptotic cell death. This could be partially rescued by the co-treatment with the ROS scavenger NAC. Taken together, our findings suggest that this iRGC model, which achieves both high yield and high purity, is suitable for investigating optic neuropathies, as well as being useful when searching for potential drugs for therapeutic treatment and/or disease prevention.

The specious art of single-cell genomics

Dimensionality reduction is standard practice for filtering noise and identifying relevant features in large-scale data analyses. In biology, single-cell genomics studies typically begin with reduction to 2 or 3 dimensions to produce "all-in-one" visuals of the data that are amenable to the human eye, and these are subsequently used for qualitative and quantitative exploratory analysis. However, there is little theoretical support for this practice, and we show that extreme dimension reduction, from hundreds or thousands of dimensions to 2, inevitably induces significant distortion of high-dimensional datasets. We therefore examine the practical implications of low-dimensional embedding of single-cell data and find that extensive distortions and inconsistent practices make such embeddings counter-productive for exploratory, biological analyses. In lieu of this, we discuss alternative approaches for conducting targeted embedding and feature exploration to enable hypothesis-driven biological discovery.

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.

Integrative Single-Cell Transcriptomics and Epigenomics Mapping of the Fetal Retina Developmental Dynamics

The underlying mechanisms that determine gene expression and chromatin accessibility in retinogenesis are poorly understood. Herein, single-cell RNA sequencing and single-cell assay for transposase-accessible chromatin sequencing are performed on human embryonic eye samples obtained 9-26 weeks after conception to explore the heterogeneity of retinal progenitor cells (RPCs) and neurogenic RPCs. The differentiation trajectory from RPCs to 7 major types of retinal cells are verified. Subsequently, diverse lineage-determining transcription factors are identified and their gene regulatory networks are refined at the transcriptomic and epigenomic levels. Treatment of retinospheres, with the inhibitor of RE1 silencing transcription factor, X5050, induces more neurogenesis with the regular arrangement, and a decrease in Müller glial cells. The signatures of major retinal cells and their correlation with pathogenic genes associated with multiple ocular diseases, including uveitis and age-related macular degeneration are also described. A framework for the integrated exploration of single-cell developmental dynamics of the human primary retina is provided.

Single-cell RNA sequencing in cornea research: Insights into limbal stem cells and their niche regulation

The corneal epithelium is composed of stratified squamous epithelial cells on the outer surface of the eye, which acts as a protective barrier and is critical for clear and stable vision. Its continuous renewal or wound healing depends on the proliferation and differentiation of limbal stem cells (LSCs), a cell population that resides at the limbus in a highly regulated niche. Dysfunction of LSCs or their niche can cause limbal stem cell deficiency, a disease that is manifested by failed epithelial wound healing or even blindness. Nevertheless, compared to stem cells in other tissues, little is known about the LSCs and their niche. With the advent of single-cell RNA sequencing, our understanding of LSC characteristics and their microenvironment has grown considerably. In this review, we summarized the current findings from single-cell studies in the field of cornea research and focused on important advancements driven by this technology, including the heterogeneity of the LSC population, novel LSC markers and regulation of the LSC niche, which will provide a reference for clinical issues such as corneal epithelial wound healing, ocular surface reconstruction and interventions for related diseases.

Single-cell RNA sequencing reveals transcriptional changes of human choroidal and retinal pigment epithelium cells during fetal development, in healthy adult and intermediate age-related macular degeneration

Age-related macular degeneration (AMD) is the most prevalent cause of blindness in the developed world. Vision loss in the advanced stages of the disease is caused by atrophy of retinal photoreceptors, overlying retinal pigment epithelium (RPE) and choroidal endothelial cells. The molecular events that underline the development of these cell types from in utero to adult as well as the progression to intermediate and advanced stages AMD are not yet fully understood. We performed single-cell RNA-sequencing (RNA-Seq) of human fetal and adult RPE-choroidal tissues, profiling in detail all the cell types and elucidating cell type-specific proliferation, differentiation and immunomodulation events that occur up to midgestation. Our data demonstrate that progression from the fetal to adult state is characterized by an increase in expression of genes involved in the oxidative stress response and detoxification from heavy metals, suggesting a better defence against oxidative stress in the adult RPE-choroid tissue. Single-cell comparative transcriptional analysis between a patient with intermediate AMD and an unaffected subject revealed a reduction in the number of RPE cells and melanocytes in the macular region of the AMD patient. Together these findings may suggest a macular loss of RPE cells and melanocytes in the AMD patients, but given the complex processing of tissues required for single-cell RNA-Seq that is prone to technical artefacts, these findings need to be validated by additional techniques in a larger number of AMD patients and controls.

Transcriptomic Analysis of the Ocular Posterior Segment Completes a Cell Atlas of the Human Eye

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here we used high-throughput single nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in the extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases - for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ∼151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the new data complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

A single-cell transcriptome atlas profiles early organogenesis in human embryos

The early window of human embryogenesis is largely a black box for developmental biologists. Here we probed the cellular diversity of 4-6 week human embryos when essentially all organs are just laid out. On the basis of over 180,000 single-cell transcriptomes, we generated a comprehensive atlas of 313 clusters in 18 developmental systems, which were annotated with a collection of ontology and markers from 157 publications. Together with spatial transcriptome on embryonic sections, we characterized the molecule and spatial architecture of previously unappreciated cell types. Combined with data from other vertebrates, the rich information shed light on spatial patterning of axes, systemic temporal regulation of developmental progression and potential human-specific regulation. Our study provides a compendium of early progenitor cells of human organs, which can serve as the root of lineage analysis in organogenesis.

Deep Sc-RNA sequencing decoding the molecular dynamic architecture of the human retina

The human retina serves as a light detector and signals transmission tissue. Advanced insights into retinal disease mechanisms and therapeutic strategies require a deep understanding of healthy retina molecular events. Here, we sequenced the mRNA of over 0.6 million single cells from human retinas across six regions at nine different ages. Sixty cell sub-types have been identified from the human mature retinas with unique markers. We revealed regional and age differences of gene expression profiles within the human retina. Cell-cell interaction analysis indicated a rich synaptic connection within the retinal cells. Gene expression regulon analysis revealed the specific expression of transcription factors and their regulated genes in human retina cell types. Some of the gene's expression, such as DKK3, are elevated in aged retinas. A further functional investigation suggested that over expression of DKK3 could impact mitochondrial stability. Overall, decoding the molecular dynamic architecture of the human retina improves our understanding of the vision system.

Zinc Supplementation Induced Transcriptional Changes in Primary Human Retinal Pigment Epithelium: A Single-Cell RNA Sequencing Study to Understand Age-Related Macular Degeneration

Zinc supplementation has been shown to be beneficial to slow the progression of age-related macular degeneration (AMD). However, the molecular mechanism underpinning this benefit is not well understood. This study used single-cell RNA sequencing to identify transcriptomic changes induced by zinc supplementation. Human primary retinal pigment epithelial (RPE) cells could mature for up to 19 weeks. After 1 or 18 weeks in culture, we supplemented the culture medium with 125 µM added zinc for one week. RPE cells developed high transepithelial electrical resistance, extensive, but variable pigmentation, and deposited sub-RPE material similar to the hallmark lesions of AMD. Unsupervised cluster analysis of the combined transcriptome of the cells isolated after 2, 9, and 19 weeks in culture showed considerable heterogeneity. Clustering based on 234 pre-selected RPE-specific genes divided the cells into two distinct clusters, we defined as more and less differentiated cells. The proportion of more differentiated cells increased with time in culture, but appreciable numbers of cells remained less differentiated even at 19 weeks. Pseudotemporal ordering identified 537 genes that could be implicated in the dynamics of RPE cell differentiation (FDR < 0.05). Zinc treatment resulted in the differential expression of 281 of these genes (FDR < 0.05). These genes were associated with several biological pathways with modulation of ID1/ID3 transcriptional regulation. Overall, zinc had a multitude of effects on the RPE transcriptome, including several genes involved in pigmentation, complement regulation, mineralization, and cholesterol metabolism processes associated with AMD.

Spatial distribution of metabolites in the retina and its relevance to studies of metabolic retinal disorders

INTRODUCTION

The primate retina has evolved regional specialisations for specific visual functions. The macula is specialised towards high acuity vision and is an area that contains an increased density of cone photoreceptors and signal processing neurons. Different regions in the retina display unique susceptibility to pathology, with many retinal diseases primarily affecting the macula.

OBJECTIVES

To better understand the properties of different retinal areas we studied the differential distribution of metabolites across the retina.

METHODS

We conducted an untargeted metabolomics analysis on full-thickness punches from three different regions (macula, temporal peri-macula and periphery) of healthy primate retina.

RESULTS

Nearly half of all metabolites identified showed differential abundance in at least one comparison between the three regions. Furthermore, mapping metabolomics results from macula-specific eye diseases onto our region-specific metabolite distributions revealed differential abundance defining systemic metabolic dysregulations that were region specific.

CONCLUSIONS

The unique metabolic phenotype of different retinal regions is likely due to the differential distribution of different cell types in these regions reflecting the specific metabolic requirements of each cell type. Our results may help to better understand the pathobiology of retinal diseases with region specificity.

Incorporating microglia-like cells in human induced pluripotent stem cell-derived retinal organoids

Microglia are the primary resident immune cells in the retina. They regulate neuronal survival and synaptic pruning making them essential for normal development. Following injury, they mediate adaptive responses and under pathological conditions they can trigger neurodegeneration exacerbating the effect of a disease. Retinal organoids derived from human induced pluripotent stem cells (hiPSCs) are increasingly being used for a range of applications, including disease modelling, development of new therapies and in the study of retinogenesis. Despite many similarities to the retinas developed in vivo, they lack some key physiological features, including immune cells. We engineered an hiPSC co-culture system containing retinal organoids and microglia-like (iMG) cells and tested their retinal invasion capacity and function. We incorporated iMG into retinal organoids at 13 weeks and tested their effect on function and development at 15 and 22 weeks of differentiation. Our key findings showed that iMG cells were able to respond to endotoxin challenge in monocultures and when co-cultured with the organoids. We show that retinal organoids developed normally and retained their ability to generate spiking activity in response to light. Thus, this new co-culture immunocompetent in vitro retinal model provides a platform with greater relevance to the in vivo human retina.

Retinal Pigment Epithelium Cell Development: Extrapolating Basic Biology to Stem Cell Research

The retinal pigment epithelium (RPE) forms an important cellular monolayer, which contributes to the normal physiology of the eye. Damage to the RPE leads to the development of degenerative diseases, such as age-related macular degeneration (AMD). Apart from acting as a physical barrier between the retina and choroidal blood vessels, the RPE is crucial in maintaining photoreceptor (PR) and visual functions. Current clinical intervention to treat early stages of AMD includes stem cell-derived RPE transplantation, which is still in its early stages of evolution. Therefore, it becomes essential to derive RPEs which are functional and exhibit features as observed in native human RPE cells. The conventional strategy is to use the knowledge obtained from developmental studies using various animal models and stem cell-based exploratory studies to understand RPE biogenies and developmental trajectory. This article emphasises such studies and aims to present a comprehensive understanding of the basic biology, including the genetics and molecular pathways of RPE development. It encompasses basic developmental biology and stem cell-based developmental studies to uncover RPE differentiation. Knowledge of the in utero developmental cues provides an inclusive methodology required for deriving RPEs using stem cells.

Optic cup morphogenesis across species and related inborn human eye defects

The vertebrate eye is shaped as a cup, a conformation that optimizes vision and is acquired early in development through a process known as optic cup morphogenesis. Imaging living, transparent teleost embryos and mammalian stem cell-derived organoids has provided insights into the rearrangements that eye progenitors undergo to adopt such a shape. Molecular and pharmacological interference with these rearrangements has further identified the underlying molecular machineries and the physical forces involved in this morphogenetic process. In this Review, we summarize the resulting scenarios and proposed models that include common and species-specific events. We further discuss how these studies and those in environmentally adapted blind species may shed light on human inborn eye malformations that result from failures in optic cup morphogenesis, including microphthalmia, anophthalmia and coloboma.

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

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

Distance covariance entropy reveals primed states and bifurcation dynamics in single-cell RNA-Seq data

Cell-fate transitions are fundamental to development and differentiation. Studying them with single-cell omic data is important to advance our understanding of the cell-fate commitment process, yet this remains challenging. Here we present a computational method called DICE, which analyzes the entropy of expression covariation patterns and which is applicable to static and dynamically changing cell populations. Using only single-cell RNA-Seq data, DICE is able to predict multipotent primed states and their regulatory factors, which we subsequently validate with single-cell epigenomic data. DICE reveals that primed states are often defined by epigenetic regulators or pioneer factors alongside lineage-specific transcription factors. In developmental time course single-cell RNA-Seq datasets, DICE can pinpoint the timing of bifurcations more precisely than lineage-trajectory inference algorithms or competing variance-based methods. In summary, by studying the dynamic changes of expression covariation entropy, DICE can help elucidate primed states and bifurcation dynamics without the need for single-cell epigenomic data.

Advances in Single-Cell Sequencing Technology and Its Application in Poultry Science

Single-cell sequencing (SCS) uses a single cell as the research material and involves three dimensions: genes, phenotypes and cell biological mechanisms. This type of research can locate target cells, analyze the dynamic changes in the target cells and the relationships between the cells, and pinpoint the molecular mechanism of cell formation. Currently, a common problem faced by animal husbandry scientists is how to apply existing science and technology to promote the production of high-quality livestock and poultry products and to breed livestock for disease resistance; this is also a bottleneck for the sustainable development of animal husbandry. In recent years, although SCS technology has been successfully applied in the fields of medicine and bioscience, its application in poultry science has been rarely reported. With the sustainable development of science and technology and the poultry industry, SCS technology has great potential in the application of poultry science (or animal husbandry). Therefore, it is necessary to review the innovation of SCS technology and its application in poultry science. This article summarizes the current main technical methods of SCS and its application in poultry, which can provide potential references for its future applications in precision breeding, disease prevention and control, immunity, and cell identification.

Approaches to benchmark and characterize in vitro human model systems

In vitro human models, such as gastruloids and organoids, are complex three-dimensional (3D) structures often consist of cells from multiple germ layers that possess some attributes of a developing embryo or organ. To use these models to interrogate human development and organogenesis, these 3D models must accurately recapitulate aspects of their in vivo counterparts. Recent advances in single-cell technologies, including sequencing and spatial approaches, have enabled efforts to better understand and directly compare organoids with native tissues. For example, single-cell genomic efforts have created cell and organ atlases that enable benchmarking of in vitro models and can also be leveraged to gain novel biological insights that can be used to further improve in vitro models. This Spotlight discusses the state of current in vitro model systems, the efforts to create large publicly available atlases of the developing human and how these data are being used to improve organoids. Limitations and perspectives on future efforts are also discussed.

Microtubule modification defects underlie cilium degeneration in cell models of retinitis pigmentosa associated with pre-mRNA splicing factor mutations

Retinitis pigmentosa (RP) is the most common cause of hereditary blindness, and may occur in isolation as a non-syndromic condition or alongside other features in a syndromic presentation. Biallelic or monoallelic mutations in one of eight genes encoding pre-mRNA splicing factors are associated with non-syndromic RP. The molecular mechanism of disease remains incompletely understood, limiting opportunities for targeted treatment. Here we use CRISPR and base edited PRPF6 and PRPF31 mutant cell lines, and publicly-available data from human PRPF31^+/− patient derived retinal organoids and PRPF31 siRNA-treated organotypic retinal cultures to confirm an enrichment of differential splicing of microtubule, centrosomal, cilium and DNA damage response pathway genes in these cells. We show that genes with microtubule/centrosome/centriole/cilium gene ontology terms are enriched for weak 3′ and 5′ splice sites, and that subtle defects in spliceosome activity predominantly affect efficiency of splicing of these exons. We suggest that the primary defect in PRPF6 or PRPF31 mutant cells is microtubule and centrosomal defects, leading to defects in cilium and mitotic spindle stability, with the latter leading to DNA damage, triggering differential splicing of DNA damage response genes to activate this pathway. Finally, we expand understanding of “splicing factor RP” by investigating the function of TTLL3, one of the most statistically differentially expressed genes in PRPF6 and PRPF31 mutant cells. We identify that TTLL3 is the only tubulin glycylase expressed in the human retina, essential for monoglycylation of microtubules of the cilium, including the retinal photoreceptor cilium, to prevent cilium degeneration and retinal degeneration. Our preliminary data suggest that rescue of tubulin glycylation through overexpression of TTLL3 is sufficient to rescue cilium number in PRPF6 and PRPF31 mutant cells, suggesting that this defect underlies the cellular defect and may represent a potential target for therapeutic intervention in this group of disorders.

Single-Cell RNA Sequencing Reveals Molecular Features of Heterogeneity in the Murine Retinal Pigment Epithelium

Transcriptomic analysis of the mammalian retinal pigment epithelium (RPE) aims to identify cellular networks that influence ocular development, maintenance, function, and disease. However, available evidence points to RPE cell heterogeneity within native tissue, which adds complexity to global transcriptomic analysis. Here, to assess cell heterogeneity, we performed single-cell RNA sequencing of RPE cells from two young adult male C57BL/6J mice. Following quality control to ensure robust transcript identification limited to cell singlets, we detected 13,858 transcripts among 2667 and 2846 RPE cells. Dimensional reduction by principal component analysis and uniform manifold approximation and projection revealed six distinct cell populations. All clusters expressed transcripts typical of RPE cells; the smallest (C1, containing 1-2% of total cells) exhibited the hallmarks of stem and/or progenitor (SP) cells. Placing C1-6 along a pseudotime axis suggested a relative decrease in melanogenesis and SP gene expression and a corresponding increase in visual cycle gene expression upon RPE maturation. K-means clustering of all detected transcripts identified additional expression patterns that may advance the understanding of RPE SP cell maintenance and the evolution of cellular metabolic networks during development. This work provides new insights into the transcriptome of the mouse RPE and a baseline for identifying experimentally induced transcriptional changes in future studies of this tissue.

Single-cell atlas of keratoconus corneas revealed aberrant transcriptional signatures and implicated mechanical stretch as a trigger for keratoconus pathogenesis

Keratoconus is a common ectatic corneal disorder in adolescents and young adults that can lead to progressive visual impairment or even legal blindness. Despite the high prevalence, its etiology is not fully understood. In this study, we performed single-cell RNA sequencing (scRNA-Seq) analysis on 39,214 cells from central corneas of patients with keratoconus and healthy individuals, to define the involvement of each cell type during disease progression. We confirmed the central role of corneal stromal cells in this disease, where dysregulation of collagen and extracellular matrix (ECM) occurred. Differential gene expression and histological analyses revealed two potential novel markers for keratoconus stromal cells, namely CTSD and CTSK. Intriguingly, we detected elevated levels of YAP1 and TEAD1, the master regulators of biomechanical homeostasis, in keratoconus stromal cells. Cyclical mechanical experiments implicated the mechanical stretch in prompting protease production in corneal stromal cells during keratoconus progression. In the epithelial cells of keratoconus corneas, we observed reduced basal cells and abnormally differentiated superficial cells, unraveling the corneal epithelial lesions that were usually neglected in clinical diagnosis. In addition, several elevated cytokines in immune cells of keratoconus samples supported the involvement of inflammatory response in the progression of keratoconus. Finally, we revealed the dysregulated cell-cell communications in keratoconus, and found that only few ligand-receptor interactions were gained but a large fraction of interactional pairs was erased in keratoconus, especially for those related to protease inhibition and anti-inflammatory process. Taken together, this study facilitates the understanding of molecular mechanisms underlying keratoconus pathogenesis, providing insights into keratoconus diagnosis and potential interventions.

Molecular profiling of stem cell-derived retinal pigment epithelial cell differentiation established for clinical translation

Human embryonic stem cell-derived retinal pigment epithelial cells (hESC-RPE) are a promising cell source to treat age-related macular degeneration (AMD). Despite several ongoing clinical studies, a detailed mapping of transient cellular states during in vitro differentiation has not been performed. Here, we conduct single-cell transcriptomic profiling of an hESC-RPE differentiation protocol that has been developed for clinical use. Differentiation progressed through a culture diversification recapitulating early embryonic development, whereby cells rapidly acquired a rostral embryo patterning signature before converging toward the RPE lineage. At intermediate steps, we identified and examined the potency of an NCAM1⁺ retinal progenitor population and showed the ability of the protocol to suppress non-RPE fates. We demonstrated that the method produces a pure RPE pool capable of maturing further after subretinal transplantation in a large-eyed animal model. Our evaluation of hESC-RPE differentiation supports the development of safe and efficient pluripotent stem cell-based therapies for AMD.

•

Transcriptional analysis of hESC-RPE differentiation benchmarked to in vivo cells

•

NCAM1 emerges as a cell-surface marker of multipotent neuroepithelial progenitors

•

hESC-RPE cells are obtained through a divergence-convergence process


•

hESC-RPE further mature in vivo upon subretinal injection into the rabbit eye

Retinal ganglion cell-specific genetic regulation in primary open-angle glaucoma

To assess the transcriptomic profile of disease-specific cell populations, fibroblasts from patients with primary open-angle glaucoma (POAG) were reprogrammed into induced pluripotent stem cells (iPSCs) before being differentiated into retinal organoids and compared with those from healthy individuals. We performed single-cell RNA sequencing of a total of 247,520 cells and identified cluster-specific molecular signatures. Comparing the gene expression profile between cases and controls, we identified novel genetic associations for this blinding disease. Expression quantitative trait mapping identified a total of 4,443 significant loci across all cell types, 312 of which are specific to the retinal ganglion cell subpopulations, which ultimately degenerate in POAG. Transcriptome-wide association analysis identified genes at loci previously associated with POAG, and analysis, conditional on disease status, implicated 97 statistically significant retinal ganglion cell-specific expression quantitative trait loci. This work highlights the power of large-scale iPSC studies to uncover context-specific profiles for a genetically complex disease.

Single-cell RNA-sequencing of mammalian brain development: insights and future directions

Understanding human brain development is of fundamental interest but is also very challenging. Single-cell RNA-sequencing studies in mammals have revealed that brain development is a highly dynamic process with tremendous, previously concealed, cellular heterogeneity. This Spotlight discusses key insights from these studies and their implications for experimental models. We survey published single-cell RNA-sequencing studies of mouse and human brain development, organized by anatomical regions and developmental time points. We highlight remaining gaps in the field, predominantly concerning human brain development. We propose future directions to fill the remaining gaps, and necessary complementary techniques to create an atlas integrated in space and time of human brain development.

Single-cell characterization of malignant phenotypes and microenvironment alteration in retinoblastoma

Retinoblastoma (RB) is the most common primary intraocular malignancy of childhood. It is known that the tumor microenvironment (TME) regulates tumorigenesis and metastasis. However, how the malignant progression in RB is determined by the heterogeneity of tumor cells and TME remains uncharacterized. Here, we conducted integrative single-cell transcriptome and whole-exome sequencing analysis of RB patients with detailed pathological and clinical measurements. By single-cell transcriptomic sequencing, we profiled around 70,000 cells from tumor samples of seven RB patients. We identified that the major cell types in RB were cone precursor-like (CP-like) and MKI67+ cone precursor (MKI67+ CP) cells. By integrating copy number variation (CNV) analysis, we found that RB samples had large clonal heterogeneity, where the malignant MKI67+ CP cells had significantly larger copy number changes. Enrichment analysis revealed that the conversion of CP-like to MKI67+ CP resulted in the loss of photoreceptor function and increased cell proliferation ability. The TME in RB was composed of tumor-associated macrophages (TAMs), astrocyte-like, and cancer-associated fibroblasts (CAFs). Particularly, during the invasion process, TAMs created an immunosuppressive environment, in which the proportion of TAMs decreased, M1-type macrophage was lost, and the TAMs-related immune functions were depressed. Finally, we identified that TAMs regulated tumor cells through GRN and MIF signaling pathways, while TAMs self-regulated through inhibition of CCL and GALECTIN signaling pathways during the invasion process. Altogether, our study creates a detailed transcriptomic map of RB with single-cell characterization of malignant phenotypes and provides novel molecular insights into the occurrence and progression of RB.

Single-Cell RNA-Seq Reveals the Cellular Diversity and Developmental Characteristics of the Retinas of an Infant and a Young Child

The human retina, located in the innermost layer of the eye, plays a decisive role in visual perception. Dissecting the heterogeneity of retinal cells is essential for understanding the mechanism of visual development. Here, we performed single-cell RNA-seq to analyze 194,967 cells from the donors of infants and young children, resulting in 17 distinct clusters representing major cell types in the retina: rod photoreceptors (PRs), cone PRs, bipolar cells (BCs), horizontal cells (HCs), amacrine cells (ACs), retinal ganglion cells (RGCs), Müller glial cells (MGs), microglia, and astrocytes (ASTs). Through reclustering, we identified known subtypes of cone PRs as well as additional unreported subpopulations and corresponding markers in BCs. Additionally, we linked inherited retinal diseases (IRDs) to certain cell subtypes or subpopulations through enrichment analysis. We next constructed extensive intercellular communication networks and identified ligand-receptor interactions that play crucial roles in regulating neural cell development and immune homeostasis in the retina. Intriguingly, we found that the status and functions of PRs changed drastically between the young children and adult retina. Overall, our study offers the first retinal cell atlas in infants and young children dissecting the heterogeneity of the retina and identifying the key molecules in the developmental process, which provides an important resource that will pave the way for research on retinal development mechanisms and advancements in regenerative medicine concerning retinal biology.

Heterogeneity of human corneal endothelium implicates lncRNA NEAT1 in Fuchs endothelial corneal dystrophy

The corneal endothelium is critical for maintaining corneal clarity by mediating hydration through barrier and pump functions. Progressive loss of corneal endothelial cells during aging has been associated with the development of Fuchs endothelial corneal dystrophy (FECD), one of the main causes of cornea-related vision loss. The mechanisms underlying FECD development remain elusive. Single-cell RNA sequencing of isolated healthy human corneas discovered 4 subpopulations of corneal endothelial cells with distinctive signatures. Unsupervised clustering analysis uncovered nuclear enriched abundant transcript 1 (NEAT1), a long non-coding RNA (lncRNA), as the top expressed gene in the C0-endothelial subpopulation, but markedly downregulated in FECD. Consistent with human corneas, a UVA-induced mouse FECD model validated the loss of NEAT1 expression. Loss of NEAT1 function by an in vivo genetic approach reproduced the exacerbated phenotype of FECD by ablating corneal endothelial cells. Conversely, gain of function by a CRISPR-activated adenoviral delivery system protected corneas from UVA-induced FECD. Our findings provide novel mechanistic insights into the development of FECD, and targeting NEAT1 offers an attractive approach for treating FECD.

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.

A Single-Cell Transcriptome Atlas of the Human Retinal Pigment Epithelium

Human retinal pigment epithelium cells are arranged in a monolayer that plays an important supporting role in the retina. Although the heterogeneity of specific retinal cells has been well studied, the diversity of hRPE cells has not been reported. Here, we performed a single-cell RNA sequencing on 9,302 hRPE cells from three donors and profiled a transcriptome atlas. Our results identified two subpopulations that exhibit substantial differences in gene expression patterns and functions. One of the clusters specifically expressed ID3, a macular retinal pigment epithelium marker. The other cluster highly expressed CRYAB, a peripheral RPE marker. Our results also showed that the genes associated with oxidative stress and endoplasmic reticulum stress were more enriched in the macular RPE. The genes related to light perception, oxidative stress and lipid metabolism were more enriched in the peripheral RPE. Additionally, we provided a map of disease-related genes in the hRPE and highlighted the importance of the macular RPE and peripheral RPE clusters P4 and P6 as potential therapeutic targets for retinal diseases. Our study provides a transcriptional landscape for the human retinal pigment epithelium that is critical to understanding retinal biology and disease.

Transcriptional profiling of the chick retina identifies down-regulation of VIP and UTS2B genes during early lens-induced myopia

Gene expression of the chick retina was examined during the early development of lens-induced myopia (LIM) using whole transcriptome sequencing. Monocular treatment of the right eyes with −10 diopter (D)...

Building the mega single-cell transcriptome ocular meta-atlas

BACKGROUND

The development of highly scalable single-cell transcriptome technology has resulted in the creation of thousands of datasets, >30 in the retina alone. Analyzing the transcriptomes between different projects is highly desirable because this would allow for better assessment of which biological effects are consistent across independent studies. However it is difficult to compare and contrast data across different projects because there are substantial batch effects from computational processing, single-cell technology utilized, and the natural biological variation. While many single-cell transcriptome-specific batch correction methods purport to remove the technical noise, it is difficult to ascertain which method functions best.

RESULTS

We developed a lightweight R package (scPOP, single-cell Pick Optimal Parameters) that brings in batch integration methods and uses a simple heuristic to balance batch merging and cell type/cluster purity. We use this package along with a Snakefile-based workflow system to demonstrate how to optimally merge 766,615 cells from 33 retina datsets and 3 species to create a massive ocular single-cell transcriptome meta-atlas.

CONCLUSIONS

This provides a model for how to efficiently create meta-atlases for tissues and cells of interest.