Anna Piskorz

Non-invasive analysis of acquired resistance to cancer therapy by sequencing of plasma DNA

Cancers acquire resistance to systemic treatment as a result of clonal evolution and selection. Repeat biopsies to study genomic evolution as a result of therapy are difficult, invasive and may be confounded by intra-tumour heterogeneity. Recent studies have shown that genomic alterations in solid cancers can be characterized by massively parallel sequencing of circulating cell-free tumour DNA released from cancer cells into plasma, representing a non-invasive liquid biopsy. Here we report sequencing of cancer exomes in serial plasma samples to track genomic evolution of metastatic cancers in response to therapy. Six patients with advanced breast, ovarian and lung cancers were followed over 1–2 years. For each case, exome sequencing was performed on 2–5 plasma samples (19 in total) spanning multiple courses of treatment, at selected time points when the allele fraction of tumour mutations in plasma was high, allowing improved sensitivity. For two cases, synchronous biopsies were also analysed, confirming genome-wide representation of the tumour genome in plasma. Quantification of allele fractions in plasma identified increased representation of mutant alleles in association with emergence of therapy resistance. These included an activating mutation in PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha) following treatment with paclitaxel; a truncating mutation in RB1 (retinoblastoma 1) following treatment with cisplatin; a truncating mutation in MED1 (mediator complex subunit 1) following treatment with tamoxifen and trastuzumab, and following subsequent treatment with lapatinib, a splicing mutation in GAS6 (growth arrest-specific 6) in the same patient; and a resistance-conferring mutation in EGFR (epidermal growth factor receptor; T790M) following treatment with gefitinib. These results establish proof of principle that exome-wide analysis of circulating tumour DNA could complement current invasive biopsy approaches to identify mutations associated with acquired drug resistance in advanced cancers. Serial analysis of cancer genomes in plasma constitutes a new paradigm for the study of clonal evolution in human cancers.

Noninvasive Identification and Monitoring of Cancer Mutations by Targeted Deep Sequencing of Plasma DNA

Sizable genomic regions were screened and low-frequency mutations were identified in circulating DNA of cancer patients using tagged-amplicon deep sequencing (TAm-Seq). Deep Sequencing Tumor DNA in Plasma Five liters of circulating blood contain millions of copies of the genome, broken into short fragments; in cancer patients, a small fraction is circulating tumor DNA (ctDNA). An even smaller number harbor mutations that affect cancer outcome. Looking for diagnostic answers in circulating DNA is a challenge, but Forshew, Murtaza, and colleagues have risen to the occasion by developing a tagged-amplicon deep sequencing (TAm-Seq) method that can amplify and sequence large genomic regions from even single copies of ctDNA. By sequencing such large regions, the authors were able to identify low-level mutations in the plasma of patients with high-grade serous ovarian carcinomas. Forshew et al. designed primers to amplify 5995 bases that covered select regions of cancer-related genes, including TP53, EGFR, BRAF, and KRAS. In plasma obtained from 38 patients with high levels of ctDNA, the authors were able to identify mutations in TP53 at allelic frequencies of 2% to 65%. In plasma samples from one patient, they also identified a de novo mutation in EGFR that had not been detected 15 months prior in the tumor mass itself. Finally, the TAm-Seq approach was used to sequence ctDNA in plasma samples collected from two women with ovarian cancer and one woman with breast cancer at different time points, tracking as many as 10 mutations in parallel. Forshew and coauthors showed that levels of mutant alleles reflected the clinical course of the disease and its treatment—for example, stabilized disease was associated with low allelic frequency, whereas patients at relapse exhibited a rise in frequency. Through several experiments, the authors were able to show that TAm-Seq is a viable method for sequencing large regions of ctDNA. Although this provides a new way to noninvasively identify gene mutations in our blood, TAm-Seq will need to achieve a more sensitive detection limit (<2% allele frequency) to identify mutations in the plasma of patients with less advanced cancers. Nevertheless, once optimized, this “liquid biopsy” approach will be amenable to personalized genomics, where the level and type of mutations in ctDNA would inform clinical decision-making on an individual basis. Plasma of cancer patients contains cell-free tumor DNA that carries information on tumor mutations and tumor burden. Individual mutations have been probed using allele-specific assays, but sequencing of entire genes to detect cancer mutations in circulating DNA has not been demonstrated. We developed a method for tagged-amplicon deep sequencing (TAm-Seq) and screened 5995 genomic bases for low-frequency mutations. Using this method, we identified cancer mutations present in circulating DNA at allele frequencies as low as 2%, with sensitivity and specificity of >97%. We identified mutations throughout the tumor suppressor gene TP53 in circulating DNA from 46 plasma samples of advanced ovarian cancer patients. We demonstrated use of TAm-Seq to noninvasively identify the origin of metastatic relapse in a patient with multiple primary tumors. In another case, we identified in plasma an EGFR mutation not found in an initial ovarian biopsy. We further used TAm-Seq to monitor tumor dynamics, and tracked 10 concomitant mutations in plasma of a metastatic breast cancer patient over 16 months. This low-cost, high-throughput method could facilitate analysis of circulating DNA as a noninvasive “liquid biopsy” for personalized cancer genomics.

