Fast and accurate mutation detection in whole genome sequences of multiple isogenic samples with IsoMut

Mutational profile of the regenerative process and de novo genome assembly of the planarian Schmidtea polychroa

Planarians are organisms with a unique capacity to regenerate any part of their body. New tissues are generated in a process that requires many swift cell divisions. How costly is this process to an animal in terms of mutational load remains unknown. Using whole genome sequencing, we defined the mutational profile of the process of regeneration in the planarian species Schmidtea polychroa. We assembled de novo the genome of S. polychroa and analyzed mutations in animals that have undergone regeneration. We observed a threefold increase in the number of mutations and an altered mutational spectrum. High allele frequencies of subclonal mutations in regenerated animals suggested that most of the cells in the regenerated animal were descendants of a small number of stem cells with high expansion potential. We provide, for the first time, the draft genome assembly of S. polychroa, an estimation of the germline mutation rate for a planarian species and the mutational spectrum of the regeneration process of a living organism.

DNA mismatch repair protects the genome from oxygen-induced replicative mutagenesis

DNA mismatch repair (MMR) corrects mismatched DNA bases arising from multiple sources including polymerase errors and base damage. By detecting spontaneous mutagenesis using whole genome sequencing of cultured MMR deficient human cell lines, we show that a primary role of MMR is the repair of oxygen-induced mismatches. We found an approximately twofold higher mutation rate in MSH6 deficient DLD-1 cells or MHL1 deficient HCT116 cells exposed to atmospheric conditions as opposed to mild hypoxia, which correlated with oxidant levels measured using electron paramagnetic resonance spectroscopy. The oxygen-induced mutations were dominated by T to C base substitutions and single T deletions found primarily on the lagging strand. A broad sequence context preference, dependence on replication timing and a lack of transcriptional strand bias further suggested that oxygen-induced mutations arise from polymerase errors rather than oxidative base damage. We defined separate low and high oxygen-specific MMR deficiency mutation signatures common to the two cell lines and showed that the effect of oxygen is observable in MMR deficient cancer genomes, where it best correlates with the contribution of mutation signature SBS21. Our results imply that MMR corrects oxygen-induced genomic mismatches introduced by a replicative process in proliferating cells.

Spontaneous mutagenesis in human cells is controlled by REV1-Polymerase ζ and PRIMPOL

Translesion DNA synthesis (TLS) facilitates replication over damaged or difficult-to-replicate templates by employing specialized DNA polymerases. We investigate the effect on spontaneous mutagenesis of three main TLS control mechanisms: REV1 and PCNA ubiquitylation that recruit TLS polymerases and PRIMPOL that creates post-replicative gaps. Using whole-genome sequencing of cultured human RPE-1 cell clones, we find that REV1 and Polymerase ζ are wholly responsible for one component of base substitution mutagenesis that resembles homologous recombination deficiency, whereas the remaining component that approximates oxidative mutagenesis is reduced in PRIMPOL-/- cells. Small deletions in short repeats appear in REV1-/-PCNAK164R/K164R double mutants, revealing an alternative TLS mechanism. Also, 500-5,000 bp deletions appear in REV1-/- and REV3L-/- mutants, and chromosomal instability is detectable in REV1-/-PRIMPOL-/- cells. Our results indicate that TLS protects the genome from deletions and large rearrangements at the expense of being responsible for the majority of spontaneous base substitutions.

Temporal convolutional network for a Fast DNA mutation detection in breast cancer data

Early detection of breast cancer can be achieved through mutation detection in DNA sequences, which can be acquired through patient blood samples. Mutation detection can be performed using alignment and machine learning techniques. However, alignment techniques require reference sequences, and machine learning techniques still cannot predict index mutation and require supporting tools. Therefore, in this research, a Temporal Convolutional Network (TCN) model was proposed to detect the type and index mutation faster and without reference sequences and supporting tools. The architecture of the proposed TCN model is specifically designed for sequential labeling tasks on DNA sequence data. This allows for the detection of the mutation type of each nucleotide in the sequence, and if the nucleotide has a mutation, the index mutation can be obtained. The proposed model also uses 2-mers and 3-mers mapping techniques to improve detection performance. Based on the tests that have been carried out, the proposed TCN model can achieve the highest F1-score of 0.9443 for COSMIC dataset and 0.9629 for RSCM dataset, Additionally, the proposed TCN model can detect index mutation six times faster than BiLSTM model. Furthermore, the proposed model can detect type and index mutations based on the patient's DNA sequence, without the need for reference sequences or other additional tools.

