Skip to main content
Download PDF

Immune characteristics of dedifferentiated retroperitoneal liposarcomas and the reliability of regional samples in evaluating their tumor immune microenvironments

World Journal of Surgical Oncology volume 22, Article number: 25 (2024) Cite this article

Abstract

Background

Tumor immunotherapy is a new treatment breakthrough for retroperitoneal liposarcoma (RPLS), which is highly invasive and has few effective treatment options other than tumor resection. However, the heterogeneity of the tumor immune microenvironment (TIME) leads to missed clinical diagnosis and inappropriate treatment. Therefore, it is crucial to evaluate whether the TIME of a certain part of the tumor reliably represents the whole tumor, particularly for very large tumors, such as RPLS.

Methods

We conducted a prospective study to evaluate the TIME in different regions of dedifferentiated RPLS (DDRPLS) by detecting the expressions of markers such as CD4+, CD8+, Foxp3+, CD20+, CD68+, LAMP3+, PD-1+ tumor-infiltrating lymphocytes (TILs), and PD-L1 in tumors and corresponding paratumor tissues via immunohistochemistry and RNA sequencing.

Results

In DDRPLS, very few TILs were observed. Differentially expressed genes were significantly enriched in cell part and cell functions, as well as the metabolic pathway and PI3K-Akt signaling pathway. In addition, for most tumors (70–80%), the TIME was similar in different tumor regions.

Conclusions

For most tumors (70–80%), the TIME in any region of the tumor reliably represents the whole tumor. DDRPLS may regulate cell functions by modulating the metabolic and PI3K-Akt signaling pathways to promote its malignant behavior.

Introduction

Retroperitoneal liposarcoma (RPLS) is the most common type of retroperitoneal tumor and is characterized by rapid growth and high invasiveness [1]. At present, there is no effective treatment of RPLS except for surgical resection of the tumor [2, 3]. Moreover, surgery often leads to trauma due to the need for extensive resection of several surrounding organs and tissues [4]. In addition, R0 resection is difficult because the tumor is extremely large and has an unclear boundary. Thus, RPLS, especially dedifferentiated RPLS (DDRPLS), has a high recurrence rate and a poor prognosis [5]. Therefore, it is essential to identify new treatment methods.

Tumor immunotherapy has demonstrated positive therapeutic effects in several tumors [6,7,8,9,10]. Therefore, tumor immunotherapy may represent a promising therapeutic approach for RPLS. However, PD-L1 expression varies across tumor regions [11,12,13], leading to erroneous classification of some PD-L1–positive tumors as PD-L1–negative tumors. In such situations, patients do not receive PD-1/PD-L1 blockers, resulting in the loss of opportunity for possible treatment. Furthermore, the characteristics of the whole tumor cannot be detected. Therefore, evaluation of the tumor immune microenvironment (TIME) before tumor immunotherapy should be based on clinical tumor specimens that accurately reflect the overall tumor characteristics. As a result, it is determined which tumor region accurately reflects the TIME characteristics of the whole tumor. Thus, in this prospective study, we analyzed the TIME characteristics of various tumor regions.

We found that there were very few tumor-infiltrating lymphocytes (TILs) in DDRPLS. Differentially expressed genes (DEGs) were significantly enriched in cell part and cell functions, as well as the metabolic and PI3K-Akt signaling pathways. For most tumors (70–80%), the TIME was similar across different tumor regions. In other words, for most tumors (70–80%), the TIME in any tumor region could reliably represent the whole tumor. DDRPLS may regulate cell functions by modulating the metabolic and PI3K-Akt signaling pathways to promote its malignant behavior.

Materials and methods

Specimen collection and clinicopathological data of patients with DDRPLS

The surgical specimens were obtained from patients with DDRPLS who were treated at Xiang'an Hospital of Xiamen University and the International Hospital of Peking University. All patients provided their written informed consent, and the study protocol was approved by the Ethics Committee of Xiang’an Hospital of Xiamen University. The study was conducted in line with the Declaration of Helsinki. DDRPLS was diagnosed on the basis of standard clinical and histological criteria. All patients underwent tumor resection. Finally, 50 tumor specimens from 10 newly resected tumors and 10 corresponding paratumor tissues were obtained from 10 patients (Fig. 1). The median tumor diameter was 22.75 cm (range, 10–29.9 cm). The clinicopathological data of patients and their tumors are presented in Table 1. Fifty specimens (T1-5) were used for immunohistochemistry (IHC) analysis and 30 specimens (T1, T3, and T5) with corresponding paratumor specimens were used for RNA sequencing (RNA-seq) (Fig. 1).

Fig. 1
figure 1

Experimental flowchart. First, the newly resected tumors were cut along their long axis. Next, specimens were obtained from T1-5 sites on the long axis (T1 is the midpoint of the long axis of a tumor; T5 is the outer edge of the long axis of the tumor; T3 is the midpoint of the line between T1 and T5; T2 is the midpoint of the line between T1 and T3; T4 is the midpoint of the line between T3 and T5), and were divided into two groups. In one group, the specimens were placed in 10% formalin for subsequent immunohistochemistry, and in another group, the specimens were placed in liquid nitrogen for subsequent RNA-seq. Finally, the IHC and RNA-seq results were quantified, and the similarity of the TIME in different tumor regions was analyzed using the Mahalanobis distance. FFPE, formalin-fixed and paraffin-embedded; IHC, Immunohistochemistry; and RNA-seq, RNA sequencing

Full size image
Table 1 Clinicopathological features of patients and their tumors
Full size table

RPLS is generally large, which makes it challenging to obtain several samples. Therefore, we hypothesized that RPLS is a round or oval tumor growing from the center to the periphery. In other words, we considered the center and periphery of RPLS to be old and new tissues, respectively. In this case, the largest difference in the TIME within the tumor tissue should be between the center and periphery; that is, if the TIME is similar between the tumor center and periphery, we can infer that the TIME of the whole RPLS is homogeneous. Therefore, we obtained samples from the abovementioned five locations.

IHC

The tissues were routinely dehydrated, embedded, and sliced. After baking the slices, the sections were successively dewaxed, hydrated, permeated, sealed, antigen-repaired, incubated with primary and secondary antibodies, color-rendered, counterstained, and sealed.

