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

Deciphering the heterogeneity of differentiating hPSC-derived corneal limbal stem cells through single-cell RNA sequencing

A comprehensive understanding of the human pluripotent stem cell (hPSC) differentiation process stands as a prerequisite for the development of hPSC-based therapeutics. In this study, single-cell RNA sequencing (scRNA-seq) was performed to decipher the heterogeneity during differentiation of three hPSC lines toward corneal limbal stem cells (LSCs). The scRNA-seq data revealed nine clusters encompassing the entire differentiation process, among which five followed the anticipated differentiation path of LSCs. The remaining four clusters were previously undescribed cell states that were annotated as either mesodermal-like or undifferentiated subpopulations, and their prevalence was hPSC line dependent. Distinct cluster-specific marker genes identified in this study were confirmed by immunofluorescence analysis and employed to purify hPSC-derived LSCs, which effectively minimized the variation in the line-dependent differentiation efficiency. In summary, scRNA-seq offered molecular insights into the heterogeneity of hPSC-LSC differentiation, allowing a data-driven strategy for consistent and robust generation of LSCs, essential for future advancement toward clinical translation.

Impact of Combined Phototherapy and Melanocyte Transplantation on Indicators of Vitiligo Activity

OBJECTIVE

To investigate the effect of phototherapy combined with melanocyte transplantation on the activity index of vitiligo.

METHODS

One hundred twenty patients with stable vitiligo were selected and divided into 2 groups: phototherapy group (n = 60) and phototherapy combined with melanocyte therapy group (n = 60). Patients' vitiligo activity scores before and 6 months after treatment, patients' skin pigmentation responses 6 months after treatment, and patients' new Koebner cases 6 months after treatment were compared. The expression of tyrosinase and Melan-A in the skin samples was analyzed by immunohistochemistry.

RESULTS

The effect of skin surface repigmentation in the observation group was better than that in the control group (p < .05). The expression of tyrosinase and Melan-A in the observation group was higher than that in the control group (p < .05), indicating that the combined treatment could enhance the function of melanocytes. After 6 months of treatment, the incidence of the Koebner phenomenon in the observation group was lower than that in the control group (p < .05).

CONCLUSION

The combination of phototherapy and melanocyte transplantation can obviously improve the activity index of vitiligo, slow down the spread of white spots, reduce the formation of new white spots, and reduce the occurrence of the Koebner phenomenon.

Decoding cellular plasticity and niche regulation of limbal stem cells during corneal wound healing

BACKGROUND

Dysfunction or deficiency of corneal epithelium results in vision impairment or blindness in severe cases. The rapid and effective regeneration of corneal epithelial cells relies on the limbal stem cells (LSCs). However, the molecular and functional responses of LSCs and their niche cells to injury remain elusive.

METHODS

Single-cell RNA sequencing was performed on corneal tissues from normal mice and corneal epithelium defect models. Bioinformatics analysis was performed to confirm the distinct characteristics and cell fates of LSCs. Knockdown of Creb5 and OSM treatment experiment were performed to determine their roles of in corneal epithelial wound healing.

RESULTS

Our data defined the molecular signatures of LSCs and reconstructed the pseudotime trajectory of corneal epithelial cells. Gene network analyses characterized transcriptional landmarks that potentially regulate LSC dynamics, and identified a transcription factor Creb5, that was expressed in LSCs and significantly upregulated after injury. Loss-of-function experiments revealed that silencing Creb5 delayed the corneal epithelial healing and LSC mobilization. Through cell-cell communication analysis, we identified 609 candidate regeneration-associated ligand-receptor interaction pairs between LSCs and distinct niche cells, and discovered a unique subset of Arg1+ macrophages infiltrated after injury, which were present as the source of Oncostatin M (OSM), an IL-6 family cytokine, that were demonstrated to effectively accelerate the corneal epithelial wound healing.

CONCLUSIONS

This research provides a valuable single-cell resource and reference for the discovery of mechanisms and potential clinical interventions aimed at ocular surface reconstruction.

Exploring single-cell RNA sequencing as a decision-making tool in the clinical management of Fuchs' endothelial corneal dystrophy

Single-cell RNA sequencing (scRNA-seq) has enabled the identification of novel gene signatures and cell heterogeneity in numerous tissues and diseases. Here we review the use of this technology for Fuchs' Endothelial Corneal Dystrophy (FECD). FECD is the most common indication for corneal endothelial transplantation worldwide. FECD is challenging to manage because it is genetically heterogenous, can be autosomal dominant or sporadic, and progress at different rates. Single-cell RNA sequencing has enabled the discovery of several FECD subtypes, each with associated gene signatures, and cell heterogeneity. Current FECD treatments are mainly surgical, with various Rho kinase (ROCK) inhibitors used to promote endothelial cell metabolism and proliferation following surgery. A range of emerging therapies for FECD including cell therapies, gene therapies, tissue engineered scaffolds, and pharmaceuticals are in preclinical and clinical trials. Unlike conventional disease management methods based on clinical presentations and family history, targeting FECD using scRNA-seq based precision-medicine has the potential to pinpoint the disease subtypes, mechanisms, stages, severities, and help clinicians in making the best decision for surgeries and the applications of therapeutics. In this review, we first discuss the feasibility and potential of using scRNA-seq in clinical diagnostics for FECD, highlight advances from the latest clinical treatments and emerging therapies for FECD, integrate scRNA-seq results and clinical notes from our FECD patients and discuss the potential of applying alternative therapies to manage these cases clinically.

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.

Glycoprotein hormone subunit alpha 2 (GPHA2): A pituitary stem cell-expressed gene associated with NOTCH2 signaling

