Integrated analysis of dermal blister fluid proteomics and genome-wide skin gene expression in systemic sclerosis: an observational study

Background Skin fibrosis is a hallmark feature of systemic sclerosis. Skin biopsy transcriptomics and blister fluid proteomics give insight into the local environment of the skin. We have integrated these modalities with the aim of developing a surrogate for the modified Rodnan skin score (mRSS), using candidate genes and proteins from the skin and blister fluid as anchors to identify key analytes in the plasma. Methods In this single-centre, prospective observational study at the Royal Free Campus, University College London, London, UK, transcriptional and proteomic analyses of blood and skin were performed in a cohort of patients with systemic sclerosis (n=52) and healthy controls (n=16). Weighted gene co-expression network analysis was used to explore the association of skin transcriptomics data, clinical traits, and blister fluid proteomic results. Candidate hub analytes were identified as those present in both blister and skin gene sets (modules), and which correlated with plasma (module membership > 0.7 and gene significance > 0.6). Hub analytes were confirmed using RNA transcript data obtained from skin biopsy samples from patients with early diffuse cutaneous systemic sclerosis at 12 months. Findings We identified three modules in the skin, and two in blister fluid, which correlated with a diagnosis of early diffuse cutaneous systemic sclerosis. From these modules, 11 key hub analytes were identified, present in both skin and blister fluid modules, whose transcript and protein levels correlated with plasma protein concentrations, mRSS, and showed statistically significant correlation on repeat transcriptomic samples taken at 12 months. Multivariate analysis identified four plasma analytes as correlates of mRSS (COL4A1, COMP, SPON1, and TNC), which can be used to differentiate disease subtype. Interpretation This unbiased approach has identified potential biological candidates that might be drivers of local skin pathogenesis in systemic sclerosis. By focusing on measurable analytes in the plasma, we generated a promising composite plasma protein biomarker that could be used for assessment of skin severity, case stratification, and as a potential outcome measure for clinical trials and practice. Once fully validated, the biomarker score could replace a clinical score such as the mRSS, which carries substantial variability.