Identification of Key Genes and FUNCTIONAL Pathway in Radioresistance of Non-Small Cell Lung Cancer

PURPOSE: For better understanding of radiotherapy resistance and its potential mechanism. METHODS: We established radioresistance cell lines of non-small cell lung cancer (NSCLC) followed by microarray analysis. 529 differentially expressed genes (DEGs) were then screened between radiation resistant cell lines compared with the sensitive cell lines. The biological functions and enrichment pathways of the above DEGs were identified using Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analyses. Gene Set Enrichment Analysis (GSEA) revealed that the radiation resistance group had the most gene sets enriched in altered immune response, such as TNF signaling pathway, when compared to the radiation sensitive group. Protein-protein interaction (PPI) network was carried out through the STRING database, and then five hub genes (CXCL10, IFIH1, DDX58, CXCL11, RSAD2) were screened by Cytoscape software. RT-PCR confirmed the expression of the above hub genes. ChIP-X Enrichment Analysis showed that STAT1 might be the transcription factor of the above hub genes. Considering that PD-L1 could be activated by STAT1 in a variety of tumors and ultimately lead to immune exhaustion, RT-PCR and Western blot verified the expression level of PD-L1. RESULTS: Five hub genes (CXCL10, IFIH1, DDX58, CXCL11, RSAD2) were screened and verified to be highly expressed in radioresistance group, STAT1 might be the transcription factor of the above hub genes. Our study found that the expression level of PD-L1 was increased after radiotherapy resistance. CONCLUSION: Although immune system activation occurs followed by radiation resistance, we hypothesized that the upregulation of PD-L1 expression caused by STAT1 activation might be one of the mechanisms of radiotherapy resistance.