Chemotherapy induces myeloid-driven spatial T-cell exhaustion in ovarian cancer.

To uncover the intricate, chemotherapy-induced spatiotemporal remodeling of the tumor microenvironment, we conducted integrative spatial and molecular characterization of 97 high-grade serous ovarian cancer (HGSC) samples collected before and after chemotherapy. Using single-cell and spatial analyses, we identify increasingly versatile immune cell states, which form spatiotemporally dynamic microcommunities at the tumor-stroma interface. We demonstrate that chemotherapy triggers spatial redistribution and exhaustion of CD8+ T cells due to prolonged antigen presentation by macrophages, both within interconnected myeloid networks termed "Myelonets" and at the tumor stroma interface. Single-cell and spatial transcriptomics identifies prominent TIGIT-NECTIN2 ligand-receptor interactions induced by chemotherapy. Using a functional patient-derived immuno-oncology platform, we show that CD8+T-cell activity can be boosted by combining immune checkpoint blockade with chemotherapy. Our discovery of chemotherapy-induced myeloid-driven spatial T-cell exhaustion paves the way for novel immunotherapeutic strategies to unleash CD8+ T-cell-mediated anti-tumor immunity in HGSC.

PIK3R1 fusion drives chemoresistance in ovarian cancer by activating ERK1/2 and inducing rod and ring-like structures.

Gene fusions are common in high-grade serous ovarian cancer (HGSC). Such genetic lesions may promote tumorigenesis, but the pathogenic mechanisms are currently poorly understood. Here, we investigated the role of a PIK3R1-CCDC178 fusion identified from a patient with advanced HGSC. We show that the fusion induces HGSC cell migration by regulating ERK1/2 and increases resistance to platinum treatment. Platinum resistance was associated with rod and ring-like cellular structure formation. These structures contained, in addition to the fusion protein, CIN85, a key regulator of PI3K-AKT-mTOR signaling. Our data suggest that the fusion-driven structure formation induces a previously unrecognized cell survival and resistance mechanism, which depends on ERK1/2-activation.

Genome-wide quantification of copy-number aberration impact on gene expression in ovarian high-grade serous carcinoma.

Copy-number alterations (CNAs) are a hallmark of cancer and can regulate cancer cell states via altered gene expression values. Herein, we have developed a copy-number impact (CNI) analysis method that quantifies the degree to which a gene expression value is impacted by CNAs and leveraged this analysis at the pathway level. Our results show that a high CNA is not necessarily reflected at the gene expression level, and our method is capable of detecting genes and pathways whose activity is strongly influenced by CNAs. Furthermore, the CNI analysis enables unbiased categorization of CNA categories, such as deletions and amplifications. We identified six CNI-driven pathways associated with poor treatment response in ovarian high-grade serous carcinoma (HGSC), which we found to be the most CNA-driven cancer across 14 cancer types. The key driver in most of these pathways was amplified wild-type KRAS, which we validated functionally using CRISPR modulation. Our results suggest that wild-type KRAS amplification is a driver of chemotherapy resistance in HGSC and may serve as a potential treatment target.

HRD related signature 3 predicts clinical outcome in advanced tubo-ovarian high-grade serous carcinoma.

OBJECTIVES: We evaluated usability of single base substitution signature 3 (Sig3) as a biomarker for homologous recombination deficiency (HRD) in tubo-ovarian high-grade serous carcinoma (HGSC). MATERIALS AND METHODS: This prospective observational trial includes 165 patients with advanced HGSC. Fresh tissue samples (n = 456) from multiple intra-abdominal areas at diagnosis and after neoadjuvant chemotherapy (NACT) were collected for whole-genome sequencing. Sig3 was assessed by fitting samples independently with COSMIC v3.2 reference signatures. An HR scar assay was applied for comparison. Progression-free survival (PFS) and overall survival (OS) were studied using Kaplan-Meier and Cox regression analysis. RESULTS: Sig3 has a bimodal distribution, eliminating the need for an arbitrary cutoff typical in HR scar tests. Sig3 could be assessed from samples with low (10%) cancer cell proportion and was consistent between multiple samples and stable during NACT. At diagnosis, 74 (45%) patients were HRD (Sig3+), while 91 (55%) were HR proficient (HRP, Sig3-). Sig3+ patients had longer PFS and OS than Sig3- patients (22 vs. 13 months and 51 vs. 34 months respectively, both p < 0.001). Sig3 successfully distinguished the poor prognostic HRP group among BRCAwt patients (PFS 19 months for Sig3+ and 13 months for Sig3- patients, p < 0.001). However, Sig3 at diagnosis did not predict chemoresponse anymore in the first relapse. The patient-level concordance between Sig3 and HR scar assay was 87%, and patients with HRD according to both tests had the longest median PFS. CONCLUSIONS: Sig3 is a prognostic marker in advanced HGSC and useful tool in patient stratification for HRD.

