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.

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.

An integrated transcriptional analysis of the developing human retina.

The scarcity of embryonic/foetal material as a resource for direct study means that there is still limited understanding of human retina development. Here, we present an integrated transcriptome analysis combined with immunohistochemistry in human eye and retinal samples from 4 to 19 post-conception weeks. This analysis reveals three developmental windows with specific gene expression patterns that informed the sequential emergence of retinal cell types and enabled identification of stage-specific cellular and biological processes, and transcriptional regulators. Each stage is characterised by a specific set of alternatively spliced transcripts that code for proteins involved in the formation of the photoreceptor connecting cilium, pre-mRNA splicing and epigenetic modifiers. Importantly, our data show that the transition from foetal to adult retina is characterised by a large increase in the percentage of mutually exclusive exons that code for proteins involved in photoreceptor maintenance. The circular RNA population is also defined and shown to increase during retinal development. Collectively, these data increase our understanding of human retinal development and the pre-mRNA splicing process, and help to identify new candidate disease genes.

IGFBPs mediate IGF-1's functions in retinal lamination and photoreceptor development during pluripotent stem cell differentiation to retinal organoids.

Development of the retina is regulated by growth factors, such as insulin-like growth factors 1 and 2 (IGF-1/2), which coordinate proliferation, differentiation, and maturation of the neuroepithelial precursors cells. In the circulation, IGF-1/2 are transported by the insulin growth factor binding proteins (IGFBPs) family members. IGFBPs can impact positively and negatively on IGF-1, by making it available or sequestering IGF-1 to or from its receptor. In this study, we investigated the expression of IGFBPs and their role in the generation of human retinal organoids from human pluripotent stem cells, showing a dynamic expression pattern suggestive of different IGFBPs being used in a stage-specific manner to mediate IGF-1 functions. Our data show that IGF-1 addition to culture media facilitated the generation of retinal organoids displaying the typical laminated structure and photoreceptor maturation. The organoids cultured in the absence of IGF-1, lacked the typical laminated structure at the early stages of differentiation and contained significantly less photoreceptors and more retinal ganglion cells at the later stages of differentiation, confirming the positive effects of IGF-1 on retinal lamination and photoreceptor development. The organoids cultured with the IGFBP inhibitor (NBI-31772) and IGF-1 showed lack of retinal lamination at the early stages of differentiation, an increased propensity to generate horizontal cells at mid-stages of differentiation and reduced photoreceptor development at the later stages of differentiation. Together these data suggest that IGFBPs enable IGF-1's role in retinal lamination and photoreceptor development in a stage-specific manner.

Transplanted pluripotent stem cell-derived photoreceptor precursors elicit conventional and unusual light responses in mice with advanced retinal degeneration.

Retinal dystrophies often lead to blindness. Developing therapeutic interventions to restore vision is therefore of paramount importance. Here we demonstrate the ability of pluripotent stem cell-derived cone precursors to engraft and restore light responses in the Pde6brd1 mouse, an end-stage photoreceptor degeneration model. Our data show that up to 1.5% of precursors integrate into the host retina, differentiate into cones, and engraft in close apposition to the host bipolar cells. Half of the transplanted mice exhibited visual behavior and of these 33% showed binocular light sensitivity. The majority of retinal ganglion cells exhibited contrast-sensitive ON, OFF or ON-OFF light responses and even motion sensitivity; however, quite a few exhibited unusual responses (eg, light-induced suppression), presumably reflecting remodeling of the neural retina. Our data indicate that despite relatively low engraftment yield, pluripotent stem cell-derived cone precursors can elicit light responsiveness even at advanced degeneration stages. Further work is needed to improve engraftment yield and counteract retinal remodeling to achieve useful clinical applications.

Developing a simple method to enhance the generation of cone and rod photoreceptors in pluripotent stem cell-derived retinal organoids.