Proportion of positive cells: We randomly obtained 10 high-power vision fields (× 400) from each slice, calculated the number of positive cells among 100 cells in each high-power vision field, and determined the average value from 10 high-power vision fields (positive index, expressed as a percentage). The cells were scored as follows: < 1%, score 0; 1–10%, score 1; 10–50%, score 2; and > 50%, score 3.

Cell density (positive cells/mm2): Ten high-power vision fields (× 400; area, 0.0788 mm2) were randomly obtained from each slice. We calculated the number of positive cells in each high-power vision field and determined the average value from 10 high-power fields. Supplementary Table S1 presents the information related to antibodies.

RNA-seq

The total RNA of the samples was extracted using a TRIzol kit (Promega, USA) in line with the manufacturer’s instructions. Residual DNA was removed by digestion with DNase I (TaKaRa Bio, Japan). A260/A280 value was measured to estimate the RNA concentration. A small quantity of RNA product was for agarose gel electrophoresis to determine the RNA quality. Next, the qualified RNAs were quantified using Qubit3.0 with the Qubit™ RNA Broad Range Assay kit (Life Technologies, Q10210). A total of 2 μg total RNA was used for stranded RNA sequencing library preparation using KC™ Stranded mRNA Library Prep Kit for Illumina® (Wuhan Seqhealth, China) in accordance with the manufacturer’s instructions. PCR products with 200–500 bps were enriched, quantified, and finally sequenced on a Novaseq 6000 sequencer (Illumina) using a PE150 model. Raw sequencing data were filtered using Trimmomatic (version 0.36). Low-quality reads were discarded, and the reads contaminated with adaptor sequences were trimmed. Reads mapped to the exon regions of each gene were counted using featureCounts (Subread-1.5.1; Bioconductor). The abundance of each immune cell type was estimated as the average expression level of the corresponding genes [14] (Supplementary Table S2).

Statistical analysis

Differences between groups, correlations between factors, and the similarity of the TIME in different tumor regions were analyzed using the Mann–Whitney U test, Spearman’s correlation coefficients, and Mahalanobis distance analysis on SPSS version 23.0 (SPSS Inc., IL, USA) and GraphPad Prism software version 8.0 (GraphPad Software Inc., San Diego, CA, USA). Statistical significance was defined as p < 0.05.

Mahalanobis distance is the covariance distance in data. It is used to effectively determine the similarity between two unknown sample sets and to effectively evaluate the outliers and similarities among individuals in a population [14]. Additionally, compared with the common Euclidean distance, Mahalanobis distance has two advantages. First, its scale remains unchanged, and thus the differences in unit measurement values of different biomarkers will not affect the analysis results. Second, it also evaluates correlations between covariates, which can capture not only the difference in a single variable but also the differences between a group of variables[15]. In summary, Mahalanobis distance analysis is an ideal method for evaluating the similarity between unknown sample sets [14].

Results

Characteristics of TIME in DDRPLS

To analyze the characteristics of the TIME in different tumor regions, we first detected the expression of immune markers including CD4+, CD8+, Foxp3+, CD20+, CD68+, LAMP3+, PD-1+ TILs, and PD-L1 in different tumor regions using IHC. As shown in Fig. 2A, the TIME was significantly heterogeneous between patients. However, the TIME was similar in different tumor regions (Fig. 2A). Furthermore, for most tumors, the distribution range of immune cells in different tumor regions was narrow (Fig. 2B, Supplementary Fig. 1), indicating similar TIME in these regions.

Fig. 2
figure 2

TIME of DDRPLS. A The TIME between different patients was heterogeneous; however, TIME in different tumor regions was similar (IHC, represented by CD4). B Immune cell constitution in DDRPLS using IHC data. For most tumors, the distribution range of immune cells was narrow for different tumor regions. Each dot denotes the immune cell density of a sample. TIME, tumor immune microenvironment; DDRPLS, dedifferentiated retroperitoneal liposarcoma; and IHC, immunohistochemistry

Full size image

In addition, very few TILs were observed in DDRPLS (Fig. 2, Supplementary Figures S1–S4), which is consistent with a previous report [16]. CD4+ T cells were the most common, followed by CD68+ macrophages and PD-1+ cells, while Foxp3+ Tregs were the least common, followed by CD8+ T cells and PD-L1+ cells (Supplementary Figures S2–S4). Interestingly, PD-1+ cells were found in all DDRPLS specimens (Table 2). However, PD-L1+ cells were very rare and not found in many specimens (Table 3). Additionally, PD-L1 expression was closely related to other TILs and was associated with increased numbers of CD8+ T cells, Foxp3+ Tregs, CD20+ B cells, CD68+ macrophages, and LAMP3+ dendritic cells (DCs) (Fig. 3A).

Table 2 PD-1 expression
Full size table
Table 3 PD-L1 expression
Full size table
Fig. 3
figure 3

PD-L1 and TLS in DDRPLS. A PD-L1 expression increases CD8+ T cells, Foxp3+ Tregs, CD20+ B cells, CD68+ macrophages, and LAMP3+ DCs. B Number or relative gene expression of TLS per sample (left: IHC data; right: RNA-seq data). Each dot denotes the number or relative gene expression of a sample. C TLS increases CD4+ T cells, CD8+ T cells, CD68+ macrophages, and PD-L1+ cells in DDRPLS. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. TLS, tertiary lymphatic structure; DDRPLS, dedifferentiated retroperitoneal liposarcoma; DCs, dendritic cells; IHC, immunohistochemistry; and RNA-seq, RNA sequencing

Full size image

The tertiary lymphoid structure (TLS) is the region of lymphocyte aggregation [17, 18], and it can provide an important functional environment for cellular and humoral immunity to maintain and promote the antitumor immune response [17,18,19]. A higher number of TLS is associated with a stronger antitumor immune response [20, 21]. Therefore, to analyze TLS in the TIME of DDRPLS, the number of TLS per tissue section was determined. As shown in Fig. 3B, 56 TLSs were found in 50 regions of 10 tumors. The highest number of TLS found in one region was 9. Interestingly, compared with tumors without TLS, the numbers of CD4+ T cells, CD8+ T cells, CD68+ macrophages, and PD-L1+ cells in the tumors with TLS were significantly increased (Fig. 3C).

Relationships between immune characteristics and clinicopathological features of DDRPLS patients

