SAPHO syndrome is a rare disease characterized by inflammatory lesions on skin and bones. Diversified manifestation and inadequate understanding of etiology has limited its diagnosis and treatment. The co-occurrence of other immune-mediated diseases strongly suggests an involvement of autoimmunity in SAPHO syndrome. However, the role of the largest population of circulating immune cells, neutrophils, is still not well explored. In this study, we performed RNA sequencing to profile the mRNA expression of neutrophils purified from peripheral blood of SAPHO patients to identify key genes associated with SAPHO syndrome, trying to find new functional molecules or biomarkers for this rare disease. A total of 442 differentially expressed genes were identified (p < 0.05, fold change > 2), in which 294 genes were upregulated and 148 genes were downregulated. Five differentially expressed genes of interest were verified by quantitative Real-Time Polymerase Chain Reaction (qRT-PCR), among which S100A12 was upregulated and positively related to high-sensitivity C-reactive protein (hsCRP), while the downregulated gene MYADM was positively related to osteocalcin. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that differentially expressed genes were enriched in "systemic lupus erythematosus" and "ECM-receptor interaction". Gene ontology (GO) enrichment showed that differentially expressed genes may participate in biological processes such as "cell migration" and "cell adhesion". In conclusion, this study provides a first insight into transcriptome characteristics of SAPHO syndrome, indicating an over-active neutrophil recruitment in patients and possibly suggesting molecular candidates for further study on diagnosis and pathology of this disease.