Cell replacement therapy is a promising treatment for irreversible retinal cell death in diverse diseases such as Stargardt's disease, age-related macular degeneration, and retinitis pigmentosa. The final impact of all retinal dystrophies is the loss of photoreceptors; hence, there is a pressing need for research into replacement. Seminal work has shown that a simple three-dimensional culture system enables differentiation of human pluripotent stem cells to retinal organoids containing large numbers of photoreceptors developing alongside retinal neurons and Müller glia cells in a laminated structure that resembles the native retina. Despite these promising developments, current protocols show different efficiencies across pluripotent stem cells and result in retinal organoids with a mixture of photoreceptor cells at varying maturation states, along with nonphotoreceptor cell types. In this study, we investigated the impact of stage-specific addition of retinoic acid (RA), 9-cis-retinal, 11-cis-retinal, levodopa (l-DOPA), triiodothyronine (T3), and γ-secretase inhibitor ((2S)-N-[(3,5-Difluorophenyl)acetyl]-l-alanyl-2-phenyl]glycine1,1-dimethylethyl ester2L [DAPT]) in the generation of cone and rod photoreceptors. Our results indicate that addition of RA + T3 during days 90 to 120 of differentiation enhanced the generation of rod and S-cone photoreceptor formation, while the combined addition of DAPT from days 28 to 42 with RA during days 30 to 120 of differentiation led to enhanced generation of L/M-cones at the expense of rods. l-DOPA when added together with RA during days 90 to 120 of differentiation also promoted the emergence of S-cones at the expense of rod photoreceptors. Collectively, these data represent an advance in our ability to direct generation of rod and cone photoreceptors in vitro.

Deconstructing Retinal Organoids: Single Cell RNA-Seq Reveals the Cellular Components of Human Pluripotent Stem Cell-Derived Retina.

The rapid improvements in single cell sequencing technologies and analyses afford greater scope for dissecting organoid cultures composed of multiple cell types and create an opportunity to interrogate these models to understand tissue biology, cellular behavior and interactions. To this end, retinal organoids generated from human embryonic stem cells (hESCs) were analyzed by single cell RNA-sequencing (scRNA-Seq) at three time points of differentiation. Combinatorial data from all time points revealed the presence of nine clusters, five of which corresponded to key retinal cell types: retinal pigment epithelium (RPE), retinal ganglion cells (RGCs), cone and rod photoreceptors, and Müller glia. The remaining four clusters expressed genes typical of mitotic cells, extracellular matrix components and those involved in homeostasis. The cell clustering analysis revealed the decreasing presence of mitotic cells and RGCs, formation of a distinct RPE cluster, the emergence of cone and rod photoreceptors from photoreceptor precursors, and an increasing number of Müller glia cells over time. Pseudo-time analysis resembled the order of cell birth during retinal development, with the mitotic cluster commencing the trajectory and the large majority of Müller glia completing the time line. Together, these data demonstrate the feasibility and potential of scRNA-Seq to dissect the inherent complexity of retinal organoids and the orderly birth of key retinal cell types. Stem Cells 2019;37:593-598.

Human-Induced Pluripotent Stem Cells Generate Light Responsive Retinal Organoids with Variable and Nutrient-Dependent Efficiency.

The availability of in vitro models of the human retina in which to perform pharmacological and toxicological studies is an urgent and unmet need. An essential step for developing in vitro models of human retina is the ability to generate laminated, physiologically functional, and light-responsive retinal organoids from renewable and patient specific sources. We investigated five different human-induced pluripotent stem cell (iPSC) lines and showed a significant variability in their efficiency to generate retinal organoids. Despite this variability, by month 5 of differentiation, all iPSC-derived retinal organoids were able to generate light responses, albeit immature, comparable to the earliest light responses recorded from the neonatal mouse retina, close to the period of eye opening. All iPSC-derived retinal organoids exhibited at this time a well-formed outer nuclear like layer containing photoreceptors with inner segments, connecting cilium, and outer like segments. The differentiation process was highly dependent on seeding cell density and nutrient availability determined by factorial experimental design. We adopted the differentiation protocol to a multiwell plate format, which enhanced generation of retinal organoids with retinal-pigmented epithelium (RPE) and improved ganglion cell development and the response to physiological stimuli. We tested the response of iPSC-derived retinal organoids to Moxifloxacin and showed that similarly to in vivo adult mouse retina, the primary affected cell types were photoreceptors. Together our data indicate that light responsive retinal organoids derived from carefully selected and differentiation efficient iPSC lines can be generated at the scale needed for pharmacology and drug screening purposes. Stem Cells 2018;36:1535-1551.

