Single-Cell RNA-Sequencing Reveals the Transcriptional Landscape of Lacrimal Gland in GVHD Mouse Model

Purpose To investigate the global transcriptional landscape of lacrimal gland cell populations in the GVHD mouse model. Methods Single-cell RNA sequencing and further bioinformatic analysis of dissociated lacrimal gland (LG) cells from the mouse model were performed. Parts of transcriptional results were confirmed by immunofluorescence staining. Results We identified 23 cell populations belonging to 11 cell types. In GVHD LG, the proportion of acinar cells, myoepithelial cells, and endothelial cells was remarkably decreased, while T cells and macrophages were significantly expanded. Gene expression analysis indicated decreased secretion function, extracellular matrix (ECM) synthesis, and increased chemokines of myoepithelial cells. A newly described epithelial population named Lrg1high epithelial cells, expressing distinct gene signatures, was exclusively identified in GVHD LG. The fibroblasts exhibited an inflammation gene pattern. The gene pattern of endothelial cells suggested an increased ability to recruit immune cells and damaged cell-cell junctions. T cells were mainly comprised of Th2 cells and effective memory CD8+ T cells. GVHD macrophages exhibited a Th2 cell-linked pattern. Conclusions This single-cell atlas uncovered alterations of proportion and gene expression patterns of cell populations and constructed cell-cell communication networks of GVHD LG. These data may provide some new insight into understanding the development of ocular GVHD.