Comprehensive investigation of the mutagenic potential of six pesticides classified by IARC as probably carcinogenic to humans.

Pesticides are significant environmental pollutants, and many of them possess mutagenic potential, which is closely linked to carcinogenesis. Here we tested the mutagenicity of all six pesticides classified probably carcinogenic (Group 2A) by the International Agency of Research on Cancer: 4,4'-DDT, captafol, dieldrin, diazinon, glyphosate and malathion. Whole genome sequencing of TK6 human lymphoblastoid cell clones following 30-day exposure at subtoxic concentrations revealed a clear mutagenic effect of treatment with captafol or malathion when added at 200 nM or 100 µM initial concentrations, respectively. Each pesticide induced a specific base substitution mutational signature: captafol increased C to A mutations primarily, while malathion induced mostly C to T mutations. 4,4'-DDT, dieldrin, diazinon and glyphosate were not mutagenic. Whereas captafol induced chromosomal instability, H2A.X phosphorylation and cell cycle arrest in G2/M phase, all indicating DNA damage, malathion did not induce DNA damage markers or cell cycle alterations despite its mutagenic effect. Hypersensitivity of REV1 and XPA mutant DT40 chicken cell lines suggests that captafol induces DNA adducts that are bypassed by translesion DNA synthesis and are targets for nucleotide excision repair. The experimentally identified mutational signatures of captafol and malathion could shed light on the mechanism of action of these compounds. The signatures are potentially suitable for detecting past exposure in tumour samples, but the reanalysis of large cancer genome databases did not reveal any evidence of captafol or malathion exposure.

CYP1A2 expression rather than genotype is associated with olanzapine concentration in psychiatric patients.

Olanzapine is a commonly prescribed atypical antipsychotic agent for treatment of patients with schizophrenia and bipolar disorders. Previous in vitro studies using human liver microsomes identified CYP1A2 and CYP2D6 enzymes being responsible for CYP-mediated metabolism of olanzapine. The present work focused on the impact of CYP1A2 and CYP2D6 genetic polymorphisms as well as of CYP1A2 metabolizing capacity influenced by non-genetic factors (sex, age, smoking) on olanzapine blood concentration in patients with psychiatric disorders (N = 139). CYP2D6 genotype-based phenotype appeared to have negligible contribution to olanzapine metabolism, whereas a dominant role of CYP1A2 in olanzapine exposure was confirmed. However, CYP1A2 expression rather than CYP1A2 genetic variability was demonstrated to be associated with olanzapine concentration in patients. Significant contribution of - 163C > A (rs762551), the most common SNP (single nucleotide polymorphism) in CYP1A2 gene, to enhanced inducibility was confirmed by an increase in CYP1A2 mRNA expression in smokers carrying - 163A, and smoking was found to have appreciable impact on olanzapine concentration normalized by the dose/bodyweight. Furthermore, patients' olanzapine exposure was in strong association with CYP1A2 expression; therefore, assaying CYP1A2 mRNA level in leukocytes can be an appropriate tool for the estimation of patients' olanzapine metabolizing capacity and may be relevant in optimizing olanzapine dosage.

Potential Association of Cytochrome P450 Copy Number Alteration in Tumour with Chemotherapy Resistance in Lung Adenocarcinoma Patients.

Resistance to anticancer agents is a major obstacle to efficacious tumour therapy and responsible for high cancer-related mortality rates. Some resistance mechanisms are associated with pharmacokinetic variability in anticancer drug exposure due to genetic polymorphisms of drug-metabolizing cytochrome P450 (CYP) enzymes, whereas variations in tumoural metabolism as a consequence of CYP copy number alterations are assumed to contribute to the selection of resistant cells. A high-throughput quantitative polymerase chain reaction (qPCR)-based method was developed for detection of CYP copy number alterations in tumours, and a scoring system improved the identification of inappropriate reference genes that underwent deletion/multiplication in tumours. The copy numbers of both the target (CYP2C8, CYP3A4) and the reference genes (ALB, B2M, BCKDHA, F5, CD36, MPO, TBP, RPPH1) established in primary lung adenocarcinoma by the qPCR-based method were congruent with those determined by next-generation sequencing (for 10 genes, slope = 0.9498, r2 = 0.72). In treatment naïve adenocarcinoma samples, the copy number multiplication of paclitaxel-metabolizing CYP2C8 and/or CYP3A4 was more prevalent in non-responder patients with progressive disease/exit than in responders with complete remission. The high-throughput qPCR-based method can become an alternative approach to next-generation sequencing in routine clinical practice, and identification of altered CYP copy numbers may provide a promising biomarker for therapy-resistant tumours.

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.

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.

A non-catalytic herpesviral protein reconfigures ERK-RSK signaling by targeting kinase docking systems in the host.

