Opening the Black Box: Spatial Transcriptomics and the Relevance of Artificial Intelligence-Detected Prognostic Regions in High-Grade Serous Carcinoma.

Image-based deep learning models are used to extract new information from standard hematoxylin and eosin pathology slides; however, biological interpretation of the features detected by artificial intelligence (AI) remains a challenge. High-grade serous carcinoma of the ovary (HGSC) is characterized by aggressive behavior and chemotherapy resistance, but also exhibits striking variability in outcome. Our understanding of this disease is limited, partly due to considerable tumor heterogeneity. We previously trained an AI model to identify HGSC tumor regions that are highly associated with outcome status but are indistinguishable by conventional morphologic methods. Here, we applied spatially resolved transcriptomics to further profile the AI-identified tumor regions in 16 patients (8 per outcome group) and identify molecular features related to disease outcome in patients who underwent primary debulking surgery and platinum-based chemotherapy. We examined formalin-fixed paraffin-embedded tissue from (1) regions identified by the AI model as highly associated with short or extended chemotherapy response, and (2) background tumor regions (not identified by the AI model as highly associated with outcome status) from the same tumors. We show that the transcriptomic profiles of AI-identified regions are more distinct than background regions from the same tumors, are superior in predicting outcome, and differ in several pathways including those associated with chemoresistance in HGSC. Further, we find that poor outcome and good outcome regions are enriched by different tumor subpopulations, suggesting distinctive interaction patterns. In summary, our work presents proof of concept that AI-guided spatial transcriptomic analysis improves recognition of biologic features relevant to patient outcomes.

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.

High-Resolution Genotyping of Formalin-Fixed Tissue Accurately Estimates Polygenic Risk Scores in Human Diseases.

