Cellular heterogeneity and immune microenvironment revealed by single-cell transcriptome in venous malformation and cavernous venous malformation

Venous malformation (VM) and cavernous venous malformation (CVM) are two types of vascular malformations. Even if the two diseases are similar in appearance and imaging, the distinct cellular components and signaling pathways between them might help distinguish the two from a molecular perspective. Here, we performed single-cell profiling of 35,245 cells from two VM samples and three CVM samples, with a focus on endothelial cells (ECs), smooth muscle cells (SMCs) and immune microenvironment (IME). Clustering analysis based on differential gene expression unveiled 11 specific cell types, and determined CVM had more SMCs. Re-clustering of ECs and SMCs indicated CVM was dominated by arterial components, while VM is dominated by venous components. Gene set variation analysis suggested the activation of inflammation-related pathways in VM ECs, and upregulation of myogenesis pathway in CVM SMCs. In IME analysis, immune cells were identified to accounted for nearly 30% of the total cell number, including macrophages, monocytes, NK cells, T cells and B cells. Notably, more macrophages and monocytes were discovered in VM, indicating innate immune responses might be more closely related to VM pathogenesis. In addition, angiogenesis pathway was highlighted among the significant pathways of macrophages & monocytes between CVM and VM. In VM, VEGFA was highly expressed in macrophages & monocytes, while its receptors were all abundantly present in ECs. The close interaction of VEGFA on macrophages with its receptors on ECs was also predicted by CellPhoneDB analysis. Our results document cellular composition, significant pathways, and critical IME in CVM and VM development.