Search

Abstract Next generation sequencing of cell-free DNA (cfDNA) is a promising method for treatment monitoring and therapy selection in metastatic breast cancer (MBC). However, distinguishing tumor-specific variants from sequencing artefacts and germline variation with low false discovery rate is challenging when using large targeted sequencing panels covering many tumor suppressor genes. To address this, we built a machine learning model to remove false positive variant calls and augmented it with additional filters to ensure selection of tumor-derived variants. We used cfDNA of 70 MBC patients profiled with both the small targeted Oncomine breast panel (Thermofisher) and the much larger Qiaseq Human Breast Cancer Panel (Qiagen). The model was trained on the panels’ common regions using Oncomine hotspot mutations as ground truth. Applied to Qiaseq data, it achieved 35% sensitivity and 36% precision, outperforming basic filtering. For 20 patients we used germline DNA to filter for somatic variants and obtained 245 variants in total, while our model found seven variants, of which six were also detected using the germline strategy. In ten tumor-free individuals, our method detected in total one (potentially germline) variant, in contrast to 521 variants detected without our model. These results indicate that our model largely detects somatic variants.

Abstract Background In breast cancer (BC), recurrent fusion genes of estrogen receptor alpha (ESR1) and AKAP12, ARMT1 and CCDC170 have been reported. In these gene fusions the ligand binding domain of ESR1 has been replaced by the transactivation domain of the fusion partner constitutively activating the receptor. As a result, these gene fusions can drive tumor growth hormone independently as been shown in preclinical models, but the clinical value of these fusions have not been reported. Here, we studied the prognostic and predictive value of different frequently reported ESR1 fusion transcripts in primary BC. Methods We evaluated 732 patients with primary BC (131 ESR1-negative and 601 ESR1-positive cases), including two ER-positive BC patient cohorts: one cohort of 322 patients with advanced disease who received first-line endocrine therapy (ET) (predictive cohort), and a second cohort of 279 patients with lymph node negative disease (LNN) who received no adjuvant systemic treatment (prognostic cohort). Fusion gene transcript levels were measured by reverse transcriptase quantitative PCR. The presence of the different fusion transcripts was associated, in uni- and multivariable Cox regression analysis taking along current clinico-pathological characteristics, to progression free survival (PFS) during first-line endocrine therapy in the predictive cohort, and disease- free survival (DFS) and overall survival (OS) in the prognostic cohort. Results The ESR1-CCDC170 fusion transcript was present in 27.6% of the ESR1-positive BC subjects and in 2.3% of the ESR1-negative cases. In the predictive cohort, none of the fusion transcripts were associated with response to first-line ET. In the prognostic cohort, the median DFS and OS were respectively 37 and 93 months for patients with an ESR1-CCDC170 exon 8 gene fusion transcript and respectively 91 and 212 months for patients without this fusion transcript. In a multivariable analysis, this ESR1-CCDC170 fusion transcript was an independent prognostic factor for DFS (HR) (95% confidence interval (CI): 1.8 (1.2–2.8), P = 0.005) and OS (HR (95% CI: 1.7 (1.1–2.7), P = 0.023). Conclusions Our study shows that in primary BC only ESR1-CCDC170 exon 8 gene fusion transcript carries prognostic value. None of the ESR1 fusion transcripts, which are considered to have constitutive ER activity, was predictive for outcome in BC with advanced disease treated with endocrine treatment.

Abstract Here, we describe a novel approach for rapid discovery of a set of tumor-specific genomic structural variants (SVs), based on a combination of low coverage cancer genome sequencing using Oxford Nanopore with an SV calling and filtering pipeline. We applied the method to tumor samples of high-grade ovarian and prostate cancer patients and validated on average ten somatic SVs per patient with breakpoint-spanning PCR mini-amplicons. These SVs could be quantified in ctDNA samples of patients with metastatic prostate cancer using a digital PCR assay. The results suggest that SV dynamics correlate with and may improve existing treatment-response biomarkers such as PSA. https://github.com/UMCUGenetics/SHARC .

Abstract Background Extracellular vesicles (EVs) are actively secreted by cells into body fluids and contain nucleic acids of the cells they originate from. The goal of this study was to detect circulating tumor-derived EVs (ctEVs) by mutant mRNA transcripts (EV-RNA) in plasma of patients with solid cancers and compare the occurrence of ctEVs with circulating tumor DNA (ctDNA) in cell-free DNA (cfDNA). Methods For this purpose, blood from 20 patients and 15 healthy blood donors (HBDs) was collected in different preservation tubes (EDTA, BCT, CellSave) and processed into plasma within 24 h from venipuncture. EVs were isolated with the ExoEasy protocol from this plasma and from conditioned medium of 6 cancer cell lines and characterized according to MISEV2018-guidelines. RNA from EVs was isolated with the ExoRNeasy protocol and evaluated for transcript expression levels of 96 genes by RT-qPCR and genotyped by digital PCR. Results Our workflow applied on cell lines revealed a high concordance between cellular mRNA and EV-RNA in expression levels as well as variant allele frequencies for PIK3CA, KRAS and BRAF. Plasma CD9-positive EV and GAPDH EV-RNA levels were significantly different between the preservation tubes. The workflow detected only ctEVs with mutant transcripts in plasma of patients with high amounts (> 20%) of circulating tumor DNA (ctDNA). Expression profiling showed that the EVs from patients resemble healthy donors more than tumor cell lines supporting that most EVs are derived from healthy tissue. Conclusions We provide a workflow for ctEV detection by spin column-based generic isolation of EVs and PCR-based measurement of gene expression and mutant transcripts in EV-RNA derived from cancer patients’ blood plasma. This workflow, however, detected tumor-specific mutations in blood less often in EV-RNA than in cfDNA.

Quantification of tumor‐specific variants (TSVs) in cell‐free DNA is rapidly evolving as a prognostic and predictive tool in patients with cancer. Currently, both variant allele frequency (VAF) and number of mutant molecules per mL plasma are used as units of measurement to report those TSVs. However, it is unknown to what extent both units of measurement agree and what are the factors underlying an existing disagreement. To study the agreement between VAF and mutant molecules in current clinical studies, we analyzed 1116 TSVs from 338 patients identified with next‐generation sequencing (NGS) or digital droplet PCR (ddPCR). On different study cohorts, a Deming regression analysis was performed and its 95% prediction interval was used as surrogate for the limits of agreement between VAF and number of mutant molecules per mL and to identify outliers. VAF and number of mutant molecules per mL plasma yielded greater agreement when using ddPCR than NGS. In case of discordance between VAF and number of mutant molecules per mL, insufficient molecular coverage in NGS and high cell‐free DNA concentration were the main responsible factors. We propose several optimization steps needed to bring monitoring of TSVs in cell‐free DNA to its full potential.

In metastatic colorectal cancer, RAS and BRAF mutations cause resistance to anti‐EGFR therapies, such as cetuximab. Heterogeneity in RAS and BRAF mutations might explain nonresponse in a subset of patients receiving cetuximab. Analyzing mutations in plasma‐derived circulating tumor DNA (ctDNA) could provide a more comprehensive overview of the mutational landscape as compared to analyses of primary and/or metastatic tumor tissue. Therefore, this prospective multicenter study followed 34 patients with metastatic colorectal cancer who were tissue‐tested as RAS wild‐type (exons 2–4) during routine work‐up and received third‐line cetuximab monotherapy. BRAF mutation status was also tested but did not exclude patients from therapy. At baseline and upon disease progression, cell‐free DNA (cfDNA) was isolated for targeted next‐generation sequencing (NGS). At 8 weeks, we determined that patients had benefited from treatment. NGS of cfDNA identified three patients with RAS mutations not detected in tumor tissue during routine work‐up. Another six patients had a BRAF or rare RAS mutation in ctDNA and/or tumor tissue. Relative to patients without mutations in RAS/BRAF, patients with mutations at baseline had shorter progression‐free survival [1.8 versus 4.9 months (P

The emerging interest in circulating tumor DNA (ctDNA) analyses for clinical trials has necessitated the development of a high‐throughput method for fast, reproducible, and efficient isolation of ctDNA. Currently, the majority of ctDNA studies use the manual QIAamp (QA) platform to isolate DNA from blood. The purpose of this study was to compare two competing automated DNA isolation platforms [Maxwell (MX) and QIAsymphony (QS)] to the current ‘gold standard’ QA to facilitate high‐throughput processing of samples in prospective trials. We obtained blood samples from healthy blood donors and metastatic cancer patients for plasma isolation. Total cell‐free DNA (cfDNA) quantity was assessed by TERT quantitative PCR. Recovery efficiency was investigated by quantitative PCR analysis of spiked‐in synthetic plant DNA. In addition, a β‐actin fragmentation assay was performed to determine the amount of contamination by genomic DNA from lysed leukocytes. ctDNA quality was assessed by digital PCR for somatic variant detection. cfDNA quantity and recovery efficiency were lowest using the MX platform, whereas QA and QS showed a comparable performance. All platforms preferentially isolated small (136 bp) DNA fragments over large (420 and 2000 bp) DNA fragments. Detection of the number variant and wild‐type molecules was most comparable between QA and QS. However, there was no significant difference in variant allele frequency comparing QS and MX to QA. In summary, we show that the QS platform has comparable performance to QA, the ‘gold standard’, and outperformed the MX platform depending on the readout used. We conclude that the QS can replace the more laborious QA platform, especially when high‐throughput cfDNA isolation is needed.

Mutations and splice variants in the estrogen receptor (ER) gene, ESR1, may yield endocrine resistance in metastatic breast cancer (MBC) patients. These putative endocrine resistance markers are likely to emerge during treatment, and therefore, its detection in liquid biopsies, such as circulating tumor cells (CTCs) and cell‐free DNA (cfDNA), is of great interest. This research aimed to determine whether ESR1 mutations and splice variants occur more frequently in CTCs of MBC patients progressing on endocrine treatment. In addition, the presence of ESR1 mutations was evaluated in matched cfDNA and compared to CTCs. CellSearch‐enriched CTC fractions (≥5/7.5 mL) of two MBC cohorts were evaluated, namely (a) patients starting first‐line endocrine therapy (n = 43, baseline cohort) and (b) patients progressing on any line of endocrine therapy (n = 40, progressing cohort). ESR1 hotspot mutations (D538G and Y537S/N/C) were evaluated in CTC‐enriched DNA using digital PCR and compared with matched cfDNA (n = 18 baseline cohort; n = 26 progressing cohort). Expression of ESR1 full‐length and 4 of its splice variants (∆5, ∆7, 36 kDa, and 46 kDa) was evaluated in CTC‐enriched mRNA. It was observed that in the CTCs, the ESR1 mutations were not enriched in the progressing cohort (8%), when compared with the baseline cohort (5%) (P = 0.66). In the cfDNA, however, ESR1 mutations were more prevalent in the progressing cohort (42%) than in the baseline cohort (11%) (P = 0.04). Three of the same mutations were observed in both CTCs and cfDNA, 1 mutation in CTCs only, and 11 in cfDNA only. Only the ∆5 ESR1 splice variant was CTC‐specific expressed, but was not enriched in the progressing cohort. In conclusion, sensitivity for detecting ESR1 mutations in CTC‐enriched fractions was lower than for cfDNA. ESR1 mutations detected in cfDNA, rarely present at the start of first‐line endocrine therapy, were enriched at progression, strongly suggesting a role in conferring endocrine resistance in MBC.

Abstract The use of blood-circulating cell-free DNA (cfDNA) as ‘liquid-biopsy’ is explored worldwide, with hopes for its potential in providing prognostic or predictive information in cancer treatment. In exploring cfDNA, valuable repositories are biobanks containing material collected over time, however these retrospective cohorts have restrictive resources. In this study, we aimed to detect tumor-specific mutations in only minute amounts of serum-derived cfDNA by using a targeted next generation sequencing (NGS) approach. In a retrospective cohort of ten metastatic breast cancer patients, we profiled DNA from primary tumor tissue (frozen), tumor-adjacent normal tissue (formalin-fixed paraffin embedded), and three consecutive serum samples (frozen). Our presented workflow includes comparisons with matched normal DNA or in silico reference DNA to discriminate germline from somatic variants, validation of variants through the detection in at least two DNA samples of an individual, and the use of public databases on variants. By our workflow, we were able to detect a total of four variants traceable as circulating tumor DNA (ctDNA) in the sera of three of the ten patients.

The aim of this study was to determine an optimal workflow to detect TP53 mutations in baseline and longitudinal serum cell free DNA (cfDNA) from high-grade serous ovarian carcinomas (HGSOC) patients and to define whether TP53 mutations are suitable as biomarker for disease. TP53 was investigated in tissue and archived serum from 20 HGSOC patients by a next-generation sequencing (NGS) workflow alone or combined with digital PCR (dPCR). AmpliSeq™-focused NGS panels and customized dPCR assays were used for tissue DNA and longitudinal cfDNAs, and Oncomine NGS panel with molecular barcoding was used for baseline cfDNAs. TP53 missense mutations were observed in 17 tissue specimens and in baseline cfDNA for 4/8 patients by AmpliSeq, 6/9 patients by Oncomine, and 4/6 patients by dPCR. Mutations in cfDNA were detected in 4/6 patients with residual disease and 3/4 patients with disease progression within six months, compared to 5/11 patients with no residual disease and 6/13 patients with progression after six months. Finally, mutations were detected at progression in 5/6 patients, but not during chemotherapy. NGS with molecular barcoding and dPCR were most optimal workflows to detect TP53 mutations in baseline and longitudinal serum cfDNA, respectively. TP53 mutations were undetectable in cfDNA during treatment but re-appeared at disease progression, illustrating its promise as a biomarker for disease monitoring.

Introduction The locoregional failure (LRF) rate in human papilloma virus (HPV)-negative oropharyngeal squamous cell carcinoma (OPSCC) remains disappointingly high and toxicity is substantial. Response prediction prior to or early during treatment would provide opportunities for personalised treatment. Currently, there are no accurate predictive models available for correct OPSCC patient selection. Apparently, the pivotal driving forces that determine how a OPSCC responds to treatment, have yet to be elucidated. Therefore, the holistiC early respOnse assessMent for oroPharyngeaL cancer paTiEnts study focuses on a holistic approach to gain insight in novel potential prognostic biomarkers, acquired before and early during treatment, to predict response to treatment in HPV-negative patients with OPSCC.Methods and analysis This single-centre prospective observational study investigates 60 HPV-negative patients with OPSCC scheduled for primary radiotherapy (RT) with cisplatin or cetuximab, according to current clinical practice. A holistic approach will be used that aims to map the macroscopic (with Intra Voxel Incoherent Motion Diffusion Kurtosis Imaging (IVIM-DKI); before, during, and 3 months after RT), microscopic (with biopsies of the primary tumour acquired before treatment and irradiated ex vivo to assess radiosensitivity), and molecular landscape (with circulating tumour DNA (ctDNA) analysed before, during and 3 months after treatment). The main end point is locoregional control (LRC) 2 years after treatment. The primary objective is to determine whether a relative change in the mean of the diffusion coefficient D (an IVIM-DKI parameter) in the primary tumour early during treatment, improves the performance of a predictive model consisting of tumour volume only, for 2 years LRC after treatment. The secondary objectives investigate the potential of other IVIM-DKI parameters, ex vivo sensitivity characteristics, ctDNA, and combinations thereof as potential novel prognostic markers.Ethics and dissemination The study was approved by the Medical Ethical Committee of Erasmus Medical Center. The main results of the trial will be presented in international meetings and medical journals.Trial registration number NL8458.