Identification of biomarkers associated with the invasion of nonfunctional pituitary neuroendocrine tumors based on the immune microenvironment.

Introduction: The invasive behavior of nonfunctioning pituitary neuroendocrine tumors (NF-PitNEts) affects complete resection and indicates a poor prognosis. Cancer immunotherapy has been experimentally used for the treatment of many tumors, including pituitary tumors. The current study aimed to screen the key immune-related genes in NF-PitNEts with invasion. Methods: We used two cohorts to explore novel biomarkers in NF-PitNEts. The immune infiltration-associated differentially expressed genes (DEGs) were obtained based on high/low immune scores, which were calculated through the ESTIMATE algorithm. The abundance of immune cells was predicted using the ImmuCellAI database. WGCNA was used to construct a coexpression network of immune cell-related genes. Random forest analysis was used to select the candidate genes associated with invasion. The expression of key genes was verified in external validation set using quantitative real-time polymerase chain reaction (qRT‒PCR). Results: The immune and invasion related DEGs was obtained based on the first dataset of NF-PitNEts (n=112). The immune cell-associated modules in NF-PitNEts were calculate by WGCNA. Random forest analysis was performed on 81 common genes intersected by immune-related genes, invasion-related genes, and module genes. Then, 20 of these genes with the highest RF score were selected to construct the invasion and immune-associated classification model. We found that this model had high prediction accuracy for tumor invasion, which had the largest area under the receiver operating characteristic curve (AUC) value in the training dataset from the first dataset (n=78), the self-test dataset from the first dataset (n=34), and the independent test dataset (n=73) (AUC=0.732/0.653/0.619). Functional enrichment analysis revealed that 8 out of the 20 genes were enriched in multiple signaling pathways. Subsequently, the 8-gene (BMP6, CIB2, FABP5, HOMER2, MAML3, NIN, PRKG2 and SIDT2) classification model was constructed and showed good efficiency in the first dataset (AUC=0.671). In addition, the expression levels of these 8 genes were verified by qRT‒PCR. Conclusion: We identified eight key genes associated with invasion and immunity in NF-PitNEts that may play a fundamental role in invasive progression and may provide novel potential immunotherapy targets for NF-PitNEts.