Exploring the cellular and molecular differences between ovarian clear cell carcinoma and high-grade serous carcinoma using single-cell RNA sequencing and GEO gene expression signatures
Cell & Bioscience volume 13, Article number: 139 (2023)
The two most prevalent subtypes of epithelial ovarian carcinoma (EOC) are ovarian clear cell carcinoma (OCCC) and high-grade serous ovarian carcinoma (HGSC). Patients with OCCC have a poor prognosis than those with HGSC due to chemoresistance, implying the need for novel treatment target. In this study, we applied single-cell RNA sequencing (scRNA-seq) together with bulk RNA-seq data from the GEO (Gene Expression Omnibus) database (the GSE189553 dataset) to characterize and compare tumor heterogeneity and cell-level evolution between OCCC and HGSC samples. To begin, we found that the smaller proportion of an epithelial OCCC cell subset in the G2/M phase might explain OCCC chemoresistance. Second, we identified a possible pathogenic OCCC epithelial cell subcluster that overexpresses LEFTY1. Third, novel biomarkers separating OCCC from HGSC were discovered and subsequently validated on a wide scale using immunohistochemistry. Amine oxidase copper containing 1 (AOC1) was preferentially expressed in OCCC over HGSC, while S100 calcium-binding protein A2 (S100A2) was detected less frequently in OCCC than in HGSC. In addition, we discovered that metabolic pathways were enriched in the epithelial compartment of the OCCC samples. In vitro experiments verified that inhibition of oxidative phosphorylation or glycolysis pathways exerted direct antitumor effects on both OCCC and HGSC cells, while targeting glutamine metabolism or ferroptosis greatly attenuated chemosensitivity only in OCCC cells. Finally, to determine whether there were any variations in immune cell subsets between OCCC and HGSC, data from scRNA-seq and mass cytometry were pooled for analysis. In summary, our work provides the first holistic insights into the cellular and molecular distinctions between OCCC and HGSC and is a valuable source for discovering new targets to leverage in clinical treatments to improve the poor prognosis of patients with OCCC.
Epithelial ovarian carcinoma (EOC) is the predominantly lethal gynecological cancer and ranks seventh among the most commonly diagnosed cancer in women worldwide . High-grade serous ovarian cancer (HGSC) is the most common EOC and accounts for as many as 60% of EOC cases. Ovarian clear cell carcinoma (OCCC) is the second most common EOC, accounting for 5–11% of EOC cases [2, 3].
The demographic, morphological, and molecular features between OCCC and HGSC differ. OCCC is more prevalent in Asians than in people of other races and is closely related to endometriosis . OCCC displays a tubulocystic architecture and clear cells, while HGSC is characterized by solid and papillary structures. OCCC is characterized by the overexpression of HNF1B [5, 6] and is typically ER- and PR-negative. However, most HGSC cells harbor TP53 and WT-1 mutations . In addition, frequent ER expression and low PR expression increase the likelihood that HGSC is hormone-dependent. HGSC and estrogen or progesterone have been linked in a number of studies [8,9,10].
Surgical procedure with subsequent chemotherapy is the conventional treatment for OCCC and HGSC. However, OCCC leads to a worse prognosis than HSGC due to OCCC chemoresistance to conventional platinum . Therefore, exploring the cellular and molecular differences between OCCC and HGSC will help in developing alternative therapies.
Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool to explore intratumor heterogeneity and evolution at the cell level. It can be used to help interpret the pathogenic mechanism and offers the possibility of personalized treatment. In the present study, we applied scRNA-seq and performed a combined analysis with bulk RNA-seq data from the GEO database to systematically characterize the cellular and molecular differences between OCCC and HGSC (Figure S1). Our study offers a valuable resource for exploring new therapeutic targets of OCCC or both clinically.
A higher proportion of LEFTY1+ and a lower proportion of the G2/M phase epithelial subset in OCCC compared with HGSC
We performed 3’ scRNA-seq with samples collected during surgery from five OCCC and five HGSC patients. Detailed clinical information was provided in Table S1. Representative images of H&E and IHC stained samples were presented in Figure S2A. After filtering low-quality cells and inferred doublets, we integrated all ten samples into a gene expression matrix. As a result, we obtained a total of 101,672 cells with 2000 genes on average in each cell. We performed unsupervised clustering analysis before performing a broad comparison of different cell types in the tumor tissues obtained from OCCC and HGSC patients. The cells were clustered into three main populations: epithelial, immune, and stromal cells (Figure S2B-C).
We then performed an unsupervised sub-clustering analysis and identified eight epithelial subpopulations (Fig. 1A). The large-scale copy number variation (CNV) analysis of malignant cells showed fewer CNVs across the whole genome in all subsets of OCCC cells than in the HGSC cell subsets (Figure S3). The proportion analysis showed that more cells were classified into Cluster 4 and fewer cells were classified into Cluster 7 in OCCC compared with HGSC (Fig. 1B).
The identification of feature markers showed that Cluster 4 expressed MUC1, FGF7, and left-right determination factor 1 (LEFTY1), and Cluster 7 expressed cell cycle-related genes (Fig. 1C). An analysis of the proportion of cells in each cell cycle phase reflected that Cluster 5 cells were arrested mostly in the S phase, while Cluster 7 cells were mainly arrested in the G2/M phase (Fig. 1D). Investigation into the regulon activity using SCENIC showed that the proliferation-related transcription factor HMGB1 was activated in Cluster 5 and Cluster 7 cells but was inactive in Cluster 4 (Fig. 1E).
To explore the development of epithelial cell subsets, we performed a pseudotime trajectory analysis using Monocle 3 based on expression and transition profiles (Fig. 1F). Interestingly, we observed that Cluster 5 cells (S phase) bifurcated into two branches leading to other branches, namely, directly to Cluster 7 or Cluster 1, when in an intermediate state to other clusters (Fig. 1F). Pathway enrichment analyses of each epithelial subtype showed that Cluster 4 was enriched with “Hedgehog signaling”, “Wnt-beta catenin signaling”, “TGF-beta signaling”, “Notch signaling” and “angiogenesis” (Fig. 1G). Taken together, based on the lower proportion of OCCC cells in Cluster 7 and higher proportions of OCCC cells in Cluster 4 compared to the proportion of HGSC cells in these clusters, as described above, we speculated that Cluster 5 OCCC cells might show an increased capacity to develop into Cluster 4 cells but a decreased capacity to develop into Cluster 7 cells compared with the HGSC cells in these clusters.
Identification and validation of distinct molecular markers in OCCC and HGSC cells discovered by combining scRNA-seq and bulk RNA-seq data from the GEO database
A differentially expressed gene (DEG) analysis of epithelial cell subpopulations between OCCC and HGSC patients showed that the expression levels of amine oxidase copper containing 1 (AOC1), left-right determination factor 1 (LEFTY1), and glutathione peroxidase 3 (GPX3) were significantly higher in OCCC samples than in HGSC samples, while the opposite was true for S100 calcium-binding protein A2 (S100A2), WAP four-disulfide core domain 2 (WFDC2), and cellular retinoic acid-binding protein 2 (CRABP2) (Fig. 2A). Further, RT‒PCR and immunoblotting supported these observations (Fig. 2B-C, Figure S4). The DEG analysis with the GSE189553 dataset validated that OCCC samples showed higher expression of AOC1, LEFTY1, and GPX3, and lower expression of CRABP2 compared to HGSC samples (Fig. 2D).
We constructed two tissue microarrays (TMAs) consisting of 128 OCCC samples and 81 HGSC samples. Representative images were illustrated in Fig. 2E and Figure S5. The IHC validation cohort confirmed that AOC1, LEFTY1, and GPX3 were more expressed in the OCCC samples, while WFDC2, CRABP2, and S100A2 were more expressed in the HGSC samples (Table 1). We analyzed the association of clinical features and IHC results with progression-free survival (PFS). The high expression of GPX3 was significantly associated with a longer PFS than low expression in OCCC samples (Fig. 2F). As shown in Table S2, in the univariate analysis, early stage, no residual tumor after primary cytoreductive surgery and high expression of GPX3 were related to a better prognosis for OCCC patients (p = 0.010, 0.013, 0.021, respectively). Early stage (p = 0.014) and expression of GPX3 (p = 0.029) were still correlated with a longer PFS in OCCC patients in the multivariate analysis. For the HGSC group, residual tumor and chemoresistance were significantly related to prognostic data (p = 0.015 and < 0.001, respectively), while chemoresistance was a susceptibility factor for recurrence (p = 0.001). Moreover, a significantly lower proportion of GPX3 expression was found in the platinum-resistant (PR) subgroup than in the platinum-sensitive (PS) subgroup of OCCC patients (Table 2).
Metabolism pathways were activated and might be promising targets for OCCC treatment
Signaling pathway enrichment of differentially expressed genes (DEGs) in the GSE189553 dataset showed that OCCC cells were enriched with “metabolic pathways” (Fig. 3A). Accordingly, signaling pathway enrichment with DEGs in Cluster 5 (dividing cells, S phase) between the OCCC and HGSC populations as distinguished with single-cell RNA seq data revealed that OCCC patient samples were enriched with genes related to “oxidative phosphorylation”, “glutathione metabolism”, “glycolysis/gluconeogenesis” and “ferroptosis” pathways (Fig. 3B).
We further treated the TOV21G cell line (an OCCC cell line) and OVCAR3 cell line (an HGSC cell line) with CB-839, liproxstatin-1, metformin, and 2-deoxy-D-glucose to explore the potential roles of glutamine metabolism, ferroptosis, oxidative phosphorylation, and the glycolysis pathway in OCCC and HGSC cells. The viability of both cell lines was inhibited by metformin or 2-deoxy-D-glucose treatment, which indicated aberrant glucose metabolism in both OCCC and HGSC cells (Fig. 3C). Cell proliferation via EdU assay and apoptosis via flow cytometry revealed that metformin inhibited cell proliferation and promoted cell apoptosis in TOV21G and OVCAR3 cells (Fig. 3D-F). 2-Deoxy-D-glucose functioned through different mechanisms to reduce cell viability. It induced anti-proliferation effects against TOV21G cells and promoted apoptosis in OVCAR3 (Fig. 3D-F).
Interestingly, we found that the inhibition of the glutathione metabolism pathway or ferroptosis by CB-839 or lip-1 significantly increased the viability of the TOV21G cell line treated with cisplatin (Fig. 4A-B). Only a mild effect was observed in the OVCAR3 cell line (Fig. 4A-B). However, this inhibitory effect was not associated with reduced cisplatin-induced apoptosis (Fig. 4C-D). These results indicated that metformin and 2-deoxy-D-glucose show good potential as therapeutic drugs. Targeting glutamine metabolism or ferroptosis greatly attenuated chemosensitivity only in OCCC cells.
Heterogeneity of immune cells in patients with OCCC or HGSC
Next, we compared the extent of tumor-infiltrating immune cells between OCCC and HGSC patients. We used ImmunCellAI to estimate and compare the immune cell abundance in the GSE189553 dataset. We found higher percentages of macrophages and CD8+ T cells in OCCC patients than in HGSC patients (Fig. 5A). In addition, an analysis of T-cell subsets revealed lower percentages of Th17 cells and naïve CD8+ T (CD8_Naive) cells in OCCC patients than in HGSC patients (Fig. 5A).
We produced subsets of immune cells from scRNA-seq data and reidentified ten immune cell populations, including CD8 T cells, CD4 T cells, NK cells, B cells, plasma B cells, dividing cells, DCs, macrophages, neutrophils and pDCs, on the basis of canonical marker levels, (Fig. 5B, Figure S6A). We found that T cells were predominantly enriched in tumor-infiltrating immune cells of HGSC and OCCC patients; however, the proportions of immune cell subsets were similar in OCCC samples compared to those in HGSC samples (Fig. 5C). We further examined the cellular subtypes in the T/NK subsets and identified eight subpopulations: CD8 CTLs, GZMKhi CD8 T cells, IFN-act CD8 T cells, NK cells, NKT/γδ T cells, GD40LG+ CD4 T cells, CXCL13hi CD4 T cells, and Tregs (Figure S6B-C). No significant differences were found in the proportions of these eight T/NK subpopulations between OCCC and HGSC samples (Fig. 5D).
We performed mass cytometry and identified nine immune cell populations, namely, in CD45+ (tumor-infiltrating lymphocytes) TILs, including B cells, CD4+ T cells, CD8+ T cells, DCs, γδ T cells, macrophages/monocytes, neutrophils, NK cells, and NKT cells using canonical marker level measurements (Fig. 5E, Figure S6D). Consistent with our scRNA-seq analysis results, T cells (CD4+ and CD8+ T) accounted for a large proportion of CD45+ TILs, and no significant difference was observed in the percentages of B cells, CD4 T cells, CD8 T cells, DCs, NK cells, macrophages and neutrophils (Fig. 5F). However, the levels of a proliferation marker (Ki-67) and inflammatory cytokine markers, including IFN-γ, IL-17 A, and TNF-α were significantly higher in CD4+ T cells but not CD8+ T cells in OCCC than in HGSC (Fig. 5G, Figure S6E).
Our study is the first to explore the cellular and molecular differences in primary tumor tissues between OCCC and HGSC using single-cell RNA sequencing and GEO gene expression signatures. We partially explained the chemoresistance of OCCC cells via a cell cycle analysis. We also identified a possible pathogenic epithelial subcluster with overexpressed LEFTY1 in OCCC samples. New biomarkers distinguishing OCCC samples from HGSC samples were identified, among which AOC1 and S100A2 were the first to be reported. Additionally, we showed that metabolic pathways were activated in OCCC cells, which indicates that they might be promising candidates for new therapeutic strategies. In addition, we characterized and compared the immune cellular profiles of tumor-infiltrating lymphocytes between OCCC and HGSC.
In the past two decades, RNA sequencing has become a ubiquitous tool used in molecular biology, and RNA sequencing results have had significant implications on the choices of specific molecular biomarkers and clinical approaches [12, 13]. Several studies that have applied scRNA-seq data to HGSC studies have suggested possible antitumor targets [14, 15], components of the tumor environment [16, 17], and chemotherapy resistance mechanisms . However, the results of scRNA-seq analysis with OCCC samples have not been reported to date.
Cell cycle phase and SCENIC analyses in this study partially explained the chemoresistance mechanism underlying OCCC oncogenesis. Recent studies discovered that cisplatin caused cell death mainly by preventing RNA transcription  and induced cells to arrest in the G2 phase instead of the S phase . The lower proportion of G2/M phase tumor cells that we found in OCCC samples might contribute to its higher platinum resistance than that of HGSC.
In this study, the dividing epithelial cell subset (Cluster 5) was more likely to develop into the LEFTY1+ subset (Cluster 4) in the OCCC group than in the HGSC group. The LEFTY1+ subset probably plays an important role in the pathogenesis of OCCC. LEFTY1 is a novel member of the TGF-beta superfamily and modulates the epithelial-mesenchymal transition and cancer stem cell properties in the context of OCCC [21, 22]. Other cancer-related pathways enriched in Cluster 4 cells, such as “Hedgehog signaling”, “Wnt-beta catenin signaling” and “Notch signaling”, had been reported in some studies on HGSC [23,24,25], but they had been rarely studied in the OCCC context. Further studies are needed to investigate whether inhibiting these signaling pathways may be an effective clinical therapy for OCCC.
In addition, fewer intensive CNVs were observed in OCCC cells than in HGSC cells. In a pancancer analysis, high CNV was associated with significantly worsened overall survival . Late-stage tumors carry distinct features ; for example, CNV burden accumulates to a higher level in the deeper lesions of squamous esophageal carcinoma . The higher CNV expression in HGSC cells may be a result of the late stages of HGSC patients compared with OCCC patients from which the samples were taken.
New diagnostic biomarkers were identified via the combined analysis of single-cell RNA sequencing and GEO gene expression data. Neither AOC1 nor S100A2 had been reported to distinguish OCCC samples from HGSC samples. AOC1, a copper-containing amine oxidase, plays different roles in different carcinomas. In gastric cancer and hepatocellular carcinoma, AOC1 functions as an oncogene, either by activating the AKT signaling pathway and EMT  or by regulating the IL-6/JAK/STAT3 pathway . However, AOC1 promotes ferroptosis and inhibits prostate cancer progression . In our study, AOC1 was a promising feature marker of OCCC cells. Additionally, platinum-sensitive patients were more likely to overexpress AOC1 than platinum-resistant patients, with a p value close to being statistically significant (p = 0.055). S100A2, a member of the largest subfamily of calcium-binding EF-hand type proteins, has been found to be related to poor prognosis in patients with ovarian cancer . Its expression is significantly higher in the advanced stage than in the early stage of an EOC . In our study, the expression of S100A2 was detected more frequently in HGSC than in OCCC samples. S100A2 may be used to differentiate OCCC from HGSC samples.
Consistent with published studies, LEFTY1 and GPX3 were also identified in this study as potential OCCC molecular biomarkers. A previous study applied shotgun proteomics and identified LEFTY1 as a specific OCCC molecular marker, showing that this protein exerted an antitumor effect by reducing cell proliferation and promoting cisplatin-induced apoptosis . Higher expression of GPX3 was discovered and found to be related to a better prognosis and chemosensitivity in patients with OCCC in our findings, which was in line with a previous study . Additionally, we found a lower proportion of primary OCCC tumor tissues of OCCC expressing WFDC2 and CRABP2 compared than was found in HGSC tissues. WAP four-disulfide core domain 2 (WFDC2/HE4) is an EOC clinical biomarker [35,36,37,38], and cellular retinoic acid-binding protein 2 (CRABP2) is upregulated in ovarian cancer and contributes to tumor growth and tumor cell migration and invasion [39,40,41].
We found that inhibiting glucose metabolism by metformin or 2-deoxy-D-glucose effectively decreased the viability of OCCC and HGSC cell lines. Metformin is a well-characterized medication used for type 2 diabetes, and increasing evidence has highlighted its potential use as an affordable, well-tolerated, and effective agent anticancer and antiaging therapy [42,43,44]. Metformin has also been demonstrated to be a novel treatment option for ovarian cancer by either preventing tumor growth or increasing chemotherapy sensitivity [45,46,47,48,49,50]. A recent Phase II clinical trial evaluating metformin therapy suggested a better-than-expected overall survival for patients with ovarian cancer . Our work first suggests an antitumor effect of metformin on OCCC patients because it suppresses tumor growth and promotes tumor cell apoptosis. Therefore, metformin can also be explored as a promising candidate for OCCC clinical treatment.
OCCC is often associated with endometriosis and is highly linked to high oxidative stress. Quenching reactive oxygen species (ROS) is expected to improve patient outcomes . Glutathione is the most abundant cellular antioxidant that can neutralize ROS. Ferroptosis is an iron-dependent cell death associated with drug resistance in many cancers [53,54,55]. In endometriosis-related ovarian cancer, such as OCCC, repeated bleeding caused by endometriosis contributes to the accumulation of iron, which stimulates the Fenton reaction, produces ROS, and further induces ferroptosis . Elevated intracellular iron levels are closely associated with ovarian cancer, and ferroptosis-inhibiting genes have been found to be related to ovarian cancer progression and chemoresistance [57,58,59,60,61]. Since our data suggest that inhibiting glutathione metabolism, which might lead to an increased ROS level, or inhibiting ROS-induced ferroptosis can greatly attenuated chemosensitivity in OCCC, it will be interesting to investigate whether and how ROS participate in OCCC chemoresistance.
Immune cells are important components in the tumor environment and exert an impact on treatment efficiency and prognosis [62,63,64,65]. Different immune cells participate in immune responses to tumors through different pathways [66, 67]. The scRNA-seq data combined with mass cytometry data in our study showed that T cells (CD4+ and CD8+ T) constitute a large proportion of CD45+ TILs in OCCC samples. The proportions of immune cell subsets did not differ. However, the results were limited by a low cell count and sequencing depth. Further study of the immune environment of OCCC and HGSC is needed.
There are, however, a number of restrictions on this study. To validate the results, we only performed TMA and in vitro cell line experiments. Further functional and mechanistic experiments are needed, in particular, murine in vivo or organoid models, to confirm our findings. Additionally, as was mentioned in the CNV section, selection bias might have occurred because to differences in the stages of the OCCC and HGSC patients.
In summary, our work provides holistic insights into the cellular and molecular differences between OCCC and HGSC samples. Specifically, we identified characteristic biomarkers, specific cell types, and metabolically related pathways involved in the tumorigenesis and drug resistance of OCCC in patients. Our study serves as a valuable resource for use in discovering new targets that improve the poor prognosis of OCCC patients.
Materials and methods
All patients underwent surgical resection at Peking Union Medical College Hospital. None of the patients had received therapy before surgery. Mixed carcinoma was ruled out in selected participants. In addition, patients with double primary carcinomas, namely, with both ovarian and uterine tumors, were excluded. The diagnosis was confirmed by 3 experienced pathologists, who evaluated tumors according to the 2020 World Health Organization classification , and staging was performed using The NCCN Guidelines for Ovarian Cancer (Version 3.2022) . Clinical data, such as residual tumor after surgery, were gathered using the criteria of previous investigations [70, 71]. Therapeutic information included chemotherapy and and other therapies were also documented . Chemoresistance was outlined as progression, persistent disease while receiving maintenance therapy, or full remission followed by relapse within six months of finishing platinum-based chemotherapy [65, 69].
Single-cell suspension preparation
Single-cell preparation for scRNA-seq and mass cytometry validation was performed according to a previous study . In brief, tumors were cut into approximately 1 mm3 pieces and digested in Hank’s balanced salt solution (HBSS) containing 1.5 mg/ml collagenase IV (Sigma‒Aldrich), 1 mg/ml hyaluronidase (Sigma‒Aldrich) and 500 µg/ml DNase I (GoldBio) while rotating at 200 rpm and 37 °C for 30 min. Cell suspensions were subsequently passed through a 70-µm cell strainer (BD, Biosciences), followed by centrifugation at 400 × g for 10 min. Red blood cells were lysed with RBC lysis buffer (Miltenyi Biotec, 130-094-183) for 2 min on ice, and dead cells were removed with a dead cell removal kit (Miltenyi Biotec, 130-090-101) following the manufacturer’s instructions.
Library construction and single-cell sequencing
Library construction was conducted using a Chromium Next GEM Single Cell 3’ GEM Library & Gel Bead Kit v3.1 (10×Genomics, USA), and sequencing was performed with an Illumina NovaSeq 6000 (Novogene, Beijing, China).
RNA-seq data processing
Raw sequencing data generated with the 10×Genomics platform were processed following the standard Chromium Cell Ranger pipeline (version 4.0.0) against the GRCh38 human reference genome. We filtered data and created Individual Seurat objects converted from single-cell counts from all samples using the Seurat analysis package in R (v4.1.0) according to a procedure described in a previous study . We then merged all samples and performed a canonical correlation analysis (CCA) to reduce dimensionality and remove the batch effect. We then performed standard cell clustering using the ScaleData function. We used the RunPCA function to calculate the principal component analysis (PCA) dimensions and the FindNeighbors and FindCluster functions for unsupervised clustering of the data. Finally, we used the RunUMAP function for cell visualization.
Cell subclustering analysis
Epithelial cells, immune cells (T cells, NK cells, myeloid cells, plasma B cells, B cells, and plasmacytoid dendritic cells (pDCs)), and stromal cells were extracted on the basis of integrated data for further subclustering. A similar process of clustering was performed for each major cell type to obtain cell subtypes. Doublet clusters were removed following the following criteria:  epithelial cells based on CD3 expression (calculated as the mean expression of CD3D, CD3E, and CD3G) > 0.1;  stromal cells based on CD3D, C1QA, EPCAM, and CD14 expression.
Differentially expressed genes and pathway enrichment
To identify differentially expressed genes (DEGs), we used the FindMarkers function with two-sided unpaired Wilcoxon tests, and p values were adjusted following the Benjamin & Hochberg protocol. DEGs were filtered with a minimum log2(fold change) of 0 and a maximum adjusted p value of 0.05. GO and KEGG pathway enrichment analyses were performed based on the DEGs using the R package ClusterProfiler (v4.0.5) .
A trajectory analysis of the epithelial cells was performed using Seurat and Monocle 3 R package software (v1.0.0) designed by Cao et al. . In brief, we used the GetAssayData function to fetch the raw expression matrix, created the CellDataSet, and used the preprocess_cds function to normalize and preprocess data. The cluster_cells and learn_graph functions were used for the trajectory inference analysis, and the results were visualized using the UMAP function.
CNV and SCENIC analysis
To identify malignant cells, we used InferCNV software (v1.2.0) to estimate the number of CNVs in each region . The CNV was calculated based on the expression level for each cell with a cutoff of 0.1. T cells were used as the reference. We used the R package pySCENIC (v0.10.0) to analyze the enrichment of transcriptome factors in each epithelial subtype . The activity of each regulon was evaluated using AUCell.
Total RNA was extracted using TRIzol reagent (Thermo Fisher Scientific) and reverse transcribed into cDNA using PrimeScript™ RT Master Mix (TAKARA, Japan). Real-time PCR was performed using Fast SYBR Green Master Mix (Thermo Fisher Scientific) and a QuantStudio™ 7 Flex Real-Time PCR System (Thermo Fisher Scientific). Gene expression was normalized to GAPDH. Tables S3 and S4 provide detailed information on the primers used.
Cells were lysed with RIPA lysis buffer (AOqing Biotechnology, Beijing, China) containing ProtLytic Protease Inhibitor Cocktail (New Cell & Molecular Biotech), and the protein concentration was determined with a BCA assay kit (Pierce Biotechnology, USA). The cell lysate was fractionated in NuPAGE™ 4 to 12% gels (Thermo Fisher Scientific) and transferred onto polyvinylidene difluoride (PVDF) membranes (Millipore, USA). The membranes were blocked and incubated with primary antibodies, including antibodies against LEFTY1, AOC1, GPX3, CRABP2, WFDC2, and S100A2, and then, the commensurate secondary antibody. The bands were visualized with enhanced chemiluminescence following the manufacturer’s instructions (Peirce).
Gene set enrichment (GSE) analysis
We obtained profiling of gene expressions in primary ovarian cancer specimens (OCCC, n = 11, HGSC, n = 8) in the GSE189553 dataset from the GEO database (http://www.ncbi.nlm.nih.gov/geo/). DEGs between OCCC specimens and HGSC specimens were identified using the limma package (v3.50.3).
Estimating the abundance of immune cell types via Immune cell abundance identifier (ImmunCellAI)
To estimate the abundance of 24 immune cell types, log2-transformed expression data from RNA-Seq results were compiled and uploaded to the web server for ImmunCellAI (http://bioinfo.life.hust.edu.cn/web/ImmuCellAI/) analysis as previously described .
Immunohistochemical staining and interpretation
The IHC validation cohort consisted of patients diagnosed with OCCC and enrolled via consecutive sampling. The patients underwent surgical resection between January 2019 and May 2022 at Peking Union Medical College Hospital with enough archived tissue for immunohistochemical testing. Finally, 128 patients with OCCC and 81 patients with HGSC were included. The clinical information was listed in Table S5. Two tissue microarrays were constructed with formalin-fixed, paraffin-embedded tissues. IHC staining was performed with DAKO Autostainer Link 48 and measured according to the standard as reported in our previous studies [78, 79]. The primary antibodies were listed in Table S6.
The location of expressed protein and the tissue control in immunohistochemical staining were also listed in Table S7 [80,81,82,83,84,85]. Composite scores (range from 0 to 12) were calculated based on the intensity and percentage of stained cells [81, 86]. Cases with composite scores of 4 points or more were termed high expression, while the others were labeled low expression.
The TOV21G cell line was purchased from ATCC (Manassas, VA, USA) and cultured in MCDB105 medium (Zhong Qiao Xin Zhou Biotechnology, Shanghai, China). The OVCAR3 cell line was purchased from the National Collection of Authenticated Cell Cultures (Shanghai, China) and cultured in RPMI-1640 medium (Zhong Qiao Xin Zhou Biotechnology, Shanghai, China) supplemented with 10% FBS (Gibco, USA) and 1% penicillin/streptomycin (Invitrogen, USA). Both cell lines were maintained in an incubator at 37 °C with 5% CO2.
Cells were seeded at a density of 5000 cells per well in 96-well plates (Corning Life Sciences, USA), grown overnight, and exposed to different concentrations of cisplatin (HY-17,394, MedChemExpress) and treated with different inhibitors, including 2-deoxy-D-glucose (HY-13,966, MedChemExpress), metformin (PHR1084, Merck), CB-839 (HY-12,248, MedChemExpress), liproxstatin-1 (HY-12,726, MedChemExpress) and ferrostatin-1 (HY-100,579, MedChemExpress), for 24 h, 48 h, 72 h, followed by 4 h of incubation with 10 µl of CCK8 reagent (Dojindo, Shanghai, China). The absorbance at 450 nm was measured using a microplate reader (Bio–Rad, USA).
Apoptosis rate determined by flow cytometry
Cells were seeded at a density of 5 × 105 cells per well in a 6-well plate, grown overnight, and treated with cisplatin combined with different inhibitors as described above. The cell apoptosis rate was assessed with a FITC Annexin V Apoptosis Detection Kit (BD, Biosciences) following the manufacturer’s instructions. The data were obtained with an Attune NxT flow cytometer (Invitrogen).
Cells were seeded at a density of 2 × 105 cells per well in a 12-well plate, grown overnight, and treated with cisplatin combined with different inhibitors as described above. Cell proliferation was determined using a Click-iT® EdU Imaging Kit (C10337, Invitrogen) according to the manufacturer’s protocol. The nucleus was stained with 1 × Hoechst 33,342 solution (5 µg/mL). Images of five randomly selected areas were taken with a fluorescence microscope (Zeiss, Germany).
GraphPad Prism software Version 8.0 (GraphPad, San Diego, CA) and SPSS Version 24.0 (SPSS Inc., Chicago, IL, USA) were used for statistical analyses and graphic presentations. The data are shown as the mean ± standard deviation (SD). Differences among multiple groups were compared by one-way ANOVA or Kruskal‒Wallis test with Bonferroni post hoc test based on a normal distribution and homogeneity of variance. Pearson’s chi-squared test and Fisher’s exact test were performed to evaluate binary variables, such as immunohistochemical results or application of chemotherapy. Univariate and multivariable analyses based on progression-free survival (PFS) were performed using Cox regression. The proportional hazards assumption was evaluated by analyzing the statistical significance of interactions between exposure and the follow-up time. A P value < 0.05 was considered statistically significant.
The accession number for the raw sequencing data has been deposited in the Genome Sequence Archive for humans (GSA-Human) under accession number HRA002536. To comply with the “Guidance of the Ministry of Science and Technology (MOST) for the Review and Approval of Human Genetic Resources”, we are required to deposit the genomic data of Chinese people under controlled access at the GSA in Beijing Institute of Genomics Data Center. To access the raw data under accession number HRA002536, please submit requests to the GSA-Human online page for this study (https://ngdc.cncb.ac.cn).
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
Torre LA, Trabert B, DeSantis CE, Miller KD, Samimi G, Runowicz CD, et al. Ovarian cancer statistics, 2018. CA Cancer J Clin. 2018;68(4):284–96.
Zhu C, Zhu J, Qian L, Liu H, Shen Z, Wu D, et al. Clinical characteristics and prognosis of ovarian clear cell carcinoma: a 10-year retrospective study. BMC Cancer. 2021;21(1):322.
Chan JK, Teoh D, Hu JM, Shin JY, Osann K, Kapp DS. Do clear cell ovarian carcinomas have poorer prognosis compared to other epithelial cell types? A study of 1411 clear cell ovarian cancers. Gynecol Oncol. 2008;109(3):370–6.
Kao YC, Lin MC, Lin WC, Jeng YM, Mao TL. Utility of hepatocyte nuclear factor-1β as a diagnostic marker in ovarian carcinomas with clear cells. Histopathology. 2012;61(5):760–8.
Tsuchiya A, Sakamoto M, Yasuda J, Chuma M, Ohta T, Ohki M, et al. Expression profiling in ovarian clear cell carcinoma: identification of hepatocyte nuclear factor-1 beta as a molecular marker and a possible molecular target for therapy of ovarian clear cell carcinoma. Am J Pathol. 2003;163(6):2503–12.
Oda K, Ikeda Y, Kashiyama T, Miyasaka A, Inaba K, Fukuda T, et al. Characterization of TP53 and PI3K signaling pathways as molecular targets in gynecologic malignancies. J Obstet Gynaecol Res. 2016;42(7):757–62.
Wu N, Zhang X, Wang Z, Zhang X, Fang C, Li H et al. Progesterone prevents HGSOC by promoting precancerous cell pyroptosis via inducing fibroblast paracrine. iScience (6107). 2023;26(4):106523.
Fraungruber P, Kaltofen T, Heublein S, Kuhn C, Mayr D, Burges A, et al. G protein-coupled estrogen receptor correlates with Dkk2 expression and has prognostic impact in Ovarian Cancer Patients. Front Endocrinol (Lausanne). 2021;12:564002.
Mungenast F, Thalhammer T. Estrogen biosynthesis and action in ovarian cancer. Front Endocrinol (Lausanne). 2014;5:192.
Itamochi H, Kigawa J, Terakawa N. Mechanisms of chemoresistance and poor prognosis in ovarian clear cell carcinoma. Cancer Sci. 2008;99(4):653–8.
Stark R, Grzelak M, Hadfield J. RNA sequencing: the teenage years. Nat Rev Genet. 2019;20(11):631–56.
Zhang J, Späth SS, Marjani SL, Zhang W, Pan X. Characterization of cancer genomic heterogeneity by next-generation sequencing advances precision medicine in cancer treatment. Precis Clin Med. 2018;1(1):29–48.
Izar B, Tirosh I, Stover EH, Wakiro I, Cuoco MS, Alter I, et al. A single-cell landscape of high-grade serous ovarian cancer. Nat Med. 2020;26(8):1271–9.
Wan C, Keany MP, Dong H, Al-Alem LF, Pandya UM, Lazo S, et al. Enhanced efficacy of simultaneous PD-1 and PD-L1 Immune Checkpoint Blockade in High-Grade Serous Ovarian Cancer. Cancer Res. 2021;81(1):158–73.
Olalekan S, Xie B, Back R, Eckart H, Basu A. Characterizing the tumor microenvironment of metastatic ovarian cancer by single-cell transcriptomics. Cell Rep. 2021;35(8):109165.
Stur E, Corvigno S, Xu M, Chen K, Tan Y, Lee S, et al. Spatially resolved transcriptomics of high-grade serous ovarian carcinoma. iScience (6107). 2022;25(3):103923.
Zhang K, Erkan EP, Jamalzadeh S, Dai J, Andersson N, Kaipio K, et al. Longitudinal single-cell RNA-seq analysis reveals stress-promoted chemoresistance in metastatic ovarian cancer. Sci Adv. 2022;8(8):eabm1831.
Forgie BN, Prakash R, Telleria CM. Revisiting the Anti-Cancer toxicity of clinically approved platinating derivatives. Int J Mol Sci. 2022;23(23).
Sorenson CM, Eastman A. Mechanism of cis-diamminedichloroplatinum(II)-induced cytotoxicity: role of G2 arrest and DNA double-strand breaks. Cancer Res. 1988;48(16):4484–8.
Akiya M, Yamazaki M, Matsumoto T, Kawashima Y, Oguri Y, Kajita S, et al. Identification of LEFTY as a molecular marker for ovarian clear cell carcinoma. Oncotarget. 2017;8(38):63646–64.
Matsumoto T, Yokoi A, Hashimura M, Oguri Y, Akiya M, Saegusa M. TGF-β-mediated LEFTY/Akt/GSK-3β/Snail axis modulates epithelial-mesenchymal transition and cancer stem cell properties in ovarian clear cell carcinomas. Mol Carcinog. 2018;57(8):957–67.
van der Ploeg P, Uittenboogaard A, Bosch SL, van Diest PJ, Wesseling-Rozendaal YJW, van de Stolpe A, et al. Signal transduction pathway activity in high-grade serous carcinoma, its precursors and fallopian tube epithelium. Gynecol Oncol. 2022;165(1):114–20.
Chehover M, Reich R, Davidson B. Expression of wnt pathway molecules is associated with disease outcome in metastatic high-grade serous carcinoma. Virchows Archiv: an international journal of pathology. 2020;477(2):249–58.
Huang CC, Cheng SH, Wu CH, Li WY, Wang JS, Kung ML, et al. Delta-like 1 homologue promotes tumorigenesis and epithelial-mesenchymal transition of ovarian high-grade serous carcinoma through activation of notch signaling. Oncogene. 2019;38(17):3201–15.
Tao Z, Wang S, Wu C, Wu T, Zhao X, Ning W et al. The repertoire of copy number alteration signatures in human cancer. Brief Bioinform. 2023;24(2).
Martínez-Jiménez F, Movasati A, Brunner SR, Nguyen L, Priestley P, Cuppen E, et al. Pan-cancer whole-genome comparison of primary and metastatic solid tumours. Nat (69504). 2023;618(7964):333–41.
Liu Z, Su R, Ahsan A, Liu C, Liao X, Tian D et al. Esophageal squamous Cancer from 4NQO-Induced mice Model: CNV alterations. Int J Mol Sci. 2022;23(22).
Xu F, Xu Y, Xiong JH, Zhang JH, Wu J, Luo J, et al. AOC1 contributes to Tumor Progression by promoting the AKT and EMT pathways in gastric Cancer. Cancer Manag Res. 2020;12:1789–98.
Ding Q, Lin D, Zhou Y, Li F, Lai J, Duan J, et al. Downregulation of amine oxidase copper containing 1 inhibits tumor progression by suppressing IL-6/JAK/STAT3 pathway activation in hepatocellular carcinoma. Oncol Lett. 2021;22(6):857.
Ding Y, Feng Y, Huang Z, Zhang Y, Li X, Liu R, et al. SOX15 transcriptionally increases the function of AOC1 to modulate ferroptosis and progression in prostate cancer. Cell Death Dis. 2022;13(8):673.
Xu HY, Song HM, Zhou Q. Comprehensive analysis of the expression and prognosis for S100 in human ovarian cancer: a STROBE study. Med (Baltim). 2020;99(47):e22777.
Lou W, Ding B, Zhong G, Du C, Fan W, Fu P. Dysregulation of pseudogene/lncRNA-hsa-mir-363-3p-SPOCK2 pathway fuels stage progression of ovarian cancer. Aging. 2019;11(23):11416–39.
Lee HJ, Do JH, Bae S, Yang S, Zhang X, Lee A, et al. Immunohistochemical evidence for the over-expression of glutathione peroxidase 3 in clear cell type ovarian adenocarcinoma. Med Oncol. 2011;28(Suppl 1):522–7.
Huang J, Chen J, Huang Q. Diagnostic value of HE4 in ovarian cancer: a meta-analysis. Eur J Obstet Gynecol Reprod Biol. 2018;231:35–42.
James NE, Gura M, Woodman M, Freiman RN, Ribeiro JR. A bioinformatic analysis of WFDC2 (HE4) expression in high grade serous ovarian cancer reveals tumor-specific changes in metabolic and extracellular matrix gene expression. Med Oncol. 2022;39(5):71.
Scaletta G, Plotti F, Luvero D, Capriglione S, Montera R, Miranda A, et al. The role of novel biomarker HE4 in the diagnosis, prognosis and follow-up of ovarian cancer: a systematic review. Expert Rev Anticancer Ther. 2017;17(9):827–39.
Dochez V, Caillon H, Vaucel E, Dimet J, Winer N, Ducarme G. Biomarkers and algorithms for diagnosis of ovarian cancer: CA125, HE4, RMI and ROMA, a review. J Ovarian Res. 2019;12(1):28.
Li N, Lin G, Zhang Y, Zhang Q, Zhang H. Exosome-related protein CRABP2 is upregulated in ovarian carcinoma and enhances cell proliferation. Discov Oncol. 2022;13(1):33.
Zhang C, Pan R, Ma S, Xu S, Wang B. Dezocine inhibits cell proliferation, migration, and invasion by targeting CRABP2 in ovarian cancer. Open Med (Wars). 2022;17(1):2052–61.
Xie T, Tan M, Gao Y, Yang H. CRABP2 accelerates epithelial mesenchymal transition in serous ovarian cancer cells by promoting TRIM16 methylation via upregulating EZH2 expression. Environ Toxicol. 2022;37(8):1957–67.
Kamarudin MNA, Sarker MMR, Zhou JR, Parhar I. Metformin in colorectal cancer: molecular mechanism, preclinical and clinical aspects. J Exp Clin Cancer Res. 2019;38(1):491.
Wheaton WW, Weinberg SE, Hamanaka RB, Soberanes S, Sullivan LB, Anso E, et al. Metformin inhibits mitochondrial complex I of cancer cells to reduce tumorigenesis. Elife. 2014;3:e02242.
Podhorecka M, Ibanez B, Dmoszyńska A. Metformin - its potential anti-cancer and anti-aging effects. Postepy Hig Med Dosw (Online). 2017;71(0):170–5.
Lengyel E, Litchfield LM, Mitra AK, Nieman KM, Mukherjee A, Zhang Y, et al. Metformin inhibits ovarian cancer growth and increases sensitivity to paclitaxel in mouse models. Am J Obstet Gynecol. 2015;212(4):479. .e1-.e10.
Dos Santos Guimarães I, Ladislau-Magescky T, Tessarollo NG, Dos Santos DZ, Gimba ERP, Sternberg C, et al. Chemosensitizing effects of metformin on cisplatin- and paclitaxel-resistant ovarian cancer cell lines. Pharmacol Rep. 2018;70(3):409–17.
Shank JJ, Yang K, Ghannam J, Cabrera L, Johnston CJ, Reynolds RK, et al. Metformin targets ovarian cancer stem cells in vitro and in vivo. Gynecol Oncol. 2012;127(2):390–7.
Yang C, Zhao N, Li D, Zou G, Chen Y. Metformin improves the sensitivity of ovarian cancer cells to chemotherapeutic agents. Oncol Lett. 2019;18(3):2404–11.
Li L, Wang L, Li J, Fan Z, Yang L, Zhang Z, et al. Metformin-Induced reduction of CD39 and CD73 blocks myeloid-derived suppressor cell activity in patients with ovarian Cancer. Cancer Res. 2018;78(7):1779–91.
Singh SK, Apata T, Singh S, McFadden M, Singh R. Clinical implication of Metformin in Relation to Diabetes Mellitus and Ovarian Cancer. Biomedicines. 2021;9(8).
Brown JR, Chan DK, Shank JJ, Griffith KA, Fan H, Szulawski R et al. Phase II clinical trial of metformin as a cancer stem cell-targeting agent in ovarian cancer. JCI Insight. 2020;5(11).
Amano T, Murakami A, Murakami T, Chano T. Antioxidants and therapeutic targets in Ovarian Clear Cell Carcinoma. Antioxid (Basel). 2021;10(2).
Zhang C, Liu X, Jin S, Chen Y, Guo R. Ferroptosis in cancer therapy: a novel approach to reversing drug resistance. Mol Cancer. 2022;21(1):47.
Chaudhary N, Choudhary BS, Shah SG, Khapare N, Dwivedi N, Gaikwad A, et al. Lipocalin 2 expression promotes tumor progression and therapy resistance by inhibiting ferroptosis in colorectal cancer. Int J Cancer. 2021;149(7):1495–511.
Yang WH, Huang Z, Wu J, Ding CC, Murphy SK, Chi JT. A TAZ-ANGPTL4-NOX2 Axis regulates ferroptotic cell death and chemoresistance in epithelial ovarian Cancer. Mol Cancer Res. 2020;18(1):79–90.
Yamaguchi K, Mandai M, Toyokuni S, Hamanishi J, Higuchi T, Takakura K, et al. Contents of endometriotic cysts, especially the high concentration of free iron, are a possible cause of carcinogenesis in the cysts through the iron-induced persistent oxidative stress. Clin cancer research: official J Am Association Cancer Res. 2008;14(1):32–40.
Zhang C, Liu N. Ferroptosis, necroptosis, and pyroptosis in the occurrence and development of ovarian cancer. Front Immunol. 2022;13:920059.
You Y, Fan Q, Huang J, Wu Y, Lin H, Zhang Q. Ferroptosis-related gene signature promotes ovarian Cancer by influencing Immune Infiltration and Invasion. J Oncol. 2021;2021:9915312.
Ma LL, Liang L, Zhou D, Wang SW. Tumor suppressor mir-424-5p abrogates ferroptosis in ovarian cancer through targeting ACSL4. NEOPLASMA (3409). 2021;68(1):165–73.
Hong T, Lei G, Chen X, Li H, Zhang X, Wu N, et al. PARP inhibition promotes ferroptosis via repressing SLC7A11 and synergizes with ferroptosis inducers in BRCA-proficient ovarian cancer. Redox Biol. 2021;42:101928.
Wang Y, Zhao G, Condello S, Huang H, Cardenas H, Tanner EJ, et al. Frizzled-7 identifies platinum-tolerant ovarian Cancer cells susceptible to ferroptosis. Cancer Res. 2021;81(2):384–99.
Sue AQR, Patel PG, Shakfa N, Nyi MN, Afriyie-Asante A, Kang EY, et al. Prognostic significance of T cells, PD-L1 immune checkpoint and tumour associated macrophages in clear cell carcinoma of the ovary. Gynecol Oncol. 2021;162(2):421–30.
Ye S, Li Q, Wu Y, Jiang W, Zhou S, Zhou X, et al. Integrative genomic and transcriptomic analysis reveals immune subtypes and prognostic markers in ovarian clear cell carcinoma. Br J Cancer. 2022;126(8):1215–23.
Matsuura H, Miyamoto M, Hada T, Ishibashi H, Iwahashi H, Kakimoto S et al. The worsening impact of programmed cell death ligand 1 in ovarian clear cell carcinomas. Arch Gynecol Obstet. 2022.
Zhu J, Wen H, Bi R, Wu Y, Wu X. Prognostic value of programmed death-ligand 1 (PD-L1) expression in ovarian clear cell carcinoma. J gynecologic Oncol. 2017;28(6):e77.
Farhood B, Najafi M, Mortezaee K. CD8(+) cytotoxic T lymphocytes in cancer immunotherapy: a review. J Cell Physiol. 2019;234(6):8509–21.
Guillerey C. NK cells in the Tumor Microenvironment. Adv Exp Med Biol. 2020;1273:69–90.
Höhn AK, Brambs CE, Hiller GGR, May D, Schmoeckel E, Horn L-C. 2020 WHO classification of female genital tumors. Geburtshilfe Frauenheilkd. 2021;81(10):1145–53.
Armstrong DK, Alvarez RD, Backes FJ, Bakkum-Gamez JN, Barroilhet L, Behbakht K, et al. NCCN Guidelines® insights: ovarian Cancer, Version 3.2022. J Natl Compr Canc Netw. 2022;20(9):972–80.
Hu D, Ma D, Zhang ZJ, Zhang Y, Huang K, Li X. Prognosis comparison between small cell carcinoma of ovary and high-grade serous ovarian cancer: a retrospective observational cohort study. Front Endocrinol (Lausanne). 2023;14:1103429.
Bekes I, Löb S, Holzheu I, Janni W, Baumann L, Wöckel A, et al. Nectin-2 in ovarian cancer: how is it expressed and what might be its functional role? Cancer Sci. 2019;110(6):1872–82.
Wu X, Liu Y, Jin S, Wang M, Jiao Y, Yang B, et al. Single-cell sequencing of immune cells from anticitrullinated peptide antibody positive and negative rheumatoid arthritis. Nat Commun. 2021;12(1):4977.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7.
Qiu X, Hill A, Packer J, Lin D, Ma YA, Trapnell C. Single-cell mRNA quantification and differential analysis with Census. Nat Methods. 2017;14(3):309–15.
Patel AP, Tirosh I, Trombetta JJ, Shalek AK, Gillespie SM, Wakimoto H, et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science. 2014;344(6190):1396–401.
Aibar S, González-Blas CB, Moerman T, Huynh-Thu VA, Imrichova H, Hulselmans G, et al. SCENIC: single-cell regulatory network inference and clustering. Nat Methods. 2017;14(11):1083–6.
Miao YR, Zhang Q, Lei Q, Luo M, Xie GY, Wang H, et al. ImmuCellAI: a Unique Method for Comprehensive T-Cell subsets abundance prediction and its application in Cancer Immunotherapy. Adv Sci (Weinh). 2020;7(7):1902880.
Gao Y, Wang P, Lu J, Pan B, Guo D, Zhang Z, et al. Diagnostic significance of parafibromin expression in parathyroid carcinoma. Hum Pathol. 2022;127:28–38.
Guo D, Zhao X, Wang A, Xie Q, Xu X, Sun J. PD-L1 expression and association with malignant behavior in pheochromocytomas/paragangliomas. Hum Pathol. 2019;86:155–62.
Owen RP, White MJ, Severson DT, Braden B, Bailey A, Goldin R, et al. Single cell RNA-seq reveals profound transcriptional similarity between Barrett’s oesophagus and oesophageal submucosal glands. Nat Commun. 2018;9(1):4261.
Liu F, Ou W, Tang W, Huang Z, Zhu Z, Ding W, et al. Increased AOC1 expression promotes Cancer Progression in Colorectal Cancer. Front Oncol. 2021;11:657210.
Yang ZL, Yang L, Zou Q, Yuan Y, Li J, Liang L, et al. Positive ALDH1A3 and negative GPX3 expressions are biomarkers for poor prognosis of gallbladder cancer. Dis Markers. 2013;35(3):163–72.
Li X, Wang C, Wang S, Hu Y, Jin S, Liu O, et al. YWHAE as an HE4 interacting protein can influence the malignant behaviour of ovarian cancer by regulating the PI3K/AKT and MAPK pathways. Cancer Cell Int. 2021;21(1):302.
Liu RZ, Li S, Garcia E, Glubrecht DD, Poon HY, Easaw JC, et al. Association between cytoplasmic CRABP2, altered retinoic acid signaling, and poor prognosis in glioblastoma. Glia. 2016;64(6):963–76.
Kumar M, Srivastava G, Kaur J, Assi J, Alyass A, Leong I, et al. Prognostic significance of cytoplasmic S100A2 overexpression in oral cancer patients. J translational Med. 2015;13:8.
He Y, Wang Y, Li P, Zhu S, Wang J, Zhang S. Identification of GPX3 epigenetically silenced by CpG methylation in human esophageal squamous cell carcinoma. Dig Dis Sci. 2011;56(3):681–8.
This study is supported by the National High Level Hospital Clinical Research Funding 2023-PUMCH-F-004 (D.G.& J.S.), 2022-PUMCH-B-062 (J.S. ), 2022-PUMCH-D-002 (D.G. ), the CAMS Innovation Fund for Medical Sciences (CIFMS) 2021-I2M-1-053 (D.G.& Z.-Y.L.) and the National Natural Science Foundation of China (No. 8197544), Qingyun Plan (PB-QY012108) (X.-Y.W.).
Ethics approval and consent to participate
This study was approved by the Institutional Review Board of Peking Union Medical College Hospital (JS-3553). Informed consent was obtained from all patients.
Consent for publication
Not applicable for that section.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
About this article
Cite this article
Guo, D., Zhang, S., Gao, Y. et al. Exploring the cellular and molecular differences between ovarian clear cell carcinoma and high-grade serous carcinoma using single-cell RNA sequencing and GEO gene expression signatures. Cell Biosci 13, 139 (2023). https://doi.org/10.1186/s13578-023-01087-3
- Single-cell RNA-seq
- Bulk RNA-seq