As shown in Supplementary Figures S5 and S6 and Supplementary Table S3, the immune characteristics of tumors were related to the clinicopathological features of DDRPLS patients. Higher numbers of CD4+ T cells, CD8+ T cells, CD20+ B cells, CD68+ macrophages, LAMP3+ DCs, PD-L1+ cells, and TLSs were detected in patients with a single tumor than in those with multiple tumors (Supplementary Figure S5). The density of CD4+ T cells (r =  − 0.5187, p = 0.0001), CD8+ T cells (r =  − 0.3221, p = 0.0225), CD68+ macrophages (r =  − 0.4097, p = 0.0031), and PD-L1+ cells (r =  − 0.4119, p = 0.003), and the number of TLSs (r =  − 0.4665, p = 0.0006), were negatively correlated with tumor size (Supplementary Figure S6). However, no significant correlation was found between immune characteristics and diagnostic status (primary or recurrent tumor) (Supplementary Table S3). In addition, higher densities of CD8+ T cells (p = 0.0193), Foxp3+ Tregs (p = 0.0029), CD68+ macrophages (p = 0.0028), and PD-L1+ cells (p = 0.0439) were detected in older patients than in younger patients. Furthermore, higher densities of CD4+ T cells (p = 0.0003), CD8+ T cells (p = 0.0188), CD20+ B cells (p = 0.0193), CD68+ macrophages (p = 0.03), and PD-L1+ cells (p = 0.0034), and a greater number of TLSs (p = 0.0044), were detected in female patients than in male patients (Supplementary Table S3).

Similarities in TIME between different regions of DDRPLS

The abovementioned results demonstrated that the TIME was similar in different regions of DDRPLS. To further confirm this finding, IHC data were analyzed to determine the Mahalanobis distance. As shown in Figure 4A, for most tumors, the TIME was similar across different tumor regions. In 7 of 10 tumors (70%; nos. 1, 2, 4, 5, 6, 9, and 10), the distribution range of Mahalanobis distance was narrow for different tumor regions, while in 3 of 10 tumors (30%; nos. 3, 7, and 8), the distribution range of Mahalanobis distance was wide for different tumor regions (Fig. 4A). In addition, in 8 of 10 tumors (80%; nos. 1, 2, 4, 5, 6, 8, 9, and 10), the Mahalanobis distance for all tumor regions was 0–22 (critical value), while it was more than 22 (critical value) in 2 of 10 tumors (20%; Nos. 3 and 7) (Fig. 4A). These data indicated that for most tumors (70–80%), the TIME was similar for different tumor regions.

Fig. 4
figure 4

Similarity of the TIME across different tumor regions using IHC data. A Mahalanobis distance of all immune markers (including eight immune cells and TLS). B Mahalanobis distance of all immune markers except one. Each dot denotes the Mahalanobis distance of a region. Dotted lines denote the critical value of Mahalanobis distance. TIME, tumor immune microenvironment; IHC, immunohistochemistry; and TLS, tertiary lymphatic structure

Full size image

To analyze the factors affecting the similarity of TIME across different tumor regions, Mahalanobis distance was calculated again after removing markers one by one. As shown in Fig. 4B, after removing CD8+ T cells, Foxp3+ Tregs, LAMP3+ DCs, and PD-1+ cells, the Mahalanobis distance for different tumor regions (no. 8) was significantly narrowed, indicating that CD8+ T cells, Foxp3+ Tregs, LAMP3+ DCs, and PD-1+ cells were the main influencing factors leading to the TIME heterogeneity in different tumor regions (no. 8). For tumor no. 3, not only the distribution range of the Mahalanobis distance for different tumor regions was wide, but also the Mahalanobis distance did not decrease to less than 22 (critical value), regardless of removing the markers (Fig. 4B). This indicated that the TIME was highly heterogeneous across the different regions of tumor no. 3. However, for tumor No. 7, the Mahalanobis distance in all regions decreased to less than 22 (critical value) after removing PD-1+ cells (Fig. 4B), indicating that PD-1+ cells were the main influencing factor leading to the heterogeneity in the TIME across different regions of tumor no. 7.

RNA-seq results showed that in 8 of 10 tumors (80%; nos. 1, 2, 4, 5, 6, 7, 8, and 10), the distribution range of the Mahalanobis distance was narrow for different tumor regions, while in 2 of 10 tumors (20%; nos. 3 and 9), the distribution range of the Mahalanobis distance was wide for different tumor regions (Fig. 5A). Furthermore, in 8 of 10 tumors (80%; nos. 1, 2, 4, 5, 6, 8, 9, and 10), the Mahalanobis distance for all tumor regions was 0–22 (critical value), and in 2 of 10 tumors (20%; nos. 3 and 7), the Mahalanobis distance for some tumor regions exceeded 22 (critical value) (Fig. 5A). These data indicated that for most tumors (80%), the TIME was similar for different tumor regions. Similarly, for tumor no. 3, the distribution range of the Mahalanobis distance in different regions was narrowed, whereas the Mahalanobis distance in all regions was decreased to less than 22 (critical value) after removing CD20+ B cells (Fig. 5B), indicating that CD20+ B cells were the main influencing factor leading to the heterogeneity in TIME among different regions of the tumor no. 3. However, for tumor No. 7, the Mahalanobis distance for some regions did not decrease to less than 22 (critical value), regardless of removing the markers (Fig. 5B), indicating that the TIME was highly heterogeneous for different regions of the tumor no. 7. In addition, for tumor no. 9, regardless of removing any of the markers, the distribution range of the Mahalanobis distance was not narrowed for different regions, although they were all within 22 (critical value) (Fig. 5B).

Fig. 5
figure 5

Similarity of the TIME across different tumor regions using RNA-seq data. A Mahalanobis distance of all immune markers (including eight immune cells and TLS). B Mahalanobis distance of all immune markers except one. Each dot denotes the Mahalanobis distance of a region. Dotted lines denote the critical value of Mahalanobis distance. TIME, tumor immune microenvironment; RNA-seq, RNA sequencing; and TLS, tertiary lymphatic structure

Full size image

