Abstract

Background Systemic sclerosis (SSc) is a rare autoimmune disease characterized by extensive skin fibrosis. There are no effective treatments due to the severity, multiorgan presentation, and variable outcomes of the disease. Here, integrated bioinformatics was employed to discover tissue-specific expressed hub genes associated with SSc, determine potential competing endogenous RNAs (ceRNA) regulatory networks, and identify potential targeted drugs. Methods In this study, four datasets of SSc were acquired. To identify the genes specific to tissues or organs, the BioGPS web database was used. For differentially expressed genes (DEGs), functional and enrichment analyses were carried out, and hub genes were screened and shown in a network of protein-protein interactions (PPI). The potential lncRNA–miRNA–mRNA ceRNA network was constructed using the online databases. The specifically expressed hub genes and ceRNA network were validated in the SSc mouse and in normal mice. We also used the receiver operating characteristic (ROC) curve to determine the diagnostic values of effective biomarkers in SSc. Finally, the Drug-Gene Interaction Database (DGIdb) identified specific medicines linked to hub genes. Results The pooled datasets identified a total of 254 DEGs. The tissue/organ-specifically expressed genes involved in this analysis are commonly found in the hematologic/immune system and bone/muscle tissue. The enrichment analysis of DEGs revealed the significant terms such as regulation of actin cytoskeleton, immune-related processes, the VEGF signaling pathway, and metabolism. Cytoscape identified six gene cluster modules and 23 hub genes. And 4 hub genes were identified, including Serpine1, CCL2, IL6, and ISG15. Consistently, the expression of Serpine1, CCL2, IL6, and ISG15 was significantly higher in the SSc mouse model than in normal mice.