Skip to main content
Download PDF

Single-cell and bulk RNA sequencing reveal ligands and receptors associated with worse overall survival in serous ovarian cancer

Cell Communication and Signaling volume 20, Article number: 176 (2022) Cite this article

Abstract

Background

Serous ovarian carcinoma is the most frequent histological subgroup of ovarian cancer and the leading cause of death among gynecologic tumors. The tumor microenvironment and cancer-associated fibroblasts (CAFs) have a critical role in the origin and progression of cancer. We comprehensively characterized the crosstalk between CAFs and ovarian cancer cells from malignant fluids to identify specific ligands and receptors mediating intercellular communications and disrupted pathways related to prognosis and therapy response.

Methods

Malignant fluids of serous ovarian cancer, including tumor-derived organoids, CAFs-enriched (eCAFs), and malignant effusion cells (no cultured) paired with normal ovarian tissues, were explored by RNA-sequencing. These data were integrated with single-cell RNA-sequencing data of ascites from ovarian cancer patients. The most relevant ligand and receptor interactions were used to identify differentially expressed genes with prognostic values in ovarian cancer.

Results

CAF ligands and epithelial cancer cell receptors were enriched for PI3K-AKT, focal adhesion, and epithelial-mesenchymal transition signaling pathways. Collagens, MIF, MDK, APP, and laminin were detected as the most significant signaling, and the top ligand-receptor interactions THBS2/THBS3 (CAFs)—CD47 (cancer cells), MDK (CAFs)—NCL/SDC2/SDC4 (cancer cells) as potential therapeutic targets. Interestingly, 34 genes encoding receptors and ligands of the PI3K pathway were associated with the outcome, response to treatment, and overall survival in ovarian cancer. Up-regulated genes from this list consistently predicted a worse overall survival (hazard ratio > 1.0 and log-rank P < 0.05) in two independent validation cohorts.

Conclusions

This study describes critical signaling pathways, ligands, and receptors involved in the communication between CAFs and cancer cells that have prognostic and therapeutic significance in ovarian cancer.

Video abstract

Background

The most common histology of malignancies of the ovary is epithelial cancer, recognized by the World Health Organization as low- and high-grade serous ovarian cancer [1, 2]. High-grade serous ovarian cancer (HGSOC) is the most prevalent and aggressive form of the disease and the leading cause of ovarian cancer deaths [1, 3]. The main therapeutic strategies are debulking surgery and platinum-based chemotherapy [4, 5]. The cytoreductive surgery is performed as primary treatment or as interval debulking surgery in conjunction with neoadjuvant chemotherapy [6]. The asymptomatic nature of ovarian cancer leads to an advanced diagnosis (FIGO stage III-IV) in 70% of cases, which results in a poor prognosis of ~ 25% 5-year survival for patients with stage IIIC-IV disease [7]. Furthermore, more than 30% of patients at advanced stages of the disease have malignant ascites or pleural effusions at diagnosis and are associated with recurrence and metastasis [8,9,10]. These fluids are drained regularly to relieve pain and discomfort and represent an inexpensive and less invasive source of metastatic tumor cells [11].

Communication and interaction between cancer cells and cancer-associated fibroblasts (CAFs) facilitate metastasis development and tumor progression [12, 13]. CAFs play a crucial role in these processes by producing components of the extracellular matrix (ECM) and secreted factors, which also influence angiogenesis, immunity, remodeling of the ECM, chemoresistance, and response to treatment (reviewed in [13, 14]). However, a comprehensive evaluation of the interactions between CAFs and cancer cells from malignant effusion fluids has not yet been determined.

Single-cell RNA sequencing (scRNA-Seq) has provided new insights into cancer biology by revealing the landscape of intercellular interactions and cell communication and improving the knowledge of tumor heterogeneity [15, 16]. In HGSOC, scRNA-Seq has been used to determine the transcriptomic heterogeneity of cancer and stromal cells from ascites and tumor tissues with a high-resolution [17,18,19,20]. Using scRNA-Seq in HGSOC, Kan et al. reported that CAFs induce epithelial-mesenchymal transition (EMT) of tumor cells via TGFβ signaling, with consequent effects on chemoresistance and metastasis [21]. A previous study showed that changes in the ECM promote an adaptive response during the co-evolution of stromal and cancer cells that drive aggressiveness in HGSOC patients [22].

In this study, we analyzed intercellular communication networks in malignant fluids from ovarian cancer patients. The inference and analysis of cell–cell communication and the signature of CAFs were obtained from the scRNA-Seq data [17] and combined with our RNA-Sequencing (RNA-Seq) data. We demonstrated that ligands and receptors of the PI3K-AKT pathway are essential mediators of CAFs communication with cancer cells. A set of these ligands and receptors genes was associated with an unfavorable prognosis. We also described promising targets for therapy.

Materials and methods

Patient samples and data collection

Patients diagnosed with serous ovarian cancer referred to the University Hospital of Southern Denmark, Vejle, DK, between February 2020 and March 2021 to drain ascites or pleural effusion were eligible to participate in the study. All patients signed an informed consent form. Malignant effusion fluids were drained and transferred to 50 mL tubes for further analysis. We also included eight normal ovarian tissue samples obtained from individuals who underwent surgery for other causes than cancer. Each sample was coded, and the personal health information was removed. Clinical and histopathological data were provided in a de-identified manner. Clinical data were obtained from patient records (Additional file 14: Table S1).

Sample handling and cell culture

Malignant effusion fluids (four ascites and four pleural effusions) were immediately processed after drainage. Fluids were centrifuged (1000 rpm/5 min), and the pellet was treated with AKC lysis buffer (Merck, Darmstadt, Germany) for 5 min at room temperature to remove erythrocytes. The pellet was resuspended DMEM/F12, GlutaMAX (Gibco, Thermo Fisher Scientific, Waltham, USA) supplemented with 10% fetal bovine serum (FBS; Biological Industries, Beit HaEmek, Israel) and 1% antibiotic–antimycotic 100x (Gibco, Thermo Fisher Scientific, Waltham, USA). Cell number, live-cell fraction, and cell size distribution were determined using the Countess II Automated Cell Counter (Thermo Fisher, Waltham, USA). Ten mL of media were used to incubate the cells during the first 12 h at 37o C with 5% CO2. After 12 h, cells attached to the flasks were cultured in two-dimensional (2D), and the supernatant was collected and centrifuged (1000 rpm/5 min) for three-dimensional (3D) cell culture. The cells were washed in medium and resuspended on ice in (1:1) cold StemPro™ hESC SFM growth full medium (which is composed by DMEM/F-12, GlutaMAX™ medium bovine serum albumin 25%, StemPro® hESC Supplement) (Thermo Fisher Scientific, Waltham, USA), added with FGF (10 µg/mL) (Thermo Fisher Scientific, Waltham, USA), 2-Mercaptoethanol (Thermo Fisher Scientific, Waltham, USA), Penicillin/Streptomycin (10,000U/ 10 mg/mL), Gentamycin, Amphotericin (2.5 µg/mL each) (Sigma-Aldrich, St. Louis, MO, USA), and Matrigel® Basement Membrane Matrix (1: 1full medium; Corning, New York, USA). We used 75 µL of the full media mixture in each well of a 24-well plate and kept it for 30 min at 37 °C, 5% CO2. After incubation, 1 mL of medium and supplements were added to each well, and the cells were incubated at 37 °C, 5% CO2. The medium was changed every 2–3 days, and 3D cells were harvested using Cell Recovery Solution (Corning, New York, USA), according to the manufacturer's protocol. Two-dimensional cell culture when confluent (5 to 10 days) was harvested with 0.5% trypsin–EDTA (Gibco, Thermo Fisher Scientific, Waltham, USA) and pellets were used for RNA isolation. Tumor-derived organoids (TDO) developed within 1–3 weeks.

Histological characterization

Primary tumor tissues, uncultured cells on day 0 (baseline), and TDO (6 to 15 days) were fixed and processed for histological analysis using automated protocols (Tissue-Teck VIP 6AI Tissue Processor, Sakura Finetek, Japan). Representative sections of primary tumors and TDO were stained with hematoxylin and eosin (H&E) and compared by a board-certified pathologist. The expression of CK7, TP53, PAX8, and calretinin was evaluated by immunohistochemistry assays in primary tumors and their derived organoids using the Benchmark Ultra automated instrument (Ventana Medical Systems, Roche, Tucson, USA) or Dako Omnis (Agilent, Santa Clara, USA) (Additional file 14: Table S2).

RNA extraction and sequencing analysis

The RNeasy mini kit (Qiagen, Valencia, CA, USA) was used to isolate total RNA from all samples (malignant effusion cells or baseline, 2D cells, TDO, and normal ovarian tissues) following the manufacturer's recommendations. RNA purity and quantity (Nanodrop spectrophotometer, Thermo Fisher Scientific, Waltham, USA) and its integrity (RNA screen tape on a 2200 TapeStation, Agilent, Santa Clara, USA) were evaluated before downstream applications.

Total RNA (100 ng) was used to prepare libraries with the Illumina Stranded Total RNA Prep Ligation with Ribo-Zero Plus kit (Illumina, San Diego, California, USA), following the manufacturer's protocol. TruSeq RNA UD indexes (Integrated DNA Technology – IDT, Coralville, Iowa, USA) were added to multiplex the samples, which were paired-end sequenced on the NovaSeq 6000 system (Illumina) using the S4 Reagent Kit v1.5 (300 cycles), according to the supplier's recommendations.