The Kaposi's sarcoma associated herpesvirus protein ORF45 binds the extracellular signal-regulated kinase (ERK) and the p90 Ribosomal S6 kinase (RSK). ORF45 was shown to be a kinase activator in cells but a kinase inhibitor in vitro, and its effects on the ERK-RSK complex are unknown. Here, we demonstrate that ORF45 binds ERK and RSK using optimized linear binding motifs. The crystal structure of the ORF45-ERK2 complex shows how kinase docking motifs recognize the activated form of ERK. The crystal structure of the ORF45-RSK2 complex reveals an AGC kinase docking system, for which we provide evidence that it is functional in the host. We find that ORF45 manipulates ERK-RSK signaling by favoring the formation of a complex, in which activated kinases are better protected from phosphatases and docking motif-independent RSK substrate phosphorylation is selectively up-regulated. As such, our data suggest that ORF45 interferes with the natural design of kinase docking systems in the host.

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.

Characterization of the Intramolecular Interactions and Regulatory Mechanisms of the Scaffold Protein Tks4.

The scaffold protein Tks4 is a member of the p47phox-related organizer superfamily. It plays a key role in cell motility by being essential for the formation of podosomes and invadopodia. In addition, Tks4 is involved in the epidermal growth factor (EGF) signaling pathway, in which EGF induces the translocation of Tks4 from the cytoplasm to the plasma membrane. The evolutionarily-related protein p47phox and Tks4 share many similarities in their N-terminal region: a phosphoinositide-binding PX domain is followed by two SH3 domains (so called "tandem SH3") and a proline-rich region (PRR). In p47phox, the PRR is followed by a relatively short, disordered C-terminal tail region containing multiple phosphorylation sites. These play a key role in the regulation of the protein. In Tks4, the PRR is followed by a third and a fourth SH3 domain connected by a long (~420 residues) unstructured region. In p47phox, the tandem SH3 domain binds the PRR while the first SH3 domain interacts with the PX domain, thereby preventing its binding to the membrane. Based on the conserved structural features of p47phox and Tks4 and the fact that an intramolecular interaction between the third SH3 and the PX domains of Tks4 has already been reported, we hypothesized that Tks4 is similarly regulated by autoinhibition. In this study, we showed, via fluorescence-based titrations, MST, ITC, and SAXS measurements, that the tandem SH3 domain of Tks4 binds the PRR and that the PX domain interacts with the third SH3 domain. We also investigated a phosphomimicking Thr-to-Glu point mutation in the PRR as a possible regulator of intramolecular interactions. Phosphatidylinositol-3-phosphate (PtdIns(3)P) was identified as the main binding partner of the PX domain via lipid-binding assays. In truncated Tks4 fragments, the presence of the tandem SH3, together with the PRR, reduced PtdIns(3)P binding, while the presence of the third SH3 domain led to complete inhibition.

Peptide Based Inhibitors of Protein Binding to the Mitogen-Activated Protein Kinase Docking Groove.

Mitogen-activated protein kinases (MAPK) are important regulatory units in cells and they take part in the regulation of many cellular functions such as cell division, differentiation or apoptosis. All MAPKs have a shallow docking groove that interacts with linear binding motifs of their substrate proteins and their regulatory proteins such as kinases, phosphatases, scaffolds. Inhibition of these protein-protein interactions may reduce or abolish the activity of the targeted kinase. Based on the wide range of their biological activity, this kind of inhibition can be useful in the treatment of many disorders like tumors, inflammation or undesired cell apoptosis. In this study a linear binding motif from the RHDF1 protein-a 15 amino acids long peptide-was selected for optimization to increase its cellular uptake but retaining its low micromolar binding affinity. First, we synthesized an octaarginine conjugate that showed efficient cellular uptake. Next, we set out to reduce the size of this construct. We were able to decrease the length of the original peptide, and to increase its cellular uptake with specific chemical modifications. These new constructs bound better to ERK2 and p38 kinases than the original peptide and they showed markedly increased cellular uptake. The new octaarginine conjugate and one of the minimized bicyclic derivatives could inhibit the phosphorylation of intracellular ERK or p38. However, the modulation of MAPK phosphorylation levels by these cell-penetrating peptides were complex, despite that in biochemical assays they all inhibited MAPK-substrate binding as well as phosphorylation. The optimized peptides depending on the applied concentration caused an expected decrease, but also some unexpected increase in MAPK phosphorylation patterns in the cell. This possibly reflects the complexity of MAPK docking groove mediated protein-protein interactions including bone fide MAPK clients such activator kinases, deactivating phosphatases or regulatory scaffolds. Thus, our findings with optimized cell-penetrating "inhibitory" peptides highlight the opportunities but also the pitfalls of docking peptide based MAPK activity regulation and call for a better quantitative understanding of MAPK mediated protein-protein interactions in cells.

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.

Evaluation and modulation of DNA lesion bypass in an SV40 large T antigen-based in vitro replication system.

