Identification and validation of TNFRSF4 as a high-profile biomarker for prognosis and immunomodulation in endometrial carcinoma
BMC Cancer volume 22, Article number: 543 (2022)
The interaction between tumor microenvironment (TME) and tumors offers various targets in mounting anti-tumor immunotherapies. However, the prognostic biomarkers in endometrial carcinoma (EC) are still limited. Here, we aimed to analyze the TME features and identify novel prognostic biomarkers for EC.
ESTIMATE, CIBERSORT, protein-protein interaction (PPI) network, univariate and multivariate Cox regression, and functional enrichment analysis were performed to identify immune- and survival-related hub genes as well as possible molecular mechanisms. The limma package and deconvolution algorithm were adopted to estimate the abundance of tumor-infiltrating immune cells (TICs) and their relationship with the target gene. In the validation section, tissue microarrays (TMAs) of EC and multiplex immunohistochemistry (m-IHC) were evaluated to validate the expression of TNFRSF4, and its correlation with immune markers, including CD4, CD8, and FOXP3. Besides, the receiver operating characteristic (ROC) curve was plotted to determine the diagnostic performance of TNFRSF4, CD4, CD8, and FOXP3 in EC.
Two genes, TNFRSF4 and S1PR4, were screened out from 386 intersection differential expression genes (DEGs) shared by ImmuneScore and StromalScore in EC. Highlighted by TNFRSF4, we found that it was not only positively correlated with the TICs (mainly CD4+ T cells, CD8+ T cells, and Tregs) but significantly related to the prognosis in patients of EC, both verified by data from The Cancer Genome Altas (TCGA)-EC database and clinical samples. At the same time, the expression trend of TNFRSF4 was further confirmed by an integrated meta-analysis based on six microarrays from the Gene Expression Omnibus database (GEO).
Collectively, TNFRSF4, a previously unrecognized key player in EC, could serve as a potential biomarker for prognosis prediction and immunomodulation of EC.
Endometrial carcinoma (EC), the second most common carcinoma of the female genital tract globally , originates from the inner epithelial lining of the uterus. The most common clinical symptoms of EC include abnormal vaginal bleeding or discharge, abdominal pain, especially for postmenopausal women. According to estimates in 2020, there will be 65,620 new cases of EC diagnosed, 12,590 deaths caused by EC in the United States , and both morbidity and mortality rates have exhibited a sustained rise, particularly in cases among women under the age of 50 years (1.3% per year from 2007 ~ 2016) . Although the prognosis of EC is well generally, however, most of the patients with grade III endometrioid endometrial carcinoma (EEC) and non-EEC (include serous clear cell, undifferentiated carcinomas and carcinosarcomas), who are confronted with enormous risks for metastasis or systemic recurrence . Once it happens, the median survival for EC patients is about 4 months. In addition, novel therapies such as optimized individual treatment based on molecular classification are not well defined in terms of surgical staging, surveil-lance scheduling and adjuvant therapy [5, 6].
The tumor microenvironment (TME) is necessary for the living of tumor cells. It, either directly or indirectly, affects tumor occurrence and development through maintaining a quiescent state of immune contexture , promoting tumor angiogenesis , changing biological features of carcinoma , or even regulating the cancer stem cell activity [10, 11]. The orchestration between tumor and TME not only engenders tumor expansion, but also increased tumor mutation probability by creates heterologous TME. Additionally, TME is inextricably linked to tumor immune suppression or activation. In addition to tumor cells, the stromal cells, chemokines, and cytokines, TME comprises innate immune cells (including macrophages, neutrophils, dendritic cells, myeloid-derived suppressor cells, and natural killer cells) and adaptive immune cells (including T cells and B cells) . Both populations of immune cells innately modulate tumor cell-intrinsic and extrinsic processes within TME. To date, some technique has been devoted to characterize the crosstalk between the tumor and its microenvironments, such as multiplex immunohistochemistry (m-IHC) and bulk tissue microarrays (TMAs) [13,14,15]. ESTIMATE (estimation of stromal and immune cells in malignant tumor tissues using expression data) is a highly rated tool developed by Kosuke Yoshihara and his co-workers for predicting tumor purity mainly based on the ssGSEA algorithm , which has helped us to estimate the abundance of tumor-infiltrating lymphocytes (TILs) and subsequently stratify patients for predicting clinical outcomes . Increasing evidence has elucidated that the difference in the efficacy of tumor immunotherapy was mainly due to the heterogeneity of TME . Even though most of these valuable methods and theories have been well applied and tested in many cancers, many works are still needed to make sense of the correlation between TME biological characteristics and the aberrant expression of immune-related genes in EC. Therefore, more reliable and effective surrogate biomarkers are left to be explored to predict the pathological behavior of EC and enable the stratification of patients according to immune-related criteria for improving prognostic and selecting appropriate adjuvant therapy to guide clinical decisions.
Here, we used the ESTIMATE algorithms and the tumor-infiltrating immune cells (TICs) profile to perform a comprehensive analysis of TME and detect related gene expression in patients with EC, through which TNFRSF4 was identified to be associated with the prognosis of EC as a crucial indicator of TME remodeling. Moreover, to support our findings, clinical specimens were applied to validate the expression of TNFRSF4 in EC and adjacent normal tissues. The correlations of TNFRSF4 with clinicopathologic features and immune-related markers (including CD4, CD8, and FOXP3) were evaluated. Overall, our results indicated that TNFRSF4 might be a biomarker of prognostic values, as a crucial role in TME of EC.
The research thoughts and structures of this study were presented in Supplementary Fig. 1.
Raw data acquisition and processing
Transcriptome RNA-seq data of 587 EC cases (552 tumor samples and 35 controls) and corresponding clinical information were obtained from the TCGA database (https://portal.gdc.cancer.gov/).
Estimation evaluation and DEGs generation
Estimate algorithm was applied to determine the immune-stromal component in TME of each sample utilizing estimate package in R software (https://r-forge.r-project.org/), the final results of which were calculated as StromalScore, ImmuneScore, and ESTIMATEScore, representing the ratio of stromal cell, immune cell, and the summation of both cells, respectively . Based on the median of the StromalScore and ImmuneScore, samples were divided into either the high or low score groups, and survival analysis was carried out by R software loaded with survival and survminer packages. Five hundred forty one tumor samples were finally analyzed with complete survival information, and the survival curve was produced by the Kaplan-Meier method and compared between groups by the log-rank test. At the same time, clinical data were acquired and analyzed to ascertain the correlation between scores and clinical parameters. Then, differentially expressed genes (DEGs) of StromalScore and ImmuneScore between the high and low score groups were screened out by limma package  with false discovery rate (FDR) < 0.05 and |log fold change (FC)| > 1.
Functional enrichment analysis
DEGs were displayed and processed as heatmaps according to the expression level using the pheatmap package for clustering (https://CRAN.R-project.org/package=pheatmap). To determine the functional enrichment for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis, R packages (clusterProfiler, enrichplot, org.Hs.eg.db, and ggplot2) were utilized to identify gene functions and achieve visualization for functional profiles, where q-value of less than 0.05 were considered significantly enriched.
PPI network construction and cox regression analysis
Protein-protein interaction (PPI) network was produced through the online analytical database STRING (https://string-db.org/cgi/input.pl) to clarify the interactive relation among genes (only were the edges and nodes kept when the confidence of interactive relationships > 0.4). To further confirm the highly interconnected genes, the MCODE plugin of Cytoscape (https://cytoscape.org/) was applied, and the interactive hub cluster was selected with 56 hub genes and 794 edges (degree cutoff = 0.2, K-core = 2). Additionally, Cox regression analysis was conducted by using R package survival. Genes with P-value < 0.05 were displayed.
Gene set enrichment analysis
H.all.v7.2.symbols Hallmark and C7.all.v7.2.symbols immunological signatures gene sets were downloaded from the Molecular Signatures Database (http://software.broadinstitute.org/gsea/msigdb/index.jsp). All gene sets were analyzed utilizing the software of GSEAS-4.0 downloaded from Broad Institute.
The limma package in R was applied to normalize the data to evaluate the proportion of TICs, and then a standardized gene expression profile was uploaded to CIBERSORT. The deconvolution algorithm was adopted to estimate the TIC abundance . Only 235 tumor samples with P-value < 0.05 were screened out by quality filtering and applied to the following analysis. Meanwhile, single sample geneset enrichment analysis (ssGSEA) was also introduced here to investigate the immune infiltration landscape of T cell subsets of EC to determine their association with the expression of TNFRSF4.
Receiver operating characteristic (ROC) curves
To determine the diagnostic values of interesting genes in EC, pROC package  was used to plot ROC curves, and the values of areas under the curves (AUCs) and corresponding 95% CI were calculated in R.
A meta-analysis based on GEO microarrays to determine the expression of TNFRSF4
The GEO database was searched to retrieve the microarray data on the expression of TNFRSF4 in EC with the retrieval strategies as follows ((‘endometrial’ OR ‘endometrium’) AND (‘cancer’ OR ‘cancerous’ OR ‘carcinoma’ OR ‘carcinomatous’ OR ‘neoplasm’ OR ‘malignancy’ OR ‘tumor’ OR ‘tumors’) AND (‘mRNA’ OR ‘RNA’ OR ‘coding RNA’). Additionally, the inclusion criteria were set to 1) patients diagnosed as EC were involved, 2) the experiments were composed of tumor samples and normal controls, and 3) the expression profiling data of TNFRSF4 was available. The expression profiling data of TNFRSF4 obtained from the queried microarrays were integrated for comparisons between the two groups. Subsequently, a comprehensive meta-analysis by the meta-package was adopted to evaluate the data from all selected microarrays, for which the forest plot presented the expression of TNFRSF4. Furthermore, the DerSimonian-Laird estimator for τ2 and I2 parameters were utilized to determine whether to adopt the corresponding model (random-effects or fixed-effects models) for the pooled estimate. The sensitivity analysis was performed to evaluate the heterogeneity among all investigated studies. Begg’s test and funnel plot were carried out to assess the potential publication bias.
Tissue microarray and ethical statement
TMAs were purchased from Shanghai Outdo Biotech Co. Ltd. (panel HUteA020CS01, HUteA045PG01, and HUteA060CS01). These TMAs, composed of 121 cores (85 carcinomas and 36 adjacent tissue), the detailed information was displayed in Supplementary Table 1. The research scheme of this study has been reviewed and approved by the Ethics Committee of Peking Union Medical College Hospital (ethics, S-K973).
Immunohistochemistry (IHC) and evaluation of immunostaining
IHC staining was performed as described previously [23, 24]. Briefly, TMAs specimens were deparaffinized, hydrated, and subjected to heat-mediated antigen retrieval with EDTA buffer at 95 °C for 10 min. After endogenous peroxidase was quenched and tissue blocked, TMAs sections were incubated with primary antibodies (listed in Supplementary Table 2) and coated at 4 °C overnight, then with biotin-conjugated secondary reagents for 30 min. At the end of the staining, whole TMAs slides were digitally scanned at × 400 using a NanoZoomer S360 (Hamamatsu, Japan) for visualization. Two experienced pathologists independently determined the semi-quantitative immunoreactive scores to evaluate the immunostaining. Given that TNFRSF4, CD8, CD4, and FOXP3 were expressed in TILs, its scoring has used the method of Erminia Massarelli  by counting positive cells in the five random square areas at 400× magnification, and the expression of each marker was recorded as the density of positive cells/mm2.
Multiplex immunohistochemistry (m-IHC) and immunofluorescence
M-IHC staining was performed on TMAs serial sections according to the manufacturer’s protocol (Opal Multiplex IHC Assay Kit; Akoya Biosciences, MA). Briefly, TMAs specimens were baked at 65 °C for 1 h, dewaxed with xylene, rehydrated with a graded series of ethanol solutions ethanol. For epitope retrieval, slides were placed in a microwave with AR6 antigen retrieval buffer for 45 s at 100% power and an additional 15 min at 20% power. After blocking, slides were incubated with the first primary antibody for 1 h at room temperature, followed by Polymer HRP Ms. + Rb and opal fluorophore working solution for 10 min at room temperature, respectively. Slides were stripped with microwave treatment before repeating until all targets of interest were detected. Apply DAPI working solution for 5 min at room temperature and covered with mounting medium, then representative images (Supplementary Fig. 2) from each sample were acquired and analyzed by the Vectra Polaris multispectral slide scanner and supporting software (Akoya Biosciences).
SPSS (Version. 22, Chicago, IL, USA) was used for statistical analyses. Shapiro-Wilk test, combined with normality plots, was applied to determine the normal distribution. For unpaired samples, comparisons between two groups were analyzed by the Mann-Whitney test or Student’s t-test. As for paired samples, the Wilcoxon rank-sum test was applied. Besides, multiple comparisons were carried out based on the Kruskal Wallis rank-sum test to compare groups’ differences. The two-sided Person chi-squared test or Spearman rank correlation test was used to evaluate the relationship between TNFRSF4 IHC scores and immune-related gene expression levels. Disease-specific survival (DSS) was used to analyze the prognostic impact of TNRSF4, which is defined as the time between the date of surgery and the date of death caused by EC. The Kaplan-Meier analysis by the log-rank test. To estimate the prognostic impact of individual continuous gene expression variables and independent prognostic factors, univariate and multivariate Cox regression analyses were performed, respectively. Statistical significance was defined when P-value was less than 0.05.
Clinicopathologic significance of the estimate scores
To reveal the relationship between the component of immune or stromal cells and the clinical outcomes, survival analysis was performed for ImmuneScore, StromalScore, and ESTIMATEScore, respectively. As presented in Fig. 1A, C, the percentage of immune components and ESTIMATEScore were positively correlated with the disease-specific survival rate despite that StromalScore made no difference (Fig. 1B). These results manifested that the immune composition in TME was a more reliable indicator of the prognosis for EC patients. We further estimate the correlation between the three scores and clinicopathologic parameters. The results demonstrated that tumors of EC patients over 60 years old had lower StromalScore than younger patients (P < 0.05) even though ImmuneScore and ESTIMATEScore did not differ in terms of age as shown in Supplementary Fig. 3A, B, C. As for the pathological grade, stage and histology, it was worth noting that our data didn’t exhibited consistently decreasing trends in the three scores as tumor progressed, and not all comparisons among groups with different clinicopathological variables were statistically significant (Supplementary Fig. 3D, E, F, G, H, I, J, K, M). These results implied that advanced age (> 60) and higher grade of EC, mainly grade III, tended to have less immune and stromal component, possibly signified poor prognosis.
DEGs shared by immunescore and stromalscore and functional enrichment
In order to explore the variation of gene profile in the wake of the alteration of the immune or stromal component in TME, we compared the expression level of genes between the high- and low-score groups. Among them, 716 DEGs (552 up-regulated and 164 down-regulated) were determined regarding ImmuneScore, 731 DEGs were screened out from StromalScore, which mainly consisted of highly expressed genes. As shown in the heatmaps, the top 50 DEGs were displayed based on ImmuneScore and StromalScore, and the gene expression pattern showed an apparent difference between groups (Fig. 1D, E). After the intersection, 386 genes were shared in common by ImmuneScore and StromalScore, with 366 up-regulated and the rest down-regulated (Fig. 1F, G). Then, functional enrichment analysis was carried out based on these overlapped genes. Our results indicated that the identified genes were chiefly enriched in the regulation of lymphocyte activation, differentiation, and proliferation in the biological process (BP). About the cellular component (CC) of GO analysis, these genes encoding proteins were the main components of the external side of the plasma membrane, immunological synapse, granule membrane, and protein complex involved in cell adhesion. As for molecular function (MF), the DEGs prevailingly enrolled in immune or cytokine receptor activity (Fig. 1H). KEGG analysis revealed that the chemokine signaling pathway and cytokine-cytokine receptor interaction were most relevant to these DEGs (Fig. 1I). It thus appeared that the overall function of these DEGs primarily focused on immune-related activities, which essentially suggested that the involvement of immune modulation was a notable feature of TME in EC.
Intersection analysis of PPI network and cox regression
To elucidate the underlying mechanism, we explored the PPI network constructed by the STRING database, and thus, the interaction between genes was identified. On this basis, the MCODE plugin of Cytoscape software was applied to seek the hub gene cluster. Finally, 56 genes were identified with 794 edges, displayed in Fig. 2A, and the bar plot in Fig. 2B showed the top 40 genes, ranked by the number of adjacent nodes. Additionally, to evaluate the prognostic impact of 386 DEGs in EC, univariate Cox regression analysis was performed. In order to present the survival-related genes visually, we are considered introducing the so-called forest plot and the results manifested that the expression level of 16 genes was highly correlated with the survival outcome (Fig. 2C). Subsequently, two genes, TNFRSF4 and S1PR4, were screened out after intersection analysis based on the above PPI network and univariate Cox regression analysis (Fig. 2D).
Prognostic significance and indicative role in TME remodeling of TNFRSF4 in EC
Since S1PR4 is not differentially expressed between paired/unpaired ECs and adjacent normal tissues (Supplementary Fig. 4A, B), TNFRSF4 was selected for follow-up studies. TNFRSF4, barely reported on EC but was a high-profile target on other cancers. TNFRSF4 was first discovered on the surface of activated CD4+ T cells, which played a vital role in immune regulation in multiple cancers as a crucial immune checkpoint. In the following analysis, our data indicated that TNFRSF4 showed higher expression levels in tumor tissues than normal tissues either in unpaired or paired samples (P < 0.001 and P = 0.028, respectively), as shown in Fig. 3A, B. To confirm the expression of this gene, six microarrays were eventually involved in the current study from the GEO database, as displayed in Supplementary Table 3. In total, 255 tumor samples and 126 normal tissues were included. Regarding the expression trends of TNFRSF4, a meta-analysis was carried out to quantify its expression according to these datasets accurately. The calculation results were presented in Fig. 4A, for which the fixed-effects model was applied given the existed homogeneity (I2 = 19%, τ2 = 0.0277, P = 0.29). Results also indicated that the expression of TNFRSF4 was dramatically augmented in the EC group (SMD = 0.34, 95% CI [0.11; 0.57], Z = 2.87, P = 0.0041). The publication bias was also assessed, as displayed in Fig. 4B. The funnel plot indicated no bias of the publications (z = 0.94, P = 0.348). Later, a sensitivity test was generated, and no items were determined to exert possible effects on the results (Fig. 4C). Taken together, it was reasonable to determine the up-regulated expression trend of TNFRSF4 in EC tissues.
Interestingly, the expression level decreased as aged when compared with the younger group (Fig. 3C). Similarly, tumors of grade II or III exhibited a consistently less amount of TNFRSF4 when compared with tumors of grade I (Fig. 3D). More importantly, high expression of TNFRSF4 was significantly associated with a better prognosis in EC based on DSS analysis (P = 0.019) (Fig. 3E).
Furthermore, on multivariate analysis, the LASSO regression model was used to reduce the potential collinearity problem and 4 variables were identified (i.e. TNFRSF4, grade, stage, histology), the optimal tuning parameter (λ, lambda.min = 0.015) was selected by cross-validation (Supplementary Fig. 5A,B,). The results suggested that there is no multi-collinearity, and the multivariable Cox regression analysis showed that positive expression of TNFRSF4 was a prognostic factor in terms of improved DSS (HR = 0.317, 95% CI 0.114–0.878, P = 0.027) independent of stage, grade, histology (Supplementary Table 4). Additionally, we also analyzed the distribution among the EC molecular subtype (including POLE, microsatellite-instable/hypermutated (MSI), Copy-Number Low and Copy-Number High) based on TCGA categories. Results show that MSI tumors exhibited higher relative frequencies of positive TNFRSF4 expression than did CN-low tumors. The distribution of TNFRSF4 expression among the four EC molecular subgroups from low to high is POLE, CN-low, CN-high and MSI (Supplementary Fig. 4C).
Then, GSEA was implemented for the high and low-expression groups of TNFRSF4. As shown in Fig. 3F, the genes in the TNFRSF4 high-expression group were chiefly enriched in immune-related activities, such as IFN-γ response, IL-2/STAT5 signaling, and IL-6/JAK/STAT3 signaling. In terms of the TNFRSF4 low-expression group, the genes were mainly involved in the G2M checkpoint, mitotic spindle, mTORC1 signaling, and protein secretion (Fig. 3G). For the C7 collection defined by MSigDB, the immunologic gene sets and multiple immune functional gene sets were enrolled in the high and low TNFRSF4 expression groups (Fig. 3H, I). These results suggested that TNFRSF4 might be a potential indicator for the status of TME.
Correlation of TNFRSF4 with the abundance of TICs
Given the above findings, in order to detect the pertinence relation of TNFRSF4 expression with the immune microenvironment, the component of TICs of each sample was further estimated using the CIBERSORT algorithm. 22 kinds of immune cell profiles in EC samples were identified, as shown in the bar plot (Fig. 5A). It was observed that immune cells in EC were mainly composed of T cells and macrophagocytes. Besides, the correlation among the immune cells was also displayed. The results showed that the proportion of CD8+ T cells was negatively correlated with the presence of CD4+ memory resting T cells and M0 macrophagocytes (correlation coefficient = − 0.55 or − 0.63, respectively). Conversely, the content of CD8+ T cells was positively related to the scale of CD4+ memory-activated T cells, as seen in Fig. 5B.
To verify the primary immune cells affected by TNFRSF4, the difference and correlation analysis were carried out, and results demonstrated that a total of 5 kinds of TICs differed between the high and low-TNFRSF4 expression groups, including CD8+ T cells, regulatory T cells, resting dendritic cells, eosinophils, and neutrophils (Fig. 6A). Furthermore, seven kinds of immune cells were correlated with the expression of TNFRSF4, as observed in Fig. 6B. In addition, four types of immune cells, including CD8+ T cells, regulatory T cells, eosinophils, and neutrophils, were identified to be vitally interconnected with TNFRSF4 (Fig. 6C). T-cell-mediated immune modulation in the immune microenvironment partly depends on the subsets of T cells. In this part, we analyzed the potential relationship between TNFRSF4 and the subsets of T cells. Results showed that TNFRSF4 was positively correlated with subsets of T cells, such as Th1 (R = 0.35), Th17 (R = 0.15) and Treg (R = 0.6), but negative correlated with Th2 (R = − 0.12) (Supplementary Fig. 4D, E). These results hinted that the expression of TNFRSF4 had much to do with the immune activity of TME in EC.
Identify the expression of TNFRSF4 and its correlation with immune-related genes
From the TCGA mRNA expression profiles, we found that the expression of TNFRSF4 was positively correlated with CD4 (R = 0.51), CD8A (R = 0.47), and FOXP3 (R = 0.59) (Fig. 7A, B, C). In IHC staining, representative images of TNFRSF4, CD4, CD8, and FOXP3 were exhibited in Fig. 7D, TNFRSF4, CD4, and CD8 proteins were over-expressed in EC tissues (mainly on tumor immune infiltrating cells) compared with paired adjacent normal tissues or even unpaired normal tissues (Fig. 7E, F). However, we did not detect any evident expression of FOXP3 at the protein level. Besides, there was a significant statistical correlation of TNFRSF4 with CD4 and CD8, except for clinicopathologic parameters (such as age and histological grade) (Supplementary Table 5). To assess the diagnostic value of TNFRSF4 in EC, the ROC diagnosis model was performed. Surprisingly, we uncovered that TNFRSF4 had a higher diagnostic significance in either the TCGA dataset (the AUC value = 0.715) (Fig. 7G) or clinical specimens (the AUC value = 0.777) (Fig. 7H), even compared with the diagnostic capability of CD4 or CD8.
Based on the results from IHC, we further carried out m-IHC to comprehensively investigate the relationship among TNFRSF4, CD4, CD8, and FOXP3 in the TME. As shown in Fig. 7I, TNFRSF4, CD4, and CD8 mainly expressed on the tumor stroma and displayed a co-localization pattern observed in both the separate and merged images. Similarly, no reliable fluorescence staining with FOXP3 was observed. Altogether, these data implied the potential role of TNFRSF4 in immune microenvironment remodeling and diagnostic performance for EC.
Given the limited understanding of the details of TME for EC, especially the complicated and volatile characteristics of the immunological microenvironment, in this study, we comprehensively analyzed the alteration of TME and the composition of TICs in EC-based on CIBERSORT algorithms. We then determined TNFRSF4 and SIRP4 as immune- and prognosis-related hub genes obtained by the intersection between Cox proportional hazards regression and PPI analysis. Finally, focused on TNFRSF4, we found that TNFRSF4 might be responsible for stably maintaining the immunomodulatory characteristics of TME, thus leading to a better prognosis, which was further validated in the TCGA dataset and clinical sample.
Exploring TME signature genes is a necessary step in uncovering intricate relationships among clinical features and molecular characteristics. As a classical tool, the ESTIMATE algorithm has been used to screen out potential novel biomarkers based on stromal and immune gene sets from each sample in the tumor expression matrix. As known, the higher score estimated in ImmuneScore or StromalScore indicated a higher proportion of corresponding immune or stromal cells in TME. Here, our data demonstrated that ImmuneScore had a more significant correlation with the overall survival rate than StromalScore and ESTIMATEScore, which suggested that TME composition was more critical in modulating cancer progression and prognosis of EC patients. To screen out novel genes related to TME in EC, 587 EC cases from the TCGA database were divided into the high and low-score groups based on the median of the StromalScore and ImmuneScore. Our results showed that 386 DEGs were shared in common by ImmuneScore and StromalScore; it suggested that these DEGs played a vital role in modulating the immunological environment and implied a series of intricate interactions among non-neoplastic cells within the TME that could impact tumor fate .
To identify potential biomarkers, the top 40 PPI genes cluster and 16 prognostic-related genes were intersected, leading to TNFRSF4 was surfaced. TNFRSF4, also known as OX40 or CD134, is a co-stimulatory molecule mainly expressed on CD4+ and CD8+ T cells. Initially, evidence was brought by Paterson, D.J , who identified TNFRSF4 as a specific marker of T-cell activation and survival when cross-linked with its ligand TNFSF4. Then emerging evidence highlighted that TNFRSF4 was a promising therapeutic target for T cell-mediated anti-tumor immunotherapy [28,29,30]. However, there are still some contradictory sides of TNFRSF4, especially since the prognostic effects of TNFRSF4 among different kinds of tumors are often quite different. These paradoxical results could be explained by heterogeneous TME and homeostatic regulation of Tregs . Encouragingly, the latest evidence reinforced the notion that immunotherapy with TNFRSF4 agonists would trigger a new trend in cancer therapies . For example, TNFRSF4 stimulation plus a CTLA-4 blockade functioned as a potential therapeutic strategy to eradicate disseminated tumors . In short, TNFRSF4 might act as a win-win path to reestablish T-cell anti-tumor activity.
Strangely, TNFRSF4, as a second-generation immune checkpoint molecule, has been poorly studied in EC. That encourages us to explore whether TNFRSF4 can serve as a reliable prognostic biomarker. Our investigation found that transcriptome signatures of immune-related activities were positively correlated with TNFRSF4 and chiefly enriched in inflammatory signaling pathways, such as IFN-γ response, which is essential for CD8+ T lymphocyte-mediated anti-tumor immunity . Immune cell infiltration has held great promise as a new biomarker of prognosis in patients with different types of solid tumors, including endometrial cancer. To further investigate the role of TNFRSF4 in TME, we firstly estimated the component of TICs in EC samples by the CIBERSORT algorithm. Consistent with previous studies, our results demonstrated that T cells and macrophagocytes were dominant in TME. The proportion of CD8+ T cells was positively related to the amount of CD4+ memory-activated T cells rather than CD4+ memory resting T cells and M0 macrophagocytes. As known, CD4+ T cell help is indispensable for sustaining CD8+ T cell function during chronic viral infection and engaged in anti-tumor activities by CD8+ cytotoxic T cells-dependent apoptosis [35, 36]. Based on our results, further confirming this correlation and evaluating the prognostic significance of TNFRSF4 in large-scale cohort studies are of great importance.
TNFRSF4 can orchestrate various subsets of T cells (including CD8+ and CD4+ T cells, Th1, Th2, Th9, Th17, Tfh, and Tregs ) to coordinate the immune microenvironment. For example, the expression of TNFRSF4 promotes CD8+ and CD4+ T cell proliferation and survival, enhance the Th1-mediated immune response, argument generation and maintenance of Th2, favor Th9 differentiation and modulate the survival and function of Th17, promote Tfh development, antagonize Treg-mediated immune-suppression through direct inhibition of FOXP3 expression [30, 38, 39]. Treg cells act as immunoregulators characterized by the expression of FOXP3, CD4, and CD25 . Under some conditions, TNFRSF4 tends to reduce the generation of Tregs by strongly antagonizing TGF-β and IL-10 and inhibiting the antigen mediated-conversion of naive CD4+ T cells developed into FOXP3+ Tregs [41, 42]. Consistent with that, we find that TNFRSF4 was negatively correlated with FOXP3 in protein levels from the patient sample but positively related to Treg cells in mRNA levels from the TCGA database. In a word, the effects of TNFRSF4 on immune modulation depend on the cellular localization and biology function, notably within the immune microenvironment.
Inspired by these results, we further verify the protein expression of TNFRSF4 in clinical specimens by IHC and m-IHC. Our result showed that TNFRSF4 was overexpressed in cancer tissues compared with adjacent normal tissues, and it was mainly located in the TILs within the tumor stroma. It is worth noting that the distinct location of TNFRSF4, as previously described, indicated different prognostic values. For example, the expression of TNFRSF4 on TILs was correlated with favorable prognosis in human cancers, including advanced gastric cancer, non-small cell lung cancer, ovarian carcinoma, malignant melanoma, and colorectal cancer [25, 43,44,45,46]. However, the expression on cancer cells was associated with a poorer outcome, including hepatocellular carcinoma and cutaneous squamous cell carcinoma [47, 48]. Furthermore, TNFRSF4-positive cells were also observed to interact with CD4+ and CD8+ T cells, possibly due to the role of TNFRSF4 in the proliferation of CD4+ or CD8+ T cells and the survival of antigen-specific memory T cells. Additionally, one unexpected result by ROC curve analysis suggested that TNFRSF4 expression could serve as a potential indicator of favorable diagnosis among patients with EC, and this diagnosis performance was better than both CD4 and CD8.
In summary, TNFRSF4 was identified as a novel molecule in EC; its elevated expression tends to signify favorable clinical outcomes and is closely related to the abundance of CD4+ and CD8+ T cells. As a second immune checkpoint, TNFRSF4 may have potential implications for the prognosis and immunomodulation of EC patients. Notably, our study suggested that further in-depth understanding the role of TNFRSF4 might provide new insights into the immunotherapy of EC.
Tumor-infiltrating immune cells
Differentially expressed genes
False discovery rate
Kyoto Encyclopedia of Genes and Genomes
Single-sample geneset enrichment analysis
Hematoxylin and eosin
Bray F, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
Lortet-Tieulent J, et al. International patterns and trends in endometrial Cancer incidence, 1978-2013. J Natl Cancer Inst. 2018;110(4):354–61.
Zong L, et al. Analysis of the immune checkpoint V-domain Ig-containing suppressor of T-cell activation (VISTA) in endometrial cancer. Mod Pathol. 2022;35(2):266–73.
Chaudhry P, Asselin E. Resistance to chemotherapy and hormone therapy in endometrial cancer. Endocr Relat Cancer. 2009;16(2):363–80.
Makker V, et al. Endometrial cancer. Nat Rev Dis Primers. 2021;7(1):88.
Phuengkham H, et al. Nanoengineered immune niches for reprogramming the immunosuppressive tumor microenvironment and enhancing Cancer immunotherapy. Adv Mater. 2019;31(34):e1803322.
De Palma M, Biziato D, Petrova TV. Microenvironmental regulation of tumour angiogenesis. Nat Rev Cancer. 2017;17(8):457–74.
Tsimberidou AM, et al. Review of precision cancer medicine: evolution of the treatment paradigm. Cancer Treat Rev. 2020;86:102019.
Prasetyanti PR, Medema JP. Intra-tumor heterogeneity from a cancer stem cell perspective. Mol Cancer. 2017;16(1):41.
Driessens G, et al. Defining the mode of tumour growth by clonal analysis. Nature. 2012;488(7412):527–30.
Hinshaw DC, Shevde LA. The tumor microenvironment innately modulates Cancer progression. Cancer Res. 2019;79(18):4557–66.
Ali HR, et al. Patterns of immune infiltration in breast Cancer and their clinical implications: a gene-expression-based retrospective study. PLoS Med. 2016;13(12):e1002194.
Newman AM, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Tan WCC, et al. Overview of multiplex immunohistochemistry/immunofluorescence techniques in the era of cancer immunotherapy. Cancer Commun (Lond). 2020;40(4):135–53.
Yoshihara K, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Bindea G, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95.
Beatty GL, Gladney WL. Immune escape mechanisms as a guide for cancer immunotherapy. Clin Cancer Res. 2015;21(4):687–92.
Ciardullo S, et al. Screening for non-alcoholic fatty liver disease in type 2 diabetes using non-invasive scores and association with diabetic complications. BMJ Open Diabetes Res Care. 2020;8(1):e000904.
Ritchie ME, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Kim K, et al. Immune gene signature delineates a subclass of papillary thyroid Cancer with unfavorable clinical outcomes. Cancers (Basel). 2018;10(12).
Robin X, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.
Zong L, et al. Expression of the immune checkpoint VISTA in breast cancer. Cancer Immunol Immunother. 2020;69(8):1437–46.
Zong L, et al. PD-L1, B7-H3 and VISTA are highly expressed in gestational trophoblastic neoplasia. Histopathology. 2019;75(3):421–30.
Massarelli E, et al. High OX-40 expression in the tumor immune infiltrate is a favorable prognostic factor of overall survival in non-small cell lung cancer. J Immunother Cancer. 2019;7(1):351.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Paterson DJ, et al. Antigens of activated rat T lymphocytes including a molecule of 50,000 Mr detected only on CD4 positive T blasts. Mol Immunol. 1987;24(12):1281–90.
Buchan SL, Rogel A, Al-Shamkhani A. The immunobiology of CD27 and OX40 and their potential as targets for cancer immunotherapy. Blood. 2018;131(1):39–48.
Jahan N, Talat H, Curry WT. Agonist OX40 immunotherapy improves survival in glioma-bearing mice and is complementary with vaccination with irradiated GM-CSF-expressing tumor cells. Neuro-Oncology. 2018;20(1):44–54.
Fu Y, et al. Therapeutic strategies for the costimulatory molecule OX40 in T-cell-mediated immunity. Acta Pharm Sin B. 2020;10(3):414–33.
Bulliard Y, et al. OX40 engagement depletes intratumoral Tregs via activating FcγRs, leading to antitumor efficacy. Immunol Cell Biol. 2014;92(6):475–80.
Cebada J, et al. OX40 agonists for cancer treatment: a patent review. Expert Opin Ther Pat. 2021;31(1):81–90.
Marabelle A, et al. Depleting tumor-specific Tregs at a single site eradicates disseminated tumors. J Clin Invest. 2013;123(6):2447–63.
Rearden R, et al. Control of B-cell lymphoma by therapeutic vaccination and acquisition of immune resistance is independent of direct tumour IFN-gamma signalling. Immunol Cell Biol. 2016;94(6):554–62.
Noël G, Langouo Fontsa M, Willard-Gallo K. The impact of tumor cell metabolism on T cell-mediated immune responses and immuno-metabolic biomarkers in cancer. Semin Cancer Biol. 2018;52(Pt 2):66–74.
Zander R, et al. CD4(+) T cell help is required for the formation of a Cytolytic CD8(+) T cell subset that protects against chronic infection and Cancer. Immunity. 2019;51(6):1028–1042.e4.
Raphael I, et al. T cell subsets and their signature cytokines in autoimmune and inflammatory diseases. Cytokine. 2015;74(1):5–17.
Valzasina B, et al. Triggering of OX40 (CD134) on CD4(+)CD25+ T cells blocks their inhibitory activity: a novel regulatory role for OX40 and its comparison with GITR. Blood. 2005;105(7):2845–51.
Kroemer A, et al. OX40 controls functionally different T cell subsets and their resistance to depletion therapy. J Immunol. 2007;179(8):5584–91.
Fehervari Z, Sakaguchi S. CD4+ Tregs and immune control. J Clin Invest. 2004;114(9):1209–17.
Gramaglia I, et al. The OX40 costimulatory receptor determines the development of CD4 memory by regulating primary clonal expansion. J Immunol. 2000;165(6):3043–50.
Nakae S, et al. Mast cells enhance T cell activation: importance of mast cell costimulatory molecules and secreted TNF. J Immunol. 2006;176(4):2238–48.
Ladányi A, et al. T-cell activation marker expression on tumor-infiltrating lymphocytes as prognostic factor in cutaneous malignant melanoma. Clin Cancer Res. 2004;10(2):521–30.
Petty JK, et al. Survival in human colorectal cancer correlates with expression of the T-cell costimulatory molecule OX-40 (CD134). Am J Surg. 2002;183(5):512–8.
Ramser M, et al. High OX40 expression in recurrent ovarian carcinoma is indicative for response to repeated chemotherapy. BMC Cancer. 2018;18(1):425.
Ohmura H, et al. OX40 and LAG3 are associated with better prognosis in advanced gastric cancer patients treated with anti-programmed death-1 antibody. Br J Cancer. 2020;122(10):1507–17.
Xie K, et al. OX40 expression in hepatocellular carcinoma is associated with a distinct immune microenvironment, specific mutation signature, and poor prognosis. Oncoimmunology. 2018;7(4):e1404214.
Speiser DE, Ho PC, Verdeil G. Regulatory circuits of T cell function in cancer. Nat Rev Immunol. 2016;16(10):599–611.
We acknowledge the TCGA database for providing their platforms and contributors for uploading their meaningful datasets. The authors would like to thank Editage (www.editage.cn) for English language polishing.
This work was supported by grants from the National Natural Science Foundation of China (Grant Nos. 81672648 and 81801433), the CAMS Initiative for Innovative Medicine (Grant Nos. 2016-I2M-1-001), the CAMS Innovation Fund for Medical Sciences (CIFMS) (Grant Nos. 2021-I2M-C&T-B-023), the CAMS Innovation Fund for Medical Sciences health (Grant Nos. 2020-PT320–003), and special funds from the CAMS Molecular Pathology Center and the Central Public Welfare Institutions of CAMS (Nos. 2017PT31008).
Ethics approval and consent to participate
This research involving human clinical specimens has been reviewed and approved by the Ethics Committee of Peking Union Medical College Hospital (ethics, S-K973). Informed consent was not required because this study was a retrospective report of cases, which was waived by the ehical committee of the Peking Union Medical College Hospital.
All methods and analysis were followed the relevant guidelines and regulations approved by the Institutional Review Board of our hospital.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The flow diagram of the research design. This flow chart presented a comprehensive bioinformatics analysis and cohort validation to screen out the putative target gene, TNFRSF4, and investigate its clinicopathologic significance in EC.
Representative images of color schemes for m-IHC. Scale bar: 100 μm.
The correlation between Estimate scores with clinicopathologic parameters. A, B, C, D, E, F, G, H, I, J, K M. Distribution of ImmuneScore, StromalScore, and ESTIMATEScore grouped by age, grade, stage and histology, respectively. Comparisons were performed by the Kruskal Wallis rank-sum test.
Characterize the expression of S1PR4, and identification of EC immune molecular subtypes. A, B. The expression of S1PR4 in unpaired and paired ECs and adjacent normal tissue based on TCGA database. C. The distribution of the expression of TNFRSF4 among the four EC molecular subtype. D. Single sample GSEA (ssGSEA) was performed using TCGA-EC gene sets. The enrichment score of ssGSEA was displayed in the heatmap. E. Validating the correlation of TNFRSF4 expression with four kinds of subsets of T cells. The correlation coefficient and P values were calculated by Spearman correlation analysis.
Independent prognostic marker selected by LASSO Cox regression analysis. A. Lasso coefficient profiles of the 5 co-variates. Each curve was on behalf of a co-variate. B. Partial likelihood deviance was calculated by the cross-validation for the best λ to determine the minimum mean cross-validated error.
Clinical baseline characteristics of all subjects of the TMAs samples.
Detailed information of primary antibody for IHC and m-IHC.
Gene Expression Omnibus datasets involved.
Results of multivariate analyses for identifying independent survival factors in patients with endometrial cancer.
Association among TNFRSF4 expression and clinicopathologic parameters, CD4 and CD8 in patients with EC of the validation cohort.
About this article
Cite this article
Ma, H., Feng, Ph., Yu, Sn. et al. Identification and validation of TNFRSF4 as a high-profile biomarker for prognosis and immunomodulation in endometrial carcinoma. BMC Cancer 22, 543 (2022). https://doi.org/10.1186/s12885-022-09654-6