Quantification of cell death and proliferation of patient-derived ovarian cancer organoids through 3D imaging and image analysis.

Patient-derived organoids (PDOs) are ideal ex vivo model systems to study cancer progression and drug resistance mechanisms. Here, we present a protocol for measuring drug efficacy in three-dimensional (3D) high-grade serous ovarian cancer PDO cultures through quantification of cytotoxicity using propidium iodide incorporation in dead cells. We also provide detailed steps to analyze proliferation of PDOs using the Ki67 biomarker. We describe steps for sample processing, immunofluorescent staining, high-throughput confocal imaging, and image-based quantification for 3D cultures. For complete details on the use and execution of this protocol, please refer to Lahtinen et al. (2023).1.

H&E image analysis pipeline for quantifying morphological features.

Detecting cell types from histopathological images is essential for various digital pathology applications. However, large number of cells in whole-slide images (WSIs) necessitates automated analysis pipelines for efficient cell type detection. Herein, we present hematoxylin and eosin (H&E) Image Processing pipeline (HEIP) for automatied analysis of scanned H&E-stained slides. HEIP is a flexible and modular open-source software that performs preprocessing, instance segmentation, and nuclei feature extraction. To evaluate the performance of HEIP, we applied it to extract cell types from ovarian high-grade serous carcinoma (HGSC) patient WSIs. HEIP showed high precision in instance segmentation, particularly for neoplastic and epithelial cells. We also show that there is a significant correlation between genomic ploidy values and morphological features, such as major axis of the nucleus.

Circulating tumor DNA-based copy-number profiles enable monitoring treatment effects during therapy in high-grade serous carcinoma.

Circulating tumor DNA (ctDNA) analysis has emerged as a promising tool for detecting and profiling longitudinal genomics changes in cancer. While copy-number alterations (CNAs) play a major role in cancers, treatment effect monitoring using copy-number profiles has received limited attention as compared to mutations. A major reason for this is the insensitivity of CNA analysis for the real-life tumor-fraction ctDNA samples. We performed copy-number analysis on 152 plasma samples obtained from 29 patients with high-grade serous ovarian cancer (HGSC) using a sequencing panel targeting over 500 genes. Twenty-one patients had temporally matched tissue and plasma sample pairs, which enabled assessing concordance with tissues sequenced with the same panel or whole-genome sequencing and to evaluate sensitivity. Our approach could detect concordant CNA profiles in most plasma samples with as low as 5% tumor content and highly amplified regions in samples with ∼1% of tumor content. Longitudinal profiles showed changes in the CNA profiles in seven out of 11 patients with high tumor-content plasma samples at relapse. These changes included focal acquired or lost copy-numbers, even though most of the genome remained stable. Two patients displayed major copy-number profile changes during therapy. Our analysis revealed ctDNA-detectable subclonal selection resulting from both surgical operations and chemotherapy. Overall, longitudinal ctDNA data showed acquired and diminished CNAs at relapse when compared to pre-treatment samples. These results highlight the importance of genomic profiling during treatment as well as underline the usability of ctDNA.

Single-cell epigenome analysis identifies molecular events controlling direct conversion of human fibroblasts to pancreatic ductal-like cells.

