Abstract: Objective To screen the genes related with leukocyte responses in mice early after burn injury by bioinformatic analysis of the gene expression profiling data. Methods The gene expression profiles were obtained from GEO (GSE7404, Mouse musculus, 25% TBSA, full-thickness) database. T test, fold changes and GO functional enrichment analysis were used to identify the differentially expressed genes (DEGs) related to leukocyte responses to burns; the interacting genes were transferred to STRING to construct the protein-protein interaction (PPI) network. Biological annotation of the sub-networks was executed using the software Cytoscape. Real-time PCR and Western blotting were used to verify the DEGs in mice. Results In mice at 1 day post-burn, a total of 658 genes were up-regulated and 1167 were down-regulated. PPI network and module analysis suggested that some of the genes (Stat1, Cdk1, Cd19, Lck and Jun) may play critical roles in the PPI network post-burn. Real-time PCR and Western blotting results in mice were consistent with those of bioinformatic analysis of Stat1, Cdk1 and Jun. Conclusion Stat1, Cdk1 and Jun might be critical players in the development of leukocyte response in mice early after burn injury. Our finding provides new insights into the pathogenesis of leukocyte response to burn injury and identifies several biomarkers as potential targets for burn treatment.