The mutagenic consequences of defective DNA repair.

Multiple separate repair mechanisms safeguard the genome against various types of DNA damage, and their failure can increase the rate of spontaneous mutagenesis. The malfunction of distinct repair mechanisms leads to genomic instability through different mutagenic processes. For example, defective mismatch repair causes high base substitution rates and microsatellite instability, whereas homologous recombination deficiency is characteristically associated with deletions and chromosome instability. This review presents a comprehensive collection of all mutagenic phenotypes associated with the loss of each DNA repair mechanism, drawing on data from a variety of model organisms and mutagenesis assays, and placing greatest emphasis on systematic analyses of human cancer datasets. We describe the latest theories on the mechanism of each mutagenic process, often explained by reliance on an alternative repair pathway or the error-prone replication of unrepaired, damaged DNA. Aided by the concept of mutational signatures, the genomic phenotypes can be used in cancer diagnosis to identify defective DNA repair pathways.

Increased frequency of CHD1 deletions in prostate cancers of African American men is associated with rapid disease progression without inducing homologous recombination deficiency.

We analyzed genomic data derived from the prostate cancer of African and European American men in order to identify differences that may contribute to racial disparity of outcome and that could also define novel therapeutic strategies. In addition to analyzing patient derived next generation sequencing data, we performed FISH based confirmatory studies of Chromodomain helicase DNA-binding protein 1 (CHD1) loss on prostate cancer tissue microarrays. We created CRISPR edited, CHD1 deficient prostate cancer cell lines for genomic, drug sensitivity and functional homologous recombination (HR) activity analysis. We found that subclonal deletion of CHD1 is nearly three times as frequent in prostate tumors of African American men than in men of European ancestry and it associates with rapid disease progression. We further showed that CHD1 deletion is not associated with homologous recombination deficiency associated mutational signatures in prostate cancer. In prostate cancer cell line models CHD1 deletion did not induce HR deficiency as detected by RAD51 foci formation assay or mutational signatures, which was consistent with the moderate increase of olaparib sensitivity. CHD1 deficient prostate cancer cells, however, showed higher sensitivity to talazoparib. CHD1 loss may contribute to worse outcome of prostate cancer in African American men. A deeper understanding of the interaction between CHD1 loss and PARP inhibitor sensitivity will be needed to determine the optimal use of targeted agents such as talazoparib in the context of castration resistant prostate cancer.

Probing the biological consequences of a previously undescribed de novo mutation of ZMYND11 in a schizophrenia patient by CRISPR genome editing and induced pluripotent stem cell based in vitro disease-modeling.

BACKGROUND: Schizophrenia (SCZ) is a severe neuropsychiatric disorder of complex, poorly understood etiology, associated with both genetic and environmental factors. De novo mutations (DNMs) represent a new source of genetic variation in SCZ, however, in most cases their biological significance remains unclear. We sought to investigate molecular disease pathways connected to DNMs in SCZ by combining human induced pluripotent stem cell (hiPSC) based disease modeling and CRISPR-based genome editing. METHODS: We selected a SCZ case-parent trio with the case individual carrying a potentially disease causing 1495C > T nonsense DNM in the zinc finger MYND domain-containing protein 11 (ZMYND11), a gene implicated in biological processes relevant for SCZ. In the patient-derived hiPSC line the mutation was corrected using CRISPR, while monoallelic or biallelic frameshift mutations were introduced into a control hiPSC line. Isogenic cell lines were differentiated into hippocampal neuronal progenitor cells (NPCs) and functionally active dentate gyrus granule cells (DGGCs). Immunofluorescence microscopy and RNA sequencing were used to test for morphological and transcriptomic differences at NPC and DGCC stages. Functionality of neurons was investigated using calcium-imaging and multi-electrode array measurements. RESULTS: Morphology in the mutant hippocampal NPCs and neurons was preserved, however, we detected significant transcriptomic and functional alterations. RNA sequencing showed massive upregulation of neuronal differentiation genes, and downregulation of cell adhesion genes. Decreased reactivity to glutamate was demonstrated by calcium-imaging. CONCLUSIONS: Our findings lend support to the involvement of glutamatergic dysregulation in the pathogenesis of SCZ. This approach represents a powerful model system for precision psychiatry and pharmacological research.

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.

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.

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.

Probing telomeric-like G4 structures with full or partial 2'-deoxy-5-hydroxyuridine substitutions.

