- Research
- Open access
- Published:
Impact of neoadjuvant chemotherapy on somatic mutation status in high-grade serous ovarian carcinoma
Journal of Ovarian Research volume 15, Article number: 50 (2022)
Abstract
Background
Patients treated with neoadjuvant chemotherapy (NACT) for advanced high-grade serous ovarian carcinoma (HGSC) have a higher rate and shorter time to platinum-resistant recurrence compared to patients treated with primary cytoreductive surgery (PCS) and adjuvant chemotherapy. The purpose of this study is to determine the impact of NACT on somatic mutation status in platinum-sensitive and resistant HGSC. Patients with advanced HGSC who had a documented response to platinum-based NACT, a banked blood sample, and a banked tumor sample before and after NACT were identified. Whole exome and/or targeted deep sequencing was performed in matched normal and pre/post-NACT tumor samples from 3 platinum-resistant and 2 platinum-sensitive patients to identify somatic non-synonymous mutations at each time point.
Results
When comparing exonic non-synonymous mutations in pre-NACT and post-NACT samples from the same patient, an average of 41% (1-68%) of genes were mutated at both time points. There were no trends detected in the mutational burden following exposure to NACT in platinum-resistant vs. platinum-sensitive cases. The majority of mutated genes were unique to each case. We identified several genes that were commonly mutated in pre-NACT samples specific to platinum-resistant (CSPG4, SLC35G5, TUBA3D) or sensitive (CYP2D6, NUTM1, DNAH5) cases. Four mutated genes emerged exclusively in the platinum-resistant cases (ADGRV1, MUC17, MUC20, PAK2) following NACT.
Conclusions
Patients with advanced HGSC present with significant intra-tumor heterogeneity. NACT significantly impacts the somatic mutation status irrespective of the time to recurrence. The mutated genes detected in chemo-naive pre-NACT tumor samples from either resistant or sensitive cases could potentially have a role in the prediction of chemotherapy response in patients scheduled to receive NACT; larger studies are required to further validate these genes.
Background
Epithelial ovarian cancer is the most lethal cancer of the female reproductive system in the developed world [1]. As 70% of patients present with advanced disease, the prognosis remains poor with high mortality and a 5-year survival of 48.6% [1]. The most common histological subtype is high-grade serous ovarian cancer (HGSC), which accounts for 90% of stage III/IV cases and 70% of ovarian cancer cases overall [2]. The standard treatment for HGSC includes primary cytoreductive surgery (PCS) and post-operative platinum- and taxane-based chemotherapy [3]. Residual disease after PCS is a major prognostic factor, with the highest survival noted in those patients who undergo a complete cytoreduction to no visible disease [4]. Patients with multiple comorbidities, poor functional status, or where the likelihood of optimal cytoreduction is low, may be considered for neoadjuvant chemotherapy (NACT) followed by interval cytoreductive surgery (ICS) as an alternative option to PCS [3]. Treatment with NACT and ICS has several advantages over PCS including less extensive surgery and reduced perioperative morbidity [5]. The rate of complete cytoreduction is higher in patients undergoing ICS compared to PCS; however, this has not been shown to provide any survival benefit [5]. This has raised concerns regarding the negative impact of NACT and its effect on platinum-resistant recurrence [6]. In addition to increasing the rate of platinum-resistant recurrence, NACT also shortens the time to platinum-resistant recurrence [6]. The response rate to subsequent chemotherapy for all patients with HGSC who develop-platinum resistance is low (15%) with a progression-free survival of 3-4 months and a median survival of less than a year [7].
Although development of platinum resistance following NACT is poorly understood, several possible mechanisms have been suggested. The introduction of NACT to a large volume of tumor pre-operatively may increase the enrichment of the platinum-resistant subclones that are already present at the time of diagnosis [8]. In contrast, this large tumor burden at presentation would be removed in patients undergoing PCS prior to treatment with adjuvant systemic chemotherapy. The limited blood supply of the large tumors may upregulate Hypoxia-inducible factor 1-alpha (HIF-1α), a hypoxia-related transcription factor linked to platinum-resistance [9]. Furthermore, operative surgical planes undergo a treatment effect, making an assessment of optimal excision at ICS challenging [10]. Inadequate tumor resection at ICS combined with NACT-induced stemness of ovarian cancer cells in the residual tumor [11] may result in a faster emergence of platinum resistance, and as a result a decreased survival in comparison to patients treated with PCS. Lastly, NACT has been shown to induce gene mutations such as TP53 and K351N, which are associated with shorter disease-free survival [12].
In this study, we examined the effects of NACT on the somatic mutation status in matched tumor samples obtained pre and post-NACT in women with HGSC to determine if NACT alters/enriches non-synonymous mutations within tumors and to identify candidate mutations that associate with platinum resistance.
Results
A total of 5 cases of stage IIIa-IIIc HGSC treated with 3-5 cycles of platinum-based NACT were selected (Table 1). This included three cases with platinum-resistant (Case 1-3) and two with platinum-sensitive disease (Case 4-5). The biochemical response in terms of reduction of Ca125 level ranged from 66 to 84% for platinum-resistant cases and 96 to 98% for platinum-sensitive cases. The average time to recurrence was 4.2 months in platinum-resistant, and 20.2 months in platinum-sensitive cases. All women with platinum-resistant HGSC died of their disease within 12.1-21.4 months from diagnosis, while both patients with platinum-sensitive HGSC were alive at the last follow-up 42.2 and 47.9 months after the diagnosis.
Whole exome sequencing was performed in buffy coat and macro-dissected pre-NACT (N = 5 samples) and post-NACT tumor (N = 6 samples). A summary of the somatic mutation counts is shown in Table 2 (exonic) and Additional File 1 (non exonic). When comparing exonic non-synonymous mutations in pre-NACT and post-NACT samples from the same patient, an average of 41% (1-68%) of genes were mutated at both time points (Table 2). There were no visible trends in increase vs. decrease in the mutational burden following exposure to NACT in resistant vs. sensitive cases.
The comparison of genes with non-synonymous exonic mutations detected at any time point (pre-NACT, post-NACT, or both) is shown in Fig. 1A. The majority of mutated genes were unique to each case, including 666/730 (91%) in Case 1, 88/112 (79%) in Case 2, 45/66 (68%) in Case 3, 125/151 (83%) in Case 4 and 183/234 (78%) in Case 5. All recurrently mutated genes (> 1 case) are listed in Additional File 2. Of note, all cases had detected mutations in MUC2, while mutations in DDX11, TP53 and TUBA3D were observed in 4 of 5 cases. Genes that were exclusively mutated in > 1 resistant or sensitive cases are highlighted in Table 3; 16 genes were exclusively mutated in > 1 resistant case and 5 genes were exclusively mutated in both sensitive cases.
The comparison of genes with non-synonymous exonic mutations detected within a given time point is shown in Fig. 1B (pre-NACT) and 1C (post-NACT). The majority of mutations were case-specific, including 85-96% among pre-NACT samples and 68-96% among post-NACT samples. However, a few genes were commonly mutated in pre-NACT samples in exclusively resistant (CSPG4, SLC35G5, TUBA3D) or sensitive (CYP2D6, NUTM1, DNAH5) cases. Among post-NACT samples, 10 genes were exclusively mutated in > 1 resistant case (ADGRV1, AOC1, CACNA1S, MTMR11, MUC17, MUC20, OR52N5, PAK2, PCDHB11, TMEM14B); no genes were exclusively mutated in both sensitive cases at this time point. As shown in Table 2 a small proportion of the same mutations were present at both time points (pre- and post-NACT) within the same case; the comparison of these stably mutated genes is shown in Fig. 1D. Of note, only 2 genes were shared between the resistant cases (SLC35G5, TUBA3D), while there were no shared genes between the sensitive cases. To compare the potential impact of NACT in our cases, we next looked at the genes that were mutated in post-NACT, but not pre-NACT samples, from the same case; this included 683/691 (99%) of post-treatment genes from Case 1, 32/65 (49%) from Case 2, 26/53 (49%) from Case 3, 56/98 (57%) from Case 4 and 25/60 (42%) from Case 5. The comparison of these emergent mutations is shown in Fig. 1E; of note, 6 genes were shared between the resistant cases (highlighted in Fig. 1F). While two of these genes showed a different pattern of mutation in sensitive cases (CACNA1S in pre-NACT only and KIR2DL3 at a low level in both pre-NACT and post-NACT), the remaining 4 genes were exclusively mutated in platinum-resistant cases (ADGRV1, MUC17, MUC20, PAK2). There were no genes with treatment-emergent mutations shared between the two sensitive cases.
The inclusion of post-NACT samples from two different tumor sites (right ovary and omentum) for platinum-resistant Case 3 allowed us to compare the impact of both sites and exposure to NACT on mutational status within the same patient (Fig. 2A; see Additional File 3 for more detailed gene list). Of note, 8 genes were commonly mutated among all samples (highlighted in Fig. 2B; CDK12, FOXJ1, HMCN1, MMRN1, MTMR11, PCDHA6, RBM12, REV3L). The same mutation was detected in all samples, with no significant changes in allele fraction following NACT. Both omental samples had detected mutations in 6 genes irrespective of treatment status, while 2 genes had the same emergent mutations in both post-NACT samples irrespective of tissue site (OR10G9, ZNF28). There were also site-specific differences in post-treatment samples, including 10 genes with emergent mutations only at the ovarian site and 14 genes with emergent mutations only at the omental site (highlighted in Fig. 2C). While five of the genes with emergent mutations showed a different pattern of mutations in sensitive cases (FRG1, KIR2DL3, MUC2, SETD8, ZNF28), the majority were exclusive to platinum-resistant cases (ARHGAP5, ARMC4, CFAP47, CLCNKA, GOLGA6L2, GPR101, KRTAP4-11, MUC17, MUC20, NUP50, OR10G9, PAK2, PCDHB11, PCDHGB6, PGAM1, PLEC, SGSM2, SPATA31D1, TBC1D3B/F, TRIM49, UGT2B11).
The results of targeted deep sequencing of 75 candidate genes selected from our exome data are summarized in Fig. 3 (see Additional File 4 for full panel gene list and selection criteria). Mutations were validated in several genes in platinum-resistant cases, while mutation of DNAH5 was confirmed in both pre-NACT and post-NACT samples from Case 4, the one platinum-sensitive case with sufficient material for targeted sequencing (Fig. 3A). Of note, emergent mutations in post-treatment samples were confirmed in seven genes in platinum-resistant Case 2 (ANKRD12, HSD17B4, KIAA1217, NBEA, SH3RF2, SH3TC2, VTI1A) and one gene in platinum-resistant Case 3 (SPATA31D1). Copy number changes were also observed in several genes prior to and/or following NACT (Fig. 3B). Of note, copy numbers of ≥7 were observed in 2 genes in Case 1 (TMED8, PTCD3), six genes in Case 2 (WSCD1, OR52N5, PRB1, KANK4, KLHL38, ATP1A2), and six genes in Case 3 (MMRN1, CDAN1, HSPB7, LRP1, MTMR11, RBM12). Furthermore, copy number loss was observed in SPATA31D1 following NACT in Case 4 only.
Discussion
Studies have demonstrated that HGSC is a multi-clonal disease with a high degree of genomic instability and extensive intra-tumor heterogeneity in functional cellularity [13], copy number alterations [14], and somatic mutation status [15]. In line with these studies, the majority of the mutated genes in our study were unique to each case. In addition, the large heterogeneity of mutation profiles was also noted within the same case between different tumour sites. Similar findings were demonstrated by Bashashati and colleagues who used whole exome, targeted deep sequencing, and copy number analysis to reveal extensive diversity in 31 spatially separated tumor samples obtained prior to treatment from 6 patients with HGSC [15]. Distinct clones were observed within and between tumor masses from the same patient, with an average of 52% (range 10-91%) of somatic mutations present in every sample [15]. Copy number and mutation profiles evolved independently. In the same study, the clonal evolution analysis in one matched primary/recurrent tumor pair revealed a much higher degree of genetic similarity than expected. Anecdotally, this case had a longer overall survival (74 months), which could suggest that clonal stability (especially following chemotherapy) results in an improved outcome [15]. We identified only 2 mutated genes which were shared between the resistant cases (SLC35G5, TUBA3D) pre-and-post-NACT, while there were no shared mutated genes between the sensitive cases. A study by Schwarz et al. used a whole-genome copy number analysis in 135 spatially and temporally separated samples from 14 patients predominantly treated by NACT [14]. Those with a high degree of clonal expansion had worse OS compared to patients with a low degree of clonal expansion, suggesting a relationship between heterogeneity and outcome. Comparison of pre-and-post-NACT samples from 10 patients revealed small alterations in heterogeneity at the CNV level. The extent of temporal heterogeneity was not predictive of survival; however, pre- and post-NACT samples were not site-matched [14]. Extensive intra- and inter- tumor heterogeneity highlight the importance of tumor sampling from many locations before making any decision about the treatment if we are to rely on the genetic information.
We found a significant impact of NACT on the somatic mutation status in all five patients, irrespective of the time to recurrence, based on the low percentage of shared mutations between pre- and post-NACT samples in all cases. However, there were no visible trends in mutational burden following the exposure to NACT in resistant vs. sensitive cases. Similar findings were demonstrated by Cooke et al. at the CNV level [8]. The authors used array comparative genomic hybridization to analyze six paired pre- and post-neoadjuvant treatment HGSC samples from the CTCR-OV01 clinical study. The results did not show extensive copy number differences, suggesting that platinum resistance develops from pre-existing minor clones.
In this study, we have identified a few commonly mutated genes in pre-NACT tumor samples from either sensitive (DNAH5, CYP2D6, NUTM1) or resistant (CSPG4, TUBA3D, SLC35G5) cases, which suggest that these genes could have a potential utility for prediction of the response to NACT. However, this would have to be validated in a larger cohort of patients.
Dynein Axonemal Heavy Chain 5 (DNAH5) encodes dynein axonemal heavy chain, which is involved in ciliary assembly and cell motility and was initially found to be frequently mutated in patients with primary ciliary dyskinesia [16]. Increased incidence of non-synonymous single-nucleotide mutations and insertions/deletions was found in DNAH2, DNAH5, and DNAH10 genes of CpG-island methylator phenotype positive clear cell renal cell carcinomas [17]. It has been also suggested that the DNAH5 gene can play a role in the development of colorectal cancer [18]. Zhu et al., using the genomic and clinical data from the TCGA project, examined the impact of the somatic mutation in DNAH genes on chemotherapy response in gastric adenocarcinoma patients treated with fluoropyrimidine and platinum [19]. The mutation rates of 13 members of the DNAH gene family were higher in chemo-sensitive than in chemo-resistant patients (mutation rate in DNAH5 gene specifically was 22.8% vs 12.5% respectively). In addition, the patients with DNAH mutations had significantly better overall survival, chemotherapy-free survival, fluoropyrimidine-free survival, and platinum-free survival. The results from TCGA analysis were validated by the authors through targeted sequencing of samples from the authors’ cohort. This confirmed association of DNAH mutations with an improved response to chemotherapy. Since the mutation in the DNAH gene affects microtubule structure and tumor cell movement, the authors speculated that it could act as a protective factor and marker of chemotherapy response. Mutation of DNAH5 was validated by targeted sequencing, including in both pre-NACT and post-NACT samples from platinum-sensitive Case 4, suggesting it is a high confidence finding.
The role of the two other mutated genes in pre-NACT tumors of sensitive cases (CYP2D6 and NUTM1) in chemotherapy response is less clear. The Cytochrome P450 Family 2 Subfamily D Member 6 (CYP2D6) gene encodes a member of the cytochrome P450 superfamily. The CYP2D6 enzyme is involved in the metabolism of 20-25% of clinically used drugs and is the main enzyme responsible for the conversion of tamoxifen into endoxifen, its most important metabolite [20]. Genetic variations in the CYP2D6 gene can cause reduced enzyme activity and influence tamoxifen metabolism. However, the association between CYP2D6 genotype and clinical outcome in patients with breast cancer treated with tamoxifen is controversial [20]. It is unclear how the mutation in CYP2D6 could impact clinical outcome after NACT in patients with ovarian cancer, as carboplatin is excreted primarily by kidneys and paclitaxel is metabolized mainly by cytochromes P450 2C8 and 3A4 [21]. Similarly, the role of NUTM1 in patients with ovarian cancer is unclear. The NUT Midline Carcinoma Family Member 1 (NUTM1) is typically expressed in normal testis and is crucial for male fertility [22]. Rearrangement in NUTM1 has been associated with poorly differentiated carcinomas with variable squamous differentiation originating in midline organs as well as in various other neoplasms including sarcoma-like tumors, poromas, and acute lymphoblastic leukemias [22].
As mentioned, three genes were mutated exclusively in pre-NACT tumor samples of resistant cases (CSPG4, TUBA3D, SLC35G5). The CSPG4 (Chondroitin sulfate proteoglycan 4) gene encodes a transmembrane protein and has been implicated in cell differentiation and migration, angiogenesis and vascularization, glial and oligodendrocyte formation and neuronal network regulation [23]. Overexpression of CSPG4 has been found in several cancer types including triple-negative breast cancer [24], sarcomas [25], and squamous cell carcinoma of the head and neck [24]. Expression of CSPG4 predicted poor survival and resistance to ionizing radiation in glioblastoma [26]. Somatic mutation in the CSPG4 gene was significantly associated with PD-L1 positivity (predictor of the response to immunotherapy with PD-1/PD-L1 inhibitors) in renal cell carcinoma tumor cells in a study by Wang et al. [27]. However, the relationship of somatic mutation in the CSPG4 gene to chemoresistance is yet to be determined. It is possible that the cells with this mutation may have been chemosensitive, since the mutation was not reflected in the surviving tumour cells.
The remaining mutated genes present in pre-NACT tumor samples from resistant cases, TUBA3D and SLC35G5, were the only two genes retained in the post-NACT tumor samples of resistant cases. The Tubulin α 3D (TUBA3D) gene encodes a member of the α tubulin family that heterodimerizes with β tubulin to form the main structural component of the microtubules. These are responsible for cellular structure, cell motility, transport and mitosis [28]. The expression of high levels of β3-tubulin was found to be associated with taxane resistance in ovarian cancer [29]. In contrast high expression of β3-tubulin was found to be associated with sensitivity, rather than resistance, to taxane-based chemotherapy in clear cell ovarian carcinoma by Aoki et al. [30]. However, there is no association of mutation in the TUBA3D gene with chemoresistance yet reported in the literature. The Solute carrier family 35 member G5 (SLC35G5) gene belongs to a solute carrier group of transporters, which transport organic or inorganic molecules across cell or organelle membranes [31]. However, its role in carcinogenesis has not been previously described.
We have identified candidate genes involved in the rapid development of platinum resistance in HGSC patients treated with NACT. Mutations in four genes emerged exclusively in the platinum-resistant cases following NACT (ADGRV1, MUC17, MUC20, PAK2). The Adhesion G Protein-Coupled Receptor V1 (ADGRV1) gene encodes an adhesion G-protein coupled receptor V1; its role in tumorigenesis is unclear. Auguste et al., [32] reported alterations in the ADGRV1 gene in two of three post-chemotherapy samples but not in the chemotherapy-naïve patients with small cell carcinoma of the ovary, hypercalcemic type. Mucin 17 (MUC17) and Mucin 20 (MUC 20) genes encode high-molecular-weight membrane glycoproteins [33, 34]. Overexpression of MUC20 was found to be a predictor of poor outcome and recurrence in colorectal cancer [35]. Low expression of MUC20 was significantly correlated with tumour regression grade in patients with esophageal squamous cell carcinoma treated with neoadjuvant cisplatin and paclitaxel chemotherapy [36]. In contrast, the blockage of MUC20 in in-vitro experiment contributed to increase sensitivity to paclitaxel [36]. Using HGSC genomics data from TCGA, Nabavi et al. identified somatic MUC17 mutation at a similar frequency of platinum-sensitive and platinum-resistant groups [37]. The P21-activated protein kinase 2 (PAK 2) is a serine/threonine kinase, which regulates cell motility/migration, and gene transcription [38]. It has been found to promote tumor cell proliferation and survival through the AKT1 and Raf–MAPK pathways [39]. Overexpression of PAK1 and PAK2 has been found in ovarian cancer cells [40]. While Gupta et al. found that elevated expression of PAK2 leads to cell proliferation and acquired chemoresistance in head and neck cancer through activation of c-Myc and PKM2 [41]. In comparison, Shuang et al., found that repression of miR-134 and leading to subsequent up-regulation of PAK2 expression contribute to paclitaxel resistance in HGSC [42]. Authors concluded that both, the miR-134 and its target PAK2 might be a potential target for therapeutic intervention of ovarian cancer paclitaxel resistance.
The additional TCGA-OV screening [43] for the commonly mutated genes found in our study, revealed only a small number of individuals with the mutations in these genes (Table 4). While the larger numbers could have been more supportive of the results of our study, it is worth noticing that the overall sample size was also limited (n = 436) and much bigger sample size would be needed for positive results.
The main limitation of our study is the small number of patients included. This is largely due to the challenge in obtaining sufficient pre and post-NACT sample from retrospectively collected cases which is largely limited by the lack of pre-NACT biopsy material. In addition, we acknowledge the lack of CRISPR-cas9 whole-genome screening in our study. While it is a powerful gene-editing tool that has been previously utilized to identify genes involved in platinum-resistance in HGSC (including but not limited to TP53, p53, ABCB1, AKT1, ERCC1, EGFR, BRCA1, PIK3CA, and MAPK1) [44, 45], we did not have access to this technology at the time of study conception. We believe that due to the small sample size, it could have been of a limited additional value. Despite the small number of cases, we were able to detect recurrently mutated genes in cases resistant to NACT (detected prior to treatment and/or following exposure to NACT). Although the impact of the identified mutations on function of the encoded protein is not known. The mutation could reflect a loss or gain of function or have no impact at all. In addition, mutations in non-exon regions could also have impact, particularly if these lie in regulatory (enhancer) regions. Study of a larger number of cases, preferably in a prospective setting, may reveal additional genes related to NACT response. Future studies should focus on further exploration of site-specific post-NACT changes in comparison to malignant ascites and/or plasma samples in platinum-resistant recurrence. Direct experimentation of the role of candidate genes and mutations in platinum resistance using patient-derived ascites may yield additional valuable information.
Conclusion
In summary, our study identified genes mutated in patients with HGSC pre-NACT that are associated with platinum resistance. Additional studies are necessary to further validate these findings and identify additional genes associated with response to chemotherapy.
Methods
Study cases and samples
HGSC patients with a documented response to platinum-based NACT who had donated a baseline blood sample to the Princess Margaret Gynecologic Oncology BioBank and had banked tumor before and after NACT were identified. Ethics approval was obtained from the University Health Network Research Ethics Board prior to case review (protocol #14-7539). Available tumor tissue blocks from all cases were subjected to secondary pathological review. Pre-NACT (at diagnosis, prior to NACT, “Pre”) and post-NACT (at surgery after 3-5 cycles of NACT, “Post”) tumor sites were matched when possible. A total of 12 tumor blocks were selected for whole exome and/or targeted sequencing, including: a pre-NACT omental biopsy and post-NACT samples from the omentum and pelvic nodule for Case 1; a pre-NACT omental biopsy and post-NACT omental tumor for Case 2; a pre-NACT omental biopsy and post-NACT samples from the right ovary and omentum for Case 3; a pre-NACT omental biopsy and a post-NACT sample from a stomach nodule for Case 4; a pre-NACT biopsy (location not specified) and a post-NACT sample from a bowel nodule for Case 5. Tumor content was enriched to > 80% cellularity through H&E-guided macrodissection of unstained sections or tissue coring prior to DNA extraction using the Qiagen QIAamp DNA FFPE extraction kit (Qiagen, Toronto, Ontario, Canada). DNA was eluted in Buffer ATE, quantified using the Qubit Fluorometer (Life Technologies, California, USA) and quality assessed (BioAnalyzer, TapeStation, qPCR). Matched buffy coat samples were included for each case as a germline control. Platinum-resistant recurrence was defined as disease recurrence within 6 months of first-line platinum-based chemotherapy, while platinum-sensitive recurrence was defined as disease recurrence 6 months and more after completion of first-line platinum-based chemotherapy [7].
Whole exome sequencing
Genomic DNA (200 ng) from the buffy coat, pre-NACT tumor, and post-NACT tumor samples was fragmented using a Covaris Focused-ultrasonicator, libraries were generated using the Agilent SureSelect XT Human All Exon v5 + UTR kit, and sequencing was performed on the Illumina HiSeq2000 using a 100-cycle paired-end protocol to achieve 50X (normal samples) or 100X (tumor samples) coverage. Genome Analysis Toolkit2 (GATK) Best Practices recommendations for sequencing analysis were followed. Briefly, fastqs were aligned to the reference human genome (hg19) using BWA-MEM (v0.7.12). BAM files were processed using several tools (Picard MarkDuplicates, GATK RealignerTargetCreator followed by IndelRealigner, GATK’s Base Quality Score Recalibration) and quality metrics were generated using Qualimap (v2.1). Candidate somatic mutations were identified using MuTect (v1.1.5), followed by filtering of VCF files (using VCFtools) to retain only high confidence somatic variants and annotation using ANNOVAR (v20160201).
Deep targeted sequencing using a custom panel
A custom SureSelect XT panel was designed using the SureDesign Tool by Agilent Technologies (Mississauga, Ontario, Canada), with 1776 included regions, a total target region size of 325.797 kbp, and a total probe size of 458.215 kbp. All exons +/− 10 bp for each gene were included, with > 99% coverage of the target region. Seventy-five genes were included, based on at least one of the following criteria: 1) presence of a mutation in only platinum-resistant (or sensitive) cases; 2) detection of a mutation in > 1 resistant case; 3) detection of a mutation only in, or at an increased percentage of variant reads in, post-NACT tumor samples from platinum-resistant cases; 4) previous implication in platinum response in ovarian or other malignancies; 5) association with PFS within The Cancer Genome Atlas (TCGA) ovarian serous carcinoma dataset.
Targeted capture and sequencing on Illumina HiSeq2000 was performed using 200 ng DNA from 9 tumor and 4 matched normal samples from Cases 1-4 (See Table 1). Samples from Case 5 were not included due to insufficient remaining DNA in the pre-NACT tumor biopsy. 5000X coverage for tumor, and 100X for normal, was achieved using the 100-cycle paired-end protocol and multiplexing. Alignment of fastqs and processing of BAM files was conducted as described for whole-exome sequencing above. MuTect (v1.1.5) was used to identify candidate somatic mutations; for targeted sequencing samples, a bed file with intervals was provided. The VCF output from this step was then filtered (using VCFtools) to retain only high confidence somatic variants. Variants were annotated using ANNOVAR (v20160201). SNVs and indels were detected using Samtools mpileup (v1.2) and VarScan (v2.3.8) and copy number variants (CNVs) were identified using VarScan2 and further verified with Sequenza (v2.1.2).
Availability of data and materials
The dataset used and/or analyzed during the current study are available from the corresponding authors on request.
References
SEER Stat Fact Sheets: Ovarian Cancer. [March 22 2021]; Available from: http://seer.cancer.gov/statfacts/html/ovary.html.
Cannistra SA. Cancer of the ovary. N Engl J Med. 2004;351(24):2519–29.
https://www.nccn.org/professionals/physician_gls/pdf/ovarian.pdf (Accessed 22 Mar 2021).
May T, Altman A, McGee J, et al. Examining Survival Outcomes of 852 Women With Advanced Ovarian Cancer: A Multi-institutional Cohort Study. Int J Gynecol Cancer. 2018;28(5):925–31. https://doi.org/10.1097/IGC.0000000000001244.
Kim SR, Kotsopoulos J, Sun P, et al. The impacts of neoadjuvant chemotherapy and of cytoreductive surgery on 10-year survival from advanced ovarian cancer [published online ahead of print, 2020 Dec 16]. Int J Gynaecol Obstet. 2020. https://doi.org/10.1002/ijgo.13542.
da Costa AA, Valadares CV, Baiocchi G, et al. Neoadjuvant Chemotherapy Followed by Interval Debulking Surgery and the Risk of Platinum Resistance in Epithelial Ovarian Cancer. Ann Surg Oncol. 2015;22(Suppl 3):S971–8. https://doi.org/10.1245/s10434-015-4623-z.
Davis A, Tinker AV, Friedlander M. "Platinum resistant" ovarian cancer: what is it, who to treat and how to measure benefit? Gynecol Oncol. 2014;133(3):624–31. https://doi.org/10.1016/j.ygyno.2014.02.038.
Cooke SL, Ng CK, Melnyk N, et al. Genomic analysis of genetic heterogeneity and evolution in high-grade serous ovarian carcinoma. Oncogene. 2010;29(35):4905–13. https://doi.org/10.1038/onc.2010.245.
Ai Z, Lu Y, Qiu S, Fan Z. Overcoming cisplatin resistance of ovarian cancer cells by targeting HIF-1-regulated cancer metabolism. Cancer Lett. 2016;373(1):36–44. https://doi.org/10.1016/j.canlet.2016.01.009.
Hynninen J, Lavonius M, Oksa S, Grénman S, Carpén O, Auranen A. Is perioperative visual estimation of intra-abdominal tumor spread reliable in ovarian cancer surgery after neoadjuvant chemotherapy? Gynecol Oncol. 2013;128(2):229–32. https://doi.org/10.1016/j.ygyno.2012.11.007.
Lim MC, Song YJ, Seo SS, Yoo CW, Kang S, Park SY. Residual cancer stem cells after interval cytoreductive surgery following neoadjuvant chemotherapy could result in poor treatment outcomes for ovarian cancer. Onkologie. 2010;33(6):324–30. https://doi.org/10.1159/000313823.
Zhang GN, Liu H, Huang JM, et al. TP53 K351N mutation-associated platinum resistance after neoadjuvant chemotherapy in patients with advanced ovarian cancer. Gynecol Oncol. 2014;132(3):752–7. https://doi.org/10.1016/j.ygyno.2014.01.028.
Stewart JM, Shaw PA, Gedye C, Bernardini MQ, Neel BG, Ailles LE. Phenotypic heterogeneity and instability of human ovarian tumor-initiating cells. Proc Natl Acad Sci U S A. 2011;108(16):6468–73.
Schwarz RF, Ng CK, Cooke SL, Newman S, Temple J, Piskorz AM, et al. Spatial and temporal heterogeneity in high-grade serous ovarian cancer: a phylogenetic analysis. PLoS Med. 2015;12(2):e1001789.
Bashashati A, Ha G, Tone A, Ding J, Prentice LM, Roth A, et al. Distinct evolutionary trajectories of primary high-grade serous ovarian cancers revealed through spatial mutational profiling. J Pathol. 2013;231(1):21–34.
Hornef N, Olbrich H, Horvath J, et al. DNAH5 mutations are a common cause of primary ciliary dyskinesia with outer dynein arm defects. Am J Respir Crit Care Med. 2006;174(2):120–6. https://doi.org/10.1164/rccm.200601-084OC.
Arai E, Gotoh M, Tian Y, et al. Alterations of the spindle checkpoint pathway in clinicopathologically aggressive CpG island methylator phenotype clear cell renal cell carcinomas. Int J Cancer. 2015;137(11):2589–606. https://doi.org/10.1002/ijc.29630.
Xiao WH, Qu XL, Li XM, et al. Identification of commonly dysregulated genes in colorectal cancer by integrating analysis of RNA-Seq data and qRT-PCR validation. Cancer Gene Ther. 2015;22(5):278–84. https://doi.org/10.1038/cgt.2015.20.
Zhu C, Yang Q, Xu J, et al. Somatic mutation of DNAH genes implicated higher chemotherapy response rate in gastric adenocarcinoma patients. J Transl Med. 2019;17(1):109. Published 2019 Apr 3. https://doi.org/10.1186/s12967-019-1867-6.
Mulder TAM, de With M, Del Re M, Danesi R, Mathijssen RHJ, van Schaik RHN. Clinical CYP2D6 Genotyping to Personalize Adjuvant Tamoxifen Treatment in ER-Positive Breast Cancer Patients: Current Status of a Controversy. Cancers (Basel). 2021;13(4):771. Published 2021 Feb 12. https://doi.org/10.3390/cancers13040771.
Spratlin J, Sawyer MB. Pharmacogenetics of paclitaxel metabolism. Crit Rev Oncol Hematol. 2007;61(3):222–9. https://doi.org/10.1016/j.critrevonc.2006.09.006.
McEvoy CR, Fox SB, Prall OWJ. Emerging entities in NUTM1-rearranged neoplasms. Genes Chromosom Cancer. 2020;59(6):375–85. https://doi.org/10.1002/gcc.22838.
Ilieva KM, Cheung A, Mele S, et al. Chondroitin Sulfate Proteoglycan 4 and Its Potential As an Antibody Immunotherapy Target across Different Tumor Types. Front Immunol. 2018;8:1911. Published 2018 Jan 10. https://doi.org/10.3389/fimmu.2017.01911.
Wang X, Osada T, Wang Y, et al. CSPG4 protein as a new target for the antibody-based immunotherapy of triple-negative breast cancer. J Natl Cancer Inst. 2010;102(19):1496–512. https://doi.org/10.1093/jnci/djq343.
Campoli M, Ferrone S, Wang X. Functional and clinical relevance of chondroitin sulfate proteoglycan 4. Adv Cancer Res. 2010;109:73–121. https://doi.org/10.1016/B978-0-12-380890-5.00003-X.
Svendsen A, Verhoeff JJ, Immervoll H, et al. Expression of the progenitor marker NG2/CSPG4 predicts poor survival and resistance to ionising radiation in glioblastoma. Acta Neuropathol. 2011;122(4):495–510. https://doi.org/10.1007/s00401-011-0867-2.
Wang J, Xi Z, Xi J, et al. Somatic mutations in renal cell carcinomas from Chinese patients revealed by whole exome sequencing. Cancer Cell Int. 2018;18:159. Published 2018 Oct 17. https://doi.org/10.1186/s12935-018-0661-5.
Berrieman HK, Lind MJ, Cawkwell L. Do beta-tubulin mutations have a role in resistance to chemotherapy? Lancet Oncol. 2004;5(3):158–64. https://doi.org/10.1016/S1470-2045(04)01411-1.
Kavallaris M, Kuo DY, Burkhart CA, et al. Taxol-resistant epithelial ovarian tumors are associated with altered expression of specific beta-tubulin isotypes. J Clin Invest. 1997;100(5):1282–93. https://doi.org/10.1172/JCI119642.
Aoki D, Oda Y, Hattori S, et al. Overexpression of class III beta-tubulin predicts good response to taxane-based chemotherapy in ovarian clear cell adenocarcinoma. Clin Cancer Res. 2009;15(4):1473–80. https://doi.org/10.1158/1078-0432.CCR-08-1274.
Song Z. Roles of the nucleotide sugar transporters (SLC35 family) in health and disease. Mol Asp Med. 2013;34(2-3):590–600. https://doi.org/10.1016/j.mam.2012.12.004.
Auguste A, Blanc-Durand F, Deloger M, et al. Small Cell Carcinoma of the Ovary, Hypercalcemic Type (SCCOHT) beyond SMARCA4 Mutations: A Comprehensive Genomic Analysis. Cells. 2020;9(6):1496. Published 2020 Jun 19. https://doi.org/10.3390/cells9061496.
Gum JR Jr, Crawley SC, Hicks JW, Szymkowski DE, Kim YS. MUC17, a novel membrane-tethered mucin. Biochem Biophys Res Commun. 2002;291(3):466–75. https://doi.org/10.1006/bbrc.2002.6475.
Higuchi T, Orita T, Nakanishi S, et al. Molecular cloning, genomic structure, and expression analysis of MUC20, a novel mucin protein, up-regulated in injured kidney. J Biol Chem. 2004;279(3):1968–79. https://doi.org/10.1074/jbc.M304558200.
Xiao X, Wang L, Wei P, et al. Role of MUC20 overexpression as a predictor of recurrence and poor outcome in colorectal cancer. J Transl Med. 2013;11:151. Published 2013 Jun 20. https://doi.org/10.1186/1479-5876-11-151.
Shen LY, Wang H, Dong B, et al. Possible prediction of the response of esophageal squamous cell carcinoma to neoadjuvant chemotherapy based on gene expression profiling. Oncotarget. 2016;7(4):4531–41. https://doi.org/10.18632/oncotarget.6554.
Nabavi S. Identifying candidate drivers of drug response in heterogeneous cancer by mining high throughput genomics data. BMC Genomics. 2016;17(1):638. Published 2016 Aug 15. https://doi.org/10.1186/s12864-016-2942-5.
Bokoch GM. Biology of the p21-activated kinases. Annu Rev Biochem. 2003;72:743–81. https://doi.org/10.1146/annurev.biochem.72.121801.161742.
Menges CW, Sementino E, Talarchek J, et al. Group I p21-activated kinases (PAKs) promote tumor cell proliferation and survival through the AKT1 and Raf-MAPK pathways. Mol Cancer Res. 2012;10(9):1178–88. https://doi.org/10.1158/1541-7786.MCR-12-0082.
Siu MK, Wong ES, Chan HY, et al. Differential expression and phosphorylation of Pak1 and Pak2 in ovarian cancer: effects on prognosis and cell invasion. Int J Cancer. 2010;127(1):21–31. https://doi.org/10.1002/ijc.25005.
Gupta A, Ajith A, Singh S, Panday RK, Samaiya A, Shukla S. PAK2-c-Myc-PKM2 axis plays an essential role in head and neck oncogenesis via regulating Warburg effect. Cell Death Dis. 2018;9(8):825. Published. https://doi.org/10.1038/s41419-018-0887-0.
Shuang T, Wang M, Shi C, Zhou Y, Wang D. Down-regulated expression of miR-134 contributes to paclitaxel resistance in human ovarian cancer cells. FEBS Lett. 2015;589(20 Pt B):3154–64. https://doi.org/10.1016/j.febslet.2015.08.047.
Li H, Li J, Gao W, Zhen C, Feng L. Systematic analysis of ovarian cancer platinum-resistance mechanisms via text mining. J Ovarian Res. 2020;13(1):27. Published 2020 Mar 11. https://doi.org/10.1186/s13048-020-00627-6.
Walton JB, Farquharson M, Mason S, et al. CRISPR/Cas9-derived models of ovarian high grade serous carcinoma targeting Brca1, Pten and Nf1, and correlation with platinum sensitivity published correction appears in Sci Rep. 2018;13;8(1):6231. Sci Rep 2017;7(1):16827 Published 2017 Dec 4. https://doi.org/10.1038/s41598-017-17119-1.
Acknowledgements
Julissa Tsao for help with DNA extraction and sequencing; Neke Ibeh for help with data analysis.
Princess Margaret Gynecologic Oncology BioBank.
Funding
Phil and Terry Browman Grant Fund (Princess Margaret Cancer Foundation).
Department of Obstetrics and Gynaecology, Mount Sinai Hospital and the University Health Network, University of Toronto.
Author information
Authors and Affiliations
Contributions
Zibi Marchocki – interpretation of results, manuscript preparation, clinical expertise. Alicia Tone – study design, case/sample selection, macrodissection, interpretation of results, manuscript preparation. Jocelyn Stewart – study design, case/sample selection, manuscript review. Blaise Clarke – review of pathology/selection of representative areas for sequencing, manuscript review. Julissa Tsao – DNA extraction, sequencing experiments, manuscript review. Carl Virtanen, Neke Ibeh, Richard de Borja – bioinformatics analysis, manuscript review. Theodore Brown – scientific expertise, manuscript review. Taymaa May – study design, interpretation of results, manuscript preparation, clinical expertise. The authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Ethics approval was obtained from the University Health Network Research Ethics Board prior to case review (protocol #14-7539).
Consent for publication
Written consent was obtained from each patient prior to collecting the blood and tissues samples.
Competing interests
None of the authors have a competing interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Marchocki, Z., Tone, A., Virtanen, C. et al. Impact of neoadjuvant chemotherapy on somatic mutation status in high-grade serous ovarian carcinoma. J Ovarian Res 15, 50 (2022). https://doi.org/10.1186/s13048-022-00983-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13048-022-00983-5