Spatial and Temporal Heterogeneity in High-Grade Serous Ovarian Cancer: A Phylogenetic Analysis

Background The major clinical challenge in the treatment of high-grade serous ovarian cancer (HGSOC) is the development of progressive resistance to platinum-based chemotherapy. The objective of this study was to determine whether intra-tumour genetic heterogeneity resulting from clonal evolution and the emergence of subclonal tumour populations in HGSOC was associated with the development of resistant disease. Methods and Findings Evolutionary inference and phylogenetic quantification of heterogeneity was performed using the MEDICC algorithm on high-resolution whole genome copy number profiles and selected genome-wide sequencing of 135 spatially and temporally separated samples from 14 patients with HGSOC who received platinum-based chemotherapy. Samples were obtained from the clinical CTCR-OV03/04 studies, and patients were enrolled between 20 July 2007 and 22 October 2009. Median follow-up of the cohort was 31 mo (interquartile range 22–46 mo), censored after 26 October 2013. Outcome measures were overall survival (OS) and progression-free survival (PFS). There were marked differences in the degree of clonal expansion (CE) between patients (median 0.74, interquartile range 0.66–1.15), and dichotimization by median CE showed worse survival in CE-high cases (PFS 12.7 versus 10.1 mo, p = 0.009; OS 42.6 versus 23.5 mo, p = 0.003). Bootstrap analysis with resampling showed that the 95% confidence intervals for the hazard ratios for PFS and OS in the CE-high group were greater than 1.0. These data support a relationship between heterogeneity and survival but do not precisely determine its effect size. Relapsed tissue was available for two patients in the CE-high group, and phylogenetic analysis showed that the prevalent clonal population at clinical recurrence arose from early divergence events. A subclonal population marked by a NF1 deletion showed a progressive increase in tumour allele fraction during chemotherapy. Conclusions This study demonstrates that quantitative measures of intra-tumour heterogeneity may have predictive value for survival after chemotherapy treatment in HGSOC. Subclonal tumour populations are present in pre-treatment biopsies in HGSOC and can undergo expansion during chemotherapy, causing clinical relapse.

Ovarian Cancer Cell Line Panel (OCCP): Clinical Importance of In Vitro Morphological Subtypes