Using Zinc Finger Nuclease Technology to Generate CRX-Reporter Human Embryonic Stem Cells as a Tool to Identify and Study the Emergence of Photoreceptors Precursors During Pluripotent Stem Cell Differentiation.

The purpose of this study was to generate human embryonic stem cell (hESC) lines harboring the green fluorescent protein (GFP) reporter at the endogenous loci of the Cone-Rod Homeobox (CRX) gene, a key transcription factor in retinal development. Zinc finger nucleases (ZFNs) designed to cleave in the 3' UTR of CRX were transfected into hESCs along with a donor construct containing homology to the target region, eGFP reporter, and a puromycin selection cassette. Following selection, polymerase chain reaction (PCR) and sequencing analysis of antibiotic resistant clones indicated targeted integration of the reporter cassette at the 3' of the CRX gene, generating a CRX-GFP fusion. Further analysis of a clone exhibiting homozygote integration of the GFP reporter was conducted suggesting genomic stability was preserved and no other copies of the targeting cassette were inserted elsewhere within the genome. This clone was selected for differentiation towards the retinal lineage. Immunocytochemistry of sections obtained from embryoid bodies and quantitative reverse transcriptase PCR of GFP positive and negative subpopulations purified by fluorescence activated cell sorting during the differentiation indicated a significant correlation between GFP and endogenous CRX expression. Furthermore, GFP expression was found in photoreceptor precursors emerging during hESC differentiation, but not in the retinal pigmented epithelium, retinal ganglion cells, or neurons of the developing inner nuclear layer. Together our data demonstrate the successful application of ZFN technology to generate CRX-GFP labeled hESC lines, which can be used to study and isolate photoreceptor precursors during hESC differentiation.

AII amacrine cells discriminate between heterocellular and homocellular locations when assembling connexin36-containing gap junctions.

Electrical synapses (gap junctions) rapidly transmit signals between neurons and are composed of connexins. In neurons, connexin36 (Cx36) is the most abundant isoform; however, the mechanisms underlying formation of Cx36-containing electrical synapses are unknown. We focus on homocellular and heterocellular gap junctions formed by an AII amacrine cell, a key interneuron found in all mammalian retinas. In mice lacking native Cx36 but expressing a variant tagged with enhanced green fluorescent protein at the C-terminus (KO-Cx36-EGFP), heterocellular gap junctions formed between AII cells and ON cone bipolar cells are fully functional, whereas homocellular gap junctions between two AII cells are not formed. A tracer injected into an AII amacrine cell spreads into ON cone bipolar cells but is excluded from other AII cells. Reconstruction of Cx36-EGFP clusters on an AII cell in the KO-Cx36-EGFP genotype confirmed that the number, but not average size, of the clusters is reduced - as expected for AII cells lacking a subset of electrical synapses. Our studies indicate that some neurons exhibit at least two discriminatory mechanisms for assembling Cx36. We suggest that employing different gap-junction-forming mechanisms could provide the means for a cell to regulate its gap junctions in a target-cell-specific manner, even if these junctions contain the same connexin.

Deciphering the spatiotemporal transcriptional and chromatin accessibility of human retinal organoid development at the single-cell level.

Molecular information on the early stages of human retinal development remains scarce due to limitations in obtaining early human eye samples. Pluripotent stem cell-derived retinal organoids (ROs) provide an unprecedented opportunity for studying early retinogenesis. Using a combination of single cell RNA-seq and spatial transcriptomics we present for the first-time a single cell spatiotemporal transcriptome of RO development. Our data demonstrate that ROs recapitulate key events of retinogenesis including optic vesicle/cup formation, presence of a putative ciliary margin zone, emergence of retinal progenitor cells and their orderly differentiation to retinal neurons. Combining the scRNA- with scATAC-seq data, we were able to reveal cell-type specific transcription factor binding motifs on accessible chromatin at each stage of organoid development, and to show that chromatin accessibility is highly correlated to the developing human retina, but with some differences in the temporal emergence and abundance of some of the retinal neurons.

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.

