# Identification of the gene signatures related to NK/T cell communication to evaluate the tumor microenvironment and prognostic outcomes of patients with prostate adenocarcinoma

## TL;DR

This study identifies gene signatures related to NK/T cell communication in prostate cancer, offering new insights into tumor behavior and potential treatment targets.

## Contribution

The study introduces a novel risk model based on NK/T cell communication genes for evaluating prostate cancer prognosis.

## Key findings

- NK/T cells showed the most significant changes in tumor groups compared to benign groups.
- The high-risk group exhibited higher immune infiltration and enrichment in immune-related pathways.
- ISYNA1 knockdown reduced prostate cancer cell proliferation, migration, and invasion.

## Abstract

Prostate adenocarcinoma (PRAD) is a leading cause of male mortality, with NK/T cell communication being key areas of the research.

The Seurat package was utilized to normalize and reduce the dimensionality of the single-cell data, and CellMarker 2.0 was employed for cell annotation. CellChat was utilized to construct the ligand-receptor interaction network of cell subsets. Differentially expressed genes (DEGs) were filtered by the limma package. Univariate Cox and the LASSO regression in the glmnet package were used to obtain biomarkers and develop a risk model. The survminer package was used to calculate the optimal threshold for dividing patients into high-risk and low-risk groups, and then Kaplan-Meier (KM) survival analysis was performed. Single-sample GSEA (ssGSEA), TIMER, and ESTIMATE packages were employed for immune infiltration analysis. Pathway analysis was conducted for the low- and high-risk groups using GSEA. Immunotherapy responses were evaluated by adopting TIDE method. Additional cellular validation (quantitative real-time PCR, CCK-8, Transwell, and scratch assay) was implemented to confirm the effects of feature genes on PRAD.

Compared with the benign group, NK/T cells were the cell type with the greatest changes in the tumor group, and their communication intensity was relatively high among all cell types. A RiskScore model was constructed as follows:
0.579*FOXS1 + 0.345*GPC6 + 0.385*ISYNA1 + 0.418*ITGAX + 0.792*MGAT4B + 0.368*PRR7 + 0.458*REXO2
. Analysis of the differences between the two risk groups showed that the level of immune infiltration was higher in the high-risk group, and it was significantly enriched in immune-correlated pathways, while the low-risk group was mainly enriched in metabolism-related pathways. TIDE analysis indicated that the high-risk group had higher immune escape potential. The cellular validation assays have revealed the higher expression of seven biomarkers in PRAD groups. Further, ISYNA1 knockdown inhibited the proliferation, migration, and invasion ability of PRAD cells.

The current research reveals key communication genes in PRAD, offering new possibilities for the exploration of new therapeutic targets.