Transcriptome analysis identifies the differentially expressed genes related to the stemness of limbal stem cells in mice

Limbal stem cells (LSCs) reside in the basal layer of limbal epithelial cells (LECs). They are crucial for maintenance of corneal epithelium homeostasis and corneal wound healing. Their stemness is determined by their gene expression pattern. Despite of several positive identifiers have been reported, the unique biomarker for LSCs still remain elusive. Differentially expressed genes (DEGs) between stem cells and differentiated cells affect the fate of stem cells via specific signaling pathway. In order to understand the DEGs in the LSCs, RNA-seq was firstly conducted using a mouse model. A total of 1907 up-regulated DEGs and 395 down-regulated DEGs were identified in the limbus (L) compared to central cornea (CC) and conjunctiva (Cj). Reliability of the expression of genes from RNA-seq analysis was evaluated by quantitative real-time polymerase chain reaction (qRT-PCR) and immunofluorescence staining. The expression pattern of putative biomarkers was considered to be age-related. In up-regulated DEGs GO analysis, 570 gene ontology (GO) terms were significantly enriched. Five groups of genes related with biological processes from these significantly enriched GO terms comprised ionic transport, regulation of tissue development, muscle contraction, visual perception, and cell adhesion, which were clustered as a weighted similar network. Whereas, in down-regulated DEGs GO analysis, 61 GO terms were significantly enriched and only one group of ATP biosynthesis and metabolic process were clustered. Furthermore, we identified 55 signaling pathways by the Kyoto Encyclopedia of Genes and Genomes (KEGG) database based on up-regulated genes and 14 KEGG pathways based on down-regulated genes. In this study, we provide a landscape of the expression of putative LSCs biomarkers and stemness-related signaling pathways in a mouse model. Our findings could aid in the identification of LSC niche factors that may be related to the stemness of the LSCs.