As shown in Supplementary Figure S7, there was a strong correlation between IHC and RNA-seq data for the evaluation of TIME. For CD4+ T cells, CD8+ T cells, CD20+ B cells, LAMP3+ DCs, and PD-1+ cells, RNA-seq was positively correlated with IHC (Supplementary Figure S7). In addition, the IHC results showed that in 7 of 10 tumors (70%; nos. 1, 2, 4, 5, 6, 9, and 10), the distribution range of the Mahalanobis distance was narrow for different tumor regions, and in 3 of 10 tumors (30%; nos. 3, 7, and 8), the distribution range of the Mahalanobis distance was wide for different tumor regions (Supplementary Figure S8A). RNA-seq results showed that in 8 of 10 tumors (80%; nos. 1, 2, 4, 5, 6, 7, 8, and 10), the distribution range of the Mahalanobis distance was narrow for different tumor regions, and in 2 of 10 tumors (20%; Nos. 3 and 9), the distribution range of the Mahalanobis distance was wide for different tumor regions(Supplementary Figure S8B). Moreover, both IHC and RNA-seq results showed that in 8 of 10 tumors (80%; nos. 1, 2, 4, 5, 6, 8, 9, and 10), the Mahalanobis distance in all tumor regions was within 22 (critical value), and in 2 of 10 tumors (20%; Nos. 3 and 7), the Mahalanobis distance for some tumor regions was more than 22 (critical value) (Supplementary Figure S8A, B). These data indicated good correlation and consistency between IHC and RNA-seq data for the evaluation of TIME in DDRPLS.

Interestingly, for regions used for RNA-seq detection (T1, T3, and T5), IHC results showed that in 9 of 10 tumors (90%; nos. 1, 2, 3, 4, 5, 6, 7, 9, and 10), the distribution range of the Mahalanobis distance was narrow for different tumor regions, and in 1 of 10 tumors (10%; no. 8), the distribution range of the Mahalanobis distance was wide for different tumor regions(Supplementary Figure S8C). RNA-seq results showed that in 8 of 10 tumors (80%; nos. 1, 2, 4, 5, 6, 7, 8, and 10), the distribution range of the Mahalanobis distance was narrow for different tumor regions, and in 2 of 10 tumors (20%; Nos. 3 and 9), the distribution range of the Mahalanobis distance was wide for different tumor regions (Supplementary Figure S8B). Moreover, both IHC and RNA-seq results showed that in 8 of 10 tumors (80%; nos. 1, 2, 4, 5, 6, 8, 9, and 10), the Mahalanobis distance in all tumor regions was within 22 (critical value), and in 2 of 10 tumors (20%; Nos. 3 and 7), the Mahalanobis distance in some tumor regions was more than 22 (critical value) (Supplementary Figure S8B, C). These data also indicated very good consistency between IHC and RNA-seq data for evaluating the TIME of DDRPLS.

Characteristics of related pathways in DDRPLS

The abovementioned results showed very good correlation and consistency between IHC and RNA-seq data for evaluating the TIME of DDRLPS. Moreover, RNA-seq data were more comprehensive and extensive than IHC data. Thus, we used RNA-seq data alone to analyze immune-related pathways in the TIME of DDRPLS, including antigen presentation, cell adhesion, co-stimulator, co-inhibitor, cytokine, receptor, and ligand [22] (Fig. 6, Supplementary Figures S9 and S10). The results showed that immune-related pathways in different tumor regions were also similar (Fig. 6, Supplementary Figure S10).

Fig. 6
figure 6

Similarity of the immune-related pathways in different tumor regions using RNA-seq data alone. Immune-related pathways in different tumor regions were similar. Each dot denotes the Mahalanobis distance of a region. Dotted lines denote the critical value of Mahalanobis distance. RNA-seq, RNA sequencing

Full size image

Gene Ontology (GO) enrichment analysis of DEGs in the tumor and corresponding paratumor tissues detected via RNA-seq showed that, compared to paratumor tissues, DEGs in DDRPLS were significantly enriched in the intracellular part, intracellular, cellular process, cell part, cell, binding, and single organism process functions, with cell part and cell functions exhibiting especially significant enrichment in almost all DDRPLS (Supplementary Figure S11). These findings suggest that DDRPLS may undergo malignant transformation primarily through the regulation of cell parts and cell functions.

In addition, further Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs in the tumor and corresponding paratumor tissues showed that, compared to paratumor tissues, DEGs in DDRPLS were mainly enriched in pathways such as ribosomes, metabolic pathway, PI3K-Akt signaling pathway, cell adhesion molecules (CAMs), phagosome, and pathways in cancer (Fig. 7). The metabolic pathway and PI3K-Akt signaling pathway exhibited especially significant enrichment in almost all DDRPLS (Fig. 7). These findings suggest that DDRPLS may regulate cell functions primarily by regulating the metabolic and PI3K-Akt signaling pathways, leading to malignant transformation.

Fig. 7
figure 7

KEGG pathway enrichment analysis of DEGs in DDRPLS. Compared to paratumor tissues, DEGs in DDRPLS were significantly enriched in several pathways, such as the ribosome, metabolic, PI3K-Akt signaling, cell adhesion molecules (CAMs), phagosome, and cancer pathways, with the most significant enrichment in the metabolic pathway and the PI3K-Akt signaling pathway. KEGG, Kyoto Encyclopedia of Genes and Genomes; DEGs, differentially expressed genes; DDRPLS, dedifferentiated retroperitoneal liposarcoma

Full size image

Discussion

In this study, we showed that in DDRPLS, the number of TILs was very low, and CD4+ T cells were the most common, while Foxp3+ Tregs were the least common. DEGs were significantly enriched in cell part and cell functions, as well as the metabolic and PI3K-Akt signaling pathways. For most tumors (70–80%), the TIME in different tumor regions was similar.

Studies evaluating the TIME have mainly focused on inter-patient or inter-tumor heterogeneity and the relationship between the TIME and clinicopathological features of patients [23,24,25,26,27,28,29,30,31,32,33]. Few studies have focused on intratumor similarity [34, 35], especially for very large tumors such as RPLS. The PD-L1 expression [11,12,13] varies across different regions of some tumors, which may lead to misclassification of some PD-L1–positive tumors as PD-L1–negative tumors, hindering possible treatment with PD-1/PD-L1. Therefore, in both clinical and experimental research, we often need to obtain tumor specimens to reliably reflect the TIME of the whole tumor, especially for large tumors, such as RPLS. Therefore, to overcome this challenge and provide guidance for clinical diagnosis and treatment, the characteristics of the TIME in different tumor regions were investigated. The results showed that for most tumors (70–80%), the TIME was similar in different tumor regions. In other words, for most tumors (70–80%), the obtained specimen can reliably represent the TIME of the whole tumor, regardless of the region it is obtained from. Therefore, for most tumors (70–80%), the TIME of any tumor region can reliably represent the whole tumor.

