Single cell transcriptomic analysis reveals cellular diversity of murine esophageal epithelium and age-associated mitochondrial dysfunction

Stratified squamous epithelium of the esophagus is comprised of basal keratinocytes that execute a terminal differentiation program in overlying suprabasal and superficial cell layers. Although morphologic progression coupled with expression of specific molecular markers has been characterized along the esophageal epithelial differentiation gradient, the molecular heterogeneity within the cell types along this trajectory has yet to be classified at the level of single cell resolution. To explore the molecular characteristics of esophageal keratinocytes along the squamous differentiation continuum, we performed single cell RNA-Sequencing transcriptomic profiling of 7,972 cells from murine esophageal epithelial sheets. We identified 8 distinct cell clusters in esophageal epithelium, unveiling an unexpected level of diversity, particularly among basal cells. We further mapped the cellular pathways and lineage trajectories within basal, suprabasal, and superficial clusters as well as within the heterogeneous basal cell populations, providing a comprehensive molecular view of esophageal epithelial cells in the context of squamous differentiation. Finally, we explored the impact of tissue aging upon esophageal epithelial cell clusters and demonstrated that mitochondrial dysfunction is a feature of aging in normal esophageal epithelium. These studies provide an unparalleled molecular perspective on murine esophageal keratinocytes that will serve as a valuable resource for dissecting cell type-specific roles in esophageal biology under conditions of homeostasis, aging, and tissue pathology.