Epithelial ovarian cancer is a highly heterogeneous disease and remains the most lethal gynaecological malignancy in the Western world. Therapeutic approaches need to account for inter-patient and intra-tumoural heterogeneity and detailed characterization of in vitro models representing the different histological and molecular ovarian cancer subtypes is critical to enable reliable preclinical testing. There are approximately 100 publicly available ovarian cancer cell lines but their cellular and molecular characteristics are largely undescribed. We have characterized 39 ovarian cancer cell lines under uniform conditions for growth characteristics, mRNA/microRNA expression, exon sequencing, drug response for clinically-relevant therapeutics and collated all available information on the original clinical features and site of origin. We tested for statistical associations between the cellular and molecular features of the lines and clinical features. Of the 39 ovarian cancer cell lines, 14 were assigned as high-grade serous, four serous-type, one low-grade serous and 20 non-serous type. Three morphological subtypes: Epithelial (n = 21), Round (n = 7) and Spindle (n = 12) were identified that showed distinct biological and molecular characteristics, including overexpression of cell movement and migration-associated genes in the Spindle subtype. Comparison with the original clinical data showed association of the spindle-like tumours with metastasis, advanced stage, suboptimal debulking and poor prognosis. In addition, the expression profiles of Spindle, Round and Epithelial morphologies clustered with the previously described C1-stromal, C5-mesenchymal and C4 ovarian subtype expression profiles respectively. Comprehensive profiling of 39 ovarian cancer cell lines under controlled, uniform conditions demonstrates clinically relevant cellular and genomic characteristics. This data provides a rational basis for selecting models to develop specific treatment approaches for histological and molecular subtypes of ovarian cancer.

Cancer predisposing BARD1 mutations in breast–ovarian cancer families

