Putative Self-Organizing Human Corneal Organoids Recapitulate Human Corneal Architecture and Cellular Diversity

Background: Corneal organoids derived from pluripotent stem cells have emerged as powerful tools for studying corneal development, disease modeling, and regenerative medicine applications. While previous protocols have successfully generated corneal tissue structures, there remains a need for three-dimensional models that recapitulate the complex cellular architecture and diversity of native human cornea. Methods: We developed a modified spontaneous three-dimensional corneal organoid model using human embryonic stem cells (hESCs) through an adapted Self-formed Ectoderm Autonomous Multi-zone (SEAM) protocol. hESCs were cultured as spheroids in ultra-low-binding plates under normoxic conditions and differentiated over 7–8 weeks. Organoids were characterized using immunofluorescence staining for corneal-specific markers and single-cell RNA sequencing to assess cellular composition and gene expression patterns. Results: Approximately 20% of organoids developed transparent regions characteristic of corneal tissue by day 30 of differentiation. Immunofluorescence analysis revealed spatially organized expression of corneal markers, including ZO-1 and E-cadherin in the outermost epithelial layers, P63α-positive putative limbal stem cells at the epithelial–stromal interface, vimentin-positive stromal cells in the interior, and laminin-1 deposition that suggests Bowman’s membrane formation. The organoids expressed cornea-specific keratins (K3, K12, and K15) and the master regulator PAX6 in appropriate cellular compartments. Single-cell RNA sequencing identified 18 distinct cell clusters, including three corneal epithelium subclusters with differential expression of MUC16, KRT12, and ΔNp63α, two stromal populations with distinct inflammatory profiles, and a corneal endothelium cluster. Transcriptomic analysis confirmed expression of key corneal genes, including AQP3, CDH1, multiple keratins, mucins, and extracellular matrix components (HAS2, CD34, CD44, COL8A1, and KERA). Conclusions: This three-dimensional spheroid-based putative corneal organoid model successfully recapitulates the multilayered architecture and cellular diversity of human cornea, including stratified epithelium, putative limbal stem cells, stroma, and endothelium in spatially appropriate arrangements. The model demonstrates molecular signatures consistent with native corneal tissue and provides a valuable platform for studying corneal development, disease mechanisms, and potential therapeutic applications. Future optimization to improve organoid formation efficiency and functional maturation will enhance the utility of this system for both basic research and translational medicine.