Skip to main content

Combined PD-1/PD-L1 and tumor-infiltrating immune cells redefined a unique molecular subtype of high-grade serous ovarian carcinoma

Abstract

Background

High-grade serous ovarian carcinoma is highly heterogeneous, and although many studies have been conducted to identify high-grade serous ovarian carcinoma molecular subtypes that are sensitive to immunotherapy, no precise molecular subtype has been proposed to date. Immune cell infiltration and immune checkpoints are highly correlated with immunotherapy. Here, we investigated immune cell infiltration and immune checkpoint values for prognosis and precise immunotherapy for high-grade serous ovarian carcinoma based on molecular subtype classification.

Results

“High antigen-presenting cells infiltration molecular subtype of high-grade serous ovarian carcinoma” was identified in immune cell infiltration profiles. Each of the three immune cell infiltration clusters (A, B, and C) demonstrated distinct immune cell characterization, with immune cell infiltration cluster C exhibiting high antigen-presenting cell infiltration, improved prognosis, and higher sensitivity to immunotherapy. Programmed death-1/programmed death ligand 1 has a prognostic and predictive role that can help classify molecular subtypes.

Conclusions

Our findings redefined a unique molecular subtype of high-grade serous ovarian carcinoma, suggesting that high-grade serous ovarian carcinoma patients with higher antigen-presenting cell infiltration and programmed death-1/programmed death ligand 1 expression can benefit from precise immunotherapy.

Peer Review reports

Background

Epithelial ovarian carcinoma (EOC) remains the dominant cause of death among gynecological malignancies [1]. High-grade serous ovarian carcinoma (HGSOC) is mainly diagnosed at advanced stages (III/IV) with a 5-year survival rate of 45.6% and accounts for 90% of EOC cases [2, 3]. High heterogeneity and resistance to platinum chemotherapy significantly contribute to the poor prognosis of HGSOC [4]. Therefore, new molecular targets for the precise treatment of HGSOC need to be explored.

Recently, precise molecular treatment has proved effective and safe against gynecological malignant tumors, such as breast ductal carcinoma and endometrial carcinoma. In particular, the precise molecular classification of tumors contributed to treatment benefits, with the 5-year disease-free survival of breast ductal carcinoma luminal A and B subtypes being significantly higher than that of triple-negative and Her-2 positive subtypes [5, 6]. A global consensus has been reached regarding the treatment of endometrial carcinoma, stating that the POLE-ultramutated subtype of endometrial carcinoma with 98% 5-year recurrence-free survival and no need for adjuvant treatment had a better prognosis than other molecular subtypes [7, 8].

Meanwhile, numerous studies have focused on exploring the substantial molecular subtypes of HGSOC. Based on the tumor microenvironment (TME) or the infiltration of carcinoma-associated fibroblasts (CAFs) with distinct prognoses, HGSOC was categorized into three distinct subtypes: non-immune group, CAF-immune subtype, and non-CAF-immune subtype [9]. By analyzing the mRNA transcriptome, HGSOC can be divided into four transcriptomic subtypes: differentiated, immunoreactive, mesenchymal, and proliferative [10]. Additionally, by integrating transcriptomics and proteomics, HGSOC can be classified into two subtypes: high sensitivity to platinum therapy and low sensitivity to platinum therapy [11]. However, regardless of the immune signatures, transcriptome, or proteomics in the TME, the robustness, availability, and variability of these subtypes remain controversial, and it is difficult to reach a consensus on the targeted treatment for the molecular subtypes of HGSOC. Therefore, identifying the potential molecular subtypes indicative of HGSOC burden and therapies is essential.

Programmed death-1 (PD-1, CD279, PDCD1) is an immune checkpoint receptor of the CD28 family and is expressed on the surface of activated T cells, B cells, and other immune cells, especially dendritic cells (DCs). PD-1 is a receptor for programmed death ligand 1 (PD-L1, CD274, B7-H1), which is expressed on the surface of tumor cells [12]. Anti-PD-1 and its ligand PD-L1 have become critical to immunotherapy strategies in advanced malignant tumors, including melanoma, non-small cell lung carcinoma, urothelial carcinoma, renal cell carcinoma, and head and neck carcinoma, and are associated with improved outcome and survival [13,14,15,16,17]. The interaction of PD-1 and PD-L1 can prevent T cell activation by downregulating the function of tumor-infiltrating immune cells (TIICs) and weakening tumor immunogenicity [18]. In a previous study, PD-L1 expression was positive in 34.9% of all EOC subtypes and 57.4% of HGSOC subtypes [19]. In addition, PD-L1 was poorly expressed in tissue biopsies but was more prevalent in tumor-infiltrating cells than in tumor cells in HGSOC tissue sections [20, 21]. Patients with PD-L1 expression may be more sensitive to immunotherapy. Therefore, tumor immunogenicity and PD-1/PD-L1 expression are essential features for the molecular classification of HGSOC.

Notably, tumor immunogenicity can be weakened by the interaction of PD-1 and PD-L1 by impairing the functions of immune cells. However, APCs, which are vital components of TIICs, can improve immunogenicity by activating naïve T cells. This suggests that patients with a “high immunogenicity molecular subtype of HGSOC” might benefit from immunotherapy. To this end, we combined PD-1/PD-L1 with TIICs to explore the potential molecular subtypes of HGSOC.

Results

PD-1 and PD-L1 were associated with improved outcomes in HGSOC

