Identifying specific TLS-associated genes as potential biomarkers for predicting prognosis and evaluating the efficacy of immunotherapy in soft tissue sarcoma

Background: The tertiary lymphatic structure (TLS) is an important component of the tumor immune microenvironment and has important significance in patient prognosis and response to immune therapy. However, the underlying mechanism of TLS in soft tissue sarcoma remains unclear.

Methods: A total of 256 RNAseq and 7 single-cell sequencing samples were collected from TCGA-SARC and GSE212527 cohorts. Based on published TLS-related gene sets, four TLS scores were established by GSVA algorithm. The immune cell infiltration was calculated via TIMER2.0 and “MCPcounter” algorithms. In addition, the univariate, LASSO, and multivariate-Cox analyses were used to select TLS-related and prognosis-significant hub genes. Single-cell sequencing dataset, clinical immunohistochemical, and cell experiments were utilized to validate the hub genes.

Results: In this study, four TLS-related scores were identified, and the total-gene TLS score more accurately reflected the infiltration level of TLS in STS. We further established two hub genes (DUSP9 and TNFSF14) prognosis markers and risk scores associated with soft tissue sarcoma prognosis and immune therapy response. Flow cytometry analysis showed that the amount of CD3, CD8, CD19, and CD11c positive immune cell infiltration in the tumor tissue dedifferentiated liposarcoma patients was significantly higher than that of liposarcoma patients. Cytological experiments showed that soft tissue sarcoma cell lines overexpressing TNFSF14 could inhibit the proliferation and migration of sarcoma cells.

Conclusion: This study systematically explored the TLS and related genes from the perspectives of bioinformatics, clinical features and cytology experiments. The total-gene TLS score, risk score and TNFSF14 hub gene may be useful biomarkers for predicting the prognosis and immunotherapy efficacy of soft tissue sarcoma.