The breast cancer susceptibility gene BARD1 (BRCA1-associated RING domain protein, MIM# 601593) acts with BRCA1 in DNA double-strand break (DSB) repair and also in apoptosis initiation. We screened 109 BRCA1/2 negative high-risk breast and/or ovarian cancer patients from North-Eastern Poland for BARD1 germline mutations using a combination of denaturing high-performance liquid chromatography and direct sequencing. We identified 16 different BARD1 sequence variants, five of which are novel. Three of them were suspected to be pathogenic, including a protein truncating nonsense mutation (c.1690C>T, p.Gln564X), a splice mutation (c.1315-2A>G) resulting in exon 5 skipping, and a silent change (c.1977A>G) which alters several exonic splicing enhancer motifs in exon 10 and results in a transcript lacking exons 2–9. Our findings suggest that BARD1 mutations may be regarded as cancer risk alleles and warrant further investigation to determine their actual contribution to non-BRCA1/2 breast and ovarian cancer families.

Optimized p53 immunohistochemistry is an accurate predictor of TP53 mutation in ovarian carcinoma

TP53 mutations are ubiquitous in high‐grade serous ovarian carcinomas (HGSOC), and the presence of TP53 mutation discriminates between high and low‐grade serous carcinomas and is now an important biomarker for clinical trials targeting mutant p53. p53 immunohistochemistry (IHC) is widely used as a surrogate for TP53 mutation but its accuracy has not been established. The objective of this study was to test whether improved methods for p53 IHC could reliably predict TP53 mutations independently identified by next generation sequencing (NGS). Four clinical p53 IHC assays and tagged‐amplicon NGS for TP53 were performed on 171 HGSOC and 80 endometrioid carcinomas (EC). p53 expression was scored as overexpression (OE), complete absence (CA), cytoplasmic (CY) or wild type (WT). p53 IHC was evaluated as a binary classifier where any abnormal staining predicted deleterious TP53 mutation and as a ternary classifier where OE, CA or WT staining predicted gain‐of‐function (GOF or nonsynonymous), loss‐of‐function (LOF including stopgain, indel, splicing) or no detectable TP53 mutations (NDM), respectively. Deleterious TP53 mutations were detected in 169/171 (99%) HGSOC and 7/80 (8.8%) EC. The overall accuracy for the best performing IHC assay for binary and ternary prediction was 0.94 and 0.91 respectively, which improved to 0.97 (sensitivity 0.96, specificity 1.00) and 0.95 after secondary analysis of discordant cases. The sensitivity for predicting LOF mutations was lower at 0.76 because p53 IHC detected mutant p53 protein in 13 HGSOC with LOF mutations. CY staining associated with LOF was seen in 4 (2.3%) of HGSOC. Optimized p53 IHC can approach 100% specificity for the presence of TP53 mutation and its high negative predictive value is clinically useful as it can exclude the possibility of a low‐grade serous tumour. 4.1% of HGSOC cases have detectable WT staining while harboring a TP53 LOF mutation, which limits sensitivity for binary prediction of mutation to 96%.

Exploratory Analysis of TP53 Mutations in Circulating Tumour DNA as Biomarkers of Treatment Response for Patients with Relapsed High-Grade Serous Ovarian Carcinoma: A Retrospective Study

Background Circulating tumour DNA (ctDNA) carrying tumour-specific sequence alterations may provide a minimally invasive means to dynamically assess tumour burden and response to treatment in cancer patients. Somatic TP53 mutations are a defining feature of high-grade serous ovarian carcinoma (HGSOC). We tested whether these mutations could be used as personalised markers to monitor tumour burden and early changes as a predictor of response and time to progression (TTP). Methods and Findings We performed a retrospective analysis of serial plasma samples collected during routine clinical visits from 40 patients with HGSOC undergoing heterogeneous standard of care treatment. Patient-specific TP53 assays were developed for 31 unique mutations identified in formalin-fixed paraffin-embedded tumour DNA from these patients. These assays were used to quantify ctDNA in 318 plasma samples using microfluidic digital PCR. The TP53 mutant allele fraction (TP53MAF) was compared to serum CA-125, the current gold-standard response marker for HGSOC in blood, as well as to disease volume on computed tomography scans by volumetric analysis. Changes after one cycle of treatment were compared with TTP. The median TP53MAF prior to treatment in 51 relapsed treatment courses was 8% (interquartile range [IQR] 1.2%–22%) compared to 0.7% (IQR 0.3%–2.0%) for seven untreated newly diagnosed stage IIIC/IV patients. TP53MAF correlated with volumetric measurements (Pearson r = 0.59, p < 0.001), and this correlation improved when patients with ascites were excluded (r = 0.82). The ratio of TP53MAF to volume of disease was higher in relapsed patients (0.04% per cm3) than in untreated patients (0.0008% per cm3, p = 0.004). In nearly all relapsed patients with disease volume > 32 cm3, ctDNA was detected at ≥20 amplifiable copies per millilitre of plasma. In 49 treatment courses for relapsed disease, pre-treatment TP53MAF concentration, but not CA-125, was associated with TTP. Response to chemotherapy was seen earlier with ctDNA, with a median time to nadir of 37 d (IQR 28–54) compared with a median time to nadir of 84 d (IQR 42–116) for CA-125. In 32 relapsed treatment courses evaluable for response after one cycle of chemotherapy, a decrease in TP53MAF of >60% was an independent predictor of TTP in multivariable analysis (hazard ratio 0.22, 95% CI 0.07–0.67, p = 0.008). Conversely, a decrease in TP53MAF of ≤60% was associated with poor response and identified cases with TTP < 6 mo with 71% sensitivity (95% CI 42%–92%) and 88% specificity (95% CI 64%–99%). Specificity was improved when patients with recent drainage of ascites were excluded. Ascites drainage led to a reduction of TP53MAF concentration. The limitations of this study include retrospective design, small sample size, and heterogeneity of treatment within the cohort. Conclusions In this retrospective study, we demonstrated that ctDNA is correlated with volume of disease at the start of treatment in women with HGSOC and that a decrease of ≤60% in TP53MAF after one cycle of chemotherapy was associated with shorter TTP. These results provide evidence that ctDNA has the potential to be a highly specific early molecular response marker in HGSOC and warrants further investigation in larger cohorts receiving uniform treatment.

A case of mast cell leukaemia with exon 9 KIT mutation and good response to imatinib

Background:  Mastocytosis is a myeloproliferative neoplasm characterized by the excessive proliferation of mast cells. Mast cell leukaemia (MCL), the aggressive form of this disease, requires cytoreductive therapy, such as cladribine, interferon‐alpha‐2b and, most recently, tyrosine kinase inhibitors – dasatinib or imatinib.

Germline whole exome sequencing and large-scale replication identifies FANCM as a likely high grade serous ovarian cancer susceptibility gene

We analyzed whole exome sequencing data in germline DNA from 412 high grade serous ovarian cancer (HGSOC) cases from The Cancer Genome Atlas Project and identified 5,517 genes harboring a predicted deleterious germline coding mutation in at least one HGSOC case. Gene-set enrichment analysis showed enrichment for genes involved in DNA repair (p = 1.8×10-3). Twelve DNA repair genes - APEX1, APLF, ATX, EME1, FANCL, FANCM, MAD2L2, PARP2, PARP3, POLN, RAD54L and SMUG1 – were prioritized for targeted sequencing in up to 3,107 HGSOC cases, 1,491 cases of other epithelial ovarian cancer (EOC) subtypes and 3,368 unaffected controls of European origin. We estimated mutation prevalence for each gene and tested for associations with disease risk. Mutations were identified in both cases and controls in all genes except MAD2L2, where we found no evidence of mutations in controls. In FANCM we observed a higher mutation frequency in HGSOC cases compared to controls (29/3,107 cases, 0.96 percent; 13/3,368 controls, 0.38 percent; P=0.008) with little evidence for association with other subtypes (6/1,491, 0.40 percent; P=0.82). The relative risk of HGSOC associated with deleterious FANCM mutations was estimated to be 2.5 (95% CI 1.3 – 5.0; P=0.006). In summary, whole exome sequencing of EOC cases with large-scale replication in case-control studies has identified FANCM as a likely novel susceptibility gene for HGSOC, with mutations associated with a moderate increase in risk. These data may have clinical implications for risk prediction and prevention approaches for high-grade serous ovarian cancer in the future and a significant impact on reducing disease mortality.We analyzed whole exome sequencing data in germline DNA from 412 high grade serous ovarian cancer (HGSOC) cases from The Cancer Genome Atlas Project and identified 5,517 genes harboring a predicted deleterious germline coding mutation in at least one HGSOC case. Gene-set enrichment analysis showed enrichment for genes involved in DNA repair (p = 1.8x10-3). Twelve DNA repair genes - APEX1, APLF, ATX, EME1, FANCL, FANCM, MAD2L2, PARP2, PARP3, POLN, RAD54L and SMUG1 - were prioritized for targeted sequencing in up to 3,107 HGSOC cases, 1,491 cases of other epithelial ovarian cancer (EOC) subtypes and 3,368 unaffected controls of European origin. We estimated mutation prevalence for each gene and tested for associations with disease risk. Mutations were identified in both cases and controls in all genes except MAD2L2, where we found no evidence of mutations in controls. In FANCM we observed a higher mutation frequency in HGSOC cases compared to controls (29/3,107 cases, 0.96 percent; 13/3,368 controls, 0.38 percent; P=0.008) with little evidence for association with other subtypes (6/1,491, 0.40 percent; P=0.82). The relative risk of HGSOC associated with deleterious FANCM mutations was estimated to be 2.5 (95% CI 1.3 - 5.0; P=0.006). In summary, whole exome sequencing of EOC cases with large-scale replication in case-control studies has identified FANCM as a likely novel susceptibility gene for HGSOC, with mutations associated with a moderate increase in risk. These data may have clinical implications for risk prediction and prevention approaches for high-grade serous ovarian cancer in the future and a significant impact on reducing disease mortality.

Selecting Short DNA Fragments In Plasma Improves Detection Of Circulating Tumour DNA

Non-invasive analysis of cancer genomes using cell-free circulating tumour DNA (ctDNA) is being widely implemented for clinical indications. The sensitivity for detecting the presence of ctDNA and genomic changes in ctDNA is limited by its low concentration compared to cell-free DNA of non-tumour origin. We studied the feasibility for enrichment of ctDNA by size selection, in plasma samples collected before and during chemotherapy treatment in 13 patients with recurrent high-grade serous ovarian cancer. We evaluated the effects using targeted and whole genome sequencing. Selecting DNA fragments between 90-150 bp before analysis yielded enrichment of mutated DNA fraction of up to 11-fold. This allowed identification of adverse copy number alterations, including MYC amplification, otherwise not observed. Size selection allows detection of tumour alterations masked by non-tumour DNA in plasma and could help overcome sensitivity limitations of liquid biopsy for applications in early diagnosis, detection of minimal residual disease, and genomic profiling.