To investigate whether PD-1/PD-L1 expression was related to clinical prognosis, clinical information was extracted from four databases, including three HGSOC databases: GSE53963 (GSE53963-HGSOC), GSE135820 (GSE135820-HGSOC), and GSE32062 (GSE32062-HGSOC) and The Cancer General Atlas (TCGA) database of OC (TCGA-OC). Based on PD-1/PD-L1 expression levels from each database, the average score was set as a cutoff value dividing patients into two groups: higher and lower expression of PD-1/PD-L1. The higher PD-1/PD-L1 expression group had significantly better overall survival (OS) in the three HGSOC databases compared with the lower expression group (p = 0.023 for PD-1 and p = 0.021 for PD-L1 in GSE5396, p = 0.016 for PD-1 and p = 0.026 for PD-L1 in GSE32062, p < 0.001 for PD-1 and p < 0.001 for PD-L1 in GSE135820, log-rank test; Fig. 1A–C). No statistically significant difference in prognosis was obtained between the higher and the lower PD-1/PD-L1 expression groups of TCGA-OV (p = 0.307 for PD-1, p = 0.065 for PD-L1, log-rank test; Fig. 1D).

Fig. 1
figure 1

Prognostic significance of PD-1/PD-L1 in HGSOC. Kaplan-Meier analysis for overall survival of patients with higher (red) or lower (blue) PD-1/PD-L1 expression in the GSE53963 (A), GSE32062 (B), GSE135820 (C), and TCGA (D) databases. The log-rank test was used to compare groups, and p < 0.05 was considered significant. HGSOC, High-grade serous ovarian carcinoma; PD-1, programmed death-1; PD-L1, programmed death ligand 1; TCGA, The Cancer General Atlas

Differentially expressed PD-1/PD-L1-related genes in HGSOC

To identify DEGs that were relevant to PD-1/PD-L1 expression, PD-1/PD-L1 was identified and screened in the GSE135820 microarray database after normalizing the chip results. Based on the PD-1/PD-L1 expression level, every sample from the microarray database was divided into two groups: higher and lower PD-1/PD-L1 expression groups. Keeping the cutoff value strict, genes with |logFC| > 0.5 of their expression between the two groups were considered DEGs and were depicted in a volcano plot. The top 20 DEGs that were conspicuously related to PD-1/PD-L1 are illustrated in a heatmap (Fig. 2A, B). Most of these DEGs were associated with the encoding of antigens on lymphatic cell membranes. The “Corrplot” R package was used to further investigate the correlations of DEGs. PD-1 was positively related to CD8A, CTLA4, and CD3D mRNA expression, and PD-L1 was significantly related to CTLA4, CD38, and CXCL9 mRNA expression (Fig. 2C, D).

Fig. 2
figure 2

DEGs were compared between relatively higher expression and lower expression groups of PD-1/PD-L1 by analyzing the microarray database GSE135820. Heatmaps show DEGs that were related to PD-1 (A) or PD-L1 (B) in the two groups. Red indicates DEGs that were positively associated with PD-1/PD-L1, white indicates DEGs that were not related to PD-1/P-L1, and blue indicates DEGs that were negatively related to PD-1/P-L1. Volcano plot demonstrates different scattering of genes between high and low expression cluster groups of PD-1 or PD-L1. Genes (depicted in green or red) showed a fold change difference of > 0.5 and p < 0.05. Correlations of DEGs that were related to PD-1 (C) or PD-L1 (D) between two groups are shown in the correlation heatmap. Red indicates a positive correlation (r = 1), white indicates no correlation (r = 0), and green indicates a negative correlation (r = −1). DEGs, differentially expressed genes; PD-1, programmed death-1; PD-L1, programmed death ligand 1

DEGs were closely associated with lymphocyte differentiation

Functional enrichment analysis of DEGs was performed to identify a heterogenic intrinsic trait of HGSOC. PD-1/PD-L1 related DEG function analysis was performed using GO, KEGG, and PPI online tools (Fig. 3). In GO analysis, lymphocyte-associated pathways, including T cell and B cell, were the most enriched in both independent experimental groups of PD-1 and PD-L1 (Fig. 3A). In KEGG analysis, DEGs were mainly enriched in proinflammatory chemokine and cytokine-associated pathways, including chemokine signaling pathway, cytokine-cytokine receptor interaction, and toll-like receptor interaction (Fig. 3B). The online database “STRING” was further used to create a visible network to elucidate the structural and functional signatures of proteins (Fig. 3C). All pathways generated by KEGG, GO, and PPI analyses demonstrated lymphocyte differentiation.

Fig. 3
figure 3

Functional analysis of DEGs associated with PD-1/PD-L1. A GO analysis of DEGs. B KEGG pathway analysis of DEGs. C PPI analysis of DEGs. The plot sizes indicate gene counts enriched in the pathway and the color depth shows the p-value from low (red) to high (blue) level. Significance was set at p < 0.05. DEGs, differentially expressed genes; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PPI, protein–protein interaction; PD-1, programmed death-1; PD-L1, programmed death ligand 1

Higher APCs infiltrating fraction in HGSOC

