A comprehensive molecular analysis of 113 primary ovarian clear cell carcinomas reveals common therapeutically significant aberrations
Diagnostic Pathology volume 18, Article number: 72 (2023)
Molecular aberrations occurring in primary ovarian clear cell carcinoma (OCCC) can be of diagnostic, predictive, and prognostic significance. However, a complex molecular study including genomic and transcriptomic analysis of large number of OCCC has been lacking.
113 pathologically confirmed primary OCCCs were analyzed using capture DNA NGS (100 cases; 727 solid cancer related genes) and RNA-Seq (105 cases; 147 genes) in order to describe spectra and frequency of genomic and transcriptomic alterations, as well as their prognostic and predictive significance.
The most frequent mutations were detected in genes ARID1A, PIK3CA, TERTp, KRAS, TP53, ATM, PPP2R1A, NF1, PTEN, and POLE (51,47,27,18,13,10,7,6,6, and 4%, respectively). TMB-High cases were detected in 9% of cases. Cases with POLEmut and/or MSI-High had better relapse-free survival. RNA-Seq revealed gene fusions in 14/105 (13%) cases, and heterogeneous expression pattern. The majority of gene fusions affected tyrosine kinase receptors (6/14; four of those were MET fusions) or DNA repair genes (2/14). Based on the mRNA expression pattern, a cluster of 12 OCCCs characterized by overexpression of tyrosine kinase receptors (TKRs) AKT3, CTNNB1, DDR2, JAK2, KIT, or PDGFRA (p < 0.00001) was identified.
The current work has elucidated the complex genomic and transcriptomic molecular hallmarks of primary OCCCs. Our results confirmed the favorable outcomes of POLEmut and MSI-High OCCC. Moreover, the molecular landscape of OCCC revealed several potential therapeutical targets. Molecular testing can provide the potential for targeted therapy in patients with recurrent or metastatic tumors.
Ovarian clear cell carcinoma (OCCC) accounts for about 10% of ovarian carcinomas. Its occurrence is geographically different – in Europe and North America it has been estimated between 5 and 13%, whereas in Asia it makes up to 25% of all ovarian carcinomas [1, 2]. OCCC differs by pathological and molecular characteristics from other ovarian carcinomas, including endometrioid (EC), high grade serous (HGSC), low grade serous (LGSC), and mucinous carcinoma (MC).
Molecular aberrations occurring in OCCC can be of predictive and prognostic significance, however, a complex molecular study including genomic and transcriptomic analysis of large number of OCCC has been missing. Mutations occurring in OCCC have been analyzed in a limited number of studies [1, 3,4,5,6,7,8,9,10,11]. Some of these studies have, however, some limitations related mostly to the limited spectrum of genes analyzed, the sensitivity of the method (low coverage of WGS or WES), or a small sample set of analyzed cases. Nevertheless, it has been shown that mutations of ARID1A, PIK3CA, and TERT promoter are a common finding, followed by KRAS, TP53, ATM, and PPP2R1A mutations. The prognostic relevance of stratifying OCCC into different molecular subtypes has been suggested [3, 4, 12]. In one study, 421 OCCC were stratified into two main subgroups. One subgroup (“classic OCCC”) included tumors with ARID1A and other common mutations (such as PIK3CA and TERT) which represented about 83% of tumors. The second subgroup (“HGSC-like”) was characterized by the TP53 mutation and showed an enriched expression of genes involved in extracellular matrix organization, mesenchymal differentiation, and immune-related pathways . However, as admitted by the authors, the main weakness of this study is that some cases in the “HGSC-like” subgroup are probably true HGSC misclassified as OCCC. Another study of 55 OCCC suggested four different molecular subgroups: “PIK3CA”, “ARID1A”, “PIK3CA-ARID1A” and “Undetermined” . Another possible approach is the stratification of OCCC into the Cancer Genome Atlas (TCGA)-based molecular subtypes used for endometrial carcinomas. Using this approach, POLE mutated (POLEmut) or mismatch repair deficient (MMR-D) cases have a better prognosis than p53 abnormal (p53abn) cases or cases with no specific molecular profile (NSMP) subgroup . Better prognosis for MMR-D OCCC has been suggested also by the results of another study .
Concerning gene rearrangements, only two studies on a limited number of OCCC have focused on this topic to date [14, 15]. However, their results are problematic, given that in one of those studies (analyzing 4 OCCC) one tumor with detected fusion CCNY::NRG4 actually represented a metastasis from primary uterine mixed clear cell and endometrioid carcinoma . The second study analyzed 20 OCCCs in which three cases with several fusions were found, but the authors admitted the possibility of a methodological problem in analyzing the data and due to this their results are equivocal at best .
Concerning transcriptome analysis, only four studies performed expression profiling using RNA sequencing (RNA-Seq). One afore-mentioned large study combined a transcriptomic study with targeted DNA sequencing . Another three studies included only a limited number of cases, namely 11, 19, and 6, respectively [16,17,18]. Two of them focused on the differences among ovarian carcinoma subtypes. The first study suggested the possible significance of expression analysis for differential diagnosis between OCCC and clear cell carcinoma of the uterus . The second study examined these differences between OCCC, HGSC, and EC, while the third study mainly focused on the differences between OCCC and HGSC [17, 18]. Moreover, two other studies which used different methodological approaches (microarray analysis) included 37 and 8 OCCCs, respectively [19, 20].
Recently, epigenome profiling in OCCC has been performed in a collaborative study including 271 cases from ten study sites and genome-wide tumor DNA methylation profiling . Their analyses supported the involvement of immune related pathways in OCCC and brought insight into epigenomic profiling. Furthermore, they revealed a higher level of chromosomal aneuploidy in OCCCs with ARID1A/PIK3CA mutation .
The primary goal of our study was to perform a comprehensive genomic and transcriptomic analysis of a well-defined sample set of 113 primary OCCCs with the aim to characterize these tumors with respect to the occurrence of molecular aberrations, as well as their prognostic and predictive value. Secondly, we focused on the possible stratification of our sample set into molecularly defined subgroups based on mutation and/or mRNA expression pattern.
The archives of the participating pathology departments were searched for cases originally diagnosed as OCCC. All cases were carefully reviewed by two pathologists (PD and MKB) and only cases meeting the strict morphological and immunohistochemical criteria were included into the study. The criteria included the morphology of clear cell carcinoma as defined elsewhere, associated with immunohistochemical profile compatible with the diagnosis, including positivity of PAX8 and at least one marker of “clear cell” differentiation (HNF1B, AMACR or napsin A) together with negativity of WT1 [22,23,24]. Finally, 120 OCCC cases were selected for DNA and RNA isolation. The clinicopathological characteristics of these cases are summarized in Table 1.
Next generation sequencing (NGS) analyses
The isolation of nucleic acids from FFPE tumor tissue for further DNA NGS and RNA-Seq analyses was performed as described previously . Out of the 120 cases included into the study, 113 OCCCs were eligible for DNA and/or RNA NGS. Samples insufficient for complex molecular analyses were excluded (20/120 DNAs; 17%, and 15/120 RNAs; 13%, respectively), 92 samples had both complete DNA and RNA NGS analyses.
DNA NGS analysis
Sequence capture NGS analysis of DNA was performed for 100 qualitatively sufficient OCCC cases in order to assess mutation pattern and frequency, tumor mutation burden (TMB), and microsatellite instability (MSI). Copy number variation (CNV) analyses were not performed because of the low DNA quality not suitable for reliable CNV testing.
The library preparation was performed using the KAPA HyperPlus kit [according to KAPA HyperCap Workflow v3.0 (Roche, Basel, Switzerland)] and a panel of hybridization probes against multiple targets of cancer relevant genes (727 genes or gene parts; 2097 kbp of the target sequence including 1708 kbp of coding regions; Roche; Supplementary Information). The prepared sample libraries were pair-end sequenced by the NextSeq 500 instrument (Illumina, San Diego, California) using the NextSeq 500/550 High Output Kit v2.5 (Illumina). The biostatistical evaluation using NextGENe software (Softgenetics) and the interpretation of DNA variants was performed as follows.
All the frameshift, no-start, and no-stop splice variants in the consensus splice sites, nonsense variants, and missense variants known as pathogenic and/or likely pathogenic (class 4/5 mutations; except for one nonsense POLE mutation which does not cause ultramutated phenotype) according to the ClinVar database were considered as deleterious. Detailed pipelines of all the NGS data analysis together with the module settings are available upon request. The analysis does not allow for the distinction between somatic and germline variants. The TP53 variants were classified according to https://p53.iarc.fr/, ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/), and https://www.cancerhotspots.org/. The size of the sequenced panel enabled us to also assess the TMB (number of mutations per 1 Mega base; mut/Mb); however, TMB was calculated only for samples with ≥ 40% tumor cells. Samples with TMB ≥ 10 mut/Mb were considered TMB-High. All synonymous and nonsynonymous variants with an allele frequency of ≥ 10% were counted. Furthermore, potential germline variants [according to databases of known germline polymorphisms including the Single Nucleotide Polymorphism database (dbSNP) and Exome Aggregation Consortium (ExAC)] and known or probable driver mutations (according to the COSMIC and ClinVar database) were determined. The resulting mutation number was normalized to 1 Mb. MSI was evaluated from NGS data using CLC Genomics Workbench software (CLC GW; Qiagen, Venlo, The Netherlands) and a module ‘Detect MSI status’ with default settings. When more than 20% of the 17 evaluated microsatellite markers were unstable, then the sample was considered microsatellite unstable (detailed analysis of the microsatellite status including a comparison of its assessment using several methods will be part of a forthcoming study). We also evaluated the hot spot variants in the TERT promoter (TERTp) that are clinically relevant (c.-124 C > T rs1561215364, c.-124 C > A, and c.–146 C > T rs1561215364).
RNA NGS analysis
Total RNA samples were processed according to the KAPA RNA HyperPrep Kit protocol (Roche; input 300 ng where available; denaturation/fragmentation 85 °C – 2 min; 11 cycles of PCR). In those samples of sufficient quality (n = 105), the target sequences were enriched by the standard KAPA HyperCap Workflow v3 (Roche) using a custom panel focused on the pan-cancer markers and potential fusion genes (147 genes; 373 kbp of the target DNA sequence; Roche; Supplementary Information). The final libraries were pair-end sequenced by the NextSeq 500 instrument using 300 cycles chemistry kits (Illumina) with a target of 10 million single reads.
The sequencing data were analyzed using the CLC GW v21.0.5. (Qiagen) by an in-house pipeline which includes targeted RNA-Seq expression analysis (RNA-Seq Analysis module) and detection of fusion genes (Detect and Refine Fusion Genes module). The bioinformatics pipeline and module settings are available upon request.
All fusions identified by the CLC GW were manually checked, filtered, and confirmed using IGV v2.11.3. (Broad Institute, Cambridge, Massachusetts). Only those fusions meeting the following criteria were considered as true fusions: (i) fusions involving protein-coding genes with standard exon-exon junctions (± 15 bp range from exon boundary), with substantial expression when compared to other samples at the respective region, (ii) ≥ 10% of reads supporting fusion presence (crossing reads) out of read counts at the respective location. Frequently repeated fusions, fusions of genes from the same gene family, or transcriptional readthroughs were excluded and considered as artefacts. The nomenclature of the detected mutations follows HGVS recommendations (https://varnomen.hgvs.org/). Each fusion that was not described in the Quiver database (http://quiver.archerdx.com/), Fusion GDB (https://ccsm.uth.edu/FusionGDB/index.html), Mitelmandatabase (https://mitelmandatabase.isb-cgc.org/) or literature was considered as novel. OncoKB (https://www.oncokb.org/) was searched for the clinical significance of the detected fusions with respect to therapeutic actionability.
Expression profiling was performed using the Heat Map for RNA-Seq module in CLC GW with Manhattan distance and Complete linkage settings.
Hierarchical clustering based on the two-dimensional heat map of normalized expression values (RPKM) of all 147 gene targets in the panel was performed using the Heat Map for RNA-Seq module in CLC GW with Manhattan distance and Complete linkage settings.
Normalization of mRNA expression was evaluated as RPKM (reads per kilobase of transcript per million reads mapped) and VCP, SF3B1, and ATP5F1B genes were used as a reference.
The software R (version 4.0.2, https://www.r-project.org/) was used to perform the statistical analyses.
Group comparisons were performed for categorical variables using Pearson chi-squared test or Fisher Exact test, based on excepted values.
Time-to-event analyses were computed using the package ‘survival’ and ‘survminer’. The associations between patient characteristics and survival were evaluated using univariable and multivariable Cox proportional hazard models (Cox PH) and described by the regression coefficient (β), hazard ratio (HR) with 95% confidence intervals (CI) and statistical significance. P-values were obtained by applying Wald statistics based on Cox PH model. Variables that showed significance in univariable model were selected for the multivariable Cox PH full model. A backward stepwise elimination was used to reach minimal adequate model. Survival curves were constructed by the Kaplan-Meier method based on the stratified risk score and compared by means of the log-rank test and likelihood ratio (LR) test.
Recurrence-free survival (RFS) was defined as the period from the date of diagnosis to the date of recurrence or death of disease, overall survival (OS) was defined as the period from the date of diagnosis to the date of recorded death. The length of the follow-up period (FU) was calculated from the date of diagnosis to the last recorded follow-up visit or death of patient. Among 87 patients included in the analysis, 21 died (24%), 10 of them from other cause than disease related.
All tests were two-sided and p-values < 0.05 were considered significant. Hierarchical clustering is based on the TMM method (trimmed mean of M values) expression data normalization, with additional calculation of TMM-adjusted log CPM (Counts per Million) for each gene with a subsequent Z-score normalization across all samples for each gene (counts for each gene are mean centered and scaled to unit variance).
Differential expression in two groups module, which is implemented in CLC GW, was used for the analysis of differential expression of individual genes included in panel between samples hierarchically clustered in group 2 versus group 1. This module is a multi-factorial statistics test based on a negative binomial Generalized Linear Model (the statistical model for this module is thoroughly described in the CLC GW manual - https://digitalinsights.qiagen.com/technical-support/manuals/).
Genomic DNA alterations
In total, the DNA NGS analysis of 100 eligible OCCC cases revealed pathogenic or likely pathogenic (class 4/5) mutations in 218 of the 727 analyzed genes (Supplementary Table 2). Out of those, for 144 genes class 4/5 mutation was detected only once in the whole sample set. The mean and median number of the detected somatic class 4/5 mutations per sample was 5.4 and 3, respectively (range 0–58). The spectrum and frequency of the mutated genes (only those affected in at least 2% of cases) in the context of clinicopathological variables and expression and fusion analysis is shown in Fig. 1. Table of genes with class 4/5 mutations in respective cases is provided in Supplementary Table 2.
The most frequently altered genes by class 4/5 mutations were ARID1A in 52/100 (51%), PIK3CA in 47/100 (47%), TERTp (27/100; 27%), KRAS 18/101 (18%), TP53 13/100 (13%), and ATM 10/100 (10%) cases. 4/100 (4%) cases showed POLE mutation, and 7/100 (7%) were MSI-High cases (one of which was also POLEmut).
The median tumor mutation burden for all 100 OCCCs was 3 mutations/megabase (TMB = 3 mut/Mb; range 0–86). Nine TMB-High OCCCs had median TMB = 25 mut/Mb (TMB range 10–86), of which 4 were POLEmut (TMB range 24–86), 4 microsatellite instable (MSI-High; TMB range 14–26), and one sample was microsatellite stable (MSS) and without POLE mutation (POLEwt).
We correlated the mutation status in frequently mutated genes with clinico-pathological characteristics (Supplementary Table 3). In our data set we observed that mutations in ARID1A or in PIK3CA correlated with a younger age of patients (p = 0.043 or p = 0.002, respectively). The same was observed for group of POLEmut and/or MSI-High patients compared with patients with POLEwt and MSS (p = 0.015). The KRAS mutations were more common in the early stages of the disease (p = 0.020). The KRAS and TP53 mutations were mutually exclusive (see Fig. 1).
RNA-Seq – fusions transcript analyses
The RNA-Seq was successful in 105 cases. Only fusions involving different protein-coding genes in ≥ 10% of reads were reported. 13 different gene fusions were identified in 14/105 (13%) of the OCCCs. Of those, six were novel (LAMB1::MET, CCDC3::AKT3, NBN::CNGB3, DHX15::ATR, LGALS3::EZH2, ARID1A::MIA2) while seven have already been described in the databases or literature (Table 2). Two fusions were recurrent, namely TFG::ADGRG7 (2/105; 2%), and MET fusions with different fusion partners (4/105; 4%). MET fusions included ST7::MET (2/105; each with different breakpoints), CAPZA2::MET, and LAMB1::MET (for fusion details see Supplementary Fig. 1 and Table 3). The 14 tumors with detected fusions did not show any specific morphologic parameters which would distinguish them from those without fusions (data not shown).
9/92 cases with complete DNA and RNA analysis had a gene fusion detected (Fig. 1). Out of those, 6 had a concurrent aberration in at least one of the considered drivers ARID1A, PIK3CA, or KRAS. One case with MET::ST7 fusion had a concurrent mutation in PPP2R1A, one case with fusion NBN::CNGB3 had a concurrent mutation in ARHGAP35, and one case with TGF::ADGRG7 fusion had a concurrent mutation in KDM5C. None of the nine cases with a gene fusion had a concurrent hot-spot alteration in the TERT promoter.
RNA-Seq – expression profile
The expression analysis of 147 genes was possible in 105 cases and revealed a heterogenous expression pattern. Unsupervised hierarchical clustering suggested two main clusters (Fig. 2). The detailed analysis of additionally normalized expression data to reference genes (VCP, SF3B1, ATP51B) comparing cluster 1 (93 cases) with cluster 2 (12 cases) is in Supplementary Table 5. On the mRNA level cluster 2 shows higher expression (≥ 2fold change) of AKT3, DDR2, CTNNB1, JAK2, KIT or PDGFRA in a majority of the samples and a lower expression of CDH1, ERBB2, ERBB3, FGFR3, HIST1H3B, HNF1B and POLQ in all samples when compared with cluster 1 (p < 0.00001).
Concerning the morphological features, cluster 2 was associated with macronucleoli (p = 0.031), increased mitotic activity (p = 0.014), and less necrosis (p < 0.001; Supplementary Table 4). We did not find any significant correlations of the expression pattern with regards to age, FIGO, or patient outcomes (Supplementary Table 3).
The follow-up data was available for 87 cases out of the 100 OCCCs. The median follow-up period was 3 years, mean ± SD was 3.5 ± 0.3 years. Cases which were POLEmut and/or MSI-High (n = 9) had median follow-up 42 months, while compared with median 33 months in cases with POLEwt and/or MSS. Based on univariate analysis, POLEmut and/or MSI-High cases had favorable RFS (p = 0.04; no case with POLEmut or MSI-High had an event compared with POLEwt/MSS cases). However, multivariable model adding age, tumor stage, radicality of surgical resection, and adjuvant therapy into consideration showed only the radicality of surgical resection as significant prognostic marker, indicating that cases with reached R0 resection margin correlated with strong decreased risk of relapse. No other stratification into subgroups based on the molecular features, including classification defined by TCGA for endometrial carcinoma (POLEmut, MSI-High, p53abn, and NSMP), classification based on the presence of ARID1A, PIK3CA, TERT, TP53 or KRAS mutation, and classification based on expression profiling, showed prognostic significance.
The molecular features of OCCC can be of prognostic and predictive significance, but the current literary data is limited and equivocal. The results of previous studies showed a wide range of mutation frequency of several genes. An explanation for these differences could include the small sample sets in some previous studies, differences in methodology, or bias caused by the inclusion of tumors of other histogenesis, especially EC and HGSC with clear cell change. In our study, only tumors meeting strict diagnostic criteria were included, which prevented bias caused by the inclusion of tumors of other histogenesis. We revealed 18 recurrently mutated targets (detected in ≥ 4 patients), namely ARID1A (52/100; 52%); PIK3CA (47/100; 47%); TERT promoter hot-spot mutations (27/100; 27%); KRAS (18/100; 18%; with one case carrying the G12C mutation); TP53 (13/100; 13%); ATM (10/100; 10%); PPP2R1A (7/100; 7%); NF1, PTEN (6/100; 6%); ARHGAP35, ARID1B, CREBBP, PIK3R1, RASA1, SHPRH (5/100; 5%); LATS1, MLL2, POLE, SETX (4/100; 4%). The majority of altered genes code for proteins which are involved in the PI3K/AKT and/or RAS/MAPK signaling pathways (KRAS, NF1, PIK3CA, PIK3R1, PTEN, RASA1), DNA repair and cell cycle regulation (ATM, LATS1, MLL2, POLE, PPP2R1A, SETX, SHPRH, TP53), and/or chromatin remodeling (ARID1A, ARID1B, CREBBP, KMT2D). Interestingly, in our study on Caucasian patients we did not detect any BRAF mutation, compared to 2/48 (4%) described in one study performed on the Japanese population . However, our finding is in concordance with the results of others [3, 4].
We identified the POLE mutation in 4/100 (4%) of cases, which is discordant to the previous 9 studies analyzing POLE mutations in OCCC (Supplementary Table 1), in which POLE mutation was detected in only 2 of 405 cases. Of note, another case in our sample set carried a nonsense POLE mutation NM_006231.2:c.3961 A > T, p.(R1321X) which is currently considered as a variant of uncertain significance (VUS). In one study, the authors classified the tumors into TCGA-based molecular subtypes, including POLEmut, MMR-D, p53abn, and NSMP. When this approach is used in our study, the results are similar for NSMP tumors (76% in both ours and their cohort), but different for other subtypes (POLEmut 4% in our study vs. 0.9% in theirs, MMR-D 7% in our study vs. 3.5%, and p53abn 13% in our study vs. 20%) . In our study, OCCC cases with POLEmut and/or MSI-High status had a better prognosis (favorable RFS, p = 0.041) compared to POLEwt /MSS primary OCCC cases using univariate analysis, which is in concordance with recently published data . However, the prognostic significance of this molecular trait was not confirmed in multivariable analysis in the context of other prognostic factors (age, stage at diagnosis, reaching R0 surgical treatment, and adjuvant therapy) where strongest prognostic factor was R0 resection margin.
Unlike the results reported in two previous studies, we did not observe worse survival in TP53 mutated OCCC (our study showed a trend towards favorable RFS and OS for TP53mut cases compared to other cases).
Another possible prognostic molecular marker suggested by one study is a mutation in the promoter region of TERT . However, in our study the presence of TERTp hot-spot mutations was not associated with worse survival.
The co-occurrence of mutations in ARID1A, PIK3CA, PIK3R1, and KRAS was common among primary OCCC. KRAS and TP53 mutations in our study were mutually exclusive (except for one MSI-High tumor). Our results are in disagreement with another study in which mutual exclusivity between somatic mutations of ARID1A, TP53, PIK3CA, and PIK3R1 was found. In this study statistically significant co-occurrence between mutations in ARID1A, PIK3CA, or TP53 and BRCA1/BRCA2 was also found . In our cohort, BRCA1 (1/100) or BRCA2 (2/100) mutations were very rare. Such discrepancies could be once again explained by the inclusion of cases of other histogenesis, as acknowledged by the authors, such as HGSC with clear cell change, which would also explain the higher frequency of TP53 mutations in their study .
In our study, 20/100 cases (20%) showed no mutation in any of the commonly mutated genes, including ARID1A, PIK3CA, TERT, KRAS, TP53, ATM, and NF1. RNA analysis of these cases showed gene fusion in 3/19 cases (16%). Of the other cases which did show mutation in some of the commonly mutated genes, gene fusion was detected less frequently (6/73; 8% cases). This OCCC group carried none or a maximum of three genomic alterations (mainly involved in DNA repair, chromatin remodeling, or PI3K/AKT signaling). Altogether, we found gene fusions in 14/105 cases. A majority of fusions were intrachromosomal rearrangements, only 4/13 resulted from interchromosomal translocation. In 5/14 (36%) cases, the 3´ gene partner was a tyrosine kinase receptor. This finding can be of clinical significance, given that several tyrosine kinase inhibitors are currently available in clinical practice. In total, we have found MET fusions in 4/105 (4%) OCCCs. Two of them were ST7::MET fusions which have been previously reported in 1 HGSC and 1 OCCC, each with different breakpoints, which according to the database OncoKB [ ; accessed September 30, 2022] were considered to be likely oncogenic . Another fusion, CAPZA2::MET has also been previously described in the Mitelman database of Chromosome Aberrations and Gene Fusions in Cancer (https://mitelmandatabase.isb-cgc.org; accessed 16, Oct 2022). Finally, we identified one novel LAMB1::MET fusion. The tyrosine kinase domain of MET was preserved in all reported fusions. Interestingly, a substantial part of the fusions affect the tyrosine kinase receptor molecule (6/13; MET fusions; ERBB4::IKZF2 and CCDC3::AKT3). Novel fusions NBN::CNGB3 and DHX15::ATR affect genes involved in DNA repair, ESR1::ARMT1 affects genes involved in DNA damage response and PI3K/AKT signaling regulation.
The mRNA expression pattern was heterogeneous, but unsupervised hierarchical clustering analysis revealed two main clusters. Cluster 1 included 93 cases with variable expression pattern, while cluster 2 included 12 cases characterized by higher expression of AKT3, DDR2, CTNNB1, JAK2, KIT or PDGFRA, and lower expression of CDH1, ERBB2, ERBB3, FGFR3, HIST1H3B (also known as H3C2), HNF1B and POLQ when compared with cluster 1 (p < 0.001). Our results are not comparable with other studies focusing on expression profiling due to the very limited overlap in genes analyzed [3, 17].
Currently, there are no targeted therapies specific for OCCC, although a number of candidate targets have recently been identified and reviewed, including targets in PI3K/AKT/mTOR pathway, and ARID1A deficiency [29, 30]. Nevertheless, some molecular aberrations such as NTRK fusions, microsatellite instability, and high tumor mutation burden have predictive value in all solid tumors. In our study, the vast majority of altered genes are involved in PI3K/AKT and/or RAS/MAPK signaling pathways, DNA repair and cell cycle regulation, chromatin remodeling, or cell adhesion. The ARID1A mutated tumors could be potentially targetable by PARP inhibitors or HDAC inhibitors [31,32,33,34]. The presence of the PIK3CA mutation might indicate susceptibility to PI3K- or mTOR-inhibitors. Targets affecting the mTOR pathway have been detected in OCCC (including PIK3CA, PTEN, FBXW7, PIK3R1, AKT3, NF1). TMB-High (detected in 9% of cases) tumors may be eligible for immune therapy. The BRCA1 and BRCA2 mutations were very rare events, as well as MSI-High cases; however, those patients could be considered for treatment by PARP inhibitors or immunotherapy. Importantly, we identified frequently altered genes coding tyrosine kinase receptors, such as MET, AKT3, DDR2, CTNNB1, KIT, or PDGFRA, which are known therapeutic targets . The MET inhibitors are currently considered as predictive biomarker in patients with lung adenocarcinoma  and c-MET has been considered as a potential therapeutic target in OCCC . Another potential target could also be the overexpression of non-receptor tyrosine kinase JAK2 .
We are aware of the limitations of our study. The main limitation is that CNV and epigenetic changes were not performed. Moreover, analysis of the outcomes can be influenced by a limited number of events in our cohort.
Knowledge about the molecular landscape of OCCC is gaining significance due to the expanding possibilities of targeted treatments available for several tumors. Our study described for the first time a complex molecular pattern in a pathologically well-defined sample set of 113 OCCC fulfilling strict inclusion criteria. Our results confirmed the favorable outcomes of POLEmut and MSI-High OCCC. Moreover, the molecular landscape of OCCC revealed several potential therapeutical targets including MET fusion, and molecular testing can provide the potential for targeted therapy in patients with recurrent or metastatic tumors.
All data supporting the findings of this study are available within the paper and its Supplementary Information. The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.
Ovarian clear cell carcinoma
High grade serous carcinoma
Low grade serous carcinoma
Formalin-fixed and paraffin-embedded
Local recurrence-free survival
Read per kilobase million
No specific molecular profile
Mismatch repair deficient
- class 4/5 mutation:
Likely pathogenic or pathogenic mutation
Friedlander ML, Russell K, Millis S, Gatalica Z, Bender R, Voss A. Molecular Profiling of Clear Cell ovarian cancers: identifying potential treatment targets for clinical trials. Int J Gynecol Cancer. 2016;26:648–54. https://doi.org/10.1097/IGC.0000000000000677.
Zhu C, Xu Z, Zhang T, Qian L, Xiao W, Wei H, et al. Updates of Pathogenesis, Diagnostic and therapeutic perspectives for ovarian Clear Cell Carcinoma. J Cancer. 2021;12:2295–316. https://doi.org/10.7150/jca.53395.
Bolton KL, Chen D, de la Corona RI, Fu Z, Murali R, M, KB, et al. Molecular subclasses of clear cell ovarian carcinoma and their impact on disease behavior and outcomes. Clin Cancer Res. 2022. https://doi.org/10.1158/1078-0432.CCR-21-3817.
Oliveira D, Schnack TH, Poulsen TS, Christiansen AP, Hogdall CK, Hogdall EV. Genomic sub-classification of Ovarian Clear Cell Carcinoma revealed by distinct mutational signatures. Cancers (Basel). 2021;13. https://doi.org/10.3390/cancers13205242.
Liu H, Zhang Z, Chen L, Pang J, Wu H, Liang Z. Next-generation sequencing reveals a very low prevalence of deleterious mutations of homologous recombination repair genes and homologous recombination Deficiency in Ovarian Clear Cell Carcinoma. Front Oncol. 2021;11:798173. https://doi.org/10.3389/fonc.2021.798173.
Shibuya Y, Tokunaga H, Saito S, Shimokawa K, Katsuoka F, Bin L, et al. Identification of somatic genetic alterations in ovarian clear cell carcinoma with next generation sequencing. Genes Chromosomes Cancer. 2018;57:51–60. https://doi.org/10.1002/gcc.22507.
Itamochi H, Oishi T, Oumi N, Takeuchi S, Yoshihara K, Mikami M, et al. Whole-genome sequencing revealed novel prognostic biomarkers and promising targets for therapy of ovarian clear cell carcinoma. Br J Cancer. 2017;117:717–24. https://doi.org/10.1038/bjc.2017.228.
Maru Y, Tanaka N, Ohira M, Itami M, Hippo Y, Nagase H. Identification of novel mutations in japanese ovarian clear cell carcinoma patients using optimized targeted NGS for clinical diagnosis. Gynecol Oncol. 2017;144:377–83. https://doi.org/10.1016/j.ygyno.2016.11.045.
Murakami R, Matsumura N, Brown JB, Higasa K, Tsutsumi T, Kamada M, et al. Exome sequencing Landscape Analysis in Ovarian Clear Cell Carcinoma Shed Light on Key chromosomal regions and mutation gene networks. Am J Pathol. 2017;187:2246–58. https://doi.org/10.1016/j.ajpath.2017.06.012.
Sugino K, Tamura R, Nakaoka H, Yachida N, Yamaguchi M, Mori Y, et al. Germline and somatic mutations of homologous recombination-associated genes in japanese ovarian cancer patients. Sci Rep. 2019;9:17808. https://doi.org/10.1038/s41598-019-54116-y.
Yin X, Bi R, Ma P, Zhang S, Zhang Y, Sun Y, et al. Multiregion whole-genome sequencing depicts intratumour heterogeneity and punctuated evolution in ovarian clear cell carcinoma. J Med Genet. 2020;57:605–9. https://doi.org/10.1136/jmedgenet-2019-106418.
Simila-Maarala J, Soovares P, Pasanen A, Ahvenainen T, Vahteristo P, Butzow R, et al. TCGA molecular classification in endometriosis-associated ovarian carcinomas: novel data on clear cell carcinoma. Gynecol Oncol. 2022;165:577–84. https://doi.org/10.1016/j.ygyno.2022.03.016.
Stewart CJ, Bowtell DD, Doherty DA, Leung YC. Long-term survival of patients with mismatch repair protein-deficient, high-stage ovarian clear cell carcinoma. Histopathology. 2017;70:309–13. https://doi.org/10.1111/his.13040.
Agostini A, Brunetti M, Davidson B, Goran Trope C, Heim S, Panagopoulos I, et al. Identification of novel cyclin gene fusion transcripts in endometrioid ovarian carcinomas. Int J Cancer. 2018;143:1379–87. https://doi.org/10.1002/ijc.31418.
Earp MA, Raghavan R, Li Q, Dai J, Winham SJ, Cunningham JM, et al. Characterization of fusion genes in common and rare epithelial ovarian cancer histologic subtypes. Oncotarget. 2017;8:46891–9. https://doi.org/10.18632/oncotarget.16781.
Alldredge J, Randall L, De Robles G, Agrawal A, Mercola D, Liu M, et al. Transcriptome analysis of ovarian and uterine clear cell malignancies. Front Oncol. 2020;10:598579. https://doi.org/10.3389/fonc.2020.598579.
Fridley BL, Dai J, Raghavan R, Li Q, Winham SJ, Hou X, et al. Transcriptomic characterization of Endometrioid, clear cell, and High-Grade Serous Epithelial Ovarian Carcinoma. Cancer Epidemiol Biomarkers Prev. 2018;27:1101–9. https://doi.org/10.1158/1055-9965.EPI-17-0728.
Nagasawa S, Ikeda K, Horie-Inoue K, Sato S, Itakura A, Takeda S, et al. Systematic identification of characteristic genes of ovarian Clear Cell Carcinoma compared with High-Grade Serous Carcinoma based on RNA-Sequencing. Int J Mol Sci. 2019;20. https://doi.org/10.3390/ijms20184330.
Stany MP, Vathipadiekal V, Ozbun L, Stone RL, Mok SC, Xue H, et al. Identification of novel therapeutic targets in microdissected clear cell ovarian cancers. PLoS ONE. 2011;6:e21121. https://doi.org/10.1371/journal.pone.0021121.
Winterhoff B, Hamidi H, Wang C, Kalli KR, Fridley BL, Dering J, et al. Molecular classification of high grade endometrioid and clear cell ovarian cancer using TCGA gene expression signatures. Gynecol Oncol. 2016;141:95–100. https://doi.org/10.1016/j.ygyno.2016.02.023.
Cunningham JM, Winham SJ, Wang C, Weiglt B, Fu Z, Armasu SM, et al. DNA methylation profiles of ovarian Clear Cell Carcinoma. Cancer Epidemiol Biomarkers Prev. 2022;31:132–41. https://doi.org/10.1158/1055-9965.EPI-21-0677.
Fadare O, Zhao C, Khabele D, Parkash V, Quick CM, Gwin K, et al. Comparative analysis of Napsin A, alpha-methylacyl-coenzyme A racemase (AMACR, P504S), and hepatocyte nuclear factor 1 beta as diagnostic markers of ovarian clear cell carcinoma: an immunohistochemical study of 279 ovarian tumours. Pathology. 2015;47:105–11. https://doi.org/10.1097/PAT.0000000000000223.
DeLair D, Oliva E, Köbel M, Macias A, Gilks CB, Soslow RA. Morphologic spectrum of immunohistochemically characterized clear cell carcinoma of the ovary: a study of 155 cases. Am J Surg Pathol. 2011;35:36–44. https://doi.org/10.1097/PAS.0b013e3181ff400e.
Köbel M, Kalloger SE, Carrick J, Huntsman D, Asad H, Oliva E, et al. A limited panel of immunomarkers can reliably distinguish between clear cell and high-grade serous carcinoma of the ovary. Am J Surg Pathol. 2009;33:14–21. https://doi.org/10.1097/PAS.0b013e3181788546.
Dundr P, Gregova M, Hojny J, Krkavcova E, Michalkova R, Nemejcova K, et al. Uterine cellular leiomyomas are characterized by common HMGA2 aberrations, followed by chromosome 1p deletion and MED12 mutation: morphological, molecular, and immunohistochemical study of 52 cases. Virchows Arch. 2022;480:281–91. https://doi.org/10.1007/s00428-021-03217-z.
Nishikimi K, Nakagawa K, Tate S, Matsuoka A, Iwamoto M, Kiyokawa T, et al. Uncommon human telomerase reverse transcriptase promoter mutations are Associated with Poor Survival in Ovarian Clear Cell Carcinoma. Am J Clin Pathol. 2018;149:352–61. https://doi.org/10.1093/ajcp/aqx166.
Chakravarty D, Gao JJ, Phillips S, Kundra R, Zhang HX, Wang JJ, et al. OncoKB: a Precision Oncology Knowledge Base. Jco Precis Oncol. 2017;1. https://doi.org/10.1200/PO.17.00011.
Lopez-Nieva P, Fernandez-Navarro P, Grana-Castro O, Andres-Leon E, Santos J, Villa-Morales M et al. Detection of novel fusion-transcripts by RNA-Seq in T-cell lymphoblastic lymphoma. Sci Rep-Uk 2019; 9. https://doi:ARTN 5179.
Khalique S, Lord CJ, Banerjee S, Natrajan R. Translational genomics of ovarian clear cell carcinoma. Semin Cancer Biol. 2020;61:121–31. https://doi.org/10.1016/j.semcancer.2019.10.025.
Amano T, Chano T, Yoshino F, Kimura F, Murakami T. Current position of the molecular therapeutic targets for ovarian Clear Cell Carcinoma: A literature review. Healthc (Basel). 2019;7. https://doi.org/10.3390/healthcare7030094.
Shen JF, Peng Y, Wei LZ, Zhang W, Yang L, Lan L, et al. ARID1A Deficiency impairs the DNA damage checkpoint and sensitizes cells to PARP inhibitors. Cancer Discov. 2015;5:752–67. https://doi.org/10.1158/2159-8290.Cd-14-0849.
Williamson CT, Miller R, Pemberton HN, Jones SE, Campbell J, Konde A, et al. ATR inhibitors as a synthetic lethal therapy for tumours deficient in ARID1A. Nat Commun. 2016;7:13837. https://doi.org/10.1038/ncomms13837.
Fukumoto T, Park PH, Wu S, Fatkhutdinov N, Karakashev S, Nacarelli T, et al. Repurposing Pan-HDAC inhibitors for ARID1A-Mutated ovarian Cancer. Cell Rep. 2018;22:3393–400. https://doi.org/10.1016/j.celrep.2018.03.019.
Bitler BG, Wu S, Park PH, Hai Y, Aird KM, Wang Y, et al. ARID1A-mutated ovarian cancers depend on HDAC6 activity. Nat Cell Biol. 2017;19:962–73. https://doi.org/10.1038/ncb3582.
Lemmon MA, Schlessinger J. Cell signaling by receptor tyrosine kinases. Cell. 2010;141:1117–34. https://doi.org/10.1016/j.cell.2010.06.011.
Plenker D, Bertrand M, de Langen AJ, Riedel R, Lorenz C, Scheel AH, et al. Structural alterations of MET trigger response to MET kinase inhibition in Lung Adenocarcinoma Patients. Clin Cancer Res. 2018;24:1337–43. https://doi.org/10.1158/1078-0432.Ccr-17-3001.
Reddy MM, Deshpande A, Sattler M. Targeting JAK2 in the therapy of myeloproliferative neoplasms. Expert Opin Ther Tar. 2012;16:313–24. https://doi.org/10.1517/14728222.2012.662956.
The authors wish to thank Mgr. Zachary Harold Kane Kendall, B.A. (Institute for History of Medicine and Foreign Languages, First Faculty of Medicine, Charles University) for the English proofreading.
This work was supported by the Ministry of Health, Czech Republic (MH CZ DRO-VFN 64165 and AZV NV19-03-00007), by Charles University (Project UNCE204065, SVV 260631), by the European Regional Development Fund (EF16_013/0001674), and by research project BBMRI_CZ LM2023033).
The authors declare that thay have no competing interests.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the General University Hospital in Prague in compliance with the Helsinki Declaration (No. 2140/19 S-IV). The Ethics Committee waived the requirement for informed consent as according to the Czech Law (Act. no. 373/11, and its amendment Act no. 202/17), it is not necessary to obtain informed consent in fully anonymized studies.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Stružinská, I., Hájková, N., Hojný, J. et al. A comprehensive molecular analysis of 113 primary ovarian clear cell carcinomas reveals common therapeutically significant aberrations. Diagn Pathol 18, 72 (2023). https://doi.org/10.1186/s13000-023-01358-0
- Capture DNA NGS
- Rare ovarian tumors
- POLE mutation