Corneal keratocytes, fibroblasts, and myofibroblasts exhibit distinct transcriptional profiles in vitro
After stromal injury to the cornea, the release of growth factors and pro-inflammatory cytokines promotes the activation of quiescent keratocytes into a migratory fibroblast and/or fibrotic myofibroblast phenotype. Persistence of the myofibroblast phenotype can lead to corneal fibrosis and scarring, which are leading causes of blindness worldwide. This study aims to establish comprehensive transcriptional profiles for cultured corneal keratocytes, fibroblasts, and myofibroblasts to gain insights into the mechanisms through which these phenotypic changes occur. Primary rabbit corneal keratocytes were cultured in either defined serum-free media (SF), fetal bovine serum (FBS) containing media, or in the presence of TGF-ß1 to induce keratocyte, fibroblast, or myofibroblast phenotypes, respectively. Bulk RNA sequencing followed by bioinformatic analyses was performed to identify significant differentially expressed genes (DEGs) and enriched biological pathways for each phenotype. Genes commonly associated with keratocytes, fibroblasts, or myofibroblasts showed high relative expression in SF, FBS, or TGF-ß1 culture conditions, respectively. Differential expression and functional analyses revealed novel DEGs for each cell type, as well as enriched pathways indicative of differences in proliferation, apoptosis, extracellular matrix (ECM) synthesis, cell-ECM interactions, cytokine signaling, and cell mechanics. Overall, these data demonstrate distinct transcriptional differences among cultured corneal keratocytes, fibroblasts, and myofibroblasts. We have identified genes and signaling pathways that may play important roles in keratocyte differentiation, including many related to mechanotransduction and ECM biology. Our findings have revealed novel molecular markers for each cell type, as well as possible targets for modulating cell behavior and promoting physiological corneal wound healing. Overall design: Comparative gene expression profiling analysis of bulk RNA-seq data from primary rabbit corneal keratocytes cultured in serum-free media, serum-containing media, or in the presence of TGF-ß1.