Samples with CIBERSORT p ≥ 0.05 were eliminated, and infiltrating proportions of 22 immune cells were depicted and calculated in each database (Table 1). The fractions of activated memory CD4 T cells, regulatory T cells, and delta gamma T cells were excluded from the subsequent analysis as they had zero proportion in normal tissue samples. The organization of LM22 infiltrating immune cells in 23 normal and 239 HGSOC samples is summarized in Fig. 4A and B. Compared to normal tissues, infiltrating immune cell patterns of HGSOC are depicted in a heatmap (Fig. 4C), and the majority of LM22 infiltrating immune cell fractions were predominant in HGSOC groups (Fig. 4D). The proportions of memory B cells, naïve CD4 T cells, regulatory T cells, gamma delta T cells, resting NK cells, M1 macrophages, activated DCs, resting DCs, and eosinophils were significantly increased in HGSOC, whereas plasma cells, resting memory CD4 cells, follicular helper T cells, activated NK cells, monocytes, M2 macrophages, resting mast cells, and activated mast cells decreased compared to normal tissues. Compared to other EOC histological subtypes, DCs, M1 macrophages, and B cells had a notably higher proportion in LM22 immune cells in HGSOC, whereas CD8 T cells showed the opposite trend (Table 1). Relationships among LM22 immune cells in normal tissues and HGSOC were investigated. For instance, M1 macrophages were highly positively associated with M0 macrophages and strongly negatively related to activated DCs in normal tissues. In the GSE135850-HGSOC database, M1 macrophages were highly positively related to resting DCs and negatively related to activated DCs in the immune phenotype profiles (Fig. 4E). Similar phenomena were found in GES32062-HGSOC, GSE53963-HGSOC, and serous OC databases (Fig. S1A, B). The relationship and patterns of immune infiltrating cells in other EOC histological types are shown in the supplemental materials (Fig. S1C–F). HGSOC has a unique immune infiltrating pattern, which is significantly different from that of other EOC subtypes.

Table 1 Proportion of immune infiltrating cells in the HGSOC (GSE135820, GSE53963, GSE32062), normal, serous, endometrioid, clear cell, and TCGA-OC databases. A difference in the proportions of immune-infiltrating cells was observed between HGSOC-135820 and the other databases. (*p < 0.05, **p < 0.01, ***p < 0.001)
Fig. 4
figure 4

Assessment of the proportions of LM22 infiltrating immune cells. A, B Bar plots exhibit the immune cell proportion differences in normal tissues and GSE135820-HGSOC. C, D Differences in LM22 immune cell fractions between HGSOC and normal tissues are demonstrated in the heat map and Violin plot. E Relevance of LM22 infiltrating immune cells in normal tissues and GSE135820-HGSOC are shown in the correlation matrix. HGSOC, High-grade serous ovarian carcinoma; LM22, leukocyte gene signature matrix

A potent immunogenicity ICI subtype of HGSOC was recognized

Based on CDF and delta curve analysis, immune cells could be classified into three optimal cluster groups (Fig. 5A). Meanwhile, 2873 HGSOC samples with CIBERSORT p < 0.05 were clustered into three subtypes: cluster 1/A (n = 1961), cluster 2/B (n = 737), and cluster 3/C (n = 175). In particular, cluster 3/C was a distinct portion (Table S4), with better outcomes, whereas ICI clusters A and B had poor outcomes (p = 0.012, log-rank test, Fig. 5B). Combined with the analysis of two immune checkpoints, PD-1 and PD-L1, in each group, ICI cluster C had a notably higher PD-1/PD-L1 expression level than ICI clusters A and B (p < 0.001, Kruskal-Wallis test, Fig. 5C). The distribution of immune cells in different cluster types is shown in the heatmap (Fig. 5D). Subsets of LM22 TIICs showed significant differences among the three cluster types (Fig. 5E); ICI cluster C was notable for a significantly higher proportion of naïve B cells, plasma cells, resting memory CD4 memory T cells, follicular helper T cells, regulatory T cells, resting NK cells, M1 macrophages, activated NK cells, activated DCs, and eosinophil infiltration. Cluster B was characterized by high levels of CD8 T cells, naïve CD4 T cells, monocytes, M1 macrophages, activated mast cells, and neutrophils. Cluster A had high levels of memory B cells, activated memory CD4 T cells, gamma delta T cells, and M2 macrophages. We concluded that ICI cluster C HGSOC patients with improved prognosis had higher APC infiltration and PD-1/PD-L1 expression levels.

Fig. 5
figure 5

Differential clinicopathological features and survival of HGSOC in cluster A/B/C subtypes in GSE135820-HGSOC. A Consensus clustering matrix for k = 3. B Kaplan-Meier analysis of OS for patients with HGSOC in three clusters. C Differences in PD-1/PD-L1expression among three ICI clusters. D Unsupervised clustering of TIICs in GSE135820-HGSOC. Rows represent TIICs, and columns represent samples. E The fraction of TIICs in three ICI clusters. (Kruskal-Wallis test, *p < 0.05, **p < 0.01, ***p < 0.001, and ns p > 0.05). HGSOC, High-grade serous ovarian carcinoma; ICI, immune cell infiltration; OS, overall survival; TIICs, tumor-infiltrating immune cells; PD-1, programmed death-1; PD-L1, programmed death ligand 1

Combined PD-1/PD-L1 with DCs identified a potent immunogenic molecular subtype of HGSOC

To investigate infiltrating immune cells that could be associated with PD-1/PD-L1 expression level in HGSOC, the relationship between PD-1/PD-L1 and immune components in the TME was visualized (Fig. 6). Immune scores and several degrees of ICI (including mainly M1 macrophages, M2 macrophages, and resting DCs) seemed relevant to PD-1/PD-L1 expression levels (Fig. 6A, B). The relationship between the three immune cells and PD-1/PD-L1 was further investigated by Pearson correlation analysis. PD-1/PD-L1 expression was positively correlated to M1 macrophages in both TCGA-OC and GSE135820-HGSOC. Moreover, PD-1/PD-L1 was positively associated with resting DCs and negatively associated with M2 macrophages in GSE135820-HGSOC (Fig. 6D). However, no solid conclusion was made regarding the correlation of PD-1/PD-L1 with resting DCs and M2 macrophages in TCGA-OC (Fig. 6C). Similar results were found in the GSE53963 database of HGSOC (Fig. S2). We constructed a three factor/two level orthogonal experimental model to recalculate the correlation between overall survival and APC infiltrating cells in 239 HGSOC patients using SPSS software (Table S5). The overall survival of 239 HGSOC patients was related to the number of M1 macrophages, dendritic cells, and B cells (p < 0.01). The overall survival of 239 HGSOC patients was not affected by the three lines of APC infiltrating cells (p > 0.05).

