Uncovering the early interactions and spatial distribution of dermal fibroblasts and immune cells in treatment-naïve diffuse cutaneous systemic sclerosis (dcSSc) patients is critical to understanding the earliest events of skin fibrosis. We generated an integrated multiomic dataset of early, treatment-naïve dcSSc skin. Skin biopsies were analyzed by single-nuclei multiome sequencing (snRNA-seq and snATAC-seq) and two different spatial transcriptomic methods to comprehensively determine the molecular changes in these individuals. We identified an immunomodulatory niche within the papillary, hypodermis, and vascular regions that are enriched for activated myeloid cells and fibroblasts characterized by expression of genes such as CXCL12, APOE, and C7. Pathway analyses showed significant enrichment of PI3K-AKT-mTOR signaling pathway expression in these cellular niches, driven by profibrotic growth factor signaling networks. Macrophage subclustering showed SSc-specific macrophage activation of the IL6-JAK-STAT signaling and the enrichment of oxidative phosphorylation pathways. Ligand-receptor analysis revealed that SSc macrophages secrete PDGF and TGF-β to activate the SSc-dominant fibroblast subclusters. Spatial transcriptomic analyses showed monocyte-derived MRC1+ macrophages express PDGF near PDGFRhighTHY1high fibroblasts. Multi-omic data integration and spatial transcriptomic neighborhood analysis revealed the co-localization of fibroblasts, macrophages, and T cells around the vasculature. These data suggest that interactions between activated immune cells and immunomodulatory fibroblasts around vascular niches are an early event in scleroderma pathogenesis.
Helen C. Jarnagin, Rezvan Parvizi, Zhiyun Gong, Rosemary Gedert, Xianying Xing, Lam (Alex) C. Tsoi, Rachael Bogle, Madeline J. Morrisson, Laurent Perreard, Patricia A. Pioli, Fred Kolling IV, Johann E. Gudjonsson, Dinesh Khanna, Michael L. Whitfield