Background High-grade serous tubo-ovarian cancer (HGSTOC) is characterised by extensive inter- and intratumour heterogeneity, resulting in persistent therapeutic resistance and poor disease outcome. Molecular subtype classification based on bulk RNA sequencing facilitates a more accurate characterisation of this heterogeneity, but the lack of strong prognostic or predictive correlations with these subtypes currently hinders their clinical implementation. Stromal admixture profoundly affects the prognostic impact of the molecular subtypes, but the contribution of stromal cells to each subtype has poorly been characterised. Increasing the transcriptomic resolution of the molecular subtypes based on single-cell RNA sequencing (scRNA-seq) may provide insights in the prognostic and predictive relevance of these subtypes. Methods We performed scRNA-seq of 18,403 cells unbiasedly collected from 7 treatment-naive HGSTOC tumours. For each phenotypic cluster of tumour or stromal cells, we identified specific transcriptomic markers. We explored which phenotypic clusters correlated with overall survival based on expression of these transcriptomic markers in microarray data of 1467 tumours. By evaluating molecular subtype signatures in single cells, we assessed to what extent a phenotypic cluster of tumour or stromal cells contributes to each molecular subtype. Results We identified 11 cancer and 32 stromal cell phenotypes in HGSTOC tumours. Of these, the relative frequency of myofibroblasts, TGF-β-driven cancer-associated fibroblasts, mesothelial cells and lymphatic endothelial cells predicted poor outcome, while plasma cells correlated with more favourable outcome. Moreover, we identified a clear cell-like transcriptomic signature in cancer cells, which correlated with worse overall survival in HGSTOC patients. Stromal cell phenotypes differed substantially between molecular subtypes. For instance, the mesenchymal, immunoreactive and differentiated signatures were characterised by specific fibroblast, immune cell and myofibroblast/mesothelial cell phenotypes, respectively. Cell phenotypes correlating with poor outcome were enriched in molecular subtypes associated with poor outcome. Conclusions We used scRNA-seq to identify stromal cell phenotypes predicting overall survival in HGSTOC patients. These stromal features explain the association of the molecular subtypes with outcome but also the latter’s weakness of clinical implementation. Stratifying patients based on marker genes specific for these phenotypes represents a promising approach to predict prognosis or response to therapy.

No abstract

Protein aggregation is an underappreciated mechanism that may contribute to the loss- and oncogenic-gain-of-function of mutant tumor suppressors such as p53 and axin. In the present study, we describe amyloid-like aggregation behaviour of the second most frequently mutated tumor suppressor in human cancer, PTEN. In silico analysis revealed a particularly high aggregation vulnerability for this protein, which was corroborated by in vitro aggregation assays. In cultured tumor cells, we found that under stress conditions, PTEN readily undergoes amyloid-like aggregation as a result of mutation. However, we also show that severe dysregulation of protein homeostasis may lead to aggregation of wild-type PTEN. These observations were supported by a small survey of patient-derived uterine tumor tissues, which found that more than 25% of tumors analyzed displayed wild-type PTEN aggregation. Finally, in an exploratory clinical study we found that PTEN aggregation status was correlated with a decline in clinical outcome. Our findings establish that the tumor suppressor PTEN is highly aggregation-prone and our work suggests that protein aggregation might be an underestimated but prevalent component of cancer cell biology.