Molecular characteristics and spatial distribution of adult human corneal cell subtypes

Cell therapy in the cornea: the emerging role of microenvironment

The cornea is an ideal testing field for cell therapies. Its highly ordered structure, where specific cell populations are sequestered in different layers, together with its accessibility, has allowed the development of the first stem cell-based therapy approved by the European Medicine Agency. Today, different techniques have been proposed for autologous and allogeneic limbal and non-limbal cell transplantation. Cell replacement has also been attempted in cases of endothelial cell decompensation as it occurs in Fuchs dystrophy: injection of cultivated allogeneic endothelial cells is now in advanced phases of clinical development. Recently, stromal substitutes have been developed with excellent integration capability and transparency. Finally, cell-derived products, such as exosomes obtained from different sources, have been investigated for the treatment of severe corneal diseases with encouraging results. Optimization of the success rate of cell therapies obviously requires high-quality cultured cells/products, but the role of the surrounding microenvironment is equally important to allow engraftment of transplanted cells, to preserve their functions and, ultimately, lead to restoration of tissue integrity and transparency of the cornea.

Single nuclei transcriptomics of the in situ human limbal stem cell niche

The corneal epithelium acts as a barrier to pathogens entering the eye; corneal epithelial cells are continuously renewed by uni-potent, quiescent limbal stem cells (LSCs) located at the limbus, where the cornea transitions to conjunctiva. There has yet to be a consensus on LSC markers and their transcriptome profile is not fully understood, which may be due to using cadaveric tissue without an intact stem cell niche for transcriptomics. In this study, we addressed this problem by using single nuclei RNA sequencing (snRNAseq) on healthy human limbal tissue that was immediately snap-frozen after excision from patients undergoing cataract surgery. We identified the quiescent LSCs as a sub-population of corneal epithelial cells with a low level of total transcript counts. Moreover, TP63, KRT15, CXCL14, and ITGβ4 were found to be highly expressed in LSCs and transiently amplifying cells (TACs), which constitute the corneal epithelial progenitor populations at the limbus. The surface markers SLC6A6 and ITGβ4 could be used to enrich human corneal epithelial cell progenitors, which were also found to specifically express the putative limbal progenitor cell markers MMP10 and AC093496.1.

A Keratin 12 Expression-Based Analysis of Stem-Precursor Cells and Differentiation in the Limbal-Corneal Epithelium Using Single-Cell RNA-Seq Data

The corneal epithelium (CE) is spread between two domains, the outer vascularized limbus and the avascular cornea proper. Epithelial cells undergo constant migration from the limbus to the vision-critical central cornea. Coordinated with this migration, the cells undergo differentiation changes where a pool of unique stem/precursor cells at the limbus yields the mature cells that reach the corneal center. Differentiation is heralded by the expression of the corneal-specific Krt12. Processing data acquired by scRNA-Seq showed that the increase in Krt12 expression occurs in four distinct steps within the limbus, plus a single continuous increase in the cornea. Differential gene analysis demonstrated that these domains reflect discreet stages of CE differentiation and yielded extensive information of the genes undergoing down- or upregulation in the sequential transition from less to more differentiate conditions. The approach allowed the identification of multiple gene cohorts, including (a) the genes which have maximal expression in the most primitive, Krt12-negative cell cohort, which is likely to include the stem/precursor cells; (b) the sets of genes that undergo continuous increase or decrease along the whole differentiation path; and (c) the genes showing maximal positive or negative correlation with the changes in Krt12.

Eyes open on stem cells

Recently, the murine cornea has reemerged as a robust stem cell (SC) model, allowing individual SC tracing in living animals. The cornea has pioneered seminal discoveries in SC biology and regenerative medicine, from the first corneal transplantation in 1905 to the identification of limbal SCs and their transplantation to successfully restore vision in the early 1990s. Recent experiments have exposed unexpected properties attributed to SCs and progenitors and revealed flexibility in the differentiation program and a key role for the SC niche. Here, we discuss the limbal SC model and its broader relevance to other tissues, disease, and therapy.

Single-cell transcriptomics of the ocular anterior segment: a comprehensive review