However, we found two tumors (nos. 3 and 7) using IHC and two tumors (nos. 3 and 7) using RNA-seq that displayed significant heterogeneity in different tumor regions, for which the TIME in any tumor region could not reliably represent the whole tumor. We thought that the heterogeneity had little to do with tumor size or type, but was related to the inherent tumor characteristics. These results are supported by a study that demonstrated that different tumor regions might undergo different stages of immune-editing of neoantigen-harboring cancer cells [36].

In addition, our results showed that the TIME of DDRPLS contained very few TILs, which is consistent with a previous report [16]. This might be related to the large volume of tumor cells, as they occupy a large space in the tumor, especially well-differentiated liposarcoma (data not shown), which might explain why many reports have shown that tumor immunotherapy was ineffective for well-differentiated liposarcoma [37, 38]. These findings also indicate that it is essential to develop strategies to increase immune cells for tumor immunotherapy of RPLS. Similarly, TILs in multiple tumors are reduced because many of these tumors are partially dedifferentiated (Table 1) with a large number of tumor cells.

Moreover, we found that the number of TLSs was associated with the number of CD4+T cells, CD8+ T cells, CD20+ B cells, and CD68+ macrophages. In other words, in tumors with a greater number of TLSs, the density of CD4+T cells, CD8+ T cells, CD20+ B cells, and CD68+ macrophages were also higher. These phenomena were mainly related to strong anti-tumor immune responses [21, 39]. In addition, compared to multiple and large tumors, the number of TLSs was greater and the density of CD4+T cells, CD8+ T cells, CD20+ B cells, and CD68+ macrophages was higher in single and small tumors. These findings indicate that, in patients with high levels of TLSs, CD4+T cells, CD8+ T cells, CD20+ B cells, and CD68+ macrophages, tumor growth may be inhibited.

Furthermore, to avoid the unreliability of a single method, such as IHC alone, we evaluated the correlation and consistency between IHC and RNA-seq for evaluating the TIME. Interestingly, our results showed a very good correlation and consistency between IHC and RNA-seq for evaluating the TIME. Moreover, compared to the IHC data, RNA-seq data were more comprehensive and extensive. These findings indicate that the use of RNA-seq data alone to evaluate the TIME is feasible and comprehensive.

Our study also had some limitations. First, the sample size was relatively small, although it was consistent with statistical principles and the sample size of many previous studies [40,41,42]. Therefore, it is necessary to expand the sample size in future studies. Second, conventional IHC cannot fully characterize the functional status of immune cells and fully detect the components of TIME, such as naive, effector, memory, or exhausted T cells, M1 or M2 polarized macrophages, and granulocytic or monocytic myeloid-derived suppressor cells. Ideal evaluation methods applicable to FFPE specimens include multiplexed fluorescent IHC [35] and mass spectrometry-based multiplexed ion beam imaging [43, 44]. However, they are not routinely used in the clinical setting and are very expensive. The correlation and consistency between RNA-seq and IHC for evaluating the TIME were very good, and RNA-seq data were very comprehensive and extensive. Therefore, RNA-seq can be used alone to evaluate the TIME in a future larger study. Third, this study suggested that DDRPLS may regulate cell functions by modulating the metabolic and PI3K-Akt signaling pathways to promote its malignant behavior. However, this needs to be further investigated in the future.

Conclusions

In this study, very few TILs were found in DDRPLS. DEGs were significantly enriched in cell part and cell functions, as well as the metabolic and PI3K-Akt signaling pathways. For most tumors (70–80%), the TIME in different tumor regions was similar. Therefore, the TIME in any tumor region could reliably represent the whole tumor. DDRPLS may regulate cell functions by modulating the metabolic and PI3K-Akt signaling pathways to promote its malignant behavior. This is important to consider when obtaining a tumor specimen to reliably represent the TIME of the whole tumor in clinical or experimental research. Furthermore, these findings may provide guidance for the clinical diagnosis and treatment, as well as for research into the mechanism of RPLS.

Availability of data and materials

The data and materials of this study are available from the corresponding author and/or the first author on reasonable request.

Abbreviations

RPLS:

Retroperitoneal liposarcoma

TIME:

Tumor immune microenvironment

DDRPLS:

Dedifferentiated retroperitoneal liposarcoma

TILs:

Tumor-infiltrating lymphocytes

IHC:

Immunohistochemistry

RNA-seq:

RNA sequencing

DEGs:

Differentially expressed genes

FFPE:

Formalin-fixed and paraffin-embedded

TLS:

Tertiary lymphoid structure

DCs:

Dendritic cells

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