CCAR2 functions downstream of the Shieldin complex to promote double-strand break end-joining

The 53BP1-RIF1 pathway restricts the resection of DNA double-strand breaks (DSBs) and promotes blunt end-ligation by non-homologous end joining (NHEJ) repair. The Shieldin complex is a downstream effector of the 53BP1-RIF1 pathway. Here, we identify a component of this pathway, CCAR2/DBC1, which is also required for restriction of DNA end-resection. CCAR2 co-immunoprecipitates with the Shieldin complex, and knockout of CCAR2 in a BRCA1-deficient cell line results in elevated DSB end-resection, RAD51 loading, and PARP inhibitor (PARPi) resistance. Knockout of CCAR2 is epistatic with knockout of other Shieldin proteins. The S1-like RNA-binding domain of CCAR2 is required for its interaction with the Shieldin complex and for suppression of DSB end-resection. CCAR2 functions downstream of the Shieldin complex, and CCAR2 knockout cells have delayed resolution of Shieldin complex foci. Forkhead-associated (FHA)-dependent targeting of CCAR2 to DSB sites re-sensitized BRCA1-/-SHLD2-/- cells to PARPi. Taken together, CCAR2 is a functional component of the 53BP1-RIF1 pathway, promotes the refill of resected DSBs, and suppresses homologous recombination.

Prospectively defined patterns of APOBEC3A mutagenesis are prevalent in human cancers

Mutational signatures defined by single base substitution (SBS) patterns in cancer have elucidated potential mutagenic processes that contribute to malignancy. Two prevalent mutational patterns in human cancers are attributed to the APOBEC3 cytidine deaminase enzymes. Among the seven human APOBEC3 proteins, APOBEC3A is a potent deaminase and proposed driver of cancer mutagenesis. In this study, we prospectively examine genome-wide aberrations by expressing human APOBEC3A in avian DT40 cells. From whole-genome sequencing, we detect hundreds to thousands of base substitutions per genome. The APOBEC3A signature includes widespread cytidine mutations and a unique insertion-deletion (indel) signature consisting largely of cytidine deletions. This multi-dimensional APOBEC3A signature is prevalent in human cancer genomes. Our data further reveal replication-associated mutations, the rate of stem-loop and clustered mutations, and deamination of methylated cytidines. This comprehensive signature of APOBEC3A mutagenesis is a tool for future studies and a potential biomarker for APOBEC3 activity in cancer.

Characterisation of the spectrum and genetic dependence of collateral mutations induced by translesion DNA synthesis

Translesion DNA synthesis (TLS) is a fundamental damage bypass pathway that utilises specialised polymerases with relaxed template specificity to achieve replication through damaged DNA. Misinsertions by low fidelity TLS polymerases may introduce additional mutations on undamaged DNA near the original lesion site, which we termed collateral mutations. In this study, we used whole genome sequencing datasets of chicken DT40 and several human cell lines to obtain evidence for collateral mutagenesis in higher eukaryotes. We found that cisplatin and UVC radiation frequently induce close mutation pairs within 25 base pairs that consist of an adduct-associated primary and a downstream collateral mutation, and genetically linked their formation to TLS activity involving PCNA ubiquitylation and polymerase κ. PCNA ubiquitylation was also indispensable for close mutation pairs observed amongst spontaneously arising base substitutions in cell lines with disrupted homologous recombination. Collateral mutation pairs were also found in melanoma genomes with evidence of UV exposure. We showed that collateral mutations frequently copy the upstream base, and extracted a base substitution signature that describes collateral mutagenesis in the presented dataset regardless of the primary mutagenic process. Using this mutation signature, we showed that collateral mutagenesis creates approximately 10–20% of non-paired substitutions as well, underscoring the importance of the process.

DNA base substitutions are the most common form of genomic mutations, formed both spontaneously and in response to environmental mutagens. One of the main mechanisms of base substitution mutagenesis is translesion synthesis, a process that relies on specialised DNA polymerases to replicate damaged DNA templates. In addition to incorrect base insertions at the site of lesions in the template, translesion polymerases may also generate ‘collateral’ mutations away from the lesion due to their lower accuracy in selecting the correct incoming nucleotide. In this study, we surveyed the whole genome sequence of experimental cell clones to examine the extent and genetic dependence of collateral mutagenesis in higher eukaryotes. Looking for close mutation pairs, we found that collateral mutations frequently occur near primary lesions generated by cisplatin or ultraviolet radiation in chicken and human cells, but are restricted to a short distance of approximately 25 base pairs. By analysing their sequence context, we showed that collateral mutations can also occur near correctly bypassed primary lesions and may be responsible for a considerable proportion of all base substitution mutations.

