Lung cancer in never smokers (LCNS) has been considered as a separate disease and the 7thcause of cancer-related death worldwide. However, limited research has focused on “female” cohorts, which have presented a higher incidence rate. In this study, the microarray data of lung cancer tissues derived from 54 female lung cancer patients, consisting of 43 nonsmokers and 11 smokers, were selected from GSE2109 dataset. A total of 249 differentially expressed genes (DEGs) including 102 up- and 147 down-regulated genes were identified and further analyzed for gene ontology (GO) terms and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment. By constructing protein–protein interaction (PPI) network and calculating key modules, 10 hub genes were screened out. The module analysis of the PPI network presented that the progression of female LCNS was significantly associated with immune response as chemokine activity and lipopolysaccharide response, and these biological processes (BP) might be mediated by chemokine signaling pathway and cytokine-cytokine receptor interaction. Then, survival analysis by Kaplan–Meier (K–M) Plotter online platform presented down-regulated gene colony stimulating factor 2 receptor beta common subunit (CSF2RB) of female LCNS might be involved in poor clinical outcome. Female LCNS with high expression of CSF2RB might be relevant with relative risk reduction of mortality, longer median survival time and higher 5-year survival rate, while female LCNS with low expression of CSF2RB might be implicated in a poor clinical outcome. In short, our results support CSF2RB to be a candidate survival predictor for female LCNS.