Elucidating the cellular and genetic composition of ocular tissues is essential for uncovering the pathophysiology of ocular diseases. Since the introduction of single-cell RNA sequencing (scRNA-seq) in 2009, vision researchers have performed extensive single-cell analyses to better understand transcriptome complexity and heterogeneity of ocular structures. This technology has revolutionized our ability to identify rare cell populations and to make cross-species comparisons of gene expression in both steady state and disease conditions. Importantly, single-cell transcriptomic analyses have enabled the identification of cell-type specific gene markers and signalling pathways between ocular cell populations. While most scRNA-seq studies have been conducted on retinal tissues, large-scale transcriptomic atlases pertaining to the ocular anterior segment have also been constructed in the past three years. This timely review provides vision researchers with an overview of scRNA-seq experimental design, technical limitations, and clinical applications in a variety of anterior segment-related ocular pathologies. We review open-access anterior segment-related scRNA-seq datasets and illustrate how scRNA-seq can be an indispensable tool for the development of targeted therapeutics.

Single-Cell Transcriptomics Reveals Cellular Heterogeneity and Complex Cell-Cell Communication Networks in the Mouse Cornea

Purpose

To generate a single-cell RNA-sequencing (scRNA-seq) map and construct cell-cell communication networks of mouse corneas.

Methods

C57BL/6 mouse corneas were dissociated to single cells and subjected to scRNA-seq. Cell populations were clustered and annotated for bioinformatic analysis using the R package "Seurat." Differential expression patterns were validated and spatially mapped with whole-mount immunofluorescence staining. Global intercellular signaling networks were constructed using CellChat.

Results

Unbiased clustering of scRNA-seq transcriptomes of 14,732 cells from 40 corneas revealed 17 cell clusters of six major cell types: nine epithelial cell, three keratocyte, two corneal endothelial cell, and one each of immune cell, vascular endothelial cell, and fibroblast clusters. The nine epithelial cell subtypes included quiescent limbal stem cells, transit-amplifying cells, and differentiated cells from corneas and two minor conjunctival epithelial clusters. CellChat analysis provided an atlas of the complex intercellular signaling communications among all cell types.

Conclusions

We constructed a complete single-cell transcriptomic map and the complex signaling cross-talk among all cell types of the cornea, which can be used as a foundation atlas for further research on the cornea. This study also deepens the understanding of the cellular heterogeneity and heterotypic cell-cell interaction within corneas.

The miR-183/96/182 cluster is a checkpoint for resident immune cells and shapes the cellular landscape of the cornea

PURPOSE

The conserved miR-183/96/182 cluster (miR-183C) regulates both corneal sensory innervation and corneal resident immune cells (CRICs). This study is to uncover its role in CRICs and in shaping the corneal cellular landscape at a single-cell (sc) level.

METHODS

Corneas of naïve, young adult [2 and 6 months old (mo)], female miR-183C knockout (KO) mice and wild-type (WT) littermates were harvested and dissociated into single cells. Dead cells were removed using a Dead Cell Removal kit. CD45+ CRICs were enriched by Magnetic Activated Cell Sorting (MACS). scRNA libraries were constructed and sequenced followed by comprehensive bioinformatic analyses.

RESULTS

The composition of major cell types of the cornea stays relatively stable in WT mice from 2 to 6 mo, however the compositions of subtypes of corneal cells shift with age. Inactivation of miR-183C disrupts the stability of the major cell-type composition and age-related transcriptomic shifts of subtypes of corneal cells. The diversity of CRICs is enhanced with age. Naïve mouse cornea contains previously-unrecognized resident fibrocytes and neutrophils. Resident macrophages (ResMφ) adopt cornea-specific function by expressing abundant extracellular matrix (ECM) and ECM organization-related genes. Naïve cornea is endowed with partially-differentiated proliferative ResMφ and contains microglia-like Mφ. Resident lymphocytes, including innate lymphoid cells (ILCs), NKT and γδT cells, are the major source of innate IL-17a. miR-183C limits the diversity and polarity of ResMφ.

CONCLUSION

miR-183C serves as a checkpoint for CRICs and imposes a global regulation of the cellular landscape of the cornea.

PLAE Web App Enables Powerful Searching and Multiple Visualizations Across One Million Unified Single-Cell Ocular Transcriptomes

Purpose

To create a high-performance reactive web application to query single-cell gene expression data across cell type, species, study, and other factors.

Methods

