“…DEGs were identified using Seurat FindAllMarkers function with the default test setting (non-parameteric Wilcoxon Rank Sum test). Expression of gene modules in the form of published signature gene sets were calculated with AddModuleScore (Seurat) taking the top DEG from telocytes (10 genes), Lo-1 FB (10 genes) and Lo-2 FB (7 genes) reported by McCarthy et al (McCarthy et al, 2020), the top 10 DEG from FB1-5, MC and SMCs reported by Hong et al (Hong et al, 2020), and the top 20 DEG from Pi16 + , Col15a1 + , Fbln1 + and Bmp4 + FB reported by Buechler et al (Buechler et al, 2021). Pearson correlations between datasets were calculated based on average expressions per cluster of overlapping variable genes and plotted with heatmap.2 (version 3.0.3) (Gregory R. Warnes, Ben Bolker, Lodewijk Bonebakker, Robert Gentleman, Wolfgang Huber, Andy Liaw, Thomas Lumley, Martin Maechler, Arni Magnusson, Steffen Moeller and Venables, 2020).…”