<figure><img src="image_1.png" /></figure>

Check for
updates

ARTICLE

Molecular Diagnostics

# Combining measures of immune infiltration shows additive effect on survival prediction in high-grade serous ovarian carcinoma

Anne Montfort<fnref n="1" />, Rowan J. Barker-Clarke²,³, Anna M. Piskorz², Anna Supernat⁴, Luiza Moore²,³, Sarwah Al-Khalidi², Steffen Böhm<fnref n="1" />, Paul Pharoah²,⁵, Jacqueline McDermott<fnref n="1" />,⁶, Frances R. Balkwill<fnref n="1" /> and James D. Brenton ID²,³,⁷

**BACKGROUND:** In colorectal and breast cancer, the density and localisation of immune infiltrates provides strong prognostic information. We asked whether similar automated quantitation and combined analysis of immune infiltrates could refine prognostic information in high-grade serous ovarian carcinoma (HGSOC) and tested associations between patterns of immune response and genomic driver alterations.

**METHODS:** Epithelium and stroma were semi-automatically segmented and the infiltration of CD45RO⁺, CD8⁺ and CD68⁺ cells was automatically quantified from images of 332 HGSOC patient tissue microarray cores.

**RESULTS:** Epithelial CD8 [p = 0.027, hazard ratio (HR) = 0.83], stromal CD68 (p = 3 × 10⁻⁴, HR = 0.44) and stromal CD45RO (p = 7 × 10⁻⁴, HR = 0.76) were positively associated with survival and remained so when averaged across the tumour and stromal compartments. Using principal component analysis, we identified optimised multiparameter survival models combining information from all immune markers (p = 0.016, HR = 0.88). There was no significant association between PTEN expression, type of TP53 mutation or presence of BRCA1/BRCA2 mutations and immune infiltrate densities or principal components.

**CONCLUSIONS:** Combining measures of immune infiltration provided improved survival modelling and evidence for the multiple effects of different immune factors on survival. The presence of stromal CD68⁺ and CD45RO⁺ populations was associated with survival, underscoring the benefits evaluating stromal immune populations may bring for prognostic immunoscores in HGSOC.

British Journal of Cancer (2020) 122:1803–1810; https://doi.org/10.1038/s41416-020-0822-x

INTRODUCTION

There is a heterogeneous immune response in the tumour microenvironment of HGSOC, but the presence of intra-epithelial CD8⁺ T cells is consistently associated with improved survival.<fnref n="1" />,² Prolonged survival is also associated with the presence of B cells and other immune cells, including CD45RO⁺ memory cells.³⁻⁶

Prognostic scoring has predominantly focused on the presence, absence or semi-quantitative analysis of immune cells in tumour epithelium.⁷ Moreover, cell densities of lymphocytes, macrophages and dendritic cells have also been shown to be prognostic in breast, ovarian and lung cancer without restricting analyses to malignant epithelial cell areas.⁸⁻¹⁰ Approaches that combine analyses of multiple immune infiltrates, such as the CD3/CD8 immunoscore in colorectal cancer,¹¹ have not yet been developed for ovarian cancer prognosis.

The notion that anti-tumoural immune responses may be altered by tumour cell-intrinsic factors is supported by several observations. Detailed temporal and spatial histological and genomic studies in a single patient with high-grade serous

ovarian carcinoma (HGSOC) showed that Wnt signalling was upregulated in a progressing tumour nodule, a phenomenon correlated with focal immunosuppression.¹² Mutations in the *TP53* and *BRCA1/BRCA2* genes as well as loss of *PTEN* expression are driver events in HGSOC development.¹³ HGSOC cases with *BRCA1* mutations have increased CD8 and CD20 intra-epithelial infiltrates,¹⁴,¹⁵ suggesting that loss of homologous recombination and DNA damage may prime immune responses. Whether different classes of *TP53* mutations alter immune infiltration in HGSOC remains unknown, but non-synonymous mutations in *TP53* may have gain-of-function or other cellular effects distinct from loss-of-function mutations.¹⁶ Mutant p53 protein may drive B cell responses and auto-antibody production.¹⁷ Loss of *PTEN* expression in melanoma was associated with both reduced T cell infiltration and resistance to immune checkpoint inhibitors,¹⁸ but these associations have not been investigated in HGSOC.

Despite the strong association of CD8⁺ infiltrate with prognosis, routine immunoscoring for ovarian cancer is not performed in the clinic and development of automatic cellular recognition tools

<footnote n="1">¹Barts Cancer Institute, Queen Mary University of London, London, UK; ²Cancer Research UK Cambridge Institute, Cambridge, UK; ³Department of Oncology, University of Cambridge, Cambridge, UK; ⁴Laboratory of Translational Oncology, Intercollegiate Faculty of Biotechnology, Medical University of Gdańsk, 80-211 Gdańsk, Poland; ⁵Department of Public Health and Primary Care, Centre for Cancer Genetic Epidemiology, University of Cambridge, Cambridge, UK; ⁶Department of Pathology, University College London Hospital, London, UK and ⁷Cambridge University Hospitals NHS Foundation Trust, Cambridge, UK</footnote>

<footnote>Correspondence: Frances R. Balkwill (f.balkwill@qmul.ac.uk) or James D. Brenton (James.Brenton@cruk.cam.ac.uk)</footnote>

<footnote>These authors contributed equally: Anne Montfort, Rowan J. Barker-Clarke</footnote>
could be beneficial for high-throughput pathology workflows. To test the hypothesis that a more integrated analysis of cytotoxic, mature and organised immune responses in the tumour micro-environment might have greater prognostic value in HGSOC, we developed new image analysis methods and workflows to test the single and integrated analysis of CD8⁺, CD45RO⁺ and CD68⁺ in the malignant epithelium and adjacent stroma. The secondary aim of this study was to test for association between quantitative measures of CD8⁺, CD45RO⁺ and CD68⁺ cells and driver genomic alterations in *BRCA1*, *BRCA2*, *PTEN* and different classes of *TP53* mutation.

# MATERIALS AND METHODS

## Patients

Samples from 570 patients from the prospective SEARCH ovarian cancer population-based study were used to construct tissue microarrays (TMAs).¹⁹ Ethical approval was granted by the Eastern Multicenter Research Ethics Committee. Among the samples from 570 patients with primary epithelial ovarian tumours, 332 were high-grade serous ovarian cancer patients. All cases underwent detailed histopathological review by a gynaecological pathologist (J.Mc.D.). Patients were staged as having localised, regional or distant disease (L/R/D).²⁰