Formalin-fixed paraffin-embedded (FFPE) tissues stored in biobanks and pathology archives are a vast but underutilized source for molecular studies on different diseases. Beyond being the "gold standard" for preservation of diagnostic human tissues, FFPE samples retain similar genetic information as matching blood samples, which could make FFPE samples an ideal resource for genomic analysis. However, research on this resource has been hindered by the perception that DNA extracted from FFPE samples is of poor quality. Here, we show that germline disease-predisposing variants and polygenic risk scores (PRS) can be identified from FFPE normal tissue (FFPE-NT) DNA with high accuracy. We optimized the performance of FFPE-NT DNA on a genome-wide array containing 657,675 variants. Via a series of testing and validation phases, we established a protocol for FFPE-NT genotyping with results comparable with blood genotyping. The median call rate of FFPE-NT samples in the validation phase was 99.85% (range 98.26%-99.94%) and median concordance with matching blood samples was 99.79% (range 98.85%-99.9%). We also demonstrated that a rare pathogenic PALB2 genetic variant predisposing to cancer can be correctly identified in FFPE-NT samples. We further imputed the FFPE-NT genotype data and calculated the FFPE-NT genome-wide PRS in 3 diseases and 4 disease risk variables. In all cases, FFPE-NT and matching blood PRS were highly concordant (all Pearson's r > 0.95). The ability to precisely genotype FFPE-NT on a genome-wide array enables translational genomics applications of archived FFPE-NT samples with the possibility to link to corresponding phenotypes and longitudinal health data.

DLL4/Notch3/WNT5B axis mediates bidirectional prometastatic crosstalk between melanoma and lymphatic endothelial cells.

Despite strong indications that interactions between melanoma and lymphatic vessels actively promote melanoma progression, the molecular mechanisms are not yet completely understood. To characterize molecular factors of this crosstalk, we established human primary lymphatic endothelial cell (LEC) cocultures with human melanoma cell lines. Here, we show that coculture with melanoma cells induced transcriptomic changes in LECs and led to multiple changes in their function. WNT5B, a paracrine signaling molecule upregulated in melanoma cells upon LEC interaction, was found to contribute to the functional changes in LECs. Moreover, WNT5B transcription was regulated by Notch3 in melanoma cells following the coculture with LECs, and Notch3 and WNT5B were coexpressed in melanoma patient primary tumor and metastasis samples. Moreover, melanoma cells derived from LEC coculture escaped efficiently from the primary site to the proximal tumor-draining lymph nodes, which was impaired upon WNT5B depletion. This supported the role of WNT5B in promoting the metastatic potential of melanoma cells through its effects on LECs. Finally, DLL4, a Notch ligand expressed in LECs, was identified as an upstream inducer of the Notch3/WNT5B axis in melanoma. This study elucidated WNT5B as a key molecular factor mediating bidirectional crosstalk between melanoma cells and lymphatic endothelium and promoting melanoma metastasis.

Characterization of Expression and Function of the Formins FHOD1, INF2, and DAAM1 in HER2-Positive Breast Cancer.

PURPOSE: Human epidermal growth factor receptor 2 (HER2)-targeted therapies, such as trastuzumab, benefit patients with HER2-positive metastatic breast cancer; however, owing to traditional pathway activation or alternative signaling, resistance persists. Given the crucial role of the formin family in shaping the actin cytoskeleton during cancer progression, these proteins may function downstream of the HER2 signaling pathway. Our aim was to uncover the potential correlations between formins and HER2 expression using a combination of public databases, immunohistochemistry, and functional in vitro assays. METHODS: Using online databases, we identified a negative prognostic correlation between specific formins mRNA expression in HER2-positive cancers. To validate these findings at the protein level, immunohistochemistry was performed on HER2 subtype breast cancer tumors to establish the links between staining patterns and clinical characteristics. We then knocked down individual or combined formins in MDA-MB-453 and SK-BR-3 cells and investigated their effects on wound healing, transwell migration, and proliferation. Furthermore, we investigated the effects of erb-b2 receptor tyrosine kinase 2 (ERBB2)/HER2 small interfering RNA (siRNA)-mediated knockdown on the PI3K/Akt and MEK/ERK1 pathways as well as on selected formins. RESULTS: Our results revealed that correlations between INF2, FHOD1, and DAAM1 mRNA expression and ERBB2 in HER2-subtype breast cancer were associated with worse outcomes. Using immunohistochemistry, we found that high FHOD1 protein expression was linked to higher histological grades and was negatively correlated with estrogen and progesterone receptor positivity. Upon formins knockdown, we observed effects on wound healing and transwell migration, with a minimal impact on proliferation, which was evident through single and combined knockdowns in both cell lines. Notably, siRNA-mediated knockdown of HER2 affected FHOD1 and INF2 expression, along with the phosphorylated Akt/MAPK states. CONCLUSION: Our study highlights the roles of FHOD1 and INF2 as downstream effectors of the HER2/Akt and HER2/MAPK pathways, suggesting that they are potential therapeutic targets in HER2-positive breast cancer.

The impact of implementing current treatment modalities and female sex on gastric cancer outcomes, 2000-2016: a longitudinal nationwide cohort study.

BACKGROUND: The implementation of current treatment modalities and their impact on nationwide gastric cancer outcomes remain poorly understood. Biological differences between females and males could impact survival. We aimed to analyze rates of gastric surgery, chemotherapy, and radiotherapy as well as changes in overall survival among gastric cancer patients diagnosed between 2000-2008 and 2009-2016, respectively, in Finland. MATERIAL AND METHODS: Data on gastric cancer patients were collected from national registries. Cox regression analysis and the Kaplan-Meier method were used to analyze differences in survival. RESULTS: We identified 9223 histologically confirmed gastric cancer patients. The rate of gastric surgery decreased from 44% (n = 2282) to 34% (n = 1368; p < 0.001). The proportion of gastric surgery patients who underwent preoperative oncological treatment increased from 0.5% (n = 12) to 16.2% (n = 222) between the calendar periods (p < 0.001) and stood at 30% in 2016. The median overall survival (OS) improved from 30 months [95% confidence interval (CI) 28-33] to 38 months (95%CI 33-42; p = 0.006) and the period 2009-2016 independently associated with a lower risk of death [hazard ratio (HR) 0.78, 95%CI 0.70-0.87] among patients who underwent gastric surgery. Females exhibited a lower risk of death (HR 0.88, 95%CI 0.81-0.97) among patients who underwent gastric surgery. CONCLUSION: Preoperative oncological treatment was gradually introduced into clinical practice and OS among gastric surgery patients improved. Moreover, female surgical patients exhibited a better survival than male patients.

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.

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.

FinnGen provides genetic insights from a well-phenotyped isolated population.

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

Serum FGF19 predicts outcomes of Kasai portoenterostomy in biliary atresia.

BACKGROUND AND AIMS: Outcomes after Kasai portoenterostomy (KPE) for biliary atresia remain highly variable for unclear reasons. As reliable early biomarkers predicting KPE outcomes are lacking, we studied the prognostic value of FGF19. APPROACH AND RESULTS: Serum and liver specimens, obtained from biliary atresia patients (N=87) at KPE or age-matched cholestatic controls (N=26) were included. Serum concentration of FGF19 and bile acids, liver mRNA expression of FGF19 , and key regulators of bile acid synthesis were related to KPE outcomes and liver histopathology. Immunohistochemistry and in situ hybridization were used for the localization of liver FGF19 expression. Serum levels (223 vs. 61 pg/mL, p <0.001) and liver mRNA expression of FGF19 were significantly increased in biliary atresia. Patients with unsuccessful KPE (419 vs. 145 pg/mL, p =0.047), and those subsequently underwent liver transplantation (410 vs. 99 pg/mL, p =0.007) had significantly increased serum, but not liver, FGF19, which localized mainly in hepatocytes. In Cox hazard modeling serum FGF19 <109 pg/mL predicted native liver survival (HR: 4.31, p <0.001) also among patients operated <60 days of age (HR: 8.77, p =0.004) or after successful KPE (HR: 6.76, p =0.01). Serum FGF19 correlated positively with increased serum primary bile acids ( R =0.41, p =0.004) and ductular reaction ( R =0.39, p =0.004). CONCLUSIONS: Increased serum FGF19 at KPE predicted inferior long-term native liver survival in biliary atresia and was associated with unsuccessful KPE, elevated serum primary bile acids, and ductular reaction.

Data-driven comorbidity analysis of 100 common disorders reveals patient subgroups with differing mortality risks and laboratory correlates.

The populational heterogeneity of a disease, in part due to comorbidity, poses several complexities. Individual comorbidity profiles, on the other hand, contain useful information to refine phenotyping, prognostication, and risk assessment, and they provide clues to underlying biology. Nevertheless, the spectrum and the implications of the diagnosis profiles remain largely uncharted. Here we mapped comorbidity patterns in 100 common diseases using 4-year retrospective data from 526,779 patients and developed an online tool to visualize the results. Our analysis exposed disease-specific patient subgroups with distinctive diagnosis patterns, survival functions, and laboratory correlates. Computational modeling and real-world data shed light on the structure, variation, and relevance of populational comorbidity patterns, paving the way for improved diagnostics, risk assessment, and individualization of care. Variation in outcomes and biological correlates of a disease emphasizes the importance of evaluating the generalizability of current treatment strategies, as well as considering the limitations that selective inclusion criteria pose on clinical trials.

Effects of Wee1 inhibitor adavosertib on patient-derived high-grade serous ovarian cancer cells are multiple and independent of homologous recombination status.

Objective: A major challenge in the treatment of platinum-resistant high-grade serous ovarian cancer (HGSOC) is lack of effective therapies. Much of ongoing research on drug candidates relies on HGSOC cell lines that are poorly documented. The goal of this study was to screen for effective, state-of-the-art drug candidates using primary HGSOC cells. In addition, our aim was to dissect the inhibitory activities of Wee1 inhibitor adavosertib on primary and conventional HGSOC cell lines. Methods: A comprehensive drug sensitivity and resistance testing (DSRT) on 306 drug compounds was performed on three patient-derived genetically unique HGSOC cell lines and two commonly used ovarian cancer cell lines. The effect of adavosertib on the cell lines was tested in several assays, including cell-cycle analysis, apoptosis induction, proliferation, wound healing, DNA damage, and effect on nuclear integrity. Results: Several compounds exerted cytotoxic activity toward all cell lines, when tested in both adherent and spheroid conditions. In further cytotoxicity tests, adavosertib exerted the most consistent cytotoxic activity. Adavosertib affected cell-cycle control in patient-derived and conventional HGSOC cells, inducing G2/M accumulation and reducing cyclin B1 levels. It induced apoptosis and inhibited proliferation and migration in all cell lines. Furthermore, the DNA damage marker γH2AX and the number of abnormal cell nuclei were clearly increased following adavosertib treatment. Based on the homologous recombination (HR) signature and functional HR assays of the cell lines, the effects of adavosertib were independent of the cells' HR status. Conclusion: Our study indicates that Wee1 inhibitor adavosertib affects several critical functions related to proliferation, cell cycle and division, apoptosis, and invasion. Importantly, the effects are consistent in all tested cell lines, including primary HGSOC cells, and independent of the HR status of the cells. Wee1 inhibition may thus provide treatment opportunities especially for patients, whose cancer has acquired resistance to platinum-based chemotherapy or PARP inhibitors.

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.

Screening for potential undiagnosed Gaucher disease patients: Utilisation of the Gaucher earlier diagnosis consensus point-scoring system (GED-C PSS) in conjunction with electronic health record data, tissue specimens, and small nucleotide polymorphism (SNP) genotype data available in Finnish biobanks.

Background: Autosomal recessive Gaucher disease (GD) is likely underdiagnosed in many countries. Because the number of diagnosed GD patients in Finland is relatively low, and the true prevalence is currently not known, it was hypothesized that undiagnosed GD patients may exist in Finland. Our previous study demonstrated the applicability of Gaucher Earlier Diagnosis Consensus point-scoring system (GED-C PSS; Mehta et al., 2019) and Finnish biobank data and specimens in the automated point scoring of large populations. An indicative point-score range for Finnish GD patients was determined, but undiagnosed patients were not identified partly due to high number of high-score subjects in combination with a lack of suitable samples for diagnostics in the assessed biobank population. The current study extended the screening to another biobank and evaluated the feasibility of utilising the automated GED-C PSS in conjunction with small nucleotide polymorphism (SNP) chip genotype data from the FinnGen study of biobank sample donors in the identification of undiagnosed GD patients in Finland. Furthermore, the applicability of FFPE tissues and DNA restoration in the next-generation sequencing (NGS) of the GBA gene were tested. Methods: Previously diagnosed Finnish GD patients eligible to the study, and up to 45,100 sample donors in Helsinki Biobank (HBB) were point scored. The GED-C point scoring, adjusted to local data, was automated, but also partly manually verified for GD patients. The SNP chip genotype data for rare GBA variants was visually assessed. FFPE tissues of GD patients were obtained from HBB and Biobank Borealis of Northern Finland (BB). Results: Three previously diagnosed GD patients and one patient previously treated for GD-related features were included. A genetic diagnosis was confirmed for the patient treated for GD-related features. The GED-C point score of the GD patients was 12.5-22.5 in the current study. The score in eight Finnish GD patients of the previous and the current study is thus 6-22.5 points per patient. In the automated point scoring of the HBB subpopulation (N ≈ 45,100), the overall scores ranged from 0 to 17.5, with 0.77% (346/45,100) of the subjects having ≥10 points. The analysis of SNP chip genotype data was able to identify the diagnosed GD patients, but potential undiagnosed patients with the GED-C score and/or the GBA genotype indicative of GD were not discovered. Restoration of the FFPE tissue DNA improved the quality of the GBA NGS, and pathogenic GBA variants were confirmed in five out of six unrestored and in all four restored FFPE DNA samples. Discussion: These findings imply that the prevalence of diagnosed patients (~1:325,000) may indeed correspond the true prevalence of GD in Finland. The SNP chip genotype data is a valuable tool that complements the screening with the GED-C PSS, especially if the genotyping pipeline is tuned for rare variants. These proof-of-concept biobank tools can be adapted to other rare genetic diseases.

Single-Cell RNA Sequencing-Based Characterization of Resident Lung Mesenchymal Stromal Cells in Bronchopulmonary Dysplasia.

Late lung development is a period of alveolar and microvascular formation, which is pivotal in ensuring sufficient and effective gas exchange. Defects in late lung development manifest in premature infants as a chronic lung disease named bronchopulmonary dysplasia (BPD). Numerous studies demonstrated the therapeutic properties of exogenous bone marrow and umbilical cord-derived mesenchymal stromal cells (MSCs) in experimental BPD. However, very little is known regarding the regenerative capacity of resident lung MSCs (L-MSCs) during normal development and in BPD. In this study we aimed to characterize the L-MSC population in homeostasis and upon injury. We used single-cell RNA sequencing (scRNA-seq) to profile in situ Ly6a+ L-MSCs in the lungs of normal and O2-exposed neonatal mice (a well-established model to mimic BPD) at 3 developmental timepoints (postnatal days 3, 7, and 14). Hyperoxia exposure increased the number and altered the expression profile of L-MSCs, particularly by increasing the expression of multiple pro-inflammatory, pro-fibrotic, and anti-angiogenic genes. In order to identify potential changes induced in the L-MSCs transcriptome by storage and culture, we profiled 15 000 Ly6a+ L-MSCs after in vitro culture. We observed great differences in expression profiles of in situ and cultured L-MSCs, particularly those derived from healthy lungs. Additionally, we have identified the location of Ly6a+/Col14a1+ L-MSCs in the developing lung and propose Serpinf1 as a novel, culture-stable marker of L-MSCs. Finally, cell communication analysis suggests inflammatory signals from immune and endothelial cells as main drivers of hyperoxia-induced changes in L-MSCs transcriptome.

Ovarian Cancers with Low CIP2A Tumor Expression Constitute an APR-246-Sensitive Disease Subtype.

Identification of ovarian cancer patient subpopulations with increased sensitivity to targeted therapies could offer significant clinical benefit. We report that 22% of the high-grade ovarian cancer tumors at diagnosis express CIP2A oncoprotein at low levels. Furthermore, regardless of their significantly lower likelihood of disease relapse after standard chemotherapy, a portion of relapsed tumors retain their CIP2A-deficient phenotype. Through a screen for therapeutics that would preferentially kill CIP2A-deficient ovarian cancer cells, we identified reactive oxygen species inducer APR-246, tested previously in ovarian cancer clinical trials. Consistent with CIP2A-deficient ovarian cancer subtype in humans, CIP2A is dispensable for development of MISIIR-Tag-driven mouse ovarian cancer tumors. Nevertheless, CIP2A-null ovarian cancer tumor cells from MISIIR-Tag mice displayed APR-246 hypersensitivity both in vitro and in vivo. Mechanistically, the lack of CIP2A expression hypersensitizes the ovarian cancer cells to APR-246 by inhibition of NF-κB activity. Accordingly, combination of APR-246 and NF-κB inhibitor compounds strongly synergized in killing of CIP2A-positive ovarian cancer cells. Collectively, the results warrant consideration of clinical testing of APR-246 for CIP2A-deficient ovarian cancer tumor subtype patients. Results also reveal CIP2A as a candidate APR-246 combination therapy target for ovarian cancer.

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.