BRCA1 deficiency specific base substitution mutagenesis is dependent on translesion synthesis and regulated by 53BP1

Defects in BRCA1, BRCA2 and other genes of the homology-dependent DNA repair (HR) pathway cause an elevated rate of mutagenesis, eliciting specific mutation patterns including COSMIC signature SBS3. Using genome sequencing of knock-out cell lines we show that Y family translesion synthesis (TLS) polymerases contribute to the spontaneous generation of base substitution and short insertion/deletion mutations in BRCA1 deficient cells, and that TLS on DNA adducts is increased in BRCA1 and BRCA2 mutants. The inactivation of 53BP1 in BRCA1 mutant cells markedly reduces TLS-specific mutagenesis, and rescues the deficiency of template switch-mediated gene conversions in the immunoglobulin V locus of BRCA1 mutant chicken DT40 cells. 53BP1 also promotes TLS in human cellular extracts in vitro. Our results show that HR deficiency-specific mutagenesis is largely caused by TLS, and suggest a function for 53BP1 in regulating the choice between TLS and error-free template switching in replicative DNA damage bypass.

Epistasis, aneuploidy, and functional mutations underlie evolution of resistance to induced microtubule depolymerization

Cells with blocked microtubule polymerization are delayed in mitosis, but eventually manage to proliferate despite substantial chromosome missegregation. While several studies have analyzed the first cell division after microtubule depolymerization, we have asked how cells cope long-term with microtubule impairment. We allowed 24 clonal populations of yeast cells with beta-tubulin mutations preventing proper microtubule polymerization, to evolve for ˜150 generations. At the end of the laboratory evolution experiment, cells had regained the ability to form microtubules and were less sensitive to microtubule-depolymerizing drugs. Whole-genome sequencing identified recurrently mutated genes, in particular for tubulins and kinesins, as well as pervasive duplication of chromosome VIII. Recreating these mutations and chromosome VIII disomy prior to evolution confirmed that they allow cells to compensate for the original mutation in beta-tubulin. Most of the identified mutations did not abolish function, but rather restored microtubule functionality. Analysis of the temporal order of resistance development in independent populations repeatedly revealed the same series of events: disomy of chromosome VIII followed by a single additional adaptive mutation in either tubulins or kinesins. Since tubulins are highly conserved among eukaryotes, our results have implications for understanding resistance to microtubule-targeting drugs widely used in cancer therapy.

Evaluating Large Spontaneous Deletions in a Bovine Cell Line Selected for Bovine Viral Diarrhea Virus Resistance

Bovine viral diarrhea virus’s (BVDV) entry into bovine cells involves attachment of virions to cellular receptors, internalization, and pH-dependent fusion with endosomal membranes. The primary host receptor for BVDV is CD46; however, the complete set of host factors required for virus entry is unknown. The Madin-Darby bovine kidney (MDBK) cell line is susceptible to BVDV infection, while a derivative cell line (CRIB) is resistant at the level of virus entry. We performed complete genome sequencing of each to identify genomic variation underlying the resistant phenotype with the aim of identifying host factors essential for BVDV entry. Three large compound deletions in the BVDV-resistant CRIB cell line were identified and predicted to disrupt the function or expression of the genes PTPN12, GRID2, and RABGAP1L. However, CRISPR/Cas9 mediated knockout of these genes, individually or in combination, in the parental MDBK cell line did not impact virus entry or replication. Therefore, resistance to BVDV in the CRIB cell line is not due to the apparent spontaneous loss of PTPN12, GRID2, or RABGAP1L gene function. Identifying the functional cause of BVDV resistance in the CRIB cell line may require more detailed comparisons of the genomes and epigenomes.

A comparative analysis of the mutagenicity of platinum-containing chemotherapeutic agents reveals direct and indirect mutagenic mechanisms