NOTCH2 is expressed in pituitary stem cells and is necessary for stem cell maintenance, proliferation, and differentiation. However, the pathways NOTCH2 engages to affect pituitary development remain unclear. In this study, we hypothesized that glycoprotein hormone subunit A2 (GPHA2), a corneal stem cell factor and ligand for the thyroid stimulating hormone receptor (TSHR), is downstream of NOTCH2 signaling. We found Gpha2 is expressed in quiescent pituitary stem cells by RNAscope in situ hybridization and scRNA seq. In Notch2 conditional knockout pituitaries, Gpha2 mRNA is reduced compared with control littermates. We then investigated the possible functions of GPHA2. Pituitaries treated with a GPHA2 peptide do not have a change in proliferation. However, in dissociated adult pituitary cells, GPHA2 increased pCREB expression and this induction was reversed by co-treatment with a TSHR inhibitor. These data suggest GPHA2 is a NOTCH2 related stem cell factor that activates TSHR signaling, potentially impacting pituitary development.

Sox9 marks limbal stem cells and is required for asymmetric cell fate switch in the corneal epithelium

Adult tissues with high cellular turnover require a balance between stem cell renewal and differentiation, yet the mechanisms underlying this equilibrium are unclear. The cornea exhibits a polarized lateral flow of progenitors from the peripheral stem cell niche to the center; attributed to differences in cellular fate. To identify genes that are critical for regulating the asymmetric fates of limbal stem cells and their transient amplified progeny in the central cornea, we utilized an in vivo cell cycle reporter to isolate proliferating basal cells across the anterior ocular surface epithelium and performed single-cell transcriptional analysis. This strategy greatly increased the resolution and revealed distinct basal cell identities with unique expression profiles of structural genes and transcription factors. We focused on Sox9; a transcription factor implicated in stem cell regulation across various organs. Sox9 was found to be differentially expressed between limbal stem cells and their progeny in the central corneal. Lineage tracing analysis confirmed that Sox9 marks long-lived limbal stem cells and conditional deletion led to abnormal differentiation and squamous metaplasia in the central cornea. These data suggest a requirement for Sox9 for the switch to asymmetric fate and commitment toward differentiation, as transient cells exit the limbal niche. By inhibiting terminal differentiation of corneal progenitors and forcing them into perpetual symmetric divisions, we replicated the Sox9 loss-of-function phenotype. Our findings reveal an essential role for Sox9 for the spatial regulation of asymmetric fate in the corneal epithelium that is required to sustain tissue homeostasis.

Integrative ATAC-seq and RNA-seq analysis of myogenic differentiation of ovine skeletal muscle satellite cell

Skeletal muscle satellite cells (SMSCs) play an important role in regulating muscle growth and regeneration. Chromatin accessibility allows physical interactions that synergistically regulate gene expression through enhancers, promoters, insulators, and chromatin binding factors. However, the chromatin accessibility altas and its regulatory role in ovine myoblast differentiation is still unclear. Therefore, ATAC-seq and RNA-seq analysis were performed on ovine SMSCs at the proliferation stage (SCG) and differentiation stage (SCD). 17,460 DARs (differential accessibility regions) and 3732 DEGs (differentially expressed genes) were identified. Based on joint analysis of ATAC-seq and RNA-seq, we revealed that PI3K-Akt, TGF-β and other signaling pathways regulated SMSCs differentiation. We identified two novel candidate genes, FZD5 and MAP2K6, which may affect the proliferation and differentiation of SMSCs. Our data identify potential cis regulatory elements of ovine SMSCs. This study can provide a reference for exploring the mechanisms of the differentiation and regeneration of SMSCs in the future.

Nanoceria-Mediated Cyclosporin A Delivery for Dry Eye Disease Management through Modulating Immune-Epithelial Crosstalk

Dry eye disease (DED) affects a substantial worldwide population with increasing frequency. Current single-targeting DED management is severely hindered by the existence of an oxidative stress-inflammation vicious cycle and complicated intercellular crosstalk within the ocular microenvironment. Here, a nanozyme-based eye drop, namely nanoceria loading cyclosporin A (Cs@P/CeO2), is developed, which possesses long-term antioxidative and anti-inflammatory capacities due to its regenerative antioxidative activity and sustained release of cyclosporin A (CsA). In vitro studies showed that the dual-functional Cs@P/CeO2 not only inhibits cellular reactive oxygen species production, sequentially maintaining mitochondrial integrity, but also downregulates inflammatory processes and repolarizes macrophages. Moreover, using flow cytometric and single-cell sequencing data, the in vivo therapeutic effect of Cs@P/CeO2 was systemically demonstrated, which rebalances the immune-epithelial communication in the corneal microenvironment with less inflammatory macrophage polarization, restrained oxidative stress, and enhanced epithelium regeneration. Collectively, our data proved that the antioxidative and anti-inflammatory Cs@P/CeO2 may provide therapeutic insights into DED management.

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.

The Potential Reversible Transition between Stem Cells and Transient-Amplifying Cells: The Limbal Epithelial Stem Cell Perspective

Stem cells (SCs) undergo asymmetric division, producing transit-amplifying cells (TACs) with increased proliferative potential that move into tissues and ultimately differentiate into a specialized cell type. Thus, TACs represent an intermediary state between stem cells and differentiated cells. In the cornea, a population of stem cells resides in the limbal region, named the limbal epithelial stem cells (LESCs). As LESCs proliferate, they generate TACs that move centripetally into the cornea and differentiate into corneal epithelial cells. Upon limbal injury, research suggests a population of progenitor-like cells that exists within the cornea can move centrifugally into the limbus, where they dedifferentiate into LESCs. Herein, we summarize recent advances made in understanding the mechanism that governs the differentiation of LESCs into TACs, and thereafter, into corneal epithelial cells. We also outline the evidence in support of the existence of progenitor-like cells in the cornea and whether TACs could represent a population of cells with progenitor-like capabilities within the cornea. Furthermore, to gain further insights into the dynamics of TACs in the cornea, we outline the most recent findings in other organ systems that support the hypothesis that TACs can dedifferentiate into SCs.

