Mutation rate and mutational signatures over disease course. A, X-axis represents samples, Y-axis represents mutations per MB. The mutation rate in pre-therapy samples is lower than that in post-therapy samples. B, X-axis represents 96 possible classes of mutations. Y-axis represents the fraction of these classes. The sample name is shown in each box.
Funding
DST | Science and Engineering Research Board (SERB)
Department of Biotechnology, Ministry of Science and Technology, India (DBT)
ARTICLE ABSTRACT
We report a deep next-generation sequencing analysis of 13 sequentially obtained tumor samples, eight sequentially obtained circulating tumor DNA (ctDNA) samples and three germline DNA samples over the life history of 3 patients with triple-negative breast cancer (TNBC), 2 of whom had germline pathogenic BRCA1 mutation, to unravel tumor evolution. Tumor tissue from all timepoints and germline DNA was subjected to whole-exome sequencing (WES), custom amplicon deep sequencing (30,000X) of a WES-derived somatic mutation panel, and SNP arrays for copy-number variation (CNV), while whole transcriptome sequencing (RNA-seq) was performed only on somatic tumor.There was enrichment of homologous recombination deficiency signature in all tumors and widespread CNV, which remained largely stable over time. Somatic tumor mutation numbers varied between patients and within each patient (range: 70–216, one outlier). There was minimal mutational overlap between patients with TP53 being the sole commonly mutated gene, but there was substantial overlap in sequential samples in each patient. Each patient's tumor contained a founding (“stem”) clone at diagnosis, which persisted over time, from which all other clones (“subclone”) were derived (“branching evolution”), which contained mutations in well-characterized cancer-related genes like PDGFRB, ARID2, TP53 (Patient_02), TP53, BRAF, BRIP1, CSF3R (Patient_04), and TP53, APC, EZH2 (Patient_07). Including stem and subclones, tumors from all patients were polyclonal at diagnosis and during disease progression. ctDNA recapitulated most tissue-derived stem clonal and subclonal mutations while detecting some additional subclonal mutations. RNA-seq revealed a stable basal-like pattern, with most highly expressed variants belonging to stem clone.
In germline BRCA1 mutated and BRCA wild-type patients, TNBC shows a branching evolutionary pattern of mutations with a single founding clone, are polyclonal throughout their disease course, and have widespread copy-number aberrations. This evolutionary pattern may be associated with treatment resistance or sensitivity and could be therapeutically exploited.