References

  1. Xue G, Wang Z, Li C, Lv A, Tian X, Wu J, Qiu H, Hao C. A novel nomogram for predicting local recurrence-free survival after surgical resection for retroperitoneal liposarcoma from a Chinese tertiary cancer center. Int J Clin Oncol. 2021;26:145–53.

    Article  PubMed  Google Scholar 

  2. Yan L, Wang Z, Cui C, Guan X, Dong B, Zhao M, Wu J, Tian X, Hao C. Comprehensive immune characterization and T-cell receptor repertoire heterogeneity of retroperitoneal liposarcoma. Cancer Sci. 2019;110:3038–48.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Yan Y, Xia S, Teng D, Hu S, Li S, Wang Y, Du X, Li R. Resection outcomes for primary and local recurrent retroperitoneal liposarcoma patients. Ann Transl Med. 2020;8:1450.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Paloyo SR, Ramirez AD, David-Paloyo FP, Dofitas RB. Wide excision of a retroperitoneal liposarcoma with en bloc ureterectomy and renal salvage by autotransplantation. Case Rep Transplant. 2019;2019:9725169.

    PubMed  PubMed Central  Google Scholar 

  5. Matsunaga T, Yamaguchi T, Kaida S, Takebayashi K, Murata S, Shimizu T, Sonoda H, Miyake T, Ueki T, Iida H, et al. Study of eight cases of retroperitoneal liposarcoma. Gan To Kagaku Ryoho. 2019;46:2104–6.

    PubMed  CAS  Google Scholar 

  6. Hao M, Hou S, Li W, Li K, Xue L, Hu Q, Zhu L, Chen Y, Sun H, Ju C, Zhang C: Combination of metabolic intervention and T cell therapy enhances solid tumor immunotherapy. Sci Transl Med 2020, 12.

  7. Perets R, Bar J, Rasco DW, Ahn MJ, Yoh K, Kim DW, Nagrial A, Satouchi M, Lee DH, Spigel DR, et al: Safety and efficacy of quavonlimab, a novel anti-CTLA-4 antibody (MK-1308), in combination with pembrolizumab in first-line advanced non-small-cell lung cancer. Ann Oncol 2020.

  8. Matull J, Livingstone E, Wetter A, Zimmer L, Zaremba A, Lahner H, Schadendorf D, Ugurel S. Durable complete response in a melanoma patient with unknown primary, associated with sequential and severe multi-organ toxicity after a single dose of CTLA-4 Plus PD-1 blockade: a case report. Front Oncol. 2020;10:592609.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Pecora A, Halpern S, Weber M, Paleoudis EG, Panush D, Patterson F, Toretsky J. Rapid and complete response to combination anti-CTLA-4 and anti-PD-1 checkpoint inhibitor therapy in a patient with stage IV refractory end-stage epithelioid sarcoma: a case report. J Immunother. 2020;43:286–90.

    Article  PubMed  Google Scholar 

  10. Keung EZ, Lazar AJ, Torres KE, Wang WL, Cormier JN, Ashleigh Guadagnolo B, Bishop AJ, Lin H, Hunt KK, Bird J, et al. Phase II study of neoadjuvant checkpoint blockade in patients with surgically resectable undifferentiated pleomorphic sarcoma and dedifferentiated liposarcoma. BMC Cancer. 2018;18:913.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Kwon MJ, Rho YS, Nam ES, Cho SJ, Park HR, Min SK, Seo J, Choe JY, Kim ES, Park B, et al. Clinical implication of programmed cell death-1 ligand-1 expression in tonsillar squamous cell carcinoma in association with intratumoral heterogeneity, human papillomavirus, and epithelial-to-mesenchymal transition. Hum Pathol. 2018;80:28–39.

    Article  PubMed  CAS  Google Scholar 

  12. Casadevall D, Clavé S, Taus Á, Hardy-Werbin M, Rocha P, Lorenzo M, Menéndez S, Salido M, Albanell J, Pijuan L, Arriola E. Heterogeneity of tumor and immune cell PD-L1 expression and lymphocyte counts in surgical NSCLC samples. Clin Lung Cancer. 2017;18:682-691.e685.

    Article  PubMed  CAS  Google Scholar 

  13. Wei Z, Fan L, Yang X, Li J, Zhan X, Ye X: PD-L1 heterogeneity in patients with non-small cell lung cancer. Asia Pac J Clin Oncol 2021.

  14. Shen YC, Hsu CL, Jeng YM, Ho MC, Ho CM, Yeh CP, Yeh CY, Hsu MC, Hu RH, Cheng AL. Reliability of a single-region sample to evaluate tumor immune microenvironment in hepatocellular carcinoma. J Hepatol. 2020;72:489–97.

    Article  PubMed  CAS  Google Scholar 

  15. Flores-Guerrero JL, Grzegorczyk MA, Connelly MA, Garcia E, Navis G, Dullaart RPF, Bakker SJL. Mahalanobis distance, a novel statistical proxy of homeostasis loss is longitudinally associated with risk of type 2 diabetes. EBioMedicine. 2021;71:103550.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Yan L, Wang Z, Cui C, Guan X, Dong B, Zhao M, Wu J, Tian X, Hao C. Comprehensive immune characterization and T-cell receptor repertoire heterogeneity of retroperitoneal liposarcoma. Cancer Sci. 2019;110:3038–48.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Schumacher TN, Thommen DS. Tertiary lymphoid structures in cancer. Science. 2022;375:eabf9419.

    Article  PubMed  CAS  Google Scholar 

  18. Domblides C, Rochefort J, Riffard C, Panouillot M, Lescaille G, Teillaud JL, Mateo V, Dieu-Nosjean MC. Tumor-associated tertiary lymphoid structures: from basic and clinical knowledge to therapeutic manipulation. Front Immunol. 2021;12:698604.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Mori T, Tanaka H, Suzuki S, Deguchi S, Yamakoshi Y, Yoshii M, Miki Y, Tamura T, Toyokawa T, Lee S, et al. Tertiary lymphoid structures show infiltration of effective tumor-resident T cells in gastric cancer. Cancer Sci. 2021;112:1746–57.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Horeweg N, Workel HH, Loiero D, Church DN, Vermij L, Léon-Castillo A, Krog RT, de Boer SM, Nout RA, Powell ME, et al. Tertiary lymphoid structures critical for prognosis in endometrial cancer patients. Nat Commun. 2022;13:1373.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Li H, Liu H, Fu H, Li J, Xu L, Wang G, Wu H. Peritumoral tertiary lymphoid structures correlate with protective immunity and improved prognosis in patients with hepatocellular carcinoma. Front Immunol. 2021;12:648812.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, Porta-Pardo E, Gao GF, Plaisier CL, Eddy JA, et al. The immune landscape of cancer. Immunity. 2018;48:812-830.e814.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Brunner SM, Itzel T, Rubner C, Kesselring R, Griesshammer E, Evert M, Teufel A, Schlitt HJ, Fichtner-Feigl S. Tumor-infiltrating B cells producing antitumor active immunoglobulins in resected HCC prolong patient survival. Oncotarget. 2017;8:71002–11.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Wang Q, Luan W, Warren L, Fiel MI, Blank S, Kadri H, Mandeli J, Hiotis SP. Prognostic role of immune cells in hepatitis B-associated hepatocellular carcinoma following surgical resection depends on their localization and tumor size. J Immunother. 2016;39:36–44.

    Article  PubMed  CAS  Google Scholar 

  25. Shu QH, Ge YS, Ma HX, Gao XQ, Pan JJ, Liu D, Xu GL, Ma JL, Jia WD. Prognostic value of polarized macrophages in patients with hepatocellular carcinoma after curative resection. J Cell Mol Med. 2016;20:1024–35.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Dong P, Ma L, Liu L, Zhao G, Zhang S, Dong L, Xue R, Chen S. CD86+/CD206+, diametrically polarized tumor-associated macrophages, predict hepatocellular carcinoma patient prognosis. Int J Mol Sci. 2016;17:320.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Ren CX, Leng RX, Fan YG, Pan HF, Li BZ, Wu CH, Wu Q, Wang NN, Xiong QR, Geng XP, Ye DQ. Intratumoral and peritumoral expression of CD68 and CD206 in hepatocellular carcinoma and their prognostic value. Oncol Rep. 2017;38:886–98.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Li X, Xing YF, Lei AH, Xiao Q, Lin ZH, Hong YF, Wu XY, Zhou J. Neutrophil count is associated with myeloid derived suppressor cell level and presents prognostic value of for hepatocellular carcinoma patients. Oncotarget. 2017;8:24380–8.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Sia D, Jiao Y, Martinez-Quetglas I, Kuchuk O, Villacorta-Martin C, Castro de Moura M, Putra J, Camprecios G, Bassaganyas L, Akers N, et al. Identification of an Immune-specific class of hepatocellular carcinoma, based on molecular features. Gastroenterology. 2017;153:812–26.

    Article  PubMed  CAS  Google Scholar 

  30. Lim CJ, Lee YH, Pan L, Lai L, Chua C, Wasser M, Lim TKH, Yeong J, Toh HC, Lee SY, et al. Multidimensional analyses reveal distinct immune microenvironment in hepatitis B virus-related hepatocellular carcinoma. Gut. 2019;68:916–27.

    Article  PubMed  CAS  Google Scholar 

  31. Zheng C, Zheng L, Yoo JK, Guo H, Zhang Y, Guo X, Kang B, Hu R, Huang JY, Zhang Q, et al. Landscape of infiltrating T cells in liver cancer revealed by single-cell sequencing. Cell. 2017;169:1342-1356.e1316.

    Article  PubMed  CAS  Google Scholar 

  32. Gabrielson A, Wu Y, Wang H, Jiang J, Kallakury B, Gatalica Z, Reddy S, Kleiner D, Fishbein T, Johnson L, et al. Intratumoral CD3 and CD8 T-cell densities associated with relapse-free survival in HCC. Cancer Immunol Res. 2016;4:419–30.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Calderaro J, Petitprez F, Becht E, Laurent A, Hirsch TZ, Rousseau B, Luciani A, Amaddeo G, Derman J, Charpy C, et al. Intra-tumoral tertiary lymphoid structures are associated with a low risk of early recurrence of hepatocellular carcinoma. J Hepatol. 2019;70:58–65.

    Article  PubMed  Google Scholar 

  34. Shi L, Zhang Y, Feng L, Wang L, Rong W, Wu F, Wu J, Zhang K, Cheng S. Multi-omics study revealing the complexity and spatial heterogeneity of tumor-infiltrating lymphocytes in primary liver carcinoma. Oncotarget. 2017;8:34844–57.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Kurebayashi Y, Ojima H, Tsujikawa H, Kubota N, Maehara J, Abe Y, Kitago M, Shinoda M, Kitagawa Y, Sakamoto M. Landscape of immune microenvironment in hepatocellular carcinoma and its additional impact on histological and molecular classification. Hepatology. 2018;68:1025–41.

    Article  PubMed  CAS  Google Scholar 

  36. Rosenthal R, Cadieux EL, Salgado R, Bakir MA, Moore DA, Hiley CT, Lund T, Tanić M, Reading JL, Joshi K, et al. Neoantigen-directed immune escape in lung cancer evolution. Nature. 2019;567:479–85.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Reichardt P. Soft tissue sarcomas, a look into the future: different treatments for different subtypes. Future Oncol. 2014;10:s19-27.

    Article  PubMed  CAS  Google Scholar 

  38. Ray-Coquard I, Serre D, Reichardt P, Martín-Broto J, Bauer S. Options for treating different soft tissue sarcoma subtypes. Future Oncol. 2018;14:25–49.

    Article  PubMed  CAS  Google Scholar 

  39. Zhou C, Liu Q, Xiang Y, Gou X, Li W. Role of the tumor immune microenvironment in tumor immunotherapy. Oncol Lett. 2022;23:53.

    Article  PubMed  Google Scholar 

  40. Yan Z, Ma J, Yao S, Yao Z, Wang H, Chu J, Zhao S, Liu Y. Anti-angiogenic agent combined with Anti-PD-1 immunotherapy showed activity in patients with classical Hodgkin lymphoma who have failed immunotherapy: a retrospective case report study. Front Immunol. 2021;12:727464.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Hu Z, Leet DE, Allesøe RL, Oliveira G, Li S, Luoma AM, Liu J, Forman J, Huang T, Iorgulescu JB, et al. Personal neoantigen vaccines induce persistent memory T cell responses and epitope spreading in patients with melanoma. Nat Med. 2021;27:515–25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Liu Y, He S, Wang XL, Peng W, Chen QY, Chi DM, Chen JR, Han BW, Lin GW, Li YQ, et al. Tumour heterogeneity and intercellular networks of nasopharyngeal carcinoma at single cell resolution. Nat Commun. 2021;12:741.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Keren L, Bosse M, Marquez D, Angoshtari R, Jain S, Varma S, Yang SR, Kurian A, Van Valen D, West R, et al. A structured tumor-immune microenvironment in triple negative breast cancer revealed by multiplexed ion beam imaging. Cell. 2018;174:1373-1387.e1319.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Steele NG, Carpenter ES, Kemp SB, Sirihorachai VR, The S, Delrosario L, Lazarus J, Amir ED, Gunchick V, Espinoza C, et al. Multimodal mapping of the tumor and peripheral blood immune landscape in human pancreatic cancer. Nat Cancer. 2020;1:1097–112.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

