Abstract: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.

Gene Expression Profile of Extracellular Matrix and Adhesion Molecules in the Human Normal Corneal Stroma

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

Proteins of the Corneal Stroma: Importance in Visual Function

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

Corneal keratocytes, fibroblasts, and myofibroblasts exhibit distinct transcriptional profiles in vitro

Patterned expression of neurotrophic factors and receptors in human limbal and corneal regions

Exploration of the Pathophysiology of High Myopia Via Proteomic Profiling of Human Corneal Stromal Lenticules.

RNA-Seq Analysis Unraveling Novel Genes and Pathways Influencing Corneal Wound Healing

Differential expression of genes in cells cultured from juxtacanalicular trabecular meshwork and Schlemm's canal

Transcriptomic landscape of quiescent and proliferating human corneal stromal fibroblasts

Fibroblastic and bone marrow-derived cellularity in the corneal stroma

Differential epithelial and stromal protein profiles in cone and non-cone regions of keratoconus corneas

Molecular Identity of Human Limbal Heterogeneity Involved in Corneal Homeostasis and Privilege.

Comprehensive Transcriptome Analysis of Patients With Keratoconus Highlights the Regulation of Immune Responses and Inflammatory Processes

Expression Profiling in Glaucomatous Human Lamina Cribrosa Cells Based on Graph-Clustering Approach.

Proteomic Characterization of Corneal Epithelial and Stromal Cell-Derived Extracellular Vesicles

Proteomic profiling of human corneal stroma from long-term contact lens wearers reveals activation of inflammatory responses

A Novel Combined Oxidative Stress and Extracellular Matrix Related Predictive Gene Signature for Keratoconus

Preferential Gene Expression in the Limbus of the Vervet Monkey.

Expression of matrix metalloproteinases and their inhibitors in corneal stromal fibroblasts and keratocytes from healthy and keratoconus corneas

Differential expression of epithelial basement membrane components nidogens and perlecan in corneal stromal cells in vitro