The Prognostic Effect of Immune Cell Infiltration Depends on Molecular Subtype in Endometrioid Ovarian Carcinomas.

Historically, ovarian carcinomas have been treated uniformly regardless of biological features of each histologic type, a trend that remains overwhelmingly common. In the last decade, multiple studies have shown that a standardized and molecular marker–based approach to determine the histologic subtypes of ovarian carcinoma can be helpful to correctly stratify patients and will thus result in improved reproducibility across biomarker and outcome studies of ovarian carcinoma ( 1–5 ). Knowledge about the importance of the immune tumor microenvironment of ovarian carcinoma has been growing in the recent years; however, studies directed to ovarian carcinomas have thus far primarily focused on high-grade serous ovarian carcinoma (HGSOC). This is justifiable given it is the most common type and responsible for the greatest mortality. Unfortunately, this has compounded poor understanding of less common histologies. Endometrioid ovarian carcinoma (ENOC) remains one of the least understood histologic subtypes. Recently, we contributed to studies involving ENOC which showed that there is room for forecasting outcome based on the presence of immune cells. There, mild but significant overall survival benefits were observed in ENOC with higher levels of intraepithelial CD8 + tumor-infiltrating lymphocytes (TIL) that also correlated with ARID1A loss ( 6, 7 ). Our group and others have taken strides in understanding the biology of ENOC, leveraging knowledge of endometrial carcinoma ( 8–11 ). ENOC accounts for approximately 10%–20% of all ovarian carcinoma cases ( 12–14 ) and are one of two endometriosis-associated ovarian carcinomas, the other being clear cell ovarian carcinoma (CCOC). ENOC arises from endometrial-type epithelium within endometriotic lesions that undergo malignant transformation ( 15, 16 ). With a shared origin of endometrial cells, ENOC tumors are histologically and molecularly similar to endometrial carcinoma, and especially endometrioid endometrial carcinomas (EEC; refs. 8–10 ). Both endometrial carcinoma and ENOC can be stratified into near-equivalent molecular subtypes: (i) an ultra-mutated and DNA polymerase epsilon (POLE)-mutant subgroup that has major defects in DNA proofreading repair ( POLE mut; ref. 17 ); (ii) a hypermutated subgroup with DNA mismatch repair deficiencies and a microsatellite instable phenotype (MMRd); (iii) a high-copy number and p53-mutant subtype (p53abn); and (iv) a low-copy number subtype, lacking the above noted features, known as no specific molecular profile (NSMP)—a potentially erroneous moniker given it references only a lack of the other three profiles. Subtypes have been identified through genomics or a small number of surrogate IHC and sequencing assays ( 8, 18–25 ). Recent successes of immunotherapy in immunologic “hot” tumor types have fueled efforts to identify ovarian carcinoma subgroups that may benefit from this new form of cancer treatment ( 26–28 ). Across many solid tumors, the density of TILs appears linked to tumor mutational burden (TMB) and the presumed accumulation of neoantigens ( 27 ). In several cancers, increasing TILs and TMB have been reported to correlate improved clinical response to immune checkpoint inhibitors (ICI; refs. 29, 30 ). Even in the absence of ICI-based therapy, multiple studies have shown that high rates of TILs correlate with improved outcomes in HGSOC and, to a lesser extent, ENOC ( 6, 7, 31, 32 ). Within both ovarian carcinoma and ENOC, both POLE mut and MMRd molecular subtypes are associated with the highest TMB ( 18, 25, 33, 34 ) and, for endometrial cancer, there is considerable evidence to support robust immune response within these ( POLE mut and MMRd) molecular subtypes ( 7, 11, 35–37 ). Our group recently reported the equivalency of molecular profiles present in both endometrial cancer and ENOC, wherein the genomic profile or surrogate assays can provide substantial prognostic information about four molecular subtypes ( 8, 18, 20, 21, 25 ). Similarities between ENOC and endometrial cancer suggest that pursuing medical management strategies for ENOC that are parallel to endometrial cancer/endometrioid endometrial cancer may be beneficial ( 38–41 ). Following this work, we now investigate the immune tumor microenvironment (iTME) in the context of molecular subtypes, to evaluate the potential additive prognostic value on overall survival (OS), disease-specific survival (DSS), and progression-free survival (PFS) in ENOC. This detailed understanding of immune infiltration in ENOC may highlight subsets for which ICIs are of potential benefit versus those unlikely to benefit. We applied multiplex immunofluorescent/IHC staining to tissue microarrays representing 210 ENOC and examined the densities of different immune cell populations across tumor epithelium and associated stroma, their relationships with molecular subtype, and resulting outcome. While molecular subtypes in ENOC and endometrial cancer are near equivalent ( 42, 43 ), each disease does have unique features and we evaluated whether iTME could provide disease or subtype-specific information to ENOC that was distinct from reports in endometrial cancer ( 23 ).