Fig. 6
figure 6

Correlation between PD-1/PD-L1 expression level and immune cell subsets. Heatmap indicates the correlation of PD-1 (A) and PD-L1 (B) with 22 immune cell subpopulations in GSE135820-HGSOC. Pearson coefficient correlation revealed the correlations between PD-1/PD-L1 expression levels and three immune cell subsets in TCGA-OC (C) and GSE135820-HGSOC (D). HGSOC, High-grade serous ovarian carcinoma; OC, ovarian carcinoma; PD-1, programmed death-1; PD-L1, programmed death ligand 1

Discussion

PD-L1 is highly expressed in several advanced carcinomas and is a favorable prognostic candidate [22,23,24]. Several studies have demonstrated that higher expression of PD-L1 indicates a favorable prognosis and restored immunogenicity [25, 26], whereas some reports showed that PD-L1 expression is not associated with OS [27, 28]. The lack of adequate clinical samples and distinct experimental approaches can explain the ambiguity of whether high PD-L1 expression in HGSOC correlates with improved prognosis. In our study, we analyzed clinical data and PD-1/PD-L1 expression in a total of 4475 samples and reported that PD-1 and PD-L1 were both associated with a better prognosis of HGSOC. Compared with previous studies, our results are reproducible and predictable as we analyzed both PD-1 and PD-L1 expression in association with HGSOC survival using data from three independent HGSOC databases (GSE53963, GSE32062, and GSE135820). The high expression of PD-1 and PD-L1 in HGSOC in our study suggests that PD-1/PD-L1 expression could have a prognostic value in HGSOC.

PD-1/PD-L1 expression in HGSOC varies across individuals. To further clarify the factors influencing PD-1/PD-L1 expression in different HGSOC patients, some DEGs, such as cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4), CXCL9, and CCL5, underwent bioinformatics analysis in our study. CTLA-4 is a common immune checkpoint inhibitor that prevents T cell activation by inhibiting the interaction of T cells with APCs [29]. T cell infiltration in the TME requires the co-expression and interaction of CCL5 and CXCL9 [30], and high expression of CXCL9 and PD-L1 is associated with longer OS [31]. The expression of DEGs regulates the infiltration ability of immune cells. Our study also showed that DEGs are associated with immune and inflammatory responses, chemotaxis and activation of lymphocytes, and cell adhesion. Collectively, these DEGs may affect the immune process by altering PD-1/PD-L1 expression, and the underlying component of immune cells could be considered a potential factor affecting PD-1/PD-L1 expression in HGSOC.

Immune cells in the TME during tumor immune response are not well understood; thus, the online tools CIBERSORT and ESTIMATE were applied to demonstrate the infiltrating patterns of immune cells and assess the value of the immune/stromal score of the immune response in various carcinomas [32]. Patients with higher immune/stromal scores had poor disease-free survival and a higher pathological T stage in prostate carcinoma [33]. Immunological components included prognostically relevant biomarkers of immune responses in HGSOC patients, and PD-1+ TIM-3+ CD8+ T cells demonstrated exhausted function and were associated with poor prognosis of HGSOC [34]. Evidence of a single EOC histological subtype regarding immune infiltration patterns has been poorly reported. A previous study used LM22 to evaluate the feasibility of leukocyte deconvolution from bulk tumors, including 547 genes that distinguish 22 human hematopoietic cell phenotypes [35]. Hence, investigating the patterns of LM22 immune cells is needed to reveal immune cell infiltration characteristics in HGSOC and other histological subtypes of EOC. Lower infiltration of T cells and M2 macrophages could be found in the TME of HGSOC. In addition, in the present study, a particularly higher APC (DCs, B cells, and M1 macrophages) infiltrating pattern in HGSOC was demonstrated compared with other histological types of EOC. Mature DCs are associated with favorable immune infiltrating patterns and improved outcomes in HGSOC patients [36]. Therefore, we clarified that HGSOC with a high APC infiltration pattern might have favorable clinical outcomes based on variations in ICI proportion among different EOC subtypes.

The ICI cluster classification method has been applied to a broad range of molecular subtypes of carcinoma because of its general usefulness in immunotherapy [37, 38]. Triple-negative breast carcinoma (TNBC) was identified in three immune molecular subtypes (immunity high, immunity medium, and immunity low) by ICI cluster classification. Identifying TNBC subtypes has potential clinical significance in TNBC treatment [39]. A special molecular subtype of HGSOC in patients with a high APC infiltrating pattern, which may be more beneficial for immune therapy, was further analyzed. HGSOC was further divided into three distinct ICI clusters based on this classification. Immune cells demonstrated distinct expression fractions of three ICI cluster groups: ICI cluster C with the highest APC infiltration had the best survival prognosis, ICI cluster B with intermediate APC infiltration had medium survival prognosis, whereas ICI cluster A with the lowest APC infiltration had poor outcomes. Interestingly, ICI cluster C, with a higher proportion of DCs, as a unique immune infiltrating subtype that differs from other ICI cluster groups, was observed in our study. DCs, which are antigen-presenting cells, can enhance tumor immunogenicity by identifying, engulfing, processing, and presenting tumor-associated antigens to activate naïve T cells [40]. Potent immunogenicity can contribute to an effective response to immunotherapy [41, 42]. Therefore, HGSOC patients in ICI cluster C with the highest APC infiltration might be more sensitive to immunotherapy. These results demonstrate the rationality of APC infiltration as the basis for immune molecular subtype classification.