Endogenous TSG-6 modulates corneal inflammation following chemical injury

PURPOSE

Tumor necrosis factor (TNF)-stimulated gene-6 (TSG-6) is upregulated in various pathophysiological contexts, where it has a diverse repertoire of immunoregulatory functions. Herein, we investigated the expression and function of TSG-6 during corneal homeostasis and after injury.

METHODS

Human corneas, eyeballs from BALB/c (TSG-6+/+), TSG-6+/- and TSG-6-/- mice, human immortalized corneal epithelial cells and murine corneal epithelial progenitor cells were prepared for immunostaining and real time PCR analysis of endogenous expression of TSG-6. Mice were subjected to unilateral corneal debridement or alkali burn (AB) injuries and wound healing assessed over time using fluorescein stain, in vivo confocal microscopy and histology.

RESULTS

TSG-6 is endogenously expressed in the human and mouse cornea and established corneal epithelial cell lines and is upregulated after injury. A loss of TSG-6 has no structural and functional effect in the cornea during homeostasis. No differences were noted in the rate of corneal epithelial wound closure between BALB/c, TSG-6+/- and TSG-6-/- mice. TSG-6-/- mice presented decreased inflammatory response within the first 24 h of injury and accelerated corneal wound healing following AB when compared to control mice.

CONCLUSION

TSG-6 is endogenously expressed in the cornea and upregulated after injury where it propagates the inflammatory response following chemical injury.

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.

SinCMat: A single-cell-based method for predicting functional maturation transcription factors

A major goal of regenerative medicine is to generate tissue-specific mature and functional cells. However, current cell engineering protocols are still unable to systematically produce fully mature functional cells. While existing computational approaches aim at predicting transcription factors (TFs) for cell differentiation/reprogramming, no method currently exists that specifically considers functional cell maturation processes. To address this challenge, here, we develop SinCMat, a single-cell RNA sequencing (RNA-seq)-based computational method for predicting cell maturation TFs. Based on a model of cell maturation, SinCMat identifies pairs of identity TFs and signal-dependent TFs that co-target genes driving functional maturation. A large-scale application of SinCMat to the Mouse Cell Atlas and Tabula Sapiens accurately recapitulates known maturation TFs and predicts novel candidates. We expect SinCMat to be an important resource, complementary to preexisting computational methods, for studies aiming at producing functionally mature cells.

Analytical Workflows for Single-Cell Multiomic Data Using the BD Rhapsody Platform

The conversion of raw sequencing reads to biologically relevant data in high-throughput single-cell RNA sequencing experiments is a complex and involved process. Drawing meaning from thousands of individual cells to provide biological insight requires ensuring not only that the data are of the highest quality but also that the signal can be separated from noise. In this article, we describe a detailed analytical workflow, including six pipelines, that allows high-quality data analysis in single-cell multiomics. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Image analysis Basic Protocol 2: Sequencing quality control and generation of a gene expression matrix Basic Protocol 3: Gene expression matrix data pre-processing and analysis Basic Protocol 4: Advanced analysis Basic Protocol 5: Conversion to flow cytometry standard (FCS) format Basic Protocol 6: Visualization using graphical interfaces.

Human conjunctiva organoids to study ocular surface homeostasis and disease

The conjunctival epithelium covering the eye contains two main cell types: mucus-producing goblet cells and water-secreting keratinocytes, which present mucins on their apical surface. Here, we describe long-term expanding organoids and air-liquid interface representing mouse and human conjunctiva. A single-cell RNA expression atlas of primary and cultured human conjunctiva reveals that keratinocytes express multiple antimicrobial peptides and identifies conjunctival tuft cells. IL-4/-13 exposure increases goblet and tuft cell differentiation and drastically modifies the conjunctiva secretome. Human NGFR+ basal cells are identified as bipotent conjunctiva stem cells. Conjunctival cultures can be infected by herpes simplex virus 1 (HSV1), human adenovirus 8 (hAdV8), and SARS-CoV-2. HSV1 infection was reversed by acyclovir addition, whereas hAdV8 infection, which lacks an approved drug therapy, was inhibited by cidofovir. We document transcriptional programs induced by HSV1 and hAdV8. Finally, conjunctival organoids can be transplanted. Together, human conjunctiva organoid cultures enable the study of conjunctival (patho)-physiology.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Differential Efficacy of Small Molecules Dynasore and Mdivi-1 for the Treatment of Dry Eye Epitheliopathy or as a Countermeasure for Nitrogen Mustard Exposure of the Ocular Surface