Platinum-based drugs are a mainstay of cancer chemotherapy. However, their mutagenic effect can increase tumour heterogeneity, contribute to the evolution of treatment resistance and also induce secondary malignancies. We coupled whole genome sequencing with phenotypic investigations on two cell line models to compare the magnitude and examine the mechanism of mutagenicity of cisplatin, carboplatin and oxaliplatin. Cisplatin induced significantly more base substitution mutations than carboplatin or oxaliplatin when used at equitoxic concentrations on human TK6 or chicken DT40 cells, and also induced the highest number of short insertions and deletions. The analysis of base substitution spectra revealed that all three tested platinum drugs elicit both a direct mutagenic effect at purine dinucleotides, and an indirect effect of accelerating endogenous mutagenic processes, whereas the direct mutagenic effect appeared to correlate with the level of DNA damage caused as assessed through histone H2AX phosphorylation and single-cell agarose gel electrophoresis, the indirect mutagenic effects were equal. The different mutagenicity and DNA-damaging effect of equitoxic platinum drug treatments suggest that DNA damage independent mechanisms significantly contribute to their cytotoxicity. Thus, the comparatively high mutagenicity of cisplatin should be taken into account in the design of chemotherapeutic regimens.

Identification of a Synthetic Lethal Relationship between Nucleotide Excision Repair Deficiency and Irofulven Sensitivity in Urothelial Cancer

PURPOSE

Cisplatin-based chemotherapy is a first-line treatment for muscle-invasive and metastatic urothelial cancer. Approximately 10% of bladder urothelial tumors have a somatic missense mutation in the nucleotide excision repair (NER) gene, ERCC2, which confers increased sensitivity to cisplatin-based chemotherapy. However, a significant subset of patients is ineligible to receive cisplatin-based therapy due to medical contraindications, and no NER-targeted approaches are available for platinum-ineligible or platinum-refractory ERCC2-mutant cases.

EXPERIMENTAL DESIGN

We used a series of NER-proficient and NER-deficient preclinical tumor models to test sensitivity to irofulven, an abandoned anticancer agent. In addition, we used available clinical and sequencing data from multiple urothelial tumor cohorts to develop and validate a composite mutational signature of ERCC2 deficiency and cisplatin sensitivity.

RESULTS

We identified a novel synthetic lethal relationship between tumor NER deficiency and sensitivity to irofulven. Irofulven specifically targets cells with inactivation of the transcription-coupled NER (TC-NER) pathway and leads to robust responses in vitro and in vivo, including in models with acquired cisplatin resistance, while having minimal effect on cells with intact NER. We also found that a composite mutational signature of ERCC2 deficiency was strongly associated with cisplatin response in patients and was also associated with cisplatin and irofulven sensitivity in preclinical models.

CONCLUSIONS

Tumor NER deficiency confers sensitivity to irofulven, a previously abandoned anticancer agent, with minimal activity in NER-proficient cells. A composite mutational signature of NER deficiency may be useful in identifying patients likely to respond to NER-targeting agents, including cisplatin and irofulven.See related commentary by Jiang and Greenberg, p. 1833.

Two main mutational processes operate in the absence of DNA mismatch repair

The analysis of tumour genome sequences has demonstrated high rates of base substitution mutagenesis upon the inactivation of DNA mismatch repair (MMR), and the resulting somatic mutations in MMR deficient tumours appear to significantly enhance the response to immune therapy. A handful of different algorithmically derived base substitution mutation signatures have been attributed to MMR deficiency in tumour somatic mutation datasets. In contrast, mutation data obtained from whole genome sequences of isogenic wild type and MMR deficient cell lines in this study, as well as from published sources, show a more uniform experimental mutation spectrum of MMR deficiency. In order to resolve this discrepancy, we reanalysed mutation data from MMR deficient tumour whole exome and whole genome sequences. We derived two base substitution signatures using non-negative matrix factorisation, which together adequately describe mutagenesis in all tumour and cell line samples. The two new signatures broadly resemble COSMIC signatures 6 and 20, but perform better than existing COSMIC signatures at identifying MMR deficient tumours in mutation signature deconstruction. We show that the contribution of the two identified signatures, one of which is dominated by C to T mutations at CpG sites, is biased by the different sequence composition of the exome and the whole genome. We further show that the identity of the inactivated MMR gene, the tissue type, the mutational burden or the patient's age does not influence the mutation spectrum, but that a tendency for a greater contribution by the CpG mutational process is observed in tumours as compared to cultured cells. Our analysis suggest that two separable mutational processes operate in the genomes of MMR deficient cells.

Correlation of homologous recombination deficiency induced mutational signatures with sensitivity to PARP inhibitors and cytotoxic agents

BACKGROUND

Homologous recombination (HR) repair deficiency arising from defects in BRCA1 or BRCA2 is associated with characteristic patterns of somatic mutations. In this genetic study, we ask whether inactivating mutations in further genes of the HR pathway or the DNA damage checkpoint also give rise to somatic mutation patterns that can be used for treatment prediction.