To further investigate the association of PD-1/PD-L1 expression with immune cells in HGSOC, three ICI clusters showed significantly distinct PD-1/PD-L1 expression. Resting DCs had a significant positive correlation with PD-1/PD-L1 expression in HGSOC. Considering that DCs can improve immunogenicity during antigen presentation, the degree of infiltration of DCs and the potency of immunogenicity can be indicated by PD-1/PD-L1 expression. In contrast, clinically relevant immunosuppression can be demonstrated by M2 macrophages in the TME [43]. In our study, PD-1/ PD-L1 expression had a strong positive correlation with M1 macrophages and a negative correlation with M2 macrophages. High M1 and low M2 macrophage infiltration of OC is associated with improved survival [44,45,46]. Hence, the intrinsic interaction between PD-1/PD-L1 and DCs could improve the prognosis of HGSOC patients. Our findings demonstrate that PD-1/PD-L1 could be a reference biomarker for molecular subtype classification of HGSOC.

Conclusions

Patients with high APC infiltration may have potent immunogenicity and a satisfactory prognosis. Comprehensive assessment of LM22 ICI models of HGSOC might facilitate the exploration of new molecular targets for precise immunotherapy. We elucidated that the “high APC infiltration molecular subtype of HGSOC” has high immunogenicity and can greatly benefit from anti-PD-1/PD-L1 immunotherapy. In addition, PD-1/PD-L1 expression level can be considered a predictor of the infiltrating level of APCs. However, the correlation between the molecular subtype of HGSOC with high APC infiltration and the expression of PD-1, PD-L1, survival, and prognosis of HGSOC patients requires further validation in independent tissue samples.

Methods

Data source and data processing