Retinal cells derived from patients with DRAM2-dependent CORD21 dystrophy exhibit key lysosomal enzyme deficiency and lysosomal content accumulation.

Biallelic mutations in DRAM2 lead to an autosomal recessive cone-rod dystrophy known as CORD21, which typically presents between the third and sixth decades of life. Although DRAM2 localizes to the lysosomes of photoreceptor and retinal pigment epithelium (RPE) cells, its specific role in retinal degeneration has not been fully elucidated. In this study, we generated and characterized retinal organoids (ROs) and RPE cells from induced pluripotent stem cells (iPSCs) derived from two CORD21 patients. Our investigation revealed that CORD21-ROs and RPE cells exhibit abnormalities in lipid metabolism, defects in autophagic flux, accumulation of aberrant lysosomal content, and reduced lysosomal enzyme activity. We identified potential interactions of DRAM2 with vesicular trafficking proteins, suggesting its involvement in this cellular process. These findings collectively suggest that DRAM2 plays a crucial role in maintaining the integrity of photoreceptors and RPE cells by regulating lysosomal function, autophagy, and potentially vesicular trafficking.

PRPF8-mediated dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects.

The carboxy-terminus of the spliceosomal protein PRPF8, which regulates the RNA helicase Brr2, is a hotspot for mutations causing retinitis pigmentosa-type 13, with unclear role in human splicing and tissue-specificity mechanism. We used patient induced pluripotent stem cells-derived cells, carrying the heterozygous PRPF8 c.6926 A > C (p.H2309P) mutation to demonstrate retinal-specific endophenotypes comprising photoreceptor loss, apical-basal polarity and ciliary defects. Comprehensive molecular, transcriptomic, and proteomic analyses revealed a role of the PRPF8/Brr2 regulation in 5'-splice site (5'SS) selection by spliceosomes, for which disruption impaired alternative splicing and weak/suboptimal 5'SS selection, and enhanced cryptic splicing, predominantly in ciliary and retinal-specific transcripts. Altered splicing efficiency, nuclear speckles organisation, and PRPF8 interaction with U6 snRNA, caused accumulation of active spliceosomes and poly(A)+ mRNAs in unique splicing clusters located at the nuclear periphery of photoreceptors. Collectively these elucidate the role of PRPF8/Brr2 regulatory mechanisms in splicing and the molecular basis of retinal disease, informing therapeutic approaches.

Ablation of retinal horizontal cells from adult mice leads to rod degeneration and remodeling in the outer retina.

In the brain, including the retina, interneurons show an enormous structural and functional diversity. Retinal horizontal cells represent a class of interneurons that form triad synapses with photoreceptors and ON bipolar cells. At this first retinal synapse, horizontal cells modulate signal transmission from photoreceptors to bipolar cells by feedback and feedforward inhibition. To test how the fully developed retina reacts to the specific loss of horizontal cells, these interneurons were specifically ablated from adult mice using the diphtheria toxin (DT)/DT-receptor system and the connexin57 promoter. Following ablation, the retinal network responded with extensive remodeling: rods retracted their axons from the outer plexiform layer and partially degenerated, whereas cones survived. Cone pedicles remained in the outer plexiform layer and preserved synaptic contacts with OFF but not with ON bipolar cells. Consistently, the retinal ON pathway was impaired, leading to reduced amplitudes and prolonged latencies in electroretinograms. However, ganglion cell responses showed only slight changes in time course, presumably because ON bipolar cells formed multiple ectopic synapses with photoreceptors, and visual performance, assessed with an optomotor system, was only mildly affected. Thus, the loss of an entire interneuron class can be largely compensated even by the adult retinal network.

Light-responsive microRNA molecules in human retinal organoids are differentially regulated by distinct wavelengths of light.