Guanine quadruplexes (G4s) are stable four-stranded secondary DNA structures held together by noncanonical G-G base tetrads. We synthesised the nucleoside analogue 2'-deoxy-5-hydroxyuridine (H) and inserted its phosphoramidite into telomeric repeat-type model oligonucleotides. Full and partial substitutions were made, replacing all guanines in all the three tetrads of a three-tier G4 structure, or only in the putative upper, central, or lower tetrads. We characterised these modified structures using CD, UV absorbance spectroscopy, native gel studies, and a capture oligo-based G4 disruption kinetic assay. The strand separation activity of BLM helicase on these substituted structures was also investigated. Two of the partially H-substituted constructs adopted G4-like structures, but displayed lower thermal stabilities compared to unsubstituted G4. The construct modified in its central tetrad remained mostly denatured, but the possibility of a special structure for the fully replaced variant remained open. H substitutions did not interfere with the G4-resolving activity of BLM helicase, but its efficiency was highly influenced by construct topology and even more by the G4 ligand PhenDC3. Our results suggest that the H modification can be incorporated into G quadruplexes, but only at certain positions to maintain G4 stability. The destabilizing effect observed for 2'-deoxy-5-hydroxyuridine indicates that the cytosine deamination product 5-hydroxyuracil and its nucleoside counterpart in RNA (5-hydroxyuridine), might also be destabilizing in cellular DNA and RNA quadruplexes. The kinetic assay employed in this study can be generally employed for a fast comparison of the stabilities of various G4s either in their free or ligand-bound states.

A fresh look at somatic mutations in cancer.

Analysis of cancer genome sequences reveals new mutational signatures.

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.

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.

BRCA1/Trp53 heterozygosity and replication stress drive esophageal cancer development in a mouse model.

BRCA1 germline mutations are associated with an increased risk of breast and ovarian cancer. Recent findings of others suggest that BRCA1 mutation carriers also bear an increased risk of esophageal and gastric cancer. Here, we employ a Brca1/Trp53 mouse model to show that unresolved replication stress (RS) in BRCA1 heterozygous cells drives esophageal tumorigenesis in a model of the human equivalent. This model employs 4-nitroquinoline-1-oxide (4NQO) as an RS-inducing agent. Upon drinking 4NQO-containing water, Brca1 heterozygous mice formed squamous cell carcinomas of the distal esophagus and forestomach at a much higher frequency and speed (∼90 to 120 d) than did wild-type (WT) mice, which remained largely tumor free. Their esophageal tissue, but not that of WT control mice, revealed evidence of overt RS as reflected by intracellular CHK1 phosphorylation and 53BP1 staining. These Brca1 mutant tumors also revealed higher genome mutation rates than those of control animals; the mutational signature SBS4, which is associated with tobacco-induced tumorigenesis; and a loss of Brca1 heterozygosity (LOH). This uniquely accelerated Brca1 tumor model is also relevant to human esophageal squamous cell carcinoma, an often lethal tumor.

A subset of lung cancer cases shows robust signs of homologous recombination deficiency associated genomic mutational signatures.

PARP inhibitors are approved for the treatment of solid tumor types that frequently harbor alterations in the key homologous recombination (HR) genes, BRCA1/2. Other tumor types, such as lung cancer, may also be HR deficient, but the frequency of such cases is less well characterized. Specific DNA aberration profiles (mutational signatures) are induced by homologous recombination deficiency (HRD) and their presence can be used to assess the presence or absence of HR deficiency in a given tumor biopsy even in the absence of an observed alteration of an HR gene. We derived various HRD-associated mutational signatures from whole-genome and whole-exome sequencing data in the lung adenocarcinoma and lung squamous carcinoma cases from TCGA, and in a patient of ours with stage IVA lung cancer with exceptionally good response to platinum-based therapy, and in lung cancer cell lines. We found that a subset of the investigated cases, both with and without biallelic loss of BRCA1 or BRCA2, showed robust signs of HR deficiency. The extreme platinum responder case also showed a robust HRD-associated genomic mutational profile. HRD-associated mutational signatures were also associated with PARP inhibitor sensitivity in lung cancer cell lines. Consequently, lung cancer cases with HRD, as identified by diagnostic mutational signatures, may benefit from PARP inhibitor therapy.

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.

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.

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.

Precision-engineered reporter cell lines reveal ABCG2 regulation in live lung cancer cells.

Expression of the ABCG2 multidrug transporter is a marker of cancer stem cells and a predictor of recurrent malignant disease. Understanding how human ABCG2 expression is modulated by pharmacotherapy is crucial in guiding therapeutic recommendations and may aid rational drug development. Genome edited reporter cells are useful in investigating gene regulation and visualizing protein activity in live cells but require precise targeting to preserve native regulatory regions. Here, we describe a fluorescent reporter assay that allows the noninvasive assessment of ABCG2 regulation in human lung adenocarcinoma cells. Using CRISPR-Cas9 gene editing coupled with homology-directed repair, we targeted an EGFP coding sequence to the translational start site of ABCG2, generating ABCG2 knock-out and in situ tagged ABCG2 reporter cells. Using the engineered cell lines, we show that ABCG2 is upregulated by a number of anti-cancer medications, HDAC inhibitors, hypoxia-mimicking agents and glucocorticoids, supporting a model in which ABCG2 is under the control of a general stress response. To our knowledge, this is the first description of a fluorescent reporter assay system designed to follow the endogenous regulation of a human ABC transporter in live cells. The information gained may guide therapy recommendations and aid rational drug design.