Cell fate can be reprogrammed by ectopic expression of lineage-specific transcription factors (TFs). However, the exact cell state transitions during transdifferentiation are still poorly understood. Here, we have generated pancreatic exocrine cells of ductal epithelial identity from human fibroblasts using a set of six TFs. We mapped the molecular determinants of lineage dynamics using a factor-indexing method based on single-nuclei multiome sequencing (FI-snMultiome-seq) that enables dissecting the role of each individual TF and pool of TFs in cell fate conversion. We show that transition from mesenchymal fibroblast identity to epithelial pancreatic exocrine fate involves two deterministic steps: an endodermal progenitor state defined by activation of HHEX with FOXA2 and SOX17 and a temporal GATA4 activation essential for the maintenance of pancreatic cell fate program. Collectively, our data suggest that transdifferentiation-although being considered a direct cell fate conversion method-occurs through transient progenitor states orchestrated by stepwise activation of distinct TFs.

ROR1-STAT3 signaling contributes to ovarian cancer intra-tumor heterogeneity.

Wnt pathway dysregulation through genetic and non-genetic alterations occurs in multiple cancers, including ovarian cancer (OC). The aberrant expression of the non-canonical Wnt signaling receptor ROR1 is thought to contribute to OC progression and drug resistance. However, the key molecular events mediated by ROR1 that are involved in OC tumorigenesis are not fully understood. Here, we show that ROR1 expression is enhanced by neoadjuvant chemotherapy, and Wnt5a binding to ROR1 can induce oncogenic signaling via AKT/ERK/STAT3 activation in OC cells. Proteomics analysis of isogenic ROR1-knockdown OC cells identified STAT3 as a downstream effector of ROR1 signaling. Transcriptomics analysis of clinical samples (n = 125) revealed that ROR1 and STAT3 are expressed at higher levels in stromal cells than in epithelial cancer cells of OC tumors, and these findings were corroborated by multiplex immunohistochemistry (mIHC) analysis of an independent OC cohort (n = 11). Our results show that ROR1 and its downstream STAT3 are co-expressed in epithelial as well as stromal cells of OC tumors, including cancer-associated fibroblasts or CAFs. Our data provides the framework to expand the clinical utility of ROR1 as a therapeutic target to overcome OC progression.

A platform for efficient establishment and drug-response profiling of high-grade serous ovarian cancer organoids.

The broad research use of organoids from high-grade serous ovarian cancer (HGSC) has been hampered by low culture success rates and limited availability of fresh tumor material. Here, we describe a method for generation and long-term expansion of HGSC organoids with efficacy markedly improved over previous reports (53% vs. 23%-38%). We established organoids from cryopreserved material, demonstrating the feasibility of using viably biobanked tissue for HGSC organoid derivation. Genomic, histologic, and single-cell transcriptomic analyses revealed that organoids recapitulated genetic and phenotypic features of original tumors. Organoid drug responses correlated with clinical treatment outcomes, although in a culture conditions-dependent manner and only in organoids maintained in human plasma-like medium (HPLM). Organoids from consenting patients are available to the research community through a public biobank and organoid genomic data are explorable through an interactive online tool. Taken together, this resource facilitates the application of HGSC organoids in basic and translational ovarian cancer research.

Evolutionary states and trajectories characterized by distinct pathways stratify patients with ovarian high grade serous carcinoma.

Ovarian high-grade serous carcinoma (HGSC) is typically diagnosed at an advanced stage, with multiple genetically heterogeneous clones existing in the tumors long before therapeutic intervention. Herein we integrate clonal composition and topology using whole-genome sequencing data from 510 samples of 148 patients with HGSC in the prospective, longitudinal, multiregion DECIDER study. Our results reveal three evolutionary states, which have distinct features in genomics, pathways, and morphological phenotypes, and significant association with treatment response. Nested pathway analysis suggests two evolutionary trajectories between the states. Experiments with five tumor organoids and three PI3K inhibitors support targeting tumors with enriched PI3K/AKT pathway with alpelisib. Heterogeneity analysis of samples from multiple anatomical sites shows that site-of-origin samples have 70% more unique clones than metastatic tumors or ascites. In conclusion, these analysis and visualization methods enable integrative tumor evolution analysis to identify patient subtypes using data from longitudinal, multiregion cohorts.

