Characterization of the heterogeneity of endothelial cells in bleomycin-induced lung fibrosis using single-cell RNA sequencing

The loss of normal alveolar capillary and deregulated angiogenesis occurs simultaneously in idiopathic pulmonary fibrosis (IPF); however the contributions of specific endothelial subpopulations in the development of pulmonary fibrosis are poorly understood. Herein, we perform single-cell RNA sequencing to characterize the heterogeneity of endothelial cells (ECs) in bleomycin (BLM)-induced lung fibrosis in rats. One subpopulation, characterized by the expression of Nos3 and Cav1, is mostly distributed in non-fibrotic lungs and also highly expresses genes related to the “response to mechanical stimulus” and “lung/heart morphogenesis” processes. Another subpopulation of ECs expanded in BLM-treated lungs, characterized by Cxcl12, is observed to be closely related to the pro-fibrotic process in the transcriptome data, such as “regulation of angiogenesis,” “collagen binding,” and “chemokine activity,” and spatially localized to BLM-induced neovascularization. Using CellPhoneDB software, we generated a complex cell–cell interaction network, which predicts the potential roles of EC subpopulations in recruiting monocytes, inducing the proliferation of fibroblasts and promoting the production and remolding of the extracellular matrix (ECM). Taken together, our data demonstrate the high degree of heterogeneity of ECs in fibrotic lung and it is proposed that the interaction between ECs, macrophages, and stromal cells contributes to pathologic IPF.