To identify the differential expressed genes of neurogenic heterotopic ossification (NHO) based on GEO database and analyze them by bioinformatics to find the disease-related genes and channels. Meth? ods: Download the gene expression microarray about NHO from GEO database. The differential genes of NHO group and healthy control group were screened by Limma package of R Software with the thresholds of adj. P-value < 0.05 and | Log 2 (fold-change) | > 2. The differential genes were enriched by GO and KEGG pathway by Matascape online tool, and the network of protein-protein interaction (PPI) of differential genes was constructed. Cytohubba and Mcode algorithms were used to find the key genes. Results: A total of 276 differential genes were screened, including 150 up-regulated genes and 126 down-regulated genes. The biological process (BP) of differentially expression genes was mainly enriched in modulation of chemical synaptic transmission, regulation of trans-synaptic signaling and ossification. Molecular function (MF) was mainly enriched in extracellular matrix structural constituent, heparin binding and glycosaminoglycan binding. The cellular component (CC) was mainly enriched in collagen trimer, collagen-containing extracellular matrix and extracellular matrix. KEEG signal pathway analysis was mainly enriched in complement and coagulation cascade, IL-17 signal pathway, Wnt signal pathway, protein digestion and absorption and so on. Ten hub genes including EGF, GPR18, ADRA2C, CXCL1, C3, CXCL6, ADRB2, PENK, CXCL2, CDH1 and four seed genes including CXCL2, PTH2, EPHB1, HTR2B were identified in PPI network. Conclusion: Through the bioinformatics analysis of the gene microarray, we found the genes and molecular regulation mechanisms that may be related to NHO.