BACKGROUND: The study of molecular mechanisms characterizing disease progression may be relevant to get insights into systemic sclerosis (SSc) pathogenesis and to intercept patients at very early stage. We aimed at investigating the proteomic profile of preclinical systemic sclerosis (PreSSc) via a discovery/validation two-step approach. METHODS: SOMAcan aptamer-based analysis was performed on a serum sample of 13 PreSSc (discovery cohort) according to 2001 LeRoy and Medsger criteria (characterized solely by Raynaud phenomenon plus a positive nailfold capillaroscopy and SSc-specific antibodies without any other sign of definite disease) and 8 healthy controls (HCs) age, gender, and ethnicity matched. Prospective data were available up to 4±0.6 years to determine the progression to definite SSc according to the EULAR/ACR 2013 classification criteria. In proteins with relative fluorescence units (RFU) > |1.5|-fold vs HCs values, univariate analysis was conducted via bootstrap aggregating models to determine the predicting accuracy (progression vs non-progression) of categorized baseline protein values. Gene Ontologies (GO terms) and Reactome terms of significant proteins at the adjusted 0.05 threshold were explored. Significant proteins from the discovery cohort were finally validated via ELISAs in an independent validation cohort of 50 PreSSc with clinical prospective data up to 5 years. Time-to-event analysis for interval-censored data was used to evaluate disease progression. RESULTS: In the discovery cohort, 286 out of 1306 proteins analyzed via SomaScan, were differentially expressed versus HCs. Ten proteins were significantly associated with disease progression; analysis through GO and Reactome showed differentially enriched pathways involving angiogenesis, endothelial cell chemotaxis, and endothelial cell chemotaxis to fibroblast growth factor (FGF). In the validation cohort, endostatin (HR=10.23, CI95=2.2-47.59, p=0.003) was strongly associated with disease progression, as well as bFGF (HR=0.84, CI95=0.709-0.996, p=0.045) and PAF-AHβ (HR=0.372, CI95=0.171-0.809, p=0.013) CONCLUSIONS: A distinct protein profile characterized PreSSc from HCs and proteins associated with hypoxia, vasculopathy, and fibrosis regulation are linked with the progression from preclinical to definite SSc. These proteins, in particular endostatin, can be regarded both as markers of severity and molecules with pathogenetic significance as well as therapeutic targets.