RNA-Seq data quality control was performed using FastQC and MultiQC [23, 24]. The reads were aligned with the Ensembl human genome assembly GRCh38 (release 99) using STAR (v.2.7.6a) [25], and the expression count matrix was generated using HTSeq [26]. The principal component plot and heatmaps were generated using the Galaxy [27] or Morpheus [28] platforms on log-transformed DESeq2-normalized counts [29]. Cluster analyses based on Euclidean distance, variance, or marker selection were also performed using Morpheus [28]. We selected the top 1,000 highly variable genes (highest standard deviation) across all samples to identify markers for malignant effusion cells, 2D cells, TDO, and normal ovarian tissues. The highly expressed genes in 2D cell cultures were inputted into the EnrichR tool [30,31,32] for a gene set enrichment analysis using the WikiPathways 2021, MSigDB Hallmark 2020, and Biological Process 2021 (GO) libraries. The top five enrichment terms (lowest adjusted P-value, Fisher's exact test) for each library were included in a final consensus list. The enrichment terms were represented in heatmaps generated by the web tool Morpheus [28].

Single-cell RNA-sequencing data processing

Droplet-based scRNA-Seq data for eight HGSOC ascites samples were downloaded from Gene Expression Omnibus[33] (GSE146026). The identification and annotation of the 9609 cells to 18 cell clusters was obtained from the original publication [17] and processed using the default parameters of Seurat V3 [34]. Cells were clustered and visualized using the t-distributed stochastic neighbor embedding (tSNE) or violin plots in Seurat. The differentially expressed genes (DEGs, positive and negative cluster markers) were defined with the 'FindAllMarkers' Seurat function. We subset the cells into four clusters containing myofibroblastic and inflammatory CAFs. This subset was then used to identify the markers among the fibroblasts clusters and visualized using the tSNE or violin plots in Seurat. The list of DEGs for each cluster was inputted into the EnrichR tool [30,31,32] for a gene set enrichment analysis using WikiPathways 2021 and MSigDB Hallmark 2020. As described above, we included the top five enrichment terms (lowest adjusted P-value, Fisher's exact test) for these libraries in a final consensus list. The heatmaps for these enrichment terms were generated with the web tool Morpheus [28]. The expression profile of 40 CAFs markers (Top 10 from each CAF cluster) was analysed in our ovarian bulk RNA-Seq data to identify samples enriched with cells expressing CAF marker genes.

Analysis of cell–cell communication

We quantitatively inferred and analyzed cellular communication networks from the HGSOC scRNA-seq data using CellChat [35] (v.1.1.3). To identify the potential ligand-receptor interactions, we loaded the Seurat object into CellChat (https://github.com/sqjin/CellChat/; accessed on 01/10/2021) using the default parameters. This analysis included 1939 human-validated interactions for paracrine/autocrine signaling (61.8%), extracellular matrix receptor interactions (ECM) (21.7%), and cell–cell contact interactions (16.5%). The aggregated cell–cell communication network was calculated by counting the interactions (represented in a circle plot), followed by the identification of cell–cell communication mediated by significant interactions (ligand-receptor pairs represented in bubble plots). We also filtered and explored specific interactions between CAFs and epithelial cancer cells using the webtool Morpheus [28]. Finally, CellChat was used to identify global communication patterns among cell types. The Cophenetic and Silhouette metrics were initially used to determine the number of incoming and outgoing communication patterns between cell groups and pathways. The signals that contribute the most to outgoing or incoming communication signaling (probability score) for each cell group were then identified and visualized using heatmaps and bar graphs.

The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) RNA-sequencing data analysis

The expression levels [log2(normalized counts + 1)] of selected ligands and receptors genes (PI3K-AKT signaling pathway) in serous ovarian cancer (The Cancer Genome Atlas, TCGA, n = 418) and healthy tissues (Genotype-Tissue Expression, GTEx, n = 88) were downloaded from UCSC Xena [36] (https://xenabrowser.net/). The ligands and receptors genes for PI3K-AKT signaling were selected for further analysis based on consistent findings of the communication patterns between CAFs and cancer cells. The DEGs were identified by Welch's t-test and presented in a box plot using the UCSC Xena. Expression profiles [log2(normalized counts + 1)] of ligands and receptors genes of the PI3K-AKT signaling pathway were further used to perform clustering analysis using Morpheus [28].

Overall survival analysis

The Kaplan–Meier plotter (http://kmplot.com/analysis) [37] was used to evaluate the effect of the expression of ligands and receptors genes (microarrays data) on patients’ overall survival. This analysis included two cohorts of patients with serous ovarian cancer, TCGA (n = 557) and GSE9891 (n = 264). The patients were divided into two groups using the autoselect best cutoff, and only JetSet best probes [38] were used for the analysis. These two groups were compared using a Kaplan–Meier survival plot. The hazard ratio with 95% confidence intervals and log-rank test P values were determined.

Results

Clinical characteristics of the serous ovarian cancer patients

Eight patients with serous ovarian cancer (FIGO stage III-IV) [39] who underwent drainage of malignant effusions (four pleural effusions and four ascites) with or without neoadjuvant treatment were included in this study (Additional file 14: Table S1). The median age of the patients was 72.5 (37–83) years old. Seven patients deceased of the disease (last follow-up February 2022). Four patients underwent debulking surgery. Six patients received neoadjuvant chemotherapy, and maintenance therapy was applied to four patients. Two patients presented previous colorectal cancer (case 4) and breast cancer (case 5) treated with surgery and surgery and chemotherapy, respectively. These malignant effusion fluids were processed, and cells were seeded in a culture medium to generate 2D cell culture and TDO. In one patient, two TDO were obtained from distinct collections and days in culture (11 and 35 days), which presented similar transcriptional profiles (Fig. 1A). The organoids derived from ovarian tumors presented different growth rates and formed structures with morphology and cohesiveness that varied from dense, low cohesive, or cystic (Additional file 1: Fig. S1), as described by Maenhoudt et al. [40].

Fig. 1
figure 1

RNA-Sequencing analysis in cells from malignant effusions (baseline), 2D, and tumor-derived organoids (TDO) compared with normal ovarian tissues. A Top 1000 highly variable genes (highest standard deviation) across all samples. Rows and columns grouped by each condition (baseline, 2D, TDO, or normal tissue) were clustered based on the Euclidean distance between normalized counts values. In one patient, two TDO were obtained from distinct collections and days in culture (11 and 35 days), and are represented by an arrow and asterisks, respectively, at the bottom of the heatmap. AS: ascites and PE: pleural effusion. B Two-dimensional (2D) cell culture system presents a set of genes enriched for epithelial-mesenchymal transition (MSigDB Hallmark 2020), collagen fibril and extracellular matrix (ECM) organization (GO – Biological Process), and Focal Adhesion and PI3K-AKT-mTOR-signaling pathways (WikiPathway 2021). C Unsupervised PCA of RNA-Seq data (normalized counts) showing the segregation of ovarian cancer 2D cells and normal tissues from the cluster generated by ovarian cancer-derived organoids and baseline

Full size image

Based on the histological analysis of the primary ovarian cancer samples, six patients were diagnosed with high-grade and two with low-grade serous carcinoma (Additional file 14: Table S1). We compared primary tumors and TDO using hematoxylin–eosin staining and immunohistochemical analyzes. The expression of CK7, PAX8, Calretinin, and TP53 in TDO was similar to their corresponding primary tumors (representative example in Additional file 2: Fig. S2).

RNA-Sequencing reveals that 2D cell culture is enriched for genes involved in the epithelial-mesenchymal transition and extracellular matrix organization

We selected the 1,000 most variable genes among all samples to capture specific biological processes and pathways associated with each condition (Fig. 1A). For most samples, cells from pleural effusions and ascites or from high-grade and low-grade tumors did not cluster closely but were more widely scattered throughout the heatmap (Fig. 1A). We found that the TDO and baseline are clustered in the PCA plot between 2D cells and normal tissues (Fig. 1C). The sequencing data of two TDO obtained from one patient were grouped together with similar expression patterns (Fig. 1A). The 2D cells enriched for a set of genes involved in the epithelial-mesenchymal transition, organization of collagen fibrils and ECM, focal adhesion and PI3K-AKT-mTOR signaling pathways (Fig. 1B, Additional file 14: Table S3). These data suggest that our 2D cells are enriched with ECM component–secreting cells and constitute a cell culture system enriched with CAFs.

Single-cell transcriptomics reveals the heterogeneity of cancer-associated fibroblasts (CAFs) in ovarian cancer ascites

We integrated our bulk RNA-Seq data with the publicly available scRNA-Seq data of malignant ascites from HGSOC patients [17]. First, we followed the single-cell pipeline described by the authors [17] to reanalyze this HGSOC scRNA-Seq dataset and generate the same 18 clusters of cells (Additional file 14: Tables S4 and S5). These clusters include ovarian cancer cells (Ep), myofibroblastic cancer-associated fibroblasts (myCAF), inflammatory cancer-associated fibroblasts (iCAF), macrophages (mac), dendritic cells (DC1), B cells (B), T cells (T), and erythrocytes (ery) (Additional file 3: Fig. S3, Additional file 14: Tables S4 and S5). Next, we confirmed the expression of canonical CAFs markers to define four clusters that include myCAF (myCAF1 and myCAF2) and iCAF (iCAF1 and iCAF2) between all cell types (Fig. 2A and 2B, Additional file 14: Tables S6 and S7). All CAFs clusters expressed ACTA2, while FAP, TGFB1, and IGF1 were expressed in iCAF1, myCAF2, and iCAF2 (Fig. 2B). The markers CXCL12 and IL6 were highly expressed in iCAF1 and iCAF2 but expressed at low levels in myCAF2 (Fig. 2B). These findings confirmed the enrichment of CAFs in our 2D system and were used to estimate the mechanisms of CAFs communication and interactions with cancer cells.

Fig. 2
figure 2

Single-cell transcriptomics re-analysis confirm the heterogeneity of cancer-associated fibroblasts (CAFs) in ovarian cancer malignant ascites. A Single-cell RNA-sequencing (scRNA-seq) data visualization of four CAFs populations using t-distributed stochastic neighbor embedding (tSNE). B Violin plots show the expression levels of canonical CAFs markers in CAFs populations, including myofibroblastic cancer-associated fibroblasts (myCAF1-2) and inflammatory cancer-associated fibroblasts (iCAF1-2). C Heatmap shows the pathway enrichment analysis (EnrichR [30, 31]) for gene sets of cluster-specific CAFs markers. Top 25 pathway terms (lowest adjusted P-value) selected from the WikiPathways 2021 and MSigDB Hallmark 2020 libraries available in EnrichR [30, 31]. D Violin plots show the expression levels of ovarian cancer CAFs markers among the four CAFs populations

Full size image

A gene set enrichment analysis using cell markers was performed to characterize pathways exclusively associated with CAFs subtypes (Fig. 2C). All CAFs were strongly enriched for MYC targets, epithelial-mesenchymal transition pathway, VEGFA-VEGFR2 signaling pathway, glycolysis, fatty acid metabolism, and mTOR1 signaling (Fig. 2C). The iCAF1 cluster featured genes related to hypoxia and TNF-alpha signaling via NF-KB, while the iCAF2 were associated with oxidative phosphorylation (Fig. 2C). myCAF2 showed a specific enrichment of genes related to mTORC1 signaling, proteasome degradation, and E2F targets (Fig. 2C), while myCAF1 presented an enrichment of genes associated with cytoplasmic ribosomal proteins.

The top DEGs capable of distinguishing the four clusters of CAFs were then identified (Fig. 2D, Additional file 14: Table S8). iCAF1 cluster presented reduced expression of RPL9, NME1-NME2, and RPS10, while IGGL5, EEF1G, and RPS26 were expressed exclusively in iCAF2 (Fig. 2D). myCAF1 exhibited reduced expression of the mitochondrial genes MT-CO3, MT-ND1, and MT-CO2, while myCAF2 was marked by increased expression of HBB, CLND4, and IFI44L (Fig. 2D).

Cancer-associated fibroblasts share common outgoing communication patterns for PI3K-AKT signaling pathway and focal adhesion in malignant effusions from ovarian cancer patients

We used CellChat [35] to explore putative cellular interactions and communication (receptor-ligand pairs) between CAFs and other cells in malignant effusion fluids of HGSOC patients evaluated by scRNA-Seq (Additional file 14: Table S9). This analysis included cellular interactions and communication based on secreted signaling (61.8%), extracellular matrix-receptor (21.7%), and cell–cell contact genes (16.5%) (Additional file 4: Fig. S4A). iCAF1, myCAF2, and iCAF2 are cells that contribute the most to outgoing signals (Fig. 3A, Additional file 4: Fig. S4B). Collagens, MIF, MDK, APP, and laminin are significantly involved in cellular interaction and communication for outgoing signals (Fig. 3A).

Fig. 3
figure 3

Cancer-associated fibroblasts (CAFs) share common outgoing communication patterns for the PI3K-AKT signaling pathway and focal adhesion in malignant effusions from ovarian cancer patients. Outgoing patterns for signaling A, cell B, and communication C of multiple cell types in malignant abdominal fluids of ovarian cancer patients obtained with CellChat [35] from scRNA-seq data [17]. The asterisks in C indicate pathways not associated with CAFs and epithelial cells communication (pattern 1, cluster 2). This communication pattern 1 (cluster 2) was used to identify pathways for gene sets of CAFs ligands D and epithelial cell receptors E for these ligands. In both D and E, the focal adhesion and PI3K-AKT-mTOR are among the top altered pathways. The top 10 pathway terms (lowest adjusted P-value) were selected from WikiPathways 2021 and MSigDB Hallmark 2020 libraries available in EnrichR [30, 31]

Full size image

We used a method of analysis of gene expression patterns available on CellChat to explore how cells and signaling pathways globally coordinate to function. First, we identified the correspondence between inferred latent communication patterns with groups of secreting cells to decipher outgoing communication patterns. Three signaling patterns were found (Additional file 5: Fig. S5A), each of them associated with specific cell types that contribute mainly to the outgoing communication: pattern #1 (CAFs), pattern #2 (immune and blood cells), and pattern #3 (epithelial cancer cells) (Fig. 3B). Then, to identify the pathways associated particularly with outgoing communication shared by the four CAFs subtypes, we selected the list of all CAFs ligands (Additional file 14: Table S9) found in the pathways most contributing to communication pattern 1 (cluster 2; Fig. 3C). These CAFs-specific ligands enriched terms related to miRNA targets in ECM and membrane receptors, PI3K-AKT signaling, focal adhesion, epithelial-mesenchymal transition, and focal adhesion-PI3K-AKT-mTOR-signaling pathway (Fig. 3D). Epithelial cancer cell receptors for these CAFs ligands highlight a similar set of epithelial-mesenchymal transition, PI3K-AKT, and focal adhesion-related pathways as essential mediators of crosstalk between CAFs and cancer cells (Fig. 3E).

Cancer-associated fibroblasts receive a specific signaling pattern associated with the differentiation pathway and EGFR tyrosine kinase inhibitor resistance

We also identified signals that most contribute to the incoming signaling to cell populations in malignant ascites (CellChat) [35]. Based on the overexpression of ligands and receptors in all cells, we verified that myCAF2 and iCAF2 are among the leading receivers of most incoming signaling, followed by B cells, cancer cells, and macrophages (Fig. 4A). Collagens, MIF, MDK, APP, and laminin were detected as the most significant incoming signaling molecules (Fig. 4A). We found three incoming signaling patterns (Fig. 4B and Additional file 5: Fig. S5B), in which all CAFs are highly associated with signaling pattern #1 (Fig. 4B). The cells ligands to CAFs (Additional file 14: Table S9) found in the communication pattern #1 (cluster 2; Fig. 4C) were enriched for neovascularization processes, differentiation pathway, EGFR tyrosine kinase inhibitor resistance, and IL-6/JAK/STAT3 signaling (cluster 2, Fig. 4D). CAFs receptors for epithelial cancer cell ligands are consistently involved in the PI3K-AKT and focal adhesion-related pathways (Fig. 4E).

Fig. 4
figure 4

Cancer-associated fibroblasts (CAFs) receive a specific communication pattern from malignant effusion cells associated with the differentiation pathway and EGFR tyrosine kinase inhibitor resistance. Incoming patterns for signaling A, cell B, and communication C of multiple cell types in malignant abdominal fluids of ovarian cancer patients obtained with CellChat [35] from scRNA-seq data [17]. The asterisks in C indicate pathways not associated with CAFs and epithelial cells communication (communication pattern 1, cluster 2). This communication pattern 1 (cluster 2) was used to identify pathways for gene sets of cells ligands D and CAFs receptors for these ligands E. In both D and E, the focal adhesion-PI3K-AKT-mTOR are among the top altered pathways. The top 10 pathway terms (lowest adjusted P-value) were selected from WikiPathways 2021 and MSigDB Hallmark 2020 libraries available in EnrichR [30, 31]

Full size image

Top interactions between cancer-associated fibroblasts and epithelial cells

The high number of ligand and receptor interactions (Additional file 6 Additional files: 78 and 9: Fig. S6A-D) indicated the complexity of the process that extended to and potentially beyond pattern #1 for outgoing and incoming signaling (cluster 2; Fig. 3C to Fig. 4C). Thus, to explore the top interactions involved in ovarian cancer, we first selected the outgoing signaling from CAFs to epithelial cancer cells with a communication probability > 0.10. Among these, we found interactions that may constitute potential targets for CAFs-based therapies, such as THBS2/THBS3 (myCAF2) and CD47 (cancer cells) or between MDK (CAFs) and SDC2/SDC4/NCL (cancer cells) (Fig. 5A).

Fig. 5
figure 5

Top interactions between cancer-associated fibroblasts (CAFs) and epithelial cells have multiple ligands and receptors associated with the PI3K-AKTsignaling pathway and focal adhesion. A Top ligand-receptor interactions (Commun. Probability > 0.10) between CAFs and epithelial cells. Heatmaps show the EnrichR pathway enrichment analysis for gene sets that define all CAFs ligands B and epithelial cell receptors C. The top 10 pathway terms (lowest adjusted P-value) were chosen from WikiPathways 2021 and MSigDB Hallmark 2020 libraries available in EnrichR [30, 31]. D PI3K-AKT ligand-receptor interactions between CAFs and epithelial cells. E Top ligand-receptors interactions (Commun. Probability > 0.10) between epithelial cells and CAFs. Heatmaps showing the EnrichR pathway enrichment analysis for gene sets that define all epithelial cells ligands F and CAFs receptors G. The top 10 pathway terms (lowest adjusted P-value) were chosen from WikiPathways 2021 and MSigDB Hallmark 2020 libraries available in EnrichR [30, 31]. Commun. Prob., communication probability

Full size image

A comprehensive evaluation of all ligands and receptors interactions between CAFs and epithelial cells (Additional file 6 Additional files: 7, 8 and 9: Fig. S6A-D) confirmed the crucial role of components of the focal adhesion-PI3K-AKT-mTOR signaling pathway (Fig. 5B–C). These interactions between CAFs and cancer cells (Additional file 6 Additional files: 78 and 9: Fig. S6A-D) were filtered, and a high-resolution map with the most relevant interactions of the PI3K-AKT pathway revealed the involvement of collagens, fibronectin, vitronectin, laminin, and osteopontin from CAFs with integrins in cancer cells (Fig. 5D).

CAFs also have an essential role in maintaining an inflammatory tumor environment by secreting cytokines, chemokines, and ECM proteins that recruit and activate different immune effector cells [41]. We used CellChat to infer and analyze CAFs-immune cells communication using single-cell data (Additional file 14: Table S9). The bi-directional crosstalk between CAFs and immune cells further confirmed that the PI3K-AKT signaling pathway is a general marker of the CAFs signaling (CAFs ligands and receptors) with other cell types (Additional file 10: Fig. S7).

MIF and MDK are essential regulators of cancer-associated fibroblasts (CAFs) by epithelial cancer cells

We speculate that the diffusion of ligands produced by epithelial cancer cells into the adjacent stroma may also stimulate CAFs. The top outgoing signaling from epithelial cancer cells to CAFs (communication probability > 0.10) showed interactions between MIF (cancer cells) and CD74 + CD44 (CAFs) or MDK (cancer cells) and NCL/SDC2/LRP1 (CAFs) (Fig. 5E). In concordance with previous results herein described, both epithelial cancer cell ligands and CAFs receptors were enriched for PI3K-AKT and focal adhesion-related pathways. These ligands were highly correlated with their receptors (and vice versa) with redundancies for these top-enriched pathways (Fig. 5F–G).

Relevance of ligands and receptors genes associated with the PI3K-AKT signaling pathway in ovarian cancer

We explored the PI3K-AKT signaling pathway using our RNA-Seq data performed in malignant fluid cells from eight HGSOC patients. This analysis included the baseline, 2D cells, and TDO combined with normal ovarian tissues (Additional file 14: Table S10). First, we confirmed that our 2D cells were eCAFs and captured the similar expression profiles found in the scRNA-Seq re-analysis. To this, we analyzed the expression profile of the top ten markers that distinguished each CAFs cluster from the other ascites cells of the scRNA-Seq re-analysis (Additional file 11: Fig. S8). The expression profile of these marker genes did not show differences between pleural effusion and ascites, nor between high- and low-grade samples (Additional file 11: Fig. S8). Next, we analyzed the expression profile of CAFs ligands targeting epithelial cells (Fig. 6A) or CAFs ligands within pattern #1 (from cluster 2, in Fig. 3C) (Additional file 12: Fig. S9), and both showed a clear enrichment of CAFs outgoing signaling. However, compared to normal tissues, epithelial cell receptors were enriched in eCAFs, TDO, and baseline (Fig. 6B). The profile of cell receptors indicated a complex regulatory interplay based on CAFs paracrine and autocrine signaling (Fig. 6B, Additional file 4: Fig. S4B). CAFs seem to act through ligands that stimulate an increase in the expression of the PI3K-AKT genes in epithelial cancer cells, including downstream effector genes such as AKT1 and MTOR (Fig. 6C). These downstream PI3K-AKT effector genes were consistently increased in baseline, TDO, and eCAF (Fig. 6C). Also, many of the PI3K-AKT ligand genes showed increased expression explicitly in eCAFs (Fig. 6D).

Fig. 6
figure 6

Establishment of in vitro system models based on tumor-derived organoids and culture enriched with cancer-associated fibroblasts (eCAFs) to validate the expression profile of ligand and receptor genes associated with the PI3K-AKT signaling pathway in ovarian cancer. Gene expression of CAFs ligands A and epithelial cells receptors (B) for these ligands. Genes were ordered based on a marker selection (signal to noise) that highlights the differences in the expression profile (normalized counts) of eCAFs A or normal tissues B compared to the other three conditions. C Top 50 highly variable genes from the PI3K-AKT signaling pathway across all samples. Genes were ordered based on a marker selection (signal to noise) that highlights the differences in the expression profile (normalized counts) between normal tissue and malignant cells. D Expression profile (normalized counts) of ligands and receptors of the PI3K-AKT signaling pathway across all samples. Rows (for ligand or receptor genes) were ordered based on a marker selection (signal to noise) that highlights the differences in the expression profile (normalized counts) of normal tissues when compared to the other three conditions

Full size image

Increased expression of genes involved in PI3K-AKT signaling between cancer-associated fibroblasts and cancer epithelial cells is associated with worse overall survival in ovarian cancer patients

The expression profile of the ligands and receptor genes of the PI3K-AKT signaling pathway in the CAFs-cancer cell crosstalk revealed genes previously associated with the outcome, response to treatment, or overall survival in ovarian cancer (Fig. 6d). This list of 34 genes includes components of ITGA and ITGB superfamily [42], laminins [43], fibroblast growth factors [44, 45], collagens [22, 46], insulin growth factors [47, 48], EGFR [49, 50], FN1 [22, 51,52,53], EPHA2 [54], TNC [55], SPP1 [55, 56], and VTN [22, 57].

We also evaluated the expression levels of these PI3K-AKT pathway ligands and receptors genes by comparing TCGA ovarian cancer cohort (n = 418) with GTEx normal samples (n = 88). Twenty genes exhibited increased and 14 decreased expression levels in ovarian cancer (Fig. 7A and Additional file 13: Fig. S10). We focused on the 20 up-regulated genes of the PI3K-AKT pathway to test their association with overall survival in two cohorts of patients (TCGA, n = 557 and GSE9891, n = 264) available on the Kaplan–Meier (KM) plotter [37]. Nine genes from the TCGA dataset and 12 genes from the GSE9891 dataset were significantly associated with a worse overall survival rate (HR > 1 and log-rank P value < 0.05) (Fig. 7b-c). Seven of these genes, COL1A1, ITGB5, COL1A2, FGFR2, FN1, IGF1, and IGF2, predicted worse overall survival in both datasets (Fig. 7B–C).

Fig. 7
figure 7

Increased expression of genes involved in PI3K-AKT signaling in cancer-associated fibroblasts and cancer epithelial cells from ovarian cancer patients is associated with worse overall survival. A Expression levels [log2(norm_count + 1)] of ligands and receptors of the PI3K-AKT signaling pathway in ovarian cystadenocarcinoma (n = 419) of The Cancer Genome Atlas (TCGA) compared to normal ovarian tissues (n = 88) tissues of Genotype-Tissue Expression (GTEx). Rows (for ligand or receptor genes) were ordered based on a marker selection (signal to noise) that highlights the differences in the expression profile [log2(norm_count + 1)] of GTEx tissues when compared to the TCGA tissues. Twenty genes exhibited increased expression levels and were selected to test their association with overall survival. B Forest plots for these twenty ligand and receptor genes of the PI3K-AKT signaling pathway in two cohorts of patients with ovarian cancer (TCGA, n = 557 and GSE9891, n = 264) available on the Kaplan–Meier (KM) plotter [37]. The hazard ratio (HR) with 95% confidence intervals (CI) determined by Cox proportional hazards model was used to evaluate the association between gene expression values and overall survival. The statistical significance was determined by a log-rank test. C Representative Kaplan–Meier curves of overall survival for patients with ovarian cancer (TCGA and GSE9891) based on the expression of IGF2, COL1A2, IGF1, and FGFR2

Full size image

Discussion

Increasing evidence indicates that communication and interactions between cancer cells and CAFs are essential determinants of tumor metastasis and progression [12, 13]. The molecular mechanisms involved in these interactions may contribute to the identification of new and more effective therapies or unveil clinically useful prognostic biomarkers for patients with ovarian cancer. We reanalyzed scRNA-seq data from malignant ascites of HGSOC patients [17] to infer and investigate cell–cell interactions and communications. Next, we integrated these gene lists with our RNA-Seq data performed in patient-derived 2D cell culture and tumor organoids. Interestingly, our 2D cells system was enriched with CAFs (eCAFs), as demonstrated by a specific transcriptional signature comprising stromal-related pathways. Our gene expression-based strategy revealed several ligand and receptor interactions between CAFs and cancer epithelial cells associated with the PI3K-AKT signaling pathway and focal adhesion. From the total of 34 ligand and receptor genes identified in these interactions, seven (COL1A1, ITGB5, COL1A2, FGFR2, FN1, IGF1, and IGF2) were consistently up-regulated and predicted worse overall survival in two additional cohorts (TCGA and GSE9891) of ovarian cancer patients.

CAFs exhibit considerable plasticity and heterogeneity [58], and different subpopulations have been reported in ovarian cancer ascites and tumor tissues. Givel et al. [59] described the stromal heterogeneity in HGSOC with four CAFs subpopulations (CAF-S1 − S4). According to the authors, the accumulation of the CAF-S1 subset was associated with an immunosuppressive tumor environment. On the other hand, Kan et al. [21] identified two subpopulations of CAFs (CAF1 and CAF2), where CAF1 was strongly associated with metastasis. Hussain et al. [60] described two CAFs states (FAP-low and FAP-high) and showed that CAFs FAP-high promote proliferation, invasion, and therapy resistance of cancer cells. Izar et al. [17] identified four CAFs sub-populations, with two expressing immune-related genes (complement factors, chemokines, and cytokines). We demonstrated that between these CAFs expressing immune-related genes (iCAF1 and iCAF2), iCAF1 featured genes related to hypoxia and TNF-alpha signaling via NF-KB, while iCAF2 were highly associated with oxidative phosphorylation (Fig. 2C). We also found that myCAF2 showed a specific enrichment of genes related to mTORC1 signaling and proteasome degradation, and E2F targets, while myCAF1 exhibited reduced expression of mitochondrial genes (Fig. 2C). Importantly, these results showed that CAFs present distinct functional states but also demonstrated that these cells share signaling patterns that can be explored for targeted therapies.

Several studies have shown that the activation of CAFs by cancer cells is mediated via secretome components associated with ECM remodeling, cytokine, chemokines, and growth factor-mediated signaling [58, 61,62,63,64]. We found that CAFs from malignant fluids of ovarian cancer receive a specific signaling pattern from cancer cells, mainly composed of ligands associated with neovascularization processes, differentiation pathway, EGFR tyrosine kinase inhibitor resistance, and IL-6/JAK/STAT3 signaling (Fig. 4D). Furthermore, we verified that MIF and MDK are essential regulators of CAFs by epithelial cancer cells targeting the heteromeric complex CD74 + CD44 and NCL/SDC2/LRP1, respectively (Fig. 5E). MIF is known to be produced by ovarian cancer cells in an autocrine manner and may promote colonization of the peritoneum and neovascularization of tumor deposits by other cytokines, chemokines, and growth factors [65]. Furthermore, MIF is secreted in ascites, and its serum levels in patients with ovarian cancer correlate with a poor prognosis [66]. Several therapeutic strategies under development use antibodies to block MIF or CD74 and thus prevent MIF signaling in different tumor types [67]. Therefore, our findings may be valuable not only for helping to explain the role of MIF in ovarian cancer but also for designing targeted therapies that act on specific molecular targets of cancer cells that may activate CAFs.

The release of factors by activated CAFs enormously impacts cancer cells in the crosstalk that occurs in ovarian cancer and is associated with several clinicopathological characteristics and disease outcomes [13, 58]. Among the top CellChat interactions, we identified MDK from CAFs targeting SDC2/SDC4/NCL of cancer cells (Fig. 5A). Thus, MDK may participate in two-way communication between CAFs and ovarian cancer cells and act in an autocrine or paracrine manner. MDK is a growth factor that acts on cancer progression and constitutes a potential therapeutic target [68]. In addition, among these top interactions, we found THBS2/THBS3 secreted by myCAF2 targeting CD47 of cancer cells (Fig. 5A). CD47 is a highly and ubiquitously expressed cell surface protein in ovarian cancer [69] that induces cancer cell growth and predicts poor prognosis [70, 71]. CD47 inhibits macrophage phagocytosis, which contributes to ovarian cancer progression [72, 73]. Consequently, CD47 has been tested in ovarian cancer as a promising CAR-T cell-based therapy [74]. A treatment that combines engineered CAR-T cells targeting CD47 and inhibits secreted THBS2/THBS3 using antibodies may constitute a valuable therapeutic strategy to inhibit ovarian cancer progression.

We show that the interplay between CAFs and epithelial cancer cells has multiple ligand and receptor interactions associated with the PI3K-AKT signaling pathway and focal adhesion (Fig. 5B–C). The activation of this pathway promotes cellular proliferation, migration, and invasion in the ovarian cancer [75]. Our RNA-seq data confirmed the relevance of these ligands and receptors of the PI3K-AKT signaling pathway (Fig. 6). CAFs induce cancer cell proliferation and metastasis by activating the PI3K-AKT-mTOR pathway in lung [76,77,78], colon [79], gastric [80], oral [81], endometrial [82], and anal [83] cancers. The CAF-derived ligand POSTN (Fig. 6A) found in our study was previously associated with chemoresistance in ovarian cancer [50]. Furthermore, the secretion of FNI by mesothelial cell-derived CAFs was also shown to decrease the platinum sensitivity of ovarian cancer cells by inducing the PI3K-AKT pathway [84]. Drugs targeting the PI3K-AKT signaling pathway, such as BEZ235, AZD5363, and NSC777213, have been tested to treat ovarian cancer with promising results [85,86,87]. A phase I trial with PI3-kinase inhibitor BKM120 in combination with PARP inhibitor olaparib is being conducted in patients with high-grade serous ovarian cancer, with evidence of clinical benefits [88]. The pharmacological and genetic inhibition of TTK, which is involved in the PI3K-AKT pathway, decreases the proliferation of ovarian cancer cells and increases their sensitivity to cisplatin by suppressing autophagy [89]. Thus, drugs targeting the crosstalk between CAFs and cancer cells by blocking PI3K-AKT ligands and receptors constitute relevant therapeutic targets.

Considering the intercommunication between CAFs and cancer cells based on components of the PI3K-AKT pathway and its impact on mechanisms that induce tumor aggressiveness or drug resistance, we evaluated whether changes in expression of these ligands and receptors genes are associated with prognosis. Several genes found in our analysis have previously been described as biomarkers of ovarian cancer or potential therapeutic targets, including genes from the ITGA and ITGB superfamily [42], laminins [43], fibroblast growth factors [44, 45, 90], collagens [22, 46], insulin growth factors [47, 48], EGFR [49, 50], FN1 [22, 51,52,53], EPHA2 [54], TNC [55], SPP1 [55, 56], and VTN [22, 57]. We demonstrated that COL1A1, COL1A2, FGFR2, FN1, IGF1, IGF2, and ITGB5 were also altered in two large cohorts of patients (TCGA and GSE9891) and are highly relevant for predicting the prognosis (Fig. 7B). These data demonstrate that CAFs induce changes in cancer cells via components of the PI3K-AKT pathway associated with ECM and adhesion. These results are aligned with findings showing that matrix adhesion is an adaptative response that drives HGSC aggressiveness by co-evolving ECM composition and sensing [22].

COL1A1 plays a central role in carboplatin resistance in ovarian cancer, acting through the ECM-receptor interaction and focal adhesion pathways [46]. In addition to contributing to changes in the PI3K-AKT pathway in ovarian cancer, these interactions may promote peritoneal metastasis by forming ascitic CAFs heterotypic aggregates with tumor cells [91]. Integrins such as ITGB1, ITGB3, ITGB6, ITGA7, and ITGB8 are also involved in cell adhesion and signaling and provide further prognostic information and druggability in HGSOC [42]. We found that the integrin ITGA3 predicted a worse prognosis in the TCGA cohort, ITGB1 in the GSE9891 cohort, while ITGB5 predicted a worse prognosis in both. Despite differences in study design and population, we and Zhu et al. [42] highlighted the prognostic value of ITGA and ITGB superfamily members in serous ovarian cancer.

RNA-Seq data provide unprecedented opportunities to explore new biological questions. Although our strategy has limitations, including the number of cases investigated, the ligand-receptor interactions were validated in two independent ovarian cancer cohorts and in a large number of cell types, strengthening the findings described here. In addition to identifying ovarian cancer biomarkers, we elaborated a detailed map of ligand-receptor interactions in ovarian cancer malignant fluids.

In conclusion, we characterized the interactions and communication between CAFs and cancer cells from ovarian cancer ascites using scRNA-Seq data. This analysis was used to compare and generate ligands and receptors from RNA-Seq data performed in patient-derived eCAFs and tumor organoids obtained from ovarian cancer malignant effusions. PI3K-AKT signaling mediating through ligands and receptors potentially constitutes major players in CAFs interactions and communication with cancer cells. We also verified that a set of ligands and receptor genes of the PI3K-AKT pathway presented potential prognostic value. The expression profile of these drivers at the single-cell level revealed molecular mechanisms and targets that may facilitate the development of therapies focusing on interactions and communication between CAFs and cancer cells.

Availability of Data and Materials

The dataset generated during the current study is available from the corresponding author upon reasonable request. Previously published droplet-based scRNA-Seq data for HGSOC ascites samples have the accession number GSE146026.

Abbreviations

CAFs:

Cancer-associated fibroblasts

eCAFs:

Patient-derived cancer-associated fibroblasts-enriched

HGSOC:

High-grade serous ovarian cancer

ECM:

Extracellular matrix

scRNA-Seq:

Single-cell RNA sequencing

EMT:

Epithelial-mesenchymal transition

RNA-Seq:

RNA-sequencing

TCGA:

The Cancer Genome Atlas

GTEx:

Genotype-Tissue Expression

2D:

Two-dimensional

3D:

Three-dimensional

TDO:

Tumor-derived organoids

H&E:

Hematoxylin and eosin

tSNE:

T-distributed stochastic neighbor embedding

DEGs:

Differentially expressed genes

Ep:

Ovarian cancer cells

myCAF:

Myofibroblastic cancer-associated fibroblasts

iCAF:

Inflammatory cancer-associated fibroblasts

mac:

Macrophages

DC1:

Dendritic cells

B:

B cells

T:

T cells

ery:

Erythrocytes

KM:

Kaplan–Meier

GEO:

Gene Expression Omnibus

References

  1. Sung H, Ferlay J, Siegel RL, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA A Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.

    Article  Google Scholar 

  2. Kurman RJ, International agency for research on cancer, world health organization, eds. WHO classification of tumours of female reproductive organs. 4th ed. International agency for research on cancer; 2014.

  3. Gershenson DM, Sun CC, Lu KH, et al. Clinical behavior of stage II-IV low-grade serous carcinoma of the ovary. Obstet Gynecol. 2006;108(2):361–8. https://doi.org/10.1097/01.AOG.0000227787.24587.d1.

    Article  PubMed  Google Scholar 

  4. du Bois A, Reuss A, Pujade-Lauraine E, Harter P, Ray-Coquard I, Pfisterer J. Role of surgical outcome as prognostic factor in advanced epithelial ovarian cancer: a combined exploratory analysis of 3 prospectively randomized phase 3 multicenter trials: by the Arbeitsgemeinschaft gynaekologische onkologie studiengruppe ovarialkarzinom (AGO-OVAR) and the groupe d’investigateurs nationaux pour les etudes des cancers de l’Ovaire (GINECO). Cancer. 2009;115(6):1234–44. https://doi.org/10.1002/cncr.24149.

    Article  CAS  PubMed  Google Scholar 

  5. Gadducci A, Guarneri V, Peccatori FA, et al. Current strategies for the targeted treatment of high-grade serous epithelial ovarian cancer and relevance of BRCA mutational status. J Ovarian Res. 2019;12(1):9. https://doi.org/10.1186/s13048-019-0484-6.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Vergote I, Tropé CG, Amant F, et al. Neoadjuvant chemotherapy or primary surgery in stage IIIC or IV ovarian cancer. N Engl J Med. 2010;363(10):943–53. https://doi.org/10.1056/NEJMoa0908806.

    Article  CAS  PubMed  Google Scholar 

  7. Lheureux S, Gourley C, Vergote I, Oza AM. Epithelial ovarian cancer. The Lancet. 2019;393(10177):1240–53. https://doi.org/10.1016/S0140-6736(18)32552-2.

    Article  Google Scholar 

  8. Kipps E, Tan DSP, Kaye SB. Meeting the challenge of ascites in ovarian cancer: new avenues for therapy and research. Nat Rev Cancer. 2013;13(4):273–82. https://doi.org/10.1038/nrc3432.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Zeppernick F, Meinhold-Heerlein I. The new FIGO staging system for ovarian, fallopian tube, and primary peritoneal cancer. Arch Gynecol Obstet. 2014;290(5):839–42. https://doi.org/10.1007/s00404-014-3364-8.

    Article  CAS  PubMed  Google Scholar 

  10. Naora H, Montell DJ. Ovarian cancer metastasis: integrating insights from disparate model organisms. Nat Rev Cancer. 2005;5(5):355–66. https://doi.org/10.1038/nrc1611.

    Article  CAS  PubMed  Google Scholar 

  11. Ford CE, Werner B, Hacker NF, Warton K. The untapped potential of ascites in ovarian cancer research and treatment. Br J Cancer. 2020;123(1):9–16. https://doi.org/10.1038/s41416-020-0875-x.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Lau TS, Chan LKY, Wong ECH, et al. A loop of cancer-stroma-cancer interaction promotes peritoneal metastasis of ovarian cancer via TNFα-TGFα-EGFR. Oncogene. 2017;36(25):3576–87. https://doi.org/10.1038/onc.2016.509.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Chen Y, McAndrews KM, Kalluri R. Clinical and therapeutic relevance of cancer-associated fibroblasts. Nat Rev Clin Oncol. 2021;18(12):792–804. https://doi.org/10.1038/s41571-021-00546-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Sahai E, Astsaturov I, Cukierman E, et al. A framework for advancing our understanding of cancer-associated fibroblasts. Nat Rev Cancer. 2020;20(3):174–86. https://doi.org/10.1038/s41568-019-0238-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Baslan T, Hicks J. Unravelling biology and shifting paradigms in cancer with single-cell sequencing. Nat Rev Cancer. 2017;17(9):557–69. https://doi.org/10.1038/nrc.2017.58.

    Article  CAS  PubMed  Google Scholar 

  16. Armingol E, Officer A, Harismendy O, Lewis NE. Deciphering cell–cell interactions and communication from gene expression. Nat Rev Genet. 2021;22(2):71–88. https://doi.org/10.1038/s41576-020-00292-x.

    Article  CAS  PubMed  Google Scholar 

  17. Izar B, Tirosh I, Stover EH, et al. A single-cell landscape of high-grade serous ovarian cancer. Nat Med. 2020;26(8):1271–9. https://doi.org/10.1038/s41591-020-0926-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Shih AJ, Menzin A, Whyte J, et al. Identification of grade and origin specific cell populations in serous epithelial ovarian cancer by single cell RNA-seq. PLoS ONE. 2018;13(11):e0206785. https://doi.org/10.1371/journal.pone.0206785.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Hu Z, Artibani M, Alsaadi A, et al. The repertoire of serous ovarian cancer non-genetic heterogeneity revealed by single-cell sequencing of normal fallopian tube epithelial cells. Cancer Cell. 2020;37(2):226-242.e7. https://doi.org/10.1016/j.ccell.2020.01.003.

    Article  CAS  PubMed  Google Scholar 

  20. Winterhoff BJ, Maile M, Mitra AK, et al. Single cell sequencing reveals heterogeneity within ovarian cancer epithelium and cancer associated stromal cells. Gynecol Oncol. 2017;144(3):598–606. https://doi.org/10.1016/j.ygyno.2017.01.015.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Kan T, Wang W, Ip PP, et al. Single-cell EMT-related transcriptional analysis revealed intra-cluster heterogeneity of tumor cell clusters in epithelial ovarian cancer ascites. Oncogene. 2020;39(21):4227–40. https://doi.org/10.1038/s41388-020-1288-2.

    Article  CAS  PubMed  Google Scholar 

  22. Pietilä EA, Gonzalez-Molina J, Moyano-Galceran L, et al. Co-evolution of matrisome and adaptive adhesion dynamics drives ovarian cancer chemoresistance. Nat Commun. 2021;12(1):3904. https://doi.org/10.1038/s41467-021-24009-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8. https://doi.org/10.1093/bioinformatics/btw354.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Wingett SW, Andrews S. FastQ screen: a tool for multi-genome mapping and quality control. F1000Res. 2018;7:1338. https://doi.org/10.12688/f1000research.15931.2.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Dobin A, Davis CA, Schlesinger F, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21. https://doi.org/10.1093/bioinformatics/bts635.

    Article  CAS  PubMed  Google Scholar 

  26. Anders S, Pyl PT, Huber W. HTSeq–a python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. https://doi.org/10.1093/bioinformatics/btu638.

    Article  CAS  PubMed  Google Scholar 

  27. Afgan E, Baker D, Batut B, et al. The galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018;46(W1):W537–44. https://doi.org/10.1093/nar/gky379.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Starruß J, de Back W, Brusch L, Deutsch A. Morpheus: a user-friendly modeling environment for multiscale and multicellular systems biology. Bioinformatics. 2014;30(9):1331–2. https://doi.org/10.1093/bioinformatics/btt772.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kuleshov MV, Jones MR, Rouillard AD, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90-97. https://doi.org/10.1093/nar/gkw377.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Chen EY, Tan CM, Kou Y, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013;14(1):128. https://doi.org/10.1186/1471-2105-14-128.

    Article  Google Scholar 

  32. Xie Z, Bailey A, Kuleshov MV, et al. Gene set knowledge discovery with enrichr. Curr Protoc. 2021;1(3):90.

    Article  Google Scholar 

  33. Edgar R. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10. https://doi.org/10.1093/nar/30.1.207.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Stuart T, Butler A, Hoffman P, et al. Comprehensive integration of single-cell data. Cell. 2019;177(7):1888-1902.e21. https://doi.org/10.1016/j.cell.2019.05.031.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Jin S, Guerrero-Juarez CF, Zhang L, et al. Inference and analysis of cell-cell communication using cell chat. Nat Commun. 2021;12(1):1088. https://doi.org/10.1038/s41467-021-21246-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Goldman MJ, Craft B, Hastie M, et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 2020;38(6):675–8. https://doi.org/10.1038/s41587-020-0546-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Győrffy B. Survival analysis across the entire transcriptome identifies biomarkers with the highest prognostic power in breast cancer. Comput Struct Biotechnol J. 2021;19:4101–9. https://doi.org/10.1016/j.csbj.2021.07.014.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Li Q, Birkbak NJ, Gyorffy B, Szallasi Z, Eklund AC. Jetset: selecting the optimal microarray probe set to represent a gene. BMC Bioinform. 2011;12(1):474. https://doi.org/10.1186/1471-2105-12-474.

    Article  Google Scholar 

  39. Prat J. FIGO committee on gynecologic oncology staging classification for cancer of the ovary fallopian tube and peritoneum. Int J Gynecol Obstet. 2014;124(1):1–5. https://doi.org/10.1016/j.ijgo.2013.10.001.

    Article  Google Scholar 

  40. Maenhoudt N, Defraye C, Boretto M, et al. Developing organoids from ovarian cancer as experimental and preclinical models. Stem Cell Reports. 2020;14(4):717–29. https://doi.org/10.1016/j.stemcr.2020.03.004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhu T, Chen R, Wang J, Yue H, Lu X, Li J. The prognostic value of ITGA and ITGB superfamily members in patients with high grade serous ovarian cancer. Cancer Cell Int. 2020;20(1):257. https://doi.org/10.1186/s12935-020-01344-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Diao B, Yang P. Comprehensive analysis of the expression and prognosis for laminin genes in ovarian cancer. Pathol Oncol Res. 2021;27:1609855. https://doi.org/10.3389/pore.2021.1609855.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Cole C, Lau S, Backen A, et al. Inhibition of FGFR2 and FGFR1 increases cisplatin sensitivity in ovarian cancer. Cancer Biol Ther. 2010;10(5):495–504. https://doi.org/10.4161/cbt.10.5.12585.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Olbrecht S, Busschaert P, Qian J, et al. High-grade serous tubo-ovarian cancer refined with single-cell RNA sequencing: specific cell subtypes influence survival and determine molecular subtype classification. Genome Med. 2021;13(1):111. https://doi.org/10.1186/s13073-021-00922-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Yang F, Zhao Z, Cai S, et al. Detailed molecular mechanism and potential drugs for COL1A1 in carboplatin-resistant ovarian cancer. Front Oncol. 2021;10: 576565. https://doi.org/10.3389/fonc.2020.576565.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Dong Y, Li J, Han F, et al. High IGF2 expression is associated with poor clinical outcome in human ovarian cancer. Oncol Rep. 2015;34(2):936–42. https://doi.org/10.3892/or.2015.4048.

    Article  CAS  PubMed  Google Scholar 

  47. Bai Y, Ren C, Wang B, et al. LncRNA MAFG-AS1 promotes the malignant phenotype of ovarian cancer by upregulating NFKB1-dependent IGF1. Cancer Gene Ther. 2021. https://doi.org/10.1038/s41417-021-00306-8.

    Article  PubMed  Google Scholar 

  48. Fujiwara S, Terai Y, Kawaguchi H, et al. GPR30 regulates the EGFR-Akt cascade and predicts lower survival in patients with ovarian cancer. J Ovarian Res. 2012;5(1):35. https://doi.org/10.1186/1757-2215-5-35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Wang K, Li D, Sun L. High levels of EGFR expression in tumor stroma are associated with aggressive clinical features in epithelial ovarian cancer. 2016: https://doi.org/10.2147/OTT.S96309

  50. Yoshihara M, Kajiyama H, Yokoi A, et al. Ovarian cancer-associated mesothelial cells induce acquired platinum-resistance in peritoneal metastasis via the FN1/Akt signaling pathway. Int J Cancer. 2020;146(8):2268–80. https://doi.org/10.1002/ijc.32854.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Liang H, Yu M, Yang R, et al. A PTAL-miR-101-FN1 axis promotes EMT and invasion-metastasis in serous ovarian cancer. Molecul Ther Oncolytic. 2020;16:53–62. https://doi.org/10.1016/j.omto.2019.12.002.

    Article  CAS  Google Scholar 

  52. Bao H, Huo Q, Yuan Q, Xu C. Fibronectin 1: a potential biomarker for ovarian cancer. Dis Mark. 2021;2021:1–11. https://doi.org/10.1155/2021/5561651.

    Article  CAS  Google Scholar 

  53. Thaker PH, Deavers M, Celestino J, et al. EphA2 expression is associated with aggressive features in ovarian carcinoma. Clin Cancer Res. 2004;10(15):5145–50. https://doi.org/10.1158/1078-0432.CCR-03-0589.

    Article  CAS  PubMed  Google Scholar 

  54. Didem T, Faruk T, Senem K, et al. Clinical significance of serum tenascin-c levels in epithelial ovarian cancer. Tumor Biol. 2014;35(7):6777–82. https://doi.org/10.1007/s13277-014-1923-z.

    Article  CAS  Google Scholar 

  55. Liu J, Meng H, Li S, et al. Identification of potential biomarkers in association with progression and prognosis in epithelial ovarian cancer by integrated bioinformatics analysis. Front Genet. 2019;10:1031. https://doi.org/10.3389/fgene.2019.01031.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Chang PY, Liao YP, Wang HC, et al. An epigenetic signature of adhesion molecules predicts poor prognosis of ovarian cancer patients. Oncotarget. 2017;8(32):53432–49. https://doi.org/10.18632/oncotarget.18515.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Wu F, Yang J, Liu J, et al. Signaling pathways in cancer-associated fibroblasts and targeted therapy for cancer. Sig Transduct Target Ther. 2021;6(1):218. https://doi.org/10.1038/s41392-021-00641-0.

    Article  CAS  Google Scholar 

  58. Givel AM, Kieffer Y, Scholer-Dahirel A, et al. miR200-regulated CXCL12β promotes fibroblast heterogeneity and immunosuppression in ovarian cancers. Nat Commun. 2018;9(1):1056. https://doi.org/10.1038/s41467-018-03348-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Hussain A, Voisin V, Poon S, et al. Distinct fibroblast functional states drive clinical outcomes in ovarian cancer and are regulated by TCF21. J Exp Med. 2020;217(8): e20191094. https://doi.org/10.1084/jem.20191094.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Wu YH, Huang YF, Chang TH, et al. COL11A1 activates cancer-associated fibroblasts by modulating TGF-β3 through the NF-κB/IGFBP2 axis in ovarian cancer cells. Oncogene. 2021;40(26):4503–19. https://doi.org/10.1038/s41388-021-01865-8.

    Article  CAS  PubMed  Google Scholar 

  61. Yeung TL, Leung CS, Wong KK, et al. TGF-β modulates ovarian cancer invasion by upregulating CAF-derived versican in the tumor microenvironment. Cancer Res. 2013;73(16):5016–28. https://doi.org/10.1158/0008-5472.CAN-13-0023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Cai J, Tang H, Xu L, et al. Fibroblasts in omentum activated by tumor cells promote ovarian cancer growth, adhesion and invasiveness. Carcinogenesis. 2012;33(1):20–9. https://doi.org/10.1093/carcin/bgr230.

    Article  CAS  PubMed  Google Scholar 

  63. Dasari S, Fang Y, Mitra AK. Cancer associated fibroblasts: naughty neighbors that drive ovarian cancer progression. Cancers (Basel). 2018;10(11):E406. https://doi.org/10.3390/cancers10110406.

    Article  CAS  Google Scholar 

  64. Hagemann T, Robinson SC, Thompson RG, Charles K, Kulbe H, Balkwill FR. Ovarian cancer cell-derived migration inhibitory factor enhances tumor growth, progression, and angiogenesis. Mol Cancer Ther. 2007;6(7):1993–2002. https://doi.org/10.1158/1535-7163.MCT-07-0118.

    Article  CAS  PubMed  Google Scholar 

  65. Krockenberger M, Kranke P, Häusler S, et al. Macrophage migration-inhibitory factor levels in serum of patients with ovarian cancer correlates with poor prognosis. Anticancer Res. 2012;32(12):5233–8.

    CAS  PubMed  Google Scholar 

  66. Penticuff JC, Woolbright BL, Sielecki TM, Weir SJ, Taylor JA. MIF family proteins in genitourinary cancer: tumorigenic roles and therapeutic potential. Nat Rev Urol. 2019;16(5):318–28. https://doi.org/10.1038/s41585-019-0171-9.

    Article  PubMed  Google Scholar 

  67. Filippou PS, Karagiannis GS, Constantinidou A. Midkine (MDK) growth factor: a key player in cancer progression and a promising therapeutic target. Oncogene. 2020;39(10):2040–54. https://doi.org/10.1038/s41388-019-1124-8.

    Article  CAS  PubMed  Google Scholar 

  68. Brightwell RM, Grzankowski KS, Lele S, et al. The CD47 “don’t eat me signal” is highly expressed in human ovarian cancer. Gynecol Oncol. 2016;143(2):393–7. https://doi.org/10.1016/j.ygyno.2016.08.325.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Wang CL, Lin MJ, Hsu CY, et al. CD47 promotes cell growth and motility in epithelial ovarian cancer. Biomed Pharmacother. 2019;119: 109105. https://doi.org/10.1016/j.biopha.2019.109105.

    Article  CAS  PubMed  Google Scholar 

  70. Li Y, Lu S, Xu Y, et al. Overexpression of CD47 predicts poor prognosis and promotes cancer cell invasion in high-grade serous ovarian carcinoma. Am J Transl Res. 2017;9(6):2901–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  71. Shimizu A, Sawada K, Kobayashi M, et al. Exosomal CD47 plays an essential role in immune evasion in ovarian cancer. Mol Cancer Res. 2021;19(9):1583–95. https://doi.org/10.1158/1541-7786.MCR-20-0956.

    Article  CAS  PubMed  Google Scholar 

  72. Liu R, Wei H, Gao P, et al. CD47 promotes ovarian cancer progression by inhibiting macrophage phagocytosis. Oncotarget. 2017;8(24):39021–32. https://doi.org/10.18632/oncotarget.16547.

    Article  PubMed  PubMed Central  Google Scholar 

  73. Shu R, Evtimov VJ, Hammett MV, et al. Engineered CAR-T cells targeting TAG-72 and CD47 in ovarian cancer. Molecul Ther Oncolytics. 2021;20:325–41. https://doi.org/10.1016/j.omto.2021.01.002.

    Article  CAS  Google Scholar 

  74. Meng Q, Xia C, Fang J, Rojanasakul Y, Jiang BH. Role of PI3K and AKT specific isoforms in ovarian cancer cell migration, invasion and proliferation through the p70S6K1 pathway. Cell Signal. 2006;18(12):2262–71. https://doi.org/10.1016/j.cellsig.2006.05.019.

    Article  CAS  PubMed  Google Scholar 

  75. Li H, Zhang Q, Wu Q, et al. Interleukin-22 secreted by cancer-associated fibroblasts regulates the proliferation and metastasis of lung cancer cells via the PI3K-Akt-mTOR signaling pathway. Am J Transl Res. 2019;11(7):4077–88.

    CAS  PubMed  PubMed Central  Google Scholar 

  76. Zhou Z, Zhou Q, Wu X, et al. VCAM-1 secreted from cancer-associated fibroblasts enhances the growth and invasion of lung cancer cells through AKT and MAPK signaling. Cancer Lett. 2020;473:62–73. https://doi.org/10.1016/j.canlet.2019.12.039.

    Article  CAS  PubMed  Google Scholar 

  77. Ying L, Zhu Z, Xu Z, et al. Cancer associated fibroblast-derived hepatocyte growth factor inhibits the paclitaxel-induced apoptosis of lung cancer A549 cells by up-regulating the PI3K/Akt and GRP78 signaling on a microfluidic platform. PLoS ONE. 2015;10(6): e0129593. https://doi.org/10.1371/journal.pone.0129593.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Ma J, Sun X, Wang Y, Chen B, Qian L, Wang Y. Fibroblast-derived CXCL12 regulates PTEN expression and is associated with the proliferation and invasion of colon cancer cells via PI3k/Akt signaling. Cell Commun Signal. 2019;17(1):119. https://doi.org/10.1186/s12964-019-0432-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Pang T, Yin X, Luo T, et al. Cancer-associated fibroblasts promote malignancy of gastric cancer cells via nodal signalling. Cell Biochem Funct. 2020;38(1):4–11. https://doi.org/10.1002/cbf.3446.

    Article  CAS  PubMed  Google Scholar 

  80. Li Y yin, Tao Y wei, Gao S, et al. Cancer-associated fibroblasts contribute to oral cancer cells proliferation and metastasis via exosome-mediated paracrine miR-34a-5p. EBioMedicine. 2018;36:209–220 https://doi.org/10.1016/j.ebiom.2018.09.006

  81. Subramaniam KS, Tham ST, Mohamed Z, Woo YL, Mat Adenan NA, Chung I. Cancer-associated fibroblasts promote proliferation of endometrial cancer cells. PLoS ONE. 2013;8(7): e68923. https://doi.org/10.1371/journal.pone.0068923.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Cacheux W, Lièvre A, Richon S, et al. Interaction between IGF2-PI3K axis and cancer-associated-fibroblasts promotes anal squamous carcinogenesis. Int J Cancer. 2019;145(7):1852–9. https://doi.org/10.1002/ijc.32178.

    Article  CAS  PubMed  Google Scholar 

  83. Chu L, Wang F, Zhang W, Li HF, Xu J, Tong XW. Periostin secreted by carcinoma-associated fibroblasts promotes ovarian cancer cell platinum resistance through the PI3K/Akt signaling pathway. Technol Cancer Res Treat. 2020;19:1533033820977535. https://doi.org/10.1177/1533033820977535.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Deng J, Bai X, Feng X, et al. Inhibition of PI3K/Akt/mTOR signaling pathway alleviates ovarian cancer chemoresistance through reversing epithelial-mesenchymal transition and decreasing cancer stem cell marker expression. BMC Cancer. 2019;19(1):618. https://doi.org/10.1186/s12885-019-5824-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Lawal B, Lo WC, Mokgautsi N, et al. A preclinical report of a cobimetinib-inspired novel anticancer small-molecule scaffold of isoflavones, NSC777213, for targeting PI3K/AKT/mTOR/MEK in multiple cancers. Am J Cancer Res. 2021;11(6):2590–617.

    CAS  PubMed  PubMed Central  Google Scholar 

  86. Fabi F, Adam P, Parent S, Tardif L, Cadrin M, Asselin E. Pharmacologic inhibition of Akt in combination with chemotherapeutic agents effectively induces apoptosis in ovarian and endometrial cancer cell lines. Mol Oncol. 2021;15(8):2106–19. https://doi.org/10.1002/1878-0261.12888.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Matulonis UA, Wulf GM, Barry WT, et al. Phase I dose escalation study of the PI3kinase pathway inhibitor BKM120 and the oral poly (ADP ribose) polymerase (PARP) inhibitor olaparib for the treatment of high-grade serous ovarian and breast cancer. Ann Oncol. 2017;28(3):512–8. https://doi.org/10.1093/annonc/mdw672.

    Article  CAS  PubMed  Google Scholar 

  88. Qi G, Ma H, Li Y, Peng J, Chen J, Kong B. TTK inhibition increases cisplatin sensitivity in high-grade serous ovarian carcinoma through the mTOR/autophagy pathway. Cell Death Dis. 2021;12(12):1135. https://doi.org/10.1038/s41419-021-04429-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Madsen CV, Steffensen KD, Olsen DA, et al. Serum platelet-derived growth factor and fibroblast growth factor in patients with benign and malignant ovarian tumors. Anticancer Res. 2012;32(9):3817–25.

    PubMed  Google Scholar 

  90. Han Q, Huang B, Huang Z, et al. Tumor cell-fibroblast heterotypic aggregates in malignant ascites of patients with ovarian cancer. Int J Mol Med. 2019;44(6):2245–55. https://doi.org/10.3892/ijmm.2019.4361.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Salgado-Albarrán M, Späth J, González-Barrios R, Baumbach J, Soto-Reyes E. CTCFL regulates the PI3K-Akt pathway and it is a target for personalized ovarian cancer therapy. npj Syst Biol Appl. 2022;8(1): 5 https://doi.org/10.1038/s41540-022-00214-z

Download references

Acknowledgements

The results shown here are in part based upon data generated by the Genotype-Tissue Expression project (GTEx) (https://gtexportal.org/) and by the TCGA Research Network (http://cancergenome.nih.gov/). The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health and NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS.

Funding

This study was supported by grants from the Region of Southern Denmark Research Fund (J.nr.: 20/14277, Efond 519), Denmark. The study was developed in collaboration with the Brazilian Federal Agency for Support and Evaluation of Graduate Education (CAPES), in the scope of the Program CAPES-PrInt (process number 88887.310463/2018–00).

Author information

Author notes

  1. Robson Francisco Carvalho and Luisa Matos do Canto these authors contributed equally

Authors and Affiliations

  1. Department of Clinical Genetics, University Hospital of Southern Denmark, 7100, Vejle, Denmark

    Robson Francisco Carvalho, Luisa Matos do Canto, Cecilie Abildgaard, Mads Malik Aagaard & Silvia Regina Rogatto

  2. Institute of Regional Health Research, University of Southern Denmark, 5230, Odense, Denmark

    Robson Francisco Carvalho, Luisa Matos do Canto, Cecilie Abildgaard & Silvia Regina Rogatto

  3. Department of Structural and Functional Biology - Institute of Bioscience, São Paulo State University (UNESP), 18.618-689, Botucatu, Brazil

    Robson Francisco Carvalho

  4. Department of Oncology, University Hospital of Southern Denmark, 7100, Vejle, Denmark

    Monica Søgaard Tronhjem, Lars Henrik Jensen & Karina Dahl Steffensen

  5. Danish Colorectal Cancer Center South, 7100, Vejle, Denmark

    Monica Søgaard Tronhjem, Lars Henrik Jensen, Karina Dahl Steffensen & Silvia Regina Rogatto

  6. Department of Clinical Pathology, University Hospital of Southern Denmark, 7100, Vejle, Denmark

    Marianne Waldstrøm

Authors
  1. Robson Francisco Carvalho
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Luisa Matos do Canto
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Cecilie Abildgaard
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Mads Malik Aagaard
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Monica Søgaard Tronhjem
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Marianne Waldstrøm
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Lars Henrik Jensen
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Karina Dahl Steffensen
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Silvia Regina Rogatto
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization and design: R.F.C. and S.R.R.; Data curation: M.M.A, R.F.C., and M.M.C.; Data acquisition: L.M.C., C.A., M.S.T., and M.A.A.; Funding acquisition: S.R.R., L.H.J, M.W., and K.D.S.; Investigation: L.M.C., C.A., M.W., and R.F.C.; Methodology: L.M.C., C.A., and R.F.C.; Supervision and resources: S.R.R.; Software: R.F.C. and M.M.A; Writing: original draft: R.F.C., L.M.C., S.R.R.; Writing: review and editing. All authors. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Robson Francisco Carvalho or Silvia Regina Rogatto.

Ethics declarations

Ethics approval and consent to participate

The Ethics Committee of the Region of Southern Denmark (S-20190102) and The Danish Data Protection Agency (Case No. 19/37168) approved the study. All patients signed a written informed consent before sample collection.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Tumor-derived organoids (TDO) from malignant effusions of serous ovarian cancer patients. Representative examples of TDO from cases 8 (A), 4 (B), and 1 (C) show differences in the growth rate and morphology. (D) Representative images of TDO morphology (dense, low cohesive, and cystic, as described by Maenhoudt et al. [40]) of individual TDO from different patients (cases 4, 7, and 6, respectively) (40x).

Additional file 2: Fig. S2.

(A) High-grade serous carcinoma showing strong and diffuse immunohistochemical expression of CK7, TP53, and PAX8. (B) Section from tumor-derived organoid day 11. The atypical cells show strong and diffuse immunohistochemical expression of CK7, TP53, and PAX, similarly to the primary tumor shown in (a). (C) Low-grade serous ovarian cancer shows a strong and diffuse immunohistochemical expression of CK7 and PAX8 and absent reaction for calretinin. (D) Tumor-derived organoids show atypical cells with strong and diffuse immunohistochemical expression of CK7 and PAX and an absence of reaction for calretinin, similarly to the primary low-grade serous ovarian cancer shown in (C).

Additional file 3: Fig. S3.

Re-analysis based on single-cell RNA-seq data from malignant ascites of eight patients with advanced high-grade serous ovarian cancer described by Izar et al. [17]. The t-SNE shows 9.609 cells analyzed in 18 clusters (colors) that include ovarian cancer cells (Ep1-5), myofibroblastic cancer-associated fibroblasts (myCAF1-2), inflammatory cancer-associated fibroblasts (iCAF1-2), macrophages (mac1-4), dendritic cells (DC1-2), B cells (B), T cells (T), and erythrocytes (ery).

Additional file 4: Fig. S4. (A)

Literature-supported ligand-receptor interactions in humans available at the CellChatDB database (http://www.cellchat.org/) [35]. The 1,939 validated interactions included paracrine/autocrine signaling (61.8%), extracellular matrix (ECM) receptor interactions (21.7%) and cell-cell contact interactions (16.5%). (B) Circos plot showing the CellChatDB inferred cell-cell communication network based on single-cell RNA-Seq data across ovarian cancer cells (Ep1-5), myofibroblastic cancer-associated fibroblasts (myCAF1-2), inflammatory cancer-associated fibroblasts (iCAF1-2), macrophages (mac1-4), dendritic cells (DC1-2), B cells (B), T cells (T), and erythrocytes (ery) from ovarian cancer patients.

Additional file 5: Fig. S5.

Cophenetic and Silhouette metrics were used in CellChatDB [35] to identify the number of outgoing (A) and incoming (B) communication patterns that cell groups and pathways coordinate to function. Both Cophenetic and Silhouette values suggest three patterns, as indicated by the first sudden drop of the measured score in this number of patterns.

Additional file 6: Fig. S6.

Ligand-receptor interactions between cancer-associated fibroblasts and epithelial cells inferred by CellChatDB [35] using single-cell RNA-seq data [17] from malignant ascites of ovarian cancer patients. (A) Inflammatory cancer-associated fibroblasts 1 (iCAF1)-epithelial cells.

Additional file 7: Fig. S6. (B)

Inflammatory cancer-associated fibroblasts 2 (iCAF2)-epithelial cells.

Additional file 8: Fig. S6. (C)

Myofibroblastic cancer-associated fibroblasts 1 (myCAF1)-epithelial cells.

Additional file 9: Fig. S6. (D)

Myofibrobastic cancer-associated fibroblast 2 (myCAF2)-epithelial cells. 

Additional file 10: Fig. S7.

Pathways associated with the signaling (ligands and receptors) from CAFs to immune cells (a and b) and from immune cells to CAFs (d and e). The top 10 pathway terms (lowest adjusted P-value) were selected from WikiPathways 2021 and MSigDB Hallmark 2020 libraries available in EnrichR [30, 31].

Additional file 11: Fig. S8.

Expression profile of 40 cancer-associated fibroblast markers (top 10 for each cluster from scRNA-Seq data) in ovarian bulk RNA-Seq. Genes were ordered based on a marker selection (signal to noise) that highlights the differences in the expression profile (normalized counts) between ovarian cancer-associated fibroblasts-enriched culture (eCAFs) and the other three conditions (Normal tissue, Normal; tumor-derived organoids, TDO; and baseline). (A) grade; (B) malignant effusions or normal tissue; and (C) experimental conditions.

Additional file 12: Fig. S9.

Expression profiling of CAFs ligands (communication pattern 1) in ovarian bulk RNA-Seq. Genes were ordered based on a marker selection (signal to noise) that highlights the differences in the expression profile (normalized counts) between ovarian cancer-associated fibroblasts-enriched culture (eCAFs) and the other three conditions (Normal tissue, Normal; tumor-derived organoids, TDO; and baseline).

Additional file 13: Fig. S10.

Box plot showing the expression levels [log2(norm_count+1)] of ligands and receptors of the PI3K-AKT signaling pathway in ovarian cystadenocarcinoma (n = 419) of The Cancer Genome Atlas (TCGA) compared to normal ovarian tissues (n = 88) tissues of Genotype-Tissue Expression (GTEx).

Additional file 14:

Supplementary Tables (S1-S10).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Carvalho, R.F., do Canto, L.M., Abildgaard, C. et al. Single-cell and bulk RNA sequencing reveal ligands and receptors associated with worse overall survival in serous ovarian cancer. Cell Commun Signal 20, 176 (2022). https://doi.org/10.1186/s12964-022-00991-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12964-022-00991-4

Keywords

Cell Communication and Signaling

ISSN: 1478-811X