RESULTS

Using whole genome sequencing of an isogenic knockout cell line panel, we find a universal HR deficiency-specific base substitution signature that is similar to COSMIC signature 3. In contrast, we detect different deletion phenotypes corresponding to specific HR mutants. The inactivation of BRCA2 or PALB2 leads to larger deletions, typically with microhomology, when compared to the disruption of BRCA1, RAD51 paralogs, or RAD54. Comparison with the deletion spectrum of Cas9 cut sites suggests that most spontaneously arising genomic deletions are not the consequence of double-strand breaks. Surprisingly, the inactivation of checkpoint kinases ATM and CHK2 has no mutagenic consequences. Analysis of tumor exomes with biallelic inactivating mutations in the investigated genes confirms the validity of the cell line models. We present a comprehensive analysis of sensitivity of the investigated mutants to 13 therapeutic agents for the purpose of correlating genomic mutagenic phenotypes with drug sensitivity.

CONCLUSION

Our results suggest that no single genomic mutational class shows perfect correlation with sensitivity to common treatments, but the contribution of COSMIC signature 3 to base substitutions, or a combined measure of different features, may be reasonably good at predicting platinum and PARP inhibitor sensitivity.

Single-Cell Omics Analyses Enabled by Microchip Technologies

Single-cell omics studies provide unique information regarding cellular heterogeneity at various levels of the molecular biology central dogma. This knowledge facilitates a deeper understanding of how underlying molecular and architectural changes alter cell behavior, development, and disease processes. The emerging microchip-based tools for single-cell omics analysis are enabling the evaluation of cellular omics with high throughput, improved sensitivity, and reduced cost. We review state-of-the-art microchip platforms for profiling genomics, epigenomics, transcriptomics, proteomics, metabolomics, and multi-omics at single-cell resolution. We also discuss the background of and challenges in the analysis of each molecular layer and integration of multiple levels of omics data, as well as how microchip-based methodologies benefit these fields. Additionally, we examine the advantages and limitations of these approaches. Looking forward, we describe additional challenges and future opportunities that will facilitate the improvement and broad adoption of single-cell omics in life science and medicine.

Long-term treatment with the PARP inhibitor niraparib does not increase the mutation load in cell line models and tumour xenografts

Background

Poly-ADP ribose polymerase (PARP) inhibitor-based cancer therapy selectively targets cells with deficient homologous recombination repair. Considering their long-term use in maintenance treatment, any potential mutagenic effect of PARP inhibitor treatment could accelerate the development of resistance or harm non-malignant somatic cells.

Methods

We tested the mutagenicity of long-term treatment with the PARP inhibitor niraparib using whole-genome sequencing of cultured cell clones and whole-exome sequencing of patient-derived breast cancer xenografts.

Results

We observed no significant increase in the number and alteration in the spectrum of base substitutions, short insertions and deletions and genomic rearrangements upon niraparib treatment of human DLD-1 colon adenocarcinoma cells, wild-type and BRCA1 mutant chicken DT40 lymphoblastoma cells and BRCA1-defective SUM149PT breast carcinoma cells, except for a minor increase in specific deletion classes. We also did not detect any contribution of in vivo niraparib treatment to subclonal mutations arising in breast cancer-derived xenografts.

Conclusions

The results suggest that long-term inhibition of DNA repair with PARP inhibitors has no or only limited mutagenic effect. Mutagenesis due to prolonged use of PARP inhibitors in cancer treatment is therefore not expected to contribute to the genetic evolution of resistance, generate significant immunogenic neoepitopes or induce secondary malignancies.

A genetic study based on PCNA-ubiquitin fusions reveals no requirement for PCNA polyubiquitylation in DNA damage tolerance