DNA damage removal by nucleotide excision repair (NER) and replicative bypass via translesion synthesis (TLS) and template switch (TSw) are important in ensuring genome stability. In this study, we tested the applicability of an SV40 large T antigen-based replication system for the simultaneous examination of these damage tolerance processes. Using both Sanger and next-generation sequencing combined with lesion-specific qPCR and replication efficiency studies, we demonstrate that this system works well for studying NER and TLS, especially its one-polymerase branch, while it is less suited to investigations of homology-related repair processes, such as TSw. Cis-syn cyclobutane pyrimidine dimer photoproducts were replicated with equal efficiency to lesion-free plasmids in vitro, and the majority of TLS on this lesion could be inhibited by a peptide (PIR) specific for the polη-PCNA interaction interface. TLS on 6-4 pyrimidine-pyrimidone photoproduct proved to be inefficient and was slightly facilitated by PIR as well as by a recombinant ubiquitin-binding zinc finger domain of polη in HeLa extract, possibly by promoting polymerase exchange. Supplementation of the extract with recombinant PCNA variants indicated the dependence of TLS on PCNA ubiquitylation. In contrast to active TLS and NER, we found no evidence of successful TSw in cellular extracts. The established methods can promote in vitro investigations of replicative DNA damage bypass.

Genome-wide alterations of uracil distribution patterns in human DNA upon chemotherapeutic treatments.

Numerous anti-cancer drugs perturb thymidylate biosynthesis and lead to genomic uracil incorporation contributing to their antiproliferative effect. Still, it is not yet characterized if uracil incorporations have any positional preference. Here, we aimed to uncover genome-wide alterations in uracil pattern upon drug treatments in human cancer cell line models derived from HCT116. We developed a straightforward U-DNA sequencing method (U-DNA-Seq) that was combined with in situ super-resolution imaging. Using a novel robust analysis pipeline, we found broad regions with elevated probability of uracil occurrence both in treated and non-treated cells. Correlation with chromatin markers and other genomic features shows that non-treated cells possess uracil in the late replicating constitutive heterochromatic regions, while drug treatment induced a shift of incorporated uracil towards segments that are normally more active/functional. Data were corroborated by colocalization studies via dSTORM microscopy. This approach can be applied to study the dynamic spatio-temporal nature of genomic uracil.

Establishment and Characterization of a Brca1-/-, p53-/- Mouse Mammary Tumor Cell Line.

Breast cancer is the most commonly occurring cancer in women and the second most common cancer overall. By the age of 80, the estimated risk for breast cancer for women with germline BRCA1 or BRCA2 mutations is around 80%. Genetically engineered BRCA1-deficient mouse models offer a unique opportunity to study the pathogenesis and therapy of triple negative breast cancer. Here we present a newly established Brca1-/-, p53-/- mouse mammary tumor cell line, designated as CST. CST shows prominent features of BRCA1-mutated triple-negative breast cancers including increased motility, high proliferation rate, genome instability and sensitivity to platinum chemotherapy and PARP inhibitors (olaparib, veliparib, rucaparib and talazoparib). Genomic instability of CST cells was confirmed by whole genome sequencing, which also revealed the presence of COSMIC (Catalogue of Somatic Mutations in Cancer) mutation signatures 3 and 8 associated with homologous recombination (HR) deficiency. In vitro sensitivity of CST cells was tested against 11 chemotherapy agents. Tumors derived from orthotopically injected CST-mCherry cells in FVB-GFP mice showed sensitivity to cisplatin, providing a new model to study the cooperation of BRCA1-KO, mCherry-positive tumor cells and the GFP-expressing stromal compartment in therapy resistance and metastasis formation. In summary, we have established CST cells as a new model recapitulating major characteristics of BRCA1-negative breast cancers.

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.

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 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.

The genome of the chicken DT40 bursal lymphoma cell line.

The chicken DT40 cell line is a widely used model system in the study of multiple cellular processes due to the efficiency of homologous gene targeting. The cell line was derived from a bursal lymphoma induced by avian leukosis virus infection. In this study we characterized the genome of the cell line using whole genome shotgun sequencing and single nucleotide polymorphism array hybridization. The results indicate that wild-type DT40 has a relatively normal karyotype, except for whole chromosome copy number gains, and no karyotype variability within stocks. In a comparison to two domestic chicken genomes and the Gallus gallus reference genome, we found no unique mutational processes shaping the DT40 genome except for a mild increase in insertion and deletion events, particularly deletions at tandem repeats. We mapped coding sequence mutations that are unique to the DT40 genome; mutations inactivating the PIK3R1 and ATRX genes likely contributed to the oncogenic transformation. In addition to a known avian leukosis virus integration in the MYC gene, we detected further integration sites that are likely to de-regulate gene expression. The new findings support the hypothesis that DT40 is a typical transformed cell line with a relatively intact genome; therefore, it is well-suited to the role of a model system for DNA repair and related processes. The sequence data generated by this study, including a searchable de novo genome assembly and annotated lists of mutated genes, will support future research using this cell line.