Article Text
Abstract
Background Ovarian clear cell carcinoma (OCCC) is a rare and chemo-resistant subtype of ovarian cancer. While immunotherapy has demonstrated effectiveness in some OCCC cases, the mechanisms for heterogeneous immunoreactivity and potential combinatory strategies remain unclear.
Methods Tumor samples from 13 patients with OCCC underwent single-cell mRNA-seq and TCR-seq to generate 1 40 683 cells transcriptome, while additionally 31 formalin-fixed paraffin-embedded samples were used for immunohistochemistry. Spatial transcriptomics of two OCCC samples and bulk RNA-seq of 58 patients were incorporated for spatial and interpatient level explorations. Serum tumor markers and radiologic images of three patients with OCCC who received combinatory VEGF and PD-1 inhibition were retrospectively analyzed.
Results OCCC exhibited a dynamic immune architecture shaped by genetic and therapeutic pressure. ARID1A mutation linked to baseline immune activation, correlated with an enrichment of neoantigen-reactive CXCL13+ CTLA4+ CD8+ T cells (p<0.001) and enhanced FASLG–FAS interactions. Recurrent OCCC was fibrotic, angiogenic, and immunosuppressive, exhibiting metabolic reprogramming towards activated activity in fatty acid metabolism. High CD36 (log-rank p=0.012, HR: 4.515) and CD47 expression (log-rank p=0.037, HR: 3.246) indicated worse progression-free survival. Treatment with bevacizumab increased intratumoral T cell infiltration and activated T cell interferon-γ signaling. Retrospective analysis of clinical cases revealed that combination therapy with anti-VEGF (vascular endothelial growth factor) and anti-PD-1 agents exerted clinical benefits in patients with OCCC with persistent, recurrent, and metastatic disease.
Conclusions ARID1A mutation correlated with OCCC baseline immune activation. Stromal reconstruction and tumor metabolic reprogramming functioned as key processes of OCCC dynamic progression. VEGF inhibition remodeled OCCC stroma, restored T cell function and potentiated immunotherapy. CD36 and CD47 might be potential therapeutic targets for recurrent OCCC.
- immunotherapy
- ovarian cancer
- tumor microenvironment - TME
- escape/evasion
- combination therapy
Data availability statement
Data are available upon reasonable request. The data generated in this study are available within the article and supplementary data files. Other data or materials are available from the corresponding author on reasonable request.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
WHAT IS ALREADY KNOWN ON THIS TOPIC
Advanced-stage ovarian clear cell carcinoma (OCCC) is chemo-insensitive and associated with unfavorable prognosis. Immunotherapy, especially when combined with other agents, emerges as a novel therapeutic approach for OCCC.
WHAT THIS STUDY ADDS
ARID1A mutation in OCCC linked to baseline immune-activation with enrichment of CXCL13+ CTLA4+ tumor-reactive T cells. OCCC dynamically manipulated the host environment through tumor-stroma interaction. Upregulation of the tumor metabolic molecules CD36 and CD47, as well as the stromal marker POSTN, were associated with unfavorable clinical outcomes. VEGF inhibition deconstructed tumor stroma and facilitated T cell activation.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
These findings uncovered a dynamic and multidimensional immune landscape in OCCC, shedding light on the impact of tumor-stroma interplay on recurrence and therapeutic resistance. VEGF inhibition combined with anti-PD1 might be a promising therapeutic strategy for OCCC.
Introduction
Ovarian cancer (OC) is a highly aggressive gynecologic malignancy, leads to approximately 200 000 women yearly deaths worldwide.1 Ovarian clear cell carcinoma (OCCC) is relatively rare, accounting for around 10% of OC, which is often diagnosed at an earlier FIGO (International Federation of Gynecology and Obstetrics) stage.2 However, of all pathologic subtypes, advanced-stage OCCC correlates with the worst prognosis,3–6 with a 3 year overall survival (OS) rate of 53.3% for stage III patients and 29.6% for stage IV patients. Moreover, recurrent OCCCs are platinum-resistant with a less than 10% objective response rate (ORR), underscoring the urgent need for novel therapeutic approaches. Unlike high-grade serous cancers (HGSCs) that are initiated by TP53 inactivation and deficiencies in homologous recombination (HR), OCCCs are driven by somatic alterations in PI3K/AKT/mTOR pathway (PIK3CA for ~50%) and SWI/SNF chromatin remodeling complex (ARID1A for ~50%, ARID1B for 6%–18%, SMARCA4 for 5%–18%).7 Somatic alterations shape the intratumoral heterogeneity, fuel tumor evolution under extrinsic stresses such as metastatic cloning and chemotherapy, and effected on the tumor microenvironment. Single-cell technology has offered an exceptional opportunity to understand the tumor plasticity and tissue architecture at a higher resolution and comprehensive view,8–10 paves the way for a deeper understanding of OCCC heterogeneity and mechanisms of therapeutic resistance.
Recent years, immunotherapy emerges as a promising therapeutic avenue for refractory OC, OCCC especially,11 with several case studies reporting a subset of cases that exhibit strong and sustained responses to immune checkpoint inhibitors (ICIs).12–14 Epithelial OC showed restrained response to single-agent PD-1/PD-L1 inhibitors, with the ORR primarily ranged around 10%, including PD-1 inhibitors pembrolizumab (8%),15 nivolumab (7.6%16–15%),17 and PD-L1 inhibitors avelumab (3.7%),18 BMS-936559 (6.0%),19 and atezolizumab (22.2%).20 Although OC might be ‘immune-cold’ and insensitive to ICIs according to above clinical trials, encouragingly, subgroup analysis implied the potential extraordinary sensitivity of OCCC to PD-1/PD-L1 inhibitors. In phase II KEYNOTE-100 study, patients with OCCC showed a trend toward greater benefits from pembrolizumab (n=19, ORR=15.8%, 95% CI: 3.4% to 39.6%).15 In phase III study of platinum-resistant OC, patients with OCCC experienced longer OS with nivolumab compared with chemotherapy (gemcitabine or pegylated liposomal doxorubicin) (HR: 0.78; 95% CI: 0.46 to 1.32).16 Moreover, a transcriptional clear cell signature was significantly higher in complete response (CR) patients with nivolumab treatment.21 Although larger-scaled clinical trials are still needed to verify above findings, great importance has been attached to clear cell pathology in the ongoing clinical studies on OC immunotherapy.
However, insufficient molecular insights into OCCC immune heterogeneity impeded the development for biomarker discovery and combinatory immunotherapy. While studies in the overall OC population have suggested a series of ICI predictors, including tumor mutational burden,22 PD-L1 expression, CD8+ T cells infiltration,15 18 transcriptional signatures of proliferation,21 glycolysis and hypoxia,22 and myeloid cell infiltration,22 little is known about the specific indicators for OCCC immunotherapy. As prevalent genomic events of OCCC, ARID1A mutations and deficient mismatch repair have been previously associated with DNA mutability and downstream immune activation.7 11Accordantly, a case study reported that patients with OCCC harboring ARID1A mutations might be particularly sensitive to immunotherapy.14 Interestingly, established evidence by immunohistochemistry (IHC) and bulk RNA sequencing have demonstrated the negative clinical correlations of tumor-infiltrating CD8+ T cells, CD4+ T cells, and CD68+ macrophages,23 24 highlighting the significant immune dysregulation within OCCC. However, little is known about the tissue-specific immunosuppression mechanisms of OCCC, as well as the dynamical processes of how OCCC respond to immunotherapy. It is essential to further explore both the tumor–host interactions specific to OCCC and the detailed functional status of intratumoral immune cells.
Recent years, the application scenarios of single-cell sequencing have been expanded from live tissues to wider-ranged samples including formalin fixed paraffin-embedded (FFPE) blocks.25 scFFPE-seq provided an opportunity to retrospectively collected clinical samples, thus facilitated investigations of rare diseases, as well as combinatory analysis of corresponding follow-up data. With the advancement in materials availability in single-cell technology, we designed this study with various source of OCCC samples, exploiting both fresh tumors and FFPE samples, to fill in the gap in translational research on OCCC immunobiology.
In the present study, we established a panoramic immune atlas of OCCC at single-cell, spatial and interpatient dimensions. Analyses of treatment-naïve tumors revealed the impact of ARID1A mutation on baseline immune activity. Longitudinal investigations in recurrent tumors demonstrated tumor-stroma interplay as acquired mechanism of immune unresponsiveness. Finally, VEGF inhibition showed the efficacy to reduce stromal density, empower host immunity and might serve as combinatory strategy to facilitate efficient immunotherapy.
Methods
Human specimen collection
All datasets used in this study are listed in online supplemental table 1. Three in-house datasets involving human subjects were enrolled in this study (dataset 1–3). The study procedures were designed according to the Declaration of Helsinki. Human specimen collection was approved by the ethic committee of Fudan University Shanghai Cancer Center (FUSCC) with the approval number No. 050432-4-2108*. Fresh tumors were collected from five patients with OCCC (dataset 1) during surgery in FUSCC from 2021 to 2023, immersed in DMEM and delivered to the laboratory on ice immediately for tissue dissociation. From 2014 to 2023, tumor samples from 39 patients (dataset 2 and 3) who were pathologically diagnosed with OCCC were collected and preserved as FFPE samples at 4℃. Those patients were followed-up for at least 1 year after surgery, and progression-free survival (PFS) was set as the endpoint of study. Tumor samples from those patients were retrospectively collected and constructed into tumor microarrays for IHC. All patients included in this study approved the utilization of corresponding data and signed written informed consent.
Supplemental material
Isolation of cancer-associated fibroblasts (CAFs)
Primary CAFs were isolated from ascites of patients with OC. Briefly, fresh ascites samples were collected from drainage bags, transported on ice and passed through a 70 µm cell strainer to remove large impurities. Subsequently, ascites was centrifuged at 800 rcf for 10 min, and resuspended with DMEM. The red blood cells were removed by Red Blood Cell Lysis Solution (Solarbio), and remaining cells were cultured with Dulbecco's modified eagle medium (DMEM) (without fetal bovine serum). Fibroblasts were separated and purified by differential adhesion method. One hour after cell seeding, the supernatant was collected into a new dish, and adherent cells were gently washed with PBS for three times. Those adherent cells were mainly fibroblasts, cultured with DMEM (10% FBS) at 37℃ and were passaged every 4–6 days according to cell density.
CAF co-culture system
Every 3 days, the culturing supernatant of fibroblasts (passage 3–5) was collected and centrifuged at 800 rcf for 3 min to remove cell pellet, adding to the tumor culture dish (supernatant:fresh medium=1:1). Tumor cell lines and cell proliferation assay were described in online supplemental data files.
Immunohistochemistry
The slides were first baked at 65℃ for 2 hours. Then, slides were immersed into three washes of dimethylbenzene, and gradient concentration ethanol (SANGON, 100%, 95%, 90%, and 85%) every 10 min in sequence for deparaffinage. Briefly, slides were immersed into antigen retrieval solution (chosen according to the manufactory’s instruction) and heated with a pressure cooker for 3 min after boiling. Thereafter, endogenous peroxidase was blocked by 0.3% H2O2 for 10 min at room temperature. Incubation of goat serum for 30 min at 37℃ was used for non-specific antigen blocking. Subsequently, slides underwent primary antibody incubation (primary antibodies: anti-rabbit human CD36 antibody, 1:400, proteintech, cat no. 18 836-1-AP; anti-mouse human POSTN antibody, 1:4000, proteintech, cat no. 66 491-1-Ig; human CD47 IHC kit, cat no. RMA-0842) overnight at 4℃ and secondary antibody (proteintech, cat no: RGAR011, RGAM001) incubation at 37℃ for 45 min. Coloration of slides by diaminobenzidine (DAB, Solarbio) lasted for 1 min and was terminated by dilated water. For counterstaining, the colored tissue sections were placed into hematoxylin for 2 min. Finally, slides were dehydrated with gradient concentration ethanol (85%, 95%, 100%) every 2 min and underwent xylene transparent for 1–2 min before mounting. Evaluation of IHC was performed by two independent oncologists who were blinded to clinical information. For evaluation of CD36 and CD47, slides with >50% strong staining (brown) tumor area were defined as ‘high expression’. For POSTN, slides with >20% strong staining (brown) stromal area were defined as ‘high expression’.
Whole-exon sequencing and identification of ARID1A mutations
Genomic DNA from tumor samples was extracted and sequenced by the Illumina NovaSeq platform. GATK4 was used to identify single nucleotide polymorphisms (SNPs) and insertions/deletions (INDELs) within the exon regions and the exon–intron boundary regions. Mutations were annotated with ANNOVAR. Two rounds of quality control were conducted to ensure the validity of mutation recognition. First, SNPs below the following thresholds were filtered out: QD <2.0, MQ <40.0, FS >60.0, SOR >3.0, MQRankSum <−12.5, and ReadPosRankSum <−8.0. Similarly, the INDEL filtering criteria were set at QD <2.0, FS >200.0, SOR >10.0, MQRankSum <−12.5, and ReadPosRankSum <−8.0. Subsequently, mutations with a variant allele frequency (VAF) <5% or mutant allele number <5 or allele sum <20 were filtered out. ARID1A mutations were identified according to OncoKB26 (https://www.oncokb.org/). Patients with ARID1A deletions and truncating mutations (loss-of function, and likely pathogenic alterations as suggested as OncoKB) were regarded as ARID1A mutant group, other patients who were not detected with those ARID1A alterations were classified as ARID1A wildtype group (online supplemental table 3).
Single-cell RNA sequencing and scFFPE-seq
Sample preparation, library construction and sequencing were performed according to 10× genomics pipeline. Fresh tumor samples were collected in the DMEM, transported on ice to the laboratory immediately during surgery after sampling. Tumors were minced into pieces as small as possible, and digested with a Tumor Dissociation Kit (Cat. 130-096-730; Miltenyi Biotec) for 30 minutes on a rotator at 37℃. Subsequently, the dissociated cells were passed through a 70 µm cell-strainer (Corning) until single cell suspensions were obtained. The red blood cells were removed by Red Blood Cell Lysis Solution (Solarbio). Single-cell RNA-seq and TCR-seq libraries were prepared using the Chromium Next GEM Single Cell 5′ Kit v2 and Chromium Single Cell Human TCR Amplification Kit from 10× Genomics. FFPE blocks were conducted as 50 µm FFPE curls, and were dissociated with the Miltenyi Biotech FFPE Tissue Dissociation Kit following the manufacturer’s instructions and underwent Chromium Fixed RNA Profiling (scFFPE-seq).25 Above libraries were sent for sequencing based on the Illumina NovaSeq 6000 platform. Single-cell quality control is described in online supplemental data files.
Identification of metaprograms (MPs) by non-negative matrix factorization (NMF) method
MPs of epithelial cells and fibroblasts were identified by the NMF method using GeneNMF package. The merged Seurat object was first split by patient, and consistent NMF MPs were computed across multiple samples with the function ‘multiNMF’. NMF MPs were annotated according to program genes and corresponding ‘hallmark’ and ‘GO-BP’ pathway activities.
Functional and pathway analysis
Functional pathway activity was analyzed with integrated approaches including AUCell, UCell, and GSVA. Transcriptional factors (TFs) activity of each cell subclusters were predicted by DecoupleR.27 Copy number alterations of epithelial cells were identified by inferCNV of the Trinity CTAT Project (https://github.com/broadinstitute/inferCNV) with T cells and macrophages as the normal controls. Metabolism28 framework was used to evaluate metabolic activity of cell clusters.
Trajectory analysis
Briefly, monocle2,29 CytoTRACE,30 and slingshot31 were used comprehensively to predict the build an evolutionary relationship and trajectory among functional cell subsets. CytoTRACE was used to evaluate cell stemness, the subcluster with the highest CytoTRACE score were presumed at the starting point of developmental route. Monocle2 and slingshot were used to infer the evolutionary trajectory in single-cell-level and subcluster-level, respectively.
Intercellular ligand–receptor interaction analysis
CellChat32 and NicheNet33 were used to speculate intercellular interaction within the OCCC tumor microenvironment. Overall interaction number and strength between cell types, as well as focused pathway activity were inferred by CellChat. Nichenet set cell type of interest as the signal receiver and tested how other cells in the tumor microenvironment impact its phenotype. Therefore, we conducted a receiver-focused ligand–receptor analysis with Nichenet to examine how TIME affected effector cells, T cells particularly.
Analysis of spatial transcriptome data
Spatial transcriptomics data were processed by Seurat v4 pipeline. SCT transform was performed and spatial-resolved gene expression was visualized by the function ‘SpatialFeaturePlot’. Spots were clustered with the resolution 0.5, and classified into ‘Epithelial-rich’, ‘Stromal-rich’, ‘Endothelial-rich’ and ‘Immune-rich’ (if presented) regions according to marker genes expression. Cell type signature scores were added to the meta data by UCell method. CellChat V.2 was used for spatial cell–cell communication analysis.
Analysis of bulk RNA sequencing cohorts
Bulk RNA-seq data of OC containing OCCC histologies were deconvoluted by cell type signature scores with ssGSEA methods. Ecotyper34 was used for immune ecotype prediction.
Statistical analysis
All graphics and statistical analysis were conducted using Prism V.8.3.0, IBM SPSS V.20 and R V.4.3.2. Detailed packages and computing frameworks are listed in online supplemental table 2. For data fit the normal distribution and homogeneity of variance, parametric tests were priorly chosen. Data that not pass the normality test and homogeneity test of variances were log-transformed or analyzed by non-parametric test as appropriate. All statistical analyses were performed two-sided. P<0.05 was regarded as statistically significant: *p<0.05, **p<0.01, ***p<0.001.
Data availability
All datasets used in this study are listed in online supplemental table 1. The following datasets were obtained from GEO (http://www.ncbi.nlm.nih.gov/geo/): GSE18488035 (dataset 5); GSE22433536 (dataset 4); GSE22687037 (dataset 6). Survival analysis on pan-cancer immunotherapy cohort was performed online using KM-plotter38 39 (https://kmplot.com/analysis/) (dataset 7). The data generated in this study are available within the article and online supplemental data files. Other data or code are available from the lead contact on reasonable request.
Results
Prospective and retrospective collection of OCCC samples for single-cell profiling
OCCC is relatively rare and poorly differentiated, which poses challenges for fresh sample collection. To provide a comprehensive immune atlas of OCCC, we integrated data from two in-house sequencing cohorts, one prospective and the other retrospective, for single-cell profiling (figure 1A, online supplemental table 1). Detailed clinical information of these patients is shown in online supplemental table 3. The prospective cohort (dataset 1) comprised three recurrent (RC) and two newly diagnosed (ND) fresh tumors obtained during first surgery and confirmed pathologically postsurgery. These samples underwent 5' single-cell RNA-seq and scTCR-seq to generate a high-quality cell type matrix.
The retrospective cohort (dataset 2) comprised eight FIGO stage III OCCC formalin-fixed paraffin-embedded (FFPE) samples, including five ND patients with more than 3 year PFS and three RC patients, for scFFPE-seq analysis. Among the five ND patients, two were identified with somatic ARID1A frameshift mutations, while the other two had wildtype ARID1A. All three RC cases underwent adjuvant chemotherapy after the initial surgery; two of them received standard TC (Taxol plus Carboplatin) chemotherapy, another one patient received a TC plus VEGFi (bevacizumab) regimen.
For IHC validation, tumor samples from 31 patients with OCCC who underwent surgery at FUSCC between 2014 and 2023 were collected (dataset 3). In addition, one spatial transcriptomics cohort consisting of two patients with OCCC (dataset 4), one external scRNA-seq cohorts of seven HGSC and five normal ovaries (dataset 5), and two bulk RNA-seq datasets (dataset 6 and 7) were obtained for spatial and interpatient level investigations.
Recurrent OCCC is fibrotic, angiogenic, and immunosuppressive
To construct a high-quality single-cell atlas of OCCC, we first analyzed fresh tumor samples (dataset 1). Following quality control, the remaining 49 228 cells were harmonized to remove batch effects and visualized using the Uniform Manifold Approximation and Projection (UMAP) method (figure 1B). The OCCC ecosystem comprised epithelial, immune, and stromal compartments, which intricately interacted. The cells were categorized into 12 major cell types (online supplemental S1A, B, online supplemental table 4). Additionally, a subset of cells overexpressing proliferation and cell cycle markers were classified as MKI67+ proliferating cells. Compared with newly diagnosed tumors, recurrent OCCC exhibited a fibrotic and angiogenic TIME. All non-immune stroma lineages including fibroblasts, endothelial cells, and pericytes were significantly enriched in recurrent cases (figure 1C). Contrastingly, immune cells were organized heterogeneously, with canonical dendritic cells (cDCs), classic monocytes, M1-polarized macrophages (online supplemental figure S2A, B), and IgM+ plasma cells (online supplemental figure S2C, D) more abundant in newly diagnosed tumors, while M2-polarized macrophages, naive B cells and IgG+ plasma cells accumulated in recurrent tumors.
The major cell types, including epithelial cells, stromal cells, myeloid cells, and T cells, were identifiable in the 91 455 FFPE-based single-cell mixture (dataset 2, figure 1D). To elucidate robust cell–cell interactions within the OCCC microenvironment, we performed CellChat analyses on both fresh tumors (online supplemental S1C) and FFPE samples (online supplemental figure S1D). Consistent findings emerged from both datasets, revealing increased outgoing secreting signaling from fibroblasts, endothelial cells, and macrophages in recurrent tumors, involving networks of periostin, complement, and galectin.
Somatic ARID1A mutation correlates with baseline immune activation
Consistent with existing data,40 41 somatic ARID1A mutations were found to be prevalent in our internal OCCC dataset (5/8 in all sequenced patients, online supplemental table 3). When stratifying treatment-naive tumors according to ARID1A status (dataset 2), increased T cells and macrophages but decreased fibroblasts were found in ARID1A mutant tumors. OCCC with somatic ARID1A mutation shaped an immune-enriched and less fibrotic TIME (figure 1E). For broader validation, Ecotyper34 was employed to deconvolute an OC bulk RNA-seq dataset (dataset 6, GSE226870). A higher proportion of ARID1A mutant tumors were classified as carcinoma ecotypes (CE) 8 and CE9, which correlated with longer survival. CE9 was also identified as a biomarker for immunotherapy, prompting further investigation into the immune landscape of ARID1A mutant OCCC.
ARID1A mutation also impacted on the cell–cell communication profile (online supplemental figure S1E). Most interestingly, the FASLG–FAS interaction between T cells and malignant cells uniquely detected in ARID1A mutant tumors. Surprisingly, the recipient of this signaling was identified as MKI67+ proliferative cells rather than the entire epithelial population, suggesting reduced proliferation as a potential mechanism to evade immune killing. Additionally, T cells in ARID1A mutant tumors exhibited activated colony stimulating factor (CSF) signaling that involved in myeloid cells recruitment, offering an explanation for the enrichment of macrophages and macrophage-derived signals.
OCCC epithelial cells are characterized by hypoxia and hypoxic metabolism
To characterize the features and heterogeneity of OCCC epithelial cells, we first focused on cells identified by fresh tumor scRNA-seq (dataset 1). Epithelial cells were observed with vibrant gene copy number variations (CNV) (online supplemental figure S3A), suggesting a malignant status. Consistent with existing experimental studies, transcription factor (TF) analysis revealed HNF1B 40 as a key governor of OCCC cells, along with other TFs including SMARCA4, SMARCA1, YAP1 and AR (online supplemental table 5). To profile the subtype-specific features, pathway activity of epithelial cells from different histologies, including OCCC, HGSC, and normal ovary, was compared using the non-parametric method AUCell (figure 2A). Notably, OCCC cells exhibited heightened activity in hypoxia and hypoxic metabolism pathways including glycolysis, adipogenesis, and fatty acid (FA) metabolism.
To dissect the key biological processes of OCCC, NMF was employed to extract transcriptomic (MPs (figure 2B). Among the 10 identified MPs, only MP2 was associated with cell cycling, while MP5, MP6, and MP9 were linked to epithelial-to-mesenchymal transition (EMT). Additionally, MP1, MP3, MP8, and MP9 indicated cellular hypoxia, and MP3, MP4, and MP9 were correlated with immune activation, again underscoring hypoxia as a central player in OCCC development.
When inquiring how genomic events effect on transcriptional programs, in scFFPE-seq dataset, ARID1A mutant tumors showed activated ROS, hypoxia, EMT, cholesterol homeostasis, and adipocyte development pathways (figure 2C). By simulating the evolutionary trajectory, ARID1A mutant tumor were positioned on a branch of epithelial cells, demonstrating sustained elevated expression of ISG15 and MXRA8 during evolution.
Recurrent tumor cells demonstrate signature of immune evasion and metabolic switch towards FA metabolism
To profile the heterogeneity of epithelial population, five functional subsets were identified from scRNA-seq dataset (figure 2D, online supplemental figure S3B). Cluster 1 (EC1) and cluster 2 (EC2), respectively, represented the majority of ND and RC tumor cells (online supplemental figure S3C). Evolutionary analysis positioned EC2 as a branch of EC1 (figure 2E), with a metabolic transition from glycolysis to active FA metabolism along the trajectory. Additionally, EC2 exhibited altered gene expression profiles, including decreased HLA-A expression and NFKB activity, increased SNAI2 activity, and loss of E-cadherin expression. By deconvoluting a bulk RNA-seq dataset (online supplemental table 6), we validated that EC1 correlated with immune activation and CD8+ T/NK cells infiltration while EC2 was associated with EMT and FA metabolism, at the interpatient level (online supplemental figure S3D, E).
After overlapping trajectory genes identified in fresh sample and FFPE samples, conservative recurrence-related processes were identified, which included FA metabolism and adipocyte development (online supplemental figure S4A, B). Specifically, two FA metabolism-related molecules, CD36 and CD47, were significantly upregulated in EC2 (figure 2F). To investigate the clinical significance of CD36 and CD47, we performed IHC staining on tumor samples from 31 patients with OCCC (dataset 3, figure 2G, online supplemental figure S4C). Interestingly, all patients with prior neoadjuvant chemotherapy showed upregulation of CD36 and CD47 (Pearson χ2 test p=0.063 and 0.019, respectively). Both CD36 (log-rank p=0.012, HR: 4.515, 95% CI: 1.249 to 16.32) and CD47 abundance (log-rank p=0.037, HR: 3.246, 95% CI: 1.015 to 10.38) correlated worse PFS, indicating their potential linkages with platinum resistance (figure 2H). However, only a borderline statistical significance (Pearson χ2 p=0.06) was observed regarding the correlation between CD36 expression and platinum-free interval (PFI), while no significant correlation was observed between the CD47 expression level and PFI group (online supplemental table 7), calling on larger-scaled cohort analysis or experimental studies to investigate the roles of FA metabolism pathway in OC platinum resistance.
GeneMANIA analysis showed that SIRPA and THBS1 were key microenvironmental regulators for CD36 and CD47 (online supplemental figure S4D), with myeloid cells and non-immune stromal cells being the major sources of SIRPA and THBS1, respectively (online supplemental figure S4E). Spatial analysis validated the colocalization of those ligand–receptor pairs in OCCC tumor (dataset 4), with CD47-THBS1 being the most prevalent pair, suggestive of the impact of non-immune stromal cells on tumor phenotype.
T cell phenotypes and neoantigen-reactive subsets in OCCC
Intratumoral T cells function as the primary anticancer effectors and are key determinants of immunotherapy sensitivity.42 43 Our analysis towards the overall OCCC ecosystem indicated comparative CD8+ T cells abundance between recurrent tumors and newly diagnosed ones (figure 1C), however, their intercellular communication patterns was altered (online supplemental figure S1C). This phenomenon prompted us to further dissect the functional continuum of CD8+ T cells. Therefore, we first classified T cells into CD4+, CD8+ and NK cells (figure 3A) and recognized six heterogenous subsets of CD8+ T cells within OCCC (figure 3B, online supplemental figure S5A).
Notably, compared with those in normal ovaries, CD8+ T cells in OCCC were more exhausted, displayed elevated expression of immune checkpoint markers such as PDCD1, LAG3, CD39, and CTLA444, alongside reduced expression of effector molecules including IFNG, GZMB, PRF1, CD107A, and KLRG1 (figure 3C). Two hyperexpanded CD8+ T cell subclusters, CXCL13+ CTLA4+ cluster 1 (enriched in ND), and GZMH+ CCL5+ cluster 2 (enriched in RC), were potentially neoantigen-responsive based on TCR repertoire analysis (figure 3D). Compared with other subsets, CXCL13+ CTLA4+ CD8+ T subset was terminally differentiated (online supplemental figure S5B), and showed a hyper-activated phenotype with mixed activation and exhaustion markers (F3E, (online supplemental figure S5C, D). Remarkably, T cells in ARID1A-mutated OCCCs exhibited higher CXCL13+ CTLA4+ CD8+ T phenotype scores (figure 3F). The presence of CXCL13+ CTLA4+ CD8+ T cells was associated with improved OS in patients receiving immunotherapy (log-rank p<0.001, HR=0.51), suggesting a potentially reversible immune dysregulation state. Ligand–receptor activity analysis uncovered fibroblasts-derived CXCL12 as a key ligand affecting CD8+ T cells (figure 3G, online supplemental figure S5E), which was previously linked to immune exclusion45–47 and tumor metastasis,48 again highlighting stromal cells as key player in mediating OCCC immunosuppression.
As immunoregulatory subsets, CD4+ T cells in OCCC also displayed dysfunctional phenotype, with naive and metabolic subsets accumulated in recurrent tumors (figure 3H–J). Five distinct subclusters were identified within the CD4+ T cell population (figure 3H and online supplemental figure S6A–D). Infiltrating CD4+ T cells in OCCC displayed reduced production of Th1 cytokines (IFNG, TNF, and IL2), yet elevated expression of regulatory markers (FOXP3 and IL2RA), dysfunctional markers (CTLA4, TIGIT, and PDCD1), and memory-related genes (CCR7, CD69, and CD45RA) (figure 3I). Compared with ND tumors, the diversity of CD4+ T cells in recurrent OCCC was found to be constrained. Two subclusters, namely the naive-like C1 (IL7R+ KLF2+) and hypermetabolic C4 (KLRB1+ LGALS3+), constituted the majority of the CD4+ T cell population in recurrent tumors.
Tumor-promoting myofibroblasts are accumulated in recurrent OCCC
Non-immune tumor stroma plays multifaced roles in tumor progression. The current study has indicated the accumulation and enhanced signaling of fibroblasts, pericytes, and endothelial cells during recurrence. We sought to identify the chief culprit throughout the heterogeneous tumor stroma.
Consistent in both fresh tumors and FFPE samples, myofibroblasts (mCAFs) were accumulated in recurrent tumors (figure 4A, B, online supplemental figure S7A–C). A unique subset of chondrocyte-like COLA2+ fibroblasts was found only in scFFPE-seq (online supplemental figure S7D). MPs of fibroblasts were identified using the NMF method from integrated scRNA-seq and scFFPE-seq data (figure 4C). MP3, linked to myogenesis and angiogenesis, was increased in recurrent OCCC and associated with mCAFs (figure 4D). Within fibroblast subsets, mCAF_C1 with the most active myogenesis process (online supplemental figure S7E) was located at the origin of the developmental trajectory (figure 4A), presenting as a stem-like recurrence-related subset. This subset was controlled by TGF-β signaling-related TFs including MYOCD, SRF, and SMAD1 (online supplemental figure S7F) and showed activated lipid metabolism (online supplemental figure S7G). Interestingly, mCAFs, especially the POSTN+ subset, were the main source of THBS1 (figure 4D), and therefore might be key regulators of tumor phenotype.
We conducted IHC staining of POSTN in OCCC tissues (figure 4E). Denser POSTN expression was observed in recurrent OCCC compared with newly diagnosed tumors (online supplemental figure S7H). High POSTN expression alone showed a borderline significant correlation with unfavorable PFS (log-rank p=0.056). However, when tumor CD36/CD47 levels were included in the analysis, patients with tumors that were double positive for CD36 or CD47 and stroma POSTN exhibited the significantly worst PFS (log-rank p=0.015) (figure 4F). Additionally, the addition of culture supernatant from cancer-associated fibroblasts isolated from surgical samples (online supplemental figure S7I) to an OCCC cell line culture system promoted tumor cell growth (figure 4G) and induced CD36 expression (online supplemental figure S7J).
Endothelial cells and pericytes form a vascular niche and correlate with unfavorable immune ecotype
Pericytes and endothelial cells are main components of tumor vascular structure. To visualize the spatial colocation of those cells in OCCC, UCell method was used to annotate each spot in spatial transcriptomic data (figure 5A and online supplemental figure S8A). The ‘Endothelial-rich’ regions showed typically codistribution of endothelial cells and pericytes, representing the vascular structure within OCCC tumors, consistent with histologic findings. Naive CD4+ T cells were primarily located around the vessels rather than in the tumor core, which partially explains their lack of expression of immune activation or regulatory molecules. The tumor stroma trapped a group of immunosuppressive cell subsets, including M2-polarized macrophages and metabolic CD4+ T cells. In contrast, tumor cores showed minimal immune infiltration. Additionally, the chemo-sensitive OCCC formed an immune-reactive interface characterized by an enrichment of DCs, M1-polarized macrophages, monocytes, apCAFs, and cluster 1 epithelial cells (online supplemental figure S8A).
Interpatient level analysis revealed the associations between cell subsets and ecosystem construction. Intriguingly, patients with OCCC with infiltrated TIME were classified into two distinct ecotypes, EC1 and EC9, which showed reverse clinical associations (online supplemental figure S8B). The most significant difference between the two subsets was the abundance of stromal cells. Among all cell subsets, B lineage was the only one that predicted a favorable carcinoma ecotype (online supplemental figure S8C, top; AUC: 0.709, p=0.011). On the contrary, pericytes (AUC: 0.101, p<0.001), mCAFs (AUC: 0.201, p<0.001), and endothelial cells (AUC: 0.283, p=0.008) were negative predictors of a good CE type.
VEGF inhibition remodels the tumor stroma and reinvigorates tumor-infiltrating T cells
Given the association with the immunosuppressive ecotype of the vascular structure, we investigated whether VEGF inhibition could improve the TIME. Comparing patients who received merely TC (Taxol plus Carboplatin) chemotherapy with those who received combinatory TC plus bevacizumab, we found that the application of bevacizumab increased the proportion of intratumoral T cells while reducing the abundance of fibroblasts (figure 5B). T cells in VEGFi-treated tumors displayed higher scores for CD8 cytotoxicity, CD8 IFN response, CD4 activation effector function, CD4 IFN response, but lower levels of CD4 naïve score,39 indicating an immune activated status (figure 5C). Cell–cell communication patterns were also altered in VEGFi-treated OCCC (figure 5D). Fibroblast-derived POSTN and SEMA3 signaling were reduced in VEGFi-treated OCCC, while CCL and galectin networks were upregulated. The patterns of VEGF signaling differed between tumors treated with TC and TC plus VEGFi. Upregulated VEGF signaling from myeloid cells and epithelial cells in VEGFi-treated OCCC partially accounted for the recurrence after adjuvant therapy.
VEGFi plus anti-PD-1 exerts clinical benefit in refractory OCCC
To evaluate the efficacy of VEGFi plus anti-PD1 in OCCC in a clinical and real-world context, we conducted a retrospective analysis of clinical cases. Clinical benefits from VEGFi plus PD1 inhibition were observed in persistent (figure 6A), recurrent (figure 6B), and metastatic (figure 6C) patients with OCCC, either combined with adjuvant chemotherapy or as maintenance treatment.
Case 1 involved a patient with ARID1A mutant OCCC with persistent disease after 2 cycles of TC adjuvant chemotherapy. The patient’s serum CA125 levels returned to the normal range after adding VEGFi plus anti-PD1 to the chemotherapy regimen (from 154 U/mL to 11.5 U/mL), and achieved radiologic partial remission (figure 6A). Case 2 was a patient with ARID1A wildtype OCCC who experienced platinum-resistant recurrence. The patient’s tumor markers remained within normal limits until 9 months after the first surgery (figure 6B). PET-CT indicated new lesions of retroperitoneal lymph node metastasis. While TC plus VEGFi treatment failed to control the elevated serum CA199, the application of GP (gemcitabine+cisplatin) + VEGFi + anti-PD1 combinatory treatment resulted in partial remission. Case 3 involved a patient with ARID1A wildtype stage IV OCCC with liver metastasis. The patient’s serum CA125 levels decreased from 227 U/mL to 27.1 U/mL after two cycles of DO (docetaxel+oxaliplatin) chemotherapy plus VEGFi and anti-PD-1 (figure 6C). While well-designed prospective clinical trials are still needed to accurately assess the rate of clinical benefit of this combinatory therapy in OCCC, the cases discussed above provide preliminary evidence of the efficacy of VEGFi plus anti-PD-1 in refractory OCCC, irrespective of ARID1A status and disease course (figure 6D).
Discussion
Advanced-stage OCCCs present a significant challenge due to their high malignancy and resistance to chemotherapy. To facilitate development of OCCC immunotherapy, we collected OCCC samples from all available sources, including fresh tumors and FFPE samples, and used single-cell technologies to establish a multidimensional and dynamic immune atlas of OCCC. OCCC exhibits distinct genomic features, with approximately 50% of cases harboring somatic ARID1A mutations.49 Consistent with previous studies in other malignancies,50 ARID1A mutation in OCCC linked to baseline immune activation and pre-existing tumor-reactive T cells.
Growing clinical trials are carrying out to systemically assess the efficacy of ICIs in OCCC, especially in combination with front-line pharmaceuticals and other immunomodulatory agents. For combinatory PD-1/CTLA4 blockade, nivolumab plus ipilimumab achieved an approximately threefold ORR (31.4%, 16/51) than nivolumab monotherapy, with clear cell histology more likely (logistic regression OR: 5.205) to respond when compared with other histologic subtypes.51 Combination therapy with novel immunomodulatory agents including IDO1 inhibitor epacadostat (ORR: 21.4%)52 are also evaluated. Interestingly, PD-1/PD-L1 blockade plus antiangiogenic agents exerted noticeable clinical effects. Pembrolizumab in combination with bevacizumab achieved 28.9% ORR53 in recurrent OC, and 47.5% when combined with oral metronomic cyclophosphamide.54 In a phase I study of the anti-PD-1/VEGFA bispecific antibody AK112, two out of three patients with platinum-resistant OCCC achieved a partial response (ORR: 66.7%).55 This ORR exceeded the average level in both patients with OC (5/19, ORR: 26.3%) and the entire solid tumor cohort (12/47, ORR: 25.5%). Emerging clinical trials are designed to test various ICI-based therapeutic combinations in recurrent and refractory OCCC. Those combinations mainly focused on the combination of PD-1/PD-L1 inhibitors and antiangiogenic drugs including Anlotinib (NCT05600998), Sintilimab (NCT04735861), Lenvatinib (NCT05296512), and AK112 (NCT06560112). To expand the therapeutic arsenal of OCCC, efforts are also made in investigating ICIs in combination with novel inhibitors, including BRCA1/2 inhibitor Pamiparib (NCT05044871), TIGIT inhibitor Etigilimab (NCT05715216), PP2A inhibitor LB-100 (NCT06065462), and HER2 ADC drug Disitamab vedotin (RC48) (NCT06540729). Particularly, one study brings AK104 (PD-1/CTLA-4 bispecific) plus chemotherapy into the context of neoadjuvant treatment (NCT06542549). Another study examining the combination of ATR inhibitor AZD6738, Olaparib and durvalumab takes ARID1A mutation into study design (NCT04065269). In the following decade, results from these trials will provide clinical-level evidences to facilitate individualized ICI-based therapies for patients with OCCC with diverse molecular events.
In this context, a major translational significance gained by our study is to biologically explain why patients with OCCC might benefit from anti-PD-1 plus VEGFi. During recurrence, OCCC tumor cells orchestrated an immunosuppressive microenvironment through dynamic reconstruction of the tumor stroma. Spatial analysis showed that the colocalization of endothelial cells and RGS5+ pericytes56 formed a typical vascular region in OCCC. Tumor cores in recurrent OCCC displayed minimal immune cell infiltration, in contrast to higher concentrations of immune cells in the stromal and perivascular regions. VEGF inhibition led to decreased stroma density, increased T cell infiltration, and reactivation of T cells towards IFNγ signaling and cytotoxicity. To obtain real-world level evidences, we reviewed a series of clinical cases. Combining VEGFi with anti-PD-1 in OCCC resulted in both biochemical and radiological response in patients with OCCC with persistent, recurrent, and metastatic disease, regardless of ARID1A status, thereby expanding the treatment options for OCCC.
Besides using VEGF inhibition to sensitize OCCC immunotherapy, our data also pinpointed new avenues for OCCC therapeutics. The longitudinal analysis of OCCC immune contexture provided an evolutionary perspective of how OCCC dynamically manipulated the host environment to evade immune killing. Recurrent tumor cells showed downregulation of HLA class I molecule but upregulation of the anti-phagocytosis receptor CD47. Antigen-presenting cells including classic DCs and IDO1+ CXCL10+ M1-polarized macrophages were dramatically reduced within recurrent tumors. Contrastingly, the non-immune stromal compartments consisted by fibroblasts, endothelial cells and pericytes were enriched in recurrent tumors and showed enhanced intercellular regulatory potential. The myofibroblast emerged as a central player in OCCC, exerting influence on both tumor phenotype through the CD47/THBS1 axis and T cell function through the CXCL12/CXCR4 pathway. These findings placed strategies to target stroma cells as future directions for OCCC treatment.
Our data also highlighted metabolic reprogramming57 as a targetable hallmark of OCCC. OCCC tumors displayed higher hypoxic metabolic activity compared with serous tumors, with recurrent OCCC tumors showing a shift from glycolysis to FA metabolism, marked by increased CD47 and CD36 expression. Additionally, CD4+ T cell in recurrent OCCC showed restricted diversity, with two main clusters representing naive and hypermetabolic states, suggesting reduced functionality due to microenvironmental hypoxia. Accordantly, a previous study suggested hypoxia transcriptional signature negatively correlated with clinical outcomes of patients with OC who received anti-PD-1 treatment.22 Whether targeting FA metabolism with statins or other novel agents could reverse chemoresistance demands further investigations, as several clinical trials are ongoing regarding the use of statins in OC (NCT06468254, NCT06468254, NCT03532139).
Finally, scFFPE-seq has revolutionized the landscape of single cell technology, and may promote the research of rare diseases. While existing studies have inferred the consistency between fresh cell scRNA-seq and scFFPE-seq using short-term preserved samples (<1 year), they also raise questions about how scFFPE performs with older FFPE blocks after years of preservation.58 Our data indicated that scFFPE-seq was capable of depicting both the cell composition and functional phenotype using FFPE blocks preserved for 2– 7 years, enabling the use of scFFPE-seq to analyze samples with long-term follow-up data. Notably, scFFPE-seq was more efficient in capturing epithelial cells and fibroblasts, and uncovered a novel subset of type II collagen secreting fibroblasts.
Certain limitations were presented in the current study. First, our general sample size is relatively small. Although we have employed various analytical methodologies and external datasets to enhance the validity and robustness of our findings, we were unable to achieve statistical significance in certain aspects of the analysis. Further research in larger patient with OCCC population is essential to deepen our understanding of the immune microenvironment in OCCC. Additionally, this study was primarily computational and correlational, leaving several biological assumptions to be experimentally and interventionally validated. Given the current lack of efficient disease models for OCCC research, the establishment of cancer models, particularly immune-competent patient-derived organoids or xenografts, could be crucial.
Taken together, this study provided a single-cell immune atlas of OCCC, assessing clinical and genomic factors in influencing OCCC immune contexture and tracing the dynamics during tumor recurrence. Although OCCC acquired chemoresistance and immunosuppression through metabolic reprogramming and stromal remodeling, VEGF inhibition reactivated anti-tumor immunity and potentiated immunotherapy.
Data availability statement
Data are available upon reasonable request. The data generated in this study are available within the article and supplementary data files. Other data or materials are available from the corresponding author on reasonable request.
Ethics statements
Patient consent for publication
Ethics approval
The study was designed and conducted according to the Declaration of Helsinki. Human specimen collection and utilization was approved by the ethic committee of Fudan University Shanghai Cancer Center (FUSCC) (No. 050432-4-2108*). Participants gave informed consent to participate in the study before taking part.
Acknowledgments
The authors wish to express their gratitude to all study participants for their cooperation and time. Figure 1A is partially modified from Servier Medical Art (http://smart.servier.com/), licensed under a Creative Common Attribution 3.0 Generic License (https://creativecommons.org/licenses/by/3.0/).
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
SX and LC contributed equally.
Contributors JiajiaL, XC and XW conceived and supervised the study. SX and LC designed the experiments, acquired and analyzed the data, and wrote the manuscript. MY, JianaL, JC, FX, MN and CL participated in sample collection, experiments conduction and revised the manuscript. All authors have read and approved the final manuscript. Dr JiajiaL is designated as the lead contact and guarantor for the authors’ contributions.
Funding This study was granted by China National Natural Science Foundation (NSFC 82272898; 82303660).
Competing interests No, there are no competing interests.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.