Our cohort consisted of a subset of previously described studies of ENOC ( n = 262; refs. 6–8 ). Tissue samples from the following Canadian and German centers were used: Department of Women's Health, Tübingen University Hospital (Tübingen, Germany); Department of Gynecology and Gynecologic Oncology, Evang. Kliniken Essen-Mitte (Essen, Germany); Medizin Campus (Bodensee, Friedrichshafen, Germany); the OVCARE gynecological tissue bank (Vancouver, Canada). Research was conducted in accordance with the Canadian Tri-Council Policy Statement on Ethical Conduct for Research Involving Humans (TCPS2, 2018) which is fully compliant with the Declaration of Helsinki. In cases where patients were noncontactable, and treated more than 5 years before the start of the study, they were included under institutionally approved waiver of consent (Tübingen University Hospital Research Ethics Board). All other patients provided informed consent. Collection and use of clinical samples, and respective clinical data, was approved through local research ethics boards. For all Vancouver cases, the date of (debulking) surgery was set as the time of origin and the follow-up period consequentially started from this date. For the German cases, the first date of histologic proof was set as start point which mostly coincided with the debulking surgery. In rare circumstances, a diagnostic surgery preceded the main debulking surgery which were only a few days apart. We used a “cohort wise” follow-up approach censoring observations on December 31st of each year following the year in which they had their surgery (year of diagnosis) to ensure random censoring and minimize ascertainment bias. Definitions of endpoints (OS, DSS, PFS) were as follows: OS: the time from surgery until the death of any cause. Patients, still alive prior to the censoring date, are censored at the time of last follow-up. DSS: the time from surgery until death due to ENOC. Patients, alive prior to the censoring date, or who died of an unrelated cause, are censored at the time of last follow-up. Unrelated causes of death were defined by the attending physician and/or rigorous chart review (reported cause of death) during clinical data abstraction. PFS: the time from surgery until there was evidence of recurrent or progressive disease (clinical evidence of recurrence, or imaging confirmation). Patients, alive and disease-free prior to the censoring date, or who died of an unrelated cause, are censored at the time of last follow-up. Unrelated causes of death were defined by the attending physician and/or rigorous chart review (reported cause of death) during clinical data abstraction. OS: the time from surgery until the death of any cause. Patients, still alive prior to the censoring date, are censored at the time of last follow-up. DSS: the time from surgery until death due to ENOC. Patients, alive prior to the censoring date, or who died of an unrelated cause, are censored at the time of last follow-up. Unrelated causes of death were defined by the attending physician and/or rigorous chart review (reported cause of death) during clinical data abstraction. PFS: the time from surgery until there was evidence of recurrent or progressive disease (clinical evidence of recurrence, or imaging confirmation). Patients, alive and disease-free prior to the censoring date, or who died of an unrelated cause, are censored at the time of last follow-up. Unrelated causes of death were defined by the attending physician and/or rigorous chart review (reported cause of death) during clinical data abstraction. ENOC molecular subtype labels and ARID1A status were taken from previous reports on this cohort (see also Fig. 1 ; refs. 7, 8 ). β-Catenin, progesterone, and estrogen receptor status was assessed according to established IHC protocols wherein a portion of the Canadian cohort has been previously reported on ( 44–46 ). Cohort classification. A, Schematic showing the samples included for the immune microenvironment analysis based on available data. B, Molecular subtype assignment of ENOC cases based on surrogate biomarkers and complete mcIHC panel data. Multicolor IHC (mcIHC) and multispectral imaging was used to identify the iTME in epithelial and stromal compartments ( 23 ). All experiments were done on tissue microarray (TMA) platforms with samples having either 2 × 0.6 mm cores, or 3 × 0.6 mm cores. The following four costaining antibody panels were used to identify T-cell lineages (panel I) CD3/CD8/PanCK, (panel II) CD25/FoxP3/CD8/PanCK; B-cell lineages (panel III) CD20/CD79a/PanCK; and macrophages (panel IV) PD-1/PD-L1/CD68/PanCK (see below for details). PanCK staining in the above-noted panels was typically used to identify epithelial/tumor versus stromal compartments ( Fig. 2 ); however, cell-specific costaining with PanCK was also used to identify populations within the panel II. For each tumor core the automated tissue segmentation and cell signal data were used to distinguish tumor and stromal compartments as well as to calculate density values for various (co-)stained cell populations. Cell phenotypes were defined by presence vs. absence of a given marker with the exception of CD25 where low and negative expression was binned (denoted as CD25 L− ). For each marker panel 9%–15% of the automated data was randomly selected for validation by manual counting. To ensure representation of both epithelial and stromal cellular compartments, TMA cores were excluded when composed of < 25% (by core area) of each epithelial and stromal tissue, or when the majority of the evaluated tissue (≥95%) was recognized as stromal. When replicate cores were interpretable the mean density value was used. Multicolor IHC. A–L, Representative images of multicolor IHC and multispectral imaging of TILs for each immune panel, including images of applied automated tissue segmentation into tissue compartments. IHC staining for the different immune marker was performed according to protocols published in Talhouk and colleagues ( 23 ). For the B-cell lineage, IHC clone SP18 (rabbit monoclonal, Abcam) for CD79a, and clone L26 (mouse monoclonal, Biocare) for CD20 were used. For the IHC T-cell lineage, clone SP7 (rabbit monoclonal, Abcam) for CD3, clone C8/144b (mouse monoclonal, Cell Marque) for CD8, clone 4C9 (mouse monoclonal, LabVision) for CD25, and clone 236A/E7 (mouse monoclonal, Abcam) for FoxP3 were applied. For the macrophage panel, IHC was performed with clone EPR4877 (rabbit monoclonal, Abcam; ref. 2 ) for PD-1, clone SP142 (rabbit monoclonal, Abcam) for PD-L1, and clone SP251 (rabbit monoclonal, Abcam) for CD68. PanCK staining was generated using clone AE1/AE3+5D3 (mouse monoclonal, Biocare). For signal detection, chromogens DAB, Warp Red, and Ferangi Blue and different fluorophores per panel (Treg:CD25–520 nm, CD8–570 nm, FoxP3–690 nm, PanCK-Coumarin; TAM: PD-1–650 nm, CD68–520, PD-L1–570 nm, PanCK-690 nm) were used. Density of immune cell infiltrates was assessed using automated image analysis and scoring software, followed by visual confirmation. In detail, stained TMA slides were scanned with Vectra Multispectral Imaging system (PerkinElmer) and InForm 2.1 software (PerkinElmer) was utilized to quantify the immune cell populations based on a minimum of three different algorithms. InForm algorithms were initially trained on ENOC and CCOC images for spatial recognition of epithelial and stromal compartments. If neither core was assessable, data were noted as “missing” for that marker. In rare circumstances when data from only one core was assessable, those were utilized. If cases were present on multiple TMAs, average values per TMA were recorded before calculating the overall case mean. Because average density values of several markers were zero-scored, a log (base 2) transformation was performed to received delta log normalized data (Supplementary Fig. S1A and S1B; refs. 23, 47 ). Similar to published approaches ( 23 ), the immune signatures were subject to unsupervised clustering using the R package diceR ( 48 ). The final consensus clustering based on an aggregating ensemble approach was established using the k-modes method. Survival benefits in association with molecular subtypes and different immune cluster were assessed using Cox proportional hazard models and Kaplan–Meier curves with the log-rank test. R packages survival and survminer were used. Multivariable analysis was performed, correcting for the known outcome modifier: age, stage, grade, residual disease, and treatment as well as ARID1A status and the immune clusters. Statistical significance was defined by P < 0.05. Anonymized data generated in this study are available upon request from the corresponding author. Some data may be restricted or subject to material/data transfer agreement to ensure compliance with patient consent.

