Abstract
Introduction: High-grade serous ovarian carcinoma (HGSOC) remains a medical challenge despite considerable improvements in the treatment. Unfortunately, over 75% of patients have already metastasized at the time of diagnosis. Advances in understanding the mechanisms underlying how ascites cause chemoresistance are urgently needed to derive novel therapeutic strategies. This study aimed to identify the molecular markers involved in drug sensitivity and highlight the use of ascites as a potential model to investigate HGSOC treatment options.
Methods: After conducting an in silico analysis, eight epithelial–mesenchymal (EM)-associated genes related to chemoresistance were identified. To evaluate differences in EM-associated genes in HGSOC samples, we analyzed ascites-derived HGSOC primary cell culture (AS), tumor (T), and peritoneal nodule (NP) samples. Moreover, in vitro experiments were employed to measure tumor cell proliferation and cell migration in AS, following treatment with doxorubicin (DOX) and cisplatin (CIS) and expression of these markers.
Results: Our results showed that AS exhibits a mesenchymal phenotype compared to tumor and peritoneal nodule samples. Moreover, DOX and CIS treatment leads to an invasive-intermediate epithelial-to-mesenchymal transition (EMT) state of the AS by different EM-associated marker expression. For instance, the treatment of AS showed that CDH1 and GATA6 decreased after CIS exposure and increased after DOX treatment. On the contrary, the expression of KRT18 has an opposite pattern.
Conclusion: Taken together, our study reports a comprehensive investigation of the EM-associated genes after drug exposure of AS. Exploring ascites and their associated cellular and soluble components is promising for understanding the HGSOC progression and treatment response at a personalized level.
1 Introduction
HGSOC is a major cause of gynecological tumor mortality, with approximately 70% out of 200,000 ovarian cancer deaths annually (), according to the GLOBOCAN 2020 report (Sung et al., 2021). Using the most recent data available, the age-standardized mortality rate of ovarian cancer in Romania was approximately 2.9 deaths per 100,000 population (), World Health Organization (2019).
HGSOC patients are typically asymptomatic in the early stages and are treated with primary debulking surgery (PDS), while chemotherapy is the preferred option for advanced-stage cancer. Carboplatin/cisplatin (CARB/CIS) and paclitaxel (PAX) are first-line chemotherapeutic agents, but the US Food and Drug Administration (FDA) has also approved DOX, gemcitabine (GEM), irinotecan (MM-398), etoposide (ETOP), oxaliplatin (OHP), and 5-fluorouracil (5-FU) (; ; ; ; Ogundipe et al., 2023). Unfortunately, most patients experience recurrence and acquire resistance to platinum-based agents (). The inhibition of apoptosis is the mechanism by which chemoresistant tumor cells are propagated (; ; Neophytou et al., 2021; ); one of the causes is the inhibition of caspase activity, which is due to the changes and mutations in signaling pathways. Mammalian cancer cells will lose their promoters and have no ability to induce apoptosis, causing resistance to cytotoxic drugs (; ; ; ).
This type of ovarian cancer is characterized by genetic mutations in tumor suppressor genes TP53, BRCA1, or BRCA2 and by a specific dissemination mechanism through the body cavity known as transcoelomic metastasis. This mechanism involves rapid growth, disruption of ovarian tumor capsules, and malignant cells spreading into the peritoneal cavity by the ascites fluid, whose role is to offer the tumor microenvironment (Tan et al., 2006; ; Shield et al., 2009; The Cancer Genome Atlas Research Network, 2011; ; Suh et al., 2014; ; ; ). Tumor formation is driven by “tumor-initiating cells” that exhibit mesenchymal and stem cell features. Activation of EMT induces HGSOC precursor lesion (secretory cell outgrowths, SCOUTS, and serous tubal intraepithelial carcinoma, STIC) by suppressing paired box protein 2 (PAX2), a key molecule in maintaining the differentiation state of oviductal epithelial cells. A critical step in the progression of HGSOC is the migration of STIC cells to the ovary (; Perets et al., 2013; ; Tone, 2017). Moreover, growth factors and hormones are secreted onto the ovarian surface to induce EMT via phosphoinositide-3-kinase/Akt (PI3K/AKT) and mitogen-activated protein kinase/extracellular signal-regulated kinase (MEK/ERK) signaling pathways (Wong and Leung, 2007; ; ). In the metastatic stage, cancer cells are released from the primary tumor directly into the peritoneal cavity and survive either as single cells or spheroids in the ascites fluid, causing the formation of peritoneal nodules (peritoneal carcinomatosis) and also metastases to distant organs (Tan et al., 2006; Shield et al., 2009; ; ). EMT upregulates α5β1 integrin, which mediates spheroid attachment to the secondary site (; ). Moreover, Rosso et al. reported higher EMT marker expression in ascites cultures than in tumor cultures, indicating its crucial role in metastatic dissemination (Rosso et al., 2017). Nevertheless, even the mesenchymal-to-epithelial transition (MET) process has been described in the metastatic cascade, as associated with epigenetic abnormalities (). At the cellular level, the partial EMT promotes ascites and metastasis formation in HGSOC ().
Along with EMT markers, the cytokeratin family (KRTs), the most abundant proteins in epithelial cells, is pivotal in maintaining keratinization and differentiation. These have been reported to conserve cell morphology, intracellular transport, and signal transduction (; Zhang et al., 2020). In many cancers, such as ovarian, breast, and lung cancers, KRTs could be used as prognosis and tumorigenesis status markers (; Shao et al., 2012; ). Their expression has been associated with a high grade of malignancy and an increase in the migratory capability of cancer cells (Zhang et al., 2020).
An accumulating number of studies have highlighted the molecular heterogeneity of ovarian cancer, suggesting the need for personalized treatment approaches, including the establishment of ascites-derived cultures (RL et al., 2014; ; Penet et al., 2018; Uno et al., 2022). The ascites contains detached cancer cells, extracellular vesicles (EVs), tumor-associated macrophages (TAMs), and host cells, together promoting proliferation, drug resistance, or metastasis (Zhang et al., 2018). Hence, ascites reveals essential information about the underlying malignancy before resection that includes molecular mechanisms and profiles (). Considering these aspects, ascites presents a chance to design a treatment plan for patients with ovarian cancer by its potential use as a liquid biopsy substrate for exploring novel therapeutic targets (; ; ).
This study aimed to display the differences between the peritoneal nodule, primary tumor, and AS, the three major types of biological material derived from patients with HGSOC. We also conducted a computational target molecule prediction involved in chemoresistance. After in vitro drug testing of ascites-derived HGSOC primary cell culture, we reported a comprehensive investigation of the EM-associated genes.
2 Materials and methods
2.1 In silico analysis of available databases
2.1.1 The acquisition of mRNA expression datasets
The transcriptome profiles and relevant clinical information on patients with ovarian serous cystadenocarcinoma (OV) have been derived from The Cancer Genome Atlas (TCGA) (https://portal.gdc.cancer.gov/), and normal human ovarian samples have been obtained from The Genotype-Tissue Expression (GTEx) (https://www.gtexportal.org/home/datasets), using the TCGAbiolinks package () in the R program (version 4.3.0). Then, both RNA sequencing (RNA-seq) data (displayed as raw counts) were combined with batch normalization using the R package “sva” (Zhang et al., 2023). TCGA-OV dataset included 421 tumor samples, and the GTEx included 108 normal ovarian samples. TCGA and GTEx expression data were normalized in transcripts per million (TPM) format. The prognostic information on TCGA-OV samples was acquired from the UCSC Xena database (https://xenabrowser.net/).
2.1.2 Differential gene expression analysis
The differentially expressed protein-coding genes (DEGs) from tumor and normal ovarian tissues were generated using the DESeq2 (version 1.40.2) package of R software (). DEGs were selected based on a false discovery rate (FDR) p. adj< 0.01 and absolute logFC (fold change) ≥2. Principal component analysis was performed to examine relationships between tumoral and normal samples using a variance stabilizing transformation function to the count data (Wu et al., 2020) (Supplementary Figure S1). Next, the EnhancedVolcano (version 1.18.0) R package was used to visualize the results of differential expression analyses (Supplementary Figure S2). GO enrichment analysis of a gene set was performed using clusterProfiler (version 4.8.2) R package, and as a result, the significantly enriched GO terms were those with adjusted p-value <0.05.
Furthermore, we used the keywords “(drug-resistant) AND (EMT) AND (keratin) AND (ovarian cancer)” to search in the GeneCards (https://www.genecards.org/) database and obtained 1,265 drug resistance-EMT-keratin-related protein-coding genes. After the published literature was reviewed, a list of eight EM-associated genes, namely, five EMT—cadherin 1 (CDH1), cadherin 2 (CDH2), epithelial cell adhesion molecule (EPCAM), vimentin (VIM), GATA-binding protein 6 (GATA6) and three keratin markers—keratin 7 (KRT7), keratin 18 (KRT18), and keratin 19 (KRT19), was further explored. To construct and analyze the protein–protein interactions of our EM-associated genes, we submitted them to STRING (version 12, https://string-db.org/), a web-based open-access software tool.
2.1.3 Validation of EM-associated genes related to chemoresistance
In addition, the mRNA expression profile from HGSOC patients and ovarian cancer cell lines treated with platinum drugs was downloaded from the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) database using the GEOquery R packages, which were GSE227100, GSE58470, and GSE98559.
These three above-mentioned datasets were, respectively, derived from Homo sapiens mRNA data, using GPL24676 Illumina NovaSeq 6000, GPL16791 Illumina HiSeq 2500,and GPL6947 Illumina HumanHT-12 V3.0 expression bead chip.
The GSE227100 included a total of 24 HGSOC patients diagnosed with FIGO Stage III/IV, from which we obtained solid tumor samples before ((pre-C/T) and after (post-C/T) completion of six cycles of CARB and Taxol combination chemotherapy. We used GSE98559 with two experimental groups, SKOV3 wild-type cells and SKOV3 cisplatin-resistant cells, each having two biological replicate samples analyzed (four samples). We also explored three ovarian carcinoma cell lines: the parental cisplatin-sensitive IGROV-1 cell line and two platinum-resistant variants (IGROV-1 CIS/IGROV-1 OHP); each cell line had three independent samples analyzed (nine samples altogether) from GSE58470. Hierarchical clustering was performed with the pheatmap (version 1.0.12) R package using “Euclidean” clustering to calculate row distances and “complete” the agglomeration method (the distance between the most distant elements in each cluster).
2.2 Patient inclusion and sample collection
The study included a total of 12 AS, 9 T, and 7 NP from 12 HGSOC patients who underwent surgical resection at the Fundeni Clinical Institute between 2019 and 2021. The study was approved by the Ethics Committee of the Fundeni Clinical Institute (52496/06.12.2018). All the experiments were conducted following the Helsinki Declaration and obeying ethical principles for medical research on human subjects. Collected clinical and pathologic data include age, TNM stage, differentiation degree, tumor size (cm), serum tumor biomarkers (CA125 (ng/mL), CA15-3 (ng/mL), CEA (ng/mL), and CA19-9 (U/mL)), and overall survival (months). Histopathological results and the tumor grade were determined by the pathologist, according to the International Federation of Gynecology and Obstetrics (FIGO) classification. The cohort had not received preoperative chemotherapy at the time of debulking surgery. NP and T tissue samples were collected during surgery—in a stabilizing solution RNAlater (Sigma, St. Louis, MO)—then cryopreserved—using the snap-frozen method—and stored at −80 °C until further analysis.
2.3 Cell culture
2.3.1 SKOV3 cells
The human ovarian cancer cell line SKOV3 was purchased from the European Collection of Authenticated Cell Cultures (ECACC). SKOV3 was cultured in McCoy’s 5A medium (modified), w: L-glutamine, w: 2.2 g/L NaHCO3 (Biochrom), supplemented with 15% FBS-fetal bovine serum (Gibco) and 1% penicillin–streptomycin (P/S), as recommended by the suppliers.
2.3.2 HGSOC-AS
HGSOC ascites was collected under sterile conditions during surgery. Cell cultures were generated from 15 mL of fresh ascites seeded in the ratio 1:1 in Dulbecco’s modified Eagle medium-GlutaMAX (Gibco) containing 4.5 g/L glucose, supplemented with 20% FBS and 1% P/S, and incubated at 37 °C and 5% CO2. After 3 days, the cells were washed with Hank’s balanced salt solution (HBSS), and fresh media were added. Cells were harvested on early passages (up to three passage stocks).
2.4 RNA isolation and quantitative real-time PCR analysis (qRT-PCR)
For the characterization of ovarian cancer samples, we obtained total RNA using approximately 50 mg of fresh frozen tissue, which was first manually homogenized using a pestle and mortar, and then, we isolated RNA from around 3 × 105 cells for therapy response testing.
Total RNA from all sample types was isolated using TRIzol (Invitrogen, Thermo Fisher Scientific), according to the manufacturer’s protocol. The quantity of the RNA was determined by the concentration and purity (A260/A280 and A260/A230), assessed by NanoDrop ND1000 (NanoDrop Technologies, Waltham, MA, United States). Total RNA quality and size distribution were analyzed by chip-based capillary electrophoresis using the Agilent 2100 Bioanalyzer with a 6000 RNA Nano Chip (Agilent Technologies, Santa Clara, CA, United States).
The EM-associated gene expression was quantified by qRT-PCR. cDNA (2000 ng) was synthesized using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Thermo Fisher Scientific). Sense and antisense primers were designed against published human sequences in Supplementary Table S1. RT-qPCR was performed using the SYBR Green PCR Master Mix (Applied Biosystems, Thermo Fisher Scientific) and 7300 Real-Time PCR System (Thermo Fisher Scientific). The resulting mRNA levels were normalized to the β-actin reference gene. Relative quantification was studied by the 2−ΔCT method (; Schmittgen and Livak, 2008).
2.5 Immunohistochemistry (IHC)
Cell blocks have been prepared from 3 × 106 AS cells and processed using the ‘cell block’ cytology technique to generate paraffin blocks. The cell block technique has been done using cell suspension materials embedded in paraffin wax, according to the adapted method described in reference (Woods and Stirling, 2019). IHC staining for AS and T was performed on 3-µm-thick section cut from the formalin-fixed cell blocks and paraffin-embedded tumor tissues. Dewaxing and rehydration were completed using the Trilogy heated buffer solution (Cell Marque, Hot Springs, AR), according to the manufacturer’s protocol. In brief, 0.03% hydrogen peroxide (H2O2) treatment was performed for 10 min, and primary antibodies (CDH1 (Cell Marque; EP700Y; 1:100), CDH2 (Thermo Fisher Scientific; MA-1–91,128; 1:500), KRT18 (Cell Marque; B22.1&B23.1; 1:250), KRT19 (Cell Marque; A-53-B/A2.26; 1:250), and KRT7 (Cell Marque; OV-TL 12/30; 1:250) were incubated overnight at 4°C. The slides were then incubated with peroxidase-labeled polymer conjugated to goat anti-rabbit IgG or goat anti-mouse IgG for 30 min. The sections were stained with DAB and counterstained with hematoxylin. Cells were evaluated with a magnification of ×40. The tissues and AS were evaluated semiquantitatively, assessing the intensity and localization, including membrane, cytoplasm, or nuclear staining.
2.6 Functional study and chemotherapy response
2.6.1 Drug treatments
AS-derived cells have been seeded at 1 × 104 cells/well in 96 well-plates with flat-bottoms for 24 h. CIS (Selleckchem, catalog no. S1166, batch no 14) was dissolved in ddH2O, and DOX (Selleckchem, Cat no. S1208, Batch no 13) was dissolved in dimethyl sulfoxide (DMSO). Compound effects were measured in a 10-point dilution series after 24 h of incubation. CIS was tested in a range of 800 to 1.562 µM, and DOX, in a range of 25 to 0.048 µM. For CIS experiments, controls consisted of ddH2O, and for DOX experiments, controls consisted of DMSO alone (maximal DMSO concentration used was 0.016%). Cell viability was determined with MTT assay (Cell Proliferation Kit, MTT, Roche), and the solubilized formazan product was spectrophotometrically quantified at 570 nm wavelength using the Sunrise Basic Tecan plate reader and Magellan V 6.5 software. Two independent experiments with three technical replicates were conducted for each tested drug. Half-maximal inhibitory concentration (IC50) values were calculated with GraphPad Prism 10 software using a dose-response curve fit model by applying the nonlinear log(inhibitor) versus response-variable slope (four parameters) equation.
2.6.2 Scratch-healing test
AS cells have been seeded at 1 × 105/well in 12-well culture plates. After 24 h, a straight scratch was carefully made across the cell monolayer using a sterile 10-µL pipette tip to create a wound, and each drug’s median IC50 (DOX = 0.6 µM, CIS = 50 µM) was added. Subsequently, the dish was placed under the microscope (IX73 Inverted Microscope, Olympus) to capture images of the scratch at two different time points: immediately after the scratch (0 h) and 24 h later (24 h). These images were analyzed using ImageJ software. The extent of wound closure was measured by quantifying the reduction in scratch width over time.
2.6.3 Apoptosis assay
Cells have been lysed with Milliplex Map Lysis Buffer 1× (Merck Millipore), supplemented with protease inhibitors and assayed using ProcartaPlex multiplex immunoassay (Thermo Fisher Scientific), according to the manufacturer’s instructions. All samples were measured in duplicate. Data were initially acquired as mean fluorescence intensity (MIF), and the ratio of fluorescence to standard magnetic microspheres was then calculated. A series of calibrators were analyzed, and standard curves and concentrations were obtained using Bio-Plex Manager Software.
2.6.4 Western blotting
AS cells were lysed in RIPA buffer (MILLIPLEX MAP Lysis Buffer 1×, Merck Millipore) with the protease inhibitor (Protease Inhibitor Cocktail 50×, Promega). A measure of 30 μg of proteins were electrophoresed on a 15% and 8% SDS-polyacrylamide gel and then transferred to a polyvinylidene difluoride (PVDF) membrane (Bio-Rad Laboratories). After the blocking step with 5% nonfat dried milk (PanReac AppliChem) diluted in PBST (1× phosphate-buffered saline with Tween 20) for 2 h at room temperature, the primary antibodies mouse anti-KRT19 (A53-B/A2.26-Ks19.1 Thermo Fisher), rabbit anti-CDH2 (D4R14, Cell Signaling), and HRP-GAPDH (PA1-987 HRP, Invitrogen) were incubated overnight at 4°C. Afterward, the membrane was washed three times with PBST and incubated with goat anti-rabbit (G21234, Invitrogen) and m-IgGk BP-HRP (SC-516102, Santa Cruz) for 1 h at room temperature and washed again in PBST buffer. The signals were detected using Pierce ECL Western blotting Substrate (Thermo Fisher Scientific) and OPTIMAX X-ray film processor (Protec). The densitometric analysis quantified protein expression by ImageJ Software.
2.7 Statistical analysis
Quantitative and semiquantitative analyses for IHC tissue markers were performed with the support of experienced pathologists (VH).
In descriptive statistics, data are presented as n (%) or median (interquartile range (IQR): Q1 and Q3). Statistical significance of univariate analysis was determined by the Mann–Whitney–Wilcoxon test with p-values calculated by the exact method and the Kruskal–Wallis test for ordinal or continuous variables with a non-normal distribution. For normal distribution (evaluated with the Shapiro test), we used an unpaired t-test with Welch’s correction and one-way analysis of variance (ANOVA) and a chi-squared test for dichotomous variables, followed by the Bonferroni post hoc test for multiple comparisons. All p-values were based on two-sided hypothesis tests, and p <0.05 was considered statistically significant.
Wilcoxon tests were used for comparisons between two independent groups. AS results were presented as the means of two independent experiments ± SEM, and the comparisons between two groups (sensitive vs resistant) were performed using Mann–Whitney t-tests. Differences were considered statistically significant at p < 0.05.
We used log-transformed TPM values to normalize the mRNA expression of EM-associated markers to obtain a normal distribution of data for all RNA-seq transcriptomic results.
A Cox proportional hazards regression was used for univariate and multivariate analyses of prognostic variables for overall survival. Multivariate survival analysis was performed for all variables that indicated significant diagnostics, based on the scaled Schoenfeld residuals, to identify independent predictors of survival. The hazard ratio (HR) and the corresponding 95% confidence interval (95% CI) were calculated. Overall survival was defined as the interval between the date of surgery and the date of death or the end of follow-ups. The R packages used in the survival analysis are gtsummary 1.7.2, survival 3.5.7, survminer 0.4.9, and flextable 0.9.4. The statistical analysis used GraphPad Prism 10 for Windows (GraphPad Software, Inc.) and R 4.3.0 software.
Figure 1 presents the workflow strategy of this study.
FIGURE 1
3 Results
3.1 Investigation of EM markers related to chemoresistance
To explore the predictive biomarkers related to chemoresistance in HGSOC, we first performed in silico analysis. First, we obtained 1,278 upregulated genes and 1,060 downregulated genes in OV compared to normal ovarian tissues using TCGA and GTEx cohorts (Supplementary Table S2). Gene Ontology (GO) analysis found that cell–cell adhesion molecular functions (Figure 2A; Supplementary Table S3) and the positive regulation of cell–cell adhesion biological function were among the top enriched (Figure 2B; Supplementary Table S4). Notably, extracellular matrix and cellular adhesion were the top 10 ranked cellular components determined by GO analysis (Figure 2C; Supplementary Table S5). The most important molecular functions related to these genes are EMT and KRTs, defined as EM-associated genes.
FIGURE 2
After the intersection of TCGA-significant DEGs and EM-associated genes related to chemoresistance from the GeneCards database (Supplementary Table S6), we obtained 263 genes (Figure 2D; Supplementary Table S7). A total of eight genes (CDH1, CDH2, VIM, GATA6, EPCAM, KRT7, KRT18, and KRT19) were correlated with chemoresistance and HGSOC after relevant literature investigation (Supplementary Table S8) and were included in further analysis.
3.2 Generation of primary cultures from HGSOC-AS
The 12 patients included in this study had a median age of 59 (interquartile range 54–65 years old). All patients recruited were diagnosed with HGSOC. A total of 67% of patients were in stage IIIC. The median overall survival in this cohort was 23 months. Clinical and pathologic parameters are summarized in Table 1.
TABLE 1
| Characteristic | N | N = 12a |
| Age (years) | 12 | |
| Median (IQR) | 59 (54, 65) | |
| Differentiation degree | 12 | |
| G1 | 1 (8.3%) | |
| G2 | 1 (8.3%) | |
| G2–G3 | 2 (17%) | |
| G3 | 8 (67%) | |
| Staging (FIGO) | 12 | |
| IIA | 1 (8.3%) | |
| IIIA2 | 1 (8.3%) | |
| IIIB | 1 (8.3%) | |
| IIIC | 8 (67%) | |
| IVB | 1 (8.3%) | |
| Tumor size (cm) | 12 | |
| Median (IQR) | 4.3 (2.5, 7.3) | |
| Overall survival (months) | 12 | |
| Median (IQR) | 23 (11, 34) | |
| CA125 (ng/mL) | 11 | |
| Median (IQR) | 445 (392, 807) | |
| CA15-3 (ng/mL) | 10 | |
| Median (IQR) | 58 (23, 125) | |
| CEA (ng/mL) | 10 | |
| Median (IQR) | 1.41 (0.60, 1.86) | |
| CA19-9 (U/mL) | 10 | |
| Median (IQR) | 14 (4, 52) |
Demographic and clinical–pathological data on HGSOC patients included in the prospective study.
Median (IQR), percentage (%).
To provide relevant molecular aspects about the HGSOC microenvironment, we established 12 HGSOC ascites-derived primary cultures, as described elsewhere (Shepherd et al., 2006; Theriault et al., 2013). Their malignancy has been confirmed through cytological detection in the Pathological Anatomy Laboratory of the Fundeni Clinical Institute. Thus, rich cell sediment was observed, represented by tumor cells, red blood cells, frequent inflammatory cellular elements, and enlarged cells isolated and grouped together, some with three-dimensional uneven appearance and enlarged slightly uneven nuclei, which are often eccentric and tachromic; the cytoplasm was basophilic with vacuolations. Immediately after isolation (0 h), ascites-derived cells revealed wide heterogeneity, including suspended cell aggregates called spheroids. After 24 h (24 h), the presence of spheroids in suspension was still detectable, but adherent cells became predominant. After first passage, at 48 hours, 48 h (P+1), the confluent AS cell monolayer illustrates a “cobblestone” phenotype. The pattern of in vitro proliferation analysis shows variations, with primary cells derived from ascites reaching confluence in 6–10 days, depending on the biological variability of each patient (Figure 3A). AS, T, and NP samples of HGSOC patients were followed by subsequent mRNA expression analyses and in vitro functional tests.
FIGURE 3
3.3 Evaluation of EM-associated markers in ascites, tumors, and peritoneal nodule
To explore the epithelial–mesenchymal phenotype of AS, T, and NP, in our cohort and SKOV3 cell line, the expression of the most used EM-associated genes (CDH1, CDH2, EPCAM, KRT7, KRT18, and KRT19) has been evaluated by qRT-PCR. The AS has a significantly lower level of CDH1 compared to NP (p = 0.0214), the same as SKOV3 to NP (p = 0.0344). The EPCAM downregulation was detected in AS cells but not in SKOV3, with p = 0.0191 compared to T and p = 0.0030 compared to NP.
On the other hand, a significantly high mRNA level of CDH2 was observed in AS, compared to T (p = 0.0011) and NP (p = 0.0084). The expression of KRT18 was considerably higher in AS compared to T (p = 0.0059), which is in convergence with KRT19 in AS with p = 0.0014 compared to T and p = 0.0301 compared to NP. For KRT7, no significant differences were noticed between the four types of samples derived from HGSOC (Figure 3B).
Spearman’s rank correlation showed that there was a strong positive relation between the levels of CDH2—KRT7 and CDH2—EPCAM in AS and T, respectively, with r > 0.6 and p < 0.05. Moreover, all investigated markers, except KRT7–CDH1, KRT19–CDH1, and KRT7–KRT19, have a significant positive correlation in T (r > 0.6, p < 0.05). In contrast, only KRT markers exhibit a strong positive correlation in NP (r > 0.6, p < 0.05) (Figure 3C). The distinct patterns of association of each EM associated marker with cytokeratin demonstrate the unique expression profile of each type of biological material derived from HGSOC, highlighting that AS has the greatest heterogeneity and a mesenchymal phenotype.
Therefore, we also validated the results of marker expressions by IHC in tumor tissues and AS, respectively. The representative images for CDH1, CDH2, KRT18, KRT19, and KRT7 are shown in Figure 3D, and the IHC scores are detailed in Supplementary Table S9.
In addition to these, we analyzed TCGA-OV and GTEx against our list of EM-associated genes. Figure 3E shows that the EM-associated markers were significantly upregulated in OV tissues compared to normal ovarian tissues.
Furthermore, we acquired and visualized the protein–protein functional associations via the STRING database (https://www.string-db.org/) (von Mering et al., 2003). Therefore, we validated the interaction between the EMT and KRT markers (PPI enrichment p-value <0.001) (Figure 3F). Results show that AS presents a mesenchymal phenotype that enhances the aggressivity of HGSOC. Thus, this model is used for further analysis to predict drug response.
3.4 The establishment of chemosensitivity in the AS model
3.4.1 Determination of drug response in the AS model
To investigate the effect of drug response on AS, we have analyzed two important HGSOC drugs that induce apoptosis by different mechanisms () because cisplatin represses cell division and tumor growth by interfering with DNA replication and causing DNA cross-linking (Siddik, 2003; ; ), while doxorubicin intercalates with DNA, inhibits topoisomerase II, and promotes reactive species accumulation (Thorn et al., 2011; ).
First, we detected the cell viability using an MTT assay on 12 AS. Both drugs’ working concentrations (0.04875–25 µM DOX and 1.562–800 µM CIS) significantly decreased viability in all cultures tested in a dose-dependent manner (Figures 4A, B). Furthermore, we established the IC50 value after 24 h drug exposure. Thus, we used each drug’s median IC50 (DOX = 0.6 µM, CIS = 50 µM) for AS cultures in subsequent analyses.
FIGURE 4
In addition, for AS, the IC50 value of drugs equal to or less than their respective median IC50 was considered a sensitive culture, while IC50 values greater than the median IC50 values were regarded as a resistant culture.
Therefore, cell migration investigated using the scratch assay demonstrates a significant decrease in the CIS-sensitive group compared to its untreated control (p = 0.0022) (Figure 4D). In contrast, migration is not affected by DOX, either in resistant or sensitive groups (Figure 4C).
3.4.2 CASP3 expression is positively correlated with AS sensitivity to cisplatin and doxorubicin
To determine drug treatment resistance, apoptosis markers, caspase-3 (CASP3), and tumor protein P53 (p53) were measured using a multiplex apoptosis assay and the Luminex platform. As shown in Figure 4E, CASP3 and p53 significantly increased after DOX (p = 0.0091 and p = 0.0108, respectively) and CIS (p = 0.0024 and p = 0.0238, respectively) treatment (Figure 4G). Interestingly CASP3 levels were elevated in both sensitivity-drug groups (Figures 4F, H). These findings suggest that DOX and CIS repress proliferation and promote apoptosis by increasing CASP3 and p53 in the sensitive-AS group. Furthermore, CIS suppresses migration, while DOX does not affect it.
3.5 Drug treatment regulates EM-associated markers in AS
In order to explore the correlation of EM-associated markers and DOX and CIS ascites culture treatment, we have evaluated the mRNA expression by qRT-PCR. The results showed that VIM (p = 0.0076) and KRT18 (p = 0.0196) were significantly upregulated in CIS-treated AS compared to untreated-AS (control), while CDH1 (p = 0.0116), CDH2 (p = 0.0005), KRT19 (p = 0.0248), and GATA6 (p = 0.0389) were decreased (Figure 5A). In addition, upon treatment with DOX, the expression of CDH1 (p = 0.0205), VIM (p = 0.0458), and GATA6 (p = 0.0329) was increased, and the expression of KRT18 (p = 0.0376) and KRT19 (p = 0.0003) was downregulated (Figure 5B). CDH2 and KRT19 expression has been assessed by Western blotting, indicating a decreased level after CIS/DOX treatment (Figure 5 C, D).
FIGURE 5
Therefore, these data demonstrate that CIS treatment induces mesenchymal phenotype by VIM expression; this mechanism has been previously associated with chemoresistance (; Sun et al., 2012; ). Interestingly, only DOX increases CDH1 significantly, possibly inducing an invasive intermediate-EMT state ().
3.6 Different chemotherapies modulate EM-associated markers in HGSOC cell cultures and patients
We have validated the mRNA expression changes of EM-associated genes in drug-resistant and drug-sensitive cells, as well as in patients.
First, we investigated GSE58470, which includes IGROV-1 (parental cell, an established cell line to study chemoresistance in ovarian cancer properly) and two platinum-resistant variants (IGROV-1 CIS and IGROV-1 OHP) (). The heatmap in Figure 6A shows the correlation of EM-markers in platinum-resistant cells compared to IGROV-1. On one hand, CDH1, EPCAM, KRT7, and KRT19 were significantly lower in platinum-resistant variants than in parental cells (p < 0.05); on the other hand, KRT18 was upregulated in IGROV-1 CIS (p = 0.0385), while KRT7 increases in IGROV-1 OHP compared to IGROV-1 (p = 0.0057) and IGROV-1 CIS (p = 0.0001) (Figure 6C).
FIGURE 6
Furthermore, we analyzed the expression of EM-associated markers in SKOV3 CIS-resistant and wild-type SKOV3 cell lines using GSE98559 (; Zhou et al., 2018). Epithelial markers (KRT18, KRT7, CDH1, and EPCAM) were alleviated in SKOV3 CIS, compared to wild-type SKOV3, while mesenchymal markers (VIM) were highly expressed, except CDH2, which decreased after CIS exposure (Figure 6B).
Finally, these EM-associated markers were investigated in patient samples from GSE227100 (Figure 7A), which contained 24 HGSOC patients before and after they were treated with six cycles of CARB and PAX combination chemotherapy (pre-C/T; post-C/T) (). The results showed that the expression of VIM and GATA6 increased significantly after chemotherapies in the late recurrence group (Figure 7B).
FIGURE 7
These findings show that mRNA changes could differentiate between administrated treatments. Indeed, CIS thereby attenuates CDH1 and KRT19 expressions in all cells treated. In contrast, DOX contributes to the induction of the intermediate-EMT state by the expression of CDH1, VIM, and GATA6 (Figure 7C).
3.7 Clinical predictive value of EM-associated markers in HGSOC patients
To evaluate the clinical prognosis of the eight EM-associated genes, we have performed uni- and multivariate Cox regression analysis in OV patients from TCGA cohort. The univariate Cox proportional hazards regression analysis confirmed the higher risk of death for patients with a more advanced stage, stage IV (HR = 2.14, p = 0.047), high mRNA expression of the KRT7 marker (HR = 1.37, p = 0.012), and with the elevated mRNA expression of the KRT19 marker (HR = 1.33, p = 0.022). In addition, primary therapy outcome status in correlation with overall survival in univariate Cox regression indicated a lower risk of death for patients with stable disease (HR = 0.49, p = 0.036), complete remission/response (HR = 0.17, p < 0.001) compared to the progressive disease group, and high mRNA expression of EPCAM (HR = 0.73, p = 0.011) (Table 2). No association was detected between overall survival and the remaining EM-associated markers.
TABLE 2
| Summary data | Univariate | Multivariate | |||||||
| Characteristic | N = 420a | N | HRb | 95% CIb | p-value | N | HRb | 95% CIb | p-value |
| Tumor stage | 417 | 335 | |||||||
| Stage I/II | 25 (6.0%) | — | — | — | — | ||||
| Stage III | 328 (79%) | 1.73 | 0.85, 3.51 | 0.13 | 2.33 | 0.92, 5.87 | 0.073 | ||
| Stage IV | 64 (15%) | 2.14 | 1.01, 4.53 | 0.047 | 1.82 | 0.69, 4.82 | 0.2 | ||
| Primary therapy outcome | 337 | 335 | |||||||
| Progressive disease | 32 (9.5%) | — | — | — | — | ||||
| Partial remission/response | 47 (14%) | 0.72 | 0.44, 1.16 | 0.17 | 0.63 | 0.38, 1.05 | 0.075 | ||
| Stable disease | 23 (6.8%) | 0.49 | 0.25, 0.96 | 0.036 | 0.46 | 0.23, 0.93 | 0.031 | ||
| Complete remission/response | 235 (70%) | 0.17 | 0.11, 0.26 | <0.001 | 0.16 | 0.10, 0.25 | <0.001 | ||
| CDH1 mRNA expression | 420 | 335 | |||||||
| Low | 210 (50%) | — | — | — | — | ||||
| High | 210 (50%) | 1.09 | 0.85, 1.39 | 0.49 | 1.13 | 0.83, 1.54 | 0.4 | ||
| CDH2 mRNA expression | 420 | 335 | |||||||
| Low | 210 (50%) | — | — | — | — | ||||
| High | 210 (50%) | 0.85 | 0.67, 1.08 | 0.19 | 0.76 | 0.57, 1.02 | 0.071 | ||
| KRT7 mRNA expression | 420 | ||||||||
| Low | 210 (50%) | — | — | ||||||
| High | 210 (50%) | 1.37 | 1.07, 1.75 | 0.012 | |||||
| KRT18 mRNA expression | 420 | 335 | |||||||
| Low | 210 (50%) | — | — | — | — | ||||
| High | 210 (50%) | 1.08 | 0.85, 1.38 | 0.52 | 1.08 | 0.81, 1.45 | 0.6 | ||
| KRT19 mRNA expression | 420 | 335 | |||||||
| Low | 210 (50%) | — | — | — | — | ||||
| High | 210 (50%) | 1.33 | 1.04, 1.70 | 0.022 | 1.37 | 1.02, 1.84 | 0.035 | ||
| EPCAM mRNA expression | 420 | 335 | |||||||
| Low | 210 (50%) | — | — | — | — | ||||
| High | 210 (50%) | 0.73 | 0.57, 0.93 | 0.011 | 0.78 | 0.57, 1.05 | 0.10 | ||
| VIM mRNA expression | 420 | 335 | |||||||
| Low | 210 (50%) | — | — | — | — | ||||
| High | 210 (50%) | 0.96 | 0.75, 1.22 | 0.73 | 0.89 | 0.66, 1.19 | 0.4 | ||
| GATA6 mRNA expression | 420 | 335 | |||||||
| Low | 210 (50%) | — | — | — | — | ||||
| High | 210 (50%) | 1.11 | 0.87, 1.42 | 0.40 | 1.07 | 0.80, 1.44 | 0.6 | ||
Univariate and multivariate Cox proportional hazards regression analyses of EM-associated gene expression markers in correlation with the overall survival of HGSOC patients (TCGA).
Percentage (%).
HR, hazard ratio; C, confidence interval.
Finally, we investigated the significant variables to describe how they correlate with overall survival. To this end, we performed a multivariate Cox regression analysis using the proportional hazards assumption for the Cox model using statistical tests and graphical diagnostics based on the scaled Schoenfeld residuals, including all variables (Supplementary Figure S3).
The multivariate Cox proportional hazards regression analysis shows that a higher risk of death is directly associated with the high mRNA expression of the KRT19 marker (HR = 1.43, p = 0.018) and that a lower risk of death is correlated with stable disease (HR = 0.45, p = 0.028) and complete remission/response (HR = 0.15, p < 0.001) compared to the progressive disease group (Table 2). After controlling for the confounding factors, multivariate Cox regression analysis reveals that primary therapy outcomes and mRNA expression of the KRT19 marker are independent variables of overall survival.
4 Discussion
Growing evidence has suggested that overexpression of EMT and KRT markers is positively correlated with the progression and occurrence of various malignant carcinomas, including HGSOC (Xiao and He, 2010; Wang and Zhou, 2011; Rosso et al., 2017; ; ; Zhang et al., 2020; Sohn et al., 2021; ; ; ). Numerous studies have been conducted to establish the EMT index and its association with patient clinicopathologic features, aiming to improve the patient prognosis. Sohn HS et al. developed an EMT-TF-based prognostic index for patients by whole-exome and RNA sequencing. According to the study, the mesenchymal type, characterized by the activation of EMT-TFs and less genomic modification, is more aggressive than the homologous recombination repair (HRR)-activated type with deficiencies in HRR genes (Sohn et al., 2021). Using single-cell sequencing, Xu et al. found that the expression of NOTCH receptor 1 (NOTCH1), SNAI2, transforming growth factor beta receptor 1 (TGFBR1), and Wnt family member 11 (WNT11) EMT-associated genes is correlated with poor patient survival. In addition, primary matrix cancer-associated fibroblasts (mCAFs) can promote EMT and cell invasion by interacting with tumoral and immune cells (Xu et al., 2022). In addition, a recent study has revealed the correlation between N6-methyladenosine (m6A) modification regulators and EMT markers in OV development (Zhang et al., 2023).
However, no independent analysis of the predictive markers related to EM-associated genes and chemoresistance in HGSOC ascites has been conducted. This study found a new gene signature to guide therapy in this pathology, based on ascites, the most accessible sample. First, as a result of our analyses, the most relevant molecules from EMT and KRT families have been investigated in all sample types that could be collected from HGSOC patients (primary tumor, metastatic peritoneal nodules, and tumor ascitic fluids) and in SKOV3, an intermediate mesenchymal ovarian cell line (Rosso et al., 2017). These data show the first EM-marker analysis in the most relevant types of samples in HGSOC and indicate that ascites has a mesenchymal phenotype. In contrast, T and NP have epithelial phenotypes and similarities in gene expression. Li Y et al. have reported that CDH1 decreases in ascites, while CDH2, VIM, and KRT-19 increase in AS, compared to tumoral tissues (Rizvi et al., 2013; ), so our data conform to the previous reports.
As discussed above, the heterogeneity of HGSOC and ascites presence affects therapy efficacy (Pogge von Strandmann et al., 2017; Zhang et al., 2018). Many studies have used cell-free ascites and BRCAwt HGSOC patients’ tissues to identify new drug sensitivity biomarkers (; ; ). Moreover, Cook DP et al. used syngeneic models and tumors derived from mice to evaluate the origin, epigenetic, and phenotypic differences in order to establish the most useful model for research and therapeutic tests ().
In addition, the feasibility of chemoresistance tests in ascites cultures has been demonstrated (). Unlike other studies, we are the first to comprehensively analyze the hub genes related to EM and drug resistance in CIS/DOX-treated and untreated AS cultures at the transcriptional level. In addition, we investigated two datasets that contain RNA-seq of drug-resistant cells (SKOV3 and IGROV-1).
Our in vitro data indicate that apoptosis markers are elevated following both treatments, while only CIS affects migration in the CIS-sensitive group. Importantly, we observed a distinct CDH1, GATA6, and KRT18 pattern in response to drug exposure. Specifically, we observed a significant increase in CDH1 and GATA6 after DOX treatment and a notable decrease after CIS treatment. In addition to this, a slight decline was observed after OHP treatment. Leung D et al. showed that CDH1 is downregulated in CARB-resistant cells, leading to enhanced cellular migration and reduced proliferation (). Moreover GATA6 depletion correlates with the downregulation of epithelial markers (; ). In contrast, KRT18 was upregulated after CIS treatment in AS and IGROV-1, while DOX treatment determined a lower gene expression. Similarly, it has been shown that KRT18 downregulation improves CIS sensitivity and EMT-independent collective migration in epithelial cancer cells (). This mechanism can explain the difference in the expression of KRT18 and CDH1 between the two types of treatments; KRT18 decreases and CDH1 increases after DOX treatment, resulting in unaffected migration.
Furthermore, we observed a significant increase in VIM and a decrease in KRT19 after both CIS and DOX treatment in AS. On the contrary, VIM, an intermediate filament protein that preserves cellular integrity, has been reported to decrease in drug-resistant ovarian cells. Its depletion leads to CIS resistance via the downregulation of cytoskeleton organization proteins, inducing the cancer stem cell phenotype and prolonged G2 arrest in drug-resistant cells (; ; ). Among these studies, Latif et al. found that VIM was upregulated at both transcriptional and post-translational levels after the CIS treatment of metastatic epithelial ovarian tumor cells, increasing stemness and drug resistance (). These discrepancies can rely on the phenotype of the cells treated. Indeed, the upregulation in VIM levels may be due to the treatment of mesenchymal cells. Therefore, designing a successful treatment strategy for advanced ovarian cancer requires investigation of the ascites microenvironment because it contains relevant molecules that lead to chemoresistance and disease progression. In addition to all these, the significant drawback of ovarian cancer cell lines is due to the accumulation of various genetic and phenotypic defects during years of culture, which no longer correctly reflect the clinical condition ().
Consistent with previous observations (Woopen et al., 2014; ; Sun et al., 2023), our analysis revealed that KRT19 and KRT7 negatively correlate with survival outcomes, while EPCAM has a positive correlation. KRT19 was reported as a potential immunotherapy target (Sun et al., 2023), and overexpression of this molecule promotes the proliferation and migration of cancer cells via Wnt/β-catenin signaling (). In contrast, downregulation of KRT19 has been reported in tumor breast tissue compared to adjacent tissue, and it has been associated with an aggressive phenotype and chemoresistance of cancer cells (; Saha et al., 2017; Saha et al., 2018).
At least, only VIM and GATA6 expressions have shown the same signature in tissue samples before/after CARB and PAX combination chemotherapy. Indeed, a lower expression has been obtained explicitly in the late-recurrence group, but a higher expression, after drug exposure.
To sum up, our results indicate that after DOX and CIS treatment, AS acquires an invasive-intermediate EMT state by the overexpression of both epithelial and mesenchymal markers. In both cases, VIM overexpression and KRT19 downregulation have been reported. Further research will investigate how these genes are involved in AS chemoresistance by expanding the sample size and conducting basic experiments. However, the ascites remains an exclusive and accessible sample from HGSOC patients to explore tumor progression and molecular pathways involved in chemoresistance. This could lead us to improve personalized treatment decisions by developing new targets to overcome drug resistance in HGSOC.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the Ethics Committee of Fundeni Clinical Institute (52496/06.12.2018). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
D-RC: conceptualization, data curation, formal analysis, investigation, methodology, validation, visualization, writing–original draft, and writing–review and editing. AS: conceptualization, data curation, formal analysis, investigation, methodology, software, validation, visualization, writing–original draft, and Writing–review and editing. A-VG: writing–original draft, writing–review and editing, conceptualization, data curation, formal analysis, investigation, validation, and visualization. DL: conceptualization, data curation, formal analysis, investigation, methodology, writing–original draft, and writing–review and editing. VH: conceptualization, funding acquisition, methodology, project administration, resources, supervision, writing–original draft, and writing–review and editing. NB: conceptualization, funding acquisition, methodology, project administration, resources, supervision, writing–original draft, and writing–review and editing. SD: conceptualization, funding acquisition, project administration, resources, supervision, visualization, writing–original draft, and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was founded by “Carol Davila” University of Medicine and Pharmacy Bucharest, Romania, through contract number 33PFE/ 30.12.2021 funded by the Ministry of Research and Innovation within PNCDI III, Program 1–Development of the National RD system, Subprogram 1.2–Institutional Performance–RDI excellence funding projects; by a grant with contract number 760009/30.12.2022, code: 14/16.11.2022, funded by Romania’s National Recovery and Resilience Plan (PNRR), Pylon III, section I5. Establishment and operationalization of Competence Centers PNRR-III-C9-2022 – I5, project “Creation, Operational and Development of the National Center of Competence in the field of Cancer”, acronym CNCC; by a grant of the Romanian Ministry of Research and Innovation, CCDI-UEFISCDI, project number PNIII-P1-1.2-PCCDI2017-0833/’’Multidisciplinary Consortium for Supporting Research Skills in Diagnosing, Treating, and Identifying Predictive Factors of Malignant Gynecologic Disorders-ONCOGIN,’’ contract number 68PCCDI/2018.3.
Acknowledgments
The authors would like to sincerely acknowledge TCGA, GTEx, GEO, and the developers of these datasets.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1363142/full#supplementary-material
References
1
AdzibolosuN.AlveroA. B.Ali-FehmiR.GogoiR.CoreyL.TedjaR.et al (2023). Immunological modifications following chemotherapy are associated with delayed recurrence of ovarian cancer. Front. Immunol.14, 1204148. 10.3389/fimmu.2023.1204148
2
AhmedN.StenversK. L. (2013). Getting to know ovarian cancer ascites: opportunities for targeted therapy-based translational research. Front. Oncol.3, 256. 10.3389/fonc.2013.00256
3
Al HabyanS.KalosC.SzymborskiJ.MccaffreyL. (2018). Multicellular detachment generates metastatic spheroids during intra-abdominal dissemination in epithelial ovarian cancer. Oncogene37, 5127–5135. 10.1038/s41388-018-0317-x
4
AlwosaibaiK.AbediniA.Al-HujailyE. M.TangY.GarsonK.CollinsO.et al (2017). PAX2 maintains the differentiation of mouse oviductal epithelium and inhibits the transition to a stem cell-like state. Oncotarget8, 76881–76897. 10.18632/oncotarget.20173
5
ArmstrongD. K.AlvarezR. D.Bakkum-GamezJ. N.BarroilhetL.BehbakhtK.BerchuckA.et al (2021). Ovarian cancer, version 2.2020, NCCN clinical practice guidelines in oncology. J. Natl. Compr. Canc Netw.19, 191–226. 10.6004/jnccn.2021.0007
6
ArrighettiN.CossaG.De CeccoL.StucchiS.CareniniN.CornaE.et al (2016). PKC-alpha modulation by miR-483-3p in platinum-resistant ovarian carcinoma cells. Toxicol. Appl. Pharmacol.310, 9–19. 10.1016/j.taap.2016.08.005
7
BacalbasaN.BalescuI.DimaS.VilcuM.BrezeanI. (2020a). Extended pelvic resections as part of secondary cytoreduction for relapsed ovarian adenocarcinoma. Vivo34, 397–400. 10.21873/invivo.11787
8
BacalbasaN.BalescuI.VilcuM.DimaS.DiaconuC.IliescuL.et al (2020b). The risk of para-aortic lymph node metastases in apparent early stage ovarian cancer. Med. Kaunas.56, 108. 10.3390/medicina56030108
9
BarbolinaM. V. (2018). Molecular mechanisms regulating organ-specific metastases in epithelial ovarian carcinoma. Cancers (Basel)10, 444. 10.3390/cancers10110444
10
BlobelG. A.MollR.FrankeW. W.Vogt-MoykopfI. (1984). Cytokeratins in normal lung and lung carcinomas. I. Adenocarcinomas, squamous cell carcinomas and cultured cell lines. Virchows Arch. B Cell Pathol. Incl. Mol. Pathol.45, 407–429. 10.1007/BF02889883
11
BowtellD. D.BohmS.AhmedA. A.AspuriaP. J.BastR. C.JR.BeralV.et al (2015). Rethinking ovarian cancer II: reducing mortality from high-grade serous ovarian cancer. Nat. Rev. Cancer15, 668–679. 10.1038/nrc4019
12
BrabletzT.KalluriR.NietoM. A.WeinbergR. A. (2018). EMT in cancer. Nat. Rev. Cancer18, 128–134. 10.1038/nrc.2017.118
13
ButtarelliM.CiucciA.PalluzziF.RaspaglioG.MarchettiC.PerroneE.et al (2022). Identification of a novel gene signature predicting response to first-line chemotherapy in BRCA wild-type high-grade serous ovarian cancer patients. J. Exp. Clin. Cancer Res.41, 50. 10.1186/s13046-022-02265-w
14
CaccuriF.SommarivaM.MarsicoS.GiordanoF.ZaniA.GiacominiA.et al (2019). Inhibition of DNA repair mechanisms and induction of apoptosis in triple negative breast cancer cells expressing the human herpesvirus 6 U94. Cancers (Basel)11, 1006. 10.3390/cancers11071006
15
Capo-ChichiC. D.RolandI. H.VanderveerL.BaoR.YamagataT.HiraiH.et al (2003). Anomalous expression of epithelial differentiation-determining GATA factors in ovarian tumorigenesis. Cancer Res.63, 4967–4977.
16
ChenE. Y.MehraK.MehradM.NingG.MironA.MutterG. L.et al (2010). Secretory cell outgrowth, PAX2 and serous carcinogenesis in the Fallopian tube. J. Pathol.222, 110–116. 10.1002/path.2739
17
ColapricoA.SilvaT. C.OlsenC.GarofanoL.CavaC.GaroliniD.et al (2016). TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res.44, e71. 10.1093/nar/gkv1507
18
CommunalL.RoyN.CahuzacM.RahimiK.KobelM.ProvencherD. M.et al (2021). A keratin 7 and E-cadherin signature is highly predictive of tubo-ovarian high-grade serous carcinoma prognosis. Int. J. Mol. Sci.22, 5325. 10.3390/ijms22105325
19
CookD. P.GalpinK. J. C.RodriguezG. M.ShakfaN.Wilson-SanchezJ.EchaibiM.et al (2023). Comparative analysis of syngeneic mouse models of high-grade serous ovarian cancer. Commun. Biol.6, 1152. 10.1038/s42003-023-05529-z
20
DeanM.DavisD. A.BurdetteJ. E. (2017). Activin A stimulates migration of the fallopian tube epithelium, an origin of high-grade serous ovarian cancer, through non-canonical signaling. Cancer Lett.391, 114–124. 10.1016/j.canlet.2017.01.011
21
Den OudenJ. E.ZamanG. J. R.DylusJ.Van DoornmalenA. M.MulderW. R.GrobbenY.et al (2020). Chemotherapy sensitivity testing on ovarian cancer cells isolated from malignant ascites. Oncotarget11, 4570–4581. 10.18632/oncotarget.27827
22
DomckeS.SinhaR.LevineD. A.SanderC.SchultzN. (2013). Evaluating cell lines as tumour models by comparison of genomic profiles. Nat. Commun.4, 2126. 10.1038/ncomms3126
23
DuanC.YuM.XuJ.LiB. Y.ZhaoY.KankalaR. K. (2023). Overcoming cancer multi-drug resistance (MDR): reasons, mechanisms, nanotherapeutic solutions, and challenges. Biomed. Pharmacother.162, 114643. 10.1016/j.biopha.2023.114643
24
FanH.AtiyaH. I.WangY.PisanicT. R.WangT. H.ShihI. M.et al (2020). Epigenomic reprogramming toward mesenchymal-epithelial transition in ovarian-cancer-associated mesenchymal stem cells drives metastasis. Cell Rep.33, 108473. 10.1016/j.celrep.2020.108473
25
FordC. E.WernerB.HackerN. F.WartonK. (2020). The untapped potential of ascites in ovarian cancer research and treatment. Br. J. Cancer123, 9–16. 10.1038/s41416-020-0875-x
26
FortierA. M.AsselinE.CadrinM. (2013). Keratin 8 and 18 loss in epithelial cancer cells increases collective cell migration and cisplatin sensitivity through claudin1 up-regulation. J. Biol. Chem.288, 11555–11571. 10.1074/jbc.M112.428920
27
GaoJ.ZhuY.NilssonM.SundfeldtK. (2014). TGF-β isoforms induce EMT independent migration of ovarian cancer cells. Cancer Cell Int.14, 72. 10.1186/s12935-014-0072-1
28
GhavamiS.HashemiM.AndeS. R.YeganehB.XiaoW.EshraghiM.et al (2009). Apoptosis and cancer: mutations within caspase genes. J. Med. Genet.46, 497–510. 10.1136/jmg.2009.066944
29
GhionescuA. V.SoropA.DimaS. O. (2023). The pivotal role of EMT-related noncoding RNAs regulatory axes in hepatocellular carcinoma. Front. Pharmacol.14, 1270425. 10.3389/fphar.2023.1270425
30
HelmC. W.StatesJ. C. (2009). Enhancing the efficacy of cisplatin in ovarian cancer treatment - could arsenic have a role. J. Ovarian Res.2, 2. 10.1186/1757-2215-2-2
31
HenriJ. L.NakhjavaniM.MccoombeS.ShigdarS. (2023). Cytotoxic effects of aptamer-doxorubicin conjugates in an ovarian cancer cell line. Biochimie204, 108–117. 10.1016/j.biochi.2022.09.005
32
HuoY.ZhengZ.ChenY.WangQ.ZhangZ.DengH. (2016). Downregulation of vimentin expression increased drug resistance in ovarian cancer cells. Oncotarget7, 45876–45888. 10.18632/oncotarget.9970
33
JacobJ. T.CoulombeP. A.KwanR.OmaryM. B. (2018). Types I and II keratin intermediate filaments. Cold Spring Harb. Perspect. Biol.10, a018275. 10.1101/cshperspect.a018275
34
JanR.ChaudhryG. E. (2019). Understanding apoptosis and apoptotic pathways targeted cancer therapeutics. Adv. Pharm. Bull.9, 205–218. 10.15171/apb.2019.024
35
JiangM.QiL.LiL.LiY. (2020). The caspase-3/GSDME signal pathway as a switch between apoptosis and pyroptosis in cancer. Cell Death Discov.6, 112. 10.1038/s41420-020-00349-0
36
JinL.HuoY.ZhengZ.JiangX.DengH.ChenY.et al (2014). Down-regulation of Ras-related protein Rab 5C-dependent endocytosis and glycolysis in cisplatin-resistant ovarian cancer cell lines. Mol. Cell Proteomics13, 3138–3151. 10.1074/mcp.M113.033217
37
JuJ. H.YangW.LeeK. M.OhS.NamK.ShimS.et al (2013). Regulation of cell proliferation and migration by keratin19-induced nuclear import of early growth response-1 in breast cancer cells. Clin. Cancer Res.19, 4335–4346. 10.1158/1078-0432.CCR-12-3295
38
KanakkantharaA.RawsonP.NorthcoteP. T.MillerJ. H. (2012). Acquired resistance to peloruside A and laulimalide is associated with downregulation of vimentin in human ovarian carcinoma cells. Pharm. Res.29, 3022–3032. 10.1007/s11095-012-0773-x
39
KerrE. H.FrederickP. J.EggerM. E.StockardC. R.SellersJ.DellamannaD.et al (2013). Lung resistance-related protein (LRP) expression in malignant ascitic cells as a prognostic marker for advanced ovarian serous carcinoma. Ann. Surg. Oncol.20, 3059–3065. 10.1245/s10434-013-2878-9
40
KimS.KimB.SongY. S. (2016). Ascites modulates cancer cell behavior, contributing to tumor heterogeneity in ovarian cancer. Cancer Sci.107, 1173–1178. 10.1111/cas.12987
41
KippsE.TanD. S.KayeS. B. (2013). Meeting the challenge of ascites in ovarian cancer: new avenues for therapy and research. Nat. Rev. Cancer13, 273–282. 10.1038/nrc3432
42
KlymenkoY.JohnsonJ.BosB.LombardR.CampbellL.LoughranE.et al (2017). Heterogeneous cadherin expression and multicellular aggregate dynamics in ovarian cancer dissemination. Neoplasia19, 549–563. 10.1016/j.neo.2017.04.002
43
KunjachanS.RychlikB.StormG.KiesslingF.LammersT. (2013). Multidrug resistance: physiological principles and nanomedical solutions. Adv. Drug Deliv. Rev.65, 1852–1865. 10.1016/j.addr.2013.09.018
44
KurmanR. J.Shih IeM. (2008). Pathogenesis of ovarian cancer: lessons from morphology and molecular biology and their clinical implications. Int. J. Gynecol. Pathol.27, 151–160. 10.1097/PGP.0b013e318161e4f5
45
LaneD.MatteI.Garde-GrangerP.LaplanteC.CarignanA.RancourtC.et al (2015). Inflammation-regulating factors in ascites as predictive biomarkers of drug resistance and progression-free survival in serous epithelial ovarian cancers. BMC Cancer15, 492. 10.1186/s12885-015-1511-7
46
LatifiA.AbubakerK.CastrechiniN.WardA. C.LiongueC.DobillF.et al (2011). Cisplatin treatment of primary and metastatic epithelial ovarian carcinomas generates residual cells with mesenchymal stem cell-like profile. J. Cell Biochem.112, 2850–2864. 10.1002/jcb.23199
47
LatifiA.LuworR. B.BilandzicM.NazaretianS.StenversK.PymanJ.et al (2012). Isolation and characterization of tumor cells from the ascites of ovarian cancer patients: molecular phenotype of chemoresistant ovarian tumors. PLoS One7, e46858. 10.1371/journal.pone.0046858
48
LengyelE. (2010). Ovarian cancer development and metastasis. Am. J. Pathol.177, 1053–1064. 10.2353/ajpath.2010.100105
49
LeungD.PriceZ. K.LokmanN. A.WangW.GoonetillekeL.KadifeE.et al (2022). Platinum-resistance in epithelial ovarian cancer: an interplay of epithelial-mesenchymal transition interlinked with reprogrammed metabolism. J. Transl. Med.20, 556. 10.1186/s12967-022-03776-y
50
LheureuxS.BraunsteinM.OzaA. M. (2019). Epithelial ovarian cancer: evolution of management in the era of precision medicine. CA Cancer J. Clin.69, 280–304. 10.3322/caac.21559
51
LiQ. Q.ChenZ. Q.CaoX. X.XuJ. D.XuJ. W.ChenY. Y.et al (2011). Involvement of NF-κB/miR-448 regulatory feedback loop in chemotherapy-induced epithelial-mesenchymal transition of breast cancer cells. Cell Death Differ.18, 16–25. 10.1038/cdd.2010.103
52
LiJ.LiuH.YuJ.YuH. (2015). Chemoresistance to doxorubicin induces epithelial-mesenchymal transition via upregulation of transforming growth factor β signaling in HCT116 colon cancer cells. Mol. Med. Rep.12, 192–198. 10.3892/mmr.2015.3356
53
LisioM. A.FuL.GoyenecheA.GaoZ. H.TelleriaC. (2019). High-grade serous ovarian cancer: basic sciences, clinical and therapeutic standpoints. Int. J. Mol. Sci.20, 952. 10.3390/ijms20040952
54
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods25, 402–408. 10.1006/meth.2001.1262
55
LiX.TangM.ZhuQ.WangX.LinY.WangX. (2020). The exosomal integrin α5β1/AEP complex derived from epithelial ovarian cancer cells promotes peritoneal metastasis through regulating mesothelial cell proliferation and migration. Cell Oncol. (Dordr)43, 263–277. 10.1007/s13402-019-00486-4
56
LiY.FeiH.LinQ.LiangF.YouY.LiM.et al (2021). ZEB2 facilitates peritoneal metastasis by regulating the invasiveness and tumorigenesis of cancer stem-like cells in high-grade serous ovarian cancers. Oncogene40, 5131–5141. 10.1038/s41388-021-01913-3
57
LoretN.DenysH.TummersP.BerxG. (2019). The role of epithelial-to-mesenchymal plasticity in ovarian cancer progression and therapy resistance. Cancers11, 838. 10.3390/cancers11060838
58
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.15, 550. 10.1186/s13059-014-0550-8
59
LuJ.KornmannM.TraubB.HuY.ZengZ. (2023). Preparation and antimicrobial activity of a film-forming polyhexamethylene biguanide teat disinfectant. Int. J. Mol. Sci.24, 17444. 10.3390/ijms242417444
60
LuQ.QuH.LouT.LiuC.ZhangZ. (2020). CK19 promotes ovarian cancer development by impacting on wnt/β-catenin pathway. Onco Targets Ther.13, 2421–2431. 10.2147/OTT.S242778
61
MachinoH.DozenA.KonakaM.KomatsuM.NakamuraK.IkawaN.et al (2023). Integrative analysis reveals early epigenetic alterations in high-grade serous ovarian carcinomas. Exp. Mol. Med.55, 2205–2219. 10.1038/s12276-023-01090-1
62
MakovecT. (2019). Cisplatin and beyond: molecular mechanisms of action and drug resistance development in cancer chemotherapy. Radiol. Oncol.53, 148–158. 10.2478/raon-2019-0018
63
MartinelliP.Carrillo-De Santa PauE.CoxT.SainzB.DusettiN.GreenhalfW.et al (2017). GATA6 regulates EMT and tumour dissemination, and is a marker of response to adjuvant chemotherapy in pancreatic cancer. Gut66, 1665–1676. 10.1136/gutjnl-2015-311256
64
MazidimoradiA.MomenimovahedZ.AllahqoliL.TiznobaikA.HajinasabN.SalehiniyaH.et al (2022). The global, regional and national epidemiology, incidence, mortality, and burden of ovarian cancer. Health Sci. Rep.5, e936. 10.1002/hsr2.936
65
MengY.ChenC. W.YungM. M. H.SunW.SunJ.LiZ.et al (2018). DUOXA1-mediated ROS production promotes cisplatin resistance by activating ATR-Chk1 pathway in ovarian cancer. Cancer Lett.428, 104–116. 10.1016/j.canlet.2018.04.029
66
MohamedM. S.BishrM. K.AlmutairiF. M.AliA. G. (2017). Inhibitors of apoptosis: clinical implications in cancer. Apoptosis22, 1487–1509. 10.1007/s10495-017-1429-4
67
MomenimovahedZ.TiznobaikA.TaheriS.SalehiniyaH. (2019). Ovarian cancer in the world: epidemiology and risk factors. Int J Womens Health11, 287–299. 10.2147/IJWH.S197604
68
MorganS.LopesF.GourleyC.AndersonR. A.SpearsN. (2013). Cisplatin and doxorubicin induce distinct mechanisms of ovarian follicle loss; imatinib provides selective protection only against cisplatin. PLoS One8, e70117. 10.1371/journal.pone.0070117
69
NeophytouC. M.TrougakosI. P.ErinN.PapageorgisP. (2021). Apoptosis deregulation and the development of cancer multi-drug resistance. Cancers (Basel)15, 4363. 10.3390/cancers13174363
70
OgundipeO. D.OlajubutuO.AdesinaS. K. (2023). Targeted drug conjugate systems for ovarian cancer chemotherapy. Biomed. Pharmacother.165, 115151. 10.1016/j.biopha.2023.115151
71
PenetM. F.KrishnamacharyB.WildesF. B.MironchikY.HungC. F.WuT. C.et al (2018). Ascites volumes and the ovarian cancer microenvironment. Front. Oncol.8, 595. 10.3389/fonc.2018.00595
72
PeretsR.WyantG. A.MutoK. W.BijronJ. G.PooleB. B.ChinK. T.et al (2013). Transformation of the fallopian tube secretory epithelium leads to high-grade serous ovarian cancer in Brca;Tp53;Pten models. Cancer Cell24, 751–765. 10.1016/j.ccr.2013.10.013
73
Pogge Von StrandmannE.ReinartzS.WagerU.MullerR. (2017). Tumor-host cell interactions in ovarian cancer: pathways to therapy failure. Trends Cancer3, 137–148. 10.1016/j.trecan.2016.12.005
74
RizviI.GurkanU. A.TasogluS.AlagicN.CelliJ. P.MensahL. B.et al (2013). Flow induces epithelial-mesenchymal transition, cellular heterogeneity and biomarker modulation in 3D ovarian cancer nodules. Proc. Natl. Acad. Sci. U. S. A.110, E1974–E1983. 10.1073/pnas.1216989110
75
RlO. D.MccormickA.MukhopadhyayA.WoodhouseL. C.MoatM.GrundyA.et al (2014). The use of ovarian cancer cells from patients undergoing surgery to generate primary cultures capable of undergoing functional analysis. PLoS One9, e90604. 10.1371/journal.pone.0090604
76
RossoM.MajemB.DevisL.LapyckyjL.BessoM. J.LlauradoM.et al (2017). E-cadherin: a determinant molecule associated with ovarian cancer progression, dissemination and aggressiveness. PLoS One12, e0184439. 10.1371/journal.pone.0184439
77
SahaS. K.ChoiH. Y.KimB. W.DayemA. A.YangG. M.KimK. S.et al (2017). KRT19 directly interacts with β-catenin/RAC1 complex to regulate NUMB-dependent NOTCH signaling pathway and breast cancer properties. Oncogene36, 332–349. 10.1038/onc.2016.221
78
SahaS. K.KimK.YangG. M.ChoiH. Y.ChoS. G. (2018). Cytokeratin 19 (KRT19) has a role in the reprogramming of cancer stem cell-like cells to less aggressive and more drug-sensitive cells. Int. J. Mol. Sci.19, 1423. 10.3390/ijms19051423
79
SchmittgenT. D.LivakK. J. (2008). Analyzing real-time PCR data by the comparative C(T) method. Nat. Protoc.3, 1101–1108. 10.1038/nprot.2008.73
80
ShaoM. M.ChanS. K.YuA. M.LamC. C.TsangJ. Y.LuiP. C.et al (2012). Keratin expression in breast cancers. Virchows Arch.461, 313–322. 10.1007/s00428-012-1289-9
81
ShepherdT. G.TheriaultB. L.CampbellE. J.NachtigalM. W. (2006). Primary culture of ovarian surface epithelial cells and ascites-derived ovarian cancer cells from patients. Nat. Protoc.1, 2643–2649. 10.1038/nprot.2006.328
82
ShieldK.AcklandM. L.AhmedN.RiceG. E. (2009). Multicellular spheroids in ovarian cancer metastases: biology and pathology. Gynecol. Oncol.113, 143–148. 10.1016/j.ygyno.2008.11.032
83
SiddikZ. H. (2003). Cisplatin: mode of cytotoxic action and molecular basis of resistance. Oncogene22, 7265–7279. 10.1038/sj.onc.1206933
84
SohnM. H.KimS. I.ShinJ. Y.KimH. S.ChungH. H.KimJ. W.et al (2021). Classification of high-grade serous ovarian carcinoma by epithelial-to-mesenchymal transition signature and homologous recombination repair genes. Genes (Basel)12, 1103. 10.3390/genes12071103
85
SuhD. H.KimH. S.KimB.SongY. S. (2014). Metabolic orchestration between cancer cells and tumor microenvironment as a co-evolutionary source of chemoresistance in ovarian cancer: a therapeutic implication. Biochem. Pharmacol.92, 43–54. 10.1016/j.bcp.2014.08.011
86
SungH.FerlayJ.SiegelR. L.LaversanneM.SoerjomataramI.JemalA.et al (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin.71, 209–249. 10.3322/caac.21660
87
SunL.YaoY.LiuB.LinZ.LinL.YangM.et al (2012). MiR-200b and miR-15b regulate chemotherapy-induced epithelial-mesenchymal transition in human tongue cancer cells by targeting BMI1. Oncogene31, 432–445. 10.1038/onc.2011.263
88
SunZ.ZhouR.DaiJ.ChenJ.LiuY.WangM.et al (2023). KRT19 is a promising prognostic biomarker and associates with immune infiltrates in serous ovarian cystadenocarcinoma. Int. J. Gen. Med.16, 4849–4862. 10.2147/IJGM.S419235
89
TanD. S.AgarwalR.KayeS. B. (2006). Mechanisms of transcoelomic metastasis in ovarian cancer. Lancet Oncol.7, 925–934. 10.1016/S1470-2045(06)70939-1
90
The Cancer Genome Atlas Research Network (2011). Integrated genomic analyses of ovarian carcinoma. Nature474, 609–615. 10.1038/nature10166
91
TheriaultB. L.PortelanceL.Mes-MassonA. M.NachtigalM. W. (2013). Establishment of primary cultures from ovarian tumor tissue and ascites fluid. Methods Mol. Biol.1049, 323–336. 10.1007/978-1-62703-547-7_24
92
ThornC. F.OshiroC.MarshS.Hernandez-BoussardT.McleodH.KleinT. E.et al (2011). Doxorubicin pathways: pharmacodynamics and adverse effects. Pharmacogenet Genomics21, 440–446. 10.1097/FPC.0b013e32833ffb56
93
ToneA. A. (2017). Taking the tube: from normal fallopian tube epithelium to ovarian high-grade serous carcinoma. Clin. Obstet. Gynecol.60, 697–710. 10.1097/GRF.0000000000000313
94
UnoK.IyoshiS.YoshiharaM.KitamiK.MogiK.FujimotoH.et al (2022). Metastatic voyage of ovarian cancer cells in ascites with the assistance of various cellular components. Int. J. Mol. Sci.23, 4383. 10.3390/ijms23084383
95
Von MeringC.HuynenM.JaeggiD.SchmidtS.BorkP.SnelB. (2003). STRING: a database of predicted functional associations between proteins. Nucleic Acids Res.31, 258–261. 10.1093/nar/gkg034
96
WangY.ZhouB. P. (2011). Epithelial-mesenchymal transition in breast cancer progression and metastasis. Chin. J. Cancer30, 603–611. 10.5732/cjc.011.10226
97
WongA. S.LeungP. C. (2007). Role of endocrine and growth factors on the ovarian surface epithelium. J. Obstet. Gynaecol. Res.33, 3–16. 10.1111/j.1447-0756.2007.00478.x
98
WoodsA. E.StirlingJ. W. (2019). Transmission electron microscopy. Bancroft's Theory Pract. Histol. Tech.
99
WoopenH.PietznerK.RichterR.FotopoulouC.JoensT.BraicuE. I.et al (2014). Overexpression of the epithelial cell adhesion molecule is associated with a more favorable prognosis and response to platinum-based chemotherapy in ovarian cancer. J. Gynecol. Oncol.25, 221–228. 10.3802/jgo.2014.25.3.221
100
WuM.ShangX.SunY.WuJ.LiuG. (2020). Integrated analysis of lymphocyte infiltration-associated lncRNA for ovarian cancer via TCGA, GTEx and GEO datasets. PeerJ8, e8961. 10.7717/peerj.8961
101
XiaoD.HeJ. (2010). Epithelial mesenchymal transition and lung cancer. J. Thorac. Dis.2, 154–159. 10.3978/j.issn.2072-1439.2010.02.03.7
102
XuJ.FangY.ChenK.LiS.TangS.RenY.et al (2022). Single-cell RNA sequencing reveals the tissue architecture in human high-grade serous ovarian cancer. Clin. Cancer Res.28, 3590–3602. 10.1158/1078-0432.CCR-22-0296
103
ZhangB.ChenF.XuQ.HanL.XuJ.GaoL.et al (2018). Revisiting ovarian cancer microenvironment: a friend or a foe?Protein Cell9, 674–692. 10.1007/s13238-017-0466-7
104
ZhangY.WangX.DuanX.DuT.ChenX. (2023). The synergistic effect of EMT regulators and m6A modification on prognosis-related immunological signatures for ovarian cancer. Sci. Rep.13, 14872. 10.1038/s41598-023-41554-y
105
ZhangZ.TuK.LiuF.LiangM.YuK.WangY.et al (2020). FoxM1 promotes the migration of ovarian cancer cell through KRT5 and KRT7. Gene757, 144947. 10.1016/j.gene.2020.144947
106
ZhouW.SunW.YungM. M. H.DaiS.CaiY.ChenC. W.et al (2018). Autocrine activation of JAK2 by IL-11 promotes platinum drug resistance. Oncogene37, 3981–3997. 10.1038/s41388-018-0238-8
Summary
Keywords
high-grade serous ovarian carcinoma, ascites-derived primary cultures, transcriptomic signature, chemoresistance, epithelial–mesenchymal transition
Citation
Constantinescu D-R, Sorop A, Ghionescu A-V, Lixandru D, Herlea V, Bacalbasa N and Dima SO (2024) EM-transcriptomic signature predicts drug response in advanced stages of high-grade serous ovarian carcinoma based on ascites-derived primary cultures. Front. Pharmacol. 15:1363142. doi: 10.3389/fphar.2024.1363142
Received
29 December 2023
Accepted
13 February 2024
Published
06 March 2024
Volume
15 - 2024
Edited by
Awanish Mishra, National Institute of Pharmaceutical Education and Research, India
Reviewed by
Hemant Kumar, National Institute of Pharmaceutical Education and Research, India
Ranjith Kumar Kankala, Huaqiao University, China
Updates
Copyright
© 2024 Constantinescu, Sorop, Ghionescu, Lixandru, Herlea, Bacalbasa and Dima.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nicolae Bacalbasa, nicolaebacalbasa@gmail.com; Vlad Herlea, herlea2002@gmail.com
† These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.