The ocular surface comprises the wet mucosal epithelia of the cornea and conjunctiva, the associated glands, and the overlying tear film. Epitheliopathy is the common pathologic outcome when the ocular surface is subjected to oxidative stress. Whether different stresses act via the same or different mechanisms is not known. Dynasore and dyngo-4a, small molecules developed to inhibit the GTPase activity of classic dynamins DNM1, DNM2, and DNM3, but not mdivi-1, a specific inhibitor of DNM1L, protect corneal epithelial cells exposed to the oxidant tert-butyl hydroperoxide (tBHP). Here we report that, while dyngo-4a is the more potent inhibitor of endocytosis, dynasore is the better cytoprotectant. Dynasore also protects corneal epithelial cells against exposure to high salt in an in vitro model of dysfunctional tears in dry eye. We now validate this finding in vivo, demonstrating that dynasore protects against epitheliopathy in a mouse model of dry eye. Knockdown of classic dynamin DNM2 was also cytoprotective against tBHP exposure, suggesting that dynasore's effect is at least partially on target. Like tBHP and high salt, exposure of corneal epithelial cells to nitrogen mustard upregulated the unfolded protein response and inflammatory markers, but dynasore did not protect against nitrogen mustard exposure. In contrast, mdivi-1 was cytoprotective. Interestingly, mdivi-1 did not inhibit the nitrogen mustard-induced expression of inflammatory cytokines. We conclude that exposure to tBHP or nitrogen mustard, two different oxidative stress agents, cause corneal epitheliopathy via different pathologic pathways. SIGNIFICANCE STATEMENT: Results presented in this paper, for the first time, implicate the dynamin DNM2 in ocular surface epitheliopathy. The findings suggest that dynasore could serve as a new topical treatment for dry eye epitheliopathy and that mdivi-1 could serve as a medical countermeasure for epitheliopathy due to nitrogen mustard exposure, with potentially increased efficacy when combined with anti-inflammatory agents and/or UPR modulators.

Transcriptional profiling specifies the pathogen-specific human host response to infectious keratitis

Purpose

Corneal infections are a leading cause of visual impairment and blindness worldwide. Here we applied high-resolution transcriptomic profiling to assess the general and pathogen-specific molecular and cellular mechanisms during human corneal infection.

Methods

Clinical diagnoses of herpes simplex virus (HSV) (n=5) and bacterial/fungal (n=5) keratitis were confirmed by histology. Healthy corneas (n=7) and keratoconus (n=4) samples served as controls. Formalin-fixed, paraffin-embedded (FFPE) human corneal specimens were analyzed using the 3' RNA sequencing method Massive Analysis of cDNA Ends (MACE RNA-seq). The cellular host response was investigated using comprehensive bioinformatic deconvolution (xCell and CYBERSORTx) analyses and by integration with published single cell RNA-seq data of the human cornea.

Results

Our analysis identified 216 and 561 genes, that were specifically overexpressed in viral or bacterial/fungal keratitis, respectively, and allowed to distinguish the two etiologies. The virus-specific host response was driven by adaptive immunity and associated molecular signaling pathways, whereas the bacterial/fungal-specific host response mainly involved innate immunity signaling pathways and cell types. We identified several genes and pathways involved in the host response to infectious keratitis, including CXCL9, CXCR3, and MMP9 for viral, and S100A8/A9, MMP9, and the IL17 pathway for bacterial/fungal keratitis.

Conclusions

High-resolution molecular profiling provides new insights into the human corneal host response to viral and bacterial/fungal infection. Pathogen-specific molecular profiles may provide the foundation for novel diagnostic biomarker and therapeutic approaches that target inflammation-induced damage to corneal host cells with the goal to improve the outcome of infectious keratitis.

The single-cell transcriptomic atlas and RORA-mediated 3D epigenomic remodeling in driving corneal epithelial differentiation

Proper differentiation of corneal epithelial cells (CECs) from limbal stem/progenitor cells (LSCs) is required for maintenance of ocular homeostasis and clear vision. Here, using a single-cell transcriptomic atlas, we delineate the comprehensive and refined molecular regulatory dynamics during human CEC development and differentiation. We find that RORA is a CEC-specific molecular switch that initiates and drives LSCs to differentiate into mature CECs by activating PITX1. RORA dictates CEC differentiation by establishing CEC-specific enhancers and chromatin interactions between CEC gene promoters and distal regulatory elements. Conversely, RORA silences LSC-specific promoters and disrupts promoter-anchored chromatin loops to turn off LSC genes. Collectively, our work provides detailed and comprehensive insights into the transcriptional dynamics and RORA-mediated epigenetic remodeling underlying human corneal epithelial differentiation.

Single-cell transcriptomic analysis of corneal organoids during development

Corneal organoids are useful tools for disease modeling and tissue transplantation; however, they have not yet been well studied during maturation. We characterized human iPSC-derived corneal organoids at 1, 2, 3, and 4 months of development using single-cell RNA sequencing to determine the cellular heterogeneity at each stage. We found pluripotent cell clusters committed to epithelial cell lineage at 1 month; early corneal epithelial, endothelial, and stromal cell markers at 2 months; keratocytes as the largest cell population at 3 months; and a large epithelial cell population at 4 months. We compared organoid to fetal corneal development at different stages and found that 4-month organoids closely resemble the corneal cellular complexity of the fetal (16 post conception week) and adult cornea. Using RNA velocity trajectory analysis, we found that less differentiated cells appear to give rise to corneal epithelial cells during development.

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.

Investigation of Potential Crucial Genes and Key Pathways in Keratoconus: An Analysis of Gene Expression Omnibus Data

Keratoconus is one of the most common causes leading to visual impairment in young adult population. The pathogenesis of keratoconus remains poorly understood. The aim of this study was to identify the potential key genes and pathways associated with keratoconus and to further analyze its molecular mechanism. Two RNA-sequencing datasets of keratoconus and paired normal corneal tissues from the Gene Expression Omnibus database were obtained. Differentially expressed genes (DEGs) were identified, and the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted. The protein-protein interaction (PPI) network of the DEGs was established, and the hub genes and significant gene modules of PPI were further constructed. Lastly, the GO and KEGG analyses of the hub gene were performed. In total, 548 common DEGs were identified. GO enrichment analysis showed that the DEGs were primarily associated with regulation of cell adhesion, the response to molecule of bacterial origin, lipopolysaccharide and biotic stimulus, collagen-containing extracellular matrix, extracellular matrix, and structure organization. KEGG pathway analysis showed that these DEGs were mainly involved in the TNF signaling pathway, IL-17 signaling pathway, Rheumatoid arthritis, Cytokine-cytokine receptor interaction. The PPI network was constructed with 146 nodes and 276 edges, and 3 significant modules are selected. Finally, top 10 hub genes were identified from the PPI network. The results revealed that extracellular matrix remodeling and immune inflammatory response could be the key links of keratoconus, TNF, IL6, IL1A, IL1B, CCL3, MMP3, MMP9, MMP1, and TGFB1 may be potential crucial genes, and TNF signaling pathway and IL-17 signaling pathway were the potential pathways accounting for pathogenesis and development of keratoconus.

