Abstract
AIMS: Endothelial cell (EC) dysfunction drives the initiation and pathogenesis of pulmonary arterial hypertension (PAH). We aimed to characterize EC dynamics in PAH at single-cell resolution. METHODS AND RESULTS: We carried out single-cell RNA sequencing (scRNA-seq) of lung ECs isolated from an EC lineage-tracing mouse model in Control and SU5416/hypoxia-induced PAH conditions. EC populations corresponding to distinct lung vessel types, including two discrete capillary populations, were identified in both Control and PAH mice. Differential gene expression analysis revealed global PAH-induced EC changes that were confirmed by bulk RNA-seq. This included upregulation of the major histocompatibility complex class II pathway, supporting a role for ECs in the inflammatory response in PAH. We also identified a PAH response specific to the second capillary EC population including upregulation of genes involved in cell death, cell motility, and angiogenesis. Interestingly, four genes with genetic variants associated with PAH were dysregulated in mouse ECs in PAH. To compare relevance across PAH models and species, we performed a detailed analysis of EC heterogeneity and response to PAH in rats and humans through whole-lung PAH scRNA-seq datasets, revealing that 51% of up-regulated mouse genes were also up-regulated in rat or human PAH. We identified promising new candidates to target endothelial dysfunction including CD74, the knockdown of which regulates EC proliferation and barrier integrity in vitro. Finally, with an in silico cell ordering approach, we identified zonation-dependent changes across the arteriovenous axis in mouse PAH and showed upregulation of the Serine/threonine-protein kinase Sgk1 at the junction between the macro- and microvasculature. CONCLUSION: This study uncovers PAH-induced EC transcriptomic changes at a high resolution, revealing novel targets for potential therapeutic candidate development.