Of 338 ENOC cases evaluated for this study, 128 were excluded based on quality assurance or missing data. A total of 210 cases were retained with full immune profiling data and, of those, 190 had available data on molecular subtype and other biomarkers ( Fig. 1A ). Clinicopathologic cohort characteristics were consistent with our prior report ( Table 1 ; ref. 8 ). The median age at diagnosis was 55 years. A total of 54.76% of cases were FIGO stage I, and 50.95% were grade 1. As expected, the molecular subtype distribution was 23/190 (12%) MMRd, 8/190 (4%) POLE mut, 17/190 (9%) p53abn, and 142/190 (75%) NSMP ( Fig. 1B ; Supplementary Table S1A). Loss of ARID1A immunoreactivity—a surrogate for loss-of-function mutations ( 49 )—affected 52 cases (25%), including 46 with complete loss (22%) and 6 with subclonal loss (Supplementary Table S1A). As expected, ARID1A loss was most prevalent in the MMRd subtype, with 14/23 (61%), but was also observed in other subtypes with 2/8 (25%) in POLE mut, 29/139 (21%) in NSMP, and 1/17 (5.9%) in the p53abn subtype. Outcomes based on clinicopathologic parameters (age, stage, grade, residual disease) and the predictive value of the molecular subtypes were not different from our previous reports (Supplementary Table S1B; Supplementary Fig. S2), and similar to endometrial carcinoma ( 8, 20 ). Significant associations were seen between molecular subtype and clinicopathologic variables of grade and stage ( P = 0.018) as well as with ARID1A status ( P < 0.001; Supplementary Table S1C). Cohort characteristics. Abbreviation: IQR, interquartile range. By analyzing stromal and epithelial compartments of each tumor separately, it was apparent that immune populations were not evenly spread among the two compartments or across molecular subtypes (Supplementary Table S2). Overall, immune infiltrates were less abundant in the T-lineage panel II in the epithelium than in the macrophage panel III. The detection of certain populations in different compartments correlated with outcomes. Stromal CD25 L− /FoxP3 − /CD8 + T cells mildly but significantly corresponded with the outcome of all survival parameters (OS: HR, 0.929; P < 0.001; DSS: 0.917; P < 0.001; PFS: 0.933; P = 0.012; Supplementary Table S3), as did stromal plasma cells (CD20 + /CD79a + ) for OS (HR, 0.964; P = 0.034) and DSS (HR, 0.952; P = 0.028). In the epithelial compartment, the density level of PD-1 + cells (HR, 0.954; P = 0.012) and FoxP3 + tumor cells (CD25 L− /FoxP3 + /CD8 − /PanCK + ) coincided with better OS (HR, 0.929; P = 0.017). In our overall cohort, we applied unsupervised clustering to the IHC dataset. The optimal number of immune clusters was determined using an unsupervised approach through the R package diceR ( 48 ). To facilitate comparison to prior reports in endometrial cancer we applied a k = 2 solution, which was determined by diceR as the optimal cluster solution for our ENOC cohort (Supplementary Fig. S1C and S1D), and the next most optimal cluster solution k = 3. These divided our cohort roughly into TIL high and TIL low (k = 2), or TIL high , TIL medium , and TIL low (k = 3) clusters, respectively. Underlying or driving components of each cluster are discussed below. The compartment-independent immune signatures/phenotypes were plotted in a grouped heat map based on cluster and molecular subtype assignments ( Fig. 3A and B ). For most markers, the stromal and epithelial values clustered in close proximity. The TIL high cluster was characterized by higher levels across most TIL subsets, in particular regulatory (CD25 + /FoxP3 + /CD8 − ) and cytotoxic T cells (CD25 L− /FoxP3 − /CD8 + and CD3 + /CD8 + ). In contrast, the TIL low cluster had lower levels of most TIL subsets and negligible B lineage (populations marked CD20 and/or CD79a) and PD-L1 + cell populations. Nonetheless, observable differences did not translate into distinct survival outcomes with k = 2 clustering ( P > 0.1; Fig. 3C ). Cluster analysis of TIL patterns for the optimal cluster solution k = 2 ( A ) and cluster alternative k = 3 ( B ). TIL high tumors showed significantly higher densities in both T-cell lineage TIL subsets. C, Forest plot showing the individual HRs of both cluster solutions in regards to OS, DSS, and PFS. The three-cluster solution presented more well-defined signatures and seemed to adhere to functional groupings; therefore, further analyses were carried out based on this. Within this, TIL high was dominated by regulatory (CD25 + /FoxP3 + ) and cytotoxic T (CD3 + /CD8 + and CD25 L− /FoxP3 − /CD8 + ) cells. The TIL medium cluster also had high levels of T cells (CD3 + populations) and macrophages (CD68 + populations) and modest levels of intraepithelial CD68 − /PD-L1 + cells but negligible levels of B-lineage cells (CD20 + and/or CD79a + ). Consequently, we deemed this the TIL B minus cluster. Finally, the TIL low cluster had a paucity of T cells, B cells, and PD-1 + and PD-L1 + populations and was notably devoid of macrophages (CD68 + ; Fig. 3B ). Intriguingly, the TIL low cluster was significantly enriched for FoxP3 + tumor cells (CD25 L− /FoxP3 + /CD8 − /PanCK + ) compared with TIL high /TIL B minus combined ( P < 0.001), whereas no such differences were found between TIL high and TIL B minus tumors ( P = 0.11). However, the k = 3 immune clusters also did not show distinct consistent survival outcomes for ENOC ( P > 0.1). ENOC typically presents at lower stage, thus we stratified analysis for low (FIGO I/II) and high (FIGO III/IV) stage. Immune clusters showed no prognostic effect in low-stage ENOC ( P > 0.1; Table 2A ). In high-stage ENOC, immune cluster groups showed a significant difference in OS (HR B minus , 4.91; 95% CI, 1.58–15.3; P < 0.05), and a trend in PFS (HR B minus , 3.57; 95% CI, 1.13–11.3; P = 0.066; Table 2A ). In multivariable analysis, adjusting for clinicopathologic and molecular parameters (age, residual disease, grade, treatment, molecular subtype, and ARID1A status), high-stage TIL B minus cases did worse than TIL high and TIL low (HR high , 0.19; 95% CI, 0.02–2.12; HR B minus , 3.23; 95% CI, 0.52–23.0—reference TIL low , respectively, both P < 0.05; data not shown). Summary of P values obtained during Kaplan–Meier survival analysis in regard to TILs and subtypes. Note: A, Effects of the immune cluster and molecular subtype on survival outcomes within ENOC. B and C, Survival analysis in regard to TIL clusters within the respective molecular subtype in all of ENOC ( B ) and only low-staged ENOC ( C ). Because of a lack of cases and events, analyses within the POLE mut subtype were inconclusive and therefore excluded. We saw significant differences in immune cell density between subtypes in the epithelial compartment for cytotoxic T-cell populations (CD3 + /CD8 + P = 0.043; and CD25 L− /FoxP3 − /CD8 + P < 0.001; Supplementary Table S2), FoxP3 + cells (CD25 L− /FoxP3 + /CD8 − P = 0.013), PD-1 + cells ( P = 0.005), and CD68 − /PD-L1 + cells ( P = 0.011). While these differences were generally driven by lower fraction of observed expression (52%–88% of detected expression) and lower infiltration (density values) in the p53abn subtype, CD68 − /PD-L1 + cells also had low density in the NSMP subtype. The CD25 L− /FoxP3 − /CD8 + population was the only one with subtype differences in both epithelial ( P < 0.001) and stromal compartments ( P = 0.043), dominantly driven by low density in p53abn (log-transformed values of −2.7 and 2.2, respectively). Stroma only infiltration was different across subtypes for CD3 + /CD8 − cells ( P = 0.048) with the highest densities in the NSMP group (Supplementary Table S2). CD68 + /PD-L1 − macrophage densities also varied by subtype in the stromal compartment ( P = 0.004) driven by high density in POLE mut/MMRd and marked lower density in p53abn. B-lymphocyte lineage populations including cells staining for CD20 and/or CD79a were the least prevalent populations with overall low densities observed across all subtypes. Although higher plasma cell density was seen in POLE mut/MMRd, and exceptionally low densities in p53abn, this was not significant (epithelial compartment P = 0.099; Supplementary Table S2). Because our previous EC study showed that molecular subtypes were superior to TIL patterns for prognostication ( 23 ), we reanalyzed ENOC cases including molecular subtype and immune cluster labels. In univariable analysis, immune clustering had no benefit, while prognostic value of the molecular subtype retained its significance in DSS and PFS ( P < 0.05; Supplementary Table S4). We next stratified the analysis by subtype. POLE mut ENOC cases ( n = 8) were omitted from this analysis due to small numbers and a lack of death or progression events. In the NSMP subtype, the TIL B minus cluster performed significantly worse for OS and DSS than the TIL high and TIL low clusters (OS: HR, 3.64; 95% CI, 1.31–10.1; DSS: HR, 4.03; 95% CI, 1.11–14.7; both P < 0.05; Table 2B ; Fig. 4 ). Benefit of immune response in NSMP ENOC. Kaplan–Meier curves of OS, DSS, and PFS in regards to TIL pattern k = 3 in all NSMP ( A – C ), in low-stage NSMP ( D – F ), and high-stage NSMP ( G – I ). We repeated this analysis further stratifying cases by low and high stage. In the MMRd and p53abn subtype, no differences were seen between TIL clusters in low-stage cases; analyses were not performed on high-stage cases due to insufficient sample size numbers ( Table 2C ). Amongst high-stage ENOC the more prevalent NSMP group showed striking differences for the k = 3 solution in OS and DSS (HR > 20; 95% CI, 0.0–0.0; both P 1.4; P ≤ 0.01), but not for high-stage restricted DSS analysis (Supplementary Table S5). Within the NSMP ENOC, the TIL B minus cluster was characterized by an absence of plasma and B cells (CD20 −/+ /CD79a + ), which was also seen with the TIL low cluster ( P > 0.3, data not shown). A retention of macrophages (CD68 + ), especially PD-L1 − macrophages, and cytotoxic T cells (CD3 + /CD8 + and CD25 − /FoxP3 − /CD8 + ) similar to the TIL high cluster was seen. Furthermore, intermediate to high levels of PD-1 + cells and regulatory T cells (CD25 + /FoxP3 + /CD8 − and CD25 − /FoxP3 + /CD8 − respectively) were observed in this NSMP TIL B minus group. Surprisingly, the highest mean density of CD3 + /CD8 − (helper) T cells was recorded in this subcohort (log transformed mean: 0.35 in TIL B minus vs. −0.67 in TIL high and −4.6 in TIL low ). This was also observed in high-stage NSMP TIL B minus cases, where additionally not even one case with infiltrating B cell was recorded. Investigating correlations with other ENOC prognostic biomarkers, we did not see associations between the TIL B minus cluster and any of progesterone receptor (PR), estrogen receptor (ER), or β-catenin positivity (all P > 0.1; Supplementary Fig. S3). However, we did note that a larger proportion of NSMP TIL B minus cases had PR-positive staining compared with the other two clusters (88.4% vs. 83.7/74% in TIL high/low ; P = 0.188; Supplementary Fig. S3). With data generally similar to trends observed in endometrial cancer we reexamined EEC from our previous dataset ( 23 ) to compare and contrast with ENOC. The EEC subtype distribution and resulting outcomes were generally similar to the full endometrial cancer cohort (Supplementary Fig. S4B–S4D; ref. 23 ). The majority of EEC cases (55.6%) were NSMP, 25.8% MMRd, 10.8% POLE mut, and 7.8% p53abn. A full side-by-side comparison of the considered immune cell populations could not be done since partially different IHC/IF panels were used between the two data sets (e.g., no assessment of FoxP3 and PanCK co-staining cells in EEC). Like ENOC, immune clustering in either k = 2 or k = 3 did not result in significant survival benefit among the full EEC cohort (OS, DSS, PFS P > 0.5; Supplementary Fig. S4). In contrast to ENOC, in the NSMP EEC subcohort we saw no survival advantage with the k = 3 immune cluster solution ( P > 0.05; Supplementary Fig. S4E–S4G).