We updated the content and structure of the underlying data (single cell Eye in a Disk [scEiaD]) and wrote the web application PLAE (https://plae.nei.nih.gov) to visualize and explore the data.

Results

The new portal provides quick visualization of over a million individual cells from vertebrate eye and body transcriptomes encompassing four species, 60 cell types, six ocular tissues, and 23 body tissues across 35 publications. To demonstrate the value of this unified pan-eye dataset, we replicated known neurogenic and cone macula markers in addition to proposing six new cone human region markers.

Conclusions

The PLAE web application offers the eye community a powerful and quick means to test hypotheses related to gene expression across a highly diverse, community-derived database.

Translational Relevance

The PLAE resource enables any researcher or clinician to study and research gene expression patterning across a wide variety of curated ocular cell types with a responsive web app.

Single-Cell RNA Sequencing: Opportunities and Challenges for Studies on Corneal Biology in Health and Disease

The structure and major cell types of the multi-layer human cornea have been extensively studied. However, various cell states in specific cell types and key genes that define the cell states are not fully understood, hindering our comprehension of corneal homeostasis, related diseases, and therapeutic discovery. Single-cell RNA sequencing is a revolutionary and powerful tool for identifying cell states within tissues such as the cornea. This review provides an overview of current single-cell RNA sequencing studies on the human cornea, highlighting similarities and differences between them, and summarizing the key genes that define corneal cell states reported in these studies. In addition, this review discusses the opportunities and challenges of using single-cell RNA sequencing to study corneal biology in health and disease.

Systematic review of SLC4A11, ZEB1, LOXHD1, and AGBL1 variants in the development of Fuchs' endothelial corneal dystrophy

Introduction

The pathogenic role of variants in TCF4 and COL8A2 in causing Fuchs' endothelial corneal dystrophy (FECD) is not controversial and has been confirmed by numerous studies. The causal role of other genes, SLC4A11, ZEB1, LOXHD1, and AGBL1, which have been reported to be associated with FECD, is more complicated and less obvious. We performed a systematic review of the variants in the above-mentioned genes in FECD cases, taking into account the currently available population frequency information, transcriptomic data, and the results of functional studies to assess their pathogenicity.

Methods

Search for articles published in 2005-2022 was performed manually between July 2022 and February 2023. We searched for original research articles in peer-reviewed journals, written in English. Variants in the genes of interest identified in patients with FECD were extracted for the analysis. We classified each presented variant by pathogenicity status according to the ACMG criteria implemented in the Varsome tool. Diagnosis, segregation data, presence of affected relatives, functional analysis results, and gene expression in the corneal endothelium were taken into account. Data on the expression of genes of interest in the corneal endothelium were extracted from articles in which transcriptome analysis was performed. The identification of at least one variant in a gene classified as pathogenic or significantly associated with FECD was required to confirm the causal role of the gene in FECD.

Results

The analysis included 34 articles with 102 unique ZEB1 variants, 20 articles with 64 SLC4A11 variants, six articles with 26 LOXHD1 variants, and five articles with four AGBL1 variants. Pathogenic status was confirmed for seven SLC4A11 variants found in FECD. No variants in ZEB1, LOXHD1, and AGBL1 genes were classified as pathogenic for FECD. According to the transcriptome data, AGBL1 and LOXHD1 were not expressed in the corneal endothelium. Functional evidence for the association of LOXHD1, and AGBL1 with FECD was conflicting.

Conclusion

Our analysis confirmed the causal role of SLC4A11 variants in the development of FECD. The causal role of ZEB1, LOXHD1, and AGBL1 variants in FECD has not been confirmed. Further evidence from familial cases and functional analysis is needed to confirm their causal roles in FECD.

Optimization of method for achieving a single-cell suspension from mouse corneas

The single-cell RNA-sequencing (scRNA-seq) technique is used to explore the biological characteristics of tissues under pathological and physiological conditions that include certain chronic eye diseases. Harvesting of single-cell suspensions is one challenge inherent to scRNA-seq procedures. This study aimed to use an optimized method to digest a whole mouse cornea to harvest single-cell suspensions. We utilized five different mouse cornea digestion methods to obtain single-cell suspensions: (1) 5 dissected mouse corneas were cut into pieces (∼0.5 mm) and digested in trypsin for 10 min, and this digestion was repeated for 10 cycles; (2) 5 dissected mouse corneas were cut into pieces and incubated with 5 mg/ml collagenase A at 37 °C for 1h and then further digested in trypsin at 37 °C for 10 min; (3) used the same approach as that used in method 2, but the second digestion step was performed in TrypLE for 20 min; (4) used the same approach as that used in method 2, but the concentration of collagenase A was 2 mg/ml and the incubation time was 2h; (5) used the same approach as that used in method 3, but the corneas were incubated in 2 mg/ml collagenase A at 37 °C for 2h. Trypan blue staining was used to calculate the cell viability and agglomeration rate. The cell types and percentages were determined using immunofluorescence staining. RNA integrity number (RIN) was measured by Agilent 2100. Method 1 showed the lowest cell yield (0.375 × 10⁵), epithelial cell percentage, and less than 70% cell viability, thus not a proper protocol. Method 2 showed the highest cell viability (over 90%), percentage of single-cell (89.53%), and high cell quantity (1.05 × 10⁵). Method 3 had a significantly lower cell viability (55.30%). Cell agglomeration rates of method 4 and 5 reached up to 20% and 13%, and with lower cell viability (72.51%, 59.87%, respectively) and decreased epithelial cell rate compared to method 2 and 3. The results suggest that method 2 (5 mg/ml collagenase A and trypsin) is a preferred protocol for digesting mouse cornea to obtain single-cell suspension which achieves the criterion of single-cell RNA sequencing.

Characterization of Porcine Ocular Surface Epithelial Microenvironment

The porcine ocular surface is used as a model of the human ocular surface; however, a detailed characterization of the porcine ocular surface has not been documented. This is due, in part, to the scarcity of antibodies produced specifically against the porcine ocular surface cell types or structures. We performed a histological and immunohistochemical investigation on frozen and formalin-fixed, paraffin-embedded ocular surface tissue from domestic pigs using a panel of 41 different antibodies related to epithelial progenitor/differentiation phenotypes, extracellular matrix and associated molecules, and various niche cell types. Our observations suggested that the Bowman's layer is not evident in the cornea; the deep invaginations of the limbal epithelium in the limbal zone are analogous to the limbal interpalisade crypts of human limbal tissue; and the presence of goblet cells in the bulbar conjunctiva. Immunohistochemistry analysis revealed that the epithelial progenitor markers cytokeratin (CK)15, CK14, p63α, and P-cadherin were expressed in both the limbal and conjunctival basal epithelium, whereas the basal cells of the limbal and conjunctival epithelium did not stain for CK3, CK12, E-cadherin, and CK13. Antibodies detecting marker proteins related to the extracellular matrix (collagen IV, Tenascin-C), cell-matrix adhesion (β-dystroglycan, integrin α3 and α6), mesenchymal cells (vimentin, CD90, CD44), neurons (neurofilament), immune cells (HLA-ABC; HLA-DR, CD1, CD4, CD14), vasculature (von Willebrand factor), and melanocytes (SRY-homeobox-10, human melanoma black-45, Tyrosinase) on the normal human ocular surface demonstrated similar immunoreactivity on the normal porcine ocular surface. Only a few antibodies (directed against N-cadherin, fibronectin, agrin, laminin α3 and α5, melan-A) appeared unreactive on porcine tissues. Our findings characterize the main immunohistochemical properties of the porcine ocular surface and provide a morphological and immunohistochemical basis useful to research using porcine models. Furthermore, the analyzed porcine ocular structures are similar to those of humans, confirming the potential usefulness of pig eyes to study ocular surface physiology and pathophysiology.

A temporal single cell transcriptome atlas of zebrafish anterior segment development

Anterior segment dysgenesis (ASD), resulting in vision impairment, stems from maldevelopment of anterior segment (AS) tissues. Incidence of ASD has been linked to malfunction of periocular mesenchyme cells (POM). POM cells specify into anterior segment mesenchyme (ASM) cells which colonize and produce AS tissues. In this study we uncover ASM developmental trajectories associated with formation of the AS. Using a transgenic line of zebrafish that fluorescently labels the ASM throughout development, Tg[foxc1b:GFP], we isolated GFP+ ASM cells at several developmental timepoints (48-144 hpf) and performed single cell RNA sequencing. Clustering analysis indicates subdifferentiation of ASM as early as 48 hpf and subsequent diversification into corneal epithelium/endothelium/stroma, or annular ligament (AL) lineages. Tracking individual clusters reveals common developmental pathways, up to 72 hpf, for the AL and corneal endothelium/stroma and distinct pathways for corneal epithelium starting at 48 hpf. Spatiotemporal validation of over 80 genes found associated with AS development demonstrates a high degree of conservation with mammalian trabecular meshwork and corneal tissues. In addition, we characterize thirteen novel genes associated with annular ligament and seven with corneal development. Overall, the data provide a molecular verification of the long-standing hypothesis that POM derived ASM give rise to AS tissues and highlight the high degree of conservation between zebrafish and mammals.

Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Mouse Corneal Epithelial Cells

Purpose

Corneal epithelial homeostasis is maintained by coordinated gene expression across distinct cell populations, but the gene regulatory programs underlying this cellular diversity remain to be characterized. Here we applied single-cell multi-omics analysis to delineate the gene regulatory profile of mouse corneal epithelial cells under normal homeostasis.

Methods

Single cells isolated from the cornea epithelium (with marginal conjunctiva) of adult mice were subjected to scRNA-seq and scATAC-seq using the 10×Genomics platform. Cell types were clustered by the graph-based visualization method uniform manifold approximation and projection and unbiased computational informatics analysis. The scRNA-seq and scATAC-seq datasets were integrated following the integration pipeline described in ArchR and Seurat.

Results

We characterized diverse corneal epithelial cell types based on gene expression signatures and chromatin accessibility. We found that cell type-specific accessibility regions were mainly located at distal regions, suggesting essential roles of distal regulatory elements in determining corneal epithelial cell diversity. Trajectory analyses revealed a continuum of cell state transition and higher coordination between transcription factor (TF) motif accessibility and gene expression during corneal epithelial cell differentiation. By integrating transcriptomic and chromatin accessibility analysis, we identified cell type-specific and shared gene regulation programs. We also uncovered critical TFs driving corneal epithelial cell differentiation, such as nuclear factor I (NFI) family members, Rarg, Elf3. We found that nuclear factor-κB (NF-κB) family members were positive TFs in limbal cells and some superficial cells, but they were involved in regulating distinct biological processes.

Conclusions

Our study presents a comprehensive gene regulatory landscape of mouse cornea epithelial cells, and provides valuable foundations for future investigation of corneal epithelial homeostasis in the context of cornea pathologies and regenerative medicine.

Cytoprotective Effects of Human Platelet Lysate during the Xeno-Free Culture of Human Donor Corneas

We evaluated the suitability of 2% human platelet lysate medium (2%HPL) as a replacement for 2% fetal bovine serum medium (2%FBS) for the xeno-free organ culture of human donor corneas. A total of 32 corneas from 16 human donors were cultured in 2%FBS for 3 days (TP1), then evaluated using phase contrast microscopy (endothelial cell density (ECD) and cell morphology). Following an additional 25-day culture period (TP2) in either 2%FBS or 2%HPL, the pairs were again compared using microscopy; then stroma and Descemet membrane/endothelium (DmE) were processed for next generation sequencing (NGS). At TP2 the ECD was higher in the 2%HPL group (2179 ± 288 cells/mm²) compared to 2%FBS (2113 ± 331 cells/mm²; p = 0.03), and endothelial cell loss was lower (ECL HPL = -0.7% vs. FBS = -3.8%; p = 0.01). There were no significant differences in cell morphology between TP1 and 2, or between 2%HPL and 2%FBS. NGS showed the differential expression of 1644 genes in endothelial cells and 217 genes in stromal cells. It was found that 2%HPL led to the upregulation of cytoprotective, anti-inflammatory and anti-fibrotic genes (HMOX1, SERPINE1, ANGPTL4, LEFTY2, GADD45B, PLIN2, PTX3, GFRA1/2), and the downregulation of pro-inflammatory/apoptotic genes (e.g., CXCL14, SIK1B, PLK5, PPP2R3B, FABP5, MAL, GATA3). 2%HPL is a suitable xeno-free substitution for 2%FBS in human cornea organ culture, inducing less ECL and producing potentially beneficial alterations in gene expression.

Dry eye disease in mice activates adaptive corneal epithelial regeneration distinct from constitutive renewal in homeostasis

Many epithelial compartments undergo constitutive renewal in homeostasis but activate unique regenerative responses following injury. The clear corneal epithelium is crucial for vision and is renewed from limbal stem cells (LSCs). Using single-cell RNA sequencing, we profiled the mouse corneal epithelium in homeostasis, aging, diabetes, and dry eye disease (DED), where tear deficiency predisposes the cornea to recurrent injury. In homeostasis, we capture the transcriptional states that accomplish continuous tissue turnover. We leverage our dataset to identify candidate genes and gene networks that characterize key stages across homeostatic renewal, including markers for LSCs. In aging and diabetes, there were only mild changes with <15 dysregulated genes. The constitutive cell types that accomplish homeostatic renewal were conserved in DED but were associated with activation of cell states that comprise "adaptive regeneration." We provide global markers that distinguish cell types in homeostatic renewal vs. adaptive regeneration and markers that specifically define DED-elicited proliferating and differentiating cell types. We validate that expression of SPARC, a marker of adaptive regeneration, is also induced in corneal epithelial wound healing and accelerates wound closure in a corneal epithelial cell scratch assay. Finally, we propose a classification system for LSC markers based on their expression fidelity in homeostasis and disease. This transcriptional dissection uncovers the dramatically altered transcriptional landscape of the corneal epithelium in DED, providing a framework and atlas for future study of these ocular surface stem cells in health and disease.

Corneal regeneration: insights in epithelial stem cell heterogeneity and dynamics

The discovery of slow-cycling cells at the corneal periphery three decades ago established the limbus as the putative corneal stem cell niche. Since then, studies have underscored the importance of the limbal stem cells in maintaining the health and function of the ocular surface. Advancements in our understanding of stem cell biology have been successfully translated into stem cell therapies for corneal diseases. Here, we review recent developments in mouse genetics, intravital imaging, and single-cell genomics that have revealed an underappreciated complexity of the limbal stem cells, from their molecular identity, function, and interactions with their niche environment. Continued efforts to elucidate stem cell dynamics of this extraordinary tissue are critical for not only understanding stem cell biology but also for advancing therapeutic innovation and development.

Transcriptomic Landscape and Functional Characterization of Human Induced Pluripotent Stem Cell-Derived Limbal Epithelial Progenitor Cells

Limbal stem cell deficiency (LSCD) is a complex, multifactorial disease affecting limbal epithelial progenitor cells (LEPC), which are essential for maintaining corneal stability and transparency. Human induced pluripotent stem cell-derived (hiPSC-) LEPC are a promising cell source for the treatment of LSCD. However, their similarity to native tissue-derived (T-) LEPC and their functional characterization has not been studied in detail. Here, we show that hiPSC-LEPC and T-LEPC have rather similar gene expression patterns, colony-forming ability, wound-healing capacity, and melanosome uptake. In addition, hiPSC-LEPC exhibited lower immunogenicity and reduced the proliferation of peripheral blood mononuclear cells compared with T-LEPC. Similarly, the hiPSC-LEPC secretome reduced the proliferation of vascular endothelial cells more than the T-LEPC secretome. Moreover, hiPSC-LEPC successfully repopulated decellularized human corneolimbal (DHC/L) scaffolds with multilayered epithelium, while basal deposition of fibrillary material was observed. These findings suggest that hiPSC-LEPC exhibited functional properties close to native LEPC and that hiPSC-LEPC-DHC/L scaffolds might be feasible for transplantation in patients suffering from LSCD in the future. Although hiPSC-LEPC-based stem cell therapy is promising, the current study also revealed new challenges, such as abnormal extracellular matrix deposition, that need to be overcome before hiPSC-LEPC-based stem cell therapies are viable.

Cellular and molecular profiles of larval and adult Xenopus corneal epithelia resolved at the single-cell level

Corneal Epithelial Stem Cells (CESCs) and their proliferative progeny, the Transit Amplifying Cells (TACs), are responsible for homeostasis and maintaining corneal transparency. Owing to our limited knowledge of cell fates and gene activity within the cornea, the search for unique markers to identify and isolate these cells remains crucial for ocular surface reconstruction. We performed single-cell RNA sequencing of corneal cells from larval and adult stages of Xenopus. Our results indicate that as the cornea develops and matures, there is an increase in cellular diversity, which is accompanied by a substantial shift in transcriptional profile, gene regulatory network and cell-cell communication dynamics. Our data also reveals several novel genes expressed in corneal cells and changes in gene expression during corneal differentiation at both developmental time-points. Importantly, we identify specific basal cell clusters in both the larval and adult cornea that comprise a relatively undifferentiated cell type and express distinct stem cell markers, which we propose are the putative larval and adult CESCs, respectively. This study offers a detailed atlas of single-cell transcriptomes in the frog cornea. In the future, this work will be useful to elucidate the function of novel genes in corneal epithelial homeostasis, wound healing and regeneration.

Single cell RNA-seq of human cornea organoids identifies cell fates of a developing immature cornea

The cornea is a protective and refractive barrier in the eye crucial for vision. Understanding the human cornea in health, disease, and cell-based treatments can be greatly advanced with cornea organoids developed in culture from induced pluripotent stem cells. While a limited number of studies have investigated the single-cell transcriptomic composition of the human cornea, its organoids have not been examined similarly. Here, we elucidated the transcriptomic cell fate map of 4-month-old human cornea organoids and human donor corneas. The organoids harbor cell clusters that resemble cells of the corneal epithelium, stroma, and endothelium, with subpopulations that capture signatures of early developmental states. Unlike the adult cornea where the largest cell population is stromal, the organoids contain large proportions of epithelial and endothelial-like cells. These corneal organoids offer a 3D model to study corneal diseases and integrated responses of different cell types.

Cell atlas of the human ocular anterior segment: Tissue-specific and shared cell types

The anterior segment of the eye consists of the cornea, iris, ciliary body, crystalline lens, and aqueous humor outflow pathways. Together, these tissues are essential for the proper functioning of the eye. Disorders of vision have been ascribed to defects in all of them; some disorders, including glaucoma and cataract, are among the most prevalent causes of blindness in the world. To characterize the cell types that compose these tissues, we generated an anterior segment cell atlas of the human eye using high-throughput single-nucleus RNA sequencing (snRNAseq). We profiled 195,248 nuclei from nondiseased anterior segment tissues of six human donors, identifying >60 cell types. Many of these cell types were discrete, whereas others, especially in the lens and cornea, formed continua corresponding to known developmental transitions that persist in adulthood. Having profiled each tissue separately, we performed an integrated analysis of the entire anterior segment, revealing that some cell types are unique to a single structure, whereas others are shared across tissues. The integrated cell atlas was then used to investigate cell type-specific expression patterns of more than 900 human ocular disease genes identified through either Mendelian inheritance patterns or genome-wide association studies.

Differentiation Trajectory of Limbal Stem and Progenitor Cells under Normal Homeostasis and upon Corneal Wounding

Limbal stem cells (LSCs) reside discretely at limbus surrounded by niche cells and progenitor cells. The aim of this study is to identify the heterogeneous cell populations at limbus under normal homeostasis and upon wounding using single-cell RNA sequencing in a mouse model. Two putative LSC types were identified which showed a differentiation trajectory into limbal progenitor cell (LPC) types under normal homeostasis and during wound healing. They were designated as “putative active LSCs” and “putative quiescent LSCs”, respectively, because the former type actively divided upon wounding while the later type stayed at a quiescent status upon wounding. The “putative quiescent LSCs” might contribute to a barrier function due to their characteristic markers regulating vascular and epithelial barrier and growth. Different types of LPCs at different proliferative statuses were identified in unwounded and wounded corneas with distinctive markers. Four maturation markers (Aldh3, Slurp1, Tkt, and Krt12) were screened out for corneal epithelium, which showed an increased expression along the differentiation trajectory during corneal epithelial maturation. In conclusion, our study identified two different types of putative LSCs and several types of putative LPCs under normal homeostasis and upon wounding, which will facilitate the understanding of corneal epithelial regeneration and wound healing.

Molecular Characteristics and Distribution of Adult Human Corneal Immune Cell Types

Background

The limbus is located at a 2-mm-wide area between the bulbar conjunctiva and the cornea and has been suggested to be the niche of corneal epithelial stem cells and immune cells. Like the skin and intestines, the cornea is also an important mucosal surface, and immune cells on the cornea play critical roles in immune surveillance to ensure barrier surface homeostasis and protection from various environmental damage and infections. Single-cell RNA sequencing (scRNA-seq) analysis of protein tyrosine phosphatase receptor type C positive (PTPRC⁺) hematopoietic cells from the corneal limbus could provide a single cell atlas of all the immune cell subsets.

Methods

We performed single-cell RNA sequencing to generate transcriptomic profile for 804 sort-purified hematopoietic cells from the corneal limbus of three healthy donors.

Results

Our analysis identified a primary transcriptomic pattern for multiple immune cell subtypes, including naive T cells, antiviral effector CD8⁺ T cells, and innate immune cells such as IDO1⁺ mature regulatory dendritic cells (mregDCs), macrophages, monocytes, and basophils in the human corneal limbus.

Conclusion

Overall, single-cell transcriptomic analysis of limbal immune cells suggested the possible contribution of these cells on the adaptive and innate immune response of the human cornea.