## Immunohistochemistry

Microarray slides composed of formalin-fixed paraffin-embedded ovarian tumour cores were dewaxed and rehydrated prior to heat-induced epitope retrieval using a pressure cooker and a citrate-based antigen unmasking solution (Vector Laboratory). Detection of CD8⁺ T cells, CD45RO⁺ memory lymphocytes and CD68⁺ macrophages was performed using the mouse anti-human CD8 (clone C8/144B, Dako), mouse anti-human CD45RO (clone UCLH, Dako) and mouse anti-human CD68 (clone M0876, Dako) antibodies, using ultrasensitive Polymer-HRP IHC Detection system (Biogenex). Immunohistochemical protocols and slide hybridisations were carried out manually. Sections were counterstained with haematoxylin and mounted with DPX mounting medium (Sigma). Previously published PTEN immunostaining data was used where high PTEN expression was considered to be positive staining and low expression to be weak, heterogeneous or negative staining, respectively.²¹

## Mutation analysis

The coding regions of *TP53* were sequenced by tagged-amplicon next-generation sequencing as previously described²² and confirmed by immunohistochemical analysis using a 4-tier core system.²³ Sequencing of germline mutations in the *BRCA1* and *BRCA2* genes was performed as previously described.²⁴

## Immune cell quantification

Stained slides were scanned using the Panoramic Slash Scanner (3D Histech). The number of CD8⁺ and CD45RO⁺ cells per mm² of epithelial and stroma areas, as well as the percentage of epithelial and stromal areas covered by CD68 staining, were digitally determined using the Tissue Studio software (Definiens™). Definiens image analysis algorithms for detection of epithelial and stromal areas were trained and the segmentation for each core was manually refined by two researchers, including a consultant gynaecological-histopathologist (J.Mc.D.). Supplementary Fig. 1 shows examples of classifications of tissue regions and cell detection and the entire data set, including these epithelial and stromal assignments, can be downloaded from the repository.

## Statistical analyses

R (version 3.5.1) was used for statistical analysis and an R markdown document containing the entire data set allowing for performing all analyses is available (https://bitbucket.org/

jamesdbrenton/search-montfort/src/master/). Quality checking for spatial bias across TMAs and effects of varying tissue area was carried out upon all cores and across TMAs using heatmaps and Shapiro-Wilk tests. All count data were transformed to log base 10 after adding a small offset to zero values. Wilcoxon's signed-rank test was used to compare the mean infiltrate between groups. Continuous data were presented as median and interquartile range (IQR) and groups were compared by the Kruskal-Wallis and pairwise Kruskal-Wallis tests. Discrete data were presented as count and percentage.

Cox proportional hazard regression analysis was applied to assess the effect of each infiltrate on overall survival. The functional form of each of the immune variables was assessed using comparison with cubic splines. The best approximations to the functional forms were carried forward for the Cox models. The clinical variables of age at diagnosis, menopause status and stage were available for the cohort and were included in the analysis. Univariable Cox regressions were used to identify best-fitting variables for the final multivariable Cox regression model. The refined model was compared with a combined multivariable Cox regression model including all immune infiltrates. Hazard ratios (HRs) refer to a single unit increase in continuous variables. The proportional hazards assumption was tested and satisfied in all cases using Schoenfeld residuals. The Kaplan-Meier analysis (with log-rank test) was applied to illustrate survival differences graphically. Two-sided *p* values <0.05 were used to indicate statistical significance. Principal component analysis (PCA) using the R package prcomp was used to extract the independent components of variance between patients. The package prcomp uses singular value decomposition and the variables were scaled to have unit variance before creating composite linear independent variables. These were then passed forward to the survival modelling. The Akaike information criterion (AIC) was used to compare the performance of survival models, which includes a penalty on the number of terms to reduce overfitting. Bonferroni *p* value corrections were carried out for all multiple testing. *P* < 0.05 was considered significant for all analyses.

# RESULTS

## Patient characteristics

Supplementary Fig. 2 shows the REMARK diagram for this study and Supplementary Table 1 shows the clinical characteristics of the 332 HGSOC patients from the study cohort. Immunohistochemical analyses on TMAs were performed to detect CD8⁺, CD45RO⁺ and CD68⁺ cells in tissue cores from primary ovarian specimens. One hundred and fifty-two HGSOC cases were available for analysis after quality assurance, data cleaning and the reduction of the data set to only cases with complete results for CD8, CD45RO and CD68 staining in both epithelium and stroma, as well as survival data.

Tagged-amplicon sequencing was performed on 248 cases and *TP53* mutation was detected in 231 samples (93%) (Supplementary Table 1). Previously published data for germline *BRCA1* and *BRCA2* mutation and PTEN expression were available for 297 and 155 cases, respectively.¹⁸,²²

## Digital pathology analysis of tumour composition and immune cell densities

Image analysis software was used to determine the area of tumour epithelium and stroma in each core (Fig. 1 and Supplementary Fig. 1). Of 964 images representing 332 HGSOC patients, 69 patients (20.8%) had images that contained malignant epithelium but no stroma; 250 patients (75.3%) had images that contained epithelium and >1% adjacent stroma and 13 patients (3.9%) had images containing no tumour epithelium (Fig. 1a). The median proportion of epithelium and stroma was 85.1% (IQR 51–100%) and 14.9% (IQR 0–49%), respectively. We expected the proportion
<figure><img src="image_2.png" /><figcaption>**Fig. 1** Proportions of immune cells vary between epithelial and stroma areas of tumours. **a** Scatter plot of the average stroma and epithelial tumour areas for each patient. **b–d** show respectively the distribution of densities of CD8<sup>+</sup>, CD45RO<sup>+</sup> and CD68<sup>+</sup> cell in epithelium and stroma. CD8<sup>+</sup> and CD45RO<sup>+</sup> densities were defined as counts per mm<sup>2</sup> and CD68<sup>+</sup> as the percentage of tissue stained for this marker. Notches on box plots extend 1.58 × IQR/sqrt(n) and approximate the 95% confidence interval for the median. Box plot whiskers extend to 1.5 × IQR.</figcaption></figure>

of tumour in a sample to be correlated with p53 mutant allele fraction, a measure of sample purity, and found them to be positively correlated ($R^2 = 0.25, p = 0.0004$).

