Single-cell spatial atlas of high-grade serous ovarian cancer unveils MHC class II as a key driver of spatial tumor ecosystems and clinical outcomes

Fernando Perez-Villatoro,Lilian van Wagensveld,Aleksandra Shabanova,Ada Junquera,Ziqi Kang,Iga Niemiec,Matias M Falco,Ella Anttila,Julia Casado,Eric Marcus,Essi Kahelin,Foteini Chamchougia,Matilda Salko,Saundarya Shah,Salvatore Russo,Jacopo Chiaro,Mikaela Gronholm,Gabe S. Sonke,Koen K Van de Vijver,Rutgerus FPM Kruitwagen,Maaike van der Aa,Anni Virtanen,Vincenzo Cerullo,Anna Vaharautio,Peter Karl Sorger,Hugo M Horlings,Anniina Farkkila

DOI: https://doi.org/10.1101/2024.11.29.626039

2024-12-03

Abstract:The tumor microenvironment (TME) is a complex network of interactions between malignant and host cells, yet its orchestration in advanced high-grade serous ovarian carcinoma (HGSC) remains poorly understood. We present a comprehensive single-cell spatial atlas of 280 metastatic HGSCs, integrating high-dimensional imaging, genomics, and transcriptomics. Using 929 single-cell maps, we identify distinct spatial domains associated with phenotypically heterogeneous cellular compositions, and demonstrate that immune cell co-infiltration at the tumor-stroma interface significantly influences clinical outcomes. To uncover the key drivers of the tumor ecosystem, we developed CEFIIRA (Cell Feature Importance Identification by RAndom forest), which identified tumor cell-intrinsic MHC-II expression as a critical predictor of prolonged survival, independent of clinicomolecular profiles. Validation with external datasets confirmed that MHC-II-expressing cancer cells drive immune infiltration and orchestrate spatial tumor-immune interactions. Our atlas offers novel insights into immune surveillance mechanisms across HGSC clinicomolecular groups, paving the way for improved therapeutic strategies and patient stratification.

Cancer Biology

What problem does this paper attempt to address?

The problems that this paper attempts to solve are: **the complex interactions in the tumor microenvironment (TME) of high - grade serous ovarian cancer (HGSC) and their impact on clinical outcomes, especially the role of immune cell co - infiltration at the tumor - stroma interface (TSI) and the key driving role of MHC class II molecules therein**.

Specifically, the research aims to:

1. **Reveal the spatial distribution of different cell types in the tumor microenvironment and their relationship with clinical prognosis**:
- By constructing a single - cell spatial atlas, the researchers analyzed 280 metastatic HGSC samples, integrating high - dimensional imaging, genomics, and transcriptomics data.
- Different spatial domains related to phenotypically heterogeneous cell composition were identified, and the significant impact of immune cell co - infiltration at the tumor - stroma interface on clinical outcomes was demonstrated.

2. **Determine the key drivers of the tumor ecosystem**:
- The researchers developed a machine - learning tool named CEFIIRA (Cell Feature Importance Identification by RAndom forest) to identify key features related to survival time.
- The results showed that the expression of MHC class II molecules intrinsic to tumor cells is a key predictor that can prolong survival independently of clinical molecular features.

3. **Verify the impact of MHC class II molecule expression on immune infiltration and clinical outcomes**:
- Using external data sets for verification, it was found that MHC class II - positive cancer cells can drive immune infiltration and coordinate tumor - immune interactions.
- The expression of MHC class II molecules is associated with more active immune cell infiltration, including CD4+ and CD8+ T cells, CD15+ myeloid cells, and M1 - like macrophages, thereby improving patient prognosis.

### Formula Examples

To ensure the correctness and readability of formulas, the following are some formula examples involved in the paper:

- **Shannon Diversity Index**: It is used to measure the heterogeneity of cell composition. The formula is:
\[
H = -\sum_{i = 1}^{n} p_i \log(p_i)
\]
where \( p_i \) is the proportion of the \( i \)-th cell type.

- **Cox Regression Model**: It is used to evaluate the impact of different features on survival time. The formula is:
\[
h(t) = h_0(t) \exp(\beta_1 X_1 + \beta_2 X_2+\cdots+\beta_p X_p)
\]
where \( h(t) \) is the hazard function, \( h_0(t) \) is the baseline hazard function, \( \beta_i \) is the regression coefficient, and \( X_i \) is the covariate.

These formulas help researchers quantitatively analyze and interpret data, further supporting the research conclusions.