Post-translational modifications of Proliferating Cell Nuclear Antigen (PCNA) play a key role in regulating the bypass of DNA lesions during DNA replication. PCNA can be monoubiquitylated at lysine 164 by the RAD6-RAD18 ubiquitin ligase complex. Through this modification, PCNA can interact with low fidelity Y family DNA polymerases to promote translesion synthesis. Monoubiquitylated PCNA can be polyubiquitylated on lysine 63 of ubiquitin by a further ubiquitin-conjugating complex. This modification promotes a template switching bypass process in yeast, while its role in higher eukaryotes is less clear. We investigated the function of PCNA ubiquitylation using a PCNA^K164R mutant DT40 chicken B lymphoblastoma cell line, which is hypersensitive to DNA damaging agents such as methyl methanesulfonate (MMS), cisplatin or ultraviolet radiation (UV) due to the loss of PCNA modifications. In the PCNA^K164R mutant we also detected cell cycle arrest following UV treatment, a reduced rate of damage bypass through translesion DNA synthesis on synthetic UV photoproducts, and an increased rate of genomic mutagenesis following MMS treatment. PCNA-ubiquitin fusion proteins have been reported to mimic endogenous PCNA ubiquitylation. We found that the stable expression of a PCNA^K164R-ubiquitin fusion protein fully or partially rescued the observed defects of the PCNA^K164R mutant. The expression of a PCNA^K164R-ubiquitin^K63R fusion protein, on which the formation of lysine 63-linked polyubiquitin chains is not possible, similarly rescued the cell cycle arrest, DNA damage sensitivity, reduction of translesion synthesis and increase of MMS-induced genomic mutagenesis. Template switching bypass was not affected by the genetic elimination of PCNA polyubiquitylation, but it was reduced in the absence of the recombination proteins BRCA1 or XRCC3. Our study found no requirement for PCNA polyubiquitylation to protect cells from replication-stalling DNA damage.

Loss of BRCA1 or BRCA2 markedly increases the rate of base substitution mutagenesis and has distinct effects on genomic deletions

Loss-of-function mutations in the BRCA1 and BRCA2 genes increase the risk of cancer. Owing to their function in homologous recombination repair, much research has focused on the unstable genomic phenotype of BRCA1/2 mutant cells manifest mainly as large-scale rearrangements. We used whole-genome sequencing of multiple isogenic chicken DT40 cell clones to precisely determine the consequences of BRCA1/2 loss on all types of genomic mutagenesis. Spontaneous base substitution mutation rates increased sevenfold upon the disruption of either BRCA1 or BRCA2, and the arising mutation spectra showed strong and specific correlation with a mutation signature associated with BRCA1/2 mutant tumours. To model endogenous alkylating damage, we determined the mutation spectrum caused by methyl methanesulfonate (MMS), and showed that MMS also induces more base substitution mutations in BRCA1/2-deficient cells. Spontaneously arising and MMS-induced insertion/deletion mutations and large rearrangements were also more common in BRCA1/2 mutant cells compared with the wild-type control. A difference in the short deletion phenotypes of BRCA1 and BRCA2 suggested distinct roles for the two proteins in the processing of DNA lesions, as BRCA2 mutants contained more short deletions, with a wider size distribution, which frequently showed microhomology near the breakpoints resembling repair by non-homologous end joining. An increased and prolonged gamma-H2AX signal in MMS-treated BRCA1/2 cells suggested an aberrant processing of stalled replication forks as the cause of increased mutagenesis. The high rate of base substitution mutagenesis demonstrated by our experiments is likely to significantly contribute to the oncogenic effect of the inactivation of BRCA1 or BRCA2.

A comprehensive survey of the mutagenic impact of common cancer cytotoxics

BACKGROUND

Genomic mutations caused by cytotoxic agents used in cancer chemotherapy may cause secondary malignancies as well as contribute to the evolution of treatment-resistant tumour cells. The stable diploid genome of the chicken DT40 lymphoblast cell line, an established DNA repair model system, is well suited to accurately assay genomic mutations.

RESULTS

We use whole genome sequencing of multiple DT40 clones to determine the mutagenic effect of eight common cytotoxics used for the treatment of millions of patients worldwide. We determine the spontaneous mutagenesis rate at 2.3 × 10(-10) per base per cell division and find that cisplatin, cyclophosphamide and etoposide induce extra base substitutions with distinct spectra. After four cycles of exposure, cisplatin induces 0.8 mutations per Mb, equivalent to the median mutational burden in common leukaemias. Cisplatin-induced mutations, including short insertions and deletions, are mainly located at sites of putative intrastrand crosslinks. We find two of the newly defined cisplatin-specific mutation types as causes of the reversion of BRCA2 mutations in emerging cisplatin-resistant tumours or cell clones. Gemcitabine, 5-fluorouracil, hydroxyurea, doxorubicin and paclitaxel have no measurable mutagenic effect. The cisplatin-induced mutation spectrum shows good correlation with cancer mutation signatures attributed to smoking and other sources of guanine-directed base damage.

CONCLUSION

This study provides support for the use of cell line mutagenesis assays to validate or predict the mutagenic effect of environmental and iatrogenic exposures. Our results suggest genetic reversion due to cisplatin-induced mutations as a distinct mechanism for developing resistance.