The RNA-seq transcriptome data and the corresponding clinical data of OC were downloaded from TCGA (https://portal.gdc.carcinoma.gov/) and Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) data portals. Gene expression data of 88 normal ovarian samples were downloaded from the Genotype-Tissue Expression (GTEx) (https://commonfund.nih.gov/GTEx/) data portal. Gene expression profiles of nine databases: GSE53963 (n = 174), GSE32060 (n = 260), GSE135820 (n = 4041), GSE17260 (n = 110), GSE73614 (n = 107), GSE73638 (n = 102), GSE68600 (n = 113), TCGA database of OC (n = 378), and GTEx accession (n = 88) were obtained. Based on different histological types of EOC in four independent databases (GSE17260, GSE73614, GSE73638, and GSE68600), three recombinational EOC histological database samples, namely serous (n = 33), endometrioid (n = 25), and clear cell OC (n = 21), were collected from the four EOC databases. The normalized data were uploaded by the authors and retrieved from the aforementioned three data portals. The “Limma” R package was used to identify differentially expressed genes (DEGs). Gene probe IDs for which an “-/NA” notation was not available were removed.

Survival analysis

Clinical survival information was collected from the GSE53963, GSE32062, GSE135820, and TCGA accessions. The clinical endpoint was OS, and this clinical information was downloaded from each dataset (Table S1). The “Survival” R package was applied to evaluate the prognostic value of PD-1 and PD-L1.

Enrichment and protein-protein network analysis

Relative DEG lists were uploaded using the online tools Database for Annotation Visualization and Integrated Discovery (DAVID, david.ncifcrf.gov/) and STRING (https://string-db.org/), and the results of Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein-protein interaction (PPI) analyses were obtained. Datasets with |log fold change | ≥0.2 and adjusted p < 0.05 were deemed meaningful enrichment pathways. The “Cluster Profiler” R package was used to combine gene annotation and gene expression analyses results [47]. The results are visualized in a circular plot.

Molecular subtypes identification and estimated fractions of the TIICs

The immune and stromal scores of each sample were assessed using the ESTIMATE algorithm of the “estimate” R package [48]. A leukocyte gene signature matrix (LM22) was designed and validated to assess the feasibility of leukocyte deconvolution from bulk tumors. It contains 547 genes that distinguish 22 human hematopoietic cell phenotypes, including seven T cell types, naïve and memory B cells, plasma cells, NK cells, and myeloid subsets [35]. The fractions of 22 immune cells from each sample were obtained by calculating the relative subsets of RNA transcripts. This approach mainly depends on the LM22. TIIC levels in EOC tissues were identified using the “CIBERSORT” R package [49]. The total proportion of 22 immune cells was one for every sample. Only samples with a CIBERSORT p < 0.05 were considered significant and subjected to subsequent analysis. The results for TIIC proportions and histological subtypes of EOC after processing and screening in each database are shown in Table S2. The consensus cluster algorithm was applied using the “ConsensuClusterPlus” R package to determine the potential clusters of TIICs [47].

Correlation between PD-1/PD-L1 and TIICs

The average proportions of LM22 immune cell subsets in each microarray dataset are shown in Table S3. PD-1/PD-L1 expression profiles in GSE135820, GSE53963, GSE32062, and TCGA-OC were obtained. The correlation between PD-1/PD-L1 and TIICs was evaluated using the R package.

Statistical analyses

OS differences were calculated using the log-rank statistic in R. Pearson correlation analysis was performed to determine the relevant LM22 immune cells in the three groups. The Wilcoxon test was used to compare immune cell infiltration levels in primary HGSOC and normal ovarian tissues. All statistical analyses were performed using R Version 4.0.2 (R-project.org), and R packages were obtained from the Bioconductor project (www.bioconductor.org). Statistical significance was set at p < 0.05.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

APCs:

Antigen-presenting cells

CAF:

Carcinoma-associated fibroblast

CDF:

Cumulative distribution functions

CTLA-4:

Cytotoxic T-lymphocyte-associated antigen 4

DC:

Dendritic cell

DEGs:

Differentially expressed genes

EC:

Endometrial carcinoma

EOC:

epithelial ovarian carcinoma

GEO:

Gene Expression Omnibus

GO:

Gene Ontology

GTEx:

Genotype-Tissue Expression

HGSOC:

High-grade serous ovarian carcinoma

ICI:

Immune cell infiltration

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LM22:

Leukocyte gene signature matrix

OC:

Ovarian carcinoma

OS:

Overall survival

PD-1:

Programmed death-1

PD-L1:

Programmed death ligand 1

PPI:

Protein-protein interaction

TCGA:

The Cancer General Atlas

TIIC:

Tumor-infiltrating immune cells

TME:

Tumor microenvironment

TNBC:

Triple-negative breast carcinoma

References

  1. 1.

    Chiofalo B, Bruni S, Certelli C, Sperduti I, Baiocco E, Vizza E. Primary debulking surgery vs. interval debulking surgery for advanced ovarian cancer: review of the literature and meta-analysis. Minerva medica. 2019;110:330–40.

    Article  Google Scholar 

  2. 2.

    Bowtell D, Böhm S, Ahmed A, Aspuria P, Bast R, Beral V, et al. Rethinking ovarian cancer II: reducing mortality from high-grade serous ovarian cancer. 2015;15:668–79.

  3. 3.

    Oronsky B, Ray C, Spira A, Trepel J, Carter C, Cottrill HJM. A brief review of the management of platinum-resistant-platinum-refractory ovarian cancer. 2017;34:103.

  4. 4.

    Matulonis UA, Sood AK, Fallowfield L, Howitt BE, Sehouli J, Karlan BY. Ovarian cancer. Nature reviews Disease primers. 2016;2:16061.

    Article  Google Scholar 

  5. 5.

    Sotiriou C, Pusztai L. Gene-expression signatures in breast cancer. The New England journal of medicine. 2009;360:790–800.

    CAS  Article  Google Scholar 

  6. 6.

    Zhong Y, Xu Y, Zhou Y, Mao F, Lin Y, Guan J, et al. Breast-conserving surgery without axillary lymph node surgery or radiotherapy is safe for HER2-positive and triple negative breast cancer patients over 70 years of age. 2020;182:117–26.

  7. 7.

    Kandoth C, Schultz N, Cherniack AD, Akbani R, Liu Y, Shen H, et al. Integrated genomic characterization of endometrial carcinoma. Nature. 2013;497:67–73.

    Article  Google Scholar 

  8. 8.

    León-Castillo A, de Boer S, Powell M, Mileshkin L, Mackay H, Leary A, et al. Molecular classification of the PORTEC-3 trial for high-risk endometrial cancer: impact on prognosis and benefit from adjuvant. Therapy. 2020;38:3388–97.

    Google Scholar 

  9. 9.

    Lu X, Ji C, Jiang L, Zhu Y, Zhou Y, Meng J, et al. Tumour microenvironment-based molecular profiling reveals ideal candidates for high-grade serous ovarian cancer immunotherapy. Cell proliferation. 2021;54:e12979.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Hu Y, Taylor-Harding B, Raz Y, Haro M, Recouvreux MS, Taylan E, et al. Are Epithelial ovarian cancers of the mesenchymal subtype actually intraperitoneal metastases to the ovary? Frontiers in cell and developmental biology. 2020;8:647.

    Article  Google Scholar 

  11. 11.

    Xie H, Xu H, Hou Y, Cai Y, Rong Z, Song W, et al. Integrative prognostic subtype discovery in high-grade serous ovarian cancer. Journal of cellular biochemistry. 2019;120:18659–66.

    CAS  Article  Google Scholar 

  12. 12.

    Bardhan K, Anagnostou T, Boussiotis VA. The PD1:PD-L1/2 Pathway from Discovery to Clinical Implementation. Frontiers in immunology. 2016;7:550.

    Article  Google Scholar 

  13. 13.

    Haslam A, Prasad V. Estimation of the percentage of US patients with cancer who are eligible for and respond to checkpoint inhibitor immunotherapy drugs. JAMA network open. 2019;2:e192535.

    Article  Google Scholar 

  14. 14.

    Carbone DP, Reck M, Paz-Ares L, Creelan B, Horn L, Steins M, et al. First-line nivolumab in stage IV or recurrent non-small-cell lung cancer. The New England journal of medicine. 2017;376:2415–26.

    CAS  Article  Google Scholar 

  15. 15.

    Schachter J, Ribas A, Long GV, Arance A, Grob JJ, Mortier L, et al. Pembrolizumab versus ipilimumab for advanced melanoma: final overall survival results of a multicentre, randomised, open-label phase 3 study (KEYNOTE-006). Lancet. 2017;390:1853–62.

    CAS  Article  Google Scholar 

  16. 16.

    Bellmunt J, de Wit R, Vaughn DJ, Fradet Y, Lee JL, Fong L, et al. Pembrolizumab as second-line therapy for advanced urothelial carcinoma. The New England journal of medicine. 2017;376:1015–26.

    CAS  Article  Google Scholar 

  17. 17.

    Motzer RJ, Escudier B, McDermott DF, George S, Hammers HJ, Srinivas S, et al. Nivolumab versus everolimus in advanced renal-cell carcinoma. The New England journal of medicine. 2015;373:1803–13.

    CAS  Article  Google Scholar 

  18. 18.

    Drakes ML, Stiff PJ. Regulation of ovarian cancer prognosis by immune cells in the tumor microenvironment. Cancers. 2018;10.

  19. 19.

    Webb JR, Milne K, Kroeger DR, Nelson BH. PD-L1 expression is associated with tumor-infiltrating T cells and favorable prognosis in high-grade serous ovarian cancer. Gynecologic oncology. 2016;141:293–302.

    CAS  Article  Google Scholar 

  20. 20.

    Lara OD, Krishnan S, Wang Z, Corvigno S, Zhong Y, Lyons Y, et al. Tumor core biopsies adequately represent immune microenvironment of high-grade serous carcinoma. Scientific reports. 2019;9:17589.

    Article  Google Scholar 

  21. 21.

    Chen H, Molberg K, Strickland AL, Castrillon DH, Carrick K, Jiang Q, et al. PD-L1 expression and CD8+ tumor-infiltrating lymphocytes in different types of tubo-ovarian carcinoma and their prognostic value in high-grade serous carcinoma. The American journal of surgical pathology. 2020;44:1050–60.

    Article  Google Scholar 

  22. 22.

    Velcheti V, Schalper K, Carvajal D, Anagnostou V, Syrigos K, Sznol M, et al. Programmed death ligand-1 expression in non-small cell lung cancer. 2014;94:107–16.

  23. 23.

    Droeser R, Hirt C, Viehl C, Frey D, Nebiker C, Huber X, et al. Clinical impact of programmed cell death ligand 1 expression in colorectal cancer. 2013;49:2233–42.

  24. 24.

    Madore J, Vilain R, Menzies A, Kakavand H, Wilmott J, Hyman J, et al. PD-L1 expression in melanoma shows marked heterogeneity within and between patients: implications for anti-PD-1/PD-L1 clinical trials. 2015;28:245–53.

  25. 25.

    Darb-Esfahani S, Kunze C, Kulbe H, Sehouli J, Wienert S, Lindner J, et al. Prognostic impact of programmed cell death-1 (PD-1) and PD-ligand 1 (PD-L1) expression in cancer cells and tumor-infiltrating lymphocytes in ovarian high grade serous carcinoma. 2016;7:1486–99.

  26. 26.

    Webb J, Milne K, Kroeger D, Nelson BJG. PD-L1 expression is associated with tumor-infiltrating T cells and favorable prognosis in high-grade serous ovarian cancer. 2016;141:293–302.

  27. 27.

    Hamanishi J, Mandai M, Iwasaki M, Okazaki T, Tanaka Y, Yamaguchi K, et al. Programmed cell death 1 ligand 1 and tumor-infiltrating CD8+ T lymphocytes are prognostic factors of human ovarian cancer. 2007;104:3360–5.

  28. 28.

    Kim H, Kim J, Lee Y, Kim S, Lee J, Nam E, et al. Expression of programmed cell death ligand 1 and immune checkpoint markers in residual tumors after neoadjuvant chemotherapy for advanced high-grade serous ovarian cancer. 2018;151:414–21.

  29. 29.

    Schneider H, Downey J, Smith A, Zinselmeyer B, Rush C, Brewer J, et al. Reversal of the TCR stop signal by CTLA-4. 2006;313:1972–5.

  30. 30.

    Dangaj D, Bruand M, Grimm A, Ronet C, Barras D, Duttagupta P, et al. Cooperation between constitutive and inducible chemokines enables T cell engraftment and immune attack in solid tumors. 2019;35:885–900.e10.

    CAS  Google Scholar 

  31. 31.

    Kubon J, Sikic D, Eckstein M, Weyerer V, Stöhr R, Neumann A, et al. Analysis of CXCL9, PD1 and PD-L1 mRNA in stage T1 non-muscle invasive bladder cancer and their association with prognosis. 2020;12.

  32. 32.

    Zhu J, Liu Y, Ao H, Liu M, Zhao M, Ma JJFio. ACK1Comprehensive Analysis of the Immune Implication of Gene in Non-small Cell Lung. Cancer. 2020;10:1132.

    Google Scholar 

  33. 33.

    Zhao H, Zhang X, Shi Z, Guo B, Zhang W, He K, et al. Identification of a prognostic signature model with tumor microenvironment for predicting disease-free survival after radical prostatectomy. 2021;12:2371–84.

  34. 34.

    Fucikova J, Rakova J, Hensler M, Kasikova L, Belicova L, Hladikova K, et al. TIM-3 dictates functional orientation of the immune infiltrate in ovarian cancer. 2019;25:4820–31.

  35. 35.

    Newman A, Liu C, Green M, Gentles A, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. 2015;12:453–7.

  36. 36.

    Truxova I, Kasikova L, Hensler M, Skapa P, Laco J, Pecen L, et al. Mature dendritic cells correlate with favorable immune infiltrate and improved prognosis in ovarian carcinoma patients. 2018;6:139.

  37. 37.

    Zhang X, Shi M, Chen T, Zhang BJMN. Characterization of the immune cell infiltration landscape in head and neck squamous cell carcinoma to aid immunotherapy. 2020;22:298–309.

  38. 38.

    Yi L, Wu G, Guo L, Zou X, Huang PJMtNa. Comprehensive analysis of the PD-L1 and immune infiltrates of mA RNA methylation regulators in head and neck squamous cell carcinoma. 2020;21, 299-314.

  39. 39.

    He Y, Jiang Z, Chen C, Wang XJJoe, CR ccr. Classification of triple-negative breast cancers based on Immunogenomic profiling. 2018;37, 327.

  40. 40.

    Banchereau J, Briere F, Caux C, Davoust J, Lebecque S, Liu Y, et al. Immunobiology of dendritic cells. 2000;18, 767-811.

  41. 41.

    Sharma P, Hu-Lieskovan S, Wargo J, Ribas AJC. Primary, adaptive, and acquired resistance to cancer immunotherapy. 2017;168, 707-723.

  42. 42.

    Olivo Pimentel V, Marcus D, van der Wiel A, Lieuwes N, Biemans R, Lieverse R, et al. Releasing the brakes of tumor immunity with anti-PD-L1 and pushing its accelerator with L19-IL2 cures poorly immunogenic tumors when combined with radiotherapy. 2021;9.

  43. 43.

    Hensler M, Kasikova L, Fiser K, Rakova J, Skapa P, Laco J, et al. M2-like macrophages dictate clinically relevant immunosuppression in metastatic ovarian cancer. 2020;8.

  44. 44.

    An Y, Yang QJLs. MiR-21 modulates the polarization of macrophages and increases the effects of M2 macrophages on promoting the chemoresistance of ovarian cancer. 2020;242, 117162.

  45. 45.

    Guo L, Zhang Y, Yang Z, Peng H, Wei R, Wang C, et al. Tunneling nanotubular expressways for ultrafast and accurate M1 macrophage delivery of anticancer drugs to metastatic ovarian carcinoma. 2019;13, 1078-1096.

  46. 46.

    Ray M, Ruffalo M, Bar-Joseph ZJG, chromosomes, cancer. Construction of integrated microRNA and mRNA immune cell signatures to predict survival of patients with breast and ovarian cancer. 2019;58, 34-42.

  47. 47.

    Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16:284–7.

    CAS  Article  Google Scholar 

  48. 48.

    Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nature communications. 2013;4:2612.

    Article  Google Scholar 

  49. 49.

    Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nature methods. 2015;12:453–7.

    CAS  Article  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This work was funded by the National Natural Science Foundation of China (Nos. 82002741, 81960465, 82060471). No funding or other benefits related to the subject of this article were received from any commercial entity.

Author information

Affiliations

Authors

Contributions

WJ designed the study; RXC, XDZ, RTF, and LT prepared the figures; PL contributed to drafting the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Wei Jia.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have 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: Figure S1:

Proportions and relevance of LM22 immune cells in other databases: (A) GSE32062-HGSOC, (B) GSE53963-HGSOC, (C) TCGA-OC, (D) Serous, (E) Endometrioid, and (F) Clear cell database. HGSOC, High-grade serous ovarian carcinoma; LM22, leukocyte gene signature matrix; TCGA, The Cancer General Atlas

Additional file 2: Figure S2:

Correlation between PD-1/PD-L1 expression and three immune cells infiltrating levels: M1 macrophages, M2 macrophages, and resting dendritic cells in (A) GSE53963-HGSOC, (B) GSE32062-HGSOC. HGSOC, High-grade serous ovarian carcinoma; PD-1, programmed death-1; PD-L1, programmed death ligand 1

Additional file 3: Table S1:

Clinical data for sample id, overall survival time, survival state, and the average expression of PD-1/PD-L1 in GSE53963, GSE32062, GSE135820, and TCGA-OV databases. PD-1, programmed death-1; PD-L1, programmed death ligand 1; TCGA, The Cancer General Atlas

Additional file 4: Table S2:

Proportions of infiltrating immune cells of each sample in GSE53963, GSE32062, GSE135820, GSE17260, GSE73614, GSE68600, GSE73638, TCGA-OC, and GETx databases. GTEx, Genotype-Tissue Expression; OC, Ovarian cancer; TCGA, The Cancer General Atlas

Additional file 5: Table S3:

Average proportions of 22 immune cells subsets in GSE135820-HGSOC, GETx-Normal, Serous, Endometrioid, Clear cell, and TCGA-OC databases. GTEx, Genotype-Tissue Expression; HGSOC, High-grade serous ovarian carcinoma; OC, Ovarian cancer; TCGA, The Cancer General Atlas

Additional file 6: Table S4:

Immune cells could be classified into three optimal cluster groups in 2873 HGSOC samples with CIBERSORT: cluster 1/A (n = 1961), cluster 2/B (n = 737), and cluster 3/C (n = 175). HGSOC, High-grade serous ovarian carcinoma

Additional file 7: Table S5:

Correlation between overall survival of 239 HGSOC patients and the content of APC infiltrating cells (M1 macrophages, dendritic cells, and B cells).

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

Verify currency and authenticity via CrossMark

Cite this article

Liu, P., Chen, R., Zhang, X. et al. Combined PD-1/PD-L1 and tumor-infiltrating immune cells redefined a unique molecular subtype of high-grade serous ovarian carcinoma. BMC Genomics 23, 51 (2022). https://doi.org/10.1186/s12864-021-08265-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-08265-y

Keywords

  • Antigen-presenting cells
  • High-grade serous ovarian carcinoma
  • Immune cell infiltration
  • PD-1
  • PD-L1