The cellular responses of corneal fibroblasts to cyclic stretching loads

Mechanical signaling plays a crucial role in maintaining extracellular matrix (ECM) homeostasis in various structures. In this study, we investigated the responses of corneal fibroblasts to cyclic stretching loads using an in vitro cell culture system. Bovine corneal fibroblasts were cultured and subjected to equibiaxial cyclic strain of 15% for 72 h at a frequency of 0.25 Hz, with bovine skin fibroblasts used as a comparison. We explored various cellular behaviors, including morphological changes, cell proliferation, and metabolism in response to mechanical stretching loads. The expression of genes, protein secretion, and enzymatic activity for several major metalloproteinases was also determined through Q-PCR, Western blot, and gel zymography. Additionally, we investigated the involvement of mitogen-activated protein kinases (MAPKs) signaling pathways in the corneal fibroblasts when subjected to mechanical stimuli. Our findings revealed that, compared to skin fibroblasts, corneal fibroblasts were reluctant to morphological changes in response to a prolonged (72 h) and high-amplitude (15% of strain) cyclic stretching load. However, cyclic stretching loads stimulated the upregulation of MMP-2 expression in corneal fibroblasts via the MAPK signaling pathways involving extracellular signal-regulated kinase and p38. Together with a lack of upregulation in type I collagen expression, our results indicate the induction of the ECM degradation process in corneal fibroblasts in response to cyclic stretching. These findings emphasize the mechanoresponsive nature of corneal fibroblasts and shed light on the potential impact of intense mechanical stress on the cornea in both normal and pathological conditions such as keratoconus, providing valuable insights for understanding corneal mechanobiology.

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.

Toward Corneal Limbus In Vitro Model: Regulation of hPSC-LSC Phenotype by Matrix Stiffness and Topography During Cell Differentiation Process

A functional limbal epithelial stem cells (LSC) niche is a vital element in the regular renewal of the corneal epithelium by LSCs and maintenance of good vision. However, little is known about its unique structure and mechanical properties on LSC regulation, creating a significant gap in development of LSC-based therapies. Herein, the effect of mechanical and architectural elements of the niche on human pluripotent derived LSCs (hPSC-LSC) phenotype and growth is investigated in vitro. Specifically, three formulations of polyacrylamide gels with different controlled stiffnesses are used for culture and characterization of hPSC-LSCs from different stages of differentiation. In addition, limbal mimicking topography in polydimethylsiloxane is utilized for culturing hPSC-LSCs at early time point of differentiation. For comparison, the expression of selected key proteins of the corneal cells is analyzed in their native environment through whole mount staining of human donor corneas. The results suggest that mechanical response and substrate preference of the cells is highly dependent on their developmental stage. In addition, data indicate that cells may carry possible mechanical memory from previous culture matrix, both highlighting the importance of mechanical design of a functional in vitro limbus model.

Liquid-biopsy proteomics combined with AI identifies cellular drivers of eye aging and disease in vivo

Single-cell analysis in living humans is essential for understanding disease mechanisms, but it is impractical in non-regenerative organs, such as the eye and brain, because tissue biopsies would cause serious damage. We resolve this problem by integrating proteomics of liquid biopsies with single-cell transcriptomics from all known ocular cell types to trace the cellular origin of 5,953 proteins detected in the aqueous humor. We identified hundreds of cell-specific protein markers, including for individual retinal cell types. Surprisingly, our results reveal that retinal degeneration occurs in Parkinson's disease, and the cells driving diabetic retinopathy switch with disease stage. Finally, we developed artificial intelligence (AI) models to assess individual cellular aging and found that many eye diseases not associated with chronological age undergo accelerated molecular aging of disease-specific cell types. Our approach, which can be applied to other organ systems, has the potential to transform molecular diagnostics and prognostics while uncovering new cellular disease and aging mechanisms.

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.

Characterization of the Limbal Epithelial Stem Cell Niche

Purpose

Limbal epithelial stem cells (LESCs) reside within a LSC niche (LSCN). We recently identified that hyaluronan (HA) is a major constituent of the LSCN, and that HA is necessary for maintaining LESCs in the "stem cell" state, both in vitro and in vivo. Herein, we characterized the LSCN to identify key components of the HA-specific LSCN.

Methods

The cornea and limbal rim were dissected from mouse corneas, subjected to mRNA extraction, and sequenced using a NextSeq 500 (Illumina) and data processed using CLC Genomics Workbench 20 (Qiagen) and the STRING database to identify key components of the LSCN. Their expression was confirmed by real-time PCR, Western blotting, and immunohistochemistry. Furthermore, the differential expression of key compounds in different corneal cell types were determined with single-cell RNA sequencing.

Results

We identified that the hyaladherins inter-alpha-inhibitor (IαI), TSG-6 and versican are highly expressed in the limbus. Specifically, HA/HC complexes are present in the LSCN, in the stroma underlying the limbal epithelium, and surrounding the limbal vasculature. For IαI, heavy chains 5 and 2 (HC5 and HC2) were found to be the most highly expressed HCs in the mouse and human limbus and were associate with HA-forming HA/HC-specific matrices.

Conclusions