Changsheng Zhou is very grateful to his wife (Ms. Fangfang Zhang) and their parents for their support and help. The authors are very grateful to Prof. Chenghua Luo, Dr. Xiaobing Chen, and Yue Hu from the International Hospital of Peking University for their collection of some specimens, and would also like to thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.

Funding

This work was supported by grants from the Natural Science Foundation of Guizhou Province (CN) (No. Qian ke he cheng guo [2019] 4444 to Changsheng Zhou), the Beijing Medical and Health Public Welfare Foundation (No. YWJKJJHKYJJ-B184054 to Changsheng Zhou), the Key Project of the Science and Technology Ministry of China (No. 2017ZXl0203206-005–002 to Wengang Li), and the Young and Middle-Aged Backbone Talent Training Program of the Fujian Health Commission (No. 2021GGB036 to Ming Li).

Author information

Authors and Affiliations

  1. School of Medicine, Xiamen University, Xiamen, Fujian, 361102, People’s Republic of China

    Changsheng Zhou, Ming Li, Yantao Ren, Fenglin Miao, Yue Wang, Ting Wu & Wengang Li

  2. Department of Hepatobiliary Surgery, Guizhou Provincial People’s Hospital, Guiyang, Guizhou, 550002, People’s Republic of China

    Changsheng Zhou & Xin Gou

  3. Department of Hepatobiliary Surgery, Xiang’an Hospital of Xiamen University, School of Medicine, Xiamen University, Xiamen, Fujian, 361102, People’s Republic of China

    Changsheng Zhou, Ming Li, Yantao Ren, Fenglin Miao, Yue Wang, Ting Wu & Wengang Li

  4. Cancer Research Center of Xiamen University, School of Medicine, Xiamen University, Xiamen, Fujian, 361102, People’s Republic of China

    Changsheng Zhou, Ming Li, Yantao Ren, Fenglin Miao, Yue Wang, Ting Wu & Wengang Li

  5. Retroperitoneal Tumor Research Center of Oncology Chapter of Chinese Medical Association, School of Medicine, Xiamen University, Xiamen, Fujian, 361102, People’s Republic of China

    Changsheng Zhou, Ming Li, Yantao Ren, Fenglin Miao, Yue Wang, Ting Wu & Wengang Li

  6. Xiamen Medicine Research Institute, Xiamen, Fujian, 361005, People’s Republic of China

    Ming Li

