Identification of key protein-coding genes and lncRNAs in spontaneous neutrophil apoptosis

Polymorphonuclear leukocytes (PMNs) are the most abundant cells of the innate immune system in humans, and spontaneous PMN apoptosis plays crucial roles in maintaining neutrophil homeostasis and resolving inflammation. However, the detailed mechanisms of spontaneous PMN apoptosis remain to be elucidated. By analysis of the public microarray dataset GSE37416, we identified a total of 3050 mRNAs and 220 long non-coding RNAs (lncRNAs) specifically expressed during PMN apoptosis in a time-dependent manner. By short time-series expression miner (STEM) analysis, Gene Ontology analysis, and lncRNA-mRNA co-expression network analyses, we identified some key molecules specifically related to PMN apoptosis. STEM analysis identified 12 gene profiles with statistically significance, including 2 associated with apoptosis. Protein-protein interaction (PPI) network analysis of the genes from 2 profiles and lncRNA-mRNA co-expression network analysis identified a 12-gene hub (including NF kappa B1 and BIRC3) associated with apoptosis, as well as 2 highly correlated lncRNAs (THAP9-AS1, and AL021707.6). We experimentally examined the expression profiles of two mRNA (NF kappa B1 and BIRC3) and two lncRNAs (THAP9-AS1 and AL021707.6) by quantitative real-time polymerase chain reaction to confirm their time-dependent expressions. These data altogether demonstrated that these genes are involved in the regulation of spontaneous neutrophil apoptosis and the corresponding gene products could also serve as potential key regulatory molecules for PMN apoptosis and/or therapeutic targets for over-reactive inflammatory response caused by the abnormality in PMN apoptosis.