The LSCN contains HA/HC complexes, which have been previously correlated with stem cell niches. The identification of HA/HC complexes in the LSCN could serve as a new therapeutic avenue for treating corneal pathology. Additionally, HA/HC complexes could be used as a substrate for culturing LESCs before LESC transplantation.

Identification of the regulatory circuit governing corneal epithelial fate determination and disease

The transparent corneal epithelium in the eye is maintained through the homeostasis regulated by limbal stem cells (LSCs), while the nontransparent epidermis relies on epidermal keratinocytes for renewal. Despite their cellular similarities, the precise cell fates of these two types of epithelial stem cells, which give rise to functionally distinct epithelia, remain unknown. We performed a multi-omics analysis of human LSCs from the cornea and keratinocytes from the epidermis and characterized their molecular signatures, highlighting their similarities and differences. Through gene regulatory network analyses, we identified shared and cell type-specific transcription factors (TFs) that define specific cell fates and established their regulatory hierarchy. Single-cell RNA-seq (scRNA-seq) analyses of the cornea and the epidermis confirmed these shared and cell type-specific TFs. Notably, the shared and LSC-specific TFs can cooperatively target genes associated with corneal opacity. Importantly, we discovered that FOSL2, a direct PAX6 target gene, is a novel candidate associated with corneal opacity, and it regulates genes implicated in corneal diseases. By characterizing molecular signatures, our study unveils the regulatory circuitry governing the LSC fate and its association with corneal opacity.

Elucidating the Corneal Endothelial Cell Proliferation Capacity through an Interspecies Transcriptome Comparison

The regenerative capacity of corneal endothelial cells (CECs) differs between species; in bigger mammals, CECs are arrested in a non-proliferative state. Damage to these cells can compromise their function causing corneal opacity. Corneal transplantation is the current treatment for the recovery of clear eyesight, but the donor tissue demand is higher than the availability and there is a need to develop novel treatments. Interestingly, rabbit CECs retain a high proliferative profile and can repopulate the endothelium. There is a lack of fundamental knowledge to explain these differences. Gaining information on their transcriptomic variances could allow the identification of CEC proliferation drivers. In this study, human, sheep, and rabbit CECs are analyzed at the transcriptomic level. To understand the differences across each species, a pipeline for the analysis of pathways with different activities is generated. The results reveal that 52 pathways have different activity when comparing species with non-proliferative CECs (human and sheep) to species with proliferative CECs (rabbit). The results show that Notch and TGF-β pathways have increased activity in species with non-proliferative CECs, which might be associated with their low proliferation. Overall, this study illustrates transcriptomic pathway-level differences that can provide leads to develop novel therapies to regenerate the corneal endothelium.

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.

Characterization of Cell Surface Glycan Profiles in Human and Mouse Corneas Using Lectin Microarrays

The advent of high-throughput sequencing technologies has facilitated the profiling of glycosylation genes at a single-cell level in complex biological systems, but the significance of these gene signatures to the composition of the glycocalyx remains ambiguous. Here, we used lectin microarrays to characterize the composition of cell surface glycans in human and mouse corneas and determine its relationship to single-cell transcriptomic data. Our results identify a series of cell surface glycan signatures that are unique to the different cell types of the human cornea and that correlate, to a certain extent, with the transcriptional expression of glycosylation genes. These include pathways involved in the biosynthesis of O-glycans in epithelial cells and core fucose on stromal and endothelial cell surfaces. Moreover, we show that human and mouse corneas display some structural differences in terms of cell surface glycan composition. These results could provide insights into the specialized function of individual cell types in the cornea and foster the identification of novel cornea-specific biomarkers.

Wnt activation as a potential therapeutic approach to treat partial limbal stem cell deficiency

Limbal epithelial stem/progenitor cells (LSCs) are adult stem cells located at the limbus, tightly regulated by their niche involving numerous signaling pathways, such as Wnt. Wnt proteins are secreted morphogens that play critical roles in embryonic development, stem cell proliferation, self-renewal, tissue regeneration, and remodeling in adults. It has been shown that a small molecule Wnt mimic could improve LSCs expansion ex vivo. Damage to the LSCs and/or their niche can lead to limbal stem cell deficiency (LSCD), a condition that can cause corneal blindness and is difficult to treat. This study explored if repopulating residual LSCs in partial LSCD through Wnt activation could be a novel therapeutic approach. To mimic LSCD due to a chemical injury, single cultured LSCs were exposed to various concentrations of sodium hydroxide. A progressive loss of the LSCs phenotype was observed: the percentage of p63bright cells and cytokeratin (K)14+ cells decreased while the percentage of K12+ increased. Wnt activation was attained by treating the LSCs with lithium chloride (LiCl) and a small-molecule Wnt mimic, respectively. After 18 h of treatment, LSCs proliferation was increased, and the LSCs phenotype was recovered, while the untreated cells did not proliferate and lost their phenotype. The percentage of p63bright cells was significantly higher in the Wnt mimic-treated cells compared with untreated cells, while the percentage of K12+ cells was significantly lower. These findings suggest that local Wnt activation may rescue LSCs upon alkaline injury.

Artificial Intelligence Models for Cell Type and Subtype Identification Based on Single-Cell RNA Sequencing Data in Vision Science

Single-cell RNA sequencing (scRNA-seq) provides a high throughput, quantitative and unbiased framework for scientists in many research fields to identify and characterize cell types within heterogeneous cell populations from various tissues. However, scRNA-seq based identification of discrete cell-types is still labor intensive and depends on prior molecular knowledge. Artificial intelligence has provided faster, more accurate, and user-friendly approaches for cell-type identification. In this review, we discuss recent advances in cell-type identification methods using artificial intelligence techniques based on single-cell and single-nucleus RNA sequencing data in vision science. The main purpose of this review paper is to assist vision scientists not only to select suitable datasets for their problems, but also to be aware of the appropriate computational tools to perform their analysis. Developing novel methods for scRNA-seq data analysis remains to be addressed in future studies.

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.