Quantitative densities of all immune populations were then generated automatically through image analysis, the quantitative assessment of CD8<sup>+</sup> T cell, CD45RO<sup>+</sup> memory lymphocyte and CD68<sup>+</sup> macrophage densities in each compartment are shown in Supplementary Fig. 1. The relationship between the fraction of tumour in a core and the density of immune infiltration in the epithelium was examined. Intra-epithelial CD8<sup>+</sup> and CD45RO<sup>+</sup> densities were weakly correlated with the purity/tumour fraction of the sample ($R^2 = 0.17, p = 0.003$ and $R^2 = 0.16, p = 0.006$), but CD68<sup>+</sup> epithelial density was not.

The three immune populations in our samples showed moderate to strong correlation between epithelium and stroma (Supplementary Fig. 3 and Supplementary Table 2). Samples with low density of stromal immune populations generally had low density of epithelial infiltrate and vice versa. The distribution of densities of immune populations within tumour epithelium and stromal areas were compared (Fig. 1b–d). The density of CD8<sup>+</sup> and CD45RO<sup>+</sup> cells were significantly higher in stroma than in tumour epithelium ($p = 0.005$ and $p = 0.004$, respectively; Welch's t test), but not significantly different for CD68<sup>+</sup> cells.

In describing the patterns of immune infiltrate, the terms immune inflamed, immune desert and immune excluded have

been used to describe varying T cell infiltration based on histological and transcriptional analyses.<sup>12,25–27</sup> Immune-inflamed and immune-desert patterns reflect high positive or negative correlations between all infiltrates, but T cell exclusion describes tumours where CD8<sup>+</sup> cells are significantly absent from tumour epithelium while still being present in the surrounding stroma.<sup>28,29</sup> Given the higher infiltration in stroma than epithelial compartments, we defined immune cell exclusion as a 10-fold difference between tumour epithelium and stromal infiltration as the standard deviation of the log 10-transformed counts was ~1. CD8<sup>+</sup> exclusion was present in 20 (10.6%) cases and 36 (20%) cases had CD45RO<sup>+</sup> exclusion. No cases had significant exclusion of CD68<sup>+</sup> infiltrate from tumour epithelium. Notably none of the cases had both CD8<sup>+</sup> and CD45RO<sup>+</sup> exclusion.

Stromal CD68<sup>+</sup> and CD45RO<sup>+</sup> densities are the strongest individual prognostic markers

Survival was modelled using Cox proportional hazards and the relationship between the immune variables and survival was found to be approximately log linear. The clinical variables accompanying the cohort were age at diagnosis, stage and menopause status and the relationship between age and survival was found to be approximately linear (see Methods, Supplementary Fig. 4 and Supplementary Table 3).
<table><caption>**Table 1.** Hazard ratios from the Cox proportional model for all infiltrates in all regions and averaged across the whole core in HGSOC, measured as log 10 (counts per mm²).</caption><thead><tr><th rowspan="2">Functional form</th><th rowspan="2">Evaluable cases</th><th rowspan="2">Tissue compartment</th><th colspan="2">Univariable</th><th colspan="2">Multivariable<sup>a</sup> (adjusted for stage)</th></tr><tr><th>HR</th><th><i>p</i> Value</th><th>HR</th><th><i>p</i> Value</th></tr></thead><tbody><tr><td>CD8<sup>+</sup></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Log 10</td><td>301</td><td>Epithelium</td><td>0.89</td><td>0.15</td><td>0.83</td><td>0.027</td></tr><tr><td>Log 10</td><td>202</td><td>Stroma</td><td>0.97</td><td>0.74</td><td>0.93</td><td>0.40</td></tr><tr><td>Log 10</td><td>315</td><td>Average</td><td>0.79</td><td>0.010</td><td>0.72</td><td>0.0006</td></tr><tr><td>CD45RO<sup>+</sup></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Log 10</td><td>290</td><td>Epithelium</td><td>0.86</td><td>0.033</td><td>0.85</td><td>0.022</td></tr><tr><td>Log 10</td><td>196</td><td>Stroma</td><td>0.76</td><td>0.001</td><td>0.76</td><td>0.0007</td></tr><tr><td>Log 10</td><td>306</td><td>Average</td><td>0.82</td><td>0.006</td><td>0.80</td><td>0.003</td></tr><tr><td>CD68<sup>+</sup></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Linear</td><td>293</td><td>Epithelium</td><td>0.99</td><td>0.67</td><td>0.99</td><td>0.43</td></tr><tr><td>Log 10</td><td>226</td><td>Stroma</td><td>0.53</td><td>0.003</td><td>0.44</td><td>0.0003</td></tr><tr><td>Log 10</td><td>308</td><td>Average</td><td>0.67</td><td>0.042</td><td>0.62</td><td>0.017</td></tr><tr><td>Stage</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>-</td><td>312</td><td>Localised</td><td>1</td><td>0</td><td>1</td><td>0</td></tr><tr><td></td><td></td><td>Regional</td><td>1.47</td><td>0.26</td><td>1.15</td><td>0.25</td></tr><tr><td></td><td></td><td>Distant</td><td>3.96</td><td>&lt;0.001</td><td>5.58</td><td>&lt;0.001</td></tr><tr><td></td><td></td><td>Unstaged</td><td>3.35</td><td>&lt;0.001</td><td>3.34</td><td>&lt;0.001</td></tr></tbody></table>

Univariable analysis showed improved survival with increasing stromal density of CD45RO⁺ [HR 0.76, 95% confidence interval (CI): 0.65–0.90, *p* = 0.001] and CD68⁺ (HR 0.53, 95% CI: 0.34–0.81, *p* = 0.003) (Table 1). Modelling each immune variable with stage showed improved predictive value for epithelial CD8⁺ density (HR = 0.83, *p* value=0.027) as well as stromal CD68⁺ and CD45RO⁺ density and epithelial CD45RO⁺ density (Table 1). Figure 2 and Supplementary Fig. 5 shows illustrative Kaplan–Meier survival curves for high and low stromal and epithelial CD68⁺, CD45RO⁺ and CD8⁺ densities.

In cores with <1% stroma, epithelial CD8⁺ malignant epithelial infiltrate remained an independent prognostic factor, but epithelial CD45RO⁺ density was not significant (Supplementary Table 4).

In clinical reporting, quantifying immune populations in exclusively tumour epithelium is technically challenging and time consuming. We tested the effects of using the average density of each marker averaged across both tumour and stromal regions from each core (Table 1). Averaging the tissue density of CD8⁺ increased the strength of the associated HR and significance of the model (HR = 0.79, *p* value = 0.010) indicating increased prognostic value over quantitation of individual epithelial and stromal infiltrates. Supplementary Fig. 5 shows illustrative Kaplan–Meier survival curves for high and low CD68⁺, CD45RO⁺ and CD8⁺ densities over combined epithelium and stroma compartments.

Multivariable Cox regression analysis including all infiltrates and stage was carried out on patients with complete data for all infiltrates (n = 152) (Table 2). In this model, only stage and CD68⁺ stromal infiltrate were significant predictors of survival.

We then refined the model by removing the least significant elements (defined as those with *p* > 0.1) (Table 2, n = 152). Interestingly, we found that the *p* values for CD68⁺ and CD45RO⁺ stromal infiltrates in the refined model become less significant and the hazard ratios are attenuated in comparison to both the univariable regression and the full model. This is likely due to the inability of Cox regression to distinguish with confidence whether stromal CD45RO⁺ or CD68⁺ density is the most significant

predictor when there is a strong correlation between all the immune variables.

Principal components of the combined immunospace describe biologically interpretable effects

As the three types of immune infiltrate vary continuously across epithelium and stroma these variables can be regarded as six dimensions of an 'immunospace' (three infiltrates, two localisations). Given the strong correlations between infiltrates (Supplementary Fig. 3), these immune variables are not independent. We used PCA to determine the independent patterns across these dimensions, using the 152 patients for whom complete data were available. PCA transformed the six correlated infiltrate variables into six independent axes with the first component containing the largest proportion of variance (60%) in the data set (Supplementary Table 5).

In principal component 1 (PC1), the weightings of all immune infiltrates are positive and similar in magnitude (Supplementary Table 5). This indicates that as one infiltrate increases so do all the others and represents the degree of coordinated immune response. The remaining PCs characterise additional patterns across immune infiltrates independent of PC1. The additive contribution of PC2 characterises negative correlation between CD8⁺ infiltrates and CD68⁺ macrophages and CD45RO⁺ memory cells. PC3 characterises additional variation where epithelial and stromal infiltrates are negatively correlated, the most positive values of PC3 correspond to high infiltration in tumour epithelium compared to stroma and the most negative values of PC3 correspond to the opposite, the aforementioned immune exclusion.

Supplementary Fig. 6 shows representative images with the largest magnitudes of PC1, PC2 and PC3 to visually illustrate the patterns described above. The variance in the remaining PCs (4–6) is smaller and less informative. Variance in PC4 is predominantly from CD45RO⁺ stromal density, PC5 is from CD45RO⁺ epithelial density and PC6 is from CD68⁺ epithelial density. Patients are plotted by their PC1 and PC2 values in Supplementary Fig. 7. Cox

<footnote><sup>a</sup>Multivariable analysis was adjusted for stage.</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2** Survival analysis of HGSOC patients relative to the density of CD68⁺ macrophages and CD45RO⁺ cells in their stroma. Kaplan-Meier survival curves using cut point of median density of stromal CD68⁺ macrophages and CD45RO⁺ memory T cells using left truncation and right censoring. Median entry to the study for all patients after diagnosis was 26.4 months. Median follow-up time from diagnosis to exit or death was 105.1 months.</figcaption></figure>

<table><caption>**Table 2.** Multivariable Cox regression hazard ratios and associated *p* values for a model with all infiltrates and stage and for a reduced model with most significant variables only (*n* = 152).</caption><thead><tr><th rowspan="2"></th><th colspan="2">Multivariable (all combined)</th><th colspan="2">Refined model</th></tr><tr><th>HR</th><th><i>p</i> Value</th><th>HR</th><th><i>p</i> Value</th></tr></thead><tbody><tr><td>CD8<sup>+</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>Epithelium</td><td>0.96</td><td>0.81</td><td>-</td><td>-</td></tr><tr><td>Stroma</td><td>1.07</td><td>0.63</td><td>-</td><td>-</td></tr><tr><td>CD45RO<sup>+</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>Epithelium</td><td>1.12</td><td>0.37</td><td>-</td><td>-</td></tr><tr><td>Stroma</td><td>0.83</td><td>0.09</td><td>0.68</td><td>0.11</td></tr><tr><td>CD68<sup>+</sup></td><td></td><td></td><td></td><td></td></tr><tr><td>Epithelium</td><td>1.16</td><td>0.63</td><td>-</td><td>-</td></tr><tr><td>Stroma</td><td>0.53</td><td>0.038</td><td>0.88</td><td>0.17</td></tr><tr><td>Stage</td><td></td><td></td><td></td><td></td></tr><tr><td>Localised</td><td>1</td><td>0</td><td>1</td><td>0</td></tr><tr><td>Regional</td><td>2.00</td><td>0.16</td><td>2.03</td><td>0.14</td></tr><tr><td>Distant</td><td>4.82</td><td>&lt;&lt;0.001</td><td>4.70</td><td>0.0001</td></tr><tr><td>Unstaged</td><td>8.15</td><td>&lt;&lt;0.001</td><td>8.25</td><td>0.0001</td></tr></tbody></table>

proportional hazard regression was used to assess whether these PCs were predictive of survival. Only PC1 was an independent predictor of survival in our cohort and was associated with improved survival (univariate: HR = 0.89, *p* value = 0.024; PC1 + stage: HR: 0.88, *p* value = 0.019) (Supplementary Table 6), reflecting the good prognosis of a strong coordinated immune response. The association of this PC with survival is also illustrated graphically using Kaplan-Meier curves in Fig. 3.

Cox regression survival models were also calculated on all combinations of PCs and stage. The AIC was used to compare the performance of these survival models and includes a penalty on the number of terms to reduce overfitting (Supplementary Fig. 8). The model combining stage, PC1 and PC5, had the best performance for predicting overall survival. The improvement with the addition of PC5 shows that the addition of this PC has a suppressor effect in the model, increasing the significance of other variables when included. This demonstrates that survival is predominantly determined by the coordinated immune response and further variation in survival from this trend can be encoded by the quantity of epithelial CD45RO⁺ infiltrate. Interestingly, the models that contained stage and either stromal CD45RO⁺ or

<figure><img src="image_4.png" /><figcaption>**Fig. 3** Survival analysis by principal component 1 (PC1). Kaplan-Meier analysis using median cut point for PC1 with left truncation for study entry variables and right censored at 120 months. Median entry to the study for all patients after diagnosis was 26.4 months. Median follow-up time from diagnosis to exit or death was 105.1 months.</figcaption></figure>

CD68⁺ infiltrate contained a similar amount of information about patient survival as the one that contained stage and PCs 1 and 5. In our cohort, the density of CD68⁺ and CD45RO⁺ stromal infiltrates are therefore the best single infiltrates for survival modelling.

Are genetic defects associated with HGSOC driving individual infiltrates or the coordinated immune response in tumours? The presence of germline *BRCA2* mutations was significantly associated with lower CD8⁺ cell density than patients with a *BRCA1* after multiple testing correction (Table 3). There was no significant association between the quantity of CD45RO⁺ or CD68⁺ infiltrate and the mutational status of either *BRCA1* or *BRCA2* genes (Table 3). No significant association was detected between *TP53* gain-of-function (GOF) and loss-of-function (LOF) mutations or PTEN expression and the densities of CD8⁺, CD45RO⁺ and CD68⁺ cells in the epithelium or stroma (Table 3). Similarly, changes in the PCs were not significantly associated with PTEN expression, *TP53* GOF or LOF mutation or germline *BRCA1/BRCA2* mutation status (Table 3).
<table><caption>**Table 3.** Association between immune infiltrate and genetic alterations in HGSOC tumours.</caption><thead><tr><th></th><th>BRCA1/BRCA2</th><th>p53 GOF/LOF</th><th>PTEN</th></tr></thead><tbody><tr><td>CD8<sup>+</sup></td><td></td><td></td><td></td></tr><tr><td>Epithelium</td><td>0.04</td><td>0.54</td><td>0.67</td></tr><tr><td>Stroma</td><td>0.17</td><td>0.29</td><td>0.07</td></tr><tr><td>Average</td><td>0.17</td><td>0.53</td><td>0.84</td></tr><tr><td>CD45RO<sup>+</sup></td><td></td><td></td><td></td></tr><tr><td>Epithelium</td><td>0.20</td><td>0.42</td><td>0.92</td></tr><tr><td>Stroma</td><td>0.46</td><td>0.84</td><td>0.82</td></tr><tr><td>Average</td><td>0.10</td><td>0.36</td><td>0.87</td></tr><tr><td>CD68<sup>+</sup></td><td></td><td></td><td></td></tr><tr><td>Epithelium</td><td>0.43</td><td>0.72</td><td>0.79</td></tr><tr><td>Stroma</td><td>0.78</td><td>0.85</td><td>0.65</td></tr><tr><td>Average</td><td>0.68</td><td>0.66</td><td>0.41</td></tr><tr><td>Principal component</td><td></td><td></td><td></td></tr><tr><td>1</td><td>0.81</td><td>0.80</td><td>0.55</td></tr><tr><td>2</td><td>0.10</td><td>0.26</td><td>0.17</td></tr><tr><td>3</td><td>0.27</td><td>0.58</td><td>0.51</td></tr><tr><td>4</td><td>0.58</td><td>0.99</td><td>0.17</td></tr><tr><td>5</td><td>0.76</td><td>0.88</td><td>0.79</td></tr><tr><td>6</td><td>0.97</td><td>0.44</td><td>0.18</td></tr></tbody></table>

P values were associated with Kruskal-Wallis test for detecting differences in mean ranks of immune infiltrate in patients grouped by mutation in BRCA1, BRCA2 or not-detected, p53 GOF or LOF and PTEN high or low in HGSOC.

## DISCUSSION

Infiltration of immune cells in HGSOC tumours has previously been correlated to prognosis for patients.¹,³⁰ However, translation of these tools to the clinic has been impeded by the lack of standardised and reliable quantitation methods. In this work, we used a semi-automated approach to quantify the CD8⁺, CD45RO⁺ and CD68⁺ infiltrates in both stromal and epithelial areas of primary HGSOC tumours. We found that immune infiltration in the tumour microenvironment was continuously distributed across a wide range and CD8⁺, CD45RO⁺ and CD68⁺ infiltrates were strongly correlated. T cell exclusion from epithelial regions has been previously reported.²⁹ We only observed epithelial exclusion of CD8⁺ or CD45RO⁺ in 10% and 20% of patients, respectively. These groups were mutually exclusive, suggesting that immune exclusion from epithelial or stromal regions, is a tumour-specific effect that is distinct from coordinated multi-infiltrate immune responses.

Our results for the positive prognostic effects of epithelial CD8⁺ T cells are consistent with the large study by Goode et al.,² who described a near log-linear relationship between the density of epithelial CD8⁺ T cells and survival. However, our observations also show that averaging CD8 density across the total core (epithelial and stromal areas) improves survival prediction. The alternative possibility that epithelial and stroma compartments have been mis-assigned is highly unlikely as these assignments were all reviewed by a gynaecological pathologist. The combination of both stroma and tumour epithelium potentially provides a better representation of the dynamics of epithelial infiltration by including cells that may be poised to infiltrate.

CD45RO is a memory marker for T and B cells. In the tumour microenvironment, we and others⁵,³¹,³² have shown that CD4⁺ and CD8⁺ T cells, myeloid cells and B cells are mainly located in the stroma of tumours, whereas CD8⁺ T cells can be found in both areas. Meta-analyses investigating the prevalence of CD45RO in

solid cancers reported that high CD45RO infiltrate was associated with better prognosis for patients.³³ We confirmed these results in HGSOC and demonstrated that the quantity of stromal and average densities of CD45RO⁺ cells over the whole core are prognostic. We also demonstrated that epithelial CD45RO⁺ density was a significant predictor of survival in cores with >1% stroma, but not in those with <1% stroma. This result emphasises the need to contextualise epithelial infiltrate and the confounding nature of tumour composition on survival impact.

The link between tumour infiltration with macrophages and patient survival is more complex. Previous studies demonstrated the positive predictive value of classically activated (M1) over tumour-promoting (M2) macrophages in ovarian cancer.⁹,³⁴ However, no significant association was found between the level of intra-epithelial CD68 infiltrate and patient survival.⁴ Our study confirms this latter result and also demonstrates a significant association between a high coverage of stromal areas with CD68⁺ macrophages and improved overall survival. This result is in line with a recent study showing a high infiltration of CD68⁺ macrophages at the invasive front of tumour sections from colorectal cancer patients is correlated to better response to chemotherapy.³⁵ Whereas tumour-associated macrophages (TAM) promote cancer progression via stimulation of multiple processes including angiogenesis, inflammation and immune escape, it now appears evident that the phenotype of TAM changes in different areas of the tumour.³⁶,³⁷ Taking into account the seemingly contradictory results obtained on this subject, additional studies aiming at deciphering the precise role distinct subpopulations of macrophages infiltrating different areas of tumours play in cancer progression is required.

All three infiltrates were also found to be predictive in multivariate analyses adjusted for tumour stage. We previously demonstrated that optimal response to neoadjuvant chemotherapy was significantly associated with a decrease in the density of FoxP3⁺ regulatory T cells in the stroma of human HGSOC tumours.³⁸ As we did not have access to treatment information for patients included in this cohort, we cannot rule it out as a potential confounding factor.

Importantly, we show that combining our six correlated types of immune infiltrates using PCA allowed us to transform inter-patient variations of the tumour immune landscape into independent, biologically interpretable PCs. We find that the main source of variation between patients is reflective of the quantity of concerted immune response. The second PC demonstrates that there is a pattern of significant variations in the CD8⁺ infiltrate of some patients that is independent and acting in opposition to the concerted immune response. The third PC measures the extent of negatively correlated epithelial and stromal infiltrates and is weighted differently by infiltrates. This supports our observation that immune exclusion is common, variable by infiltrate and shows that this is an additional effect to that of the concerted immune response.

The main limitation of our analysis of the 'immunospace' and its PCs is that we could only apply it to images that contain both epithelium and stroma. This reduced the size of the cohort we could analyse. With foresight, sampling could be designed to include both regions in TMAs and sectioning and imputation could be carried out for missing data. It is worth noting that being able to combine multiple correlated infiltrates and elucidate patterns and sources of variation will become even more important as many more (30+) immune markers are combined on single sections using imaging methods such as Hyperion. This method is particularly useful in that it also allows us to measure the strength of immune patterns occurring across multiple infiltrates. Analysing PCs also avoids some statistical issues that are associated with frequent co-correlation between different immune populations. Including correlated variables in a typical multivariable Cox regression survival model not only reduces
interpretability but can also cause model results to be variable under bootstrapping as seen in the refined model of our cohort.

We found that the first PC was independently associated with survival. This result is consistent with a previous study showing that the simultaneous infiltration of tumours by different subsets of leucocytes (e.g. T cells, plasma cells, B cells), likely reflecting the establishment of a concerted immune response, gives a survival advantage to ovarian cancer patients.⁵ We used the PCs to model survival and found that a model including PC1 and PC5 was as predictive as individually modelled stromal CD45RO and CD68 infiltrates. Therefore, CD68 or CD45RO markers analysed in epithelium-adjacent stroma constitute the best single prognostic markers in our cohort.

In human lung adenocarcinoma, Mansuet-Lupo et al.³⁹ showed that oncogenic mutations in the *TP53* gene were positively correlated with CD8⁺ infiltrate. In their study, intra-epithelial CD8⁺ T cell numbers and *TP53* status were both linked to prognosis with patients harbouring tumours with low CD8⁺ infiltrate and non-disruptive *TP53* mutation (associated with GOF) being linked to poorer survival. In glioblastoma, *TP53* mutations were also associated with increased immune infiltrate.⁴⁰ In the present work, we integrated histological and genomic features to evaluate whether genetic alterations could be linked to different quantities of immune infiltration. We found no correlation between the nature of *TP53* mutations and the amount of CD8⁺ T cells, CD45RO⁺ memory lymphocytes and CD68⁺ macrophages in the stromal and malignant areas of 197 HGSOC tumours or the values of the PCs we derived. Overall, the exact nature of the relationship between subtypes of *TP53* mutations and immune infiltrate is still unclear and likely to vary across different cancers, as evidenced by two studies showing *TP53* GOF mutations to be associated with pro-tumour effects related to inflammation in glioblastoma⁴⁰ and immunosuppression in lung adenocarcinoma.⁴¹

PTEN, another essential tumour suppressor, regulates the production of immunosuppressive cytokines by melanoma cells.⁴² In human melanoma, *PTEN* deletion was correlated with a decrease in infiltrating CD8⁺ T cells.¹⁸ In HGSOC, however, we observed no correlation between the expression level of PTEN and the number and/or localisation of CD8⁺, CD45RO⁺ and CD68⁺ leucocytes in the tumour microenvironment or the PCs we derived. The discrepancy between the results obtained with HGSOC and melanoma tumours might lie in the fact that in our study loss of PTEN expression was only partial with tumours divided into expressing high or low levels of PTEN.

In contradiction to a previous work,⁴³ we did not observe a significant correlation between the number of intra-epithelial CD8⁺ T cells and mutations in the *BRCA1* gene. As there were only 18 patients with a *BRCA1* mutation in this cohort, this is most likely related to a broad confidence interval as compared with the null population. We did however find a significant difference in intra-epithelial CD8⁺ T cell infiltration between patients with *BRCA1* and *BRCA2* mutations, with *BRCA1* patients having significantly higher infiltration than *BRCA2* cases. The discrepancy could also be linked to the automated and continuous quantitation method used in our study to quantify the immune infiltrate of tumours, as scoring in the previous study was done manually and in a stratified manner.⁴³ Furthermore, we did not assess the methylation in the gene coding region for *BRCA1* and restricted our analysis to the identification of germline mutations in the *BRCA1* and *BRCA2* genes.

Our work shows that averaging immune infiltration over the whole tissue core could be as useful as current methods and potentially provides more prognostic information. Methods that average counts of immune cells across tissues are also simpler to implement than epithelial:stromal segmentation methods. Nonetheless, careful sampling of the tumour microenvironment with inclusion of both epithelium and stroma remains very important as both areas have independent prognostic significance.

Our results showing strong positive prognostic significance of stromal CD68⁺ infiltrate in HGSOC tumours warrants further investigation into the role and properties of CD68⁺ macrophages in HGSOC and also may caution against the use of new potential macrophage-depleting therapies.

## ACKNOWLEDGEMENTS

We are grateful to the patients who consented to donate tumour samples and to Dr. Filippe Correia Martins for helpful discussions about the PTEN analysis.

## AUTHOR CONTRIBUTIONS

Experiments and data analysis: A.M., R.J.B.C., S.B., J.McD., P.P., J.D.B., F.R.B., A.M.P., A.S., L.M., S.A.-K.; TMA cohort: P.P.; writing paper: J.D.B., F.R.B., A.M., R.J.B.C.; designed study and acquired funding: F.R.B. and J.D.B.

## ADDITIONAL INFORMATION

**Ethics approval and consent to participate** Ethical approval for the collection of tumour biopsies from ovarian cancer patients was granted by the Eastern Multicenter Research Ethics Committee. All patients enrolled in this study provided informed consent for the use of their tissue samples for research purposes. This study was performed in accordance with the Declaration of Helsinki.

**Data availability** The R markdown document containing the entire data set and allowing for reproducing all analyses performed in this manuscript is available online (https://bitbucket.org/jamesdbrenton/search-montfort/src/master/).

**Competing interests** P.P. is a member of the British Journal Cancer Editorial Board. All other authors declare no competing interests.

**Funding information** We acknowledge funding and support from Cancer Research UK (grant numbers A16354, A25714, A15973, A15601, A22905, A17197), ERC Advanced Grant ERC322566, the Queen Mary University of London, the National Institute for Health Research Cambridge, the National Cancer Research Network and the Cancer Research UK Experimental Cancer Medicine Centres and Hutchison Whampoa Limited. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41416-020-0822-x.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## REFERENCES

1. Zhang, L., Conejo-Garcia, J. R., Katsaros, D., Gimotty, P. A., Massobrio, M., Regnani, G. et al. Intratumoral T cells, recurrence, and survival in epithelial ovarian cancer. *N. Engl. J. Med.* **348**, 203–213 (2003).

2. Goode, E. L., Block, M. S., Kalli, K. R., Vierkant, R. A., Chen, W., Fogarty, Z. C. et al. Dose-response association of CD8⁺ tumor-infiltrating lymphocytes and survival time in high-grade serous ovarian cancer. *JAMA Oncol.* **3**, e173290 (2017).

3. Zhang, Z., Huang, J., Zhang, C., Yang, H., Qiu, H., Li, J. et al. Infiltration of dendritic cells and T lymphocytes predicts favorable outcome in epithelial ovarian cancer. *Cancer Gene Ther.* **22**, 198–206 (2015).

4. Milne, K., Köbel, M., Kalloger, S. E., Barnes, R. O., Gao, D., Gilks, C. B. et al. Systematic analysis of immune infiltrates in high-grade serous ovarian cancer reveals CD20, FoxP3 and TIA-1 as positive prognostic factors. *PLoS ONE* **4**, e6412 (2009).

5. Kroeger, D. R., Milne, K. & Nelson, B. H. Tumor-infiltrating plasma cells are associated with tertiary lymphoid structures, cytolytic T-cell responses, and superior prognosis in ovarian cancer. *Clin. Cancer Res.* **22**, 3005–3015 (2016).

6. Leffers, N., Gooden, M. J., de Jong, R. A., Hoogeboom, B. N., ten Hoor, K. A., Hollema, H. et al. Prognostic significance of tumor-infiltrating T-lymphocytes in primary and metastatic lesions of advanced stage ovarian cancer. *Cancer Immunol. Immunother.* **58**, 449–459 (2009).

7. Fridman, W. H., Zitvogel, L., Sautés-Fridman, C. & Kroemer, G. The immune context in cancer prognosis and treatment. *Nat. Rev. Clin. Oncol.* **14**, 717–734 (2017).
8. Ali, H. R., Dariush, A., Thomas, J., Provenzano, E., Dunn, J., Hiller, L. et al. Lymphocyte density determined by computational pathology validated as a predictor of response to neoadjuvant chemotherapy in breast cancer: secondary analysis of the ARTemis trial. *Ann. Oncol.* **28**, 1832–1835 (2017).

9. Zhang, M., He, Y., Sun, X., Li, Q., Wang, W., Zhao, A. et al. A high M1/M2 ratio of tumor-associated macrophages is associated with extended survival in ovarian cancer patients. *J. Ovarian Res.* **7**, 19 (2014).

10. Dieu-Nosjean, M. C., Antoine, M., Danel, C., Heudes, D., Wislez, M., Poulot, V. et al. Long-term survival for patients with non-small-cell lung cancer with intratumoral lymphoid structures. *J. Clin. Oncol.* **26**, 4410–4417 (2008).

11. Pagès, F., Kirilovsky, A., Mlecnik, B., Asslaber, M., Tosolini, M., Bindea, G. et al. In situ cytotoxic and memory T cells predict outcome in patients with early-stage colorectal cancer. *J. Clin. Oncol.* **27**, 5944–5951 (2009).

12. Jiménez-Sánchez, A., Memon, D., Pourpe, S., Veeraraghavan, H., Li, Y., Vargas, H. A. et al. Heterogeneous Tumor-Immune Microenvironments among Differentially Growing Metastases in an Ovarian Cancer Patient. *Cell* **170**, 927–938.e920 (2017).

13. Bowtell, D. D., Böhm, S., Ahmed, A. A., Aspuria, P. J., Bast, R. C., Beral, V. et al. Rethinking ovarian cancer II: reducing mortality from high-grade serous ovarian cancer. *Nat. Rev. Cancer* **15**, 668–679 (2015).

14. Dai, Y., Sun, C., Feng, Y., Jia, Q. & Zhu, B. Potent immunogenicity in BRCA1-mutated patients with high-grade serous ovarian carcinoma. *J. Cell Mol. Med.* https://doi.org/10.1111/jcmm.13678 (2018).

15. McAlpine, J. N., Porter, H., Köbel, M., Nelson, B. H., Prentice, L. M., Kalloger, S. E. et al. BRCA1 and BRCA2 mutations correlate with TP53 abnormalities and presence of immune cell infiltrates in ovarian high-grade serous carcinoma. *Mod. Pathol.* **25**, 740–750 (2012).

16. Zhou, X., Hao, Q. & Lu, H. Mutant p53 in cancer therapy—the barrier or the path. *J. Mol. Cell. Biol.* **11**, 293–305 (2019).

17. Kumar, S., Mohan, A. & Guleria, R. Prognostic implications of circulating anti-p53 antibodies in lung cancer—a review. *Eur. J. Cancer Care (Engl.)* **18**, 248–254 (2009).

18. Peng, W., Chen, J. Q., Liu, C., Malu, S., Creasy, C., Tetzlaff, M. T. et al. Loss of PTEN promotes resistance to T cell-mediated immunotherapy. *Cancer Discov.* **6**, 202–216 (2016).

19. Song, H., Ramus, S. J., Quaye, L., DiCioccio, R. A., Tyrer, J., Lomas, E. et al. Common variants in mismatch repair genes and risk of invasive ovarian cancer. *Carcinogenesis* **27**, 2235–2242 (2006).

20. Walters, S., Maringe, C., Butler, J., Brierley, J. D., Rachet, B. & Coleman, M. P. Comparability of stage data in cancer registries in six countries: lessons from the International Cancer Benchmarking Partnership. *Int. J. Cancer* **132**, 676–685 (2013).

21. Martins, F. C., Santiago, I., Trinh, A., Xian, J., Guo, A., Sayal, K. et al. Combined image and genomic analysis of high-grade serous ovarian cancer reveals PTEN loss as a common driver event and prognostic classifier. *Genome Biol.* **15**, 526 (2014).

22. Piskorz, A. M., Ennis, D., Macintyre, G., Goranova, T. E., Eldridge, M., Segui-Gracia, N. et al. Methanol-based fixation is superior to buffered formalin for next-generation sequencing of DNA from clinical cancer samples. *Ann. Oncol.* **27**, 532–539 (2016).

23. Köbel, M., Piskorz, A. M., Lee, S., Lui, S., LePage, C., Marass, F. et al. Optimized p53 immunohistochemistry is an accurate predictor of. *J. Pathol. Clin. Res.* **2**, 247–258 (2016).

24. Song, H., Cicek, M. S., Dicks, E., Harrington, P., Ramus, S. J., Cunningham, J. M. et al. The contribution of deleterious germline mutations in BRCA1, BRCA2 and the mismatch repair genes to ovarian cancer in the population. *Hum. Mol. Genet.* **23**, 4703–4709 (2014).

25. Ros, X. R. & Vermeulen, L. Turning cold tumors hot by blocking TGF-β. *Trends Cancer* **4**, 335–337 (2018).

26. Haanen, J. B. A. G. Converting cold into hot tumors by combining immunotherapies. *Cell* **170**, 1055–1056 (2017).

27. Chen, D. S. & Mellman, I. Elements of cancer immunity and the cancer-immune set point. *Nature* **541**, 321–330 (2017).

28. Peranzoni, E., Lemoine, J., Vimeux, L., Feuillet, V., Barrin, S., Kantari-Mimoun, C. et al. Macrophages impede CD8 T cells from reaching tumor cells and limit the efficacy of anti-PD-1 treatment. *Proc. Natl Acad. Sci. USA* **115**, E4041–E4050 (2018).

29. Joyce, J. A. & Fearon, D. T. T cell exclusion, immune privilege, and the tumor microenvironment. *Science* **348**, 74–80 (2015).

30. Sato, E., Olson, S. H., Ahn, J., Bundy, B., Nishikawa, H., Qian, F. et al. Intraepithelial CD8⁺ tumor-infiltrating lymphocytes and a high CD8+/regulatory T cell ratio are associated with favorable prognosis in ovarian cancer. *Proc. Natl Acad. Sci. USA* **102**, 18538–18543 (2005).

31. Montfort, A., Pearce, O., Maniati, E., Vincent, B. G., Bixby, L., Böhm, S. et al. A strong B-cell response is part of the immune landscape in human high-grade serous ovarian metastases. *Clin. Cancer Res.* **23**, 250–262 (2017).

32. Bindea, G., Mlecnik, B., Tosolini, M., Kirilovsky, A., Waldner, M., Obenauf, A. C. et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. *Immunity* **39**, 782–795 (2013).

33. Hu, G. & Wang, S. Tumor-infiltrating CD45RO. *Sci. Rep.* **7**, 10376 (2017).

34. Zhang, Q. W., Liu, L., Gong, C. Y., Shi, H. S., Zeng, Y. H., Wang, X. Z. et al. Prognostic significance of tumor-associated macrophages in solid tumor: a meta-analysis of the literature. *PLoS ONE* **7**, e50946 (2012).

35. Malesci, A., Bianchi, P., Celesti, G., Basso, G., Marchesi, F., Grizzi, F. et al. Tumor-associated macrophages and response to 5-fluorouracil adjuvant therapy in stage III colorectal cancer. *Oncoimmunology* **6**, e1342918 (2017).

36. Pinto, M. L., Rios, E., Duraes, C., Ribeiro, R., Machado, J. C., Mantovani, A. et al. The two faces of tumor-associated macrophages and their clinical significance in colorectal cancer. *Front Immunol.* **10**, 1875 (2019).

37. Yang, M., McKay, D., Pollard, J. W. & Lewis, C. E. Diverse functions of macrophages in different tumor microenvironments. *Cancer Res.* **78**, 5492–5503 (2018).

38. Bohm, S., Montfort, A., Pearce, O. M., Topping, J., Chakravarty, P., Everitt, G. L. et al. Neoadjuvant chemotherapy modulates the immune microenvironment in metastases of tubo-ovarian high-grade serous carcinoma. *Clin. Cancer Res.* **22**, 3025–3036 (2016).

39. Mansuet-Lupo, A., Alifano, M., Pécuchet, N., Biton, J., Becht, E., Goc, J. et al. Intratumoral immune cell densities are associated with lung adenocarcinoma gene alterations. *Am. J. Respir. Crit. Care Med.* **194**, 1403–1412 (2016).

40. Ham, S. W., Jeon, H. Y., Jin, X., Kim, E. J., Kim, J. K., Shin, Y. J. et al. TP53 gain-of-function mutation promotes inflammation in glioblastoma. *Cell Death Differ.* **26**, 409–425 (2019).

41. Biton, J., Mansuet-Lupo, A., Pécuchet, N., Alifano, M., Ouakrim, H., Arrondeau, J. et al. TP53, STK11, and EGFR mutations predict tumor immune profile and the response to anti-PD-1 in lung adenocarcinoma. *Clin. Cancer Res.* **24**, 5710–5723 (2018).

42. Dong, Y., Richards, J. A., Gupta, R., Aung, P. P., Emley, A., Kluger, Y. et al. PTEN functions as a melanoma tumor suppressor by promoting host immune response. *Oncogene* **33**, 4632–4642 (2014).

43. Clarke, B., Tinker, A. V., Lee, C. H., Subramanian, S., van de Rijn, M., Turbin, D. et al. Intraepithelial T cells and prognosis in ovarian carcinoma: novel associations with stage, tumor type, and BRCA1 loss. *Mod. Pathol.* **22**, 393–402 (2009).

<figure><img src="image_5.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.