A comprehensive survey of immune cytolytic activity-associated gene co-expression networks across 17 tumor and normal tissue types

Cytolytic immune activity in solid tissue can be quantified by transcript levels of two genes, GZMA and PRF1, which is named the CYT score. A previous study has investigated the molecular and genetic properties of tumors associated CYT, but a systematic exploration of how co-expression networks across different tumors are shaped by anti-tumor immunity is lacking. Here, we examined the connectivity and biological themes of CYT-associated modules in gene co-expression networks of 14 tumor and 3 matched normal tissues constructed from the RNA-Seq data of the “The Cancer Genome Atlas” project. We first found that tumors networks have more diverse CYT-correlated modules than normal networks. We next identified and investigated tissue-specific CYT-associated modules across 14 tumor types. Finally, a common CYT-associated network across 14 tumor types was constructed. Two common modules have mixed signs of correlation with CYT in different tumors. Given the tumors and normal tissues surveyed, our study presents a systematic view of the regulation of cytolytic immune activity across multiple tumor tissues.