Differential roles of FOXC2 in the trabecular meshwork and Schlemm's canal in glaucomatous pathology

Impaired development and maintenance of Schlemm's canal (SC) are associated with perturbed aqueous humor outflow and intraocular pressure. The angiopoietin (ANGPT)/TIE2 signaling pathway regulates SC development and maintenance, whereas the molecular mechanisms of crosstalk between SC and the neural crest (NC)-derived neighboring tissue, the trabecular meshwork (TM), are poorly understood. Here, we show NC-specific forkhead box (Fox)c2 deletion in mice results in impaired SC morphogenesis, loss of SC identity, and elevated intraocular pressure. Visible-light optical coherence tomography analysis further demonstrated functional impairment of the SC in response to changes in intraocular pressure in NC-Foxc2 -/- mice, suggesting altered TM biomechanics. Single-cell RNA-sequencing analysis identified that this phenotype is predominately characterized by transcriptional changes associated with extracellular matrix organization and stiffness in TM cell clusters, including increased matrix metalloproteinase expression, which can cleave the TIE2 ectodomain to produce soluble TIE2. Moreover, endothelial-specific Foxc2 deletion impaired SC morphogenesis because of reduced TIE2 expression, which was rescued by deleting the TIE2 phosphatase VE-PTP. Thus, Foxc2 is critical in maintaining SC identity and morphogenesis via TM-SC crosstalk.

Posterior corneoscleral limbus: Architecture, stem cells, and clinical implications

The limbus is a transition from the cornea to conjunctiva and sclera. In human eyes, this thin strip has a rich variation of tissue structures and composition, typifying a change from scleral irregularity and opacity to corneal regularity and transparency; a variation from richly vascularized conjunctiva and sclera to avascular cornea; the neural passage and drainage of aqueous humor. The limbal stroma is enriched with circular fibres running parallel to the corneal circumference, giving its unique role in absorbing small pressure changes to maintain corneal curvature and refractivity. It contains specific niches housing different types of stem cells for the corneal epithelium, stromal keratocytes, corneal endothelium, and trabecular meshwork. This truly reflects the important roles of the limbus in ocular physiology, and the limbal functionality is crucial for corneal health and the entire visual system. Since the anterior limbus containing epithelial structures and limbal epithelial stem cells has been extensively reviewed, this article is focused on the posterior limbus. We have discussed the structural organization and cellular components of the region beneath the limbal epithelium, the characteristics of stem cell types: namely corneal stromal stem cells, endothelial progenitors and trabecular meshwork stem cells, and recent advances leading to the emergence of potential cell therapy options to replenish their respective mature cell types and to correct defects causing corneal abnormalities. We have reviewed different clinical disorders associated with defects of the posterior limbus and summarized the available preclinical and clinical evidence about the developing topic of cell-based therapy for corneal disorders.

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.

Cultivated autologous limbal epithelial cell (CALEC) transplantation: Development of manufacturing process and clinical evaluation of feasibility and safety

To treat unilateral limbal stem cell (LSC) deficiency, we developed cultivated autologous limbal epithelial cells (CALEC) using an innovative xenobiotic-free, serum-free, antibiotic-free, two-step manufacturing process for LSC isolation and expansion onto human amniotic membrane with rigorous quality control in a good manufacturing practices facility. Limbal biopsies were used to generate CALEC constructs, and final grafts were evaluated by noninvasive scanning microscopy and tested for viability and sterility. Cultivated cells maintained epithelial cell phenotype with colony-forming and proliferative capacities. Analysis of LSC biomarkers showed preservation of "stemness." After preclinical development, a phase 1 clinical trial enrolled five patients with unilateral LSC deficiency. Four of these patients received CALEC transplants, establishing preliminary feasibility. Clinical case histories are reported, with no primary safety events. On the basis of these results, a second recruitment phase of the trial was opened to provide longer term safety and efficacy data on more patients.

Colocalization of corneal resistance factor GWAS loci with GTEx e/sQTLs highlights plausible candidate causal genes for keratoconus postnatal corneal stroma weakening

Background: Genome-wide association studies (GWAS) for corneal resistance factor (CRF) have identified 100s of loci and proved useful to uncover genetic determinants for keratoconus, a corneal ectasia of early-adulthood onset and common indication of corneal transplantation. In the current absence of studies to probe the impact of candidate causal variants in the cornea, we aimed to fill some of this knowledge gap by leveraging tissue-shared genetic effects. Methods: 181 CRF signals were examined for evidence of colocalization with genetic signals affecting steady-state gene transcription and splicing in adult, non-eye, tissues of the Genotype-Tissue Expression (GTEx) project. Expression of candidate causal genes thus nominated was evaluated in single cell transcriptomes from adult cornea, limbus and conjunctiva. Fine-mapping and colocalization of CRF and keratoconus GWAS signals was also deployed to support their sharing causal variants. Results and discussion: 26.5% of CRF causal signals colocalized with GTEx v8 signals and nominated genes enriched in genes with high and specific expression in corneal stromal cells amongst tissues examined. Enrichment analyses carried out with nearest genes to all 181 CRF GWAS signals indicated that stromal cells of the limbus could be susceptible to signals that did not colocalize with GTEx's. These cells might not be well represented in GTEx and/or the genetic associations might have context specific effects. The causal signals shared with GTEx provide new insights into mediation of CRF genetic effects, including modulation of splicing events. Functionally relevant roles for several implicated genes' products in providing tensile strength, mechano-sensing and signaling make the corresponding genes and regulatory variants prime candidates to be validated and their roles and effects across tissues elucidated. Colocalization of CRF and keratoconus GWAS signals strengthened support for shared causal variants but also highlighted many ways into which likely true shared signals could be missed when using readily available GWAS summary statistics.