In this study, we evaluated the iTME across the four molecular subtypes of endometrioid ovarian carcinoma. Immune cell infiltrates showed three distinct patterns: a “TIL low ” or cold cluster, a “TIL high ” or hot cluster with a combined T- and B-cell response, and a “TIL B minus ” cluster that was distinguished by T cells with negligible B and plasma cells. TIL levels were higher in high-stage ENOC. Across the overall ENOC cohort, molecular subtyping presented superior prognostic significance compared to immune cell infiltrates alone. Nonetheless, after molecular subtype stratification, immune infiltrates showed prognostic value in the NSMP subtype, with the TIL B minus cluster showing significantly worse survival than the TIL high and TIL low clusters. We examined a cohort of 210 ENOC which were previously characterized by pathology review and IHC-based confirmation of the histotype ( 8 ). A total of 190 cases had full molecular subtype annotation; the prevalence of these subtypes was similar to other reports, as were the overall clinicopathological features (age, stage, grade, residual disease) and ARID1A loss enriched in the MMRd subtype ( 7, 8, 11 ). On average, when evaluating the different compartment restricted immune cell populations, an association with outcome benefits was seen across the full ENOC cohort. High level of cytotoxic CD8 + (CD25 L− /FoxP3 − /CD8 + ) T cells or plasma cells (CD20 + /CD79a + ) in the stroma were linked to better survival outcomes, consistent with HGSOC ( 50, 51 ). However, given a relatively weak HR for this observation, including all ENOC, it is likely a diluted effect from the majority NSMP subset. It should also be noted our cohort was smaller than the previous reports and the described effect was subtle in ENOC ( 6, 7 ). Clustering of the immune cell densities revealed three dominant immune phenotypes (TIL high , TIL B minus , and TIL low ); however, examination of ENOC as a whole suggested these phenotypes had no significant outcome benefits. High-stage ENOC outcomes appeared to be more influenced by immune infiltration and clusters with generally higher levels of TILs had better OS with supporting trends in DSS and/or PFS, although this may be confounded with molecular subtype ( 7, 8 ). While a general trend of high-moderate-low overall immune infiltration could be observed across these clusters, FoxP3 + tumor cells showed a distinct pattern. FoxP3 is typically used as a marker of regulatory T cells; however, our and others’ analyses indicate FoxP3 can also be expressed by tumor cells ( 52, 53 ). Although it has been suggested that epithelial FoxP3 expression influences carcinogenesis and outcome in other cancers ( 54, 55 ), no effect on outcomes was observed here in ENOC. Further studies into the nature of FoxP3 + tumor cells may still be warranted to better describe this biologically distinct group with otherwise low levels of immune infiltration. Molecular subtype influenced the immune cell density across both the epithelial and stromal compartments. Consistent with findings in endometrial cancer ( 23 ), molecular subtypes with higher mutational burden tended to have a greater level of immune cell infiltration. That said, a wide range of immune cell densities was visible within each subtype, and no immune cluster group was entirely confined to any one molecular subtype. POLE mut and MMRd subgroups were enriched for TIL high tumors and had the highest levels of intraepithelial PD-L1 + cells. Consistent with this, PD-L1 positivity has been associated with improved outcomes in HGSOC and endometrial cancer ( 56 , 57 ). High infiltration of PD-1 + cells was also observed across epithelial compartments of all molecular subtypes, except for p53abn cases. In p53abn ENOC an absence and low density of cytotoxic T (CD8 + ) cells was recorded in stromal and epithelial compartments. In fact, only a few CD3 + /CD8 − lymphocytes were detected in p53abn stroma, while they were most prevalent in the stromal compartment of NSMP tumors (p53 wild-type). This may be consistent with recent reports suggesting that perturbed p53 leads to an immune suppressive and tumor promoting environment, manifested as reduced cytotoxic cells and increased inhibitory myeloid cells ( 58, 59 ). Overlaying molecular subtype and immune density clustering did not add additional prognostic value over molecular subtype alone. However, restricted analysis of NSMP ENOC suggested the TIL B minus cluster performed significantly worse than others. Likewise, superior performance of TIL high tumors was most evident in high-stage NSMP, however, this observation was based on a small number of cases ( n = 20) and has to be seen in more of an exploratory sense. The TIL B minus NSMP cluster was characterized by relatively high levels of PD-L1 − macrophages, cytotoxic T cells, CD8 − (helper) T cells and FoxP3-expressing (regulatory) T cells in epithelial compartments. Indeed, the levels of those cell types in TIL B minus were not significantly different from the TIL high cluster; however, it remains unknown whether these immune infiltrates might be exhausted or dysfunctional within their iTME. Likewise, TIL B minus tumors almost entirely lacked B or plasma cells. In other cancers, high levels of B cells, especially in tertiary lymphoid structures, can modulate the overall iTME and positively influence outcome and prognosis ( 51, 60 ). One could additionally speculate that the lack of B cells in TIL B minus tumors shifts T cells towards a more suppressive phenotype, which aligns with accumulating evidence that B-lineage cells have a net positive influence on antitumor immunity ( 61 ). Because our IHC/IF data were based on TMAs rather than whole tissue sections, it is possible that we have underestimated B and plasma cell responses involving tertiary lymphoid structures. Nonetheless, high levels of cytotoxic and Th cells PD-L1 + macrophages, suggest PD-1/PD-L1-directed ICI therapy may be effective ( 62 ). Similarly, cases with a high level of regulatory T cells may also benefit from ICIs as activating these populations has been shown to recruit cancer antigen–specific effector cells and potentiate immune therapy responses ( 63 ). These observations should be correlated to data from ICI trials including ENOC and/or in functional, immune-competent, models. This is among the first reports of TIL low (ovarian) tumors being associated with more favorable prognosis than TIL B minus tumors, which exhibit a medium level of immune infiltration. This could be attributed to a low-stage cohort, typical of ENOC. Indeed, while immune cell cluster distributions were similar in stage I and II cases, an enrichment of TIL high tumors was seen stage III and IV cases. However, we do want to point out that an effect overestimation based on small sample size in the sub analyses may be present. Given the known similarity between endometrial cancer, and especially EEC, with ENOC we were unsurprised to find similar trends related to immune cell infiltration ( 23 ). However, as there were differences in the immune cell panels examined between endometrial cancer and ENOC cohorts some subtlety between the uterine versus ovarian environments may have been missed. In fact, we were not able to recapitulate a poor-outcome TIL B minus cluster in EEC and it is unclear if this is a result of ovarian-specific microenvironment or panel differences between studies. Nonetheless our data are still supportive of future ENOC trials being combined with endometrial cancer /EEC. This may be particularly beneficial for NSMP tumors which are not well characterized in either the ovarian or uterine presentation. While NSMP tumors are proportionally less common among endometrial cancer/EEC compared with ENOC, endometrial cancer overall is a much more common cancer. Large-scale studies of endometrial cancer suggest there are a number of mutations prevalent across NSMP tumors, however, those are less uniform than the other subtypes such that NSMP is generally regarded as a subtype by omission: lacking POLE mutation, MMRd/MSI, and high levels of DNA copy-number alterations (p53-abnormality). The same is true for their ENOC counterparts, and both may benefit from further stratification. The value of biomarkers that add prognostic (and potentially predictive) information to NSMP tumors is high. Our study suggests immune infiltration may be a valid and independent stratification mechanism in NSMP. However, as the usage of multiplex IHC and/or IF is limited in current clinical practice, simpler signatures and surrogate biomarker investigations are still needed to make these findings clinically applicable. In ENOC and endometrial cancer, the presence of hormone receptors, especially PR, has been associated with favorable outcome ( 45, 64 ). Indeed investigation of hormone-based therapies is a prime candidate for NSMP endometrial cancer trials ( 65 ). In our NSMP ENOC cohorts, both PR and ER showed prognostic significance. However, no correlation was seen between the immune clusters and positivity for either hormone receptor. This suggests that the prognostic effect of the immune cluster in NSMP ENOC may be independent of PR/ER status. Our study is the first multifaceted iTME evaluation in ENOC in the context of molecular subtypes. This highlights the value of considering molecular subtype stratification in prognostication prior to targeted therapeutic considerations or other stratification. Consequently, for prognostication of NSMP cases, we suggest a tiered stratification system first by molecular subtype and second by immune response may be beneficial. Our data also propose that the NSMP subtype warrants consideration for immunotherapy clinical trials.

Supplementary Figures S1-S4 Supplementary Tables S1-S5