Cells in the human retina must rapidly adapt to constantly changing visual stimuli. This fast adaptation to varying levels and wavelengths of light helps to regulate circadian rhythms and allows for adaptation to high levels of illumination, thereby enabling the rest of the visual system to remain responsive. It has been shown that retinal microRNA (miRNA) molecules play a key role in regulating these processes. However, despite extensive research using various model organisms, light-regulated miRNAs in human retinal cells remain unknown. Here, we aim to characterize these miRNAs. We generated light-responsive human retinal organoids that express miRNA families and clusters typically found in the retina. Using an in-house developed photostimulation device, we identified a subset of light-regulated miRNAs. Importantly, we found that these miRNAs are differentially regulated by distinct wavelengths of light and have a rapid turnover, highlighting the dynamic and adaptive nature of the human retina.

pRB-Depleted Pluripotent Stem Cell Retinal Organoids Recapitulate Cell State Transitions of Retinoblastoma Development and Suggest an Important Role for pRB in Retinal Cell Differentiation.

Retinoblastoma (Rb) is a childhood cancer of the developing retina, accounting for up to 17% of all tumors in infancy. To gain insights into the transcriptional events of cell state transitions during Rb development, we established 2 disease models via retinal organoid differentiation of a pRB (retinoblastoma protein)-depleted human embryonic stem cell line (RB1-null hESCs) and a pRB patient-specific induced pluripotent (iPSC) line harboring a RB1 biallelic mutation (c.2082delC). Both models were characterized by pRB depletion and accumulation of retinal progenitor cells at the expense of amacrine, horizontal and retinal ganglion cells, which suggests an important role for pRB in differentiation of these cell lineages. Importantly, a significant increase in the fraction of proliferating cone precursors (RXRγ+Ki67+) was observed in both pRB-depleted organoid models, which were defined as Rb-like clusters by single-cell RNA-Seq analysis. The pRB-depleted retinal organoids displayed similar features to Rb tumors, including mitochondrial cristae aberrations and rosette-like structures, and were able to undergo cell growth in an anchorage-independent manner, indicative of cell transformation in vitro. In both models, the Rb cones expressed retinal ganglion and horizontal cell markers, a novel finding, which could help to better characterize these tumors with possible therapeutic implications. Application of Melphalan, Topotecan, and TW-37 led to a significant reduction in the fraction of Rb proliferating cone precursors, validating the suitability of these in vitro models for testing novel therapeutics for Rb.

Human Retinal Organoids Provide a Suitable Tool for Toxicological Investigations: A Comprehensive Validation Using Drugs and Compounds Affecting the Retina.

Retinal drug toxicity screening is essential for the development of safe treatment strategies for a large number of diseases. To this end, retinal organoids derived from human pluripotent stem cells (hPSCs) provide a suitable screening platform due to their similarity to the human retina and the ease of generation in large-scale formats. In this study, two hPSC cell lines were differentiated to retinal organoids, which comprised all key retinal cell types in multiple nuclear and synaptic layers. Single-cell RNA-Seq of retinal organoids indicated the maintenance of retinal ganglion cells and development of bipolar cells: both cell types segregated into several subtypes. Ketorolac, digoxin, thioridazine, sildenafil, ethanol, and methanol were selected as key compounds to screen on retinal organoids because of their well-known retinal toxicity profile described in the literature. Exposure of the hPSC-derived retinal organoids to digoxin, thioridazine, and sildenafil resulted in photoreceptor cell death, while digoxin and thioridazine additionally affected all other cell types, including Müller glia cells. All drug treatments caused activation of astrocytes, indicated by dendrites sprouting into neuroepithelium. The ability to respond to light was preserved in organoids although the number of responsive retinal ganglion cells decreased after drug exposure. These data indicate similar drug effects in organoids to those reported in in vivo models and/or in humans, thus providing the first robust experimental evidence of their suitability for toxicological studies.

Co-expression of SARS-CoV-2 entry genes in the superficial adult human conjunctival, limbal and corneal epithelium suggests an additional route of entry via the ocular surface.