Toxicity and therapy outcome associations in LIG3, SLCO1B3, ABCB1, OPRM1 and GSTP1 in high-grade serous ovarian cancer.

Adverse effects are the major limiting factors in combinatorial chemotherapies. To identify genetic associations in ovarian cancer chemotherapy-induced toxicities and therapy outcomes, we examined a cohort of 101 patients receiving carboplatin-paclitaxel treatment with advanced high-grade serous ovarian cancers. Based on literature and database searches, we selected 19 candidate polymorphisms, designed a multiplex single nucleotide polymorphism-genotyping assay and applied Cox regression analysis, case-control association statistics and the log-rank Mantel-Cox test. In the Cox regression analysis, the SLCO1B3 rs1052536 AA-genotype was associated with a reduced risk of any severe toxicity (hazard ratio = 0.35, p = 0.023). In chi-square allelic test, the LIG3 rs1052536 T-allele was associated with an increased risk of neuropathy (odds ratio [OR] = 2.79, p = 0.031) and GSTP1 rs1695 G allele with a poorer response in the first-line chemotherapy (OR = 2.65, p = 0.026). In Kaplan-Meier survival analysis, ABCB1 rs2032582 TT-genotype was associated with shorter overall survival (uncorrected p = 0.025) and OPRM1 rs544093 GG and GT genotypes with shorter platinum-free interval (uncorrected p = 0.027) and progression-free survival (uncorrected p = 0.012). Results suggest that SLCO1B3 and LIG3 variants are associated with the risk of adverse effects in patients receiving carboplatin-paclitaxel treatment, the GSTP1 variant may affect the treatment response and ABCB1 and OPRM1 variants may influence the prognosis.

Functional Homologous Recombination Assay on FFPE Specimens of Advanced High-Grade Serous Ovarian Cancer Predicts Clinical Outcomes.

PURPOSE: Deficiency in homologous recombination (HR) repair of DNA damage is characteristic of many high-grade serous ovarian cancers (HGSC). It is imperative to identify patients with homologous recombination-deficient (HRD) tumors as they are most likely to benefit from platinum-based chemotherapy and PARP inhibitors (PARPi). Existing methods measure historical, not necessarily current HRD and/or require high tumor cell content, which is not achievable for many patients. We set out to develop a clinically feasible assay for identifying functionally HRD tumors that can predict clinical outcomes. EXPERIMENTAL DESIGN: We quantified RAD51, a key HR protein, in immunostained formalin-fixed, paraffin-embedded (FFPE) tumor samples obtained from chemotherapy-naïve and neoadjuvant chemotherapy (NACT)-treated HGSC patients. We defined cutoffs for functional HRD separately for these sample types, classified the patients accordingly as HRD or HR-proficient, and analyzed correlations with clinical outcomes. From the same specimens, genomics-based HRD estimates (HR gene mutations, genomic signatures, and genomic scars) were also determined, and compared with functional HR (fHR) status. RESULTS: fHR status significantly predicted several clinical outcomes, including progression-free survival (PFS) and overall survival (OS), when determined from chemo-naïve (PFS, P < 0.0001; OS, P < 0.0001) as well as NACT-treated (PFS, P < 0.0001; OS, P = 0.0033) tumor specimens. The fHR test also identified as HRD those PARPi-at-recurrence-treated patients with longer OS (P = 0.0188). CONCLUSIONS: We developed an fHR assay performed on routine FFPE specimens, obtained from either chemo-naïve or NACT-treated HGSC patients, that can significantly predict real-world platinum-based chemotherapy and PARPi response. See related commentary by Garg and Oza, p. 2957.

Kinase Inhibitors in the Treatment of Ovarian Cancer: Current State and Future Promises.