Identification of BST2 as a conjunctival epithelial stem/progenitor cell marker

The conjunctival epithelium consists of conjunctival epithelial cells and goblet cells derived from conjunctival epithelial stem/progenitor cells. However, the source of these cells is not well known because no specific markers for conjunctival epithelial stem/progenitor cells have been discovered. Therefore, to identify conjunctival epithelial stem/progenitor cell markers, we performed single-cell RNA sequencing of a conjunctival epithelial cell population derived from human-induced pluripotent stem cells (hiPSCs). The following conjunctival epithelial markers were identified: BST2, SLC2A3, AGR2, TMEM54, OLR1, and TRIM29. Notably, BST2 was strongly positive in the basal conjunctival epithelium, which is thought to be rich in stem/progenitor cells. Moreover, BST2 was able to sort conjunctival epithelial stem/progenitor cells from hiPSC-derived ocular surface epithelial cell populations. BST2-positive cells were highly proliferative and capable of successfully generating conjunctival epithelial sheets containing goblet cells. In conclusion, BST2 has been identified as a specific marker of conjunctival epithelial stem/progenitor cells.

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.

Sequence variants contributing to dysregulated inflammatory responses across keratoconic cone surface in adolescent patients with keratoconus

Background

Keratoconus (KTCN) is the most common corneal ectasia resulting in a conical shape of the cornea. Here, genomic variation in the corneal epithelium (CE) across the keratoconic cone surface in patients with KTCN and its relevance in the functioning of the immune system were assessed.

Methods

Samples from four unrelated adolescent patients with KTCN and two control individuals were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated towards the whole-genome sequencing (WGS) study embracing a total of 18 experimental samples. The coding and non-coding sequence variation, including structural variation, was assessed and then evaluated together with the previously reported transcriptomic outcomes for the same CE samples and full-thickness corneas.

Results

First, pathway enrichment analysis of genes with identified coding variants pointed to "Antigen presentation" and "Interferon alpha/beta signaling" as the most overrepresented pathways, indicating the involvement of inflammatory responses in KTCN. Both coding and non-coding sequence variants were found in genes (or in their close proximity) linked to the previously revealed KTCN-specific cellular components, namely, "Actin cytoskeleton", "Extracellular matrix", "Collagen-containing extracellular matrix", "Focal adhesion", "Hippo signaling pathway", and "Wnt signaling" pathways. No genomic heterogeneity across the corneal surface was found comparing the assessed topographic regions. Thirty-five chromosomal regions enriched in both coding and non-coding KTCN-specific sequence variants were revealed, with a most representative 5q locus previously recognized as involved in KTCN.

Conclusion

The identified genomic features indicate the involvement of innate and adaptive immune system responses in KTCN pathogenesis.

Corneal epithelial development and homeostasis

The corneal epithelium (CE), the most anterior cellular structure of the eye, is a self-renewing stratified squamous tissue that protects the rest of the eye from external elements. Each cell in this exquisite three-dimensional structure needs to have proper polarity and positional awareness for the CE to serve as a transparent, refractive, and protective tissue. Recent studies have begun to elucidate the molecular and cellular events involved in the embryonic development, post-natal maturation, and homeostasis of the CE, and how they are regulated by a well-coordinated network of transcription factors. This review summarizes the status of related knowledge and aims to provide insight into the pathophysiology of disorders caused by disruption of CE development, and/or homeostasis.

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.

Consensus label propagation with graph convolutional networks for single-cell RNA sequencing cell type annotation

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) data, annotated by cell type, is useful in a variety of downstream biological applications, such as profiling gene expression at the single-cell level. However, manually assigning these annotations with known marker genes is both time-consuming and subjective.

RESULTS

We present a Graph Convolutional Network (GCN)-based approach to automate the annotation process. Our process builds upon existing labeling approaches, using state-of-the-art tools to find cells with highly confident label assignments through consensus and spreading these confident labels with a semi-supervised GCN. Using simulated data and two scRNA-seq datasets from different tissues, we show that our method improves accuracy over a simple consensus algorithm and the average of the underlying tools. We also compare our method to a nonparametric neighbor majority approach, showing comparable results. We then demonstrate that our GCN method allows for feature interpretation, identifying important genes for cell type classification. We present our completed pipeline, written in PyTorch, as an end-to-end tool for automating and interpreting the classification of scRNA-seq data.

AVAILABILITY AND IMPLEMENTATION

Our code for conducting the experiments in this paper and using our model is available at https://github.com/lewinsohndp/scSHARP.

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.

Dendrimers in Corneal Drug Delivery: Recent Developments and Translational Opportunities

Dendrimers are biocompatible organic nanomaterials with unique physicochemical properties, making them the focus of recent research in drug delivery. The cornea of the human eye presents a challenge for drug transit due to its inherently impenetrable nature, requiring nanocarrier-mediated targeted drug delivery. This review intends to examine recent advancements in the use of dendrimers for corneal drug delivery, including their properties and their potential for treating various ocular diseases. The review will also highlight the benefit of the novel technologies that have been developed and applied in the field, such as corneal targeting, drug release kinetics, treatments for dry eye disease, antibacterial drug delivery, corneal inflammation, and corneal tissue engineering. The review seeks to provide a comprehensive overview of the current state of research in this field, along with the translational developments in the field of dendrimer-based therapeutics and imaging agents and inspire the potential for future developments and translational opportunities in dendrimers based corneal drug delivery.