Abstract
The association between immune checkpoint inhibitors (ICIs) and immune gene networks in squamous lung cancer (LUSC) and lung adenocarcinoma (LUAD) was studied. Immune gene networks were constructed using RNA-seq data from the gene expression omnibus (GEO) database. Datasets with more than 10 samples of normal control and tumor tissues were selected; of these, GSE87340, GSE120622, and GSE111907 were suitable for analysis. Gene set enrichment for pathway analysis was performed. For immune gene network construction, 998 unique immune genes were selected from 21 pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG). Gene function annotation was performed based on the KEGG, Gene Ontology, and Reactome databases. Tumor tissues showed decreased coagulation, hematopoiesis, and innate immune pathways, whereas complement- and coagulation-related genes were prominent in the tumor immune gene network. The average numbers of neighbors, clustering coefficients, network diameters, path lengths, densities, and heterogeneities were highest for normal tissue, followed by LUAD and LUSC. Decreased coagulation genes, which were prominent in tumor immune networks, imply functional attenuation. LUAD was deviated from normal tissue, based on network parameters. Tumor tissues showed decreased immune function, and the deviation of LUSC from normal tissue might explain LUSC's better therapeutic response to ICI treatment.