Ovarian cancer is the deadliest gynecological cancer, the high-grade serous ovarian carcinoma (HGSC) being its most common and most aggressive form. Despite the latest therapeutical advancements following the introduction of vascular endothelial growth factor receptor (VEGFR) targeting angiogenesis inhibitors and poly-ADP-ribose-polymerase (PARP) inhibitors to supplement the standard platinum- and taxane-based chemotherapy, the expected overall survival of HGSC patients has not improved significantly from the five-year rate of 42%. This calls for the development and testing of more efficient treatment options. Many oncogenic kinase-signaling pathways are dysregulated in HGSC. Since small-molecule kinase inhibitors have revolutionized the treatment of many solid cancers due to the generality of the increased activation of protein kinases in carcinomas, it is reasonable to evaluate their potential against HGSC. Here, we present the latest concluded and on-going clinical trials on kinase inhibitors in HGSC, as well as the recent work concerning ovarian cancer patient organoids and xenograft models. We discuss the potential of kinase inhibitors as personalized treatments, which would require comprehensive assessment of the biological mechanisms underlying tumor spread and chemoresistance in individual patients, and their connection to tumor genome and transcriptome to establish identifiable subgroups of patients who are most likely to benefit from a given therapy.

Optimized detection of homologous recombination deficiency improves the prediction of clinical outcomes in cancer.

Homologous recombination DNA-repair deficiency (HRD) is a common driver of genomic instability and confers a therapeutic vulnerability in cancer. The accurate detection of somatic allelic imbalances (AIs) has been limited by methods focused on BRCA1/2 mutations and using mixtures of cancer types. Using pan-cancer data, we revealed distinct patterns of AIs in high-grade serous ovarian cancer (HGSC). We used machine learning and statistics to generate improved criteria to identify HRD in HGSC (ovaHRDscar). ovaHRDscar significantly predicted clinical outcomes in three independent patient cohorts with higher precision than previous methods. Characterization of 98 spatiotemporally distinct metastatic samples revealed low intra-patient variation and indicated the primary tumor as the preferred site for clinical sampling in HGSC. Further, our approach improved the prediction of clinical outcomes in triple-negative breast cancer (tnbcHRDscar), validated in two independent patient cohorts. In conclusion, our tumor-specific, systematic approach has the potential to improve patient selection for HR-targeted therapies.

A Subset of Secreted Proteins in Ascites Can Predict Platinum-Free Interval in Ovarian Cancer.

The time between the last cycle of chemotherapy and recurrence, the platinum-free interval (PFI), predicts overall survival in high-grade serous ovarian cancer (HGSOC). To identify secreted proteins associated with a shorter PFI, we utilized machine learning to predict the PFI from ascites composition. Ascites from stage III/IV HGSOC patients treated with neoadjuvant chemotherapy (NACT) or primary debulking surgery (PDS) were screened for secreted proteins and Lasso regression models were built to predict the PFI. Through regularization techniques, the number of analytes used in each model was reduced; to minimize overfitting, we utilized an analysis of model robustness. This resulted in models with 26 analytes and a root-mean-square error (RMSE) of 19 days for the NACT cohort and 16 analytes and an RMSE of 7 days for the PDS cohort. High concentrations of MMP-2 and EMMPRIN correlated with a shorter PFI in the NACT patients, whereas high concentrations of uPA Urokinase and MMP-3 correlated with a shorter PFI in PDS patients. Our results suggest that the analysis of ascites may be useful for outcome prediction and identified factors in the tumor microenvironment that may lead to worse outcomes. Our approach to tuning for model stability, rather than only model accuracy, may be applicable to other biomarker discovery tasks.

Multiomics characterization implicates PTK7 in ovarian cancer EMT and cell plasticity and offers strategies for therapeutic intervention.