Authors
  1. Changsheng Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ming Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yantao Ren
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Fenglin Miao
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yue Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Ting Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Xin Gou
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Wengang Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

C. Z contributed to the design and implementation of the experiment; the acquisition, analysis, and interpretation of data; and writing the manuscript. Y. R, F. M and Y. W contributed to the design and implementation of the experiment and the acquisition and analysis of data. M. L and T. W contributed to the design and implementation of the experiment; the acquisition and analysis of data; and revised the paper critically for important intellectual content. X. G and W. L contributed to the experiment design and data interpretation and revised the paper critically for important intellectual content.

Corresponding authors

Correspondence to Ming Li, Ting Wu, Xin Gou or Wengang Li.

Ethics declarations

Ethics approval and consent to participate

The details are shown in the materials and methods section.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table S1.

Antibody information. Supplementary Table S2. All immune markers corresponding genes. Supplementary Table S3. Correlation between immune markers and clinicopathological features of patients.

Additional file 2: Supplementary Figure 1.

Immune cell constitution in DDRPLS using RNA-seq data. For most tumors, the distribution range of immune cells was narrow for different tumor regions. Each dot denotes the relative gene expression of a sample. DDRPLS, Dedifferentiated retroperitoneal liposarcoma; RNA-seq, RNA sequencing. Supplementary Figure 2. Immune cell constitution in the same region of DDRPLS. Very few immune cells were observed in DDRPLS. A: Immune cells in the same region of DDRPLS (IHC, represented by T3). B: Immune cell density in the same region of DDRPLS (represented by T3). Each dot denotes a tumor. C: Immune cell proportion in the same region of DDRPLS (represented by T3). Each dot denotes a tumor. DDRPLS, Dedifferentiated retroperitoneal liposarcoma; IHC, Immunohistochemistry. Supplementary Figure 3. Representative IHC images of immune cells in DDRPLS. Very few immune cells were observed in DDRPLS. IHC, Immunohistochemistry; DDRPLS, Dedifferentiated retroperitoneal liposarcoma. Supplementary Figure 4. Immune cell constitution in the same region of DDRPLS. Very few immune cells were observed in DDRPLS. A: Immune cell density in the same region of DDRPLS (represented by T3). B: Immune cell proportion in the same region of DDRPLS (represented by T3). DDRPLS, Dedifferentiated retroperitoneal liposarcoma. Supplementary Figure 5. Relationship between immune markers and tumor type in DDRPLS. Fewer immune markers were detected in multiple tumorsa single tumor. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. DDRPLS, Dedifferentiated retroperitoneal liposarcoma; TLS, Tertiary lymphatic structure. Supplementary Figure 6. Relationship between immune markers and tumor size in DDRPLS. Immune markers were negatively correlated with tumor size in DDRPLS. DDRPLS, Dedifferentiated retroperitoneal liposarcoma; TLS, Tertiary lymphatic structure. Supplementary Figure 7. Correlation between IHC and RNA-seq  data for the evaluation of TIME. RNA-seq was positively correlated with IHC data for the evaluation of TIME. IHC, Immunohistochemistry; RNA-seq, RNA sequencing; and TIME, Tumor immune microenvironment. Supplementary Figure 8. Consistency between IHC and RNA-seq for evaluating TIME. IHC and RNA-seq data demonstrated consistent results for the evaluation of TIME. A: Mahalanobis distance of all regions (T1-5) used for IHC. B: Mahalanobis distance of all regions used for RNA-seq. C: Mahalanobis distance of regions (T1, T3, and T5) used for RNA-seq (detected by IHC). Each dot denotes the Mahalanobis distance of a region. Dotted lines denote the critical value of Mahalanobis distance. IHC, Immunohistochemistry; RNA-seq, RNA sequencing; and TIME, Tumor immune microenvironment. Supplementary Figure 9. Relative expression of immune-related pathways in DDRPLS. The relative expression of immune-related pathways in DDRPLS was low (represented by T3). DDRPLS, Dedifferentiated retroperitoneal liposarcoma. Supplementary Figure 10. Similarity of cytokines in different tumor regions of using RNA-seq data alone. Cytokines in different tumor regions were similar. Each dot denotes the Mahalanobis distance of a region. Dotted lines denote the critical value of Mahalanobis distance. RNA-seq, RNA sequencing. Supplementary Figure 11. GO enrichment analysis of DEGs in DDRPLS. Compared withto paratumor tissues, DEGs in DDRPLS were significantly enriched in the intracellular part, intracellular, cellular process, cell part, cell, binding, and single organism process functions, with the most significant enrichment in the cell part and cell functions. GO, Gene Ontology; DEGs, Differentially expressed genes; DDRPLS, Dedifferentiated retroperitoneal liposarcoma.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, C., Li, M., Ren, Y. et al. Immune characteristics of dedifferentiated retroperitoneal liposarcomas and the reliability of regional samples in evaluating their tumor immune microenvironments. World J Surg Onc 22, 25 (2024). https://doi.org/10.1186/s12957-023-03295-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12957-023-03295-0

Keywords

World Journal of Surgical Oncology

ISSN: 1477-7819

Contact us