PURPOSE: The high infection rate of SARS-CoV-2 necessitates the need for multiple studies identifying the molecular mechanisms that facilitate the viral entry and propagation. Currently the potential extra-respiratory transmission routes of SARS-CoV-2 remain unclear. METHODS: Using single-cell RNA Seq and ATAC-Seq datasets and immunohistochemical analysis, we investigated SARS-CoV-2 tropism in the embryonic, fetal and adult human ocular surface. RESULTS: The co-expression of ACE2 receptor and entry protease TMPRSS2 was detected in the human adult conjunctival, limbal and corneal epithelium, but not in the embryonic and fetal ocular surface up to 21 post conception weeks. These expression patterns were corroborated by the single cell ATAC-Seq data, which revealed a permissive chromatin in ACE2 and TMPRSS2 loci in the adult conjunctival, limbal and corneal epithelium. Co-expression of ACE2 and TMPRSS2 was strongly detected in the superficial limbal, corneal and conjunctival epithelium, implicating these as target entry cells for SARS-CoV-2 in the ocular surface. Strikingly, we also identified the key pro-inflammatory signals TNF, NFKß and IFNG as upstream regulators of the transcriptional profile of ACE2+TMPRSS2+ cells in the superficial conjunctival epithelium, suggesting that SARS-CoV-2 may utilise inflammatory driven upregulation of ACE2 and TMPRSS2 expression to enhance infection in ocular surface. CONCLUSIONS: Together our data indicate that the human ocular surface epithelium provides an additional entry portal for SARS-CoV-2, which may exploit inflammatory driven upregulation of ACE2 and TMPRSS2 entry factors to enhance infection.

A single cell atlas of human cornea that defines its development, limbal progenitor cells and their interactions with the immune cells.

PURPOSE: Single cell (sc) analyses of key embryonic, fetal and adult stages were performed to generate a comprehensive single cell atlas of all the corneal and adjacent conjunctival cell types from development to adulthood. METHODS: Four human adult and seventeen embryonic and fetal corneas from 10 to 21 post conception week (PCW) specimens were dissociated to single cells and subjected to scRNA- and/or ATAC-Seq using the 10x Genomics platform. These were embedded using Uniform Manifold Approximation and Projection (UMAP) and clustered using Seurat graph-based clustering. Cluster identification was performed based on marker gene expression, bioinformatic data mining and immunofluorescence (IF) analysis. RNA interference, IF, colony forming efficiency and clonal assays were performed on cultured limbal epithelial cells (LECs). RESULTS: scRNA-Seq analysis of 21,343 cells from four adult human corneas and adjacent conjunctivas revealed the presence of 21 cell clusters, representing the progenitor and differentiated cells in all layers of cornea and conjunctiva as well as immune cells, melanocytes, fibroblasts, and blood/lymphatic vessels. A small cell cluster with high expression of limbal progenitor cell (LPC) markers was identified and shown via pseudotime analysis to give rise to five other cell types representing all the subtypes of differentiated limbal and corneal epithelial cells. A novel putative LPCs surface marker, GPHA2, expressed on the surface of 0.41% ± 0.21 of the cultured LECs, was identified, based on predominant expression in the limbal crypts of adult and developing cornea and RNAi validation in cultured LECs. Combining scRNA- and ATAC-Seq analyses, we identified multiple upstream regulators for LPCs and demonstrated a close interaction between the immune cells and limbal progenitor cells. RNA-Seq analysis indicated the loss of GPHA2 expression and acquisition of proliferative limbal basal epithelial cell markers during ex vivo LEC expansion, independently of the culture method used. Extending the single cell analyses to keratoconus, we were able to reveal activation of collagenase in the corneal stroma and a reduced pool of limbal suprabasal cells as two key changes underlying the disease phenotype. Single cell RNA-Seq of 89,897 cells obtained from embryonic and fetal cornea indicated that during development, the conjunctival epithelium is the first to be specified from the ocular surface epithelium, followed by the corneal epithelium and the establishment of LPCs, which predate the formation of limbal niche by a few weeks. CONCLUSIONS: Our scRNA-and ATAC-Seq data of developing and adult cornea in steady state and disease conditions provide a unique resource for defining genes/pathways that can lead to improvement in ex vivo LPCs expansion, stem cell differentiation methods and better understanding and treatment of ocular surface disorders.