Most patients with ovarian cancer (OC) are diagnosed at a late stage when there are very few therapeutic options and a poor prognosis. This is due to the lack of clearly defined underlying mechanisms or an oncogenic addiction that can be targeted pharmacologically, unlike other types of cancer. Here, we identified protein tyrosine kinase 7 (PTK7) as a potential new therapeutic target in OC following a multiomics approach using genetic and pharmacological interventions. We performed proteomics analyses upon PTK7 knockdown in OC cells and identified novel downstream effectors such as synuclein-γ (SNCG), SALL2, and PP1γ, and these findings were corroborated in ex vivo primary samples using PTK7 monoclonal antibody cofetuzumab. Our phosphoproteomics analyses demonstrated that PTK7 modulates cell adhesion and Rho-GTPase signaling to sustain epithelial-mesenchymal transition (EMT) and cell plasticity, which was confirmed by high-content image analysis of 3D models. Furthermore, using high-throughput drug sensitivity testing (525 drugs) we show that targeting PTK7 exhibited synergistic activity with chemotherapeutic agent paclitaxel, CHK1/2 inhibitor prexasertib, and PLK1 inhibitor GSK461364, among others, in OC cells and ex vivo primary samples. Taken together, our study provides unique insight into the function of PTK7, which helps to define its role in mediating aberrant Wnt signaling in ovarian cancer.

POIBM: batch correction of heterogeneous RNA-seq datasets through latent sample matching.

MOTIVATION: RNA sequencing and other high-throughput technologies are essential in understanding complex diseases, such as cancers, but are susceptible to technical factors manifesting as patterns in the measurements. These batch patterns hinder the discovery of biologically relevant patterns. Unbiased batch effect correction in heterogeneous populations currently requires special experimental designs or phenotypic labels, which are not readily available for patient samples in existing datasets. RESULTS: We present POIBM, an RNA-seq batch correction method, which learns virtual reference samples directly from the data. We use a breast cancer cell line dataset to show that POIBM exceeds or matches the performance of previous methods, while being blind to the phenotypes. Further, we analyze The Cancer Genome Atlas RNA-seq data to show that batch effects plague many cancer types; POIBM effectively discovers the true replicates in stomach adenocarcinoma; and integrating the corrected data in endometrial carcinoma improves cancer subtyping. AVAILABILITY AND IMPLEMENTATION: https://bitbucket.org/anthakki/poibm/ (archived at https://doi.org/10.5281/zenodo.6122436). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Longitudinal single-cell RNA-seq analysis reveals stress-promoted chemoresistance in metastatic ovarian cancer.

Chemotherapy resistance is a critical contributor to cancer mortality and thus an urgent unmet challenge in oncology. To characterize chemotherapy resistance processes in high-grade serous ovarian cancer, we prospectively collected tissue samples before and after chemotherapy and analyzed their transcriptomic profiles at a single-cell resolution. After removing patient-specific signals by a novel analysis approach, PRIMUS, we found a consistent increase in stress-associated cell state during chemotherapy, which was validated by RNA in situ hybridization and bulk RNA sequencing. The stress-associated state exists before chemotherapy, is subclonally enriched during the treatment, and associates with poor progression-free survival. Co-occurrence with an inflammatory cancer-associated fibroblast subtype in tumors implies that chemotherapy is associated with stress response in both cancer cells and stroma, driving a paracrine feed-forward loop. In summary, we have found a resistant state that integrates stromal signaling and subclonal evolution and offers targets to overcome chemotherapy resistance.

QuantISH: RNA in situ hybridization image analysis framework for quantifying cell type-specific target RNA expression and variability.

RNA in situ hybridization (RNA-ISH) is a powerful spatial transcriptomics technology to characterize target RNA abundance and localization in individual cells. This allows analysis of tumor heterogeneity and expression localization, which are not readily obtainable through transcriptomic data analysis. RNA-ISH experiments produce large amounts of data and there is a need for automated analysis methods. Here we present QuantISH, a comprehensive open-source RNA-ISH image analysis pipeline that quantifies marker expressions in individual carcinoma, immune, and stromal cells on chromogenic or fluorescent in situ hybridization images. QuantISH is designed to be modular and can be adapted to various image and sample types and staining protocols. We show that in chromogenic RNA in situ hybridization images of high-grade serous carcinoma (HGSC) QuantISH cancer cell classification has high precision, and signal expression quantification is in line with visual assessment. We further demonstrate the power of QuantISH by showing that CCNE1 average expression and DDIT3 expression variability, as captured by the variability factor developed herein, act as candidate biomarkers in HGSC. Altogether, our results demonstrate that QuantISH can quantify RNA expression levels and their variability in carcinoma cells, and thus paves the way to utilize RNA-ISH technology.